APP下载

渗流圆的等效原则

2018-12-18李传亮邓鹏朱苏阳刘东华

新疆石油地质 2018年6期
关键词:油井渗流水平井

李传亮,邓鹏,朱苏阳,刘东华

(1.西南石油大学 石油与天然气工程学院,成都 610599;2.中海石油(中国)有限公司 湛江分公司,广东 湛江 524057)

油气水在地层中的流动与电流在导体中的流动类似,描述渗流的达西定律与描述电流的欧姆定律也具有相同的形式,根据电学原理研究渗流力学问题,即水电相似原理[1]。等值渗流阻力法是渗流力学和油藏工程的常用方法[2-3],是根据水电相似原理建立起来的一种等效方法,是欧姆定律在渗流力学中的应用。由于该方法简单且直观,因此,在石油工程领域得到了广泛应用,以解决复杂的渗流问题[4-7]。等值渗流阻力法包括线性流和径向流两种基本流动形态的等效处理。线性流的等效处理较为简单,实际应用中基本没有分歧。但是,对于径向流的等效处理却一直存在不同观点,主要是对渗流圆的等效原则有不同的理解,等效后的渗流圆有大、小之分,因此,等效后的渗流阻力和油井产量公式便存在差别,计算的油井产量也不同。为更好地使用等值渗流阻力法解决渗流力学问题,必须对渗流圆的等效原则加以明确,只有统一了等效原则,才能得到相同的等效结果。

1 径向流与线性流

若一口直井位于圆形地层的中心(图1),如不考虑地层伤害,油井的产量(地层的流量)公式[8]为

图1 圆形地层直井渗流示意

径向流的所有流线,都沿着圆形地层的径向线指向圆心(图1),也就是一个渗流圆。把径向流的产量公式(1)式写成欧姆定律的形式,则为

由(1)式和(2)式,得径向渗流阻力:

若地层为一线性地层,流线为一组平行线(图2),此即线性流,地层的流量[8]为

图2 线性地层渗流示意

把(4)式写成欧姆定律的形式,则为

由(4)式和(5)式,得线性流的渗流阻力为

2 简单渗流区域形状

实际地层通常不是圆形或线性地层,而是不规则的形状,流线也往往不是平行线和径向线,更多时候是弯曲的曲线,这种情况下的产量公式很难给出。但是,前人想到了一种近似处理方法,即把不规则的渗流区域形状等效成圆形或线性地层,这样就可以方便地给出油井的产量公式。

图3为一个正方形地层,正方形的边长为a,渗流区域的面积为a2,直井位于正方形的中心,若把它等效成图1的渗流圆,有两种等效处理方法。

图3 正方形渗流区域

第一种等效处理方法就是面积等效,让正方形与圆形的面积相等,则等效渗流圆的半径为

第二种等效处理方法就是周长等效,让正方形与圆形的周长相等,则等效渗流圆的半径为

很显然,用周长等效法得到的渗流圆是大渗流圆,用面积等效法得到的渗流圆是小渗流圆,因此,计算的产量自然有所不同。

对于复杂的任意渗流区域形状,前人通过严格的解析方法,给出了一个统一的等效公式[9],即

对于正方形地层,CA=30.9,代入(9)式,得re=0.57a,与面积等效法得到的小渗流圆十分接近,即小渗流圆在理论上更合理一些,而大渗流圆则高估了渗流阻力,低估了油井产量。

3 复杂渗流区域形状

封闭边界的任意泄油区域形状的等效问题,通过(9)式基本上都能解决。但是,若地层存在部分开放边界或定压边界,(9)式就无能为力了,这时必须进行分区等效。

3.1 线性流—径向流组合

排状注水开发井网的复杂渗流,可以分解成多个渗流单元,每一个单元的渗流都相同,研究其中的一个单元即可。图4为排状注水开发井网的基本渗流单元,单元两侧的封闭边界为开发井网的分流线,两端的定压边界代表注水井排。由其流线分布可以看出,远离油井区域的流线近乎平行,可以用线性流进行等效;而近井区的流线向井底汇聚,可以用径向流进行等效。因此,整个流动可分解为两个线性流和一个径向流。

径向流的渗流圆可以采用两种方法加以确定,一种方法是让渗流圆的直径与地层的宽度相等(图5),即直径等效法,等效后的渗流圆为大渗流圆,渗流圆半径为

图4 排状注水开发井网基本渗流单元

图5 排状注水开发井网直径等效法等效渗流单元(大渗流圆)

渗流圆的渗流阻力为

线性流由左右两个完全相等的部分并联组成,其渗流阻力为

由于地层的长度远大于渗流圆的半径,因此(12)式忽略了渗流圆对线性流的影响。地层总的渗流阻力为

按照欧姆定律,油井的产量公式为

另外一种处理方法,是让渗流圆的周长与地层的宽度相等(图6),即周长等效法,等效后的渗流圆为小渗流圆,渗流圆半径为

图6 排状注水开发井网周长等效法等效渗流单元(小渗流圆)

油井的产量公式为

(14)式与(16)式的唯一不同就是分母上差了一个π,大渗流圆的渗流阻力偏大,油井产量就相应地偏小。实际工作中大都采用周长等效法,但是,哪一种更为合理,需要用后面的径向流—径向流组合加以证明。

3.2 径向流—径向流组合

图7 圆形地层环状井网

在圆形地层中布置一个环形井排,井排半径为rp,油井数为n(图7),环形井排的供给半径为re.求解环形井排的生产问题,与求解扇形地层等分线上一口直井的生产问题完全相同。

图7中扇形地层中的渗流包括两个径向流,一个是从供给边界到生产井排的径向流(外部径向流),另一个是从生产井排到油井的径向流(内部径向流)。外部径向流不需要等效处理,其渗流阻力称作外阻,外阻的大小为

内部径向流有两种等效处理方法,一种方法是让渗流圆的直径与扇形地层井排的弧长相等,即直径等效法,等效后的渗流圆为大渗流圆,渗流圆半径为

内部渗流圆的阻力称作内阻,内阻的大小为

渗流总阻力为

于是,按照欧姆定律得油井的产量公式为

另一种处理方法是让渗流圆的周长与扇形地层井排的弧长相等,即周长等效法,等效后的渗流圆为小渗流圆,渗流圆半径为

油井的产量公式为

(21)式与(23)式的唯一不同,就是分母上差了一个π.(23)式在渗流力学中利用保角变换进行了严格的理论证明[10],因此具有一定的理论合理性,而(21)式则没有经过理论上的证明。换言之,小渗流圆比大渗流圆更加合理。

4 水平井的情况

图8 圆形地层水平井渗流

若在圆形地层中心打一口水平井,则地层中的流动既不是径向流,也不是线性流,而具有复杂的流动形态(图8)。

图10 水平井近井区的径向流

求解水平井的渗流,可以把图8中的复杂渗流简化成远井区的径向流(外部渗流圆,图9)和近井区的径向流(内部渗流圆,图10)两部分,地层的渗流阻力也由外阻和内阻两部分组成。

关于图10中内部渗流圆的等效处理结果也有大渗流圆和小渗流圆之分。大渗流圆的等效方法是让渗流圆的直径等于地层的厚度(直径等效法),小渗流圆的等效方法是让渗流圆的周长等于地层的厚度(周长等效法)。这两种处理方法与图5和图6中直井的等效处理方法类似,只是把地层的宽度变成了图10中的地层厚度。

许多水平井的产量公式都采用了小渗流圆的等效处理方法,如Borisov公式、Giger-Resis-Jourdan公式、Renard-Dupuy公式。

Borisov公式[11]:

采用大渗流圆进行等效处理的水平井产量公式较少,主要有Joshi公式和陈元千公式。Joshi公式[14]:

(27)式与(26)式只在分母上相差一个π.

陈元千公式[15]:

5 计算分析

5.1 水平井

若地层渗透率为0.05 D,地层厚度为10 m,生产压差为1 MPa,流体黏度为1 mPa·s,地层泄油半径为300 m,水平段长度为300 m,油井完井半径为0.1 m.把有关参数代入Renard-Dupuy公式得油井的产量为183.34m3/d,代入Joshi公式得油井的产量为178.74m3/d.

很显然,Joshi公式低估了油井的产量,但二者的差别并不大,仅为2.57%.

实际上,Joshi公式与Renard-Dupuy公式只是在分母上相差一个π,而且π又位于对数符号内,加上水平井的水平段长度远大于地层厚度,因此,不同等效方法对水平井的产量计算结果影响甚微。因此,尽管Joshi公式缺少理论上的合理性,但却一直保留了下来。

5.2 直井

若地层长度为600 m,地层宽度为300 m,其他参数同前。把有关参数代入(14)式,得到的油井产量为25.96 m3/d;代入(16)式,得油井的产量为29.15 m3/d,二者相差12.7%.

由此可见,不同等效方法对直井产量的影响较大,采用大渗流圆的(14)式明显低估了油井的产量,实际工作中应避免使用。

6 结论

(1)等值渗流阻力法有两种渗流圆的等效方法,分别等效出了大渗流圆和小渗流圆。大渗流圆高估了渗流阻力,而低估了油井的产量。

(2)由于水平井长度对内阻的削弱作用,不同等效方法对水平井产量的计算结果影响较小。但是,不同等效方法对直井产量的计算结果影响却较大。

(3)由于大渗流圆的等效方法缺少理论上的合理性,而且又对直井产量的计算产生了明显的影响,建议今后废弃大渗流圆的等效方法,在应用等值渗流阻力法时尽量采用小渗流圆的等效方法。

符号注释

a——正方形边长,m;

A——正方形泄油区域面积,m2;

b——地层长度,m;

CA——Dietz形状因子,不同泄油区域形状取不同数值;

d——地层宽度,m;

h——储集层厚度,m;

k——储集层渗透率,D;

L——水平井的水平段长度,m;

n——环形井排上的井数,口;

q——油井产量(地下),m3/ks;

re——油井泄油半径,m;

rei——内部渗流圆的半径,m;

rp——环形井排半径,m;

rw——油井完井半径,m;

Rl——线性地层的渗流阻力,mPa·s/(D·m);

Rr——径向渗流阻力,mPa·s/(D·m);

Rre——外部径向流的渗流阻力,mPa·s/(D·m);

Rri——内部径向流的渗流阻力,mPa·s/(D·m);

Δp——生产压差,MPa;

γ——常数,取1.781…;

μ——流体黏度,mPa·s.

猜你喜欢

油井渗流水平井
油井遭袭
长河坝左岸地下厂房渗流场研究及防渗优化
考虑各向异性渗流的重力坝深层抗滑稳定分析
新型油井水泥消泡剂的研制
一种油井水泥用抗分散絮凝剂
浅论水平井压裂酸化技术的改造
陆梁油田油井H2S综合治理技术分析
一种计算水平井产能的新方法
水平井段超临界CO2携岩数值模拟
热采水平井加热半径计算新模型