1. 三维地震数据体加载与地震相分析工区建立
为什么工区建立是地震相分析的第一道质量关口
地震相识别与河道砂体预测的所有后续工作——频谱分析、层位追踪、属性提取、相分类——都建立在同一个三维数据体与同一套空间参照之上。如果SEG-Y加载时坐标基准、面元网格或深度/时间域参数设置错误，后续所有属性图都会出现系统性偏移或拉伸，且很难在后期被发现。因此本例的核心不是"把文件读进来"，而是让数据体的几何信息与地质空间严格对齐，并留下可核查的完整性证据。
SEG-Y文件本身只携带道头信息（Inline/Crossline号、CDP坐标、时间采样数等），并不强制规定坐标系统。加载工具通常需要人工指定：坐标单位（米/英尺）、坐标参考系（如CGCS2000高斯投影带、UTM带）、面元尺寸、Inline/Crossline增量方向。这些参数一旦与采集设计不符，工区就会出现镜像、旋转或错位。
加载前的参数准备：从采集设计书提取关键量
在打开加载对话框之前，先把采集设计中的四个量抄到一张纸上，它们决定后续所有空间运算：
- 面元尺寸：常见 12.5 m × 12.5 m、25 m × 25 m。它决定 Inline/Xline 增量对应的实际距离。
- 坐标系统与投影带：国内陆上三维多用 CGCS2000 高斯-克吕格投影，需确认中央经线与带号；若用 UTM 则确认带号与北/南半球。
- Inline/Xline 增长方向：道头中 Inline 号增大对应北向还是南向，Xline 号增大对应东向还是西向。方向搞反会导致工区镜像。
- 时间采样与记录长度：采样间隔（1/2/4 ms）与记录长度（如 6000 ms）必须与道头一致，否则时间轴会被拉伸。
以某陆上三维为例：面元 25 m × 25 m，CGCS2000 3°带中央经线 105°E，Inline 号向北增大，Xline 号向东增大，采样 2 ms、记录长度 4000 ms。这些值在加载时逐项填入，而不是依赖软件自动识别——自动识别在道头不规范时经常给出错误结果。
加载操作的关键设置与常见陷阱
不同软件（如 Petrel、Kingdom、GeoEast、OpenDtect）的加载入口不同，但需要设置的参数类别一致。以下按参数类别说明设置逻辑与判断依据：
- 文件格式识别：确认是标准 SEG-Y（IBM 浮点或 IEEE 浮点）。若道头采样数为 0 或异常，需手动指定采样数与采样间隔，不能信任卷头。
- 道头字节映射：Inline 与 Xline 号在道头中的字节位置因采集系统而异（常见 189–192、193–196，也有 9–12、21–24 等）。加载前用十六进制查看器或软件的道头浏览功能确认实际字节位置，再建立映射。
- 坐标缩放因子：道头中的 CDP 坐标常以整数存储并带缩放因子（如坐标值 × 0.01 得米）。若忽略缩放，工区坐标会差两个数量级。
- 坐标系统指定：选择投影坐标系而非地理坐标系，输入中央经线与假东/假北。若数据已带投影坐标，不要再做二次投影。
- 工区范围裁剪：按目标区 Inline/Xline 范围设定，保留一定外扩边（如各方向外扩 20 道）以便后续层位追踪有边界控制。
最常见的三个错误：① 把 Xline 增量方向设反，导致河道走向镜像；② 坐标缩放因子漏乘，工区落在错误位置；③ 采样间隔按卷头默认值而非道头实际值，导致时间轴比例错误。这三类错误在时间切片上往往表现为"看起来正常"，直到与井位对比才暴露。
完整性核查：用可观察的证据确认工区可用
加载完成后不要直接进入解释，先做四项核查，每项都有明确的通过判据：
| 核查项 | 操作方法 | 通过判据 | 异常处理 |
|---|---|---|---|
| 道头一致性 | 统计 Inline/Xline 号范围与步长，检查是否有跳号或重复 | 步长恒为 1，无重复道，空道数为 0 | 跳号处标记为数据缺口，裁剪工区避开 |
| 空间参照 | 将工区边界与已知井位、老图件叠加 | 井位落在预期位置，边界方向与采集设计一致 | 检查投影带、缩放因子、增量方向 |
| 时间域参数 | 查看道头采样数与采样间隔，抽查初至位置 | 采样 2 ms、记录 4000 ms，初至在合理时间 | 手动重设采样参数后重新加载 |
| 可视化抽检 | 抽 3 条 Inline、3 条 Xline、3 张时间切片 | 剖面连续、无条带、振幅量级正常 | 条带多为道缺失或增益问题，回查道头 |
四项全部通过后，工区才具备进入频谱分析与层位追踪的条件。建议把核查结果（Inline/Xline 范围、面元、采样参数、坐标系统、缺口位置）记录在工区说明中，作为后续所有属性提取的空间基准。
工区建立阶段留下的参数记录，是后续河道砂体预测出现"属性图与井不符"时最快的排查入口。先核对工区参数，再怀疑属性算法，能省去大量返工。