瞬变电磁场直接时域数值分析.pdfVIP

第 45 卷 第 2 期 地 球 物 理 学 报 Vol. 45 , No. 2 　　　　　 　　　　　　 　2002 年 3 月 CHINESE JOURNAL OF GEOPHYSICS 　 Mar. , 20 0 2 ( ) [文章编号] 　0001 - 5733 2002 02 - 0275 - 10 [ 中图分类号] 　P631 瞬变电磁场的直接时域数值分析 1) 2) 1) 闫　述 　陈明生 　傅君眉 1) 西安交通大学电子与信息工程学院微波与光通信研究所 ,西安　710049 2) 煤炭科学研究总院西安分院 ,西安　710054 [摘 　要] 　为了深入了解瞬变电磁场的勘探原理 ,直接在时间域对负阶跃脉冲激发的二维瞬 态场进行了数值分析. 采用的方法是从反映电磁场基本规律的麦克斯韦方程组出发 ,导出时 域电场的齐次扩散方程 ,对所研究的空间区域作差分离散 ,源作为初始条件加入 ,利用准静态 近似处理空中边界 ,然后进行时间的逐步递推 , 由此展现瞬变电磁场在地下扩散随时间发展 的全过程. 通过模拟计算不同时刻瞬态电场在地下的分布形态及地面上感生电动势相应的 变化 ,揭示了低阻异常体对感应涡流的聚集作用 ,低阻覆盖层对瞬变场扩散的减速作用 ,及瞬 变场的延时效应. 因此 ,瞬变电磁法对低阻体是敏感的 ,有上覆低阻层时探测同样的深度需 要较长的时间 ,而延时效应瞬变场的晚期时段可反映埋藏较浅的异常体. [ 关键词] 　瞬变电磁场 ,有限差分 ,扩散 ,“烟圈”概念,延时. 1 　引　　言 对于瞬变电磁场响应特征的分析 ,大多是将相应情况下的频率域结果转换到时间域 , 而且在地面上接收装置中的响应. 如水平层状大地上的瞬变响应[1 —3 ] 、导电半空间中或 层状大地中异常体响应的分析[4 —6 ] ,大都采取了将频率域结果经傅里叶变换转换到时间 域的方式. 虽然也有直接在时间域用积分方程的形式进行分析的[7 ,8 ] ,仅根据地面上的响 应对瞬变电磁场的探测原理作出物理解释 ,分析瞬变场方法的探测深度、探测精度、分辨 率等是不够的 ,而且积分方程法不具有一般性 ,每次计算不得不对具体的几何边界和介质 特性进行再推导. 微分形式的方法主要是有限差分法和有限元法[9 ,10 ] ,在时间域 ,有限差 分的方法更为方便[11 —13 ] . Oristaglio 等[14 ] 和 Wang[15 ] 等提出了地球物理勘探中瞬变电磁场 的时域有限差分算法 ,即在每一时间步计算网格空间各点的场量 ,随着时间步的推进 ,就 能直接模拟电磁波的传播及其与地质体的相互作用过程 ,使电磁场的时域特性被直接反 映出来. 这一特点使它能直接给出非常丰富的电磁场问题的时域信息 ,给复杂的物理过 程描绘出清晰的物理图像. 这样 ,直接在时间域对瞬变电磁场的信号进行分析 ,可对瞬变 [ 收稿日期] 　2000 - 11 - 16 收到 ,2001 - 07 - 31 收到修定稿. ( ) [ 基金项目] 　国家自然科学基金项目 . [ 作者简介] 　闫　述 ,女 ,1953 年生 ,1978 年毕业于南开大学物理系 ,高级工程师. 现为西安交通大学在读博士生 ,主 要从事电磁法勘探的研究工作. Email :chenms @pub. xaonline. com ( )