4. 沿层地震属性提取与属性优选
为什么“多提属性”不等于“提对属性”
沿层属性提取的关键不是把软件里能算的属性全部算一遍，而是围绕河道砂体这一目标地质体，构造一个“层位控制 + 时窗控制 + 属性族覆盖 + 敏感性筛选”的闭环。河道砂体在三维地震上的响应通常表现为：振幅类属性对砂泥岩波阻抗差异敏感，频率类属性对薄层调谐与厚度变化敏感，相位/连续性类属性对河道边界和地层终止关系敏感。因此本例沿目的层顶界（T2）建立提取基准，分别提取振幅、频率、相位三大属性族，再用相关性矩阵剔除冗余、用交会图挑选互补组合。
提取参数设定：层位、时窗与采样
本例使用一个 3D 地震数据体（示例参数：主频约 32 Hz，频带 8–55 Hz，面元 25 m×25 m，目的层 T2 顶界解释层位已通过闭合差检查）。沿层属性提取的核心参数如下：
- 基准层位：T2 顶界解释层位（作为沿层提取的参考面）。
- 提取时窗：T2 顶界向下 0～40 ms（约 1.2 个主周期），覆盖河道砂体主要调谐响应段；同时保留一个 ±10 ms 的“层位容差窗”用于稳健性检查。
- 属性族：振幅类（均方根振幅 RMS、最大绝对振幅、平均能量）、频率类（瞬时频率、谱分解 25 Hz 分量能量）、相位/几何类（瞬时相位、相干/相似性）。
- 输出网格：与地震面元一致，25 m×25 m，便于后续成图与聚类。
属性提取与相关性分析流程
下图展示了从层位输入到属性组合输出的完整数据流：层位与时窗共同约束属性提取，多属性进入相关性矩阵与交会图两级筛选，最终输出对河道砂体敏感且彼此互补的属性组合。
相关性分析与冗余剔除
对提取出的 7 个属性在目的层平面网格上逐点采样，计算两两 Pearson 相关系数 r。示例结果（教学示例数据，非实测）如下表，用于说明筛选逻辑：
| 属性对 | 相关系数 r | 判断 | 处理 |
|---|---|---|---|
| RMS 振幅 — 最大绝对振幅 | 0.94 | 高度冗余 | 保留 RMS 振幅 |
| RMS 振幅 — 平均能量 | 0.91 | 高度冗余 | 剔除平均能量 |
| 瞬时频率 — 谱分解 25 Hz 能量 | 0.72 | 中等相关 | 二者择一，保留瞬时频率 |
| RMS 振幅 — 瞬时频率 | 0.28 | 弱相关 | 互补，保留组合 |
| RMS 振幅 — 相干 | −0.35 | 弱负相关 | 互补，保留组合 |
| 瞬时相位 — 相干 | 0.61 | 中等相关 | 相干对边界更直接，保留相干 |
筛选原则：同一属性族内 |r| > 0.85 视为冗余，只保留物理意义最直接的一个；跨属性族之间 |r| < 0.5 视为互补，优先进入组合。据此得到候选组合：RMS 振幅 + 瞬时频率 + 相干。
属性交会图判读
将候选属性两两做交会图，横轴为 RMS 振幅、纵轴为瞬时频率，按已知井点（示例：2 口井钻遇河道砂体、1 口井钻遇河道边缘）着色。判读要点：
- 河道砂体井点若在交会图上聚成可分离的簇，说明该属性对具有判别力；
- 若砂体点与泥岩背景点大面积重叠，说明该属性对不敏感，应换属性或增加第三维；
- 交会图上的线性趋势若过强，提示两属性仍存在残余相关，组合信息量下降。
示例判读：RMS 振幅与瞬时频率交会图上，河道砂体点集中在“高振幅—中低频”象限，与背景泥岩点分离度较好；再加入相干作为第三维后，河道边界点进一步被区分。因此确定最终组合为 RMS 振幅 + 瞬时频率 + 相干，作为下一阶段地震相聚类划分的输入。
可复用的属性优选检查清单
| 检查项 | 通过标准 | 不通过时的动作 |
|---|---|---|
| 层位闭合差 | 目的层闭合差在允许范围内 | 重新追踪或人工修正层位 |
| 时窗长度 | 覆盖 1～1.5 个主周期 | 依据频谱分析调整时窗 |
| 族内相关性 | |r| ≤ 0.85 | 剔除冗余属性 |
| 跨族互补性 | |r| < 0.5 | 更换属性或增加属性维 |
| 交会图分离度 | 砂体点与背景点可区分 | 换属性对或检查井震标定 |
| 层位容差稳健性 | ±10 ms 内平面形态稳定 | 复查层位与信噪比 |
完成以上步骤后，输出的属性组合即可直接用于后续地震相单元划分；若交会图分离度仍不理想，应优先回到层位解释与频谱分析环节排查，而非盲目增加属性数量。