3. 常用频谱分解方法原理对比

频谱分解的本质，是把一维地震道从“时间—振幅”域映射到“时间—频率—振幅”的联合域，从而在时间轴上定位不同厚度的地质体调谐响应。不同方法的核心差异，在于它们如何构造时频原子、如何假设信号模型，以及由此带来的时频分辨率折中。理解这一点，比记住各方法的公式更重要，因为方法选择直接决定了薄层识别能力和计算成本。

五种方法的基本原理

短时傅里叶变换（STFT）采用固定长度的窗函数沿时间轴滑动，对每个窗内信号做傅里叶变换。窗长一旦确定，时间分辨率与频率分辨率就被同时锁定：短窗时间定位好但频率分辨率差，长窗反之。其实现简单、相位保真，适合做快速扫描和频谱属性提取，但对非平稳信号的适应性最弱。

连续小波变换（CWT）用可伸缩、可平移的小波基替代固定窗，尺度大对应低频、尺度小对应高频，因而具有“低频高频率分辨率、高频高时间分辨率”的多分辨率特性。它对突变和薄层边界敏感，但尺度与频率的对应关系依赖小波选取，且冗余度高、计算量较大。

S变换（ST）可视为STFT与小波变换的折中：其窗长随频率变化，低频用长窗、高频用短窗，同时保持与傅里叶谱的直接对应关系，因此频率轴具有明确的物理意义，相位谱也可直接用于反演。它在地震信号处理中常被用于高精度时频分析和衰减估计，但计算量高于STFT。

匹配追踪（MP）不再使用固定基，而是从过完备原子库中逐次贪婪地选取与残差最匹配的时频原子，用少量原子稀疏表示信号。它能自适应地刻画地震子波和薄层调谐，时频聚焦性最好，但计算成本高，且原子库设计、迭代次数和停止准则对结果影响显著。

同步挤压变换（SST）在小波或S变换基础上，对时频面上的能量沿频率方向做重分配，把发散的能量“挤压”到瞬时频率附近，从而在保持可逆性的同时显著提高频率分辨率。它适合处理非平稳、多分量信号，但对噪声敏感，且需要选择合适的瞬时频率估计方式。

分辨率、效率与适用场景的取舍

从时频聚焦性看，匹配追踪和同步挤压最强，S变换次之，CWT再次，STFT最弱；从计算效率看则基本相反。实际生产中，若目标是快速生成分频振幅体做河道扫描，STFT或CWT已足够；若目标是精细刻画薄储层调谐厚度、做定量厚度估算，则更倾向S变换或匹配追踪；若信号非平稳性强、需要高分辨频率切片，同步挤压更有优势。

以一个教学示例说明参数影响：设采样率1 ms、主频35 Hz、目的层顶底反射时差约12 ms（对应调谐厚度约17 m，速度3000 m/s，属合理假设）。用STFT取64 ms窗时，频率分辨率约15.6 Hz，12 ms的双反射难以分开；改用CWT的Morlet小波（中心频率1 Hz、带宽参数6），在35 Hz附近有效时窗约30 ms，可分辨该调谐；若用S变换，其35 Hz处窗长约28 ms，效果接近但频率轴更易与井中频谱标定对应。这说明窗长/尺度参数必须与目标厚度匹配，而非固定取值。

实现要点与结果判读

工程实现上，STFT需注意窗函数类型与补零长度，避免频谱泄漏；CWT需明确尺度与频率的换算关系（如 f = fc / (a·dt)），否则频率轴不可解释；S变换要保证频率采样与时间采样匹配，防止高频段过采样；匹配追踪需设定原子库（如Ricker子波主频与相位网格）、最大迭代次数和残差阈值，通常以残差能量下降至初始的5%～10%作为停止条件；同步挤压需先得到小波或S变换时频面，再沿频率方向做瞬时频率重分配。

判读时不要只看单张频率切片，应结合分频振幅的横向连续性、频率随厚度的变化趋势以及井震标定共同判断。例如在35 Hz切片上出现强振幅且横向连续，而在25 Hz切片上减弱，往往指示该处存在与调谐厚度匹配的薄层。方法本身没有绝对优劣，关键是让时频原子的尺度与目标地质尺度对齐。