Windows环境下基于Qt的SEGY格式地震数据可视化.docVIP

W indows环境下基于 Qt的SEGY格式地震 W indows环境下基于 Qt的 SEGY格式地震数据可视化 张 婉 王 磊 (“地下信息探测技术与仪器 ”教育部重点实验室 中国地质大学 ·北京 100083) 摘 要 SEGY格式的地震数据以工作站存储方式存储 ,在非工作站环境下无法直接查看数据文件中 的道头信息及采样数据 。针对这个问题 ,笔者在 W indow s环境下以 Q t为开发平台实现了 SEGY格式 地震数据可视化 ,阐述了地震数据的组成结构以及相关信息 ,并给出了具体应用实例 。 关键词 Q t; SEGY格式 ;可视化 中最广泛使用的 SEGY格式数据 。 SEGY地震数据格 式是美国 SEG学会 (勘探地球物理学家协会 )制定并 推荐的一种交换格式标准 ,目的是为了使不同的地球 物理公司采集的地震数据能够共享和交流 ,节约计算 机资源 ,此种地震数据格式目前正广泛应用于石油勘 探工业界 。 SEGY格式的地震数据文件由多炮地震记录以道 序方式编排而成 ,每炮记录由文件头和各地震道数据 体组成 ,地震道数据体又分为道头和采样值两部分 。 标准的 SEGY格式包括 3200 字节的字符串头段 , 400 字节的二进制头段和 240字节的道头字 。其中 3200字 节的字符串头段与 400字节的二进制头段合起来统称 为地震数据文件的文件头 ,总共占用 3600字节 。地震 道数据从第 241个字节依次开始存放 (如图 1所示 ) 。 0 引 言 随着地震勘探技术的不断提高 ,野外采集的地震 勘探数据可以通过采集仪器直接转换成以道序方式 排列的 SEGY数据格式 。一般地震资料采集器使用 的都是 U n ix工作站 ,而且 SEGY格式的地震数据是 以二进制形式存储的 ,对于一般用户而言 ,不能直接 观察到文件中的道头信息和振幅数据值 。但在实际 工作中 ,对地震数据进行处理之前 ,地球物理人员是 非常希望了解有关道头信息的 ,因为运用软件进行数 据处理的第一步就是加载数据 ,数据加载的正确与否 关系到后续的一系列工作 ,查看相关的道头信息及原 始数据就是为了确保数据加载的正确性 。基于此目 的 ,作者在微机的 M ic ro soft W indow s环境下以 Q t为 开发平台 ,实现了对 SEGY格式地震数据的读取及可 视化 。 1 SEGY地震数据格式 图 1 SEGY带结构示意图 文中用于实现读取及可视化的数据是地震勘探 第一作者简介 张婉 ( 1982 - ) ,女 ,中国地质大学 (北京 )地球物理与信息技术学院硕士研究生 ,主要研究方向为地震勘探 。 换成 ASC II码 ,要达 换成 ASC II码 ,要达到供用户文本阅读的效果 ,还需将 得到的 ASC II码转换为字符串形式 。在此调用自定义 函数 : SEGYedbcdic2 string ( _SEGYebcdic ebcdic) ,函数 中使 用 Q t 下 的 文 本 流 类 QTextStream , 调 用 形 式 为 QTextStream ts ( str, IO _W riteOnly ) ,以写入的形式构 造一个通过内部设备操作 Unicode QString的 str文本 流 。由于 SEGY数据中各个头段所占用的字节数固定 不变 ,在对数据读取过程中 ,可以利用 QFile 类中的 readB lock函数 ,调用形式为 : Q _LONG QFile: : readB lock ( char3 p , Q _ULONG len ) ,从打开的文件中读取最多 len字节到 p 并且返回实际读取得字节数 ,实现了对数 据的批量处理 ,提高了数据的读取效率 。因此对后来 的各头段的读取 ,皆采取此种方法 。 3. 2 400字节二进制头段 400字节二进制头段 ,以 16 位或 32 位二进制数 记录了整个 SEGY数据文件的基本情况 ,如采样点 数 、采样率 、数据存放类型等重要信息 。原始的 SEGY 格式地震数据是以工作站中的存储格式存储的 ,它与 微机上的存放格式不同 ,这是由于各计算机生产厂商 之间技术相互保密造成的 。工作站中的二进制数据 是高字节在前 、低字节在后 ,即 B igEnd ian; 而微机中 的则是低字节在前 、高字节在后 ,即 L ittleEnd ian。所 以要在微机的 w indow s环境下实现对 SEGY格式的 数据进行读取 ,第一步就是将工作站上的 SEGY数据 格式转换成微机能识别的 SEGY数据格式 ,也就是在 读取过程中进行高低位字节的转换 。根据存储信息 所占字节数不同 ,给出了具体的实现代码 ,其中 i2 [ ] 为存储 16位数据的数组 ,