郭文硯 陳清通 賈新果 王志杰

（1.煤炭科學(xué)技術(shù)研究院有限公司，北京 100013；2.煤炭資源高效開(kāi)采與潔凈利用國(guó)家重點(diǎn)實(shí)驗(yàn)室（煤炭科學(xué)研究總院），北京 100013）

礦產(chǎn)資源開(kāi)發(fā)利用造成了大面積土地沉陷，嚴(yán)重影響了礦業(yè)城市的經(jīng)濟(jì)建設(shè)發(fā)展［1］。隨著城鎮(zhèn)化和現(xiàn)代化建設(shè)進(jìn)程加快，土地資源日益緊缺，開(kāi)發(fā)沉陷區(qū)土地作為建設(shè)場(chǎng)地是實(shí)現(xiàn)土地利用的有效途徑［2］。目前對(duì)于煤礦沉陷區(qū)建設(shè)場(chǎng)地穩(wěn)定性評(píng)價(jià)研究較多。張俊英［3］研究了建筑荷載對(duì)采空區(qū)地基穩(wěn)定性的影響；孫慶先等［4］研究了沉陷區(qū)天然氣管道、光伏建設(shè)穩(wěn)定性評(píng)價(jià)方法；賈新果［5］研究了老采空區(qū)殘余沉降預(yù)計(jì)模型；杜明澤等［6］研究了閉坑礦井場(chǎng)地地下水污染影響；劉洋等［7］研究了基于AHP-模糊評(píng)價(jià)的鐵礦采空區(qū)失穩(wěn)預(yù)測(cè)模型。鐵礦礦體因地質(zhì)條件和開(kāi)采方法不同于煤礦，其采空區(qū)穩(wěn)定性分析和地表變形存在一定的差異，不同學(xué)者也進(jìn)行了相關(guān)研究［8-10］。但對(duì)于鐵礦采空區(qū)地表建筑建設(shè)適宜性評(píng)價(jià)的案例相對(duì)較少。

張馬屯鐵礦位于濟(jì)南市歷城區(qū)，1977年投產(chǎn)，2017年8月關(guān)閉。由于濟(jì)南園區(qū)建設(shè)，部分規(guī)劃用地位于張馬屯鐵礦采空區(qū)上方（以下簡(jiǎn)稱(chēng)為“建設(shè)場(chǎng)地”）。該礦礦房采空區(qū)數(shù)量多、分布范圍廣、埋深差異大，采空區(qū)進(jìn)行全尾砂膠結(jié)充填后仍存在剩余空區(qū)，使得地面建筑建設(shè)存在安全隱患。本研究通過(guò)統(tǒng)計(jì)分析、理論分析、經(jīng)驗(yàn)類(lèi)比和現(xiàn)場(chǎng)實(shí)測(cè)，分析剩余空區(qū)頂板穩(wěn)定性，劃分地表不同程度殘余變形區(qū)域，研究建筑物載荷影響深度與覆巖破壞高度的相對(duì)關(guān)系，綜合評(píng)價(jià)地表建筑建設(shè)適宜性，并提出建設(shè)場(chǎng)地建設(shè)安全保障措施，為類(lèi)似金屬礦山沉陷區(qū)土地建筑復(fù)墾提供借鑒。

1 采礦地質(zhì)條件與采空區(qū)分布

張馬屯鐵礦建設(shè)場(chǎng)地下方為第四系土體，厚度為35～55 m，上部為黏性土結(jié)構(gòu)，下部為黏性土與砂性土互層，土體工程地質(zhì)性質(zhì)良好。礦體頂、底板巖體以閃長(zhǎng)巖和蝕變閃長(zhǎng)巖為主，部分為大理巖和石灰?guī)r。主要開(kāi)采Ⅰ號(hào)礦體，采用分段礦房采礦法和全尾砂膠結(jié)充填采空區(qū)，礦體厚度為1.52～70.73 m，階段高12 m，礦房長(zhǎng)度一般為50 m，礦房與礦柱寬度為16 m。根據(jù)相關(guān)調(diào)查，在-200～-360 m水平上共分布有82個(gè)礦房采空區(qū)，如圖1所示。由于采空區(qū)充填可能未及時(shí)或者受地質(zhì)環(huán)境條件的改變和開(kāi)采爆破作業(yè)影響，采空區(qū)仍然有剩余空區(qū)與裂隙存在，通過(guò)瞬變電磁物探也驗(yàn)證了存在一定范圍的空區(qū)。本研究統(tǒng)計(jì)了礦房采出礦石量、面積、充填量和充填率等參數(shù)（表1），對(duì)礦房剩余空區(qū)高度進(jìn)行了估算，可知采空區(qū)平均剩余空區(qū)高度為4.8 m，大于10 m的礦房有14個(gè)，最大為Z3602礦房22.9 m；平均充填率85%，大于85%的有41個(gè)，整體充填效果較好。

根據(jù)采空區(qū)分布，可以將采空區(qū)分為小片采空區(qū)（1～4個(gè)礦房）、連片采空區(qū)（5～7個(gè)礦房）和上下兩層采空區(qū)3種類(lèi)型。小片采空區(qū)跨度最大為64 m，剩余最大空區(qū)高度13.2 m（Z2801～Z2803和Z2807）；連片采空區(qū)最大跨度為103 m，剩余空區(qū)最大高度22.9 m（Z3361～Z3363和Z3601～Z3603）；兩層采空區(qū)間隔最小為8 m（Z3242與Z2761）、最大跨度56 m，下層剩余空區(qū)最大高度19.6 m。連片采空區(qū)和兩層采空區(qū)的穩(wěn)定性相對(duì)較差。

2 采空區(qū)地表建筑建設(shè)適宜性評(píng)價(jià)

2.1 剩余空區(qū)頂板穩(wěn)定性分析

鐵礦采空區(qū)頂板巖體強(qiáng)度相對(duì)較高，有一定的穩(wěn)定性，可通過(guò)頂板安全厚度計(jì)算方法，根據(jù)采空區(qū)形態(tài)、頂板力學(xué)參數(shù)計(jì)算頂板穩(wěn)定所需的安全厚度，進(jìn)而分析空區(qū)穩(wěn)定性［11-12］。

2.1.1 厚跨比法

厚跨比法常用于穩(wěn)定圍巖，根據(jù)近似的水平投影跨度和頂部最薄處厚度，以厚跨比大于0.5作為安全厚度評(píng)價(jià)依據(jù)。則頂板最小安全厚度可按式（1）計(jì)算，

式中，H0為采空區(qū)頂板最小安全厚度，m；L為采空區(qū)跨度，m；K為安全系數(shù)，取3。連片采空區(qū)最大跨度取103 m，計(jì)算可得最小安全厚度為154.5 m。

2.1.2 普氏拱理論計(jì)算法

該理論認(rèn)為采空區(qū)頂板將形成拋物線形的拱帶，空區(qū)上部巖體質(zhì)量由拱承擔(dān)。采用壓力拱和2倍安全系數(shù)來(lái)計(jì)算頂板安全厚度，公式為

式中，φ為巖石內(nèi)摩擦角，（°）；H為剩余空區(qū)最大高度，m；f為普氏系數(shù)，取f=R/ 10；R為巖石的單軸抗壓強(qiáng)度，MPa。采空區(qū)最大跨度取103 m，空區(qū)最大剩余高度取22.9 m，內(nèi)摩擦角取40°，單軸抗壓強(qiáng)度取閃長(zhǎng)巖和大理巖濕式條件下的較小值78.5 MPa，計(jì)算可得頂板最小安全厚度為15.8 m。

2.1.3 結(jié)構(gòu)力學(xué)梁理論計(jì)算法

假定采空區(qū)頂板巖體是一個(gè)兩端固定的平板梁結(jié)構(gòu)，按梁板受彎，以巖層的抗拉強(qiáng)度作為控制指標(biāo)，安全厚度可按式（3）計(jì)算，

式中，b為頂板單位計(jì)算寬度，取1 m；γ為上覆巖體容重，kN/m3；σt為上覆巖體抗拉強(qiáng)度，MPa；P為建筑載荷附加應(yīng)力，kPa，取值為0。巖體抗拉強(qiáng)度可基于Hoek-Brown準(zhǔn)則和GSI指標(biāo)的力學(xué)參數(shù)折算方法來(lái)估算［13］，公式為：

式中，σci為巖塊的單軸抗壓強(qiáng)度，MPa；mb、mi、s為巖體的Hoek-Brown常數(shù)。

頂板巖體以閃長(zhǎng)巖為主，根據(jù)巖體質(zhì)量和經(jīng)驗(yàn)常數(shù)，GSI取65，mi取30，mb取 8.595，s取0.097，單軸抗壓強(qiáng)度取78.5 MPa，可得巖體抗拉強(qiáng)度為886 kPa。采空區(qū)跨度取103 m，巖體容重取27 kN/m3，計(jì)算可得頂板最小安全厚度為161.6 m。

2.1.4 空?qǐng)霰确?/p>

假定空區(qū)頂板為一塊嵌固梁板，受均衡連續(xù)載荷作用，當(dāng)采空區(qū)長(zhǎng)度與跨度之比大于2時(shí)，頂板最小安全厚度可進(jìn)行如下計(jì)算：

式中，P1為爆破載荷，kPa，取0。

代入式（3）中相應(yīng)參數(shù)，計(jì)算可得頂板最小安全厚度為82.3 m。

綜合以上4種計(jì)算方法，連片采空區(qū)剩余空區(qū)頂板安全厚度計(jì)算的較大值為161.6 m，而頂板閃長(zhǎng)巖厚度大于190 m，因此剩余空區(qū)頂板相對(duì)穩(wěn)定，不會(huì)發(fā)生整體塌陷。對(duì)于兩層采空區(qū)，采空區(qū)跨度為16～56 m，采空區(qū)之間頂板厚度為8～24 m，分析可知兩采空區(qū)之間的頂板穩(wěn)定性較差，可能發(fā)生垮落，但由于存在兩層采空區(qū)的范圍較小，垮落巖體碎脹可充填部分剩余空間。按上述理論方法分析可知，上層采空區(qū)剩余空區(qū)頂板相對(duì)穩(wěn)定，也不會(huì)發(fā)生整體塌陷。

2.2 建設(shè)場(chǎng)地地表殘余移動(dòng)變形計(jì)算

概率積分法具有參數(shù)容易確定、實(shí)用性強(qiáng)等優(yōu)點(diǎn)，是地表移動(dòng)變形計(jì)算最常用的方法之一。一般認(rèn)為，地表殘余移動(dòng)變形與移動(dòng)期內(nèi)移動(dòng)變形規(guī)律一致［5］，因此，張馬屯鐵礦剩余空區(qū)地表殘余移動(dòng)變形可以類(lèi)比概率積分法計(jì)算。根據(jù)礦房尺寸、剩余空洞高度、埋深等條件，沉陷預(yù)計(jì)按照最危險(xiǎn)情況下所有礦房空洞被垮塌壓實(shí)考慮，結(jié)合覆巖巖性、開(kāi)采時(shí)間與房柱式開(kāi)采方法原理，單個(gè)礦房采空區(qū)殘余下沉系數(shù)取0.2，連片采空區(qū)適當(dāng)增大下沉系數(shù)。根據(jù)移動(dòng)變形軟件（CESSDG系統(tǒng)）計(jì)算結(jié)果，地表殘余下沉為0～213 mm，傾斜量為3.4～+2.5 mm/m，水平變形為-3.5～+3.6 mm/m。參考“三下采煤”規(guī)范［14］中磚混結(jié)構(gòu)建筑物損壞等級(jí)劃分標(biāo)準(zhǔn)，將地表下沉小于10 mm的區(qū)域劃分為無(wú)影響區(qū)（A區(qū)），將水平變形不大于2.0 mm/m或傾斜變形不大于3 mm/m的區(qū)域劃分為輕微變形區(qū)（B區(qū)），將水平變形大于2.0 mm/m且傾斜變形大于3 mm/m的區(qū)域劃分為Ⅱ級(jí)輕度損壞變形區(qū)（C區(qū)），則A區(qū)面積占建設(shè)場(chǎng)地8.6%，B區(qū)占75.8%，C區(qū)占15.6%，變形分區(qū)見(jiàn)圖2。因此，建設(shè)場(chǎng)地內(nèi)大部分區(qū)域地表建筑屬于I級(jí)輕微損壞變形，少部分區(qū)域受到Ⅱ級(jí)輕度損壞變形。

2.3 覆巖破壞高度與建筑載荷影響深度的關(guān)系

礦體開(kāi)采后，頂板一定范圍內(nèi)巖層發(fā)生垮落、斷裂破壞，巖體強(qiáng)度降低，存在裂縫和離層，如果建筑荷載附加應(yīng)力影響到該破壞區(qū)域，將造成巖體進(jìn)一步變形和沉降，可能導(dǎo)致建筑物地基破壞。礦房開(kāi)采相比于煤層長(zhǎng)壁開(kāi)采的范圍小，頂板巖體強(qiáng)度好，其覆巖破壞范圍通常小于長(zhǎng)壁開(kāi)采，但考慮部分剩余空區(qū)高度較大，因此鐵礦開(kāi)采覆巖破壞高度計(jì)算可類(lèi)比“三下采煤”規(guī)范［14］中厚煤層分層開(kāi)采堅(jiān)硬頂板條件下的“兩帶”高度經(jīng)驗(yàn)公式，

式中，Hli為導(dǎo)水裂縫帶高度，m；∑M為累計(jì)采厚，m。

選取連片采空區(qū)進(jìn)行計(jì)算，平均剩余空洞高度取9.5 m，代入上式中計(jì)算并取較大值為102.5 m。

根據(jù)《建筑地基基礎(chǔ)設(shè)計(jì)規(guī)范》（GB 50007—2011）［15］，存在采空區(qū)時(shí)，以地基中建筑載荷附加應(yīng)力等于相應(yīng)位置處地基土層自重應(yīng)力10%時(shí)的深度作為載荷影響深度。取擬建建筑物載荷為120、240、360 kPa，單棟建筑平面尺寸為60 m×20 m，第四系土層厚度取40 m，容重取20 kN/m3，基巖層容重取25 kN/m3，基礎(chǔ)埋深為5 m，可計(jì)算出不同建筑載荷的影響深度為29、39、44 m（圖3）。

要保證建設(shè)場(chǎng)地穩(wěn)定性，應(yīng)使建筑載荷附加應(yīng)力不對(duì)采空區(qū)破壞巖體產(chǎn)生影響，則采空區(qū)的臨界開(kāi)采深度（H臨）應(yīng)滿足：

式中，hz為載荷影響深度，取44 m；a為保護(hù)層厚度，取10 m，則臨界開(kāi)采深度應(yīng)大于156.5 m。

建設(shè)場(chǎng)地地表標(biāo)高+35 m、開(kāi)采標(biāo)高-200～-360 m、礦房高度小于60 m，則實(shí)際采深大于臨界開(kāi)采深度，因此地表新增360 kPa建筑荷載對(duì)采空區(qū)穩(wěn)定性無(wú)影響。

2.4 地面沉降監(jiān)測(cè)分析

張馬屯鐵礦地表布設(shè)有沉降監(jiān)測(cè)點(diǎn)，在建設(shè)場(chǎng)地內(nèi)部和附近的部分監(jiān)測(cè)點(diǎn)位置見(jiàn)圖1，監(jiān)測(cè)時(shí)間從2015年2月到2019年10月，測(cè)點(diǎn)累計(jì)沉降曲線見(jiàn)圖4?？梢钥闯?，J5測(cè)點(diǎn)最大累計(jì)沉降為11 mm，主要由于靠近連片采空區(qū)，沉降較大。該礦2017年8月31日停產(chǎn)，從2018年5月之后，測(cè)點(diǎn)沉降趨于穩(wěn)定。地表整體沉降較小，采空區(qū)穩(wěn)定性較好。

結(jié)合地表沉降監(jiān)測(cè)和移動(dòng)變形預(yù)計(jì)結(jié)果，該區(qū)域地表最大傾斜值為3.4 mm/m?？紤]建筑自身荷載作用，根據(jù)《建筑地基基礎(chǔ)設(shè)計(jì)規(guī)范》（GB 50007—2011）［15］規(guī)定，當(dāng)多層和高層建筑高度小于24 m時(shí)，建筑整體傾斜允許值為4 mm/m，建筑高度為24～60 m時(shí)，傾斜允許值為3 mm/m。因此，根據(jù)預(yù)測(cè)的地表最大傾斜值，建議該區(qū)域建筑高度限制在36 m以下，不適宜建設(shè)超高層建筑。

2.5 地表建筑建設(shè)適宜性評(píng)價(jià)

綜上分析，采空區(qū)剩余空區(qū)頂板穩(wěn)定，不會(huì)發(fā)生大面積塌陷；建設(shè)場(chǎng)地地表75.8%區(qū)域建筑受到I級(jí)輕微損壞變形，15.6%區(qū)域受到Ⅱ級(jí)輕度損壞變形，危險(xiǎn)程度較低；地表新增360 kPa建筑荷載對(duì)采空區(qū)穩(wěn)定性無(wú)影響；通過(guò)地面沉降長(zhǎng)期監(jiān)測(cè)，地表整體沉降較小，采空區(qū)穩(wěn)定性較好。因此，建設(shè)場(chǎng)地進(jìn)行一般工程建設(shè)（建筑高度36 m以下）的適宜性較好。

3 建設(shè)場(chǎng)地安全保障措施

（1）奧陶紀(jì)石灰?guī)r裂隙含水層為礦井主要充水因素，礦區(qū)采取注漿大帷幕對(duì)該含水層進(jìn)行封堵，有效降低了礦井涌水量。但長(zhǎng)期滲透使帷幕注漿結(jié)石體受到侵蝕，帷幕失效可能導(dǎo)致礦井涌水量增加，涌水侵蝕會(huì)降低礦柱及頂板巖體強(qiáng)度，加速采空區(qū)垮塌，影響采空區(qū)穩(wěn)定性。因此，建設(shè)過(guò)程中需進(jìn)一步排查不良封閉鉆孔，加強(qiáng)對(duì)帷幕內(nèi)外水文觀測(cè)孔水位及井下涌水量在線監(jiān)測(cè)。

（2）建設(shè)場(chǎng)地不適宜建設(shè)（超）高層建筑，建議該區(qū)域建筑高度不超過(guò)36 m，并增加基礎(chǔ)抗變形能力，建議在連片采空區(qū)和剩余空區(qū)高度較大的礦房及短礦柱上方地面增加地表及建筑物沉降觀測(cè)點(diǎn)，定期進(jìn)行觀測(cè)。

（3）對(duì)于下方繼續(xù)使用的巷道，應(yīng)加強(qiáng)頂板監(jiān)測(cè)及巷道維修保護(hù)，若巷道不再繼續(xù)使用，建議對(duì)廢棄巷道進(jìn)行充填支護(hù)，避免巷道變形影響采空區(qū)穩(wěn)定。

4 結(jié) 論

（1）綜合比較厚跨比、普氏拱理論、結(jié)構(gòu)力學(xué)梁理論和空?qǐng)霰?種頂板安全厚度計(jì)算結(jié)果和剩余空區(qū)頂板實(shí)際厚度，認(rèn)為張馬屯鐵礦剩余空區(qū)頂板厚度滿足安全要求，不會(huì)發(fā)生整體塌陷。

（2）通過(guò)剩余空區(qū)殘余變形預(yù)計(jì)，建設(shè)場(chǎng)地地表75.8%區(qū)域內(nèi)建筑受到I級(jí)輕微損壞變形，15.6%區(qū)域內(nèi)建筑受到Ⅱ級(jí)輕度損壞變形，危險(xiǎn)程度較低?，F(xiàn)場(chǎng)觀測(cè)地表5 a累計(jì)最大沉降11 mm，剩余空區(qū)穩(wěn)定性較好。

（3）通過(guò)載荷影響深度分析，建設(shè)場(chǎng)地地表新增360 kPa建筑荷載對(duì)剩余空區(qū)穩(wěn)定性無(wú)影響，場(chǎng)地進(jìn)行工程建設(shè)（建筑高度36 m以下）的適宜性較好。