3D full-time anisotropic TEM modelling using a mixed BDF2/SAI method. (28th December 2021)

• Record Type: Journal Article
• Authors: Huang, Kui; Lu, Kailiang; Zhou, Jianmei; Li, Xiu; Meng, Lifei
• Is Part Of: Journal of geophysics and engineering. Volume 18:Number 6(2021)
• Page Start: 995
• Page End: 1006
• Publication Date: 2021-12-28
• Subjects: TEM -- numerical modelling -- matrix exponential -- BDF2 algorithm -- anisotropy; Geophysics -- Periodicals; Prospecting -- Geophysical methods -- Periodicals; Engineering -- Periodicals; 622.1505
• DOI: 10.1093/jge/gxab066
• Languages: English
• ISSNs: 1742-2132

Abstract

Abstract: Transient electromagnetic (TEM) data are affected by resistivity anisotropy, which should be considered in 3D modelling. The influence of anisotropy on full-time response is the main focus of this research. For spatial discretisation of an anisotropic model, the mimetic finite volume approach was applied. The accuracy of the shift-and-invert (SAI) Krylov subspace approach and the two-step backward differentiation formula (BDF2) for modelling 3D full-time electromagnetic data has been demonstrated. However, both algorithms require time-consuming calculations. The SAI technique requires a number of projection subspace constructions, whereas the BDF2 algorithm necessitates numerous coefficient matrix decompositions. We proposed a novel mixed BDF2/SAI algorithm in this paper, which combines the advantages of the two algorithms. The on-time response is computed using BDF2, while the off-time response is computed using the SAI-Krylov subspace method. The forward results of a 1D model with a half-sine waveform demonstrated that the new algorithm is accurate and faster than both the BDF2 algorithm and the SAI algorithm. During the full-time period, the forward results of a 3D anisotropic model with half-sine waveform show that abnormal responses can be induced. It was shown that the relative abnormal of ${{{\bf b}}_{\boldsymbol{z}}}$ is higher during the on-time period, while the relative abnormal of $\partial {{{\bf b}}_{\boldsymbol{z}}}/\partial t$ is higher during the off-time period. Furthermore, the change in relative anomaly is more obvious as the anisotropic block rotates around the x -axis. And the larger the rotation angle, the larger the relative anomaly.