3D full-time anisotropic TEM modelling using a mixed BDF2/SAI method

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.