ARTICLE · INTELLIGENCE

战地情报 · 详情页

来自尧图项目组的一线实战观察与深度解析

瞬变电磁数据反演全流程:从bin导出到IX1Dv3一维反演实战

瞬变电磁数据反演全流程:从bin导出到IX1Dv3一维反演实战 简介《软件培训讲义.pptx》是一份面向瞬变电磁法TEM从业者和软件初学者的操作培训材料系统梳理了瞬变电磁法基本原理、工作方式、应用场景以及IX1Dv3软件从数据导入、工区建立、数据编辑、初始模型建立、人机交互反演到成果保存的完整流程。讲义对Terratem采集数据的bin格式转换为USF格式作了重点讲解并给出主界面File、Open、Export等菜单路径可直接对照开展实际数据处理。资源包共1个文件文件类型为pptx演示文稿整体大小约4.19MBPPT形式便于随堂讲解和后期修订。目前已有63人学习。对于正在使用IX1Dv3处理TEM测深数据、希望补齐软件操作细节的地质勘查、沉积盆地研究和航空电磁填图相关技术人员这份讲义兼具原理说明和分步操作指引能帮助快速上手并减少格式转换与反演流程中的常见错误。1. 瞬变电磁数据反演为什么野外拿到 bin 数据还只是开始做过 TEM 的人都知道野外拉回来的 bin 文件只是半成品。断电后二次场电压随时间衰减的曲线并不直接等于地下电阻率要把不同延时道的响应“翻译”成随深度变化的电阻率模型核心动作是一维反演。IX1Dv3 的瞬变电磁模块恰好就是干这个的它把 Terratem 实测数据、工区管理、数据编辑和反演计算串成一条完整链路而这条链路的第一道门槛是数据格式转换。这篇笔记基于一份软件培训讲义整理把从 bin 导出 USF、建工区、编辑衰减曲线、跑反演到保存成果的全流程拆开讲延时怎么选、单位怎么避坑、最小二乘和平滑模型反演各自用在哪以及我实际处理数据时踩过的翻车点一次说清楚。适合正在跑 TEM 测深、手里有 IX1Dv3 但还没跑通全流程的从业者。2. IX1Dv3 与 TEM从二次场衰减到电阻率模型这套软件到底解决了什么2.1 瞬变电磁法的“一次场-二次场”逻辑先聊原理想明白了反演参数就不容易乱。瞬变电磁法以不接地回线磁偶源或接地线源电偶源作为场源发射电流脉冲然后断电后测量地下介质被激发出来的二次涡流场。这个二次场的本质是电磁感应定律断电瞬间地下感生涡流涡流随时间扩散并衰减接收线圈里看到的是一条电压随时间下降的响应曲线。这条曲线本身不是电阻率剖面。早、中、晚期物理含义完全不同早期二次场对应频率域中的高频成分衰减快、趋肤深度小反映的是浅部地电信息晚期对应低频成分衰减慢、趋肤深度大携带的是深部信息。所以一条完整衰减曲线里实际压缩了从浅到深的一组地电响应一维反演做的就是把这条曲线拆回成“层电阻率层厚度”的模型。理解这个“一次场-二次场”的逻辑后面设初始模型时才不会被“看起来拟合得很好但层参数完全不合理”这种结果骗过去。2.2 三种工作方式与其适用边界讲义把 TEM 的工作方式分成三类实际处理数据时选错类型会直接影响反演时的几何参数设置。第一类是同点类型包括同一回线、重叠回线和中心回线。同点组合下接收装置与探测目标体耦合最好异常响应强、形态简单、分辨能力强是精细测深的首选反演时几何参数基本不用额外调。第二类是分离回线类型收发距拉开适合普查或航空电磁填图但浅部分辨率会牺牲一点反演时收发距必须写进工区。第三类是大定源组合常见于矿产普查阶段覆盖范围大但数据质量更依赖场源位置与测点关系的准确记录。选工作方式和选反演算法是联动的同一回线数据跑最小二乘反演没问题分离回线的数据如果收发距没写对早期道反演出的浅层电阻率会明显偏掉。这里有个容易被忽略的细节如果现场用了重叠回线接收回线和发射回线的尺寸在反演模型里是两个参数导入时一定要确认工区里的回线边长和实际布设一致差一米早期道的响应就可能差出百分之十几。2.3 IX1Dv3 的模块组成与标准流程IX1Dv3 不是单为 TEM 写的它是一维直流电阻法DC、激发极化法IP、大地电磁法MT、电磁法EM测深正演和反演软件包。其中 TEM 模块的核心能力有四个多种数据格式接口、便捷的数据编辑功能、两种基本 1D 反演方法最小二乘反演、平滑模型反演、带约束的人机交互解释功能。这些能力对应到实际工作里就是“能接不同仪器的数据能手动剔坏道能跑反演能加约束调模型”。软件里还有两个容易忽视但很实用的功能一是对多地层做平顺模拟二是对反向分层模型做等价性分析。等价性分析我留到最后一章单独说反演结果不稳定时它比反复调参数管用得多。软件的标准操作流程是七步我按培训讲义整理成表格步骤动作产出物这一步常见问题1数据导入测点衰减曲线进入数据库格式不兼容、单位错2建立工区测点坐标、高程、剖面高程缺失、测点排序乱3编辑原始数据干净可用的衰减道误删有效道4建立初始模型层数、厚度、电阻率初值层数过多不收敛5人机交互反演计算拟合后的地电模型参数约束不当6修正模型符合地质认识的模型只看拟合差忽略地质合理性7保存反演成果模型文件加图件只存图不存模型这个顺序基本就是后面各章展开的顺序。前两步决定数据能不能被正确读进去中间三步决定反演能不能收敛最后一步决定成果能不能复用。很多人拿到软件先点反演结果数据单位、测点顺序都是乱的反演跑完才发现剖面完全没法看回头重做一遍反而更慢。我现在的习惯是老老实实按这个七步走每一步的产出物都确认过再往下进行。再补一点选型理由。能处理 TEM 的软件不止一个IX1Dv3 在一维测深批量处理这个场景里最突出的优势是数据管理大量测深数据放在同一个数据库里支持导入剖面、从导入的深度数据创建剖面图人机交互时能一条曲线一条曲线地盯着看。代价是它对输入格式要求严格Terratem 的 bin 文件不能直接读必须先导出成 USF。这就是下一章要解决的问题。3. 从 bin 到 USFTerratem 导出、工区建立与剖面的完整路径3.1 为什么仪器里的 bin 不能直接拖进 IX1Dv3Terratem 实测得到的原始数据格式是 bin。bin 是仪器自己的存储格式里面除了衰减曲线还带着仪器配置、标定参数和原始采集状态字段定义只有自家采集软件认识。IX1Dv3 不读 bin讲义里说得非常直白Terratem 能导出的数据格式里只有 USF 格式的文件能够用 IX1Dv3 直接处理。USF 是通用格式文件。所谓“通用”指的是多个软件约定用这种格式交换测深数据但不同仪器导出的 USF 字段并不完全一致所以导入前最好先人工确认内容。我一般会把 USF 当纯文本查看先确认测点号、延时道、电压值这三类信息都在再进 IX1Dv3。这个确认动作看着多此一举实际上能拦住后面一大半“测点消失”和“曲线断道”的问题。3.2 Export 操作路径与延时选择格式转换在 Terratem 采集软件里完成不进 IX1Dv3。具体路径分四步第一步在采集软件主界面选择 File点击 Open。默认弹出的文件过滤器是 *.bin先点 Select directory选择实测数据的存储路径。第二步在文件列表里选中要转换的 bin 文件点 OK。此时会弹出一个对话框要求选择延时delay time。第三步延时选择完成后回到主界面依次点击 File 和 Export导出 USF 格式文件。这里延时选择是整个导出过程最关键的参数不是随便选一个就完事。延时对应的是断电后开始采样的时间窗口选早了会把关断电流段的畸变数据带进去早期道电压异常饱和反演时浅层电阻率直接失真选晚了又会丢掉浅部地电信息浅层分辨能力下降。讲义里只提示了“选择延时”这个动作没展开讲怎么选我的经验是如果野外采集时用了多组延时导出时要把同一个测点的所有延时道都选上如果一组延时里明显包含关断电流段的畸变道就要在后续编辑里单独剔除。分离回线数据还要额外确认收发距与延时道一一对应不然反演时早期道对不上几何关系。提示导出的 USF 文件名建议直接用测点号命名。我见过有人导出后文件名是一串乱码式编号进 IX1Dv3 后测点识别混乱最后只能对照采集手簿手工重命名白白浪费一下午。3.3 导出后先验证USF 文件里到底有什么导出之后不要急着开 IX1Dv3先花一分钟确认文件内容。这里给一个最简单的检查方式# 读取 USF 文件前 30 行确认测点编号、延时道数量和电压列 with open(site_001.usf, r, encodingutf-8, errorsignore) as f: for idx, line in enumerate(f): if idx 30: break print(idx, line.strip()[:120])这段脚本的逻辑很简单把 USF 当纯文本逐行读出来打印头部 30 行人工确认文件里至少包含测点号、延时道列表和数据列。参数说明site_001.usf 换成实际导出文件的路径errorsignore 是为了防止个别仪器导出的文件里夹带非 UTF-8 字符导致读入中断strip() 和切片 [:120] 是为了避免打印超长行刷屏。正常情况头部 30 行里应该能看到测点号的变化。如果 30 行里全是同一个测点号说明这个文件只装了一个测点批量导入时要核对文件名与测点的映射关系。如果延时道数量与采集手簿对不上比如野外记了 28 道、文件里只有 20 道就别往下走了回到采集软件重新导出。这个检查习惯能省掉后续工区里很多莫名其妙的问题。3.4 建工区导入原始数据、高程与生成剖面USF 就绪之后进 IX1Dv3开始建工区。建工区不是简单的新建文件它决定了后期剖面的横纵坐标是否可信。先把 USF 数据导入数据库设定工区名称与坐标系。这里要特别注意单位问题不同仪器默认单位可能不一样米和英尺混用时会直接造成剖面横向尺度失真。导入完成后检查测点坐标的极值范围如果测区长度实际只有几百米坐标却显示几千基本就是单位没换算回到导入设置里重新选。高程数据建议单独导入。野外采集时测点高程如果没有写进 USF工区里的测点会落在同一水平面上生成的剖面图是平的地形与反演深度的关系就丢了。讲义里把“导入原始数据”和“导入高程数据”并列列出不是摆设。实际操作时我一般准备一个 csv 文件包含测点号和高程两列导入后相邻测点高差是否符合地形趋势一眼就能看出来。生成剖面是在数据确认无误之后再做。生成之前先看测点排序方式默认排序可能按采集时间或文件名而不是按测线里程排序一旦错了剖面图上测点首尾连接完全乱掉横向对比反演成果就无从谈起。生成剖面后回看一遍测点分布发现有跳跃的测点就回到数据列表修正排序再重新生成。4. 数据编辑与一维反演最小二乘、平滑模型与人机交互约束怎么配4.1 编辑原始数据先把坏道挑出去进了工区第一步不是建模型而是编辑原始数据。野外实测的衰减曲线总会有不干净的早期道受关断电流畸变影响、晚期道落在噪声底之上、个别测点出现明显飞点。IX1Dv3 的数据编辑模块支持直接查看单条衰减曲线逐道检查数值。我的编辑习惯是按“道”而不是按“整条曲线”来操作。先看全部测点在同一延时道的数值分布如果某一道在多个测点上一致性地跳高或跳低这通常是系统性问题比如延时选错或仪器道设置有问题要回到导出环节重新处理如果只是个别测点的某一道异常剔除该点就行。编辑时最重要的原则是少删只删确定异常的道不要因为曲线“看着不顺”就大段删删多了浅部信息就没了后期反演的浅层分辨率会明显变差。操作完成后最好把编辑结果存成一个独立的编辑记录。IX1Dv3 能保存编辑后的数据状态但不会自动记录你删了哪些道。我一般会另存一份原始未编辑的备份文件防止后面发现误删了有效道却没法恢复。4.2 最小二乘反演与平滑模型反演选型与参数IX1Dv3 提供两种基本 1D 反演方法最小二乘反演和平滑模型反演。名字看着像实际上模型取向完全不同选错方法比参数设错更致命。最小二乘反演适合层状地层模型。它的输出是有限层数、每层有明确电阻率和厚度适合地层分层明显、有钻孔或测井资料约束的场景。初始模型一般从两层半空间开始跑收敛后再逐步加层。这里容易犯的错误是一上来就设八层、十层——层数越多反演自由度越大结果发散得越快拟合差看起来低但模型形态完全失真。平滑模型反演适合电阻率渐变或没有明显分层界面的场景。它不强行分层而是让电阻率随深度平顺变化输出是一条连续曲线对应讲义里“多地层平顺模拟”的能力。它的优势是对初始模型不敏感、不容易在层边界处抖动代价是深度分辨率比层状模型低不能当作精确的分层依据。两种方法的差异整理成表项目最小二乘反演平滑模型反演输出形态层状模型层厚度加层电阻率连续电阻率-深度曲线适用场景分层明显、有先验信息电阻率渐变、无明确分层初始模型重要性高初值离谱容易不收敛低对初值不敏感迭代收敛节奏十到二十次内看拟合差变化五到十次看曲线形态稳定成果用途分层解释、钻孔对比电性趋势判断、填图迭代次数和收敛判据也要说清楚。一维反演的拟合差是观测曲线与正演曲线的均方根差。不要只看拟合差数字要同时看拟合曲线的形态早期道拟合好了但晚期道整体偏移说明浅部模型对了、深部没约束住需要调整深部初值或加约束如果拟合差一直降不下来先检查数据编辑是否干净再检查初始模型。4.3 建初始模型与人机交互约束初始模型设置是反演里最容易被当成玄学的环节实际有章可循。常见做法是两层半空间起步第一层电阻率用早期道视电阻率的平均值厚度用估算的探测深度除以层数深部电阻率用晚期道的平均值。先跑一遍模型收敛了再逐步加层。这个流程的好处是每一步都知道自己在加什么而不是一开始就把模型空间铺满让反演自己去瞎撞。约束是人机交互的关键。已知的层参数可以直接固定比如有测井资料确认浅层是高阻层就锁定这一层的电阻率范围与厚度范围不确定的参数给一个宽松范围让反演自行调整。IX1Dv3 的 TEM 模块带有约束功能用途就在这——把地质先验信息变成反演的边界而不是让反演在完全无约束的参数空间里搜索。修正模型这一步容易被忽略。反演计算出结果后要回到模型里逐层看合理性电阻率有没有出现负值、厚度有没有突破工区范围的物理边界、深部模型是不是在晚期道噪声范围内跳变。如果这些方面有问题就调整约束范围重新计算。反演的目的不是拟合差最小而是得到一个既拟合数据又符合地质认识的模型这两点缺一个成果都站不住。5. 常见问题与排查六个让反演翻车的实测坑5.1 测点数量与实测对不上现象导入工区后只显示 80 个测点野外明明采了 90 个。原因导出 USF 时选错了目录或者选中了部分文件。Terratem 的 Open 对话框默认显示 *.bin如果存储路径下混着多个测线的数据批量转换时容易漏选。更隐蔽的情况是 bin 文件里有部分测点没有正确写入延时道采集软件导出时直接跳过但软件不报错。解决回到采集软件重新确认导出范围选中整条测线确认每一个测点都出现在文件列表里再 Export。如果测点数量还是对不上检查是否有多组延时文件拆成多个 USF需要合并导入。5.2 剖面图的地形是平的现象生成剖面后所有测点高程完全一致剖面图上地形线是一条直线。原因USF 文件里没有高程列或者高程列全部为零。Terratem 在实测时需要接入高程数据才会写进 bin 文件野外接测距仪或 RTK 时如果没接好高程字段就是空的导出 USF 时自然也没有高程信息。解决单独导入高程数据。准备一个包含测点号和高程两列的 csv 文件在 IX1Dv3 的工区管理里执行导入高程操作导入后检查相邻测点高差与实际地形是否吻合。我当时处理某矿区数据时就是漏了这一步反演跑完才发现整条剖面地形是平的地形的低阻覆盖层厚度变化全都没反映出来只能重做一遍。5.3 反演到第五次迭代时拟合差反而上升现象迭代刚开始拟合差下降到第五次左右不降反升继续迭代也越来越差。原因初始模型电阻率初值离真实值太远正演响应在早期道严重溢出。常见于把浅层初始电阻率设得过高而实际浅层是低阻覆盖层正演出来的早期道电压比实测大几个量级梯度计算被带偏。解决回到初始模型把前两层电阻率改成早期道平均视电阻率重新启动反演。如果还是发散再检查数据编辑里是否保留了关断电流畸变的早期道这类坏道会让反演朝错误方向修正。5.4 初始模型层数一加多反演结果就乱现象两层模型跑得好加到五层后浅层出现极高或极低电阻率的怪异薄层拟合差虽然低了但模型完全不可信。原因层数增加导致自由度变多相邻层的电阻率与厚度互相补偿出现等值性范围内的虚假分层。这是 TEM 一维反演的固有问题不是软件 bug。解决减少层数或对可疑薄层的电阻率加约束范围。运行等价性分析如果目标层的厚度和电阻率可以在很大范围内变化而不改变拟合差说明该层是虚假自由度造成的要么合并到相邻层要么在报告里标注为不可靠层。5.5 成果只保存了图件找不到模型文件现象反演结果导出了 JPG 图隔几天想重新编辑模型发现项目里只剩图片反演模型文件不在。原因只执行了绘图导出没有单独保存反演模型。IX1Dv3 的“保存反演成果”和“导出图片”是两个动作导出了图片不等于模型数据被持久化。解决每次反演确认合理后立即用软件内保存功能把反演模型写入工区数据库并另外备份一份模型文件。从那以后我每个测点反演完都强制先保存模型再出图顺序反了就很可能白跑。5.6 同一测点重复观测两次反演结果差异大现象同一点位两次采集曲线形态相似但反演出来的深部电阻率差出一个数量级。原因深层信息在晚期道信噪比低反演对深部参数的分辨能力天然不足同时深部电阻率与浅层厚度之间存在等价性问题即“薄层低阻加厚层高阻”组合和“厚层稍高阻”组合可以产生几乎一样的响应。两次观测的微小噪声差异就会被反演放大到深部模型上。解决对深部参数加宽约束或固定范围用平滑模型反演替代层状模型反演先看整体趋势再分层。再把等价性分析结果附在报告里明确标注哪些层位是唯一解、哪些层位是范围解。我一般会在同一测区选 10% 的测点做重复观测专门用来评估深部反演结果的稳定性重复观测都稳的层位才敢写进解释结论。6. 进阶用多地层平顺模拟与等价性分析复核反演成果6.1 怎么用等价性分析判断结果可靠性反演收敛不等于结果唯一。TEM 一维反演最常见的坑是不同层序模型能产生几乎相同的响应——浅层低阻薄层加高阻基底和单一中电阻率厚层衰减曲线可能高度重合。这是电磁测深的固有局限硬件再好的仪器也躲不开。IX1Dv3 的人机交互解释功能里反向分层模型等价性分析就是专门对应这个问题的。实际操作时先跑完一层状模型反演然后打开等价性分析在拟合差允许的范围内观察各层电阻率和厚度的波动区间。波动区间小的层是唯一解可以放心写进报告波动区间大的层是范围解只能给出上下限不能给出精确值。判断标准可以归纳为层电阻率波动小于半个数量级、层厚度波动小于一倍的可信电阻率波动超过一个数量级或厚度可以在好几个区间游走的不可信。6.2 一个收尾习惯跨测点成图对比我现在的收尾动作有两步。第一步是在工区里生成完整的电阻率-深度剖面图按测线里程排列所有测点从横向上检查相邻测点反演结果的连续性。野外实测的衰减曲线再怎么有噪声相邻测点的地电模型也不该出现剧烈跳变如果某个测点的反演电阻率和两侧差出一个数量级先怀疑这个测点的数据质量或反演约束参数而不是急着解释成“局部异常”。第二步是每完成一条剖面就把层状模型反演和平滑模型反演的结果放在一起对比。层状模型给出分层框架平滑模型给出电性渐变趋势两者一致的部分是可靠成果不一致的部分大概率落在等价性模糊区需要在报告里如实说明。这套流程看着笨却是实打实的后悔药。从那以后我每次拿到一批 TEM 测深数据都强制自己走一遍“USF 验证、工区单位检查、按道编辑、初始模型从两层起步、收敛后做等价性分析、跨测点成图对比”的完整链路跑完的结果基本都能经得起复查。希望这几条实测经验帮到你少走点弯路。本文还有配套的精品资源点击获取
RELATED READING

延伸阅读

更多一线实战笔记与深度复盘,助您持续精进