Triggering types, characteristics and disaster mechanism of seismic loess landslides
-
摘要: 基于现场调查勘探、无人机航测和大型振动台试验，系统研究了黄土地震滑坡的触发类型、特征与成灾机制。结果表明，黄土地震滑坡在空间分布、单体规模、致灾范围、平面形态、地形水文条件、地震动强度、土层厚度、与发震断层关系等方面具有显著特征。黄土地震滑坡从触发机理角度可划分为剪切型滑坡、液化型滑坡、震陷型滑坡三种类型。剪切型滑坡根据滑动面地层岩性可进一步分为黄土层内滑坡、黄土与泥岩接触面滑坡、切入基岩的切层滑坡；液化型滑坡可根据液化层位置划分为底部液化滑移型、表层液化泥流型、底-表层联合液化滑流型等。震陷型滑坡可根据坡体破坏形式细分为陷滑型、崩滑型等两种。本文可为黄土地震滑坡的风险评估与防控提供科学依据。Abstract: Based on the field investigation and exploration, unmanned aerial survey and large-scale shaking table tests, the triggering types, characteristics and disaster-generating mechanism of seismic loess landslides are systematically studied. The results show that the earthquake-induced loess landslides have their distinctive characteristics in spatial distribution, single size, influencing area, plane modality, topographical and hydrological conditions, seismic intensity, deposit thickness and relations to seismic faults. They can be classified into three types from the perspective of the triggering mechanism: shear landslides, liquefaction landslides and seismic subsidence landslides. The shear landslides can be further classified according to the lithology of the sliding surface strata into three types: landslides within a loess layer, landslides on the contact surface between loess and mudstone, and landslides cutting into bedrock. The liquefaction landslides can be divided according to the location of the liquefaction layer into three types: deep liquefaction sliding type, surface liquefaction mudflow, and combined deep-surface liquefaction type. The seismic subsidence landslides can be divided into two types of landslides, subsidence slide and avalanche slide, according to the damage modes caused by seismic subsidence. This study may provide a scientific basis for the risk assessment, prevention and control of loess seismic landslides.
-
0. 引言
岩石材料在地质构造及矿物成形过程中,矿物晶体结构一般具有较明显的非均质性特点,在成分组成、几何尺寸及力学强度[1]等方面具有较大差异。在外力作用下,岩石矿物可能发生晶体尺度局部应力集中,并导致破坏的发生。岩石破坏后,其宏观裂纹在微观尺度上主要分为两大类：①沿矿物晶体边界发展的晶间裂纹；②矿物晶体内部发展的晶内裂纹（图1）。
近年来,针对岩石微观结构非均质性对宏观力学特性的影响,国内外学者对非均质性岩石进行了较多的试验研究[5-7]。Tugrul等[8]通过试验及数据相关性分析发现,矿物平均晶体尺寸及矿物组成对岩石单轴压缩强度（UCS）有较大的影响；Ündül[9]通过研究矿物组分与岩石强度之间的关系指出,斜长石矿物含量的增加将导致岩石UCS增加。此外,岩石特征应力与微观矿物结构也具有一定的关系。Ündül等[10]研究表明岩石起裂应力
受矿物成分影响较大,随斜长石矿物含量增加, 呈增加的变化趋势。Cowie等[2]对花岗岩特征应力进行统计并拟合发现,云母矿物含量与 及 具有负相关的关系,且随岩石中最小矿物尺寸的增加, 与 的比值呈增加的变化趋势。然而,Eberhardt等[11]对比3种不同粒径花岗岩单轴压缩试验结果指出,岩石中矿物晶体的尺寸对 影响不大,但随矿物晶体尺寸增加, 呈减小的变化趋势。目前,尽管国内外学者对非均质性岩石宏观力学特性进行了较多试验研究,但天然岩石材料微观矿物结构的随机性使得试验结果表现出较大的离散性,故通过试验手段定量化研究非均质岩石宏观力学仍有一定的困难。 Cundall等[12]提出的离散单元法（discrete element method, DEM）可在有效研究岩石宏观力学特性的前提下,从微观尺度对岩石内部裂纹的起裂、扩展及聚合规律进行研究。尽管颗粒流软件（particle flow code, PFC）先后开发的两种团簇模型cluster和clump模型[13-15]可有效克服传统圆盘颗粒所带来的摩擦角及强度值较小问题,但上述模型在研究矿物晶体结构对宏观力学特性影响方面仍具有不足,如cluster模型不能有效控制颗粒簇尺寸,clump模型不能模拟晶体内破裂等；经过多年发展,PFC中矿物晶体模型（PFC-GBM）的提出使得通过岩石内部矿物晶体组成、尺寸等信息进行定量化晶体建模成为了可能,且通过对矿物晶体赋予不同的微观参数即可定量化模拟非均质性岩石宏观力学特性[16-18]。周喻等[19]利用上述模型从细观角度揭示了岩石加载条件下的破裂机制及强度特性。Bewick等[20]研究了非均质岩石直剪条件下宏观力学特性及微裂纹扩展规律。Hofmann等[21]和Peng等[22]研究了花岗岩单轴及三轴压缩条件下岩石内部微裂纹的变化规律。综上,采用GBM对非均质性岩石进行矿物晶体尺度模拟是一种十分有效的方法。
本文结合花岗岩室内试验,采用颗粒流矿物晶体模型（PFC-GBM）对花岗岩矿物晶体结构进行建模,通过对比试验与模拟声发射结果并分析晶体尺度微裂纹发展规律研究了单轴压缩条件下非均质岩石微观力学特性变化规律。同时,改变岩石矿物组成进行了一系列单轴压缩试验模拟,分析并讨论了矿物结构非均质性对岩石宏观力学特性的影响及变化原因。
1. PFC-GBM简介及多尺度断裂机理
1.1 PFC-GBM简介
Potyondy[23]将Voronoi多边形网格引入PFC软件,通过在矿物晶体内赋予平行黏结模型(Parallel bonded model, PBM),矿物晶体间赋予光滑节理模型(Smooth joint model,SJM)可实现非均质岩石内微观矿物晶体建模,其建模流程主要包括以下3个步骤（图2）：①生成大尺寸颗粒以控制矿物晶体尺寸及含量（图2（a））；②连接接触多边形质心形成矿物晶体边界（图2（b）,2（c））；③在晶体边界内填充小尺寸颗粒并进行分组形成矿物晶体模型（图2（d））。
尽管上述方法的提出使得GBM在PFC软件中实现流程化,但并不能定量化控制岩石的矿物结构,尤其晶体组成、尺寸大小。因此,本文采用其他方法对岩石微观结构进行定量化建模,并将模型数据与PFC软件相结合,成功实现定量化研究岩石非均质性对宏观力学特性的影响。
1.2 多尺度断裂机理
根据断裂机理不同,GBM颗粒间接触破坏所产生的晶体尺度裂纹（晶内/晶间裂纹）可进一步分为拉伸裂纹和剪切裂纹。当矿物晶体内颗粒法向拉应力大于拉伸强度时,颗粒间接触将发生拉伸破坏并产生晶内拉伸裂纹；当矿物晶体内颗粒切向剪切应力大于剪切强度时,颗粒间接触将发生剪切破坏并产生晶内剪切裂纹。同理,矿物晶体间拉伸及剪切裂纹断裂原理与上述相同。
2. 花岗岩模型建立及微观参数确定
2.1 花岗岩模型建立
图3为花岗岩室内试验示意图,试样采用尺寸为50 mm×100 mm的标准圆柱试样。单轴压缩试验在中国科学院武汉岩土力学研究所MTS815.04型电液伺服岩石试验系统上进行,轴向加载速率设置为0.001 mm/s。
试验过程中,采用轴向引伸计和环向应变规实时测量岩石的轴向应变及环向应变,同时采用声发射装置对岩石破裂过程中声发射特征参数进行监测与采集。
结合室内试验,本文通过Neper软件对花岗岩非均质结构进行建模。建模流程如图4,主要分为3个步骤,首先在50 mm×100 mm的模型区域内生成直径大小服从均匀分布的圆盘颗粒,颗粒半径最小为0.15 mm,且最大与最小半径的比值为1.66（步骤1）；随后在Neper软件中生成矿物晶体多边形,各矿物晶体尺寸使用软件内相关命令设置为服从正态分布（步骤2）；最后将矿物晶体信息以相应的格式导入PFC软件,根据步骤2所生成的矿物晶体种类将步骤1所生成的颗粒分组以形成GBM模型（步骤3）。
2.2 微观参数确定
结合文献[21, 22, 24, 25]中GBM参数校正方法及模拟效果,模型中矿物晶体力学特性（弹性模量、微观拉伸强度及黏结强度等）设置满足以下关系：石英矿物>长石矿物>云母矿物>其他矿物。另外,由于光滑节理模型并不能较好地模拟实际矿物晶体间凹凸不平的接触面及矿物界面间的内锁问题[26-27],本文采用Li等[16]的方法对矿物晶体内及晶体间接触赋予平行黏结模型（PBM）,且晶体间微观参数采用表2公式进行折减以简化GBM的参数校正过程。使用“试错法”将微观参数进行校正后,得模拟与试验力学特性比对表见表1。
表 2 模型微观参数表Table 2. Micro-parameters of grain-based model参数 石英矿物 长石矿物 云母矿物 其他矿物 矿物晶体间参数 体积含量/% 24 49 24 3 矿物晶体尺寸d/mm 2.43±0.25 3.42±0.25 1.98±0.25 1.17±0.25 颗粒最小半径Rmin/mm 0.15 颗粒半径比Rmax/Rmin 1.66 密度/(kg·m-3) 2650 2600 3050 1650 线性弹性模量E/GPa 51 42 33 24 0.2×2E1E2/(E1+E2) 摩擦系数μ 0.35 0.41 0.33 0.44 0.8 线性刚度比kn/ks 1.1 1.3 1.7 3.7 1.3 平行弹性模量pb_emod/GPa 51 42 33 24 0.2×2E1E2/(E1+E2) 黏结刚度比kn/ks 1.1 1.3 1.7 3.7 1.3 拉伸强度pb_ten/MPa 126±16 105±16 98±13 77±9 21±0 黏结强度pb_coh/MPa 196±42 162±28 146±22 105±0 50±0 摩擦角φ/(°) 19.5 22.4 17.3 23.7 38.6 注： 表中E1,E2分别为晶间两侧矿物颗粒的接触模量值。表 1 宏观力学特性比对表Table 1. Comparison of mechanical properties宏观力学特性 试验结果 模拟结果 误差/% UCSσc/MPa 172 171 0.6 拉伸强度σt/MPa 13.6 14.3 5.1 弹性模量E/GPa 52.4 55.4 5.7 泊松比 0.229 0.228 0.44 峰值应变εp/% 0.36 0.34 5.6 由表可知,使用矿物晶体模型所得宏观力学特性与试验所得结果误差较小,故该模型可在保证有效模拟宏观力学特性的前提下,进一步研究岩石非均质性对宏观力学特性的影响,表2为该模型微观参数表。
3. 矿物晶体模型模拟结果
3.1 应力-应变曲线
图5为花岗岩试验及模拟所得应力-应变曲线。由图5（a）可知,岩石单轴压缩试验力学响应包括以下5个阶段：（Ⅰ）裂纹闭合阶段；（Ⅱ）弹性变形阶段；（Ⅲ）裂纹稳定扩展阶段；（Ⅳ）裂纹不稳定扩展阶段；（Ⅴ）峰后阶段。通过分析裂纹体积应变及试样体积应变的变化关系,可得花岗岩试样特征应力,即闭合应力
、起裂应力 、损伤应力 分别为44.8,71及133.2 MPa。 由于没有考虑材料内部的初始缺陷,模拟结果并不能有效地模拟岩石裂纹闭合阶段,故本文模拟所得力学响应存在4个阶段：（Ⅱ）弹性变形阶段；（Ⅲ）裂纹稳定扩展阶段；（Ⅳ）裂纹不稳定扩展阶段；（Ⅴ）峰后阶段。结合文献[28～30]判断起裂应力及损伤应力的方法,当AE事件率开始变化时,此时所对应的轴向应力可认为是起裂应力
；随着轴向加载的进行,AE事件率开始进入稳定变化段,此时岩石内部微裂纹开始稳定起裂并扩展；当AE事件率变化速率快速变化且在后一段加载过程中波动性较大时,可认为此临界应力为损伤应力 ,此时岩石内部微裂纹开始不稳定扩展及聚合。本文通过上述方法对特征应力进行判断,可得模拟过程中岩石 及 分别为79.2,143.3 MPa。该结果与试验所得特征应力相差在10 MPa以内,差异性较小。 3.2 声发射演化规律
试验过程中,通过对声发射特征参数进行空间定位处理可反演岩石真实破裂过程。模拟过程中,本文借鉴文献[31, 32]中声发射事件的确定方法将一个微裂纹的产生对应于1次声发射事件,通过对岩石内微裂纹空间位置进行监测并记录,可得到岩石的破裂过程。结果表明,模拟与试验所得声发射事件具有相似的变化规律（图6）。当轴向应力达到岩石起裂应力时,试样进入裂纹稳定扩展阶段,此时微裂纹开始起裂且在空间上呈随机分布状态；当轴向应力达到岩石损伤应力时,试样内微裂纹进一步起裂并发生非稳定扩展,此时微裂纹发生聚集成簇的现象；一旦轴向应力超过岩石峰值应力,微裂纹聚集成簇的程度增加进而导致岩石宏观破裂面的产生。从模拟结果可知,岩石发生破坏时其理想声发射空间定位结果应主要聚集在岩石宏观破坏面附近；相比而言,尽管声发射试验受多方面因素影响,如试验环境、岩石非均质结构等,定位的精确性及规律性并不是很完美,但其结果仍可提高对岩石破裂机理的了解。
![]()
图 6 模拟与试验累计声发射规律Figure 6. Accumulative acoustic emission properties of simulation and experiment3.3 晶体尺度裂纹演化规律
图7为岩石加载过程中晶体尺度裂纹演化规律。当轴向应力超过起裂应力后,岩石内部裂纹起裂主要以矿物晶间裂纹为主,且破坏模式大多为拉伸型破坏,这与Kranz[33]发现材料拉伸强度普遍小于剪切强度的结论所一致。随轴向加载的进行,晶内拉伸裂纹数量将超过晶间拉伸裂纹并占主导地位。当试样发生破坏时,晶内裂纹与总裂纹数的占比约为60.95%,拉伸裂纹与总裂纹数的占比约为93.87%。
为验证上述模型模拟所得晶体尺度裂纹发展规律,将试验后花岗岩试样分别进行偏光显微镜及电子显微镜观察。通过在试样表面多处进行观察发现,花岗岩试样受力破坏后所产生的宏观裂纹在矿物晶体尺度主要以不同矿物内部裂纹为主,少部分为矿物晶间裂纹,且该试验结果与文献[34]所得结论一致,故GBM模型可较真实地模拟岩石外力作用下晶体尺度裂纹演化规律。花岗岩代表性微裂纹分析结果如图8所示,图中Qz为石英矿物,Bt为云母矿物,Pl为斜长石,Kf为钾长石。
图9为岩石加载过程中矿物晶体内裂纹演化规律。加载初期,长石矿物晶体内裂纹首先发生起裂,随后石英、云母等矿物内裂纹依次开始起裂并扩展；当试样发生破坏时,矿物内裂纹数量关系如下：长石内裂纹>云母内裂纹>石英内裂纹>其他矿物内裂纹。由于长石矿物在此模型中的含量较高且力学特性相对石英矿物较差,故其内部裂纹最先发生起裂且数量始终占主导地位；相比云母矿物与石英矿物,尽管两者含量相近,但云母矿物相比石英矿物具有较差的力学性质,故岩石发生破坏时,云母内裂纹数量多于石英内裂纹。
在此,为进一步研究矿物晶体裂纹空间演化特征,图9中各加载点（A,B,C,D点）处晶体尺度裂纹分布图如图10所示。由图可知,当轴向应力达到A点时,模型中晶体尺度裂纹主要为矿物晶间裂纹,且裂纹随机分布在模型中；当轴向应力达到B点时,晶体尺度裂纹开始发生聚集现象,此时宏观破裂面开始形成。随着微裂纹的聚集程度加深,岩石宏观破裂面的发展及贯通将导致模型发生破坏。另外,观察宏观破裂面处晶体内裂纹分布发现,长石及云母矿物内裂纹占主导地位,石英矿物内裂纹次之,故在此可推断岩石发生破坏时微裂纹扩展及聚合主要以矿物晶体内裂纹为主,尤其长石及云母矿物内裂纹。
3.4 不同长石矿物含量岩石力学特性
在上述模型基础上,为研究不同矿物组成岩石宏观力学特性,本文在控制石英及其他矿物成分含量不变的条件下,通过Neper软件改变岩石中长石含量w以10%的梯度从20%变化至60%,生成5组不同的GBM模型,且每组模型进行3次共计15次数值模拟。模型如图11所示,黑色多边形为云母矿物,灰色多边形为长石矿物,白色多边形为石英矿物,蓝色多边形为其他矿物。
图12为不同长石含量岩石峰值应力及特征应力变化规律。由图12（a）可知,模拟所得峰值应力整体随长石矿物含量增加呈增加的变化趋势；当长石含量从20%增加至60%时,模拟所得峰值应力从168 MPa变化至178 MPa,但其变化幅度与文献[2]中统计拟合结果相比较小,这可能与模拟及文献[2]中岩石微观结构差异较大具有一定的联系。另外,随长石矿物含量增加,岩石损伤应力呈轻微增加的变化趋势,而起裂应力受长石矿物含量影响较小（图12（b））。
![]()
图 12 不同长石含量模型峰值应力及特征应力变化曲线Figure 12. Curves of peak and characteristic stresses of models with different feldspar contentsÜndül[9]通过试验指出,斜长石矿物与矿物基质含量占比的增加将增大单轴压缩强度,而矿物含量与起裂应力之间并没有较明确的关系,故此模拟所得结果与试验结果相吻合。
现从微观角度对不同长石含量岩石单轴压缩条件下晶体尺度裂纹变化规律进行分析,以进一步研究并讨论峰值应力变化的原因。将上述模型破坏时矿物晶体内裂纹数量取平均值并与总裂纹数进行正则化处理后,得不同长石含量GBM模型矿物晶体内裂纹占比变化规律如图13。
![]()
图 13 不同长石含量模型矿物内裂纹变化规律Figure 13. Variations of intragranular cracks of models with different feldspar contents由图13可知,随长石矿物含量增加,岩石云母矿物内裂纹占比呈递减的变化趋势,而长石矿物内裂纹占比呈递增的变化趋势,石英及其他矿物内裂纹占比则变化不大,故长石矿物含量的升高将导致岩石内部裂纹形式从以云母内裂纹为主转变为以长石内裂纹为主,模型破坏时相应的晶体尺度裂纹分布见图14。随长石矿物含量增加,岩石宏观破裂面处云母内裂纹数量减少,而长石内裂纹数量增加；相比云母矿物而言,长石矿物由于具有更好的力学性质且内部裂纹发生扩展所需要的能量越多,故随着长石含量的增加,岩石内部裂纹形式从以云母内裂纹为主转变为以长石内裂纹为主可能导致岩石峰值应力整体呈增加变化趋势。
![]()
图 14 不同长石含量模型晶体尺度裂纹分布图Figure 14. Distributions of grain-scale cracks of models with different feldspar contents非均质岩石力学特性的数值方法。4. 结论
本文结合花岗岩室内单轴试验并采用矿物晶体模型（GBM）对岩石矿物非均质结构进行了离散元模拟,改变岩石矿物组成研究了岩石微观非均质性对宏观力学特性的影响,并从矿物晶体尺度对其变化原因进行了较为合理的解释。主要结论如下：
（1）矿物晶体模型（GBM）可在考虑岩石矿物晶体结构及强度非均质性的前提下,对其应力-应变曲线、声发射特性进行有效模拟,并可研究矿物晶体尺度微裂纹起裂、扩展及聚合规律,是一种有效研究
（2）岩石单轴压缩过程中,微裂纹起裂首先以晶间裂纹为主,随后以晶内裂纹为主,且主要破坏机理为拉伸型破坏。岩石破坏时,拉伸裂纹及晶内裂纹与总裂纹的占比分别约为93.87%,60.95%。
（3）岩石单轴压缩过程中,长石矿物内裂纹首先起裂,随后云母、石英等矿物内裂纹依次发生起裂并扩展。岩石破坏时,微裂纹聚集将导致宏观破裂面的形成,且宏观破裂面的形成主要与晶体内裂纹扩展及聚合有关。
（4）随长石矿物含量增加,岩石峰值应力及损伤应力整体呈增加的变化趋势,而起裂应力变化并不明显。从矿物晶体尺度裂纹变化规律可推断,岩石峰值应力增大的原因可能与长石内裂纹占比增加、云母内裂纹占比减少有关。
-
![]()
图 1 1920年海原大地震地震滑坡和灾害分布图（据Close等[4], 本文清绘）
Figure 1. Distribution of landslides and disasters induced by Haiyuan Ms 8.5 earthquake in 1920 (Upton Close, 1922, clearly drawn by the authors)
![]()
图 2 1920年海原MS8.5级地震引发的大规模黄土滑坡
Figure 2. Large-scale loess landslides induced by Haiyuan MS 8.5 earthquake in 1920
![]()
图 4 1920年海原8.5地震静宁孙家沟黄土层内剪切型滑坡
Figure 4. Shear landslide within loess strata at Sunjiagou in Jingning county induced by Haiyuan MS 8.5 earthquake
![]()
![]()
图 7 2013年岷县-漳县MS 6.6级地震永光村西底-表联合液化型黄土流滑
Figure 7. Mud flow and landslide at west Yongguang village induced by liquefaction in both bottom and surface layers during the Minxian-Zhangxian MS 6.6 earthquake in 2013
![]()
图 10 剪切型滑坡振动台试验宏观破坏特征
Figure 10. Macroscopic failure characteristics of shear landslidebased on shaking table tests
![]()
图 11 黄土-基岩接触面剪切型滑坡振动台试验模型与传感器布设
Figure 11. Shaking table test model and layout of sensors for loess-bedrock interface shear landslide
![]()
图 12 黄土-基岩接触面剪切滑坡不同位置加速度放大效应
Figure 12. Acceleration amplification effects of shear landslide on loess-bedrock interface at different locations
![]()
图 13 黄土-基岩滑坡破坏时接触面和上覆黄土土压力分布
Figure 13. Soil pressure distribution of loess-bedrock interface and overlying loess during shear landslide failure
![]()
图 14 黄土-基岩接触面滑坡振动台试验宏观破坏特征
Figure 14. Macroscopic failure characteristics of loess-bedrock shear landslide based on shaking table tests
![]()
图 15 底部液化型黄土滑移振动台试验模型与传感器布设
Figure 15. Shaking-table test model for liquefaction at bottom-triggered sliding of loess deposit and layout of sensors
![]()
图 16 液化型黄土滑坡饱和土层孔隙水压力增长趋势
Figure 16. Increasing trend of pore water pressure in saturated soil of liquefaction loess landslide base on shaking table tests
![]()
图 17 底部饱和液化型黄土滑坡振动台试验宏观破坏特征
Figure 17. Failure characteristics of liquefaction at bottom-triggered sliding of loess deposit
![]()
图 18 表层液化流滑型滑坡试验模型尺寸与传感器布设
Figure 18. Sizes of shaking-table test model and layout of sensors for surface liquefaction-triggered sliding flow of loess slope
![]()
图 19 表层饱和液化泥流型黄土斜坡流滑振动台试验宏观破坏特征
Figure 19. Macroscopic failure characteristics of liquefaction- triggered mud flow of loess slope with saturated surface based on shaking table tests
![]()
图 20 坡角至坡顶不同部位急剧增长的孔隙水压力
Figure 20. Dramatic increase of pore water pressures at different locations from slope foot to slope crest
![]()
图 21 黄土震陷性滑坡振动台试验模型与传感器布设
Figure 21. Shaking-table test model and layout of sensors for loess landslides due to seismic subsidence
![]()
图 22 震陷型黄土滑坡振动台试验宏观破坏特征
Figure 22. Macroscopic failure characteristics of seismic subsidence- triggered loess landslide in shaking table tests
![]()
图 23 基于应变观测的黄土斜坡震陷变形失稳过程
Figure 23. Seismic subsidence-triggered deformation and instability process of loess slope based on strain observation
表 1 黄土地区9次大震和强震诱发的滑坡灾害
Table 1 Landslides induced by 9 great and strong earthquakes in Loess Plateau
地震名称 震级 震中烈度 断裂性质 遇难人数 地震岩土灾害 1303年洪洞地震 8.0 Ⅺ 右旋走滑 20万 滑坡、液化滑移 1556年华县地震 Ⅺ 正断 83万 密集滑坡、震陷 1654年天水南地震 8.0 Ⅺ 左旋走滑兼正断 3.1万 密集滑坡 1695年临汾地震 Ⅹ 正断 5.3万 滑坡、液化滑移 1718年通渭地震 7.5 Ⅹ 逆冲 4万 密集滑坡 1879年武都地震 8.0 XI 逆冲兼右旋走滑 2.2万 滑坡、崩塌 1920年海原地震 8.5 Ⅻ 左旋走滑 27万 密集滑坡、液化滑移 1927年古浪地震 8.0 Ⅺ 西段逆冲，东段正断 4万 密集滑坡 2013年岷—漳地震 6.6 Ⅷ 逆冲为主，兼具左旋走滑 12 滑坡、液化泥流 
下载: 导出CSV
表 2 底部饱和液化滑移型模型地层参数
Table 2 Stratum parameters of bottom saturated liquefaction slip model
地层土样 取样深度/m 干密度/(g·cm-3) 含水率w/% 孔隙比 饱和度Sr /% 颗粒组成/% 黏粒 粉粒 砂粒 上部非饱和黄土层（Q3、Q4） 6.5 1.35 5.51 1.12 13.36 13.15 77.99 8.86 第一古土壤层（Q3） 13.5 1.68 11.14 0.79 38.08 16.97 81.01 2.02 饱和砂质黄土层（Q3） 16.5 1.65 27.69 1.09 75.40 12.20 59.66 28.14 第二古土壤层（Q2） 35.5 1.72 26.82 1.00 72.82 13.22 66.21 20.57
下载: 导出CSV
-
[1] 王兰民, 蒲小武, 陈金昌. 黄土高原地震诱发滑坡分布特征与灾害风险[J]. 城市与减灾, 2019(3): 33-40. https://www.cnki.com.cn/Article/CJFDTOTAL-CSJZ201903009.htm WANG Lanmin, PU Xiaowu, CHEN Jinchang. Distribution characteristics and disaster risk of earthquake-induced landslides in loess plateau[J]. City and Disaster Reduction, 2019(3): 33-40. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-CSJZ201903009.htm
[2] 赵晋泉, 张大卫, 高树义, 等. 1303年山西洪洞8级大地震郇堡地滑之研究[J]. 山西地震, 2003(3): 16-22. https://www.cnki.com.cn/Article/CJFDTOTAL-SXDZ200303006.htm ZHAO Jinquan, ZHANG Dawei, GAO Shuyi, et al. Huanbu ground slide, the relic of 1303 Hongtong, Shanxi, earthquake of M8[J]. Earthquake Research in Shanxi, 2003(3): 16-22. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SXDZ200303006.htm
[3] 陈国顺, 王国强. 临汾1695年八级地震及震后地表垂直运动[J]. 山西地震, 1984(4): 31-36. https://www.cnki.com.cn/Article/CJFDTOTAL-SXDZ198404007.htm CHEN Guoshun, WANG Guoqiang. The vertical movement of ground surface after the 1695 M8 earthquake in Linfen[J]. Earthquake Research in Shanxi, 1984(4): 31-36. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SXDZ198404007.htm
[4] CLOSE U, MCCORMICK E. Where the mountains walked[J]. The National Geographic Magazine, 1922, LI(5): 443-464.
[5] 张晓超, 裴向军, 张茂省, 等. 强震触发黄土滑坡流滑机理的试验研究: 以宁夏党家岔滑坡为例[J]. 工程地质学报, 2018, 26(5): 1219-1226. https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ201805013.htm ZHANG Xiaochao, PEI Xiangjun, ZHANG Maosheng, et al. Experimental study on mechanism of flow slide of loess landslides triggered by strong earthquake—a case study in Dangjiacha, Ningxia Province[J]. Journal of Engineering Geology, 2018, 26(5): 1219-1226. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ201805013.htm
[6] 王家鼎, 张倬元. 典型高速黄土滑坡群的系统工程地质研究[M]. 成都: 四川科学技术出版社, 1999. WANG Jiading, ZHANG Zhuoyuan. Systematic Engineering Geological Study on Typical High-Speed Loess Landslide Group[M]. Chengdu: Sichuan Scientific & Technical Publishers, 1999. (in Chinese)
[7] 王家鼎, 张倬元. 地震诱发高速黄土滑坡的机理研究[J]. 岩土工程学报, 1999, 21(6): 670-674. http://www.cgejournal.com/cn/article/id/10420 WANG Jiading, ZHANG Zhuoyuan. A study on the mechanism of high speed loess landslide induced by earthquake[J]. Chinese Journal of Geotechnical Engineering, 1999, 21(6): 670-674. (in Chinese) http://www.cgejournal.com/cn/article/id/10420
[8] WANG G H, SASSA K. Pore-pressure generation and movement of rainfall-induced landslides: effects of grain size and fine-particle content[J]. Engineering Geology, 2003, 69(1/2): 109-125.
[9] 张振中. 黄土地震灾害预测[M]. 北京: 地震出版社, 1999. ZHANG Zhenzhong. Earthquake Disaster Prediction of Loess[M]. Beijing: Seismological Press, 1999. (in Chinese)
[10] 卢育霞, 石玉成, 陈永明, 等. 地震诱发黄土滑坡的滑距估测[J]. 西北地震学报, 2006(3): 248-251. https://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ200603009.htm LU Yuxia, SHI Yucheng, CHEN Yongming, et al. Slippage estimation of the loess landslide triggered by earthquake[J]. Northwestern Seismological Journal, 2006(3): 248-251. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ200603009.htm
[11] 王鼐, 王兰民, 王谦, 等. 黄土高原地震作用下黄土滑坡滑距预测方法[J]. 地震工程学报, 2016, 38(4): 533-540. https://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ201604007.htm WANG Nai, WANG Lanmin, WANG Qian, et al. Forecasting method for sliding distance of seismic landslides on the loess plateau, China[J]. China Earthquake Engineering Journal, 2016, 38(4): 533-540. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ201604007.htm
[12] 王兰民. 黄土地层大规模地震液化滑移的机理与风险评估[J]. 岩土工程学报, 2020, 42(1): 1-19. doi: 10.11779/CJGE202001001 WANG Lanmin. Mechanism and risk evaluation of sliding flow triggered by liquefaction of loess deposit during earthquakes[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(1): 1-19. (in Chinese) doi: 10.11779/CJGE202001001
[13] 国家地震局兰州地震研究所, 宁夏回族自治区地震队. 一九二○年海原大地震[M]. 北京: 地震出版社, 1980. Lanzhou Institute of Seismology, State Seismological Bureau, Seismological Team of Ningxia Hui Autonomous Region. Haiyuan Earthquake in 1920[M]. Beijing: Seismological Press, 1980. (in Chinese)
[14] 王兰民. 黄土动力学[M]. 北京: 地震出版社, 2003. WANG Lanmin. Loess Dynamics[M]. Beijing: Seismological Press, 2003. (in Chinese)
[15] 白铭学, 王增光. 宁夏黄土高原区低角度滑移研究报告[R]. 宁夏: 宁夏回族自治区地震局, 1987. BAI Mingxue, WANG Zengguang. Study on Low Angle Slip in Loess Plateau of Ningxia[R]. Ningxia: Earthquake Agency of Ningxia Hui Autonomous Region, 1987. (in Chinese)
[16] 王兰民, 吴志坚. 岷县漳县6.6级地震震害特征及其启示[J]. 地震工程学报, 2013, 35(3): 401-412. https://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ201303002.htm WANG Lanmin, WU Zhijian. Earthquake damage characteristics of the Minxian-Zhangxian MS6.6 earthquake and its lessons[J]. China Earthquake Engineering Journal, 2013, 35(3): 401-412. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ201303002.htm
-
期刊类型引用(3)
1. 郑刚，甄洁，程雪松，邓旭，宋许根. 基坑坑底抗隆起稳定性研究现状与展望. 建筑结构. 2024(20): 129-142 . 
百度学术
2. 肖伟，路轶，杨严禄，肖桃李，郝勇. 考虑尺寸效应的狭窄型基坑稳定性分析. 科学技术与工程. 2022(10): 4105-4112 . 百度学术
3. 李林，黄奕鹏，苏宁川. 狭窄基坑抗隆起稳定性计算方法的改进研究. 地下空间与工程学报. 2022(S1): 74-81 . 百度学术
其他类型引用(12)
-
其他相关附件
计量
- 文章访问数: 785
- HTML全文浏览量: 77
- PDF下载量: 287
- 被引次数: 15
