Three-Dimensional Forward Modeling of Transient Electromagnetic Method Considering Induced Polarization Effect Based on Spectral Element Method

The transient electromagnetic method (TEM) is widely used in the exploration of mineral, petroleum, and geothermal resources due to its sensitivity to low-resistivity bodies, limited site constraints, and strong resistance to interference. In practical applications, the TEM often uses a long wire source instead of an idealized horizontal electric dipole (HED) source as the excitation source. This is due to the complex external conditions and the relatively large distance between the receiving zone and the transmitter source. Compared to the HED, the long wire source can provide a larger excitation current, generating stronger signals to meet the requirements of a higher signal-to-noise ratio or deeper exploration. It also produces longer-duration signals, thereby providing better resolution. Additionally, for the interpretation of TEM data, three-dimensional forward modeling plays a crucial role. However, the mature traditional TEM forward method is based on a simple, sometimes inappropriate model, as it is well established that the induced polarization (IP) effect is widely present in the deep earth, especially in oil and gas reservoirs. The presence of the IP effect results in negative responses in field data that do not conform to the traditional theoretical decay law of TEM, which can significantly impact data processing and inversion results. To address this issue, a TEM forward modeling method considering the IP effect based on the spectral element method (SEM) has been developed in this study. Firstly, starting from the Helmholtz equation satisfied by the time domain electric field, we introduce the Debye model with polarization information into the forward modeling by utilizing the differential form of Ohm's law. As a result, we derive the boundary value problem for the time domain electric field that considers the induced polarization effect. Using Gauss-Lobatto-Legendre (GLL) polynomials as the basis functions, the SEM is employed to discretize the governing equations at each time step and obtain spectral element discretization equations. Then, temporal discretization equations are derived using the second-order backward Euler formula, and the linear system of equations is solved using the Pardiso direct solver. Finally, the electromagnetic responses at any time channel are obtained via SEM interpolation and numerical integration, thereby achieving three-dimensional TEM forward modeling considering the IP effect. The results indicate that this method can effectively reflect the spatial distribution of polarizable subsurface media. It provides valuable references for studying the polarization parameters of subsurface media and performing a three-dimensional inversion of TEM data considering the induced polarization effect.