[发明专利]基于地震振幅差的反演方法及系统有效
申请号: | 201910960266.3 | 申请日: | 2019-10-10 |
公开(公告)号: | CN112649844B | 公开(公告)日: | 2023-08-15 |
发明(设计)人: | 周单;唐金良 | 申请(专利权)人: | 中国石油化工股份有限公司;中国石油化工股份有限公司石油物探技术研究院 |
主分类号: | G01V1/28 | 分类号: | G01V1/28 |
代理公司: | 北京思创毕升专利事务所 11218 | 代理人: | 孙向民;廉莉莉 |
地址: | 100027 北*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 基于 地震 振幅 反演 方法 系统 | ||
1.一种基于地震振幅差的反演方法,其特征在于,包括:
针对时移地震前后两次的地震数据分别进行归一化处理,获得归一化新地震数据与归一化旧地震数据;
计算所述归一化新地震数据与所述归一化旧地震数据的振幅差值;
根据所述振幅差值建立关于反射系数的目标函数,计算所述目标函数趋于0时的反射系数;
根据所述反射系数,针对每一道的振幅差值进行波阻抗反演,获得每一道的反演结果;
其中,所述目标函数为:
J=(W(t)R(t)-seis(t)*)T(W(t)R(t)-seis(t)*)+λ(R(t)TR(t)-C) (3)
其中,J为目标函数,R(t)为反射系数,W(t)为给定子波,seis(t)*为振幅差值,λ为阻尼系数,C为常数;
其中,通过公式(4)进行波阻抗反演:
其中,Zn+1为波阻抗反演,Z0为该道第一个振幅差值,N为每个地震道数值的总个数。
2.根据权利要求1所述的基于地震振幅差的反演方法,其中,所述归一化新地震数据为:
其中,为第t个归一化新地震数据,seis(t)new为第t个新地震数据,seis(min)new为新地震数据的最小值,seis(max)new为新地震数据的最大值。
3.根据权利要求1所述的基于地震振幅差的反演方法,其中,所述归一化旧地震数据为:
其中,为第t个归一化旧地震数据,seis(t)old为第t个旧地震数据,seis(min)old为旧地震数据的最小值,seis(max)old为旧地震数据的最大值。
4.一种基于地震振幅差的反演系统,其特征在于,该系统包括:
存储器,存储有计算机可执行指令;
处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:
针对时移地震前后两次的地震数据分别进行归一化处理,获得归一化新地震数据与归一化旧地震数据;
计算所述归一化新地震数据与所述归一化旧地震数据的振幅差值;
根据所述振幅差值建立关于反射系数的目标函数,计算所述目标函数趋于0时的反射系数;
根据所述反射系数,针对每一道的振幅差值进行波阻抗反演,获得每一道的反演结果;
其中,所述目标函数为:
J=(W(t)R(t)-seis(t)*)T(W(t)R(t)-seis(t)*)+λ(R(t)TR(t)-C) (3)
其中,J为目标函数,R(t)为反射系数,W(t)为给定子波,seis(t)*为振幅差值,λ为阻尼系数,C为常数;
其中,通过公式(4)进行波阻抗反演:
其中,Zn+1为波阻抗反演,Z0为该道第一个振幅差值,N为每个地震道数值的总个数。
5.根据权利要求4所述的基于地震振幅差的反演系统,其中,所述归一化新地震数据为:
其中,为第t个归一化新地震数据,seis(t)new为第t个新地震数据,seis(min)new为新地震数据的最小值,seis(max)new为新地震数据的最大值。
6.根据权利要求4所述的基于地震振幅差的反演系统,其中,所述归一化旧地震数据为:
其中,为第t个归一化旧地震数据,seis(t)old为第t个旧地震数据,seis(min)old为旧地震数据的最小值,seis(max)old为旧地震数据的最大值。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国石油化工股份有限公司;中国石油化工股份有限公司石油物探技术研究院,未经中国石油化工股份有限公司;中国石油化工股份有限公司石油物探技术研究院许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201910960266.3/1.html,转载请声明来源钻瓜专利网。
- 上一篇:一种人体结构化的方法
- 下一篇:一种仪表控制系统软件升级的方法、装置和汽车