勘探地球物理层析成像软件系统及其应用-物探与化探.PDF

第 22 卷第 2 期 物 　探 　与 　化 　探 Vol . 22 ,No . 2 1998 年 4 月 GEOPH YSICAL GEOCHEM ICAL EXPL ORA TION Apr . ,1998 勘探地球物理层析成像软件系统及其应用 朱介寿 　严忠琼 　曹俊兴 　张雪梅 (成都理工学院 ,成都 　6 10059) 摘 　要 　本文对人工场源激发的地震波 、声波 、电磁波及电阻率分布而重建地下物理参数图像的 理论 、正反演算法 、计算机软件 、实验室超声波模型实验 、电阻率实验等做了大量的研究 ,完成了由 Windows95 支持的勘探地球物理层析成像软件系统的研制 ,并介绍了三个勘探应用实例 。 关键词 　井间地球物理层析成像 ,软件系统 ,勘探应用 。 0 　引言 近年来 , 由于专门用于地球物理层析成像的观测仪器及计算技术的进步 ,勘探地球物理层 析成像技术及应用范围已取得了很大的进展 。首先 ,在井间地震波层析技术方面 , 已经有多种 激发源 ,利用三分量井中检波器或压电灵敏元件做接收器 。美国、加拿大等国在井间采用可控 压电式震源及三分量检波器 ,观测 1 000～5 000 Hz 的高频地震波 ,采用多次叠加及对地震波 进行分离处理 ,地下层间分辨率已提高到 1 m 。我国研制的第四代电磁波透射仪 ,可发射 0 . 2 ～32 M Hz 的强信号 ,使井间电磁波层析成像效果有很大提高 。由日本发展的井中及地面高 密度电阻率剖面观测系统 ,为井间、坑道间或地井 、地面间的电阻率层析成像开拓了良好的应 用前景 。 1 　井间地球物理层析成像理论与算法 1. 1 　井间地震走时层析成像 ( ) 基于惠更斯 Huygens 原理与网络理论的最小走时最短路径射线追踪为基础的走时层析 成像正演理论及算法 ,能模拟任意复杂介质射线 ,保证阴影区亦有射线通过 。该方法计算速度 快 、收敛稳定 、分辨率高 ,是 目前最先进的射线追踪算法 。发展了层析成像的反演理论及多种 快速稳定的重建算法 ,对迭代求解算法的稳定性 、收敛速度 、分辨能力进行了比较和研究 。 无论介质有多复杂 , 由点源激发的地震波的传播完全可以用 Huygens 原理做可靠的描述 。 Huygens 原理的射线追踪有两种实现方法 ,一是基于网络理论的最短路径算法 〔11〕,一是基于 动力学的波阵面算法 。这两种算法都能模拟直达波 、折射波 、反射波 、散射及绕射波 ,而且一次 计算即可得到一个共激发点记录的全部走时 ,计算效率很高 。当网格剖分密度适当时 ,所达到 的精度也相当高 。其中建立在网络理论基础上而寻求最短路径的方法 ,是 目前追踪不均匀介 质中真实射线较好的方法 ,很适合于层析成像问题的大量高精度射线追踪计算 。本软件系统 采用基于 Huygens 原理的网络追踪 ———最短路径射线追踪算法作为正反演计算方法 。 ( ) 最短路径算法的基本思路如下 :设 t i 为激发点 S 到第 i 个接收点的最小走时 1997 年 8 月 5 日收稿, 同年 10 月 15 日收修改稿 。 1 期 朱介寿等 : 勘探地球物理层析成像软件系统及其应用 ·9 1 · ( ) ( ) ( ) t i = min 〔t j + dij 〕　　　i , j ∈N 1 j ≠i 其中 N 为包括任意 n 个结点及弧段的集, 而 j 表示 i