[发明专利]一种基于希尔伯特变换的振幅增益方法及处理终端有效
申请号: | 201810574209.7 | 申请日: | 2018-06-06 |
公开(公告)号: | CN109031414B | 公开(公告)日: | 2020-03-27 |
发明(设计)人: | 刘玉萍;张宝金 | 申请(专利权)人: | 广州海洋地质调查局 |
主分类号: | G01V1/30 | 分类号: | G01V1/30 |
代理公司: | 广州君咨知识产权代理有限公司 44437 | 代理人: | 江超 |
地址: | 510000 广东*** | 国省代码: | 广东;44 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 基于 希尔伯特 变换 振幅 增益 方法 处理 终端 | ||
1.一种基于希尔伯特变换的振幅增益方法,其特征在于:包括以下步骤:
步骤S1:获取原始地震资料,并对地震资料进行预处理,获得处理后的叠前道集;
步骤S2:提取出叠前道集的包络振幅和相位值,根据公式①获得地震波的包络振幅:
式中,A(t)为地震波的包络振幅,x(t)为实际地震记录道,从叠前道集中直接获得,为x(t)的正交道,由希尔伯特变换求得,
根据公式②或③获得地震波的相位值:
x(t)=A(t)cosθ(t)------②
其中,θ(t)为地震波的相位值;
步骤S3:分离出包络振幅的主分量和奇异值分量,首先,求得步骤S2中的每道包络振幅A(t)的平均振幅Aave值,然后,根据平均振幅Aave值将步骤S2中获得的A(t)值分离出主分量B(t)值:
包络振幅A(t)中小于Aave的值保持不变,大于Aave的值则用Aave代替,从而获得主分量B(t)值,根据公式④计算出包络振幅的奇异值分量C(t):
C(t)=A(t)-B(t)------④
步骤S4:计算权重w值,权重w值由公式⑤获得:
步骤S5:构建新的包络振幅Anew,其由公式⑥获得:
Anew(t)=B(t)+w*C(t)------⑥
最终,根据新构建的包络振幅Anew来获得重新构建后的地震道集,用xnew表示,并由公式⑦获得:
xnew(t)=Anew(t)cosθ(t)------⑦
其中,θ(t)值由步骤S2中计算获得。
2.根据权利要求1所述的基于希尔伯特变换的振幅增益方法,其特征在于:所述预处理包括球面扩散补偿、去噪音、涌浪静校正、去气泡、多次波压制和信号增强中的一种或多种。
3.一种处理终端,其特征在于:包括,
存储器,用于存储程序指令;
处理器,用于运行所述程序指令,以执行如下步骤:
步骤S1:获取原始地震资料,并对地震资料进行预处理,获得处理后的叠前道集;
步骤S2:提取出叠前道集的包络振幅和相位值,根据公式①获得地震波的包络振幅:
式中,A(t)为地震波的包络振幅,x(t)为实际地震记录道,可从叠前道集中直接获得,为x(t)的正交道,由希尔伯特变换求得,
根据公式②或③获得地震波的相位值:
x(t)=A(t)cosθ(t)------②
其中,θ(t)为地震波的相位值;
步骤S3:分离出包络振幅的主分量和奇异值分量,首先,求得步骤S2中的每道包络振幅A(t)的平均振幅Aave值,然后,根据平均振幅Aave值将步骤S2中获得的A(t)值分离出主分量B(t)值:
包络振幅A(t)中小于Aave的值保持不变,大于Aave的值则用Aave代替,从而获得主分量B(t)值,根据公式④计算出包络振幅的奇异值分量C(t):
C(t)=A(t)-B(t)------④
步骤S4:计算权重w值,权重w值由公式⑤获得:
步骤S5:构建新的包络振幅Anew,其由公式⑥获得:
Anew(t)=B(t)+w*C(t)------⑥
最终,根据新构建的包络振幅Anew来获得重新构建后的地震道集,用xnew表示,并由公式⑦获得:
xnew(t)=Anew(t)cosθ(t)------⑦
其中,θ(t)值由步骤S2中计算获得。
4.根据权利要求3所述的处理终端,其特征在于:所述预处理包括球面扩散补偿、去噪音、涌浪静校正、去气泡、多次波压制和信号增强中的一种或多种。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于广州海洋地质调查局,未经广州海洋地质调查局许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201810574209.7/1.html,转载请声明来源钻瓜专利网。