[发明专利]一种预测地震波中的多次波和一次波信号的方法

权利要求书

1.一种预测地震波中的多次波信号的方法，包括以下步骤：

1)用如下式(9)的目标函数优化滤波器系数a(i)：

max J(a(i))＝[E{G(ε₀(n))}-E{G(v(n))}]²，i＝1，2，…，m    (9)；

2)将步骤1)中得到的优化滤波器系数a(i)代入下式(1)，得到多次波信号

式(9)和式(1)中，ε₀(n)是误差信号ε(n)的归一化信号，v(n)是标准高斯信号，G(·)是非线性函数，J(a(i))是用来衡量ε₀(n)的非高斯性程度的变量，q是地震观测信号的预测步长，i是滤波器系数的序号，n是地震观测信号的序号，m是滤波器的阶数，y(n-q-i+1)是地震观测信号；

所述G(·)为下式(10)中的函数：

G(x)＝-exp(-0.5x²)    (10)。

2.如权利要求1所述的方法，其特征在于：所述用式(9)的目标函数优化滤波器系数a(i)的方法为迭代求解如下方程组：

R(i,s)=Σn=q+mN[exp(-0.5ϵ02(n))y(n-i-q+1)y(n-s-q+1)],]]>i＝1，…，m，s＝1，…，m(12)；

b(i)=Σn=q+mN[exp(-0.5ϵ02(n))y(n-i-q+1)y(n)],]]>i＝1，…，m    (13)；

式(11)、式(12)和式(13)中，m为滤波器的阶数，N为地震观测信号y的长度，q为地震观测信号的预测步长。

3.一种预测地震波中的一次波信号的方法，是将多次波信号代入下式(2)得到一次波信号ε(n)；

ϵ(n)=y(n)-y^(n)=y(n)-Σi=1ma(i)y(n-q-i+1)---(2);]]>

式(2)中，所述y(n)为地震观测信号；所述是通过权利要求1或2所述的方法计算得到的。