劉王政　張云鵬　劉愛興,2　鄭文龍

(1．華北理工大學(xué)礦業(yè)工程學(xué)院；2．徐州鐵礦集團(tuán)有限公司；3.河北鋼鐵集團(tuán)司家營鐵礦)

?

崩落圍巖處理采空區(qū)對尾礦庫的影響研究

劉王政1張云鵬1劉愛興1,2鄭文龍3

(1．華北理工大學(xué)礦業(yè)工程學(xué)院；2．徐州鐵礦集團(tuán)有限公司；3.河北鋼鐵集團(tuán)司家營鐵礦)

摘要某鐵礦地下采空區(qū)的地表為尾礦庫的一部分，為確保尾礦庫穩(wěn)定，采用崩落充填法處理后運用“三帶”理論、巖石移動帶理論及數(shù)值模擬對采空區(qū)進(jìn)行研究。結(jié)果表明，崩落處理采空區(qū)后，地表水平位移和豎直位移均處在標(biāo)準(zhǔn)范圍內(nèi)，地表沉降裂隙帶距尾礦庫最小距離為127 m，采空區(qū)處理后能夠保證地表尾礦庫穩(wěn)定。但隨著該礦深部礦體的進(jìn)一步開采，應(yīng)加強(qiáng)對地表裂隙帶范圍的監(jiān)控，以保證尾礦庫的安全。

關(guān)鍵詞崩落圍巖采空區(qū)尾礦庫穩(wěn)定性分析

隨著時間的推移，滯留的采空區(qū)將誘使上覆圍巖發(fā)生沉降，對深部礦體的安全開采產(chǎn)生嚴(yán)重影響。上覆圍巖沉降發(fā)展到地表形成一定區(qū)域的巖石移動帶，在移動帶內(nèi)建構(gòu)物的穩(wěn)定性將受到嚴(yán)重影響，給礦山及附近居民的生活帶來巨大的安全隱患。

國內(nèi)外許多專家研究了礦山地下開采引起地表沉陷及對地表建筑物的影響，并取得了較多成果。前蘇聯(lián)專家阿維爾申運用塑性理論，得出礦體采空后引起地表下沉盆地剖面方程，形成了水平移動與地面傾斜成正比的結(jié)論[1]；梁力等人利用ANSYS有限元分析軟件，對尾礦壩穩(wěn)定性進(jìn)行了數(shù)值模擬分析[2]；李宏儒等人采用有效應(yīng)力算法對壩體邊坡的穩(wěn)定性進(jìn)行了數(shù)值模擬分析[3]；李國政等運用結(jié)構(gòu)可靠度指標(biāo)，對壩體穩(wěn)定性進(jìn)行了評價研究[4]。然而，崩落法處理后采空區(qū)對地表尾礦庫穩(wěn)定性的影響研究相對較少。

以某地下鐵礦為例，該礦采區(qū)共有3條礦體：由北向南分別為Fe1#、Fe2#和Fe3#礦脈，其中Fe1#正在開采，220 m水平以上Fe2#和Fe3#礦脈已采完，形成采空區(qū)，采用崩落法治理采空區(qū)并回收礦柱。在采區(qū)上盤存在即將廢棄的尾礦庫，庫區(qū)面積約6.25萬m2，庫容約250萬m3。隨著礦床開采深度的逐漸加大，采空區(qū)及其周圍可能產(chǎn)生下沉或塌陷，并誘發(fā)裂隙。由于尾礦庫的一部分位于采空區(qū)巖石移動范圍內(nèi)，本文采用三帶理論、巖石移動帶分析法和數(shù)值模擬，對崩落圍巖處理的采空區(qū)進(jìn)行研究分析，以確保礦山正常生產(chǎn)及尾礦庫的穩(wěn)定。

1采空區(qū)對地表影響評價方法

1.1三帶理論分析方法

礦床開采后，空區(qū)圍巖及上覆巖層在應(yīng)力重新分配的過程中運動。為研究其移動破壞規(guī)律，由下至上依次將巖層移動劃分為冒落帶、裂隙帶、彎曲下沉帶，如圖1所示。通過計算分析得出地表建構(gòu)物所處的巖層移動帶，并判斷其穩(wěn)定性。

圖1　采空區(qū)上覆巖層三帶示意

(1)冒落帶?？諈^(qū)頂板所屬巖層發(fā)生彎曲、變形以及斷裂、冒落等破壞。冒落帶高度為：

(1)

式中,m為采高，m；kp為巖體松散系數(shù)，為1.25～1.50。

(2)裂隙帶。巖層沿層理方向發(fā)生裂隙進(jìn)而產(chǎn)生破壞，下沉量及變形量大，位于兩帶之間。

(3)彎曲下沉帶。位于裂隙帶上方至地表，該范圍內(nèi)巖層受自重影響，發(fā)生層面法向彎曲下沉且不斷裂。

1.2巖石移動帶分析方法

地下采空區(qū)促使巖體原有應(yīng)力的平衡發(fā)生改變，在新應(yīng)力平衡過程中，采空區(qū)周圍巖體發(fā)生的形變和移動會逐漸延展至地表，使地表形成一定范圍的移動帶和塌陷帶，空區(qū)底部與移動帶邊緣連線和水平面所形成的夾角叫做移動角。在移動帶以內(nèi)的地表建筑物及施工工程，都可能會受到巖體移動的影響而發(fā)生破壞，如圖2所示。通過計算巖石移動角和塌落角，判斷地表建構(gòu)物是否處于移動帶內(nèi)。

1.3巖石移動帶分析方法

隨著計算機(jī)技術(shù)的發(fā)展，數(shù)值模擬方法在地質(zhì)環(huán)境復(fù)雜、不利于現(xiàn)場試驗的地下工程研究中得到廣泛應(yīng)用。在已知礦巖體物理參數(shù)的條件下，通過建立數(shù)值模型對采空區(qū)及上覆圍巖進(jìn)行模擬研究，分析地表形變量及應(yīng)力分布圖，進(jìn)而得出空區(qū)的穩(wěn)定性及對地表建構(gòu)物的影響程度。

圖2　巖石移動帶與塌陷帶范圍

2穩(wěn)定性綜合分析

2.1三帶理論分析

采空區(qū)崩落的圍巖為散體，壓實后會在頂部出現(xiàn)少量空區(qū)，每個階段采空區(qū)崩落的巖石充滿率按95%計算。尾礦庫位于勘探線3-3和3-7附近，故選取這兩個剖面為研究對象，結(jié)合礦山開采現(xiàn)狀及巖性，取巖體松散系數(shù)為1.3，裂隙帶高度取采高的25倍，計算各勘探線巖體的冒落帶和裂隙帶高度。如表1所示。

表1　各勘探線冒落帶高度計算結(jié)果

由表1可知，在3-3線剖面上，礦體Fe1#、Fe2#、Fe3#采出后，崩落采空區(qū)圍巖形成的冒落帶高度均小于采高，但裂隙帶高度均大于采高，地表產(chǎn)生裂隙；在3-7線剖面上，F(xiàn)e1#及Fe3#采空區(qū)處理后形成的冒落帶和裂隙帶高度均小于采高，對尾礦庫穩(wěn)定性無影響，F(xiàn)e2#采空區(qū)處理后形成的冒落帶高度小于采高，裂隙帶高度超過采高，地表存在裂隙。綜合分析地表裂隙并不發(fā)育，對地表尾礦庫穩(wěn)定性影響不大，但隨時間推移裂隙將進(jìn)一步發(fā)育，可能對尾礦庫產(chǎn)生影響，應(yīng)加強(qiáng)對地表監(jiān)測，及時消除尾礦庫的安全隱患。

2.2巖石移動帶理論分析

采用工程類比法，選取上下盤的巖石移動角65°，崩落角70°，據(jù)此圈定3-3線及3-7勘探線剖面的巖石移動范圍和巖石崩落范圍，如圖3、圖4所示。

圖3　3-3線剖面移動帶

圖4　3-7線剖面移動帶

由圖3、圖4可知，按選定角度圈定的巖石移動范圍和崩落范圍，一部分尾礦庫位于圈定的范圍內(nèi)，表明尾礦庫存在安全隱患，但對尾礦庫的排水涵管無影響。由于該方法沒有考慮礦體埋深、礦體厚度等因素，因此具有一定的局限性。

2.3“三帶”和巖石移動帶理論綜合分析

“三帶”理論可以分析出采空區(qū)的變形破壞能否延展至地表，但對地表的影響范圍不能確定。巖石移動帶理論能夠彌補(bǔ)“三帶”理論的不足，便于確定地表影響范圍。本次研究基于“三帶”理論，采用巖石移動帶評價方法對采空區(qū)進(jìn)行綜合評價，最終得出采空區(qū)對地表的影響范圍，結(jié)果如圖5、圖6所示。

圖5　3-3線剖面移動帶

由圖5、圖6可知，3-3及3-7剖面有部分采空區(qū)的裂隙帶發(fā)育至地表，距尾礦庫最小距離分別為127.43 m和196.06 m，所有采空區(qū)冒落帶高度均小于采高，故采空區(qū)對尾礦庫穩(wěn)定性基本無影響。根據(jù)巖石移動理論計算分析處理后，采空區(qū)地表裂隙帶范圍如圖7所示。尾礦庫位于崩落法處理采空區(qū)地表裂隙帶范圍之外，尾礦庫的穩(wěn)定性能夠得到保障。

圖6　3-7線剖面移動帶

圖7　巖石移動界限

2.4數(shù)值模擬綜合分析

基于有限元ANSYS軟件，建立采空區(qū)數(shù)值模型，分別對采空區(qū)處理前后的地表水平及豎直位移進(jìn)行模擬研究，分析尾礦庫最大位移變化，最終圈定巖石移動對尾礦庫的影響范圍。

2.4.1模型建立及巖石力學(xué)參數(shù)

針對3-3及3-7兩個剖面，分別建立(長×高)566 m×500 m、626 m×504 m的兩個巖體數(shù)值模型，選取plane 82模型單元，模型左右兩側(cè)施加水平位移約束，底部施加垂直位移約束，地表為自由面。礦巖力學(xué)參數(shù)如表2。

表2　礦巖力學(xué)參數(shù)

2.4.2采空區(qū)未處理時數(shù)值模型

當(dāng)采空區(qū)未處理時，通過施加重力場，對兩剖面采空區(qū)進(jìn)行數(shù)值模擬，分析采空區(qū)對地表水平位移影響，如圖8～9所示。

圖8　3-3采空區(qū)未處理時位移云圖

圖9　3-7采空區(qū)未處理時位移云圖

依據(jù)采空區(qū)建筑地基穩(wěn)定性評價標(biāo)準(zhǔn)，水平方向位移變化量為0～14 mm、垂直方向位移變化量為0～40 mm時，采空區(qū)對地表無影響。由圖8～圖9可知，兩個剖面水平位移分別為2 mm和4 mm，未超過水平位移變化量；兩個剖面的豎直位移為分別52 mm和60 mm，均大于豎直位移變化量。因此，采空區(qū)未處理不利于地表尾礦庫穩(wěn)定。

2.4.3崩落處理采空區(qū)數(shù)值模型

崩落充填采空區(qū)后，通過數(shù)值模擬研究，地表水平及豎直位移較采空區(qū)未處理時變化較大，結(jié)果如圖10～11所示。

圖10　采空區(qū)處理后3-3位移云圖

由圖10～11可知，剖面3-3和3-7的水平位移分別為6 mm及5 mm，豎直位移分別為30 mm和20 mm，采空區(qū)地表影響范圍距尾礦庫最小距離分別為57 m和68 m。根據(jù)評價標(biāo)準(zhǔn)沉降量，剖面3-3和 3-7采空區(qū)對地表尾礦庫無影響，崩落處理采空區(qū)，有效的控制了空區(qū)對地表的影響范圍，確保了尾礦庫穩(wěn)定。

圖11　采空區(qū)處理后3-7位移云圖

3結(jié)論

針對該礦空區(qū)現(xiàn)狀及尾礦庫與采空區(qū)的空間位置關(guān)系，選定3-3和3-7兩個典型剖面，分別運用巖石移動帶理論、“三帶”理論、綜合兩種方法及數(shù)值模擬對采空區(qū)進(jìn)行分析。研究表明，采用崩落法處理采空區(qū)后，地表裂隙帶范圍減小，水平和數(shù)值位移沉降量降至安全范圍內(nèi)，對地表影響很小，采空區(qū)地表破壞點距離尾礦庫最小距離為127 m，確保了尾礦庫穩(wěn)定性。但隨著深部礦體的進(jìn)一步開采，裂隙還會逐漸發(fā)展，應(yīng)加強(qiáng)監(jiān)測，確保安全生產(chǎn)。

參考文獻(xiàn)

[1]孫超,薄景山,孫有為.采空區(qū)沉陷研究歷史及現(xiàn)狀[J].防災(zāi)科技學(xué)院學(xué)報,2008(4):128-131.

[2]梁力,李明,王偉,等.尾礦庫壩體穩(wěn)定性數(shù)值分析方法[J].中國安全生產(chǎn)科學(xué)技術(shù),2007,3(5):11-15.

[3]李宏儒,胡再強(qiáng)，陳存禮，等.金堆城尾礦壩加高方案數(shù)值模擬及穩(wěn)定性分析[J].巖土力學(xué),2008，29(4):1138-1142.

[4]李國政,李培良,徐宏達(dá).基于結(jié)構(gòu)可靠度指標(biāo)的尾礦庫壩體穩(wěn)定性分析[J].黃金，2005，26(6)：48-50.

(收稿日期2015-10-27)

Study of the Influence of Using Caving Rock to Deal with Goaf to Tailings

Liu Wangzheng1Zhang Yunpeng1Liu Aixing1,2Zheng Wenlong3

(1.School of Mining Engineering,North China University of Science and Technology;2.Xuzhou Iron Group Co.,Ltd.; 3.Sijiaying Iron Mine, Heibei Iron & Steel Group Co.,Ltd.)

AbstractPart of the tailings is existed at the surface of a underground iron mine goaf, in order to ensure the tailings stability, the caving filling method is adopted to deal with the goaf, the “three zones”theory, rock displacement belt theory and numerical simulation method are adopted to analyze the tailings stability. The results show that, after dealing with the goaf by caving filling method, the surface horizontal displacement variation and vertical displacement variation are within the standard range, the minimum distance of surface subsidence fissure belt to tailings is 127 m, the goaf treatment could ensure the surface tailings stability, with the further mining of the deep ore-body of the mine, it is necessary to strengthen the monitoring of surface crack zone so as to ensure the safety of surface tailings.

KeywordsCaving rock, Goaf, Tailings, Stability analysis

劉王政(1988—)，男，碩士研究生，063009 河北省唐山市新華西道46號。