地震数据压缩重构的正则化与零范数稀疏最优化方法  

Seismic data restoration based on compressive sensing using the regularization and zero-norm sparse optimization

作　　者：曹静杰[1,2,3] 王彦飞[2] 杨长春[2] 

机构地区：[1]石家庄经济学院,石家庄050031 [2]油气资源研究重点实验室,中国科学院地质与地球物理研究所,北京100029 [3]中国科学院研究生院,北京100049

出　　处：《地球物理学报》2012年第2期596-607,共12页Chinese Journal of Geophysics

基　　金：国家自然科学基金项目(10871191;40974075);中国科学院知识创新工程重要方向性项目(KZCX2-YW-QN107);石家庄经济学院博士科研启动基金联合资助

摘　　要：地震数据重构问题是一个病态的反演问题.本文基于地震数据在curvelet域的稀疏性,将地震数据重构变为一个稀疏优化问题,构造0范数的逼近函数作为目标函数,提出了一种投影梯度求解算法.本文还运用最近提出的分段随机采样方式进行采样,该采样方式能够有效地控制采样间隔并且保持采样的随机性.地震数值模拟表明,基于0范数逼近的投影梯度法计算效率有明显的提高;分段随机采样方式比随机欠采样有更加稳定的重构结果.Seismic data restoration is an ill-posed inverse problem. Based on the sparseness of seismic data in the curvelet domain, this problem can be transformed into a sparse optimization problem. This paper proposes to use the approximation of zero-norm as the objective function and develop a projected gradient method to solve the corresponding minimization problem. We also employ a recently proposed piecewise random sampling method which can both control the sampling gap and keep the randomness of sampling. Numerical results show that the projectedgradient method can reduce the amount of computation greatly, and the restoration based on the piecewise random sampling are better than that of random sub-sampling.

关 键 词：波场重构 CURVELET变换 压缩传感 0范数逼近 反问题 不适定性 稀疏优化 

分 类 号：P631[天文地球—地质矿产勘探]