人工湖下連采連充煤炭開采地表移動變形預測研究

張國建 孟浩 熊威 白濤 孟現(xiàn)臣 王軍 呂曉

摘要：為研究人工湖下煤炭連采連充開采后的地表移動變形規(guī)律，對連采連充工作面進行了室內力學試驗和現(xiàn)場取芯力學試驗，驗證充填體的可行性；基于等價采高概率積分法，對連采連充工作面進行地表沉陷預測；通過數(shù)值模擬計算，分析導水裂縫帶發(fā)育高度，并與概率積分法結果進行對比。結果表明：充填體的強度5.063 MPa超過設計強度2.0 MPa，能夠確保礦井的安全開采；連采連充工作面開采后，地表傾斜值極值0.3 mm/m，地表水平變形極值-0.2 mm/m，小于磚混結構建筑I級損壞范圍，周圍地表沉陷平緩，無安全隱患；導水裂縫帶發(fā)育高度約49.7 m，距離隔水層約160.3 m，水下采煤安全，F(xiàn)LAC3D數(shù)值模擬和概率積分法結果比較接近，驗證連采連充技術可有效減緩地表移動變形。

關鍵詞：連采連充；充填體；地表移動變形預測；概率積分法；數(shù)值模擬

中圖分類號：TD823?? 文獻標志碼：A?? 文章編號：1002-4026（2023）05-0033-11

Predicting surface movement and deformation for continuous mining and continuous backfilling under an artificial lake

ZHANG Guojian1，2a，3，MENG Hao1*，XIONG Wei4，BAI Tao5，MENG Xianchen5，WANG Jun2b，LV Xiao2b

（1.School of Environment and Spatial Informatics， China University of Mining and Technology， Xuzhou 221116， China;2. a.School of Surveying and Geo-Informatics; b. School of Civil Engineering， Shandong Jianzhu University，Jinan 250101， China; 3. Technology Research Institute of Shandong Energy Group， Jinan 250101， China;4. Shandong Road and Bridge Engineering Design Consulting Institute， Jinan 250000， China;5. Inner Mongolia Yuxing Mining Institute， Chifeng 016064， China）

Abstract∶To investigate surface movement and deformation characteristics due to continuous mining and continuous backfilling （CMCB）of coal under artificial lakes， laboratory and field coring mechanical tests were conducted on the CMCB area to verify the feasibility of the filling body. Based on the equivalent mining height probability integration method， the surface subsidence of the CMCB area was predicted. The height of the water-conducting fracture zone was analyzed using numerical simulation， and its results were compared with those of the probability integration method. The results show that the strength of the filling body is 5.063 MPa， which is higher than the designed strength of 2.0 MPa， ensuring safe mining.Owing to continuous mining and backfilling in the area， the maximum inclination value of the surface was 0.3 mm/m and the maximum horizontal deformation value of the surface was -0.2 mm/m， respectively， which is less than the range of grade Ⅰ damage to brick and concrete structures. The surrounding surface subsidence was gentle， and there was no safety hazard. The height of the water-conducting fracture zone was about 49.7 m，and the distance from the waterproof layer was about 160.3 m， indicating the safety of underwater coal mining. Results of the FLAC3D numerical simulation and probability integration method were close， thereby verifying that the CMCB technology can effectively slow down surface movement and deformation.

Key words∶continuous mining and continuous backfilling;backfill; surface movement and deformation prediction; probability integration method; numerical simulation

在河流、湖泊、池塘、水庫等水體下方進行的采煤工作叫做水體下采煤。在我國，水體下采煤曾引起多次突水安全事故，造成嚴重的經(jīng)濟損失和人員傷亡。充填開采是進行水體下安全采煤的有效防治措施，是煤礦綠色開采的研究方向之一[1-3]。當采空區(qū)采用充填法管理頂板時，可有效抑制采煤引起的導水裂縫帶發(fā)育，減輕工作面支撐壓力產(chǎn)生的礦壓顯現(xiàn)，減緩地表沉陷，實現(xiàn)水體下安全采煤。因此，研究水體下充填采煤技術具有重要的理論價值和研究意義。

近些年來，許多學者對采煤引起地表沉陷的控制方法進行研究。在充填開采方面，王炯等[4]根據(jù)實測數(shù)據(jù)資料對花園煤礦充填開采地表沉陷預計參數(shù)進行計算反演，研究出充填開采巖層移動變形特征。郭廣禮等[5]通過相似材料模擬和數(shù)值模擬實驗探究了充填開采的地表沉陷規(guī)律，發(fā)現(xiàn)充填開采巖層下沉緩慢，上覆巖層破壞程度較低，僅在采空區(qū)附近發(fā)育一定范圍的裂縫帶，且發(fā)育高度有限。在條帶開采方面，王冰等[6]采用數(shù)值模擬的方法，提出煤層錯動系數(shù)用于解釋條帶開采的地表移動變形、上覆巖層的應力變化規(guī)律。殷和健等[7]采用理論分析和數(shù)值模擬，設計出采留比的優(yōu)化設計方案，研究在特定地質采礦條件下的條帶開采最佳采寬和留寬。但是利用充填開采管理頂板時，需要專門的充填材料及充填系統(tǒng)，這會額外增加開采成本；而條帶開采往往會造成煤炭資源的浪費，煤炭采出率低，工作效率低。針對上述問題，國內有學者提出了全負壓式連采連充采煤技術[8]。經(jīng)過多年來的理論研究和現(xiàn)場工程應用，不斷改進充填采煤工藝、研制相應的采充系統(tǒng)，形成了比較完善的煤礦連采連充式充填采煤技術。

本研究驗證了赤峰市某礦連采連充技術充填體的可行性。采用等價采高概率積分法預測地表移動變形。同時，使用FLAC3D軟件模擬了煤礦連采連充數(shù)值模型，分析了導水裂縫帶高度，并將其與概率積分法結果進行對比。通過綜合實驗和數(shù)值模擬方法，研究了人工湖下煤炭連采連充后地表移動變形，提供了地表沉陷預測，驗證了技術的安全性和有效性。

1 煤礦概況

1.1 采區(qū)概況

該煤礦位于內蒙古自治區(qū)鄂托克旗境內，行政區(qū)劃歸赤峰市鄂托克旗棋盤井鎮(zhèn)鎮(zhèn)政府管轄。礦區(qū)屬溫帶半荒漠大陸性氣候區(qū)，日照豐富，太陽輻射強烈，年總降水量小于200 mm，蒸發(fā)量大。礦區(qū)地勢北高南低。圖1為礦區(qū)綜合水文柱狀圖，下二疊系中下統(tǒng)和上石炭中下統(tǒng)為本區(qū)主要含煤地層，含2～16號共13個煤層。主要由灰色中粒砂巖，雜砂泥巖，細砂巖，灰色中粒、細粒砂巖，灰黑色砂泥巖，灰色粗粒、細粒石英砂巖，泥巖和煤層組成。

1.2 連采連充工作面

連采連充工作面的開采厚度4.5～6.6 m，平均6 m，煤層傾角在6°～14°，平均8°。井下標高+1 134 m，地面標高+1 314 m，位于棋盤井鎮(zhèn)鎮(zhèn)政府以西500 m處，工作面周圍地形多為丘陵，地面無大型建筑。工作面由相鄰布置的若干采場支巷和保護煤柱組成，間隔5 m分別布置采場支巷和保護煤柱，形成奇數(shù)巷和偶數(shù)巷的相鄰對稱的布置局面。連采連充采場可選擇長壁式工作面形式布設回風巷、運輸巷及開切眼，使工作面整體上形成負壓通風模式。在實際生產(chǎn)開采作業(yè)過程中，采用大型綜采設備進行開采，無軌化輪胎煤車牽引運輸煤炭。開采之前，將采場支巷劃分成2～5個開采階段，為確保充填和開采始終處于穩(wěn)定且獨立的作業(yè)空間內，采場支巷在同一開采階段需要彼此獨立，并保持一定距離的間隔，采場支巷在同一開采階段需采用連續(xù)采煤設備依次進行開采，并對采煤結束的采場支巷立即密閉并迅速充填，使得工作面呈現(xiàn)出連采連充的作業(yè)模式，直至所有采場支巷開采并充填完成[9-11]。

2 充填體物理力學性能

2.1 室內力學試驗

矸石、水泥、粉煤灰、水的質量分數(shù)分別為：77.5%、7.0%、5.0%、10.5%，選取立方體28 d試塊取標準巖芯后在實驗室做了圍壓為1、3 MPa三軸壓縮試驗，試驗結果見圖2，試驗強度極限、殘余強度結果如表1所示。

2.2 現(xiàn)場取芯力學試驗

本次試驗以該礦的連采連充工作面作為工程實踐背景，旨在研究充填體的力學性能。因此在現(xiàn)場分別鉆取了10、30、60 m深度的鉆孔已獲取充填體試樣。在取芯后，立即進行了密封保存，以最大程度地防止風化并保持原貌。由于本節(jié)主要考察充填體的強度，因此主要完成了與單軸抗壓強度相關的試驗數(shù)據(jù)，即從現(xiàn)場取得的部分充填體試塊，加工成標準充填體試樣進行單軸壓縮試驗，圖3所示為充填體試塊及標準試件，充填體單軸壓縮試驗結果見表2。

根據(jù)《煤礦采空區(qū)建（構）筑物地基處理技術規(guī)范》（GB 51180—2016）[12]中要求，綜合考慮煤礦地質采礦條件，要求充填體的最終單軸抗壓強度不低于2.0 MPa?，F(xiàn)場取芯力學試驗中，對連采連充工作面5# 10 m、9# 10 m、5# 30 m、9# 30 m、5# 60 m、9# 60 m進行了現(xiàn)場取芯，并進行了單軸壓縮試驗，結果表明：現(xiàn)場取樣充填體單軸抗壓強度平均為5.093 MPa。通過不同位置、不同區(qū)段的充填體取芯試驗表明，煤矸石通過連采連充工藝充填后，充填體的強度超過設計強度2.0 MPa，能夠確保礦井的安全開采。

3 連采連充地表移動變形預測

3.1 預測方法

本文采用等價采高概率積分法進行變形預測。在采用傳統(tǒng)的完全垮落法處理采空區(qū)時，上覆巖層下沉和地表移動變形程度主要是由煤層實際采高決定。在相同地質采礦條件下，對礦區(qū)進行充填開采，充填采空區(qū)的充填體會占據(jù)上覆巖層一部分有效下沉空間，從而控制上覆巖層的下沉和移動變形，使得上覆巖層下沉和移動變形的主控因素變成了由實際采高減去充填體高度后的剩余高度，這被稱為等價采高[13]，如圖4所示。

假設巖層和巖體介質是由眾多微小的顆粒單元組成的介質，這些顆粒單元之間沒有聯(lián)系，保持著獨立性。大量單元顆粒介質之間的相互移動被視為一個隨機過程。在水平煤層條件下，單元開采導致單元下沉盆地的剖面形態(tài)呈正態(tài)分布概率密度函數(shù)。這一方法被稱為概率積分法。通過查詢相關資料，我們發(fā)現(xiàn)所研究煤礦開采引起的沉陷與概率積分法預測模型是相符的。

利用基于等價采高的概率積分法，可以預測充填開采引起的地表移動變形，預測公式如下：

式中：w（x）、i（x）、k（x）、u（x）、ε（x）分別表示為地表下沉、地表傾斜、地表曲率、地表水平移動、地表水平變形，r表示為主要影響半徑，W0表示為最大下沉值，b表示為水平移動系數(shù)。

根據(jù)等價采高原理，可以利用等價采高代替實際采高，然后應用概率積分法進行地表沉陷預測。這種方法被稱為基于等價采高的概率積分法，其預測參數(shù)可以通過對實測地表移動變形數(shù)據(jù)進行擬合計算來確定。

3.2 預測參數(shù)選取

礦井在二采區(qū)東翼9#煤層充填開采過程中均設置了地面移動觀測站，故本次預測參數(shù)的選取參考二采區(qū)東翼9#煤層充填開采時的觀測數(shù)據(jù)。通過對觀測的下沉曲線進行擬合計算，其擬合結果見擬合曲線函數(shù)圖5，相關系數(shù)為0.965 25，滿足精度要求，得到充填開采概率積分法預計參數(shù)：下沉系數(shù)q=0.5、水平移動系數(shù)b=0.2、主要影響角正切tan β=1.8、開采影響傳播角θ=87°、拐點偏移距S=0 m。

3.3 地表移動變形預測

基于等價采高概率積分法，采用表4充填開采預測參數(shù)，使用山東科技大學特殊開采研究所設計的巖移數(shù)據(jù)處理系統(tǒng)進行計算，得到連采連充地表移動變形數(shù)據(jù)，利用Surfer軟件繪制地表移動變形等值線圖。連采連充工作面充填開采后地表下沉及水平移動等值線圖如圖6所示。

連采連充工作面充填開采后，地表最大下沉值34.0 mm；走向方向地表傾斜極值0.2 mm/m，傾向方向地表傾斜極值0.3 mm/m；走向方向地表水平移動極值7.0 mm，傾向方向地表水平移動極值11.0 mm；走向方向地表水平變形極值-0.1 mm/m，傾向方向地表水平變形極值-0.2 mm/m。遠小于《建筑物、水體、鐵路及主要井巷煤柱留設與壓煤開采指南》[14]中的磚混結構建筑Ⅰ級損壞范圍（水平變形<2.0 mm/m；i<3.0 mm/m）。連采連充工作面充填開采后，觀光塔區(qū)域地表沉陷平緩，無安全隱患。

4 數(shù)值模擬

4.1 數(shù)值模型建立

根據(jù)連采連充工作面的地質采礦條件，對巖層的厚度進行近似取整，傾角簡化為0°，簡化后的模型沿地表方向垂直從上往下依次為：表土層、中粒砂巖、雜砂質泥巖、中粗粒砂巖、砂泥巖、細粒砂巖、煤層、砂質泥巖[15]。其模型巖層物理力學參數(shù)見表3。

為研究連采連充工作面充填開采條件下的地表下沉和地表移動變形規(guī)律，采用莫爾－庫倫模型，建立數(shù)值模擬模型，模型的尺寸為1 200 m×1 200 m×218 m（長×寬×高），模型共劃分為187 200個單元、307 461個節(jié)點，模型從左到右開采，為減小邊界效應影響，左側留設200 m，一次開采170 m，共充填開采4次，推進長度800 m，模型底部邊緣沿y軸方向位移是恒定，兩側邊界邊緣沿x軸方向位移是恒定，而模型頂部不設置位移約束條件，為自由表面[16]。建立的數(shù)值模型見圖7所示。[]

4.2 數(shù)值模擬結果分析

經(jīng)過圖8和圖9的數(shù)值模擬結果以及塑性區(qū)導水裂縫帶分析，隨著充填開采次數(shù)的增加，上覆巖層的塑性區(qū)破壞程度逐漸加劇，導水裂縫帶的發(fā)育高度與之呈正相關。在工作面開采完畢后，導水裂縫帶的發(fā)育高度約為49.7 m。此后高度基本穩(wěn)定，由此得出FLAC3D模擬的導水裂縫帶最大發(fā)育高度為49.7 m。

選取模擬數(shù)值49.7 m作為該連采連充工作面的導水裂縫帶發(fā)育高度，人工湖處的煤層埋藏深度在204～210 m左右，最小的防水安全煤柱的高度為154.3 m，符合防水安全煤巖柱的高度要求，導水裂縫帶不會波及人工湖。數(shù)值模擬塑性區(qū)顯示，煤層頂板以彎曲下沉為主，不存在垮落帶，上覆巖層僅產(chǎn)生少量斷裂隙，主要集中在煤層上方0～3 m，故人工湖不會對工作面開采產(chǎn)生影響，湖水不會流入礦井下，可實現(xiàn)水下安全開采[17]。

利用FLAC3D軟件建立數(shù)值模型，提取在連采連充工作面充填開采后的地表下沉值和地表水平移動值。圖10是利用Origin軟件繪制出的連采連充工作面開采后的地表下沉和地表水平移動曲線圖。本研究比較了兩種不同的方法：一種是基于FLAC3D軟件的數(shù)值模擬，另一種是基于等價采高概率積分法的分析。對這兩種方法的相對誤差進行對比，旨在評估它們在解決采煤引起的地表沉陷問題時的一致性和準確性。具體結果見表4。

由表4可知，用FLAC3D軟件模擬的連采連充工作面開采后產(chǎn)生的地表最大下沉值是28.742 2 mm，地表最大水平移動是9.701 6 mm，而利用等價采高概率積分法求的地表最大下沉值是34 mm，地表最大水平移動是11 mm。兩種方法地表最大下沉值相差5.257 8 mm，相對誤差約為16.8%，地表最大水平移動相差1.298 4 mm，相對誤差約12.5%；[JP]利用概率積分法求的地表下沉和水平移動均大于數(shù)值模擬結果，其產(chǎn)生誤差的原因可能是利用FLAC3D進行數(shù)值模擬巖層時，其復雜巖體的裂縫、節(jié)理等因素未充分考慮、未能著重考慮時間效應等非采動因素造成的沉陷，造成數(shù)值模擬結果偏小。假如考慮上述因素，相對誤差將進一步減小，概率積分法和數(shù)值模擬的結果也可以進一步驗證出連采連充技術可以有效減緩地表移動變形。

4.3 地表沉陷機理解釋

如圖11（a）所示，常規(guī)垮落法對頂板進行管理，在工作面推進過程中，由于上覆巖層壓力增加和巖層自身重力影響，頂板會因為開采空間的增大，產(chǎn)生向下的彎曲變形。當頂板上覆巖層極限荷載壓力小于其上覆巖層壓力時，會使頂板產(chǎn)生破碎和斷裂現(xiàn)象，形成許多的碎塊巖石，雜亂地堆積在采空區(qū)內，形成垮落帶。由于開采空間的進一步增大，更會出現(xiàn)離層和斷裂問題，形成斷裂帶。水下采煤時，垮落帶和斷裂帶合稱導水裂縫帶。這個現(xiàn)象會在開采過程中循環(huán)發(fā)生，直至上覆巖層和垮落的破碎巖石重新達到應力平衡后結束，在斷裂帶之上形成彎曲帶[18]。

如圖11（b）所示，工作面采用連采連充技術時，充填體主要由開采出的煤矸石和其他廢料組成。采空區(qū)的充填體占據(jù)了上覆巖層的一部分有效下沉空間。上覆巖層在下沉過程中很快接觸到充填體，在充填體的支撐作用之下，整個上覆巖層的下沉速度將逐漸變得緩慢，在充填體逐漸被壓實后，其承載應力迅速增加、壓縮變形逐漸減小并最后趨向于穩(wěn)定。在充填體支撐作用限制下，上覆巖層垮落空間極小或無垮落空間，僅發(fā)育斷裂帶和彎曲帶，導水裂縫帶發(fā)育高度低，有利于水下安全采煤[19]。因此在連采連充工作面開采時，上覆巖層沉陷緩慢，這就是連采連充技術可以有效減緩地表移動變形的原因。

5 結論

本文以赤峰市某礦連采連充技術為例，通過充填體力學實驗，采用等價采高概率積分法對連采連充工作面進行地表沉陷預測，結合數(shù)值模擬實驗，深入研究了連采連充技術減緩地表沉陷的機理，得到以下主要結論：

（1）通過室內力學試驗和現(xiàn)場取芯力學試驗，充填體單軸抗壓強度平均為5.093 MPa，煤矸石通過連采連充工藝充填后，充填體的強度超過設計強度2.0 MPa，能夠確保礦井的安全開采。

（2）根據(jù)煤礦二采區(qū)東翼9#煤層實測數(shù)據(jù)，擬合反演出預測參數(shù)，基于等價采高概率積分法，使用巖移數(shù)據(jù)處理系統(tǒng)進行計算，對連采連充工作面開采后進行地表沉陷預測，地表傾斜值極值0.3 mm/m，地表水平變形極值-0.2 mm/m，小于磚混結構建筑Ⅰ級損壞范圍。工作面充填開采后，地面沉陷平緩，無安全隱患。

（3）利用FLAC3D軟件，模擬連采連充工作面開采，分析塑性區(qū)導水裂縫帶發(fā)育高度，可實現(xiàn)水下安全采煤，模擬結果與概率積分法結果對比，驗證連采連充技術可以有效減緩地表移動變形。

參考文獻：

［1］李濤. 豐堠溝水體下采煤安全開采方案研究與探討[J]. 煤炭與化工， 2021， 44（4）： 16-20. DOI： 10.19286/j.cnki.cci.2021.04.005.

[2]馮帆. 淺埋煤層水體下采煤作業(yè)防治水分析[J]. 能源與節(jié)能， 2019（3）： 25-26. DOI： 10.16643/j.cnki.14-1360/td.2019.03.011.

[3]黃玉誠， 董羽， 段仲捷， 等. 似膏體充填水體下采煤技術實踐[J]. 中國礦業(yè)， 2013， 22（7）： 80-82. DOI： 10.3969/j.issn.1004-4051.2013.07.021.

[4]王炯， 郭廣禮， 朱曉峻， 等. 條帶開采與固體充填開采地表沉陷規(guī)律研究[J]. 金屬礦山， 2015（5）： 166-170. DOI： 10.3969/j.issn.1001-1250.2015.05.037.

[5]郭廣禮， 郭凱凱， 張國建， 等. 深部帶狀充填開采復合承載體變形特征研究[J]. 采礦與安全工程學報， 2020， 37（1）： 101-109. DOI： 10.13545/j.cnki.jmse.2020.01.011.

[6]王冰， 郭廣禮， 朱曉峻， 等. 多煤層開采條帶錯動程度對巖層控制的影響[J]. 金屬礦山， 2015（11）：153-157.

[7]殷和健， 査劍鋒， 仲崇武， 等. 固體充填條帶開采采寬與留寬優(yōu)化設計[J]. 煤炭工程， 2019， 51（6）： 34-38. DOI： 10.11799/ce201906007.

[8]李永亮， 路彬， 楊仁樹， 等. 煤礦連采連充式膠結充填采煤技術與典型工程案例[J]. 煤炭學報， 2022， 47（3）： 1055-1071. DOI： 10.13225/j.cnki.jccs.2021.1612.

[9]于濤， 王壯， 劉洋， 等. 長城五礦連采連充自流充填工藝設計[J]. 煤炭工程， 2022， 54（3）：17-21. DOI： 10.11799/ce202203004.

[10]常西坤， 王明國， 王中青， 等. 西部生態(tài)脆弱區(qū)村莊下固廢綠色充填開采技術應用研究[J]. 煤炭工程， 2022， 54（7）：1-6.[JP] DOI： 10.11799/ce202207001.

[11]李浩， 劉音， 王凱， 等. 建筑垃圾骨料-粗粉煤灰基膠結料充填體力學性能試驗研究[J]. 煤礦安全， 2019， 50（12）： 60-63. DOI： [JP]10.13347/j.cnki.mkaq.2019.12.014.

[12]中華人民共和國住房和城鄉(xiāng)建設部.煤礦采空區(qū)建（構）筑物地基處理技術規(guī)范：GB 51180—2016[S].北京：中國計劃出版社，2017.

[13]郭廣禮， 查劍鋒. 礦山開采沉陷學[M]. 徐州： 中國礦業(yè)大學出版社， 2020.

[14]胡炳南， 張華興， 申寶宏. 建筑物、水體、鐵路及主要井巷煤柱留設與壓煤開采指南[M]. 北京： 煤炭工業(yè)出版社， 2017.

[15]GUO G L ， ZHU X J ， ZHA J F ， et al. Subsidence prediction method based on equivalent mining height theory for solid backfilling mining[J].Transactions of Nonferrous Metals Society of China， 2014， 24（10）：3302-3308.

[16]張國建. 巨厚弱膠結覆巖深部開采巖層運動規(guī)律及區(qū)域性控制研究[D]. 徐州： 中國礦業(yè)大學， 2020.

[17]欒元重， 于水， 胡軍偉， 等. 充填開采覆巖導水裂縫帶高度發(fā)育規(guī)律研究[J]. 金屬礦山， 2021（6）： 216-220. DOI： 10.19614/j.cnki.jsks.202106029.

[18]王啟春， 郭廣禮. 村莊下厚煤層綜合機械化矸石充填開采地表沉陷與變形分析[J]. 煤礦安全， 2020， 51（1）： 222-228. DOI： 10.13347/j.cnki.mkaq.2020.01.049.

[19]劉瑞斌. 固體充填綜采地表沉陷預測方法研究[D]. 焦作： 河南理工大學，2014.