[发明专利]叠前混合非线性反演方法及计算机可读存储介质有效
| 申请号： | 201710464861.9 | 申请日： | 2017-06-19 |
| 公开（公告）号： | CN109143346B | 公开（公告）日： | 2020-10-20 |
| 发明（设计）人： | 张远银;刘喜武;刘宇巍;霍志周;刘志远;刘炯;钱恪然 | 申请（专利权）人： | 中国石油化工股份有限公司;中国石油化工股份有限公司石油勘探开发研究院 |
| 主分类号： | G01V1/30 | 分类号： | G01V1/30 |
| 代理公司： | 北京思创毕升专利事务所 11218 | 代理人： | 孙向民;廉莉莉 |
| 地址： | 100728 北*** | 国省代码： | 北京;11 |
| 权利要求书： | 查看更多 | 说明书： | 查看更多 |
| 摘要： | |||
| 搜索关键词： | 混合 非线性 反演 方法 计算机 可读 存储 介质 | ||
1.一种叠前混合非线性反演方法，包括：
对叠前地震数据进行叠前线性反演，获得横波反射系数Rs的初始值与横纵波速度比Rk的初始值；
基于横波反射系数Rs、横纵波速度比Rk、密度梯度RD，构造组合弹性参量KS；
将叠前纵波反射系数公式表示为关于所述密度梯度RD、纵波反射系数RP和所述组合弹性参量KS的线性方程，对所述叠前纵波反射系数公式进行线性求解，获得所述密度梯度RD、纵波反射系数RP和所述组合弹性参量KS；
构造关于所述横波反射系数Rs和横纵波速度比Rk的目标函数，基于所述横波反射系数Rs的初始值与横纵波速度比Rk的初始值对所述目标函数进行非线性求解，进而获得最终的横波反射系数与最终的横纵波速度比；
其中，所述组合弹性参量KS为：
其中，所述目标函数为公式(4)或公式(5)：
其中，F(Rs,Rk)表示目标函数。
2.根据权利要求1所述的叠前混合非线性反演方法，其中，所述叠前纵波反射系数公式为：
其中，Rpp(θ)表示叠前纵波反射系数，θ表示反射角与入射角的平均值。
3.根据权利要求1所述的叠前混合非线性反演方法，其中，通过标准粒子群算法对所述目标函数进行非线性求解。
4.一种计算机可读存储介质，其上存储有计算机程序，其中，所述程序被处理器执行时实现以下步骤：
对叠前地震数据进行叠前线性反演，获得横波反射系数Rs的初始值与横纵波速度比Rk的初始值；
基于横波反射系数Rs、横纵波速度比Rk、密度梯度RD，构造组合弹性参量KS；
将叠前纵波反射系数公式表示为关于所述密度梯度RD、纵波反射系数RP和所述组合弹性参量KS的线性方程，对所述叠前纵波反射系数公式进行线性求解，获得所述密度梯度RD、纵波反射系数RP和所述组合弹性参量KS；
构造关于所述横波反射系数Rs和横纵波速度比Rk的目标函数，基于所述横波反射系数Rs的初始值与横纵波速度比Rk的初始值对所述目标函数进行非线性求解，进而获得最终的横波反射系数与最终的横纵波速度比；
其中，所述组合弹性参量KS为：
其中，所述目标函数为公式(4)或公式(5)：
其中，F(Rs,Rk)表示目标函数。
5.根据权利要求4所述的计算机可读存储介质，其中，所述叠前纵波反射系数公式为：
其中，Rpp(θ)表示叠前纵波反射系数，θ表示反射角与入射角的平均值。
该专利技术资料仅供研究查看技术是否侵权等信息，商用须获得专利权人授权。该专利全部权利属于中国石油化工股份有限公司;中国石油化工股份有限公司石油勘探开发研究院，未经中国石油化工股份有限公司;中国石油化工股份有限公司石油勘探开发研究院许可，擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作，请联系【客服】
本文链接：http://www.vipzhuanli.com/pat/books/201710464861.9/1.html，转载请声明来源钻瓜专利网。
