
金属矿地球物理勘探技术与设备:回顾与进展
吕庆田1,2,张晓培3,汤井田4,金胜5,梁连仲6,牛建军3,王绪本7,林品荣2,姚长利5,高文利2,顾建松8,韩立国9,蔡耀泽10,张金昌11,刘宝林12,赵金花1
1中国地质科学院矿产资源研究所
2 中国地质科学院地球物理地球化学勘查研究所
3 吉林大学建设工程学院
4 中南大学地球科学与信息物理学院
5 中国地质大学(北京)地球物理与信息技术学院
6 北京奥地探测仪器有限公司
7 成都理工大学地球物理学院
8 中国船舶重工集团公司第七一五研究所
9 吉林大学地球探测科学与技术学院
10 中地装(重庆)地质仪器有限公司
11 中国地质科学院勘探技术研究所
12 中国地质大学(北京)工程技术学院
作者介绍:吕庆田,研究员,博士生导师,863计划重大项目首席专家,主要从事矿产勘查技术方法及应用研究。
地球物理探测和钻探工作,是深部矿产资源勘查必不可缺的技术方法手段。2023年,《深部矿产资源勘探技术与装备》项目荣获自然资源科学技术奖特等奖。
长期以来,我国地球物理勘查技术和仪器严重依赖国外进口,同时我国也在自主研发,努力赶超国外先进技术水平。近年来不断涌现的勘查新技术和新装备成果,有力的支持了新一轮找矿突破战略行动,为国家矿产资源安全保障提供了探测技术支撑。
“十二五”国家高技术研究发展计划(863计划)资源环境技术领域,在已有研究成果基础上,设立了“深部矿产资源勘探技术”重大项目,目标是深部500-2000m金属矿精细探测,提高深部矿产资源探测的深度、精度、分辨率和抗干扰能力,共设置了15个研发课题,涉及4个方面:
课题15:金属矿地震处理、解释新技术与软件系统研发。
以上15个课题涵盖了金属矿勘查所需的重要探测技术方法和装备。经多年攻关研究,突破了高精度微重力传感器、铯光泵磁场传感器、宽带感应式电磁传感器等10余项关键技术;研发、完善和升级了地面高精度数字重力仪、质子磁力仪、大功率伪随机广域电磁探测系统、分布式多参数电磁探测系统等18套勘探地球物理仪器设备;创新和完善了20余项勘探地球物理数据处理、正反演方法,研发和完善了2套适合金属矿数据处理及解释的大型软件系统和8套其他专用软件系统,大幅度提升了我国地球物理勘探技术水平。
“十二五”国家实施的“深部矿产资源勘探技术”重大项目,不仅对以往勘查技术研究成果进行了全面总结或升级,创新研发取得了一大批技术与装备突破成果,更是提前吹响了我国新时代地质矿产勘查技术与装备自主研发的号角,为近年来很多勘查技术创新成果提供了研究方向和技术基础,为有效应对国外技术封锁和关键矿产管制提供了有力支撑。
本文回顾了金属矿勘探技术发展历程,介绍了金属矿勘探技术新进展,分析了技术突破挑战和下一步研发方向。文中内容极其丰富,全面展示了地球物理重、磁、电、震和钻探、井中物探等勘查技术新进展。
人类利用矿产资源的历史已经有几千年,然而,发明和利用地球物理勘探技术进行探矿的历史只有100多年。近百年的统计数据告诉我们,近地表发现矿床的概率越来越小,矿产勘探逐渐向深部(500~2000m)和“难进入”地区转移。随着找矿难度越来越大,勘探技术越加重要。有资料表明,近20年全球范围重要矿床的发现都与地球物理技术有关。从历史上看,矿床发现数量与勘探技术的创新密切相关,每次新技术的出现都促进大批矿床的发现。比如,20世纪50年代到70年代中期(1950—1975),由于激发极化(IP)法、航空电磁法,γ能谱和质子磁力仪的应用,促进全球一大批矿床的发现,占同期发现矿床的69%;而1970至1975出现新一轮矿床发现高峰则是由于勘探数据处理、集成分析和反演成像技术的进步。因此,随着勘探深度越来越深,勘探地球物理技术的重要性不言而喻。遗憾的是对于金属矿来说,没有一种技术像反射地震在沉积盆地中寻找油气那样有效,必须使用各种重磁、电法和电磁探测综合技术圈定深部找矿靶区,通过钻探实现最终找矿突破。
近30年,国际勘探地球物理技术突飞猛进,代表性技术包括:航空重力、重力梯度测量系统,直升机吊舱式时间域电磁探测系统,地面大深度3D分布式DCIP探测系统,电磁模拟反演技术,多参数重、磁处理及解释工具以及3D物性反演技术;金属矿地震技术及应用;井中地球物理勘查技术,包括井中地震VSP、井中重力等技术;地质与地球物理数据集成建模技术等等,这些技术的创新极大提高了矿产勘查的深度、精度和效率,带来了新一轮的勘探技术“革命”。本文的目的不是介绍国际勘探地球物理技术的发展情况,这里不多费笔墨,感兴趣的读者可参阅Nabighian和Asten(2002)、Nabighian等(2005a,b)、Vallee等(2011)、Hansen (2001)、Malehmir等(2012)和Patrick(2017, 2018)等综述性文章,了解国际勘探技术的最新进展的详细情况。
我国地球物理勘探技术已有80多年的历史,经过几代人的不懈努力,通过引进、吸收和创新,建立了较为完善的勘探地球物理技术体系,为国家矿产、能源的发现做出了巨大贡献。然而,我国勘查地球物理技术与国外先进水平相比差距较大。航空地球物理勘探技术除了航磁、航放技术外,航空重力、重力梯度,航空电磁探测技术等目前还形成不了实质的勘探能力,一些技术仍属于空白;地面地球物理勘探技术在传感器核心技术、全三维数据采集以及三维数据处理解释等方面还有进一步提升空间;井中(包括测井技术)和井间勘探技术在测量参数、发射功率以及解释技术方面存在一定差距,井中重力测量尚属于空白。
为了加快我国勘探地球物理技术的发展,满足未来深部矿产勘探的需求,“十二五”期间,国家863计划资源环境技术领域设立了“深部矿产资源勘探技术”重大项目。项目在前人研究的基础上,瞄准国际矿产勘查技术前沿,开展核心技术攻关,研制大深度实用化仪器装备,发展深部探测方法技术,创新和完善我国深部矿产资源勘探技术体系,缩小与国外的技术差距,初步实现从“跟跑”到“并跑”的技术跨越。为了便于读者全面了解我国地球物理技术勘探的发展情况和重大项目取得的成果,本文首先回顾我国地球物理勘探技术的基础和发展历史,然后重点介绍863重大项目在核心技术、仪器设备和方法技术方面取得的主要成果,最后对我国下一步地球物理勘探技术的发展方向提出作者个人的认识和建议。
我国的矿产地球物理勘探最早可以追索到1930年,以李四光在《地质论评》上发表的《扭转天平之理论》为起算点,至今已有近90年的历史。新中国成立以来,勘探地球物理进入了良性发展时期,建立机构、组建队伍、开拓应用,先后经历了大发展时期(1949—1961)、调整提高阶段(1962—1978)、全面发展阶段(1979—1990)和改革发展阶段(1991—2000)(夏国治等,2004)。早期的物探仪器主要从原苏联和东欧进口,从原苏联引进了地震仪、重力仪、电法仪和测井仪;从东德进口了大量磁力仪,从匈牙利引进了大地电磁仪、扭秤和测井仪;从瑞典进口重力仪等。在大量引进国外仪器设备的同时,我国从50年代着手开始国产仪器的研制,并开展相关方法技术的研究。直到今天,在自主研发的道路上前进,国产仪器设备、自主研发的软件系统为我国的矿产勘探做出了重要贡献。下面从重、磁、电、震、井中物探及测井、钻探等六个方面进行回顾,对取得的成果和存在的问题进行总结。
重力测量仪根据测量原理可分为相对重力测量与绝对重力测量两种。20世纪60年代,西安石油仪器厂、北京地质仪器厂研制并生产了一些低精度金属弹簧、石英弹簧重力仪。直到1975年,北京地质仪器厂生产的ZSM 3型重力仪投入批量生产,开始了重力仪国产化的新阶段。到20世纪80年代末,北京地质仪器厂先后生产了ZSM 3、ZSM 4和ZSM 5三种型号的石英弹簧重力仪,ZSM 4重力仪器量程≥5000 mGal,读数精度为0.01 mGal,观测精度≤±0.03 mGal,混合零漂≤±0.1 mGal/h(陆其鹄等, 2007;吴天彪,2007)。到80年代末,已经累计生产千余台。虽然,这些重力仪只达到了中等精度,但基本上满足了我国区域重力测量和军方的需要,也使我国成为国际上少数可以批量生产重力仪的国家。此外,中国地震局也研制了DZW型相对重力仪,主要用于对地球重力固体潮的长周期连续性观测。
我国是最早进行绝对重力仪研制的国家之一,从20世纪60年代开始至今,中国计量科学研究院先后研制了NIM-1、NIM-2、NIM-3A型等绝对重力测量仪器。清华大学研制了T-1、T-1A、T-2型绝对重力仪等。此外,中国地震局研制了新型激光干涉绝对重力仪,上述仪器精度可达5~10 μGal左右。
超导重力仪是利用超导原理和技术以测量相对重力变化的新型仪器(郑家祺等,1988)。它利用超导体在超导转变之后的无限导电性和完全抗磁性建立超导磁悬浮系统,以代替常规重力仪的弹簧支撑机构。目前国内在超导重力仪研制方面还处于起步阶段(刘向东等,2019)。重力测量的另一个研究前沿是高精度原子干涉重力梯度仪,国内在此方面起步较晚,主要有中船重工集团第七一七研究所、中国科学院武汉数学与物理研究所、华中科技大学、浙江大学等单位开展了相关技术研究,并取得了较大进展(魏学通,2017)。总体而言,由于仪器系统的技术综合性强、研制难度大,绝对重力仪仍有一些关键核心技术问题有待突破。已有的研究成果主要侧重于实验室基础物理理论分析与实验验证,且由于目前体积较大,可移动性差,在矿产资源勘查中的应用还有较长的路要走。
磁力仪的研发经历了从机械式到质子、磁通门,再到光泵、超导和冷原子磁力仪的发展历程。北京地质仪器厂从1958年开始一直到1991年,连续生产了11种型号的机械式磁力仪,目前机械式磁力仪已逐渐被各类新型电子磁力仪取代。质子磁力仪是我国物探领域应用最多的磁力仪器。北京地质仪器厂从80年代开始不断研制完善了CZM系列质子磁力仪。重庆奔腾数控技术研究所、核工业京核鑫隆公司、廊坊瑞星仪器公司等单位也在引进学习国外仪器的基础上,推出了WCZ系列、G856F型、PM-1A型等多种质子磁力仪,国产质子磁力仪灵敏度都在0.1 nT左右(吴天彪,2007)。磁通门磁力仪在矿产勘查中使用很广泛, 国内研发生产单位也很多。如北京地质仪器厂先后推出了CCM-4、CTM/BS-1、CTM-DT06、CTSD-1、CGM-02D、CSM-05等多种型号的磁通门磁力仪。中国科学院地球物理研究所也研制出了CTM-302型三分量高分辨率磁通门磁力仪(刘士杰等,1990),中国地震局地球物理研究所还研制了磁通门磁力仪野外台阵观测系统(王晓美等,2012)。上述磁通门磁力仪的分辨率也都在0.1 nT左右。基于磁通门技术的磁梯度测量仪器,我国于1975年投产,到1996 年更新换代,灵敏度达到0.1 nT/0.5 m(吴天彪,2007)。
光泵磁力仪、超导磁力仪和原子磁力仪是具有更高精度的磁力探测仪器。20世纪70 年代,北京地质仪器厂曾开始研制台站式铯光泵磁力仪,沈阳仪器仪表研究所也曾着手研制便携式半导体氦光泵磁力仪,由于条件限制,当时的研究产品未能推广使用(吴天彪,2007)。从90 年代中期开始,中国国土资源航空物探遥感中心研制出HC-95型手持式氦光泵磁力仪(灵敏度0.02 nT)和HC-2000 型航空氦光泵磁力仪(灵敏度0.003 nT)。中船重工集团第七一五研究所也研制出了GB系列、RS系列等氦光泵磁力仪。上述磁力仪主要用于海洋磁测和航空磁测,在地面磁测中的应用相对较少。总体而言,我国在光泵磁力仪灵敏度、稳定性和实用性方面接近国外水平,但尚有差距,相关产品在矿产勘查领域发挥了重要作用。关于超导磁力仪,目前我国还处在发展完善之中,已开展的研究主要包括:中国科学院物理研究所研究的高温DC-SQUID平面梯度计(郎佩琳等,2004)、吉林大学研究的高温 RF-SQUID 梯度计样机(赵静等,2011)、中国科学院上海微系统与信息技术研究所研究的低温 SQUID磁强计与梯度计、北京大学研究的超导量子干涉仪磁强计(王赤军,2009)等。中国地质科学院地球物理地球化学勘查研究所从1989 年开始进行高温超导磁强计在地球物理领域的应用研究,目前研制了三分量高温超导磁强计,分辨率达到97 fT/√(Hz),带宽0~30 kHz,可初步应用于电磁法勘探(陈晓东等,2002;陆其鹄等,2009)。除了上述类型的磁力仪,一些单位开展了基于冷原子干涉原理的磁力仪研制(又称全光学磁力仪),国内仍处于起步阶段,但进展很快。
综上,我国可用于实际矿产资源勘查的磁场测量仪器基本可以满足需求,但是精度更高、更稳定、梯度容限更大、量程范围更宽的铯光泵测量技术还有待进一步完善和提升,钾光泵技术还是空白。在磁场测量参数多样化、仪器的设计工艺、稳定性等方面与国外尚有一定差距。总体上看,磁测仪器的发展趋势是:高采样率、高分辨率、高精度、大梯度容限、宽的温度适用范围、大屏幕彩显、大容量存储,以及集GPS于一体的智能化和轻便化仪器。
在重力数据处理与解释方面,精细处理方法不断创新,不仅包括干扰消除、异常分离、延拓、求导等常规转换,也包括小波、分形、人工神经网络、遗传算法等各种非线性方法以及计算机可视化技术。20世纪90年代以来,国内各科研单位开始研究一系列的重磁、重磁电或重磁震数据联合反演和建模方法技术,解决了许多靠单一方法无法解决的问题,使资料解释可靠性、精度有了明显提高,形成了一些有影响力的方法和处理软件系统。
重磁勘探方法由于数据获取方便、测量精度高,被广泛应用于金属矿勘查等领域,近年相应的数据处理与解释技术得到快速发展。在重磁异常弱信号的提取、不同深度异常的分离、低纬度变倾角化磁极,以及位场曲面延拓、复杂条件下高精度三维重磁场多参数联合反演等方面得到快速发展,目前我国已初步形成了多种软件系统,实现了数据处理智能化、人机交互可视化、成果图绘制自动化等功能(王万银等,1991;田黔宁等,2001;管志宁等,2002;李华东等,2015;黄朝等,2015;Wang et al。, 2017c)。在重磁地质应用技术方面,岩性填图技术(严加永等,2014;郭冬等,2014)、基于离散反演和广义反演的三维地质建模技术快速发展(Lü et al。,2013; 祁光等,2012, 2014),为矿产勘查提供更多关于结构、岩性和蚀变等信息,极大提高了利用重磁数据解决地质矿产问题的能力。
电法勘探方法众多,主要可分为基于几何测深原理的直流电法(DCIP)、基于频率测深原理的大地电磁法(MT/AMT/CSAMT),以及基于时间域测深原理的瞬变电磁法(TEM)等。
在直流电法仪器研制方面,北京地质仪器厂、重庆地质仪器厂、重庆奔腾数控技术研究所等单位具有悠久的历史。20世纪70—80年代,北京地质仪器厂研制出了DDJ系列多功能电法仪和DWJ型微机激电仪;重庆地质仪器厂研制出了DZD系列多功能直流电法仪、DDC型电阻率仪及DJF型大功率激电仪。随后,重庆奔腾数控技术研究所研制出了WDJD系列、WDDS系列、WDA型等直流激电仪。上述仪器系统在我国早期的金属矿产勘查中发挥了重要作用。此外,北京中西远大科技有限公司、北京大地华龙科技公司、桔灯公司和中国地质大学(北京)等也都推出了各具特色的电法测量系统。中南大学何继善院士等在总结国际上传统的变频激电和奇次谐波激电法的基础上提出了双频激电和伪随机多频激电法,并研发了相应的仪器系统(何继善,1978;柳建新等,2001;陈儒军等,2003)。随着电子技术的进步,以提高采集效率和分辨率为核心的高密度电法仪器迅速发展,比如,北京地质仪器厂研制的DDJ系列、重庆地质仪器厂研制的DUK系列和重庆奔腾数控技术研究所研制的WGMD系列高密度电法仪等,还有吉林大学骄鹏公司、西安澳立华公司、北京大地华龙科技有限公司等单位都开发出了各具特色的高密度电法仪(李志武等,2004;李晓斌等,2008;何刚等,2014)。这些国产品牌的直流探测仪器在工程地质、矿产地质、环境地质等领域发挥了重要作用。从目前国际上直流电法的发展态势看,全波形、大规模三维分布式采集技术和相应的反演解释技术将是未来重要的发展方向。
基于几何测深原理的直流电法受发射极距、发射功率的影响,探测深度有限。借助于天然场或人工场源的感应类电磁法(MT/TEM/CSAMT)具有更大的穿透深度。我国在大地电磁仪研制方面具有较早的历史。原中国科学院兰州地球物理研究所、国家地震局地质研究所及石油部地球物理勘探局仪器厂等,在20世纪60年代中期便开始了大地电磁观测仪器的研制,但因高频特性不好、温漂大、移动不便、仪器动态范围小、数采技术等原因,并未得到进一步推广应用。80年代后,中国科学院地球物理研究所、长春地质学院、中国地质大学(北京)等相继研发出了低频数字大地电磁测深仪、GEM-1宽频数字大地电磁测深仪和长周期大地电磁测深仪等,一些主要技术指标已经达到当时国际同类仪器水平(张秀成,1989;王家映,1997;巩秀钢等,2014)。在瞬变电磁仪研制方面我国也具有较早的历史。原地矿部物化探研究所、原中南工业大学、北京矿产地质研究院等较早开展了TEM仪器的研发。目前,技术相对成熟、且有一定市场的瞬变电磁仪器系统有:重庆地质仪器厂GATEM系列、重庆奔腾数控研究所WTEM系列、北京矿产地质研究院TEMS系列、骄鹏科技公司MDTEM64系列、中国地质科学院地球物理地球化学勘查研究所IGGETEM系统、中国科学院地质与地球物理研究所CASTEM系统、吉林大学ATEM 系列、中国地质大学(武汉)CUGTEM系列和中南大学HPTEM系列等(王华军和梁庆九,2005;嵇艳鞠等,2005;周平和施俊法,2007;付志红等,2008;武欣等,2016;席振铢等,2016)。
在可控源电磁法和多功能电磁法方面,早期几乎全部被国外仪器垄断。21世纪初开始,在国家“十一五”863计划、科技部重大仪器专项,国家公益性行业专项和中国科学院战略先导性专项等项目支持下,我国在仪器设备自主研发方面取得跨越式发展,研制出了适用于三维分布式探测的可控源电磁法及多功能电磁法仪器,并在深部矿产勘查中得到较好的应用。比如中南大学研制的JSGY广域电磁探测系统、中国地质科学院地球物理地球化学勘查研究所研制的DEM多功能电磁探测系统、中国科学院地质与地球物理研究所研制的SEP地面电磁探测系统、吉林大学研制的JLEM大深度分布式电磁探测系统、中船重工第七二二研究所研制的CEMT-03大地电磁探测系统、骄鹏科技公司研制的E60EM-3D多功能电磁法系统、中石油集团东方地球物理勘探公司研制的TFEM时频电磁仪和中国地震局等单位研制的大功率极低频/超低频(SLF/ELF)电磁接收机等(何继善,2010;蒋奇云,2010;林品荣等,2010;王兰炜等,2010;张文秀等,2012;底青云等,2013;吕雪峰等,2015)。在“十二五”863计划支持下,我国在已有的研发基础上,完成了电磁探测技术与装备的升级和性能的全面提升,实现了更大发射功率、更大探测深度、更强抗干扰能力、更多采集参数的功能。这些进展将在本文第二部分进行详细介绍。
感应类电磁探测技术的核心是电磁场传感器。我国市场上使用的高、中、低频段的高灵敏度电磁传感器(包括线圈感应式和超导式)几乎全部被国外垄断,部分低频产品对我国禁止出口。由于磁性材料、工艺等技术的原因,很长一段时间我国还不能生产高灵敏度磁场传感器,尤其是低频磁场传感器,更缺乏对高灵敏度传感器的检测、标定技术和装置,长期以来一直制约着我国电磁探测技术的发展。通过“十二五”的努力,多家单位在高灵敏度、宽频电磁传感器研制方面取得突破性进展,总体达到、局部超过国外主流仪器的技术指标。比如,中国科学院电子所研发的MT 磁场传感器工作频率范围达到10-3Hz~1 kHz,0.1 Hz时噪声水平达到5 pT/√Hz(王言章等,2007; 巨汉基等,2010;陈兴朋等,2012;朱万华等,2015;赵毅等,2016; Du et al。, 2018)。在863计划支持下,重大项目实现了高灵敏度电磁传感器关键部件研制的重要进展,并在检测、标定技术和装置方面取得了零的突破,具体内容见本文第二部分。
近年,虽然我国在电磁勘探核心技术、仪器研制方面取得显著进展,但受仪器的工业化、商品化设计、元器件质量等方面的影响,仪器在稳定性、实用性和可操作性方面与国外设备相比仍有较大差距。占据我国矿产勘查领域的主流仪器仍是V8、GDP32、EH4、GMS等国外仪器。因此,早日实现国产电磁法仪器的实用化和产品化,完全替代国外产品,还需要研发机构和企业紧密联合、共同努力。
随着电磁探测技术在矿产资源勘查、工程勘察中越来越广泛的应用,相应的数据处理及解释技术得到极大发展。在正反演方面,有限差分、自适应非结构化网格有限元、边界元、积分方程、有限体积等高精度数值模拟方法已广泛用于正反演计算;一维/二维反演技术已经非常成熟,三维带地形直流电法、激电法、可控源电磁法、大地电磁法等反演算法也已经实现,并得到较好的应用;电性各向异性的研究也取得系列成果,多参数互约束联合反演技术研究已经展开,并取得进展。中国地质大学(北京)、中南大学、吉林大学、中国地质科学院地球物理地球化学勘查研究所等单位研发了适用于起伏地形的TDIP/SIP、CSAMT、AMT/MT反演解释软件和TEM定量解释软件,极大提高了电与电磁探测技术的实用性。在电磁抗干扰数据处理技术方面,除了常规的时序叠加、频谱分析、数字滤波等技术之外,相关辨识、稳健统计、数学形态滤波、经验模态分解、小波分析、分形技术等手段也被应用于海量电磁数据的自动化处理,为进一步实现电磁三维高精度探测打下了坚实基础。
地震仪器按照数据传输方式可分为三种:有缆遥测地震仪、节点地震仪和无缆遥测地震仪。无缆遥测地震仪由于受到数据传输效率、速度和稳定性等影响,市场应用不太普及,但前景看好。按照震源性质分为人工源和天然源地震仪,前者更多用于勘探,后者多用于探测和监测,近年出现的节点地震仪架起了人工源和天然源的桥梁,即可同时接收人工和天然源信号。从硬件角度,并没有针对金属矿勘探研制的地震采集设备,从实际应用看,金属矿地震勘探多数仍使用油气勘探的有缆数字地震采集设备,比如,法国Sercel公司的428、408、388等;美国I/O公司的System II、System 2000、Image; 加拿大Geo-X公司的ARAM-24、ARAM-ARIES等。我国一直在致力开发具有自主知识产权的有缆遥测地震仪(万道地震仪)以及配套的人工震源和检波器,经过多年的努力,近年中国石油集团公司东方地球物理公司、中国科学院地质与地球物理研究所等单位在数字地震仪器、MEMS检波器、可控震源等方面取得突出进展,国产地震仪器取代进口指日可待。。
国内有多家单位已经研发出多型号的节点地震仪,比如,重庆地质仪器厂研制的EPS便携式数字地震仪、CZS一体化宽频数字地震仪;吉林大学地球探测科学与技术学院在863计划、国家公益性行业专项支持下,研发了无缆遥测地震仪;中国科学院地质与地球物理研究所等单位研发出了无缆存储式数字地震仪等(孙明等,2001;孙明和林君,2002;杨泓渊等,2009;王肃静等,2015;赵金龙等,2016;田入运等,2017),这些仪器已经开展了前期实验,取得一定效果,但距大规模产品化使用尚有一定距离。国内研制的无缆节点式地震采集仪器在电源管理、制造工艺和数据回收系统等方面与国外仪器尚有一定差距(郭建和刘光鼎,2009)。
目前,国内外还缺少针对金属矿地震勘探的采集系统,大多数仪器还是用在油气勘探、煤田和工程勘探上。金属矿产勘探具有地质、地形条件复杂、难“进入”和有效信号弱等特点,客观上要求无缆、轻便、多道、大动态范围的仪器设备,适应金属矿勘探特殊的地表和地下地质条件。
我国金属矿地震数据处理与解释技术整体仍处在发展阶段。由于金属矿区成层性(连续性)差、波阻差异较小、反射信号弱等特点,客观上要求数据处理技术的创新。但是目前还没有看到有针对金属矿产特殊的处理解释技术,总体上还是沿用油气地震勘探的技术,只是在个别处理环节针对金属矿区和硬岩地区的特点进行重点处理,如去噪、均衡、偏移等。近十多年来,也有一些专家针对金属矿产勘查的特点,探索了金属矿散射波的模拟和成像技术,但是更多的还是停留在理论层面。比如,国内的孙明和林君(2001)依据微扰理论,进行了金属矿地震散射波场的数值模拟研究。实验结果表明:可通过地震波散射响应的强弱推断矿体,散射波相干性的好坏与杂乱散射体的不均匀性有关。刘学伟课题组对散射波的基础理论、物理模拟和成像剖面进行了研究,通过大量地质模型(单点、多点、层状、单体、多体)的地震散射波场的数值模拟和正演研究,对单炮记录和叠加偏移剖面上散射波的波组特征与其所反映的地质模型之间的对应关系提出了新的认识(尹军杰等,2005;李灿苹等,2005;勾丽敏等,2007;李敏锋等,2007)。
金属矿区的地震剖面解释一般不能沿用盆地反射地震的解释思路。一般要根据地震反射特征,比如:反射纹理、强弱变化、空间关系、密集程度、反射“亮点”等,结合矿区地质特征、成矿和控矿模式等进行解释。在这方面一些专家做了有意义的尝试。Li和Eaton(2005)在东天山土屋铜矿成功地运用地震方法发现大约 1 km深的与斑岩铜矿有关的岩体和岩体的内部结构;吕庆田等(2004)在铜陵狮子山矿田利用反射地震发现了主要赋矿层“五通组”的深部分布形态;在庐枞的罗河铁矿深部发现了新的反射界面,后经钻探验证,在 1600 m深处发现了小包庄铁矿(吕庆田等,2010)。
虽然金属矿地震借鉴已经成熟的石油地震技术,但仍然面临很多技术上的挑战,如复杂构造、低信噪比、不连续反射和强散射等,要解决这些问题,必须发展新的采集和处理解释技术。
我国是开展地下电磁波勘探较早的国家之一,自20世纪60年代初地矿部物化探研究所研制出第一台仪器并开展工作至今,全国先后有几十个科研单位、高等院校和生产单位开展了地下电磁波层析技术的应用和研究。研究内容涉及理论分析、模型实验、仪器制造、层析算法、软件 研制、解释和应用诸方面。先后研制和生产了具有代表性的 JWT-1、JWT-2、JWQ-3、WKT1-4、JW-4、JW-5 等多种型号的仪器,开发了多种行之有效的层析算法,在矿产勘查、工程地质、环境地质和考古等众多领域中发挥了重要作用。
最近二三十年,中国地质科学院地球物理地球化学勘查研究所、重庆地质仪器厂、重庆奔腾数控技术研究所、北京地质仪器厂、上海地学仪器研究所等单位先后研制出地下电磁波CT、声波CT、井中磁力仪、井中激电测量仪、井中TEM测量仪等多种类型的仪器产品。但整体而言,我国的井中物探装备研发一直处于相对跟踪状态,随着深部找矿工作的深入开展,已有仪器暴露出的问题越来越多(周平等,2009)。目前,我国井中物探技术还面临着测量参数少、发射功率和接收灵敏度低等问题。
在综合地球物理测井仪器方面,代表性的产品有重庆地质仪器厂的JGS型、上海地学仪器研究所的JHQ-2D型、渭南煤矿专用设备厂的TYSC-QB型、北京地质研究院的HD-40型等数字测井系统。但上述仪器在装备性能和批量配套生产方面与西方公司仍有一定差距。中国石油测井有限公司等单位曾开展了高性能测井技术与装备的研发工作,但主要针对油气测井,在固体矿产资源勘查中很少应用。
目前,金属矿小口径测井设备存在测井深度浅、测量参数少、缺乏必要的检测标定装置(参数井)以及现场实时处理解释软件等问题,与国外同类产品存在一定差距(周平等,2009)。因此,研发深孔小口径、多参数、轻便自动化、稳定可靠的综合测井仪器是目前深部资源勘探的一项重要内容。
近20年来,我国在地下物探及测井数据处理解释方面取得了重要进展。国内学者先后开展了井-地电磁三维正反演研究,开发出了井-地电阻率、激电相位三维正反演软件和可视化CT处理与解释系统,可实现地下电磁波、声波数据与钻孔资料的实时交互处理(Wang et al。, 2006; 武军杰等,2017)。中国地质科学院地球物理地球化学勘查研究所等单位还研发了一套包含地下电磁波和声波透视层析成像的地下物探综合工作站;吉林大学、中国地质大学(北京)、中国地质大学(武汉)、成都理工大学、核工业203研究所等单位先后设计实现了测井资料处理解释软件系统(邵才瑞等,2005;原福堂等,2005;邹长春等,2010)。在地质解释方面,中国地质调查局及所属单位先后完成了《地-井脉冲瞬变电磁法物理模拟曲线图册》、《井-地大功率充电法数值和物理模拟曲线图册》,制定了井中磁测、井中激发极化法、钻孔电磁波法、金属矿地球物理测井等相关技术规程,为地下物探异常的推断解释提供了有实用价值的参考资料。上述成果为今后逐步实现地下多物理场多参数多成像系统的综合集成打下了良好的基础。
20世纪 50 年代前,我国的岩心钻探设备主要靠进口。50 年代中期国内开始仿制,60 年代中国地质科学院勘探技术研究所研制出我国第一台岩心钻机;70 年代起,为了适应金刚石取心钻进的需要,开始研制高转速立轴式钻机及与之配套的泥浆泵,成功研制了XY 系列钻机与BW系列泥浆泵(甄玉娜和王均,2001)。在全液压顶驱式钻机研制方面,“八五”期间成功研制了我国第一台车装钻机,主要供300 m以内金属矿勘探空气反循环岩心取样钻进使用。“十五”中期,又通过引进、消化吸收国外先进技术,研发出我国第一台拖装式千米全液压顶驱钻机样机(文田,2007;石艺,2008; 张金昌,2016)。
“十二五”之前国内研制的全液压动力头地质岩心钻机主要有: 连云港黄海机械厂的HYDX-4、HYDX-5、HYDX-6型,钻深能力分别为700 m、1300 m 和1600 m。山东省地质探矿机械厂的XD-3、XD-5型,钻深能力分别为700 m 和1200 m。还有其他一些厂家研制的不同型号的全液压动力头钻机,钻探深度多在2000 m以内。2007年立项的“十一五”国家“863”计划重点科技项目“2000 m 地质岩心钻探关键技术与装备”于2010年底结题,完成了优质、高效深孔钻进工艺研究以及2000 m全液压岩心钻机及配套设备的研制。进入“十二五”以来中国地质科学院探矿工艺研究所联合中国地质装备集团公司研制成功3000 m电动直驱顶驱钻机,中国地质调查局勘探技术研究所研制了3500 m地质岩心钻机、400 m轻便钻机、浅层取样钻机和600 m反循环钻机等(张金昌等,2012;张金昌,2016)。
我国拥有居国际领先地位的液动锤冲击回转钻进技术,该技术将轴向高频冲击载荷施加到钻头上,使钻进效率较纯回转钻进方法大大提高,在东海中国大陆科学钻探工程中,创造了5158 m的液动锤工作深度之最(王人杰和苏长寿,1999)。我国绳索取心钻进技术起步于20世纪80年代初,但因钻杆强度问题存在较大的阻力(孙建华等,2011)。近年,针对钻杆强度薄弱的问题,一些研究单位与施工单位选择了优质管材和加大钻杆壁厚的方法,但钻进孔径59 mm钻孔的小口径绳索取心钻杆强度仍不能很好地解决,因此,难以进行深孔钻进。小口径定向钻探技术达到了相当的水准,80年代末在铜陵冬瓜山铜矿床于一个主干孔内成功施工了6个全方位分布的定向分枝孔,平均深度达856 m,钻孔中靶精度在1.29~5.86 m。在组合钻具的研究上也做了大量的工作。80年代初、中期,我国扩孔翼式不提钻换钻头钻具(钻孔直径75 mm)、绳索取心液动锤钻具(最小钻孔直径59 mm)相继问世,尽管存在较多的问题,但也指明了一定的完善途径与发展前景。2001年,中国大陆科学钻探工程启动,使我国首创的大口径岩心钻探组合工具的研究与应用达到了前所未有的高度。螺杆钻+液动锤提钻取心钻进方法取得了举世瞩目的辉煌成果;螺杆钻+绳索取心钻具、液动锤+绳索取心钻具(绳冲二合一)、螺杆钻+液动锤+绳索取心钻具(三合一)样机入井实钻试验,均有一定的成果,为研究小口径组合钻具积累了可借鉴的经验。
进入“十二五”之后,我国在钻探技术和装备研制方面不断取得突破,本文将在第二部分进行详细介绍。
为加快我国矿产勘查技术的发展,“十二五”期间,科技部设立了“深部矿产资源勘探技术”重大项目,该项目根据金属矿勘查需要多方法综合探测的特点,从重磁探测技术、电法及电磁探测技术、金属矿地震探测技术、钻探及井中探测技术四个方面部署研发任务(图1)。其中,重磁探测技术包括:全数字地面重力勘探仪器、高精度地面量子磁测技术及仪器,以及地球物理多参数约束反演技术与软件系统研发等内容;电法及电磁探测技术包括:高灵敏度、宽频带电磁传感器技术,大功率伪随机广域电磁探测技术与装备,大功率井-地电磁成像系统,分布式多参数高密度电磁探测系统,长周期分布式大地电磁观测系统等研发内容;金属矿地震探测技术包括:轻便分布式遥测地震勘探系统和金属矿地震处理、解释新技术与软件系统研发;钻探及井中探测(测井)技术包括:金属矿小口径地下物探技术与设备、大深度小口径多参数测井技术及仪器研制,4000 m地质岩心钻探成套技术装备和自动化智能化岩心钻探技术与设备研制等内容。通过上述技术方法的研发,为500~2000 m深部金属矿精细探测提供技术支撑。
图1 863重大项目“深部矿产资源勘探技术”研究框架及课题设置示意图
下文按照四个方面的任务部署,分别介绍四个方面取得的成果。首先是进展概述,介绍每个方面在关键技术突破、仪器研制、方法创新和软件研发取得的总体进展情况;然后是代表性成果,选择2~3个“亮点”成果,介绍关键技术突破的思路,或仪器的研制方案和软件系统框架与功能。
关键技术研发方面,项目突破了高稳定性熔融石英弹簧装置、高分辨率差动电容测微(位移)技术、高精度倾斜测量技术和高稳定性恒温与测温技术等核心技术,大幅度降低了零漂,测量灵敏度达到国外同类技术(CG-5)的水平(耿启立,2016;秦佩等,2019);攻克了铯光泵磁场传感器核心技术,突破了铯灯室制作、铯探头设计、宽带信号处理检测系统等核心技术(张谨等,2015;黄成功等,2019)。
在仪器研制方面,成功研制了地面高精度数字重力仪、地面高精度绝对重力仪,其中地面高精度数字重力仪实现了商品化,实现了批量生产能力(耿启立,2016;黄雯迪等,2017);研制出了高精度直流激发型质子磁力仪、动态激发核磁共振磁力仪、DGB-8 型数字式氦光泵磁力仪、铯光泵磁力仪(样机)等4种高灵敏度磁力仪(张谨等,2015;贺宁波,2016;黄成功等,2019)。高精度直流激发型质子磁力仪CZM-863T 是对前几代质子磁力仪性能指标的全面提升。通过磁传感器的优化设计,降低了噪声干扰,增大了梯度容限。对旋进信号进行数字化处理和精确测频,实现磁场值的高精度测量,分辨率达到0.05 nT。动态激发核磁共振磁力仪突破了多项技术难点,填补了国内在动态激发核磁共振磁力仪领域的技术空白,测程范围为20000~120000 nT,静态噪声0.01 nT。DGB-8 型数字式氦光泵磁力仪是基于氦光泵技术的新一代数字化磁力仪,用于总场、梯度的连续测量,测程范围为10000~100000 nT,静态噪声0.005 nT。CZM-863T和DGB-8磁力仪设备集成了GPS 定位模块,系统软件功能齐全,操作简便,适用于地面、低空磁法勘探。铯光泵磁力仪测量技术实现零的突破。
在方法创新和软件系统研发方面,开发出了先进的重磁数据处理、反演解释技术和软件系统。实现了复杂地表条件下的重、磁三维约束反演、重-震匹配三维反演、直流电阻率三维反演、激发极化法三维反演、地磁与井中磁测联合多参量三维反演、重力及其梯度数据三维联合反演等反演解释技术(李泽林等,2019;刘强等,2019;于平等,2019;张刚和吕庆田,2019);优化了包括化极求导在内的十八种常规方法、曲化平、低纬度化极、场源边界识别等数据处理技术,通过软件集成开发,研发出了一个用于重磁电数据处理与正反演解释的综合软件(刘彩云等,2015;英高海等,2016;王泰涵等,2019;汤井田等,2019;郑元满等,2019)。并在下列关键技术和算法上实现创新:
使用小波压缩、等效正演计算、数据自适应采样、CPU/GPU并行计算等手段,实现了大数据重磁三维反演计算,使得大数据三维反演达到实用化;
提出重磁三维稀疏反演算法,有效地利用已知的物性信息,实现了深度分辨率高、具有尖锐边界的反演效果,有效提高了反演的分辨率;
提出了重-震匹配反演算法,通过引入重-震匹配项,提高了反演结果的横向和纵向分辨率。提出了重磁场及其梯度数据三维联合反演算法,有效地利用已知的多参量重磁场数据,改善反演结果的深度分辨率。
高精度重力传感器是高精度重力仪的核心部件(图2a),主要由石英弹簧系统、位移传感电路、恒温测温电路、倾斜测量设备、真空仓等组成。它通过作用在重荷(动极板)上的重力、石英弹簧的弹力、反馈静电力间的平衡原理,实现地面重力变化的精确测定(图2b)。高精度重力仪在我国具有较长的研发历史,原北京地质仪器厂曾经成功研制出ZSM-3、4、5型中精度重力仪,但由于材料、微位移测量、温控/温补技术不过关,一直达不到高精度测量标准。
图2 熔融石英重力传感器(a)、石英重力传感器原理示意图(b)和ZSM-6高精度数字重力仪(c)
石英扭丝导电膜制作工艺。对于直径仅有几十微米的石英扭丝,采用了特殊的技术和工艺实现其导电膜的制作,使其具有良好的导电性和附着力;
石英部件退火工艺。用于高精度重力传感器的石英部件,采用特殊的工艺进行退火处理,消除了石英部件内应力;
石英弹簧的灵敏度调整工艺。经过特殊的抗压精密调整,精细调整石英弹簧的粗细,保证其弹力系数符合设计要求;
石英传感器参数选择。精细结构参数和电参数的设计,确定测微传感器基础电容值和测程范围;
三维限位调整。通过限位间隙百纳米级三维限位精密调整,减小弹性后效对重力仪读数的影响,提高仪器使用的安全性。
高精度恒温测温技术。通过使用双层恒温电路,并通过PID闭环控温电路控制其温度,实用精密温度差动全桥测量电路检测残余温度变化,精确测出温度变化后设计软件算法进行自动温度补偿,最后确保测温精度0。00001 ℃,重力温度补偿精度为1微伽/0。00001 ℃。
新研制的ZSM-6地面高精度数字重力仪,各项技术指标基本达到国际先进水平(图2c,表1),实现了温度、倾斜、漂移和固体潮自动补偿改正,重力值测量过程全自动化,填补了我国高精度数字重力仪的空白。目前已实现批量生产,可以满足地质科学研究、矿产资源勘探、工程勘查、地质灾害调查和国防建设等领域的需求。
表1 ZSM-6数字重力仪与CG-5重力仪指标对比
铯光泵磁传感器为单光系结构,包含:铯灯室、光学镜片、光敏传感器(光电转换)、圆偏振片、铯灯室加热和恒温部件等。铯灯、铯吸收室为玻璃封装结构,内部充入高纯度的铯原子(Cs133)金属,铯吸收室被加热线圈加热汽化并控制45℃恒温。铯灯工作时被功率模块激励发光输出894.3 nm的泵浦光,铯原子气室内铯原子在光泵作用下产生磁共振。通过研究铯灯、铯室制作工艺以及铯灯室制作系统,设计铯光泵传感器探头,研制信号处理电路、恒温控制电路等,获取光磁共振信号,并输出磁场信号,完成了铯光泵磁力仪原理样机研制,实现了我国铯光泵磁力仪研制零的突破,样机技术指标基本达到设计指标(表2)。
表2 铯光泵磁力仪(样机)的技术指标
铯灯和铯吸收室的制作工艺:通对铯原子气室缓冲气体最佳压力计算和试验分析、仿真计算,铯灯室样品的磁共振线宽测试和检测分析等工作,完成关键部件铯灯、铯室的设计和样品制作;
铯光泵磁力仪探头的方向误差:从铯光泵的物理基础、铯的磁共振信号、“分裂波束”设计入手,分析铯光泵磁力仪的Sx信号及角相关性,并在探头光学设计中采用分裂波束设计,解决了铯光泵磁力仪探头的方向误差问题;
铯光泵磁力仪信号检测技术:通过激励、恒温控制、磁共振信号调理、90°宽带移相、拉莫尔计数器电路的设计,实现对铯光泵传感器输出的磁共振信号的检测。
地球物理数据转换与反演是获得地质认识的基础,它通常通过软件系统来实现,在勘探地球物理学科中有极其重要的地位。项目以现代反演解释技术为重点,在前期研发的基础上,实现了多种有效的重磁三维反演方法技术,解决了起伏地形、大数据、有剩磁影响等复杂情况下的三维反演关键技术问题。形成了一套交互性强、方法先进,集数据管理、处理、地质建模和反演的多功能一体化的软件系统。
软件系统架构采用自下而上分层设计策略,包括:底层基础系统层、软件框架系统层、基础功能模块层、专业处理模块层、软件交互操作界面层和系统集成层。软件系统采用模块化设计,模块之间通过相应接口进行调用,使得各个模块具有良好的封装性和独立性。根据功能划分子系统,一共划分为5个子系统,即:重磁常规处理子系统、辅助处理工具子系统、图形图像绘制子系统、交互处理与解释子系统和多参数反演子系统。各系统功能:
① 重磁常规处理子系统包括:化极、延拓、滤波、边界增强等处理。开发实现了21个处理模块:磁异常化极、磁异常三分量转换、磁异常任意方向分量转换、磁源重力异常、重力插值切割、向上延拓、向下延拓、一阶总导数模、一阶水平导数模、一阶任意方向导数模、二阶任意方向导数模、补偿圆滑滤波、正则化滤波、深度信息分离提取、构造增强、趋势分析、窗口滑动平均法光滑、圆周平均法光滑、等效源曲化平、等效源化极、全张量梯度数据组合分析识别场源边界(图3)。
图3 地球物理多参数反演解释软件系统:重磁数据处理和可视化界面
② 辅助处理工具子系统包括:网格数据基本信息查看工具、网格数据任意剖面数据提取工具、三维数据任意切片提取工具等。
③ 图形图像绘制子系统包括:二三维数据成图可视化功能,我们开发实现了剖面曲线图、平剖图、等值线图、颜色栅格图、立体晕渲图、三维曲面图、三维等值面显示图、三维切片显示图、三维阈值显示图。
④ 交互处理与解释子系统包括:二度半模块交互建模处理解释系统、三维模型交互建模处理解释系统。
⑤ 多参数反演子系统包括:三维重磁光滑反演、三维重磁聚焦反演、三维重磁稀疏反演、三维重磁场及梯度数据联合反演、重震匹配反演、井地磁多参量联合反演、二维电阻率反演、三维电阻率反演、二维激发极化法反演、三维激发极化法反演等。
关键技术研发方面,项目通过磁芯的导磁材料的组份配比、制备方法、热处理工艺参数研究,以及感应式电磁传感器的制造及装配工艺研究,攻克了高灵敏度宽频感应式电磁传感器技术,研制出了适用于瞬变电磁法的TEMC104型电磁传感器,适用于CSAMT/MT的MC-30、MC-50型电磁传感器,在一致性、稳定性和重复性上超过国外同类产品。针对我国缺乏感应式电磁传感器检测和标定设备及规程的现状,项目研发出了适用于高灵敏度感应式电磁传感器检测系统,并制定了相应的检测和标定技术规程,为进一步提升电磁传感器的技术性能提供了研究平台(汪杰等,2017;王勇等,2019)。
采用基于传统双磁芯差分探头结构和检测线圈与反馈线圈共用的技术方案,通过磁通门建模和电磁场仿真等手段论证了磁通门结构、尺寸和性能指标,提出制作方案,成功研发出高灵敏度三轴分离式磁通门传感器(LMT-FS02),经过第三方机构检测主要指标达到国际同类产品(英国Bartington公司的Mag-03传感器)的水平,频响范围DC-1 kHz时,单轴噪声水平约为6pT/√(Hz)@1Hz,三轴正交误差<0.5%(高嵩等,2019)。
为实现硬岩表面电场的有效观测,借鉴于医疗设备中生物电信号观测的耦合式电极原理,成功研制出电容耦合式电场传感器,有效解决了井地电磁方法在坑道内进行电磁场观测的特殊需求(Chen et al。, 2017; 陈凯等,2019)。
在仪器研制方面,研制、完善和发展了分布式多参数高密度电磁探测系统、大功率伪随机广域电磁探测系统、坑-井-地电磁成像系统和长周期分布式大地电磁观测系统,具体技术指标达到或超过国外同类产品,前三者实现小批量生产和初步的产品化。分布式多功能电磁法仪器、伪随机广域电磁法仪器已在10多个地区进行了应用,取得了良好的地质效果。相比常规电磁仪器,具有抗干扰能力强、勘探深度大、分辨率高、工作效率高等优点。三套设备的研发成功,将逐渐改变我国多功能电法仪器依赖国外的现状,提高国产电法/电磁法勘探的效率,满足不同勘探目标对多功能电磁勘探方法(DCIP/SIP/CSAMT/AMT/MT/LMT)的需求。
地球物理电磁法理论、方法取得创新性进展,形成了以全域电阻率、多场源激励和信噪分离为核心的时空阵列(广域)电磁法技术体系。包括:提出了多场源激励-阵列分布式接收的三维电磁法勘探方案;实现了电性源和磁性源所有电磁场分量及其比值的广域视电阻率定义及稳定计算,为矢量和张量广域电磁法测量与应用奠定了基础;提出了基于磁场旋度的视电阻率定义;研究了多种强干扰情况下的信噪分离算法和评价方法,包括基于压缩感知的稀疏表示与正交匹配追踪的强干扰分离算法,基于FFT-EEMD及移不变稀疏编码(SISC)、基于周期信号特征及逆FFT的电磁法信噪分离方法,以及各种噪声评价和有效信号筛选方法。建立了多场源输入-多时刻激励-多站道输出的电磁勘探模型,实现了基于场源的信号分离,不仅可压制输入端强相关干扰及其他噪声,还可以统一大地电磁法、可控源音频大地电磁法和广域电磁法的野外采集和数据处理工作,初步形成了时空阵列广域电磁法理论技术体系(周聪等,2019)。
在人工源电磁法正反演方面,提出并实现了基于有限元-无限元耦合的人工源电磁法三维正演与反演、带任意地形的基于非结构化网格和虚拟场的体-面积分方程法正演,以及任意复杂模型的波数空间混合域人工源电磁法三维快速正反演(汤井田等,2014;Ren et al。, 2018a; Zhou et al。, 2018; 原源等,2019)。在此基础上,开发了基于工区管理的广域电磁法数据处理、三维正反演和可视化交互解释软件平台,极大地提升了广域电磁法资料处理与解释水平。
在常规电磁法正反演方面,提出了多种大深度电磁探测三维正反演方案,实现了基于工区的可视化三维反演与交互解释软件方案(郭泽秋和董浩,2019;李勇等,2019;罗威等,2019; 王宁等,2019)。研究了基于FFT的快速2.5D正反演;提出了坑-地一体化带地形三维反演算法;还研发了其他电磁软件系统,界面友好、交互性较好的分布式长周期大地电磁数据处理软件(LMTPro-V1.0)和长周期大地电磁二、三维正反演解释软件LMT3D-V1.0等。
电磁传感器是电磁探测仪器的核心部件,我国电磁勘探仪器设备严重依赖国外进口,特别是高性能电磁传感器更是为国外少数几家公司垄断,一些低频传感器对我国采取严格的出口管制。根据国内电磁探测设备的需要,项目开展了瞬变电磁(TEM)、音频电磁(CSAMT/AMT)和大地电磁法(MT)电磁传感器的研制。
感应式电磁传感器一般包括:磁芯、感应线圈、标定线圈、前置放大电路、电反馈电路、屏蔽层、支撑及封装件等部分组成。项目从磁芯材料、装配工艺、检测技术等方面开展研制工作,通过大量的对比试验,优选出适合制作高性能感应式电磁传感器磁芯的导磁材料的组份配比、工艺制备方法、热处理工艺参数,研制了一套高性能电磁传感器的制造及装配工艺(图4),使得研制的高性能电磁传感器在一致性、稳定性和重复性上达到,或超过国外同类产品的技术水平。
图4 全数字化一体式设计,具温度补偿、标定与数据采集功能的MT传感器
磁芯软磁材料制造工艺创新:加入由某种合金材料制成的磁芯,再开展组份配比、热处理工艺参数、制备工艺方法以及残余应力消除方法等方面的研究,通过磁性参数的对比分析,形成一整套制造高灵敏度电磁传感器磁芯的软磁材料的制造工艺参数;
装配工艺创新:通过封装材料选型、屏蔽材料选型以及部件之间连接方式的研究,避免磁芯受力,以降低传感器使用的时效性,研发一套从材料、加工以及装配的完整的制造工艺流程;
磁屏蔽技术及温控系统:如何屏蔽外界电磁环境噪声是电磁检测和标定系统的关键。项目按照GB/T 50719-2011的技术标准建造了电磁屏蔽室,用于抑制1 KHz~10 GHz的射频干扰。在标定仓内部设计有四层屏蔽体,其中最内层采用1J85坡膜合金作为屏蔽材料,用于抑制低频电磁干扰。
TEM(TEMC104)传感器主要技术指标:频带范围为1 Hz~250 kHz,灵敏度为 5.952 mV/nT@1 kHz,静态噪声水平 9.44×10-4nT/√(Hz)@1 kHz;
CSAMT/AMT(MC-30)传感器的主要技术指标:频带范围为0。01 Hz~10 kHz,静态噪声水平在1 Hz时为0.3 pT/√(Hz),10Hz时为0.03 pT/√(Hz),100 Hz时为10 fT/√(Hz),1kHz时为3 fT/√(Hz),灵敏度为215 mV/nT;
MT(MC-50)传感器主要技术指标:频带范围为 0.0001 Hz~360 Hz,静态噪声水平在0.01 Hz时为10 pT/√(Hz),1Hz时为0.05pT/√(Hz),10Hz时为10fT/√(Hz),100Hz时为0.5fT/√(Hz),灵敏度为1359 mV/nT。
研制的三种类型的电磁传感器与国外同类型传感器进行了室内及野外的对比试验,从对比试验的结果看,项目研制的三种类型的传感器技术性能已经达到国外同类型产品的水平。
在高灵敏度传感器研发过程中,还建立了一套电磁传感器的测评系统和测评方法。测评系统主要包括:磁屏蔽室、温度控制系统、交变磁场复现装置、标定仓及测控系统等部分。其中测控系统分别设置有8个高速通道和8个低速通道,可以实现被检测传感器不同频带范围内数据的采集工作。由于系统内部设置有温度控制系统,可以研究环境温度对电磁传感器技术性能的影响效应。
感应式电磁传感器检测与标定技术关键是磁屏蔽技术及温控系统,按照 GB/T 50719-2011 的技术标准项目建造了电磁屏蔽室,用于抑制1 kHz~10 GHz的射频干扰。在标定仓内部设计有四层屏蔽体,最内层采用某种坡膜合金作为屏蔽材料,用于抑制低频电磁干扰。温控系统采用水浴循环方法实现标定仓内范围为-40 ℃~+80 ℃温度的调节,该方式的最大优点在于能够在标定仓内部创建均匀的温度场。
该系统可以实现MT、CSAMT/AMT以及TEM电磁传感器的主要技术指标的测评工作。为实现测试过程的标准化,项目还制定了传感器主要技术指标的测试流程以及数据处理方法等,为感应式电磁传感技术指标的测评提供重要的技术支撑。
广域电磁探测系统由大功率发射系统和电磁信号采集系统构成。主要包括:大功率伪随机电流信号发送机、伪随机信号控制器、多分量电磁数据采集站和数据采集与控制软件等。发电机组交流输出经调压整流后,在GPS与伪随机信号控制器的控制下向接地导线或不接地回线发送伪随机电流,并全波形记录电流数据;在数据采集软件控制下,二通道或五通道数据采集站全波形、高精度、大动态范围地记录大地对发送电流的电磁场响应信号。通过对电流数据与接收的电磁数据进行去噪等处理、广域电阻率计算和反演,获得地下一定深度视电阻率。
伪随机信号控制器:伪随机信号控制器是系统的重要部件,要求具有人机交互功能、数字信号传输与处理功能以及少量的模拟量处理能力。控制器硬件结构以A2F500为核心,外围再增加电源电路、USB通信、LAN通信、信号驱动电路,以及人机交互的LCD显示屏、键盘接口电路,和LED指示灯电路。A2F500内部集成了cortexM3硬核和50万门的FPGA,具有安全性、稳定性、低功耗等特点。
200 kW大功率电磁信号发送机:发送机总体结构包括主回路和控制回路。主回路将750 V 50 Hz的三相交流电通过输入过流保护和主接触器输送到由二极管全桥构成的三相整流器,再送往滤波组件变为1000 V的平滑直流电。直流高压输往IGBT全桥逆变模块,变换为含有信号波形的输出功率信号,通过输出电缆将此信号送到发送电极AB端(图5)。控制回路包括控制电源系统、控制开关、监控电路、显示电路,以及IGBT驱动电路等单元。
图5 200 kW电磁信号发送机总体结构
叠层母排技术:为了减弱寄生电感对功率开关IGBT的危害,减小发射机体积,增强发射机的长时间工作可靠性与稳定性,发射机中将主回路,包括整流、滤波、逆变等组件,采用叠层母排技术,将其融为一个整体。叠层母排具有低成本、安全、可靠和节省空间等优点;叠层母排本身的电感和阻抗比较低,与传统电缆相比,具有更好的载流能力,温度升高较慢。
外观结构设计:野外面临灰尘、振动、潮湿等复杂环境。外观和内部结构设计上充分考虑了以下因素:内部模块承重、散热、防震、防水、防潮、防盐雾等,按照相关标准考虑安全性和电磁兼容性,以及便于零件安装和维修检测等。
电磁数据采集站:数据采集站接收电场信号或磁场信号,经过超低噪声前置放大、低通滤波、程控放大等信号处理,然后进行24位ADC数字化;采用高可靠FPGA系统控制ADC的转换,并将多道电磁数据存储或传送给ARM处理器进行存储。采集站硬件电路设计包括三部分:模拟电路、数字电路和GPS电源三大部分,每部分零部件的选型充分考虑了先进性、可靠性和适用性。
经第三方测试,伪随机广域电磁发射机和接收机技术指标如下:
发射机技术指标:最高发射电压:1025 V;最大发射电流:205 A;最大发射功率:201 kW;测试期间同步误差变化范围:18~25 ns;发送机可靠性:满负荷工作5h无异常;频率稳定度Fs :2.12×10-8。
接收机主要技术指标:噪声:0.1 Hz~10 Hz各通道最大RMS值为0.208 μV;动态范围:最小值为122.3 dB;一致性:误差均方值为:0.17%;输入阻抗:最小值为5.84 MΩ;精度:各个通道精度最大值为:0.39%(RMS);接收机带宽:0.01~28 kHz。50 Hz陷波指标:最大陷波深度为:-66.73 dB;环境温度:-40 ℃,且在常温与-40 ℃环境温度下一致性误差最大为1.04%。
项目成功研制出实用化、产品化的200 kW广域电磁探测系统(图6),技术指标总体达到和优于国外同类型产品,实现了小批量生产能力。并在多个矿区进行了应用,取得了较好应用效果。
大功率坑-井-地电磁成像系统由发射系统、接收系统及其辅助设备组成。发射系统主要包括:大功率发电机、大功率可调变压整流器、大功率电磁发射机,以及发射控制器等。接收系统根据接收的位置可以分为地面、坑道和井中单元、发射电流记录仪等。地面单元主要包括数据记录仪、磁传感器、不极化电极等。坑道和井中接收单元分别根据空间限制采用不同的传感器和配套的设备。系统可以利用已有的坑道、矿井、井孔空间,实现地面和地下准三维观测,获取地下介质导电性、极化特性等多种参数,并进行综合地质解释。系统突破了传统地表观测或者井中观测方式,支持地面发射,井中、或坑道、或地面接收方式;井中发射,井中、或坑道、或地面接收等多种工作模式。发射系统支持的发送制式包括:可控源音频大地电磁测深法、时间域激电法、频率域激电法、多频、伪随机等,同时也支持频率域任意频率的发射。
①发射机与发射控制器:发射机主要包括:全桥整流、支撑电容、大功率逆变模块、散热系统、状态监测系统等(图7a)。外部的高压通过高压输入接口进入发射机,经过全桥整流单元和支撑电容后给大功率逆变模块提供电能,然后通过高压输出接口将逆变产生的多制式矩形人工源信号发送出去;状态监测系统采集发射机内部的温度、电流和电压信号,并通过RS485总线传送到发射控制器,显示发射机工作状态,同时控制大功率逆变模块的信号发送。
图7 井-地电磁系统发射机和发射控制器原理结构(a)、(b)及实物图(c)、(d)
发射机控制器核心主控采用STM32F103单片机与复杂可编程逻辑器件(CPLD)组合而成,主要电路包括时钟、实时时钟(RTC)、无线通信、电源管理、数据存储和驱动放大电路,以及与发射机通信接口等。发射机通过无线Wi-Fi与上位机进行数据交互,主控通过读取上位机的控制指令,将指令中的频率、发送时间长度、开启发射时间等信息提取出来,并用于发射机的发射配置(图7b)。
② 电流记录仪:电流记录仪主要包括电流传感器、ADC电路、主控FGPA+ARM核心电路、时钟、电源以及外围接口电路等。研制的关键技术是电流幅值的精度、时间精度、全波形记录等。采用与数据记录仪相同的软硬件平台,保证了较高的数据吞吐率;在电流传感器、通道动态范围方面,实现了电流记录器的低噪声及大动态范围特性;高精度时间方案与数据记录仪方案一致。
③ 地面-坑道电磁接收机:接收机内部主要有采集电路、电池组、机壳、接插件等组成,采集电路包括前端接口板、磁场通道板、电场通道板、FPGA时钟逻辑板、ARM控制板、时钟板与电源板等7部分。外围接口包括GPS天线、3通道磁场、3通道电场、LED指示灯、软开关按钮和电池输出、输入口等。
④ 井中电磁接收机:井中接收机包括井口控制单元、绞车和电缆、井中探管等部分组成。井口控制单元实现绞车电缆收放控制、电缆长度记录、与井中探管通讯、数据显示、计算、存储等工作;井中探管实现四分量电磁场信号的传感、采集、存储、传输。
坑道环境下进行磁场信号观测时,由于空间受限、磁棒方向难以确定,采用微型三分量音频磁传感器。它将三轴感应线圈集成至一个边长为30 cm的立方箱体,三轴线圈及放大电路保持相对独立,三轴线圈相比磁棒现场作业有更好的正交度。
经第三方测试,大功率井-地电磁成像系统的技术指标如下:
发射机技术指标:最大发射电流:60 A;最大发射功率:48 kW;发射机时间同步精度:初始误差:±0.1 μs,时钟漂移:20 μs/hr;发送机可靠性:满负荷工作8h无异常;发射波形频率:10 kHz~0.01 Hz(IP:100~0.01 Hz)。</