UNIVERSIDAD NACIONAL DE SALTA

REFRACCIÓN ESCUELA DE GEOLOGÍA

Tiempo de la Refracción CATEDRA de GEOFÍSICA

UNIDAD 4 : Sísmica de Refracción : Método de

refracción.- Caso de una o más capas horizontales.-
Trayectorias y relaciones tiempo-distancia.-
Dromocronas.- Determinación de espesores.- Capas
inclinadas.- velocidad aparente.- Dromocrona.- Perfil
y contraperfil.- Interpretación: fallas, escalones,
domos salinos, etc.

Prof. Néstor Vitulli

Geólogo
 UNIVERSIDAD NACIONAL DE SALTA
REFRACCIÓN ESCUELA DE GEOLOGÍA

Tiempo de la Refracción CATEDRA de GEOFÍSICA

En el estudio de la demostración de la ley de Snell, cuando la perturbación se realiza por el contacto entre V1–V2

sen i 
v1
c AB BC CD
v2 T ABCD    por la geometría de la figura
v1 v2 v1
x AB  CD
A x1
Dx x2 Dx D
2 AB BC
n n T ABCD 
v1

v2
i ic  i i12  i
i
z c ic ic ic
ic
z BC  x  2 z tg i
Dx Dx v1 AB 
cos i
B E F C v2
r=90o 2z x  2 z tg i v1 tg i  sen i
Tx   sen i 
x cos i . v1 v2 v2 cos i

Tx 
x

2z

2 z sen i v 2
1  x  2 z  2 z sen i
Tx
x

2z

. 1  sen2 i  
x 2 z cos i

Tx 
x 2 z.cos i

v2 cos i v
1
cos i v v
2 1 v 2 v1
cos i cos i v
1 v2 cos i v
1 v 2 v1
v2 v1

1
Esta expresión en un plano T - x corresponde a la ecuación de una recta con pendiente v
2
y ordenada al origen (término independiente) 2 z .cos i
x
v1 
T d v1
Si la perturbación viaja de A a D sin mediar refracción, lo hará con V1 y el tiempo que tarda será:
1 Prof. Néstor Vitulli
y está es la ecuación de una recta con pendiente v 1 Geólogo
 UNIVERSIDAD NACIONAL DE SALTA
REFRACCIÓN ESCUELA DE GEOLOGÍA

Dromocrona CATEDRA de GEOFÍSICA

xi xi 2 z.cos i
Ti   
v1 v 2 v1
también podemos despejar z
 1 1  v1 xi  v1  1
z  xi .    .  . 1   .
 v 1 v 2  2 cos i 2  v 2  cos i

si cos i  1  sen2 i
 v1 
xi v 2 . v 2  v1 xi
1
1 
xi  v2 v 2  v1
z .  .  . 
2 2
2 1 2  2 2 v 2  v1. v 2  v1
 v1 
1   v2 2 1
v v
 v2
v2  v12 v 2  v1. v 2  v1 xi v 2  v1 xi v 2  v 1
xi xi
z .  . z .
2 v 2  v1v 2  v1 2 v 2  v1. v 2  v1 2 v 2  v1 2 v 2  v 1

xi se denomina distancia de quiebre o intersección de las dromocronas de las ondas directa y refractada
La denominación de dromocrona se aplica a las gráficas que relacionan tiempo de recorrido con la distancia en
superficie entre detector y fuente. Cuando se hable de dromocronas verticales estas expresan la relación de los tiempos
de recorrido con la profundidad (distancia vertical de la fuente al detector ).
Prof. Néstor Vitulli
Geólogo
 UNIVERSIDAD NACIONAL DE SALTA
REFRACCION ESCUELA DE GEOLOGÍA

Dos capas paralelas a sup CATEDRA de GEOFÍSICA

Trayectorias de rayos en refracción
(caso de dos capas horizontales)
v3> v2> v1 v1  v 2  v 3
x
A PE x1 Dx x2
R Analizaremos la trayectoria para la perturbación
n n que transcurre por la interfase v 2   v 3
z1 a’ v1
a y arriba a un detector a distancia x de la fuente.
Sabemos que: sen a v 1 v1
z2
B
b
n E
v2 
sen g v 3
g  90 sen a 
b v3

además sen a v1
C g r=90o D v3 
sen b v 2
i1-3 = a i2-3 = b
El tiempo de recorrido de A a R, pasando por B, C, D, E, será:

AB BC CD DE ER z1 z2 CD  x  2 z1 tg a  2 z 2 tg b
TABCDER      AB  ER 
cos a
BC  DE 
cos b
y
v1 v 2 v3 v 2 v1
2 z1 2z 2 x  2z1 tg i 13  2z 2 tg i 23
llamamos a  i13 b  i 23 a ' i12 TA  R    
v1 cos i 13 v1 cos i 23 v3
x 2 z1 2 z1tgi13 2z2 2 z 2tgi23 sen i 13 tg i 23  sen i 23
reordenando TABCDER      tg i 13  cos i 23
v3 v1 cos i13 v3 v 2 cos i 23 v3 cos i 13

x 2 z1  sen i 13  2 z 2  sen i 23 
TA  R   1 v1  1 v 2 Prof. Néstor Vitulli
Remplazando y factor común v 3 v 1 cos i 13  v 3  v 2 cos i 23  v3  Geólogo
 UNIVERSIDAD NACIONAL DE SALTA
REFRACCION ESCUELA DE GEOLOGÍA

Dos capas paralelas a sup CATEDRA de GEOFÍSICA

,

sen i 12 v 1 sen i 13 v 1 v2 v2 v2 v1
 y 
sen i 23 v 2
 sen i 23 
sen i 13 y entonces sen i 13   sen i 13 
1 v2 v1 v1 v 3 v3
y si tuviera más capas (n) podría escribirse una generalidad
v1
sen i 1n 
vn
Retomando formula del tiempo
v1 v2
TA  R ,  sen i 13 y  sen i 23 Tx  TA  R 
x

2 z1
 
1  sen 2 i 13 
2z 2
 
1  sen 2 i 23 
v3 v3 v 3 v1 cos i 13 v 2 cos i 23

Para dos capas x 2 z1 cos i 13 2 z 2 cos i 23 x 2 z.cos i

TA  R    Recordando una Tx  
v3 v1 v2 capa v2 v1
Determinación de espesores.-

Si hacemos x=0 como z1 se puede calcular con los arribos que se propagan por la interface v1  v2
se puede despejar z2, escribiendo entonces:
 2 z 1 cos i 13  v 2
2 z1 2z 2 z 2   T 3  
T 3  cos i 13  cos i 23  v1  2 cos i 23
v1 v2  z1 
El valor de z1 era T 2 v1 v1 v v2 z 2   T 3   k 23
z1   T 2. k 12 haciendo k 12  k 13 
1
k 23   k 13 
2 cos i 12 2 cos i 12 2 cos i 13 2 cos i 23
Prof. Néstor Vitulli
Geólogo
 UNIVERSIDAD NACIONAL DE SALTA
REFRACCION ESCUELA DE GEOLOGÍA

Dos capas paralelas a sup CATEDRA de GEOFÍSICA

T(seg) 1/v3
Tx
Tx

T2 v1
1/v2 z1=
2 cos i1-2
T3 k1-2
z1 .
z2 = (T3 - k
T2 k1-3 ) 2-3
1/v1 v2
v1
k1-3= k2-3=
2 cos i1-3 2 cos i2-3
O xi x (m)

esta es una expresión que aparece en las planillas de cálculos en el campo, para varias capas horizontales.
Por ejemplo para 3 (tres) capas la formula sería:

z1 
T 2 v1
 T 2. k 12  z1   z1 z2 
z 2   T 3   k 23 z 3   T 4    k 34
2 cos i 12  k 13   k 14 k 24 

Prof. Néstor Vitulli

Geólogo
 UNIVERSIDAD NACIONAL DE SALTA
REFRACCION ESCUELA DE GEOLOGÍA

Planilla cálculo de refracción sencilla CATEDRA de GEOFÍSICA

Prof. Néstor Vitulli

Geólogo
 UNIVERSIDAD NACIONAL DE SALTA
REFRACCION ESCUELA DE GEOLOGÍA

Capa inclinada CATEDRA de GEOFÍSICA

Trayectorias de refracción, una capa inclinada
Caso de una capaperfilinclinada. v1
y contraperfil
sen a  a  i12
x
A R
v 2 AB BC CR
 TABCR   
n  v1 v2 v1
D
z AB z
n E i 12 = ac  CR Z
i 12
v1 cos i 12v1 
Z AE = x cos  v 1 cos i 12v 1
i 12
G 
RE = x sen  BC x cos  z tg i 12  Z tg i 12
B 
v1 v2 v2
C
F
2zcos i 12 x v2
Tbx = + v1
sen (i 12) z x cos  z tg i 12  Z tg i 12 Z
v1
Txb  TA  R   
2Zcos i 12 x
cos i 12v1 v2 cos i 12v1
Tsx = + sen (i 12)
v1 v1
z z tg i 12 Z Z tg i 12 x cos
sen i 12 v Txb  TA  R     
tg i 12  sen i 12 
1
cos i 12 v 1 v 2 cos i 12 v 1 v 2 v2
cos i 12 v2
Reemplazado tg por sen/cos y multiplicando y dividiendo por V1 en 2do y 4to término nos queda:
x cos
Txb 
z
cos i 12v1

1  sen 2 i 12 
Z
cos i 12v1
 
1  sen 2 i 12 
v2

z cos i 12 Z cos i 12 x cos Zz
Txb 
v1

v1

v2
Reordeno términos y Z  z  x sen  Ya que sen  
x
reemplazo a Z por

x sen cos i 2 z cos i 12 x  v 1  v1

Tbx 
x cos  V1

z cos i
12 
z cos i
12  12 T bx    cos   sen  cos i 12  sen i 12
v 1 v 1  v2  v 2
v2 V1 v1 v1 v1
2 z cos i 12 x
Tbx   sen i 12    Prof. Néstor Vitulli
Geólogo
v1 v1
 UNIVERSIDAD NACIONAL DE SALTA
REFRACCION ESCUELA DE GEOLOGÍA

Capas inclinada CATEDRA de GEOFÍSICA

Explotamos en R hacia A
Z cos i 12 x cos z cos i 12
Tsx    reemplazamos Z= Z  x sen 
v1 v2 v1
2 Z cos i 12 x sen  cos i 12 x cos  sen i 12 2 Z cos i 12 x
Tsx      sen  cos i 12  cos  sen i 12 
v1 v1 v1 v1 v1
2 Z cos i 12 x 2 z cos i 12 x
Tsx   sen i 12    y Tbx   sen i 12   
v1 v1 v1 v1
Si derivamos ambas expresiones
con respecto a x, tendremos:
Tsx 1 sen i 12    Tbx 1 sen i 12   
   
dx vaps v1 dx vapb v1
de aquí podemos escribir
v1
 sen i 12   ;  arcsen
v1
 i 12  
vapb vapb
v1
 sen i 12   ;  arcsen
v1
 i 12  

vaps vaps
sumando y restando las últimas expresiones obtenemos los valores de los ángulos i12
1 v1 v1  1 v1 v1 
i 12   arcsen  arcsen     arcsen  arcsen 
2 vapb vaps  2 vapb vaps 
Los valores de v1 y vapb y vaps los obtenemos de los registros o gráficas dromocrónicas. Así conocido i 12 puedo calcular v2 de la
expresión: v1
sen i 12  v2  v1 / sen i 12
v2 conocido v2 ó directamente i12 y calcular z y Z. Prof. Néstor Vitulli
Geólogo
 UNIVERSIDAD NACIONAL DE SALTA
REFRACCION ESCUELA DE GEOLOGÍA

Primeros arribos CATEDRA de GEOFÍSICA

Tiempos
(segundos)

Prof. Néstor Vitulli

Geólogo
 UNIVERSIDAD NACIONAL DE SALTA
REFRACCION ESCUELA DE GEOLOGÍA

Interpretación CATEDRA de GEOFÍSICA

Recordar
Limite de viajes de frentes de ondas directas, refracción de segunda y tercera capa
y la gráfica dromocrónica correspondiente.
1
v3

1
v2

1
v1
3 km 6 km

Ds
s 1 2 Ds1 3 4
seg
v1

Ds2 v2

Ds3 v3
Prof. Néstor Vitulli
Curva de contacto de frentes Geólogo
 UNIVERSIDAD NACIONAL DE SALTA
REFRACCION ESCUELA DE GEOLOGÍA

Interpretación CATEDRA de GEOFÍSICA

tiempo (seg)
Caso de contactos subverticales cubiertos. discordancias

Al registrar hacia la derecha en las condiciones del esquema la

dromocrona que se obtiene presenta la forma normal de una
0,8 -
secuencia de capas horizontales. La registración del contra-perfil
correspondiente producirá una dromocrona bien diferente y
permitirá estimar la anomalía. 1 2
1000 m/s distancia (km)

2000 m/s 4000 m/s

Refracción a través de un domo y de un pliegue.

-

tiempo (seg)
2

Los esquemas ilustran los

casos de perfiles en un
tiempo (seg)

/s)
-
2 sentido a través de un domo 370
0 (m

1
v2 y de una estructura plegada, 1 -
la respuesta es similar con 4900
m/s
1 - 1
v2 4500
m/s
la presencia de adelantos m/
s)
0(
típicos en las dromocronas 370
2 4
1
v1 2 4
distancia (km) en las distancias corres- distancia (km)
v1 (m/s) pondientes a la anomalía en 3700 (m/s)

v2 (m/s) el subsuelo. 4900 m/s

4500 m/s

Prof. Néstor Vitulli

Geólogo
 UNIVERSIDAD NACIONAL DE SALTA
REFRACCION ESCUELA DE GEOLOGÍA

Interpretación CATEDRA de GEOFÍSICA

Casos de un escalón con rechazo moderado.

1
v2 1
v2
Los gráficos de las figuras corres-

tiempo (seg)
ponden a las dromocronas y es- Dt
Dt
quema de subsuelo para fuente 1 1
sobre la parte baja y alta del esca- 6- v2 v2
lón de poco rechazo. Se dan las 4-
expresiones para la estimación 1
v1
del rechazo 2-

d2-d1 = Dt V1 cos i
20k m 20k m

-
v1 d1 v1
d2 d1 d2
v2 v2
- .
d2 d1 = D tv 1 cos i

Dromocronas en presencia de escalones (fallas) de gran rechazo.

tiempo (seg)
En las figuras se muestran las dromocronas 1

/s
v2

m
y esquemas de subsuelo correspondientes a

0
0
17
1
escalones de gran rechazo para registración v1 1
con la fuente sobre la parte alta del escalón y 0,6-
v2

/s
m
00
x2
con la fuente sobre la parte baja. Se aprecia
17
e
la necesidad de registrar perfiles opuestos y x2
se da la expresión para calcular la distancia x1 x1
10 20km 10 20km
del escalón a la fuente en la primera 3km i
situación.
d
/s
1700 m v1 v2

e = x2- d tg i 5000 m/s 5000 m/s

Prof. Néstor Vitulli

Geólogo
 UNIVERSIDAD NACIONAL DE SALTA
REFRACCION ESCUELA DE GEOLOGÍA

Interpretación CATEDRA de GEOFÍSICA

Diagrama que muestra las posiciones de las pendientes (inversas de las velocida-
des) para capa horizontal (paralela a la superficie) y capas buzamiento arriba y
abajo.

El tiempo de interceptación es
1 1 aproximadamente el mismo pero
tiempo

v1 v2d
1 necesitamos la registración en
v2h
1 contraperfil para poder determinar
T
v2u
el ángulo de buzamiento y el ángulo
crítico i, y de aquí la pendiente
distancia
correspondiente a la velocidad real
xu xd del medio subyacente a la interface.
v1 v2u

v2h

v2 v2d

Prof. Néstor Vitulli

Geólogo
 UNIVERSIDAD NACIONAL DE SALTA
REFRACCION ESCUELA DE GEOLOGÍA

Interpretación fallas típicas CATEDRA de GEOFÍSICA

Interferencia de las dromocrona en perfil y contraperfil por un escalón de falla.

TTOTAL TTOTAL

rayos emergentes
S1 S2
d1
d2

frentes de ondas Prof. Néstor Vitulli

afectados Geólogo
 UNIVERSIDAD NACIONAL DE SALTA
Próxima Unidad REFLEXIÓN ESCUELA DE GEOLOGÍA

CATEDRA de GEOFÍSICA

FIN

Prof. Néstor Vitulli

Geólogo