CN107045142B - 基于压缩感知的小波域地震数据实时压缩与高精度重构方法 - Google Patents
基于压缩感知的小波域地震数据实时压缩与高精度重构方法 Download PDFInfo
- Publication number
- CN107045142B CN107045142B CN201710181390.0A CN201710181390A CN107045142B CN 107045142 B CN107045142 B CN 107045142B CN 201710181390 A CN201710181390 A CN 201710181390A CN 107045142 B CN107045142 B CN 107045142B
- Authority
- CN
- China
- Prior art keywords
- data
- algorithm
- compression
- wavelet
- calculation matrix
- 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.)
- Active
Links
- 238000007906 compression Methods 0.000 title claims abstract description 42
- 230000006835 compression Effects 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 21
- 239000011159 matrix material Substances 0.000 claims abstract description 25
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 238000010276 construction Methods 0.000 claims abstract description 6
- 238000005070 sampling Methods 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 8
- 238000013461 design Methods 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 6
- 238000012544 monitoring process Methods 0.000 claims description 3
- 238000013144 data compression Methods 0.000 abstract description 11
- 238000002474 experimental method Methods 0.000 abstract description 2
- 238000012545 processing Methods 0.000 description 7
- 230000008901 benefit Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 4
- 238000000638 solvent extraction Methods 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241000260732 Buteo regalis Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000006837 decompression Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M7/00—Conversion of a code where information is represented by a given sequence or number of digits to a code where the same, similar or subset of information is represented by a different sequence or number of digits
- H03M7/30—Compression; Expansion; Suppression of unnecessary data, e.g. redundancy reduction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
- G01V2210/23—Wavelet filtering
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Theoretical Computer Science (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种基于压缩感知的小波域地震数据实时压缩与高精度重构方法,包括以下步骤:首先对微震信号在小波域稀疏表示;然后利用Logistic混沌序列构造混沌伯努利测量矩阵(CBMM),并用测量矩阵对稀疏表示的微震信号压缩观测;最后采用贝叶斯小波树结构压缩感知重构算法(BTSWCS),恢复出完整原始数据。本发明经实际对比实验结果表明:使用本方法对总采样点为28的数据压缩,压缩时间可缩短至10‑5s,也就是说,若地震仪采样率是1KSPS,CBMM测量矩阵对它采集的0.25s的数据基本可以实现实时压缩。低信噪比情况下,本发明重构算法使PSNR值至少提升5dB,相对贪婪算法,明显提高了峰值信噪比,降低了重构误差。
Description
技术领域
本发明涉及一种微地震数据的处理方法,特别涉及一种基于压缩感知的小波域地震数据实时压缩与高精度重构方法。
背景技术
数据实时传输是制约无缆存储式地震仪发展的最主要因素。随着地震勘探的不断深入,地震勘探正向多维、多分量和高分辨等方向发展,使得地震勘探数据持续膨胀,给无缆地震仪的采样和无线传输速度、存储器的存储容量以及计算机的处理速度造成极大压力。通过对地震数据进行合理的压缩处理,可以提高实时处理速度,改善地震仪无线通讯数据传输性能。因此,缩减压缩时间、提高重构精度成为迫切需要研究的课题。
CN104749633A公开了一种《高效高精度地震勘探数据压缩方法》对地球物理勘探仪器中采集的信号经24位模数转换后的数据压缩1/3的数据量,把24位3字节数据压缩为16位2字节数据。首先,设24位模数转换器输出的数据为D24;接着,设经压缩后的数据格式为4位阶码E4+16位尾数M16;然后把24位数据在保留符号位的情况下转换为4位阶码E4加16位尾数M16,阶码E4的取值为0~8;最后把4个采样点数据组合成10个字节,即4个样点的阶码占用2个字节,4个尾数占用8个字节。该方法具有压缩误差小,能够保持数据原貌,实现高效高精度压缩的优点。
CN104378118A公开了一种《高效自适应地震数据流无损压缩及解压缩方法》,通过使用编码方式将采样数据原始的24位3字节形式自适应的压缩为1字节或2字节或3字节或者和少量数据转化为4字节;在数据压缩前,首先根据原始数据的数值大小进行所需字节的判定,0~63及-64~-1区间的数据压缩后为1字节,64~8191及-8192~-65区间的数据压缩后为2字节,8192~104875及-104876~-8193区间的数据压缩方法操作后为3字节,占用字节与原来相同,除以上数据范围外,其他24位有符号整型数据可以表示的其他整型数在经压缩算法操作后需要使用4个字节表示。该方法具有大量节省存储空间,显著提高数据传输效率的优点。
2011年石水勇公开了《基于提升小波变换的地震勘探数据有损压缩研究》,以小波变换为基础,研究了不同小波转化为提升方案的原理,用于地震勘探数据提升小波正逆变换。对基于小波变换的三种典型的有损压缩方法一EZW、SPIHT、EBCOT方法,研究了其各自的原理和优缺点,从而选择出SPIHT算法作为地震勘探数据的编码方法。对于地震勘探原始数据和重构数据,设计并使用5种不同的评价标准对重构数据进行分析。2015年徐锋涛公开了《基于EZW算法的地震数据压缩》,对地震信号进行二维小波变换;利用系数之间的相关性对变换后的系数采用EZW方法重新组织;采用自适应算术编码对量化后的系数进行无损熵编码。
虽然上述现有技术可用于微地震信号压缩重构,但由于地球物理数据的特殊性,数据采集时要求尽量保持数据原貌,但这些数据压缩方式均不能用于数据采集阶段,必须对完整的地震数据进行操作,而不能对一个或少量几个采样数据处理,无法对数据流实时压缩,并且解码纯粹是编码的逆运算,压缩倍数是以损失信号精度为代价,无法做到降低压缩比的同时,减少重构误差。
发明内容
本发明目的在于针对上述现有技术的不足,提供一种基于压缩感知的小波域地震数据实时压缩与高精度重构方法。
由压缩感知理论可知:若信号是稀疏的或者是在某个变换域内稀疏,那么信号可以由测量矩阵将原始信号变换为低维信号,再通过重构算法优化求解得到原始信号。本发明首先,在小波域对微震数据稀疏表示;然后,构造测量矩阵对稀疏系数边采集边压缩;最后,设计重构算法高精度恢复出完整数据。本发明与常规压缩编解码算法的区别在于:不是先采集后压缩,而是边采集变压缩;解码不受编码的制约,而是独立的重构算法恢复出采集的完整数据。
本发明的目的是通过以下技术方案实现的:
一种基于压缩感知的小波域地震数据实时压缩与高精度重构方法,包括以下步骤:
A、输入微地震监测实际数据x;
B、对输入数据用小波基Ψ稀疏表示,得到稀疏系数θ=Ψx;
C、构造混沌伯努利测量矩阵(CBMM)Φ对数据进行压缩;
D、利用测量矩阵对稀疏系数观测,得到压缩后的数据y=Φθ+n,n表示噪声;
E、设计贝叶斯小波树结构压缩感知重构算法(BTSWCS)恢复信号;
F、用重构算法求解完整数据的稀疏系数
G、对求得的稀疏系数反变换,得到完整的微震数据
H、与多级树集合分裂(Set Partitioning in Hierarchical Trees,SPIHT)算法,贪婪算法中的正交匹配追踪(Orthogonal Matching Pursuit,OMP)算法、广义匹配追踪(Generalized Orthogonal Matching Pursuit,GOMP)算法、压缩采样匹配追踪(Compressive Sampling Matching Pursuit,CoSaMP)算法作比较,分析试验结果,对测量矩阵和重构算法进行评价。
优选地,步骤B,所述小波基Ψ为bior4.4小波基。
与现有技术相比,本发明的有益效果在于:本发明基于压缩感知理论,采用构造测量矩阵的方式,对稀疏表示后的微震信号边采样边压缩,有效地解决了常规压缩方式需要先采集再压缩的问题。同时,通过设计重构算法对观测数据高精度恢复重构,且具有良好的抗噪性能。通过对比实验可以看出,本发明实现实时压缩,对低信噪比数据也能高精度重构,可以广泛地应用于地球物理勘探仪器和其他场合,有利于改善地震仪无线通讯数据传输性能、降低地球物理仪器野外部件的数据存储空间、提高数据传输效率的优点。
附图说明
图1为基于压缩感知的小波域地震数据实时压缩与高精度重构流程图;
图2为中国西北部某区的地震勘探实际数据图;
图3a为当数据长度为214时,分别用SPIHT编码算法和CBMM压缩所有时间柱状图;
图3b为当数据长度为212时,分别用SPIHT编码算法和CBMM压缩所有时间柱状图;
图3c为当数据长度为210时,分别用SPIHT编码算法和CBMM测量矩阵压缩所有时间柱状图;
图3d为当数据长度为28时,分别用SPIHT编码算法和CBMM测量矩阵压缩所有时间柱状图;
图4为BTSWCS算法精确重构所需的MCMC迭代次数曲线图;
图5a为压缩比0.2时,贪婪算法重构地震记录图;
图5b为压缩比0.2时,BTSWCS算法重构地震记录图;
图6a为实际数据小波变换产生的小波系数图;
图6b为BTSWCS算法重构出的小波系数图;
图7a1为实际数据第2道记录的时域波形图;
图7a2为BTSWCS算重构的实际数据第2道记录的时域波形图;
图7a3为实际数据第2道记录的频谱图;
图7a4为BTSWCS算重构的实际数据第2道记录的频谱图;
图7b1为实际数据第15道记录的时域波形图;
图7b2为BTSWCS算重构的实际数据第15道记录的时域波形图;
图7b3为实际数据第15道记录的频谱图;
图7b4为BTSWCS算重构的实际数据第15道记录的频谱图;
图8a为贪婪算法、SPIHT解码算法和BTSWCS算法重构PSNR折线图;
图8b为贪婪算法、SPIHT解码算法和BTSWCS算法重构RMSE折线图。
具体实施方式
下面结合附图对本发明进行详细说明。
如图1所示,一种基于压缩感知的小波域地震数据实时压缩与高精度重构方法,包括以下步骤:
A、输入微地震监测实际数据x;
B、对输入数据用小波基Ψ稀疏表示,得到稀疏系数θ=Ψx;
C、构造混沌伯努利测量矩阵(CBMM)Φ;
D、利用测量矩阵对稀疏系数观测,得到压缩后的数据y=Φθ+n,n表示噪声;
E、设计贝叶斯小波树结构压缩感知重构算法(BTSWCS);
F、用重构算法求解完整数据的稀疏系数
G、对求得的稀疏系数反变换,得到完整的微震数据
H、与多级树集合分裂(Set Partitioning in Hierarchical Trees,SPIHT)算法,贪婪算法中的正交匹配追踪(Orthogonal Matching Pursuit,OMP)算法、广义匹配追踪(Generalized Orthogonal Matching Pursuit,GOMP)算法、压缩采样匹配追踪(Compressive Sampling Matching Pursuit,CoSaMP)算法进行比较,分析试验结果,对测量矩阵和重构算法进行评价,评价标准是峰值信噪比(PSNR)和均方根误差(RMSE)。
上述步骤中有三个关键技术:稀疏表示、构造测量矩阵、设计重构算法,下面对它们的实现原理介绍。
1、本发明微震数据用小波基稀疏表示原理:
微震数据压缩方面小波分析是其理论基础,小波理论提供的自适应时频分析法能自动调节时-频窗,以适应实际分析的需要。由于地震信号通常表现出分频带性强、幅度变化大、在二维空间上相关性高等特性,且小波基与地震信号具有相似性,所以对地震信号进行小波变换后,得到的小波系数具有很好的稀疏性(可压缩性)。对一维信号(单道微震信号)的离散小波多分辨率分析,Mallet算法描述公式为:
式中,Cj,k、Dj,k分别是原始信号在尺度2-j下分解得到的低频分量和高频分量,h(m-2k)、g(m-2k)分别是低通和高通滤波器的系数。相应的重建公式如下:
对二维地震信号进行分解与重构可以分别对行和列进行一维小波变换得到,其中行号和列号分别对应地震剖面的时间和道号。根据小波基性质和算法仿真效果,算法选用“bior4.4”对地震信号稀疏表示。
2、本发明混沌伯努利测量矩阵(CBMM)的构造原理:
已知Logistic混沌系统如式(3)所示:
该Logistic系统在参数λ∈[1.872,2.0]区间的条件下,当初值u0=0.23,0.37,0.7时,其Lyapunov指数均大于0,此时的Logistic系统是混沌系统。
由该Logistic系统产生的随机序列{un}经过式(4)符号函数映射成序列{an},其中:
当λ=2.0时,由该Logistic系统产生的序列{an}满足Bernoulli分布,同时也满足RIP性质,则由Logistic系统产生的序列{an}可作为CS的测量矩阵。
本发明构造测量矩阵算法步骤为:
步骤1,根据Logistic映射扩频序列的相关分布,取λ=2.0,u0=0.37,通过该Logistic混沌系统来产生混沌序列{un},序列长度为M×N-1。
步骤2,将步骤1生成的混沌序列{un}通过式子(4)符号函数映射成序列{an}。
步骤3,将步骤2生成的序列{an}取N长截断形成M×N维测量矩阵Φ。
3、本发明所用贝叶斯小波树结构压缩感知重构算法(BTSWCS)设计原理:
为克服低信噪比下,常用贪婪算法重构性能差的问题,本发明考虑的是贝叶斯框架下的压缩感知恢复算法,利用参数的先验概率和样本信息所满足的后验概率,继而获得总体的概率密度函数。
BTSWCS算法步骤为:
步骤1,假设贝叶斯先验模型,设定稀疏系数θ中每个元素都服从钉铆先验分布
步骤2,对先验模型参数后验估计,本发明用MCMC推理,该方法是通过足够的样本数来估计后验分布,通过迭代每个随机变量(模型参数和中间变量)来收集这些样本,迭代的随机变量是从它的条件后验分布给出的所有其他随机变量的最新值。
p(αn|-)=Gamma(a0+N/2,b0+(y-Φθ)T(y-Φθ)/2) (14)
本发明实验所用实际数据如图2所示,为证明用CBMM算法压缩要比常用SPIHT编码算法更具时效性,选用长度分别为214、212、210和28的数据做压缩比为0.2、0.25、0.3、0.35、0.4、0.45、0.5和0.55的处理,记录CBMM算法和SPIHT编码算法完成压缩所用时间。
因为测量矩阵是由随机数组成的,所以结果取50次试验的平均值,如图3a、3b、3c和3d所示。从柱状对比图可以清晰看出:压缩相同数据,达到相同压缩比情况下,CBMM算法要比SPIHT编码算法所用时间更短。特别是,随着待压缩数据量的减少,CBMM压缩的优势更加突出,对长度为28的数据用CBMM矩阵观测,所用压缩时间更是缩短至10-5s级;也就是说,若地震仪采样率是1kHz,CBMM测量矩阵对它采集的0.25s的数据基本可以实现实时压缩。
下位机通过CBMM对地震数据边采集变压缩,上位机通过BTSWCS重构算法恢复出完整记录。本发明实验选用的实际数据(图2)总长214,在压缩比分别为0.2、0.25、0.3、0.35、0.4、0.45、0.5和0.55的情况下,BTSWCS算法精确重构所需的MCMC迭代次数曲线如图4所示。为了验证BTSWCS重构算法在低信噪比下的优越性,首先对图2所示的实际数据用CBMM进行压缩比为0.2的压缩,然后用贪婪算法恢复(OMP),结果如图5a所示。实际数据在压缩比为0.2的情况下,根据图4,BTSWCS算法200次迭代用重构出完整数据,结果如图5b所示。可以很直观的看出,用BTSWCS算法重构效果明显比贪婪算法好,婪算法对本身含有噪声和含有测量噪声的数据进行重构处理时,重构性能明显下降,不能达到重构精度要求。图6和图7是对图5b重构结果的细节进行评价。图6a是实际数据小波系数,图6b是重构出的小波系数;图7a1是实际数据第2道记录的时域波形,图7a2是重构的第2道记录的时域波形,图7a3是实际数据第2道记录的频谱图,图7a4是重构的第2道记录的频谱图;图7b1是实际数据第15道记录的时域波形,图7b2是重构的第15道记录的时域波形,图7b3是实际数据第15道记录的频谱图,图7b4是重构的第15道记录的频谱图。
用贪婪算法、SPIHT解码和BTSWCS算法对压缩的实际数据恢复重构,结果取30次实验平均值,三种算法重构性能比较如图8所示。图8a是三类算法重构实际数据时的PSNR对比曲线,对实际数据的处理,贪婪算法整体峰值信噪比小于65dB,SPIHT解码和BTSWCS算法重构的PSNR值均大于贪婪算法重构。尤其是BTSWCS算法,它的PSNR至少比贪婪算法高出5dB,随着压缩比或迭代次数的增加,PSNR值会更高。图8b是三类算法重构实际数据时的RMSE对比曲线,贪婪算法在重构实际数据时,重构误差均大于0.1,重构误差很大,SPIHT解码和BTSWCS算法要比贪婪算法重构误差小很多,特别是BTSWCS,随着压缩比或迭代次数的增加,误差基本趋近于0。
Claims (2)
1.一种基于压缩感知的小波域地震数据实时压缩与高精度重构方法,其特征在于,包括以下步骤:
A、输入微地震监测实际数据x;
B、对输入数据用小波基Ψ稀疏表示,得到稀疏系数θ=Ψx;
C、构造混沌伯努利测量矩阵(CBMM)Φ对数据进行压缩;
D、利用测量矩阵对稀疏系数观测,得到压缩后的数据y=Φθ+n,n表示噪声;
E、设计贝叶斯小波树结构压缩感知重构算法(BTSWCS)恢复信号;
F、用重构算法求解完整数据的稀疏系数
G、对求得的稀疏系数反变换,得到完整的微震数据
H、与多级树集合分裂算法、贪婪算法中的正交匹配追踪算法、广义匹配追踪算法、压缩采样匹配追踪算法作比较,分析试验结果,对测量矩阵和重构算法进行评价。
2.根据权利要求1所述的一种基于压缩感知的小波域地震数据实时压缩与高精度重构方法,其特征在于:步骤B,所述小波基Ψ为bior4.4小波基。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710181390.0A CN107045142B (zh) | 2017-07-04 | 2017-07-04 | 基于压缩感知的小波域地震数据实时压缩与高精度重构方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710181390.0A CN107045142B (zh) | 2017-07-04 | 2017-07-04 | 基于压缩感知的小波域地震数据实时压缩与高精度重构方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107045142A CN107045142A (zh) | 2017-08-15 |
CN107045142B true CN107045142B (zh) | 2019-04-30 |
Family
ID=59545001
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710181390.0A Active CN107045142B (zh) | 2017-07-04 | 2017-07-04 | 基于压缩感知的小波域地震数据实时压缩与高精度重构方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107045142B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107561979B (zh) * | 2017-08-23 | 2020-02-11 | 厦门大学 | 一种面向冲击监测的数字异步压缩采样系统及方法 |
CN108076508B (zh) * | 2017-12-21 | 2021-03-26 | 儒安物联科技集团有限公司 | 从时空双维度扩展无线传感网的低功耗方法 |
CN107942377A (zh) * | 2018-01-05 | 2018-04-20 | 吉林大学 | 一种地震数据压缩与重构方法 |
CN110160642B (zh) * | 2019-05-13 | 2020-07-17 | 浙江大学 | 一种小样本条件下的螺旋桨噪声线谱重构方法 |
CN111329468A (zh) * | 2020-02-10 | 2020-06-26 | 郭洪军 | 一种心内科实时监测系统及其监测方法 |
CN111313988B (zh) * | 2020-02-22 | 2021-06-08 | 吉林大学 | 基于鞅变换模型实现urllc的免授权noma方法 |
CN112054806B (zh) * | 2020-09-15 | 2023-07-25 | 西南交通大学 | 一种基于二维自适应阈值小波的地铁塞拉门数据压缩方法 |
CN115379021A (zh) * | 2022-07-08 | 2022-11-22 | 金陵科技学院 | 基于椭圆曲线伪随机序列的煤矿微震数据压缩采集方法 |
CN115201904B (zh) * | 2022-07-18 | 2023-03-03 | 北京石油化工学院 | 基于边缘智能的微地震数据压缩和事件检测方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6895336B2 (en) * | 2002-09-12 | 2005-05-17 | Kelman Technologies, Inc. | Method of using matrix rank reduction to remove random noise from seismic data |
US7647182B2 (en) * | 2004-07-15 | 2010-01-12 | Baker Hughes Incorporated | Apparent dip angle calculation and image compression based on region of interest |
US8547789B2 (en) * | 2007-12-27 | 2013-10-01 | Schlumberger Technology Corporation | Method for gas zone detection using sonic wave attributes |
CN106199698A (zh) * | 2016-06-29 | 2016-12-07 | 中国石油化工股份有限公司 | 基于多次波信息的频率域地震数据重构方法 |
CN106371138B (zh) * | 2016-08-17 | 2018-10-16 | 中国石油天然气集团公司 | 地震数据重建方法和装置 |
CN106597541A (zh) * | 2017-02-22 | 2017-04-26 | 中国石油大学(华东) | 基于Shearlet变换的地震数据重构方法 |
-
2017
- 2017-07-04 CN CN201710181390.0A patent/CN107045142B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107045142A (zh) | 2017-08-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107045142B (zh) | 基于压缩感知的小波域地震数据实时压缩与高精度重构方法 | |
CN103595414B (zh) | 一种稀疏采样与信号压缩感知重构方法 | |
Goyal et al. | Quantized overcomplete expansions in ir/sup n: analysis, synthesis, and algorithms | |
CN105181122B (zh) | 机械振动信号数据压缩采集方法 | |
CN107192878A (zh) | 一种基于压缩感知的电力系统谐波检测方法及装置 | |
CN107942377A (zh) | 一种地震数据压缩与重构方法 | |
US6594394B1 (en) | Fast compression and transmission of seismic data | |
CN112134567B (zh) | 一种吸收光谱数据的自适应实时压缩方法与系统 | |
CN108989817A (zh) | 一种基于参考帧错位预测的雷达数据压缩方法 | |
Khelifati et al. | CORAD: correlation-aware compression of massive time series using sparse dictionary coding | |
CN103560795A (zh) | 一种心电信号的压缩方法、传输方法、重构方法及装置 | |
Donoho et al. | High-performance seismic trace compression | |
Kumar et al. | Computational efficient method for ECG signal compression based on modified SPIHT technique | |
Bai et al. | High-efficiency observations: compressive sensing and recovery of seismic waveform data | |
CN102291152B (zh) | 基于奇异值分解的压缩感知含噪信号重构系统 | |
US11327187B2 (en) | Functional quantization based data compression in seismic acquisition | |
CN103281087A (zh) | 基于多级压缩感知的信号采样系统及方法 | |
CN206727986U (zh) | 基于fpga的地震数据压缩装置 | |
CN113362405A (zh) | 一种基于StOMP压缩感知图像重构系统的构建方法 | |
Abdelwahab et al. | LiDAR data compression challenges and difficulties | |
Reddy et al. | A fast curvelet transform image compression algorithm using with modified SPIHT | |
Wakin | The geometry of low-dimensional signal models | |
Zheng et al. | A fast compression algorithm for seismic data from non-cable seismographs | |
CN103427845A (zh) | 基于二维分块dct变换的电力系统谐波数据压缩与重构方法 | |
Pascal et al. | Evaluation of EMG signals compression by JPEG 2000 called 1D |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |