用户名: 密码: 验证码:
强干扰环境下电磁探测技术研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
在矿区周围开展“就矿找矿”是解决目前能源危机的有效办法。利用电磁探测方法进行的“就矿找矿”工作,由于受地形因素及各种人文强干扰影响,对仪器设备的要求较高,且数据受多种环境噪声影响质量较差,需要进行相关处理后方可使用。
     本文从仪器系统的角度出发,介绍了电磁探测采用的仪器设备,包括发射系统、接收系统的软硬件结构。主要对发射天线的理论模型进行仿真验证,得到不同条件下的电阻率值与收发距、工作频率的关系,为野外实验装置的布局提供理论指导。为满足接收机操作要求,设计了基于LabVIEW的软件操作界面,完成仪器控制、显示、信号预处理等功能。磁传感器方面,对于影响传感器动态灵敏度的参数进行分析,并提出相应的参数修正方案。为了得到更精确的传感器实测灵敏度曲线,在分析了标定磁场分布的不均匀性之后,采用均值法对其进行处理,提高了标定精度。针对电磁数据的噪声干扰问题,提出采用FastICA算法对采集数据进行处理。为验证该方法的有效性,在分析强干扰探测环境如矿区的噪声特点之后,模拟仿真各种噪声,建立含噪电磁数据模型,作为FastICA的处理对象,并建立了相应的数据质量评价标准,从量化的角度评价去噪效果。分别对仿真数据模型和实测数据进行去噪处理,处理后数据的信噪比有明显改善,达到去噪目的。
“Ore-prospecting nearby the existing deposits” is an effective wayto solve energy crisis. However, the working environment of mine areais complex, so it puts forwards higher request to the detection equipmentand the data processing method. From the view of ground and deviceplaced, a loop antenna which can be approximated to a horizontalmagnetic dipole has been designed and it has been demonstrated intheory; In order to obtain a higher sensitivity of magnetic sensor, theparameter analysis and correction has been done; In the last, for theelectromagnetic noise question, the FastICA algorithm has been putforward to separate signal from data with noise, and then thehigher-quality electromagnetic data will be gotten.
     The main contents of this paper can be summarized as followings:
     1) The theories argument of transmitting antenna used inelectromagnetic detection has been done. In order to make up for the lack of natural magnetic field signal, the high-frequency electromagneticwave has been produced by artificial method. The transmitting device isdesigned as a couple of loop antennas which are mutually perpendicular,and they are perpendicular to the ground while working. The signalswith14different frequency components which is between1kHz and80kHz were generated by selecting different combination of capacitors. Insome certain condition, the transmitting antenna can be approximated toa horizontal magnetic dipole. Through theoretical deduction andexperimental simulation, the electromagnetic wave generated by thehorizontal magnetic dipole is deemed to satisfy the plane wave in somesituations. The simulation result of resistivity under different conditionsgives the the relationship between distance, frequency and targetresistivity, and it provides a theoretic guide for the device placed. Thereceiver surface has been designed using LabVIEW graphicprogramming platform, and it provides a friendly man-machine interface.Using the interface, the experimental data can be observed at any time,and the function such as data pre-processing can be realized togetherwith some simple algorithms.
     2) The non-uniformity of calibration magnetic field which is used inthe sensor sensitivity test is demonstrated in theory, and thenon-uniformity has been processed using average method. Theenvironment for magnetic sensor calibration requires no electromagneticinterference, and the calibration magnetic field must be stationary anduniform. However, the ideal situation cannot be achieved in practice.This paper analyzed the distribution of the magnetic field generated bydifferent equipment and decided to use long straight solenoid as themagnetic field generating equipment. The magnetic field generated bythe long straight solenoid can be deemed as stationary and uniform. Inthe past calibration experiments, the inner magnetic field of the solenoidis regarded as uniformly distributed, and the magnetic field around itscentral axis was used as the estimation of its inner magnetic field. Thisestimation caused the sensitivity of magnetic sensor to be lower than itsactual value so that the deviation is introduced. There is no numericalsolution of the non-axis magnetic formula, so the high-precisionGaussian quadrature formula is proposed and the result is consideredaccurate. The Simulation results show that the inner magnetic field of solenoid has little variation in landscape orientation. However, in thevertical orientation, the magnetic field is homogenous in the middle butvaries a lot in both ends. The intensity of magnetic field at both endsdrops to50%of that in the middle. In order to reduce the impact of thenon-uniformity and increase the accuracy, the average method isproposed to deal with this question, and the length of step is minimizedas possible to improve the accuracy. The experimental sensitivityprocessed by this method has improved than before.
     3) Various factors influencing sensor theoretical sensitivity has beenanalysed in detail, and the impact of the magnetic core’s averagemagnetic permeability is mainly discussed. Moreover, a correctionprogram for parameters in theoretical sensitivity formula of magneticsensor is also put forward. Since the magnetic core of magnetic sensorwas made from permalloy, and permalloy has a highmagnetic permeability under low frequency but themagnetic permeability drops under high frequency due to the increase inmagnetic hysteresis loss. In the broadband measurment, the effective magnetic permeability should be a function to the frequencybut not a fixed value. Unfortunately, the magnetic permeability ofmagnetic core material has a complicated relationship to the frequency,and it is also influenced by its shape. Due to the complexity, there is noquantitative function to calculate the relationship between them.Through the comparison with the experiment, the functionμ app (f)=105
     0.005f2+0.0005f+105μapp低is put forward to calculate therelationship between the magnetic permeability of magnetic core and itsoperating frequency. The compare between theoretical sensitivity andexperimental sensitivity showed that the corrected theoretical sensitivityis coincide with the experimental sensitivity.
     4) After analyzing the feature of the electromagnetic signal andnoise, the noisy electromagnetic data model has been built. The workingenvironment of mine area is complicated and there are various noises.The correct method to deal with the noise and gain high quality sensingdata were adopt only when the feature and classification of the noisewere recognized. Some common noise source in mine area has been analyzed, and time-domain waveform, bandwidth and impact on datainterpretation of the noise has also been introduced in detailsrespectively. The simulated data of different noise is produced after theanalysis. The simulated noise data were added to the original signal sothat the target wave for verifying different noise-cancellation methodwas generated.
     5) The FastICA algorithm is used to process the simulated data andthe experimental data of the electromagnetic detection. After analyzingthe feature of signal and noise of mine area, it was decided to adoptFastICA algorithm to process the data with noise. In order to evaluatethe denoising effectiveness of the algorithm objectively, two qualityevaluate parameters such as SNR and coherency were introduced. Firstly,the noisy electromagnetic data model was processed using FastICAalgorithm, and we found from the comparison of wave in bothtime-domain and frequency-domain before and after denoising that thesignal can be clearly distinguish in the denoised data while it is totallysubmerged in the data before denoising. The SNR of processed signal was much higher than the original data, so that it is proved this methodis very practical in processing the simulated data of mine area. In orderto test the method’s application effectivity for the experimental data, theFastICA algorithm is used to denoise the received pseudo-randomthree-frequency signal and the pseudo-random five-frequency signal.The processing result shows that this method is capable for extractingmeaningful signal from noisey data, and the coherency of electric signaland magnetic signal is improved obviously.
     The main innovation works of this paper is as followings:
     1) The circuits model and physical model of the loop antenna whichis perpendicular have been built. The spatial distribution ofelectromagnetic field produced by antenna under the different geologicalmodel has been gained while the antenna is approximated to a horizontalmagnetic dipole. Combined with the plane wave detection theory, therelationship between the resistivity of conductor and the distance as wellas the working frequency was obtained, which provided a theoreticalfoundation for antenna design.
     2) The average method for processing the calibration magnetic fieldof long straight solenoid is introduced, so the precision of inductivemagnetic sensor’s sensitivity has been improved. Before this paper,almost every magnetic sensor was calibrated with the inner magneticfield of long straight solenoid deeming as homogenous distributed. Andthe intensity of solenoid’s central axis was deemed as the intensity of thewhole magnetic field. This approximation causes the experimentalsensitivity to be lower than its actual value. This paper took thenon-homogenization into consideration and proved the existence of thisnon-homogenization through theoretical deduction and numericalcalculation. Using averaging method to deal with it made the resultsmore close to the actual value.
     3) A function model of the relationship between themagnetic permeability of magnetic core and its working frequency wasbuilt, and this model resolved the problem that there is no formulabetween the magnetic permeability of magnetic core and its workingfrequency, and then the theoretical model of the magnetic sensor has been corrected. The effective magnetic permeability of magnetic corewhich is made from magnetic materials is not only related to the shapeof core but is also affected by the working frequency. Theeffective magnetic permeability of magnetic cores made from permalloymaterial is high while working under the low frequency, and it dropsalong with the increase of working frequency. So far, there is no specificformula to characterize the law of the change. This paper reviewedrelated documents and considered the comparison between magneticsensor’s actual measured sensitivity and its theoretical sensitivity. Thispaper refered to related documents and considered the comparisonbetween experimental sensitivity and its theoretical sensitivity. And thena low-pass filter’s transfer function was utilized as the initial mathematicmodel for effective magnetic permeability of magnetic core and itsworking frequency. The parameters of the initial mathematic model weremodified to achieve a more precise function expression after manyexperiments. The result of the research has modified the formulation tothe magnetic sensor’s sensitivity.
     4) The noisy electromagnetic data model of mine area has beenbuilt, and an electromagnetic processing method based on FastICAalgorithm for mine area detection is so also proposed.that a standardelectromagnetic database of mine area for signal processing is provided.Using this method, the quality of data has improved as well as theutilization of the data has increased. The working environment of minearea is complicated and there are various noises. It will bring moreconvenience for the electromagnetic data processing of mine area tobuild a standard noisy electromagnetic data model of mine area afteranalyzing the features of the noise and simulating the noise wave.
     The production of this paper has provided theory foundation for the“Ore-prospecting nearby the existing deposits” of mine areaelectromagnetic detection, and and it also provided a reference forinstrument design.
引文
[1]张立公.新一轮金属矿的勘查与开发[J].资源调查与环境,2002,23(3):196-199.
    [2]傅雷,仲冰.中国矿产资源现状与思考[J].资源与产业,2008,10(1):83-86.
    [3]刘江天,刘荫椿,付乐.就矿找矿理论浅析[J].黄金,2009,30(2):10-13.
    [4]郭明凤.浅析如何就矿找矿[J].科学之友,2012,2:36-37.
    [5]吴大忱.“就矿找矿”论[J].地质学刊.,2009,33(2):222-224.
    [6]曹晓明,周贤旭,钟浩.“就矿找矿”的认识与实践[J].东华理工大学学报(自然科学版),2011,34(1):51-56.
    [7]余志庆.谈“就矿找矿”问题[J].中国地质,1990,11:27.
    [8] A.A.考夫曼, G.V.凯勒.大地电磁测深法[M].北京:地震出版社,1987.
    [9] Séran, H.C, P. Fergeaua. An optimized low-frequency three-axissearch coil magnetometer for space research[J]. Review of ScientificInstruments,2005(76): p.1-10.
    [10]曲放,翁爱华.瞬变电磁法及其应用研究[J].光机电信息,2011,28(12):82-88
    [11]裴婧.广域电磁仪接收机硬件设计与实现[D].中南大学,2010.
    [12] U.Matzander. ADU-07Operating Manual. Bransschweiger,Germany: Metronix,2008.
    [13]林品荣,赵子言.分布式被动源电磁法系统及其应用[J]地震地质,2001,23(002):138-142.
    [14]魏文博.我国大地电磁测深新进展及瞻望[J].地球物理学进展,2002,17(002):245-254.
    [15]程德福,王君,李秀平等.混场源电磁法仪器研制进展[J].地球物理学进展,2004,19(4):778-781.
    [16]邓明,魏文博.海底天然大地电磁场的探测[J].测控技术,2003,22(1):5-8.
    [17]陆其鹄,孙进忠.发展技术,增强基础-地球物理仪器国际化与地球物理技术在工程上的应用研讨会综述[J].地球物理学进展,2004,19(4):726-729.
    [18]陆其鹄,彭克中,易碧金.我国地球物理仪器的发展[J].地球物理学进展,2007,22(4):1332-1337.
    [19]滕吉文.中国地球物理仪器和实验设备研究与研制的发展与导向[J].地球物理学进展,2005,20(2):276-281.
    [20]滕吉文.中国地球物理仪器的研制和产业化评述[J].石油物探,2006,45(3):209-216.
    [21]陈儒军.伪随机多频电磁法观测系统研究[D].中南大学,2003.
    [22]李文彬.大地电磁测深系统磁传感器的研究与实现[D].中国舰船研究所,2012.
    [23]章克来,朱海明.微弱信号检测技术[J].航空电子技术.2009,40(2):30-36.
    [24]蒋洪波,陈超,秦其明.天然源超低频频谱的曲波分解与分析[J].光谱学与光谱分析,2012,32(2):472-475.
    [25]方晓梅,昌彦君.自适应滤波器用于瞬变电磁法工频干扰压制研究[J].中国地球物理学会第十八届年会,2002:164.
    [26]车宏伟.基于主成分的时间域航空电磁数据神经网络反演方法研究[J].吉林大学,2012.
    [27] Zhu Kai-Guang, Ma Ming-Yao, Che Hong-Wei,etc. PC-based artificial neural network inversion for airborne time-domainelectromagnetic data[J]. APPLIED GEOPHYSICS,2012,9(1):1-8.
    [28]李培军,金慧然,宋本钦.基于独立成分分析的超低频电磁探测信号滤波[J].北京大学学报(自然科学版),2009,45(4):621-627.
    [29]陈乐寿.大地电磁测深资料处理与解释[M].石油工业出版社,1989.
    [30] Louis Cagniard. Basic theory of the magneto-telluric method ofgeophysical prospecting[J]. Geophysics,1953,18(3):605-635.
    [31] William H. Hayt and Jr. John A. Buck. Engineeringelectromagnetics(Sixth Edition)[M]. McGraw-Hill,2001.
    [32] David J. Griffiths. Introduction to electrodynamics(Third Edition)
    [M]. Person Education,2005.
    [33]陈重,崔正勤.电磁场理论基础[M].北京:北京理工大学出版社,2003.
    [34]王言章.混场源电磁探测关键技术研究[D].吉林大学,2010.
    [35] Jutten C, Herault J. Independent component analysis versus PCA.In: Proceeding of European Signal Processing Conf,1988.287-314.
    [36]Comon P. Independent component analysis-A new concept. SignalProcessing,1994,36(3):287-314.
    [37] Jut t en C, Her aul t J. Blind separat ion of sources, part I: Anadaptive algorithm based on neu romimetic architecture. S ignalProcessing,1991,24(1):1-10.
    [38] Oja E. The nonlinear PCA learing rule in independent componentanalysis. Neur ocomputing,1997,17(1):25-46.
    [39] Hyvarinen A. Survey on independent component analysis. NeuralComputing Surveys,1999,2(4):94-128.
    [40] Hyvarinen A. Fast and robust fixed-point algorithms f or independent component analysis. IEEETrans. on Neural Networks,1999,10(3):626-634.
    [41]赵峰.独立成分分析改进算法与仿真研究[D].大连交通大学,2010.
    [42]张建明,林亚平,吴宏斌等.独立成分分析的研究进展[J].系统仿真学报,2006,18(4):992-997.
    [43] Hyvarinen A, Oja E. Independent component analysis: a tutorial.Neural Networks,2000,13(4-5):411~430.
    [44]姚志湘.过程系统状态空间统计模型和化工过程监测[D].华南理工大学,2003.
    [45]曾生根.快速独立分量分析方法及其在图像分析中的若干应用研究[D].南京理工大学,2004.
    [46]马建.基于独立成分分析的说话人识别技术研究[D].电子科技大学,2004.
    [47]徐丽琴.典型ICA算法在盲源分离中的实现[J].科技信息,2010,26:519-521.
    [48]毕杨.基于快速独立分离分析的盲源分离算法研究及应用[D].西安理工大学,2007.
    [49]冶继民.盲信号分离fast-ICA算法研究[D].西安电子科技大学,2011.
    [50]吴逊.基于独立成分分析的特征提取方法研究[D].厦门大学,2007.
    [51] Giannakopoulos X, Karhunen J, Oja E. Experimental comparisonof neural ICA algorithms. In: Proc. Int. Conf.on Artif icial NeuralNetworks (ICANN'98), Skvde, Sweden,1998.651~656.
    [52]邹大力.基于计算智能的结构损伤识别研究[D].大连理工大学,2005.
    [53] Kaufman A.A. and Keller G.V. Frequency and transientsoundings[M]. Elsevier Science Publish,1983.
    [54]张秀成.大地电磁测深测量及仪器[M].地质出版社,1989.
    [55]魏文博.我国大地电磁测深新进展及瞻望[J].地球物理学进展,2002,17(2):245-254.
    [56]王家映.我国大地电磁测深研究新进展[J].地球物理学报,1997,40(增刊):206-216.
    [57]刘国栋,陈乐寿.大地电磁测深研究[M].北京:地震出版社,1984.
    [58]万云霞,程德福,卢浩等.混场源电磁探测仪器发射天线设计[J].吉林大学学报(工学版),2012,42(6):1432-1436.
    [59]米萨克N纳比吉安.勘查地球物理[M].电磁法,第一卷,理论,1991.
    [60]吴旭东,侯文生,郑小林等.磁偶极子的定位模型及实验验证[J].仪器仪表学报,2008,29(2):326-329.
    [61]任来平,赵俊生,侯世喜.磁偶极子磁场空间分布模式[J].海洋测绘,2002,22(2):18-21.
    [62] Misac N. Nabighian. Electromagnetic methods in appliedgeophysics[M]. Volume1, Theory,1991.
    [63] Chave A.D. Numerical integration of related hankel transforms byquadrature and continued fraction expansion[J]. Geophysics,1983,48(12):1671-1686.
    [64] F.N. Kong. Hankel transform filters for dipole antenna radiation in aconductive medium[J]. Geophysical Prospecting,2007,55(1):83-89.
    [65]江少恩.基于快速汉克尔变换数值求解波动方程[J].数值计算与计算机应用,1997,3:161-166.
    [66]耿书宇.EM-II型混场源电磁法发射机的研制[D].吉林大学,2008.
    [67]万云霞,程德福,王言章等.混场源电磁接收系统同步数据采集技术[J].吉林大学学报(信息科学版),2009,27(2):162-166.
    [68]王剑.基于LabVIEW的数据采集及分析系统的开发[D].哈尔滨理工大学,2004.
    [69]徐超. LabVIEW在实时测控系统中的应用研究[D].重庆大学,2005.
    [70]孙春龙.基于LabVIEW多通道数据采集分析系统开发[D].武汉大学,2004.
    [71]吕向锋,高洪林,马亮等.基于LabVIEW串口通信的研究[J].国外电子测量技术,2009,28(12):27-31.
    [72]许美玲,李春茂,朱俊峰等.基于LabVIEW串口通讯的多路数据采集系统[J].机械与电子,2010,7(1):172-174.
    [73]王珏,郑建勇.基于LabVIEW的串口调试系统设计[J].江西科学,2007,25(6):755-757.
    [74]陈勖,李尔园. GPS接收机工作原理及发展现状[J].全球定位系统,2010,6:51-55.
    [75]刘娣,薄煜明,赵高鹏. GPS软件接收机信号的快速捕获与跟踪[J].计算机工程,2010,36(1):239-241.
    [76]王尔申. GPS接收机及其应用的关键技术研究[D].大连海事大学,2009
    [77]孙召红.井底钻压扭传感系统及标定装置设计研究[D].中国石油大学,2010.
    [78]蒋安林,周穗华,张晓兵等.感应式磁传感器的设计[J].航电技术,2011.31(12):35-37.
    [79]谢石林,陈胜来,张希农等.传感器标定的神经网络杂交建模方法[J].机械工程学报,2010.46(22):6-15.
    [80]俞富.电路噪声的成因与抑制措施[J].现代电影技术,2006,10:36-38.
    [81] Reza Moghimi.电路噪声分析以及优化[J].电子与电脑,2009,10:43-44.
    [82]邵英秋.宽频带感应式磁传感器的研制[D].吉林大学,2012.
    [83]王岗岭,王坚宾.传感器噪声及其抑制方法[J].传感器世界,2007,2:28-31.
    [84]朱妍妍,段传波,柳至和等.城市轨道交通噪声试验研究[C].中国环境科学学会学术年会论文集,2010:3740-3743.
    [85]谢明,刘湘京,窦燕生.工业噪声环境影响预测方法研究[J].中国卫生工程学,2005,4(3):134-136.
    [86]孙远涛.建筑施工噪声烦恼度阈限值研究[D].长安大学,2008.
    [87] H.C Seran and P. Fergeau. An optimized low-frequency three-axissearch coil magnetometer for space research[J]. Review ofScientific Instruments,2005(76):1-10.
    [88]李久林.电能计量箱柜门禁装置[J].内江电业局.
    [89]岳高翔,马新生.简易磁屏蔽工装的设计[C].第二届民用飞机制造技术及装备高层论坛资料汇编(论文集),2010.
    [90]汝强,胡社军,胡显奇等.电磁屏蔽理论及屏蔽材料的制备[J].包装工程,2004,25(4):21-24.
    [91]赵灵智,胡社军,何琴玉等.电磁屏蔽材料的屏蔽原理与研究现状[J].包装工程,2006,27(2):1-4.
    [92]刘虹让,刘景岩,马全喜.通电圆线圈轴线上磁场分布分析[J].物理通报,,2012,3:19-20.
    [93]李君良,吕玉祥,魏循等.亥姆霍兹线圈的磁场及其测量[J].太原理工大学学报,2005,36(5):564-566.
    [94]刘耀康.导出圆电流的磁感应强度的简便方法[J].大学物理,2007,26(7):32-33.
    [95]刘耀康.螺线管磁场的数值解[J].高师理科学刊,2008,28(1):60-64.
    [96]徐士良.计算机常用方法[M].清华大学出版社,1995.
    [97]郝殿华.不能忘记“就矿找矿”方法[N].中国矿业报,2004-10-30.
    [98]刘国仁.就矿找矿的理论及新疆阿尔泰的找矿实践[J].新疆有色金属,2009,增刊:84-87.
    [99]何玉辉.天然电磁场测深[J].物探化探计算技术,1999,21(2):97-115.
    [100]张木生.大地电磁测深仪器中的几个关键技术[M].地球物理专辑(第二集,电磁测深).地质出版社,1989.
    [101]陈乐寿,王光愕.大地电磁测深法[M].地质出版社,1990
    [102] AA考夫曼,凯勒.大地电磁测深法[M].地质出版社,1987
    [103]荒木忠次,高杉真司,筝原信男.MT法概论[J].地热技术,1982,7(3):1-52.
    [104]蒋洪波,陈超,秦其明.天然源超低频频谱的曲波分解与分析[J].光谱学与光谱分析,2012,,32(2):472-475.
    [105]杨生.大地电磁测深法环境噪声抑制研究及其应用[D].中南大学,2004.
    [106]孙洁,晋光文,白登海等.大地电磁测深资料的噪声干扰[J].物探与化探,2000,24(2):119-127.
    [107]朱威,范翠松,姚大为等.矿集区大地电磁噪声场源分析及噪声特点[J].物探与化探,2011,35(5):659-662.
    [108]范翠松.矿集区强干扰大地电磁噪声特点及去噪方法研究[D].吉林大学,2009.
    [109] Gamble T D,Goubau W M,Clarke J. Error Analysis for RemoteReference magnetotelluric[J]. Geophysics,1979,44.
    [110] Mallat S G. A theory for multire solution singalde composition:the wavele trepresentation[J]. Transactionson Pattern Analysisand Machine Intelligence,1989,11(7).
    [111] Tradand D O,Travassos J M.Wavelet filtering of magnetotblluricdata[J]. Geophysics,2000,65.
    [112]王宏亮.基于中压电力线载波的通信技术研究[D].哈尔滨工业大学,2006.
    [113]李桐林,刘福春,韩英杰等.50万伏超高压输电线的电磁噪声的研究[J].长春科技大学学报,2000,30(1):80-83.
    [114]徐志敏,汤井田,强建科.矿集区大地电磁强干扰类型分析[J].物探与化探,2012,36(2):214-219.
    [115]汤井田,徐志敏,肖晓等.庐枞矿集区大地电磁测深强噪声的影响规律[J].地球物理学报,2012,55(12):4147-4159.
    [116] Foster,M.R.and Hermance,The Coefficient of coherence:Itsestimation and use in Geophsical data processing. Geophysics,1978,32(4):743-760.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700