王 可，葉曉冬，潘志成，閆茂林，溫一波

（1.云南大學(xué)建筑與規(guī)劃學(xué)院，云南昆明650500；2.云南大學(xué)資源環(huán)境與地球科學(xué)學(xué)院，云南昆明650500）

小灣滑坡地震穩(wěn)定性研究

王 可1，葉曉冬1，潘志成1，閆茂林1，溫一波2

（1.云南大學(xué)建筑與規(guī)劃學(xué)院，云南昆明650500；2.云南大學(xué)資源環(huán)境與地球科學(xué)學(xué)院，云南昆明650500）

滑坡災(zāi)害破壞性大、治理難度高，特別是地震作用下的穩(wěn)定性問題十分復(fù)雜，是土木工程界的一大難題。以小灣滑坡為研究對象，運用有限元法計算分析小灣滑坡在小震、中震、大震作用下的穩(wěn)定性，得到小灣滑坡在地震作用下的穩(wěn)定系數(shù)時程曲線，并以此為依據(jù)對小灣滑坡的地震穩(wěn)定性進(jìn)行評價，研究了小灣滑坡穩(wěn)定系數(shù)與地震強(qiáng)度的關(guān)系。

小灣滑坡；地震穩(wěn)定性；有限元法

我國地震災(zāi)害頻繁，并且具有震源淺、強(qiáng)度高、破壞性大的特點，而滑坡災(zāi)害也是常見的地震次生災(zāi)害之一［1］。王秀英等［2］運用隸屬函數(shù)法對模型累積位移隸屬度進(jìn)行計算，并以此判定滑坡的穩(wěn)定程度。許沖等［3］將四川省蘆山縣太平鎮(zhèn)東北部一個矩形區(qū)作為研究區(qū)，統(tǒng)計分析了地震滑坡的密度與地層巖性、地形、地震強(qiáng)度的關(guān)系。方薇等［4］在不平衡推力法中引入比例因子以確定條間力傾角，建立力的平衡方程與力矩平衡方程，計算滑坡地震穩(wěn)定系數(shù)并驗證了其合理性。顧成壯等［5］研究了滑坡的啟動機(jī)制，從而推導(dǎo)出滑坡啟動速度的公式。毛彥龍等［6］提出累進(jìn)破壞效應(yīng)、啟程加速效應(yīng)和啟動效應(yīng)導(dǎo)致滑坡發(fā)生的觀點。陳曉利等［7］提出以研究區(qū)域地形、地質(zhì)、地貌等為影響因子，綜合評價滑坡災(zāi)害危險等級的方法，并以此對汶川地震進(jìn)行了地震危險等級劃分并證明其合理性。本文在前人研究的基礎(chǔ)上，增添新資料，以小灣滑坡為實例，運用有限元數(shù)值模擬法分析小灣滑坡在小震、中震、大震作用下的穩(wěn)定性。根據(jù)計算結(jié)果，對小灣滑坡地震穩(wěn)定性進(jìn)行了研究分析。

1 資料與方法

1.1 小灣滑坡

小灣滑坡地處云南省臨滄市鳳慶縣，根據(jù)《建筑抗震設(shè)計規(guī)范》［8］（GB50011-2010）規(guī)定：鳳慶縣屬抗震設(shè)防烈度Ⅷ度區(qū)第二組，設(shè)計基本地震加速度值為0.20 g，并規(guī)定抗震設(shè)防烈度為Ⅷ度的地區(qū)，多遇地震的地震加速度時程最大值為0.07 g，罕遇地震的地震加速度時程最大值為0.4 g。小灣滑坡區(qū)域范圍圖見圖1，圖2為小灣滑坡的剖面示意簡圖。

圖1 小灣滑坡區(qū)域范圍圖

圖2 小灣滑坡剖面示意簡圖

《云南省鳳慶縣小灣鎮(zhèn)政府駐地滑坡勘察及防治工程可行性研究報告》中，小灣滑坡長約140 m，平均寬200 m，平面形態(tài)近半圓型，主滑方向78°，滑塊平均厚約6 m，體積約19.2×104m3，屬于淺層中型滑坡。土體類型主要是粉質(zhì)黏土和弱風(fēng)化花崗片麻巖。

小灣滑坡動力分析土體參數(shù)取值如表1所示。

表1 小灣滑坡動力分析參數(shù)取值表

粉質(zhì)黏土中，土體應(yīng)力狀態(tài)是最大動剪切模量的重要影響因素，不能取某一常數(shù)值［10-11］，可以根據(jù)土體的最大埋深、超固結(jié)比、孔隙比等估計最大動剪切模量與Y向有效應(yīng)力的關(guān)系［12］，估算得到的關(guān)系曲線如圖3所示。

圖3 粉質(zhì)黏土最大動剪切模量-Y向有效應(yīng)力估算曲線

1.2 相關(guān)方法

1.2.1 分析流程

分析小灣滑坡地震穩(wěn)定性的流程如圖4所示。

圖4 地震穩(wěn)定性分析流程圖

本次研究涉及到GeoStudio軟件中兩大板塊——QUAKE/W板塊和SLOPE/W板塊，在研究中應(yīng)注意邊界條件的設(shè)置，在初始應(yīng)力狀態(tài)分析時，將滑坡模型兩側(cè)設(shè)置為X向約束，下端設(shè)置為X/Y向約束；而在進(jìn)行動力分析時，兩側(cè)約束為Y方向約束。因為滑坡在地震作用下的破壞表現(xiàn)為滑坡體的水平移動，而水平地震作用對滑坡有著較大的破壞［14］。

本次滑坡地震穩(wěn)定性研究的方法是基于有限元的傳遞系數(shù)法，GeoStudio提供了多種網(wǎng)格劃分的方法，本次研究網(wǎng)格劃分統(tǒng)一采用Quads&Triangles方式進(jìn)行劃分。單元的尺寸根據(jù)滑動面的大小獨立設(shè)置，不采用統(tǒng)一的單元尺寸［15］。

1.2.2 地震加速峰值的選取

在對小灣滑坡進(jìn)行動力分析時，輸入的峰值加速度分別按大震、中震、小震的加速度時程最大值取值。峰值加速度取值方法見表2、表3。

表2 抗震設(shè)防烈度與設(shè)計基本加速度值

表3 所用地震加速度時程的最大值

1.2.3 地震波的選取

本次對小灣滑坡的地震穩(wěn)定性的研究選用唐山波進(jìn)行分析。

對9例蘇醒后顱腦損傷昏迷病人進(jìn)行訪談，發(fā)現(xiàn)病人心理體驗主要表現(xiàn)在病人對生命價值認(rèn)同、對健康狀況存在期望、擔(dān)憂工作生活改變、希望建立社會支持系統(tǒng)等方面，情緒體驗可進(jìn)一步歸納為2種類型：積極情緒體驗、消極情緒體驗。

唐山波持續(xù)時間40 s，時間間隔為0.02 s，峰值加速度為65.94 cm/s2，出現(xiàn)在第19.36 s，其加速度時程曲線如圖5所示。

圖5 唐山波加速度時程曲線

2 滑坡數(shù)值模擬結(jié)果與分析

2.1 唐山波動力分析

唐山波為入射地震波，等比例改變其加速度時程曲線中加速度的大小，使其加速度峰值為0.07 g、0.2 g、0.4 g，模擬分析小灣滑坡在小震、中震、大震環(huán)境下的穩(wěn)定性問題，計算得到小灣滑坡穩(wěn)定系數(shù)時程曲線見圖6～圖8。

圖6 小震穩(wěn)定系數(shù)時程曲線（0.07 g）

圖7 中震穩(wěn)定系數(shù)時程曲線（0.2 g）

圖8 大震穩(wěn)定系數(shù)時程曲線（0.4 g）

以唐山波為入射波時，小灣滑坡穩(wěn)定系數(shù)時程曲線關(guān)鍵數(shù)據(jù)匯總表見表4。

表4 小灣滑坡穩(wěn)定系數(shù)曲線關(guān)鍵數(shù)據(jù)匯總表

以唐山波為入射波，分析小灣滑坡在地震作用下的穩(wěn)定性時：

（1）小震影響下，小灣滑坡在第19.38 s出現(xiàn)首個小于1的穩(wěn)定系數(shù)0.9956，說明滑坡在第19.38 s發(fā)生滑動。最小穩(wěn)定系數(shù)0.9556表示滑坡在小震作用下的穩(wěn)定性較為合理［16］。該值比自然穩(wěn)定系數(shù)降低了22%。

（2）中震影響下，小灣滑坡在第0.22 s出現(xiàn)首個小于1的穩(wěn)定系數(shù)0.8851，說明滑坡在第0.22 s已經(jīng)開始滑動；出現(xiàn)在第19.38 s的最小穩(wěn)定系數(shù)0.6998僅能從理論上表示滑坡的穩(wěn)定程度。該值比自然穩(wěn)定系數(shù)降低了45%。

（3）大震影響下，小灣滑坡在第0.22 s出現(xiàn)首個小于1的穩(wěn)定系數(shù)0.6711，說明滑坡在第0.22 s已經(jīng)開始滑動；出現(xiàn)在第19.38 s的最小穩(wěn)定系數(shù)0.4707僅能從理論上表示滑坡的穩(wěn)定程度。該值比自然穩(wěn)定系數(shù)降低了63%。

小震、中震、大震時滑坡穩(wěn)定系數(shù)最小的時刻均為第19.38 s，比地震加速度峰值出現(xiàn)的時刻晚0.02 s。

2.2 數(shù)值模擬合理性評價

為了驗證數(shù)值模擬結(jié)果的合理性，現(xiàn)運用Janbu法計算滑坡在小震、中震、大震作用下的穩(wěn)定系數(shù)，與數(shù)值模擬的結(jié)果形成對比，如表5所示。

表5 小灣滑坡數(shù)值分析精度評價表

表5顯示，Janbu法的計算結(jié)果和數(shù)值模擬的結(jié)果較為接近?？烧J(rèn)為本次數(shù)值模擬結(jié)果基本合理［16］。

3 結(jié) 論

本次研究以唐山波作為入射波，運用有限元法對小灣滑坡在小震、中震、大震影響下的穩(wěn)定性進(jìn)行了分析，得到了一些初步的結(jié)論：

（1）穩(wěn)定系數(shù)時程曲線最低點所對應(yīng)的時間點，并不會因緊隨地震強(qiáng)度的改變而改變，與地震波加速度峰值對應(yīng)的時刻并不重合，可能會提前或者推遲。

（2）小灣滑坡在小震影響下可以保持整體穩(wěn)定性，在中震和大震作用下則可能發(fā)生劇烈的滑動。

［1］ 黃潤秋.20世紀(jì)以來中國的大型滑坡及其發(fā)生機(jī)制［J］.巖石力學(xué)與工程學(xué)報，2007，26（3）：434－454.

［2］ 王秀英，聶高眾，王 松.利用模糊數(shù)學(xué)方法建立汶川地震滑坡災(zāi)害評判標(biāo)準(zhǔn)［J］.巖土力學(xué)，2011，32（2）：403－410.

［3］ 許 沖，肖建章.蘆山地震滑坡空間分布分析［J］.地震地質(zhì)，2013，35（2）：436－451.

［4］ 方 薇，孫正兵，何曉霞，等.考慮整體力矩平衡的改進(jìn)不平衡推力法［J］.四川大學(xué)學(xué)報（自然科學(xué)版），2012，44（3）：28－32.

［5］ 顧成壯，胡卸文，羅 剛，等.地震滑坡啟程動力學(xué)機(jī)理［J］.西南交通大學(xué)報，2012，47（4）：567－572.

［6］ 毛彥龍，胡廣韜，趙法鎖，等.地震動觸發(fā)滑坡體滑動的機(jī)理［J］.地球科學(xué)與環(huán)境學(xué)報，1998，20（4）：45－48.

［7］ 陳曉利，冉洪流，王明明.潛在地震滑坡危險區(qū)區(qū)劃方法［J］.地球物理學(xué)報，2012，55（4）：1267－1277.

［8］ 中華人民共和國住房和城鄉(xiāng)建設(shè)部.建筑抗震設(shè)計規(guī)范：GB50011-2010［S］.北京：中國建筑工業(yè)出版社，2010.

［9］ 《工程地質(zhì)手冊》編寫委員會.工程地質(zhì)手冊［M］.3版.北京：中國建筑工業(yè)出版社，1992.

［10］ 齊文浩，薄景山.土層地震反應(yīng)等效線性化方法綜述［J］.世界地震工程，2007，23（4）：221－226.

［11］ 陳 堯.二種土層地震反應(yīng)模擬方法的比較研究和昆明市地鐵三號線場址土層地震反應(yīng)分析［D］.昆明：云南大學(xué)，2012.

［12］ 溫一波，趙洞明，陳 堯.呈貢新區(qū)剪切波特征與分布［J］.云南大學(xué)學(xué)報（自然科學(xué)版），2012，34（S2）：221－226.

［13］ 谷德振著.巖體工程地質(zhì)力學(xué)基礎(chǔ)［M］.北京：科學(xué)出版社，1979.

［14］ 湯 卓，謝建斌，卞榮森，等.材料非飽和對尾礦庫三維滲流影響分析［J］.水利與建筑工程學(xué)報，2016，14（1）：15－20.

［15］ 錢玉林，劉長林，呂 媛.基于三維數(shù)值分析的土釘支護(hù)結(jié)構(gòu)研究［J］.水利與建筑工程學(xué)報，2016，14（2）：55－60.

［16］ 姚晨輝，魏江波，任光明，等.考慮巖土體參數(shù)變異性的滑坡破壞概率分析［J］.水利與建筑工程學(xué)報，2016，14（2）：100－103.

Analysis of Xiaowan Landslide Stability in Earthquake

WANG Ke1，YE Xiaodong1，PAN Zhicheng1，YAN Maolin1，WEN Yibo2
（1.College of Architecture and Planning，Yunnan University，Kunming，Yunnan 650500，China；2.College of Resources Environment and Earth Sciences，Yunnan University，Kunming，Yunnan 650500，China）

Landslides as a kind of geological disaster have great destructive power and is very difficult to control.The stability of the landslide under the seismic action is very complicated.It is a multi－discipline research issue and still a big problem for civil engineering.Taking Xiaowan landslide as research objects this paper analyzed its stability under the action of small，medium and strong earthquakes by using finite element method，then the time－h(huán)istories curves of landslide stability coefficient under the influence of the earthquake was acquired.Finally，the stability of Xiaowan landslide in earthquake was analyzed combining with the data，the relationship between stability coefficient and seismic intensity was also analyzed.

Xiaowan landslides；tearthquake stability；finite element method

P315.9

A

1672—1144（2016）05—0177—04

10.3969/j.issn.1672-1144.2016.05.034

2016－06－17

2016－07－13

云南省自然科學(xué)基金資助項目（2009CD001）

王 可（1993—），女，山東菏澤人，碩士研究生，研究方向為地震引發(fā)滑坡。E－mail：610762362@qq.com

溫一波（1962—），男，四川成都人，教授，主要從事地震衰減特性研究，震波分析，地震勘測等方面的工作。E－mail：1213201943@qq.com