[发明专利]一种油气藏型地下储气库多周期注采动态库容计算方法

权利要求书

1.一种油气藏型地下储气库多周期注采动态库容计算方法，其特征在于，包括以下步骤：

获取不同物性实验岩样的岩石常压孔隙度、不同实验净应力和所述不同实验净应力对应的覆压孔隙度；

基于所述不同物性实验岩样的岩石常压孔隙度、所述不同实验净应力和所述不同实验净应力对应的覆压孔隙度建立所述实验岩样的不可逆孔隙度损失率模型，包括，基于所述岩石不同净应力及其对应的覆压孔隙度拟合覆压孔隙度比值与无因次净应力关系，得到第一拟合参数；

基于所述岩石覆压孔隙度拟合孔隙度应力敏感指数与常压孔隙度关系，得到第二拟合参数；

基于所述第一拟合参数和所述第二拟合参数建立所述实验岩样的不可逆孔隙度损失率模型；

基于所述不同物性实验岩样的岩石常压孔隙度和不同实验岩样的初始覆压孔隙度建立岩石初始覆压孔隙度计算模型，包括，建立岩石初始覆压孔隙度与岩石常压孔隙度之间的转换关系式，根据岩石初始覆压孔隙度与岩石常压孔隙度转换关系式，通过常压孔隙度计算得到岩石的初始覆压孔隙度模型；

基于所述不可逆孔隙度损失率模型和所述岩石初始覆压孔隙度计算模型建立初始净应力下的岩石覆压孔隙度计算模型，计算模型为：

式中：为i号实验岩样在第j个储气库运行周期内初始净应力下的覆压孔隙度；

基于所述初始净应力下的岩石覆压孔隙度计算模型建立储气库注采动态库容计算模型，包括，整理不可逆孔隙度损失率计算模型、岩石初始覆压孔隙度与岩石常压孔隙度转换关系以及初始净应力下的覆压孔隙度计算模型，根据岩石室内实验常压孔隙度，计算得到任意储气库运行周期初始净应力下的覆压孔隙度，考虑孔隙度应力敏感性和孔隙度交变应力敏感性的储气库多周期注采动态库容评价模型为：

式中：G为储气库库容，m3；A_g为储气库储层面积，m2；h为储气库储层厚度，m；S_wc为束缚水饱和度；B_gi为天然气体积系数；

将储气单元平均参数输入所述储气库注采动态库容计算模型得到研究目标储气库的储层孔隙度和储气库运行周期的储气库库容。

2.如权利要求1所述的油气藏型地下储气库多周期注采动态库容计算方法，其特征在于，

所述获取不同物性实验岩样的岩石常压孔隙度、所述不同实验净应力和所述不同实验净应力对应的覆压孔隙度的具体方式为：

对实验岩样开展常压孔隙度测试实验，得到常压孔隙度；

对所述实验岩样开展多周期循环载荷覆压孔隙体积实验，得到覆压孔隙体积；

基于所述覆压孔隙体积计算所述实验岩样在不同净应力状态下对应的岩石常压覆压孔隙度。

3.如权利要求1所述的油气藏型地下储气库多周期注采动态库容计算方法，其特征在于，

所述基于所述岩石覆压孔隙度拟合覆压孔隙度比值与无因次净应力关系，得到第一拟合参数的具体方式为：

将所述岩石覆压孔隙度除以初始覆压孔隙度，得到覆压孔隙度比值；

将净应力除以初始净应力得到无因次净应力；

基于所述覆压孔隙度比值和所述无因次净应力，得到第一拟合参数。

4.如权利要求3所述的油气藏型地下储气库多周期注采动态库容计算方法，其特征在于，

所述基于所述岩石常压覆压孔隙度拟合孔隙度应力敏感指数与常压孔隙度关系，得到第二拟合参数的具体方式为：

采用幂函数回归分析建立所述净应力增大过程应力敏感指数与常压孔隙度的关系模型，得到净应力增大过程应力敏感指数关系模型；

将所述净应力减小过程的系数进行算术平均处理，建立净应力减小过程应力敏感指数的关系模型，得到净应力减小过程应力敏感指数关系模型；

基于所述净应力增大过程应力敏感指数关系模型和所述净应力减小过程应力敏感指数关系模型得到第二拟合参数。