[发明专利]一种页岩气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反演目标函数:贝叶斯公式表示为:<mrow><mi>P</mi><mrow><mo>(</mo><mi>B</mi><mo>|</mo><mi>A</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>p</mi><mrow><mo>(</mo><mi>B</mi><mo>)</mo></mrow><mi>P</mi><mrow><mo>(</mo><mi>A</mi><mo>|</mo><mi>B</mi><mo>)</mo></mrow></mrow><mrow><mo>&Integral;</mo><mi>P</mi><mrow><mo>(</mo><mi>B</mi><mo>)</mo></mrow><mi>P</mi><mrow><mo>(</mo><mi>A</mi><mo>|</mo><mi>B</mi><mo>)</mo></mrow><mi>d</mi><mi>B</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo></mrow>式(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)中,舍弃常数α,建立最终反演的目标函数,即:<mrow><mtable><mtr><mtd><mrow><mo>&lsqb;</mo><mi>&phi;</mi><mo>,</mo><msub><mi>V</mi><mrow><mi>s</mi><mi>h</mi></mrow></msub><mo>,</mo><mi>T</mi><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mi>arg</mi><mi> </mi><mi>M</mi><mi>a</mi><mi>x</mi><mi>P</mi><mfenced open = "(" close = ")"><mtable><mtr><mtd><mrow><mi>P</mi><mrow><mo>(</mo><mo>&lsqb;</mo><mi>&phi;</mi><mo>,</mo><msub><mi>V</mi><mrow><mi>s</mi><mi>h</mi></mrow></msub><mo>,</mo><mi>T</mi><mo>&rsqb;</mo><mo>)</mo></mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>EI</mi><mn>1</mn></msub><mo>|</mo><mo>&lsqb;</mo><mi>&phi;</mi><mo>,</mo><msub><mi>V</mi><mrow><mi>s</mi><mi>h</mi></mrow></msub><mo>,</mo><mi>T</mi><mo>&rsqb;</mo><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>EI</mi><mn>2</mn></msub><mo>|</mo><mo>&lsqb;</mo><mi>&phi;</mi><mo>,</mo><msub><mi>V</mi><mrow><mi>s</mi><mi>h</mi></mrow></msub><mo>,</mo><mi>T</mi><mo>&rsqb;</mo><mo>)</mo></mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>EI</mi><mn>3</mn></msub><mo>|</mo><mo>&lsqb;</mo><mi>&phi;</mi><mo>,</mo><msub><mi>V</mi><mrow><mi>s</mi><mi>h</mi></mrow></msub><mo>,</mo><mi>T</mi><mo>&rsqb;</mo><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>7</mn><mo>)</mo></mrow><mo>,</mo></mrow>式(1‑7)中,P([φ,Vsh,T])为储层TOC的先验分布,P(EI1|[φ,Vsh,T])P(EI2|[φ,Vsh,T])P(EI3|[φ,Vsh,T])为联系先验分布和后验概率分布的似然函数;步骤二,基于弹性阻抗的页岩TOC叠前反演:根据对测井数据的统计分析而建立所述储层TOC的先验分布:假定测井资料中所提供的各储层物性参数均服从混合高斯分布,并且每个混合高斯分布都是由N个高斯分量通过加权平均构成,利用EM算法求出混合高斯分布中每一个高斯分量中的各项统计参数:均值、方差及权重,然后根据各项统计参数分别建立N个高斯分布,并按照它们各自的权重进行加权平均,最终得出储层TOC的混合高斯分布,其表达式如下:<mrow><mi>P</mi><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>&alpha;</mi><mi>k</mi></msub><mi>N</mi><mrow><mo>(</mo><mi>T</mi><mo>;</mo><msubsup><mi>&mu;</mi><mi>T</mi><mi>k</mi></msubsup><mo>,</mo><msubsup><mi>&Sigma;</mi><mi>T</mi><mi>k</mi></msubsup><mo>)</mo></mrow><mo</div> </div> <div class="b20"></div> <div class="down-box" id="down-box"> <div class="msg" style="display: block;"> <span>下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。</span> </div> <div class="btns"> <span class="btn paydown">免登录下载</span><a href="/login.html?p=8684656C99F5006F2432FC0618ACFC4928C45BC7BD2736E9" class="btn green" target="_blank">普通用户下载</a><a href="http://yh.vipzhuanli.com/member/service/pay-vip.html?p=v1" target="_blank" class="btn red">升级VIP会员,免费下载</a> </div> </div> <div class="warning"> <p>该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国石油大学(华东),未经中国石油大学(华东)许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【<a href="https://wpa1.qq.com/l11yQAzu?_type=wpa&qidian=true">客服</a>】</p> <p>本文链接:http://www.vipzhuanli.com/patent/201610707298.9/,转载请声明来源钻瓜专利网。</p> </div> <ul class="clear_div other_o"><li class="prev">上一篇:<a href="/patent/201610718098.3/" title="一种基于偏振干涉仪的光谱复原修正方法">一种基于偏振干涉仪的光谱复原修正方法</a></li><li class="next">下一篇:<a href="/patent/201610727565.9/" title="一种汽车内外饰件超市化设计平台及方法">一种汽车内外饰件超市化设计平台及方法</a></li></ul> <div class="oth-box"> <dl class="d_th"><dd><span>同类专利</span></dd><dt class="th_a"></dt></dl> <dl class="d_th" style="padding-top:15px;"><dd><span>专利分类</span></dd></dl> <div class="ps_c"> <div><a href="/ipc/G/" target="_blank" title="物理">G 物理</a></div><a class="ml1" href="/ipc/G01/" target="_blank" title="测量;测试">G01 测量;测试</a><br/><a class="ml2" href="/ipc/G01V/" target="_blank" title="地球物理;重力测量;物质或物体的探测;示踪物">G01V 地球物理;重力测量;物质或物体的探测;示踪物</a><br/><a class="ml3" href="/pat/ipc/G01V1/00/" target="_blank" title="地震学;地震或声学的勘探或探测">G01V1-00 地震学;地震或声学的勘探或探测</a><br/><a class="ml3" href="/pat/ipc/G01V1/02/" target="_blank" title=".产生地震能">G01V1-02 .产生地震能</a><br/><a class="ml3" href="/pat/ipc/G01V1/16/" target="_blank" title=".地震信号的接收元件;接收元件的配置或改进">G01V1-16 .地震信号的接收元件;接收元件的配置或改进</a><br/><a class="ml3" href="/pat/ipc/G01V1/22/" target="_blank" title=".将地震信号送至记录或处理设备">G01V1-22 .将地震信号送至记录或处理设备</a><br/><a class="ml3" href="/pat/ipc/G01V1/24/" target="_blank" title=".地震数据的记录">G01V1-24 .地震数据的记录</a><br/><a class="ml3" href="/pat/ipc/G01V1/28/" target="_blank" title=".地震数据的处理,例如,分析、用于解释、用于校正">G01V1-28 .地震数据的处理,例如,分析、用于解释、用于校正</a><br/> </div> </div> </div> <div class="content-r"> <div class="btns content-list" id="downdd"> <div class="header"> <div class="header-title"><a >专利文件下载</a></div> <hr /> </div> <span class="btn paydown">免登录下载</span><a href="/login.html?p=8684656C99F5006F2432FC0618ACFC4928C45BC7BD2736E9" class="btn green" target="_blank">普通用户下载</a><a href="http://yh.vipzhuanli.com/member/service/pay-vip.html?p=v1" target="_blank" class="btn red">升级VIP会员,免费下载</a> </div> <div class="content-list"> <div class="header"> <div class="header-title"><a href="/patent/list.html?kw=%e5%8f%8d%e6%bc%94 ">反演 相关专利</a></div> <hr /> </div> <ul> <li><a href="/patent/201510276042.2/">一种瞬变电磁电阻率的反演方法与装置</a></li> <li><a href="/patent/201610429648.X/">截断时窗的低通滤波多尺度全波形反演方法</a></li> <li><a href="/patent/201610694535.2/">电法与地震同步联合反演方法及系统</a></li> <li><a href="/patent/201710525723.7/">一种多模型协同的水质参数浓度反演方法和装置</a></li> <li><a href="/patent/201710527454.8/">一种地震数据反演方法及系统</a></li> <li><a href="/patent/201811140059.5/">一种分频迭代地震反演方法及系统</a></li> <li><a href="/patent/201910272992.6/">一种GNSS-R综合海洋环境要素反演装置以及方法</a></li> <li><a href="/patent/202010037302.1/">一种基于大地电磁和直流电阻率数据的联合反演方法</a></li> <li><a href="/patent/202010774885.6/">层状介质电阻率各向异性海洋可控源电磁快速反演方法</a></li> <li><a href="/patent/202011245522.X/">基于乘积函数的海洋可控源电磁与大地电磁联合反演方法</a></li> </ul> </div> <div class="content-list"> <div class="header"> <div class="header-title"><a href="/patent/list.html?kw=%e9%9a%8f%e6%9c%ba%e6%8a%bd%e6%a0%b7 ">随机抽样 相关专利</a></div> <hr /> </div> <ul> <li><a href="/patent/201220344707.0/">一种用于随机抽样的骰子及随机抽样装置</a></li> <li><a href="/patent/201410416213.2/">一种计及径流预报不确定性的水库调度风险评估方法</a></li> <li><a href="/patent/201510698056.3/">一种基于互补随机抽样的无透镜衍射成像方法</a></li> <li><a href="/patent/201610623109.X/">一种马赛克瓷砖随机抽样检测装置</a></li> <li><a href="/patent/201611203570.6/">一种基于改进的水塘抽样的大数据集多遍随机抽样方法</a></li> <li><a href="/patent/201620826898.2/">一种马赛克瓷砖随机抽样检测装置</a></li> <li><a href="/patent/201621289034.8/">一种化妆品合格检测用随机抽样设备</a></li> <li><a href="/patent/201710545574.0/">一种基于空间邻接指数的遥感分类野外调查样方抽选方法</a></li> <li><a href="/patent/202011119733.9/">一种发动机寿命预测方法、存储介质及计算设备</a></li> <li><a href="/patent/89215308.3/">随机抽样器</a></li> </ul> </div> <div class="content-list"> <div class="header"> <div class="header-title"><a href="/patent/list.html?kw=%e5%85%88%e9%aa%8c%e5%88%86%e5%b8%83 ">先验分布 相关专利</a></div> <hr /> </div> <ul> <li><a href="/patent/200510011550.4/">基于贝叶斯方法的软件运行剖面获取方法</a></li> <li><a href="/patent/201110257955.1/">一种基于混合贝叶斯先验分布的可靠性验证测试方法</a></li> <li><a href="/patent/201210125683.4/">一种贝叶斯软件可靠性验证测试方法及其计算机辅助工具</a></li> <li><a href="/patent/201380045315.1/">用于根据通过随机反演生成的后验分布确定信息价值度量的系统和方法</a></li> <li><a href="/patent/201610029926.2/">基于可靠度先验信息融合的冷备系统可靠度估计方法</a></li> <li><a href="/patent/201710110476.4/">一种商品评价得分的处理方法及系统</a></li> <li><a href="/patent/201810770361.2/">视频预测方法和装置、电子设备及车辆</a></li> <li><a href="/patent/201810812432.0/">视频预测方法、视频预测装置、电子设备和车辆</a></li> <li><a href="/patent/201911029097.8/">一种器官图像分割方法和装置</a></li> <li><a href="/patent/202011116184.X/">一种X射线图像线性重建方法</a></li> </ul> </div> <div class="content-list"> <div class="header"> <div class="header-title"><a href="/patent/list.html?kw=%e9%a1%b5%e5%b2%a9%e6%b0%94 ">页岩气 相关专利</a></div> <hr /> </div> <ul> <li><a href="/patent/201210227088.1/">页岩气增产方法及页岩气增产设备</a></li> <li><a href="/patent/201410015264.4/">页岩生烃、吸附及解吸模拟系统和模拟方法</a></li> <li><a href="/patent/201410597334.1/">一种页岩气扩散能力检测方法、装置及系统</a></li> <li><a href="/patent/201420019460.4/">页岩生烃、吸附及解吸模拟系统</a></li> <li><a href="/patent/201710811197.0/">一种页岩气富集成藏分类方法</a></li> <li><a href="/patent/201711384255.2/">获取页岩气的类型的方法和装置</a></li> <li><a href="/patent/201810321479.7/">一种页岩含气量的获取方法</a></li> <li><a href="/patent/201810414945.6/">页岩散失气总含量的确定方法、装置及存储介质</a></li> <li><a href="/patent/202010006082.6/">一种页岩气目标保存条件定量评价方法</a></li> <li><a href="/patent/202011402428.0/">一种页岩气资源评价方法及相关装置</a></li> </ul> </div> </div> </div> </div> <input type="hidden" id="hid_id" /> <script type="text/javascript"> /* <![CDATA[ */ var pat_ajax_url = "/down/check.html"; var wppay_ajax_url = "/pay/down"; var pnum = "201610707298.9"; var openNo = "CN106324674B"; var op = "20180612"; var y = "2018"; /* */
×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

tel code back_top