Abstract：2.5D DC resistivity modeling was realized by adaptive finite-element numerical simulation algorithm. This algorithm uses unstructured triangular element that is easy to simulate complex topography and inclined interface to subdivide mesh, and utilizes dual weighted posteriori error estimation guide mesh refining process. In the instance model analysis, the authors calculated the DC resistivity response of the layered model and the vertical dike model, which were compared with their analytical solutions respectively. The results show that the numerical solution for this algorithm has high accuracy,and the relative error of the numerical solution is less than 5 percent. Finally, the authors calculated 2.5D geoelectric model apparent resistivity anomaly of undulate topography and corrected the undulate topography by using the comparative method. The contrast between topographic correction results and level terrain results show that the comparative method is effective in eliminating the terrain effect and giving prominence to local geological anomalies.
严波, 刘颖, 叶益信. 基于对偶加权后验误差估计的2.5维直流电阻率自适应有限元正演[J]. 物探与化探, 2014, 38(1): 145-150.
YAN Bo, LIU Ying, YE Yi-xin. 2.5D DIRECT CURRENT RESISTIVITY ADAPTIVE FINITE-ELEMENT NUMERICAL MODELING BASED ON DUAL WEIGHTED POSTERIORI ERROR ESTIMATION. Geophysical and Geochemical Exploration, 2014, 38(1): 145-150.
[1] 徐世浙.地球物理中的有限单元法[M].北京:科学出版社,1994:183-194. [2] 熊斌,阮百尧,罗延钟.复杂地形条件下直流电阻率异常三维数值模拟研究[J].地质与勘探,2003,39(4):60-64. [3] 强建科,罗延钟.三维地形直流电阻率有限元法模拟[J].地球物理学报,2007,50(5):1606-1613. [4] Sasaki Y. 3D resistivity inversion using the finite-element method[J].Geophysics, 1994, 59:1839-1848. [5] Tang J T, Wang F Y, Ren Z Y. 2.5D DC resistivity modeling by adaptive finite-element method with unstructured triangulation[J].Geophs, 2010, 53(3):708-716. [6] 王飞燕. 基于非结构化网格的2.5D直流电祖率法自适应有限元数值模拟.长沙:中南大学,2009. [7] Ovall J S. Duality-based adaptive refinement for elliptic PDEs.San Diego:University of California, 2004. [8] Ovall J S. Asymptotically exact functional error estimators based on super convergent gradient recovery[J].Numeriche Mathematik,2006,102:543-558. [9] 罗延钟,张贵青.电子计算机在电法勘探中的应用[M].武汉:武汉地质学院出版社,1987:67-69. [10] 阮百尧,徐世浙.电导率分块线性变化二维地电断面电阻率测深有限元数值模拟[J].地球科学:中国地质大学学报,1998,23(3):303-307. [11] Li Y G, Kerry K. 2D marine controlled-source electromagnetic modeling:Part 1-An adaptive finite-element algorithm[J].Geophysics,2007,72(2):51-62. [12] Li Y G, Josef P. Adaptive finite element modeling of two-dimensional magnetotelluric fields in general anisotropic media[J].Geophys J Int, 2008,175:942-954. [13] Bank R E, Xu J C.Asymptotically exact a posteriori error estimators, Part Ⅱ:General unstructured grids[J].SIAM Journal on numerical analysis, 2003,41:2313-2332. [14] Babuska I, Ihlenburg F, Stroubousli T, et al. A posteriori error estimation for finite element solution of Helmholtz equation.ⅡEstimation of the pollution error[J].International Journal for Numerical Methods in Engineering, 1997, 41:3883-3900. [15] Oden J, Feng Y, Prughomme S. Local and pollution error estimation for stokesian flows[J].International Journal for Numerical Methods in Fluids, 1998, 27:33-39. [16] 刘兴国.电法勘探原理与方法[M].北京:地质出版社,2003:48-52. [17] Pidlisecky A, Knight R. FW2_5D:A MATLAB 2.5D electrical resistivity modeling code.Computers & Geosciences, 2008, 34:1645-803. [18] Telford W, Geldart L P, Sheriff R E. Applied geophysics second edition[M].Cambridge University Press, 1990:559-561. [19] 刘海飞,阮百尧,吕玉增.直流激电测深二维反演的若干问题研究[J].物探与化探,2007,31(1):47-50. [20] 段本春, 朱永盛.最优化波数在点源二维有限元正演计算中的应用[J].物探与化探,1994,18(3):186-191. [21] 李金铭.地电场与电法勘探[M].北京:地质出版社,2005:185-190.
