Acoustic VTI reverse time migration based on an improved source wavefield storage strategy
Ying Shi 1 3 Xiuzheng Fang 2 Weihong Wang 1 Xuan Ke 11 Northeast Petroleum University, School of Earth Science, Daqing, Heilongjiang 163318, China.
2 The Unconventional Natural Gas Institute, China University of Petroleum, Beijing 102249, China.
3 Corresponding author. Email: shiying@nepu.edu.cn
Exploration Geophysics 49(6) 891-897 https://doi.org/10.1071/EG17018
Submitted: 24 January 2017 Accepted: 23 January 2018 Published: 16 April 2018
Abstract
Advances in computational capabilities as well as ongoing improvements in storage strategies have made reverse time migration (RTM) a feasible method for capturing images of complex structures. However, large storage requirements still restrict RTM applications, especially in anisotropic media. Utilising a first-order quasi-P-wave equation in vertically transversely isotropic (VTI) media, we investigate anisotropy and deduce an RTM equation for a staggered-grid high-order finite difference (FD) scheme incorporating a perfectly matched layer (PML) boundary in this study. We also develop an improved source wavefield storage strategy via a PML boundary method for VTI medium RTM using graphic processing unit (GPU) accelerated computation. Our proposed method significantly reduces the total volume of data storage required for conventional RTM while increasing calculation time by just a small amount. Checkpoints can be set based on GPU memory size, leading to the generation of high precision and high efficiency subsurface images. We carried out a series of numerical tests on simple anisotropic media and complex Hess 2D VTI models to verify the effectiveness of our proposed method.
Key words: GPU, PML, reverse time migration, storage, VTI.
References
Alkhalifah, T., 1998, Acoustic approximations for processing in transversely isotropic media: Geophysics, 63, 623–631| Acoustic approximations for processing in transversely isotropic media:Crossref | GoogleScholarGoogle Scholar |
Alkhalifah, T., 2000, An acoustic wave equation for anisotropic media: Geophysics, 65, 1239–1250
| An acoustic wave equation for anisotropic media:Crossref | GoogleScholarGoogle Scholar |
Baysal, E., Kosloff, D., and Sherwood, J., 1983, Reverse time migration: Geophysics, 48, 1514–1524
| Reverse time migration:Crossref | GoogleScholarGoogle Scholar |
Clapp, R. G., 2009, Reverse time migration with random boundaries: SEG Technical Program Expanded Abstracts, 2809–2813.
Collino, F., and Tsogka, C., 2001, Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media: Geophysics, 66, 294–307
Du, X., Fletcher, R. P., and Fowler, P. J., 2008, A new pseudo-acoustic wave equation for VTI media: 70th EAGE Conference and Exhibition, Extended Abstracts, H033.
Dussaud, E., Symes, W. W., Williamson, P., Lemaistre, L., Singer, P., Denel, B., and Cherrett, A., 2008, Computational strategies for reverse-time migration: SEG Technical Program Expanded Abstracts, 2267–2271.
Duveneck, E., Milcik, P., Bakker, P. M., and Perkins, C., 2008, Acoustic VTI wave equations and their application for anisotropic reverse-time migration: SEG Technical Program Expanded Abstracts, 2186–2190.
Feng, B., and Wang, H., 2012, Reverse time migration with source wavefield reconstruction strategy: Journal of Geophysics and Engineering, 9, 69–74
| Reverse time migration with source wavefield reconstruction strategy:Crossref | GoogleScholarGoogle Scholar |
Han, F., and Sun, S. Z., 2015, GPU acceleration of amplitude-preserved Q compensation prestack time migration: Computers & Geosciences, 82, 214–224
| GPU acceleration of amplitude-preserved Q compensation prestack time migration:Crossref | GoogleScholarGoogle Scholar |
Hestholm, S., 2007, Acoustic VTI modeling using high-order finite-differences: SEG Technical Program Expanded Abstracts, 139–143.
Hestholm, S., 2009, Acoustic VTI modeling using high-order finite-differences: Geophysics, 74, T67–T73
| Acoustic VTI modeling using high-order finite-differences:Crossref | GoogleScholarGoogle Scholar |
Liu, Y., and Sen, M. K., 2010, Acoustic VTI modeling with a time-space domain dispersion-relation-based finite-difference scheme: Geophysics, 75, A11–A17
| Acoustic VTI modeling with a time-space domain dispersion-relation-based finite-difference scheme:Crossref | GoogleScholarGoogle Scholar |
McMechan, G., 1983, Migration by extrapolation of time-dependent boundary values: Geophysical Prospecting, 31, 413–420
| Migration by extrapolation of time-dependent boundary values:Crossref | GoogleScholarGoogle Scholar |
Pestana, R. C., Ursin, B., and Stoffa, P. L., 2012, Rapid expansion and pseudo spectral implementation for reverse time migration in VTI media: Journal of Geophysics and Engineering, 9, 291–301
| Rapid expansion and pseudo spectral implementation for reverse time migration in VTI media:Crossref | GoogleScholarGoogle Scholar |
Schiemenz, A., and Igel, H., 2013, Accelerated 3-D full-waveform inversion using simultaneously encoded sources in the time domain: application of Valhall ocean-bottom cable data: Geophysical Journal International, 195, 1970–1988
| Accelerated 3-D full-waveform inversion using simultaneously encoded sources in the time domain: application of Valhall ocean-bottom cable data:Crossref | GoogleScholarGoogle Scholar |
Shin, J., Ha, W., Jun, H., Min, D. J., and Shin, C., 2014, 3D Laplace-domain full waveform inversion using a single GPU card: Computers & Geosciences, 67, 1–13
| 3D Laplace-domain full waveform inversion using a single GPU card:Crossref | GoogleScholarGoogle Scholar |
Sun, W., and Fu, L., 2013, Two effective approaches to reduce data storage in reverse time migration: Computers & Geosciences, 56, 69–75
| Two effective approaches to reduce data storage in reverse time migration:Crossref | GoogleScholarGoogle Scholar |
Symes, W. W., 2007, Reverse time migration with optimal checkpointing: Geophysics, 72, SM213–SM221
| Reverse time migration with optimal checkpointing:Crossref | GoogleScholarGoogle Scholar |
Thomsen, L., 1986, Weak elastic anisotropy: Geophysics, 51, 1954–1966
| Weak elastic anisotropy:Crossref | GoogleScholarGoogle Scholar |
Virieux, J., 1984, SH-wave propagation in heterogeneous media: Velocity-stress finite-difference method: Geophysics, 49, 1933–1942
| SH-wave propagation in heterogeneous media: Velocity-stress finite-difference method:Crossref | GoogleScholarGoogle Scholar |
Virieux, J., 1986, P-SV wave propagation in heterogeneous media: Velocity-stress finite-difference method: Geophysics, 51, 889–901
| P-SV wave propagation in heterogeneous media: Velocity-stress finite-difference method:Crossref | GoogleScholarGoogle Scholar |
Wang, S. D., 2003, Absorbing boundary condition for acoustic wave equation by perfectly matched layer: Oil Geophysical Prospecting, 38, 31–34
Whitmore, D., 1983, Iterative depth migration by backward time propagation: SEG Technical Program Expanded Abstracts, 382–385.
Xie, W., Yang, D., Liu, F., and Li, J., 2014, Reverse-time migration in acoustic VTI media using a high-order stereo operator: Geophysics, 79, WA3–WA11
| Reverse-time migration in acoustic VTI media using a high-order stereo operator:Crossref | GoogleScholarGoogle Scholar |
Yan, H., and Liu, Y., 2013, Acoustic VTI modeling and pre-stack reverse-time migration based on time-space domain staggered-grid finite-difference method: Journal of Applied Geophysics, 90, 41–52
| Acoustic VTI modeling and pre-stack reverse-time migration based on time-space domain staggered-grid finite-difference method:Crossref | GoogleScholarGoogle Scholar |
Zhou, H., Zhang, G., and Bloor, R., 2006, An anisotropic acoustic wave equation for VTI media: 68th EAGE Conference and Exhibition incorporating EUROPEC, Extended Abstracts, 1171–1174.