1、第 卷 第 期年 月地球物理学报 ，王一博 ，郑忆康 ，薛清峰等 基于 变换的全波场分离逆时偏移成像 地球物理学报 ，（）：，： ， ， （ ），（）：，：基于变换的全波场分离逆时偏移成像王一博，郑忆康，薛清峰，常旭， ， 中国科学院地质与地球物理研究所，北京 ， ， ， 摘要逆时偏移方法利用双程波算子模拟波场的正向和反向传播 ，通常采用互相关成像条件获得偏移剖面 ，是一种高精度的成像方法 但是传统的互相关成像条件会在偏移结果中产生低频噪声 ；此外 ，如果偏移速度中存在剧烈速度变化还可能进一步产生偏移假象 为了提高逆时偏移的成像质量 ，可在成像过程中先对震源波场和检波点波场分别进行波场分离 ，
2、然后选择合适的波场成分进行互相关成像 本文基于 变换 ，推导了可在偏移过程中进行上下行和左右行波场分离的高效波场分离公式以及相应的成像条件 ，结合 合成数据 ，给出了不同波场成分的互相关成像结果 数值算例结果表明 ，采用本文提出的高效波场分离算法以及合理的波场成分互相关成像条件可以获得高信噪比的成像结果 关键词逆时偏移 ；互相关成像条件 ；波场分离 ；变换： 中图分类号 收稿日期 ，收修定稿基金项目国家自然科学基金项目 （）资助 作者简介王一博 ，研究员 ，博士生导师 ，主要从事勘探地球物理和储层地质力学方面的研究工作 ： ， ， ， ， ， ， ， ， ， （） ， ， ， ， ； ； ；
3、期 王一博等 ：基于 变换的全波场分离逆时偏移成像 引言深度域偏移成像方法是油气勘探领域的一项关键技术 ，传统的偏移成像方法 ，如偏移 、单程波偏移等 ，存在一定局限性 ，无法满足复杂地区 ，复杂地表情况下的高精度成像要求基于射线理论的积分偏移法 （，； ，），求取的是波动方程高频近似解 ，对复杂构造中的高陡倾角地层不能很好地成像基于单程波算子的波动方程偏移法 ，例如相移法（，）、分步傅里叶法 （ ，）、傅里叶有限差分法 （ ，）等 ，尽管有算法稳定 、背景噪声小的优点 ，但在横向速度变化特别大的时候成像精度有限在世纪年代 ，一些学者提出了逆时偏移成像方法及实现策略（ ，；，；，），但是受当时
4、的计算能力限制 ，该方法并没有得到广泛应用近年来 ，随着计算能力的飞速发展 ，逆时偏移以其能对高陡倾角成像 ，并能够处理地震数据中多种成分 （反射波 ，绕射波 ，棱柱波等 ）的优点 ，在学术界和工业界都成为了研究热点在逆时偏移中 ，首先需要求解双程波动方程 ，构建震源波场和检波点波场 ，然后利用合适的成像条件生成成像结果一般我们采用的是有限差分法来求解常密度声波方程 ，在时间空间域进行波场的正向和反向传播最常用的成像条件是经典的互相关成像条件 （，），即震源波场和检波点波场在所有时间采样点上进行零延迟互相关后再进行叠加这种成像方法在计算和存储能力满足需求的情况下比较容易实现 ，但其主要问题是成
5、像结果中存在高振幅的低频噪声这种低频噪声是由波动方程双向传播的特点所引起的 ，会严重干扰成像质量为了避免在反射界面处产生强噪声 ，可以在波场模拟中进行针对性处理等 （）提出在波场模拟中采用无反射的双程波动方程 ，等（）提出对偏移速度进行充分光滑以降低反射波能量此外 ，等 （）提出采用滤波方法来消除低频噪声根据低频噪声的产生原理 ，如果能在成像时分离不同方向的波场能量 ，则可以从本质上解决低频噪声问题但常规的基于波场分离的逆时偏移成像方法 ，通常需要存储所有时刻的波场值并进行 傅 里 叶 变 换 ，需 要 较 大 的 存 储 量 和 计 算 量和（）利用矢量确定波场传播方向并据此来进行波场分离成
6、像等（）提出了一种基于变换的波场分离成像条件 ，可以有效去除高振幅的低频噪声等（）进一步指出 ，在偏移速度模型存在速度剧烈变化区域时 ，常规逆时偏移方法会产成成像假象为此他们利用变换构建了一个波场分离算子 ，来实现下行震源波场和上行检波点波场的互相关成像 ，并有效去除了成像假象 ，但他们只是探讨了成像过程中的上下行波场分离 ，针对的成像问题有一定局限性为此 ，本文首先分析了常规逆时偏移过程中噪声和假象的产生原理 ，然后在上下行波场分离成像条件的基础上 ，构建了基于变换的全波场分离算子 ，推导了上左 、上右 、下左以及下右四个方向的波场分离算子 、相应的互相关成像条件以及具体实现策略 ，最后针对
7、 模型给出了种波场分离成像条件对应的逆时偏移成像结果数值计算结果表明本文方法不仅可以有效地压制常规逆时偏移成像方法的噪声和假象 ，还可以根据成像目标的地质情况获取最佳成像结果 上下行波场分离及成像条件常规逆时偏移采用的是不进行波场分离的互相关成像条件 （，），其表达式为（）（，）（，）， （）这里 ，是地下某一点的空间位置 ，表示时间 ，是最大波场传播时间 ，一般选择为地震数据的最大记录时间（，）是正向传播的震源波场 ，（，）是反向传播的检波点波场这两个波场可以通过求解双程波动方程获得 ，例如采用有限差分法求解常密度声波方程为了避免存储整个正传或反传波场 ，逆时偏移可以采用以下数值步骤实现 ：
8、（）对于每一炮 ，正演震源波场 ，记录所有边界点的波场值和最后时刻的波场快照 ；（）反向传播检波点波场 ，同时根据存储的边界点波场值和最后时刻的波场快照按时间逆推震源波场 ；（）在每个时间采样点应用零延迟互相关成像条件 ，即对检波点波场和震源波场进行点乘 ，然后叠加所有时间采样点和所有炮的结果 ，获得偏移剖面逆时偏移采用双程波动方程模拟波场传播过程 ，可以准确地生成沿各个方向传播的波组成分即使是在其他成像方法中难以处理的回折波或棱柱地 球 物 理 学 报 （）卷 波 ，逆时偏移也可以对他们正确成像但是常规互相关成像条件会在成像结果中生成强振幅的低频噪声 ，尤其当偏移速度中存在强反射界面时 ，噪
9、声更为明显为了提高成像精度 ，有必要改进常规成像条件假设地下介质存在水平反射层 ，当波场遇到该反射层时 ，由于双程波动方程的特点 ，入射波场会分为向下透射和向上反射两部分 ，可以用公式表示为（，）（，）（，）， （）（，）（，）（，）， （）这里的下标和分别表示上行波和下行波 ，把公式 （）和 （）代入互相关成像条件 （）中 ，可将成像结果进一步表达为如下形式 ：（）（，）（，）（，）（，）（，）（，）（，）（，）（）等 （）指出正是最后两项 ，和，生成了逆时偏移中强振幅的低频噪声这两部分表示的是波场中沿相同方向传播的能量 ，不应该包括在互相关成像条件中为此 ，他们提出应当在公式 （）中去掉最
10、后两项等 （）进一步指出 ，在偏移速度中存在速度变化特别剧烈的区域时 ，公式 （）中的第二项将产生假象并严重干扰成像质量 ，应该在互相关成像条件中去除 ，则公式 （）所示的成像条件可进一步表示为（）（，）（，）（）对于公式 （）所示的波场分离成像条件 ，其实现难点在于如何快速有效地进行波场分离以二维地震数据的逆时偏移为例 ，传统的波场分离方法是将震源与检波点波场全部存储 ，然后变换到频率波数域 ，采用公式 （）至公式 （）实现波场分离 ，获得上行波和下行波 ：（，）（，）， ， ，（）（，）（，）， ， ，（）（，）（，）， ， ，（）（，）（，）， ， ，（）这里的（，）对应震源波场的二维傅
11、里叶变换 ，（，）对应检波点波场的二维傅里叶变换公式（）至公式 （）所示的分离方法需要较大的数据存储量 ，在实际生产中存在一定问题为此等 （）提出了一种利用变换的高效波场分离算法对于一个给定的函数（），变换具有如下性质 ： （ ）（） （ ），（）这里的表示傅里叶变换算子 ，表示变换算子 ，下标表示变换作用的变量 ，表示符号函数 ，定义为（），烅烄烆（）可以构建波场分离算子： （ ）（） （）（）该算子的傅里叶变换具有如下性质 ： （ ）（），（）该算子既可以作用于时间变量也可以作用于空间坐标或，则公式 （）中的成像条件可以表示为（ ，）（） （， ） （， ）（）经推导可将公式 （）进一步表
12、达为 （ ，）（）（）（） （ ）（）（）， （）这里的代表震源波场 ，代表检波点波场 ，满足，其中是基于波动方程计算获得的格林函数 ，表示实际地震记录 ，可证明（）（），即（）表示的是一个地震记录经过变换之后正演模拟得到的波场 全波场分离及成像条件在某些复杂情况下 ，仅使用上下行波的波场分离成像条件是不够的 ，需要进一步获取精细成像结果 ，例如在地下存在高陡倾角构造时 ，就需要将上下行波场进一步分离为左行与右行波场震源与检波点波场按照左上 、左下 、右上与右下四个方向分解后可以表示为（，）（，）（，）（，）（，），（）（，）（，）（，）（，）（，），（）期 王一博等 ：基于 变换的全波场分离
13、逆时偏移成像这里的下标和分别表示左行波和右行波把公式 （）和公式 （）代入互相关成像条件 （）中 ，将获得个成像结果可以根据地层倾角特点对个成像结果进行选择性组合 ，以获取相对合理的成像结果 ，例如对位于震源和检波点下方左侧的垂直反射层成像时 ，相对合理的成像结果可能是（）（，）（，）（）以下推导全波场分离公式 ，首先定义一个新的算子，其表达式为 （ ）（） （）， （）将该算子作用于空间变量，可得 （ ）（），烅烄烆（）则 左 右 行 波 的 波 场 分 离 算 子 可 以 表 示 为 （， ） 和 （， ）可推导获得左下 、左上 、右下 、右上四种类型波场成分表达式为（，） （， ）（）（
14、）（），（）（，） （， ）（）（）（），（）（，） （， ）（）（）（），（）（，） （， ）（）（）（）（）实际计算时 ，只需在常规逆时偏移基础上额外计算一个变换后的地震记录的正演波场 ，对其沿不同坐标方向进行变换并加以组合 ，就可以获得四个方向的波场分离结果整个计算过程不需要存储波场快照和进行傅里叶变换 ，具有很高的计算效率根据上述波场分离公式 ，可进一步推导相应的波场分离互相关成像条件 ，例如左下行震源波场和右上行检波点波场互相关成像条件为（）（，）（，） （， ） （， ）（）值得注意的是由于算子作用的震源波场或者检波点波场在频率域内对于负的均为，我们只需在时间变量上对震源波场或者检
15、波点波场作用一次算子，则公式 （）可以简化为（） （， ） （， ）（）将波场分离算子的定义公式 （）和 （）代入公式（），可得左下行震源波场与右上行检波点波场互相关成像条件的最终表达式 ：（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）在对波场进行左右行和上下行分离后 ，我们一共可以得到种互相关成像结果 ，每一种成像公式的表达形式均与公式 （）类似我们认为其中有四种成像结果是 模型相对合理的成像结果 ，第一种即为公式 （）所示 ，其余三种成像公式经推导具有如下表达式第二种是左下行震源波场和左上行检波点波场的互相关 ：（）（，）（
16、，） （， ） （， ）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）地 球 物 理 学 报 （）卷 第三种是右下行震源波场和右上行检波点波场的互相关 ：（）（，）（，） （， ） （， ），（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）第四种是右下行震源波场和左上行检波点波场互相关 ：（）（，）（，） （， ） （， ）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）（）当炮点和检波点都在地表 ，地下只有水平反射层
17、时 ，仅有下行震源波场和上行检波点波场的互相关是正确的成像结果但如果地下反射层的倾角接近时 ，就应该采用左行震源波场与右行检波点波场的互相关成像结果 ，或者采用右行震源波场与左行检波点波场的互相关成像结果我们需要针对待处理工区地质结构的特点 ，对各种成像结果进行选择性组合 ，以获取最佳成像效果 数值算例我们采用一个简单模型来验证前文推导的波场分离公式的有效性图是一个双层速度模型 ，该模型的水平方向有个网格点 ，垂直方向有个网格点 ，网格点间距为震源位置在地表水平坐标处图是时刻的震源波场快照图至图是震源波场的四个方向分离结果 ，如图所示 ，波场的左下 、左上 、右下与右上四个不同方向的波场成分均
18、得到了有效分离我们采用 数据来验证波场分离成像条件的有效性如图所示 ，该模型的水平方向有个网格点 ，垂直方向有个网格点 ，网格点间距为该数据有炮 ，炮点间距是，单炮最大道数为道 ，炮点间距为，检波点间距为，炮点和检波点深度均为，炮图 双层速度模型 图 的震源波场快照 期 王一博等 ：基于 变换的全波场分离逆时偏移成像图 四个方向的波场分离结果 ：（）左下 ；（）左上 ；（）右下 ；（）右上 （）；（）；（）；（）图 （） 速度模型；（） 偏移速度模型 （） ；（） 地 球 物 理 学 报 （）卷 图 （） 数据的常规逆时偏移剖面；（）对图 （）所示的常规逆时偏移剖面进行滤波的结果 （） ；（）
19、 （）图 左下类型的震源波场与分离后的两个检波点波场进行互相关成像的结果 （）；（） 期 王一博等 ：基于 变换的全波场分离逆时偏移成像图 模型对应的四种相对合理的成像结果（）；（）；（）；（） 地 球 物 理 学 报 （）卷 图 左上类型的震源波场与分离后的四个检波点波场进行互相关成像的结果（）；（）；（）；（） 期 王一博等 ：基于 变换的全波场分离逆时偏移成像图 右上类型的震源波场与分离后的四个检波点波场进行互相关成像的结果（）；（）；（）；（） 地 球 物 理 学 报 （）卷 图 右下类型的震源波场与分离后的两个检波点波场进行互相关成像的结果（）；（） 记录的时间长度为，时间采样间隔为
20、我们采用自由表面多次波去除之后的数据进行逆时偏移图是 数据对应的偏移速度模型 ，图是采用常规互相关成像条件的逆时偏移结果 ，图是图经过滤波之后的结果图至图分别是采用公式 （）至公式 （）所示的四种波场分离互相关成像条件获得的逆时偏移成像结果图至图是 模型对应的四个相对合理的成像结果 ，但图与图的盐丘上部均存在一些噪声 ，这是由于盐丘上边界为倾斜边界 ，而图与图对应的成像公式中存在水平方向同向传播的震源与检波点波场 ，这两个同向传播波场的互相关导致了盐丘上部出现低频噪声图至图分别是左上类型的震源波场与分离后的四种检波点波场采用互相关成像条件的逆时偏移结果图与图分别是左下类型的震源波场与左下 、右
21、下两种检波点波场采用互相关成像条件的逆时偏移结果图至图分别是右上类型的震源波场与分离后的四种检波点波场采用互相关成像条件的逆时偏移结果图与图分别是右下类型的震源波场与左下 、右下两种检波点波场采用互相关成像条件的逆时偏移结果图至图所示的种成像结果大都是偏移噪声和假象 ，主要是由同方向传播的震源与检波点波场互相关产生的 ，如，等 结论逆时偏移是地震勘探领域一种重要的成像方法 ，它成像精度高并且适用于复杂地质构造成像 ，但常规逆时偏移结果存在低频噪声 ，而且在偏移速度梯度较大时还可能存在偏移假象本文在上下行波场分离互相关成像条件的基础上 ，推导了基于变换的全波场分离公式以及相应的互相关成像条件 ，并针对 模型给出了种采用全波场分离互相关成像条件的偏移结果 ，可以有效地分离常规逆时偏移方法产生的噪声和假象本文建立的全波场分离逆时偏移成像方法计算效率高 ，易于实现 ，在实际应用中可以利用本文方法根据目标地层的特期 王一博等 ：基于 变换的全波场分离逆时偏移成像点高效地获取最优成像结果 ，这对于逆时偏移的实际应用有一定意义 ， ， ，（）： ， ， ，（）： ，（）： ， ， ， ： ，（）： ，（）： ， ，（）： ， ， ， ，（）： ， ， ，（）： ，（）： ， ，（）： ， ， ， ，（）： ， ，（）： ， ，（）：，： ， ： ， ，：（本文编辑胡素芳 ）
