CN107977732B - 边水气藏水平井的见水时间预测方法及装置 - Google Patents
边水气藏水平井的见水时间预测方法及装置 Download PDFInfo
- Publication number
- CN107977732B CN107977732B CN201710998036.7A CN201710998036A CN107977732B CN 107977732 B CN107977732 B CN 107977732B CN 201710998036 A CN201710998036 A CN 201710998036A CN 107977732 B CN107977732 B CN 107977732B
- Authority
- CN
- China
- Prior art keywords
- water
- gas
- reservoir
- target well
- breakthrough time
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 282
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000035699 permeability Effects 0.000 claims abstract description 42
- 238000004519 manufacturing process Methods 0.000 claims abstract description 39
- 239000012071 phase Substances 0.000 claims description 106
- 239000008346 aqueous phase Substances 0.000 claims description 23
- 230000005484 gravity Effects 0.000 claims description 17
- 230000033001 locomotion Effects 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 14
- 230000001133 acceleration Effects 0.000 claims description 12
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 239000007789 gas Substances 0.000 description 127
- 238000010586 diagram Methods 0.000 description 13
- 238000003860 storage Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 239000012530 fluid Substances 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 239000002343 natural gas well Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000010924 continuous production Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Agronomy & Crop Science (AREA)
- General Health & Medical Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Primary Health Care (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Husbandry (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Solid-Sorbent Or Filter-Aiding Compositions (AREA)
- Battery Electrode And Active Subsutance (AREA)
Abstract
本申请实施例提供了一种边水气藏水平井的见水时间预测方法及装置，该方法包括：确定目标井的特征参数，所述特征参数包括水相渗透率、气相渗透率、水相粘度、气相粘度、水相密度、气相密度、储层孔隙度、含水饱和度、储层生产压差、气体产量、气体体积系数、初始气水边界与所述目标井的距离、储层厚度、所述目标井的水平段长度和储层倾角；根据所述特征参数确定水舌点移动速度；根据所述水舌点移动速度确定所述目标井的见水时间。本申请实施例可提高边水气藏水平井的见水时间预测的准确度和通用性。
Description
技术领域
本申请涉及天然气井的见水风险预测技术领域，尤其是涉及一种边水气藏水平井的见水时间预测方法及装置。
背景技术
见水时间预测是气藏开发的重要问题，国内外学者做了大量研究，运用各种方法推导出不同见水时间预测模型，其中对含底水气藏见水时间的研究较多，而对边水气藏见水时间的研究相对较少。
目前已有运用多孔介质流体渗流理论推导出边水气藏见水时间预测模型，其预测是基于直井进行的推导，而目前水平井在气藏开发中的应用日益增多，因此，若使用上述模型对水平井进行计算，必然会造成较大误差。因此，现有天然气井的见水时间预测方案的准确度和适用性亟待提高。
发明内容
本申请实施例的目的在于提供一种边水气藏水平井的见水时间预测方法及装置，以提高边水气藏水平井的见水时间预测的准确度和通用性。
为达到上述目的，一方面，本申请实施例提供了一种边水气藏水平井的见水时间预测方法，包括：
确定目标井的特征参数，所述特征参数包括水相渗透率、气相渗透率、水相粘度、气相粘度、水相密度、气相密度、储层孔隙度、含水饱和度、储层生产压差、气体产量、气体体积系数、初始气水边界与所述目标井的距离、储层厚度、所述目标井的水平段长度和储层倾角；
根据所述特征参数确定水舌点移动速度；
根据所述水舌点移动速度确定所述目标井的见水时间。
本申请实施例的边水气藏水平井的见水时间预测方法，所述根据所述特征参数确定水舌点移动速度，包括：
根据公式
确定水舌点移动速度；
其中，Viw为水舌点移动速度，Kw为水相渗透率，Kg为气相渗透率，μw为水相粘度，μg为气相粘度，ρw为水相密度，ρg为气相密度，φ为储层孔隙度，Sw为含水饱和度，ΔP为储层生产压差，q为气体产量，Bg为气体体积系数，L为初始气水边界与所述目标井的距离，h为储层厚度，Lh为所述目标井的水平段长度，α为储层倾角，g为重力加速度。
本申请实施例的边水气藏水平井的见水时间预测方法，所述根据所述水舌点移动速度确定所述目标井的见水时间，包括：
根据公式
确定目标井的见水时间；
其中，tbt(0,t)为目标井的见水时间，L1为水舌点移动距离，Viw为水舌点移动速度，Kw为水相渗透率，Kg为气相渗透率，μw为水相粘度，μg为气相粘度，ρw为水相密度，ρg为气相密度，φ为储层孔隙度，Sw为含水饱和度，ΔP为储层生产压差，q为气体产量，Bg为气体体积系数，L为初始气水边界与所述目标井的距离，h为储层厚度，Lh为所述目标井的水平段长度，α为储层倾角，g为重力加速度。
本申请实施例的边水气藏水平井的见水时间预测方法，在在地层均质含倾角情况下，所述目标井的见水时间通过以下公式确定：
本申请实施例的边水气藏水平井的见水时间预测方法，所述储层生产压差通过以下公式确定：
其中，ΔP为储层生产压差，ΔPw为水区压差，ΔPg为气区压差，Kw为水相渗透率，Kg为气相渗透率，μw为水相粘度，μg为气相粘度，ρw为水相密度，ρg为气相密度，q为气体产量，Bg为气体体积系数，L为初始气水边界与所述目标井的距离，L1为水舌点移动距离，h为储层厚度，Lh为所述目标井的水平段长度，α为储层倾角，g为重力加速度。
另一方面，本申请实施例还提供了一种边水气藏水平井的见水时间预测装置，包括：
特征参数确定模块，用于确定目标井的特征参数，所述特征参数包括水相渗透率、气相渗透率、水相粘度、气相粘度、水相密度、气相密度、储层孔隙度、含水饱和度、储层生产压差、气体产量、气体体积系数、初始气水边界与所述目标井的距离、储层厚度、所述目标井的水平段长度和储层倾角；
移动速度确定模块，用于根据所述特征参数确定水舌点移动速度；
见水时间确定模块，用于根据所述水舌点移动速度确定所述目标井的见水时间。
本申请实施例的边水气藏水平井的见水时间预测装置，所述根据所述特征参数确定水舌点移动速度，包括：
根据公式
确定水舌点移动速度；
其中，Viw为水舌点移动速度，Kw为水相渗透率，Kg为气相渗透率，μw为水相粘度，μg为气相粘度，ρw为水相密度，ρg为气相密度，φ为储层孔隙度，Sw为含水饱和度，ΔP为储层生产压差，q为气体产量，Bg为气体体积系数，L为初始气水边界与所述目标井的距离，h为储层厚度，Lh为所述目标井的水平段长度，α为储层倾角，g为重力加速度。
本申请实施例的边水气藏水平井的见水时间预测装置，所述根据所述水舌点移动速度确定所述目标井的见水时间，包括：
根据公式
确定目标井的见水时间；
其中，tbt(0,t)为目标井的见水时间，L1为水舌点移动距离，Viw为水舌点移动速度，Kw为水相渗透率，Kg为气相渗透率，μw为水相粘度，μg为气相粘度，ρw为水相密度，ρg为气相密度，φ为储层孔隙度，Sw为含水饱和度，ΔP为储层生产压差，q为气体产量，Bg为气体体积系数，L为初始气水边界与所述目标井的距离，h为储层厚度，Lh为所述目标井的水平段长度，α为储层倾角，g为重力加速度。
本申请实施例的边水气藏水平井的见水时间预测装置，在在地层均质含倾角情况下，所述目标井的见水时间通过以下公式确定：
本申请实施例的边水气藏水平井的见水时间预测装置，所述储层生产压差通过以下公式确定：
其中，ΔP为储层生产压差，ΔPw为水区压差，ΔPg为气区压差，Kw为水相渗透率，Kg为气相渗透率，μw为水相粘度，μg为气相粘度，ρw为水相密度，ρg为气相密度，q为气体产量，Bg为气体体积系数，L为初始气水边界与所述目标井的距离，L1为水舌点移动距离，h为储层厚度，Lh为所述目标井的水平段长度，α为储层倾角，g为重力加速度。
另一方面，本申请实施例还提供了另一种边水气藏水平井的见水时间预测装置，包括存储器、处理器、以及存储在所述存储器上的计算机程序，所述计算机程序被所述处理器运行时执行如下步骤：
确定目标井的特征参数，所述特征参数包括水相渗透率、气相渗透率、水相粘度、气相粘度、水相密度、气相密度、储层孔隙度、含水饱和度、储层生产压差、气体产量、气体体积系数、初始气水边界与所述目标井的距离、储层厚度、所述目标井的水平段长度和储层倾角；
根据所述特征参数确定水舌点移动速度；
根据所述水舌点移动速度确定所述目标井的见水时间。
由以上本申请实施例提供的技术方案可见，本申请实施例在确定目标井的包含储层倾角和水平段长度等在内的特征参数的基础上，根据该特征参数确定水舌点移动速度，然后基于水舌点移动速度确定目标井的见水时间。由于本申请实施例充分考虑储层倾角和水平段长度等参数对见水时间的影响，从而使得本申请实施例预测出的边水气藏水平井的见水时间更加准确，并且，在储层倾角为零时，本申请实施例还可以用来计算无倾斜均质地层水平井的见水时间，从而使得本申请实施能适用于不同生产条件下水平井的见水时间预测，因而本申请实施提高了边水气藏水平井的见水时间预测的通用性。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本申请一实施例中边水气藏水平井的见水时间预测方法的流程图；
图2为本申请一实施例中倾斜边水气藏水平井的生产示意图；
图3为本申请一实施例中见水时间与水平段长度的关系示意图；
图4为本申请一实施例中见水时间与储层倾角的关系示意图；
图5为本申请一实施例中边水气藏水平井的见水时间预测装置的结构框图；
图6为本申请另一实施例中边水气藏水平井的见水时间预测装置的结构框图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
为便于本领域技术人员更加清楚的理解本申请，先对本申请实施例进行原理性描述。
图2为倾斜边水气藏水平井的生产示意图，由图2可知，某一带倾角的边水气藏，其中一口水平井为C，可以将不规则气水界面近似为一条直线。储层倾角为α，储层厚度为h，初始气水边界与水平井C相距为L，初始气水边界平面径向流动，时间t后渗流到A点，移动距离为L1，与水平井C的水平段距离为L2，水的渗流速度为Vw，水区压差为ΔPw，气的渗流速度为Vg，气区压差为ΔPg。
在满足倾斜地层气水运动变化的条件下，为简化模型推导过程，作以下假设：①倾斜储层等厚且均质；②气水渗流过程中，忽略毛管力的影响；③流体微可压缩，气水密度和黏度为常数；④气相、水相均满足达西定律。
根据气水两相渗流理论，考虑储层倾角的气水两相渗流公式为：
式中，Kw、Kg分别为水相和气相的渗透率，单位为μm2；μw、μg分别为水相和气相的粘度，单位为mPa·s；ρw、ρg分别为液体和气体的密度，单位为g/cm3；g为重力加速度，单位为m/s2。
水平井C的水平段沿轴向的过流截面积可表示为：
A＝Lhh (3)
式中，Lh为水平井C的水平段长度，单位为m。
根据流体力学基本理论，气体、液体渗流速度可表示为：
式中，q为气体产量，单位为m3/d；Bg为气体体积系数，无量纲。
联立式(1)、式(3)与式(5)，可求得水区压差为：
联立式(2)、式(3)与式(4)，可求得气区压差为：
联立式(6)和式(7)，储层生产压差可表示为：
由于L＝L1+L2，代入式(8)，进而可得：
由于水舌点移动速度与水相在储层孔隙中的流速相同，考虑孔隙度和含水饱和度的影响，则水舌点移动速度可表示为：
式中，Viw为修正后的水舌点移动速度，φ为储层孔隙度，Sw为含水饱和度。
将式(1)代入式(10)，可得：
联立式(8)和式(11)，可求得水舌点移动速度为：
经过整理可得：
令常数部分
则水舌点移动速度可表示为：
依据气水界面在孔隙介质中移动距离和移动时间的关系，初始时刻t＝0，气井见水时，t＝tbt，移动距离为L，则
将式(16)代入式(17)，经整理，气藏水平井见水时间可表示为
因此，在地层均质含倾角情况下，若需求得边水气藏水平井的见水时间，将式(14)、式(15)中的代入式(18)即可；若使α＝0，便可得到无倾斜均质地层水平井的计算模型。
由以上可知，影响边水气藏水平井见水时间的因素很多，例如水平井的水平段长度、储层倾斜角、边水距气井长度、气水流度比、气井产量和储层厚度等。下面就水平段长度和储层倾角这两个因素为例，对水平井边水气藏见水时间的影响程度进行分析。
一、水平段长度的影响
气体流速Vg和水平段长度Lh为负相关，由式(14)、式(15)可知，水平段长度Lh影响参数a和b，水平段越长，气藏的储量越大，在水平井产量一定的情况下流速就会降低，故舌进移动速度减小，见水时间增大。如图3所示，见水时间与水平段长度呈线性关系，水平井见水时间随水平段长度的增长而增长，水平段长度每增加10m，见水时间增加230d(即230天)。
二、储层倾角的影响
储层倾角主要是通过重力作用影响见水时间，重力是气水在多孔介质中流动的阻力，储层倾角越大，重力作用越大，导致流速就会降低，故舌进移动速度减小，见水时间增大。由图4可知，见水时间随着储层倾角增大而增大，但斜率先平缓后陡峭，两者之间呈现复杂非线性关系。当α＜5°时，重力作用不显著，见水时间增长缓慢；当α＞5°时，重力作用明显，见水时间增长迅速，每增加10°见水时间约为20d，在α＝90°时见水时间最大，此时变为底水气藏。
参考图1所示，在以上描述基础上，本申请实施例的边水气藏水平井的见水时间预测方法可以包括以下步骤：
S101、确定目标井的特征参数，所述特征参数包括水相渗透率、气相渗透率、水相粘度、气相粘度、水相密度、气相密度、储层孔隙度、含水饱和度、储层生产压差、气体产量、气体体积系数、初始气水边界与所述目标井的距离、储层厚度、所述目标井的水平段长度和储层倾角；
S102、根据所述特征参数确定水舌点移动速度。
S103、根据所述水舌点移动速度确定所述目标井的见水时间。
上述各步骤的细节请参见上述原理性描述，在此不再赘述。此外，虽然上文描述的过程流程包括以特定顺序出现的多个操作，但是，应当清楚了解，这些过程可以包括更多或更少的操作，这些操作可以顺序执行或并行执行(例如使用并行处理器或多线程环境)。
由此可见，本申请实施例以多孔介质流体渗流理论为基础，利用水平井渗流模型得到压力场分布，通过水质点舌进移动速度反演见水时间，推导出了考虑储层倾角和水平段长度的见水时间预测新模型，使得预测更为准确，且该模型在特定的条件下(即储层倾角为零)可转变为无倾斜均质地层水平井的计算模型，从而使得该见水时间预测新模型能适用于不同生产条件下水平井的见水时间预测，从而提高了其通用性。
下面介绍本申请一示例性实施例，具体如下：
西部某一边水气藏水平井A井和B井分别属于该气田的两个区块，两口井的现场基础数据如下表1和表2所示。
表1气井A现场基础数据表
表2气井B现场基础数据表
采用本申请实施例与现有技术的几个常用计算模型对现场基础数据进行计算和精度分析，其对应的计算结果如表3所示。
表3水平井边水气藏见水时间预测对比
由表3可以看出，王会强模型的预测精度最低，相对误差高达33％和36％，其主要原因是未考虑储层倾角和水平段长度的影响，由于储层倾角和水平段长度均会导致气体流速减小，舌进速度减小，见水时间增大，因此该模型计算时间较实际时间偏小；而李涛模型相对前者来说，预测精度略有提高，但相对误差仍很大，均超过27％，其主要原因是考虑了水平段长度但未考虑储层倾角的影响，因此该模型计算时间较实际时间偏小，且误差较大；而用本申请实施例所计算的结果相对误差很小，均低于10％，比现有技术的几种常用模型的计算精度都要高，从而与现场工况吻合度较高，进而验证了本申请实施例可准确预测边水气藏水平井的见水时间，有效指导现场持续生产。
参考图5所示，本申请实施例的边水气藏水平井的见水时间预测装置可以包括：
特征参数确定模块51，可以用于确定目标井的特征参数，所述特征参数包括水相渗透率、气相渗透率、水相粘度、气相粘度、水相密度、气相密度、储层孔隙度、含水饱和度、储层生产压差、气体产量、气体体积系数、初始气水边界与所述目标井的距离、储层厚度、所述目标井的水平段长度和储层倾角；
移动速度确定模块52，可以用于根据所述特征参数确定水舌点移动速度；
见水时间确定模块53，可以用于根据所述水舌点移动速度确定所述目标井的见水时间。
参考图6所示，本申请实施例的边水气藏水平井的见水时间预测装置，包括存储器、处理器、以及存储在所述存储器上的计算机程序，所述计算机程序被所述处理器运行时执行如下步骤：
确定目标井的特征参数，所述特征参数包括水相渗透率、气相渗透率、水相粘度、气相粘度、水相密度、气相密度、储层孔隙度、含水饱和度、储层生产压差、气体产量、气体体积系数、初始气水边界与所述目标井的距离、储层厚度、所述目标井的水平段长度和储层倾角；
根据所述特征参数确定水舌点移动速度；
根据所述水舌点移动速度确定所述目标井的见水时间。
为了描述的方便，描述以上装置时以功能分为各种单元分别描述。当然，在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白，本申请的实施例可提供为方法、系统或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。
Claims (7)
1.一种边水气藏水平井的见水时间预测方法，其特征在于，包括：
确定目标井的特征参数，所述特征参数包括水相渗透率、气相渗透率、水相粘度、气相粘度、水相密度、气相密度、储层孔隙度、含水饱和度、储层生产压差、气体产量、气体体积系数、初始气水边界与所述目标井的距离、储层厚度、所述目标井的水平段长度和储层倾角；
根据所述特征参数确定水舌点移动速度；
根据所述水舌点移动速度确定所述目标井的见水时间；
其中，所述根据所述特征参数确定水舌点移动速度，包括：
根据公式
确定水舌点移动速度；
所述根据所述水舌点移动速度确定所述目标井的见水时间，包括：
根据公式
确定目标井的见水时间；
其中，tbt(0,t)为目标井的见水时间，L1为水舌点移动距离，Viw为水舌点移动速度，Kw为水相渗透率，Kg为气相渗透率，μw为水相粘度，μg为气相粘度，ρw为水相密度，ρg为气相密度，φ为储层孔隙度，Sw为含水饱和度，ΔP为储层生产压差，q为气体产量，Bg为气体体积系数，L为初始气水边界与所述目标井的距离，h为储层厚度，Lh为所述目标井的水平段长度，α为储层倾角，g为重力加速度。
2.如权利要求1所述的边水气藏水平井的见水时间预测方法，其特征在于，在地层均质含倾角情况下，所述目标井的见水时间通过以下公式确定：
3.如权利要求1或2所述的边水气藏水平井的见水时间预测方法，其特征在于，所述储层生产压差通过以下公式确定：
其中，ΔP为储层生产压差，ΔPw为水区压差，ΔPg为气区压差，Kw为水相渗透率，Kg为气相渗透率，μw为水相粘度，μg为气相粘度，ρw为水相密度，ρg为气相密度，q为气体产量，Bg为气体体积系数，L为初始气水边界与所述目标井的距离，L1为水舌点移动距离，h为储层厚度，Lh为所述目标井的水平段长度，α为储层倾角，g为重力加速度。
4.一种边水气藏水平井的见水时间预测装置，其特征在于，包括：
特征参数确定模块，用于确定目标井的特征参数，所述特征参数包括水相渗透率、气相渗透率、水相粘度、气相粘度、水相密度、气相密度、储层孔隙度、含水饱和度、储层生产压差、气体产量、气体体积系数、初始气水边界与所述目标井的距离、储层厚度、所述目标井的水平段长度和储层倾角；
移动速度确定模块，用于根据所述特征参数确定水舌点移动速度；
见水时间确定模块，用于根据所述水舌点移动速度确定所述目标井的见水时间；
其中，所述根据所述特征参数确定水舌点移动速度，包括：
根据公式
确定水舌点移动速度；
所述根据所述水舌点移动速度确定所述目标井的见水时间，包括：
根据公式
确定目标井的见水时间；
其中，tbt(0,t)为目标井的见水时间，L1为水舌点移动距离，Viw为水舌点移动速度，Kw为水相渗透率，Kg为气相渗透率，μw为水相粘度，μg为气相粘度，ρw为水相密度，ρg为气相密度，φ为储层孔隙度，Sw为含水饱和度，ΔP为储层生产压差，q为气体产量，Bg为气体体积系数，L为初始气水边界与所述目标井的距离，h为储层厚度，Lh为所述目标井的水平段长度，α为储层倾角，g为重力加速度。
5.如权利要求4所述的边水气藏水平井的见水时间预测装置，其特征在于，在地层均质含倾角情况下，所述目标井的见水时间通过以下公式确定：
6.如权利要求4或5所述的边水气藏水平井的见水时间预测装置，其特征在于，所述储层生产压差通过以下公式确定：
其中，ΔP为储层生产压差，ΔPw为水区压差，ΔPg为气区压差，Kw为水相渗透率，Kg为气相渗透率，μw为水相粘度，μg为气相粘度，ρw为水相密度，ρg为气相密度，q为气体产量，Bg为气体体积系数，L为初始气水边界与所述目标井的距离，L1为水舌点移动距离，h为储层厚度，Lh为所述目标井的水平段长度，α为储层倾角，g为重力加速度。
7.一种边水气藏水平井的见水时间预测装置，包括存储器、处理器、以及存储在所述存储器上的计算机程序，其特征在于，所述计算机程序被所述处理器运行时执行如下步骤：
确定目标井的特征参数，所述特征参数包括水相渗透率、气相渗透率、水相粘度、气相粘度、水相密度、气相密度、储层孔隙度、含水饱和度、储层生产压差、气体产量、气体体积系数、初始气水边界与所述目标井的距离、储层厚度、所述目标井的水平段长度和储层倾角；
根据所述特征参数确定水舌点移动速度；
根据所述水舌点移动速度确定所述目标井的见水时间；
其中，所述根据所述特征参数确定水舌点移动速度，包括：
根据公式
确定水舌点移动速度；
所述根据所述水舌点移动速度确定所述目标井的见水时间，包括：
根据公式
确定目标井的见水时间；
其中，tbt(0,t)为目标井的见水时间，L1为水舌点移动距离，Viw为水舌点移动速度，Kw为水相渗透率，Kg为气相渗透率，μw为水相粘度，μg为气相粘度，ρw为水相密度，ρg为气相密度，φ为储层孔隙度，Sw为含水饱和度，ΔP为储层生产压差，q为气体产量，Bg为气体体积系数，L为初始气水边界与所述目标井的距离，h为储层厚度，Lh为所述目标井的水平段长度，α为储层倾角，g为重力加速度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710998036.7A CN107977732B (zh) | 2017-10-24 | 2017-10-24 | 边水气藏水平井的见水时间预测方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710998036.7A CN107977732B (zh) | 2017-10-24 | 2017-10-24 | 边水气藏水平井的见水时间预测方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107977732A CN107977732A (zh) | 2018-05-01 |
| CN107977732B true CN107977732B (zh) | 2021-11-02 |
Family
ID=62012590
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710998036.7A Active CN107977732B (zh) | 2017-10-24 | 2017-10-24 | 边水气藏水平井的见水时间预测方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107977732B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108830410A (zh) * | 2018-06-05 | 2018-11-16 | 中国石油天然气集团有限公司 | 低渗透底水气藏的见水时间预测方法及装置 |
| CN112836438B (zh) * | 2019-11-25 | 2022-08-30 | 中国石油化工股份有限公司 | 一种非均质边水气藏气井见水时间预测方法及装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105631529A (zh) * | 2015-10-28 | 2016-06-01 | 中国石油化工股份有限公司 | 一种边水气藏见水时间预测方法 |
| WO2017139271A2 (en) * | 2016-02-08 | 2017-08-17 | Rs Energy Group Topco, Inc. | Method for estimating oil/gas production using statistical learning models |
-
2017
- 2017-10-24 CN CN201710998036.7A patent/CN107977732B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105631529A (zh) * | 2015-10-28 | 2016-06-01 | 中国石油化工股份有限公司 | 一种边水气藏见水时间预测方法 |
| WO2017139271A2 (en) * | 2016-02-08 | 2017-08-17 | Rs Energy Group Topco, Inc. | Method for estimating oil/gas production using statistical learning models |
Non-Patent Citations (2)
| Title |
|---|
| A new algorithm to quantify water coning in hydraulically fractured tight oil reservoirs;Shams Kalam;Sami Alnuaim;Muzammil Hussain Rammay;《2015 International Field Exploration and Development Conference (IFEDC 2015)》;IET;20160328;全文 * |
| 考虑储层倾角和水侵的边水气藏见水时间预测研究;李元生; 杨志兴; 藤赛男; 廖恒杰; 陈自立;《石油钻探技术》;20170125;91-96 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107977732A (zh) | 2018-05-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106991223B (zh) | 确定多层油藏水驱波及系数的方法及装置 | |
| Balmforth et al. | Shallow viscoplastic flow on an inclined plane | |
| CN109034447B (zh) | 一种水驱油藏含水上升率变化的预测方法及装置 | |
| CN105569641B (zh) | 一种改善近井压降的水平井瞬态压力快速计算模型建立方法 | |
| CN107451311B (zh) | 一种油水相对渗透率曲线计算方法及装置 | |
| CN107977732B (zh) | 边水气藏水平井的见水时间预测方法及装置 | |
| Christie et al. | Fast procedure for upscaling compositional simulation | |
| CN107630684B (zh) | 连续油管气举排水采气动态过程模拟方法及装置 | |
| Sinha et al. | Effective rheology of bubbles moving in a capillary tube | |
| Ramon et al. | On the hydrodynamic interaction between a particle and a permeable surface | |
| Zhao et al. | Spontaneous imbibition in a square tube with corner films: Theoretical model and numerical simulation | |
| CN108681793B (zh) | 深层油藏采油指数预测方法及装置 | |
| CN106960264B (zh) | 低渗透砂砾岩地层压力预测方法及装置 | |
| Todd | Unsteady flow in porous media by means of a Hele‐Shaw viscous fluid model | |
| CN104360412A (zh) | 致密深盆气成藏预测方法和装置 | |
| Beresnev et al. | Theory of breakup of core fluids surrounded by a wetting annulus in sinusoidally constricted capillary channels | |
| CN108952676A (zh) | 一种页岩气藏非均质性评价方法及其装置 | |
| CN112796745A (zh) | 一种蒸汽驱油藏生产动态数据的确定方法和装置 | |
| CN107975364B (zh) | 一种确定气井极限排砂时机的方法及装置 | |
| CN114021494A (zh) | 基于页岩纳米孔中两相流体渗流阻力的开发方案确定方法 | |
| Leung | A fast convolution method for implementing single-porosity finite/infinite aquifer models for water-influx calculations | |
| CN113392372B (zh) | 气井全生命周期无阻流量确定方法及系统 | |
| CN112112639B (zh) | 一种凝析气藏循环注气条件下地层压力确定方法及系统 | |
| CN110969273B (zh) | 生产井产量预测方法及装置 | |
| CN108920781B (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 | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 100007 No. 9 North Main Street, Dongcheng District, Beijing, Dongzhimen Patentee after: CHINA NATIONAL PETROLEUM Corp. Patentee after: CNPC ENGINEERING TECHNOLOGY R & D Co.,Ltd. Patentee after: CNPC JIANGHAN MACHINERY RESEARCH INSTITUTE Co.,Ltd. Address before: 100007 No. 9 North Main Street, Dongcheng District, Beijing, Dongzhimen Patentee before: China National Petroleum Corp. Patentee before: CNPC Drilling Engineering Technology Research Institute Co.,Ltd. Patentee before: Hubei PetroChina Jianghan Machinery Research Institute Co.,Ltd. Address after: 100007 No. 9 North Main Street, Dongcheng District, Beijing, Dongzhimen Patentee after: China National Petroleum Corp. Patentee after: CNPC Drilling Engineering Technology Research Institute Co.,Ltd. Patentee after: Hubei PetroChina Jianghan Machinery Research Institute Co.,Ltd. Address before: 100007 No. 9 North Main Street, Dongcheng District, Beijing, Dongzhimen Patentee before: China National Petroleum Corp. Patentee before: DRILLING RESEARCH INSTITUTE OF CHINA NATIONAL PETROLEUM Co. Patentee before: JIANGHAN MACHINERY Research Institute DR CNPC |
|
| CP01 | Change in the name or title of a patent holder | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20221229 Address after: 100007 No. 9 North Main Street, Dongcheng District, Beijing, Dongzhimen Patentee after: CHINA NATIONAL PETROLEUM Corp. Patentee after: CNPC ENGINEERING TECHNOLOGY R & D Co.,Ltd. Patentee after: CNPC JIANGHAN MACHINERY RESEARCH INSTITUTE Co.,Ltd. Patentee after: BEIJING PETROLEUM MACHINERY Co.,Ltd. Address before: 100007 No. 9 North Main Street, Dongcheng District, Beijing, Dongzhimen Patentee before: CHINA NATIONAL PETROLEUM Corp. Patentee before: CNPC ENGINEERING TECHNOLOGY R & D Co.,Ltd. Patentee before: CNPC JIANGHAN MACHINERY RESEARCH INSTITUTE Co.,Ltd. |
|
| TR01 | Transfer of patent right |