CN113236207B - 一种强非均质性储层中产水气井的定产量递减预测方法 - Google Patents

一种强非均质性储层中产水气井的定产量递减预测方法 Download PDF

Info

Publication number
CN113236207B
CN113236207B CN202110787098.XA CN202110787098A CN113236207B CN 113236207 B CN113236207 B CN 113236207B CN 202110787098 A CN202110787098 A CN 202110787098A CN 113236207 B CN113236207 B CN 113236207B
Authority
CN
China
Prior art keywords
water
gas well
pressure
well
producing gas
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110787098.XA
Other languages
English (en)
Other versions
CN113236207A (zh
Inventor
谭晓华
崔苗逢
李晓平
漆麟
万翠蓉
孟展
丁磊
彭港珍
邓永建
韩晓冰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN202110787098.XA priority Critical patent/CN113236207B/zh
Publication of CN113236207A publication Critical patent/CN113236207A/zh
Application granted granted Critical
Publication of CN113236207B publication Critical patent/CN113236207B/zh
Priority to US17/812,106 priority patent/US20230017069A1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
    • E21B43/20—Displacing by water
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • E21B49/08—Obtaining fluid samples or testing fluids, in boreholes or wells
    • E21B49/087—Well testing, e.g. testing for reservoir productivity or formation parameters
    • E21B49/0875—Well testing, e.g. testing for reservoir productivity or formation parameters determining specific fluid parameters
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B2200/00—Special features related to earth drilling for obtaining oil, gas or water
    • E21B2200/20—Computer models or simulations, e.g. for reservoirs under production, drill bits
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B2200/00—Special features related to earth drilling for obtaining oil, gas or water
    • E21B2200/22—Fuzzy logic, artificial intelligence, neural networks or the like
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00—Survey of boreholes or wells
    • E21B47/06—Measuring temperature or pressure
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00—Survey of boreholes or wells
    • E21B47/06—Measuring temperature or pressure
    • E21B47/07—Temperature

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
  • Measuring Fluid Pressure (AREA)

Abstract

本发明为一种强非均质性储层中产水气井的定产量递减预测方法，它主要是收集目标产水气井相关数据，拟合得水驱常数与水侵常数，采用Blasingame图版法拟合动态储量，根据双重介质模型拟合弹性储容比与窜流系数，计算储层的非均质系数，得到稳产末期的井底流压，以时间1天为迭代步长，计算目标产水气井定量生产新一天的地层压力，直到该地层压力小于等于稳产末期地层压力结束迭代，绘制目标产水气井定产量递减预测曲线。本发明能结合试井分析定量评价储层非均质性，针对强非均质性产水气井的定量生产，对其进行生产预测，得出气井稳产年限，实现强非均质储层产水气井定产量递减预测。

Description

一种强非均质性储层中产水气井的定产量递减预测方法
技术领域
本发明属于气藏工程开发领域，具体涉及一种强非均质性储层中产水气井的定产量递减预测方法。
背景技术
天然气是较为安全的燃气之一，采用天然气作为能源，不仅可减少煤和石油的用量，从而大大改善环境污染问题；还能能减少二氧化碳、二氧化硫和粉尘的排放量，有助于减少酸雨形成，从根本上改善环境质量，因此天然气的开发显得尤为重要。但是对于强非均质性储层产水气井，储层的强非均质性以及产水都会对引起气井生产上的变化，导致气井的递减不同于常规气井。因此如何评价储层不均匀以及产水时气井的递减趋势，在气藏开发中是一个难题。
目前，已有发明专利CN201410638125.7《一种油气井产量递减分析方法和系统》提供了一种可应用于页岩气井和其他类型油气井产量分析和动态评价的油气井产量递减分析方法和系统，但是此方法为预测气井的产量递减，但实际气井生产过程往往为定产量生产，且并未针对非均质以及产水气井进行生产预测；发明专利CN201310314083.7《一种缝洞型碳酸盐岩气藏动态分析方法及系统》可以针对缝洞型非均质气藏进行定产量生产动态预测，并不能针对产水气井进行生产动态预测。因此，为了更好的形成一种强非均质性储层产水气井定产量递减预测方法，本发明针对强非均质性储层产水气井，进行了定产量生产的动态预测。
发明内容
本发明目的是针对强非均质产水气藏，建立定产量生产情况下气井的生产动态预测，形成强非均质性储层产水气井定产量递减预测方法，为气藏开发奠定基础。
本发明所采用的技术方案是：
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，具体步骤为：第一、基于目标产水气井压力恢复试井得到的压力随试井的变化数据，采用双重介质模型进行数据拟合，得到弹性储容比ω与窜流系数λ；第二、根据拟合得到的弹性储容比ω与窜流系数λ，代入
Figure 100002_DEST_PATH_IMAGE002
计算储层非均质系数D，其中，α为形状因子，由目标产水气井所在储层中的岩心取心获得，单位为m-2；r w为井筒半径，单位为m；λ为窜流系数，无单位；ω为弹性储容比，无单位；D为储层非均质系数，无单位；
S500、根据收集到的气体相对密度γ g、原始地层压力p i、点测静压数据p，采用D-A-K法求得原始地层压力下的偏差因子z i、点测静压下的偏差因子z；
S600、结合水封气物质平衡公式
Figure 100002_DEST_PATH_IMAGE004
，采用牛顿迭代法计算得出水侵常数C，其中，p为点测静压数据，单位为MPa；z为点测静压下的偏差因子，无单位；p i为原始地层压力，单位为MPa；z i为原始地层压力下的偏差因子，无单位；D为储层非均质系数，无单位；R为储量的采出程度，无单位；C为水侵常数，无单位，具体步骤为：第一、由水封气物质平衡公式得到水侵常数C为未知量的公式
Figure 100002_DEST_PATH_IMAGE006
，其中，f（C）为表征水侵常数C的公式，无单位；第二、根据f（C）对水侵常数C进行求导，得到
Figure 100002_DEST_PATH_IMAGE008
，其中f′（C）为f（C）对水侵常数C求导后的公式，无单位；第三、设定水侵常数C为1，带入f（C）与f′（C），通过C减去f（C）与f′（C）的比值计算得出新的水侵常数C 1；第四、计算C与C 1的绝对差，若C与C 1的绝对差小于0.00001，则C 1为所求目标产水气井的水侵常数；若C与C 1的绝对差大于0.00001，则将C替换为C 1，重新带入f（C）与f′（C），计算得出全新的水侵常数C 1，直至C与C 1的绝对差小于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天得出目标产水气井的稳产年限。
上述一种强非均质性储层中产水气井的定产量递减预测方法中，其特征在于，Blasingame图版法为采用RTA软件，输入目标产水气井的生产数据、原始地层压力、地层温度、井筒产层中部深度、井筒半径，在理论曲线图版上对实际生成曲线进行拟合，再由RTA软件自动计算出目标产水气井的动态储量。
上述一种强非均质性储层中产水气井的定产量递减预测方法中，其特征在于，D-A-K法为根据气样的相对密度γ g采用经验公式
Figure 100002_DEST_PATH_IMAGE010
与
Figure DEST_PATH_IMAGE012
计算拟临界压力p pc与拟临界温度T pc，再根据某一地层压力p k与某一地层温度T k采用
Figure DEST_PATH_IMAGE014
与
Figure DEST_PATH_IMAGE016
计算拟对比压力p pr与拟对比温度T pr，采用
Figure DEST_PATH_IMAGE018
、
Figure DEST_PATH_IMAGE020
与
Figure DEST_PATH_IMAGE022
这三个公式联立迭代计算出偏差因子，其中，γ g为气样的相对密度，无单位；p pc为拟临界压力，单位为MPa；T pc为拟临界温度，单位为K；p k为某一地层压力，单位为MPa；T k为某一地层温度，单位为K；p pr为拟对比压力，单位为MPa；T pr为拟对比温度，单位为K；ρ pr为拟对比密度，无单位；z k为某一地层压力对应的偏差因子，无单位；F（ρ pr）为表征拟对比密度的公式，无单位；A 1=0.3265，无单位；A 2=-1.0700，无单位；A 3=-0.5339，无单位；A 4=0.01569，无单位；A 5=-0.05165，无单位；A 6=0.5475，无单位；A 7=-0.7361，无单位；A 8=0.1844，无单位；A 9=0.1056，无单位；A 10=0.6134，无单位；A 11=0.7210，无单位；F′（ρ pr）为F（ρ pr）对ρ pr求导后的公式，无单位。
上述一种强非均质性储层中产水气井的定产量递减预测方法中，其特征在于，Hagedom-Brown法为
Figure DEST_PATH_IMAGE024
、
Figure DEST_PATH_IMAGE026
、
Figure DEST_PATH_IMAGE028
、
Figure DEST_PATH_IMAGE030
、
Figure DEST_PATH_IMAGE032
、
Figure DEST_PATH_IMAGE034
与
Figure DEST_PATH_IMAGE036
结合计算井底流压，其中，Δp为井管道压力增加量，单位为MPa；ΔH为井管道深度增加量，单位为m；ρ m为气水混合物密度，单位为kg/m3；g为重力加速度，单位为m/s2；f m为两相摩阻系数，无单位；G mA为单位井管道横截面积下的混合物质量流量，kg/s/m2；r w为井筒半径，单位为m；ρ w为目标产水气井水密度，由物性分析可知，单位为kg/m3；ρ g为目标产水气井气密度，由物性分析可知，单位为kg/m3；H L为持液率，无单位；e为管壁绝对粗糙度，管壁分析可知，单位为m；N Re为两相雷诺数，无单位；q g为日产气量，单位为m3；q w为日产水量，单位为m3；V sl为液相表观速度，单位为m/s；V sg为气相表观速度，单位为m/s；μ w为水黏度，由物性分析可知，单位为mPa·s；μ g为气黏度，由物性分析可知，单位为mPa·s。
上述一种强非均质性储层中产水气井的定产量递减预测方法中，其特征在于，一点法公式为
Figure DEST_PATH_IMAGE038
，其中，q g为日产气量，单位为m3；q AOF为无阻流量，单位为m3；p min为稳产末期地层压力，单位为MPa；p wfmin为稳产末期的井底流压，单位为MPa；
本发明的优点：能结合试井分析定量评价储层非均质性，针对强非均质性产水气井的定量生产，对其进行生产预测，得出气井稳产年限，实现强非均质储层产水气井定产量递减预测。
附图说明
在附图中：
图1是一种强非均质性储层中产水气井的定产量递减预测方法步骤图。
图2是某强非均质产水气井甲型水驱曲线图。
图3是某强非均质产水气井Blasingame图版拟合图。
图4是某强非均质产水气井双重介质拟合图。
图5是某强非均质产水气井定产量递减预测曲线。
具体实施方式
下面结合附图对本发明做进一步说明。
本发明提供了一种强非均质性储层中产水气井的定产量递减预测方法，图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，具体步骤为：第一、基于目标产水气井压力恢复试井得到的压力随试井的变化数据，采用双重介质模型进行数据拟合，得到弹性储容比ω与窜流系数λ；第二、根据拟合得到的弹性储容比ω与窜流系数λ，代入
Figure 650829DEST_PATH_IMAGE002
计算储层非均质系数D，其中，α为形状因子，由目标产水气井所在储层中的岩心取心获得，单位为m-2；r w为井筒半径，单位为m；λ为窜流系数，无单位；ω为弹性储容比，无单位；D为储层非均质系数，无单位；
S500、根据收集到的气体相对密度γ g、原始地层压力p i、点测静压数据p，采用D-A-K法求得原始地层压力下的偏差因子z i、点测静压下的偏差因子z；
S600、结合水封气物质平衡公式
Figure 567970DEST_PATH_IMAGE004
，采用牛顿迭代法计算得出水侵常数C，其中，p为点测静压数据，单位为MPa；z为点测静压下的偏差因子，无单位；p i为原始地层压力，单位为MPa；z i为原始地层压力下的偏差因子，无单位；D为储层非均质系数，无单位；R为储量的采出程度，无单位；C为水侵常数，无单位，具体步骤为：第一、由水封气物质平衡公式得到水侵常数C为未知量的公式
Figure 410024DEST_PATH_IMAGE006
，其中，f（C）为表征水侵常数C的公式，无单位；第二、根据f（C）对水侵常数C进行求导，得到
Figure 867550DEST_PATH_IMAGE008
，其中f′（C）为f（C）对水侵常数C求导后的公式，无单位；第三、设定水侵常数C为1，带入f（C）与f′（C），通过C减去f（C）与f′（C）的比值计算得出新的水侵常数C 1；第四、计算C与C 1的绝对差，若C与C 1的绝对差小于0.00001，则C 1为所求目标产水气井的水侵常数；若C与C 1的绝对差大于0.00001，则将C替换为C 1，重新带入f（C）与f′（C），计算得出全新的水侵常数C 1，直至C与C 1的绝对差小于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天得出目标产水气井的稳产年限。
进一步的，所述一种强非均质性储层中产水气井的定产量递减预测方法中，Blasingame图版法为采用RTA软件，输入目标产水气井的生产数据、原始地层压力、地层温度、井筒产层中部深度、井筒半径，在理论曲线图版上对实际生成曲线进行拟合，再由RTA软件自动计算出目标产水气井的动态储量。
进一步的，所述一种强非均质性储层中产水气井的定产量递减预测方法中，D-A-K法为根据气样的相对密度γ g采用经验公式
Figure 245704DEST_PATH_IMAGE010
与
Figure 130484DEST_PATH_IMAGE012
计算拟临界压力p pc与拟临界温度T pc，再根据某一地层压力p k与某一地层温度T k采用
Figure 397517DEST_PATH_IMAGE014
与
Figure 658734DEST_PATH_IMAGE016
计算拟对比压力p pr与拟对比温度T pr，采用
Figure 389929DEST_PATH_IMAGE018
、
Figure 944145DEST_PATH_IMAGE020
与
Figure 698475DEST_PATH_IMAGE022
这三个公式联立迭代计算出偏差因子，其中，γ g为气样的相对密度，无单位；p pc为拟临界压力，单位为MPa；T pc为拟临界温度，单位为K；p k为某一地层压力，单位为MPa；T k为某一地层温度，单位为K；p pr为拟对比压力，单位为MPa；T pr为拟对比温度，单位为K；ρ pr为拟对比密度，无单位；z k为某一地层压力对应的偏差因子，无单位；F（ρ pr）为表征拟对比密度的公式，无单位；A 1=0.3265，无单位；A 2=-1.0700，无单位；A 3=-0.5339，无单位；A 4=0.01569，无单位；A 5=-0.05165，无单位；A 6=0.5475，无单位；A 7=-0.7361，无单位；A 8=0.1844，无单位；A 9=0.1056，无单位；A 10=0.6134，无单位；A 11=0.7210，无单位；F′（ρ pr）为F（ρ pr）对ρ pr求导后的公式，无单位。
进一步的，所述一种强非均质性储层中产水气井的定产量递减预测方法中，Hagedom-Brown法为
Figure 497803DEST_PATH_IMAGE024
、
Figure 83506DEST_PATH_IMAGE026
、
Figure 513350DEST_PATH_IMAGE028
、
Figure 817292DEST_PATH_IMAGE030
、
Figure 718514DEST_PATH_IMAGE032
、
Figure 424302DEST_PATH_IMAGE034
与
Figure 759469DEST_PATH_IMAGE036
结合计算井底流压，其中，Δp为井管道压力增加量，单位为MPa；ΔH为井管道深度增加量，单位为m；ρ m为气水混合物密度，单位为kg/m3；g为重力加速度，单位为m/s2；f m为两相摩阻系数，无单位；G mA为单位井管道横截面积下的混合物质量流量，kg/s/m2；r w为井筒半径，单位为m；ρ w为目标产水气井水密度，由物性分析可知，单位为kg/m3；ρ g为目标产水气井气密度，由物性分析可知，单位为kg/m3；H L为持液率，无单位；e为管壁绝对粗糙度，管壁分析可知，单位为m；N Re为两相雷诺数，无单位；q g为日产气量，单位为m3；q w为日产水量，单位为m3；V sl为液相表观速度，单位为m/s；V sg为气相表观速度，单位为m/s；μ w为水黏度，由物性分析可知，单位为mPa·s；μ g为气黏度，由物性分析可知，单位为mPa·s。
进一步的，所述一种强非均质性储层中产水气井的定产量递减预测方法中，一点法公式为
Figure 550707DEST_PATH_IMAGE038
，其中，q g为日产气量，单位为m3；q AOF为无阻流量，单位为m3；p min为稳产末期地层压力，单位为MPa；p wfmin为稳产末期的井底流压，单位为MPa。
利用一种强非均质性储层中产水气井的定产量递减预测方法的步骤说明，以某一强非均质产水气井为例，进行气井的定产量生产动态预测，以确定气井的稳产年限。
收集该强非均质产水气井生产数据、物性分析数据、储层数据，基于甲型水驱公式拟合得到水驱常数a与b，得出a=4.948，b=0.000000046，如图2所示；采用Blasingame图版法拟合该强非均质产水气井动态储量，为518000000m3，如图3所示；再基于恢复试井，采用双重介质模型进行数据拟合，得到弹性储容比为0.223，窜流系数为0.00000111，计算出储层的非均质系数为3.4842，如图4所示；采用D-A-K法求得原始地层压力下的偏差因子为1.71；结合水封气物质平衡公式，采用牛顿迭代法计算得出水侵常数为2。再采用Hagedom-Brown法，求取井口输压下的井底流压为36.918MPa，根据一点法公式计算稳产末期的井底流压下的稳产末期地层压力为44.862MPa；通过时间步长的迭代，得出井底流压随时间的变化曲线，最终得到目标产水气井定产量递减预测曲线，如图5所示；根据定产量递减预测曲线，得出该强非均质产水气井的稳产年限为1.115年。
与现有气井定产量预测方法相比，本发明具有以下有益效果：能结合试井分析定量评价储层非均质性，针对强非均质性产水气井的定量生产，对其进行生产预测，得出气井稳产年限，实现强非均质储层产水气井定产量递减预测。
最后所应说明的是：以上实施例仅用以说明而非限制本发明的技术方案，尽管参照上述实施例对本发明进行了详细说明，本领域的普通技术人员应该理解：依然可以对本发明进行修改或者等同替换，而不脱离本发明的精神和范围的任何修改或局部替换，其均应涵盖在本发明的权利要求范围当中。

Claims (3)

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，具体步骤为：第一、基于目标产水气井压力恢复试井得到的压力随试井的变化数据，采用双重介质模型进行数据拟合，得到弹性储容比ω与窜流系数λ；第二、根据拟合得到的弹性储容比ω与窜流系数λ，代入
Figure DEST_PATH_IMAGE002
计算储层非均质系数D，其中，α为形状因子，由目标产水气井所在储层中的岩心取心获得，单位为m-2；r w为井筒半径，单位为m；λ为窜流系数，无单位；ω为弹性储容比，无单位；D为储层非均质系数，无单位；
S500、根据收集到的气体相对密度γ g、原始地层压力p i、点测静压数据p，采用D-A-K法求得原始地层压力下的偏差因子z i、点测静压下的偏差因子z；
S600、结合水封气物质平衡公式
Figure DEST_PATH_IMAGE004
，采用牛顿迭代法计算得出水侵常数C，其中，p为点测静压数据，单位为MPa；z为点测静压下的偏差因子，无单位；p i为原始地层压力，单位为MPa；z i为原始地层压力下的偏差因子，无单位；D为储层非均质系数，无单位；R为储量的采出程度，无单位；C为水侵常数，无单位，具体步骤为：第一、由水封气物质平衡公式得到水侵常数C为未知量的公式
Figure DEST_PATH_IMAGE006
，其中，f（C）为表征水侵常数C的公式，无单位；第二、根据f（C）对水侵常数C进行求导，得到
Figure DEST_PATH_IMAGE008
，其中f′（C）为f（C）对水侵常数C求导后的公式，无单位；第三、设定水侵常数C为1，带入f（C）与f′（C），通过C减去f（C）与f′（C）的比值计算得出新的水侵常数C 1；第四、计算C与C 1的绝对差，若C与C 1的绝对差小于0.00001，则C 1为所求目标产水气井的水侵常数；若C与C 1的绝对差大于0.00001，则将C替换为C 1，重新带入f（C）与f′（C），计算得出全新的水侵常数C 1，直至C与C 1的绝对差小于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天得出目标产水气井的稳产年限。
2.根据权利要求1所述的一种强非均质性储层中产水气井的定产量递减预测方法，其特征在于，步骤S300中所述Blasingame图版法为采用RTA软件，输入目标产水气井的生产数据、原始地层压力、地层温度、井筒产层中部深度、井筒半径，在理论曲线图版上对实际生成曲线进行拟合，再由RTA软件自动计算出目标产水气井的动态储量。
3.根据权利要求1所述的一种强非均质性储层中产水气井的定产量递减预测方法，其特征在于，步骤S700中所述一点法公式为
Figure DEST_PATH_IMAGE010
，其中，q g为日产气量，单位为m3；q AOF为无阻流量，单位为m3；p min为稳产末期地层压力，单位为MPa；p wfmin为稳产末期的井底流压，单位为MPa。
CN202110787098.XA 2021-07-13 2021-07-13 一种强非均质性储层中产水气井的定产量递减预测方法 Active CN113236207B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202110787098.XA CN113236207B (zh) 2021-07-13 2021-07-13 一种强非均质性储层中产水气井的定产量递减预测方法
US17/812,106 US20230017069A1 (en) 2021-07-13 2022-07-12 Prediction method for constant production decline of water-producing gas well in highly heterogeneous reservoir

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110787098.XA CN113236207B (zh) 2021-07-13 2021-07-13 一种强非均质性储层中产水气井的定产量递减预测方法

Publications (2)

Publication Number Publication Date
CN113236207A CN113236207A (zh) 2021-08-10
CN113236207B true CN113236207B (zh) 2021-09-10

Family

ID=77135301

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110787098.XA Active CN113236207B (zh) 2021-07-13 2021-07-13 一种强非均质性储层中产水气井的定产量递减预测方法

Country Status (2)

Country Link
US (1) US20230017069A1 (zh)
CN (1) CN113236207B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114742330B (zh) * 2022-06-13 2022-09-13 西南石油大学 一种高含硫有水气藏水封气量预测方法
CN115935121B (zh) * 2023-02-23 2023-05-12 中国石油大学(华东) 一种单相排水期最大日产水量确定方法及系统
CN117648523B (zh) * 2024-01-29 2024-04-05 成都英沃信科技有限公司 一种有水气藏动态储量及水侵常数计算方法

Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104453876A (zh) * 2014-11-03 2015-03-25 中国石油天然气股份有限公司 致密油气储层水平井油气产量的预测方法及预测装置
CN104832156A (zh) * 2015-05-05 2015-08-12 中国石油天然气股份有限公司 一种预估气井产量的方法
CN105134191A (zh) * 2015-08-25 2015-12-09 中国石油天然气股份有限公司 致密油油井储量的评价方法
CN105528648A (zh) * 2015-11-30 2016-04-27 中国石油天然气股份有限公司 缝洞单元的生产动态预测方法及装置
CN106194154A (zh) * 2016-07-11 2016-12-07 西南石油大学 一种非常规油气藏中长期产能预测方法
CN107423844A (zh) * 2017-06-06 2017-12-01 西南石油大学 一种预测页岩气/致密气井可采储量的新方法
CN107435528A (zh) * 2016-05-20 2017-12-05 中国石油天然气股份有限公司 火山岩气藏气井配产的方法
CN107944599A (zh) * 2017-10-31 2018-04-20 中国石油天然气股份有限公司 油气水平井产量的预测方法
CN108280534A (zh) * 2017-12-25 2018-07-13 中国石油天然气股份有限公司 一种气井产量递减率预测方法
CN108694254A (zh) * 2017-04-06 2018-10-23 中国石油化工股份有限公司 一种变产变压生产气井产量经验递减曲线分析方法
CN108804819A (zh) * 2018-06-10 2018-11-13 西南石油大学 一种低渗气藏动态储量评价方法
CN109033541A (zh) * 2018-07-02 2018-12-18 中国地质大学（北京） 一种基于eur的压后页岩气藏非均质性评价方法
CN110242291A (zh) * 2019-07-02 2019-09-17 中国石油化工股份有限公司 一种超强非均质油气储层非均质性表征方法
CN110334431A (zh) * 2019-07-02 2019-10-15 西南石油大学 一种低渗透致密气藏单井控制储量计算及剩余气分析方法
CN110552682A (zh) * 2019-09-16 2019-12-10 西南石油大学 一种碳酸盐岩三重介质气藏直井产量递减分析方法
CN111091293A (zh) * 2019-12-18 2020-05-01 中国石油天然气股份有限公司 一种油藏开发动态预警方法
CN111502652A (zh) * 2020-07-01 2020-08-07 西南石油大学 一种三孔介质气藏水平井产量递减及生产动态预测方法
CN112360422A (zh) * 2020-12-08 2021-02-12 西南石油大学 一种页岩气藏压裂水平井产量预测方法及系统

Patent Citations (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104453876A (zh) * 2014-11-03 2015-03-25 中国石油天然气股份有限公司 致密油气储层水平井油气产量的预测方法及预测装置
CN104832156A (zh) * 2015-05-05 2015-08-12 中国石油天然气股份有限公司 一种预估气井产量的方法
CN105134191A (zh) * 2015-08-25 2015-12-09 中国石油天然气股份有限公司 致密油油井储量的评价方法
CN105528648A (zh) * 2015-11-30 2016-04-27 中国石油天然气股份有限公司 缝洞单元的生产动态预测方法及装置
CN107435528A (zh) * 2016-05-20 2017-12-05 中国石油天然气股份有限公司 火山岩气藏气井配产的方法
CN106194154A (zh) * 2016-07-11 2016-12-07 西南石油大学 一种非常规油气藏中长期产能预测方法
CN108694254A (zh) * 2017-04-06 2018-10-23 中国石油化工股份有限公司 一种变产变压生产气井产量经验递减曲线分析方法
CN107423844A (zh) * 2017-06-06 2017-12-01 西南石油大学 一种预测页岩气/致密气井可采储量的新方法
CN107944599A (zh) * 2017-10-31 2018-04-20 中国石油天然气股份有限公司 油气水平井产量的预测方法
CN108280534A (zh) * 2017-12-25 2018-07-13 中国石油天然气股份有限公司 一种气井产量递减率预测方法
CN108804819A (zh) * 2018-06-10 2018-11-13 西南石油大学 一种低渗气藏动态储量评价方法
CN109033541A (zh) * 2018-07-02 2018-12-18 中国地质大学（北京） 一种基于eur的压后页岩气藏非均质性评价方法
CN110242291A (zh) * 2019-07-02 2019-09-17 中国石油化工股份有限公司 一种超强非均质油气储层非均质性表征方法
CN110334431A (zh) * 2019-07-02 2019-10-15 西南石油大学 一种低渗透致密气藏单井控制储量计算及剩余气分析方法
CN110552682A (zh) * 2019-09-16 2019-12-10 西南石油大学 一种碳酸盐岩三重介质气藏直井产量递减分析方法
CN111091293A (zh) * 2019-12-18 2020-05-01 中国石油天然气股份有限公司 一种油藏开发动态预警方法
CN111502652A (zh) * 2020-07-01 2020-08-07 西南石油大学 一种三孔介质气藏水平井产量递减及生产动态预测方法
CN112360422A (zh) * 2020-12-08 2021-02-12 西南石油大学 一种页岩气藏压裂水平井产量预测方法及系统

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
Rapid and accurate evaluation of reserves in different types of shale-gas wells: Production-decline analysis;Ke Wang 等;《International Journal of Coal Geology》;20200229;1-12页 *
三重介质复合气藏水平井不稳定产量递减动态分析;黄雨,李晓平,谭晓华;《天然气地球科学》;20181231;1190-1197页 *
低渗透压敏油藏产量递减动态预测方法研究;杨勇;《大庆石油地质与开发》;20081231;64-67页 *
计算天然气偏差因子的DAK方法的修正;胡建国，郭分乔，许进进;《石油与天然气地质》;20130228;120-123页 *
页岩气井压力递减分析新方法;彭朝阳,李井亮,韩永胜,王林,黄小青,吕一;《天然气勘探与开发》;20200625;104-109页 *
页岩气多级压裂水平井不稳定产量递减探讨;谢维扬,李晓平,张烈辉,王俊超,程子洋,袁淋;《天然气地球科学》;20151231;384-390页 *
高压产水气井井底压力计算方法研究;李玲、王波;《内蒙古石油化工》;20111231;224-225页 *

Also Published As

Publication number Publication date
US20230017069A1 (en) 2023-01-19
CN113236207A (zh) 2021-08-10

Similar Documents

Publication Publication Date Title
CN113236207B (zh) 一种强非均质性储层中产水气井的定产量递减预测方法
CN107563899A (zh) 油气井产能预测方法及装置
CN103592687B (zh) 煤岩吸附气含量的定量计算方法
CN111927444A (zh) 一种枯竭油气藏储气库注气能力评价方法
CN115587674B (zh) 油藏改建储气库扩容达产过程气井动态产能预测方法
CN107130955B (zh) 井底流压的确定方法及储集体天然能量的确定方法
Yadua et al. Performance of a gas-lifted oil production well at steady state
CN112943229A (zh) 一种储气库盖层突破压力连续预测方法
CN115345090A (zh) 一种欠饱和煤层气储层压降漏斗动态传播的计算方法
CN112746836B (zh) 基于层间干扰的油井各层产量计算方法
CN116341299A (zh) 一种低渗油藏co2驱气窜识别的方法
Jiang et al. Reserves estimation for geopressured gas reservoirs
CN111287741B (zh) 一种致密油藏体积压裂改造区渗透率的快速计算方法
CN111219180B (zh) 在海上低渗透油田中筛选适合注减氧空气开发试验区的方法
JOEL et al. a Novel Approach To Forecasting Production Rate of Dry Gas Wells
Mutalik et al. Production data analysis of horizontal CBM wells in Arkoma Basin
CN115492573B (zh) 一种柱塞气举井地层流入动态确定方法
Yao et al. Case study on diagnosis and identify the degree of bottom hole liquid accumulation in double-branch horizontal wells in PCOC
CN114547850B (zh) 一种基于多元回归的气井早期采收率计算方法
CN115749702B (zh) 一种页岩油天然气吞吐注气量优化设计方法
Iktissanov et al. Updating of limited bottomhole pressures for carbonate and sandstone reservoirs operated by Tatneft PJSC (Russian)
CN116822111A (zh) 一种高含硫气藏动态储量的计算方法
Levitina Calculation of Reservoir Characteristics by the Pressure Curve at the Inflow Using Reference Curves
CN117874988A (zh) 基于生产数据求油井单井控制储量及动态地层压力的方法
Djossou Feasibility Analysis of Ethane and Flue Gas Injection for Enhanced Oil Recovery with a Niger Delta Case Study

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant