裂缝型底水气藏水侵动态分析方法

熊钰 杨水清 乐宏 唐建荣 余翔

1.西南石油大学石油工程学院 2.中国石油西南油气田公司蜀南气矿

裂缝型底水气藏水侵动态分析方法

熊钰1杨水清2乐宏2唐建荣2余翔2

1.西南石油大学石油工程学院 2.中国石油西南油气田公司蜀南气矿

对有水气藏水侵动态和储量核实的研究是气藏开发动态分析的基础工作,而针对裂缝型底水气藏相关分析方法的研究还较少,如何利用底水气藏的开发动态数据来研究水体性质及水侵动态,对底水气藏的开发后期实施有针对性的排水采气措施具有重要的实际意义。在参照底水气藏数学模型的基础上,结合针对边水气藏的A IF模型,运用水体影响函数理论,建立了底水驱A IF和压力动态分析模型,得到水体影响下的裂缝型底水气藏水侵动态分析方法,并引入非线性最优化法以获得气藏的水侵动态。利用Visual Basic语言编制气藏水侵动态分析程序,并依据一个实际典型裂缝型底水气藏——川南威远气田的生产数据进行了水侵动态验证分析。实例研究结果表明,该方法可以较为准确获得裂缝型底水气藏动态储量、水体性质和水侵动态,是裂缝型底水气藏水侵动态分析的好方法。

底水气藏 水侵 动态分析 底水驱 裂缝(岩石) 数学模型 威远气田

0 引言

在有水气藏的水侵动态分析研究中,针对边水气藏的水体影响函数法[1-2]对于边水气藏来说,拟合效果比较好,拟合压力的误差较小,能够比较真实客观地反映边水气藏的压力动态变化,得到较为准确的气藏动态储量和水侵动态数据。但是对于具有大裂缝构造或是底水入侵的气藏来讲,边水气藏水体影响函数的压力拟合和储量计算结果会产生较大误差,影响拟合和计算的精度。针对裂缝型底水气藏,有必要建立一种新的适用模型结合水体影响函数法来计算其动态储量以及在开发中气藏随时间变化的水侵量和水侵速度,进而得到裂缝型底水气藏的水侵动态特征。

1 底水驱模型的建立

由于裂缝型底水气藏的水侵多表现为边底水复合侵入的方式,笔者的研究对象均是外边界封闭的气藏,故可将其简化为一种底水模型来研究气水界面处压力动态表征。

对具有底水的裂缝型气藏可简化为由含气区和含水区构成,其底水驱模型[3-4]如图1所示,气层和水层都是圆柱形的,气层位于水层的上方,并假定水体区域的厚度(ha)和半径(ra)都是常数,气层的半径为 rR,水流入气层的界面半径r为0～rR,气层和水层的其他外边界是封闭的,于是可以作如下假设:

图1 裂缝型底水气藏底水驱模型示意图

1)水层孔隙度和厚度为定值;

2)水层渗透率是均质的,但具有方向性,水平方向和垂直方向的渗透率分别是 KH和 KV;

3)岩石和水的压缩系数为常数;

4)气层和水层是同心的,水层半径为 rR,气层半径为ra;

5)气水界面处压力已知; 6)水层岩石完全饱和水。

根据已建立的裂缝型底水气藏底水驱模型,考虑气水界面处的多级压力降,采用径向流坐标系中二维扩散方程表示水流从水层流向气层,定义渗流方程为:

其中:Δp(t,r,z)=pai-p(t,r,z)

初始和边界条件为:

使用拉普拉斯变换和变量分离对式(1)求解,此外使用压力平均法求取气水界面处的压力。

则在拉普拉斯空间,由定水侵速度在气水界面处产生的压力降为:

其中

s是与真实时间相对应的拉普拉斯变量;λn是由J1(λnra)=0解出的特征值。

2 压力动态描述

在研究中,处理底水影响下气藏的压力动态特征沿用水体影响函数法的思路[1],由杜哈美原理有:

根据水体影响函数法,有

对比式(4)、式(5),可知水体影响函数为:

由于事实上Δpxu(t)是式(3)在实空间下的解析解,因此反演式(3)可得:

当k的截断项定为1,这样 F(t)简化为:

其中

为了采用最优化的方法求解压力差,可将式(8)改写为:

故对于裂缝型底水气藏,气藏压力和函数 F(t)的关系可表示为:

3 水体影响函数的最优问题及其解

式(10)为一非线性方程,可构成一个最优化问题(即为非线性最优化拟合问题)。当引入牛顿迭代和最速下降原理相结合可采用列维贝格算法计算偏差函数,当最优拟合的偏差函数达到最小以后的解为最优解[1]。最优化问题为:

在实际应用中,式(11)中的变量水侵速度 ∑qv是利用两个生产时间间隔中的平均水侵速度来替代,即

4 实例分析

威远气田地处四川盆地南部威远、荣县和资中3县境内,1964年发现埋深3 000 m的震旦系灯影组气藏,是一个大型致密白云岩裂缝—孔洞型底水气藏[5]。其震旦系顶界含气范围为216 km2,含气高度为244 m,原始气水界面海拔为-2 434 m,原始地层压力为29.533 MPa,地层温度为120℃,基质平均孔隙度在2%左右,渗透率介于0.001×10-3～0.05×10-3μm2[6]。

通过统计和整理威远气田1965年至2007年的生产数据和平均地层静压测试数据,选取威远气田15个基本生产动态数据点作为拟合目标(见表1)。

表1 威远气田生产动态数据点表

威远气田底水模型建立所需的所有参数均可从《威远气藏描述暨地质模型研究报告》中获得,所涉及的参数整理可得:气层半径(rR)为8 292 m(由含气面积216 km2计算得到),水层半径(ra)为10 202 m(由含水面积327 km2计算得到),水层厚度(ha)为132 m,水层平面渗透率(KH)为6×10-3μm2,水层垂向渗透率(KV)为4×10-3μm2,水层平均孔隙度(φa)为0.02,水的地下黏度(μw)为0.277 65 m Pa·s,岩石有效压缩系数(cea)为14.23×10-4M Pa-1(由 Hall图版法算出)。

应用本方法对威远气田生产动态数据进行计算。通过计算得到储量与偏差函数的关系,当满足偏差函数达到最小时(对应的偏差函数值为0.587 6)的解为最优解。最优拟合动态地质储量为380×108m3、采出程度为37.86%。

根据计算得到的储量值,可得到威远气田最优拟合数据及压力拟合情况(图2、表2)。

图2 威远气田压力拟合情况图

表2 威远气田最优拟合数据表

图3 威远气田水侵速度和累计水侵量对比图

图3是在满足最优拟合偏差函数最小时,在水侵动态分析和生产动态数据的基础上得到的水侵速度和累积水侵量与生产时间的关系。水侵速度呈台阶式上升,这和开采阶段密切相关。

研究威远气田各个生产时段水侵速度与生产曲线的对比(图4)可知,威远气田生产初期的采气速度较大,整体水侵速度较小,在开始提高采气速度后水侵逐渐加剧。当气田开采进入中期,采气速度开始逐渐减小时,产水速度开始逐渐增大,并呈台阶式波动,同时水侵速度也呈现出一种阶梯式的递增,表明威远气田受水侵的影响是在生产中期逐渐加大的,采用气田内排水的方式加剧了裂缝性水窜;随着气藏开发进入生产末期,气田内生产井大量排水形成了水侵通道,水侵速度到达最大值,造成许多气井在短时间内大量见水,但由于底水水体的弹性能量大量释放以后,水侵能量减弱,水侵速度已经开始表现出逐渐减小的趋势。

图4 威远气田采气、产水速度与水侵速度对比示意图

另一方面,从函数 F(t)随时间的变化关系曲线(图5)可知,在威远气田开井投产初期,由于采气速度高、底水沿裂缝窜进侵入气藏、水体能量表现为逐渐增大;在进入生产中期时,水体的弹性能量开始逐渐减弱,水体的膨胀程度开始逐渐减小,进入生产末期后,使函数 F(t)形成一种持续下降甚至为负值,表明随着水体能量的大量释放以及受威远气田大面积关井影响,部分地层水有可能出现沿大裂缝系统回退的现象,说明整个水体能量表现为具有可排性,如果在合适的位置打排水井大量排水,可使气藏的采收率得到大幅度地提高。

图5 威远气田水体影响函数 F(t)随时间的变化关系图

5 结论

1)建立了裂缝型底水气藏的底水驱模型和水侵动态分析方法,利用威远气田的生产数据成功地分析了水侵动态。

2)改进的底水气藏水体影响函数法是裂缝型底水气藏水侵动态分析的有效方法。

3)底水气藏由于构造封闭造成的构造局部封存水体的水体影响函数同样为抛物线形。

4)利用该模型可以较好地解决裂缝型底水气藏水侵动态分析中存在的问题,得到的分析结果符合实际生产情况。

5)威远气田地层水具有可排性,合适的排水井大量排水,可大大提高气藏的采收率。

符 号 说 明

Δp为压降,M Pa;pai为水层原始地层压力,M Pa;Δpv为变压降,M Pa;η为水层扩散常数;r为径向空间变量,m;z为垂向空间变量,m;KV为水层垂向渗透率,10-3μm2;KH为水层水平方向渗透率,10-3μm2;t为时间,mon;φa为水层孔隙度,小数; μw为地下水黏度,m Pa·s;cea为水体有效压缩系数,MPa-1;ha为水体厚度,m;qx为气水界面处定流量,m3/mon;qv为气水界面处变流量,m3/mon;rR为气层半径,m;ra为水层半径,m;We为累计水侵量,m3。

[1]熊钰,胡述清,曲林,等.非均质气藏局部封存水体性质及水侵动态分析方法研究[J].天然气工业,2004,24(2):78-81.

[2]张新征,张烈辉,李玉林.预测裂缝型有水气藏早期水侵动态的新方法[J].西南石油大学学报,2007,29(5):82-85.

[3]YILD IZ T,KHOSRAV I A.An analytical bottom waterdrive aquifermodel fo rmaterial balance analysis[J].Journal of SPE Reservoir Evaluation&Engineering,2007,10(6): 618-628.

[4]张建业,熊钰,黄兰,等.底水驱解析新模型[J].国外油田工程,2007,23(6):13-16.

[5]张虎权,卫平生,张景廉,等.也谈威远气田的气源[J].天然气工业,2005,25(7):4-7.

[6]冯异勇,贺胜宁.裂缝性底水气藏气井水侵动态研究[J].天然气工业,1998,18(3):40-44.

A new method of water influx performance analysis on fracture gas reservoirs with bottom water

Xiong Yu1,Yang Shuiqing2,Le Hong2,Tang Jianrong2,Yu Xiang2
(1.College of Petroleum Engineering,Southw est Petroleum University,Chengdu,Sichuan 610050,China; 2.South Gas Field,PetroChina Southw est Oil&Gasf ield Com pany,L uzhou,Sichuan 646001,China)

NATUR.GAS IND.VOLUM E 30,ISSUE 1,pp.61-64,1/25/2010.(ISSN 1000-0976;In Chinese)

The studying of water invasion behavior and reserves substantiation is fundamentalwork for exp loitation perfo rmance analysis of aquifer gas reservoirs.However,there is still not enough related research on the fracture gas reservoirs w ith bottom water. Therefo re,using the p roduction data to study the water aquifer character and water invasion performance is very impo rtant for perfo rming themeasures of drainage gas recovery especially w hen the life of gas reservoir is coming to the declining period.Thus this paper p resents a new kind of aquifer influence function(A IF)of bottom aquifer driving gas reservoirs,the model of w hich is built up,in combination w ith the descrip tion of the p ressure performance characteristics,to get a new method of water influx performance analysison fracture gas reservoirs w ith bottom water and analyze the water invasion performance by the nonlinear op timization law. A s a case study,this app roach was used in the Weiyuan gas reservoir w hich is a classic fracture gas reservoir w ith bo ttom water. The results have indicated that this app roach can not only calculate the reserves exactly,but also analyze the aquifer behavio r and water invasion perfo rmance for those fracture gas reservoirs w ith bottom water.

gas reservoir w ith bottom water,water invasion,perfo rmance analysis,bottom aquifer drive,fracture,mathematical model,Weiyuan gas reservoir

四川省天然气开采重点实验室基金项目“裂缝性非均质气藏水侵动态预测方法”(编号:GPL04-07)。

熊钰,教授,1968年生;从事油气藏工程和流体相态方面的教学和研究工作。地址:(610500)四川省成都市新都区新都大道8号。电话:(028)83033082。E-mail:xiongyu-swpi@126.com

熊钰等.裂缝型底水气藏水侵动态分析方法.天然气工业,2010,30(1):61-64.

DO I:10.3787/j.issn.1000-0976.2010.01.017

(修改回稿日期 2009-11-16 编辑 韩晓渝)

DO I:10.3787/j.issn.1000-0976.2010.001.017

Xiong Yu,p rofesso r,was bo rn in 1968.He ismainly engaged in teaching and research of reservoir engineering and fluid phases.

Add:No.8,Xindu Rd.,Xindu District,Chengdu,Sichuan 610500,P.R.China

Tel:+86-28-8303 3082E-mail:xiongyusw pi@126.com