低滲水驅油藏水平井見水后產能研究新方法

袁　淋，李曉平，李松岑

（1.西南石油大學油氣藏地質及開發工程國家重點實驗室，成都610500；2.中國石化中原油田普光分公司，四川達州636155）

油氣田開發

低滲水驅油藏水平井見水后產能研究新方法

袁淋1，李曉平1，李松岑2

（1.西南石油大學油氣藏地質及開發工程國家重點實驗室，成都610500；2.中國石化中原油田普光分公司，四川達州636155）

在利用水平井開發低滲水驅油藏的過程中，油井見水是制約其產量的一個重要因素?；贘oshi水平井產能分析理論，將水平井三維滲流問題簡化為2個二維平面滲流問題，考慮啟動壓力梯度，利用保角變換方法分別獲得了2個二維平面內的產量公式，再結合水驅油藏中非活塞式驅替過程油水兩相滲流特征，最終求得了低滲水驅油藏中水平井產油量計算新公式。通過實例計算與對比發現，本文公式與水平井單相經典產能公式計算結果雖然均比實際產油量大，但是本文公式計算結果與實際產油量相對誤差最小，為2.94%，從而證明本文公式準確性較高、實用性較好。此項研究為低滲水驅油藏水平井見水后產量預測提供了新的思路。

低滲儲層；水驅開發；水平井；油水兩相；油井產能

0　引言

在利用水平井開發低滲水驅油藏過程中，由于黏度差異、毛細管現象以及油層非均質性等因素的影響，驅替方式通常為非活塞式驅替，且在驅替相與被驅替相之間存在一油水兩相滲流的過渡帶，加之啟動壓力梯度的影響，使得低滲水驅油藏中水平井周圍的滲流更加復雜，因而用常規水平井單相產能公式無法準確預測其產能［1-5］。國內外學者［6-14］將理論與實驗相結合，研究了低滲油藏中油水兩相滲流規律，但對于水驅油藏水平井見水后產能的研究尚存在不足。筆者認為水平井見水后外部滲流場為驅替相的單相滲流，內部滲流場為油水兩相滲流，進而在考慮啟動壓力梯度的基礎上，求得低滲水驅油藏水平井見水后產油量預測新公式，并通過實例證明新公式的準確性以及產水對水平井產能的影響，以期為低滲水驅油藏水平井合理工作制度的研究提供理論依據。

1　產能公式推導

流體在水平井周圍地層中的滲流為橢球形三維滲流。對于該三維滲流問題的處理，目前最常用的方法是Joshi［3，15］提出的將三維滲流簡化為水平平面與垂直平面2個二維滲流。在水驅油藏見水后，水平平面內的滲流為水相的滲流，而垂直平面內的滲流為油水兩相的滲流。在2個二維滲流平面內分別求取產能公式，再利用等值滲流阻力法與水電相似原理便可獲得水平井產能公式。

1.1水平平面內的產能公式

流體在水平平面內的滲流為橢圓形供給邊界中存在一排液坑道的滲流問題［5］，引入儒柯夫斯基變換，取保角變換函數，即

式中：L為水平段長度，m；Z與ξ分別為儒柯夫斯基變換前后2個復平面內的變量，Z=x+iy，ξ=u+iv。

通過此變換，將Z1平面內長半軸為a、短半軸為b的橢圓滲流場轉換成Z2平面內半徑為（a+b）/（L/2）的圓形區域，并將線段（-L/2，0）到（L/2，0）映射成單位圓周，這就等價于半徑為（a+b）/（L/2）的泄油區域內一口半徑為1的直井的滲流問題（圖1）。

經過儒柯夫斯基變換后，啟動壓力梯度變為

式中：λ′為保角變換后的啟動壓力梯度，MPa/m；re為泄油半徑，m；λ為啟動壓力梯度，MPa/m；a為橢圓形泄油區域長半軸長度，m；b為橢圓形泄油區域短半軸長度，m。

圖1　水平平面內保角變換示意圖Fig.1　The sketch map of conformal mapping in horizontal plane

水平平面內的滲流為單相水的流動，考慮啟動壓力梯度的運動方程為［16］

式中：v為油藏中任一點的滲流速度，m/d；kh為油藏水平方向滲透率，mD；μw為地層水黏度，mPa·s；p為油藏中任一點的壓力，MPa；r為油藏中任一點的徑向距離，m。

分離變量，將式（3）在r區間［1，（a+b）/（L/2）］上積分，化簡得到水平平面內的壓降為

式中：pe為油藏驅動壓力，MPa；pF為內部滲流場驅動壓力，MPa；h為油藏厚度，m；qL為日產液量，m3/d。

1.2垂直平面內的產能公式

垂直平面內的滲流可看成上下封閉邊界中一點匯的滲流問題，引入儒柯夫斯基變換，取保角變換函數

通過此變換，將Z1平面的帶型區域（-0.5h＜y＜0.5h）轉換成Z2平面的單位圓形區域，Z1平面上的匯點（0，0）變成Z2平面上的圓心（0，0）（圖2）。

圖2　垂直平面內保角變換示意圖Fig.2　The sketch map of conformal mapping in vertical plane

Z1平面上的井筒半徑rw在Z2平面上變成ρw

式中：rw為水平井井筒半徑，m；ρw為保角變換后的水平井井筒半徑，m。

啟動壓力梯度變為

垂直平面內為油水兩相滲流，考慮啟動壓力梯度的油水兩相運動方程為［17］

式中：kv為油藏垂直方向滲透率，mD；krw為水相相對滲透率；kro為油相相對滲透率；μo為原油黏度，mPa·s。

式中：pwf為水平井井底流壓，MPa。

在垂直平面內，油藏的含水飽和度分布滿足徑向Buckley-Levereet方程［16］：

式中：φ為儲層孔隙度，%；ΣQ1為累計產液量，m3；rf為垂直平面內供給半徑，m，文中為1 m；fw為油藏中任一點的分流率，%；Sw為油藏中任一點的含水飽和度，%。

當r為井筒半徑時，則此時位于油水驅替前緣，rf=1，則式（10）可以寫為

式中：Swf為前緣含水飽和度，%。

綜合式（10）和（11）可得

式（12）等號兩邊同時微分得

將式（13）代入式（9），并利用分流率的表達式

則式（9）可以變化為

1.3水平井產能公式

結合式（4）與（15），可得到水平井的產液量公式，即

利用產液量與產油量之間的關系，將式（16）轉化為產油量公式

其中

式中：qo為日產油量，m3/d。

A的求取方法有2種：一是利用數值積分法求解，根據梯形公式或辛普森公式進行計算；二是采用圖形擬合法求解。考慮到數值積分法計算結果精度不高，因此采用圖形擬合法求解，具體步驟如下：①根據兩相滲透率實驗測試數據做油水兩相相對滲透率曲線；②根據含水率計算公式，結合相滲曲線，做分流率隨含水飽和度變化的關系曲線，再過（Swr，0）點做該曲線的切線，切點對應的含水飽和度即為水驅前緣含水飽和度；③利用分流率曲線，計算求得分流率導數隨含水飽和度變化的關系曲線，擬合得到分流率導數與含水飽和度的表達式；④根據第①～②步分流率計算結果以及第③步分流率導數與含水飽和度擬合關系式，計算得到B與fw′（Sw）關系曲線，進而擬合得到二者之間關系式，將變量fw′（Sw）在區間［0，fw′（Swf）］上積分得到A的值。

2　實例計算與分析

某低滲水驅油藏處于開發中后期，日產液量較高，但油井含水率較高，產油量低，其中一水平井參數如表1所列。

表1　某低滲水驅油藏一水平井基本參數Table1　The parameters of a horizontal well in one low permeability water flooding reservoir

若要計算水平井產液量以及產油量，則必須先計算出A的值。根據所給相滲測試資料，做該區塊油水兩相相對滲透率曲線（圖3）。

圖3　油水兩相相對滲透率曲線Fig.3　The relative permeability curve of oil-water two phase

根據該相滲曲線，利用式（16）分流率計算公式獲得分流率曲線（圖4）。

圖4　分流率隨含水飽和度變化曲線Fig.4　The change of shunt index with water saturation

過（0.3，0）點做分流率曲線的切線，切點對應的含水飽和度為47%，即前緣含水飽和度，對應的前緣分流率約為0.7。根據該分流率曲線，通過計算可以得到分流率導數曲線（圖5）。

圖5　分流率導數曲線Fig.5　The derivative curve of shunt index

通過擬合發現，采用六次多項式擬合，相關系數為0.975 1，準確性較高，則得到分流率導數與含水飽和度的關系式為

結合式（19）、相滲曲線以及分流率曲線，做B與fw′（Sw）關系曲線（圖6）。

圖6　B與fw′（Sw）的擬合關系曲線Fig.6　The fitting curves between B and fw′（Sw）

為了提高擬合精度，可以進行分段擬合，從0增加到fw′（Sw）max為第一段，fw′（Sw）max遞減到fw′（Swf）為第二段。由擬合結果可以看出，相關系數為1，擬合程度較高，可以得到關系式為

fw′（Sw）max的值通過分流率導數曲線可以得到，約為5.5，對式（20）進行積分得到A值為84.305 2。將表1中的參數代入式（16）和（17），進而求得產液量為30.628 m3/d，產油量為8.06 m3/d。

將所給已知參數代入考慮啟動壓力梯度的水平井單相經典產能公式進行計算，并將計算結果（表2）與實際產油量7.83 m3/d進行對比。

表2　水平井單相經典產能公式與本文公式產油量計算結果對比Table2　Comparison of results calculated by new formula and single-phase productivity formula of horizontal well

從表2可以看出，本文公式與水平井單相經典產能公式計算結果雖然均比實際產油量偏高，但本文公式計算結果與實際產油量相對誤差最小，為2.94%。這是因為水平井見水后，井筒周圍滲流場由純油區變為油水兩相區，滲流阻力增加，進而導致產油量比水平井單相經典產能公式計算結果偏低，說明水平井井底見水對其產油量具有較大的影響。因此，水驅油藏開發過程中應采取合理控水措施，以延緩水的突破，延長無水采油期。

3　結論

（1）在考慮啟動壓力梯度的基礎上建立了低滲水驅油藏滲流數學模型，求得了水平井產液量以及產油量公式，實例計算表明本文公式與水平井單相經典產能公式計算結果雖然均比實際產油量偏高，但是本文公式計算結果更加接近實際產油量。

（2）在水驅油藏中，由于井底水體的突破，使得井筒附近的滲流變為油水兩相滲流，并且含水飽和度沿徑向呈非均勻分布，而遠井地帶仍為單相水的滲流。

（References）：

［1］Borisov J P.Oil production using horizontal and multiple deviation wells［M］.Oklahoma：The R&D Library Translation，1984.

［2］Giger F M.Horizontal wells production techniques in heterogenous reservior［R］.SPE 13710，1985.

［3］Joshi S D.Augmentation of well production using slant and horizontal wells［R］.SPE 15375，1986.

［4］Renard G，Dupuy J M.Formation damage effects on horizontal well flow efficiency［R］.SPE 19414，1991.

［5］陳元千.水平井產量公式的推導與對比［J］.新疆石油地質，2008，29（1）：68-71.

Chen Yuanqian.Derivation and correlation of production rate formula for horizontal well［J］.Xinjiang Petroleum Geology，2008，29（1）：68-71.

［6］Wu Yushu.Theoretical studies on non-Newtonian and Newtonian fluid flow through porous Media［D］.Berkeley：University of California，1990.

［7］孫衛，何娟.姬塬延安組儲層水驅油效率及影響因素［J］.石油與天然氣地質，1999，20（1）：26-29.

Sun Wei，He Juan.Displacing efficiency and affecting factors of reservoirs in Yan'an formation，Jiyuan region［J］.Oil&Gas Geology，1999，20（1）：26-29.

［8］宋付權，劉慈群.含啟動壓力梯度油藏的兩相滲流分析［J］.石油大學學報：自然科學版，1999，23（3）：47-50.

Song Fuquan，Liu Ciqun.Two phase flow analysis in oil reservoir considering start-up pressure gradient［J］.Journal of the University of Petroleum，China：Edition of Natural Science，1999，23（3）：47-50.

［9］宋付權，劉慈群，張慧生.低滲透油藏中水平井兩相滲流分析［J］.水動力學研究與進展，2002，17（2）：181-188.

Song Fuquan，Liu Ciqun，Zhang Huisheng.Analysis of two-phase percolation of a horizontal well in low-permeability porous media［J］.Journal of Hydrodynamics，2002，17（2）：181-188.

［10］姚約東，葛家理，李相方.低滲透油藏油水兩相滲流研究［J］.石油大學學報：自然科學版，2005，29（2）：52-56.

Yao Yuedong，Ge Jiali，Li Xiangfang.Oil-water two-phase fluid flow in low permeability reservoir［J］.Journal of the University of Petroleum，China：Edition of Natural Science，2005，29（2）：52-56.

［11］劉曉彥，謝吉兵，廖建波，等.靖安油田ZJ2井區侏羅系油藏水驅狀況研究與調整對策［J］.巖性油氣藏，2007，19（1）：124-129.

Liu Xiaoyan，Xie Jiebin，Liao Jianbo，et al.Strategy and waterflood condition research on Jurassic reservoir of ZJ2 well-area，Jing'an Oilfield［J］.Lithologic Reservoirs，2007，19（1）：124-129.

［12］田敏，李曉平，李楠，等.低滲透油藏油水兩相非線性滲流規律［J］.特種油氣藏，2011，18（6）：92-95.

Tian Min，Li Xiaoping，Li Nan，et al.Theoretical study on oil-water two-phase non-linear flow in low permeability reservoirs［J］.Special Oil and Gas Reservoirs，2011，18（6）：92-95.

［13］何文祥，楊億前，馬超亞.特低滲透率儲層水驅油規律實驗研究［J］.巖性油氣藏，2010，22（4）：109-111.

He Wenxiang，Yang Yiqian，Ma Chaoya.Expermental study on waterflooding in ultra-low permeability reservoirs［J］.Lithologic Reservoirs，2010，22（4）：109-111.

［14］車洪昌，任耀宇，劉漢平，等.龍虎泡油田活性水驅油室內實驗研究［J］.巖性油氣藏，2011，23（2）：128-132.

Che Hongchang，Ren Yaoyu，Liu Hanping，et al.Laboratory study on oil displacement with active water in Longhupao Oilfield［J］. Lithologic Reservoirs，2011，23（2）：128-132.

［15］Joshi S D.A Review of horizontal well and drain hole technology［R］.SPE 16868，1988：339-352.

［16］李曉平.地下油氣滲流力學［M］.北京：石油工業出版社，2008：136-137.

Li Xiaoping.The underground oil and gas percolation mechanics［M］.Beijing：Petroleum Industry Press，2008：136-137.

［17］計秉玉.油水兩相徑向穩定滲流條件下壓力分布及產量變化特征［J］.大慶石油地質與開發，1995，14（3）：44-46.

Ji Bingyu.Under condition of radial stable seepage flow with two phase of oil and water［J］.Petroleum Geology and Development in Daqing，1995，14（3）：44-46.

（本文編輯：于惠宇）

New method for productivity study of horizontal well after water breakthrough in low permeability water flooding reservoir

YUAN Lin1，LI Xiaoping1，LI Songcen2
（1.State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation，Southwest Petroleum University，Chengdu 610500，China；2.Gas Production Plant，Puguang Branch，Zhongyuan Oilfield Company，Sinopec，Dazhou 636155，Sichuan，China）

During the development of low permeability water flooding reservoir with horizontal well，water breakthrough is an important factor to restrict the production.Based on the Joshi theory of productivity of horizontal well，this paper simplified3Dflowofhorizontalwellintotwo2Dflows；consideringthestartuppressuregradient，usedconformal mapping method to obtain productivity formulas in two 2D flow planes respectively；combined with the flow characteristics of oil-water two phase during non-piston-like displacement in water flooding reservoir，eventually obtained a new production formula of horizontal well in low-permeability water flooding reservoir.Through case study and comparison，the results calculated by new formula and single-phase productivity formula of horizontal well are more than the actual output，but the result calculated by new formula has a small relative error with actual output，only of 2.94%.It demonstrates that the new formula has a high degree of accuracy and is very practical.The study will provide a new method for predicting the production after water breakthrough of horizontal well in low permeability water flooding reservoir.

low permeability reservoir；developing with water flooding；horizontal well；oil-water two phase；productivityofoil well

TE328

A

1673-8926（2015）03-0127-06

2014-10-30；

2014-12-04

國家杰出青年科學基金項目“油氣滲流力學”（編號：51125019）資助

袁淋（1990-），男，西南石油大學在讀碩士研究生，研究方向為油氣藏工程與滲流力學。地址：（610500）四川省成都市新都區西南石油大學油氣藏地質及開發工程國家重點實驗室B403室。E-mail：yuanlin343@163.com

李曉平（1963-），男，教授，博士生導師，主要從事滲流力學、試井分析及油氣藏工程領域的教學和科研工作。E-mail：nclxphm@126.com。