一种储气库合理采气速度计算方法

技术特征：
1.一种储气库合理采气速度计算方法，其特征在于，步骤如下：s1、根据现场情况给定井口压力p
wh
，在该井口压力条件下，任意给定偏差系数初值z0、合理产气速度初值q0、对应的井底流压初值p
wf0
；s2、在给定的q0、p
wf0
状态，计算以下各项参数：井筒静气柱平均温度无因次指数s、摩阻系数f、气流雷诺数re、井筒流速v、天然气体积系数b
g
、井筒静气柱平均压力天然气黏度μ
g
、天然气密度ρ
g
、天然气偏差系数z；s3、利用计算得到的参数和试井资料代入到以下两个公式中计算产气能力q
g
和对应的井底流压p
wf
：：其中，式中，p
r
为地层压力，e为绝对粗糙度，d为油管内径，m；a和b为二项式产能方程系数；s4、对比q
g
和q0，p
wf
和p
wf0
，判断q
g
与q0之间的绝对差值、p
wf
与p
wf0
之间的绝对差值是否同时满足小于0.00001，若不满足，则将q
g
和p
wf
代入步骤s2中重新计算各项参数；然后进行步骤s3计算出新的q
g1
和p
wf1
；再对比q
g1
和q
g
，p
wf1
和p
wf
，判断q
g1
与q
g
之间的绝对差值、p
wf1
与p
wf
之间的绝对差值是否同时满足小于0.00001；如此重复迭代，直至两个绝对差值均小于0.00001为止；s5、输出最后一次迭代计算出的q
gn
和p
wfn
，即井口压力p
wh
下的合理产气速度以及其对应的井底流压；s6、改变给定的井口压力p
wh
，重复步骤s1-s5，计算得到不同井口压力p
wh
条件下的合理产气速度以及其对应的井底流压；s7、根据地层压力p
r
与临界出砂压差δp计算得出该地层压力下的临界出砂流量q
sc
：s8、根据s7中计算的临界出砂流量q
sc
，结合筛选条件：合理产气速度不超过临界出砂流量q
sc
，从步骤s5和s6计算的合理产气速度中筛选出最佳的合理采气速度。2.如权利要求1所述的储气库合理采气速度计算方法，其特征在于，步骤s2中，井筒静
气柱平均温度计算公式如下：式中，t
wf
，t
wh
、分别为气井井底温度、井口温度、井筒静气柱平均温度，k。3.如权利要求2所述的储气库合理采气速度计算方法，其特征在于，步骤s2中，无因次指数s计算公式为：式中，γ
g
为天然气相对密度，无因次；h为井口到气层中部深度，m；为井筒静气柱平均偏差系数。4.如权利要求1所述的储气库合理采气速度计算方法，其特征在于，步骤s2中，井筒流速v计算公式为：速v计算公式为：式中，q
gc
为标准状态条件下的天然气流量，m3/d；b
g
为天然气体积系数，小数；m为井筒截面积，m2；z为天然气偏差系数；为井筒静气柱平均温度；为井筒静气柱平均压力。5.如权利要求1所述的储气库合理采气速度计算方法，其特征在于，步骤s2中，井筒静气柱平均压力计算公式为：式中，p
wf
、p
wh
分别为气井井底流压、井口静压，mpa。6.如权利要求1所述的储气库合理采气速度计算方法，其特征在于，步骤s2中，天然气偏差系数z按照cranmer方法迭代计算。7.如权利要求1所述的储气库合理采气速度计算方法，其特征在于，步骤s2中，天然气黏度μ
g
按照lee-gonzalez-eakin半经验法计算。8.如权利要求1所述的储气库合理采气速度计算方法，其特征在于，步骤s2中，天然气密度ρ
g
计算公式为：式中，γ
g
为天然气相对密度，无因次；为井筒静气柱平均偏差系数；为井筒静气柱平均温度；为井筒静气柱平均压力。9.如权利要求1所述的储气库合理采气速度计算方法，其特征在于，还包括步骤s9，根据s7中所计算出的临界出砂流量q
sc
与地层压力p
r
计算出相应的地层压力与流量下的井口压力p
whsc
：
式中，为井筒静气柱平均温度、s为无因次指数、f为摩阻系数、为井筒静气柱平均压力、p
wf
为井底流压、e为绝对粗糙度、d为油管内径，m、z为天然气偏差系数。

技术总结
本发明公开了一种储气库合理采气速度计算方法，步骤S1、给定井口压力P


技术研发人员：胡书勇 董旭 廖伟 张赟新 张佳轶 张季 刘鑫哲 胡洋
受保护的技术使用者：西南石油大学
技术研发日：2022.10.25
技术公布日：2023/1/13