高密度电阻率法的维反演.pdfVIP

第39卷增刊 地质与勘探 V01．39 2003年10月 GEOLOGYANDPROSPECTING October，2003 高密度电阻率法的2．5维反演 罗延钟，万 乐，董浩斌，熊 彬，王传雷 (中国地质大学，武汉430047) [摘要]讨论了电阻率法2．5维正演和反演的算法，并在此基础上编制了高密度电阻率法2．5维 反演程序，该程序可用于8种常用电极装置观测结果的反演。对理论和实测数据的反演结果表明，采用 的算法正确，程序运行稳定，反演效果很好。 [关键词]高密度电阻率法反演 [中图分类号]P631．3 2，…，Ⅳ)就是反演要确定的未知参数。 O 前言 1．1反演线性方程组的建立 高密度电阻率法在我国工程勘查中已获得相当 假设在实际工作中测到一组(共M个)实测视 广泛的应用，但目前国内的高密度电阻率法应用软 电阻率数据，称为实测数据矢量{n}， (1) 件，基本上只能做数据处理和成图等定性解释，尚无 {Ps}={风。，风：，…，风i，…，风肘}1 具有自主知识产权的定量解释商业软件。笔者在早 反演的目的是通过这组实测数据求出地下岩石 年研制“电法勘探工作站”的基础上，开发了适用于 多种电极装置的高密度电阻率法2．5维定量解释软 {仃}， 件，经试算效果甚佳。本文将介绍该软件及其试算 {仃}-{∥，，or：，…，q，…，orⅣ}7 (2) 情况。 某一电极距在某一测点的实测视电阻率风i是 l 算法原理 假设地电构造是二维的，观测剖面垂直于构造 可以给定；{6仃}为模型矢量的改正量，是未知的。 走向。反演的任务是由上述剖面的视电阻率观测数 只要求到了改正量，模型矢量也就求得。 据，定量推断该剖面下的地电断面——确定地下岩 将风。({or})作多变量台劳级数展开： 石的真电阻率分布。设在高密度电阻率法观测剖面 上布置的总电极数为N。，相邻电极间距为o，电极装 J21 oor二 置的最小和最大间隔系数为Ⅳ胁和Ⅳ。。。。将反演区 ‘呐+R (3) 沿横向剖分为步长为。的虬个节点，分别与观测剖 尺是高次项，将其略去，上式就是一个线性方 面上的各个电极位置相重合；沿纵向按同样步长。 程。其中，Pfi是与实测视电阻率风。相对应的、对 初始模型矢量{oro}，用有限单元法计算出的视电阻 剖分为(Ⅳ。。。+1)／2+1个节点。这样，反演区就被 剖分为N=(Nd一1)×(N…+1)／2个正方形单元 率理论值；塑掣是第i个视电阻率对第J个 组成的长方形网格。由其分别向两侧和下方延伸出 doi 15个步长依次按1．3倍增大的矩形单元，连同反演模型参数的一阶偏导数理论值。对于肘个实测视 区共同组成有限单元法正演计算的网格。 电阻率数据，我们有 将反演区内Ⅳ个正方形单元分别视为电性均 匀，它们的电导率Ori或电阻率Pi(Or，=1／pi，i=1， [收稿日期]2