|
|
Seismic data reconstruction based on segmented random sampling and MCA |
WANG De-Ying1( ), ZHANG Kai2, LI Zhen-Chun2, ZHANG Yi-Kui3, XU Xin1 |
1. Research Institute of Petroleum Exploration & Development-Northwest (NWGI),PetroChina,Lanzhou 730030,China 2. School of Geosciences,China University of Petroleum,Qingdao 266580,China 3. Wuhua Energy Technology Co.,Ltd,Xi'an 710067,China |
|
|
Abstract Data reconstruction is a critical preliminary work in the processing of seismic data.Compressed sensing (CS) has been well applied in data reconstruction.The key to CS is random sampling,which converts the mutual coherent alias caused by regular under-sampling into lower-amplitude incoherent noise. But traditional sampling methods lack constraints on sampling points, resulting in excessive noise interference. The segmented random sampling (SRS) method can effectively control the distance between sampling points. Furthermore, a single mathematical transformation will lead to incomplete sparse representation and impact data reconstruction. The morphological component analysis (MCA) can decompose a signal into several components with outstanding morphological features to approximate the complex internal structure of data. A new dictionary combination (Shearlet+DCT) has been found under the MCA framework, and the block coordinate relaxation (BCR) algorithm has been used to get the optimal solution to obtain desired reconstruction results. Tests of real data have proven that the proposed method can produce good effects when used to reconstruct the SRS data.
|
Received: 25 November 2021
Published: 03 January 2023
|
|
|
|
|
|
Schematic diagram of segmented random sampling
|
|
Sampled data and spectrum a—synthetic seismic record;b—spectrum of synthetic seismic record wavenumber;c—synthetic seismic data after regular sampling;d—spectrum of synthetic seismic data after regular sampling wavenumber;e—synthetic seismic data after random sampling;f—spectrum of synthetic seismic data after random sampling wavenumber;g—synthesized seismic record after SRS;h—spectrum of synthesized seismic record after SRS
|
|
Original data, sampling matrix and missing data a—original data;b—Random sampling matrix;c—Jitter sampling matrix;d—segmented random sampling matrix;e—missing data after random sampling;f—missing data after Jitter sampling;g—missing data after segmented random sampling
|
|
Reconstruction results and difference profiles of DCT dictionary a—DCT reconstruction result of random sampling;b—DCT reconstruction result of Jitter sampling;c—DCT reconstruction result of segmented random sampling;d—DCT reconstruction result of segmented random sampling;e—DCT difference profile of Jitter sampling;f—DCT difference profile of segmented random sampling
|
|
Reconstruction results and differences profile of Shearlet dictionary a—Shearlet reconstruction result of random sampling;b—Shearlet reconstruction result of Jitter sampling;c—Shearlet reconstruction result of segmented random sampling;d—Shearlet reconstruction result of segmented random sampling;e—Shearlet difference profile of Jitter sampling;f—Shearlet difference profile of segmented random sampling
|
|
Reconstruction results and difference profiles of Dictionary Combination a—Shearlet+DCT reconstruction result of random sampling;b—Shearlet+DCT reconstruction result of Jitter sampling;c—Shearlet+DCT reconstruction result of segmented random sampling;d—Shearlet+DCT reconstruction result of segmented random sampling;e—Shearlet+DCT difference profile of Jitter sampling;f—Shearlet+DCT difference profile of segmented random sampling
|
| | 缺失数据 | DCT | Shearlet | DCT+Shearlet | 完全随机采样 | SNR/dB | 2.9288 | -0.5554 | 0.5201 | 2.6251 | | PSNR/dB | 28.5367 | 24.5367 | 25.6122 | 27.7172 | Jitter随机采样 | SNR/dB | 3.3736 | -0.3296 | 7.2277 | 10.3023 | | PSNR/dB | 28.4657 | 24.7625 | 30.0910 | 31.3944 | 分段随机采样 | SNR/dB | 3.0158 | -0.5368 | 8.6044 | 12.0367 | | PSNR/dB | 28.1079 | 25.5553 | 30.6965 | 32.1288 |
|
Evaluation parameters of reconstruction results
|
|
Lagrangian multiplier and number of iterations and reconstruction signal-to-noise ratio
|
|
Under sampling of noise synthetic seismic data a—original noise data;b—50% missing data
|
|
Three reconstruction results a—DCT reconstruction result;b—Shearlet reconstruction result;c—Shearlet+DCT reconstruction result
|
原始数据 | 欠采样数据 | 方法 | 重建数据 | 变换值 | SNR=3.0102dB PSNR=19.8125dB | SNR=2.246dB PSNR=17.4515dB | DCT | SNR=13.1747dB | +10.1645 | PSNR=29.7691dB | Shearlet | SNR=15.6185dB | +12.6083 | PSNR=30.0910dB | Shearlet+DCT | SNR=17.1428dB | +14.1326 | PSNR=33.945dB |
|
Evaluation parameters of reconstruction results
|
|
Original profile
|
|
Shearlet+DCT reconstruction of 11 consecutive missing a—11 consecutive missing images;b—reconstruction result;c—absolute error
|
|
Shearlet+DCT reconstruction of 22 consecutive missing a—22 consecutive missing images;b—reconstruction result;c—absolute error
|
[1] |
张华, 陈小宏. 基于Jitter采样和曲波变换的三维地震数据重建[J]. 地球物理学报, 2013, 56(5):1637-1649.
|
[1] |
Zhang H, Chen X H. Seismic data reconstruction based on jittered sampling and curvelet transform[J]. Chinese J. Geophys., 2013, 56(5):1637-1649.
|
[2] |
唐刚. 基于压缩感知和稀疏表示的地震数据重建与去噪[D]. 北京: 清华大学, 2010.
|
[2] |
Tang G. Seismic data reconstruction and denoising based on compressive sensing and sparse representation[D]. Beijing: Tsinghua University, 2010.
|
[3] |
Leneman O. Random sampling of random processes:Impulse response[J]. Information and Control, 1966, 9(2):347-363.
|
[4] |
Hennenfent G, Herrmann F J. Simply denoise:Wavefield Reconstruction via jittered undersampling[J]. Geophysics, 2008, 73(3): V19-V28.
|
[5] |
Herrmann F J, Wang D, Hennenfent G, et al. Curvelet-based seismic data processing: A multiscale and nonlinear approach[J]. Geophysics, 2008, 73(1): A1-A5.
|
[6] |
Mosher C C. Generalized windowed transforms for seismic processing and imaging[C]// 82nd Annual International Meeting Expanded Abstracts,SEG, 2012.
|
[7] |
Yang P, Fomel S. Seislet-based morphological component analysis using scale-dependent exponential shrinkage[J]. Journal of Applied Geophysics, 2015, 118:66-74.
|
[8] |
Men Z, Ning H X, Zhang M G, et al. A method and application of irregular geometry design based on compressive sensing[C]// SEG Technical Program Expanded Abstracts, 2019.
|
[9] |
Sardy S, Bruce A G, Tseng P. Block coordinate relaxation methods for nonparametric wavelet denoising[J]. Journal of Computational and Graphical Statistics, 2000, 9(2):361-379.
|
[10] |
刘成明, 王德利, 王通, 等. 基于Shearlet变换的地震随机噪声压制[J]. 石油学报, 2014, 35(4):692-699.
|
[10] |
Liu C M, Wang D L, Wang T, et al. Random seismic noise attenuation based on the Shearlet transform[J]. Acta Petrolei Sinica, 2014, 35(4): 692-699.
|
[11] |
李海山, 吴国忱, 印兴耀. 形态分量分析在地震数据重建中的应用[J]. 石油地球物理勘探, 2012, 47(2):236-243.
|
[11] |
Li H S, Wu G C, Yin X Y. Morphological component analysis in seismic data reconstruction[J]. Oil Geophysical Prospecting, 2012, 47(2):236-243.
|
[12] |
周亚同, 刘志峰, 张志伟. 形态分量分析框架下基于DCT与曲波字典组合的地震信号重建[J]. 石油物探, 2015, 54(5):560-568.
|
[12] |
Zhou Y T, Liu Z F, Zhang Z W. Seismic signal reconstruction under the morphological component analysis framework combined with discrete cosine transform (DCT) and curvelet dictionary[J]. Geophysical Prospecting for Petroleum, 2015, 54(5):560-568.
|
[13] |
张良, 韩立国, 许德鑫, 等. 基于压缩感知技术的Shearlet变换重建地震数据[J]. 石油地球物理勘探, 2017, 52(2):220-225.
|
[13] |
Zhang L, Han L G, Xu D X, et al. Seismic data reconstruction with Shearlet transform based on compressed sensing technology[J]. Oil Geophysical Prospecting, 2017, 52(2): 220-225.
|
[14] |
徐卫, 张华, 张落毅. 基于复值曲波变换的地震数据重建方法[J]. 物探与化探, 2016, 40(4):750-756.
|
[14] |
Xu W, Zhang H, Zhang L Y. A study of seismic data reconstruction based on complex-valued curvelet transform[J]. Geophysical and Geochemical Exploration, 2016, 40(4):750-756.
|
[15] |
石战战, 夏艳晴, 周怀来, 等. 一种基于L1-L1范数稀疏表示的地震反演方法[J]. 物探与化探, 2019, 43(4):851-858.
|
[15] |
Shi Z Z, Xia Y Q, Zhou H L, et al. Seismic reflectivity inversion based on L1-L1-norm sparse representation[J]. Geophysical and Geochemical Exploration, 2019, 43(4): 851-858.
|
[16] |
孔旭, 密文天, 莫雄, 等. 基于MRAS证据权重法的湖南怀化地区金矿成矿预测[J]. 物探与化探, 2016, 40(3):467-474.
|
[16] |
Kong X, Mi W T, Mo X, et al. Metallogenic prediction of gold deposits with weighting of evidence based on MRAS in Huaihua area,Hunan Province[J]. Geophysical and Geochemical Exploration, 2016, 40(3): 467-474.
|
[17] |
Wu R S, Geng Y, Ye L. Preliminary study on Dreamlet based compressive sensing data recovery[C]// SEG Technical Program Expanded Abstracts, 2013.
|
[18] |
何真, 曹思远, 郝婳婕, 等. 基于自适应K-SVD的能量泄漏恢复研究[J]. 物探与化探, 2020, 44(2):362-371.
|
[18] |
He Z, Cao S Y, Hao H J, et al. Research on energy leakage recovery of adaptive K-SVD[J]. Geophysical and Geochemical Exploration, 2020, 44(2):362-371.
|
[19] |
张凯, 张医奎, 李振春, 等. MCA框架下Shearlet和DCT字典组合地震数据重建[J]. 石油地球物理勘探, 2019, 54(5):12.
|
[19] |
Zhang K, Zhang Y K, Li Z C, et al. Seismic data reconstruction method combined with Discrete Cosine Transform and Shearlet dictionary under Morphological Component Analysis framework[J]. Oil Geophysical Prospecting, 2019, 54(5):12.
|
[20] |
Li X, Guo M J, Li W H, et al. Sparsity promoting reconstruction with compressively acquired land data[C]// SEG Technical Program Expanded Abstracts, 2019.
|
[21] |
Zwartjes P, Gisolf A. Fourier reconstruction with sparse inversion[J]. Geophysical Prospecting, 2007, 5(2):199-221.
|
[22] |
郭奇, 曾昭发, 于晨霞, 等. 基于高精度字典学习算法的地震随机噪声压制[J]. 物探与化探, 2017, 41(5):907-913.
|
[22] |
Guo Q, Zeng Z F, Yu C X, et al. Seismic random noise suppression based on the high-precision dictionary learning algorithm[J]. Geophysical and Geochemical Exploration, 2017, 41(5):907-913.
|
[23] |
Kumar R, Wason H, Herrmann F J. Source separation for simultaneous towed-streamer marine acquisition—A compressed sensing approach[J]. Geophysics, 2015, 80(6):WD73-WD88.
|
[24] |
Mosher C C, Kaplan S T, Janiszewski F D. Non-uniform optimal sampling for seismic survey design[C]// 74th EAGE Conference and Exhibition, 2012.
|
[25] |
Neelamani R, Baumstein A, Gillard D, et al. Coherent and random noise attenuation using the curvelet transform[J]. The Leading Edge, 2008, 27:240-248.
|
[26] |
罗勇, 毛海波, 杨晓海, 等. 基于双重稀疏表示的地震资料随机噪声衰减方法[J]. 物探与化探, 2018, 42(3):608-615.
|
[26] |
Luo Y, Mao H B, Yang X H, et al. Seismic random seismic noise attenuation method on basis of the double sparse representation[J]. Geophysical and Geochemical Exploration, 2018, 42(3):608-615.
|
[1] |
Fei-Long YANG, Dai YU, Yuan SUN. Fresnel beam migration method based on compressed sensing under complex topographic conditions[J]. Geophysical and Geochemical Exploration, 2019, 43(6): 1285-1290. |
[2] |
Yong LUO, Hai-Bo MAO, Xiao-Hai YANG, Wen-Jie LI, Wen-Chao CHEN. Seismic random seismic noise attenuation method on basis of the double sparse representation[J]. Geophysical and Geochemical Exploration, 2018, 42(3): 608-615. |
|
|
|
|