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

黄俊革; 阮百尧; 鲍光淑

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

黄俊革^{1, 2,},
阮百尧²,
鲍光淑¹

1.
中南大学信息物理学院, 湖南长沙 410083
2.
桂林工学院资源与环境工程系, 广西桂林 541004

基金项目:

国家自然科学基金 40064001

详细信息

作者简介:
黄俊革(1966-), 男, 中南大学应用地球物理专业博士研究生. E-mail: hjunge@163.com

中图分类号: P624
计量
- 文章访问数: 3621
- HTML全文浏览量: 426
- PDF下载量: 11
- 被引次数: 0
出版历程
- 收稿日期: 2002-06-08
- 刊出日期: 2003-05-25

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

1.
Institute of Information and Physics, Central South University, Changsha 410083, China
2.
Department of Resources and Environment Engineering, Guilin Institute of Technology, Guilin 541004, China

摘要

摘要: 用有限单元法求解三维地电断面激发极化法正演模拟算法.首先给出了三维构造中点源电场异常电位的边值问题与变分问题, 简化无穷远边界上的边界条件以提高计算速度及计算精度.以此为基础计算视电阻率对模型电阻率的偏导数矩阵, 并进行三维地电断面激发极化法正演模拟, 与等效电阻率法相比节约了计算时间.对几例较典型的地电模型进行计算, 结果表明本方法是正确可行的.
- 三维地电断面 /
- 激发极化法 /
- 有限单元法 /
- 偏导数矩阵
Abstract: 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.
- 3-D geoelectric section /
- IP /
- finite element method /
- partial derivatives

HTML全文

图 1 极化率正演(a₁, b₁, c₁) 与等效电阻率法计算结果(a₂, b₂, c₂) 对比(等值线单位为%)

Fig. 1. Comparison of IP modeling and equivalent resistivity method

下载: 全尺寸图片幻灯片

参考文献(7)

[1]	Seigel H O. Mathematical formulation and type curves for induced polarization[J]. Geophysics, 1959, 24: 547- 565. doi: 10.1190/1.1438625
[2]	Tripp A C, Hohmann G W, Swift C M, et al. Two-dimensional resistivity inversion[J]. Geophysics, 1984, 49: 1708- 1717. doi: 10.1190/1.1441578
[3]	阮百尧. 视电阻率对模型电阻率的偏导数矩阵计算方法[J]. 地质与勘探, 2001, 37: 39- 41. https://www.cnki.com.cn/Article/CJFDTOTAL-DZKT200106009.htm RUAN B Y. A generation method of the partial derivatives of the apparent resistivity with respect to the model resistivity parameter[J]. Geology and Prospecting, 2001, 37: 39- 41. https://www.cnki.com.cn/Article/CJFDTOTAL-DZKT200106009.htm
[4]	徐世浙. 地球物理中的有限单元法[M]. 北京: 科学出版社, 1994. 178- 188. XU S Z. Finite element method in geoghysics[M]. Beijing: Science Press, 1994. 178- 188.
[5]	阮百尧, 熊彬, 徐世浙. 三维地电断面电阻率测深有限元数值模拟[J]. 地球科学——中国地质大学学报, 2001, 26(1): 73- 77. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200101013.htm RUAN B Y, XIONG B, XU S Z. Finite element method for modeling resistivity sounding on 3-D geoelectric section [J]. Earth Science— Journal of China University of Geosciences, 2001, 26(1): 73- 77. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX200101013.htm
[6]	黄俊革, 阮百尧, 鲍光淑. 齐次边界条件下三维地电断面电阻率有限元数值模拟法[J]. 桂林工学院学报, 2002, 22: 11- 14. https://www.cnki.com.cn/Article/CJFDTOTAL-GLGX200201002.htm HUANG J G, RUAN B Y, BAO G S. FEM under quantic-boundary condition for modeling resistivity on 3-D geoelectric section[J]. Journal of Guilin Institute of Technology, 2002, 22: 11- 14. https://www.cnki.com.cn/Article/CJFDTOTAL-GLGX200201002.htm
[7]	阮百尧, 村上裕, 徐世浙. 电阻率/激发极化法数据的二维反演程序[J]. 物探化探计算技术, 1999, 21: 116- 125. RUAN B Y, Yutaka Murakami, XU S Z. 2-D inversion programs of induced polarization data[J]. Geophysical and Geochemical Exploration, 1999, 21: 116- 125.