摘 要

　　　　　　——以川西頁巖氣藏某水平井為例摘　要：頁巖儲層需要水力壓裂才能獲得理想的產(chǎn)能，壓裂時在追求較大改造體積的同時也應(yīng)注重形成與儲層相匹配的縫網(wǎng)導流

　　　　　　——以川西頁巖氣藏某水平井為例

摘　要：頁巖儲層需要水力壓裂才能獲得理想的產(chǎn)能，壓裂時在追求較大改造體積的同時也應(yīng)注重形成與儲層相匹配的縫網(wǎng)導流能力，以提高改造體積滲透率。基于Warren Root模型，將壓裂后形成的縫網(wǎng)考慮為高滲透帶，利用等效滲流理論建立了等效高滲透帶模型，在地質(zhì)模型中利用體積及等效滲透率對壓裂縫網(wǎng)進行表征，通過產(chǎn)能模擬并借助凈現(xiàn)值理論對高滲透帶長度、等效滲透率等參數(shù)進行優(yōu)選，并以優(yōu)選的高滲透帶參數(shù)為目標，結(jié)合縫網(wǎng)模擬便可得到目標條件下的最優(yōu)施工參數(shù)。針對川西頁巖氣藏某水平井使用該方法得到最優(yōu)高滲透帶長度為200～220m，最優(yōu)等效滲透率為4～5mD。結(jié)合縫網(wǎng)模擬得到目標條件下的施工參數(shù)為：總液量為l600m³，總砂量為53m³，平均砂比為l0％，最高砂比為28％，施工排量為10m³／min。該設(shè)計為頁巖氣儲層改造作業(yè)提供了技術(shù)支撐。

關(guān)鍵詞：頁巖  水平井  水力壓裂  縫網(wǎng)  壓裂設(shè)計Warren-Root模型  等效滲流理論  高滲透帶  四川盆地

An optimal design of multi-stage fracking for horizontal shale gas wells：A case study from the western Sichuan Basin

Abstract：Hydraulic fracturing treatment is a major technique stimulating the field of a gas shale reservoir．While a large stimulated volume is pursued，the fracture network conductivity should be controlled so well to match the reservoir that the permeabilitv of the stimulated volume is improved．Based on the Warren-Root model and with the fracture network upon fracturing taken as the highpermeability zone，the equivalent high permeability zone model was established by use of the equivalent seePage theory；secondly，the network was characterized in terms of volume and equivalent permeability in the geological model；thirdly，through performing capacity simulation while utilizing the net present value theory，such parameters were optimized as the length and equivalent permeability of the high permeability zone；and on this basis，optimal construction parameters wete obtained under the objective condition combining the network simulation．In a case study from the western Sichuan Basin，the optimized high permeabilitv zone length ranged from 200 to 220m，and the optimized equivalent permeability from 4 to 5mD．Combined with the network simulation，the construction parameters under the objective condition are：the total fluid volume of l600m³，the total sand volume of 53m³，the average sand rate of l0％，the maximum sand rate of 28％，and construction displacement of l0m³／min．This design provides a technical support for stimulating shale gas reservoirs．

Keywords：shale gas，horizontal well，multi-stage fracking，network， fracturing design，Warren Root model， equivalent seepage theory，high permeability zone，Sichuan Basin

頁巖油氣資源對于緩解國內(nèi)能源供需矛盾、保障國家能源安全具有重要作用^[1-2]。頁巖儲層具有極低的滲透率和孔隙度，其最主要的儲集和運移通道來自一系列的天然裂縫網(wǎng)絡(luò)，因此需要通過水力壓裂打碎有效儲集體，溝通裂縫網(wǎng)絡(luò)，創(chuàng)造油氣運移通道，增大裂縫面與基質(zhì)接觸面積，從而實現(xiàn)基質(zhì)中油氣從任意方向向裂縫最短距離滲流。一直以來體積壓裂以增加改造體積為主要目標，忽略了其與縫網(wǎng)導流能力的匹配。Warpinski、Ge等^[4-5]研究認為，頁巖壓裂在實現(xiàn)較大改造體積的同時，還需獲得一定的改造帶內(nèi)縫網(wǎng)導流能力，以提高改造帶內(nèi)整體滲透率。張小濤等^[6]認為，頁巖壓裂縫只有具有一定導流能力時才可能成為對氣藏開發(fā)有效的裂縫。于榮澤等^[7]認為頁巖壓裂誘導縫存在著最優(yōu)導流能力。Mayerhofer等^[8]研究認為縫網(wǎng)導流能力介于0.015～1.5D·cm時才能夠獲得滿意的產(chǎn)能。因此，獲得與縫網(wǎng)相適應(yīng)的裂縫導流能力對于頁巖壓裂同樣重要。頁巖壓裂設(shè)計不同于常規(guī)雙翼裂縫^[9]，目前主要依賴于產(chǎn)能預(yù)測來優(yōu)化縫網(wǎng)幾何參數(shù)以及縫網(wǎng)導流能力，產(chǎn)能預(yù)測主要有數(shù)學解析模型和使用離散網(wǎng)絡(luò)模型(DFN)和雙重孔隙介質(zhì)模型進行數(shù)值模擬等方法^[10]，解析模型準確性差，數(shù)值模擬過程復雜，參數(shù)精度要求較高，不便于現(xiàn)場應(yīng)用。因此，綜合考慮改造體積以及縫網(wǎng)導流能力與儲層的匹配，建立一種簡單實用的設(shè)計方法對于指導頁巖壓裂設(shè)計具有重要意義。

基于Warren-Root模型將壓裂后形成的縫網(wǎng)考慮為高滲透帶，利用等效滲流理論建立了等效高滲透帶模型，在地質(zhì)模型中利用體積及等效滲透率對壓裂縫網(wǎng)進行表征，通過產(chǎn)能模擬并借助凈現(xiàn)值理論對高滲透帶長度、等效滲透率等參數(shù)進行優(yōu)選，并以優(yōu)選的高滲透帶參數(shù)為目標，結(jié)合縫網(wǎng)模擬得到目標條件下的最優(yōu)施工參數(shù)，建立了頁巖水平壓裂優(yōu)化的新方法。針對川西頁巖氣藏某水平井應(yīng)用該方法進行了優(yōu)化沒計，取得了良好的效果。

1　川西頁巖氣儲層特征

四川盆地上三疊統(tǒng)須家河組五段是川西坳陷最主要的烴源巖發(fā)育層段，地層孔隙度為3％左右，滲透率為0.01mD左右，為典型的低孔隙度、低滲透率頁巖氣藏，需進行縫網(wǎng)壓裂才能實現(xiàn)有效增產(chǎn)。該層段脆性指數(shù)為54，水平段應(yīng)力差異系數(shù)為0.24，且發(fā)育大量天然裂縫，具備壓裂形成復雜裂縫的條件^[11-12]。前期微地震監(jiān)測表明該層段壓裂形成了復雜裂縫，結(jié)果如圖l所示。

 

2　縫網(wǎng)壓裂優(yōu)化設(shè)計步驟

針對川西頁巖壓裂形成復雜裂縫的情況建立了優(yōu)化縫網(wǎng)參數(shù)及施工參數(shù)的設(shè)計步驟：

1)利用等效滲流原理將縫網(wǎng)帶等效為高滲透帶，利用體積及等效滲透率在地質(zhì)模型中對縫網(wǎng)進行表征，建立高滲透帶參數(shù)優(yōu)化模型。

2)通過產(chǎn)能模擬并借助凈現(xiàn)值理論優(yōu)化高滲透帶長度、等效滲透率等參數(shù)。

3)以優(yōu)選的高滲透帶參數(shù)為目標，利用壓裂設(shè)計軟件模擬縫網(wǎng)參數(shù)，通過參數(shù)反演并與最優(yōu)參數(shù)對比，達到日標時的施工參數(shù)即為優(yōu)化得到的施工參數(shù)。

3　壓裂優(yōu)化模型

3．1　等效高滲透帶模型的建立

頁巖在壓裂后以射孔段為中心形成復雜裂縫網(wǎng)絡(luò)，每一簇復雜裂縫可視為沿井筒方向分布的高滲透帶，根據(jù)等效滲流理論，建立了將縫網(wǎng)帶等效為高滲透帶的等效高滲透帶模型。

模型假設(shè)條件：①壓裂形成的裂縫網(wǎng)絡(luò)使用warren-Root模型進行表征；②支撐劑均勻鋪置在支撐裂縫網(wǎng)絡(luò)中，縫網(wǎng)體積不等于支撐劑體積；③高滲透帶中向井筒滲流等效為高滲透帶中基質(zhì)向井筒滲流和縫網(wǎng)向井筒滲流兩部分；④高滲透帶的滲流為穩(wěn)定滲流。

高滲透帶滲流單元由基質(zhì)及其附近的裂縫空間組成如圖2所示。其中基質(zhì)內(nèi)的滲流可表示為：

 

 

式中qml為基質(zhì)中流體流量，m³；K。為基質(zhì)滲透率，mD；Am為基質(zhì)截面積，m²；p1為入口端壓力，MPa；p2為出口端壓力，MPa；m為流體黏度，mPa·s；Lm為單元體內(nèi)基質(zhì)體長度，m。

裂縫內(nèi)的滲流可表示為：

 

式中qf1為單元體內(nèi)裂縫流量，m³；Kf為裂縫滲透率，mD；Af為基質(zhì)截面積，m²；wf為裂縫寬度，m。

因此對于整個裂縫網(wǎng)絡(luò)，有

 

式中`K為等效滲透率，mD；qt為總流量，m³。

根據(jù)等效滲流理論，有

qt＝qm+qf 　    　　　         　　　 (6)

結(jié)合式(3)～(6)得到： 

 

 

對于單元體而言Lm≈wf，因此式(7)、式(8)可寫成：

 

式(9)為高滲透帶等效滲透率的計算表達式，其為縫網(wǎng)體積、縫網(wǎng)滲透率、基質(zhì)體積、基質(zhì)滲透率的加權(quán)平均。通過此關(guān)系式將改造體積與縫網(wǎng)導流能力的優(yōu)化轉(zhuǎn)變?yōu)楦邼B透帶體積與等效滲透率的優(yōu)化。

3．2　地質(zhì)模型的建立

考慮吸附氣解析作用，根據(jù)地質(zhì)特征建立了的壓裂水平井地質(zhì)模型(圖3)，模型長l200 ，寬800m，厚80m，基質(zhì)滲透率為0.01mD，水平段長為767m，單井控制面積為96×l0⁴m²。以川西頁巖氣藏某水平井為例(圖2)，將9段裂縫網(wǎng)絡(luò)等效為9段高滲透帶植入地質(zhì)模型。

 

4　高滲透帶參數(shù)優(yōu)化

4．1　縫網(wǎng)長度優(yōu)化

結(jié)合經(jīng)濟評價模型，通過模擬高滲透帶長度為80～260m下的產(chǎn)量，利用凈現(xiàn)值法優(yōu)選裂縫長度，凈現(xiàn)值表達式為：

 

式中NPV為凈現(xiàn)值，元；RF為壓裂井現(xiàn)值，元；RO為未壓裂井現(xiàn)值，元；CF為壓裂施工成本，元；n為年限；(VF)j為壓裂井第j年度總收入，元；i為貼現(xiàn)率；(VO)j為未壓裂井第j年度總收入，元。

通過模擬不同縫網(wǎng)長度下的壓裂井產(chǎn)能，并分析凈現(xiàn)值與縫網(wǎng)長度的關(guān)系，其模擬結(jié)果如圖4、5所示。隨著裂縫長度的增加，累計產(chǎn)量與凈現(xiàn)值在200～220m均出現(xiàn)拐點。因此優(yōu)選的縫網(wǎng)長度為200～220m，此長度下地質(zhì)模型中的改造體積為128×10⁴m³。

 

 

4．2　等效滲透率優(yōu)化

改變地質(zhì)模型中高滲透帶滲透率，模擬高滲透帶滲透率分別為2.5、3.0、3.5、4.0、5.0、6.0、7.0mD時的產(chǎn)量，分析凈現(xiàn)值與高滲透帶等效滲透率的關(guān)系，其模擬結(jié)果見圖6、7。隨著高滲透帶滲透率的增加產(chǎn)量不斷增加，當高滲透帶等效滲透率增加到4mD后，產(chǎn)量增幅及凈現(xiàn)值曲線均趨于平穩(wěn)，因此優(yōu)選的高滲透帶滲透率為4～5mD。

 

 

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

以優(yōu)化得到的高滲透帶長度、滲透率為目標，借助頁巖壓裂設(shè)計軟件對不同液量、不同砂量下的壓裂網(wǎng)絡(luò)裂縫參數(shù)進行模擬。模擬發(fā)現(xiàn)：改造體積與液量成正相關(guān)關(guān)系，縫網(wǎng)長度與前置液成正相關(guān)關(guān)系，它們均與砂量無關(guān)，給定液量下裂縫導流能力與加砂量成正相關(guān)關(guān)系。最終在總液量為l600m³、總砂量為53m³下得到縫長202m，改造體積102×10⁴m³ (其小于地質(zhì)模型中改造體積的原因在于模擬得到的裂縫網(wǎng)絡(luò)為橢球形，而地質(zhì)模型中的等效高滲透帶為長方體，后續(xù)研究中將對地質(zhì)模型進一步改進)，縫網(wǎng)平均導流能力為0.039D·cm，平均縫網(wǎng)寬度為0.37mm，結(jié)合縫網(wǎng)參數(shù)計算得到等效滲透率為4.49mD，滿足設(shè)計要求。模擬得到的網(wǎng)絡(luò)裂縫參數(shù)如下：裂縫網(wǎng)絡(luò)長404m，裂縫網(wǎng)絡(luò)寬40m，裂縫網(wǎng)絡(luò)高92m，平均DFN寬度為0.37mm，改造體積(SRV)1025500m³，DFN體積為43lm³，DFN平均導流能力為0.0397D·cm，等效滲透率為4.49mD。網(wǎng)絡(luò)裂縫形態(tài)如圖8、9所示。

 

 

根據(jù)模擬結(jié)果推薦的施工參數(shù)為：總液量l600m³，總砂量53m³，平均砂比10％，最高砂比28％，施工排量10m³／min。

6　結(jié)論

1)基于Warren Root模型將壓裂后形成的縫網(wǎng)考慮為高滲透帶，利用等效滲流理論建立了頁巖水平井縫網(wǎng)壓裂優(yōu)化設(shè)計的新方法。

2)針對川西頁巖氣藏某水平井使用該方法得到高滲透帶長度為200～220m，最優(yōu)等效滲透率為4～5mD。施工參數(shù)為：總液量l600m³，總砂量53m³，平均砂比10％，最高砂比28％，施工排量為10m³／min。

3)建立的方法能夠有效解決復雜縫網(wǎng)參數(shù)優(yōu)化問題，對類似儲層的縫網(wǎng)壓裂具有一定的指導意義。

 

參考文獻

[1]胡文瑞，鮑敬偉．探索中國式的頁巖氣發(fā)展之路[J]．天然氣工業(yè)，2013，33(1)：1-7．

HU Wenrui，BAO Jingwei．To explore the way of Chinesestyle shale gas development[J]．Natural Gas Industry，2013，33(1)：1-7．

[2]張金川，林臘梅，李玉喜，等．頁巖油分類與評價[J]．地學前緣，2012，1 9(5)：322-331．

ZHANG Jinchuan，LIN Lamei，LI Yuxi，et al．Classificarion and evaluation of shale oil[J]．Earth Science Frontiers，2012，19(5)：322-331．

[3]吳奇，胥云，王騰飛，等．增產(chǎn)改造理念的重大變革——體積改造技術(shù)概論[J]．天然氣工業(yè)，2011，31(4)：7-12．

WU Qi，XU Yun，WANG Tengfei，et al．The revolution of reservoir stimulation：An introduction of volume fracturing[J]．Natural Gas Industry，2011，3l(4)：7-12．

[4] WARPINSKI N R，MAYERHOFER M J，VINCENT M C，et al．Stimulating unconventional reservoirs：Maximizing network growth while optimizing fracture conductivity[C]//paper ll4173 presented at the SPE Unconventional Res ervoirs Conference，10-12 February 2008，Keystone，Colorado，USA．New York：SPE，2008．

[5]GE J，GHASSEMI A．Permeability enhancement in shale gas reservoirs after stimulation by hydraulic fracturing[C]//paper ll-514 presented at the 45^th U．S．Rock Mechanics／Geomechanics Symposium，26-29 June 2011，San Francis co，California，USA．New York：ARMA，2011．

[6]張小濤，吳建發(fā)，馮曦，等．頁巖氣藏水平井分段壓裂滲流特征數(shù)值模擬[J]．天然氣工業(yè)，2013，33(3)：47-52．

ZHANG Xiaotao，WU Jianfa，FENG Xi，et al．Numerical simulation of seepage flow characteristics of multi stage fracturing(MSF)in horizontal shale gas wells[J]．Natural Gas Industry，2013，33(3)：47-52．

[7]于榮澤，張曉偉，卞亞南，等．頁巖氣藏流動機理與產(chǎn)能影響因素分析[J]．天然氣工業(yè)，2012，32(9)：10-15．

YU Rongze，ZHANG Xiaowei，BIAN Yanan，et al．Flowmechanism of shale gas reservoirs and influential factors of their productivity[J]．Natural Gas Industry，2012，32(9)：10-15．

[8]MAYERFIOFER M J，LOLON E P，YOUNGBLOOD I E，et al．Integration of microseismic fracture mapping resuits with numerical fracture network production modeling in the Barnett Shale[C]//paper l02-103 presented at the SPE Annual Technical Conference and Exhibition,24-27 September 2006，San Antonio，Tcxas，USA．New York：SPE，2006．

[9]郭建春，梁豪，趙志紅．基于最優(yōu)支撐劑指數(shù)法優(yōu)化低滲氣藏裂縫參數(shù)[J]．西南石油大學學報：自然科學版，2013，35(1)：93-98．

GUO Jianchun，LIANG Hao，ZHA()Zhihong．Optimizing the fracture parameters of low permeability gas reservoirs[J]．Journal of Southwest Petroleum University：Science& Technology Edition，2013，35(1)：93-98．

[10]程遠方，董丙響，時賢，等．頁巖氣藏三孔雙滲模型的滲流機理[J]．天然氣工業(yè)，2012，32(9)：44-47．

CFIENG Yuanfang，DONG Bingxiang，SHI Xian，et al．Seepage mechanism of a triple-porosity／dual-permeability model for shale gas reservoirs EJ]．Natural Gas Industry，2012，32(9)：44-47．

[11]趙金洲，任嵐，胡永全．頁巖儲層壓裂縫成網(wǎng)延伸的受控因素分析[J]．西南石油大學學報：自然科學版，2013，35(1)：1-9．

ZHAO Jinzhou，REN Lan，HU Yongquan．Controlling factors of hydraulic fractures extending into network in shale formations[J]．Journal of Southwest Petroleum University：Science&Technology Edition，2013，35(1)：1-9．

[12]RICKMAN R，MULLEN M，PETRE E，et al．A Dractical use of shale petrophysics for stimulation design optimization：All shale plays are not clones of the Barnett Shale[C]//paper ll5-258 presented at the SPE Annual Technical Conference and Exhibition，21-24 September 2008，Denver，Colorad0，USA．New York：SPE，2008．

 

本文作者：郭建春  梁豪  趙志紅  王興文  林立世

作者單位：“油氣藏地質(zhì)及開發(fā)工程”國家重點實驗室·西南石油大學

  中國石化集團兩南油氣分公司工程技術(shù)研究院