Semblance analysis

Last updated March 06, 2024

Semblance analysis is a process used in the refinement and study of seismic data. The use of this technique along with other methods makes it possible to greatly increase the resolution of the data despite the presence of background noise. The new data received following the semblance analysis is usually easier to interpret when trying to deduce the underground structure of an area. Weighted semblance can be used for increasing the resolution of traditional semblance or make traditional semblance capable of analyzing more complicated seismic data.^[1]^[2]^[3]

History

Semblance analysis is a technique that first began to be developed and used in the late 1960s. Prior to the discovery of this method, identifying the main reflections produced by the many layers under the ground was fairly difficult. The primary reflections of these layers were often obscured by the background noise as well as noise from the many secondary reflections that are produced. The use of semblance analysis allows for the removal of the extra noise and leaves only the primary reflection.

Process

Uncorrected data showing hyperbolic curve. KR-example.png — Uncorrected data showing hyperbolic curve.

Semblance analysis allows for the refinement of seismic data. This is done by developing a velocity spectra ^[4] display to determine the velocity through different layers at depth.^[5] The easiest way to accomplish this is by recording the normal incidence path (NIP). The NIP is where you have the shot and the geophone in the same location and the path taken by the recorded sound waves is perpendicular to the boundaries between the layers.^[6] This path represents the shortest amount of time that can be taken to reach a layer and return. With this information it becomes fairly easy to calculate the velocity of the waves as they travel through each layer by using the equation for the root mean square velocity starting with the top layer and working downward.

V_{\mathrm {rms} }={\sqrt {\frac {\sum t_{i}\cdot {V_{i}}^{2}}{\sum t_{i}}}}

Once all of the velocities for the layers are known then it is possible to calculate the time needed for the wave to travel the distance down to the midpoint between each geophone and the shot point for each of the layers. As the geophones are farther away from the shot, the more the time taken for the wave to travel there increases, this forms a hyperbola in a graph of time vs. distance. The velocity data is used to correct the curves of the hyperbolas and create a flat line where all points are at an equal depth. The final step for the semblance analysis is to sum all of the data that has been corrected for velocity. This is done with the use of a computer filter to sum together all of the events that the traces share, then remove the ones they don't. The result is a single data set that has all of the primary peaks strongly displayed with most of the noise removed.^[7]

Problems

While this technique can be very useful in the analysis there are several situations in which it will not work. Semblance analysis will not work properly when the offset from the shot is greater than the depth of the reflecting layers because the data no longer has a hyperbolic pattern. To correct this it is necessary to use more complex equations that model non-hyperbolic moveout.^[8] Also in situations where there is large offset there can also be polarity reversals with moveout then the data will be heavily distorted. To make moveout analysis suitable for data with polarity reversals a method known as AK semblance developed. This method first worked only for 2D models but has since been further refined for 3D as well.^[9]

Related Research Articles

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.

Noise reduction is the process of removing noise from a signal. Noise reduction techniques exist for audio and images. Noise reduction algorithms may distort the signal to some degree. Noise rejection is the ability of a circuit to isolate an undesired signal component from the desired signal component, as with common-mode rejection ratio.

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.

In geophysics and reflection seismology, amplitude versus offset (AVO) or amplitude variation with offset is the general term for referring to the dependency of the seismic attribute, amplitude, with the distance between the source and receiver. AVO analysis is a technique that geophysicists can execute on seismic data to determine a rock's fluid content, porosity, density or seismic velocity, shear wave information, fluid indicators.

In geophysics, vertical seismic profile (VSP) is a technique of seismic measurements used for correlation with surface seismic data. The defining characteristic of a VSP is that either the energy source, or the detectors are in a borehole. In the most common type of VSP, hydrophones, or more often geophones or accelerometers, in the borehole record reflected seismic energy originating from a seismic source at the surface.

Exploration geophysics is an applied branch of geophysics and economic geology, which uses physical methods at the surface of the Earth, such as seismic, gravitational, magnetic, electrical and electromagnetic, 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. It can also be used to detect the presence of unexploded ordnance.

Seismic anisotropy is the directional dependence of the velocity of seismic waves in a medium (rock) within the Earth.

Geophysical imaging is a minimally destructive geophysical technique that investigates the subsurface of a terrestrial planet. Geophysical imaging is a noninvasive imaging technique with a high parametrical and spatio-temporal resolution. It can be used to model a surface or object understudy in 2D or 3D as well as monitor changes.

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.

In reflection seismology, normal moveout (NMO) describes the effect that the distance between a seismic source and a receiver has on the arrival time of a reflection in the form of an increase of time with offset. The relationship between arrival time and offset is hyperbolic and it is the principal criterion that a geophysicist uses to decide whether an event is a reflection or not. It is distinguished from dip moveout (DMO), the systematic change in arrival time due to a dipping layer.

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 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. 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. Synthetic seismograms are generated using specialized geophysical software.

In reflection seismology, a seismic attribute is a quantity extracted or derived from seismic data that can be analysed in order to enhance information that might be more subtle in a traditional seismic image, leading to a better geological or geophysical interpretation of the data. Examples of seismic attributes can include measured time, amplitude, frequency and attenuation, in addition to combinations of these. Most seismic attributes are post-stack, but those that use CMP gathers, such as amplitude versus offset (AVO), must be analysed pre-stack. They can be measured along a single seismic trace or across multiple traces within a defined window.

In seismology, first-break picking is the detecting or picking the onset arrivals of refracted signals from all the signals received by receiver arrays and produced by a particular source signal generation. It is also called first arrival picking or first break detection. First-break picking can be done automatically, manually or as a combination of both. With the development of computer science and the size of seismic surveys, automatic picking is often preferred.

Interferometry examines the general interference phenomena between pairs of signals in order to gain useful information about the subsurface. Seismic interferometry (SI) utilizes the crosscorrelation of signal pairs to reconstruct the impulse response of a given media. Papers by Keiiti Aki (1957), Géza Kunetz, and Jon Claerbout (1968) helped develop the technique for seismic applications and provided the framework upon which modern theory is based.

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.

During seismic exploration, P-waves penetrate down into the earth. Due to mode conversion, a P-wave can reflect upwards as an S-wave when it hits an interface. Other P-wave to S-wave (P-S) conversions can occur, but the down-up conversion is the primary focus. Unlike P-waves, converted shear waves are largely unaffected by fluids. By analyzing the original and converted waves, seismologists obtain additional subsurface information, especially due to (1) differential velocity (V_P/V_S), (2) asymmetry in the waves' angles of incidence and reflection and (3) amplitude variations.

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.

Reinhard Bortfeld was a German geophysicist.

Seismic data acquisition is the first of the three distinct stages of seismic exploration, the other two being seismic data processing and seismic interpretation. Seismic acquisition requires the use of a seismic source at specified locations for a seismic survey, and the energy that travels within the subsurface as seismic waves generated by the source gets recorded at specified locations on the surface by what is known as receivers.

References

↑ Chen, Yangkang; Liu, Tingting; Chen, Xiaohong (July–August 2015). "Velocity analysis using similarity-weighted semblance". Geophysics. 80 (4): A75–A82. Bibcode:2015Geop...80A..75C. doi:10.1190/geo2014-0618.1.
↑ Gan, Shuwei; Wang, Shoudong; Chen, Yangkang; Qu, Shan; Zu, Shaohuan (2016). "Velocity analysis of simultaneous-source data using high-resolution semblance—coping with the strong noise". Geophysical Journal International. 204 (2): 768–779. doi: 10.1093/gji/ggv484 .
↑ Ebrahimi, Saleh; Kahoo, Amin; Chen, Yangkang; Porsani, Milton (March–April 2017). "A high-resolution weighted AB semblance for dealing with amplitude-variation-with-offset phenomenon". Geophysics. 82 (2): V85–V93. Bibcode:2017Geop...82...85E. doi:10.1190/geo2016-0047.1.
↑ "Velocity Spectra". Stanford Exploration Project. Stanford University. 2000.
↑ M. TURHAN TANER; FULTON KOEHLERS; Society of Exploration Geophysicists (December 1969). "VELOCITY SPECTRA-DIGITAL COMPUTER DERIVATION AND APPLICATIONS OF VELOCITY FUNCTIONS" (PDF). GEOPHYSICS. p. 859.
↑ Schlumberger Limited (2020). "Oilfield Glossary: Normal Incidence".
↑ Taner, Turhan; Fulton Koehler (December 1969). "Velocity Spectra-Digital Computer Derivation and Application of Velocity Functions". Geophysics. 34 (6): 859–881. doi:10.1190/1.1440058.
↑ Sergey Fomel; Alexey Stovas (2010). "Generalized nonhyperbolic moveout approximation". Geophysics. Society of Exploration Geophysicists. 75 (2): U9–U18. Bibcode:2010Geop...75....9F. doi:10.1190/1.3334323. eISSN 1942-2156. ISSN 0016-8033 – via SEG Library.
↑ Yan, Jia; Ilya Tsvankin (March–April 2008). "AVO-sensitive semblance analysis for wide-azimuth data". Geophysics. 73 (2): U1–U11. doi:10.1190/1.2834115.

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.

[vscan-1] Chen, Yangkang; Liu, Tingting; Chen, Xiaohong (July–August 2015). "Velocity analysis using similarity-weighted semblance". Geophysics. 80 (4): A75–A82. Bibcode:2015Geop...80A..75C. doi:10.1190/geo2014-0618.1.

[svscan-2] Gan, Shuwei; Wang, Shoudong; Chen, Yangkang; Qu, Shan; Zu, Shaohuan (2016). "Velocity analysis of simultaneous-source data using high-resolution semblance—coping with the strong noise". Geophysical Journal International. 204 (2): 768–779. doi: 10.1093/gji/ggv484 .

[abvscan-3] Ebrahimi, Saleh; Kahoo, Amin; Chen, Yangkang; Porsani, Milton (March–April 2017). "A high-resolution weighted AB semblance for dealing with amplitude-variation-with-offset phenomenon". Geophysics. 82 (2): V85–V93. Bibcode:2017Geop...82...85E. doi:10.1190/geo2016-0047.1.

[4] "Velocity Spectra". Stanford Exploration Project. Stanford University. 2000.

[5] M. TURHAN TANER; FULTON KOEHLERS; Society of Exploration Geophysicists (December 1969). "VELOCITY SPECTRA-DIGITAL COMPUTER DERIVATION AND APPLICATIONS OF VELOCITY FUNCTIONS" (PDF). GEOPHYSICS. p. 859.

[6] Schlumberger Limited (2020). "Oilfield Glossary: Normal Incidence".

[7] Taner, Turhan; Fulton Koehler (December 1969). "Velocity Spectra-Digital Computer Derivation and Application of Velocity Functions". Geophysics. 34 (6): 859–881. doi:10.1190/1.1440058.

[8] Sergey Fomel; Alexey Stovas (2010). "Generalized nonhyperbolic moveout approximation". Geophysics. Society of Exploration Geophysicists. 75 (2): U9–U18. Bibcode:2010Geop...75....9F. doi:10.1190/1.3334323. eISSN 1942-2156. ISSN 0016-8033 – via SEG Library.

[9] Yan, Jia; Ilya Tsvankin (March–April 2008). "AVO-sensitive semblance analysis for wide-azimuth data". Geophysics. 73 (2): U1–U11. doi:10.1190/1.2834115.

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]