李騰飛　陳洪濤　王瑞青

(中交公路規(guī)劃設(shè)計(jì)院有限公司　北京　100088)

?

湖北宜昌鹽池河滑坡成因機(jī)理分析*

李騰飛陳洪濤王瑞青

(中交公路規(guī)劃設(shè)計(jì)院有限公司北京100088)

1980年6月3日湖北宜昌鹽池河發(fā)生大規(guī)模山體滑坡地質(zhì)災(zāi)害，滑坡體總體積約130×104m3，導(dǎo)致284人死亡。通過(guò)現(xiàn)場(chǎng)調(diào)查和分析后認(rèn)為，地下采礦是鹽池河滑坡的重要原因。針對(duì)已有對(duì)鹽池河滑坡研究中存在不足，考慮到隨機(jī)介質(zhì)理論在巖層移動(dòng)變形分析中的應(yīng)用，采用隨機(jī)介質(zhì)理論方法對(duì)地下采礦作用下鹽池河滑坡體移動(dòng)變形規(guī)律進(jìn)行研究。并結(jié)合邊坡穩(wěn)定性分析和數(shù)值模擬方法，對(duì)采動(dòng)作用下坡體移動(dòng)變形規(guī)律進(jìn)行分析、驗(yàn)證。結(jié)果表明，采礦誘發(fā)的坡體不均勻沉降是鹽池河滑坡的重要原因。采動(dòng)變形導(dǎo)致山體中Ⅰ號(hào)裂縫形成、擴(kuò)展，降低了山體內(nèi)側(cè)對(duì)滑動(dòng)山體的側(cè)向阻滑作用，并加劇地表降雨入滲趨勢(shì)，導(dǎo)致滑帶力學(xué)參數(shù)劣化，最終導(dǎo)致坡體失穩(wěn)。

鹽池河滑坡地下采礦隨機(jī)介質(zhì)理論機(jī)理分析

0　引　言

隨著人類地下采礦活動(dòng)增加，由此引發(fā)的山體滑坡自然災(zāi)害等也時(shí)有發(fā)生，典型案例如發(fā)生在湖北的宜昌鹽池河滑坡(姚寶魁等，1988)。該滑坡災(zāi)害發(fā)生于1980年6月3日，滑坡總體積超過(guò)1.3×106m3，滑坡體從高差約200m的陡崖垮落，形成面積約500m×478m堆積區(qū)并阻塞下方河道，最大堆積厚度近40m。滑坡摧毀該磷礦整個(gè)工礦區(qū)，導(dǎo)致284人死亡和2500萬(wàn)元財(cái)產(chǎn)損失，并誘發(fā)1.6級(jí)地震(黃潤(rùn)秋，2007；連志鵬等，2013)。

鹽池河滑坡位于黃陵背斜東北翼，區(qū)域地質(zhì)構(gòu)造為單斜構(gòu)造，巖層走向總體上呈南北走向，傾向東，巖層傾角近水平約14°，地層巖性為中、厚層塊狀白云巖。受構(gòu)造、地下采礦等作用影響，山體發(fā)育4組主要結(jié)構(gòu)面，幾組陡傾構(gòu)造結(jié)構(gòu)面和巖體中夾含的若干軟弱夾層、巖層面一起，把山體巖體切割呈“積木塊”狀。發(fā)生滑動(dòng)的斜坡上部近于直立，坡向朝北，東、西兩側(cè)也為臨空面，坡體整體上三面臨空(李騰飛，2014)。

鹽池河磷礦區(qū)建于1969年，自1978年冬開始上方陡崖多次發(fā)生小規(guī)模落石和崩塌，1979年7月大規(guī)模爆破房間礦柱后，陡崖處陸續(xù)出現(xiàn)10條地表裂縫，并逐漸發(fā)生變形。1980年6月2～3日，山體內(nèi)側(cè)Ⅰ號(hào)裂縫下沉位移速度達(dá)1m·d-1以上，前緣不斷發(fā)生落石和小型崩滑，山體最終發(fā)生整體性滑坡。磷礦層(開采厚度為2m)賦存在滑坡山體下部，礦體產(chǎn)狀近水平分布，與巖層產(chǎn)狀一致，埋深約300m。截至山體滑坡發(fā)生時(shí)，地下累計(jì)采空區(qū)面積達(dá)6.4×104m2。

圖1　宜昌鹽池河滑坡Fig.1　The Yanchihe landslide in Yichang city

圖2　主要結(jié)構(gòu)面和地層、坡面空間關(guān)系Fig.2　The spatial relationships of joint sets，sratum and slope surface

針對(duì)鹽池河滑坡的失穩(wěn)機(jī)理，姚寶魁等(1988)采用有限元方法進(jìn)行過(guò)分析，主要研究由于地下開采引起的山體內(nèi)的應(yīng)力狀態(tài)變化，結(jié)果表明采礦引起的拉應(yīng)力及其導(dǎo)致山體裂縫產(chǎn)生是坡體失穩(wěn)的主要原因。連志鵬等(2013)以磷礦采礦深厚比(H/M)為評(píng)價(jià)依據(jù)，采用GIS方法對(duì)磷礦采空區(qū)穩(wěn)定性進(jìn)行評(píng)價(jià)。已有研究拓寬了對(duì)鹽池河滑坡失穩(wěn)機(jī)理的認(rèn)識(shí)，但對(duì)磷礦采礦如何影響坡體穩(wěn)定的分析還不夠詳細(xì)。本文以鹽池河滑坡為例，在現(xiàn)場(chǎng)調(diào)查的基礎(chǔ)上，從地下采礦引起上覆山體移動(dòng)變形的角度，通過(guò)隨機(jī)介質(zhì)理論分析和數(shù)值模擬，對(duì)鹽池河滑坡的成因機(jī)理進(jìn)行深入分析，總結(jié)地下采礦條件下山體的移動(dòng)變形規(guī)律，為我國(guó)西南地區(qū)類似滑坡的研究和科學(xué)減災(zāi)提供科學(xué)經(jīng)驗(yàn)和理論依據(jù)。

1　采動(dòng)山體移動(dòng)機(jī)理研究

研究手段采用隨機(jī)介質(zhì)理論方法，該方法由波蘭學(xué)者李特維尼申(Litwinszyn，1957)建立，目的是計(jì)算地表移動(dòng)，我國(guó)學(xué)者劉寶琛等(劉寶琛等，1965)將其應(yīng)用于“三下”開采評(píng)價(jià)中。該方法主要優(yōu)勢(shì)在于構(gòu)建起地下采礦與上覆巖體變形之間的聯(lián)系。認(rèn)為地下采礦引起的上覆巖體移動(dòng)變形，由單元體開采引起的巖體內(nèi)部單元下沉曲線(何國(guó)清等，1994；胡斌等，2012)累積組成。地下開采下沉和水平移動(dòng)變形的公式如下：

(1)

(2)

(3)

公式中，m為礦層開采厚度(m)；k為下沉系數(shù)；s為積分變量，為開采坐標(biāo)；Ha為最小開采深度(m)；α為計(jì)算剖面的礦層傾角(°)；θ為地下開采影響傳播角(°)。

圖3　鹽池河滑坡堆積區(qū)Fig.3　The Yanchihe landslide area1.地下采空區(qū)界限；2.分段號(hào)；3.采礦及放頂日期；4.裂縫編號(hào)及其生產(chǎn)的日期；5.崩坍山體界限；6.崩坍物堆積范圍

對(duì)于鹽池河滑坡體，隨機(jī)介質(zhì)理論方法的適用性體現(xiàn)為：基于山體千米尺寸，考慮兩組近于正交的陡傾結(jié)構(gòu)面和近水平層面的影響(圖2)。塊狀巖體結(jié)構(gòu)在一定程度上可認(rèn)為是均質(zhì)的松散介質(zhì)，這與該方法的假設(shè)條件是相吻合的。

為簡(jiǎn)化分析，兩個(gè)剖面滑動(dòng)方向AA′和崩塌剖面BB′地下采礦都按照4個(gè)開采步長(zhǎng)(step)，礦柱寬度設(shè)置為8m(根據(jù)現(xiàn)場(chǎng)調(diào)查)，斜坡下磷礦層厚度平均約2.5m。

表1　計(jì)算參數(shù)

Table 1　The calculation parameters

分析剖面Section采礦步長(zhǎng)step/m礦層厚度m/m最小開采深度Ha/m巖層視傾角α/(°)AA'BB'87.3955.52.52.529.4965912.830剖面Section水平移動(dòng)系數(shù)B(z)地表下沉系數(shù)k主要影響角β(z)/(°)影響傳播角θ/(°)AA'BB'—0.40.650.6557578090

1.1滑動(dòng)方向移動(dòng)分析

圖4為山體滑動(dòng)方向剖面圖，分4布開采。數(shù)據(jù)記錄點(diǎn)1-5沉降變形(圖5)。分析可知，數(shù)據(jù)記錄點(diǎn)垂直沉降最大值出現(xiàn)在最靠近采空區(qū)的5點(diǎn)，最大沉降值到達(dá)0.9m。特別是后兩步開采完成后沉降值急劇增加，這主要是因?yàn)楹髢刹介_采位于變形坡體正下方。采動(dòng)變形導(dǎo)致AA′剖面的巖層面傾角降低，減緩了坡體的下滑趨勢(shì)。

圖4　山體滑動(dòng)方向剖面Fig.4　The section of landslide direction

圖5　山體滑動(dòng)方向數(shù)據(jù)點(diǎn)沉降變形Fig.5　The settlements of data point landslide direction section

圖6　山體崩塌方向剖面Fig.6　The section of collapse direction

圖7　山體崩塌方向數(shù)據(jù)點(diǎn)沉降變形Fig.7　The settlements of data points in collapse direction section

圖8　山體崩塌方向數(shù)據(jù)點(diǎn)水平移動(dòng)Fig.8　The horizontal deformation of data points in collapse direction section

1.2崩塌方向移動(dòng)分析

圖6為山體崩塌方向剖面圖。由圖7分析可知，地下開采作用導(dǎo)致的外側(cè)變形沉降趨勢(shì)更明顯。坡體內(nèi)外側(cè)采動(dòng)沉降差異，導(dǎo)致鹽池河山體表現(xiàn)出向外側(cè)崩塌變形的趨勢(shì)。山體內(nèi)側(cè)分布Ⅰ號(hào)裂縫，其所在位置數(shù)據(jù)點(diǎn)1-5水平移動(dòng)情況(圖8)，坡頂?shù)乃阶冃沃底畲?，這說(shuō)明坡體水平移動(dòng)導(dǎo)致山體產(chǎn)生張拉作用，這與實(shí)際山體中內(nèi)側(cè)Ⅰ號(hào)裂縫相吻合。

綜上所述，地下開采誘發(fā)的山體變形導(dǎo)致山體內(nèi)側(cè)Ⅰ號(hào)裂縫形成，使得潛在滑坡體實(shí)際上被切割呈四面臨空的狀態(tài)。

對(duì)于鹽池河滑坡體穩(wěn)定性的影響，地下采礦主要體現(xiàn)在以下幾個(gè)方面。作用①：走向NW66°剖面即AA′剖面采動(dòng)變形使得巖層面傾角降低，最大沉降差值達(dá)0.54m，導(dǎo)致滑動(dòng)面角度降低0.28°，這反而對(duì)坡體穩(wěn)定有利。作用②：走向NS 剖面即BB′剖面采動(dòng)變形使得山體南側(cè)Ⅰ號(hào)裂縫產(chǎn)生，導(dǎo)致滑坡體處于實(shí)際上的“四面臨空”狀態(tài)，減弱了山體內(nèi)側(cè)巖體對(duì)滑動(dòng)體的側(cè)向阻擋，對(duì)坡體穩(wěn)定不利。作用③：內(nèi)側(cè)Ⅰ號(hào)張拉裂縫產(chǎn)生，客觀上加劇了降雨入滲作用。而降雨入滲導(dǎo)致滑帶軟化則直接降低了坡體的穩(wěn)定性。

2　邊坡穩(wěn)定性分析

根據(jù)上述分析，斜坡體存在東側(cè)、北側(cè)和西側(cè)三面臨空的地形條件，內(nèi)側(cè)即南側(cè)Ⅰ號(hào)裂縫產(chǎn)生，組合作用下斜坡體實(shí)際上表現(xiàn)為四面臨空，滑坡體此時(shí)與斜面體系類似。影響斜面體系穩(wěn)定有以下兩個(gè)原因：斜面角度和斜面力學(xué)參數(shù)，這兩個(gè)原因恰與上述地下采礦影響作用①和作用③向?qū)?yīng)。將斜坡體簡(jiǎn)化為斜面體系，作用①和作用③所引起采動(dòng)山體安全系數(shù)變化程度即為其影響。

表2　采動(dòng)影響分析

Table 2　The analysis of undermining influence

影響作用作用①作用②作用③表現(xiàn)形式降低滑動(dòng)面傾角Ⅰ號(hào)裂縫產(chǎn)生→降低山體內(nèi)側(cè)阻滑力Ⅰ號(hào)裂縫產(chǎn)生→降雨入滲滑帶劣化對(duì)坡體穩(wěn)定影響有利不利不利

對(duì)于鹽池河滑坡體，地下采礦作用誘發(fā)的不均勻?qū)е翧A′剖面不均勻沉降，最大沉降值達(dá)0.54m，從而使視傾角由11.94°變?yōu)?1.70°，進(jìn)而導(dǎo)致對(duì)應(yīng)的真傾角(滑動(dòng)方向角)由14°變?yōu)?3.72°，坡體安全系數(shù)提高1.9%，即地下采礦影響作用①對(duì)山體穩(wěn)定的有利影響。

地下采礦影響作用③具體體現(xiàn)為滑帶遇水前后力學(xué)性質(zhì)參數(shù)的變化。內(nèi)側(cè)Ⅰ號(hào)裂縫擴(kuò)張導(dǎo)致的滑帶遇水飽和后的滑坡體安全系數(shù)降低17%，即地下采礦影響作用③對(duì)山體穩(wěn)定的不利影響。

3　數(shù)值模擬分析

對(duì)鹽池河采動(dòng)滑坡體進(jìn)行數(shù)值模擬分析，采用基于連續(xù)介質(zhì)離散元分析，該方法理論依據(jù)是動(dòng)態(tài)松弛法求解。結(jié)合現(xiàn)場(chǎng)調(diào)查，確定模型大小和采空區(qū)位置，共39.4萬(wàn)單元。根據(jù)實(shí)際沿中央斜井自西向東開采情況，模型中斜井的兩側(cè)位置設(shè)置5m寬的安全礦柱，礦房也按每隔4m的間距設(shè)置3m×5m礦柱。根據(jù)實(shí)際調(diào)查，模型中地下采礦層厚度為3m，全部21個(gè)采空區(qū)簡(jiǎn)化為3步完成開采，按照第1步A1(1～7)→第2步A2(8～16)→第3步A3(17～21)的順序累計(jì)開采。在模型中設(shè)置了若干數(shù)據(jù)記錄點(diǎn)，以便于分析地下采礦作用對(duì)山體模型的影響。

圖9　鹽池河斜坡地質(zhì)模型Fig.9　The geological model of Yanchihe landslide

圖10　模型中開采順序Fig.10　The mining sequence in the geological model

圖11　開采完成后坡體失穩(wěn)情況(放大5倍)Fig.11　The undermining-induced landslide(×5)

圖12　滑帶沉降變形Fig.12　The settlements of slip zone

鹽池河滑坡體表現(xiàn)為整體滑動(dòng)特征，失穩(wěn)前坡體前緣不斷發(fā)生滾石和小型崩塌。圖11 為地下采礦后坡體的移動(dòng)變形情況，分析可知隨地下采礦的不斷推進(jìn)，原山體內(nèi)側(cè)裂縫逐漸擴(kuò)張，陡崖外側(cè)變形破壞增加，形成整體破壞。

采礦作用下滑坡體底面軟弱層帶沉降變形(圖12)，隨開采面積增加，沉降變形不斷增加，靠近陡崖外側(cè)沉降量最大，最外側(cè)沉降接近1m。圖13為地下采礦作用下滑坡體內(nèi)側(cè)裂縫水平移動(dòng)變形，靠近山頂處水平變形最大。

圖13　內(nèi)側(cè)裂縫水平移動(dòng)變形Fig.13　The horizontal deformation of inside crack

4　結(jié)　論

針對(duì)鹽池河滑坡體的失穩(wěn)機(jī)制，分別采用隨機(jī)介質(zhì)理論方法和數(shù)值模擬進(jìn)行分析。雖然隨機(jī)介質(zhì)方法和數(shù)值計(jì)算得出的采動(dòng)變形量有所不同，但兩種方法得出的坡體移動(dòng)變形規(guī)律是相同的。地下開采磷礦誘發(fā)上覆巖層發(fā)生不均勻沉降，導(dǎo)致山體有向外側(cè)發(fā)生傾倒的趨勢(shì)，這也與原山體中Ⅰ號(hào)裂縫形成、擴(kuò)展相吻合。采動(dòng)作用下坡體內(nèi)側(cè)Ⅰ號(hào)裂縫逐漸擴(kuò)張的直接結(jié)果是減少了山體內(nèi)側(cè)對(duì)滑動(dòng)山體的側(cè)向阻滑作用，會(huì)對(duì)坡體穩(wěn)定性產(chǎn)生不利影響。坡體內(nèi)側(cè)Ⅰ號(hào)裂縫逐漸擴(kuò)張的間接結(jié)果則是加劇了地表降雨入滲趨勢(shì)，從而導(dǎo)致滑帶力學(xué)參數(shù)的劣化，進(jìn)一步降低坡體的安全系數(shù)。

He G Q，Yang L.1994.Mine subsidence induced by undermining[M].Xuzhou，China University of Mining and Technology Press：118～147．

Hu B，Liu Y L，Tang H M，et al.2012.Research on ground subsidence due to tunnel excavation in Huquan-Mingdu section of Wuhan subway[J].Chinese Journal of Rock Mechanics and Engineering，31(5)：908～913.

Huang R Q.2007.Large-scale landslides and their sliding mechanism in China since the 20th century[J].Chinese Journal of Rock Mechanics and Engineering，26(3)：433～454．

Li T F.2014.Research on the Mechanism of Mountainous Landslide Geohazrds Induced by Underground Mining[Ph.D.Thesis][D].Beijing：University of Chinese Academy of Sciences．

Lian Z P，Tan J M，Li J F.2013.Stability evaluation of Yanchihe phosphate mine in Yuan’an，Yichang，China[J].Geology and Mineral Resources of South China，29(1)：60～65．

Litwinszyn J.1957.Fundamental principles of the mechanics of stochastic medium[C]∥Proceedings of the 3rd Conference Theo.Appl.Mech.Bongalore，India：[s.n.]：50～61．

Liu B C，Liao G H.1965.Basic law of surface movement in coal mine[M].Beijing：China Industry Press：50～78.

Yao Ba K，Sun Y K.1988.Yanchi River phosphate landslides and its collapse failure mechanism in Yichang[C].∥Ground rock Engineering Committee，CSRME & Engineering geology Committee，GSC.Typical landslides in China，Beijing：Science Press：89～98．

何國(guó)清，楊倫.1994.礦山開采沉陷學(xué)[M].徐州：中國(guó)礦業(yè)大學(xué)出版社：118～147．

胡斌，劉永林，唐輝明，等.2012.武漢地鐵虎泉-名都區(qū)間隧道開挖引起的地表沉降研究[J].巖石力學(xué)與工程學(xué)報(bào)，31(5)：908～913.

黃潤(rùn)秋.2007.20世紀(jì)以來(lái)中國(guó)的大型滑坡及其發(fā)生機(jī)制[J].巖石力學(xué)與工程學(xué)報(bào)，26(3)：433～454．

李騰飛.2014.地下采礦誘發(fā)山體崩滑地質(zhì)災(zāi)害機(jī)理研究[D].北京：中國(guó)科學(xué)院大學(xué)．

連志鵬，譚建民，李景富.2013.湖北遠(yuǎn)安鹽池河磷礦開采區(qū)穩(wěn)定性評(píng)價(jià)[J].華南地質(zhì)與礦產(chǎn)，29(1)：60～65．

劉寶琛，廖國(guó)華.1965.煤礦地表移動(dòng)的基本規(guī)律[M].北京：中國(guó)工業(yè)出版社：50～78．

姚寶魁，孫玉科.1988.宜昌鹽池河磷礦山崩及其崩塌破壞機(jī)制[C].∥中國(guó)巖石力學(xué)與工程學(xué)會(huì)地面巖石工程專業(yè)委員會(huì)&中國(guó)地質(zhì)學(xué)會(huì)工程地質(zhì)專業(yè)委員會(huì).中國(guó)典型滑坡，北京：科學(xué)出版社：89～98．

FORMATION MECHANISM OF YANCHIHE LANDSLIDE IN YICHANG CITY，HUBEI PROVINCE

LI TengfeiCHEN HongtaoWANG Ruiqing

(CCCC Highway Consultants CO.，Ltd.，Beijing100088)

A large scale of landslide disaster occurred on June 3，1980 at Yanchihe in Yichang city of Hubei province.Its volume is 130×104m3.The landslide killed 284 people.Based on the site investigation and analysis，underground mining is an important reason of Yanchihe landslide.Considering the limitations of previous research，and the paper uses the stochastic medium theory that has been used widely in ground surface movement and deformation in diggings.The slope movement and deformation principle resulted from underground mining are illuminated with the stochastic medium theory.Combined with slope stability analysis and numerical simulation，deformation of slope movement induced by underground mining is analysed and verified.The conclusions are as follow：The underground mining caused differential settlements.The underground mining deformation exacerbated the formation of crack inside the slope，providing a channel for rainfall infiltration.Rainfall infiltration in crack caused the slip zone softening，causing the landslide．

Yanchihe landslide，Underground mining，Stochastic medium theory，F(xiàn)ormation analysis

10.13544/j.cnki.jeg.2016.04.012

2015-07-11;

2015-12-31.

李騰飛(1985-)，男，博士，工程師.主要從事高速公路工程地質(zhì)和路基路面方面的研究.Email: Robustfei@126.comcn

P642.22

A