如何使用segy数据绘制地震剖面

SEG-Y简介

SEG-Y格式是由SEG （Society of Exploration Geophysicists）提出的标准磁带数据格式之一，它是石油勘探行业地震数据的最为普遍的格式之一。

画图

Seismic Unix

首先需要安装Seismic Unix，安装方法（王建提供）：

#Install required packages
For fedora
yum install gcc gcc* libx* freeglut-devel mu
yum install gcc-gfortran
yum install xorg-x11-server-devel libXt-devel
 
for Ubuntu:
 sudo apt-get install build-essential
  sudo apt-get install libx11-dev
  sudo apt-get install libxt-dev
  sudo apt-get install freeglut3-dev
  sudo apt-get install libxmu-dev
  sudo apt-get install libxi-dev
  sudo apt-get install gfortran
 
#install the Seismic Unix
mkdir -p /home/wangj/programs/cwp
cd ~/programs/cwp
wget ftp://ftp.cwp.mines.edu/pub/cwpcodes/cwp_su_all_43R3.tgz
tar -zxvf cwp_su_all_43R3.tgz
 
#write the following two lines into the file ​~/.bashrc​
export CWPROOT=:/home/wangj/programs/cwp
export PATH=$PATH: /home/wangj/programs/cwp/bin
 
source ~/.bashrc
 
edit the /src/Makefile: CWPROOT =/home/wangj/programs/cwp
 
#compile
cd $CWPROOT/src
make install
make xtinstall
make finstall 
make mglinstall
make xminstall  (optional)
make sfinstall 
make utils
 
#Testing the install
suplane | suximage title="My First Plot"
suplane | suxwigb

注意：需要将路径中的/home/wangj/替换成自己的路径。
若wget不能下载成功，可以去GitHub下载安装包。

一个例子的结果如下：

注意：此脚本绘制的剖面是以Offset排列的。

Matlab

使用Matlab处理sgy数据主要需要使用SegyMAT程序包，可以在GitHub下载程序包。
使用如下脚本即可读取sgy文件并绘制剖面。

clear; close all;
[seismic_data]=ReadSegyFast(datafile);
% wiggle(seismic_data);
imagesc(seismic_data);
xlabel('TraceCount','FontSize',15,'FontWeight', 'bold');
ylabel('SampleCount','FontSize',15,'FontWeight', 'bold');

一个例子的结果如下：

注意：此脚本绘制的剖面是以Trace排列的。

SeiSee

SeiSee是一个功能强大的应用程序，以一种快速的方法来可视化SEG-Y和CST格式的地震数据。只需将sgy文件导入程序，调整参数即可，由于参数调整都是可视化的，因此使用起来比较简单。遗憾的是该程序只有Windows版本。

一个例子的结果如下：