航空瞬变电磁全时域全空域快速成像.pdf

第 卷第期 年月物探化探计算技术 收稿日期：基金项目：陕西省教育厅科学研究计划项目（）第一作者：郑建波（），女，硕士，讲师，主要从事电磁勘探方法研究，：。文章编号：（）航空瞬变电磁全时域全空域快速成像郑建波，李美艳（西安外事学院工学院，西安 ）摘要：航空瞬变电磁法以其快速高效的优势已获得广泛应用，然而航空瞬变电磁采样密集，数据量巨大。为了实现航空瞬变电磁观测数据的快速解释，开展航空瞬变电磁快速成像算法研究。利用反函数定理建立成像迭代格式，并在成像过程中考虑发射源高度和观测时间的影响，进而实现航空瞬变电磁全时域全空域快速成像。这里首先给出了全时域全空域成像理论和实现方法，然后对航空瞬变电磁装置进行了介绍，最后建立了典型采空区模型，利用开发的航空瞬变电磁全时域全空域视电阻率成像方法进行数据成像，成像结果与真实模型基本一致，验证了方法的有效性。关键词：航空瞬变电磁；全时域；全空域；成像中图分类号：文献标志码：犇 犗 犐：引言航空瞬变电磁法将电磁探测装置搭载在飞机上，通过发射线圈进行大功率发射，激发地下感应涡流，通过观测断电间歇感应涡流产生的二次磁场的空间和时间分布特征，对地下目标进行快速电磁探测。由于该方法具有快速、高效、受地形影响小等特点，已被广泛应用于矿产勘查、环境地质调查等诸多领域。由于航空瞬变电磁法在探测过程中进行连续采样，获得海量观测数据，因此具有较高的分辨率。然而，海量数据量也给航空瞬变电磁数据解释带来了巨大困难。为了能够实现航空瞬变电磁数据的快速解释，成像方法成为了首选。航空电磁在数据采集过程中，飞行高度会随着地形起伏发生变化，观测时间跨度大，这给传统的成像方法提出了挑战。为此，笔者提出了一种航空瞬变电磁全时域全空域快速成像方法，进而建立了一套不受观测装置、观测时间和飞行高度限制的通用成像系统。全时域全空域视电阻率成像思想最早由张莹莹等 提出。该方法起初被应用于电性源地空瞬变电磁数据电阻率成像中，通过对偏移距、发射源位置、电流方向等参数的影响进行分析，建立视电阻率迭代格式实现了全时域、全空域电阻率成像。赵越等 将该方法应用于大定源回线地空瞬变电磁电阻率成像中；武军杰等 利用全时域、全空域视电阻率成像方法实现了电性源地井瞬变电磁数据解释；曾友强等 进一步将其应用于金矿水害探测当中。笔者将全时域、全空域电阻率成像思想引入到航空瞬变电磁数据成像领域，建立一套通用的航空瞬变电磁快速成像方法。电阻率成像理论 航空瞬变电磁探测原理分析航空瞬变电磁系统利用飞机作为移动搭载平台，利用多匝线圈作为发射源，向地下发送一次电磁激发信号，然后在一次场断电间歇测量由地下介质中感应涡流产生的二次磁场。当地下存在电性异常体时，观测的二次磁场将发生畸变，因此通过分析二次磁场的空间和时间分布规律就可以实现对地下电性结构的探测。以地下含水采空区探测为例，由于采空区充水后，相对周围岩石介质表现为低阻特征，即采空区与周围大地介质间存在明显的电性差异，这种巨大的电性差异，将导致观测的二次磁场在空间上表现为正异常，在时间上表现为慢衰减特征，这就是利用航空电磁法探测地下结构的物理机理。利用视电阻率成像方法，便可以将二次磁场的空间和时间分布特征转化为地下结构的电性分布特征，进而揭示地下目标体的三维空间分布。时间域航空电磁数值模拟方法均匀大地表面，一维频率域航空电磁响应由下式给出。犎狕（，犺，狕，）犿 狕犚犚犜（狕）（）式中：为大地电阻率；犺为飞行高度；狕为接收点坐标；为角频率；犿为发射磁矩；狕犺狕；犚狉狕槡；狉为接收线圈半径；犜（狕）是阶贝塞尔积分。利 用 汉 克 尔 变 换 技 术 可 以 完 成 式（）中犜（狕）的计算。这里采用阶跃电流发射波形，其满足：犐（狋）犐狋狋（）即在时刻之前供稳定电流，在时刻突然断电。我们采用正 弦变化 方法将 频 率 域 航 空 电 磁 响 应犎狕（，犺，狕，）转换到间域，获得二次场犅狕（，犺，狕，狋）狋槡 犎狕（，犺，狕，）犻 槡犑（狋）（）其中：表示取虚部；犑为 阶贝塞尔函数，再次利用汉克尔变换即可完成时间域航空电磁阶跃响应计算。全时域全空域视电阻率成像方法当航空电磁系统数据采集结束后，可以获得准确的装置参数和观测信息，包括飞行高度犺，接收点坐标狕和观测时间狋，此时观测数据犅狕（，犺，狕，狋）成为了一个只与地下电阻率有关的物理量，现将其改写为犅（）。通过大量数值实验可以发现犅（）响应是关于电阻率的单调函数。根据反函数定理可知，必然存在一个电阻率值唯一对应着一个犅（）值。为此，我们给出一个电阻率初值，然后在点对犅（）进行级数展开犅（）犅（）犅（）（）！犅（）（）狀！犅狀（）（）狀（）为了建立迭代格式，我们忽略高阶项，只保留其线性主部，建立视电阻率迭代格式为式（）。犻 犻 犻（）其中犻是电阻率修正量，由式（）给出。犅（，犮，狋）犅（犻）犅（犻）（）利用式（）进行多次迭代后，当满足式（）条件时，迭代结束。犅 犅（犻）犅 （）其中：犅 是实际观测数据；犅（犻）是电阻率为犻的半空间模型航空电磁理论响应，犅（犻）由式（）计算；是迭代终止误差。由于在式（）中，同时考虑了飞行高度和观测时间的影响，因此基于上述迭代格式的视电阻率成像结果将不再受到飞行高度和观测时间的影响，实现全时域全空域视电阻率成像，获得更加精确的地下电性分布信息。航空瞬变电磁装置航空瞬变电磁探测装置如图所示，主要包括飞行器装置（图中的直升机）、吊舱、供电系统、发射系统（包括发射机和发射线圈）、接收系统（包括接收机和接收线圈）、辅助系统（提供飞行时必要的辅助信息，包括高度、经纬度）等，飞行器下方挂着吊舱，吊舱 承 载 着 发 射 和 接 收 线 圈。供 电 系 统、接收机及辅助系统这三部分放在机舱内，发射机可置于机舱 或 者 吊 舱内。对 于 航 空 瞬 变 电 磁 系 统来说，发射系统至关重要，笔者重点介绍发射系统关键技术。发射系统 电源取自直升机的直流发电机，而发射系统最终的发射波形是如图所示的双极性梯形波，即要实现从直流到交流的转换，所以主电路采用如图所示的全桥逆变器。在图中，为全控型电力电子器件，实际工作时，和 为一组、和 为一组，同一组开关的导通和关断时间是一致的，通过交替控制两组开关的导通和关断时期郑建波，等：航空瞬变电磁全时域全空域快速成像图航空瞬变电磁装置 图双极性梯形波电流示意图 图航空瞬变电磁发射桥路拓扑图 刻，实现输出波形为双极性梯形波。全控器件的控制脉冲基于（脉冲宽度调制）原理实现，具体过程如下：根据输出双极性梯形波的电流幅值及每一部分的持续时间，依据电阻和电感的伏安特性方程计算出每一段所对应的输出电压，然后以此输出电压波形作为调制波，以三角波作为载波，准确计算出每个开关器件的导通和关断时刻。图不同角度视图的采空区模型 理论模型验证笔者以采空区探测为例，进行理论模型验证。在煤矿采煤结束后，通常会撤去支护，留下一个空腔，此时地下水会沿着地下裂隙进入。此外，随着时间增加，由于重力作用，地下空腔还会出现垮塌，进而形成采塌陷空区。在没有采空的地方，地下岩石比较致密，导电性较差，表现为高阻特性。而采空区通常充水，导电性较好，表现为低阻特性。为了模拟采空区特征，设计了如图所示的采空区模型。一个电阻率为 的低阻椭球体采空区埋藏在 的高阻围岩中。地下采空区的长轴半径为 ，短轴半径为 的椭球体，中心埋深为 。采用航空电磁装置进行探测，其发射线圈为一边长为 的正 边形，线圈匝数为匝，发射如图所示的双极性梯形波，峰值发射电流为 。单个梯形波供电时长犜，其中上升沿时间为 ，平顶沿时间为 ，下降沿时间为 。两个双极性梯形波之间的断电时长为 ，因此单个波形总时长犜 。在每个梯形波断电后观测地下产生的纯二次场。进行面积性测量，点距和线距均为 ，整个测区共有 个测点，观测断电后 对数等间隔的 个时间道的数据。利用齐彦福 物探化探计算技术 卷图全域视电阻率切片图 图全域视电阻率三维透视图 在 年开发的有限元正演程序模拟理论模型航空电磁观测数据。在获得正演数据后，采用笔者开发的航空电磁全域视电阻率成像方法进行数据成像。对于所有测点、所有时间道的初始模型均采用 的均匀半空间模型。该数值算例在 机上进行，主频 ，计算机内存 。图和图展示了视电阻率成像结果，可以看出，全域视电阻率成像技术能够清晰反映出地下电性分布情况。从图可以看出，在椭球体中心位置的视电阻率最低。随着距离增加，视电阻率逐渐升高，逐渐趋于背景围岩电阻率。对比图和图，我们可以看到笔者开发的任意发射波形全域视电阻率可以很好地反应出地下电性信息，准确地恢复了地下良导椭球体的形状、位置和尺寸，这为后续地质解释提供了基础，本次成像共耗时 。结论笔者针对航空电磁采样密集，数据量巨大的特点，成功建立了一套航空瞬变电磁全时域全空域快速成像算法。利用时间域航空电磁观测响应与地下介质电阻率间的单调关系，利用反函数定理建立了视电阻率迭代格式，从而实现了考虑装置形式和观测时间的航空瞬变电磁全时域全空域快速成像。通过将成像方法应用于理论采空区模型，结果显示成像结果可以很好地反映地下采空区的真实分布，从而验证了算法的可靠性。参考文献：宿晓明，李世华黑龙江小兴安岭丽林清水地区利用航空瞬变电磁法找矿前景分析世界有色金属，（）：，（）：（）彭莉红，江民忠，程莎莎，等航空瞬变电磁法在新疆阿期郑建波，等：航空瞬变电磁全时域全空域快速成像奇山地区铅锌多金属找矿中的应用研究地质与勘探，（）：，（）：（）何怡原，梁盛军，连晟基于航空瞬变电磁法的宝清地区含水体分布特征研究地球物理学进展，（）：，（）：（）彭莉红，张策，骆燕，等航空瞬变电磁法在塔源二支线铅锌铜矿床的应用研究物探化探计算技术，（）：，（）：（）赵思为，尹小康，张振雄，等半航空瞬变电磁法在某滑坡探测中的应用研究物探化探计算技术，（）：，（）：（），（）：，（）：，（）：，（）：张莹莹，李貅，姚伟华，等多辐射场源地空瞬变电磁法多分量全域视电阻率定义地球物理学报，（）：，（）：（）赵越，李貅，王鹏大回线源瞬变电磁全域视电阻率定义地球物理学进展，（）：，（）：（）武军杰，李貅，智庆全等电性源地井瞬变电磁全域视电阻率 定 义 地 球 物 理 学 报，（）：，（）：（）曾友强，谢国胜，智庆全，等小回线源 法全域视电阻率成像在某金矿水害中的应用地球物理学进展，（）：，（）：（）殷长春，黄威，贲放时间域航空电磁系统瞬变全时响应正演模拟地球物理学报，（）：，（）：（）姜健时间域航空电磁法双波形发射电路研制长春：吉林大学，：，（）李琳琳半航空瞬变电磁发射机关键技术研究成物探化探计算技术 卷都：成都理工大学，：（）齐彦福，殷长春，刘云鹤，等基于瞬时电流脉冲的三维时间域航空电磁全波形正演模拟地球物理学报，（）：，（）：（）犃 犲 狉 犻 犪 犾 狋 狉 犪 狀 狊 犻 犲 狀 狋 犲 犾 犲 犮 狋 狉 狅 犿 犪 犵 狀 犲 狋 犻 犮犪 犾 犾 狋 犻 犿 犲 犱 狅 犿 犪 犻 狀犪 狀 犱犪 犾 犾 狊 狆 犪 狋 犻 犪 犾 犳 犪 狊 狋 犻 犿 犪 犵 犻 狀 犵 ，（，）犃 犫 狊 狋 狉 犪 犮 狋：，犓 犲 狔 狑 狅 狉 犱 狊：；期郑建波，等：航空瞬变电磁全时域全空域快速成像