申请/专利权人：成都英沃信科技有限公司

申请日：2024-01-29

公开（公告）日：2024-04-05

公开（公告）号：CN117648523B

主分类号：G06F17/11

分类号：G06F17/11;G06Q50/02

优先权：

专利状态码：有效-授权

法律状态：2024.04.05#授权;2024.03.22#实质审查的生效;2024.03.05#公开

摘要：本发明涉及一种有水气藏动态储量及水侵常数计算方法，其技术方案是考虑已有常规水驱特征曲线高含水期与特高含水期拟合不准确，推导新型水驱特征曲线方程，给定动态储量G和水侵常数B的上下界，通过调整动态储量G和水侵常数B来拟合实际的水气比，得到动态储量G和水侵常数B，并绘制水气比拟合曲线图和气水相对渗透率拟合曲线图。本发明特点在于能反映高含水期水驱特征曲线上翘的现象，拟合更加准确，可用于计算动态储量以及水侵常数。

主权项：1.一种有水气藏动态储量及水侵常数计算方法，其特征在于，该方法包括以下步骤：S100、推导新型水驱特征曲线方程，步骤如下：S101、将存水体积系数的定义方程和方程依次代入存水体积与天然气地下体积之比方程推导含水饱和度方程： 式中，B为水侵常数，定义为，单位为无量纲；Sw为含水饱和度，单位为无量纲；Swi为原始条件下含水饱和度，单位为无量纲；R为采出程度，单位为无量纲；Bgi为原始条件下天然气体积系数，单位为无量纲；We为水侵量，单位为方；Wp为累产水，单位为方；为存水体积系数，单位为无量纲；Bw是地层水体积系数，单位为无量纲；G是动态储量，单位为万方；RB为采出程度的水侵常数次方，根据定义有采出程度的水侵常数次方RB等于存水体积系数；S102、将S101中推导出的含水饱和度方程和气水相对渗透率比值方程依次代入水气比方程推导新型水驱特征曲线方程： 其中，Bg为天然气的体积系数，定义为，单位为无量纲；psc为地面压力，单位为MPa；p为地层压力，单位为MPa；T为地层温度，单位为K；Tsc为地面温度，单位为K；Z为地层条件下的气体偏差因子，单位为无量纲；Zsc为地面条件下的气体偏差因子，单位为无量纲；Sw为含水饱和度，单位为无量纲；Swi为原始条件下含水饱和度，单位为无量纲；R为采出程度，单位为无量纲；WGR是水气比，单位为方万方；是天然气黏度，单位为mPa·s；是地层水黏度，单位为mPa·s；Bw是地层水体积系数，单位为无量纲；Krw是岩心水相相对渗透率，单位为无量纲；Krg是岩心气相相对渗透率，单位为无量纲；n、b为与储层和流体物性有关的常数，由气水相对渗透率曲线拟合得到，单位为无量纲；RB为采出程度的水侵常数次方，根据定义有采出程度的水侵常数次方RB等于存水体积系数；S200、准备静态资料和生产历史数据，静态资料包括气水相对渗透率曲线、原始含水饱和度、天然气黏度、地层水黏度、地层温度、地面温度、地面压力；生产历史数据包括日产气量、日产水量、地层压力；S300、根据S102中的新型水驱特征曲线方程，进行生产历史拟合的思路为：S301、给定动态储量G，水侵常数B两个参数的上下界限，其中动态储量G的上界为20000万方，动态储量G的下界为1000万方，水侵常数B的上界为20，水侵常数B的下界为1，上下界限内随机给定初值；采出程度R由累产Gp除以动态储量G计算，其中的动态储量G和水侵常数B均为可变参数；S302、结合S301中的动态储量G和水侵常数B代入新型水驱特征曲线方程可以算出预测的水气比；其中，在计算地层条件下的天然气体积系数Bg时，地层条件下的气体偏差因子Z由地层压力根据DAK法计算；S303、再重复S301~S302的过程，以生产数据中的一天为一个周期，计算得到整个生产阶段的水气比；S304、在参数的上下限内调整参数，采用自动拟合方法拟合实际的水气比，寻求预测值与实际值的最佳拟合，水气比拟合的收敛条件表示为： 其中，E为偏差值；m为迭代次数；lnWGRiG,B为预测水气比取对数的值，单位为方万方；lnWGRi为实际的水气比取对数的值，单位为方万方；S305、对有水气藏生产历史进行拟合，拟合成功后，停止调整动态储量G、水侵常数B两个参数，以拟合最优的动态储量G、水侵常数B为计算结果。

全文数据：

权利要求：

百度查询： 成都英沃信科技有限公司 一种有水气藏动态储量及水侵常数计算方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD