地球科学与环境学报
JOURNAL OF EARTH SCIENCES AND ENVIRONMENT
2005 Vol.27 No.3 P.63-67


基于序列二次规划方法的高密度电阻率反演

伍开江  周启友 

摘 要:针对高密度电阻率成像法反演计算中存在的困难,探讨了基于序列二次规划(SQP)方法建立可行的电阻率反演算法的问题.在对三维点电源二维地电体电位场模型有限元法正演计算的基础上,建立了基于SQP方法的高密度电阻率成像法反演算法.根据勘测中测量电极间距的不同,提出按不同的电极间距分别建立优化模型进行参数优化.通过对模拟电阻率模型和实际观测资料的反演计算,表明该算法进行高密度电阻率反演是可行的,其具有对初始模型无特殊要求、收敛速度快的特点.这一反演算法可望在高密度电阻率成像的反演计算中得到应用.
关键词:电阻率;反演;有限元法;参数优化;SQP方法
分类号:P631.3+22  文献标识码:A

文章编号:1672-6561(2005)03-0063-05

High density electrical resistivity inversion based on sequential quadratic programming method

WU Kai-jiang  ZHOU Qi-you 

基金项目:国家自然科学基金项目(40371021)
作者简介:伍开江(1979-),男,四川双流人,硕士研究生,从事高密度电阻率成像法的应用研究.
作者单位:伍开江(南京大学,地球科学系,江苏,南京,210093) 
     周启友(南京大学,地球科学系,江苏,南京,210093) 

参考文献:

[1]Daily W, Ramirez A, LaBrecque D, et al. Electrical resistivity tomography of vadose zone water movement[J]. Water Resour Res, 1992,28(5): 1424~1442.
[2]Ramirer A, Daily W, Labrecque D, et al. Monitoring an underground steam injection process using electrical resistance tomography[J]. Water Resour Res, 1993,29(1): 73~87.
[3]Binley A, Henry-Poulter S, Shaw B. Examination of solute transport in an undisturbed soil column using electrical resistance tomography[J]. Water Resour Res, 1996,32 (4):763~769.
[4]Guo Y Q, Batoz J L, Naceur H, et al. Recent developments on the analysis and optimum design of sheet metal forming parts using a simplified inverse approach[J]. Computers and Structures, 2000,78: 133~148.
[5]Kyoung K S, Lee W I. Solution of inverse heat conduction problems using maximum entropy method[J]. International Journal of Heat and Mass Transfer Volume, 2002,45: 381~391.
[6]Villafafila A, Mujtaba I M. Fresh water by reverse osmosis based desalination: simulation and optimization[J]. Desalination, 2003,155: 1~13.
[7]Tripp A C, Hohmann G W, Swift C M Jr. Two-dimensional resistivity inversion[J]. Geophysics, 1984,49(10): 1708~1717.
[8]Sasaki Y. Two-dimensional joint inversion of magnetotelluric and dipole-dipole resistivity data[J]. Geophysics, 1989,54(2): 254~262.
[9]Sasaki Y. 3-D resistivity inversion using the finite-element method[J]. Geophysics, 1994,59 (11): 1839~1848.
[10]Shima H. Two-dimensional automatic resistivity inversion technique using alpha centers[J]. Geophysics, 1990,55(6):682~694.
[11]Zhang J, Mackie R L,Madden T R. 3-D resistivity forward modeling and inversion using conjugate gradients[J]. Geophysics, 1995,60(5): 1313~1325.
[12]Zhou Q Y, Shimada J, Sato A. Three-dimensional soil resistivity inversion using patching method[J]. Jour Japan Soc Eng Geol, 1999,39(6): 524~532.
[13]杨冰.实用最优化方法及计算机程序[M].哈尔滨:哈尔滨船舶工程学院出版社,1994.
[14]袁亚湘,孙文瑜.最优化理论与方法[M].北京:科学出版社,1997.

收稿日期:2004年9月26日

出版日期:2005年9月1日