News and Announcements
ARTICLE
Reflectivity wavefield inversion of synthetic data from elastic and viscoelastic anisotropic flat-layered, shale-bounded fractured sandstone reservoir models
SYLVIA PACHECO,
GEORGE A. MCMechAN
Show Less
Center for Lithospheric Studies, The University of Texas at Dallas, 800 W. Campbell Road, Richardson, TX 75080-3021, U.S.A.,
JSE 2013, 22(3), 223–250;
Submitted: 23 December 2012 | Accepted: 2 May 2013 | Published: 1 July 2013
© 2013 by the Authors. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution -Noncommercial 4.0 International License (CC-by the license) ( https://creativecommons.org/licenses/by-nc/4.0/ )
Abstract
Pacheco, S. and McMechan, G.A., 2013. Reflectivity wavefield inversion of synthetic data from elastic and viscoelastic anisotropic flat-layered, shale-bounded fractured sandstone reservoir models. Journal of Seismic Exploration, 22: 223-250. Full-wavefield inversion based on frequency-domain reflectivity modeling is developed and applied, in the 7-p domain, to synthetic data for flat-layered viscoelastic, anisotropic models. P- and S-wave quality factors and HTI and VTI anisotropy are extracted, by inversion, from synthetic data for viscoelastic, fractured sand reservoir models. Ignoring viscoelasticity or anisotropy results in errors and increased uncertainty in the estimates of the other seismic parameters, as the parameterization is incomplete. Increasing the number of layers or parameters, or adding noise, makes the inversion more difficult to converge to the correct solution; parameter correlations confirm this observation, showing more inter-layer correlations. The noise is mainly left in the residual wavefield. Better results are obtained with layer stripping than without, as the solution space has fewer local minima when the model contains fewer parameters (fewer degrees of freedom). Multiple layer stripping runs, solving for parameters adjacent pairs of layers together, improves the inversion performance and the accuracy of the results.
Keywords
wavefield seismic inversion
sandstones
reflectivity
7-p
viscoelastic
anisotropic
References
- Bakker, P.M., 1995. About the completeness of the classification of cases of elliptic anisotropy.Proc. Roy. Soc. London A, 451: 367-373.
- Bakulin, A., Woodward, M., Nichols, D., Osypov, K. and Zdraveva, O., 2010. Building tiltedtransversely isotropic depth models using localized anisotropic tomography with wellinformation. Geophysics, 75(4): D27-D36.
- Baraka-Lokmane, S., 2002. Hydraulic versus pneumatic measurements of fractured sandstonepermeability. J. Petrol. Sci. Eng., 36: 183-192.
- Barbosa, B., Costa, J., Gomes, E. and Schleicher, J., 2008. Resolution analysis forstereotomography in media with elliptic and analeptic anisotropy. Geophysics, 73(4): R49-R58.
- Chang, H. and McMechan, G.A., 2009a. 3D, 3-C full wavefield elastic inversion for anisotropicparameters: A feasibility study with synthetic data. Geophysics, 74(6): WCC159-WCC175.
- Chang, H. and McMechan, G.A., 2009b. Synthetic data tests of 3D full-wavefield inversion for
- P-wave anisotropic parameter estimation in flat layered VTI, HTI and orthorhombic media.J. Seismic Explor., 18: 249-270.
- Golikov, P. and Stovas, A., 2012. Traveltime parameters in a tilted elliptical anisotropic medium.Geophys. Prosp., 60: 433-443.
- Grechka, V., 2009. Nonuniqueness of traveltime inversion in elliptically anisotropic media.
- Expanded Abstr., 79th Ann. Internat. SEG Mtg., Houston, 28: 196-200.
- Guiton, M., Sassi, W., Leroy, Y. and Gauthier, B., 2003. Mechanical constraints on the chronologyof fracture activation in folded Devonian sandstone of the western Moroccan Anti-Atlas. J.Struct. Geol., 25: 1317-1330.
- Hennings, P., Olson, J. and Thompson, L., 2000. Combining outcrop data and threedimensionalstructural models to characterize fractured reservoirs: an example from Wyoming. AAPGBull., 84: 830-849.250 PACHECO & MCMECHAN
- Hu, Y. and McMechan, G.A., 2010. Sensitivity of 3-D, 3-C finite-difference elastic seismic datato inclusion parameters in HTI and TTI media with high inclusion density. Geophysics,75(2): T49-T61.
- Kelly, J., Parnell, J. and Chen, H.H., 2000. Application of fluid inclusions to studies of fracturedsandstone reservoirs. J. Geochem. Explor., 69-70: 705-709.
- Khazanehdari, J. and McCann, C., 2005. Acoustic and petrophysical relationships in low-shalesandstone reservoir rocks. Geophys. Prosp., 53: 447-461.
- Kleyn, A.H., 1956. On seismic wave propagation in anisotropic media with applications in the Betunarea, south Sumatra. Geophys. Prosp., 4: 56-69.
- Koesoemadinata, A.P. and McMechan, G.A., 2001. Empirical estimation of viscoelastic seismicparameters from petrophysical properties of sandstone. Geophysics, 66: 1457-1470.
- Koesoemadinata, A.P. and McMechan, G.A., 2003. Petro-seismic inversion for sandstoneproperties. Geophysics, 68: 1611-1625.
- Maultzsch, S., Chapman, M., Liu, E. and Li, X.Y., 2003. Modelling frequency-dependent seismicanisotropy in fluid-saturated rock with aligned fractures: implication of fracture sizeestimation from anisotropic measurements. Geophys. Prosp., 51: 381-392.
- Martinez, R.D. and McMechan, G.A., 1991a. 7-p seismic data for viscoelastic media. Part 1:Modelling. Geophys. Prosp., 39: 141-156.
- Martinez, R.D. and McMechan, G.A., 1991b. 7-p seismic data for viscoelastic media. Part 2:
- Linearized inversion. Geophys. Prosp., 39: 157-182.
- Martinez, R.D., 1993. Wave propagation effects on amplitude variation with offset measurements:A study. Geophysics, 58: 534-543.
- McCaulay, A.D., 1986. Plane-layer prestack inversion in the presence of surface reverberation.Geophysics, 51: 1789-1800.
- Mora, P., 1989. Inversion = migration + tomography. Geophysics, 54: 1575-1586.
- Odling, N., Gillespie, P., Bourgine, B., Castaing, C., Chils, J., Christensen, N., Fillion, E.,
- Genter, A., Olsen, C. Thrane, L., Trice, R., Aarseth, E., Walsh, J. and Watterson, J.,
- Variations in fracture system geometry and their implications for fluid flow infractured hydrocarbon reservoirs. Petrol. Geosc., 5: 373-384.
- Ohlsen, F. and Macbeth, C., 1999. Elliptical anisotropy - regression or advance? Expanded Abstr.,69th Ann. Internat. SEG Mtg., Houston: 1600-1603.
- Ortega, O. and Marrett, R., 2000. Prediction of macrofracture properties using microfractureinformation, Mesaverde Group sandstones, San Juan basin, New Mexico. J. Struct. Geology,22: 571-588.
- Ozdenvar, T. and McMechan, G.A., 1997. Algorithm for staggered grid computations forporoelastic, elastic, acoustic, and scalar wave equations. Geophys. Prosp., 45: 403-420.
- Pan, G.S., Phinney, R.A. and Odom, R.I., 1988. Full-waveform inversion of plane-waveseismograms in stratified acoustic media: Theory and feasibility. Geophysics, 53: 21-31.
- Pica, A., Diet, J.P. and Tarantola, A., 1990. Nonlinear inversion of seismic reflection data in alaterally invariant medium. Geophysics, 55(3): 284-292.
- Schwab, F.A. and Knopoff, L., 1972. Fast Surface waves and free mode computations. Meth.Comput. Physics, 11: 87-180.
- Stoffa, P.L. and Sen, M.K., 1992. Seismic waveform inversion using global optimization. J. SeismicExplor., 1: 9-27.
- Thomsen, L., 2002. Understanding Seismic Anisotropy in Exploration and Exploitation. SEGDistinguished Instructor Series, no. 5.
- Tiwari, U.K. and McMechan, G.A., 2007a. Effects of incomplete parameterization on inversion offull wavefield viscoelastic seismic data for petrophysical reservoir properties. Geophysics,72: 09-17.
- Tiwari, U.K. and McMechan, G.A., 2007b. Estimation of effective pressure and water saturationby viscoelastic inversion of synthetic time-lapse seismic data for a gas sandstone reservoir.J. Seismic Explor., 16: 57-68.
- Xu, T., McMechan, G.A. and Sun, R., 1995. 3-D prestack full-wavefield inversion. Geophysics,60: 1805-1818.
