基于QWE法的汉克尔变换在电磁法数值计算中的应用  

Hankel transform based on QWE in the electromagnetic numerical computation

作　　者：胡俊华[1] 昌彦君[1] 左乾坤[1] 孙卫斌 武光海[3] 

机构地区：[1]中国地质大学(武汉)地球物理与空间信息学院,地球内部多尺度成像湖北省重点实验室,湖北武汉430074 [2]东方地球物理公司,河北涿州072750 [3]国家海洋局第二研究所,浙江杭州310000

出　　处：《石油地球物理勘探》2015年第2期370-375,7-8,共6页Oil Geophysical Prospecting

基　　金：“大型油气田及煤层气开发科技重大专项(2011ZX05019-007)”;“国际海域资源调查与开发‘十二五’课题(DY125-11-R)”联合资助

摘　　要：在电磁模拟的数值计算中,需要有高性能的数值滤波方法计算汉克尔积分,以获得良好的计算精度和运算速度。将直接数值积分方法即基于Shanks变换的正交插值求积方法(QWE)应用于两个实际问题:重叠回线装置的一维瞬变电磁响应和电偶源电磁张量格林函数的计算。与快速汉克尔变换(FHT)计算结果相比较,对于重叠回线装置,应用QWE法可获得更准确的垂直磁场高频段虚部响应和早中期瞬变响应,效率接近FHT法。对于电磁张量格林函数的计算,选用5节点和100个求积段数仅有0.1%的误差,在相同精度下QWE法比FHT法节省一半计算时间,显示QWE法在电磁模拟的数值计算中具有潜在的优势和应用价值。A high-performance numerical filtering method is usually needed to calculate Hankel integral in electromagnetic numerical simulations in order to obtain good calculation accuracy and computation speed.We propose in this paper a direct numerical integration method,quadrature-with-extrapolation(QWE)based on Shanks transform,and use it to compute one dimensional transient electromagnetic response of coincident-loop configuration and electromagnetic tensor Green's function of an electric dipole source.For the former,the QWE method can obtain more accurate imaginary components of vertical magnetic field at high frequencies and lead to more precise transient response in early and mediate times correspondingly in comparison with the fast Hankel transform(FHT)method,while its efficiency does not differ with FHT.For the latter,the QWE method with 5nodes and 100 intervals has only 0.1% error,and saves half the time of the FHT computation under the same precision.These results show the QWE algorithm has the potential advantage and application value in the numerical computation of electromagnetic simulation.

关 键 词：快速汉克尔变换 外推式正交插值求积 Shanks变换 瞬变电磁法 张量格林函数 

分 类 号：P631[天文地球—地质矿产勘探]