专利名称：一种电法勘探装置的制作方法
技术领域：
本实用新型涉及电法勘探领域，具体而言，涉及一种电法勘探装置和系统。
背景技术：
在寻找金属矿产方面电法勘探(特别是激电法)的应用非常广泛，但是由于生产矿山及其附近存在大量的高压输电线和各种用电设施，这将产生极其强烈的电磁干扰信号；地下开采坑道内大量的金属管道、车辆轨道和废石渣都给电法的供电和电场接收造成很大的困难。所有这一切，都影响了电法勘探方法的应用效果，因此电法勘探在解决危机矿山问题时存在严重的噪声干扰问题。生产矿山附近电磁干扰尽管很大，但一般是随机性的干扰，而相关辨识技术是一种可以有效地去除随机噪声干扰的系统辨识方法，在机械、自动化、仪器仪表、认知科学、生物信息学等领域都已经得到了广泛的应用，但是在地球物理探测方面却没有得到足够的重视和应用。目前，电法勘探中一般是通过堆叠、滤波、远参考等方法去除噪声，但是在矿山电法勘探中这些方法用于去除噪声就显得无能为力，因此勘察效果差，直接带来人力财力的巨大浪费。而且，为了提高信噪比，一个常用的方法就是提高发送电源的功率，这必然造成发送机异常笨重，这给山区的勘探工作带来巨大的困难。另一方面，激电法是电法勘探中用于找矿的最有效的手段，特别是复电阻率法，能够提供较多的电性参数，探测精度较高。复电阻率法是指通过测量得到复电阻率谱也就是频率响应的一种电法勘探方法，但是需要测量多个频点的电场响应才能得到一条频率响应曲线，而且低频的时候需要的时间很长，因此探测效率很低，这是造成复电阻率法效果虽好却不能得到广泛运用的原因。在电法勘探中如果采用相关辨识技术，那么一次测量就能得到一条频率响应曲线，可以有效解决复电阻率法的效率低下问题，并且具有复电阻率法探测精度高的优点，因此这是一种新的时间域测量的复电阻率法。图1(a)和(b)示出了现有技术中复电阻率法的观测装置和等效电路。传统的复电阻率法是通过在相当宽的超低频段上观测视复电阻率的幅度谱和相位谱或实分量和虚分量谱，以研究地下介质情况。如图1(a)和(b)所示，发送电流源AB每次发送一个频率的正弦波或者方波，接收端子MN测得对应的电压信号，多次测量得到一个频段内的多个频率的幅度和相位特性，得到复电阻率谱，再利用科尔-科尔模型进行拟合和反演，得到直流电阻率、充电率、时间常数和频率相关系数等电性参数进行地球物理解释。这种方法能提供比较丰富的激电信息，但需在许多频率上做观测才能获得较完整的频谱，所以生产效率低
实用新型内容
[0010]本实用新型的目的是利用相关辨识技术去除随机噪声；在时间域测量一次就可以得到一条频率响应曲线从而提高探测效率；另外，通过采用双科尔-科尔模型进行解释，不仅可以去除电磁耦合效应的影响，还可以得到电磁参数，探测精度提高。根据本实用新型的一个方面，提供一种用于电法勘探的发送机，其特征在于所述发送机包括全球定位系统GPS (Global Positioning System)�？�(1)、伪随机信号发生器(2)、驱动及保护�？�(3)、智能功率�？镮PMdntelligent Power Modules) (5)、温度控制�？棰群凸┑缒？�(7)，和发送电极AB (6)；其中，GPS�？棰磐ü谙蛭彼婊藕欧⑸�(2)传递授时同步信息，完成数据信息的存储，并控制驱动及保护�？�(3)；驱动及保护�？�(3)驱动IPM(5)工作，保护电路用于保护整个系统的安全；温度控制系统(4)监测发射机内部的温度，当超过预定值后复位机器或关闭机器；IPM(5)将输出信号转换成用于电法勘探的电流信号并向发送电极AB(6)输送电流；供电�？�(7)为板卡提供低压电源，并为IPM (5)提供所需高压。优选地，伪随机信号发生器( 进一步根据需要输出不同阶数、不同频率的m序列或者逆重复m序列，配置LCD用于显示发射机信号类型、参数及状态。根据本实用新型的另一方面，提供一种用于电法勘探的接收机，其特征在于所述接收机包括GPS�？�(8)、伪随机信号发生器(9)、接收电极丽(10)、前置放大�？�(11)、滤波及陷波�？�(12)、精密放大�？�(13)、自相关器(14)、互相关器(15)、第一 FFT变换器(16)、第二FFT变换器(17)、频率响应生成器(18)、复电阻率谱生成器(19)、拟合器Q0)和反演器01)。其中，GPS�？�(8)通过串口向伪随机信号发生器(9)传递授时同步信息，实现精准同步并且和发送机的伪随机序列保持精确同相，伪随机信号发生器(9)根据需要输出与发送机相同阶数、相同频率的m序列；接收电极丽(10)接收通过地下介质传来的电位差，然后经过前置放大�？�(11)将接收到的电磁信号进行初步的放大，并经过滤波及陷波�？�(12)滤除高频噪声及系统的工频噪声，再通过精密放大�？�(13)将滤除噪声的信号进行低噪精密放大；自相关器(14)实现伪随机信号的自相关运算；互相关器(15)实现伪随机信号和接收信号的互相关运算。第一 FFT变换器(16)、第二 FFT变换器(17)分别对自相关结果和互相关结果进行时频变换，频率响应生成器(18)运算得到去除观测系统影响后的大地系统频率响应，大地系统频率响应经过复电阻率谱生成器(19)得到大地系统的复电阻率谱，拟合器00)依据双科尔-科尔模型拟合获得视谱参数，反演器进行联合反演获得真谱参数；从而获得地球物理参数信息，并进行实时数据图形显示及存储。根据本实用新型的再一方面，提供一种电法勘探系统，其特征在于包含上述发送机以及接收机。与现有技术相比，根据本实用新型的技术方案能够带来以下的技术效果1、能够去除接收信号中叠加的随机干扰。根据维纳-霍夫方程，用输入的伪随机信号对输出信号进行互相关运算，由于输入信号和输出信号中叠加的随机噪声不相关，所以随机噪声被消除。2、能够去除观测系统对测量结果的影响。对于一个实际系统，电极、导线、接地阻抗及测量装置都会对测量结果产生影响而降低辨识的准确度，利用反卷积运算，这个影响
5很容易被去除。所述反卷积运算可以在频域实现。3、能够得到目标极化体的真谱参数和几何分布。单输入单输出辨识得到的冲激响应和频率响应反映的是发送接收电极之间勘探体积内地下介质的整体效应，采用多极距的高密度几何测深方法，能够得到目标极化体的真谱参数和几何分布，具有水平和纵向分辨力高、反映异常幅度大等优点，特别是在复杂地形条件下地质勘查效果较好。4、能够去除电磁耦合效应的影响。激电法中如何去除电磁耦合效应的影响一直以来都是人们亟待解决的问题。利用双科尔-科尔模型的复电阻率法不仅可以得到激电效应的4个参数，还可以得到关于电磁耦合效应的3个参数，因此不仅可以分离电磁效应对激电参数的影响，利用电磁参数还可以得到更多的信息，因此探测精度大大提高。
图1(a)示出了复电阻率法的观测装置；图1 (b)示出了复电阻率法的等效电路；图2示出了根据本实用新型的电法勘探方法的辨识模型；图3示出了待辨识系统的冲激响应(以一个典型的二阶系统为例)；图4 (a)示出了一个周期的m序列；图4(b)示出了一个周期的逆重复m序列；图4(c)示出了高斯白噪声的示意图；图5示出了利用m序列辨识得到的冲激响应和真实冲激响应对比示意图；图6示出了利用逆重复m序列辨识得到的冲激响应和真实冲激响应对比示意图；图7示出了根据本实用新型的电法勘探方法的辨识模型中考虑观测系统的影响的情形；图8示出了待辨识系统冲激响应反卷积前后的对比示意图；图9示出了多极距观测系统的示意图；图10示出了双科尔-科尔模型幅度谱辨识前后的对比示意图；图11示出了双科尔-科尔模型相位谱辨识前后的对比示意图；图12示出了双科尔-科尔模型拟合的幅度/相位示意图；图13示出了双科尔-科尔模型拟合的实部/虚部示意图；图14示出了根据本实用新型的电法勘探装置中发送机的示意图；图15示出了根据本实用新型的电法勘探装置中接收机的示意图。
具体实施方式
图2示出了根据本实用新型的电法勘探方法的辨识模型。它的基本原理是维纳-霍夫方程。它的核心是采用伪随机序列(包括m序列和逆重复m序列)对接收到的时间序列进行相关处理，得到冲激响应，把冲激响应变到频域就是频率响应。系统输入为v(t)，输出u(t)也就是测量到的信号，h(t)是待辨识系统的冲激响应，r(t)是输入v(t)经过被测系统后的响应，外部的随机噪声干扰η (t)，其中输入ν (t)、输出u(t)是已知的，干扰n(t)为不可测的、未知的。系统输入输出信号ν (t)、u (t)有如下关系[0043]u(t) = r(t)+n(t) = v(t)*h(t)+n(t)(1)式中h(t)为被测系统的时域冲激响应，*为卷积运算。式(1)两边同时用输入信号ν (t)进行互相关运算Rvu (t) = Rvv(t)*h(t)+Rvn(t)(2)式中Rvu(t)表示输入和输出的互相关，Rvv (t)表示输入的自相关，RTO(t)表示干扰n(t)和输入信号的互相关。由于外部的随机噪声干扰和输入信号不相关，因此(2)式中的Rvn(t)为0。可以看到，互相关运算的结果是消去了随机噪声的干扰。因此得到Rvu (t) = Rvv(t)*h(t)(3)这就是著名的维纳-霍夫方程。因此，我们只要求出输入的自相关函数Rvv (t)和输入输出的互相关Rvu(t)，就可以通过解卷积的方法得到被测系统的时域冲激响应h(t)。变换到频域就是频率响应。这样一次发送和接收就得到了一条频率响应曲线，包括幅频响应和相频响应，大大提高了探测效率。但是如果输入是一般信号，那么求解方程C3)是比较困难的。为此我们要寻找一些特殊的信号作为输入信号以简化求解过程。而白噪声信号的自相关函数具有特殊的形式，其自相关函数是一个δ函数。当被辨识系统的输入为白噪声时，只要求出输入输出信号之间的互相关函数，就可以相除直接得到系统的冲激响应。但白噪声只是数学上的一个抽象，工程上不容易实现，更不具有可重复性。工程中常使用伪随机信号作为输入信号。它具有近似白噪声的性质，所以可以保证良好的辨识精度，而且在工程上容易实现。辨识通常采用的伪随机信号是m序列或者逆重复m序列。m序列是指最大长度移位寄存器序列。逆重复m序列也称L序列，是对m序列扩展一倍，然后隔位取反后产生的，由于在一个周期内均值为0，因此在进行系统相关辨识时具有比m序列更优良的特性。附图3是一个待辨识系统的冲激响应h(t)，以一个典型的二阶系统为例。附图4(a)是一个周期，长度为31的m序列，图4(b)是一个周期的逆重复m序列，图4(c)是叠加的噪声n(t)。附图5是利用m序列进行相关辨识的结果，可以看到，经过相关辨识恢复的冲激响应基本和真实的冲激响应重合，噪声的影响变得很小。附图6是利用逆重复m序列进行相关辨识的结果，辨识效果比m序列更好。另外，在本实用新型的一个优选实施例中，为了去除电极、导线、接地阻抗及测量装置的影响，在接收机的数据处理过程中进行了反卷积运算，其原理如下将电极、导线、接地阻抗及测量装置的影响归结为一个系统，它的冲激响应为hs(t)，大地系统的冲激响应记为 ⑴，如附图7所示。那么⑴式变成u(t) = v(t)*he(t)*hs(t)+n(t)(4)(4)式两边用输入伪随机序列进行互相关运算，就可以消去随机噪声的影响，得到Rvu (t) = Rvv(t)*he(t)*hs(t)(5)因此，先求输入和输出信号的互相关rTO(t)以及输入信号的自相关Rvv(t)，反卷积运算后得到观测系统和大地系统总的冲激响应he (t) *hs (t)，然后再通过反卷积运算去除观测系统hs(t)的影响得到大地系统的冲激响应、(t)。从附图8可以看出，去除观测系统的影响后的冲激响应更接近待辨识系统的冲激响应。也可以把卷积和反卷积运算转化成在频域或者复频域相乘或者相除的运算。因此(5)式变为Ρ (ω) = Ρνν(ω) · He (ω) · Η3(ω)(6)这样就得到待辨识系统的频率响应He(O)=(7)
Pvv(CO)-Hs(CO)在发送机附近相关辨识得到的结果就认为是观测系统的冲激响应、(0，因为在发送机附近，可以忽略大地系统的影响，认为只有观测系统的影响。把它变到频域即为Η3(ω)0在地球物理中用到的是复电阻率谱P (ω)。它和频率响应的关系是ρ ( ω ) = K · He ( ω )(8)这里K称为装置系数。如图1 (a)所示，A、B是一对发送电极，M、N是一对接收电极，那么装置系数的大小取决于A、M、N、B四电极的相对位置
Τζ ____(Q)“ J___J___
'am IN IM BM式中AM、an、BN、BM为电极间的距离。得到复电阻率谱以后，利用复电阻率法来进行地球物理解释。通常复电阻率法采用单科尔-科尔模型进行拟合
1
Ι + υωτΥ 这样拟合得到4个激电参数P ρ m、τ、c。由于复电阻率谱中含有电磁耦合效应 ρ(ω) = p0h-m
1-·
(10)
的影响，所以这样得到的激电参数是不够准确的。 而利用双科尔-科尔模型(式(11)所示)进行拟合不仅可以得到4个激电参数P O^m1, C1, T1，还可以得到关于电磁耦合效应的3个参数m2、c2、τ2。
1ρ(ω) = P0 <{ 1 - W1
1
.ι+(Χ
l-m2
1--
\ + {jOT2)
c2
(11)
这7个参数反映的是电极排列勘探体积内的总体信息，称为视谱参数，它们的物理意义是①Ptl-表示零频视电阻率反映电极排列勘探体积内的平均电阻率。②叫-激电效应的视充电率激电效应强度参数(％ )，与电极排列勘探体积内的可极化物质体积含量正相关。③T1-激电效应的视时间常数激电效应特征参数(秒)，与电极排列勘探体积内的可极化物质的粒度大小等结构信息相关。④C1-激电效应的视频率相关系数激电效应过程参数(无量纲)，与电极排列勘探体积内的可极化物质的激电效应类型以及极化物质混合分布均勻性相关。⑤m2-电磁效应的视充电率。⑥τ 2-电磁效应的视时间常数。⑦C2-电磁效应的视频率相关系数。因此该方法不仅可以分离激电谱和电磁耦合谱，利用得到的电磁参数还能得到剩
8余电磁效应(REM)参数和电磁视电阻率P ω等，因此可以提高探测精度。在一个测点得到的是视谱参数，我们需要的是目标极化体(也就是需探明的矿体或者油气藏等)的真科尔-科尔参数和它的几何分布(包括深度、厚度、产状等信息)，所以需要2个以上的测点的信息。采取多极距观测，信息量就极大丰富了。因此在本实用新型的一个优选实施例中，采用了几何测深的方法——多极距观测，可以得到不同深度的异常体的信息。这种工作方式如图9所示，AB是发送电极，M1N1是第一对接收电极，N1M3 (即为M2N2)是第二对接收电极，M3N3是第三对接收电极。同线排列的多个接收器同时接收发送电偶极源的电场响应，收发距按一定的间距增加。由于不同的收发距探测的深度不同，收发距越大反映的深度越大，因此在一个供电点一次供电，就可以接收到不同深度的多个频率的电场响应。沿剖面按一定的点距多次移动观测系统，就可观测到一条剖面上不同收发距的多个频率的径向电场响应，获取丰富的电场信息。可以采取以下两种方法由多极距观测结果得到目标极化体的真谱参数和几何分布。①对多极距的测量结果也就是多个测点的视谱进行联合反演，可以求出目标极化体的几何分布及真谱参数。先对每个测点的复电阻率谱(包括振幅和相位)用常规复电阻率法反演分离出极化体和围岩的视谱，再对分离出的多个不同极距的视相位谱联合反演得到目标极化体的真谱参数和几何分布。②利用分离出的电磁效应的参数可以计算剩余电磁效应(REM)参数和电磁视电阻率2个电磁参数，它们比常规视电阻率能更灵敏地反应地下导电性异常。这样反演每一个测点的复电阻率谱可获得六个参数四个激电参数，两个电磁参数。多极距观测结果用六个参数的拟断面图表示，可以反映地电构造沿剖面和随深度的变化。本领域的技术人员应当理解，在图1、图9中，发送端不一定用电偶极源，也可以使用磁性源线圈等任意一种电法勘探用主动源。接收端也不一定用偶极子接收，可以使用磁棒等任何一种电法勘探用接收器。另外，发送接收不一定采用如图9的多极距观测装置，可以是这种同线排列的，也可以是轴向测量，可以是这种一个发多个收的，也可以一个发一个收，或者多发多收。另外，对得到的频率响应，除了用复电阻率法来解释，还可以用直流电阻率法、激电法、电磁法来解释。在频率接近0时，频率响应曲线在该点的值就是该测点的直流视电阻率，得到多个测点的值后可以用直流电阻率法反演。在高频(例如大于10咕2)，可以求取一些电磁参数，采用电磁法的资料处理和反演方法。在低频段则可以得到一些激电参数，例如视频散率、视频散率道间变化率、道间低频复电阻率的变化率等参数。另外，也可以直接对时间域的冲激响应进行地球物理解释。在考虑复电阻率法时，也不一定是这种乘积形式的双科尔-科尔模型，也可以是和形式的双科尔-科尔模型，还可以是Dias模型、Brown模型或者其他模型。与现有技术相比较，根据本实用新型的技术方案的不同之处至少在于1、采用伪随机序列(包括m序列和逆重复m序列)对电法接收机接收到的时间序列进行互相关运算，得到系统的冲激响应和频率响应；2、采用反卷积去除观测系统的影响；3、采用多极距的观测方式得到目标极化体的真谱参数和几何分布；4、对得到的频率响应采用复电阻率法进行反演解释。采用双科尔-科尔模型可以得到电磁参数，探测更精细。为了对本实用新型技术方案的技术效果有更好的了解，以下通过一个仿真实例进行说明。由于典型二阶系统的特性是我们所熟知的，辨识效果也容易比对，所以下面以对一个二阶系统的辨识结果为例。假设一个待辨识二阶系统冲激响应为h(t)=a · ebt · Sin(Ct)，取a = 11. 547，b = -5, c = 8. 66，如图3所示。该二阶系统截止频率为2Hz，它的冲激响应在Is后基本趋于零，所以系统的过渡时间约为Is。根据参数选择公式(12)和(13)，伪随机序列的码片宽度可以选择0.02s，而序列的周期长度为63即可。
权利要求1.一种用于电法勘探的接收机，其特征在于所述接收机包括GPS�？�(8)、伪随机信号发生器(9)、接收电极MN(IO)、前置放大模块(11)、滤波及陷波�？�(12)、精密放大模块(13)、自相关器(14)、互相关器(15)、第一 FFT变换器(16)、第二 FFT变换器(17)、频率响应生成器(18)、复电阻率谱生成器(19)、拟合器00)和反演器；其中，GPS�？�(8)通过串口向伪随机信号发生器(9)传递授时同步信息，实现精准同步并且和发送机的伪随机序列保持精确同相，伪随机信号发生器(9)根据需要输出与发送机相同阶数、相同频率的m序列；接收电极丽(10)接收通过地下介质传来的电位差，然后经过前置放大模块(11)将接收到的电磁信号进行初步的放大，并经过滤波及陷波�？�(12)滤除高频噪声及系统的工频噪声，再通过精密放大模块(13)将滤除噪声的信号进行低噪精密放大；自相关器(14)实现伪随机信号的自相关运算；互相关器(1 实现伪随机信号和接收信号的互相关运算；第一 FFT变换器(16)、第二 FFT变换器(17)分别对自相关结果和互相关结果进行时频变换，频率响应生成器(18)运算得到去除观测系统影响后的大地系统频率响应，大地系统频率响应经过复电阻率谱生成器(19)得到大地系统的复电阻率谱，拟合器OO)依据双科尔-科尔模型拟合获得视谱参数，反演器进行联合反演获得真谱参数；从而获得地球物理参数信息，并进行实时数据图形显示及存储。
2.一种电法勘探系统，其特征在于通过发送机发送信号以及接收机对接收到的信号进行处理，所述发送机包括全球定位系统GPS (Global Positioning System)模块(1)、伪随机信号发生器(2)、驱动及保护�？�(3)、智能功率�？镮PMdntelligent Power Modules) (5)、温度控制�？�(4)和供电�？�(7)，和发送电极AB (6)；其中，GPS�？�(1)通过串口向伪随机信号发生器(2)传递授时同步信息，完成数据信息的存储，并控制驱动及保护模块(3)；伪随机信号发生器(2)根据需要输出不同阶数、不同时钟频率的m序列或者逆重复m序列等伪随机序列，配置LCD用于显示发射机信号类型、参数及状态。驱动及保护�？镃3)驱动IPM(5)工作，保护电路用于保护整个系统的安全；温度控制系统(4)监测发射机内部的温度，当超过预定值后复位机器或关闭机器；IPM(5)将输出信号转换成用于电法勘探的电流信号并向发送电极AB(6)输送电流；供电�？�(7)为板卡提供低压电源，并为IPM(5)提供所需高压；所述接收机包括GPS�？�(8)、伪随机信号发生器(9)、接收电极MN(IO)、前置放大�？�(11)、滤波及陷波�？�(12)、精密放大�？�(13)、自相关器(14)、互相关器(15)、第一 FFT变换器(16)、第二 FFT变换器(17)、频率响应生成器(18)、复电阻率谱生成器(19)、拟合器(20)和反演器；其中，GPS�？�(8)通过串口向伪随机信号发生器(9)传递授时同步信息，实现精准同步并且和发送机的伪随机序列保持精确同相，伪随机信号发生器(9)根据需要输出与发送机相同阶数、相同频率的m序列；接收电极丽(10)接收通过地下介质传来的电位差，然后经过前置放大�？�(11)将接收到的电磁信号进行初步的放大，并经过滤波及陷波�？�(12)滤除高频噪声及系统的工频噪声，再通过精密放大�？�(13)将滤除噪声的信号进行低噪精密放大；自相关器(14)实现伪随机信号的自相关运算；互相关器(15)实现伪随机信号和接收信号的互相关运算；第一 FFT变换器(16)、第二 FFT变换器(17)分别对自相关结果和互相关结果进行时频变换，频率响应生成器(18)运算得到去除观测系统影响后的大地系统频率响应，大地系统频率响应经过复电阻率谱生成器(19)得到大地系统的复电阻率谱，拟合器00)依据双科尔-科尔模型拟合获得视谱参数，反演器进行联合反演获得真谱参数；从而获得地球物理参数信息，并进行实时数据图形显示及存储。
专利摘要本实用新型涉及一种电法勘探装置，其通过发送机发送信号以及接收机对接收到的信号进行处理。其中发送机包括全球定位系统GPS(GlobalPositioning System)�？�(1)、伪随机信号发生器(2)、驱动及保护�？�(3)、智能功率�？镮PM(Intelligent Power Modules)(5)、温度控制�？�(4)和供电模块(7)，和发送电极AB(6)；接收机包括GPS�？�(8)、伪随机信号发生器(9)、接收电极MN(10)、前置放大�？�(11)、滤波及陷波�？�(12)、精密放大�？�(13)、自相关器(14)、互相关器(15)、第一FFT变换器(16)、第二FFT变换器(17)、频率响应生成器(18)和复电阻率谱生成器(19)、拟合器(20)和反演器(21)。
文档编号G01V3/00GK202330736SQ201120453809
公开日2012年7月11日 申请日期2011年11月16日 优先权日2011年11月16日
发明者李梅, 杜守德, 王海宾, 王超, 范东琦, 荣攀, 邢玉川, 魏文博 申请人:中国地质大学(北京)
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
