Three-dimensional direct current resistivity forward modeling based on the hp-adaptive finite element method

Ning Zhao,Yahang Shen,Ce Qin,Xuben Wang

DOI: https://doi.org/10.1016/j.cageo.2020.104566

2020-11-01

Abstract:<p>In recent years, the adaptive finite element method has been used for direct current resistivity forward modeling to improve the accuracy of numerical solutions. The accuracy of adaptive finite element solutions is mainly affected by two factors: the cell size (<em>h</em>) and the order of the shape function (<em>p</em>). To further improve the accuracy of the adaptive finite element solutions and keep low computational costs, we present a <em>hp</em>-adaptive finite element algorithm, combining <em>h</em>-adaptive and <em>p</em>-adaptive, for three-dimensional direct current forward modeling. The adaptive mesh refinement is guided by a posteriori error estimator and a smoothness estimator of the solution. In the adaptive process, for a cell with a large error, if the smoothness of the solution is low, the cell will be refined; otherwise, the order of the shape function on the cell will be increased by one order. To obtain high-precision finite element solutions for complex models with topographies the unstructured grids are adopted. The octree-based mesh refinement method is used to refine the cells, and the shape functions with arbitrary orders in three-dimensional space are generated by using tensor products of one-dimensional polynomials. A two-layer model is used to verify the correctness of the algorithm. The comparison between the convergence rate of global, <em>h</em>-adaptive, and <em>hp</em>-adaptive refinement algorithms indicates that our algorithm can provide the most accurate solution with the lowest computation costs. It shows exponentially convergence rate. Finally, a topography model with an abnormal body and a complex abnormal model are used to verify the robustness of our algorithm.</p>

geosciences, multidisciplinary,computer science, interdisciplinary applications

What problem does this paper attempt to address?

The problem this paper attempts to address is: how to improve the accuracy of 3D direct current resistivity forward modeling while maintaining low computational costs. Specifically, the authors propose an hp-adaptive finite element algorithm that combines h-adaptive and p-adaptive methods to enhance the accuracy of numerical solutions. Traditional h-adaptive methods mainly improve accuracy by refining the mesh, but this approach is not effective when the solution is sufficiently smooth in certain regions. On the other hand, p-adaptive methods improve accuracy by increasing the order of the shape functions, but they are limited in regions where the solution is not smooth. Therefore, the hp-adaptive method proposed in this paper can select the appropriate strategy based on the smoothness of the solution, achieving higher accuracy and lower computational costs.

To verify the effectiveness of the algorithm, the authors conducted the following work:
1. **Correctness Verification**: The algorithm's correctness was verified using a two-layer model, and the results showed that the numerical solution gradually converged to the analytical solution, with a significant reduction in relative error.
2. **Convergence Rate Comparison**: The convergence rates of global, h-adaptive, and hp-adaptive algorithms were compared using a uniform half-space model. The results indicated that the hp-adaptive algorithm had a higher convergence rate and lower computational cost at high degrees of freedom.
3. **Terrain Model Verification**: The robustness of the algorithm was verified using a terrain model containing an anomalous body, and the results showed that the algorithm could handle complex terrain models well.
4. **Complex Anomaly Model Testing**: The effectiveness of the algorithm was further tested using a complex anomaly model, and the results showed that the numerical solution gradually converged to the analytical solution, with a significant reduction in the maximum relative error.

In summary, this paper aims to resolve the contradiction between accuracy and computational cost in 3D direct current resistivity forward modeling by proposing a new hp-adaptive finite element algorithm, providing a more efficient and accurate numerical simulation tool for fields such as geological exploration.