CN105422070B - 优化特低渗透非均质气藏压裂水平井裂缝位置的方法 - Google Patents
优化特低渗透非均质气藏压裂水平井裂缝位置的方法 Download PDFInfo
- Publication number
- CN105422070B CN105422070B CN201510889494.8A CN201510889494A CN105422070B CN 105422070 B CN105422070 B CN 105422070B CN 201510889494 A CN201510889494 A CN 201510889494A CN 105422070 B CN105422070 B CN 105422070B
- Authority
- CN
- China
- Prior art keywords
- unit
- flow
- pressure
- crack
- permeable
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000035699 permeability Effects 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 title claims abstract description 35
- 239000012530 fluid Substances 0.000 claims abstract description 24
- 238000004088 simulation Methods 0.000 claims abstract description 13
- 230000008878 coupling Effects 0.000 claims abstract description 4
- 238000010168 coupling process Methods 0.000 claims abstract description 4
- 238000005859 coupling reaction Methods 0.000 claims abstract description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 36
- 238000004519 manufacturing process Methods 0.000 claims description 20
- 230000008569 process Effects 0.000 claims description 8
- 239000000758 substrate Substances 0.000 claims description 6
- 230000008595 infiltration Effects 0.000 claims description 5
- 238000001764 infiltration Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 238000009738 saturating Methods 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 6
- 239000007789 gas Substances 0.000 description 58
- 238000010586 diagram Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000005336 cracking Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000006399 behavior Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007717 exclusion Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Separation Using Semi-Permeable Membranes (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本申请公开了一种优化特低渗透非均质气藏压裂水平井裂缝位置的方法,主要包括以下步骤:收集储层、流体、水平井筒的基本参数;将非均质气藏沿水平井筒的长度方向划分成至少两个渗流带,渗流带含有人工裂缝;建立每一个渗流带的渗流单元模拟模型;建立至少两个渗流单元的渗流模拟模型;耦合各个渗流单元的渗流模拟模型,计算非均质气藏压裂水平井的产量;对比所述产量,选择特低渗透非均质气藏压裂水平产量最高的裂缝位置。本申请的技术方案充分考虑了储层非均质性和启动压力梯度的影响,从而克服了现有技术中只能实现对均质气藏、不考虑启动压力的压裂水平井裂缝位置选择,从而提高特低渗透非均质气藏水平井压裂改造的有效性和效果。
Description
技术领域
本申请属于油气田开发领域,具体地说,涉及一种优化特低渗透非均质气藏压裂水平井裂缝位置的方法。
背景技术
特低渗透气藏突出表现为储层渗透率平面非均质性和存在启动压力梯度的特点,通过对非均质气藏水平井分段压裂可以显著提高产量。在对特低渗透气藏压裂水平井裂缝位置优选时,必须综合考虑渗透率非均质性和启动压力梯度等因素的影响。目前针对压裂水平井产量预测主要有电模拟、保角变换、叠加原理、复变函数理论等方法。
上述方法均没有同时考虑特低渗透气藏非均质性以及特低渗透气藏流体流动时存在启动压力梯度的综合影响。
发明内容
有鉴于此,本申请所要解决的技术问题是没有同时考虑特低渗透气藏非均质性以及特低渗透气藏流体流动时存在启动压力梯度的综合影响。
为了解决上述技术问题,本申请公开了一种优化特低渗透非均质气藏压裂水平井裂缝位置的方法,主要包括以下步骤:
1)收集储层、流体、水平井筒的基本参数;
2)将非均质气藏沿是所述水平井筒的长度方向划分成至少两个渗流带,所述渗流带含有人工裂缝;
3)建立每一个所述渗流带的渗流单元模拟模型;
4)建立至少两个所述渗流单元的渗流模拟模型;
5)耦合各个所述渗流单元的渗流的模拟模型,计算非均质气藏压裂水平井的产量;
6)根据产量,优化特低渗透非均质气藏压裂水平井的裂缝位置。
进一步的,所述步骤1)中,所述储层的基本参数包括:原地应力方向、储层厚度、储层温度、孔隙度、渗透率;所述流体的基本参数包括:气体粘度、气体临界压力、气体偏差因子、气体临界温度;所述水平井筒的基本参数包括:水平井方位、水平井井筒长度。
进一步的,所述步骤2)具体为:根据储层在水平井筒长度方向上渗透率的差异,将渗透率相同的储层划分成同一个含有人工裂缝渗透带的渗流带。
进一步的,所述渗流单元的基本流动过程包括:基质线性流、当量井径向流和裂缝内的线性渗流;以第1个条带基质中的第1个渗流单元为例计算渗流单元中各流动过程的阻力。
进一步的,所述步骤3)具体包括:
a、根据关系式:计算Ⅰ区线性流的阻力,式中:Ruij为
第i渗透带中第j个渗流单元Ⅰ区基质线性流阻力,单位为MPa2/((m3/d)·(mPa·s));T为气
藏温度,单位为K;Tsc为标准大气温度,单位为K;psc为标准大气压,单位为MPa;lij为第i渗透
带中第j条裂缝与上一条裂缝渗流距离,单位为m;xf为裂缝半长,单位为m;ki为第i渗透带基
质渗透率,单位为mD;h为储层厚度,单位为m;
b、根据关系式:计算Ⅱ区径向流阻力;式
中:Rnij为第i渗透带中第j个渗流单元当量井径向流阻力,单位为MPa2/((m3/d)·(mPa·
s));kfij为裂缝渗透率,单位为Dc;wfij为裂缝宽度,单位为cm;
c、根据所述的Ⅰ区线性流的阻力和Ⅱ区径向流阻力,建立所述渗流带的渗流单元的模拟模型。
进一步的,所述步骤a中,所述关系式中,Ruij=ψuij-
ψdij,
根据气藏启动压力梯度下的基质线性流动流量与拟压力关系,得出公式:
根据拟压力定义,气体流过距离d时消耗的启动压力为Gij·d,将Gij·d定义为启动压力pB;因此可做以下转化:
式中:Ψuij、Ψdij分别为Ⅰ区基质线性流外拟压力和内拟压力,ΨBuij-ΨBdij为流体流过长度为dij Ⅰ区基质的启动拟压差,单位为MPa2/(mPa·s);Gij为Ⅰ区基质的启动压力梯度,单位为MPa/m;dij为Ⅰ区基质线性流长度,单位为m;qij为第i渗透带中第j条裂缝流量,单位为m3/d。
进一步的,所述步骤b中:
根据同一个渗流单元内,I区和II区边界位置处的压力相等,Ⅰ区线性流的内压即为Ⅱ区径向流的外压,则径向流流量与压差关系式为:
式中:Ψwfij为Ⅱ区径向流内拟压力,单位为MPa2/(mPa·s);rwij为当量井半径,单位为m。
所述当量直井的半径rw通过当量井径模型求取,具体为:
通过建立裂缝产量和直井产量相等的关系式从而得到当量井径与裂缝半长之间关系式,下面分别建立求取当量直径和裂缝的产量:
直井产量:求取定压供给边界条带气藏中一口直井A1的稳态产量,根据镜像反映原理,将其转化为一排生产井和一排注气井的产量稳态进行求解;
根据等值渗流阻力法,得一口井产量:
有限导流裂缝产量:
式中:Ψ为裂缝方向上x点处拟压力,单位为MPa2/(mPa·s)。
解出裂缝内平均压力近似为裂缝内压力,流体从基质流到裂缝内满足线性流动规律,线性流外拟压力为Ψuij,求出基质线性流流量,再由质量守恒得裂缝产量:
对比当量直井的产量式(6)与有限导流裂缝的产量式(8)得当量井半径公式:
式中,rw为当量井半径,单位为m;
将式(8)代入式(5)可以得到渗流单元Ⅱ区当量井附近等效径向流阻力为:
进一步的,所述步骤4)具体为:
在求得各区流动阻力之后,根据不同渗流区域间满足压力连续、流量相等的基本原理,利用等值渗流阻力法即求得压裂水平井产能公式。
进一步的,所述步骤5)具体为:流体由两边向中间流动,先分析两边渗透带裂缝产能方程再分析中间裂缝产能方程,相邻裂缝间用等值渗流阻力法连接各阻力得到个方程组成的线性方程组;
计算第1渗透带中裂缝产量,公式为:
计算第3渗透带中裂缝产量,公式为:
由于第2渗透带最外两个渗流单元分别与1、3渗透带相连,则其Ⅰ区基质线性流阻力为:
假设中间分流裂缝为第2渗透带中第m条裂缝所等效的井排,该井排为左右两流量汇合处,所以可列出N2+1个方程:
式中:q2m为左边流向分流井排的流量,单位为m3/d;q2(m+1)为右边流向分流井排的流量,单位为m3/d;分流裂缝的实际流量为两者之和,单位为m3/d;
根据拟压力定义,拟启动拟压差可采用压力平方差表示:
同理可得任意一点压力p的拟压力函数表达式:
根据气藏的对称性,压裂水平井总的产量为:
式中,Q为压裂水平井总产量,单位为m3/d。
进一步的,所述步骤6)具体包括:
a、设置非均质储层不同的裂缝位置方案;
b、对比计算不同裂缝位置方案的产量;
c、以压裂后的产量最优选选择裂缝位置。
与现有技术相比,本申请可以获得包括以下技术效果:
充分考虑了储层非均质性和启动压力梯度的影响,从而克服了现有技术中只能实现对均质气藏、不考虑启动压力的压裂水平井裂缝位置选择,从而提高特低渗透非均质气藏水平井压裂改造的有效性和效果。
当然,实施本申请的任一产品必不一定需要同时达到以上所述的所有技术效果。
附图说明
此处所说明的附图用来提供对本申请的进一步理解,构成本申请的一部分,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1是本申请实施例的非均质储层压裂水平井示意图;
图2A是本申请实施例的渗流单元实际流动的等效示意图;
图2B是本申请实施例的渗流单元等效流动的等效示意图;
图3A是本申请实施例的定压封闭边界镜像反映示意图;
图3B是本申请实施例的定压封闭边界镜像反映后的电路图;
图4是本申请实施例的气藏压裂水平井等效电路图;
图5是本申请实施例的不同渗透率组合下水平井产量曲线图;
图6是本申请实施例的不同裂缝位置布置示意图;
图1中1表示定压边界、2表示封闭边界、3表示水平井筒、4表示裂缝;图6中3表示水平井筒、4表示裂缝。
具体实施方式
以下将配合附图及实施例来详细说明本申请的实施方式,藉此对本申请如何应用技术手段来解决技术问题并达成技术功效的实现过程能充分理解并据以实施。
本申请公开了一种优化特低渗透非均质气藏压裂水平井裂缝位置的方法,主要包括以下步骤:
1)收集储层、流体、水平井筒的基本参数;
2)将非均质气藏沿是所述水平井筒的长度方向划分成至少两个渗流带,所述渗流带含有人工裂缝;
3)建立每一个所述渗流带的渗流单元模拟模型;
4)建立至少两个所述渗流单元的渗流模拟模型;
5)耦合各个所述渗流单元的渗流的模拟模型,计算非均质气藏压裂水平井的产量;
6)根据产量,优化特低渗透非均质气藏压裂水平井的裂缝位置。
进一步的,所述步骤1)中,所述储层的基本参数包括:原地应力方向、储层厚度、储层温度、孔隙度、渗透率;所述流体的基本参数包括:气体粘度、气体临界压力、气体偏差因子、气体临界温度;所述水平井筒的基本参数包括:水平井方位、水平井井筒长度。
如图1所示,所述步骤2)具体为:根据储层在水平井筒长度方向上渗透率的差异,将渗透率相同的储层划分成同一个含有人工裂缝渗透带的渗流带。
如图2A和图2B所示,压裂水平井的渗流可以看作是由很多相似流动过程的单元组成的,简称为渗流单元,每个渗流单元含有基质线性流、当量井径向流和裂缝内的线性渗流三个基本流动过程,本实施例以第1个条带基质中的第1个渗流单元为例计算渗流单元中各流动过程的阻力;
所述步骤3)具体包括:
a、计算Ⅰ区线性流的阻力:
考虑气藏启动压力梯度下的基质线性流动流量与拟压力关系式为:
上述(1)式中积分项难以确定,可采用拟压力定义。气体流过距离d时消耗的启动压力为Gij·d,将Gij·d定义为启动压力PB。因此可做以下转化:
式中:Ψuij、Ψdij分别为Ⅰ区基质线性流外拟压力和内拟压力,ΨBuij-ΨBuij为流体流过长度为dij Ⅰ区基质的启动拟压差,单位为MPa2/(mPa·s);Gij为Ⅰ区基质的启动压力梯度,单位为MPa/m;dij为Ⅰ区基质线性流长度,单位为m;qij为流量,单位为m3/d;xf为裂缝半长,单位为m;ki为第i渗透带基质渗透率,单位为mD;h为储层厚度,单位为m;T为气藏温度,单位为K;Tsc为标准大气温度,单位为K;psc为标准大气压,单位为MPa。
由于当量井附近存在等效径向泄气区域,则Ⅰ区线性流长度dij小于裂缝间距lij,径向流泄气周长可近似为裂缝半长,泄气面积近似为裂缝的直线源到裂缝间的矩形面积,则Ⅰ区线性流长度为:
所以Ⅰ区线性流的阻力为:
式中:Ruij为第i渗透带中第j个渗流单元Ⅰ区基质线性流阻力,单位为MPa2/((m3/d)·(mPa·s));
b、计算Ⅱ区径向流阻力:
根据同一个渗流单元内,I区和II区边界位置处的压力相等,Ⅰ区线性流的内压即为Ⅱ区径向流的外压,则径向流流量与压差关系式为:
式中:Ψwfij为Ⅱ区径向流内拟压力,单位为MPa2/(mPa·s);rwij为当量井半径,单位为m;
利用当量井径模型求取当量直井的半径rw,其基本思想就是通过建立裂缝产量和直井产量相等的关系式从而得到当量井径与裂缝半长之间关系式。下面分别建立求取当量直径和裂缝的产量。
如图3A和图3B所示,直井产量:为了求取定压供给边界条带气藏中一口直井A1的稳态产量,根据镜像反映原理,可将其转化为一排生产井和一排注气井的产量稳态进行求解;
利用等值渗流阻力法求得一口井产量:
有限导流裂缝产量:采用下述常微分方程描述有限导流裂缝的产量计算:
式中:式中:Ψ为裂缝方向上x点处拟压力,单位为MPa2/(mPa·s);kfij为裂缝渗透率,单位为Dc;wfij为裂缝宽度,单位为cm。
求解出裂缝内平均压力近似为裂缝内压力,考虑流体从基质流到裂缝内满足线性流动规律,线性流外拟压力为Ψuij,求出基质线性流流量,再由质量守恒得裂缝产量:
对比当量直井的产量式(6)与有限导流裂缝的产量式(8)得当量井半径公式:
将式(8)代入式(5)可以得到渗流单元Ⅱ区当量井附近等效径向流阻力为:
式中:Rnij为第i渗透带中第j个渗流单元当量井径向流阻力,单位为MPa2/((m3/d)·(mPa·s))。
进一步的,所述步骤4)具体包括:
计算物理模型:
如图4所示,在求得各区流动阻力之后,根据不同渗流区域间满足压力连续、流量相等的基本原理,利用等值渗流阻力法即求得压裂水平井产能公式;这里与气藏有三个不同的渗透率区域分布为例,第i渗透带压开的Ni条裂缝,忽略水平井筒压降的影响,由不同渗流单元的组合分布特征得到气藏等效电路图;
流体由两边向中间流动,先分析两边渗透带裂缝产能方程再分析中间裂缝产能方程,相邻裂缝间用等值渗流阻力法连接各阻力得到个方程组成的线性方程组。
进一步的,所述步骤4)具体包括:
计算产量:
计算第1渗透带中裂缝产量,公式为:
计算第3渗透带中裂缝产量,公式为:
由于第2渗透带最外两个渗流单元分别与1、3渗透带相连,则其Ⅰ区基质线性流阻力为:
假设中间分流裂缝为第2渗透带中第m条裂缝所等效的井排,该井排为左右两流量汇合处,所以可列出N2+1个方程:
式中:q2m为左边流向分流井排的流量,单位为m3/d;q2(m+1)为右边流向分流井排的流量,单位为m3/d,分流裂缝的实际流量为两者之和,单位为m3/d;;
根据拟压力定义,拟启动拟压差可采用压力平方差表示:
同理可得任意一点压力p的拟压力函数表达式:
根据气藏的对称性,压裂水平井总的产量为:
式中,Q为压裂水平井总产量,单位为m3/d。
进一步的,所述步骤6)具体包括:
a、设置非均质储层不同的裂缝位置方案;
b、对比计算不同裂缝位置方案的产量;
c、以压裂后的产量最优选选择裂缝位置。
实施例
本实施例应用优化特低渗透非均质气藏压裂水平井裂缝位置的方法,具体如下:
某非均质气藏水平井基础参数如下:非均质气藏长为800m,宽为120m,厚度为15m,两边边界压力为25MPa,水平井筒内压力22MPa,两边渗透带(第1和3渗透带)长度均为300m、基质渗透率均为0.6mD、均压裂2条均匀分布的垂直裂缝,中间渗透带(第2渗透带)长为200m、基质渗透率为1.2mD、压开1条垂直裂缝,所有裂缝半长为60m、裂缝导流能力为30D·cm,天然气粘度为0.0215mPa·s,偏差系数为0.9218。
表1为三个渗透带不同基质渗透率组合,为了对比均质气藏和非均质气藏压裂水平井产量计算结果的差异,同时考虑组合II、组合III另外两种储层渗透率参数的组合,其它参数同上。
表1
根据基质启动压力梯度和渗透率的关系:
式中:A为气藏常数;n取值为1。
如图5所示,现将A依次取值为0、0.05、0.1、0.15、0.2,得到与不同渗透带基质渗透率组合相应的启动压力梯度的组合,在不同气藏常数A下压裂水平井的日产量。
由图5可知:随着启动压力梯度增大(A值增大),同一组合下压裂水平井产量逐渐减少;不同组合减小趋势不一,且两边渗透率低、中间渗透率高的组合(Ⅰ)的下降趋势大于两边渗透率高、中间渗透率低的组合(III)压裂水平井产量;在同一启动压力梯度下,两边渗透率低、中间渗透率高的气藏(Ⅰ)压裂水平井产量低于两边渗透率高、中间渗透率低的气藏(III)压裂水平井产量,因此对于非均质气藏类型Ⅰ,若采用基质渗透率加权平均的方法处理为均质气藏II,则预测产能明显偏高;对于非均质气藏类型III,若采用基质渗透率加权平均的方法处理为均质气藏II,则预测产能明显偏低。
如图6所示,进一步的对比了该气藏压裂水平井3形成5条裂缝4,在不同裂缝4位置布置方案下的压裂水平井3产量。
表2为不同裂缝位置分布下的压裂水平井产量,注:d1-2表示第1条和第2条裂缝之间的距离;d2-3表示第2条和第3条裂缝之间的距离;d3-4表示第3条和第4条裂缝之间的距离;d4-5表示第4条和第5条裂缝之间的距离。由表2可知,采用方案I布缝的压裂水平井产量最高,因此优选方案I进行布缝。
表2
本申请充分考虑了储层非均质性和启动压力梯度的影响,从而克服了现有技术中只能实现对均质气藏、不考虑启动压力的压裂水平井裂缝位置选择,从而提高特低渗透非均质气藏水平井压裂改造的有效性和效果。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的商品或者***不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种商品或者***所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的商品或者***中还存在另外的相同要素。
上述说明示出并描述了本申请的若干优选实施例,但如前所述,应当理解本申请并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述申请构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本申请的精神和范围,则都应在本申请所附权利要求的保护范围内。
Claims (7)
1.一种优化特低渗透非均质气藏压裂水平井裂缝位置的方法,其特征在于,主要包括以下步骤:
1)收集储层、流体、水平井筒的基本参数;
2)将非均质气藏沿所述水平井筒的长度方向划分成至少两个渗流带,所述渗流带含有人工裂缝;
3)建立每一个所述渗流带的渗流单元的渗流模拟模型;
4)建立至少两个所述渗流单元的渗流模拟模型;
5)耦合各个所述渗流单元的渗流模拟模型,计算非均质气藏压裂水平井的产量;
6)根据产量,优化特低渗透非均质气藏压裂水平井的裂缝位置;
所述步骤3)具体包括:
a、根据关系式:计算Ⅰ区线性流的阻力,式中:Ruij为第i渗透带中第j个渗流单元Ⅰ区基质线性流阻力,单位为MPa2/((m3/d)·(mPa·s));T为气藏温度,单位为K;Tsc为标准大气温度,单位为K;psc为标准大气压,单位为MPa;lij为第i渗透带中第j条裂缝与上一条裂缝渗流距离,单位为m;xf为裂缝半长,单位为m;ki为第i渗透带基质渗透率,单位为mD;h为储层厚度,单位为m;
b、根据关系式:计算Ⅱ区径向流阻力;
式中:Rnij为第i渗透带中第j个渗流单元当量井径向流阻力,单位为MPa2/((m3/d)·(mPa·s));kfij为裂缝渗透率,单位为Dc;wfij为裂缝宽度,单位为cm;
c、根据所述的Ⅰ区线性流的阻力和Ⅱ区径向流阻力,建立所述渗流带的渗流单元的模拟模型。
2.如权利要求1所述的方法,其特征在于,所述步骤1)中,所述储层的基本参数包括:原地应力方向、储层厚度、储层温度、孔隙度、渗透率;所述流体的基本参数包括:气体粘度、气体临界压力、气体偏差因子、气体临界温度;所述水平井筒的基本参数包括:水平井方位、水平井井筒长度。
3.如权利要求2所述的方法,其特征在于,所述步骤2)具体为:根据储层在水平井筒长度方向上渗透率的差异,将渗透率相同的储层划分成同一个含有人工裂缝渗透带的渗流带。
4.如权利要求3所述的方法,其特征在于,所述渗流单元的基本流动过程包括:基质线性流、当量井径向流和裂缝内的线性渗流;以第1个条带基质中的第1个渗流单元为例计算渗流单元中各流动过程的阻力。
5.如权利要求4所述的方法,其特征在于,所述步骤a中,所述关系式中,Ruij=ψuij-ψdij,
根据气藏启动压力梯度下的基质线性流动流量与拟压力关系,得出公式:
根据拟压力定义,气体流过距离d时消耗的启动压力为Gij·d,将Gij·d定义为启动压力PB;因此可做以下转化:
式中:Ψuij、Ψdij分别为Ⅰ区基质线性流外拟压力和内拟压力,Ψ Buij - Ψ Bdij 为流体流过长度为dij的Ⅰ区基质的启动拟压差,单位为MPa2/(mPa·s);Gij为Ⅰ区基质的启动压力梯度,单位为MPa/m;dij为Ⅰ区基质线性流长度,单位为m;qij为第i渗透带中第j条裂缝流量,单位为m3/d;xf为裂缝半长,单位为m;ki为第i渗透带基质渗透率,单位为mD;h为储层厚度,单位为m;T为气藏温度,单位为K;Tsc为标准大气温度,单位为K;psc为标准大气压,单位为MPa。
6.如权利要求5所述的方法,其特征在于,所述步骤b中:
根据同一个渗流单元内,I区和II区边界位置处的压力相等,Ⅰ区线性流的内压即为Ⅱ区径向流的外压,则径向流流量与压差关系式为:
式中:Ψwfij为Ⅱ区径向流内拟压力,单位为MPa2/(mPa·s);rw为当量直井半径,单位为m;
所述当量直井的半径rw通过当量井径模型求取,具体为:
通过建立裂缝产量和直井产量相等的关系式从而得到当量井径与裂缝半长之间关系式,下面分别建立求取当量直径和裂缝的产量:
直井产量:求取定压供给边界条带气藏中一口直井A1的稳态产量,根据镜像反映原理,将其转化为一排生产井和一排注气井的产量稳态进行求解;
根据等值渗流阻力法,得一口井产量:
有限导流裂缝产量:
式中:Ψ为裂缝方向上x点处拟压力,单位为MPa2/(mPa·s);kfij为裂缝渗透率,单位为Dc;wfij为裂缝宽度,单位为cm;
解出裂缝内平均压力近似为裂缝内压力,流体从基质流到裂缝内满足线性流动规律,线性流外拟压力为Ψuij,求出基质线性流流量,再由质量守恒得裂缝产量:
对比当量直井的产量式(6)与有限导流裂缝的产量式(8)得当量井半径公式:
式中,rw为当量直井半径,单位为m;
将式(8)代入式(5)可以得到渗流单元Ⅱ区当量井附近等效径向流阻力为:
7.如权利要求6所述的方法,其特征在于,所述步骤4)具体为:
在求得各区流动阻力之后,根据不同渗流区域间满足压力连续、流量相等的基本原理,利用等值渗流阻力法即求得压裂水平井产能公式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510889494.8A CN105422070B (zh) | 2015-12-07 | 2015-12-07 | 优化特低渗透非均质气藏压裂水平井裂缝位置的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510889494.8A CN105422070B (zh) | 2015-12-07 | 2015-12-07 | 优化特低渗透非均质气藏压裂水平井裂缝位置的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105422070A CN105422070A (zh) | 2016-03-23 |
CN105422070B true CN105422070B (zh) | 2018-11-23 |
Family
ID=55500565
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510889494.8A Active CN105422070B (zh) | 2015-12-07 | 2015-12-07 | 优化特低渗透非均质气藏压裂水平井裂缝位置的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105422070B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108256208B (zh) * | 2018-01-15 | 2021-07-27 | 中国石油化工股份有限公司 | 大斜度拟水平井分段压裂产能计算方法 |
CN109033674B (zh) * | 2018-08-09 | 2022-06-10 | 西南石油大学 | 一种靶向压裂酸化井裂缝参数优化方法 |
CN109236287A (zh) * | 2018-11-08 | 2019-01-18 | 中国石油大学(北京) | 浅层低渗油藏一层多水平缝压裂裂缝参数确定方法及装置 |
CN112765534A (zh) * | 2021-01-05 | 2021-05-07 | 中国石油天然气股份有限公司 | 二项式产能方程的应用方法、装置、终端及存储介质 |
CN116029232B (zh) * | 2023-03-28 | 2023-07-28 | 西南石油大学 | 一种底水气藏复杂井型水侵前沿预测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103256035A (zh) * | 2013-05-22 | 2013-08-21 | 中国石化集团华北石油局 | 一种致密气田水平井压裂缝地质设计方法 |
CN104500017A (zh) * | 2014-12-12 | 2015-04-08 | 中国石油天然气集团公司 | 一种优化水平井分段压裂位置的方法 |
CN104594872A (zh) * | 2015-01-04 | 2015-05-06 | 西南石油大学 | 一种优化致密气藏压裂水平井裂缝导流能力的方法 |
CN105019876A (zh) * | 2014-04-24 | 2015-11-04 | 中国石油化工股份有限公司 | 分段压裂水平井注水开发裂缝间距及井距确定方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2852881A1 (en) * | 2011-06-24 | 2012-12-27 | Board Of Regents, The University Of Texas System | Method for determining spacing of hydraulic fractures in a rock formation |
-
2015
- 2015-12-07 CN CN201510889494.8A patent/CN105422070B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103256035A (zh) * | 2013-05-22 | 2013-08-21 | 中国石化集团华北石油局 | 一种致密气田水平井压裂缝地质设计方法 |
CN105019876A (zh) * | 2014-04-24 | 2015-11-04 | 中国石油化工股份有限公司 | 分段压裂水平井注水开发裂缝间距及井距确定方法 |
CN104500017A (zh) * | 2014-12-12 | 2015-04-08 | 中国石油天然气集团公司 | 一种优化水平井分段压裂位置的方法 |
CN104594872A (zh) * | 2015-01-04 | 2015-05-06 | 西南石油大学 | 一种优化致密气藏压裂水平井裂缝导流能力的方法 |
Non-Patent Citations (1)
Title |
---|
"新马滚动区水平井分段压裂优化设计研究";于东;《中国优秀硕士学位论文全文数据库工程科技I辑》;20150915;第1页、第34页至第61页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105422070A (zh) | 2016-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105422070B (zh) | 优化特低渗透非均质气藏压裂水平井裂缝位置的方法 | |
CN104594872B (zh) | 一种优化致密气藏压裂水平井裂缝导流能力的方法 | |
CN105422071B (zh) | 评价低渗透非均质气藏压裂水平井裂缝参数合理性的方法 | |
CN105507893A (zh) | 预测页岩储层体积改造后生产效果的方法 | |
CN105484741A (zh) | 低渗透非均质应力敏感储层压裂水平井产量预测方法 | |
CN109408859B (zh) | 页岩气藏压裂水平井二维三重介质数值模型建立方法 | |
CN105386751B (zh) | 一种基于油藏渗流模型的水平井测井产能预测方法 | |
CN106321051B (zh) | 一种用于优化多段压裂水平井网络裂缝参数的方法 | |
RU2455471C1 (ru) | Система разработки монолитного малопродуктивного зонально-неоднородного нефтяного пласта | |
Rhein et al. | Channel fracturing in horizontal wellbores: the new edge of stimulation techniques in the eagle ford formation | |
CN109826609B (zh) | 一种针对多隔夹层油气藏水力压裂的支撑剂精准置放方法 | |
CN105350961A (zh) | 低渗透非均质应力敏感储层体积压裂水平井产量预测方法 | |
CN112434426A (zh) | 页岩气多级压裂水平井台阶梯度压降开发方法及装置 | |
CN106593400A (zh) | 一种超低渗致密油藏体积压裂定向井布井方法 | |
CN109146185B (zh) | 一种注采井网差异化水驱优化调整方法 | |
Economides et al. | Fracturing horizontal transverse, horizontal longitudinal and vertical wells: criteria for decision | |
CN113836695B (zh) | 一种基于无网格连接元的油藏数值模拟方法 | |
CN109101773A (zh) | 致密油水平井立体缝网簇网压裂优化方法 | |
CN106469333B (zh) | 一种低渗稠油油藏水平井热采压力分布预测方法 | |
CA2744749A1 (en) | Basal planer gravity drainage | |
CN105545275B (zh) | 气田气井部署方法和装置 | |
CN108005626B (zh) | 一种基于热管技术的天然气水合物开采装置及方法 | |
CN105205318A (zh) | 确定多层多段水平裂缝采油井的总产量的方法和装置 | |
WO2012145286A3 (en) | Well system with lateral main bore and strategically disposed lateral bores and method of forming | |
CN106096858B (zh) | 一种浅层稠油油藏地下原油粘度获取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |