一种各向异性介质大地电磁无网格数值模拟方法

摘要

本发明针对地下空间介质中广泛存在的各向异性问题，提出了一种基于离散节点构造形函数的一种大地电磁无网格数值模拟方法。本发明假设各向异性介质的电性主轴一个垂直于层面方面，另外一个平行于层面方向，构造了各向异性介质模型的电导率张量模型。从各向异性大地电磁边值问题出发，推导了对应无网格径向基点插值法的等价线性方程组，形函数采用复合2次径向基函数（MQ-RBF）构造，通过Krylov子空间的                                               预处理和正则化拟残量极小方法（QMR）方法实现了大型稀疏线性方程组高效，精确的求解。本发明构造的形函数光滑稳定，可以实现电磁法物性参数复杂分布，高精度，自适应的数值模拟。

主权项

一种各向异性介质大地电磁无网格数值模拟方法，其特征在于：包括如下的步骤：1)从各向异性大地电磁边值问题出发，构造等价泛函，推导对应无网格径向基点插值法的等价线性方程组；2)读取当前模型参数，包括频率参数，节点坐标，背景单元，支持域，形状参数，极化模式等；3)对当前极化模式进行判断，若为TE极化模式，则计算区域包含空气层，若为TM极化模式，则不含空气层；4)对所有背景网格进行循环，对背景网格的所有高斯积分点循环，搜索该高斯积分点支持域内的有效节点，计算支持域内节点处的形函数，求取系数矩阵和右端项；5)加载本质边界条件，求解线性方程组，得到各个节点的场值；6)由视电阻率计算公式，代入场值求取地面处各个方向的视电阻率和相位。