磁偶源25维瞬变电磁场全空间FDTD数值模拟-物探与化探.PDFVIP

　第３２卷第３期 物　探　与　化　探 Ｖｏｌ．３２，Ｎｏ．３　 　２００８年６月 ＧＥＯＰＨＹＳＩＣＡＬ＆ＧＥＯＣＨＥＭＩＣＡＬＥＸＰＬＯＲＡＴＩＯＮ Ｊｕｎ．，２００８　 磁偶源２．５维瞬变电磁场全空间ＦＤＴＤ数值模拟 杨海燕，岳建华 （中国矿业大学 资源与地球科学学院，江苏徐州　２２１１１６） 摘 要：为解决模拟瞬变电磁场在３维介质中传播出现的计算量大、对计算机要求高等问题，使用了３维磁偶极源 来模拟２维均匀介质和层状介质中方形和薄板低阻体的全空间响应特征。模拟中使用了时域有限差分法，采用非 均匀网格和时间步长，分析了瞬变电磁场在均匀介质和层状介质中的传播规律，以及巷道和交界面对场的影响特 征，为解释全空间均匀介质和层状介质ＴＥＭ异常提供参考。 关键词：磁偶源；瞬变电磁场；时域有限差分法；模拟计算 中图分类号：Ｐ６３１　　　文献标识码：Ａ　　　文章编号：１０００－８９１８（２００８）０３－０３２６－０５ 　　瞬变电磁法的三维正演模拟一直是人们关注的 Ｅ（ｒ，ｔ）为电场强度，Ｅ（ｒ，ｔ）为磁场强度，、 式中， μσ 焦点，但由于计算机内存、ＣＰＵ等条件的限制，模拟 （ｒ）分别为磁导率和电导率。 难度非常大。多次的模拟经验表明，当计算的数据 采用非均匀网格将计算区域剖分成若干个长方 文件的大小超过内存后，计算速度变得非常慢，而 形单元，从而将求解无限大区域转化为求解各小长 ２．５维（二维模型三维源）模拟则可以进行。 方形单元。源附近网格较密，网格步长按等比级数 笔者从麦克斯韦方程组出发，使用时间域有限 递增。为保证ＦＤＴＤ计算稳定性，时间离散的步长 差分法（ＦＤＴＤ），推导出无源瞬变电磁场的二维有 （ｔ）与空间离散步长（ｘ、ｙ、ｚ）间应满足一定的 Δ Δ Δ Δ 限差分格式，并采用均匀全空间中阶跃波脉冲信号 关系。分析表明，时间步长可选为电磁波传播一个 ［５］ 下激发的磁偶极源（三维）作为初始源。为了模拟 空间步长所需时间的一半 。图１为在直角坐标系 较大区域，使用了非均匀网格对场域进行剖分。另 下场域的网格剖分模型，式（１）两边对图１中的矩 一方面，考虑到解的收敛性，初始时间步长最大只能 形ＡＢＣＤ积分得到 设置为 －８ １０ ｓ，而要得到毫秒级的模拟结果，迭代次 Ｈ ２  Ｈｄｘｄｚ＝ μσ ｄｘｄｚ。 （２） ５   ｔ 数将超过１０次，因此在模拟中也使用了非均匀时 ＡＢＣＤ ＡＢＣＤ  间步长，从而使工作量相对减小。文中仅模拟了均 匀介质和层状介质中巷道底板岩层内部方形和直立 板状低阻异常体的响应特征，绘制出了不同时刻磁 场强度的垂直分量Ｈ及其对时间的偏导 Ｈ／ｔ的   ｚ ｚ 等值线。 １　差分方程与源 １．１　差分方程的推