4、地震勘探技术原理：反射波法、地震波传播、分辨率与信噪比等核心概念

地震勘探是石油天然气勘探中最主要的地球物理方法，其核心逻辑是通过人工激发弹性波，观测波在地下介质中的传播规律，进而推断地下岩层的几何形态与物性参数。在技术方案设计中，理解反射波法的物理机制、地震波传播规律、分辨率与信噪比的权衡关系，是后续开展采集参数论证、处理流程设计以及解释方案制定的基础。本章围绕这四个核心概念展开，帮助方案设计人员建立从物理原理到工程参数映射的完整认知链条。

4.1 反射波法的物理基础与观测系统

反射波法的本质是利用人工震源激发弹性波，当波传播至波阻抗（速度与密度的乘积）差异界面时，一部分能量反射回地表，被检波器记录。反射波的到达时间取决于反射界面的深度和上覆介质的平均速度，反射波的振幅则与界面两侧波阻抗差的大小直接相关。在方案设计中，需要明确反射波法的两个基本假设：一是地下介质呈层状或准层状结构，二是反射界面两侧存在可检测的波阻抗差异。

核心原理

反射波法的理论基础是斯奈尔定律和佐普里兹方程。对于垂直入射情况，反射系数 R 可简化为：

R = (ρ₂V₂ - ρ₁V₁) / (ρ₂V₂ + ρ₁V₁)

其中 ρ 为密度，V 为纵波速度。反射系数越大，反射能量越强。当 R 为正值时，反射波与入射波同相位；R 为负值时，反射波发生 180° 相位反转。这一特性在解释中常用于识别含气砂岩（速度降低导致负反射系数）。

观测系统设计要点

观测系统定义了炮点与检波点的空间关系，直接影响地下覆盖次数和成像质量。设计时需考虑以下参数：

• 道间距（ΔX）：需满足空间采样定理，避免空间假频。一般要求 ΔX ≤ V_min / (2·f_max)，其中 V_min 为浅层最低速度，f_max 为有效波最高频率。
• 排列长度：决定最大炮检距，需兼顾浅层折射波与深层反射波的接收。通常最大炮检距应达到目标深度的 0.8~1.2 倍。
• 覆盖次数：通过多次叠加提高信噪比，常规二维勘探为 30~60 次，三维勘探为 60~120 次。
• 方位角分布：三维采集中，宽方位角有利于各向异性分析和裂缝预测，但成本显著增加。

4.2 地震波传播规律与速度场

地震波在地下传播过程中，经历波前扩散、吸收衰减、反射、透射、折射和绕射等多种物理过程。理解这些传播规律是进行时深转换、速度分析和偏移成像的前提。在方案设计中，速度模型的精度直接决定了构造成像的可靠性和深度预测的准确性。

波传播的主要效应

• 波前扩散：球面波传播时能量密度随距离平方衰减，导致深层反射信号弱。处理中需进行几何扩散补偿。
• 吸收衰减：地层对高频成分的吸收远大于低频，导致深层信号主频降低、波形变宽。用品质因子 Q 描述，Q 值越小衰减越强。
• 薄层调谐：当薄层厚度小于 λ/4 时，顶底反射发生干涉，振幅随厚度变化呈调谐曲线。这是储层厚度预测的物理基础。
• 多次波：自由表面多次波和层间多次波会干扰一次反射波，尤其在海底或强反射界面发育区。方案设计中需评估多次波压制难度。

速度分析的关键作用

速度是地震勘探中最重要的参数之一。均方根速度（V_rms）通过 Dix 公式可转换为层速度：

V_int,n = sqrt[(V_rms,n²·t_n - V_rms,n-1²·t_n-1) / (t_n - t_n-1)]

在方案设计中，速度分析的精度受以下因素影响：

1. 覆盖次数：覆盖次数越高，速度谱能量越聚焦，拾取精度越高。
2. 最大炮检距：足够的炮检距才能观察到双曲线时距曲线的曲率变化。
3. 信噪比：低信噪比区域速度谱能量团发散，需要人工干预。
4. 各向异性：页岩或裂缝发育区存在速度各向异性，需引入各向异性参数（ε、δ、γ）进行修正。

4.3 分辨率：纵向与横向的物理极限

分辨率是衡量地震数据能够分辨最小地质体尺度的能力。纵向分辨率（时间方向）决定了能够识别的最小薄层厚度，横向分辨率（空间方向）决定了能够分辨的最小地质体宽度。两者共同约束了储层描述和构造解释的精度。

纵向分辨率

纵向分辨率由地震子波的波长决定，通常以 λ/4 作为可分辨的极限。λ = V / f，其中 V 为层速度，f 为主频。提高纵向分辨率的途径包括：

• 提高主频：通过优化激发参数（小药量、浅井深）或使用可控震源的高频扫描信号。
• 拓宽频带：低频和高频成分的保留有助于压缩子波旁瓣，提高分辨率。
• 反褶积处理：压缩子波长度，但需注意信噪比损失。

横向分辨率

横向分辨率由第一菲涅尔带半径决定：R_f = sqrt(λ·h / 2)，其中 h 为目标深度。偏移处理可以将菲涅尔带压缩至约 λ/2 的尺度，但偏移速度误差会直接影响横向聚焦效果。在方案设计中，道间距和方位角分布决定了偏移后的横向分辨率上限。

4.4 信噪比：数据质量的核心指标

信噪比（SNR）定义为有效信号能量与噪声能量的比值。地震噪声分为相干噪声（多次波、面波、折射波）和随机噪声（环境噪声、仪器噪声）。在方案设计中，信噪比是决定数据可用性的第一要素，也是处理流程中各项去噪参数选择的依据。

噪声类型与压制策略

• 面波：低频、低速、强能量，沿地表传播。通过检波器组合（组内距）和频率滤波压制。
• 多次波：周期性出现，速度低于一次波。通过叠加（速度差异）或预测反褶积压制。
• 折射波：在高速层顶界面滑行，先于反射波到达。通过切除（mute）处理。
• 随机噪声：无规律，通过覆盖次数叠加和去噪算法（如 f-x 反褶积）压制。

信噪比与覆盖次数的关系

理论上次 N 次覆盖叠加可使随机噪声信噪比提高 sqrt(N) 倍。例如，从 30 次提高到 60 次覆盖，信噪比提升约 1.41 倍。但实际提升受限于噪声相关性和静校正精度。方案设计中，应根据目标信噪比反推所需覆盖次数：

N_required = (SNR_target / SNR_single)²

其中 SNR_single 为单道信噪比，通常通过微测井或试验线数据估算。

4.5 案例拆解：某碎屑岩储层三维地震采集方案设计

背景与目标：某油田目标储层为埋深 2500~3000 m 的河道砂体，单砂体厚度 5~15 m，横向宽度 100~300 m。地质任务要求识别 5 m 以上的薄砂层，并预测含油气性。工区地表为农田，地形平坦，环境噪声中等。

初始参数与问题：初版方案采用道距 25 m、覆盖次数 60 次、最大炮检距 4000 m。经论证发现：目标层主频约 35 Hz，层速度约 3500 m/s，纵向分辨率极限 λ/4 ≈ 25 m，无法满足 5 m 薄层识别需求。

分析过程：

1. 要识别 5 m 薄层，纵向分辨率需达到 5 m，即 λ/4 ≤ 5 m，λ ≤ 20 m，主频需 ≥ 175 Hz。但该深度地层吸收严重，175 Hz 信号无法到达目标层。
2. 退而求其次，利用调谐振幅属性预测薄层厚度。调谐厚度约 λ/4 = 25 m，5 m 砂体处于调谐曲线上升段，可通过振幅标定预测厚度。
3. 为满足振幅保真度要求，需保证目标层信噪比 ≥ 5:1。根据试验线数据，单道信噪比约 0.5，所需覆盖次数 N = (5/0.5)² = 100 次。
4. 横向分辨率要求：第一菲涅尔带半径 R_f = sqrt(λ·h/2) = sqrt(100×2750/2) ≈ 370 m。偏移后横向分辨率约 λ/2 = 50 m，可分辨 100 m 宽的砂体。

方案调整：

• 道间距从 25 m 缩小至 15 m，确保空间采样密度。
• 覆盖次数从 60 次提高至 96 次（通过增加横向接收线数实现）。
• 采用可控震源低频+高频组合扫描（2~96 Hz），拓宽频带。
• 增加微测井控制点，提高近地表静校正精度。

结果与复盘：实施后目标层主频达到 45 Hz，纵向分辨率提升至约 19 m，调谐振幅属性成功预测了 5~15 m 砂体分布，钻井符合率 82%。复盘经验：薄储层识别不一定非要追求"直接分辨"，利用调谐效应和振幅属性同样可达到地质目的，但必须保证振幅保真度，这对采集方案提出了更高要求。