[发明专利]去除地震数据中交流电干扰信号的方法

权利要求书

1.一种去除地震数据中交流电干扰信号的方法，所述方法包括： 

（1）获取至少一道地震数据，对每道地震数据进行傅立叶变换，得到至少一道频率域地震数据； 

（2）在设定的频率范围内，从每道频率域地震数据中选取一个最大振幅值的频率域地震数据，将选取的频率域地震数据的振幅值与预设振幅值做比较，将振幅值大于等于预设振幅值的频率域地震数据组成新的频率域地震数据，对新的频率域地震数据进行傅里叶反变换，得到时间域的观测矢量； 

（3）确定最佳分离矩阵，利用分解模型计算交流电干扰信号； 

（4）从时间域的观测矢量中减去交流电干扰信号，得到有效信号。 

2.如权利要求1所述的方法，其中，所述设定的频率范围为48～52Hz。 

3.如权利要求1所述的方法，其中，确定最佳分离矩阵的步骤为， 

（a）对时间域的观测矢量进行预处理； 

（b）利用预处理后的观测矢量，得到分离矩阵； 

（c）对分离矩阵进行归一化处理，得到最佳分离矩阵。 

4.如权利要求3所述的方法，其中，对观测矢量进行预处理的方法为中心化和白化。 

5.如权利要求3所述的方法，其中，步骤（b）包括：利用下面的公式获得J(y)为最大值时的分离矩阵w， 

J(y)∝[E{g(y)}-E{g(y_gauss)}]²

其中，y为随机量，w为分离矩阵，T表示转置，为预处理后的观测矢量，J(y)表示随机量y的负熵，y_gauss为均值为零，且与随机量y具有相同方差的高斯分布的随机量，g(y)和g(y_gauss)为非线性函数， a为常数，E{g(y)}为g(y)的均值运算，E{g(y)}为g(y_gauss)的均值运算。 

6.如权利要求3所述的方法，其中，步骤（c）包括： 

（ⅰ）利用下面的公式获得分离矩阵w的新值w₊， 

其中，w₊为分离矩阵w的新值，g(y)为非线性函数， g′(y)为g(y)的导数，为的均值运算，E{g′(y)}为g′(y)的均值运算； 

（ⅱ）对分离矩阵的新值w₊进行归一化处理，得到最佳分离矩阵w₀，归一化处理的计算公式为， 

其中，‖w₊‖为分离矩阵的新值w₊的模。 

7.如权利要求1所述的方法，其中，所述分解模型为， 

s₂=w₀X， 

其中，X为观测矢量，w₀为最佳分离矩阵，s₂为交流电干扰信号。