膠結充填體下采場結構參數優化研究

李龍福 江東平

(1.中鋼集團馬鞍山礦山研究院有限公司，安徽 馬鞍山 243000；2.金屬礦山安全與健康國家重點實驗室，安徽 馬鞍山 243000)

膠結充填體下采場結構參數優化研究

李龍福1,2江東平1,2

(1.中鋼集團馬鞍山礦山研究院有限公司，安徽 馬鞍山 243000；2.金屬礦山安全與健康國家重點實驗室，安徽 馬鞍山 243000)

要在大體積膠結充填體下對礦體進行安全高效地開采，合理的采場結構是采礦設計首先需要考慮的問題。在綜合分析某地下礦山地質概況、開采技術條件現狀的基礎上，采用經驗法與理論計算法對采場結構參數進行計算，提出了4種采場結構參數初步方案，建立對應的三維數值計算模型，將圍巖的變形、應力分布及塑性區大小等作為衡量采場結構參數優劣的指標。根據模擬結果對比分析了各方案在礦房回采后采場圍巖的力學響應情況，得到了采場處于最有利力學狀態時的結構參數，結合各方案盤區礦石理論回采率，進而對采場結構方案進行了優選，即當頂柱厚度為6 m，礦房跨度38 m，礦柱寬14 m時，采場結構穩定性處于最有利狀態并且具有較高的理論回采率。優化結果可為后續開采設計提供參考。

膠結充填體 采場結構參數 數值優化 圍巖力學響應

地下礦山要安全高效地生產，采場結構參數是否合理是其開采前提[1]，而對于某一特定條件下的礦體開采，其結構參數的確定應視具體情況具體分析。目前國內外對采場結構參數優化選擇的方法主要有經驗類比法、半經驗半理論法、不確定性方法、數值模擬法、模糊數學理論、人工神經網絡理論等。而隨著計算機科學技術的發展及數值理論分析的日趨成熟，各種功能強大的數值模擬軟件得到迅猛發展，使得三維數值模擬計算成為采場結構參數優化的有效手段之一[2-4]。其中，三維數值模擬可以定量地計算和分析回采過程中圍巖的應力、位移和塑性區的分布狀況，分析它們在回采過程中的動態變化，從而對采場圍巖的穩定性狀態做出判斷[5-6]。

本研究應用FLAC3D三維有限差分軟件，結合某地下礦山開采現狀，建立相關的三維數值模型，通過多方案的計算和分析，比較在不同方案下圍巖的力學響應特征，從而確定出安全合理的采場結構參數。

1 工程概況

某礦山前期由于資金、技術人員等方面實力比較薄弱，所以采用落后、短期投入的開采技術。由于上述原因，目前在礦區地表、井下形成了大量的采坑(約92萬m3)和未充填的采空區(約12.57萬m3)。這些空區受礦體自然形態及開采條件的影響，空間形態分布各異，而且由于巖體的力學特性呈各向異性，導致巖體內天然應力場發生改變，應力場重新調整，在一些部位發生應力集中、局部釋放，最終導致圍巖的失穩破壞。2010年礦山對前期已經回采的3個中段中殘留的大量礦柱進行回收及對采空區進行膠結充填處理，由此在已經回采中段中會留下大體積的充填體，這對于下層礦體的安全開采提出了新的要求，即合理布置采場結構，優化選擇采場結構參數。

本研究選擇礦區南礦段規模最大的5#礦體作為研究對象，礦體分布于南礦段的西南部，地表控制礦體長219 m，坑道工程控制礦體長319 m，垂深已控制180 m，在-180 m中段礦體厚2.07～27.79 m，平均厚15.02 m，變化系數0.67，礦體形態復雜程度中等。在地表處傾向東，至地下漸變為總體傾向西，傾向250°～335°，傾角52°～81°，走向8°～46°。

礦體圍巖主要為礫巖，其次為花崗巖和炭質頁巖。頂柱為礫巖時，巖石呈塊狀，結構穩定，取樣深度180 m，其飽和抗壓強度為59.5～146 MPa，穩定性較好；頂柱為炭質頁巖時，巖體呈層狀，結構不穩定，其飽和抗壓強度為1.06 MPa，穩定性一般；頂柱為蝕變黑云母斜長花崗巖時，綠泥石化較強烈，小節理裂隙發育，巖體聯結能力弱，巖石無水狀態下較堅硬，遇水很快變軟，力學性能強度較低，會發生崩塌，穩定性較差。本研究涉及的相關材料的物理力學參數如表1。

表1 材料物理力學參數

2 采場結構參數的初步確定

-180 m中段以下礦體的開采方案為空場嗣后充填采礦法，即先采礦房，回采完一個中段后全部充填，之后再回收礦柱。在此過程中采場的穩定情況大致可歸納為以下4種類型[7]：礦柱穩定、頂柱巖層穩定；礦柱穩定、頂柱巖層局部破壞；礦柱失穩、頂柱巖層完整；礦柱失穩、頂柱巖層破壞。而與此有關因素主要包括頂柱厚度、礦房長度、礦柱寬度和礦柱面積比率及巖體強度等，因此，主要對頂柱厚度、礦房跨度、礦柱寬度進行優化選擇，其初始值確定如下。

(1)頂柱厚度。頂柱厚度的確定采用經驗類比法。根據經驗，在空場法開采中，中厚—厚礦體的開采，頂柱厚度一般取3～6 m，5#礦體平均厚度為15.02 m，屬于厚礦體，因此，頂柱厚度分別取3 m、4 m、5 m、6 m，并以此作為其他參數確定的依據。

(2)礦房跨度。由于礦體為全厚開采，所以其寬度為礦體的厚度，取平均值。采場極限跨度L可通過“巖梁”理論[8]來計算：

L=4H2Rs/3q(H-2),

(1)

式中，Rs為巖體的抗剪強度，Pa；q為上覆巖體的均布載荷，kN/m2；H為巖梁的厚度，m。因涉及的對象為-180 m、-220 m中段結構參數，中段高差40 m，應力增幅為1.08 MPa，遠小于巖體的強度，因此，將兩中段的參數作統一考慮，即上覆巖體的均布載荷取兩者的平均值。代入相應參數，計算結果為44 m、39 m、41 m、44 m。

另外，針對頂柱圍巖的力學分布結構，考慮頂柱上覆巖層的壓力、巖體的力學特征等與礦房跨度的關系，得出空場法中礦房跨度b的理論計算公式[9]：

(2)

式中，λ為側壓力系數；γ為礦石的容重，kN/m3；H為開采深度(頂柱距地表深度)，m；α1為頂柱的許可抗拉強度，Pa。代入相應參數，計算結果為32 m。

綜合考慮各種因素，取兩者的平均值，分別為38 m、35 m、36 m、38 m。

(3)礦柱寬度。根據Bieniawski[10]的研究，礦柱寬度與安全系數幾乎成正比關系，寬度越大，安全系數越高，進而可以確定出礦柱的合理寬度，且存在如下關系：

(3)

式中，K為礦柱安全系數；σc為礦巖的強度，MPa；h為礦房高度，m；β為常數，取值根據礦柱的寬高比值而定；Lp為礦柱截面長，m；L0為礦柱間隔的長度，m；Wo為礦房寬度，m；Wp為礦柱寬度，m；γ為上覆巖層的容重，N/m3。根據Bieniawski等的實驗研究成果，當礦柱的寬高比大于5時，β=1.4；而當礦柱的寬高比小于5時，β=1。代入參數，得安全系數與礦柱寬度的關系如圖1所示。

圖1 安全系數與礦柱寬度的關系

考慮到現場實際條件的影響，一般安全系數取值在1.2以上，對應于圖1所示關系，礦柱寬度應分別大于等于14 m、12 m、13 m、14 m。

綜合以上分析，可得4個結構參數組合方案，如表2所示。

表2 采場結構參數初步方案

3 模型的建立及開挖方案的確定

3.1 三維模型的建立

考慮到計算量和計算效果，取-180 m中段控礦范圍的一半作為計算區域，即盤區的一半160 m，中段高40 m，礦體厚度取平均值15 m，礦體傾角取平均值65°，共可布置3個礦房。

根據彈塑性理論，在地下工程中，巖體的開挖僅對距開挖中心點3～5倍跨度范圍內的圍巖產生影響，因此，計算模型范圍：垂直礦體走向為x軸，取210 m；沿礦體走向為y軸，取500 m；垂直方向為z軸，為進一步優化計算模型，將垂直方向從+995 m水平取至+790 m水平(地表標高為+1 190 m)。得出模型長×寬×高為500 m×210 m×205 m。

計算邊界采用位移-應力混合約束，將模型4個側面法線方向的水平位移進行固定(ux=0,uy=0)，對模型底部平面位移固定(ux=0,uy=0,uz=0)，上部邊界為上覆巖體的自重應力，σz=-2.57 MPa。采用Mohr-Coulomb屈服準則。

3.2 開挖方案的確定

根據礦山采礦方法中采場結構布置，模擬計算方案分3步完成：礦房1回采→礦房2回采→礦房3回采，如圖2所示。

圖2 采場結構示意

3.3 模擬結果分析

由于研究礦體埋深較淺，原巖應力場只考慮自重應力場的作用，水平應力根據泊松效應進行計算。選取礦房全部回采后，采場中圍巖變形、應力分布、盤區礦石理論回采率及塑性區大小作為衡量采場結構參數優劣的指標。具體見表3所示。

表3 不同方案模擬結果對比

(1)計算結果表明，在礦房回采結束后，最大垂直位移、最大水平位移均發生在采場的頂板中。從表3的數據可知，當頂柱厚度由3 m增加到6 m時，頂柱中的位移顯著減小，相對于方案一減小幅度分別為7.4%、13.5%、19.0%；5.7%、15.9%、24.0%。說明頂柱在控制采場變形方面起著重要作用。

在上部充填體中，如圖3所示，礦房回采后位移在垂直方向呈現拱形分布，礦房正上方位移比周邊大，由中心向外逐漸減小；并且從云圖中發現，當頂柱厚度在5 m以上時，礦房正上方的位移等值云圖中較大變形區沒有發生疊加現象，說明在該厚度下，頂柱可以較好地控制位移的傳遞，阻止變形在充填體中的擴展，維護采場穩定。

圖3 垂直方向位移云圖

(2)隨著礦房的回采，采場圍巖中出現了不同程度的應力集中現象，以礦房底柱中最為顯著，礦柱與頂板中次之，如圖4所示。最大主應力值相差不大，最大值出現在方案一中，為-10.17 MPa，遠小于巖體的抗壓強度，因此，采場處于較小的壓應力場中。從充填體中應力分布情況分析，頂柱厚度不同，充填體的受力狀態有差異，如圖4所示，隨頂柱厚度的增加，圍巖與充填體的受力更趨一致，說明當頂柱達到一定厚度時，可以把礦柱的支撐效果幾乎均勻地分配到整條巖梁上，從而表現出整體變形模式，而不會出現由于頂柱強度不足而出現的分段彎曲，表現在礦柱正上方撓度小，其他地方如礦房中部撓度大等特征。因此，要使充填體與圍巖作為一個整體受力，頂柱的厚度是一個重要的影響因素。

圖4 最大主應力云圖

對比最小主應力值可以發現，在方案二中出現120 kPa大小的拉應力，雖小于圍巖的抗拉強度，但巖體作為一種特殊的脆性材料，且其中存在大量的微裂隙、不連續面，這些結構弱面在很小的拉應力作用下就發生斷裂，因此，一般情況下要盡可能的避免出現拉應力。其余方案中均沒有出現拉應力，但從最小主應力值發現，隨頂柱厚度的增加，最小主應力的絕對值逐漸增大，說明圍巖的受力條件更好。

(3)由圖5可知，采場下盤圍巖與礦房底柱中分別發生大面積拉伸破壞和剪切破壞，方案三、方案四的破壞范圍稍小于方案一、方案二；在頂柱中，方案三、方案四中只有局部小范圍的拉伸破壞，而在方案一與方案二中除了有拉伸破壞，還出現剪切破壞，說明頂柱中的剪應力較大，局部超過巖體的剪切強度；在礦柱中塑性區深度幾乎相同。

圖5 塑性區分布云圖

整體而言，4種方案的塑性區情況類似，在礦柱與頂柱中均沒有出現貫通性的塑性破壞區，說明在該工況下采場都是安全的，但從下盤圍巖與底柱塑性區范圍大小而言，方案三、方案四要優于方案一與方案二。

經過綜合對比分析，考慮盤區礦石理論回采率，當頂柱厚度為6 m，礦房跨度38 m，礦柱寬14 m時，采場結構穩定性處于最有利狀態并且具有較高回采率，所以將方案四中的參數確定為最優采場結構參數。

4 結 論

(1)利用數值模擬手段，對采用經驗與理論計算相結合確定的采場結構參數初始值進行計算，對比分析了礦房回采后圍巖的變形、應力及塑性區分布情況，考慮盤區礦石回采率，對結構參數進行了優化，即當頂柱厚6 m，礦房跨度38 m，礦柱寬14 m時，采場穩定性處于最有利狀態并且具有較高回采率。

(2)采場垂直方向位移呈現拱形分布，位移由中心向外逐漸減小；當頂柱厚度在5 m以上時，礦房正上方的大變形區沒有發生疊加現象，說明在該厚度下，頂柱可以較好地控制位移的傳遞，阻止變形在充填體中的擴展。垂直與水平方向的最大變形均分布在頂柱中，且從變形的幅值可以看出，頂柱在控制采場變形方面起著重要作用。

(3)充填體的受力情況與頂柱厚度有關，當頂柱達到一定厚度時，可以把礦柱的支撐效果幾乎均勻地分配到整條巖梁上，從而使充填體與圍巖表現出整體變形模式。

(4)雖同為理論計算下的安全參數，但在方案二中卻出現拉應力，顯示出經驗與理論計算等傳統方法的局限性，而三維數值模擬卻可以較好地解決這一問題，且隨著數值仿真技術的發展，將在地下工程結構優化中發揮出更大的作用。

[1] 王新民，李潔慧，張欽禮，等．基于FAHP 的采場結構參數優化研究[J]．中國礦業大學學報，2010，39(2)：163-168． Wang Xinmin，Li Jiehui，Zhang Qinli，et al.Optimizing mining stope structural parameters using a FAHP[J].Journal of China University of Mining & Technology，2010，39(2)：163-168．

[2] 王新民，曹 剛，張欽禮，等.康家灣礦深部難采礦體采場穩定性及結構參數優化研究[J].河南理工大學學報：自然科學版,2007,26(6):634-640． Wang Xinmin，Cao Gang，Zhang Qinli，et al.Research on stability and optimum sizes of deep stope with weak roof in Kangjiawan Mine[J].Journal of Henan Polytechnic University：Natural Science，2007，26(6)：634-640．

[3] 羅周全，管佳林，馮富康，等.盤區隔離礦柱采場結構參數數值優化[J].采礦與安全工程學報,2012,29(2):261-264． Luo Zhouquan，Guan Jialin，Feng Fukang，et al.Stope structural parameters of panel isolation pillar numerical optimization[J].Journal of Mining & Safety Engineering，2012，29(2):261-264．

[4] 崔少東,劉保國.石人溝鐵礦北區采場結構參數的優化[J].北京交通大學學報,2009,33(1):131-134． Cui Shaodong，Liu Baoguo.Stope structure parameters.optimization of North Shirengou Iron Mine[J].Journal of Beijing Jiaotong University，2009，33(1):131-134．

[5] 李佳洋.貢北金礦破碎頂柱下緩傾斜薄至中厚礦體安全開采技術研究[D].長沙:中南大學,2011． Li Jiayang.Study of Safety Mining Technology for Mining Gently Inclined Thin to Medium Thick Ore under Broken Roof in Gongbei Gold Mine[D].Changsha：Central South University,2011．

[6] Hart R，Cundall P A，Lemos J.Formulations of three-dimensional distinct element model PAPTⅡ:Mechanical calculation of a system composed of many polyhedral blokes[J].International Journal of Rock Mechanics and Mining Sciences，1988(3)：117-121．

[7] 黃英華，徐必跟，唐紹輝.房柱法開采礦山采空區失穩模式及機理[J].礦業研究與開發,2009,29(4):24-26． Huang Yinghua，Xu Bigen，Tang Shaohui.Study on the damage patterns and mechanism of mined-out area in mines using room-and-pillar mining method[J].Mining Research & Development，2009，29(4)：24-26．

[8] 柴修偉，張電吉，周 麟，等.磷塊巖礦柱置換安全穩定性研究[J].金屬礦山,2011(7):8-11． Chai Xiuwei，Zhang Dianji，Zhou Lin，et al Research on safety and stability of phosphorite ore pillar replacement[J].Metal Mine，2011(7)：8-11．

[9] 趙 奎，胡慧明，王曉軍，等.某金礦房柱法采場人工礦柱參數選取[J].采礦技術,2011,11(2):15-17． Zhao Kui，Hu Huiming，Wang Xiaojun，et al.Artificial column parameters selected of room and pillar mining in a gold mine[J].Mining Technology，2011，11(2)：15-17．

[10] 廖文景.石膏礦采空區積水對礦柱穩定性的影響分析[J].采礦技術,2009，9(3):52-54． Liao Wenjing.Effect of gypsum mine goaf hydrops on pillar stability analysis[J].Mining Technology，2009，9(3)：52-54

(責任編輯 石海林)

StudyonOptimizationofStopeStructuralParametersunderCementedFill

Li Longfu1,2Jiang Dongping1,2

(1.SinosteelMaanshanInstituteofMiningResearchCo．，Ltd.，Maanshan243000，China;2.StateKeyLaboratoryofSafetyandHealthforMetalMine，Maanshan243000，China)

To safely and efficiently mining the ore under large cemented filling body，reasonable stope structure is the first issue that need to consider in mining design.On the basis of comprehensive analysis on geological situation and current mining technical condition of some underground mine，stope structure parameters were calculated by using empirical method with theoretical calculation，and four kinds of the preliminary scheme about the stope structure were put forward.Then，the corresponding 3-dimensional numerical model is established，where the deformation，stress distribution of surrounding rock and the plastic zone size were taken as indicators to measure the performance of stope structure parameters.Based on the simulation results，the mechanical response of the various programs in the stope wall rock after room mining are comprehensively analyzed，obtaining the structural parameters under the most favorable mechanical state of the stope.Combining with the theoretical recovery rate of each panel，the programs of stope structure were optimized，that is，when the thickness of top column is at 6 m，the stope span at 38 m and the pillar width at 14 m，stope structural stability keeps in the most favorable state and the stope has a high theoretical recovery rate.The optimized results could provide reference for the follow-up mining design.

Cemented fill，Stope structural parameter，Numerical optimization，Mechanical response of surrounding rock

2014-08-04

李龍福(1986—)，男，碩士研究生。

TD853

A

1001-1250(2014)-12-035-05