TEMPROS瞬变电磁反演软件说明书

概述
TEMPROS是为TEM数据预处理制作的软件，主要适用于加拿大POTEM系列仪器和美国Zonge公司的GDP系列仪器采集的TEM数据处理。

可对Gx7文件（POTEM 系列）和avg文件(GDP系列仪器），以及用户自定义数据进行可视化数据编辑，其效果直观、快速有效。

反演方法有线框的一维快速反演，大定源的一维反演。

本软件由编程者与北京欧华联科技有限责任公司联合出品，版权归睿德所有。

一、Protem数据导入
1、在“数据导入”菜单下的“Protem->dat”，可接读Protem57的数据
“*.Gx7”文件，直接存入“*.dat”文件里。

2、*.dat文件的数据格式
点号 时间（ms） 感应电动势[uV/(A*m*m)]
160.0 0.09000 990.60000
160.0 0.11310 823.60000
160.0 0.14250 806.60000
160.0 0.18060 841.40000
160.0 0.23000 825.00000
160.0 0.29250 828.40000
160.0 0.37000 616.20000
160.0 0.46630 371.00000
160.0 0.58500 247.40000
160.0 0.72500 146.58000
160.0 0.88130 96.00000
160.0 1.06900 59.12000
160.0 1.31300 35.82000
160.0 1.61900 21.56000
160.0 2.00600 12.66200
160.0 2.50600 7.44600
160.0 3.14400 4.29200
160.0 3.95700 2.28000
160.0 4.99400 1.23080
160.0 6.31300 0.73880
160.0 7.99400 0.24440
二、曲线显示
1、打开文件，可选择*.dat文件后，点击“曲线显示”后，无数据编辑功能。

2、点击滑动键，变化不同点的衰减曲线，可显示当前点点号。

三、选择点号
通过按滑键，选择点的数据。

数据编辑手段
1、单击鼠标左键，改变单个频点的数值。

2、曲线的整体移动方法是对准需要调到的位置，点击鼠标右键，最好鼠标也放在在任意频率轴上。

四、预处理
1、预处理方法选择：数据编辑、数据截断和后延改正。

1.1 数据编辑
1.1.1 打开文件，可选择*.dat文件，进行数据编辑，此时无曲线显示功能；
1.1.2 选择“数据编辑”，进入数据编辑。

通过按滑块选择某点的曲线，通过鼠标点击曲线单点修该数值；
1.1.3 保存编辑数据
1.1.4 其它编辑功能
（1）当前是编辑单个点衰减曲线，可以点击“三点曲线”，可以参考当前点前后两个点曲线，“单点曲线”和“三点曲线”可以切换运用；
（2）如果测量观测数据为负数时，可以点击“单点取负值”对当前点所有时间道的衰减曲线数据取负值，也可根据情况对所有测点进行“所有点取负值”处理。

（3）对数据编辑后，在退出编辑时，提示保存数据，以保存数据。

2、数据剪裁
在扰大时，通过数据剪裁，剔除早期道或晚期道的数据，打开*.dat文
件后，
显示第一道时间和最后道时间，可以输入合适的时间，保存这个时间段
的数据。

五、二维小波变换，区分正常场、提取异常场
1、输入参数选择，通常和外延两种小波运算方法从TEM数据中区分纯异常场和正常场，图形显示原始数据的剖面曲线、衰减曲线，异常场和正常场的剖面曲线、衰减曲线。

2、处理结果以剖面、测点格式保存。

六、回线反演
回线反演包括回线数据转换、回线反演、结果显示、绘制surfer等值线图等功能。

1、回线数据转换
将dat数据转换成Teminv反演的数据格式。

读入*.dat 文件，生出*.inp文件。

一般经验选择：陷波＝0.01，是高频平均，数据的接收线圈面积已归一，保存*.inp 文件。

2、回线反演为中心回线反演
2.1 快速反演法，输入*.inp文件，电流已归一为1，需输入发射线圈X边和Y边长度、线圈匝数。

2.2 反演计算生成 *.iat和 *.res文件。

3、 反演成果显示
反演结束后，上图自动退出，这时选择结果显示，出现反演电阻率断面图。

4、绘制surfer等值线图
可调用surfer8，选择*.iat或*.res反演结果，直接生产surfer等值线图。

七、定源一维反演
定源反演法，读入参数文件*.fix，读入数据文件*.dat，改变各参数后进行反演，反演生成*.inv文件。

其步骤有：
1、层数、最大迭代次数；
2、回线X方向和Y方向边的长度；
3、电性层参数输入，输入层数后，显示如数个电性层参数，用鼠标左键双击后，输入数字，确定回车键，可保存参数文件（*.fix）内，以便在下次反演前可直接读入参数文件；
4、若存有fix文件，可以直接读入参数文件；
5、反演需读入数据文件*.dat文件，可直接生成反演结果文件名*.inv，若另存其它文件名，可以选“反演另存”，输入其它文件名；
6、反演
7、以下为形成反演图件的文件数据和surfer电阻率等值线图生成
（1）在生成做图数据之前，需输入成图迭代次数、距离(m)/点间隔、绘制最深深度，以及是否带有地形数据，以便将深度转为地形高程；
（2）读入的“*.inv”文件，生成第n次迭代数据，文件为*n.TXT。

（3）最后，确定“生成surfer图”键，生成surfer图。

八、实测资料试算结果
为了验证本研究的方法技术及软件的有效性，我们选取北方某铜镍矿一条实测的TEM 剖面进行了处理，得出分离后获得的异常剖面，见图2、图3所示。

2.1 工区地质情况
测区主要地层是前寒武系下统强变质岩系和前寒武系上统的浅变质岩系。

矿床是与超镁铁岩有关的硫化铜镍矿床。

岩体直接围岩主要为蛇纹石化大理岩、黑云母片麻岩、二云石英片岩和条带状混合岩等。

岩体由三次侵位的超基性岩叠加形成。

早期岩石以中细粒含二辉橄榄岩为主；中期岩石以中粗粒二辉橄榄岩为主；晚期岩石则为纯橄榄岩。

三期岩石相互间呈侵入接触，是深部分异、先后侵位的产物。

岩体主要为纯橄榄岩－二辉橄榄岩－斜长二辉橄榄岩。

2.2 电性特征
超基性岩与围岩的电阻率一般为几十Ω·ｍ，有时增至为几百Ω·ｍ，片麻岩有时呈较低值。

海绵状及块状矿石的电阻率较低，
而星点状矿石与超基性岩和岩体围岩的电
阻率差不多是同等级的。

总的来说，区内各种岩石的电阻率差别不大，采用普通的电阻率法在本区无明显的效果。

(μV/A)ε/I ε/I
2.3 资料解释：
从42线地质断面图(图4) 中可以看出：地层南倾，220点以北出露的地层岩性为蛇纹大理岩和大理岩，640点以南出露的岩性为条带状混合岩，从220点至640点间为二辉橄榄岩，其下部的矿体发育于岩体的下侧部,以星点状贫矿为主,其次为海绵损铁状富矿和块状特富矿。

图2的多测道感应电动势/电流(I/ε)剖面曲线反映了地质体的横向变化，曲线在200～600点间的数值高于两侧，与北侧差异较小，与南侧有较大幅度差异。

曲线中段的中值响应反映了中部相对为低阻的二辉橄榄岩，在200点北侧的相对于中段较弱的响应反映了北端为高阻的的大理岩，600点以南的低值响应是南侧条带状混合岩的反应。

在300点晚期道曲线有明显的高阻异常，是深部低阻体的反应，钻孔未揭露到此异常。

在600～800点深部的海绵损铁状富矿和块状特富矿在晚期道有较弱反应，但不明显。

经二维小波变换提取的TEM异常剖面如图3中的曲线，对岩性在横向变化有明显的差异。

在600点处呈现低值响应，正是二辉橄榄岩和条带状
点号
混合岩接触带的反应。

在600～800点处的晚期道出现了局部高值异常响应，很好地反映了相对低电阻的海绵损铁状富矿和块状特富矿。

图5为TEMPROS 快速反演结果，在008点深部的低阻与海绵损铁状富矿和块状特富矿相吻合。

图5 反演电阻率断面图。