考慮煤層傾角和鉆孔傾角的鉆屑量理論研究

周 洋，梁 冰，石占山，孫維吉，劉 學(xué)

(1.遼寧工程技術(shù)大學(xué) 礦業(yè)學(xué)院，遼寧 阜新 123000；2.遼寧工程技術(shù)大學(xué) 力學(xué)與工程學(xué)院，遼寧 阜新123000；3.煤炭科學(xué)技術(shù)研究院有限公司 安全分院，北京 100013)

煤與瓦斯突出是煤礦中發(fā)生的極為復(fù)雜的煤與瓦斯涌出現(xiàn)象,是嚴(yán)重威脅煤礦安全生產(chǎn)的主要災(zāi)害之一[1-3]。國內(nèi)外廣大科技工作者雖已對煤與瓦斯突出問題進(jìn)行了較深入的認(rèn)識和研究,但深部開采情況下地質(zhì)條件復(fù)雜、受高地應(yīng)力及其他各種因素影響,煤與瓦斯突出災(zāi)害問題還未根本解決。由于鉆屑量是煤體應(yīng)力的綜合反映，是煤礦煤與瓦斯突出預(yù)測和預(yù)防的主要指標(biāo)之一，因此,研究影響鉆屑量及相關(guān)規(guī)律的因素對煤與瓦斯突出的預(yù)測和預(yù)防具有重要意義[4-7]。

對于鉆屑法的研究,國內(nèi)外學(xué)者很早便開展了大量的理論和實驗研究工作。勃羅納[8]將鉆屑排放分為兩部分,其中一部分與圓柱形煤體的鉆孔直徑相同,另一部分是鉆孔內(nèi)壁的位移產(chǎn)生量；王佑安[9]通過Mohr-Coulomb準(zhǔn)則分析了利用鉆屑量如何反映煤體應(yīng)力的問題；潘一山[10]通過考慮膨脹因子等因素,得到鉆屑量與煤體應(yīng)力間的關(guān)系；王振等[11]提出了應(yīng)力梯度下鉆屑量的增量變化指數(shù)；曲效成等[12]確定了鉆屑量和鉆孔應(yīng)力之間的聯(lián)系;董杰等[13]分析了煤巖物理和力學(xué)性質(zhì),壓力和測量誤差對鉆屑量的影響；陳峰等[14]研究了最大鉆屑量對沖擊地壓的影響；王志明等[15]指出圍巖的應(yīng)力狀態(tài)可以通過鉆屑法確定。

我國的大部分煤層都有傾角，隨著開采深度的增加,水平地應(yīng)力的作用也越大,同時在打鉆時鉆孔傾角也會對水平地應(yīng)力的計算產(chǎn)生一定影響。因此在考慮煤層傾角與鉆孔傾角因素下，通過對有效地應(yīng)力的計算修正,從而建立新理論公式并對鉆屑量與煤層傾角和鉆孔傾角之間的關(guān)系進(jìn)行分析，根據(jù)平煤有限公司八礦的鉆屑量實測值，對鉆屑量的新理論公式進(jìn)行了對比驗證,為準(zhǔn)確預(yù)測地應(yīng)力提供參考。

1 考慮傾角的鉆屑量公式研究

從理論分析角度，深部開采條件下，鉆孔傾角與煤層傾角對水平地應(yīng)力的作用不能忽略，引入鉆孔傾角與煤層傾角推導(dǎo)得到新鉆屑量公式。

1.1 模型的建立

采用Mohr-Coulomb準(zhǔn)則進(jìn)行分析，建立了鉆孔圍巖應(yīng)力分布力學(xué)模型，如圖1所示。其中，σv為垂直地應(yīng)力；σh為最小水平地應(yīng)力；r為鉆頭半徑；R0為破碎圈半徑；R為非彈性變形區(qū)半徑。

Ⅰ—破碎區(qū)；Ⅱ—塑性軟化區(qū)；Ⅲ—彈性區(qū)。

1.2 新鉆屑量公式推導(dǎo)

鉆孔鉆屑量主要分為靜態(tài)部分和動態(tài)部分鉆屑量。靜態(tài)鉆屑量即鉆孔實芯部分G1；動態(tài)鉆屑量由鉆孔彈性變形區(qū)部分G2，破碎區(qū)與彈性區(qū)交界處部分G3，煤體擴(kuò)容部分G4組成。總鉆屑量為：

G=G1+G2+G3+G4

(1)

鉆孔煤體實芯部分G1表示為：

G1=πr2ρ0

(2)

式中：r為鉆頭半徑，mm；ρ0為原始煤體密度，g/cm3。

鉆孔彈性變形部分G2表示為：

(3)

式中：μ為煤體泊松比；E為煤體彈性模量，MPa；σe為有效應(yīng)力，MPa。

采用易俊等[16]試驗推導(dǎo)的有效應(yīng)力σe計算公式:

(4)

式中:σv[sin2(90°-α)+Kcos2(90°-α)]-1為引用汪大海等[17]推導(dǎo)的正應(yīng)力計算公式；σH為最大水平地應(yīng)力，MPa；α為鉆孔傾角，(°)；K為煤層側(cè)壓系數(shù);p1為鉆孔非彈性區(qū)瓦斯壓力，MPa；p2為鉆孔形成后周圍的瓦斯壓力，取0.1 MPa。

煤層的水平應(yīng)力主要由上覆巖層壓力、構(gòu)造運(yùn)動的作用引起，而煤層傾角和鉆孔傾角的存在也影響了水平地應(yīng)力與垂直地應(yīng)力的計算,筆者引用樓一珊[18]推導(dǎo)的煤層側(cè)壓系數(shù)計算公式:

K=KαKβ

(5)

式中:Kα為鉆孔傾角影響系數(shù),Kα=(α+β)/(α-β)；Kβ為煤層傾角影響系數(shù),Kβ=β/(α-β)；β為煤層傾角，(°)。

考慮煤層傾角與鉆孔傾角的水平地應(yīng)力計算模型為:

(6)

G3為破碎區(qū)與彈性區(qū)交界處因彈性卸載產(chǎn)生的附加鉆屑量：

(7)

(8)

式中：K1為破碎區(qū)煤體的三軸殘余強(qiáng)度系數(shù)；R0為破碎圈半徑，mm；σ0為煤體殘余強(qiáng)度，MPa。

G4為鉆孔形成后，煤體擴(kuò)容鉆屑量。煤體應(yīng)變軟化本構(gòu)方程為：

(9)

式中：m為煤體塑性軟化系數(shù)；σc為單軸抗壓強(qiáng)度，MPa；εc為單軸抗壓強(qiáng)度相應(yīng)的應(yīng)變。

依據(jù) Mohr-Coulomb 準(zhǔn)則，建立彈塑性平衡方程為：

(10)

式中：σθ為徑向應(yīng)力，MPa；σr為切向應(yīng)力，MPa。

屈服條件為：

σθ=qσr+σc

(11)

邊界條件：當(dāng)r=R0時，σr=0；當(dāng)r=R時，得出非彈性變形區(qū)半徑R的解析式為：

(12)

當(dāng)m=0時，式(12)即為理想Kastner公式。在不考慮擴(kuò)容情況時，鉆孔非彈性部分與彈性部分交界處的徑向位移uR為：

Kcos2(90°-α))-1-σc]}

(13)

孔內(nèi)壁徑向位移由體積不變條件得：

(14)

考慮煤體非彈性變形擴(kuò)容時的孔內(nèi)壁徑向位移ur為：

(15)

綜上可得，擴(kuò)容單位孔深鉆屑量G4推導(dǎo)為：

(16)

由以上推導(dǎo)可得新鉆屑量G為：

(17)

式(17)即為考慮鉆孔傾角和煤層傾角情況下新鉆屑量理論公式。

1.3 新鉆屑量公式影響因素分析

從式(17)可以看出，煤層傾角和鉆孔傾角對鉆屑量有一定的影響。以典型的高瓦斯突出礦井——平煤股份有限公司八礦為例，取屑鉆孔直徑為 94 mm，通過對新鉆屑量公式進(jìn)行分析，并將表1中的物理力學(xué)參數(shù)代入到式(17)中，得出在考慮鉆孔傾角與煤層傾角影響下的新鉆屑量變化規(guī)律。

我國煤礦井下條件較為復(fù)雜，煤層的賦存條件不同[21]，大多數(shù)煤礦煤層不同程度地存在著煤層傾角,并且隨著開采深度的加大，由于重力勢能作用，對地應(yīng)力的影響越大，而打鉆時鉆孔傾角也影響著鉆孔受到的垂直應(yīng)力和水平應(yīng)力的大小，在兩者共同作用下鉆屑量的大小必然也會受到影響。由式(4)～(6)可知，在考慮煤層傾角與鉆孔傾角因素下與鉆屑量的關(guān)系如圖2～4所示。

圖2 鉆屑量與煤層傾角和鉆孔傾角間的關(guān)系

圖3 煤層傾角0°時鉆屑量與鉆孔傾角變化關(guān)系

圖4 鉆屑量減小比率隨鉆孔傾角的變化關(guān)系

從圖2～4可知，鉆屑量與煤層傾角和鉆孔傾角均呈冪指數(shù)遞減關(guān)系，且遞減幅度先大再逐漸減小后再次逐漸增大。煤層傾角與鉆孔傾角對鉆屑量的計算有較大影響,當(dāng)煤層傾角較小時影響較小,煤層傾角較大時影響較大；當(dāng)鉆孔傾角較小時影響較小,鉆孔傾角較大時影響較大。主要原因是大多數(shù)煤層是各向異性體或橫觀各向同性體,其各方向上的彈性參數(shù)是不相同的,當(dāng)煤層具有一定傾角時,這種差別更為明顯。而鉆孔具有一定傾角時，所受水平地應(yīng)力與垂直地應(yīng)力均會不同程度的變化，變化幅度隨傾角增大而增大。因此，對于大傾角煤層、深部開采高地應(yīng)力、高瓦斯壓力礦井，考慮煤層傾角與鉆孔傾角因素計算鉆屑量預(yù)測煤與瓦斯突出是十分有必要的。

2 現(xiàn)場對比驗證

在平煤股份有限公司八礦己15-21030風(fēng)巷進(jìn)行鉆屑量實測研究，對比驗證新鉆屑量理論公式準(zhǔn)確性。平煤股份有限公司八礦己15-21030風(fēng)巷現(xiàn)場基本情況如下：地面標(biāo)高+75～+90 m，己15煤層區(qū)段標(biāo)高-728～-840 m，埋深803～930 m；工作面按綜采工作面布置，采高3.5 m，平均采長205 m，可采走向長1 568 m，儲量147萬t；煤層厚度一般在3.0～4.4 m內(nèi)，平均厚度為3.5 m，煤層傾角一般為5°～15°；該工作面正常涌水量為2～3 m3/h，最大涌水量為 5 m3/h，屬突出危險工作面。己15-21030風(fēng)巷鉆孔鉆屑量G(采用稱重法)測試：鉆孔采用ZYWL-4000SY煤礦用雙履帶式全液壓鉆機(jī)鉆進(jìn)，搭配直徑94 mm螺旋鉆桿，煤層實際傾角為14°，鉆孔開孔高度1.2 m，每個鉆孔實際鉆進(jìn)15 m，取前10 m鉆屑量進(jìn)行測試，現(xiàn)場在同一區(qū)域排除其他影響因素的情況下共計測試5個鉆孔，鉆孔間距2 m。鉆屑量現(xiàn)場實測值如表 2 所示。

表2 鉆屑量實測值

將表1中平煤八礦的煤體力學(xué)參數(shù)代入到新鉆屑量公式中，計算出鉆屑量理論值，并將其與現(xiàn)場實測值進(jìn)行誤差對比分析，分析結(jié)果如表3所示。

表3 鉆屑量理論值與實測值誤差分析

從表3中可知，原鉆屑量理論值誤差隨著傾角的變小逐漸減小，而在考慮了煤層傾角與鉆孔傾角影響下的新鉆屑量理論值誤差波動較小，為3.9%。

3 結(jié)論

1)在煤層傾角一定時，鉆屑量隨鉆孔傾角增大呈冪指數(shù)遞減關(guān)系，減小的比率呈“V”形變化趨勢；在鉆孔傾角一定時，鉆屑量隨煤層傾角增大呈冪指數(shù)遞減關(guān)系，減小的比率也呈“V”形變化趨勢。

2)工程實例中考慮煤層傾角與鉆孔傾角影響因素，計算得到的鉆屑量與現(xiàn)場實測值平均誤差為3.9%，比原鉆屑量理論值平均誤差(5.5%)減小1.6%。說明在考慮煤層傾角與鉆孔傾角影響因素下建立的新鉆屑量公式是有效合理的，更適用于礦山深部開采與現(xiàn)場實際情況，能為預(yù)測與防治煤與瓦斯突出提供參考。