波动方程CT技术在地震数据处理中的应用.pdfVIP

维普资讯 第 34卷 第 2期 地 球 物 理 学 报 Vo1．34．No．2 1991年 3 月 ACTA GEOPHYSICA SIN ICA M 4r． 1991 波动方程 CT技术在地震数据 处理中的应用 田 钢 唱鹤鸣 何樵登 (长春地质学院，长春 13ooz6) 摘 要 本文从畸变的 Born近似的散扰技求出发 ，给 出了利用广义 Radon变换和 Fourier积分 算子的理论反演介质间断性的原理． 将声学的广义 Radon变换与经典 RIdon 变换进行类 比，近似地导出了声学广义 Radon变换的反演公式． 本文对于反射地震学的情况，提出了一种拟线性化方法，考虑了成象点的一次散射场，从 某种程度上减少了 Born近似对弱散射的苛求． 利用同一模型的理论记录和物理实验记录的反演计算结果对提出的方法进行了验证，并 讨论了进一步提高成象精度的方法． 关键词 计算机层析技术，广义 Radon变换 ，拟线性化 ，奇异性，超声物理模型． 一 、 引 言 在医学上应用的 CT (ComputerTomograph)技术 ，以经典的 Radon 变换 (简称 CRT)为基础 ．但在其他领域 ，特别是 由于地球物理学中问题的特殊性和复杂性，远没 有达到它在医学界的应用效果，其重要原因就是 CRT的限制性条件问题． 广义 Radon 变换 (简称 GRT)的引入使得克服这一局限性成为可能，而现代数学的 Fourier积分算 子(简称 FIO)和拟微分算子的理论为解决这一问题提供了新的工具 ． 波动方程偏移技术蓬勃发展以来，常规的衍射叠加 已变得越来越不引人注意 ，原因是 它的理论基础建立在几何地震学之上．近十几年中波动方程偏移已成为地震数据处理中 必不可少的手段．但 由于引入了一个爆炸反射面模型的假说，其明确的意义却很难解释 ， 所 以对波动方程偏移来说 ，还存在着其 自身理论难 以克服的弊端． 相反，常规的衍射叠加提供了一个对地下点成象的直接过程．它建立在声学理论基 础上的理论 问题，近来已由新的反演GRT理论所证明． 本文研究的是基于 GRT理论的 波动方程 CT技术 ，它既有别于医学上基于 CRT理论的 CT 技术，也有别于天然地震等 领域的基于射线走时反演的射线 CT技术． 该方法有严密的数学基础、明确的物理意义 与直观 的几何解释 ，三者达到 比较完美的统一． 本文 1988年 l0月2§日收到，1989年 10月3日收到修改稿． 维普资讯 2期 田 钢等：波动方程 CT技术在地震数据处理中的应用 二、反演问题的线性化 频率域的声波方程 审 (，J，m)+ 彳09z (，J，)一一8(x—J)， () 这里 (，J，m)是总的声压场，J是源的位置， 为频率，C()为变化的声速 ，假设 密度是常数．设 C。()是初始速度的估计 ，这样可 以将变化的速度 c()分解成 已知和 未知的两项 c ()一 c ()+Kx)， (2) 未知的扰动 f被称为介质的(声)散射势．背景速度 Co()可以随深度 或水平方 向 变化． 取 G·(，y，m)为背景介质的格林函数，故有 vGo(x,y，∞)+ G0(，y，∞)一 一 —y)， (3) 在此基础上 ，(1)式可改写为一积分方程 (y，J，∞)一G0(