致密儲(chǔ)層CO2壓裂裂縫擴(kuò)展規(guī)律數(shù)值模擬

張健 張國(guó)祥 鄒雨時(shí) 趙文韜 王金意

1 .中國(guó)華能集團(tuán)清潔能源技術(shù)研究院有限公司；2.中國(guó)石油大學(xué)（北京）

目前致密儲(chǔ)層壓裂多采用滑溜水、線性膠等低黏度水基壓裂液，通過(guò)高排量和大液量注入提高儲(chǔ)層改造體積［1］。但水基壓裂液不適用于水敏性、水相圈閉地層，且存在耗水量巨大、化學(xué)試劑污染環(huán)境等問(wèn)題［2］。CO2壓裂具有節(jié)約水資源、降低地層傷害、提高采收率及實(shí)現(xiàn)CO2埋存等優(yōu)點(diǎn)，可避免水基壓裂液的環(huán)境問(wèn)題，具有很好的應(yīng)用前景［3-4］。

國(guó)內(nèi)外學(xué)者對(duì)CO2干法壓裂的研究主要集中于設(shè)備、液體等方面［5］，針對(duì)其造縫機(jī)理研究較為不足。King（1983）［6］研究了低滲透油氣田的 CO2干法壓裂，指出運(yùn)用CO2干法壓裂可消除滲透率和毛細(xì)管力的損害；Settari等（1986）［7］建立了考慮支撐劑沉降和側(cè)向移動(dòng)影響的CO2干法壓裂數(shù)值模型，并通過(guò)礦場(chǎng)實(shí)例驗(yàn)證了其準(zhǔn)確性；王斌等（2004）［8］在擬三維模型基礎(chǔ)上引入簡(jiǎn)化的傳熱模型，建立了CO2干法壓裂裂縫擴(kuò)展數(shù)值模型；蒲春生等（2007）［9］基于斷裂力學(xué)理論，研究了氣體壓裂的裂縫起裂過(guò)程以及增產(chǎn)機(jī)理；孫小輝等（2015）［10］研究了CO2壓裂的溫度場(chǎng)解析模型，但未考慮裂縫擴(kuò)展；陳立強(qiáng)等（2015）［11］研究了 CO2壓裂起裂壓力模型，指出黏度、加壓速率對(duì)起裂壓力影響顯著。

綜上所述，國(guó)內(nèi)外針對(duì)CO2干法壓裂相關(guān)的裂縫擴(kuò)展數(shù)值模擬沒(méi)有考慮壓裂過(guò)程的熱流固耦合問(wèn)題，無(wú)法準(zhǔn)確認(rèn)識(shí)CO2干法壓裂裂縫擴(kuò)展規(guī)律。為此，建立了CO2壓裂熱流固耦合模型，采用有限元法求解，并通過(guò)與解析解對(duì)比驗(yàn)證模型可靠性。最后采用本文方法對(duì)吉林油田致密砂巖CO2干法壓裂的主要參數(shù)進(jìn)行了敏感性分析，并給出了考慮熱流固耦合的CO2干法壓裂的裂縫擴(kuò)展規(guī)律，可為CO2干法壓裂提供設(shè)計(jì)依據(jù)。

1 數(shù)學(xué)模型

CO2壓裂裂縫擴(kuò)展為熱流固耦合過(guò)程，包括注入流體在裂縫和巖石基質(zhì)中的流動(dòng)、巖石變形和CO2壓裂液與地層的傳熱［12］。數(shù)學(xué)模型即針對(duì)這3個(gè)物理過(guò)程建立。注入流體在裂縫和基質(zhì)的流動(dòng)考慮為達(dá)西滲流過(guò)程［13］；壓裂過(guò)程注入排量恒定，注入流體井底溫度恒定［14］；考慮注入CO2為液態(tài)，忽略溫度變化對(duì)流體黏度、相態(tài)的影響［15］。本文應(yīng)力符號(hào)采用彈性力學(xué)的規(guī)定，即取拉應(yīng)力為正、壓應(yīng)力為負(fù)。

1.1 流體流動(dòng)方程

水力壓裂過(guò)程中，裂縫入口不斷有液體注入，從而引起裂縫擴(kuò)展和流體向地層的滲流（濾失）。CO2黏度較低，沿裂縫向地層的濾失較為嚴(yán)重。注入CO2向地層的流動(dòng)方程為［16］

式中，p為流體壓力，MPa；η為壓力傳導(dǎo)系數(shù)，m2/s；t為注入時(shí)間，s。

式（1）描述了流體在地層的濾失，因此可用于分析CO2壓裂過(guò)程流體向地層的壓力擴(kuò)散過(guò)程。

1.2 熱傳導(dǎo)方程

由于注入流體CO2溫度低于地層溫度，水力裂縫擴(kuò)展過(guò)程存在注入流體與地層的熱交換問(wèn)題。根據(jù)傅里葉定律和能量守恒原理，可得熱傳導(dǎo)控制方程為［17］

式中，T為溫度，℃；c為溫度傳導(dǎo)系數(shù)，m2/s。

在CO2干法壓裂過(guò)程中，由于CO2注入流體與地層存在溫度差，會(huì)產(chǎn)生熱應(yīng)力（溫變應(yīng)力）［18］。根據(jù)溫度控制方程，可得到地層任一點(diǎn)的溫度變化，從而可以計(jì)算地層溫度變化引起的熱應(yīng)力?？紤]地層溫度差引起的變形符合線性規(guī)律，則溫度場(chǎng)變化引起的熱應(yīng)力為［18］

式中，K為體積模量，MPa；ΔT為溫度差，℃；β為熱膨脹系數(shù)，℃-1；δij為克羅內(nèi)克函數(shù)。

1.3 巖石變形方程

各向同性介質(zhì)中熱-孔彈性平衡方程為［19］

結(jié)合巖石變形幾何方程，得到固體位移場(chǎng)控制方程為

式中，σij為應(yīng)力分量，MPa；εij、εkk為應(yīng)變分量；G為剪切模量，MPa；v為泊松比；α為Biot系數(shù)。

1.4 擴(kuò)展條件

采用最大拉應(yīng)力準(zhǔn)則判斷巖石是否發(fā)生破壞。當(dāng)計(jì)算單元的應(yīng)力大于抗拉強(qiáng)度時(shí)，巖石發(fā)生破壞。最大拉應(yīng)力準(zhǔn)則為

式中，σ1為最大主應(yīng)力，MPa；σt為巖石抗拉強(qiáng)度，MPa。

1.5 邊界和初始條件

（1）流動(dòng)邊界：外邊界封閉，內(nèi)邊界（注入位置）恒流量。

（2）熱傳導(dǎo)邊界：外邊界絕熱，內(nèi)邊界（注入位置）恒溫。

（3）應(yīng)力邊界：外邊界為遠(yuǎn)場(chǎng)最大最小主應(yīng)力，內(nèi)邊界為注入流體壓力。

（4）初始條件：地層流體壓力為原始地層壓力，溫度為原始地層溫度，初始位移為0。

1.6 模型求解

采用有限元方法求解固體、流體以及熱傳導(dǎo)方程。方程聯(lián)立求解較為復(fù)雜，本文采用弱耦合方法求解［20］，即首先求解流動(dòng)方程和熱傳導(dǎo)方程，再根據(jù)壓力場(chǎng)和溫度場(chǎng)計(jì)算應(yīng)力變化。弱耦合方法在保證一定準(zhǔn)確性的基礎(chǔ)上，減少了大量迭代計(jì)算帶來(lái)的時(shí)間成本問(wèn)題。計(jì)算流程為：首先對(duì)模型設(shè)置初始值，然后依次求解滲流方程、熱傳導(dǎo)方程、固體變形方程，然后判斷單元是否破壞。如果是，則改變單元的物性參數(shù)，如滲透率k、彈性模量E等，并令時(shí)間增加預(yù)設(shè)的增量dt，并再次計(jì)算滲流、熱傳導(dǎo)方程、固體變形方程；否則，則直接令時(shí)間增加預(yù)設(shè)的增量dt，并再次計(jì)算滲流方程、熱傳導(dǎo)方程、固體變形方程。循環(huán)該過(guò)程，直至達(dá)到注入結(jié)束時(shí)間。程序采用Fortran語(yǔ)言編制，具體計(jì)算流程如圖1所示。

圖1 程序計(jì)算框圖Fig. 1 Block diagram of program calculation

2 模型驗(yàn)證

幾何模型如圖2所示。模型X方向長(zhǎng)度為200 m，Y方向長(zhǎng)度為150 m，共劃分為60000個(gè)三角形單元，30351個(gè)節(jié)點(diǎn)。

圖2 幾何模型Fig. 2 Geometric model

為了驗(yàn)證模型的準(zhǔn)確性，將數(shù)值模型模擬結(jié)果與KGD結(jié)果［21］進(jìn)行對(duì)比。模型驗(yàn)證中參數(shù)取值：黏度為 0.2 mPa·s，排量為 0.05 m3/s，彈性模量為 40 GPa，縫高為50 m。圖3為模型計(jì)算裂縫擴(kuò)展動(dòng)態(tài)與解析結(jié)果對(duì)比，可以看出，模型在不同時(shí)刻的縫長(zhǎng)和縫寬計(jì)算結(jié)果與解析解誤差均小于5%，滿足工程計(jì)算精度，表明所建模型可用于裂縫擴(kuò)展的數(shù)值模擬分析。

圖3 所建模型與KGD模型結(jié)果對(duì)比Fig. 3 Result comparison between thermo-hydro-mechanical coupling model and KGD model

3 結(jié)果分析

為分析CO2壓裂裂縫擴(kuò)展規(guī)律，基于第1部分建立的數(shù)值模型對(duì)地層滲透率、注入液體黏度、排量和地層溫度等參數(shù)進(jìn)行敏感性分析。模擬參數(shù)（表1）主要基于吉林油田致密砂巖儲(chǔ)層X(jué)井，滲透率范圍為 0.001～0.100 mD，注入黏度為 0.02～1.00 mPa·s，最大施工排量為5 m3/min，注入液溫度為地表溫度。

表1 輸入?yún)?shù)取值依據(jù)Table 1 Valuing basis of input parameters

3.1 滲透率的影響

圖4為不同滲透率條件下井底壓力變化。圖4顯示，k=0.001 mD時(shí)，地層破裂前增壓速率較大，在約100 s時(shí)刻發(fā)生破裂；k=0.01 mD時(shí)，地層破裂前增壓速率減小，破裂時(shí)刻延后至180 s；k=0.03 mD時(shí)，地層破裂前增壓速率大幅減小，破裂時(shí)刻為780 s；而k=0.1 mD時(shí)，地層破裂前增壓速率很小，注入1000 s內(nèi)并沒(méi)有破裂壓力顯示。注入時(shí)間越長(zhǎng)說(shuō)明注入液量越多，因此，地層滲透率越高，發(fā)生起裂需要流體越多。

圖4 不同滲透率條件的井底壓力隨時(shí)間變化Fig. 4 Change of bottom hole pressure over the time at different permeability

圖5為不同滲透率條件下縫長(zhǎng)和縫寬變化，可以看出，k=0.1 mD時(shí)，裂縫長(zhǎng)度和寬度接近為0，表明裂縫未發(fā)生擴(kuò)展，與井底壓力動(dòng)態(tài)吻合。對(duì)比不同滲透率的縫長(zhǎng)和縫寬變化特征可知，地層滲透率越小，裂縫擴(kuò)展長(zhǎng)度越大，縫寬越大。

圖5 不同滲透率下縫長(zhǎng)和縫寬變化Fig. 5 Change of fracture length and wideth at different permeability

圖6為注入500 s時(shí)刻地層壓力分布。對(duì)比不同滲透率的壓力分布可知，滲透率越大，縫長(zhǎng)越短，沿縫長(zhǎng)的壓力分布范圍越小，而垂直于縫長(zhǎng)方向的壓力擴(kuò)散范圍越大。地層滲透率是影響壓力擴(kuò)散的重要參數(shù)，滲透率越大，注入流體沿地層的擴(kuò)散越快，因此，高滲透地層裂縫起裂時(shí)間會(huì)延后，濾失量多，進(jìn)而縫長(zhǎng)、縫寬減小。

圖6 不同滲透率下的壓力分布Fig. 6 Pressure distribution at different permeability

3.2 黏度的影響

圖7為不同黏度下的井底壓力變化，可以看出，黏度為0.1～1 mPa·s時(shí)，地層破裂前增壓速率基本不變；黏度為0.02 mPa·s，增壓速率略微減小。由此可知，對(duì)于致密儲(chǔ)層，改變注入液體（CO2壓裂液）黏度，不會(huì)顯著影響裂縫起裂時(shí)間。該結(jié)果與Zhang等［22］考慮液體黏度的裂縫起裂數(shù)值模擬結(jié)果一致。

圖7 不同黏度條件的井底壓力隨時(shí)間變化Fig. 7 Change of bottom hole pressure over the time at different viscosity

圖8為不同滲透率條件下縫長(zhǎng)和縫寬變化。對(duì)于CO2壓裂，其他條件相同情況下，注入液體黏度越低，裂縫擴(kuò)展縫長(zhǎng)越短。因此，為解決低黏度液體改造效率過(guò)低問(wèn)題，可適當(dāng)使用增稠劑增加液體黏度，從而提高液體效率和改造效果低黏度CO2濾失速度較大，壓力擴(kuò)散較為顯著，縫內(nèi)壓力耗散較快，因而裂縫長(zhǎng)度較短，而增大黏度會(huì)減小濾失速度，減緩縫內(nèi)壓力耗散，從而增大裂縫長(zhǎng)度。

圖8 不同黏度下縫長(zhǎng)和縫寬變化Fig. 8 Change of fracture length and wideth at different viscosity

3.3 排量的影響

排量對(duì)起裂壓力的影響如圖9所示，可以看出，當(dāng)排量為0.5 m3/min時(shí)，井底增壓速率較小，裂縫在注入約800 s時(shí)才開(kāi)始起裂；隨著排量逐漸增加為1 m3/min、2.5 m3/min、5 m3/min 時(shí)，注入端達(dá)增壓速率依次增大。當(dāng)排量為5 m3/min時(shí)，井底壓力迅速上升，并在注入約50 s后達(dá)到破裂壓力，裂縫起裂與延伸速度較快。綜合不同排量下井底壓力隨時(shí)間變化情況可知，排量越大，裂縫起裂時(shí)間越短。同時(shí)，不同排量下，起裂壓力均為38 MPa，表明注入排量不影響起裂壓力大小。這一現(xiàn)象與起裂壓力理論模型結(jié)果一致，井筒起裂壓力受地應(yīng)力、井筒形狀和地層巖石力學(xué)參數(shù)影響，與注入?yún)?shù)無(wú)關(guān)［23］。

圖9 不同排量下井底壓力隨時(shí)間變化Fig. 9 Change of bottom hole pressure over the time at different displacement

圖10為不同排量下縫長(zhǎng)和縫寬擴(kuò)展動(dòng)態(tài)。結(jié)果顯示，排量為0.5 m3/min時(shí)，縫長(zhǎng)擴(kuò)展緩慢，最大縫長(zhǎng)僅為10 m左右，隨著排量增大，縫長(zhǎng)和縫寬擴(kuò)展顯著，如排量為5 m3/min時(shí)，注入1000 s時(shí)縫長(zhǎng)達(dá)到160 m，縫寬達(dá)到8.2 mm。

圖10 不同排量下縫長(zhǎng)和縫寬變化Fig. 10 Change of fracture length and wideth at different displacement

圖11為不同排量下壓力分布。結(jié)果顯示，增大排量可導(dǎo)致壓力沿縫長(zhǎng)擴(kuò)展，而注入排量較低時(shí)，壓力主要沿垂直于縫長(zhǎng)方向擴(kuò)散。因此，CO2壓裂設(shè)計(jì)應(yīng)在考慮設(shè)備能力的基礎(chǔ)上，盡量提高注入排量，提高裂縫擴(kuò)展長(zhǎng)度和改造體積。在表1參數(shù)條件下，適合采用高于2.5 m3/min的排量進(jìn)行壓裂，低于該值時(shí)地層破裂前增壓速率低，地層不易破裂，形成的裂縫較短且窄，不利于加砂作業(yè)。目前CO2干法壓裂的砂比一般為1%～4%［24］，壓裂后支撐裂縫導(dǎo)流能力較低，導(dǎo)致產(chǎn)能遞減較快。

圖11 不同排量下的壓力分布Fig. 11 Pressure distribution at different displacement

3.4 地層溫度的影響

通過(guò)模擬不同地層溫度下CO2壓裂裂縫擴(kuò)展動(dòng)態(tài)，對(duì)比分析熱應(yīng)力對(duì)裂縫擴(kuò)展的作用。圖12為不同地層溫度條件下井底壓力變化動(dòng)態(tài)。結(jié)果顯示，不考慮溫度場(chǎng)作用時(shí)，地層破裂壓力為39.2 MPa；地層溫度50 ℃時(shí)，破裂壓力為38.3 MPa；地層溫度100℃時(shí)，破裂壓力為36.2 MPa；地層溫度150 ℃時(shí)，破裂壓力為34.5 MPa。因此，相同情況下，CO2壓裂液與地層溫差越大，發(fā)生熱交換量越多，地層熱應(yīng)力越大，破裂壓力越小，同時(shí)破裂時(shí)間也有一定程度減小。

圖12 不同地層溫度下注入壓力隨時(shí)間變化Fig. 12 Change of injection pressure over the time at different reservoir temperature

圖13為不同溫度條件下縫長(zhǎng)和縫寬變化。增大地層溫度與注入液溫差會(huì)促進(jìn)裂縫延伸，提高裂縫寬度。這是因?yàn)闇囟炔钤酱?，地層產(chǎn)生熱應(yīng)力越大，而熱應(yīng)力是拉伸應(yīng)力，有利于減小裂縫延伸阻力，促進(jìn)裂縫延伸。該現(xiàn)象與地?zé)峋畨毫阎袩釕?yīng)力現(xiàn)場(chǎng)相類似［25］，因此，CO2壓裂設(shè)計(jì)不能忽略地層熱應(yīng)力的造縫作用。

圖13 不同地層溫度下縫長(zhǎng)和縫寬變化Fig. 13 Change of fracture length and wideth at different reservoir temperature

圖14為不同地層溫度下，注入CO2導(dǎo)致的地層溫度場(chǎng)分布。

圖14 不同地層初始溫度的溫度場(chǎng)分布Fig. 14 Temperature field distribution correspondind to different initial reservoir temperature

注入CO2導(dǎo)致裂縫附近存在低溫區(qū)，但由于CO2熱傳導(dǎo)范圍較小，低溫區(qū)局限于垂直裂縫距離5 m之內(nèi)。當(dāng)考慮熱應(yīng)力作用時(shí)，裂縫擴(kuò)展縫長(zhǎng)增加。地層溫度越高，注入流體與地層溫度差越大，裂縫長(zhǎng)度增加越大。

4 結(jié)論

（1）CO2干法壓裂設(shè)計(jì)需要考慮注入CO2與地層的熱傳導(dǎo)和熱應(yīng)力作用。地層與注入液的溫度差越大，熱應(yīng)力越大，地層發(fā)生破裂的起裂壓力越小，裂縫延伸越容易。

（2）地層滲透率、注入流體黏度和排量是影響CO2干法壓裂流體濾失的關(guān)鍵因素。滲透率較高儲(chǔ)層（k＞0.01 mD）不適合低黏度CO2壓裂；CO2低黏度特性導(dǎo)致濾失速度較大，壓力擴(kuò)散顯著，裂縫長(zhǎng)度較短、寬度較小。

（3）CO2干法壓裂地層起裂壓力大小不受排量和黏度顯著影響，但排量和黏度會(huì)改變起裂時(shí)間和延伸壓力。