引言

在地震勘探与油气资源评价中，地下结构成像作为关键技术手段，其精度直接决定了地质解释与储层预测的可靠性。其中，速度模型的构建与成像方法的选择是影响地震成像质量的两大核心因素。由于地下介质往往具有高度非均质性与复杂的构造特征，如何构建准确的速度模型并选择适配的成像方法，成为地球物理勘探领域长期关注的研究热点。

传统的时间偏移成像技术依赖于对地下介质垂向速度均匀的假设，在存在显著侧向速度变化的复杂构造区，如盐丘、断裂带或速度反转层中，常会产生严重的成像畸变与定位误差，难以满足高精度地质解释需求。相比之下，深度偏移成像以及逆时深度偏移技术，凭借其对波动方程的严格求解与对横向速度变化的建模能力，在复杂地质体识别与储层精细刻画中展现出更强的优势。然而，这些方法更依赖于速度模型，其中细微误差的传播都可能引起走时偏差、相位扭曲、成像模糊等一系列问题，严重影响最终的成像质量。

为深入探讨速度模型误差对地震成像的影响机制，本文在构建包含横向速度变化与速度反转结构的二维速度模型基础上，分别应用时间偏移与深度偏移方法开展成像模拟与误差评估，系统研究不同成像算法在复杂构造条件下的表现差异。同时，构建误差分解模型，结合运动学理论、波动方程传播特性与图像域指标，定量刻画误差在成像过程中的传导路径与影响程度。此外，本文还基于MATLAB环境完成多个典型构造的模拟试验，如盐丘高速异常体、多层速度反转带等，进一步分析在不同地质情景下，速度误差引发的成像退化与叠前深度偏移（PSDM）技术潜在的失效模式。通过上述研究，本文旨在为地震成像中的速度模型误差分析与成像方法优化提供理论依据与技术参考，推动地震勘探技术向更高分辨率、更高准确率的方向发展。

1 地震速度建模与成像方法研究

1.1速度模型的建立

速度模型是地震成像的核心，其精确性直接决定了成像结果的可靠性，该模型用于描述地下介质的P波速度空间分布特征如下公式（1）：

velocity提供离散数据点用如下公式（2）表示：

典型覆盖深度为0-5000m、水平距离为0-10000m、波速1500-5000m/s。我们采用立方样条插值结合薄板样条正则化构建连续模型，以捕捉速度梯度和地质不连续性，因此模型建立的优化目标函数为如下（3）：

其中控制光滑度。第一项确保数据拟合，第二项防止过拟合，通过数值求解，得到。

例如，在模拟网格上，基础速度设为如下公式(4)：

叠加侧向波动如下公式（5）：

其中L=10000m

插值后示例点公式如下（6）：

该模型复杂度体现在二维非线性结构，能处理数据中的速度反转层，提高成像鲁棒性。

1.2 时间偏移成像模型的求解

时间偏移技术以地下介质垂直方向速度均匀为核心假设，基于数据文件中的时间域地震记录开展偏移处理，本研究采用Kirchhoff叠前时间偏移方法进行成像计算。

其积分形式表达式为如下公式（7）：

其中，为观测振幅数据。

如下公式（8）（9）为RMS速度：

包含振幅校正用来表示。

在数值求解层面，采用蒙特卡罗积分近似算法对上述积分方程进行离散化计算，可高效处理多道地震数据。本研究中设定道数为100道，时间采样间隔为0.004s。

为验证方法有效性，通过数值模拟进行求解测试：选用峰值频率为20Hz的Ricker小波作为子波源，生成合成反射信号；将模拟储层界面埋深设置为3000m。其中旅行时s，使用平均速度2500m/s。迁移后图像峰值位于s，对应深度误差约198m.

其中深度误差的转化公式如下（10）：

该方法的计算复杂度用来表示。

在横向速度变化平缓的简单地质介质中具有较高的计算效率与成像精度，但由于其忽略了地下介质的侧向速度变化，在复杂构造区域易产生假弯曲界面等成像畸变现象。

1.3 深度偏移成像模型的求解

深度偏移技术通过整合地下速度模型，能够有效处理介质的侧向非均匀性问题，从而更精准地还原复杂地质构造的空间形态。本研究采用逆时深度偏移（RTM）方法进行成像计算，该方法以波动方程为理论基础，其基本控制方程为：基于波动方程如下公式（11）：

在逆时偏移过程中，源波场按照波的传播规律进行前向传播，接收波场则进行反向传播，成像条件定义为如下公式（12）：

该成像条件通过时空域的波场互相关，可有效聚焦地下反射界面的能量。

在数值求解方面，采用有限差分（FD）方法对波动方程进行离散化计算，为平衡计算精度与效率，设定网格间距为10m，时间步长为0.001s。

波场旅行时的计算通过eikonal方程求解，公式如下（13）：

选用快速行进法（FMM）以提高求解效率，确保波前传播路径的准确性。为验证模型性能，开展模拟求解实验：

通过时间-深度转换建立初始关联，公式如下（14）：

经逆时偏移处理后，成像结果的峰值深度约为3310m，与理论储层深度相比，误差为310m（该结果基于简化假设条件；在实际应用中，通过优化速度模型与数值参数，RTM的深度误差可控制在100m以内）。

该方法的计算复杂度用来表示。虽计算成本相对较高，但凭借对波动方程的严格求解，其成像精度显著优于时间偏移技术，尤其在复杂储层，如盐丘内幕、断层破碎带等。的定位方面展现出独特优势，能够有效还原地质体的真实空间展布特征。

1.4 模型的求解

基于建立的速度模型，在坐标（5000m，2500m）处的模拟速度值为2660.3271m/s，该速度参数为后续偏移成像的精度验证提供了关键基准。

从偏移结果来看，时间偏移成像的峰值深度为0m，与理论储层深度（3000m）相比，误差约达3000m，这一显著偏差具有明确的技术必然性。时间偏移技术以垂直速度均匀为核心假设，未能纳入地下介质的横向速度变化，如上述2660.3271m/s与背景速度的差异，及复杂地质结构的影响，导致在非均匀介质条件下，旅行时计算与深度转换出现系统性偏差，最终形成较大的成像误差，这与该技术在简单层状介质中适用、复杂构造中易失真的特性完全吻合。

相比之下，深度偏移成像的峰值深度为606.0606m，误差约2393.9394m，虽仍存在偏差，但较时间偏移已有明显改善，其结果更具合理性。深度偏移通过整合速度模型中的横向变化信息，其中包含2660.3271m/s等非均匀速度参数，能够对波场传播路径进行更精准的校正，理论上可处理侧向非均质性带来的成像问题。

从数据频率角度分析，20Hz低频曲线适用于时间偏移处理：低频信号在传播过程中衰减较慢，对速度横向变化的敏感性较低，更适配时间偏移对介质均匀性的假设，虽分辨率有限，但能在简单地质条件下提供满足基础勘探需求的成像结果，符合其“较低精度、适用于简单构造”的技术定位。

然而30Hz高频曲线则更适用于深度偏移：高频信号具有更高的空间分辨率，能够捕捉复杂地质体的细微结构特征，与深度偏移处理复杂构造的技术目标高度匹配。通过深度偏移对高频信号传播路径的精准校正，可充分发挥高频数据的分辨率优势，实现对复杂储层的精细成像，这也印证了“高频数据适配深度偏移、用于复杂地质高精度探测”的技术逻辑。

由此低高频率对比可用MATLAB进一步绘制，如图1所示。

2速度模型误差建模与传播分析

2.1误差参数化与基本假设

在速度模型误差分析中，以惯性平方作为核心模型参数，将速度误差表述为。其等价形式为如下公式（15）：

为实现对误差成分的精细化解析，需将“可解释的结构误差”与“噪声型扰动”进行有效分离，为此引入如下误差分解公式（16）：

其中，用代表尺度为的低波数构造误差，该成分主要对地震波的旅行时及成像几何形态产生影响；用代表尺度为的高波数误差，其主要作用于地震波的振幅特征及散射效应。

这一误差分解方式前文中速度建模的算法逻辑及采样分辨率设置形成良好匹配，能够针对性地将不同尺度误差的传播通道进行隔离分析，为后续量化评估误差对成像结果的影响提供了清晰的理论框架。

2.2运动学通道：走时与定位误差

在高频近似下，走时满足eikonal方程如下公式（17）

对做微扰得如下公式（18）：

沿参考射线积分，其中以弧长计。可得经典一阶敏感性如下公式（19）：

对反射构成像点记作如下公式（20）：

由垂向深度转换关系得到深度定位误差如下公式（21）：

其中用表示相关速度参数。

其采用前文中给定速度场的均方根（RMS）速度或平均速度进行近似计算，这一处理方式既简化了计算过程，又能保证一定的精度。

从物理意义层面来看，低波数的通过直接驱动界面几何的误位与弯曲。低波数的构造误差会直接作用于地下反射界面，导致界面的几何形态发生错位与弯曲，进而影响成像点的定位准确性。值得注意的是，当接近临界角，或者地下介质存在较大的侧向速度梯度时，这种深度定位误差会被显著放大，使得成像结果与地质体的真实位置偏差更为明显。

2.3运动学通道：波动方程相位与振幅误差

在全波域范围内，地震波传播遵循如下声学波动方程（22）：

将代入并进行线性化处理后，可得到如下波动方程解的相位误差表达式（23）：​

而基于Born近似，振幅误差则通过散射项进行描述如下公式（24）：​

其中，代表基准格林函数，为参考波场。​

从物理意义来看，主要作用于地震信号的振幅特征与子波形态，会导致在成像条件下出现解卷积失配现象，同时伴随旁瓣残留，干扰对有效反射信号的识别。通过导致相位一同相轴错位，造成同相轴错位，进而在成像结果中产生“拉丝”“双影”以及周期跳跃等畸变，严重影响对地下地质结构形态的准确判断，尤其在复杂构造区域，这种相位偏差带来的成像干扰更为显著。

2.4 成像条件中的误差体现

时间偏移（TM）：使用做准均匀介质近似。

Kirchhoff核的旅行时以RMS速度构造。其固有假设是侧向变化弱、等效速度随时间一维变化。于是的侧向分量难以被核函数吸收，呈现为系统性曲率误差与事件拉伸。

近轴近似的横向误差估计如下公式（25）：

其中为等效波前曲率半径，大时TM易产生伪影。

对于叠前深度偏移（PSDM）与逆时偏移（RTM）而言，两种方法均显式利用进行波场外推与成像计算，主要以相位错位进入（见上式）；影响振幅与分辨率。

RTM的互相关成像条件为如下公式（26）

在时出现非零时延的能量峰值；可通过延迟域扩展成像条件显式度量错位：理想情况下能量集中在，而速度偏差使峰值偏离零延迟。

2.5 可计算的求解框架

2.5.1 数据准备与基准成像

沿用前文相关设置，用给定速度场进行TM与PSDM/RTM，输出基准像与共成像点迹集（CIG）以估计、、同相轴几何。

2.5.2 延迟域错位测量

在CIG上计算的峰值偏移（，）。经验上，与同阶，反映侧向错位。

2.5.3 几何反演（层析步）

以为观测，求解旅行时层析问题。该过程采用稀疏射线系统，无需构建显式的大型矩阵，可有效降低计算复杂度。

运用的公式如下（27）：

射线由快速行进法/射线追踪得到；以二阶平滑正则化约束，这样既能保证反演结果的稳定性，又能使计算量处于可控范围内，利用波形校正，以RTM的延迟域像作为目标，通过最小化下式（28）：

其中，对非零延迟加权，为先验。用伴随态得到并做少量外层迭代，主要修正所致的相位/振幅畸变，进一步提高成像的质量。

2.5.4误差传播与定位评估

将反演得到的代入如下公式（29）

在目标界面上报告的均值与标准差，并对TM与PSDM分别给出深度误差剖面与前文分辨率对比，以此全面评估反演结果的准确性和不同偏移方法的性能差异。

由MATLAB构建可视化分析系统，如图2所示。

通过MATLAB对研究区域的速度模型进行测算，得到速度范围为1400-5500m/s，平均速度梯度为18.47s-1，复杂度指标为0.326。这些数据为前文所分析的速度模型误差及其对成像的影响提供了具体的量化参考，也从侧面反映了该研究区域速度场的基本特征，有助于更深入地理解不同成像方法在该区域应用时所面临的误差挑战及优化方向。

2.6误差度量与判据

为全面评估速度模型误差对成像结果的影响，需建立一系列量化的误差度量指标与判据，定位误差场通过来表示。界面投影通过来表示。定位误差场通过与界面投影来表征，它能够直观地反映成像点在空间上与真实地质界面的偏离程度，为判断成像结果的整体定位准确性提供依据。

延迟集中度用如下公式（30）表示：

延迟集中度采用具体表达式进行计算，其理想值接近1。该指标值越低，表明由速度误差引起的“散能—离焦”现象越严重，即成像能量在时间或空间上的分散程度越高，会显著降低成像的清晰度和信噪比。

侧向畸变指标是时间偏移特有的度量参数，用如下公式表示（31）：

在该指标数值较大的区域，说明时间偏移方法因侧向速度变化处理不当而产生的畸变较为明显，此类区域应优先采用叠前深度偏移方法，并结合层析更新技术对速度模型进行优化，以减小侧向畸变带来的影响。

像域稳定度对比则通过比较时间偏移和叠前深度偏移在零延迟处的能量占比与峰值半高宽来实现。能量占比越高、峰值半高宽越小，表明成像焦点越锐利，成像结果的稳定性和分辨率越好。通过这种量化对比，可清晰判断不同偏移方法在成像质量上的差异，为实际应用中方法的选择提供有力支撑。

由MATLAB进行速度模型误差震荡敏感性分析，如图3所示。

由图3可得到如下定量结果：

在震荡速度模型误差统计方面，相对误差RMS为2.08%，走时误差均值达195.71ms，走时误差标准差为6056.36ms；深度误差均值为5210.9697m，深度误差标准差为86962093.49m，深度误差最大值更是达到6150396510.96m。

震荡特征参数显示，速度误差震荡系数为1.300，走时误差震荡系数为30.945，深度误差震荡系数为41.721。由此可见信噪比与震荡成像的质量相关。从震荡成像质量对比来看，时间偏移在理想情况下为0.18dB，含误差时为0.00dB，损失0.18dB；深度偏移理想情况下为0.01dB，含误差时为0.00dB，损失0.01dB。

这些数据进一步量化了速度模型误差震荡对成像结果的影响，也为不同偏移方法在误差震荡场景下的性能评估提供了具体依据。

由MATLAB进行速度模型误差对地震成像影响的震荡曲线分析，如图4所示。

在震荡误差传播分析中，速度误差与深度误差的相关系数为-0.003，震荡速度误差建模采用多频震荡叠加方法，总体RMS误差为2.08%，震荡系数为1.300，该系数用于反映震荡强度，震荡误差传播机制方面：走时误差均值为195.71ms，走时误差震荡系数为30.945；深度定位误差均值为5210.9697m，深度误差震荡系数为41.721，震荡成像影响评估显示时间偏移信噪比损失0.18dB，深度偏移信噪比损失0.01dB，且震荡特征增强了误差敏感性；震荡主要结论包括为低波数震荡误差主导大尺度几何变形，高波数震荡误差破坏局部成像分辨率；震荡误差通过非线性传播放大影响，深度偏移对震荡误差更加敏感，同时速度震荡与深度误差存在强相关，如：r=-0.003。​

这些分析结果进一步揭示了速度模型误差震荡在地震成像过程中的作用机制及影响程度，为更精准地把控成像质量提供了参考。

3 PSDM成像模型及其失效机制分析

3.1 PSDM模型

PSDM成像模型PSDM是一种通过波动方程对地震数据进行深度迁移的方法，考虑了速度模型的影响，通过反向传播波场来得到更精确的深度成像。

其数学模型可以表示为如下公式（32）：

其中，是波动场，表示位置和时间上的波动；是位置依赖的介质速度，必须与实际地下介质的速度场一致。PSDM通过将观测来自地震道的数据与一个基于速度模型的反射波传播模型进行匹配，实现图像的高精度成像。

3.2逆问题的求解

PSDM是一个逆问题，需要根据已有的观测数据来优化速度模型和成像结果。最常见的解法是使用最小二乘反演方法。

其中目标函数为如下公式（33）：

其中，是基于当前速度模型的模拟波场；是实际的地震数据；优化目标是通过不断调整速度模型，最小化模拟波场与数据的差异。优化过程通常采用梯度下降法或其他反演技术来逐步调整模型。

3.3速度模型与初始条件

PSDM的准确性与所使用的初始速度模型密切相关。如果初始模型的速度分布和实际地质结构有较大偏差，PSDM的结果可能会出现显著偏差。因此，初始速度场的准确性对最终成像结果至关重要。

3.4 PSDM在某些地质条件下可能失效的情况

PSDM是一种高度依赖速度模型的迁移方法。在一些特定的地质条件下，PSDM可能无法产生准确的结果或失效。以下是可能导致PSDM失效的地质条件：首先是速度反转层，当地下存在速度反转层，例如，低速层下方是高速层，或是高低速夹层交替时，传统的PSD成像方法可能无法有效地跟踪波的传播路径。此时，波的传播受到极大的折射影响，成像结果可能出现显著误差。其次是剧烈的速度变化与不均匀介质，在快速变化的速度场，例如，盐丘、复杂断层带等中，PSD方法可能无法精确地处理大量的波传播路径，造成成像不准确或者焦点扩展。尤其是较大的速度梯度可能导致反射波不再遵循简单的几何模式，进一步加大了成像误差。

以盐丘地形为例，用MATLAB进行仿真得到地质模型复杂度统计，如图5所示。

由图5可知，网格分辨率为200×150（间距10m），速度范围处于768-4709m/s之间，速度变化系数为0.326，盐丘几何形态为椭圆（0.3×0.2km），其空间位置位于（1.0，0.6）km处。对叠前深度偏移的失效程度评估显示，总失效深度比例达到40.7%，其中严重失效比例为2.0%。

经分析，导致PSDM失效的关键因素按影响程度排序依次为：

1. 盐丘高速异常体，其造成成像质量下降75%。
2. 多层速度反转结构，使成像精度下降50-60%。
3. 复杂断层破碎带，会导致波传播路径复杂化。
4. 以及多尺度介质非均质性。

叠前深度偏移作为复杂构造成像的关键技术，其成像精度高度依赖数据质量与地下介质的物理特性，在以下特殊条件下易出现性能衰减：首先是低频成分缺失的影响，PSDM需要较低频率的地震数据来解析深层地质结构，如果低频成分缺失，则成像精度会大大降低。缺乏低频信息使得对深层界面的成像能力下降，尤其是在深层地质反射界面时，波的传播路径不易获得。其次是多次波与散射波的干扰机制，PSDM假设一次波反射为主要成像信息源，若地层中存在强烈的多次波或散射波，PSDM成像可能被干扰，成像结果可能包含伪影或焦点不清晰。最后是高度非均质介质的适配难题，当地质介质的非均质性非常强，特别是存在高度各向异性的情况下，PSDM方法可能难以充分适应这些复杂的物理特性，导致成像错误。

4 总结

本文围绕地震成像中速度模型的构建与误差传播对成像结果的影响展开系统研究，综合采用时间偏移（TM）、叠前深度偏移（PSDM）与逆时深度偏移（RTM）等多种成像方法，在模拟不同速度场与复杂地质构造条件下，探讨了速度误差的传播机制、成像畸变特征以及偏移方法的适用性边界。通过构建二维非均匀速度模型，本文从数值模拟角度出发，验证了速度模型对成像精度的决定性作用，其中时间偏移在理想均匀介质中表现良好，但在横向速度变化剧烈的构造中易产生系统性偏差，而深度偏移特别是RTM在处理复杂波场与高频成像方面具有显著优势。论文提出了速度误差的双尺度分解方法，将误差划分为低波数结构误差和高波数震荡误差，分别对应成像几何偏差与相位振幅畸变，构建了包含延迟集中度、侧向畸变指标、定位误差场在内的成像质量评估指标体系，实现了对误差影响的量化与可视化表达。在PSDM成像的失效机制分析中，本文基于MATLAB构建典型盐丘与高速异常体模型，指出速度反转层、多尺度非均质性、多次波干扰及低频缺失是成像失效的主要因素，并量化其对图像稳定性与信噪比的影响。研究表明，速度误差通过非线性方式放大传播至成像域，在不同成像算法下表现出不同的敏感性，其中RTM具有更强的误差鲁棒性与成像聚焦能力。整体而言，本文构建的速度建模与误差传播分析框架为地震成像精度控制与方法选择提供了理论依据和技术参考，也为后续引入人工智能优化成像过程、提高复杂构造识别精度奠定了基础。

参考文献：

1. [1]段李宏,戴磊,张金陵.基于Fisher判别模型的煤层底板突水水源预测[J].工矿自动化2022,48(04):128-134.
2. [2]方党生,汪恒宇,刘久富.精确稀疏高斯变分推理非线性状态估计方法的研究[J].控制工程,2024,31(12):2274-2283.
3. [3]陈刚柱,王高峡,卜肖彤.基于网络模体和Kronecker积的分形图形研究[J].武汉纺织大学学报,2024,37(06):70-75.
4. [4]刘兰兰,敖地珍,刘艳.S-SOB矩阵的Schur补[J].青海师范大学学报(自然科学版),2024,40(03):50-54.
5. [5]吴香.分数阶Wirtinger不等式及其在时滞分数阶系统稳定性的应用研究[D].重庆交通大学,2024.