紅嶺鉛鋅礦采場結構參數優化研究

張建偉 趙江波 姚景亮

收稿日期：2023-12-06； 修回日期：2024-02-18

基金項目：國家自然科學基金面上項目（51974043）

作者簡介：張建偉（1980—），男，工程師，從事金屬礦充填開采相關研究工作；E-mail：1009589285@qq.com

摘要：為確定紅嶺鉛鋅礦分段鑿巖階段空場嗣后充填采礦法的合理采場結構參數，采用數值模擬的方式對頂板厚度分別為3 m、5 m、7 m和10 m，采場長度分別為20 m、30 m、40 m和60 m的16種工況進行了研究，對比了頂板厚度為3 m情況下各工況的應力、位移和塑性區分布。結果表明：紅嶺鉛鋅礦分段鑿巖階段空場嗣后充填采礦法采場頂板厚度3 m，采場長度35 m時，采場穩定性較好。

關鍵詞：分段鑿巖；階段空場嗣后充填；采場結構參數；頂板厚度；采場長度；采區穩定性

中圖分類號：TD853.34????????? 文章編號：1001-1277（2024）06-0006-05

文獻標志碼：Adoi：10.11792/hj20240602

引? 言

隨著人類社會不斷發展，人們對礦產資源的需求不斷增加，采礦活動逐漸由淺部走向深部［1-3］；在建設綠色礦山、推動綠色發展的背景下，充填采礦法以充填體支撐頂板、回收利用礦山固體廢物等優點逐漸成為硬巖礦山地下開采的首選方式［3-6］。

赤峰山金紅嶺有色礦業有限責任公司（下稱“紅嶺鉛鋅礦”）位于內蒙古自治區巴林左旗烏蘭達壩蘇木。紅嶺鉛鋅礦床是以鉛、鋅為主的矽卡巖型多金屬礦床，產于燕山晚期花崗巖與二迭系大石寨組地層的外接觸帶上，出露地層有上古生界二迭系大石寨組、中生界侏羅系滿克頭顎博組及新生界第四系松散沉積物。該區域內構造交錯復雜，圍巖較為破碎，955～855 m采場留有大量礦柱，下部中段礦體沒有開采，導致大量礦體積壓，且該區域內礦體距離地表高差僅有150 m左右，礦體埋深較淺，礦體開采易造成上盤圍巖不穩固。

對于紅嶺鉛鋅礦南西翼礦區955～855 m采場，擬采用分段鑿巖階段空場嗣后充填采礦法進行開采，在前人研究的基礎上［7-11］，采用有限元法分析不同工況下采場及圍巖的穩定性，結合數值模擬的方式，應用基于樣本點構造方法的先進結構優化方法對采場主要參數進行選擇，在確保回采安全、保證最大生產效率，為礦山創造最大經濟效益的前提下，確定最優的頂板厚度和采場結構參數。

ABAQUS作為通用的模擬工具，除了能解決大量結構（應力/位移）問題，還可以在非線性分析中，自動選擇相應載荷增量和收斂限度，通過準確的定義參數能夠很好地控制數值計算結果，保證分析過程中能夠得到有效的精確解，為紅嶺鉛鋅礦分段鑿巖階段空場嗣后充填采礦法采場結構參數研究提供了強有力的計算工具。

1? 數值模擬材料及參數

1.1? 計算模型

紅嶺鉛鋅礦分段鑿巖階段空場嗣后充填采礦法標準方案如圖1所示。分礦房、礦柱進行回采，開采時先采礦柱再采礦房。根據紅嶺鉛鋅礦現場情況，分別對頂板厚度為3 m、5 m、7 m和10 m，采場長度為20 m、30 m、40 m和60 m的16種工況進行模擬，計算工況如表1所示，其三維模型如圖2所示。

1.2? 材料模型及參數

根據現場巖體力學調查及充填報告中實驗室測試結果，折減后圍巖、礦體及廢石充填體的物理力學參數如表2所示。

2? 開采模擬及結果

2.1? 采場圍巖應力

從已有的力學參數及巖體質量調查來看，上盤圍巖力學性能低于下盤圍巖，且礦體呈70°左右傾斜，因此礦體上盤的圍巖更易表現出不穩定。由于巖石為脆性材料，一般采用第一主應力破壞強度理論，通過上盤的最大主應力圖可以看出拉應力比較大的區域，在主應力最大值附近設置觀察點，對比觀察點處的拉應力值與巖體抗拉強度分析可能出現拉應力破壞的區域。

本文針對頂板厚度為3 m時的采場（工況1、工況5、工況9、工況13）進行分析，通過有限元分析法來研究采場及上盤圍巖的穩定性，不同采場長度情況下的最大主應力剖面圖如圖3所示。

2024年第6期/第45卷? 礦業工程礦業工程? 黃? 金

分析不同情況下最大主應力剖面圖：在工況1中，采場及其上盤圍巖最大主拉應力超過巖體抗拉強度的區域在垂直方向上呈點狀分布；在工況5中，采場及其上盤圍巖最大主拉應力大于巖體抗拉強度的區域在水平方向和垂直方向上均存在，分別呈線性分布；在工況9中，采場上盤圍巖最大主拉應力超過巖體抗拉強度的區域幾乎貫通，且呈片狀分布；在工況13中，采場上盤圍巖最大主拉應力超過巖體抗拉強度的區域已經貫通，且呈扇形分布。在輸出結果中繼續選取如圖4所示的觀察路徑。

對比觀察路徑上各節點的應力狀況，查看當頂板厚度為3 m時，在不同采場長度情況下各點處的拉應力值，將結果與巖體的抗拉強度對比。觀察路徑1，2上相應觀察點最大主應力與觀察點沿路徑距離的關系曲線如圖5、圖6所示。

工況1中：垂直采場長度方向上，采場及其上盤圍巖有26 m范圍內的巖體處于拉應力區域內。其中，距采場底面9.67～12.79 m的巖體拉應力大于巖體的抗拉強度，該區域內巖體受拉破壞。沿采場長度方向上，分別在沿采場長度1/3和2/3處有最大拉應力，值為1.16 MPa。

工況5中：垂直采場長度方向上，采場及其上盤圍巖有24 m范圍內的巖體處于拉應力區域內。其中，距采場底面9.66～12.98 m的巖體拉應力大于巖體的抗拉強度，該區域內巖體受拉破壞。沿采場長度方向上，分別在沿采場長度7/25和18/25處有最大拉應力，值為0.87 MPa。

工況9中：垂直采場長度方向上，采場及其上盤圍巖有23 m范圍內的巖體處于拉應力區域內。其中，距采場底面10.41～13.74 m的巖體拉應力大于巖體的抗拉強度，該區域內巖體受拉破壞。沿采場長度方向上，分別在沿采場長度1/5和4/5處有最大拉應力，值為0.71 MPa。

工況13中：垂直采場長度方向上，采場及其上盤圍巖有25.5 m范圍內的巖體處于拉應力區域內。其中，距采場底面7.30～10.63 m的巖體拉應力大于巖體的抗拉強度，該區域內巖體受拉破壞。沿采場長度方向上，分別在沿采場長度3/40和37/40處有最大拉應力，值為1.12 MPa。采場邊緣位置豎直方向的應力狀態如表3所示。

從以上幾點可以看出，不同長度的采場及其上盤圍巖都會在距采場底面10 m左右產生受拉破壞，因此在開采該部分區域時應該注意加強支護，保證安全。

2.2? 采場圍巖位移

數值計算中，將不產生塑性貫通區的頂板厚度定為各采場的安全厚度，不同頂板厚度下采場上盤圍巖最大位移如圖7所示。由圖7可知：隨采場長度增大，圍巖的位移逐漸增大，采場長度為20～40 m時，其上盤圍巖最大位移相差不到10 mm，而采場長度為60 m時，其上盤圍巖最大位移可達52.74 mm；不同頂板厚度對于采場上盤圍巖的位移幾乎沒有影響。因此，考慮到礦山經濟效益，建議紅嶺鉛鋅礦南西翼礦區頂板厚度為3 m。

2.3? 采場圍巖塑性區變化

ABAQUS有限元分析時，判斷產生破壞的標準之一是產生塑性區的貫通。沿采場長度方向，拉應力狀態區的應力值都接近或大于巖體的抗拉強度，對比采場上盤圍巖的塑性區發展圖（如圖8所示），將采場上盤圍巖從中間剖開，觀察圍巖內部的塑性區發展。

工況1中，采場及其上盤圍巖屈服區呈點狀及線狀分布，但上下兩屈服區未形成貫通；工況2中，采場及其上盤圍巖屈服區呈線狀分布；工況3中，采場及其上盤圍巖在垂直方向上屈服區仍呈線狀分布，且面積進一步擴大；工況4中，采場及其上盤圍巖屈服區相連成片，這些區域表明巖石由于本身的性質或外部因素的影響而產生了破壞，該現象表明，該區域存在著一定的不穩定性。同時，可以觀察到，在不同采場長度下，其上盤圍巖的塑性貫通區與底部可能產生拉應力破壞的區域完全一致。因此，在頂板厚度3 m的情況下進行開采時，采場底部區域也會產生沿采場長度方向的破壞。

分析不同工況的頂板塑性區發展情況，結果如表4所示。

其中，采場長度為30 m及40 m時，破壞向內部延伸的最小，且二者破壞程度相近。因此，選取30 m或40 m的采場長度是合理的，考慮到開采效率，建議紅嶺鉛鋅礦南西翼礦區采場長度為35 m。

3? 結? 論

1） 采場長度為20～60 m時，采場及上盤圍巖最大主拉應力值接近或超過巖體抗拉應力的區域均位于距采場底面10 m左右的區域內。

2） 隨采場長度增大，圍巖的位移逐漸增大，采場長度為20～40 m時，采場及其上盤圍巖最大位移相差不到10 mm，不同頂板厚度對于采場上盤圍巖的位移幾乎沒有影響。

3） 采場長度為30 m及40 m時，破壞向內部延伸的最小，且二者破壞程度相近。

綜上所述，紅嶺鉛鋅礦南西翼礦區采用分段鑿巖階段空場嗣后充填采礦法的采場建議頂板厚度為3 m，采場長度為35 m，同時在開采距采場底面10 m左右的區域時應加強支護，保證作業安全。

［參 考 文 獻］

［1］? PENG K，LIU Z M，ZHANG Y L，et al.Determination of isolation layer thickness for undersea mine based on differential cubature solution to irregular Mindlin plate［J］.Journal of Central South University，2017，24（3）：708-719.

［2］? 彭康，李夕兵，彭述權，等.基于響應面法的海下框架式采場結構優化選擇［J］.中南大學學報（自然科學版），2011，42（8）：2 417-2 422.

［3］? 古德生，李夕兵.有色金屬深井采礦研究現狀與科學前沿［C］∥中國有色金屬學會.中國有色金屬學會第5屆學術年會論文集.北京：中國有色金屬學會，2003：1-5.

［4］? 周逸文，張濤，段隆臣，等.我國礦山采空區綜合治理研究綜述［J］.安全與環境工程，2022，29（4）：220-230.

［5］? 侯俊.某礦山浮選尾礦充填材料特性試驗研究［J］.黃金，2023，44（8）：26-29，34.

［6］? 梅世洲.充填采礦工藝特點及發展趨勢分析［J］.現代礦業，2016，32（3）：253-254.

［7］? 余昕，周家祥，宋衛東，等.分段鑿巖階段空場嗣后充填采礦工藝及爆破設計［J］.金屬礦山，2018（8）：7-12.

［8］? 汪海萍，宋衛東，張興才，等.大冶鐵礦分段鑿巖階段空場嗣后膠結充填法工程實踐［J］.礦業研究與開發，2015，35（5）：26-29.

［9］? 孟慶學，宋德林，李朝民.中深孔階段礦房采礦法在急傾斜薄礦體中的應用［J］.黃金，2022，43（3）：36-39.

［10］? 楊春城.某礦分段采場上盤圍巖穩定性分析及工程應用［J］.金屬礦山，2019（4）：20-25.

［11］? 楊春城.某礦分段充填法采場結構參數優化及工程應用［J］.礦業研究與開發，2019，39（4）：1-5.

Study on the optimization of stope structural parameters of Hongling Lead-Zinc Mine

Zhang Jianwei，Zhao Jiangbo，Yao Jingliang

（Chifeng Shanjin Hongling Nonferrous Mining Co.，Ltd.）

Abstract：To determine the reasonable stope structural parameters for the sublevel drilling stage open stope and subsequent filling mining method in Hongling Lead-Zinc Mine，numerical simulation was employed to study 16 different scenarios with roof thicknesses of 3 m，5 m，7 m，and 10 m，and stope lengths of 20 m，30 m，40 m，and 60 m，respectively.The stress，displacement，and plastic zone distribution of each scenario with a roof thickness of 3 m were compared.The results indicate that the stability of the stope in Hongling Lead-Zinc Mine is fairly good when the roof thickness is 3 m and the stope length is 35 m for the sublevel drilling stage open stope and subsequent filling mining method.

Keywords：sublevel drilling;stage open stope and subsequent filling mining;stope structural parameters;roof thickness;stope length;mining area stability