CN111810119A - 一种高压碳酸盐岩有水气藏气井产能计算方法 - Google Patents

一种高压碳酸盐岩有水气藏气井产能计算方法 Download PDF

Info

Publication number
CN111810119A
CN111810119A CN202010705350.3A CN202010705350A CN111810119A CN 111810119 A CN111810119 A CN 111810119A CN 202010705350 A CN202010705350 A CN 202010705350A CN 111810119 A CN111810119 A CN 111810119A
Authority
CN
China
Prior art keywords
gas
water
well
pressure
coefficient
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
Application number
CN202010705350.3A
Other languages
English (en)
Other versions
CN111810119B (zh
Inventor
庞进
刘洪�
欧霖
卢灿洋
田园
卢宇
马誉畅
张旭
于希南
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University of Science and Technology
Original Assignee
Chongqing University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing University of Science and Technology filed Critical Chongqing University of Science and Technology
Priority to CN202010705350.3A priority Critical patent/CN111810119B/zh
Publication of CN111810119A publication Critical patent/CN111810119A/zh
Application granted granted Critical
Publication of CN111810119B publication Critical patent/CN111810119B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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/04—Measuring depth or liquid level
    • 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/08—Measuring diameters or related dimensions at the borehole
    • 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

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Geology (AREA)
  • Geophysics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Business, Economics & Management (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Marketing (AREA)
  • Agronomy & Crop Science (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Animal Husbandry (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明提供一种高压碳酸盐岩有水气藏气井产能计算方法,包括以下步骤，第一步：建立高压碳酸盐岩有水气藏气井产能模型；第二步：求解产能模型，求出达西系数和非达西系数的表达式；第三步：根据岩心渗透率应力敏感实验数据回归确定应力敏感系数b；第四步：根据测井原始饱和度和生产数据计算目前含水饱和度；第五步：回归气水相渗曲线，确定气水相渗曲线方程，根据不同储层的相对渗透率曲线方程确定气的相对渗透率Krg和水的相对渗透率Krw；第六步：计算达非达西系数A和达西系数B，确定产能方程表达式。本发明的有益效果是,采用本发明提供的技术方案可以准确计算高压有水气藏不同介质储层气井在应力敏感，不同产水方式下的气井产能。

Description

一种高压碳酸盐岩有水气藏气井产能计算方法
技术领域
本发明涉及油气开发技术领域，具体涉及一种高压碳酸盐岩有水气藏气井产能计算方法。
背景技术
随着我国深层海相碳酸盐岩气藏不断加大勘探开发力度，在这类气藏开发过程中，由于地层异常高压、普遍存在边底水等地质特征，给这类气藏的产能评价带来了挑战，常规的气井产能评价方法例如二项式、指数式等均为地层在正常压力系统，无应力敏感条件下，同时无边底水、无层内孔隙可动水条件下建立起来的，无论是从理论原理还是现场实施评价效果来看，已不再适用于高压有水气藏气井产能的评价。急需一种能够反映高压碳酸盐岩气藏应力敏感、边底水和孔隙可动水产出、不同类型储层介质的产能评价方法，指导这类气藏的高效开发。
发明内容
本发明要解决的技术问题是，提供一种能够反映高压碳酸盐岩气藏应力敏感、边底水和孔隙可动水产出、不同类型储层介质的产能评价方法，指导高压有水气藏的高效开发。
本发明提供一种高压碳酸盐岩有水气藏气井产能计算方法,包括以下步骤，
第一步：建立高压碳酸盐岩有水气藏气井产能模型；
第二步：求解产能模型，求出达西系数和非达西系数的表达式；
第三步：根据岩心渗透率应力敏感实验数据回归确定应力敏感系数b；
第四步：计算目前含水饱和度；
第五步：回归气水相渗曲线，确定气水相渗曲线方程，根据不同储层的相对渗透率曲线方程确定气的相对渗透率Krg和水的相对渗透率Krw；；
第六步：计算达非达西系数A和达西系数B，确定产能方程表达式。
进一步的，所述步骤第一步包括，
收集气井的储层厚度、井半径、日产气量、日产水量、地层压力、井底流压、气体密度、水密度数据，确定产能模型假设条件，建立高压有水气藏气井产能数学模型。
进一步的，包括以下步骤，
所述确定产能模型假设条件包括，
①水平均质等厚无限大圆形气水同层的储层，中心一ロ井；
②气水彼此不互溶；
③产层全部打开，流体径向流入井内；
④地层流体可压缩；
⑤流体粘度为常数，考虑气水两相高速非达西渗流，但不考虑启动压力梯度；
⑥忽略重力和毛管力的影响；
⑦流体为等温流动；
所述建立高压有水气藏气井产能数学模型包括，
①气相流动方程：
Figure BDA0002594525380000021
②水相流动方程：
Figure BDA0002594525380000022
③气水两相拟压力函数：
Figure BDA0002594525380000023
④储层应力敏感方程：
Figure BDA0002594525380000024
⑤边界条件：r＝rw,p＝pwf；r＝re,p＝pe；
式中,b为敏感系数,无因次；
pi是原始地层压力，单位为MPa；
ki是原始地层压力下对应的渗透率，无因次；
krw为水相相对渗透率、krg为气相相对渗透率，无因次；
pw为水相压力、pg为气相压力，单位为MPa；
vw为水相的速度、vg为气相的速度，单位为m/s；
μw为水相的粘度、μg为气相的粘度，单位为mPa·s；
βw为水相的速度系数、βg为气相的速度系数，
Figure BDA0002594525380000031
Figure BDA0002594525380000032
δ＝7.75×10-6；
ρw为水相密度、ρg为气相的密度，单位为kg/m3；
Kg为水相渗透率、kw为气相渗透率，单位为10-3μm2；
mg为气的质量流量、mw为水的质量流量、ml为气水和的质量流量，单位为kg/s；
ρsc为标准状态下气体的密度，单位为kg/m3；
qsc为标准状态下气体的体积流量，单位为m3/d；
α为水气质量比，kg/kg，根据生产数据得到。
h为油层厚度，单位为m，根据测井解释得到；
re为气藏半径，单位为m，根据不稳定生产曲线拟合得到；
rw为井眼半径，单位为m，钻井实录；
pwf为井底流压，单位为MPa，实测或者计算得到；
pe为地层压力，单位为MPa，实测或者计算得到。
进一步的，所述步骤第二步具体包括，
确定达非达西系数A和达西系数B的计算公式：
Figure BDA0002594525380000033
写为：
Figure BDA0002594525380000041
或
Figure BDA0002594525380000042
其中：
Figure BDA0002594525380000043
式中，A为产能方程达西系数，无因次；
B为产能方程非达西系数，无因次；
S为表皮系数，无因次；由试井解释得到。
Figure BDA0002594525380000044
为气体平均粘度，单位为mPa·s；
Figure BDA0002594525380000049
为气体平均偏差因子，无因次；
式中：
Figure BDA0002594525380000047
和
Figure BDA0002594525380000048
可以通过查天然气高压物性曲线得到。
进一步的，所述第三步具体包括，
根据储层类型的岩心渗透率应力敏感实验数据，利用应力敏感方程
Figure BDA0002594525380000045
确定应力敏感系数b。
进一步的，所述第四步具体包括，
若边底水未侵入气井，气井主要产出孔隙水，则
根据公式
Figure BDA0002594525380000046
计算Sw，
式中Swi为原始含水饱和度，无因次；
Sw为目前含水饱和度，无因次；
WGR为标状下水气体积比，单位为m3/104m3；
Bw为地层水体积系数，无因次；
式中，原始含水饱和度Swi可由测井解释得到；
如边底水已侵入气井，气井产水主要为边底水，则：
根据公式
Figure BDA0002594525380000051
计算Sw，
式中We为累积水侵量，单位为m3；
Wp为累积产水量，单位为m3；根据生产数据得到；
Gp为累积产气量，单位为m3；根据生产数据得到；
Bgi为原始地层压力下的天然气体积系数，无因次；
Bg为天然气体积系数，无因次；
式中，累积水侵量可由水驱气藏物质平衡方程计算得到。
进一步的，所述第五步具体包括，
回归气水相渗曲线，确定气水相渗曲线方程
Figure BDA0002594525380000052
Figure BDA0002594525380000053
将目前含水饱和度Sw代入气水相渗曲线方程，计算得到Krg和Krw。
进一步的，所述第六步具体包括，
将应力敏感系数b，相对渗透率Krg和Krw代入达非达西系数A和达西系数B的计算公式，计算A和B，确定气井的产能方程表达式，计算气井无阻流量。
本发明的有益效果是,采用本发明提供的技术方案可以准确计算高压有水气藏不同介质储层气井在应力敏感，不同产水方式下的气井产能。为高压碳酸盐岩有水气藏产能评价提供更科学的计算方法。
附图说明
图1为本发明一实施例流程图。
图2为本发明一实施例应力敏感实验数据回归曲线。
图3为本发明一实施例相对渗透率曲线。
具体实施方式
本发明的发明构思是，通过分析高压碳酸盐岩有水气藏地质和生产特征，建立高压碳酸盐岩有水气藏气井产能计算数学模型，并对模型进行近似求解。在该产能模型中考虑两个因素：一是不同储层介质(基质、裂缝、孔洞以及各种组合介质)的应力敏感性；二是气井两种可能的产水类型：边底水和孔隙水。分别通过岩心渗透率应力敏感测试数据回归分析和相对渗透率曲线的回归分析来描述上述两个因素，并将所描述的数学公式耦合到气井产能计算模型中，用于准确计算高压有水气藏不同介质储层气井在应力敏感，不同产水方式下的气井产能。为高压碳酸盐岩有水气藏产能评价提供更科学的计算方法。
实施例1
如图1所示，本发明提供一种高压碳酸盐岩有水气藏气井产能计算方法，包括以下七个操作步骤。
第一步：收集准备气井生产完井和生产数据，建立高压碳酸盐岩有水气藏气井产能模型。
收集气井的储层厚度、井半径、日产气量、日产水量、地层压力、井底流压、气体密度、水密度等基础数据，提出产能模型假设条件，建立高压有水气藏气井产能数学模型。
Ⅰ、模型假设条件
①水平均质等厚无限大圆形气水同层的储层，中心一ロ井；
②气水彼此不互溶；
③产层全部打开，流体径向流入井内；
④地层流体可压缩；
⑤流体粘度为常数，考虑气水两相高速非达西渗流，但不考虑启动压力梯度；
⑥忽略重力和毛管力的影响；
⑦流体为等温流动。
II、数学模型构成
①气相流动方程：
Figure BDA0002594525380000071
②水相流动方程：
Figure BDA0002594525380000072
③气水两相拟压力函数：
Figure BDA0002594525380000073
④储层应力敏感方程：
Figure BDA0002594525380000074
⑤边界条件：r＝rw,p＝pwf；r＝re,p＝pe；
式中krw、krg为水相、气相相对渗透率，无因次；
pw、pg为水相、气相压力，MPa；
vw、vg为水相、气相的速度，m/s；
μw、μg为水相、气相的粘度，mPa·s；
βw、βg为水相、气相的速度系数，
Figure BDA0002594525380000075
Figure BDA0002594525380000076
δ＝7.75×10-6；
ρw、ρg为水相、气相的密度，kg/m3；
Kg、kw为水相、气相渗透率，10-3μm2；
mg、mw、ml为气、水、气水和的质量流量，kg/s；
ρsc为标准状态下气体的密度，kg/m3；
qsc为标准状态下气体的体积流量，m3/d；
α为水气质量比，kg/kg，根据生产数据得到。
h为油层厚度，m，根据测井解释得到；
re为气藏半径，m，根据不稳定生产曲线拟合得到；
rw为井眼半径，m，钻井实录；
pwf为井底流压，MPa，实测或者计算得到；
pe为地层压力，MPa，实测或者计算得到；
第二步：采用近似方法求解产能模型，并引入表皮系数S，求出达西系数和非达西系数的表达式。
采用近似方法计算获取近似解，并考虑气井不完善性，引入表皮系数S，求解得到：
Figure BDA0002594525380000081
写为：
Figure BDA0002594525380000082
或
Figure BDA0002594525380000083
其中：
Figure BDA0002594525380000084
Figure BDA0002594525380000085
式中A为产能方程达西系数，无因次；
B为产能方程非达西系数，无因次；
S为表皮系数，无因次；由试井解释得到。
Figure BDA0002594525380000086
为气体平均粘度，mPa·s；
Figure BDA0002594525380000087
为气体平均偏差因子，无因次。
式中：
Figure BDA0002594525380000088
和Z可以通过查天然气高压物性曲线得到。
第三步：根据岩心渗透率应力敏感实验数据回归确定应力敏感系数b
根据不同储层类型的岩心渗透率应力敏感实验数据，利用应力敏感方程
Figure BDA0002594525380000089
(b为待定系数)确定应力敏感系数b。
第四步：根据气井产水情况计算目前含水饱和度。
根据气井实际产水情况，分为两种情况：
1、边底水未侵入气井，气井主要产出孔隙水时：
根据公式
Figure BDA0002594525380000091
计算Sw，
式中Swi为原始含水饱和度，无因次；
Sw为目前含水饱和度，无因次；
WGR为标状下水气体积比，m3/104m3；
Bw为地层水体积系数，无因次；
式中，原始含水饱和度Swi可由测井解释得到。
2、边底水已侵入气井，气井产水主要为边底水时：
根据公式
Figure BDA0002594525380000092
计算Sw，
式中We为累积水侵量，m3；
Wp为累积产水量，m3；根据生产数据得到。
Gp为累积产气量，m3；根据生产数据得到。
式中，累积水侵量可由水驱气藏物质平衡方程计算得到。
第五步：回归气水相渗曲线，确定气水相渗曲线方程，根据不同储层的相对渗透率曲线方程确定气和水的相对渗透率Krg和Krw
回归气水相渗曲线，确定气水相渗曲线方程
Figure BDA0002594525380000093
Figure BDA0002594525380000094
将目前含水饱和度Sw代入气水相渗曲线方程，计算得到目前含水饱和度下的相对渗透率Krg和Krw。
第六步：计算达非达西系数A和达西系数B，确定产能方程表达式
将应力敏感系数b，相对渗透率Krg和Krw代入达非达西系数A和达西系数B的计算公式，计算A和B，确定气井的产能方程表达式，计算气井无阻流量。
本发明第一步中，建立高压碳酸盐岩有水气藏气井产能模型的过程可以省去，直接利用本发明中给出的达非达西系数A和达西系数B计算，即在基础资料收集好后从第三步开始工作。
本发明适用于孔隙、裂缝、孔洞单一储层介质以及组合而成的多重介质，其应力敏感系数b和气水相对渗透率Krg和Krw取决于所选的储层介质类型。
本发明计算的气井产能方程可用于气井无阻流量计算、气井优化配产等。
实施例2
本实施例以四川某高压碳酸盐岩有水气藏气井为例对本发明的技术方案进行说明。
1、基础资料收集
该井地层厚度h＝49.6m，渗透率Ki＝2.2×10-3μm2，泄流半径re＝908m,表皮系数S＝0.35,原始含水饱和度Swi＝0.26；目前地层压力pe＝48.05MPa；目前井底流压pwf＝42.06MPa；目前地层压力天然气体积系数Bg＝0.0038,原始地层压力天然气体积系数Bgi＝0.0028天然气密度ρg＝187.7408Kg/m3,天然气偏差因子
Figure BDA0002594525380000101
天然气粘度
Figure BDA0002594525380000102
日产气量131998m3/d，日产水量63.28m3/d；累积产气量Gp＝8.59×108m3，累积产水量Wp＝2.47×104m3，累积水侵量We＝96.55×104m3；井半径0.1m，
2、应力敏感系数确定
如图2应力敏感实验数据回归曲线所示，根据该井所在工区的应力敏感实验数据回归和地质认识，该井储层为孔洞型，应力敏感系数b＝0.0289。
3、目前含水饱和度计算
该井边底水已侵入的情况，根据公式
Figure BDA0002594525380000103
计算目前含水饱和度Sw＝0.66
4、目前含水饱和度下的气水相对渗透率
如图3相对渗透率曲线所示
该井储层为孔洞型，回归相对渗透率曲线为：
Figure BDA0002594525380000104
Sw＝0.66代入其中得到目前含水饱和度下的气水相对渗透率分别为krg＝0.153、krw＝0.068。
5、将b＝0.0289、krg＝0.153、krw＝0.068以及其他生产和完井数据代入公式(1)和(2)中，计算达非达西系数A和达西系数B，分别为A＝2.17×10-8、B＝0.0181。该井的产能方程为：pe 2-pwf 2＝6.2×10-10qs 2 c+5.17×10-4qsc，气井无阻流量为155.7×104m3/d。
本发明的有益效果是，采用本发明提供的技术方案可以准确计算高压有水气藏不同介质储层气井在应力敏感，不同产水方式下的气井产能。为高压碳酸盐岩有水气藏产能评价提供更科学的计算方法。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (8)

1.一种高压碳酸盐岩有水气藏气井产能计算方法,其特征在于，包括以下步骤，
第一步：建立高压碳酸盐岩有水气藏气井产能模型；
第二步：求解产能模型，求出达西系数和非达西系数的表达式；
第三步：根据岩心渗透率应力敏感实验数据回归确定应力敏感系数b；
第四步：根据测井原始饱和度和生产数据计算目前含水饱和度；
第五步：回归气水相渗曲线，确定气水相渗曲线方程，根据不同储层的相对渗透率曲线方程确定气的相对渗透率Krg和水的相对渗透率Krw；
第六步：计算达非达西系数A和达西系数B，确定产能方程表达式。
2.如权利要求1所述的一种高压碳酸盐岩有水气藏气井产能计算方法,其特征在于，
所述步骤第一步包括：
收集气井的储层厚度、井半径、日产气量、日产水量、地层压力、井底流压、气体密度、水密度数据，确定产能模型假设条件，建立高压有水气藏气井产能数学模型。
3.如权利要求2所述的一种高压碳酸盐岩有水气藏气井产能计算方法,其特征在于，
所述确定产能模型假设条件包括：
①水平均质等厚无限大圆形气水同层的储层，中心一ロ井；
②气水彼此不互溶；
③产层全部打开，流体径向流入井内；
④地层流体可压缩；
⑤流体粘度为常数，考虑气水两相高速非达西渗流，但不考虑启动压力梯度；
⑥忽略重力和毛管力的影响；
⑦流体为等温流动；
所述建立高压有水气藏气井产能数学模型包括，
①气相流动方程：
Figure FDA0002594525370000021
②水相流动方程：
Figure FDA0002594525370000022
③气水两相拟压力函数：
Figure FDA0002594525370000023
④储层应力敏感方程：
Figure FDA0002594525370000024
⑤边界条件：r＝rw,p＝pwf；r＝re,p＝pe；
式中b为敏感系数,无因次；
pi是原始地层压力，单位为MPa；
ki是原始地层压力下对应的渗透率，无因次；
krw为水相相对渗透率、krg为气相相对渗透率，无因次；
pw为水相压力、pg为气相压力，单位为MPa；
vw为水相的速度、vg为气相的速度，单位为m/s；
μw为水相的粘度、μg为气相的粘度，单位为mPa·s；
βw为水相的速度系数、βg为气相的速度系数，
Figure FDA0002594525370000025
δ＝7.75×10-6；
ρw为水相密度、ρg为气相的密度，单位为kg/m3；
Kg为水相渗透率、kw为气相渗透率，单位为10-3μm2；
mg为气的质量流量、mw为水的质量流量、ml为气水和的质量流量，单位为kg/s；
ρsc为标准状态下气体的密度，单位为kg/m3；
qsc为标准状态下气体的体积流量，单位为m3/d；
α为水气质量比，kg/kg，根据生产数据得到。
h为油层厚度，单位为m，根据测井解释得到；
re为气藏半径，单位为m，根据不稳定生产曲线拟合得到；
rw为井眼半径，单位为m，钻井实录；
pwf为井底流压，单位为MPa，实测或者计算得到；
pe为地层压力，单位为MPa，实测或者计算得到。
4.如权利要求1所述的一种高压碳酸盐岩有水气藏气井产能计算方法,其特征在于，
所述步骤第二步具体包括，
确定达非达西系数A和达西系数B的计算公式：
Figure FDA0002594525370000031
写为：
Figure FDA0002594525370000032
或
Figure FDA0002594525370000033
其中：
Figure FDA0002594525370000034
Figure FDA0002594525370000035
式中A为产能方程达西系数，无因次；
B为产能方程非达西系数，无因次；
S为表皮系数，无因次；由试井解释得到。
Figure FDA0002594525370000036
为气体平均粘度，单位为mPa·s；
Figure FDA0002594525370000037
为气体平均偏差因子，无因次；
式中，
Figure FDA0002594525370000038
和
Figure FDA0002594525370000039
可以通过查天然气高压物性曲线得到。
5.如权利要求1所述的一种高压碳酸盐岩有水气藏气井产能计算方法,其特征在于，所述第三步具体包括：
根据储层类型的岩心渗透率应力敏感实验数据，利用应力敏感方程
Figure FDA0002594525370000041
确定应力敏感系数b。
6.如权利要求1所述的一种高压碳酸盐岩有水气藏气井产能计算方法,其特征在于，所述第四步具体包括，
若边底水未侵入气井，气井主要产出孔隙水，则
根据公式
Figure FDA0002594525370000042
计算Sw，
式中Swi为原始含水饱和度，无因次；
Sw为目前含水饱和度，无因次；
WGR为标状下水气体积比，单位为m3/104m3；
Bw为地层水体积系数，无因次；
式中，原始含水饱和度Swi可由测井解释得到；
如边底水已侵入气井，气井产水主要为边底水，则：
根据公式
Figure FDA0002594525370000043
计算Sw，
式中We为累积水侵量，单位为m3；
Wp为累积产水量，单位为m3；根据生产数据得到；
Gp为累积产气量，单位为m3；根据生产数据得到；
Bgi为原始地层压力下的天然气体积系数，无因次；
Bg为天然气体积系数，无因次；
式中，累积水侵量可由水驱气藏物质平衡方程计算得到。
7.如权利要求1所述的一种高压碳酸盐岩有水气藏气井产能计算方法,其特征在于，所述第五步具体包括，
回归气水相渗曲线，确定气水相渗曲线方程
Figure FDA0002594525370000051
mg、ng、mw、nw为回归待定系数，将目前含水饱和度Sw代入气水相渗曲线方程，计算得到Krg和Krw。
8.如权利要求1所述的一种高压碳酸盐岩有水气藏气井产能计算方法,其特征在于，所述第六步具体包括，
将应力敏感系数b，相对渗透率Krg和Krw代入达非达西系数A和达西系数B的计算公式，计算A和B，确定气井的产能方程表达式，计算气井无阻流量。
CN202010705350.3A 2020-07-21 2020-07-21 一种高压碳酸盐岩有水气藏气井产能计算方法 Active CN111810119B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010705350.3A CN111810119B (zh) 2020-07-21 2020-07-21 一种高压碳酸盐岩有水气藏气井产能计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010705350.3A CN111810119B (zh) 2020-07-21 2020-07-21 一种高压碳酸盐岩有水气藏气井产能计算方法

Publications (2)

Publication Number Publication Date
CN111810119A true CN111810119A (zh) 2020-10-23
CN111810119B CN111810119B (zh) 2022-06-28

Family

ID=72861780

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010705350.3A Active CN111810119B (zh) 2020-07-21 2020-07-21 一种高压碳酸盐岩有水气藏气井产能计算方法

Country Status (1)

Country Link
CN (1) CN111810119B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112360449A (zh) * 2020-12-11 2021-02-12 西南石油大学 一种应用于低渗透气藏和气井的可采气量预测方法
CN112392473A (zh) * 2020-11-17 2021-02-23 中国石油天然气股份有限公司 一种低渗透气藏型储气库气井注采能力评价方法
CN114756826A (zh) * 2022-06-16 2022-07-15 中海油田服务股份有限公司 气井产能预测方法、装置、电子设备和存储介质
CN115060875A (zh) * 2022-06-17 2022-09-16 大连理工大学 一种基于达西定律的水合物储层生产压力区间的确定方法
CN115853503A (zh) * 2022-12-16 2023-03-28 陕西延长石油（集团）有限责任公司 一种基于考虑储层应力敏感计算气井地层压力的方法
CN116245423A (zh) * 2023-03-31 2023-06-09 西南石油大学 一种考虑相渗影响的有水气藏开发效果评价方法
CN118008259A (zh) * 2024-02-05 2024-05-10 重庆科技大学 一种考虑应力敏感的气井产能试井解释方法
CN115060875B (zh) * 2022-06-17 2024-10-22 大连理工大学 一种基于达西定律的水合物储层生产压力区间的确定方法

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1486536A (en) * 1975-01-20 1977-09-21 Permeator Corp Method of maintaining and stabilizing the permeability of hydrocarbon reservoir rock
US20040060351A1 (en) * 2002-09-30 2004-04-01 Gunter William Daniel Process for predicting porosity and permeability of a coal bed
US20120158378A1 (en) * 2010-12-15 2012-06-21 Guillaume Enchery Method of developing a petroleum reservoir from a flow model calibrated through pressure and saturation map scaling
CN103046925A (zh) * 2012-12-31 2013-04-17 中国石油天然气股份有限公司 基于二项式的凝析气藏绝对无阻流量的获取方法及系统
CN103806904A (zh) * 2012-11-12 2014-05-21 中国石油天然气集团公司 砂岩储层非均匀损害深度的计算方法
CN104863560A (zh) * 2015-03-09 2015-08-26 东方宝麟科技发展(北京)有限公司 用于页岩气开采的宽网压裂方法
CN105243182A (zh) * 2015-09-08 2016-01-13 中国石油天然气股份有限公司 致密油压裂水平井的动态储量计算方法
CN106529184A (zh) * 2016-11-24 2017-03-22 重庆科技学院 倾斜有水气藏产水气井产能的计算方法
CN107301306A (zh) * 2017-08-04 2017-10-27 中国石油化工股份有限公司华北油气分公司 用于致密砂岩气藏压裂水平井的动态无阻流量预测方法
CN108133087A (zh) * 2017-12-11 2018-06-08 常州大学 一种气-水两相渗流应力敏感储层原始渗透率反演方法
CN108520143A (zh) * 2018-04-08 2018-09-11 中海石油(中国)有限公司 一种注气开发油藏气油比上升率表征方法
CN109025942A (zh) * 2018-08-09 2018-12-18 西南石油大学 一种致密气藏斜井压裂不规则多裂缝的产量计算方法
CN109488276A (zh) * 2019-01-16 2019-03-19 重庆科技学院 经水力压裂改造的产水页岩气井页岩气产量预测方法
CN110658100A (zh) * 2019-10-12 2020-01-07 重庆科技学院 一种气相阈压梯度的实验测试系统、方法和数据处理方法
CN110685651A (zh) * 2019-10-14 2020-01-14 重庆科技学院 一种多层合采气井产量劈分方法及系统

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1486536A (en) * 1975-01-20 1977-09-21 Permeator Corp Method of maintaining and stabilizing the permeability of hydrocarbon reservoir rock
US20040060351A1 (en) * 2002-09-30 2004-04-01 Gunter William Daniel Process for predicting porosity and permeability of a coal bed
US20120158378A1 (en) * 2010-12-15 2012-06-21 Guillaume Enchery Method of developing a petroleum reservoir from a flow model calibrated through pressure and saturation map scaling
CN103806904A (zh) * 2012-11-12 2014-05-21 中国石油天然气集团公司 砂岩储层非均匀损害深度的计算方法
CN103046925A (zh) * 2012-12-31 2013-04-17 中国石油天然气股份有限公司 基于二项式的凝析气藏绝对无阻流量的获取方法及系统
CN104863560A (zh) * 2015-03-09 2015-08-26 东方宝麟科技发展(北京)有限公司 用于页岩气开采的宽网压裂方法
CN105243182A (zh) * 2015-09-08 2016-01-13 中国石油天然气股份有限公司 致密油压裂水平井的动态储量计算方法
CN106529184A (zh) * 2016-11-24 2017-03-22 重庆科技学院 倾斜有水气藏产水气井产能的计算方法
CN107301306A (zh) * 2017-08-04 2017-10-27 中国石油化工股份有限公司华北油气分公司 用于致密砂岩气藏压裂水平井的动态无阻流量预测方法
CN108133087A (zh) * 2017-12-11 2018-06-08 常州大学 一种气-水两相渗流应力敏感储层原始渗透率反演方法
CN108520143A (zh) * 2018-04-08 2018-09-11 中海石油(中国)有限公司 一种注气开发油藏气油比上升率表征方法
CN109025942A (zh) * 2018-08-09 2018-12-18 西南石油大学 一种致密气藏斜井压裂不规则多裂缝的产量计算方法
CN109488276A (zh) * 2019-01-16 2019-03-19 重庆科技学院 经水力压裂改造的产水页岩气井页岩气产量预测方法
CN110658100A (zh) * 2019-10-12 2020-01-07 重庆科技学院 一种气相阈压梯度的实验测试系统、方法和数据处理方法
CN110685651A (zh) * 2019-10-14 2020-01-14 重庆科技学院 一种多层合采气井产量劈分方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
宋力等: "致密砂岩气藏压裂井两相流流入动态研究", 《天然气勘探与开发》 *
黄万书等: "压敏效应对低渗透砂岩气藏产能的影响", 《天然气技术》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112392473A (zh) * 2020-11-17 2021-02-23 中国石油天然气股份有限公司 一种低渗透气藏型储气库气井注采能力评价方法
CN112392473B (zh) * 2020-11-17 2023-11-28 中国石油天然气股份有限公司 一种低渗透气藏型储气库气井注采能力评价方法
CN112360449A (zh) * 2020-12-11 2021-02-12 西南石油大学 一种应用于低渗透气藏和气井的可采气量预测方法
CN114756826A (zh) * 2022-06-16 2022-07-15 中海油田服务股份有限公司 气井产能预测方法、装置、电子设备和存储介质
CN115060875A (zh) * 2022-06-17 2022-09-16 大连理工大学 一种基于达西定律的水合物储层生产压力区间的确定方法
CN115060875B (zh) * 2022-06-17 2024-10-22 大连理工大学 一种基于达西定律的水合物储层生产压力区间的确定方法
CN115853503A (zh) * 2022-12-16 2023-03-28 陕西延长石油（集团）有限责任公司 一种基于考虑储层应力敏感计算气井地层压力的方法
CN116245423A (zh) * 2023-03-31 2023-06-09 西南石油大学 一种考虑相渗影响的有水气藏开发效果评价方法
CN116245423B (zh) * 2023-03-31 2023-11-17 西南石油大学 一种考虑相渗影响的有水气藏开发效果评价方法
CN118008259A (zh) * 2024-02-05 2024-05-10 重庆科技大学 一种考虑应力敏感的气井产能试井解释方法

Also Published As

Publication number Publication date
CN111810119B (zh) 2022-06-28

Similar Documents

Publication Publication Date Title
CN111810119B (zh) 一种高压碳酸盐岩有水气藏气井产能计算方法
US11054540B2 (en) Computer implemented method for measurement of hydrocarbon content of tight gas reservoirs
CN106285646B (zh) 基于多信息融合的钻井漏失层位识别方法
CN104912550A (zh) 一种核磁共振测井资料定量计算储层产液剖面方法
US7849736B2 (en) Method for calculating the ratio of relative permeabilities of formation fluids and wettability of a formation downhole, and a formation testing tool to implement the same
US20240151139A1 (en) Systems and Methods for Identifying Two or More Charges into Reservoir Using Downhole Fluid Analysis
RU2601344C2 (ru) Устройство, система и способ для отбора проб пластовой среды
WO2006120366A1 (en) Methods for analysis of pressure response in underground formations
US10378349B2 (en) Methods of plotting advanced logging information
CN109555515A (zh) 地层坍塌压力确定方法和装置
CN110552694B (zh) 一种考虑多因素影响的泥质白云岩油藏油井产能评价方法
Oliveira et al. Deep Transient Testing Digital Products Create Novel Real-Time Reservoir Insights
KR101800796B1 (ko) 극저투수율암석의 투수율 측정 장치 및 그 방법
Liu et al. Simultaneous interpretation of three-phase relative permeability and capillary pressure for a tight carbonate reservoir from wireline formation testing
CN112036097B (zh) 一种水锁气井的产能计算方法
Badry et al. New wireline formation tester techniques and applications
CN109538199A (zh) 一种煤系地层含气量评价方法、装置及电子设备
US9013954B2 (en) Method for determining properties of a formation
CN105758780A (zh) 一种低渗透气藏非均质复合压力衰竭程度测试方法
Jianxin et al. Productivity Evaluation of Pronounced Heterogeneous Gas Reservoir Drilled at High Overbalance
Zweers et al. Relative Permeability Measurements on Heterogeneous Samples. A Pragmatic Approach
Hawker Direct gas in mud measurement at the well site
Blinov et al. Exploration and Formation Testing of Low Permeability Clastic Reservoir in Hot Slim-Diameter Open-Hole Well
WO2023064325A1 (en) Combination of a surface well testing facility and a cable formation tester with an active circulation system for obtaining inflow and measuring formation fluid parameters on the surface
CN114461970A (zh) 一种利用mdt电缆地层测试资料计算储层含油饱和度的方法

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
EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20201023

Assignee: Guangxi ronghua Ship Technology Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2023980053987

Denomination of invention: A Method for Calculating Gas Well Productivity in High Pressure Carbonate Water Gas Reservoirs

Granted publication date: 20220628

License type: Common License

Record date: 20231227

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20201023

Assignee: GUANGZHOU JUMING AUTOMATIC CONTROL TECHNOLOGY CO.,LTD.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2023980054328

Denomination of invention: A Method for Calculating Gas Well Productivity in High Pressure Carbonate Water Gas Reservoirs

Granted publication date: 20220628

License type: Common License

Record date: 20240103

Application publication date: 20201023

Assignee: GUANGZHOU TIANYI TECHNOLOGY Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2023980054325

Denomination of invention: A Method for Calculating Gas Well Productivity in High Pressure Carbonate Water Gas Reservoirs

Granted publication date: 20220628

License type: Common License

Record date: 20240103

Application publication date: 20201023

Assignee: Guangzhou Xiangyun Technology Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2023980054323

Denomination of invention: A Method for Calculating Gas Well Productivity in High Pressure Carbonate Water Gas Reservoirs

Granted publication date: 20220628

License type: Common License

Record date: 20240103

Application publication date: 20201023

Assignee: Laishi (Guangzhou) Digital Technology Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2023980054322

Denomination of invention: A Method for Calculating Gas Well Productivity in High Pressure Carbonate Water Gas Reservoirs

Granted publication date: 20220628

License type: Common License

Record date: 20240103

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20201023

Assignee: Shandong Kaicheng Geoprofessions#Geological engineering Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2024980000294

Denomination of invention: A Method for Calculating Gas Well Productivity in High Pressure Carbonate Water Gas Reservoirs

Granted publication date: 20220628

License type: Common License

Record date: 20240110

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20201023

Assignee: Zhongke Liben (Chongqing) Technology Development Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2024980000644

Denomination of invention: A Method for Calculating Gas Well Productivity in High Pressure Carbonate Water Gas Reservoirs

Granted publication date: 20220628

License type: Common License

Record date: 20240119

Application publication date: 20201023

Assignee: Yuao Holdings Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2024980000642

Denomination of invention: A Method for Calculating Gas Well Productivity in High Pressure Carbonate Water Gas Reservoirs

Granted publication date: 20220628

License type: Common License

Record date: 20240119

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20201023

Assignee: FOSHAN TENGPU INDUSTRIAL DESIGN CO.,LTD.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2024980003021

Denomination of invention: A Method for Calculating Gas Well Productivity in High Pressure Carbonate Water Gas Reservoirs

Granted publication date: 20220628

License type: Common License

Record date: 20240322

Application publication date: 20201023

Assignee: FOSHAN YIQING TECHNOLOGY Co.,Ltd.

Assignor: Chongqing University of Science & Technology

Contract record no.: X2024980003019

Denomination of invention: A Method for Calculating Gas Well Productivity in High Pressure Carbonate Water Gas Reservoirs

Granted publication date: 20220628

License type: Common License

Record date: 20240322