一种电磁信号强噪声去除方法技术

本发明专利技术提供的一种电磁信号强噪声去除方法，包括以下步骤：(1)设置两套发射频率表；(2)采集站接收两套发射信号；(3)去除数据中的非平稳随机噪声，即去除随机噪声E

A method of removing strong noise of electromagnetic signal

【技术实现步骤摘要】
一种电磁信号强噪声去除方法
本专利技术属于一种电磁信号强噪声去除方法，具体为电磁勘探信号处理

技术介绍
大地电磁测深法(MT)是上个世纪50年代初被提出以来的，该方法已广泛应用于石油、金属矿、地下水等资源的勘探中。但是由于MT方法利用天然源信号，因此其抗干扰能力弱，测量精度比较低。为了克服MT方法的这些缺点，地球物理学家提出了采用人工发射源作为信号的可控源电磁法(CSEM)，这些方法包括可控源音频大地电磁法CSAMT，极低频探地电磁法WEM等。这些方法的显著特点是采用人工可以控制的发射源发射不同频率的电磁信号，接收到的场值具有较高的信噪比，提高了野外工作效率。CSEM法虽然具有相对较高的信噪比，但其探测应用仍然受到了野外电磁干扰的影响。电磁噪声主要是由人工活动产生或者人工产生的，包括随机噪声及平稳噪声，这些噪声的频率分布范围广，并且随着国民经济和现代化建设的发展，这种干扰日趋严重。在强电磁干扰区，电磁法往往无法施工应用，即使勉强工作，也会得到错误的探测结果，可能导致巨大的经济损失。噪声的处理与压制问题已经成为CSEM探测方法得以成功实施的关键所在。现有技术设置一套频率表进行电磁信号的发射，信号接收端通过增加观测时长，对同频点进行多次观测，通过将采集得到的多次观测数据进行叠加处理，实现噪声压制的目的。叠加处理的具体实施方式为对观测到的时间序列数据分段分别进行处理，并得到同频点对应的多个电磁场幅值，对多个数据进行叠加、平均、筛选等。比如对于1Hz的信号，观测10分钟可得到含有600个周期信号的时间序列数据。如果对这些数据分10段进行处理可以得到10个1Hz的场值，对这10个数据进行筛选及叠加平均可以得到最终的1Hz对应的场值结果。现有技术仅能够去除信号中随机非平稳噪声的干扰。实际上，野外测量中遇到的噪声类型包括随机噪声与平稳噪声。随机噪声包括高精度传感器仪器系统的振动噪声及测区附近的各种用电设备开启关闭及运行产生的脉冲干扰等。而平稳噪声主要来源于周围环境中的高压电缆、测区存在的工频用电器以及工频地电流等，较强的稳定的50Hz及谐波干扰是CSEM方法中的主要强干扰来源。对于这些广泛存在的稳定电磁噪声，现有技术无能为力。尤其随着工业化程度的增强，这些平稳噪声的干扰越来越严重，直接导致CSEM方法探测能力的失效。总之，现有野外施工与处理技术忽略了平稳噪声的影响，采用单次发射与接收端多次叠加的方式，进行信号的接收与处理，将处理后的数据直接作为有效信号使用。而实际上信号里面包含了平稳噪声及有效噪声，这种忽略平稳噪声的方式在环境噪声较大时会导致处理结果误差大，甚至可能产生错误的探测结果。
技术实现思路
本专利技术提出了一种去除电磁噪声的方法，尤其对广泛存在的平稳强噪声更加有效。根据本专利技术提出的去噪方法，可计算得到高精度的有效信号，去除背景平稳噪声，实现对强电磁干扰的高效压制与去除，大幅提高CSEM方法的数据质量，解决此类方法在高干扰区施工的难题。具体的技术方案为：一种电磁信号强噪声去除方法，包括以下步骤：(1)设置两套发射频率表按照现有的CSEM野外施工流程，完成发射机及发射天线电极的布设；设置完成第一套发射机发射频率表Table1；设置完成第二套发射机的发射频率表Table2，其中发射频点序列F2＝F1，各频点的发射持续时间T2＝T1，各频点的发射电流大小I2＝a×I1，其中系数a≠1；(2)采集站接收两套发射信号在噪声干扰区完成阵列式电磁信号接收机的布设，准备接收发射机发射的电磁信号；发射机首先按照频率表Table1发射信号，各接收机完成频率表Table1的电磁信号采集；发射机按照频率表Table2发射信号，阵列式接收机完成频率表Table2的电磁信号采集；可控源电磁法在野外观测到的场值包含需要的有效数据、平稳噪声及非平稳随机噪声三部分，如下公式所示：其中Ei1和Ei2分别为接收机第1次和第2次观测到的第i个频点的电场或磁场幅值，E有效为我们需要的真正的携带地下信息的数据；设第二次的发射电流大小为第一次发射电流大小的a倍，则Ei2有效＝a×Ei1有效；Ei平稳噪声是稳定存在的环境噪声，其值在不同的观测时间是不变的；(3)去除数据中的非平稳随机噪声去除随机噪声Ei1随机噪声和Ei2随机噪声，得到下式(4)计算得到真实信号及环境噪声各测点对各频点进行计算，得到此测点此频点的有效信号及平稳噪声大小；根据电磁探测原理，因为第二次发射电流为第一次发射电流的a倍，则第二次接收的有效信号也应为第一次有效信号的a倍：得到：其中Ei1及Ei2为由接收机采集得到，Ei有效及Ei平稳噪声为未知数，可以求解得到。本专利技术的原理：电场的有效数据可以通过如下公式正演模拟计算得到其中Ex为电场的x方向分量，发射电矩P＝I×L，其中I为发射电流，L为发射天线的长度。ρ为大地视电阻率(是未知的，是在实际工作中我们想要得到的)，r为接收测点位置到发射源中心的距离，x和y为接收测点位置的坐标；k为波数，ω＝2πf，f为发射频率；由公式1可知，某接收位置采集到的电磁场幅值是与发射源的电矩P成正比的。在野外天线长度不变的条件下，可通过改变发射电流的大小来改变发射电矩P的大小，从而改变电场的大小。基于此原理，本专利技术通过改变发射电流的大小，进行两次信号发射并采集相应的信号的方式，建立起包含两个未知数的两个方程，从而可得到真实的有效场值以及环境平稳噪声。CSEM法近年得到了广泛的应用，但是其一直受到环境电磁噪声影响的困扰。电磁噪声包括平稳噪声与随机噪声，现有的方法主要采用数据叠加等方式进行噪声压制，此种方法仅对于随机噪声具有效果，但是无法去除广泛存在的平稳强干扰信号，这会导致CSEM法在强干扰区探测失败。本专利技术提出了一种采用不同发射电流进行两次发射的方式，可以利用采集到的两套数据处理得到有效信号，高效压制环境噪声。本专利技术仅通过两次发射即可去除长期制约着电磁法应用的强噪声干扰，得到测点处的有效电磁场，解决此重要问题。本专利技术提供的一种电磁信号强噪声去除方法，不仅野外施工操作简单方便，而且对于去除强干扰噪声十分有效，可以得到非常干净的有效信号，大幅提高CSEM的数据质量。由于CSEM方法的探测精度随着干扰的减小而逐渐提高，本专利技术的数据去噪方法可为提高CSEM法的探测精度提供重要保证。具体实施方式结合实施例说明本专利技术的具体技术方案。本实例中设大地为电阻率为1000Ωm的均匀大地，发射极长度为2km，接收点距离发射机中点为10km。本实例中设发射机共发射24个频率点，分别为(单位Hz)：1，2，4，7680，5120，3840，2560，1920，1280，1024，711.11，512，355.56，256，177.78，128，88.89，64，44.44，32，22.22，16本文档来自技高网...

【技术保护点】
1.一种电磁信号强噪声去除方法，其特征在于，包括以下步骤：/n(1)设置两套发射频率表/n按照现有的CSEM野外施工流程，完成发射机及发射天线电极的布设；设置完成第一套发射机发射频率表Table1；设置完成第二套发射机的发射频率表Table2，其中发射频点序列F2＝F1，各频点的发射持续时间T2＝T1，各频点的发射电流大小I2＝a×I1，其中系数a≠1；/n(2)采集站接收两套发射信号/n在噪声干扰区完成阵列式电磁信号接收机的布设，准备接收发射机发射的电磁信号；发射机首先按照频率表Table1发射信号，各接收机完成频率表Table1的电磁信号采集；发射机按照频率表Table2发射信号，阵列式接收机完成频率表Table2的电磁信号采集；/n可控源电磁法在野外观测到的场值包含需要的有效数据、平稳噪声及非平稳随机噪声三部分，如下公式所示：/n

【技术特征摘要】
1.一种电磁信号强噪声去除方法，其特征在于，包括以下步骤：
(1)设置两套发射频率表
按照现有的CSEM野外施工流程，完成发射机及发射天线电极的布设；设置完成第一套发射机发射频率表Table1；设置完成第二套发射机的发射频率表Table2，其中发射频点序列F2＝F1，各频点的发射持续时间T2＝T1，各频点的发射电流大小I2＝a×I1，其中系数a≠1；
(2)采集站接收两套发射信号
在噪声干扰区完成阵列式电磁信号接收机的布设，准备接收发射机发射的电磁信号；发射机首先按照频率表Table1发射信号，各接收机完成频率表Table1的电磁信号采集；发射机按照频率表Table2发射信号，阵列式接收机完成频率表Table2的电磁信号采集；
可控源电磁法在野外观测到的场值包含需要的有效数据、平稳噪声及非平稳随机噪声三部分，如下公式所示：


...

【专利技术属性】
技术研发人员：付长民，底青云，王妙月，雷达，王若，安志国，王啸天，
申请(专利权)人：中国科学院地质与地球物理研究所，
类型：发明
国别省市：北京;11