化曉鋒

(南水北調(diào)東線山東干線有限責(zé)任公司，山東 濟南 250109)

自由垮落法管理頂板開采方式開采地下煤炭資源必然會引起較大范圍的地面沉降，由此會對沉降范圍內(nèi)地表建(構(gòu))筑物的安全穩(wěn)定性及環(huán)境產(chǎn)生一定的不良影響[1-4]。關(guān)于煤層開采造成的開采沉陷預(yù)測，諸多研究者運用不同手段進行了分析研究工作，如現(xiàn)場監(jiān)測[5-6]、概率積分法[7-8]、數(shù)值模擬[9-10]及相似模型等[11-12]。其中，概率積分法以其計算形式簡單、理論性強和結(jié)果可靠等諸多優(yōu)點目前已成為國內(nèi)最為常用的煤礦開采地表移動量化預(yù)測方法[7-8,13-14]。為此，本文以南水北調(diào)東線二級泵站廠區(qū)平臺南部存在高莊煤礦開采為研究背景，采用概率積分法分析開采沉陷對廠區(qū)平臺的影響并進行預(yù)測，并與現(xiàn)場地表布設(shè)的沉降監(jiān)測結(jié)果進行對比，探究高莊煤礦開采沉陷可能影響泵站安全穩(wěn)定性問題。

1 工程概況

南水北調(diào)東線工程從江蘇揚州附近長江下游干流抽引長江水，進入南四湖下級湖，經(jīng)二級壩泵站提水進入上級湖[15-16]。二級壩泵站是南水北調(diào)東線工程的第10級抽水梯級泵站，工程場地位于蘇魯兩省邊界處的南四湖，是南水北調(diào)東線工程的重要調(diào)蓄水庫和蘇魯兩省分水地點。二級壩泵站為大(1)型泵站，是南水北調(diào)東線工程實現(xiàn)梯級調(diào)水目標的樞紐工程，泵站工程包括主副廠房、變電所、進水閘、引水渠、出水渠、導(dǎo)流渠等一級主體建筑物及二級壩公路橋、引水渠交通橋等三級附屬構(gòu)筑物。泵站廠區(qū)位于南四湖上、下級湖之間的二級壩南側(cè)。

南四湖跨江蘇和山東兩省，湖下煤炭資源豐富，南水北調(diào)東線一期工程共有24座泵站工程，調(diào)查發(fā)現(xiàn)僅二級壩泵站工程受棗礦集團高莊煤礦采煤沉陷影響較大[17-18]。泵站工程場地南側(cè)臨近高莊煤礦開采的西十一采區(qū)，高莊煤礦靠近二級壩泵站廠區(qū)平臺最近的工作面為3上1107工作面。根據(jù)高莊煤礦西十一采區(qū)3上煤回采規(guī)劃，引水渠段大部分處于3上煤規(guī)劃采空范圍，北側(cè)走向回采邊界與泵站廠區(qū)平臺間距最小約350m(見圖1)。根據(jù)2013—2017年廠區(qū)平臺、引水渠、主副廠房等沉降監(jiān)測資料分析，主、副廠房和廠區(qū)平臺受采煤影響較小，但引水渠南部沉降較大，引水渠較大監(jiān)測范圍及南部的交通橋因開采沉陷已沉入水下[19]?？琢钗涞萚20]結(jié)合2013—2019年的變形觀測資料分析了采煤沉陷對二級壩泵站引水渠過流能力影響，認為現(xiàn)狀引水渠過流能力滿足設(shè)計輸水功能要求。

圖1 泵站廠區(qū)平臺與開采工作面相對位置模型

二級壩泵站樞紐工程下面賦存的煤炭資源屬于華北型沉積煤層，為全隱蔽的石炭～二疊紀煤田。對泵站廠區(qū)平臺影響的南部高莊煤礦主采山西組含煤地層，即3煤層，在該區(qū)域大部分又分為3上和3下兩煤層，部分地段合并為一層(3煤)，基本全區(qū)可采。據(jù)初步統(tǒng)計，3上煤層平均厚度5.0m，3下煤層平均厚度4.0m，目前為止3下煤層仍未開采。為與地表監(jiān)測數(shù)據(jù)進行對比，在此主要討論高莊煤礦開采西十一采區(qū)3上煤層后所引起的沉降波及泵站廠區(qū)平臺。

2 基于概率積分法的采動地表變形對廠區(qū)平臺的沉陷分析

根據(jù)泵站廠區(qū)平臺及南部已開采的高莊煤礦西十一采區(qū)工作面位置關(guān)系，結(jié)合該區(qū)域含煤地層結(jié)構(gòu)賦存特征，采用概率積分法對西十一采區(qū)開采沉降及波及范圍進行了計算。首先對西十一采區(qū)進行單元網(wǎng)格劃分，確定相對坐標系，然后代入基于MATLAB基礎(chǔ)上的自行編寫的概率積分法應(yīng)用軟件進行各塊段開采3上煤層所引起的開采沉降的計算，最后進行疊加獲得整個區(qū)域開采3上煤層所引起的西十一采區(qū)沉降及波及范圍。在此基礎(chǔ)上，考慮殘余變形，又計算出3上煤層開采完畢后的殘余沉降量。計算的主要結(jié)果如表1和圖2所示。

表1 泵站廠區(qū)平臺4個永久坐標點沉降量數(shù)據(jù) mm

圖2 3上煤層采后廠區(qū)平臺場地采動沉降變形示意

采用概率積分法計算的開采3上煤層后廠區(qū)泵站平臺4個角點采動沉降量和殘余沉降量及總沉降量對比如表1所示。發(fā)現(xiàn)泵站廠區(qū)平臺從北到南受開采影響沉降逐漸增大。東南部的17號點位置受影響最大，總沉降量達13.6mm。但從數(shù)值的絕對大小來說，該值仍很小，基本接近文獻[13]中規(guī)定的開采后不受影響的穩(wěn)定邊界區(qū)域。

高莊煤礦西十一采區(qū)3上煤層全部開采完畢后所引起的沉降分布特征如圖2a所示。由此可見，最靠近泵站廠區(qū)平臺的3上1107工作面開采完畢后移動邊界線距離泵站廠區(qū)平臺最近的17號點約110m，沉降邊界已越過該點，進入平臺約30m，但這時影響已很小。

根據(jù)泵站廠區(qū)平臺4個角點間距，結(jié)合概率積分法計算預(yù)測的泵站廠區(qū)平臺四角標識點位置的沉降變形累積值(見表1)，結(jié)合4個角點間距分別計算了開采3上煤層后廠區(qū)平臺采動變形梯度(見圖2b)。預(yù)測的泵站廠區(qū)平臺場地在走向、傾向方向的采動變形梯度均在0.01%以下，量化反映了采動引起的場地不均勻沉降程度非常小。

3 理論計算與監(jiān)測對比分析

為保證二級壩泵站安全運行，組織實施了安全監(jiān)測網(wǎng)的布設(shè)與監(jiān)測。南水北調(diào)東線山東干線有限責(zé)任公司于2012年委托山東省水利勘測設(shè)計院進行了主要位置的相關(guān)變形監(jiān)測[21]。為此，筆者對這些監(jiān)測資料進行了系統(tǒng)整理與分析，根據(jù)廠區(qū)泵站平臺主要監(jiān)測點所在建筑物分布特征，組成了集概率積分法計算點和已有歷年監(jiān)測點于一體的泵站廠區(qū)平臺綜合對比圖(見圖3)。

圖3 廠區(qū)平臺剖面線(A—A′，B—B′)及場地主要監(jiān)測點布置平面

為便于分析，建立了廠區(qū)泵站平臺中沿著進水渠和出水渠兩側(cè)如圖3所示A—A′和B—B′兩主要剖面線。根據(jù)圖2所示廠區(qū)平臺場地采動沉降變形坡度，可獲得整體上沿著A—A′和B—B′兩主要剖面線的開采3上煤層沉陷引起的沉降坡降分別為0.024 5，0.024 8mm/m。結(jié)合廠區(qū)平臺主要建筑物進水閘、進水池和主、副廠房沿著A—A′與B—B′兩主要剖面線布置的主要監(jiān)測點各自與A點或B點間距，可初步獲得各監(jiān)測點位置處因開采沉陷所引起的總沉降量(見表2)。由此可獲得廠區(qū)平臺主要建筑物3上煤層開采完畢后的最終總沉降量大致范圍為：進水閘橋附近11.63～12.12mm，進水池附近8.40～8.52mm，主、副廠房附近7.74～8.33mm，出水池附近7.44～7.56mm。對比廠區(qū)平臺沿兩剖面所布設(shè)的地表變形監(jiān)測數(shù)據(jù)來看，12組數(shù)據(jù)基于概率積分法理論計算的總沉降量和實際監(jiān)測的總沉降量基本一致，但各監(jiān)測點對應(yīng)的累計總沉降量均小于開采3上煤層計算的總沉降量。初步分析認為有3方面主要原因：①監(jiān)測數(shù)據(jù)在時間上存在不統(tǒng)一性，監(jiān)測點的監(jiān)測時間有的早有的較晚，導(dǎo)致各監(jiān)測數(shù)據(jù)總沉降量偏小，而理論計算是按預(yù)先設(shè)計開采時間持續(xù)計算的，是統(tǒng)一的；②各監(jiān)測點在埋設(shè)前所在位置可能已發(fā)生開采沉陷，這樣就無法監(jiān)測到前期沉降量，而概率積分法計算是從一開始有沉降就考慮的；③基于概率積分法計算的開采引起總沉降量包括殘余變形，這是需要很長時間才可能達到的結(jié)果。

表2 泵站廠區(qū)平臺主要監(jiān)測點開采3上煤層與監(jiān)測的沉降量部分數(shù)據(jù)對比

主、副廠房主要測點2013年2月至2020年5月沉降量觀測累計沉降量和月變化量歷時變化特征曲線如圖4所示，可明顯看出，6個監(jiān)測點的累計沉降歷時變化曲線基本一致，整體上在波動中基本無大變化，最大沉降量約8.7mm(JH5，2014年2月)(見圖4a)；月沉降量一開始變化較大，最大變化約7.36mm(JH9，2016年1月)(見圖4b)。后期連續(xù)6個月基本不變，根據(jù)文獻[13]相關(guān)規(guī)定，主、副廠房沉降基本穩(wěn)定，開采沉陷對主、副廠房的影響基本停止或沒受開采沉陷影響。

圖4 主、副廠房測點2013年2月至2020年5月沉降量變化特征曲線

限于篇幅，現(xiàn)場對主要建筑物同時進行了水平位移監(jiān)測，如通過對泵站主、副廠房變形監(jiān)測點水平位移累積量和月變化進行分析，發(fā)現(xiàn)這兩個數(shù)據(jù)在各監(jiān)測點對應(yīng)值也均很小，各期坐標差基本處于測量誤差范圍內(nèi)，主、副廠房基本無水平位移產(chǎn)生；在垂直方向上，JH6，JH9，JH10 3點先有小幅度上升，后又緩慢沉降，年度累計沉降量變化不大，各點變化情況總體保持一致，變化幅度均勻且在允許范圍內(nèi)，說明泵站主、副廠房處于穩(wěn)定狀態(tài)。

4 結(jié)語

1)系統(tǒng)分析了二級壩泵站廠區(qū)平臺工程及場地地質(zhì)與煤礦開采條件相互位置關(guān)系，提出了廠區(qū)平臺南部的高莊煤礦開采西十一采區(qū)對泵站安全運行可能存在一定影響，建立了開采影響范圍的理論計算模型。

2)采用概率積分法計算了開采3上煤層對泵站廠區(qū)平臺四角位置的沉降變形累積值和采動變形梯度，量化反映了僅開采3上煤層引起的場地不均勻沉降程度對廠區(qū)泵站平臺各主要設(shè)備的影響。

3)沿著廠區(qū)平臺進水池到出水池方向設(shè)定2條監(jiān)測線(A—A′，B—B′)，結(jié)合監(jiān)測點布置情況，依次計算了12個監(jiān)測點因開采3上煤層所引起的總沉降量，與對應(yīng)的監(jiān)測總沉降量對比分析，發(fā)現(xiàn)兩者沉降量基本一致，監(jiān)測值對應(yīng)小于理論計算值，并結(jié)合主、副廠房監(jiān)測值認為泵站廠區(qū)平臺目前沉降基本穩(wěn)定。