[发明专利]一种消除地震数据随机噪声的方法和系统

权利要求书

1.一种消除地震数据随机噪声的方法，其特征在于，包括：

采集含有随机噪声的地震数据；

对采集的所述地震数据进行曲波变换；

根据曲波变换所得到的系数确定阈值序列；

定义参数P_i，P_i＝log(||x_i||₁)+log(||r_i||₂)；其中，x_i表示对所述地震数据在曲波域中进行硬阈值运算得到的去噪结果，r_i表示在曲波域采用阈值序列对所述地震数据进行硬阈值去噪后的结果与含噪声的所述地震数据之间的残差，i＝1,2,……N，N表示序列长度；

确定参数P_i为最大值时，所对应的阈值；

根据所述阈值，求取得到去噪后的地震数据。

2.根据权利要求1所述的一种消除地震数据随机噪声的方法，其特征在于，所述阈值序列为一个阈值从大到小逐渐降低的阈值序列，τ₁＞τ₂＞τ₃＞…＞τ_N；其中，N表示序列长度。

3.根据权利要求2所述的一种消除地震数据随机噪声的方法，其特征在于，所述阈值序列中的第一个阈值为所述曲波变换得到系数的绝对值的最大值。

4.根据权利要求3所述的一种消除地震数据随机噪声的方法，其特征在于，所述阈值序列中的最后一个阈值为一个大于0小于1的正实数。

5.根据权利要求1所述的一种消除地震数据随机噪声的方法，其特征在于，x_i为r_i为r_i＝C^*x_i-d_noise；

其中，C表示曲波变换，d_noise为含有随机噪声的地震数据，C^*表示曲波变换的逆变换，TH_τ(·)表示在曲波域进行硬阈值运算，τ表示阈值的大小，τ_i为所述阈值序列，i＝1,2，……N，N表示序列长度。

6.根据权利要求1所述的一种消除地震数据随机噪声的方法，其特征在于，所述根据所述阈值，求取得到去噪后的地震数据，具体为：其中，d_denoise表示去噪后的地震数据，τ_n表示参数P_i为最大值时对应的阈值，C表示曲波变换，C^*表示曲波变换的逆变换。

7.一种消除地震数据随机噪声的系统，其特征在于，包括：

采集模块，用于采集含有随机噪声的地震数据；

曲波变换模块，用于对采集的所述地震数据进行曲波变换；

阈值序列确定模块，用于根据曲波变换所得到的系数确定阈值序列；

参数定义模块，用于定义参数P_i，P_i＝log(||x_i||₁)+log(||r_i||₂)；其中，x_i表示对所述地震数据在曲波域中进行硬阈值运算得到的去噪结果，r_i表示在曲波域采用阈值序列对所述地震数据进行硬阈值去噪后的结果与含噪声的所述地震数据之间的残差，r_i＝C^*x_i-d_noise，i＝1,2,……N，N表示序列长度，C表示曲波变换，d_noise为含有随机噪声的地震数据，C*表示曲波变换的逆变换，TH_τ(·)表示在曲波域进行硬阈值运算，τ表示阈值的大小，τ_i为所述阈值序列；

对应阈值确定模块，用于确定参数P_i为最大值时，所对应的阈值；

求取模块，用于根据所述阈值，求取得到去噪后的地震数据。