伪随机信号电法仪的抗干扰参数_罗延钟.pdf

2023 年第 38 卷 第1期2023，38(1):0474-0483地球物理学进展Progress in Geophysicshttp:/wwwprogeophyscnISSN 1004-2903CN 11-2982/P罗延钟，陆占国，王寒冰，等2023 伪随机信号电法仪的抗干扰参数 地球物理学进展，38(1):0474-0483，doi:106038/pg2023FF0063LUO YanZhong，LU ZhanGuo，WANG HanBing，et al 2023 Anti-interference parameters of instrument for electrical prospecting using pseudo randomsignal Progress in Geophysics(in Chinese)，38(1):0474-0483，doi:106038/pg2023FF0063伪随机信号电法仪的抗干扰参数Anti-interference parameters of instrument for electrical prospecting using pseudorandom signal罗延钟1，陆占国2，王寒冰2，孙国良3，高隆钦2LUO YanZhong1，LU ZhanGuo2，WANG HanBing2，SUN GuoLiang3，GAO LongQin2收稿日期2021-05-12;修回日期2022-06-27投稿网址http:/www progeophys cn第一作者简介罗延钟，男，1935 年生，1957 年毕业于北京地质学院，教授，长期从事电法勘查教学和科学研究 E-mail:yluo cugb edu cn1 中国地质大学，武汉4300742 北京桔灯地球物理勘探有限公司，北京1022003 北京航空航天大学，北京1000831 China University of Geosciences，Wuhan 430074，China2 Beijing Orangelamp Geophysical Exploration Company Limited，Beijing 102200，China3 Beijing University of Aeronautics and Astronautics，Beijing 100083，China摘要基于伪随机信号系统辨识的电法勘查仪器多试图借助于相关处理 计算大地(系统)输入信号(供电电流 I(t)和输出信号(电位差 U(t)的互相关 I，U(t)和自相关 I，I(t)以期提高获取大地阻抗的瞬变过程(冲激响应和阶跃响应)和/或频谱(频率响应)的抗干扰性能 笔者先前的研究结果表明，仅在偏移时间 t 小于位宽 t，特别是在零偏移时间，互相关I，U(t)和自相关 I，I(t)函数值较大时才有较强的抗干扰能力;而时间偏移 t 大于位宽 t 时，相关函数值非常小，易受干扰而畸变 利用全周期相关函数值计算大地(系统)时间响应和频率响应的伪随机相关辨识的抗干扰性，并不理想 在电磁干扰较大或信噪比较低(近于或小于 0 dB)时，不能取得可靠的时间响应和频率响应数据 基于上述情况，本文提出仅利用零偏移时间的互相关和自相关函数值 I，U(0)和 I，I(0)，获取具有高抗干扰性能的导电和激电参数 文中给出了这些抗干扰参数的具体算法，并列举仿真模拟和野外观测数据，证明它们的实际抗干扰能力关键词电法勘查;观测仪器;伪随机信号;抗干扰参数中图分类号P631文献标识码Adoi:10 6038/pg2023FF0063AbstractThe instruments for electrical and electromagneticprospecting using pseudo random signal have tried toImprove anti-interferenceabilityforgetingtheearthfrequency response or/and time response by means ofcorrelation processing to calculate the auto-correlationII(t)and cross-correlation IU(t)of the input signal(strength of supply current I(t)and output signal(potential difference U(t)Previous research by theauthers suggests that just in the condition of the offset timeis litter than the bit width of the pseudo random signal，especially at the zero offset time，the function value ofauto-correlation I，I(t)and cross-correlation I，U(t)arelarge，it has strong anti-interference ability;and when theoffset time great than the bit width，the function value ofauto-correlation I，I(t)and cross-correlation I，U(t)arevery small，susceptible to interference and distortion Theanti-jamming capability of the“pseudo random correlationidentification”，which using the correlation function valueof the whole cycle，to get the earth frequency response or/and time response，is clearly not ideal In the condition ofhigh electromagnetic interference or low SN，cant getreliable data of earth frequency response or/and timeresponse Based on above-mentioned condition，this paperproposes to use only the function value of auto-correlationI，I(0)and cross-correlation I，U(0)at the zero offsettime，to obtain conductive and IP parameters with highanti-interference performance The concrete algorithm ofthese parameters and some simulation and field observationdata，which proved actual anti-jamming ability of theseparameters are given in this paperKeywordsElectrical prospecting;Instrument;Pseudorandom signal;Anti-interference parameter2023，38(1)罗延钟，等:伪随机信号电法仪的抗干扰参数(www progeophys cn)0引言基于伪随机信号系统辨识的电法勘查仪器多试图借助于相关处理 计算大地(系统)输入信号(供电电流 I(t)和输出信号(电位差 U(t)的互相关 I，U(t)和自相关 I，I(t)以期提高获取大地阻抗的瞬变过程(冲激响应和阶跃响应)和/或频谱(频率响应)，即“相关辨识”的抗干扰性 能(Duncan et al，1980;徐建华，1981;李白南，1987;Chu，1990;Wright et al，2002;Zepernick and Filger，2005;Hobbs et al，2006;Ziolkowski et al，2007;李梅等，2008;汤井田和罗维斌，2008;刘景贤，2010;刘金琨等，2013;淳少恒等，2014;张宇等，2014;罗延钟等，2015，2016)近期笔者(罗延钟等，2020)通过仿真模拟和野外试验，初步探索了“相关辨识”的抗干扰性能 得到的结果是，无论是否存在干扰(包括理论的随机干扰和实测的电磁干扰)，相关谱分析算法计算的阻抗频谱 Z2(i)和直接谱分析算法计算的阻抗频谱 ZI，U(i)，都几乎完全重合 这表明，相关处理并未明显增加伪随机信号观测数据计算大地阻抗频谱的抗干扰能力 究其原因是，仅在偏移时间t 小于伪随机信号位宽 t，特别是在零偏移时间，互相关 I，U(t)和自相关 I，I(t)函数值较大时才有较强的抗干扰能力;而时间偏移 t 大于位宽 t 时，相关函数值非常小，易受干扰而畸变;致使利用全周期相关函数值计算大地(系统)时间响应和频率响应的伪随机相关辨识的抗干扰性，并不理想 在电磁干扰较大或信噪比较低(近于或小于 0 dB)时，不能取得可靠的时间响应和频率响应数据 本文将在上述研究成果的基础上，提出仅利用零偏移时间的互相关和自相关函数值 I，U(0)和 I，I(0)，获取具有高抗干扰性能的导电和激电参数;并用仿真模拟和野外试验结果，验证这些参数的实际抗干扰效果1伪随机信号电法仪观测的“视电阻”受前文(罗延钟等，2020)“零偏移时间的互相关和自相关函数值 I，U(0)和 I，I(0)具有高抗干扰性能”的启发，我们计算两者的比值 该比值具有电阻量纲，我们定义为“视电阻”s，用数学公式表示为:s=I，U(0)I，I(0)(1)为了检验“视电阻”s的抗干扰性能，我们仍旧利用先前(罗延钟等，2016，2020)使用的电阻和电容电路(图 1)作仿真模拟研究 图中，发送机向电路供入按 m 序列波形变化的供电电流 Im，其在标准电阻 0上的电位差用接收机 1 记录下来作为电路的输入信号 I(t);供电电流流过电阻 和串联的两组并联电阻 1与电容 C1和电阻 2与电容 C2，其上的电位差 U(t)由接收机 2 记录下来作为电路的输出 各电阻和电容取值为 =350，1=50，2=100，C1=0 1 F 和 C2=10 F 对应的并联电阻 1与电容 C1和电阻 2与电容 C2的时间常数分别为1=1C1=5 s 和 2=2C2=1 ms，用以模拟低频段的激电效应和高频段的电磁效应 仿真模拟可选用不同的 m 序列参数(阶数 n 和位宽 t)，分别模拟低频段，中频段和高频段的响应 本文重点选用较大的位宽 t 以检验低频段探测激电效应的抗干扰性图 1用作计算机仿真模拟的模型 C 电路图=350，1=50，2=100，C1=01 F，C2=10 FFig 1The circuit diagram for computer simulation=350，1=50，2=100，C1=01 F，C2=10 F图 2对图 1 所示 C 电路，采用 7 阶和位宽为 2 s 的m 序列激发计算的阻抗振幅频谱曲线Fig 2The impedance amplitude spectrum curvecalculated for the C circuit shown in Fig 1574地球物理学进展www progeophys cn2023，38(1)图2 示出了对图1 所示 C 电路，采用7 阶和位宽为 2 s 的 m 序列，计算的阻抗振幅频谱理论曲线最高频率只有 10 Hz，频谱曲线主要反映激电效应在频率 0.0110 Hz 频段上，阻抗振幅随频率增高，大约由 500 下降到 450，对应的“极化率”大约为(500 450)/500=10%表 1 示出了采用不同阶数和位宽的 m 序列，在添加和不添加干扰条件下，计算的零时间偏移互相关 I，U(0)，零时间偏移自相关 I，I(0)和进而计算的视电阻 s=I，U(0)/I，I(0)对比图 2 和表 1 可以看出:(1)视电阻值 s为某一频率 fs的大地(系统)阻抗振幅值(2)视电阻值 s与所采用 m 序列的阶数 n 和位宽 t 有关:阶数 n 以及特别是位宽 t越大，对应的频率 fs越低，即视电阻 s接近于大地的低频阻抗值;反之，阶数 n 以及特别是位宽 t 越小，对应的频率 fs越高，即视电阻 s接近于大地的高频阻抗值(3)若在一个测区采用统一的阶数 n和位宽 t，则用一定电极装置测得的视电阻 s可反映测区大地的导电性;将视电阻 s乘以装置系数K，可计算电阻率法的常用参数视电阻率:s=K s(2)2伪随机信号电法仪观测的激电参数若在同一测点上，采用一定电极装置对两个相同阶数、不同大小位宽的 m 序列作观测，则可获得两个分别代表较高频率 fGs和较低频率 fDs的视电阻Gs和 Ds;由此，可进而计算反映激电效应的参数 视频散率(百分频率效应 PFE):Ps=Ds GsGs100%(3)表 1对不同阶数和位宽计算的视电阻Table 1The apprarent resistance calculated for different parameters n and t阶数 n位宽 t/s零偏移互相关 I，U/VA零偏移自相关 I，I/A2视电阻 s=I，U/I，I/是否包含干扰?5104771831477183无干扰510477E+021477189随机干扰 10%5104629630998981463435实测干扰 10%5104791160998981479605实测干扰 10%90.02451 560999996451562随机干扰 10%110.024499761449976无干扰110.02450 23145023随机干扰 10%110.024502811450281实测干扰 10%110.02449.0211449.021实测干扰 10%130.02450.0171450.017随机干扰 10%150.024499521449952随机干扰 10%110.000163565511356551无干扰110.000163565361356536随机干扰 1%110.000163564091356409随机干扰 10%110.000163551371355137随机干扰 100%110.000023499511349951无干扰6742023，38(1)罗延钟，等:伪随机信号电法仪的抗干扰参数(www progeophys cn)表 2 给出了包含不同随机干扰时，双 m 序列视电阻和视频散率的仿真模拟结果 可以看出，即使随机干扰方差达到 100%时，视电阻和视频散率所受的影响仍很小;尤其是在两个 m 序列的位宽相差较大时，几乎看不出随机干扰对视电阻和视频散率的影响表 2包含不同随机干扰时，双 m 序列视电阻和视频散率的仿真模拟结果Table 2Apparent resistances and apparent frequencydispersivities calculated by the computer simulation inthe case of different random disturbances(for dual m sequence pattern)模型:C 电路，极化率 m=10%，时间常数 =05 sm 序列:阶数 n=3，第一位宽 t1=10 s，第二位宽 t2=02 s随机干扰方差/%第一视电阻 Ds第二视电阻 Gs视频散率 PS/%0.00 515135E+030472257E+038 685100 515185E+030 472298E+0386861000 515644E+030472674E+038 696m 序列:阶数 n=7，第一位宽 t1=2 s，第二位宽 t2=02 s随机干扰方差/%第一视电阻 Ds第二视电阻 Gs视频散率 PS/%0.00 482958E+030453087E+036 382100 482206E+030 452253E+0364111000 475455E+030444751E+036 673同样的原理，只用一个 m 序列作观测，也可将观测数据划分为相对较高频率 fGs和较低频率 fDs的两个部分，并分别计算其对应的视电阻 Gs和 Ds，进而由式(3)计算视频散率 PS 图 3 给出了用单个 m序列作观测，计算视频散率的图示 将供电电流和电位差每个脉冲在时间上均分为两个部分，然后分别计算每个脉冲前段和后段的电流和电位差的零时间偏移自相关 I，I(0)和互相关 I，U(0)，进而计算视电阻 前段计算的为较高频视电阻 Gs，后段计算的为较低频视电阻 Ds，由此可由式(3)计算视频散率Ps=Ds GsGs100%表 3 列出了对图 1 所示 C 电路模型作仿真模拟，计算用单 m 序列供电，极化率 m=10%，时间常数 =0 5 s，阶数 n=3，位宽 t=2 s，观测和计算抗图 3单个 m 序列作观测，计算视频散率的图示Fig 3A graphical representation for calculating the apparentfrequency dispersivities(for single m sequence pattern)干扰参数的结果 计算数据表明，在加入方差 20%的随机干扰情况下，抗干扰参数视电阻和视频散率几乎不受影响表 3单个 m 序列的抗干扰参数仿真模拟结果Table 3Apparent resistances and apparent frequencydispersivities calculated by the computer simulation inthe case of different random disturbances(for single m sequence pattern)噪声/%视电阻 s/视频散率 PS/%04495364930651044960649299820449702492441为初步验证视电阻和视频散率的抗干扰能力，我们在北京市郊的空闲场地进行了现场观测试验进行试验的场地存在较强的工业游散电流干扰 分别供电和观测一个和两个 m 序列周期，第一天作“基本观测”，第二天作“检查观测”观测装置:对称四极装置;m 序列参数:阶数 n=7，位宽 t=10 s，观测结果示于表 4，可以看到，无论是作一个 m 序列周期，还是两个周期观测，基本观测和检查观测的视电阻和视频散率，都非常接近，相差小于千分之五(0 5%)这都初步验证了这两个参数的高抗干扰能力774地球物理学进展www progeophys cn2023，38(1)表 4单个 m 序列观测的抗干扰参数场地试验结果Table 4Apparent resistances and apparent frequencydispersivities of ground test results(for singlem sequence pattern)观测周期基本观测视电阻 s/检查观测视电阻 s/基本观测视频散率 PS/%检查观测视频散率 PS/%1326.058325 3882 464242 457522325422325 5482 468332 44616图 4河北某老矿山电法野外采集的原始数据(黑色)和智能去噪后数据(白色)Fig 4The field-observed original data(black)and intelligent de-noised data(bright color)at an old mine in Hebei Province3野外试验结果为检验上述两个抗干扰参数实际应用效果，我们应用研制的伪随机信号激电仪 Abollo IP 在几个金属矿区开展了野外试验 Abollo IP 是我们研制的第一款伪随机信号电法仪，观测参数是视电阻率 s(m)和视频散率 Ps(%)下面报告三个矿区的试验结果(1)河北某老矿山野外试验效果(王寒冰等，2019)该矿山 2010 年曾利用传统激电仪开展激电测深工作 由于严重电磁干扰(见图 4)，观测数据凌乱，不可靠 反演的充电率断面图(图 5a)指示极化体深部向西倾斜 据此施工的钻孔 ZK002 未见矿;而东侧施工的探槽 TC003 和平洞 CM001 打到了矿体(参见图 5d)2018 年 5 月应矿山勘探队邀请，采用 Abollo-IP伪随机信号激电仪再次开展激电测深工作 因为Abollo-IP 仪器观测灵敏度高，对弱信号的识别能力很强，所以选用“中梯测深”装置测量 7 天内，高效地完成了 3 条线、48 个测深点的野外观测工作 图5b 和图 5c 示出了与 2010 年同一条剖面上，对采用Abollo-IP 观测的抗干扰参数视电阻率和视频散率，作反演获得的电阻率和频散率断面图 其上的低电阻率和高频散率异常与探矿工程探明的矿体(图5d)吻合的很好 这很好的验证了伪随机信号电法仪(激电仪)的新参数 视电阻率和视频散率的高抗干扰性能(2)内蒙新巴尔虎右旗某铅锌银矿试验效果该矿处于勘探阶段，电磁干扰并不很强 为检验Abollo-IP 的抗干扰能力，我们采用不同供电电流强度进行试验观测，考查不同信噪比情况下，观测数据的可靠性 为对比，我们在完全相同条件下，采用当前国内最受欢迎的传统激电仪 从美国进口的GDD 时间域激电仪与 Abollo-IP 激电仪同时作观测，对比两种仪器的观测结果 GDD 时间域激电仪采用后延触发，观测视电阻率(m)和视充电率(mV/V，即);而 Abollo-IP 激电仪对应的观测参数是视电阻率(m)和视频散率(%)试验采用中间梯度装置，供电电极距 AB=840 m，测量电极距 MN=测点距=20 m，共 30 个测点，分布于 AB 中间 120 720 m我们选用该矿第 88 号勘探线作为试验剖面 图6 示出了该勘探线的地质剖面图，图中的深绿色线条指示为钻探控制的铅锌银矿脉图 7a 示出了查干布拉根铅锌银矿 88 线分别由Abollo-IP 和 GDD 激电仪在不同供电电流强度条件8742023，38(1)罗延钟，等:伪随机信号电法仪的抗干扰参数(www progeophys cn)下，实测的视电阻率剖面曲线图 可以看出，Abollo-IP 在不同供电电流观测的视电阻率曲线都非常接近 相比之下，不同供电电流 GDD 观测的视电阻率曲线虽然形状基本相同，但数值有明显差异;而且，在最小供电电流(I=7 mA)时，无法获得观测值相对来说，图 7b 所示 Abollo-IP 和 GDD 观测的激电参数，精确度不如图 7a 所示视电阻率的精度高 在大电流时(I=1 A)，两种仪器都能测得形状接近的圆滑激电异常;当供电电流下降到 I=70 mA时，两种仪器测得的激电参数剖面曲线都围绕大电流的圆滑异常有一定程度的上下“蹦跳”;在更小电流时(I=20 mA)，GDD 无法取得激电参数读数，而Abollo-IP 仍能测得视频散率激电异常，虽然在个别测点上围绕大电流圆滑异常的“蹦跳”较大查干布拉根铅锌银矿 88 线的野外试验结果表明，相对于当今最受欢迎的传统激电仪 GDD 而言，Abollo-IP 伪随机信号激电仪，在大供电电流(高信噪比或弱干扰)条件下，能测得精度更高的视电阻率和观测精度相当的激电参数;而在较小供电电流(低信噪比或强干扰)，GDD 不能观测到视电阻率和974地球物理学进展www progeophys cn2023，38(1)图 5河北某老矿山激电测深反演断面图(a)传统激电仪数据反演的充电率;(b)Abollo IP 数据反演的电阻率;(c)Abollo IP 数据反演的频散率;(d)探矿工程资料绘制的地质剖面图Fig 5IP sounding inversion results at an old mine in Hebei Province(a)Charging rate retrieved from traditional IP instrument data;(b)esistivity inversion of Apollo IP data;(c)Dispersion rate of Apollo IPdata inversion;(d)Geological profile drawn from prospecting engineering data图 6查干布拉根铅锌银矿 88 线地质剖面图Fig 6Geologic profile of Chaganbulageng lead-zinc mine 88 line0842023，38(1)罗延钟，等:伪随机信号电法仪的抗干扰参数(www progeophys cn)图 7查干布拉根铅锌银矿 88 线激电中梯试验剖面曲线(a)视电阻率试验曲线;(b)视激电参数试验曲线Fig 7Profile curves of IP intermediate gradient on the Chaganbulageng lead-zinc mine 88 line(a)Apparent resistivity test curve;(b)Apparent IP parameter test curve图 8安徽铜陵冬瓜山铜矿实测电磁噪声 时间特性(a)小比例尺时间轴;(b)中比例尺时间轴;(c)大比例尺时间轴Fig 8Electromagnetic noise observed at Dongguashan Copper Mine in Anhui Province(a)Small scale timeline;(b)Medium scale timeline;(c)Large scale timeline激电参数的条件下，Abollo-IP 仍能获得有效观测值这再次验证了伪随机信号电法仪(激电仪)的新参数 视电阻率和视频散率的高抗干扰性能(3)安徽铜陵冬瓜山铜矿试验效果冬瓜山铜矿是一个老矿山，地下正在采矿，地面为矿山职工及家属生活区 电磁干扰非常严重 以往各种电法勘查方法在此试验，都因电磁干扰无法获得可信的观测数据 我们在此采用 Abollo-IP 伪随机信号激电仪作激电测深，检验新参数的抗干扰能力为了解冬瓜山铜矿区的电磁干扰情况，我们在测深点上不向地下供电的情况下，采用 80 m 测量电极距(MN=80 m)观测纯电磁干扰信号 图8 示出了实测的电磁噪声时间特性 其中，图 8a 为小比例尺时间轴的图形，表现为随时间急剧震荡，幅度达到 100mV/80 m 适当增大时间轴比例尺后，由图 8b 可看出，电磁噪声的震荡周期大致为 20 ms 及其整倍数，即频率大致为 50 Hz 及其高次谐波，这表明电磁噪声主要是工业游散电流干扰 利用简单的数据处理剔除工频干扰后，采用更大时间轴比例尺绘制的电磁噪声时间特性图形示于图 8c 可以看出，剔除工频干扰后剩余的干扰幅度大为降低，大体包含两个部分，一部分是幅度很小(不超过 1 mV/80 m)的随机干扰，另一部分为幅度较大(n mV/80 m)的跃变干扰 跃变干扰应该是地下矿井或坑道中运行的直流电器(特别是电气机车)开启、关断或变速所引起的184地球物理学进展www progeophys cn2023，38(1)激电测深采用对称四极装置，测量电极距固定为 MN=80 m，供电电极距 AB 取值为 n 160 m，n=1，2，321 显然，供电电极距 AB 随着 n 增大而增大，电位差观测信号(因而信噪比)随之降低 我们采用在不同观测日期进行“独立检查观测”，以评判观测数据的可靠性表 5 列出了冬瓜山第 400 号测深点，在较小供电电极距 AB=1600 m 和较大电极距 AB=3360 m时，基本观测和两次检查观测测得的视电阻 s和视频散率 Ps 可以看出，较小电极距(AB=1600 m)的信噪比较高，基本观测和检查观测值非常接近，即观测精度很高;而较大电极距(AB=3360 m)的信噪比较低，基本观测和检查观测值相差很大，即观测精度很低 这些观测结果表明，对于冬瓜山铜矿这类电磁干扰，在大电极距、低信噪比条件下，本文提出的“抗干扰参数”也不能获得可靠的观测结果表 5冬瓜山测深点 400 基本观测和检查观测结果Table 5Observed results of the basic observation and observations made for examination atsounding point number 400 at Dongguashan Copper Mine in Anhui Province观测模式AB=1600 mAB=3360 m视电阻 s/视频散率 Ps/%视电阻 s/视频散率 Ps/%基本观测0 244496 1050 192843 100259654 1060615633 10检查观测 10 244135 1050 196598 100293438 1060208383 10检查观测 20 243943 1050 195705 100323255 1060449890图 9冬瓜山铜矿实测供电电流强度和电位差波形图(片段)(a)AB=1600 m;(b)AB=3360 mFig 9Measured oscillographs of the supply current and the potential difference at Dongguashan Copper Mine in Anhui Province(a)AB=1600 m;(b)AB=3360 m为了探索为什么本文提出的“抗干扰参数”未取得有效的抗干扰效果，我们研究了冬瓜山实测的供电电流强度 I 和电位差 U 的“波形图”(图 9)由图可见，供电电流波形(图中红线所示)基本正常:正反向大体对称;只是供电电流强度随着时间增长缓慢减小 采用较小供电电极距观测的电位差信噪比较强时(图9a)，电位差波形(图中蓝线所示)大体上和供电电流波形相同;而采用较大供电电极距观测的电位差信噪比较弱时(图 9b)，图中蓝线所示电位差波形很复杂:除较大幅度的随机干扰外，还有幅度更大的跃变干扰 后者是引起很大观测误差的主要因素，至今还无法克服4结论本文讨论伪随机信号电法仪的抗干扰观测参数 仿真模拟、现场试验观测和野外试验结果表明，2842023，38(1)罗延钟，等:伪随机信号电法仪的抗干扰参数(www progeophys cn)基于供电电流和电位差时间序列在零偏移时间的互相关和自相关函数值比值 I，U(0)/I，I(0)，计算的导电性参数视电阻 s或视电阻率 s和激电参数视频散率 Ps，对于随机干扰、50 Hz 及其谐波的工频干扰和其他低频干扰等，具有很强的抗干扰能力 但是，对于与大功率直流电器(特别是矿山和坑道电气机车)开启、关断或变速相关的“跃变干扰”，仅在信噪比较高时，视电阻 s或视电阻率 s和激电参数视频散率 Ps才有较好的抗干扰能力;而在信噪比较低时，抗干扰效果欠佳 对于“跃变干扰”的抗干扰问题，尚待进一步研究解决致谢感谢审稿专家提出的修改意见和编辑部的大力支持!eferencesChu W T 1990 Impulse-response and reverberation-decay measurementsmade by using a periodic pseudorandom sequence Appl Acoust，29(3):193-205Chun S H，Chen J，Geng M H 2014 eview of the pseudo-random msequence and its application in electrical prospecting of explorationgeophysics Progress in Geophysics(in Chinese)，29(1):439-446，doi:10 6038/pg20140164Duncan P M，Hwang A，Edwards N，et al 1980 The developmentand applications of a wide band electromagnetic sounding systemusing a pseudo-noise source Geophysics，45(8):1276-1296Hobbs B，Ziolkowski A，Wright D 2006 Multi-transient electromagnetics(MTEM)-controlled source equipment for subsurface resistivityinvestigation/18th IAGA WG1 2 Workshop on ElectromagneticInduction in the Earth El Vendrell，Spain，17-23Li B N 1987 Pseudo andom Signal and Correlation Identification(inChinese)Beijing:Science PressLi M，Wei W B，Deng M 2008 Pseudo-random sequence in electricalprospecting Chinese Journal of Scientific Instrument(in Chinese)，29(4S):233-235Liu J K，Shen X，Zhao L 2013 Theory of System Identification andMATLAB Simunation(in Chinese)Beijing:Electronic Industry PressLiu J X 2010 Study on the m-sequence correlation identification and itsapplication to induced polarization prospecting for the method ofsuppressing noiseMasters thesis(in Chinese)Beijing:ChinaUniversity of GeosciencesLuo Y Z，Lu Z G，Sun G L，et al 2015 New generation of instrumentsfor electrical and electromagnetic prospecting instruments usingpseudo random signal Progress in Geophysics(in Chinese)，30(1):411-415，doi:106038/pg20150159Luo Y Z，Lu Z G，Sun G L，et al 2016 Algorithms of second solutionfor instruments of electrical and electromagnetic prospecting usingpseudo random signal Progress in Geophysics(in Chinese)，31(6):2604-2608，doi:10 6038/pg20160634Luo Y Z，Lu Z G，Sun G L，et al 2020 ecognition of anti-interferenceperforman