Two-dimensional forward modeling of magnetotelluric using a Wavelet-Galerkin method
Magnetotelluric method is widely used in the exploration of deep geological structure,mineral resources,hydrogeology and so on because of its large exploration depth and low cost.Efficient forward modeling is the key to fast and stable magnetotelluric inversion.Therefore,based on the properties of Daubechies wavelet and the scale function of Daubechies wavelet as the basis function of Galerkin method,a two-dimensional forward algorithm for magnetotellurics based on wavelet is constructed in this paper.Firstly,based on the partial differential equations satisfying the magnetotelluric field,the wavelet Galerkin method equations for two-dimensional forward magnetotelluric field are derived in detail,and the bistable conjugate gradient algorithm of incomplete LU decomposition is used to solve the equations,which improves the efficiency of solution.Finally,a typical geoelectric model is designed for forward calculation and compared with the finite element method.The comparison results show that the wavelet Galerkin method is feasible and effective for two-dimensional magnetotelluric forward modeling.The algorithm presented in this paper has good universality and can be extended to other electromagnetic numerical simulation problems.
NETL
NSTL
万方数据
