专栏首页DSPSTACK.COM地震信号系列完结篇-反卷积方法
原创

地震信号系列完结篇-反卷积方法

前言

本篇将详细地讲解地震信号中用到的反卷积方法。反卷积方法的作用在文章 地震信号的一些基本概念 中已经阐述过,简单的说就是:在压缩原信号的同时,对频谱进行补偿(反卷积的输出信号)。而在地震信号处理中,除了前面的优势,该方法还可以使反射系数成为尖脉冲,从而提高地震记录的分辨率。

地震记录:就是由仪器在地表所测的地震数据 反射系数:表征的是地球的每一层对地震波的反射特征,相当于一个系统或是滤波器,可看作为滤波器系数 地震记录的分辨率:地震记录的目的是想把地球的每一层进行区分和识别,而区分的细度就是所谓的分辨率,所能识别的底层越多,则分辨率越高,反之越低。

地震记录定义(输出信号)

我们首先给出地震记录的定义:

其中

:地震记录,
:地震子波,
:反射系数序列

地震子波包括震源子波,记录滤波器,地表反射和检波器的组合响应。上述定义省去了噪音的影响。

地震子波:与相关子波(相关序列)是不同的概念,它不是相关序列

反卷积要处理的问题

在地震勘探中,当地下层位较密集时,它们之间的反射层将难以区分,见下图:

前三个反射层间距很接近。

所以希望反射层的反射系数为一系列的尖脉冲,这样就能从地震记录中把各反射层区分开来,从而提高了地震记录的分辨率。如何将非尖脉冲的信号压缩成尖脉冲的反射系数,就是反卷积要处理的问题和重点。

反卷积方法

首先我们使用地震记录定义

对其进行时域到频域的变换(傅里叶变换),得到

其中

为地震记录
的频谱;

显然,

转到时域为:

其中

根据上式和卷积的同元性质可知:

其中

地震反子波。

根据已知地震子波

,即可求出地震反子波
,而根据卷积模型

就可求出反射系数
上述过程是反卷积的全过程,特别地,它是反卷积的一个特例,即反滤波,所谓反滤波是因为滤波算子
是地震子波的逆。

地震子波最小相位的假设:由于在求取反卷积滤波算子时,地震子波的频谱作为分母出现,很明显为了使反滤波算子稳定,地震子波必须是最小相位的,因为最小相位的零点都位于单位圆内,也就是作为分母出现的极点

扩展知识

上一节的过程看似顺理成章,其实隐含了一个最优方法求解滤波算子的过程。我们从最初的地震记录卷积模型得到:

时,即为上一节的式(1)。

其实

是一个期望输出,而实际输出为
,那么使用最小二乘原理的目标误差函数为

所以最理想的情况就是误差函数为0,即

此外,为了反射系数达到更好的波形,我们可以选择俞式子波作为期望输出(该知识也超出本篇范围)。

小结

本篇既是地震信号补偿频谱的基础也是对提高反射坡面识别率的重要基础,地震信号的处理技术很多都基于反卷积知识,当然还有最小相位,混合相位和最大相位的区分,还有期望输出子波的选择等等,它们都是以此为基础来对地震信号进行处理的。

(如有错误,请指正)

原创声明,本文系作者授权云+社区发表,未经许可,不得转载。

如有侵权,请联系 yunjia_community@tencent.com 删除。

我来说两句

0 条评论
登录 后参与评论

相关文章

  • Savitsky-Golay 平滑算法

    通过多项式对移动窗口内的数据进行多项式最小二乘拟合,算出窗口内中心点关于其周围点的加权平均和。

    DSPSTACK.COM
  • 地震信号-相关子波零相位化

    子波的振幅不变,改变相位谱的滤波器,以达到子波形状改变的过程称为子波整形或整形滤波

    DSPSTACK.COM
  • 地震信号的一些基本概念

    可控震源是指通过一个与大地紧密耦合的振动平板,以反作用方式向地下传送一组连续振动的弹性波信号(又称扫描信号),再经过对地面接收到的反射波信号的处理和辨识,用于解...

    DSPSTACK.COM
  • Outlook自动转发

    林万程
  • 使用SAP云平台portal service之前,需要做好哪些准备

    Jerry Wang
  • mac软件这样使用,你知道几个。(技术创作101训练营)

    听了腾讯云社区的授课,加上工作之后,将会使用到windows笔记本作为办公电脑。在此就想着把平时在mac上使用的软件整理一下,能帮到自己的软件分享给大家。

    崔笑颜
  • Spark-ALS 分布式实现详解

    本文重点阐述了 Spark ML 库中 ALS 的实现,要看懂以上计算流程,请务必结合源代码理解,在实际源码实现中,使用了很多优化技巧。

    涂小刚
  • SAP SRM ABAP Webdynpro和CFCA usb key集成的一个原型开发

    Jerry Wang
  • OBS推流

    治电小白菜
  • 区块链浅谈

    本文是笔者作为小白,参阅了区块链相关书籍资料后了,整理的随手笔记,在项目团队做了分享,比较初级简洁,希望对第一次接触区块链的小伙伴有所帮助。

    MT

扫码关注云+社区

领取腾讯云代金券