摘" 要：厚大礦體資源的開采面臨施工難、采動(dòng)地壓大的難題，而采用以廂式充填法為主的分層條帶充填開采在杉樹埡磷礦實(shí)現(xiàn)順利應(yīng)用。為研究厚大礦體充填條帶開采巖層損傷與地表沉降，采用地質(zhì)雷達(dá)的方法，監(jiān)測得到厚大磷礦層分層條帶充填開采產(chǎn)生的頂板損傷深度在8.5 m左右，采動(dòng)影響范圍較低，廂式充填采礦法在厚大礦體開采中應(yīng)用效果較好。厚大礦體的礦柱在5～10 m范圍內(nèi)存在明顯的巖層裂隙，而充填體的裂隙節(jié)理主要分布在17 m范圍內(nèi)，超過17 m深度的充填體裂隙數(shù)目明顯降低，表現(xiàn)出完整性的提升。通過DS-InSAR地表位移監(jiān)測方法，研究得到2019—2023年，401段和402段的地表沉陷最大值為47 mm，且總體上變化不均勻而出現(xiàn)明顯波動(dòng)性，表明厚大礦體開采并未對地表沉陷造成明顯影響。

關(guān)鍵詞：磷礦；地質(zhì)雷達(dá)；采動(dòng)應(yīng)力；DS-InSAR；地表位移

中圖分類號：TD326" " " "文獻(xiàn)標(biāo)志碼：A" " " " " 文章編號：2095-2945（2025）05-0112-05

Abstract： The mining of thick and large ore body resources faces difficulties in construction and high mining ground pressure. However， layered strip filling mining based on box filling method has been successfully applied in Shanshuya Phosphate Mine. In order to study the rock formation damage and surface subsidence in thick and large ore bodies filled strip mining， the geological radar method was used to monitor and determine that the depth of roof damage caused by layered strip filling mining in thick and large phosphorus ore layers was about 8.5 m， and the mining influence range was relatively low. The box-filling mining method has a good application effect in the mining of thick and large ore bodies. There are obvious rock fractures in the range of 5 to 10 m in the pillars of thick and large ore bodies， while the crack joints of the filling body are mainly distributed in the range of 17 m. The number of cracks in the filling body beyond 17 m is significantly reduced， showing an improvement in integrity. Through the DS-InSAR ground displacement monitoring method， it is found that from 2019 to 2023， the maximum value of surface subsidence in sections 401 and 402 is 47 mm， and the overall change is uneven and there is obvious fluctuation， indicating that the mining of thick and large ore bodies has not significantly affected surface subsidence.

Keywords： phosphate ore; geological radar; mining stress; DS-InSAR; ground displacement

磷礦資源是世界上重要的戰(zhàn)略資源，屬于不可再生的資源，具有不可代替性[1]。推行安全高效的磷礦開采技術(shù)，努力提高資源的回采率和開發(fā)利用率，是當(dāng)前我國資源開發(fā)利用的原則之一[2]。在目前磷礦開采過程中，相較于其他礦體，受限于自身所固有的開采技術(shù)條件及礦體的賦存條件和其產(chǎn)狀的特殊性，在開采過程中往往會(huì)使巖層結(jié)構(gòu)發(fā)生破壞損傷從而引起覆巖破斷導(dǎo)致大量節(jié)理裂隙產(chǎn)生，延伸至地表后發(fā)生地表沉降，嚴(yán)重影響生態(tài)環(huán)境保護(hù)[3-6]。因此，在開采過程中探索安全可靠、高效合理的監(jiān)測技術(shù)，對于預(yù)防潛在的地質(zhì)災(zāi)害及生態(tài)環(huán)境修復(fù)具有重要意義。

隨著測繪技術(shù)的快速發(fā)展，國內(nèi)學(xué)者針對采礦引起的巖層損傷及地表沉陷進(jìn)行了諸多深入的研究[7-9]。敖嫩等[10]采用現(xiàn)場實(shí)地調(diào)查與統(tǒng)計(jì)分析的手段，分析了鄂爾多斯地區(qū)地裂縫形成的原因，提出了解決地裂縫問題的防治方法；朱鵬程[11]基于無人機(jī)測量技術(shù)提出了地表水平移動(dòng)變形的監(jiān)測方法，在南屯煤礦應(yīng)用效果良好；周安民[12]針對開采沉陷數(shù)據(jù)采集時(shí)序問題，提出了基于濾波算法+遙感測繪技術(shù)的新型監(jiān)測方法，準(zhǔn)確率可達(dá)95%以上；張麗媛等[13]采用無人機(jī)和InSAR多源手段聯(lián)合應(yīng)用的方法，實(shí)現(xiàn)了對礦區(qū)開采沉陷的有效監(jiān)測；楊中輝等[14]運(yùn)用模糊層次分析法對黃土溝壑區(qū)開采沉陷進(jìn)行了綜合評價(jià)，具有較高推廣價(jià)值；此外還有學(xué)者利用SBAS-InSAR技術(shù)對開采沉陷進(jìn)行了分析研究[15-16]；李濤等[17]構(gòu)建了基于DBD-Net的InSAR開采沉陷盆地監(jiān)測網(wǎng)絡(luò)；惠甜甜等[18]將SBAS-InSAR技術(shù)和LSTM-GS模型相結(jié)合獲取了時(shí)序累積沉降值。

綜上所述，學(xué)者們對于開采引起的覆巖損傷與地表沉陷監(jiān)測已取得了諸多成果，但對于不同礦區(qū)由于其地質(zhì)條件的特殊性、復(fù)雜性，在監(jiān)測方面缺乏聯(lián)合有效的手段對其進(jìn)行分析研究。本文以杉樹埡磷礦充填條帶開采為工程背景，通過聯(lián)合地質(zhì)雷達(dá)井下監(jiān)測和DS-InSAR地表沉陷監(jiān)測，為礦區(qū)厚大磷礦體充填開采產(chǎn)生的巖層損傷與地表沉降進(jìn)行透明化、精細(xì)化分析，并為相關(guān)厚大礦體開采工程提供參考與指導(dǎo)。

1" 工程概況

1.1" 厚大礦體開采概況

杉樹埡磷礦的開采面臨著礦床厚度增大的挑戰(zhàn)，使得傳統(tǒng)房柱式開采方法的應(yīng)用受到了限制。為此，該礦進(jìn)行開采方法的革新，針對厚大礦體采用分層充填條帶開采方法。圖1為該礦正在開采的四采區(qū)平面圖。

1.2" 分層充填條帶開采

針對四采區(qū)水平或緩傾斜厚大礦體，采用分層充填條帶開采，即廂式充填采礦法，其采礦作業(yè)流程如圖2所示。

針對四采區(qū)厚大磷礦層，先在礦塊兩邊和中間分別施工上、下部通道，條帶礦房回采工藝采用分層開采，即先由上部通道進(jìn)行切頂，再由下部通道進(jìn)行降低完成全礦層回采。礦塊內(nèi)膠結(jié)充填條帶和干式充填條帶間隔布置，分2步回采，先回采膠結(jié)充填條帶，待回采完所有膠結(jié)充填條帶，進(jìn)行膠結(jié)充填，充填體達(dá)到一定強(qiáng)度后再回采干式充填條帶，待回采完干式充填條帶，進(jìn)行干式充填。采用鑿巖臺車鑿巖，每次回采步距2.5～3 m，回采崩下礦石由鏟運(yùn)機(jī)搬運(yùn)到上、下部通道附近的集中裝載點(diǎn)裝入坑內(nèi)卡車，經(jīng)由中段運(yùn)輸巷道運(yùn)至采區(qū)溜井，由此完成厚大礦體的開采和采空區(qū)的充填作業(yè)。

2" 巖層采動(dòng)損傷地質(zhì)雷達(dá)監(jiān)測

2.1" 地質(zhì)雷達(dá)監(jiān)測原理

隨著開采礦體的厚度增大，其開采對巖層的損傷范圍需要進(jìn)一步監(jiān)測，尤其是厚大礦體開采后頂板的巖層損傷范圍。同時(shí)，由于采用充填體來置換條帶的礦石，充填體的支撐性能需要進(jìn)行監(jiān)測，可對礦柱支撐幫部和充填體支撐幫部進(jìn)行對比監(jiān)測，評估置換前后的性能。此處，采用基于高頻電磁波技術(shù)的地質(zhì)雷達(dá)，包括括激發(fā)裝置100 MHz天線以及觸摸屏數(shù)字視頻記錄器（DVL）。

2.2" 頂板巖層采動(dòng)損傷監(jiān)測

在厚大礦體開采區(qū)域進(jìn)行測點(diǎn)的選擇，選取至少3處頂板相較平整的區(qū)域進(jìn)行布置監(jiān)測。需保障測點(diǎn)范圍處于受礦房開采的擾動(dòng)區(qū)域內(nèi)，且盡量選擇采場的中部區(qū)域。綜合地質(zhì)雷達(dá)頂板損傷區(qū)監(jiān)測結(jié)果，分析評估厚大礦體的開采擾動(dòng)對于頂板的損傷范圍，監(jiān)測的現(xiàn)場圖與結(jié)果如圖3所示。

通過對頂板30 m深度范圍進(jìn)行地質(zhì)雷達(dá)探測，測線0～3 m部分由于頂板周圍有金屬，數(shù)據(jù)受干擾程度較大；通過觀察檢測結(jié)果發(fā)現(xiàn)巖層裂隙集中在4～8.5 m處；超過9 m的巖層基本完整性和致密性較好，未受到采動(dòng)損傷。由此說明，厚大礦體開采產(chǎn)生的頂板損傷深度在8.5 m左右。

2.3" 巷幫巖層、充填體采動(dòng)損傷監(jiān)測

條帶開采充填開采前期是需要礦柱支撐頂板，而后期會(huì)由充填體置換礦柱，由充填體來支撐頂板。因此，對比礦柱和充填體在采動(dòng)影響下的損傷范圍將利于分析和評價(jià)充填體的性能。在現(xiàn)場由地質(zhì)雷達(dá)探測30 m范圍的礦柱和充填體，如圖4所示。

礦柱受采動(dòng)損傷的范圍明顯分布在10 m內(nèi)，其淺部同樣受到支護(hù)體、金屬網(wǎng)等干擾，監(jiān)測到的裂隙分布不明顯，結(jié)果顯示5～10 m范圍內(nèi)存在明顯的巖層裂隙，而超過10 m深度的巖層基本無明顯的裂隙帶，巖層的致密性未受到影響。充填體的整個(gè)監(jiān)測范圍的致密性相較原生礦柱明顯較差，且裂隙節(jié)理主要分布在17 m范圍內(nèi)，越淺部的區(qū)域裂隙數(shù)目越多，而超過17 m深度的充填體裂隙數(shù)目明顯降低，表現(xiàn)出完整性的提升。通過兩者采動(dòng)損傷范圍監(jiān)測結(jié)果的直觀對比，說明充填體相較原生礦柱受到厚大礦體開采損傷的程度高，且充填體邊界區(qū)域（17 m之內(nèi)）的支撐性能不能完全比肩礦柱，深部區(qū)域的充填體性能顯著提升，可以實(shí)現(xiàn)較好的支撐效果。

3" 厚大礦體開采地表沉陷監(jiān)測

3.1" DS-InSAR地表沉陷監(jiān)測原理

由于時(shí)序InSAR技術(shù)在變形計(jì)算上應(yīng)用穩(wěn)定相位的散射體，實(shí)現(xiàn)了高相干性，這可應(yīng)對常規(guī)技術(shù)的時(shí)空相干性缺失及延遲的問題，但其在礦山采空區(qū)地表中植被較少的低相干地區(qū)監(jiān)測點(diǎn)密度低，在應(yīng)用上遇到阻礙。DS-InSAR技術(shù)可對其缺陷進(jìn)行彌補(bǔ)，它應(yīng)用分布式散射體，并集成同質(zhì)像元識別和相位優(yōu)化，可較好地應(yīng)用于礦上開采后采空區(qū)地表沉陷量的監(jiān)測，其生成時(shí)間序列形變圖的步驟原理如圖5所示。

3.2" 開采沉陷監(jiān)測區(qū)域選擇

董家河河床兩岸下伏存在四采區(qū)401和402 兩個(gè)塊段，2019—2023年兩個(gè)塊段進(jìn)行著厚層礦體高強(qiáng)度開采、充填作業(yè)。通過井下平面圖的位置信息，采用DS-InSAR調(diào)取董家河附近401和402 兩個(gè)塊段地表區(qū)域近5年的天通衛(wèi)星數(shù)據(jù)，即2019—2023年的地表位移數(shù)據(jù)。對董家河沿岸采區(qū)進(jìn)行井下采掘工程平面圖、衛(wèi)星地表地勢圖和累計(jì)變形數(shù)據(jù)圖進(jìn)行對照，如圖6所示。

將礦區(qū)四采區(qū)的401分段和402分段在衛(wèi)星云圖與地表沉陷監(jiān)測數(shù)據(jù)圖上進(jìn)行標(biāo)注，可以得到401段和402段區(qū)域的地表沉陷數(shù)據(jù)。在圖6的監(jiān)測區(qū)域累計(jì)變形量圖中，出現(xiàn)變形量的點(diǎn)均會(huì)被展示，包括抬升量和下沉量。監(jiān)測可得，2019—2023年，兩個(gè)區(qū)域地表的局部抬升與局部下沉并存，而重點(diǎn)關(guān)注的是發(fā)生下沉的區(qū)域。

3.3" 區(qū)域地表位移監(jiān)測分析

調(diào)取四采區(qū)401分段區(qū)域2019—2023年的累計(jì)變形量，其包含抬升和沉降量。通過結(jié)果對比分析出最大位移沉降數(shù)據(jù)點(diǎn)，并對數(shù)據(jù)進(jìn)行按月份的統(tǒng)計(jì)，得到位移量時(shí)間變化過程，如圖7所示。

401分段的地表區(qū)域內(nèi)，西側(cè)存在抬升區(qū)域和位移不變區(qū)域，而中部和東側(cè)均為發(fā)生沉陷區(qū)。根據(jù)數(shù)據(jù)的定位，得到最大沉降點(diǎn)的最大沉降量為47 mm，平均值為24 mm，總體上呈現(xiàn)緩慢波動(dòng)下沉的情況，且趨勢逐漸變緩。最大沉降點(diǎn)的位置在401分段對應(yīng)地表的東北區(qū)域，為河流沿岸公路邊坡側(cè)，而此處地下的礦體并未進(jìn)行大規(guī)模充分的開采，其沉降的原因與降雨沖刷和滾石滑坡等情況相關(guān)。

同樣調(diào)取四采區(qū)402分段區(qū)域在2019—2023年的累計(jì)變形量，按月統(tǒng)計(jì)和分析出最大位移沉降數(shù)據(jù)點(diǎn)的變化過程，如圖8所示。

402分段的地表區(qū)域內(nèi)，西側(cè)、北側(cè)和東南側(cè)存在明顯的地表抬升，而中部及南側(cè)存在沉陷區(qū)和位移不變區(qū)域。根據(jù)數(shù)據(jù)的定位，最大沉降點(diǎn)的變化波動(dòng)較大，正位移與負(fù)位移交替，最大沉降量為32 mm，平均值為15 mm。由此說明，此處地表為山體，其表面的抬升與降低交替存在，說明其變化與地下采礦活動(dòng)的相關(guān)性不充分。

4" 結(jié)論

1）地質(zhì)雷達(dá)監(jiān)測得到厚大磷礦層分層條帶充填開采產(chǎn)生的頂板損傷深度在8.5 m左右，采動(dòng)影響范圍較低，廂式充填采礦法在厚大礦體開采中應(yīng)用效果較好。

2）厚大礦體的礦柱在5～10 m范圍內(nèi)存在明顯的巖層裂隙，而充填體的裂隙節(jié)理主要分布在17 m范圍內(nèi)，超過17 m深度的充填體裂隙數(shù)目明顯降低，表現(xiàn)出完整性的提升。

3）401分段和402分段的地表沉陷最大值為47 mm，且總體上變化不均勻而出現(xiàn)明顯波動(dòng)性，表明厚大礦體開采并未對地表沉陷造成明顯影響。

參考文獻(xiàn)：

[1] 賈晗，劉軍省，焦森，等.長江經(jīng)濟(jì)帶磷礦資源開發(fā)與生態(tài)保護(hù)現(xiàn)狀分析及對策建議[J].中國礦業(yè)，2021，30（2）：67-72.

[2] 高林，蘆慶和，梁維，等.上向分段充填開采覆巖移動(dòng)變形及影響范圍預(yù)測——以貴州某磷礦為例[J].化工礦物與加工，2023，52（9）：59-65.

[3] 陳冉麗，王紅軍，刁鑫鵬，等.InSAR形變監(jiān)測在開采沉陷損害鑒定及穩(wěn)定性評價(jià)中的應(yīng)用[J].測繪通報(bào)，2023（12）：106-111.

[4] 武毅藝，何滿潮，李輝，等.深部多形式應(yīng)用交岔點(diǎn)圍巖失穩(wěn)機(jī)理及能量演化研究（英文）[J].Journal of Central South University， 2024，31（3）：890-911.

[5] 余華婷，徐連滿，李瑋哲，等.地表開采沉陷監(jiān)測技術(shù)研究[J].科技創(chuàng)新與應(yīng)用，2023，13（20）：183-186.

[6] 何滿潮，武毅藝，高玉兵，等.深部采礦巖石力學(xué)進(jìn)展[J].煤炭學(xué)報(bào)，2024，49（1）：75-99.

[7] 張明媚，薛永安.多煤層重復(fù)開采地表下沉系數(shù)演變研究[J].煤炭工程，2024，56（2）：152-156.

[8] 楊澤發(fā)，馬澤林，喬思宇.基于InSAR技術(shù)的開采沉陷影響邊界劃定方法[J].金屬礦山，2023（1）：119-125.

[9] 夏元平.基于InSAR/GIS的礦區(qū)地下非法開采監(jiān)測關(guān)鍵技術(shù)研究[J].測繪學(xué)報(bào)，2024，53（2）：394.

[10] 敖嫩，王漢元，劉宏磊，等.鄂爾多斯高強(qiáng)度煤炭開采區(qū)地裂縫發(fā)育特征[J].能源與環(huán)保，2023，45（3）：251-256，262.

[11] 朱鵬程.基于無人機(jī)的礦區(qū)地表水平移動(dòng)監(jiān)測方法[J].地理空間信息，2024，22（3）：16-19，23.

[12] 周安民.基于遙感測繪技術(shù)的礦山開采沉陷監(jiān)測方法研究[J].采礦技術(shù)，2022，22（6）：69-72.

[13] 張麗媛，張兆江，劉闊，等.基于多源數(shù)據(jù)的礦區(qū)形變監(jiān)測研究——以邯鄲萬年礦為例[J].測繪與空間地理信息，2023，46（9）：31-35.

[14] 楊中輝，李培現(xiàn)，郝登程，等.黃土溝壑區(qū)煤礦開采地面沉陷災(zāi)害危險(xiǎn)性評價(jià)[J].能源與環(huán)保，2024，46（3）：239-244.

[15] 李博，張雨欣，白喜成，等.紅柳林煤礦采空區(qū)地面沉陷InSAR監(jiān)測與分析[J].中國礦業(yè)，2023，32（8）：55-62.

[16] 周健，周大偉，安士凱，等.基于SBAS-InSAR的西部開采沉陷超前影響角反演研究[J].金屬礦山，2023（4）：168-173.

[17] 李濤，鄒英杰，范洪冬，等.基于DBD-Net的InSAR礦區(qū)開采沉陷盆地檢測方法[J].煤礦安全，2024，55（4）：177-186.

[18] 惠甜甜，劉長星，王圣杰，等.基于SBAS-InSAR技術(shù)和LSTM-GS模型的礦區(qū)開采沉陷預(yù)測[J].地理空間信息，2022，20（8）：13-17.