CN109611088B - 一种诊断边水早期水侵的试井分析方法 - Google Patents
一种诊断边水早期水侵的试井分析方法 Download PDFInfo
- Publication number
- CN109611088B CN109611088B CN201811534935.2A CN201811534935A CN109611088B CN 109611088 B CN109611088 B CN 109611088B CN 201811534935 A CN201811534935 A CN 201811534935A CN 109611088 B CN109611088 B CN 109611088B
- Authority
- CN
- China
- Prior art keywords
- water
- pressure
- sector
- well
- area
- 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.)
- Expired - Fee Related
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 202
- 230000009545 invasion Effects 0.000 title claims abstract description 114
- 238000012360 testing method Methods 0.000 title claims abstract description 70
- 238000004458 analytical method Methods 0.000 title claims abstract description 26
- 238000003745 diagnosis Methods 0.000 claims abstract description 20
- 238000013178 mathematical model Methods 0.000 claims abstract description 16
- 230000009466 transformation Effects 0.000 claims abstract description 14
- 230000008878 coupling Effects 0.000 claims abstract description 8
- 238000010168 coupling process Methods 0.000 claims abstract description 8
- 238000005859 coupling reaction Methods 0.000 claims abstract description 8
- 238000000034 method Methods 0.000 claims abstract description 8
- 238000004519 manufacturing process Methods 0.000 claims description 21
- 239000002131 composite material Substances 0.000 claims description 15
- 230000035699 permeability Effects 0.000 claims description 11
- 230000000704 physical effect Effects 0.000 claims description 11
- 230000014509 gene expression Effects 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 4
- 238000003860 storage Methods 0.000 claims description 4
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 230000002500 effect on skin Effects 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 abstract description 3
- 230000008859 change Effects 0.000 abstract description 3
- 238000011161 development Methods 0.000 description 8
- 238000009792 diffusion process Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 230000001186 cumulative effect Effects 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 208000012661 Dyskinesia Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000002349 well water Substances 0.000 description 1
- 235000020681 well water Nutrition 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
- 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
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Geology (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Geochemistry & Mineralogy (AREA)
- Environmental & Geological Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Databases & Information Systems (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Software Systems (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Geophysics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
- Examining Or Testing Airtightness (AREA)
Abstract
本发明公开了一种诊断边水早期水侵的试井分析方法，包括以下步骤：采用扇区化方式将气井流动区域进行区域划分并建立多区渗流数学模型；将所述多区渗流数学模型的参数并进行无因次化；对所述多区渗流数学模型作Laplace变换；通过多区耦合建立单井试井解释模型，计算获得边水舌进水侵条件下的理论解释图版；绘制压降典型曲线和双对数压力诊断诊断曲线；对比分析上述几种曲线诊断是否存在早期水侵。本发明的方法在气井见水之前，通过试井手段判断水侵的存在、估计水侵前缘的距离和水侵的区域范围，本发明的双对数压力诊断诊断曲线能够反映气井外围地层流动能力的变化，水区的流度远低于气区，通过多次试井对比分析就能够实现早期水侵诊断。
Description
技术领域
本发明涉及边水气藏开发技术领域，特别涉及一种诊断边水早期水侵的试井分析方法。
背景技术
气藏中直接与气层连通的地下水叫做气层水，如果气层厚度不大或构造较陡时气充满圈闭的高部位，而且水环绕在气藏的周缘，这种水称为边水，该气藏则形成了边水气藏。在边水气藏的开发过程中，边水侵入会造成气井出水，不仅增加气藏的开发、开采难度，而且还会造成气井产能损失，降低气藏采收率，影响气藏开发效益。因此，需要对边水气藏的水侵动态进行判断和诊断，通过试井手段判断水侵的存在、估计水侵前缘距离和水侵范围，为边水气藏制定科学、合理的开发方案提供理论支持，提高气藏开发效益。
发明内容
本发明的目的在于克服现有技术的不足，提供一种诊断边水早期水侵的试井分析方法。
本发明的目的是通过以下技术方案来实现的：
一种诊断边水早期水侵的试井分析方法，包括以下几个步骤：
S1、采用扇区化方式将气井流动区域划分为气藏流动扇区和水侵影响扇区，所述水侵影响扇区包括边水舌进扇区和近井气藏区，以水侵前缘处的压力与气藏流动扇区压力相同、压力梯度相同为耦合条件，在气藏流动扇区和水侵影响扇区在井点的压力相同、流量之和等于井的总流量的条件下，建立多区渗流数学模型，所述多区渗流数学模型包括单扇区流动模型和复合扇区流动模型；
S2、将所述多区渗流数学模型的参数并进行无因次化后，定义所述气井流动区域的物性：区域导压系数比η、区域流度比γ、边水舌进扇区系数β、无因次半径、无因次井筒半径、无因次水侵前缘半径、无因次时间和无因次压力；
S3、对所述单扇区流动模型作Laplace变换，获得气藏流动扇区的Laplace空间无因次压力解；对所述复合扇区流动模型作Laplace变换，再根据界面流量、界面压力、内边界定流量产出和外边界无限大的条件整理得到所述近井气藏区1在井筒处的定产压降Laplace空间无因次压力解；
S4、通过多区耦合建立单井试井解释模型，计算获得边水舌进水侵条件下的理论解释图版，将近井气藏区、气藏流动扇区的定产压力解释转换为定压流量解，通过两区产出流量求和获得完整单井模型的定压流量，将Laplace空间井筒压力进行数值反演获得实空间的井底压力；
S5、根据步骤S2-S4中的边水舌进扇区系数β绘制压降典型曲线；假设气区和水区的渗透率k相同，确定气藏流动扇区与边水舌进扇区的扩散系数比η、流度比λ，绘制出第一试井典型曲线，当边水舌进扇区系数β增加与舌进半径R减小同时发生时绘制出第二试井典型曲线；
S6、对比分析所述压降典型曲线和双对数压力诊断诊断曲线诊断是否存在早期水侵，所述双对数压力诊断诊断曲线包括第一试井典型曲线和第二试井典型曲线，具体分析过程如下：
分析压降典型曲线：压力导数线上翘，表明压降波及到水侵区域，水侵区域存在导致双对数压力诊断诊断曲线中、后期出现压力导数驼峰，扇区系数β增加对应边水舌进扇区的范围扩大，水侵影响阶段的压降则显著增加；
分析第一试井典型曲线：舌进半径R越大，水侵驼峰出现越晚，扇区系数β越大，水侵区域越大，水侵驼峰越高；
分析第二试井典型曲线：水侵驼峰前移和驼峰高度增加，表明水侵区域越大。
进一步地，所述气井流动区域的各物性的表达式为：
区域导压系数比η
区域流度比γ
无因次半径rD＝r/L
无因次井筒半径RwD＝Rw/L
无因次水侵前缘半径RD＝R/L
无因次时间
无因次压力
边水舌进扇区系数β
其中，k为渗透率、μ为粘度、φ为孔隙度、Ct为综合压缩系数、L为参考长度、为舌进半径R、ref为参考物性、α为边水舌进扇区的夹角。
进一步地，所述单扇区流动模型用以下方程式表示：
气藏流动扇区的控制方程
内边界定流量产出
外边界无限
初始条件
结合初始条件对(1)～(3)式作Laplace变换，获得气藏流动扇区3的Laplace压力解为：
其中，K0()为第二类变形贝塞尔函数，z为Laplace变量，
进一步地，所述复合扇区流动模型用以下方程式表示：
初始条件
初始条件
界面流量相等
界面压力相等
外边界无限大
内边界定流量产出
结合初始条件对(6)～(13)式作Laplace变换，得(6)和(7)式的通解分别为：
其中，A1、B1、A2、B2为常数，I0()为第一类变形贝塞尔函数，K0()为第二类变形贝塞尔函数，z为Laplace变量，
进一步地，所述获得近井气藏区在井筒处的定产压降Laplace空间无因次压力解的求解方程组为：
求解方程组，解得系数A1、B1、B2后，获得近井气藏区在井筒RwD处的定产压降Laplace空间无因次压力
为：
进一步地，所述步骤S4还包括以下步骤：
将完整单井模型的定压流量反转为完整单井模型的定产压力解，再进一步迭加含井储和表皮效应，获得的井筒压力。
本发明的有益效果是：
1)本发明采用扇区化方式建立气区以及水侵复合区的径向渗流数学模型，能正确解释评价水侵区域参数，利用模型解释出水侵的前缘位置及水侵区域范围的大小，为后续边水气藏制定科学、合理的开发方案提供了理论判据，提高气藏开发效益。
2)本发明的方法在气井见水之前，通过试井手段判断水侵的存在、估计水侵前缘的距离和水侵的区域范围，本发明的双对数压力诊断诊断曲线能够反映气井外围地层流动能力的变化，水区的流度远低于气区，通过多次试井对比分析就能够实现早期水侵诊断。
3)通过本发明的诊断边水早期水侵的试井分析方法，前后两次测试的解释参数，前后测试的早期动态参数、中期的渗透率和表皮系数解释值基本相同，准确度高。
附图说明
图1为本发明的诊断边水早期水侵的试井分析方法的流程框图；
图2为本发明的气井流动区扇区划分示意图；
图3为本发明的压降典型曲线图；
图4-5为本发明的第一试井典型曲线图；
图6为本发明的第二试井典型曲线图；
图7为本发明的第1次压恢测试的双对数压力诊断曲线图；
图8为本发明的第2次压恢测试的双对数压力诊断曲线图；
图9为本发明应用分析实例的m18井的累产动态图；
图10为本发明应用分析实例的m18井第一次双对数压力拟合曲线图；
图11为本发明应用分析实例的m18第二次井双对数压力拟合曲线图；
图中，1-近井气藏区，2-水侵区，3-气藏流动扇区。
具体实施方式
下面将结合实施例，对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域技术人员在没有付出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
参阅图1-11，本发明提供一种技术方案：
早期水侵的试井分析方法目标是在气井见水之前，通过试井手段判断水侵的存在、估计水侵前缘距离和水侵范围。试井曲线能够反映气井外围地层流动能力的变化，水区的流度远低于气区，通过多次试井的对比分析实现早期水侵诊断，但是常规的径向复合试井解释模型不适应边水水侵情况，不能正确解释评价水侵区域参数。本发明提出了一种描述水侵舌进的多区复合试井解释新模型，采用扇区化方式建立气区以及水侵复合区的径向渗流数学模型，考虑气水区域的流度、导压系数差异，通过多区耦合形成单井试井解释模型，计算获得边水舌进水侵条件下的理论解释图版，水侵区的存在导致双对数压力诊断诊断曲线中后期出现压力导数驼峰，利用新模型可解释出水侵的前缘位置及水侵区域范围的大小。
首先，如图2所示，假设条件设置如下：
将气井流动区域划分为气藏流动扇区3和水侵影响扇区，其中水侵影响扇区考虑水侵舌进，又划分为近井气藏区1和水侵区2，定义水侵前缘到井点的半径为水侵半径即气区半径R，水侵前缘即图1中近井气藏区1的弧形段。
假设在以下条件下建立本发明中的建立多区渗流数学模型：
在水侵前缘处的压力与气藏流动扇区3处压力相同、压力梯度相同为耦合条件；
气藏流动扇区3和水侵影响扇区在井点的压力相同，两个区域的流量贡献不同，流量之和等于井的总流量。
求解思路：在水侵前缘处的两个区，即近井气藏区1和水侵区2压力相同、压力梯度相同为耦合条件，两个区即气藏流动扇区3和水侵影响扇区在井点的压力一致作为多区耦合条件，两个区即气藏流动扇区3和水侵影响扇区流量之和等于井的总流量为内边界条件，建立和求解多区渗流数学模型。
如图1所示，本发明的一种诊断边水早期水侵的试井分析方法，具体操作包括以下几个步骤：
第一步：采用扇区化方式将气井流动区域划分为气藏流动扇区3和水侵影响扇区，所述水侵影响扇区包括边水舌进扇区2和近井气藏区1，以水侵前缘处的压力与气藏流动扇区3压力相同、压力梯度相同为耦合条件，在气藏流动扇区3和水侵影响扇区在井点的压力相同、流量之和等于井的总流量的条件下，建立多区渗流数学模型，所述多区渗流数学模型包括单扇区流动模型和复合扇区流动模型；
第二步：将所述多区渗流数学模型的参数并进行无因次化后，定义所述气井流动区域的物性：区域导压系数比η、区域流度比γ、边水舌进扇区系数β、无因次半径、无因次井筒半径、无因次水侵前缘半径、无因次时间和无因次压力；
进一步地，所述气井流动区域的各物性的表达式为：
区域导压系数比η
区域流度比γ
无因次半径rD＝r/L
无因次井筒半径RwD＝Rw/L
无因次水侵前缘半径RD＝R/L
无因次时间
无因次压力
边水舌进扇区系数β
其中，k为渗透率、μ为粘度、φ为孔隙度、Ct为综合压缩系数、L为参考长度、为舌进半径R、ref为参考物性、α为边水舌进扇区2的夹角。
第三步：对所述单扇区流动模型作Laplace变换，获得气藏流动扇区3的Laplace空间无因次压力解；对所述复合扇区流动模型作Laplace变换，再根据界面流量、界面压力、内边界定流量产出和外边界无限大的条件整理得到所述近井气藏区1在井筒处的定产压降Laplace空间无因次压力解；
进一步地，所述单扇区流动模型用以下方程式表示：
气藏流动扇区的控制方程
内边界定流量产出
外边界无限
初始条件
结合初始条件对(1)～(3)式作Laplace变换，获得气藏流动扇区3的Laplace压力解为：
其中，K0()为第二类变形贝塞尔函数，z为Laplace变量，
进一步地，所述复合扇区流动模型用以下方程式表示：
初始条件
初始条件
界面流量相等
界面压力相等
外边界无限大
内边界定流量产出
结合初始条件对(6)～(13)式作Laplace变换，得(6)和(7)式的通解分别为：
其中，A1、B1、A2、B2为常数，I0()为第一类变形贝塞尔函数，K0()为第二类变形贝塞尔函数，z为Laplace变量，
进一步地，所述获得近井气藏区(1)在井筒处的定产压降Laplace空间无因次压力解的求解方程组为：
求解方程组，解得系数A1、B1、B2后，获得近井气藏区1在井筒RwD处的定产压降Laplace空间无因次压力
为：
第四步：通过多区耦合建立单井试井解释模型，计算获得边水舌进水侵条件下的理论解释图版，将近井气藏区1、气藏流动扇区3的定产压力解释转换为定压流量解，通过两区产出流量求和获得完整单井模型的定压流量，再反转为完整单井模型的定产压力解，进一步迭加含井储和表皮效应，获得的井筒压力，最后将Laplace空间井筒压力进行数值反演获得实空间的井底压力PwD。
第五步：根据第二步到第五步中的边水舌进扇区系数β绘制压降典型曲线；假设气区和水区的渗透率k相同，确定气藏流动扇区3与边水舌进扇区2的扩散系数比η、流度比λ，绘制出第一试井典型曲线，当边水舌进扇区系数β增加与舌进半径R减小同时发生时绘制出第二试井典型曲线；所述第一试井典型曲线和第二试井典型曲线合称为双对数压力诊断诊断曲线。
如图3所示，边水舌进扇区系数β从0增加到1时的压降典型曲线计算结果表示在图3中，当压降波及到水侵
区时，水侵带的低渗特性导致压力导数线上翘；扇区系数β为0对应无水侵区的均质模型响应，扇区系数β为1对应四周水侵的常规2区径向复合模型响应；扇区系数β增加对应水侵舌进区范围的扩大，水侵影响阶段的压降则显著增加。
根据某气藏的物性参数见表1.1，假设气区、水区的渗透率k相同，确定气藏流动扇区3与水侵区的扩散系数比η、水气流度比λ，计算出第一试井典型曲线如图图4-5所示，边水舌进扇区2的存在导致压力导数线在中期出现第二个驼峰，如图4所示，水侵半径即气区半径R越大，水侵驼峰出现越晚；如图5所示，水侵区域越大即扇区系数β越大，水侵驼峰越高。计算发现水侵驼峰高度主要受扇区系数β、水气流度比γ影响，受扩散系数比η的影响较小。
表1.1某气藏的物性参数
按照井点周围射线状的流线分布形态，边水水侵区域逐步扩大并且水侵前缘越来越接近井点，当水侵扇区系数β增加与舌进半径即水侵半径R减小同时发生时的第二试井典型曲线，如图6所示，水侵舌进导致水侵驼峰前移和驼峰高度增加。通过上述2次实测试井曲线的对比，可以判断是否存在上述特征，作为早期水侵的试井诊断判据。
总的来说，对比分析所述压降典型曲线和双对数压力诊断诊断曲线诊断是否存在早期水侵，具体分析过程如下：
分析压降典型曲线：压力导数线上翘，表明压降波及到水侵区域，水侵区域存在导致双对数压力诊断诊断曲线中、后期出现压力导数驼峰，扇区系数β增加对应边水舌进扇区2的范围扩大，水侵影响阶段的压降则显著增加；
分析第一试井典型曲线：舌进半径R越大，水侵驼峰出现越晚，扇区系数β越大，水侵区域越大，水侵驼峰越高；
分析第二试井典型曲线：水侵驼峰前移和驼峰高度增加，表明水侵区域越大。
应用分析实例：
m18井第1次压恢测试时关井前累产气5984万方，第2次压恢测试时关井前累产气14128万方，两次测试施工成功、仪器正常，双对数压力诊断曲线上的反映井储控制早期动态相似，但反映地层响应的中期形态差异较大，如图7-8所示，测试期间周围邻井的工作制度无异常，可排除邻井的干扰异动，测试井本身处于边水附近，怀疑是边水舌进水侵：
第1次测试：图7的诊断曲线中期反映出“似复合”气藏或舌进水侵气藏特征，导数线上翘；
第2次测试：图8的诊断曲线中期反映出“似均质”气藏特征，当舌进水侵影响与早期井储影响叠加在一起时也反映出类似动态；
由m18井的累产动态图，如图9上可见，在累产气15298万方开始该井产水量快速上升，证实有边水突破。
采用舌进水侵模型解释，第一次和第二次测试的双对数压力拟合曲线分别见图10和图11，前后两次测试的解释参数对比见表1.2，前后测试的早期动态参数、中期的渗透率和表皮系数解释值基本相同，水侵影响主要表现在舌进水侵的前缘半径减大幅度小，反映出水侵推进，控制舌进区形态的扇区系数的略降，反映出水侵舌进区前缘更尖。
表1.2两次测试的解释结果对比
以上所述仅是本发明的优选实施方式，应当理解本发明并非局限于本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、修改和环境，并能够在本文所述构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附权利要求的保护范围内。
Claims (3)
1.一种诊断边水早期水侵的试井分析方法，其特征在于，包括以下步骤：
步骤S1、采用扇区化方式将气井流动区域划分为气藏流动扇区(3)和水侵影响扇区，所述水侵影响扇区包括边水舌进扇区(2)和近井气藏区(1)，以水侵前缘处的压力与气藏流动扇区(3)压力相同、压力梯度相同为耦合条件，在气藏流动扇区(3)和水侵影响扇区在井点的压力相同、流量之和等于井的总流量的条件下，建立多区渗流数学模型，所述多区渗流数学模型包括单扇区流动模型和复合扇区流动模型；
步骤S2、将所述多区渗流数学模型的参数并进行无因次化后，定义所述气井流动区域的物性：区域导压系数比η、区域流度比λ、边水舌进扇区系数β、无因次半径、无因次井筒半径、无因次水侵前缘半径、无因次时间和无因次压力；
步骤S3、对所述单扇区流动模型作Laplace变换，获得气藏流动扇区(3)的Laplace空间无因次压力解；对所述复合扇区流动模型作Laplace变换，再根据界面流量、界面压力、内边界定流量产出和外边界无限大的条件整理得到所述近井气藏区(1)在井筒处的定产压降Laplace空间无因次压力解；
步骤S4、通过多区耦合建立单井试井解释模型，计算获得边水舌进水侵条件下的理论解释图版，将近井气藏区(1)、气藏流动扇区(3)的定产压力解释转换为定压流量解，通过两区产出流量求和获得完整单井模型的定压流量，将Laplace空间井筒压力进行数值反演获得实空间的井底压力；
步骤S5、根据步骤S2-S4中的边水舌进扇区系数β绘制压降典型曲线；假设气区和水区的渗透率k相同，确定气藏流动扇区(3)与边水舌进扇区(2)的区域导压系数比η、区域流度比λ，绘制出第一试井典型曲线，当边水舌进扇区系数β增加与舌进半径R减小同时发生时绘制出第二试井典型曲线；
步骤S6、对比分析所述压降典型曲线和双对数压力诊断曲线诊断是否存在早期水侵，所述双对数压力诊断曲线包括第一试井典型曲线和第二试井典型曲线，具体分析过程如下：
分析压降典型曲线：压力导数线上翘，表明压降波及到水侵区域，水侵区域存在导致双对数压力诊断曲线中、后期出现压力导数驼峰，边水舌进扇区系数β增加对应边水舌进扇区(2)的范围扩大，水侵影响阶段的压降则显著增加；
分析第一试井典型曲线：舌进半径R越大，水侵驼峰出现越晚，边水舌进扇区系数β越大，水侵区域越大，水侵驼峰越高；
分析第二试井典型曲线：水侵驼峰前移和驼峰高度增加，表明水侵区域越大；
气井流动区域的各物性的表达式为：
区域导压系数比
区域流度比
无因次半径rD＝r/L
无因次井筒半径RwD＝Rw/L
无因次水侵前缘半径RD＝R/L
无因次时间
无因次压力
边水舌进扇区系数
其中，k为渗透率、μ为粘度、φ为孔隙度、Ct为综合压缩系数、L为参考长度、为舌进半径R、ref为参考物性、α为边水舌进扇区(2)的夹角；
单扇区流动模型用以下方程式表示：
气藏流动扇区的控制方程
内边界定流量产出
外边界无限
初始条件
结合初始条件对(1)～(3)式作Laplace变换，获得气藏流动扇区(3)的Laplace压力解为：
其中，K0(rDu3)为第二类变形贝塞尔函数，z为Laplace变量，
复合扇区流动模型用以下方程式表示：
初始条件
初始条件
界面流量相等
界面压力相等
外边界无限大
内边界定流量产出
结合初始条件对(6)～(13)式作Laplace变换，得(6)和(7)式的通解分别为：
其中，A1、B1、A2、B2为常数，I0(rDu1)、I0(rDu2)分别为第一类变形贝塞尔函数，K0(rDu1)、K0(rDu2)分别为第二类变形贝塞尔函数，
z为Laplace变量。
2.根据权利要求1所述的诊断边水早期水侵的试井分析方法，其特征在于：所述获得近井气藏区(1)在井筒处的定产压降Laplace空间无因次压力解的求解方程组为：
求解方程组，解得系数A1、B1、B2后，获得近井气藏区(1)在井筒RwD处的定产压降Laplace空间无因次压力
为：
3.根据权利要求1所述的诊断边水早期水侵的试井分析方法，其特征在于：所述步骤S4还包括以下步骤：
将完整单井模型的定压流量反转为完整单井模型的定产压力解，再进一步叠加含井储和表皮效应，获得的井筒压力。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811534935.2A CN109611088B (zh) | 2018-12-14 | 2018-12-14 | 一种诊断边水早期水侵的试井分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811534935.2A CN109611088B (zh) | 2018-12-14 | 2018-12-14 | 一种诊断边水早期水侵的试井分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109611088A CN109611088A (zh) | 2019-04-12 |
| CN109611088B true CN109611088B (zh) | 2022-04-22 |
Family
ID=66009278
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811534935.2A Expired - Fee Related CN109611088B (zh) | 2018-12-14 | 2018-12-14 | 一种诊断边水早期水侵的试井分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109611088B (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5090238A (en) * | 1990-09-27 | 1992-02-25 | Santa Fe Energy Resources, Inc. | Oil well production testing |
| AU2006235887A1 (en) * | 2000-02-22 | 2006-11-23 | Schlumberger Technology B.V. | Integrated reservoir optimization |
| CN105443120A (zh) * | 2015-10-26 | 2016-03-30 | 中国石油化工股份有限公司 | 一种海相整装气田边水水侵早期特征分析方法 |
| CN106503407A (zh) * | 2016-12-05 | 2017-03-15 | 中海石油(中国)有限公司 | 存在部分连通断层的线性水侵油藏的试井分析方法及装置 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8132453B2 (en) * | 2005-05-10 | 2012-03-13 | Schlumberger Technology Corporation | Method for analysis of pressure response in underground formations |
| CN105569641B (zh) * | 2015-12-17 | 2021-09-28 | 西南石油大学 | 一种改善近井压降的水平井瞬态压力快速计算模型建立方法 |
| US11280179B2 (en) * | 2017-01-30 | 2022-03-22 | Schlumberger Technology Corporation | Evaluation of pressure-transient behavior of wells |
| CN107563027B (zh) * | 2017-08-21 | 2019-07-23 | 西南石油大学 | 用于分析多缝洞单元串联结构的试井解释模型与方法 |
| CN107676086B (zh) * | 2017-10-27 | 2020-10-09 | 中国石油天然气股份有限公司 | 一种分散边底水型气藏水侵阶段的判别方法及其装置 |
-
2018
- 2018-12-14 CN CN201811534935.2A patent/CN109611088B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5090238A (en) * | 1990-09-27 | 1992-02-25 | Santa Fe Energy Resources, Inc. | Oil well production testing |
| AU2006235887A1 (en) * | 2000-02-22 | 2006-11-23 | Schlumberger Technology B.V. | Integrated reservoir optimization |
| CN105443120A (zh) * | 2015-10-26 | 2016-03-30 | 中国石油化工股份有限公司 | 一种海相整装气田边水水侵早期特征分析方法 |
| CN106503407A (zh) * | 2016-12-05 | 2017-03-15 | 中海石油(中国)有限公司 | 存在部分连通断层的线性水侵油藏的试井分析方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109611088A (zh) | 2019-04-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107563027B (zh) | 用于分析多缝洞单元串联结构的试井解释模型与方法 | |
| Prijambodo et al. | Well test analysis for wells producing layered reservoirs with crossflow | |
| Bourdet et al. | Use of pressure derivative in well-test interpretation | |
| CN111425193B (zh) | 一种基于聚类分析测井岩石物理相划分的储层可压性评价方法 | |
| CN102953726A (zh) | 一种水驱油田优势通道识别方法及装置 | |
| CN110043254B (zh) | 一种基于电缆地层测试资料地层有效渗透率的获取方法 | |
| CN105298479A (zh) | 压裂直井产油(气)位置的诊断方法及其系统 | |
| CN109057786A (zh) | 一种用于非常规油气藏可采储量估算的自适应模型 | |
| CN112377178B (zh) | 一种边水水侵前缘的定量诊断方法 | |
| CN111339481B (zh) | 一种煤层气井的试井分析方法 | |
| CN110821473B (zh) | 一种井筒管柱损伤评价方法 | |
| CN109611088B (zh) | 一种诊断边水早期水侵的试井分析方法 | |
| Aibazarov et al. | Well Spacing Verification At Gas Condensate Field Using Deconvolution Driven Long-Term Pressure and Rate Analysis | |
| CN109710881B (zh) | 一种非均质井区井间干扰快速解析计算方法 | |
| CN110486008B (zh) | 一种径向复合油藏的参数解释方法及系统 | |
| CN109184676B (zh) | 一种页岩气藏有效改造体积评价方法 | |
| CN112036097A (zh) | 一种水锁气井的产能计算方法 | |
| Zheng et al. | Extraction of interference from long-term transient pressure using multi-well deconvolution algorithm for well test analysis | |
| CN110552682A (zh) | 一种碳酸盐岩三重介质气藏直井产量递减分析方法 | |
| CN105550438A (zh) | 鱼骨状多分支水平井产油（气）位置的诊断方法及其系统 | |
| CN112700039B (zh) | 一种火电厂负荷运行数据的稳态检测与提取方法 | |
| Ogbe et al. | Pulse testing with wellbore storage and skin effects | |
| Ballin et al. | New reservoir dynamic connectivity measurement for efficient well placement strategy analysis under depletion | |
| CN118008259A (zh) | 一种考虑应力敏感的气井产能试井解释方法 | |
| Ahmed et al. | Production logging as an integral part of horizontal-well transient-pressure test |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220422 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |