[发明专利]一种地震噪声去除方法及装置

摘要

本发明属于地震资料处理领域，具体属于一种地震噪声去除方法及装置。其中通过针对层状介质层间多次波预测算子的快速计算，得到层间多次波预测的炮集数据。基于预测算子的层间多次波预测技术可以一次预测所有层相关的层间多次波，利用层状介质的层间多次波和下覆层介质在层状介质中产生的层间多次波计算得到地震噪声，并在炮集数据中减去所计算的层间多次波炮集数据从而得到纯净的炮集数据。本发明的有益效果在于能够大幅度的减小计算代价。

主权项

一种地震噪声去除方法，其特征在于包括，步骤101，对炮集数据进行分离，分为层状介质数据Slayer和下覆层数据Syu；步骤102，对炮集数据中的层状介质数据Slayer进行τ‑p变换得到STlayer，利用计算STlayer每个P值对应的地震道的第n个地震数据的预测算子，其中，An1、An2分别为第n个地震数据的预测算子的实部和虚部，frn、fin分别为P值对应的地震道中第1到第n个地震数据的傅里叶变换的实部和虚部，n取1至N的正整数，N为STlayer在时间方向地震数据的个数，ωn为圆频率；步骤103，利用Anmultiple＝An(ωn)‑An‑1(ωn‑1)求取第n个地震数据的多次波预测算子Anmultiple，对Anmultiple做反傅里叶变换，得到时间域多次波预测算子anmultiple，令anmultiple(n)＝0,得到anint；其中，anmultiple为第n个地震数据的多次波预测算子,anint为第n个地震数据的层间多次波预测算子；步骤104，按照公式预测所有Slayer相关的τ‑p域层间多次波；其中，di为第i个地震数据,ε为子波长度，*表示数学上的褶积运算；步骤105，计算层状介质数据中所有的层间多次波预测算子序列并求和记为alayer，对下覆层数据Syu进行τ‑p变换得到STyu，按照公式Myu＝STyu*alayer得到下覆层介质在层状介质中产生的层间多次波Myu；步骤106，将Mlayer和Myu进行相加，得到τ‑p域预测的层间多次波炮集数据，然后进行τ‑p反变换，得到X‑T域层间多次波的炮集数据M；步骤107，在炮集数据中减去计算得到的层间多次波的炮集数据M，得到不含层间多次波的炮集数据。