姜 永 別旭偉 劉洪洲 王 迪 吳浩君

(中海石油(中國)有限公司天津分公司 天津 300459)

雙重介質(zhì)油藏注水井試井解釋模型的建立及應(yīng)用*

姜 永 別旭偉 劉洪洲 王 迪 吳浩君

(中海石油(中國)有限公司天津分公司 天津 300459)

姜永,別旭偉,劉洪洲，等.雙重介質(zhì)油藏注水井試井解釋模型的建立及應(yīng)用[J].中國海上油氣，2017,29(4)：98-103.

JIANG Yong,BIE Xuwei,LIU Hongzhou,et al.Establishment and application of a well test interpretation model for injection wells in dual porosity reservoir [J].China Offshore Oil and Gas,2017,29(4):98-103.

針對雙重介質(zhì)油藏注水開發(fā)特征，基于Buckley-Leverett飽和度分布方程建立了注水井壓降試井解釋模型，通過Laplace變換和Stehfest數(shù)值反演對模型進(jìn)行了求解，利用無因次壓力降落解獲得了注水井的壓力降落典型圖版。分析圖版可知，兩相區(qū)對注水井壓力響應(yīng)特征較大，油水黏度差異會導(dǎo)致雙對數(shù)壓力導(dǎo)數(shù)曲線上翹，并且兩者差異越大，兩相區(qū)雙對數(shù)壓力導(dǎo)數(shù)曲線上翹幅度越大，因此在試井解釋時應(yīng)區(qū)別儲層物性變化或不滲透邊界造成的壓力導(dǎo)數(shù)曲線上翹。選取JZ油田1口注水井進(jìn)行了實例應(yīng)用，結(jié)果表明本文解釋圖版擬合結(jié)果與現(xiàn)場測試數(shù)據(jù)吻合較好，并且解釋所得的儲層參數(shù)合理。本文建立的試井解釋模型可以計算注水前緣的位置，分析注水井周圍地層信息及邊界情況，對評價雙重介質(zhì)油藏注水井的動態(tài)及注水方案的設(shè)計具有重要的指導(dǎo)意義。

雙重介質(zhì)油藏；注水井；試井；解釋模型；壓力降落典型圖版

注水井試井可以獲取油相和水相滲透率、注水井周圍污染情況、注水前緣距離、油藏邊界等信息，為注水井措施提供依據(jù)。關(guān)于注水井試井解釋方法，國內(nèi)外許多學(xué)者開展了相關(guān)研究[1-4]。于九政 等[5]針對低滲透油藏注水開發(fā)特點，在考慮啟動壓力梯度方向性和時效性的情況下，提出了低滲儲層注水井井下關(guān)井試井壓力響應(yīng)機(jī)理，建立了低滲儲層井下關(guān)井條件下注水井的試井解釋模型。劉佳潔 等[6]針對注水井試井未考慮含水率、水驅(qū)前緣等參數(shù)的問題，根據(jù)質(zhì)量守恒原理，結(jié)合注水井巖心歸一化相滲曲線，在油水兩相區(qū)運用油水兩相流滲流理論，建立了三區(qū)流體復(fù)合注水井試井模型。李乃華[7]建立了多層均質(zhì)定壓邊界灌注試井?dāng)?shù)學(xué)模型，針對適于現(xiàn)場應(yīng)用的折線型灌注函數(shù)，獲得了層間及井底壓降的實空間精確解析表達(dá)式。張艷玉 等[8]假設(shè)測試時間內(nèi)地層壓力均勻分布，油藏為兩區(qū)復(fù)合儲層模型，不考慮每個區(qū)中飽和度梯度，建立了單一介質(zhì)注水井壓降試井解釋模型。但目前有關(guān)雙重介質(zhì)油藏注水井壓降試井解釋方法的研究較少。

為了更好地掌握雙重介質(zhì)地層滲流動態(tài)信息，為注水井的增注措施提供依據(jù)，本文考慮地層中飽和度梯度和注水階段對壓力降落試井的影響，建立了雙重介質(zhì)油藏注水井壓降試井解釋模型，并進(jìn)行了實例應(yīng)用，結(jié)果表明本文解釋圖版擬合結(jié)果與現(xiàn)場測試數(shù)據(jù)吻合較好，解釋所得的儲層參數(shù)合理可靠。

1 模型的建立及求解

1.1 注水階段的物理模型

注水井物理模型見圖1。

圖1 注水井物理模型示意圖Fig .1 Schematic of water injection well model

為研究方便，對模型作如下基本假設(shè)：

1) 地層水平等厚，油藏由裂縫系統(tǒng)和基質(zhì)系統(tǒng)組成；

2) 注水井以恒定速率注水；

3) 流體和巖石微可壓縮，且壓縮系數(shù)為常數(shù)；

4) 流體在地層等溫流動，滲流滿足達(dá)西定律；

5) 忽略重力和毛細(xì)管力影響；

6) 考慮井儲和表皮系數(shù)影響；

7) 水驅(qū)油為非活塞式的，裂縫系統(tǒng)和基質(zhì)系統(tǒng)飽和度分布滿足Buckley-Leverett方程。

1.2 注水井試井解釋模型

注水階段的物理模型假設(shè)考慮了油藏內(nèi)油水飽和度的分布，為方便求解，將油水兩相區(qū)剖分為N個圓環(huán)(圖2)，即將連續(xù)的飽和度分布用階梯型分布的飽和度代替，當(dāng)圓環(huán)個數(shù)趨近于無限大時，階梯型飽和度分布趨于連續(xù)的飽和度分布(圖3)。

圖2 油水兩相區(qū)剖分為N個圓環(huán)的模型示意圖Fig .2 Schematic of the post split model with Ncircular rings

圖3 飽和度分布示意圖Fig .3 Schematic of saturation distribution

剖分后每一圓環(huán)內(nèi)的飽和度為一定值，即可認(rèn)為每一圓環(huán)內(nèi)的水驅(qū)油是活塞式的，其飽和度分布可由Buckley-Leverett方程確定，即

(1)

每一圓環(huán)內(nèi)的飽和度為

(2)

(3)

(4)

將每個圓環(huán)內(nèi)的Swi代入Buckley-Leverett方程，則可計算出每一圓環(huán)無因次半徑rDi。

對于每一圓環(huán)，其擴(kuò)散方程為

裂縫系統(tǒng)

(5)

基質(zhì)系統(tǒng)

(6)

初始條件

(7)

內(nèi)邊界條件

(8)

(9)

外邊界條件

1) 無限大外邊界

(10)

2) 定壓外邊界

(11)

3) 封閉外邊界

(12)

銜接條件

(13)

(14)

1.3 數(shù)學(xué)模型的求解

定義拉普拉斯變換

(15)

則數(shù)學(xué)模型可化為

裂縫系統(tǒng)

(16)

基質(zhì)系統(tǒng)

(17)

整理得

(18)

內(nèi)邊界條件

(19)

(20)

銜接條件

(21)

(22)

外邊界條件

1) 定壓外邊界

(23)

2) 封閉外邊界

(24)

3) 無限大外邊界

(25)

對于上述數(shù)學(xué)模型，其方程為一組零階的虛宗量Bessel方程，通解為

(27)

zi=(ln2/t)i

(28)

(29)

2 典型試井曲線圖版分析

通過求解數(shù)學(xué)模型可以獲得無因次井底壓力pwD與無因次時間tD的關(guān)系。假設(shè)關(guān)井測試階段水驅(qū)前緣不再發(fā)生變化，利用疊加原理得到關(guān)井階段的壓力降落解為

(30)

利用無因次壓力降落解可獲得注水井的壓力降落典型圖版，如圖4所示，計算參數(shù)為CD=50、S=5、λ=0.000 01、ω=0.02。根據(jù)導(dǎo)數(shù)曲線特征可劃分7個流動段,分別對應(yīng)井筒儲集效應(yīng)階段、表皮效應(yīng)階段、基質(zhì)向裂縫竄流階段、水區(qū)徑向流階段、油水兩相區(qū)響應(yīng)階段、總系統(tǒng)徑向流段和邊界響應(yīng)階段。從圖4可以看出，在竄流階段結(jié)束后出現(xiàn)第1條水平直線段，表征的是井周圍水區(qū)的響應(yīng)；之后由于油水黏度的差異，壓力導(dǎo)數(shù)曲線呈逐漸上翹的特征，反映的是油水兩相區(qū)的響應(yīng)，并且油水黏度差異越大，兩相區(qū)雙對數(shù)曲線上翹幅度越大(圖5)。另外，水驅(qū)前緣距離越遠(yuǎn)，油水兩相區(qū)時間越長，雙對數(shù)曲線上到達(dá)前緣壓力響應(yīng)時間越晚(圖6)，但不影響由于流體性質(zhì)差異導(dǎo)致的曲線上翹幅度。因此，在試井解釋時應(yīng)區(qū)別儲層物性變化或不滲透邊界造成的壓力導(dǎo)數(shù)曲線上翹，在到達(dá)水驅(qū)前緣后導(dǎo)數(shù)曲線上出現(xiàn)一個凸起，接著導(dǎo)數(shù)曲線出現(xiàn)第2條水平直線段，反映油區(qū)的壓力響應(yīng)特征。

圖4 雙重介質(zhì)油藏注水井壓力降落典型曲線Fig .4 Typical pressure drop curves of an injection well in fractured reservoir

圖5 油水黏度差異對雙重介質(zhì)油藏注水井壓力降落 典型曲線的影響Fig .5 Effect of oil viscosity difference on typical pressure drop curves of an injection well in fractured reservoir

圖6 水驅(qū)前緣對雙重介質(zhì)油藏注水井壓力降落 典型曲線的影響Fig .6 Effect of water flooding front on typical pressure drop curves of an injection well in fractured reservoir

3 實例應(yīng)用

JZ油田一口注水井X井，在進(jìn)行壓力降落測試前平均日注水256 m3，注水2個月累注水1.54×104m3。該注水井及油藏基本參數(shù)見表1，地層流體的相對滲透率曲線及地層分流率曲線見圖7。

該注水井實測壓力數(shù)據(jù)與本文典型圖版擬合情況見圖8，可以看出，兩者吻合較好。根據(jù)壓力擬合可計算出殘余油飽和度下的水相滲透率，根據(jù)時間擬合可計算出前緣的位置[8]。該注水井最終解釋參數(shù)值見表2，可以看出，儲層物性較好，注水井存在一定程度污染，水驅(qū)波及距離達(dá)到183 m。

表1 注水井X井及油藏基本參數(shù)Table 1 Basic parameters of the injection Well X and the reservoir

圖7 注水井X井相對滲透率曲線和地層分流率曲線Fig .7 Relative permeability curve and water diversion rate curve of the injection Well X

圖8 注水井X井壓降曲線與典型圖版擬合圖Fig .8 Fitting chart of pressure drop curves of the injection Well X with typical curves表2 注水井X井試井解釋結(jié)果Table 2 Well test interpretation results of the injection Well X

參數(shù)數(shù)值參數(shù)數(shù)值地層滲透率/mD382井儲系數(shù)/(m3·MPa-1)0.15油相滲透率/mD163竄流系數(shù)0.0003水相滲透率/mD206彈性儲能比0.16前緣距離/m183地層壓力/MPa15.0表皮系數(shù)6.3

4 結(jié)論

1) 基于Buckley-Leverett飽和度分布方程建立了雙重介質(zhì)油藏注水井壓降試井解釋模型，通過Laplace變換和Stehfest數(shù)值反演對模型進(jìn)行了求解，得到了注水井的壓力降落典型圖版。分析圖版可知，兩相區(qū)對注水井壓力響應(yīng)特征較大，油水黏度差異會導(dǎo)致雙對數(shù)曲線上翹，并且油水黏度差異越大，兩相區(qū)雙對數(shù)曲線上翹幅度越大，因此在試井解釋時應(yīng)區(qū)別儲層物性變化或不滲透邊界造成的壓力導(dǎo)數(shù)曲線上翹。

2) 利用本文建立的試井解釋模型可以計算注水前緣的位置，分析注水井周圍地層信息及邊界情況。實例應(yīng)用表明，本文解釋圖版擬合結(jié)果與現(xiàn)場測試數(shù)據(jù)吻合較好，解釋所得的儲層參數(shù)合理可靠。

符號注釋

Sw—含水飽和度，%；

rD(Sw)—Sw所對應(yīng)的無因次距離；

q—注水量，m3/d；

tp—注水時間，h；

h—油藏厚度，m；

rw—井半徑，m；

Sor—殘余油飽和度，%；

SwF—水驅(qū)油前緣飽和度，%；

Kf—裂縫滲透率，D；

K—地層平均滲透率，D；

λ—竄流系數(shù)；

φf—裂縫孔隙度,%；

φ—地層孔隙度,%；

Ctf—裂縫系統(tǒng)壓縮系數(shù),1/MPa；

Ct—地層綜合壓縮系數(shù),1/MPa；

μ—流體黏度,mPa·s；

re—油藏半徑，m；

B—體積系數(shù),m3/m3；

t—時間，h；

S—表皮系數(shù),無因次；

CD—無因次井筒儲存系數(shù)；

pi—原始地層壓力,MPa；

z—拉普拉斯變量；

I0、K0—第1類、第2類零階修正貝塞爾函數(shù)；

Ai、Bi—待定常數(shù)；

rDf—水驅(qū)油前緣無因次距離。

[1] 尹洪軍,付春權(quán),萬貴春，等.多重復(fù)合低滲透油藏注水井壓力動態(tài)分析[J].石油學(xué)報,2002,23(4):70-74.YINHongjun,FUChunquan,WANGuichun,etal.Pressuretransientanalysisofwaterinjectionwellinmultiplecompositelowpermeabilityreservoir[J].ActaPetroleiSinica,2002,23(4):70-74.

[2] 巨亞鋒，于九政，郭方元，等.低滲兩區(qū)復(fù)合油藏注水井試井解釋模型[J].石油天然氣學(xué)報(江漢石油學(xué)院學(xué)報),2012,34(7):109-113.JU Yafeng，YU Jiuzheng，GUO Fangyuan，et al.Well-test interpretation model for water injection wells in two-regional low permeability composite reservoirs [J].Journal of Oil and Gas Technology(J.JPI),2012,34(7):109-113.

[3] 尹洪軍，曾昭英，劉義坤，等.合理利用典型曲線分析注水井試井資料的方法[J].大慶石油學(xué)院學(xué)報,1993,17(2):112-114.YIN Hongjun,ZENG Zhaoying，LIU Yikun，et al.A method for analyzing well test data of water injection wells by using type curves reasonably [J].Journal of Daqing Petroleum Institute,1993,17(2):112-114.

[4] 劉佳潔,孟英峰,李皋，等.基于油水兩相滲流的地層流體復(fù)合注水井試井模型[J].石油天然氣學(xué)報(江漢石油學(xué)院學(xué)報),2014,36(2):128-132.LIU Jiajie,MENG Yingfeng，LI Gao，et al.Composite well test model based on oil-water two-phase flow for injection wells[J].Journal of Oil and Gas Technology(J.JPI),2014,36(2):128-132.

[5] 于九政，晏耿成，巨亞鋒，等.低滲油藏注水井試井壓力響應(yīng)機(jī)理及解釋模型[J].斷塊油氣田,2013,20(1):92-95.YU Jiuzheng，YAN Gengcheng，JU Yafeng，et al.Response mechanism and interpretation model of water injection well testing pressure in low permeability reservoir [J].Fault-Block Oil and Gas Field,2013,20(1):92-95.

[6] 劉佳潔,孟英峰,李皋，等.三區(qū)流體復(fù)合注水井試井模型研究[J].石油天然氣學(xué)報(江漢石油學(xué)院學(xué)報),2013,35(2):120-125.LIU Jiajie,MENG Yingfeng,LI Gao，et al.Three-zone fluid composite well test model of injection wells[J].Journal of Oil and Gas Technology(J.JPI),2013,35(2):120-125.

[7] 李乃華.多層定壓邊界灌注井?dāng)?shù)學(xué)模型[J].河北工業(yè)大學(xué)學(xué)報,2006,35(5):32-35.LI Naihua.Mathematical model of the injection well testing for layers with constant pressure[J].Journal of Hebei University of Technology,2006,35(5):32-35.

[8] 張艷玉,姚軍.現(xiàn)代試井解釋原理與方法[M].北京:中國石油大學(xué)出版社,2006.

[9] 賈永祿,聶仁仕,王永恒，等.二次梯度非線性滲流問題求解[J].石油天然氣學(xué)報(江漢石油學(xué)院學(xué)報),2008,30(4):119-123.JIA Yonglu,NIE Renshi,WANG Yongheng,et al.Solution for the problem of quadratic gradient nonlinear percolation[J].Journal of Oil and Gas Technology(J.JPI),2008,30(4):119-123.

[10] 劉永良,劉彬，劉大偉，等.三重介質(zhì)氣藏單井生產(chǎn)動態(tài)分析模型研究[J].石油天然氣學(xué)報(江漢石油學(xué)院學(xué)報)，2011,33(2):123-127.LIU Yongliang,LIU Bin，LIU Dawei，et al.Analytic model of single well production performance in triple porosity medium reservoirs [J].Journal of Oil and Gas Technology(J.JPI),2011,33(2):123-127.

[11] 胡建國.圓形定壓邊界油藏試井分析方法[J].中國海上油氣(地質(zhì)),1997,11(5):355-363.HU Jianguo.A method of well-test analysis for a circular reservoir with a constant pressure outer boundary [J].China Offshore Oil and Oas(Geology),1997,11(5):355-363.

[12] 廖新維,陳欽雷.均質(zhì)油藏斜井試井模型研究[J].中國海上油氣(地質(zhì)),1998,12(5):332-338.LIAO Xinwei,CHEN Qinlei.Discussion on well test model of homogeneous reservoir in slanted well[J].China Offshore Oil and Gas(Geology),1998,12(5):332-338.

[13] 賈永祿,趙必榮.拉普拉斯變換及數(shù)值反演在試井分析中的應(yīng)用[J].天然氣工業(yè),1992,12(1):60-64.JIA Yonglu,ZHAO Birong.Application of Laplace transform and numerical inversion to well test analysis [J].Natural Gas Industry,1992,12(1):60-64.

(編輯：楊 濱)

Establishment and application of a well test interpretation model for injection wells in dual porosity reservoir

JIANG Yong BIE Xuwei LIU Hongzhou WANG Di WU Haojun

(TianjinBranchofCNOOCLtd.,Tianjin300459,China)

Based on water flooding characteristics of dual porosity reservoir and Buckley-Leverett equation, the injection pressure drop testing interpretation model is established.Laplace transform and Stehfest numerical inversion are used to solve the model and a pressure drop typical plate of injection well is obtained with dimensionless pressure drop solution.Results show that oil-water two-phase region has large pressure response characteristics, oil-water viscosity difference causes double logarithmic curve upward.The larger the difference of oil-water viscosity, the more two-phase double logarithmic curve upturns.So the upward curve caused by reservoir change or impermeability boundary should be distinguished.The application in an injection well of JZ oilfield shows that the interpretation results agree with field test data.The proposed model can be used to calculate the position of water flood front, and to analyze the formation around the injection well and boundary condition, thus guiding water injection well evaluation and the design of injection scheme.

dual porosity reservoir; water injection well; well test; interpretation model; typical chart of pressure drop

*國家自然科學(xué)基金重點項目“儲層裂縫形成機(jī)理(編號：40772089)”部分研究成果。

姜永，男，工程師，碩士，主要從事試井、油藏工程方面的研究工作。地址：天津市濱海新區(qū)海川路2121號渤海石油管理局B座(郵編：300459)。E-mail：jiangyong198786858@163.com。

1673-1506(2017)04-0098-06

10.11935/j.issn.1673-1506.2017.04.012

TE344

A

2017-01-12 改回日期：2017-03-31