F-X域波场自适应插值法深度偏移.pdfVIP

F-X 域波场自适应插值法深度偏移1 1 1，2 1 刘文革 ，贺振华 ，黄德济 1 成都理工大学信息工程学院，四川成都 （610059 ） 2 成都理工大学“油气藏地质与开发工程” 国家重点实验室，四川成都（610059) E-mail：lwg71309@163.com 摘 要：本文通过在频率空间域波场的自适应插值，实现了复杂构造的快速准确的深度偏移 成像，并在二维空间进行了 SEG/EAGE 模型数据的深度偏移试处理。试验结果与分步傅立 叶（SSF ）、傅立叶有限差分（FFD ）深度偏移法相比较，F-X 域波场自适应插值法深度偏移 算法稳定，并能很好的处理速度场横向变化，优于SSF、FFD 法深度偏移的结果。因该法可 根据速度模型的复杂程度，自动优化选取参考速度的数量，提高了计算效率，所以可适用于 不同地质情况的深度偏移，不失为一种通用的深度成像方法。…… 实验结果证明了方法的有效性。 关键词：深度偏移 自适应插值分步傅立叶法 傅立叶有限差分法 1．引言 随着物探研究方法和计算机技术的飞速发展，许多新技术正不断应用于实践。在油气勘 探领域，复杂区域的勘探难度不断增加，为了实现精确勘探，减少风险，高精度成像的深度 偏移技术已被广泛用于复杂构造成像。目前的主流深度偏移方法是Kirchhoff积分法，但基于 波动方程波场延拓的各种方法正成为当今研究的热点。如分步傅立叶法（SSF）、傅立叶有 限差分法（FFD ）、广义屏法、频率空间域法等。这类方法具有理论完整、没有高频近似、 振幅保真等诸多优点[1] 。但其存在问题是当处理横向速度变化时，一般采用舍弃慢度场摄动 的二阶项[2] [3] 或用有限差分法对误差项进行校正 的方法，从而会造成成像精度不高或频散严 [4] 重、算法不稳定等现象。Jeno. Gazdag 1984年提出的相移加插值（PSPI ）的波场延拓方法 ， 可避免上述不利因素，但因计算效率低而制约了该方法的发展。通过深入研究，偏移实现过 程中在其理论基础上进行了方法改进，提高了计算效率，并编制了相应的软件，在模型试算 中取得了高质量的成像效果。 2 ．方法原理 相移法偏移[5] （Stolt，1978）是在频率波数域求解波动方程，由于借助快速傅立叶变换， 计算效率非常高，并且计算结果准确而无条件稳定。其理论基础是假设介质是水平层状的， 速度只随深度而发生变化。在垂向上把速度划分成若干单元，在每一单元深度内速度可视为 常数，此时波动方程简化为常系数偏微分方程。当使用波动方程进行波场外推时，可以在频 率波数域通过旋转相位角得以实现。相移法波场外推算子简单，并能适合任意地层倾角，因 而成为频率波数域和频率空间域各种偏移算法发展的基础，其主要缺陷是不能处理速度的横 向变化。 若地下介质的速度场既是横向变化的，又是垂向变化的。如何处理速度的横向变化，成 为偏移算法实现的关键[6] 。在深入研究各类偏移算法的基础之上，本文以PSPI 方法为基础， 通过改进而形成的频率空间域波场自适应插值方法能有效解决速度的横向变化。 偏移中参考速度的选取和插值方法是自适应插值法实现过程的两个关键环节： 2.1 参考速度的选取 1本课题得到高等学校博士学科点专项科研基金（项目编号：20040616001）资助 - 1 - 根据Gazdag 的研究结果，当相邻参考速度v1 ，v2 （v1v2 ）满足v21.5v1 时，频散方 程在大传播角度时中可对垂向波数kz 做近似处理 ω 2 1/ 2 kz (1−K