CN107725035B - 一种用于非均匀产液水平井判别来水方向的压力监测方法 - Google Patents
一种用于非均匀产液水平井判别来水方向的压力监测方法 Download PDFInfo
- Publication number
- CN107725035B CN107725035B CN201710720472.8A CN201710720472A CN107725035B CN 107725035 B CN107725035 B CN 107725035B CN 201710720472 A CN201710720472 A CN 201710720472A CN 107725035 B CN107725035 B CN 107725035B
- Authority
- CN
- China
- Prior art keywords
- liquid production
- horizontal well
- well
- uniform liquid
- water injection
- 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
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 189
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 187
- 239000007788 liquid Substances 0.000 title claims abstract description 134
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000012544 monitoring process Methods 0.000 title claims abstract description 37
- 238000002347 injection Methods 0.000 claims abstract description 162
- 239000007924 injection Substances 0.000 claims abstract description 162
- 230000008859 change Effects 0.000 claims abstract description 29
- 238000004458 analytical method Methods 0.000 claims abstract description 27
- 238000012360 testing method Methods 0.000 claims abstract description 11
- 239000012530 fluid Substances 0.000 claims description 26
- 230000035699 permeability Effects 0.000 claims description 21
- 239000000243 solution Substances 0.000 claims description 19
- 230000015572 biosynthetic process Effects 0.000 claims description 18
- 238000004891 communication Methods 0.000 claims description 14
- 238000003860 storage Methods 0.000 claims description 7
- 238000009826 distribution Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 238000009795 derivation Methods 0.000 claims description 4
- 238000013178 mathematical model Methods 0.000 claims description 3
- 230000002093 peripheral effect Effects 0.000 claims description 3
- 238000011084 recovery Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 2
- 238000007906 compression Methods 0.000 claims description 2
- 150000001875 compounds Chemical class 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000003129 oil well Substances 0.000 description 3
- 230000007423 decrease Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 235000020681 well water Nutrition 0.000 description 2
- 239000002349 well water Substances 0.000 description 2
- YBJHBAHKTGYVGT-ZKWXMUAHSA-N (+)-Biotin Chemical compound N1C(=O)N[C@@H]2[C@H](CCCCC(=O)O)SC[C@@H]21 YBJHBAHKTGYVGT-ZKWXMUAHSA-N 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- FEPMHVLSLDOMQC-UHFFFAOYSA-N virginiamycin-S1 Natural products CC1OC(=O)C(C=2C=CC=CC=2)NC(=O)C2CC(=O)CCN2C(=O)C(CC=2C=CC=CC=2)N(C)C(=O)C2CCCN2C(=O)C(CC)NC(=O)C1NC(=O)C1=NC=CC=C1O FEPMHVLSLDOMQC-UHFFFAOYSA-N 0.000 description 1
Images
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
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
-
- 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
- E21B47/00—Survey of boreholes or wells
- E21B47/10—Locating fluid leaks, intrusions or movements
-
- 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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
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)
- Geophysics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明提供了一种一种用于非均匀产液水平井判别来水方向的压力监测方法，包括以下步骤：监测非均匀产液水平井井底压力变化数据；考虑周围注水井工作制度的变化，设计注水井与非均匀产液水平井之间的干扰试井；建立注水井干扰下的非均匀产液水平井井底压力分析图版；利用得到的非均匀产液水平井井底压力分析图版对测得的井底压力数据进行拟合分析，得到不同注水井改变注水量下的储层参数，通过分析比较不同注水井工作制度下拟合得到的参数，判别来水方向。本发明由于考虑了注水井注水量的变化，能够准确地判别非均匀产液水平井来水方向，同时利用该方法进行非均匀产液水平井监测压力解释更符合实际情况，提高了压力监测解释的精准度。
Description
技术领域
本发明涉及油气开采领域，尤其涉及非均匀产液水平井判别来水方向的压力监测方法的领域。
背景技术
直井注水、水平井采油的联合注采方式在低渗油藏的开发中成为一种普遍现象，一般情况下水平井并非整个水平段都参与生产，即水平井是非均匀产液，且水平井产量递减很快，一旦见水含水率迅速上升，通过水平井井地电位法、水平井分段产液识别的试井分析方法等可以有效地识别非均匀产液水平井的出水位置。
在一个注采单元中，一口水平井对应多口注水井，虽然可以确定非均匀产液水平井的出水位置和出水量，但是现有的压力分析方法不能准确地判别这些水来源于哪一口注水井。
例如申请号为CN201621039603.3的中国专利，提供了一种判断水平井出水位置与来水方向的工艺管柱，包括油管和套管，油管末端设有导向堵丝，油管上对应各层段均从前至后依次设有水平封隔器、扶正器和压力计，其中，最后一个层段的扶正器和压力计之间设有智能开关器，油管内设有和各层段压力计连接的电缆，所述电缆与地面的上位机连接传输各层段压力计数据。该专利根据油井段数，利用封隔器将各层段卡封，各卡封段下入压力计，井下分段实时采集，判断油井出水位置；通过改变周围注水井工作制度，测试各层段压力变化情况，识别来水方向。但是上述专利只提供了工艺管柱，未提供判别来水方向的具体方法，而运用常规方法，现有的关于非均匀产液水平井的模型和方法没有考虑注水井的干扰影响，对于生产后期实测压力导数曲线下降的特征拟合效果不好。
发明内容
为了克服现有技术中的问题，本发明通过在直井注水、水平井采油的注采井网中，设计非均匀产液水平井与周围注水井的干扰试井，建立了一种新型的用于非均匀产液水平井判别来水方向的压力监测方法，能够确定井间连通情况，判别来水方向，对于低渗油藏调整注水井工作制度，提高油井产量，具有重要的意义。
本发明解决前述技术问题所采用的技术方案是：一种用于非均匀产液水平井判别来水方向的压力监测方法，包括以下步骤：
①监测非均匀产液水平井井底压力变化数据；
②考虑周围注水井工作制度的变化，设计注水井与非均匀产液水平井之间的干扰试井；
③建立注水井干扰下的非均匀产液水平井井底压力分析图版；
④利用得到的非均匀产液水平井井底压力分析图版对测得的井底压力数据进行拟合分析，得到不同注水井改变注水量下的储层参数，通过分析比较不同注水井工作制度下拟合得到的参数，确定井间连通情况，判别来水方向。
优选的是，步骤①中，非均匀产液水平井井底压力变化数据由压力计测得，所述压力计下入位置为非均匀产液水平井造斜位置。
上述任一方案优选的是，步骤②中，所述注水井工作制度的变化具体包括改变非均匀产液水平井周围注水井的注入量。
上述任一方案优选的是，所述改变非均匀产液水平井周围注水井的注入量包括使非均匀产液水平井以某一产量生产，改变其周围某一口注水井的注入量，保持其他注水井的注入量不变，测取非均匀产液水平井的井底压力；依次改变该非均匀产液水平井周围其他注水井的注入量，同时保持其他注水井的注入量不变，测取非均匀产液水平井在注水井不同工作制度下的井底压力。注入量的变化理论上是任意的，范围不限，适用于实际生产中所有的变化情况，开始是多少和变化后是多少可根据实际情况而定。
上述任一方案优选的是，步骤③中，进一步包括根据测井和完井地质资料，设置模型的初始参数，根据初始参数，计算出井底压力随生产时间的变化值，绘制典型图版。
上述任一方案优选的是，所述初始参数通过在地质资料直接获取，包括产层厚度、孔隙度、地层流体粘度、综合压缩系数、体积系数、井半径、非均匀产液水平井产量、产液段数、各段产液量、原始地层压力。
上述任一方案优选的是，井底压力随生产时间的变化值即非均匀产液水平井井底压力降，包括非均匀产液水平井生产时造成的井底压力降和由于周围注水井注水在非均匀产液水平井井底造成的压力恢复。
上述任一方案优选的是，利用叠加原理建立数学模型，所述井底压力随生产时间的变化值满足下列关系式，
式中，
△p——井底压力差，atm；
pi——原始地层压力，atm；
pwf——井底压力，atm；
Kh——水平渗透率，D；
Kv——垂直渗透率，D；
φ——孔隙度，小数；
μ——地层流体粘度，cP；
Ct——综合压缩系数，atm-1；
η——导压系数，D·atm/cP；
B——体积系数，无因次；
qwi——第i个生产段的产量，cm3/s；
Lwi——第i个生产段的长度，cm；
h——产层厚度，cm；
rw——井半径，cm；(xwi,ywi,zwi)——第i个生产段中心点的坐标位置；
Gx，Gyz——x，yz方向上的源函数；
qj——第j个注水井的注入量，cm3/s；
(xj,yj,zj)——第j个注水井的井底坐标位置；
t——时间，s；
τ——时间变量，s。
上述任一方案优选的是，所述关系式的推导过程为
(1)求得非均匀产液水平井单独生产时的压力解
引入
则水平井三维渗流偏微分方程可化为如下标准形式：
考虑水平井为圆柱源，利用瞬时源解和Newman乘积法，得到第i个产液段的地层压力分布公式：
其中，
考虑各个产液段不同的表皮系数Si后，N个产液段水平井地层压力为：
(2)求得周围注水井注水时的压力解
每一口注水井的渗流微分方程为：
利用瞬时源解和叠加原理求得M口注水井地层压力为：.
(3)求得注水井干扰下的非均匀产液水平井的井底压力解利用叠加原理，求得注水井干扰下的非均匀产液水平井的井底压力为：
上述推导过程中，瞬时源解具体为：
瞬时点源压力表达式：
无限大平面直线源压力解：
无限大平面条带源压力解：
条带形地层中直线源压力解：
Newman乘积法：
该方法最初用于求解三维热传导问题，由于热传导问题与渗流问题在数学上具有某些共性，Newman乘积法也可被应用于多维渗流问题的瞬时源函数，即一个多维油藏可看做二至三个一维“油藏”的相交，或一维与二维油藏的相交，其多维源函数是多个一维源函数的乘积。相交以后，几个一维共有的流动区域构成一个多维流动区域；几个一维共有的源汇分布区构成一个多维的源汇分布区；边界条件则由一维边界条件的适当组合给出。
上述任一方案优选的是，所述步骤④进一步包括，比较和判断实测井底压力随时间的变化与计算理论值是否相等，若实测井底压力差随时间的变化关系与计算理论值相等，则将初始参数作为非均匀产液水平井的实际拟合参数；若不等，则修正拟合参数，直至两者在一个可接受的误差范围之内，得到非均匀产液水平井实际储层参数，进而确定井间连通情况，判别来水方向。拟合得到的导压系数、渗透率等参数比较大，说明井间连通情况较好，拟合得到的导压系数、渗透率等参数比较小，说明井间连通情况较差。参数值的大小即可反映井间连通情况；通过比较不同注水井工作制度下拟合得到的参数大小来判别来水方向。例如当改变某口井的注水量时，通过拟合得到的渗透率、导压系数等参数比改变其它井的注水量时拟合的结果大，则可判别该水平井的产水主要来自该注水井。
上述任一方案优选的是，所述储层参数包括储层渗透率、导压系数、井筒存储系数、表皮系数。
上述任一方案优选的是，步骤④中，所述对测得的井底压力数据进行拟合分析包括利用步骤③得到的非均匀产液水平井井底压力分析图版对测得的井底压力数据进行手动拟合分析或者使用计算机编程进行自动拟合分析。
上述任一方案优选的是，所述拟合得到的参数包括储层渗透率、导压系数、流动系数、各产液段长度和出水量、表皮系数、井筒存储系数。
本发明的方法简单易懂，易于操作，提高了生产效率。由于考虑了注水井注水量的变化，因此，采用该方法能够准确地判别非均匀产液水平井来水方向，同时利用该方法进行非均匀产液水平井监测压力解释更符合实际情况，提高了压力监测解释的精准度，为低渗油藏非均匀产液水平井及其周围注水井工作制度的调整提供了理论依据。具体体现在：
(1)提出了一种用于非均匀产液水平井判别来水方向的压力监测方法。
(2)考虑了周围注水井工作制度的变化，更符合实际。
(3)解释了生产后期实测压力导数曲线下降的问题。现有的关于非均匀产液水平井随时间的压力变化模型和方法没有考虑注水井的干扰影响，对于生产后期实测压力导数曲线下降的特征拟合效果不好；而本方法考虑了注水井干扰，提出的新方法能够和实测数据进行很好的拟合，因此可以解释生产后期实测压力导数曲线下降是由于注水井的干扰影响。
(4)在考虑注水井干扰的情况下，现有的技术无法实现非均匀产液水平井压力动态的拟合，本发明的方法能拟合好非均匀产液水平井压力动态，使理论曲线和实测数据能够近似重合在一起。
附图说明
图1为本发明的一种用于非均匀产液水平井判别来水方向的压力监测方法的一优选实施例的流程图。
图2为本发明的一种用于非均匀产液水平井判别来水方向的压力监测方法的一优选实施例的压力计下入位置示意图。
图3为本发明的一种用于非均匀产液水平井判别来水方向的压力监测方法的一优选实施例的注水井干扰下的非均匀产液水平井井底压力分析图版。
图4为本发明的一种用于非均匀产液水平井判别来水方向的压力监测方法的一优选实施例的非均匀产液水平井试井模型(单井)与实测数据拟合示意图。
图5为本发明图4所示实施例的采用本发明方法计算的理论曲线与实测数据拟合示意图。
具体实施方式
为了更进一步了解本发明的发明内容，下面将结合具体实施例对本发明作更为详细的描述，实施例只对本发明具有示例性作用，而不具有任何限制性的作用；任何本领域技术人员在本发明的基础上作出的非实质性修改，都应属于本发明保护的范围。
实施例1
如图1所示，在直井注水、非均匀产液水平井采油的注采井网中(例如五点法井网，七点法井网)，本实施例以五点法井网为例，实施一种用于非均匀产液水平井判别来水方向的压力监测方法，包括以下步骤：
①监测非均匀产液水平井井底压力变化数据；
②考虑周围注水井工作制度的变化，设计注水井与非均匀产液水平井之间的干扰试井；
③建立注水井干扰下的非均匀产液水平井井底压力分析图版；
④利用得到的非均匀产液水平井井底压力分析图版对测得的井底压力数据进行拟合分析，得到不同注水井改变注水量下的储层参数，通过分析比较不同注水井工作制度下拟合得到的参数，确定井间连通情况，判别来水方向。
如图2所示，在本实施例中，步骤①中，非均匀产液水平井井底压力变化数据由压力计测得，所述压力计下入位置为非均匀产液水平井造斜位置。
在本实施例中，步骤②中，所述注水井工作制度的变化具体包括改变非均匀产液水平井周围注水井的注入量。
在本实施例中，所述改变非均匀产液水平井周围注水井的注入量包括使非均匀产液水平井以某一产量生产，五点法井网中注水井分别编号A、B、C、D，改变其周围A注水井的注入量，保持其他B、C、D三口注水井的注入量不变，测取非均匀产液水平井的井底压力；然后依次改变注水井B、C、D的注入量，同时保持其他三口注水井的注入量不变，测取非均匀产液水平井在注水井不同工作制度下的井底压力。
在本实施例中，步骤③中，进一步包括根据测井和完井地质资料，设置模型的初始参数，根据初始参数，计算出井底压力随生产时间的变化值，绘制典型图版。如图3所示，本实施例中的图版是在五点法井网中，根据油田实际情况，首先给定一些基础参数值，初始参数的给定可根据油藏实际资料给出，也可先采用常规水平井模型拟合，将其结果作为初始参数，然后改变四口注水井注水量，利用本实施例中的下述井底压力随生产时间的变化公式计算出6种不同注水量下的理论井底压力和压力导数，将其绘制在双对数坐标系中。
在本实施例中，所述初始参数通过在地质资料直接获取，包括产层厚度、孔隙度、地层流体粘度、综合压缩系数、体积系数、井半径、非均匀产液水平井产量、产液段数、各段产液量、原始地层压力。
在本实施例中，井底压力随生产时间的变化值即非均匀产液水平井井底压力降，包括非均匀产液水平井生产时造成的井底压力降和由于周围注水井注水在非均匀产液水平井井底造成的压力恢复。
在本实施例中，利用叠加原理建立数学模型，所述井底压力随生产时间的变化值满足下列关系式，
式中，
△p——井底压力差，atm；
pi——原始地层压力，atm；
pwf——井底压力，atm；
Kh——水平渗透率，D；
Kv——垂直渗透率，D；
φ——孔隙度，小数；
μ——地层流体粘度，cP；
Ct——综合压缩系数，atm-1；
η——导压系数，D·atm/cP；
B——体积系数，无因次；
qwi——第i个生产段的产量，cm3/s；
Lwi——第i个生产段的长度，cm；
h——产层厚度，cm；
rw——井半径，cm；(xwi,ywi,zwi)——第i个生产段中心点的坐标位置；
Gx，Gyz——x，yz方向上的源函数；
qj——第j个注水井的注入量，cm3/s；
(xj,yj,zj)——第j个注水井的井底坐标位置；
t——时间，s；
τ——时间变量，s。
在本实施例中，所述关系式的推导过程为
(1)求得非均匀产液水平井单独生产时的压力解
引入
则水平井三维渗流偏微分方程可化为如下标准形式：
考虑水平井为圆柱源，利用瞬时源解和Newman乘积法，得到第i个产液段的地层压力分布公式：
其中，
考虑各个产液段不同的表皮系数Si后，N个产液段水平井地层压力为：
(2)求得周围注水井注水时的压力解
每一口注水井的渗流微分方程为：
利用瞬时源解和叠加原理求得M口注水井地层压力为：.
(3)求得注水井干扰下的非均匀产液水平井的井底压力解
利用叠加原理，求得注水井干扰下的非均匀产液水平井的井底压力为：
在本实施例中，所述步骤④进一步包括，比较和判断实测井底压力随时间的变化与计算理论值是否相等，若实测井底压力差随时间的变化关系与计算理论值相等，则将初始参数作为非均匀产液水平井的实际拟合参数；若不等，则修正拟合参数，直至两者在一个可接受的误差范围之内，得到非均匀产液水平井实际储层参数，进而确定井间连通情况，判别来水方向。拟合得到的导压系数、渗透率等参数比较大，说明井间连通情况较好，拟合得到的导压系数、渗透率等参数比较小，说明井间连通情况较差。参数值的大小即可反映井间连通情况；通过比较不同注水井工作制度下拟合得到的参数大小来判别来水方向。
在本实施例中，所述储层参数包括储层渗透率、导压系数、井筒存储系数、表皮系数。
在本实施例中，步骤④中，所述对测得的井底压力数据进行拟合分析包括利用步骤③得到的非均匀产液水平井井底压力分析图版对测得的井底压力数据进行手动拟合分析或者使用计算机编程进行自动拟合分析。
在本实施例中，所述拟合得到的参数包括储层渗透率、导压系数、流动系数、各产液段长度和出水量、表皮系数、井筒存储系数。
本发明的方法简单易懂，易于操作，提高了生产效率。由于考虑了注水井注水量的变化，因此，采用该方法能够准确地判别非均匀产液水平井来水方向，同时利用该方法进行非均匀产液水平井监测压力解释更符合实际情况，提高了压力监测解释的精准度，为低渗油藏非均匀产液水平井及其周围注水井工作制度的调整提供了理论依据。
实施例2
如图1所示，一种用于非均匀产液水平井判别来水方向的压力监测方法，包括以下步骤：
①监测非均匀产液水平井井底压力变化数据；
②考虑周围注水井工作制度的变化，设计注水井与非均匀产液水平井之间的干扰试井；
③建立注水井干扰下的非均匀产液水平井井底压力分析图版；
④利用得到的非均匀产液水平井井底压力分析图版对测得的井底压力数据进行拟合分析，得到不同注水井改变注水量下的储层参数，通过分析比较不同注水井工作制度下拟合得到的参数，确定井间连通情况，判别来水方向。
在本实施例中，优选油田一口非均匀产液水平井，该非均匀产液水平井投产初期产液量为10m3/d。该水平井周围有4口注水井，分别编号注水井1、注水井2、注水井3、注水井4，其注水量分别为20、25、30、15m3/d，与非均匀产液水平井的距离分别为420、446、462、450m。
非均匀产液水平井基本数据如表1所示。
表1某致密油田一口非均匀产液水平井基础参数
按照上述方法，首先，将压力计下入非均匀产液水平井造斜位置，实时监测非均匀产液水平井开井生产时井底压力变化值。记录不同注水井工作制度下的非均匀产液水平井测压数据，将测得(或折算)的井底压力差与时间以及井底压力差与时间的导数曲线绘制在双对数坐标系下。
其次，给定一个初始参数值，利用下述公式进行计算，得到双对数坐标系下的压力及压力导数图版。所述公式为：
式中，
△p——井底压力差，atm；
pi——原始地层压力，atm；
pwf——井底压力，atm；
Kh——水平渗透率，D；
Kv——垂直渗透率，D；
φ——孔隙度，小数；
μ——地层流体粘度，cP；
Ct——综合压缩系数，atm-1；
η——导压系数，D·atm/cP；
B——体积系数，无因次；
qwi——第i个生产段的产量，cm3/s；
Lwi——第i个生产段的长度，cm；
h——产层厚度，cm；
rw——井半径，cm；(xwi，ywi，zwi)——第i个生产段中心点的坐标位置；
Gx，Gyz——x，yz方向上的源函数；
qj——第j个注水井的注入量，cm3/s；
(xj，yj，zj)——第j个注水井的井底坐标位置；
t——时间，s；
τ——时间变量，s。
然后，对比理论值与实际测试值之间的差异，若两者基本重合，说明预设的参数值合理；若两者差距较大，需继续给定初始值，不断调整，直至两者在合理的误差范围之内。在上述注水井注水量下，首先将实测数据与非均匀产液水平井单井试井模型的理论曲线进行拟合，如图4所示，曲线拟合效果较差，后期压力导数下降阶段拟合不上。然后将实测数据与本专利中方法计算的理论曲线进行拟合，如图5所示，曲线拟合很好，拟合结果为：井储系数C＝0.95m3/MPa，表皮系数S＝0.16，有效长度65.5m，水平渗透率Kh＝4.5×10-3μm2，各水平段产液量分别为3.46、6.54m3/d，流动系数
导压系数η＝2678.31×10-3μm2·MPa/(mPa·s)。
最后，判别来水方向。为判别非均匀产液水平井的来水方向，对其周围4口注水井的注水量实施动态调整。调整方式为：依次改变4口井的注水量(在某口井注水量改变的同时，其它3口井保持原来的注水量不变)，每口注水井的日注水量都在原来的基础上增加10m3。测取非均匀产液水平井在4种不同注水井工作制度下的井底压力，将4种情形下的实测数据与理论曲线进行拟合，得到解释参数，如表2所示。
表2 4种注水井工作制度下的拟合参数结果
从表2中可以看出，改变注水井2的注水量时，拟合得到的渗透率、流动系数和导压系数最大，所以可以判断此非均匀产液水平井的产水主要来源于注水井2。
尽管具体地参考其优选实施例来示出并描述了本发明，但本领域的技术人员可以理解，可以作出形式和细节上的各种改变而不脱离所附权利要求书中所述的本发明的范围。以上结合本发明的具体实施例做了详细描述，但并非是对本发明的限制。凡是依据本发明的技术实质对以上实施例所做的任何简单修改，均仍属于本发明技术方案的范围。
Claims (11)
1.一种用于非均匀产液水平井判别来水方向的压力监测方法，包括以下步骤：
①监测非均匀产液水平井井底压力变化数据；
②考虑周围注水井工作制度的变化，设计注水井与非均匀产液水平井之间的干扰试井；
③建立注水井干扰下的非均匀产液水平井井底压力分析图版；
④利用得到的非均匀产液水平井井底压力分析图版对测得的井底压力数据进行拟合分析，得到不同注水井改变注水量下的储层参数，通过分析比较不同注水井工作制度下拟合得到的参数，确定井间连通情况，判别来水方向；
步骤③中，进一步包括根据测井和完井地质资料，设置模型的初始参数，根据初始参数，计算出井底压力随生产时间的变化值，绘制典型图版；
井底压力随生产时间的变化值即非均匀产液水平井井底压力降，包括非均匀产液水平井生产时造成的井底压力降和由于周围注水井注水在非均匀产液水平井井底造成的压力恢复；
利用叠加原理建立数学模型，所述井底压力随生产时间的变化值满足下列关系式，
式中，△p——井底压力差，atm；
pi——原始地层压力，atm；
pwf——井底压力，atm；
Kh——水平渗透率，D；
Kv——垂直渗透率，D；
φ——孔隙度，小数；
μ——地层流体粘度，cP；
Ct——综合压缩系数，atm-1；
η——导压系数，D·atm/cP；
B——体积系数，无因次；
qwi——第i个生产段的产量，cm3/s；
Lwi——第i个生产段的长度，cm；
h——产层厚度，cm；
rw——井半径，cm；(xwi,ywi,zwi)——第i个生产段中心点的坐标位置；
Gx，Gyz——x，yz方向上的源函数；
qj——第j个注水井的注入量，cm3/s；
(xj,yj,zj)——第j个注水井的井底坐标位置；
t——时间，s；
τ——时间变量，s；
Si——第i产液段的表皮系数；
N——产液段数；
M——注水井口数。
2.根据权利要求1所述的用于非均匀产液水平井判别来水方向的压力监测方法，其特征在于：步骤①中，非均匀产液水平井井底压力变化数据由压力计测得，所述压力计下入位置为非均匀产液水平井造斜位置。
3.根据权利要求1所述的用于非均匀产液水平井判别来水方向的压力监测方法，其特征在于：步骤②中，所述注水井工作制度的变化具体包括改变非均匀产液水平井周围注水井的注入量。
4.根据权利要求3所述的用于非均匀产液水平井判别来水方向的压力监测方法，其特征在于：所述改变非均匀产液水平井周围注水井的注入量包括使非均匀产液水平井以某一产量生产，改变其周围某一口注水井的注入量，保持其他注水井的注入量不变，测取非均匀产液水平井的井底压力；依次改变该非均匀产液水平井周围其他注水井的注入量，同时保持其他注水井的注入量不变，测取非均匀产液水平井在注水井不同工作制度下的井底压力。
5.根据权利要求1所述的用于非均匀产液水平井判别来水方向的压力监测方法，其特征在于：所述初始参数通过在地质资料直接获取，包括产层厚度、孔隙度、地层流体粘度、综合压缩系数、体积系数、井半径、非均匀产液水平井产量、产液段数、各段产液量、原始地层压力。
6.根据权利要求1所述的用于非均匀产液水平井判别来水方向的压力监测方法，其特征在于：所述关系式中，
7.根据权利要求1所述的用于非均匀产液水平井判别来水方向的压力监测方法，其特征在于：所述关系式的推导过程为
(1)求得非均匀产液水平井单独生产时的压力解
引入
则水平井三维渗流偏微分方程可化为如下标准形式：
考虑水平井为圆柱源，利用瞬时源解和Newman乘积法，得到第i个产液段的地层压力分布公式：
其中，
考虑各个产液段不同的表皮系数Si后，N个产液段水平井地层压力为：
(2)求得周围注水井注水时的压力解
每一口注水井的渗流微分方程为：
利用瞬时源解和叠加原理求得M口注水井地层压力为：.
(3)求得注水井干扰下的非均匀产液水平井的井底压力解
利用叠加原理，求得注水井干扰下的非均匀产液水平井的井底压力为：
8.根据权利要求1所述的用于非均匀产液水平井判别来水方向的压力监测方法，其特征在于：所述步骤④进一步包括，比较和判断实测井底压力随时间的变化与计算理论值是否相等，若实测井底压力差随时间的变化关系与计算理论值相等，则将初始参数作为非均匀产液水平井的实际拟合参数；若不等，则修正拟合参数，直至两者在一个可接受的误差范围之内，得到非均匀产液水平井实际储层参数，进而确定井间连通情况，判别来水方向。
9.根据权利要求8所述的用于非均匀产液水平井判别来水方向的压力监测方法，其特征在于：所述储层参数包括储层渗透率、导压系数、井筒存储系数、表皮系数。
10.根据权利要求1所述的用于非均匀产液水平井判别来水方向的压力监测方法，其特征在于：步骤④中，所述对测得的井底压力数据进行拟合分析包括利用步骤③得到的非均匀产液水平井井底压力分析图版对测得的井底压力数据进行手动拟合分析或者使用计算机编程进行自动拟合分析。
11.根据权利要求10所述的用于非均匀产液水平井判别来水方向的压力监测方法，其特征在于：所述拟合得到的参数包括储层渗透率、导压系数、流动系数、各产液段长度和出水量、表皮系数、井筒存储系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710720472.8A CN107725035B (zh) | 2017-08-21 | 2017-08-21 | 一种用于非均匀产液水平井判别来水方向的压力监测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710720472.8A CN107725035B (zh) | 2017-08-21 | 2017-08-21 | 一种用于非均匀产液水平井判别来水方向的压力监测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107725035A CN107725035A (zh) | 2018-02-23 |
| CN107725035B true CN107725035B (zh) | 2020-08-28 |
Family
ID=61205222
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710720472.8A Active CN107725035B (zh) | 2017-08-21 | 2017-08-21 | 一种用于非均匀产液水平井判别来水方向的压力监测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107725035B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111178785B (zh) * | 2020-01-06 | 2023-01-24 | 中国海洋石油集团有限公司 | 一种排状井网下水平井变流量油水两相压力动态评价方法 |
| CN115680609B (zh) * | 2021-07-27 | 2024-06-18 | 中国石油天然气股份有限公司 | 一种基于注采井组试井资料确立连通关系的方法 |
| CN115929279B (zh) * | 2022-10-18 | 2023-08-08 | 常州大学 | 一种注采平衡下判断干扰强度的图版分析法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101050699B (zh) * | 2006-12-11 | 2011-01-26 | 大庆油田有限责任公司 | 一种利用同层垂向干扰试井获取厚油层垂向渗透率的方法及测试管柱 |
| US9063250B2 (en) * | 2009-08-18 | 2015-06-23 | Schlumberger Technology Corporation | Interference testing while drilling |
| CN206071468U (zh) * | 2016-09-06 | 2017-04-05 | 中国石油天然气股份有限公司 | 一种判断水平井出水位置与来水方向的工艺管柱 |
-
2017
- 2017-08-21 CN CN201710720472.8A patent/CN107725035B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN107725035A (zh) | 2018-02-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107725034B (zh) | 一种用于多级压裂水平井判别来水方向的压力监测方法 | |
| CN106837297B (zh) | 一种识别井间连通性及油水动态预测的方法 | |
| CA2974979C (en) | Multi-phase polymer apparent viscosity determination in polymer coreflood simulation study workflow | |
| CN109522634B (zh) | 一种致密气多段体积压裂水平井数值分析方法 | |
| CN105822298B (zh) | 基于产气指数的页岩气层绝对无阻流量的获取方法 | |
| CN104504604B (zh) | 一种定性气井井筒积液的方法 | |
| CN107503739B (zh) | 一种用于水平井判别来水方向的压力监测方法 | |
| CN107725035B (zh) | 一种用于非均匀产液水平井判别来水方向的压力监测方法 | |
| CN105464644B (zh) | 一种井口压力监测方法及控制系统 | |
| CN108661631B (zh) | 一种产量预测方法 | |
| CN103032064A (zh) | 一种钻井过程中气侵位置检测的方法及装置 | |
| DK179910B1 (en) | Method for estimating well productivity along a section of a well bore | |
| CN103510940B (zh) | 机械采油井工况综合诊断分析方法及装置 | |
| CN110847894B (zh) | 一种井下节流气井流压的确定方法 | |
| CN105089566A (zh) | 一种气井系统配产方法 | |
| CN109441415B (zh) | 基于邻井干扰的聚合物驱油藏测试井的试井解释方法 | |
| CN108614902A (zh) | 一种生产制度变化的页岩气井生产数据分析方法 | |
| Ismadi et al. | The use of combined static-and dynamic-material-balance methods with real-time surveillance data in volumetric gas reservoirs | |
| CN110735633A (zh) | 一种低渗碳酸盐岩气藏气井井筒积液早期判识方法 | |
| CN108106687B (zh) | 一种含软夹层的基岩地下水流网探究方法及双胶囊止水器 | |
| CN112796725B (zh) | 一种分段压裂页岩气井压裂段产气贡献率确定方法及系统 | |
| CN111950112A (zh) | 一种适用于底部封闭的碳酸盐岩储层动态分析方法 | |
| Kumar et al. | Machine learning applications for a qualitative evaluation of the fracture network in the Wolfcamp shale using tracer and completion data | |
| CN112836438B (zh) | 一种非均质边水气藏气井见水时间预测方法及装置 | |
| CN103334740B (zh) | 考虑启动压力梯度的确定泄油前缘的方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |