Three-dimensional Forward Calculation for Loop Source Transient Electromagnetic Method Based on Electric Field Helmholtz Equation
Li Jian-Hui,Hu Xiang-Yun,Zeng Si-Hong,Lu Jin-Ge,Huo Guang-Pu,Han Bo,Peng Rong-Hua
DOI: https://doi.org/10.6038/cjg20131228
2013-01-01
Abstract:Forward calculation is the base of parameters selection in field work, data processing, and interpretation; particularly the accurate, stable and efficient three-dimensional forward solver is very important. In this paper, based on the strategy that solving the Laplace-domain electrical field firstly, then acquiring the time-domain magnetic field by Gaver-Stehfest algorithm, we employed staggered finite difference (SFD) and finite volume (FV) to discretize the electrical field Helmholtz equation for three-dimensional forward calculation of loop source transient electromagnetic method. Comparing with the numerical solutions of integral equation (IE), finite difference in time domain (FDTD), vector finite element method (VFEM) and spectral Lanczos decomposition method (SLDM) for low resistivity body in homogeneous half-space, we validated the SFD and FV algorithms. Because the unknown electrical field is arranged on the edge of brick for SFD, FV and VFEM based on brick element, so these three numerical algorithms could utilize the same method for coding the unknown electrical field, computing the background field and post processing. However, comparing with the VFEM, the coefficient matrixes produced by SED and FV are sparser, so these two algorithms are more efficient. For low resistivity horizontal sheet model, we use SFD, FV and 1D analytical method to compute the induced electromagnetic force (EMF). The results show that SFD has higher accuracy than FV, and the contours of EMF for these three algorithms coincide with each other very well although the range of EMF has fine distinctions.