Dynamic characteristics and fragmentation evolution of columnar rockfall: A case study of the Zengziyan rockfall in Chongqing, China
-
摘要: 柱状岩体崩塌具有分布范围广、破坏能力强、影响范围大的特点。2004年8月12号,重庆甑子岩W12危岩体发生崩塌,崩塌体运动距离约600 m,形成显著超前空气冲击效应,激起浮尘高度约150 m。文章基于MatDEM离散元软件对甑子岩崩塌动力特征与破碎规律进行了研究,建立了按照实际节理分布的崩塌模型,实现了崩塌全过程的模拟,并结合影像资料验证了模型的有效性,在此基础上对MatDEM进行二次开发,统计分析了崩塌过程中岩块粒径演化规律,确定了崩塌过程中的四个显著颗粒破碎时刻,分别为崩塌源区底部岩体受压破碎、中上部岩体撞击低速三角区、中部岩体撞击斜坡地面与上部岩体撞击斜坡地面。引入分形维数与双参数Weibull分布模型分析了崩塌前后颗粒破碎规律,结果显示崩塌后颗粒破碎明显,细粒颗粒占比显著增加。文章为岩体崩塌的动力特征与破碎规律的研究提供了依据。Abstract: Columnar rockfall occur widely in China and cause extreme damage to man-made structures and facilities. On 12 August 2004, Zengziyan W12 perilous rock collapsed and ran out a horizontal distance of 600 m. The rockfall caused airblast with a height of 150 m. In this paper, the dynamic characteristics and fragmentation evolution is studied based on the discrete element software MatDEM. The simulation of the rockfall process is achieved based on the discrete element model with real distribution of joints. The simulation is valid by comparing the simulation and video data. Through deeper development on MatDEM, the rock size evolution is analyzed, it was found that there were four significant breakage moments during the rockfall. The moments are rock compression breakage at the bottom of the rockfall, collision breakage of middle-upper rock mass with low-velocity triangle area, collision of middle and upper part of rock mass with ground, respectively. The fragmentation evolution of the pre- and post- rockfall is studied based on the fractal theory and dual parameters Weibull distribution. The results show that the proportion of fine particles increase obviously after the rockfall. This paper offers the basis of dynamic characteristics and fragmentation evolution of other rockfalls.
-
0. 引言
库岸滑坡是深切割高山峡谷型库岸常见的破坏形式,多集中分布于我国西南山区[1-2]。库岸滑坡由地表内外营力相互作用而形成,人类工程建设及库区水位变化则使其演化特征更为突出,通常表现为库区蓄水之后库岸下缘坡体岩层软化,从而引起上部库岸的形变破坏[3-4]。库岸滑坡受多方面因素影响,具有成因复杂、类型多样和危害巨大等特点[5-7]。在水电站库区,蓄水和泄洪等因素导致的水位变化直接影响库岸滑坡的稳定性,库岸滑坡一旦失稳会诱发一系列次生灾害,破坏区域生态系统,毁坏库区坝体和发电设施[8],严重威胁库区上下游居民生命财产安全。因此,对水电站库岸滑坡进行变形监测具有重要意义。
常规监测手段已难以识别和监测大面积的库岸滑坡形变,相对传统的精密水准测量、全球卫星导航系统(Global Navigation Satellite System, GNSS)和光学遥感技术,合成孔径雷达干涉测量(Interferometric Synthetic Aperture Radar,InSAR)技术因其大范围、全天时、全天候、高精度和高分辨率等特点,已成功应用于库岸滑坡灾害识别监测分析中[9-11]。国内外学者利用InSAR技术在库岸滑坡灾害的应用方面做了大量研究,康亚等[8]利用三类InSAR产品和DEM数据对金沙江流域乌东德水电站段进行滑坡早期识别,成功探测到多处未知和已知的滑坡体;徐帅等利用墨兰指数对SBAS-InSAR技术获取的形变点进行空间域异常值分析和聚类处理,成功识别出三峡库区巫山—奉节段高概率潜在滑坡范围[12];王振林等[13]利用SBAS-InSAR技术提取雅砻江流域锦屏一级水电站库区左岸边坡的形变特征信息,推断出大幅水位上升是诱发滑坡复活的主要因素;朱同同等[14]结合时序InSAR技术和GPS观测值分析了降雨和蓄水对三峡库区树坪滑坡变形的影响;史绪国等[15]联合分布式目标与点目标的时序InSAR技术对三峡库区藕塘滑坡进行稳定性监测;Tantianuparp等[16]联合多种SAR数据和PS-InSAR技术对三峡巴东进行滑坡探测,并将PS点时序形变与水位时间变化进行初步相关性分析;Zhou等[17]利用时序InSAR技术发现三峡库区木鱼堡滑坡变形主要发生在水库涨落期和高水位期;Liu等[18]利用SBAS-InSAR技术对三峡巴东地区进行滑坡探测并分析季节性滑坡运动与水位变化之间的相关性,上述研究证明了时序InSAR技术在库岸滑坡监测中的可靠性,可以对水电站库岸滑坡变形进行有效分析。白鹤滩水电站地处四川和云南交界,自2021年4月开始蓄水,库区水位由660 m升至825 m,上升幅度达165 m;水电站运行期间最低水位765 m,最高水位825 m,升降水位差60 m,最大库容达256×108 m3[19-20]。库区地形起伏较大、断裂构造发育,加之蓄水引起的水位变化直接影响库岸潜在滑坡的变形趋势,对水电站基础设施和上下游居民生命财产安全造成潜在威胁[21-22]。因此,亟需对白鹤滩水电站库岸潜在滑坡进行变形分析。
文章联合2019年7月3日至2021年7月28日的150景升降轨Sentinel-1 SAR数据集,采用SBAS-InSAR技术获取白鹤滩水电站库区雷达视线方向(Line of sight,LOS)形变时间序列,在分析地表形变时间演化规律和空间分布特征的基础上,结合无人机野外调查,分析白鹤滩水电站库岸潜在滑坡的变形特征,重点研究蓄水因素对库岸潜在滑坡变形趋势的影响。
1. SBAS-InSAR技术
小基线集InSAR(Small Baseline Subset InSAR,SBAS-InSAR)技术最早由Berardino和Lanari等[23]提出,该方法通过组合数据的方式获得一系列短空间基线差分干涉图,这些差分干涉图能较好地克服空间失相关现象。SBAS-InSAR技术利用奇异值分解(SVD)法求解形变速率,将被较大空间基线分开的孤立SAR数据进行连接,进一步提高观测数据的时间采样率[23]。该方法可以有效减弱大气效应,降低相位噪声和误差[24],其基本原理及流程如下:
假定已获取覆盖同一区域的按时间序列排序的
$ N + 1 $ 幅SAR影像:$$ T = {\left[ {{T_0},{T_1}, \cdots ,{T_N}} \right]^{\rm{T}}} $$ (1) 根据干涉组合规则,生成M幅干涉图且M应当满足:
$$ \frac{{N + 1}}{2} \leqslant M \leqslant \frac{{N\left( {N + 1} \right)}}{2} $$ (2) 假设以
$ {t_0} $ 作为影像获取起始时刻且$ {t_0} $ 时刻影像覆盖区域位移为0,则在去除轨道误差、平地效应及地形相位的影响后,第$ i\left( {1 \leqslant i \leqslant M} \right) $ 幅影像某像素的干涉相位可表示为:$$ \Delta {\varphi _i} = {\varphi _{{t_1}}} - {\varphi _{{t_2}}} \approx \Delta {\varphi _{{i_{\rm{def}}}}} + \Delta {\varphi _{{i_{\rm{topo}}}}} + \Delta {\varphi _{{i_{\rm{atm}}}}} + \Delta {\varphi _{{i_{\rm{noise}}}}} $$ (3) $$ \left\{ \begin{gathered} \Delta {\varphi _{{i_{\rm{def}}}}}\left( {x,r} \right) = \frac{{4\text{π} }}{\lambda }\left[ {d\left( {{t_2}} \right) - d\left( {{t_2}} \right)} \right],i = 1,2, \cdots ,m \\ \Delta {\varphi _{{i_{\rm{topo}}}}}\left( {x,r} \right) = \frac{{4\text{π} }}{\lambda } \cdot \frac{{{B_ \bot }\Delta h}}{{r{\sin}\theta }} \\ \Delta {\varphi _{{i_{\rm{atm}}}}}\left( {x,r} \right) = {\varphi _{\rm{atm}}}\left( {{t_2}} \right) - \varphi \left( {{t_1}} \right) \\ \end{gathered} \right. $$ (4) 式中:
$ \Delta {\varphi _{{i_{\rm{def}}}}} $ ——斜距向形变产生的相位;$ \Delta {\varphi _{{i_{\rm{topo}}}}} $ ——地形相位;$ \Delta {\varphi _{{i_{\rm{atm}}}}} $ ——大气延迟引起的相位;$ \Delta {\varphi _{{i_{\rm{noise}}}}} $ ——相干噪声造成的相位。利用最小二乘或者奇异值分解(SVD)对m个解缠相位进行三维时空相位解缠即可获得不同SAR时刻对应的时序形变速率。
2. 研究区概况和研究数据
本文以四川省与云南省交界白鹤滩水电站库区作为研究区域,如图1所示。研究区长约30.38 km,宽约11.65 km,总面积353.93 km2,地处横断山脉东北部、青藏高原东南边缘,区域内断裂构造发育,构造运动强烈,河谷深切,山体陡峻,地震频发[25-27]。最高海拔3556 m,最低海拔520 m,高差达3036 m,地势陡峭,致使该区存在大量滑坡、崩塌和泥石流等地质灾害隐患。
形变监测数据选用从欧州航天局(European Space Agency, ESA)免费下载的150景C波段Sentinel-1雷达影像(其中升轨数据50景,降轨数据100景并在每个时间点上下两景拼接),升降轨数据覆盖区域如图2所示。时间跨度为2019年7月3日至2021年7月28日,极化方式为VV,成像方式为IW,数据参数如表1所示。为提高影像轨道精度,引入POD精密定轨星历数据。使用日本宇宙航空研究开发机构(Japan Aerospace Exploration Agency,JAXA)发布的ALOS WORLD 3D 30 m空间分辨率的数字高程模型(Digital Elevation Model,DEM),用于去除地形相位影响,如图3所示。
表 1 Sentinel-1A数据参数Table 1. Sentinel-1A data parameters轨道方向 成像模式 极化方式 波长 波段 入射角/(°) 升轨 IW VV 5.63 C 39.44 降轨 IW VV 5.63 C 39.28 3. SBAS-InSAR技术数据处理
采用SBAS-InSAR技术,选取经镶嵌、配准和裁剪后的100景Sentinel-1A斜距单视复数(Single Look Complex,SLC)影像(升降轨数据各50景),根据时间基线和垂直基线最优原则,升轨和降轨数据分别以日期为20191216和20200204的影像作为超级主影像。设置时间基线阈值180d,空间基线为临界基线阈值的50%,共生成654和888对干涉像对。为抑制斑点噪声,设置多视数为1∶4,采用Minimum Cost Flow 解缠方法和Goldstein滤波方法进行干涉处理,将组合干涉对经过配准,调整删除不理想的数据后生成干涉图,研究区部分较理想的干涉图如图4所示。
经过轨道精炼和重去平,利用最小二乘法和奇异值矩阵分解进行形变反演,然后估算和去除大气相位,得到研究区时间序列形变信息,对时序信息地理编码后获取研究区2019年7月3日至2021年7月28日LOS方向的形变结果。如图5所示,形变速率为正值表示靠近卫星,负值表示远离卫星。对比图5(a)、(b)研究区形变结果可知,降轨数据集探测的形变信息较为丰富,主要集中在库区西岸,最大LOS向形变速率−61.425 mm/a;升轨数据集仅在库区东岸部分区域形变较为明显,最大LOS向形变速率为91.426 mm/a。升降轨数据集形变信息不一致的原因是白鹤滩水电站库区两岸地形起伏较大,山势陡峭险峻,而升轨数据飞行方向大致沿东南向西北,雷达视线方向位于右侧,降轨数据则与之相反,故利用InSAR探测形变过程中阴影、叠掩和透视收缩等几何畸变现象严重。
4. 试验结果与分析
4.1 研究区库岸典型潜在滑坡选取
对升轨和降轨数据获取的研究区形变结果进行综合解译,升轨数据库岸形变区域解译结果如图6所示,共选取库岸形变较大区域4处。结合无人机野外调查结果,发现典型潜在滑坡2处,分别用H1和H2表示;非滑坡形变区2处,分别用X1和X2表示,升轨数据详细解译结果如表2所示。
表 2 升轨数据库岸形变区域解译结果列表Table 2. List of interpretation results of shore deformation region in orbit lifting database编号 形变区域名称 最大形变速率/(mm·a−1) 形变区域类别 H1 观音岩 19.846 潜在滑坡 H2 鱼坝 18.537 潜在滑坡 X1 六城村 76.259 非滑坡形变 X2 半坡 55.947 非滑坡形变 降轨数据库岸形变区域解译结果如图7所示,共选取库岸形变较大区域6处。结合无人机野外调查结果,发现典型潜在滑坡4处,分别用H3至H6编号;非滑坡形变区2处,分别用X3和X4表示,降轨数据详细解译结果如表3所示。
表 3 降轨数据库岸形变区域解译结果列表Table 3. List of interpretation results of shore deformation region in orbit descent database编号 形变区域名称 最大形变速率/(mm·a−1) 形变区域类别 H3 观音岩 10.726 潜在滑坡 H4 清水沟 17.605 潜在滑坡 H5 鱼坝 19.326 潜在滑坡 H6 大湾子 15.888 潜在滑坡 X3 六城村 48.871 非滑坡形变 X4 半坡 61.425 非滑坡形变 对比升轨和降轨数据解译结果可以看出,非滑坡形变区X1、X2与X3、X4分别相同,潜在滑坡H1、H2与H3、H5相互对应。另外,降轨数据还解译出除上述区域以外的潜在滑坡H4和H6,同一时间段不同轨道SAR数据集探测的形变结果能够相互对应,从侧面验证了本文InSAR结果的准确性,但受时间、空间失相干因素和几何畸变影响,升降轨形变信息有所差异,说明升降轨结合的方式能够有效弥补仅利用单一轨道识别结果不全面、不准确的缺陷,提升库岸潜在滑坡灾害识别和监测的准确性和有效性。
4.2 库岸典型潜在滑坡变形分析
结合4.1节升降轨数据集库岸潜在滑坡解译结果,本文选取H1、H2、H4和H6四处典型潜在滑坡进行变形分析,分别在各滑坡形变结果中选取特征点,引入研究区降雨数据,绘制特征点在蓄水前后的时序形变曲线,并结合无人机野外调查结果分析库岸典型潜在滑坡的变形特征。
H1滑坡地处观音岩,位于沿江公路东岸,滑坡形变速率如图8(a)所示。滑坡整体形变速率范围为−10.726~15.433 mm/a,分别选取滑坡体上缘和下缘特征点A、B与降雨数据构建时序形变曲线如图8(b)所示,特征点A和B时序形变速率波动趋势大致相同,每年雨季形变速率较旱季明显增加。2019年10月—2020年5月形变速率减小,2021年4月后形变速率增大,同比增加约16 mm/a。
经实地勘察,该滑坡坡体上缘为自然坡体,坡体下缘已进行边坡加固,故在2019年10月至2020年5月间B点较A点形变速率变化相对稳定。受降雨因素影响,坡体在雨季滑动速率增大。2021年4月至5月,降雨量几乎为零,水电站蓄水导致库区水位上升,库岸下缘受到江水侵蚀改变坡体上下缘间的平衡关系,使该滑坡体形变量增大。
H2滑坡地处鱼坝村,金沙江支流末端。滑坡形变速率如图9(a)所示,形变较大值处于坡体中上部,形变范围在−19.326~8.254 mm/a。在坡体两侧分别选择特征点C、D结合降雨数据构建时序形变曲线见图9(b),特征点C和D形变速率变化趋势基本一致,与降雨数据呈现一定相关性,雨旱两季形变速率差异较小。2020年1月至10月间,形变速率逐渐减小,2021年1月后形变速率振荡变化,2021年4月之后,形变速率增加值超过10 mm/a。
经野外实地调查,H2滑坡滑面自上而下呈倒“V”字形。由于隧道工程尚未完工,附近仍伴有部分工程活动,故在2020年雨季坡体滑动速率对降雨因素响应较弱,表现为形变速率逐渐减小。从图9(b)可以看出,2021年4月以后,特征点C和D形变速率相对之前有所增加,此时降雨量较小,说明蓄水导致的库区水位抬升也对远离河道的坡体产生影响。
H4滑坡体地处清水沟,位于库区西岸。滑坡形变速率如图10(a)所示,滑坡整体形变范围为−17.605~9.012 mm/a,形变速率较大区域位于坡体中部。由图10(b)特征点与降雨数据构建的时序形变曲线可知,特征点E呈振荡变化趋势,雨旱两季形变速率差异明显。2020年8月后形变速率急剧增大,2021年4月之后形变速率相比同期增加约17 mm/a。
通过野外调查可知,H4滑坡属于临江大型冲沟,沟面呈褶皱形态,目前尚未发育为真正意义的滑坡。图10(b)时序形变曲线在2020年雨季后呈梯度下降趋势,主要原因是降水冲刷沟壑表面使冲沟坡面向下滑动。2021年4月之后相比同期形变速率明显增加,此时受降雨影响微弱,说明该滑坡体对水位变化有较强响应,原本裸露的坡体下缘遭受江水侵蚀,下缘坡体在动水压力作用下土壤结构趋向松散状态,上缘冲沟体失稳,自然产生向下形变。
H6滑坡位于库区西岸大湾子隧道,滑动面处于隧道临江一侧,其形变速率如图11所示,整体形变速率为−15.888~16.326 mm/a, 选取滑坡体中部特征点F与降雨数据建立时序形变曲线(图11),特征点F形变速率整体呈波动趋势,在2020年雨季形变速率较旱季增速明显。2021年4月之后,形变速率缓慢增大,较同期增加约16 mm/a。
经野外实地考察,发现H6滑坡已发育且有部分滑动痕迹,在坡体顶端还发育有一定程度的裂缝(图11),所以在雨季降水冲刷坡面且沿裂缝渗入坡体改变其土体应力结构,使坡体产生较大形变。由图11可以看出,2021年4—5月间,降雨量几乎为零,但形变速率变化明显,说明该坡体对水位抬升具有较强响应,降水沿裂缝进入坡体内部,促进了破裂面的贯通,而水位抬升致使坡体下缘和滑动面软化,降低其抗剪强度,降雨和水位抬升的共同作用可能使H6滑坡进一步发育,后续应当对该滑坡进行重点监测。
5. 结论
本文联合升降轨Sentinel-1 SAR数据,采用SBAS-InSAR技术并结合无人机野外调查数据,分析白鹤滩水电站库岸潜在滑坡的变形特征,得到以下结论:
(1)白鹤滩水电站库区LOS方向形变速率为−90.959~91.426 mm/a,受蓄水因素影响,各库岸典型潜在滑坡形变速率明显加快,蓄水前后形变平均增速达10 mm/a以上;
(2)白鹤滩水电站库岸潜在滑坡对水位变化具有较强响应,蓄水量增加是当前库岸潜在滑坡发育的关键性诱因,水位抬升之后潜在滑坡形变速率变化明显,在降雨和蓄水等因素共同作用下,白鹤滩水电站库岸潜在滑坡存在失稳风险;
(3)降轨数据集探测的形变信息较为丰富,主要集中在库区西岸,而升轨数据集仅在库区东岸部分区域形变较为明显,故联合升降轨SAR数据能有效克服仅利用单一轨道导致的几何畸变等问题,使水电站库岸潜在滑坡变形监测更加准确、全面。
-
表 1 模型材料力学参数
Table 1 Mechanical parameters of the Zengziyan rockfall model
材料属性 符号与单位 源区上部岩体 源区下部岩体 节理 密度 ρ/(kg·m−3) 2700 2700 2600 弹性模量 E/GPa 65 58 8 泊松比 v 0.20 0.17 0.14 抗拉强度 Cu/MPa 5.3 1.4 0.9 抗压强度 Tu/MPa 51 13 3 内摩擦系数 μi 0.87 0.7 0.7 单元直径 di/m 0.4±0.08 0.4±0.08 0.4±0.08 -
[1] CROSTA G B,IMPOSIMATO S,RODDEMAN D. Numerical modeling of 2-D granular step collapse on erodible and nonerodible surface[J]. Journal of Geophysical Research,2009,114(F3):F03020.
[2] UTILI S,ZHAO T,HOULSBY G T. 3D DEM investigation of granular column collapse:evaluation of debris motion and its destructive power[J]. Engineering Geology,2015,186:3 − 16. DOI: 10.1016/j.enggeo.2014.08.018
[3] ZHOU Y Y,SHI Z M,ZHANG Q Z,et al. 3D DEM investigation on the morphology and structure of landslide dams formed by dry granular flows[J]. Engineering Geology,2019,258:105151. DOI: 10.1016/j.enggeo.2019.105151
[4] HUANG B L,WANG J,ZHANG Q,et al. Energy conversion and deposition behaviour in gravitational collapse of granular columns[J]. Journal of Mountain Science,2020,17(1):216 − 229. DOI: 10.1007/s11629-019-5602-9
[5] 陈智强,李渝生. 重庆市南川甑子岩危岩形成演化机制分析及防治措施探讨[J]. 中国地质灾害与防治学报,2004,15(1):78 − 81. [CHEN Zhiqiang,LI Yusheng. Analysis on formation and development mechanism and discussion on prevention measures for Zenziyan dangerous rock mass in Chongqing[J]. The Chinese Journal of Geological Hazard and Control,2004,15(1):78 − 81. (in Chinese with English abstract) DOI: 10.3969/j.issn.1003-8035.2004.01.017 [6] 贺凯. 塔柱状岩体崩塌机理研究[D]. 西安: 长安大学, 2015 HE Kai. Research on collapse mechanism of tower rock[D]. Xi’an: Changan University, 2015. (in Chinese with English abstract)
[7] 贺凯,殷跃平,李滨,等. 塔柱状岩体崩塌运动特征分析[J]. 工程地质学报,2015,23(1):86 − 92. [HE Kai,YIN Yueping,LI Bin,et al. Video imaged based analysis of motion characteristic for tower rock collapse[J]. Journal of Engineering Geology,2015,23(1):86 − 92. (in Chinese with English abstract) DOI: 10.13544/j.cnki.jeg.2015.01.013 [8] 贺凯,殷跃平,冯振,等. 重庆南川甑子岩-二垭岩危岩带特征及其稳定性分析[J]. 中国地质灾害与防治学报,2015,26(1):16 − 22. [HE Kai,YIN Yueping,FENG Zhen,et al. Analysis of characteristics and stability for Zengziyan-Eryayan unstable rocks belt in Nanchuan County Chongqing[J]. The Chinese Journal of Geological Hazard and Control,2015,26(1):16 − 22. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2015.01.003 [9] 冯振,陈云霞,李滨,等. 重庆南川甑子岩山体崩塌机制研究[J]. 水文地质工程地质,2016,43(1):50 − 56. [FENG Zhen,CHEN Yunxia,LI Bin,et al. Failure mechanism on the Zengziyan collapse in Nanchuan of Chongqing[J]. Hydrogeology & Engineering Geology,2016,43(1):50 − 56. (in Chinese with English abstract) DOI: 10.16030/j.cnki.issn.1000-3665.2016.01.08 [10] 孙敬辉,石豫川. 重庆甑子岩崩塌落石动力学特征及危险性分区[J]. 中国地质灾害与防治学报,2019,30(3):6 − 11. [SUN Jinghui,SHI Yuchuan. Dynamics and hazard zoning of collapse and rockfall in Zengziyan,Chongqing[J]. The Chinese Journal of Geological Hazard and Control,2019,30(3):6 − 11. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2019.03.02 [11] 任幼蓉,陈鹏,张军,等. 重庆南川市甑子岩W12#危岩崩塌预警分析[J]. 中国地质灾害与防治学报,2005,16(2):28 − 31. [REN Yourong,CHEN Peng,ZHANG Jun,et al. Early-warning analysis on the rockfall for Zenziyan W12# dangerous rock mass in Nanchuan City of Chongqing[J]. The Chinese Journal of Geological Hazard and Control,2005,16(2):28 − 31. (in Chinese with English abstract) DOI: 10.3969/j.issn.1003-8035.2005.02.006 [12] 刘春,范宣梅,朱晨光,等. 三维大规模滑坡离散元建模与模拟研究—以茂县新磨村滑坡为例[J]. 工程地质学报,2019,27(6):1362 − 1370. [LIU Chun,FAN Xuanmei,ZHU Chenguang,et al. Discrete element modeling and simulation of 3-dimensional large-scale landslide-taking xinmocun landslide as an example[J]. Journal of Engineering Geology,2019,27(6):1362 − 1370. (in Chinese with English abstract) DOI: 10.13544/j.cnki.jeg.2018-234 [13] 奚悦. 基于离散单元法的岩石颗粒破碎研究[D]. 上海: 上海交通大学, 2016 XI Yue. Study on the rock particle crushing using discrete element method[D]. Shanghai: Shanghai Jiao Tong University, 2016. (in Chinese with English abstract)
[14] HOU T X,XU Q,ZHOU J W. Size distribution,morphology and fractal characteristics of brittle rock fragmentations by the impact loading effect[J]. Acta Mechanica,2015,226(11):3623 − 3637. DOI: 10.1007/s00707-015-1409-0
[15] TURCOTTE D L. Fractals and fragmentation[J]. Journal of Geophysical Research,1986,91(B2):1921. DOI: 10.1029/JB091iB02p01921
[16] 王健,黄波林,张全,等. 碎裂化柱状危岩体崩塌-堆积特征概化模型研究[J]. 水利水电技术,2020,51(2):136 − 143. [WANG Jian,HUANG Bolin,ZHANG Quan,et al. Study on generalized model of collapse-deposit characteristics of cataclastic and columnar dangerous rock mass[J]. Water Resources and Hydropower Engineering,2020,51(2):136 − 143. (in Chinese with English abstract) DOI: 10.13928/j.cnki.wrahe.2020.02.016 [17] LUBE G, HUPPERT H E, SPARKS R S J, et al. Collapses of two-dimensional granular columns[J]. Physical Review E, Statistical, Nonlinear, and Soft Matter Physics, 2005, 72(4 Pt 1): 041301.
[18] 郝明辉,许强,杨兴国,等. 高速滑坡-碎屑流颗粒反序试验及其成因机制探讨[J]. 岩石力学与工程学报,2015,34(3):472 − 479. [HAO Minghui,XU Qiang,YANG Xingguo,et al. Physical modeling tests on inverse grading of particles in high speed landslide debris[J]. Chinese Journal of Rock Mechanics and Engineering,2015,34(3):472 − 479. (in Chinese with English abstract) [19] KERMANI E,QIU T,LI T B. Simulation of collapse of granular columns using the discrete element method[J]. International Journal of Geomechanics,2015,15(6):4015004. DOI: 10.1061/(ASCE)GM.1943-5622.0000467
[20] 陈小婷,黄波林. FEM/DEM法在典型柱状危岩体破坏过程数值分析中的应用[J]. 水文地质工程地质,2018,45(4):137 − 141. [CHEN Xiaoting,HUANG Bolin. Application of the FEM/DEM method to numerical analyses of the failure process of representative pillar-shape dangerous rockmass[J]. Hydrogeology & Engineering Geology,2018,45(4):137 − 141. (in Chinese with English abstract) DOI: 10.16030/j.cnki.issn.1000-3665.2018.04.20 [21] 纳曼·麦麦提,米红林. 基于离散元的危岩群体崩塌影响因素分析[J]. 人民长江,2021,52(2):99 − 104. [NAMAN·Maimaiti,MI Honglin. Numerical simulation on collapsing influence factors of perilous rock groups based on discrete element method[J]. Yangtze River,2021,52(2):99 − 104. (in Chinese with English abstract) DOI: 10.16232/j.cnki.1001-4179.2021.02.016 [22] 张家勇,邹银先,杨大山. 基于PFC3D的鱼鳅坡滑坡运动过程分析[J]. 中国地质灾害与防治学报,2021,32(4):33 − 39. [ZHANG Jiayong,ZOU Yinxian,YANG Dashan. Analysis of Yuqiupo landslide motion process based on PFC3D[J]. The Chinese Journal of Geological Hazard and Control,2021,32(4):33 − 39. (in Chinese with English abstract) [23] 陶志刚,张海江,尹利洁,等. 基于FDEM的戒台寺古滑体开裂破坏过程数值模拟[J]. 水文地质工程地质,2017,44(3):105 − 112. [TAO Zhigang,ZHANG Haijiang,YIN Lijie,et al. Numerical modeling of cracking for the Jietai temple ancient landslide with the combined finite-discrete element method[J]. Hydrogeology & Engineering Geology,2017,44(3):105 − 112. (in Chinese with English abstract) DOI: 10.16030/j.cnki.issn.1000-3665.2017.03.16