Meshfree magnetotelluric modelling
Meshfree magnetotelluric modelling
复制标题
DOI：
10.1093/gji/ggu207
复制
发表时间：
2014-08
影响因子：
2.8
通讯作者：
J. Wittke;B. Tezkan
中科院分区：
文献类型：
--
作者：
J. Wittke;B. Tezkan
与有限元或有限差分等传统数值方法相比，我们提出了一种二维大地电磁正演数值模拟的新方法。用于求解偏微分方程的方法基于无网格技术，不需要底层网格或网格。我们使用无网格局部 Petrov-Galerkin (MLPG) 方法与径向基函数相结合来模拟给定电导率模型对平面波源的响应。我们将无网格解决方案与已知的模拟程序和简单的解析解决方案进行了比较。此外，我们从实施和稳定性方面讨论了新的大地电磁建模方法。首先，我们用简单的半空间电导率模型研究了新方法的收敛性和离散化误差。然后，我们将简单的二维电导率模型的无网格模拟结果与众所周知的有限元程序的结果进行比较。最后，我们提供了使用无网格方法计算的平滑电导率模型。即使节点随机分布，建模结果也与有限元方法获得的结果非常吻合。
We present a new approach for 2-D magnetotelluric forward numerical modelling in contrast to traditional numerical methods like finite elements or finite differences. The method used for solving the partial differential equations is based on a mesh-free technique which does not need an underlaying mesh or grid. We use the Meshless Local Petrov–Galerkin (MLPG) method in combination with radial basis functions to simulate the response of a given conductivity model to a plane-wave source.We compare the mesh-free solution with known simulation programs and simple analytical solutions. Furthermore, we discuss the new magnetotelluric modelling method in terms of implementation and stability. First, we study the convergence and discretization errors of the new method with a simple half-space conductivity model. Then we compare our mesh-free simulation results with simple 2-D conductivity models with the results of a well-known finite element program. In the end, we provide a smooth conductivity model calculated with the mesh-free approach. The modelling results, even with randomly distributed nodes, are in a good agreement with those obtained by the finite element method.