论文部分内容阅读
本研究提出一种数值方法,计算频率域中由偶极源在二维地质结构体上产生的电磁响应,偶极源在地面激发,假定电磁参数(σ、ε、μ)沿走向方向没有变化,仅在垂直走向的平面变化,此类问题即为偶极源二维半问题。人工源电磁法(CSEM)作为勘探地球物理主要方法之一,可以很快采集到高质量的电磁数据,是一种有效的勘探手段。然而,人工源电磁法采用分离一次场和二次场的常规做法,对复杂地质结构是很难做到的,因此其对复杂的二维三维结构体很难模拟。另外,在实际工作中,地形的影响往往也要考虑,这更增加了模拟的困难。为了解决上述困难,通过采用一种狄拉克δ函数代表源,则不需要分离一次场和二次场;采用等参有限元法可以很好地模拟复杂模型和地形影响;求积分时通过构造合适的高阶插值函数,从而提高求解精度。整个做法的思路是:从麦克斯韦方程组出发,对电场磁场强度各个分量沿走向方向做傅立叶变换,将空间域变换到波数域,用有限元法对研究区域进行剖分,运用等参有限元法模拟实际地形影响,然后根据微分方程和边界条件,采用Galerkin方法推出有限元方程,求波数域电磁场各分量的弱解(数值解),最后通过傅立叶反变换得到频率域的解。最后用模型对这种方法进行验证。