技术特征：
1.一种大地电磁强噪声时段压制方法，其特征在于，包括：获取大地电磁数据阵列；所述大地电磁数据阵列包括5道数据，第1道数据为南北向的电道测量值，第2道数据为东西向的电道测量值，第3道数据为南北向的磁道测量值，第4道数据为东西向的磁道测量值，第5道数据为垂向磁道测量值；基于所述第1道数据和所述第2道数据构造复数道序列；对所述复数道序列进行hilbert变换，并获取所述复数道序列的奇部序列和偶部序列；对所述奇部序列和所述偶部序列进行短时傅里叶变换，获取瞬时频谱；基于所述瞬时频谱计算极化椭圆率、极化方向以及正交电磁分量的相位差；基于所述极化椭圆率、所述极化方向以及所述正交电磁分量的相位差，确定极化性质时频谱；所述极化性质时频谱包括极化椭圆率时频谱、极化方向时频谱以及正交电磁场分量相位差时频谱；基于所述极化性质时频谱确定强噪声时段并进行删除，得到去噪后的大地电磁数据阵列。2.根据权利要求1所述的大地电磁强噪声时段压制方法，其特征在于，所述复数道序列的构造公式如下：c(t)＝e
x
(t)+ie
y
(t)其中，c(t)为t时间点的复数道序列，i为虚数单位，e
x
为第1道数据，e
y
为第2道数据。3.根据权利要求1所述的大地电磁强噪声时段压制方法，其特征在于，所述瞬时频谱的计算公式如下：计算公式如下：其中，为偶部序列c
+
对应的瞬时频谱，为偶部序列c
－
对应的瞬时频谱，t＝1,2,3,...,n，为不同时间点，m为瞬时频谱的时间因子，n为瞬时频谱的频率因子，g*为窗函数g的共轭，i为虚数单位。4.根据权利要求3所述的大地电磁强噪声时段压制方法，其特征在于，所述极化椭圆率ε的计算公式如下：5.根据权利要求3所述的大地电磁强噪声时段压制方法，其特征在于，所述极化方向θ的计算公式如下：6.根据权利要求3所述的大地电磁强噪声时段压制方法，其特征在于，所述正交电磁分量的相位差δφ的计算公式如下：
7.一种大地电磁强噪声时段压制系统，其特征在于，包括：大地电磁数据阵列获取模块，用于获取大地电磁数据阵列；所述大地电磁数据阵列包括5道数据，第1道数据为南北向的电道测量值，第2道数据为东西向的电道测量值，第3道数据为南北向的磁道测量值，第4道数据为东西向的磁道测量值，第5道数据为垂向磁道测量值；复数道序列构造模块，用于基于所述第1道数据和所述第2道数据构造复数道序列；奇部序列和偶部序列获取模块，用于对所述复数道序列进行hilbert变换，并获取所述复数道序列的奇部序列和偶部序列；瞬时频谱获取模块，用于对所述奇部序列和所述偶部序列进行短时傅里叶变换，获取瞬时频谱；计算模块，用于基于所述瞬时频谱计算极化椭圆率、极化方向以及正交电磁分量的相位差；极化性质时频谱确定模块，用于基于所述极化椭圆率、所述极化方向以及所述正交电磁分量的相位差，确定极化性质时频谱；所述极化性质时频谱包括极化椭圆率时频谱、极化方向时频谱以及正交电磁场分量相位差时频谱；去噪模块，用于基于所述极化性质时频谱确定强噪声时段并进行删除，得到去噪后的大地电磁数据阵列。
技术总结
本发明公开了一种大地电磁强噪声时段压制方法及系统。该方法包括：获取大地电磁数据阵列；基于大地电磁数据阵列中的第1道数据和第2道数据构造复数道序列；对复数道序列进行Hilbert变换，并获取复数道序列的奇部序列和偶部序列；对奇部序列和偶部序列进行短时傅里叶变换，获取瞬时频谱；基于瞬时频谱计算极化椭圆率、极化方向以及正交电磁分量的相位差；基于极化椭圆率、极化方向以及正交电磁分量的相位差，确定极化性质时频谱；基于极化性质时频谱确定强噪声时段，并进行删除得到去噪后的大地电磁数据阵列。本发明能够针对性地去除大地电磁时间序列中的噪声干扰，同时适用性很强，适用于绝大多数的人文噪声压制且不受噪声复杂波形的制约。复杂波形的制约。复杂波形的制约。
技术研发人员：康建强 韩江涛 刘立家 郭振宇
受保护的技术使用者：吉林大学
技术研发日：2022.06.29
技术公布日：2022/9/30