Analysis of the landslide movement characteristics in Biyi Village, Ganluo County, Sichuan Province, through SPH numerical simulation
-
摘要:
四川甘洛县比依市村受连续强降雨影响下形成滑坡灾害,造成大量财产损失。文章以该滑坡为例,通过工程地质条件探究了比依市村滑坡产生的原因,采用光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法对滑坡全过程进行真三维数值模拟,对滑坡关键位置进行位移、速度监测,结果表明:(1)不良地质和长历时、高强度的持续降雨是滑坡的主要诱发因素;(2)SPH方法得到的位移、速度等运动特征参数与实际勘察效果基本一致,具有较好的可靠性;(3)该滑坡是由于滑坡后缘推动滑坡中部导致整体下错,前缘受阻形成鼓胀区,最终破坏模式表现为推移式破坏;(4)比依市村滑坡仍处于蠕变阶段,利用SPH模拟大变形的优势可以为此类灾害的机理研究和工程防治提供新的研究思路。
Abstract:Due to continuous heavy rainfall, landslides were triggered in Biyi Village, Ganluo County, Sichuan Province, resulting in significant property losses. This paper explores the causes of the landslide in Biyi Village by analyzing the local engineering geological conditions, and conducts a true three-dimensional numerical simulations of the entire landslide process using the smoothed particle hydrodynamics (SPH) method. Displacement and velocity monitoring were performed at key positions of the landslide. The results indicate that: (1) unfavorable geological conditions and long-duration, high-intensity continuous rainfall were the main triggering factors of the landslide. (2) The SPH method obtained displacement, velocity, and other motion-related parameters that closely aligned with the actual survey results, demonstrating its good reliability and accuracy. (3) The landslide was caused by the sliding of the rear edge of the landslide, pushing the middle part of the landslide, resulting in overall downward movement and blockage at the front edge, ultimately leading to the formation of a bulging zone. The final landslide failure mode was characterized as a translational slide. (4) The Biyi Village landslide is still in the creep stage. Utilizing the capabilities of the SPH method in simulating large deformations opens up new avenues for research in understanding the mechanisms and engineering strategies for the prevention of similar disasters.
-
0. 引言
地裂缝为常见的地质灾害之一,一般指在内外应力协同作用下,岩土体的连续性被破坏,并发生变形破裂,进而在浅表层形成的具一定规模的变形带或裂缝的灾害现象,其变形破坏常呈现为水平剪切、垂直剪切与拉张破坏[1 − 3]。地裂缝在许多国家都有广泛发育,在美国、中国、墨西哥、澳大利亚以及埃塞俄比亚、肯尼亚等地均出现了大规模的地裂缝灾害,对当地居民生产生活造成严重影响,产生了难以计量的破坏和损失[4 − 10]。
调查发现,自上世纪60年代起,北京平原区勘察到地裂缝超过40条,至今已发育有7个地裂缝带,其中以高丽营、顺义和通州地区所发育的地裂缝最为典型,对周围建筑和公共设施造成了严重破坏,引起了人们的高度关注[11]。在北京地裂缝的发育特征和成因机理研究方面:贾三满等认为高丽营地裂缝主要受断裂的控制,为基底断裂活动延伸至地表导致的,是断裂构造和地面沉降共同作用的结果[12];彭建兵等通过对顺义区土沟村地裂缝进行研究,认为该地裂缝与下伏断裂相连,并呈现出同沉积断层的活动特征[13];刘德成等认为通州地裂缝是非构造性地裂缝,主要受地貌形态所控制,并分布于河流两岸,顺河堤方向展布[14];姜媛等在对大量资料分析的基础上,得到了高丽营地裂缝和当地地面沉降在时间、空间和活动方式上的关联,认为二者在加剧活动时间上具有同步性,在活动特征上具有相似性,并指出地面沉降在断裂两侧引起的不均匀沉降是高丽营地裂缝形成的主要因素[15];卢全中等对北七家至高丽营段的地裂缝进行研究,确定其上、下盘影响带总宽度为43 m,地裂缝对房屋与道路的破坏形式有3种[16];刘方翠等认为是区域构造活动及差异沉降变形机制的共同作用导致了北京平原区顺义地裂缝、高丽营地裂缝、羊房地裂缝和北小营地裂缝的快速形成发育[17];周永恒等通过地球物理探测技术与联合钻孔剖面方法揭示了黄庄-高丽营断裂(房山-涞水段)具典型的正断活动特征[18];关金环等发现北京首都国际机场发生的差异性沉降主要是由顺义-良乡断裂所发育的顺义地裂缝造成的[19];孟振江等以北京地裂缝为原型研究了耦合型地裂缝的发育活动特征和成因机制,揭示了不同位错量和水位下降量引发的地层位移场和应力场的变化特征[20];任雅哲等依据断层位错理论定量分析了顺义断裂蠕滑活动对机场地裂缝形成的影响,认为顺义断裂上盘中、深层地下水抽采引起的地面差异沉降是导致机场地裂缝加剧的主要因素[21]。以上研究在一定程度上揭示了北京地裂缝的成因与影响因素,为后续的深入研究提供了重要的研究基础和理论参考。
作为特殊的工程场地,地裂缝场地对动力的响应明显且产生的影响较大。相关研究表明,傅里叶频率与场地条件、土层结构完整性及力学性质等紧密相关,而地脉动频率及卓越周期都可以反映出试验场地土体结构的变化[22 − 23]。在利用地脉动测试方法研究方面:薛捷等通过地脉动测试的方法对西安地裂缝场地进行了研究,从而为西安市抗震设防工作提供理论依据[24];慕焕东等使用FLAC3D数值分析软件对西安地裂缝场地开展研究,发现地裂缝场地具有明显的放大效应,且上盘效应显著[25];张磊刚等应用地脉动测试手段对西安F6地裂缝场地进行研究,探明了场地因地裂缝而产生的放大效应和峰值变化规律[26];王晗等通过地脉动现场测试方法,得到了西安地区典型地裂缝的场地放大效应和变化规律[27];崔思颖等选取河北平原具有代表性的地裂缝进行地脉动测试揭示了动力响应特征[28];Chang J和Xuan Y等利用地脉动测试技术获得了山西临汾和大同两处盆地的典型地裂缝场地的数据并开展研究,研究结果表明可以通过地脉动信号来分析地裂缝场地的动力响应特征,并揭示了其变化规律与放大效应[29 − 32]。
在对北京地区的地裂缝研究方面,前人所做研究大多集中于对地裂缝破坏现象的描述和产状勘察,以及从发育场地地貌和地质构造角度对其成因机理进行研究,在利用地脉动测试北京地裂缝发育场地的动力响应研究方面还属于空白,因该方法的应用还未被普及,加上北京多个地裂缝段现场的干扰较大,故尚未开展专门的地脉动现场测试[33]。鉴于此,本文的研究对象选取为北京通州区宋庄地裂缝场地,采用现场踏勘、地脉动测试和动力分析等方法,开展针对性的宋庄地裂缝场地的动力响应研究,以探究地裂缝的存在对场地动力响应的“峰值”及其“卓越频率”所造成的影响,并在此基础上揭示地裂缝发育场地的动力响应规律,以期为宋庄地区场地工程建设和防震减灾工程的设计和规划提供理论依据,同时为其他地裂缝发育区开展此类研究提供借鉴与参考。
1. 北京地裂缝分布及活动特征
北京平原区地处太行山隆起与燕山隆起和华北断陷盆地相交汇的地带,自新生代以来,该地区主要处于伸展变形的构造环境,地质构造复杂且独特,断裂构造发育。近年来,随着人类工程活动的日益加剧,引发了多起地裂缝灾害[33]。
通州区地裂缝于上世纪60年代开始零星出现,唐山大地震发生后,郎府、西集、宋庄等地都出现了大量的地裂缝,其中以宋庄地裂缝最为典型。宋庄镇地处北京迭断陷与大兴迭隆起相交界的区域,其下伏具有活动迹象的早、中更新世南苑-通县断裂,该断裂带长约110 km,走向NE35~50°,倾向NW,倾角60°~75°,NW盘下降,SE盘上升,为正断层(图1)。断层西北侧震旦亚界埋深400~
1000 余米,上覆有巨厚的新近系、第四系沉积,而东南侧寒武系及中上元古界埋深仅60~300米,上覆有第四系及少量新近系沉积;对断层附近钻孔岩芯观察,具有明显挤压现象,应属压型断层。宋庄地裂缝于2010年后活动尤为剧烈[34]。经野外调查发现,宋庄地裂缝沿南苑-通县断裂展布,南起小中河在地表呈近直线延伸,经双埠头村、沟渠村、大庞村至平家疃村后向北东方向延伸,长度约8.7 km,整体走向NE31°,倾向NW。该地裂缝以水平拉张破坏为主,地表水平拉张量在2.15~56.7 cm之间,垂直位错量较小,地裂缝影响变形带最宽处达到400 m(图2)。宋庄地裂缝的地层剖面如图3示,地表以下由粉砂、细砂、粉土、粉质黏土互层组成[35]。在浅表层,地裂缝近直立,地裂缝以偏拉张开裂为主,张开量为1~15 cm;在地下埋深50 m以下的地层开始逐渐出现位错现象,埋深53.2 m处的钻孔岩芯显示,在浅黑灰色砂质粘土层中,发育不连续构造面,构造面倾角64°,且两侧地层颜色不同。总体上分析发现地层上下盘位错量随埋深逐渐增大,地裂缝NW盘地层相对下降,西北侧第四系地层沉积厚度大于东南侧。在宋庄地裂缝影响范围内,致灾最为严重的是双埠头村至大庞村段,据统计,2016年、2018年和2022年地裂缝造成的直接经济损失总值分别为
4623.24 万元、7033.49 万元和10751.14 万元,可见地裂缝的影响和破坏在不断加剧[36]。2. 地脉动测试
2.1 地脉动基本原理
地脉动通常是指地表面振动周期在0.05~10s、振幅在千分之几至几微米的微幅振动,是由环境振动在地球表面所产生的复杂随机振动,包括地球内部应力、地震等自然因素及交通等人类活动所引起的振动,并经过不同传波介质和不同场地环境后综合作用于地球表面的结果。目前地脉动研究大致分为两类:一类是常时微动,一般是指周期小于1s的微动,侧重其振幅、周期性与观测点场地分类和振动特性研究;另一类是长波微动,一般是指周期大于1s的微动,研究的重点是由微动提取面波,进而解释推断观测点地下横波波速1s的常时微动,可用以研究地裂缝场地动力响应特性。
2.2 地脉动测试
地脉动测试要求的频段范围在0.1~15Hz之间,本次进行地脉动测试所使用的设备是高灵敏度伺服型速度网络地震仪,可满足测试精度要求,主要规格参数见表1。
表 1 选取地震仪的主要技术参数Table 1. Main technical parameters of the Seismic Monitor仪器型号 CV-374AV 动态范围 136 db 灵敏度 1000 mv/(cm/s)通频带 0.1~100 Hz 测量范围 ±0.02 m/s 线性 0.03% 采样频率 100 Hz AD分辨频率 24 bit 高灵敏度伺服型速度网络地震仪尺寸为180*120*100 mm,灵敏度较高,在进行地脉动测试时,严格按照操作规程测试,确保测试数据的精确性。为避开不利环境因素和人类活动对地脉动监测结果的影响,选择天气条件良好情况下进行地脉动测试,测试工作选择在凌晨1点至5点期间。
选择在北京通州区宋庄镇发育较为典型的双埠头村地裂缝段进行地脉动测试,现场布设测点具体原则为:在与主裂缝走向垂直的方向布设2条测线,结合地裂缝的影响带宽度范围,布设长度为60 m的测线,其中,上、下盘长度各设30 m。将地裂缝作为测线的中点,向远离地裂缝的方向进行测点的布设。以上盘为例,首先,在距地裂缝1.5 m处布设第一个测点,再沿测线以3 m为间距布设测点,共布设5个;接着按每隔5 m布设一个测点,共布设3个。因此,所布设测点与地裂缝间的水平距离由近到远设成1.5、3、6、9、12、15、20、25和30 m。下盘测点与上盘对称布设。
针对地脉动观测各个测点的数据,在稳定波段内以10秒为单位进行波段的截取并通过Pwave32软件将其转换导出为加速度时程曲线,之后应用SeismoSignal软件对数据进行预处理,以得到反应谱曲线、傅里叶谱曲线以及Arias烈度曲线,通过频谱特征分析,进而获得地裂缝场地地脉动响应的加速度幅值、卓越频率、影响范围等场地动力响应特征;通过利用3种从不同角度出发的动力分析方法,对采集并经处理的地脉动数据进行分析,进一步提炼地脉动产生的场地动力特性。
3. 测试结果分析
3.1 傅里叶分析
通过傅里叶谱分析,可以将地脉动加速度时程曲线转换成傅里叶谱曲线,从而得到场地卓越频率和傅里叶谱峰值,通过对不同测点卓越频率和峰值的对比,可得到场地不同区域地脉动响应的分布特征,从而获得地裂缝场地动力响应特征。
3.1.1 双埠头村地裂缝傅里叶分析
图4、图5为双埠头村地裂缝测线上盘和下盘在三个方向上的傅里叶谱图(其中X、Y、Z分别为垂直地裂缝走向方向、平行地裂缝走向方向及和垂直地面方向,下文同),图中各测点的傅里叶谱峰值显著,呈现出多峰的特点。主频带宽度较大,而各测点的卓越频率相对集中,因此其波速的变化没有明显的规律。傅里叶谱图反映出的多峰、谱面积大以及谐波丰富等特征,这与双埠头村地裂缝发育场地中较为松软的粘性土层结构特征相一致。通过对双埠头村地裂缝上下盘的傅里叶谱统计分析,得到该区域地裂缝场地的卓越频率在2.10~4.30 Hz之间浮动。而场地卓越频率的大小和测点距地裂缝的水平距离间的关系并不显著,均值约为3.30 Hz。
图6、图7为双埠头村地裂缝各测点的傅里叶卓越峰值图。从图中可以看出,在距离地裂缝区域较近的位置,地脉动响应由于地裂缝的存在而出现了显著的放大效应,傅里叶谱峰值与该距离呈负相关关系,距地裂缝的距离越大,峰值越低,在距离较远处曲线趋于平缓,其临界值约在距地裂缝12 m处,傅里叶谱峰值大致保持稳定不变。在距离地裂缝12 m以内的区域,振动的幅值和能量密度相对较高,而在距离地裂缝12米以外的区域,地表振动的幅值和能量密度则相对较低,这种现象表明靠近地裂缝的区域地表振动会比远离地裂缝的区域更强烈。
3.1.2 沟渠村地裂缝傅里叶分析
图8、图9分别为沟渠村地裂缝各测点X、Y、Z三个方向上傅里叶谱峰值与距地裂缝的距离变化关系的散点及曲线图,发现距离地裂缝越远的位置则傅里叶谱峰值越小,直至在较远的位置,曲线逐渐接近平缓;最终在距地裂缝约12 m处,曲线基本保持不变,傅里叶谱峰值趋于一致。通过统计分析,得出沟渠村地裂缝的场地卓越频率在2.30~4.30 Hz范围内变化,卓越频率的高低与测点所在的位置之间没有明显关系,其均值约为3.40 Hz。
3.2 反应谱分析
通过反应谱分析,可以得到特征周期和相应的反应谱峰值,通过对不同测点特征周期和反应谱峰值的对比,可以得到场地不同区域地脉动响应的分布特征,从而获得地裂缝场地动力响应特征。
3.2.1 双埠头村地裂缝反应谱分析
经数据分析得到,双埠头村地裂缝各测点在同一方向上的反应谱谱型均为“单峰型”,且具有主峰突出、频带窄和谱面积小的特点;反应谱卓越周期集中在0.1 s~0.2 s之间,部分卓越周期可达0.3s,此外卓越周期值与地裂缝和测点间的水平距离无明显关系(图10、11)。
由图12、图13可得,在距离地裂缝较近处,因地裂缝的存在地脉动响应表现出显著的放大效应,随着与地裂缝距离的增大,地脉动的响应加速度基本不变,在距地裂缝12m处趋于平稳,反应谱峰值基本一致。
3.2.2 沟渠村地裂缝反应谱分析
由图14、图15可以看出,沟渠村地裂缝各测点在X、Y、Z三个方向的上反应谱加速峰值随与地裂缝间的垂直距离变化,距地裂缝较近的区域的地脉动响应表现出显著的放大效应,且距离地裂缝越远,反应谱峰值越小,直至曲线接近稳定不变,反应谱峰值基本一致。
通过对比双埠头村地裂缝和沟渠村地裂缝反应谱曲线,发现布设在宋庄地裂缝上下盘的测点,在与地裂缝同等距离时,上盘反应谱的峰值均大于下盘,这表明地裂缝的存在使得地裂缝两侧场地的地脉动响应发生了差异变化,存在显著“上盘效应”,且上盘峰值大于下盘。
3.3 Arias烈度分析结果
Arias烈度是一种目前较为成熟的评判场地地震动响应强弱、反映场地抗震能力和设防标准的方法,通过分析Arias烈度曲线的变化特征进而了解地裂缝场地不同位置的抗震特性,可从场地抗震烈度方面揭示地裂缝对场地的放大效应。
3.3.1 双埠头村地裂缝Arias烈度分析
图16、图17分别为双埠头村地裂缝上盘与下盘不同测点的Arias烈度曲线图,图中注释栏“MP”表示监测点。结果表明,Arias烈度呈现出随时间逐渐增长的趋势,但各个测点的增长幅度有所不同,且有着不同的最大值。
为反映地裂缝对场地不同位置烈度值产生的影响,将双埠头村地裂缝各测点Arias烈度的峰值沿垂直地裂缝走向方向展开,得到上下盘各测点在X、Y、Z三个方向上的Arias烈度峰值随与地裂缝距离变化的散点数据曲线图(图18、19)。可见,场地不同位置的烈度值受地裂缝的影响程度也不同。Arias烈度在距离地裂缝较近的区域,其峰值较大,且该烈度值随测点与地裂缝间的距离增大而减小并逐渐趋于稳定,其临界距离为12 m,而临界区域以外的范围受放大效应的影响较小。
4. 结论
(1)宋庄地裂缝延伸8.7 km,地表破坏以拉张开裂为主,水平拉张量在2.15~56.7 cm之间;其中双埠头村地裂缝的傅里叶卓越频率主要集中在2.10~4.30 Hz之间,反应谱的卓越周期集中在0.1~0.2 s之间;沟渠村地裂缝的傅里叶卓越频率集中在2.30~4.30 Hz之间,反映了地裂缝不同区段活动性的变化。
(2)通过反应谱分析法、傅里叶谱分析法分析得到,测点的卓越频率无明显的规律性变化;Arias烈度峰值在靠近地裂缝处较大,且随着距离的增大而逐渐减小,在距地裂缝较远的区域趋于平缓。
(3)宋庄地裂缝场地动力响应在反应谱、傅里叶谱和Arias烈度上都存在明显的“放大效应”,即靠近地裂缝的区域大于远离地裂缝的区域,且该临界距离为距地裂缝12 m。
(4)宋庄地裂缝场地有着较为显著的“上盘效应”,即距地裂缝相同距离情况下,位于上盘的傅里叶卓越频率、反应谱卓越周期和Arias烈度峰值总大于下盘。
-
表 1 岩土体参数
Table 1 Rock and soil mass parameters
名称 密度/(kg·m−3) 弹性模量/MPa 泊松比 黏聚力/kPa 内摩擦角/(°) 滑床 1846 20000 0.25 26 28 滑体 1846 80 0.35 14 18 滑带 - - - 7.51 10.74 -
[1] 郑琅,张欣,王立娟. 四川省甘洛县山体滑坡应急调查与成因机制分析[J]. 人民长江,2022,53(8):117 − 122. [ZHENG Lang,ZHANG Xin,WANG Lijuan. Emergency investigation and formation mechanism of landslide in Ganluo County,Sichuan Province[J]. Yangtze River,2022,53(8):117 − 122. (in Chinese with English abstract)] ZHENG Lang, ZHANG Xin, WANG Lijuan . Emergency investigation and formation mechanism of landslide in Ganluo County, Sichuan Province[J]. Yangtze River,2022 ,53 (8 ):117 −122 . (in Chinese with English abstract)[2] 范刚,聂锐华,陈骎,等. 2020年甘洛县黑西洛沟山洪-滑坡-堰塞湖灾害链成灾机理分析[C]//. 中国水利学会2021学术年会论文集第二分册,2021:291 − 296. [FAN Gang,NIE Ruihua,CHEN Qin,et al. An analysis of the disaster mechanism of the Mountain flood,landslides and barrier Lake disaster chain in Hexiluogou,Ganluo County,in 2020[C]//China water conservancy society 2021 academic conference proceedings of the second volume,2021:291 − 296. (in Chinese)] FAN Gang, NIE Ruihua, CHEN Qin, et al. An analysis of the disaster mechanism of the Mountain flood, landslides and barrier Lake disaster chain in Hexiluogou, Ganluo County, in 2020[C]//China water conservancy society 2021 academic conference proceedings of the second volume, 2021: 291 − 296. (in Chinese)
[3] ZHOU Jinxing,ZHU Chunyun,ZHENG Jingming,et al. Landslide disaster in the loess area of China[J]. Journal of Forestry Research,2002,13(2):157 − 161. DOI: 10.1007/BF02857244
[4] HUANG Da,MENG Qiujie,SONG Yixiang,et al. Dynamic process analysis of the Niumiangou landslide triggered by the Wenchuan earthquake using the finite difference method and a modified discontinuous deformation analysis[J]. Journal of Mountain Science,2021,18(4):1034 − 1048. DOI: 10.1007/s11629-020-6188-y
[5] LIU Bo, HU Xiewen, HE Kun, et al. The starting mechanism and movement process of the coseismic rockslide:A case study of the Laoyingyan rockslide induced by the “5•12” Wenchuan earthquake[J]. Journal of Mountain Science,2020,17(5):1188 − 1205.
[6] BERNANDER S,KULLINGSJÖ A,GYLLAND A S,et al. Downhill progressive landslides in long natural slopes:triggering agents and landslide phases modeled with a finite difference method[J]. Canadian Geotechnical Journal,2016,53(10):1565 − 1582. DOI: 10.1139/cgj-2015-0651
[7] 邹德兵,丁刚,熊瑶,等. 贵州响水水库下汤章滑坡稳定性研究[J]. 水利水电技术,2020,51(增刊1):250 − 255. DOI:10.13928/j.cnki.wrahe. [ZOU Debing,DING Gang,XIONG Yao,et al. Study on the stability of Xiatangzhang landslide in Xiangshui Reservoir,Guizhou Province [J]. Water Resources and Hydropower Engineering,2020(sup 1):250 − 255. (in Chinese with English abstract)] ZOU Debing, DING Gang, XIONG Yao, et al. Study on the stability of Xiatangzhang landslide in Xiangshui Reservoir, Guizhou Province [J]. Water Resources and Hydropower Engineering, 2020(sup 1): 250 − 255. (in Chinese with English abstract) DOI: 10.13928/j.cnki.wrahe.2020.S1.044
[8] 王伟,王卫,戴雄辉. 四川美姑拉马阿觉滑坡复活特征与影响因素分析[J]. 中国地质灾害与防治学报,2022,33(4):9 − 17. [WANG Wei,WANG Wei,DAI Xionghui. Analysis of reactivated characteristics and influencing factors of the Lamajue landslide in Meigu County of Sichuan Province[J]. The Chinese Journal of Geological Hazard and Control,2022,33(4):9 − 17. (in Chinese with English abstract)] WANG Wei, WANG Wei, DAI Xionghui . Analysis of reactivated characteristics and influencing factors of the Lamajue landslide in Meigu County of Sichuan Province[J]. The Chinese Journal of Geological Hazard and Control,2022 ,33 (4 ):9 −17 . (in Chinese with English abstract)[9] 张家勇,邹银先,杨大山. 基于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)] 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)[10] 王明辉, 曹熙平, 谯立家. 危岩体精细调查与崩塌过程三维场景模拟——以西南某水电站高边坡为例[J]. 中国地质灾害与防治学报,2023,34(6):86 − 96. [WANG Minghui, CAO Xiping, QIAO Lijia. Comprehensive analysis of hazardous rock mass and simulation of potential rockfall processes using 3D terrain model: A case study of the high cut slope near damsite of a hydropower station in southern China[J]. The Chinese Journal of Geological Hazard and Control,2023,34(6):86 − 96. (in Chinese with English abstract)] DOI: 10.16031/j.cnki.issn.1003-8035.2021.01.03 WANG Minghui, CAO Xiping, QIAO Lijia . Comprehensive analysis of hazardous rock mass and simulation of potential rockfall processes using 3D terrain model: A case study of the high cut slope near damsite of a hydropower station in southern China[J]. The Chinese Journal of Geological Hazard and Control,2023 ,34 (6 ):86 −96 . (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2021.01.03[11] 张龙,唐辉明,熊承仁,等. 鸡尾山高速远程滑坡运动过程PFC3D模拟[J]. 岩石力学与工程学报,2012,31(增刊1):2601 − 2611. [ZHANG Long,TANG Huiming,XIONG Chengren,et al. Movement process simulation of high-speed long-distance Jiweishan landslide with PFC3D[J]. Chinese Journal of Rock Mechanics and Engineering,2012,31(Sup 1):2601 − 2611. (in Chinese with English abstract)] ZHANG Long, TANG Huiming, XIONG Chengren, et al . Movement process simulation of high-speed long-distance Jiweishan landslide with PFC3D[J]. Chinese Journal of Rock Mechanics and Engineering,2012 ,31 (Sup 1 ):2601 −2611 . (in Chinese with English abstract)[12] 丛凯, 魏洁, 杨亚兵, 等. 基于坡表变形分析与降雨响应模拟的立节北山滑坡运动特征[J]. 地质科技通报,2022,41(6):54 − 65. [CONG Kai, WEI Jie, YANG Yabing, et al. Investigation of the kinematic characteristic of Lijie Beishan landslide through surface displacement monitoring and rainfall response numerical simulation[J]. Bulletin of Geological Science and Technology,2022,41(6):54 − 65. (in Chinese with English abstract)] CONG Kai, WEI Jie, YANG Yabing, et al . Investigation of the kinematic characteristic of Lijie Beishan landslide through surface displacement monitoring and rainfall response numerical simulation[J]. Bulletin of Geological Science and Technology,2022 ,41 (6 ):54 −65 . (in Chinese with English abstract)[13] 廖德武, 郑冰, 杜艳松, 等. 兴仁“6•10” 彭家洞高速滑坡运动特征与形成机理[J]. 地质科技通报,2022,41(6):66 − 76. [LIAO Dewu, ZHENG Bing, DU Yansong, et al. Movement characteristics and formation mechanism of the “6•10” Pengjiadong high speed landslide in Xingren[J]. Bulletin of Geological Science and Technology,2022,41(6):66 − 76. (in Chinese with English abstract)] LIAO Dewu, ZHENG Bing, DU Yansong, et al . Movement characteristics and formation mechanism of the “6•10” Pengjiadong high speed landslide in Xingren[J]. Bulletin of Geological Science and Technology,2022 ,41 (6 ):66 −76 . (in Chinese with English abstract)[14] 张巍,史卜涛,施斌,等. 土质滑坡运动全过程物质点法模拟及其应用[J]. 工程地质学报,2017,25(3):815 − 823. [ZHANG Wei,SHI Butao,SHI Bin,et al. Material point method for Run-out process simulation of soil landslides and application[J]. Journal of Engineering Geology,2017,25(3):815 − 823. (in Chinese with English abstract)] DOI: 10.13544/j.cnki.jeg.2017.03.029 ZHANG Wei, SHI Butao, SHI Bin, et al . Material point method for Run-out process simulation of soil landslides and application[J]. Journal of Engineering Geology,2017 ,25 (3 ):815 −823 . (in Chinese with English abstract) DOI: 10.13544/j.cnki.jeg.2017.03.029[15] 王升,曾鹏,李天斌,等. 土质滑坡失稳、运动及冲击压力物质点法模拟研究[J]. 工程地质学报,2022,30(4):1362 − 1370. [WANG Sheng,ZENG Peng,LI Tianbin,et al. Initiation,movement and impact simulation of soil land-slide with material point method[J]. Journal of Engineering Geology,2022,30(4):1362 − 1370. (in Chinese with English abstract)] WANG Sheng, ZENG Peng, LI Tianbin, et al . Initiation, movement and impact simulation of soil land-slide with material point method[J]. Journal of Engineering Geology,2022 ,30 (4 ):1362 −1370 . (in Chinese with English abstract)[16] 胡嫚,谢谟文,王立伟. 基于弹塑性土体本构模型的滑坡运动过程SPH模拟[J]. 岩土工程学报,2016,38(1):58 − 67. [HU Man,XIE Mowen,WANG Liwei. SPH simulations of post-failure flow of landslides using elastic-plastic soil constitutive model[J]. Chinese Journal of Geotechnical Engineering,2016,38(1):58 − 67. (in Chinese with English abstract)] DOI: 10.11779/CJGE201601005 HU Man, XIE Mowen, WANG Liwei . SPH simulations of post-failure flow of landslides using elastic-plastic soil constitutive model[J]. Chinese Journal of Geotechnical Engineering,2016 ,38 (1 ):58 −67 . (in Chinese with English abstract) DOI: 10.11779/CJGE201601005[17] DAI Zili. 3D numerical modeling using smoothed particle hydrodynamics of flow-like landslide propagation triggered by the 2008 Wenchuan earthquake[J]. Engineering Geology,2014,180:21 − 33. DOI: 10.1016/j.enggeo.2014.03.018
[18] HUANG Yu,ZHANG Weijie,XU Qiang,et al. Run-out analysis of flow-like landslides triggered by the Ms 8.0 2008 Wenchuan earthquake using smoothed particle hydrodynamics[J]. Landslides,2012,9(2):275 − 283. DOI: 10.1007/s10346-011-0285-5
[19] PRAVEEN KUMAR R,DODAGOUDAR G R. Meshfree analysis of two-dimensional contaminant transport through unsaturated porous media using EFGM[J]. International Journal for Numerical Methods in Biomedical Engineering,2010,26(12):1797 − 1816. DOI: 10.1002/cnm.1266
[20] GOODIN C, PRIDDY J D. Comparison of SPH simulations and cone index tests for cohesive soils[J]. Journal of Terramechanics,2016,66:49 − 57.
[21] LIANG Dongfang,HE Xuzhen. A comparison of conventional and shear-rate dependent Mohr-Coulomb models for simulating landslides[J]. Journal of Mountain Science,2014,11(6):1478 − 1490. DOI: 10.1007/s11629-014-3041-1
[22] BUI H H, SAKO K, FUKAGAWA R. Numerical simulation of soil–water interaction using smoothed particle hydrodynamics (SPH) method[J]. Journal of Terramechanics,2007,44(5):339 − 346. DOI: 10.1016/j.jterra.2007.10.003
[23] GUO W B,SUN Z W,WANG C S. Slope stability numerical calculation using geotechnical strength yield criterion and sph instability criterion analysis method. Fresenius environmental bulletin. 2020,29(7A):5739 − 5745.
[24] WANG Haibin, YAN Fei, ZHANG Liwei, et al. Mechanism and flow process of debris avalanche in mining waste dump based on improved SPH simulation[J]. Engineering Failure Analysis,2022,138:106345.
[25] FÁVERO NETO A H,ASKARINEJAD A,SPRINGMAN S M,et al. Simulation of debris flow on an instrumented test slope using an updated Lagrangian continuum particle method[J]. Acta Geotechnica,2020,15(10):2757 − 2777. DOI: 10.1007/s11440-020-00957-1
[26] DAI Zili, HUANG Yu, CHENG Hualin, et al. 3D numerical modeling using smoothed particle hydrodynamics of flow-like landslide propagation triggered by the 2008 Wenchuan earthquake[J]. Engineering Geology,2014,180:21 − 33.
[27] 费康,张建伟. ABAQUS在岩土工程中的应用[M]. 北京:中国水利水电出版社,2010. [FEI Kang,ZHANG Jianwei. Application of ABAQUS in geotechnical engineering[M]. Beijing:China Water & Power Press,2010. (in Chinese)] FEI Kang, ZHANG Jianwei. Application of ABAQUS in geotechnical engineering[M]. Beijing: China Water & Power Press, 2010. (in Chinese)