Total Pageviews

Saturday, 19 March 2011

EMP

Electromagnetic Pulse Energy Source

Introduction
EMP sources based on the electrodynamic and electromagnetic principles. Charge in the cylinder lifts the plate in EMP and released as soon as the it discharges and a pulse is generated of high energy . Look at the amplitude spectrum of the EMP source, the frequency content for dynamite is taken  upto 125 Hz, but any amplitude is hardly seen after 70 Hz. But the amplitude spectrum you can very well see it goes upto 90 hz. So its quite comparable. 


Migration Methods

Seismic Data Migration Methods


1. Kirchoff Migration
2. Gaussian Beam migration

Gaussian beam migration (Hill, 1990, 2001) is an elegant, accurate, and efficient depth migration method. It has the ability to image complicated geologic structures with fidelity exceeding that of single-arrival Kirchhoff migration and approaching that of wave-equation migration. In fact, its accuracy can exceed that of most wave-equation migrations in imaging very steep dips, especially in three dimensions and especially in the presence of anisotropy.
The success of Gaussian beam migration has sparked interest in beam migrations, some based on Gaussian beams, others not. Some variations of beam migration have as their sole purpose the speedup of Kirchhoff migration (Sun et al., 2000); others use wavefield extrapolation to build migration Green’s functions (Brandsberg-Dahl, 2003); still others use wavelet decomposition of seismic wavefields to emulate the simultaneous space/wavenumber localization property of Gaussian beam migration (Wu et al., 2002).



3. Reverse Time migration





Friday, 18 March 2011

GeGraphix & Kingdom Comparison

GEOGRAPHIX SEISMIC INTERPRETATION SOFTWARE SUITE

In 1984 David Armiitage founded GeoGraphix with a goal of simplifying and automating standard geology tasks using developing computer technologies. In the very next year, he launched the industry's first automated base map generation and printing solution for the PC. In 1987 he Launched IsoMap, the industry's first PC-based surface modeling system. 1989 Released the original GeoGraphix Exploration System (GES, Package included SeisBase, Wellbase and IsoMap.)  on Windows 3.0. In June 1995 GeoGraphix was acquired by Landmark to address the technology needs of consultants and small- to mid-sized operators. It was then in July 1996 Halliburton acquired Landmark (and GeoGraphix) and later they established strategic alliance with ESRI, the industry leader in mapping and map object technologies(ESRI 16-bit technologies incorporated into GeoGraphix platform). By this time in 1998 geographix database was migrated from a flat file database system to Sybase relational ODBC database.


Then came the Geographix Discovery Suite in Aug-1999, the first suite of geoscience applications built on a shared data management and common visualization foundation that provided dynamic integration between geological and geophysical data. GeoGraphix Discovery provided unprecedented integrated interpretation capabilities on the Windows NT® desktop. All GeoGraphix geoscience interpretation solutions-SeisVision™, GESXplorer™, and PRIZM™ - were incorporated into this first GeoGraphix Discovery suite.

In 2000 Established first reseller partner in China-GNT International and later in Nov 2001 they realeased significant performance and capability enhancement to Discovery suite. Applications performed as much as 2 to 10 times faster and included 90+ enhancements. In June 2002 Wellbase was named a laureate to the 2002 Computerworld Honors Collection and then in 2005 Established GeoGraphix regional champions network.Later in Feb 2005 Halliburton acquired smartSECTION® geologic software business from A2D Technologies, significantly expanding the well log correlation and subsurface interpretation capabilities in the GeoGraphix suite of software. 
In Nov 2005 - Nov. Released GeoGraphix Discovery on OpenWorks® Database, interpreters having their project in geographix can migrate it to openworks and make use of the landmark's advance interpretation technology . 

In 2007 - 2008 established business relationship with WhiteStar Data and Spatial Energy. It was then in Oct 2009 Haliburton Launched GeoGraphix Discovery 3D software, one of the first oil and gas solutions to leverage state-of-the-art interactive visualization tools from Microsoft's Xbox gaming console and the 64-bit capabilities and performance of the new Microsoft Windows 7 operating system.

But things were different in the 21st century, many small software company grew up and gave the industry such advancement that previously used to be with unix based seismic interpretation softwares. SMT came up in a larger way with its Kingdom advance seismic interpretation suite. 

Not later than in July 2010, Landmark and LMKR signed an agreement transferring the management, sales, marketing and product development of GeoGraphix to LMKR.  Landmark remains vested in the success of the GeoGraphix business.  Customer support and training will continue to be offered through the Landmark organization. With Kingoms Advance came advance horizon auto-picker, horizon mapping, flex gridding, fulat mapping, fault modeling, atrribute analysis, post processing attributes, AVO anaysis, PostSTK Inversion, Geomodeling(Jewel suite application), Data migration(project data/ASCII data), and 3D data vizualization while interpretation, network security and project data transfer to different software through openspirit tunnel technology. Still there are lot areas which needs to make them more stronger. kingdom is stronger in geophysics and week in geology. They are rather moving towards modeling world like Petrel.
Petrel is the leading in Geomodeling. Geologist like petrel in many ways for its ability in geomodeling flexiblity. 
Like for inversion and geostatistical mapping its always been Hampson Russel at the top. Kriging is the most famous geostatistical application, which more accurate method for time to deppth/porosity/permeability mapping. Using minimum number of wells one can easily convert time to depth maps.

In 2010 GeoGraphix released version 5000.0.2.0.  New technologies released include streaming GIS, Geomodeling While Drilling (smartSTRAT), Advanced 3D Visualization and Optimized Well and Field Planning workflows. When Shale Exploration was capturing the major market in India, then geographix came back with a breakthrough of technology geosteering while drilling for horizontal wells for shale gas. Its a breakthrough where u model the path of the well, through the horizon, so that it doesn't misses the horizon and passes along the top of the horizon.
Rest everything remaining same as before seisvision, Geoatlas, Prizm, coordinate manger, well base, but with introduction of tree in modules like as in kingdoms


Geographix 

 Hardware Recommendation :
a)      Windows7/10 - 64 bit OS, Ram 8gb, cache memory > 2gb,  recommended for Discovery 3D  
b)      Windows7/10 - 32 bit OS, Ram 4gb
c)       DirectX  10  
d)      V-Ram ~  1gb
e)      .Net  Framework 3.5 sp1

OS: Discovery 3D viewer runs in Windows7/10.

Database: Oracle client/Sybase.

Discovery 3D is much advance 3D viewer: multiple surface views and blending them
 ( no interpretation done through 3D Viewer).

Seisvision; Well base; Seisbase; Geoatlas; Prizm is same as before, with few more interactive tools inbuilt.

Arc GIS capability: (transparency tool added)
a)      Live streaming satellite image (require url)
b)      Geo-Tiff images import

Smart Section

a)      Well correlation
b)      SmartSRAT using frame builder geomodeling : Model the surface
c)       Geosteering-along LWD data Shale formation
d)      Mircoseismic data analysis: manage and plan wells.

Field planning for shale exploration – GeoAtlas- Wellbase- SSL Google Earth; Discovery 3D-TrackPlanner Xpress

Attribute
a)      Seisvision has curvature attribute
b)      Semblance coherency attribute in PStax

Importing well data from excel file: SSL


Version

GeoGraphix 

Kingdom 2017 
Database:
Database: Oracle(Track Planner); Sybase® SQL®
SQL server / Express 2012 or 2014 or 2017 (Auto/Manual) /Oracle
OS:
·   Windows XP, but needs Windows 7 for Discovery 3D View Advance visualization / windows 10
·   Windows 7/10
Data Management
·    Import/export isn’t so flexible
·    Project backup has to taken before leaving the desk
·    It can handle data of different projection into a single Project.
·         Stable Database
·         Very flexible import/export.
·         Portable project - copy paste.
·         It can’t handle data with different projection system into a single project, have to have other method to get through
3D Visualization
Yes (Require Windows7 and high VRam), But No interpretation can be done, only view
Very fast interactive interpretation can be done
Mapping
·    GeoAtlas is excellent
·    Streaming Satellite Images(require URL)
·    Attaching Geo-referenced image
·         No separate mapping Tool; inbuilt in basemap
·         Webmap services over internet streaming geo-referenced images
Microseismic
Yes
Yes
Geomodeling
·   Geomodeling While drilling: Geosteering of horizontal wells
·      3D Geomodeling
·      Geosteering available
Field Planning
Yes, Map fracture Ellipse by padding required specially for shale Gas to optimize field production
Not yet available
Attribute
Curvature attribute and Coherency; 
·      50+ Attribute Geometrical, Wavelet, Instantaneous & Spectral Decomposition (Gabor-Morlet Method). 
·      Generates math based attribute 
Volumetrix & Reserve Estimation
·      Volumetrics can be done
·      Deterministic Analysis and no probabilistic, analysis can be done
·      Multi-layer Volumetrix
·      Probabilistic Analysis(Monte Carlo Simulation)
·      Economic Analysis
Inversion
No
Yest ; Works with gather also for
AVO; 1DFM; Colored Inversion
Synthetic
Yes
Yes; Advance GeoSyn module available
Petrophysical
Deterministic Interpretation
Available
Post Stack processing
Yes
Yes
Velocity Modeling
No
Yes
Direct Data Transfer 
·      SeisXchange;
·      WellXchange:
·      GridXchange to openworks
·      Open spirit; 
·      Landmark connect ; 
·      Kingdom Petrel Gateway available to transfer data to petrel;
·      Data can be directly transfered on the fly from mapping  software like Petrosys and pushed back to Petrel
A comparison based on industry based seismic interpretation software


Thursday, 17 March 2011

Inversion techniques

Seismic Data Inversion


I. Post-Stack Inversion


In the classic post-stack domain, analysis of 2D or 3D post-stack seismic volumes is carried on to produce an acoustic impedance volume.

Seismic inversion may be based on the following methods:
  1. Model-Based Inversion
  2. Sparse-Spike Inversion-
  3. Bandlimited - recursive-fast track Inversion
  4. Colored Inversion
  5. Neural Network Inversion
The data input to a post-stack STRATA project typically consist of a set of wells, with sonic and density logs, a series of interpreted horizons, and a seismic volume.



Current fast track methods for band-limited inversion to relative impedance are usually direct unconstrained transforms of the seismic data, such as phase rotation, trace integration and recursive inversion, and as such are prone to error because no account is taken of the seismic wavelet or calibration made to the earth.  More sophisticated techniques such as sparse-spike inversion do take account of these factors, but require specialist skills to implement correctly. 

Coloured Inversion implicitly takes into accounts for the seismic wavelet, is consistent with log data, yet is easy and fast to implement. Additionally it benchmarks well against unconstrained Sparse-Spike inversion. Coloured Inversion ( CI ) invert data within hours and establish a base-case against which all subsequently more sophisticated techniques must be judged.


II. Prestack Inversion

In the pre-stack domain, STRATA analyzes angle gathers or angle stacks to produce volumes of acoustic impedance, shear impedance and density.

This can be especially useful for analyzing data with AVO anomalies. The seismic input for this option is typically angle gathers of NMO corrected data. Tools such as Super Gather, Trim Statics, Noise Attenuation, and Angle Gather are used as the input data to enhance the AVO inversion process.

Tuesday, 15 March 2011

AWD

Introduction

Dynamite was the only source used until 1953 when the weight dropping Thumper technique was introduced. Thumper’s advanced technology what we call AWD Automatic weight drop seismic energy source. But  Dynamite has serious drawback with respect to AWD handling safety, high drilling cost, logistical problems (villages, buildings and pipelines). Whereas Vibroseis requires many number in serial and a big crew running all the time, which also has a logistical issue(roads has to be conditioned more than that of AWD), whereas AWD can be mounted on a farm tractor or a heavy Russian Truck with minor modifications and easily implemented within a short time. Its problems are lesser in field with a single crew running for all mechanical issue, can be fixed with no time.

Basis for choosing a seismic energy source
  1. Penetration of the required depth
  2. Bandwidth for the required resolution
  3. Signal to Noise ratio
  4. Environment
  5. Availability and cost


The AWDis a vehicle mounted, high performance gas-charged accelerated (impact) weight drop energy source designed for use on seismic exploration surveys. Applications include: shallow and deep refraction surveys, 2D and 3D seismic reflection surveys, Vertical Seismic Profiling (VSP), and downhole seismic or LVL surveys. Where applicable, the AWD is an environmentally friendly, economical, and efficient alternative to using dynamite and vibrator seismic energy sources. For many 2D and 3D seismic surveys that use explosives or vibrators as the primary energy source in only a portion of the survey area, the Digipulse can also be used as an “infill” source to obtain near offset data in cultural zones that may prohibit drilling and the use of explosives and large vibrators. Deeper multi-layer refraction and LVL, and VSP surveys can be conducted without the requirement for costly shot hole drilling and dynamite, or vibrator seismic sources.


AWD Operation is Estreme conditions and in a very difficult terrain
AWD Operation in Hilly terrain mounted on a Heavy Truck

Theory of Operation - The Physics of AWD


Unlike the explosive energy source, which is placed in shot holes below the ground, the AWD weight  drop is a surface “impact” type seismic energy source. The energy produced by the AWD (approximately 70 to 100 K-joules) is derived from a large hardened steel hammer mass (1180 lbs) that impacts a ground coupled base plate. The most important part is the coupling and the source synchronization with the instrument. The synchronization can be done through radio signal( minimum delay) or through wireline. A sensor attached to the base plate connected to a geophone sensing the earths vibration, is used for synchronizing with the instrument. Mechanically, the AWD uses a hydraulic system to lift the hammer mass to a “mass loaded” position. In the loaded position, a nitrogen gas charged cylinder and piston  assembly applies a downward force on the hammer mass. The approximate pressure applied to the hammer mass is a minimum 1,000 psi, depending on the model. Testings are carried on to adjust the Nitrogen pressure, which gives the best subsurface information. When released, the hammer mass is propelled at high velocity to impact the base plate. 

Conclusions
The concept of using weight drop, or accelerated weight drop technology, is not new to the seismic exploration industry. Several forms of weight drop systems have been in use for seismic exploration during the last 40 years. Currently, the AWD system is the only latest version of high performance accelerated weight drop systems in production. Depending on surface soil conditions, the source signature characteristics and energy levels of the AWD are usually sufficient for acquiring up to 4-6 second records (two-way travel time). Current users of the AWD report reaching target depths of up to 10,000 feet under ideal conditions. Typical seismic target objectives range from 2,000 to 15,000 feet.
AWD can deliver sufficient energy to reach shallow geologic target objectives, the AWD  1180 gas charged system is a more economical alternative in almost every aspect. A new AWD with a 3,000 lb mass mounted on an Minibuggy platform is now available. These larger Digipulse systems will deliver more than twice the downgoing energy of the Digipulse AWD units, and will likely reach the deeper seismic targets while still offering the same economical and production efficiency benefits.
The Digipulse AWD is also a more suitable source alternative in areas with sensitive cultural and environmental issues. In some cases, the AWD can be combined with  dynamite or vibrator operations to acquire only nearoffset data within areas that restrict or prohibit using energy sources that are hazardous or damaging to the environment. Moreover, the production rate in terms of completing a survey can be much higher, and the overall operating costs associated with the AWD is much less than with seismic surveys using dynamite or Vibroseis.













Discussion- Seismic Data Processing

Seismic Data

Typical Raw Seismic Record where headwaves-refractions, reflections, airwaves, Surface waves, Ground-rolls are identified to undertand the different types of noise with reflection
Seismic Record delineating reflection and noise

Refraction Headwaves- Blue, Noise(Ground Roll, P & S wave wave)-Red,  Reflection-dark Blue, Surface Wave -Black, Air



Seismic Data Processing 
Alteration of seismic data to suppress noise Enhance signal and migrate seismic events to the appropriate location in space.
Processing steps typically include analysis of velocities and frequencies, static corrections, deconvolution, normal moveout, dip moveout, stacking, and migration, which can be performed before or after stacking. Seismic processing facilitates better interpretation because subsurface structures and reflection geometries are more apparent. 

Brute Stack
Field Data primarily after being loading in the processing software, a single velocity function is applied to NMO in the gather and after CDP stacking, we get the brute stack.






Data Reduction:
  • Demultiplex
  • Gain Recovery
  • Editing: for bad are NTBC 
  • Migration: Summing (Vertical stack)
  • Correlation 
  • Gain Function
Geometry Correction:
  • CDP Sorting
  • Receiver and Shot (Datum) static
  • Uphole Static
  • NMO Correction
  • Residual Static
1. Create Geometry 
Create geometry for 2D/3D using OB Log informations. Shot Point, Station Number, X, Y, Shot Elevation, Shot Static, Offset and Skid in Shot spread sheet. Receiver Spread starting according to the shot and spread geometry(End-on/Slip Spread), Rec Elevation, Receiver static. 
Binning the data with bin interval, max and min bin size

2. Merge Geometry
Merge the geometry with the raw SegY

3. Editing
Editing for dead/bad traces


4. Static corrections
Often called statics, a bulk shift of a seismic trace in time during seismic processing. A common static correction is the weathering correction, which compensates for a layer of low seismic velocity material near the surface of the Earth. Other corrections compensate for differences in topography and differences in the elevations of sources and receivers. 

A major goal of reflection processing is to provide reflectivity images of the correct reflector geometry, which can be thwarted by the statics problem. The statics problem is defined to be static time shifts introduced into the traces by, e.g., near-surface velocity anomalies and/or topography, which distort the true geometry of deep reflectors.
An undulating topography can produce moveout delays in the CSG that can also be interpreted as undulations in the reflector, even if the reflector is flat. After an elevation statics correction (i.e., time shifts applied to traces) the data appear to have been collected on a flat datum plane. Static shifts can also be introduced by near-surface velocity anomalies which usually delay the traveltimes, resulting in reflections having a non-hyperbolic moveout curve.

Land seismic data are usually recorded over an irregular surface, and static correction has long been a problem. There are two kinds of statics, long-wavelength (large time shift) and short-wavelength (small time shift) statics. They are solved in two separate stages. The first stage is datum correction, used to correct long-wavelength statics, by which simple vertical time shifts are applied to common-shot gathers and common-receiver gathers to datum from an irregular earth surface to a planar surface. This planar surface is either above, below, or through all the shot and receiver elevations, depending on the requirements of a survey area. The details of this datuming procedure can be found in seismic data processing manuals. Next, the sorted common midpoint (CMP) gathers are normal moveout (NMO) corrected and stacked to generate the zero-offset section. Then statistical residual static corrections are applied to the zero-offset stack section to flatten a user-predefined main reflection event, thus correcting the short-wavelength statics. Throughout the years, even though many methods have been developed to deal with this problem, static correction still remains a difficult unsolved problem. 





Differential weathering correction
A type of static correction that compensates for delays in seismic reflection or refraction times from one point to another, such as among geophone groups in a survey. These delays can be induced by low-velocity layers such as the weathered layer near the Earth's surface. 

Elevation correction
Any compensating factor used to bring measurements to a common datum or reference plane. In gravity surveying, elevation corrections include the Bouguer and free-air corrections. Seismic data undergo a static correction to reduce the effects of topography and low-velocity zones near the Earth's surface. Well log headers include the elevation of the drilling rig's kelly bushing and, for onshore locations, the height of the location above sea level, so that well log depths can be corrected to sea level. 

Kelly bushing
An adapter that serves to connect the rotary table to the kelly. The kelly bushing has an inside diameter profile that matches that of the kelly, usually square or hexagonal. It is connected to the rotary table by four large steel pins that fit into mating holes in the rotary table. The rotary motion from the rotary table is transmitted to the bushing through the pins, and then to the kelly itself through the square or hexagonal flat surfaces between the kelly and the kelly bushing. The kelly then turns the entire drillstring because it is screwed into the top of the drillstring itself. Depth measurements are commonly referenced to the KB, such as 8327 ft KB, meaning 8327 feet below the kelly bushing. 

Photograph of kelly bushing
Kelly bushing. The kelly bushing connects the kelly to the rotary table. This hexagonal kelly turns the entire drillstring


Bouguer correction
The adjustment to a measurement of gravitational acceleration to account for elevation and the density of rock between the measurement station and a reference level. It can be expressed mathematically as the product of the density of the rock, the height relative to sea level or another reference, and a constant, in units of mGal:

Bouguer correction

Strictly interpreted, the Bouguer correction is added to the known value of gravity at the reference station to predict the value of gravity at the measurement level. The difference between the actual value and the predicted value is the gravity anomaly, which results from differences in density between the actual Earth and reference model anywhere below the measurement station. 

Free-air correction
In gravity surveying, a correction of 0.3086 mGal/m [0.09406 mGal/ft] added to a measurement to compensate for the change in the gravitational field with height above sea level, assuming there is only air between the measurement station and sea level. 

First break
The earliest arrival of energy propagated from the energy source at the surface to the geophone in the wellbore in vertical seismic profiles and check-shot surveys, or the first indication of seismic energy on a trace. On land, first breaks commonly represent the base of weathering and are useful in making static corrections.

TAR- True Amplitude Recovery
Steps in seismic processing to compensate for attenuation, spherical divergence and other effects by adjusting the amplitude of the data. The goal of TAR is to get the data to a state where the reflector amplitudes relate directly to the change in rock properties giving rise to them. 

II - Data Enhancement

  •         Mute
  •        Bandpass Filter
  •         Notch Filter
  •        Deconvolve
  •        2-D (F-K) Filter
  •       Stack
  •       Trace Equalize(AGC)

FK Filtering
FK analysis is done for delineating noise level from the data, os it depends on the stack, to decide when to usek FK filter. Events with a slow moveout velocity (1000-4000 ft/s) are typically unwanted noise such as surface waves. Their velocity is usually well separated from the apparent velocity of deep reflections so we can filter them out in a domain in which they are well separated from one another, namely the FK (i.e., frequency-wavenumber) domain. Thus muting one line from the other based on their different slopes (i.e., velocities) is trivial in the FK domain.

Stack after FK Analysis 


Deconvolution
The recorded seismic signal may be considered as the convolution of the source signal with the instruments, the geophones, and the response of the earth. The earth response includes some undesirable effects, such as reverberation, attenuation, and ghosting. The objective of deconvolution is to estimate these effects as linear filters, and then design and apply inverse filters. 
A step in seismic signal processing to recover high frequencies, attenuate multiples, equalize amplitudes, produce a zero-phase wavelet or for other purposes that generally affect the waveshape. Deconvolution, or inverse filtering, can improve seismic data that were adversely affected by filtering, or convolution that occurs naturally as seismic energy is filtered by the Earth. Deconvolution can also be performed on other types of data, such as gravity, magnetic or well log data. 

Decon with PD 4 and OL 200

Normal Moveout
The effect of the separation between receiver and source on the arrival time of a reflection that does not dip, abbreviated NMO. A reflection typically arrives first at the receiver nearest the source. The offset between the source and other receivers induces a delay in the arrival time of a reflection from a horizontal surface at depth. A plot of arrival times versus offset has a hyperbolic shape.

Diagram of common midpoint and moveout
The traces from different source-receiver pairs that share a midpoint, such as receiver 6 (R6), can be corrected during seismic processing to remove the effects of different source-receiver offsets, called normal moveout or NMO. After NMO corrections, the traces can be stacked to improve the signal-to-noise ratio.

Slant stack
A process used in seismic processing to stack, or sum, traces by shifting traces in time in proportion to their offset. This technique is useful in areas of dipping reflectors.
Slant stack is a transformation of the offset axis. It is like steering a beam of seismic waves. It is used as a part of a migration method. Mathematically, the slant-stack concept is found in the Radon [1917] transformation.


Linear moveout converts the task of identifying tangencies to constructed parallel lines to the task of locating the tops of convex events.
The slant-stack idea resembles the Snell trace method of organizing data around emergent angle. The Snell trace idea selects data based on a hypothetical velocity predicting the local stepout pdt/dx. Slant stack does not predict the stepout, but extracts it by filtering. Thus slant stack does its job correctly whether or not the velocity is known. When the velocity of the medium is known, slant stack enables immediate downward continuation even when mixed apparent velocities are present as with diffractions and multiple reflections.

Travel-time curves for a data gather on a multilayer earth model of constant velocity before and after slant stacking.
In other words, slant stacking takes us from two dimensions to one, but a $\sqrt{t}$ remains to correct the conical wavefront of three dimensions to the plane wave of two.

Stack
A processed seismic record that contains traces that have been added together from different records to reduce noise and improve overall data quality. The number of traces that have been added together during stacking is called the fold.

To sum traces to improve the signal-to-noise ratio, reduce noise and improve seismic data quality. Traces from different shot records with a common reflection point, such as common midpoint (CMP) data, are stacked to form a single trace during seismic processing. Stacking reduces the amount of data by a factor called the fold. 

Velocity
The rate at which a wave travels through a medium (a scalar) or the rate at which a body is displaced in a given direction (a vector), commonly symbolized by v. Unlike the physicist's definition of velocity as a vector, its usage in geophysics is as a property of a medium-distance divided by traveltime. Velocity can be determined from laboratory measurements, acoustic logs, vertical seismic profiles or from velocity analysis of seismic data. Velocity can vary vertically, laterally and azimuthally in anisotropic media such as rocks, and tends to increase with depth in the Earth because compaction reduces porosity. Velocity also varies as a function of how it is derived from the data. For example, the stacking velocity derived from normal moveout measurements of common depth point gathers differs from the average velocity measured vertically from a check-shot or vertical seismic profile (VSP). Velocity would be the same only in a constant velocity (homogeneous) medium. 


NMO Stretching problem(Red) and static problem(blue)


25 % Strech mute

Velocity Analysis  
1Km-V - 2D Supergather- BandPass - Vel Analysis precompute   - Vel Analysis.
First PassOnce flattened, the traces in a NMO-corrected gather can be stacked together for  constructive reinforcement of the reflection events.

Determining the stacking/NMO velocities from the CMG. Here, CMG amplitudes are summed together along different hyperbolas described by the traveltime equation 

t(x) =  
               

V2 =

The summation is carried out for different values of (t0, VNMO), and the result is contoured in (t0,VNMO) space (see the RHS of Figure 1.17). The correct values of (t0, VNMO) will describe an hyperbola that coincides with the reflection event to give a large summation value; otherwise the summation values (or semblance) are small. The summation value is called a semblance value because we really sum over a thick hyperbolic line, weight the summation with a normalization factor, and insure positive summation values by taking the absolute value of the sum. Identifying the bullseyes on the RHS of this figure, and connecting lines between these bullseyes describes the optimal vNMO vs t0curve. The vNMO can also be thought of as the stacking velocity for a layered medium, and is considered to be good estimate of vRMS.

Apparent velocity
The speed of a wavefront in a certain direction, typically measured along a line of receivers and symbolized by va. Apparent velocity and velocity are related by the cosine of the angle at which the wavefront approaches the receivers:



Apparent velocity


Velocity correction
A change made in seismic data to present reflectors realistically. Velocity corrections typically require that assumptions be made about the seismic velocities of the rocks or sediments through which seismic waves pass. 


Stacking velocity
The distance-time relationship determined from analysis of normal moveout (NMO) measurements from common depth point gathers of seismic data. The stacking velocity is used to correct the arrivaltimes of events in the traces for their varying offsets prior to summing, or stacking, the traces to improve the signal-to-noise ratio of the data. 

Interval velocity
The velocity, typically P-wave velocity, of a specific layer or layers of rock, symbolized by vint and commonly calculated from acoustic logs or from the change in stacking velocity between seismic events on a common midpoint gather. 


Dix
An equation used to calculate the interval velocity within a series of flat, parallel layers, named for American geophysicist C. Hewitt Dix (1905 to 1984). Sheriff (1991) cautions that the equation is misused in situations that do not match Dix's assumptions. The equation is as follows:
Dix formula

Average velocity
In geophysics, the depth divided by the traveltime of a wave to that depth. Average velocity is commonly calculated by assuming a vertical path, parallel layers and straight raypaths, conditions that are quite idealized compared to those actually found in the Earth. 


Wavelet extraction
A step in seismic processing to determine the shape of the wavelet, also known as the embedded wavelet, that would be produced by a wave train impinging upon an interface with a positive reflectioncoefficient. Wavelets may also be extracted by using a model for the reflections in a seismic trace, such as a synthetic seismogram. A wavelet is generated by deconvolving the trace with the set of reflection coefficients of the synthetic seismogram, a process also known as deterministic wavelet extraction. Wavelets may be extracted without a model for the reflections by generating a power spectrum of the data. By making certain assumptions, such as that the power spectrum contains information about the wavelet (and not the geology) and that the wavelet is of a certain phase (minimum, zero), a wavelet may be generated. This is also called statistical wavelet extraction. A particular processing approach to establishing the embedded wavelet is to compare the processed seismic response with the response measured by a vertical seismic profile (VSP) or generated synthetically through a synthetic seismogram in which the embedded wavelet is known. The wavelet can also be extracted through the autocorrelationof the seismic trace, in which case the phase of the wavelet has to be assumed. 

Residual Statics :
Static shifts introduced by topographic variations fall under the class of field statics, and those due to near-surface lithological variations that occur within a cable length fall under the class of residual statics. Correcting for static shifts in the traces can make a significant difference in the quality of a migrated or stacked image. It is easy to determine elevation static corrections, but not so easy to find the residual static corrections. One means is to determine the near-surface velocity distribution by refraction tomography.

It is applied calculating residual shift for consecutive traces in a gather. For a time gate, and respective picked time shift,  it  creates a delta-T matrix and then solves the matrix to give the final output.

DMO
A seismic processing operation to correct for the fact that, for dipping reflections, the traces of a CMP gather do not involve a common reflection point. DMO effectively corrects for the reflection-point smear that results when reflectors dip. 

Events with various dips stack with the same velocity after DMO.









III IMAGING 
Post-Stack Migrate  MIGRAITON AFTER STACK( After twice Residual static DMO is applied twice to get to raw stack and then DMO pics applied to the raw stack to get final post stack)


Depth Conversion  - VELOCTY MODEL SENSETIVE  - VERY DELICATE TO HANDLE
Pre-Stack Migrate  - EXPENSIVE - MIGRATION BEFORE STACK- VELOCITY ANALYSIS OF GATHERS - Migration is done after Residual Static correction on gathers.
  
Migration
The movement of hydrocarbons from their source into reservoir rocks. The movement of newly generated hydrocarbons out of their source rock is primary migration, also called expulsion. The further movement of the hydrocarbons into reservoir rock in a hydrocarbon trap or other area of accumulation is secondary migration. Migration typically occurs from a structurally low area to a higher area because of the relative buoyancy of hydrocarbons in comparison to the surrounding rock. Migration can be local or can occur along distances of hundreds of kilometers in large sedimentary basins, and is critical to the formation of a viable petroleum system. 

A step in seismic processing in which reflections in seismic data are moved to their correct locations in the x-y-time space of seismic data, including two-way traveltime and position relative to shotpoints. Migration improves seismic interpretation and mapping because the locations of geological structures, especially faults, are more accurate in migrated seismic data. Proper migration collapses diffractions from secondary sources such as reflector terminations against faults and corrects bow ties to form synclines. There are numerous methods of migration, such as dip moveout (DMO), frequency domain, ray-trace and wave-equation migration. 


Diagram of a datum and hydrophone array
During seismic processing, migration adjusts the location of events in seismic traces to compensate for dipping reflectors.

Time migration
A migration technique for processing seismic data in areas where lateral velocity changes are not too severe, but structures are complex. Time migration has the effect of moving dipping events on a surface seismic line from apparent locations to their true locations in time. The resulting image is shown in terms of traveltime rather than depth, and must then be converted to depth with an accurate velocity model to be compared to well logs.

Kirchhoff migration
A method of seismic migration that uses the integral form (Kirchhoff equation) of the wave equation. All methods of seismic migration involve the backpropagation (or continuation) of the seismic wavefield from the region where it was measured (Earth's surface or along a borehole) into the region to be imaged. In Kirchhoff migration, this is done by using the Kirchhoff integral representation of a field at a given point as a (weighted) superposition of waves propagating from adjacent points and times. Continuation of the wavefield requires a background model of  seismic velocity, which is usually a model of constant or smoothly varying velocity. Because of the integral form of Kirchhoff migration, its implementation reduces to stacking the data along curves that trace the arrival time of energy scattered by image points in the earth. 

smile
A concave-upward, semicircular event in seismic data that has the appearance of a smile and can be caused by poor data migration or migration of noise.

PSTM:
Prestack Kirchoff Time Migration Dips upto 90 deg. This Algorithm  is applied on common offset non- NMO Corrected gather(through trace binning) , uses both vertically and lateral RMS/stacking velocity from floating datum to final datum(Flat datum).
Maximum Frequency: Less than Nyquist frequency will speed up the computation and above it should be band filtered.
Migration Aperture: Horizontal resolution is twice the vertical resolution so depending upon the depth of interest the apperture test is done to take care of the width of the horizontal distance that energy can migrate.Tapering should be applied. Migration can be split into four conceptual pieces. As a rule of thumb, these four pieces will help us understand what migration is and how it naturally completes the imaging process.
1.The first piece is called normal moveout(NMO). When the world is flat, NMO corrects for the fact that the source and receiver are not coincident, but it cannot do so when the reflections come from dipping horizons.
2.The second piece of migration corrects for dip. Historically, this second piece was called dip-moveout(DMO), but, in the cases of interest here, it happens within the migration methodology itself.
3.The third migration piece shifts events on each moveout-corrected offset to its true subsurface position.
4.The fourth and final piece sums (stacks) all the redundant traces into the final image.
time Migration doesn’t accounts for ray bending, but  interval velocity in depth migration takes care of that.
  • Kirchoff
  • FK
  • Phase Shift
  • Finite Difference
  • Reverse Time


FX: implemented with spatially variant convolutional filters and are often grouped with finite difference technique12. DAS -Q Compensation - Spectral Whitening - Filter test RADON transforms data to time-moveout space for picking a mute to attenuate energy with undesirable moveout. Radon Filter is commonly used for suppression of multiples. The normal technique is to model multiples and subtract them from the input seismic trace data. Radon filters are then applied to pass the primary energy. However, in practice this tends to produce an artificial and wormy appearing result.