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

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

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

中圖分類號：TD823?? 文獻(xiàn)標(biāo)志碼：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ìn)行的采煤工作叫做水體下采煤。在我國，水體下采煤曾引起多次突水安全事故，造成嚴(yán)重的經(jīng)濟(jì)損失和人員傷亡。充填開采是進(jìn)行水體下安全采煤的有效防治措施，是煤礦綠色開采的研究方向之一[1-3]。當(dāng)采空區(qū)采用充填法管理頂板時，可有效抑制采煤引起的導(dǎo)水裂縫帶發(fā)育，減輕工作面支撐壓力產(chǎn)生的礦壓顯現(xiàn)，減緩地表沉陷，實現(xiàn)水體下安全采煤。因此，研究水體下充填采煤技術(shù)具有重要的理論價值和研究意義。

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

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

1 煤礦概況

1.1 采區(qū)概況

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

1.2 連采連充工作面

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

2 充填體物理力學(xué)性能

2.1 室內(nèi)力學(xué)試驗

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

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

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

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

3 連采連充地表移動變形預(yù)測

3.1 預(yù)測方法

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

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

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

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

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

3.2 預(yù)測參數(shù)選取

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

3.3 地表移動變形預(yù)測

基于等價采高概率積分法，采用表4充填開采預(yù)測參數(shù)，使用山東科技大學(xué)特殊開采研究所設(shè)計的巖移數(shù)據(jù)處理系統(tǒng)進(jìn)行計算，得到連采連充地表移動變形數(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。遠(yuǎn)小于《建筑物、水體、鐵路及主要井巷煤柱留設(shè)與壓煤開采指南》[14]中的磚混結(jié)構(gòu)建筑Ⅰ級損壞范圍（水平變形<2.0 mm/m；i<3.0 mm/m）。連采連充工作面充填開采后，觀光塔區(qū)域地表沉陷平緩，無安全隱患。

4 數(shù)值模擬

4.1 數(shù)值模型建立

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

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

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

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

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

利用FLAC3D軟件建立數(shù)值模型，提取在連采連充工作面充填開采后的地表下沉值和地表水平移動值。圖10是利用Origin軟件繪制出的連采連充工作面開采后的地表下沉和地表水平移動曲線圖。本研究比較了兩種不同的方法：一種是基于FLAC3D軟件的數(shù)值模擬，另一種是基于等價采高概率積分法的分析。對這兩種方法的相對誤差進(jìn)行對比，旨在評估它們在解決采煤引起的地表沉陷問題時的一致性和準(zhǔn)確性。具體結(jié)果見表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ù)值模擬結(jié)果，其產(chǎn)生誤差的原因可能是利用FLAC3D進(jìn)行數(shù)值模擬巖層時，其復(fù)雜巖體的裂縫、節(jié)理等因素未充分考慮、未能著重考慮時間效應(yīng)等非采動因素造成的沉陷，造成數(shù)值模擬結(jié)果偏小。假如考慮上述因素，相對誤差將進(jìn)一步減小，概率積分法和數(shù)值模擬的結(jié)果也可以進(jìn)一步驗證出連采連充技術(shù)可以有效減緩地表移動變形。

4.3 地表沉陷機(jī)理解釋

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

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

5 結(jié)論

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

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

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

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

參考文獻(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]黃玉誠， 董羽， 段仲捷， 等. 似膏體充填水體下采煤技術(shù)實踐[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]郭廣禮， 郭凱凱， 張國建， 等. 深部帶狀充填開采復(fù)合承載體變形特征研究[J]. 采礦與安全工程學(xué)報， 2020， 37（1）： 101-109. DOI： 10.13545/j.cnki.jmse.2020.01.011.

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

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

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

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

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

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

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

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

[14]胡炳南， 張華興， 申寶宏. 建筑物、水體、鐵路及主要井巷煤柱留設(shè)與壓煤開采指南[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]張國建. 巨厚弱膠結(jié)覆巖深部開采巖層運(yùn)動規(guī)律及區(qū)域性控制研究[D]. 徐州： 中國礦業(yè)大學(xué)， 2020.

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

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

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