基于加权平均法的高密度电法数据预处理--《陕西省地球物理文集（八）：深部矿井灾害源探测实践》2008年
基于加权平均法的高密度电法数据预处理
【摘要】：应用加权平均的方法对高密度电法的某些突变数据进行预处理,对处理后的高密度电法资料进行重构,便可得到消除了突变后的数据。经过验证,该方法应用于高密度电法资料的预处理对提高资料质量有较好的效果。
【作者单位】：
【关键词】：
【分类号】：P631.34【正文快照】：
高密度电阻率法是近年来应用较为广泛的地球物理勘探方法之一,与常规电阻率法相比,高密度电阻率法在野外数据采集过程中可组合使用多达十余种装置形式,具有观测精度高、数据采集量大、地质信息丰富、分辨率高、生产效率高等特点。在实际探测工作中,由于人为错误、外界强干扰、
欢迎：、、)
支持CAJ、PDF文件格式，仅支持PDF格式
【参考文献】
中国期刊全文数据库
熊光楚;[J];物探化探计算技术;2000年02期
刘荣昌,薛隆泉,刘宏昭,王慧武,王庆祝;[J];西安理工大学学报;2002年04期
【共引文献】
中国期刊全文数据库
熊光楚;[J];物探化探计算技术;2000年03期
吴国平,彭放,徐忠祥;[J];物探化探计算技术;2001年04期
刘海飞;张赛民;阮百尧;;[J];物探化探计算技术;2006年03期
中国博士学位论文全文数据库
杨福俊;[D];大连理工大学;2002年
钟清;[D];中国地质大学（北京）;2006年
刘晓光;[D];上海交通大学;2008年
【二级参考文献】
中国期刊全文数据库
邵梦林,陆正立,梁锦文;[J];桂林冶金地质学院学报;1982年02期
管志宁;安玉林;;[J];物化探计算技术;1984年01期
【相似文献】
中国期刊全文数据库
王宇玺;肖宏跃;雷宛;贺智娟;;[J];工程勘察;2009年11期
祝杰;亢会明;杜毅;雷宛;;[J];勘察科学技术;2011年02期
罗高玲;陈鹏飞;罗起;;[J];中国科技信息;2011年16期
祝杰;杜毅;亢会明;雷宛;张银松;;[J];工程勘察;2011年10期
侯烈忠,秋兴国,罗奕;[J];煤田地质与勘探;1997年04期
王士鹏;[J];水文地质工程地质;2000年01期
王运生,王旭明,谢向文,孙雅芳,冷元宝,郭玉松;[J];勘察科学技术;2001年01期
胡勘,赵贵民;[J];贵州地质;2002年04期
龚道平;[J];华东地质学院学报;2002年03期
龚道平;[J];勘察科学技术;2002年05期
中国重要会议论文全文数据库
马利平;徐晓英;;[A];山东地球物理六十年[C];2009年
董浩斌;王传雷;;[A];资源、环境与工程电法勘探方法技术交流会论文集[C];2003年
王建平;;[A];中国地球物理.2003——中国地球物理学会第十九届年会论文集[C];2003年
吕玉增;阮百尧;;[A];勘探地球物理2005学术交流会论文集[C];2005年
黄宴辉;王诗东;熊壮;朱世山;;[A];中国地球物理2013——第二十四专题论文集[C];2013年
祁民;张宝林;梁光河;陈友明;管刚;;[A];中国科学院地质与地球物理研究所2006年论文摘要集[C];2007年
祁民;张宝林;梁光河;;[A];中国科学院地质与地球物理研究所2006年论文摘要集[C];2007年
董浩斌;王传雷;唐筱眸;;[A];2001年中国地球物理学会年刊——中国地球物理学会第十七届年会论文集[C];2001年
刘晓东;张虎生;邢应太;;[A];资源、环境与工程电法勘探方法技术交流会论文集[C];2003年
寇玉才;王俊祥;程洪明;;[A];勘探地球物理2005学术交流会论文集[C];2005年
中国重要报纸全文数据库
邱宏喜 张志臣;[N];地质勘查导报;2007年
中国硕士学位论文全文数据库
梁冰;[D];吉林大学;2010年
雷世红;[D];河海大学;2005年
钟韬;[D];成都理工大学;2008年
张碧勇;[D];吉林大学;2013年
郭清石;[D];西南交通大学;2013年
蔡斌;[D];吉林大学;2011年
刘海生;[D];太原理工大学;2006年
何晓辉;[D];郑州大学;2012年
侯智超;[D];中南大学;2013年
张永伟;[D];华东交通大学;2014年
&快捷付款方式
&订购知网充值卡
400-819-9993
《中国学术期刊（光盘版）》电子杂志社有限公司
同方知网数字出版技术股份有限公司
地址：北京清华大学 84-48信箱 大众知识服务
出版物经营许可证 新出发京批字第直0595号
订购热线：400-819-82499
服务热线：010--
在线咨询：
传真：010-
京公网安备75号一种集中式高密度电法仪的数据采集控制方法
申请号：.1 申请日：
摘要：本发明提供一种集中式高密度电法仪的数据采集控制方法，包括：步骤S1，获取特殊场合下集中电缆和电极的布置图；步骤S2，设计特殊场合下供电电极A、B和测量电极M、N的测量方式；步骤S3，设计供电电极A、B所在的电极编号和测量电极M、N所在的电极编号；步骤S4，创建文本格式的文本装置文件，文本装置文件包括文件头和根据供电电极A、B的电极编号、测量电极M、N的电极编号以及测点序号填写的电极排列方式；步骤S5，将创建的文本装置文件存储到高密度电法仪测量程序所在目录。根据本发明获得的任意装置的文本装置文件，可以按用户需求设置任意装置任意电极排列方式，极大地拓展了集中式高密度电法仪的适用范围。
地址： 518000 广东省深圳市南山区高新区南区虚拟大学园楼A208室
发明(设计)人：
主分类号：
&实质审查的生效IPC(主分类):G01V
3/00申请日:
注：本法律状态信息仅供参考，即时准确的法律状态信息须到国家知识产权局办理专利登记簿副本。
&一种集中式高密度电法仪的数据采集控制方法，包括：步骤S1，获取特殊场合下集中电缆和电极的布置图；步骤S2，设计所述特殊场合下供电电极A、B和测量电极M、N的测量方式；步骤S3，设计所述供电电极A、B所在的电极编号和所述测量电极M、N所在的电极编号；步骤S4，创建文本格式的文本装置文件，所述文本装置文件包括文件头和根据所述供电电极A、B的电极编号、所述测量电极M、N的电极编号以及测点序号填写的电极排列方式；步骤S5，将创建的所述文本装置文件存储到高密度电法仪测量程序所在目录。
专利代理机构
进入国家日期抗干扰编码电法仪的实现及应用
&&&&2014, Vol. 29 Issue (2): 944-951
罗先中, 李达为, 彭芳苹, 罗林涛, 赵金水. 2014,抗干扰编码电法仪的实现及应用[J]. 地球物理学进展,
29(2): 944-951, doi: 10.6038/pg&&
LUO Xian-zhong, LI Da-wei, PENG Fang-ping, LUO Lin-tao, ZHAO Jin-shui.2014, Implementation and applications of an coded electrical instrument with anti-interference ability. Progress in Geophysics,
29(2): 944-951, doi: 10.6038/pg &&
抗干扰编码电法仪的实现及应用
罗先中1,2,3, 李达为4, 彭芳苹4, 罗林涛4, 赵金水4&&&&
1. 油气资源与勘探技术教育部重点实验室（长江大学）, 武汉 430010;2. 中国地质调查局发展研究中心, 北京 100037;3. 国土资源部矿产勘查技术指导中心, 北京 100120;4. 北京勘察技术工程有限公司, 北京 100085
基金项目：国家自然科学基金面项目“伪随机扩频脉冲在地球介质中的传播特性研究”（）和国土资源部公益性行业科研专项基金项目“深井物探和抗干扰电法技术研究与应用示范”（）资助.
作者简介：罗先中，男，1968年出生，1991年毕业于中国地质大学（武汉），教授级高级工程师，长期从事电磁法理论、抗干扰电法仪器等研究.（E-mail：）
摘要：为了提高电法或电磁法仪器的抗干扰能力，采用伪随机编码的脉冲序列作为激励源，应用通讯领域的伪随机编码相关检测技术和信号过采样等先进技术，开发出具有强抗干扰能力的电法勘探硬件采集系统和资料处理软件系统.通过仿真检测和在不同勘探领域的实测试验，证明了该方法的可行性与有效性.实现的抗干扰编码电法仪达到了预期效果，展示了优越的抗干扰能力，可以用于城市、矿山、电站和铁路周边等强干扰地区的电法勘探.
伪随机编码&&&&
相关检测&&&&
过采样&&&&
抗干扰&&&&
电法仪器&&&&
Implementation and applications of an coded electrical instrument with anti-interference ability
LUO Xian-zhong1,2,3, LI Da-wei4, PENG Fang-ping4, LUO Lin-tao4, ZHAO Jin-shui4&&&&
1. Key Laboratory of Exploration Technologies for Oil and Gas Resources, Yangtzi University, Wuhan 430010, C2. Development and Research Center, China Geological Survey, Beijing 100037, C3. Mineral Exploration and Technical Guidance Center, Ministry of Land and Resources, Beijing 100120, C4. Beijing Explo-Tech Engineering Co., Ltd., Beijing 100085, China
Abstract: To improve the anti-interference ability of the electrical and electromagnetic survey instrument, a coded pseudo-random pulse sequence is used as the driving source, and a set of data acquisition hardware system for electric or electromagnetic survey has been developed. The algorithms of data processing and software implementation are based on the technologies of correlation detection of pseudo-random sequences and over sampling which have been used successfully in communication for anti-interference improvement. The effectiveness of the hardware and software system is proved through a series of laboratory simulation test and field tests of various survey sites, and superior anti-interference ability has been exhibited through comparison tests. It is suitable for surveys in strong anti-interference environment, such as mining area, near power station or railway, urban and suburb area, etc.
Key words:
pseudo-random coding&&&&
correlation superimposing&&&&
over sampling&&&&
anti-interference&&&&
electrical prospecting instrument&&&&
作为地球物理勘探的一种重要方法，电法勘主要针对勘探目标的电性特征，利用已知的可控(可测)电磁信号或电流在勘探目标所在介质区域产生电(磁)场，通过测量勘探目标的对激励电(磁)场的响应，分析和了解勘探目标的导电(导磁、极化)性质.矿山开采需要了解地层及其富(通)水性、构造破碎带、采空区、塌陷区，或在其深部和周边寻找新的矿产资源量；新建铁路线需要路基勘查；油气勘探中需要通过电法勘探获取油气储层的电性参数以评价储层的含油气性，油田开发压裂中需要检测压裂效果；钻井过程中需要查证漏水漏油程度.这些最需要勘查的地方，往往存在严重的电磁干扰，然而现有的电法(电磁法)、磁法均难以获得令人满意的结果.如何减少甚至克服电磁干扰，提高电(磁)法仪器的抗干扰能力已成为电法勘探方法亟待解决的重要技术难题.
地球物理仪器的研制目前越来越受到我国学界的重视和国家层面的关怀，并开始取得成效(；；；；；；).伪随机编码技术是广泛应用于航天遥测和现代通讯领域中的一项抗电磁干扰技术().国内外地球物理界一些有识之士早就意识到该技术在提高地球物理电法仪器抗干扰能力方面的应用潜力().何继善院士有专著论述，并推出了伪随机信号激电仪()；谢丽等人也对扩频、相关检测、解扩技术在瞬变电磁法中的应用进行过研究与探讨(胡文宝等，年Larose等人对国际上伪随机波场相关技术的地球物理中的应用进行了综述报道( ).但现代通信中的抗干扰技术绝非仅靠某种单一的方法，而需要多方面技术的综合应用( ).本文重点讨论伪随机码扩频发送、相关检测、过采样等技术在抗干扰电法仪的硬件系统和处理算法中的实现及压制电磁干扰的应用效果.
相关叠加抑噪提高信噪比原理示意图
Fig. 1 Schematic diagram of noise depression by correlation superimposing
1 编码电法仪的抗干扰原理
编码电法仪设计的基本思路是：向地下发射伪随机编码的脉冲信号(；)，利用相关叠加抑噪技术()压制干扰，实现微弱信号的检测和参数提取；采用过采样技术降低系统的量化噪声，提高数据采集的精度.
1.1 相关叠加抑噪
假定勘探目标及其所在介质区域为线性时不变系统，则电法勘探采集数据的数学过程可表示为
其中：y(t)是仪器获得的测量信号，即接收电极MN之间的测量信号；G(t－σ)是大地的响应函数；u(σ)是激励源信号(伪随机时间序列)；e(t)是干扰信号.
将仪器获得的测量信号y(t)与激励源信号u(σ)进行褶积，得到
由于使用的激励源信号是伪随机时间序列，为脉冲函数，因此，上式又可表示为
式中： G(T)为响应函数的信号；n为伪随机时间序列的项次数.由于干扰信号与伪随机码的相关性很小，所以可以近似成
式中： 为各项伪随机信号激励对应的噪声.这样当n足够大时，有用信号与噪声的放大倍数之比为：
(2n-1)/2n，所以通过改变n获得足够信噪比特性，从而获得更接近于真值的大地响应函数.
1.2 过采样技术
当信号的采样频率高于奈奎 斯特(Nyquist)采样频率时，称为过采样.除应用相关叠加抑噪技术提升信噪比外，编码电法仪还采用了过采样技术，过采样技术通过数字平均来减小折合到输入端的噪声，提高信噪比，从而提高分辨率，提高数据采集的精度(；).
在一个过采样率为k的过采样系统中，采样将通过过采样率k来抽取.经公式(5)运算，提高后的效果见表 1.采样率每增加1倍，通过过采样可以增加信噪比3 db.
式中：SNR为信躁比，N一个采样周期，db为信躁比单位：分贝.
表1(Table 1)
表 1 过采样方法提高信噪比的效果
Table 1 The effect of improving signal-to-noise by over sampling
过采样率系数信噪比提高扩展的分辨率
表 1 过采样方法提高信噪比的效果
Table 1 The effect of improving signal-to-noise by over sampling
2 MATLAB系统仿真结果
当前的电法勘探仪器在噪声功率为信号功率2倍的前提下，几乎不可能进行有效观测.而在进行仿真设计时，采用噪声功率为信号功率2倍的信号源进行信号源模拟(；；).
图 2为利用MATLAB软件在无干扰情况下对于伪随机信号进行相关检测的仿真结果，信噪比提高明显，快速傅里叶变换(FFT)后得到的幅频曲线和相频曲线都呈现很明显的周期变化，因此可以作为标准曲线来进行对比分析.
无干扰情况下对于伪随机信号进行
相关检测的仿真结果
(a)FFT后幅频曲线；(b)进行相关检测后的信号；
(c)FFT后相频曲线；(d)无干扰情况下发出的信号.
Fig. 2 The simulation results of pseudo-r and om signal
by related detection method with no interference
(a)Amplitude frequency curve after FFT;(b)The signal after
related (c)Phase frequency curve
after FFT;(d)The signal with no interference.
图 3是在强噪声(噪声强度为信号强度的2倍)情况下的仿真结果，观测波形被严重干扰，进行相关检测后信号被检出，但是噪声水平很高.在强干扰条件下仅仅采用相关检测技术不能完全解决干扰问题.
强干扰(噪声强度是信号强度的2倍)情况下
对于伪随机信号进行相关检测的仿真结果
(a)FFT后幅频曲线；(b)相关检测后的信号；(c)FFT后
相频曲线；(d)强干扰信号与原信号叠加情况下的信号.
Fig. 3 The simulation results of pseudo-r and om signal
by related detection method with strong interference
(The noise intensity is two times as the signal intensity)
(a)Amplitude frequency curve after FFT;(b)The signal after
related (c)Phase frequency curve after
FFT;(d)The signal produced by the overlay of strong
interference signal and the original signal.
图 4(过采样系数k=100)和图 5(过采样系数k=1000)显示的是对于伪随机信号进行过采样后，再进行相关检测的仿真结果.可以看出压制干扰信号的效果明显，FFT后得到的幅频曲线和相频曲线，呈现很明显的周期变化.与无干扰情况下标准曲线进行对比发现：曲线形态一致，还原了信号本象，抗干扰效果明显.
强干扰(噪声强度是信号强度的2倍)情况下
对于伪随机信号进行过采样(过采样率为100)和
相关检测后的仿真结果
(a)FFT后幅频曲线；(b)过采样(过采样率为100)和
相关检测后的信号；(c)FFT后相频曲线；
(d)强干扰信号与原信号叠加情况下的信号.
Fig. 4 The simulation results of pseudo-r and om signal
after over sampling(over sampling rate is 100) and
related detection with strong interference(The noise
intensity is two times as the signal intensity)
(a)Amplitude frequency curve after FFT;(b)The signal
after over sampling(over sampling rate is 100) and related
(c)Phase frequency curve after FFT;
(d)The signal produced by the overlay of strong
interference signal and the original signal.
强干扰(噪声强度是信号强度的2倍)情况下
对于伪随机信号进行过采样(过采样率为1000)
和相关检测后的仿真结果
(a)FFT后幅频曲线；(b)过采样(过采样率为1000)和
相关检测后的信号；(c)FFT后相频曲线；
(d)强干扰信号与原信号叠加情况下的信号.
Fig. 5 The simulation results of pseudo-r and om
signal after over sampling(over sampling rate is 1000)
and related detection with strong interference
(The noise intensity is two times as the signal intensity)
(a)Amplitude frequency curve after FFT;(b)The signal
after over sampling(over sampling rate is 1000) and related
(c)Phase frequency curve after FFT;
(d)The signal produced by the overlay of strong
interference signal and the original signal.
综上仿真结果可以看出，抗干扰编码电法仪通过采用伪随机编码信号的相关叠加，压制带内干扰，同时利用过采样技术压制带外干扰，二者结合使用，能够有效压制强干扰.
3 系统的实现
编码电法仪通过发射机向地下发送具有伪随机码特性的电流波形作为激励信号，接收机接收由激励信号产生的电压信号.利用发送的伪随机编码电流序列对接收信号做互相关处理得到大地响应函数，然后对发送的激励信号和提取的大地响应函数做快速傅立叶变换(FFT)，分别得到其幅频和相频特征曲线.用大地响应的相频特征曲线的各频点值减去发送的激励信号的相频特征曲线相应频点值，得到各个频点的绝对相位值φs(；)，用大地响应的幅频特征曲线各频点值除以发送的激励信号的幅频特征曲线对应频点的值，再乘以装置系数K即得到视电阻率值为
在硬件实现上，电流波形伪随机编码的码长可调、本源多项式可选择、码速率可调、起始相位可调，所发送的激励信号是从初相开始，按照设定的码长、本源多项式、码速率和起始相位将具有伪随机特性的电流波形依次发送，以2n-1为周期，n是一个大于0小于20的整数.
3.1 发射机设计
伪随机编码电法仪发送机由以下8个部分组成，如图 6所示：大规模可编程逻辑器，时钟同步电路，逻辑时序发生器，自采样处理电路，直流电源，驱动电路及保护电路，控制计算机.
发射机结构组成框图
Fig. 6 Block diagram of transmitter structure
大规模可编程逻辑芯片是发射机的数学运算和逻辑控制的重要单元，它要完成伪随机码的形成、供电数据的采样、滤波、相关计算、时钟同步、数据后处理计算等功能.其主要功能框图如图 7所示.
发射机大规模可编程逻辑芯片功能框图
Fig. 7 Functional block diagram of the large-scale programmable logic chips of transmitter
电流信号通过AD采样电路变成数字信号，数字信号通过零点漂移校正，温度系数校正，卷积处理自相关形成脉冲过渡函数.脉冲过渡函数再经过傅立叶变换形成电流信号的幅频和相频曲线.
同步信号是GPS模块输出的秒脉冲信号，秒脉冲信号经过本地的数字锁相环把本地时钟同GPS时钟同步.
发射机的主要功能如下：
(1)发射机发送的伪随机码波形可以按照不同的码长、本源多项式、码速率和初相任意选择确定，可以针对不同地区的干扰情况选择合适的伪随机码波形取得满意的信噪比，从而能很好压制生产矿山工业电磁干扰，可在生产矿山深边部开展勘探.
(2)同一地区工作的多台测量装置可以选择不同的码长和本源多项式工作，解决了电法勘探仪器互为干扰源从而不能同时工作的问题.
(3)发射机和接收机精确同步，保证了接收机和发射机起始相位相同，可实现多个频率的振幅、绝对相位、相对相位谱精确测量，可根据所测得的谱特征准确判别地下异常源性质，进行找矿.
(4)测量装置的抗干扰设计为降低主动源的功率提供了可能，实现了轻便化和可一机发送多机接收工作方式，能满足快速大面积勘查，特别是降低了在西部荒漠和高山地区的劳动强度，同时提高了工作效率.
3.2 接收机设计
接收机的原理框图如图 8所示，它主要由滤波及保护电路、放大及AD采样电路、数据处理及滤波电路、时钟同步电路、数据显示及控制计算机组成.
接收机结构组成框图
Fig. 8 Block diagram of receiver structure
接收机大规模可编程逻辑芯片功能框图
Fig. 9 Functional block diagram of the large-scale programmable logic chips of receiver
接收机的主要功能如下：
(1)滤波及保护电路主要是完成信号的滤波及内部电路的保护功能，滤波电路是由高通滤波器、低通滤波器、陷波滤波器综合构成的，保护电路主要是钳制输入端的信号不大于输入电路的最大电压值以保护内部电路的安全和可靠运行.
(2)放大及AD采样电路的主要是放大弱小输入信号、完成模数转换的功能.这部分电路是这块板的模拟电路部分的核心，它要求低噪声、高带宽等，以完成前端弱小信号的采样，为后续处理提供坚实的基础.
(3)数据处理及滤波电路主要功能是以FPGA及周边电路组成，它要完成数据同步采样、线性滤波、数字抗混滤波、低通滤波、采样控制、模式控制、通讯等功能，是接收机的核心.
(4)时钟同步电路主要是完成同发射机的时钟同步功能，同步可以使用GPS同步也可以使用线同步方式，时钟同步电路的精度决定了接收机的激发极化参数的精度，所以这部分电路的设计也十分重要.
(5)数据显示及储存计算机的功能是显示、存储、控制接收机硬件电路的模式选择、参数选择、采样数据显示、输出数据的显示及存储.
电流信号通过AD采样电路变成数字信号，数字信号通过零点漂移校正，温度系数校正，通过卷积处理自相关形成脉冲过渡函数.脉冲过渡函数再经过傅立叶变换形成电流信号的幅频和相频曲线.
同步信号是GPS模块或发射及输出的秒脉冲信号，秒脉冲信号经过本地的数字锁相环把本地时钟同GPS时钟或发射机同步.
3.3 阻容网络验证
为了检验仪器测量参数的正确性，根据大地响应特征而设计的阻容网络试验电路，如图 10所示.
图 10图 10
阻容网络试验电路图
Fig. 10 Test circuit diagram of resistance and capacitance network
图 10中R1是发射机内部的采样电阻，R2是模拟接收极MN之间的电阻，C是模拟MN之间的电容.选定电容C，根据上图计算出MN两点间各个频点的幅值和相位，与编码电法仪实际测的参数对比，用以验证仪器定义电阻率、频散率、相位的正确性.根据电路图，解析解推导公式为
S是我们需要的脉冲响应函数.图 11和图 12是编码电法仪所测结果与理论解析解的比较图，其中电容选择了1000 μF、2200 μF、4700 μF、10000 μF四种情况.通过电路试验，说明了编码电法仪在原理电路上测量的参数和理论计算的参数是一致的.
图 11图 11
幅频曲线实测值与理论值对比图
Fig. 11 Comparison diagram of measured and theoretical
values of amplitude frequency curve
图 12图 12
相频曲线实测值与理论值对比图
Fig. 12 Comparison diagram of measured and theoretical values of phase frequency curve
4 应用效果分析
4.1 基础调查—内蒙古四子王旗电阻率标准场建立
试验区行政区划属内蒙古自治区四子王旗，区内无城市或县城，距离达尔汗茂明安联合旗8 km、距离四子王旗20 km，试验区内有国道通过，另外还有较多的苏木级别公路相通，交通比较便利.本项目共计设计电测深点945个，要求在2月完成，时间紧任务重，常规仪器在同一工区工作互相影响，一个台班一天6个点算，需要半年完成.布设抗干扰仪器可以上多个台班选择不同的多项式发射伪随机码互相之间不影响.该仪器的投产，能适应大规模、会战式工作开展，规模效益明显.
图 13为AB=920 m时的视电阻率平面等值线图(1:10万)，可以看出工区电性可以分为三个区：西北部电阻率低值区，中部电阻率高值区，南部电阻率高值区，中部高值区和南部高值区被断裂破碎带分隔，与地质构造图 14相吻合.
图 13图 13
AB=920 m时电阻率平面等值线图
Fig. 13 Resistivity plane contour map(AB=920 m)
图 14图 14
地质构造图
Fig. 14 Geology structural map
4.2 城市活断裂勘查—秦皇岛活断裂勘查
金梦海湾位于秦皇岛市海港区西部滨海地带，东起闻涛园西侧排洪河，西至滨海大道山东堡立交桥，北至京山铁路港口连接线，南至渤海潮汐线，工作区干扰严重，施工现场存在大型机械，火车高压线，居民生活用电、潮汐等强干扰.北京地震局布设浅层地震勘探和常规电测深工作，均由于干扰未能完成.编码电法仪采用电阻率对称四极测深装置，查明了断裂位置.
图 15图 15
金梦海湾10线综合剖面图
Fig. 15 Comprehensive cross-section of the tenth line in gold dream gulf
5.1&&&&抗干扰编码电法仪通过伪随机编码信号的相关检测理论压制带内干扰，利用过采样技术压制带外干扰，二者合一，理论上证明其能有效压制强噪声干扰.
5.2&&&&完成了仪器的硬件和软件系统的设计，通过阻容网络实验与理论计算结果对比，验证了仪器的抗干扰性与所取得参数的正确性.
5.3&&&&该技术使得仪器体积大大减小，提高了仪器的便携性能，为野外高效施工提供了便利.
5.4&&&&几个具体工程实例检验了仪器的抗干扰特性.抗干扰编码电法仪能够压制电磁干扰，可检测微弱的有用信号，增加了勘探深度；且无需大功率、大电流发射，仪器更加简便，大幅度减少仪器制作及工程探测施工成本.
5.5&&&&抗干扰编码电法仪投入工程应用，必将产生显著的经济与社会效益.
Barbier M G, Viallix J R. 1973. Sosie: A new tool for marine seismology[J].
Cha G M, Xiong X Z. 1990. Spread Spectrum Communication (in Chinese)[M]. Xi’an: Xi’an Electronic and Science University Press.
Duncan P M, Hwang A, Edwards R N, et al. 1980. The development and applications of a wide band electromagnetic sounding system using a pseudo-noise source[J].
Guo H Y, Bi H J. 2004. Application of MATLAB in generation and simulation of PN code [J]. Computer Simulation (in Chinese), 21(3): 149-152.
He J S. 2010. Wide Field Electromagnetic Sounding and Pseudo-random Electrical Prospecting Method (in Chinese)[M]. Beijing: Higher Education Press.
Larose E, Margerin L, Derode A, et al. 2006. Correlation of random wavefields: An interdisciplinary review[J].
Li Y Y, Chen J, Zhang Y M, et al. 2011. Method for SNR estimation based on autocorrelation function [J]. Journal of Military Communications Technology (in Chinese), 32(2): 18-21.
Li Z L. 2006. Modern wireless and mobile communication technology (in Chinese)[M]. Beijing: Science Press.
Liu P,Wang Y F,Yang M M,et al.2013.Seismic data decomposition using sparse Gaussian beams. Chinese J.Geophys..
Liu X J, Xu Y Y, Ji C H. 2003. Application of Phase-coding Signal in LPI Radar [J]. Modern Radar (in Chinese), 25(7): 11-13.
Luo X Z, Zhao B R, Hu W B. 2011. KGR-1 anti-jamming power law instrument development results (in Chinese)[C]. The 27th Chinese Geophysical Society Annual Meeting Changsha, Hunan.
Ma L, Lu H Y. 2009. The application of oversampling techniques in DSP [J]. Journal of Shenyang Institute of Engineering: Natural Science (in Chinese), 5(1): 81-83.
Wan L,Liu T T,Lin J,et al.2013.Joint inversion of MRS and TEM data based on adaptive genetic algorithm. Chinese J.Geophys.(in Chinese),56(11):.
Wang Q,Liu J,Yu S B,et al.2013.Study on influence and correction of coil attitude and bird swing for the fixed-wing time-domain electromagnetic system. Chinese), Chinese J.Geophys.(in Chinese),56(11):.
Wang Y D, Wang J X. 2010. MATLAB-based simulation on frequency-hopping spread aspectrum communication system [J].
Wu M J, Du T C. 2002. The pseudonoise sequence code and the birth for computer [J].
Xiao L Z,Xie Q M,Xie R H,et al.2013.Noise reduction for NMR logging with regularization-heursure algorithm. Chinese), Chinese J.Geophys..
Xie L, Cao Y, Xiang F, et al. 2007. Spread spectrum technology in the application of transient electromagnetic sounding (in Chinese)[C]. Jingzhou: The 8th Chinese International Earth’s Electromagnetic Science Seminar.
Zhang Q S,Deng M,Liu N,et al.2013.Development of the expendable current profiler.Chinese J.Geophys.(in Chinese),56(11):.
Zeng X W, Liu N A, Sun X P. 2004. Extended Spectrum Communications and Multiple Access Technology (in Chinese)[M]. Xi’an: Xi’an Electronic and Science University Press.
Zhao C L, Zhang C. 2008. To improve the measure dynamic range by the sampling technology and FFT algorithm [J]. Science and Technology Information (in Chinese), (17): 424-425.
Zhu W H,Di Q Y,Liu L S,et al?.2013.Development of search coil magnetometer based on magnetic flux negative feedback structure. Chinese J.Geophys.(in Chinese),56(11):.
查光明, 熊贤祚. 1990. 扩频通信[M]. 西安: 西安电子科技大学出版社.
郭海燕, 毕红军. 2004. MATLAB 在伪随机码的生成及仿真中的应用[J].
何继善. 2010. 广域电磁法和伪随机信号电法[M]. 北京: 高等教育出版.
李银永, 陈瑾, 张玉明, 等. 2011. 基于自相关函数的信噪比估计方法[J].
李仲令. 2006. 现代无线与移动通信技术[M]. 北京: 科学出版社.
刘鹏,王彦飞,杨明名,等,2013.地震数据的稀疏高斯束分解方法..
刘晓娟, 徐永元, 姬长华. 2003. 相位编码信号在LPI雷达中的应用[J].
罗先中, 赵壁如, 胡文宝. 2011. KGR-1抗干扰电法仪器研制成果[C]. 中国地球物理学会第二十七届年会..
马玲, 卢海英. 2009. 过采样技术在DSP中的应用[J]. 沈阳工程学院学报: 自然科学版, 5(1): 81-83.
万玲,林婷婷,林君,等,2013.基于自适应遗传算法的MRS-TEM联合反演方法研究.地球物理学报,56(11):.
王琦,林君,于生宝,等,2013.固定翼航空电磁系统的线圈姿态及吊舱摆动影响研究与校正..
王玉德, 王金新. 2010. 基于MATLAB的跳频扩频通信系统的仿真研究[J].
吴明捷, 杜天苍. 2002. 伪随机码及计算机的产生[J].
肖立志,谢庆明,谢然红,等,2013.核磁共振测井的正则化-启发式阈值降噪研究.地球物理学报,56(11):.
谢丽, 曹杨, 向峰, 等. 2007. 扩频技术在瞬变电磁测深中的应用初探[C]. 荆州: 第8届中国国际地球电磁学讨论会.
曾兴雯, 刘乃安, 孙献璞. 2004. 扩展频谱通信及其多址技术[M]. 西安: 西安电子科技大学出版社.
张启升,邓明,刘宁,等,2013.投弃式海流电场剖面仪研制.地球物理学报,56(11):.
赵成岭, 张超. 2008. 采用过采样技术和FFT算法提高测量动态范围[J].
朱万华,底青云,刘雷松,等.2013.基于磁通负反馈结构的高灵敏度感应式磁场传感器研制..