1.一种自然递减气井的目标产气量预测方法,其特征在于:方法如下:步骤1:计算拟压力ψ(p)及压力转换表皮的初始值S p (i);通过下式(2)式中:ψ(p)为压力p对应的拟压力,MPa 2 /mPa·s;a为μ(p)Z(p)-p 2 曲线采用线性拟合得到的第一系数,无量纲;b为μ(p)Z(p)-p 2 曲线采用线性拟合得到的第二系数,无量纲;p为压力,MPa;p 0 为参考压力,MPa;μ(p)为压力p对应的气体粘度,mPa·s;Z(p)为压力p对应的偏差系数,无量纲;当p取值p i 、p wf (1)、p wf * 时,计算得到拟原始地层压力ψ(p i )、p wf (1)时的拟井底流压ψ(p wf (1))、p wf * 时的拟井底流压ψ(p wf * );计算压力转换表皮的初始值S p (i);通过下式(4)式中:S p (i)为压力转换表皮的初始值,无量纲;p sc 为标准压力,MPa;Z sc 为标准偏差系数,无量纲;T sc 为标准温度,K;T为储层温度,K;K为储层渗透率,10 -3 μm 2 ;h为储层厚度,m;ψ(p i )为拟原始地层压力,MPa 2 /mPa·s;p i 为原始地层压力,MPa;ψ(p wf (1))为p wf (1)时的拟井底流压,MPa 2 /mPa·s;p wf (1)为t=1时的井底流压,MPa;q sc (1)为t=1时的井口产气量,m 3 /d;计算得到压力转换表皮的初始值S p (i);步骤2:计算递减阶段开始时间t * ;(1)设置t取值范围为[1,t T ];其中,t为生产时间,d;t T 为目标生产时间,d;t * 为递减阶段开始时间,d;(2)采用物质平衡计算平均地层压力 生产时间t的取值从1逐一增大到目标生产时间t T ,通过下式(6)式中: 为t时的平均地层压力,MPa;Z r (t)为 对应的偏差系数,无量纲;Z i 为p i 对应的偏差系数,无量纲;G为气井控制储量,m 3 ;q sc * 为稳产阶段的产气量,m 3 /d;计算得到对应的 通过下式(11)式中:A为 曲线采用线性拟合得到的第一系数,无量纲;B为 曲线采用线性拟合得到的第二系数,无量纲;计算得到对应的平均地层压力 绘制物质平衡方程的 曲线;(3)采用产能方程计算平均地层压力 生产时间t的取值从1逐一增大到目标生产时间t T ,通过下式(7)式中: 为t时的拟平均地层压力,MPa 2 /mPa·s;ψ(p wf * )为p wf * 时的拟井底流压,MPa 2 /mPa·s;p wf * 为递减阶段的井底流压,MPa;t D 为无因次生产时间,t D =t/1,t D >0,无量纲;计算得到对应时的拟平均地层压力 令(3)式中的 则(3)式中的 通过下式(3)计算得到对应的平均地层压力 绘制产能方程的 曲线;(4)获取递减阶段开始时间t *将物质平衡方程的 曲线与产能方程的 曲线绘制到一张图上,两条曲线的交点对应的横坐标即为气井的递减阶段开始时间t * ;步骤3:判断目标生产时间t T 与递减阶段开始时间t * 的关系;若目标生产时间t T <递减阶段开始时间t * ,则该时刻目标生产时间t T 对应的目标产量q sc (t T )=稳产阶段的产气量q sc * ,计算结束。
2.根据权利要求1所述一种自然递减气井的目标产气量预测方法,其特征在于:若目标生产时间t T ≥递减阶段开始时间t * ,目标生产时间t T 位于递减阶段;还包括:步骤4:计算递减阶段的产气量q sc (t);(1)令生产时间t=递减阶段开始时间t * ;(2)计算t时的产气量q sc (t);①假设产气量的迭代初值q sc (t) B ,迭代初值q sc (t) B 的取值范围为[1,q sc * ];②通过下式(8)式中:q sc (t) B 为递减阶段产气量的迭代初值,m 3 /d;G p (t-1)为t-1时的累产气量,m 3 ;计算得到对应的 通过公式(11)计算得到对应的平均地层压力 ③令公式(2)中的p取值 通过公式(2),计算得到对应的拟平均地层压力 ④将拟平均地层压力 代入下式(9)式中:q sc (t) E 为递减阶段产气量的迭代终值,m 3 /d;计算得到产气量的迭代终值q sc (t) E ;⑤判断迭代初值q sc (t) B 与迭代终值q sc (t) E 之间的迭代误差是否满足下式(10)|q sc (t) B -q sc (t) E |/q sc (t) E ≤ε(q sc ) (10)式中:|q sc (t) B -q sc (t) E |/q sc (t) E 为迭代误差,无量纲;ε(q sc )为迭代计算的误差上限,0≤ε(q sc )≤10%,无量纲;若满足,令q sc (t)=q sc (t) E ,直接执行⑥;若不满足,令q sc (t) B =q sc (t) E ,重复①到⑤继续迭代;直至迭代误差满足公式(10),输出q sc (t)=q sc (t) E ,执行⑥;⑥判断生产时间t与目标生产时间t T 的关系;若生产时间t=目标生产时间t T ,输出q sc (t T )=q sc (t),计算结束;若生产时间t<目标生产时间t T ,执行(3);(3)将t增大1后作为新的t,即令t=t+1,通过步骤(2)进行迭代,计算该时刻的产气量q sc (t),直至该时刻的产气量q sc (t)满足迭代误差要求,同时满足生产时间t=目标生产时间t T ,输出q sc (t T )=q sc (t),计算结束。