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

权利要求书

1.一种地震噪声去除方法，其特征在于包括，

步骤101，对炮集数据进行分离，分为层状介质数据S_layer和下覆层数据S_yu；

步骤102，对炮集数据中的层状介质数据S_layer进行τ-p变换得到ST_layer，利用An1(ωn)An2(ωn)=cosh(fr2+fi2)+fr·sinh(fr2+fi2)fr2+fi2fi·sinh(fr2+fi2)fr2+fi2]]>计算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，得到不含层间多次波的炮集数据。