基于对偶加权后验误差估计的2.5维直流电阻率自适应有限元正演  

2.5D DIRECT CURRENT RESISTIVITY ADAPTIVE FINITE-ELEMENT NUMERICAL MODELING BASED ON DUAL WEIGHTED POSTERIORI ERROR ESTIMATION

作　　者：严波[1] 刘颖[1] 叶益信[1] 

机构地区：[1]中国海洋大学海底科学与探测技术教育部重点实验室,山东青岛266100

出　　处：《物探与化探》2014年第1期145-150,共6页Geophysical and Geochemical Exploration

基　　金：国家自然科学基金(41204055;41130420)

摘　　要：采用易于模拟复杂地形起伏和倾斜界面的非结构三角单元剖分网格,并利用对偶加权后验误差估计指导网格自动细化过程,实现了2.5维直流电阻率法自适应有限元数值模拟。在实例模型分析中,分别计算了层状模型和垂直岩脉模型的直流电阻率响应,并与其解析解进行了比较。对比结果表明,该算法所得数值解精度很高,解的相对误差小于0.5%。最后,计算了起伏地形2.5维地电模型视电阻率异常,并利用比较法进行了地形改正。地形改正结果与水平地形时的结果对比表明,比较法可以较好地消除地形影响,突出局部地质体异常。2.5D DC resistivity modeling was realized by adaptive finite-element numerical simulation algorithm. This algorithm uses un-structured 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 to-pography by using the comparative method. The contrast between topographic correction results and level terrain results show that the com-parative method is effective in eliminating the terrain effect and giving prominence to local geological anomalies.

关 键 词：直流电阻率 自适应有限元法 地形改正 2 5维 后验误差估计 

分 类 号：P631[天文地球—地质矿产勘探]