[发明专利]一种基于站间传递函数的大地电磁噪声压制方法

摘要

本发明专利提供一种基于站间传递函数的大地电磁噪声压制方法，该方法通过估算同步测点之间电磁场的传递函数，合成本地测点时间序列信号，从而去除噪声，其实现步骤为：1.对多个大地电磁测点时间序列数据进行同步；2.估算站间电磁场之间的传递函数；3.利用参考测点数据和站间传递函数，合成本地测点的时间序列信号；4.识别本地测点噪声所在的位置；5.用合成信号替换噪声，得到去噪后的新数据。通过上述实现步骤，本发明可以有效压制强相关噪声，抑制飞点，减小估计误差。

主权项

1.一种基于站间传递函数的大地电磁噪声压制方法，其特征在于，能够利用站间传递函数对本地测点进行噪声压制，包括如下步骤：(1)对多个大地电磁测点时间序列数据进行同步；(2)估算本地测点与参考测点电磁场之间的站间传递函数，站间传递函数包括：本地测点水平电场与参考测点水平磁场之间的传递函数，本地测点水平磁场与参考测点水平磁场之间的传递函数，本地测点水平电场与参考测点水平电场之间的传递函数，本地测点水平磁场与参考测点水平电场之间的传递函数，本地测点垂直磁场与参考测点水平磁场之间的传递函数；(3)利用站间传递函数和参考测点数据，合成本地测点的时间序列信号：先对参考测点时间序列数据进行快速傅里叶变换，得到参考测点电场和磁场的频谱，再利用参考测点电场或磁场的频谱，结合站间传递函数，合成本地测点电场和磁场的频谱，对本地测点电场频谱和磁场频谱分别进行傅里叶反变换，得到本地测点电场和磁场的时间序列合成信号；(4)识别本地测点噪声所在的位置；(5)用合成信号替换噪声数据，得到去噪后的新数据，达到去噪的目的，替换时在窗口的左边界和右边界进行加权处理，避免引入阶跃噪声，具体采用如下形式：左边界：右边界：其中d_new表示本地测点的新数据，d_obs表示本地测点的观测数据，d_syn表示本地测点的合成数据，L为窗口边界的长度，L小于窗口长度N，为加权的权重。