[发明专利]一种用于地震勘探仪的数据传输无损压缩算法

说明书

技术领域

本发明涉及采集大量波形并实时传输的仪器设备，特别是石油地震勘探仪的数据传输无损压缩算法。

背景技术

地震勘探仪是对地震波进行采集、传输并记录的仪器，可用于陆地、海洋和井中地震勘探。以陆地地震勘探仪为例其组成结构如图1所示，包括：1)中央处理单元，由主机、显示器、电源、磁盘、打印机等组成；2)交叉站单元M个，管理M条大线；3)采集站单元M*N，组成阵列，每个采集站可采集1～4个通道地震波数据。系统工作时，成千上万个采集通道的地震波数据由各个采集站向本条大线的交叉站发送，经过交叉站上传给中央处理单元。目前地震波采样均选用24位高精度AD，采样频率逐步提高，道间距逐渐缩小，所需传输的数据量巨大且呈增加态势。

地震勘探仪面临的问题：不管是放炮方式还是连续监测采集方式，所有的数据必须在规定时间内上传到中央处理单元进行分析和记录，以保证数据的正确性和安全性，否则会造成巨大损失。面对大数据量的实时数据传输，带宽的瓶颈直接制约着单条大线的带道能力和道间距、整个系统的带道能力，从而影响施工效率、勘探分辨能力等。压缩可以从传输的源头减少数据量，提升系统能力。

发明内容

本发明旨在根据地震勘探波形的特点，在每个采集站内部采用一种特殊方法对AD采样的24位原始数据进行无损压缩，从而减小所需上传的数据量，最终在中央处理单元进行解压缩，完整还原每个采样点的24位原始数据。

算法的提出和实现。通过对大量地震波形的研究，发现波形是连续变化的，而且采样频率越高所显示的波形连续性越好，在大部分范围内相邻三个采样点的幅值越接近于线性。三个采样点间的幅值越接近线性，预测值与真实值的差值ΔV越小，所占用的有效位数也越少。因此提出了用预测差值的方法对采样点24位原始数据进行替换。设一帧包含多个采样点，构造指示数据位数段放在帧的前部，用以指示每个ΔV数据的有效位数，这在解压缩时是必要的。算法实现分两个步骤完成：

步骤1计算预测差值方法中的差值ΔV，确定差值ΔV的有效位数，将差值ΔV作为上传数据。

步骤2编制和处理指示数据位数段。求取一帧中每个采样点差值ΔV的有效数据位数作为指示数据位数段的基础数据。先对指示数据位数段的基础数据采用游程编码处理，然后再用霍夫曼编码做进一步压缩处理，以缩短指示数据位数段的长度。

本发明在算法的实际效果方面，由于根据数据特点采用了多次压缩处理，压缩能力强于通用压缩软件WinRAR，一般地震勘探波形压缩率可达到50％(一般WinRAR压缩率为75％)，在井中波形的压缩率可达到30％(一般WinRAR压缩率为50％)，大大减少了所需上传的数据量，相当于提高了系统的整体带宽和带道能力；在算法的复杂度和时效性方面，需要的信息量少、节省内存、方法简单，对采集单元的微处理器任务量或实时性没有显著影响；另外不需要额外硬件支持，造价低廉，易于实现。

附图说明

图1为地震勘探仪组成框图；图2为预测差值示意图：横轴为采样时刻，纵轴为采样点的幅值；图3为指示数据位数信息图；图4为游程编码示意图。

具体实施方式

以下结合附图对本发明的原理和特征进行描述，所举实例仅用于解释本发明，并非用于限定本专利。

预测差值定义。预测值就是采用当前采样点之前相邻的两个采样点真实值对本采样点进行预测所求得的值。预测差值就是预测值与当前采样点的真实值之差，用ΔV表示。因此在一帧采样数据中除了前2个采样点上传原始值，其他采样点上传ΔV即可。

预测差值方法中的差值ΔV计算。如图2所示，在系统正常采集状态下，设已获得相邻3个采样点的时刻和幅值分别是(t1,a)、(t2,b)、(t3,c)。由于均为等间隔采样，因此采样间隔为t3-t2＝t2-t1。依据采样点的幅值v与时间t的线性方程式v＝kt+C进行线性预测，将(t1,a)、(t2,b)代入，可求得在t3时刻对应的预测值是c＇＝2b-a。预测差值ΔV＝c-c＇＝c-2b+a。在解压缩时将求得原始值c＝ΔV+c＇。在压缩计算过程中，预测值的获得并不需要复杂的运算，甚至用加减法和移位就可以完成，这也是应用线性预测的好处之一。