2016_WPI_23_时间空间域全波形反演.docVIP

时间-空间域全波形反演 任浩然，王华忠 同济大学海洋与地球科学学院波现象与反演成像研究组，上海，200092 摘 要 地震波反演成像的理论研究地震波场在地球介质中的传播现象，从地震数据直接反求物性参数。而全波形反演是利用全波场进行地震反演的方法。在反演框架中，正演方法是反演成败的基础。基于二次型泛函，非线性寻优可以采用梯度导引类的解法，也可以采用牛顿类的解法。我们主要讨论梯度导引类方法求解时间-空间域全波形反演的方案。基于模型试验，讨论了全波形反演方法的优缺点，并就全波形反演的实用化提出了相应策略。 [关键词]：地震波反演，全波形反演，特征波场，非线性反演 1 引言 从Tarantola(1984)建立以波动方程为基础的地震反演方法以后的近三十年来，地震波反演方法已经建立了较为完善的理论体系。这些方法可以分别从不同的模型参数化、不同的正演方法和不同的反演策略来区分。其中，介质模型的参数化与正演方法相互对应。对模型有什么样的认识，才能建立什么样的正演方程,有什么样的正演方程才会有什么样的反演策略。 逆散射层析技术基于波场的线性近似，将地下的成像点看作是互不相干的散射点。传统上的散射层析只针对均匀背景模型（Devaney，1982；Wu和 Toks?z，1987）。后来发展的广义散射层析基于De Wolf近似(De Wolf，1985)，这种对散射波场进行多次前向散射和单次后向散射的近似使得广义散射层析在非均匀背景上进行反演(朱小三，2010)。基于旅行时的层析反演方法可以分为：利用初至波或折射波的近地表速度反演(刘玉柱，2007)，利用反射波的深层速度反演（Billette和Lambare，1998），利用井间透射波旅行时反演井旁精细速度反演（Harris等，1995）等。基于射线理论的各种反演方法本质上还是模拟波场的旅行时信息，因此可以反演影响旅行时的速度等参数。然而，勘探工业发展的现实表明，仅仅研究速度参数是不够的。后期的地震解释要求偏移成像是一个保振幅的，即按照偏移或者反演估计出的反射系数随角度的变化应当能够反应真实的地下物性。这就要求对模型的参数化要更加贴近物理现实、正问题过程要更加真实的描述地震波场的传播过程。另外，近年来随着计算机技术的发展，基于全波动方程的逆时偏移技术已经得到工业界大规模应用。因此，有理由相信，本世纪第二个十年将会是基于波动方程的反演方法大规模应用的阶段。 当前，基于声波全波方程的反演还没能得到很好的应用，制约的因素主要有：如何估计一个较为精确的初始地球物理模型；如何构造更好的泛函以减少各种因素造成的振幅误差敏感性；正则化的问题如何处理，重建地震数据的低频信息。当然，对实际地震数据的FWI，仍然存在诸多制约因素，比如低信噪比、低频缺失、初始模型偏差太大、各向异性、粘性衰减，这些因素都有可能造成寻优过程的不收敛。在FWI的实用化过程中，面临的主要问题可以归纳为： （1）建立一个较为精确的初始地球物理模型，体现对模型的某种假设（地质的或地球物理的）； （2）构建一个更好的泛函，以减少FWI对各种因素造成的振幅误差的敏感性，提高多参数反演中的稳定性； （3）构造适当的正则化参数，突出地质和偏移成像的先验信息参与度； （4）重建地震数据中的低频信息。 从计算方法上，随着过去三十年中地球物理学家对利用迭代策略的全波形反演的研究,全波形反演的实现方式上主要有梯度导引类方法（如，Tape等，2007）和牛顿类方法（如高斯-牛顿法）。 本文，首先从梯度导引类的全波形反演方法出发,分析声波方程下，全波形反演的梯度和步长的求取，并讨论迭代的优化算法。基于优化的梯度类全波形反演策略进行程序实现，在简单模型上进行数值实验。针对全波形反演相关的一些问题进行了分析。最后，对全波形反演的实用化给出了我们的策略。 2 全波形反演理论基础 地震反演问题本质上是一个非线性问题。但人们在研究非线性问题时,常常在局部甚至全局线性化,以简化问题的求解过程。地震波场的正传播过程可以记为， (1) 这里，为地震地球物理参数矢量，如速度、密度、弹性参数等，为观测到的地震数据, 描述了依赖于的地震波场正传播过程，它表达了地震波场传播的系统。在此系统下，研究特定模型中地震波的传播过程以得到数据的过程为地震正演。反之，由数据反推模型参数的过程为地震反演。反演的过程可以记为， (2) 这里，即描述了利用各种数学工具进行地震反演的过程。相对地，地震反演问题比正演问题要复杂，这是由于地震正传播得到的数据只把模型的部分信息带入了数据中。而且，模型和正问题的传播过程都只是实际情况的近似。是依赖于模型参数的非线性函数。定义误差泛函为， (3) 这里，上标表示矩阵转置。定义Fréchet微商为， (4) 在非线性问题中，Fréchet微商是依赖