基于經濟指標的下向進路斷面安全系數法分析

紀永剛 彭 康 趙 杰 李紀玉(.錫林郭勒盟山金阿爾哈達礦業有限公司,內蒙古 東烏珠穆沁旗 026300；2.重慶大學煤礦災害動力學與控制國家重點實驗室,重慶 400044； 3.重慶大學資源與環境科學學院，重慶 400044)

基于經濟指標的下向進路斷面安全系數法分析

紀永剛1彭 康2,3趙 杰1李紀玉1
(1.錫林郭勒盟山金阿爾哈達礦業有限公司,內蒙古 東烏珠穆沁旗 026300；2.重慶大學煤礦災害動力學與控制國家重點實驗室,重慶 400044； 3.重慶大學資源與環境科學學院，重慶 400044)

為實現下向進路安全高效經濟開采，采用彈性力學理論建立了下向進路薄板力學模型，得到充填體失穩破壞抗拉強度判據，在此基礎上建立了下向進路P(h/M)經濟模型。以山東某金礦為例分析了安全系數為1.8時，不同進路寬度、高度與下向進路充填經濟指標的關系，研究表明：經濟指標與進路寬度呈直線增長關系，與進路高度呈冪函數降低關系，且M≤1.5 m時經濟指標大幅下降，M>1.5 m時經濟指標緩慢下降。進而研究了進路半寬為1.3 m到1.7 m，高度為2.8 m到3.6 m的經濟指標，建議將傳統安全系數分析結果(l=1.5 m、M=3.6 m、h=1.49 m)改為經濟指標分析結果(l=1.4 m、M=3.6 m、h=1.40 m)，此時充填成本降低1.64元/m3。由此得出，該方法具有很好的理論及現實意義。

薄板模型 下向進路 安全系數 經濟指標

下向進路膠結充填采礦法是回采高品位、貴重及稀有金屬、礦巖極破碎礦床的主要方法[1]，隨著無軌自行設備在地下金屬礦山的推廣應用，下向進路膠結充填采礦法以其回收率高、貧化損失小、回采強度大等顯著優勢，進入現代化采礦技術之列[2-3]，而安全高效經濟開采一直是礦山追求的目標。目前針對下向進路的穩定性研究主要有力學分析[4]、數值模擬[5]、安全系數法[6]和可靠度分析[7]，對下向進路經濟研究主要從尋求廉價的充填骨料和水泥代替材料、優化充填料漿配合比、改變充填方式與工藝等方面進行[8-10]，未見下向進路穩定條件下經濟最優分析方法。

本研究在建立下向進路薄板力學模型的基礎上，提出了下向進路經濟評價模型，結合山東某金礦對下向進路寬度、高度對經濟指標的影響進行研究，并設計不同進路斷面得到其經濟指標，進而對基于經濟指標安全系數法與傳統安全系數法結果比較，提出最優進路斷面，降低充填成本。

1 下向進路薄板力學模型

根據進路開采中承載層受力特征，作如下假設：

(1)進路承載層為連續、均質、各向同性的，符合彈性力學假設條件的彈性板；

(2)礦體和承載層在屈服破壞之前為線彈性體，其本構方程為σ=Eε；

(3)承載層厚度h與承載層水平方向上的最小尺寸L的比值：h/L≤1/5；

(4)承載層上受均布荷載q；

(5)開采中進路長度都遠大于進路寬度。

由此可將進路側幫視為彈性基礎，承載層視為彈性基礎之上由彈性介質組成的薄板來研究[11-12](如圖1)，并由文獻[13]得到：

圖1 下向進路薄板力學模型及受力分析

(1)當-L≤x≤l上，xc=0處彎矩有極值(不考慮彎矩作用方向)：

(3)

其中，

式中，Ej為進路側幫基礎的彈性模量；EL為頂板的彈性模量；μ為充填體的泊松比；h為頂板厚度，一般h≤1.5 m；l為進路寬度之半，一般l≥1.25 m；M為進路高度。

(2)當x>l時(x<-l與其對稱，不再重述)，

處，彎矩有極大值：

(4)

由此得到下向進路最大拉應力計算式(5)，并以此作為判斷充填體穩定性的判據：若承載層下表面拉應力超過其抗拉強度，則認為充填體發生失穩破壞[8-9]。

σtmax=6Mmax/h2.

(5)

2 下向進路經濟評價模型

下向進路充填體由承載層、補口層和接頂層構成，設承載層充填成本a元/m3，補口層和接頂層充填成本b元/m3，由此得到單位體積充填料漿成本P的計算為

(6)

承載層、補口層和接頂層單位體積充填成本通常視為定值，且a>b>0，可知充填成本P為h/M的函數，且P(h/M)隨h/M的減小而減小，故h/M取值越小越經濟。

由此得到下向進路參數滿足充填體抗拉強度判據式(1)，且式(6)中h/M越小，則方案越優，從而實現礦山安全經濟開采。

3 工程實例

3.1 工程背景

山東某金礦屬急傾斜中厚極破碎礦體，礦山采用C料尾砂下向進路充填采礦法，進路寬度L=3.0 m，進路高度M=3.0 m，承載層厚度h=1.0 m。進路沿礦體走向時由上盤向下盤連續回采，垂直礦體走向時隔一采一，因此礦體回采中“硬支薄板”、“軟硬支混合”、“軟支薄板”結構均有，由文獻[8]知“軟支薄板”結構最易破壞，故可只對“軟支薄板”結構進行分析。由C料尾砂配比實驗得到承載層、補口層和接頂層力學參數如表1所示。

表1 充填體力學參數Table 1 Mechanical parameters of backfill

3.2 傳統安全系數法分析

下向進路承載層的穩定是進路回采的關鍵，因此承載層所受拉應力需留有一定的安全儲備，為此引入下向進路穩定性的安全系數判別式：

η=σtmax/σt,

(7)

式中,η為安全系數；σtmax為承載層所能承受最大拉應力；σt為承載層內實際最大拉應力。

根據地下礦山安全開采經驗，安全系數應在1.6～1.8以上[14]，故可取η≥1.8作為穩定條件評價標準。由現場觀察知承載層的穩定性較好，普通充填體穩定性較差，故每個分層所受荷載為承載層和普通充填體的自重，即：

q=γ1h+γ2(M-h).

(8)

令

η=σtmax/σt=1.8

其中，σt=0.42MPa，代入人工假頂力學參數，運用Mathematica編程得到進路高度M、半寬l與所需承載層厚度h關系(如圖2)，可知所需承載層厚度隨進路寬度和高度的增加而增加。為進一步研究進路寬度和高度對承載承載層厚度的影響，采用控制變量法得到圖3和圖4。

圖2 進路高度及半寬與所需承載層厚度關系

圖3 M=3.0 m時不同進路寬度下承載層厚度與安全系數關系

該金礦現有進路l=1.5 m、M=3.0 m，對承載層厚度與安全系數作非線性擬合，得到擬合關系為

圖4 l=1.5 m時不同進路高度下承載層厚度與安全系數關系

(9)

由此式易得現有承載層厚度h=1.0 m，安全系數η僅為1.27，遠低于要求值1.8，需將承載層厚度提高至1.27 m(η=1.8)，才可保證礦山安全開采。

進一步分析知進路寬度對下向進路的穩定性較大，而進路高度對下向進路的穩定性影響很小。建議保持進路寬度不變，而將進路高度調整為3.6 m，再擬合計算得到承載層厚度為1.49 m。

3.3 基于經濟指標的安全系數法分析

由上述分析知傳統安全系數法常通過確定的斷面計算得到滿足安全系數下的最低承載層厚度，或者先控制某個變量，再確定研究另外2個變量最優關系，而此時得到的進路參數并不一定是最經濟的，為此引入經濟指標進行分析。

3.3.1 進路斷面對經濟指標的影響

令η=1.8，運用 Mathematica 編程得到進路高度M、半寬l與經濟評價指標P(h/M)關系圖(如圖5)，由圖5可知經濟指標隨進路高度增加而減小，隨進路寬度增加而增加。

圖5 進路高度及半寬與經濟評價指標關系

為進一步得到進路高度、寬度對經濟指標的影響大小，采用控制變量法對其關系進行擬合。

(1) 當M=3.0 m時，經濟指標P(h/M)：

P(h/M)=0.261l+0.032.

(10)

可知經濟指標與進路寬度呈直線增長關系(如圖6所示)，且由系數0.261知，當進路寬度增加時，經濟指標增長較快，因此在滿足生產的條件下進路寬度需盡可能的小為好。

圖6 M=3.0 m時進路半寬與經濟指標關系曲線

(2) 當l=1.5 m時，經濟指標P(h/M)：

P(h/M)=0.574 5M0.294 1.

(11)

可知經濟指標與進路高度近似呈冪函數關系(如圖7所示)，經濟指標隨進路高度增加而減小，且M≤1.5 m時經濟指標大幅下降，M>1.5 m時經濟指標下降緩慢，因此進路高度應在大于1.5 m范圍內，且在滿足限制條件下需盡可能的大為好。

圖7 l=1.5 m時進路高度與經濟指標關系曲線

3.3.2 進路斷面優選方案設計

根據該金礦已有生產經驗，選取進路半寬從1.3 m到1.7 m，進路高度M從2.8 m到3.6 m，共25種不同的斷面組合，研究其在安全系數為1.8時的經濟指標見表2。

由表2，可得原有進路斷面l=1.5 m、M=3.0 m對應的經濟指標P(h/M)為0.425，傳統安全系數法分析結果l=1.5 m、M=3.6 m,對應的經濟指標P(h/M)為0.413。根據采礦工藝可將進路半寬調為l=1.4 m、M=3.6 m，此時對應的經濟指標為0.389。

表2 不同進路斷面下的經濟指標Table 2 The economic index at different drift sections

由C料尾砂膠結充填工藝知，承載層充填成本a=125.8元/m3，補口層和接頂層充填成本b=57.5元/m3，因此基于經濟指標分析結果較傳統安全系數法分析結果充填成本降低1.64元/m3，按該金礦現有生產能力23萬t/a(即8.58萬m3/a)，年可節約充填成本14.07萬元。

4 結 論

(1)建立了下向進路薄板力學模型和經濟評價模型，得到承載層在滿足最大拉應力判據時，經濟指標P(h/M)越小，方案越優，從而實現礦山安全經濟開采。

(2)傳統安全系數法分析知現有進路l=1.5 m，M=3.0 m，h=1.0 m，安全系數η僅為1.27，需提高至h=1.27 m(η=1.8)，才可保證礦山安全開采。進路寬度對穩定性影響較大，而高度影響很小，建議將進路參數調整為l=1.5 m、M=3.6 m、h=1.49 m。

(3)基于經濟指標安全系數法分析得到經濟指標與進路寬度呈直線增長關系，與進路高度呈冪函數關系，且M≤1.5 m時經濟指標大幅下降，M>1.5 m時經濟指標下降緩慢，通過優選方案得到選用進路斷面l=1.4 m，M=3.6 m，h=1.40 m，該結果較傳統安全系數法分析結果充填成本降低1.64元/m3。

[1] 楊紅偉,趙 彬,蒲成志,等.下向進路人工假頂失穩機理研究[J].現代礦業,2009(12):13-15. Yang Hongwei，Zhao Bin，Pu Chengzhi，et al.Research on instability mechanism of artificial roof of underhand drift[J].Modern Mining,2009(12):13-15.

[2] Souza E D E,Degagne D，Archibald J F.Minefill applications，practices and trends in Canadian mines[C]∥Proceedings of the 7th International Symposium on Mining With Backfill[S.l.]:The Society for Mining，Metallurgy and Exploration，2001:311-319.

[3] O′Toole D.The basics of mine back-fill[J].Engineering and Mining Journal,2004,205(12):27-29.

[4] 何 書，趙 奎.下向進路側幫可動塊體的穩定性力學分析[J].金屬礦山,2012(12):5-8. He Shu，Zhao Kui.Mechanical analysis of stability of movable block body underhand cut-and-fill drift side[J].Metal Mine，2012(12):5-8.

[5] 劉 冬，金長宇，夏自峰，等.焦家金礦破碎礦體采場進路的跨度優化研究[J].金屬礦山，2014(2):18-21. Liu Dong，Jin Changyu，Xia Zifeng，et al.Research on the stope span optimization of Jiaojia Gold Mine with fractured rock mass[J].Metal Mine，2014(2):18-21.

[6] 趙 彬.焦家金礦尾砂固結材料配比試驗及工藝改造方案研究[D].長沙:中南大學,2009. Zhao Bin.Study of Tailing-cemented Materials Proportion and Back-filling Technology Transformation [D].Changsha:Central South University,2009.

[7] 韓 斌.金川二礦區充填體可靠度分析與1#礦體回采地壓控制優化研究[D].長沙:中南大學,2004. Han Bin.Study on Backfill Reliability Analysis and Ground Control Optimization on 1#Ore-body in Jinchuan No.2 Mine Area[D].Changsha:Central South University,2004.

[8] 孫恒虎，劉文永，黃玉誠，等.高水固結充填采礦[M].北京:機械工業出版社,1998:313-317. Sun Henghu，Liu Wenyong，Huang Yucheng，et al.Mining Technology Backfilled by High-water Solidified Material[M].Beijing:China Machine Press,1988:313-317.

[9] Chadwick J.Benefits of back-fill[J].Mining Magazine,2004(7 ):26-28.

[10] Swan G，Brummer R.Backfill design for deep under-hand drift-and-fill mining[C]∥Stone D.Proceedings of the 7th International Symposium on Mining with Backfill.Littleton:SME，2001:359-368.

[11] 韓 斌,張升學,鄧 建,等.基于可靠度理論的下向進路充填體強度確定方法[J].中國礦業大學學報,2006，35(3):372-376. Han Bin，Zhang Shengxue，Deng Jian，et al.Determining method of backfill strength for underhand drift cut and fill stopping based on reliability theory[J].Journal of China University of Mining & Technology,2006,35(3):372-376.

[12] 江文武,徐國元,中國生.下向分層進路式膠結充填體頂板穩定性分析[J].江西有色金屬,2008,21(4):12-15. Jiang Wenwu，Xu Guoyuan，Zhong Guosheng.Analysis on stability of cemented fills ceiling of underhand delamination heading[J].Jiangxi Nonferrous Metals,2008,21(4):12-15.

[13] 黃玉誠,孫恒虎.高水固結充填下向進路假頂參數設計方法[J].有色金屬:礦山部分,2001(2):15-17. Huang Yucheng，Sun Henghu.An underhand drift roof design method with high water filling consolidation[J].Nonferrous Metal:Mine Section,2001(2):15-17.

[14] 陳玉民,李夕兵.海底大型金屬礦床安全高效開采技術[M].北京:冶金工業出版社,2013:196. Chen Yumin，Li Xibing.Research on Technology of Efficient and Safe Mining in Large-undersea Metal Deposit[M].Beijing：Metallurgical Industry Press,2013:196.

(責任編輯 徐志宏)

Safety Coefficient Analysis in Underhand Drift Section Based on Economic Index

Ji Yonggang1Peng Kang2,3Zhao Jie1Li Jiyu1
(1.XilinguoleAerhadaMiningCo.，Ltd.，ShandongGoldGroup，EastUjimqinBanner026300，China;2.StateKeyLaboratoryofCoalMineDisasterDynamicsandControl，ChongqingUniversity，Chongqing400044，China;3.CollegeofResourcesandEnvironmentalScience，ChongqingUniversity，Chongqing400044，China)

To achieve an efficient and safe operation in Underhand drift stoping，the mechanical model of thin plate with downward drift was built by using elasticity mechanics and a criterion for tensile strength failure of filling body was obtained.Then the economic model of downward drift stopingP(h/M) was built.In the case of a gold mine in Shandong，the relationship between height and width of different drifts and economical index of downward drift filling was analyzed when safety coefficient is 1.8.Research shows that the relationship between economical index and drift width is in linear growth，while the relationship between economical index and drift height is a negative power function.WhenM≤1.5 m，economical index drops significantly.While economical index falls slowly asM>1.5 m.Therefore，the economical index is analyzed when a half width of drift ranges between 1.3 m and 1.7 m，and height ranges between 2.8 m and 3.6 m.It is recommended that the analysis results of safety factor，l=1.5 m，M=3.6 m，h=1.49 m，should be substituted by analysis results of economical indexesl=1.4 m，M=3.6 m，h=1.40 m，for the filling cost has decreased by 1.64 yuan per cubic meter.Therefore，this research has great theoretic value and realism significance.

Thin plate model，Downward drift，Safety coefficient，Economic index

2014-11-04

紀永剛(1970—)，男，高級工程師。 通訊作者 彭 康(1986—)，男，講師，博士。

TD853

A

1001-1250(2015)-01-034-05