抗干扰电法仪试用申请

1、伪随机信号电法仪（KGR）推广方案国土资源部矿产勘查技术指导中心2015年1月8日目录一、伪随机信号电法仪理论介绍3二、KGR-1电法仪应用的几个实例4三、KGR-1抗干扰电法仪的不足及解决方案4四、KGR-1抗干扰电法仪推广方案4五、申报书内容4六、联系方式4一、伪随机信号电法仪理论介绍电法勘查是以地球物质电性差异为基础，通过观测和研究地球电流场，推断地下电性分布，进而勘查地下地质结构和寻找地下资源的一类地质勘查方法。大部分电法勘查方法观测和研究的地球电流场，是勘查人员主动供电在地中产生的电流场。这类电法勘查方法称为“主动源电法勘查方法”。1. 电法勘查和系统辨识主动源电法勘查方法可概括为图

2、1(a)所示观测系统¾¾由发送系统的电流源供给一定波形和强度的电流I，通过发送传感器（不接地回线或者接地电极等）在地下产生一定的电流场；接收系统通过接收传感器（不接地回线或者接地电极等）获得反映地下电流场的感应电动势或电位差U，并由观测仪器（比如，电位仪）记录下来。图1(b)示出了上述观测系统的等效电路¾¾其输入端（左端，发送系统）和输出端（右端，接收系统）之间，代表大地系统。 发送系统 接收系统电流源 电位仪 电流源 电位仪 地球 I I U I 系统 U 发送传感器 地面 接收传感器（a） (b) 图1. 主动源电法勘查观测系统（a）及其等效电路（b

3、）目前，主动源电法勘查有两种观测方式¾¾时间域观测和频率域观测。时间域观测方式通常用一定占空比的正负方波供电，观测供电期间电流场电位差的建场过程或/和断电后的放电过程U(t)，获得归一化电位差或阻抗的“阶跃响应”Z(t) =U(t)/ I0 。时间域激电法和瞬变电磁法等属于这类方法。频率域观测方式通常依次用不同频率的连续正负方波供电，观测基波或谐波频率电流场的复电位差U(i)，获得归一化复电位差或复阻抗的“频谱” Z(i)=U(i) / I0。频谱激电法和各种频率测深法(CSAMT及广域电磁法等)属于这类方法。按上述两种观测方式建立的现有电法勘查仪器，都存在抗干扰能力差和野

4、外观测效率低的严重问题。解决这些问题的出路是，利用现代系统科学关于“系统辨识”的新成就，建立全新的电法勘查仪器。图1（b）示出了主动源电法勘查观测系统，按系统科学简化的等效电路。左端发送机通过发送传感器向大地供入的电流强度I(t)，作为被研究的“地球系统”的“输入”，右端接收机通过接收传感器观测到的电位差U(t)，则是“地球系统”的“输出”。大地系统通常可以看成是“线性”和“时不变”的。系统科学理论表明【2】，对于一个线性、时不变系统，任何时间序列输入信号I(t)与输出信号U (t)之间的关系，可以通过如下卷积积分（Wiener-Hopf方程）给出： （1）式中，为系统的冲激响应，即供电电流I

5、(t)为(t)函数（单位脉冲）时，大地的时间响应（输出）；符号*表示卷积， s为积分变量（具时间量纲）。原则上，可以通过输入信号I(t)与输出信号U (t)的反卷积，计算出冲激响应，进而通过对积分计算得到系统的阶跃响应： （2）对作傅氏变换可得到系统的频谱（传输函数）¾¾对于时间因子，可写出 （3）系统科学中，基于公式（1）由实测系统的输入和输出，确定系统的冲激响应（进而可获得系统的时间响应和频率响应），称为“系统辨识”。所以，从系统科学的观点看，电法勘查就是对地球系统的系统辨识。2. 伪随机二进制输入信号 图2. 周期为T0，位数N=15的最长伪随机二进制信号序列（m序列

6、）【1】。I(t)+I0-I0 虽然，对于任意给定的输入信号I(t)和相应的输出信号U (t)，都可以通过反卷积计算出冲激响应；但除计算繁琐外，还经常由于观测误差和干扰导致计算失败。所以，人们总在寻找适当的输入信号波形，以获得高精度的冲激响应。在系统辨识的诸多方法中，有一种特殊和十分有效的方法，这就是“伪随机二进制输入信号（PRBS）”的系统辨识方法1,2,3。它采用按“最大长度伪随机二进制序列”（m序列） 伪随机信号（序列）是指，特征与随机信号的特征充分相似的确定性信号。有各种各样的伪随机信号，其中，最大长度伪随机二进制序列（m序列）是同时具有平衡标准、游程标准和相关性标准的最佳周期自相关函

7、数序列。在n阶GF(p)域序列中，总数（即“位数”N）为pn-1的序列，定义为“最大长度序列”或“m序列”。因此，n阶二进制m序列的位数N=2n-1。所谓“n阶”是指序列中连续出现某一数的最大次数为n。变化的电流源I(t)供电，供电电流强度I(t)是宽度随机变化，而幅度保持为常数I0的周期性出现的一系列正、负方波脉冲（见图2）。其特点是：（1）电流强度只取值为+ I0或I0。（2）而连续取+ I0或I0的时间长度K*t是随机变化的（t为常数，称为m序列的“位宽”或“钟脉冲周期”；K为整数，称为“蝉联数”）。（3）经过时间T0后，供电电流波形重复出现。周期T0=N*t（N为整数，称为m序列的“位

8、数”或“码数”）。（4）一个周期内，取值+ I0的各（正向）脉冲的蝉联数之和，与取值- I0的各（反向）脉冲的蝉联数之和相差不超过1。（5）取，则在一个周期内，蝉联数为n的非零（或者为+ I0，或者为-I0）方波脉冲出现，而且只出现一次（n称为m序列的“阶数”）；而蝉联数为K（1K<n）的方波脉冲（包括正的和负的方波脉冲）出现的机率为。上述伪随机二进制信号序列，可以用专用设备（线性反馈移位寄存器LFSR）产生，也可以由数字计算机（运行线性递推式）简单地产生【2,3】。3. 伪随机二进制输入信号系统辨识的基本算法前已述及，大地电性系统可视为是线性和时不变系统，在此条件下，供电电流时间序列I

9、(t)和实测电位差时间序列U(t)之间，有形如（1）式的卷积关系。 另一方面，在供电电流I(t)是以T0为周期的最大长度伪随机二进制信号序列（m序列）时，I(t)的周期自相关函数（PACF）RI,I （4）图3. m序列的周期自相关函数【1】图3给出了（4）式的图形表示，它表明m序列的自相关函数是一个三角形脉冲；它具有周期性，且其周期与m序列相同；除在=0点附近外，的数值很小。由（4）式可进一步写出近似式：RI,I()/t （5）式中 或 （6）为“狄拉克函数”或称“脉冲函数”。其一个重要性质是，对于任何一个连续函数f(s)，在包含s =的区间Q上的积分 （7）以及，对f(s)=1的特例，可写

10、出 当m序列电流幅值I0=1和位宽t « 周期T0时，（5）式足够近似成立。 进一步，考查供电电流I(t)与实测电位差U(t)的周期互相关函数（PCCF）RI,U： （8）将（1）式代入（8）式，并考虑到当t<0时，可得 （9）考虑到（5）式，上式可进一步写成： （10）若选择m序列的周期T0足够大，使在< T0时冲激响应已衰减殆尽，则（10）式可写成： （11）这说明，适当选择m序列的周期T0和位宽t（也就是足够大的位数N），观测和计算供电电流输入I(t)和电位差输出U(t)的互相关函数，就可近似得到大地的冲激响应。文献【2】导出了近似性更好的表达式（未考虑相邻周期的影

11、响）： （11.a）式中， （11.b）S为自相关函数图形中的三角形面积： （11.c）在获得大地的冲激响应后，可进一步按（2）和（3）式，算出大地的阶跃响应和频谱（传输函数）；进而，换算出常规谱激电法和电磁测深法的观测数据¾¾由阶跃相应可算出常规时间域观测的充、放电过程；而由频谱（传输函数）Z (i) 乘以装置系数K，可得复电阻率法的视复电阻率s(i)： （12）或者计算各种频率测深方法（可控源音频大地电磁法，广域电磁法或频率测深法等）的视电阻率。可见，基于上述伪随机信号系统辨识原理，可组成全新的主动源电法勘查仪器¾¾伪随机信号电法仪器。这种仪器原则上

12、一个周期的观测便可获得常规电法勘查宽频带范围内的丰富信息¾¾时间响应和视电阻率频谱。应该指出，上述信息（首先是冲激响应）是通过互相关、积分和傅氏变换获得的，这些计算本身都是数字滤波运算，可以压制干扰，能获得很高的信噪比。此外，还可以通过“过采样”（在每一个位宽t上分布成百上千个采样点，以这些采样时间记录的数据的平均值，作为该位宽的观测值）和在多个周期上进行观测（相当于“多次叠加”），进一步压制干扰，提高观测精度。这就是说，伪随机信号电法仪器既能高效率获取大地的时间响应和频谱的丰富信息，又有很强的抗干扰能力。图4 伪随机信号电法仪的原理框图4. 伪随机信号电法仪的原理框图伪随

13、机信号电法仪的原理框图见图4。由最长二进制伪随机信号（m序列）发生器控制的供电电流I(t)，供入大地系统。其输出电位差U(t)进入相关器，与来自m序列发生器、经时间偏移的供电电流信号I(t -)计算互相关函数，按理有。对其分别做时间积分和傅氏变换，最后输出阶跃响应和频谱。5. 结论伪随机信号电法仪是基于系统科学关于“系统辨识”最新成果建立的，不同于现行电法勘查仪器，全新探测理念的主动源电法勘查仪器，其主要优点是：（1）一次观测即可获得大地的时间响应（瞬变响应）和频率响应（频谱），可用于各种主动源电法勘查方法。（2）具有极强的抗干扰能力。（3）发送机和多道接收机分置，不需太大发送功率，观测系统非

14、常轻便。 说明：本章的内容由中国地质大学（武汉）罗延钟教授提供参考文献1. 徐建华编著，状态估计与系统辨识，【M】北京：科学出版社，1981.2. 李白南，伪随机信号及相关辨识，【M】北京：科学出版社，1987.3. Aans-Jutgen Zepernick and Adolf Filger, 伪随机信号处理¾¾理论与应用，【M】北京：电子工业出版社，2007. 4. Duncan P.M. et al, The development and applications of a wide band electromagnetic sounding system usin

15、g a pseudo-noisesource, Geophysics, 45 (8) 1276-1296, 1980.5. 赵璧如, 赵 健, 张洪魁等，PS100 型I P 到端可控源高精度大地电测仪系统CDMA 技术首次在地电阻率测量中的应用，地球物理学进展，21(2)675-682,2006.6. Bruce Hobbs, Anton Ziolkowski , David Wright, Multi-Transient Electromagnetics (MTEM) controlled source equipment for subsurface resistivity invest

16、igation, 18th IAGA WG 1.2 Workshop on Electromagnetic Induction in the Earth 4/4 El Vendrell, Spain, September 17-23, 2006.二、KGR-1电法仪应用的几个实例1.水槽实验 设计如图所示的水槽，在水中放一直立铜板，利用KGR-1A仪器和多种常规电法仪进行中间梯度测量，将得到的结果进行对比。图5 水槽实验示意图图6 KGR-1和FX-1电法仪干扰前后相位对比曲线图7 KGR-1和SQ-3C电法仪干扰前后Fs对比曲线图8 KGR-1和DZD-6电法仪干扰前后Fs对比曲线水槽试验结

17、果表明：（1）在无干扰环境下，KGR-1与常规仪器在视电阻率、视频散率、视相位三个参数测量结果上吻合较好，精度指标相近，验证了仪器的正确性。（2）在加干扰后，所测三个参数的均方差，KGR-1明显优于常规仪器，即抗干扰能力更强。 2. 延庆石槽铜矿实验北京市延庆县石槽村存在一已知铜矿体，利用KGR-1A仪器与重庆仪器厂的DJS-8大功率激电仪在同一条剖面上进行中间梯度测量，将测得的视电阻率和极化参数进行对比验证。图9 KGR-1和DJS-8视电阻率对比曲线图9 KGR时间常数与DJS-8视极化率对比曲线通过上述对比结果可以看出，视电阻率曲线重合，极化率与时间常数曲线形态一致。并且，由已知的地质和

18、物探资料可知，该测线的40号点左右存在矿体，两种仪器在矿体上都有明显反应。3.山西某矽卡岩型银铜矿矿区位于军用雷达站旁边，电磁干扰较强，常规激电仪无法获得稳定数据。通过KGR-1A测深剖面工作，获得高精度稳定的电法数据，下图为电阻率反演图及时间常数断面图，与已知接触面及矿体有较好的对应关系。图10 上图为电阻率反演图，下图为时间常数断面图4.江西某矽卡型铜矿在江西某矿区利用KGR-1仪器与GDD仪器在同一条剖面上进行激电剖面试验，激电剖面与测深反演断面对比结果显示，KGR-1与GDD的视电阻率和极化参数的对应关系良好（GDD是32个测点，KGR是30个测点，右侧对齐）。且对已知矿体有明显的反演

19、（见矿孔位置在540点左右，斜孔，见矿深部150米左右）。图11 电阻率反演断面对比图图12 充电率（时间常数）反演断面对比图三、KGR-1抗干扰电法仪的不足及解决方案该仪器的较强抗干扰能力，通俗讲属于效率换精度，勘查区干扰水平越大，发射信号的码长就需要设置的足够长，增加了单点数据的采集时间，大大降低了工作效率。因此国内最常用的对称四极测深装置，应用该种仪器，导致加长了工作时间，增加野外开支。解决办法：用多台套接收机的办法换效率，因此，测深工作必须采用三极测深装置（国外普遍采用的装置）。比如：我们在江西某铜矿区，采用三极装置，利用了15台接收机（共30个物理点），码长为6的情况下，2天即可完成。效率大大高于对称四极，完成同样工作量的情况下，节省了时间和经费。只是，投入仪器较多，增加了仪器消耗成本！四、KGR-1抗干扰电法仪推广方案1.目标通过在整装勘查区或矿集区组织KGR-1新技术的示范推广，进一步挖掘仪器的不足并不断完善，加快推进KGR-1抗干扰电法仪的实用化进程，促进科技成果转化为生产力，提升找矿突破战略行动科技引领作用。2. 申请单位的条件及义务（1）资质要求：承担整装勘查区或者矿集区2015年度IP工作；具备相应地球物理资质，或项目负责人具有丰富的野外电法施工经验。（2）遵循自愿原则，各勘查单位按照附件要求填写申请书并报送部矿产勘查技术指导中心，由技术指导中心