论文部分内容阅读
大地电磁正演领域,目前流行的数值计算方法包括有限差分法和有限元法,该类数值计算方法都是基于网格插值求解计算域物理场,对复杂地电模型的网格剖分比较繁琐,计算精度严重依赖于单元剖分的形状和大小。为克服这些困难,无网格法作为有网格法的一种重要补充和发展具有重要的研究意义和实用价值。本文重点研究基于全局弱式无网格法的大地电磁二维正演问题。 首先从大地电磁场二维问题所满足的Maxwell方程组出发,利用变分问题与边值问题的等价性,对大地电磁二维问题满足的微分方程进行推导并得到其相应的变分形式。利用一系列离散场节点对求解域进行离散,采用移动最小二乘法(MLS)构造积分点支持域形函数。根据Galerkin法原理,推导出与大地电磁法二维边值问题对应的无网格全局弱式矩阵表达式。采用含背景网格的高斯积分将计算域离散化,求得最终总体矩阵的离散表达式。 后续通过多个一维层状模型数值模拟计算,研究了相关重要参数,如高斯点数目,无量纲支持域尺寸等对计算精度及计算效率的影响,并确定其参数最优值,选取高斯点数量为2~3、无量纲支持域尺寸a为1.1~1.3。确定最优参数值后,分别用有限元法和全局弱式无网格法计算几个典型一维模型和二维模型,通过对比计算结果及计算精度,验证本文中提出的全局弱式无网格法的有效性及在处理复杂模型上的优越性。 最后,根据冬瓜山铜矿实际情况设计了二维地质模型,并用本文提出的全局弱式无网格法对其进行数值模拟计算,其计算结果对低阻矿体实际形态反应较好,进一步证明了无网格算法在地球物理电磁波数值模拟中具有重要的意义。