Please wait a minute...
E-mail Alert Rss
 
物探与化探  2022, Vol. 46 Issue (5): 1214-1224    DOI: 10.11720/wtyht.2022.1458
  方法研究·信息处理·仪器研制 本期目录 | 过刊浏览 | 高级检索 |
基于分段采样和MCA的地震数据重建
王德英1(), 张凯2, 李振春2, 张医奎3, 许鑫1
1.中国石油勘探开发研究院西北分院,甘肃 兰州 730030
2.中国石油大学(华东) 地球科学与技术学院,山东 青岛 266580
3.物华能源科技有限公司,陕西 西安 710067
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
全文: PDF(9370 KB)   HTML
输出: BibTeX | EndNote (RIS)      
摘要 

数据重建是地震资料处理中一项重要的前期工作。压缩感知(compress sensing, CS)已经在数据重建领域取得了很好的应用。CS的关键是采样的随机性,随机采样将常规欠采样引起的互相干假频转化为较低能量的不相干噪声。一方面,传统的随机采样方法缺乏对采样点的约束,导致产生过多的噪声干扰,分段随机采样可有效地控制采样点之间的距离。另一方面,单一的数学变换会导致信号的不完全稀疏表达,影响数据重建效果,形态分量分析(morphological component analysis, MCA)将信号分解成几个具有显著特征的成分以逼近数据复杂的内部结构。本文在MCA框架下找到了一个新的字典组合(Shearlet+DCT),并使用块坐标松弛(block coordinate relaxation,BCR)算法得到最优解,从而获得理想重构结果。对实际资料的实验表明,该方法在重建分段随机采样数据时具有较好效果。

服务
把本文推荐给朋友
加入引用管理器
E-mail Alert
RSS
作者相关文章
王德英
张凯
李振春
张医奎
许鑫
关键词 压缩感知分段随机采样形态分量分析数据重建    
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.

Key wordscompressed sensing    segmented random sampling    morphological component analysis    data reconstruction
收稿日期: 2021-11-25      修回日期: 2022-05-18      出版日期: 2022-10-20
ZTFLH:  P631.4  
基金资助:中国石油勘探与生产分公司科技项目“薄储层全频处理方法研究与目标精细刻画技术攻关试验”(2022KT1503)
作者简介: 王德英(1997-),男,2019年毕业于防灾科技学院(现应急管理大学)并获得地球物理学理学学士学位,同年进入中国石油大学(华东)就读硕士研究生,主要从事压缩感知在地震勘探中的应用、速度建模、全波形反演等方面的研究工作,现在中国石油勘探开发研究院西北分院从事高分辨率处理和机器学习等方面的研究工作。Email:498987895@qq.com
引用本文:   
王德英, 张凯, 李振春, 张医奎, 许鑫. 基于分段采样和MCA的地震数据重建[J]. 物探与化探, 2022, 46(5): 1214-1224.
WANG De-Ying, ZHANG Kai, LI Zhen-Chun, ZHANG Yi-Kui, XU Xin. Seismic data reconstruction based on segmented random sampling and MCA. Geophysical and Geochemical Exploration, 2022, 46(5): 1214-1224.
链接本文:  
https://www.wutanyuhuatan.com/CN/10.11720/wtyht.2022.1458      或      https://www.wutanyuhuatan.com/CN/Y2022/V46/I5/1214
Fig.1  分段随机采样示意
Fig.2  欠采样数据及其频谱
a—合成地震记录;b—合成地震记录频谱;c—规则欠采样数据;d—规则欠采样数据频谱;e—完全随机采样数据;f—完全随机采样数据频谱;g—分段随机采样数据;h—分段随机采样数据频谱
Fig.3  原始数据、采样矩阵和缺失数据
a—原始数据;b—随机采样矩阵;c—Jitter采样矩阵;d—分段随机采样矩阵;e—完全随机采样数据;f—Jitter采样数据;g—分段随机采样数据
Fig.4  DCT字典重建数据与差值剖面
a—随机采样的DCT重建结果;b—Jitter采样的DCT重建结果;c—分段随机采样的DCT重建结果;d—随机采样的DCT差值剖面;e—Jitter采样的DCT差值剖面;f—分段随机采样的DCT差值剖面
Fig.5  Shearlet字典重建数据与差值剖面
a—随机采样的Shearlet重建结果;b—Jitter采样的Shearlet重建结果;c—分段随机采样的Shearlet重建结果;d—随机采样的Shearlet差值剖面;e—Jitter采样的Shearlet差值剖面;f—分段随机采样的Shearlet差值剖面
Fig.6  字典组合重建结果及差值剖面
a—随机采样的Shearlet+DCT重建结果;b—Jitter采样的Shearlet+DCT重建结果;c—分段随机采样的Shearlet+DCT重建结果;d—随机采样的Shearlet+DCT差值剖面;e—Jitter采样的Shearlet+DCT差值剖面;f—分段随机采样的Shearlet+DCT差值剖面
缺失数据 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
Table 1  重建结果评价参数
Fig.7  拉格朗日乘子及迭代次数与重建信噪比
Fig.8  含噪合成地震数据欠采样
a—原始含噪数据;b—分段随机采样50%数据
Fig.9  三种重建结果
a—DCT重建结果;b—Shearlet重建结果;c—Shearlet+DCT重建结果
原始数据 欠采样数据 方法 重建数据 变换值
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
Table 2  重建结果评价参数
Fig.10  原始剖面
Fig.11  连续缺失11道Shearlet+DCT重建
a—连续缺失11道图像;b—重建结果;c—绝对误差
Fig.12  连续缺失22道Shearlet+DCT重建
a—连续缺失22道图像;b—重建结果;c—绝对误差
[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.
doi: 10.1016/S0019-9958(66)80002-5
[4] Hennenfent G, Herrmann F J. Simply denoise:Wavefield Reconstruction via jittered undersampling[J]. Geophysics, 2008, 73(3): V19-V28.
doi: 10.1190/1.2841038
[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.
doi: 10.1190/1.2799517
[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.
doi: 10.1016/j.jappgeo.2015.04.003
[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.
doi: 10.7623/syxb201404009
[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.
doi: 10.7623/syxb201404009
[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.
doi: 10.1190/geo2015-0108.1
[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.
doi: 10.1190/1.2840373
[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] 杨飞龙, 俞岱, 孙渊. 基于压缩感知的复杂地表菲涅尔束偏移[J]. 物探与化探, 2019, 43(6): 1285-1290.
[2] 徐卫, 张华, 张落毅. 基于复值曲波变换的地震数据重建方法[J]. 物探与化探, 2016, 40(4): 750-756.
[3] 路交通, 曹思远, 董建华, 张. 基于稀疏变换的地震数据重构方法[J]. 物探与化探, 2013, 37(1): 175-179.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
京ICP备05055290号-3
版权所有 © 2021《物探与化探》编辑部
通讯地址:北京市学院路29号航遥中心 邮编:100083
电话:010-62060192;62060193 E-mail:whtbjb@sina.com