吕玉增-电阻率极化率2.5维正反演及地-井方位IP三维人机交互解释1

电阻率/极化率2.5维正反演地-井方位IP三维人机交互解释吕玉增桂林理工大学中国矿业大学（北京）2010.12

直流激电正反演研究现状应用中的热点问题直流激电2.5维正反演软件介绍地井方位激电数据的人机交互解释软件介绍汇报提纲一直流激电正反演研究现状国外上世纪70、80年代基本完成了相关的研究已经完成一维、二维和三维正反演系统研究2.5维、三维反演基本成熟，已达到应用阶段复杂地形等条件下三维正反演及其应用井中、坑道等不同环境工作方法、解释和应用坑道超前探测等方面研究

研究现状

目前研究热点二应用中的热点问题

中梯数据的反演问题常规对称四极测深数据反演问题直流激电深部找矿问题观测（测线、测点）问题等中梯数据的反演（1）激电中梯扫面数据3DIPPs改正效率中梯数据的反演（1）激电中梯（不对称四极）剖面数据2D（2.5D）问题可以是任意组合：

A1-B1，A1-B2，A2-B1，A2-B2，……

＋观测电极的任意组合数据量大，信息量更丰富常规对称四极测深数据反演问题（2）1.极距的变化（高密度－均匀间隔）（对称四极－近似对数等间隔）常规对称四极测深数据反演问题（2）2.网格的剖分依照高密度数据等网格剖分反演：一个测点，观测10个极距的对称四极2.5维反演相当于40个电极的高密度数据反演，如果有10几，几十个测点的对称四极数据反演计算时间太长，不实用；常规对称四极测深数据反演问题（2）2.网格的剖分（不对称网格剖分+插值技术）以10个对称四极测点，每个测点10个极距计，等网格剖分水平方向需要10×4×10＝400个用以上方法剖分，仅需要20+10=30个激电深部勘探问题（3）

1.地面激电深度勘探？（直流激电是几何测深，地面物探工作要想勘探深度大势必增大供电极距。）

2.开展井中、坑道激电等地下物探勘查应该是深部勘查的一个方向地井井地井井

观测（测线、测点）问题（4）测线二维构造

1.测线问题

观测（测线、测点）问题（4）

2.测点问题不够灵活－仪器设备、软件解释等限制

不等间隔网格剖分技术变不规则剖分为规则剖分

不等间隔剖分技术＋插值技术

2D

3D

任意3D不规则网观测和解释三直流激电2.5维正反演软件介绍

正反演软件是在桂林理工大学阮百尧教授开发的电阻率/极化率数据2.5维正反演程序的基础上研制开发的。2.5D正演软件适合于目前高密度电法和常规对称四极等不同装置下、任意复杂地形和2.5维地电体的正演模拟。2.5D反演软件适用于目前高密度电法、常规对称四极等不同装置的直流激电、双频仪等的实测数据（视电阻率和视极化率）反演解释。也适用于水平地形和起伏地形下的数据解释，且不需要对地形影响进行校正。正演算法采用有线单元法，三角单元剖分、电位和电导率连续变化等先进技术，正演精度高。反演采用了约束最小二乘最优化反演算法，加入了最光滑、最简单模型、背景电性参数等约束条件，舍得反演效果能更接近实际情况。三直流激电2.5维正反演软件介绍

1.2.5D正演算法异常电位的边值问题付氏变换——将点源二维变为波数域中的二维如果地下介质的电性分布有走向，取走向方向为z,对电位在z方向进行傅立叶变换对若干个k，解出波数域中点电源总电位V，用傅氏反变换，可以求出主剖面上的三维空间中的v。傅氏反变换的公式是：异常电位的边值问题波数域中异常电位的微分方程异常电位的边值问题傅氏反变换的kI和gi

（徐世浙，1998）

12345kI0.002680.024190.096210.339911.26051gi0.005650.024060.078490.274091.09590

异常电位的变分问题波数域中异常电位的微分方程异常电位的变分问题微分方程转变为变分问题

有限元法解变分问题解线性方程组得节点上波数域异常电位异常电位的变分问题波数域中异常电位的微分方程付氏逆变换加上正常电位，可得总电位v。按观测方式计算视电阻率。

由上述方法计算不同波数时各节点上的波数域异常电位值，然后用下式付氏反变换得各节点空间域异常电位。2.5D正演模拟软件介绍2.5D正演