天然裂縫對(duì)頁(yè)巖儲(chǔ)層水力裂縫擴(kuò)展影響數(shù)值模擬研究

沈永星，馮增朝，周 動(dòng),2 ，馬玉林

(1.太原理工大學(xué) 原位改性采礦教育部重點(diǎn)實(shí)驗(yàn)室， 山西 太原 030024；2.太原理工大學(xué) 安全與應(yīng)急管理工程學(xué)院，山西 太原 030024；3.遼寧工程技術(shù)大學(xué) 力學(xué)與工程學(xué)院，遼寧 阜新 123000)

0 引 言

水力壓裂作為一種非常規(guī)油氣藏的增產(chǎn)技術(shù)，已經(jīng)廣泛應(yīng)用在頁(yè)巖氣開(kāi)采中[1]。頁(yè)巖儲(chǔ)層為低滲透儲(chǔ)層，通過(guò)水力壓裂產(chǎn)生大量裂縫，以改造頁(yè)巖儲(chǔ)層滲透性實(shí)現(xiàn)增產(chǎn)[2]。地質(zhì)調(diào)查與現(xiàn)場(chǎng)研究發(fā)現(xiàn)頁(yè)巖儲(chǔ)層中存在較多的天然裂縫，會(huì)影響水力壓裂裂縫的擴(kuò)展路徑，繼而對(duì)裂縫內(nèi)壓裂液和支撐劑的流動(dòng)產(chǎn)生影響[3-5]。因此考慮天然裂縫影響下的頁(yè)巖水力壓裂數(shù)值模擬工作對(duì)頁(yè)巖氣的增產(chǎn)開(kāi)發(fā)具有重要意義[6]。

許多學(xué)者通過(guò)試驗(yàn)手段研究了壓裂裂縫的擴(kuò)展規(guī)律[7-8]。文獻(xiàn)[9-10]通過(guò)室內(nèi)試驗(yàn)和理論分析總結(jié)了水力裂縫與天然裂縫相互作用準(zhǔn)則，提出了水力裂縫貫穿天然裂縫的條件。認(rèn)為水力裂縫與天然裂縫相交作用表現(xiàn)為：水力裂縫的穿越、天然裂縫的張開(kāi)與水力裂縫的阻止3種行為。朱寶存等[11]利用有限元方法分析了地應(yīng)力、天然裂縫對(duì)煤巖水力壓裂起裂壓力的影響。KHOEI等[12]基于擴(kuò)展有限元方法利用力學(xué)響應(yīng)模擬了水力裂縫與天然裂縫交點(diǎn)的關(guān)系。DAHI 等[13]提出了一種基于豐富節(jié)點(diǎn)的數(shù)值模型，研究了不滲透介質(zhì)中水力裂縫被天然裂縫阻止情況下的裂縫交叉點(diǎn)關(guān)系。KESHAVARZI等[14]假設(shè)水力裂縫中保持恒定、均勻的凈壓力，研究了水力裂縫與天然裂縫之間的相互作用。KHOEI等[15]利用一個(gè)交叉富集函數(shù)來(lái)表示裂縫之間的交點(diǎn)，模擬了流體驅(qū)動(dòng)裂縫與摩擦性天然斷層之間的相互作用,利用罰函數(shù)法表示斷層表面之間的摩擦接觸行為，但假設(shè)裂縫內(nèi)的流體壓力是恒定的。文獻(xiàn)[16-17]利用擴(kuò)展有限元方法采用倫肖-波拉德判據(jù)，確定水力裂縫與天然裂縫交點(diǎn)的穿越或阻止行為。研究了摩擦型天然裂縫與膠結(jié)型天然裂縫的差異，發(fā)現(xiàn)含摩擦型天然裂縫的頁(yè)巖儲(chǔ)層的水力壓裂裂縫網(wǎng)絡(luò)更加復(fù)雜。

以上研究為水力裂縫與天然裂縫相交作用提供了理論基礎(chǔ)，但在計(jì)算過(guò)程為簡(jiǎn)化計(jì)算往往忽略裂縫內(nèi)的流體流動(dòng)，只是施加均勻水壓，這與實(shí)際情況不符。同時(shí)對(duì)于水力裂縫與天然裂縫相交后的裂縫擴(kuò)展變化規(guī)律研究較少，天然裂縫對(duì)水力壓裂得到壓裂裂縫影響規(guī)律不明晰。因此筆者在考慮巖石儲(chǔ)層基質(zhì)變形、裂縫內(nèi)的流體流動(dòng)以及天然裂縫與水力裂縫的相交作用情況下，建立了二維水力壓裂裂縫相交擴(kuò)展模型。對(duì)水平應(yīng)力差、天然裂縫相交角與天然裂縫長(zhǎng)度影響裂縫相交前后的水力壓裂裂縫擴(kuò)展問(wèn)題進(jìn)行了研究。

1 水力壓裂裂縫相交擴(kuò)展模型

圖1 水力壓裂裂縫相交擴(kuò)展模型示意

1.1 巖石基質(zhì)變形方程

水力壓裂模型中的巖石基質(zhì)在地應(yīng)力和水力裂縫內(nèi)的流體壓力共同作用下發(fā)生變形，基于線彈性理論方程得巖石基質(zhì)變形控制方程為

?·σ+tb=0

(1)

式中：?為L(zhǎng)aplace算子；σ為應(yīng)力張量；tb為體力張量。

對(duì)應(yīng)的邊界條件如式(2)所示，即

(2)

1.2 水力裂縫內(nèi)滲流方程

假設(shè)水力裂縫內(nèi)流體是不可壓縮牛頓流體，水力裂縫內(nèi)流體流動(dòng)方程由潤(rùn)滑理論定律可得[18]

(3)

式中：q為流體流量；w為水力裂縫張開(kāi)寬度；μ為流體黏度。

設(shè)t為時(shí)間，則水力裂縫內(nèi)流體應(yīng)滿(mǎn)足局部質(zhì)量守恒定律，則有

(4)

將式(3)代入式(4)可得

(5)

同時(shí)滿(mǎn)足整體質(zhì)量守恒方程式(6),即

(6)

壓裂液注入點(diǎn)為定流量條件，注入速率為Q0。假設(shè)裂尖寬度為0，流體無(wú)滯留，則邊界條件為

(7)

1.3 裂縫相交與擴(kuò)展準(zhǔn)則

當(dāng)水力裂縫與天然裂縫發(fā)生相交作用時(shí)，水力裂縫尖端與天然裂縫相交如圖2所示。其中β為水力裂縫與天然裂縫相交角，σH和σh分別表示最大水平地應(yīng)力和最小水平地應(yīng)力，τβ和σβ,n為作用于天然裂縫表面剪切力和法向力。

圖2 水力裂縫與天然裂縫相交應(yīng)力場(chǎng)示意

當(dāng)最大主應(yīng)力σ1大于巖石抗拉強(qiáng)度T0，且天然裂縫摩擦面不發(fā)生滑移，水力裂縫將穿越天然裂縫，在天然裂縫另一側(cè)巖石基質(zhì)中產(chǎn)生一個(gè)新裂縫，其臨界條件[10]為

σ1=T0

(8)

|τβ|

(9)

式中：τβ和σβ,n為作用于天然裂縫表面剪切應(yīng)力和法向應(yīng)力；C為天然裂縫界面黏結(jié)力；μf為天然裂縫摩擦因數(shù)。

裂縫擴(kuò)展行為采用最大周向應(yīng)力準(zhǔn)則進(jìn)行判斷[19]。先利用相互積分法計(jì)算的裂縫尖端應(yīng)力強(qiáng)度因子KⅠ與KⅡ，通過(guò)比較破裂方向等效應(yīng)力強(qiáng)度因子Ke與巖體的斷裂韌度KIC，即可獲得對(duì)應(yīng)的裂縫是否擴(kuò)展。其臨界條件為

(11)

裂縫起裂角θ計(jì)算公式為

(12)

水力壓裂裂縫相交擴(kuò)展模型的求解是考慮巖石基質(zhì)變形與裂縫內(nèi)流體流動(dòng)的非線性流固耦合問(wèn)題。筆者以巖石基質(zhì)變形位移和裂縫內(nèi)流體壓力作為基本變量，采用擴(kuò)展有限元方法和有限元方法分別對(duì)巖石基質(zhì)變形方程和裂縫內(nèi)滲流方程進(jìn)行離散，利用Newton-Raphson方法全耦合求解模型，編制了數(shù)值程序，進(jìn)行數(shù)值模擬研究。

2 模型驗(yàn)證

為驗(yàn)證水力壓裂裂縫相交擴(kuò)展模型的正確性，分別與經(jīng)典KGD模型解析解和Blanton室內(nèi)水力壓裂試驗(yàn)結(jié)果進(jìn)行了對(duì)比。

2.1 與KGD模型對(duì)比

首先對(duì)無(wú)天然裂縫存在時(shí)的二維水力壓裂裂縫擴(kuò)展數(shù)值結(jié)果與經(jīng)典KGD模型解析解[20]比較，數(shù)值計(jì)算輸入?yún)?shù)如下：

巖石基質(zhì)彈性模量E/GPa20巖石基質(zhì)泊松比ν0.2巖石基質(zhì)斷裂韌度KIC/(MPa·m1/2)0.5壓裂液注入速率Q0/(m2·s-1)0.001壓裂液黏度μ/(Pa·s)0.1

注水壓力隨時(shí)間變化關(guān)系如圖3所示，可知數(shù)值計(jì)算結(jié)果與經(jīng)典KGD模型解析解基本吻合，兩者誤差為1.1%，初步驗(yàn)證了模型的正確性。

圖3 注水壓力變化曲線

2.2 與Blanton試驗(yàn)對(duì)比

BLANTON等[9]在室內(nèi)進(jìn)行了三軸應(yīng)力狀態(tài)下預(yù)壓裂材料水力壓裂試驗(yàn)，研究了水力裂縫與天然裂縫在不同水平應(yīng)力差和裂縫相交角下擴(kuò)展情況。為驗(yàn)證本模型考慮水力裂縫與天然裂縫相交作用時(shí)的正確性，進(jìn)行了6種水平應(yīng)力差和4種相交角度組合下的水力壓裂數(shù)值試驗(yàn)對(duì)比Blanton室內(nèi)試驗(yàn)，其中水平應(yīng)力差Δσ=σH-σh。計(jì)算模型如圖4所示，區(qū)域模型具體參數(shù)如下：

巖石基質(zhì)彈性模量E/GPa15巖石基質(zhì)泊松比ν0.25巖石基質(zhì)斷裂韌度KIC/(MPa·m1/2)0.176巖石基質(zhì)抗拉強(qiáng)度T0/MPa3.1天然裂縫摩擦因數(shù)μf0.75天然裂縫黏結(jié)力C/MPa0壓裂液注入速率Q0/(m2·s-1)0.001壓裂液黏度μ/(Pa·s)0.1

圖4 水力裂縫與天然裂縫相交計(jì)算模型

相交角等于60°時(shí)，應(yīng)力差分別為2 MPa和15 MPa情況下的水力裂縫與天然裂縫相交擴(kuò)展情況(圖5)?？芍趹?yīng)力差為2 MPa時(shí)，水力裂縫與天然裂縫相交后，天然裂縫張開(kāi)，水力裂縫沿天然裂縫軌跡擴(kuò)展。此時(shí)天然裂縫段的寬度沿相交點(diǎn)上下出現(xiàn)差異性變化，天然裂縫上段受水力裂縫產(chǎn)生的誘導(dǎo)應(yīng)力影響，其寬度小于天然裂縫下段寬度。當(dāng)應(yīng)力差為15 MPa時(shí)，水力裂縫直接穿越天然裂縫繼續(xù)傳播。11組數(shù)值試驗(yàn)結(jié)果與Blanton試驗(yàn)的對(duì)比結(jié)果見(jiàn)表1，除了一種結(jié)果外，其余都與Blanton試驗(yàn)結(jié)果相同，由此驗(yàn)證了筆者建立的水力壓裂裂縫相交擴(kuò)展模型的有效性。

表1 數(shù)值模型與Blanton試驗(yàn)對(duì)比結(jié)果

圖5 不同水平應(yīng)力差下水力裂縫與天然裂縫相交結(jié)果

3 裂縫相交擴(kuò)展數(shù)值模擬

3.1 水平應(yīng)力差對(duì)裂縫擴(kuò)展的影響

為研究水平應(yīng)力差Δσ對(duì)水力裂縫與天然裂縫相交后裂縫擴(kuò)展的影響，進(jìn)行了當(dāng)水平應(yīng)力差Δσ分別為0、1、2、3 MPa時(shí)的水力壓裂數(shù)值試驗(yàn)。計(jì)算模型，為25 m×25 m的正方形區(qū)域。模型四周受最小水平主應(yīng)力σh和最大水平主應(yīng)力σH作用，沿模型中央點(diǎn)建立xoy坐標(biāo)系(圖6)。初始水力裂縫半

圖6 計(jì)算模型1示意

長(zhǎng)a=2 m，壓裂液注入速率為Q0。2條天然裂縫沿y軸布置在水力裂縫兩側(cè)，其中心點(diǎn)坐標(biāo)(x1、y1)和(x2、y2)分別為(-4，0)和(4，0)，長(zhǎng)度b1=b2=4 m，裂縫相交角β1=β2=90°，模型具體參數(shù)取值如下：

最大水平主應(yīng)力σH/MPa4,5,6,7最小水平主應(yīng)力σh/MPa4巖石基質(zhì)彈性模量E/MPa20巖石基質(zhì)泊松比ν0.25巖石基質(zhì)斷裂韌度KIC/(MPa·m1/2)0.5巖石基質(zhì)抗拉強(qiáng)度T0/MPa2天然裂縫摩擦因數(shù)μf0.69天然裂縫黏結(jié)力C/MPa1壓裂液注入速率Q0/(m2·s-1)0.001壓裂液黏度μ/(Pa·s)0.1

σh=4 MPa、σH=4 MPa,應(yīng)力差為0時(shí)前3個(gè)壓裂步的水力壓裂裂縫擴(kuò)展應(yīng)力云圖如圖7所示。隨著壓裂液的注入，當(dāng)t=2.14 s時(shí)水力裂縫開(kāi)始在水壓和地應(yīng)力共同作用下，沿最大主應(yīng)力方向(水平方向)向兩翼擴(kuò)展，如圖7a、圖7d、圖7g所示。

t=8.35 s時(shí)水力裂縫與天然裂縫相交，此時(shí)儲(chǔ)層應(yīng)力狀態(tài)發(fā)生明顯的變化。如圖7b、圖7d、圖7h所示，x方向應(yīng)力減小，y方向應(yīng)力和剪應(yīng)力增大(本文應(yīng)力為正值時(shí)表示拉應(yīng)力)。水力裂縫兩側(cè)的壓應(yīng)力減小，水力裂縫與天然裂縫相交點(diǎn)處出現(xiàn)x方向應(yīng)力、y方向應(yīng)力與剪應(yīng)力集中區(qū)，同時(shí)天然裂縫兩尖端均有較大的y方向應(yīng)力與剪應(yīng)力產(chǎn)生。圖7c、圖7f、圖7i為水力裂縫與天然裂縫相交后沿天然裂縫兩端擴(kuò)展情況，可知水力裂縫貫通天然裂縫后沿天然裂縫兩端繼續(xù)沿儲(chǔ)層基質(zhì)傳播，裂縫相交點(diǎn)和裂縫尖端均有較大的剪應(yīng)力產(chǎn)生。

4種水平應(yīng)力差水平下水力壓裂裂縫最終擴(kuò)展路徑分布如圖8a所示，A-F為Δσ=3 MPa時(shí)每一破裂步的裂縫擴(kuò)展位置，A點(diǎn)為第1次裂縫擴(kuò)展的裂尖位置。隨著水平應(yīng)力差的減小，與天然裂縫相交后裂縫沿天然裂縫兩端繼續(xù)擴(kuò)展的裂縫路徑由水平方向逐漸向垂直方向傾斜，向最大主應(yīng)力方向偏轉(zhuǎn)。

由圖8b可知當(dāng)水平應(yīng)力差Δσ分別為0、1、2和3 MPa時(shí)，裂縫起裂壓力p0分別為5.93、6.93、7.94和8.93 MPa，可知裂縫的起裂壓力隨水平應(yīng)力差線性增加。水力裂縫未與天然裂縫相交前，注水點(diǎn)壓力從裂縫起裂A點(diǎn)開(kāi)始逐漸減小，水力裂縫與天然裂縫在B點(diǎn)相交,但此時(shí)水力裂縫中的壓裂液還未與天然裂縫直接溝通，C點(diǎn)為水力裂縫溝通天然裂縫并沿天然裂縫兩端傳播得到的裂縫裂尖，此時(shí)注水點(diǎn)凈水壓力快速降低了0.46 MPa。這是由于天然裂縫黏結(jié)力較低，當(dāng)水力裂縫與天然裂縫相交，發(fā)生天然裂縫張開(kāi)行為時(shí)，壓裂液快速貫通天然裂縫，使得注水水壓出現(xiàn)驟降現(xiàn)象。

圖8 水平應(yīng)力差對(duì)裂縫擴(kuò)展以及注水壓力的影響

3.2 天然裂縫相交角與長(zhǎng)度對(duì)裂縫擴(kuò)展的影響

為研究天然裂縫相交角對(duì)水力壓裂裂縫擴(kuò)展的影響，模擬計(jì)算了在最大水平主應(yīng)力σH=4 MPa和最小水平主應(yīng)力σh=2 MPa時(shí)，天然裂縫相交角β分別為45°、60°和75°的水力壓裂裂縫擴(kuò)展數(shù)值實(shí)驗(yàn)。為簡(jiǎn)化計(jì)算，本節(jié)分析水力裂縫單翼擴(kuò)展與單天然裂縫相交擴(kuò)展的對(duì)稱(chēng)模型如圖9所示。其中模型原點(diǎn)位于左下角，注水點(diǎn)坐標(biāo)(0，12.5)水力裂縫沿水平方向x布置，其長(zhǎng)度a=2.5 m天然裂縫長(zhǎng)度b=4 m，天然裂縫中點(diǎn)坐標(biāo)為(8.5,12.5)，模型其他參數(shù)見(jiàn)上文。

圖9 計(jì)算模型2示意

為更加清楚描述天然裂縫相交角大小與天然裂縫長(zhǎng)度對(duì)水力壓裂效果的影響。對(duì)于本文中二維水力壓裂裂縫擴(kuò)展模型，引入壓裂裂縫面積Sf，對(duì)裂縫寬度和長(zhǎng)度進(jìn)行量化描述。其中壓裂裂縫面積Sf計(jì)算公式為

(13)

式中：s=0處為注水點(diǎn)位置，s=l為壓裂裂縫末端尖位置;wi為第i個(gè)裂縫計(jì)算點(diǎn)的裂縫寬度。

由圖10a可知天然裂縫與水力裂縫相交后天然裂縫上端裂尖停止向巖石基質(zhì)傳播，僅沿天然裂縫下端繼續(xù)傳播?？芍烊涣芽p上段靠近水力裂縫，受水力裂縫產(chǎn)生的誘導(dǎo)應(yīng)力影響作用明顯，使得天然裂縫上端裂尖停止向巖石基質(zhì)傳播。而天然裂縫下段離水力裂縫較遠(yuǎn)，受其產(chǎn)生的誘導(dǎo)應(yīng)力影響較小，天然裂縫下端裂尖可繼續(xù)擴(kuò)展。

裂縫寬度與長(zhǎng)度變化規(guī)律如圖10b所示，可知：①水力裂縫A點(diǎn)與天然裂縫相交B點(diǎn)，裂縫寬度逐漸減小，3種相交角度下的裂縫寬度減小趨勢(shì)相同，相交角越大AB段裂縫寬度越大。②天然裂縫BD段，裂縫寬度出現(xiàn)不連續(xù)變化，從相交點(diǎn)B到天然裂縫端D寬度先增大后減小。由4.1節(jié)可知，水力裂縫與天然裂縫相交點(diǎn)周?chē)霈F(xiàn)應(yīng)力集中區(qū)，因此天然裂縫段寬度分布受到限制，使得裂縫相交處的裂縫寬度出現(xiàn)間斷變化，容易造成壓裂液流動(dòng)困難，支撐劑堵塞等情況。由圖10b可知裂縫相交角越大限制作用越強(qiáng)，間斷現(xiàn)象越明顯。③D點(diǎn)之后為沿天然裂縫下端擴(kuò)展的裂縫，其寬度分布先增大后逐漸減小為0，相交角越小裂縫寬度越大。3種相交角45°、60°和75°最終得到的壓裂裂縫長(zhǎng)分別為24.15、23.46和22.12 m，壓裂面積Sf分別為0.312 7、0.301 7和0.283 6 m2。可知隨著天然裂縫相交角度的增大，水力壓裂總裂縫長(zhǎng)度與壓裂面積都隨之減小。

圖10 天然裂縫相交角的影響

同時(shí)研究了天然裂縫長(zhǎng)度對(duì)裂縫擴(kuò)展的影響，分別計(jì)算了天然裂縫長(zhǎng)度為b=6、8 和10 m三種情況下水力壓裂裂縫擴(kuò)展情況，天然裂縫相交角設(shè)置β為60°，其他模型參數(shù)取值與上一算例相同。從圖11a可知3種天然裂縫長(zhǎng)度下，水力裂縫與天然裂縫相交后裂縫擴(kuò)展路徑相似，最終獲得的裂縫長(zhǎng)度分別為25.77、25.78、25.78 m，差別很小，即天然裂縫長(zhǎng)度對(duì)水力壓裂總裂縫長(zhǎng)度影響不大。但天然裂縫長(zhǎng)度影響最終的裂縫寬度分布，天然裂縫長(zhǎng)度的增加，裂縫相交前AB段裂縫寬度逐漸減小，3種天然裂縫長(zhǎng)度下的寬度減小趨勢(shì)相同。同時(shí)隨著天然裂縫的增大天然裂縫BD段寬度限制作用逐漸增大，D點(diǎn)之后的裂縫寬度越小。由圖11b可知3種長(zhǎng)度下最終的壓裂裂縫面積Sf分別為：0.442 5、0.353 1和0.306 5 m2。隨著天然裂縫長(zhǎng)度的增加，壓裂裂縫面積逐漸減小。

圖11 天然裂縫長(zhǎng)度的影響

4 結(jié) 論

在考慮巖石儲(chǔ)層基質(zhì)變形，裂縫內(nèi)的流體流動(dòng)以及天然裂縫與水力裂縫的相交作用情況下建立了二維水力壓裂裂縫相交擴(kuò)展模型，研究分析了水平應(yīng)力差、天然裂縫相交角和天然裂縫長(zhǎng)度對(duì)水力壓裂過(guò)程中發(fā)生水力裂縫與天然裂縫相交擴(kuò)展時(shí)裂縫擴(kuò)展變化的規(guī)律，得出以下結(jié)論：

1)隨著水平應(yīng)力差的線性增加，裂縫起裂注水壓力呈線性增大趨勢(shì)。水平應(yīng)力差對(duì)水力裂縫與天然裂縫相交后的裂縫擴(kuò)展路徑有顯著影響，裂縫逐漸向最大主應(yīng)力方向偏轉(zhuǎn)。

2)水力裂縫與天然裂縫相交時(shí)，注水點(diǎn)水壓出現(xiàn)驟降現(xiàn)象，裂縫相交點(diǎn)處易出現(xiàn)應(yīng)力集中區(qū)，天然裂縫段寬度受到限制作用，裂縫相交處寬度出現(xiàn)間斷現(xiàn)象。

3)天然裂縫相交角影響壓裂裂縫長(zhǎng)度、寬度和面積。天然裂縫相交角越大，裂縫相交點(diǎn)前裂縫寬度越大，而裂縫相交點(diǎn)后的裂縫寬度變的越小，裂縫相交處寬度間斷現(xiàn)象越明顯。最終的壓裂裂縫長(zhǎng)度和壓裂裂縫面積也越小。

4)天然裂縫長(zhǎng)度對(duì)壓裂裂縫長(zhǎng)度影響較小，但對(duì)壓裂裂縫的寬度和面積影響明顯。天然裂縫長(zhǎng)度越長(zhǎng)，裂縫相交點(diǎn)前裂縫寬度越小，裂縫相交點(diǎn)后的裂縫寬度也越小，壓裂裂縫面積越小。