CN113946929A - 产水气井产能评价方法、装置、电子设备及介质 - Google Patents
产水气井产能评价方法、装置、电子设备及介质 Download PDFInfo
- Publication number
- CN113946929A CN113946929A CN202010681022.4A CN202010681022A CN113946929A CN 113946929 A CN113946929 A CN 113946929A CN 202010681022 A CN202010681022 A CN 202010681022A CN 113946929 A CN113946929 A CN 113946929A
- Authority
- CN
- China
- Prior art keywords
- gas well
- gas
- water
- productivity
- corrected
- 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.)
- Pending
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 143
- 238000011156 evaluation Methods 0.000 title claims description 8
- 230000035699 permeability Effects 0.000 claims abstract description 72
- 238000004519 manufacturing process Methods 0.000 claims abstract description 48
- 238000000034 method Methods 0.000 claims abstract description 37
- 239000007788 liquid Substances 0.000 claims abstract description 24
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 238000010606 normalization Methods 0.000 claims abstract description 17
- 239000012071 phase Substances 0.000 claims description 71
- 230000009545 invasion Effects 0.000 claims description 15
- 210000002615 epidermis Anatomy 0.000 claims description 14
- 230000006378 damage Effects 0.000 claims description 12
- 210000003491 skin Anatomy 0.000 claims description 11
- 238000003860 storage Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 8
- 230000015572 biosynthetic process Effects 0.000 claims description 7
- 239000000126 substance Substances 0.000 claims description 7
- 208000027418 Wounds and injury Diseases 0.000 claims description 6
- 239000008346 aqueous phase Substances 0.000 claims description 6
- 208000014674 injury Diseases 0.000 claims description 6
- 239000012530 fluid Substances 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 15
- 238000011161 development Methods 0.000 abstract description 5
- 235000020681 well water Nutrition 0.000 abstract description 3
- 239000002349 well water Substances 0.000 abstract description 3
- 230000008859 change Effects 0.000 description 14
- 238000012360 testing method Methods 0.000 description 8
- 239000002131 composite material Substances 0.000 description 6
- 239000007791 liquid phase Substances 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000012937 correction Methods 0.000 description 4
- 230000000704 physical effect Effects 0.000 description 4
- 102000003712 Complement factor B Human genes 0.000 description 3
- 108090000056 Complement factor B Proteins 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000005514 two-phase flow Effects 0.000 description 2
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Images
Classifications
-
- 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
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- Geometry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Separation Using Semi-Permeable Membranes (AREA)
Abstract
公开了一种产水气井产能评价方法、装置、电子设备及介质。该方法可以包括：建立气井产能方程；针对相渗曲线进行归一化处理，获得气相相对渗透率的回归关系式；根据气井生产液气比与回归关系式，获得WGR关系式；根据WGR关系式，校正气井产能方程中的层流系数与紊流系数；根据气井产能方程与校正后的层流系数、校正后的紊流系数，计算产水气井的无阻流量。本发明考虑了气井见水后紊流效应对气井产能的影响，对气井产能方程的紊流系数进行修正，修正后的产能方程计算的见水气井无阻流量考虑了紊流效应的影响，计算结果更加符合实际情况，能够为气田开发提供指导。
Description
技术领域
本发明涉及有水气藏勘探开发领域，更具体地，涉及一种产水气井产能评价方法、装置、电子设备及介质。
背景技术
对边、底水气藏而言，在气田开发的中后期普遍存在见水现象。气井产水对气藏的正常生产有很大的影响，储层中水相的存在可能使近井地带气相相对渗透率降低，进而降低气井的产能。气井见水后的产能评价对指导气田开发方案调整、提高剩余气动用程度和气田最终采收率具有至关重要的意义。
在前人的研究中，逐渐发展出了一些采用不同手段得到的产水气井产能评价方法，包括：
在印尼Arun气田的双对数压恢试井曲线中发现了液相伤害的复合模型特征，因此将水侵带液相伤害等效为近井表皮的复合模型，通过引入水相伤害表皮因子，将气藏作为伤害区-未伤害区复合模型来研究液相聚集对产能的影响。此后，部分国内学者也采用类似方法对气水两相复合模型进行了研究。
采用Hawkins表皮表征近井地带水侵带来的储层损害。综合来看，将水侵带液相伤害等效为近井表皮的复合模型这种方法的优点是计算简单，缺点则是计算精度较差，且无法考虑水侵区的气液两相流动的动态变化对气井产能的影响。
考虑水气比的拟压力产能方程方法，主要是采用气水两相拟压力研究气水两相气井产能。还可以同时考虑启动压力梯度、应力敏感和滑脱效应等影响。近年来又有学者将气、水的层流系数和紊流系数引入气水两相拟压力的定义式中，考虑气藏压力变化对气、水高压物性的影响。虽然气水两相拟压力方法计算结果较水侵带液相伤害复合模型更为准确，但该方法需要持续监测井底流压的变化，而现场生产条件往往无法满足，因而导致实际操作性低。
相对渗透率曲线修正法是将气井生产液气比WGR通过分流率的定义引入气相相对渗透率的表达式中，再将气相相对渗透率引入气井二项式产能方程中，分别校正见水气井二项式产能方程的层流系数A与紊流系数B。关于紊流系数B，前人的理论和实验研究表明：气体的紊流因子β对含水饱和度的变化非常敏感，随着含水饱和度的上升，紊流因子变大，当含水饱和度从40％上升到70％时，紊流因子增加8倍。部分国外学者通过实验研究得到气体渗透率和紊流因子的经验关系式，这些实验和关系式都表明当气体渗透率降低时，紊流因子会显著升高。还有实验表明含水饱和度的上升不但影响紊流因子，同时也会增加气相的密度和黏度。而在校正见水气井的产能方程系数时，有学者认为只有层流系数A与气相相对渗透率有关，而紊流系数B为定值，因此只对层流系数A进行了校正，得到了修正后的产水气井二项式产能方程，这一结果是基于紊流系数不受相渗和含水饱和度影响这一假设而得出的，因而没有考虑紊流系数的变化。
现有技术表明，相对渗透率曲线修正法将液气比动态变化与相渗曲线结合，有效利用了见水前期测试成果，在计算简便的基础上同时能考虑水侵区的气液两相流动的动态变化对气井产能的影响，但该方法目前尚未考虑气井见水后紊流效应对气井产能变化的影响。当液气比升高时，含水饱和度的增加不但影响到气体相对渗透率的变化(层流系数A的变化)，也会影响到气体的紊流因子和物性(层流系数B的变化)，而这种影响又是不可忽视的。因此，有必要开发一种产水气井产能评价方法、装置、电子设备及介质。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种产水气井产能评价方法、装置、电子设备及介质，其考虑了气井见水后紊流效应对气井产能的影响，对气井产能方程的紊流系数进行修正，修正后的产能方程计算的见水气井无阻流量考虑了紊流效应的影响，计算结果更加符合实际情况，能够为气田开发提供指导。
第一方面，本公开实施例提供了一种产水气井产能评价方法，包括：
建立气井产能方程；
针对相渗曲线进行归一化处理，获得气相相对渗透率的回归关系式；
根据气井生产液气比与所述回归关系式，获得WGR关系式；
根据所述WGR关系式，校正所述气井产能方程中的层流系数与紊流系数；
根据所述气井产能方程与校正后的层流系数、校正后的紊流系数，计算产水气井的无阻流量。
优选地，所述气井产能方程为：
其中，PR为地层压力，Pwf为井底流压，q为气井产量，A为无水采气期层流系数，B为无水采气期紊流系数。
优选地，气相相对渗透率的回归关系式为：
其中，a、b为回归系数，Krg(Sw)为见水后任意含水饱和度下的气相相对渗透率，Krw(Sw)为见水后任意含水饱和度下的水相相对渗透率。
优选地，气井生产液气比为：
其中，Bg为气体体积系数，μg为气相黏度，μw为水相黏度，WGR为气井生产液气比，Rcw凝析液气比。
优选地，WGR关系式为：
其中，
为关于WGR的气相相对渗透率。
优选地，通过公式(5)校正层流系数：
其中，A’为校正后的层流系数，Krg(Swi)为束缚水饱和度下的气相相对渗透率，St为气井见水后的总表皮，St＝S+Sdam，S为见水前表皮，Sdam为见水引起的附加表皮，
rdam为水侵带伤害半径。
优选地，通过公式(6)校正紊流系数：
其中，B’为校正后的紊流系数。
作为本公开实施例的一种具体实现方式，
第二方面，本公开实施例还提供了一种产水气井产能评价装置，包括：
方程建立模块，建立气井产能方程；
归一化模块，针对相渗曲线进行归一化处理，获得气相相对渗透率的回归关系式；
WGR计算模块，根据气井生产液气比与所述回归关系式，获得WGR关系式；
校正模块，根据所述WGR关系式，校正所述气井产能方程中的层流系数与紊流系数；
流量计算模块，根据所述气井产能方程与校正后的层流系数、校正后的紊流系数，计算产水气井的无阻流量。
优选地，所述气井产能方程为：
其中，PR为地层压力，Pwf为井底流压，q为气井产量，A为无水采气期层流系数，B为无水采气期紊流系数。
优选地，气相相对渗透率的回归关系式为：
其中，a、b为回归系数，Krg(Sw)为见水后任意含水饱和度下的气相相对渗透率，Krw(Sw)为见水后任意含水饱和度下的水相相对渗透率。
优选地，气井生产液气比为：
其中，Bg为气体体积系数，μg为气相黏度，μw为水相黏度，WGR为气井生产液气比，Rcw凝析液气比。
优选地，WGR关系式为：
其中，
为关于WGR的气相相对渗透率。
优选地，通过公式(5)校正层流系数：
其中，A’为校正后的层流系数，Krg(Swi)为束缚水饱和度下的气相相对渗透率，St为气井见水后的总表皮，St＝S+Sdam，S为见水前表皮，Sdam为见水引起的附加表皮，
rdam为水侵带伤害半径。
优选地，通过公式(6)校正紊流系数：
其中，B’为校正后的紊流系数。
第三方面，本公开实施例还提供了一种电子设备，该电子设备包括：
存储器，存储有可执行指令；
处理器，所述处理器运行所述存储器中的所述可执行指令，以实现所述的产水气井产能评价方法。
第四方面，本公开实施例还提供了一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的产水气井产能评价方法。
其有益效果在于：
(1)通过在气井二项式产能方程中引入气井水侵后气相相对渗透率对产能方程中紊流因子进行校正，考虑了紊流效应对产水气井产能的影响，进一步提高了产水气井产能预测精度；
(2)无需持续监测井底流压，条件简单，可操作性较高；
(3)实际计算在简单方便的同时能保持较高的计算精度，矿场应用效果较好。
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
图1示出了根据本发明的一个实施例的产水气井产能评价方法的步骤的流程图。
图2示出了根据本发明的一个实施例的气相相对渗透率随液气比的变化关系的示意图。
图3示出了根据本发明的一个实施例的不同水气比下P105-2井的IPR曲线的示意图。
图4示出了根据本发明的一个实施例的考虑与不考虑紊流系数时P105-2井的产能变化曲线的对比图。
图5示出了根据本发明的一个实施例的一种产水气井产能评价装置的框图。
附图标记说明：
201、方程建立模块；202、归一化模块；203、WGR计算模块；204、校正模块；205、流量计算模块。
具体实施方式
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。
气井见水后，近井地带气相相对渗透率急剧降低，气体紊流因子发生改变，从而产生气体紊流效应。紊流效应的产生会使气井产能进一步降低，紊流效应越强，气井产能降低幅度越大。本发明根据气水两相相渗曲线测试数据，将考虑气井水侵的相对渗透率引入到气井二项式产能方程中，分别对气井的层流系数A和紊流系数B进行校正，进而得到考虑紊流效应影响的产水气井二项式产能方程，计算产水气井无阻流量。该方法适用前提是气井见水前做过产能试井测试(已知见水前二项式产能方程)、压恢测试(已知见水前气井表皮和泄气半径)，且有产气层岩心相渗实验数据。
本发明提供一种产水气井产能评价方法，包括：
建立气井产能方程；在一个示例中，气井产能方程为：
其中，PR为地层压力，Pwf为井底流压，q为气井产量，A为无水采气期层流系数，B为无水采气期紊流系数。
具体地，建立气井产能方程为公式(1)。
针对相渗曲线进行归一化处理，获得气相相对渗透率的回归关系式；在一个示例中，气相相对渗透率的回归关系式为：
其中，a、b为回归系数，Krg(Sw)为见水后任意含水饱和度下的气相相对渗透率，Krw(Sw)为见水后任意含水饱和度下的水相相对渗透率。
具体地，将气井生产层气水两相相渗曲线进行归一化处理，得到任意含水饱和度下的气相相对渗透率表达式为公式(2)。
根据气井生产液气比与回归关系式，获得WGR关系式；在一个示例中，气井生产液气比为：
其中，Bg为气体体积系数，μg为气相黏度，μw为水相黏度，WGR为气井生产液气比，Rcw凝析液气比。
在一个示例中，WGR关系式为：
其中，
为关于WGR的气相相对渗透率。
具体地，气井生产液气比为公式(3)，将公式(3)带入公式(2)中，得到WGR关系式为为公式(4)。
根据WGR关系式，校正气井产能方程中的层流系数与紊流系数；在一个示例中，通过公式(5)校正层流系数：
其中，A’为校正后的层流系数，Krg(Swi)为束缚水饱和度下的气相相对渗透率，St为气井见水后的总表皮，St＝S+Sdam，S为见水前表皮，Sdam为见水引起的附加表皮，
rdam为水侵带伤害半径。
在一个示例中，通过公式(6)校正紊流系数：
其中，B’为校正后的紊流系数。
具体地，根据WGR关系式，通过公式(5)校正层流系数，通过公式(6)校正紊流系数。
根据气井产能方程与校正后的层流系数、校正后的紊流系数，计算产水气井的无阻流量。
具体地，将校正后的层流系数与紊流系数带入公式(1)中，计算产水气井的无阻流量。
本发明还提供一种产水气井产能评价装置，包括：
方程建立模块，建立气井产能方程；在一个示例中，气井产能方程为：
其中，PR为地层压力，Pwf为井底流压，q为气井产量，A为无水采气期层流系数，B为无水采气期紊流系数。
具体地，建立气井产能方程为公式(1)。
归一化模块，针对相渗曲线进行归一化处理，获得气相相对渗透率的回归关系式；在一个示例中，气相相对渗透率的回归关系式为：
其中，a、b为回归系数，Krg(Sw)为见水后任意含水饱和度下的气相相对渗透率，Krw(Sw)为见水后任意含水饱和度下的水相相对渗透率。
具体地，将气井生产层气水两相相渗曲线进行归一化处理，得到任意含水饱和度下的气相相对渗透率表达式为公式(2)。
WGR计算模块，根据气井生产液气比与回归关系式，获得WGR关系式；在一个示例中，气井生产液气比为：
其中，Bg为气体体积系数，μg为气相黏度，μw为水相黏度，WGR为气井生产液气比，Rcw凝析液气比。
在一个示例中，WGR关系式为：
其中，
为关于WGR的气相相对渗透率。
具体地，气井生产液气比为公式(3)，将公式(3)带入公式(2)中，得到WGR关系式为为公式(4)。
校正模块，根据WGR关系式，校正气井产能方程中的层流系数与紊流系数；在一个示例中，通过公式(5)校正层流系数：
其中，A’为校正后的层流系数，Krg(Swi)为束缚水饱和度下的气相相对渗透率，St为气井见水后的总表皮，St＝S+Sdam，S为见水前表皮，Sdam为见水引起的附加表皮，
rdam为水侵带伤害半径。
在一个示例中，通过公式(6)校正紊流系数：
其中，B’为校正后的紊流系数。
具体地，根据WGR关系式，通过公式(5)校正层流系数，通过公式(6)校正紊流系数。
流量计算模块，根据气井产能方程与校正后的层流系数、校正后的紊流系数，计算产水气井的无阻流量。
具体地，将校正后的层流系数与紊流系数带入公式(1)中，计算产水气井的无阻流量。
本发明还提供一种电子设备，电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述的产水气井产能评价方法。
本发明还提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述的产水气井产能评价方法。
为便于理解本发明实施例的方案及其效果，以下给出四个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
实施例1
图1示出了根据本发明的一个实施例的产水气井产能评价方法的步骤的流程图。
如图1所示，该产水气井产能评价方法包括：步骤101，建立气井产能方程；步骤102，针对相渗曲线进行归一化处理，获得气相相对渗透率的回归关系式；步骤103，根据气井生产液气比与回归关系式，获得WGR关系式；步骤104，根据WGR关系式，校正气井产能方程中的层流系数与紊流系数；步骤105，根据气井产能方程与校正后的层流系数、校正后的紊流系数，计算产水气井的无阻流量。
以一口实际生产气井A井为例，A井于2010年1月8日投产，油管内径0.0759m，初期放喷测试无阻流量326×104m3/天。2013年1月气井见水后液气比开始迅速上升，2014年10月关井，关井前水气比4m3/104m3。气井见水前在2012年7月实施了4个工作制度下的产能测试，得到A井二项式产能方程为
泄气半径646m，表皮系数-1.91。相渗曲线归一化处理后得到的关系式为
气井的高压物性参数如表1所示。
表1
根据气井的高压物性参数，计算不同水气比、不同压力下的气相相对渗透率和水侵带附加表皮，计算结果如表2所示。
表2
图2示出了根据本发明的一个实施例的气相相对渗透率随液气比的变化关系的示意图。
图3示出了根据本发明的一个实施例的不同水气比下P105-2井的IPR曲线的示意图。
图4示出了根据本发明的一个实施例的考虑与不考虑紊流系数时P105-2井的产能变化曲线的对比图。
气相相对渗透率随液气比的变化关系如图2所示，利用公式(5)和公式(6)分别校正气井二项式产能方程的系数A和系数B为A’和B’，得到新的考虑紊流效应的见水气井产能方程，进而计算见水气井无阻流量结果如图3、图4和表3所示。
表3
结果表明考虑紊流因子B变化时，无阻流量比不考虑紊流系数B变化时减少50万方/天左右，差值随液气比增大而减小；随着液气比增加，若不考虑紊流系数B的变化，计算无阻流量误差不可忽略，为35％左右。
实施例2
图5示出了根据本发明的一个实施例的一种产水气井产能评价装置的框图。
如图5所示，该产水气井产能评价装置，包括：
方程建立模块201，建立气井产能方程；
归一化模块201，针对相渗曲线进行归一化处理，获得气相相对渗透率的回归关系式；
WGR计算模块203，根据气井生产液气比与回归关系式，获得WGR关系式；
校正模块204，根据WGR关系式，校正气井产能方程中的层流系数与紊流系数；
流量计算模块205，根据气井产能方程与校正后的层流系数、校正后的紊流系数，计算产水气井的无阻流量。
作为可选方案，气井产能方程为：
其中，PR为地层压力，Pwf为井底流压，q为气井产量，A为无水采气期层流系数，B为无水采气期紊流系数。
作为可选方案，气相相对渗透率的回归关系式为：
其中，a、b为回归系数，Krg(Sw)为见水后任意含水饱和度下的气相相对渗透率，Krw(Sw)为见水后任意含水饱和度下的水相相对渗透率。
作为可选方案，气井生产液气比为：
其中，Bg为气体体积系数，μg为气相黏度，μw为水相黏度，WGR为气井生产液气比，Rcw凝析液气比。
作为可选方案，WGR关系式为：
其中，
为关于WGR的气相相对渗透率。
作为可选方案，通过公式(5)校正层流系数：
其中，A’为校正后的层流系数，Krg(Swi)为束缚水饱和度下的气相相对渗透率，St为气井见水后的总表皮，St＝S+Sdam，S为见水前表皮，Sdam为见水引起的附加表皮，
rdam为水侵带伤害半径。
作为可选方案，通过公式(6)校正紊流系数：
其中，B’为校正后的紊流系数。
实施例3
本公开提供一种电子设备包括，该电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述产水气井产能评价方法。
根据本公开实施例的电子设备包括存储器和处理器。
该存储器用于存储非暂时性计算机可读指令。具体地，存储器可以包括一个或多个计算机程序产品，该计算机程序产品可以包括各种形式的计算机可读存储介质，例如易失性存储器和/或非易失性存储器。该易失性存储器例如可以包括随机存取存储器(RAM)和/或高速缓冲存储器(cache)等。该非易失性存储器例如可以包括只读存储器(ROM)、硬盘、闪存等。
该处理器可以是中央处理单元(CPU)或者具有数据处理能力和/或指令执行能力的其它形式的处理单元，并且可以控制电子设备中的其它组件以执行期望的功能。在本公开的一个实施例中，该处理器用于运行该存储器中存储的该计算机可读指令。
本领域技术人员应能理解，为了解决如何获得良好用户体验效果的技术问题，本实施例中也可以包括诸如通信总线、接口等公知的结构，这些公知的结构也应包含在本公开的保护范围之内。
有关本实施例的详细说明可以参考前述各实施例中的相应说明，在此不再赘述。
实施例4
本公开实施例提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的产水气井产能评价方法。
根据本公开实施例的计算机可读存储介质，其上存储有非暂时性计算机可读指令。当该非暂时性计算机可读指令由处理器运行时，执行前述的本公开各实施例方法的全部或部分步骤。
上述计算机可读存储介质包括但不限于：光存储介质(例如：CD－ROM和DVD)、磁光存储介质(例如：MO)、磁存储介质(例如：磁带或移动硬盘)、具有内置的可重写非易失性存储器的媒体(例如：存储卡)和具有内置ROM的媒体(例如：ROM盒)。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (10)
1.一种产水气井产能评价方法，其特征在于，包括：
建立气井产能方程；
针对相渗曲线进行归一化处理，获得气相相对渗透率的回归关系式；
根据气井生产液气比与所述回归关系式，获得WGR关系式；
根据所述WGR关系式，校正所述气井产能方程中的层流系数与紊流系数；
根据所述气井产能方程与校正后的层流系数、校正后的紊流系数，计算产水气井的无阻流量。
2.根据权利要求1所述的产水气井产能评价方法，其中，所述气井产能方程为：
其中，PR为地层压力，Pwf为井底流压，q为气井产量，A为无水采气期层流系数，B为无水采气期紊流系数。
3.根据权利要求1所述的产水气井产能评价方法，其中，气相相对渗透率的回归关系式为：
其中，a、b为回归系数，Krg(Sw)为见水后任意含水饱和度下的气相相对渗透率，Krw(Sw)为见水后任意含水饱和度下的水相相对渗透率。
4.根据权利要求3所述的产水气井产能评价方法，其中，气井生产液气比为：
其中，Bg为气体体积系数，μg为气相黏度，μw为水相黏度，WGR为气井生产液气比，Rcw凝析液气比。
5.根据权利要求4所述的产水气井产能评价方法，其中，WGR关系式为：
其中，
为关于WGR的气相相对渗透率。
6.根据权利要求5所述的产水气井产能评价方法，其中，通过公式(5)校正层流系数：
其中，A’为校正后的层流系数，Krg(Swi)为束缚水饱和度下的气相相对渗透率，St为气井见水后的总表皮，St＝S+Sdam，S为见水前表皮，Sdam为见水引起的附加表皮，
rdam为水侵带伤害半径。
7.根据权利要求5所述的产水气井产能评价方法，其中，通过公式(6)校正紊流系数：
其中，B’为校正后的紊流系数。
8.一种产水气井产能评价装置，其特征在于，包括：
方程建立模块，建立气井产能方程；
归一化模块，针对相渗曲线进行归一化处理，获得气相相对渗透率的回归关系式；
WGR计算模块，根据气井生产液气比与所述回归关系式，获得WGR关系式；
校正模块，根据所述WGR关系式，校正所述气井产能方程中的层流系数与紊流系数；
流量计算模块，根据所述气井产能方程与校正后的层流系数、校正后的紊流系数，计算产水气井的无阻流量。
9.一种电子设备，其特征在于，所述电子设备包括：
存储器，存储有可执行指令；
处理器，所述处理器运行所述存储器中的所述可执行指令，以实现权利要求1-7中任一项所述的产水气井产能评价方法。
10.一种计算机可读存储介质，其特征在于，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现权利要求1-7中任一项所述的产水气井产能评价方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010681022.4A CN113946929A (zh) | 2020-07-15 | 2020-07-15 | 产水气井产能评价方法、装置、电子设备及介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010681022.4A CN113946929A (zh) | 2020-07-15 | 2020-07-15 | 产水气井产能评价方法、装置、电子设备及介质 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113946929A true CN113946929A (zh) | 2022-01-18 |
Family
ID=79325865
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010681022.4A Pending CN113946929A (zh) | 2020-07-15 | 2020-07-15 | 产水气井产能评价方法、装置、电子设备及介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113946929A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114991724A (zh) * | 2022-06-17 | 2022-09-02 | 中海石油(中国)有限公司 | 一种致密气井产能预测方法及系统 |
-
2020
- 2020-07-15 CN CN202010681022.4A patent/CN113946929A/zh active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114991724A (zh) * | 2022-06-17 | 2022-09-02 | 中海石油(中国)有限公司 | 一种致密气井产能预测方法及系统 |
| CN114991724B (zh) * | 2022-06-17 | 2024-01-02 | 中海石油(中国)有限公司 | 一种致密气井产能预测方法及系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10344584B2 (en) | Systems and methods for transient-pressure testing of water injection wells to determine reservoir damages | |
| Can et al. | Simple tools for forecasting waterflood performance | |
| CN105893679B (zh) | 低产水平井续流修正试井解释方法 | |
| CN107526891B (zh) | 一种聚合物驱大孔道油藏试井分析方法 | |
| US10858936B2 (en) | Determining geologic formation permeability | |
| CN104632158B (zh) | 二氧化碳混相驱条件下油井受效阶段划分及判别方法 | |
| CN106250984A (zh) | 油井的油水关系模式的判断方法及装置 | |
| US9471730B2 (en) | Generalized inflow performance model for oil wells of any inclined angle and a computer-implemented method thereof | |
| CN113946929A (zh) | 产水气井产能评价方法、装置、电子设备及介质 | |
| Gong et al. | Modeling of liquid injectivity in surfactant-alternating-gas foam enhanced oil recovery | |
| Ding et al. | Near-wellbore formation damage effects on well performance: A comparison between underbalanced and overbalanced drilling | |
| CN112036097B (zh) | 一种水锁气井的产能计算方法 | |
| CN111577264A (zh) | 裂缝孔隙型油藏水平井产能预测方法及装置 | |
| Jaripatke et al. | Eagle ford completions optimization-an operator's approach | |
| Pandey et al. | Vertical growth of hydraulic fractures in layered formations | |
| CN107956472B (zh) | 裂缝性碳酸盐岩储集层完井方式确定方法、装置及介质 | |
| CN113622908B (zh) | 一种水侵气藏废弃地层压力的确定方法 | |
| Jinasena et al. | Modeling and analysis of fluid flow through a non-prismatic open channel with application to drilling | |
| CN111582532A (zh) | 应力敏感性油藏水平井产液能力预测方法及装置 | |
| CN112145162B (zh) | 水油体积比计算方法及系统 | |
| CN113417632B (zh) | 基于压裂后压力确定煤储层原始地层压力的方法及装置 | |
| Ji et al. | Methods for Modelling Static Fractures in Reservoir Simulation | |
| CN115828787B (zh) | 一种改进的mdh试井分析方法 | |
| CN113006752B (zh) | 注聚压力预测的方法和装置 | |
| US20230259670A1 (en) | Permeability modeling in a reservoir simulation model using dynamic pressure transient analysis |
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 |