BS模式SEG-Y格式地震数据读取与演示.docVIP

BS模式SEG-Y格式地震数据读取与演示 　　【摘要】本文详细介绍SEG-Y文件的组成，在对SEG-Y详细研究的基础上，研究B/S模式下，用户远程对服务器提供的SEG-Y资源文件进行读取与演示。研究主要基于.NET平台，使用ASP.NET动态网站开发技术，代码主要由C#编写，最终实现用户对SEG-Y的读取，完成地震剖面图的绘制。 　　【关键词】B/S；剖面图；远程演示；SEG-Y 　　中图分类号：P315 文献标识码：A 文章编号： 　　1绪论 　　SEG-Y地震数据是地球物理领域运用很广泛的一种数据存储格式，本文从理论和实践两个角度出发，通过设计和实现一个基于.NET的SEG-Y浏览系统，深入探讨和了解B/S模式结构，ASP.NET在动态网站建设方面的优势和SEG-Y文件的分析读取与演示。 　　2开发技术 　　本文主要使用ASP.NET动态网页开发技术，网页后台程序由C#语言编写，前台使用HTML+CSS实现，客户端动态数据验证使用js。 　　3SEG-Y文件分析 　　3.1SEG-Y文件组成 　　3.1.1组成简介 　　SEG-Y文件由文件头文件和数据体两部分组成。文件头记录了此文件的各项参数，为文件的进一步读取给出必要信息。数据体则是文件的主体部分，它记录了主要的地震数据，是绘图的依据。 　　3.1.2文件头文件介绍与分析 　　文件头总长度为3600字节，分两部分:第一部分为字符型文件头，文件头第二部分为二进制文件头。 　　3.1.3数据体部分介绍与分析 　　数据体由多个数据道组成,每道数据分两部分：道头和采样数据。 　　3.2SEG-Y文件处理及实现 　　3.2.1SEG-Y文件字符型文件头（Textual File Header）处理及实现 　　字符型文件头数据由EBCDIC代码组成，所以，必须将其转换为ASCII代码才可以正确显示。具体方法是建立两种代码关系对应表，每次读入一个字节EBCDIC代码，循环查找EBCDIC代码表，找出它在EBCDIC表中的位置，在返回ASCII代码表中对应位置的ASCII代码。 　　3.2.2SEG-Y文件二进制文件头（Binary File Header）处理及实现 　　二进制文件头由16位和32位二进制数据构成，但在读取其中的各参数值时就要注意，SEG-Y文件以工作站的存储格式存储。工作站中的二进制数据是高字节在前，低字节在后面；而微机中的数据则是低字节在前，高字节在后。 　　我们采用先通过FileStream将400个字节的文件信息读入内存，400字节数据按从低地址到高地址的顺序存放在一个数组中，然后根据数据的存储位数，决定转换字节数，如读取16位数据，先读取第一个字节放入一个16位临时整型数中，将整型数据左移8位，然后，再读入下一字节，将其放入另一个16位临时整型数中，最后，将临时16位数按位或，得到16位微机格式的地震数据。 　　3.2.3SEG-Y文件道头文件（Trace Header）的处理及实现 　　道头与二进制文件头数据存储相同，所以，可以采用相同的数据处理方式，这里不再赘述。 　　3.2.4 SEG-Y文件道数据（Trace Data）的处理及实现 　　道数据格式与二进制头文件的控制位有关，长度由道采样点数及采样点格式共同决定。二进制头文件3225-3226字节标志了采样点数据（道数据）采用的数据格式，采样点数由道头115-116字节给出，若二进制文件头3503-3504字节标志为0，则该文件道长一致，可以根据二进制文件头3221-3222字节得出，无需分别读出每道采样点数。若采样点使用32位数据记录（如：IBM浮点型，IEEE浮点型，4字节整型），道长则为采样点数×4字节，若采样点使用16位数据记录（如：2字节整型），道长则为采样点数×2字节。 　　道数据依然是以工作站存储方式存储数据，所以，首先要进行高低位转换，然后根据数据类型，转换为微机存储格式以便读出。 　　SEG-Y道数据主要用IBM浮点型数据记录，后来也常用IEEE浮点型数据记录。这里主要介绍IBM浮点数的存储格式和转换方法。IBM浮点数分三部分：符号位(sign)占1 位, 指数部分(e)占7位, 尾数部分(m)占24位。 　　所得数值result = (-1)^sign * ( m / (2^24 ) )*16 ^ ( e^ - 64 )。 　　读取IBM浮点数具体方法为，先顺序读入道数据至数组，逐个读出4字节数据，每个数据先高低位换位（方法同二进制文件头处理方法）。将32位数右移31位，所剩为符号位，取出保存；将32位数与0x7f000000按位与，获得浮点数指数部分，保存；将32位数与0x00ffffff按位与，获得浮点数位数部分，保存