基于解析时间波场外推与波场分解的逆时偏移方法研究  

Reverse time migration using analytical time wavefield extrapolation and separation

作　　者：胡江涛[1] 王华忠[1] 

机构地区：[1]同济大学海洋与地球科学学院波现象与反演成像研究组,上海200092

出　　处：《地球物理学报》2015年第8期2886-2895,共10页Chinese Journal of Geophysics

基　　金：973项目(2011CB201002);国家自然科学基金(41374117);国家重大专项(2011ZX05005-005-008HZ;2011ZX05006-002;2011ZX05023)资助

摘　　要：双程波方程逆时深度偏移是复杂介质高精度成像的有效技术,但其结果中通常包含成像方法引起的噪音和假象,一般的滤波方法会破坏成像剖面上的振幅,其中的假象也会给后续地质解释带来困扰.将波场进行方向分解然后实现入射波与反射波的相关成像能够有效地消除这类成像噪音,并提高逆时偏移成像质量.波传播方向的分解通常在频率波数域实现,它会占用大量的存储和计算资源,不便于在沿时间外推的逆时深度偏移中应用.本文提出解析时间波场外推方法,可以在时间外推的每个时间片上实现波传播方向的显式分解,逆时深度偏移中利用分解后的炮检波场进行对应的相关运算,实现成像噪音和成像信号的分离.在模型和实际数据上的测试表明,相比于常规互相关逆时偏移成像结果,本文方法能够有效地消除低频成像噪音和特殊地质构造导致的成像假象.Two-way wave equation based reverse time migration is powerful in imaging complex area.However,its result usually contains imaging noise which poses challenge to seismic interpretation.The imaging amplitude would be distorted by simple noise filter. Wavefield separation imaging condition is an effective way to suppress the imaging noise and enhance the imaging quality.Conventionally,the wavefield separation is achieved in frequency wavenumber domain which cannot be conveniently implemented in time domain reverse time migration.We propose the analytical wavefield extrapolation method which can efficiently separate wavefield into its directional components during the time domain extrapolation.We use this method in reverse time migration and separate both the source and receiver wavefield.And the imaging noise and signal can be separated after applying the imaging condition to the separated wavefield.The wavefield propagation direction is usually defined by temporal and spatial Fourier transform.In the Fourier domain,the wavefield propagation direction is defined by the sign of frequency and spatial wavenumber.In order to efficiently apply the wavefield separation imagingcondition in reverse time migration,we extend the analytical time signal into the wavefield and call it analytical time wavefield.Since the analytical time wavefield contains only the positive frequency component,the wavefield propagation direction can be defined by the sign of spatial wavenumber.To avoid the I/O cost in generating the analytical wavefield at every imaging point,we propose an analytical wavefield propagation equation based on the linear relation between the source term and wavefield in wave equation.We solve the proposed equation by finite difference method.Then we separate the source and receiver wavefield into their up and down going components and apply the imaging condition to the separated wavefield. Four imaging components（i.e.,imaging component from up going source and down going receiver wavefield,imaging component from do

关 键 词：解析时间波场外推 波场分解 噪音消除 逆时深度偏移 

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