Typical 3-D processing sequence

Excerpted from "3-D Seismic Interpretation" by M. Bacon, R. Simm and T. Redshaw
====================================================================

1. Reformat
2. Designature
3. Resample from 2 to 4 ms
4. Low cut filter (5/12 minimum phase filter)
5. Remove bad traces
6. Merge navigation with seismic headers
7. Spatial resampling from 12.5 to 25 m groups. Normal moveout (NMO) correction, K-filter, trace drop, inverse NMO
8. Spherical divergence gain corrections
9. Deconvolution before stack
10. Shot interpolation to double fold in CMP gathers
11. Radon demultiple on interpolated gathers
12. High-frequency noise removal
13. Drop of interpolated traces
14. Flex binning increasing from 37.5 m on near to 50 m on far offsets
15. Sort to common offset
16. Dip moveout (including approximated NMO) halving the number of offset planes on output
17. Pre-stack time migration using constant velocity 1600 m/s
18. Inverse NMO
19. Re-pick velocities (0.5 km grid)
20. NMO

Processing hereafter continuing on three volumes:
21. Stack to generate three volumes: near offset stack, far offset stack and full offset stack
22. 3-D constant velocity inverse time migration
23. Bulk static (gun and cable correction)
24. K-notch filter to remove pattern caused by the acquisition
25. FXY deconvolution
26. FXY interpolation to 12.5 m X 12.5 m bin grid
27. Pre-migration data conditioning (e.g. amplitude balance, edge tapers, etc.)
28. One pass steep dip 3-D time migration (using time and spatially varying velocity field)
29. Zero-phase conversion by matching to wells
30. Spectral equalisation
31. Bandpass filter
32. Residual scaling

Seismic Noise Suppression

Suppression

In seismic acquisition and processing, the attenuation of amplitudes to reduce the effects of noise or to prevent overload from the high energy of first breaks.

1. Alias filter: A filter, or a set of limits used to eliminate unwanted portions of the spectra of the seismic data, to remove frequencies that might cause aliasing during the process of sampling an analog signal during acquisition or when the sample rate of digital data is being decreased during seismic processing. (Synonyms: antialias filter)
2. Coherence filtering: A technique for removing noise and emphasizing coherent events from multiple channels of seismic data.
3. Mute: To remove the contribution of selected seismic traces in a stack to minimize air waves, ground roll and other early-arriving noise. Low-frequency traces and long-offset traces are typical targets for muting.
4. Tail mute: A cutoff in time, offset or both that has the effect of eliminating some types of noise from seismic data. A tail mute can be used to exclude slow surface waves such as ground roll.

Seismic Noise and Distortion

Noise
Anything other than desired signal. Noise includes disturbances in seismic data caused by any unwanted seismic energy, such as shot generation ground roll, surface waves, multiples, effects of weather and human activity, or random occurrences in the Earth. Noise can be minimized by using source and receiver arrays, generating minimal noise during acquisition and by filtering and stacking data during processing.

1. Aliasing: The distortion of frequency introduced by inadequately sampling a signal, which results in ambiguity between signal and noise. Aliasing can be avoided by sampling at least twice the highest frequency of the waveform or by filtering frequencies above the Nyquist frequency, the highest frequency that can be defined accurately by that sampling interval.
2. Bubble effect: Bubble pulses or bubble noise that affect data quality. In marine seismic acquisition, the gas bubble produced by an air gun oscillates and generates subsequent pulses that cause source-generated noise. Careful use of multiple air guns can cause destructive interference of bubble pulses and alleviate the bubble effect. A cage, or a steel enclosure surrounding a seismic source, can be used to dissipate energy and reduce the bubble effect.
3. Coherent noise: Undesirable seismic energy that shows consistent phase from trace to trace, such as ground roll and multiples.
4. Cultural noise: Undesirable energy, or noise, generated by human activity, such as automobile traffic that interferes with seismic surveying, or electrical power lines or the steel in pipelines that can adversely affect electromagnetic methods.
5. Ground roll: A type of coherent noise generated by a surface wave, typically a low-velocity, low-frequency, high-amplitude Rayleigh wave. Ground roll can obscure signal and degrade overall data quality, but can be alleviated through careful selection of source and geophone arrays, filters and stacking parameters.
6. Surface wave: A wave that propagates at the interface between two media as opposed to through a medium. A surface wave can travel at the interface between the Earth and air, or the Earth and water. Love waves and Rayleigh waves are surface waves.
7. Multiple reflection: Multiply reflected seismic energy, or any event in seismic data that has incurred more than one reflection in its travel path. Depending on their time delay from the primary events with which they are associated, multiples are characterized as short-path or peg-leg, implying that they interfere with the primary reflection, or long-path, where they appear as separate events. Multiples from the water bottom (the interface of the base of water and the rock or sediment beneath it) and the air-water interface are common in marine seismic data, and are suppressed by seismic processing. (Synonyms: secondary reflection)

-------------------------------------------------------------------------------------------------

Multiply-reflected seismic energy occurs in several ways but is typically removed by seismic processing. Long-path multiples appear as distinct events and generally originate deep in the subsurface. Short-path multiples are added to primary reflections and tend to come from shallow subsurface phenomena.

Diagram of ghosts and multiple reflections

Multiply-reflected seismic energy from the water bottom is common in marine seismic data, but, like many multiples, seismic processing attempts to minimize its presence.

Diagram of multiple reflections

Multiple can be divided into three major types:

1. Short-path multiple: Multiply-reflected seismic energy with a shorter travel path than long-path multiples. Short-path multiples tend to come from shallow subsurface phenomena or highly cyclical sedimentation and arrive soon after, and sometimes very near, the primary reflections. Short-path multiples are less obvious than most long-path multiples and are less easily removed by seismic processing.
2. Long-path multiple: A type of multiply-reflected seismic energy that appears as an event. Long-path multiples generate distinct events because their travel path is much longer than primary reflections giving rise to them. They typically can be removed by seismic processing. Long-path multiples originate deep in the subsurface and generate distinct events because their travel path is longer than primary reflections from the same depth.
   
   
3. Peg-leg multiple: A type of short-path multiple, or multiply-reflected seismic energy, having an asymmetric path. Short-path multiples are added to primary reflections, tend to come from shallow subsurface phenomena and highly cyclical deposition, and can be suppressed by seismic processing. In some cases, the period of the peg-leg multiple is so brief that it interferes with primary reflections, and its interference causes a loss of high frequencies in the wavelet.

Comparison of forward modeling and inversion

Seismic inversion begins with a seismic trace and outputs an acoustic impedance model.

A mathematical process by which data are used to generate a model that is consistent with the data, the process of solving the inverse problem. In seismology, surface seismic data, vertical seismic profiles and well log data can be used to perform inversion, the result of which is a model of Earth layers and their thickness, density and P- and S-wave velocities. Successful seismic inversion usually requires a high signal-to-noise ratio and a large bandwidth.

Wave

Wave:
A periodic vibrational disturbance in which energy is propagated through or on the surface of a medium without translation of the material. Waves can be differentiated by their frequency, amplitude, wavelength and speed of propagation.

Body wave:
A wave that propagates through a medium rather than along an interface. P-waves and S-waves are examples of body waves.

1. P-wave: An elastic body wave or sound wave in which particles oscillate in the direction the wave propagates. P-waves are the waves studied in conventional seismic data. P-waves incident on an interface at other than normal incidence can produce reflected and transmitted S-waves, in that case known as converted waves. (Synonyms: acoustic wave, compressional wave, dilatational wave)
2. S-wave: An elastic body wave in which particles oscillate perpendicular to the direction in which the wave propagates. S-waves are generated by most land seismic sources, but not by air guns. P-waves that impinge on an interface at non-normal incidence can produce S-waves, which in that case are known as converted waves. S-waves can likewise be converted to P-waves. S-waves, or shear waves, travel more slowly than P-waves and cannot travel through fluids because fluids do not support shear. Recording of S-waves requires receivers coupled to the solid Earth. Interpretation of S-waves can allow determination of rock properties such as fracture density and orientation, Poisson's ratio and rock type by crossplotting P-wave and S-wave velocities, and by other techniques. (Synonyms: shear wave, tangential wave)
3. SH-wave: A shear wave that is polarized so that its particle motion and direction of propagation are contained in a horizontal plane.
4. SV-wave: A shear wave that is polarized so that its particle motion and direction of propagation occur in a vertical plane.
5. Converted wave: A seismic wave that changes from a P-wave to an S-wave, or vice versa, when it encounters an interface.

Surface wave:
A wave that propagates at the interface between two media as opposed to through a medium. A surface wave can travel at the interface between the Earth and air, or the Earth and water. Love waves and Rayleigh waves are surface waves.

1. Love wave: A type of surface wave in which particles oscillate horizontally and perpendicularly to the direction of wave propagation.
2. Rayleigh wave: A type of surface wave in which particles move in an elliptical path within the vertical plane containing the direction of wave propagation. At the top of the elliptical path, particles travel opposite to the direction of propagation, and at the bottom of the path they travel in the direction of propagation. Because Rayleigh waves are dispersive, with different wavelengths traveling at different velocities, they are useful in evaluation of velocity variation with depth. Rayleigh waves make up most of the energy recorded as ground roll.
   
3. Stoneley wave: A type of large-amplitude interface, or surface, wave generated by a sonic tool in a borehole. Stoneley waves can propagate along a solid-fluid interface, such as along the walls of a fluid-filled borehole and are the main low-frequency component of signal generated by sonic sources in boreholes. Analysis of Stoneley waves can allow estimation of the locations of fractures and permeability of the formation. Stoneley waves are a major source of noise in vertical seismic profiles.
4. Tube wave: (1) A Stoneley wave that occurs at the low frequencies of seismic data. (2) An interface wave that occurs in cased wellbores when a Rayleigh wave encounters a wellbore and perturbs the fluid in the wellbore. The tube wave travels down the wellbore along the interface between the fluid in the wellbore and the wall of the wellbore. A tube wave suffers little energy loss and typically retains a very high amplitude which interferes with reflected arrivals occurring later in time on vertical seismic profile (VSP) data. Because the tube wave is coupled to the formation through which it is traveling, it can perturb the formation across open fractures intersecting the borehole. This squeezing effect can generate secondary tube waves which travel both up and down from the fracture location. Such events can be diagnostic of the presence of open fractures and their amplitude related qualitatively to the length and width, e.g., volume of the fluid-filled fracture space. This effect is generally seen only in shallow formations where the overburden pressure is lower.
5. Dispersion: A type of distortion of a wave train in which the velocity of the wave varies with frequency. Surface waves and electromagnetic body waves typically exhibit dispersion, whereas P-waves in most rocks show little change in velocity with frequency.
6. Ground roll: A type of coherent noise generated by a surface wave, typically a low-velocity, low-frequency, high-amplitude Rayleigh wave. Ground roll can obscure signal and degrade overall data quality, but can be alleviated through careful selection of source and geophone arrays, filters and stacking parameters.

Sources of free exploration-geophysical software

• Madagascar (A GPL software package with a separated header / data format very similar to classic SEPlib, but rewritten from scratch. Madagascar includes an environment for creating and testing reproducible software experiments.)
• CWP/SU Seismic Un*x Home Page (A free seismic processing system that uses an expanded SEGY-like data format.)
• qiWorkbench (A processing environment originally from BHP Billiton that includes the bhpViewer as one component.)
• Software from the Stanford Exploration Project (SEPlib, a data-cube-based seismic processing system, along with other miscellaneous software.)
• The FreeUSP home page (An Amoco/BP-heritage processing system that uses an expanded SEGY-like data format.)
• The Data Dictionary System (DDS) home page (An Amoco/BP-heritage processing system supporting high-dimensional flexible-format data files.)
• ARCO heritage source code (The Arco SPARC processing system, and some reflectivity modeling codes.)
• Computers and Geosciences journal software archive

干扰波的类型和特点

根据干扰波的出现规律，可以分为规则干扰和无规则干扰（随机干扰）两大类。

无规则干扰主要指没有一定频率，也没有一定的传播方向的波。在记录上形成杂乱无章的干扰背景。

规则干扰波是指一定的主频和一定视速度的干扰波。例如面波、声波、浅层折射波、侧面波等。它们的主要特点如下:

1. 面波(surface wave)：地震勘探中遇到的面波，它的特点是频率低，一般为几Hz到30Hz；速度低，一般为100米/秒~1000米/秒，以200米/秒～500米/秒最为常见。面波时距曲线是直线，因此在小排列(100~150米)的波形记录上面波同相轴是直的。面波随着传播距离的增大，振动延续时间也越长，形成“扫帚状”，即发生频散。随着远离爆炸点，面波的强相位逐渐向后或向前转换。强相位追踪不远，一般只有几十米。面波能量的强弱与激发岩性、激发深度以及表层地震地质条件有关。这是因为在淤泥、厚黄土及沙漠地区，由于对有效波的能量强烈吸收，有效波能量减弱，面波能量相对增强；在疏松的低速岩层中激发或所用炸药量过大造成激发频率降低，使面波能量相对加强；爆炸井较深时面波减弱，井较浅时面波增强。妥善选择激发条件和组合是克服面波的主要办法。

2. 声波(acoustic wave)：在坑中，浅水池中、河中和干井中爆炸，都会出现强烈声波。声波是空气中传播的弹性波，速度为340米/秒左右，比较稳定，频率较高，延续时间较短，呈窄带出现。为了避免声波干扰，应尽量不在浅水及浅井中放炮，尽量采用井中爆炸，并用埋井的方法以增强有效波的能量和防止声波干扰。在山区工作时，有时还会碰到多次声波的干扰。

3. 浅层折射波：当表层存在高速层，或第四系下面的老地层埋藏浅，可能观测到同相轴为直线的浅层折射波。

4. 侧面波：在地表条件比较复杂的地区进行地震勘探工作时，还会出现一种被称为侧面波的干扰波。例如在黄土高原地区，由于水系切割，形成谷沟交错的复杂地形。黄土高原的侧面是沟，原和沟的相对高差达几百米，在原与沟的交界为陡峻的黄土与空气的接触面，形成一个强波阻抗分界面，因而地震波激发后，传播到黄土边沿，被反射回来，记录上可能出现来自不同方向的具有不同视速度的干扰波。这种干扰波是一种侧面波。

还有一类是次生的低速干扰和次生的高速干扰。它和随机干扰不同，但与上面谈的规则干扰相比又显得更复杂一些，所以单独列为一类。这类干扰在频率域中与有效波是不能分离的，在视速度和视波长域中次生高速干扰又和反射有效波部分重合。它们在记录上出现的位置可以占全部记录的任意一个角落。它们是地表附近各种地物障碍物（沟、坝、公路、树木、电杆、房屋建筑物以及小山包等等）以及地表岩性不均匀性所造成的，是由反射波到达地面后，使地面产生振动，地面上任何不均匀性和地物障碍物受激发而等于地面做“敲击”动作。于是在近处产生次生的直达波和面波，在远处产生次生的折射波。这些波种类繁多，个数太多且来自于不同方向，因此就造成了极大的复杂性。

为了识别和区分这些干扰波与有效波，可以抓住它们两者之间的一个主要差异，即干扰波的最大真速度和有效波的视速度范围不同。干扰波是沿着地表附近传播的；有效波是从地下垂直来到地面的，这是它们的区别。当然，当次生干扰源相对于测线的位置不同时，次生干扰波到达测线上的视速度可以等于或大于其真速度以至无穷大，但它总会在地面的某个方向上暴露出它传播的真速度。

最后，对于干扰波的类型作一个小结：干扰波可以分为规则干扰和随机干扰两大类，规则干扰又可分为（1）沿着水平方向传播的，如面波和车辆等引起的干扰；沿着垂直方向传播的，如多次波；（2）具有重复性的，如面波；不具有重复性的，如一些人为因素产生的干扰。随机干扰也可以分为能重复出现的，如地表不均匀性引起的散射；不重复出现的，如风吹草动等自然因素引起的随机干扰。