Methodology of waveform inversion for acoustic orthorhombic media

Wang, H. and Tsvankin, I., 2018. Methodology of waveform inversion for acoustic orthorhombic media. Journal of Seismic Exploration, 27: 201-226. Three-dimensional seismic waveform inversion (WI) for anisotropic media is highly challenging due to its computational cost and trade-offs between multiple model parameters. Here, we develop a methodology of 3D WI for orthorhombic media in the acoustic approximation using two mixed-domain modeling algorithms, one of which is based on low-rank decomposition and the other on the generalized pseudospectral method. Numerical testing shows that both techniques produce kinematically accurate compressional wavefields with an acceptable computational cost. To take advantage of the superior stability and accuracy of the low-rank-decomposition-based method, it is employed to simulate both the forward and adjoint wavefields. The gradient of the data- difference objective function, however, is more convenient to obtain from the wave equations derived with the pseudospectral method. The inversion is conducted with a limited-memory version of the quasi-Newton optimization algorithm. Under the assumption that the symmetry-plane orientation is known, we invert wide-azimuth data for all six parameters of acoustic orthorhombic media. The performance of the developed wavefield-extrapolation and gradient-computation algorithms is evaluated for a medium with Gaussian anomalies in each parameter. Then we apply the method to a modified SEG/EAGE overthrust model to demonstrate the feasibility of waveform inversion and illustrate parameter trade-offs for structurally complicated orthorhombic media.
- Albertin, U., Shen, P., Sekar, A., Johnsen, T., Wu, C., Nihei, K. and Bube, K., 2016. 3D
- orthorhombic elastic full-waveform inversion in the reflection domain from
- hydrophone data. Expanded Abstr., 86th Ann. Internat. SEG Mtg., Dallas: 1094-
- Alkhalifah, T., 1998. Acoustic approximations for processing in transversely isotropic
- media. Geophysics, 63: 623-631.
- Alkhalifah, T., 2000. An acoustic wave equation for anisotropic media. Geophysics, 65:
- 1239-1250.
- Alkhalifah, T., Masmoudi, N. and Oh, J.-W., 2016. A recipe for practical full-waveform
- inversion in orthorhombic anisotropy. Expanded Abstr., 86th Ann. Internat. SEG
- Mtg., Dallas: 317-321.
- Aminzadeh, F., Brac, J. and Kunz, T., 1997. 3D salt and overthrust model. Modeling Series
- I, SEG, Tulsa, OK.
- Baydin, V., Charara, M. and Dovgilovich, L., 2015. Elastic orthorhombic full-waveform
- inversion for 3D VSP case. Extended Abstr., 77th EAGE Conf., Madrid: 1-5.
- Benson, S.J. and Moré, J.J., 2001. A limited memory variable metric method in
- subspaces and bound constrained optimization problems. Technical Report,
- ANL/ACS-P909-0901, Argonne National Labor.
- Birdus, S., Sun, J, Sun, W., Xie, Y., Gazzoli, M., Andreolli, M. and Ursulic, A., 2012.
- Multi-azimuth PSDM processing in the presence of orthorhombic anisotropy — a
- case history offshore north west Australia. Expanded Abstr., 82nd Ann. Internat.
- SEG Mtg., Las Vegas: 1-5.
- Chen, P., 2011. Full-wave seismic data assimilation: theoretical background and recent
- advances. Pure Appl. Geophys., 168: 1527-1552.
- Du, X., Fowler, P.J. and Fletcher, R.P., 2014. Recursive integral time-extrapolation
- methods for waves: A comparative review. Geophysics, 79(1): T9-T26.
- Duveneck, E. and Bakker, P.M., 2011. Stable P-wave modeling for reverse-time migration
- in tilted TI media. Geophysics, 76(2): S65—S75.
- Etgen, J.T. and Dellinger, J., 1989. Accurate wave equation modeling. Expanded Abstr.,
- 59th Ann. Internat. SEG Mtg., Dallas: 494-497.
- Fichtner, A., 2010. Full Seismic Waveform Modelling and Inversion. Springer Verlag,
- Berlin.
- Fichtner, A., Bunge, H.-P. and Igel, H.,2006a. The adjoint method in seismology: I. Theory.
- hys. Earth Planet. Inter., 157: 86-104.
- Fichtner, A., Bunge, H.-P. and Igel, H., 2006b. The adjoint method in seismology: II.
- Applications: traveltimes and sensitivity functionals. Phys. Earth Planet. Inter.,
- 157: 105-123.
- Fletcher, R.P., Du, X. and Fowler, P., 2008. A new pseudo-acoustic wave equation for TI
- media. Expanded Abstr., 78th Ann. Internat. SEG Mtg., Las Vegas: 2082-2086.
- Fletcher, R.P., Du, X. and Fowler, P.J., 2009. Reverse time migration in tilted transversely
- isotropic (TTD media. Geophysics, 74(6): WCA179-WCA187.
- Fomel, S., Sava, P., Vlad, I, Liu, Y. and Bashkardin, V., 2013a. Madagascar: Open-source
- software project for multidimensional data analysis and _ reproducible
- computational experiments. J. Open Res. Softw., 1: e8.
- Fomel, S., Ying, L. and Song, X., 2013b. Seismic wave extrapolation using lowrank symbol
- approximation: Geophys. Prosp., 61: 526-536.
- Fowler, P., 2015. Some pitfalls in orthorhombic parameter estimation. Expanded Abstr.,
- 85th Ann. Internat. SEG Mtg., New Orleans: 493-497.
- Fowler, P.J., Du, X. and Fletcher, R.P., 2010. Coupled equations for reverse time
- migration in transversely isotropic media. Geophysics, 75(1): S$11-S22.
- Fowler, P.J. and King, R., 2011. Modeling and reverse time migration of orthorhombic
- pseudoacoustic P-waves Expanded Abstr., 81st Ann. Internat. SEG Mtg., San
- Antonio: 190-195.
- Fowler, P.J. and Lapilli, C., 2012. Generalized pseudospectral methods for orthorhombic
- modeling and reverse-time migration. Expanded Abstr., 82nd Ann. Internat. SEG
- Mtg., Las Vegas: 1-5.
- Gholami, Y., Brossier, R., Operto, S., Prieux, V., Ribodetti, A. and Virieux, J., 2013.
- Which parameterization is suitable for acoustic vertical transverse isotropic full
- waveform inversion? Part 2: Synthetic and real data case studies from Valhall.
- Geophysics, 78(2): R107-R124.
- Grechka, V., Zhang, L. and Rector, J.W., 2004. Shear waves in acoustic anisotropic
- media. Geophysics, 69: 576-582.
- Kamath, N. and Tsvankin, I., 2013. Full-waveform inversion of multicomponent data for
- horizontally layered VTI media. Geophysics, 78(5): WC113-WC121.
- Kamath, N. and Tsvankin, I., 2016. Elastic full-waveform inversion for VTI media:
- Methodology and sensitivity analysis. Geophysics, 81(2): C53-C68.
- Karazincir, M. and Orumwense, R., 2014. Tilted orthorhombic velocity model building and
- imaging of Zamzama gas field with full-azimuth land data. The Leading Edge, 33:
- 1024-1028.
- Kolda, T.G., O’Leary, D.P. and Nazareth, L., 1998. BFGS with update skipping and varying
- memory. SIAM J. Optimizat., 8: 1060-1083.
- Kosloff, D.D. and Baysal, E., 1982. Forward modeling by a Fourier method. Geophysics,
- 47: 1402-1412.
- Lailly, P., 1983. The seismic inverse problem as a sequence of before stack migrations.
- Conf. Inver. Scatter.: Theory Applicat, Soc. Industr. Appl. Mathemat.,
- Philadelphia, PA: 206-220.
- Liu, D.C. and Nocedal, J., 1989. On the limited memory BFGS method for large scale
- optimization. Mathemat. Programm., 45: 503-528.
- Liu, Q. and Tromp, J. 2008. Finite-frequency sensitivity kernels for global seismic wave
- propagation based upon adjoint methods. Geophys. J. Internat., 174: 265-286.
- Luo, J. and Wu, R.-S., 2015. Seismic envelope inversion: reduction of local minima and
- noise resistance. Geophys. Prosp., 63: 597-614.
- Masmoudi, N. and Alkhalifah, T., 2016. A new parametrization for waveform inversion in
- acoustic orthorhombic media. Geophysics, 81(4): R157-R171.
- Mathewson, J., Bachrach, R., Ortin, M., Decker, M., Kainkaryam, S., Cegna, A., Paramo, P.,
- Vincent, K. and Kommedal, J., 2015. Offshore Trinidad tilted orthorhombic
- imaging of full-azimuth OBC data. Expanded Abstr., 85th Ann. Internat. SEG
- Mtg., New Orleans: 356-360.
- Nocedal, J., 1992. Theory of algorithms for unconstrained optimization. Acta Numer., 1:
- 199-242.
- Nocedal, J. and Wright, S., 2006. Numerical Optimization. Springer Science & Business
- Media.
- Oh, J.-W. and Alkhalifah, T., 2016. 3D elastic-orthorhombic anisotropic full-waveform
- inversion. Application to field OBC data. Expanded Abstr., 86th Ann. Internat.
- SEG Mtg., Dallas: 1206-1210.
- Pestana, R.C. and Stoffa, P.L., 2010. Time evolution of the wave equation using rapid
- expansion method: Geophysics, 75(4): T121-T131.
- Plessix, R.-E., 2006. A review of the adjoint-state method for computing the gradient of a
- functional with geophysical applications. Geophys. J. Internat., 167: 495-503.
- Royle, G., 2011. Viscoelastic orthorhombic full wavefield inversion: Development of
- multiparameter inversion methods. Expanded Abstr., 81st Ann. Internat. SEG
- Mtg., San Antonio: 4329-4333.
- Song, X. and Alkhalifah, T., 2013. Modeling of pseudo-acoustic P-waves in
- orthorhombic media with low-rank approximation. Geophysics, 78(4): C33-C40.
- Tal-Ezer, H., 1986. Spectral methods in time for hyperbolic equations. SIAM J. Numer.
- Analys., 23: 11-26.
- Tal-Ezer, H., Kosloff, D. and Koren, Z., 1987. An accurate scheme for seismic forward
- modeling. Geophys. Prosp., 35: 479-490.
- Tarantola, A., 1984. Inversion of seismic reflection data in the acoustic approximation.
- Geophysics, 49: 1259-1266.
- Tarantola, A., 1988. Theoretical background for the inversion of seismic waveforms
- including elasticity and attenuation. Pure Appl. Geophys., 128: 365-399.
- Thiébaut, E., 2002. Optimization issues in blind deconvolution algorithms.
- Astronomical Telescopes and Instrumentation, Internat. Soc. Optics Photonics,
- Thomas, M., Mothi, S. and McGill, P., 2012. Improving subsalt images using tilted-
- orthorhombic RTM in Green Canyon, Gulf of Mexico. Expanded Abstr., 82nd
- Ann. Internat. SEG Mtg., Las Vegas: 1-5.
- Tromp, J., Tape, C. and Liu, Q., 2005. Seismic tomography, adjoint methods, time
- reversal and banana-doughnut kernels. Geophys. J. Internat., 160: 195-216.
- Tsvankin, I., 1997. Anisotropic parameters and P-wave velocity for orthorhombic media.
- Geophysics, 62: 1292-1309.
- Tsvankin, I., 2012. Seismic Signatures and Analysis of Reflection Data in Anisotropic
- Media. SEG, Tulsa, OK.
- Tsvankin, I. and Grechka, V., 2011. Seismology of Azimuthally Anisotropic Media and
- Seismic Fracture Characterization. SEG, Tulsa, OK.
- Wang, S. and Wilkinson, D., 2012. Imaging for unconventional resource plays using an
- orthorhombic velocity model. Expanded Abstr., 82nd Ann. Internat. SEG Mtg.,
- Las Vegas:, 1-5.
- Wu, R.-S., Luo, J. and Wu, B., 2014. Seismic envelope inversion and modulation signal
- model. Geophysics, 79(3): WA13-WA24.
- Xie, Y., Zhou, B., Zhou, J., Hu, J., Xu, L., Wu, X., Lin, N., Loh, F.C., Liu, L. and Wang,
- Z., 2017. Orthorhombic full-waveform inversion for imaging the Luda field using
- wide-azimuth ocean-bottom-cable data. The Leading Edge, 36: 75-80.
- Zhang, HJ. and Zhang, Y., 2011. Reverse time migration in vertical and tilted
- orthorhombic media. Expanded Abstr., 81st Ann. Internat. SEG Mtg., San
- Antonio: 185-189.
- Zhang, Y. and Zhang, G., 2009. One-step extrapolation method for reverse time
- migration. Geophysics, 74(4): A29-A33.
- Zhou, H., Zhang, G. and Bloor, R., 2006a. An anisotropic acoustic wave equation for
- modeling and migration in 2D TTI media. Expanded Abstr., 76th Ann. Internat.
- SEG Mtg., New Orleans: 194-198.
- Zhou, H., Zhang, G. and Bloor, R., 2006b. An anisotropic acoustic wave equation for VTI
- media. Extended Abstr., 68th EAGE Conf., Vienna: 1-5.