一种气藏储量计算方法、装置、系统以及存储介质与流程

1.本发明主要涉及油气开发技术领域，具体涉及一种气藏储量计算方法、装置、系统以及存储介质。


背景技术：

2.地下储气库大多为枯竭型、弱边水气藏改建，在注采运行过程中，气体和水体始终接触，周期运行形成的气水过渡带上两相交互驱替的渗流特征，对储气库的注采能力和储气库储量计算造成了一定影响。目前气藏储量计算方法已较为成熟。sh i等采用物质平衡法拟合动态生产数据对气藏储量进行了计算。luo等建立了超深天然裂缝性致密砂岩气藏储量计算方法，并分析了影响动态储量评价主要影响因素。okoro等利用递减曲线对气藏储量进行了计算。
3.但现有气藏储量计算方法不能完全适用于产水储气库，la i等基于物质平衡原理，考虑了co2注采周期性，建立了枯竭湿/干气藏co2储量计算方法。但是，此方法采用的地层压力数据，在实际生产中很难获取。程时清等将气藏物质平衡方程与二项式产能方程相结合，联立求解气井井底流压，与实际生产过程中的井底流压曲线进行拟合，计算单井动态储量。此方法避免求解地层压力，充分考虑了井底压力的变化，但是忽略了产水情形，而且拟合方法比较耗时。徐兵祥等建立了产水气井物质平衡方程，并描述了水气两相物质平衡方程求解气井动态储量的流程。但是此方法只考虑了气藏产水情况，不能完全适用于多注采周期产水储气库动态储量评价。因此，如何结合产水和储气库注采周期特性计算更加准确的储量是一个亟需解决的难题。


技术实现要素：

4.本发明所要解决的技术问题是针对现有技术的不足，提供一种气藏储量计算方法、装置、系统以及存储介质。
5.本发明解决上述技术问题的技术方案如下：一种气藏储量计算方法，包括如下步骤：
6.s1：获得初始动态储量值和产能方程系数，并导入储气库气井生产数据以及储层参数；
7.s2：通过所述储气库气井生产数据、所述储层参数、所述初始动态储量值以及所述产能方程系数分析井口压力，得到预测井口压力；
8.s3：通过所述预测井口压力和所述储气库气井生产数据计算井口压力的误差值，根据计算结果得到目标动态储量值，并将所述目标动态储量值作为气藏储量计算结果。
9.本发明解决上述技术问题的另一技术方案如下：一种气藏储量计算装置，包括：
10.数据获得模块，用于获得初始动态储量值和产能方程系数，并导入储气库气井生产数据以及储层参数；
11.压力分析模块，用于通过所述储气库气井生产数据、所述储层参数、所述初始动态
储量值以及所述产能方程系数分析井口压力，得到预测井口压力；
12.计算结果获得模块，用于通过所述预测井口压力和所述储气库气井生产数据计算井口压力的误差值，根据计算结果得到目标动态储量值，并将所述目标动态储量值作为气藏储量计算结果。
13.基于上述一种气藏储量计算方法，本发明还提供一种气藏储量计算系统。
14.本发明解决上述技术问题的另一技术方案如下：一种气藏储量计算系统，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，当所述处理器执行所述计算机程序时，实现如上所述的气藏储量计算方法。
15.基于上述一种气藏储量计算方法，本发明还提供一种计算机可读存储介质。
16.本发明解决上述技术问题的另一技术方案如下：一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，当所述计算机程序被处理器执行时，实现如上所述的气藏储量计算方法。
17.本发明的有益效果是：通过储气库气井生产数据、储层参数、初始动态储量值以及产能方程系数分析井口压力得到预测井口压力，通过预测井口压力和储气库气井生产数据计算井口压力的误差值，根据计算结果得到目标动态储量值，并将目标动态储量值作为气藏储量计算结果，降低了对储气库动态储量预测的影响，解决了近似偏差因子带来的误差，具有较快的拟合速度，提高了动态储量的预测精度。
附图说明
18.图1为本发明实施例提供的气藏储量计算方法的流程示意图；
19.图2为本发明实施例提供的气藏储量计算装置的模块框图。
具体实施方式
20.以下结合附图对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。
21.图1为本发明实施例提供的一种气藏储量计算方法的流程示意图。
22.如图1所示，一种气藏储量计算方法，包括如下步骤：
23.s1：获得初始动态储量值和产能方程系数，并导入储气库气井生产数据以及储层参数；
24.s2：通过所述储气库气井生产数据、所述储层参数、所述初始动态储量值以及所述产能方程系数分析井口压力，得到预测井口压力；
25.s3：通过所述预测井口压力和所述储气库气井生产数据计算井口压力的误差值，根据计算结果得到目标动态储量值，并将所述目标动态储量值作为气藏储量计算结果。
26.应理解地，获取所述储气库气井生产数据及所述储层参数；设置所述初始动态储量值和所述产能方程系数。
27.上述实施例中，通过储气库气井生产数据、储层参数、初始动态储量值以及产能方程系数分析井口压力得到预测井口压力，通过预测井口压力和储气库气井生产数据计算井口压力的误差值，根据计算结果得到目标动态储量值，并将目标动态储量值作为气藏储量计算结果，降低了对储气库动态储量预测的影响，解决了近似偏差因子带来的误差，具有较
快的拟合速度，提高了动态储量的预测精度。
28.可选地，作为本发明的一个实施例，所述s2的过程包括：
29.通过所述储气库气井生产数据、所述储层参数、所述初始动态储量值以及所述产能方程系数计算井底压力，得到井底压力；
30.计算所述储气库气井生产数据与所述井底压力的井口压力，得到预测井口压力。
31.应理解地，地下储气库大多由枯竭型、弱边水气藏改建，与气藏不同的是储气库存在即注气又采气复杂生产模型，如果直接利用气藏储量评价方法直接用于储气库，计算结果可能存在误差。在注采运行过程中，气体和水体始终接触，周期运行形成的气水过渡带上两相交互驱替的渗流特征，对储气库的注采能力和储气库储量评价造成一定影响。因此，需要建立产水物质平衡方程。
32.上述实施例中，通过储气库气井生产数据、储层参数、初始动态储量值以及产能方程系数计算井底压力得到井底压力，计算储气库气井生产数据与井底压力的井口压力得到预测井口压力，降低了对储气库动态储量预测的影响，解决了近似偏差因子带来的误差，具有较快的拟合速度，提高了动态储量的预测精度。
33.可选地，作为本发明的一个实施例，所述储气库气井生产数据包括累积产水量、累积注气量、当前地层压力下偏差系数和日产气量，所述储层参数包括原始地层压力，所述产能方程系数包括指数产能方程系数和渗流指数，
34.所述通过所述储气库气井生产数据、所述储层参数、所述初始动态储量值以及所述产能方程系数计算井底压力，得到井底压力的过程包括：
35.通过第一式、所述累积产水量、所述累积注气量、所述日产气量、所述原始地层压力、所述初始动态储量值、所述渗流指数以及所述指数产能方程系数计算井底压力，得到井底压力，所述第一式为：
[0036][0037]
其中，p
wf
为井底压力，pi为原始地层压力，zi为预设地层压力下偏差系数，zr为当前地层压力下偏差系数，g
p
为累积产气量，g
in
为累积注气量，g为初始动态储量值，bw为水体积系数，w
p
为累积产水量，b
gi
为气体体积系数，qg为日产气量，c为指数产能方程系数，n为渗流指数。
[0038]
应理解地，利用储气库气井基本储层参数(即所述累积产水量、所述累积注气量、所述日产气量、所述原始地层压力、所述初始动态储量值、所述渗流指数以及所述指数产能方程系数)，计算所述井底压力，公式如下：
[0039]
[0040]
或
[0041][0042]
式中，w
p
为累积产水量，m3；b
gi
为气体在原始条件下的体积系数，m3/m3；bg为气体在某时刻的体积系数，m3/m3；bw为水在某时刻的体积系数，m3/m3；g
p
为累积产气量，104m3；g
in
为累积注气量，104m3；g为气井动态储量，104m3；pi为原始地层压力，mpa；zr为地层压力下偏差系数；zi为原始地层压力下偏差系数；a为二项式产能方程层流系数，mpa/(104m3·
d-1
)；b为二项式产能方程紊流系数，mpa/(104m3·
d-1
)2；c为指数产能方程系数，104(m3·
d-1
)/mpa
2n
；n为渗流指数；qg为日产气量，m3；p
wf
为井底流压，mpa。
[0043]
具体地，首先根据物质平衡原理，将储气库累积注气量和累积产气量因素考虑其中，得到气体物质平衡方程，公式如下：
[0044][0045]
和水方程：
[0046][0047]
式中，w
p
为累积产水量，m3；b
gi
为气体在原始条件下的体积系数，m3/m3；bg为气体在某时刻的体积系数，m3/m3；bw为水在某时刻的体积系数，m3/m3；s
wi
为原始含水饱和度；sw为某时刻含水饱和度；v为总的孔隙体积，104m3；g
p
为累积产气量，104m3；g
in
为累积注气量，104m3；g为气井动态储量，104m3。
[0048]
由式(2)得，公式如下：
[0049][0050]
将式(3)带入式(1)中得，公式如下：
[0051][0052]
其中，
[0053][0054][0055]
将公式(5)和(6)带入公式(4)中，公式(4)可化简为：
[0056][0057]
式中，pr为地层压力，mpa；pi为原始地层压力，mpa；zr为地层压力下偏差系数，zi为原始地层压力下偏差系数。
[0058]
s2：建立物质平衡-产能方程
[0059]
地层压力在实际施工过程中通常很难测量，为了避免用到地层压力数据，由物质平衡方程(7)导出如下关系式：
[0060][0061]
则地层压力可由式(8)右侧求解得到。但是，公式(8)中zr和bw是与pr相关的未知数。为了解决此问题，采用交替迭代算法对地层压力pr、zr和bw进行求解。具体求解步骤如下：
[0062]
step1：初始化p
r0
；
[0063]
step2：根据p
r0
，求解zr和bw；
[0064]
step3：将求解出的zr和bw带入公式(8)，计算pr；
[0065]
step4：求pr和p
r0
之间的绝对误差err0＝|p
r-p
r0
|；
[0066]
step5：如果err0＜err，退出循环；如果err0＞err，将pr赋值给p
r0
，重复step2，step3，step4。
[0067]
通过以上交替迭代算法即可求出zr和bw，然后将式(8)代入产能方程中。其中产能方程可分为二项式产能方程，公式如下：
[0068][0069]
和指数式产能方程，公式如下：
[0070][0071]
带入后得到如下公式：
[0072][0073]
或
[0074]
[0075]
其中，a为二项式产能方程层流系数，mpa/(104m3·
d-1
)；b为二项式产能方程紊流系数，mpa/(104m3·
d-1
)2；c为指数产能方程系数，104(m3·
d-1
)/mpa
2n
；n为渗流指数；qg为日产气量，m3；p
wf
为井底流压，mpa。
[0076]
上述实施例中，通过第一式、累积产水量、累积注气量、日产气量、原始地层压力、初始动态储量值、渗流指数以及指数产能方程系数计算井底压力得到井底压力,无需使用地层压力数据，降低了数据获得的难度，具有较快的拟合速度，提高了动态储量的预测精度。
[0077]
可选地，作为本发明的一个实施例，所述储气库气井生产数据包括日产气量，所述计算所述储气库气井生产数据与所述井底压力的井口压力，得到预测井口压力的过程包括：
[0078]
通过第二式计算与所述井底压力的井口压力，得到预测井口压力，所述第二式为：
[0079][0080]
其中，
[0081]
其中，ph为预测井口压力，p
wf
为井底压力，d为井筒内径，qg为日产气量，为井筒平均温度，s为无因次量，f为摩阻系数，为井筒平均偏差系数，rg为井筒内气体相对密度，h为井筒长度，re为雷诺数，ε为井筒内管道粗糙度。
[0082]
应理解地，利用井底压力，计算井口压力：
[0083][0084]
其中，
[0085][0086]
式中，ph为井口压力，mpa；为井筒平均温度，k；为井筒的平均偏差系数；h为井筒长度，m；rg为井筒内气体的相对密度；d为井筒内径，m；
[0087][0088]
式中，re为雷诺数；ε为注采井井筒内管道的粗糙度，m。
[0089]
具体地，得到井底流压，然后通过井口压力计算公式，将井底流压折算成井口压力，公式如下：
[0090][0091]
其中
[0092][0093]
ph为井口压力，mpa；为井筒平均温度，k；为井筒的平均偏差系数；h为井筒长度，m；rg为井筒内气体的相对密度；d为井筒内径，m；
[0094][0095]
式中，re为雷诺数；ε为注采井井筒内管道的粗糙度，m。
[0096]
式(11)和(12)中g
p
，g
in
，pi，w
p
和qg为已知量。zi可由ti和pi通过standing-katz图版计算得到。b
gi
可由地层温度t、zi和pi计算得到。然后，由式(13)得到的井口压力与实际井口压力拟合，可计算出公式(11)或(12)中的未知参数g，a和b或c和n。此方法即能评价产水储气库动态储量，又能确定气井二项式或指数式产能方程系数。
[0097]
上述实施例中，计算储气库气井生产数据与井底压力的井口压力得到预测井口压力，既能评价产水储气库动态储量，又能确定气井二项式或指数式产能方程系数，提高了动态储量的预测精度。
[0098]
可选地，作为本发明的一个实施例，
[0099]
所述储气库气井生产数据包括实际井口压力，所述s3中，通过所述预测井口压力和所述储气库气井生产数据计算井口压力的误差值，根据计算结果得到目标动态储量值的过程包括：
[0100]
s31：计算所述预测井口压力与所述实际井口压力之差，得到井口压力误差值；
[0101]
s32：判断所述井口压力误差值是否小于预设误差值，若否，则执行s33；若是，则将所述初始动态储量值作为目标动态储量值；
[0102]
s33：根据所述井口压力误差值分别对所述初始动态储量值和所述产能方程系数进行更新，更新后返回步骤s1。
[0103]
优选地，所述预设误差值可以为0.0001。
[0104]
具体地，求解计算的井口压力(即所所述预测井口压力)和实际生产井口压力(即所述实际井口压力)的绝对误差err(即所述井口压力误差)；当err《0.0001时，输出动态储量值(即所述初始动态储量值)和所述产能方程系数；否则，更新动态储量值(即所述初始动态储量值)和所述产能方程系数。
[0105]
上述实施例中，通过预测井口压力和储气库气井生产数据计算井口压力的误差值，根据计算结果得到目标动态储量值，降低了对储气库动态储量预测的影响，解决了近似偏差因子带来的误差，具有较快的拟合速度，提高了动态储量的预测精度。
[0106]
可选地，作为本发明的另一个实施例，物质平衡-产能方程是气井储量计算常用方
法之一，但是在求解过程中近似偏差因子和忽略气井产水影响，显著增加储量计算误差。此外，该方法还存在压力拟合过程耗时的缺点。为解决以上问题，本发明被提出。具体地，产水量这一因素和注采气累产量被引入物质平衡方程来降低其对储气库动态储量预测的影响。然后，交替迭代算法被用来求解物质平衡方程中的地层压力和偏差因子用于解决近似偏差因子带来的误差。最后，粒子群优化算法被用来拟合实际生产数据计算储气库动态储量，具有较快的拟合速度。实验表明了所提出方法可以有效提高拟合速度，并且对产水储气库动态储量评价更加合理和准确。
[0107]
可选地，作为本发明的另一个实施例，本发明将采用粒子群优化算法对所述预测井口压力和实际生产数据(即所述实际井口压力)进行拟合，得到储气库动态储量预测值。
[0108]
粒子群优化算法通过设计一种无质量的粒子来模拟鸟群中的鸟，粒子仅具有两个属性：速度v和位置x，速度代表移动的快慢，位置代表移动的方向。找到最优个体极值作为整个粒子群的当前全局最优解g
best
，粒子群中的所有粒子根据自己找到的当前个体极值p
best
和整个粒子群共享的当前全局最优解g
best
来调整自己的速度和位置，具体步骤如下描述：
[0109]
step1：初始化粒子群：设置最大速度区间，防止超出最大的区间。位置信息即为整个搜索空间，我们在速度区间和搜索空间上随机初始化速度和位置，设置群体规模。
[0110]
step2：个体极值与全局最优解：个体极值为每个粒子找到的历史上最优的位置信息，并从这些个体历史最优解中找到一个全局最优解，并与历史最优解比较，选出最佳的作为当前的历史最优解。
[0111]
step3：更新速度和位置：
[0112]vid
＝ω
×vid
+c1r1(p
id-x
id
)+c2r2(p
gd-x
id
)
ꢀꢀ
(16)
[0113]
x
id
＝x
id
+v
id
ꢀꢀ
(17)
[0114]
其中，ω为惯性权重，非负数，调节解空间的搜索范围；c1，c2为加速度常数，调节学习最大步长；r1，r2为[0,1]范围内的随机数。
[0115]
step4：终止条件：如果相邻两代之间的偏差在一个指定的范围内或达到虽大迭代次数即停止。
[0116]
此时种群的全局最优粒子，即为计算井口压力与实测井口压力达到最优拟合的动态储量和产能方程系数。
[0117]
可选地，作为本发明的另一个实施例，本发明首先，产水量和储气库注气和采气累积量被引入物质平衡方程来降低其对储气库动态储量预测的影响。然后，交替迭代算法被用来求解物质平衡方程中的地层压力和偏差因子解决近似偏差因子带来的误差。最后，粒子群优化算法被用来拟合实际生产数据计算储气库动态储量，具有较快的拟合速度。
[0118]
图2为本发明实施例提供的一种气藏储量计算装置的模块框图。
[0119]
可选地，作为本发明的另一个实施例，如图2所示，一种气藏储量计算装置，包括：
[0120]
数据获得模块，用于获得初始动态储量值和产能方程系数，并导入储气库气井生产数据以及储层参数；
[0121]
压力分析模块，用于通过所述储气库气井生产数据、所述储层参数、所述初始动态储量值以及所述产能方程系数分析井口压力，得到预测井口压力；
[0122]
计算结果获得模块，用于通过所述预测井口压力和所述储气库气井生产数据计算
井口压力的误差值，根据计算结果得到目标动态储量值，并将所述目标动态储量值作为气藏储量计算结果。
[0123]
可选地，作为本发明的一个实施例，所述压力分析模块具体用于：
[0124]
通过所述储气库气井生产数据、所述储层参数、所述初始动态储量值以及所述产能方程系数计算井底压力，得到井底压力；
[0125]
计算所述储气库气井生产数据与所述井底压力的井口压力，得到预测井口压力。
[0126]
可选地，作为本发明的一个实施例，所述储气库气井生产数据包括累积产水量、累积注气量、当前地层压力下偏差系数和日产气量，所述储层参数包括原始地层压力，所述产能方程系数包括指数产能方程系数和渗流指数，
[0127]
所述压力分析模块具体用于：
[0128]
通过第一式、所述累积产水量、所述累积注气量、所述日产气量、所述原始地层压力、所述初始动态储量值、所述渗流指数以及所述指数产能方程系数计算井底压力，得到井底压力，所述第一式为：
[0129][0130]
其中，p
wf
为井底压力，pi为原始地层压力，zi为预设地层压力下偏差系数，zr为当前地层压力下偏差系数，g
p
为累积产气量，g
in
为累积注气量，g为初始动态储量值，bw为水体积系数，w
p
为累积产水量，b
gi
为气体体积系数，qg为日产气量，c为指数产能方程系数，n为渗流指数。
[0131]
可选地，本发明的另一个实施例提供一种气藏储量计算系统，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，当所述处理器执行所述计算机程序时，实现如上所述的气藏储量计算方法。该系统可为计算机等系统。
[0132]
可选地，本发明的另一个实施例提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，当所述计算机程序被处理器执行时，实现如上所述的气藏储量计算方法。
[0133]
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。
[0134]
所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，上述描述的装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
[0135]
在本技术所提供的几个实施例中，应该理解到，所揭露的装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，单元的划分，仅仅为
一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。
[0136]
作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本发明实施例方案的目的。
[0137]
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以是两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
[0138]
集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分，或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例方法的全部或部分步骤。而前述的存储介质包括：u盘、移动硬盘、只读存储器(rom，read-on ly memory)、随机存取存储器(ram，random access memory)、磁碟或者光盘等各种可以存储程序代码的介质。
[0139]
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

技术特征：
1.一种气藏储量计算方法，其特征在于，包括如下步骤：s1：获得初始动态储量值和产能方程系数，并导入储气库气井生产数据以及储层参数；s2：通过所述储气库气井生产数据、所述储层参数、所述初始动态储量值以及所述产能方程系数分析井口压力，得到预测井口压力；s3：通过所述预测井口压力和所述储气库气井生产数据计算井口压力的误差值，根据计算结果得到目标动态储量值，并将所述目标动态储量值作为气藏储量计算结果。2.根据权利要求1所述的气藏储量计算方法，其特征在于，所述s2的过程包括：通过所述储气库气井生产数据、所述储层参数、所述初始动态储量值以及所述产能方程系数计算井底压力，得到井底压力；计算所述储气库气井生产数据与所述井底压力的井口压力，得到预测井口压力。3.根据权利要求2所述的气藏储量计算方法，其特征在于，所述储气库气井生产数据包括累积产水量、累积注气量、当前地层压力下偏差系数和日产气量，所述储层参数包括原始地层压力，所述产能方程系数包括指数产能方程系数和渗流指数，所述通过所述储气库气井生产数据、所述储层参数、所述初始动态储量值以及所述产能方程系数计算井底压力，得到井底压力的过程包括：通过第一式、所述累积产水量、所述累积注气量、所述日产气量、所述原始地层压力、所述初始动态储量值、所述渗流指数以及所述指数产能方程系数计算井底压力，得到井底压力，所述第一式为：其中，p
wf
为井底压力，p
i
为原始地层压力，z
i
为预设地层压力下偏差系数，z
r
为当前地层压力下偏差系数，g
p
为累积产气量，g
in
为累积注气量，g为初始动态储量值，b
w
为水体积系数，w
p
为累积产水量，b
gi
为气体体积系数，q
g
为日产气量，c为指数产能方程系数，n为渗流指数。4.根据权利要求2所述的气藏储量计算方法，其特征在于，所述储气库气井生产数据包括日产气量，所述计算所述储气库气井生产数据与所述井底压力的井口压力，得到预测井口压力的过程包括：通过第二式计算与所述井底压力的井口压力，得到预测井口压力，所述第二式为：其中，
其中，p
h
为预测井口压力，p
wf
为井底压力，d为井筒内径，q
g
为日产气量，为井筒平均温度，s为无因次量，f为摩阻系数，为井筒平均偏差系数，r
g
为井筒内气体相对密度，h为井筒长度，re为雷诺数，ε为井筒内管道粗糙度。5.根据权利要求1所述的气藏储量计算方法，其特征在于，所述储气库气井生产数据包括实际井口压力，所述s3中，通过所述预测井口压力和所述储气库气井生产数据计算井口压力的误差值，根据计算结果得到目标动态储量值的过程包括：s31：计算所述预测井口压力与所述实际井口压力之差，得到井口压力误差值；s32：判断所述井口压力误差值是否小于预设误差值，若否，则执行s33；若是，则将所述初始动态储量值作为目标动态储量值；s33：根据所述井口压力误差值分别对所述初始动态储量值和所述产能方程系数进行更新，更新后返回步骤s1。6.一种气藏储量计算装置，其特征在于，包括：数据获得模块，用于获得初始动态储量值和产能方程系数，并导入储气库气井生产数据以及储层参数；压力分析模块，用于通过所述储气库气井生产数据、所述储层参数、所述初始动态储量值以及所述产能方程系数分析井口压力，得到预测井口压力；计算结果获得模块，用于通过所述预测井口压力和所述储气库气井生产数据计算井口压力的误差值，根据计算结果得到目标动态储量值，并将所述目标动态储量值作为气藏储量计算结果。7.根据权利要求6所述的气藏储量计算装置，其特征在于，所述压力分析模块具体用于：通过所述储气库气井生产数据、所述储层参数、所述初始动态储量值以及所述产能方程系数计算井底压力，得到井底压力；计算所述储气库气井生产数据与所述井底压力的井口压力，得到预测井口压力。8.根据权利要求7所述的气藏储量计算装置，其特征在于，所述储气库气井生产数据包括累积产水量、累积注气量、当前地层压力下偏差系数和日产气量，所述储层参数包括原始地层压力，所述产能方程系数包括指数产能方程系数和渗流指数，所述压力分析模块具体用于：通过第一式、所述累积产水量、所述累积注气量、所述日产气量、所述原始地层压力、所述初始动态储量值、所述渗流指数以及所述指数产能方程系数计算井底压力，得到井底压力，所述第一式为：其中，p
wf
为井底压力，p
i
为原始地层压力，z
i
为预设地层压力下偏差系数，z
r
为当前地层
压力下偏差系数，g
p
为累积产气量，g
in
为累积注气量，g为初始动态储量值，b
w
为水体积系数，w
p
为累积产水量，b
gi
为气体体积系数，q
g
为日产气量，c为指数产能方程系数，n为渗流指数。9.一种气藏储量计算系统，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征在于，当所述处理器执行所述计算机程序时，实现如权利要求1至5任一项所述的气藏储量计算方法。10.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，当所述计算机程序被处理器执行时，实现如权利要求1至5任一项所述的气藏储量计算方法。

技术总结
本发明提供一种气藏储量计算方法、装置、系统以及存储介质，属于油气开发领域，方法包括：获得初始动态储量值和产能方程系数，并导入储气库气井生产数据以及储层参数；通过储气库气井生产数据、储层参数、初始动态储量值以及产能方程系数分析井口压力得到预测井口压力；通过预测井口压力和储气库气井生产数据计算井口压力的误差值，根据计算结果得到目标动态储量值，并将目标动态储量值作为气藏储量计算结果。本发明降低了对储气库动态储量预测的影响，解决了近似偏差因子带来的误差，具有较快的拟合速度，提高了动态储量的预测精度。提高了动态储量的预测精度。提高了动态储量的预测精度。


技术研发人员：任众鑫 王海峰 王多才 赵廉斌 刘峻峰 张宏 张栩赫 鲁俊 韩娜 段冲
受保护的技术使用者：国家石油天然气管网集团有限公司
技术研发日：2023.05.23
技术公布日：2023/8/14