A strong shielding removal method of reflection coefficient domain based on compressed sensing with L2 norm constraint along layer
ZHANG Junhua1,2, WANG Jing3, WANG Yanguang4, LIU Libin4, LI Hongmei4, WANG Xi'an1,2
1. School of Geosciences, China University of Petroleum (East China), Qingdao, Shandong 266580, China; 2. Key labortory of Deep Oil and Gas, China University of Petroleum (East China), Qingdao, Shandong 266580, China; 3. Exploration and Development Research Institute, Shengli Oilfield Branch Company, SINOPEC, Dongying, Shandong 257015, China; 4. Geophysical Research Institute, Shengli Oilfield Branch Company, SINOPEC, Dongying, Shandong 257022, China
Abstract:Conventional compressed sensing methods are based on the sparse inversion of reflection coefficients to improve the resolution of seismic data. However, if the strong shielding layer information and the weak reservoir information are superimposed, the reservoir cannot be effectively predicted, because the information of the weak reflection layer will be lost when the strong shielding layer is removed. For this reason, taking advantage of the strong sparse representation ability of compressed sensing and the high resolution of the reflection coefficient domain without wavelet overlap, this paper proposes a strong shielding removal method of reflection coefficient domain based on compressed sensing with L2 norm constraint along the layer. This method is based on the theory of compressed sensing. First, according to the sparse characteristics of the reflection coefficient in the time domain, the strong shielding layer and the reservoir are separated using the information along the layer. Then the sparse inversion is performed, and finally the reflection coefficient after strong shielding removal is convolved with the original wavelet to obtain the high-resolution results in the absence of a strong shielding layer. The advantage of the method is that the high-resolution reflection coefficient can separate the information of the strong shielding layer and the reservoir, which is beneficial to extract and remove the strong shielding layer. The model tests and actual seismic data application show that the proposed method can accurately separate the reflection information of the strong shielding layer from the weak reflection information of the reservoir and thus improve the accuracy of reservoir identification. On the time-frequency slices after strong shielding layer removal, the weak energy of the reservoir can be seen, and low-frequency accompanying phenomena appear. On the energy half-time attribute slices along the layer after strong shielding removal, the attribute has a good correlation with the beach-bar sand reservoir, which can effectively identify favorable reservoir areas.
张军华, 王静, 王延光, 刘立彬, 李红梅, 王喜安. 基于压缩感知的反射系数域沿层L2范数约束去强屏蔽方法[J]. 石油地球物理勘探, 2022, 57(2): 405-413.
ZHANG Junhua, WANG Jing, WANG Yanguang, LIU Libin, LI Hongmei, WANG Xi'an. A strong shielding removal method of reflection coefficient domain based on compressed sensing with L2 norm constraint along layer. Oil Geophysical Prospecting, 2022, 57(2): 405-413.
赵铭海.地震相似背景分离技术在东营凹陷的应用[J].油气地质与采收率,2004,11(4):25-27.ZHAO Minghai.Application of separation technology in seismic similar background to Dongying sag[J].Petroleum Geology and Recovery Efficiency,2004,11(4):25-27.
[3]
李军亮.相似背景目标分离技术在滩坝砂储层预测中的应用——以济阳坳陷东营凹陷沙四段为例[J].工程地球物理学报,2008,5(3):332-336.LI Junliang.The application of separation technology in seismic semblance background to reservoir prediction of beach-bar sandbody:Case of Es4 sandbody in the area of Dongying sag of Jiyang depression[J].Chinese Journal of Engineering Geophysics,2008,5(3):332-336.
[4]
赵爽,李仲东,许红梅,等.多子波分解技术检测含煤砂岩储层[J].天然气工业,2007,27(9):44-47.ZHAO Shuang,LI Zhongdong,XU Hongmei,et al.Using multiple wavelet decomposition technique to detect the sandstone reservoir with coal layer[J].Na-tural Gas Industry,2007,27(9):44-47.
[5]
张军华,刘振,刘炳杨,等.强屏蔽层下弱反射储层特征分析及识别方法[J].特种油气藏,2012,19(1):23-26.ZHANG Junhua,LIU Zhen,LIU Bingyang,et al.Analysis and identification of reservoir characteristics of weak reflectors under strong shielding layer[J].Special Oil&Gas Reservoirs,2012,19(1):23-26.
[6]
向龙斌,李东安.基于多子波分解和旋回分析法去火成岩地震反射屏蔽[J].地质科技情报,2016,35(3):210-215.XIANG Longbin,LI Dong'an.Eliminating seismic shielding of igneous rock using wavelet decomposition and seismic cycle analysis[J].Geological Science and Technology Information,2016,35(3):210-215.
[7]
金成志,秦月霜.利用长、短旋回波形分析法去除地震强屏蔽[J].石油地球物理勘探,2017,52(5):1042-1048,1111.JIN Chengzhi and QIN Yueshuang.Seismic strong shield removal based on the long and short cycle ana-lysis[J].Oil Geophysical Prospecting,2017,52(5):1042-1048,1111.
[8]
江馀,张军华,韩宏伟,等.优化变分模态分解方法消除强反射影响——以东营凹陷沙四段滩坝砂目标处理为例[J].石油地球物理勘探,2020,55(1):147-152,166.JIANG Yu,ZHANG Junhua,HAN Hongwei,et al.Elimination of strong reflection influence based on optimized variational mode decomposition method:a case study of the target processing of beach bar sand of Es4 in Dongying Sag[J].Oil Geophysical Prospecting,2020,55(1):147-152,166.
[9]
WANG Y H.Seismic time-frequency spectral decomposition by matching pursuit[J].Geophysics,2007,72(1):V13-V20.
[10]
WANG Y H.Multichannel matching pursuit for seismic trace decomposition[J].Geophysics,2010,75(4):V61-V66.
张在金,张军华,李军,等.煤系地层地震强反射剥离方法研究及低频伴影分析[J].石油地球物理勘探,2016,51(2):376-383.ZHANG Zaijin,ZHANG Junhua,LI Jun,et al.A method for stripping coal seam strong reflection and low-frequency shadow analysis[J].Oil Geophysical Prospecting,2016,51(2):376-383.
[13]
LEVY S,FULLAGAR P K.Reconstruction of a sparse spike train from a portion of its spectrum and application to high-resolution deconvolution[J].Geophysics,1981,46(9):1235-1243.
[14]
CHOPRA S,CASTAGNA J,PORTNIAGUINE O.Thin-bed reflectivity inversion[C].SEG Technical Program Expanded Abstracts,2006,25:2057-2061.
[15]
PURYEAR C I,CASTAGNA J P.An algorithm for calculation of bed thickness and reflection coefficients from amplitude spectrum[J].SEG Technical Program Expanded Abstracts,2006,25:1767-1770.
[16]
张锐,张奎.基于压缩感知理论的深度域薄层反演[C].SPG/SEG北京2016国际地球物理会议电子论文集,2016,528-530.ZHANG Rui,ZHANG Kui.Compressed sensing inversion for thin-bed resolution in depth domain[C].Proceedings of SPG/SEG Beijing 2016 International Geophysical Conference,Beijing,China,2016,528-530.
[17]
陈祖庆,王静波.基于压缩感知的稀疏脉冲反射系数谱反演方法研究[J].石油物探,2015,54(4):459-466.CHEN Zuqing,WANG Jingbo.A spectral inversion method of sparse-spike reflection coefficients based on compressed sensing[J].Geophysical Prospecting for Petroleum,2015,54(4):459-466.
[18]
顾雯,章雄,徐敏,等.强屏蔽下薄储层高精度预测研究——以松辽盆地三肇凹陷为例[J].石油物探,2017,56(3):439-448.GU Wen,ZHANG Xiong,XU Min,et al.High precision prediction of thin reservoir under strong shielding effect and its application:a case study from Sanzhao Depression,Songliao Basin[J].Geophysical Prospecting for Petroleum,2017,56(3):439-448.
[19]
宋维琪,吴彩端.利用压缩感知方法提高地震资料分辨率[J].石油地球物理勘探,2017,52(2):214-219.SONG Weiqi,WU Caiduan.Seismic data resolution improvement based on compressed sensing[J].Oil Geophysical Prospecting,2017,52(2):214-219.
[20]
夏红敏,刘兰锋,张显辉,等.地震数据谱反演压缩感知算法实现及应用[J].石油地球物理勘探,2021,56(2):295-301.XIA Hongwei,LIU Lanfeng,ZHANG Xianhui,et al.Implementation and application of compressed sensing algorithm for seismic spectrum inversion[J].Oil Geophysical Prospecting,2021,56(2):295-301.
[21]
DONOHO D L.Compressed sensing[J].IEEE Transactions on Information Theory,2006,52(4):1289-1306.
[22]
CANDES E J,ROMBERG J,TAO T.Robust uncertainty principles:exact signal reconstruction from highly incomplete frequency information[J].IEEE Transactions on Information Theory,2006,52(2):489-509.
[23]
BARANIUK R G.Compressive sensing[J].IEEE Signal Processing Magazine,2007,24(4):118-121.
[24]
CANDES E J,WAKIN M B.An introduction to compressive sampling[J].IEEE Signal Processing Magazine,2008,25(2):21-30.
[25]
BECK A,TEBOULLE M.A fast iterative shrinkage-thresholding method for linear inverse problem[J].SIAM Journal on Imaging Sciences,2009,2(1):183-202.