地震作用下陈庄黄土-泥岩接...面滑坡失稳机理及稳定性研究_马紫娟.pdf

第 卷第期 年月地震工程学报 ，收稿日期：基金项目：甘肃省地震局地震科技发展基金项目（，）第一作者简介：马紫娟（），女，硕士，高级工程师，主要从事地震地质灾害方面的研究。：。通信作者：张有龙（），男，硕士，高级工程师，主要从事地震地质和工程地震方面的研究和应用。：。马紫娟，张有龙，刘小丰，等 地震作用下陈庄黄土泥岩接触面滑坡失稳机理及稳定性研究 地震工程学报，（）：，（）：地震作用下陈庄黄土泥岩接触面滑坡失稳机理及稳定性研究马紫娟，张有龙，刘小丰，郑龙，高中南（甘肃省地震局，甘肃 兰州 ）摘要：以甘肃礼县陈庄黄土泥岩接触面滑坡为研究对象，在对该滑坡的地震地质环境、基本特征等调查的基础上，运用数值分析对陈庄黄土泥岩接触面滑坡在各人工计算场地基岩地震动时程作用下的失稳机理及稳定性进行动态模拟研究。结果表明，陈庄滑坡在地震诱发下变形破坏模式为受拉和受剪的复合破坏模式。失稳机理过程可概括为：地震作用底部岩体受张拉潜在滑面主滑段上部剪应力增大剪应变集中带出现剪应变集中带延伸扩展贯通滑坡失稳；滑坡稳定性方面，随着输入地震动的增大，陈庄黄土泥岩接触面滑坡稳定性降低，其中在 年超越概率的基岩地震动时程作用下滑坡处于不稳定状态。关键词：黄土泥岩接触面滑坡；陈庄滑坡；失稳机理；稳定性中图分类号：文献标志码：文章编号：（）：，（，）：，：，：；引言黄土新近纪红色泥岩接触面滑坡是黄土地区分布最广、发生最频繁的滑坡类型之一。此类滑坡由于“双层异质”的结构，在地震等外因作用下破坏力强，灾害损失惨重，影响力巨大。目前尽管众多学者对黄土滑坡机理做了大量的研究，但针对黄土泥岩接触面滑坡失稳机理的研究还较少：文宝萍等提出了牵引、累进式变形和沿接触面的解体式块体滑移并伴随沿古土壤层的局部滑移的失稳方式；李媛等提出天水市孟家山滑坡为蠕动拉裂滑带扩展剪出口形成突滑的变形破坏机理；石瑞红应用工程地质成因法对黄土泥岩（红层）接触面滑坡的变形特征和机理进行了探讨；郑程波分析了坡角、降雨以及人为活动引起的黄土泥岩接触面滑移机理；马健应用 分析了池沟滑坡稳定性，探讨了其形成演化机理；朱建东等分析了两种典型降雨模式下黄土泥岩接触面滑坡的失稳模式等。上述研究大都是在降雨、开挖等诱因下的机理，地震作用下的研究目前主要集中在黄土震陷和液化失稳机理方面：如陈永明等运用黄土震陷对兰州永登地震产生的黄土滑坡机理进行了解释；石玉成等通过黄土的微结构来研究震陷机理，解释黄土地震滑坡的形成机理；王兰民用土动力学方法探讨了黄土层的液化机理。但应用数值模拟动态来研究黄土泥岩接触面滑坡在地震作用下的失稳机理却不多。为此本文以甘肃礼县陈庄黄土泥岩接触面滑坡为例，对其地震地质条件、地形地貌和地层结构等进行调查，在此基础上对该滑坡在地震作用下的逐步变形失稳过程进行动态数值模拟，分析其在不同人工计算场地基岩地震动时程工况下的瞬时动态应力、变形和稳定系数等参数，从而展开对黄土泥岩接触面滑坡在地震作用下的失稳机理和稳定性研究。滑坡体地震地质条件陈庄滑坡所处礼县位于青藏高原东缘西礼盆地东部，即西汉水及其支流河谷地带，大地构造单元为青藏高原东北缘西秦岭北部褶皱带。新构造运动强烈，境内及周边 范围内发育了规模不等，性质各异的 条断裂，其中洮坪盐官断裂为礼县罗家堡断裂的一部分，是晚更新世全新世活动断裂，该断裂带上地震时有发生，震级可达级，其中 年天水南级地震就发生在该断裂带上。目前研究区及外围 范围内共发生中强地震 次，其中级地震 次，级地震次，级地震次，这些地震的烈度几乎都大于度。研究认为，产生滑坡的地震烈度黄土区约为度，石质山区度，可见在礼县能产生这类滑坡的地震是比较常见的。滑坡基本特征滑坡概况陈庄滑坡（图）位于礼县城关镇石碑村陈庄组与土山村之间。该滑坡所处位置宏观地貌为中山地貌，微观地貌为陡坡，滑坡位于斜坡的下部，河流的左岸，为顺向坡，形状为凹形。图陈庄滑坡全貌 该滑 坡 属 特 大 型 滑 坡，长 约 ，宽 约 ，厚约，体积约 ，前后缘相对高差约 。滑坡主滑方向为 ，坡度 。滑坡体内发育了多个次级滑坡和潜在滑坡，其中沿主滑方向分布了个滑坡。滑坡的形貌特征（）平面形态该滑坡整体呈似簸箕形的不规则平面形态，滑第 卷 第期马紫娟，等：地震作用下陈庄黄土泥岩接触面滑坡失稳机理及稳定性研究坡侧壁以东西部山脊为界，两侧山脊向北延伸并与北侧的陡坡相连形成平均高差为 （最大不足）的黄土陡壁为滑坡后缘，前缘以西汉水河谷西北侧田地村庄为界。（）后壁和侧壁该滑坡后壁大致呈长弧形圈椅状。除后部突出山脊处的后壁略缓外，整体较陡，后壁倾角约为 ，局部近直立，高 （图）。后壁整体为黄土覆盖，底部局部出露红层。此外后壁表面有较明显的滑动擦痕，局部有规模较小的崩塌和裂缝（图）。图滑坡后壁特征 图滑坡后壁的裂缝和局部崩塌 滑坡两侧壁处分布着经过长时间风化等作用后不明显的雁形排列羽状张裂缝。（）前部及剪出口滑坡前部剪出口已被农耕修整，黄土堆积物覆盖未找到明显的剪出口，滑坡舌呈小范围“平原”向谷地外延伸的形态，滑坡前缘后部局部可见由次级滑坡滑动造成的陡坎（图）。滑坡微地貌特征（）平台及陡坎滑坡后壁前存在平台，未被夷为平地，靠近侧壁和后壁处较陡。平台上覆黄土疏松，节理裂隙较发育。在滑坡后缘和次级滑坡的后壁上发育有高约 的陡坎，倾角 ，物质组成为黄土和红层，其中后部有红层出露（图）。图滑坡前缘特征 图次级滑坡后壁的陡坎 （）冲沟坡体中冲沟发育，长 不等，形状为“”型（图），沟口处于滑坡前部，深度随两侧地形不等，两侧有小崩塌，沟内黄土层垂直节理和裂缝发育。局部沟底黄土与棕红色泥岩接触面清晰，为厚约 的褐红色滑带土。部分沟中有径流冲刷侵蚀而成的黄土桥连接冲沟两侧。图坡体上的冲沟 （）落水洞与裂缝落水洞与裂缝主要分布在后壁下方、平台后部、陡坎前部和冲沟尾部。落水洞呈圆形，大小深度各不相等，直径，深度（图）。此外坡体内发育多条纵向裂缝，排列位置大致呈放射状，贯通性好，长度在 间，张开度 地震工程学报 年间，裂缝内均有少量风化土充填。图坡体上的落水洞和裂缝 滑坡地层结构该滑坡上覆第四系风成黄土，下伏第三系红色泥岩（图）。图陈庄滑坡剖面图 滑体为第四系风成黄土，结构松散零乱，滑体中碎石含量约，碎石大小为 （长）（宽）（高）。垂直节理发育，承载力较差，属低承载力、高压缩性土。滑面为弧形，后缘较陡，埋深约。滑带土属软岩滑带土，为粉质黏土，其中黏粒含量较高，厚度约，湿，呈软塑状。滑床为第三系泥岩，局部为砂质泥岩，基岩产状为 ，棕红色，表面风化破碎，风化带厚度约，具层理，成岩作用较弱，属极软岩，在原状结构条件下，力学性质较差。滑坡失稳机理及稳定性分析模型建立地震作用考虑为在基岩底部输入水平地震加速度时程（不考虑竖向地震力），并认为在地震到达时地震加速度同时作用在模型的每一个单元节点上。模拟时首先采用线弹性土体模型进行静力计算获得初始应力，将其作为动力分析的初始条件，采用等效线性材料模式进行有限元动力分析，得到不同时刻地震力作用下的应力、应变；同时采用 法和 永久位移法 得出滑坡沿指定滑面的安全系数和永久位移时程曲线。选取陈庄滑坡典型剖面（图）作为研究对象，综合考虑计算速度与精度要求，共划分 个有限元网格单元，个节点。边界条件采用位移约束：静力计算时，左、右边界约束水平向位移，底部双向位移约束；动力计算时，左、右边界为竖向约束，底部边界保持双向约束不变。计算参数地层力学参数方面，通过室内试验和工程类比，该滑坡体主要组成地层的力学参数指标列于表。表滑坡体岩土层物理力学参数 岩性天然重度（）黏聚强度 摩擦角（）弹性模量 初始剪切模量 泊松比阻尼比滑体风成黄土 滑带土粉质黏土 滑床泥岩 地震参数设置方面，本次模拟选取 礼县县城及盐官镇地震小区划项目报告所提供的礼县县城 年超越概率、的人工计算场地基岩地震动时程分别作为输入地震动（简称为工况一、工况 二、工 况 三），其 地 震 动 加 速 度 峰 值 分 别 为 、.、.。在保证计算精度的情况下，为方便求解，本次模拟选取 的数据，定义时间步长为。数值模拟采用的地震波时程曲线见图。结果分析（）应力和变形分析以工况三为例研究滑坡在地震作用下的应力和变形特征，将滑坡在天然自重作用下的变形归零后，施加地震动加速度时程曲线。地震作用后，滑坡体应力最大值和整个应力场第 卷 第期马紫娟，等：地震作用下陈庄黄土泥岩接触面滑坡失稳机理及稳定性研究礼县县城及盐官镇地震小区划项目报告兰州：甘肃省地震工程研究院，（图）均发生了变化，滑坡体应力较震前有明显的增加，并出现了应力集中带。坡面上，由于底层岩体受到张拉力的作用，沿着坡面一定深度范围内的应力值减小。潜在滑面上，由于主滑段上向拉应力集中带导致其上部剪应力有明显的增加。最大剪应变方面（图），在地震作用下，坡体产生了剪切变形。随着地震动的作用，沿滑带中上部出现剪应变集中带，相对地震前应变值明显增大，图地震波时程曲线 图 滑坡在工况三下的应力分布云图（单位：）（：）图 滑坡体在工况三各时段最大剪应变云图 贯通率也显著增大，地震初期（、）时伴随着剪应变值在潜在滑带和滑体内增大的同时，剪应变集中带 分 布 范 围 沿 着 主 滑 段 延 伸，地 震 作 用 后 期（、）剪应变值继续在滑面和滑体内增大；主滑段最大剪应变集中带继续扩展，同时下部抗滑段上也出现了剪应变集中带，最后基本贯通了整个滑面。从位移分布（图 ）可以看出，滑坡的位移地震工程学报 年最大值为 ，水平向位移分量最大值为 ，出现在坡面中后部地形陡缓过渡处，即次级滑面和前缘潜在滑面后缘处。此外滑坡的后缘位移较前缘大，这与滑坡在变形破坏中后缘易出现拉裂缝变形现象相吻合。下覆基岩整体稳定，位移变化不大，上覆滑体位移较下伏基岩位移大一些，潜在滑面上主滑段和牵引段基本都出现了一定的位移，可见该滑坡会发生整体滑动。图 滑坡在工况三下的位移云图（）（）图 滑坡在工况三下的向位移云图（）（）（）稳定性评价陈庄滑坡在地震动作用下沿现场勘探确定的推测主滑面和前缘潜在滑面（图）的安全系数时程曲线和 永久位移时程曲线见图 。滑坡在工况一：年超越概率，加速度峰值 为 地震动作用下（图），沿主滑面的图 滑坡在工况一下的稳定性参数时程曲线 安全系数最小值为 ，发生在时，永久位移为；沿前缘潜在滑面的安全系数最小值为 ，发生在 时，永久位移为。可见滑坡在工况一下虽按最小安全系数判定处于基本稳定状态，但沿主滑面和前缘潜在滑面都没有滑动，是稳定的。滑坡在工况二：年超越概率，加速度峰第 卷 第期马紫娟，等：地震作用下陈庄黄土泥岩接触面滑坡失稳机理及稳定性研究值为 的地震动荷载作用下（图），沿主滑面的安全系数最小值为 ，发生在 时，永久位移为；沿前缘潜在滑面的安全系数最小值为 ，发 生 在 时，永 久 位 移 为 。可见该滑坡在工况二下沿主滑面按最小安全系数判定处于欠稳定状态，没发生位移；沿前缘潜在滑面按最小安全系数判定处于不稳定状态，发生了很小的位移。图 滑坡在工况二下的稳定性参数时程曲线 滑坡在工况三：年超越概率，加速度峰值为 的地震动荷载作用下（图），沿主滑面的安全系数最小值为 ，发生在 时，永久位移为 ；沿前缘潜在滑面的安全系数最小值为 ，发生在时，永久位移为 。可见滑坡在工况三下沿主滑面和前缘滑面都处于不稳定状态，发生了滑动。在地震荷载作用下，选用整个地震历时中的最小安全系数评价滑坡的稳定性没有考虑地震荷载作用的过程和边坡破坏的积累性是不合理的。本文采用文献 提出的最小平均安全系数 式（）作为评价指标：（）（）式中：为边坡在静力作用下的安全系数；为地震作用下安全系数随震动过程而波动的最小值。采用上述方法沿指定滑动面下所计算的最小平均安全系数列于表，与上述 永久位移结合最小安全系数判定的稳定性结论基本一致。同时随着输入地震动的增大，最小平均安全系数逐渐减小，滑坡稳定性降低，其中在 年超越概率的基岩地震动时程作用下陈庄滑坡处于不稳定状态。结论本文以礼县陈庄黄土泥岩接触面滑坡为例，在对其发育特征调查的基础上，对其在地震作用下的动力响应进行了数值分析，从而对地震作用下黄土泥岩接触