M´etodos Geof´ısicos: Parte II M´etodo S´ısmico de Refracci´on M´etodo S´ısmico de Reflexi´on Jairo Torres
[email protected] F´ısico, Universidad Industrial de Santander Candidato a Mag´ıster en Ciencias de la Tierra, Universidad EAFIT Escuea de Geolog´ıa Universidad Industrial de Santander
20 de febrero de 2012
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
1 / 65
Contenido de la Presentaci´on 1
Libro de Referencia y Material para el II Corte
2
Interpretaci´on de S´ısmica de Refracci´ on Interpretaci´on de Modelos Una Interface Horizontal Simple Una Interface Horizontal Simple Interfaces Horizontales Varias Interfaces Inclinada Simple
3
S´ısmica de Refracci´on: Adquisici´ on y Procesamiento Conceptos B´asicos
4
Tipos de Velocidad
5
Tipos de Velocidad
6
Tipos de Velocidad Interpretaci´on Tectonica de Perfiles de S´ısmica de Reflexi´on Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
2 / 65
Libro de Referencia y Material para el II Corte
Fecha II Parcial Marzo 09 de 2012 Capitulo 5 de M.Slawisnki (Simetr´ıas Is´ otropa, VTI, HTI, TTI) Capitulo 3, Capitulo 4, Capitulo 5 y 6 Whole Earth Geophysics (Robert. J Lillie) Articulos de Referencia: I I I I I I
Levantamiento de S´ısmica de Pozos, M´as alla del Perfil Vertical. Elevaci´ on de estandar en los datos s´ısmicos Inversi´ on de Datos S´ısmicos: Lectura entre L´ıneas La recompensa que yace bajo la Sal. La Naturaleza de los Yacimientos Naturalmente Fracturados La era de las im´agenes en escala de Profundidad
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
3 / 65
Interpretaci´on de S´ısmica de Refracci´on
La refracci´on de ondas s´ısmicas provee contrastes en la composici´on de las estructruas de ciertas parte del subsuelo por que las ondas refractadas encuentran cambios en la velocidad s´ısmica. Los cambios en la velocidad estan relacionados a cambios en los modulos de rigidez, volumetricos y la densidad. El m´etodo de refracci´ on es el mas usado en los casos que exista un incremento abrupto en la velocidad con la profundidad, por que la onda P refractada eventualmente adelante a las otras ondas.
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
4 / 65
Interpretaci´on de S´ısmica de Refracci´on
2 problemas tratados efectivamente por la s´ısmica de refracci´on. Crustal Thickness. (frontera entre la corteza y el manto ) Contraste fuerte en velocidades entre material blando y material compuesto por Roca firme.
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
5 / 65
Interpretaci´on de S´ısmica de Refracci´on
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
6 / 65
Interpretaci´on de S´ısmica de Refracci´on
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
7 / 65
Interpretaci´on de S´ısmica de Refracci´on
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
8 / 65
Interpretaci´on de S´ısmica de Refracci´on Interpretaci´ on de Modelos
Una interface horizontal Varias Interfaces Horizontales Una Interface Inclinada Para estos modelos podremos deducir el problema directo considerando Modelo directo, anticipando curvas de travel-time o por medio de inversi´ on de campos de velocidad.
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
9 / 65
Interpretaci´on de S´ısmica de Refracci´on Una Interface Horizontal Simple
Para un modelo de una sola interface horizontal, el tiempo de arribo de la onda compresional de refracci´ on cr´ıtica esta dada por la expresi´on t = t1 +
x v2
Donde: t1 : es el intercepto en el eje vertical (tiempo) t1 =
2h cos(θc ) v1
θc , el ´angulo cr´ıtico. Distancia critica xc = 2h tan(θc ) h, es la profundidad v1 , la velocidad de la onda P en la capa superior v2 , la velocidad de la onda P en la capa inferior xcr , la distancia que corta la onda de refracc´ on con la onda directa r v2 + v1 xcr = 2h v2 − v1 Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
10 / 65
Interpretaci´on de S´ısmica de Refracci´on Una Interface Horizontal Simple
Podemos en base a las curvas interpretadas buscar los valores, usando las pendientes en las curvas para la onda directa y refractada:
Onda Directa mwd =
1 v1
→
v1 =
1 mwd
v1 =
1 mhw
Onda Refractada mhw =
1 v1
→
θc y h θc = sin−1 (v1 /v2 )
Jairo Torres (Geo UIS)
t1 =
2h cos(θc ) → v1 M´ etodos Geof´ısicos
h
20 de febrero de 2012
11 / 65
Interpretaci´on de S´ısmica de Refracci´on Una Interface Horizontal Simple
Comparaci´on de tiempos de intercepto t1 y distancia xcr para diferentes espesores de la corteza.
a) Espesores en regiones tipicas de corteza continental y b) Corteza oce´anica (zonas de rift continental; m´argenes continental pasivos y cuencas oce´anicas) Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
12 / 65
Interpretaci´on de S´ısmica de Refracci´on Una Interface Horizontal Simple
Comparaci´on de tiempos de intercepto t1 y distancia xcr para diferentes espesores de la corteza.
c) Espesores en regiones tipicas debajo de monta˜ nas. Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
13 / 65
Interpretaci´on de S´ısmica de Refracci´on Una Interface Horizontal Simple
Depth to Bedrock
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
14 / 65
Interpretaci´on de S´ısmica de Refracci´on Una Interface Horizontal Simple
Depth to Bedrock
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
15 / 65
Interpretaci´on de S´ısmica de Refracci´on: Varias Int Horizontales Vamos a considerar el caso de varias interfaces horizontales en donde la velocidad de profundidad de las capas siempre es mayor entonces v1 < v2 < v3 < v4 El ´angulo para la primera refracci´ on critica es: v1 sin(θ1,2 ) = v2 El ´angulo para la segunda refracci´ on critica es: v1 sin(θ1,4 ) = v4 v2 sin(θ2,4 ) = v4 v3 sin(θ3,4 ) = v4 Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
16 / 65
Interpretaci´on de S´ısmica de Refracci´on: Varias Int Horizontales El ´angulo para la tercera refracci´ on critica es: v1 sin(θ1,3 ) = v3 v2 sin(θ2,3 ) = v3 Donde θ1,2 , θ2,3 y θ3,4 son los ´angulos criticos para la segunda, tercera y cuarta capa.
Las velocidades aparentes registradas a lo largo de la superifice de la tierraa a
(son inversos de las pendientes observadas en la curva de travel-time:
vap2 =
v1 = v2 , sin(θ1,2 )
vap3 =
v1 = v3 , sin(θ1,3 )
vap4 =
v1 = v4 sin(θ1,4 )
vap2 , vap3 y vap4 . Asociadas a las velocidades aparente de refracci´on de la segunda, tercer y cuarta capa Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
17 / 65
Interpretaci´on de S´ısmica de Refracci´on: Varias Int Horizontales
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
18 / 65
Interpretaci´on de S´ısmica de Refracci´on:Varias Int Horizontales
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
19 / 65
Generalizaci´on de la Refracci´on de un modelo de varias interfaces
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
20 / 65
Generalizaci´on de la Refracci´on de un modelo de varias interfaces
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
21 / 65
Generalizaci´on de la Refracci´on de un modelo de varias interfaces
Tiempo de Arribo para el Caso General tn = tn−1 +
x vn
Donde: tn , es el tiempo de propagaci´ on hacia abajo desde la fuente, a lo largo de la capas superior n tn−1 , es el intercepto con el eje t para la refracci´ on en la capa n x, es la distancia de la fuente al receptor y vn , la velocidad en la capa n
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
22 / 65
Generalizaci´on de la Refracci´on de un modelo de varias interfaces
tn−1 tn−1 =
n X 2hi cos(θi,n ) i=1
vn
Distancia Cr´ıtica xc xci =
n X
2hi tan(θi,n )
i=1
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
23 / 65
Refracci´on de la Corteza Oce´anica Velocidad: I I I I I
de arribo de la onda directa en el agua ≈ 1.5km/s Refracci´ on de la Capa 1 ≈ 2.0km/s Refracci´ on de la Capa 2 ≈ 5.07 ± 0.63km/s Refracci´ on de la Capa 3 ≈ 6.69 ± 0.26km/s Refracci´ on de la Capa 4 ≈ 8.13 ± 0.24km/s
Espesores I I I I
Agua ≈ 4.5 km Capa 1 ≈ 0.5 km Capa 2 ≈ 1.71 ± 0.75 km Capa 3 ≈ 4.86 ± 1.42 km
Secuencia tipica de Ophiolite (veremos observaciones de cuenca oce´anica) Ver la figura basada en la inversi´ on de algunos perfiles de refracci´on
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
24 / 65
Refracci´on de la Corteza Oce´anica
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
25 / 65
Refracci´on de la Corteza Oce´anica
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
26 / 65
Una Interface Inclinada Principio de Reciprocidad: El tiempo de propagaci´on de una onda s´ısmica del punto A al B, es el mismo que del punto B hacia A. Esto es una consecuencia del principio de Fermat, o del recorrido del tiempo m´ınimo El principio de Reciprocidad nos dice que : tAB = tBA
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
27 / 65
Una Interface Inclinada
v2d , es la velocidad aparente (hacia abajo de la interface inclinada) v2u , es la velocidad aparente (hacia arriba de la interface inclinada) La velocidad aparente: v1 v2d = sin(θc + α) v1 Donde: v2u = sin(θc − α)
Jairo Torres (Geo UIS)
v1 , velocidad de la capa. θc , ´angulo cr´ıtico α, inclinaci´on de la interface
M´ etodos Geof´ısicos
20 de febrero de 2012
28 / 65
Una Interface Inclinada
v2u simultaneamente:
−1
sin α= El ´angulo de inclinaci´on de la interface y el ´angulo cr´ıtico puede ser determinado resolviendo v2d y
Jairo Torres (Geo UIS)
sin
M´ etodos Geof´ısicos
v1 v2d
−1
v1 v2u
v1 v2u
− sin 2
−1
θc =
v1 v2d
+ sin
−1
2
20 de febrero de 2012
29 / 65
Una Interface Inclinada por la lectura de la velocidad sobre la capa (v1 ) como la inversi´on de la pendiente del arribo directo, y si resolvemos para la velocidad de refracci´ on de la capa (v2 ) de la ley de snell:
v1 = La velocidad verdadera puede ser determinada de la curva travel-time
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
pendiente v2 =
1 arribo
directo
v1 sin(θc )
20 de febrero de 2012
30 / 65
Una Interface Inclinada Los tiempos de propagaci´on considerando la distancia horizontal (x): x td = t1d + v2d x tu = t1u + v2u Donde: 2hd cos(θc ) t1d = v1 2hu cos(θc ) t1u = v1 hd , es la distancia perpendicular a la interface cuando se esta disparando hacia abajo v1 t1d hd = 2 cos θc hu , es la distancia perpendicular a la interface cuando se esta disparando hacia arriba v1 t1u hu = 2 cos θc Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
31 / 65
Una Interface Inclinada
Las profundidades verticales de la interface respecto a los puntos A y B son respectivamente
Jairo Torres (Geo UIS)
zd =
hd cos α
zu =
hu cos α
M´ etodos Geof´ısicos
20 de febrero de 2012
32 / 65
Problemas del M´etodo S´ısmico La se˜ nal que se recibe es muy peque˜ na para distinguirla del ruido. El m´etodo del CMP y el Apilado (stack) resuelven en gran medida este problema. Solo sobreviven las frecuencias bajas, lo que implica que no se pueden dsitinguir ecos de capas delgadas. Hasta cierto punto la deconvoluci´on ayuda. La se˜ nal rebota varias veces haciendo que lleguen varios ecos de la misma capa. T´ecnicas de filtrado en diversos espacios disminuyen este problema. Las capas pueden estar inlinadas, lo que dificulta saber a que profundidad se gener´o el eco. La correcci´ on DMO y la migraci´on ayuda dan soluci´on a esta dificultad.
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
33 / 65
El espectro de frecuencias s´ısmico
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
34 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento Recordando
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
35 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento a
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
36 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento a
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
37 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
38 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
39 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
40 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
41 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
42 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
43 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
44 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
45 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
46 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
47 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
48 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
49 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
50 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
51 / 65
Tipos de Velocidad
Velocidad Promedio Velocidad RMS Velocidad de Apilado Velocidad de Intervalo
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
52 / 65
Tipos de Velocidad Velocidad Promedio Es la distancia en la cual se ecuentra la interface z, dividido por el tiempo de propagaci´on en una sola direcci´on vav =
z t
Donde: z, es la profundidad de la capa. t, es el tiempo de propagaci´on de la onda (de ida) Para el caso de varias interfaces: PN
vav = Pi=1 N
vi t i
i=1 ti
Donde: vi , es la velocidad en la i -esima capa ti , es el tiempo de propagaci´on (ida-vuelta) en la i -esima capa Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
53 / 65
Tipos de Velocidad
Velocidad RMS: v u N u 1 X t vrms = vi2 ti (0) ttot (0) i=1
Velocidad de Apilado: Es la velocidad con la cual se corrigen mejor los eventos en el
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
CMP Gather para el Normal Moveout (NMO) Velocidad de Intervalo : Es la velocidad promedio del material entre dos interface: vint =
∆Z ∆t
20 de febrero de 2012
54 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
55 / 65
S´ısmica de Refracci´on: Adquisici´on y Procesamiento
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
56 / 65
Interpretaci´on Tectonica de S´ısmica de Reflexi´on Interface Inclinada
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
57 / 65
Interpretaci´on Tectonica de S´ısmica de Reflexi´on Secci´ on migrada en tiempo, Rayos Normales
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
58 / 65
Interpretaci´on Tectonica de S´ısmica de Reflexi´on Migraci´ on en tiempo a migraci´ on en profundidad
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
59 / 65
Interpretaci´on Tectonica de S´ısmica de Reflexi´on Eventos No migrados y migrados
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
60 / 65
Interpretaci´on Tectonica de S´ısmica de Reflexi´on Estructuras de Interfaces inclinadas
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
61 / 65
Interpretaci´on Tectonica de S´ısmica de Reflexi´on Expresi´ on de una punto difractor
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
62 / 65
Interpretaci´on Tectonica de S´ısmica de Reflexi´on Expresi´ on s´ısmica de un anticlinal
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
63 / 65
Interpretaci´on Tectonica de S´ısmica de Reflexi´on Expresi´ on s´ısmica de un sinclinal
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
64 / 65
Interpretaci´on Tectonica de S´ısmica de Reflexi´on ondulating or hummocky surface (l´ıneas solidas porciones de un sinclinal y l´ınea punteada porciones de un anitclinal
Jairo Torres (Geo UIS)
M´ etodos Geof´ısicos
20 de febrero de 2012
65 / 65