ustc-seismic unix期末SU大作业（2）速度模型构建与射线追踪实战解析

1. SEG/EAGE盐丘模型与SU环境搭建

SEG/EAGE盐丘模型是国际地震勘探领域的标准测试数据集，模拟了墨西哥湾典型盐丘构造。这个模型包含676×210个网格点，横向和纵向网格间距均为6米，速度范围从1500m/s（海水层）到4500m/s（盐丘体）。我首次接触这个模型时，发现其复杂速度结构能完美测试算法的稳定性——盐丘体与周围地层的强烈速度反差会引发射线路径的剧烈弯曲。

数据准备实战：
将提供的seg676x210_0001.dat文件放入工作目录后，建议先用SU的xdip检查数据头：

sugethw < seg676x210_0001.dat key=ns,n1,n2,dx,dz

确认输出应为n1=210 n2=676 dx=6 dz=6。常见坑点是Windows系统下载的dat文件可能有字节序问题，遇到乱码时用dos2unix转换格式。

2. 速度模型可视化技巧

2.1 SU基础绘图

用ximage快速查看模型全貌：

ximage < seg676x210_0001.dat n1=210 n2=676 width=800 height=300 title="Velocity Model" legend=1 &

关键参数：

• n1：纵向采样点数（深度方向）
• n2：横向采样点数
• legend=1：显示色标栏

地质解译要点：

• 深度400-800m处出现>4000m/s高速体（白色区域），对应盐丘顶面
• 盐丘底部可见向下穿刺的"根状"构造
• 1100m附近水平断层表现为速度突变带

2.2 增强显示方案

对于复杂构造，建议使用xwigb突出细节：

xwigb < seg676x210_0001.dat n1=210 n2=676 width=800 height=400 title="Enhanced View" &

通过调整perc=99参数可以压制异常值，使盐丘边界更清晰。我曾通过对比发现，设置perc=98能更好显示盐丘与围岩的接触关系。

3. 射线追踪核心算法

3.1 最短路径法实现

在SU中可通过suray模块实现：

suray vfile=seg676x210_0001.dat xs=500 zs=10 nr=300 \ rx=500 rz=2000 dxrec=5 method=1 > rays.dat

参数解析：

• xs/zs：震源位置（单位：网格点）
• nr：接收点数量
• method=1：选择Dijkstra算法

调试经验：
当射线在盐丘边界出现"锯齿状"路径时，说明网格间距过大。此时应减小dxrec或对速度模型进行平滑处理（如使用sufilter）。

3.2 波前构建法进阶

对于多路径情况，推荐使用波前追踪：

suwfront < seg676x210_0001.dat xs=500 zs=10 \ nt=1000 dt=0.001 > wavefronts.dat

输出文件中每道代表不同时刻的波前位置。通过supsimage可生成时序动画：

supsimage < wavefronts.dat n1=210 n2=676 \ title="Wavefront Propagation" > wavefronts.ps

4. 照明分析与偏移成像

4.1 炮点照明计算

使用suillum计算照明能量：

suillum < rays.dat nz=210 nx=676 dz=6 dx=6 \ illum=illum.dat

生成的illum.dat可用ximage可视化，暖色区域表示照明充分。实际项目中，我发现盐丘侧翼常出现照明阴影区，这正是叠前深度偏移的难点区域。

4.2 偏移成像验证

选择相移法偏移进行验证：

supspmig < shot_gather.su vfile=seg676x210_0001.dat \ > migrated.su

对比照明分析结果，偏移剖面中信噪比低的区域应与照明阴影区吻合。这个交叉验证步骤在我的作业中帮助发现了速度模型构建的误差。

5. 性能优化策略

5.1 并行计算方案

对于5炮数据的批量处理，建议使用GNU Parallel：

parallel -j 4 "suray vfile=seg676x210_0001.dat xs={} zs=10 nr=300" ::: 100 300 500 700 900

-j 4表示使用4个CPU核心。实测表明，在Intel i7处理器上可将总计算时间从45分钟缩短到12分钟。

5.2 内存管理技巧

大型模型计算时可能遇到内存溢出，通过ulimit -v限制内存使用：

ulimit -v 4000000 # 限制为4GB suray vfile=large_model.dat ...

另一种方案是将模型分块处理，再用suappend合并结果。

这次作业让我深刻体会到，速度模型的质量直接决定成像效果。有个深夜调试的案例：当盐丘顶部速度值偏差5%时，偏移剖面中出现了明显的构造假象。这也印证了业界那句老话："Garbage in, garbage out"。建议同学们在提交报告前，务必用不同颜色方案反复检查速度模型的合理性。