FDTD／MOMTD混合方法在探地雷达散射问题中的应用 - Journal of.PDF

2007年2月 重庆大学学报(自然科学版) Feb．2oo7 第3O卷第2期 Journal of Chongqing University(Natural Science Edition) Vo1．3O No．2 文章编号：1000-582X(2007)02-0050-03 FDTD／MOMTD混合方法在探地雷达散射问题中的应用 周 超，俞集辉 (重庆大学电气工程学院，重庆 400030) 摘 要：为了分析和研究探地雷达的地下散射问题，构建了时域有限差分法和时域矩量法相结合的 混合方法，并采用该混合方法进行仿真，当散射良导体距离地面高度不同时，计算得出探地雷达接收天 线中心馈电点的电流波形，通过定性和定量分析证明该混合方法所得的结果的正确性． 关键词：时域有限差分法；时域矩量法；混合方法 中图分类号：TN959 文献标识码：A 探地雷达(Ground Penetrating Radar，简称GPR)方 × ： + ． (2) 法是一种用于确定地下介质分布的广谱(1MHz一 1GHz)电磁技术．利用探地雷达的一个天线发射电磁 FDTD在空间和时间上采取交替抽样的离散方 波，另一个天线接受来自地下介质的反射波．考虑到电 式[3】，每一个 (或日)场分量周围有4个 (或日)场 磁波在介质中传播时，其路径、电磁场量强度及其波形 分量环绕，形成了将麦克斯韦旋度方程转化为含时间 将随所通过介质的电性质和几何形态而变化．于是，可 变量的一组差分方程，并在时间轴上逐步推进地求解 根据接收到电磁波的旅行时间(亦称双程走时)、幅度 空间分布的电磁场． 与波形等信息，来分析、推测地下介质的分布和结构特 用MOMTD处理电场积分方程 征．由此可反映出探地雷达的广阔应用前景． (r,t)：一 一 (叫)， (3) 探地雷达的应用，其关键点是复杂的地下散射电 式中： 磁场的数值仿真，它已受到人们广泛的关注．考虑到探 地雷达发射和吸收的都是瞬态脉冲波，若采用时域方 A(r,t)= 等dV ， (4) 法分析探地雷达应当比采用频域方法会显得方便有 效．时域有限差分法 (FDTD)对于复杂情况的电磁 (r,t)= I． 等dV ． (5) 散射问题有着广泛的应用，而时域矩量法-3 (MO- 在一维情况下，对式(3)一(5)进行离散化 ， MTD)对于处理脉冲天线辐射问题有着良好的表现． 从而逐步求解出各时间步上各个离散点的时变电流密 文献[5]的作者充分利用它们的长处，把这2种方法 度．，和时变电荷密度P，再据此求解空间任意点的 结合起来形成了其混合方法．文献[6]的作者应用其 和日．其中 混合方法计算基站天线的吸收电磁能量的比值，得到 日(r，‘)= ×A(r，‘)． (6) 了很好的效果．该论文将其混合方法用以求解探地雷 达的地下散射问题． 若使用这2种方法考虑在时间上同步，在计算空间上 相吻合H-2]．先利用MOMTD式(3)一(5)，以及由 1 FDTD／MOMTD混合方法的构成思路 FDTD方法解得的散射体在天线上的散射电场，根据 对于麦克斯韦旋度方程