CN102721982A - 基于扩频编码技术的地面电磁勘探方法及其探测系统 - Google Patents

Abstract

一种基于扩频编码技术的地面电磁勘探方法及其探测系统，通过发送电极向地下供入以一定的编码序列跳变的电流信号作为人工激励源，接收机接收电磁场响应信息，同时同步记录发送的电流信号和接收机所在位置的地电响应信息，经过相关辨识方法，得到包含有地电阻率的大地系统响应，最终得出地电阻率分布特征。应用本发明勘探方法的勘探系统，由发送机发送整个预设频段范围的单极性或双极性的正弦波或方波信号，发送序列按预先设定好的频率图案跳变，由相关辨识检测方法去除不相关噪声。

Description

基于扩频编码技术的地面电磁勘探方法及其探测系统

技术领域

本发明涉及一种地球物理勘探新方法，特别涉及一种扩频编码技术的地面电磁勘探方法及探测系统。

背景技术

频率域电磁法勘探为地球物理勘探的重要手段。频率域电磁法的原理是根据电磁信号在地下的穿透深度与地下介质的导电性质以及所用电磁波的频率有密切的关系。当地下介质的电导率不变时，低频的电磁波穿透的更深，了解到地下信息更多，高频则相反。频率域电法在寻找地下金属矿藏、地质灾害探测、水资源分布、煤炭、油气等资源勘查方面都有着十分广泛的应用。

目前比较常用的频率域电法主要有：变频法、奇次谐波法、双频电流法、可控源音频大地电磁法（CSAMT）、aⁿ序列伪随机法等。传统的变频激电法需要逐次改变频率，一个一个频率的供电和测量，所供电流的波形很难保持完全不变，接收时所受的干扰程度也不同，精度和效率低；奇次谐波法比起变频法有所进步，能够一次供电进行多个频率谐波的测量，然而奇次谐波法谐波强度随谐波次数成反比衰减，且谐波频率分布不均匀；双频激电法将两种频率的电流合成，同时供入地下，一次性同时接收两种频率的信号，其成功克服了变频激电的不足，与奇次谐波法相比，它的两个频率成分强度完全相等。但对于需要提取多个频率地电信息的频谱激电和电磁法来说，仍然难以一次获得较为完整的激电频谱；可控源音频大地电磁法克服了大地电磁（MT）法场源随机性的缺点，但其沿用MT法卡尼亚公式要求在“远区”进行测量，达不到远区公式不成立。

aⁿ序列伪随机法实现了工作效率高、观测精度高、仪器轻便、观测参数丰富等优点。aⁿ序列伪随机信号电法是根据伪随机信号编码的数学原理，用-1，0和1三个码元分别表示电流I=-I，0，I，将n个不同频率的电流，组合为含有n个主要频率成分的合成电流，同时供入地下。一次观测可以从地下提取n个不同频率的响应。这种方法由于探测信号含有多个频率成分，可同时测量，一定程度上提高了工作效率。以2n伪随机信号为例，其场源电流中含有的主频率按2n步进，主频率的个数为3、5、7…157…等。根据频率测深的原理，探测信号频率的不连续性会造成地质探测深度的不连续，分辨率较低，进而影响地质探测的效果；另外，当外界信号与探测信号中某一频率重合或接近时，极易影响整个波形的探测结果。为此，在采用这种方法进行探测时，人们不得不加大探测发送机功率，以提高接收信号的信噪比。但功率的加大造成其体积和重量的成倍的增长，不利于野外和山区的勘测。

对于频率域电法而言，我们关注的是大地对供电波形的响应，理想频率域场源信号应该具有频率范围宽、丰富的频谱信息、频谱精细化程度高、主频能量分布均匀、信噪比高等特点。

就发送而言，目前的人工场源电磁探测方法发展的主要瓶颈为发送序列和发送功率两方面，主要表现在发送序列的单一化和能量分布不均。已有的电磁探测法存在着发送序列复杂度和频谱细化程度较低，信号功率谱密度低于噪声功率谱密度、频谱密度与主频能量矛盾等问题，制约着电磁探测法的发展，使得单次观测探测效率低，信噪比的提高依赖于系统能量的提高；就信号检测方法而言，目前的电磁测探方法为了提高信噪比，均采用了多次叠加的观测方法，抗干扰能力低，需要多次重复测量来获得可靠信号。随着对勘探深度和纵向分辨率要求的提高，以及现代社会人文干扰日益加剧，对电磁勘探检测方法也提出了更高的要求。

英国爱丁堡大学Hobbs等人提出了用伪随机编码代替阶跃电流，国内学者也有关于m序列及逆重复m序列伪随机电磁法的研究。但此处提到的m序列及逆重复m序列电磁法所使用的序列复杂度不够，频谱均是在频带内等间距分布，其伪随机谱线的疏密对辨识准确度影响很大，所使用的伪码序列的限制使得信号功率谱密度低于噪声功率谱密度；且对于m序列和逆m序列，半周期数T的选取对相关辨识算法精度影响很大，T和信号幅度I越大，误差越小，T越大，谱线越密，但T太大，信号频带变宽，主频能量变小，各种突变干扰对辨识脉冲响应的不良影响也越大，欲提高系统的信噪比，就需要加大系统功率。即所提出的方法，由于受到扩展信号频谱的限制，仍然在走增加信号功率，减少噪声，提高信噪比的路。

就信号检测方法而言，目前的电磁测探方法为了提高信噪比，均采用了多次叠加的观测方法，抗干扰能力低，需要很多次重复测量来获得可靠信号。

发明内容

为了克服上述现有方法的不足，本发明提出一种扩频编码技术的地面电磁勘探方法及其探测系统。本发明采用的编码序列频率范围宽、频谱信息丰富、频谱精细化程度高、主频能量分布均匀，可克服发送序列的单一化和能量分布不均，单次观测探测效率低，信噪比的提高依赖于系统能量的提高等问题。本发明同步记录发送的电流信号和接收机所在位置的地电响应信息，采用相关辨识计算，无用信号的干扰码序列不起作用。这样总能避开噪声干扰，提高探测的分辨率、抗多频干扰能力，提高探测效率。

本发明通过发送电极向地下供入以扩频编码序列跳变的电流信号作为人工激励源，接收机接收电磁场响应信息，同时配合GPS同步技术记录发送的电流信号和接收机所在位置的地电响应信息，经过相关辨识方法，得到大地的频谱响应，计算出视电阻率、视极化率、视复电阻率等地球物理信息。本发明为地面电磁探测的精细化、高分辨率测深提供新的思路，同时拓宽了扩频技术的应用领域。随着新方法的不断进步和完善，预测本发明将在资源探测领域发挥重要作用。

本发明方法的步骤如下：

（1）发送机通过发送电极向地下供入以扩频编码序列跳变的电流信号作为人工激励源，同时记录发送电流信号的频率、幅度、相位等相关信息参数；

（2）在距离人工激励场源数米到数十公里范围内接收电极处和场源附近的电磁场响应信息；

（3）发送信号的采集与接收信号的采集严格同步，同步方式采用GPS同步来实现；

（4）利用相关算法，将接收端接收到的信号与发送序列做相关运算，去除不相关噪声，可提取大地系统冲激响应；

（5）大地系统冲激响应为地电阻率、收发距及频率或时间的函数，可求得相关地球物理信息。

应用本发明方法的勘探系统主要包括发送机与接收机。

所述的发送机用于发送扩频编码信号。所述的发送机包括供电电源、快熔、升压电路、整流滤波、主控单元FPGA，即现场可编程门阵列、隔离驱动单元、保护单元、逆变单元、电压电流检测模块、A/D转换模块、输入控制模块、显示存储模块、GPS同步模块以及隔离电源组；所述的主控单元FPGA包括编码序列频率-波形合成单元、数据采集单元、USB接口程序；发送机的连接方式为：供电电源的输出端连接快熔后输入给升压电路，升压电路输出作为整流滤波单元的输入，整流滤波单元的输出电压输入给逆变单元；主控单元FPGA的编码序列频率—波形合成单元输出跳频编码序列，输入给隔离驱动电路，隔离驱动电路的输出为逆变单元开关管的输入驱动信号，主控单元FPGA与控制保护单元相连，保护单元与逆变单元相连接，用于对逆变单元进行保护；逆变单元与电压电流检测单元相连接，电压电流检测的输出作为A/D转换的输入，而后输出给主控单元FPGA的信号同步采集；信号同步采集单元采集到的信号通过USB接口程序实现数据的显示和存储，主控单元FPGA与显示存储单元相连接；输入控制端的输出信号输入给主控单元FPGA，主控单元FPGA通过解读控制信号来进行系统同步时间和输出波形的设定；GPS同步单元与主控单元双向连接，用于进行同步时间信息的提取和处理。

本发明中采用的信号为：发送频率为0.001Hz~65535Hz之间按编码图案变化，发送波形为单极性或双极性的正弦波或方波。发送频率范围和波形形式可预先设定。

所述的接收机包括阻抗匹配电路、陷波电路、高低通滤波电路、信号放大电路，A/D转换、高速信号采集、大容量数据存储电路、GPS同步单元及人机界面；接收机通过接收电极接收包含地电信息的响应信号，电极送来的信号通过阻抗匹配电路进行阻抗变化，然后通过陷波电路，陷波电路的输入端与阻抗匹配电路的输出端相连，陷波电路的输出端与高低通滤波电路的输入端相连，高低通滤波电路的输出端与信号放大电路的输入端相连，信号放大电路的输出端再将信号传递给A/D转换后进行高速数据采集，数据采集的输出通过基于FPGA的USB接口程序输送给存储单元；接收到信号的传递顺序为：由阻抗匹配电路至陷波电路至高低通滤波电路、信号放大电路至A/D转换采集电路至存储电路。

本发明采用的发送信号为跳频编码序列波形电流，跳频编码序列波形电流的跳频编码载波频率按照某种跳频图案伪随机序列在0.001Hz~65535Hz范围内跳变。跳频图案由伪随机码控制，使载频的跳变具有均匀分布的性质。这样发送电流既有随机波形电流的特征，又能重复产生。跳频编码发送信号功率谱密度高于噪声功率谱密度；频率跳变量使得各子频率相互不重叠，保持了较高的频带利用率；同时跳频编码的发送方式决定其具有能量比较集中的优点。均符合一个最优电磁探测发送序列的要求。

本发明所述的发送序列由编码序列发生器和直接数字频率合成技术，即DDS配合实现。跳频编码序列波形产生过程为：由线性反馈移位寄存器产生控制字，输出的控制字经过控制字调整模块转化为DDS的频率控制字，经过相位累加器相位累加后作为一个查表地址输送给地址调整模块，地址调整模块根据输出波形选择信号，调整地址值对应于波形存储器（ROM）中不同波形的存储区域，最后根据查表的地址值输出所需数字波形。

由编码序列发生器和直接数字频率合成技术配合实现本发明所述的gold序列和其他跳频序列。发送序列的产生由线性反馈移位寄存器来实现，线性反馈移位寄存器是目前FPGA常用的一种实现伪随机序列的方法，这种方法结构简单、易于实现、而且所产生的伪随机序列具有周期长、随机性好的特点。线性反馈移位寄存器的输出序列具有周期性。定义所述的n级线性反馈移位寄存器模块的生成多项式为：

$p\left(x\right)=\underset{i}{\overset{n}{\mathrm{\Σ}}}{g}_i{x}^i=g_n{x}^n+g_{n-1}{x}^{n-1}+g_{n-2}{x}^{n-2}+\·\·\·+g_2{x}^2+g_1{x}^1+g_0---\left(1\right)$

式中，p（x）为关于变量x的多项式；g_i（i=0,1,…，n）取二进制的0或1，而g_i=1表示在产生的序列中被连接参与反馈，g_i=0表示被断开不参与反馈；n为线性反馈移位寄存器的级数；线性反馈移位寄存器产生的周期最长的一种序列，又称最长线性移位寄存器序列，即为m序列。线性反馈移位寄存器的输出序列长度N和线性反馈移位寄存器的级数n的关系为N=2ⁿ-1，n为大于0的整数。将两组m序列的移位寄存器串联或并联后就可得到Gold序列。序列均满足上述对发送波形的需求，发送机按所述跳频编码序列进行发送。

本发明所述跳频编码技术的地面电磁勘探接收方法，是在距离人工激励源数米到数十公里的范围内放置接收机进行数据接收，同时实时同步记录发送的电流信号和接收机所在位置的地电响应信息，同步精度影响到后期数据处理结果的准确性，可采用授时精度达到纳秒级的GPS模块来实现同步。

本发明所述接收机接收到的信息的处理方法为：通过大地系统对发送信号的频谱响应接收相关信息，求解大地系统的冲激响应来反演出相关地球物理参数。跳频编码地面电磁探测方法携带信息的机制与传统跳频有根本的区别。接收单元记录的电磁场响应不仅有大地系统响应，还有接收系统自身的系统响应。收发距为r处记录的大地电磁系统响应u(r,t)，即接收端大地电磁系统响应为接收系统传递函数、大地系统冲激响应及场源编码序列三者的卷积与噪声信号的和。

u(r,t)=y(r,t)+n(r,t)=g_r(t)*g_e(r,t)*f(t)+n(r,t)        （2）

式中，g_r(t)是接收系统的传递函数，g_e(t)是未知的大地系统冲激响应，n(r,t)是不相关噪声，f(t)是场源扩频编码调制后的编码序列。由式（2）可知，如果想要得到大地系统的冲激响应，首先要已知接收系统冲激响应g_r(t)。接收系统冲激响应可通过直接记录发送电流，与同步记录的场源附近电磁场响应卷积得到。

假设所有电磁场接收单元有相同的系统特性，对式（2）求输入与输出信号的互相关，由于噪声信号与激励信号不相关，因而其互相关等于零，随后进行傅立叶变换，求得大地频率响应特性G_e(jω)为：

$G_e\left(\mathrm{j\ω}\right)=\left(\frac{S_{\mathrm{uf}}\left(\mathrm{j\ω}\right)}{S_f\left(\mathrm{j\ω}\right)}\right)/G_r\left(\mathrm{j\ω}\right)---\left(3\right)$

其中S_uf(jω)和S_f(jω)分别是u(r,t)与f(t)的互相关R_uf(τ)和f(t)的自相关R_f(τ)的傅立叶变换，进行幅度和相位分析可估计出地电阻率随频率的变化。通过式（3）进行运算可分离大地冲激响应和观测系统冲激响应，大地冲激响应是地电阻率、收发距及频率或时间的函数。根据大地频率响应特性可得到不同探测深度条件下的地球物理参数信息。

本发明不同于目前普遍应用的以方波为场源的频率域电法勘探，而采用按跳频图案变化的单极性或双极性正弦波和方波进行探测，其频率范围宽、频谱精细化程度高，使得信号包含丰富的地球物理信息，提高地层分辨率和探测效率；而发送序列的复杂度高的特点使得序列有着良好的自相关和互相关特性，可尽可能减小信号与干扰信号的相关性，有效避开环境温度和湿度、外界地电和天电干扰等随时间而变化的因素对观测结果的影响，同时提高了观测效率；跳频编码发送序列发射能量集中的特点，可实现小能量条件下高信噪比探测。在信号检测方面，本发明在信号接收端进行相关接收，利用跳频编码序列进行相关运算，求解包含地球物理信息的大地系统冲激响应，避免了不相关信号码序列的干扰。这样采用多种技术手段来优化传统发送接收方法，为地面电磁探测的精细化、高分辨率测深，及仪器小型化发展提供了新的思路。

附图说明

图1本发明扩频编码电法勘探方法原理框图；

图2本发明发送机原理框图；

图3本发明接收机原理框图；

图4基于FPGA实现的跳频序列波形合成单元原理图；

图5GPS同步系统原理图；

图6跳频编码序列方波信号示意图；

图7跳频编码序列正弦波信号示意图。

具体实施方式

下面结合附图和具体实施方式对本发明做进一步说明。

所述的扩频编码电法勘探新方法原理框图如图1所示。由扩频编码序列发生单元为主要组成的发送单元发送信号，通过大地系统和观测系统的响应，叠加噪声信号后得到记录的大地电磁系统响应。而接收系统冲激响应可通过直接记录发送电流，与同步记录的场源附近电磁场响应卷积得到。通过接收端记录的观测结果对发送的编码序列做相关运算，最终得出大地系统的冲激响应，大地冲激响应是地电阻率、收发距及频率或时间的函数。根据大地频率响应特性可得到不同探测深度条件下的地球物理参数信息。其中发送机单元发送序列和接收机单元接收信号要通过同步跟踪，目的为保持收发同步和采集同步。

应用本发明跳频编码技术的电法勘测系统由发送机和接收机两部分组成，发送信号为按跳频图案变化的单极性或双极性方波及正弦波，编码序列为最大周期为N=2ⁿ-1的m序列通过串联或并联合成的2ⁿ+1个Gold序列和正交或准正交的其他扩频编码序列，发送机和接收机之间采用GPS进行采样同步和收发同步控制。

所述的发送机原理框图如图2所示。所述的发送机包括供电电源、快熔、升压电路、整流滤波、主控单元FPGA，即现场可编程门阵列、隔离驱动单元、保护单元、逆变单元、电压电流检测模块、A/D转换模块、输入控制模块、显示存储模块、GPS同步模块以及隔离电源组。所述的主控单元FPGA包括编码序列频率波形合成单元、数据采集单元、USB接口程序；发送机的连接方式为：供电电源的输出端连接快熔后输入给升压电路，升压电路输出作为整流滤波单元的输入，整流滤波单元的输出电压输入给逆变单元；主控单元FPGA的编码序列频率—波形合成单元输出跳频编码序列，输入给隔离驱动电路，隔离驱动电路的输出为逆变单元开关管的输入驱动信号，主控单元FPGA与控制保护单元相连，保护单元与逆变单元相连接，用于对逆变单元进行保护；逆变单元与电压电流检测单元相连接，电压电流检测的输出作为A/D转换的输入，而后输出给主控单元的信号同步采集；信号同步采集单元采集到的信号通过USB接口程序实现数据的显示和存储，主控单元FPGA与显示存储单元相连接；输入控制端的输出信号输入给主控单元FPGA，主控单元FPGA通过解读控制信号来进行系统同步时间和输出波形的设定；GPS同步单元与主控单元双向连接，用于进行同步时间信息的提取和处理。

供电电源由蓄电池或小型发电机供入、经过升压和整流滤波，提供给发送机。在发送机中，主控单元采用FPGA作为控制芯片，其主要实现的功能有编码序列频率—波形合成的实现、同步信息提取、USB传输接口程序及输入控制和显示存储的实现，编码序列频率—波形合成单元用于合成跳频编码序列图案，输出单极性或双极性的正弦波或方波跳频编码信号，同时可输出特定的单一频率方波或正弦波信号；合成的编码序列信号与隔离驱动电路驱动相连，输出给逆变器单元，逆变单元对波形进行功率放大后通过发送电极供入大地。保护单元对发送机进行过压、欠压、过流保护。电流电压检测检测逆变单元的工作状态，而后与A/D转换相连接，所得的数字信号通过高速数据采集单元传送给主控单元FPGA，再通过USB接口传输给上位机显示及存储。发送机同步单元为全球定位系统GPS来实现，通过主控单元FPGA对GPS模块接收到的世界时信息进行提取，利用世界时秒脉冲进行运算，最终实现同步使能信号的输出，GPS同步使能信号传递给主控单元FPGA来实现发送系统A/D采集及发送的同步控制。

所述的接收机原理如图3所示。所述的接收机包括阻抗匹配电路、陷波电路、高低通滤波电路、信号放大电路，A/D转换、高速信号采集、大容量数据存储电路、GPS同步单元及人机界面；接收机通过接收电极M、N接收包含地电信息的响应信号，电极送来的信号通过阻抗匹配电路进行阻抗变化，然后通过陷波电路，陷波电路的输入端与阻抗匹配电路的输出端相连，陷波电路的输出端与高低通滤波电路的输入端相连，高低通滤波电路的输出端与信号放大电路的输入端相连，信号放大电路的输出端再将信号传递给A/D转换后进行高速数据采集，数据采集的输出通过基于FPGA的USB接口程序输送给存储单元；接收到信号的传递顺序为：由阻抗匹配电路至陷波电路至高低通滤波电路、信号放大电路至A/D转换采集电路至存储电路。阻抗匹配电路用于对大地进行阻抗匹配，由于各个电极具有较高的接地电阻且接地状况各不相同，相应的阻抗变换电路可减弱不同接地电阻所带来的影响；陷波器单元的设计是为了应对50Hz、150Hz谐波工业游散电流的严重干扰；信号滤波及放大单元对接收到的信号进行进一步的调理；A/D转换、高速信号采集及大容量数据存储单元进行存储；FPGA对GPS模块接收到的世界时信息进行提取，利用世界时秒脉冲进行运算，最终实现同步使能信号的输出，GPS同步使能信号传递给FPGA来实现发送系统A/D采集及发送的同步控制。

所述的勘探系统采用跳频编码序列的方波和正弦波信号。在本实施例中采用FPGA作为主控单元。图4所示为基于FPGA实现的跳频序列波形合成原理图，整个跳频编码序列单极性或双极性方波和正弦波及单频方波和正弦波均由FPGA来实现，采用的方法为编码序列发生器结合直接数字频率合成技术。FPGA产生编码序列波形的内部模块有：编码序列发生器、单频控制字发生器、波形模式选择模块、控制字调整模块、相位累加器、地址调整模块、波形存储器。连接方式包括外部连接和内部连接。外部连接为：外部控制端smode连接波形模式选择模块实现对发送波形模式的选择，包括编码序列输出模式和单频输出模式；外部控制端swave实现对发送波形形式的选择，包括方波波形或正弦波波形。内部连接顺序为：编码序列发生器及单频控制字发生器、波形模式选择模块、控制字调整模块、相位累加器、地址调整模块、波形存储器（ROM）。其中，由编码序列发生器或单频控制字发生器产生控制字，经波形模式选择模块后再传递给控制字调整模块，转化为DDS的频率控制字，经过相位累加器相位累加后作为一个查表地址输送给地址调整模块，地址调整模块根据输出波形选择信号调整地址值对应于波形存储器中不容波形的存储区域，最后根据查表的地址值输出所需数字波形。编码序列发生器产生的跳频编码图案由线性反馈移位寄存器产生。

所述勘探系统的GPS同步单元原理图如图5所示。主要由GPS同步模块、FPGA主控芯片来实现。连接方式为：GPS通过天线接收卫星传送的世界时和秒脉冲1PPS信号，GPS模块将接收到的信息传送给FPGA控制器，同时外部对FPGA输入设定时间，通过FPGA内部逻辑运算输出使能信号，使能信号传送给编码序列频率—波形合成单元和数据采集单元，实现收发同步和采集同步。FPGA内部模块包括同步时间设置模块、协调世界时（UTC）提取模块、减法器模块、计数器模块及比较器模块。时间提取模块用于提取GPS模块接收到的协调世界时信息，同步时间设置模块及UTC时间提取模块将信息传送给减法器进行减法运算，运算结果输出给比较器模块，计数器模块通过GPS提取的秒脉冲信号来计数。当1PPS的计数值小于Δt时，比较器模块输出低电平，A/D转换数据接收模块不工作，当1PPS的计数值等于Δt时，比较器模块输出高电平使能信号。

图6、图7所示分别为跳频编码序列方波信号和跳频编码序列正弦波信号。

所述的勘探系统的信号采集存储电路需要高速、高分辨率的A/D采样，建议A/D转换器在16位以上，数据存储需要大容量的存储器，存储长时间段的包括发送信号和接收信号在内的数据。以100kHz数据采集为例，至少包括发送和接收信号在内的两路信号采集10个小时的数据为：16×100000×60×60×10=5.76×10¹⁰（bit）≈6.71（G），所以建议存储器容量≥10G。

所述的接收机接收到的信号输送到计算机进行数据处理，进行发送信号与接收信号之间的相关运算便可去除干扰信号的影响，得到有用的地质探测数据，进而根据常规计算就可以得到大地的视电阻率信息。相关性运算的运算方法如发明内容所述。

Claims (9)

1.一种基于扩频编码技术的地面电磁勘探方法，其特征在于，所述的方法是通过发送电极向地下供入跳频编码序列波形电流信号作为人工激励源，同时同步记录发送的电流信号和接收机所在位置的地电响应信息，经过相关辨识方法，得到大地的频谱响应，计算地球物理信息。

2.根据权利要求1所述的基于扩频编码技术的地面电磁勘探方法，其特征在于，所述跳频编码序列波形电流信号为随跳频图案跳变的编码信号，为单极性或双极性的正弦波或方波信号；所述信号频谱分布无限接近白噪声频谱特性；所述电流波形信号的编码序列频率覆盖整个频率段。

3.根据权利要求1或2所述的基于扩频编码技术的地面电磁勘探方法，其特征在于，所发送的电流波形编码序列信号包括由m序列通过串联或并联合成的2ⁿ+1个Gold序列，其中m序列为由n级线性反馈移位寄存器产生的最长序列；线性反馈移位寄存器的输出序列具有周期性；定义所述的n级线性反馈移位寄存器模块的生成多项式为：

$p\left(x\right)=\underset{i}{\overset{n}{\mathrm{\Σ}}}{g}_i{x}^i=g_n{x}^n+g_{n-1}{x}^{n-1}+g_{n-2}{x}^{n-2}+\·\·\·+g_2{x}^2+g_1{x}^1+g_0$

式中，p（x）为关于变量x的多项式；g_i（i=0,1,…，n）取二进制的0或1，而g_i=1表示在产生的序列中被连接参与反馈，g_i=0表示被断开不参与反馈；n为线性反馈移位寄存器的级数；线性反馈移位寄存器的输出序列长度N和线性反馈移位寄存器的级数n的关系为N=2ⁿ-1，n为大于0的整数。

4.根据权利要求1所述的基于扩频编码技术的地面电磁勘探方法，其特征在于：所述编码序列单极性、双极性的正弦波和方波电流的最低频率为0.001Hz，最高频率为65535Hz，发送频率为0.001Hz~65535Hz之间按编码图案变化，且发送频率范围可预先设定。

5.根据权利要求1所述的基于扩频编码技术的地面电磁勘探方法，其特征在于：在所述的勘探过程中，实时同步记录发送的电流信号和接收机所在位置的地电响应信息，实时同步记录的发送电流信号用于相关辨识算法，以求解大地电阻率信息。

6.根据权利要求1所述的基于扩频编码技术的地面电磁勘探方法，其特征在于：所述勘探系统发送机产生的跳频编码序列电流信号通过发送电极供入地下，接收机在接收端接收经大地响应后携带地电信息的信号，同时对发送的电流信号和接收机所在位置记录的地电响应信息进行GPS同步采集和存储，经过相关辨识方法，对接收信号和发送信号进行相关性运算，收发距为r处记录的地电磁系统响应u(r,t)：u(r,t)=y(r,t)+n(r,t)=g_r(t)*g_e(r,t)*f(t)+n(r,t)

式中，g_r(t)是接收系统的传递函数，g_e(t)是未知的大地系统冲激响应，n(r,t)是不相关噪声，f(t)是场源扩频编码调制后的编码序列；若得到大地系统的冲激响应，首先要已知接收系统冲激响应g_r(t)；接收系统冲激响应通过直接记录发送电流，与同步记录的场源附近电磁场响应卷积得到；假设所有电磁场接收单元有相同的系统特性，求输入与输出信号的互相关，由于噪声信号与发送信号不相关，因而其互相关等于零，随后进行傅立叶变换，求得大地频率响应特性G_e(jω)为：

$G_e\left(\mathrm{j\ω}\right)=\left(\frac{S_{\mathrm{uf}}\left(\mathrm{j\ω}\right)}{S_f\left(\mathrm{j\ω}\right)}\right)/G_r\left(\mathrm{j\ω}\right)$

其中S_uf(jω)和S_f(jω)分别是u(r,t)与f(t)的互相关R_uf(τ)和f(t)的自相关R_f(τ)的傅立叶变换，进行幅度和相位分析估计出地电阻率随频率的变化；通过大地频率响应特性进行运算，分离大地冲激响应和观测系统冲激响应，大地冲激响应是地电阻率、收发距及频率或时间的函数；根据大地频率响应特性得到不同探测深度条件下的地球物理参数信息。

7.一种应用权利要求1所述方法的勘测系统，包括发送机与接收机，其特征在于：所述发送机用于发送单极性或双极性的正弦波或方波跳频编码信号，通过电极供入地下；所述的发送机包括供电电源、快熔、升压电路、整流滤波、主控单元FPGA，即为现场可编程门阵列、隔离驱动电路、保护单元、逆变单元、电压电流检测模块、A/D转换模块、输入控制端、显示存储模块、GPS同步模块以及隔离电源组；所述的现场可编程门阵列FPGA主控单元包括编码序列频率—波形合成单元、数据采集单元、USB接口程序；编码序列频率—波形合成单元用于合成跳频编码序列图案，输出单极性或双极性的正弦波或方波跳频编码信号，同时输出单一频率方波或正弦波信号；由编码序列频率—波形合成单元合成的编码序列信号通过隔离驱动电路驱动逆变器单元对波形进行功率放大后，通过发送电极供入大地；A/D转换及数据采集单元对发送的电流信号进行同步采集和存储；发送及采集同步单元采用全球定位系统GPS，通过输入控制端设置系统启动时间，依靠GPS的世界时和秒脉冲对系统进行收发同步和采样同步；保护单元对发送机进行过压、欠压、过流保护；

所述的接收机通过接收电极来接收反应地电信息的响应信号；所述的接收机主要包括阻抗匹配电路、陷波电路、高低通滤波电路、信号放大电路，A/D转换、高速信号采集、大容量数据存储电路，GPS同步单元及人机界面。

8.根据权利要求7所述的勘测系统，其特征在于，所述的发送机包括供电电源、快熔、升压电路、整流滤波、主控单元FPGA，即现场可编程门阵列、隔离驱动电路、保护单元、逆变单元、电压电流检测模块、A/D转换模块、输入控制端、显示存储模块、GPS同步模块以及隔离电源组；所述的主控单元FPGA包括编码序列频率—波形合成单元、数据采集单元、USB接口程序；发送机的连接方式为：供电电源的输出端连接快熔后输入给升压电路，升压电路输出作为整流滤波单元的输入，整流滤波单元的输出电压输入给逆变单元；主控单元FPGA的编码序列频率—波形合成单元输出跳频编码序列，输入给隔离驱动电路，隔离驱动电路的输出为逆变单元开关管的输入驱动信号，主控单元FPGA与控制保护单元相连，保护单元与逆变单元相连接，用于对逆变单元进行保护；逆变单元与电压电流检测单元相连接，电压电流检测的输出作为A/D转换的输入，而后输出给主控单元FPGA的信号同步采集；信号同步采集单元采集到的信号通过USB接口程序实现数据的显示和存储，主控单元FPGA与显示存储单元相连接；输入控制端的输出信号输入给主控单元FPGA，主控单元FPGA通过解读控制信号来进行系统同步时间和输出波形的设定；GPS同步单元与主控单元双向连接，用于进行同步时间信息的提取和处理。

9.根据权利要求7所述的勘测系统，其特征在于，所述的接收机包括阻抗匹配电路、陷波电路、高低通滤波电路、信号放大电路，A/D转换、高速信号采集、大容量数据存储电路、GPS同步单元及人机界面；接收机通过接收电极接收包含地电信息的响应信号，电极送来的信号通过阻抗匹配电路进行阻抗变化，然后通过陷波电路，陷波电路的输入端与阻抗匹配电路的输出端相连，陷波电路的输出端与高低通滤波电路的输入端相连，高低通滤波电路的输出端与信号放大电路的输入端相连，信号放大电路的输出端再将信号传递给A/D转换后进行高速数据采集，数据采集的输出通过基于FPGA的USB接口程序输送给存储单元；接收到信号的传递顺序为：由阻抗匹配电路至陷波电路至高低通滤波电路、信号放大电路至A/D转换采集电路至存储电路。

Images