下载文档原格式
第28卷 第4期物探化探计算技术 2006年11月收稿日期:2005-04-29文章编号:1001—1749(2006)04—0337—05可控源音频大地电磁测深法应用实例黄力军,陆桂福,刘瑞德(中国地质科学院地球物理地球化学勘查研究所,廊坊065000)摘 要:目前对电法的探测深度要求越来越高,常规电法有时难以满足实际工作需要。
与常规电法相比,可控源音频大地电磁测深法具有探测深度大、设备相对轻便和横向分辨率高等特点,已广泛应用于地热资源、水文地质、有色金属矿产和工程地质勘查工作之中,并取得了令人十分满意的实际应用效果。
这里介绍的是一些具有代表性可控源音频大地电磁测深应用实例和应用成果,这些实例基本函概了目前电法的主要应用范畴,常规电法难以胜任的地质调查工作。
实例中包括断陷盆地型地热资源和深部基岩裂隙水调查,危机矿山深部找矿和复杂地形条件下铁路隧道工程勘查。
关键词:电磁测深;地热资源;水文地质;有色金属矿产;工程地质;调查中图分类号:P 63113+25 文献标识码:A0 前言可控源音频大地电磁测深是在音频大地电磁测深(AMT )基础上,发展起来的一种人工源频率测深方法。
上世纪五十年代,在卡尼亚(L Cagn 2iard )论文[1]的基础上,发展形成了基于观测天然场大地电场和磁场正交分量,计算视电阻率的大地电磁测深方法。
在音频(n ×10-1Hz ～n ×103Hz )范围内,大地电磁场相对较弱,人文干扰较大。
为了克服上述困难,七十年代初,D W Strang way 教授和他的学生M A Goldstein [2,3]提出沿用AMT 的测量方式,观测人工供电产生的音频电磁场。
由于所观测的电磁场频率、场强和方向均由人工控制,其观测的方式与AMT 相同,所以称其为可控源音频大地电磁测深(简称CS AMT 法)。
CS AMT 法可以采用磁性源或电性源二种人工场源,目前主要应用电性源可控源音频大地电磁测深法。
可控源音频大地电磁测深在地热资源中的应用可控源音频大地电磁测深在地热资源中的应用引言多年实践表明,地热资源是一种十分宝贵的综合性矿产资源,其功能多,用途广。
地热资源的综合开发利用,在社会、经济和环境效益均很显著,在发展国民经济中已显示出越来越重要的作用。
一般来说地热资源埋藏深,开采风险大。
为了减少开采风险,提高效率,开发地热资源必须进行地质调查,地球物理勘查是地热资源勘查的重要方法之一。
随着地热资源开发的难度越来越大,深度也越来越深,这就要求我们寻找更有效方法。
我们通过CSAMT法在某地的深部地热勘查中取得较好的地质效果。
1 方法简介CSAMT法是上世纪八十年代兴起的一种物探方法。
该方法根据电磁感应的趋肤效应,高频的电磁场穿透深度浅、低频电磁穿透深度深的原因,随着频率的改变,探测深度也随之改变[1]。
从电磁波的趋肤效应理论分析可得到趋肤深度公式：H≈356■（1）式中：H为探测深度,？籽为电阻率,f为频率。
2 技术措施CSAMT法具有信噪比高,观测信号强,设备较为轻便,生产效率高等优点。
这符合我们进行地热勘查的要求。
本次工作所使用的V8型接收机系统出队前经过了认真检查,符合有关规定,均可在野外施工使用。
在野外数据采集前与外业工作结束后均按技术要求对设备进行了标定。
经检查,标定的相位谱和电阻率谱曲线规则、光滑,符合《CSAMT法技术规定》要求。
①供电系统：电源为30kega;,个别点位因地表条件较差,接地电阻在4000Ω左右,采集频率同发射频率。
主要采集技术参数为：每测点按设计工作频率表依次扫频观测,每个频率叠加次数按采集时间满足1分钟或叠加60次。
自动记录、存储。
③同步系统：发射与接收采用GPS卫星时间同步,精度为0.1ns。
④观测方式：CSAMT采用赤道偶极排列。
3 实例分析该工区平行布置了4条剖面,各剖面异常形态与电性结构具相似性。
4条剖面较好地反映了不同地电特征的岩性界面。
可控源音频大地电磁法应用于地热资源勘查可控源音频大地电磁法应用于地热资源勘查(3)图4 花家地(J-76)孔旁测深柱状图理论上一维层状模型反演的前提条件是:地电模型基本满足一维地电条件，初始参数正确，反演解释能够符合实际情况。
在实际工作中，有时地质断面不能满足一维条件，这时的解很难符合实际情况。
在CSAMT没有二维反演解释方法的情况下，采用上述两种反演相结合的解释方法，对提高解释准确性会有帮助。
在CSAMT测深曲线解释中，K型和H型曲线电性分层相对较为容易，A型和Q型曲线解释相对较为困难，这是指划分电性层的难易而言。
无论是什么类型曲线将电性分层进行地质分层解释时的难度都是比较大，特别是新区难度更大。
如怀柔庙城地热井和大兴庞各庄地热井均为明显的K型曲线，电性分层没有问题，但是在地质解释上都出了大问题。
实钻表明前者将蓟县系铁岭组解释为奥陶系、寒武系地层(误差达1520m，286%);后者将蓟县系雾迷山组地层解释为寒武系(误差1800m，450%)。
前者是附近有一眼旧井的误导，后者是属于不够慎重(前人只推有Qn，这次又加上?，错上加错)。
产生问题的根源是物探的地质解释多解性。
O、?及Jx都可以是中高阻电性层。
解释为其中的哪一个一定要慎重，充分搜集已知孔的资料(准确的)，由已知推未知，做出正确的判断。
这两口地热井实施结果表明:地质分层出了问题(张冠李戴)，不仅仅是热储层埋深误差大小问题，更为严重的是推断的盖层条件(地层和厚度)与实际出入很大，直接影响温度预测，出现温度风险，给工作造成了很大的被动。
目前庞各庄地热井采取的补救措施是将热储目的层由蓟县系雾迷山组改为长城系高于庄组。
从北京统计的一部分地热井的测井电阻率看，电性层与地层不是完全一一对应的。
有时在同一个时代地层有两个电性层或多个电性层，如K或J中的电性层(有时是系中组的反映)，也1有一个电性层是由跨两个时代地层的同种岩性组成(如?+Q)。
因此，地质分层解释是一个1cnj非常细致和复杂的工作，一定要倍加慎重。
大地电磁测深数据处理及应用研究的开题报告
一、研究背景和研究意义
大地电磁测深技术在石油勘探、地质灾害预测、水文地质和环境地
球物理等领域已经得到了广泛的应用。
该技术可以获得地下多层结构的
电阻率信息，是研究地下细节结构和探测深部资源的重要手段。
近年来，随着该技术的发展，采集到的数据量也越来越大，数据处理和解释的难
度也随之增加。
因此，对大地电磁测深数据的处理和应用研究具有重要
的现实意义和科学价值。
二、研究内容和研究方法
本研究主要包括以下内容：
1. 大地电磁测深数据预处理：包括数据质量控制、数据去噪和数据
重采样等预处理步骤，以提高数据质量和处理效率。
2. 电阻率反演和成像：采用多种反演算法对大地电磁测深数据进行
电阻率反演和成像，生成地下电阻率分布图像。
3. 基于电阻率图像的地质解释和资源定位：将电阻率图像与地质信
息结合，进行地质解释和资源定位研究。
本研究主要采用数值模拟和实测数据分析相结合的方法，通过MATLAB编程实现数据处理和反演算法。
三、研究进展和计划
目前，已完成大地电磁测深数据的预处理工作，包括数据质量控制
和噪声去除等方面。
下一步计划是进行电阻率反演和成像，并将其与地
质信息相结合，进行地质解释和资源定位研究。
同时，对比分析不同反
演算法的效果，寻求最优解，并对研究结果进行验证和评估。
四、结论
本研究将有助于进一步提高大地电磁测深技术的应用效率和数据处理精度。
同时，该研究成果可以为石油勘探、地质灾害预测、水文地质和环境地球物理等领域的研究提供有价值的参考和支持。
1．大地电磁测深应用范围、面临问题及解决方法大地电磁测深（MT ）已成功的应用于油气田勘探、矿产资源勘查、地热资源调查、工程勘察、地壳和上地幔深部地质构造的研究中，取得了明显的地质效果，已成为这些领域的主要应用手段之一。
但由于天然电磁场某些频段振幅弱（图1）、某些地区工频信号及谐波干扰又很强（图2），导致信噪比低，使某些频点的视电阻率和相位值失真。
另一方面，在传统的MT 方法中都是在频带的每个量级中取8-12个目标频率计算视电阻率和相位值，然后将其连成曲线。
由于频点稀可能使某些薄层被圆滑掉了，降低了MT 法的分辨能力。
德国Metronix 公司研发的GMS-07e 综合大地电磁仪所携带的Mapros 数据处理软件成功的解决了上述的两个难题，显著提高了MT 法的应用效果。
`图2测点频谱叠加图50Hz 工频信号 工频谐波干扰 图1 测点频谱叠加图振幅较弱，数据易受干扰2．抑制工频信号及谐波干扰德国metronix 公司研发的综合大地电磁仪GMS-07e 是目前最先进的频率域电磁法勘探仪器，它采集的是天然电磁场时间域信号，而反映地下电性结构的视电阻率和相位值是在频率域中求得的，因此需通过傅里叶变换将时间域信号转换成频率域信号。
若数据采集时存在工频干扰，它也会记录到采集的数据中。
由于工频信号不满足平面电磁波垂直入射的条件，导致视电阻率和相位值失真。
为了抑制工频干扰，在Mapros 预处理软件中可调节目标频率的窗口宽度，避开工频信号及其谐波对数据的干扰，可极大程度的提高数据信噪比，从而提高观测精度。
Mapros 不仅可以通过筛选时间序列、改变目标频率窗口宽度避开噪声频率，而且也可以自定义输出频点个数提高大地电磁法的分辨率。
现举例如下：为了比较目标频率窗口宽度的选择对数据质量的影响，在其他处理参数相同的前提下，对同一频谱数据采用不同窗口宽度计算视电阻率和相位值，其中窗口宽度公式为2f t C r /π，f t 为目标频率，C r 为parzen 半径。
可控源音频大地电磁测深法在双尖山矿区勘探中的应用1. 引言1.1 背景介绍双尖山矿区位于中国华北地区，地处河北省境内，是一个潜在的矿产资源富集区。
该区域地质构造复杂，受多次构造运动的影响，形成了多个矿体赋存的地质背景。
由于地下深部构造复杂，传统的地质勘探方法在该区域已经难以满足勘探需求。
为了更好地解决双尖山矿区的地质勘探难题，研究人员开始尝试应用可控源音频大地电磁测深法。
这是一种以高频电磁波为信号源的深部地球物理勘探方法，能够有效地穿透地下覆盖层，获取更加准确的地下构造信息。
通过对矿区进行可控源音频大地电磁测深法的应用实践，可以更好地揭示地下矿体的位置和规模，为矿产资源的开发提供科学依据。
本研究旨在探讨可控源音频大地电磁测深法在双尖山矿区的应用效果，为该区域的地质勘探工作提供有力支撑，为资源的开发利用提供科学依据。
通过本研究，希望能够为类似地质背景下的矿产资源勘探提供新思路和方法。
1.2 研究意义双尖山矿区是一个重要的矿业资源区域，地下矿产资源潜力巨大。
由于地下地质情况较为复杂，传统的勘探方法往往难以满足勘探需求。
研究如何更有效地进行矿区勘探具有重要的意义。
通过本次研究，我们将不仅可以验证可控源音频大地电磁测深法在双尖山矿区的适用性和有效性，也有望为矿区勘探提供新的思路和方法。
研究结果将对该地区的矿产资源开发和利用具有积极的促进作用，为地质勘探技术的发展和矿产资源的综合利用提供有力支撑。
本研究具有重要的理论和实践意义。
1.3 研究目的研究目的是为了深入探究可控源音频大地电磁测深法在双尖山矿区勘探中的应用效果和优势，验证其在矿区勘探中的可行性和有效性。
通过对矿区内地下结构和矿藏分布进行精准探测，为矿区的资源开发和利用提供科学依据和技术支持。
通过研究可控源音频大地电磁测深法的应用实践，总结其在实际勘探中遇到的问题和挑战，为未来在矿区勘探中更好地应用该技术提供经验和参考。
最终的目的是全面评估该方法在双尖山矿区勘探中的效益和应用价值，为矿区勘探工作提供科学依据和技术支持。
能可控源音频大地电磁法在深部矿体定位中的应用电法在金属矿产物探勘查工作中一直占有重要地位，随着找矿工作的进一步深入，要求电法在寻找隐伏金属矿床方面发挥更大作用。
实际应用结果表明，在地质调查工作中，合理地使用电法，在确定成矿有利地段和隐伏矿定位预测等方面起到其它方法无法替代的重要作用。
这里的电法定位预测实际上是指在详查找矿阶段，采用以电测深为主的方法手段对矿（化）体进行二维或三维空间定位，根据其产状、规模和延伸（深）情况，结合已知地质和其它资料进行合理地质解释的一种找矿方法技术。
本文主要以某典型矿床的勘查研究成果为基础，介绍可控源音频大地电磁测深方法在寻找隐伏矿方面的定位预测作用和实际应用效果。
标签：可控源音频大地电磁法;黄狮涝金矿;寻找隐伏金属矿;定位预测一、前言为了进一步对某典型矿床进行深部矿体定位预测，采用可控源音频大地电磁法（CSAMT），在研究工作区53、56勘探线进行了CSAMT测量，目标是大致查明实测剖面内异常地质结构或地质体（矿体）的电性特征，预测隐伏―深部矿的赋存有利部位（靶区），为深部找矿和工程验证钻孔提供依据。
文中着重介绍可控源音频大地电磁测深的实际应用效果。
二、岩石电性特征和解释依据测线附近无物性资料参考，只能结合邻近区域的岩石物性资料进行分析。
我们可以通过对电磁测深数据的分析，推断测线下的地层和岩矿石分布。
断层判断依据：由于断层造成了两边地层的错动，因此可以通过电性不连续，来判断断层的存在。
再者，断层活动使断层附近可能存在破碎带，这样较陡立的低阻带也可能是断层存在的表现。
由于地质构造是复杂的，可能有多种多样的表现形式;不同岩性的电阻率也多有重叠部分。
因此电磁测深的推断解释，往往存在多解性。
故此在进行推断解释时，必须结合已知的地质和物探资料进行分析类比，反复推敲。
三、方法技术1、方法原理可控源音频大地电磁测深法（简称CSAMT法）是以有限长接地电偶极子为场源，在距偶极中心一定距离处同时观测电、磁场参数的一种电磁测深方法。
音频大地电磁测深法在某隧道中的应用摘要:本文主要讲述了音频大地电磁测深在某隧道的实际应用，从野外的施工到室内的数据处理和资料成图解释。
依据音频大地电磁法（AMT）的带地形二维非线性共轭梯度（NLCG）反演结果，把电阻率异常分为吁类、郁类、芋类及域类，分别对应极破碎岩体（吁级围岩）、破碎岩体（郁级围岩）、较破碎岩体（芋级围岩）和完整岩体（域级围岩）。
结合围岩分级结果，查明了隧道洞深附近岩性、断层破碎带及岩溶位置，对岩体完整性进行了评估。
Abstract院This article mainly introduces the application of audio-frequency magnetotellurics in a tunnel, from wild construction toindoor data processing and data mapping. According to the inversion results of the topography based two-dimensional nonlinear conjugategradient (NLCG) of audio-frequency magnetotelluric (AMT), the resistivity anomaly is divided into Class吁, 郁, 芋and 域, corresponding toextremely broken rock mass (Class吁surrounding rock), fractured rock mass (Class 郁surrounding rock), relatively broken rock mass (Class芋surrounding rock ) and complete rock mass (Class域surrounding rock). In combination with the surrounding rock classification result, thetunnel nearby lithology, fault fracture zone and the karst location are found out, and the integrity of rock mass is evaluated.关键词院音频大地电磁；隧道；电阻率异常；围岩分级Key words院audio-frquency magnetotelluric；tunnel；resistivity anomaly；classification of surrounding rock中图分类号院P631；U45 文献标识码院A 文章编号院1006-4311（2014）29-0116-030 引言音频大地电磁测深法（AMT）勘测精度高，仪器轻便，费用低，所以在工程建设中应用广泛。
大地电磁资料处理中有效视电阻率的利用李爱勇;柳建新;杨生【摘要】有效视电阻率与大地电磁响应阻抗张量矩阵的模相对应,是坐标旋转不变量,在一维条件时有效视电阻率等于常规视电阻率ρa,在二维条件时等于ρTE和ρTM的几何平均,具有降维特征.因静态效应不影响有效视电阻率的曲线形态,所以可用平移法进行静校正,由于与旋转轴无关,对有效视电阻率进行二维反演可避免因极化模式判别不准确而带来的反演结果偏差.经系统的讨论,认为在大地电磁测深资料处理解释中,可充分利用有效视电阻率发挥作用.%Based on the systematical discussion of the effective apparent resistivity properties, the effective apparent resistivity is an invariable under the coordinate rotation, which is correspondent to the modulus of magnetotelluric response impedance tensor matrix. Under the 1-D condition, It is equal to the normal apparent resistivity ρa, and under 2-D condition, it is the average of the ρTE andρTM, which has the dimension reduction property. The parallel moving algorithm was used to do static correction, as the static effect can not affect its curl shape. At the same time, it is irrelative to coordinate rotation, so the inversed deviation because of the incorrect mode discrimination, can be avoided in the effective apparent resistivity 2-D inversion. So the utilization of effective apparent resistivity in processing and interpretation of the magnetotelluric sounding data was recommended.【期刊名称】《物探化探计算技术》【年(卷),期】2011(033)005【总页数】5页(P496-500)【关键词】大地电磁测深;有效视电阻率;旋转不变量;降维特征【作者】李爱勇;柳建新;杨生【作者单位】中南大学信息物理工程学院,湖南长沙410083;江苏省有色金属华东地质勘查局八一四队,江苏镇江212005;中南大学信息物理工程学院,湖南长沙410083;中南大学信息物理工程学院,湖南长沙410083;有色金属矿产地质调查中心,北京100073【正文语种】中文【中图分类】P631.3+25大地电磁测深法(简称MT)野外数据采集，随着仪器硬件系统性能的提高，设计思路的完善，特别是将GPS应用到采集信号的时间同步上，使远参考技术从设想变为现实，大大改善了野外数据的采集质量［1］。
EMTF在宽频大地电磁测深数据处理中的应用徐坤;叶高峰;魏文博;郭泽秋【摘要】The problems as to how to suppress various kinds of disturbance efficiently, raise signal-to-noise ratio and guarantee the data quality have become the key issues in the collection and handling of MT data. To get more accurate results, the authors introduced the EMTF system ( single, reference and MMT) and used this system to process some of the MT data obtained from North China, and also compiled Matlab program so as to rewrite parameters easily. In comparison with SSMT-2000, the EMTF system can suppress the noise in most cases and greatly improve the efficiency of data processing.%如何有效地压制各种电磁噪声，提高信噪比，保证大地电磁测深数据质量，始终是大地电磁资料采集与处理的核心问题。
为寻求更好的处理效果，使用EMTF系统（单点，远参考以及多站同步MMT）对华北多个宽频MT测点数据进行了处理，并编写Matlab接口来方便地改写处理参数。
与SSMT-2000对比，结果表明，EMTF系统在多数情况下可以压制噪声，并获得较好的处理效果，提高数据处理的效率。
【期刊名称】《物探与化探》【年(卷),期】2014(000)005【总页数】5页(P1055-1059)【关键词】大地电磁测深;EMTF;SSMT2000;Robust;MMT【作者】徐坤;叶高峰;魏文博;郭泽秋【作者单位】中国地质大学地下信息探测技术与仪器教育部重点实验室，北京100083; 中国地质大学地球物理与信息技术学院，北京 100083;中国地质大学地下信息探测技术与仪器教育部重点实验室，北京 100083; 中国地质大学地球物理与信息技术学院，北京 100083; 地质过程与矿产资源国家重点实验室，北京100083;中国地质大学地下信息探测技术与仪器教育部重点实验室，北京 100083; 中国地质大学地球物理与信息技术学院，北京 100083; 地质过程与矿产资源国家重点实验室，北京 100083;中国地质大学地下信息探测技术与仪器教育部重点实验室，北京 100083; 中国地质大学地球物理与信息技术学院，北京 100083【正文语种】中文【中图分类】P631大地电磁测深法（MT）是研究地质构造的一种地球物理方法，它以天然交变电磁场为场源，野外施工简捷，成本低廉，具有探测深度大、不受高阻层屏蔽影响、对低阻有较高的分辨能力等优点［1］，已成为油气田和煤田普查勘探、岩石圈结构探测、地下水和地热资源寻找、工程勘探和地震预报等领域的一种重要手段，尤其在地震勘探方法难以施工或无法取得有效信息的地区，大地电磁是一种很好的替代方法。
但是，在大地电磁测深阻抗张量计算过程中，当自功率谱中包含噪声时，计算的阻抗有误差；互功率谱［2］、Robust［3－6］等方法在压制电磁噪声、改善数据质量等方面有一定的效果，但仅能消除不相关噪声的干扰，对于相关噪声的影响是无能为力的。
T．D．Gamble等人提出的远参考道［7］方法不仅可以消除不相关噪声的影响，还可以同时消除本地的相关噪音，数据处理质量大大提高。
但是，远参考处理并不是处处有效的，对于互参考的两点，当干扰均较强时，处理效果就会不尽人意［8］；并且在进行大面积的勘探时，往往会有多站点同时采集，并不能保证至少有一个是无干扰的测点，而互做远参考的方法无疑降低了计算的效率。
EMTF［9－10］由 Gary Egbert开发，涵盖了单点、远参考以及多站同步处理。
它的提出基于多元统计思想，消除了各道中“飞点”对整体数据的影响，进而达到压制相干噪声的效果。
笔者使用EMTF系统并编写Matlab接口对华北多个MT测点进行单点、远参考、多点同步处理，并和国内广泛使用的与宽频仪器（V5⁃2000）配套的 SSMT⁃2000处理系统进行对比，结果表明，EMTF处理系统压制了噪声，结果真实可靠。
1 EMTF系统简介1．1 EMTF 处理步骤在EMTF中，转换函数（TF）估计程序的使用可以分为三个（或两个）步骤。
（1）重定数据格式。
通常，数据文件有不同的格式来源，所以要想开发出一个普适于所有可能数据格式的处理软件是很困难的。
对时间序列的数据文件，采用一种虽然并不常见，然而却是标准的空间高效格式来表示。
一旦数据转换成这种格式，那么，在第二阶段——傅里叶变换阶段中，dnff程序就很容易运行了。
EMTF系统含有能将两种数据格式分别转化成这种标准格式的转换程序rfemi和rfasc。
程序rfemi是针对EMI MT⁃1系统的，而程序rfasc则是针对普通的ASCII文件的。
对于几种特定的数据格式，dnff程序可以直接在原始数据上运行而不需要对他们进行格式转换。
这些特定的格式包括：普通的 ASCII格式，EMI MT⁃24格式以及mini⁃SEED文件格式。
（2）对数据段的开窗和傅里叶变换。
这是由程序dnff来完成的，输入该程序，包括一个数据文件（既可以是标准的二进制格式也可以是普通的ASCII格式）、一个描述系统参数的文件及一些控制程序选项，还包括一个可选的数据段剔除指令文件。
该程序的输出是一个由按频率顺序排列的傅里叶系数（FCs）组成的文件。
程序dnff的关键特征包括一个内嵌分样表和一些由不同测点得到的矫正时间窗特征，这些特征允许在不同的起始时间或采样率的情况下对独立采集点进行远参考及多点处理。
输出的傅里叶系数文件在随后所有的频域处理程序中都要用到，这些程序包括单点及远参考转换函数估计程序tranmt和多点处理程序multmtrn。
（3）对转换函数的robustROBUST估计。
这是由程序tranmt来完成的（或者是多点处理程序mult⁃mtrn来完成）。
该程序的输入是有序的傅里叶系数文件以及一些控制程序选项的文件。
1．2 EMTF 目录结构EMTF目录结构见图1，说明如下。
图1 EMTF目录RF／＝＝＞数据格式重定程序RFEMI＝＝＞将EMI MT－1格式转换成标准的二进制格式RFASC＝＝＞将普通的ASCII格式转换为标准的二进制格式BIN2ASC＝＝＞将标准的二进制数据文件转换成普通的ASCII格式FCOR＝＝＞由一个标准的系统参数（．sp）文件算出一个滤波器／系统响应校正检索表（这种计算还需要传感器文件）C／＝＝＞数据清除（平滑）程序D／＝＝＞分样／快速傅里叶变换程序T／＝＝＞单点以及远参考ROBUST转换函数估计程序MMT／＝＝＞多点ROBUST转换函数估计程序include＝＝＞包含文件，这些文件为好几个目录下的程序所共用matlab／＝＝＞一些matlab绘图及后期处理程序test／＝＝＞一些用于测试的简单（人造的）数据样例DATA／＝＝＞包含两个压缩的普通 ASCII文件 test1．asc．Z和test2．asc．Z，以及一个时钟文件SP／＝＝＞系统参数文件样例CF／＝＝＞配置文件样例，以及测试运行需要的所有文件FC／＝＝＞空目录，该目录用于存放得到的傅里叶系数，这些系数后面带有说明MT＝＝＞空目录，该目录用于存放程序tranmt的输出文件MMT＝＝＞空目录，该目录用于存放程序multmtrn的输出文件sensors＝＝＞空目录，该目录用于存放系统参数检索表（如果该检索表存在），或者用一个叫sensors的链接取代它，该链接用来定位系统检索表的存放路径results＝＝＞当测试的每一步完成之后出现在 D，FC，MT和MMT目录中的结果bin＝＝＞空目录，该目录中将安置用于测试的可执行对象EMI－MT1＝＝＞用于EMI－MT1数据格式重定的测试数据doc／＝＝＞程序文档TEX＝＝＞程序文档的latex格式源文件ps＝＝＞postscript格式的程序文档文件从上可以看出，EMTF系统需要手动修改许多参数配置文件，也不能实现多测点数据的批处理功能，这样处理过程会非常繁琐，数据处理效率非常低。
为此，笔者尝试编写Matlab接口来整合参数配置文件及可执行程序，视窗式的操作实现了参数配置文件的修改和核心程序的后台调用，并可绘制响应曲线以供分析，提高了数据处理的效率。
2 接口处理实测数据流程本次华北数据的野外采集使用凤凰公司的MTU⁃5A／P（V5⁃2000 格式）仪器，分时、分段采集，得到TBL参数文件和TS∗数据文件，多站点每天同时采集，并保证具有8 h以上的共同数据段。
编写的Matlab接口如图2所示。
其数据处理流程如下。
（1）转换TS∗文件的格式为ASC文件。
（2）去除记录信息，截取采集的数据，得到数据ASC文件以及时钟clk文件。
由于该程序仅仅要求是时间序列文件，因此必须将数据ASC文件中的记录参数去除，拼接成完整的时间序列（对不同采样率文件分别截取）。
为此编写了Matlab接口来实现这一功能，并得到时钟clk文件。
对于远参考和多点同步处理需要截取共同时间段采集的数据。
以下是clk文件示例。
16．0 ＜＝＝＝采样率，单位用秒。
85 07 20 0 0 0 ＜＝＝＝时钟设置：年，月，日，时，钟，秒。
85 07 20 0 0 0 ＜＝＝＝通用时钟零时间：年，月，日，时，分，秒。
图2 matlab接口图示（3）傅里叶变换。
这一过程包括了格值转换、数据预白、级联分样和开窗，调用了EMTF中的dnff程序块。
总的来说，开窗和傅里叶变换中开窗长度和级联分样非常重要，选取得不好，频点的拼接就会有重叠或者是缺少某些频率段。
其基本规律是“开窗”长度大，分样水平级数低，越能得到响应的低频段信息，反之，能得到高频段的信息。
通过实验，得到了一组稳定的处理参数。
其中，TS3对应高频，开窗长度选择512，分样级数为2；TS4对应中频，开窗长度选择256，分样级数为2；TS5对应着低频，采取分两段处理的办法，第一次得到相对高频的数据，开窗长度选择256，分样级数为4，第二次是得到相对低频的数据，开窗长度选择4 096，分样级数是5。
编写的matlab接口实现了这一功能。
