[发明专利]多波动态匹配方法及系统有效
申请号: | 201810630263.9 | 申请日: | 2018-06-19 |
公开(公告)号: | CN110618448B | 公开(公告)日: | 2021-02-09 |
发明(设计)人: | 刘韬;陈天胜;刘春园;季玉新 | 申请(专利权)人: | 中国石油化工股份有限公司;中国石油化工股份有限公司石油勘探开发研究院 |
主分类号: | G01V1/28 | 分类号: | G01V1/28;G01V1/30 |
代理公司: | 北京思创毕升专利事务所 11218 | 代理人: | 孙向民;廉莉莉 |
地址: | 100027 北*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 波动 匹配 方法 系统 | ||
1.一种多波动态匹配方法,其特征在于,包括:
步骤1:输入初始纵横波速比,进而获得初始纵横波速比变化率;
步骤2:通过所述初始纵横波速比进行时间匹配,获得纵波和匹配后的转换波;
步骤3:分别对所述纵波和所述转换波进行AVO反演,根据反演属性值计算基于反演的纵横波速比变化率;
步骤4:将所述基于反演的纵横波速比变化率与所述初始纵横波速比变化率做差,获得纵横波速比变化率差值:
步骤5:根据所述纵横波速比变化率差值,调整所述初始纵横波速比,重复步骤1-4,直到所述纵横波速比变化率差值最小,将此时的转换波作为最终转换波;
其中,所述基于反演的纵横波速比变化率通过公式(3)计算:
其中,为基于反演的纵横波速比变化率,n为时间采样点,γ为纵横波速比,A、B为纵波属性参数,D为转换波属性参数;
其中,所述初始纵横波速比变化率根据公式(6)计算:
其中,n为时间采样点,γt为初始纵横波速比。
2.根据权利要求1所述的多波动态匹配方法,其中,所述反演属性值包括纵波属性参数和转换波属性参数,其中,
所述纵波属性参数的值通过对所述纵波进行AVO反演获得;
所述转换波属性参数的值通过对所述转换波进行AVO反演获得。
3.根据权利要求2所述的多波动态匹配方法,其中,通过公式(1)对所述纵波进行AVO反演,获得纵波属性参数的值:
Rpp(θ)≈A+Bsin2θ+C(tan2θ-sin2θ) (1)
其中,Rpp(θ)为纵波反射系数,θ为角度,A、B、C为纵波属性参数。
4.根据权利要求2所述的多波动态匹配方法,其中,通过公式(2)对所述纵波进行AVO反演,获得转换波属性参数的值:
Rps(θ)≈Dsinθ+Esin3θ (2)
其中,Rps(θ)为转换波反射系数,D、E为转换波属性参数,θ为角度。
5.根据权利要求1所述的多波动态匹配方法,其中,所述纵横波速比通过公式(4)计算:
6.根据权利要求1所述的多波动态匹配方法,其中,所述纵横波速比差值根据公式(5)计算:
其中,△F为纵横波速比变化率差值,为初始纵横波速比变化率。
7.一种多波动态匹配系统,其上存储有计算机程序,其特征在于,所述程序被处理器执行时实现以下步骤:
步骤1:输入初始纵横波速比,进而获得初始纵横波速比变化率;
步骤2:通过所述初始纵横波速比进行时间匹配,获得纵波和匹配后的转换波;
步骤3:分别对所述纵波和所述转换波进行AVO反演,根据反演属性值计算基于反演的纵横波速比变化率;
步骤4:将所述基于反演的纵横波速比变化率与所述初始纵横波速比变化率做差,获得纵横波速比变化率差值:
步骤5:根据所述纵横波速比变化率差值,调整所述初始纵横波速比,重复步骤1-4,直到所述纵横波速比变化率差值最小,将此时的转换波作为最终转换波;
其中,所述基于反演的纵横波速比变化率通过公式(3)计算:
其中,为基于反演的纵横波速比变化率,n为时间采样点,γ为纵横波速比,A、B为纵波属性参数,D为转换波属性参数;
其中,所述初始纵横波速比变化率根据公式(6)计算:
其中,n为时间采样点,γt为初始纵横波速比。
8.根据权利要求7所述的多波动态匹配系统,其中,所述反演属性值包括纵波属性参数和转换波属性参数,其中,
所述纵波属性参数的值通过对所述纵波进行AVO反演获得;
所述转换波属性参数的值通过对所述转换波进行AVO反演获得。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国石油化工股份有限公司;中国石油化工股份有限公司石油勘探开发研究院,未经中国石油化工股份有限公司;中国石油化工股份有限公司石油勘探开发研究院许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201810630263.9/1.html,转载请声明来源钻瓜专利网。