伪随机编码源激发下的时域电磁信号合成

王若， 王妙月， 底青云， 薛国强， 殷长春， 雷达

1 中国科学院页岩气与地质工程重点实验室，中国科学院地质与地球物理研究所， 北京 100029 2 吉林大学地球探测科学与技术学院， 长春 130026



伪随机编码源激发下的时域电磁信号合成

王若1， 王妙月1， 底青云1， 薛国强1， 殷长春2， 雷达1

1 中国科学院页岩气与地质工程重点实验室，中国科学院地质与地球物理研究所， 北京 100029 2 吉林大学地球探测科学与技术学院， 长春 130026

将伪随机编码技术引入到人工源电磁法后，可以通过加大发射功率以及应用后续的相关处理技术来达到压制噪声、加大探测深度及提高分辨率的目的，因此引起了越来越多学者的关注及研究，但大多数研究集中在资料处理的相关技术上，对模拟电磁信号关注较少.然而，资料处理工作大多是从电磁信号出发的，模拟伪随机编码源激发下的电磁信号不但可以为资料处理环节提供理论数据，而且可以为检测资料处理的效果提供中间结果，因此，模拟电磁信号工作必不可少.本文根据获得接收信号的物理过程来实现伪随机编码源激发下的电磁信号合成.首先用解析公式获得特定地电结构的大地频率域响应，然后通过余弦变换得到时间域阶跃响应，接下来用阶跃响应的时间导数得到大地脉冲响应，通过将大地脉冲响应与伪随机编码源的褶积得到理想接收信号，最后，用低通滤波器来模拟发射设备和接收设备的频带限制，将之和噪声一起加到理想接收信号上，最终模拟出仿真的合成信号.通过和野外实际接收信号对比发现本文合成信号仿真度较高，可以服务于后续的数据处理环节.

伪随机编码技术； 时间域； 大地脉冲响应； 接收信号； 合成

1 引言

近年来，传统的电磁方法在能源、矿产资源、水资源、环境地质及工程地质勘察中都得到了广泛应用且发挥了重要作用(底青云等，2006；Xue et al., 2014).随着国民经济的发展，需要探测埋深更深、精度更高的目标体时，这些方法的应用受到了一定的限制，亟待开发新的技术来提高探测精度与探测深度.

近年来，一些新技术新方法被陆续提出来，英国爱丁堡大学的Wright(2002)提出了多道瞬变电磁法的概念和探测油气目标体的处理技术，何继善(2010)提出了可提高分辨率及加大勘探深度的广域电磁法及相关技术.

在这些新方法中，都引入了伪随机编码发射技术，其全新的处理方式及应用效果引起各个研究机构的关注.近几年，中国科学院、中南大学、中国地质大学、国土资源部以及相关单位对这种发射技术(赵璧如等，2006；罗先中等，2014；罗延钟等，2015)及后续的处理工作都进行了不同程度的研究.汤井田和罗维斌(2008)分析了逆重复m序列的波形特征、相关函数及频谱特性，讨论了对各种干扰的压制方法及测量精度，薛国强等(2015)分析了多通道瞬变电磁法(MTEM)伪随机编码发射的响应特性和方法的关键技术，齐彦福等(2015)利用方波响应移位叠加和电流导数与阶跃响应褶积两种方法实现了理论m序列和实际发射波形的全时正演模拟，武欣等(2015)研究了m序列伪随机编码源电磁响应的精细辨识问题.

学者多把研究重点放在相关辨识及对噪声的压制上，对相关辨识的方法理论做了充分的探讨，但在模拟电磁信号方面却着墨较少.然而，采用伪随机编码源的资料处理工作大多是从电磁信号出发的，模拟伪随机编码源激发下的电磁信号不但可以为资料处理环节提供理论数据，而且可以为检测资料处理的质量提供中间结果，因此，模拟电磁信号工作是必不可少的.

本文根据获得接收信号的物理过程来实现伪随机编码源激发下的电磁信号合成.首先简述伪随机源激发下的电磁信号合成基本理论，然后根据合成公式中的大地脉冲响应、伪随机编码源、仪器带限影响及噪声逐一进行研究，最后给出信号合成结果，并与一个野外实测信号作对比，来验证信号合成的正确性.

2 基本理论

伪随机源激发下的电磁信号合成工作是指给定一个m序列伪随机编码源，在编码源的激发下，模拟大地及发射、接收系统和噪声的共同响应，继而合成一个仿真的野外观测信号的过程.

在野外工作中，发射机通过两个接地电极向地下供入伪随机交变电流信号，在离发射电极一定距离处用不极化电极排列接收相邻电极间的电位差，根据这种发射接收装置，发射机发射的电磁信号在大地中传播，然后通过电极与导线、接收系统接收.在勘探工作中，大地是探测目标，经过大地的信号为有用信号.对有用信号产生影响的因素主要为噪声，除此之外，电极附近的极化效应、接收器的频率响应限制(即带限)、电线的感抗等都对有用信号或多或少产生影响.在这些影响因素中，除噪声外，每一部分的影响都相当于在有用信号上加载了一个滤波器，最后接收到的信号是对有用信号经过多次滤波后的综合信号.因此，模拟电磁信号工作应当考虑各种因素的影响，然而，现实中有些因素的影响不易模拟或相对较小，本文只考虑噪声及仪器带限影响，并且将发射和接收系统的带限影响综合考虑.

接收信号的合成公式可以用(1)式表示(Ziolkowski et al., 2007)：

v(t)=i(t)*g(t)*r(t)+n(t)，

(1)

其中，t表示接收时间，v(t)表示接收机在接收点接收的信号，i(t)表示发射系统向大地发射的电流，g(t)表示经由大地滤波后在接收点处的大地脉冲响应，r(t)表示仪器频率响应.n(t)为接收点处的噪声.

从公式(1)可以看出，只要获得了模型地电结构的大地脉冲响应g(t)、伪随机发射源信号i(t)、仪器带限响应r(t)及噪声n(t)，就可合成出仿真的电磁信号v(t).下面对公式(1)中各项的实现过程逐一阐述.

3 大地脉冲响应

理想情况下，向地下输入脉冲电源，在不计及各种影响因素时，接收机测量的输出即大地脉冲响应.但在野外工作中，一般采用方波作为输入.在数值模拟时，先在频率域获得不同频率的电磁响应，然后通过余弦变换转换到时间域，得到大地阶跃响应，再通过计算阶跃响应的时间导数来得到大地脉冲响应.

3.1 频率域响应

本文用均匀半空间地电模型来模拟地表电场信号的合成过程.根据文献(Wright，2004)，只有与源同方向的电场可提供有用信息，其他方向的场不能提供更多的信息，因此本文只用与源同方向的电场来合成所需信号.假设源的方向为x方向，则频率域中同方向的电场可写为(朴化荣，1990)：

[3cos2φ-2+(1+k1r)e-k1r]，

(2)

用于频率域到时间域变换的方法较多，如正弦变换或余弦变换(Anderson,1982；Newman et al., 1986; Li and Constable,2010)、折线法(Li et al., 2011)，这三种方法同属于傅里叶逆变换.在这三种方法中，正弦变换与余弦变换需要计算较宽频带内的频率响应，其频带宽度与时间点有关；折线法相对来说计算速度快，但相对粗糙.除此之外，还有基于拉普拉斯逆变换的系列算法(Li and Farquharson,2015)可以实现信号的频时域转换，这类算法的优势在于可以用较少的滤波系数得到较精确的结果，但其对计算机精度的要求太高(4倍精度).本文最终选择了余弦变换作为把信号从频率域转换到时间域的方法.

余弦变换的表达式为(朴化荣，1990)：

(3)

其中，f(t)为时间域的场，F(ω)为频率域的场.

参考Anderson (1982)的方法，选用函数对

(4)

来求余弦变换的滤波权系数，其中，a为大于0的实数.利用滤波系数可将(3)式写为：

}/t，

(5)

其中Ci为余弦变换的滤波权系数，Ai-x为移动的横坐标，N1和N2的大小由积分时间自动调节，本文所用的余弦变换子程序为公开的代码(方文藻等，1993).

3.3 大地脉冲响应

殷长春等(2013)利用阶跃响应与电流的时间导数褶积代替脉冲响应与电流褶积的方法成功实现了任意波形全时响应的计算，避免了记录时间t→0时大地脉冲奇异值的出现，但只考虑了电流变化感应出的电磁场，并没有考虑电流自身产生的直流场，所以还需再将直流场加进感应的场中.本文采用Strack(1992)的提法，直接计算阶跃响应的时间导数来得到大地脉冲响应.t→0处的值用第一个阶跃值与该处时间的比值来代替，通过将起始时间设计得较小(如0.01 ms)来减小t→0时的误差.

设计一个电阻率是100 Ωm的均匀半空间模型，假设源的方向为x方向，源的长度为200 m，计算点位于源的延长线上且离源中心点500 m处，电场的频率域响应曲线如图1a所示.图1a显示了从10-7Hz到107Hz的计算结果，这个频带宽度可满足1 s内所有时间点上进行余弦变换时所需的频率范围.从图1a可以看出，在高频和在低频一定范围内，电场曲线保持水平，不受频率的影响.通过余弦变换将频率域结果转换到时间域，得到如图1b所示的阶跃响应，将其与解析公式计算结果对比，发现二者的吻合程度非常高，验证了余弦变换过程的正确性.采用3.3节中提到的阶跃响应的时间导数来计算脉冲响应，计算结果如图1c所示，对图中所示的x方向电场来说，曲线中出现了一个峰值，这个峰值对后续的处理非常重要，因不是本文的研究内容，所以不进行阐述.这里为了突出峰值的形态，对曲线进行了截断.所用的时间域解析公式来自文献(朴化荣，1990).

图1 均匀半空间模型的频率域与时间域响应(a) 频率域大地响应曲线； (b) 时间域大地阶跃响应； (c) 时间域大地脉冲响应.Fig.1 Responses of homogenous half space in frequency domain and time domain(a) Response curve in frequency domain; (b) Step response of earth in time domain; (c) Impulse response of earth in time domain.

4 m序列伪随机源

4.1 m序列伪随机编码

伪随机序列是由移位寄存器产生的， 根据文献(林可祥和汪一飞，1978)，可通过查表的方式得到不同级移位寄存器产生的伪随机序列的本原多项式F(x).通过F(x)与序列多项式G(x)的关系式G(x)=1/F(x)，用多项式的长除算法可计算出相应的序列多项式G(x)，从G(x)中可直接获得m序列伪随机编码.

序列多项式G(x)最终可写为下面的表达形式：

(6)

其中⊕表示模二加法运算符，x为移位运算符，其上标的数字表示移位的位数，例如：akx2表示将ak位移两位，an所组成的序列{an}=a0,a1,a2,a3,…即为m序列伪随机编码.(5)式中的N为伪随机码一个周期的码元个数，N=2r-1，r为移位寄存器级数.

图2是4级和6级的m序列伪随机编码一个周期的波形，它是通过文献(林可祥和汪一飞，1978)本原多项式列表中4级和6级第一个本原多项式得到的m序列的编码波形，从图中看出其随机性随着移位寄存器级数的增多而愈发明显.

图2 m序列伪随机编码一个周期的波形(a) 4级移位寄存器产生的波形； (b) 6级移位寄存器产生的波形.Fig.2 Waveforms of pseudo random binary sequence(a) Generated by 4-rank shift register; (b) Generated by 6-rank shift register.

4.2 m序列伪随机编码源的频谱

对m序列伪随机编码应用快速傅氏变换，获得了其作为源时相应的频谱.通过对比不同级数的频谱幅值，可了解不同级数的伪随机信号所包含的频率成分及幅频特性.本文给出了4、6、8、10、12级伪随机编码源的频谱，其中与伪随机编码相关的参数见表1.为了对比不同源谱的宽度，将所有伪随机编码周期与采样率都定为0.25 s及16384 Hz.伪随机编码源的频谱如图3所示.

表1 伪随机编码相关参数

从图3可以看出，不同级移位寄存器得到的频谱形态不同，级数越高，主瓣的宽度越宽，当移位寄存器的级数r=12时，从0 Hz到8192 Hz范围内的频谱丰富，包络线的形态接近水平，接近于脉冲源的频谱，由此源激发得到的大地响应会接近于脉冲源

图3 不同阶m序列伪随机编码源发射波形的频谱特性(a) 4级移位寄存器； (b) 6级移位寄存器； (c) 8级移位寄存器； (d) 10级移位寄存器； (e) 12级移位寄存器.Fig.3 Spectral properties of pseudo-random binary sequence generated by shift registers of different ranks(a) 4-rank; (b) 6-rank; (c) 8-rank; (d) 10-rank; (e) 12-rank.

激发得到的响应.图3中从10级到12级伪随机编码源的频谱外包络线不光滑，这是采样率不够导致的.若提高采样率，则能得到外包络线较光滑、幅值较均匀的谱图.从图3还可以看出，随着寄存器级数的增加，主瓣的幅值也随之下降，说明信号的能量降低.对地球物理勘探来说，信号能量降低会导致电磁信号的穿透深度变浅.换句话说，虽然主瓣频率成分提高，分辨率也提高，但信号强度因能量分散而降低，导致探测深度浅.若想解决这个问题，需在加大发射功率的同时增加垂直叠加的时间，而且在后续处理采用类地震水平叠加的技术.

5 其他影响因素

5.1 仪器频带限制的影响

发射机和接收器自身的频带宽度是有限的，作为理论研究，可以用有一定带宽的低通滤波器来模拟仪器的频带宽度(即带限问题).

低通滤波器的设计方法有多种，不同方法得到的滤波器的旁瓣幅值大小以及从主瓣到旁瓣的过渡带宽度不同.理想滤波器要求过渡带尽可能窄，旁瓣幅值尽可能小，以减小能量的泄露，使能量尽可能多地保留在主瓣内.本文选用海明窗函数法作为滤波器的设计方法，因为海明窗可将99.963%的能量集中在窗谱的主瓣内，旁瓣的峰值小于主瓣峰值的1%(程佩青，2007)，符合理想滤波器的要求.

海明窗函数为：

(7)

借鉴频率域电磁仪器的频带宽度，设计一个低通滤波器，其通带为8192Hz，阻带为9600Hz.理想的低通滤波器的频带如图4a所示，在通带内，信号可以通过，通带外是阻带，阻带内的信号全部被滤除.理想的低通滤波器在时间域内如图4c中的黑实线所示(实际为无限信号，为了和加窗后的信号相比，只显示了其中的一部分).图4b是海明窗在频率域的表现形式，用海明窗对图4c中的信号进行截断，截断后的信号如图4c中的离散点所示，将加窗后的时间域滤波器变换到频率域，如图4d所示，相对于理想低通滤波器的频带(图4a)，图4d出现了通带与阻带之间的过渡带，但较窄，阻带出现了小幅

图4 低通滤波器的特性(a) 理想的低通滤波器的谱; (b) 海明窗; (c) 理想的低通滤波器与经过海明窗截断后的滤波器对比; (d) 加窗后的信号对应的低通滤波器的谱.Fig.4 Properties of low-pass filter(a) Ideal spectrum; (b) Hamming window; (c) Comparison of ideal low-pass filter and filter cut-off by Hamming window; (d) Spectrum of low-pass filter after cutting.

振荡(因其相对于主瓣来说，能量较小，振荡不明显，若放大该部分，可以看到振荡现象)，说明所设计的滤波器过渡带较窄，阻带能量泄漏很小，是一个性能优良的滤波器.

5.2 噪声的加入

当用电磁法在野外工作时，常遇到的噪声类型有两种，一是白噪声，二是50 Hz的工频噪声.

(1) 白噪声

本文用正态分布随机数来模拟噪声：

(8)

其中，n1为所求的正态分布随机数，μ为正态分布的均值，本文为观测资料，σ2为正态分布的方差，在本文中σ为观测值与噪声水平之积，Ni为0到1之间均匀分布的随机数，m为随机数的个数.

(2) 50 Hz噪声

用一个频率是50 Hz的正弦波来模拟工频噪声：

n2=Asin(100πt),

(9)

其中n2为噪声值，A为50 Hz工频噪声的振幅.

6 伪随机源激发下的电磁信号合成

6.1 信号合成

通过(1)式将以上内容进行组合，便可分别得到理想接收信号、加带限的信号和带噪声的信号.

以4阶m序列伪随机编码源作为激发源(见图5a)，将源和均匀半空间大地脉冲响应(见图1c)做褶积，得到电场的理想信号，如图5b所示.将理想信号与低通滤波器(仪器带限)做褶积，便可得到考虑了带限影响的电场信号，如图5c所示.从这些图可以看出，不加带限前，信号比较平滑，加了带限后，信号出现了一些毛刺状的干扰.

对加了带限的信号分别加上了5%及10%的白噪声，如图6a和6b所示，噪声叠加在原来的信号之上，引起信号的畸变，当噪声达到10%时，形成的干扰已非常明显.图6c为幅值是4倍于有用信号的50 Hz干扰信号，干扰信号和有用信号叠加后的信号如图6d所示，可以看出，有用信号湮没在了50 Hz干扰信号中.

6.2 合成信号与实测信号的对比

中国科学院项目组2015年在河北张北地区进行了野外试验.用一个长度是300 m的接地电偶极发射一个12阶的伪随机编码源，采用轴向装置，用一系列极距是60 m的电偶极子一字排开，在离发射源不同距离处的接收端接收耦合了大地响应及各种影响的信号.为了和野外信号对比，在室内也进行了模拟.根据测区已知资料，在室内用电阻率是100 Ωm的均匀半空间模型来近似地下介质，发射源也用12级移位寄存器生成的m序列伪随机编码序列，编码的部分形态如图7a所示.采用了和野外一致的码元长度及采样率(码元长度为1/1024 s，采样率为16 kHz，发射一个周期的信号大约需要4 s)，因信号太长，为了突出细节，我们仅显示了其中的一小部分.值得注意的是，室内的编码序列和室外的编码序列是不一致的，因为对于任何级的移位寄存器来说，可以产生很多种形态编码，级数越多，能产生的编码形态也越多，所以并未刻意和野外用相同的编码形式.模拟时，分别在偏移距由近及远的位置接收信号，300 m和2700 m处的脉冲响应分别如图7b和图7c所示，离源较近的接收点脉冲响应的信号很强，并且和离源较远时接收到的脉冲响应形态有差异，这是因为离源近，信号接近饱和，离源越远，得到的脉冲响应越平缓.

图5 在4阶伪随机源激发下的接收信号(a) 4阶伪随机源; (b) 不考虑带限的电场理想信号; (c) 考虑带限的电场信号.Fig.5 Signal excited by 4-rank pseudo-random binary sequence(a) Source; (b) Ideal signal; (c) Signal considering band limit.

图6 在4阶伪随机源激发下考虑噪声的信号(a) 加5%白噪声; (b) 加10%白噪声; (c) 50 Hz干扰噪声; (d) 加50 Hz干扰噪声.Fig.6 Signal with different noise type(a) Adding 5% white noise; (b) Adding 10% white noise; (c) 50 Hz industry noise; (d) Adding 50 Hz industrial noise.

图7 模拟信号和野外实测信号的对比(a) 12级m序列伪随机信号部分波形; (b) 偏移距是300 m时的电场脉冲响应; (c) 偏移距是2700 m时的电场脉冲响应; (d) 偏移距是300 m时的电场计算信号; (e) 偏移距是300 m时的电场实测信号.Fig.7 Comparison of simulated signal to survey signal(a) Source; (b) Impulse response at offset 300 m; (c) Impulse response at offset 2700 m; (d) Simulated signal at offset 300 m; (e) Measured signal at offset 300 m.

由于本文所依据的理论是电偶极子理论，当发射偶极与接收偶极子距离较近时，对发射偶极源进行了离散叠加.对300 m的收发距，将偶极源划分成一系列5 m长的小偶极子，离散后的小偶极子是300 m收发距的1/60，满足偶极子理论，同时大大减小了线间的感应耦合.

我们计算了部分点的接收信号，并和野外接收信号做了对比，发现二者的相似程度较高，这里只展示偏移距是300 m时的合成信号，如图7c所示，通过和野外300 m处接收信号的屏幕截图7e对比，可以看出信号的相似度较高.二者不尽相同的原因有三个：一是所用的编码形式不同，二是地下介质的电阻率不同，这两个在前面也有所描述，还有一个原因，就是在野外实测中所受的干扰更多，如电极附近的极化效应，接收电线的感抗，而在数值模拟中只是模拟了较为典型的干扰.即便如此，二者较高的相似度说明了本文信号合成过程是正确的.

7 结论

本文根据伪随机编码源信号的传播物理过程对接收信号进行了合成，通过研究，发现大地脉冲响应离源较近时会出现饱和现象，是一条单调衰减曲线，随着测点远离源，峰值逐渐出现.

本文通过计算阶跃响应的时间导数来得到大地脉冲响应，t→0处的值用第一个阶跃值与该处时间的比值来代替，通过将起始时间设计得较小(如0.01 ms)来减小t→0时的误差，通过数值模拟发现，这种近似方法可行.

当用移位寄存器得到m序列伪随机编码时，移位寄存器级数越多，其频谱主瓣的宽度越宽，当移位寄存器的级数r=12时，从0～8192 Hz范围内的各个频率点上的幅值外包络线形态接近水平，接近于脉冲响应的频谱，由此源激发而得到大地响应会接近于脉冲源激发得到的响应，从而有较高的电性结构分辨能力.然而随着寄存器级数的增加，主瓣的各个频率信号的幅值也随之下降，影响了勘探深度的提高，为此需通过增加源的发射时间(便于加大垂直叠加次数)和增加拟地震水平叠加次数来弥补因增加频率成分而导致的信号强度损失，以达到既提高分辨率，又提高探测深度的双重目标.

由本文的方法得到的大地脉冲响应与理论值吻合，合成信号与野外信号非常接近，因此，本文方法所产生的中间结果及理论信号合理可靠，可用在后续的数据处理中，并可用于检验后续数据处理方法的可靠性和有效性.

Anderson W L. 1982. Calculation of Transient Soundings for a Coincident Loop System. Denver, Colo.: U. S. Geological Survey.Cheng P Q. 2007. Digital Signal Processor Tutorial. 3rd ed. (in Chinese). Beijing: Tsinghua University Press.

Di Q Y, Wang G J, An Z G, et al. 2006. Geophysical exploration of a long deep tunnel on the west route of South to North Water Diversion Project.ChineseJ.Geophys. (in Chinese), 49(6): 1836-1842.Fang W Z, Li Y G, Li X. 1993. Theory of Transient Electromagnetic Sounding (in Chinese). Xi′an: Northwestern Polytechnical University Press. He J S. 2010. Full Domain Electromagnetic Method and Pseudo-Random Signal Electric Method (in Chinese). Beijing: China Higher Education Press.

Li J H, Farquharson C G. 2015. Two effective inverse Laplace transform algorithms for computing time-domain electromagnetic responses.∥ SEG Technical Program Expanded Abstracts. SEG: 957-962.Li J H, Zhu Z Q, Liu S C, et al. 2011. 3D numerical simulation for the transient electromagnetic field excited by the central loop based on the vector finite-element method.JournalofGeophysicsandEngineering, 8(4): 560-567.Li Y G, Constable S. 2010. Transient electromagnetic in shallow water: insights from 1D modeling.ChineseJ.Geophys., 53(3): 737-742, doi: 10.3969/j.issn.0001-5733.2010.03.029.

Lin K X, Wang Y F. 1978. Theory and Application of Pseudo-Random Code (in Chinese). Beijing: Posts and Telecom Press.

Luo X Z, Li D W, Peng F P, et al. 2014. Implementation and applications of an coded electrical instrument with anti-interference ability.ProgressinGeophysics(in Chinese), 29(2): 944-951, doi: 10.6038/pg20140263.

Luo Y Z, Lu Z G, Sun G L, et al. 2015. New generation of instruments for electrical and electromagnetic prospecting-instruments using pseudo random signal.ProgressinGeophysics(in Chinese), 30(1): 411-415, doi: 10.6038/pg20150159.Newman G A, Hohmann G W, Anderson W L. 1986. Transient electromagnetic response of a three-dimensional body in a layered earth.Geophysics, 51(8): 1608-1627.

Piao H R. 1990. Theory of Electromagnetic Sounding (in Chinese). Beijing: Geological Publishing House.

Qi Y F, Yin C C, Wang R, et al. 2015. Multi-transient EM full-time forward modeling and inversion of m-sequences.ChineseJ.Geophys. (in Chinese), 58(7): 2566-2577, doi: 10.6038/cjg20150731.Strack K M. 1992. Exploration with Deep Transient Electromagnetics. Netherlands: Elsevier Science Publishers.Tang J T, Luo W B. 2008. Pseudo-random electromagnetic exploration based on invert-repeated m-Sequence correlation identification.ChineseJ.Geophys. (in Chinese), 51(4): 1226-1233.

Wright D, Ziolkowski A, Hobbs B A. 2002. Hydrocarbon detection and monitoring with a multicomponent transient electromagnetic (MTEM) survey.TheLeadingEdge, 21(9): 852-864. Wright D A. 2004. Detection of hydrocarbons and their movement in a reservoir using time-lapse multichannel transient electromagnetic (MTEM) data [Ph. D. thesis]. Edinburgh: University of Edinburgh.Wu X, Xue G Q, Di Q Y, et al. 2015. Accurate identification for the electromagnetic impulse response of the earth with pseudo random coded waveforms.ChineseJ.Geophys. (in Chinese), 58(8): 2792-2802, doi: 10.6038/cjg20150814.

Xue G Q, Gelius L J, Sakyi P A, et al. 2014. Discovery of a hidden BIF deposit in Anhui province, China by integrated geological and geophysical investigations.OreGeologyReview, 63: 470-477.

Xue G Q, Yan S, Di Q Y, et al. 2015. Technical analysis of multi-transient electromagnetic method.JournalofEarthSciencesandEnvironment(in Chinese), 37(1): 94-100.

Yin C C, Huang W, Ben F. 2013. The full-time electromagnetic modeling for time-domain airborne electromagnetic systems.ChineseJ.Geophys. (in Chinese), 56(9): 3153-3162, doi: 10.6038/cjg20130928.

Zhao B R, Zhao J, Zhang H K, et al. 2006. The PS100 high precision earth electricity instrument system (IP to IP) with controllable source-Application of CDMA technology to the measurement of earth resistivity for the first time.ProgressinGeophysics(in Chinese), 21(2): 675-682.

Ziolkowski A, Hobbs B A, Wright D. 2007. Multitransient electromagnetic demonstration survey in France.Geophysics, 72(4): F197-F209.

附中文参考文献

程佩青. 2007. 数字信号处理教程(第三版). 北京: 清华大学出版社.

底青云, 王光杰, 安志国等. 2006. 南水北调西线千米深长隧洞围岩构造地球物理勘探. 地球物理学报, 49(6): 1836-1842.

方文藻, 李予国, 李貅. 1993. 瞬变电磁测深法原理. 西安: 西北工业大学出版社.

何继善. 2010. 广域电磁法和伪随机信号电法. 北京: 高等教育出版社.

林可祥, 汪一飞. 1978. 伪随机码的原理与应用. 北京: 人民邮电出版社.

罗先中, 李达为, 彭芳苹等. 2014. 抗干扰编码电法仪的实现及应用. 地球物理学进展, 29(2): 944-951, doi: 10.6038/pg20140263.

罗延钟, 陆占国, 孙国良等. 2015. 新一代电法勘查仪器——伪随机信号电法仪. 地球物理学进展, 30(1): 411-415, doi: 10.6038/pg20150159.

朴化荣. 1990. 电磁测深法原理. 北京: 地质出版社.

齐彦福, 殷长春, 王若等. 2015. 多通道瞬变电磁m序列全时正演模拟与反演. 地球物理学报, 58(7): 2566-2577, doi: 10.6038/cjg20150731.

汤井田, 罗维斌. 2008. 基于相关辨识的逆重复m序列伪随机电磁法. 地球物理学报, 51(4): 1226-1233.

武欣, 薛国强, 底青云等. 2015. 伪随机编码源电磁响应的精细辨识. 地球物理学报, 58(8): 2792-2802, doi: 10.6038/cjg20150814.

薛国强, 闫述, 底青云等. 2015. 多道瞬变电磁法(MTEM)技术分析. 地球科学与环境学报, 37(1): 94-100.

殷长春, 黄威, 贲放. 2013. 时间域航空电磁系统瞬变全时响应正演模拟. 地球物理学报, 56(9): 3153-3162, doi: 10.6038/cjg20130928.

赵璧如, 赵健, 张洪魁等. 2006. PS100型IP到端可控源高精度大地电测仪系统——CDMA技术首次在地电阻率测量中的应用. 地球物理学进展, 21(2): 675-682.

(本文编辑 何燕)

Signal synthesis of time-domain electromagnetic response excited by a pseudo random binary sequence source

WANG Ruo1, WANG Miao-Yue1, DI Qing-Yun1, XUE Guo-Qiang1, YIN Chang-Chun2, LEI Da1

1KeyLaboratoryofShaleGasandEngineeringGeology,InstituteofGeologyandGeophysics,ChineseAcademyofSciences,Beijing100029,China2CollegeofGeo-explorationScienceandTechnology,JilinUniversity,Changchun130026,China

After the pseudo-random binary sequence source is introduced into the artificial electromagnetic method, the exploration depth and the resolution are improved, because in this approach correlation technique can be used in the data processing to suppress noise. Most researchers didn’t focus on the synthetic signal rather than data processing. However, some data processing techniques have to use the electromagnetic signal, therefore it is indispensable to rely on synthesis of electromagnetic signal. In this paper, the electromagnetic signal is synthesized according to the propagating routine of the electromagnetic field. Firstly, electromagnetic response of a simulated geo-electric model in the frequency domain is calculated according to analytic formula. Secondly, the step electromagnetic response in the time domain is obtained by applying cosine transform to the response in the frequency domain, and then impulse response of the earth is obtained from time derivative of the step response. Through convoluting the impulse response of the earth and with PRBS source, the ideal electromagnetic signal can be composed. Finally, a low-pass filter is used to simulate the synthetic system response of the transmitter and receiver because all of the instruments are band-limited. Then the simulated signal is obtained by convoluting a low pass filter to the ideal signal added with noise. Comparing with actual field received signal proves that the composed signal is correct, and the signal of high simulation can be used to serve the data processing technique.

Pseudo randomly coded technology; Time domain; Impulse response; Received signal; Compose

10.6038/cjg20161203.

国家重大科研装备研制项目“深部资源探测核心装备研发”(ZDYZ2012-1-05-04)，国家重点研发计划(2016YFC0600507)，国家自然科学基金(41174111)资助.

王若，女，博士，主要从事电磁法正反演研究.E-mail：wangruo@mail.iggcas.ac.cn

10.6038/cjg20161203

P631

2016-02-01，2016-11-14收修定稿

王若， 王妙月， 底青云等. 2016. 伪随机编码源激发下的时域电磁信号合成. 地球物理学报，59(12):4414-4423,

Wang R, Wang M Y, Di Q Y, et al. 2016. Signal synthesis of time-domain electromagnetic response excited by a pseudo random binary sequence source.ChineseJ.Geophys. (in Chinese),59(12):4414-4423,doi:10.6038/cjg20161203.