卞曉冰，蔣廷學，蘇 瑗，王海濤，劉致嶼

(1.中石化石油工程技術研究院，北京 100101;2.中國石油大學，北京 102249)

引 言

頁巖氣藏具有儲層巖性復雜、基質(zhì)滲透性差、頁理層理發(fā)育等特征，壓裂后形成的網(wǎng)絡裂縫體系是其主要的滲流通道[1-2]。頁巖氣水平井單井壓裂液用量一般為2×104～4×104m3，大量液體在壓裂施工結(jié)束后滯留于裂縫及地層中[3-5]。因此，建立頁巖氣水平井壓裂后排采模型的關鍵是對裂縫形態(tài)的描述和壓裂后地層壓力、飽和度的初始化。目前國內(nèi)外對于頁巖氣水平井壓裂后排采尚沒有進行系統(tǒng)的研究[6-8]，工藝方面主要集中在利用壓裂液和支撐劑示蹤分析壓裂液的返排及井間干擾[9-11];機理方面主要通過數(shù)值模擬手段建立三維氣水兩相模型，考慮了天然裂縫體系、毛管力等單因素的影響[12-13]。本文以川東南一口頁巖氣評價井為例，在模擬壓裂注入過程的基礎上，通過精細描述裂縫體系形態(tài)，研究其對壓裂后排采的影響及主控因素，并以此為基礎進行生產(chǎn)歷史擬合及排采預測。

1 模型的建立

為了盡可能準確地模擬頁巖氣水平井多段壓裂施工后的排采過程，利用Eclipse從頁巖地質(zhì)建立、裂縫形態(tài)描述、壓力場及飽和度場初始化3方面進行精細建模:①采用雙重介質(zhì)模型給頁巖的基質(zhì)系統(tǒng)和天然裂縫系統(tǒng)分別賦值，包括孔隙度、滲透率、含水飽和度和相對滲透率等，應用Langmuir曲線表征頁巖氣的解吸附過程，同時考慮氣體從基質(zhì)向微裂縫的擴散;②基于主、次裂縫是否連通，初步把人工裂縫分為單一縫、復雜縫、網(wǎng)絡縫3種類型，縫高剖面在縱向上逐漸變化;③模擬壓裂井注入過程，井口壓力為施工壓力，注入液量為施工總液量，以此來初始化壓裂后的壓力場和飽和度場，盡可能真實模擬壓裂施工后的狀態(tài)。

以川東南某頁巖氣評價井為例，模型所采用基礎數(shù)據(jù)如下:氣藏大小為2 000 m×1 500 m×60 m，網(wǎng)格數(shù)量為218×75×10個，水平段長為1 260 m，地層壓力為30 MPa，總含氣量為3.91 m3/t(其中吸附氣含量占30%)，基質(zhì)滲透率為10-7μm2，孔隙度為4.6%，共壓裂22段，每段壓裂2簇，簇間距為20 m，每段注入量為1 900 m3，單一裂縫半長為300 m，裂縫導流能力為2 μm2·cm，井底流壓為22 MPa。

2 裂縫參數(shù)影響因素分析

建立單段壓裂模型，以研究不發(fā)生縫間干擾情況下不同裂縫參數(shù)，包括裂縫半長、縫高剖面和裂縫形態(tài)對頁巖氣水平井壓后累計產(chǎn)氣量與返排率的影響。

2.1 裂縫半長的影響

模擬裂縫半長為100～400 m，600 d時的累計產(chǎn)氣量與返排率曲線如圖1所示。由模擬結(jié)果可知:200 m裂縫半長是累計產(chǎn)量和返排率變化的拐點;累計產(chǎn)氣量隨著裂縫半長的增加而增加，但其增幅逐漸減緩;而隨著裂縫半長的增加，返排率逐漸遞減。

2.2 縫高剖面的影響

圖1 600d時不同裂縫半長對應的累計產(chǎn)氣量與返排率

頁巖氣水平井壓裂后，真實縫高剖面在縱向上是變化的，模擬了4種縫高剖面(圖2)，有效裂縫面積依次減小。

圖2 不同縫高剖面模擬圖

表1 600 d時不同縫高剖面對應的累計產(chǎn)氣量與返排率

4種縫高剖面對應600 d時的累計產(chǎn)氣量與返排率曲線如表1所示。由模擬結(jié)果可知:縫高剖面越小累計產(chǎn)氣量越低。但縫高剖面對返排率的影響更為復雜:由于壓裂液基于裂縫面向地層中濾失，因此在相同的總壓裂液量條件下，縫高剖面越小，滯留在近井筒地層中的壓裂液也越多，而近井帶的高壓聚集有利于壓后返排(返排率大小順序為剖面3、剖面2、剖面1);但縫高剖面太小(剖面4裂縫面積僅為剖面1的2/5)會導致壓裂液向垂直于裂縫面更遠處的地層中過度濾失，從而導致返排率下降。

2.3 裂縫形態(tài)的影響

模擬裂縫形態(tài)為單一縫、復雜縫和網(wǎng)絡縫對生產(chǎn)效果的影響(圖3)。其中紅色區(qū)域代表主裂縫，綠色區(qū)域代表次裂縫，次裂縫導流能力設置為主裂縫的1/5。在復雜裂縫中，次裂縫連通了2條主裂縫;而網(wǎng)絡縫中各條主、次裂縫均是連通的。

圖3 不同裂縫形態(tài)模擬圖

3種裂縫形態(tài)對應600 d時的累計產(chǎn)氣量與返排率曲線如表2所示。由模擬結(jié)果可知:裂縫形態(tài)越復雜，累計產(chǎn)氣量越大，但返排率較低。以網(wǎng)絡縫為例，由于次裂縫的導流能力較低，大量存在于其中的壓裂液是網(wǎng)絡縫返排率較低的原因之一。因此，相互連通的網(wǎng)絡縫可極大提高產(chǎn)氣量，但增加了壓裂液的返排難度。

表2 600 d時不同裂縫形態(tài)對應的累計產(chǎn)氣量與返排率

3 裂縫參數(shù)主控因素分析

在考慮段間干擾和縫間干擾的情況下，為了全面并定性地分析各裂縫參數(shù)對頁巖氣水平井壓后排采的影響程度，對裂縫半長、縫高剖面和裂縫形態(tài)各取3個值，共設計27種方案進行數(shù)值模擬(表3)。

表3 裂縫參數(shù)取值

27個方案600 d的累計產(chǎn)氣量和返排率值如表4所示，方差分析結(jié)果見表5。綜合評價結(jié)果顯示，縫高剖面和裂縫形態(tài)是影響壓后產(chǎn)氣量的顯著因素，而裂縫形態(tài)對返排率影響最大，但裂縫參數(shù)對返排率均無顯著影響。因此，在對于頁巖氣水平井進行壓裂設計時，可以采取變排量、膠液前置造縫、混合加砂等措施，保證縫高在縱向上充分延伸并盡可能形成復雜裂縫體系，以提高壓裂后生產(chǎn)效果。而對于已施工井，頁巖氣水平井壓裂后裂縫參數(shù)的準確判斷則是建立可靠排采模型的基礎。

表4 600d時27個方案對應的累計產(chǎn)氣量與返排率

表5 裂縫參數(shù)對累計產(chǎn)氣量及返排率影響程度顯著性對比

4 實例應用

川東南某頁巖氣水平井投產(chǎn)后自噴產(chǎn)氣、水量少，后下電泵抽汲。根據(jù)施工曲線擬合結(jié)果及G函數(shù)曲線特征，診斷裂縫半長為250～300 m;縫高剖面為楔形。應用本文裂縫參數(shù)模擬方法建立了該井壓裂后排采模型，并進行生產(chǎn)歷史擬合(圖4)?？梢?，該模型能反應實際井的生產(chǎn)動態(tài)，同時也驗證了本文所建立排采模型方法的可靠性。

圖4 某頁巖氣井生產(chǎn)歷史擬合結(jié)果

對于下泵階段，可根據(jù)氣相流入流出曲線的生產(chǎn)協(xié)調(diào)點獲得優(yōu)化的井底流壓[14]，最優(yōu)值為16.4 MPa，對應產(chǎn)氣量為2.2×104m3/d。以井底流壓為基礎對擬合后的模型進行排采預測，與真實排采情況對比，壓后3個月的累計產(chǎn)量增加了34.7%，返排率增加了2.5個百分點。此優(yōu)化結(jié)果對該井下一步措施調(diào)整具有一定的指導意義。

5 結(jié)論

(1)建立了頁巖氣水平井多段壓裂排采模型，模型考慮了頁巖地質(zhì)特征、裂縫參數(shù)特征和壓力場及飽和度場的初始化，盡可能準確地模擬了頁巖氣水平井多段壓裂排采過程。

(2)裂縫參數(shù)對頁巖氣水平井壓后排采影響較大，主控因素分析結(jié)果顯示，縫高剖面和裂縫形態(tài)是影響壓后產(chǎn)氣量的顯著因素，但裂縫參數(shù)對返排率無顯著影響。

(3)對川東南某頁巖氣水平井的歷史擬合結(jié)果表明，本文所建立的排采模型能真實反應實際井的生產(chǎn)動態(tài)，初步優(yōu)化了下泵排采的生產(chǎn)制度，優(yōu)化結(jié)果對該井下一步措施調(diào)整具有一定的指導意義。

[1]李慶輝，陳勉，金衍，等.壓裂參數(shù)對水平頁巖氣井經(jīng)濟效果的影響[J].特種油氣藏，2013，20(1):146-150.

[2]單衍勝，張金川，李曉光，等.遼河東部凸起上古生界海陸過渡相頁巖氣富集條件[J].大慶石油地質(zhì)與開發(fā)，2012，31(3):13-18.

[3]陳尚斌，朱炎銘，王紅巖，等.中國頁巖氣研究現(xiàn)狀與發(fā)展趨勢[J]. 石油學報，2010，31(4):689-694.

[4]薛承瑾.頁巖氣壓裂技術現(xiàn)狀及發(fā)展建議[J].石油鉆探技術，2011，39(3):24-29.

[5]孫海成，湯達禎，蔣廷學，等.頁巖氣儲層壓裂改造技術[J]. 油氣地質(zhì)與采收率，2011，18(4):90-93.

[6]錢旭瑞，劉廣忠，唐佳，等.頁巖氣井產(chǎn)能影響因素分析[J]. 特種油氣藏，2012，19(3):81-83.

[7]孫海成，湯達禎，蔣廷學.頁巖氣儲層裂縫系統(tǒng)影響產(chǎn)量的數(shù)值模擬研究[J].石油鉆探技術，2011，39(5):63-67.

[8]畢曉明，邵銳，唐亞會，等.頁巖氣滲流機理對氣藏開采的影響[J].大慶石油地質(zhì)與開發(fā)，2011，30(2):80-83.

[9]Clarkson C R.Modeling 2-phase flowback of multifractured horizontal wells completed in shale[C].SPE162593，2012:28-30.

[10]孫海，姚軍，孫致學，等.頁巖氣數(shù)值模擬技術進展及展望[J]. 油氣地質(zhì)與采收率，2012，19(1):46-49.

[11]William G E，Stuart J E，Thomas E E.New single-well standalone gas lift process facilitates Barnett Shale fracture-treatment flowback[C].SPE118876，2010:4-8.

[12]Cheng Y.Impact of water dynamics in fractures on the performance of hydraulically fractured wells in gasshale reservoirs[C].SPE127863，2012:10-12.

[13]James W C.Flowback performance in intensely naturally fractured shale gas reservoirs[C].SPE131785，2010:23-25.

[14]蔣廷學，卞曉冰，王海濤，等.頁巖氣水平井分段壓裂排采規(guī)律研究[J].石油鉆探技術，2013，41(5):21-25.