Fórmulas de Geofísica aplicada Prospección sísmica

PROSPECCIÓN SÍSMICA

Ecuaciones de sísmica de refracción

Las siguientes ecuaciones sólo son válidas si las velocidades aumentan hacia abajo.
x 2 ⋅ e ⋅ cos ic
Caso de dos capas horizontales (una interfase horizontal): t= + , donde t es
V2 V1
el tiempo, x la distancia, e el espesor de la capa superior, V1 y V2 la velocidad de las ondas P
en las capas superior e inferior respectivamente, e ic el ángulo crítico. Esa ecuación
relaciona distancias con tiempos de llegada en función del ángulo crítico. Las velocidades
son el inverso de las pendientes de las dromocronas. Se usa más en la forma siguiente, que
está sólo en función de las velocidades:
x 2 ⋅ e ⋅ V2 − V1 2 ⋅ e ⋅ V22 − V12
2 2

t= + . Para x = 0 , el tiempo de intersección es: t i =

V2 V1 ⋅ V2 V1 ⋅ V2
Se usan además las siguientes fórmulas para calcular el espesor (e) del lecho superior:
t V1 ⋅ V2 x V − V1
e= i ⋅ , en función del tiempo de intersección (ti), y e = cr ⋅ 2 , en función
2 V22 − V12 2 V2 + V1
de la distancia de cruce (xcr).
x 2 ⋅ e1 ⋅ V3 − V1 2 ⋅ e2 ⋅ V32 − V22
2 2

Caso de tres capas horizontales: t = + + , donde t es el

V3 V1 ⋅ V3 V2 ⋅ V3
tiempo, x la distancia, e1 y e2 los espesores de la capa superior e intermedia, y V1, V2 y V3 la
velocidad de las ondas P en las capas superior, media e inferior respectivamente. Las
velocidades son el inverso de las pendientes de las dromocronas.
Para x = 0 , el tiempo de intersección de la primera interfase es ti1, y se calcula por la fórmula para
el caso de dos capas. Para la segunda interfase, el tiempo de intersección es:
2 ⋅ e1 ⋅ V32 − V12 2 ⋅ e2 ⋅ V32 − V22
ti2 = +
V1 ⋅ V3 V2 ⋅ V3
El espesor e1 se calcula por la fórmula para dos capas. El espesor e2 es:
1 ⎛⎜ 2 ⋅ e1 ⋅ V32 − V12 ⎞⎟ ⎛⎜ V2 ⋅ V3 ⎞⎟
e2 = ⋅ t i 2 − ⋅
2 ⎜ V1 ⋅ V3 ⎟ ⎜ V 2 −V 2 ⎟
⎝ ⎠ ⎝ 3 2 ⎠

Caso de cuatro capas horizontales:

x 2 ⋅ e1 ⋅ V4 − V1
2 2
2 ⋅ e2 ⋅ V42 − V22 2 ⋅ e3 ⋅ V42 − V32
t= + + +
V4 V1 ⋅ V4 V2 ⋅ V4 V3 ⋅ V4
Las velocidades se obtienen de la forma habitual, invirtiendo las pendientes. Los espesores e1 y e2
se obtienen por las fórmulas de los casos anteriores. El tiempo de intersección para la tercera
2 ⋅ e1 ⋅ V42 − V12 2 ⋅ e 2 ⋅ V42 − V22 2 ⋅ e3 ⋅ V42 − V32
interfase es: ti3 = + + ,
V1 ⋅ V4 V2 ⋅ V4 V3 ⋅ V 4
y el espesor e3 se calcula despejando de esa fórmula:
1 ⎛ 2 ⋅ e1 ⋅ V42 − V12 2 ⋅ e2 ⋅ V42 − V22 ⎞⎟ ⎛⎜ V3 ⋅ V4 ⎞⎟
e3 = ⋅ ⎜ t i 3 − − ⋅
2 ⎜ V1 ⋅ V4 V2 ⋅ V4 ⎟ ⎜ V 2 −V 2 ⎟
⎝ ⎠ ⎝ 4 3 ⎠

Pueden darse ecuaciones para más capas, pero no tiene sentido porque 4 ya suelen ser muy
difíciles de diferenciar en un perfil sísmico de refracción de una exploración local.

9
Fórmulas de Geofísica aplicada Prospección sísmica

Caso de dos capas inclinadas: en este caso, se trabaja con dos disparos, uno en cada extremo del
dispositivo. Se definen dos pendientes (md y mu) y dos tiempos de intersección (tid y tiu) para
las dromocronas correspondientes a la interfase. Se usan los subíndices d y u para cuando
las ondas de Mintrop recorren la interfase hacia abajo (md, td) y hacia arriba (mu, tu).
Los tiempos hacia abajo y hacia arriba, en función de la distancia, son:
x 2 ⋅ e A ⋅ cos ic x 2 ⋅ e D ⋅ cos ic
t d = ⋅ sen (ic + β ) + y t u = ⋅ sen (ic − β ) + ,
V1 V1 V1 V1
donde x es la distancia, eA y eD los espesores de la capa superior pendiente arriba y
pendiente abajo respectivamente, V1 la velocidad de las ondas P en la capa superior, ic el
ángulo crítico y β el buzamiento de las capas.
sen (ic + β ) sen (ic − β )
Las pendientes para la onda Mintrop son: md = y mu =
V1 V1
2 ⋅ e A ⋅ cos ic 2 ⋅ e D ⋅ cos ic
Y los tiempos de intersección son: t i ⋅d = y t iu =
V1 V1
Esas ecuaciones se combinan para obtener las más empleadas:
ic = ⋅ (arc sen ( md ⋅ V1 ) + arc sen (mu ⋅ V1 ) )
1
2
β = ⋅ (arc sen (m d ⋅ V1 ) − arc sen (mu ⋅ V1 ) )
1
2
2 ⋅ cos β 2
La velocidad V2 es: V2 = y, si el buzamiento es pequeño, V2 ≅
( m d + mu ) ( m d + mu )
Salto de una falla (s): se parte de dos capas horizontales desplazadas por una falla vertical. Se
definen dos tiempos de intersección, para las dromocronas del labio levantado y hundido
V1 ⋅ V2
respectivamente (ti1 y ti2). El salto viene dado por: s = (t i 2 − t i1 ) ⋅
V22 − V12
Medio con gradiente continuo de velocidad: la velocidad aumenta con la profundidad según
VP = V0 + k ⋅ z , donde k es una constante y V0 la velocidad para z = 0 ; k vale entre 0,5 y 1
s-1 en capas sedimentarias y sus unidades son km/s cada kilómetro, es decir, km s-1/km = s-1.

Relación entre la distancia que recorre un rayo y la profundidad a la que penetra, en un medio con
gradiente continuo de velocidad:
2
V0 ⎞ ⎡ 1 − p 2 ⋅ V02 ⎤
2
1 ⎛
= ⎜ z + ⎟ + ⎢ x − ⎥ , donde V0 es la velocidad de las ondas P
k 2 ⋅ p2 ⎝ k ⎠ ⎢ k⋅p ⎥⎦
⎣
en la superficie, k la constante que regula el incremento de velocidad con la profundidad, x
la proyección sobre la horizontal de la distancia recorrida por el rayo, z la profundidad para
esa distancia y p el parámetro del rayo.
1
Es la ecuación de una circunferencia cuyo radio es , y cuyo centro está sobre una línea
k⋅p
V
horizontal a una distancia 0 por encima de la superficie. La distancia entre la fuente y la
k
1 − p 2 ⋅ V02
proyección sobre la horizontal del centro de la circunferencia viene dada por
k⋅p

10
Fórmulas de Geofísica aplicada Prospección sísmica

Correcciones en sísmica de refracción

Tiempo de retardo (ΔtR): es la diferencia entre el tiempo empleado por la onda refractada en llegar
a un punto separado una distancia x de la fuente, y el tiempo que hubiera empleado viajando
x
horizontalmente esa misma distancia por el lecho de velocidad mayor (V2): Δt R = t −
V2
El tiempo de retardo tiene dos componentes, una relacionada con el extremo del disparo
(ΔtR1)y otra con el del receptor (ΔtR2): Δt R = Δt R1 + Δt R 2
En el caso de una interfase horizontal, las dos componentes son iguales: Δt R1 = Δt R 2 , y por
2 ⋅ e ⋅ V22 − V12
tanto: Δt R = Δt R1 + Δt R 2 = 2 ⋅ Δt R1 =
V1 ⋅ V2
En el caso general, los espesores son distintos bajo el disparo (e1) y bajo el receptor (e2).
e1 ⋅ V22 − V12 e2 ⋅ V22 − V12
Entonces: Δt R1 = , y Δt R 2 =
V1 ⋅ V2 V1 ⋅ V2
Normalmente se realiza una aproximación, válida cuando los buzamientos no son muy
e ⋅ V22 − V12
grandes, y se supone que Δt R1 ≅ Δt R 2 = , donde e es el espesor medio.
V1 ⋅ V2
Corrección topográfica (ΔtT): se trata de "re-posicionar" tanto la fuente como el receptor en un
plano horizontal de referencia ("datum"). Es decir, la corrección topográfica es el tiempo
que hay que sumar o restar al tiempo total para obtener el tiempo que habría tardado la onda
refractada si la fuente y el receptor estuvieran a la misma cota:
(h1 − p + h2 − 2 ⋅ d ) ⋅ V22 − V12
Δt T = , donde h1 es la cota de la
V1 ⋅ V2
embocadura del sondeo desde el que se hace el disparo, p su profundidad, h2 la cota del
receptor, y d la cota del "datum". La corrección topográfica se resta al tiempo observado si
el "datum" está por debajo del disparo y del receptor. Eso es lo normal, es decir, el "datum"
se escoge por debajo. La resta equivale a bajar toda la dromocrona un tiempo igual a ΔtT.
Corrección de alteración (Δta): se trata también de "re-posicionar" tanto la fuente como el receptor
en un plano horizontal de referencia ("datum"), pero ahora hay que tener en cuenta la
existencia de un lecho alterado de baja velocidad (V0) y espesor a.
Si no hay diferencias de cota entre el disparo y el receptor, la corrección de alteración es:
a ⋅ V22 − V02
Δt a = , donde a es el espesor del lecho alterado por debajo del receptor.
V0 ⋅ V 2
Pero lo normal es que haya diferencias de cota y, además, el disparo se realiza en un sondeo,
que atraviesa el lecho alterado, de forma que la explosión se hace en la roca no alterada.
Entonces, la corrección de alteración es:
V22 − V12 a ⋅ V22 − V02
Δt a = (h1 − p + h2 − 2 ⋅ d − a ) ⋅ + , donde h1 es
V1 ⋅ V2 V0 ⋅ V 2
la cota de la embocadura del sondeo desde el que se hace el disparo, p su profundidad, h2 la
cota del receptor, d la cota del "datum" y a el espesor del lecho alterado por debajo del
receptor. Como normalmente el sondeo atraviesa todo el lecho alterado en el lado del
disparo, su espesor ahí no se tiene en cuenta. La corrección de alteración se resta al tiempo
observado si el "datum" está por debajo del disparo y del receptor, que es lo normal. La resta
equivale a bajar toda la dromocrona un tiempo igual a Δta.

11
Fórmulas de Geofísica aplicada Prospección sísmica

Exploración cortical y sísmica de gran ángulo

Tiempo reducido (tr): los perfiles de sísmica de gran ángulo, o DSS ("deep seismic sounding") se
x
suelen proyectar en un gráfico distancia-tiempo reducido (o normalizado): t r = t − ,
Vr
donde Vr es la velocidad de reducción. En corteza continental se emplea Vr = 6 km s-1 y
para investigar el manto, Vr = 8 km s-1. Eso para las ondas P. Cuando se usan ondas S, se
aplica a esas velocidades una reducción de un factor 1 3 . Una onda refractada que viaje a
velocidad igual a Vr, dará llegadas que se alinearán según una recta horizontal en el gráfico
distancia-tiempo reducido. Si su velocidad es menor que Vr, la recta tendrá su pendiente
inclinada hacia el disparo, y si es mayor, la pendiente se inclinará en sentido opuesto. Las
llegadas que se alinean según curvas son reflexiones de alto ángulo, o refracciones en un
lecho con un gradiente de velocidades.
x
Para calcular el tiempo real a partir del reducido, no hay más que despejar: t = tr +
Vr

Ecuaciones de las diferentes llegadas en una agrupación de disparo en sísmica de reflexión

x
Onda directa: t= , donde t es el tiempo, x la distancia y V1 la velocidad del primer lecho.
V1
x 2 ⋅ e ⋅ V2 − V1
2 2

Refracciones en interfases horizontales: t = + , donde e es el espesor del

V2 V1 ⋅ V2
primer lecho y V1 y V2 las velocidades encima y debajo de la interfase.
2
2 ⎛ x⎞
Reflexiones en interfases horizontales: t = ⋅ z2 + ⎜ ⎟ , donde z es la profundidad y V1 la
V1 ⎝2⎠
velocidad por encima de la interfase.

Resolución en sísmica de reflexión

Resolución vertical: para que las capas sean detectables en sísmica de reflexión, su espesor (e)
debe guardar cierta relación con la longitud de onda (λ) dominante. El límite de
detectabilidad está alrededor de e = λ 30 .
Cuando e = λ 4 , se le llama espesor de sintonización, y entonces, techo y muro de la capa
vienen señalados por un surco y un pico consecutivos (en pulsos cero-fase).
Resolución se define como el espesor para el cual pueden diferenciarse techo y muro de una
capa, y es cualquier espesor igual o superior al espesor de sintonización.
Resolución horizontal: viene determinada por la anchura de la zona de Fresnel, que es el área de
un reflector que produce una interferencia constructiva de las ondas reflejadas.
Este área es un círculo en un reflector horizontal, y su radio (rf) es el del círculo intersectado
por el frente de onda en la superficie horizontal, cuando la onda que va 1/4 de λ por detrás
de llega a esa superficie.
V t
Su valor es rf = ⋅ , donde V es la velocidad media, t el tiempo de ida y vuelta, o
2 f
TWTT ("two-way travel time"), y f la frecuencia dominante.

________________________________________________________________________________
José Ramón Martínez Catalán. Octubre de 2007

12