反射菲涅尔体地震层析成像和其应用研究.pdfVIP

040 SPG/SEG 深圳 2011 国际地球物理会议论文集 反射菲涅尔体地震层析成像及其应用研究 * 杨积忠 刘玉柱 董良国 (同济大学海洋与地球科学学院) 摘要 根据地震波传播的有限频理论，对于某个特定震相的观测信息，不仅射线路径上的点对该信息 具有影响，射线邻域上的其它点对该接收信息也具有影响，这种影响可以用核函数来表达。本文在前 人研究的基础上，对均匀介质情况下二维反射菲涅尔体走时层析核函数进行了深入研究，结合前人提 出的带限透射菲涅尔体边界的确定方法，将反射菲涅尔体走时层析成像理论应用于理论模型实验。反 演结果表明，反射菲涅尔体走时层析成像方法比传统的反射射线层析成像理论具有更高的反演精度和 纵向分辨率。为了解决反射层析成像中速度、深度的耦合问题，本文进一步提出了利用零偏或近偏叠 加剖面，采用速度、深度交替迭代，逐层反演的层析成像策略。理论模型实验同样验证了该策略的有 效性。 关键词 菲涅尔体 核函数 反射波 层析成像 射线 述等诸多实际问题，严重制约了波动方程层析 1 引言 在实际地震反演中的应用。 基于无限高频假设的射线理论，由于其求 根据地震波传播的有限频理论，地震信号 解速度快、算法实现相对简单，已被广泛应用 具有一定的带宽。对于某种特定震相的观测信 于 地 震 的 正 反 演 问 题 中 （ Sheriff, 1982; 息，不仅射线路径上的点对该信息具有影响， Dziewonski, 1984; Nolet, 1987; Pulliam 等, 1993; 射线邻域上的其它点对该接收信息也具有影 BilletteLambare, 1998; Aki Richards, 响，这种影响可以用核函数来表达（Liu 等, 2002 ）。但由于它无法精确描述地震波的传播规 2009 ）。理论研究表明，影响地震波传播的主要 律，导致射线层析成像反演分辨率低，在绕射 区域集中于射线邻域的第一菲涅尔体内。为了 效应比较明显的介质中的基于射线理论的偏移 简单起见，下文将第一菲涅尔体简称为菲涅尔 成像结果也不准确。 体。 基 于 波动 理 论 的波 动 方 程层 析 成 像 自从Slaney 等（1984）对均匀介质中Born （ Devaney, 1984; WuToksöz, 1987; 近似与 Rytov 近似散射层析进行了深入的分析 PrattGoulty, 1991; Woodward, 1992; 对比研究后，大量学者对 Fréchet 核函数与菲 SchusterAksel, 1993 ），从理论上讲比射线层 涅尔体层析成像方法进行了研究。WuToksöz 析具有更高的反演分辨率，但它需要计算地震 （1987）、Woodward （1992）、SniederLomax 波的正向传播以及剩余波场的反向传播，还要 （1996）、Marquering 等（1998）等在理论上 计算二者的相关沿时间的积分，计算效率非常 建立了 Rytov 波场与 Born 散射场之间的关系， 低，再加上地震子波提取困难、实际地震资料 并导出了相位与振幅扰动与 Born 散射场之间 的信噪比较低以及实际地震波传播难以准确描 的 定 量 表 达 式 。 Marquering （1999 ）、 8