[发明专利]一种瞬时相位解缠绕方法、全波形反演方法及计算机设备有效
申请号: | 201910982006.6 | 申请日: | 2019-10-16 |
公开(公告)号: | CN110687594B | 公开(公告)日: | 2020-11-24 |
发明(设计)人: | 张凯;蔺玉曌;李振春;向晨;许鑫 | 申请(专利权)人: | 中国石油大学(华东) |
主分类号: | G01V1/28 | 分类号: | G01V1/28;G06F17/11 |
代理公司: | 北京汇泽知识产权代理有限公司 11228 | 代理人: | 刘淑敏 |
地址: | 266580 山东省青岛市西海岸*** | 国省代码: | 山东;37 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 瞬时 相位 缠绕 方法 波形 反演 计算机 设备 | ||
1.一种全波形反演方法,其特征在于,包括:
步骤1,基于初始的模型参数按照确定地震子波对应的正演地震信号其中,F为正演算子,是空间坐标,t是地震数据时间,i的初始值为0;
步骤2,确定正演地震信号的瞬时相位其中,其中,为的希尔伯特变换;
步骤3,按照对正演地震信号的瞬时相位进行解缠绕,得到解缠绕后的正演瞬时相位其中为的包络,C为积分路径,为相位梯度,为的初始相位;
步骤4,按照对地震子波对应的观测地震信号的瞬时相位进行解缠绕,得到解缠绕后的观测瞬时相位其中,为的包络,为相位梯度,为的初始相位;
步骤5,确定基于所述正演瞬时相位和所述观测瞬时相位确定的目标函数的梯度其中,其中,F-1为正演算子F的转置,为伴随震源,且
步骤6,按照对模型参数进行迭代,得到迭代后的模型参数并按照i=i+1更新迭代次数i,其中,α为步长;
步骤7,判断迭代次数i是否达到预设次数,如果否,将所述迭代的模型参数作为所述步骤1中的初始的模型参数返回所述步骤1;如果是,将迭代后的模型参数作为模型参数
2.根据权利要求1所述的全波形反演方法,其特征在于,所述模型参数为速度。
3.根据权利要求1所述的全波形反演方法,其特征在于,所述步长α通过抛物线插值方法确定。
4.一种瞬时相位解缠绕方法,其特征在于,包括:
确定正演地震信号的瞬时相位其中,是空间坐标,t是地震数据时间,为的希尔伯特变换;
按照对信号的瞬时相位进行解缠绕,得到解缠绕后的瞬时相位其中为的包络,C为积分路径,为相位梯度,为的初始相位。
5.一种计算机设备,其特征在于,包括:
存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序;
所述计算机程序被所述处理器执行时实现如权利要求1至3中任一项或4所述的方法的步骤。
6.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有全波形反演程序,所述全波形反演程序被处理器执行时实现如权利要求1至3中任一项所述的全波形反演方法的步骤。
7.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有瞬时相位解缠绕程序,所述瞬时相位解缠绕程序被处理器执行时实现如权利要求4所述的瞬时相位解缠绕方法的步骤。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国石油大学(华东),未经中国石油大学(华东)许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201910982006.6/1.html,转载请声明来源钻瓜专利网。