煤田地质与勘探

作者：梅竹虚, 原源, 周峰, 周聪, 张林成,

关键词：射频大地电磁法, 位移电流, 有限元正演, 并行计算, Julia语言,

射频大地电磁法(Radio-magnetotelluric, RMT)是浅地表电磁勘探的重要手段之一。由于勘探频段为10~300 kHz，其电磁场传播受地下介质介电常数的影响较大。传统准静态条件下的电磁响应严重制约RMT正演模拟精度，并进一步影响反演成像分辨率。针对这一问题，提出了一种基于Julia并行加速的全电流RMT电磁响应数值模拟方法，利用Julia的分布式计算将各频点的计算发送到不同进程进行求解，从而达到加速计算的目的，同时在计算中考虑位移电流的影响，提升正演的模拟精度。通过计算几种典型高阻/高介电模型的RMT响应，分析并总结了位移电流对射频段电磁场视电阻率及相位响应的影响规律。数值模拟结果表明：当浅部存在高阻覆盖层时，基于准静态假设条件计算的RMT视电阻率和相位响应偏高，且频率越高、覆盖层电阻率越大，响应偏差越大；对于煤炭采空区模型，RMT法能有效反映异常体位置，但忽略位移电流会在采空区及其附近引起较大计算误差；起伏地形的算例表明地形会覆盖地下异常体的RMT数值响应，尤其是地形拐角处；2种不同规模的并行算例对比证明了并行算法的高效性，且随着求解问题规模增大，并行算法效率也随之增大。本文研究提高了RMT的正演速度与精度，为后续快速反演算法的实现奠定基础，数值算例表明RMT法在煤田采空区等实际勘探中具有较好的应用前景。

来源：2023年第3期

《煤田地质与勘探》期刊编辑部

查看煤田地质与勘探杂志2023年第3期