曲俊霖,王繼仁,韓新平,南存全
(遼寧工程技術(shù)大學(xué) 礦業(yè)學(xué)院,遼寧 阜新 123000)
螺旋鉆機(jī)作為一種可以高效安全地回采露天礦最終境界下端幫壓煤的礦山機(jī)械,已經(jīng)在德國(guó)和澳大利亞有著廣泛且成功的應(yīng)用[1]。螺旋鉆機(jī)運(yùn)用螺旋鉆桿開采出1排或多排的水平鉆孔,從端幫的鄰空面回采壓煤且無需支護(hù),適用于使用單一的露天或井工開采方法難以回采或者經(jīng)濟(jì)效益不理想的情況。目前國(guó)內(nèi)學(xué)者對(duì)螺旋鉆機(jī)開采端幫煤炭的研究多集中在對(duì)多層螺旋鉆孔的采場(chǎng)穩(wěn)定性研究,且研究方法大多為利用數(shù)值模擬軟件分析比較不同方案中參數(shù)的合理性。曾錢幫等[2]運(yùn)用數(shù)值模擬的方法,研究了螺旋鉆機(jī)不同的鉆孔開采順序?qū)Σ蓤?chǎng)穩(wěn)定性的影響。研究表明:層間不同開采順序?qū)γ褐€(wěn)定性影響較大,而同層內(nèi)的鉆孔開挖順序?qū)γ褐€(wěn)定性影響不明顯;黃侃等[3]運(yùn)用數(shù)值模擬方法對(duì)螺旋鉆機(jī)開采下不同水平孔間煤柱寬度及垂向隔層的不同排布方式對(duì)采場(chǎng)穩(wěn)定性的影響進(jìn)行了研究。研究表明:由于采場(chǎng)間垂向隔層為煤柱提供了側(cè)向應(yīng)力,所以在煤柱優(yōu)化設(shè)計(jì)時(shí)需要考慮采場(chǎng)間的垂向隔層厚度;程國(guó)明等[4-5]對(duì)特厚煤層條件下采用螺旋鉆機(jī)開采細(xì)長(zhǎng)窄煤柱時(shí),煤柱的破壞過程進(jìn)行了數(shù)值模擬。研究表明:螺旋鉆機(jī)開采所引起的煤柱內(nèi)應(yīng)力分布以及破壞方式與條帶式分布不同。
鑒于以往鮮有對(duì)單層螺旋鉆孔參數(shù)設(shè)計(jì)的研究,且未兼顧壓力拱這一采場(chǎng)中的“小結(jié)構(gòu)”對(duì)圍巖內(nèi)應(yīng)力分布狀態(tài)、關(guān)鍵層以及永久隔離煤柱內(nèi)應(yīng)力集中程度的影響。筆者以礦山實(shí)例為背景對(duì)單層螺旋鉆孔設(shè)計(jì)參數(shù)進(jìn)行研究,基于壓力拱理論分析數(shù)值模擬結(jié)果,選出最優(yōu)開采方案,希望為螺旋鉆機(jī)開采邊坡煤設(shè)計(jì)方案提供依據(jù)。
螺旋鉆機(jī)開采單層鉆孔與條帶開采在開采形式上相近,都可以簡(jiǎn)單概括為采出一個(gè)條帶后留設(shè)一個(gè)條帶,如此往復(fù)。留設(shè)的條帶煤柱能承載上覆巖層的載荷,保證地表移動(dòng)和變形的均勻且可控。2種開采方式在開采過程中,相鄰采空區(qū)間煤柱內(nèi)應(yīng)力變化曲線有如下相同點(diǎn):①在未進(jìn)行開采時(shí),煤體受上覆巖層均布荷載作用;②當(dāng)煤柱一側(cè)開采結(jié)束后,煤柱一定范圍內(nèi)形成支承壓力帶及局部塑性區(qū);③煤柱兩側(cè)均開采結(jié)束后,若煤柱能保持穩(wěn)定狀態(tài),則煤柱所受垂直應(yīng)力呈馬鞍型分布,煤柱兩側(cè)存在一定寬度的塑性破壞區(qū),煤柱中出現(xiàn)彈性核區(qū);④在周圍其他采硐開采的影響作用下,煤柱兩側(cè)塑性破壞區(qū)逐漸擴(kuò)展,彈性核區(qū)中心應(yīng)力達(dá)到煤柱極限強(qiáng)度,核區(qū)應(yīng)力變化曲線進(jìn)入平臺(tái)期,此時(shí)是煤柱失穩(wěn)破壞的臨界狀態(tài),若彈性核區(qū)中心應(yīng)力稍有上升煤柱將迅速失穩(wěn)[6]。突變理論可以指導(dǎo)條帶煤柱的參數(shù)設(shè)計(jì)[7-8],即建立支承煤柱的尖點(diǎn)突變方案,導(dǎo)出極限狀態(tài)下非彈性區(qū)在煤柱中的占比,判定煤柱寬度是否合理?;谝陨舷嗤c(diǎn),突變理論也適用于螺旋鉆機(jī)回采露天礦端幫壓煤的煤柱參數(shù)設(shè)計(jì),但需要確定一種定義相鄰鉆孔間煤柱有效寬度的方法。不同于矩形斷面硐室,圓形鉆孔內(nèi)最大切向應(yīng)力發(fā)生在孔的兩幫中點(diǎn)和頂?shù)椎闹胁?。鉆孔上形成的壓力拱將拱頂上方壓力傳遞至兩側(cè)拱腳處,由于圍巖在發(fā)生破壞時(shí)首先沿著壓力拱塌落,而鉆孔拱頂本身就接近最佳硐形,所以鉆孔周圍煤巖破壞主要發(fā)生在鉆孔兩側(cè),參數(shù)設(shè)計(jì)時(shí)無需考慮鉆孔拱頂破壞問題。
2.1.1 孔間煤柱有效寬度
規(guī)則煤柱的有效寬度容易確定,但對(duì)于相鄰圓形鉆孔間這種不規(guī)則煤柱的有效寬度筆者提出一種方法:相鄰螺旋鉆孔間煤柱可以看作無梁樓蓋設(shè)計(jì)中上下端帶有斜柱帽的房柱。無梁樓蓋是一種建筑結(jié)構(gòu),其特點(diǎn)為無主梁和次梁的板柱結(jié)構(gòu)體系。柱帽具有實(shí)際作用:較寬的斜柱帽能減小頂板柱帽及跨中位置處的彎矩[9],房柱對(duì)頂板中彎矩的影響主要來自斜柱帽,中柱則負(fù)責(zé)承擔(dān)載荷,這與鉆孔周圍的應(yīng)力分布規(guī)律相吻合。由于煤柱強(qiáng)度隨煤柱最大內(nèi)切圓直徑的增加而增強(qiáng)[10-11],所以對(duì)于相鄰鉆孔間煤柱,其有效寬度應(yīng)為中柱最大內(nèi)切圓直徑。如圖1所示,其中a為煤柱的有效寬度(簡(jiǎn)稱為孔間煤柱寬度)。
圖1 相鄰鉆孔間煤柱有效寬度Fig.1 Effective width of pillar between adjacent boreholes
2.1.2 采場(chǎng)穩(wěn)定性判定標(biāo)準(zhǔn)
王東等[7]基于尖點(diǎn)突變理論得到了突變分叉集方程:
(1)
式中:Δ為系統(tǒng)分叉集控制方程;Ws為煤柱寬度;Wm為采硐寬度;Y為煤柱屈服區(qū)寬度;E為煤柱的彈性模量;H為上覆巖層的平均厚度;γ為容重;h為煤層厚度;u0為峰值載荷下煤柱的變形。
式(1)表明:當(dāng)Δ<0時(shí),系統(tǒng)會(huì)發(fā)生突變,即支承煤柱屈服區(qū)寬度2Y與煤柱寬度Ws之比大于0.88時(shí),煤柱將發(fā)生突變失穩(wěn)。根據(jù)有效支承區(qū)域理論,即煤柱支承著自身上部和間隔相鄰煤柱平分的上覆巖層重力[12]。孟達(dá)等[13]對(duì)房柱式開采采場(chǎng)穩(wěn)定性進(jìn)行了相似材料模擬試驗(yàn)并得出結(jié)論:當(dāng)采場(chǎng)中間煤柱失穩(wěn)垮落瞬間,相鄰煤柱所承載的有效支承區(qū)域載荷增加1倍,每次垮落都導(dǎo)致基本頂承載的載荷成倍增加,而且相鄰煤柱由破壞前的軸心受壓轉(zhuǎn)化為偏心受壓,逐步從小偏心受壓到大偏心受壓,從而導(dǎo)致基本頂連續(xù)破壞,產(chǎn)生連鎖反應(yīng),最終導(dǎo)致整體切冒。所以確保采場(chǎng)中點(diǎn)處煤柱的穩(wěn)定是保障鉆孔群及邊坡穩(wěn)定的重要前提。選取采場(chǎng)中點(diǎn)處煤柱為監(jiān)測(cè)煤柱(若采場(chǎng)中孔間煤柱數(shù)量為奇數(shù),則選取中點(diǎn)處煤柱;若采場(chǎng)中孔間煤柱數(shù)量為偶數(shù),則選取中點(diǎn)處2根煤柱),以監(jiān)測(cè)煤柱中屈服區(qū)寬度與煤柱有效寬度的比值是否大于0.88為判定采場(chǎng)穩(wěn)定性的標(biāo)準(zhǔn)。
ANWAR[14]研究認(rèn)為當(dāng)相鄰隧道間距較小時(shí),2個(gè)隧道上方形成1個(gè)壓力拱;當(dāng)間距較大時(shí),各自形成獨(dú)立的壓力拱,如圖2所示。隨著相鄰采空區(qū)數(shù)量的增多,原本各自獨(dú)立的壓力拱形成1個(gè)大壓力拱,拱腳前移,拱頂高度增加,圍巖內(nèi)應(yīng)力集中程度增加。當(dāng)應(yīng)力集中程度大于頂板巖層的承載能力,上覆巖層沿壓力拱輪廓發(fā)生破壞。同時(shí),根據(jù)巖層控制的關(guān)鍵層理論[15],主關(guān)鍵層對(duì)直至地表的巖層活動(dòng)起控制作用。所以同一水平上相鄰鉆孔數(shù)量不能過多。
圖2 相鄰隧道壓力拱Fig.2 Adjacent tunnel pressure arch
螺旋鉆機(jī)開采邊坡煤鉆孔排布參數(shù)可以看作多孔之間相互影響問題??字苓叺膽?yīng)力集中系數(shù)K可以表示為
K=σ′/σ
(2)
其中:σ′為孔間煤柱內(nèi)部應(yīng)力峰值;σ為原巖應(yīng)力。應(yīng)力集中系數(shù)一方面隨著D與B之比增大而增大(D為孔徑,B為孔周邊的間距),另一方面又受同一水平上孔的數(shù)量影響。顯然,孔的數(shù)量越多,孔周邊的應(yīng)力集中系數(shù)越大[15]。
而合理寬度的永久隔離煤柱可以有效地減小相鄰鉆孔群之間的相互影響。當(dāng)鉆孔上獨(dú)立的壓力拱形成一個(gè)更大的壓力拱時(shí),永久隔離煤柱實(shí)際上為壓力拱的拱腳處,拱腳處的應(yīng)力狀態(tài)受壓力拱傳遞的不穩(wěn)定載荷的集中程度影響。在埋深、鉆孔直徑不變的情況下,鉆孔數(shù)量愈多、相鄰鉆孔間煤柱寬度愈小,拱腳處的應(yīng)力集中系數(shù)愈大。隔離煤柱增加了相鄰鉆孔群間距,有效地分隔了相鄰工作面垮落空間的橫向貫通和縱向擴(kuò)展,隔開了相鄰工作面上覆巖層沉降曲線[16]。定義孔周邊應(yīng)力集中系數(shù)K在永久隔離煤柱中點(diǎn)處等于1時(shí),永久隔離煤柱寬度為最優(yōu)值。
某露天礦開采后,在邊坡下有大量的壓煤,應(yīng)用其他采煤方法很難將這部分煤采出,故采用螺旋鉆機(jī)將部分壓煤采出。煤層坡面角60°,南北走向傾角較小,可認(rèn)為是近水平煤層。邊坡上覆巖層最小厚度20 m,最大厚度130 m。煤層平均厚度5 m,頂板為粉砂巖、細(xì)砂巖,底板以泥巖為主,煤層全區(qū)發(fā)育,煤層結(jié)構(gòu)簡(jiǎn)單。螺旋鉆機(jī)對(duì)端幫煤層露頭進(jìn)行單水平順序鉆孔回采,鉆頭最大直徑1.2 m,鉆孔與鉆孔之間留設(shè)隔離煤柱,采出一系列鉆孔后留設(shè)一段永久隔離煤柱。
壓力拱高度與關(guān)鍵層高度有如下關(guān)系:壓力拱發(fā)育至關(guān)鍵層后若關(guān)鍵層破斷則基于破斷位置重新起拱,若不破斷則拱頂位于關(guān)鍵層附近[17]。所以,同一水平相鄰鉆孔數(shù)量應(yīng)保證在壓力拱發(fā)育至關(guān)鍵層后,拱內(nèi)圍巖應(yīng)力集中小于頂板圍巖的承載能力,避免拱內(nèi)圍巖破壞導(dǎo)致關(guān)鍵層下圍巖離層引起關(guān)鍵層破斷。
“正交設(shè)計(jì)法”是研究與處理多因素試驗(yàn)的一種科學(xué)的試驗(yàn)方法。經(jīng)分析選取相鄰鉆孔間煤柱寬度、相鄰鉆孔數(shù)量、永久隔離煤柱寬度為3個(gè)主要因素,每個(gè)因素取3個(gè)水平,共設(shè)置9次數(shù)值模擬試驗(yàn),繪制L9(34)正交表(表1)。鉆孔埋深為固定值,不作為正交試驗(yàn)中的因素,試驗(yàn)方案見表2。
表1 單層鉆孔排布各因素水平Table 1 Single layer drilling arrangement of each factor level
表2 試驗(yàn)方案Table 2 Experimental scheme
數(shù)值模擬計(jì)算主要借助有限差分軟件。考慮到鉆孔長(zhǎng)度可達(dá)60 m,并可能進(jìn)一步增加,且在數(shù)值計(jì)算的過程中常把該問題視為平面問題求解,所以在數(shù)值模擬的過程中將模型簡(jiǎn)化為三維切片。計(jì)算應(yīng)用摩爾-庫(kù)侖彈塑性模型,結(jié)合采場(chǎng)的實(shí)際情況,模型位置選取距鉆孔底部1 m處厚度1 m的切片。根據(jù)圣維南原理,局部效應(yīng)的影響可以由足夠大的圍巖范圍消除,取模型高度為煤層厚度的6~7倍。最終模型的尺寸為:120 m×1 m×98 m,共1 623 400個(gè)單元。左右和下邊界施加位移約束,上邊界沒有施加約束。模型與應(yīng)力監(jiān)測(cè)路徑如圖3所示。
圖3 數(shù)值模擬模型及應(yīng)力監(jiān)測(cè)路徑Fig.3 Numerical simulation model and stress monitoring path
通過對(duì)取自采場(chǎng)的煤巖樣進(jìn)行單軸抗壓、抗拉強(qiáng)度試驗(yàn)以及單軸壓縮變形試驗(yàn)并經(jīng)胡克布朗準(zhǔn)則折減后獲得了模擬所需要的采礦巖石力學(xué)參數(shù),見表3。
表3 煤巖樣物理力學(xué)試驗(yàn)結(jié)果Table 3 Results of coal physical and mechanical test
由于模擬結(jié)果數(shù)據(jù)較多,無法詳盡列出,僅選取9組方案中拱頂路徑上的原巖應(yīng)力與開挖后水平方向和垂直方向應(yīng)力變化曲線,如圖4、圖5所示,拱腰路徑上應(yīng)力集中系數(shù)見表4。
表4 拱腰路徑上應(yīng)力集中系數(shù)Table 4 Coefficient of stress concentration on arch waist path
圖4 拱頂路徑1、2、3上水平應(yīng)力變化Fig.4 Horizontal stress variation on vault paths 1,2 and 3
圖5 拱頂路徑1、2、3上垂直應(yīng)力變化Fig.5 Vertical stress variation on vault path 1,2 and 3
隨著開采進(jìn)程的推進(jìn),鉆孔開采對(duì)拱頂監(jiān)測(cè)路徑上應(yīng)力的影響主要反映在:路徑1、路徑3上原本隨埋深增加而增加的水平應(yīng)力與豎直應(yīng)力,即拱頂監(jiān)測(cè)路徑上的切向應(yīng)力和徑向應(yīng)力均出現(xiàn)拐點(diǎn)。具體現(xiàn)象為:拱頂監(jiān)測(cè)路徑1、路徑3上,鉆孔周圍的切向應(yīng)力在開挖后沿拱頂監(jiān)測(cè)路徑方向先增大,出現(xiàn)切向應(yīng)力壓緊區(qū),后下降并回歸至原巖應(yīng)力水平;沿拱頂監(jiān)測(cè)路徑方向,徑向應(yīng)力先增大后減小,在各條拱頂監(jiān)測(cè)路徑上,除方案8、方案9外,方案1—方案7中徑向應(yīng)力始終低于原巖應(yīng)力水平,說明圍巖中出現(xiàn)松動(dòng)區(qū)域。
拱腰應(yīng)力監(jiān)測(cè)路徑上的應(yīng)力集中系數(shù)表明,保持鉆孔數(shù)量不變,孔間煤柱寬度的增加會(huì)降低拱腳處的應(yīng)力集中程度;保持孔間煤柱寬度不變,相鄰鉆孔數(shù)量的增加會(huì)提高拱腳處應(yīng)力集中程度;永久隔離煤柱能有效地減小拱腳處應(yīng)力集中程度和相鄰鉆孔群之間的影響。
1)對(duì)于拱頂內(nèi)邊界和外邊界的確定主要選取如下2個(gè)準(zhǔn)則:①根據(jù)切向應(yīng)力判定準(zhǔn)則,將切向應(yīng)力σxx由最大值減小到原巖應(yīng)力水平的點(diǎn)作為壓力拱的外邊界。該準(zhǔn)則反映出鉆孔的開挖對(duì)周圍煤巖的影響范圍。②根據(jù)拱體內(nèi)部主應(yīng)力極值判定準(zhǔn)則:鉆孔開挖后切向應(yīng)力由低于原巖應(yīng)力增加到最大值后又逐漸恢復(fù)到原巖應(yīng)力水平,將切向應(yīng)力增加到最大值的點(diǎn)作為壓力拱的內(nèi)邊界。
由圖5可知,在拱頂應(yīng)力監(jiān)測(cè)路徑2上,除方案7外,剩余8組方案中的切向應(yīng)力均在距鉆孔拱頂9 m處增加到最大值,在14 m處由最大值回歸至原巖應(yīng)力水平,即壓力拱內(nèi)邊界高度9 m,外邊界高度14 m;而方案7中切向應(yīng)力最大值位于距鉆孔拱頂18 m處,距鉆孔拱頂25 m處仍未回歸至原巖應(yīng)力水平。在監(jiān)測(cè)路徑1和路徑3上,方案1—方案6中切向應(yīng)力在7 m處增加到最大值,在9 m處回歸至原巖應(yīng)力水平,即方案1—方案6壓力拱內(nèi)邊界高度7 m,外邊界高度9 m;方案8、方案9中切向應(yīng)力距鉆孔拱頂6 m處達(dá)到最大值,在14 m附近回歸至原巖應(yīng)力水平,即壓力拱內(nèi)邊界6 m,外邊界14 m。壓力拱將上覆巖體載荷轉(zhuǎn)化為拱的切向應(yīng)力,并將其傳遞至兩側(cè)永久隔離煤柱,保證回采作業(yè)處于應(yīng)力降低區(qū)。根據(jù)采場(chǎng)中巖層賦存條件,距鉆孔拱頂9~25 m處的粗砂巖層符合關(guān)鍵層標(biāo)準(zhǔn)。除方案7外的所有方案中壓力拱外邊界均位于關(guān)鍵層附近。由圖5可知,在監(jiān)測(cè)路徑1和路徑3上,僅有方案8、方案9中徑向應(yīng)力在關(guān)鍵層附近回歸至原巖應(yīng)力水平,表明未出現(xiàn)松動(dòng)區(qū)。在監(jiān)測(cè)路徑2上,方案9中徑向應(yīng)力與原巖應(yīng)力差值最小,可知方案9中關(guān)鍵層附近松動(dòng)區(qū)域最小,關(guān)鍵層附近圍巖較穩(wěn)定。
2)9組方案的塑性區(qū)分布如圖6所示,由于有限差分軟件具有動(dòng)態(tài)松弛特性,只有處于平衡狀態(tài)的單元才有意義,所以將shear-p設(shè)為空單元。方案1、方案4、方案7中監(jiān)測(cè)煤柱內(nèi)塑性區(qū)寬度與煤柱的有效寬度比值為93%~100%,遠(yuǎn)大于臨界值88%,表明煤柱內(nèi)塑性區(qū)已貫通,將發(fā)生突變失穩(wěn)。方案2、方案5、方案8中監(jiān)測(cè)煤柱內(nèi)塑性區(qū)寬度與其有效支承寬度比值為55%~60%,未達(dá)到臨界值,煤柱保持穩(wěn)定。方案3、方案6、方案9中監(jiān)測(cè)煤柱兩幫處存在小范圍剪切破裂區(qū),塑性區(qū)寬度與煤柱有效寬度比值為36%~40%,遠(yuǎn)小于臨界值,煤柱保持穩(wěn)定。
3)對(duì)于拱腰處的應(yīng)力監(jiān)測(cè)路徑,徑向應(yīng)力σ′zz=σxx,切向應(yīng)力σ′xx=σzz。拱腰路徑上的壓力拱內(nèi)邊界同樣適合內(nèi)部主應(yīng)力極值確定準(zhǔn)則,但不對(duì)其進(jìn)行研究。拱腰路徑上壓力拱的外邊界可以用該路徑上的切向應(yīng)力集中系數(shù)確定開挖對(duì)周圍煤巖內(nèi)應(yīng)力大小的影響范圍,故選取拱腰路徑上切向應(yīng)力集中系數(shù)等于1處為拱腰路徑上壓力拱的外邊界。
由表4可知,在鉆孔直徑不變的情況下,孔間煤柱寬度的減小和鉆孔數(shù)量的增加都會(huì)使拱腳應(yīng)力集中系數(shù)增大。在拱腰應(yīng)力監(jiān)測(cè)路徑上,切向應(yīng)力自鉆孔群兩側(cè)至永久隔離煤柱中點(diǎn)方向先增大后減小,但始終大于原巖應(yīng)力,說明在拱腰路徑上未出現(xiàn)松動(dòng)區(qū)域。除方案7、方案8外,剩余方案中沿拱腰應(yīng)力監(jiān)測(cè)路徑上的切向應(yīng)力在永久隔離煤柱中點(diǎn)處基本回歸至原巖應(yīng)力水平,表明這7組方案中的永久隔離煤柱能有效地消除相鄰鉆孔的影響。
1)將壓力拱理論、正交試驗(yàn)法與數(shù)值模擬相結(jié)合,基于壓力拱理論分析正交試驗(yàn)所得方案在開采后的數(shù)值模擬結(jié)果,優(yōu)選出合理單層鉆孔排布方案,從而使參數(shù)設(shè)計(jì)更科學(xué),降低開采作業(yè)風(fēng)險(xiǎn)。
2)結(jié)合建筑學(xué)中無梁樓蓋設(shè)計(jì)理論,提出一種定義相鄰圓形鉆孔間煤柱有效寬度的方法,即將煤柱分解成由斜柱帽與中柱組成的房柱,該方法適合螺旋鉆機(jī)開采鉆孔參數(shù)設(shè)計(jì)研究。
3)經(jīng)分析,不同的鉆孔排布參數(shù)對(duì)開采后圍巖內(nèi)的應(yīng)力分布有不同程度的影響;相鄰鉆孔數(shù)量25個(gè),相鄰鉆孔間煤柱寬度1 m,永久隔離煤柱寬度10 m,關(guān)鍵層附近松動(dòng)區(qū)域最小、監(jiān)測(cè)煤柱穩(wěn)定、永久隔離煤柱寬度合理,優(yōu)于其他方案。但缺少實(shí)際開挖所得數(shù)據(jù),未對(duì)最優(yōu)方案的合理性進(jìn)行科學(xué)驗(yàn)證。