CN110219624A - 岩石孔隙收缩和束缚水膨胀条件下确定水驱气藏参数方法 - Google Patents
岩石孔隙收缩和束缚水膨胀条件下确定水驱气藏参数方法 Download PDFInfo
- Publication number
- CN110219624A CN110219624A CN201910317922.8A CN201910317922A CN110219624A CN 110219624 A CN110219624 A CN 110219624A CN 201910317922 A CN201910317922 A CN 201910317922A CN 110219624 A CN110219624 A CN 110219624A
- Authority
- CN
- China
- Prior art keywords
- gas reservoir
- water
- measured
- invasion
- enchroachment
- 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 153
- 238000000034 method Methods 0.000 title claims abstract description 17
- 230000009545 invasion Effects 0.000 claims abstract description 38
- 239000000463 material Substances 0.000 claims abstract description 13
- 238000004519 manufacturing process Methods 0.000 claims description 20
- 230000004941 influx Effects 0.000 claims description 10
- 230000001186 cumulative effect Effects 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000008602 contraction Effects 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 abstract 1
- 239000007789 gas Substances 0.000 description 101
- 238000010586 diagram Methods 0.000 description 5
- 102000003712 Complement factor B Human genes 0.000 description 3
- 108090000056 Complement factor B Proteins 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 3
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 239000008239 natural water Substances 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 241000282994 Cervidae Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000003860 storage 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
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
- E21B43/20—Displacing by water
-
- 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
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/25—Methods for stimulating production
- E21B43/26—Methods for stimulating production by forming crevices or fractures
-
- 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)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了岩石孔隙收缩和束缚水膨胀条件下确定水驱气藏参数方法,包括:S1、将考虑岩石孔隙收缩和束缚水膨胀条件的水驱气藏物质平衡方程转化为变形后的物质平衡方程;S2、将待测气藏水侵量计算公式代入变形后的物质平衡方程,得到水侵物质平衡指示曲线方程;S3、设待测气藏的水体倍数为10,在此水体倍数下绘制水侵物质平衡指示曲线;S4、设待测气藏的水体倍数分别为20、30、50、150,重复S3步骤,得到不同水体倍数下的水侵物质平衡指示曲线;S5、相关系数最大的指示曲线所对应的水体倍数为待测水驱气藏水体倍数,曲线的斜率为待测气藏水侵常数,曲线的截距为待测气藏动态储量。本发明原理可靠,操作简便,为气藏开发提供数据支撑和理论依据。
Description
技术领域
本发明属于气藏开发领域,具体涉及一种在岩石孔隙收缩和束缚水膨胀条件下确定水驱气藏水体倍数、水侵常数和动态储量的方法。
背景技术
我国天然气资源丰富,其中水驱气藏是我国气藏的重要组成部分。水驱气藏边底水构成水驱气藏水体,研究并准确确定水驱气藏水体大小对气藏的生产有重要意义,对储层气井产量的计算、动态分析、数值模拟非常重要。对水驱气藏水体大小的研究还有利于气藏中后期及时调整生产方案,提高水驱气藏的采收率。
目前计算水驱气藏水体大小的方法中,最优拟合法(廖恒杰,鹿克峰等.采用最优化拟合法计算水驱气藏动储量及水体大小[J].重庆科技学院学报.2018,20(2):59-62)使用Fetkovich拟稳态方法计算水驱气藏的水侵量,根据水侵量大小进而分析水体大小。异常高压凝析气田算法虽然公式推导严谨(吴克柳,李相方等.异常高压凝析气藏水侵量及水体大小计算方法[J].中国矿业大学学报.2013,42(1):105-111),可以用来计算水体大小,但其适用范围小,仅仅适用于异常高压凝析气藏,而且适用条件苛刻,必须是当压力波及整个天然水域,参与流动的天然水域范围固定时计算误差才会缩小,不适合应用于实际矿场。
发明内容
本发明的目的在于提供岩石孔隙收缩和束缚水膨胀条件下确定水驱气藏参数方法,该方法原理可靠,操作简便,能够准确合理地确定水驱气藏水体倍数、水侵常数和动态储量,为气藏开发提供数据支撑和理论依据。
为达到以上技术目的,本发明采用以下技术方案。
岩石孔隙收缩和束缚水膨胀条件下确定水驱气藏参数方法,依次包括以下步骤:
S1、将考虑岩石孔隙收缩和束缚水膨胀条件的水驱气藏物质平衡方程(李传亮.油藏工程原理第三版.石油工业出版社,2017年):
GBgi=(G-Gp)Bg+We-WpBw+CcGBgiΔP
转化为变形后的物质平衡方程:
GpBg+WpBw=G(CcBgiΔP+Bg-Bgi)+We
其中Gp为待测气藏的累积产气量,m3;
Bgi为待测气藏的气体原始体积系数;
Bg为待测气藏的气体体积系数;
We为待测气藏的水侵量,m3;
Wp为待测气藏的累积产水量,m3;
Bw为待测气藏产出水的体积系数;
G为待测气藏的动态储量,m3;
Cc为待测气藏岩石孔隙收缩和束缚水膨胀综合系数,Mpa-1;
ΔP为待测气藏初始压力与当前压力的压力差,Mpa;
S2、将待测气藏水侵量We的计算公式
代入变形后的物质平衡方程,得到水侵物质平衡指示曲线方程:
特殊地:当i=1时,ΔP1=(P0-P1)/2
其中We为待测气藏水侵量,m3;
B为待测气藏水侵常数,m3/MPa;
P0为待测气藏原始压力,MPa;
Pi为待测气藏i时刻的压力,MPa;
tn为气藏初始时刻到n时刻所用的生产时间,天;
ti-1为气藏初始时刻到i-1时刻所用的生产时间,天;
(tn-ti-1)D为与tn-ti-1时间段有关的无因次时间;
q(tn-ti-1)D为与tn-ti-1时间段有关的无因次时间所对应的无因次流量;
S3、设待测气藏的水体倍数为10,在此水体倍数下计算各个地层压力测点所对应的(GpBg+WpBw)/(CcBgiΔP+Bg-Bgi)值和Σ(Δpi)q(tn-ti-1)D/(CcBgiΔP+Bg-Bgi)值,绘制水侵物质平衡指示曲线;
S4、设待测气藏的水体倍数分别为20、30、50、150,重复S3步骤,得到不同水体倍数下的水侵物质平衡指示曲线;
S5、确定每条水驱气藏水侵物质平衡指示曲线的相关系数,相关系数最大的指示曲线所对应的水体倍数即为待测水驱气藏的水体倍数,曲线的斜率为待测气藏水侵常数,曲线的截距为待测气藏动态储量。
所述步骤S1中,通过下式计算待测气藏岩石孔隙收缩和束缚水膨胀综合系数Cc:
其中Cp为待测气藏的岩石孔隙收缩系数,MPa-1;
Cw为待测气藏的束缚水膨胀系数,MPa-1;
Swc为待测气藏的束缚水饱和度。
所述步骤S2中,通过van Everdingen-Hurst水侵计算公式(Klins,M.A.,Bouchard,A.J.,Cable,C.L.A Polynomial Approach to the van Everdingen-HurstDimensionless Variables for Water Encroachment,Society of PetroleumEngineers,1988.2)计算待测气藏水侵量We:
特殊地:当i=1时,ΔP1=(P0-P1)/2
其中B为待测气藏水侵常数,m3/Mpa;
P0为待测气藏原始压力,Mpa;
Pi为待测气藏i时刻的压力,Mpa;
tn为气藏初始时刻到n时刻所用的生产时间,天;
ti-1为气藏初始时刻到i-1时刻所用的生产时间,天;
(tn-ti-1)D为与tn-ti-1时间段有关的无因次时间;
q(tn-ti-1)D为与tn-ti-1时间段有关的无因次时间所对应的无因次流量。
所述步骤S3中,待测气藏的累积产气量Gp,待测气藏的累积产水量Wp,待测气藏的气体原始体积系数Bgi和待测气藏的气体体积系数Bg,待测气藏产出水的体积系数Bw,待测气藏原始压力P0和不同时刻地层压力Pi均从油田生产资料中获得。
通过待测气藏不同时刻地层压力,进而求得待测气藏初始压力与当前压力的压力差ΔP;通过待测气藏从初始时刻到n时刻和i-1时刻的生产时间tn和ti-1,可计算与tn-ti-1时间段有关的无因次时间(tn-ti-1)D与tn-ti-1时间段有关的无因次时间所对应的无因次流量q(tn-ti-1)D。
与现有技术相比,本发明具有以下有益效果:
本发明利用待测气藏数据和水驱气藏物质平衡水侵指示曲线可以确定出实际气藏水体的大小、水侵常数和动态储量。该方法获取数据快速简单,利用多组数据得到不同水侵指示曲线,根据得到的相关系数,得到最优结果。本发明提供的方法,提高了水体大小、水侵常数和动态储量确定的真实性和准确性,为气藏开发领域提供基础数据支撑和理论依据。
附图说明
图1为实施例中10倍水体倍数下拟合得到的水侵指示曲线示意图。
图2为实施例中20倍水体倍数下拟合得到的水侵指示曲线示意图。
图3为实施例中30倍水体倍数下拟合得到的水侵指示曲线示意图。
图4为实施例中50倍水体倍数下拟合得到的水侵指示曲线示意图。
图5为实施例中150倍水体倍数下拟合得到的水侵指示曲线示意图。
具体实施方式
下面根据附图和实例进一步说明本发明,以便于本技术领域的技术人员理解本发明。但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,均在保护之列。
实施例1
岩石孔隙收缩和束缚水膨胀条件下确定水驱气藏参数方法,依次包括以下步骤:
S1、把考虑岩石孔隙收缩和束缚水膨胀的水驱气藏物质平衡方程转化为变形后的物质平衡方程;
S2、将Van Everdingen-Hurst水侵计算公式代入变形后的物质平衡方程,得到水侵物质平衡指示曲线方程;
S3、从油田生产资料获得待测气藏的累积产气量Gp,待测气藏的累积产水量Wp,待测气藏的气体原始体积系数Bgi和待测气藏的气体体积系数Bg,待测气藏原始压力P0和不同时刻地层压力Pi;通过待测气藏的岩石孔隙收缩系数Cp、待测气藏的束缚水膨胀系数Cw和待测气藏束缚水饱和度Swc,可求取待测气藏岩石孔隙收缩和束缚水膨胀综合系数Cc;通过待测气藏不同时刻地层压力,进而求得待测气藏初始压力与当前压力的压力差ΔP;通过待测气藏从初始时刻到n时刻和i-1时刻的生产时间,可计算与tn-ti-1时间段有关的无因次时间(tn-ti-1)D、与tn-ti-1时间段有关的无因次时间所对应的无因次流量q(tn-ti-1)D。
设置水体倍数为10,在此水体倍数下计算各个地层压力测点所对应的(GpBg+WpBw)/(CcBgiΔP+Bg-Bgi)值和Σ(Δpi)q(tn-ti-1)D/(CcBgiΔP+Bg-Bgi)值,以Σ(Δpi)q(tn-ti-1)D/(CcBgiΔP+Bg-Bgi)为横坐标,(GpBg+WpBw)/(CcBgiΔP+Bg-Bgi)为纵坐标,建立直角坐标系,绘制水侵物质平衡指示曲线;
S4、分别设置水体倍数为20、30、50、150,重复S3步骤,得到多条水驱气藏水侵指示曲线;
S5、确定每条水驱气藏水侵指示曲线的相关系数,相关系数最大的水驱气藏水侵指示曲线所对应的水体倍数即为实际待测水驱气藏水体倍数,曲线的斜率为待测气藏水侵常数,曲线的截距为待测气藏动态储量。
本实施例随机选取了一组实验数据在图1到图5中展示。其中图1、图2、图3、图4、图5分别是在10倍、20倍、30倍、50倍、150倍水体大小下拟合得到的水侵指示曲线图。不同水体倍数下水侵指示曲线相关系数R2不同,对比不同相关系数,图1的水侵指示曲线的相关系数最高,即实际气藏的水体倍数为10,在10倍水体大小条件下即可求得待测气藏的水侵常数为1.3273×105m3/Mpa,待测气藏的动态储量为3.22×108m3。
Claims (2)
1.岩石孔隙收缩和束缚水膨胀条件下确定水驱气藏参数方法,依次包括以下步骤:
S1、将考虑岩石孔隙收缩和束缚水膨胀条件的水驱气藏物质平衡方程:
GBgi=(G-Gp)Bg+We-WpBw+CcGBgiΔP
转化为变形后的物质平衡方程:
GpBg+WpBw=G(CcBgiΔP+Bg-Bgi)+We
其中Gp为待测气藏的累积产气量,m3;
Bgi为待测气藏的气体原始体积系数;
Bg为待测气藏的气体体积系数;
We为待测气藏的水侵量,m3;
Wp为待测气藏的累积产水量,m3;
Bw为待测气藏产出水的体积系数;
G为待测气藏的动态储量,m3;
Cc为待测气藏岩石孔隙收缩和束缚水膨胀综合系数,Mpa-1;
△P为待测气藏初始压力与当前压力的压力差,Mpa;
S2、将待测气藏水侵量We的计算公式
代入变形后的物质平衡方程,得到水侵物质平衡指示曲线方程:
当i=1时,△P1=(P0-P1)/2
其中We为待测气藏水侵量,m3;
B为待测气藏水侵常数,m3/MPa;
P0为待测气藏原始压力,MPa;
Pi为待测气藏i时刻的压力,MPa;
tn为气藏初始时刻到n时刻所用的生产时间,天;
ti-1为气藏初始时刻到i-1时刻所用的生产时间,天;
(tn-ti-1)D为与tn-ti-1时间段有关的无因次时间;
q(tn-ti-1)D为与tn-ti-1时间段有关的无因次时间所对应的无因次流量;
S3、设待测气藏的水体倍数为10,在此水体倍数下计算各个地层压力测点所对应的(GpBg+WpBw)/(CcBgiΔP+Bg-Bgi)值和Σ(Δpi)q(tn-ti-1)D/(CcBgiΔP+Bg-Bgi)值,绘制水侵物质平衡指示曲线;
S4、设待测气藏的水体倍数分别为20、30、50、150,重复S3步骤,得到不同水体倍数下的水侵物质平衡指示曲线;
S5、确定每条水驱气藏水侵物质平衡指示曲线的相关系数,相关系数最大的指示曲线所对应的水体倍数即为待测水驱气藏的水体倍数,曲线的斜率为待测气藏水侵常数,曲线的截距为待测气藏动态储量。
2.如权利要求1所述的岩石孔隙收缩和束缚水膨胀条件下确定水驱气藏参数方法,其特征在于,所述步骤S1中,通过下式计算待测气藏岩石孔隙收缩和束缚水膨胀综合系数Cc:
其中Cp为待测气藏的岩石孔隙收缩系数,MPa-1;
Cw为待测气藏的束缚水膨胀系数,MPa-1;
Swc为待测气藏的束缚水饱和度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910317922.8A CN110219624A (zh) | 2019-04-19 | 2019-04-19 | 岩石孔隙收缩和束缚水膨胀条件下确定水驱气藏参数方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910317922.8A CN110219624A (zh) | 2019-04-19 | 2019-04-19 | 岩石孔隙收缩和束缚水膨胀条件下确定水驱气藏参数方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110219624A true CN110219624A (zh) | 2019-09-10 |
Family
ID=67819810
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910317922.8A Pending CN110219624A (zh) | 2019-04-19 | 2019-04-19 | 岩石孔隙收缩和束缚水膨胀条件下确定水驱气藏参数方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110219624A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110688764A (zh) * | 2019-10-09 | 2020-01-14 | 西南石油大学 | 一种水驱气藏水体倍数的计算方法 |
CN111507537A (zh) * | 2020-04-26 | 2020-08-07 | 中国石油大学(北京) | 一种水驱气藏的储量和水侵量的预测方法及装置 |
CN116006132A (zh) * | 2023-01-29 | 2023-04-25 | 西南石油大学 | 一种水驱气藏考虑气中凝析水影响的水侵量计算方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105443120A (zh) * | 2015-10-26 | 2016-03-30 | 中国石油化工股份有限公司 | 一种海相整装气田边水水侵早期特征分析方法 |
CN109339775A (zh) * | 2018-10-25 | 2019-02-15 | 西南石油大学 | 一种确定水驱气藏水体大小的方法 |
-
2019
- 2019-04-19 CN CN201910317922.8A patent/CN110219624A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105443120A (zh) * | 2015-10-26 | 2016-03-30 | 中国石油化工股份有限公司 | 一种海相整装气田边水水侵早期特征分析方法 |
CN109339775A (zh) * | 2018-10-25 | 2019-02-15 | 西南石油大学 | 一种确定水驱气藏水体大小的方法 |
Non-Patent Citations (1)
Title |
---|
童孝华 等: "《油气藏工程基础》", 31 March 1996, 石油工业出版社 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110688764A (zh) * | 2019-10-09 | 2020-01-14 | 西南石油大学 | 一种水驱气藏水体倍数的计算方法 |
CN111507537A (zh) * | 2020-04-26 | 2020-08-07 | 中国石油大学(北京) | 一种水驱气藏的储量和水侵量的预测方法及装置 |
CN111507537B (zh) * | 2020-04-26 | 2023-06-20 | 中国石油大学(北京) | 一种水驱气藏的储量和水侵量的预测方法及装置 |
CN116006132A (zh) * | 2023-01-29 | 2023-04-25 | 西南石油大学 | 一种水驱气藏考虑气中凝析水影响的水侵量计算方法 |
CN116006132B (zh) * | 2023-01-29 | 2024-05-14 | 西南石油大学 | 一种水驱气藏考虑气中凝析水影响的水侵量计算方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110110435B (zh) | 一种基于广义管流渗流耦合的流动模拟及瞬变井分析方法 | |
CN110219624A (zh) | 岩石孔隙收缩和束缚水膨胀条件下确定水驱气藏参数方法 | |
CN103912269B (zh) | 页岩气储层地层破裂压力梯度录井确定方法 | |
CN104453804A (zh) | 一种气驱油藏开发动态监测与评价方法 | |
CN107435528B (zh) | 火山岩气藏气井配产的方法 | |
CN109339775A (zh) | 一种确定水驱气藏水体大小的方法 | |
CN111502652B (zh) | 一种三孔介质气藏水平井产量递减及生产动态预测方法 | |
CN112412411B (zh) | 一种压裂井多井系统试井分析方法、装置及存储介质 | |
CN106761613B (zh) | 一种二氧化碳驱替前缘的试井确定方法 | |
Zhou et al. | Investigation of microscopic pore structure and permeability prediction in sand-conglomerate reservoirs | |
CN104237107A (zh) | 地层中低渗透率储层的视渗透率解释方法及系统 | |
CN111241460B (zh) | 一种复杂致密储层孔隙度计算方法 | |
CN116205163B (zh) | 一种天然气水合物藏数值模拟中的地层稳定性分析方法 | |
CN109655394B (zh) | 一种孔喉结构参数约束下的核磁共振t2谱渗透率计算方法 | |
CN115994500A (zh) | 评估页岩气井压裂裂缝导流能力动态变化的方法及系统 | |
CN103590813B (zh) | 两层合采气井动态储量的测试方法、测试管柱及测试装置 | |
CN113622908B (zh) | 一种水侵气藏废弃地层压力的确定方法 | |
CN107288633A (zh) | 具有重力分异现象的凝析气藏地质储量计算方法 | |
CN113486538B (zh) | 一种非常规油气井产量预测及压裂效果评价方法 | |
CN111305806A (zh) | 自支撑裂缝导流能力的分析方法及装置 | |
CN110095584A (zh) | 一种储层油水饱和度校正方法 | |
CN112561223B (zh) | 油气井增产潜力的定量评价方法 | |
CN114493212A (zh) | 凝析气藏动态储量的评估方法 | |
CN113356843A (zh) | 针对地层的井壁稳定性分析方法、装置、介质及设备 | |
Bhatti et al. | Generation of Saturation Functions from Historic SCAL Data for Simulation of a Super-Giant Carbonate Reservoir |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190910 |