Numerical Simulation of a Three-Dimensional Strongly Magnetic Material Tensor in Wave Vector Domain
Li Kun, Shi Hui, Wu Yu-cheng*, Kang Chen, Zhao Dong-dong
1. School of Earth Sciences and Technology, Southwest Petroleum University, Chengdu 610500, China
2. State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Chengdu 610500, China.
3. Sinopec Geophysical Corporation Nanfang Branch, Chengdu 610200, Sichuan ,China
4. Exploration and Development Research Institute, PetroChina Southwest Oil & Gasfield Company, Chengdu, 610041,China
5. School of Life and Environmental Sciences, Guilin University of Electronic Technology, Guilin 541000, China.
Abstract:
Demagnetization effects are typically ignored in numerical simulations of magnetic material tensors. However, as the magnetization susceptibility increases, demagnetization effects become more pronounced, thereby gradually decreasing the accuracy of magnetic material tensor simulations. An iterative calculation method based on the wave vector domain is proposed to address this issue in the context of strongly magnetic materials. This method uses the integral equations satisfied by the magnetic potential, employs a twodimensional Fourier transform to convert the three-dimensional integral problem into a one-dimensional problem, and utilizes the shape function method to solve this one-dimensional problem. An iterative tightbinding operator is introduced to ensure the convergence of the algorithm. High-precision numerical simulations of strongly magnetic material tensors are achieved on the basis of this foundation. Comparisons with analytical solutions for strongly magnetic spherical and spherical shell models demonstrate that the proposed method is applicable to a wide range of magnetization susceptibilities. Under normal conditions of strong magnetization (χ = 1SI), the method converges rapidly and closely matches the analytical solutions. At higher magnetization levels (1SI < χ ≤ 100SI), it stabilizes after a fi nite number of iterations while maintaining high computational accuracy. At extremely strong magnetization levels (100SI < χ ≤ 500SI), the convergence speed and computational accuracy of the method decrease relatively, but it still satisfi es the requirements for computational accuracy.