评地球物理反演的发展趋向_杨文采

评地球物理反演的发展趋向

杨文采

(中国地质科学院 地质研究所,北京100037)

摘 要:自20世纪70年代以来,地球物理反演方法已走向成熟,这些方法包括地震走时反演与波场反演、层析成像反演、线性化迭代反演、仿真淬火法反演、遗传算法及联合反演方法等。地震反演的目标不仅仅是波速模型,而且是地球内部的精细组构和动力学过程,并向密度、各向异性及粘度等更多参数成像发展。波动方程线性化迭代反演基于开放物理系统状态发生相变的原理,要进一步改善模型参数化的技术和迭代过程的自组织;仿真淬火法与遗传算法基于自然过程的指数率或生物演化的优生率,可以相互结合以提高解估计的分辨率与置信度;联合反演要结合岩石物理性质的统计规律,才能取得兼容地质与综合方法的应用效果。地球动力学中的反问题不仅涉及偏微分方程系数项的求解,而且涉及初始条件或初始边界的求解,对地球动力学作用过程研究特别重要。关键词:地球物理;反演;非线性;联合反演;地球精细结构

中图分类号:P3 1 文献标识码:A 文章编号:10052321(2002)04038908

收稿日期:2002

03

04;修订日期:2002

0521

基金项目:国家重大科学工程资助项目 中国大陆科学深钻工程 ;国家自然科学基金资助项目(40174013)

作者简介:杨文采(1942 ),男,博士,研究员,博士生导师,地球物理学专业。

30多年前,在板块构造学说与计算机革命的双重推动下,基于现代数学物理理论的现代地球物理反演方法诞生了

[1,2]

。在21世纪到来之际,全球地

球动力学发展和能源资源勘探对地球物理反演提出了更高的要求,信息产业的迅速发展又为地球物理反演方法的发展提供了更好的条件。作为这个领域中的前沿,自20世纪90年代以来非线性地球物理反演已经取得一些可喜的进展[3~15]

。本文在回顾线性化迭代反演法、仿真淬火法、遗传算法与联合反演方法的基础上,从历史沿革、自然哲学与偏微分方程反问题的角度展望了今后非线性地球物理反演方法研究的一些方向。

1 地震走时反演与层析成像

地球物理反演大约是从1907年发表地震走时反演论文开始的,这个领域历经百年风尘,取得了丰硕成果。目前用于地球广角折/反射资料反演波速结构的软件,如RAYINV 和SEIS88已经非常流行,以至于非反演专业人员也可以应用这些软件反演地球波速结构。然而,不同的人用这些软件反演出的地球波速结构有很大的差异,而且用这些软件

不可能反演出复杂的地球波速结构。图1给出反演地球波速结构的典型剖面。图1a 是20世纪80年代我们常见的由广角折射走时反演出的地球波速结构剖面,基于层状介质模型,这种反演包含的地球信息不丰富。图1b 为过温哥华岛的地壳纵波波速模型,它基于纵横向交叉分块的波速模型,给出了太平洋板块俯冲的详细信息

[16]

。图1c 为过温哥华

岛的地壳横波波速与电阻率联合反演模型,它在P 波速度结构的基础上细化了波速模型,给出了太平洋板块俯冲及科迪勒拉造山带内幕更为详细的信息。由此可见,为取得更详细的地球信息我们不能依赖于现成软件,而要发展基于复杂模型的反演方法。为此首先要发展基于复杂模型的正演方法,克服复杂模型引起的射线焦散及波场奇点等困难。

以走时反演为基础的地震层析成像自80年代以来取得了巨大成功,尤其是提供了全地幔P 波、横波与面波波速结构图像,扩大了地球科学家的眼界,这些成果已为大家熟知[17]。近年来,地球物理界已不满足于波速成像,开始向密度、各向异性及粘度等更多参数成像发展,促进人们对深部岩性及流变性质的了解和相关地球动力学问题的解决[18~

20]

。图2是过格鲁吉亚西天山造山带岩石圈

层析成像的一个例子。通过层析成像反演计算出P 波波速剖面(图2a),横波波速剖面(图2b),密度反演剖面(图2c),通过综合分析取得了比较可靠的推测岩性剖面(图2d)。由此例可见,当前地球物理反

第9卷第4期2002年10月

地学前缘(中国地质大学,北京)

Earth Scien ce Frontiers (China University of Geosciences,Beij ing)Vol.9No.4Oct.2002

演的目标不仅仅是波速模型,而是地球内部的精细

组构和动力学过程,

这方面以后还要讨论。

图1 由地震走时反演取得地球波速结构的比较。(a)80年代由广角折射走时反演出的地球波速结构;

(b)过温哥华岛的地壳纵波波速模型及(c)横波速度/电阻率模型

[16]

Fi g 1 Comparison of reconstructed velocity models obtained via traveltime inversions:(a)old GGT produced from w ide angle refraction data;(b)P w ave velocity profi le and (c)S wave velocity and resistivity profile across the Vancouver island

[

16]

2 过格鲁吉亚西天山造山带岩石圈的P 波波速剖面(a),横波波速剖面(b),密度反演剖面(c)与推测岩性剖面(d)

[20]

Fi g.2 A profile acros s w est T i an S han orogene:(a)P wave velocity,(b)S w ave veloci ty,(c)density sections and

(d)inferred petrology section based on sections (a-c)

[20]

岩性编号:1 酸性岩;2 中性岩;3 基性岩;4 超基性岩;5 韧性流变层

390

杨文采/地学前缘(Ear th Science F ront iers)2002,9(4)