局部倾角约束最小二乘偏移方法研究
|
摘要
随着石油勘探难度的进一步加大,地震数据往往存在采样不规则、地震道缺失等现象,如果不对其进行处理,会对后续的地震成像产生影响,引入成像噪音.针对这一问题,一般是通过地震道插值或数据规则化对叠前数据进行处理,然后采用常规的偏移方法进行成像,本文则是将地震成像看作最小二乘反演问题,在共成像点道集引入平滑算子,在共偏移距/角度道集引入平面波构造算子(PWC)进行约束,通过预条件共轭梯度法使得反偏移后数据与输入数据之间的误差达到最小,最终得到信噪比更高、振幅属性更为可靠的成像结果.理论模型和实际资料处理表明,本文方法不仅可以有效压制数据不规则对成像产生的噪音,而且具有更高的成像精度.
As the difficulty of oil exploration increases,the phenomenon of irregular sampling and missing traces often exits in seismic data,which will introduce imaging noise without special data processing.In order to solve the problem,the conventional method is implementing seismic trace interpolation or data regularization to pre-stack data before imaging with conventional migration.In this paper,seismic imaging is considered as least-squares inverse problem,with constraints of smoothing operator in common image gathers and plane-wave constructor in common offset/angle gathers,to obtain an artifact-reduced seismic image by iteratively minimizing the difference between de-migrated data and input data with preconditioning conjugate gradient method.Experimental results on theoretical model and seismic field data show that the proposed method can suppress the imaging noise introduced by data irregularity thus providing a more accurate image.
As the difficulty of oil exploration increases,the phenomenon of irregular sampling and missing traces often exits in seismic data,which will introduce imaging noise without special data processing.In order to solve the problem,the conventional method is implementing seismic trace interpolation or data regularization to pre-stack data before imaging with conventional migration.In this paper,seismic imaging is considered as least-squares inverse problem,with constraints of smoothing operator in common image gathers and plane-wave constructor in common offset/angle gathers,to obtain an artifact-reduced seismic image by iteratively minimizing the difference between de-migrated data and input data with preconditioning conjugate gradient method.Experimental results on theoretical model and seismic field data show that the proposed method can suppress the imaging noise introduced by data irregularity thus providing a more accurate image.
引文
[1]Claerbout J F.Towards a unified theory of reflectormapping.Geophysics,1971,36(3):467-481.
[2]Lailly P.The seismic inverse problem as a sequence of beforestack migration.∥Bednar J ed.Conference on InverseScattering:Theory and Applications,SIAM,ExpandedAbstracts,1983:206-220.
[3]Tarantola A.Inversion of seismic reflection data in theacoustic approximation.Geophysics,1984,49(8):1259-1266.
[4]Kuehl H,Sacchi W D.Least-squares wave-equation migrationfor AVP/AVA inversion.Geophysics,2003,68(1):262-273.
[5]Clapp M L.Imaging under salt:Illumination compensationby regularized inversion[Ph.D.thesis].Stanford:StanfordUniv.,2005.
[6]Valenciano A A,Biondi B,Guitton A.Target-oriented wave-equation inversion.Geophysics,2006,71(4):A35-A38.
[7]Prucha M,Biondi B.Subsalt event regularization withsteering filters.∥72nd Annual International Meeting,SEG,Expanded Abstracts,2002:824-827.
[8]Tang Y X.Target-oriented wave-equation least-squaresmigration/inversion with phase-encoded Hessian.Geophysics,2009,74(6):WCA95-WCA107.
[9]Wang J,Sacchi M D.Structure constrained least-squaresmigration.∥79th Ann.Internat Mtg.,Soc.Expl.Geophys.Expanded Abstracts,2009:2763-2767.
[10]Symes W W.Approximate linearized inversion by optimalscaling of prestack depth migration.Geophysics,2008,73(2):R23-R35.
[11]Rickett J E.Illumination-based normalization for wave-equation depth migration.Geophysics,2003,68(4):1371-1379.
[12]Aoki N,Schuster G T.Fast least-squares migration with adeblurring filter.Geophysics,2009,74(6):WCA83-WCA93.
[13]Herrmann F J,Brown C R,Erlangga Y A,et al.Curvelet-based migration preconditioning and scaling.Geophysics,2009,74(4):A41-A46
[14]Stolt R H,Benson A.Seismic Migration:Theory andPractice.Amsterdam:Geophysical Press,1986.
[15]Claerbout J F.Earth Soundings Analysis-Processing VersusInversion.Blackwell Scientific Publications,1992.
[16]Hestenes M,Steifel E.Methods of conjugate gradients forsolving linear systems.Journal of Research of the NationalInstitute of Standards,1952,49(6):409-436.
[17]Fomel S B,Guitton A.Regularizing seismic inverse problemsby model reparameterization using plane-wave construction.Geophysics,2006,71(5):A43-A47.
[18]Fomel S B.Applications of plane-wave destruction filters.Geophysics,2002,67(6):1946-1960.
[19]Fomel S B.Shaping regularization in geophysical estimationproblems.∥75th Annual International Meeting,SEG ExpandedAbstracts,2005:1673-1676.
[20]Fomel S B.Three Dimension Seismic Data Regularization.[Ph.D.thesis].Stanford:Stanford Univ.,2001.
[2]Lailly P.The seismic inverse problem as a sequence of beforestack migration.∥Bednar J ed.Conference on InverseScattering:Theory and Applications,SIAM,ExpandedAbstracts,1983:206-220.
[3]Tarantola A.Inversion of seismic reflection data in theacoustic approximation.Geophysics,1984,49(8):1259-1266.
[4]Kuehl H,Sacchi W D.Least-squares wave-equation migrationfor AVP/AVA inversion.Geophysics,2003,68(1):262-273.
[5]Clapp M L.Imaging under salt:Illumination compensationby regularized inversion[Ph.D.thesis].Stanford:StanfordUniv.,2005.
[6]Valenciano A A,Biondi B,Guitton A.Target-oriented wave-equation inversion.Geophysics,2006,71(4):A35-A38.
[7]Prucha M,Biondi B.Subsalt event regularization withsteering filters.∥72nd Annual International Meeting,SEG,Expanded Abstracts,2002:824-827.
[8]Tang Y X.Target-oriented wave-equation least-squaresmigration/inversion with phase-encoded Hessian.Geophysics,2009,74(6):WCA95-WCA107.
[9]Wang J,Sacchi M D.Structure constrained least-squaresmigration.∥79th Ann.Internat Mtg.,Soc.Expl.Geophys.Expanded Abstracts,2009:2763-2767.
[10]Symes W W.Approximate linearized inversion by optimalscaling of prestack depth migration.Geophysics,2008,73(2):R23-R35.
[11]Rickett J E.Illumination-based normalization for wave-equation depth migration.Geophysics,2003,68(4):1371-1379.
[12]Aoki N,Schuster G T.Fast least-squares migration with adeblurring filter.Geophysics,2009,74(6):WCA83-WCA93.
[13]Herrmann F J,Brown C R,Erlangga Y A,et al.Curvelet-based migration preconditioning and scaling.Geophysics,2009,74(4):A41-A46
[14]Stolt R H,Benson A.Seismic Migration:Theory andPractice.Amsterdam:Geophysical Press,1986.
[15]Claerbout J F.Earth Soundings Analysis-Processing VersusInversion.Blackwell Scientific Publications,1992.
[16]Hestenes M,Steifel E.Methods of conjugate gradients forsolving linear systems.Journal of Research of the NationalInstitute of Standards,1952,49(6):409-436.
[17]Fomel S B,Guitton A.Regularizing seismic inverse problemsby model reparameterization using plane-wave construction.Geophysics,2006,71(5):A43-A47.
[18]Fomel S B.Applications of plane-wave destruction filters.Geophysics,2002,67(6):1946-1960.
[19]Fomel S B.Shaping regularization in geophysical estimationproblems.∥75th Annual International Meeting,SEG ExpandedAbstracts,2005:1673-1676.
[20]Fomel S B.Three Dimension Seismic Data Regularization.[Ph.D.thesis].Stanford:Stanford Univ.,2001.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心