CN109033541B - 一种基于eur的压后页岩气藏非均质性评价方法 - Google Patents

一种基于eur的压后页岩气藏非均质性评价方法 Download PDF

Info

Publication number
CN109033541B
CN109033541B CN201810709858.3A CN201810709858A CN109033541B CN 109033541 B CN109033541 B CN 109033541B CN 201810709858 A CN201810709858 A CN 201810709858A CN 109033541 B CN109033541 B CN 109033541B
Authority
CN
China
Prior art keywords
eur
well
gas reservoir
shale gas
heterogeneous
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.)
Active
Application number
CN201810709858.3A
Other languages
English (en)
Other versions
CN109033541A (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.)
China University of Geosciences Beijing
Original Assignee
China University of Geosciences Beijing
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 China University of Geosciences Beijing filed Critical China University of Geosciences Beijing
Priority to CN201810709858.3A priority Critical patent/CN109033541B/zh
Publication of CN109033541A publication Critical patent/CN109033541A/zh
Application granted granted Critical
Publication of CN109033541B publication Critical patent/CN109033541B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

本发明提出了一种基于EUR(最终可采储量)的压后页岩气藏非均质性评价方法。该方法包括：获取单井EUR；计算每口气井的EUR累积贡献百分数和井数累积贡献百分数，绘制EUR洛伦兹曲线；通过所述单井EUR洛伦兹曲线计算非均质系数，建立页岩气藏非均质评价标准，定量评价页岩气藏非均质程度。该评价方法利用了EUR动态参数结合洛伦兹曲线，能够克服由于压后页岩气藏大量支撑剂的泵入改变储层物性空间分布带来的影响，能够很好的解决静态参数在评价压后页岩气藏非均质性不准确的问题，提高压后页岩气藏非均质性评价结果的精度，适合矿场应用。

Description

一种基于EUR的压后页岩气藏非均质性评价方法
技术领域
本发明涉及页岩气开采技术领域，更为具体地，涉及一种基于EUR的压后页岩气藏非均质性评价方法。
背景技术
气藏非均质性是指储层岩性、物性或含油气性等在三维空间分布及其变化的差异性，分为宏观非均质性和微观非均质性，而宏观非均质性又包括层间、平面及层内非均质性。气藏非均质性是影响地下流体运移及油气采收率的主要因素之一，对于强非均质储层，钻遇低效井的概率往往比较高，增加了勘探与开发决策的难度。由此，准确认识气藏的非均质特征，有利于勘探开发部署，降低开发风险，提高开发效益。
常规油气藏非均质性评价主要以储层的非均质性评价为主，一般采用渗透率的变异系数、突进系数和渗透率级差等静态参数描述储层非均质特征，上述参数值越大，则气藏非均质性越强。页岩气藏由于采用水平井多段压裂，大量支撑剂泵入地层形成人工裂缝网络，改变了页岩气藏的原始物性分布，使得页岩气藏具备人工气藏的特征，静态指标的评价只能适用于原始页岩气藏，无法反映出页岩压后形成的人工气藏特征。因此，基于常规储层的静态参数(渗透率)非均质性评价方法不适用于页岩气藏。
近些年，国内外学者在页岩气藏非均质性评价方法进行了积极的探索。一些新的静态参数，如应力、裂缝、有机碳含量和岩石脆性等静态参数被用来研究页岩气藏非均质性。另外，一些学者提出采用气藏动态参数评价气藏非均质性，如初期峰值产量、阶段产气量(高峰月产量或首年产气量)等。但是目前采用的气井初期产量及阶段累积产量受生产条件的影响较大。因而不管是静态参数或者目前所使用的动态参数在评级页岩气藏非均质性时都存在一定的不适性，急需形成一种评价压后页岩气藏非均质程度的方法。
发明内容
本发明的目的在于克服上述现有的页岩气藏非均质评价方法的不足，提供一种基于EUR(Estimated Ultimate Recovery)的压后页岩气藏非均质性评价方法。
本发明是通过下述技术方案实现：一种基于EUR的压后页岩气藏非均质性评价方法，包括以下步骤：
(1)通过递减分析法获取页岩气藏单井EUR；
建立递减分析模型中产量和累积产量表达式：
Figure BDA0001716397690000021
Figure BDA0001716397690000022
Figure BDA0001716397690000023
在双对数坐标中q/Gp～t呈线性关系，利用生产数据作图拟合确定参数a和m；
通过整体生产数据确定初期产量qi，设：
Figure BDA0001716397690000024
考虑井筒储集，产量关系可以表示为：
q＝qit(a，m)+q∞
在直角坐标系中做出q～t(a,m)的关系，确定初期产量qi和附加产量值q∞；
确定废弃产量和对应的时间，通过对应的时间预测最终可采储量：
Figure BDA0001716397690000025
Figure BDA0001716397690000026
其中：qi为初期产量，104m3/d；Gp为累积产量，104m3；q为产量，104m3/d；t为生产时间，d；D为递减率，无因次；a、m为与递减率有关的参数，无因次；q∞为考虑井筒储集的附加产量，104m3/d。
(2)计算每口气井EUR贡献值累积百分数和井数贡献值累积百分数，绘制单井EUR洛伦兹曲线；
①按单井EUR值递增序列进行排序，并对气井从1开始依次编号，按照井编号递增顺序，计算每口气井EUR贡献值累积百分数；
②按照井编号增序的顺序，计算井数贡献值累积百分数；
③分别以每口气井EUR贡献值累积百分数为横坐标、井数贡献值累积百分数为纵坐标绘制单井EUR洛伦兹曲线。
所述每口气井EUR贡献值累积百分数和所述井数贡献值累积百分数的计算按下式：
Figure BDA0001716397690000031
Figure BDA0001716397690000032
Figure BDA0001716397690000033
其中：i为气井编号；n为总井数；EURi为第i口井对应的EUR值，104m3；Li为第i口井EUR对整体贡献百分数，％；Lj为第j口井EUR对整体贡献百分数，％；Xi为1～i井EUR贡献值累积百分数，％；Yi为1～i井的井数贡献值累积百分数，％。
(3)通过所述单井EUR洛伦兹曲线计算非均质系数，建立页岩气藏非均质评价标准，定量评价页岩气藏非均质程度。
所述的非均质系数为原洛伦兹曲线中基尼系数模型，非均质系数是通过曲线拟合积分法来计算的。
所述非均质系数的计算公式为：
K＝S1/(S1+S2)
其中：K为非均质系数，K表示实际气藏EUR洛伦兹曲线偏离理想气藏EUR曲线的面积与理想气藏EUR曲线面积的比值；S1为EUR洛伦兹曲线与45°倾斜直线所围成的面积，表示实际气藏EUR洛伦兹曲线偏离理想气藏EUR曲线的程度；S2为EUR洛伦兹曲线与纵坐标轴和上坐标轴围成的面积。
K＝0，表示气藏完全均质，K＝1，表示气藏具极度非均质特点，一般页岩气藏非均质系数位于0～1之间。
所述曲线拟合积分法采用ExpAssoc方程，所述ExpAssoc方程表达式为：
Figure BDA0001716397690000034
通过对所述单井EUR洛伦兹曲线进行曲线拟合得到拟合参数值，代入积分公式求得面积，将所述面积值代入所述非均质系数的计算公式，进而计算得到所述非均质系数。
所述积分公式为：
Figure BDA0001716397690000041
其中：y为井数贡献值累积百分数，％；x为EUR贡献值累积百分数，％；y0、A1、t1、A2、t2为方程中的拟合参数；y0为偏移量，％；A1、A2为幅度，无因次；t1、t2为宽度，无因次。
所述建立页岩气藏非均质评价标准时，参考致密气藏非均质评价标准，利用所述非均质系数，定量评价页岩气藏非均质程度。
本发明的有益效果是：提供一种定量评价压后页岩气藏非均质性的方法，该方法解决了页岩气藏由于压裂施工，大量支撑剂的泵入改变了储层物性的空间分布，现有参数不能准确评价页岩气藏非均质性，导致页岩气产区低效井比例大，开发效益差的问题。
附图说明
图1为基于EUR的压后页岩气藏非均质性评价方法流程示意图；
图2为本发明实例提供的某页岩气产区4509口气井EUR分布曲线；
图3为通过本发明绘制出的单井EUR洛伦兹曲线。
具体实施方式
下面结合实施例，对本发明作进一步详细说明，但本发明的实施方式不仅限于此。
图1示出了根据本发明的基于EUR的压后页岩气藏非均质性评价方法的流程。
如图1所示，本发明提供的一种基于EUR的压后页岩气藏非均质性评价方法，包括以下步骤：
步骤S1：获取页岩气藏单井EUR；
步骤S2：计算每口气井EUR贡献值累积百分数和井数贡献值累积百分数，绘制单井EUR洛伦兹曲线；
步骤S3：通过所述单井EUR洛伦兹曲线计算非均质系数，建立页岩气藏非均质评价标准，定量评价页岩气藏非均质程度。
为了便于理解，本发明以一个具体实例进行说明。
以我国西部某盆地某页岩气藏为例，采用递减分析法获取页岩气藏产区4509口气单井EUR。
图2为采用递减分析法获取的该页岩气藏产区4509口气井单井EUR分布曲线。
计算该页岩气藏每口气井EUR贡献值累积百分数和井数贡献值累积百分数，绘制单井EUR洛伦兹曲线。
按照井编号递增顺序，计算每口气井EUR贡献值累积百分数：
Figure BDA0001716397690000051
Figure BDA0001716397690000052
式中：i为气井编号；n为总井数；EURi为第i口井对应的EUR值，104m3；Li为第i口井EUR对整体贡献百分数，％；Lj为第j口井EUR对整体贡献百分数，％；
Xi为1～i井EUR贡献值累积百分数，％。
按照井编号增序的顺序，计算井数贡献值累积百分数：
Figure BDA0001716397690000053
式中：Yi为1～i井的井数贡献值累积百分数，％。
分别以每口气井EUR贡献值累积百分数Xi为横坐标、井数贡献值累积百分数Yi为纵坐标绘制EUR洛伦兹曲线，图3为区块单井EUR洛伦兹曲线。
通过计算该页岩气藏单井EUR洛伦兹曲线计算非均质系数，建立页岩气藏非均质评价标准，定量评价该页岩气藏非均质程度。
所述的非均质系数为原洛伦兹曲线中基尼系数模型，非均质系数是通过曲线拟合积分法来计算的。
非均质系数的计算公式为：K＝S1/(S1+S2)。
所述曲线拟合积分法采用ExpAssoc方程，所述ExpAssoc方程表达式为：
Figure BDA0001716397690000054
通过对所述单井EUR洛伦兹曲线进行曲线拟合得到拟合参数y0＝2.49、A1＝111.92、t1＝71.85、A2＝13.18、t2＝2.62，代入积分公式求得面积S1＝0.1433和S2＝0.3567
所述积分公式为：
Figure BDA0001716397690000061
其中：y为井数贡献值累积百分数，％；x为EUR贡献值累积百分数，％；y0、A1、t1、A2、t2为方程中的拟合参数；y0为偏移量，％；A1、A2为幅度，无因次；t1、t2为宽度，无因次。
将所述面积值S1＝0.1433和S2＝0.3567代入所述非均质系数的计算公式，进而计算非均质系数K＝0.287。
所述建立页岩气藏非均质评价标准时，参考致密气藏非均质评价标准，建立的页岩气藏非均质评价标准如表1所示，利用所述非均质系数K＝0.287，可以得到区块整体具弱非均质特征。
表1页岩气藏EUR非均质性评价标准
Figure BDA0001716397690000062
以上所述，仅是本发明的较佳实施例并非对本发明对任何形式上的限制，凡是依据本发明的技术、方法实质上以上实施例所作的任何简单修改、等同变化，均落入本发明的保护范围之内。

Claims (6)

1.一种基于EUR的压后页岩气藏非均质性评价方法，其特征在于，包括以下步骤：
步骤S1：获取页岩气藏单井EUR；所述页岩气藏单井EUR是通过递减分析方法计算得到的；所述递减分析方法主要包括以下步骤：
①建立递减分析模型中产量和累积产量表达式：
Figure FDA0002382537740000011
Figure FDA0002382537740000012
Figure FDA0002382537740000013
②在双对数坐标中q/Gp～t呈线性关系，利用生产数据作图拟合确定参数a和m；
③通过整体生产数据确定初期产量qi，设：
Figure FDA0002382537740000014
考虑井筒储集，产量关系可以表示为：
q＝qit(a，m)+q∞
在直角坐标系中做出q～t(a,m)的关系，确定初期产量qi和附加产量值q∞；
④确定废弃产量和对应的时间，通过对应的时间预测最终可采储量：
Figure FDA0002382537740000015
Figure FDA0002382537740000016
其中：qi为初期产量，104m3/d；Gp为累积产量，104m3；q为产量，104m3/d；t为生产时间，d；D为递减率，无因次；a、m为与递减率有关的参数，无因次；q∞为考虑井筒储集的附加产量，104m3；
步骤S2：计算每口气井EUR贡献值累积百分数和井数贡献值累积百分数，绘制单井EUR洛伦兹曲线；所述每口气井EUR贡献值累积百分数的计算、所述井数贡献值累积百分数的计算和所述单井EUR洛伦兹曲线的绘制方法包括：
①按单井EUR值递增序列进行排序，并对气井从1开始依次编号，按照井编号递增顺序，计算每口气井EUR贡献值累积百分数；
②按照井编号增序的顺序，计算井数贡献值累积百分数；
③分别以每口气井EUR贡献值累积百分数为横坐标、井数贡献值累积百分数为纵坐标绘制单井EUR洛伦兹曲线；
所述每口气井EUR贡献值累积百分数和所述井数贡献值累积百分数的计算按下式：
Figure FDA0002382537740000021
Figure FDA0002382537740000022
Figure FDA0002382537740000023
其中：i为气井编号；n为总井数；EURi为第i口井对应的EUR值，104m3；Li为第i口井EUR对整体贡献百分数，％；Lj为第j口井EUR对整体贡献百分数，％；Xi为1～i井EUR贡献值累积百分数，％；Yi为1～i井的井数贡献值累积百分数，％；
步骤S3：通过所述单井EUR洛伦兹曲线计算非均质系数，建立页岩气藏非均质评价标准，定量评价页岩气藏非均质程度。
2.根据权利要求1所述的基于EUR的压后页岩气藏非均质性评价方法，其特征在于：所述步骤S3中所述的非均质系数的计算公式为：
K＝S1/(S1+S2)。
其中：K为非均质系数，K表示实际气藏EUR洛伦兹曲线偏离理想气藏EUR曲线的面积与理想气藏EUR曲线面积的比值；S1为EUR洛伦兹曲线与45°倾斜直线所围成的面积，表示实际气藏EUR洛伦兹曲线偏离理想气藏EUR曲线的程度；S2为EUR洛伦兹曲线与纵坐标轴和上坐标轴围成的面积。
3.根据权利要求2所述的基于EUR的压后页岩气藏非均质性评价方法，其特征在于：所述的非均质系数为原洛伦兹曲线中基尼系数模型，非均质系数是通过曲线拟合积分法来计算的。
4.根据权利要求3所述的基于EUR的压后页岩气藏非均质性评价方法，其特征在于：所述曲线拟合积分法采用ExpAssoc方程，通过对所述单井EUR洛伦兹曲线进行曲线拟合得到拟合参数值，代入积分公式求得面积，将所述面积值代入所述非均质系数的计算公式，进而计算非均质系数。
5.根据权利要求4中所述的基于EUR的压后页岩气藏非均质性评价方法，其特征在于：所述的ExpAssoc方程表达式为：
Figure FDA0002382537740000031
所述积分公式为：
Figure FDA0002382537740000032
其中：y为井数贡献值累积百分数，％；x为EUR贡献值累积百分数，％；y0、A1、t1、A2、t2为方程中的拟合参数；y0为偏移量，％；A1、A2为幅度，无因次；t1、t2为宽度，无因次。
6.根据权利要求1中所述的基于EUR的压后页岩气藏非均质性评价方法，其特征在于：所述建立页岩气藏非均质评价标准时，参考致密气藏非均质评价标准，利用所述非均质系数，定量评价页岩气藏非均质程度。
CN201810709858.3A 2018-07-02 2018-07-02 一种基于eur的压后页岩气藏非均质性评价方法 Active CN109033541B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810709858.3A CN109033541B (zh) 2018-07-02 2018-07-02 一种基于eur的压后页岩气藏非均质性评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810709858.3A CN109033541B (zh) 2018-07-02 2018-07-02 一种基于eur的压后页岩气藏非均质性评价方法

Publications (2)

Publication Number Publication Date
CN109033541A CN109033541A (zh) 2018-12-18
CN109033541B true CN109033541B (zh) 2020-07-28

Family

ID=65521322

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810709858.3A Active CN109033541B (zh) 2018-07-02 2018-07-02 一种基于eur的压后页岩气藏非均质性评价方法

Country Status (1)

Country Link
CN (1) CN109033541B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111258894B (zh) * 2020-01-13 2023-08-08 中信银行股份有限公司 软件风险的评价方法、装置、存储介质及电子设备
CN112184915B (zh) * 2020-10-10 2022-11-11 东南大学 一种建筑冷热需求时空不同步特性的定量描述方法
CN113236207B (zh) * 2021-07-13 2021-09-10 西南石油大学 一种强非均质性储层中产水气井的定产量递减预测方法
CN113837890B (zh) * 2021-08-09 2024-07-12 中国石油天然气股份有限公司 一种基于数据统计的页岩气水平井井型优选方法
CN113642246B (zh) * 2021-08-24 2024-04-30 中国石油化工股份有限公司 一种页岩气水平井布井参数优化方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103075142A (zh) * 2013-01-14 2013-05-01 中国石油大学(华东) 一种注水开发油田堵水油井选井方法
CN107622165A (zh) * 2017-09-25 2018-01-23 西南石油大学 一种页岩气水平井重复压裂产能计算方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2951720B1 (en) * 2013-01-31 2018-05-02 Betazi, LLC Production analysis and/or forecasting methods, apparatus, and systems
CN104750896B (zh) * 2013-12-31 2016-07-06 中国石油化工股份有限公司 一种缝洞型碳酸盐岩油藏数值模拟方法
CN104405374B (zh) * 2014-10-16 2017-12-05 中国石油天然气股份有限公司 一种致密气藏储层应力敏感性的测量方法
CN105525909A (zh) * 2015-11-19 2016-04-27 薛云飞 分析油藏储层非均质性质的方法
CN105507893A (zh) * 2015-12-07 2016-04-20 西南石油大学 预测页岩储层体积改造后生产效果的方法
CN105528522A (zh) * 2015-12-15 2016-04-27 中国石油天然气股份有限公司 一种基于空间网格计算连续性油气藏资源量的方法及装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103075142A (zh) * 2013-01-14 2013-05-01 中国石油大学(华东) 一种注水开发油田堵水油井选井方法
CN107622165A (zh) * 2017-09-25 2018-01-23 西南石油大学 一种页岩气水平井重复压裂产能计算方法

Also Published As

Publication number Publication date
CN109033541A (zh) 2018-12-18

Similar Documents

Publication Publication Date Title
CN109033541B (zh) 一种基于eur的压后页岩气藏非均质性评价方法
CN112561144B (zh) 一种致密油压裂水平井产能主控因素评判与产能预测方法
CN109594968B (zh) 一种页岩气多段压裂水平井压后裂缝参数评价方法及***
CN104899411B (zh) 一种储层产能预测模型建立方法和***
CN111581786A (zh) 用于分析缝洞串联模式双孔复合储层的试井解释模型
CN107387051B (zh) 低渗透非均质油藏多段压裂水平井重复压裂选井的方法
CN109933745A (zh) 一种基于模糊判断的水合物钻采风险智能判断方法
CN109632459B (zh) 一种页岩可压裂性评价方法
CN106150477A (zh) 一种确定单井控制储量的方法
CN111456709B (zh) 一种基于测井曲线的水平井多级压裂分段分簇方法
CN108960651A (zh) 一种致密油气藏多级压裂水平井完井效率的综合评价方法
CN113283108B (zh) 一种定量评价页岩油储层可压裂性的方法及***
CN110608023B (zh) 稠油分层注汽的适应性界限分析评价方法
CN104695950A (zh) 火山岩油藏产能预测方法
CN108894777A (zh) 一种分层压裂多层合采油气藏储层物性及裂缝特性参数的确定方法
CN108661631A (zh) 一种页岩气井产量预测方法
CN111927417A (zh) 一种页岩气分段压裂水平井组储量动用状况评价方法
CN112196513A (zh) 基于水平井轨迹评价的龙马溪组页岩气井产能预测方法
CN113361111B (zh) 一种低渗透油藏储层特征模型的划分方法
CN108664679A (zh) 一种油气井生产数据分析方法
CN112211626B (zh) 一种非均质气藏气井产能试井测试类型的优选方法
CN111950112A (zh) 一种适用于底部封闭的碳酸盐岩储层动态分析方法
CN115434698B (zh) 基于光电吸收截面指数的地层钙质胶结物含量计算方法
CN108664677A (zh) 一种油气井生产数据分析方法
CN111523696A (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