Abstract：
Prestack reverse time migration (PSTM) is a common imaging method; however low-frequency noises reduce the structural imaging precision. Thus, the suppression of migration noises must be considered. The generation mechanism of low-frequency noises is analyzed and the up-, down-, left-, and right-going waves are separated using the Poynting vector of the acoustic wave equation. The computational complexity and memory capacitance of the proposed method are far smaller than that required when using the conventional separation algorithm of 2D Fourier transform. The normalized wavefield separation cross-correlation imaging condition is used to suppress low-frequency noises in reverse time migration and improve the imaging precision. Numerical experiments using the Marmousi model are performed and the results show that the up-, down-, left-, and right-going waves are well separated in the continuation of the wavefield using the Poynting vector. We compared the imaging results with the conventional method, Laplacian filtering, and wavefield separation with the 2D Fourier transform. The comparison shows that the migration noises are well suppressed using the normalized wavefield separation cross-correlation imaging condition and higher precision imaging results are obtained.
CHEN Ting,HE Bing-Shou. A normalized wavefield separation cross-correlation imaging condition for reverse time migration based on Poynting vector[J]. APPLIED GEOPHYSICS, 2014, 11(2): 158-166.
[1]
Baysal, E., Kosloff, D. D., and Sherwood, J. W. C., 1983, Reverse time migration: Geophysics, 48(11), 1514-1524.
[2]
Berenger, J. P., 1994, A perfectly matched layer for the absorption of electromagnetic waves: Journal of Computational Physics, 114(2), 185-200.
[3]
Claerbout, J. F., 1971, Toward a unified theory of reflector mapping: Geophysics, 36(3), 467-481.
[4]
Claerbout, J. F., and Doherty, S. M., 1972, Downward continuation of moveout-corrected seismograms: Geophysics, 37(5), 741-768.
[5]
Clapp, R. G., 2008, Reverse time migration: Saving the boundaries: 136, 136-144.
[6]
Clapp, R. G., 2009, Reverse time migration with random boundaries: 79th Annual International Meeting., SEG Expanded Abstracts, 2809-2813.
[7]
Dong, L. G., Ma, Z. T., Cao, J. Z., et al., 2000, A staggered-grid high-order difference method of one-order elastic wave equation: Chinese Journal of Geophysics(in Chinese), 43(3), 411-419.
[8]
Deng, S. Q., 2012, Study on numerical simulation of whole-space elastic wave and reverse time migration imaging method: PhD Thesis, China University of Mining and Technology, Xuzhou.
[9]
Jin, S., Jiang, F., Ren, Y., 2010, Comparison of Isotropic VTI And TTI Reverse Time Migration: an Experiment On BP Anisotropic Benchmark Dataset: 80th Annual International Meeting., SEG Expanded Abstracts, 3198-3202.
[10]
Li, B., Liu, H. W., Liu, G. F., et al., 2010, Computational strategy of seismic pre-stack reverse time migration on CPU/GPU: Chinese Journal of Geophysics (in Chinese), 53(12), 2938-2943.
[11]
Liu, F. Q., Zhang, G. Q., and Morton, S. A., 2011, An effective imaging condition for reverse-time migration using wavefield decomposition: Geophysics, 76(1), S29-S39.
[12]
Liu, Y., and Sen M.K., 2012, A hybrid absorbing boundary condition for elastic staggered-grid modeling: Geophysical Prospecting, 60(6), 1114 - 1132.
[13]
Loewenthal, D., Stoffa, P. L., and- Faria, E. L., 1987, Suppressing the unwanted reflections of the full wave equation: Geophysics, 52(7), 1007-1012.
[14]
Loewenthal, D., and Mufti, I., 1983. Reversed time migration in spatial frequency domain: Geophysics, 48(5), 627-635.
[15]
McMechan, G, A., 1983, Migration by extrapolation of time-dependent boundary values: Geophysical Prospecting, 31(2), 413-420.
[16]
Mu, Y. G., and Pei, Z. L., 2005, Seismic numerical modeling for 3-D complex media: Petroleum Industry Press, China, 33-34.
[17]
Poynting, J. H., 1884, On the transfer of energy in the electromagnetic field: Philosophical Transactions of the Royal Society of London, 175, 343-361.
[18]
Schleicher, J., Costa, J., and Novais, A., 2008, A comparison of imaging conditions for wave-equation shot-profile migration: Geophysics, 73(6), S219-S227.
[19]
Sun, R., McMechan, G. A., Lee, C. S., et al., 2006, Prestack scalar reverse-time depth migration of 3D elastic seismic data: Geophysics, 71(5), S199-S207.
[20]
Symes, W. W., 2007, Reverse time migration with optimal checkpointing: Geophysics, 72(5), SM213-SM221.
[21]
Tang, C., Wang, D. L., 2012, Reverse time migration with source wavefield reconstruction and wavefield decomposition: Global Geology(in Chinese), 31(4), 803-812.
[22]
Whitmore, N. D., 1983, Iterative depth migration by backward time propagation: 53th Annual International Meeting., SEG Expanded Abstracts, 382-385.
[23]
Yan, H.Y., and Liu, Y., 2013, Acoustic prestack reverse time migration using the adaptive high-order finite-difference method in time-space domain: Chinese Journal of Geophysics (in Chinese), 56(3), 971 - 984.
[24]
Yoon, K., Marfurt, K. J., 2006, Reverse-time migration using the Poynting vector: Exploration Geophysics, 37(1), 102-107.
[25]
Zhang, J.H, and Yao, Z.X., 2013, Optimized finite-difference operator for broadband seismic wave modeling: Geophysics, 78(1), A13 - A18.
[26]
Zhang, Y., Sun, J., 2008, Practical issues of reverse time migration: true amplitude gathers, noise removal and harmonic-source encoding: 70th EAGE Conference & Exhibition, Rome, Italy.
[27]
Zhang, Y., Zhang, H., Zhang, G., 2011, A stable TTI reverse time migration and its implementation: Geophysics, 76(3), WA3-WA11.
2014, Vol. 11 
