3D surface-related multiple elimination based on sparse inversion
-
摘要: 三维表面多次波压制是海洋地震资料预处理中的重要研究课题,基于波动理论的三维表面多次波压制方法(3D SRME)是数据驱动的方法,理论上来说,可有效压制复杂构造地震数据表面多次波.但该方法因对原始地震数据采集要求高而很难在实际资料处理中广泛应用.本文基于贡献道集的概念,将稀疏反演方法引入到表面多次波压制中,应用稀疏反演代替横测线积分求和,无需对横测线进行大规模重建,进而完成三维表面多次波预测,这样可有效解决实际三维地震数据横测线方向稀疏的问题.基于纵测线多次波积分道集为抛物线的假设,为保证预测后三维表面多次波和全三维数据预测的多次波在运动学和动力学特征上基本一致,文中对预测数据实施基于稳相原理的相位校正.理论模型和实际数据的测试结果表明,本文基于稀疏反演三维表面多次波压制方法可在横测线稀疏的情况下,有效压制三维复杂介质地震资料中的表面多次波,从而更好地提高海洋地震资料的信噪比,为高分辨率地震成像提供可靠的预处理数据保障.Abstract: Three dimensional surface-related multiple elimination (SRME) is one of the important topics in the processing of seismic data from marine exploration,theoretically, the data-driven SRME based on wave-equation, can suppress all surface-related multiples from complex structure, both in 2D and 3D sense. But actually, because of the high requirement for seismic data acquisition, it is usually difficult to apply 3D SRME for field data demultiple processing.3D multiple suppression approach based on sparse inversion is analyzed.We classify multiple suppression method into two categories, filter and SRME method respectively. For the seismic data from complex geological structure, filter approach doesn't work well. However, the data-driven SRME approach based on wave equation, can suppress multiple better, which has no requirement for velocity information.For SRME, the full wavefield data requirement is an important disadvantage, which cannot be meet for almost all marine field data. Therefore the data reconstruction is necessary for traditional multiple suppression using SRME.For current marine acquisition geometries, the data is densely sampled in the inline direction, but very sparsely in the crossline direction, we introduce contribution gather concept, and calculate sparsely sampled crossline multiple contribution by means of sparse inversion algorithm.Therefore, the data reconstruction is unnecessary before demultiple,comparing with traditional 3D SRME algorithm, which decrease storage cost greatly. Based on the assumption that the crossline time-distance curves are hyperbolic or parabolic, i.e., after integrated along inline direction after finished the first step processing of contribution gathers, we applied the phase correction algorithm which based on the principle of stationary phase approximation, predicted multiples by the proposed method and result by the full 3D SRME method are basically the same on kinematics and dynamics characteristics. Multiple prediction and adaptive subtraction are two crucial steps for 3D multiple suppression using SRME. We simulate complex model data to test the proposed 3D multiple prediction algorithm, and 10 thousands shot records are modelled, the shot interval and trace interval are both 25 m. The result comparison show the 3D multiple prediction approach can predict the multiple's amplitude and phase correctly, and also, the subtraction result is superior than 2D algorithm. The horizontal four-layered media is also designed to test the sparse inversion 3D multiple suppression algorithm, there are 3136 shot records in total, the trace and shot intervals are both 25 m, the line interval is 75 m. The single trace and common-offset result show that the proposed approach can predict the multiple's amplitude and traveltime correctly, and they are very close to the full data circumstance. The test on field data from some area in China show that the proposed sparse inversion method is applicable and effective, where the trace interval is 12.5 m, shot interval is 50 m, and the line interval is 100 m. After summation along inline direction, the partial integration data is transformed to Radon domain using apex-shifted Radon approach based on the assumption of hyperbolic or parabolic events, stacked the Radon imaging, and also applied the phase correction to the sparse inversion solution, the predicted multiple is acquired.After theoretical investigation and data tests, we have the following conclusions, (1) the proposed method is suitable for simple and complex 3D model data, (2) for real seismic data, the inline direction reconstruction is needed, (3) because only forward Radon transform, not inverse Radon transform is used,phase correction is demanded, and (4) the algorithm do not rely on the assumption of full data using sparse inversion approach, and also, widen the application extent of 3D multiple suppression.The result show that the proposed 3D multiple suppression algorithm can improve the S/N ratio during the course of preprocessing, and provide the high quality data for the subsequent high resolution imaging.
-
Key words:
- Sparse inversion /
- Surface-related multiple /
- Contribution gather /
- Multiple prediction
-
[1] Abma R, Kabir N, Matson K H, et al. 2005. Comparisons of adaptive subtraction methods for multiple attenuation. The Leading Edge, 24(3): 277-280.
[2] Berkhout A J. 1980. Seismic Migration: Imaging of Acoustic Energy by Wave Field Extrapolation. Part A: Theoretical Aspects. Elsevier.
[3] Berkhout A J, Verschuur D J. 1997. Estimation of multiple scattering by iterative inversion, part I: Theoretical considerations. Geophysics, 62(5): 1586-1595.
[4] Bisley R, MacNeill M. 2008. Application of true-azimuth 3D SRME in the Northwest Shelf, Australia. //78th Annual International Meeting, SEG, Expanded Abstracts, 2436-2440.
[5] Bleistein N. 1984. Mathematical Methods for Wave Phenomena. Academic Press.
[6] Dragoset B, Moore I, Kostov C. 2006. The impact of field-survey characteristics on surface-related multiple attenuation. Geophysical Prospecting, 54(6): 781-791.
[7] Dragoset B, Verschuur D J, Moore I, et al. 2010. A perspective on 3D surface-related multiple elimination. Geophysics, 75(5): 75A245-75A261.
[8] Dragoset W H, Jeričević Ž. 1998. Some remarks on surface multiple attenuation. Geophysics, 63(2): 772-789.
[9] Hokstad K, Sollie R. 2003. 3D surface-related multiple elimination using parabolic sparse inversion. //73rd Annual International Meeting, SEG, Expanded Abstracts, 1961-1964.
[10] Jakubowicz H. 1998. Wave equation prediction and removal of interbed multiples. //68th Annual International Meeting, SEG, Expanded Abstracts, 1527-1530.
[11] Kabir M N, Verschuur D J. 1995. Restoration of missing offsets by parabolic Radon transform. Geophysical Prospecting, 43(3): 347-368.
[12] Kelamis P G, Verschuur D J. 2000. Surface-related multiple elimination on land seismic data—Strategies via case studies. Geophysics, 65(3): 719-734.
[13] Li X, Hu T Y. 2009. Surface-related multiple removal with inverse scattering series method. Chinese J. Geophy. (in Chinese), 52(6): 1633-1640, doi: 10.3969/j.issn.0001-5733.2009.06.026.
[14] Matson K H, Abma R. 2005. Fast 3D surface-related multiple elimination using azimuth moveout for multiples. //75th Annual International Meeting, SEG, Expanded Abstracts, 2064-2067.
[15] Ross W S, Yu Y, Gasparotto F A. 1999. Travel time prediction and suppression of 3-D multiples. Geophysics, 64(1): 261-277.
[16] Sacchi M D, Ulrych T J. 1995. High-resolution velocity gathers and offset space reconstruction. Geopyhsics, 60(4): 1169-1177.
[17] Shi Y, Wang W H, Li Y, et al. 2013. 3D surface-related multiple prediction approach investigation based on wave equation. Chinese J. Geophys. (in Chinese), 56(6): 2023-2032, doi: 10.6038/cjg20130623.
[18] Sun Y. 1999. Anti-aliasing multiple prediction beyond 2-D. //69th Annual International Meeting SEG, Expanded Abstracts, 1338-1341
[19] Trad D, Ulrych T, Sacchi M. 2003. Latest views of the sparse Radon transform. Geophysics, 68(1): 386-399.
[20] van Dedem E J, Verschuur D J. 2001. 3D surface multiple prediction using sparse inversion. //71st Annual International Meeting, SEG, Expanded Abstracts, 1285-1288.
[21] van Dedem E J, Verschuur D J. 2005. 3D surface-related multiple prediction: A sparse inversion approach. Geophysics, 70(3): V31-V43.
[22] Verschuur D J. 1991. Surface-related multiple elimination, an inversion approach. Netherlands: Delft University of Technology.
[23] Verschuur D J, Berkhout A J. 1997. Estimation of multiple scattering by iterative inversion, Part II: Practical aspects and examples. Geophysics, 62(5): 1596-1611.
[24] Verschuur D J, Berkhout A J, Wapenaar C P A. 1992. Adaptive surface-related multiple elimination. Geophysics, 57(9): 1166-1177.
[25] Wang B L, Sacchi M D, Yin X Y, et al. 2014. Multiple attenuation based on amplitude preserving Radon transform. Chinese J. Geophys. (in Chinese), 57(6): 1924-1933, doi: 10.6038/cjg20140623.
[26] Wang Y H. 2003. Sparseness-constrained least-squares inversion: Application to seismic wave reconstruction. Geophysics, 68(5): 1633-1638.
[27] Wapenaar C P A. 1992. The stationary phase method and its application in seismics and acoustics. Lecture notes: Akoestische Beeldvorming, Delft University of Technology.
计量
- 文章访问数: 1904
- PDF下载数: 2003
- 施引文献: 0