非线性地球物理反演方法：回顾与展望-评地球物理反演的发.pdfVIP

地学前缘 Vol.10, 2003 评地球物理反演的发展趋向* 杨文采 中国地质科学院地质研究所（北京，100037） [摘要] 自20 世纪 70 年代以来，地球物理反演方法已走向成熟，这些方法包括地震走时反演与 波场反演、层析成像反演、线性化迭代反演、仿真淬火法反演、遗传算法及联合反演方法等。地 震反演的目标不仅仅是波速模型, 而是地球内部的精细组构和动力学过程, 并向密度、各向异性 及粘度等更多参数成像发展。波动方程线性化迭代反演基于开放物理系统状态发生相变的原理， 要进一步改善模型参数化的技术和迭代过程的自组织；仿真淬火法与遗传算法基于自然过程的指 数率或生物演化的优生率，可以相互结合以提高解估计的分辨率与置信度; 联合反演要结合岩石 物理性质的统计规律，才能取得兼容地质与综合方法的应用效果。地球动力学中的反问题不仅涉 及偏微分方程系数项的求解，而且涉及初始条件或初始边界的求解，对地球动力学作用过程研究 特别重要。 [关键词] ：地球物理，反演，非线性，联合反演，地球精细结构 三十多年前，在板块构造学说与计算机革命的双重推动下，基于现代数学物理理 论的现代地球物理反演方法诞生了[1-2] 。在21 世纪到来之际，全球地球动力学发展和 能源资源勘探对地球物理反演提出了更高的要求，信息产业的迅速发展又为地球物 理反演方法的发展提供了更好的条件。作为这个领域中的前沿，自 20 世纪 90 年代 [3-15] 。本文在回顾线性化迭代反演 以来非线性地球物理反演已经取得一些可喜的进展 法、仿真淬火法、遗传算法与联合反演方法的基础上，从历史沿革、自然哲学与偏 微分方程反问题的角度展望了今后非线性地球物理反演方法研究的一些方向。 1 地震走时反演与层析成像 地球物理反演大约是从 1907 年发表地震走时反演论文开始的, 这个领域经百年 *本课题由国家重大科学工程项目“中国大陆科学深钻工程”与国家自然科学基金(共同资助 1 风尘, 取得了丰硕成果。目前用于地球广角/ 折反射资料反演波速结构的软件,如 RAYINV 和 SEIS88 已经非常流行, 以至于非反演专业人员也可以应用这些软件反 演地球波速结构。然而, 不同的人用这些软件反演出的地球波速结构有很大的差异, 而且用这些软件不可能反演出复杂的地球波速结构。图 1 给出反演地球波速结构的 图 1. 由地震走时反演取得地球波速结构的比较。(a) 80 年代由广角折射走时反演出的地球波 速结构; (b) 过温哥华岛的地壳纵波波速模型及(c)横波速度/ 电阻率模型。 Fig.1 Comparison of reconstructed velocity models obtained via traveltime inversions: (a) old GGT produced from wide-angle refraction data; (b) P-wave velocity profile and (c) S-wave velocity and resistivity profile across the Vancouver island (After Ron Clowes et al [16]) 典型剖面。图 1(a) 是 80 年代我们常见的由广角折射走时反演出的地球波速结构剖面, 基于层状介质模型, 这种反演包含的地球信息不丰富。图 1(b) 为过温哥华岛的地壳 纵波波速模型,它基于纵横向交叉分块的波速模型, 给出了太平洋板块俯冲的详细信 息 [16] 。图 1(c) 为过温哥华岛的地壳横波波速与电阻率联合反演模型,它在P波速度结 构的基础上细化了波速模型, 给出了太平洋板块俯冲及科迪勒拉造山带内幕更为详 细的信息。由此可见, 为取得更详细的地球信息我们不能依赖于现成软件, 而要发展 基于复杂模型的反演方法。为此首先要发展基于复杂模型的正演方法,克服复杂模型