ustc-seismic unix期末SU大作业（2）绘制速度模型

防盗几行，下面就是正文干货：qaq码字不易孩子赚点饭钱qaq

每组的数据和相关资料都单独放在同一个压缩包里。整套数据来自于SEG/EAGE salt model，是地震资料处理中的重要标准模型之一。数据处理的工具建议使用Seismic Unix（SU），这是科罗拉多矿业学院开发的开源处理软件，功能丰富、灵活简单，在教育界获得了广泛使用。另一个选择是较新的开源软件Madagascar，处理功能更加完善，但软件安装稍复杂。在开始作业之前，你需要有一台安装了Linux(Unix)系统的电脑、搜索下载并安装好SU或Madagascar。

请完成以下作业内容，并制作成8-10分钟的ppt报告：

1. 研究数据，明确数据道头，画出速度模型及采集系统。
2. 画出给定的5炮数据中的1炮。根据采集系统、速度模型，分析判断数据中的主要同相轴。简单描述数据的信噪比和分辨率。
3. 对5炮分别进行射线追踪，了解并分析各炮对地下区域的照明，借此明确偏移成像的优势区域和挑战区域。输出射线追踪的相关结果。注意：相邻射线的出射角度间距不能太小，以免黑作一团。要保证能看出射线分布的密度差异。
4. 对5炮数据分别进行偏移。可任选一种偏移方法。描述你所选择的偏移方法的特点，并用偏移结果加以验证。请输出各个单炮的偏移结果，分析判断是否与射线追踪的结果相呼应。
5. 对5个单炮成像结果进行叠加，从而完成一个叠前偏移流程。输出最后的成像结果，分析这一成像结果的优劣。
6. 给出上述各步骤的计算机时间，分析影响效率的主要因素。
7. 如果你使用更多炮，将得到更好的成像结果。鼓励有余力的同学使用不止5炮进行成像。除上述内容外，也鼓励在报告中进行其他力所能及的处理和分析。

1. 速度模型：seg676x210_0001.dat

横向网格点数nx=676

纵向网格点数nz=210

网格间距dx=dz=6m

Ricker子波主频f0=20Hz

激发时间延迟ts=0.1s

炮点深度zs=8

横向炮点位置：炮集数据名字的后四位，如data0001_0005.dat，位置为5。

注：

对于上述各步骤，要在ppt报告中以截图形式给出你使用的核心命令行或 shell 语句等，勿仅展示结果。

暂定考试前一周进行课堂汇报。

速度模型：seg676x210_0001.dat

python matlab su都能画

这里先展示su的画法

python matlab的画法后续写

ximage<seg676x210_0001.dat n1=210 n2=676 title="VM" legend=1&

分析：

随地层深度增加波速变大，正常波速从1500m/s递增至2500m/s,中间凹陷可能是背斜构造；

在深度400-800m且横向位置1000m-3000m范围有高速（>4000m/s）异常体,可能是盐丘；

下部相对低速的异常可能与油气储藏有关;

有大量方向与高速体两侧近一致的断层，1100m深度左右有水平断层 (dv/dz~100/(s*m) )

这速度模型图片太大了插入不进来qaq

后续我想办法加进来

反正就是一个盐丘