一种地震噪声去除方法及装置

摘要

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

主权项

一种地震噪声去除方法，其特征在于包括，步骤101，对炮集数据进行分离，分为层状介质数据S_layer和下覆层数据S_yu；步骤102，对炮集数据中的层状介质数据S_layer进行τ‑p变换得到ST_layer，利用$\left(\begin{array}{c}{A}_n^1\left({\omega }_n\right)\\ A_n^2\left({\omega }_n\right)\end{array}\right)=\left\{\begin{array}{c}\cosh \left(\sqrt{{\mathrm{fr}}^2+{\mathrm{fi}}^2}\right)+\mathrm{fr}·\frac{\sinh \left(\sqrt{{\mathrm{fr}}^2+{\mathrm{fi}}^2}\right)}{\sqrt{{\mathrm{fr}}^2+{\mathrm{fi}}^2}}\\ \mathrm{fi}·\frac{\sinh \left(\sqrt{{\mathrm{fr}}^2+{\mathrm{fi}}^2}\right)}{\sqrt{{\mathrm{fr}}^2+{\mathrm{fi}}^2}}\end{array}\right.$计算ST_layer每个P值对应的地震道的第n个地震数据的预测算子，其中，A_n¹、A_n²分别为第n个地震数据的预测算子的实部和虚部，fr、fi分别为P值对应的地震道中第1到第n个地震数据的傅里叶变换的实部和虚部，n取1至N的正整数，N为ST_layer在时间方向地震数据的个数；步骤103，利用A_n^multiple＝A_n(ω_n)‑A_n‑1(ω_n‑1)求取第n个地震数据的多次波预测算子A_n^multiple，对A_n^multiple做反傅里叶变换，得到时间域多次波预测算子a_n^multiple，令a_n^multiple(n)＝0,得到a_n^int；其中，a_n^multiple为第n个地震数据的多次波预测算子,a_n^int为第n个地震数据的层间多次波预测算子；步骤104，按照公式预测所有S_layer相关的τ‑p域层间多次波；其中，d_i为第i个地震数据,ε为子波长度，*表示数学上的褶积运算；步骤105，计算层状介质数据中所有的层间多次波预测算子序列并求和记为a_layer，对下覆层数据S_yu进行τ‑p变换得到ST_yu，按照公式M_yu＝ST_yu*a_layer得到下覆层介质在层状介质中产生的层间多次波M_yu；步骤106，将M_layer和M_yu进行相加，得到τ‑p域预测的层间多次波炮集数据，然后进行τ‑p反变换，得到X‑T域层间多次波的炮集数据M；步骤107，在炮集数据中减去计算得到的层间多次波M，得到不含层间多次波的炮集数据。