CN111810140A - 一种单采气井的最终累计采气量预测方法 - Google Patents

一种单采气井的最终累计采气量预测方法 Download PDF

Info

Publication number
CN111810140A
CN111810140A CN202010750697.XA CN202010750697A CN111810140A CN 111810140 A CN111810140 A CN 111810140A CN 202010750697 A CN202010750697 A CN 202010750697A CN 111810140 A CN111810140 A CN 111810140A
Authority
CN
China
Prior art keywords
gas production
final accumulated
well
coefficient
regression
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
CN202010750697.XA
Other languages
English (en)
Other versions
CN111810140B (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.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202010750697.XA priority Critical patent/CN111810140B/zh
Publication of CN111810140A publication Critical patent/CN111810140A/zh
Application granted granted Critical
Publication of CN111810140B publication Critical patent/CN111810140B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing 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
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Geology (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • General Business, Economics & Management (AREA)
  • Marine Sciences & Fisheries (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Health & Medical Sciences (AREA)
  • Animal Husbandry (AREA)
  • Agronomy & Crop Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明了公开了一种单采气井的最终累计采气量预测方法,通过静态地质条件、压裂改造情况,预测单采气井最终累计采气量,减少了传统预测最终累计采气量方法中所需参数,降低了参数获取过程中可能造成的误差,通过反映储层条件的地质参数及反映储层改造效果的压裂施工参数,简单、快速多元线性回归出最终累计采气量预测公式,节约开发成本,促进了气井全生命周期高效管理及低产低效井挖潜方式优化,确保气田经济效益开发。

Description

一种单采气井的最终累计采气量预测方法
技术领域
本发明涉及油气田勘探开发领域,具体为一种单采气井的最终累计采气量预测方法。
背景技术
目前,单井最终累计采气量多是利用RTA等模拟软件,以历史生产数据为基础,通过设置原始井口及地层压力、静态储层参数及生产工艺参数等,选择Blasingame、Fetkovich等产量分析模型或Flowing Material Balance(流动物质平衡法)方法进行预测。这些方法需要在软件中设置反映油气藏初始状态的参数(压力、流体性质、孔隙度、渗透率、含油气饱和度等)和开采方式的参数(油管特征、套管特征等),部分参数获取难度大、耗费时间长,不同模型模拟结果有时差别较大。
储层条件是一口气井最终累积产气量的基础,压裂改造效果一定程度上决定了采出程度,基于此,通过反映储层条件的静态地质参数及压裂改造效果的压裂施工参数预测一口气井最终累计采气量,促进气井全生命周期高效管理及低产低效井挖潜方式优化,确保气田经济效益开发,目前此方面研究相对较少。
发明内容
本发明的提供了一种单采气井的最终累计采气量预测方法,目的是为了解决现有预测方法中,部分参数获取难度大、耗费时间长的问题。
本发明的目的是通过以下技术手段实现的,一种单采气井的最终累计采气量预测方法,包括以下步骤:
第一步,获取单试井预测最终累计采气量Gp;
第二步,获取已知单试井射孔段所在位置气层储层的声波时差AC,其中声波时差AC作为回归参数X1;
第三步,得到射孔段气层储层地层系数和储能系数,将地层系数和储能系数分别作为另外两个静态地质回归参数X2和X3;
第四步,获取射孔段压裂改造段的破裂压力、压裂液入地总量、平均砂比和反排率这四个压裂施工参数,分别作为参与回归参数X4、X5、X6和X7;
第五步,将静态地质、压裂施工七个回归参数X1、X2、X3、X4、X5、X6、X7与最终累计采气量Gp分别绘制交汇图;
第六步,以最终累计采气量Gp为因变量,以静态地质、压裂施工七个回归参数X1、X2、X3、X4、X5、X6、X7为自变量,进行多元线性回归,确定各回归参数的权重系数及回归常数c,X1、X2、X3、X4、X5、X6和X7的权重系数依次为a1、a2、a3、a4、a5、a6和a7
第七步,获得待测井的七个回归参数,根据第六步中的权重系数和回归常数c得到待测井的最终累计采气量Gp。
所述第三步中,地层系数和储能系数的获取方法为,获取已知单试井射孔段所在位置气层储层的测井解释气层厚度H、孔隙度Φ、渗透率K、含气饱和度Sg,则地层系数为K*H,储能系数为H*Φ*Sg。
所述第五步中,绘制交汇图后还剔除奇异值。
所述第七步中,待测井的最终累计采气量Gp计算公式为,Gp=a1*X1+a2*X2+a3*X3+a4*X4+a5*X5+a6*X6+a7*X7+c,
Gp:最终累计采气量,单位104m3
a1、a2、a3、a4、a5、a6、a7:权重系数;
X1:声波时差,单位μs/m
X2:地层系数H*K,其中气层厚度H单位m,渗透率K单位10-3μm2
X3:储能系数H*Φ*Sg,其中孔隙度Φ、含气饱和度Sg单位%
X4:破裂压力,单位MPa
X5:入地总量,单位m3
X6:平均砂比,单位%
X7:反排率,单位%
C: 回归常数。
本发明的有益效果在于:根据静态地质条件和压裂改造情况,对单采气井的最终累计采气量进行预测,相比现有技术,参数获取难度小,耗费时间少,且预测误差低,从而达到气井全生命周期科学高效管理及低产低效井挖潜方式优化,确保气田经济效益开发。
附图说明
图1为本发明单采气井预测最终累计采气量与常规方法预测最终累采气量对比图;
以下将结合附图及实施例对本发明做进一步详细说明。
具体实施方式
【实施例1】
一种单采气井的最终累计采气量预测方法,包括以下步骤:
第一步,获取单试井预测最终累计采气量Gp;
第二步,获取已知单试井射孔段所在位置气层储层的声波时差AC,其中声波时差AC作为回归参数X1;
第三步,得到射孔段气层储层地层系数和储能系数,将地层系数和储能系数分别作为另外两个静态地质回归参数X2和X3;
第四步,获取射孔段压裂改造段的破裂压力、压裂液入地总量、平均砂比和反排率这四个压裂施工参数,分别作为参与回归参数X4、X5、X6和X7;
第五步,将静态地质、压裂施工七个回归参数X1、X2、X3、X4、X5、X6、X7与最终累计采气量Gp分别绘制交汇图;
第六步,以最终累计采气量Gp为因变量,以静态地质、压裂施工七个回归参数X1、X2、X3、X4、X5、X6、X7为自变量,进行多元线性回归,确定各回归参数的权重系数及回归常数c,X1、X2、X3、X4、X5、X6和X7的权重系数依次为a1、a2、a3、a4、a5、a6和a7
第七步,获得待测井的七个回归参数,根据第六步中的权重系数和回归常数c得到待测井的最终累计采气量Gp。
所述第三步中,地层系数和储能系数的获取方法为,获取已知单试井射孔段所在位置气层储层的测井解释气层厚度H、孔隙度Φ、渗透率K、含气饱和度Sg,则地层系数为K*H,储能系数为H*Φ*Sg。
所述第五步中,绘制交汇图后还剔除奇异值。
所述第七步中,待测井的最终累计采气量Gp计算公式为,Gp=a1*X1+a2*X2+a3*X3+a4*X4+a5*X5+a6*X6+a7*X7+c,
Gp:最终累计采气量,单位104m3
a1、a2、a3、a4、a5、a6、a7:权重系数;
X1:声波时差,单位μs/m
X2:地层系数H*K,其中气层厚度H单位m,渗透率K单位10-3μm2
X3:储能系数H*Φ*Sg,其中孔隙度Φ、含气饱和度Sg单位%
X4:破裂压力,单位MPa
X5:入地总量,单位m3
X6:平均砂比,单位%
X7:反排率,单位%
C: 回归常数。
综上所述,具体工作流程为,获取已经开采的单试井的最终累计采气量Gp,测井解释气层厚度H、孔隙度Φ、渗透率K、含气饱和度Sg、声波时差AC,气层储层地层系数K*H、储能系数H*Φ*Sg、破裂压力、压裂液入地总量、平均砂比和反排率。以已经开采的单试井的最终累计采气量Gp为因变量,以已开采井的X1:声波时差;X2:地层系数H*K;X3:储能系数H*Φ*Sg;X4:破裂压力;X5:入地总量;X6:平均砂比;X7:反排率为自变量,进行多元线性回归,得到各个回归参数对应的权重系数,即a1、a2、a3、a4、a5、a6和a7,以及回归常数c,由此建立最终累计采气量预测公式Gp=a1*X1+a2*X2+a3*X3+a4*X4+a5*X5+a6*X6+a7*X7+c。
获取待测井的X1:声波时差;X2:地层系数H*K;X3:储能系数H*Φ*Sg;X4:破裂压力;X5:入地总量;X6:平均砂比;X7:反排率,带入到最终累计采气量预测公式中,根据已开采井得到的权重系数,即a1、a2、a3、a4、a5、a6和a7,以及回归常数c,即可得到待测井的最终累计采气量Gp。
减少了传统预测最终累计采气量方法中所需参数,降低了参数获取过程中可能造成的误差,通过反映储层条件的地质参数及反映储层改造效果的压裂施工参数,简单、快速多元线性回归出最终累计采气量预测公式,节约开发成本,促进了气井全生命周期高效管理及低产低效井挖潜方式优化,确保气田经济效益开发。
如图1所示,本发明在S气田应用50口井,预测无阻流量与现场测试无阻流量平均误差小于20%,满足现场生产需求。
【实施例2】
在实施例1的基础上,S气田SX1井,验证预测最终累计采气量计算结果。
1)分别获取20口单采气井常规方法预测最终累计采气量Gp。
2)获取20口井射孔段所在位置气层储层对应的测井解释气层厚度H、孔隙度Φ、渗透率K、含气饱和度Sg、声波时差AC,其中AC直接作为其中一个回归参数X1。
3)计算20口井射孔段气层储层地层系数K*H、储能系数H*Φ*Sg,将地层系数、储能系数分别作为另外两个静态地质回归参数X2、X3;
4)获取20口井射孔段压裂改造段的破裂压力、压裂液入地总量、平均砂比、反排率4个压裂施工参数,分别作为参与回归参数X4、X5、X6、X7;
5)将静态地质、压裂施工七个回归参数X1、X2、X3、X4、X5、X6、X7与试气无阻流量Gp分别绘制交汇图,剔除奇异值;
6)以20口井预测最终累采气量Gp为因变量,以静态地质、压裂施工七个回归参数X1、X2、X3、X4、X5、X6、X7为自变量,进行多元线性回归,得出各参数的权重系数分别为a1=76.7、a2=239.7、a3=3822.1、a4=-5.3、a5=12.4、a6=235.6、a7=46.9,公式常数c=-26583;
7)获取SX1井,即未开采井射孔段气层储层测井解释气层厚度H=5.8m、孔隙度Φ=7.18%、渗透率K=0.461×10-3μm2、含气饱和度Sg=70.1%,计算地层系数X2=2.67、储能系数X3=0.29;获取声波时差X1=218.5μs/m、破裂压力X4=40.9MPa、压裂液入地总量X5=363.1m3、平均砂比X6=16.8%、反排率X7=72.9%;
8)利用多元线性回归公式Gp=α1*X1+α2*X2+α3*X3+α4*X4+α5*X5+α6*X6+α7*X7+c,计算出SX1井最终累产气量3444.5×104m3,SX1井常规方法预测最终累产气量3935.4×104m3/d,误差12.5%,小于20%,满足生产实际需求。
【实施例3】
在实施例2的基础上,权重系数沿用实例2中的值,即a1=76.7、a2=239.7、a3=3822.1、a4=-5.3、a5=12.4、a6=235.6、a7=46.9,公式常数c=-26583;
2)获取SX2井射孔段气层储层测井解释气层厚度H=8.5m、孔隙度Φ=7.4%、渗透率K=0.653×10-3μm2、含气饱和度Sg=52.6%,计算地层系数X2=5.55、储能系数X3=0.33;获取声波时差X1=237.7μs/m、破裂压力X4=20.4MPa、压裂液入地总量X5=483.6m3、平均砂比X6=16.0%、反排率X7=50.0%;
3)利用多元线性回归公式Gp=α1*X1+α2*X2+α3*X3+α4*X4+α5*X5+α6*X6+α7*X7+c,计算出SX2井最终累产气量5878.1×104m3,SX2井常规方法预测最终累产气量5278.9×104m3/d,误差11.4%,小于20%,满足生产实际需求。
如图1所示,本发明所采用的预测方法,减少了传统预测最终累计采气量方法中所需参数,降低了参数获取过程中可能造成的误差,通过反映储层条件的地质参数及反映储层改造效果的压裂施工参数,简单、快速多元线性回归出最终累计采气量预测公式,节约开发成本,促进了气井全生命周期高效管理及低产低效井挖潜方式优化,确保气田经济效益开发。

Claims (4)

1.一种单采气井的最终累计采气量预测方法,其特征在于,包括以下步骤:
第一步,获取单试井预测最终累计采气量Gp;
第二步,获取已知单试井射孔段所在位置气层储层的声波时差AC,其中声波时差AC作为回归参数X1;
第三步,得到射孔段气层储层地层系数和储能系数,将地层系数和储能系数分别作为另外两个静态地质回归参数X2和X3;
第四步,获取射孔段压裂改造段的破裂压力、压裂液入地总量、平均砂比和反排率这四个压裂施工参数,分别作为参与回归参数X4、X5、X6和X7;
第五步,将静态地质、压裂施工七个回归参数X1、X2、X3、X4、X5、X6、X7与最终累计采气量Gp分别绘制交汇图;
第六步,以最终累计采气量Gp为因变量,以静态地质、压裂施工七个回归参数X1、X2、X3、X4、X5、X6、X7为自变量,进行多元线性回归,确定各回归参数的权重系数及回归常数c,X1、X2、X3、X4、X5、X6和X7的权重系数依次为a1、a2、a3、a4、a5、a6和a7
第七步,获得待测井的七个回归参数,根据第六步中的权重系数和回归常数c得到待测井的最终累计采气量Gp。
2.根据权利要求1所述的一种单采气井的最终累计采气量预测方法,其特征在于:所述第三步中,地层系数和储能系数的获取方法为,获取已知单试井射孔段所在位置气层储层的测井解释气层厚度H、孔隙度Φ、渗透率K、含气饱和度Sg,则地层系数为K*H,储能系数为H*Φ*Sg。
3.根据权利要求1所述的一种单采气井的最终累计采气量预测方法,其特征在于:所述第五步中,绘制交汇图后还剔除奇异值。
4.根据权利要求1所述的一种单采气井的最终累计采气量预测方法,其特征在于:所述第七步中,待测井的最终累计采气量Gp计算公式为,Gp=a1*X1+a2*X2+a3*X3+a4*X4+a5*X5+a6*X6+a7*X7+c,
Gp:最终累计采气量,单位104m3
a1、a2、a3、a4、a5、a6、a7:权重系数;
X1:声波时差,单位μs/m
X2:地层系数H*K,其中气层厚度H单位m,渗透率K单位10-3μm2
X3:储能系数H*Φ*Sg,其中孔隙度Φ、含气饱和度Sg单位%
X4:破裂压力,单位MPa
X5:入地总量,单位m3
X6:平均砂比,单位%
X7:反排率,单位%
C: 回归常数。
CN202010750697.XA 2020-07-30 2020-07-30 一种单采气井的最终累计采气量预测方法 Active CN111810140B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010750697.XA CN111810140B (zh) 2020-07-30 2020-07-30 一种单采气井的最终累计采气量预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010750697.XA CN111810140B (zh) 2020-07-30 2020-07-30 一种单采气井的最终累计采气量预测方法

Publications (2)

Publication Number Publication Date
CN111810140A true CN111810140A (zh) 2020-10-23
CN111810140B CN111810140B (zh) 2023-07-25

Family

ID=72863366

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010750697.XA Active CN111810140B (zh) 2020-07-30 2020-07-30 一种单采气井的最终累计采气量预测方法

Country Status (1)

Country Link
CN (1) CN111810140B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112360449A (zh) * 2020-12-11 2021-02-12 西南石油大学 一种应用于低渗透气藏和气井的可采气量预测方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102536200A (zh) * 2012-02-17 2012-07-04 中国石油化工股份有限公司 致密碳酸盐岩气层初始产能预测方法
CN104832166A (zh) * 2015-03-20 2015-08-12 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 一种页岩气水平井初期产能预测方法
CN107590550A (zh) * 2017-07-26 2018-01-16 长江大学 有关超低渗透储层油田产能评价及预测的方法
CN109958430A (zh) * 2019-03-10 2019-07-02 东北石油大学 复杂致密气藏产能预测方法
US20200211126A1 (en) * 2018-12-29 2020-07-02 Petrochina Company Limited Prediction method for shale oil and gas sweet spot region, computer device and computer readable storage medium

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102536200A (zh) * 2012-02-17 2012-07-04 中国石油化工股份有限公司 致密碳酸盐岩气层初始产能预测方法
CN104832166A (zh) * 2015-03-20 2015-08-12 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 一种页岩气水平井初期产能预测方法
CN107590550A (zh) * 2017-07-26 2018-01-16 长江大学 有关超低渗透储层油田产能评价及预测的方法
US20200211126A1 (en) * 2018-12-29 2020-07-02 Petrochina Company Limited Prediction method for shale oil and gas sweet spot region, computer device and computer readable storage medium
CN109958430A (zh) * 2019-03-10 2019-07-02 东北石油大学 复杂致密气藏产能预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李跃刚等: "基于气水相对渗透率曲线的产水气井开采效果评价――以苏里格气田致密砂岩气藏为例", 《天然气工业》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112360449A (zh) * 2020-12-11 2021-02-12 西南石油大学 一种应用于低渗透气藏和气井的可采气量预测方法

Also Published As

Publication number Publication date
CN111810140B (zh) 2023-07-25

Similar Documents

Publication Publication Date Title
CN109594968B (zh) 一种页岩气多段压裂水平井压后裂缝参数评价方法及系统
RU2663011C2 (ru) Система и способ выполнения операции интенсификации скважины
US9715026B2 (en) System and method for performing microseismic fracture operations
CN108894777B (zh) 一种分层压裂多层合采油气藏储层物性及裂缝特性参数的确定方法
CN108009705A (zh) 一种基于支持向量机技术的页岩储层可压性评价方法
CN112502701B (zh) 一种低渗透储层综合地质-工程的分类评价方法
CN102418511B (zh) 低渗透储层井下关井水井压降试井分析方法
US20060015310A1 (en) Method for simulation modeling of well fracturing
CN112036098A (zh) 一种深层油气藏水力裂缝扩展数值模拟的方法
WO2015095557A1 (en) Method for determining hydraulic fracture orientation and dimension
CN103857876A (zh) 用于执行井眼压裂作业的系统和方法
CN111456709B (zh) 一种基于测井曲线的水平井多级压裂分段分簇方法
CN108952657B (zh) 一种水平井平台压裂裂缝长度确定方法及装置
CN112016212A (zh) 一种基于渗流控制方程的储层纵向非均质性评价方法
Usui et al. Effect of poroelasticity on hydraulic fracture interactions
CN115510778A (zh) 一种陆相页岩储层无限级压裂工艺优化方法及系统
CN109426689B (zh) 水平井压裂裂缝的评价方法及系统
CN111810140A (zh) 一种单采气井的最终累计采气量预测方法
CN115809536A (zh) 一种页岩气井多段压裂改造的评价方法
CN112100707A (zh) 一种穿层压裂图版构建方法
CN115270533A (zh) 重复压裂设计方法、装置、存储介质及电子设备
CN112068198B (zh) 基于地震波全波形特征的裂缝破裂尺度的描述方法
RU2709260C1 (ru) Способ повышения эффективности разработки слабопроницаемых нефтяных залежей
CN109339760B (zh) 一种水平井一段多簇压裂裂缝条数诊断方法
CN107704646B (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