[发明专利]一种页岩气TOC叠前地震反演预测方法有效
申请号: | 201610707298.9 | 申请日: | 2016-08-23 |
公开(公告)号: | CN106324674B | 公开(公告)日: | 2018-06-12 |
发明(设计)人: | 宗兆云;印兴耀;吴国忱;王泓鉴 | 申请(专利权)人: | 中国石油大学(华东) |
主分类号: | G01V1/30 | 分类号: | G01V1/30 |
代理公司: | 北京汇泽知识产权代理有限公司 11228 | 代理人: | 胡小龙 |
地址: | 266580 山东省*** | 国省代码: | 山东;37 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明涉及非常规油气地震勘探领域,具体公开了一种页岩气TOC叠前地震反演预测方法,该方法包括以下步骤:步骤一,建立页岩储层TOC反演目标函数;步骤二,基于弹性阻抗的页岩TOC叠前反演:根据对测井数据的统计分析而建立所述储层TOC的先验分布,通过蒙特卡洛仿真模拟技术分别对已建立的先验分布做随机抽样,最终获得储层TOC的随机样本空间分布,估计储层TOC后验概率的最大值,此最大值所在的位置处所对应的TOC值就是反演的最终结果。本发明综合应用了贝叶斯理论、统计岩石物理模型、蒙特卡洛随机抽样技术等理论,能够同时反演几种物性参数,消除了单独反演某一种参数时受到其他参数限制的影响,进而增强了反演的可信度。 | ||
搜索关键词: | 反演 随机抽样 先验分布 页岩气 页岩 蒙特卡洛仿真 岩石物理模型 油气地震勘探 参数限制 测井数据 弹性阻抗 后验概率 空间分布 模拟技术 目标函数 随机样本 物性参数 最终结果 地震 统计分析 贝叶斯 可信度 预测 应用 统计 | ||
【主权项】:
一种页岩气TOC叠前地震反演预测方法,其特征在于,其包括以下步骤:步骤一,建立页岩储层TOC反演目标函数:贝叶斯公式表示为: 式(1‑1)中,A代表观测样本信息,B则表示待估计参数,p(B)代表待估计参数B的先验分布,P(B|A)代表待估计参数B的后验分布,P(A|B)是联系随机事件A和随机事件B的似然函数,∫P(B)P(A|B)dB是观测样本的全概率,设∫P(B)P(A|B)dB是一个常数,记为将式(1‑1)应用于储层TOC反演过程中目标函数的推导得到:P([φ,Vsh,T]|[EI1,EI2,EI3])=α×P([φ,Vsh,T])P([EI1,EI2,EI3]|[φ,Vsh,T]) (1‑4),式(1‑4)中,φ,Vsh,T依次表示孔隙度、泥质含量以及TOC,EI1,EI2,EI3分别表示三个不同角度的弹性阻抗,P(·)代表概率密度函数,对三个不同角度的弹性阻抗做去相关处理,使三者彼此之间相互独立,对式(1‑4)做变形得到下式:P([EI1,EI2,EI3]|[φ,Vsh,T])=P(EI1|[φ,Vsh,T])P(EI2|[φ,Vsh,T])P(EI3|[φ,Vsh,T]) (1‑5),式(1‑5)中,EI1,EI2,EI3分别为经去相关处理后的弹性阻抗参数,寻找后验分布中最大后验概率的所在位置,该位置所对应的参数值即为储层TOC的最终反演结果:[φ,Vsh,T]=arg MaxP([φ,Vsh,T]|[EI1,EI2,EI3]) (1‑6),将式(1‑4)和式(1‑5)代入式(1‑6)中,舍弃常数α,建立最终反演的目标函数,即: 式(1‑7)中,P([φ,Vsh,T])为储层TOC的先验分布,P(EI1|[φ,Vsh,T])P(EI2|[φ,Vsh,T])P(EI3|[φ,Vsh,T])为联系先验分布和后验概率分布的似然函数;步骤二,基于弹性阻抗的页岩TOC叠前反演:根据对测井数据的统计分析而建立所述储层TOC的先验分布:假定测井资料中所提供的各储层物性参数均服从混合高斯分布,并且每个混合高斯分布都是由N个高斯分量通过加权平均构成,利用EM算法求出混合高斯分布中每一个高斯分量中的各项统计参数:均值、方差及权重,然后根据各项统计参数分别建立N个高斯分布,并按照它们各自的权重进行加权平均,最终得出储层TOC的混合高斯分布,其表达式如下:
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国石油大学(华东),未经中国石油大学(华东)许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201610707298.9/,转载请声明来源钻瓜专利网。