粘滞声波波动方程偏移方法的应用与研究论文.pdfVIP

粘滞声波波动方程偏移方法的应用与研究 匡斌左建军季玉新陈世军 摘要地下介质的粘弹性吸收和频散，使常规的声波偏移方法不能使地下构造真实形态正 确成像，成像振幅也不能正确地反映地下介质的岩性和物性变化。本文将地震波的吸收和频散考 虑到偏移延拓算子中去，在波场延拓的同时，补偿地震披的吸收和频散。通过线性最优化和非线 性最优化问题的联合求解来获得粘滞声波延拓算子，保证了粘滞声波延拓算子的稳定性和计算精 度。通过设计合理的并行算法，获得理想的并行实现，从而在实用过程中解决了本算法计算量大 的技术难题。通过对数值模型和实际地震资料的处理表明，与传统的声波偏移技术相比，粘滞声 波偏移技术能够获得更真实、更精确的地下成像效果。 一、引 言 经典的弹性理论建立在精确遵守虎克定律的基础之上。弹性波动理论适用于对弹性介质 中传播的地震波的刻画和描述，但是无论是对地下岩心的实验室测试，还是油田的勘探开发实 践，都已证实地下介质或多或少地表现出非完全弹性，即粘弹性行为(衰减和频散)。特别是当 地层中含有油气时，地震波的粘弹性行为表现的尤为明显。粘弹性介质造成地震波能量衰减、 频率发散、相位畸变，并使得一些基于弹性波理论的地震处理和解释技术面I临踏进“陷阱”的危 险。 随着偏移技术的发展和地震勘探精度的不断提高，人们对偏移技术的要求由原来的构造 形态精确成像，提高到通过偏移不仅要得到精确的地下反射图象，而且成像振幅要相对正确地 反映地下构造的物性和岩性变化。AVO技术作为检测地下岩性属性最主要的工具，在许多地 区取得成功应用实例的同时，也在许多地区出现了AVO分析的“陷阱”，特别是当构造相对复 杂，绕射能量较发育时，AVO分析的可信度大打折扣。通过研究，人们发现这与反射和绕射能 量未能通过偏移很好地重新正确分配有很大关系。因而最近几年，保持振幅偏移(Hubral， 地震波振幅影响较大的因素主要有散射、透射和粘弹性吸收，因而真正的保持振幅偏移必须考 虑地震渡在粘弹性介质中所经历的衰减和频散影响。 最早的粘弹性补偿是以一维的形式出现的，这就是目前工业界仍在广泛使用的反Q滤波 不同形式的反Q滤波方法。但是，粘弹性地震波的衰减不仅与传播时间有关，而且与子渡的 频率、速度场的分布、射线路径、地下的构造展布等都有关系，仅仅利用反Q滤波对粘弹性吸 (1995)在频率空间域利用显式褶积滤波器完成了粘滞声波的叠前深度偏移。并就偏移结果对 宏Q模型误差的敏感性进行了数值分析。 一30一 利用频率、速度、吸收系数确定粘滞声波延拓算子的问题是典型的非线性最优化问题，本文 争j用线性最优化问题与非线性最优化同题的结合有效地提高了粘滞声波延拓算子的计算精度和 计算效率；为提高三维偏移的计算效率，增强本方法在实际地震资料处理中的可应用性，本文在 干Ⅱ用查表法实现显示偏移的基础上采用McClellan变换应用一维粘滞声波延拓算子实现三维波 场的向下延拓，从而很大程度地提高了三维偏移的计算效率。通过设计合理的并行算法，在微机 群上获得了理想的并行实现，从而解决了本算法在实用化过程中计算量大的技术难题。 在对粘滞声波延拓算子的误差进行分析的基础上，利用本方法分别对数值模型和实际地 震资料进行了处理，与声波偏移的结果相比，粘滞声波偏移取得了更好的成像结果。 二、基本原理和实现方法 标准线性体粘滞声波方程表示为 ’ ’ 1 M 、1(1)7 az2 av2 (z， 曼!醴铲+鱼!丝铲十鱼!丝虫堕d划+———彘=0az2 Y，2，∞)和一 式中M——复模量； p——密度。 定义复速度为 口(z，Y，#，∞)=~，丽 (2) 则式(2)的解析解为 “(∞，矗。，是，，z十△z)=e±√-‘向l‘‘：“；)Azu(坩，虹，矗y，2) (3) 定义复波数女(m)为 矗(【t．)=∞／口(∞)=矗R(∞)+旋f