朱博良,王延壽,王俊卿
(甘肅省有色金屬地質(zhì)勘查局 天水礦產(chǎn)勘查院,甘肅 天水 741020)
我國黃土分布面積之廣、厚度之大,堪稱世界之最。在偌大的黃土高原,黃土邊坡既是人類生存的重要環(huán)境,又是工程建設(shè)的重要載體。氣候變化、環(huán)境破壞和人類活動,導(dǎo)致黃土邊坡出現(xiàn)了滑坡、崩塌等地質(zhì)災(zāi)害。近年來,隨著“西部大開發(fā)”和“一帶一路”經(jīng)濟(jì)建設(shè),黃土高原地區(qū)大量基礎(chǔ)性建設(shè)工程上馬,不可避免地開挖坡腳,形成了較多的人工邊坡,使得黃土邊坡的穩(wěn)定性問題更為突出。
導(dǎo)致黃土邊坡失穩(wěn)的因素很多,黃土的微結(jié)構(gòu)特征、邊坡高度、降雨、地震、凍融、灌溉、開挖坡腳等都是黃土邊坡失穩(wěn)的重要因素,尤其是較大的坡度,為滑坡和崩塌的產(chǎn)生提供了必要條件。本研究以位于黃土高原地區(qū)的天水市麥積區(qū)某黃土邊坡為例,設(shè)計了不同坡度的76例黃土邊坡模型,采用基于極限平衡條分法的瑞典法、簡化畢肖普法和簡化簡布法在天然狀態(tài)下進(jìn)行邊坡穩(wěn)定性系數(shù)計算,討論黃土邊坡穩(wěn)定性系數(shù)隨坡度的變化規(guī)律,并建立預(yù)測模型,用強度折減法驗證模型的合理性。
本研究以位于黃土高原地區(qū)的天水市麥積區(qū)某黃土邊坡作為研究對象。該邊坡坡體由第四系上更新統(tǒng)馬蘭黃土構(gòu)成,坡高6.0 m,天然容重18.5 kN/m3,飽和容重 19.0 kN/m3,黏聚力20.0 kPa,內(nèi)摩擦角24.0°,坡度15°~90°,剖面示意見圖1。
保持計算模型中其他尺寸不變,只以間隔1°變換坡度,得到坡度介于15°~90°的76例計算模型,利用理正巖土邊坡穩(wěn)定性分析系統(tǒng),采用瑞典法、簡化畢肖普法和簡化簡布法對不同坡度的76例黃土邊坡模型進(jìn)行穩(wěn)定性系數(shù)計算,研究邊坡穩(wěn)定性系數(shù)與坡度的關(guān)系。
圖1 計算模型示意(單位:m)
建立黃土邊坡穩(wěn)定性系數(shù)與坡度的關(guān)系式為
y=kf(x)
(1)
式中:y為除坡度外其他影響穩(wěn)定性系數(shù)的因子為恒量時的邊坡穩(wěn)定性系數(shù);x為坡度,(°);k為常數(shù)。
采用瑞典法、簡化畢肖普法和簡化簡布法對不同坡度的76例黃土邊坡模型進(jìn)行穩(wěn)定性系數(shù)計算,從穩(wěn)定性系數(shù)隨坡度的變化關(guān)系(圖2)中可以看出,采用這3種方法計算得到的穩(wěn)定性系數(shù)均隨坡度的增大而減小,即:坡度15°時穩(wěn)定性系數(shù)為最大,分別為3.464、3.706、3.718;坡度90°時穩(wěn)定性系數(shù)為最小,分別為1.340、1.244、1.228。具體來說,穩(wěn)定性系數(shù)隨坡度變化的過程可分為3個階段:當(dāng)坡度為15°~55°時,穩(wěn)定性系數(shù)隨坡度增加下降非常明顯,稱為陡降段;當(dāng)坡度>66°時,穩(wěn)定性系數(shù)隨坡度增加變化不明顯,趨于平緩,稱為平緩段;當(dāng)坡度介于55°~66°時,穩(wěn)定性系數(shù)隨坡度增加變化較明顯,緩慢下降,稱為過渡段。
圖2 穩(wěn)定性系數(shù)與坡度關(guān)系曲線
根據(jù)文獻(xiàn)[1]的研究成果,采用瑞典法、簡化畢肖普法和簡化簡布法計算得到的穩(wěn)定性系數(shù)略小于強度折減法計算結(jié)果,這是因為強度折減法考慮了土體的非線性本構(gòu)關(guān)系,以及應(yīng)力與應(yīng)變的關(guān)系。綜合考慮,在采用瑞典法、簡化畢肖普法和簡化簡布法計算結(jié)果不一致時應(yīng)選用最大值。結(jié)合圖2,當(dāng)坡度為15~55°時,應(yīng)選用簡化畢肖普法或簡化簡布法;當(dāng)坡度≥55°時,應(yīng)選用瑞典法。
用穩(wěn)定性系數(shù)與坡度的關(guān)系建立黃土邊坡穩(wěn)定性系數(shù)預(yù)測模型。從3種方法得到的穩(wěn)定性系數(shù)隨坡度變化的擬合曲線(圖3)可以看出,3條曲線的擬合度均極高,決定系數(shù)分別為0.999 3、0.998 5、0.998 0,且各曲線均呈現(xiàn)穩(wěn)定性系數(shù)隨坡度的增加而降低的變化趨勢,其函數(shù)關(guān)系式分別為
圖3 穩(wěn)定性系數(shù)隨坡度變化的擬合曲線
y=14.525x-0.534
(2)
y=19.809x-0.612
(3)
y=20.291x-0.619
(4)
根據(jù)上文研究成果,當(dāng)坡度≥55°時,可選用式(2)進(jìn)行預(yù)測;當(dāng)坡度<55°時,可選用式(3)或式(4)進(jìn)行預(yù)測。考慮到用簡化畢肖普法和簡化簡布法計算得到的穩(wěn)定性系數(shù)基本重合,將二者取平均值,可得
y= 20.048x-0.616
(5)
將式(2)、(5)聯(lián)合,得到黃土邊坡穩(wěn)定性系數(shù)隨坡度變化的預(yù)測模型為
(6)
根據(jù)上式,不同坡度可選用相應(yīng)的函數(shù)式進(jìn)行計算。當(dāng)黃土邊坡其他參數(shù)與本算例參數(shù)基本一致且已知坡度時,可利用預(yù)測模型進(jìn)行計算。下面擬通過實例來驗證該預(yù)測模型的合理性。
利用文獻(xiàn)[1]中某黃土邊坡工程實例對預(yù)測模型的合理性進(jìn)行驗證。該黃土邊坡其他參數(shù)與本算例參數(shù)基本一致,且已知邊坡坡度為62°。用邁達(dá)斯軟件中巖土與隧道分析系統(tǒng)(GTS)和理正巖土邊坡穩(wěn)定性分析系統(tǒng)兩種方法共同驗證。由基于有限元法的邁達(dá)斯GTS求得的穩(wěn)定性系數(shù)為1.825,用理正系統(tǒng)基于極限平衡條分法的瑞典法、簡化畢肖普法和簡化簡布法得到的穩(wěn)定性系數(shù)分別為1.601、1.600、1.599,平均值為1.600。強度折減法計算結(jié)果稍大,而極限平衡法計算結(jié)果略偏于安全。根據(jù)式(6),當(dāng)坡度x=62°時,得到的穩(wěn)定性系數(shù)為1.603,數(shù)值介于上述兩種方法計算的結(jié)果之間,證明該預(yù)測模型具有一定的合理性。
通過研究76例不同坡度的黃土邊坡模型穩(wěn)定性系數(shù)隨坡度的變化規(guī)律及其預(yù)測模型,得到以下結(jié)論。
(1)黃土邊坡穩(wěn)定性系數(shù)隨坡度的增大而減小,該過程大致可分為3個階段,即坡度<55°時的陡降段、坡度>66°時的平緩段、坡度介于55~66°時的過渡段。
(2)預(yù)測模型可用冪函數(shù)表達(dá),擬合曲線的決定系數(shù)較高,通過驗證,該預(yù)測模型具有合理性,可用于計算具有相似條件的不同坡度的黃土邊坡穩(wěn)定性系數(shù)。