Research on 3D modelling of induced polarization in polarizability anisotropic media

Abstract

To enhance the 3D numerical simulation of the induced polarization method within anisotropic media, our study employs the 2D Fourier transform technique. This technique is utilized to convert the 3D integral of the abnormal potential from the space domain into a 1D integral in the wave number domain. Subsequently, we apply the shape function integration method, which is founded on quadratic interpolation, to resolve the 1D integral equation effectively. This methodology significantly decreases the necessary computational resources and storage while simultaneously harnessing the high efficiency and accuracy of the 1D shape function integration method, as well as the high efficiency of the fast Fourier transform, optimizing the numerical simulation process of the induced polarization method. We validate the accuracy of our algorithmic approach using an equivalent uniform layered model. Furthermore, by employing the sphere model, we conduct a comparison of computation time with the finite element method, thereby demonstrating high efficiency of the proposed algorithm. Utilizing the OpenMP parallel algorithm, we confirm that the proposed algorithm has a high degree of parallelism. We also analyse the differences in the equivalent apparent resistivity and apparent polarizability for various electrical parameters, using a prismatic model as the basis for our analysis. Our results clearly indicate that the anisotropy of the polarizability exerts substantial influence on the observe data. Consequently, the implications of polarizability anisotropy are deemed critical and not be disregarded in the field detection applications.