en>Wcherowi: /* References */ alphabetize

2013-10-06T19:14:08Z

References: alphabetize

New page

'''''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 is possible to greatly increase the resolution of the data despite the presence of background noise. This new data is usually easier to interpret when trying to deduce the underground structure of an area.

== 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 however allows for the removal of the extra noise and leaves only the primary reflection.

== Process ==

[[File:KR-example.png|thumb|Uncorrected data showing hyperbolic curve.]]
[[File:Corrected wigb.png|thumb|Corrected data.]]

Semblance analysis allows for the refinement of seismic data. This is done by developing a velocity spectra display to determine the velocity through different layers at depth. The easiest way to accomplish this is by recording the normal incidence path (NIP). The NIP is the where you have the shot and the geophone in the same location and the path taken by the recorded sound waves in perpendicular to the boundaries between the layers. 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 the using the equation for the root mean square velocity starting with the top layer and working downward.

:<math>V_\mathrm{rms} = \sqrt { \frac { \sum t_i \cdot {V_i}^2 }{\sum t_i} }</math>

Once all of the velocities for the layers are known then it is possible to calculate the amount of time that is needed for the wave to travel the distance down to the midpoint between the each geophone and the shot point for each of the layers. As the geophones are farther away from the shot the time taken for the wave to travel there increases which 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 which will sum together all of the events that the traces share and remove the ones that they do not. The result is a single data set that has all of the primary peaks strongly displayed with most of the noise removed. <ref>{{cite journal|last=Taner|first=Turhan|coauthors=Fulton Koehler|title=Velocity Spectra-Digital Computer Derivation and Application of Velocity Functions|journal=Geophysics|year=1969|month=December|volume=34|issue=6|pages=859–881}}</ref>

==Problems==

While this technique can be very useful in the analysis there are several situations in which the process 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. Also in situations where there is large offset there can also be polarity reverses 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. <ref>{{cite journal|last=Yan|first=Jia|coauthors=Ilya Tsvankin|title=AVO-sensitive semblance analysis for wide-azimuth data|journal=Geophysics|year=2008|month=March-April|volume=73|issue=2|pages=U1–U11}}</ref>

==Code of program==
This code is for a program in perl that allows for the determination of multiple moveout velocities for the correction of seismic data.

#!/bin/sh

set -x

# nmo_test.sh

# Nov. 8, 2012

# Program to test nmo's

# several constant velocity moveouts are tested

# starting at 600 m/s and ending at 1000 m/s

# STEP 1: Data is sorted by cdp and offset

#STEP 1A: Data is windowed

# STEP 2: DATA is moved out

# STEP 3: data is filtered

# STEP 4: data is gained

# STEP 5: data is displayed

# Juan M. Lorenzo

# set up working directories

SU_DIR='Directory'

this_file='inpute_file'

counter=0

vel_start=7500

vel_last=35000

vel_inc=5000

first_cmp=1

last_cmp=1

for ((vel=$vel_start; vel<=$vel_last; vel=$vel+$vel_inc))

do

echo $vel

susort <$SU_DIR/$this_file.su cdp offset \

|

suwind key=cdp min=$first_cmp max=$last_cmp \

| \

sunmo vnmo=$vel \

| \

sufilter f=0,3,400,600 \

| \

sugain agc=1 wagc=0.1 \

| \

suximage \

\

xbox=$[$counter*200] ybox=0 wbox=200 hbox=600 \

title="$vel m/s" & \

counter=$[$counter+1]

done

==References==
{{Reflist}}

[[Category:Seismology]]

Linearised polynomial - Revision history

en>Wcherowi: /* References */ alphabetize