Stable Calculation of Responses for a Plane Electromagnetic Wave in an Anisotropic Layered Media for the Modeling of Magnetotelluric Data
0
조회
0
다운로드
0
인용
0
참고문헌
키워드
초록
Research in multi-dimensional magnetotelluric (MT) modeling and inversion, particularly within anisotropic media, is gaining substantial interest. One-dimensional (1-D) modeling is utilized as boundary conditions necessary for accurate 2-D or 3-D anisotropic modeling. Based on Maxwell’s equations and interface boundary conditions, we developed a new algorithm to calculate electromag- netic fields for both 1-D isotropic and anisotropic layered models. However, this algorithm encounters stability issues at considerable depths or within exceptionally thick layers. To address these challenges, we introduced two stabilization strategies. The first strategy alters the exponential terms in the plane-wave solutions to mitigate instabilities, while the second evaluates each layer’s contribution to numerical instability, replacing problematic thick layers with an infinite half-space. Through numerical tests on a four-layer anisotropic model, these strategies proved effective, ensuring stable computations of electromagnetic fields at any depth and in configurations with thick layers, thereby enhancing the model’s applicability for in-depth anisotropic analyses.
Abstract
Research in multi-dimensional magnetotelluric (MT) modeling and inversion, particularly within anisotropic media, is gaining substantial interest. One-dimensional (1-D) modeling is utilized as boundary conditions necessary for accurate 2-D or 3-D anisotropic modeling. Based on Maxwell’s equations and interface boundary conditions, we developed a new algorithm to calculate electromag- netic fields for both 1-D isotropic and anisotropic layered models. However, this algorithm encounters stability issues at considerable depths or within exceptionally thick layers. To address these challenges, we introduced two stabilization strategies. The first strategy alters the exponential terms in the plane-wave solutions to mitigate instabilities, while the second evaluates each layer’s contribution to numerical instability, replacing problematic thick layers with an infinite half-space. Through numerical tests on a four-layer anisotropic model, these strategies proved effective, ensuring stable computations of electromagnetic fields at any depth and in configurations with thick layers, thereby enhancing the model’s applicability for in-depth anisotropic analyses.
참고문헌 (13)
- Bai, N., "3D edge-based and nodal finite element modeling of magnetotelluric in general anisotropic media", Computers & Geosciences, 158, pp.104975~ (2022)
- Cao, X., "3D MT anisotropic inversion based on unstructured finiteelement method", Journal of Environmental and Engineering Geophysics, 26(1), pp.49~60 (2021)
- Key, K., "MARE2DEM: a 2-D inversion code for controlled-source electromagnetic and magnetotelluric data", Geophysical Journal International, 207(1), pp.571~588 (2016)
- Li, Y., "Adaptive finite element modeling of two-dimensional magnetotelluric fields in general anisotropic media", Geophysical Journal International, 175(3), pp.942~954 (2008)
- Liu, Y., "Adaptive finite element modeling of three-dimensional magnetotelluric fields in general anisotropic media", Journal of Applied Geophysics, 151, pp.113~124 (2018)
- O’Brien, D. P., "Electromagnetic fields in an n-layer anisotropic half-space", Geophysics, 32(4), pp.668~677 (1967)
- Pek, J., "Magnetotelluric impedances and parametric sensitivities for 1-D anisotropic layered media", Computers & Geosciences, 28(8), pp.939~950 (2002)
- Reddy, I. K., "Magnetotelluric effect of dipping anisotropies", Geophysical Prospecting, 19(1), pp.84~97 (1971)
- Rong, Z., "Three-dimensional magnetotelluric inversion for arbitrarily anisotropic earth using unstructured tetrahedral discretization", Journal of Geophysical Research: Solid Earth, 127(8), pp.e2021JB023778~ (2022)
- Usui, Y., "3-D inversion of magnetotelluric data using unstructured tetrahedral elements: applicability to data affected by topography", Geophysical Journal International, 202(2), pp.828~849 (2015)
- Usui, Y., "Marine magnetotelluric inversion with an unstructured tetrahedral mesh", Geophysical Journal International, 214(2), pp.952~974 (2018)
- Wannamaker, P. E., "Anisotropy versus heterogeneity in continental solid earth electromagnetic studies: fundamental response characteristics and implications for physicochemical state", Surveys in Geophysics, 26, pp.733~765 (2005)
- Xiao, T., "Three-dimensional magnetotelluric modeling in anisotropic media using edge-based finite element method", Journal of Applied Geophysics, 149, pp.1~9 (2018)