丁文潔，馬祥配，彭永和，別鳳華，楊超，解露茜

(山東省第八地質(zhì)礦產(chǎn)勘查院，山東地礦局有色金屬礦找礦與資源評(píng)價(jià)重點(diǎn)實(shí)驗(yàn)室，日照地質(zhì)地理信息大數(shù)據(jù)研究院，日照市土地質(zhì)量評(píng)價(jià)與污染修復(fù)重點(diǎn)實(shí)驗(yàn)室，山東 日照 276826)

0 引言

隨著我國(guó)經(jīng)濟(jì)建設(shè)的加快，基礎(chǔ)設(shè)施建設(shè)不斷加強(qiáng)。在公路建設(shè)中由于受道路運(yùn)輸?shù)葪l件的限制，不可避免的會(huì)產(chǎn)生大量棄土棄渣，因而，在公路建設(shè)時(shí)必然要設(shè)置一定的棄土場(chǎng)來(lái)滿足棄土處置的需要[1]。

公路棄土場(chǎng)由于其欠固結(jié)、多孔隙、非飽和等結(jié)構(gòu)特點(diǎn)[2]，在降雨及動(dòng)力荷載的影響下極易發(fā)生失穩(wěn)滑坡[3]。研究棄土場(chǎng)邊坡的穩(wěn)定性，可以盡可能減小棄土場(chǎng)的潛在破壞性，對(duì)避免由于棄渣產(chǎn)生次生災(zāi)害、減少水土流失具有重要意義。在棄土場(chǎng)的選址、支擋措施等設(shè)計(jì)中就可有針對(duì)性地采取相應(yīng)的措施，達(dá)到安全、經(jīng)濟(jì)和有效的目的。

目前，我國(guó)針對(duì)公路棄土場(chǎng)開(kāi)展的研究很少。譚鵬[4]分析了棄土場(chǎng)失穩(wěn)破壞機(jī)理并研究提出了棄土場(chǎng)穩(wěn)定性分析方法，沈明祥[5]通過(guò)現(xiàn)場(chǎng)調(diào)查和計(jì)算分析評(píng)估了貴州省六盤水至威寧高速公路棄土場(chǎng)穩(wěn)定性，鄭鶴丹[6]運(yùn)用Geo Studio軟件對(duì)棄土體邊坡穩(wěn)定性進(jìn)行了評(píng)估，分析正常工況和非常工況條件下棄土體邊坡的穩(wěn)定性，李研[7]對(duì)高速公路棄土場(chǎng)的處置方案進(jìn)行了研究。以上文獻(xiàn)從分析機(jī)理、現(xiàn)場(chǎng)調(diào)查、計(jì)算分析、風(fēng)險(xiǎn)評(píng)價(jià)、處置方案設(shè)計(jì)等方面研究了高速公路棄土場(chǎng)的相關(guān)問(wèn)題，但均缺乏穩(wěn)定影響因素敏感性的分析。因此，在實(shí)際工作中，有必要結(jié)合具體工程實(shí)例對(duì)高速公路棄土場(chǎng)進(jìn)行穩(wěn)定影響因素敏感性分析，預(yù)防地質(zhì)災(zāi)害的發(fā)生。

公路棄土場(chǎng)在填筑過(guò)程中，需要根據(jù)棄土場(chǎng)的用途、穩(wěn)定性，以及綜合經(jīng)濟(jì)性等多種因素確定壓實(shí)控制標(biāo)準(zhǔn)。本文以日照市嵐山區(qū)某棄土場(chǎng)為例，通過(guò)MIDAS GTS NX進(jìn)行數(shù)據(jù)模擬，在不同壓實(shí)控制標(biāo)準(zhǔn)的情況下，歸類分析高度和坡度對(duì)邊坡穩(wěn)定性系數(shù)的變化規(guī)律，從而找到邊坡失穩(wěn)的敏感因素。

1 敏感性分析原理

影響邊坡穩(wěn)定的因素有坡角、坡高、黏聚力、內(nèi)摩擦角、含水率、重度以及地震、降雨等，各因素對(duì)邊坡失穩(wěn)的破壞程度是不同的，因此,通過(guò)邊坡穩(wěn)定影響因素的敏感性分析，可以找出影響棄土場(chǎng)邊坡穩(wěn)定的敏感因素，分析影響邊坡穩(wěn)定性各因素與邊坡安全系數(shù)之間的相關(guān)性，即分析各因素的變化對(duì)于邊坡安全系數(shù)的影響程度[8]。

敏感性分析是系統(tǒng)分析中對(duì)系統(tǒng)穩(wěn)定性的一種不確定性分析方法[9-10]。

本文通過(guò)MIDAS GTS NX進(jìn)行模擬，計(jì)算出不同條件下的安全系數(shù)，然后選取任意一個(gè)影響因素在可能范圍內(nèi)變化，且保持其他影響因素不變來(lái)計(jì)算安全系數(shù)，最后計(jì)算各影響因素的敏感性。

2 棄土場(chǎng)現(xiàn)狀

日照市嵐山區(qū)某棄土場(chǎng)，坡腳距實(shí)驗(yàn)小學(xué)邊坡坡頂約10m。整個(gè)棄土場(chǎng)坡面平緩，坡面未見(jiàn)滑坡，坡頂未見(jiàn)裂縫、坍塌現(xiàn)象，整個(gè)坡面完整性好。

棄土場(chǎng)堆土巖性為素填土，最大厚度約6.0m，中密，濕，成分以回填風(fēng)化巖碎屑為主；其下為殘積成因砂質(zhì)黏性土，硬塑狀；底部為晉寧期花崗片麻巖，強(qiáng)風(fēng)化，變晶結(jié)構(gòu)，片麻狀構(gòu)造，主要成分為石英、長(zhǎng)石及其他暗色礦物，巖體呈致密的砂礫狀及角礫狀(圖1)。坡頂植被不發(fā)育。

1—填土；2—強(qiáng)風(fēng)化花崗巖圖1 棄土場(chǎng)剖面圖

3 建立基本模型

3.1 建模

本文根據(jù)實(shí)際的邊坡坡面形態(tài)、各巖層巖性及其物理力學(xué)特征構(gòu)建模型[12-14]。采用MIDAS-GTSNX軟件中的有限元強(qiáng)度折減法(SRM)[15-18]算出邊坡的等效應(yīng)變塑性區(qū)及安全系數(shù)，以此分析邊坡穩(wěn)定性。

模型分為兩個(gè)部分，分別是棄土場(chǎng)基底區(qū)A和棄土區(qū)B(圖2、圖3、圖4)，圖中單位為m，計(jì)算時(shí)線網(wǎng)格控制間距0.5m。

約束條件:模型左右兩側(cè)為水平約束,底部固定,上部及坡面為自由邊界。

本構(gòu)模型:根據(jù)土體的應(yīng)力應(yīng)變特性,選取本構(gòu)模型為Mohr-Coulomb模型。

初始應(yīng)力:僅考慮自重應(yīng)力產(chǎn)生的應(yīng)力場(chǎng)。

計(jì)算收斂標(biāo)準(zhǔn):收斂誤差小于定值10-6,最大迭代次數(shù)1000次，計(jì)算終止。重力加速度選取9.80665m/s2，幾何誤差0.0001。

圖2 一級(jí)邊坡二維模型(A)二級(jí)邊坡二維模型(B)

圖3 一級(jí)邊坡剖面等效應(yīng)變?cè)茍D

圖4 二級(jí)邊坡剖面等效應(yīng)變?cè)茍D

3.2 軟土強(qiáng)度參數(shù)試驗(yàn)

堆放于棄土場(chǎng)的土主要為坡積粉質(zhì)黏土，開(kāi)展了直剪試驗(yàn)[19]，得到了不同壓實(shí)度下的黏聚力、內(nèi)摩擦角值和容重值(表1)。

表1 直剪試驗(yàn)結(jié)果一覽表

4 模擬計(jì)算

根據(jù)前人[20]對(duì)工程邊坡進(jìn)行過(guò)的研究，整理總結(jié)并通過(guò)現(xiàn)場(chǎng)工程地質(zhì)勘查分析，影響邊坡穩(wěn)定性的因素有很多[21]，主要影響因素有坡腳角度、邊坡高度、巖(土)體的抗剪強(qiáng)度性指標(biāo)內(nèi)摩擦角和巖(土)體的黏聚力[22]。壓實(shí)度決定了巖(土)體的抗剪強(qiáng)度性指標(biāo)內(nèi)摩擦角和巖(土)體的黏聚力，本文主要分析在不同壓實(shí)控制標(biāo)準(zhǔn)的情況下，研究坡腳角度、邊坡高度對(duì)邊坡穩(wěn)定性系數(shù)的影響。

本文采用控制變量法[23]來(lái)探究邊坡穩(wěn)定性，運(yùn)用MIDAS GTS NX進(jìn)行安全系數(shù)計(jì)算。假設(shè)該邊坡為均質(zhì)土層，在不同壓實(shí)控制標(biāo)準(zhǔn)的情況下，選取不同的坡角β、坡高H，來(lái)獲得不同的安全系數(shù)，并對(duì)獲得的穩(wěn)定系數(shù)進(jìn)行分析。

(1)壓實(shí)度為90%時(shí)，黏聚力為28kPa，內(nèi)摩擦角為14.6°，坡高分別選取5m、6m、7m、8m、9m、10m，坡比分別選取1∶0.5、1∶0.75、1∶1.00、1∶1.25、1∶1.50，得到相應(yīng)的安全系數(shù),詳見(jiàn)表2、表3。

表2 90%壓實(shí)度下不同坡高、坡比對(duì)應(yīng)的安全系數(shù)(一級(jí)坡)

表3 90%壓實(shí)度下不同坡高、坡比對(duì)應(yīng)的安全系數(shù)(二級(jí)坡)

(2)壓實(shí)度為92%時(shí)，黏聚力為33kPa，內(nèi)摩擦角為15.3°，坡高分別選取5m、6m、7m、8m、9m、10m，坡比分別選取1∶0.5、1∶0.75、1∶1.00、1∶1.25、1∶1.50，得到相應(yīng)的安全系數(shù),詳見(jiàn)表4、表5。

表4 92%壓實(shí)度下不同坡高、坡比對(duì)應(yīng)的安全系數(shù)(一級(jí)坡)

表5 92%壓實(shí)度下不同坡高、坡比對(duì)應(yīng)的安全系數(shù)(二級(jí)坡)

(3)壓實(shí)度為95%時(shí)，黏聚力為41kPa，內(nèi)摩擦角為15.5°，坡高分別選取5m、6m、7m、8m、9m、10m，坡比分別選取1∶0.5、1∶0.75、1∶1.00、1∶1.25、1∶1.50，得到相應(yīng)的安全系數(shù),詳見(jiàn)表6、表7。

表6 95%壓實(shí)度下不同坡高、坡比對(duì)應(yīng)的安全系數(shù)(一級(jí)坡)

表7 95%壓實(shí)度下不同坡高、坡比對(duì)應(yīng)的安全系數(shù)(二級(jí)坡)

5 邊坡穩(wěn)定性影響因素的敏感分析

公路棄土場(chǎng)在填筑過(guò)程中，已經(jīng)根據(jù)實(shí)際需要確定了壓實(shí)控制標(biāo)準(zhǔn)。由于不同壓實(shí)控制標(biāo)準(zhǔn)決定了巖(土)體的內(nèi)摩擦角和黏聚力，所以本次工作選取在不同壓實(shí)控制標(biāo)準(zhǔn)的情況下，改變坡高H、坡角β，進(jìn)行影響因素的敏感分析。

5.1 邊坡高度的敏感性分析

利用數(shù)值模擬的結(jié)果(表2—表7),繪制坡高H與安全系數(shù)的關(guān)系圖(圖5)。通過(guò)圖5知,在相同的坡比下,坡高與安全系數(shù)呈向上凹的拋物線關(guān)系,且隨著坡高的增加,邊坡安全系數(shù)在不斷減小。坡高在5～7m范圍逐漸增大時(shí)，邊坡安全系數(shù)降低很快；在坡高在7～10m之間增大時(shí)，安全系數(shù)降低速率相較于穩(wěn)定性系數(shù)在5～7m范圍內(nèi)的變化要小。

a—90%壓實(shí)度一級(jí)坡坡高與安全系數(shù)關(guān)系圖;b—90%壓實(shí)度二級(jí)坡坡高與安全系數(shù)關(guān)系圖;c—92%壓實(shí)度一級(jí)坡坡高與安全系數(shù)關(guān)系圖;d—92%壓實(shí)度二級(jí)坡坡高與安全系數(shù)關(guān)系圖;e—95%壓實(shí)度一級(jí)坡坡高與安全系數(shù)關(guān)系圖;f—95%壓實(shí)度二級(jí)坡坡高與安全系數(shù)關(guān)系圖圖5 不同類別下坡高與安全系數(shù)關(guān)系圖

根據(jù)敏感度原理，計(jì)算出不同類別下棄土場(chǎng)坡高的敏感度，例如：90%壓實(shí)度一級(jí)坡，坡比為1∶0.5，當(dāng)高度從7m減少到5m(降低約為28.57%)時(shí)，安全系數(shù)由1.42增加至1.89(增加約為33.10%)，棄土場(chǎng)坡高的敏感度約為1.158，而高度從7m增加到10m(增加約為42.86%)時(shí)，安全系數(shù)由1.42減小至1.08(降低約為23.94%)，棄土場(chǎng)坡高的的敏感度約為0.559。各個(gè)類別下的敏感度計(jì)算如表8。

(1)棄土場(chǎng)邊坡高度下降敏感度數(shù)值遠(yuǎn)大于邊坡高度上升敏感度數(shù)值，由此可以得出：棄土場(chǎng)邊坡高度下降對(duì)安全系數(shù)影響比棄土場(chǎng)邊坡高度增加對(duì)安全系數(shù)影響敏感。

(2)隨著坡比的減小，無(wú)論邊坡高度下降還是增加，坡高對(duì)安全系數(shù)的影響有減小的趨勢(shì),即說(shuō)明隨著坡比的減小,邊坡高度的變化對(duì)穩(wěn)定性的影響在減小。

(3)相同條件下，壓實(shí)度越大，無(wú)論坡高減小還是增加，坡高對(duì)安全系數(shù)的敏感度有下降的趨勢(shì)，說(shuō)明壓實(shí)度越大，邊坡高度對(duì)安全系數(shù)的影響在減小。

表8 棄土場(chǎng)坡高敏感度統(tǒng)計(jì)表

5.2 邊坡坡比的敏感性分析

利用數(shù)值模擬的結(jié)果(表2—表7),繪制坡比與安全系數(shù)的關(guān)系圖(圖6)。通過(guò)圖6知,在相同的高度下,坡比與安全系數(shù)呈一定的線性關(guān)系,且隨著坡比由1∶0.5變化到1∶1.5,邊坡安全系數(shù)在不斷增大。

(1)棄土場(chǎng)坡比減小敏感度數(shù)值大于坡比增加敏感度數(shù)值，由此可以得出，棄土場(chǎng)坡比減小對(duì)安全系數(shù)影響比棄土場(chǎng)坡比增加對(duì)安全系數(shù)影響大。

(2)隨著高度的增加，無(wú)論坡比增加還是減小，坡比對(duì)安全系數(shù)的敏感度呈增加的趨勢(shì)，說(shuō)明隨著高度的增加，坡比對(duì)安全系數(shù)的影響在增大。

根據(jù)敏感度原理，計(jì)算出不同類別下棄土場(chǎng)坡比的敏感度，例如：90%壓實(shí)度一級(jí)坡，坡高為5m，當(dāng)坡比從1∶1.00增加到1∶0.50(增加約為100%)時(shí)，安全系數(shù)由2.40減少至1.89(減少約為21.25%)，棄土場(chǎng)坡高的敏感度約為0.213，而坡比從1∶1.00減小到1∶1.50(減小約為33.33%)時(shí)，安全系數(shù)由2.40增加至2.84(增加約為18.33%)，棄土場(chǎng)坡高的敏感度約為0.550。各個(gè)類別下的敏感度計(jì)算如表9。

(3)在相同高度下，壓實(shí)度越大，無(wú)論坡比增加還是減小，坡比對(duì)安全系數(shù)的敏感度減小，說(shuō)明相同高度下，壓實(shí)度越大，坡比對(duì)安全系數(shù)的影響減小。

a—90%壓實(shí)度一級(jí)坡坡比與安全系數(shù)關(guān)系圖;b—90%壓實(shí)度二級(jí)坡坡比與安全系數(shù)關(guān)系圖;c—92%壓實(shí)度一級(jí)坡坡比與安全系數(shù)關(guān)系圖;d—92%壓實(shí)度二級(jí)坡坡比與安全系數(shù)關(guān)系圖;e—95%壓實(shí)度一級(jí)坡坡比與安全系數(shù)關(guān)系圖;f—95%壓實(shí)度二級(jí)坡坡比與安全系數(shù)關(guān)系圖圖6 不同類別下坡比與安全系數(shù)關(guān)系圖

表9 棄土場(chǎng)坡比敏感度統(tǒng)計(jì)表

6 結(jié)論

利用MIDAS GTS NX軟件,對(duì)壓實(shí)度分別為90%、92%、95%，高度為5～10m,以及坡角為1∶1.50～1∶0.50的一級(jí)坡和二級(jí)坡進(jìn)行了模擬,系統(tǒng)的對(duì)邊坡穩(wěn)定性進(jìn)行了分析。利用控制變量法,分析了在坡高和坡比單獨(dú)變化下邊坡的穩(wěn)定性。利用敏感度,分析了在不同壓實(shí)控制標(biāo)準(zhǔn)下,坡高和坡比對(duì)安全系數(shù)的影響。以下為本文的主要結(jié)論:

(1)棄土場(chǎng)邊坡高度下降比棄土場(chǎng)邊坡高度增加對(duì)安全系數(shù)影響敏感；棄土場(chǎng)坡比減小比棄土場(chǎng)坡比增加對(duì)安全系數(shù)影響敏感。

(2)隨著坡比的減小,邊坡高度的變化對(duì)穩(wěn)定性的影響在不斷減??；隨著高度的增加，坡比的變化對(duì)安全系數(shù)的影響增大。

(3)相同條件下，壓實(shí)度越大，邊坡高度和坡比對(duì)安全系數(shù)的影響越小。