基于非均一场源的球坐标大地电磁模拟方法

大地电磁测深理论与数据处理解释均假定平面电磁波垂直入射地下空间,但随着研究尺度的逐渐扩大,使得因地球弧度产生的影响难以忽略.此时,传统笛卡尔坐标体系及平面波场源不再适用于大尺度的大地电磁数据正反演解释.本文提出并实现了一种基于球坐标系的大地电磁交错网格有限差分三维正演算法,并对电场进行极向-环向分解,结合球谐函数和贝塞尔函数构建了可取代平面波的场源模型.首先利用经纬度信息构建三维地电模型,将场源设置于模型空间正上方,然后通过直接求解球坐标系下麦克斯韦方程来获得大地电磁响应.在此基础上,本文设计了球坐标下具有不同分辨率的多个三维地电模型,阐述了由球体模型到笛卡尔模型的转换方法,详细对比了两种坐标体系在计算效率、所求得的电场和视电阻率方面的差异.结果表明二者差异度主要与电性横向分布和地图投影方法有关,与周期并不存在明显的单调递增关系.