Henggao’s HomePage
Henggao’s HomePage/
💥
Res Seismic/
🧠
Res SU (2):SU测试
Search
🧠

Res SU (2):SU测试

1. SU读取Segy2. .画剖面图2.1.用.su画变密度剖面图2.2.用.su画变面积剖面图2.3.画自定义时窗的地震剖面图3. 文件读取——segywrite（su→segy）3.1 将Segy文件转为su文件3.2 将su数据转换成segy3.3 su文件转为bin3.4 将binary文件转为su/segy文件相关参数绘制频谱曲线

1. SU读取Segy

读.segy生成.su
# segyread tape=文件名.segy verbose=1 ns=采样点数 conv=1>文件名.su
root@seismic:/opt/cwp# segyread tape=test.segy verbose=1 | segyclean > test.su
  • filename是需要转换的文件名的名称
  • verbose 表示是否显示读取数据的每一步，=0表示不显示，=1表示显示
  • endian 表示(字，字节等的)存储机制， endian =1表示big_endian, =0表示little_endian，一般Intel等x86架构的PC机上都是little endian
notion image
 

2. .画剖面图

2.1.用.su画变密度剖面图

# suximage < 文件名.su legend=1 d2=tr.d2 d1=tr.d1 title="图名" cmap=rgb1 &
suximage < test.su legend=1 d2=tr.d2 d1=tr.d1 title="test_su" cmap=rgb1 &
【注】：cmap=rgb1时为红白蓝、默认灰黑白。
n1：时间序列的采样个数，n1=tr.ns，ns是n1的别称
n2：横向道数，n2=tr.ntr，ntr是n2的别称
d1：时间采样间隔，单位：sec，也可以修改为meter
d2：横向采样间隔，单位自定义，可以为cmp，可以为distance
f1：时间轴的起始坐标，从上方开始
f2：空间轴的起始坐标，从左侧开始
notion image

2.2.用.su画变面积剖面图

# suxwigb < 文件名.su d2=tr.d2 d1=tr.d1 title="图名" &
suxwigb < test.su d2=tr.d2 d1=tr.d1 title="test_su2" &
【注】道数较多时用变密度图
notion image

2.3.画自定义时窗的地震剖面图

# suximage < 文件名.su legend=1 d2=tr.d2 d1=tr.d1 title="图名" cmap=rgb1 x1beg=纵坐标起点 x1end=纵坐标结束点 x2be
suximage < test.su legend=1 d2=tr.d2 d1=tr.d1 title="test_su3" cmap=rgb1 x1beg=50 x1end=300 x2beg=200 x2end=600 &
notion image
用SeismicUnix（SU）读segy文件画地震剖面_seismicunix代码详解_吃最辣的锅的博客-CSDN博客
1、读.segy生成.susegyread tape=文件名.segy verbose=1 ns=采样点数 conv=1&gt;文件名.su2、用.su画变密度剖面图（红白蓝）suximage &lt; 文件名.su legend=1 d2=tr.d2 d1=tr.d1 title=“图名” cmap=rgb1 &amp;3.自定义时窗时：suximage &lt; 文件名.su legend=1 d2=tr.d2 d1=tr.d1 title=“图名” cmap=rgb1 x1beg=纵坐标起点
用SeismicUnix（SU）读segy文件画地震剖面_seismicunix代码详解_吃最辣的锅的博客-CSDN博客
https://blog.csdn.net/qq_37611077/article/details/106730278?spm=1001.2101.3001.6650.1&utm_medium=distribute.pc_relevant.none-task-blog-2~default~CTRLIST~Rate-1-106730278-blog-83148657.235^v38^pc_relevant_sort_base1&depth_1-utm_source=distribute.pc_relevant.none-task-blog-2~default~CTRLIST~Rate-1-106730278-blog-83148657.235^v38^pc_relevant_sort_base1&utm_relevant_index=2
 

3. 文件读取——segywrite（su→segy）

3.1 将Segy文件转为su文件

segyread tape='data.sgy' verbose=1 endian=0>data.su
// 其中，filename是需要转换的文件名的名称
// verbose 表示是否显示读取数据的每一步，=0表示不显示，=1表示显示
// endian =1表示服务器类型, =0表示PC机格式

3.2 将su数据转换成segy

首先需要segyhdrs来生成两个头文件
// 1) 生成头文件
segyhdrs< data.su
//执行完该命令，会在当前的文件夹下出现headr和binary的文件
//  2）格式转换
segywrite< data.su tape=data.segy verbose=1 endian=0 
// endian =1表示服务器类型, =0表示PC机格式
notion image

3.3 su文件转为bin

su转binary直接脱掉su道头即可
sustrip < in.su > ou.su

3.4 将binary文件转为su/segy文件

首先，确定好binary文件的快慢轴、大小端以及数据类型。然后，给binary文件增加道头 suaddhead；
suaddhead < indata.bin  ns=3001 > oudata.su
由于suaddhead可以增加的参数有限，在此基础上继续增加其他的道头参数 sushw。比如
sushw < indata.su  key=dt a=2 > oudata.su

相关参数

SU(Seismic Unix)画图函数supsimage使用方法_yangshun_cug的博客-CSDN博客
一、简介SU是基于类Unix平台的开源地震处理软件，不仅具有强大的数据处理能力，还能画出美观的地震剖面图。这里介绍一个重要的画图函数supsimage。supsimage是一个二维绘图函数，类似matlab的imagesc，是基于linux系统的画图函数psimage再开发的。二、supsimage函数的参数介绍我们先看一下supsimage函数的参数。$ supsimage SUPSIMAGE - PostScript IMAGE plot of a segy data
SU(Seismic Unix)画图函数supsimage使用方法_yangshun_cug的博客-CSDN博客
https://blog.csdn.net/weixin_39509073/article/details/116995598?spm=1001.2101.3001.6650.16&utm_medium=distribute.pc_relevant.none-task-blog-2~default~CTRLIST~Rate-16-116995598-blog-83148657.235^v38^pc_relevant_sort_base1&depth_1-utm_source=distribute.pc_relevant.none-task-blog-2~default~CTRLIST~Rate-16-116995598-blog-83148657.235^v38^pc_relevant_sort_base1&utm_relevant_index=20
 

绘制频谱曲线

suwind<data.su key=tracl max=60 min= 60>singleline.su
 
Copyright 2022 Henggao Geng