2. 地震相干体基础原理
相干体技术的本质，是把三维地震数据中每个采样点周围的一小段地震道集合，压缩成一个描述“这些道彼此有多像”的数值。断层、河道边缘、地层尖灭等地质不连续处，相邻地震道的波形、振幅、相位会发生突变，相似性下降，相干值随之降低。因此理解相干体，关键要抓住三件事：用什么数学量度量相似、这个量如何由协方差矩阵与特征值得到、以及计算流程中哪些参数会显著改变结果。
相似性度量：从道间差到多道统计
最直观的相似性度量是两道之间的归一化互相关。设两道离散序列为 x(n) 与 y(n)，在时窗长度 N 内，零延迟相关系数为：
r = Σ x(n)·y(n) / sqrt( Σ x(n)² · Σ y(n)² )它取值在 −1 到 1 之间，1 表示波形完全一致。但实际相干体要处理的是分析窗内 J 道、每道 N 个采样点的道集，两道相关无法直接推广。工程上常用两类多道度量：一类基于道间能量差异（如相似系数），一类基于多道协方差矩阵的特征值结构。前者计算快、物理直观，后者对噪声与道数更稳健，也是 C1、C2、C3 等经典相干算法的共同数学基础。
协方差矩阵与特征值分解
把分析窗内的地震道排列成矩阵 D，行对应时间采样、列对应地震道，元素 dnj 为第 j 道第 n 个采样点的振幅（通常先去均值）。协方差矩阵定义为：
C = DᵀD / N （J×J 矩阵，J 为参与计算的道数）矩阵 C 的对角元素是各道能量，非对角元素是道间零延迟互相关。对 C 做特征值分解 C = VΛVᵀ，得到 J 个非负特征值 λ₁ ≥ λ₂ ≥ … ≥ λJ。特征值的分布直接反映道集的相干程度：
- 完全相干：各道波形只差一个比例因子，协方差矩阵秩为 1，只有 λ₁ 显著，其余接近 0。
- 完全不相干：各道相互独立，特征值近似相等，能量均匀分散到所有特征方向上。
- 实际地质体：介于两者之间，断层处 λ₁ 占比下降，小特征值抬升。
由此得到最常用的特征值相干（C3 类）公式：
Ec = λ₁ / Σ λⱼ （j = 1…J）该值在 0～1 之间，越接近 1 越相干。另一种形式利用特征值离散度，对弱相干更敏感。需要强调的是，特征值分解对道数 J 敏感：J 越大，矩阵维度越高，统计更稳，但横向平均效应增强，小断层容易被抹平。
计算流程与关键参数
实际计算按“逐点滑动、逐窗统计”的方式推进：对数据体中每个采样点，以其为中心取一个 inline×crossline 的道集窗口，再沿时间方向取长度 N 的时窗，构成 N×J 矩阵，去均值后求协方差矩阵、分解特征值、按公式算出该点的相干值，最后把结果写回同位置，形成相干数据体。整个过程可理解为一次三维滑动窗口滤波，只是滤波核不是加权求和，而是特征值统计。
参数选择直接决定成像效果，工程上有几条经验值得记住。时窗 N 通常取目的层主周期的 1～2 倍：以主频 30 Hz、采样率 2 ms 为例，主周期约 33 ms，时窗可取 32～64 ms（16～32 个采样点）。时窗过短，随机噪声使小特征值抬升，相干值整体偏低且抖动；时窗过长，纵向分辨率下降，薄层与相邻断层容易被合并。道数 J 常用 3×3=9 道或 5×5=25 道，前者保细节、后者抗噪但横向平滑明显，识别小断层时宜用较小孔径。若目的层倾角较大，应启用倾角/方位扫描，沿同相轴方向取窗，否则倾斜同相轴本身就会被误判为低相干。此外，计算前应做振幅均衡与去噪，避免强能量异常道主导协方差矩阵。
需要区分的是，相干值低并不等于一定是断层。地层倾角突变、采集脚印、强噪声、火山岩等特殊岩性体都可能造成低相干。因此相干体解释必须结合振幅、相位、倾角等多种属性以及区域地质认识综合判断，这也是后续断裂组合解析要解决的问题。