[发明专利]基于地空分解策略的大地电磁正演方法及装置、存储介质

权利要求书

1.一种基于地空分解策略的大地电磁正演方法，其特征在于：包括如下步骤：

步骤1：构建求解区域三维MT的地电几何模型；

步骤2：将所述地电几何模型离散为一系列不重叠的四面体单元；

步骤3：获取基于地空分解策略的A-Φ-Ψ耦合势三维MT正演满足的积分弱解形式，并采用节点基函数和四面体单元的节点对所述积分弱解形式进行基于节点离散的转换得到稀疏线性方程组，其中，A表示磁场矢量位、Φ表示电场标量位，Ψ表示磁场标量位，节点为离散四面体单元之间的连接枢纽；

步骤4：加载入射场置于稀疏线性方程组的右端项并求解所述稀疏线性方程得到求解区域中四面体单元各个节点的磁场矢量位A、电场标量位Φ，磁场标量位Ψ，其中，空气区域内节点的磁场矢量位A、电场标量位Φ为0，地下区域内节点的磁场标量位Ψ为0；

步骤5：根据所述磁场矢量位A、电场标量位Φ，磁场标量位Ψ计算出观测点的电场、磁场一个或多个参数的组合。

2.根据权利要求1所述的方法，其特征在于：步骤3中所述积分弱解形式如下：

式中，V表示测试函数，E表示电场，B为磁场,Ω表示求解区域，表示求解区域的边界，Ω₀为空气空间区域，Ω₁为地下空间区域，Γ₀表示空气层的上顶面，Γ₁表示地下区域的下底面，Γ为地下区域和空气区域的分界面，n表示边界外法向量，dv、ds分别表示体积微分元和面积微分元，v表示体积、s表示面积，表示阻抗率，表示导纳率，μ为磁导率。

3.根据权利要求1所述的方法，其特征在于：步骤3中所述稀疏线性方程组如下：

式中，均为稀疏的子矩阵，均为地空分界面的面积分项形成的子矩阵,均表示零矩阵，B₀,B₁,B₂,B₃,B₄均为稀疏线性方程组的右端项，X₀、X₁、X₂、X₃、X₄分别为稀疏线性方程组中的未知数矩阵，且由求解区域内所有节点的磁场矢量位A、电场标量位Φ，磁场标量位Ψ构成。