Relaxation algorithms for matrix completion, with applications to seismic travel-time data interpolation

被引:0
作者
Baraldi, Robert [1 ]
Ulberg, Carl [2 ]
Kumar, Rajiv [3 ]
Creager, Kenneth [2 ]
Aravkin, Aleksandr [1 ]
机构
[1] Univ Washington, Dept Appl Math, Seattle, WA 98195 USA
[2] Univ Washington, Dept Earth & Space Sci, Seattle, WA 98195 USA
[3] DownUnder GeoSolut, Perth, WA, Australia
基金
美国国家科学基金会;
关键词
interpolation; constrained optimization; regional seismology; CONVERGENCE; RECOVERY;
D O I
10.1088/1361-6420/ab3204
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
Travel time tomography is used to infer the underlying three-dimensional wavespeed structure of the Earth by fitting seismic travel time data collected at surface stations. Data interpolation and denoising techniques are important pre-processing steps that use prior knowledge about the data, including parsimony in the frequency and wavelet domains, low-rank structure of matricizations, and local smoothness. We show how local smoothness structure can be combined with low rank constraints using level-set optimization formulations, and develop a new relaxation algorithm that can efficiently solve these joint problems. In the seismology setting, we use the approach to interpolate missing stations and de-noise observed stations. The new approach is competitive with alternative algorithms, and offers new functionality to interpolate observed data using both smoothness and low rank structure in the presence of data fitting constraints.
引用
收藏
页数:24
相关论文
共 29 条
[1]  
[Anonymous], 1997, NUMERICAL LINEAR ALG
[2]  
[Anonymous], 2000, Numerical Optimization
[3]   FAST METHODS FOR DENOISING MATRIX COMPLETION FORMULATIONS, WITH APPLICATIONS TO ROBUST SEISMIC DATA INTERPOLATION [J].
Aravkin, Aleksandr ;
Kumar, Rajiv ;
Mansour, Hassan ;
Recht, Ben ;
Herrmann, Felix J. .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2014, 36 (05) :S237-S266
[4]   Level-set methods for convex optimization [J].
Aravkin, Aleksandr Y. ;
Burke, James V. ;
Drusvyatskiy, Dmitry ;
Friedlander, Michael P. ;
Roy, Scott .
MATHEMATICAL PROGRAMMING, 2019, 174 (1-2) :359-390
[5]   A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems [J].
Beck, Amir ;
Teboulle, Marc .
SIAM JOURNAL ON IMAGING SCIENCES, 2009, 2 (01) :183-202
[6]   Near-optimal signal recovery from random projections: Universal encoding strategies? [J].
Candes, Emmanuel J. ;
Tao, Terence .
IEEE TRANSACTIONS ON INFORMATION THEORY, 2006, 52 (12) :5406-5425
[7]   Matrix Completion With Noise [J].
Candes, Emmanuel J. ;
Plan, Yaniv .
PROCEEDINGS OF THE IEEE, 2010, 98 (06) :925-936
[8]   Exact Matrix Completion via Convex Optimization [J].
Candes, Emmanuel J. ;
Recht, Benjamin .
FOUNDATIONS OF COMPUTATIONAL MATHEMATICS, 2009, 9 (06) :717-772
[9]  
Davis D, 2016, SCI COMPUT, P115, DOI 10.1007/978-3-319-41589-5_4
[10]   Simply denoise: Wavefield reconstruction via jittered undersampling [J].
Hennenfent, Gilles ;
Herrmann, Felix J. .
GEOPHYSICS, 2008, 73 (03) :V19-V28