【seismic unix中如何给su文件添加cdp属性】

添加CDP属性到正演数据的步骤

在Seismic Unix（SU）中，正演数据通常没有道集（CDP）属性，需手动添加。以下方法适用于为模拟数据分配CDP号。

生成或准备数据
确保正演数据已加载或生成，例如使用suplane或suwavelet等命令生成模拟数据。示例生成测试数据：

suplane > synthetic.su

使用suaddhead添加CDP号
通过suaddhead直接为每道分配CDP号。以下命令为100道数据分配连续的CDP号（从1开始）：

suaddhead < synthetic.su ns=100 cdp=1,100 > output.su

参数说明：

• ns=100：设置道数为100（与数据一致）。
• cdp=1,100：CDP号从1线性递增到100。

按需设置非连续CDP号
若需非连续编号，先创建CDP列表文件（如cdp_list.txt），再通过supermute处理：

# 创建CDP列表（示例每道递增2）
seq 1 2 200 > cdp_list.txt
# 应用CDP列表
supermute < synthetic.su cdp=cdp_list.txt > output.su

验证结果
使用sugethw检查添加的属性：

sugethw cdp < output.su

或通过suximage可视化确认数据与CDP号对应关系。

其他相关操作

批量处理实际数据
若正演数据来自外部文件（如SEGY），先转换为SU格式：

segyread tape=input.segy > synthetic.su

再执行上述CDP添加步骤。

结合采集几何信息
若需更复杂的CDP分配（如基于炮检距），需使用sushw或自定义脚本处理头段字段（如scalco、gx等）。

注意：确保CDP号与数据道数匹配，避免溢出或重复。对于复杂场景，建议参考SU官方文档或结合编程（如Python）预处理头段信息。

suchw <inputdata key1=cdp key2=gx key3=sx b=1 c=1 d=2 >outdata 

可以看出cdp已经添加到trace中，更具sx和gx计算得出正确的结果。