下载文档原格式
基于VC++的SEGY数据格式地震剖面图绘制刘继承;姚亮【期刊名称】《自动化技术与应用》【年(卷),期】2012(031)003【摘要】地震数据显示是地震资料解释的基础,将数据转换为图形的形式在屏幕上显示出来可以使其更加直观、易于理解.本文在分析了SEGY格式地震数据文件结构的基础上,利用VC++实现了地震数据的准确读取,完成了波形、面积填充和变密度这几种常见地震数据剖面显示方式的绘制.对绘制的原理和过程进行了介绍,并展示了绘制结果.%The seismic data presentation is the foundation of seismic data interpretation, converting the seismic data into graphics or images and showing them on the screen can make it more intuitively and understandable. Using Visual C++ 6.0 as development platform, seismic data exeat reading is realized, several common seismic profile displaying patterns are completed such as the waveform graph,area filling and variable density displaying based on the analysis of the the structure of SEGY file. The principles and process of drawing are introduced. At last the results of drawing are also demonstrated in the thesis.【总页数】5页(P24-27,39)【作者】刘继承;姚亮【作者单位】东北石油大学电气信息工程学院,黑龙江大庆163318;东北石油大学电气信息工程学院,黑龙江大庆163318【正文语种】中文【中图分类】TP391.411【相关文献】1.M-TEM数据拟地震剖面图绘制技术与实现 [J], 胡建平;赵昆明2.地震剖面图Wiggle曲线绘制与正振幅填充算法 [J], 刘彩云;陈忠;熊杰3.基于GPU的地震剖面图形快速绘制算法 [J], 邓博文;刘春松;吴凡贤;姚兴苗4.M-TEM 数据拟地震剖面图绘制技术与实现 [J], 赵昆明;5.基于VC的地震剖面图的绘制 [J], 孟祥琦; 王文涛; 刘腾; 干大勇; 范岩俊因版权原因，仅展示原文概要，查看原文内容请购买。
segy处理资料文件格式说明4 处理资料文件格式说明：4.1 SEG-Y 记录格式(标准)(1)卷头: 3600字节，(a) ASCII 区域: 3200字节(40条记录 x 80 字节/每条记录，特性：EBCDIC字符集，参数卡，需要转换为ASCII码后才能显示)。
(b) 二进制数区域: 400字节(3201~3600，数据类型：32位、16位的整型；特性：二进制头，记录数据体信息。
)。
3213~3214 字节—每个记录的数据道数(每炮道数或总道数)。
3217~3218 字节—采样间隔(μs)。
3221~3222 字节—样点数/每道(道长)。
3225~3226 字节—数据样值格式码1-浮点；3255~3256 字节—计量系统：1-米， 2-英尺。
3261~3262* 字节—文件中的道数(总道数)。
3269~3270* 字节—数据域(性质)：0-时域，1-振幅，2-相位谱“ * “号字为非标准定义。
(2)道记录块:(a) 道头字区: 含: 60个字/4字节整或120个字/2字节整，共240个字节，按二进制格式存放。
•SEG—Y格式道头说明：字号(4字节) 字号(2字节) 字节号内容说明1 1—2 1—4 一条测线中的道顺序号，如果一条测线有若干卷磁带，顺序号连续递增。
2 3—4 5—8 在本卷磁带中的道顺序号。
每卷磁带的道顺序号从l开始。
3 5—6 9—12 * 原始的野外记录号（炮号）。
4 7—8 13—16 在原始野外记录中的道号。
5 9—10 17—20 测线内炮点桩号(在同一个地面点有多于一个记录时使用)。
6 11—12 21—24 CMP号（或CDP号）。
（弯线=共反射面元号）7 13—14 25—28 在CMP道集中的道号(在每个CMP道集中道号从1开始)。
8—1 15 29—30* 道识别码：l=地震数据； 4=爆炸信号； 7＝计时信号；2=死道； 5=井口道；8＝水断信号；3=无效道（空道）；6＝扫描道；9…N=选择使用(N=32767)8—2 16 31—32 构成该道的垂直叠加道数(1是一道；2是两道相加;…)9—l 17 33—34 构成该道的水平叠加道数(1是一道; 2是两道叠加;…)9—2 18 35—36 数据类型：1=生产； 2=试验10 19—20 37—40 从炮点到接收点的距离(如果排列与激发前进方向相反取负值) （分米）。
地球物理数据：SEG_Y 格式BY : 张新君SEG-Y格式是地球物理探测协会开发的几种基于标准磁带格式的数据格式。
在探测和生产工业中，它被最广泛的应用于地震数据中。
然而，它是创建于1973年，其中存在着很多‘旧’的特性。
SEG-Y是为存储在IBM主机上的IBM 9磁道磁带中的单行地震数据而设计的。
现代大多数的SEG-Y变种主要是为了突破这些局限性。
SEG-Y中的一些过时的特性包括：●EBCDIC（扩充的二进制编码的十进制交换码）描述头（而不是现在流行的ASCII）●IBM 浮点指针数据（而不是现在流行的IEEE）(注: 标准总线接口; 电气和电子工程师协会)●单行存储(而不是现在普遍的3D形式)官方标准的SEG-Y格式由以下部分组成：●一个3200字节的EBCDIC描述卷头记录●一个400字节的二进制卷头记录●磁道记录由以下部分组成：⏹一个240字节的二进制磁道头⏹磁道记录谈及早期，还有许多标准形式的变种。
SEG-Y EBCDIC Textual File Reel HeaderEBCDIC卷头相当于40个IBM打孔卡。
官方设计的这些80个字节的卡片是如下所示的内容的EBCDIC等价物：12345678901234567890123456789012345678901234567890123456789012345678901234567890C 1 CLIENT COMPANY CREW NOC 2 LINE AREA MAP IDC 3 REEL NO DAY-START OF REEL YEAR OBSERVERC 4 INSTRUMENT: MFG MODEL SERIAL NOC 5 DATA TRACES/RECORD AUXILIARY TRACES/RECORD CDP FOLDC 6 SAMPLE INTERVAL SAMPLES/TRACE BITS/IN BYTES/SAMPLEC 7 RECORDING FORMAT FORMAT THIS REEL MEASUREMENT SYSTEMC 8 SAMPLE CODE: FLOATING PT FIXED PT FIXED PT-GAIN CORRELATEDC 9 GAIN TYPE: FIXED BINARY FLOATING POINT OTHERC10 FILTERS: ALIAS HZ NOTCH HZ BAND - HZ SLOPE - DB/OCTC11 SOURCE: TYPE NUMBER/POINT POINT INTERVALC12 PATTERN: LENGTH WIDTHC13 SWEEP: START HZ END HZ LENGTH MS CHANNEL NO TYPEC14 TAPER: START LENGTH MS END LENGTH MS TYPEC15 SPREAD: OFFSET MAX DISTANCE GROUP INTERVALC16 GEOPHONES: PER GROUP SPACING FREQUENCY MFG MODELC17 PATTERN: LENGTH WIDTHC18 TRACES SORTED BY: RECORD CDP OTHERC19 AMPLITUDE RECOVERY: NONE SPHERICAL DIV AGC OTHERC20 MAP PROJECTION ZONE ID COORDINATE UNITSC21 PROCESSING:C22 PROCESSING:C23C24C25C26C27C28C29C30C31C32C33C34C35C36C37C38C39C40 END EBCDIC卡中空白的空间用空白填充。
怎样保存SEGY数据供工程项目使用？（含程序）在《地震数据处理实战入门》的第五课《怎样修改保存二维和三维地震数据？》中，给大家详细讲解了怎样保存SEGY数据。
很多同学使用Python程序成功构建了新的SEGY数据，但是将数据投入工程项目使用时（如做高分辨率、断层识别等），出现了问题。
比如在商业地质勘探软件Petrel中，导入我们新建的数据资料，报了这样的错误：这是什么原因呢？01 从SEGY文件结构找原因。
SEGY地震数据一般以地震道为单位进行组织，采用SEG-Y文件格式存储。
SEG-Y格式是由SEG （Society of Exploration Geophysicists）提出的标准磁带数据格式之一，它是石油勘探行业地震数据的最为普遍的格式之一。
SEGY的数据结构比较复杂，详细介绍可以参考另外一门课《深度学习地震去噪实战》的第4节课《地震仿真噪声实战》。
这里简单介绍一下：标准SEG-Y文件一般包括三部分。
第一部分是EBCDIC卷头说明。
第二部分是二进制文件头。
第三部分是实际的地震道。
回顾一下当时保存地震数据的程序：可以看到，在上面的程序中，新保存的地震数据里面没有卷头说明和文件头的信息，只保存了地震道的数据。
这样的数据使用简单的画图程序打开没有问题，但是用到专业的处理软件做后面复杂的工作就有问题了。
02 一种简单的解决方案。
找到了问题原因，要解决起来就有了思路，就是要构造新数据的卷头说明和文件头信息。
如果大家只是对原始数据做了数值上的处理，没有改变数据尺寸，这还比较好办。
就是把原始数据的卷头说明和文件头传给新数据即可。
主要程序可以这样写：这里有几个点要说明：一是专业软件读取地震头时主要关注的是数据的尺寸，所以如果没有对数据尺寸进行修改，就可以直接复制原始数据的卷头说明和文件头信息。
二是src.text[0]读出了原始数据的卷头说明。
三是src.header读出了所有地震道的文件头信息。
03 更复杂的解决方案。
第一章帮助工具1. suhelp 显示可执行的程序和Shell脚本。
2. suname 列出SU中各项命令的名字和简短描述，以及编码的地址。
3. sudoc 得到编码的DOC列表，列出SU中各条目的在线文档。
4. sufind 在自述文档中得到信息，使用给定的字符串查找SU命令。
5. Demo演示程序：l SU软件包中有一套Shell脚本演示程序，放在目录$CWPROOT/src/demos下。
$CWPROOT/src/demos/README文件是这些演示程序的说明书。
l Making Data 演示程序显示使用susynlv程序制作合成记录炮集和共偏移距道集的基础内容。
应当特别注意演示中好的标注风格。
l Filter/Sufilter 使用实际数据处理例子演示说明消除地滚波和初至。
l Deconvolution 使用supef和其它工具简单合成脉冲道集来实例说明去混响和脉冲反褶积处理。
演示程序包括使用loops系统检验滤波参数影响的命令。
6. sukeyword 列出SU道头中的关键字 sukeyword –o7. SU基本用法是通过创建Shell脚本来实现相应的数据处理。
su/examples目录下有很多这样的例子。
第二章核心SU程序一、从磁带上读写数据下面的程序对于地球物理应用中特定的数据输入和输出任务是有用的，对内部SU数据格式也一样。
BHEDTOPAR–把二进制磁带HEaDer文件转换成PAR文件格式DT1TOSU –把Sensors & Software X.dtl GPR的地质雷达数据转成SU格式SEGDREAD –读取SEG-D磁带SEGYCLEAN - zero out unassigned portion of headerSEGYREAD–读SEG-Y磁带SEGYHDRS–为segywrite构造SEG-Y文件的ascii和二进制头文件SEGYWRITE –写SEG-Y磁带SETBHED - 设置一个SEGY二进制磁带HEaDer文件的道头字SUADDHEAD–为裸道加上头文件并设置tracl和ns道头字SUSTRIP–从道中去掉SEGY头文件SUPASTE –为已存在的数据加上已存在的SEGY头文件下面的程序可以用于一般的数据输入，输出和数据格式转换，在磁带读写中也可以使用。
SEGY格式SEGY格式是地震勘探中最常用的数据格式，所以了解SEGY格式、学会读取SEGY格式数据是非常必要的。
现将SEGY格式说明如下。
1、 SEGY格式的一般情况每个数据占4个字节（既每个数据由32位2进制数字组成）；一个SEGY数据文件的结构如下图：每个数据的4个字节的摆放顺序是：低位在前，高位在后。
如有一个十进制数据一千五百二十一，在SEGY格式中表示为：1251。
当然，SEGY格式是二进制的，这里用十进制为例，仅仅为了说明而已。
所以在读取SEGY格式的步骤有两个，Step1：读取一个32位的数据；Step2：互换该数据的第一个字节和第四个字节，互换该数据的第二个字节和第三个字节。
这时得到的数据才是确切的数据。
2、 SEG-Y 格式道头说明字(32位) 字节号说明1 1-4* 一条测线中的道顺序号。
如果一条测线有若干卷带，顺序号连续递增。
2 5-8 在本卷磁带中的道顺序号。
每卷带的道顺序号从1开始。
3 9-12* 原始的野外记录号。
4 13-16* 在原始野外记录中的道号。
5 17-20 震源点号(在同一个地面点有多于一个记录时使用)。
6 21-24 CMP号。
7 25-28 在CMP道集中的道号(在每个CMP道集中道号从1开始)。
8-1 29-30* 道识别码：1=地震数据；4=时断；7=记时；2=死道；5=井口时间；8=水断；3=DUMMY；6=扫描道；9…N=选择使用(N=32767)8-2 31-32 产生这一道的垂直叠加道数(1是一道； 2是两道相加；…)。
9-1 33-34 产生这一道的水平叠加道数(1是一道； 2是两道叠加；…)。
9-2 35-36 数据类型：1=生产；2=试验。
10 37-40 炮检距(如果是相反向激发为负值)。
11 41-44 接收点高程。
高于海平而的高程为正，低于海平面为负。
12 45-48 炮点的地面高程。
13 49-52 炮点低于地面的深度(正数)(井深)。
