[发明专利]基于成分相位反相的拖曳瞬变电磁数据工频噪声抑制方法

权利要求书

1.一种基于成分相位反相的拖曳瞬变电磁数据工频噪声抑制方法，其特征在于，包括以下步骤：

S1、获取环境的工频噪声基频频率；

S2、根据环境的工频噪声基频频率确定发射频率；

S3、根据步骤S2中的确定的发射频率计算一个测点内有效采集数据的周期个数；

S4、提取第一个测点的所有有效采集数据，提取正负周期反相数据，相减获得一组数据；

S5、将步骤S4提取的所有数据做叠加平均处理得到第一个测点预处理后的原始数据；

S6、其他测点按照步骤S4和步骤S5进行处理，获得拖曳式数据所有预处理后的原始数据；

步骤S2中，根据探测需求，设计确定发射频率：根据步骤S1取得的数据，确定n的取值，并确定发射频率，确定原则需满足：其中f_tr表示发射频率，f_h表示环境的工频噪声频率，n表示确定发射频率相关参数，采用双极性方波发射，设定占空比为25％，数据连续不间断采集。

2.按照权利要求1所述的方法，其特征在于，步骤S3中具体为，根据预设横向分辨率L与航空飞行器或拖曳车的移动速度v，确定一组内有效采集数据的周期个数M：其中f_tr表示发射频率，即每M个周期的数据等效为一个测点的数据。

3.按照权利要求1所述的方法，其特征在于，步骤S4具体包括：取第一个测点的所有数据和为第一个测点第i周期的正周期数据，是由于第一个测点第i周期的正周期二次场信号第一个测点第i周期的正周期工频噪声和随机噪声构成，表示第一个测点第i+n周期的负周期信号数据，是由于第一个测点第i+n周期的负周期二次场信号第一个测点第i+n周期的负周期工频噪声和随机噪声构成，其中，信号反相，工频噪声同相，进一步做差得到消噪后的数据：

处理后数据y_i(t)滤除了工频噪声数据。

4.按照权利要求1所述的方法，其特征在于，步骤S5中第一个测点的所有数据进行叠加平均处理得到第一个测点预处理后的原始数据包括采用如下的公式：

其中y_i为滤除了工频噪声数据，M为一组内有效采集数据的周期个数。