基于SEG—Y格式地震勘探数据无失真压缩研究.docVIP

基于SEG—Y格式地震勘探数据无失真压缩研究 　　【摘 要】随着地震勘探的深入，地震勘探数据量也越来越大，那地震数据压缩是解决地震数据传输和存储问题的关键。在本文中研究了SEG-Y数据文件的特点，并分析了道头数据各道之间的统计规律，实现对道头数据无失真压缩。 　　【关键词】SEG-Y；统计规律；无失真；数据压缩 　　1.研究背景概述 　　由于地震勘探数据以及确定的深海的石油和天然气资源逐年增加，如何有效的传输和访问这些数据就是目前需要解决的问题。大量数据压缩是信息收集和处理技术的一个关键技术[1][2]。 　　2.SEG-Y格式结构 　　SEG（Society of Exploration Geop- 　　hysicists）在不同的时期发布了多种地震数据交换标准，在本文中主要用的是1975年版的Y格式（即SEG-Y rev 0）。 　　（ni为第i道中的样点数，一般情况下n1=n2=…=nm=n） 　　如图2-1所示，卷头数据块和地震道数据两大部分是SEG-Y数据文件主要的组成部分。这其中，卷头数据块共有3600字节，它又可以分为两部分别为3200字节的EBCDIC卡和400字节的二进制编码头部分；而m个地震道数据块（Trace data block）就组成地震道数据，而每个地震道数据块就有240字节的道头数据和多个样点数据。在本文我们根据样点数据的类型不同，有的每个样点数据由4字节（浮点型和定点4字节型）构成，有的由2字节（定点2字节型）构成[3][4]。 　　3.SEG-Y格式道头数据无失真压缩 　　前面分析了SEG-Y格式地震数据头段部分的特点，在此基础上将对道头数据无失真压缩做出具体分析并在VC++开发环境下用C语言编程实现。 　　3.1 道头数据的读取 　　通过对赵宇、王志良《石油工业计算机应用》在辽河取得的[5]多个SEG-Y地震数据文件进行分析，发现每一道头内的86个数据项之间存在着相关性，例如，第1～第4字节的测线中的道序号和第5～第8字节的卷中的道序号一般是相同的。相邻地震道的道头数据的对应数据项之间的强相关性是道头数据有冗余的主要原因。我们设SEG-Y数据文件中共记录m个地震道，道头数据中的86个数据项用Xi，1，Xi，2，……，Xi，86表示，其中i=1，2，……，m， 　　假设：X0，1=X0，2=……=X0，86=0和X-1，1=X-1，2=……=X-1，86=0，定义??头数据项的二阶差分为： 　　△2Xi，j=Xi，j-2Xi，j-1+Xi-2，j（i=1，2，……，m；j=1，2，……，86） （3-1） 　　通过对记录的10个SEG-Y数据文件[5]的道头数据进行分析，发现绝大多数道头数据项的二阶差分为零，并且从第3道至第m道的道头数据项的二阶差分全为零的数据项也很多，在86个数据项中平均只有7.9个数据项的道头数据项的二阶差分全为零。也就说二阶差分不全为零的数据项仅占7.9/86=9.18，这一发现是能有效地压缩道头数据的依据。 　　3.2 道头数据的压缩 　　通过前一节对读取出来的道头数据的分析发现，大部分数据项的二阶差分值都为0，因而可以只存储前2道数据和二阶差分值就可以将道头数据完全恢复。具体做法如下： 　　（1）首先需要计算道头的每个数据项二阶差分，按照（3-1）式的方法计算。对第i个数据项（i=1，2，…，86），判断Δ2x3，i=Δ2x4，i=……Δ2xm，i=0是否成立。若成立，则将标志位flag[i]置为0；否则，置flag[i]为1。对每个flag[i]（i=1，2，……，86）： 　　（a）若flag[i]为1，则将二阶差分值Δ2x3，i，Δ2x4，i……Δ2xm，i全部存储在另外一个文件中。 　　（b）若flag[i]为0，则不用存储二阶差分值到文件中。最后存储在文件中的是86个标志位flag[i]（i=1，2，……，86），前2道的道头数据和flag[i]为1时的二阶差分值Δ2x3，i，Δ2x4，i……Δ2xm，i。 　　对压缩倍数进行优化： 　　由于flag[i]（i=1，2，……，86）的值全是0和1，故可以将每8位用一个unsigned char型变量来表示，然后再存储到文件中，这样可以进一步减小文件的大小使压缩倍数进一步提高。 　　由于道头数据有可能超过65535，因而如要直接存储前2道数据，必须用4个字节来存储。为了使文件得到进一步的压缩，可以先计算前2道数据之间的一阶差分值，然后存储第1道数据和第2道数据与第1道数据之间的一阶差分值。而前2道数据之间的一阶差分值不会大于65535，因而一阶差分值只需要用2个字节来存储，即用short int类型的变量来表示。 　　（2）其次对于道头数据第一次压缩后产生的