SOME PROBLEMS CONCERNING 2D INVERSION OF DIRECT CURRENT IP SOUNDING DATA
LIU Hai-fei1, RUAN Bai-yao2, LÜ Yu-zeng1,2
1. School of Info-Physics and Geomatics Engineering Central South University, Changsha 410083, China;
2. Department of Resources and Environment Engineering Guilin Institute of Technology, Guilin 541004, China
At present, the researches on 2D inversion of IP data mainly aim at the HD (high density) IP data, which is mainly due to the equivalent distances between neighboring poles during data collection. Under such a condition, the division of finite element grid is convenient and the forward and inverse operation is easy. In forward operation, the finite element grid is automatically divided according to the practical data, and dual-grid system is used. The finite element modeling with linear variation of conductivity within each block is done, and the error of modeling result can be modified. In inversion, the partial derivatives of the apparent resistivity with respect to model resistivity are calculated by means of the relationship between the source and receiver locations and Broyden's method. And it is solved by using the conjugate gradient method, in which the calculation precision is almost the same as that in the SVD method. Through inversing the modeling data and observed data, it is shown that the algorithms of this paper are workable.
刘海飞, 阮百尧, 吕玉增. 直流激电测深二维反演的若干问题研究[J]. 物探与化探, 2007, 31(1): 47-50,54.
LIU Hai-fei, RUAN Bai-yao, LÜ Yu-zeng. SOME PROBLEMS CONCERNING 2D INVERSION OF DIRECT CURRENT IP SOUNDING DATA. Geophysical and Geochemical Exploration, 2007, 31(1): 47-50,54.
[1] Tripp A C,Hohmann G W,Swift C M.Two-dimensional resistivity inversion[J].Geophysics,1984,49:1708.
[2] Yutaka Sasaki.3D Resistivity inversion using the finite-element method[J]. Geophysics,1994,59:1839.
[3] 罗延钟.电子计算机在电法勘探中的应用[M].武汉:武汉地质学院出版社,1987.
[4] 阮百尧,村上裕,徐世浙.电阻率/激发极化率数据的二维反演程序[J].物探化探计算技术,1999,21(5).
[5] 阮百尧.三角单元剖分电导率分块连续变化点源二维电场有限元数值模拟[J].广西科学,2001,8(2).
[6] 吴小平,汪彤彤.电阻率三维复杂结构的快速反演[J].煤田地质与勘探,2001,29(2).
[7] 阮百尧,徐世浙.二维直流电阻率测深曲线的快速反演[J].物探与化探,1996,20(6).
[8] Torres-Verdin,Druskin C,Fang V.A dual-grid nonlear inversion technique with applications to the interpretation of dc resistivity data[J].Geophysics,2000,65:1733.
[9] 阮百尧.视电阻率对模型电阻率的偏导数矩阵计算方法[J].地质与勘探,2001,37(6).
[10] Loke M H,Barker R D.Rapid least-squares inversion of apparent resistivity pseudosections by a quasi-Newton method[J].Geophysical prospecting,1996,44:131.
[11] 李家康,姜兰茵,刘宇.共轭梯度法和奇异值分解法与降元迭代法的对比讨论[J].地球物理学报,1994,37,(Ⅰ).
[12] 周竹生,赵荷晴.广义共轭梯度算法[J].物探与化探,1996,20(5).
[13] 吴小平,汪彤彤.利用共轭梯度方法的电阻率三维反演若干问题研究[J].地震地质,2001,23(2).
[14] 周竹生,何继善,赵荷晴.利用广义共轭梯度算法求解地震道反演问题[J].石油地球物理勘探,1998,33(4).