三维电磁响应积分方程法数值模拟.pdfVIP

第 30 卷第 5 期 中南工业大学学报 Vol. 30 　No. 5 　　　　　　　　　　　　 　　　　　　　　　　　　　　 1 9 9 9 年 1 0 月 J . CENT. SOUTH UNIV. TECHNOL. Oct. 　1 9 9 9 三维电磁响应积分方程法数值模拟 鲍光淑 　张碧星 　敬荣中 　毛先进 ( 中南工业大学资源环境与建筑工程学院 ,长沙 ,410083) 摘 　要 　从电磁基本理论出发 ,深入研究了均匀导电半空间频域三维电磁散射问题 , 由张量格林函数求解 Fredholm 方程 ,将异常体划分为一系列立方单元 ,并将该积分方程转化为矩阵方程 ,从而计算出空间任意点的电场和磁场 ; 此外 ,还给出了计算实例 ,并与国外结果进行分析对比 ,证明作者运用的计算方法结果正确 ,算法快速. 关键词 　三维电磁响应 ; 张量格林函数 ; 积分方程法 ; 数值模拟 分类号 　P318 ; P319 　　在实际应用中 ,无论是频率域电磁法还是时间 与国外存在较大差距. 为了提高电磁法解释水平 ,推 域电磁法 ,都获得了比较好的地质效果 ,引起国内外 动我国电磁法的应用发展 ,必须加强这方面的应用 学者的普遍重视 ,得到愈来愈广泛的应用. 目前, 国 基础理论研究. 为此 ,作者进行了均匀导电半空间频 内外在电磁法的仪器设备方面不断推出新产品,发 率域三维电磁响应的数值模拟. 展很快 ,然而 , 电磁法的解释技术相对落后,制约了 1 　基本理论 电磁法的发展. 野外地质情况复杂 ,在用电磁法解决 各类地质问题时 ,迫切需要二维和三维正、反演软件 　　所有电磁现象均可用麦克斯韦方程描述. 在一 作为资料解释的重要工具. 而高维电磁响应的数值 ( 5 ) 定频率范围内 10 Hz ,可忽略位移电流. 因此 ,麦 模拟和反演成像是一项非常复杂的工作 ,也是当前 克斯韦尔方程变成一组扩散方程 : 勘探地球物理领域前沿的研究课题. 尤其是三维体 2 H E - uσ = 0 t 电磁响应的数值模拟,更是迫切需要解决的问题 ,它 2 H 是反演成像的基础 , 同时对深入理解所测量到的三 H - μσ = 0 t 维体电磁场响应的物理性质、场的分布规律和特点 , 2 μσω E - i E = 0 是非常重要的 ,是解释实际资料的理论依据.