[发明专利]一种大地电磁测深资料时间域压噪的方法及系统

权利要求书

1.一种大地电磁测深资料时间域压噪的方法，其特征在于，包括以下步骤：

读取时间序列的电磁场数据，

将所述时间序列的电磁场数据分隔成若干波段，具体包括，

步骤S11、将时间序列的电磁场数据进行隔点抽样，得到对应的奇数序列和偶数序列；

步骤S12、在所述奇数序列和偶数序列中，若存在其下一顺序位置有对应奇数序列项的偶数序列项，则将该偶数序列项减去与该偶数序列项下一顺序位置对应的奇数序列项的一半以及与该偶数序列项处于相同顺序位置的奇数序列项的一半，得到变换后的偶数序列项；

若存在其下一顺序位置没有对应奇数序列项的偶数序列项，将该偶数序列项减去与其处于相同顺序位置的奇数序列项，得到变换后的偶数序列项；所有变化后的偶数序列项形成了变换后的偶数序列；

在所述奇数序列和偶数序列中，若存在与偶数序列项具有相同顺序位置的奇数序列项，则将处于相同顺序位置的奇数序列项和变换后的偶数序列项对应相加，得到变换后的奇数序列项；

若存在与偶数序列项没有相同顺序位置的奇数序列项，则对该奇数序列项加上与该奇数序列项前一顺序位置对应的变换后的偶数序列项，得到变换后的奇数序列项；所有变换后的奇数序列项形成了变换后的奇数序列；

步骤S13、对变换后的奇数序列再次进行隔点抽样，得到对应的奇数序列和偶数序列；

步骤S14、执行步骤S12；

步骤S15、重复执行步骤S13～S14，直至到达预设隔点抽样次数；形成的所有奇数序列和偶数序列构成了若干波段，形成电磁场数据波段；

获取噪音的频率范围，将落在噪音频率范围内的电磁场数据波段对应的数据设置为零；

将所有电磁场数据波段重构为一个完整的电磁场数据。

2.根据权利要求1所述的大地电磁测深资料时间域压噪的方法，其特征在于，所述将所有电磁场数据波段重构为一个的完整电磁场数据，具体包括：

将最后得到的奇数序列和偶数序列执行下述步骤S21，

步骤S21、若存在与偶数序列项具有相同顺序位置的奇数序列项，则将该奇数序列项减去对应顺序位置的偶数序列项，得到重构后奇数序列项；

若存在与偶数序列项没有相同顺序位置的奇数序列项，则对该奇数序列项减去与该奇数序列项前一顺序位置对应的偶数序列项，得到重构后的奇数序列项；所有重构后的奇数序列项形成了重构后的奇数序列；

若存在其下一顺序位置有对应奇数序列项的偶数序列项，则将该偶数序列项加上与该偶数序列项下一顺序位置对应的重构后的奇数序列项的一半以及该偶数序列项相同顺序位置的重构后的奇数序列项的一半，得到重构后偶数序列项；

若存在其下一顺序位置没有对应奇数序列项的偶数序列项，将偶数序列项加上对应顺序位置的重构后的奇数序列项，得到重构后的偶数序列项；所有重构后的偶数序列项形成了重构后偶数序列；

步骤S22、将重构后的奇数序列和偶数序列，按照隔点抽样的逆过程进行组合，得到重构后的序列；

步骤S23、将重构后的序列作为一个奇数序列，将所述奇数序列以及与其对应的偶数序列，执行步骤S21；

步骤S24，重复执行步骤S22～S23，直至步骤S22中得到的重构后的序列没有对应的偶数序列为止；所述得到的重构后的序列即为一个完整的电磁场数据。

3.根据权利要求1所述的大地电磁测深资料时间域压噪的方法，其特征在于，所述获取噪音的频率范围，具体包括：依据电磁测点的视电阻曲线、相位曲线或者信噪比曲线获取噪音的频率范围。