摘要
采用三角网格剖分的有限元法,研究了2.5维航空瞬变电磁法正演模拟问题。利用时频变换数值方法将时间域电磁场转换到拉氏域,再利用傅里叶变换将三维问题降维变为2.5维问题,然后由有限元法求解得到拉氏域二维电磁场,逆拉氏变换后得到时间域航空瞬变响应。为了回避正演模拟中总感应磁场在场源处的奇异性问题,采用异常场算法,场源响应通过在微分方程中施加背景电磁场实现。由于瞬变电磁信号具有较大的动态范围,而且需要经过两次正、逆拉氏变换和傅里叶变换,每个环节的计算精度和速度要严格控制在较高的水平上,否则积累误差会非常大。模型计算表明均匀大地和层状大地模型解析解与数值解吻合很好。这证明该算法是正确可行的,可作为研究二维复杂地质体的方法手段。
Based on the triangular mesh dissection of the finite element method, this paper presents the 2.5D airborne transient electromagnetic method forward modeling. Firstly the time domain electromagnetic partial equations are converted into laplace domain with the numerical timefrequency transform algorithm, then the 2.5D expressions are obtained from 3D case by adopting the fourier transform. The airborne data are acquired from the results of the digital inverse laplace transform of the 2D electric and magnetic components in laplacian domain, which is computed using finite element method. Instead of the conventional algorithm, the anomalous field method is employed throughout the finite element process, so that the source singularity of total field is avoided. the source impact is embodied in the differencial equation by appling the background electromagnetic field item. The computing accuracy and efficiency would be managed rigorously due to the wide dynamic range of the transient singnal, as well as the twotime inverse laplace and fourier transforms. Otherwise the accumulative error will soar acceptably. The synthetic model experiments shows that the numeric modeling solutions match the analytical results of homogeneous and layered erath responses well, which also proves the effectiveness of the algorithm.
出处
《物探与化探》
CAS
CSCD
2015年第5期1059-1062,共4页
Geophysical and Geochemical Exploration
基金
国家自然科学基金项目(41174104
41472301)
关键词
航空瞬变电磁
2.5维正演
有限单元法
异常场法
响应特征
ATEM
2.5-D forward modeling
finite element method
anomalous field algorithm
response characteristic
