共享知识 分享快乐罗延钟1，万乐1.中国地质大学数学系，武汉4300742.中国科学院高能物理研究所，北京100049：介绍在地质大学电法科研组“高密度电阻率法正、反演软件”的基础上，研制高密度电2.5维电阻率成像系统；讨论不平地形条件下高密度电阻率法2.5维电阻率成像的算法。
：高密度；电阻率法；成像0 前言目前，国内外高密度电阻率法地下电阻率成像主要应用Luke M.研制的软件系统[1－2]。
生成的图像甚至和导电体的倾斜产状相反[3]。
本电法科研组研制2.5维电阻率成像系统克服了前人系统的缺陷，在应用中取得了较好的2005年4月优选出中国地质大学电法科研组的“高密度电阻率法正、反演软件”设在高密度电阻率法观测剖面上。
反演基础的正演计算假设地面起伏和地电构造都是1 算法概述高密度电阻率法在预先选定的地面观测剖面和测点上，同时布置几十乃至上百个电极，2.5维。
2.5维。
反演线性方程组的建立59637050）；“十五”国家科技攻关项目（2004BA523B ）附件二维和具有相同走向的，且观测剖面垂直于构造走向。
成像任务是由观测剖面上高密度电阻率法的视电阻率观测数据，反演该剖面下地电断面[3，6]。
地下电场成三维分布，故称所论地下电阻率成像是2.5维。
因为假定地电结构是二维的；而电极装置使用电极故称所论地下电阻率成像。
1.1.1 原理和方法为检验成像系统的实用性，用该系统对浙江某地探测岩溶和秦始皇陵考古探测的高密度电阻率法实测数据作了成像。
用Surfer 软件绘制的成像结果示于图3和图4（略）。
图3显著地指示出地下低电阻率岩溶的位置。
而秦皇陵的成像结果与洛阳铲钻探结果吻合得很高密度电阻率法在预先选定的地面观测剖面和测点上，将这些电极组合成指定的电极装置和电极距[4－5]。
假设地面起伏和地电构造都是二维和具有相同走向的，且观测剖面垂直于构造走向。
成像任务是由观测剖面上高密度电阻率法的视电阻率观测数据，反演该剖面下地电断面，绘制地下岩石的真电阻率分布图像。
⎥⎥⎦⎤⎢⎢⎣⎡==∑j jp ij ip ip Ιt w f t x f y )()]([ （1）)()()()(t w t w t w t E ηij ij pij p Λ+=∂∂∑（2）同时布置几十乃至上百个电极，然后用多芯电缆将它们连接到特制的电极转换装置。
成像任务是由观测剖面上高密度电阻率法的视电阻率观测数据，反演该剖面下地电断面，绘制地下岩石的真电阻率分布图像[6，8，10－14]。
假设地面起伏和地电构造都是二维和具有相同走向的，且观测剖面垂直于构造走向。
成像任务是由观测剖面上高密度电阻率法的视电阻率观测数据，反演该剖面下地电断面，绘制地下岩石的真电阻率分布图像。
成像任务是由观测剖面上高密度电阻率法的视电阻率观测数据，反演该剖面下地电断面，绘制地下岩石的真电阻率分布图像。
假设地面起伏和地电构造都是二维和具有相同走向的，且观测剖面垂直于构造走向。
成像任务是由观测剖面上高密度电阻率法的视电阻率观测数据，反演该剖面下地电断面，绘制地下岩石的真电阻率分布图像。
成像任务是由观测剖面上高密度电阻率法的视电阻率观测数据，反演该剖面下地电断面，绘制地下岩石的真电阻率分布图像。
假设地面起伏和地电构造都是二维和具有相同走向的，且观测剖面垂直于构造走向。
成像任务是由观测剖面上高密度电阻率法的视电阻率观测数据，反演该剖面下地电断面，绘制地下岩石的真电阻率分布图像。
图1 典型BP网络示意图高密度电阻率法在预先选定的地面观测剖面和测点上，同时布置几十乃至上百个假设地面起伏和地电构造都是二维和具有相同走向的，且观测剖面垂直于构造走向。
成像任务是由观测剖面上高密度电阻率法的视电阻率观测数据，绘制地下岩石的真电阻率分布图像。
高密度电阻率法在预先选定的地面观测剖面和测点上，同时布置几十乃至上百个电极，然后用多芯电缆将它们连接到特制的电极转换装置；可根据操作员的指令，将这些电极组合成指定的电极装置和电极距。
假设地面起伏和地电构造都是二维和具有相同走向的，且观测剖面垂直于构造走向。
成像任务是由观测剖面上高密度电阻率法的视电阻率观测数据，反演该剖面下地电断面，绘制地下岩石的真电阻率分布图像。
4 结束语本文介绍了中国地质大学电法科研组研制的不平地形条件下高密度电阻率法2.5维成像系统，该系统可用于九种电极装置观测数据的成像；系统的输入文件简单；系统采用的Quick-Win界面美观，用户友好性强，运行操作简便；成像结果输出形式多样。
参考文献[1] 成谷, 马在田, 耿建华,等. 地震层析成像发展回顾[J]. 勘探地球物理进展, 2002, 25(3): 4-7.[3] Aki K, Lee W H K. Determination of three-dimensional velocity anomalies under a seismic arrayusing first P arrival times from local earthquakes[J]. Geophys Res, 1976, 81(23): 4381-4399.structure of the lithosphere Geophys[J]. Journal of Geophysical Research, 1977, 82: 277-96.and a reexamination of the correction of ellipticity Geophys[J]. J R Astron Soc, 1976, 44: 7-17.李娟. 首都圈地区Pn和PmP波层析成像研究[D].中国地震局地球物理研究所博士学位论文, 2003.[7] 杨晓春, 李小凡, 张美根. 地震波反演方法研究的某些进展及其数学基础[J]. 地球物理学进展,2001, 16(4): 96-109.[8] 吴振利, 李家虝, 阮爱国. 地震层析成像及其国内研究进展[J]. 东海海洋, 2003, 21(3): 54-64.[9] Zhao Dapen. New advance of seismic tomography and its application to subduction zones andearthquake fault zones[J]: A review, The Island Arc. 2001, 10: 68-84.共享知识 分享快乐[10] 郭飚. 流动宽频带地震台阵: 地震定位与地层层析成像[D]. 中国地震局地质研究所硕士论文,2003.[11] Thurber C, Aki K. Three-dimensional seismic imaging[J]. Annual Review of Earth and PlanetarySciences, 1987, 15: 115-139.[12] Tarantola A, Valette B. Inverse problems = quest for information[J]. Geophys, 1982, 50:159-170.[13] Tarantola A, Nercessian A. 1984. Three-dimensional inversion without blocks[J]. Geophys JR. astr. Soc. 1982, 76: 299-306.[16] 王夫运, 张先康, 陈棋福, 等. 北京地区上地壳三维细结构层析成像[J]. 地球物理学报, 2005,48(2): 359-366.[15] Burmakov, J A, Treussov A V, Vinnik L P. Determination of three dimensional velocity structurefrom observations of refracted body waves[J]. Geophys J R astr Soc, 1984, 79: 285-292.[16] 朱露培, 曾融生, 刘福田. 一种新的三维速度结构反演模型参数化方法[J]. 地球物理学报, 1990,33(1): 34-43.[17] 罗延钟, 万乐, 董浩斌,等. 高密度电阻率法的2.5维反演[J]. 地质与勘探, 2003, 39(1): 107-113。
[18] 熊彬，张桂青. 电子计算机在电法勘探中的应用[M]. 武汉: 武汉地质学院出版社，1987. [19] Tripp A C, Hohmann G W, Swift C. G. Two-dimensional resistivity inversion: Geophysics[J].1984, 49: 1708-1717. [20] Jupp, D, Jonnie L B, Vozoff K. Two-dimensional magnetotelluric inversion: Geophys[J]. J Roy. Astr, 1976, 50: 333-352. [21] Qin F, Olsen K B, Cai W et al. Finite-difference solution of the eikonal equation alongexpanding wavefronts[J]. Geophysics, 1992, 57(3): 478-487.[22] Benxi Ke, Bo Zhao, Jiaming Cai et al. 2-D finite element acoustic wave modeling includingrugged tomography[A]. Socv Expl Geophys 71th Ann. Internat. Mtg[C]. Tulsa, USA. 2001. [23] Kirkpatrck S, Gelatt C, Jr D, et al. Optimazation by simulated annealing[J]. Science, 1983,220: 671-680.[24] 成谷, 马在田, 耿建华,等. 地震层析成像发展回顾[J]. 勘探地球物理进展, 2002, 25(3).High Density Resistivity ImageLuo Yan-zhong ，(China University of Geosciences, Wuhan 430074, China )Abstract ：The algorithms and software system of 2.5-D resistively image for High DensityResistively in the condition of rolling terrain are discussed. This software system can be used fornine common pole-arrays. The image results for modeling and measured data indicate that the algorithms are exactness, the program is stabilization and image effects are fine.Key words ：High Density Resistively ；2.5 dimension ；Resistively image盛年不重来，一日难再晨。