張愛卿， 李金云， 吳愛祥， 王貽明

（1.北華航天工業(yè)學院 建筑工程學院，河北 廊坊 065000； 2.北京科技大學 膏體充填采礦技術(shù)研究中心，北京 100083）

隨著礦產(chǎn)資源日益枯竭以及礦山生產(chǎn)的迫切需求，礦柱的安全高效回采已成為許多礦山企業(yè)面臨的難題。 礦房內(nèi)礦石開采后采空區(qū)殘留礦柱起到支撐空區(qū)頂板的作用，導致礦柱應力集中，尤其是一些特殊工況下，比如該中段的上下中段礦房均已開采、只有上下礦段對應的礦柱存在時，礦柱的回采形式就成為保證礦柱安全高效回采的關(guān)鍵。 合理的礦柱回采方式對于提高礦床資源采出率、實現(xiàn)礦山可持續(xù)發(fā)展具有重要意義［1－5］。

近年來，專家學者對礦柱的回采形式進行了大量研究［6－9］，但都是針對嗣后充填礦柱連續(xù)回采展開，對于多礦柱同時回采的研究報道相對較少。

本文建立了3個礦柱同時回采模型，總結(jié)文獻研究成果和多年現(xiàn)場實踐經(jīng)驗，確定影響充填體側(cè)向位移的主要因素為礦柱一次回采高度、礦房寬度、礦柱跨度和充填體強度，并結(jié)合某鉛鋅礦實際開采情況，分析了礦柱一次回采高度、礦柱跨度、礦房寬度及膠結(jié)充填體強度對充填體側(cè)向位移的影響規(guī)律，最終確定了該鉛鋅礦適宜的礦柱一次回采高度和充填體強度值。 研究成果對其他類似金屬礦山進行礦柱回采具有一定借鑒意義。

1 嗣后充填多礦柱同時回采的有限元模型

嗣后充填多礦柱回采示意圖如圖1 所示，圖中h為礦柱一次回采高度，b為礦柱跨度，B為礦房寬度。本文提及的充填體側(cè)向位移，對于1＃礦柱左側(cè)的充填體僅考慮其靠近空區(qū)的一側(cè)，對于2＃礦柱左側(cè)的充填體考慮兩側(cè)的位移值，由于該模擬是對稱的，分析僅針對模型左側(cè)，同時，定義1＃礦柱左側(cè)充填體的側(cè)向位移值為單側(cè)位移，2＃礦柱兩側(cè)充填體的側(cè)向位移值為雙側(cè)位移。

圖1 嗣后充填多礦柱回采示意圖

根據(jù)圖1 建立模型并劃分網(wǎng)格，以垂直礦體走向為X軸方向，沿礦體走向為Y軸方向，模型的高為Z軸方向。 運用有限元軟件FLAC3D中Mohr-Coulomb 模型進行彈塑性分析，由于礦山賦存礦體和圍巖巖體是一個復雜地質(zhì)構(gòu)造體，為便于計算且快速收斂，對巖體介質(zhì)性質(zhì)及計算模型等做以下假設(shè)：

1） 將礦體和巖體視為各向同性、連續(xù)均質(zhì)的力學介質(zhì)體；

2） 忽略不良地質(zhì)構(gòu)造的影響，如斷層、節(jié)理裂隙等；

3） 計算過程只在靜荷載下進行，不考慮地下水、地震及爆破振動等動力對空區(qū)穩(wěn)定性的影響。

考慮設(shè)計開采范圍周圍巖體受采動影響范圍的大小，將模型尺寸擬定為長410 m、寬210 m、高600 m，建立的三維數(shù)值模型如圖2 所示。 初始模型采用四面體單元剖分，重點剖分礦柱部分，精度以確保無畸變單元為原則并在局部適當加密。 數(shù)值模擬離散后的模型六面體單元數(shù)共140 000個網(wǎng)格。

圖2 FLAC3D數(shù)值模擬模型

數(shù)值計算采用文獻［10］中的動態(tài)強度折減法，僅對充填體的內(nèi)摩擦角和黏聚力值按照一定比例折減，分析多礦柱回采后兩側(cè)充填體水平位移場的變化規(guī)律。

2 工程實例

某鉛鋅礦礦體采用淺孔留礦法進行開采，礦房垂直礦體走向布置，目前開采工作面在860 m 中段，垂直距離地表深度235 m，開采順序為：先采礦房，形成礦柱支撐采空區(qū)，對空區(qū)充填后再回采礦柱（見圖1）。礦房結(jié)構(gòu)參數(shù)為長35 m、寬15 m、高50 m，礦柱結(jié)構(gòu)參數(shù)為長35 m、寬9 m、高50 m。

參考凡口鉛鋅礦及國內(nèi)其他金屬礦山有關(guān)嗣后充填礦柱回采時選取的參數(shù)，假設(shè)礦柱跨度分別取值6 m、9 m 和12 m，礦房寬度分別取值12 m、15 m 和20 m，充填體強度分別取值1.5 MPa、2 MPa 和2.5 MPa，礦柱一次回采高度分別取值4 m、6 m 和8 m。

2.1 礦柱及充填體物理力學參數(shù)

某鉛鋅礦上盤、下盤巖性均為安山質(zhì)凝灰?guī)r，礦體樣本為含銀鉛鋅原生硫化礦礦石，通過現(xiàn)場取芯測試其物理力學參數(shù)，結(jié)果見表1。

表1 某鉛鋅礦巖石物理力學參數(shù)

通過對礦區(qū)內(nèi)外部開采技術(shù)條件進行綜合技術(shù)經(jīng)濟分析，確定礦山的充填方式為全尾砂膠結(jié)充填。 全尾砂從選廠取出運送到實驗室，采用PC32.5 級復合硅酸鹽水泥，根據(jù)影響充填體強度的主要因素、坍落度實驗觀察到的料漿流動情況，設(shè)計灰砂比1 ∶4、1 ∶6、1 ∶8、1 ∶10、1 ∶12 共5 組，料漿濃度分別取76%、74%、72%、70%共4 組，合計20 組不同材料配比實驗。 選取充填體28 d 抗壓強度為1.5 MPa、2 MPa 和2.5 MPa 的物理力學參數(shù)，見表2。

表2 充填體物理力學參數(shù)

2.2 充填體側(cè)向位移影響因素分析

2.2.1 一次回采高度的影響

假定礦柱高度50 m、跨度9 m、礦房寬度15 m、充填體強度2 MPa，以一次回采高度h＝4 m、6 m 和8 m為例，研究了一次回采高度對充填體側(cè)向位移的影響，選取礦柱回采后每隔1 m 所對應的充填體側(cè)向位移數(shù)值繪制成圖，如圖3 所示。 從圖3可以看出，充填體雙側(cè)側(cè)向位移較單側(cè)側(cè)向位移值大，不論礦柱回采高度是否變化，充填體側(cè)向位移變化均呈凸拋物線形分布，由礦柱靠近中央位置向兩端逐漸減小。 充填體側(cè)向位移隨著礦柱一次回采高度增加而增大。 對于全尾膠結(jié)充填體來說，產(chǎn)生的側(cè)向位移值越小越好。 因此，一次回采高度4 m 是安全的回采高度。

圖3 礦柱一次回采高度對充填體側(cè)向位移的影響

2.2.2 礦柱跨度的影響

假定礦柱高度50 m、一次回采高度6 m、礦房寬度15 m、充填體強度2 MPa，以礦柱跨度b＝6 m、9 m 和12 m 為例，研究了礦柱跨度對充填體側(cè)向位移的影響，結(jié)果見圖4。 從圖4可以看出，充填體雙側(cè)側(cè)向位移較單側(cè)側(cè)向位移值大，充填體側(cè)向位移呈凸拋物線形分布，由礦柱靠近中央位置向兩端逐漸減小。 充填體側(cè)向位移值隨著礦柱跨度增加而增大，呈正相關(guān)變化規(guī)律。

圖4 礦柱跨度對充填體側(cè)向位移的影響

2.2.3 礦房寬度的影響

假定礦柱高度50 m、一次回采高度6 m、礦柱跨度9 m、充填體強度2 MPa，以礦房寬度B＝12 m、15 m 和20 m 為例，研究了礦房寬度對充填體側(cè)向位移的影響，結(jié)果如圖5 所示。 從圖5可以看出，充填體雙側(cè)側(cè)向位移較單側(cè)側(cè)向位移值大，充填體側(cè)向位移呈凸拋物線形分布，由礦柱靠近中央位置向兩端逐漸減小。充填體側(cè)向位移值隨著礦房寬度增加而減小，呈負相關(guān)變化規(guī)律。

圖5 礦房寬度對充填體側(cè)向位移的影響

2.2.4 充填體強度的影響

假定礦柱高度50 m、一次回采高度6 m、礦房寬度15 m，礦柱跨度9 m，以充填體強度P＝1.5 MPa、2 MPa和2.5 MPa 為例，研究了充填體強度對充填體側(cè)向位移的影響，結(jié)果如圖6 所示。 從圖6可以看出，充填體雙側(cè)側(cè)向位移較單側(cè)側(cè)向位移值大，充填體側(cè)向位移呈凸拋物線形分布，由礦柱靠近中央位置向兩端逐漸減小。 充填體側(cè)向位移值隨著充填體強度增加而減小，呈負相關(guān)變化規(guī)律。

圖6 充填體強度對充填體側(cè)向位移的影響

2.2.5 極差分析結(jié)果

對上述4 組數(shù)據(jù)分別求極差，一次回采高度、礦房寬度、礦柱跨度和充填體強度的極差分別為0.75、0.9、0.85 和0.55。 根據(jù)極差分析判斷上述4個影響因素對充填體側(cè)向位移產(chǎn)生的敏感程度由大到小為：礦房寬度、礦柱跨度、一次回采高度和充填體強度，如果將充填體作為礦柱來考慮的話，充填體側(cè)向位移敏感程度的排序與文獻［7］研究結(jié)論一致。

2.3 工程實例中P 和h 的確定

大量工程實踐表明，充填料漿沉縮后導致充填體接頂困難［11－12］，因此，采用高水自立強度公式來確定充填體強度P［9］：

式中P為充填體自立強度，MPa；L為礦房長度，m；H為礦房高度，m；W為礦房寬度，m；γ為充填體密度，g/cm3；C為充填體黏結(jié)強度，MPa；φ為充填體內(nèi)摩擦角，（°）；α為充填材料坍塌角，（°），α＝45°＋φ/2；k為側(cè)壓系數(shù)；φ1為圍巖摩擦角，（°）；C1為圍巖黏結(jié)強度，MPa。

將礦房結(jié)構(gòu)參數(shù)和表2 數(shù)據(jù)代入式（1）中，得出充填體強度2.5 MPa 時，可以滿足充填體自立要求。

假定3個礦柱同時一次回采高度分別為4 m、6 m和8 m，通過數(shù)值計算，得到充填體側(cè)向位移（雙側(cè)）隨一次回采高度變化規(guī)律如圖7 所示。 從圖7可以看出，充填體側(cè)向位移與礦柱一次回采高度呈正相關(guān)。一次回采高度4 m 時，充填體側(cè)向位移最小，從礦柱回采安全角度考慮，確定某鉛鋅礦的一次回采高度為4 m。

圖7 不同一次回采高度下充填體雙側(cè)側(cè)向位移變化規(guī)律

3 結(jié) 論

1） 3個礦柱同時回采時，充填體雙側(cè)側(cè)向位移值大于單側(cè)側(cè)向位移值，充填體側(cè)向位移最大值出現(xiàn)在礦柱高度中間位置，由礦柱靠近中央位置向兩端逐漸減小。 最大側(cè)向位移值隨著一次回采高度增大而增加，隨著礦柱跨度增加而增大，與礦房寬度和充填體強度呈負相關(guān)。

2） 極差分析結(jié)果表明，一次回采高度、礦柱跨度、礦房寬度和充填體強度4個影響因素對充填體側(cè)向位移產(chǎn)生的敏感程度由大到小排序為：礦房寬度、礦柱跨度、一次回采高度和充填體強度。

3） 對于某鉛鋅礦，當?shù)V房、礦柱結(jié)構(gòu)參數(shù)固定時，利用高水自立強度公式計算出充填體自立強度值為2.5 MPa，從礦柱回采安全性考慮，最終確定3個礦柱同時回采時一次回采高度為4 m。