[发明专利]一种强非均质性储层中产水气井的定产量递减预测方法

权利要求书

1.一种强非均质性储层中产水气井的定产量递减预测方法，其特征在于，该方法包括以下步骤：

S100、收集目标产水气井的原始地层压力p_i、井口输压p_t、点测静压数据p_j与点测静压相对应的累产气量G_pj、地层温度T_i、井口温度t、井筒产层中部深度h、井筒半径r_w、无阻流量q_AOF、目前的累产气量G_p、累产水量W_p、日产气量q_g、日产水量q_w、气样的相对密度γ_g、氮气摩尔分数y_N2、二氧化碳摩尔分数y_CO2、硫化氢摩尔分数y_H2S，水样的相对密度γ_w和氯化钠摩尔分数y_NaCl；

S200、基于每天的累产水量和每天的累产气量，得出水驱常数a、水驱常数b，并得到目标产水气井的甲型水驱公式；

S300、采用Blasingame图版法拟合得出目标产水气井的动态储量G，由点测静压相对应的累产气量除以目标产水气井的动态储量得出点测静压相对应的储量采出程度R_j；

S400、收集目标产水气井压力恢复试井数据，进行压力恢复试井分析，计算得出目标产水气井所在储层的非均质系数D，具体步骤为：第一、基于目标产水气井压力恢复试井得到的压力随试井的变化数据，采用双重介质模型进行数据拟合，得到弹性储容比ω与窜流系数λ；第二、根据拟合得到的弹性储容比ω与窜流系数λ，代入计算储层非均质系数D，其中，α为形状因子，由目标产水气井所在储层中的岩心取心获得，单位为m^-2；r_w为井筒半径，单位为m；λ为窜流系数，无单位；ω为弹性储容比，无单位；D为储层非均质系数，无单位；

S500、根据收集到的气体相对密度γ_g、原始地层压力p_i、点测静压数据p，采用D-A-K法求得原始地层压力下的偏差因子z_i、点测静压下的偏差因子z；

S600、结合水封气物质平衡公式，采用牛顿迭代法计算得出水侵常数C，其中，p为点测静压数据，单位为MPa；z为点测静压下的偏差因子，无单位；p_i为原始地层压力，单位为MPa；z_i为原始地层压力下的偏差因子，无单位；D为储层非均质系数，无单位；R为储量的采出程度，无单位；C为水侵常数，无单位，具体步骤为：第一、由水封气物质平衡公式得到水侵常数C为未知量的公式，其中，f（C）为表征水侵常数C的公式，无单位；第二、根据f（C）对水侵常数C进行求导，得到，其中f′（C）为f（C）对水侵常数C求导后的公式，无单位；第三、设定水侵常数C为1，带入f（C）与f′（C），通过C减去f（C）与f′（C）的比值计算得出新的水侵常数C₁；第四、计算C与C₁的绝对差，若C与C₁的绝对差小于0.00001，则C₁为所求目标产水气井的水侵常数；若C与C₁的绝对差大于0.00001，则将C替换为C₁，重新带入f（C）与f′（C），计算得出全新的水侵常数C₁，直至C与C₁的绝对差小于0.00001，最终得出目标气井的水侵常数C；

S700、对目标产水气井进行定产量的递减预测，得出目标产水气井定产量的稳产年限，具体步骤为：第一、采用Hagedom-Brown法，将原始地层压力p_i、井口输压p_t、地层温度T_i、井口温度t、井筒产层中部深度h、井筒半径r_w、日产气量q_g、日产水量q_w、气样的相对密度γ_g、氮气摩尔分数y_N2、二氧化碳摩尔分数y_CO2、硫化氢摩尔分数y_H2S，水样的相对密度γ_w和氯化钠摩尔分数y_NaCl带入，求得井口输压下的井底流压p_wfmin，即稳产末期的井底流压p_wfmin；第二、根据一点法公式计算稳产末期的井底流压p_wfmin下的稳产末期地层压力p_min；第三、由目标产水气井目前的累产气量除以目标产水气井的动态储量得出储量的采出程度R，结合水封气物质平衡公式与D-A-K法得出目前地层压力p与目前地层压力对应的压缩因子z；第四、目标产水气井以q_g定量生产，以时间1天为迭代步长，叠加G_p得到新一天的累产气量，代入目标产水气井的甲型水驱公式计算新一天的累产水量，结合水封气物质平衡公式与D-A-K法得到新一天的地层压力，直到新一天的地层压力小于等于稳产末期地层压力p_min结束迭代，代入一点法公式反算井底流压，绘制井底流压随时间的变化曲线，得到目标产水气井定产量递减预测曲线；第五、根据迭代结束时的时间除以365天得出目标产水气井的稳产年限。