[发明专利]一种三维各向异性射频大地电磁自适应有限元正演方法

摘要

本发明公开了一种三维各向异性射频大地电磁自适应有限元正演方法，包括：构建表示求解区域的三维地电几何模型，并获取不同区域的控制参数；将三维地电几何模型剖分为若干个互不相交的四面体单元；获取三维各向异性射频大地电磁对应的积分弱形式，并构建稀疏线性方程组，再在地表引入两种不同极化方向的入射场获得两种稀疏线性方程组中的不同右端项，并求解所述两种稀疏线性方程组得到两种极化方向上的每条棱边上的近似电场切向分量；最后基于两种极化方向上每条棱边的近似电场切向分量获取测点处的电场、磁场、阻抗张量、磁倾子矢量、视电阻率、相位响应中一个或多个的参数组合。该方法基于自适应有限元法实现3D RMT各向异性正演模拟。

主权项

1.一种三维各向异性射频大地电磁自适应有限元正演方法，其特征在于：包括如下步骤：S1：构建表示求解区域的三维地电几何模型，并获取三维地电几何模型中不同区域的控制参数；所述控制参数为各个区域的电导率主值参数、相对介电常数主值参数、对应的主轴坐标系与观测坐标系之间的旋转角；S2：将三维地电几何模型剖分为若干个互不相交的四面体单元；其中，一个四面体单元仅属于三维地电几何模型中的一个区域；S3：根据麦克斯韦方程组和Garlerkin有限元原理获取三维各向异性射频大地电磁对应的积分弱形式，并基于所述积分弱形式、1st order Nédélec基函数、四面体单元每条棱边的近似电场切向分量构建稀疏线性方程组，S4：在地表引入两种不同极化方向的入射场得到所述稀疏线性方程组中的两种不同右端项，并求解所述两种含不同边界条件值的稀疏线性方程组获得两种极化方向上分别对应的每条棱边上的近似电场切向分量；S5：基于两种极化方向上每条棱边的近似电场切向分量获取测点处的目标参量，所述目标参量为：电场、磁场、阻抗张量、磁倾子矢量、视电阻率、相位响应中一个或多个的组合。