三维地电断面激发极化法有限元数值模拟  

Finite Element Method for IP Modeling on 3-D Geoelectric Section

作　　者：黄俊革 阮百尧[2] 鲍光淑[1] 

机构地区：[1]中南大学信息物理学院,湖南长沙410083 [2]桂林工学院资源与环境工程系,广西桂林541004

出　　处：《地球科学（中国地质大学学报）》2003年第3期323-326,共4页Earth Science-Journal of China University of Geosciences

基　　金：国家自然科学基金 (No .40 0 6 40 0 1);广西自然科学基金匹配项目;教育部优秀青年教师基金项目;广西培养新世纪学术和技术带头人专项基金 .

摘　　要：用有限单元法求解三维地电断面激发极化法正演模拟算法 .首先给出了三维构造中点源电场异常电位的边值问题与变分问题 ,简化无穷远边界上的边界条件以提高计算速度及计算精度 .以此为基础计算视电阻率对模型电阻率的偏导数矩阵 ,并进行三维地电断面激发极化法正演模拟 ,与等效电阻率法相比节约了计算时间 .对几例较典型的地电模型进行计算 ,结果表明本方法是正确可行的 .FEM is used to modeling IP sounding on 3-D geoelectric section. First, this paper presents the 3-D boundary value and variation equation of point-source electrical field on a 3-D structure. Then, the boundary condition for the infinite boundary is simplified to improve calculation speed and precision. On this basis, the partial derivatives matrix for the ratio between the apparent resistivity and the model resistivity parameters is calculated, and at the same time, the IP on 3-D geoelectric section is positively simulated. Compared with equivalent resistivity method, this simulation saves much calculation time. The test of several relative typical geoelectric models shows that this modeling method is correct and feasible.

关 键 词：有限单元法 三维地电断面 激发极化法 正演 三维构造 电位 视电阻率 偏导数矩阵 

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