張照偉

（中煤天津地下工程智能研究院有限公司，天津 300121）

隨著煤炭和石油等資源開采深度的不斷拓展，常規(guī)能源的開發(fā)程度隨采深逐漸受限，煤層氣和頁(yè)巖油氣等非常規(guī)油氣資源開發(fā)利用具備廣闊前景，成為當(dāng)前國(guó)內(nèi)外學(xué)者研究的重點(diǎn)[1]。煤層氣等資源的開采主要針對(duì)煤巖儲(chǔ)層裂紋網(wǎng)絡(luò)改造，主要通過(guò)水力壓裂工藝形成較為復(fù)雜的裂紋網(wǎng)絡(luò)體系，從而提高煤層氣資源的抽采效率[2-4]。

實(shí)際煤巖儲(chǔ)層地質(zhì)結(jié)構(gòu)復(fù)雜，特別是復(fù)合體層狀儲(chǔ)層具有縱向巖性復(fù)雜、層間物性差異較大，具有非均質(zhì)性強(qiáng)且天然裂隙發(fā)育等特點(diǎn)，需考慮水力裂縫在不同巖性介質(zhì)中擴(kuò)展特征的差異，采用穿層壓裂技術(shù)綜合開采煤層氣資源具有重要意義[5]。由于復(fù)合體層間巖石力學(xué)性質(zhì)、地應(yīng)力、結(jié)構(gòu)弱面及壓裂工藝參數(shù)的不同，水壓裂縫呈現(xiàn)的形態(tài)差異較大[6-7]；針對(duì)水壓裂縫在多巖性復(fù)合體儲(chǔ)層中的延伸擴(kuò)展方面，國(guó)內(nèi)外學(xué)者進(jìn)行了大量的研究。TAN 等[8-9]對(duì)比研究了煤巖、砂巖及頁(yè)巖不同組合條件下，水壓裂縫從煤巖層頂?shù)装迤鹆鸭翱缃缑娲訑U(kuò)展規(guī)律，分析了結(jié)構(gòu)面膠結(jié)強(qiáng)度對(duì)水壓裂縫穿層擴(kuò)展的影響作用，掌握了水壓裂縫垂向非對(duì)稱延伸特征；LIU 等[10]在考慮分層介質(zhì)中方位角、不同井斜角及射孔參數(shù)的情況下，對(duì)水壓裂的起裂擴(kuò)展規(guī)律進(jìn)行了研究；馬衍坤等[11-12]研究了跨界面水力壓裂增滲作用機(jī)制，指出水力壓裂通過(guò)引起巖體破碎膨脹而使得煤層卸壓增透，提高煤層透氣性；XING 等[13]研究了裂縫垂直擴(kuò)展行為，在考慮層間應(yīng)力差、界面膠結(jié)強(qiáng)度、縫內(nèi)凈壓力及垂向應(yīng)力差的影響下，構(gòu)建了多參數(shù)影響的控制模型；武鵬飛等[14]開展了大尺寸水力壓裂試驗(yàn)，發(fā)現(xiàn)水壓裂縫跨界面擴(kuò)展時(shí)，容易在界面處形成偏轉(zhuǎn)型、貫穿型、止裂型裂縫形態(tài)；付世豪等[15]開展了真三軸水力壓裂物理模擬試驗(yàn)，研究了多因素對(duì)水力裂縫垂向擴(kuò)展行為的影響，發(fā)現(xiàn)了水力裂縫在縱向上呈現(xiàn)非對(duì)稱、非平面擴(kuò)展特征，依據(jù)水力裂縫與巖性界面不同的作用方式，水力裂縫具有停止、轉(zhuǎn)向、分叉、穿透等多種復(fù)雜擴(kuò)展模式；GUO 等[16]通過(guò)ABAQUS 數(shù)值模擬軟件對(duì)四川盆地砂巖-頁(yè)巖相互組合的儲(chǔ)層中裂縫與界面之間的相互作用行為進(jìn)行了計(jì)算，發(fā)現(xiàn)提高壓裂液注入壓力和減小壓裂液黏度更有利于使得水力裂縫克服不同巖性儲(chǔ)層的阻力穿過(guò)界面形成穿層裂縫；ZOU 等[17]通過(guò)實(shí)驗(yàn)室物理試驗(yàn)和數(shù)值模擬結(jié)合的方式，對(duì)滲透性層理和水力裂縫之間的相互作用進(jìn)行了研究，發(fā)現(xiàn)處于較低壓裂液黏度時(shí)在天然層理面的濾失作用會(huì)增強(qiáng)，大排量注入會(huì)有利于水力裂縫穿過(guò)層里面擴(kuò)展；LI 等[18]通過(guò)滲流-損傷-應(yīng)力耦合程序?qū)Ρ』映练e巖組合儲(chǔ)層進(jìn)行了水力壓裂模擬，發(fā)現(xiàn)水力裂縫在巖石交界面出的轉(zhuǎn)折程度增加會(huì)減小滲流通道的開度。

綜上，目前煤巖復(fù)合體水力壓裂的研究多集中在煤巖組分和壓裂參數(shù)等方面，水力裂縫在不同組合巖體條件下擴(kuò)展的完整過(guò)程尚不明確，考慮煤巖復(fù)合體水力裂縫擴(kuò)展過(guò)程中與煤巖交界面以及層理等天然裂縫的復(fù)雜相互作用行為還需要深入研究。為此，基于黏聚力單元法模擬煤巖內(nèi)部基質(zhì)縫網(wǎng)結(jié)構(gòu)，結(jié)合ABAQUS 數(shù)值模擬軟件開展煤巖復(fù)合體水力壓裂數(shù)值計(jì)算，針對(duì)應(yīng)力場(chǎng)、煤巖交界面強(qiáng)度和煤層層理傾角等因素對(duì)煤巖復(fù)合體水力壓裂裂紋擴(kuò)展規(guī)律的影響進(jìn)行研究，以期為煤巖復(fù)合體儲(chǔ)層裂紋網(wǎng)絡(luò)改造提供參考理論基礎(chǔ)。

1 黏聚力單元法

1.1 黏聚力單元簡(jiǎn)介

黏聚力指的是在同種物質(zhì)內(nèi)部相鄰各部分之間的相互吸引力，這種相互吸引力是同種物質(zhì)分子之間存在分子力的微觀表現(xiàn)，只有在物質(zhì)內(nèi)分子間的距離十分接近時(shí)，一般為小于10-6cm 時(shí)，這種黏聚力才顯現(xiàn)出來(lái)[19]。在水力壓裂裂縫擴(kuò)展模擬過(guò)程中，黏聚力單元本身并不代表任何材料，只是黏聚力單元中存在的黏聚力可以抵抗裂縫尖分離時(shí)產(chǎn)生的拉應(yīng)力，正是因?yàn)轲ぞ哿卧哂写颂卣鳎攀沟迷谒毫蚜芽p擴(kuò)展模擬過(guò)程中，基于黏聚力單元的各類模型得到了廣泛的應(yīng)用。

黏聚力單元嵌入在煤巖體實(shí)體單元內(nèi)，當(dāng)煤巖體發(fā)生破壞發(fā)生時(shí)，黏聚力單元張開從而引起裂縫的起裂及擴(kuò)展。因此模擬水力裂縫的擴(kuò)展主要是依賴于黏聚力單元的嵌入，而水力裂紋如何擴(kuò)展，這取決于黏聚力單元的網(wǎng)格劃分，其網(wǎng)格精度對(duì)數(shù)值模擬結(jié)果有著重要的影響[20]。

1.2 力學(xué)損傷破壞準(zhǔn)則

在黏聚力單元破壞理論中，裂縫擴(kuò)展過(guò)程是裂縫尖端克服黏聚力作用并發(fā)生分離斷裂，裂縫起裂與擴(kuò)展由牽引-分離準(zhǔn)則控制。

1）初始損傷。裂縫牽引力的拉伸組分主要由裂縫流體壓力引起，其剪切組分由地應(yīng)力差和煤巖內(nèi)部天然裂縫的存在引起的局部剪應(yīng)力場(chǎng)所誘導(dǎo)。對(duì)于拉-剪混合破壞模式的裂縫，其初始損傷位移并不為定值，而是與具體荷載條件下的張拉-剪切混合方式及混合比例有關(guān)。采用拉-剪混合模式的二階應(yīng)力準(zhǔn)則預(yù)測(cè)初始損傷，如式（1）[21]：

式中：σt、σs為拉伸強(qiáng)度和剪切強(qiáng)度，Pa；tn、ts為法向牽引力和切向牽引力，Pa。

2）損傷演化。通過(guò)引入損傷因子D來(lái)表征黏聚力單元的損傷程度，初始損傷發(fā)生后，D從0單調(diào)遞增至1。則由損傷引起的應(yīng)力變化可表示為：

式中：tn為描述拉-剪混合模式下?lián)p傷演化；、為當(dāng)前分離位移在線彈性法則下的應(yīng)力組分。

1.3 流體流動(dòng)準(zhǔn)則

假設(shè)縫內(nèi)壓裂液選取為不可壓縮的牛頓流體，符合牛頓黏性定律：

縫內(nèi)流體切向流動(dòng)方程為：

式中：qt為單位面積壓裂液流速，m/s；w為裂縫寬度，m；?pw裂縫方向的壓力梯度，Pa/m；可以理解為滲透性或流動(dòng)阻力。

縫內(nèi)向煤巖多孔基質(zhì)濾失行為的法向流定義如下：

式中：qa、qb為單元頂面和底面的法向?yàn)V失流動(dòng)速率，m/s；ca、cb為頂?shù)酌鏋V失系數(shù)，m/（Pa·s）；pm為 單元中面裂隙流壓力，Pa；pa、pb為單元頂?shù)酌嫣幍目紫秹毫Γ琍a。

水力壓裂黏聚力模型如圖1，模型中的裂縫分為黏性斷裂區(qū)與黏性非斷裂區(qū)。

圖1 水力壓裂黏聚力模型Fig.1 Hydraulic fracturing cohesion model diagram

2 煤巖復(fù)合體水力壓裂數(shù)值模型及試驗(yàn)方案

1）煤巖復(fù)合體水力壓裂數(shù)值模型。煤巖復(fù)合體水力壓裂數(shù)值計(jì)算模型如圖2，整個(gè)模型尺寸與壓裂區(qū)尺寸之比為5∶1，模型共劃分46 857 個(gè)單元，其中孔隙流體-應(yīng)力耦合平面應(yīng)變單元（類型CPE4P）為34 533 個(gè)，黏聚力孔壓?jiǎn)卧愋虲OH2D4P）為12 324 個(gè)。模型邊界節(jié)點(diǎn)自由度采用固支約束，即約束x方向和y方向自由度，xy平面內(nèi)的轉(zhuǎn)動(dòng)自由度不被約束。模型物理力學(xué)參數(shù)見(jiàn)表1。

表1 數(shù)值計(jì)算模型主要物理力學(xué)參數(shù)表Table 1 Table of physical and mechanical parameters of numerical calculation model

圖2 煤巖復(fù)合體水力壓裂數(shù)值模擬模型Fig.2 Numerical simulation model of hydraulic fracturing of coal rock complex

2）試驗(yàn)方案。將應(yīng)力場(chǎng)、煤巖交界面強(qiáng)度和煤層層理作為單因素控制變量來(lái)研究不同條件下水壓裂紋擴(kuò)展規(guī)律，x水平方向和y垂直方向應(yīng)力場(chǎng)（σx，σy）設(shè)置4 組，分別為（5 MPa，2.5 MPa）、（5 MPa，5 MPa）、（5 MPa，7.5 MPa）、（5 MPa，10 MPa）、（5 MPa，12.5 MPa）、（5 MPa，15 MPa）；煤巖交界面彈性模量分別為1.5、2.5、5 GPa；煤體層理傾角（層理與x水平方向）分別設(shè)置為0°、30°、60°；壓裂液排量設(shè)置為0.006 m3/s。

3 數(shù)值模擬結(jié)果

3.1 應(yīng)力場(chǎng)對(duì)煤巖復(fù)合體水力壓裂的影響

為了分析不同應(yīng)力條件對(duì)煤巖復(fù)合體水力壓裂裂紋擴(kuò)展規(guī)律的影響，設(shè)置x方向應(yīng)力σx為5 MPa 不變，只改變y方向應(yīng)力σy，間隔為2.5 MPa，進(jìn)行對(duì)比水力壓裂模擬試驗(yàn)研究。裂紋擴(kuò)展過(guò)程中孔隙壓力數(shù)值模擬結(jié)果如圖3，水力壓裂裂紋張開度結(jié)果對(duì)比云圖如4，不同應(yīng)力條件下注入壓力全程變化對(duì)比曲線如圖5。

圖3 不同應(yīng)力條件下孔隙壓力數(shù)值模擬結(jié)果對(duì)比云圖Fig.3 Comparison of numerical simulation results of pore pressure under different stress conditions

由圖3 和圖4 可知：應(yīng)力對(duì)水力裂紋的擴(kuò)展過(guò)程和形態(tài)有重要影響作用。在應(yīng)力場(chǎng)（σx，σy）=（5 MPa，2.5 MPa）的條件下，水壓裂縫在上覆巖層中起裂并沿著最大主應(yīng)力方向擴(kuò)展，水力裂縫在上煤巖交界面水平平行擴(kuò)展并未沿著煤巖交界弱面擴(kuò)展，由于上覆巖力學(xué)特性較強(qiáng)且層脆性較大，所以水壓裂縫張開度整體較小，裂縫傳播速度較快；在應(yīng)力場(chǎng)（σx，σy）=（5 MPa，5 MPa）的條件下，水壓裂縫在上覆巖層中起裂，到達(dá)上煤巖交界面后形成了2 條主裂縫，其中1 條主裂縫直接穿過(guò)煤巖交界面到達(dá)煤層中并持續(xù)在煤體里邊擴(kuò)展，另1 條主裂縫在達(dá)到上煤巖交界面后沿著上煤巖交界面上擴(kuò)展延伸，在煤巖交界面處裂縫張開度最大；在應(yīng)力場(chǎng)（σx，σy）=（5 MPa，7.5/10/12.5/15 MPa）的條件下，水壓裂縫在上覆巖層中起裂后，隨著應(yīng)力差的逐漸增大，水力壓裂裂縫受應(yīng)力場(chǎng)的控制作用越來(lái)越顯著，沿著最大主應(yīng)力直接穿過(guò)煤巖交界面到達(dá)煤層后持續(xù)沿著此方向擴(kuò)展。但隨著應(yīng)力差越來(lái)越大，水壓裂縫在上覆巖層中擴(kuò)展的距離越來(lái)越短，穿過(guò)交界面后在煤層中擴(kuò)展越來(lái)越遠(yuǎn)。

圖4 不同應(yīng)力條件下裂紋張開度數(shù)值模擬結(jié)果對(duì)比云圖Fig.4 Comparison of numerical simulation results of crack opening under different stress conditions

由圖5 可知：σx應(yīng)力為5 MPa，σy應(yīng)力依次為2.5、5、7.5、10、12.5、15 MPa 條件下的首次壓裂壓力分別為21.14、21.41、23.77、23.83、22.65、22.51 MPa；當(dāng)應(yīng)力差較低時(shí)，應(yīng)力場(chǎng)難以控制水力壓裂裂縫的擴(kuò)展方向，只能在巖體內(nèi)部產(chǎn)生擴(kuò)展裂縫；隨著應(yīng)力差增大，水力壓裂裂縫會(huì)沿著最大主應(yīng)力方向，對(duì)煤巖復(fù)合體產(chǎn)生不同程度的孔隙壓力，導(dǎo)致煤巖復(fù)合體內(nèi)部產(chǎn)生水壓裂縫；后期σy`應(yīng)力為7.5 MPa 和10 MPa 注入壓力平穩(wěn)時(shí)相對(duì)較高，σy應(yīng)力為12.5 MPa 和15 MPa注入壓力相對(duì)較小一些，是因?yàn)樾纬闪芽p的大小不同差異導(dǎo)致卸壓不同。

圖5 不同應(yīng)力條件下注入壓力全程變化對(duì)比曲線Fig.5 Comparison curves of injection pressure change under different stress conditions

3.2 交界面強(qiáng)度對(duì)煤巖復(fù)合體水力壓裂的影響

為分析交界面強(qiáng)度對(duì)煤巖復(fù)合體水力壓裂裂紋擴(kuò)展規(guī)律的影響，設(shè)置應(yīng)力場(chǎng)在（σx，σy）=（5 MPa，10 MPa）條件下，改變煤巖交界面的彈性模量進(jìn)行對(duì)比水力壓裂模擬試驗(yàn)研究。不同交界面彈性模量條件下水力壓裂模擬結(jié)果對(duì)比云圖如圖6。

圖6 不同交界面彈層模量條件下水力壓裂模擬結(jié)果對(duì)比云圖Fig.6 Comparison cloud charts of hydraulic fracturing simulation results under different interface elastic layer modulus conditions

如圖6（a），煤巖交界面強(qiáng)度1.5 GPa，水力壓裂裂縫從巖層擴(kuò)展并貫穿交界面，但是由于交界面強(qiáng)度較弱難以形成水壓聚集達(dá)到煤層的開裂條件，裂縫最終沿著煤層上表面向兩側(cè)擴(kuò)展，此條件下水壓裂縫不能誘導(dǎo)煤層出現(xiàn)裂縫。

如圖6（b），煤巖交界面強(qiáng)度2 GPa 和2.5 GPa水壓裂縫能夠穿越交界面并最終在煤層中擴(kuò)展延伸，但是此時(shí)交界面強(qiáng)度較弱水壓裂縫不僅穿越煤巖交界面并誘導(dǎo)煤層起裂擴(kuò)展并且在煤層上表面也沿交界面擴(kuò)展一定距離。

如圖6（c），當(dāng)交界面強(qiáng)度超過(guò)煤層時(shí)由于煤巖交界面強(qiáng)度較大破裂所需壓力大于煤層，所以基本上不存在沿煤層上表面擴(kuò)展的現(xiàn)象，裂縫直接穿越交界面并在煤層中擴(kuò)展。

從水壓裂縫在煤層中擴(kuò)展的長(zhǎng)度和縫寬來(lái)講，相同應(yīng)力條件下交界面強(qiáng)度越大在煤層中擴(kuò)展的長(zhǎng)度越長(zhǎng)，同時(shí)界面強(qiáng)度越大，水壓裂縫在煤層中的縫寬越寬。特別的，當(dāng)交界面強(qiáng)度很大超過(guò)煤層時(shí)由于煤巖交界面強(qiáng)度大破裂所需破裂壓力大于煤層，所以基本上不存在沿煤層上表面擴(kuò)展的現(xiàn)象，裂縫直接穿越交界面并在煤層中擴(kuò)展。同時(shí)從水壓裂縫在煤層中擴(kuò)展的長(zhǎng)度和縫寬來(lái)講，相同應(yīng)力條件下交界面強(qiáng)度越大在煤層中擴(kuò)展的長(zhǎng)度越長(zhǎng)，同時(shí)界面強(qiáng)度越大，水壓裂縫在煤層中的縫寬越寬。

結(jié)合應(yīng)力分析：煤巖復(fù)合體水力壓裂裂縫穿過(guò)交界面并在煤層中擴(kuò)展存在閾值，主應(yīng)力差，交界面強(qiáng)度對(duì)此閾值存在重要的影響。主應(yīng)力差相同時(shí)，煤巖交界面強(qiáng)度越大越容易形成水壓聚集達(dá)到煤層破裂擴(kuò)展的閾值，同時(shí)界面強(qiáng)度越大水力壓裂裂縫穿越煤巖交界面后更容易在煤層中擴(kuò)展。

3.3 交界面強(qiáng)度對(duì)煤巖復(fù)合體水力壓裂的影響

為分析煤層層理傾角對(duì)煤巖復(fù)合體水力壓裂裂縫擴(kuò)展規(guī)律的影響，在（σx，σy）=（5 MPa，10 MPa）應(yīng)力場(chǎng)條件下，改變煤層層理傾角進(jìn)行對(duì)比水力壓裂模擬試驗(yàn)研究。不同煤層層理傾角條件下水力壓裂模擬結(jié)果對(duì)比云圖如圖7。

圖7 不同煤層層理傾角條件下水力壓裂模擬結(jié)果對(duì)比云圖Fig.7 Comparison cloud charts of hydraulic fracturing simulation results under different bedding angles of coal seams

由圖7 可知：層理傾角為0°和30°條件下，主裂縫均穿過(guò)層理面并在層理處形成張開裂縫，其中層理傾角30°條件下在層理面擴(kuò)展程度較為明顯；層理傾角為60°條件下未能穿過(guò)層理面，但在層理面出現(xiàn)明顯的裂縫張開滑移；層理傾角為0°時(shí)水力裂縫可以穿過(guò)煤巖交界面并在煤層中穿過(guò)層理擴(kuò)展，在主裂縫與層理遭遇時(shí)形成一定數(shù)量的分值裂縫；層理傾角為30°時(shí)水力裂縫在穿過(guò)交界面時(shí)先在交界面中擴(kuò)展一定距離后沿多個(gè)層里面延伸，最終貫穿層理面；層理傾角為60°時(shí)，主裂縫穿越交界面后在煤層中受層理約束沿層理擴(kuò)展未能貫穿層理面壓裂的作用范圍僅在兩相鄰層理之內(nèi)的煤層中，說(shuō)明此條件下層理弱面對(duì)裂縫擴(kuò)展方向的控制作用強(qiáng)于地應(yīng)力。

4 結(jié)語(yǔ)

1）煤巖復(fù)合體水壓裂縫在上覆巖層中起裂后，隨著σy與σx應(yīng)力差的逐漸增大，水力壓裂裂縫受應(yīng)力場(chǎng)的控制作用越來(lái)越顯著，沿最大主應(yīng)力方向持續(xù)擴(kuò)展。在應(yīng)力場(chǎng)（σx，σy）=（5 MPa，5 MPa）的條件下，水壓裂縫在上覆巖層中起裂后在煤巖交界面和煤層中擴(kuò)展形成了的2 條主水力裂縫。

2）煤巖交界面強(qiáng)度較小時(shí)，在交界面中難以形成水壓積聚達(dá)到煤層的開裂條件，裂縫最終沿著煤層上表面向兩側(cè)擴(kuò)展，不能誘導(dǎo)煤層出現(xiàn)裂縫。煤巖交界面強(qiáng)度較大時(shí)水壓裂縫能夠穿越交界面并最終在煤層中擴(kuò)展延伸，能夠在煤巖交界面中擴(kuò)展一定距離并完成憋壓，當(dāng)壓力達(dá)到煤層開裂條件時(shí)水壓裂縫穿越煤巖交界面并誘導(dǎo)煤層起裂擴(kuò)展；煤巖交界面強(qiáng)度很大超過(guò)煤層強(qiáng)度時(shí)，水壓裂縫基本上不沿煤層上表面擴(kuò)展，而是裂縫直接穿越交界面并在煤層中擴(kuò)展。同時(shí)從水壓裂縫在煤層中擴(kuò)展的長(zhǎng)度和縫寬來(lái)講，相同應(yīng)力條件下交界面強(qiáng)度越大在煤層中擴(kuò)展的長(zhǎng)度越長(zhǎng)，同時(shí)界面強(qiáng)度越大，水壓裂縫在煤層中的縫寬越寬。

3）層理傾角為0°時(shí)水力裂縫可以穿過(guò)煤巖交界面并在煤層中穿過(guò)層理擴(kuò)展，在主裂縫與層理遭遇時(shí)形成一定數(shù)量的分值裂縫；層理傾角為30°時(shí)水力裂縫在穿過(guò)交界面時(shí)先在交界面中擴(kuò)展一定距離后沿多個(gè)層里面延伸，最終貫穿層理面；層理傾角為60°時(shí)，主裂縫穿越交界面后在煤層中受層理約束沿層理擴(kuò)展未能貫穿層理面壓裂的作用范圍僅在兩相鄰層理之內(nèi)的煤層中，說(shuō)明此條件下層理弱面對(duì)裂縫擴(kuò)展方向的控制作用強(qiáng)于地應(yīng)力。