CN114970153B - 一种油气藏型地下储气库多周期注采动态库容计算方法 - Google Patents
一种油气藏型地下储气库多周期注采动态库容计算方法 Download PDFInfo
- Publication number
- CN114970153B CN114970153B CN202210580896.XA CN202210580896A CN114970153B CN 114970153 B CN114970153 B CN 114970153B CN 202210580896 A CN202210580896 A CN 202210580896A CN 114970153 B CN114970153 B CN 114970153B
- Authority
- CN
- China
- Prior art keywords
- porosity
- rock
- net
- stress
- overburden
- 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
Links
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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/70—Combining sequestration of CO2 and exploitation of hydrocarbons by injecting CO2 or carbonated water in oil wells
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)
- Filling Or Discharging Of Gas Storage Vessels (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及数据计算技术领域,具体涉及一种油气藏型地下储气库多周期注采动态库容计算方法,包括获取不同物性实验岩样的岩石常压孔隙度、不同实验净应力和不同实验净应力对应的覆压孔隙度;建立实验岩样的不可逆孔隙度损失率模型和建立岩石初始覆压孔隙度计算模型,然后基于上述两个模型建立初始净应力下的岩石覆压孔隙度计算模型;基于初始净应力下的岩石覆压孔隙度计算模型和岩石初始覆压孔隙度计算模型建立储气库注采动态库容计算模型;将储气单元平均参数输入储气库注采动态库容计算模型得到特定储层孔隙度和特定储气库运行周期的储气库库容,解决了因缺乏储气库动态库容评价方法而导致无法准确制定储气库运行方案的问题。
Description
技术领域
本发明涉及数据计算技术领域,尤其涉及一种油气藏型地下储气库多周期注采动态库容计算方法。
背景技术
油气藏型地下储气库周期性注采运行特征,致使储层岩石的净应力周期性变化,从而产生交变应力敏感,严重损害其岩石物性,导致储气库运行过程中储气库库容大幅减小,影响储气库的正常运行。因此,理清油气藏型储气库多周期注采过程中的库容动态变化规律是准确预测和把握储气库运行动态,保证储气库安全高效运行的重要前提。
目前,现有技术公开的计算储气库库容的方法有基于储层参数的静态法和基于动态监测资料的动态法,但目前缺乏储气库动态库容评价方法,会导致在储气库运行过程中无法准确把握库容动态变化情况,进而造成无法准确制定储气库运行方案的问题。
发明内容
本发明的目的在于提供一种油气藏型地下储气库多周期注采动态库容计算方法,旨在解决因缺乏储气库动态库容评价方法而导致无法准确制定储气库运行方案的问题。
为实现上述目的,本发明提供了一种油气藏型地下储气库多周期注采动态库容计算方法,包括以下步骤:
获取不同物性实验岩样的岩石常压孔隙度、不同实验净应力和所述不同实验净应力对应的覆压孔隙度;
基于所述不同物性实验岩样的岩石常压孔隙度、所述不同实验净应力和所述不同实验净应力对应的覆压孔隙度建立所述实验岩样的不可逆孔隙度损失率模型;
基于所述不同物性实验岩样的岩石常压孔隙度和所述不同实验岩样的初始覆压孔隙度建立岩石初始覆压孔隙度计算模型;
基于所述不可逆孔隙度损失率模型和所述岩石初始覆压孔隙度计算模型建立初始净应力下的岩石覆压孔隙度计算模型;
基于所述初始净应力下的岩石覆压孔隙度计算模型建立储气库注采动态库容计算模型;
将储气单元平均参数输入所述储气库注采动态库容计算模型得到特定储层孔隙度和特定储气库运行周期的储气库库容。
其中,所述获取不同物性实验岩样的岩石常压孔隙度、所述不同实验净应力和所述不同实验净应力对应的覆压孔隙度的具体方式为:
对实验岩样开展常压孔隙度测试实验,得到常压孔隙度;
对所述实验岩样开展多周期循环载荷覆压孔隙体积实验,得到覆压孔隙体积;
基于所述覆压孔隙体积计算所述实验岩样在不同净应力状态下对应的岩石常压覆压孔隙度。
其中,所述基于所述不同物性实验岩样的岩石常压孔隙度、所述不同实验净应力和所述不同实验净应力对应的覆压孔隙度建立所述实验岩样的不可逆孔隙度损失率模型的具体方式为:
基于所述岩石不同净应力及其对应的覆压孔隙度拟合覆压孔隙度比值与无因次净应力关系,得到第一拟合参数;
基于所述岩石覆压孔隙度拟合孔隙度应力敏感指数与常压孔隙度关系,得到第二拟合参数;
基于所述第一拟合参数和所述第二拟合参数建立所述实验岩样的不可逆孔隙度损失率模型。
其中,所述基于所述岩石覆压孔隙度拟合覆压孔隙度比值与无因次净应力关系,得到第一拟合参数的具体方式为:
将所述岩石覆压孔隙度除以初始覆压孔隙度,得到覆压孔隙度比值;
将净应力除以初始净应力得到无因次净应力;
基于所述覆压孔隙度比值和所述无因次净应力,得到第一拟合参数。
其中,所述基于所述岩石常压覆压孔隙度拟合孔隙度应力敏感指数与常压孔隙度关系,得到第二拟合参数的具体方式为:
采用幂函数回归分析建立所述净应力增大过程应力敏感指数与常压孔隙度的关系模型,得到净应力增大过程应力敏感指数关系模型;
将所述净应力减小过程的系数进行算术平均处理,建立净应力减小过程应力敏感指数的关系模型,得到净应力减小过程应力敏感指数关系模型;
基于所述净应力增大过程应力敏感指数关系模型和所述净应力减小过程应力敏感指数关系模型得到第二拟合参数。
本发明的一种油气藏型地下储气库多周期注采动态库容计算方法,通过包括获取不同物性实验岩样的岩石常压孔隙度、不同实验净应力和所述不同实验净应力对应的覆压孔隙度;基于所述不同物性实验岩样的岩石常压孔隙度、所述不同实验净应力和所述不同实验净应力对应的覆压孔隙度建立所述实验岩样的不可逆孔隙度损失率模型;基于所述不同物性实验岩样的岩石常压孔隙度、所述不同实验净应力和所述不同实验净应力对应的覆压孔隙度建立初始覆压孔隙度计算模型;基于所述不可逆孔隙度损失率模型和所述岩石初始覆压孔隙度计算模型建立初始净应力下的岩石覆压孔隙度计算模型;基于所述初始净应力下的岩石覆压孔隙度计算模型计算储气库运行周期初始覆压孔隙度;基于所述初始净应力下的岩石覆压孔隙度计算模型建立储气库注采动态库容计算模型;将储气单元平均参数输入所述储气库注采动态库容计算模型得到特定储层孔隙度和特定储气库运行周期的储气库库容,解决了因缺乏储气库动态库容评价方法而导致无法准确制方案的问题。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明提供的一种油气藏型地下储气库多周期注采动态库容计算方法的流程图。
图2是覆压孔隙度与净应力关系曲线图。
图3是覆压孔隙度比值与无因次净应力关系曲线图。
图4是孔隙度应力敏感指数与常压孔隙度关系图。
图5是系数c和d差值与净应力变化过程关系图。
图6是孔隙度应力损失率与孔隙度关系图。
图7是储层A库容随储气库运行周期变化曲线图。
图8是储层B库容随储气库运行周期变化曲线图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
请参阅图1至图8,本发明提供一种油气藏型地下储气库多周期注采动态库容计算方法,包括以下步骤:
S1获取不同物性实验岩样的岩石常压孔隙度、不同实验净应力和所述不同实验净应力对应的覆压孔隙度;
具体方式为:
S11对实验岩样开展常压孔隙度测试实验,得到常压孔隙度;
S12对所述实验岩样开展多周期循环载荷覆压孔隙体积实验,得到覆压孔隙体积;
具体的,开展不同物性实验岩样多周期循环载荷覆压孔隙体积实验:根据现有实验测试系统和实验测试方法,开展不同物性实验岩样多周期循环载荷覆压孔隙体积实验,获取不同物性实验岩样在不同净应力下的覆压孔隙体积。
S12基于所述覆压孔隙体积计算所述实验岩样在不同净应力状态下对应的岩石常压覆压孔隙度;
具体的,计算不同物性实验岩样各净应力状态下覆压孔隙度:
根据实验测定的不同物性实验岩样各净应力状态下覆压孔隙体积,通过下式计算不同物性实验岩样各净应力状态下覆压孔隙度:
式中:为i号实验岩样在第j个净应力变化周期内受第k个实验设计净应力σk状态下的覆压孔隙体积;/>为i号实验岩样外观体积;/>为i号实验岩样在第j个净应力变化周期内受第k个实验设计净应力σk状态下的覆压孔隙度;T为净应力变化周期。
S2基于所述不同物性实验岩样的岩石常压孔隙度、所述不同实验净应力和所述不同实验净应力对应的覆压孔隙度建立所述实验岩样的不可逆孔隙度损失率模型;
具体方式为:
S21基于所述岩石不同净应力及其对应的覆压孔隙度拟合覆压孔隙度比值与无因次净应力关系,得到第一拟合参数;
具体方式为:
S211将所述岩石覆压孔隙度除以初始覆压孔隙度,得到覆压孔隙度比值;
具体的,初始覆压孔隙度为第1个净应力变化过程初始净应力下的覆压孔隙度。
S212将净应力除以初始净应力得到无因次净应力;
S213基于所述覆压孔隙度比值和所述无因次净应力,得到第一拟合参数。
具体的,绘制覆压孔隙度比值与无因次净应力关系曲线,通过下式拟合覆压孔隙度比值与无因次净应力关系,获取拟合参数。
式中:为i号实验岩样初始覆压孔隙度(第1个净应力变化过程初始净应力下的覆压孔隙度);σi为第i个实验设计净应力;σmin为第1个净应力变化过程初始净应力;ai|jD为i号实验岩样拟合乘幂系数,1-ai|jD也称为不可逆孔隙度损失率;bi|jT为i号实验岩样拟合幂指数相反数,也称为孔隙度应力敏感指数;D为储气库运行周期。
S22基于所述岩石覆压孔隙度拟合孔隙度应力敏感指数与常压孔隙度关系,得到第二拟合参数;
具体方式为:
S221采用幂函数回归分析建立所述净应力增大过程应力敏感指数与常压孔隙度的关系模型,得到净应力增大过程应力敏感指数关系模型;
具体的,绘制应力敏感指数与孔隙度关系曲线,通过下式拟合孔隙度应力敏感指数与常压孔隙度关系,获取应力敏感指数与孔隙度拟合参数。
式中:ci|jT、di|jT均为i号实验岩样拟合方程回归系数;φi为i号实验岩样的常压孔隙度,%。
净应力增大时,随着净应力变化过程数增多,净应力对岩样造成了不同程度的塑性形变,因而应力敏感指数随净应力变化过程数的增多而逐渐降低,随着净应力变化过程数增多,岩样致密程度增大,岩样越来越难以被压缩,塑性形变越来越小,应力敏感指数随净应力变化过程数的增多而降低的趋势逐渐减缓;净应力下降时,仅弹性变形可以恢复,而各净应力上升过程的应力变化相同,各净应力上升过程的弹性变形规律基本相同,净应力下降时弹性变形恢复规律也基本相同,因而孔隙度应力敏感指数基本不随净应力变化过程数的增多而改变。
基于上述分析可以得出:对于净应力增大过程,系数ci|jT和di|jT随净应力变化过程数的增多而减小;对于净应力减小过程,系数ci|jT和di|jT随净应力变化过程数的增多基本不发生变化。
因此,采用幂函数回归分析建立净应力增大过程应力敏感指数与常压孔隙度的关系模型:
式中:为i号实验岩样净应力增大过程的孔隙度应力敏感指数;为i号实验岩样净应力增大过程的方程回归系数ci|jT随净应力变化周期变化的幂回归函数;/>为净应力增大过程的方程回归系数di|jT随净应力变化周期变化的幂回归函数;D为储气库运行周期。
S222将所述净应力减小过程的系数进行算术平均处理,建立净应力减小过程应力敏感指数的关系模型,得到净应力减小过程应力敏感指数关系模型;
具体的,将净应力减小过程的系数ci|jT和di|jT进行算术平均处理,建立净应力减小过程应力敏感指数关系模型:
式中:为i号实验岩样净应力减小过程的孔隙度应力敏感指数;/>为i号实验岩样净应力减小过程的方程回归系数/>的算术平均值;/>为i号实验岩样净应力减小过程的方程回归系数/>的算术平均值。
S223基于所述净应力增大过程应力敏感指数关系模型和所述净应力减小过程应力敏感指数关系模型得到第二拟合参数。
S23基于所述第一拟合参数和所述第二拟合参数建立所述实验岩样的不可逆孔隙度损失率模型。
具体的,建立实验岩样不可逆孔隙度损失率模型:
各净应力变化过程覆压孔隙度比值与无因次净应力关系曲线是连续的,净应力增大过程的终点即是下一个净应力减小过程的起点,净应力减小过程的终点又是下一个净应力增大过程的起点。第(2j-1)个净应力变化过程为净应力增大过程,第2j个净应力变化过程为净应力减小过程。根据覆压孔隙度比值与无因次净应力关系曲线的连续性,在最大无因次净应力下,净应力增大过程与下一个净应力减小过程的覆压孔隙度比值相等,即:
式中:为i号实验岩样第(2j-1)个净应力变化过程对应的孔隙度保持率;ai|2jD为i号实验岩样第2j个净应力变化过程对应的孔隙度保持率;bi|(2j-1)D为i号实验岩样第(2j-1)个净应力变化过程的孔隙度应力敏感指数;bi|2jD为i号实验岩样第2j个净应力变化过程的孔隙度应力敏感指数;/>为最大无因次净应力。
则第2j个净应力变化过程对应的孔隙度保持率为:
第1个净应力变化过程对应的孔隙度保持率始终为1,根据覆压孔隙度比值与无因次净应力关系曲线的连续性,可得第2j个净应力变化过程对应的孔隙度保持率为:
则第j个储气库运行周期的不可逆损失率为:
式中:为i号实验岩样第j个储气库运行周期的不可逆孔隙度损失率。
至此,根据不可逆孔隙度损失率计算模型和净应力增大过程的孔隙度应力敏感指数计算模型,可计算得到岩石相对于初始覆压孔隙度的不可逆孔隙度损失率。
S3基于所述不同物性实验岩样的岩石常压孔隙度和所述不同实验岩样的初始覆压孔隙度建立岩石初始覆压孔隙度计算模型;
具体的,拟合初始覆压孔隙度与常压孔隙度关系(初始覆压孔隙度计算模型);
目前,油气田储量计算和储气库库容计算一般采用的孔隙度为室内实验测试所得岩样孔隙度的算术平均或加权平均值,而室内实验多为常压孔渗实验,测试结果与实际储层条件下的孔隙度差距较大。因此,需要建立岩石初始覆压孔隙度(第1个净应力变化过程初始净应力下的覆压孔隙度)与岩石常压孔隙度之间的转换关系式:
根据岩石初始覆压孔隙度(第1个净应力变化过程初始净应力下的覆压孔隙度)与岩石常压孔隙度转换关系式,即可通过室内实验常压孔隙度计算得到储气库储集层岩石的初始覆压孔隙度。
S4基于所述不可逆孔隙度损失率模型和所述岩石初始覆压孔隙度计算模型建立初始净应力下的岩石覆压孔隙度计算模型;
储集层岩石在净应力增大时存在着一定的塑性变形,随着净应力变化周期数的增多,岩石累积塑性形变量增大,不可逆孔隙度损失率增大,储集层岩石不仅存在孔隙度应力敏感性,还存在一定的孔隙度交变应力敏感性。
无因次净应力为1时(也即初始净应力下)的覆压孔隙度计算模型为:
式中:为i号实验岩样在第j个储气库运行周期内无因次净应力为1时(也即初始净应力下)的覆压孔隙度。
S5基于所述初始净应力下的岩石覆压孔隙度计算模型建立储气库注采动态库容计算模型;
具体的,整理不可逆孔隙度损失率计算模型、岩石初始覆压孔隙度与岩石常压孔隙度转换关系以及无因次净应力为1时的覆压孔隙度计算模型,可得:
根据岩石室内实验常压孔隙度,采用上式可计算得到任意储气库运行周期初始净应力下的覆压孔隙度。
考虑孔隙度应力敏感性和孔隙度交变应力敏感性的储气库多周期注采动态库容评价模型为:
式中:G为储气库库容,m3;Ag为储气库储层面积,m2;h为储气库储层厚度,m;Swc为束缚水饱和度;Bgi为天然气体积系数。
S6将储气单元平均参数输入所述储气库注采动态库容计算模型得到特定储层孔隙度和特定储气库运行周期的储气库库容。
具体的,根据研究目标地下储气库的储气单元平均参数,采用建立的多周期注采动态库容计算模型,计算储气库在未来多个注采周期内的动态库容,研究其变化规律。
实施案例:
根据建立的油气藏型地下储气库多周期注采动态库容评价方法,评价了中国中部地区XX气藏型地下储气库在运行过程中动态库容的变化情况。
(1)开展不同物性实验岩样多周期循环载荷覆压孔隙体积实验
选取了6块标准岩心样品(表1)。采用岩石周期载荷覆压孔隙体积实验测试系统及测试方法,测定了6块岩样每块岩样7个净应力变化过程的覆压孔隙体积。
表1实验岩心基础数据表
(2)计算不同物性实验岩样各净应力状态下覆压孔隙度
根据实验测试结果,采用下式计算了覆压孔隙度,绘制覆压孔隙度与净应力关系曲线(图2)。
(3)拟合覆压孔隙度比值与无因次净应力关系
根据实验测试结果,将覆压孔隙度除以初始覆压孔隙度(第1个净应力变化过程初始净应力下的覆压孔隙度)得到覆压孔隙度比值,将净应力除以初始净应力得到无因次净应力,绘制覆压孔隙度比值与无因次净应力关系曲线(图3)。
通过下式拟合覆压孔隙度比值与无因次净应力关系,获取拟合参数(表2~表3)。
表2不可逆孔隙度损失率计算结果
表3孔隙度应力敏感指数回归分析结果
(4)拟合孔隙度应力敏感指数与常压孔隙度关系
绘制应力敏感指数与常压孔隙度关系曲线(图4),获取应力敏感指数与常压孔隙度拟合参数(表4)。
表4方程系数c和d回归分析结果统计表
对于净应力增大过程,系数ci|jT和di|jT随净应力变化过程数的增多而减小并趋于净应力减小过程的ci|jT和di|jT,则净应力增大过程系数ci|jT和di|jT与净应力减小过程系数ci|jT和di|jT的差值会随净应力变化过程数的增多而减小并趋于0。绘制系数ci|jT和di|jT差值与净应力变化过程关系曲线(图5),采用幂函数回归分析建立净应力增大过程应力敏感指数与常压孔隙度的关系模型:
对于净应力减小过程,将净应力减小过程的系数ci|jT和di|jT进行算术平均处理,建立净应力减小过程应力敏感指数关系模型:
(5)建立实验岩样不可逆孔隙度损失率模型
基于应力敏感指数关系模型建立岩样孔隙度保持率计算模型:
则第j个储气库运行周期的不可逆损失率为:
(6)拟合初始覆压孔隙度与常压孔隙度关系(初始覆压孔隙度计算模型)
作用在岩石上的净应力从0增大到初始净应力,岩石被压缩,岩石孔隙度减小,引入岩石孔隙度应力损失率,定义式为:
式中:为孔隙度应力损失率。
根据XX储气库6块岩心样品覆压孔隙度测试结果(图1),依据上式计算岩样孔隙度应力损失率(表5),根据计算结果分析孔隙度应力损失率与常压孔隙度间成幂函数关系(图6)。
表5孔隙度应力损失率计算结果
采用幂函数分析建立孔隙度损失率与常压孔隙度关系式为:
结合岩石孔隙度应力损失率定义式与孔隙度损失率与常压孔隙度关系式可得初始覆压孔隙度计算式为:
(7)计算初始净应力下岩样覆压孔隙度(初始净应力下的岩石覆压孔隙度计算模型)
根据无因次净应力为1时的覆压孔隙度计算模型,结合不可逆孔隙度损失率计算模型、岩石初始覆压孔隙度与岩石常压孔隙度转换关系,可得:
(8)建立多周期注采动态库容计算模型
建立考虑孔隙度应力敏感性和孔隙度交变应力敏感性的储气库多周期注采动态库容计算模型为:
(9)计算油气藏型地下储气库多周期注采动态库容
根据表6所示的XX储气库储层平均参数,采用考虑孔隙度应力敏感性和孔隙度交变应力敏感性的储气库多周期注采动态库容评价模型计算考虑孔隙度应力敏感和交变应力敏感综合影响的第30个储气库运行周期XX储气库储层A和储层B的库容变化情况,根据计算结果绘制库容随储气库运行周期变化曲线(图7~图8)。
表6 XX储气库库容计算单元基本参数
XX储气库储层A的原始库容为105.85×108m3,储层B的原始库容为8.46×108m3,总原始库容为114.31×108m3。当XX储气库运行至第30个周期时,储层A和储层B库容分别为94.40×108m3和7.47×108m3,总库容为101.87×108m3,总库容损失将达到12.44×108m3,总库容损失率达10.89%。孔隙度应力敏感性和周期应力敏感性对储气库库容影响较大,要准确预测储气库运行动态,需要考虑孔隙度应力敏感性和周期应力敏感性的影响。
以上所揭露的仅为本发明一种油气藏型地下储气库多周期注采动态库容计算方法较佳实施例而已,当然不能以此来限定本发明之权利范围,本领域普通技术人员可以理解实现上述实施例的全部或部分流程,并依本发明权利要求所作的等同变化,仍属于发明所涵盖的范围。
Claims (4)
1.一种油气藏型地下储气库多周期注采动态库容计算方法,其特征在于,包括以下步骤:
获取不同物性实验岩样的岩石常压孔隙度、不同实验净应力和所述不同实验净应力对应的覆压孔隙度;
基于所述不同物性实验岩样的岩石常压孔隙度、所述不同实验净应力和所述不同实验净应力对应的覆压孔隙度建立所述实验岩样的不可逆孔隙度损失率模型,包括,基于所述岩石不同净应力及其对应的覆压孔隙度拟合覆压孔隙度比值与无因次净应力关系,得到第一拟合参数;
基于所述岩石覆压孔隙度拟合孔隙度应力敏感指数与常压孔隙度关系,得到第二拟合参数;
基于所述第一拟合参数和所述第二拟合参数建立所述实验岩样的不可逆孔隙度损失率模型;
基于所述不同物性实验岩样的岩石常压孔隙度和不同实验岩样的初始覆压孔隙度建立岩石初始覆压孔隙度计算模型,包括,建立岩石初始覆压孔隙度与岩石常压孔隙度之间的转换关系式,根据岩石初始覆压孔隙度与岩石常压孔隙度转换关系式,通过常压孔隙度计算得到岩石的初始覆压孔隙度模型;
基于所述不可逆孔隙度损失率模型和所述岩石初始覆压孔隙度计算模型建立初始净应力下的岩石覆压孔隙度计算模型,计算模型为:
式中:为i号实验岩样在第j个储气库运行周期内初始净应力下的覆压孔隙度;
基于所述初始净应力下的岩石覆压孔隙度计算模型建立储气库注采动态库容计算模型,包括,整理不可逆孔隙度损失率计算模型、岩石初始覆压孔隙度与岩石常压孔隙度转换关系以及初始净应力下的覆压孔隙度计算模型,根据岩石室内实验常压孔隙度,计算得到任意储气库运行周期初始净应力下的覆压孔隙度,考虑孔隙度应力敏感性和孔隙度交变应力敏感性的储气库多周期注采动态库容评价模型为:
式中:G为储气库库容,m3;Ag为储气库储层面积,m2;h为储气库储层厚度,m;Swc为束缚水饱和度;Bgi为天然气体积系数;
将储气单元平均参数输入所述储气库注采动态库容计算模型得到研究目标储气库的储层孔隙度和储气库运行周期的储气库库容。
2.如权利要求1所述的油气藏型地下储气库多周期注采动态库容计算方法,其特征在于,
所述获取不同物性实验岩样的岩石常压孔隙度、所述不同实验净应力和所述不同实验净应力对应的覆压孔隙度的具体方式为:
对实验岩样开展常压孔隙度测试实验,得到常压孔隙度;
对所述实验岩样开展多周期循环载荷覆压孔隙体积实验,得到覆压孔隙体积;
基于所述覆压孔隙体积计算所述实验岩样在不同净应力状态下对应的岩石常压覆压孔隙度。
3.如权利要求1所述的油气藏型地下储气库多周期注采动态库容计算方法,其特征在于,
所述基于所述岩石覆压孔隙度拟合覆压孔隙度比值与无因次净应力关系,得到第一拟合参数的具体方式为:
将所述岩石覆压孔隙度除以初始覆压孔隙度,得到覆压孔隙度比值;
将净应力除以初始净应力得到无因次净应力;
基于所述覆压孔隙度比值和所述无因次净应力,得到第一拟合参数。
4.如权利要求3所述的油气藏型地下储气库多周期注采动态库容计算方法,其特征在于,
所述基于所述岩石常压覆压孔隙度拟合孔隙度应力敏感指数与常压孔隙度关系,得到第二拟合参数的具体方式为:
采用幂函数回归分析建立所述净应力增大过程应力敏感指数与常压孔隙度的关系模型,得到净应力增大过程应力敏感指数关系模型;
将所述净应力减小过程的系数进行算术平均处理,建立净应力减小过程应力敏感指数的关系模型,得到净应力减小过程应力敏感指数关系模型;
基于所述净应力增大过程应力敏感指数关系模型和所述净应力减小过程应力敏感指数关系模型得到第二拟合参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210580896.XA CN114970153B (zh) | 2022-05-25 | 2022-05-25 | 一种油气藏型地下储气库多周期注采动态库容计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210580896.XA CN114970153B (zh) | 2022-05-25 | 2022-05-25 | 一种油气藏型地下储气库多周期注采动态库容计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114970153A CN114970153A (zh) | 2022-08-30 |
CN114970153B true CN114970153B (zh) | 2023-07-25 |
Family
ID=82955816
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210580896.XA Active CN114970153B (zh) | 2022-05-25 | 2022-05-25 | 一种油气藏型地下储气库多周期注采动态库容计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114970153B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011133885A1 (en) * | 2010-04-23 | 2011-10-27 | The Board Of Regents Of The University Of Oklahoma | Total storage capacity and total porosity of porous media |
CN113011736A (zh) * | 2021-03-17 | 2021-06-22 | 西南石油大学 | 基于动态分析的获取油气井储量起算标准的方法及装置 |
CN114372352A (zh) * | 2021-12-23 | 2022-04-19 | 成都理工大学 | 渗流-温度双场耦合数值模拟预测复杂断块油藏储气库调峰能力方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2558838C1 (ru) * | 2014-07-02 | 2015-08-10 | Общество с ограниченной ответственностью "Научно-исследовательский институт природных газов и газовых технологий-Газпром ВНИИГАЗ" | Способ моделирования и оценки активного объема подземного хранилища газа в водоносных трещиновато-поровых структурах |
CN111287740B (zh) * | 2020-04-09 | 2023-06-30 | 中国石油天然气集团有限公司 | 基于真实应变下的异常高压气藏动态储量计算方法 |
CN112255159B (zh) * | 2020-10-10 | 2021-08-10 | 西南石油大学 | 一种非均质页岩储层水化作用能力的综合评价方法 |
CN114441401B (zh) * | 2020-10-30 | 2023-12-26 | 中国石油天然气股份有限公司 | 覆压孔隙度的确定方法、装置、计算机设备及存储介质 |
CN112763140B (zh) * | 2020-12-23 | 2021-08-27 | 重庆科技学院 | 一种油藏型储气库盖层的动态密封性评价方法 |
CN113008752B (zh) * | 2021-02-01 | 2022-02-15 | 重庆科技学院 | 一种油藏型储气库库容计算的有效孔隙体积确定方法 |
-
2022
- 2022-05-25 CN CN202210580896.XA patent/CN114970153B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011133885A1 (en) * | 2010-04-23 | 2011-10-27 | The Board Of Regents Of The University Of Oklahoma | Total storage capacity and total porosity of porous media |
CN113011736A (zh) * | 2021-03-17 | 2021-06-22 | 西南石油大学 | 基于动态分析的获取油气井储量起算标准的方法及装置 |
CN114372352A (zh) * | 2021-12-23 | 2022-04-19 | 成都理工大学 | 渗流-温度双场耦合数值模拟预测复杂断块油藏储气库调峰能力方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114970153A (zh) | 2022-08-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104141490B (zh) | 有水气藏单井水侵情况判断及气井产量控制方法和装置 | |
CN109374497B (zh) | 一种岩石微观孔隙结构测试方法 | |
CN103046914B (zh) | 一种低渗气藏水平井分段压裂效果判断方法 | |
CN109057786B (zh) | 用于非常规油气藏可采储量估算的自适应模型的确定方法 | |
CN103900942A (zh) | 基于高压压汞分析的储层微观孔喉参数的连续表征方法 | |
CN104612635A (zh) | 煤层群瓦斯联合抽采达标预评判方法 | |
RU2015113307A (ru) | Система и способ для определения количества удерживаемого углеводородного флюида | |
CN114279842B (zh) | 一种确定岩石裂纹起裂应力和损伤应力的方法及系统 | |
CN110878690B (zh) | 基于碳同位素演化确定页岩气井生产递减曲线的方法 | |
CN110909413A (zh) | 一种混凝土面板堆石坝水平位移预测精确方法及系统 | |
CN111157346A (zh) | 地下水封储油库水幕系统设计及修正的分析系统及方法 | |
CN112528455A (zh) | 吸附不饱和煤层气井储层压力与动态储量计算方法及系统 | |
CN114970153B (zh) | 一种油气藏型地下储气库多周期注采动态库容计算方法 | |
CN114580261A (zh) | 可采储量预测模型建立方法、可采储量预测方法和装置 | |
CN103675236A (zh) | 一种测定致密油临界充注孔喉半径门限的方法及装置 | |
CN112784454A (zh) | 基于蒙特卡罗法的地热资源量计算方法、系统、介质及终端 | |
CN110633848B (zh) | 煤层气井稳产时间的预测方法及装置 | |
CN110857625B (zh) | 碳酸盐岩有封闭水体缝洞单元储量获取方法、装置及设备 | |
CN112377184A (zh) | 串珠状缝洞型碳酸盐岩储层物性参数分析方法及装置 | |
CN110032777B (zh) | 一种基于不确定性分析的钻井液密度安全窗口估计方法 | |
CN112834404B (zh) | 页岩储层的产气能力的确定方法和装置 | |
CN113627069A (zh) | 针对缝洞型油藏油气井的试井动态产量评价方法及系统 | |
CN117189038A (zh) | 一种新型基于稳产期的致密气藏气井合理配产方法及系统 | |
CN114109324B (zh) | 关井时间的确定方法及装置 | |
CN113266328A (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 |