Procesado de Datos Sísmicos PDF

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 125

PROCESADO DE

DATOS SSMICOS

PROYECTO BRIGADA ESCUELA


DE GEOCIENCIAS
PEMEX UANL
BRIGADA 4
Elabor: Ing. Luis Mario Palacios Reyes

Revis: Ing. Luis Mario Palacios Reyes

Aprob: Hctor de Len Gmez

Fecha: 20.03.2006

Versin No.1

Pginas:125

Edicin No. 1

Procesado de datos ssmicos

TEMARIO
1. Fundamentos de Sistemas Lineales
2. Ruidos Ssmicos
3. Amplitudes Ssmicas
4. Filtros digitales
5. Auto-correlacin/ correlacin cruzada
6. Filtro optimo de Wiener
7. Modelo Convolucional
8. Correcciones Estticas
9. Correccin Dinmica
10. Anlisis de velocidad
11. Apilamiento de trazas
12. Velocidades Ssmicas
13. Factores que afectan las velocidades
14. Concepto de Migracin Ssmica
15. Seccin en tiempo
16. Seccin en profundidad

Procesado de datos ssmicos

PROCESADO DE DATOS SSMICOS


Objetivo: Uso eficiente de parmetros de proceso, para
mejorar la informacin ssmica adquirida en el campo y
proporcionar imgenes de calidad para su interpretacin.

Procesado de datos ssmicos

1.- FUNDAMENTOS DE SISTEMAS LINEALES

La comunicacin es el
intercambio de informacin

distancias muy cortas


distancias muy largas

(Medio de transporte)
Seal
Seal
Informacin

Informacin
4

Procesado de datos ssmicos

1.- FUNDAMENTOS DE SISTEMAS LINEALES


La seal puede tener diferentes formas y transportar diversos tipos
de informacin simultneamente
Onda acstica generada por las cuerdas vocales del humano
conduce palabras (informacin)
Onda electromagntica lleva informacin desde un transmisor a un
receptor de radio o televisor
Una seal puede ser funcin de una sola variable continua, como el
tiempo; o de dos o mas variables continuas (x, y, t) donde x e y son
coordenadas espaciales; o bien puede ser funcin de una o mas
variables discretas.
5

Procesado de datos ssmicos

1.- FUNDAMENTOS DE SISTEMAS LINEALES


Representacin

Procesado de datos

Manipulacin
Transformacin

seales y la
informacin
que llevan

Ejemplos:
Realzar la seal reduciendo el ruido o algunas otras
interferencias
Extraccin de cierta informacin (palabras contenidas en una
seal de dialogo, la identidad de una persona en una
fotografa o la clasificacin de un blanco en una seal de
radar)
6

Procesado de datos ssmicos

1.- FUNDAMENTOS DE SISTEMAS LINEALES


Telecomunicaciones

Modulacin de amplitud

Sistemas de
transmisin de
seales a travs de:

Tcnicas de control
automtico

Frecuencia
Cdigos de impulsos (variacin de
intensidad, de longitud o duracin, de
posicin en el tiempo, de nmero, etc.)

La energa se transmite en esquemas de


baja potencia con alto contenido de
informacin.
7

Procesado de datos ssmicos

1.- FUNDAMENTOS DE SISTEMAS LINEALES


Una seal es una funcin de una o ms variables, que transportan
informacin acerca de la naturaleza de un fenmeno fsico.
Seal de una sola variable, se le denomina seal unidimensional.
Cuando depende de dos o ms variables, la seal es
multidimensional.
Un sistema se define como una entidad que manipula una o ms
seales para llevar a cabo una funcin, produciendo de ese modo,
nuevas seales.

Seal de entrada

Sistema

Seal de salida
8

Procesado de datos ssmicos

1.- FUNDAMENTOS DE SISTEMAS LINEALES


SISTEMA, es una interconexin de operaciones que transforma una
seal de entrada, en una seal de salida con propiedades diferentes
a las de la seal de entrada.
CONTINUO

x(t)

DISCRETO

y(t)

x[n]

y(t) = H{x(t)}

y[n]

y[n] = H{x[n]}

x[n]

Sk

y[n-k]
9

Procesado de datos ssmicos

1.- FUNDAMENTOS DE SISTEMAS LINEALES

Problema de control y posicin (Chi - Tsong Chen)

10

Procesado de datos ssmicos

1.- FUNDAMENTOS DE SISTEMAS LINEALES


Problema de control de temperatura (Chi - Tsong Chen)

11

Procesado de datos ssmicos

1.- FUNDAMENTOS DE SISTEMAS LINEALES


Un sistema es lineal si cumple con las propiedades de
Homogeneidad y superposicin y adems es invariante en el tiempo.

Homogeneidad: Para una seal de entrada x[n]

ax[n]

ay[n]

Superposicin: Si la seal de entrada es:


n

xt ai xi t
i 1

12

Procesado de datos ssmicos

1.- FUNDAMENTOS DE SISTEMAS LINEALES

La seal de salida ser

y t H xt H
ai xi t

i 1

De acuerdo al principio de superposicin se puede expresar como:

y t

a y t
i

i 1

Dado que:

yi t H xi t y t

a H x t
i

i 1

13

Procesado de datos ssmicos

1.- FUNDAMENTOS DE SISTEMAS LINEALES


Invariancia con el tiempo. Un sistema es invariante en el tiempo
cuando a un retraso o corrimiento en tiempo de la seal de entrada
tambin hay el mismo retraso o corrimiento en tiempo de la seal de
salida. De otro modo el sistema es variante con el tiempo.

x(t)

x(t-to)

Sto

x(t)

yi(t)

Sto

yo(t)

14

Procesado de datos ssmicos

1.- FUNDAMENTOS DE SISTEMAS LINEALES


Si y(t) = H{x(t)} y x(t-to) = Sto{x(t)} entonces:

yi(t) = H{x(t-to)}= H{Sto{x(t)}= HSto{x(t)}

yo(t) = Sto{y(t)}=Sto{H{x(t)}}= StoH{x(t)}

Si yi(t) = yo(t) para toda seal x(t)


entonces el sistema es invariante en el tiempo si HSto= StoH
15

Procesado de datos ssmicos

1.- FUNDAMENTOS DE SISTEMAS LINEALES


Estabilidad. El sistema es estable de entrada acotada salida
acotada (BIBO) si y solo si toda entrada acotada origina una salida
acotada. La salida de este sistema no diverge si la entrada no
diverge. El operador H es estable si la seal de salida y(t) satisface
la condicin:

yt My t
Siempre que:

xt Mx t

16

Programa general

1.- FUNDAMENTOS DE SISTEMAS LINEALES


Memoria: El sistema posee memoria si la seal de salida depende
de los valores anteriores de la seal de entrada.

y[n] = 1/3 (x[n]+x[n-1]+x[n-2])

Un sistema no tiene memoria si su seal de salida depende solo del


valor presente de la seal de entrada.

y[n] = x2[n]

17

Procesado de datos ssmicos

1.- FUNDAMENTOS DE SISTEMAS LINEALES


Causalidad. El sistema es causal si el valor presente de la seal de
salida depende solo de los valores presente y/o pasado de la seal
de entrada. El sistema es no causal si depende de valores futuros de
la seal de entrada.

y[n] = 1/3 (x[n]+x[n-1]+x[n-2])

Causal

y[n] = 1/3 (x[n+1]+x[n]+x[n-1])

No causal

18

Procesado de datos ssmicos

1.- FUNDAMENTOS DE SISTEMAS LINEALES


Invertibilidad. El sistema es invertible si la entrada puede
recuperarse de la salida.

x[n]

y[n]

H-1

x[n]

H-1{y(t)} = H-1{H{x(t)}}=H-1H{x(t)}

donde H-1H= I es un operador identidad

19

Procesado de datos ssmicos

1.- FUNDAMENTOS DE SISTEMAS LINEALES

Seal de
Entrada

Funcin de
Transferencia

Seal de
Salida

. ... ... . ... ........ .. . .. . . .


.
. .. . .. . . . . .
.
.
*
* *
*
*
* *
*
+ +
+ + +
+
+
+ +
+ + +++ +

20
14

Procesado de datos ssmicos

1.- FUNDAMENTOS DE SISTEMAS LINEALES


Cintas
de
campo

Demultiplexado
Prueba de
Filtros
Filtrado
Prueba de
Deconvolucin

Anlisis de
velocidades

Prueba de
Filtros

Deconvolucin

Correcciones
Estticas
Correcciones
Dinmicas
Apilamiento

Filtrado y
normalizacin

Migracin

Filtrado y
normalizacin

21
21

Procesado de datos ssmicos

2.- RUIDOS SSMICOS


Calidad de los registros ssmicos depende del rea
Seal es todo evento de un registro ssmico que puede proporcionar
informacin
Ruido es el caso contrario e incluye eventos coherentes que
interfieren con las observaciones y medidas de las seales
Razn seal/ruido (S/R), es la relacin entre la energa de la seal y
la energa total del ruido en una regin del registro
Ruido coherente, puede aparecer en algunas trazas de manera
continua; puede ser horizontal y vertical
Ruido incoherente, no se repite en ninguna traza; se denomina
como ruido aleatorio lo que implica que no es predecible pero tiene
ciertas propiedades estadsticas.
22

Procesado de datos ssmicos

2.- RUIDOS SSMICOS


El ruido puede ser repetible y no repetible; en resumen se puede
decir que el ruido tiene tres caractersticas: Coherencia, direccin y
repetibilidad.

Caractersticas de los eventos ssmicos (Telford, Geldart, Sheriff & Keys, 1976)

23

Procesado de datos ssmicos

2.- RUIDOS SSMICOS


SEAL UTIL:

REFLEJOS
ONDAS DIRECTAS

ONDAS REFRACTADAS
ONDAS DIFRACTADAS

RUIDO:
PROPAGACIN VERTICAL

FUENTE DE ENERGA
PROPAGACIN HORIZONTAL

INSTRUMENTALES

AMBIENTAL
24

Procesado de datos ssmicos

2.- RUIDOS SSMICOS


RUIDO:

FUENTE DE ENERGA

ONDAS DE AIRE
ONDAS SUPERFICIALES

PROPAGACIN VERTICAL

RAYLEIGH
LOVE
STOHELEY
HIDRODINMICA
OTRAS ONDAS GUIADAS

ONDAS S

DIRECTA
REFRACTADA

ONDAS SECUNDARIAS

DIFRACCIONES LATERALES
EVENTOS LATERALES
REFRACCIONES
REFLEJADAS

MLTIPLES

FANTASMAS
REFRACCIONES
MLTIPLES
MLTIPLES INTERNOS

PROPAGACIN HORIZONTAL
DIFRACCIONES

25

Procesado de datos ssmicos

2.- RUIDOS SSMICOS


RUIDO:

INSTRUMENTAL
INDUCCIONES
RUIDO AMBIENTAL
EQUIPO DEFECTUOSO
FALSOS CONTACTOS
EQUIPO PERIFRICO DEFECTUOSO
MAL PLANTADOS DE DETECTORES

AMBIENTAL
TRAFICO
VIENTO
LLUVIA
OLEAJE
ANIMALES
MAQUINARIA
CORRIENTE ALTERNA

26

Procesado de datos ssmicos

2.- RUIDOS SSMICOS


Anlsis de ruido o Walkaway

Fuente: Vibroseis; geofnos espaciados 1.5 m con offset al primer gefono de 427 m. Eventos identificados: Refraccin de la
base de la capa intemperizada con 1890 m/s groundroll con 533 m/s y 625 m/s; onda de aire con 332 m/s y refraccin con
3140 m/s (Sheriff, 1973)
27

Procesado de datos ssmicos

2.- RUIDOS SSMICOS


Diferentes seales registradas

Offset (m)
0

a = onda directa; V = 650 m/s


b = refraccin en la base de la capa intemperizada;
VH = 1640 m/s

0.2

c = refraccin de una capa plana; VR = 4920 m/s

500

1000

a
b

0.4

d = reflexin de refractor en c; Vav = 1640 m/s

f = reflexin de una capa plana; Vav = 2300 m/s


g = reflexin de una capa inclinada; Vav = 2630 m/s
h = mltiple de d;
i = mltiple de e;
j = ground roll; V = 575 m/s
k = onda de aire; V = 330 m/s

d
Tiempo de arribo (s)

e = reflexin de una capa plana; Vav = 1970 m/s

1500

0.6

h
0.8

1.0

k
1.2

a
j
m
f

1.4

l = refraccin reflejada en la interrupcin del


refractor de c a lo largo de la lnea
1.6

m = refraccin reflejada en la interrupcin del


refractor de c en la parte transversal de la lnea

l
g
i

1.8

28

Procesado de datos ssmicos

2.- RUIDOS SSMICOS


Tiro ssmico con ruido

29

Procesado de datos ssmicos

2.- RUIDOS SSMICOS


Tipos de mltiples

30

Procesado de datos ssmicos

2.- RUIDOS SSMICOS


Seccin ssmica con mltiples

31

Procesado de datos ssmicos

2.- RUIDOS SSMICOS


Tiro de Dinamita en un tendido simtrico de 120 canales

32

Procesado de datos ssmicos

2.- RUIDOS SSMICOS


Tiro de Vibroseis en un tendido simtrico de 100 canales

Barrido de 10 s con canales auxiliares a


la izquierda que contienen el barrido de
referencia.

Punto de tiro con correlacin cruzada


para obtener los 3 s de informacin.

33

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


Causas de ruido ssmico
Efecto del cable
Ruido
ambiental

Sismgrafo

Acoplamiento de gefonos
Respuesta en frecuencia de gefonos
Efecto del arreglo
Ground Roll

Punto de
tiro

Fantasma en
superficie libre

Capa de baja
velocidad
Refracciones

Efecto de la fuente

Divergencia
esfrica

Perdida por
transmisin

Refracciones

Factor Q

Dispersin

Mltiples de
perodo corto
Mltiples de
perodo largo

34

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


Factores que afectan la amplitud ssmica

Superficie terrestre

Absorcin
Entrada

Zona de
explosin

Coeficientes de
reflexin y refraccin

Receptor

Salida

Divergencia esfrica

35

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


Factores que afectan la amplitud ssmica

36

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


Factores que modifican la amplitud ssmica
La Dispersin es el decremento de la Densidad e Intensidad de Energa por
el cambio de frecuencia debido a condiciones Geolgicas (contacto de las
partculas del medio transmisor) y/o Geomtricas (por divergencia esfrica)
Densidad de Energa
Es la energa por volumen unitario en las proximidades de un punto, para
una onda armnica P esfrica, la energa cintica (Ec) por unidad de volumen
(v) est dada por:

Ec 1
1
2
u 2 A2 sen 2 t
v
2
2
1 22 A2 2 2 f 2 2
donde: es la frecuencia angular, es la fase, t el tiempo, es la densidad;
f es la frecuencia y A es la amplitud

37

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


Factores que modifican la amplitud ssmica
Intensidad de Energa

Se define como la cantidad de energa que fluye a travs de un rea unitaria


normal a la direccin de propagacin de la onda en unidad de tiempo.

I EV 1 2V 2 A2
donde: V es la velocidad
El flujo de energa por unidad de tiempo en direccin radial, es el producto de
la intensidad y el rea:

I1 A1 I 2 A2
Las reas A1 y A2 son proporcionales al cuadrado de su radio, entonces:

I 2 I1 A1

A2 r1 r2 E2 E1
2

38

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


Factores que modifican la amplitud ssmica
Ya que E es proporcional a I ,
esto significa que la expansin
geomtrica provoca que la
intensidad y densidad de energa
de la onda esfrica disminuyan
con respecto al cuadrado de la
distancia a la fuente (Newman,
1973). A este fenmeno se le
denomina, divergencia esfrica.
El frente de onda es esfrico si la
velocidad ssmica del subsuelo
es constante, es decir el rayo es
recto.

39
45

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


Factores que modifican la amplitud ssmica
Absorcin. La disminucin de la densidad de energa de una onda, se debe
tambin a otros dos mecanismos: la absorcin y la particin de energa en
diferentes interfases. Cuando la onda pasa a travs de un medio; la energa
elstica asociada con el movimiento de onda se absorbe gradualmente por el
medio, transformndose en calor hasta desaparecer el movimiento de onda.
La perdida de energa por absorcin es exponencial respecto a la distancia,
para ondas elsticas en las rocas.

I I 0e x
Donde I e I0 son los valores de intensidad en dos puntos separados una
distancia x y el valor es el coeficiente de absorcin.

En laboratorio se encontr que el coeficiente de absorcin es proporcional


a la frecuencia f , y representa un mecanismo para explicar la perdida de
altas frecuencias con la distancia.
40

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


Factores que modifican la amplitud ssmica
Las perdidas por divergencia son mas importantes que las perdidas por
absorcin en bajas frecuencias y distancias cortas.

Perdida de energa por absorcin y divergencia ( 0.25 db/l y


V = 2 000 m/seg
Frecuencia

Absorcin

Divergencia

Distancia desde el punto de tiro xs

1200 m

2200 m

4200 m

8200 m

1 Hz

0.55 db

1.1 db

2.2 db

4.4 db

1.6

3.3

6.6

13

10

5.5

11

22

44

30

16

33

66

130

100

55

110

220

440

todas

16

21

26

32
41

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


Factores que modifican la amplitud ssmica
Dispersin y velocidad de grupo. La velocidad que aparece en la ecuacin
de onda se refiere a la velocidad de fase, ya que es la distancia que recorre
un punto de fase constante por unidad de tiempo (pico o valle). Pero no es
necesariamente la misma velocidad con que se transmite la energa, a la cual
se le denomina velocidad de grupo y se denota por U como se ilustra en la
figura.

42
45

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


Factores que modifican la amplitud ssmica
Si la velocidad es la misma en todas las frecuencias el pulso permanecer
con la misma velocidad de fase y grupo. Si la velocidad cambia con la
frecuencia, el pulso tambin cambia cuando viaja y la velocidad de grupo es
diferente de la velocidad de fase en un medio dispersivo.

dV
dV
U V l
V
dl
d
Donde V, l, , dV/dl y dV/d son valores promedio en el rango de
frecuencias del que esta formado el pulso.

43

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


Factores que modifican la amplitud ssmica
Si la velocidad es la misma en todas las frecuencias el pulso permanecer
con la misma velocidad de fase y grupo. Si la velocidad cambia con la
frecuencia, el pulso tambin cambia cuando viaja y la velocidad de grupo es
diferente de la velocidad de fase en un medio dispersivo.

dV
dV
U V l
V
dl
d
Donde V, l, , dV/dl y dV/d son valores promedio en el rango de
frecuencias del que esta formado el pulso.

44

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


Comparacin entre velocidad de grupo y de fase

Velocidad de grupo

Velocidad de fase

x
U
t g

x
V
t p

45

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


Ley de Reflexin

i r

i r
V1
V2

Ley de Refraccin

seni senR

V1
V2
El ngulo i es crtico cuando
el ngulo R mide 90

V1
senic
V2

i
V1

V2
R
46

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


La amplitud incidente es la suma de las amplitudes reflejada y
transmitida

i r t
Las amplitudes de las ondas reflejadas y transmitidas varan con el
ngulo de incidencia, dadas por las ecuaciones de Zoeppritz.
Ai

Ar

V1
V2
R

At
47

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


Para un ngulo de incidencia igual acero, las amplitudes se
relacionan en coeficientes de reflexin y transmisin:
Coeficiente de reflexin:

A1 Z 2 Z1

A0 Z 2 Z1

Coeficiente de transmisin:

A2
2Z1
T

A0 Z 2 Z1
Z es la impedancia acstica (onda P) de la capa, y est dada por Z=V,
dnde el V es la velocidad de la onda y la densidad.
Las mismas formulas aplican a las ondas S de incidencia normal.

48

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


Coeficientes de reflexin y transmisin
Si

A1 A0 R

entonces

A2 1 R A1

La amplitud de la onda que se transmite en el subsuelo va perdiendo parte


de su energa por reflexin en cada interfase.
R1(1 R02)

Si A0 = 1
1

R0

R1(1 R0)(1+ R0)


R0

(1-R0)

(1-R0)R1
R1

(1-R0) - R1(1-R1)
R2
49

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


FACTORES DE ATENUACIN DE LA SEAL SSMICA

AFECTAN LA
AMPLITUD

AFECTAN LA
AMPLITUD, LA
FRECUENCIA Y LA
FASE

PERDIDAS POR
TRANSMISIN
PERDIDAS POR DIVERGENCIA
ESFRICA

ABSORCIN
INELSTICA
PERDIDAS POR
DISPERSIN
50

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


FACTORES DE ATENUACIN DE LA SEAL SSMICA

Amplitud promedio de la seal registrada


Amplitud promedio de los primeros quiebres
Prdidas por reflexin
Prdidas por 1divergencia esfrica

cft
Atenuacin e

Reflexin fuerte

Amplitud del promedio de la seal


ssmica

Reflexin dbil

TIEMPO
51

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


FACTORES DE ATENUACIN DE LA SEAL SSMICA

A
M
P
L
I
T
U
D

Amplitud real de la seal registrada


(traza ssmica)

TIEMPO

52

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


NORMALIZACIN EXPONENCIAL
B
A
M
P
L
I
T
U
D

1 ct
e
t

TIEMPO

G
A
N
A
N
C
I
A

1 ct
te
B
Curva de ganancia

TIEMPO
53

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


NORMALIZACIN RMS
Traza original

x 2x

Factores centrales
Por ventana

H1

H2

4x

H3

4x

H4

Traza de salida

RMS 2
Factor
RMS1
RMS 2 Deseado

RMS1 Existente
54

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


NORMALIZACIN DE AMPLITUDES

55

Procesado de datos ssmicos

3.- AMPLITUD SSMICA


NORMALIZACIN RMS

56

Procesado de datos ssmicos

4.- FILTROS DIGITALES


ESPECTROS DE AMPLITUD Y FASE

57

Procesado de datos ssmicos

4.- FILTROS DIGITALES


ESPECTROS DE AMPLITUD Y FASE

58

Procesado de datos ssmicos

4.- FILTROS DIGITALES


ESPECTROS DE AMPLITUD Y FASE

59

Procesado de datos ssmicos

4.- FILTROS DIGITALES


SEAL DIGITAL Y FRECUENCIA DE NYQUIST

60

Procesado de datos ssmicos

4.- FILTROS DIGITALES


INTERVALO DE MUESTREO

61

Procesado de datos ssmicos

4.- FILTROS DIGITALES


INTERVALO DE MUESTREO

62

Procesado de datos ssmicos

4.- FILTROS DIGITALES


ESPECTROS DE AMPLITUD Y FASE

63

Procesado de datos ssmicos

4.- FILTROS DIGITALES


ESPECTROS DE AMPLITUD Y FASE

64

Procesado de datos ssmicos

4.- FILTROS DIGITALES


ESPECTROS DE AMPLITUD Y FASE

65

Procesado de datos ssmicos

4.- FILTROS DIGITALES


ESPECTROS DE AMPLITUD Y FASE

66

Procesado de datos ssmicos

4.- FILTROS DIGITALES


ESPECTROS DE AMPLITUD Y FASE

67

Procesado de datos ssmicos

4.- FILTROS DIGITALES


ESPECTROS DE AMPLITUD Y FASE

68

Procesado de datos ssmicos

4.- FILTROS DIGITALES


FILTRO ANLOGO

69

Procesado de datos ssmicos

4.- FILTROS DIGITALES


FILTRO DIGITAL

70

Procesado de datos ssmicos

4.- FILTROS DIGITALES


CONVOLUCIN

71

Procesado de datos ssmicos

4.- FILTROS DIGITALES


CONVOLUCIN

C1 = A1B1
C2 = A1B2 + A2B1
C3 = A1B3 + A2B2 + A3B1

F = (C1 + C2 + . CZ)

C4 = A1B4 + A2B3 + A3B2 + A4B1

F = (A1 + A2 + AN)(B1 + B2 + . + BM)

F = (A1 + A2 + AN) (Respuesta al Impulso)

.
CZ = ANBM

72

Procesado de datos ssmicos

4.- FILTROS DIGITALES


FILTROS DIGITALES EN TIEMPO Y FRECUENCIA

73

Procesado de datos ssmicos

4.- FILTROS DIGITALES


DIFERENTES FILTROS DIGITALES

74

Procesado de datos ssmicos

4.- FILTROS DIGITALES


DIFERENTES FILTROS DIGITALES

75

Procesado de datos ssmicos

4.- FILTROS DIGITALES


FILTRADO PASA BANDA PARA ATENUAR RUIDO INCOHERENTE

Anlisis de Fourier
para determinar el
nivel de ruido
incoherente, y
establecer el rango de
frecuencias de la seal
ssmica de reflexin.

76

Procesado de datos ssmicos

4.- FILTROS DIGITALES


FILTRADO PASA BANDA PARA ATENUAR RUIDO INCOHERENTE

Aplicar filtro pasa banda para eliminar


ruido incoherente y
mejorar la seal.

77

Procesado de datos ssmicos

4.- FILTROS DIGITALES


DIFERENTES BANDAS DE FRECUENCIA

78

Procesado de datos ssmicos

4.- FILTROS DIGITALES


EJEMPLOS DE FILTRADO

79

Procesado de datos ssmicos

4.- FILTROS DIGITALES


EJEMPLOS DE FILTRADO

80

Procesado de datos ssmicos

4.- FILTROS DIGITALES


FILTRO f k PARA RUIDO COHERENTE
El Filtro f-k,
consiste en
transformar
los datos del
dominio x-t
(distancia vs
tiempo) de los
sismogramas
al dominio f-k
(frecuencia vs
nmero de
onda), donde
son fciles de
distinguir los
ruidos
coherentes.
81

Procesado de datos ssmicos

4.- FILTROS DIGITALES


FILTRO f k PARA RUIDO COHERENTE

Aplicacin del
Filtro f-k, se
limita un
abanico similar
al pasa bandas
que elimina los
ruidos
coherentes.

82

Procesado de datos ssmicos

5.- AUTOCORRELACIN/CORRELACIN CRUZADA


Esperanza
La esperanza E se refiere a un promedio estadstico que aplicado a una
ondcula proporciona su energa. Es la autocorrelacin de la ondcula sin
desplazamiento t 0.

E b bt2 b02 b12 b22 ... energa


2
t

t 0

Para una serie de tiempo se obtiene la potencia con la autocorrelacin sin


desplazamiento t = 0.

1
E u lm
u2T u2T 1 .... uT2 potencia
T 2T 1
2
t

83

Procesado de datos ssmicos

5.- AUTOCORRELACIN/CORRELACIN CRUZADA


AUTOCORRELACIN Y CORRELACIN CRUZADA
La autocorrelacin para todos los valores de desplazamiento t proporciona
informacin de la energa o potencia en diferentes componentes de
frecuencia.

bb t bt t bt Ebt t bt
uu t Eut t ut

para la ondcula

para la serie de tiempo

La correlacin cruzada es una forma de conocer la similitud entre dos seales


y a diferencia de la autocorrelacin, no es conmutativa y depende de cual
seal est fija. Para las seales xt y yt la correlacin cruzada es:

xy t Ext t yt
84

Procesado de datos ssmicos

6.- FILTRO PTIMO DE WIENER


FILTRO DE MNIMOS CUADRADOS
Sea ft, t=0,1,2,(n-1) la respuesta al impulso del filtro deseado (secuencia
finita). La salida ct para una secuencia de entrada estocastica {xt} es:
n 1

ct f s xt s t 0,1,2,......
s 0

El diseo requiere obtener un {ft} tal que {ct} es el mnimo cuadrado de la


seal deseada {dt}. Para ello se minimiza la energa o potencia Jn del error (et
= dt ct), donde n es la longitud del filtro.

J n E e E dt ct .
2
t

85

Procesado de datos ssmicos

6.- FILTRO PTIMO DE WIENER


FILTRO DE MNIMOS CUADRADOS
Minimizando Jn, respecto a cada coeficiente del filtro desconocido de manera
que :

et
J n
E 2et
0 k 0,1,2,..., (n 1)

f k
f k
Lo que se reduce a una ecuacin matricial:

R xxf rxd

Donde Rxx es la matriz de correlaciones y rxd el vector columna de la


correlacin cruzada entre la seal deseada y la seal de entrada:

f R xx1rxd

El filtro de Wiener f es ptimo porque el error


de mnimos cuadrados entre la salida actual y
la deseada es mnimo:
86

Procesado de datos ssmicos

6.- FILTRO PTIMO DE WIENER


FLUJO DEL DISEO DEL FILTRO PTIMO DE WIENER

Entrada

Salida deseada

Autocorrelacin

Correlacin cruzada

Filtro de Wiener
*

Salida actual

87

Procesado de datos ssmicos

6.- FILTRO PTIMO DE WIENER


FLUJO DE INTERRELACIN ENTRE VARIOS FILTROS DE
DECONVOLUCIN
Optimum Wiener
Filters

Least
Squares
Inverse Filter

Desired Output
Inverse
Filter
Zero-Delay
Spike

Time-Advanced
Version of Input with
Prediction Lag, a

Unit Prediction
Lag (a = 1)

Spiking
Deconvolution

Predictive
Deconvolution

Any Other Form

Wavelet Shaping

88

Procesado de datos ssmicos

6.- FILTRO PTIMO DE WIENER


ELIMINACIN DE MLTIPLES

89

Procesado de datos ssmicos

6.- FILTRO PTIMO DE WIENER


DECONVOLUCIN
El proceso de Deconvolucin, consiste en regenerar la forma de onda
emitida al subsuelo, con la finalidad de atenuar todos los efectos que sufre la
energa ssmica al paso por el subsuelo.

90

Procesado de datos ssmicos

6.- FILTRO PTIMO DE WIENER


DECONVOLUCIN
Deconvolucin Spiking

reverberacin de
energa

91

Procesado de datos ssmicos

6.- FILTRO PTIMO DE WIENER


DECONVOLUCIN
Deconvolucin Spiking que suprime la reverberacin de energa que
acompaa las reflexiones de los gathers de tiro comn en la lmina anterior

92

Procesado de datos ssmicos

7.- MODELO CONVOLUCIONAL


SERIE REFLECTIVA
Al conjunto de coeficientes de reflexin R, se le conoce como la Serie
Reflectiva propia del subsuelo.

Los coeficientes de reflexin


pueden ser graficados como
una funcin de tiempo creando
un series de tiempo llamada
Serie Reflectiva

93

Procesado de datos ssmicos

7.- MODELO CONVOLUCIONAL


TRAZA SSMICA
La traza ssmica es generada por la convolucin de la serie reflectiva del
subsuelo con la ondcula generada por la fuente de energa en superficie.

94

Procesado de datos ssmicos

7.- MODELO CONVOLUCIONAL


GEOMETRA DE ADQUISICIN
Debido al gran volumen de datos grabados en campo, se pueden organizar
en conjuntos de trazas o gathers.
Un gather de tiro comn es una
coleccin de trazas ssmicas
grabada por varios receptores de un
solo tiro. sta es la configuracin en
que los datos ssmicos son
adquiridos en el campo.

Un gather de receptor comn es


una coleccin de trazas ssmicas que
corresponden a varios tiros grabados
en un solo receptor.

95

Procesado de datos ssmicos

7.- MODELO CONVOLUCIONAL


GEOMETRA DE ADQUISICIN

Un gather punto medio comn


(CMP) es una coleccin de trazas en
que el tiro y el receptor estn
distribudos simtricamente sobre la
misma posicin del punto medio.

El punto de profundidad comn


(CDP) es el punto en una interface
horizontal plana donde se generan
todas las reflexiones en un CMP.

96

Procesado de datos ssmicos

8.- CORRECIONES ESTTICAS


CORRECCIN ESTTICA
Las correcciones verticales invariantes en tiempo, aplicadas para colocar
fuentes y detectores sobre un mismo plano horizontal imaginario (Datum), se
le conoce como Correccin Esttica.
Tiene como objetivo cancelar los efectos ssmicamente indeseables de la
porcin superficial de la Tierra, que afectan a todas las ondas que viajan a
travs de los estratos superficiales.

97

Procesado de datos ssmicos

8.- CORRECIONES ESTTICAS


CORRECCIN ESTTICA
Estos efectos son debidos a variaciones de elevaciones, variaciones errticas
del estrato superficial (capa de intemperismo), sean de espesor o de cambios
laterales de velocidad.

Cest.

Elev . fuente Elev .NR


V

Elev .receptor Elev NR


V

Elev . fuente Esp.w Elev .w Elev .NR

Cest.

Vw
V1

Elev .receptor Esp.w Elev .w Elev .NR

Vw
V1

98

Procesado de datos ssmicos

8.- CORRECIONES ESTTICAS


CORRECCIN ESTTICA

99

Procesado de datos ssmicos

9.- CORRECCIN DINMICA


CORRECCIN DINMICA NMO
La Correccin Dinmica consiste en llevar todas las trayectorias oblicuas,
de una familia CMP, a la vertical, o sea, suponer que la fuente y el receptor
se encuentran en la misma estacin.

100

Procesado de datos ssmicos

9.- CORRECCIN DINMICA


CORRECCIN DINMICA NMO
La Correccin Dinmica est definida por:

T
T
CD t TX T0 2 X 0
2 2

del triangulo rectngulo se tiene:

donde VC es considerada VRMS

2VC

T0 TX

2 2
2

x

VC

x
t TX T0 T02
VC

T02 TX2

T0

101

Procesado de datos ssmicos

9.- CORRECCIN DINMICA


CORRECCIN DINMICA NMO
La Correccin Dinmica requiere conocer estadsticamente las velocidades
de correccin.

102

Procesado de datos ssmicos

9.- CORRECCIN DINMICA


CORRECCIN DINMICA NMO
Un Gather CDP sin y con Correccin Dinmica

103

Procesado de datos ssmicos

10.- ANLISIS DE VELOCIDAD


ANLISIS DE VELOCIDAD PARA CORRECCIN DINMICA
Efecto de la aplicacin de velocidades incorrectas

104

Procesado de datos ssmicos

10.- ANLISIS DE VELOCIDAD


ANLISIS DE VELOCIDAD PARA CORRECCIN DINMICA
La correccin dinmica como proceso para la eliminacin de eventos
mltiples en los gathers CDPs

105

Procesado de datos ssmicos

11.- APILAMIENTO DE TRAZAS


APILAMIENTO DE TRAZAS
El proceso de Apilamiento consiste en sumar algebraicamente las diferentes
trazas individuales que conforman un mismo gather CDP

106

Procesado de datos ssmicos

11.- APILAMIENTO DE TRAZAS


APILAMIENTO DE TRAZAS
Seccin ssmica Apilada

107

Procesado de datos ssmicos

12.- VELOCIDADES SSMICAS


CLASIFICACIN TIPO DE MEDICIN DE VELOCIDADES

Estudios
convencionales
de pozos

Disparo en pozos

Mtodo X2 T2

Mtodo T T

Registro Snico
(BHC)

108

Procesado de datos ssmicos

12.- VELOCIDADES SSMICAS


ESTUDIOS CONVENCIONALES DE POZOS
Disparo en pozos

Z
t to 2
2
z

Donde:
t0, es el tiempo observado

t, es el tiempo de trnsito vertical

z
V
t
V, es la velocidad promedio

109

Procesado de datos ssmicos

12.- VELOCIDADES SSMICAS


ESTUDIOS CONVENCIONALES DE POZOS
Disparo en pozos

La velocidad promedio en el intervalo


(zm y zn), para dos disparos

zm zn
Vi
tm tn
Desventajas:
- Caro
- Posibles daos al pozo
110

Procesado de datos ssmicos

12.- VELOCIDADES SSMICAS


ESTUDIOS CONVENCIONALES DE POZOS
Registro Snico (BHC)

Fuente 1: Receptor 2, Receptor 4


Fuente 2: Receptor 3, Receptor 1
La distancia entre receptores es de 2 pies
t1 t1

1
t

t2 t2

t1 t2 s
(
)
pie
2

(Velocidad de onda P)

111

Procesado de datos ssmicos

12.- VELOCIDADES SSMICAS


MTODO X2 - T2
Se usa cuando no se tienen muchas
distancias offset (x).

Se utilizan offsets largos,


basados en el arreglo descrito por Dix.
2
x
t 2 2 t02
Vrms

112

Procesado de datos ssmicos

12.- VELOCIDADES SSMICAS


MTODO X2 - T2
Un estudio de este tipo nos proporciona velocidades precisas si:
1) Si se tiene un nmero moderado de reflexiones
2) Se aplican correcciones superficiales
3) El trabajo de campo y la interpretacin se hacen cuidadosamente.
Una vez que se han determinado
las velocidades (Vrms) sucesivas
a dos reflectores paralelos en
base a la ecuacin:

2
x
t
2
0

n 1

Se encuentran las velocidades


de intervalo usando la frmula de
Dix.

2
L

(V
Vn2

2
Vrms

t V t )
2
U

tn
113

Procesado de datos ssmicos

12.- VELOCIDADES SSMICAS


MTODO T - T
Este mtodo se basa en el Normal Moveout (NMO)

El NMO se puede rescribir como:

x
V
2t0tn

Para calcular tn en un tendido simtrico se emplea la siguiente


frmula:

tn

1
t1 t0 tl t0 1 t1 tl t0
2
2
114

Procesado de datos ssmicos

13.- FACTORES QUE AFECTAN LAS VELOCIDADES


ISOTROPIA Y ANISOTROPIA

115

Procesado de datos ssmicos

13.- FACTORES QUE AFECTAN LAS VELOCIDADES


MEDIDAS DE VELOCIDAD EN DIFERENTES ROCAS (Birch (1966)

116

Procesado de datos ssmicos

13.- FACTORES QUE AFECTAN LAS VELOCIDADES


RELACIN VELOCIDAD DENSIDAD EN DIFERENTES ROCAS
(Gardner et al., (1974)

117

Procesado de datos ssmicos

13.- FACTORES QUE AFECTAN LAS VELOCIDADES


RELACIN VELOCIDAD PROFUNDIDAD DE SEPULTAMIENTO
(Jankowsky., (1970)

118

Procesado de datos ssmicos

14.- CONCEPTO DE MIGRACIN SSMICA


POSICIONAMIENTO REAL DE LOS REFLECTORES
Debido a que las diferentes trayectorias de las ondas se reflejan en incidencia
normal se realiza el proceso de Migracin para llevar los eventos a su
verdadera posicin en el subsuelo

119

Procesado de datos ssmicos

14.- CONCEPTO DE MIGRACIN SSMICA


MTODO GEOMTRICO DE MIGRACIN

120

Procesado de datos ssmicos

14.- CONCEPTO DE MIGRACIN SSMICA


EJEMPLOS DE PROCESO DE MIGRACIN

Seccin Apilada

Seccin Migrada

121

Procesado de datos ssmicos

14.- CONCEPTO DE MIGRACIN SSMICA


ESPECTROS DE AMPLITUD Y FASE

Seccin Final Normal

Seccin Final Migrada

122

Procesado de datos ssmicos

15.- SECCIN EN TIEMPO


SECCIN FINAL EN TIEMPO

123

Procesado de datos ssmicos

16.- SECCIN EN PROFUNDIDAD


SECCIN FINAL EN PROFUNDIDAD

124

Procesado de datos ssmicos

FINALMENTE
EL Software de procesado ssmico, contiene programas y/o
subrutinas que realizan un proceso especfico

Normalmente un buen paquete de programas contienen alrededor de 250 programas


125

También podría gustarte