摘要
由何继善院士发明并命名的伪随机多频信号带宽可控,其能量集中在数个对数坐标上呈等间距分布的频率上,是电法勘探的理想场源。以此为激励信号的伪随机多频激电法勘探,不仅具有时间域激电法的激电异常形态明显、在单次测量中获取丰富信息等特点,而且具有频率域激电法的抗干扰能力强、仪器轻便等优点。
笔者基于对各种激电法充分认识的基础上,分析研究了各种激电参数的作用与特性,并对前人在研究频率域激电法中特别是在研究影响岩矿石激电相位或时间常数的因素上存在的问题,提出了自己的见解。认为比较岩矿石中导电矿物的含量、成分、结构构造对时间常数或绝对相位的影响程度是不合理的,因这种影响因素的划分是重叠的,从而使测试结果不具备可比性。
笔者借鉴时间域激电法对充电率的定义,提出了频率域激电法中的充电率的概念,并对其具体含义作了简单的分析。
基于对双频激电法的认识,笔者借鉴幅频率的定义,提出了频率域激电法的另一个激电参数:相频率,并详细分析了时间常数、充电率、频率相关系数对相频率的影响情况,从而得出相频率在识别不同极化体的时间常数方面有其独特的优势,从而有助于激电异常源性质的区分。
在对激电法的综合分析的基础上,在对幅频率、激电相位、激电相对相位等参数频谱特性的详细剖析上,特别是在频率域激电法充电率与相频率激电参数的研究上,笔者提出了伪随机激电多参数谱法的概念,并对其基本作用与意义作了有益的概括与分析。
基于对双频激电法与相位激电法的综合分析基础上,提出了双频激电相位法,运用该方法既能测量激电异常体的激电效应强弱,还能反映异常源的时间常数大小,进而区分激电异常源性质,而工作量与双频激电法一样。
首次详细研究了伪随机多频激电相对相位法在区分激电异常源性质上的基本原理,为伪随机多频激电相对相位法区分激电异常源提供理论上的支持。
基于对伪随机激电多参数谱法的定义,笔者对均匀介质中的三维立方极化体进行了多参数谱(相位谱、相对相位谱、相频率谱、幅频率谱、电阻率谱)的数值模拟,并分析了各种参数谱对区分激电异常源性质的作用。
针对伪随机激电法原始数据存在着数据格式不统一、电磁耦合、数据拼接、数据插值、人文干扰等问题或要求,为剔除各种人文干扰引起的数据畸变,作者研究了伪随机激电法数据处理。为提高解释系统开发的效率,还探讨了VB与FORTRAN混合语言编程在系统开发方面的应用。为使解释结果更直接、更形象地反映处理结果,笔者在分析各种绘图软件的基础上,对采用VB与OpenGL相结合的方法编写更方便解释系统应用的绘图功能模块作了尝试。
上述对激电法多参数谱、数据处理、绘图功能的研究是为将来对电阻率法、激电法、电磁法等进行联合反演所作的前期研究与尝试,也是将来伪随机多功能电法系统的数据解释系统的一部分。
最后笔者对在研究中未能解决的问题与下一步计划研究的问题作了总结。如COLE-COLE模型的修改问题、混合矿体的区分问题、激电法专家解释系统开发问题、激电法中电极效应校正问题、伪随机多参数谱数据解释研究问题等。
The controlled multi-frequency data bandwidth pseudorandomequipment, invented by and named after academician He Jishan, caninduce an ideal field source of electrical survey, whose energy wascentralized on several equal-spacing dispersed frequency on the logcoordinate. The pseudorandom multi-frequency detection based on theexited signal, not only has the properties of obvious anomaly, as in thetime domain IP method, and rich information in the single survey, butalso has the properties of strong anti-interruption ability as the frequencydomain IP method and portability. Based on the IP method study history,This essay mainly compared the merits and shortcomings among thevariable-frequency method, phase IP method, frequency spectrum methodand dual-frequency IP, all which belong to frequency domain IP.
Throuth the thorough recognization of different IP methods, thevarious IP parameter function and property was studied and the proposalwas promoted on the existed problems in the frequency domain IPmethod, especially about the lithology IP phase and the time constant.The comparison of the conductive mineral content, component andstructure influence degree in the ore on the time constant and absolutephase is unreasonable, because the devision is overlapped.
The frequency domain IP charge ratio was put forward depended onthe definition of time domain IP charge ratio, then the concrete meaningwas analyzed.
Depended on the dual-frequency IP recognition and the definition ofM-frequency, the another IP parameter phase frequency of the frequencydomain IP method was put forward. The influence of time constant,charge ratio and frequency correlation parameter on the phase frequencewas analyzed detailedly. Then the unique advantage on recognition ofdifferent polarization body time constant was concluded, which facilitatedthe property distiction of IP anomaly resouce.
From the analysis on frequency spectrum property of magnitudefrequency, IP phase and IP relative phase, especially the research of thefrequency domain IP method charge ratio and phase-frequency IP parameters, the notion of pseudorandom IP multi-parameters spectrummethod was put forward as well as the fuction and meaning.
In the article, the author advanced and studied the dual-frequency IPphase method, on the basis of dual-frequency IP and phase IP method.The method can measure the IP effect, also reflect the time constant value,so it can distinguish the source of anomaly.
It is the first time to research the keystone of pseudo-randommulti-frequency relative phase method on the division source of anomaly,and the conclusions can supply the method in property distinguish intheory.
In order to analyze the effect of different parameter spectrum on theproperty distiction of IP anomaly resouce, the 3D cubic polarization bodywas simulated numerically and physically by multi-parameter spectrumsuch as: phase spectrum, relative phase spectrum, phase frequencyspectrum and magnitude frequency spectrum, depended on the definitionof pseudorandom IP multi-parameters spectrum method.
Depended on the problems, such as: un-union data format,electromagnetic coupling, data linkage, data interpolation and humaninterruption, the author studyed Pseudorandom IP method datainterpretation system in order to eliminate the mal-data caused bydifferent human interruption. In order to promote the interpretationsystem exploitation effective, the synthesized language program of VBand FORTRAN was discussed. Depend on the analysis of variousdrawing software, some try has done on the drawing module compiled bythe VB and OpenGL method, in order to make the interpretation resultmore direct and visualized.
The study of data processing and drawing is mainly about someex-research and trials on the united inversion of the resistivity method, IPmethod, electromagnetic method and so on for future. Also it will becomethe part of the processing and interpretation system of the pseudorandommulti-function electrical system.
Finally the unsolved and further study problem was concluded. Suchas the COLE-COLE model modification problem, mixture ore-bodydiscrimination problem, IP expert interpretation exploitation problem, electrode effect correction problem and multi-parameter spectruminterpretation problem of the Pseudorandom.
引文
[1] J.S萨姆纳(美)著.陈文华,译.地球物理勘探的激发极化原理.北京:地质出版社.1981
[2] 何继善.双频激电法.北京:高等教育出版社,2005.
[3] 傅良魁主编.激发极化法.北京:地质出版社,1982
[4] 何继善等著.双频道激电法研究.长沙:湖南科学术出版社,1989
[5] 罗延钟,张桂青.频率域激电法原理.北京:地质出版社,1988
[6] 刘崧.谱激电法.武汉:中国地质大学出版社,1998
[7] 慕文斋,龙占奎,刘静.相位激电法中相角测量方法纵横谈.维普资讯
[8] 何继善.伪随机三频电法研究.中国有色金属学报,1994,4(1):1~7
[9] 何继善.2~n系列伪随机信号及应用.见:中国地球物理学会年刊.北京:地震出版社,1998.199
[10] 陈儒军.伪随机多频电磁法观测系统:[博士学位论文].长沙:中南大学,2003
[11] 何继善等.伪随机电磁法及多功能仪器研究报告(内部交流).长沙:中南大学地球物理勘查新技术研究所,1999
[12] 吴汉荣.激电相对相位测量及其地质效果[A].吴汉荣,谢婷婷,等.勘查地球物理勘查地球化学文集第3集—激发极化法专集[C].北京:地质出版社,1985.295~300
[13] 陈儒军,何继善,白宜诚,汤井田.多频激电相对相位谱研究.中南大学学报(自然科学版).2004.35(1):106~111
[14] 张友山,何继善.多频组合波激电法的去耦研究.中南矿冶学院学报,1994,25(5):113—115
[15] 张宪润,陈儒军.激电相对相位法区分矿与非矿异常的成功实例.物探与化探,1998,22(4):251~254
[16] 张友山,何继善.三频激电相对相位观测法.中南矿冶学院学报,1994,25(4):417~421.
[17] 张友山,何继善.伪随机三频波激电法.中南大学学报,1995,26(2):157~161
[18] 何继善.多参数双频激电研究.地球与空间科学观测技术进展.北京:地震出版社,1995,413~421
[19] 陈儒军等.高精度多频电法数据采集系统.物探与化探,2003,27(5):375~378
[20] 何继善,柳建新.伪随机多频相位法应用简介.中国有色金属学报,2002,12(2):374—376
[21] 鲍光淑,何继善.多参数频域激电观测系统的研究.地学仪器,1995,(2):5~ 10
[22] 何继善,白宜诚,浦慧如等.伪随机电磁法及多功能仪器研究科研报告[R].长沙:中南大学地球物理勘查新技术研究所,1999
[23] 何继善.电法勘探的发展和展望.地球物理学报,1997,40(增):308~316
[24] 柳建新.广西泗顶铅锌矿及外围地球物理研究.长沙:中南大学地球物理勘察新技术研究所,2003
[25] 何继善,柳建新.双频激电法矿产普查工作手册.长沙:中南大学,2002
[26] 柳建新,何继善,等.一种区分矿与非矿的有效方法-伪随机多频相位法原理及其应用.中国地质,2001,28(9):41~46
[27] 何继善.提取和利用EM效应和IP效应的方波相干法.见:中国地球物理学会年刊.北京:地震出版社,1994
[28] B.A.柯马罗夫.激发极化法勘探,阎立光等译.北京:地质出版社,1983
[29] 何继善,李大庆,汤井田.频谱激电非线性效应的理论模型.地球物理学报1995.38(5):662~670
[30] I. M., Spectral IP parameters derived from time-domain measurements, in Ward, S. H., Ed., Induced polarization: Soc. Of Expl. Geophys., 1990, 57~78
[31] Lewis, R. J. G. and Bishop, J. R., Target discrimination with time-domain spectral IP inversion, 57th Ann. Internat. Mtg: Soc. Of Expl. Geophys., 198T, Session: MIN1. 6.
[32] Seigel, H. O., Vanhala, H. and Sheard, S. N., Some case histories of source discrimination using time-domain spectral IP. Geophysics, 1997, 62: 1394~1408
[33] Seigel, H. O., Ehrat, R. and Brcic, I., A microprocessor based multichannel time-domain IP receiver, in Ward, S. H., Ed., Induced polarization: Soc. Of Expl. Geophys., 1990, 91~103
[34] Seigel, Harold O., The magnetic induced polarization(MIP) method. Geophysics, 1974, 39: 321~339
[35] 周安昌,杨冠鼎.激发极化相对相位测量在两个矿床上的实验结果[A].吴汉荣,谢婷婷,等.勘查地球物理勘查地球化学文集第3集—激发极化法专集[C].北京:地质出版社,1985:291~294
[36] 何继善.电法勘探的发展与展望.地球物理学报.1997,40:308~315
[37] 张赛珍,等译.地球物理勘探专辑第7辑—电法勘探(激发极化法).北京:中国工业出版社.1965
[38] 张赛珍等.我国几个金属矿区岩(矿)石的低频相位频率特性及影响因素.地理物理学报,1984,27(2)
[39] Pelton, W. H., Ward, S. H. , Hallof, P. G. , Sill, W. R. , and Nelson, P. H. Mineral discrimination and removal of inductive coupling with multi-frequency IP. Geophysics, 1978, 43(3): 588~609
[40] 日丹诺夫M C.电法勘探.潘玉玲,王守坦译.武汉:中国地质大学出版社,1990.397~405
[41] 中国地质大学(武汉)应用地球物理系组编.应用地球物理学进展.武汉:中国地质大学出版社,1996
[42] 王书民.激电相位的一种近似解释方法与应用效果.物探与化探.2004,28(4):317~319
[43] 王书民,雷达.相位激电法(偶极—偶极)单频电磁耦合校正方法.物探与化探.2002,26(1):57~63
[44] 张美玲,王艳,孙宏智.利用相位差与电阻率曲线反演介电常数的方法研究.大庆石油地质与开发.2000,19(3):46~48
[45] 王自力,张赛珍.实测视复电阻率振幅谱换算相位的可行性.物探与化探,1989,13(6):453~457
[46] 罗延钟,吴之训,谱激电中频率相关系数的应用,地球物理学报,1992,35(4)
[47] 张赛珍,王庆乙,罗延钟.中国电法勘探发展概况.地球物理学报,1994,37(增):408~424
[48] 罗延钟,熊宗厚,吴之训,崔先文.频谱激电的异常规律.勘查地球物理勘查地球化学文集(20),罗延钟,史保连,朴化荣主编.北京:地质出版社,1996.
[49] 吴之训,崔先文,罗延钟.频谱激电法在金属矿和油气藏勘查中的应用.勘查地球物理勘查地球化学文集(20),罗延钟,史保连,朴化荣主编.北京:地质出版社,1996
[50] 罗延钟,孟永良,熊宗厚.频谱激电法的模拟准则.勘查地球物理勘查地球化学文集(20),罗延钟,史保连,朴化荣主编.北京:地质出版社,1996
[51] 林长佑,武玉霞,罗东山,等.频率测深中单分量视电阻率—相位联合反演.煤田地质与勘探.1993,21(3):46~52
[52] 黄力军,刘瑞德,陆桂福.相位激电测深在有色金属矿产勘查中的应用.物探化探计算技术.2006,28(1):22~25
[53] 罗延钟.频谱激电法及其在近几年的发展,勘查地球物理勘查地球化学文集(20),罗延钟,史保连,朴化荣主编.北京:地质出版社,1996
[54] 鲍光淑,张碧星,敬荣中,等.三维电磁响应积分方程法数值模拟.中南工业大学学报.1999,30(5):472~474
[55] 张碧星,鲁来玉,鲍光淑,等.电张量格林并矢与三维电磁模拟.中南工业大学学报.1999,30(6):575~577
[56] 夏训银,张宪润,张碧星.均匀半空间三维体的三频激电相对相位法积分方 程模拟.桂林工学院学报.2004,24(4):422~425
[57] 夏训银.相对相位法数值模拟和物理模拟研究:[硕士学位论文].长沙:中南大学,2000
[58] 朴化荣,薛爱民,金东和,等.积分方程求解三度极化体的激发极化效应.物化探计算技术.1985,17(3)
[59] 雷照银,马信山.均匀半空间导体的并矢格林函数解析式.地球物理学报.1997,40(2)
[60] 米萨克,纳比吉安.电磁法(理论部分).赵经祥等译.北京:地质出版社,1992
[61] 唐章洪,薛赛男,冯峰等编著.Visual Fortran程序设计.北京:人民邮电出版社.2000
[62] 郑咸义编.数值计算方法与FORTRAN语言.北京:电子工业出版社.1985
[63] 钱乐秋,陈增荣,朱扬勇编著.FORTRAN语言程序设计与数值计算.上海:复旦大学出版社.1997
[64] 徐士良编.FORTRAN常用算法程序集.北京:清华大学出版社.1992
[65] 廖秀英.基于激发极化数据处理的系统研究:[硕士学位论文].长沙:中南大学,2006
[66] Zonge, K.L. and Wynn, I.C. Recent advances and applications in complex resistivity measurement, Geophsics, 1975, 40(5)
[67] ZONGE INC. GDP MANUAL. 2003
[68] 何继善,汤井田,王少武.直接同时提取激电效应的方波相干法.地球物理学报.1994年,2:516~523
[69] 刘海飞.高密度数据处理方法研究:[硕士学位论文].长沙:中南大学,2004
[70] 罗延钟,张桂青编著.电子计算机在电法勘探中的应用.武汉:武汉地质学院出版社.1987
[71] 李庆扬,王能超,易大义.数值分析,武汉:华中理工大学出版社,1982
[72] 李棋,物探数值方法导论,地质出版社,1991
[73] 李庆扬等,数值分析,华中理工大学出版社,1982
[74] 胡家赣,线性代数方程组的迭代解法,科学出版社,1991
[75] 汤井田,宋守根,何继善.小波变换与阻抗函数的奇性分析及其应用(摘要).中南矿冶学院学报.1994,25(5):29
[76] 宋守根,汤井田,何继善.小波变换与电磁法勘探中高维地质体的特征刻划.中国有色金属学报.1993,3(4):1
[77] 杨福生.小波变换的工程分析与应用.北京:科学出版社,2000:178~179
[78] 宋国乡,甘小冰.数值泛函及小波分析初步.郑州:河南科学出版社,1993
[79] 邓东素,彭立中.小波分析.数学进展,1991,Vol.20(3),294~310
[80] 孙延奎.小波分析及其应用.机械工业出版社,2005
[81] 解海军,陈明生,阎述.利用小波分析压制静态效应.煤田地质与勘探.1998,,26(4)
[82] 秦前清,杨宗凯.实用小波分析.西安:西安电子科技大学出版社.1995
[83] 侯遵泽,杨文采.小波分析应用研究.物探化探计算技术.1995,17(3)
[84] Alexandrescu M Gibert D. Detection of geomagnetic jerks using wavelet analysis. Journal of Geophysical Research. 1995;(B7): 12557~12572
[85] Mallat S Hwang WL. Singularity detection and processing with wavets IEEE trans. On Information Theory, 1992; 38: 617~643
[86] Daubechies I. Orthonormal bases of compactly supported wavelets. Communicationson Pure and AppliedMathematics, 1988; XLI: 909~996
[87] 宋守根,汤井田,何继善.小波分析与电磁测深中静态效应的识别、分离及压制.地球物理学报.1995,38(1)
[88] 林春丽,王克成,黄轶.小波变换与傅立叶变换在信号消噪中的对比研究华北科技学院学报.2005,12(4)
[89] 柳建新,韩世礼,马捷,等.小波分析在地震资料去噪中的应用.物探化探计算技术.2005,3
[90] 牛全营,王长丰,杨宗岭.小波分析在信噪分离中的应用原理.焦作大学学报.2006,03
[91] 韩世礼.时间域多参数联合测深法的应用研究:[硕士学位论文].长沙:中南大学,2006
[92] 林彦杰.基于MATLAB的FIR数字滤波器的设计.电子工程师.2005,01
[93] 董绍平,陈世耕,王洋.数字信号处理基础(修订版).哈尔滨:哈尔滨工业大学出版社,1996.
[94] 刘乃琦编著.IBM PC混合语言编程技术.北京:电子工业出版社.1999
[95] 张超锋,黄艽.Visual Basic 6.0与Visual Fortran的混合编程.研究生论丛.2003,No.2:69~72
[96] 中国科学院数学研究所主编.BASIC语言常用算法程序汇编.北京:中国铁道出版社,1981
[97] Francesco Bslena著.Visual Basic 6编程技术大全,翔实翻译组译,北京:中国机械工业出版社,1981
[98] 李兰友,万振凯,等编著.Visual Basic 6图像处理开发与实例.北京:电子工业出版社.2000
[99] 黄嘉辉,黄悦珊编著.Visual Basic与Windows APt程序设计高手.北京:清华大学出版社.2001
[100] John Paul Mueller著.VisualBasicCOM+编程宝典.刘建军,卢刚,宋江洪等译.北京:电子工业出版社.2001
[101] 汤洪志,周华中,胡文平,周亚中.VB与Visual Fortran混合语言编程及其在物探化探数值计算中的应用.物化探计算技术.2002,24(1):71~76
[102] 张学胜.用VB和Fortran混合编程开发科学计算软件.计算机应用.2003,23:11~12
[103] 朱从旭,邓宏贵.用VB和Fortran混合编程开发科学计算与作图软件.计算机应用.2000,20(4):65~67
[104] 王峰.用VISUAL BASIC实现OpenGL三维地形显示.测绘与信息工程.2000.13(3):13~16
[105] 肖晓玲等.VC与Fortran混合编程及其在大地电磁测深中的应用.物探化探计算技术.2000(2),22(1)
[106] 白世彪,陈晔,王建.等值线绘图软件SURFER7中九种插值法介绍.物探化探计算技术.2002,24(2):157~162
[107] 张二勇,李云峰,王玮.Surfer软件绘图接口的开发及应用.地下水.2005,27(3):212~214
[108] the golden software Inc. the manual of surfer 8.0.2001
[109] 刘川宁.GRAPHER、Surfer32与Core1DRAW10在物探成图中的应用.矿产与地质.2003,17(1),62~64
[110] 武汉中地信息工程有限公司.MAPGIS操作手册,2005
[111] 贾志刚编著.OpenGL编程:入门与提高.北京:中国环境科学出版社.1999
[112] 和平鸽工作室编著.OpenGL高级编程与可视化系统开发(系统开发篇).北京:中国水利水电出版社.2002
[113] The Rockware Inc. the manual of the rockworks 2004. 2004;
[114] 藤蔓等.最新石油地下数据形象化软件Rock Works2004.石油工业计算机应用.2004,12(3),23~25
[115] 童艳光,刘春明等.ROCKWORKS三维成图在物探找矿工作中的应用研究.物化探计算技术,2006(待刊)
[116] 郝永鸿.地质成图方法研究.黑龙江国土资源.2005,5,42~43
[117] 唐令广.CAD成图软件在矿山生产中的应用.矿业工程.2004,2(2),56~58
[118] 高淑芳.地质找矿过程中如何提高物探工作质量.地球物理学进展.2005,20(1),97~103