Taller geofísica - Refracción sísmica
Profesor: Gaspar Monsalve
Por:
Luisa Fernanda Huertas Marulanda CC1018512801
Manuela Herrera Ochoa CC 1035919014
Fecha de entrega: 23 de marzo
Universidad Nacional de Colombia
Sede Medellín
1. Se tiene un registro de sismogramas correspondiente a una unidad de refracción sísmica de
12 geófonos (imagen abajo). Asuma que los geófonos están separados entre sí 1.5 m, y que
existe esa misma separación entre el sitio de disparo y el primer geófono. La escala de tiempo
está dada en segundos.
(a) Recolecte los tiempos de llegada de la primera onda P y represéntelos en un gráfico t vs x
(b) ¿Cuántas capas se identifican?, Halle las velocidades y espesores de cada capa suponiendo una
configuración horizontal
t (s) x (m)
1.5 0.017 x (m) contra t (s)
3 0.021
0,05
4.5 0.025
6 0.028 0,04
7.5 0.029 0,03
x (m)
9 0.031
0,02
10.5 0.034
12 0.037 0,01
13.5 0.039 0
15 0.041 0 5 10 15 20
16.5 0.044 t (s)
18 0.047
En la gráfica podemos observar claramente que tenemos 2 capas, la segunda con un aumento en
su velocidad respecto a la primera, dado que podemos obtener los datos de pendiente e
intercepto por medio de la gráfica, esto lo podemos usar para hallar la velocidad y el espesor de
las capas.
• Tramo 1
Capa 1
y = 0,0025x + 0,0135
0,03
0,025
0,02
x (m)
0,015
0,01
0,005
0
0 1 2 3 4 5 6 7
t (s)
Dado que tenemos la ecuación de la recta podemos despejar la velocidad pues la pendiente es la
inverso de esta:
1 1
𝑚= 𝑉1 = 𝑉1 = 400 𝑚/𝑠
𝑉1 0,0025
• Tramo 2
Capa 2
y = 0,0017x + 0,0161
0,05
0,045
0,04
0,035
0,03
x (m)
0,025
0,02
0,015
0,01
0,005
0
7 9 11 13 15 17 19
t (s)
A continuación, despejamos a su vez la velocidad de la capa 2.
1 1
𝑚= 𝑉1 = 𝑉2 = 588,23𝑚/𝑠
𝑉1 0,0017
Dado que ya tenemos las velocidades de las capas podemos hallar los ángulos por medio de la ley
de Snell:
𝑉1 𝑉1
𝑆𝑒𝑛 (𝑖𝑐) = 𝑖𝑐 = 𝑎𝑟𝑐𝑠𝑒𝑛 ( ) 𝑖𝑐 = 42,84
𝑉2 𝑉2
Ahora por medio del valor del intercepto t0 podremos hallar los espesores de las capas de la
mediante la siguiente fórmula:
2 ∗ ℎ1 ∗ cos(𝑖𝑐)
𝑡0 =
𝑉1
• Espesor para la capa 1
t0 en la capa 1 tiene un valor de 0,0135, de este modo despejamos el espesor h1:
𝑡0 ∗ 𝑉1 0,0135 ∗ 400
ℎ1 = ℎ1 = ℎ1 = 3,68 𝑚
2 ∗ cos (𝑖𝑐) 2 ∗ cos (42,82)
• Espesor para la capa 2
t0 en la capa 2 tiene un valor de 0,0161, de este modo despejamos el espesor h1:
𝑡0 ∗ 𝑉1 0,0161 ∗ 588,23
ℎ2 = ℎ2 = ℎ2 = 6,455 𝑚
2 ∗ cos (𝑖𝑐) 2 ∗ cos (42,82)
Los datos recolectados en este estudio fueron:
Capa 1 Capa 2
Velocidad m/s 400 588,23
Espesor m 3,68 6,455
2.
3. Para la siguiente configuración del subsuelo elabore el gráfico de dromocrónicas (gráficos t
vs x en ambos sentidos) con disparos desde ambos extremos de la línea: Se tienen 3 medios
formando una estructura de capas paralelas, con velocidad constante, y contactos buzando
11° en la dirección del experimento. La longitud total del arreglo es de 200m. V1 = 750 m/s, V2
= 1400 m/s, V3 = 2100 m/s, h1 = 8 m, h2 = 12 m. ¿Hay capas ocultas en el registro?
Para poder realizar el gráfico dromocrónicas será necesario obtener los tiempos de onda de la
refracción directa, primera refracción y segunda refracción.
En principio realizaremos el viaje de izquierda a derecha, empezando así a calcular los tiempos de
la onda directa, primera refracción y segunda refracción, para esto tomaremos primeramente el
buzamiento a favor.
Los datos que podemos hallar de inmediato o que fueron suministrados son:
𝑉1 𝑉1 𝑉2
𝑖11 = 𝑖12 = 𝑖13 =
𝑉2 𝑉3 𝑉3
Incidencia Radianes Grados
V1 = 750 m/s L = 200 m
i11 0,56535 32,39
V2 = 1400 m/s h1 = 8 m
i12 0,36521 20,92
i22 0,66667 41,9 V3 = 2100 m/s h2 = 12m
Buzamiento a favor: Para el buzamiento a favor hallamos los datos de pendiente y de intercepto
para así entonces realizar el respectivo diagrama:
• Tramo 1
1 1
𝑃𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒 = =
𝑉1 750
𝑃𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒 = 0,0013
• Tramo 2
𝑠𝑒𝑛(𝑖11 + 𝜃) 𝑠𝑒𝑛(32 + 11) 2 ∗ ℎ1 ∗ cos (𝑖11 )
𝑃𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒 = = 𝐼𝑛𝑡𝑒𝑟𝑐𝑒𝑝𝑡𝑜 =
𝑉1 750 𝑉1
𝑃𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒 = 0,0009 𝐼𝑛𝑡𝑒𝑟𝑐𝑒𝑝𝑡𝑜 = 0,018
• Tramo 3
𝑠𝑒𝑛(𝑖12 + 𝜃) 𝑠𝑒𝑛(20,9 + 11)
𝑃𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒 = = = 0,0007
𝑉1 750
2 ∗ ℎ1 ∗ cos (20,9) 2 ∗ ℎ2 ∗ cos (41,9)
𝐼𝑛𝑡𝑒𝑟𝑐𝑒𝑝𝑡𝑜 = + = 0,034
𝑉1 𝑉2
Para el buzamiento en contra debemos encontrar en principio h´1, para hallarla usaremos la
siguiente ecuación:
ℎ´1 = ℎ1 + 𝐿 ∗ 𝑆𝑒𝑛𝜃 ℎ´1 = 8 + 200 ∗ 𝑆𝑒𝑛(11º) ℎ´1 = 46,16
Con esta información podemos calcular las pendientes e intercepto de las capas de la :
• Tramo 1
1 1
𝑃𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒 = =
𝑉1 750
𝑃𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒 = 0,0013
• Tramo 2.
𝑆𝑒𝑛(𝑖11 − 𝜃)
𝑃𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒 = = 0,00048
𝑉1
2 ∗ ℎ´1 ∗ cos (𝑖11 )
𝐼𝑛𝑡𝑒𝑟𝑐𝑒𝑝𝑡𝑜 = = 0,1039
𝑉1
• Tramo 3
𝑠𝑒𝑛 (𝑖12 − 𝜃)
𝑃𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒 = = 0,00022
𝑉1
2 ∗ ℎ1´ ∗ cos (𝑖12 ) 2 ∗ ℎ2 ∗ cos (𝑖22 )
𝐼𝑛𝑡𝑒𝑟𝑐𝑒𝑝𝑡𝑜 = + = 0,12775
𝑉1 𝑉2
Luego de recolectar todos los datos, las ecuaciones entonces a favor y en contra del buzamiento
serían:
A favor del buzamiento:
- Onda directa= 0,0013x
- Primera onda refractada:0,0009x+0,018
- Segunda onda refractada: 0,0007x+0,034
En contra del buzamiento
- Onda directa = 0,0013x
- Primera onda refractada: 0,00048x+0,1039
- Segunda onda refractada: 0,00022x+0,12775
En base a los datos recolectados es posible realizar el gráfico de dromocrónicas: Profe este no es el
gráfico, se me dañó el archivo de Excel y no lo pude recuperar a tiempo perdón
4.
5.
Se evidencia un contacto vertical en el punto rojo, gracias a que la pendiente tiene un aumento en
vez de una disminución, en los otros dos puntos se muestra una disminución de la pendiente la
cual se debe a que es una refracción crítica; en el contra perfil no se evidencia el contacto vertical
porque las velocidades van de menor a mayor. Puede ser que el perfil horizontal esté buzando un
poco gracias a que las líneas de la primera pendiente tienen una línea mas larga que la otra lo que
quiere decir que demoro más tiempo un lado del otro lo que puede inferir que el estrato
horizontal está buzando hacia la derecha
Distancia(m) T Perfil ms Distancia(m) T contra-p
0 0 0 51
10 5 10 49
20 10 20 47
30 14 30 45
40 17 40 43
50 20 50 41
60 23 60 39
70 26 70 38
80 29 80 35
90 32 90 30
100 35 100 25
110 38 110 20
120 41 120 15
130 44 130 10
140 47 140 5
150 51 150 0
A)
Dist Distancia (m) T Perfil ms Distancia(m ) T contra-p
Perfil 0 0 0 51
60 10 5 10 47
40
20 10 20 44
20
30 14 30 41
0
0 50 100 150 200
40 17 40 38
50 20 50 35
60 23 60 32
70 26 70 29
80 29 80 26
90 32 90 23
T contra-p 100 35 100 20
60
110 38 110 17
120 41 120 14
40
130 44 130 10
20
140 47 140 5
0
150 51 150 0
0 50 100 150 200
Perfil y contra perfil
60
40
20
0
0 50 100 150 200
B) Hallar la velocidad de la onda directa
Tomamos la primera pendiente y para hallar V1
V1=1/m1
V1=1/0,5
V1=2
C) Hallar las velocidades aparentes (Inversas de las pendientes) de la capa que esta debajo, en perfil y
contra perfil.
Vap2=1/m2 60
y = 0,3033x + 4,7802
40
20
0
0 50 100 150 200
Vap2=1/0,3033
Vap2= 3,29
Vap3=1/m3 60
50
40
y = -0,195x + 50,911
30
20
10
0
0 50 100
Vap3= 1/0,19
Vap3= 5.3
D)Estime el buzamiento de la frontera en la dirección del perfil m2=sen(i11+a)/V1
Viaje en el buzamiento a favor.
m3=sen(i11-a)/V1 Viaje en el buzamiento en contra.
Despejamos he igualamos las expresiones i11=30°
a=7.4°
E) Hallar la V2 por ley de Snell
sen(i11)=v1/v2
V2=V1/sen(i11) V2=4
