唐煊赫 朱海燕 李奎東

1. “油氣藏地質(zhì)及開發(fā)工程”國(guó)家重點(diǎn)實(shí)驗(yàn)室·成都理工大學(xué) 2. 成都理工大學(xué)能源學(xué)院

3. 中國(guó)石化江漢油田分公司石油工程技術(shù)研究院

0 引言

四川盆地及周緣頁(yè)巖氣資源分布廣泛，其中上奧陶統(tǒng)五峰組—下志留統(tǒng)龍馬溪組頁(yè)巖，為目前頁(yè)巖氣開發(fā)的主要目的層系，具有富有機(jī)質(zhì)頁(yè)巖厚度大、有機(jī)質(zhì)豐度高、有機(jī)質(zhì)熱演化程度高、脆性礦物含量高、地層壓力系數(shù)高、單井產(chǎn)量高、天然裂縫發(fā)育等特點(diǎn)，是我國(guó)頁(yè)巖氣商業(yè)開發(fā)的典范。但該地區(qū)頁(yè)巖氣儲(chǔ)層開發(fā)過程中存在天然裂縫對(duì)該類儲(chǔ)層水力裂縫網(wǎng)絡(luò)擴(kuò)展的影響尚不明確的問題[1-3]。

不僅如此，涪陵頁(yè)巖氣五峰組—龍馬溪組平均埋深超過3 000 m，部分已達(dá)5 000 m，呈現(xiàn)地質(zhì)構(gòu)造復(fù)雜、頁(yè)理發(fā)育、高應(yīng)力和高強(qiáng)度等典型特征，致使儲(chǔ)層施工壓力高達(dá)80～110 MPa。根據(jù)涪陵頁(yè)巖氣田近3年在一期中淺層（垂深2 000～3 500 m）的開發(fā)實(shí)踐發(fā)現(xiàn)：頁(yè)巖氣水平井分段多簇體積壓裂過程中，儲(chǔ)層破裂壓力和裂縫延伸壓力高，常遇到施工壓力超限儲(chǔ)層壓不開、近井地帶裂縫形態(tài)復(fù)雜而引起砂堵的惡性事故。隨著我國(guó)頁(yè)巖氣田逐漸向深層邁進(jìn)（垂深大于3 500 m），地應(yīng)力和巖石強(qiáng)度增大，施工壓力升高，頁(yè)巖儲(chǔ)層近井地帶的多裂縫起裂和擴(kuò)展問題已成為工程技術(shù)人員和科研工作者無(wú)法回避的難題[4-5]。

近年來，國(guó)內(nèi)外學(xué)者針對(duì)頁(yè)巖氣體積壓裂復(fù)雜裂縫的起裂和擴(kuò)展行為開展了大量的室內(nèi)實(shí)驗(yàn)和數(shù)值模擬研究，而數(shù)值模擬研究一般被視作研究頁(yè)巖縫網(wǎng)形成規(guī)律的重要手段，其中常用的復(fù)雜裂縫擴(kuò)展數(shù)值模擬方法包括：離散裂縫網(wǎng)絡(luò)（DFN）[6]、邊界元法（BEM）[7-9]、位移不連續(xù)法（DDA）[10-12]、離散元法（DEM）[13-16]、擴(kuò)展有限元法（XFEM）[17-22]、有限元法（FEM）[23-25]及其發(fā)展方法[26-28]等，均廣泛應(yīng)用于復(fù)雜裂縫擴(kuò)展模擬研究中。目前，相比于其他方法，F(xiàn)EM及其發(fā)展方法（如結(jié)合離散元法或相場(chǎng)法等）模擬裂縫性儲(chǔ)層水力裂縫擴(kuò)展時(shí)優(yōu)勢(shì)顯著：①相比于DFN，能夠考慮天然裂縫、材料各向異性、孔隙壓力等對(duì)裂縫網(wǎng)絡(luò)的影響；②相較于DDM和DDA，適應(yīng)非均質(zhì)、本構(gòu)模型復(fù)雜的材料；③相較于DEM，能直接使用彈性參數(shù)和斷裂力學(xué)參數(shù)；④相較于XFEM，能夠解決大量裂縫的交叉和分岔，而無(wú)需處理復(fù)雜的形函數(shù)；⑤能夠?qū)崿F(xiàn)基質(zhì)孔隙內(nèi)部及其與裂縫系統(tǒng)之間的流固耦合。

在裂縫損傷起裂與擴(kuò)展的表征模型中，黏聚力模型避免了裂尖因子應(yīng)力奇異性導(dǎo)致的不收斂問題，因而得到了廣泛應(yīng)用于裂縫擴(kuò)展過程中的應(yīng)力—應(yīng)變狀態(tài)表征[29]。同時(shí)，水力壓裂過程中，巖石破裂及擴(kuò)展主要依靠注入的壓裂液作外部載荷，需要在黏聚力單元中額外靠攏流體對(duì)裂縫變形的作用，一般將水力裂縫擴(kuò)展區(qū)塊稱作孔壓黏聚力區(qū)域[30]。裂縫擴(kuò)展過程中的損傷力學(xué)行為描述存在不同的線性或非線性模型，但針對(duì)四川盆地中深層頁(yè)巖裂縫擴(kuò)展損傷尚缺乏深入的認(rèn)識(shí)。同時(shí)，國(guó)內(nèi)外在宏觀大型天然裂縫網(wǎng)絡(luò)對(duì)頁(yè)巖水力裂縫擴(kuò)展影響機(jī)理方面缺乏深入研究。

筆者針對(duì)天然裂縫發(fā)育的四川盆地頁(yè)巖氣儲(chǔ)層水力壓裂復(fù)雜裂縫擴(kuò)展問題，基于有限元—離散裂縫網(wǎng)絡(luò)（以下簡(jiǎn)稱FEM-DFN）建立了復(fù)雜裂縫擴(kuò)展的滲流—應(yīng)力—損傷多場(chǎng)耦合模型，并利用解析模型、現(xiàn)場(chǎng)施工注入曲線和現(xiàn)場(chǎng)微地震監(jiān)測(cè)結(jié)果，驗(yàn)證了模型的準(zhǔn)確性和適應(yīng)性。進(jìn)一步以涪陵中深層頁(yè)巖氣藏為目標(biāo)對(duì)象，研究了裂縫頁(yè)巖儲(chǔ)層水力裂縫與天然裂縫的排斥、吸引、分岔和相交競(jìng)爭(zhēng)擴(kuò)展規(guī)律，探討了包括天然裂縫走向、天然裂縫面摩擦系數(shù)、水平應(yīng)力差、施工排量和壓裂液黏度對(duì)水力—天然裂縫交錯(cuò)機(jī)制的影響機(jī)理，分析了多簇水力裂縫在天然裂縫網(wǎng)絡(luò)中的擴(kuò)展規(guī)律。

1 FEM-DFN復(fù)雜裂縫擴(kuò)展理論

FEM-DFN是一類考慮地層天然裂縫網(wǎng)絡(luò)的有限元方法。本文所述頁(yè)巖多裂縫交錯(cuò)擴(kuò)展FEM-DFN方法假設(shè)裂縫沿著網(wǎng)格邊界擴(kuò)展，采用有限單元對(duì)模型進(jìn)行離散，并在所有塊體單元之間插入表征節(jié)理的裂縫單元，從而實(shí)現(xiàn)對(duì)復(fù)雜裂縫網(wǎng)絡(luò)任意擴(kuò)展的模擬。

1.1 巖石介質(zhì)離散與裂縫單元插入

FEM-DFN裂縫擴(kuò)展數(shù)值建模的首要任務(wù)是如何對(duì)連續(xù)介質(zhì)巖體中產(chǎn)生不連續(xù)面（即裂縫）的數(shù)值實(shí)現(xiàn)。在塊體單元之間插入內(nèi)聚力單元，如圖1所示，利用FEM-DFN模擬裂縫擴(kuò)展是將所有表征巖石的塊體（solid）單元離散；各離散塊體單元之間通過插入黏聚力單元表征塊體單元的相互作用以及裂縫變形擴(kuò)展。黏聚力單元以塊體單元的相鄰邊（面）為幾何邊界進(jìn)行構(gòu)建，同時(shí)在上下邊（面）中間插入一層連續(xù)的節(jié)點(diǎn)，用于表征流體（壓裂液）在裂縫內(nèi)的流動(dòng)。黏聚力單元分為完整未損傷單元（未受到裂縫擴(kuò)展影響）、損傷單元（裂縫正在破裂）、完全失效單元（上下巖體已經(jīng)完全斷脫，裂縫面之間無(wú)相互作用力）[31]。

圖1 基于FEM-DFN的網(wǎng)格離散與水力裂縫擴(kuò)展過程圖

對(duì)于脆性較強(qiáng)的頁(yè)巖，其水力壓裂裂縫擴(kuò)展過程包括3個(gè)連續(xù)且互相耦合的過程：①流體（壓裂液）在波及某一未損傷黏聚力單元時(shí)，由于各單元間流動(dòng)層節(jié)點(diǎn)的連續(xù)性，該波及單元中必然至少有1個(gè)流動(dòng)節(jié)點(diǎn)壓力不為0，隨著液體的不斷進(jìn)入，該流動(dòng)節(jié)點(diǎn)壓力不斷上升；②節(jié)點(diǎn)壓力作用在該單元上下邊（面）上，使得裂縫起裂并隨著壓力增大而發(fā)生擴(kuò)展，同時(shí)產(chǎn)生黏聚力，與該黏聚力單元相鄰的2個(gè)巖石固體單元受到施加在裂縫面上的流體壓力以及裂縫擴(kuò)展過程中產(chǎn)生的黏聚力作用，產(chǎn)生（彈性）應(yīng)變；③當(dāng)裂縫上下邊（面）之間的位移達(dá)到某一臨界值，裂縫面之間的黏聚力完全消失，裂縫完全損傷，隨著壓裂液的不斷注入和流動(dòng)，致使受到壓裂液波及進(jìn)而發(fā)生擴(kuò)展的單元增加，裂縫向前擴(kuò)展。

1.2 控制方程

與前文所述FEM-DFN水力裂縫擴(kuò)展計(jì)算過程相對(duì)應(yīng)，分別建立縫內(nèi)流體流動(dòng)、表征裂縫的黏聚力單元變形以及巖石基質(zhì)變形的控制方程。本文中相關(guān)公式單位均基于SI國(guó)際單位制。

1.2.1 縫內(nèi)流體流動(dòng)

本文假設(shè)裂縫內(nèi)流體為不可壓縮牛頓流體，其縫內(nèi)壓裂液流量可采用立方定律描述：

式中q表示縫內(nèi)流量；s表示裂縫到縫口的距離；t表示注液時(shí)間；w表示縫寬；μ表示壓裂液黏度；?p表示縫內(nèi)流體沿縫長(zhǎng)方向壓力梯度。

考慮裂縫內(nèi)流體的質(zhì)量平衡，本文通過潤(rùn)滑方程來描述流體流動(dòng)：

同時(shí)，考慮裂縫壁面的滲透性，由于縫內(nèi)流體與地層孔壓之間存在壓差，因而任意單元內(nèi)水力裂縫壁面向地層滲流可表示為：

式中qleakoff表示縫內(nèi)流體濾失量；cleakoff表示壁面濾失系數(shù)；pf表示縫內(nèi)流體壓力；pp表示儲(chǔ)層孔隙壓力。

通過給定初始注入?yún)?shù)和邊界條件，即可聯(lián)立求解式（1）～（3）。其中，壓裂液流動(dòng)邊界為裂縫起裂位置（注入點(diǎn)）和裂縫尖端，因而邊界條件可表示為：

式中qcrack(t)表示裂縫在任意時(shí)間t的流量；q0(t)表示起裂位置（注入點(diǎn)）在任意時(shí)間t的流量；wtip(t)表示裂縫尖端在任意時(shí)間t的縫寬；qtip(t)表示裂縫尖端在任意時(shí)間t的流量。

1.2.2 裂縫擴(kuò)展非線性變形

根據(jù)斷裂力學(xué)理論[32]，巖體損傷破裂過程中，其裂縫區(qū)域從未損傷到完全損傷可大致分為：完整巖體區(qū)、彈性損傷區(qū)、損傷演化區(qū)、應(yīng)力釋放區(qū)。而根據(jù)Hillerborg理論，裂縫擴(kuò)展區(qū)域稱之為破裂過程區(qū)（Fracture Process Zone，簡(jiǎn)稱FPZ），可以劃分為應(yīng)力釋放區(qū)和非線性損傷區(qū)（將處于破裂過程中的損傷演化區(qū)和彈性損傷區(qū)進(jìn)行整合），非線性損傷區(qū)的損傷應(yīng)力稱為黏聚力。傳統(tǒng)FEM-CZM（Cohesive Zone Model，即黏聚力模型，簡(jiǎn)稱CZM）方法通過引入Cohesive黏聚力區(qū)域模型形成CFEM方法，該方法假設(shè)任意2個(gè)巖石基質(zhì)單元之間的互作用遵循牽引—分離準(zhǔn)則（Traction-separation Law），并以此來描述整個(gè)巖體的損傷及裂縫擴(kuò)展[33]，避免了裂縫尖端的應(yīng)力奇點(diǎn)。FEM-DFN模型采用CZM即黏聚力模型，來描述裂縫破裂損傷過程[34]。而黏聚力模型中最為關(guān)鍵的是損傷黏聚力與裂縫變形之間的關(guān)系。而針對(duì)頁(yè)巖破裂過程中的應(yīng)力—應(yīng)變關(guān)系目前尚存在較大爭(zhēng)議。Zheng等[35]認(rèn)為頁(yè)巖牽引—分離準(zhǔn)則應(yīng)服從非線性關(guān)系，且黏聚力隨著損傷應(yīng)變的增大而不斷線性降低，直至完全破裂失效。然而，從大多數(shù)材料（特別是巖石）的損傷特性來看，裂縫在起裂之前一般會(huì)經(jīng)歷一定的硬化階段，即在損傷初期，黏聚力隨著損傷應(yīng)變的增大而增大，并在達(dá)到起裂應(yīng)變時(shí)才開始發(fā)生實(shí)際的破裂[36]。Li等[37]模擬了層理性頁(yè)巖的裂縫擴(kuò)展，并認(rèn)為其黏聚力—應(yīng)變關(guān)系服從雙線性牽引—分離準(zhǔn)則。而另一方面，根據(jù)Mahabadi等[38]對(duì)包含Ⅰ型、Ⅱ型和Ⅰ-Ⅱ混合型裂縫損傷模式的研究，認(rèn)為在利用FEM-DFN模型描述巖石裂縫擴(kuò)展時(shí)，非線性牽引—分離準(zhǔn)則（特別是損傷屈服階段）能夠較好地通過抗拉強(qiáng)度測(cè)試和單軸壓縮測(cè)試進(jìn)行驗(yàn)證[39]。Dahi等[40]提出了線性硬化—非線性損傷屈服相結(jié)合的黏聚力模型，較好地體現(xiàn)了雙線性和非線性牽引—分離準(zhǔn)則的特征。因此，筆者將頁(yè)巖的牽引—分離準(zhǔn)則描述為裂縫破裂應(yīng)力—應(yīng)變應(yīng)經(jīng)歷3個(gè)階段，即：初始線性硬化、非線性損傷屈服、完全失效。

筆者提出的非線性牽引—分離準(zhǔn)則，除了將線性硬化—非線性損傷屈服相結(jié)合外，同時(shí)還包括對(duì)Ⅰ型、Ⅱ型和Ⅰ-Ⅱ混合型裂縫損傷模式的表征（圖2）。對(duì)于Ⅰ型裂縫（即拉伸或壓縮），如圖2-a和公式（6）所示，其初始硬化過程采用Hillerborg等[31]提出的線性硬化，當(dāng)裂縫面法向黏聚力隨著應(yīng)變線性增大到頁(yè)巖的抗拉強(qiáng)度（σt）時(shí)，此時(shí)的裂縫應(yīng)變?yōu)槌跏紦p傷應(yīng)變（δn0），巖體開始發(fā)生損傷屈服。而損傷屈服階段則采用Munjiza[41]的非線性模型進(jìn)行描述，但與Munjiza[41]的不同之處在于，本文模型假設(shè)2個(gè)裂縫面間法向位移增長(zhǎng)到裂縫失效臨界法向位移（δnf）后，裂縫面之間將不存在法向黏聚力。

圖2 非線性牽引—分離準(zhǔn)則圖

式中δn表示法向位移；表示起裂時(shí)的法向位移；σt0表示巖石抗拉強(qiáng)度；δnf表示裂縫失效臨界法向位移。

對(duì)于Ⅱ型裂縫（即剪切），本文采用Ida[42]提出的滑移—弱化（slip-weakening）模型，如圖2-b和公式（7）所示，與Ⅰ型破裂模式類似，初始剪切黏聚力也隨著切向應(yīng)變?cè)龃蠖€性增大，直至達(dá)到頁(yè)巖的剪切強(qiáng)度（τs）即開始剪切損傷屈服。由于開始損傷后，切向黏聚力始終與剪切強(qiáng)度有關(guān)，因此損傷屈服過程也服從Mohr–Coulomb準(zhǔn)則，如公式（8）所示。當(dāng)巖石完全切向失效后，裂縫面之間完全斷脫，與天然裂縫面一致，采用Coulomb剪切定律[43]描述裂縫面之間的摩擦效應(yīng)。

式中δs表示切向位移；δs0表示起裂切向位移；τs表示巖石剪切強(qiáng)度；δsf表示裂縫失效臨界法向位移；τf表示裂縫切向失效后的摩擦阻力。

其中，抗剪強(qiáng)度τs為：

式中c表示巖石的內(nèi)聚力；?i表示巖石的內(nèi)摩擦角；μf表示天然裂縫面的摩擦系數(shù)。

摩擦阻力τf為：

巖石的裂縫擴(kuò)展除需要考慮為Ⅰ型裂縫或Ⅱ型裂縫外，還涉及Ⅰ-Ⅱ型混合裂縫模式，且損傷演化階段更傾向于非線性軟化，混合模式判斷準(zhǔn)則如圖2-c所示。

由于Ⅱ型裂縫模式的剪切強(qiáng)度計(jì)算與法向應(yīng)力相關(guān)，本文中混合裂縫模式通過考慮將法向位移和切向位移進(jìn)行耦合分析[44]，如式（10）所示：

完全失效臨界法向/切向位移是準(zhǔn)確約束圖2-a、b所示曲線的關(guān)鍵，為此，本文引入斷裂能[43]來表征裂縫的完全失效；同時(shí)，引入二次能量準(zhǔn)則[45]來約束應(yīng)變能的增加，作為圖2-c所示混合模式裂縫完全失效的判斷準(zhǔn)則。

在裂縫擴(kuò)展過程中，巖石除了受到原有地應(yīng)力的作用外，其外表面（即裂縫面）同時(shí)還受到裂縫擴(kuò)展過程中產(chǎn)生的黏聚力和流體施加的壓力，因而將發(fā)生相應(yīng)的變形。本文假設(shè)巖石在裂縫擴(kuò)展過程中為各向同性彈性變形，并將巖石基質(zhì)考慮為多孔介質(zhì)，并服從虎克定律。另一方面，由于頁(yè)巖氣藏屬于特低滲透致密巖層，壓裂液的濾失通道主要是張開的天然裂縫或水力裂縫[46-48]，且本文針對(duì)壓裂裂縫擴(kuò)展及延伸的施工時(shí)間尺度范圍內(nèi)開展研究，因而壓裂液向地層濾失導(dǎo)致的地層孔隙壓力變化可以忽略。模型初始狀態(tài)下無(wú)剪切應(yīng)變，且外邊界條件為位移邊界。

1.2.3 FEM-DFN在時(shí)間及空間上的離散

在任意嵌入了水力裂縫的巖體區(qū)域內(nèi)，必然包含巖體受到外載荷的外邊界，同時(shí)又會(huì)有水力裂縫在巖體內(nèi)部形成弱面，流體壓力施加在區(qū)域內(nèi)部裂縫表面形成內(nèi)載荷的內(nèi)邊界。同時(shí)，由于縫內(nèi)流體流動(dòng)是平板流動(dòng)，流體壓力和流量均為與時(shí)間有關(guān)的變量。因而，在不考慮模型所受體力的情況下，本文模型中縫內(nèi)流動(dòng)、巖石破裂、巖石變形均是與時(shí)間有關(guān)的瞬態(tài)過程，通過三角形固體單元和無(wú)厚度的黏聚力平面單元進(jìn)行離散，并從動(dòng)力學(xué)角度建立整體有限元方程：

對(duì)上述方程采用顯式積分法，對(duì)巖石變形和裂縫擴(kuò)展方程采用向前時(shí)間差分。

式中M表示集中質(zhì)量矩陣；C表示阻尼矩陣；K表示剛度矩陣；x表示全局節(jié)點(diǎn)位移；Fext表示區(qū)域所受外載荷；Fc表示裂縫擴(kuò)展過程中的黏聚力；Fint表示區(qū)域內(nèi)部裂縫表面所受載荷；H表示集中容積矩陣；p表示節(jié)點(diǎn)壓力向量；Iint和Iext分別表示區(qū)域內(nèi)部和外部流量載荷。

2 模型驗(yàn)證

2.1 室內(nèi)實(shí)驗(yàn)驗(yàn)證

為了測(cè)試模型在頁(yè)巖儲(chǔ)層中的適用性，引入現(xiàn)有公開文獻(xiàn)中的頁(yè)巖室內(nèi)水力壓裂大型物理模擬實(shí)驗(yàn)數(shù)據(jù)，并建立相同幾何及力學(xué)的FDEM模型，模擬實(shí)驗(yàn)載荷及注入條件下的裂縫擴(kuò)展情況，對(duì)比室內(nèi)實(shí)驗(yàn)結(jié)果。本文參考文獻(xiàn)[49]為針對(duì)包含層理和天然裂縫的頁(yè)巖水力壓裂大型物理模擬實(shí)驗(yàn)。實(shí)驗(yàn)巖樣尺寸為300 mm×300 mm×300 mm，在平行于層理方向上：彈性模量17.2 GPa，泊松比0.183，抗拉強(qiáng)度1.35 MPa，內(nèi)聚力6.15 MPa；而垂直于層理方向上：彈性模量17.3 GPa，泊松比0.183，抗拉強(qiáng)度3.85 MPa，內(nèi)聚力11.41 MPa。并選取該文獻(xiàn)中變流量注入實(shí)驗(yàn)組14號(hào)：注入流量為50 cm3/min，注入液體黏度為65.0 mPa·s。注入方式為階梯式變流量注入，注入速率為0.05×(n+1) mL/s，三向地應(yīng)力分別為：σV＝ 25 MPa，σH＝ 15 MPa，σh＝ 12 MPa。

實(shí)驗(yàn)結(jié)束后分別在平行和垂直于層理方向?qū)r樣進(jìn)行CT切片掃描，如圖3-a所示。同時(shí)，分別從CT1切面和CT2切面建立FDEM水力壓裂裂縫擴(kuò)展模型（圖3-b～c），在CT1切面上建立裂縫擴(kuò)展幾何模型，通過幾何反演得到巖樣的層理分布，并嵌入模型中。其中，天然裂縫無(wú)膠結(jié)，層理為無(wú)膠結(jié)或弱膠結(jié)。

圖3 14號(hào)巖樣實(shí)驗(yàn)后切面觀測(cè)面示意圖

從圖4、5可以看出，水力裂縫在擴(kuò)展過程中基本上沿著最大（水平）主應(yīng)力方向擴(kuò)展，在與無(wú)膠結(jié)的天然裂縫相遇時(shí)，僅僅將誘發(fā)了天然裂縫發(fā)生剪切滑移，然后就直接穿透裂縫。

圖4 CT1視角下的模型對(duì)比驗(yàn)證圖

2.2 現(xiàn)場(chǎng)施工驗(yàn)證

通過涪陵頁(yè)巖氣水平井壓裂施工曲線驗(yàn)證計(jì)算得到的破裂壓力是否合理、通過微地震監(jiān)測(cè)解釋結(jié)果來驗(yàn)證計(jì)算得到儲(chǔ)層改造范圍是否合理。

圖5 CT2視角下的模型對(duì)比驗(yàn)證圖

2.2.1 基于儲(chǔ)層地質(zhì)特征與工程條件的基準(zhǔn)模型

四川盆地涪陵頁(yè)巖氣某典型水平井區(qū)的天然裂縫主要以高角度裂縫（裂縫傾角均大于87°，本文統(tǒng)一處理為90°）和水平層理為主，而層理主要為充填裂縫，因此，本文主要考慮水力裂縫在擴(kuò)展過程中主要與高角度非膠結(jié)天然裂縫互相溝通形成復(fù)雜裂縫網(wǎng)絡(luò)。根據(jù)該區(qū)塊的前期露頭觀測(cè)和成像測(cè)井分析結(jié)果顯示[50]，該區(qū)塊頁(yè)巖儲(chǔ)層天然裂縫走向在3個(gè)角度左右較為集中，分別為：215°（占比約38%）、289°（占比約32%）、360°（占比約28%），如圖6所示。同時(shí)，結(jié)合該井所處區(qū)域的地震數(shù)據(jù)體反演統(tǒng)計(jì)與分析結(jié)果可知，建模井段附近裂縫的平面分布和縫長(zhǎng)則服從對(duì)數(shù)正態(tài)分布[51-54]。

圖6 高角度天然裂縫走向與分布圖

該頁(yè)巖氣水平井水力壓裂采用分段實(shí)施壓裂施工，單段內(nèi)進(jìn)行了以3簇為主的多簇射孔。因此，在表1的地質(zhì)力學(xué)和注入?yún)?shù)基礎(chǔ)上，選擇其中某3簇射孔壓裂段，建立近井地帶多簇同步壓裂裂縫擴(kuò)展模型。如圖7所示，模型的幾何尺寸為200 m×200 m，并沿井筒方向設(shè)置3個(gè)射孔簇作為注入點(diǎn)，各簇間距30 m。同時(shí)，在模型中以縫面密度約為0.007 5/m2，嵌入該儲(chǔ)層3種不同走向天然裂縫，天然裂縫長(zhǎng)度分布范圍為15～30 m，并結(jié)合該區(qū)域螞蟻體數(shù)據(jù)進(jìn)行修正。

圖7 裂縫性地層多簇壓裂幾何與DFN圖

表1 模型的地質(zhì)力學(xué)參數(shù)表

由前期鉆完井資料可知，該區(qū)塊儲(chǔ)層的最大水平主應(yīng)力方向大約為西—東向，最小水平主應(yīng)力方向大約為北—南向，水平井筒沿著最小水平主應(yīng)力方向分布，因此，設(shè)置模型的x方向?yàn)樽畲笏街鲬?yīng)力方向，y方向?yàn)樽钚∷街鲬?yīng)力方向，水平井筒沿著模型y軸分布，射孔點(diǎn)（壓裂液注入點(diǎn)）處于井筒中央，模型幾何邊界為不可滲透邊界。同時(shí)，根據(jù)前期地質(zhì)力學(xué)資料，該區(qū)塊地質(zhì)力學(xué)參數(shù)如表1所示。

2.2.2 微地震監(jiān)測(cè)結(jié)果驗(yàn)證

裂縫擴(kuò)展模擬的最終形態(tài)如圖8-a所示，圖中裂縫包括天然裂縫網(wǎng)絡(luò)、張性水力裂縫和剪切水力裂縫，反映水力裂縫擴(kuò)展過程中與天然裂縫溝通形成的改造區(qū)域形態(tài)和大致范圍，Lf表示總裂縫長(zhǎng)度，Lt表示拉伸裂縫長(zhǎng)度，Ls表示剪切裂縫長(zhǎng)度；而圖8-b為該井段壓裂過程中井下微地震監(jiān)測(cè)反演得到的微地震事件點(diǎn)，一定程度上反映實(shí)際壓裂改造區(qū)域形態(tài)和大致范圍。對(duì)比圖8-a、b可知，模擬裂縫形態(tài)和現(xiàn)場(chǎng)監(jiān)測(cè)微地震事件點(diǎn)的形態(tài)都呈現(xiàn)出左側(cè)裂縫在區(qū)域長(zhǎng)度和寬度上均比右側(cè)裂縫擴(kuò)展更為充分、左側(cè)改造區(qū)各簇?cái)U(kuò)展相對(duì)均勻、右側(cè)改造區(qū)上下兩簇改造范圍大于中間簇等相似特征。同時(shí)，模擬裂縫改造左側(cè)區(qū)域長(zhǎng)度約87 m，右側(cè)約65 m；而現(xiàn)場(chǎng)監(jiān)測(cè)結(jié)果左側(cè)區(qū)域長(zhǎng)度約80 m，右側(cè)約68 m。說明模擬裂縫能夠較好地反映實(shí)際壓裂改造范圍。

圖8 基于現(xiàn)場(chǎng)微地震監(jiān)測(cè)的多簇裂縫擴(kuò)展模擬結(jié)果驗(yàn)證圖

3 復(fù)雜裂縫擴(kuò)展規(guī)律

3.1 水力裂縫與天然裂縫互作用機(jī)理

為了研究水力壓裂裂縫在裂縫性頁(yè)巖儲(chǔ)層中的擴(kuò)展規(guī)律，必須首先明確在該類地層條件下，水力裂縫與天然裂縫的互作用機(jī)理。在圖7模型基礎(chǔ)上進(jìn)行簡(jiǎn)化，僅在考慮一條天然裂縫的情況下，探討水力裂縫在擴(kuò)展過程中主要與高角度非膠結(jié)天然裂縫互作用機(jī)理。如圖9所示，將走向?yàn)?15°（逼近角55°，為該地層天然裂縫走向分布的峰值）的天然裂縫嵌入到水力—天然裂縫互作用模型中。

模型計(jì)算結(jié)果如圖9～11所示，包括水力裂縫擴(kuò)展過程及其巖石變形，以及注入壓力隨時(shí)間變化的時(shí)程曲線。水力裂縫在含有天然裂縫的地層中擴(kuò)展，大致可以描述為如下幾個(gè)過程。

圖9 水力裂縫與單條天然裂縫互作用數(shù)值模型圖

1）水力裂縫擴(kuò)展過程中，注入壓力隨著液體的注入不斷升高，直到在注入時(shí)間71 s時(shí)達(dá)到破裂壓力開始起裂，并迅速降低至閉合壓力，如圖10中第Ⅰ階段所示；縫內(nèi)壓力隨著壓裂液的注入而不斷增大，并通過施加在裂縫表面使得裂縫兩側(cè)地層受到擠壓而發(fā)生變形，同時(shí)波及前方巖石使其發(fā)生斷裂向前擴(kuò)展。

圖10 注入壓力隨模擬時(shí)間變化曲線圖

2）注入時(shí)間為130 s時(shí)，天然裂縫起裂，對(duì)應(yīng)壓力時(shí)程曲線中出現(xiàn)的壓力突降，如圖10中第Ⅱ階段所示；當(dāng)水力裂縫與天然裂縫交匯時(shí)，如圖10所示，天然裂縫開始起裂，且水力裂縫進(jìn)入天然裂縫南側(cè)分支裂縫使其優(yōu)先開裂，且天然裂縫北側(cè)分支發(fā)生一定程度的剪切。

3）當(dāng)水力裂縫結(jié)束在天然裂縫中的擴(kuò)展重新轉(zhuǎn)向時(shí)，需要克服更大的應(yīng)力，注入壓力逐漸上升，如圖10中第Ⅲ階段所示；水力裂縫沿著天然裂縫南側(cè)分支不斷擴(kuò)展，且由于水力裂縫內(nèi)壓力主要用于使得天然裂縫向前擴(kuò)展，因而水力主裂縫的東側(cè)分支暫停向前擴(kuò)展；在天然裂縫內(nèi)擴(kuò)展時(shí)，南側(cè)分支兩側(cè)巖體所受壓應(yīng)力存在差異：①在靠近注入點(diǎn)一側(cè)巖體同時(shí)受到原水力主裂縫擴(kuò)展時(shí)對(duì)其的擠壓和天然裂縫擴(kuò)展時(shí)對(duì)其的擠壓（雙向擠壓）；②在靠近模型邊界一側(cè)的巖體僅受到天然裂縫擴(kuò)展對(duì)其的擠壓（單向擠壓），且由于裂尖的低應(yīng)力區(qū)疊加作用，使得該側(cè)的應(yīng)力遠(yuǎn)小于靠近模型邊界一側(cè)（圖11）。

圖11 水力—天然裂縫互作用過程中的相遇與重新轉(zhuǎn)向擴(kuò)展圖

4）當(dāng)右側(cè)水力裂縫擴(kuò)展至模型的不可滲透邊界時(shí)，注入壓力僅用于擴(kuò)展左側(cè)水力裂縫的張性破裂，因此會(huì)經(jīng)歷一段穩(wěn)定壓力階段，直到左側(cè)裂縫也擴(kuò)展到不可滲透邊界，液體不斷注入導(dǎo)致已擴(kuò)展裂縫內(nèi)的凈壓力不斷上升，所示圖10中第Ⅳ階段所示；在擴(kuò)展到其南側(cè)分支結(jié)束后，由于天然裂縫兩側(cè)的應(yīng)力差和最大水平主應(yīng)力的作用，迫使其重新轉(zhuǎn)向最大水平主應(yīng)力方向擴(kuò)展。此時(shí)，裂縫兩側(cè)所受應(yīng)力狀態(tài)相近，縫內(nèi)壓力迫使兩側(cè)裂縫均勻向前擴(kuò)展。當(dāng)裂縫擴(kuò)展至模型的不滲透邊界，縫內(nèi)凈壓力升高，導(dǎo)致天然裂縫北側(cè)分支的破裂行為從剪切變?yōu)槔臁?/p>

3.2 施工參數(shù)與地質(zhì)力學(xué)參數(shù)的影響

地質(zhì)力學(xué)特征與壓裂施工參數(shù)均會(huì)影響頁(yè)巖氣儲(chǔ)層水力壓裂裂縫擴(kuò)展，需要基于現(xiàn)場(chǎng)施工實(shí)際情況（以表1中參數(shù)為基礎(chǔ)），開展注入排量、壓裂液黏度、水平應(yīng)力差、裂縫面摩擦系數(shù)和天然裂縫走向等不同參數(shù)條件下，天然—水力裂縫互作用機(jī)理如表2所示。

表2 注入?yún)?shù)與地質(zhì)力學(xué)參數(shù)表

通過對(duì)不同的影響因素進(jìn)行模擬計(jì)算，發(fā)現(xiàn)在頁(yè)巖水力裂縫擴(kuò)展過程中，其與天然裂縫的交錯(cuò)機(jī)制主要為沿著天然裂縫擴(kuò)展到其末端再進(jìn)行轉(zhuǎn)向和直接穿透天然裂縫2種交錯(cuò)機(jī)制。因此，針對(duì)不同地質(zhì)力學(xué)或注入?yún)?shù)等因素下天然裂縫對(duì)水力裂縫擴(kuò)展影響機(jī)理的關(guān)注重點(diǎn)應(yīng)該著眼于：隨著影響因素的增大或減小，水力—天然裂縫交錯(cuò)機(jī)制在何時(shí)發(fā)生轉(zhuǎn)變，即沿天然裂縫擴(kuò)展和穿透天然裂縫之間的閾值。

3.2.1 排量、黏度、應(yīng)力差、摩擦系數(shù)的影響

圖12為注入壓力隨注入排量、壓裂液黏度、水平應(yīng)力差和裂縫面摩擦系數(shù)的變化情況，并結(jié)合對(duì)應(yīng)參數(shù)下模擬得到的水力—天然裂縫互作用情況進(jìn)行討論。

如圖12-a所示，在注入排量較低時(shí)，水力—天然裂縫交錯(cuò)機(jī)制為沿著天然裂縫誘導(dǎo)轉(zhuǎn)向。隨著注入排量的不斷增大，注入壓力不斷上升，這是由于天然裂縫為無(wú)膠結(jié)裂縫，且天然裂縫表征黏聚力單元的破裂力學(xué)參數(shù)小于巖石的破裂參數(shù)，壓裂液不斷向天然裂縫內(nèi)濾失，因此穿透天然裂縫所需要克服的阻力大于壓開天然裂縫的阻力。隨著排量的不斷上升，水力裂縫兩側(cè)巖石受到垂直于裂縫面方向擠壓，從而產(chǎn)生橫向膨脹變形，進(jìn)而使得天然裂縫面法向壓力增大；在排量達(dá)到26 m3/min時(shí)，水力裂縫壓開天然裂縫所受阻力大于穿透天然裂縫所受阻力，水力裂縫穿透天然裂縫，此時(shí)注入壓力出現(xiàn)陡增現(xiàn)象。

圖12 不同因素對(duì)水力—天然裂縫交錯(cuò)機(jī)制的影響圖

如圖12-b所示，在壓裂液黏度較低時(shí)，水力—天然裂縫交錯(cuò)機(jī)制為沿著天然裂縫誘導(dǎo)轉(zhuǎn)向。隨著黏度的不斷增大，注入壓力不斷上升；在黏度達(dá)到40 mPa·s時(shí)，水力裂縫穿透天然裂縫，且此時(shí)注入壓力陡增，這是由于較大的黏度將使得剪切應(yīng)力增大。同時(shí)，隨著黏度增大，壓裂液向天然裂縫濾失減少，裂縫內(nèi)流體壓力增加變慢，剪切滑移阻力增大，使得水力裂縫在遭遇天然裂縫時(shí)，更容易穿過天然裂縫。

如圖12-c所示，當(dāng)水平應(yīng)力差較低時(shí)，水力裂縫擴(kuò)展區(qū)域附近的巖石受擠壓導(dǎo)致應(yīng)力差進(jìn)一步減小，水力裂縫一旦與天然裂縫相遇，水力裂縫被誘導(dǎo)沿著天然裂縫擴(kuò)展。隨著水平應(yīng)力差的不斷增大，注入壓力不斷下降（裂縫起裂損傷的臨界壓力降低）；在水平應(yīng)力差達(dá)到18 MPa時(shí)，水力裂縫穿透天然裂縫，且此時(shí)注入壓力下降速度增大，這是由于水平應(yīng)力差越大，水力裂縫壁面受到的相對(duì)阻力越小，越有利于向前擴(kuò)展，且模型中是通過使最小水平主應(yīng)力降低更多來增大應(yīng)力差，也有利于裂縫向前擴(kuò)展。

如圖12-d所示，在天然裂縫面摩擦系數(shù)較低時(shí)，水力—天然裂縫交錯(cuò)機(jī)制為沿著天然裂縫誘導(dǎo)轉(zhuǎn)向。隨著摩擦系數(shù)的不斷增大，注入壓力不斷上升；且在摩擦系數(shù)達(dá)到0.8時(shí)，水力裂縫將直接穿透天然裂縫，且此時(shí)注入壓力上升速度增大，這是由于摩擦系數(shù)越大，水力裂縫在與天然裂縫相遇時(shí)，雖然容易受其誘導(dǎo)擴(kuò)展，但是較大的摩擦系數(shù)增大了其剪切一側(cè)裂縫的難度，此時(shí)裂尖壓力不斷上升，直接將裂縫前方的黏聚力單元撕開。因此，準(zhǔn)確測(cè)量天然裂縫面的摩擦系數(shù)對(duì)評(píng)價(jià)裂縫性儲(chǔ)層壓裂裂縫復(fù)雜度至關(guān)重要。

3.2.2 天然裂縫走向的影響

圖13-a為天然裂縫與最大水平主應(yīng)力夾角（即逼近角）19°時(shí)的裂縫擴(kuò)展情況。與逼近角55°時(shí)不同，水力裂縫并非在與天然裂縫交匯出進(jìn)入天然裂縫，而是受到天然裂縫誘導(dǎo)提前轉(zhuǎn)向匯入天然裂縫；在與天然裂縫交錯(cuò)后，其擴(kuò)展規(guī)律則與逼近角55°時(shí)類似。

圖13 與最大水平主應(yīng)力夾角19°和90°時(shí)的裂縫擴(kuò)展情況圖

圖13-b為天然裂縫與最大水平主應(yīng)力夾角90°時(shí)的裂縫擴(kuò)展情況。與逼近角為19°和55°時(shí)均不同，水力裂縫在與天然裂縫交匯后沿著天然裂縫兩側(cè)分支均發(fā)生剪切作用，這是由于最大水平主應(yīng)力垂直于裂縫面，裂縫面受到壓應(yīng)力最大，裂縫發(fā)生拉伸變形的難度最高；但當(dāng)左側(cè)裂縫擴(kuò)展至不可滲透邊界時(shí)，縫內(nèi)凈壓力隨著液體的泵入不斷上升，使得水力—天然裂縫交匯處壓力能夠克服最大水平主應(yīng)力，迫使原本為剪切狀態(tài)的天然裂縫發(fā)生拉伸斷裂。

3.3 多簇壓裂復(fù)雜裂縫擴(kuò)展規(guī)律

基于水力—天然裂縫互作用機(jī)理認(rèn)識(shí)，針對(duì)圖7所示模型的實(shí)際多簇射孔壓裂裂縫擴(kuò)展規(guī)律模擬（圖14），分別為模型位移云圖和裂縫幾何形態(tài)圖。裂縫幾何形態(tài)圖中區(qū)分了張性裂縫和剪切裂縫。在裂縫擴(kuò)展初期，由于第1簇和第3簇附近大量天然裂縫的存在，該兩簇水力裂縫一開始即被天然裂縫誘導(dǎo)，而第2簇裂縫則沿著最大水平主應(yīng)力方向擴(kuò)展。此外，由于水力裂縫擴(kuò)展過程中產(chǎn)生的誘導(dǎo)應(yīng)力場(chǎng)，導(dǎo)致附近的巖石發(fā)生變形錯(cuò)動(dòng)，進(jìn)而使得部分未連通的天然裂縫也發(fā)生了一定程度的剪切。

圖14 裂縫性地層井周多簇裂縫擴(kuò)展圖

如圖14-a～d所示，在密集分布的天然裂縫作用下，第3簇右翼裂縫被誘導(dǎo)至與第2簇右翼裂縫相交匯，導(dǎo)致第3簇右翼裂縫沿最大水平主應(yīng)力擴(kuò)展范圍僅10.0 m。同時(shí)，由于兩簇裂縫的交匯形成一條優(yōu)勢(shì)主裂縫，使得該縫內(nèi)流量增大，更加易于穿透天然裂縫形成張性水力裂縫，進(jìn)而導(dǎo)致剪切裂縫占比降低。此外，由于第1簇和第3簇裂縫擴(kuò)展過程中導(dǎo)致的儲(chǔ)層巖石變形以及誘導(dǎo)應(yīng)力，使得第2簇裂縫兩側(cè)巖石變形增大，縫寬受到限制，遠(yuǎn)小于另外兩簇裂縫。因此，簇間距優(yōu)化重點(diǎn)考慮以下兩點(diǎn)：①避免兩側(cè)裂縫對(duì)中間裂縫的限制，盡量使得各簇裂縫均勻擴(kuò)展；②在高密度天然裂縫地層，避免各簇水力主裂縫過早地交匯形成優(yōu)勢(shì)主裂縫，降低了壓裂裂縫復(fù)雜度。

隨著注液時(shí)間的增加，由圖14-e～f可知，第1簇裂縫向右擴(kuò)展至密集天然裂縫帶后，在天然裂縫的強(qiáng)誘導(dǎo)作用下，大量壓裂液向天然裂縫帶濾失，導(dǎo)致該區(qū)域形成剪切裂縫帶并限制了第1簇裂縫擴(kuò)展；由于天然—水力裂縫夾角作用，當(dāng)?shù)?簇右翼裂縫與第2簇右翼裂縫發(fā)生交匯后形成的裂縫，同時(shí)被第2簇右側(cè)和兩簇連通的兩個(gè)90°夾角天然裂縫所遮擋，使得右側(cè)裂縫難以向前擴(kuò)展。此時(shí)，第2簇和第3簇內(nèi)的縫內(nèi)凈壓力主要用于其左翼裂縫的擴(kuò)展，使得第2簇水力裂縫能夠連續(xù)穿透多條天然裂縫。

總的來說，相比于不含天然裂縫的儲(chǔ)層中以地應(yīng)力作為裂縫擴(kuò)展主控因素，裂縫性頁(yè)巖儲(chǔ)層中的天然裂縫（包括天然裂縫的產(chǎn)狀、發(fā)育程度、裂縫面等）對(duì)水力裂縫擴(kuò)展的影響比地應(yīng)力更為顯著。

4 結(jié)論

1）當(dāng)水力裂縫擴(kuò)展至天然裂縫時(shí)，水力—天然裂縫交匯點(diǎn)兩側(cè)裂縫所受正應(yīng)力之差決定了兩側(cè)裂縫的開啟方式。頁(yè)巖儲(chǔ)層水力裂縫與天然裂縫交錯(cuò)機(jī)制主要為兩類：被天然裂縫誘導(dǎo)并沿天然裂縫擴(kuò)展、直接穿透天然裂縫。不同的影響因素對(duì)交錯(cuò)機(jī)制轉(zhuǎn)變的存在一定的閾值，在本文算例中分別對(duì)應(yīng)是：當(dāng)排量超過26 m3/min時(shí)，或壓裂液黏度超過40 mPa·s時(shí)，或應(yīng)力差超過18 MPa時(shí)，或裂縫面摩擦系數(shù)超過0.8時(shí)，水力裂縫將直接穿透天然裂縫。

2）天然裂縫對(duì)水力裂縫的誘導(dǎo)作用隨著兩者夾角的變化而變化，在夾角較小的情況下，水力裂縫容易在最大水平主應(yīng)力方向交匯點(diǎn)之前提前被天然裂縫誘導(dǎo)；而當(dāng)夾角增大時(shí)則在交匯點(diǎn)處進(jìn)入天然裂縫擴(kuò)展，且表現(xiàn)為一側(cè)拉張裂縫一側(cè)剪切裂縫；當(dāng)天然—水力裂縫相互垂直時(shí)，水力裂縫擴(kuò)展至天然裂縫處，容易使天然裂縫兩側(cè)均發(fā)生剪切，并限制水力裂縫在天然裂縫內(nèi)的擴(kuò)展。

3）對(duì)于裂縫發(fā)育頁(yè)巖儲(chǔ)層的多簇水力裂縫擴(kuò)展，相鄰兩簇水力裂縫容易在天然裂縫的誘導(dǎo)作用下相互吸引，進(jìn)而交匯成一條優(yōu)勢(shì)水力主裂縫并向前擴(kuò)展。因此，天然裂縫的分布狀態(tài)是決定水力壓裂裂縫復(fù)雜度的關(guān)鍵性因素。

4）本文的研究結(jié)論可為中深層頁(yè)巖氣壓裂裂縫分析和施工參數(shù)優(yōu)化提供理論與技術(shù)支撐。