伪随机信号电法仪的抗干扰参数_罗延钟.pdf
《伪随机信号电法仪的抗干扰参数_罗延钟.pdf》由会员分享,可在线阅读,更多相关《伪随机信号电法仪的抗干扰参数_罗延钟.pdf(10页珍藏版)》请在咨信网上搜索。
1、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 pr
2、ospecting 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
3、 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
4、 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
5、)函数值较大时才有较强的抗干扰能力;而时间偏移 t 大于位宽 t 时,相关函数值非常小,易受干扰而畸变 利用全周期相关函数值计算大地(系统)时间响应和频率响应的伪随机相关辨识的抗干扰性,并不理想 在电磁干扰较大或信噪比较低(近于或小于 0 dB)时,不能取得可靠的时间响应和频率响应数据 基于上述情况,本文提出仅利用零偏移时间的互相关和自相关函数值 I,U(0)和 I,I(0),获取具有高抗干扰性能的导电和激电参数 文中给出了这些抗干扰参数的具体算法,并列举仿真模拟和野外观测数据,证明它们的实际抗干扰能力关键词电法勘查;观测仪器;伪随机信号;抗干扰参数中图分类号P631文献标识码Adoi:10
6、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-co
7、rrelation 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
8、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-jammin
9、g 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 freque
10、ncy 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 ofthe
11、se 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引言基于伪随机信号系统辨识的电法勘查
12、仪器多试图借助于相关处理 计算大地(系统)输入信号(供电电流 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;淳少恒等,20
13、14;张宇等,2014;罗延钟等,2015,2016)近期笔者(罗延钟等,2020)通过仿真模拟和野外试验,初步探索了“相关辨识”的抗干扰性能 得到的结果是,无论是否存在干扰(包括理论的随机干扰和实测的电磁干扰),相关谱分析算法计算的阻抗频谱 Z2(i)和直接谱分析算法计算的阻抗频谱 ZI,U(i),都几乎完全重合 这表明,相关处理并未明显增加伪随机信号观测数据计算大地阻抗频谱的抗干扰能力 究其原因是,仅在偏移时间t 小于伪随机信号位宽 t,特别是在零偏移时间,互相关 I,U(t)和自相关 I,I(t)函数值较大时才有较强的抗干扰能力;而时间偏移 t 大于位宽 t 时,相关函数值非常小,易受干
14、扰而畸变;致使利用全周期相关函数值计算大地(系统)时间响应和频率响应的伪随机相关辨识的抗干扰性,并不理想 在电磁干扰较大或信噪比较低(近于或小于 0 dB)时,不能取得可靠的时间响应和频率响应数据 本文将在上述研究成果的基础上,提出仅利用零偏移时间的互相关和自相关函数值 I,U(0)和 I,I(0),获取具有高抗干扰性能的导电和激电参数;并用仿真模拟和野外试验结果,验证这些参数的实际抗干扰效果1伪随机信号电法仪观测的“视电阻”受前文(罗延钟等,2020)“零偏移时间的互相关和自相关函数值 I,U(0)和 I,I(0)具有高抗干扰性能”的启发,我们计算两者的比值 该比值具有电阻量纲,我们定义为“
15、视电阻”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与
16、电容 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
17、的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 示出了采用不
18、同阶数和位宽的 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可反
19、映测区大地的导电性;将视电阻 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零偏移
20、互相关 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.017145
21、0.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%时,视电阻和视频散
22、率所受的影响仍很小;尤其是在两个 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,
23、第一位宽 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+0304447
24、51E+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 电路模型作仿真模拟,
25、计算用单 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
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 随机 信号 电法仪 抗干扰 参数 罗延钟
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。