无网格局部Petrov-Galerkin法大地电磁二维正演

Magnetotelluric two-dimensional forward modeling using meshless local Petrov-Galerkin method

将无网格局部Petrov-Galerkin算法用于大地电磁二维正演。介绍了该方法的基本原理;从大地电磁二维边值问题出发,利用子域法详细推导了与之对应的局部Petrov-Galerkin弱式方程,并用高斯积分法将其离散化。论述了无网格局部Petrov-Galerkin法较无单元Galerkin法及有限元法的优缺点,最后通过二维模型的计算验证了算法的有效性。

This Paper devotes meshless local Petrov-Galerkin method to magnetotelluric two-dimensional forwardmodeling. Its basic principle was introduced. From magnetolelluric 2-D boundary values, using subdomain method,the corresponding weak equation of local Petrov-Galerkin is derived and discreted by gauss integral method. Theadvantages and disadvantages of meshless local Petrov-Galerkin method are discussed in comparison withelement-free Galerkin method and finite element method. At last, the validity of the algorithm is verified bytwo-dimensional model calculations.