[发明专利]一种集合卡尔曼滤波方法及装置在审
| 申请号： | 201710517113.2 | 申请日： | 2017-06-29 |
| 公开（公告）号： | CN109214013A | 公开（公告）日： | 2019-01-15 |
| 发明（设计）人： | 周单 | 申请（专利权）人： | 中国石油化工股份有限公司;中国石油化工股份有限公司石油物探技术研究院 |
| 主分类号： | G06F17/50 | 分类号： | G06F17/50;G06N7/00 |
| 代理公司： | 北京思创毕升专利事务所 11218 | 代理人： | 孙向民;廉莉莉 |
| 地址： | 100728 北*** | 国省代码： | 北京;11 |
| 权利要求书： | 查看更多 | 说明书： | 查看更多 |
| 摘要： | |||
| 搜索关键词： | 静态模型 集合 校正 预测状态向量 卡尔曼滤波 地震响应 动态参数 静态参数 初始状态向量 地震数据 非均质性 高斯模拟 观测数据 滤波结果 模型实验 状态向量 同化 抽取 吸收 | ||
1.一种集合卡尔曼滤波方法，其特征在于，所述方法包括以下步骤：
1)建立储层静态模型，并利用序贯高斯模拟方法产生多个储层静态模型集合；
2)基于所述多个储层静态模型集合得到储层静态参数Rs、初始时刻动态参数Rd、地震差异数据ΔSe组成初始状态向量Sk,j；
3)基于所述多个储层静态模型集合和所述初始状态向量Sk,j得到所述多个储层静态模型集合的各个实现对应的动态参数和地震响应差异；
4)基于储层静态参数以及多个储层静态模型集合的各个实现对应的动态参数和地震响应差异预测状态向量；
5)以预测状态向量同化吸收观测数据进行储层静态模型的校正，并得到校正状态向量；
6)抽取校正后的各个实现对应的储层静态参数，求取平均值作为滤波结果。
2.根据权利要求1所述的集合卡尔曼滤波方法，其特征在于，所述多个储层静态模型集合的各个实现对应的动态参数通过对所述多个储层静态模型集合进行数值模拟而得到。
3.根据权利要求1所述的集合卡尔曼滤波方法，其特征在于，所述多个储层静态模型集合的各个实现对应的地震响应差异通过对岩石物理模型计算岩石物理参数的变化而得到。
4.根据权利要求1所述的集合卡尔曼滤波方法，其特征在于，步骤5)中所使用的校正公式为：
式中，Kgk、为卡尔曼滤波因子。
5.根据权利要求4所述的集合卡尔曼滤波方法，其特征在于，卡尔曼滤波因子Kgk、通过以下公式得到：
式中，观测因子Hk＝[0|I]，其中I单位阵，I的维数为观测数据个数，Ne为集合中实现的数目。
6.一种集合卡尔曼滤波装置，其特征在于，所述装置包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现以下步骤：
1)建立储层静态模型，并利用序贯高斯模拟方法产生多个储层静态模型集合；
2)基于所述多个储层静态模型集合得到储层静态参数Rs、初始时刻动态参数Rd、地震差异数据ΔSe组成初始状态向量Sk,j；
3)基于所述多个储层静态模型集合和所述初始状态向量Sk,j得到所述多个储层静态模型集合的各个实现对应的动态参数和地震响应差异；
4)基于储层静态参数以及多个储层静态模型集合的各个实现对应的动态参数和地震响应差异预测状态向量；
5)以预测状态向量同化吸收观测数据进行储层静态模型的校正，并得到校正状态向量；
6)抽取校正后的各个实现对应的储层静态参数，求取平均值作为滤波结果。
7.根据权利要求6所述的集合卡尔曼滤波装置，其特征在于，所述多个储层静态模型集合的各个实现对应的动态参数通过对所述多个储层静态模型集合进行数值模拟而得到。
8.根据权利要求6所述的集合卡尔曼滤波装置，其特征在于，所述多个储层静态模型集合的各个实现对应的地震响应差异通过对岩石物理模型计算岩石物理参数的变化而得到。
9.根据权利要求6所述的集合卡尔曼滤波装置，其特征在于，步骤5)中所使用的校正公式为：
式中，Kgk、为卡尔曼滤波因子。
10.根据权利要求9所述的集合卡尔曼滤波装置，其特征在于，卡尔曼滤波因子Kgk、通过以下公式得到：
式中，观测因子Hk＝[0|I]，其中I单位阵，I的维数为观测数据个数，Ne为集合中实现的数目。
该专利技术资料仅供研究查看技术是否侵权等信息，商用须获得专利权人授权。该专利全部权利属于中国石油化工股份有限公司;中国石油化工股份有限公司石油物探技术研究院，未经中国石油化工股份有限公司;中国石油化工股份有限公司石油物探技术研究院许可，擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作，请联系【客服】
本文链接：http://www.vipzhuanli.com/pat/books/201710517113.2/1.html，转载请声明来源钻瓜专利网。
