4. 频谱分解关键参数选择与质控
频谱分解的结果是否可信，往往不取决于算法本身有多先进，而取决于四个参数是否与目标地质体匹配：时窗长度、小波基类型、尺度（频率）范围、频率采样间隔。同一份叠后数据，用 40 ms 窗和 120 ms 窗做短时傅里叶变换，得到的调谐体可能给出完全相反的甜点结论。本节围绕这四个参数的物理含义、选择依据和参数扫描质控流程展开，目标是让读者能独立完成一次可复现的参数试验并判断结果是否可用。
一、时窗长度：时间分辨率与频率分辨率的直接交换
时窗长度是频谱分解中最敏感的参数。短时傅里叶变换（STFT）中，时窗 Δt 与频率分辨率 Δf 满足不确定性关系 Δt·Δf ≥ 常数（高斯窗取等号时约为 0.5）。实践含义很直接：40 ms 窗能分辨约 25 Hz 的频率差，但时间定位好，适合追踪薄层（小于调谐厚度）的横向变化；120 ms 窗频率分辨率提高到约 8 Hz，但时间上会把相邻两个砂体“糊”在一起。
选择时窗的经验起点是目标储层时间厚度的 1.5～3 倍。例如目标砂体厚 20 ms（双程时间），可从 40 ms、60 ms、80 ms 三档开始扫描。需要注意两个陷阱：其一，时窗必须以分析点为中心对称，若用“向前开窗”会在同相轴处引入相位漂移，表现为低频能量沿层位系统性偏移；其二，时窗跨越强反射界面（如煤层顶、不整合面）时，强界面的旁瓣会污染目标频带，此时应缩短时窗或改用匹配追踪类方法。
二、小波基类型：决定时频能量的“形状”
常用小波基分三类。Morlet 小波在频率域近似高斯，时频聚焦性好，适合做连续小波变换（CWT）的频谱分解，参数为中心频率 f₀ 与带宽 σ，通常取 f₀/σ ≈ 1 保证近似零均值。Ricker 子波是地震资料自身的子波形态，用其做匹配追踪或谱反演时物理意义最直接，但频率域旁瓣较多。Gabor 窗（高斯窗 STFT）实现简单、可逆性好，是商业软件默认选项。
选择原则：做地层厚度与横向连续性分析优先 Morlet；做与井震标定直接挂钩的振幅解释优先 Ricker；需要快速批量处理、强调可逆重构时用 Gabor。小波基一旦确定，整批数据必须统一，否则不同工区或不同批次的调谐体无法对比。
三、尺度范围与频率采样：把计算量花在有效频带上
尺度范围对应分解的频率区间。地震数据有效频带通常为 8～80 Hz，但目标储层的调谐频率往往集中在 20～50 Hz。若把尺度范围设成 5～100 Hz 且采样过密，计算量成倍增加而信息量不增。推荐做法是先对目的层段做振幅谱统计，取能量高于峰值 10% 的频带作为分解范围，再向外各留 5 Hz 余量。
频率采样间隔 Δf 需满足不发生假频：对 CWT，尺度采样应保证相邻尺度的中心频率比不超过 1.05～1.1；对 STFT，Δf 应小于时窗频率分辨率的 1/2。以 60 ms 窗为例，Δf 理论分辨率约 16 Hz，采样间隔取 4～8 Hz 即可，取 1 Hz 只会带来冗余。下表给出一组可直接参考的示例参数（属教学假设，非实测）。
| 参数 | 薄层追踪场景 | 厚层边界刻画场景 | 质控要点 |
|---|---|---|---|
| 时窗长度 | 40～60 ms | 100～140 ms | 为目标厚度 1.5～3 倍，对称开窗 |
| 小波基 | Morlet（f₀/σ≈1） | Gabor 高斯窗 | 全工区统一 |
| 频率范围 | 20～60 Hz | 10～50 Hz | 覆盖有效频带 10% 能量区间 |
| 频率采样 | 4 Hz | 6 Hz | 小于分辨率 1/2，避免冗余 |
四、参数扫描与质控流程
参数不能凭经验一次定死，应通过扫描建立“参数—响应”关系。推荐流程如下：
- 选井控点：选 2～3 口有完整测井曲线的井，标定已知砂体的顶底时间与厚度，作为真值参照。
- 单参数扫描：固定其余三个参数，只变时窗（40/60/80/100 ms），观察过井地震道的分频能量是否在已知砂体位置出现单峰、峰值频率是否随厚度单调变化。
- 假频检查：把频率采样间隔减半重算，若调谐体切片出现明显差异，说明原采样不足；若几乎不变，说明采样已足够。
- 模糊检查：在已知两套薄互层之间取一道，看分频剖面能否分开两个能量团；分不开即时窗过长或频带过窄。
- 盲井验证：留一口井不参与标定，用选定参数预测砂体厚度趋势，与实钻对比，误差在可接受范围才推广到全区。
五、质控记录与可复现性
参数试验必须留痕，否则换人重做无法复现。建议每个工区保存一份参数卡，至少包含：数据版本、时窗长度及开窗方式、小波基及参数、频率范围与采样间隔、扫描档位、盲井编号与预测误差、最终采用参数及理由。下面给出一个可直接填写的示例结构（数值为教学假设）。
工区: A区三维叠后
数据版本: PSTM_2024_v3
时窗: 60 ms, 对称开窗
小波基: Morlet, f0=1.0, sigma=1.0
频率范围: 20-60 Hz, 采样 4 Hz
扫描档位: 时窗 40/60/80/100 ms; Δf 2/4/8 Hz
盲井: Well-07, 预测厚度误差 12%
采用参数: 60 ms / Morlet / 20-60 Hz / 4 Hz
质控结论: 通过, 假频与模糊检查均合格最后强调一点：参数选择没有全局最优，只有针对目标地质体尺度与资料品质的局部最优。当目标从薄层追踪切换到厚层边界刻画时，必须重新扫描时窗与频带，不能沿用上一轮参数。把参数扫描、假频检查、模糊检查和盲井验证固化为流程，才能避免因参数不当产生的假频条带与模糊能量团被误读为储层甜点。