范　磊，劉　聰

(1.長安大學(xué) 環(huán)境科學(xué)與工程學(xué)院，陜西 西安 710054；2.西安市城市規(guī)劃設(shè)計(jì)研究院，陜西 西安 710082)

2008年汶川8級(jí)地震的重災(zāi)區(qū)14個(gè)縣(市)主要公路全部癱瘓，主要鐵路、通信、水庫和城鎮(zhèn)遭受嚴(yán)重破壞，特別是北川縣城和汶川的映秀鎮(zhèn)成為一片廢墟，是繼1976年唐山7.8級(jí)地震之后我國又一個(gè)死亡和損失巨大的毀滅性地震[1-2]。以往國內(nèi)外的震害研究表明：破壞性特大地震，在距震中300 km內(nèi)和地震烈度VII度區(qū)內(nèi)均能造成大量崩塌滑坡[3]；汶川地震對(duì)地質(zhì)環(huán)境造成的影響前所未有，也給中國地學(xué)界提出了一系列重大問題[4]。專家們認(rèn)為,地震不僅激活了部分已存在的崩塌滑坡體，還產(chǎn)生了大量新的地質(zhì)災(zāi)害，對(duì)地質(zhì)環(huán)境、生態(tài)環(huán)境和地下水環(huán)境等都帶來了一定的影響[1-4]。地震影響區(qū)的地質(zhì)環(huán)境將隨當(dāng)?shù)亟?jīng)受的地震動(dòng)烈度、地質(zhì)環(huán)境條件組合、地表位錯(cuò)量和外在因素干擾等而發(fā)生一定的變化，這種影響是十分罕見的。同時(shí)，地震對(duì)地質(zhì)環(huán)境的影響時(shí)間滯后持續(xù)效應(yīng)十分明顯[5]。

西安市位于關(guān)中盆地，南靠秦嶺，以秦嶺山脈主脊為界，東有黃土堆積臺(tái)塬地貌，如白鹿塬等，西以太白山及青化臺(tái)塬為界，南北最大縱距約100 km，東西最大橫距約204 km，總面積9 983 km2。距離汶川地震中央斷裂帶北段約350 km，屬于汶川地震波及區(qū)，受汶川地震影響，西安全市發(fā)生多起地質(zhì)災(zāi)害，對(duì)交通、旅游帶來嚴(yán)重影響[6]。西安市地質(zhì)災(zāi)害主要為受斷層、地下水開采等因素誘發(fā)的地裂縫與地面沉降問題和受河流切割、地貌抬升影響產(chǎn)生的崩塌、滑坡和泥石流問題。特殊的區(qū)域特征，使得西安地區(qū)地質(zhì)災(zāi)害類型多、分布廣，西安地區(qū)屬于汶川地震Ⅴ～Ⅵ度烈度區(qū)，這些地質(zhì)災(zāi)害在汶川地震中均受到一定程度的影響。因此，客觀評(píng)價(jià)汶川地震對(duì)西安地質(zhì)環(huán)境影響的程度，有效減輕地質(zhì)災(zāi)害對(duì)社會(huì)的不利影響，對(duì)西安社會(huì)經(jīng)濟(jì)可持續(xù)發(fā)展具有十分重要的現(xiàn)實(shí)意義。

地震誘發(fā)的崩塌滑坡研究在地震次生山地災(zāi)害最為廣泛和深入，已成為滑坡研究的一個(gè)重要領(lǐng)域。Keffer對(duì)1989年加利福尼亞6.9級(jí)地震誘發(fā)的滑坡研究發(fā)現(xiàn)，其發(fā)生與距斷裂帶距離呈很強(qiáng)的反相關(guān)，而與巖石類型和軟硬度的相關(guān)性較弱，這表明滑坡與巖石特性的關(guān)系較為復(fù)雜[7]。Chigira 和Yagi于2006年對(duì)日本2004年發(fā)生的Mid Niigta縣6.8級(jí)地震引發(fā)的滑坡地貌和地質(zhì)特征進(jìn)行了分析，發(fā)現(xiàn)深層滑坡中老滑坡和凹形坡受地震影響較大，此外平面滑坡面對(duì)形成災(zāi)難性的滑坡敏感性較強(qiáng)[8]。針對(duì)滑坡區(qū)域性評(píng)價(jià)方面，目前較為流行的是利用NEWMARK模型，對(duì)區(qū)域性永久位移進(jìn)行評(píng)價(jià)，獲取區(qū)域的永久位移，評(píng)價(jià)區(qū)域邊坡的穩(wěn)定性及地震對(duì)區(qū)域長期的影響作用[9-11]。對(duì)于單體滑坡而言，利用數(shù)值模擬評(píng)價(jià)其對(duì)不同烈度下的穩(wěn)定性是當(dāng)前較為常用的方法[12-13]。

本文根據(jù)西安地區(qū)地質(zhì)環(huán)境條件和汶川地震斷裂破壞情況，選取西安地區(qū)典型滑坡體耶柿村滑坡對(duì)其在地震作用下的穩(wěn)定性進(jìn)行評(píng)價(jià)，研究地震對(duì)西安地區(qū)邊坡的長期穩(wěn)定性的影響。

表1　耶柿村滑坡土層物理力學(xué)性質(zhì)

1　滑坡體概況及研究方法

1.1　耶柿村滑坡體概況

耶柿村滑坡位于西安市長安區(qū)魏寨鄉(xiāng)耶柿村，109°10′53″E，34°6′10″N，距西安市中心約30 km。該滑坡東西長1.0 km左右，南北寬約0.6 km，體積約900萬m3，屬大型黃土滑坡(圖1)。該滑坡總體呈上小下大階梯狀，滑坡后壁呈圈椅狀，坡度約46°。該滑坡由老滑坡和滑體上的次級(jí)復(fù)活體組成。滑坡體按其性質(zhì)和特征差異從上至下可分為2層：第一層為黃土狀土(粉質(zhì)粘土)，組成滑體土和滑帶土；第二層為第三系粘土巖，組成滑床。土層的物理力學(xué)特性見表1。

根據(jù)野外調(diào)查結(jié)果，耶柿村滑坡是一個(gè)發(fā)生在黃土塬邊的老滑坡，由于滑體地形較為平緩，有利于降水的入滲，且滑體的前部有耶柿村分布，農(nóng)田灌溉，滑體前部削坡建房和修路等人類工程活動(dòng)，使得老滑坡在強(qiáng)降雨時(shí)發(fā)生復(fù)活。根據(jù)記載和調(diào)查，該滑坡體在1976年和2003年暴雨后均有大量裂縫產(chǎn)生，并發(fā)現(xiàn)前緣緩慢前移現(xiàn)象。該滑坡體處于汶川地震Ⅴ度烈度區(qū)，汶川地震對(duì)該滑坡的穩(wěn)定性造城一定影響，地震后在滑坡體后緣出現(xiàn)裂縫。

圖1　耶柿村滑坡區(qū)域位置與遙感影像圖

1.2　GEOSTUDIO數(shù)值分析

本文采用GEOSTUDIO中的SLOPE/w和QUAKE/w兩個(gè)模塊，對(duì)滑坡體進(jìn)行有限元數(shù)值模擬。QUAKE/W模塊用于處理分析系統(tǒng)的動(dòng)態(tài)部分，可用來分析震動(dòng)前后滑坡體內(nèi)部應(yīng)力應(yīng)變的動(dòng)態(tài)變化特征。當(dāng)QUAKE/W計(jì)算完成后，SLOPE/W模塊可調(diào)用Quake/w模塊的計(jì)算結(jié)果來分析地震動(dòng)過程中邊坡的穩(wěn)定性，尋找最危險(xiǎn)滑面并求得到安全系數(shù)。

1.2.1計(jì)算剖面選取

計(jì)算剖面的選取將與滑坡主滑方向一致且剖面具有代表性，能反映滑坡的變形、破壞特征，選取的剖面數(shù)據(jù)基本可靠，能較真實(shí)反映巖土體的情況。根據(jù)以上剖面選取原則，本次穩(wěn)定性定量評(píng)價(jià)選取主剖面線進(jìn)行力學(xué)計(jì)算，計(jì)算剖面圖見圖2。

圖2　耶柿村滑坡力學(xué)計(jì)算剖面圖

1.2.2計(jì)算模型與參數(shù)選擇

在有限元計(jì)算中，網(wǎng)格劃分密度應(yīng)滿足計(jì)算精度的要求，因此對(duì)單元網(wǎng)格進(jìn)行劃分，既要有較高的精度，又要能模擬任意幾何形狀，符合實(shí)際情況。數(shù)值模型受周邊巖體的約束，故模型采用位移約束，底部位移不變，即約束X、Y兩個(gè)方向位移，左右邊界豎向位移不變，只約束Y方向的位移，坡面無約束，為自由平面。坐標(biāo)軸選擇豎直水平方向，Y軸為豎直，向上為正，X軸為水平，指向臨空面為正。數(shù)值模擬物理力學(xué)參數(shù)按表1取值。

數(shù)值模擬過程中采用理想彈塑性模型，屈服準(zhǔn)則采用Mohr-Coufomb強(qiáng)度準(zhǔn)則，屈服函數(shù)選用下列莫爾屈服準(zhǔn)則驗(yàn)算地震荷載作用下邊坡的穩(wěn)定性[14]。

(1)

式中：F為邊坡穩(wěn)定性強(qiáng)度；φ為內(nèi)摩擦角，θσ為應(yīng)力第二洛德角一般取值范圍在- 30°≤θσ≤30°之間。I1為應(yīng)力張量的第一不變量:

I1=σx+σy+σz。

(2)

J2為應(yīng)力偏量的第二不變量：

(3)

地震波施加在垂直邊坡側(cè)的水平方向。如果地震荷載持續(xù)時(shí)間過長，將導(dǎo)致計(jì)算的工程量非常大、而且計(jì)算結(jié)果的后處理也非常麻煩，根據(jù)汶川地震特征及其在西安地區(qū)的影響，本次模擬將地震持時(shí)間定為50 s，定義時(shí)間步長為0.2 s，需計(jì)算100步。

2　結(jié)果分析

根據(jù)汶川地震烈度分布圖，西安地區(qū)主要位于汶川地震Ⅴ～Ⅵ度烈度區(qū)域，因此，為了更好的評(píng)價(jià)地震對(duì)該滑坡的作用過程，主要采用了Ⅵ～Ⅷ度對(duì)應(yīng)的地震加速度進(jìn)行分析。在采用時(shí)程分析法對(duì)邊坡進(jìn)行地震反應(yīng)分析時(shí)，需要輸入地震動(dòng)加速度時(shí)程?？紤]滑坡體地質(zhì)條件，采用修正后的汶川地震波作為動(dòng)力響應(yīng)的輸入波，提取前10 s加速度數(shù)據(jù)進(jìn)行分析，并認(rèn)為在地震到達(dá)時(shí)地震加速度同時(shí)施加到邊坡每個(gè)有限元網(wǎng)格的節(jié)點(diǎn)上。由于地震傳播方向的不規(guī)則性，為了簡(jiǎn)化，本文假定地震加速度作用與坡體臨空面的水平方向一致(圖3)。

圖3　汶川地震波波形圖(假設(shè)加速度峰值為0.1)

首先利用GEOSTUDIO中的SLOPE/w模塊，對(duì)滑坡進(jìn)行自然狀態(tài)下的穩(wěn)定性進(jìn)行分析，邊坡因自重應(yīng)力形成的水平和豎向應(yīng)力分布較均勻，整體處于受壓狀態(tài)，基覆界面處水平應(yīng)力最大值一般在為0.2～0.6 MPa。而剪應(yīng)力整體受力均勻，表層堆積物其值一般在0.05～0.2 MPa。該滑坡在自然狀態(tài)下的穩(wěn)定系數(shù)為1.31，說明該滑坡在自然狀態(tài)下，其整體穩(wěn)定性較好，發(fā)生整體失穩(wěn)的可能性極小(圖4)。

由于滑坡體的阻尼作用，坡體對(duì)地震動(dòng)力響應(yīng)相對(duì)地震波作用有一定的滯后，而坡體各節(jié)點(diǎn)的加速度在地震作用的25 s前后達(dá)到最大，故提取t=24 s時(shí)滑坡在加速度峰值為0.08 g(相當(dāng)于Ⅵ度烈度區(qū))、0.10 g(相當(dāng)于Ⅶ度烈度區(qū))和0.20 g(相當(dāng)于Ⅷ度烈度區(qū))的地震動(dòng)荷載作用下的位移、應(yīng)力和剪應(yīng)變?cè)隽糠植继卣鳌?/p>

地震動(dòng)荷載作用下坡體的位移、應(yīng)力和剪應(yīng)變?cè)隽糠植继卣鳎瑥睦碚撋虾芎玫卦佻F(xiàn)了邊坡從靜態(tài)進(jìn)入動(dòng)態(tài)所發(fā)生的變化，充分說明地震對(duì)邊坡穩(wěn)定性的影響極其大。

坡體在加速度峰值為0.08 g的地震動(dòng)荷載作用后，最大水平位移出現(xiàn)在坡體上部，為6.55 cm，水平總應(yīng)力表現(xiàn)沿坡體深度呈明顯增大趨勢(shì)，但局部不均勻，基覆界面處水平應(yīng)力最大值為2～2.5 MPa，滑坡體基覆界面處的最大剪應(yīng)力最大為0.6 MPa。相應(yīng)剪應(yīng)變主要發(fā)生滑坡體后緣下部滑體與基巖交界內(nèi)，最大剪應(yīng)變值為1.2×10-3(圖5～圖6)，坡體整體處于穩(wěn)定狀態(tài)，計(jì)算得出此時(shí)的安全系數(shù)為1.22，由此可見在此地震動(dòng)荷載作用下，魏寨鄉(xiāng)耶柿村滑坡整體上處于較高的穩(wěn)定狀態(tài)。

圖4　滑體天然狀態(tài)下自重應(yīng)力水平x向及最大剪應(yīng)力云圖

圖5　地震動(dòng)荷載作用下滑坡體剖面水平位移云圖

圖6　地震動(dòng)荷載作用下滑坡體剖面圖水平應(yīng)力云圖

圖7　地震動(dòng)荷載作用下剖面最大剪應(yīng)力云圖

圖8　地震動(dòng)荷載作用下剖面剪應(yīng)變?cè)茍D

圖9　地震動(dòng)荷載作用下邊坡穩(wěn)定性系數(shù)

在加速度峰值為0.10 g的地震動(dòng)荷載作用下(相當(dāng)于Ⅶ度烈度區(qū))，從水平位移云圖看，位移較大的地方主要發(fā)生在滑體表層，滑坡后緣基巖部位變形不明顯，坡體表層覆蓋物出現(xiàn)最大水平位移為8.4 cm。由圖7可以看出，水平總應(yīng)力與剪應(yīng)力隨著峰值加速度增加增加，地震時(shí)坡體水平應(yīng)力總體表現(xiàn)沿坡體深度呈明顯增大趨勢(shì)，但局部不均勻，水平總應(yīng)力最大值出現(xiàn)在基覆界面處，為0.2～1.6 MPa；而坡體最大剪應(yīng)力出現(xiàn)在坡體與基巖交界內(nèi)，值為0.8 MPa(圖7)。由圖可以看出，一條剪應(yīng)變集中帶已明顯出現(xiàn)，但剪應(yīng)變僅出現(xiàn)在滑坡體后緣及坡腳處，且擴(kuò)展范圍較小，最大剪應(yīng)變值為1.8×10-3。由此可知，坡體在此地震破壞作用下具有小規(guī)模失穩(wěn)的可能，主要發(fā)生在坡體前緣，計(jì)算得出此時(shí)的安全系數(shù)為1.12。由此可見，在此地震動(dòng)荷載作用下，魏寨鄉(xiāng)耶柿村滑坡整體上處于穩(wěn)定狀態(tài)，只有局部發(fā)生失穩(wěn)的可能。

在加速度峰值為0.20 g的地震動(dòng)荷載作用下(相當(dāng)于Ⅷ度烈度區(qū))，隨著動(dòng)峰加速度的增加，水平位移呈明顯上升，最大水平位移為13.1 cm，主要發(fā)生在滑體表層(圖8)。由于基巖傾角相對(duì)較緩，上覆堆積體也較薄，此部分堆積體在地震動(dòng)荷載作用下，剪應(yīng)變集中帶依然沒有貫通，最大剪應(yīng)變也達(dá)到2.9×10-3。但是上部坡體的水平應(yīng)力與剪應(yīng)變均具有明顯增大的趨勢(shì)，水平總應(yīng)力最大值出現(xiàn)在基覆界面處，為 2.8 MPa，坡體剪應(yīng)力值增大為1.3 MPa；由此可見，在此地震動(dòng)荷載作用下，滑坡上部堆積體邊坡位移較大，滑動(dòng)面具有整體貫通的趨勢(shì)，計(jì)算得出此時(shí)的安全系數(shù)為1.03。在該震動(dòng)作用下，坡體處于欠穩(wěn)定的狀態(tài)，具有發(fā)生整體失穩(wěn)的可能。

通過對(duì)比分析該滑坡在不同狀態(tài)下的穩(wěn)定系數(shù)可以看出，在天然狀態(tài)下該滑坡整體穩(wěn)定性較好，很難發(fā)生整體失穩(wěn)，但是隨著地震峰值加速度的增加，該滑坡的穩(wěn)定系數(shù)出現(xiàn)急劇下降的趨勢(shì)，從天然狀態(tài)下的1.31下降到1.03的欠穩(wěn)定狀態(tài)(圖9)。尤其是隨著地震峰值加速度的增加，邊坡存在部分失穩(wěn)的可能。

3　耶柿村滑坡防治措施

根據(jù)上述數(shù)值模擬結(jié)果，該滑坡體在地震動(dòng)作用下均發(fā)生位移，尤其是當(dāng)?shù)卣鹆叶却笥冖葧r(shí)，該滑坡體具有復(fù)活失穩(wěn)的可能，由此提出該滑坡體的防治措施。

(1)監(jiān)測(cè)措施：首先采取監(jiān)測(cè)措施，以便隨時(shí)掌握滑坡的變形情況，提出預(yù)警預(yù)報(bào)，并適時(shí)啟動(dòng)“耶柿村滑坡防、搶、撤方案”或采取工程治理措施。

(2)避讓措施：由于耶柿村滑坡穩(wěn)定性較差，威脅人口多，危險(xiǎn)性大。為了確保耶柿村人民生命財(cái)產(chǎn)安全，最大限度地減少人員傷亡和經(jīng)濟(jì)損失，應(yīng)根據(jù)經(jīng)濟(jì)條件和危險(xiǎn)性分區(qū)逐批將村民進(jìn)行搬遷；在提前預(yù)警和突發(fā)情況下按“耶柿村滑坡防、搶、撤方案”進(jìn)行緊急撤離。

(3)工程措施：針對(duì)該滑坡的特點(diǎn)，為防止滑坡進(jìn)一步活動(dòng)，首先應(yīng)采取簡(jiǎn)易工程措施，如截排水溝，對(duì)裂縫加以填埋等簡(jiǎn)易工程措施。為了確保居住在滑坡區(qū)的人民生命財(cái)產(chǎn)安全，建議對(duì)該滑坡進(jìn)行詳細(xì)的勘察工作，進(jìn)而進(jìn)行滑坡防治工程設(shè)計(jì)，對(duì)滑坡進(jìn)行相應(yīng)的工程治理。

4　結(jié)論

以魏寨鄉(xiāng)耶柿村滑坡為例，根據(jù)工程地質(zhì)條件、構(gòu)造發(fā)育特征，通過運(yùn)用數(shù)值模擬動(dòng)力分析方法，分析邊坡的動(dòng)力變化特性與動(dòng)力響應(yīng)規(guī)律。利用Geostudio計(jì)算安全系數(shù)分析穩(wěn)定性。隨著地震動(dòng)峰值加速度的增大邊坡穩(wěn)定性逐漸降低，安全系數(shù)降低到1.06。計(jì)算結(jié)果表明，魏寨鄉(xiāng)耶柿村滑坡僅在大于Ⅷ度地震烈度作用下，滑坡體處于欠穩(wěn)定狀態(tài)，具有整體發(fā)生失穩(wěn)的可能；在Ⅶ度地震烈度作用下，在滑坡體前緣具有發(fā)生小規(guī)模失穩(wěn)的可能；在Ⅵ度地震烈度作用下，該滑坡體整體穩(wěn)定性系數(shù)為1.22，穩(wěn)定性較高。但是該滑坡體在地震作用下，邊坡內(nèi)部均出現(xiàn)位移，對(duì)坡體內(nèi)部具有一定的損傷，由此提出了相應(yīng)的防治對(duì)策。本文只是數(shù)值模擬，具體的滑坡狀態(tài)和防治措施，還需進(jìn)一步深入研究。