Abstract
A finite-difference approach with non-uniform meshes was presented for simulating magnetotelluric responses in 2D structures. We presented the calculation formula of this scheme from the boundary value problem of electric field and magnetic field, and compared finite-difference solutions with finite-element numerical results and analytical solutions of a 1D model. First, a homogeneous half-space model was tested and the finite-difference approach can provide very good accuracy for 2D magnetotelluric modeling. Then we compared them to the analytical solutions for the two-layered geo-electric model; the relative errors of the apparent resistivity and the impedance phase were both increased when the frequency was increased. To conclude, we compare our finite-difference simulation results with COMMEMI 2D-0 model with the finite-element solutions. Both results are in close agreement to each other. These comparisons can confirm the validity and reliability of our finite-difference algorithm. Moreover, a future project will extend the 2D structures to 3D, where non-uniform meshes should perform especially well.
Highlights
The magnetotelluric method is a passive electromagnetic exploration technique that measures orthogonal components of the electric and magnetic fields on the Earth’s surface [1]
The two-layered model is used as a test example for the comparison of finite difference (FD) numerical solution and analytical solution
We presented the calculation formula of this approach from the boundary value problem of electromagnetic fields
Summary
The magnetotelluric method is a passive electromagnetic exploration technique that measures orthogonal components of the electric and magnetic fields on the Earth’s surface [1]. The magnetotelluric numerical modeling aims to solve the boundary value problem derived from frequency-domain Maxwell’s equations and calculate the spatial distribution of electric and magnetic fields in the subsurface for a given conductivity distribution and a range of frequencies. Numerical modeling approaches such as finite difference (FD), finite element (FE), have been developed and applied as the process of forward modeling for 2D magnetotelluric regularized inversion [4,5,6,7,8]. To verify the accuracy of the FD forward algorithm, the resulting numerical was compared to both an analytical solutions and the FE numerical solutions
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.