Loading [MathJax]/jax/output/HTML-CSS/jax.js
  岩性油气藏  2024, Vol. 36 Issue (4): 178-188       PDF    
×
边水气藏水侵动态分析方法及水侵主控因素
秦正山1, 何勇明1, 丁洋洋1, 李柏宏2, 孙双双1    
1. 成都理工大学 能源学院, 成都 610000;
2. 中国石油长庆油田公司 第六采油厂, 西安 710000
摘要: 基于经典分流理论与气水渗流规律,推导了考虑气体高速非达西渗流的边水气藏分流方程,并根据四川盆地某边水气藏多口产水气井的基础资料,分析了水侵动态特征参数的变化规律及主控因素,探讨了延缓气藏水侵的技术思路与对策。研究结果表明:①新的改进分流方程能够确定不同开发时刻气藏(或气井)边水的整体推进情况及见边水时间,相较于传统的基于达西定律的分流方程,该计算结果更加可靠。②边水气藏水侵动态受制于多种因素的综合影响,其中储层渗透率对其影响最为显著,其次为非达西流系数、相对渗透率及孔隙度,而有效厚度、供气边界及水侵流量的影响程度较小。③充分发掘物性较均一的低渗致密段储层的开发潜力是提升边水气藏开发效果的关键,制定合理的气藏采气强度是主动控水、稳水的重要手段。
关键词: 分流理论    水侵动态    非达西流    熵权法    物性    有效厚度    供气边界    水侵流量    低渗致密储层    采气强度    边水气藏    
Water invasion performance and main controlling factors for edge-water gas reservoirs
QIN Zhengshan1, HE Yongming1, DING Yangyang1, LI Baihong2, SUN Shuangshuang1    
1. School of Energy, Chengdu University of Technology, Chengdu 610000, China;
2. No. 6 Oil Production Plant, PetroChina Changqing Oilfield Company, Xi'an 710000, China
Abstract: Based on classical fractional flow theory and gas-water percolation law, the fractional flow equation for edge-water gas reservoirs considering high-speed non-Darcy flow of gas phase was derived. The characteristic parameters and controlling factors of water invasion performance were analyzed based on the data of multiple water-producing gas wells in an edge-water gas reservoir in Sichuan Basin, and the technical approaches and strategies for delaying water invasion in gas reservoirs were discussed. The results show that: (1)The derived fractional flow equation can determine the overall advancement and breakthrough time of edge water in the reservoir(or gas well)at different development times, and the calculation results are more reliable than that by the traditional Darcy flow model.(2)The water invasion performance of the edge-water gas reservoirs is affected by a combination of multiple factors, the reservoir permeability has the most significant influence, followed by non-Darcy flow coefficient, relative permeability and porosity, while effective thickness, gas supply boundary and water invasion flow rate are less influential.(3)Fully exploiting the development potential of low-permeability and tight reservoirs with relatively uniform physical properties is the key to enhancing the development effectiveness of edge-water gas reservoirs, and formulating a reasonable gas recovery intensity is an important means to actively control and stabilize water.
Key words: fractional flow theory    water invasion performance    non-Darcy flow    entropy method    physical properties    effective thickness    gas supply boundary    water invasion flow rate    low-permeability and tight reservoir    gas recovery intensity    edge-water gas reservoir    
0 引言

对于水驱气藏而言,在其开发过程中随着气藏压力的下降,边(底)水逐渐推进至井底,导致气井见水,产水量迅速增加,产气量大幅下降,甚至部分气井出现水淹的状况,从而影响了气藏开发效果[1-2]。水驱气藏在我国气藏类型中占有相当大的比例(约占65%),正确认识水驱气藏的水侵动态变化规律及控制因素,是成功开发该类气藏的关键[3-4]。目前,采用水驱特征曲线[5-6]、物质平衡法[7-8]、经验公式[9]等方法研究水驱气藏各参数的动态大小及变化规律已有较多研究成果,但这些方法仅能笼统分析水侵量、累计产气量和累计产水量随时间的变化[10-13],难以体现复杂的地下渗流及驱替特征,也无法了解不同时刻地层水的推进情况及水侵动态的控制机制。一些研究人员[14-15]通过气藏数值模拟手段,虽然可以分析不同开发阶段的水线推进情况来了解气藏的水侵动态变化,但是存在复杂渗流机理认识不深入的缺点,分析结果存在较多模糊性。

基于分流理论,提出一种新的改进分流方程。首先推导考虑气体高速非达西渗流的分流方程,以四川盆地某边水气藏为例,分析典型因素对边水气藏动态变化的影响规律;其次,结合正交实验设计与熵权法,明确边水气藏动态变化的主控因素;然后,利用新的改进分流方程和传统的基于达西定律的分流方程,分别计算多口产水气井的理论见边水时间,并将计算结果分别与矿场统计值进行对比与分析;最后,基于研究区水侵主控因素的分析结果,探讨延缓边水气藏水侵的技术思路与对策,以期为水侵动态分析、治水稳产方案的编制提供一定参考。

1 边水气藏分流方程推导 1.1 经典分流理论的适用性

Buckley等[16]基于一维水驱油线性系统,提出了经典分流理论(图 1a),描述了油水两相非混相驱替过程,其中油水两相流体渗流过程满足达西定律。经典分流理论包含2个重要方程,即分流率方程和等饱和度推进方程。

下载原图 图 1 基本几何流动模型示意图 Fig. 1 Schematic diagram of basic geometric flow models

驱替相分流率方程[12]

fD=qDqT=qDqD+qo=vDvD+vo=11+KoμDKDμo (1)

式中:fD为驱替相分流率;qDqoqT分别为驱替相、被驱替相及两相总体积流量,m3/s;vDvo分别为驱替相与被驱替相渗流速度,m/s;KDKo分别为驱替相与被驱替相有效(相)渗透率,mD;μDμo分别为驱替相与被驱替相黏度,mPa·s。

等饱和度推进方程(或推进速度方程)[16]

Δu=qTφA(dfDdSD)Δθ (2)

(uθ)SD=qTφA(fDSD)θ (3)

式中:u为驱替相沿着渗流方向的推进距离,m;SD为驱替相饱和度,%;θ为时间,ks;A为渗流截面积,m2φ为孔隙度,%。

根据水体与天然气的分布可将水驱气藏分为边水气藏和底水气藏。边水气藏的流体性质和渗流规律与水驱油藏具有明显差异:①水驱油藏与边水气藏的高压物性参数不同,导致二者的驱替特征不同;②对于水驱油藏,渗流速度通常达不到紊流发生的速度界限,而对于边水气藏,当其整体采气强度或气井采气量提高时,气井容易在近井地带产生高速非达西效应,随之产生由流量变化引起的附加压降,流体的渗流规律偏离达西定律;③边水气藏径向系统下的基本几何流动模型与水驱油藏不同,径向系统下边水气藏流体的流动方向由供气边界流向井筒,基本几何流动模型为“点汇”模型(图 1b),而常规水驱油径向基本几何流动模型为“点源”模型,即流体的流动方向由注入端向外延伸(图 1c)。以上原因导致常规水驱油的分流方程难以应用于边水气藏。

1.2 改进的分流率方程

基于达西定律的水相运动方程[15-17](不考虑毛管力、重力)为

vw=kkrwμwdpwdr (4)

式中:vw为水相渗流速度,m/s;k为绝对渗透率,mD;krw为水相相对渗透率;pw为水相流体压力,MPa;μw为水相黏度,mPa·s;r为径向距离,m。

考虑到径向地层气井高速生产过程中容易产生高速非达西效应,根据Forchheimer[18]提出的二次方程,对于气相有

dpgdr=μgvgkkrg+βρgvg2 (5)

对Forchheimer二次方程求根,物理上方程仅存在唯一解,可将式(5)进一步表示为

vg=μgkkrg+[μ2gk2k2rg+4ρgβdpgdr]1/22ρgβ (6)

式中:vg为气相渗流速度,m/s;μg为气相黏度,mPa·s;krg为气相相对渗透率;ρg为气相密度,kg/m3β为非达西流系数,m-1pg为气相流体压力,MPa。

Evans等[19]的研究表明,如果利用单相条件下的非达西流系数与绝对渗透率、孔隙度构建的经验方程,来预测气水两相渗流系统下的非达西流系数,会导致预测结果与实验结果之间出现数量级的偏差。由此提出采用气相有效渗透率来替换绝对渗透率,采用气相饱和度来修正孔隙度,所获经验关系更能反映气水两相渗流系统下的渗流规律及由惯性效应主导的物理规律。Evans等[19]、LIU等[20]将两相渗流系统下的非达西流系数称为“有效非达西流系数”(本文用“βeff”表示),关系式为

βeff=Cβ(kkrg)5/4[φ(SgSgr)]]3/4 (7)

式中:βeff为有效非达西流系数,m-1Cβ为非达西流常数,m3/2,砂岩储层取值高于碳酸盐岩储层1~2个数量级[20-21],SI单位制下取值3.2×10-7~3.2×10-5 m3/2SgSgr分别为含气饱和度和残余气饱和度,%。

计算气体高速非达西渗流条件下的各相分流率时,首先需要获得各相渗流速度。假设不考虑毛管压力,气相与水相具有相同的压力梯度,式(5)可变形为

vg=kkrgμgdp dr(1+βeffρgkkrgvgμg)1 (8)

δ=βeffρgkkrgμg,χ=(1+δvg)1 (9)

联立式(4)、式(6)、式(7)—(9),考虑气相高速非达西渗流的分流率方程可写为

fhw=vwvw+vg=11+krgkrwμwμgχ (10)

fhg=vgvw+vg=11+krwkrgμgμwχ1 (11)

式中:fwhfgh分别表示考虑气体高速非达西渗流的水相与气相分流率。

1.3 饱和度推进方程、水侵前缘及见边水时间

图 2可看出,在边水侵入的初始时刻(t=0),地层总体积流量(Qt)等于水相体积流量(Qw);在边水逐渐侵入过程中,相邻微元体中的流体在某一时刻(t或某一时间段Δt内)满足质量守恒定律。

下载原图 图 2 边水气藏径向地层微元体示意图 re:供给边界;rw:井筒半径;Δr:微元体径向距离变化量;h:储层有效厚度。 Fig. 2 Schematic diagram of micro-unit body of radial formation in edge-water gas reservoirs

微元体满足如下质量守恒关系:

[(Qwρw)r+Δr(Qwρw)r]Δt=πhφ{[re(r+Δr)]2(rer)2}[(Swρw)t+Δt(Swρw)t] (12)

式中:Qw为水相体积流量,m3/s;Sw为含水饱和度,%;re为供给边界,m;r为某一特定径向位置或推进距离,m;Δr为微元体径向距离变化量,m;ρw为水相密度,kg/m3t为边水推进时间,s;Δt为边水推进Δr所需时间,s;h为储层有效厚度,m。

对式(12)化简可得

[(Qwρw)r+Δr(Qwρw)r]Δt=πhφ(Δr22reΔr+2rΔr)[(Swρw)t+Δt(Swρw)t] (13)

对于微元体,由于Δr → 0,Δt → 0,有下式满足近似相等关系:

Δr22reΔr+2rΔr2rΔr2reΔr (14)

将式(14)右项代入式(13)。假设水相密度不变,根据偏微分的定义,式(13)可写为

2πhφ(rwr)(Sw)t=(Qw)r (15)

根据分流率基本定义,式(15)可表示为

2πhφ(rer)(Sw)t=(Qtfhw)r (16)

式中:Qt为总体积流量,m3/s。

对于气井某一特定生产时刻,式(16)可进一步化简为

2πhφ(rer)Qt(Sw)t=fhwr (17)

由于水相分流率仅是含水饱和度的函数[15-16],对上述偏微分方程利用链式求导法展开,可得到边水气藏径向地层中两相流动的渗流微分方程(或驱替方程):

2πhφ(rer)QtSwt+fhwSwSwr=2πhφ(rer)QtSwt+fhwSwr=0 (18)

式中:fwh′为水相分流率导数。

含水饱和度Sw是时间t和位置r的二元函数,即Sw=frt)。选择地层流体等饱和度作为特征线,在特征线上,下式成立[15-16]

drdt=SwtSwr (19)

将式(19)带入式(18),整理得到径向地层任意水相等饱和度推进速度方程为

d[re2(rer)2]dt=fhwQtπhφ (20)

对式(20)关于时间t和径向距离r积分,得

re2(rer)2=Qtfhwπhφ(tt0) (21)

式中:t0为水侵初始时刻,s。

显然,式(21)仅存在唯一解满足物理规律,即边水气藏任意水相等饱和度面推进距离为

r=re[re2Q1fhwπhφ(tt0)]1/2 (22)

由于水侵流量在气藏开发过程中并非常数,平均水侵速度vt可采用数值积分方法计算,即采用物质平衡法获得不同阶段的水侵流量Qti(t),计算Qti(t)与dt围成的面积,并除以累计生产时间t,表达式为

¯vt=t0Qti(t)dtAt=t0Qti(t)dt2πreht (23)

式中:vt为平均水侵速度,m/s;Qti(t)为气井不同阶段生产阶段的水侵流量,m3/s。

联立式(22)、式(23),当t0=0时,有

r=re[r2et0Qti(t)dtπhφfhw]1/2=re[r2e2¯vtretφfhw]1/2 (24)

则井筒距水侵前缘的距离rw-f

rwf=[r2et0Qti(t)dtπhφfhw]1/2=[r2e2¯vtretφfhw]1/2 (25)

式中:rw-f为井筒距水侵前缘的距离,m。

通过以上推导,可得到边水气藏(或气井)径向地层任意水相等饱和度的线移动速度及推进方程[式(20)、式(22)或式(24)],也可得到井筒距水侵前缘的距离[式(25)]。

当气相区增加的水量等于边水侵入量时有:

r10(SwSwi)πhφd[r2e(rer)2]=Qt(tt0) (26)

式中:rf为水侵前缘推进距离,m。

整理式(26)为

r0(SwSwi)(2re2r)πhφdr=Qt(tt0) (27)

对式(21)求微分,可得

(2re2r)dr=Qt(tt0)πhφdfhw (28)

将式(28)带入式(27),可得

(SwSwi)fhw|rf0rf0 dfhw=1 (29)

式中:Swi为初始地层含水饱和度,%;fhwf,fhw分别为水侵前缘含水饱和度(Swf)对应的分流率及分流率导数。

进一步积分,得到边水气藏(或气井)水侵前缘饱和度方程为

fhw=fhwfSwfSwi (30)

式(30)的几何意义是由初始含水饱和度位置作水相分流率曲线的切线,切线的斜率值即fhw,切点位置对应的含水饱和度为水侵前缘饱和度。

水侵前缘到达井筒时,rf=re-rwrw2re2,由式(22)可得径向地层气井见水时间tbt

tbt=t0+π(re2rw2)hφQtfhwft0+πre2hφQtfhw (31)

式中:tbt为径向地层气井见水时间,s;rw为井筒半径,m。

2 水侵动态影响因素分析与水侵对策研究 2.1 水侵动态影响因素分析

边水气藏水侵动态变化的影响因素较多,选取部分典型因素(如非达西流常数、润湿性、水侵流量)分析其对水侵动态参数的影响规律。

当非达西流常数(Cβ)取值不同时,有效非达西流系数(βeff)差距较大。在相同的水侵流量下(Qt=0.40 m3/s),对Cβ取不同数值,计算川中L02-X1井在稳定水侵流量下的水相分流率与水侵前缘推进情况,并与达西流条件下的曲线进行对比(图 3图 4)。由图 3图 4可知,在水侵流量相同的情况下,非达西流常数Cβ越大,惯性效应越严重,气藏的“非活塞”驱替特征就越显著;水侵前缘饱和度越小,水侵速度越快,水侵前缘推进距离就越远。

下载原图 图 3 不同非达西流常数下L02-X1井水相分流率曲线 Fig. 3 Water fractional flow curves for well L02-X1 under different non-Darcy flow constants
下载原图 图 4 不同非达西常数下L02-X1井水侵前缘推进曲线(生产时间为3年) Fig. 4 Water invasion front edge curves for well L02-X1 under different non-Darcy flow constants

对于某一特定气藏,润湿性/饱和历史对气水相渗曲线也具有显著影响,进而影响气藏水侵动态[18-21]。计算不同润湿性(驱替与渗吸)条件下的分流率曲线与前缘推进曲线如5图 6所示。

下载原图 图 5 不同润湿性条件下L02-X1井水相分流率曲线 Fig. 5 Water fractional flow curves for well L02-X1 under different wettability
下载原图 图 6 不同润湿性条件下L02-X1井水侵前缘推进曲线 Fig. 6 Water invasion front edge curves for well L02-X1 under different wettability

在不同润湿性条件下,多相流体的运动机制会发生变化。当储层润湿相饱和度较高时,非润湿相逐渐取代润湿相,驱替作用占主导;当储层润湿相饱和度较低时,润湿相逐渐取代非润湿相,渗吸作用占主导。由图 5图 6可看出,不同润湿性条件下水相分流率曲线、水侵前缘推进曲线具有明显差异,对气藏的水侵动态具有一定影响。

以L01-1井为例,计算不同水侵流量下水相分流率及水侵前缘推进曲线如图 7图 8所示。

下载原图 图 7 不同水侵流量下L01-1井水相分流率曲线 Fig. 7 Water phase fractional flow curves for well L01-1 under different water invasionflow rates
下载原图 图 8 不同水侵流量下L01-1井水侵前缘推进曲线(生产时间为3年) Fig. 8 Water invasion front edge curves for well L01-1 under different water invasion flow rates

图 7图 8可看出,水侵流量越大,水侵前缘饱和度越小,边水推进速度就越快。当水侵流量达到2.5 m3/s时,边水早已推进至井筒,且该水侵流量下的水侵前缘饱和度(18.82%)仅略大于初始地层含水饱和度,反映了更为显著的“非活塞”驱替特征(图 8)。

2.2 水侵动态影响因素敏感性分析

边水气藏水侵动态的影响因素较多,通过边水气藏水侵动态因素敏感性分析,可进一步明确边水气藏动态变化的主控因素,探究延缓水侵的技术思路与对策。

2.2.1 水侵动态影响因素与评价指标的选择及量化

选择储层孔隙度、渗透率、水侵流量、有效厚度、供气半径、非达西流常数、相渗曲线特征参数等作为水侵动态的影响因素,将水侵前缘推进速度与见边水时间作为水侵特征的评价指标。由于相渗曲线特征参数较多,若将表征同一类因素的所有特征参数纳入,会因多因素分析而难以形成规律性认识,由此提出采用相渗曲线“共渗三角区面积占比”来表征各类储层渗流能力的大小,并作为反映相渗特征的唯一因素纳入量表(表 1)。“共渗三角区面积占比”为相对渗透率曲线与饱和度轴围成的双曲边三角形闭合面积占总相渗曲线区域面积的比例,可结合改进Corey幂律模型[19],利用定积分计算,表达式为

SΔ=SewSwkrw(Sw)d(Sw)+1SgrScwkrg(Sw)d(Sw) (32)

下载CSV 表 1 熵权法计算各因素影响权重结果 Table 1 Weight results of various factors calculated by entropy method

式中:SΔ为共渗三角区面积占比;Swc为共渗点水相饱和度,%;krwSw),krgSw)为水相与气相相对渗透率与含水饱和度函数。

表 1可以看出,储层渗透率在整个多元复杂系统中提供的信息最为关键,其对水侵动态的影响最为显著,其次为非达西流常数(或非达西流系数)、共渗三角区面积(或相对渗透率)、孔隙度,而有效厚度、供气边界、水侵流量的影响程度相同且较低。

水侵特征评价指标(见边水时间与水侵前缘推进速度)相较于各影响因素,其熵值低,在反映水侵动态的系统中具有较大信息量,在系统中权重系数总和高达49.18%,水侵特征评价指标与影响因素在整个系统中的地位同等重要,说明选用的评价指标可作为水侵变化的重要衡量指标。

依据岩石类型、孔隙特征及储层物性,可将研究区储层分为Ⅰ类(k ≥ 20 mD,φ≥12%)、Ⅱ类(0.1 mD ≤k<20 mD,6%≤φ<12%),Ⅲ类(0.01 mD≤k<0.1 mD,2%≤φ<6%))、Ⅳ类(k<0.01 mD,φ<2%)。由图 9可看出,从Ⅰ类储层到Ⅳ类储层,随着渗透率与孔隙度的减小,束缚水饱和度增大,气水两相区范围变窄,等渗点变低,等渗三角区面积变小,两相渗流能力变弱。

下载原图 图 9 气相、水相相对渗透率曲线“共渗三角区”及其储层分类特征示意图 Fig. 9 Schematic diagram of"co-permeable triangle area"of gas-water relative permeability curves and reservoir classification characteristics
2.2.2 正交设计

基于研究区静态、动态资料,确定影响因素水平,设计L32.4.7/9正交表(32组算例、4个水平、7个影响因素、2个指标),制定标准正交表。

各影响因素的水平取值为:渗透率(40.00 mD,30.00 mD,10.00 mD,0.05 mD);孔隙度(24%,18%,9%,4%);有效厚度(100 m,75 m,50 m,25 m);供气半径(2 000 m,1 500 m,1 000 m,500 m);Qt(1 m3/s,0.75 m3/s,0.5 m3/s,0.25 m3/s);非达西流常数(0.2×10-5 m3/2,3.2×10-5 m3/2,6.2×10-5 m3/2,9.2×10-5 m3/2)。对于砂岩储层,水相与气相的Corey拟合指数nwng均取值为1.5~4.0[19]。基于数值优化结果,制定相渗曲线综合表征参数(共渗三角区面积占比SΔ%)的水平为6.26%(ng=nw=1.5),3.75%(ng=nw=2.0),1.66%(ng=2.0,nw=4.0)和0.58%(ng=nw=4.0)。

2.3 延缓气藏水侵的技术思路与对策

储层渗透率、非达西流常数(或非达西流系数)、相对渗透率(共渗三角区面积)均是反映储层与流体渗流特性的重要影响因素,且各影响因素之间存在一定的相关性[15, 18-19, 22-24],而有效厚度、供气边界、水侵流量是反映储层与水体性质的宏观度量参数,并非地层流体渗流的决定性因素。

储层渗透率对边水气藏水侵动态的影响最为显著,渗透率发生变化时,非达西流常数与相对渗透率均会发生变化。通常,渗透率越大,边水推进速度越快,非均质性越强,边水非活塞驱替特征就越明显,采收率也就越低。因此,开展低渗致密有效储层识别与筛选,并结合水体分布及其性质认识,充分发掘物性分布相对均匀的低渗致密段储层的开发潜力是提升边水气藏开发效果的关键。

水侵流量尽管对水侵动态无直接相关关系,但是保持合理的储层压力可降低水侵流量,促进水线均匀推进,减少残余气局部富集,从而提高边水气藏最终采收率。因此,从生产角度,制定合理的气藏采气强度(或气井工作制度)是控制水侵速度及控水稳水的重要手段。

3 实例计算

四川盆地某边水气藏构造上为一近东西向的低幅度背斜,南北两翼稍陡,整体以薄层为主,平面连续性好,气藏温度为85 ℃,压力为34 MPa,储层孔隙度为4.0%~22.7%,渗透率为0.001~28.210 mD,非均质性强,气层有效厚度约30 m,天然气黏度为0.027 8 mPa·s,地层水黏度为0.357 1mPa·s,天然气密度为239.98 kg/m3,地层水密度为984.69 kg/m3,目前该气藏大量气井已见水,甚至部分气井出现水淹的状况。

利用新的改进分流方程对该气藏进行分析,首先统计该气藏10口气井的静态、动态资料,采用物质平衡法[2, 7, 10, 25],并结合数值积分法,计算各井的水侵流量;然后,利用Evans公式[15]与改进Corey模型[19],采用数值迭代法计算各气井的分流参数及各井的理论见边水时间,并将新的改进分流方程和传统分流方程的理论见边水时间分别与矿场统计的见边水时间进行对比,进而分析水侵动态指标的差异。

3.1 计算水线推进情况与见边水时间

以L02-X1井为例,该井投产后产量相对稳定,2019年7月产气量迅速下降,产水量迅速上升。井筒测试表明,该气井井筒积液逐渐增多,在已有的工作制度下气井难以携液生产。2021年3月气井停产,通过优化气井工作制度,开展排水采气、油管解堵等措施,气井产能得到一定程度的恢复。目前,该井油压由初期的5.8 MPa下降至3.2 MPa,下降了44.82%。该气井生产曲线与产出水矿化度变化如图 10图 11所示。

下载原图 图 10 L02-X1井生产曲线 Fig. 10 Gas production curves of well L02-X1
下载原图 图 11 L02-X1井产水前后矿化度对比 Fig. 11 Comparison of mineralization before and after water production in well L02-X1

边水气藏开发初期产水主要为凝析水与原生地层水,气井见边水后,产水量急剧增大,产气量急剧减小,L02-X1井符合典型的边水气藏气井生产动态特征[2, 4]图 10)。由于地层水中通常具有多种离子(如氯离子、钠离子、钙离子、锶离子等),其矿化度远大于凝析水中的矿化度。根据产出水矿化度分析资料(图 11)可知,随着气藏不断被开采,产出水矿化度具有明显的上升趋势,见边水时间与生产曲线产水量迅速增加时间基本吻合。利用新的改进分流方程,计算L02-X1井理论见边水时间为1 100 d,矿场统计见边水时间为1 125 d(表 2),二者差距较小。

下载CSV 表 2 气井、地层基础参数与见边水时间预测误差对比 Table 2 Basic parameters of gas wells and formation and prediction errors for edge water breakthrough time
3.2 见边水时间计算值与矿场统计值对比

理论见边水时间与矿场统计见边水时间的平均相对误差(Mean Relative Error,MRE%)的计算公式为

MRE(%)=100NiNii=1|tcalc treal |treal  (33)

式中:Ni为统计井数量,口;i为第i个统计数据;tcalctreal分别为计算的和矿场统计的见边水时间,d。

通过计算10口井的见边水时间(表 2)可知,新的改进分流方程理论见边水时间与矿场统计见边水时间的相对误差为-17.52%~15.82%,MRE%为9.77%。传统分流方程理论见边水时间与矿场统计见边水时间的相对误差为47.32%~163.62%,MRE%高达78.05%。因此,新提出的改进分流方程能较好地反映地层流体的渗流规律,但该分流方程对于常规储层的适应性更好,而对于非常规储层,由于存在特殊渗流机理,适应性仍需进一步研究。

4 结论

(1)通过扩展经典分流理论,推导了考虑气体高速非达西渗流的边水气藏分流方程,定量研究了不同开发时刻产水气井的边水推进状况及水侵动态参数变的化规律。新的改进分流方程相较于传统的基于达西定律的分流方程,计算结果更加可靠。

(2)渗透率对边水气藏水侵动态的影响最为显著,其次为非达西流常数、相对渗透率(共渗三角区面积)、孔隙度,而有效厚度、供气边界、水侵流量的影响程度均较低。

(3)充分发掘物性分布相对均匀的低渗致密段储层的开发潜力是提升边水气藏开发效果的关键,从生产角度,制定合理的气藏采气强度(或气井工作制度)是主动控水、稳水的重要手段。

参考文献
[1]
LEE J, WATTENBARGAR R A. Gas Reservoir Engineering. Houston: Society of Petroleum Engineers, 1996.
[2]
邓成刚, 李江涛, 柴小颖, 等. 涩北气田弱水驱气藏水侵早期识别方法. 岩性油气藏, 2020, 32(1): 128-134.
DENG Chenggang, LI Jiangtao, CHAI Xiaoying, et al. Early identification methods of water invasion in weak water drive gas reservoirs in Sebei gas field, Qaidam Basin. Lithologic Reservoirs, 2020, 32(1): 128-134. DOI:10.12108/yxyqc.20200114
[3]
杨宇, 孙晗森, 彭小东, 等. 气藏动态储量计算原理. 北京: 科学出版社, 2016.
YANG Yu, SUN Hansen, PENG Xiaodong, et al. Principles of gas reservoirs dynamic reserve calculation. Beijing: Science Press, 2016.
[4]
廖恒杰, 杨志兴, 李元生, 等. 西湖凹陷气藏出水规律及控制因素. 岩性油气藏, 2017, 29(6): 135-141.
LIAO Hengjie, YANG Zhixing, LI Yuansheng, et al. Water production laws and controlling factors of gas reservoir in Xihu Sag. Lithologic Reservoirs, 2017, 29(6): 135-141. DOI:10.3969/j.issn.1673-8926.2017.06.017
[5]
苟燕, 刘华勋, 高树生, 等. 水驱气藏的水驱特征曲线与应用效果分析. 石油钻采工艺, 2013, 35(3): 63-65.
GOU Yan, LIU Huaxun, GAO Shusheng, et al. Research and application effect analysis on water-drive characteristic curve for water-drive gas reservoir. Oil Drilling & Production Technology, 2013, 35(3): 63-65.
[6]
胡俊坤, 李晓平, 宋代诗雨. 水驱气藏动态储量计算新方法. 天然气地球科学, 2013, 24(3): 628-632.
HU Junkun, LI Xiaoping, SONG Daishiyu. A new method for dynamic reserves evaluation of water-drive gas reservoir. Natural Gas Geoscience, 2013, 24(3): 628-632.
[7]
陈恒, 杜建芬, 郭平, 等. 裂缝型凝析气藏的动态储量和水侵量计算研究. 岩性油气藏, 2012, 24(1): 117-120.
CHEN Heng, DU Jianfen, GUO Ping, et al. Study on calculation of dynamic reserves and water influx in fractured condensate gas reservoir. Lithologic Reservoirs, 2012, 24(1): 117-120. DOI:10.3969/j.issn.1673-8926.2012.01.023
[8]
谭先红, 梁斌, 王帅, 等. 一种低渗储层凝析气藏气井产能评价方法研究. 油气藏评价与开发, 2021, 11(5): 724-729.
TAN Xianhong, LIANG Bin, WANG Shuai, et al. A productivity evaluation method of gas wells in condensate gas reservoirs with low permeability. Reservoir Evaluation and Development, 2021, 11(5): 724-729.
[9]
ZHANG Jiqun, DENG Baorong, HU Changjun, et al. Computation method for water influx in different layers of natural edge water. Petroleum Exploration and Development, 2016, 43(5): 825-831. DOI:10.1016/S1876-3804(16)30098-2
[10]
邓成刚, 孙勇, 曹继华, 等. 柴达木盆地涩北气田地质储量和水侵量计算. 岩性油气藏, 2012, 24(2): 98-101.
DENG Chenggang, SUN Yong, CAO Jihua, et al. Calculation of original gas in place and water influx in the Sebei gas field, Qaidam Basin. Lithologic Reservoirs, 2012, 24(2): 98-101. DOI:10.3969/j.issn.1673-8926.2012.02.019
[11]
蒋琪, 王中元, 宫宇宁, 等. 顶水水侵对SAGD开采效率的影响及顶水下窜速率预测. 特种油气藏, 2022, 29(2): 72-76.
JIANG Qi, WANG Zhongyuan, GONG Yuning, et al. Influence of top water intrusion on SAGD production efficiency and prediction of downward channeling rate of top water. Special Oil & Gas Reservoirs, 2022, 29(2): 72-76.
[12]
何云峰, 杨小腾. 活跃边水气藏水侵系数与稳产期关系研究. 油气藏评价与开发, 2021, 11(1): 124-128.
HE Yunfeng, YANG Xiaoteng. Relation between water invasion coefficient and stable production period in gas reservoirs with active edge water. Reservoir Evaluation and Development, 2021, 11(1): 124-128.
[13]
张冰岩, 陈小凡, 乐平. 水侵缝洞型碳酸盐岩底水油藏弹性驱动单元开采研究. 油气藏评价与开发, 2020, 10(2): 71-75.
ZHANG Bingyan, CHEN Xiaofan, YUE Ping. Research on unit mining by elastic drive of fractured-vuggy carbonate reservoir with bottom water by water intrusion. Reservoir Evaluation and Development, 2020, 10(2): 71-75.
[14]
AMIRSARDARI M, ALALI N, AFSARI K. Modeling aquifer flow behavior in low-dip edge-water drive gas reservoirs. Natural Resources Research, 2021, 30(1): 479-493. DOI:10.1007/s11053-020-09741-0
[15]
李传亮. 油藏工程原理. 北京: 石油工业出版社, 2017.
LI Chuanliang. Principles of oil reservoir engineering. Beijing: Petroleum Industry Press, 2017.
[16]
BUCKLEY S E, LEVERETT M C. Mechanism of fluid displacement in sands. Transactions of the AIME, 1942, 146(1): 107-116. DOI:10.2118/942107-G
[17]
DAKE L P. Fundamentals of reservoir engineering. Amsterdam: Elsevier BV, 1983.
[18]
FORCHHEIMER P H. Wasserbewegung durch boden. Zeitschrift des Vereines Deutscher Ingenieure, 1901, 45(50): 1782-1788.
[19]
EVANS E V, EVANS R D. The influence of an immobile or mobile saturation upon Non-Darcy compressible flow of real gases in propped fractures. Journal of Petroleum Technology, 1988, 40(10): 1343-1351. DOI:10.2118/15066-PA
[20]
LIU X, CIVAN F, EVANS R D. Correlation of the non-Darcy flow coefficient. Journal of Canadian Petroleum Technology, 1995, 34(10): 50-64.
[21]
HOCKING R R. Methods and applications of linear models: Regression and the analysis of variance. Ishpeming: John Wiley & Sons, 2013.
[22]
GAO Chenglu, LI Shucai, WANG Jing, et al. The risk assessment of tunnels based on grey correlation and entropy weight method. Geotechnical and Geological Engineering, 2018, 36(3): 1621-1631. DOI:10.1007/s10706-017-0415-5
[23]
GHANBARNEZHAD-MOGHANLOO R, LAKE L W. Applying fractional-flow theory under the loss of miscibility. SPE Journal, 2012, 17(3): 661-670. DOI:10.2118/129966-PA
[24]
BROOKS R H. Hydraulic properties of porous media. Fort Collins: Colorado State University, 1964.
[25]
陈军, 秦柯, 任洪伟, 等. 利用气藏生产指示曲线计算凝析气藏水侵量. 岩性油气藏, 2015, 27(2): 103-108.
CHEN Jun, QIN Ke, REN Hongwei, et al. Estimation of water influx in condensate gas pool by means of productivity index curve. Lithologic Reservoirs, 2015, 27(2): 103-108. DOI:10.3969/j.issn.1673-8926.2015.02.016
边水气藏水侵动态分析方法及水侵主控因素
秦正山, 何勇明, 丁洋洋, 李柏宏, ...