龐立寧，付書俊，蘇 波

（1.天地科技股份有限公司 開采設(shè)計(jì)事業(yè)部，北京 100013；2.中煤科工開采研究院有限公司，北京 100013；3.華陽(yáng)新材料科技集團(tuán)有限公司 智能礦山事業(yè)部，山西 陽(yáng)泉 045000）

煤層大直徑鉆孔和頂板預(yù)裂孔卸壓技術(shù)是沖擊地壓礦井最常用的防沖手段之一[1-3]。國(guó)內(nèi)外學(xué)者對(duì)卸壓鉆孔防沖機(jī)理進(jìn)行過不少研究，如賈傳洋[4]認(rèn)為卸壓鉆孔防沖機(jī)制在于鉆孔周邊煤巖體裂紋擴(kuò)展貫通導(dǎo)致應(yīng)力釋放；李金奎[5]等通過數(shù)值模擬發(fā)現(xiàn)大直徑卸壓鉆孔使得煤體應(yīng)力峰值向深部轉(zhuǎn)移從而達(dá)到卸壓效果；蘭永偉[6]通過數(shù)值模擬得出鉆孔卸壓作用根本原因是其周圍煤體的弱化；吳鑫等[7]通過3DEC數(shù)值模擬手段得出鉆孔卸壓使得巷幫應(yīng)力峰值向圍巖深部轉(zhuǎn)移，且在鉆孔末端形成應(yīng)力集中；閆永敢[8]研究發(fā)現(xiàn)，鉆孔周圍非彈性區(qū)范圍大小與卸壓效果成正比。但大多都從應(yīng)力轉(zhuǎn)移等卸壓機(jī)制方面研究，從改變巷道圍巖結(jié)構(gòu)、增大動(dòng)載荷應(yīng)力波傳遞耗散角度研究卸壓鉆孔防沖機(jī)理的文獻(xiàn)較少，為此，結(jié)合葫蘆素煤礦現(xiàn)場(chǎng)實(shí)際，進(jìn)行煤層大直徑鉆孔和頂板預(yù)裂鉆孔防沖機(jī)理研究，制定煤層大直徑鉆孔和頂板壓裂孔卸壓技術(shù)方案，在21103工作面回風(fēng)巷進(jìn)行了工業(yè)性試驗(yàn)，取得了好的效果。

1 煤層大直徑鉆孔防沖機(jī)理

1.1 巷道圍巖力學(xué)模型

巷道開挖前，煤層上覆巖層重力均勻分布在煤體上，即準(zhǔn)巷道圍巖均處在原巖應(yīng)力狀態(tài)。

煤層巷道開挖后，巷道圍巖應(yīng)力重新分布，根據(jù)切應(yīng)力大小，可分為減壓區(qū)、增壓區(qū)、穩(wěn)壓區(qū)[9]，另外，以支承壓力峰值為界，可分為極限平衡區(qū)和彈性區(qū)。巷道圍巖不同位置應(yīng)力不同，但滿足下式：

式中：σx、εx分別為距離巷幫x位置處的垂直應(yīng)力和應(yīng)變；E為彈性模量；ρ為上覆巖層的密度；H為埋深；L為距巷幫的距離。

不同深度巷道圍巖的力學(xué)性能不同，可由不同特性的力學(xué)模型來(lái)模擬。根據(jù)巖石應(yīng)力應(yīng)變?nèi)糖€，巖石峰后區(qū)破壞階段呈現(xiàn)出應(yīng)變軟化特性[10]，即隨著應(yīng)力增大，應(yīng)變?cè)鲩L(zhǎng)的速率加快，這與巷道圍巖極限平衡區(qū)巖體力學(xué)性能較類似，可用彈黏體元件串聯(lián)模擬。巷道圍巖彈性區(qū)可用完全彈性體元件模擬。因此，巷道圍巖力學(xué)模型如圖1。

圖1 巷道開挖后圍巖力學(xué)模型Fig.1 Mechanical model of surrounding rock after roadway excavation

1.2 大直徑鉆孔對(duì)巷道圍巖結(jié)構(gòu)影響

幫部大直徑卸壓鉆孔的實(shí)施為頂板巖層的下沉或變形提供了空間，同時(shí)幫部大直徑卸壓區(qū)的存在也改變了巷道圍巖的力學(xué)結(jié)構(gòu)模型。當(dāng)支承壓力峰值以外人為產(chǎn)生卸壓區(qū)后，相當(dāng)于將該區(qū)域完全彈性體元件結(jié)構(gòu)模型變?yōu)閺楌んw串聯(lián)元件模型，大直徑卸壓鉆孔實(shí)施后巷道圍巖力學(xué)模型如圖2。

圖2 大直徑卸壓鉆孔實(shí)施后巷道圍巖力學(xué)模型Fig.2 Mechanical model of roadway surrounding rock after large diameter pressure relief drilling

大直徑卸壓鉆孔的存在，使得頂板下沉或變形主要以壓縮卸壓孔空間（黏體元件）方式實(shí)現(xiàn)，此時(shí)煤體內(nèi)產(chǎn)生的應(yīng)變大都被“黏體元件”吸收，因此不會(huì)導(dǎo)致應(yīng)力的增大。當(dāng)卸壓孔空間被壓實(shí)（可理解為黏體元件中阻尼器行程已滿）以后，彈黏體串聯(lián)元件中的彈性元件開始發(fā)揮作用，隨著卸壓孔圍巖內(nèi)彈性形變的增加（可理解為完全彈性體彈簧行程增加），垂直應(yīng)力慢慢恢復(fù)，這也是大直徑卸壓鉆孔存在時(shí)效性的原因。

2 頂板預(yù)裂鉆孔防沖機(jī)理

2.1 應(yīng)力波衰減機(jī)制

沖擊地壓?jiǎn)?dòng)理論中煤巖系統(tǒng)集中靜載荷是內(nèi)因，系統(tǒng)外集中動(dòng)載荷是外因，主要起到誘發(fā)沖擊地壓的作用[11]。集中動(dòng)載荷的產(chǎn)生主要來(lái)源于采場(chǎng)周圍上覆頂板巖層的斷裂或者采掘爆破等人為擾動(dòng)[12-13]，其通過應(yīng)力波形式作用于集中靜載荷區(qū)。

采場(chǎng)頂板斷裂等產(chǎn)生的應(yīng)力波在煤巖介質(zhì)中傳播時(shí)會(huì)發(fā)生衰減，其衰減主導(dǎo)機(jī)制[14]主要包含摩擦效應(yīng)（基質(zhì)裂隙相對(duì)運(yùn)動(dòng)時(shí)接觸面摩擦生熱消耗應(yīng)力波攜帶的彈性能）、散射效應(yīng)（組成煤巖體介質(zhì)的晶粒對(duì)應(yīng)力波散射）、流動(dòng)效應(yīng)（孔隙介質(zhì)骨架與流體之間相對(duì)運(yùn)動(dòng)引發(fā)內(nèi)摩擦耗能）、裂紋擴(kuò)展效應(yīng)（應(yīng)力波攜帶能量轉(zhuǎn)化為裂紋擴(kuò)展表面能）。

巖石物理學(xué)中Hudson模型描述煤體中摩擦效應(yīng)導(dǎo)致的應(yīng)力波衰減系數(shù)α，可以用下式表達(dá)[15-16]：

式中：V為煤巖體傳播介質(zhì)總體積（數(shù)值與應(yīng)力波產(chǎn)生點(diǎn)和沖擊啟動(dòng)點(diǎn)距離正相關(guān)）；l、N分別為裂隙長(zhǎng)度和數(shù)量；cd、cs分別為煤巖體系統(tǒng)中縱波和橫波波速；A、B分別為cd/cs的5次方函數(shù)；U1、U3的取值取決于煤體含水情況；ω為角頻率。

2.2 頂板預(yù)裂孔防沖機(jī)理

通過應(yīng)力波衰減表達(dá)式可以看出，煤巖體系統(tǒng)中裂隙長(zhǎng)度和數(shù)量與應(yīng)力波衰減系數(shù)正相關(guān)。煤巖體中大直徑卸壓鉆孔的施工不僅可以轉(zhuǎn)移集中靜載荷，還可以增大卸壓區(qū)附近裂隙發(fā)育情況，在一定程度上衰減了頂板斷裂或采掘爆破等產(chǎn)生的應(yīng)力波。通過施工頂板孔切縫、壓裂或爆破，可以在頂板內(nèi)人為制造很多裂縫，這些裂縫的存在耗散了很大一部分應(yīng)力波攜帶的能量。

在集中動(dòng)載荷產(chǎn)生位置和靜載荷集中區(qū)距離一定的條件下，應(yīng)力波的作用效果與其持續(xù)時(shí)間和應(yīng)力波幅值呈正相關(guān)，而其持續(xù)時(shí)間正比于頂板斷裂尺度，這也是采場(chǎng)堅(jiān)硬厚頂板條件比松軟薄頂板條件更容易發(fā)生沖擊地壓的原因；其應(yīng)力波幅值正比于頂板或者煤體等傳播介質(zhì)的強(qiáng)度，這也是采場(chǎng)頂板或煤體強(qiáng)度越高越容易發(fā)生沖擊地壓的原因。頂板壓裂孔可以原位改變頂板巖層的物理力學(xué)性質(zhì)，降低頂板巖層的完整性和強(qiáng)度，減小沖擊地壓危險(xiǎn)。

頂板預(yù)裂孔除了可以衰減頂板破斷產(chǎn)生的動(dòng)載荷應(yīng)力波外，還可以切斷采空區(qū)側(cè)向懸頂，減少懸露頂板作用在煤柱上的靜荷載；這2個(gè)方面均是頂板預(yù)裂孔可防治沖擊地壓的原因。

3 現(xiàn)場(chǎng)案例

3.1 工程概況

21103工作面為葫蘆素煤礦首采區(qū)第2個(gè)回采工作面，煤層厚度平均2.54 m，鑒定結(jié)果表明，葫蘆素煤礦2-1煤具有強(qiáng)沖擊傾向性，頂板具有弱沖擊傾向性，底板無(wú)沖擊傾向性。頂板巖性多以粉砂巖、砂質(zhì)泥巖為主，底板巖性多為砂質(zhì)泥巖及粉砂巖。煤層頂?shù)装逄卣饕姳?。

表1 煤層頂?shù)装逄卣鞅鞹able 1 Characteristics of coal seam roof and floor

21103工作面回采過程中發(fā)生過4次沖擊地壓顯現(xiàn)，分別于2017年11月18日、2018年2月6日、2018年2月15日、2018年3月8日。其中2017年11月18日沖擊顯現(xiàn)最為嚴(yán)重，“11.18”沖擊位置如圖3。21103回風(fēng)巷沖擊顯現(xiàn)長(zhǎng)度60 m，沖擊瞬間造成煤柱側(cè)幫鼓起，最大變形為0.7 m，導(dǎo)致機(jī)尾無(wú)安全出口，電纜落地，共破壞錨桿30根，錨索15根，損壞單體10根。

圖3“11.18”沖擊位置Fig.3“11.18”rock burst position

3.2 21103工作面回風(fēng)巷沖擊地壓發(fā)生機(jī)理

21103工作面回風(fēng)巷圍巖不同位置靜載荷集中可以用函數(shù)J（σx，εx）表示，不同位置動(dòng)載荷集中可以用函數(shù)D（σx，εx）表示，構(gòu)建一個(gè)范函數(shù)F：

式中：a為靜載荷集中對(duì)沖擊地壓?jiǎn)?dòng)的影響因子；b為動(dòng)載荷對(duì)沖擊地壓?jiǎn)?dòng)的影響因子。

當(dāng)b＝0時(shí)，式（3）代表沖擊地壓類型為純靜載型沖擊地壓；式（3）的表達(dá)適用于不同的沖擊地壓理論：當(dāng)采用強(qiáng)度理論時(shí)，函數(shù)J（σx，εx）和函數(shù)D（σx，εx）即為σx和εx的強(qiáng)度函數(shù)；當(dāng)采用能量理論時(shí)，函數(shù)J（σx，εx）和函數(shù)D（σx，εx）即為σx和εx的能量函數(shù)?？傊?，在諸多靜載荷集中和動(dòng)載荷集中函數(shù)組成的范函數(shù)集合中，定然存在若干種不同的組合使得沖擊啟動(dòng)發(fā)生。

沖擊地壓?jiǎn)?dòng)理論指出[2]，煤巖系統(tǒng)內(nèi)集中靜載荷達(dá)到臨界狀態(tài)是沖擊地壓發(fā)生的必要條件。同時(shí)，沖擊地壓和突出的統(tǒng)一失穩(wěn)理論[17]指出沖擊地壓和突出發(fā)生的先決條件是煤巖體局部進(jìn)入峰值強(qiáng)度后的變形區(qū)，即沖擊地壓和突出都是發(fā)生在煤巖體內(nèi)的高應(yīng)力區(qū)。據(jù)此，可認(rèn)為沖擊地壓?jiǎn)?dòng)地點(diǎn)為煤巖體內(nèi)最高靜載集中應(yīng)力區(qū)。

葫蘆素煤礦埋深大、煤體堅(jiān)硬，這使得21103工作面回風(fēng)巷兩幫容易積聚大量彈性能，這是沖擊啟動(dòng)發(fā)生的靜荷載來(lái)源，且煤體具有強(qiáng)沖擊傾向性，在煤巖體系統(tǒng)內(nèi)已形成了沖擊啟動(dòng)體；臨空煤柱及采空區(qū)上方的厚硬砂巖頂板密度大且不易垮落，在工作面回采過程中，其大面積突然斷裂或上覆高位堅(jiān)硬頂板瞬間破斷是沖擊啟動(dòng)的動(dòng)荷載來(lái)源，臨空巷道動(dòng)靜載誘發(fā)沖擊啟動(dòng)示意圖如圖4。

圖4 臨空巷道動(dòng)靜載誘發(fā)沖擊啟動(dòng)示意圖Fig.4 Shock start induced by dynamic and static load of free entry

最高靜載集中區(qū)的靜應(yīng)力和動(dòng)載荷應(yīng)力波產(chǎn)生的動(dòng)態(tài)應(yīng)力疊加形成的應(yīng)變能超過煤巖體動(dòng)力破壞所消耗的能量時(shí)，發(fā)生沖擊地壓?jiǎn)?dòng)。即滿足：

式中：EJΩ為最高靜載集中區(qū)的靜應(yīng)力；EDΩ為動(dòng)載荷應(yīng)力波產(chǎn)生的動(dòng)態(tài)應(yīng)力；Efmin為煤巖體動(dòng)力破壞所消耗的能量。

煤巖體發(fā)生沖擊地壓等動(dòng)力破壞時(shí)釋放的能量往往遠(yuǎn)大于誘發(fā)能量，根據(jù)巖體破壞的最小能量原理[18]，無(wú)論煤巖體是在三向受力雙向受力還是單向受力狀態(tài)下破壞，其破壞真正消耗的能量總是一定的，即巖體破壞最小能量EC：

式中：σC為單軸抗拉強(qiáng)度；τC為抗剪切強(qiáng)度；E為彈性模量；G為剪切模量。

根據(jù)能量守恒，在煤巖體破壞消耗能量一定的條件下，煤巖體自身積聚的彈性能越大，沖擊地壓時(shí)伴隨的動(dòng)能就越大，該動(dòng)能來(lái)源于煤巖體破壞彈性余能WC：

式中：Ei為i維受力狀態(tài)下煤巖體中儲(chǔ)存的彈性能，由源于靜載荷集中彈性能EJΩ和動(dòng)應(yīng)力彈性能EDΩ構(gòu)成。

根據(jù)縱波衰減Hudson模型中應(yīng)力波衰減系數(shù)和動(dòng)態(tài)應(yīng)力關(guān)系式[14]可推導(dǎo)出動(dòng)態(tài)應(yīng)力表達(dá)式：

式中：σ（t）為動(dòng)態(tài)應(yīng)力；K、M為煤巖體本身物理力學(xué)性質(zhì)有關(guān)的常數(shù)；α為應(yīng)力波衰減系數(shù)；φ為煤巖體動(dòng)力破壞所消耗的能量；c為材料常數(shù)。

動(dòng)載荷應(yīng)力波傳遞至最高靜載集中應(yīng)力區(qū)時(shí)，單位體積內(nèi)動(dòng)態(tài)應(yīng)變能Vε為：

最高靜載應(yīng)力集中區(qū)單位體積內(nèi)積聚的彈性應(yīng)變能[19]VJε為：

式中：σ1、σ2、σ3均為主應(yīng)力；μ為泊松比。

由式（6）、式（8）、式（9）沖擊地壓發(fā)生時(shí)，沖擊啟動(dòng)區(qū)煤巖體破壞后彈性余能為：

該部分彈性余能以波的形式向周圍傳播[20-21]，除去耗散在21103工作面回風(fēng)巷圍巖介質(zhì)中的部分能量外，剩下的轉(zhuǎn)化成動(dòng)能釋放到巷道空間，造成煤柱側(cè)瞬間幫鼓、錨桿錨索崩落等沖擊顯現(xiàn)。

3.3 頂板壓裂孔和大直徑鉆孔卸壓技術(shù)

21103工作面回風(fēng)巷"鉆、切、壓"施工布置圖如圖5。

圖5 21103工作面回風(fēng)巷“鉆、切、壓”施工布置圖Fig.5“Drilling,cutting and pressing”construction layout of return air roadway in 21103 working face

針對(duì)堅(jiān)硬煤體，在21103工作面回風(fēng)巷兩幫施工了大直徑鉆孔，鉆孔直徑150 mm，孔間距1 m。

針對(duì)堅(jiān)硬頂板，采用高壓水泵配合專用高壓切縫鉆桿、高壓封隔器，進(jìn)行超前預(yù)裂，預(yù)裂孔每組2孔。煤柱幫側(cè)鉆孔方位角270°、仰角60°、孔徑56 mm、孔間距15 m、孔深24 m；回采幫側(cè)鉆孔方位角135°、仰角45°、孔徑56 mm、孔間距15 m、孔深30 m。

3.4 大直徑鉆孔和頂板壓裂孔防沖效果監(jiān)測(cè)

3.4.1 大直徑鉆孔防沖效果監(jiān)測(cè)

21103工作面回風(fēng)巷幫部大直徑卸壓前后微震及沖擊事件的時(shí)間序列如圖6。

圖6 21103工作面回風(fēng)巷幫部大直徑卸壓前后微震及沖擊事件的時(shí)間序列Fig.6 Time series of microseismic and shock events before and after large diameter pressure relief in return air side of 21103 working face

回風(fēng)巷幫部大直徑卸壓鉆孔實(shí)施前，21103工作面回風(fēng)巷104J以上能量事件共計(jì)發(fā)生過4次，其中2月6日和2月15日的2起4次方能量事件誘發(fā)了沖擊地壓，誘發(fā)沖擊比例達(dá)到50%，巷道最大破壞長(zhǎng)度達(dá)48 m。實(shí)施煤層大直徑卸壓鉆孔后，共發(fā)生104J以上能量事件9次，僅有1起4次方能量事件誘發(fā)了輕微沖擊，誘沖比例為11%。上述分析表明，煤層大直徑卸壓鉆孔的施工雖然沒有遏制高能量事件頻發(fā)態(tài)勢(shì)，但可以明顯的降低高能量事件的誘發(fā)沖擊地壓比例，其原因是幫部煤層大直徑卸壓鉆孔不僅使得集中應(yīng)力向圍巖深部轉(zhuǎn)移，而且增加了鉆孔附近裂隙發(fā)育情況，在巷道兩幫形成一定范圍的卸壓保護(hù)帶，一定程度上衰減了頂板斷裂引發(fā)的動(dòng)態(tài)應(yīng)力波，耗散了應(yīng)力波攜帶的彈性應(yīng)變能量，從而避免了沖擊地壓的啟動(dòng)和顯現(xiàn)。

3.4.2 頂板壓裂孔防沖效果監(jiān)測(cè)

利用頂板壓裂孔不同施工階段的微震數(shù)據(jù)來(lái)宏觀分析工作面總體預(yù)裂效果，第1階段21103綜采工作面從距離切眼753～1 322 m段回風(fēng)巷頂板采取頂板水力壓裂斷頂技術(shù)；第2階段從距離切眼1 322～1 912 m段回風(fēng)巷頂板均未采取頂板水力壓裂斷頂技術(shù)。

在第1階段21103工作面共監(jiān)測(cè)到有效微震事件7 740次，1次方以下、1～2次方、2～3次方、3～4次方、4次方及5次方能級(jí)事件所占比例分別為74.32%、19.53%、4.63%、1.37%及0.16%；在第2階段21103工作面共監(jiān)測(cè)到有效微震事件6 631次，1次方以下、1次方、2次方、3次方、4次方及5次方能級(jí)事件所占比例分別為59.60%、28.50%、9.24%、2.37%及0.27%；2個(gè)階段微震事件統(tǒng)計(jì)表見表2。

表2 2個(gè)階段微震事件統(tǒng)計(jì)表Table 2 Statistical table of microseismic events in two stages

為避免傳感器挪移或故障等情況造成小能量微震事件數(shù)量波動(dòng)，僅對(duì)2次方（即100 J）以上微震事件進(jìn)行統(tǒng)計(jì)對(duì)比分析，可見第2階段的微震事件在各能量級(jí)區(qū)間內(nèi)均顯著高于第1階段。結(jié)果表明：工作面在未進(jìn)行水力壓裂區(qū)域推采時(shí)，易引起能量的集中釋放，引發(fā)大能量事件發(fā)生；工作面在回風(fēng)巷均進(jìn)行了水力壓裂的區(qū)域推采時(shí)，工作面的安全性，明顯高于未進(jìn)行水力壓裂施工的區(qū)域。

4 結(jié)語(yǔ)

1）構(gòu)建了巷道圍巖力學(xué)結(jié)構(gòu)模型，圍巖深部彈性區(qū)用完全彈性體元件模擬，圍巖極限平衡區(qū)可用彈黏體元件串聯(lián)模擬。

2）大直徑鉆孔改變了巷道圍巖的力學(xué)結(jié)構(gòu)模型，支承壓力峰值以外人為產(chǎn)生卸壓區(qū)使得完全彈性體元件結(jié)構(gòu)模型變?yōu)閺楌んw串聯(lián)元件模型，黏體元件可以吸收部分應(yīng)變，而應(yīng)力不增加，相當(dāng)于在巷道周圍形成卸壓保護(hù)帶，該保護(hù)帶也可以衰減動(dòng)荷載應(yīng)力波的傳播。

3）動(dòng)荷載應(yīng)力波在預(yù)裂后頂板內(nèi)傳播時(shí)會(huì)以摩擦效應(yīng)、散射效應(yīng)、流動(dòng)效應(yīng)、裂紋擴(kuò)展效應(yīng)等機(jī)制發(fā)生衰減；頂板壓裂孔可以原位改變頂板巖層的物理力學(xué)性質(zhì)，降低頂板巖層的完整性和強(qiáng)度，減小沖擊地壓危險(xiǎn)。

4）對(duì)比葫蘆素煤礦21103工作面回風(fēng)巷煤層大直徑鉆孔和頂板壓裂孔施工前后的沖擊地壓顯現(xiàn)情況和微震能量事件，發(fā)現(xiàn)21103工作面回風(fēng)巷卸壓工程方案取得了好的效果，有效地減少了沖擊地壓顯現(xiàn)的發(fā)生。