Difference between revisions of "Procesamiento sísmico"

From Sinfronteras
Jump to: navigation, search
(Secciones Sísmicas)
 
(50 intermediate revisions by the same user not shown)
Line 1: Line 1:
 +
[[File:Aplicaciones_Geofisicas_de_los_filtros_digitales-Mario_Caicedo-Milagrosa_Aldana.pdf]] -- [https://docplayer.es/16950716-Aplicaciones-geofisicas-de-los-filtros-digitales.html It's been downloaded from this link] [Mario Caicedo & Milagrosa Aldana]
 +
 +
 +
This is a good explanation about the Relations and differences between time-series analysis and statistical signal processing? https://stats.stackexchange.com/questions/52270/relations-and-differences-between-time-series-analysis-and-statistical-signal-pr#:~:text=2%20Answers&text=As%20a%20signal%20is%20by,significant%20overlap%20between%20the%20two.
 +
 +
 +
<br />
 
==Convetional Seismic data processing==
 
==Convetional Seismic data processing==
  
===Yilmaz (2001)===
+
 
 +
<br />
 +
===Yilmaz - 2001===
 
{| class="wikitable"
 
{| class="wikitable"
|- style="margin: auto; width: 100%; vertical-align:center;"
+
|- style="vertical-align: top;  background: white;"
| rowspan="6" | '''1 Pre-processing'''
+
| rowspan="6" style="width: 100pt; text-align: center" |  
| '''Demultiplexing'''
+
'''1 Pre-processing'''
 +
| style="width: 150pt" |'''Demultiplexing'''
 
| Orders seismic data by seismic trace instead of being ordered by time. Multiplexed recording is still found in older data or in current data shots with old recording instruments.
 
| Orders seismic data by seismic trace instead of being ordered by time. Multiplexed recording is still found in older data or in current data shots with old recording instruments.
 
| '''Démultiplexage:'''
 
| '''Démultiplexage:'''
|-
+
|- style="vertical-align: top;  background: white;"
 
| '''Reformatting'''
 
| '''Reformatting'''
 
| '''(SEG-Y/SEG-D)'''
 
| '''(SEG-Y/SEG-D)'''
 
Puts the data in a format that is understandable by computers (e.g. SEG-Y format, SEG-D format, etc.).
 
Puts the data in a format that is understandable by computers (e.g. SEG-Y format, SEG-D format, etc.).
 
| '''Reformatage:'''
 
| '''Reformatage:'''
|-
+
|- style="vertical-align: top;  background: white;"
 
| '''Seismic data edition'''
 
| '''Seismic data edition'''
 
| Searchin for noisy, monofrequency and incorrect polarities traces.
 
| Searchin for noisy, monofrequency and incorrect polarities traces.
 
Examines seismic data for bad, noisy, and\or monofrequency traces, and deletes them, and also spots any traces with incorrect polarities and corrects them.
 
Examines seismic data for bad, noisy, and\or monofrequency traces, and deletes them, and also spots any traces with incorrect polarities and corrects them.
 
| '''Édition:'''
 
| '''Édition:'''
|-
+
|- style="vertical-align: top;  background: white;"
 
| '''Geometrical spreading correction'''
 
| '''Geometrical spreading correction'''
 
| It is also called amplitude correction. Because the signal becomes weaker as it moves away from its source (i.e. its amplitude decreases with time), amplitude correction is needed to count for this lose, and to make amplitudes stronger.
 
| It is also called amplitude correction. Because the signal becomes weaker as it moves away from its source (i.e. its amplitude decreases with time), amplitude correction is needed to count for this lose, and to make amplitudes stronger.
 
| '''Correction de la dispersion géométrique:'''
 
| '''Correction de la dispersion géométrique:'''
|-
+
|- style="vertical-align: top;  background: white;"
 
| '''Set-up of field geometry - Geometry QC'''
 
| '''Set-up of field geometry - Geometry QC'''
 
| Incorporates field geometry with seismic data processing. Coordinates, shot\ receiver locations and spacing must be entered to the system carefully and precisely because processes like CMP sorting, for example, highly depend on this.
 
| Incorporates field geometry with seismic data processing. Coordinates, shot\ receiver locations and spacing must be entered to the system carefully and precisely because processes like CMP sorting, for example, highly depend on this.
 
http://radexpro.com/scope/qc/
 
http://radexpro.com/scope/qc/
 
| '''Géométrie de terrain:'''
 
| '''Géométrie de terrain:'''
|-
+
|- style="vertical-align: top;  background: white;"
 
| '''Application of field statics corrections'''
 
| '''Application of field statics corrections'''
 
| This is needed for land data collected on non-flat areas. It reduces the travel times into a common datum level. The datum could be the sea level or any other local datum.
 
| This is needed for land data collected on non-flat areas. It reduces the travel times into a common datum level. The datum could be the sea level or any other local datum.
 
| '''Corrections statiques:'''
 
| '''Corrections statiques:'''
|-
+
|- style="vertical-align: top;  background: white;"
| '''2'''
+
| style="text-align: center" | '''2'''
 
| '''Deconvolution and trace balancing'''
 
| '''Deconvolution and trace balancing'''
| Deconvolution makes the signal look better by increasing the temporal resolution and removing echoes. Trace balancing makes the amplitude uniform.
+
|
 +
<div class="mw-collapsible mw-collapsed">
 +
'''Deconvolution''' makes the signal look better by increasing the temporal resolution and removing echoes:
 +
<div class="mw-collapsible-content">
 +
https://www.youtube.com/watch?v=kzNlXZ-8tTs
 +
 
 +
The seismic wave can be mathematically represented by what is called the '''Convolutional model for the seismic trace'''. This is a mathematical model that states that the Seismic trace is the result of the convolution between the Seismic Wavelet (The initial wave produced by the seismic source <math>W(t)</math>) and the Reflectivity function.
 +
 
 +
The reverse calculation, which is the process of calculating the Reflectivity function from the trace is known as "Deconvolution". So, when we apply Deconvolution to real seismic data, we are looking to extract the reflectivity function from the seismic trace, and the reflectivity function should be a better representation of the reflectors in the subsurface. It should highlight structures better and increase the temporal resolution.
 +
 
 +
[[File:The_convolutional_model_for_the_seismic_trace1.png|700px|thumb|center|Taken from https://www.youtube.com/watch?v=kzNlXZ-8tTs]]
 +
</div>
 +
</div>
 +
 
 +
 
 +
'''Trace balancing''' makes the amplitude uniform.
 
| '''Déconvolution et normalisation des traces:'''
 
| '''Déconvolution et normalisation des traces:'''
|-
+
|- style="vertical-align: top;  background: white;"
| '''3'''
+
| style="text-align: center" | '''3'''
 
| '''CMP sorting'''
 
| '''CMP sorting'''
 
| Sorts the data into CMP gathers, so it can be corrected for the NMO and stacked after that (see step 7).
 
| Sorts the data into CMP gathers, so it can be corrected for the NMO and stacked after that (see step 7).
 
| '''Triage par CMP:'''
 
| '''Triage par CMP:'''
|-
+
|- style="vertical-align: top;  background: white;"
| '''4'''
+
| style="text-align: center" | '''4'''
 
| '''Velocity analysis'''
 
| '''Velocity analysis'''
 
| Gives info about velocities in the subsurface layers. It finds the stacking velocity (very close in value to the RMS velocity) that best fits our data (This step may be delayed after step 5).
 
| Gives info about velocities in the subsurface layers. It finds the stacking velocity (very close in value to the RMS velocity) that best fits our data (This step may be delayed after step 5).
 
| '''Analyse de vitesse:'''
 
| '''Analyse de vitesse:'''
|-
+
|- style="vertical-align: top;  background: white;"
| '''5'''
+
| style="text-align: center" | '''5'''
 
| '''Residual statics corrections'''
 
| '''Residual statics corrections'''
 
| Counts for near-surface velocity variations that causes some static and dynamic distortion problems.
 
| Counts for near-surface velocity variations that causes some static and dynamic distortion problems.
 
| '''Corrections statiques résiduelles:'''
 
| '''Corrections statiques résiduelles:'''
|-
+
|- style="vertical-align: top;  background: white;"
| '''6'''
+
| style="text-align: center" | '''6'''
 
| '''Velocity analysis'''
 
| '''Velocity analysis'''
 
| (see step 4). If we have residual statics problems, we do velocity analysis after we count for them.
 
| (see step 4). If we have residual statics problems, we do velocity analysis after we count for them.
 
| '''Analyse de vitesse:'''
 
| '''Analyse de vitesse:'''
|-
+
|- style="vertical-align: top;  background: white;"
| '''7'''
+
| style="text-align: center" | '''7'''
 
| '''NMO Correction'''
 
| '''NMO Correction'''
 
| Counts for the increase in travel time with increasing offset distance. This increase makes flat reflectors look dipping, and makes dipping reflectors look even more dipping. The amount of correction needed decreases with depth, so that shallower reflectors get more “stretched” than do the deeper ones.
 
| Counts for the increase in travel time with increasing offset distance. This increase makes flat reflectors look dipping, and makes dipping reflectors look even more dipping. The amount of correction needed decreases with depth, so that shallower reflectors get more “stretched” than do the deeper ones.
 
Muting: It is just a fancy word for deleting a part of a trace. Muting a whole trace is called “killing”.  
 
Muting: It is just a fancy word for deleting a part of a trace. Muting a whole trace is called “killing”.  
 
| '''Correction NMO:'''
 
| '''Correction NMO:'''
|-
+
|- style="vertical-align: top;  background: white;"
 
|
 
|
 
|  
 
|  
 
|
 
|
 
| '''Correction DMO:'''
 
| '''Correction DMO:'''
|-
+
|- style="vertical-align: top;  background: white;"
 
|
 
|
 
|
 
|
 
|
 
|
 
| '''Correction NMO inverse:'''
 
| '''Correction NMO inverse:'''
|-
+
|- style="vertical-align: top;  background: white;"
 
|
 
|
 
|
 
|
 
|
 
|
 
| '''Analyse de vitesse:'''
 
| '''Analyse de vitesse:'''
|-
+
|- style="vertical-align: top;  background: white;"
 
|  
 
|  
 
|  
 
|  
 
|  
 
|  
 
| '''Correction NMO:'''
 
| '''Correction NMO:'''
|-
+
|- style="vertical-align: top;  background: white;"
| '''8'''
+
| style="text-align: center" | '''8'''
 
| '''Stacking'''
 
| '''Stacking'''
 
| After sorting the data into CMP gathers and applying the NMO correction, reflectors line up nicely and hence their stacking gives a stronger signal. Multiples\ random noises do not line up. Stacking increases S/N ratio by decreasing multiples\ random noise from the data which enhances the overall quality. It also reduces the seismic data volume to the plane of the seismic section
 
| After sorting the data into CMP gathers and applying the NMO correction, reflectors line up nicely and hence their stacking gives a stronger signal. Multiples\ random noises do not line up. Stacking increases S/N ratio by decreasing multiples\ random noise from the data which enhances the overall quality. It also reduces the seismic data volume to the plane of the seismic section
 
| '''Sommation:'''
 
| '''Sommation:'''
|-
+
|- style="vertical-align: top;  background: white;"
| '''9'''
+
| style="text-align: center" | '''9'''
 
|  
 
|  
 
|
 
|
 
| '''Déconvolution:'''
 
| '''Déconvolution:'''
|-
+
|- style="vertical-align: top;  background: white;"
| '''10'''
+
| style="text-align: center" | '''10'''
 
|  
 
|  
 
|
 
|
 
| '''Blanchiment de spectre à temps variable:'''
 
| '''Blanchiment de spectre à temps variable:'''
|-
+
|- style="vertical-align: top;  background: white;"
| '''11'''
+
| style="text-align: center" | '''11'''
| '''Time- variant band-pass filtering'''
+
| '''Time-variant band-pass filtering'''
| Filters unwanted “signals” based on their frequencies. For example, ground roll has a lower frequency (and higher amplitude) compared to the rest of the section, so we can filter it out based on that fact.
+
| Filters unwanted "signals" based on their frequencies. For example, ground roll has a lower frequency (and higher amplitude) compared to the rest of the section, so we can filter it out based on that fact.  
 +
 
 +
 
 +
[noise reduction] [time/frequency domain transforms]
 +
 
 +
 
 +
https://www.glossary.oilfield.slb.com/en/Terms/g/ground_roll.aspx
 +
 
 +
https://wiki.seg.org/wiki/Ground_roll
 
| '''Filtre à temps variable:'''
 
| '''Filtre à temps variable:'''
|-
+
|- style="vertical-align: top;  background: white;"
| '''12'''
+
| style="text-align: center" | '''12'''
 
| '''Migration'''
 
| '''Migration'''
 
| So far, each trace is plotted under its CMP location which puts reflectors in the wrong subsurface location. Migration process moves those reflectors into their true subsurface locations, and that improves lateral resolution. Migration also collapses diffractions into identifiable points on the seismic section.
 
| So far, each trace is plotted under its CMP location which puts reflectors in the wrong subsurface location. Migration process moves those reflectors into their true subsurface locations, and that improves lateral resolution. Migration also collapses diffractions into identifiable points on the seismic section.
 
| '''Migration:'''
 
| '''Migration:'''
|-
+
|- style="vertical-align: top;  background: white;"
| '''13'''
+
| style="text-align: center" | '''13'''
 
| '''Gain recovery'''
 
| '''Gain recovery'''
 
| Seismic energy gets lost in many different ways (e.g. scattering, frication, etc.). Gain is the inverse function of energy loss. It is very hard to predict the attenuation function because primary reflections, multiple reflections, and random noise have different decay functions. Instead, we examine (test) different gain values, and see which value gives a better looking data.  
 
| Seismic energy gets lost in many different ways (e.g. scattering, frication, etc.). Gain is the inverse function of energy loss. It is very hard to predict the attenuation function because primary reflections, multiple reflections, and random noise have different decay functions. Instead, we examine (test) different gain values, and see which value gives a better looking data.  
Line 114: Line 147:
 
|}
 
|}
  
==what is different between VSP and check- shot survey? please explain me...==
+
 
 +
<br />
 +
==What is different between VSP and check-shot survey==
 
https://www.linkedin.com/grp/post/1781702-189955245
 
https://www.linkedin.com/grp/post/1781702-189955245
  
==Secciones Sísmicas==
+
 
===Migración pre-apilamiento en tiempo (Kirchhoff)===
+
<br />
 +
==Seismic sections==
 +
 
 +
 
 +
<br />
 +
===Migracion pre-apilamiento en tiempo - Kirchhoff===
 
[[File:Inline_1450_PreSTM.jpg |900px | thumb | center |Inline 1450 (Kirchhoff PreSTM).]]
 
[[File:Inline_1450_PreSTM.jpg |900px | thumb | center |Inline 1450 (Kirchhoff PreSTM).]]
  
 
[[File:Inline_1450_PreSTM-acercamiento.jpg |900px | thumb | center |Acercamiento sobre la inline 1450 (Kirchhoff PreSTM).]]
 
[[File:Inline_1450_PreSTM-acercamiento.jpg |900px | thumb | center |Acercamiento sobre la inline 1450 (Kirchhoff PreSTM).]]
  
===Migración pre-apilamiento en profundidad (Kirchhoff)===
+
 
 +
<br />
 +
===Migracion pre-apilamiento en profundidad - Kirchhoff===
 
[[File:Inline_1450_PreSDM-acercamiento.jpg |900px | thumb | center |PreSDM empleando distintos campos de velocidad (inline 1450). Las dos primeras, de izquierda a derecha, se corresponden con el modelo generado en 2008 y con el actual. La tercera ha sido migrada empleando un modelo similar al actual pero disminuyendo la velocidad de la capa de baja velocidad hasta 2100 ''m/s'', y la cuarta es el mismo caso anterior pero con una velocidad de 1900 ''m/s''.]]
 
[[File:Inline_1450_PreSDM-acercamiento.jpg |900px | thumb | center |PreSDM empleando distintos campos de velocidad (inline 1450). Las dos primeras, de izquierda a derecha, se corresponden con el modelo generado en 2008 y con el actual. La tercera ha sido migrada empleando un modelo similar al actual pero disminuyendo la velocidad de la capa de baja velocidad hasta 2100 ''m/s'', y la cuarta es el mismo caso anterior pero con una velocidad de 1900 ''m/s''.]]
 +
 +
 +
<br />

Latest revision as of 19:54, 23 February 2021

File:Aplicaciones Geofisicas de los filtros digitales-Mario Caicedo-Milagrosa Aldana.pdf -- It's been downloaded from this link [Mario Caicedo & Milagrosa Aldana]


This is a good explanation about the Relations and differences between time-series analysis and statistical signal processing? https://stats.stackexchange.com/questions/52270/relations-and-differences-between-time-series-analysis-and-statistical-signal-pr#:~:text=2%20Answers&text=As%20a%20signal%20is%20by,significant%20overlap%20between%20the%20two.



Convetional Seismic data processing


Yilmaz - 2001

1 Pre-processing

Demultiplexing Orders seismic data by seismic trace instead of being ordered by time. Multiplexed recording is still found in older data or in current data shots with old recording instruments. Démultiplexage:
Reformatting (SEG-Y/SEG-D)

Puts the data in a format that is understandable by computers (e.g. SEG-Y format, SEG-D format, etc.).

Reformatage:
Seismic data edition Searchin for noisy, monofrequency and incorrect polarities traces.

Examines seismic data for bad, noisy, and\or monofrequency traces, and deletes them, and also spots any traces with incorrect polarities and corrects them.

Édition:
Geometrical spreading correction It is also called amplitude correction. Because the signal becomes weaker as it moves away from its source (i.e. its amplitude decreases with time), amplitude correction is needed to count for this lose, and to make amplitudes stronger. Correction de la dispersion géométrique:
Set-up of field geometry - Geometry QC Incorporates field geometry with seismic data processing. Coordinates, shot\ receiver locations and spacing must be entered to the system carefully and precisely because processes like CMP sorting, for example, highly depend on this.

http://radexpro.com/scope/qc/

Géométrie de terrain:
Application of field statics corrections This is needed for land data collected on non-flat areas. It reduces the travel times into a common datum level. The datum could be the sea level or any other local datum. Corrections statiques:
2 Deconvolution and trace balancing

Deconvolution makes the signal look better by increasing the temporal resolution and removing echoes:

https://www.youtube.com/watch?v=kzNlXZ-8tTs

The seismic wave can be mathematically represented by what is called the Convolutional model for the seismic trace. This is a mathematical model that states that the Seismic trace is the result of the convolution between the Seismic Wavelet (The initial wave produced by the seismic source ) and the Reflectivity function.

The reverse calculation, which is the process of calculating the Reflectivity function from the trace is known as "Deconvolution". So, when we apply Deconvolution to real seismic data, we are looking to extract the reflectivity function from the seismic trace, and the reflectivity function should be a better representation of the reflectors in the subsurface. It should highlight structures better and increase the temporal resolution.


Trace balancing makes the amplitude uniform.

Déconvolution et normalisation des traces:
3 CMP sorting Sorts the data into CMP gathers, so it can be corrected for the NMO and stacked after that (see step 7). Triage par CMP:
4 Velocity analysis Gives info about velocities in the subsurface layers. It finds the stacking velocity (very close in value to the RMS velocity) that best fits our data (This step may be delayed after step 5). Analyse de vitesse:
5 Residual statics corrections Counts for near-surface velocity variations that causes some static and dynamic distortion problems. Corrections statiques résiduelles:
6 Velocity analysis (see step 4). If we have residual statics problems, we do velocity analysis after we count for them. Analyse de vitesse:
7 NMO Correction Counts for the increase in travel time with increasing offset distance. This increase makes flat reflectors look dipping, and makes dipping reflectors look even more dipping. The amount of correction needed decreases with depth, so that shallower reflectors get more “stretched” than do the deeper ones.

Muting: It is just a fancy word for deleting a part of a trace. Muting a whole trace is called “killing”.

Correction NMO:
Correction DMO:
Correction NMO inverse:
Analyse de vitesse:
Correction NMO:
8 Stacking After sorting the data into CMP gathers and applying the NMO correction, reflectors line up nicely and hence their stacking gives a stronger signal. Multiples\ random noises do not line up. Stacking increases S/N ratio by decreasing multiples\ random noise from the data which enhances the overall quality. It also reduces the seismic data volume to the plane of the seismic section Sommation:
9 Déconvolution:
10 Blanchiment de spectre à temps variable:
11 Time-variant band-pass filtering Filters unwanted "signals" based on their frequencies. For example, ground roll has a lower frequency (and higher amplitude) compared to the rest of the section, so we can filter it out based on that fact.


[noise reduction] [time/frequency domain transforms]


https://www.glossary.oilfield.slb.com/en/Terms/g/ground_roll.aspx

https://wiki.seg.org/wiki/Ground_roll

Filtre à temps variable:
12 Migration So far, each trace is plotted under its CMP location which puts reflectors in the wrong subsurface location. Migration process moves those reflectors into their true subsurface locations, and that improves lateral resolution. Migration also collapses diffractions into identifiable points on the seismic section. Migration:
13 Gain recovery Seismic energy gets lost in many different ways (e.g. scattering, frication, etc.). Gain is the inverse function of energy loss. It is very hard to predict the attenuation function because primary reflections, multiple reflections, and random noise have different decay functions. Instead, we examine (test) different gain values, and see which value gives a better looking data. Gain:



What is different between VSP and check-shot survey

https://www.linkedin.com/grp/post/1781702-189955245



Seismic sections


Migracion pre-apilamiento en tiempo - Kirchhoff

Inline 1450 (Kirchhoff PreSTM).
Acercamiento sobre la inline 1450 (Kirchhoff PreSTM).



Migracion pre-apilamiento en profundidad - Kirchhoff

PreSDM empleando distintos campos de velocidad (inline 1450). Las dos primeras, de izquierda a derecha, se corresponden con el modelo generado en 2008 y con el actual. La tercera ha sido migrada empleando un modelo similar al actual pero disminuyendo la velocidad de la capa de baja velocidad hasta 2100 m/s, y la cuarta es el mismo caso anterior pero con una velocidad de 1900 m/s.