CN102522999B - 一种三维声波测井数据的实时压缩方法 - Google Patents
一种三维声波测井数据的实时压缩方法 Download PDFInfo
- Publication number
- CN102522999B CN102522999B CN201110408530.6A CN201110408530A CN102522999B CN 102522999 B CN102522999 B CN 102522999B CN 201110408530 A CN201110408530 A CN 201110408530A CN 102522999 B CN102522999 B CN 102522999B
- Authority
- CN
- China
- Prior art keywords
- data
- mapping
- enum
- compressed
- encoded
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Compression, Expansion, Code Conversion, And Decoders (AREA)
Abstract
本发明公开了一种三维声波测井数据的实时压缩方法,通过对原始声波数据进行拆分为高n1位保留数据和低n2位待压缩数据,对待压缩数据进行分割长度为m的待压缩数据;然后按组进行对待压缩数据进行DCT变换,得到映射数据,对映射数据进行量化,即数值较大的数据保留较多的位数,而数值较小的数据则保留较少的位数,实现数据的实时压缩,对量化后的数据进行编码;最后将保留数据数组中的高n1位保留数据与编码后的数据进行合并,得到压缩数据,从而减少数据存储量和传输量,提高测井效率。
Description
技术领域
本发明属于数据压缩技术领域,更为具体地讲,涉及一种三维声波测井数据的实时压缩方法。
背景技术
三维声波测井是在正交偶极声波测井技术基础上发展起来的新一代声波测井技术,其测量原理是利用目前所有的声波,即单极、偶极及斯通利波测量模式对各种频带的波形进行综合测量以获取地层的三维声波特性,即纵波时差、横波及斯通利波在井筒轴向、径向和周向的变化,对地层特性的方位性提供完整的描述。通过该技术可以清晰地对均质地层和非均质地层的各向异性及各向异性形成的各种机理进行分析。
图1是三维声波测井仪器结构示意图。
如图1所示,三维声波测井仪器采用三个单极发射换能器和两个正交的偶极发射换能器发射声波信号,采用十三组接收换能器R1~R13的接收阵列接收声波信号。其中,三个单极发射器分别是上单极发射器、下单极发射器和远单极发射器,上单极发射器和下单极发射器位于接收器阵列两端,而远单极发射器和两个X、Y正交偶极发射器位于仪器下部的较远处。三个单极发射器用于产生不同源距的纵波、横波和斯通利波,X、Y偶极发射器用于产生弯曲波。通过不同组合模式可获得所需要的多种模式波形数据。
图2是三维声波测井仪的信号采集示意图。
如图2所示,采集板在DSP的控制下进行信号采集,每组接收阵列包括8个接收传感器,共需采集104通道的声波信号。
三维声波测井仪器在工作过程中,需要在不同井眼深度进行多次信号采集,每次采集含多个采样周期,各个周期分别采集不同模式的波形数据,每个采样周期,每通道需采样500点数据,因此,测井数据量非常大。
由于三维声波测井数据量非常大,井下数据传输率受限制,给数据的实时存储和传输带来了困难。所以,必须选用合适的数据压缩算法,对测井数据进行实时压缩处理,减少数据量,提高测井效率
发明内容
本发明的目的在于克服三维声波测井数据大,在井下数据传输受限的情况下,实时存储和传输困难的缺陷,提供一种三维声波测井数据的实时压缩方法,以实现井下仪器在大量采集信号的同时对测井数据进行实时、有效地压缩,减少数据存储量和传输量,提高测井效率。
为实现上述发明目的,本发明三维声波测井数据的实时压缩方法,其特征在于,包括以下步骤:
(1)、原始声波数据的预处理
对输入缓冲区的n位原始声波数据拆分为高n1位保留数据和低n2位待压缩数据,并分别存储到保留数据数组和待压缩数据数组中;
从待压缩数组中读入待压缩数据,将其分割为一组一组的长度为m的待压缩数据;
(2)、待压缩数据的压缩
2.1)、DCT映射变换
按组对待压缩数据进行DCT正变换,将待压缩数据的能量尽量集中在少量的系数上,从而最大限度地去除压缩数据中的相关性,得到映射数据;
2.2)、动态量化和数据编码
采用动态量化的方式对映射数据进行量化,对数值较大的数据保留较多的位数,而对数值较小的数据则保留较少的位数;对量化后的数据进行编码;
(3)、数据合并与存储
将保留数据数组中的高n1位保留数据与编码后的数据进行合并,得到压缩数据,并存储到输出缓冲区。
本发明的发明目的是这样实现的:
本发明三维声波测井数据的实时压缩方法,通过对原始声波数据进行拆分为高n1位保留数据和低n2位待压缩数据,对待压缩数据进行分割长度为m的待压缩数据;然后按组进行对待压缩数据进行DCT变换,得到映射数据,对映射数据进行量化,即数值较大的数据保留较多的位数,而数值较小的数据则保留较少的位数,实现数据的实时压缩,对量化后的数据进行编码;最后将保留数据数组中的高n1位保留数据与编码后的数据进行合并,得到压缩数据,从而减少数据存储量和传输量,提高测井效率。
附图说明
图1是三维声波测井仪器结构示意图;
图2是三维声波测井仪的信号采集示意图;
图3是三维声波测井数据的实时压缩方法一具体实施方式流程图;
图4是原始声波数据的预处理示意图;
图5是待压缩数据的DCT映射变换、映射数据的量化、编码示意图;
图6是量化数据的编码示意图;
图7是数据合并与存储示意图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述,以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是,在以下的描述中,当已知功能和设计的详细描述也许会淡化本发明的主要内容时,这些描述在这里将被忽略。
实施例
图3是三维声波测井数据的实时压缩方法一具体实施方式流程图。
1、预处理
在本实施例中,如图3所示,首先是对原始声波数据进行预处理,其预处理的过程如图4所示,在本实施例中,采样得到的有效声波数据为24位,为便于处理,在输入缓冲区中保存为32位的原始声波数据,其最高8位全为零。如图4所示,将输入缓冲区中的n=32位的原始声波数据拆分为高n1=16位保留数据和低n2=16位待压缩数据,并分别存储到保留数据数组和待压缩数据数组中。在后续操作中,只对所有待压缩数组中的数据,即原始声波数据的低16位数据进行压缩,待压缩完毕后,再将保留数据和压缩数据整理为完整的输出数据。
将待压缩数据的长度进行补零使之成为m=8的整数倍,然后分割为一组一组的长度为m=8的待压缩数据。
在实际的应用中,采集数据长度为500,则需要将长度500的待压缩数据进行补零修改为8的整数倍,即504,然后进行分割得到63组待压缩数据。
2、数据压缩
待压缩数据的DCT映射变换、映射数据的量化、编码示意图如图5所示。
2.1、DCT映射变换
读入一组长度为8的待压缩数据,进行一维DCT映射变换,变换公式为:Y=CX,其中,X=(x1,...,x8)T是由8个待压缩数据构成的8×1的列向量,是8×8的变换矩阵,Y是由变换后得到的8个映射数据构成的8×1的列向量。
在本实施例中,取N=8,此时的C称为余弦基本函数。我们可以根据C的定义式计算出这个函数的所有值,得到64个浮点系数,从而实现DCT映射变换。矩阵C是归一化正交矩阵,因此DCT是正交变换,如果输入数据的长度(实际应用中为504)大于N,那么可以将数据以N为单位分段,各段数据都应用C进行变换也是不相关的,也就是说,由余弦基本函数计算得到的64个浮点系数可用于所有输入数据的映射变换,在整个处理过程中无需更换变换矩阵,这一点可以大大加速计算过程。
在本实施例中,采用德州仪器公司(TI)的型号为TMS320F2812的DSP芯片来实现数据压缩。事先将计算好的变换矩阵的64个浮点系数转换为Q15格式(TI公司DSP的定点数据格式)的定点系数,保存在TMS320F2812的系数缓冲区中,这样可以避免浮点运算,而且能通过查表的方式简化操作,节约处理时间。此外,在进行DCT映射变换时,使用汇编语言进行程序设计,以乘累加指令MAC为核心,采用循环寻址方式访问输入缓冲区、系数缓冲区和输出缓冲区,这样可以在一个指令周期内完成一次乘加运算,充分发挥DSP的性能,快速实现向量点积运算,完成DCT映射变换。同时,在具体实施过程中,将实现DCT映射变换的代码搬移到TMS320F2812的内部RAM中运行,提高程序运行速度,满足三维声波测井的实时数据处理要求。
2.2、动态量化和数据编码
DCT映射变换本身不能实现数据压缩,数据压缩是由量化和编码实现的。为了说明量化和编码的实现过程,假设有如下的一组经DCT映射变换后得到的映射数据,如表1所示。映射数据用16位的无符号数表示,其符号位单独保存在符号位数组中。从表中可以看到,该组映射数据的上半部分,即前四个映射数据的数值较大,而下半部分,即后四个映射数据数值较小。
在本实施例中,采用动态量化的方式进行对映射数据进行量化:每次量化前检测每组映射数据的最大值,根据该数值确定该组数据的有效位,按保留位数从最高有效位开始保留,得到量化数据,并记录下量化数据最低位在映射数据中的位数,即量化位数,这样便于解压缩数据时,数据的恢复。
检测每组向量的最大值,根据最大值的有效位数确定量化的程度,有效位数越多则量化程度越低,反之量化程度越高。这样能最大程度地保存原始数据的信息。
由此,每组映射数据的量化位数随该组映射数据的有效位的位数动态变化。为了最大程度地保留映射数据的信息,减小压缩失真度,在本实施例中,量化时,每组映射数据的前四个映射数据保留有效位的最高7位,后四个映射数据保留有效位的最高6位,对于舍弃部分的按照过半进位的原则进位。从表中可以看出,量化后的数据总共只有“4×7+4×6”位,数据量为原来映射数据的40.6%。
表1
为了恢复原始数据,需要将这些量化后的数据进行编码。编码规则为:第一个编码字节(8位)的低4位存放量化位数,即后四个映射数据的量化位数。在本实施例中,如表1、图6所示,该组映射数据中最大数据(0x194A)的有效位数为13(前3位均为0,是无效数据),对于前四个映射数据,保留前7位而舍弃后6位,对于后四个映射数据,保留前6位而舍弃后7位,所以该组数据的量化位数为7(0111B)。第一个编码字节的第5位保存第一个映射数据的符号位(0),第一个编码字节最高位3位和第二个编码字节的最低4位保存第一个量化数据(1100101B);第二个编码字节的第5位保存第二个映射数据的符号位(0),第二个编码字节最高位3位和第三个编码字节的最低4位保存第二个量化数据(0011110B),依次类推,将前四个量化数据进行保存;
第五个编码字节的第5位保存第五个映射数据的符号位,第五个编码字节最高3位和第六个编码字节的最低3位保存第五个量化数据(000001B);第六个编码字节的第4位保存第六个映射数据的符号位,第六个编码字节最高位4位和第七个编码字节的最低2位保存第六个量化数据(000001B);第七个编码字节的第3位保存第七个映射数据的符号位,第七个编码字节最高位5位和第八个编码字节的最低1位保存第七个量化数据(000000B);第八个编码字节的第2位保存第八个映射数据的符号位,第八个编码字节最高位六位保存第八个量化数据(000000B)。
具体如图6所示。按照此原则,编码后数据的码字长度为“8×8”位,由于编码前映射数据码字的长度为“8×16”位,编码实现了50%的固定压缩比。
在本实施例中,由于采用了动态量化的方式对DCT映射变换后的波形数据进行量化,有效降低了数据压缩的失真度,仅在原始波形幅度较大的地方有相对较大的误差,完全满足三维声波测井的数据失真度指标要求。
3、数据合并与存储
将保留数据数组中的高n1=8位保留数据与8位压缩数据(由原始声波数据的低16位数据压缩得到)组合为16位的字数据并保存到输出缓冲区中,如图7所示。由于输入缓冲区中的原始声波数据是32位的双字数据,最终实现的声波数据固定压缩比为50%。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (4)
1.一种三维声波测井数据的实时压缩方法,其特征在于,包括以下步骤:
(1)、原始声波数据的预处理
对输入缓冲区的n位原始声波数据拆分为高n1位保留数据和低n2位待压缩数据,并分别存储到保留数据数组和待压缩数据数组中;
从待压缩数据数组中读入待压缩数据,将其分割为一组一组的长度为m的待压缩数据;
(2)、待压缩数据的压缩
2.1)、DCT映射变换
按组对待压缩数据进行DCT正变换,将待压缩数据的能量尽量集中在少量的系数上,从而最大限度地去除压缩数据中的相关性,得到映射数据;
2.2)、动态量化和数据编码
采用动态量化的方式对映射数据进行量化,对数值较大的数据保留较多的位数,而对数值较小的数据则保留较少的位数;对量化后的数据进行编码;
(3)、数据合并与存储
将保留数据数组中的高n1位保留数据与编码后的数据进行合并,得到压缩数据,并存储到输出缓冲区。
2.根据权利要求1所述的三维声波测井数据的实时压缩方法,其特征在于,量化时,检测每组映射数据的最大值,根据该数值确定该组数据的有效位,按保留位数从最高有效位开始保留,舍弃部分的按照过半进位的原则进位,得到量化数据,并记录下量化数据最低位在映射数据中的位数,即量化位数。
3.根据权利要求2所述的三维声波测井数据的实时压缩方法,其特征在于,所述的原始声波数据为32位,拆分为高16位保留数据和低16位待压缩数据,分组长度m为8;
量化时,对于前四个映射数据,保留前7位而舍弃后6位,对于后四个映射数据,保留前6位而舍弃后7位,以后四个量化数据最低位在映射数据中的位数作为量化位数。
4.根据权利要求3所述的三维声波测井数据的实时压缩方法,其特征在于,所述的编码为:第一个编码字节的最低4位保存量化位数,第一个编码字节的第5位保存第一个映射数据的符号位,第一个编码字节最高位3位和第二个编码字节的最低4位保存第一个量化数据;第二个编码字节的第5位保存第二个映射数据的符号位,第二个编码字节最高位3位和第三个编码字节的最低4位保存第二个量化数据,依次类推,将前四个量化数据进行保存;
第五个编码字的第5位保存第五个映射数据的符号位,第五个编码字节最高3位和第六个编码字节的最低3位保存第五个量化数据;第六个编码字的第4位保存第六个映射数据的符号位,第六个编码字节最高位4位和第七个编码字节的最低2位保存第六个量化数据;第七个编码字的第3位保存第七个映射数据的符号位,第七个编码字节最高位5位和第八个编码字节的最低1位保存第七个量化数据;第八个编码字的第2位保存第八个映射数据的符号位,第八个编码字节最高位六位保存第八个量化数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110408530.6A CN102522999B (zh) | 2011-12-09 | 2011-12-09 | 一种三维声波测井数据的实时压缩方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110408530.6A CN102522999B (zh) | 2011-12-09 | 2011-12-09 | 一种三维声波测井数据的实时压缩方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102522999A CN102522999A (zh) | 2012-06-27 |
CN102522999B true CN102522999B (zh) | 2014-06-25 |
Family
ID=46293813
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110408530.6A Expired - Fee Related CN102522999B (zh) | 2011-12-09 | 2011-12-09 | 一种三维声波测井数据的实时压缩方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102522999B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103498658B (zh) * | 2013-08-27 | 2019-03-19 | 中国石油化工集团公司 | 一种伽马测井数据处理方法 |
CN103711474B (zh) * | 2013-12-19 | 2016-08-17 | 天津大学 | 一种正交偶极子声电组合测井仪器 |
CN103840839B (zh) * | 2014-03-21 | 2017-06-27 | 中国科学院声学研究所 | 一种井下声波成像测井数据的实时压缩方法 |
CN105227503B (zh) * | 2015-09-08 | 2019-01-18 | 北京航空航天大学 | 一种基于无线随钻测量***的井下信源信道联合编码方法 |
CN106527975A (zh) * | 2016-10-14 | 2017-03-22 | 惠州Tcl移动通信有限公司 | 一种基于移动终端的文件数据存储及还原处理方法及*** |
CN106772567A (zh) * | 2016-12-23 | 2017-05-31 | 中国石油大学(华东) | 一种用于地震勘探仪的数据传输无损压缩算法 |
CN106685764B (zh) * | 2016-12-30 | 2020-10-23 | 新控技术(北京)有限公司 | 一种基于数据采集传输的通讯***及通讯方法 |
CN106845534B (zh) * | 2017-01-05 | 2021-05-28 | 北京大学深圳研究生院 | 一种基于数据压缩表示的支持向量机快速实现方法和装置 |
CN108494408B (zh) * | 2018-03-14 | 2021-07-13 | 电子科技大学 | 基于哈希字典的随钻密度测井仪井下高速实时压缩方法 |
CN108494409B (zh) * | 2018-03-14 | 2021-07-13 | 电子科技大学 | 基于小字典的随钻中子测井仪井下高速实时压缩方法 |
CN109831209B (zh) * | 2018-12-29 | 2023-06-20 | 长江大学 | 一种声波数据记录方法及装置 |
CN109600202B (zh) * | 2019-01-14 | 2021-06-01 | 同济大学 | 一种基于网络的远程ais接口适配装置 |
CN109905131B (zh) * | 2019-02-27 | 2021-01-29 | 中国科学院半导体研究所 | 一种脉冲序列压缩编码方法和装置 |
CN110120852A (zh) * | 2019-05-09 | 2019-08-13 | 上海格联医疗科技有限公司 | 一种手持超声设备无线传输提高带宽的数据处理方法 |
CN114151066B (zh) * | 2021-10-09 | 2023-04-25 | 电子科技大学 | 一种超声Lamb波测井井壁声学界面逆时偏移成像方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7251376B2 (en) * | 2003-08-29 | 2007-07-31 | Canadian Space Agency | Data compression engines and real-time wideband compressor for multi-dimensional data |
US7433527B2 (en) * | 2004-02-26 | 2008-10-07 | Mitsubishi Denki Kabushiki Kaisha | Time series data dimensional compression apparatus |
-
2011
- 2011-12-09 CN CN201110408530.6A patent/CN102522999B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102522999A (zh) | 2012-06-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102522999B (zh) | 一种三维声波测井数据的实时压缩方法 | |
CN103256044B (zh) | 一种随钻声波信号处理装置 | |
CN105181122B (zh) | 机械振动信号数据压缩采集方法 | |
CN107045142B (zh) | 基于压缩感知的小波域地震数据实时压缩与高精度重构方法 | |
CA2817632C (en) | System and method for computations utilizing optimized earth model representations | |
CN107942377A (zh) | 一种地震数据压缩与重构方法 | |
US6567781B1 (en) | Method and apparatus for compressing audio data using a dynamical system having a multi-state dynamical rule set and associated transform basis function | |
CN103742131A (zh) | 随钻声波井下信号采集与处理***的时差实时提取方法 | |
CN108989817A (zh) | 一种基于参考帧错位预测的雷达数据压缩方法 | |
CN103840839B (zh) | 一种井下声波成像测井数据的实时压缩方法 | |
CN105206277A (zh) | 基于单比特压缩感知的语音压缩方法 | |
CN101290774B (zh) | 音频编码和解码*** | |
CN104077479B (zh) | 一种基于守恒迎风格式获取参量阵声场空间分布的方法 | |
CN206727986U (zh) | 基于fpga的地震数据压缩装置 | |
CN105354867A (zh) | 自适应冗余字典压缩感知的高光谱图像压缩算法研究 | |
CN103776907A (zh) | 基于sinc插值的超声相控阵接收信号精延时方法 | |
CN101025919A (zh) | 音频解码中的合成子带滤波方法和合成子带滤波器 | |
CN101350199A (zh) | 音频编码器及音频编码方法 | |
CN104104390A (zh) | 一种信号压缩方法和重构方法及相关装置和*** | |
CN103985100A (zh) | 一种基于自适应观测组合优化的分块压缩感知方法 | |
CN104008237A (zh) | 一种计算溢流环换能器阵声辐射特性的方法 | |
CN114885036B (zh) | 一种探地雷达数据的实时有损压缩方法及*** | |
CN103366752B (zh) | 生成用于编码信息信号的候选码矢的方法和设备 | |
CN106572093B (zh) | 一种无线传感器阵列数据压缩方法及*** | |
CN115329260A (zh) | 一种基于偏移对角矩阵的Transformer加速器 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140625 Termination date: 20161209 |
|
CF01 | Termination of patent right due to non-payment of annual fee |