鄭松青

（中國石化石油勘探開發(fā)研究院，北京 100083）

非均勻分段壓裂水平井產(chǎn)能計算

鄭松青

（中國石化石油勘探開發(fā)研究院，北京 100083）

為考慮非均勻分布裂縫對水平井產(chǎn)能的影響，在保角變換的基礎(chǔ)上，引入局部坐標(biāo)系，建立非均勻分段壓裂水平井的產(chǎn)能計算模型，考慮裂縫不同長度、不同間距、同水平井呈一定夾角及裂縫間的相互干擾對產(chǎn)能的影響，對紅河地區(qū)3口分段壓裂水平井產(chǎn)能進(jìn)行計算，并同郎兆新等建立的模型進(jìn)行對比，結(jié)果顯示模型具有更小的誤差，但對輸入?yún)?shù)要求更高．

保角變換；局部坐標(biāo)系；多段壓裂；水平井產(chǎn)能

0 引言

水平井分段壓裂是低滲致密油氣藏、頁巖油氣藏、稠油油藏等復(fù)雜油氣藏高效開發(fā)的主要技術(shù)手段之一［1—3］．產(chǎn)能計算是水平井分段壓裂優(yōu)化設(shè)計的基礎(chǔ)．壓裂水平井產(chǎn)能計算方法有兩類，一是數(shù)值模擬方法［4—6］，二是解析、半解析方法［7—19］．在數(shù)值模擬方法中，早期一般對裂縫做等效處理［5］，隨著PEBI網(wǎng)格的引入，在一定程度上提高裂縫描述和模擬的精度［6］．解析、半解析方法也可分兩類，一是非穩(wěn)態(tài)模型［7—9］，二是穩(wěn)態(tài)模型［10—19］．同非穩(wěn)態(tài)模型相比，穩(wěn)態(tài)模型計算簡單，應(yīng)用方便，因而得到更多關(guān)注．Norris S O［10］、Rajagopal R［11］、Hegre T M［12］、Guo Boyun［13—14］等通過不同的方式建立計算模型．郎兆新［15］、范子菲［16］、丁一萍［17—18］等也建立一系列模型．其中郎兆新應(yīng)用位勢理論和疊加原理建立的計算模型得到廣泛應(yīng)用．在該模型中，假設(shè)裂縫為垂直裂縫、等間距分布、高度等于儲層的厚度，且產(chǎn)量相等．寧正福等［19］在此基礎(chǔ)上進(jìn)行改進(jìn)，考慮裂縫間的相互干擾，但仍然假設(shè)裂縫同井筒垂直、等間距分布．受儲層條件和工藝水平的影響，這些假設(shè)過于理想化．

筆者在保角變換的基礎(chǔ)上，引入局部坐標(biāo)系，建立考慮裂縫非均勻分布的水平井產(chǎn)能計算模型．

1 單條裂縫在平面上的勢分布

單條裂縫的勢分布函數(shù)是分段壓裂水平井產(chǎn)能計算的基礎(chǔ)．對于有限長度裂縫的勢分布函數(shù)，郎兆新等并未給出推導(dǎo)過程［15］，程林松利用保角變換研究勢的分布函數(shù)，但在后續(xù)推導(dǎo)中，將橢圓等勢線簡化為圓形等勢線［20］．文中采用保角變換，重新對其進(jìn)行推導(dǎo)．假設(shè)一條裂縫長度為2L，壓開整個儲層，忽略裂縫內(nèi)部的壓力損失，將三維流動簡化為二維流動．裂縫在Z平面坐標(biāo)系中的位置見圖1．

采用保角變換的方法研究其勢分布，取變換函數(shù)，即

式中：z＝x+i y，Z為平面上一點(diǎn)；x＝u+i v；W為平面上一點(diǎn)．

通過保角變換，將Z平面上半平面變換為W 平面上帶寬為π的右邊無限大平面，將Z平面下半平面變換為W 平面上帶寬為π的左邊無限大平面，將長度為2L的裂縫變換為寬度為π的排液坑道．于是，Z平面上向裂縫的流動經(jīng)過變換變成W 平面上的坑道流（見圖2，Φ為勢分布）．坑道流勢分布簡單，可據(jù)此解決裂縫流動的問題．

圖1 裂縫分布示意Fig．1 Diagram of a fracture in Z plane

圖2 裂縫流動保角變換示意Fig．2 Diagram of conformal transformation

根據(jù)雙曲函數(shù)的性質(zhì)，對式（1）實(shí)部和虛部分離后，可得

根據(jù)式（2），可得u同x，y的關(guān)系為

如果裂縫的產(chǎn)量為q，儲層厚度為h，對于W 平面上的坑道流，勢的分布函數(shù)Ф為

將式（3）代入式（4），可得單條裂縫在Z平面上勢的分布為

式中：（x，y），（u，v）分別為Z平面和W平面上點(diǎn)的坐標(biāo)；C為常數(shù)．

式（5）勢函數(shù)假設(shè)裂縫中點(diǎn)在坐標(biāo)原點(diǎn)，裂縫位于x軸．如果裂縫的中點(diǎn)在（x0，y0），且同y軸夾角為θ，就不能直接使用式（5）．為解決這個問題，引入局部坐標(biāo)系（見圖3），在局部坐標(biāo)系下變量為（x′，y′）滿足式，根據(jù)全局坐標(biāo)系同局部坐標(biāo)系的關(guān)系，可得平面上任一條裂縫的勢分布為

圖3 任意位置裂縫局部坐標(biāo)系與全局坐標(biāo)系關(guān)系示意Fig．3 Diagram of local and global coordinate system

2 產(chǎn)能計算模型

一口水平井壓裂多條裂縫，不考慮未壓裂區(qū)域水平井產(chǎn)能，多條裂縫產(chǎn)能可通過勢的疊加原理求得．根據(jù)勢的疊加原理，如果存在N條裂縫，空間某點(diǎn)（x，y）的勢為

式中：x′i＝（x—xi0）sinθi—（y—yi0）cosθi；y′i＝（x—xi0）cosθi—（y—yi0）sinθi；qi為第i條裂縫的產(chǎn)量；（xi0，yi0）為第i條裂縫的中點(diǎn)坐標(biāo)；θi為第i條裂縫同y軸（水平井）的夾角；Li為第i條裂縫的半長．

令

將式（7）簡化為

假設(shè)第j條裂縫同水平井井筒交點(diǎn)坐標(biāo)為（0，yj），該處的壓力為井底流壓pwf，在邊界（xe，ye）處壓力為pe．根據(jù)勢同壓力的關(guān)系，有

將式（9）代入式（10），可得

式中：μ為原油黏度；K為儲層滲透率；pj為第j條裂縫的壓力（井底流壓）．

則水平井的產(chǎn)能Q為

3 計算實(shí)例

分別利用文中建立的非均勻模型和郎兆新等建立的模型，計算紅河油田37井區(qū)3口分段壓裂水平井產(chǎn)能，并與實(shí)際產(chǎn)量進(jìn)行對比（見表1）．其中HH37P24井的基本參數(shù)見表2～3．

表1 2種模型計算結(jié)果與實(shí)際結(jié)果對比Table 1 2Results comparison between calculation and real data t·d—1

表2 HH37P24井基本參數(shù)Table 2 Well parameters of HH37P24well

表3 HH37P24井壓裂參數(shù)Table 3 Fractures parameters of HH37P24well

由表3可知，受到地應(yīng)力及施工工藝、地層物性等影響，壓裂裂縫同水平井完全垂直，長度一致是很困難的．文中非均勻模型考慮裂縫同水平井夾角、裂縫長度不一致等因素影響，得出的結(jié)果同實(shí)際相比存在較大誤差，但比等效模型結(jié)果誤差要?。饕蚴橇芽p間存在干擾，裂縫長度、相對位置、間距等不同，干擾程度不同，整口井的產(chǎn)能也會出現(xiàn)差異．

同郎兆新等建立的模型相比，文中模型雖然誤差較小，但對輸入?yún)?shù)要求更高，這些參數(shù)在現(xiàn)場獲取相對困難．

4 結(jié)論

（1）在保角變換的基礎(chǔ)上，引入局部坐標(biāo)系，推導(dǎo)非均勻分段壓裂水平井產(chǎn)能公式．該公式考慮裂縫不等長、非均勻分布及同水平井不垂直等因素的影響，更符合實(shí)際情況．

（2）同郎兆新等建立的模型進(jìn)行對比，文中模型具有更小誤差，但對參數(shù)要求更高．

［1］ 劉振宇，劉洋，賀麗艷，等．人工壓裂水平井研究綜述［J］．大慶石油學(xué)院學(xué)報，2002，26（4）：96—99．Liu Zhenyu，Liu Yang，He Liyan，et al．Overview of the research of hydraulically fractured horizontal wells［J］．Journal of Daqing Petroleum Institute，2002，26（4）：96—99．

［2］ 王樹平，史云清，嚴(yán)謹(jǐn)，等．考慮裂縫干擾的氣藏壓裂水平井產(chǎn)能預(yù)測模型［J］．大慶石油學(xué)院學(xué)報，2012，36（3）：67—71．Wang Shuping，Shi Yunqing，Yan Jin，et al．A new productivity prediction model of fracturing horizontal well in gas reservoir taking into account the mutual interference between the fracture［J］．Journal of Daqing Petroleum Institute，2012，36（3）：67—71．

［3］ 任俊杰，郭平，王德龍，等．頁巖氣藏壓裂水平井產(chǎn)能模型及影響因素［J］．東北石油大學(xué)學(xué)報，2012，36（6）：76—81．Ren Junjie，Guo Ping，Wang Delong，et al．Productivity model and influence factors for fracturing horizontal wells in shale gas［J］．Journal of Daqing Petroleum Institute，2012，36（6）：76—81．

［4］ Machiel B，Loyd E，Igor Kuvshinov．The potential of multiple fractured horizontal wells in layered reservoirs［R］．SPE 102633，2006．

［5］ 程林松，皮建，廉培慶，等．裂縫性油藏水平井產(chǎn)能計算方法［J］．計算物理，2011，28（2）：230—236．Cheng Linsong，Pi Jian，Lian Peiqing，et al．A computational method for productivity of horizontal well in naturally fractured reser—voirs［J］．Chinese Journal of Computational Physics，2011，28（2）：230—236．

［6］ 吳軍來，劉月田．基于PEBI網(wǎng)格加密的水平井分段壓裂模擬［J］．大慶石油學(xué)院學(xué)報，2010，34（6）：53—57，63．Wu Junlai，Liu Yuetian．Simulation of staged fracturing for horizontal wells based on PEBI grid refinements［J］．Journal of Daqing Petroleum Institute，2010，34（6）：53—57，63．

［7］ 徐嚴(yán)波，齊桃，楊鳳波，等．壓裂后水平井產(chǎn)能預(yù)測新模型［J］．石油學(xué)報，2006，27（1）：89—91，96．Xu Yanbo，Qi Tao，Yang Fengbo，et al．New model for productivity test of horizontal well after hydraulic fracturing［J］．Acta Petro—lei Sinica，2008，15（2）：64—68．

［8］ 閃從新，李曉平，秦海菲，等．多裂縫水平井非穩(wěn)態(tài)產(chǎn)能模型及計算方法研究［J］．油氣井測試，2008，17（6）：5—7．Shan Congxin，Li Xiaoping，Qin Haifei，et al．Study on the unsteady productivity model of horizontal well with multiple transversal fracture and calculation method［J］．Well Testing，2008，17（6）：5—7．

［9］ 廉培慶，程林松，曹仁義，等．低滲透油藏壓裂水平井井筒與油藏耦合的非穩(wěn)態(tài)模型［J］．計算物理，2010，27（2）：203—210．Lian Peiqing，Cheng Linsong，Cao Renyi，et al．A coupling model of low permeability reservoir and fractured horizontal wellbore in nonsteady state［J］．Chinese Journal of Computational Physics，2010，27（2）：203—210．

［10］ Norris S O，Hunt J L，Soliman M Y，et al．Predicting horizontal well performance：A review of current technology［R］．SPE 21793，1991．

［11］ Rajagopal R，Sada D，Joshi．Productivity of multiple drainholes or fractured horizontal wells［R］．SPE 21243，1993．

［12］ Hegre T M，Larsen L．Productivity of multi—fractured horizontal wells［R］．SPE 28845，1994．

［13］ Guo Boyun，Schechter D S．A simple and acturate mathematical model for predicting productivity of multifractured horizontal wells ［R］．SPE 38104，1997．

［14］ Guo Boyun，Yu Xiance．A Simple and accurate mathematical model for predicting productivity of multi—fractured horizontal wells ［R］．SPE 114452，2008．

［15］ 郎兆新，張麗華，程林松．壓裂水平井產(chǎn)能研究［J］．石油大學(xué)學(xué)報：自然科學(xué)版，1994，18（2）：43—46．Lang Zhaoxin，Zhang Lihua，Cheng Linsong．Investigation on productivity of fractured horizontal well［J］．Journal of the University of Petroleum：Natural Sciences Edition，1994，18（2）：43—46．

［16］ 范子菲，方宏長．裂縫性油藏水平井穩(wěn)態(tài)解產(chǎn)能公式研究［J］．石油勘探與開發(fā)，1996，23（3）：52—57．Fan Zifei，F(xiàn)ang Hongchang．A steady solution formula of horizontal well productivity in a fractured reservoir［J］．Petroleum Explo—ration and Development，1996，23（3）：52—57．

［17］ 丁一萍．壓裂水平井產(chǎn)能研究［D］．北京：中國地質(zhì)大學(xué)（北京），2006．Ding Yiping．Analysis for productivity of horizontal fractured well［D］．Beijing：China University of Geosciences（Beijing），2006．

［18］ 丁一萍，王曉冬，邢靜，等．一種壓裂水平井產(chǎn)能計算方法［J］．特種油氣藏，2008，15（2）：64—68．Ding Yiping，Wang Xiaodong，Xing Jing，et al．A method of productivity calculation for fractured horizontal well［J］．Special Oil ＆Gas Reservoirs，2008，15（2）：64—68．

［19］ 寧正福，韓樹剛，程林松，等．低滲透油氣藏壓裂水平井產(chǎn)能計算方法［J］．石油學(xué)報，2002（2）：68—71．Ning Zhengfu，Han Shugang，Cheng Linsong，et al．Productivity calculation method of fractured horizontal wells in low permeability oil or gas field［J］．Acta Petrolei Sinica，2002（2）：68—71．

［20］ 程林松．高等滲流力學(xué)［M］．北京：石油工業(yè)出版社，2011：55—56．Cheng Linsong．Advanced seepage mechanics［M］．Beijing：Petroleum Industry Press，2011：55—56．

TE357．11

A

2095—4107（2013）04—0055—05

DOI 10．3969／j．issn．2095—4107．2013．04．008

2013—07—01；編輯：關(guān)開澄

中石化科技攻關(guān)項(xiàng)目（P12099）

鄭松青（1982—），男，博士，工程師，主要從事油藏工程方面的研究．