3D Anisotropic Forward Modeling for Direct-Current Resistivity Method Based on Multiscale Finite-Element Algorithm
YanFu Qi,Liang Wu,HuaKe Cao,Xiu Li,NaiQuan Sun,ZhiPeng Qi,JianMei Zhou,WenTao Liu
DOI: https://doi.org/10.6038/cjg2023r0077
2024-01-01
Abstract:To handle the forward modeling problem of DC resistivity over anisotropic earth with large-area complex structures,the traditional numerical simulation methods need to generate fine mesh in large area,which leads to a huge amount of calculation.In order to improve the computational efficiency,we develop a 3D forward modeling method for DC resistivity over anisotropic earth based on multiscale finite-element algorithm.By solving the boundary value problem of local partial differential equations in parallel on multiscale meshes,we construct the multiscale basis functions to reflect small-scale electrical information to large scale.Then,by solving the original current conduction problem on a coarse mesh,the forward modeling is accelerated.Aiming at the resonance effect caused by the similarity between the grid size and the heterogeneous-target size,we adopt the oversampling technique to suppress the resonance error and improve its calculation accuracy.We first establish an anisotropic layered model,and verify the accuracy of our algorithm by comparing with the analytical solution.On this basis,we establish two 3D anisotropic resistivity models with complex structures.By comparing their forward modeling results obtained by our multiscale method with those obtained by the traditional finite-element method on fine mesh,we find that the proposed algorithm has higher computational efficiency while ensuring the calculation accuracy.Finally,we further analyze the influence of oversampling technique and coarse/fine ratio of the multiscale mesh on both the numerical accuracy and computational efficiency.It can provide guidance for setting optimal parameters of multiscale finite-element algorithm.