发明名称 一种预测地震波中的多次波和一次波信号的方法
摘要 本发明公开了一种预测地震波中的多次波和一次波信号的方法,属于地震观测信号处理技术领域。本发明所提供的预测地震波中的多次波信号的方法,包括以下步骤:1)用式(9)的目标函数优化滤波器系数a(i);2)将步骤1)中得到的优化滤波器系数a(i)代入式(1),得到多次波信号。用地震观测信号减去所述多次波的信号还可以得到地震波中的一次波信号。本发明克服了传统的预测方法需要假设预测多次波和一次波是正交的限制,更有效地消除了地震信号中的多次波。<img file="200810117757.3_ab_0.GIF" wi="242" he="56" />(1);maxJ(a(i))=[E{G(ε<sub>0</sub>(n))}-E{G(v(n))}]<sup>2</sup>,i=1,2,…,m,(9)。
申请公布号 CN101334482B 申请公布日期 2011.05.25
申请号 CN200810117757.3 申请日期 2008.08.05
申请人 中国海洋石油总公司;中海石油研究中心;清华大学 发明人 赵伟;陆文凯
分类号 G01V1/28(2006.01)I 主分类号 G01V1/28(2006.01)I
代理机构 北京纪凯知识产权代理有限公司 11245 代理人 关畅;任凤华
主权项 1.一种预测地震波中的多次波信号的方法,包括以下步骤:1)用如下式(9)的目标函数优化滤波器系数a(i):max J(a(i))=[E{G(ε<sub>0</sub>(n))}-E{G(v(n))}]<sup>2</sup>,i=1,2,…,m    (9);2)将步骤1)中得到的优化滤波器系数a(i)代入下式(1),得到多次波信号<img file="FSB00000322104300011.GIF" wi="155" he="65" /><img file="FSB00000322104300012.GIF" wi="1643" he="117" />式(9)和式(1)中,ε<sub>0</sub>(n)是误差信号ε(n)的归一化信号,v(n)是标准高斯信号,G(·)是非线性函数,J(a(i))是用来衡量ε<sub>0</sub>(n)的非高斯性程度的变量,q是地震观测信号的预测步长,i是滤波器系数的序号,n是地震观测信号的序号,m是滤波器的阶数,y(n-q-i+1)是地震观测信号;所述G(·)为下式(10)中的函数:G(x)=-exp(-0.5x<sup>2</sup>)    (10)。
地址 100027 北京市东城区东直门外小街6号