Sparse-spike deconvolution based on adaptive step FISTA algorithm
PAN Shulin1, YAN Ke2, LI Lingyun3, JIANG Congyuan4, SHI Linguang3
1. School of Earth Science and Technology, Southwest Petroleum University, Chengdu, Sichuan 610500, China;
2. Southwest Oil & Gas Field Company, PetroChina, Chengdu, Sichuan 610051, China;
3. Research Institute of Geophysics, Shengli Oilfield Branch Co., SINOPEC, Dongying, Shandong 257000, China;
4. Department of Electrical and Electronic Engineering, Sichuan Vocational and Technical College, Suining, Sichuan 629000, China
Abstract:The optimal internal gradient step L can be searched by linear methods in the fast iterative shrinkage-thresholding algorithm (FISTA).But the step L can only be increased,which seriously affects the convergence of FISTA algorithm.In this paper,a sparse-spike deconvolution based on adaptive step FISTA algorithm is proposed.On the basis of FISTA algorithm,the step L is reduced appropriately before each iteration,and then the step L is found by a linear search to achieve the adaptive adjustment of L.In order to realize the algorithm theoretical converge,the auxiliary sequence of the conventional FISTA algorithm is modified by combining two consecutive L.Finally,the algorithm is converged.Theoretical model and field data tests show that the proposed method has better convergence and better anti-noise ability than the conventional FISTA algorithm,and obtain ideal inversion results of different signal-to-noise ratio data.
宋维琪,吴彩端.利用压缩感知方法提高地震资料分辨率[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.
[3]
王宇,韩立国,周家雄,等.L1-L2范数联合约束稀疏脉冲反演的应用[J].地球科学——中国地质大学学报,2009,34(5):835-840.WANG Yu,HAN Liguo,ZHOU Jiaxiong,et al.A-pplication of combined norm constrained sparseness spike inverse[J].Earth Science:Journal of China University of Geosciences,2009,34(5):835-840.
[4]
曹静杰.基于广义高斯分布和非凸Lp范数正则化的地震稀疏盲反褶积[J].石油地球物理勘探,2016,51(3):428-433.CAO Jingjie.Seismic sparse blind deconvolution based on generalized Gaussian distribution and non-convex Lp norm regularization[J].Oil Geophysical Prospecting,2016,51(3):428-433.
[5]
Chen S S,Donoho D L,Saunders M A.Atomic Decomposition by Basis Pursuit[M].Society for Industrial and Applied Mathematics,1998.
[6]
Nesterov Y.A method of solving a convex progra-mming problem with convergence rate O(1/k2)[J].Soviet Mathematics Doklady,1983,27(2):372-376.
[7]
Nesterov Y.Introductory Lectures on Convex Optimization[M].Kluwer Academic Publishers,2014.
[8]
Nesterov Y.Gradient methods for minimizing com-posite objective function[J].Core Discussion Papers,2007,140(1):125-161.
[9]
Beck A,Teboulle M.A fast iterative shrinkage-thre-sholding algorithm for linear inverse problems[J].Siam Journal on Imaging Sciences,2009,2(1):183-202.
Becker S R,Candès E J,Grant M C.Templates for convex cone problems with applications to sparse signal recovery[J].Mathematical Programming Computation,2011,3(3):165-218.
[13]
Zhang Z,Saligrama V.Rapid:Rapidly accelerated proxi-mal gradient algorithms for convex minimization[C].IEEE International Conference on Acoustics,Speech and Signal Processing,2015,doi:10.1109/ICASSP.2015.7178681.
[14]
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(46):1235-1243.
[15]
Mehrotra S.On the implementation of a primal-dual interior point method[J].Siam Journal on Optimization,1992,2(4):575-601.
[16]
Parikh N,Boyd S.Proximal algorithms[J].Foundations & Trends in Optimization,2013,1(3):121-231.
[17]
Chambolle A,Dossal C.On the convergence of the iterates of the "fast iterative shrinkage thresholding algorithm"[J].Journal of Optimization Theory and Applications,2015,166(3):968-982.
[18]
Levitin E S,Polyak B T.Constrained minimization methods[J].USSR Computational Mathematics & Mathematical Physics,1966,6(5):1-50.
[19]
蔡志东,刘聪伟,王勇,等.井地联合地震数据反褶积[J].石油地球物理勘探,2017,52(1):8-12.CAI Zhidong,LIU Congwei,WANG Yong,et al.Seismic data deconvolution with VSP operator[J].Oil Geophysical Prospecting,2017,52(1):8-12.
[20]
吕铁良,王永刚,谢万学,等.稀疏脉冲反演技术在井间地震反演中的应用[J].石油物探,2007,46(1):58-63.LYU Tieliang,WANG Yonggang,XIE Wanxue,et al.Application of sparse pulse inversion in crosswell seismic data[J].Geophysical Prospecting for Petroleum,2007,46(1):58-63.