2.5维电磁响应的有限元模拟与波数取值研究
利用有限元法实现了任意方向偶极子源在二维介质中频率域电磁响应的数值模拟,研究了波数取值对模拟结果的影响.通过对构造走向的Fourier变换,将全三维电磁问题,转化为一系列二维问题,并在波数域求解,极大地减小了计算工作量,导出了波数域耦合适用于二维电性介质中任何方向电或磁偶极子响应计算的电磁场方程.针对每个给定的波数,上述耦合电磁场方程用等参有限元方法在x-z平面内求解.采用Fourier逆变换,将波数域解积分,得到空间域电磁场.针对电磁模拟计算中,源点的奇异性,采用具有一定面积的伪δ函数表达源电流分布,使数值解精度得以提高.另外,采用等参有限元,使地下复杂地质体得到准确表达.利用不同波数值对均匀介质与层状介质的模拟结果与解析解的对比,验证了算法的正确性与精度.利用层状介质模型的解析解与数值计算结果的对比,分析了波数的优化取值范围及取值点数对数值模拟结果的影响,考察了算法对非均匀介质的适应性.
孙文博,SUN Wen-bo(中国石油大学,资源与信息学院,地球物理系,北京,102249)
刊 名： 物探化探计算技术 ISTIC 英文刊名： COMPUTING TECHNIQUES FOR GEOPHYSICAL AND GEOCHEMICAL EXPLORATION 年，卷(期)： 2008 30(2) 分类号： P631.3+25 关键词： 2.5维 波数范围 电磁响应模拟 等参有限元 伪δ函数