第五章 波动方程法保幅偏移模板.ppt

图5-29 对图5-28给出的分离P波和S波进行偏移得出的成像剖面 §5.3 基于弹性波波动方程的保幅叠前深度偏移 经典的弹性理论是建立在精确遵守虎克定律的基础之上的，弹性波动理论适用于对弹性介质中传播的地震波的刻画和描述。但是无论是对地下岩芯的实验室测试，还是油田的勘探开发实践，都已经证实地下介质或多或少地表现出非完全弹性，即粘弹性行为（衰减和频散）。特别是当地层中含有油气时，地震波的粘弹性行为表现的尤为明显。粘弹性介质造成地震波能量衰减、频率发散和相位畸变，并使得一些基于弹性波理论的地震处理和解释技术面临踏进“陷阱”的危险，因而要实现保幅偏移必须考虑地震波在粘弹性介质中的衰减和频散影响。 最早的粘弹性补偿是以一维的形式出现的，这就是目前工业界仍在广泛使用的反Q滤波技术，Bickel（1985）、Gelius（1987）、Varela（1993）、Hargreves（1994）和Saatcilar（1996）分别完成了不同形式的反Q滤波。Dai和West（1994）利用伪谱法完成了粘滞声波的叠后深度偏移，将粘弹性补偿由一维发展到了二维；Mittet和Sollie（1995）在频-空域利用显式褶积滤波器完成了粘滞声波的叠前深度偏移，并就偏移结果对宏Q模型误差的敏感性做了数值分析。 利用频率、速度和吸收系数确定粘滞声波延拓算子的问题是典型的非线性最优化问题，通过线性最优化问题与非线性最优化问题的结合能有效地提高粘滞声波延拓算子的计算精度和计算效率，在确定粘滞声波延拓算子的基础上结合McClellan变换可以实现三维波场的偏移成像。下面对其基本原理及模型试算作简单的讨论。 §5.4 基于粘滞声波波动方程的保幅叠前深度偏移 一．基本原理 标准线性体粘滞声波方程表示为 其中M为复模量， 为密度。定义复速度 ，将（5-9）式变换到频率-波数域可得其通解为： （5-9） （5-10） 定义复波数， 则上、下行波场的相移算子分别表示为： 采用线性最优与非线性最优相结合的方法设计最佳的褶积延拓算子，该算子是频率、速度和粘滞性品质因子Q的函数，该算子结合McClellan变换（Hale, 1991）可以实现三维保幅叠前深度偏移，由于其计算量较大的缘故，可以通过并行算法来实现。 （5-11） §5.4 基于粘滞声波波动方程的保幅叠前深度偏移 二．模型试算 模型一 图5-30是一组粘滞声波叠后时间偏移理论模型实验。图5-30a是设计的地下构造模型，在800毫秒处有两套薄层，与下伏地层成角度不整合接触。上覆的水平盖层与下伏的隆起形成油气圈闭，圈闭之内模拟油气亮点反射。图中标出了主要的物性参数。图5-30b是图5-30a的粘滞声波自激自收响应，由于粘滞性吸收，图5-30b的频率和能量明显降低，反射相位有所畸变。图5-30c是对粘滞声波模拟的自激自收响应进行声波偏移后的结果，能量、频率和地层间的接触关系都与图5-30a有较大的误差，特别是两组倾斜反射由于衰减和频散没有得到有效的校正和补偿，致使波形畸变、能量分布不均、与上下地层的接触关系不准。图5-30d是利用粘滞声波偏移对粘滞声波自激自收响应进行偏移后的结果，可以看出：粘滞声波偏移能很好地补偿和校正地震波在传播过程中的衰减、频散和相位畸变，使反射和绕射能量高精度的聚焦和归位，从而实现真正的保幅偏移。 §5.4 基于粘滞声波波动方程的保幅叠前深度偏移 图5-30a 时间域地下构造模型 图5-30b 时间域地下构造模型的粘滞声波自激自收响应 §5.4 基于粘滞声波波动方程的保幅叠前深度偏移 图5-30c 对粘滞声波自激自收响应做声波偏移后的结果 图5-30d 对粘滞声波自激自收响应做粘滞声波偏移后的结果 §5.4 基于粘滞声波波动方程的保幅叠前深度偏移 模型二 图5-31是一组粘滞声波叠前深度偏移理论模型实验。，图5-31a是地下构造模型，图中标出了主要的物性参数。设计的道间距为20米，炮间距为40米，最小偏移距为20米，道数为150道，单边放炮，大号接收，采样率为4毫秒，记录长度为3秒，震源子波采用主频为22赫兹的雷克子波。图5-31b是CDP80处激发得到的粘滞声波单炮记录。图5-31c是利用声波叠前深度偏移对粘滞声波合成炮记录进行偏移后的结果，由于在偏移延拓过程中未能消除粘滞性吸收和频散的影响，反射波的能量和反射波的归位存在较大的误差。图5-31d是对粘滞声波合成炮记录做粘滞声波叠前深度偏移的结果，可以看到反射波的能量、频率和相位都得到了很好的恢复。 §5.4 基于粘滞声波波动方程的保幅叠前深度偏移 图5-31a 中央凹陷地质模型 图5-31c 声波叠前深度偏移的