电磁波场数值计算中时间域有限差分法与时间域伪谱法的对比研究

现代地质 ›› 2012, Vol. 26 ›› Issue (6): 1193-1198.

电磁波场数值计算中时间域有限差分法与时间域伪谱法的对比研究

林树海，王伟利

中国地质大学地球物理与信息技术学院, 北京100083

收稿日期:2012-06-03 修回日期:2012-08-11 出版日期:2012-10-24 发布日期:2012-10-31
作者简介:林树海，男，副教授，1964年出生，地球物理学专业，主要从事电磁法勘探与地球物理测井方面的研究工作。Email: lin_shuhai@163.com。

Study on the Differences Between Finite Difference Time Domain Method and Pseudo-spectral Time Domain Algorithm in Numerical Calculation of Electromagnetic Wave Fields

LIN Shu-Hai, WANG Wei-Li

School of Geophysics and Information Technology, China University of Geosciences, Beijing100083, China

Received:2012-06-03 Revised:2012-08-11 Online:2012-10-24 Published:2012-10-31

摘要/Abstract

摘要：

为提高电磁波场数值计算精度，对时间域有限差分法与时间域伪谱法进行了对比研究。时间域有限差分法是一种目前流行的电磁场时域数值计算方法，已被广泛应用于求解与时间有关的偏导数方程。对于大规模数值计算，时间域有限差分法需要较多的内存空间。时间域伪谱法基于时间域有限差分法，该方法使用快速傅里叶度换来计算麦克斯韦方程中的空间导数。由于傅里叶变换的准确性，时间域伪谱法使计算精度提高，数值计算时所需的格子数大大减少，这极大地节省了计算机内存空间，适合于大规模正反演问题的数值计算。

关键词: 数值计算, 时间域有限差分法, 时间域伪谱法, 快速傅里叶变换, 精度

Abstract:

In order to improve the resolution of numerical calculation of electromagnetic wave fields, the differences between finite difference time domain (FDTD) method and psedo-spactral time domain (PSTD) algorithm were studied. The FDTD method is a common numerical calculation method of electromagnetic wave fields, and it is widely applied to the solution of the partial derivative equations related to time.More memory is needed using FDTD method in large-scale numerical calculation. The PSTD algorithm is based on FDTD method and employs fast Fourier transform (FFT) to calculate the spatial derivatives in Maxwell's equations. Due to accuracy by FFT, the calculation precision is improved and grid numbers needed is decreased greatly by using PSTD algorithm. The PSTD algorithm has the advantage of saving computational memory space and is suitable for large-scale numerical calculation of forward modeling and inversion.

Key words: numerical calculation, finite difference time domain method, pseudospectral time domain algorithm, fast Fourier transform algorithm, precision

林树海, 王伟利. 电磁波场数值计算中时间域有限差分法与时间域伪谱法的对比研究[J]. 现代地质, 2012, 26(6): 1193-1198.

LIN Shu-Hai, WANG Wei-Li. Study on the Differences Between Finite Difference Time Domain Method and Pseudo-spectral Time Domain Algorithm in Numerical Calculation of Electromagnetic Wave Fields[J]. Geoscience, 2012, 26(6): 1193-1198.

[1]	蒙嘉琪, 王志猛, 贾三石, 付建飞, 张岩松. 基于地球物理技术圈定沉积变质型富铁矿体: 辽宁鞍山齐大山铁矿例析[J]. 现代地质, 2024, 38(01): 87-97.
[2]	闫迪, 付建飞, 贾三石, 程浩, 王志猛. 鞍山—本溪地区眼前山铁矿床深(边)部错断隐伏矿体磁法和瞬变电磁法联合定位及预测[J]. 现代地质, 2024, 38(01): 35-45.
[3]	陈格格, 高志前, 焦存礼, 胡宗全, 袁钰轩, 卫端, 翟昕箐, 畅哲. 塔里木盆地北部下寒武统肖尔布拉克组海进-海退转换背景下的微生物礁滩发育模式[J]. 现代地质, 2023, 37(05): 1169-1181.
[4]	胡子奇, 张德贤, 刘磊. 束斑直径和能量密度对锆石U-Pb定年准确度的影响研究[J]. 现代地质, 2023, 37(03): 722-732.
[5]	李兆令. 河北小山火山的地磁场分布特征及其地质意义[J]. 现代地质, 2021, 35(05): 1480-1486.
[6]	芮志锋, 林畅松, 郭佳, 杜家元, 丁琳, 李潇. 珠江口盆地惠州地区珠江组砂体上倾尖灭的地质-地球物理“逐级预测”方法[J]. 现代地质, 2019, 33(06): 1229-1240.
[7]	韩剑发，王清龙，陈军，林畅松，宋玉斌，徐俊博，李浩，薛学亚，王宇盟. 塔里木盆地西北缘中—下奥陶统碳酸盐岩层序结构和沉积微相分布[J]. 现代地质, 2015, 29(3): 599-608.
[8]	左璠璠，林畅松，高达，赵建华，夏世强，李虹. 塔中地区西北部良里塔格组沉积特征及其演化规律[J]. 现代地质, 2014, 28(5): 1008-1016.
[9]	张启升, 邓明, 杨鹍鹏, 陈凯, 王猛, 孙志成. SoPC技术在高精度地电数据采集中的应用[J]. 现代地质, 2012, 26(6): 1306-1311.
[10]	贾冀辉, 牛滨华, 孙春岩, 郭继亮, 孙逊, 刘思思. 水合物密度的数值计算及储层识别[J]. 现代地质, 2011, 25(6): 1180-1186.
[11]	吴怀春, 王成善, 张世红, 杨天水, 万晓樵. “地时”（Earthtime）研究计划：“深时”（Deep Time）记录的定年精度与时间分辨率[J]. 现代地质, 2011, 25(3): 419-428.
[12]	刘思思, 牛滨华, 孙春岩, 杨奎. 饱水白云岩临界点、骨架和流体弹性参数的数值计算[J]. 现代地质, 2011, 25(1): 129-136.
[13]	蔡振忠, 张丽娟, 林畅松, 雷燕平, 苏真真, 司宝玲. 库车坳陷却勒—羊塔井区库姆格列木群底砂岩段高精度层序对比及油气开发意义[J]. 现代地质, 2010, 24(4): 662-668.
[14]	刘景彦林畅松常艳艳高达张宇. 贺兰山—桌子山地区石炭—二叠系陆表海高精度层序地层模拟与聚煤分析[J]. 现代地质, 2009, 23(4): 625-630.
[15]	张英志林畅松马利民刘景彦张宇梁秀丽. 大庆油田葡一油层组高精度层序格架内储集砂体的沉积构成及其对油田开发的意义[J]. 现代地质, 2006, 20(2): 225-231.