2. 频谱分解的物理基础与地震信号特性

频谱分解之所以能用于储层预测,根本原因在于地震波在地下传播时,不同频率成分的衰减规律不同,而薄层的调谐效应又会让特定频率的振幅随厚度发生系统性变化。理解这两条物理主线,才能把频谱分解结果从“彩色图件”读成“厚度与物性信息”。本节围绕频率衰减、薄层调谐与频谱响应模式展开,并给出可操作的判读方法。

一、频率衰减:高频比低频“走得更累”

地震波在黏弹性介质中传播时,质点的非弹性形变会把弹性能转化为热能,表现为振幅随传播距离衰减,且衰减量与频率近似成正比。常用品质因子 Q 描述:

A(f, t) = A0(f) · exp(-π f t / Q)

其中 f 为频率,t 为传播时间(或双程时),Q 越小衰减越强。对同一目的层,高频成分(如 45 Hz)经过上覆地层后的相对损失明显大于低频成分(如 15 Hz)。这带来两个直接后果:一是深层、高吸收区的地震资料主频偏低、频带变窄;二是含气储层常因流体黏滞耗散产生“高频衰减异常”,在频谱上表现为高频能量相对低频能量下降,即低频阴影或高频衰减梯度异常。需要强调,衰减异常是“相对量”,必须与相邻道、相邻时窗或区域背景对比,单点绝对能量没有判读意义。

实际处理中,可沿目的层顶底开时窗做短时傅里叶变换或小波变换,计算各频率的能量或振幅谱,再构造衰减指示参数,例如高频段与低频段能量比、谱质心频率、高频衰减梯度等。示例参数(合理假设,非实测):目的层埋深 2800 m,时窗 20 ms,频带 5–60 Hz,若某道 40–55 Hz 能量与 10–25 Hz 能量之比相对背景道下降约 30%,且谱质心向低频偏移 5–8 Hz,可作为含气可能性的辅助证据,但仍需结合振幅、AVO 与地质认识综合判断。

二、薄层调谐:厚度决定“哪个频率最亮”

当储层厚度小于地震子波的四分之一波长时,顶底反射相互干涉,复合振幅不再与厚度线性相关,而呈周期性调谐。对厚度为 h、层速度为 v 的薄层,其调谐频率近似满足:

f_tune ≈ v / (4h)   (垂直入射、顶底极性相反时的主调谐近似)

例如层速度 3000 m/s、厚度 15 m,则 f_tune ≈ 50 Hz;厚度增至 30 m,调谐频率降至约 25 Hz。也就是说,薄层越薄,其调谐能量越向高频移动;越厚,能量越向低频移动。这正是频谱分解识别厚度的物理依据:在分频振幅切片上,不同厚度区域会在不同频率上出现“亮点”,把各频带振幅按频率排列,就能得到厚度变化的频谱响应模式。

需要注意,调谐频率还受子波频带、顶底反射系数极性与大小、层内多次波影响。若顶底反射系数同号(如低速层夹在高速层之间但波阻抗变化方向一致),干涉规律会改变,调谐关系需重新标定。因此实际应用前应先用井资料做正演,建立本区“厚度—频率—振幅”量板,而不是直接套用公式。

薄层厚度—调谐频率—分频振幅响应关系 地质模型(层速度 v=3000 m/s) 上覆层(高速) 储层 h=15 m 下伏层(高速) f_tune ≈ 50 Hz(高频亮) 上覆层(高速) 储层 h=30 m 下伏层(高速) f_tune ≈ 25 Hz(低频亮) 上覆层(高速) 储层 h=60 m 下伏层(高速) 顶底可分辨,调谐减弱 正演 分频振幅谱(示意) 频率 / Hz 振幅 0 15 30 45 60 75 h=15 m 峰值≈50 Hz h=30 m 峰值≈25 Hz h=60 m 顶底分离、双峰 薄层:能量向高频集中 中厚层:能量向低频移动 厚层:调谐减弱、双峰

三、不同厚度与速度对比下的频谱响应模式

把上述物理关系落到实际资料上,可归纳出三类典型响应。第一类是薄层高频亮:厚度约 8–20 m、层速度 2800–3500 m/s 时,调谐频率落在 35–60 Hz,高频分频体上出现与储层展布一致的强振幅条带,低频体上则相对平淡。第二类是中厚层低频亮:厚度 25–40 m 时调谐频率降至 20–30 Hz,低频切片更能刻画储层边界,高频切片反而因顶底干涉减弱而能量分散。第三类是厚层顶底分离:厚度超过约 1/4 波长(本例约 40–50 m)后,顶底反射可分辨,频谱上出现双峰或宽峰,调谐不再主导,此时应回到常规振幅与反演手段。

速度对比同样关键。储层与围岩的波阻抗差决定反射系数大小,进而决定调谐振幅的强弱。若储层为低速含气砂岩(如 2400 m/s)夹在 3200 m/s 泥岩中,阻抗差大、反射强,调谐特征清晰;若为致密砂岩与泥岩速度接近(如 3000 m/s 对 3100 m/s),反射系数小,调谐振幅弱,频谱分解的敏感性明显下降。因此解释前应先统计目的层段的波阻抗差与反射系数,评估该方法在本区的适用性,而不是对所有目标一律套用。

判读要点:分频振幅的“亮点”必须与调谐频率对应。若某高频切片亮点在更低频切片同样强、且形态不变,更可能是岩性界面或强反射的宽带响应,而非薄层调谐,应结合正演量板排除。

四、从物理到操作:建立本区调谐量板

实际工作中,建议按以下顺序把物理基础转化为可判读的依据。首先从井出发,提取目的层厚度、层速度与密度,计算波阻抗与反射系数;其次用井旁地震道提取子波,做不同厚度(如 5、10、15、20、30、40、60 m)的一维正演,得到各厚度的振幅谱;然后把正演谱的峰值频率与厚度拟合成关系曲线,形成“厚度—调谐频率”量板;最后将实际分频数据与量板比对,把频率域亮点翻译为厚度区间,并用井点标定误差。示例量板(合理假设):厚度 10 m 对应约 70 Hz、15 m 约 50 Hz、20 m 约 38 Hz、30 m 约 25 Hz、40 m 约 19 Hz,实际值需按本区速度与子波重新计算。

需要提醒的是,频谱分解得到的是“调谐厚度”而非真实厚度,且受噪声、子波非稳态、地层吸收与处理流程影响。合理解释应坚持三条原则:一是用井标定建立量板,不凭经验直接读厚度;二是多频带联合分析,观察峰值频率的横向迁移而非单频切片;三是与振幅、AVO、反演及沉积相认识交叉验证。把频率衰减与薄层调谐这两条物理主线贯穿到参数选择、正演标定与结果判读中,频谱分解才能真正成为储层甜点预测的定量工具,而不只是显示手段。