Multi-component seismic-resolution analysis using finite-difference acquisition modelling*
Shaun Strong 1 2 Steve Hearn 1 21 Velseis Pty Ltd and University of Queensland, Brisbane, Australia.
2 Corresponding authors. Email: sstrong@velseis.com; shearn@velseis.com
Exploration Geophysics 39(4) 189-197 https://doi.org/10.1071/EG08023
Submitted: 7 April 2008 Published: 15 December 2008
Abstract
Various rules-of-thumb (e.g. Fresnel radius, Rayleigh limit) are commonly used to predict seismic resolution, based on the dominant frequency on the image. However, seismic resolution ultimately depends on more fundamental parameters including survey design, source bandwidth, geology, and data processing. A more instructive analysis is possible via numerical modelling of the acquisition process. Here we demonstrate the improved insight available with this approach, using examples taken from the petroleum and coal sectors.
We use viscoelastic finite-difference modelling to simulate 2D multi-component acquisition sequences. The ability to allow for anelastic attenuation is important as it permits a more realistic comparison of the resolution achievable on P-wave and converted-wave (PS) imagery.
An examination of vertical resolution for a wedge model on a petroleum scale indicates that processed P-wave sections have poorer resolution (62 m) than predicted by the Widess (20 m) and Rayleigh (40 m) resolution limits. For this model the vertical resolution for the PS data is comparable to that of the P-wave data. This is in agreement with the theoretical relative-resolution relationship.
A second example examines detection of lens-like features at petroleum depth. The resolving ability on the P-wave imagery is broadly consistent with analytical predictions appropriate to migrated data (100 m laterally and 40 m vertically). Again PS resolution is comparable to P resolution.
Analysis of a typical coal target suggests that barren-zones of width 5–10 m can be resolved. The interplay of wavelength and attenuation is such that the PS image is likely to exhibit comparable, or slightly reduced, lateral resolution, provided statics are not a problem. Resolution can be downgraded significantly if statics are more severe, and in practice this is likely to have greater impact on the PS image.
Realistic numerical modelling, simulating the full acquisition and processing sequence, leads to a more pragmatic understanding of seismic resolution issues. It is a valuable tool for survey planning and image interpretation.
Key words: modelling, multi-component, resolution.
Acknowledgments
This research makes use of the processing framework provided by Seismic Unix from Colorado School of Mines. We thank Randall Taylor and Binzhong Zhou for constructive reviews, and Lindsay Thomas for editorial suggestions.
Chen, J., and Schuster, G. T., 1999, Resolution limits of migrated images: Geophysics 64, 1046–1053.
| Crossref | GoogleScholarGoogle Scholar |
Okaya, D. A., 1995, Spectral properties of the earth’s contribution to seismic resolution: Geophysics 60, 241–251.
| Crossref | GoogleScholarGoogle Scholar |
Toksöz, M. N., Johnston, D. H., and Timur, A., 1979, Attenuation of seismic waves in dry and saturated rocks: 1. Laboratory measurements: Geophysics 44, 681–690.
| Crossref | GoogleScholarGoogle Scholar |
Widess, M. B., 1973, How thin is a thin bed? Geophysics 38, 1176–1180.
| Crossref | GoogleScholarGoogle Scholar | we can define the relative resolution as
A ratio less than unity implies that the PS image has better resolution, and vice versa. Velseis (2003) provides additional comment on this topic, including practical examples.