作者简介:陈民锋(1971―),男,中国石油大学(北京)副教授、博士. 研究方向:油气田开发系统理论方法和提高采收率. E-mail:cmfllp96@126.com
中文责编:晨 兮; 英文责编:天 澜
Chen Minfeng, Wang Zhaoqi, Zhang Qichen, and Chen LuCollege Petroleum Engineering, China University of Petroleum(Beijing), Beijing 102249, P.R.China
oilfield development; low permeability reservoirs; threshold pressure gradient; pressure and velocity field; effective producing area; injector producer distance
DOI: 10.3724/SP.J.1249.2017.01091
低渗透储层具有明显的启动压力梯度,导致注采井间的渗流场分布与常规条件存在一定差异.基于典型注采条件下的渗流特点和驱替过程分析,在产量方程中引入启动压力梯度,构建出考虑启动压力影响下一注一采定产渗流模型.通过推导得出注采井间压力场、速度场分布的数学表达式,并建立计算注采井间有效驱替面积、有效注采井距的方法,得到不同开发条件下井间渗流场分布和注采单元中有效驱替范围变化的规律.结合典型低渗透油藏参数进行的研究表明,根据有效驱替要求确定有效注采井距,并以此为基础部署开发井网,能在更大程度上有效驱动注采单元内的流体,进而提高低渗透油藏的开发效果.
Due to the existence of significant threshold pressure gradient, the flow in low permeability reservoirs is low-speed non-Darcy flow. Under these conditions, interwell flow filed has significant difference compared to the fields of conventional reservoirs. In view of the flow characteristics and displacement process analysis of a typical water flooding regime, threshold pressure gradient is introduced into the basic flow equation and combined with the condition that both production and injection wells are working in a constant rate. A mathematical model is established for one injector producer pair in infinite formation. According to the initial and boundary conditions, interwell pressure and velocity field distribution function are deduced from the model. Based on the mathematical description, interwell pressure and velocity field distribution and variation could be defined. The analysis shows that for the waterflood development of low permeability oilfield, only by taking effective flooding strategy, can the oil between the injection and production wells be effectively displaced to a greater extent and can the development effectiveness of low permeability reservoirs be improved.
与常规油藏相比,低渗透储层具有较为明显的启动压力梯度,在水驱开发过程中,流体在多孔介质中的渗流偏离达西定律,导致在注采井间部分区域内,若驱替压力梯度小于启动压力梯度,流体不能流动,直接影响了低渗透油藏储量的水驱动用效率[1- 4].
为保证注采井间的注采压差大于流体流动所需的启动压差,确定开发井距通常采用的方法是,在给定最大压差条件下,计算井间驱替压力梯度等于启动压力梯度时对应的极限注采井距.但这样确定的开发井距在实际生产中存在一定的问题,一方面,先利用势叠加原理求得压力分布,再减去启动压力梯度得到的计算公式,没有考虑启动压力梯度对渗流过程的影响,是一个等效处理手段; 另一方面,在这种极限注采井距下建立的注采压差,只能使主流线附近区域内的流体能够流动(即注采连线的“线驱动”),而整个注采单元大部分区域的流体由于受到的注采压差较小,而不能得到有效动用[5-12].
启动压力梯度存在下的注采井间渗流场分布规律,与经典渗流理论得到的渗流场存在较大的差异,为了在注采井间有更多区域内的流体能够被驱替(即注采井间“面驱动”),需要从基本渗流方程入手,研究在不同条件下注采井间压力场和速度场分布及变化规律.
本研究基于典型注采条件下的渗流特点和驱替过程分析,在产量方程中引入启动压力梯度,结合油井定产生产和注入井定量注入的条件,构建出无限大的地层中考虑启动压力影响的一注一采渗流模型,通过数学推导得出注采井间速度场和压力场分布的数学描述.基于典型注采条件下渗流场的分布特点,提出实现有效驱替的条件,通过研究注采单元中渗流场分布规律及驱替压力梯度与有效驱替范围的关系,建立确定低渗透油田水驱开发有效驱替注采井距的计算方法.
不同形式的注采井网,由于其注采井数比、井距、排距和位置等方面的不同,其注采井间的压力场、速度场分布规律不尽相同; 但压力场是由各个单井(注入、采出)的势分布迭加、变换而成,而速度场是根据压力场分布进行相应的计算得到; 因此,一注一采条件下的压力场、速度场就是典型的注采井间渗流场,是研究不同注采形式井网下注采井间渗流场的基础,反映了渗流场分布的最基本变化规律.因此,基于一注一采条件下的渗流问题,在基本渗流方程中引入启动压力梯度,结合油井定量生产和注入井定量注入的条件,构建出无限大地层中考虑启动压力影响的一注一采渗流模型,作为定量分析注采井间渗流场分布及变化规律的基础.
研究的目标油藏类型是孔隙型的低渗透砂岩油藏,在注水开采的过程中,渗吸现象(渗析驱替)及其作用很小,一般可以忽略,因此,在相应的注采渗流方程中,都省略渗吸项.
考虑启动压力梯度,渗流速度[13-15]应满足
v={0, |(p)/(r)|≤G
-k/μ(p)/(r)(1+G/(p/r)), |(p)/(r)|>G(1)
其中, p为半径r处的地层压力; G为启动压力梯度; k为储层渗透率; μ为地层原油黏度.
设无限大地层中存在等产量的两口井, A点注水井产量为-q, B点生产井产量为q; 流体为单相不可压缩液体.见图1.
根据假设条件,则无限大地层等产量一注一采考虑启动压力时的基本微分方程[13-15]为
{Z1/r(d)/(dr)[r((dp)/(dr)-G)]=0
r1((dp)/(dr))|r</sub>1→r</sub>w=(qμ)/(2πkh)
r2((dp)/(dr))|r</sub>2→r</sub>w=(-qμ)/(2πkh)(2)
其边界条件为
{p(r)|r=r</sub>1=pwf
p(r)|r=r</sub>2=pwi(3)
其中,两井相距为2d; 井筒半径为rw; r1和r2分别为任意点M到注入井、生产井的距离; 生产井井底流压为pwf; 注水井注入压力为pwi; h为储层有效厚度.
将一注一采同时生产的问题,先分解为无限大地层一口生产井和一口注入井单独工作,然后再进行速度的“矢量”叠加.
无限大地层中一口生产井单独生产时,任一点流体的速度大小v生为
v生=q/(2π((x-d)2+y2)1/2)(4)
生产井引起在x方向上的分速度v生x为
v生x=q/(2π((x-d)2+y2)1/2)·(d-x)/(((x-d)2+y2)1/2)(5)
同理,无限大地层中一口注入井单独生产时,任一点流体的速度v注为
v注=q/(2π((x+d)2+y2)1/2)(6)
注入井流体在x方向上的分速度v注x为
v注x=q/(2π((x+d)2+y2)1/2)·(x+d)/(((x+d)2+y2)1/2)(7)
而式(7)中的井口产量q, 在考虑启动压力的情况下可表示为
q=Q/h·(πk(pwi-pwf-G·2d))/(μln(2d/rw))(8)
将式(5)和式(7)两个解叠加,可得一注一采情况下地层任一点在x方向上的合速度为
vx=(k(pwi-pwf-G·2d))/(2μln(2d/rw))·(x+d)/((x+d)2+y2)+
(k(pwi-pwf-G·2d))/(2μln(2d/rw))·(x-d)/((x-d)2+y2)(9)
令式(9)中y=0, 可得在注采主流线(x轴)上的速度分布表达式为
vx|y=0=(k(pwi-pwf-G·2d))/(2μln(2d/rw))·
(1/(x+d)+1/(d-x))(10)
令式(9)中x=0, 可得在注采中间线(y轴)上的速度分布vx为
vx|x=0=(k(pwi-pwf-G·2d))/(2μln(2d/rw))·((2d)/(d2+y2))(11)
结合式(1),对式(10)进行积分,可得主流线上压力分布的情况为
p|y=0=-(pwi-pwf-G·2d)/(2ln(2d/rw))·[ln(x+d)-
ln(d-x)]-Gx+C(12)
利用注采井间中心点(x=0,y=0)处压力大小来确定式(12)的常数C
p|x=0, y=0=-(pwi-pwf-G·2d)/(2ln(2d/rw))·[ln(d)-
ln(d)]+C=C=(pwi+pwf)/2(13)
由此可得主流线上压力分布表达式为
p|y=0=-(pwi-pwf-G·2d)/(2ln(2d/rw))·[ln(x+d)-
ln(d-x)]-Gx+(pwi+pwf)/2(14)
同理,结合式(1),对式(9)进行积分,可得一注一采条件下井间压力分布表达式为
p|y=0=-(pwi-pwf-G·2d)/(2ln(2d/rw))·
{1/2ln[(x+d)2+y2]-
1/2ln[(d-x)2+y2]}-Gx+
(pwi+pwf)/2(15)
结合式(4)和式(6),通过矢量相加, v=v生+v注,可得一注一采条件下井间速度大小v分布为
v=q/(π)d/(((x-d)2+y2)1/2·((x+d)2+y2)1/2)(16)
在考虑启动压力梯度的影响时,注采井间压力场、速度场的分布规律,与经典渗流理论的结果存在一定差别.当给定注采压差和启动压力梯度等条件,在注采井中间线(y轴)上存在一点M0(x=0,y=b0), 该点压力梯度等于启动压力梯度,则该点所处流线与对称流线(关于x轴对称),所包含区域内的压力梯度均大于启动压力梯度,其中的流体都能够流动,称为可驱动区域.在图1中,方框内区域为注采井控单元,阴影区域为其中的可流动区域.
基于以上分析,可建立考虑启动压力影响、在不同生产条件下确定油藏注采单元进行有效驱替规律的步骤如下:
① 已知储层启动压力梯度G和注采压差△p;
② 由式(10),在注采主流线中心处(x=0), 当速度值等于0时,可求得注采井间主流线上能够驱动(线驱动)的极限注采井距Dc;
③ 对于如图1的典型注采单元,图1中井控单元面积为SA=4bd; 在以下分析中,井控单元长度为2d, 宽度取值2b=2d.
① 在极限注采井距内,取不同的井距进行计算分析,即取值0<2d≤Dc;
② 由式(10),可先得出此时主流线中心处(x=0)的速度vm;
③ 在可动区域内,不同位置流动速度差别较大,为保证有效开发,要求可动边界上的速度大于0.由式(11),当其值等于θ×vm时(θ为外流动边界速度的倍数, 0<θ≤1), 可得到相应条件下注采单元中能够有效驱替的边界be; 则相应包含的区域,即为井控单元中的有效驱替区域.
① 当边界流速不为0时,其边界包含的面积即为注采井控单元中有效驱替面积.结合式(16),有效驱替面积可表达为
Se=∫x=dx=-d∫y=bey=-bev(x, y)dxdy=∫x=dx=-d∫y=bey=-beq/(π)·
d/(((x-d)2+y2)1/2((x+d)2+y2)1/2)dxdy(17)
② 注采井控单元有效驱替面积系数为
ωe=(Se)/(SA)(18)
一般0<ωe≤1. 有效驱替面积系数ωe越大,说明在井控单元中,得到有效驱替的范围越大.
对于给定的油藏条件和生产条件,可以作出不同注采压差、启动压力梯度下,可驱动面积系数、有效驱替面积系数与注采井距的关系,为实际油藏开发井网的部署,提供技术支持.
利用建立的方法进行相关理论研究.基本参数为:井半径为0.1 m,注采井距变化为±200 m; 启动压力梯度取值范围为0~0.05 MPa/m,注采压差△p取值范围为5~20 MPa.
基于式(15)和式(16),得到一注一采情况下,分别考虑有无启动压力梯度时的压力场、速度场分布,见图2和图3.
由图2和图3可见,由于启动压力梯度直接影响着流体的渗流过程,因此,有无启动压力梯度时的压力场分布存在一定的差别.在无启动压力梯度时(G=0)的压力场,与经典渗流理论结果相同, y轴线是等压线, x轴是流线,整个水动力学场关于y轴对称.在存在启动压力梯度时(G≠0)的压力场,在相同的注采压差下,由于井间的压力消耗增加,使得等压线的分布形态出现较大变化,等压线在注采井间的变化幅度大、变化频率快; 而速度场也关于x轴和y轴对称,但流速大于0的范围明显变小.
对于不同启动压力梯度和注采压差条件,首先计算出相应的极限注采井距,再确定有效驱替面积系数的变化规律.根据式(10)计算不同条件下极限注采井距,结果见图4.
由图4可见,随着启动压力梯度的增大,极限注采井距呈幂函数逐渐减小; 极限井距是注采井网部署时的最大限制井距,在相同的启动压力梯度条件下,增大注采压差可以在一定程度上增大极限注采井距.
图4 不同条件下极限注采井距变化曲线
Fig.4 Variation curves of limit injector producer distance under different conditions
取启动压力梯度G为0.01~0.05 MPa/m、注采压差△p为5~20 MPa,有效驱替面积系数与注采井距的变化关系,结果见图5.
由图5可见,在极限井距内,随着注采井距的减小,有效驱替面积系数逐渐增大; 注采压差越大、启动压力梯度越小,井控单元内有效流动区域的面积越大.
在实际开发中,希望在注采井控单元中能够有效流动的区域越大越好,一般有效驱替面积系数ωe>0.5, 此时其对应的开发井距可称为有效注采井距.可以得出不同启动压力梯度、不同注采压差下,有效注采井距的变化如图6.
图6 不同条件下有效注采井距变化曲线
Fig.6 Variation curves of effective injector producer distance under different conditions
由图6可见,随着启动压力梯度的增大,油藏有效注采井距逐渐减小; 而随着注采压差的增大,有效注采井距逐渐增大.
设系数λ为有效注采井距与极限注采井距的比值,其值大小与启动压力梯度和生产压差呈非线性关系.对比图4和图6,在相同启动压力梯度G下,系数λ随注采压差△p的增大而增大.
系数λ与启动压力梯度G和注采压差△p的变化关系,可表示为
λ=0.533 4+e(-2.108+11.72G+0.0125 7△p)(19)
在G为0.01~0.05 MPa/m、 △p为5~20 MPa时,有效注采井距与极限注采井距的比值约在0.7~0.8.由于极限注采井距的计算相对简单,在实际操作中可先求出油藏的极限注采井距,然后利用关系式(19)求出λ, 进而得到对应的有效注采井距,以此来指导低渗透油藏的开发部署.
综上研究可知:
1)考虑启动压力梯度后,使注采井间渗流阻力增加,导致注采井间的压力场、速度场分布规律发生明显变化.
2)启动压力梯度对水驱开发的影响,主要体现在增加渗流阻力、降低流体在介质中的流动能力,使得水驱有效动用的范围明显减小; 为满足注采单元有效驱替范围的要求,有效注采井距一般都小于其极限注采井距,两者的比值与启动压力梯度和生产压差呈非线性关系.
3)本研究基于典型注采条件下渗流场的分布特点,提出实现有效驱替的要求,进而建立了计算低渗透水驱开发有效注采井距的方法; 计算分析表明,本方法具有很好的操作性,可以指导此类油藏注采井距的合理部署.
深圳大学学报理工版
JOURNAL OF SHENZHEN UNIVERSITY SCIENCE AND ENGINEERING
(1984年创刊 双月刊)
主 管 深圳大学
主 办 深圳大学
编辑出版 深圳大学学报理工版编辑部
主 编 阮双琛
国内发行 深圳市邮电局
国外发行 中国国际图书贸易集团有限公司(北京399信箱)
地 址 北京东黄城根北街16号
邮 编 100717
电 话 0755-26732266
0755-26538306
Email journal@szu.edu.cn
标准刊号 ISSN 1000-2618
CN 44-1401/N