2.5D induced polarization forward modeling using the adaptive finite-element method

The conventional finite-element (FE) method often uses a structured mesh, which is designed according to the user’s experience, and it is not sufficiently accurate and flexible to accommodate complex structures such as dipping interfaces and rough topography. We present an adaptive FE method for 2.5...

Full description

Saved in:
Bibliographic Details
Published in:Applied geophysics 2014-12, Vol.11 (4), p.500-507
Main Authors: Ye, Yi-Xin, Li, Yu-Guo, Deng, Ju-Zhi, Li, Ze-Lin
Format: Article
Language:English
Subjects:
Adaptive algorithms
Algorithms
Earth and Environmental Science
Earth Sciences
Finite element analysis
Finite element method
Geophysics
Geophysics/Geodesy
Geotechnical Engineering & Applied Earth Sciences
Induced polarization
Mathematical analysis
Mathematical models
Mesh generation
Receivers
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:The conventional finite-element (FE) method often uses a structured mesh, which is designed according to the user’s experience, and it is not sufficiently accurate and flexible to accommodate complex structures such as dipping interfaces and rough topography. We present an adaptive FE method for 2.5D forward modeling of induced polarization (IP). In the presented method, an unstructured triangulation mesh that allows for local mesh refinement and flexible description of arbitrary model geometries is used. Furthermore, the mesh refinement process is guided by dual error estimate weighting to bias the refinement towards elements that affect the solution at the receiver locations. After the final mesh is generated, the Jacobian matrix is used to obtain the IP response on 2D structure models. We validate the adaptive FE algorithm using a vertical contact model. The validation shows that the elements near the receivers are highly refined and the average relative error of the potentials converges to 0.4 % and 1.2 % for the IP response. This suggests that the numerical solution of the adaptive FE algorithm converges to an accurate solution with the refined mesh. Finally, the accuracy and flexibility of the adaptive FE procedure are also validated using more complex models.
ISSN:1672-7975
1993-0658
DOI:10.1007/s11770-014-0455-z