Procesado de Datos Sísmicos PDF
Procesado de Datos Sísmicos PDF
Procesado de Datos Sísmicos PDF
DATOS SSMICOS
Fecha: 20.03.2006
Versin No.1
Pginas:125
Edicin No. 1
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
La comunicacin es el
intercambio de informacin
(Medio de transporte)
Seal
Seal
Informacin
Informacin
4
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
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.)
Seal de entrada
Sistema
Seal de salida
8
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
10
11
ax[n]
ay[n]
xt ai xi t
i 1
12
y t H xt H
ai xi t
i 1
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
x(t)
x(t-to)
Sto
x(t)
yi(t)
Sto
yo(t)
14
yt My t
Siempre que:
xt Mx t
16
Programa general
y[n] = x2[n]
17
Causal
No causal
18
x[n]
y[n]
H-1
x[n]
H-1{y(t)} = H-1{H{x(t)}}=H-1H{x(t)}
19
Seal de
Entrada
Funcin de
Transferencia
Seal de
Salida
20
14
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
Caractersticas de los eventos ssmicos (Telford, Geldart, Sheriff & Keys, 1976)
23
REFLEJOS
ONDAS DIRECTAS
ONDAS REFRACTADAS
ONDAS DIFRACTADAS
RUIDO:
PROPAGACIN VERTICAL
FUENTE DE ENERGA
PROPAGACIN HORIZONTAL
INSTRUMENTALES
AMBIENTAL
24
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
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
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
Offset (m)
0
0.2
500
1000
a
b
0.4
d
Tiempo de arribo (s)
1500
0.6
h
0.8
1.0
k
1.2
a
j
m
f
1.4
l
g
i
1.8
28
29
30
31
32
33
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
Superficie terrestre
Absorcin
Entrada
Zona de
explosin
Coeficientes de
reflexin y refraccin
Receptor
Salida
Divergencia esfrica
35
36
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
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
39
45
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.
Absorcin
Divergencia
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
42
45
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
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
Velocidad de grupo
Velocidad de fase
x
U
t g
x
V
t p
45
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
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
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
A1 A0 R
entonces
A2 1 R A1
Si A0 = 1
1
R0
(1-R0)
(1-R0)R1
R1
(1-R0) - R1(1-R1)
R2
49
AFECTAN LA
AMPLITUD
AFECTAN LA
AMPLITUD, LA
FRECUENCIA Y LA
FASE
PERDIDAS POR
TRANSMISIN
PERDIDAS POR DIVERGENCIA
ESFRICA
ABSORCIN
INELSTICA
PERDIDAS POR
DISPERSIN
50
cft
Atenuacin e
Reflexin fuerte
Reflexin dbil
TIEMPO
51
A
M
P
L
I
T
U
D
TIEMPO
52
1 ct
e
t
TIEMPO
G
A
N
A
N
C
I
A
1 ct
te
B
Curva de ganancia
TIEMPO
53
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
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
C1 = A1B1
C2 = A1B2 + A2B1
C3 = A1B3 + A2B2 + A3B1
F = (C1 + C2 + . CZ)
.
CZ = ANBM
72
73
74
75
Anlisis de Fourier
para determinar el
nivel de ruido
incoherente, y
establecer el rango de
frecuencias de la seal
ssmica de reflexin.
76
77
78
79
80
Aplicacin del
Filtro f-k, se
limita un
abanico similar
al pasa bandas
que elimina los
ruidos
coherentes.
82
t 0
1
E u lm
u2T u2T 1 .... uT2 potencia
T 2T 1
2
t
83
bb t bt t bt Ebt t bt
uu t Eut t ut
para la ondcula
xy t Ext t yt
84
ct f s xt s t 0,1,2,......
s 0
J n E e E dt ct .
2
t
85
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
f R xx1rxd
Entrada
Salida deseada
Autocorrelacin
Correlacin cruzada
Filtro de Wiener
*
Salida actual
87
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
Wavelet Shaping
88
89
90
reverberacin de
energa
91
92
93
94
95
96
97
Cest.
Cest.
Vw
V1
Vw
V1
98
99
100
T
T
CD t TX T0 2 X 0
2 2
2VC
T0 TX
2 2
2
x
VC
x
t TX T0 T02
VC
T02 TX2
T0
101
102
103
104
105
106
107
Estudios
convencionales
de pozos
Disparo en pozos
Mtodo X2 T2
Mtodo T T
Registro Snico
(BHC)
108
Z
t to 2
2
z
Donde:
t0, es el tiempo observado
z
V
t
V, es la velocidad promedio
109
zm zn
Vi
tm tn
Desventajas:
- Caro
- Posibles daos al pozo
110
1
t
t2 t2
t1 t2 s
(
)
pie
2
(Velocidad de onda P)
111
112
2
x
t
2
0
n 1
2
L
(V
Vn2
2
Vrms
t V t )
2
U
tn
113
x
V
2t0tn
tn
1
t1 t0 tl t0 1 t1 tl t0
2
2
114
115
116
117
118
119
120
Seccin Apilada
Seccin Migrada
121
122
123
124
FINALMENTE
EL Software de procesado ssmico, contiene programas y/o
subrutinas que realizan un proceso especfico