基于GeoStudio灌木修復(fù)下尾礦壩穩(wěn)定性分析

侯永莉，曹明杰，郝 喆,*，周素航

(1.遼寧有色勘察研究院，遼寧 沈陽(yáng) 110013；2.遼寧大學(xué) 環(huán)境學(xué)院，遼寧 沈陽(yáng) 110036)

0 引言

隨著對(duì)礦產(chǎn)資源需求的日益增加，礦山開(kāi)采逐漸過(guò)度，據(jù)統(tǒng)計(jì)我國(guó)每年尾礦產(chǎn)量超過(guò)6億噸，堆存量已超過(guò)200億噸，現(xiàn)存尾礦超過(guò)2萬(wàn)余座[1-4]。隨著礦山生產(chǎn)的進(jìn)行，尾礦容量隨之增大，進(jìn)而導(dǎo)致尾礦庫(kù)穩(wěn)定性極差，尾礦庫(kù)一旦失穩(wěn)發(fā)生滑坡將對(duì)礦區(qū)及周圍居民造成極大損害，造成難以估量的損失[5-8]。因此尾礦庫(kù)穩(wěn)定性及其安全系數(shù)成為眾多學(xué)者研究的課題，國(guó)內(nèi)外眾多學(xué)者對(duì)尾礦壩穩(wěn)定性及根系固坡效果做了大量研究。李強(qiáng)等研究了基于流固耦合理論尾礦壩失穩(wěn)特性及穩(wěn)定性分析[9]；陳建宏等研究了基于 PCA-BP 神經(jīng)網(wǎng)絡(luò)的尾礦庫(kù)壩體穩(wěn)定性分析，并得出結(jié)論可以使用 PCA-BP 神經(jīng)網(wǎng)絡(luò)模型分析壩體穩(wěn)定性[10]。董紅娟等研究基于3DEC的尾礦壩邊坡監(jiān)測(cè)[11]。李嫣綺等基于強(qiáng)度折減法對(duì)尾礦壩穩(wěn)定性數(shù)值進(jìn)行模擬分析，發(fā)現(xiàn)分析結(jié)果與瑞典圓弧法計(jì)算結(jié)果一致[12]。張年龍使用GPS技術(shù)對(duì)尾礦庫(kù)進(jìn)行在線監(jiān)測(cè)[13]。但鮮有學(xué)者研究植物生長(zhǎng)年限對(duì)尾礦壩安全系數(shù)的影響，因此研究開(kāi)展基于GeoStudio灌木修復(fù)下尾礦壩穩(wěn)定性分析，從尾礦壩采集不同年限的修復(fù)植物及根系土，進(jìn)行根土復(fù)合剪切實(shí)驗(yàn)，進(jìn)一步建立尾礦庫(kù)模型，開(kāi)展GeoStudio模擬安全系數(shù)計(jì)算。本研究旨在為提高尾礦壩安全系數(shù)提供理論依據(jù)。

1 研究方法

1.1 采樣點(diǎn)的確定

尾礦砂與修復(fù)灌木均采集自本溪歪頭山尾礦庫(kù)尾礦壩外坡，歪頭山尾礦庫(kù)共有6級(jí)臺(tái)階，修復(fù)灌木為沙棘，本次實(shí)驗(yàn)采樣根據(jù)歪頭山尾礦庫(kù)生態(tài)修復(fù)的實(shí)際效果，確定對(duì)修復(fù)年限分別為1 a、4 a、10 a的臺(tái)階作為實(shí)驗(yàn)對(duì)象進(jìn)行采樣。

1.2 采樣方式

修復(fù)植物沙棘需連根取出，并保證根部完整，沙棘取出后迅速將根部包裹以防水分蒸干[14]，所取沙棘為其所在臺(tái)階的平均大小，不能過(guò)大或過(guò)小。對(duì)應(yīng)在所取沙棘處的根部取擾動(dòng)土樣，使用鋁盒裝好后編號(hào)帶回實(shí)驗(yàn)室。沙棘與根際土每個(gè)采樣點(diǎn)取3個(gè)重復(fù)，采集樣品見(jiàn)圖1。

1.3 實(shí)驗(yàn)參數(shù)的測(cè)定

取樣現(xiàn)場(chǎng)進(jìn)行灌水試驗(yàn)，測(cè)定土壤天然密度，利用環(huán)刀取原狀樣，測(cè)定含水率和孔隙比，用來(lái)計(jì)算土的天然重度。

進(jìn)行根土復(fù)合剪切實(shí)驗(yàn)，測(cè)量加根后尾礦砂的內(nèi)聚力(c)與內(nèi)摩擦角(φ)，不同修復(fù)年限的尾礦砂加根量，根據(jù)對(duì)應(yīng)沙棘的生物量確定，加根量分別設(shè)置為 2.0 g(1 a)、3.5 g(4 a)、7.6 g(10 a)，實(shí)驗(yàn)結(jié)果如表1所示。

表1 根土復(fù)合剪切實(shí)驗(yàn)結(jié)果

2 模擬計(jì)算結(jié)果及分析

2.1 分析方法

(1)尾礦庫(kù)模型建立：對(duì)尾礦壩建立幾何模型，主要對(duì)壩體剖面進(jìn)行研究。分析尾礦砂粒級(jí)分布特征，分層特征相近歸并為同一層，進(jìn)而建立計(jì)算模型與地質(zhì)結(jié)構(gòu)分層圖。針對(duì)平面應(yīng)變問(wèn)題進(jìn)行分析，進(jìn)一步將計(jì)算模型轉(zhuǎn)變?yōu)槎S數(shù)值模型，并統(tǒng)一采用三角形單元，使用網(wǎng)格檢驗(yàn)的方法優(yōu)化網(wǎng)格。

(2)修復(fù)尾礦土修復(fù)模型建立：測(cè)量不同修復(fù)年限臺(tái)階植被修復(fù)深度。通過(guò)剪切實(shí)驗(yàn)測(cè)量根-土復(fù)合土的相關(guān)力學(xué)指標(biāo)，包括內(nèi)聚力與內(nèi)摩擦角等；測(cè)量植物生長(zhǎng)狀態(tài)等相關(guān)指標(biāo)，包括植株高度、主根長(zhǎng)等指標(biāo)。據(jù)此建立表層生態(tài)修復(fù)區(qū)的根-土復(fù)合土模型和修復(fù)植物結(jié)構(gòu)模型。復(fù)合土定義為實(shí)體有限單元、主根為錨固單元、上部修復(fù)植物為點(diǎn)荷載、草本植物為分布荷載。

(3)滲流數(shù)值模擬計(jì)算：① 正常水位浸潤(rùn)線：初期壩為堆石透水壩，并選擇適合的滲透參數(shù)，將堆積壩面作為自由透水面，施加水頭邊界后進(jìn)行相關(guān)計(jì)算。對(duì)比、擬合實(shí)測(cè)浸潤(rùn)線與計(jì)算浸潤(rùn)線，進(jìn)而驗(yàn)證、反演滲流參數(shù)。如實(shí)測(cè)浸潤(rùn)線與計(jì)算浸潤(rùn)線擬合不好，需重新選擇參數(shù)。② 洪水位浸潤(rùn)線：施加最小干灘長(zhǎng)度所對(duì)應(yīng)水位的水頭邊界，其他邊界條件不變的條件下進(jìn)行洪水位浸潤(rùn)線計(jì)算。③ 計(jì)算正常水位與洪水位浸潤(rùn)線下的尾礦壩體孔隙水壓力和滲流速度，分析其對(duì)壩體穩(wěn)定性的影響。

(4)靜力分析：模型底部加固定約束，左右兩側(cè)加水平約束，選擇合適參數(shù)，分別對(duì)正常運(yùn)行工況和洪水運(yùn)行工況進(jìn)行靜力穩(wěn)定性計(jì)算。

(5)壩體穩(wěn)定性計(jì)算：采用畢肖普法進(jìn)行計(jì)算，采用《尾礦設(shè)施設(shè)計(jì)規(guī)范》(GB 50863—2013)的規(guī)定的最小安全系數(shù)。

2.2 模型建立

尾礦庫(kù)邊坡穩(wěn)定性分析屬于典型的平面應(yīng)變問(wèn)題，首先確定模型高度，剖面長(zhǎng)度，坡面附加荷載等條件，模擬區(qū)域以(0，0)點(diǎn)為起點(diǎn)，在Y軸方向上確定模型高度的點(diǎn)，剖面長(zhǎng)度在X軸方向上確定，連線建立邊坡模型，見(jiàn)圖2，按彈塑性模型考慮，采用Mohr-Coulomb準(zhǔn)則[15]。

設(shè)置3種試驗(yàn)?zāi)Ｐ徒l件，分別為壩高、坡比、修復(fù)時(shí)間。壩高選取15、30、45、60、80 m；坡比選取1∶3、1∶4、1∶5；修復(fù)時(shí)間選取1 a、4 a、10 a，部分模擬計(jì)算結(jié)果見(jiàn)圖3a～3d(全部計(jì)算結(jié)果見(jiàn)圖4～13)。

2.3 模擬計(jì)算結(jié)果及分析

2.3.1 修復(fù)時(shí)間分析穩(wěn)定性

上述模型分析均采用畢肖普法計(jì)算。模擬計(jì)算結(jié)果見(jiàn)圖4～8。由圖可知，在壩高一定的條件下，安全系數(shù)隨著坡比的降低呈現(xiàn)逐漸升高的趨勢(shì)。坡比為1∶3、1∶4的處理組，其不同修復(fù)年限的安全系數(shù)無(wú)顯著性差異(P>0.05)，坡比為1∶5的處理組則呈現(xiàn)一定差異性(P<0.05)，這說(shuō)明坡比越大，安全系數(shù)越小，修復(fù)植物根系的固土作用越弱。在壩高和坡比一定時(shí)，總體規(guī)律為修復(fù)時(shí)間久的模型安全系數(shù)高：修復(fù)10 a>修復(fù)4 a>修復(fù)1 a，且兩種修復(fù)時(shí)間之間增加幅度較小，差值大部分在0.02～0.04之間波動(dòng)，說(shuō)明植物生態(tài)修復(fù)的年限對(duì)尾礦壩體穩(wěn)定性提高效果一般。在坡比與修復(fù)年限相同的條件下，隨著尾礦壩的增高，安全系數(shù)呈降低趨勢(shì)。無(wú)論坡比與壩高的參數(shù)如何設(shè)定，4 a與1 a相比差異性均不顯著(P>0.05)，個(gè)別還出現(xiàn)了1a的安全系數(shù)大于、等于4a的情況，這可能由于尾礦基質(zhì)營(yíng)養(yǎng)成分差，導(dǎo)致沙棘生長(zhǎng)緩慢，從而降低了安全系數(shù)。所有的模擬計(jì)算結(jié)果中，尾礦壩安全系數(shù)，在壩高為15 m、坡比為1∶5、修復(fù)10 a的條件下，達(dá)到最高值2.84，分別比相同條件下1 a、4 a安全系數(shù)提高了1.45倍、1.42倍，同時(shí)是最小安全系數(shù)的4.8倍，此條件下穩(wěn)定性達(dá)到最大。

2.3.2 坡比分析穩(wěn)定性

模擬計(jì)算結(jié)果見(jiàn)圖9～13。在壩高、修復(fù)年限相同的條件下，隨著坡比的減小其安全系數(shù)呈增大趨勢(shì)，安全系數(shù)差值集中在0.2～0.8之間，安全系數(shù)計(jì)算結(jié)果差異均顯著(P>0.05)，說(shuō)明坡比大小是影響尾礦壩安全系數(shù)的主要因素。依據(jù)畢肖普法計(jì)算，3等尾礦壩最小安全按系數(shù)為1.30，4等和5等尾礦壩最小安全系數(shù)為1.25，由坡比為1∶3的模型可知，當(dāng)壩高為45、60、80 m時(shí)，其安全系數(shù)均小于1.25，故為危庫(kù)或險(xiǎn)庫(kù)；坡比1∶4，壩高為15、30、45 m為安全性能良好的尾礦壩，其余模型為危庫(kù)或險(xiǎn)庫(kù)；坡比1∶5的模型安全系數(shù)均大于規(guī)范要求的最小值，為安全性較高的尾礦壩。上述分析說(shuō)明坡比對(duì)尾礦庫(kù)邊坡穩(wěn)定性影響顯著。

3 結(jié)論

(1)植物修復(fù)尾礦壩模型安全系數(shù)的總體規(guī)律為：修復(fù)10 a>修復(fù)4 a>修復(fù)1 a，差值大部分在0.02～0.04之間波動(dòng)。修復(fù)4 a與修復(fù)1 a效果差異并不顯著。只有當(dāng)坡比為1∶5時(shí)，修復(fù)10 a與修復(fù)1 a、4 a尾礦壩安全系數(shù)才有差異性。表明隨著坡比的增大，可能會(huì)減弱植物根系固土的效果。

(2)在坡比與修復(fù)年限相同的條件下，隨著尾礦壩的增高，其安全系數(shù)呈降低趨勢(shì)。尾礦壩安全系數(shù)，在壩高為15 m、坡比為1∶5、修復(fù)10 a的條件下，達(dá)到最高值2.84，分別比相同條件下1 a、4 a安全系數(shù)提高了1.45倍、1.42倍，同時(shí)是最小安全系數(shù)的4.8倍。

(3)在壩高、修復(fù)年限相同的條件下，隨著坡比的減小其安全系數(shù)呈增大趨勢(shì)，安全系數(shù)差值集中在0.2～0.8之間，安全系數(shù)計(jì)算結(jié)果間差異均顯著(P>0.05)，說(shuō)明坡比大小是影響尾礦壩安全系數(shù)的主要因素。

(4)坡比大小是影響尾礦壩安全系數(shù)的主要因素，尾礦壩修復(fù)植物可以提高安全系數(shù)，但在修復(fù)植物生長(zhǎng)年限短的情況下效果不明顯，且坡比對(duì)植物提高尾礦壩安全系數(shù)效果的影響極大。