Synthetic seismogram

Last updated

A synthetic seismogram is the result of forward modelling the seismic response of an input earth model, which is defined in terms of 1D, 2D or 3D variations in physical properties. In hydrocarbon exploration this is used to provide a 'tie' between changes in rock properties in a borehole and seismic reflection data at the same location. It can also be used either to test possible interpretation models for 2D and 3D seismic data or to model the response of the predicted geology as an aid to planning a seismic reflection survey. In the processing of wide-angle reflection and refraction (WARR) data, synthetic seismograms are used to further constrain the results of seismic tomography. [1] In earthquake seismology, synthetic seismograms are used either to match the predicted effects of a particular earthquake source fault model with observed seismometer records or to help constrain the Earth's velocity structure. [2] Synthetic seismograms are generated using specialized geophysical software.

Contents

1D synthetics

Seismic reflection data are initially only available in the time domain. In order that the geology encountered in a borehole can be tied to the seismic data, a 1D synthetic seismogram is generated. This is important in identifying the origin of seismic reflections seen on the seismic data. Density and velocity data are routinely measured down the borehole using wireline logging tools. These logs provide data with a sampling interval much smaller than the vertical resolution of the seismic data. The logs are therefore often averaged over intervals to produce what is known as a 'blocked-log'. [3] This information is then used to calculate the variation in acoustic impedance down the well bore using the Zoeppritz equations. [4] This acoustic impedance log is combined with the velocity data to generate a reflection coefficient series in time. This series is convolved with a seismic wavelet to produce the synthetic seismogram. The input seismic wavelet is chosen to match as closely as possible to that produced during the original seismic acquisition, paying particular attention to phase and frequency content.

1.5D seismic modelling

The convolutional 1D modelling produces seismograms containing approximations of primary reflections only. For more accurate modelling involving multiple reflections, head waves, guided waves and surface waves, as well as transmission effects and geometrical spreading, full waveform modelling is required. For 1D elastic models the most accurate approach to full waveform modelling is known as the reflectivity method. [5] This method is based on the integral transform approach, whereby the wave field (cylindrical or spherical wave) is represented by a sum (integral) of time-harmonic plane waves. [6] The reflection and transmission coefficients for individual plane waves propagating in a stack of layers can be computed analytically using a variety of methods, such as matrix propagator, [7] [8] [9] [10] [11] global matrix [12] or invariant embedding. [13] This group of methods is called 1.5D because the earth is represented by a 1D model (flat layers), while wave propagation is considered either in 2D (cylindrical waves) or 3D (spherical waves).

2D synthetic seismic modeling

A similar approach can be used to examine the seismic response of a 2D geological cross-section. This can be used to look at such things as the resolution of thin beds or the different responses of various fluids, e.g. oil, gas or brine in a potential reservoir sand. [14] It may also be used to test out different geometries of structures such as salt diapirs, to see which gives the best match to the original seismic data. A cross-section is built with density and seismic velocities assigned to each of the individual layers. These can be either constant within a layer or varying in a systematic fashion across the model both horizontally and vertically. The software program then runs a synthetic acquisition across the model to produce a set of 'shot gathers' that can be processed as if they were real seismic data to produce a synthetic 2D seismic section. The synthetic record is generated using either a ray-tracing algorithm or some form of full waveform modelling, depending on the purpose of the modelling. Ray-tracing is quick and sufficient for testing the illumination of the structure, [15] but full waveform modelling will be necessary to accurately model the amplitude response. [16]

3D synthetic seismic modelling

The approach can be further expanded to model the response of a 3D geological model. This is used to reduce the uncertainty in interpretation by modelling the response of the 3D model to a synthetic seismic acquisition that matches as closely as possible to that actually used in acquiring the data that has been interpreted. [17] The synthetic seismic data is then processed using the same sequence as that used for the original data. This method can be used to model both 2D and 3D seismic data that has been acquired over the area of the geological model. During the planning of a seismic survey, 3D modelling can be used to test the effect of variation in seismic acquisition parameters, such as the shooting direction or the maximum offset between source and receiver, on the imaging of a particular geological structure. [18] [19]

WARR data modelling

Wide Aperture Reflection and Refraction (WARR) models' initial processing is normally carried out using a tomographic approach in which the time of observed first arrivals is matched by varying the velocity structure of the subsurface. The model can be further refined using forward modelling to generate synthetic seismograms for individual shot gathers. [1]

Earthquake modelling

Source modelling

In areas that have a well understood velocity structure it is possible to use synthetic seismograms to test out the estimated source parameters of an earthquake. Parameters such as the fault plane, slip vector and rupture velocity can be varied to produce synthetic seismic responses at individual seismometers for comparison with the observed seismograms. [20]

Velocity modelling

For seismic events of known type and location, it is possible to obtain detailed information about the Earth's structure, at various scales, by modelling the teleseismic response of the event. [2]

Related Research Articles

<span class="mw-page-title-main">Seismology</span> Scientific study of earthquakes and propagation of elastic waves through a planet

Seismology is the scientific study of earthquakes and the propagation of elastic waves through the Earth or through other planet-like bodies. It also includes studies of earthquake environmental effects such as tsunamis as well as diverse seismic sources such as volcanic, tectonic, glacial, fluvial, oceanic, atmospheric, and artificial processes such as explosions. A related field that uses geology to infer information regarding past earthquakes is paleoseismology. A recording of Earth motion as a function of time is called a seismogram. A seismologist is a scientist who does research in seismology.

<span class="mw-page-title-main">Seismic wave</span> Seismic, volcanic, or explosive energy that travels through Earths layers

A seismic wave is a wave of acoustic energy that travels through the Earth. It can result from an earthquake, volcanic eruption, magma movement, a large landslide, and a large man-made explosion that produces low-frequency acoustic energy. Seismic waves are studied by seismologists, who record the waves using seismometers, hydrophones, or accelerometers. Seismic waves are distinguished from seismic noise, which is persistent low-amplitude vibration arising from a variety of natural and anthropogenic sources.

<span class="mw-page-title-main">Seismogram</span> Graph output by a seismograph

A seismogram is a graph output by a seismograph. It is a record of the ground motion at a measuring station as a function of time. Seismograms typically record motions in three cartesian axes, with the z axis perpendicular to the Earth's surface and the x- and y- axes parallel to the surface. The energy measured in a seismogram may result from an earthquake or from some other source, such as an explosion. Seismograms can record many things, and record many little waves, called microseisms. These tiny microseisms can be caused by heavy traffic near the seismograph, waves hitting a beach, the wind, and any number of other ordinary things that cause some shaking of the seismograph.

Seismic tomography or seismotomography is a technique for imaging the subsurface of the Earth with seismic waves produced by earthquakes or explosions. P-, S-, and surface waves can be used for tomographic models of different resolutions based on seismic wavelength, wave source distance, and the seismograph array coverage. The data received at seismometers are used to solve an inverse problem, wherein the locations of reflection and refraction of the wave paths are determined. This solution can be used to create 3D images of velocity anomalies which may be interpreted as structural, thermal, or compositional variations. Geoscientists use these images to better understand core, mantle, and plate tectonic processes.

<span class="mw-page-title-main">Reflection seismology</span> Explore subsurface properties with seismology

Reflection seismology is a method of exploration geophysics that uses the principles of seismology to estimate the properties of the Earth's subsurface from reflected seismic waves. The method requires a controlled seismic source of energy, such as dynamite or Tovex blast, a specialized air gun or a seismic vibrator. Reflection seismology is similar to sonar and echolocation. This article is about surface seismic surveys; for vertical seismic profiles, see VSP.

Exploration geophysics is an applied branch of geophysics and economic geology, which uses physical methods, such as seismic, gravitational, magnetic, electrical and electromagnetic at the surface of the Earth to measure the physical properties of the subsurface, along with the anomalies in those properties. It is most often used to detect or infer the presence and position of economically useful geological deposits, such as ore minerals; fossil fuels and other hydrocarbons; geothermal reservoirs; and groundwater reservoirs.

<span class="mw-page-title-main">Earthscope</span> Earth science program exploring the structure of the North American continent

Earthscope was an National Science Foundation (NSF) funded earth science program that, from 2003-2018, used geological and geophysical techniques to explore the structure and evolution of the North American continent and to understand the processes controlling earthquakes and volcanoes. The project had three components: USArray, the Plate Boundary Observatory, and the San Andreas Fault Observatory at Depth. Organizations associated with the project included UNAVCO, the Incorporated Research Institutions for Seismology (IRIS), Stanford University, the United States Geological Survey (USGS) and National Aeronautics and Space Administration (NASA). Several international organizations also contributed to the initiative. EarthScope data are publicly accessible.

Seismic anisotropy is a term used in seismology to describe the directional dependence of the velocity of seismic waves in a medium (rock) within the Earth.

<span class="mw-page-title-main">Geodynamics</span> Study of dynamics of the Earth

Geodynamics is a subfield of geophysics dealing with dynamics of the Earth. It applies physics, chemistry and mathematics to the understanding of how mantle convection leads to plate tectonics and geologic phenomena such as seafloor spreading, mountain building, volcanoes, earthquakes, faulting. It also attempts to probe the internal activity by measuring magnetic fields, gravity, and seismic waves, as well as the mineralogy of rocks and their isotopic composition. Methods of geodynamics are also applied to exploration of other planets.

Seismic migration is the process by which seismic events are geometrically re-located in either space or time to the location the event occurred in the subsurface rather than the location that it was recorded at the surface, thereby creating a more accurate image of the subsurface. This process is necessary to overcome the limitations of geophysical methods imposed by areas of complex geology, such as: faults, salt bodies, folding, etc.

Sonic logging is a well logging tool that provides a formation’s interval transit time, designated as , which is a measure of a how fast elastic seismic compressional and shear waves travel through the formations. Geologically, this capacity varies with many things including lithology and rock textures, most notably decreasing with an increasing effective porosity and increasing with an increasing effective confining stress. This means that a sonic log can be used to calculate the porosity, confining stress, or pore pressure of a formation if the seismic velocity of the rock matrix, , and pore fluid, , are known, which is very useful for hydrocarbon exploration.

The receiver function technique is a way to image the structure of the Earth and its internal boundaries by using the information from teleseismic earthquakes recorded at a three-component seismograph.

In geophysics, seismic inversion is the process of transforming seismic reflection data into a quantitative rock-property description of a reservoir. Seismic inversion may be pre- or post-stack, deterministic, random or geostatistical; it typically includes other reservoir measurements such as well logs and cores.

A seismic metamaterial, is a metamaterial that is designed to counteract the adverse effects of seismic waves on artificial structures, which exist on or near the surface of the earth. Current designs of seismic metamaterials utilize configurations of boreholes, trees or proposed underground resonators to act as a large scale material. Experiments have observed both reflections and bandgap attenuation from artificially induced seismic waves. These are the first experiments to verify that seismic metamaterials can be measured for frequencies below 100 Hz, where damage from Rayleigh waves is the most harmful to artificial structures.

<span class="mw-page-title-main">Surface wave inversion</span>

Seismic inversion involves the set of methods which seismologists use to infer properties through physical measurements. Surface-wave inversion is the method by which elastic properties, density, and thickness of layers in the subsurface are obtained through analysis of surface-wave dispersion. The entire inversion process requires the gathering of seismic data, the creation of dispersion curves, and finally the inference of subsurface properties.

<span class="mw-page-title-main">Near-surface geophysics</span> Geophysics of first tens of meters below surface

Near-surface geophysics is the use of geophysical methods to investigate small-scale features in the shallow subsurface. It is closely related to applied geophysics or exploration geophysics. Methods used include seismic refraction and reflection, gravity, magnetic, electric, and electromagnetic methods. Many of these methods were developed for oil and mineral exploration but are now used for a great variety of applications, including archaeology, environmental science, forensic science, military intelligence, geotechnical investigation, treasure hunting, and hydrogeology. In addition to the practical applications, near-surface geophysics includes the study of biogeochemical cycles.

In reflection seismology, the anelastic attenuation factor, often expressed as seismic quality factor or Q, quantifies the effects of anelastic attenuation on the seismic wavelet caused by fluid movement and grain boundary friction. As a seismic wave propagates through a medium, the elastic energy associated with the wave is gradually absorbed by the medium, eventually ending up as heat energy. This is known as absorption and will eventually cause the total disappearance of the seismic wave.

<span class="mw-page-title-main">Linear seismic inversion</span> Interpretation of seismic data using linear model

Inverse modeling is a mathematical technique where the objective is to determine the physical properties of the subsurface of an earth region that has produced a given seismogram. Cooke and Schneider (1983) defined it as calculation of the earth's structure and physical parameters from some set of observed seismic data. The underlying assumption in this method is that the collected seismic data are from an earth structure that matches the cross-section computed from the inversion algorithm. Some common earth properties that are inverted for include acoustic velocity, formation and fluid densities, acoustic impedance, Poisson's ratio, formation compressibility, shear rigidity, porosity, and fluid saturation.

Michel Campillo is a French geophysicist and seismologist. He was born in 1957 and grew up in Chambon sur Lignon.

Seismic wide-angle reflection and refraction is a technique used in geophysical investigations of Earth's crust and upper mantle. It allows the development of a detailed model of seismic velocities beneath Earth's surface well beyond the reach of exploration boreholes. The velocities can then be used, often in combination with the interpretation of standard seismic reflection data and gravity data, to interpret the geology of the subsurface.

References

  1. 1 2 Makris, J., Egloff, F. & Rihm, R. 1999. WARRP (Wide Aperture Reflection and Refraction Profiling): The principle of successful data acquisition where conventional seismic fails, SEG 1999 Expanded Abstracts
  2. 1 2 Helmberger, D.V. 1974, Understanding Seismograms by Constructing Numerical Models, Engineering and Science, 38, 26–29.
  3. ↑ Goldberg, D., Wilkens, R.H. & Moos, D. 1987. Seismic modeling of diagenetic effects in Cenozoic marine sediments at Deep Sea Drilling Project sites 612 and 613, DSDP Initial report on Leg 95, 23
  4. ↑ OBartels, T., Krastel, S., and Spiess, V., 2007. Correlation of high-resolution seismic data with ODP Leg 208 borehole measurements. In Kroon, D., Zachos, J.C., and Richter, C. (Eds.), Proc. ODP, Sci. Results, 208: College Station, TX (Ocean Drilling Program), 1–27
  5. ↑ Fuchs, K., and G. Muller, 1971, Computation of synthetic seismograms with the reflectivity method and comparison with observations, Geophys. J. R. Astron. Soc, 23, 417.
  6. ↑ Aki, K. and Richards, R.G., Quantitative Seismology, Theory and Methods, Vol. I, W. H. Freeman, 1980.
  7. ↑ Thomson, W.T., 1950, Transmission of elastic waves through a stratified solid material, Journal of Applied Physics, 21, 89–93.
  8. ↑ Haskell, N. A., The dispersion of surface waves in multilayered media, , Bulletin of the Seismological Society of America, 43, 17–34,1953.
  9. ↑ Dunkin, I.W., 1965, Computation of model solutions in layered elastic media at high frequencies, Bulletin of the Seismological Society of America, 55, 335–358.
  10. ↑ Thrower, E.N., The computation of the dispersion of elastic waves in layered media, Journal of Sound and Vibration, 2, 210–226.
  11. ↑ Molotkov L.A., 1984, Matrix method in the theory of wave propagation in layered elastic and fluid media, Nauka (in Russian).
  12. ↑ Schmidt, H and Tango., 1986, Efficient global matrix approach to the computation of synthetic seismograms, Geophysical Journal of the Royal Astronomic Society, 84, pp 331–359.
  13. ↑ Kennett, B. L. N., 1985, Seismic wave propagation in stratified media, Cambridge University Press.
  14. ↑ Hodgetts, D. & Howell, J.A. 2000. Synthetic seismic modelling of a large-scale geological cross-section from the Book Cliffs, Utah, USA, Petroleum Geoscience, 6, 221–229.
  15. ↑ Graham, S., Lawton, D. & Spratt, D. 2005. Sub-thrust imaging:modelling example from the Cusiana oilfield, Llanos Basin, Colombia, CSEG National Convention, Abstract.
  16. ↑ Li, Y., Downton, J. & Xu, Y. 2004. AVO Modeling in Seismic Processingand Interpretation II. Methodologies, CSEG Recorder, January, 38–44.
  17. ↑ Gawith, D.E. & Gutteridge, P.A. 1996. Seismic validation of reservoir simulation using a shared earth model, Petroleum Geoscience, 2, 97–103.
  18. ↑ Gjøystdal, H., Iversen, E., Lecomte, I., Kaschwich, T., Drottning, Å. and Mispel, J. 2007. Improved applicability of ray tracing in seismic acquisition, imaging, and interpretation, Geophysics, 72, 261–271.
  19. ↑ Ray, A., Pfau, G. & CHen, R. 2004. Importance of ray-trace modeling in the discovery of Thunder Horse North Field, Gulf of Mexico, The Leading Edge, 23, 68–70.
  20. ↑ Cotton, F. & Campillo, M. 1994. Application of seismogram synthesis to the study of earthquake source from strong motion records, Annali di Geofisica, 37, 1539–1564.
This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.