L1范数约束被动源数据稀疏反演一次波估计
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
对于被动源地震数据,运用常规的互相关算法得到的虚拟炮记录中,不仅含有一次波反射信息,还包括了表面相关多次波.然而,通过传统的被动源数据稀疏反演一次波估计(EPSI)方法,可以求得只含有一次波,不含表面相关多次波的虚拟炮记录.本文改进了传统的被动源数据稀疏反演一次波估计问题的求解方法,将被动源稀疏反演一次波估计求解问题转化为双凸L1范数约束的最优化求解问题,避免了在传统的稀疏反演一次波估计过程中用时窗防止反演陷入局部最优化的情况.在L1范数约束最优化的求解过程中,又结合了2DCurvelet变换和小波变换,在2DCurvelet-wavelet域中,数据变得更加稀疏,从而使求得的结果更加准确,成像质量得到了改善.通过简单模型和复杂模型,验证了本文提出方法的有效性.
This work has improved the original algorithm to estimate primaries by sparse inversion of passive-source seismic data through replacing the original algorithm to solve the convex optimization problem with L1-norm constraint.It avoids using a time-window to prevent the inversion from into local optimization situations when estimating primaries by sparse inversion.Moreover,during the solving of the optimization problem with L1-norm constraint,2DCurvelet transform and wavelet transform are used at the same time.In 2D Curvelet and wavelet domains,the data become more sparse,then the results obtained are more accurate and the quality of imaging is improved.First,the method of convex optimization problem with L1-norm constraint is introduced to solve the problem of estimating primaries by sparse inversion of passive-source seismic data,instead of the steepest descent method under L0-norm constraint.Second,2DCurvelet transformand wavelet transform are combined during the sparse inversion.In the 2D Curvelet-wavelet domain,the data become more sparse.Comparing with 3DCurvelet transform,the velocity of 2D Curvelet-wavelet transform is improved.Third,a simple model and a complex model are used to simulate the passive seismic data.The method of convex optimization problem with L1-norm constraint and that combined with 2D Curvelet transform and wavelet transform are used to estimate primaries from the passive seismic data,respectively.At last,comparison with the results obtained by the traditional LSQR algorithm illustrates that the method proposed is feasible and effective.The method of estimating primaries by sparse inversion can directly estimate primaries from the passive seismic data,and obtain the virtual-shot gathers which are free of the surface-related multiples.Under the assumption that the data is sparse,this work uses the method of convex optimization problem with L1-norm constraint to replace the traditional one to estimate primaries,which avoids using a time-window to prevent the inversion from into local optimization situations during sparse inversion,and improves the precision of the primaries estimated.It also suppresses artificial influence and improves the imaging quality.Comparing with the result obtained by lsqr algorithm shows the accuracy and superiority of the convex optimization problem with L1-norm constraint.During sparse inversion by L1-norm inversion,2DCurvelet transform and wavelet transform are combined to make the data more sparse,and improve the precision of primaries estimated.At the same time,the artificial influence suppression is improved further.Comparing with the traditional method to estimate primaries,the convex optimization problem with L1-norm constraint can avoid using a time-window to prevent the inversion from into local optimization situations,and the primaries estimated become more accurate.2DCurvelet transform and wavelet transform introduced make the data sparse,and improve the precision of primaries estimated.
引文
Baardman R H,Verschuur D J,Van Borselen R G,et al.2010.Estimation of primaries by sparse inversion using dual-sensordata.80th Annual International Meeting,SEG,Expanded Abstracts,3468-3472.
    Claerbout J F.1968.Synthesis of a layered medium from its acoustictransmission response.Geophysics,33(2):264-269.
    Feng F,Wang D L,Zhu H,et al.2013.Estimating primaries bysparse inversion of the 3Dcurvelet transform and the L1-normconstraint.Applied Geophysics(in Chinese),10(2):201-209.
    Lin T T Y,Herrmann F J.2011.Estimating primaries by sparse inversion in a curvelet-like representation domain.73rd EAGEConference and Exhibition,EAGE,Extended Abstracts,H043.
    Lin T T Y,Herrmann F J.2013.Robust estimation of primaries bysparse inversion via one-norm minimization.Geophysics,78(3):R133-R150.
    Van Den Berg E,Friedlander M P.2008.Probing the paretofrontier for basis pursuit solution.SIAM Journal on ScientificComputing,31(2):890-912.
    Van Groenestijn G J A,Verschuur D J.2009a.Estimatingprimaries by sparse inversion and application to near-offset datareconstruction.Geophysics,74(3):A23-A28.
    Van Groenestijn G J A,Verschuur D J.2009b.Estimation ofprimaries and near-offset reconstruction by sparse inversion:marine data application.Geophysics,74(6):R119-R128.
    Van Groenestijn G J A,Verschuur D J.2010.Estimation ofprimaries by sparse inversion from passive seismic data.Geophysics,75(4):SA61-SA69.
    Van Groenestijn G J A,Ross W.2011.Primary estimation on OBCdata by sparse inversion.81th Annual International Meeting,SEG,Expanded Abstracts,3531-3535.
    Vasconcelos I,Snieder R.2008a.Interferometry by deconvolution,Part 1—Theory for acoustic waves and numerical examples.Geophysics,73(3):S115-S128.
    Vasconcelos I,Snieder R.2008b.Interferometry by deconvolution:Part 2—Theory for elastic waves and application to drill-bitseismic imaging.Geophysics,73(3):S129-S141.
    Wapenaar K,van der Neut J,Ruigrok E.2008.Passive seismicinterferometry by multidimensional deconvolution.Geophysics,73(6):A51-A56.
    Wapenaar K,van der Neut J,Ruigrok E,et al.2011.Seismicinterferometry by cross-correlation and by multi-dimensionaldeconvolution:a systematic comparison.Geophysical JournalInternational,185(3):1335-1364.
    Ying L X,Demanet L,Candes E J.2005.3-D discrete curvelettransform.Proceedings SPIE 5914,Wavelets XI.San Diego:SPIE,344-354.
    Zhu H,Wang D L,Shi Z A,et al.2012.Passive seismic imaging ofseismic interferometry.Progress in Geophysics(in Chinese),27(2):496-502,doi:10.6038/j.issn.1004-2903.2012.02.012.
    冯飞,王德利,朱恒等.2013.三维曲波变换L1范数约束稀疏反演一次波估计方法研究.应用地球物理,10(2):201-209.
    朱恒,王德利,时志安等.2012.地震干涉技术被动源地震成像.地球物理学进展,27(2):496-502,doi:10.6038/j.issn.1004-2903.2012.02.012.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心