王海文,張 寶,顧生春,李志軍,吳亞飛
(1.西部黃金克拉瑪依哈圖金礦有限責任公司, 新疆 克拉瑪依市 834025;2.長沙礦山研究院有限責任公司, 湖南 長沙 410012)
新疆某金礦賦存于上石炭統(tǒng)含碳泥質灰?guī)r夾薄層灰?guī)r條帶中,礦體嚴格受此層位控制,礦體與圍巖界線清楚。礦體由條帶狀、透鏡狀或薄層狀組成,呈層狀產出。產狀與圍巖基本一致,總體走向近東西向。該脈礦體中等穩(wěn)固、上盤圍巖不穩(wěn)固、下盤圍巖中等穩(wěn)固,礦體傾角65°,平均品位6.2 g/t,厚度達2~8 m,屬于急傾斜高品位中等穩(wěn)固中厚礦體。礦石松散系數(shù)為1.6,體重為2.85 t/m3,巖石體重為2.70 t/m3。
礦山采用平硐+豎井聯(lián)合開拓。對于該高品質礦體,目前采用上向進路充填法進行開采。設計采場長度為50 m,寬度為礦體厚度,高度為中段高度50 m,分段高度10 m,分層高度2.5 m。該方法回采貧損率<8%。
由于采場上盤破碎,給礦山安全管理帶來了極大的挑戰(zhàn)。為了預防采礦過程上盤圍巖剝落,必須了解上下盤圍巖巖體質量等級。通過對該金礦上盤巖體質量進行調查,發(fā)現(xiàn)其巖芯破碎呈短柱狀、碎塊狀為主,裂隙極發(fā)育,常見綠泥石化,淺部巖石中發(fā)育有風化裂隙,巖體質量指標<0.12,抗壓強度<30 MPa,巖體較破碎,厚度巨大,巖體質量分級為Ⅳ類,RQD分級、RMR分級、Q系統(tǒng)分級分別為 32,30,2.9,上盤巖體質量總體為差巖體。設計對開挖后進路上盤加以支護,采用螺紋鋼錨桿+噴漿+雙筋條進行聯(lián)合支護,以確保在回采期間上盤的穩(wěn)定性。即首先對進路上盤及頂板進行噴漿(2~3 cm)處理,然后對采場上盤圍巖進行螺紋鋼錨桿(全長錨固)+雙筋條支護,再進行一次噴漿支護(2~3 cm),服務周期約1個月。
為了獲得科學合理的采場破碎上盤支護參數(shù),根據(jù)現(xiàn)場情況及類似巖性支護參數(shù),設計兩種支護方案:第一種方案的錨桿網度為2 m×2 m,錨桿長度為2 m,噴射混凝土厚度為10 cm;第二種方案的錨桿網度為1 m×1 m,錨桿長度為3.5 m,噴射混凝土厚度為4~6 cm。
本文將采用數(shù)值仿真模擬方法,從理論上對這兩種支護參數(shù)進行對比研究,優(yōu)選出一種支護參數(shù)方案,以獲得更適合礦山實際情況的安全、可靠、經濟的支護方案。
1.2.1 支護結構模型
數(shù)值分析軟件采用巖土計算專業(yè) FLAC3D。在建立了采場數(shù)值模型后,通過其FISH程序命令語言添加錨桿和錨索(cable)單元、噴射混凝土(line二步驟)單元,錨桿和噴射混凝土單元需在分層礦體開挖后一排排安裝在上盤,首先噴漿,然后安裝錨桿和雙筋條(局部采用鋼絲網),最后再噴射混凝土。
結合礦山所使用的支護形式(樹脂錨桿,噴漿支護)及其具體的支護參數(shù),參考類似礦山同類型支護的參數(shù),得出礦山支護數(shù)值模擬分析的FLAC3D參數(shù)如表1所示。同時,設置錨桿的預應力為6×104Pa。
表1 支護結構單元參數(shù)
1.2.2 開采步驟數(shù)值模擬
根據(jù)上向進路充填法回采技術特點,為盡可能模擬礦體回采-支護-充填過程中應力、應變的疊加演變,共設計 18次“采-支-充”循環(huán)進行模擬計算。自下向上依次劃分為 18個分層,每個分層高為2.5 m。
(1)模擬采場拉底時,一次性完成“開挖1”達到控頂高度,隨后采用錨桿和混凝土支護上盤圍巖,最后采用高配比混凝土“充填1”。
(2)完成“充填1”后,繼續(xù)開采2,即完成“開挖2”;隨后采用錨桿和混凝土支護上盤圍巖,最后采用高配比混凝土“充填 2”形成本階段的假底。
(3)依次循環(huán)向上發(fā)展直至本階段回采結束。同時,在進行“采-支-充”循環(huán)模擬分析時,記錄每一分層開挖平衡后的最大壓應力、最大拉應力、最大垂向位移和塑性區(qū)變化情況。
因本次需重點分析穩(wěn)定性的對象為上盤的圍巖,為證明錨桿支護的效果和選擇最佳的支護參數(shù),本節(jié)對4 m、6 m厚礦體開采時在3種支護參數(shù):不支護、2 m×2 m網度支護、1 m×1 m網度支護條件下,探討上盤的位移、應力和塑性區(qū)情況。
(1)上盤應力演變。開采4 m厚礦體完成整體開挖后,最大主應力即最大拉應力出現(xiàn)在采場的底板和底板的位置,且隨著網度的增加,頂?shù)装宄霈F(xiàn)的拉應力集中區(qū)變大,但拉應力變小,由于數(shù)值不大,對上盤的穩(wěn)定性影響不大。而隨著支護網度的增加,礦體的上盤和下盤的應力集中區(qū)域面積逐漸減小,數(shù)值逐漸增加,應力集中現(xiàn)象明顯。由圖1,圖2可以看出,對于不支護或2 m網度的支護結構,隨著回采的進行,上盤拉應力逐漸增加,但對于1 m網度的支護結構,上盤拉應力卻在逐漸減小,開挖完成后圍巖內的最大應力為0.15 MPa,說明這種支護結構能有效釋放上盤內的拉應力。但壓應力卻較另外兩種方案有所增加,且在上盤內造成了十分明顯的壓應力集中區(qū)域,但盡管壓應力有集中現(xiàn)象,但絕對數(shù)值較小,在8.7~9.2 MPa之間,對上盤的穩(wěn)定性構不成威脅,可忽略。
圖1 不同支護方式開采過程中最大拉應力演變
圖2 不同支護方式開采過程中最大壓應力演變
(2)上盤垂直位移演變。由圖 3可知,不支護和2 m網度的支護結構上盤中的最大位移基本相同,為28~36 mm,說明2 m的支護結構不能有效控制上盤變形。1 m網度的支護結構上盤下沉位移保持在11.6 mm左右,控制了采場周邊圍巖的變形,維持了上盤的穩(wěn)定性。
圖3 不同支護方式開采過程中最大垂向位移演變
(3)上盤圍巖塑性區(qū)。根據(jù)數(shù)值模擬結果,隨著采場不斷回采,支護結構可以有效的減小圍巖的塑性區(qū),但2 m網度的上盤圍巖還是出現(xiàn)了成片的塑性區(qū),這代表此種支護結構不能最大限度保持上盤圍巖的穩(wěn)定性,而1 m網度的支護結構上盤圍巖內基本沒有塑性區(qū),在下盤內同樣沒有明顯的塑性區(qū),這代表圍巖的穩(wěn)定性良好,小網度支護結構可以有效控制上盤圍巖的變形,使之保持穩(wěn)定。
(1)上盤應力演變。開采6 m厚礦體完成整體開挖后,最大主應力即最大拉應力。由圖 4、圖5可以看出,對于不支護或2 m網度的支護結構,隨著回采的進行,上盤拉應力逐漸增加,但對于1 m網度的支護結構,上盤拉應力卻在逐漸減小,開挖完成后圍巖內的最大應力小于0.2 MPa,這說明這種支護結構能有效釋放上盤拉應力。但壓應力卻較另外兩種方案有所增加,基本保持在10 MPa左右,且在上盤內造成了十分明顯的壓應力集中區(qū)域。盡管壓應力有集中現(xiàn)象,但絕對數(shù)值較小,對上盤的穩(wěn)定性構不成威脅。
(2)上盤豎直位移演變。從圖 6可知,不支護和2 m網度的支護結構上盤中的最大位移基本相同,為35~42 mm,說明2 m的支護結構不能有效控制上盤變形。1 m網度的支護結構上盤下沉位移,保持在16 mm左右,控制了采場周邊圍巖的變形,維持了上盤的穩(wěn)定性。
圖4 不同支護方式開采過程中最大拉應力演變
圖5 不同支護方式開采過程中最大壓應力演變
圖6 不同支護方式開采過程中最大垂向位移演變
(3)上盤圍巖塑性區(qū)。根據(jù)數(shù)值模擬結果,隨著采場不斷回采,支護結構可以有效減小圍巖的塑性區(qū),但2 m網度的上盤圍巖還是出現(xiàn)了成片的塑性區(qū),這代表此種支護結構不能最大限度保持上盤圍巖的穩(wěn)定性,而1 m網度的支護結構上盤圍巖內基本沒有塑性區(qū),在下盤內同樣沒有明顯的塑性區(qū),這代表圍巖的穩(wěn)定性良好,小網度支護結構可以有效控制上盤圍巖變形,使之保持穩(wěn)定。
(1)對于不支護或2 m網度的支護結構,隨著回采的進行,上盤內拉應力逐漸增加。但對于1m網度的支護結構,上盤內的拉應力會逐漸減小,此種支護結構能釋放上盤內拉應力,同時壓應力卻較另外兩種方案有所增加,上盤內有明顯的應力集中區(qū),但絕對數(shù)值較小,約為9 MPa,不會對上盤的穩(wěn)定性構造成影響。
(2)不支護和2 m網度的支護結構上盤中的最大位移基本相同,平均為35 mm左右,基本為1 m網度支護結構的上盤位移的2倍,2 m的支護結構不能有效控制上盤變形。1 m網度的支護結構能有效減小上盤圍巖的下沉位移,控制了圍巖變形,保持上盤的穩(wěn)定性。
(3)支護結構可以有效的控制圍巖的塑性區(qū),但2 m網度的上盤圍巖內會出現(xiàn)成片的塑性區(qū),表明此種支護結構不能保持上盤圍巖的穩(wěn)定性。而 1 m網度的支護結構上下盤圍巖內均沒有塑性區(qū),只有充填體內出現(xiàn)零星塑性區(qū),對上盤穩(wěn)定性影響不大。這表明1 m網度的噴漿+錨桿+雙筋條支護結構可以有效地控制上盤圍巖變形,保持其穩(wěn)定。
因此,為了保證上盤圍巖的穩(wěn)定性,確保采場安全高效回采,推薦上向進路充填法采場上盤圍巖的支護參數(shù)為:錨桿支護網度1 m×1 m,錨桿長度3.5 m,噴射混凝土厚度4~6 cm。