地震segy格式介绍
SEG—Y地震数据格式解析及转换方法
ＳＥＧ — ｆ ｒ ａｔｉ ｔ ｎ ｒ ｅ ｏｒ ｎｇ ｆ ｍ ａ ｅ ｓ ｉ ｔ ｏ ｔｔ ｔ ｄ ｂｙ ｔ Ｙ ｏ ｍ ｓａ ｓ ａ ｄａ ｄ ｒ ｃ ｄｉ ｏｒ ｔｏｆｓ ｉｍ ｃｄａ ａｃ ｎｓ ｉｕ ｅ ｈｅＳＥＧ ｎｓｉｕｔ ．Ｔ ｈｅＳＥＧ－ ｄａ ａ Ｉ ｔｔ ｅ Ｙ ｔ ｇ ｎｅ ａ ｅ ｔ ｉｔｙ ｉ ａ ｃ ｄ ｎｃ ｉｈ ｔ ｓ ｓ ａ ｄａ ｄ ｒ ａｌｄ ｔ ｔ ｎｄ ｒ ＥＧ－ ｄ ｔ ． Ｂｕｔｏｔ ｅ ｅ ｒ ｔ ｄ ｓ ｒｃ ｌ ｎ ｃ ｏｒ ａ ｅ ｗ ｔ ｈｉ ｔ ｎ ｒ ａ ｅ ｃ ｌｅ ｈｅ ｓ ａ ａ ｄ Ｓ Ｙ ａ ａ ｈ ｒ ＳＥＧ－ ｄａ ａ ｔ ｔ Ｙ ｔ ｈａ
的转换 过 程 。
释 系统一 般都 支持 Ｓ Ｇ Ｙ 地 震数 据 格 式 的 加 载与 Ｅ —
输 出 ， 同系统 之 间可 以借 助 Ｓ Ｇ Ｙ格 式 的地震 数 不 Ｅ — 据来共 享 和交流 成果 。
本 文 介绍 的 就是 对 标 准 和非 标 准 Ｓ Ｇ Ｙ 格 式 Ｅ — 进行灵 活解 析 ， 以及 与不 同数据 格式 之 间方便 、 直观 地转换 方 法 。
据格 式 。系统 一 般是 通 过输 入 程 序先 把 Ｓ Ｇ— 数 Ｅ Ｙ
０ 引言
Ｓ Ｇ— Ｅ Ｙ格 式不仅 是地 震勘 探 采集 的数 据格 式 ， 而且 也是 一种 通用 的地 震 数 据 交换 格 式 。处 理 、 解
据转 成 系统 内部格 式 数 据之 后 再 作 处 理 ； 通 过输 或 出程 序把 内部格 式数 据输 出成 Ｓ Ｇ Ｙ数 据 ， Ｅ — 以便输 入到其 它 系统 做 进 一 步处 理 或 解 释 。这 里 的输 入 、 输 出过程 就是 Ｓ Ｇ Ｙ 数 据 与 另一 种 数 据 格 式之 间 Ｅ —
SEGY格式说明
SEG Y数据格式译稿(1)(2007-11-13 22:41:45)标签：IT/科技分类：专业SEG Y修订版1数据交换格式1SEG技术标准委员会2版本1.0，2002年5月12001,勘探地球物理学家学会，版权所有2编者：Michael W. Norris 和Alan K. Faichney目录1.简介2.概述2.1. 不变的条目2.2. 修订版0到修订版1的变化2.3. 注释2.4. 监管机构2.5. 致谢3.SEG Y文件结构3.1. 记录介质3.2. 文件结构3.3. 数字格式3.4. 变道长3.5. 坐标4.原文文件头5.二进制文件头6.扩展原文文件头6.1. 扩展原文文件头结构6.2. 结尾文本段6.3. 文本段示例7.数据道7.1. 道头7.2. 道头数据附录A. 写SEG Y数据到磁盘文件附录B. SEG Y磁带标签附录C. 磁带上的SEG Y文件块附录D. 扩展原文文本段D-1. 位置数据D-1.1 位置数据文本段D-1.2 位置数据文本段示例D-2. 面元网格定义D-2.1 面元网格定义文本段D-2.2 面元网格定义文本段示例 D-3. 资料地理范围和覆盖区域D-3.1 资料地理范围文本段D-3.2 资料地理范围文本段示例D-3.3 覆盖区域文本段D-3.4 覆盖区域文本段示例此例基于图3D-4. 数据取样测量单位D-4.1 数据取样测量单位文本段D-4.2 数据取样测量单位文本段示例D-5. 处理历史D-5.1 处理历史文本段D-5.2 处理历史文本段示例D-6. 震源类型/方位D-6.1 震源类型/方位文本段D-6.2 震源类型/方位文本段示例D-7. 震源测量单位D-7.1 震源测量单位文本段D-7.2 震源测量单位文本段示例附录E. 文字格式数据附录F. EBCDIC码和ASCII码附录G. 参考文献图片图1. 带N个扩展原文文件头记录和M道记录的SEGY文件字节流结构图2. 面元网格定义图3. 地震调查的数据范围和覆盖区域变化表格表1. 原文文件头表2. 二进制文件头表3. 道头表4. SEG Y磁带标签表5. 位置数据文本段表6. 面元网格定义文本段表7. 资料地理范围文本段表8. 覆盖区域文本段表9. 数据取样测量单位文本段表10. 处理历史文本段表11. 震源类型/方位文本段表12. 震源测量单位文本段表13. IBM 3270 字符集参考CH 10，GA27－2837－9，1987年4月1．简介最早的SEG Y数据交换格式（修订版0，参考页数45）自1975年出版以来在地球物理行业得到广泛的使用。
地震数据文件格式
炮点的野外一次静校正值（ms）
接收点的野外一次静校正值（ms）
总野外一次静校正量 (若未用静校时为零，ms)
延迟时间—A，以ms表示 时间延迟—B，以ms表示 延迟记录时间，以ms表示 起始切除时间（ms）
结束切除时间（ms） 本道的采样点数
本道的采样间隔，以us表示
野外仪器的增益类型
73—88字节中坐标的比 例因子=1，土10， 土100，土1000 土10000。
如果为正，乘以因子； 如果为负，则除以因子
接收点坐标—Y（分米）
（如果坐标单位是弧度·秒 ；X值代表径度，Y值代表纬度；正值代表格林 威治子午线东或者赤道北的秒数。负值则为西或者南的秒数）
坐标单位； 1=长度(米或者英尺)； 2=弧度·秒 接收点下风化层速度 （低速带速度，m/S） 接收点下次风化层速度 （降速带速度，M/S） 震源处的井口时间（ms）
trace number of sweep channel (扫描辅助道数) sweep trace taper length at start if tapers.
sweep trace taper at the end (扫描类型码) sweep trace taper type code: (扫描斜坡类型码) 1 = linear; 2 = cos-squared ;3 = other
接收点的地面高程。 高于海平面为正， 低于海平面为负(cm) 炮点的地面高程(cm)
炮井深度(正数，cm) 接收点基准面高程(cm)
炮点基准面高程(cm) 炮点的水深(cm)
接收点的水深(cm)
炮点坐标—X（分米） 炮点坐标—Y（分米） 接收点坐标—X（分米）
41一68字节中高程 和深度的比例因子=l， 土10，土100，土1000 或者 土10000。 如果为正，乘以因子； 如果为负，则除以因子
地震资料的处理体会Y
地震资料的处理及特殊处理方法一、地震资料的纪录格式SEG-D、SEG-Y格式最为常用：Y格式的道序如下：这里的间隙是物理间隙，不是指字节其中卷头由2部分组成，第一部分40张卡片，每张80列，即80字节，这样第一部分共计3200字节；第二部分为二进制卷头编码数据，共400字节长。
卡片上纪录着工区名、测线号、采样间隔、纪录长度等信息。
可打印出来。
二进制卷头（400字节）纪录信息如下：1-4 byte:作业鉴别号5-8 byte ：测线号9-12 byte：磁带卷号13-14 byte：每个纪录内得道数17-18 byte：采样间隔21-22 byte：每道的样点个数25-26 byte：数据样点格式1、浮点数（4字节）2、定点数（4字节）3、定点数（2字节）4、带增益码的定点数（4字节）27-28 byte：CDP号29-30 byte：道选排方式1、炮序2、CDP选排3、单次覆盖4、水平叠加55-56 byte：测量系统1、公制，m2、英制，in在卷头之后为地震道，每一个地震道由道头加数据两部分组成。
道头和数据之间没有间隙，道头的长度为240字节长，道头后面跟着数据，一般为32位浮点数。
240个字节安排如下：1-4字节：该道在测线的顺序号4-8字节：该道在一卷磁带内的顺序号9-12字节：原始野外纪录号（炮号）13-16字节：该道在原始炮中的道号21-24字节：cdp号29-30字节：道鉴别编码1、地震道号2、死道3、全零道4-8辅助道115-116字节：本道采样点数117-118字节：采样间隔二、反Q滤波由于大地的滤波作用，地层可作为一个滤波体系。
从能量衰减意义上看，Q 值就是存储能量与消耗能量的比率，它代表了岩石的一种固有属性。
地震波传播过程中，频率愈高，时间愈长，衰减愈快。
反Q滤波就是对地震纪录进行频率、振幅、相位的补偿。
三、提高剖面信噪比方法1、多项式拟合技术地震纪录在单炮纪录和水平叠加剖面上都表现为能量渐变和波形渐变，随深度的的增加，道与道之间的相关性越来越好。
地震记录格式说明
§4 处理资料文件格式说明：4.1 SEG-Y 记录格式(标准)(1)卷头: 3600字节(a)ASCII 区域: 3200字节(40条记录x 80 字节/每条记录)。
(b)二进制数区域: 400字节(3201~3600)。
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 从炮点到接收点的距离(如果排列与激发前进方向相反取负值) （分米）。
SEGY道头详细说明
101-102
检波器组的校正量（毫秒）
103-104
应用的总静校正量（毫秒）（如没有应用静校正量为零）
105-106
延迟时间A――以毫秒表示的240字节道识别头的结束和时间断点之间的时间。当时间断点出现在头之后，该值为正；当时间断点出现在头之前，该值为负。时间断点是最初脉冲，它由辅助道记录或由其他记录系统指定。
3=2字节，两互补整数。4=4字节带增益定点数（过时，不再使用）。5=4字节IEEE浮点数。6=现在没有使用。7=现在没有使用。8=1字节，两互补整数。
3227-3228
道集覆盖次数――每个数据集的期望数据道数（例如CMP覆盖次数）。强烈推荐所有类型的数据使用
3229-3230
道分选码（即集合类型）：-1=其他（应在用户扩展文件头文本段中解释）0=未知。1=同记录（未分选）。2=CDP道集。3=单次覆盖连续剖面。4=水平叠加。5=共炮点。6=共接收点。7=共偏移距。8=共中心点。9=共转换点。强烈推荐所有类型的数据使用。
脉冲极化码：1=压力增大或检波器向上运动在磁带上记作负数。2=压力减小或检波器向下运动在磁带上记作正数。
3259-3260
可控源极化码：地震信号滞后引导信号：1=337.5°-22.5°。2=22.5°-67.5°。3=67.5°-112.5°。4=112.5°-157.5°。5=157.5°-202.5°。6=202.5°-247.5°。7=247.5°-292.5°。8=292.5°-337.5°。
77-80
震源坐标――Y
81-84
检波器组坐标――X
85-88
检波器组坐标――Y
89-90
坐标单位：1=长度（米或英尺）。2=弧度秒。3=小数度。4=度，分，秒（DMS）。
地震数据格式知多少
地震数据格式这么多，你晕了吗？王清振每个物探人员都要和地震数据打交道，但是很多同学常常被眼花缭乱的地震数据格式迷得七荤八素。
SEGA、SEGB、SEGC、SEGD、SEGY、SU，工作站格式、微机格式、整形、浮点型、大头(Big-Endian)、小头(Little-Endian)、IBM、IEEE等等，这些名称连起来可以绕脑门好几圈，直接把你绕蒙圈。
今天就为大家抽丝剥茧，聊聊地震数据格式那些事。
SEGA、SEGB、SEGC、SEGD、SEGYSEGA、SEGB、SEGC、SEGD、SEGY这几种格式一看就是亲兄弟，它们系出名门，都是美国勘探地球物理学会(SEG)推荐的几种数字磁带记录格式。
如果地球物理界也是一个江湖，那么SEG就是该江湖中的武林大会了，该机构推荐的格式基本就是江湖标准，应用那是相当的广泛。
SEG-A和SEG-B记录格式简称A格式和B格式，是SEG1967年推荐的两种数字磁带记录格式。
分别适用于21轨一英寸磁带和九轨半英寸磁带。
A格式目前早已随着一英寸磁带的淘汰而被淘汰。
SEG-C是1972年推荐的一种九轨半英寸磁带记录格式，它以时序方式对数据多路编排记带，与B格式的区别在于所记录的数据是以32位浮点的IBM格式记录的。
SEGD和SEGY都是1975年推出的新的数据格式，也是目前应用最多的两种，SEGD在野外采集时用的比较多，载体多为磁带，SEGY多用在室内数据传递，载体以磁盘居多。
其实对我们大部分从业人员来说，ABCD都只是传说，SEGY才是我们熟悉的家常菜。
SU格式SU全称Seismic Unix，是科罗拉多矿院CWP（Center for Wave Phenomena）实验室开发的一套开源地球物理数据处理系统。
CWP在江湖上具有响当当的名声，地位比肩少林武当，门下人才辈出，John W. Stockwell Jr., Jack K. Cohen, Einar Kjartansson, and Joshua (Shuki) Ronen等人也因为在SU中做出的贡献与2002年被SEG授予特殊贡献奖, 目前Dave Hale正带领着团队持续维护并发展着SU。
C# 读取sgy代码
sgyrinfo.ntrpr = rp.ReadUInt16();
input.Seek(2, SeekOrigin.Current);
sgyrinfo.hdt = rp.ReadUInt16();
temp = data32[1];
data32[1] = data32[2];
data32[2] = temp;
}
public unsafe void convert(int* a)
{
temp = data32[1];
data32[1] = data32[2];
data32[2] = temp;
}
// 写文件
public unsafe void write(string path,float **data,ushort interval,ushort tracenum,ushort samplenum)
if (testformat >= sgyrinfo.hns)
isWS = true; //为真则是在工作站下的写的
input.Seek(2, SeekOrigin.Current);
sgyrinfo.format = rp.ReadInt16();
FileStream input = new FileStream(path, FileMode.Open);
BinaryReader rp = new BinaryReader(input);
//获取卷头信息&&判断文件是PC机写的还是工作站写的
读写SGY格式的c#代码
首先说一下，什么是segy？
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
地震segy格式介绍
地震SEG-Y格式
SEGY格式是地震勘探中最常用的数据格式，所以了解SEGY格式、学会读取SEGY格式数据是非常必要的。
现将SEGY格式说明如下。
1、 SEGY格式的一般情况
每个数据占4个字节（既每个数据由32位2进制数字组成）；
每个数据的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 炮点低于地面的深度(正数)(井深)。
14 53-56 接收点的基准面高程。
15 57-60 炮点的基准面高程。
16 61-64 炮点的水深。
17 65-68 接收点的水深。
18-1 69-70 对41-68字节中的所有高程和深度应用了此因子给出真值。
比例因子=1，±10，±100，±1000或者±10000。
如果为正，乘以因子；如果为负，则除以因子。
18-2 71-72 对73-88字节中的所有坐标应用了此因子给出真值。
比例因子=1，±10，±100，±1000或者±10000。
如果为正，乘以因子；如果为负，则除以因子(在GRISYS中为10)。
19 73-76 炮点坐标--X |- 如果坐标单位是弧度的秒，X值代表
20 77-80 炮点坐标--Y | 径度，Y值代表纬度。
正值代表格林
21 81-84 检波点坐标--X | 威治子午线东或者赤道北的秒数。
负
22 85-88 检波点坐标--Y |- 值则为西或者南的秒数
23-1 89-90 坐标单位；1=长度(米或者英尺)；2=弧度的秒。
23-2 91-92 风化层速度。
24-1 93-94 降速层速度。
24-2 95-96 震源处的井口时间。
25-1 97-98 接收点处的井口时间。
25-2 99-100 炮点的静校正。
26-1 101-102 接收点的静校正。
26-2 103-104 应用的总静校正量(如果没有应用静校正为零)。
27-1 105-106 延迟时间-A，以ms 表示。
240字节的道标识的结束和时间信号之间的时间。
如果时间信号出现在道头结束之前为正。
如果时间信号出现在道头结束之后为负。
时间信号就是起始脉冲，它记录在辅助道上或者由记录系统指定。
27-2 107-108 时间延迟-B，以ms 表示。
为时间信号和能量起爆之间的时间。
可正可负。
28-1 109-110 时间延迟时间，以ms 表示。
能量源的起爆时间和开始记录数据样点之间的时间(深水时，数据记录不从时间零开始。
)
28-2 111-112 起始切除时间。
29-1 113-114 结束切除时间。
29-2 115-116* 本道的采样点数。
30-1 117-118* 本道的采样间隔，以ms 表示。
30-2 119-120 野外仪器的增益类型：
1=固定增益；2=二进制增益；3=浮点增益；4…N=选择使用。
31-1 121-122 仪器增益常数。
31-2 123-124 仪器起始增益(DB)。
32-1 125-126 相关码：1=没有相关；2=相关。
32-2 127-128 起始扫描频率(HZ)。
33-1 129-130 结束扫描频率(HZ)。
33-2 131-132 扫描长度，以ms 表示。
34-1 133-134 扫描类型：1=线性；2=抛物线；3=指数；4=其他34-2 135-136 扫描道起始斜坡长度，以ms 表示。
35-1 137-138 扫描道终了斜坡长度，以ms 表示。
35-2 139-140 斜坡为型：1=线性；2= ；3=其他。
36-1 141-142 滤假频的频率(如果使用)。
36-2 143-144 滤假频的斜率。
37-1 145-146 陷波滤波器频率。
37-2 147-148 陷波斜率。
38-1 149-150 低截频率(如果使用)。
38-2 151-152 高截频率(如果使用)。
39-1 153-154 低截频率的斜率。
39-2 155-156 高截频率的斜率。
40-1 157-158 数据记录的年。
40-2 159-160 数据记录的日。
41-1 161-162 小时(24时制)。
41-2 163-164 分。
42-1 165-166 秒。
42-2 167-168 时间代码：1=当地时间；2=格式威治时间；3=其他43-1 169-170 道加权因子。
(最小有效位定义为2**(-n), n=0，1，2， (32767)
43-2 171-172 覆盖开关位置1的检波器道号。
44-1 173-174 在原始野外记录中道号1的检波点号。
44-2 175-176 在原始野外记录中最后一道的检波点号。
45-1 177-178 缺口大小(缺少的检波点总数)。
45-2 179-180 在测线的开始或者结束处的斜坡位置：
1=在后面；2=在前面。
46－1 181－182 数据道数。
46－2 183－184 未用。
47 185－188 炮点位置剩余静校正量。
48 189－192 检波点位置剩余静校正量。
49-1 193-194 应用后的CMP基准面静校正量。
49-2 195-196 检波点位置站号。
50 197-198 三维线号。
199-200 未用。
51 201-204 CMP点X坐标。
52 205-208 CMP点Y坐标。
53 209-212 统一坐标系投影方式。
54 213-216 标识号(IDENT)。
55-1 217-218 野外文件号。
55-2 219-220 炮点位置站号。
56-1 221-222 CMP位置高程。
56-2 223-224 CMP位置站号。
57-(60-1) 未用。
60-2 239-240 二维、三维识别符：0－二维；1－三维。
说明：带 * 的字节的信息必须记录。