[发明专利]去除地震信号中的随机噪声的方法和装置

权利要求书

1.一种去除地震信号中的随机噪声的方法，包括：

对地震信号进行多尺度Contourlet变换，得到在多尺度多方向的Contourlet系数，Contourlet系数c_l,k(m)表示为：

其中，l表示尺度编号，k表示方向编号，0≤k＜2^l-1；m、n表示在尺度l的二维离散希尔伯特空间L²(Z²)的位置信息，Z²表示二维离散信号的整数网格，L²表示平方可积空间；b_l(n)表示在尺度l上进行DFB分解前的高频子带信号；S_l,k表示在尺度l、方向k的采样系数，h_l,k(S_l,km-n)表示在尺度l的DFB分解中在方向k的楔形滤波器系数，是二维离散希尔伯特空间L²(Z²)中的一组基；

得到在每一尺度的每个方向，基于在该尺度该方向上的强有效信号对应的第一子带系数P1(c_l,k,m)：

其中H表示预设的硬阈值；

基于第一子带系数P1(c_l,k,m)重构所述地震信号的高频子带信号。

2.根据权利要求1所述的去除地震信号中的随机噪声的方法，其中，基于第一子带系数P1(c_l,k,m)重构所述地震信号的高频子带信号还包括：

每一尺度的每一方向包括多个道，每一道包括多个点，每个点对应于一个Contourlet系数c_l,k(m)，相应地，每个点对应于一个第一子带系数P1(c_l,k,m)，在每一尺度的每个方向，以互相关值为判断标准对各个点对应的第一子带系数P1(c_l,k,m)进行倾角扫描，以及基于每个点在各个倾角方向上的互相关值判断该点是否是强有效信号；

在每一尺度的每个方向，基于在该尺度该方向上的强有效信号对应的第一子带系数P1(c_l,k,m)确定在该尺度该方向的软阈值；

得到第二子带系数P2(c_l,k,m)：

其中F_l,k表示在尺度l、方向k的软阈值；

基于第二子带系数P2(c_l,k,m)重构所述地震信号的高频子带信号。

3.根据权利要求2所述的去除地震信号中的随机噪声的方法，其中，基于每个点在各个倾角方向上的互相关值判断该点是否是强有效信号包括：

针对每个点，得到其在各个倾角方向上的互相关值中的最大值和第二大值，如果该最大值与该第二大值间的差大于互相关阈值，则判断该点是强有效信号。

4.根据权利要求2所述的去除地震信号中的随机噪声的方法，其中，基于在该尺度该方向上的强有效信号对应的第一子带系数P1(c_l,k,m)确定在该尺度该方向的软阈值包括：

基于在该尺度该方向上的各个强有效信号对应的各个第一子带系数P1(c_l,k,m)的绝对值中的最小值确定在该尺度该方向的软阈值。

5.根据权利要求4所述的去除地震信号中的随机噪声的方法，其中，所述软阈值为在该尺度该方向上的各个强有效信号对应的各个第一子带系数P1(c_l,k,m)的绝对值中的最小值的10％。