抗干扰电法仪试用申请

伪随机信号电法仪（KGR）推广方案国土资源部矿产勘查技术指导中心2015年1月8日目录一、伪随机信号电法仪理论介绍..............................................3二、KGR-1电法仪应用的几个实例.......................................11三、KGR-1抗干扰电法仪的不足及解决方案.......................15四、KGR-1抗干扰电法仪推广方案.......................................16五、申报书内容........................................................................17六、联系方式............................................................................17一、伪随机信号电法仪理论介绍电法勘查是以地球物质电性差异为基础，通过观测和研究地球电流场，推断地下电性分布，进而勘查地下地质结构和寻找地下资源的一类地质勘查方法。大部分电法勘查方法观测和研究的地球电流场，是勘查人员主动供电在地中产生的电流场。这类电法勘查方法称为“主动源电法勘查方法”。1.电法勘查和系统辨识主动源电法勘查方法可概括为图1(a)所示观测系统由发送系统的电流源供给一定波形和强度的电流I，通过发送传感器（不接地回线或者接地电极等）在地下产生一定的电流场；接收系统通过接收传感器（不接地回线或者接地电极等）获得反映地下电流场的感应电动势或电位差ΔU，并由观测仪器（比如，电位仪）记录下来。图1(b)示出了上述观测系统的等效电路其输入端（左端，发送系统）和输出端（右端，接收系统）之间，代表大地系统。目前，主动源电法勘查有两种观测方式时间域观测和频率域观测。时间域观测方式通常用一定占空比的正负方波供电，观测供电期间电流场电位差的建场过程或/和断电后的放电过程ΔU(t)，获得归一化电位差或阻抗的“阶跃响应”Z(t)=ΔU(t)/I0。时间域激电法和瞬变电磁法等属于这类方法。频率域观测方式通常依次用不同频率的连续正负方波供电，观测基波或谐波频率电流场的复电位差ΔU(iω)，获得归一化复电位差或复阻抗的“频谱”Z(i发送系统接收系统电流源电位仪电流源电位仪地球IIΔUI系统ΔU)(tZ发送传感器地面接收传感器（a）(b)图1.主动源电法勘查观测系统（a）及其等效电路（b）ω)=ΔU(iω)/I0。频谱激电法和各种频率测深法(CSAMT及广域电磁法等)属于这类方法。按上述两种观测方式建立的现有电法勘查仪器，都存在抗干扰能力差和野外观测效率低的严重问题。解决这些问题的出路是，利用现代系统科学关于“系统辨识”的新成就，建立全新的电法勘查仪器。图1（b）示出了主动源电法勘查观测系统，按系统科学简化的等效电路。左端发送机通过发送传感器向大地供入的电流强度I(t)，作为被研究的“地球系统”的“输入”，右端接收机通过接收传感器观测到的电位差ΔU(t)，则是“地球系统”的“输出”。大地系统通常可以看成是“线性”和“时不变”的。系统科学理论表明【2】，对于一个线性、时不变系统，任何时间序列输入信号I(t)与输出信号ΔU(t)之间的关系，可以通过如下卷积积分（Wiener-Hopf方程）给出：)()()()()(tItZdsstIsZtU（1）式中，)(tZ为系统的冲激响应，即供电电流I(t)为δ(t)函数（单位脉冲）时，大地的时间响应（输出）；符号*表示卷积，s为积分变量（具时间量纲）。原则上，可以通过输入信号I(t)与输出信号ΔU(t)的反卷积，计算出冲激响应)(tZ，进而通过对)(tZ积分计算得到系统的阶跃响应)(tZ：tdttZtZ0)()(（2）对)(tZ作傅氏变换可得到系统的频谱（传输函数）)(iZ对于tie时间因子，可写出dtetZiZti)()(（3）系统科学中，基于公式（1）由实测系统的输入和输出，确定系统的冲激响应（进而可获得系统的时间响应和频率响应），称为“系统辨识”。所以，从系统科学的观点看，电法勘查就是对地球系统的系统辨识。2.伪随机二进制输入信号虽然，对于任意给定的输入信号I(t)和相应的输出信号ΔU(t)，都可以通过反卷积计算出冲激响应)(tZ；但除计算繁琐外，还经常由于观测误差和干扰导致计算失败。所以，人们总在寻找适当的输入信号波形，以获得高精度的冲激响应)(tZ。在系统辨识的诸多方法中，有一种特殊和十分有效的方法，这就是“伪随机二进制输入信号（PRBS）”的系统辨识方法[1,2,3]。它采用按“最大长度伪随机二进制序列”（m序列）1变化的电流源I(t)供电，供电电流强度I(t)是宽度随机变化，而幅度保持为常数I0的周期性出现的一系列正、负方波脉冲（见图2）。其特点是：（1）电流强度只取值为+I0或－I0。（2）而连续取+I0或－I0的时间长度K*Δt是随机变化的（Δt为常数，称为m序列的“位宽”或“钟脉冲周期”；K为整数，称为“蝉联数”）。（3）经过时间T0后，供电电流波形重复出现。周期T0=N*Δt（N为整数，称为m序列的“位数”或“码数”）。（4）一个周期内，取值+I0的各（正向）脉冲的蝉联数之和，与取值-I0的各（反向）脉冲的蝉联数之和相差不超过1。（5）取12nN，则在一个周期内，蝉联数为n的非零（或者为+I0，或者为-I0）方波脉冲出现，而且只出现一次（n称为m序列的“阶数”）；而1伪随机信号（序列）是指，特征与随机信号的特征充分相似的确定性信号。有各种各样的伪随机信号，其中，最大长度伪随机二进制序列（m序列）是同时具有平衡标准、游程标准和相关性标准的最佳周期自相关函数序列。在n阶GF(p)域序列中，总数（即“位数”N）为pn-1的序列，定义为“最大长度序列”或“m序列”。因此，n阶二进制m序列的位数N=2n-1。所谓“n阶”是指序列中连续出现某一数的最大次数为n。图2.周期为T0，位数N=15的最长伪随机二进制信号序列（m序列）【1】。I(t)+I0-I0蝉联数为K（1≤Kn）的方波脉冲（包括正的和负的方波脉冲）出现的机率为K2/1。上述伪随机二进制信号序列，可以用专用设备（线性反馈移位寄存器LFSR）产生，也可以由数字计算机（运行线性递推式）简单地产生【2,3】。3.伪随机二进制输入信号系统辨识的基本算法前已述及，大地电性系统可视为是线性和时不变系统，在此条件下，供电电流时间序列I(t)和实测电位差时间序列ΔU(t)之间，有形如（1）式的卷积关系。另一方面，在供电电流I(t)是以T0为周期的最大长度伪随机二进制信号序列（m序列）时，I(t)的周期自相关函数（PACF）RI,I10,)11(121,100,)()(1)(0ttNNtNTIIdttItITR（4）图3给出了（4）式的图形表示，它表明m序列的自相关函数)(,IIR是一个三角形脉冲；它具有周期性，且其周期与m序列相同；除在τ=0点附近外，)(,IIR的数值很小。由（4）式可进一步写出近似式：00,)()(nIInTR（5）式中0100)(或sss10)(（6）图3.m序列的周期自相关函数)(,IIR【1】RI,I(τ)τ/Δt为“狄拉克δ函数”或称“脉冲函数”。其一个重要性质是，对于任何一个连续函数f(s)，在包含s=τ的区间Q上的积分Qfdsssf)()()(（7）以及，对f(s)=1的特例，可写出Qdssτ1)(当m序列电流幅值I0=1和位宽Δt«周期T0时，（5）式足够近似成立。进一步，考查供电电流I(t)与实测电位差ΔU(t)的周期互相关函数（PCCF）RI,U：dttUtITRTUI)()(1)(000,（8）将（1）式代入（8）式，并考虑到当t0时，0)()(tItZ，可得dssRsZdsdtstItITsZdtdsstIsZtITRIITTUI)()()()(1)()()()(1)(,0000000,00（9）考虑到（5）式，上式可进一步写成：......)2()()(......])2()()([)()(00000,TZTZZdssTsTssZRUI（10）若选择m序列的周期T0足够大，使在τT0时冲激响应)(Z已衰减殆尽，则（10）式可写成：)()(,ZRUI（11）这说明，适当选择m序列的周期T0和位宽Δt（也就是足够大的位数N），观测和计算供电电流输入I(t)和电位差输出ΔU(t)的互相关函数)(,UIR，就可近似得到大地的冲激响应)(Z。文献【2】导出了近似性更好的表达式（未考虑相邻周期的影响）：StRRZtStttRRZtStRRZIUUIIUUIIUUI/)]()0([2)0(0,)21(21/)]()([)(,/)]()([)(,,,,,,（11.a）式中，000,1TIUdtttItUTtR（11.b）S为自相关函数图形中的三角形面积：tNS)11(（11.c）在获得大地的冲激响应)(Z后，可进一步按（2）和（3）式，算出大地的阶跃响应)(tZ和频谱（传输函数）)(iZ；进而，换算出常规谱激电法和电磁测深法的观测数据由阶跃相应)(tZ可算出常规时间域观测的充、放电过程；而由频谱（传输函数）Z(iω)乘以装置系数K，可得复电阻率法的视复电阻率ρs(iω)：)()(iZKis（12）或者计算各种频率测深方法（可控源音频大地电磁法，广域电磁法或频率测深法等）的视电阻率。可见，基于上述伪随机信号系统辨识原理，可组成全新的主动源电法勘查仪器伪随机信号电法仪器。这种仪器原则上一个周期的观测便可获得常规电法勘查宽频带范围内的丰富信息时间响应和视电阻率频谱。应该指出，上述信息（首先是冲激响应)(tZ）是通过互相关、积分和傅氏变换获得的，这些计算本身都是数字滤波运算，可以压制干扰，能获得很高的信噪比。此外，还可以通过“过采样”（在每一个位宽Δt上分布成百上千个采样点，以这些采样时间记录的数据的平均值，作为该位宽的观测值）和在多个周期上进行观测（相当于“多次叠加”），进一步压制干扰，提高观测精度。这就是说，伪随机信号电法仪器既能高效率获取大地的时间响应和频谱的丰富信息，又有很强的抗干扰能力。4.伪随机信号电法仪的原理框图伪随机信号电法仪的原理框图见图4。由最长二进制伪随机信号（m序列）发生器控制的供电电流I(t)，供入大地系统。其输出电位差ΔU(t)进入相关器，与来自m序列发生器、经时间偏移的供电电流信号I(t-τ)计算互相关函数)(,UIR，按理有)()(,ZRUI。对其分别做时间积分tdttZtZ0)()(和傅氏变换dtetZiZti)()(，最后输出阶跃响应)(tZ和频谱)(iZ。5.结论伪随机信号电法仪是基于系统科学关于“系统辨识”最新成果建立的，不同于现行电法勘查仪器，全新探测理念的主动源电法勘查仪器，其主要优点是：（1）一次观测即可获得大地的时间响应（瞬变响应）和频率响应（频谱），可用于各种主动源电法勘查方法。（2）具有极强的抗干扰能力。（3）发送机和多道接收机分置，不需太大发送功率，观测系统非常轻便。说明：本章的内容由中国地质大学（武汉）罗延钟教授提供参考文献1.徐建华编著，状态估计与系统辨识，【M】北京：科学出版社，1981.2.李白南，伪随机信号及相关辨识，【M】北京：科学出版社，1987.3.Aans-Jutgen