地震作用下陳莊黃土－泥巖接觸面滑坡失穩機理及穩定性研究

馬紫娟　張有龍　劉小豐　鄭龍　高中南

摘要： 以甘肅禮縣陳莊黃土－泥巖接觸面滑坡為研究對象，在對該滑坡的地震地質環境、基本特征等調查的基礎上，運用數值分析對陳莊黃土－泥巖接觸面滑坡在各人工計算場地基巖地震動時程作用下的失穩機理及穩定性進行動態模擬研究。結果表明，陳莊滑坡在地震誘發下變形破壞模式為受拉和受剪的復合破壞模式。失穩機理過程可概括為：地震作用→底部巖體受張拉→潛在滑面主滑段上部剪應力增大→剪應變集中帶出現→剪應變集中帶延伸擴展貫通→滑坡失穩;滑坡穩定性方面，隨著輸入地震動的增大，陳莊黃土－泥巖接觸面滑坡穩定性降低，其中在50年超越概率2%的基巖地震動時程作用下滑坡處于不穩定狀態。

關鍵詞： 黃土－泥巖接觸面滑坡; 陳莊滑坡; 失穩機理; 穩定性

中圖分類號： P642.22文獻標志碼：A 文章編號： 1000－0844（2023）04-0810-09

DOI：10.20000/j.1000－0844.20201009001

Instability mechanism and stability analysis of Chenzhuangloess－mudstone interface landslide under seismic action

MA Zijuan， ZHANG Youlong， LIU Xiaofeng， ZHENG Long， GAO Zhongnan

Abstract：? The loess－mudstone interface landslide in Chenzhuang Village， Lixian County， Gansu Province， was investigated in this study. A dynamic numerical simulation of the instability mechanism and stability of the landslide subjected to ground motion at various artificial calculation sites was conducted based on an analysis of the seismic geological environment and basic characteristics of the landslide. Results show that the deformation and failure of the Chenzhuang landslide are associated with a composite failure mode of tension and shear when subjected to ground motion. The process of instability can be summarized as follows： ground motion due to an earthquake → tension at the bottom of rock mass → increase in shear stress at the upper part of the main slip section of the potential slip surface → emergence of a shear strain concentration zone → extension and penetration of the shear strain concentration zone → landslide instability. The stability of the Chenzhuang loess－mudstone interface landslide decreases with an increase in the ground motion， and the landslide is in an unstable state when subjected to bedrock ground motion with a exceeding probability of 2% in 50 years.

Keywords： loess－mudstone interface landslide; Chenzhuang landslide; instability mechanism; stability

0 引言

黃土－新近紀紅色泥巖接觸面滑坡是黃土地區分布最廣、發生最頻繁的滑坡類型之一［1］。此類滑坡由于 “雙層異質”的結構，在地震等外因作用下破壞力強，災害損失慘重，影響力巨大。目前盡管眾多學者對黃土滑坡機理做了大量的研究，但針對黃土－泥巖接觸面滑坡失穩機理的研究還較少：文寶萍等［1］提出了牽引、累進式變形和沿接觸面的解體式塊體滑移并伴隨沿古土壤層的局部滑移的失穩方式;李媛等［2］提出天水市孟家山滑坡為蠕動－拉裂－滑帶擴展－剪出口形成－突滑的變形破壞機理;石瑞紅［3］應用工程地質成因法對黃土－泥巖（紅層）接觸面滑坡的變形特征和機理進行了探討;鄭程波［4］分析了坡角、降雨以及人為活動引起的黃土泥巖接觸面滑移機理;馬健［5］應用Geostudio分析了池溝滑坡穩定性，探討了其形成演化機理;朱建東等［6］分析了兩種典型降雨模式下黃土－泥巖接觸面滑坡的失穩模式等。上述研究大都是在降雨、開挖等誘因下的機理，地震作用下的研究目前主要集中在黃土震陷和液化失穩機理方面：如陳永明等［7］運用黃土震陷對蘭州永登地震產生的黃土滑坡機理進行了解釋;石玉成等［8］通過黃土的微結構來研究震陷機理，解釋黃土地震滑坡的形成機理;王蘭民［9］用土動力學方法探討了黃土層的液化機理。但應用數值模擬動態來研究黃土－泥巖接觸面滑坡在地震作用下的失穩機理卻不多。

為此本文以甘肅禮縣陳莊黃土－泥巖接觸面滑坡為例，對其地震地質條件、地形地貌和地層結構等進行調查，在此基礎上對該滑坡在地震作用下的逐步變形失穩過程進行動態數值模擬，分析其在不同人工計算場地基巖地震動時程工況下的瞬時動態應力、變形和穩定系數等參數，從而展開對黃土－泥巖接觸面滑坡在地震作用下的失穩機理和穩定性研究。

1 滑坡體地震地質條件

陳莊滑坡所處禮縣位于青藏高原東緣西—禮盆地東部，即西漢水及其支流河谷地帶，大地構造單元為青藏高原東北緣西秦嶺北部褶皺帶。新構造運動強烈，境內及周邊25 km范圍內發育了規模不等，性質各異的15條斷裂，其中洮坪—鹽官斷裂為禮縣—羅家堡斷裂的一部分，是晚更新世－全新世活動斷裂，該斷裂帶上地震時有發生，震級可達6～8級，其中1654年天水南8級地震就發生在該斷裂帶上。目前研究區及外圍150 km范圍內共發生中強地震92次，其中6.0～6.9級地震20次，7.0～7.9級地震5次，8.0～8.9級地震2次，這些地震的烈度幾乎都大于Ⅵ度。研究認為，產生滑坡的地震烈度黃土區約為Ⅵ度，石質山區Ⅶ度［10］，可見在禮縣能產生這類滑坡的地震是比較常見的。

2 滑坡基本特征

2.1 滑坡概況

陳莊滑坡（圖1）位于禮縣城關鎮石碑村陳莊組與土山村之間。該滑坡所處位置宏觀地貌為中山地貌，微觀地貌為陡坡，滑坡位于斜坡的下部，河流的左岸，為順向坡，形狀為凹形。

該滑坡屬特大型滑坡，長約1 250 m，寬約1 300 m，厚約30 m，體積約4.875×107 m3，前后緣相對高差約444 m。滑坡主滑方向為155°，坡度30°。滑坡體內發育了多個次級滑坡和潛在滑坡，其中沿主滑方向分布了2個滑坡。

2.2 滑坡的形貌特征

（1） 平面形態

該滑坡整體呈似簸箕形的不規則平面形態，滑坡側壁以東西部山脊為界，兩側山脊向北延伸并與北側的陡坡相連形成平均高差為15 m（最大不足20 m）的黃土陡壁為滑坡后緣，前緣以西漢水河谷西北側田地村莊為界。

（2） 后壁和側壁

該滑坡后壁大致呈長弧形圈椅狀。除后部突出山脊處的后壁略緩外，整體較陡，后壁傾角約為60°～80°，局部近直立，高10～20 m（圖2）。后壁整體為黃土覆蓋，底部局部出露紅層。此外后壁表面有較明顯的滑動擦痕，局部有規模較小的崩塌和裂縫（圖3）。

滑坡兩側壁處分布著經過長時間風化等作用后不明顯的雁形排列羽狀張裂縫。

（3） 前部及剪出口

滑坡前部剪出口已被農耕修整，黃土堆積物覆蓋未找到明顯的剪出口，滑坡舌呈小范圍“平原”向谷地外延伸的形態，滑坡前緣后部局部可見由次級滑坡滑動造成的陡坎（圖4）。

2.3 滑坡微地貌特征

（1） 平臺及陡坎

滑坡后壁前存在平臺，未被夷為平地，靠近側壁和后壁處較陡。平臺上覆黃土疏松，節理裂隙較發育。

在滑坡后緣和次級滑坡的后壁上發育有高約5～15 m的陡坎，傾角50°～70°，物質組成為黃土和紅層，其中后部有紅層出露（圖5）。

（2） 沖溝

坡體中沖溝發育，長20～300 m不等，形狀為“V”型（圖6），溝口處于滑坡前部，深度隨兩側地形不等，兩側有小崩塌，溝內黃土層垂直節理和裂縫發育。局部溝底黃土與棕紅色泥巖接觸面清晰，為厚約10～25 cm的褐紅色滑帶土。部分溝中有徑流沖刷侵蝕而成的黃土橋連接沖溝兩側。

（3） 落水洞與裂縫

落水洞與裂縫主要分布在后壁下方、平臺后部、陡坎前部和沖溝尾部。落水洞呈圓形，大小深度各不相等，直徑0.5～1.5 m，深度1～3 m（圖7）。此外坡體內發育多條縱向裂縫，排列位置大致呈放射狀，貫通性好，長度在5～20 m間，張開度4～30 cm間，裂縫內均有少量風化土充填。

2.4 滑坡地層結構

該滑坡上覆第四系風成黃土，下伏第三系紅色泥巖（圖8）。

滑體為第四系風成黃土，結構松散零亂，滑體中碎石含量約5%，碎石大小為30 cm（長）×20 cm（寬）×10 cm（高）。垂直節理發育，承載力較差，屬低承載力、高壓縮性土。

滑面為弧形，后緣較陡，埋深約30 m。滑帶土屬軟巖滑帶土，為粉質黏土，其中黏粒含量較高，厚度約0.25 m，濕，呈軟塑狀。

滑床為第三系泥巖，局部為砂質泥巖，基巖產狀為190°∠5°，棕紅色，表面風化破碎，風化帶厚度約0.3 m，具層理，成巖作用較弱，屬極軟巖，在原狀結構條件下，力學性質較差。

3 滑坡失穩機理及穩定性分析

3.1 模型建立

地震作用考慮為在基巖底部輸入水平地震加速度時程（不考慮豎向地震力），并認為在地震到達時地震加速度同時作用在模型的每一個單元節點上。

模擬時首先采用線彈性土體模型進行靜力計算獲得初始應力，將其作為動力分析的初始條件，采用等效線性材料模式進行有限元動力分析，得到不同時刻地震力作用下的應力、應變;同時采用Morgenstern－Price法和Newmark永久位移法［11］得出滑坡沿指定滑面的安全系數和永久位移時程曲線。

選取陳莊滑坡典型剖面（圖8）作為研究對象，綜合考慮計算速度與精度要求，共劃分4 773個有限元網格單元，4 879個節點。邊界條件采用位移約束：靜力計算時，左、右邊界約束水平向位移，底部雙向位移約束;動力計算時，左、右邊界為豎向約束，底部邊界保持雙向約束不變。

3.2 計算參數

地層力學參數方面，通過室內試驗和工程類比，該滑坡體主要組成地層的力學參數指標列于表1。

地震參數設置方面，本次模擬選取《禮縣縣城及鹽官鎮地震小區劃項目報告》禮縣縣城及鹽官鎮地震小區劃項目報告.蘭州：甘肅省地震工程研究院，2013.所提供的禮縣縣城50年超越概率63%、10%、2%的人工計算場地基巖地震動時程分別作為輸入地震動（簡稱為工況一、工況二、工況三），其地震動加速度峰值分別為0.069? 97g、0.287 14g、0.518 72g。在保證計算精度的情況下，為方便求解，本次模擬選取0～20 s的數據，定義時間步長為0.02 s。數值模擬采用的地震波時程曲線見圖9。

3.3 結果分析

（1） 應力和變形分析

以工況三為例研究滑坡在地震作用下的應力和變形特征，將滑坡在天然自重作用下的變形歸零后，施加地震動加速度時程曲線。

地震作用后，滑坡體應力最大值和整個應力場（圖10）均發生了變化，滑坡體應力較震前有明顯的增加，并出現了應力集中帶。坡面上，由于底層巖體受到張拉力的作用，沿著坡面一定深度范圍內的應力值減小。潛在滑面上，由于主滑段上X向拉應力集中帶導致其上部剪應力有明顯的增加。

最大剪應變方面（圖11），在地震作用下，坡體產生了剪切變形。隨著地震動的作用，沿滑帶中上部出現剪應變集中帶，相對地震前應變值明顯增大，貫通率也顯著增大，地震初期（5 s、10 s）時伴隨著剪應變值在潛在滑帶和滑體內增大的同時，剪應變集中帶分布范圍沿著主滑段延伸，地震作用后期（15 s、20 s）剪應變值繼續在滑面和滑體內增大;主滑段最大剪應變集中帶繼續擴展，同時下部抗滑段上也出現了剪應變集中帶，最后基本貫通了整個滑面。

從位移分布（圖12～13）可以看出，滑坡的位移最大值為3.059 m，水平向位移分量最大值為2.726 m，出現在坡面中后部地形陡緩過渡處，即次級滑面和前緣潛在滑面后緣處。此外滑坡的后緣位移較前緣大，這與滑坡在變形破壞中后緣易出現拉裂縫變形現象相吻合。下覆基巖整體穩定，位移變化不大，上覆滑體位移較下伏基巖位移大一些，潛在滑面上主滑段和牽引段基本都出現了一定的位移，可見該滑坡會發生整體滑動。

（2） 穩定性評價

陳莊滑坡在地震動作用下沿現場勘探確定的推測主滑面和前緣潛在滑面（圖8）的安全系數時程曲線和Newmark永久位移時程曲線見圖14～16。

滑坡在工況一：50年超越概率63%，加速度峰值為0.069 97g地震動作用下（圖14），沿主滑面的安全系數最小值為1.139 984，發生在4.8 s時，永久位移為0;沿前緣潛在滑面的安全系數最小值為1.132 197，發生在16.8 s時，永久位移為0。可見滑坡在工況一下雖按最小安全系數判定處于基本穩定狀態，但沿主滑面和前緣潛在滑面都沒有滑動，是穩定的。

滑坡在工況二：50年超越概率10%，加速度峰值為0.287 14g的地震動荷載作用下（圖15），沿主滑面的安全系數最小值為1.028 769，發生在18.6 s時，永久位移為0;沿前緣潛在滑面的安全系數最小值為0.979 481，發生在7.4 s時，永久位移為0.000 91 m。可見該滑坡在工況二下沿主滑面按最小安全系數判定處于欠穩定狀態，沒發生位移;沿前緣潛在滑面按最小安全系數判定處于不穩定狀態，發生了很小的位移。

滑坡在工況三：50年超越概率2%，加速度峰值為0.518 72g的地震動荷載作用下（圖16），沿主滑面的安全系數最小值為0.692 95，發生在18.6 s時，永久位移為0.681 55 m;沿前緣潛在滑面的安全系數最小值為0.740 833，發生在7.4 s時，永久位移為1.189 13 m。可見滑坡在工況三下沿主滑面和前緣滑面都處于不穩定狀態，發生了滑動。

在地震荷載作用下，選用整個地震歷時中的最小安全系數評價滑坡的穩定性沒有考慮地震荷載作用的過程和邊坡破壞的積累性是不合理的［12］。本文采用文獻［12］提出的最小平均安全系數［式（1）］作為評價指標：

Fsmin=Fso-0.65（Fso-Fsmin）（1）

式中：Fso為邊坡在靜力作用下的安全系數;Fsmin為地震作用下安全系數隨震動過程而波動的最小值。

采用上述方法沿指定滑動面下所計算的最小平均安全系數列于表2，與上述Newmark永久位移結合最小安全系數判定的穩定性結論基本一致。同時隨著輸入地震動的增大，最小平均安全系數逐漸減小，滑坡穩定性降低，其中在50年超越概率2%的基巖地震動時程作用下陳莊滑坡處于不穩定狀態。

4 結論

本文以禮縣陳莊黃土－泥巖接觸面滑坡為例，在對其發育特征調查的基礎上，對其在地震作用下的動力響應進行了數值分析，從而對地震作用下黃土－泥巖接觸面滑坡的失穩機理和穩定性進行了動態分析研究。結果表明：

（1） 陳莊黃土－泥巖接觸面滑坡長約1 250 m，寬約1 300 m，高約330 m，厚約30 m，體積約4.875×107 m3，坡度為30°，屬巨型滑坡。該滑坡呈現出如下特征：簸箕形的平面形態，陡峭的后壁，兩側分布著雁形排列的羽狀張裂縫，前緣大量黃土堆積，坡內發育著眾多沖溝，落水洞及裂縫，黃土上覆于第三系泥巖的地層結構。這些特征為陳莊黃土－泥巖接觸面滑坡失穩機理的探討提供了研究基礎。

（2） 陳莊滑坡在地震誘發作用下失穩機理過程為：地震作用下，首先由于底層巖體的張拉作用，導致潛在滑面主滑段中上部剪應力較地震前有明顯的增加，出現了剪應變集中帶;隨著地震動的作用，剪應變集中帶從主滑段中上部向兩端延伸擴展，直至在潛在滑面中形成一個貫通的滑動區域導致滑坡失穩破壞，是受拉和受剪的復合破壞作用，即地震作用→底部巖體受張拉→潛在滑面主滑段上部剪應力增大→剪應變集中帶出現→剪應變集中帶延伸擴展貫通→滑坡失穩。

（3） 隨著輸入地震動的增大，滑坡穩定性降低：陳莊滑坡在50年超越概率63%的基巖地震動作用下沿主滑面處于穩定狀態，沿前緣潛在滑面處于基本穩定狀態;在50年超越概率10%的基巖地震動作用下沿主滑面處于基本穩定狀態，沿前緣潛在滑面處于欠穩定狀態;在50年超越概率2%的基巖地震動作用下沿主滑面和前緣潛在滑面都處于不穩定狀態。

參考文獻（References）

［1］文寶萍，王思敬，王恩志，等.黃土—紅層接觸面滑坡的變形特征［J］.地質學報，2005，79（1）：144.

WEN Baoping，WANG Sijing，WANG Enzhi，et al.Deformation characteristics of landslide at loess－red bed interface［J］.Acta Geologica Sinica，2005，79（1）：144.

［2］李媛，吳奇.孟家山黃土－紅層接觸面滑坡破壞機理研究［J］.水文地質工程地質，2001，28（1）：52－54.

LI Yuan，WU Qi.The study of deformation mechanism on loess－red mudstone sliding surface in Mengjiashan landslide［J］.Hydrogeology and Engineering Geology，2001，28（1）：52－54.

［3］石瑞紅.黃土－紅層接觸面滑坡形成演化機理［D］.蘭州：蘭州大學，2013.

SHI Ruihong.Formation and evolution mechanism of the loess－red beds interface landslide［D］.Lanzhou：Lanzhou University，2013.

［4］鄭程波.內蒙古自治區黃土泥巖接觸面滑移的數值模擬［D］.呼和浩特：內蒙古工業大學，2013.

ZHENG Chengbo.The numerical simulation of sliding about contact surface between loess and mudstone in Inner Mongolia Autonomous Region［D］.Hohhot：Inner Mongolia University of Tehchnology，2013.

［5］馬健.池溝滑坡成因機理與穩定性分析［D］.蘭州：蘭州大學，2014.

MA Jian.Analysis on genetic mechanism and stability of the Chigou landslide［D］.Lanzhou：Lanzhou University，2014.

［6］朱建東，鄢好，李紹紅，等.黃土－泥巖接觸面滑坡的兩種雨型模型試驗［J］.工程地質學報，2019，27（3）：623－631.

ZHU Jiandong，YAN Hao，LI Shaohong，et al.Laboratory model experiment of landslides along loess－mudstone interface induced by rainfall patterns［J］.Journal of Engineering Geology，2019，27（3）：623－631.

［7］陳永明，石玉成，徐暉平，等.1995年永登地震黃土震陷變形特征及其形成機理［J］.西北地震學報，2000，22（4）：465－470，475.

CHEN Yongming，SHI Yucheng，XU Huiping，et al.Deformation characteristic of seidmic subsidence and its formation mechanism of loess of the Yongdeng earthquake in 1995［J］.Northwestern Seismological Journal，2000，22（4）：465－470，475.

［8］石玉成，李蘭.黃土震陷變形特征的細觀分析［J］.巖石力學與工程學報，2003，22（增刊2）：2829－2833.

SHI Yucheng，LI Lan.Micro－analysis on seismic subsidence characteristics of loess［J］.Chinese Journal of Rock Mechanics and Engineering，2003，22（Suppl02）：2829－2833.

［9］王蘭民.黃土動力學［M］.北京：地震出版社，2003.

WANG Lanmin.Loess dynamics［M］.Beijing：Seismological Press，2003.

［10］吳錦奎.甘肅省禮縣地質災害發育特征及危險性評價［D］.蘭州：蘭州大學，2002.

WU Jinkui.Development characteristics and risk assessment of geological disasters in Lixian County，Gansu Province［D］.Lanzhou：Lanzhou University，2002.

［11］NEWMARK N M.Effects of earthquakes on dams and embankments［J］.Géotechnique，1965，15（2）：139－160.

［12］劉漢龍，費康，高玉峰.邊坡地震穩定性時程分析方法［J］.巖土力學，2003，24（4）：553－556，560.

LIU Hanlong，FEI Kang，GAO Yufeng.Time history analysis method of slope seismic stability［J］.Rock and Soil Mechanics，2003，24（4）：553－556，560.

（本文編輯：張向紅）

收稿日期：2020-10-09

基金項目：甘肅省地震局地震科技發展基金項目（2021M7，2016M04）

第一作者簡介：馬紫娟（1983- ），女，碩士，高級工程師，主要從事地震地質災害方面的研究。E－mail：65052042@qq.com。

通信作者：張有龍（1968-），男，碩士，高級工程師，主要從事地震地質和工程地震方面的研究和應用。E－mail：Ljxx1968＠163.com。