大紅山銅鐵礦地壓活動規律數值模擬分析

朱國輝,胡平安

(1.長沙有色冶金設計研究院,湖南長沙 410011;2.湖南有色金屬研究院,湖南長沙 410015)

·采 選·

大紅山銅鐵礦地壓活動規律數值模擬分析

朱國輝1,胡平安2

(1.長沙有色冶金設計研究院,湖南長沙 410011;2.湖南有色金屬研究院,湖南長沙 410015)

頻繁的地壓活動對大紅山銅鐵礦安全生產構成了嚴重威脅。通過礦區地壓演化過程的數值模擬分析,掌握了礦區地壓活動產生的原因及規律,根據模擬結果提出了相應的治理措施。研究結果對控制礦區地壓及安全生產具指導作用。

地壓活動;數值模擬;塑性區;拉應力區

受工程地質條件復雜及多中段高強度開采影響,大冶大紅山銅鐵礦區發生了較明顯的地壓活動,主要表現在:采場頂、底板圍巖發生大面積跨塌;部分礦(間)柱產生張裂縫及滑移,采場內大量充填料沿張裂縫泄漏;穿脈巷道出現冒頂、片幫及底鼓,地壓活動出現從局部向全礦區發展的趨勢,對礦山安全生產構成了嚴重威脅。

為了解地壓活動產生原因,準確預測礦區地壓活動發展趨勢,開展了礦區地壓活動規律數值模擬分析,結論與建議對礦區地壓控制及安全生產具一定指導作用。

1 礦區概況及開采技術條件

1.1 礦區概況

礦區三面環湖,南面居民房屋、學校、工礦企業密集,武-九鐵路、大-金省道緊繞礦區而過,地表不允許塌陷。

礦區-200 m以上民采空區積水嚴重,現已停止采礦活動,-200～-220 m標高間設20 m厚防水隔離礦柱,-220 m以下屬大冶大紅山銅鐵礦開采范圍,共設-270 m、-320 m、-370 m、-420 m 4個中段。兩步驟回采,一步驟礦房采用淺孔留礦法開采,嗣后分級尾砂膠結充填;二步驟礦(間)柱采用上向低分段空場法開采,嗣后分級尾砂充填。為滿足開采規模要求,礦區實行多中段生產。因某方面原因,已造成-320 m、-370 m中段部分礦房、礦柱完全錯位。

1.2 開采技術條件

礦床屬中酸性花崗閃長斑巖與嘉陵江組灰巖接觸交代作用形成的高中溫熱液矽卡巖型銅鐵礦床。礦體受舌狀體或捕虜體和反“S”型接觸帶構造控制。主礦體走向長900 m,厚4～100 m,傾角5°～85°,傾向延伸 50～550 m,平均品位 Cu:1.37%,TFe: 37.26%,Au:0.46 g/t,Ag:7.07 g/t。

礦體頂、底盤圍巖為花崗閃長巖、大理巖及矽卡巖。花崗閃長巖堅硬穩固,大理巖強度高,但層理面發育,矽卡巖與礦巖接觸帶內巖體松軟、破碎。其中礦體與下盤花崗閃長巖接觸帶一旦揭露遇水即發生膨脹崩解,強度迅速降低。礦體與上盤大理巖接觸帶部位巖體強度相對較低,但接觸帶影響寬度較小,且結合較緊密。

礦區水文地質條件復雜。直接頂板大理巖溶發育。強巖溶發育帶底部平均標高-72.47 m,弱巖溶發育帶底部平均標高-326.84 m,弱帶以下巖溶基本消失。礦區實施帷幕注漿防治水工程后,幕內外水位差達36.48～36.91 m,堵水率60%左右。

2 數值分析模型的創建

2.1 模型的建立

數值模擬的可靠性一定程度上取決于所選取的計算模型。因礦區工程地質條件復雜,三維建模繁雜,且根據工程實踐,同等條件下,二維數值模擬結果較三維數值模擬結果在應力和位移方面高約10%,此對工程巖體穩定性評價是有利的,故本次分析采用二維模型。

數值模擬采用有限元分析軟件。礦房垂直投影尺寸12.5 m×46 m,高50 m。為更真實地反映礦區地壓演變過程,避免應力過度集中,礦房頂、底部邊角進行了圓角處理。

需要指出的是:數值模擬無論采用二維還是三維模型,均經過了較多簡化,與實際情況存在一定差異,計算結果雖不能完全準確地反應客觀實際,但用于宏觀定性分析有較好的參考作用。就本研究而言,數值模擬結果在一定程度上反映了采礦活動對礦巖及充填體穩定性的影響,初步揭示了礦區地壓活動規律特征,對礦區合理回采順序的選擇具一定指導作用。

2.2 材料參數

數值計算結果的可靠性還取決于巖體力學參數選取的合理性,受客觀條件及技術難度影響,進行現場大型原位巖體力學參數試驗十分困難。一般工程則以實驗室巖塊試驗為基礎,同時考慮巖體的結構效應,結合工程實踐對巖塊力學參數進行適當修正得到數值計算所用材料參數,詳見表1。

表1 礦巖物理力學參數表

2.3 原巖應力場

根據中科院武漢巖土力學研究所提供的礦區原巖應力測試結果表明:礦區最大主應力為水平方向的構造應力,其次為垂直方向的自重應力,二者比值0.46～1.36。本研究計算中,垂直方向應力按自重應力計算。其中-320 m水平垂直方向自重應力約為9.50 MPa,水平方向的應力約為13 MPa。

2.4 礦巖破壞準則

圍巖破壞按兩種情況考慮:

1.張拉破壞。處于空區頂部巖體,受自重作用而產生受拉破壞。受拉應力作用巖體破壞方式按最大拉應力準則處理,即:

2.壓剪破壞。開挖作用影響,空區周邊會產生壓應力集中。處于壓應力集中區的巖體,可能產生壓剪破壞。壓剪破壞采用德魯克-普拉格(Drucker -Prager)準則。即:

式中I1為第一應力不變量;J2為第二偏應力不變量;α、β為正常數,它可由巖塊的c,φ求得:

3 模擬結果與分析

3.1 -270 m中段礦房開采時

1.圍巖中僅局部產生較大位移。最大位移產生在端部礦(間)柱中央,方向近似水平,指向采空區。

2.礦房開采后,礦(間)柱中水平方向應力得到解除,垂直方向應力變為最大主應力。在端部礦房中,出現較大范圍的垂直方向的拉應力,中間礦(間)柱中,出現量值較小的水平方向的拉應力(見圖1)。

3.礦房頂、底部出現局部塑性區,范圍較小。

圖1 -270 m中段礦(間)柱中拉應力情況

3.2 充填-270 m中段礦房、開采-320 m中段礦房時

1.-270 m中段礦房采空區充填后,該中段礦(間)柱垂直方向應力值有所降低;受充填體作用,礦(間)柱中的水平拉應力消失。

2.-320 m中段礦房開采所形成的拉應力區形狀與上中段礦房開采時所形成的形狀基本一致,主要發生在礦(間)柱中。端部礦房拉應力方向垂直,量值最大。中間礦房拉應力方向水平,量值較小(見圖2)。

圖2 -320中段開采時礦(間)柱中拉應力情況

3.礦房頂部出現較小范圍的塑性區。

4.充填體中的應力為壓應力,但量值較小。

3.3 不充填-320 m中段礦房直接開采-370 m礦房時

1.由于上下礦房、礦(間)柱不對應,開采-370 m水平礦房之后,-320 m礦房與-370 m礦房之間形成較大范圍塑性區,塑性區貫通了整個頂底柱(見圖3)。

圖3 -320 m中段貫通的塑性區

2.除上部礦房開采時在礦(間)柱中產生拉應力區之外,-370 m中段礦房頂部也出現了拉應力區,而且拉應力量值明顯增大。

3.4 充填所有礦房后開采-270 m中段礦(間)柱時

1.當開采礦(間)柱時,主要荷載還是由遠處其余礦(間)柱承擔,充填體只承擔了較小一部分荷載(見圖4)。

圖4 礦房充填體中應力分布

2.從計算結果來看,在礦(間)柱開采時,充填體產生了較大范圍的塑性區,該塑性區水平方向上完全貫通了整個充填體(見圖5)。

3.礦(間)柱開采時,其空區頂部產生了較大范圍的拉應力區,可能影響-200～-220 m間隔離礦柱的穩定(見圖6)。

3.5 所有礦房充填后開采-320 m中段礦(間)柱時

1.開采-320 m中段礦(間)柱時塑性區范圍及應力大小與開采-270 m中段礦(間)柱情形基本一致。

圖5 充填體在水平方向產生貫通性塑性區

圖6 空區頂部(隔離礦柱)出現拉應力區

2.在空區頂部,出現更大范圍的拉應力區。對應的頂部礦(間)柱中垂直方向的壓應力得到了部分釋放(見圖7)。

圖7 空區頂部出現的拉應力區

4 結論與建議

1.一步驟礦房開采建議按從上至下單中段或相鄰上、下兩個中段同時生產;合理布置采場,確保礦房、礦柱上、下對應;堅持“強采強出強充”原則,盡量縮短采場回采周期;保證礦房充填體強度,確保二步驟礦(間)柱回采安全。

2.礦房開采過程中,礦柱經歷了“由垂直方向出現拉應力到水平方向應力解除再到垂直方向應力增大”的應力演化過程,該過程使礦柱中的結構面經歷了“垂直方向張開至垂直方向閉合進而沿結構面滑移”的一個過程,這是礦柱中部急傾斜結構面產生破壞的根本原因。

3.礦柱中部的主應力方向產生了互換。礦柱兩側礦房開采后,垂直方向壓應力增大,水平方向壓應力基本為零,此時礦柱中結構面很易產生受壓破壞。

4.基于結論與建議2、3方面原因,對長大礦柱的回采,可考慮在其長軸方向上將礦房分解為2～3個采場,以改善礦柱的受力狀況,減少礦房開采對礦柱的破壞。

5.上、下礦柱不對應造成上、下相鄰礦房之間礦柱產生貫通性塑性區以及對應礦柱垂直方向應力部分解除是大紅山礦區產生威脅性地壓活動的兩個重要誘導因素。

6.礦房充填體對改善礦柱的應力分布起到了重要作用。礦房開采后,宜盡快充填,盡量地降低暴露時間;充填體要有足夠強度,應能完全承擔荷載。

7.-200～-220 m防水隔離礦柱下礦柱回采會在隔離礦柱中形成拉應力區,對隔離礦柱的穩定產生不利影響。因此,在礦區上部民采空區積水威脅解除之前,建議各中段礦柱的回采順序采用先采下部礦柱后采上部礦柱的上行式順序,且礦柱回采時,應加強兩側礦房充填體的保護及全礦區地壓監測。

[1] 蔡美峰.金屬礦山采礦設計優化與地壓控制—惠顧論與實踐[M].北京:科學出版社,2001.

[2] 胡飛宇.安慶銅礦特大型采場充填體穩定性數值模擬研究[J].有色金屬(礦山部分),2006,58(3):14-17.

[3] 李愛兵,周先明.安慶銅礦高階段回采充填體-礦體-巖體穩定性的有限元分析[J].礦業研究與開發,2000,20(1):19-21.

[4] 謝茂輝,伍佑倫,馮聲浪,等.留礦法嗣后充填開采中的地壓活動與控制[J].礦業研究與開發,2006,26(B11):16-29.

[5] 楊彪,羅周全,劉曉明,等.凡口鉛鋅礦礦柱回采穩定性數值模擬研究[J].礦業研究與開發,2007,27(5):26-28.

[6] 秦艷華,王曉軍,鐘春暉,等.礦柱回采對采空區穩定性影響數值模擬研究[J].銅業工程,2007,3:8-10.

Numerical Simulation Analysis of the Rule of Geostatic Activity in Dahongshan Copper-iron Mine

ZHU Guo-hui1,HU Ping-an2

(1.Changsha Engineering and Research Institute of Nonferrous Metallurgy,Changsha410011,China;2.Hunan Research Institute of Nonferrous Metals,Changsha410015,China)

Frequent geostatic activity has been a serious threat to production safety in the Dahongshan copper-iron mine.Through numerical simulation analysis of the evolution process of geostatic activity in the mining area,the reason and rule of the geostatic activity has been taken,and corresponding control measures has been brought up according to the simulation results.Research results have a guiding role to control geostatic activity and safety production.

geostatic activity;numerical simulation;plastic zone;tensile stress area

TD167

A

1003-5540(2011)05-0001-04

朱國輝(1965-),男,高級工程師,主要從事采礦工程設計工作。

2011-08-20