[发明专利]混叠噪声压制方法及装置

权利要求书

1.一种混叠噪声压制方法，其特征在于，包括：

对混叠数据进行动校正；

采用指定步长对动校正后的数据进行中值滤波；

对中值滤波后的数据进行反动校正，获得第一数据；

基于复曲波域的阈值去噪方法从残余数据中提取有效信号，获得第二数据；所述残余数据为所述混叠数据减去所述第一数据后的数据；

将所述第一数据和所述第二数据相加，获得有效数据；

确定所述混叠数据的第一伪分离数据，并根据预设的混叠算子获取所述有效数据的第二伪分离数据；

确定所述第一伪分离数据与所述第二伪分离数据的差值；

以所述差值作为新的混叠数据迭代计算新的差值，并将所述有效数据与所述新的差值相加获得新的有效数据，依此递推直至当前迭代出的有效数据满足预设信噪比为止；

所述确定所述混叠数据的第一伪分离数据，包括：

根据公式m'＝Γ^Hd_bl获取所述混叠数据的第一伪分离数据；

其中，m'为第一伪分离数据，d_bl为混叠数据，Γ^H为混叠算子Γ的共轭转置；

所述根据预设的混叠算子获取所述有效数据的第二伪分离数据，包括：

根据公式m＝Γ^HΓm获取所述有效数据的第二伪分离数据；

其中，m为第二伪分离数据，m为有效数据,Γ为混叠算子，Γ^H为混叠算子Γ的共轭转置。

2.如权利要求1所述的混叠噪声压制方法，其特征在于，在迭代过程中，所述中值滤波中的指定步长及所述复曲波域的阈值去噪方法中的阈值逐渐减小。

3.如权利要求1所述的混叠噪声压制方法，其特征在于，所述采用指定步长对动校正后的数据进行中值滤波，包括：

采用指定步长并根据公式对动校正后的数据进行中值滤波；

其中，u_m为数据集的中值，u_i和u_j分别为数据集中的第i个和第j个数据，W为指定的步长，p为定义距离时所用的模。

4.如权利要求1所述的混叠噪声压制方法，其特征在于，所述基于复曲波域的阈值去噪方法从残余数据中提取有效信号，获得第二数据，包括：

根据公式m_n+1＝C-¹T_λC[m_n+γΓ^H(d_bl-Γm_n)]从残余数据中提取有效信号，获得第二数据；

其中，C和C-¹表示复曲波变换的正反变换对，m_n表示第n次迭代结果，Γ为混叠算子，Γ^H为混叠算子Γ的共轭转置，γ为标量系数，T_λ为阈值函数，λ为所取阈值。

5.如权利要求4所述的混叠噪声压制方法，其特征在于，所述阈值函数包括：

T_λ(x,λ)＝sign(x)·max(0,|x|-λ)；

其中，T_λ(x,λ)为阈值函数，sign(x)为符号函数，x为残余数据在复曲波域的变换系数。

6.如权利要求1所述的混叠噪声压制方法，其特征在于，所述有效数据的信噪比通过以下公式确定：

其中，SNR_n为有效数据m_n的信噪比，d_bl为混叠数据，Γ为混叠算子，Γ^H为混叠算子Γ的共轭转置。