邊坡地震作用下穩定性研究評價

宋常旭

摘 要：地震對于邊坡作用是造成邊坡失穩的主要原因。文章對地震作用下邊坡穩定性分析的多種方法進行了分析總結，根據現有的理論方法，總體上分成了確定性方法和不確定方法兩個大類。其中確定性方法是當前主流的研究方法，并且從定性方法和定量方法兩種手段進行了相應的介紹、分析、評述，分析了各種方式的特點及適用性，從而提出了現階段研究上存在的不足，希望可以對未來研究發展的方向起到一定幫助。

關鍵詞：地震 邊坡穩定性 評價方法 確定性

中圖分類號：P31 文獻標識碼：A 文章編號：1674-098X（2016）05（a）-0020-04

地震是造成人員傷亡損失的嚴重的自然災害，地震作用下的邊坡穩定性更是巖土工程、地質工程、地震工程關注的重要問題。對于地震作用下邊坡穩定性的分析評價方法的有很多分類，比如：祁生文[1]將邊坡地震穩定性評價方法歸納為：擬靜力法、有限滑動位移法、Makdisi-Seed簡化分析法、剪切楔法、概率分析方法以及數值方法。劉紅帥等[2]從是否考慮邊坡巖土體參數的不確定性角度，將巖土邊坡地震穩定性評價方法分為確定性方法和概率分析方法兩大類。趙體[3]從評價原理及手段的角度出發全面考慮邊坡地震穩定性分析的評價方法，將邊坡地震穩定性分析方法大致分為統計分析、物理試驗和計算分析3類。以上分類都全面地涵蓋了國內外地震作用下邊坡穩定性的分析方法，但是對于多學科交叉的方法，并沒有太多涉及，因此，該文在已有的分析方法上，對地震作用下邊坡穩定性的分析方法中為確定性方法進行了分類、評價。

1 定性分析

定性分析是一種使用簡單方便，可以快速評價地震作用下邊坡穩定性的方法。它根據地震滑坡理論，對已發生地震的地區震后調查的大量結果分析，對震災區域潛在滑坡體做出快速、客觀、準確地判斷，其中，工程地質調查法和資料統計為兩種常用方法。

辛鴻博等[4]通過對我國地震觸發的285例典型邊坡崩滑實例進行分析研究，提出了以場地烈度、震級、震中距為標準的地震邊坡崩滑初判準則。丁彥慧等[5]通過對我國地震引起的194例地震崩滑實例的研究分析，建立了適用于我國西部地區的分級判別地震崩滑的預測方法。建立了地震崩滑最大震中距與地震震級間存在的對應關系，提出了一種判定邊坡穩定性綜合指標法，并且得出當邊坡的綜合指標值小于穩定邊坡臨界值時，邊坡穩定，否則邊坡失穩。王余慶等[6]通過大量的已有地震滑坡資料分析，在第三級預測法中進行改進，在綜合考慮各種條件的共同作用的基礎上，確定以巖性構造、坡角、坡高、地震、降雨為主要影響因子，提出了各因子分級賦值的原則和方案，建立了可對不同因子進行區分和量化的含5個因子的綜合指標法計算式，得出邊坡危害指數，并結合實際地震進行了對比驗證。黃潤秋等[7]通過“5.12”汶川地震災后地質災害的現場調查和遙感解譯手段獲得11 308處地質災害點，對地震地質災害的分布與距發震斷裂距離坡度、高程、巖性等因素的關系進行統計分析，指出了這些地震地質災害的發生存在明顯的上盤效應；距離發震斷裂越近地質災害發育密度越高；大多數的崩滑災害多發生在河谷岸坡度較陡，巖體卸荷強烈，地震響應最為突出的部位。王秀英等[8]通過收集整理汶川地震記錄的近3 000個崩滑點數據進行了地震誘發滑坡與地震動峰值加速度之間的對應關系研究，發現震動峰值加速度與地震誘發崩滑之間存在非常明顯的正相關性。不同地質區域對應斜坡臨界加速度有所不同，大多數在0.05～0.15 g之間變化，峰值加速度如果超過局部場地斜坡臨界加速度后，更加可能誘發滑坡的發生。在龍門山震區存在0.2 g的峰值加速度分界線，高于臨界值地震滑坡造成的地質災害嚴重。

定性分析法是地震邊坡穩定性判斷分析最直接的一種方法，能準確、快速判斷地震作用下邊坡的總體規律，對地震造成的邊坡失穩，以及工程抗震提供了有力依據。但是資料統計需要大量的數據以及資料，對于時間比較遠的地震資料，往往記錄的較少，這樣對統計結果造成一定影響。工程地質調查法在復雜的邊坡工程中有著明顯的優勢，但缺點也是顯而易見的，如需要更多的人員的現場調查，或者更多采集數據的手段做支撐，耗費大量人力物力，做一項完整詳細的調查，需要很長的時間。現階段，技術手段相比以往已經很先進，但還不能足夠快速地進行評定。

2 定量分析

定量分析的方法在地震作用下的邊坡穩定性評價中主要是通過一定的數據收集，進行一系列復雜分析計算，進而得到一些盡可能精確的數據來對邊坡的穩定性進行評價，給出一個較為準確的結論，對工程和實際研究有一定的指導作用。其中對于邊坡地震作用下穩定性計算中，國內外常用的方法為計算邊坡的地震安全系數或是最大滑動位移。

2.1 安全系數法

考慮地震安全系數判斷或評價邊坡穩定性有多種計算方式，其中擬靜力法是最常用的一種方式。

1950年Terzaghi首次將擬靜力法用于計算地震作用下邊坡穩定性問題中，此方法是基于靜力計算，將地震作用力簡化為水平和豎直兩個方向的不變的加速度作用，作用于滑坡體的重心，根據極限平衡理論建立靜力平衡方程和力矩平衡方程，求解邊坡的穩定性安全系數。擬靜力法概念清晰，計算工作量小，因此得到了廣泛應用，并且被編入到相關規范中。對于擬靜力法，關鍵的問題是地震系數的選取。Terzaghi[9]曾經指出，一般性地震情況，水平向系數取0.1；破壞性地震情況，系數取0.2；災難性地震情況，系數取0.5。國內方面，呂擎峰等[10]認為傳統的擬靜力法采用豎向條分法進行受力分析不能準確合理地計算地震慣性力矩，會產生誤差。他們在通過計算豎向力時采用豎向條分，計算水平地震力對滑動力矩和抗滑力矩的作用時采用水平條分的方法計算安全系數，但是這樣得到的安全系數小于傳統擬靜力法計算得到的安全系數。張建海等[11]提出了采用剛體彈簧元計算邊坡等在地震波作用下的動力安全系數，此方法得出的安全系數是隨地震波的作用而發生變化，從而更深刻地反映了動力現象本質。

擬靜力法建立在假定理想的剛體基礎上，實際的邊坡并非剛性體，同時也很難考慮地震的特性以及邊坡的材料的動力特性，因而，擬靜力方法不能特別準確地評價邊坡地震作用下的穩定性。

通過數值模擬計算安全系數的方法已經被國內外廣泛應用。常用的方法的有有限元法、有限差分法、離散元法。薄景山等[12]使用有限元分析得到邊坡動力作用下的應力場，再分別對每一時刻下確定的折線滑裂面和圓弧滑裂面上的抗滑力和下滑力進行積分，得到安全系數時程曲線，來確定為最小安全系數。但是用最小安全系數的方法過于保守，當小于最小安全系數時，邊坡不一定徹底破壞。因此，劉漢龍等[13]認為用整個地震歷時中的最小安全系數評價邊坡的地震穩定性不太合適，提出了用最小平均安全系數作為評價指標，即利用安全系數最大振幅的0.65倍作為平均振幅來反映安全系數隨地震波動變化的過程。吳兆營等[14]采用有限元計算某巖體邊坡的動力安全系數時程，并給出時間加權平均安全系數。此外，馬芳芳[15]提出以時間長度上安全系數的平均值作為評價指標，在可能滑裂區內搜索最小安全系數均值及對應的最危險滑裂面，考慮了地震荷載對邊坡失穩的時間作用。鄧東平等[16]基于滑動面搜索新方法對地震作用下邊坡穩定性擬靜力分析，通過算例對比分析，表明此方法計算得的安全系數略小于簡化畢肖普、楊布圓弧法和極限分析上限法，得到的臨界滑動面與臨界圓弧滑動面也頗為相似，而在非均質土坡中表現為非圓弧型。

數值模擬的方法相比擬靜力方法的優勢顯而易見，它能解決復雜條件下邊坡動力問題。然而，也對計算機性能和軟件硬件條件提出了更高的要求。此外，對于數值軟件分析是所用到的模型、參數、邊界條件、本構關系、加載方式等合理地選取，如何更加準確合理地使用，在未來的研究中還需要更加深入地分析研究。

目前，安全系數法是國內外使用最多的用來評價邊坡地震作用下穩定性的方法，但是仍然存在如何確定最小安全系數這個難題。當達到最小安全系數時，邊坡仍然有可能并未完全滑落，這對實際工程和災后的判斷產生了重要影響。因此，如何綜合考慮邊坡地震作用下穩定性的各方面因素，提出一種合理的最小安全系數，并且是研究與工程實際聯系起來，更加合理地應用到實際的工程中，是未來研究中的核心。

2.2 永久位移法

邊坡地震作用下穩定性分析另外一種重要的定量分析是根據永久位移分析邊坡的穩定性。這是一種有別于安全系數法判別穩定性的方法，主要有Newmark滑塊位移法和能量分析法兩大類。

Newmark滑塊分析法是邊坡動力穩定性分析的常用方法之一。它根據地震作用下永久位移的標準來評價邊坡的地震穩定性。Newmark法最關鍵的是屈服地震系數。Bray等[17]根據Newmark滑塊機理，考慮了初始屈服地震系數（Ky）、初始基本周期（Ts）和當退化周期為1.5Ts時地面運動的譜加速度的影響，提出了一種簡單的預測邊坡地震作用下永久位移的模型。祁生文[18]提出了一種考慮了結構面有起伏的條件下，結構面退化的巖質邊坡的地震永久位移算法。

Newmark滑塊位移法只考慮了地震作用過程中所產生的永久位移，而忽略了計算地震作用后的永久位移。董建華等[19]提出了永久位移由地震作用過程中的位移和地震發生后位移兩部分組成，并且建立地震作用后永久位移計算模型并求解，計算結果表明兩個位移量的比值很大，因此震后位移不能忽略。通過試驗模型研究位移變形是一種最為直觀的方式。王思敬[20]將有限滑動位移法思路引入到動力條件下邊坡穩定性分析中，采用大型振動臺試驗驗證了基巖振動加速度可以超過極限加速度值而不一定導致動力失穩，但可能產生一定的最終位移。

國內外許多學者根據滑塊位移法思想，通過數值模擬的方式研究邊坡穩定性。陶連金等[21]應用動力離散元計算了某節理巖體邊坡的永久位移，并對其進行了穩定性評價。劉春玲等[22]利用FLAC3D對某一實際邊坡進行了動力分析，討論了利用FLAC3D進行邊坡動力分析時如何對參數進行選擇，結果表明，在地震作用下，邊坡會發生一定的永久位移。

以往的Newmark滑塊位移法計算條件簡單，并沒有考慮豎向地震動和孔隙水壓力條件。劉忠玉等[23]選用A型曲線對地震荷載作用下飽和黃土的孔壓變化規律進行了分析發現：孔壓累積會降低屈服地震系數，且與常規的Newmark法分析的進行了比較，結果顯示考慮孔壓累積時的永久位移計算值要大于常規的Newmark法分析的結果。黃建梁等[24]提出并建立了根據水平和豎直地震加速度時程估計坡體穩定的加速度、速度和位移時程的方法；討論了地震動加速度時程的確定問題，評價了地震過程中坡體抗滑強度的衰減問題和孔隙的動態響應問題以及坡體地震穩定性的評定問題，根據實例得出了豎直地震動對地震位移計算的影響很小，在位移計算中可忽略豎直地震動的影響的結論。此外，李紅軍[25]等利用改進的Newmark滑塊位移法，采用時程豎向加速度計算屈服角加速度，運用擬靜力極限平衡分析和地震動力反應分析計算土工結構地震永久變形，根據算例表明：算屈服角加速度計算采用時程豎向加速度得到的滑動位移位于豎向加速度假定恒定向上或向下之間。

永久位移的方法很好的反應了地震作用下邊坡累計效應，用宏觀的位移量來表示邊坡變形的情況，也是不同于安全系數法這樣的理論計算方法，此外，該方法在對地震動、邊坡動力特性有著準確認識和對地震反應有著合理模擬。因此，近些年受到了國內外的專家學者的廣泛關注和深入研究。但是，永久位移的方法還存在很多的問題。在Newmark方法中，地震運動是三維的，但在進行邊坡地震永久位移計算時常采用一維水平地震動輸入；在計算永久位移時，對巖土體材料、特性以及參數的變化的研究還不夠明確；Newmark法計算出的為一個平均位移，這樣造成了無法確定邊坡潛在滑移面各點的分布規律。近些年，雖有一些專家學者進行過一些大型振動臺或離心機試驗，但用于邊坡地震穩定評價還不夠明確。基于能量分析法判斷永久位移，雖有相關研究，但數量還是很少，很難應用到實際工程或是穩定性評價。

3 結語與展望

該文就邊坡地震作用下穩定性分析問題從過往的研究、發展和分析方法方面進行了總結，并進行了分類。雖然國內外的專家學者在邊坡穩定性分析評價方面已經取得了大量的富有成效的成果以及結論，但是仍在存在一定的不足以及未來發展的需要關注的方面。

（1）定性分析方法依靠工程地質調查和資料統計手段來進行分析評價，但對資料與數據的數量上來說，還是相對的比較缺乏，但是今后如何建立數據庫，以及將已有數據和資料與更多新的方法理論聯系起來，用今后的工程實例來驗證，可能未來發展的新方向。

（2）根據安全系數法確定穩定性方法是國內外最為常用的方法，在對于地震系數的研究上已經取得了一定的效果，但在地震動峰值加速度的放大系數的研究還不夠深入，還不能準確的給出快速準確的給出地震系數。

（3）計算機軟件的發展，依據相關軟件分析提出的方法也是屢見不鮮，但如何將在地震作用下產生邊坡失穩的各種因素，比如地震動參數，巖土體等各種參數都考慮全面，綜合起來，得出更為準確的安全系數，并且用已發生的實例進行詳細分析對比，從而應用到未來的地質災害中或是工程案例中。

（4）豎向地震動已經被認為對邊坡穩定性等影響起到了一定的作用，不少專家學者已經有相關研究，但是對于豎向地震動在地震作用下如何產生作用，以及對邊坡穩定性影響所占的比重是多少的相關研究卻很少。

參考文獻

[1]祁生文.邊坡動力響應分析及應用研究[D].北京：中國科學院地質與地球物理研究所，2002.

[2]趙春明.巖土邊坡地震穩定性評價方法研究進展[J].防災科技學院學報，2007（3）：20-27.

[3]趙體，李麒，白田子.邊坡地震穩定性評價方法研究進展[J].華北科技學院學報，2014（10）：49-54.

[4]辛鴻博，王余慶.巖土邊坡地震崩滑及其初判準則[J].巖土工程學報，1999（5）：591-594.

[5]丁彥慧，王余慶，孫進忠，等.地震崩滑預測方法及其工程應用研究[J].工程地質學報，2000（4）：475-480.

[6]王余慶，辛鴻博，高艷平，等.預測巖土邊坡地震崩滑的綜合指標法研究[J].巖土工程學報，2001（3）：311-314.

[7]黃潤秋，李為樂.“5.12”汶川大地震觸發地質災害的發育分布規律研究[J].巖石力學與工程學報，2008（12）：2585-2592.

[8]王秀英，聶高眾，王登偉.汶川地震誘發滑坡與地震動峰值加速度對應關系研究[J].巖石力學與工程學報，2010（1）：82-89.

[9]Cogen william M. Mechanics of Landslides[D].Cailfornia Institite of Technology，2015.

[10]呂擎峰，殷宗澤，王叔華，等.擬靜力法邊坡穩定分析的改進[J].巖土力學，2005（S1）：35-38.

[11]張建海，范景偉，何江達.用剛體彈簧元法求解邊坡、壩基動力安全系數[J].巖石力學與工程學報，1999（4）：19-23.

[12]薄景山，徐國棟，景立平.土邊坡地震反應及其動力穩定性分析[J].地震工程與工程振動，2001（2）：116-120.

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

[14]吳兆營，薄景山，劉紅帥，等.巖體邊坡地震穩定性動安全系數分析方法[J].防災減災工程學報，2004（3）：237-241.

[15]馬芳芳.基于地震動力時程反應的有限元邊坡穩定性分析[D].大連理工大學，2005.

[16]鄧東平，李亮.基于滑動面搜索新方法對地震作用下邊坡穩定性擬靜力分析[J].巖石力學與工程學報，2012（1）：86-98

[17]Jonathan D. Bray， Thaleia Travasarou. Pseudostatic Coefficient for Use in Simplified Seismic Slope Stability Evaluation [J]. Journal of Geotechnical and Geoenvironmental Engineering，2009，135（9）：1336-1340.

[18]祁生文.考慮結構面退化的巖質邊坡地震永久位移研究[J].巖土工程學報，2007（3）：452-457.

[19]董建華，朱彥鵬.地震作用下土釘支護邊坡永久位移計算方法研究[J].工程力學，2011（10）：101-110.

[20]王思敬，張菊明.邊坡巖體滑動穩定的動力學分析[J].地 質科學，1982（2）：162-170.

[21]陶連金，蘇生瑞，張倬元.節理巖體邊坡的動力穩定性分析[J].工程地質學報，2001（1）：32-38.

[22]劉春玲，祁生文，童立強，等.利用FLAC～（3D）分析某邊坡地震穩定性[J].巖石力學與工程學報，2004，23（16）：2730–2733.

[23]劉忠玉，魏建東.飽和黃土邊坡的動力穩定性分析[J].巖土力學，2005（2）：198-202.

[24]黃建梁，王威中，薛宏交.坡體地震穩定性的動態分析[J].地震工程與工程振動，1997（4）：114-123.

[25]李紅軍，遲世春，林皋.平均屈服加速度的Newmark滑塊位移法[J].哈爾濱工業大學學報，2009（10）：100-104.