[发明专利]一种改进的油气藏高效自动历史拟合方法有效
申请号: | 202110770767.2 | 申请日: | 2021-07-08 |
公开(公告)号: | CN113486556B | 公开(公告)日: | 2022-03-11 |
发明(设计)人: | 张钊 | 申请(专利权)人: | 山东大学 |
主分类号: | G06F30/23 | 分类号: | G06F30/23;G06Q50/02;G06F119/14 |
代理公司: | 北京卓胜佰达知识产权代理有限公司 16026 | 代理人: | 张串串 |
地址: | 250100 *** | 国省代码: | 山东;37 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 改进 油气藏 高效 自动 历史 拟合 方法 | ||
1.一种改进的油气藏高效自动历史拟合方法,其特征在于,包括以下步骤:
S1、获取井点硬数据,并对所述井点硬数据进行处理,得到若干个先验地质模型,将若干个先验地质模型看作一个集合;
S2、基于所述先验地质模型使用快速行进法和差分进化算法进行拟合计算:
S2.1、获取目标函数值,所述目标函数值为预测数据和观测数据之间误差的最小值,所述目标函数值包括第一目标函数值、第二目标函数值,将第一目标函数值与第二目标函数值进行比较,若第一目标函数值小于第二目标函数值,则将第二目标函数值对应的先验地质模型替换为新地质模型,将所述新地质模型放到所述集合中,组成一个新的集合,从而完成一次进化;
所述第一目标函数值为新地质模型的目标函数值,所述第二目标函数值为所述集合中的目标函数最大值;
S2.2、重复所述S2.1直至所述新的集合中的所有地质模型对应的目标函数收敛;
所述新地质模型的获取过程如下:
A、对每个先验地质模型进行网格划分得到若干个网格单元,每个先验地质模型的网格单元一一对应,在集合中随机选择两个先验地质模型,计算两个先验地质模型每个对应网格单元之间差值,将所述差值加到随机选择第三个先验地质模型的对应网格单元上,从而得到一个所述新地质模型;
对新地质模型,再执行以下步骤:
计算网格单元的扩散飞行时间的传播速度f,
其中k为渗透率,为孔隙度,μ是动力粘度,ct是总压缩系数;
利用快速行进法求解程函方程得到网格单元上τ的值;
计算所述目标函数值y:
其中,Nt为时间步数,pi是第i个时间步的井点压力模拟值,是第i个时间步的井点压力观测值,也就是拟合目标。
2.如权利要求1所述的改进的油气藏高效自动历史拟合方法,其特征在于,S1中的所述处理过程为:使用序贯高斯模拟、直接序贯模拟、序贯高斯协模拟和直接序贯协模拟其中的一种,对所述井点硬数据进行处理,得到若干个先验地质模型,若干个先验地质模型均不相同。
3.如权利要求2所述的改进的油气藏高效自动历史拟合方法,其特征在于,
利用欧拉法求解:
得到井点处压力随时间的变化,其中p为井点处压力模拟值,qw为井点处流量,ct是总压缩系数,t为网格单元实际飞行时间,V(t)为t时刻的泄油体积。
4.如权利要求3所述的改进的油气藏高效自动历史拟合方法,其特征在于,所述t时刻的泄油体积V(t)通过累加小于等于t时刻的网格单元体积的得到。
5.如权利要求4所述的改进的油气藏高效自动历史拟合方法,其特征在于,基于网格单元的扩散飞行时间τ计算得到所述网格单元实际飞行时间t:
其中,β为常系数。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于山东大学,未经山东大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202110770767.2/1.html,转载请声明来源钻瓜专利网。