一种地震数据压缩方法、压缩存储结构及随机访问方法
【技术领域】
[0001]本发明涉及一种石油勘探数据处理技术，尤其是一种支持随机读取的地震数据压缩方法、压缩存储结构及随机访问方法。
【背景技术】
[0002]地震勘探已经进入全方位、宽频带、数字化的时代。随着地震勘探向多维度，多分量、高精度方向发展，勘探数据呈爆炸性增长。尤其是海量勘探地震数据的出现，给现代计算机存储和处理技术带来严峻的挑战。由于海量勘探地震数据的存储和处理都需要占用极大的硬件资源，因此如何对地震数据进行压缩以解决海量地震数据的存储问题，一直是业内技术人员的研究课题。
[0003]现有的地震数据压缩方法主要分为无损压缩和有损压缩两大类型。无损数据压缩方法主要去除原始数据中的冗余，利用尽可能少的码字来描述原始数据。有损数据压缩方法主要是利用某种正交变换(如小波变换，离散余弦变换等)将原始数据变换到频率域数据，并对频率域数据进行量化，然后对量化后的数据进行熵编码。
[0004]目前，无论是无损数据压缩还是有损数据压缩，在对数据进行访问和处理时，都必须先经历数据的解压缩还原过程。S卩，先将压缩后的数据全部还原并存储在磁盘上，然后才能对数据进行访问和处理。由于不能直接对压缩数据进行随机读取，因此势必会占用大量的存储空间，以及增加运行时间，降低工作效率。
【发明内容】
[0005]针对上述问题，本发明提出了一种新的可以支持随机访问的地震数据压缩方法、压缩存储结构，以及相应的随机访问方法。
[0006]本发明提供一种地震数据压缩方法，用于压缩包含若干地震道数据块的地震数据，所述压缩方法包括以下步骤:
[0007]道头无损压缩步骤，对每一个地震道数据块的道头数据去除冗余字段；
[0008]样点有损压缩步骤，对每一个地震道数据块求取地震道特征值以及各样点数据的特征向量值，利用地震道特征值以及各样点数据的特征向量值表示样点数据。
[0009]进一步地，上述道头无损压缩步骤包括:
[0010]利用一道头描述信息记录对道头数据中的各数据项的属性定义以及道头描述信息本身所占用的空间大小。
[0011]根据本发明的实施例，上述数据项的属性包括字面意义、数据类型和字节长度。
[0012]进一步地，上述样点有损压缩步骤包括:
[0013]读取一地震道数据块中的样点数据；
[0014]计算相邻两个样点数据之间的差值，取其中绝对值最大的差值记作样点最大差的绝对值；
[0015]根据样点最大差的绝对值和指定的压缩级别，按照下式求取该地震道数据块的地震道特征值Bin:
[0016]Bin = 2 X Dmax/Short_Leve I
[0017]式中，Dmax是样点最大差的绝对值，Short_Level是与压缩级别有关的参数；
[0018]根据地震道特征值，求取每一个样点数据的特征向量值；
[0019]依次输出压缩后的道头数据、第一个样点数据、样点最大差的绝对值以及其余样点数据的特征向量值。
[0020]根据本发明的实施例，通过以下步骤求取样点数据的特征向量值:
[0021]计算样点数据的浮点特征向量FTi = Di/Bin,
[0022]式中，Di为相邻两个样点之间的差值；
[0023]对FTi取整获得特征向量值Ti = Round (FTi)，
[0024]式中，Round函数是四舍五入取整函数。
[0025]根据本发明的实施例，上述样点数据可以是4字节浮点型数据，样点数据的特征向量值可以是2字节整形数据。
[0026]相应地，与压缩级别有关的参数Short_Level可以设置为65535。
[0027]此外，本发明还提供一种地震数据压缩存储结构，其特征在于，包括:
[0028]卷头；
[0029]文件头；
[0030]道头描述信息，用于记录对道头数据中各数据项的属性定义以及道头描述信息本身所占用的空间大小；
[0031]若干地震道数据块，每一地震道数据块中包括:
[0032]压缩的道头数据；
[0033]第一个样点数据；
[0034]样点最大差的绝对值；
[0035]其余样点数据的特征向量值。
[0036]根据本发明的实施例，上述数据项的属性包括字面意义、数据类型和字节长度。
[0037]本发明还提供一种对采用上述压缩存储结构的地震数据的随机访问方法，其包括以下步骤:
[0038]分别获取卷头、文件头和道头描述信息所占用的空间大小；
[0039]通过下式确定第M道地震道数据块的起始位置:
[0040]Start Pos = VheaderLength+FHeaderLength+DesLength+TRlengthX (M_l)。
[0041]式中，VheaderLength、FHeaderLength、DesLength以及 TRlength 分别是卷头、文件头、道头描述信息以及一个地震道数据块所占用的空间大小。
[0042]与现有技术相比，本发明利用地震数据的特点，分别对地震道头数据进行无损压缩和对地震样点数据进行有损压缩，去掉冗余的道头字段和压缩样点数据，通过此方法压缩的地震数据能够支持后期的随机定位，实时在内存中解压缩还原，能够直接为实际生产所用。
[0043]本发明的其它特征和优点将在随后的说明书中阐述，并且部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
【附图说明】
[0044]附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例共同用于解释本发明，并不构成对本发明的限制。在附图中:
[0045]图1是现有技术中地震数据以SEGY结构存储的示意图；
[0046]图2是本发明实施例中采用的地震数据压缩方法的工作流程图；
[0047]图3是本发明实施例中压缩后的地震数据存储结构的示意图；
[0048]图4是本发明实施例中道头描述信息的内容格式的示意图；
[0049]图5是本发明实施例中采用的地震数据解压缩方法的工作流程图。
【具体实施方式】
[0050]下面结合附图和实施例对本发明的技术方案进行详细地描述。
[0051]实施例一
[0052]图1是现有技术中地震数据以SEGY结构存储的示意图。SEGY是目前地震勘探中最常用的数据存储形式。一份SEGY数据文件通常包括卷头、文件头和地震道数据三大部分。按照目前惯例，卷头数据一般为3200字节，文件头一般为400字节，地震道数据分为若干地震道数据块(分别记为第I道?第M道)。每一个地震道数据块又分为240字节的道头数据以及若干样点数据(分别记为样点I?样点N)。在本实施例中，样点数据为4字节浮点型。
[0053]针对上述地震数据的特点，本发明提出一种新的地震数据压缩方法，对每一个地震道数据块中的道头数据和样点数据分别进行压缩，并将压缩后的地震数据以一种新的结构存储，以便后期不通过解压缩就能够读取地震数据。下面结合实施例详细地说明该压缩方法的工作原理和工作流程(参见图2)，以及压缩后的地震数据存储结构(参见图3)。
[0054]如图2所示，该压缩方法可以大致分为道头无损压缩和样点有损压缩两大步骤。
[0055]S110、对每一个地震道数据块的道头数据去除冗余字段。
[0056]一般而言，地震道数据块的道头数据中总有部分数据字段是空值或者无效值，其比例有时甚至高达70%。因此有必要去除道头数据中毫无意义的冗余字段，也即对道头数据进行无损压缩。在本实施例中，可以通过人工干预的方式保留有意义的数据项，并重新定义各数据项的属性。具体地，如图3所示，优选在文件头与地震道数据之间增加一道头描述信息，用于记录对道头数据中的各数据项的属性定义(如字面意义、数据类型和字节长度的定义)，以及记录道头描述信息本身所占用的空间大小(DesLength)。
[0057]图4显示了本实施例中道头描述信息(TraceHeaderDescript1n)的一段内容。其中，<HeaderItem>表示道头数据中的一个数据项，<descript1n>XX<descript1n>定义该数据项的字面意义，<type>XX<type>定义该数据项的数据类型，<length>XX<length>定义该数据项的字节长度。道头描述信息本身所占用的空间大小DesLength =这段内容的字节长度+4字节。这是因为本实施例中采用4个字节保存DesLength参数。
[0058]S120、对每一个地震道数据块求取地震道特征值以及各样点数据的特征向量值，利用地震道特征值以及各样点数据的特征向量值表示样点数据。
[0059]具体地，对于每一个地震道数据块都执行下面的步骤。
[0060]S120.1、读取该地震道数据块中的样点数据。
[0061]S120.2、计算相邻两个样点数据之间的差值，并将其中绝对值最大的差值记为样点最