王 博，姜在炳，杜天林，賈立龍，牟全斌，楊建超，陳崇楓，李浩哲

(1. 中煤科工集團(tuán)西安研究院有限公司，陜西西安 710054；2. 西安科技大學(xué)，陜西西安 710054)

井下瓦斯(煤層氣)抽采被認(rèn)為是煤礦安全的一項(xiàng)重要技術(shù)措施，但常規(guī)的井下鉆孔瓦斯抽采半徑小，瓦斯?jié)舛刃。椴尚Ч?，不能滿足發(fā)電的要求，大部分井下瓦斯均被直接排放到大氣中，根據(jù)相關(guān)研究，1m3甲烷的溫室效應(yīng)是二氧化碳的20余倍[1]，甲烷的排放造成了對大氣和生態(tài)環(huán)境的嚴(yán)重污染[2-3]。利用地面煤層氣開發(fā)技術(shù)對煤儲層及其頂?shù)装暹M(jìn)行壓裂改造可以提高瓦斯?jié)舛?，達(dá)到利用瓦斯利用的預(yù)期效果[4-6]。由于高瓦斯地區(qū)煤層多為碎軟煤層，其結(jié)構(gòu)破碎，在煤層內(nèi)定向鉆孔易塌孔，施工難度大，因此選擇在緊鄰煤層頂板中施工大直徑定向長鉆孔，然后在頂板內(nèi)向下定向射孔，實(shí)施分段壓裂，能夠?qū)崿F(xiàn)瓦斯高效抽采。煤層頂板水平井分段壓裂模式抽采的瓦斯?jié)舛雀?，基本可?00%利用，有利于大幅降低甲烷排放量，增加清潔能源供給，保護(hù)生態(tài)環(huán)境，未來的社會(huì)效益和環(huán)境效益巨大，是實(shí)現(xiàn)“雙碳”目標(biāo)的可靠路徑。煤層頂板分段水力壓裂技術(shù)目前已在淮北廬嶺礦區(qū)、山西保德及陜西韓城的礦區(qū)進(jìn)行了工程試驗(yàn)，取得了較好的效果[7-8]。以往研究更多針對壓裂縫在垂向上能否穿過頂板進(jìn)入煤層開展，并未對裂縫穿層后的延伸規(guī)律深入研究；分段優(yōu)化設(shè)計(jì)在頁巖氣開發(fā)等方面應(yīng)用較多，在瓦斯抽采方面分段參數(shù)如何設(shè)置的問題還有待探討。文章以淮北礦區(qū)祁東煤礦為例，通過煤層頂板水平長鉆孔裂縫穿層延展特性以及壓裂縫長度和壓裂縫間距的優(yōu)化的研究，闡述了煤層頂板水平井分段壓裂瓦斯治理模式，可為區(qū)域化瓦斯治理的工程布置方式提供理論依據(jù)。

1 工作區(qū)概況

祁東煤礦位于宿南向斜的東南端，其構(gòu)造形態(tài)基本為一走向近東西、向北緩傾的單斜構(gòu)造。含煤地層為二疊系山西組(P1s)以及下石河子組(P1x)，主要可采煤層為61、71、82和9煤層。

本次以71煤層為研究對象，71煤層為本礦主采煤層，煤厚2.61～3.91m，平均3.50m，煤層厚度變化較大，煤層結(jié)構(gòu)較復(fù)雜，多含1～2層泥巖夾矸，夾矸厚0.28～1.04m，平均0.51m。煤層埋深680～710 m。直接頂為深灰色泥巖、粉砂巖；老頂為細(xì)粒砂巖，黑灰色，薄層狀，水平層理。直接底為深灰色泥巖，塊狀，含較多植物化石碎片及植物根莖化石。通過交叉偶極子聲波測井方法，解釋得到71煤層垂直應(yīng)力平均為17.0 MPa，最大水平主應(yīng)力為13.0 MPa，最小水平主應(yīng)力值為5.0 MPa。71煤層頂板垂直應(yīng)力為17.0 MPa，最大主應(yīng)力為13.0MPa，最小主應(yīng)力值為6.0 MPa。頂板及煤層的垂向應(yīng)力均大于水平應(yīng)力，水力壓裂過程中煤層將形成垂直裂縫;頂?shù)装宓淖钚∷街鲬?yīng)力大于煤層的水平主應(yīng)力，能夠使頂板中形成的垂直裂縫擴(kuò)展延伸到煤層，提高煤層的增產(chǎn)效果。

煤層頂板水平井分段壓裂技術(shù)工藝，即將水平井設(shè)置在緊鄰煤層的頂板巖層中，可以規(guī)避在碎軟低滲煤層中鉆水平井時(shí)易垮孔、固井質(zhì)量差的問題，并實(shí)施跨煤巖界面穿層壓裂溝通井筒與下部煤層，利用煤層頂板相對較強(qiáng)的脆性造長縫，提高壓裂改造效果，實(shí)現(xiàn)瓦斯高效抽采。煤層頂板水平井分段壓裂模如圖1所示。

圖1 煤層頂板水平井分段壓裂模式示意圖Figure 1 Schematic diagram of staged fracturing mode of horizontal well in the roof of coal seam

2 裂縫穿層特性研究

水力壓裂穿過煤巖層界面及在煤巖體內(nèi)的擴(kuò)展特性是施工效果的關(guān)鍵因素[9-12]。為了研究煤層頂板壓裂水力裂縫的穿層特性，利用數(shù)值模擬軟件建立了基于頂板水平井穿層壓裂數(shù)值模型[13-15]，相比于常規(guī)有限元方法，擴(kuò)展有限元方法模塊具有計(jì)算結(jié)果精度高和計(jì)算量小的優(yōu)點(diǎn)，可對水力壓裂過程中多種儲層物性參數(shù)及壓裂施工參數(shù)進(jìn)行模擬分析，且裂縫形態(tài)逼真，裂縫面凹凸程度清晰，結(jié)果準(zhǔn)確[16-18]。

2.1 模型參數(shù)設(shè)置

通過定義裂縫起裂準(zhǔn)則和損傷演化準(zhǔn)則，建立應(yīng)力-滲流-損傷耦合關(guān)系，模擬計(jì)算水力壓裂裂縫擴(kuò)展規(guī)律[19-21]。模型從上至下設(shè)計(jì)與實(shí)際地層接近的地層模型，如圖2所示，劃分為4層，即細(xì)粒砂巖5 m，砂質(zhì)泥巖3 m，煤層3.5 m，泥巖5 m，水平井的位置部署在距離煤層頂部1 m的砂質(zhì)泥巖頂板內(nèi)。因?yàn)閴毫涯Ｐ蜑閷ΨQ方式，因此本次模擬選取裂縫的單翼進(jìn)行研究[22-24]。模擬的模型長寬高設(shè)定為100 m×40 m×16 m，模型劃分形成12 600個(gè)單元。模型中插入采用擴(kuò)展有限元法，應(yīng)用應(yīng)力-滲流-損傷耦合理論模型來模擬水壓致裂裂縫擴(kuò)展過程，提取了不同時(shí)刻形成的裂縫形態(tài)和應(yīng)力分布狀態(tài)圖，如圖3所示，分析了頂板水平井穿層壓裂裂縫動(dòng)態(tài)擴(kuò)展過程。

圖2 水力壓裂數(shù)值模擬模型Figure 2 Numerical simulation model of hydrofracturing

a.5min裂縫形態(tài)

b.10min裂縫形態(tài)

c.15min裂縫形態(tài)

d.20min裂縫形態(tài)圖3 煤層頂板壓裂裂縫延展圖Figure 3 Extension diagram of fracturing fractures in coal seam roof

水平井壓裂施工排量為5 m3/min，壓裂液黏度0.96 Pa·s；巖石力學(xué)參數(shù)、孔隙率、滲透率、地應(yīng)力等參數(shù)采用了淮北礦區(qū)祁東煤礦實(shí)際測試數(shù)據(jù)(表1)。

表1 頂板壓裂模型計(jì)算參數(shù)

2.2 數(shù)值模擬結(jié)果分析

從圖3中可以看出，裂縫擴(kuò)展是一個(gè)漸進(jìn)過程，壓裂初始階段裂縫主要在煤層頂板泥巖內(nèi)擴(kuò)展，此時(shí)裂縫會(huì)沿著縱向高度和橫向長度兩個(gè)方向同時(shí)擴(kuò)展，擴(kuò)展速度較快。當(dāng)裂縫高度擴(kuò)展到上界面時(shí)，裂縫內(nèi)水壓難以壓開頂板細(xì)粒砂巖，裂縫開始向下傳遞，在定向射孔的誘導(dǎo)作用下，裂縫會(huì)向下穿越界面進(jìn)入煤層把煤層全部壓開，與此同時(shí)，水力裂縫繼續(xù)擴(kuò)展，在長度方向上裂縫擴(kuò)展增加較快，同時(shí)由于裂縫在頂板泥巖內(nèi)的擴(kuò)展速度要大于裂縫在煤層的擴(kuò)展速度，橫向上裂縫在頂板巖層內(nèi)快速延伸，從而牽引煤層中的裂縫在橫向上快速延伸。煤層的朔性比頂板泥巖大，可以形成比頂板泥巖更為寬泛的壓裂縫。

2.3 施工參數(shù)優(yōu)化

基于實(shí)際地層參數(shù)，采用MFracSuite三維壓裂系統(tǒng)軟件，進(jìn)一步研究壓裂施工排量對煤層頂板水平井壓裂裂縫擴(kuò)展形態(tài)的影響規(guī)律[25]。數(shù)值模擬時(shí)，固定其他參數(shù)不變，設(shè)定水平井與軟煤層的距離為1.0 m，分別模擬施工排量為3.0 m3/min、4.0 m3/min、5.0 m3/min、6.0 m3/min時(shí)裂縫形態(tài)。模擬結(jié)果如圖4所示：煤層頂板水平井壓裂都能把下伏煤層全部壓開，且隨著壓裂施工排量的增加，縫長和裂縫高度都有增加，但是縫高的增加更加明顯。對于水力壓裂而言，裂縫縫高盡量控制在煤層有效厚度內(nèi)，如果縫高超過儲層有效厚度，將會(huì)造成支撐劑的無效支撐，甚至?xí)捍└魧?，影響壓裂效果，降低增產(chǎn)效果。因此，根據(jù)祁東煤礦地質(zhì)特征，施工排量控制在5m3/min左右較合理。

a.注入排量為 3m3 / min 時(shí)的裂縫延伸形態(tài)

b.注入排量為 4m3 / min 時(shí)的裂縫延伸形態(tài)

c.注入排量為 5m3 / min 時(shí)的裂縫延伸形態(tài)

d.注入排量為 6m3 / min 時(shí)的裂縫延伸形態(tài)圖4 不同注入排量下的裂縫延伸形態(tài)Figure 4 Fracture extension morphology under different injection displacement

3 工程設(shè)計(jì)優(yōu)化研究

上述壓裂縫擴(kuò)展研究結(jié)果表明，緊鄰煤層頂板巖石的裂縫能夠延展至煤層，且在頂板巖石的撕裂作用下，煤層內(nèi)的裂縫能夠延展一定的距離。但是如何設(shè)置裂縫的形態(tài)才能得到最優(yōu)的抽采效果，是區(qū)域化瓦斯治理必須解決的問題。因此在考慮裂縫穿層機(jī)理的基礎(chǔ)上，研究裂縫長度以及裂縫間距等設(shè)計(jì)參數(shù)，對瓦斯區(qū)域治理具有重要實(shí)用價(jià)值[26-28]。

以下應(yīng)用Eclipse產(chǎn)能模擬軟件，采用數(shù)值模擬的方法，從裂縫參數(shù)對累計(jì)產(chǎn)氣量的影響方面，提出最優(yōu)裂縫設(shè)計(jì)參數(shù)。

3.1 產(chǎn)能模擬參數(shù)設(shè)置

本次模擬采用雙孔模型，同時(shí)考慮壓裂裂縫，形成煤基質(zhì)、水力壓裂裂縫的復(fù)合模型。模擬氣藏區(qū)域?yàn)?600 m×600 m，高度3.50 m，模型x方向網(wǎng)格設(shè)置60個(gè)，y方向網(wǎng)格設(shè)置60個(gè)，z方向的網(wǎng)格數(shù)為4個(gè)，在各壓裂段網(wǎng)格進(jìn)行加密。對于壓裂裂縫段，利用裂縫“等效導(dǎo)流能力法”對裂縫網(wǎng)格進(jìn)行處理。產(chǎn)能模擬的儲層參數(shù)參考擬進(jìn)行工程試驗(yàn)的淮北礦區(qū)祁東煤礦，參數(shù)設(shè)置見表2所示。

表2 瓦斯抽采產(chǎn)能模擬參數(shù)

裂縫的縫長、裂縫的間距的選取范圍是在參考現(xiàn)有施工條件及臨近礦井實(shí)際生產(chǎn)數(shù)據(jù)的基礎(chǔ)上，結(jié)合以往相關(guān)專家學(xué)者的研究成果合理選取的。根據(jù)礦區(qū)附近蘆嶺井田應(yīng)用煤層頂板水平井壓裂技術(shù)后，日最高產(chǎn)氣量為10 000 m3進(jìn)行反演，可知壓裂后裂縫導(dǎo)流能力約為100 mD，因此本次設(shè)置裂縫導(dǎo)流能力為100 mD；根據(jù)許耀波等人研究，壓裂段數(shù)過密，段間距過小會(huì)產(chǎn)生縫間應(yīng)力干擾[10]，段間距過大則瓦斯的抽采不充分，因此設(shè)定壓裂縫間距大于等于50 m，即單孔壓裂段在6～10段。在目前壓裂技術(shù)條件下，壓裂施工排量有限，據(jù)已有地面煤層氣井的微地震監(jiān)測顯示，壓裂縫的長度一般小于160 m，因此本次設(shè)置裂縫長度取值范圍在40～140 m。

3.2 裂縫縫長的選擇

裂縫長度是影響水平井產(chǎn)能的一項(xiàng)重要因素，本次設(shè)置裂縫導(dǎo)流能力為100 mD，單孔壓裂段數(shù)為6段，改變壓裂裂縫長度(40 m、60 m、80 m、100 m、120 m及140 m)，模擬不同方案4 a瓦斯含量的變化情況，如圖5所示。分析不同方案的累計(jì)產(chǎn)氣量隨抽采時(shí)間的變化規(guī)律，并根據(jù)4 a的累積量計(jì)算不同裂縫數(shù)目下累計(jì)產(chǎn)氣量增幅的同比增長率及環(huán)比增長率，分別如圖6和圖7所示。

由模擬結(jié)果可以看到，隨著壓裂裂縫長度的增加，煤層4 a累計(jì)產(chǎn)氣量是逐漸增加的，但增加幅度逐漸減慢，這是由于雖然較長的裂縫建立了井眼到遠(yuǎn)端的高滲流通道，但解吸速度較慢，無法提供足夠的氣量向通道中運(yùn)送，從而導(dǎo)致每米裂縫貢獻(xiàn)值下降，從累計(jì)產(chǎn)氣來看，當(dāng)裂縫長度大于80 m后，累計(jì)產(chǎn)氣量變化幅度趨于平緩，單條裂縫對產(chǎn)能的貢獻(xiàn)逐漸變小，結(jié)合井下施工條件及安全考慮，由于井下孔徑的限制，壓裂排量一般小于5 m3/min，平均砂比控制在10%以下，因此優(yōu)選取裂縫長度80 m。

a.縫長40m抽采4a瓦斯剩余量

b.縫長80m抽采4a瓦斯剩余量

c.縫長120m抽采4a瓦斯剩余量

d.縫長160m抽采4a瓦斯剩余量圖5 4a瓦斯剩余含量Figure 5 Residual gas content in 4a

圖6 累計(jì)產(chǎn)氣量隨抽采時(shí)間的變化Figure 6 Variation graph of cumulative gas production with drainge time

圖7 不同的裂縫長度對累計(jì)產(chǎn)氣量的影響Figure 7 Effect of different fracture length on cumulative gas production

3.3 壓裂段數(shù)的選擇

對于水平井長為500 m，縫長以80 m的等間距裂縫進(jìn)行模擬，壓裂段數(shù)目分別為6段、7段、8段、9段和10段(裂縫的間距分別為84 m、72 m、63 m、55 m和50 m)，模擬不同方案4 a瓦斯含量的變化情況，如圖8所示。分析了不同方案的累計(jì)產(chǎn)氣量隨抽采時(shí)間的變化規(guī)律，并根據(jù)4 a的累積量計(jì)算不同裂縫數(shù)目下累計(jì)產(chǎn)氣量增幅的同比增長率及環(huán)比增長率，分別如圖9和圖10所示。

a.單孔分6段抽采4a瓦斯剩余量

b.單孔分7段抽采4a瓦斯剩余量

c.單孔分8段抽采4a瓦斯剩余量

d.單孔分9段抽采4a瓦斯剩余量圖8 4a瓦斯含量變化情況Figure 8 Variation of gas content in 4a

圖9 累計(jì)產(chǎn)氣量隨抽采時(shí)間的變化Figure 9 Variation of cumulative gas production with the change of drainge time

圖10 不同的壓裂段數(shù)對累計(jì)產(chǎn)氣量的影響Figure 10 Effect of different fracturing stages on cumulative gas production

從圖9和圖10可以看出隨著裂縫數(shù)目的增加，增長率值不斷減小，即總產(chǎn)氣量的增加速度逐漸減慢。裂縫數(shù)目的增加促進(jìn)了產(chǎn)量的增長，但對一定長度水平井其裂縫的間距必就減小，相鄰裂縫由于流動(dòng)相互干擾，其流動(dòng)區(qū)域變小，抑制了產(chǎn)量的增長，所以產(chǎn)量增長呈現(xiàn)出遞減的規(guī)律。若同比增長率10 %時(shí)決定壓裂更多裂縫有較好的經(jīng)濟(jì)效益，則本井例的壓裂段劃分為8段，即裂縫間距約63 m較合適。

4 工程驗(yàn)證

祁東煤礦71號煤層埋深在680～710 m，平均煤層厚度3.50m，空氣干燥基氣含量為10.85m3/t。水平井層位布置在71號煤層頂板較堅(jiān)硬的巖層中，為了解決定向鉆進(jìn)中的技術(shù)難題，施工期間采用邊界探測儀技術(shù)，即利用隨鉆邊界探測儀實(shí)時(shí)探測井眼軌跡一定范圍的煤、巖層電阻率和自然伽馬等參數(shù)，根據(jù)71煤層與其頂板泥巖、粉砂巖的電阻率和自然伽馬差異，來判斷井眼所在位置及其與煤層頂面的距離。據(jù)此采用地質(zhì)導(dǎo)向配套工具及時(shí)調(diào)整鉆進(jìn)方向，以最大程度地確保水平段井眼軌跡在煤層頂面0.5～2m范圍內(nèi)，滿足了后續(xù)射孔、壓裂的要求，使壓裂裂縫有效溝通煤層，形成煤層氣的滲流通道。水平井水平段長度為720 m，將水平井分割成10 段進(jìn)行壓裂改造，各段間距在60～80 m，壓裂注入排量為5～10 m3/min，平均每段加砂77 m3，平均每段注入壓裂液938 m3。在壓裂施工過程中，通過微地震裂縫監(jiān)測技術(shù)對頂板水平井壓裂裂縫形態(tài)進(jìn)行監(jiān)測，結(jié)果顯示，第一段主裂縫長度80～120 m，裂縫高度10.3m，裂縫方位角357.3°，證明了頂板裂縫能夠下穿進(jìn)入煤層，且能夠延展較長的距離，實(shí)現(xiàn)對煤層的強(qiáng)化改造的目的。

水平井產(chǎn)氣最高突破5 000 m3/d，截至2021年11月，2a累計(jì)產(chǎn)氣量約87萬m3，平均單段日產(chǎn)氣量約為121 m3。

5 結(jié)論

1)采用數(shù)值模擬的方法模擬了煤巖界面處裂縫從頂板起裂后的穿層及延展規(guī)律，認(rèn)為在煤層頂板巖石的裂縫能夠延展至煤層，且在頂板巖石的撕裂作用下，煤層內(nèi)的裂縫能夠快速延展且煤層的裂縫寬度大于頂板巖石的裂縫寬度。

2)在礦井工程地質(zhì)及煤儲層參數(shù)一定的條件下，通過瓦斯抽采的效率，確定了分段壓裂的裂縫長度，以及壓裂縫間距，認(rèn)為壓裂裂縫的縫長為80 m，壓裂縫的間距為63 m時(shí)抽采效果。