1.一种基于混合网格的大地电磁高阶有限元正演模拟方法,其特征在于:所述模拟方法如下:采用Delaunay三角算法将地球表面自由划分为三角形网格,在三角形网格的垂直方向上延伸三角形网格来形成三棱柱单元;与三棱柱单元相邻的区域被离散成四面体单元,从而通过三角形平面将三棱柱单元的上下界面耦合在起来,四面体单元的大小基于耦合连接的三角形的大小,并逐渐向远离该三角形一侧增长;利用构建的混合网格离散地下电性模型,通过高阶形函数进行有限元的单元分析,离散大地电磁波动方程,形成大地电磁高阶有限元系统;求解该混合网格有限元系统,获得大地电磁响应信号。
2.根据权利要求1所述的一种基于混合网格的大地电磁高阶有限元正演模拟方法,其特征在于:所述三角形网格延伸的距离是趋附深度的1到2倍,每个细分层的垂直尺寸应小于趋附深度的1/3。
3.根据权利要求1所述的一种基于混合网格的大地电磁高阶有限元正演模拟方法,其特征在于:采用二阶形函数离散三棱柱单元,所述三棱柱单元的形函数由矢量形函数和标量节点形函数组成,二阶三棱柱单元有15个节点形函数;
其中三棱柱单元6个顶点处的节点形函数为:
三角形各边中点处的节点形函数为:
矩形各边中点处的节点形函数为:
矢量形状函数是从节点形状函数派生出来的;它们的关系为:
其中边i‑j与节点j相关联。