CN108180013A - 一种边水气藏水侵方向判别方法 - Google Patents
一种边水气藏水侵方向判别方法 Download PDFInfo
- Publication number
- CN108180013A CN108180013A CN201711474029.3A CN201711474029A CN108180013A CN 108180013 A CN108180013 A CN 108180013A CN 201711474029 A CN201711474029 A CN 201711474029A CN 108180013 A CN108180013 A CN 108180013A
- Authority
- CN
- China
- Prior art keywords
- water
- invasion
- gas well
- target
- permeability
- 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.)
- Granted
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 268
- 230000009545 invasion Effects 0.000 title claims abstract description 127
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000035699 permeability Effects 0.000 claims abstract description 107
- 230000001186 cumulative effect Effects 0.000 claims description 22
- 235000020681 well water Nutrition 0.000 claims description 9
- 239000002349 well water Substances 0.000 claims description 9
- 230000001174 ascending effect Effects 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims 2
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 230000010354 integration Effects 0.000 abstract description 3
- 238000011084 recovery Methods 0.000 abstract description 2
- 238000004458 analytical method Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000008595 infiltration Effects 0.000 description 3
- 238000001764 infiltration Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000002596 correlated effect Effects 0.000 description 2
- 208000004130 Blepharoptosis Diseases 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 206010015995 Eyelid ptosis Diseases 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000000875 corresponding effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000008398 formation water Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000010297 mechanical methods and process Methods 0.000 description 1
- 230000005226 mechanical processes and functions Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 201000003004 ptosis Diseases 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000009738 saturating Methods 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
- 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)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Filling Or Discharging Of Gas Storage Vessels (AREA)
Abstract
本发明公开了一种边水气藏水侵方向判别方法，属于边水气藏开发技术领域。该方法包括：获取气藏区域内多个气井的参数值，参数包括平均渗透率，渗透率变异系数，避水高度，以及井底到气水边界的距离，之后对参数值进行归一化处理；根据经归一化处理后的参数值，按照预设规则获取目标气井的水侵风险系数；根据目标气井的水侵风险系数获取气藏区域内的水侵风险系数等值线分布图，并根据水侵风险系数等值线分布图获取目标气井的水侵方向。本发明所提供的边水气藏水侵方向判别方法，实现了水侵方向的定量预测，且综合考虑各个因素，更为全面、科学。同时适用于边水气藏开发的早期阶段，为合理制定开采方案提供理论支持。
Description
技术领域
本发明涉及边水气藏开发技术领域，特别涉及一种边水气藏水侵方向判别方法。
背景技术
气藏中直接与气层连通的地下水叫做气层水。如果气层厚度不大，或构造较陡时气充满圈闭的高部位，而且水环绕在气藏的周缘，这种水称为边水，该气藏则形成了边水气藏。在边水气藏的开发过程中，边水侵入会造成气井出水，不仅增加气藏的开发、开采难度，而且还会造成气井产能损失，降低气藏采收率，影响气藏开发效益。因此，需要对边水气藏的水侵动态进行判断，特别是对水侵方向进行判别。通过判断水侵方向可为边水气藏制定科学、合理的开发方案提供理论支持，提高气藏开发效益。
现有技术中提供了一种判别方法，具体通过综合分析构造、储层及气井所处位置可以实现边水气藏水侵方向的初步判断。
发明人在实现本发明的过程中，发现现有技术中至少存在以下问题：
现有技术所提供的方法忽略了储层非均质性以及射孔段对水侵方向的影响，且判定结果不够科学全面。
发明内容
为了解决现有技术中水侵方向判断方法不够科学、且适用性低的问题，本发明实施例提供了一种边水气藏水侵方向判别方法，该技术方案具体如下：
一种边水气藏水侵方向判别方法，包括：
获取气藏区域内多个气井的参数值，参数包括平均渗透率，渗透率变异系数，避水高度，以及井底到气水边界的距离，并对所述参数值进行归一化处理；
根据经所述归一化处理后的参数值，按照预设规则获取所述目标气井的水侵风险系数，并根据所述目标气井的水侵风险系数获取所述气藏区域内的水侵风险系数等值线分布图；
根据所述水侵风险系数等值线分布图获取所述目标气井的水侵方向。
可选地，获取平均渗透率值包括：在所述目标气井的深度方向上，选定n个采样点，相邻所述采样点间隔预设距离，其中n为大于或者等于1的整数；
获取每个所述采样点的渗透率值，并根据如下公式(1)获取所述平均渗透率值：
其中，Kx是所述目标气井的渗透率，10-3μm2；
m是所述目标气井中第m个采样点；
n是所述目标气井中所述采样点的数量，个；
Kx(m)是所述目标气井中第m个采样点处的渗透率，10-3μm2。
可选地，所述平均渗透率的归一化处理，包括：对所述目标气井的平均渗透率进行归一化处理：
其中，x是所述边水气藏中第x口气井；
K′x是所述目标气井的归一化处理后平均渗透率，无量纲；
i是所述边水气藏中所述气井的数量，口。
可选地，获取所述目标气井的所述渗透率变异系数包括：
根据所述渗透率值，对所述目标气井的n个所述采样点进行降序或者升序排列，并获取所述采样点排列后的序数；
获取目标采样点的渗透率累计百分比，获取所述目标采样点的序数累计百分比，以所述渗透率累计百分比为纵坐标，以所述序数累计百分比为横坐标建立坐标系，在所述坐标系中绘制所述渗透率累计百分比的洛伦兹曲线；
获取所述洛伦兹曲线与直线y＝x所围区域的面积，获取直线y＝x与坐标轴所围区域的面积，通过如下公式(3)获取所述渗透率变异系数：
其中，VKx是所述目标气井的渗透率变异系数，无量纲；
Sa是所述曲线与直线y＝x所围区域的面积，cm2；
Sb是直线y＝x与所述坐标轴所围区域的面积，cm2。
可选地，所述渗透率变异系数的归一化处理，包括：并根据如下公式(4)对所述目标气井的渗透率变异系数进行归一化处理：
其中，x是所述边水气藏中第x口气井；
VK′x是所述目标气井的归一化处理后渗透率变异系数，无量纲；
i是所述边水气藏中所述气井的数量，口。
可选地，获取所述避水高度，包括：
当射孔段底界位于气水界面之上时，所述避水高度为所述射孔段底界海拔与所述气水界面海拔的差值；
当所述射孔段底界海拔位于所述气水界面海拔之下时，所述避水高度的值为0。
可选地，所述避水高度的归一化处理，包括：根据如下公式(5)对所述目标气井的避水高度进行归一化处理：
其中，x是所述边水气藏中第x口气井；
Hx是所述目标气井的避水高度，m；
H′x是所述目标气井的归一化处理后避水高度，无量纲；
i是所述边水气藏中所述气井的数量，口。
可选地，所述井底到气水边界的距离的归一化处理，包括：
获取所述边水气藏区域内所有所述气井的井底到气水边界的距离，并根据如下公式(6)对所述井底到气水边界的距离进行归一化处理：
其中，x是所述边水气藏中第x口气井；
Lx是所述目标气井的井底到气水边界的距离，m；
L′x是所述目标气井的归一化处理后避水高度，无量纲；
i是所述边水气藏中所述气井的数量，口。
可选地，所述根据所述气井经过所述归一化处理后的参数值获取预设气井的水侵风险系数，包括：
获取所述水侵风险系数与所述参数的关系表达式；
根据所述关系表达式，以及经过所述归一化处理后的参数值获取所述目标气井的所述水侵风险系数。
可选地，获取所述水侵风险系数与所述参数的关系表达式，包括：
获取如下公式(7)：
Wx＝(-a×H′x-b×L′x+c×Vk′x+d×k′x)×100％ (7)
其中，Wx是目标气井水侵风险系数，无量纲；
H′x是所述目标气井经过归一化处理后的避水高度，无量纲；
L′x是所述目标气井经过归一化处理后的井底到气水边界的距离，无量纲；
Vk′x是所述目标气井经过归一化处理后的渗透率变异系数，无量纲；
k′x是所述目标气井经过归一化处理后的平均渗透率，无量纲；
a是H′x的权重系数，b是L＇x的权重系数，
c是Vk＇x的权重系数，d是k＇x的权重系数；
根据所述参数对水侵风险影响的重要程度，获取所述公式(7)中所述参数的权重系数。
可选地，所述获取所述公式(7)中所述参数的权重系数，包括：
以四个所述参数中任意一个为目标参数，对比所述目标参数与所有所述参数对所述水侵风向影响的重要程度，并确定所述目标参数的得分；
若所述目标参数对所述水侵风向的影响程度高，则所述目标参数得第一预设分值；
若所述目标参数与任一其他所述参数对所述水侵风向的影响程度相当，则所述目标参数得第二预设分值；
若所述目标参数对所述水侵风向的影响程度低，则所述目标参数得第三预设分值；
所述目标参数的得分为所述目标参数与每个所述参数对比后所得分数之和；
所述目标参数的所述权重系数通过如下公式(8)获取：
其中，A是所述目标参数的权重系数，无量纲；
Aw是所述目标参数的得分，无量纲。
本发明实施例提供的技术方案带来的有益效果是：
在本发明实施例所提供的水侵方向判别方法中，通过综合考虑气井的平均渗透率、渗透率变异系数、避水高度以及井底到气水边界的距离四大参数对水侵方向的影响。同时，本方法对上述四个参数的取值进行归一化处理，克服了不同参数间的数量级差，能够更佳科学地反馈四个参数对水侵方向的影响。此外，采用水侵风险系数等值线分布图来判别水侵方向，其效果直观、分析准确。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的边水气藏水侵方向判断流程示意图；
图2是本发明实施例提供的获取储层渗透率变异系数的洛伦兹曲线图；
图3是本发明实施例提供的某边水气藏水侵风险系数等值线平面分布图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
本发明实施例提供了一种边水气藏水侵方向判别方法，该方法具体包括：
S1、获取气藏区域内多个气井的参数值，这些参数包括平均渗透率，渗透率变异系数，避水高度，以及井底到气水边界的距离，并对参数值进行归一化处理；
S2、根据经归一化处理后的参数值，按照预设规则获取目标气井的水侵风险系数；
S3、根据目标气井的水侵风险系数获取气藏区域内的水侵风险系数等值线分布图，并根据水侵风险系数等值线分布图获取目标气井的水侵方向。
其中，气井的渗透率对水侵方向的影响显著，因此考量气井渗透率对水侵方向的影响非常必要。且由于在不同深度方向上，气井渗透率差异较大，因此在本发明实施例所提供的方法中，采用气井不同深度处渗透率的算数平均值，即平均渗透率，作为考量因素。
气藏区域内储层的非均质性对于水侵方向也具有显著影响，非均质性越强，优先出现水侵。在本发明实施例所提供的方法中，采用渗透率变异系数来表征储层非均质性，使得通过该判别方法得到的判别结果更为科学、客观。
避水高度体现的是气井射孔段底界海拔与气水界面海拔的高度差。将避水高度作为本判别方法的考量参数之一，体现出了射孔段对水侵方向的影响。
此外，还需说明的是，在本发明实施例中，井底到气水边界的距离是指，在边水气藏平面图上气井井底到气水边界的最短距离。
本发明实施例对于上述参数值的获取方法不做具体限定，例如可通过直接测量获取，或者借助相关软件计算获取。
且在获取水侵风向系数等值线分布图后，可根据水侵风险系数等值线的走向趋势来判定目标气井的水侵方向。
在本发明实施例所提供的水侵方向判别方法中，通过综合考虑气井的平均渗透率、渗透率变异系数、避水高度以及井底到气水边界的距离四大参数对水侵方向的影响。同时，本方法对上述四个参数的取值进行归一化处理，克服了不同参数间的数量级差，能够更加科学地反馈四个参数对水侵方向的影响。此外，采用水侵风险系数等值线分布图来判别水侵方向，其效果直观、分析准确。
以下将对步骤S1的具体过程进行阐述。
第一方面，在步骤S1中，获取目标气井的平均渗透率时，需要取气井不同深度处渗透率的算数平均值，具体按照以下步骤进行：
S111、在目标气井的深度方向上，选定n个采样点，相邻采样点间隔预设距离，其中n为大于或者等于1的整数；
S112、获取每个采样点的渗透率值，并根据如下公式(1)获取平均渗透率值：
其中，Kx是目标气井的平均渗透率，10-3μm2；
m是目标气井中第m个采样点；
n是目标气井中采样点的数量，个；
Kx(m)是目标气井中第m个采样点处的渗透率，10-3μm2。
通常在确定采样点时，在气井的深度方向上间距0.125m取一个采样点，对于采样点的数量n本发明实施例不做具体限定，例如400个、500个、600个等。
进一步地，按照如下方法对目标气井的平均渗透率进行归一化处理。具体包括：根据如下公式(2)对目标气井的平均渗透率进行归一化处理：
其中，x是边水气藏中第x口气井；
K＇x是目标气井的归一化处理后平均渗透率，无量纲；
i是边水气藏中气井的数量，口。
第二方面，在步骤S1中，获取目标气井的渗透率变异系数时，采用洛伦兹曲线法，具体按照以下步骤进行：
S121、根据渗透率值的大小，对目标气井的n个采样点进行降序或升序排列，并获取采样点排列后的序数。
S122、获取目标采样点的渗透率累计百分比，获取目标采样点的序数累计百分比，以渗透率累计百分比为纵坐标，以序数累计百分比为横坐标建立坐标系，在坐标系中绘制渗透率累计百分比的洛伦兹曲线；
在步骤S122中需要对采样点的渗透率累计百分比，以及序数累计百分比进行解释说明：
其中渗透率累计百分比为：按照步骤S121中的排序，第y个采样点的渗透率累计百分比为第1个至第y个采样点的渗透率之和除以所有采样点的渗透率之总和。序数累计百分比为第1个至第y个采样点的序数之和与所有采样点序数的总和的比值。示例地，参见表1，其中渗透率按照降序排列。
表1采样点渗透率取值表
其中，假设总采样点为6个，以第5个采样点为例，其渗透率累计百分比的取值等于第1个至第5个采样点的渗透率之和除以6个采样点渗透率之总和，即其渗透率累计百分比的取值等于136.668/160.189，为85.32％。
其序数累计百分比的取值等于15/21，为71.43％。
此处还需说明，由于在步骤S122中，需以渗透率累计百分比为纵坐标，以序数累计百分比为横坐标来绘制洛伦兹曲线，因此在取值时，增加序数0，其对应的渗透率为0％。
S123、获取洛伦兹曲线与直线y＝x所围区域的面积，获取直线y＝x与坐标轴所围区域的面积，通过如下公式(3)获取渗透率变异系数：
其中，VKx是目标气井的渗透率变异系数，无量纲；
Sa是曲线与直线y＝x所围区域的面积，cm2；
Sb是直线y＝x与坐标轴所围区域的面积，cm2。
更直观地，以图2为例，其中Sa位于y＝x的上方，Sb位于y＝x的下方。
进一步地，按照如下方法对目标气井的渗透率变异系数进行归一化处理，具体包括：根据如下公式(4)对目标气井的渗透率变异系数进行归一化处理：
其中，x是边水气藏中第x口气井；
VK′x是目标气井的归一化处理后渗透率变异系数，无量纲；
i是边水气藏中气井的数量，口。
通过获取气井的渗透率变异系数，体现了边水气藏区域的储层非均质性对水侵风险的影响，进而为科学判别气井水侵方向提供更有力支撑。
第三方面，在步骤S1中，获取目标气井的避水高度具体通过如下步骤进行：
步骤S131、当射孔段底界位于气水界面之上时，避水高度为射孔段底界海拔与气水界面海拔的差值；
步骤S132、当射孔段底界海拔位于气水界面海拔之下时，避水高度的值为0。
进一步地，按照如下方法对目标气井的避水高度进行归一化处理。具体包括：根据如下公式(5)对目标气井的避水高度进行归一化处理：
其中，x是边水气藏中第x口气井；
Hx是目标气井的避水高度，m；
H′x是目标气井的归一化处理后避水高度，无量纲；
i是边水气藏中气井的数量，口。
第四方面，在步骤S1中，对目标气井的井底到气水边界的距离的归一化处理按照如下步骤进行：
根据如下公式(6)对井底到气水边界的距离进行归一化处理：
其中，x是边水气藏中第x口气井；
Lx是目标气井的井底到气水边界的距离，m；
L′x是目标气井的归一化处理后避水高度，无量纲；
i是边水气藏中气井的数量，口。
上述对步骤S1中的技术方案作出了详细说明，下面将具体阐述步骤S2。
在步骤S2中，首先，获取水侵风险系数与参数的关系表达式。之后，根据关系表达式，以及经过归一化处理后的参数值获取目标气井的水侵风险系数。
其中，在获取水侵风险系数与参数的关系表达式时，具体按照如下步骤进行：
步骤S21、获取如下公式(7)：
Wx＝(-a×H′x-b×L′x+c×Vk′x+d×k′x)×100％ (7)
其中，Wx是目标气井水侵风险系数，无量纲；
H′x是目标气井经过归一化处理后的避水高度，无量纲；
L′x是目标气井经过归一化处理后的井底到气水边界的距离，无量纲；
Vk′x是目标气井经过归一化处理后的渗透率变异系数，无量纲；
k＇x是目标气井经过归一化处理后的平均渗透率，无量纲；
a是H＇x的权重系数，b是L＇x的权重系数，
c是Vk＇x的权重系数，d是k＇x的权重系数。
其中，目标气井的水侵风险系数越大，目标气井出现水侵的可能性就越大。并且，避水高度、井底到气水边界的距离与水侵风险系数成负相关，渗透率变异系数、平均渗透率与水侵风险系数成正相关。
步骤S22、根据参数对水侵风险影响的重要程度，获取公式(7)中参数的权重系数。
此处需要说明的是，不同地质类型的气藏，上述四个参数对水侵风险的影响重要程度顺序不同。举例来说，在孔隙型储层的气藏中，四个参数对水侵风险的影响重要程度排序为：避水高度＝井底到气水边界的距离＞渗透率变异系数＞平均渗透率；在裂缝-孔隙型储层的气藏中，四个参数对水侵风险的影响重要程度排序为：渗透率变异系数＞平均渗透率＞避水高度＝井底到气水边界的距离。
其中，在步骤S22中，公式(7)中参数的权重系数根据以下步骤确定。
步骤S221、以四个参数中任意一个为目标参数，对比目标参数与所有参数对水侵风向影响的重要程度，并确定目标参数的得分。
其中，若目标参数对水侵风险的影响程度高，则目标参数得第一预设分值。若目标参数与任一其他参数对水侵风险的影响程度相当，则目标参数得第二预设分值；此时包括了目标参数与其自身相对比的情况。若目标参数对水侵风险的影响程度低，则目标参数得第三预设分值。且目标参数的得分为该目标参数与每个参数对比后所得分数之和。
具体地，对于第一预设分值、第二预设分值以及第三预设分值的取值不做具体限定。示例地，与其他目标参数相比，当目标参数对水侵风险程度高时，第一预设分值为1分；当对比的两个参数对水侵风险影响程度相当时，第二预设分值为0.5分；当目标参数对水侵风险程度低时，第三预设分值为0。其中，当目标参数与其自身进行对比时，取第二预设分值。
在实际操作过程中，可采用列表计算的方法获取目标参数的得分。
步骤S222、目标参数的权重系数通过如下公式(8)获取：
其中，A是目标参数的权重系数，无量纲；Aw是所述目标参数的得分，无量纲。
在获取各个参数的权重系数后，根据目标气井的经过归一化处理后的参数值，以及公式(7)能够获取目标气井的水侵风险系数。
在获取目标气井的水侵风险系数之后，继续进行步骤S3。具体根据边水气藏区域内气井的水侵风险系数来获取气藏区域内的水侵风险系数等值线分布图，并根据水侵风险系数等值线分布图获取目标气井的水侵方向。
在获取水侵风险系数等值线分布图时，可采用克里金插值算法生成，此时需要获取至少三个气井的水侵风险参数。之后根据水侵风险参数等值线的走向与排布，可以直接判别目标气井的水侵方向。
本发明实施例所提供的边水气藏水侵方向判别方法，实现了水侵方向的定量预测，且综合考虑各方因素，更为全面、科学。同时适用于气藏开发的早期阶段，为合理制定开采方案提供理论支持。
下面将结合具体实施例来论述本发明实施例所提供的水侵方向判别方法的科学有效性。
实施例
下面将利用本发明实施例提供的水侵方向判别方法来判断上述气藏区域内各个气井的水侵方向。
步骤S1、获取气藏区域内目标气井的参数值，参数包括平均渗透率，渗透率变异系数，避水高度，以及井底到气水边界的距离，之后对参数值进行归一化处理；
按照步骤S111～S112、以及步骤S121～S123获取气藏区域内各个气井的平均渗透率、渗透率变异系数、以及进行归一化处理后的平均渗透率和渗透率变异系数，具体表2所示。
表2单井储层渗透率平均值及渗透率变异系数计算数据表
按照步骤S131～S132，以及第四方面所述方法获取避水高度、井底到气水边界的距离，以及归一化处理后避水高度、井底到气水边界的距离的值，具体如表3所示。
表3单井避水高度及其井底距原始气水边界的距离计算数据表
步骤S2、根据经归一化处理后的参数值，按照预设规则获取目标气井的水侵风险系数。
步骤S21、确定如下公式：
Wx＝(-a×H′x-b×L′x+c×Vk′x+d×k′x)×100％
其中，Wx是目标气井水侵风险系数，无量纲；
H′x是目标气井经过归一化处理后的避水高度，无量纲；
L′x是目标气井经过归一化处理后的井底到气水边界的距离，无量纲；
Vk′x是目标气井经过归一化处理后的渗透率变异系数，无量纲；
k′x是目标气井经过归一化处理后的平均渗透率，无量纲；
a是H′x的权重系数，b是L′x的权重系数，
c是Vk′x的权重系数，d是k′x的权重系数。
该边水气藏储层为孔隙型储层，则四个参数：平均渗透率(Kx)、渗透率变异系数(Vkx)、避水高度(Hx)、井底到气水边界的距离(Lx)的重要性排序为Hx＝Lx>Vkx>Kx。
步骤S22、确定上述四个参数的权重系数a、b、c、d，具体如表4所示，其中，以纵向与横向相比确定得分。
表4水侵风险系数表达式中权重系数a、b、b、d的权重计算数据表
| 序号 | 比较指标 | Hx′ | Lx′ | VKx′ | Kx′ | 求和 | 权重系数 | 权重 |
| 1 | Hx′ | 0.5 | 0.5 | 1 | 1 | 3 | a | 0.375 |
| 2 | Lx′ | 0.5 | 0.5 | 1 | 1 | 3 | b | 0.375 |
| 3 | Vkx′ | 0 | 0 | 0.5 | 1 | 1.5 | c | 0.1875 |
| 4 | Kx′ | 0 | 0 | 0 | 0.5 | 0.5 | d | 0.0625 |
如此，确定a、b、c、d的取值分别为0.375、0.375、0.1875、0.0625，从而确定该边水气藏水侵风险系数Wx数学表达式：
Wx＝(-0.375×H′x-0.375×L′x+0.1875×Vk′x+0.0625×k′x)×100％
在确定了风险系数Wx的数学表达式之后，获取目标气井的水侵风险系数，如表5所示。
表5气藏水侵风险系数计算数据表
| 序号 | 井名 | Hx′ | Lx′ | Kx′ | Vkx′ | Wx |
| 1 | H5 | 0 | 0.0106 | 0.0262 | 0.0754 | 1.18 |
| 2 | H1-X1 | 0.0493 | 0.0223 | 0.0512 | 0.1029 | -0.44 |
| 3 | H4-X3 | 0.0532 | 0.0267 | 0.0083 | 0.0914 | -1.23 |
| 4 | H4-X1 | 0.0691 | 0.0411 | 0.1422 | 0.0994 | -1.38 |
| 5 | H1-X2 | 0.0617 | 0.0436 | 0.0042 | 0.0903 | -2.23 |
| 6 | H10 | 0.1096 | 0.1126 | 0.3965 | 0.0926 | -4.12 |
| 7 | H1 | 0.0909 | 0.0977 | 0.2475 | 0.0629 | -4.35 |
| 8 | H4-2 | 0.1342 | 0.0858 | 0.0804 | 0.0937 | -5.99 |
| 9 | H8 | 0.1309 | 0.1703 | 0.0144 | 0.0857 | -9.60 |
| 10 | H4 | 0.1471 | 0.1883 | 0.0269 | 0.1086 | -10.37 |
| 11 | H4-X4 | 0.1540 | 0.2010 | 0.0023 | 0.0971 | -11.47 |
步骤S3、根据目标气井的水侵风险系数获取气藏区域内的水侵风险系数等值线分布图，并根据水侵风险系数等值线分布图获取目标气井的水侵方向。
据表5中计算的11口单井水侵风险系数，采用克里金插值算法得到气藏水侵风险系数等值线平面分布图，如图3所示。
根据水侵风险系数的地质意义，地层水首先沿水侵风险系数较大的井区侵入。因此，由图3可判别，气藏北翼H4-X3井区地层水是由气藏东北侧侵入的，气藏南翼H1-X1井区地层水是由气藏南侧侵入的。
由气藏实际出水井情况看，气藏北翼低部位的H5井测试气水同产，水侵风险系数为+1.18％，气藏南翼边部的H1-X1井于2010年产出地层水，水侵风险系数为-0.44％，气藏北翼边部H4-X3井于2012年产出地层水，水侵风险系数为-1.23％，出水井水侵风险系数大小与气藏实际出水顺序是一致的。说明本发明实施例所提供的水侵方向判别方法切实可行。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (11)
1.一种边水气藏水侵方向判别方法，其特征在于，包括：
获取气藏区域内多个气井的参数值，参数包括平均渗透率，渗透率变异系数，避水高度，以及井底到气水边界的距离，之后对所述参数值进行归一化处理；
根据经所述归一化处理后的参数值，按照预设规则获取所述目标气井的水侵风险系数；
根据所述目标气井的水侵风险系数获取所述气藏区域内的水侵风险系数等值线分布图，并根据所述水侵风险系数等值线分布图获取所述目标气井的水侵方向。
2.根据权利要求1所述的判别方法，其特征在于，获取平均渗透率值包括：
在所述目标气井的深度方向上，选定n个采样点，相邻所述采样点间隔预设距离，其中n为大于或者等于1的整数；
获取每个所述采样点的渗透率，并根据如下公式(1)获取所述平均渗透率值：
其中，Kx是所述目标气井的平均渗透率，10-3μm2；
m是所述目标气井中第m个采样点；
n是所述目标气井中所述采样点的数量，个；
Kx(m)是所述目标气井中第m个采样点处的渗透率，10-3μm2。
3.根据权利要求2所述的判别方法，其特征在于，所述平均渗透率的归一化处理，包括：
根据如下公式(2)对所述目标气井的平均渗透率进行归一化处理：
其中，x是所述边水气藏中第x口气井；
K′x是所述目标气井的归一化处理后平均渗透率，无量纲；
i是所述边水气藏中所述气井的数量，口。
4.根据权利要求2所述的判别方法，其特征在于，获取目标气井的所述渗透率变异系数值，包括：
根据所述渗透率值，对所述目标气井的n个所述采样点进行降序或者升序排列，并获取所述采样点排列后的序数；
获取所述目标采样点的序数累计百分比，获取目标采样点的渗透率累计百分比，以所述渗透率累计百分比为纵坐标，以所述序数累计百分比为横坐标建立坐标系，在所述坐标系中绘制所述渗透率累计百分比的洛伦兹曲线；
获取所述洛伦兹曲线与直线y＝x所围区域的面积，获取直线y＝x与坐标轴所围区域的面积，通过如下公式(3)获取所述渗透率变异系数：
其中，VKx是所述目标气井的渗透率变异系数，无量纲；
Sa是所述曲线与直线y＝x所围区域的面积，cm2；
Sb是直线y＝x与所述坐标轴所围区域的面积，cm2。
5.根据权利要求4所述的判别方法，其特征在于，所述渗透率变异系数的归一化处理，包括：
根据如下公式(4)对所述目标气井的渗透率变异系数进行归一化处理：
其中，x是所述边水气藏中第x口气井；
VK′x是所述目标气井的归一化处理后渗透率变异系数，无量纲；
i是所述边水气藏中所述气井的数量，口。
6.根据权利要求1所述的判别方法，其特征在于，获取避水高度值，包括：
当射孔段底界位于气水界面之上时，所述避水高度为所述射孔段底界海拔与所述气水界面海拔的差值；
当所述射孔段底界海拔位于所述气水界面海拔之下时，所述避水高度的值为0。
7.根据权利要求6所述的判别方法，其特征在于，所述避水高度的归一化处理，包括：
根据如下公式(5)对所述目标气井的避水高度进行归一化处理：
其中，x是所述边水气藏中第x口气井；
Hx是所述目标气井的避水高度，m；
H′x是所述目标气井的归一化处理后避水高度，无量纲；
i是所述边水气藏中所述气井的数量，口。
8.根据权利要求1所述的判别方法，其特征在于，所述井底到气水边界的距离的归一化处理，包括：
获取所述边水气藏区域内所有所述气井的井底到气水边界的距离，并根据如下公式(6)对所述井底到气水边界的距离进行归一化处理：
其中，x是所述边水气藏中第x口气井；
Lx是所述目标气井的井底到气水边界的距离，m；
L′x是所述目标气井的归一化处理后避水高度，无量纲；
i是所述边水气藏中所述气井的数量，口。
9.根据权利要求1所述的判别方法，其特征在于，所述根据所述气井经过所述归一化处理后的参数值获取预设气井的水侵风险系数，包括：
获取所述水侵风险系数与所述参数的关系表达式；
根据所述关系表达式，以及经过所述归一化处理后的参数值获取所述目标气井的所述水侵风险系数。
10.根据权利要求9所述的判别方法，其特征在于，获取所述水侵风险系数与所述参数的关系表达式，包括：
获取如下公式(7)：
Wx＝(-a×H′x-b×L′x+c×Vk′x+d×k′x)×100％ (7)
其中，Wx是目标气井水侵风险系数，无量纲；
H′x是所述目标气井经过归一化处理后的避水高度，无量纲；
L′x是所述目标气井经过归一化处理后的井底到气水边界的距离，无量纲；
Vk′x是所述目标气井经过归一化处理后的渗透率变异系数，无量纲；
k′x是所述目标气井经过归一化处理后的平均渗透率，无量纲；
a是H′x的权重系数，b是L′x的权重系数，
c是Vk′x的权重系数，d是k′x的权重系数；
根据所述参数对水侵风险影响的重要程度，获取所述公式(7)中所述参数的权重系数。
11.根据权利要求10所述的判别方法，其特征在于，所述获取所述公式(7)中所述参数的权重系数，包括：
以四个所述参数中任意一个为目标参数，对比所述目标参数与所有所述参数对所述水侵方向影响的重要程度，并确定所述目标参数的得分；
若所述目标参数对所述水侵风向的影响程度高，则所述目标参数得第一预设分值；
若所述目标参数与任一其他所述参数对所述水侵风向的影响程度相当，则所述目标参数得第二预设分值；
若所述目标参数对所述水侵风向的影响程度低，则所述目标参数得第三预设分值；
所述目标参数的得分为所述目标参数与每个所述参数对比后所得分数之和；所述目标参数的所述权重系数通过如下公式(8)获取：
其中，A是所述目标参数的权重系数，无量纲；
Aw是所述目标参数的得分，无量纲。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711474029.3A CN108180013B (zh) | 2017-12-29 | 2017-12-29 | 一种边水气藏水侵方向判别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711474029.3A CN108180013B (zh) | 2017-12-29 | 2017-12-29 | 一种边水气藏水侵方向判别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108180013A true CN108180013A (zh) | 2018-06-19 |
| CN108180013B CN108180013B (zh) | 2021-04-30 |
Family
ID=62549085
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711474029.3A Active CN108180013B (zh) | 2017-12-29 | 2017-12-29 | 一种边水气藏水侵方向判别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108180013B (zh) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109102538A (zh) * | 2018-07-17 | 2018-12-28 | 成都信息工程大学 | 利用等值线数据提取槽脊特征点与槽脊特征线的方法、天气图分析方法 |
| CN111927411A (zh) * | 2020-09-24 | 2020-11-13 | 西南石油大学 | 一种智能化有水气藏水侵跟踪和预警方法 |
| CN112377178A (zh) * | 2020-11-24 | 2021-02-19 | 中国石油天然气股份有限公司 | 一种边水水侵前缘的定量诊断方法 |
| CN112709547A (zh) * | 2019-10-24 | 2021-04-27 | 中国石油化工股份有限公司 | 一种产水气井堵水时机的判别方法 |
| CN112983401A (zh) * | 2021-04-30 | 2021-06-18 | 西南石油大学 | 一种边水气藏水侵边界计算方法 |
| CN113268846A (zh) * | 2020-02-17 | 2021-08-17 | 中国石油天然气股份有限公司 | 确定有水气藏治水对策的方法、装置、终端和存储介质 |
| CN113705878A (zh) * | 2021-08-24 | 2021-11-26 | 中海石油(中国)有限公司深圳分公司 | 水平井出水量的确定方法、装置、计算机设备及存储介质 |
| CN115577213A (zh) * | 2022-11-11 | 2023-01-06 | 成都理工大学 | 一种水侵方向预测方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2405933C1 (ru) * | 2009-04-27 | 2010-12-10 | Игорь Анатольевич Чернобровкин | Способ исследования газовых и газоконденсатных скважин |
| CN102041995A (zh) * | 2010-12-02 | 2011-05-04 | 中国海洋石油总公司 | 复杂油藏水淹状况监测系统 |
| CN104100245A (zh) * | 2013-04-08 | 2014-10-15 | 中国石油化工股份有限公司 | 人工边水驱断块油藏筛选评价方法 |
| CN104459823A (zh) * | 2014-12-17 | 2015-03-25 | 中冶集团武汉勘察研究院有限公司 | 一种综合物探测试海水侵入界面方法 |
| CN105443120A (zh) * | 2015-10-26 | 2016-03-30 | 中国石油化工股份有限公司 | 一种海相整装气田边水水侵早期特征分析方法 |
| CN105631529A (zh) * | 2015-10-28 | 2016-06-01 | 中国石油化工股份有限公司 | 一种边水气藏见水时间预测方法 |
-
2017
- 2017-12-29 CN CN201711474029.3A patent/CN108180013B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2405933C1 (ru) * | 2009-04-27 | 2010-12-10 | Игорь Анатольевич Чернобровкин | Способ исследования газовых и газоконденсатных скважин |
| CN102041995A (zh) * | 2010-12-02 | 2011-05-04 | 中国海洋石油总公司 | 复杂油藏水淹状况监测系统 |
| CN104100245A (zh) * | 2013-04-08 | 2014-10-15 | 中国石油化工股份有限公司 | 人工边水驱断块油藏筛选评价方法 |
| CN104459823A (zh) * | 2014-12-17 | 2015-03-25 | 中冶集团武汉勘察研究院有限公司 | 一种综合物探测试海水侵入界面方法 |
| CN105443120A (zh) * | 2015-10-26 | 2016-03-30 | 中国石油化工股份有限公司 | 一种海相整装气田边水水侵早期特征分析方法 |
| CN105631529A (zh) * | 2015-10-28 | 2016-06-01 | 中国石油化工股份有限公司 | 一种边水气藏见水时间预测方法 |
Non-Patent Citations (1)
| Title |
|---|
| 冯曦等: "有效治理气藏开发过程中水侵影响的问题及认识", 《开发工程》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109102538B (zh) * | 2018-07-17 | 2020-09-29 | 成都信息工程大学 | 利用等值线数据提取槽脊特征点与槽脊特征线的方法、天气图分析方法 |
| CN109102538A (zh) * | 2018-07-17 | 2018-12-28 | 成都信息工程大学 | 利用等值线数据提取槽脊特征点与槽脊特征线的方法、天气图分析方法 |
| CN112709547A (zh) * | 2019-10-24 | 2021-04-27 | 中国石油化工股份有限公司 | 一种产水气井堵水时机的判别方法 |
| CN113268846A (zh) * | 2020-02-17 | 2021-08-17 | 中国石油天然气股份有限公司 | 确定有水气藏治水对策的方法、装置、终端和存储介质 |
| CN113268846B (zh) * | 2020-02-17 | 2022-11-04 | 中国石油天然气股份有限公司 | 确定有水气藏治水对策的方法、装置、终端和存储介质 |
| CN111927411A (zh) * | 2020-09-24 | 2020-11-13 | 西南石油大学 | 一种智能化有水气藏水侵跟踪和预警方法 |
| CN112377178A (zh) * | 2020-11-24 | 2021-02-19 | 中国石油天然气股份有限公司 | 一种边水水侵前缘的定量诊断方法 |
| CN112377178B (zh) * | 2020-11-24 | 2024-04-30 | 中国石油天然气股份有限公司 | 一种边水水侵前缘的定量诊断方法 |
| CN112983401A (zh) * | 2021-04-30 | 2021-06-18 | 西南石油大学 | 一种边水气藏水侵边界计算方法 |
| CN113705878A (zh) * | 2021-08-24 | 2021-11-26 | 中海石油(中国)有限公司深圳分公司 | 水平井出水量的确定方法、装置、计算机设备及存储介质 |
| CN113705878B (zh) * | 2021-08-24 | 2024-02-09 | 中海石油(中国)有限公司深圳分公司 | 水平井出水量的确定方法、装置、计算机设备及存储介质 |
| CN115577213A (zh) * | 2022-11-11 | 2023-01-06 | 成都理工大学 | 一种水侵方向预测方法 |
| CN115577213B (zh) * | 2022-11-11 | 2023-03-24 | 成都理工大学 | 一种水侵方向预测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108180013B (zh) | 2021-04-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108180013A (zh) | 一种边水气藏水侵方向判别方法 | |
| CN110363325B (zh) | 复杂断块油藏多目标注采调控优化方法 | |
| CN107203651A (zh) | 一种辫状河致密砂岩气藏储层构型解剖方法 | |
| CN104020509B (zh) | 基于Bayes判别分析的砂质滩坝沉积微相测井识别方法 | |
| CN107965318A (zh) | 一种火山岩油藏有效储层定量分类的方法 | |
| CN104732000B (zh) | 地下策略优化 | |
| CN111352162A (zh) | 页岩地层元素相的划分方法和元素相评价方法 | |
| CN108518219B (zh) | 一种基于天然能量开发油田规模提液选井的评价方法 | |
| CN103853922B (zh) | 一种油气田相似性的可视化定量评价方法 | |
| WO2019214262A1 (zh) | 确定控制陆源分散有机质分布三角洲体系外缘坡度的方法 | |
| CN103604916A (zh) | 一种连续型致密砂岩气藏含气范围的测量方法及系统 | |
| JP7137889B2 (ja) | 深層－超深層海相カーボナタイトのオイルプール賦存深度下限の予測方法 | |
| CN109298464B (zh) | 致密砂岩储层成岩相测井识别方法及装置 | |
| CN105888659A (zh) | 确定岩性油气藏成藏概率的方法及装置 | |
| CN109555518B (zh) | 一种基于聚类和判别算法的冲积扇单井构型识别方法 | |
| CN110991933A (zh) | 一种山地页岩气资源评价方法及系统 | |
| CN104965237A (zh) | 一种定量优选非常规油气有利区的方法 | |
| CN112610205B (zh) | 一种强非均质储层气井早期产能评价方法 | |
| CN117115376A (zh) | 一种识别与绘制gde范围的方法、系统、设备及介质 | |
| Li et al. | Quantification of river network types based on hierarchical structures | |
| CN112989257B (zh) | 海相页岩油气储层的产气量测定方法 | |
| Shirakawa et al. | Two new deer truffle species, Elaphomyces marmoratus and Elaphomyces fuscus spp. nov., from a secondary forest in Japan | |
| WO2018087988A1 (ja) | 金属鉱床の探鉱方法、資源開発方法、採鉱方法、二次硫化銅の産出方法、資源産出方法、鉱山開発方法およびボーリング方法 | |
| CN106845685B (zh) | 一种降低水平井试油压裂作业成本的参数优化方法 | |
| CN104699977B (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 |