High-order FDTD method for multi-offset electromagnetic numerical imaging
WANG Ke1,2
1. College of Electrical and Information Engineering, Guangxi University of Science and Technology, Liuzhou 545006, China; 2. Hunan Railway Professional Technology College, Zhuzhou 412001, China
Abstract：Numerical simulation is an important part of data processing and inversion interpretation. In order to improve the efficiency and the accuracy of electromagnetic wave imaging technology in engineering survey, it is important to understand the propagation characteristics of electromagnetic waves by forward modeling. In this study, the finite-difference time-domain FDTD (2, 4) was used to simulate the numerical imaging of electromagnetic waves in three kinds of geoelectric models, i.e., karst caves, faults and boulders. The multi-offset mode was used to acquire the electromagnetic wave gathers. Based on the analysis of the electromagnetic wave wavefield and the waveform profile, the authors summarized the propagation law of electromagnetic waves in different geological bodies and provided reliable basis for practical engineering investigation. The experimental results show that FDTD (2, 4) can simulate the numerical imaging of electromagnetic waves in complex media with high accuracy, and also prove that the multi-offset acquisition mode has high efficiency and flexibility.
[1] 赵连锋,朱介寿,严忠琼,等.岩溶勘察层析成像试验:成像效果分析[J].物探与化探计算技术,2002,24(4):332-336. [2] 陈昌彦,沈小克,苏兆锋,等.电磁波层析成像技术在复杂地质边坡工程勘察中的应用研究[J].地球物理学进展,2012,27(2):0796-0803. [3] 冯锐.电磁波层析成象[J].地震学报,1997,19(5):524-534. [4] 张曼舳,师学明.电磁波层析成像技术进展[J].工程地球物理学报,2009,6(4):418-425. [5] 刘昌铨,孙振国,李长法.层析成像技术在核反应堆地基勘察中的应用[J].物探与化探,2000,24(1):69-75. [6] YEE K S. Numerical solution of initial boundary value problems involving Maxwell��s equations in isotropic media [J]. IEEE Transactionon Antennasand Propagation, 1966, 14(4):302. [7] 林树海,赵立英.井问电磁波层析成像中的高精度时间域正演计算[J].中国地质大学学报:地球科学,2007,32(4):469 473. [8] Fang J. Time domain infinite difference computation for Maxwell’s equations [D]. Berkeley:Univ California, 1989. [9] Lan K, Liu Y W, Lin W G. A higher order (2,4) scheme for reducing dispersion in FDTD algorithm [J]. IEEE Trans Electromagnetic Compatibility, 1999, 41(2):160-165. [10] 吴丰收,曾昭发,黄玲,等.探地雷达信号的高阶时间域有限差分模拟[J].物探与化探计算技术,2009,32(4):308-313. [11] 孔繁敏,李康,郭毅峰.高阶FDTD法分析电-大尺寸光波导器件[J].光电子·激光,2004,15(6):672-6789. [12] 王健.高阶FDTD方法及其在散射问题的应用[D].安徽:安徽大学,2007. [13] 曾昭发,刘四新,王者江,等.探地雷达方法原理及应用[M].北京:科学出版社,2006. [14] 冯德山,戴前伟,何继善,等.探地雷达GPR正演模拟的时域有限差分实现(英文)[J]. 地球物理学进展,2006,02:630-636. [15] 葛德彪,阎玉波.电磁场时域有限差分法[M].西安:西安电子科技大学出版社,2002.
