Three-Dimensional DC Anisotropic Resistivity Modeling Using a Method in the Mixed Space-Wavenumber Domain

A newly developed method in the mixed space-wavenumber domain is proposed for 3D DC resistivity modeling in anisotropic structures. Unlike the traditional method for numerical simulation of DC resistivity, this method performs a 2D Fourier transform directly on the partial differential equation of DC abnormal potential along the horizontal direction, so that the horizontal direction is converted into the wavenumber domain while the vertical direction remains in the spatial domain. In this way, the 3D partial differential equation in the spatial domain is transformed into a 1D ordinary differential equation in different wave numbers. Therefore, a large-scale 3D numerical simulation problem is decomposed into many 1D numerical simulation problems. The 1D finite element method is used to solve the equations, and a contraction operator is used for iterative calculation to obtain a highly accurate solution. The accuracy of the algorithm is verified by comparison with the goal-oriented adaptive finite element method. Compared with the COMSOL Multiphysics commercial finite element simulation software, the algorithm in this paper has high computational efficiency. To identify the electrical anisotropy, we propose a new measurement and then use it to study the response characteristics of a half-space model with two prisms. The test results show that this measurement can accurately reflect the underground electrical anisotropy. Finally, the code is applied to the numerical simulation of a vertical-borehole survey and shows that the significant information on the hydraulically induced fracturing system can be reflected well by the electric anisotropy.