PROSPECCIÓN GEOFÍSICA

PROSPECCIÓN SÍSMICA

TÉCNICAS DE REFRACCIÓN

Jesús Díaz Curiel

 ÍNDICE

™ INTRODUCCIÓN
– Velocidades de propagación sísmica
– Relación de la velocidad con módulos dinámicos
– Fenómenos de propagación en una interfase
– Instrumentación sísmica
– Fuentes sísmicas
– Tipos de geófonos y su disposición
– Técnicas sísmicas
™ SÍSMICA DE REFRACCIÓN
– Bases de la técnica de refracción
– Estratificación plana no horizontal
– Aplicaciones de la sísmica de refracción
– Resolución y limitaciones
– Realización de perfiles de refracción
– Instrumentación de refracción
– Metodología de trabajo en campo
– Interpretación de perfiles de refracción
™ SÍSMICA EN SONDEO
– Up-hole
– Cross-hole
– Tomografía sísmica
 INTRODUCCIÓN

™ Base - El método sísmico de prospección se basa en la distinta

velocidad de propagación de las ondas sísmicas (elásticas) a
través de determinados medios y la correspondiente variación
de trayectoria en las superficies de separación entre los
medios de distinta velocidad.

™ Proceso genérico - Mediante la generación de una onda

sísmica y el estudio de los tiempos de llegada de las ondas
resultantes a un cierto número de puntos de medida, se
pueden deducir las posiciones y disposición de esas
interfases y la velocidad de propagación en cada uno de los
medios de distinta velocidad.
 ONDAS SÍSMICAS

™ CONCEPTO: Onda sísmica es la propagación de una deformación a

través de un medio elástico
™ TIPOS DE ONDAS
¾ Volumétricas
‰ Longitudinales (Compresionales)
‰ Transversales (cizalla) - medios sólidos

¾ Superficiales (Love y Rayleigh)

 VELOCIDADES DE PROPAGACIÓN SÍSMICA

™ FACTORES DE DEPENDENCIA:
Porosidad y fluido que rellena los poros. Cuando el
fluido que rellena los poros es líquido se tiene que:

1 = 1-φ + φ
Vr Vm Vf
– Grado de alteración (reduce la velocidad)
– Anisotropía: V varía según la dirección
– Compactación y presión de confinamiento
– Contenido en arcillas para medios sedimentarios
 VELOCIDADES DE PROPAGACIÓN DE ONDAS P
(m/s)

VARIACIÓN MEDIA VARIACIÓN MEDIA

SEDIMENTOS 300 – 1200 750 SAL GEMA 4000 – 6000 5000
ROCAS
1200 – 4000 2500 YESO 3000 – 4000 3500
SEDIMENTARIAS
ROCAS
4500 – 6800 5500 ANHIDRITA 4500 – 6500 5500
METAMÓRFICAS
ROCAS VOLCÁNICAS 3000 – 6500 4500 GRANITO 5600 – 6400 6000
ROCAS PLUTÓNICAS 5000 – 7500 6000 PIZARRAS 2700 – 4800 3750
SUELOS 100 – 500 300 CALIZAS 2000 – 5000 3500
ARENAS SECAS 100 – 600 350 DOLOMÍAS 3500 – 5500 4000
ARENAS HÚMEDAS 200 – 800 500 METANO 400 – 600 500
LIMOS 300 – 900 600 PETRÓLEO 900 – 1400 1200
MORRENAS 1500 – 2700 2000 LIGNITO 1000 – 2200 2000
ARENISCAS 1500 – 4000 2500 HULLA 1800 – 3000 2500
ARCILLAS 1200 – 2500 2000 ANTRACITA 2500 – 3400 3000
MARGAS 2500 – 3500 3000 AIRE 310 – 360 340
ROCAS ALTERADAS 1900 – 4000 3000 AGUA 1350 – 1650 1500
ACERO 5300 – 6000 HIELO 3100 – 4200 3600
 VELOCIDAD Y MÓDULOS DINÁMICOS

K= P
™ MODULO DE ∆V
COMPRESION 4
K = δ·VP2 - δ·VS2
3
deformació n transversal unitaria
Coeficiente de Poisson µ =
deformació n longitudinal unitaria

™ MÓDULO DE
POISSON
 VELOCIDAD Y MÓDULOS DINÁMICOS

Fuerza por unidad de área .

™ MÓDULO DE YOUNG, E = –––––––––––––––––––––––––––––––––––––
Cambio de longitud por unidad de longitud

F
F
™ MÓDULO DE E= A
YOUNG
L ∆L
L L
ED = 2·δ· VS2 (1 + µ)

F
F
G= A
™ MÓDULO DE
∆x
RIGIDEZ L L
G = δ· VS2
 PROPAGACIÓN DE ONDAS SÍSMICAS

™ EN LA PROPAGACIÓN DE UNA ONDA A TRAVÉS DE LAS

SUCESIVAS INTERFASES SE PRODUCEN DISTINTOS CAMBIOS DE
TRAYERCORIA SIMULTÁNEAMENTE
FENÓMENOS DE PROPAGACIÓN EN UNA INTERFASE

™ INTERFASE - SUPERFICIE DE SEPARACIÓN ENTRE MEDIOS DE

DISTINTA VELOCIDAD
™ FRENTES DE ONDA - PRINCIPIO DE HUYGENS:
™ RAYO (PERPENDICULAR AL FRENTE DE ONDAS)
™ TRANSMISIÓN Y REFLEXIÓN SIMULTÁNEAMENTE

Rayo Rayo Angulo de

Incidente Reflejado Rayo Incidencia
Incidente Critico Rayo
i c Reflejado
r
V1
V1

Rayo V2 Rayo
V2 Refractado Refractado
90`
t Critico

™ LEY DE SNELL
sen α2 = v 2
sen α1 v1
 INSTRUMENTACIÓN SÍSMICA

™ GEÓFONOS y CABLES

™ FUENTE DE ONDAS ELÁSTICAS

™ DETECTOR DEL MOMENTO DEL

IMPACTO (DISPARADOR)

™ SISMÓGRAFO
 TIPOS DE GEÓFONOS Y SU DISPOSICIÓN

Las ondas sísmicas producen

desplazamientos en dirección
arbitraria los geófonos deben
posicionarse para registrar la
mayor componente que
produzcan (vertical u horizontal).

™GEÓFONOS UNIDIRECCIONALES
¾ Miden movimientos horizontales o los verticales.
¾ Frecuencia propia (5 a 50 Hz).

™DISPOSICIONES TRIAXIALES
¾ Tres geófonos colocados perpendicularmente entre sí.
 TIPOS DE GEÓFONOS Y SU DISPOSICIÓN

Las ondas sísmicas producen desplazamientos en dirección arbitraria ?

los geófonos deben posicionarse para registrar la mayor componente que
produzcan (vertical u horizontal).

™ GEÓFONOS UNIDIRECCIONALES
– Miden movimientos horizontales o los verticales.
– Frecuencia propia (5 a 50 Hz).
™ DISPOSICIONES TRIAXIALES
– Tres geófonos colocados perpendicularmente entre sí.

TIPO DE ESTUDIO
Onda de estudio Superficie Pared horizontal Pared vertical
p directa horiz / long horiz / long vertical
p refracción crítica vertical horiz / trans horiz / normal
s directa y refractada horiz / trans vertical horiz / paralelo
TIPO Y DISPOSICIÓN DE GEÓFONO
ASPECTO DE TRAZAS SÍSMICAS

TIEMPOS CRECIENTES HACIA ABAJO

 FUENTES DE ENERGÍA SÍSMICA

™ EXPLOSIVAS
• GOMAS (DINAMITA)
• GELES
• MEZCLA DE GASES
™ IMPLOSIVAS
• VACÍO ENTRE PLACAS
™ SECUENCIALES DE PERCUSIÓN
• CAÍDA LIBRE DE PESO
• CAÍDA IMPULSADA (MUELLES, AIRE COMPRIMIDO, ...)
™ SECUENCIALES ESPECIALES
• CARTUCHOS
• ARCO ELECTRICO (SPARKER)
™ MARINAS (ACUÁTICAS)
• LIBERACIÓN DE AIRE COMPRIMIDO (AIRGUN)
• EMISIÓN DE AGUA A PRESIÓN
• VAPOR DE AGUA (STEAMGUN, VAPORCHOC)
™ VIBRATORIAS
• MONOFRECUANCIA (VIBROSEIS, MINISEIS)
• MULTIFRECUENCIA
 TIPOS DE DISPARADORES

™ DETECCIÓN DE TIEMPO CERO

– Relés de contacto (mecánicos

y electrónicos) para fuentes IEEE - 488 RS 232 C

impulsivas
VIDEO SYSTEM TRIGGER

GEOPH SWITCH

– Bobinas intermedias en los

cables de los explosores POWER 5A PRINTER

™ Geófonos de tiempo cero

(trigger) próximos al punto de
golpeo.
GND

– Ajuste de sensibilidad.
– Errores en apilamiento
 TÉCNICAS DE SÍSMICA DE REFRACCIÓN

SI EN UNA INTERFASE ENTRE DOS MEDIOS

V2 > V1 Ö sen α2 > sen α1 Ö α2 > α1
SI CONSIDERAMOS UN FRENTE DE ONDAS ESFÉRICO, UNO DE LOS
RAYOS INCIDIRÁ CON UN ÁNGULO αC TAL QUE:
V1
α2 = 90° → sen α2 = 1 → sen α1 = → α1 = Angulo crítico
V2
0 XK

PT G1 G2 G3 G4 G5

Angulo de
Rayo Incidencia
Incidente Critico
c
Z

V
1 r1 r1 r1 r1

V 90 o Rayo
2 Refractado
Critico
 REPRESENTACIÓN DE LA DROMOCRONA

TIEMPOS CRECIENTES HACIA ARRIBA

26 26

24 DROMOCRONA 24

22 22

20 20
18

Tiempo (ms)
Tiempo (ms)

18
16 16

14 14
12 12

10 10

8 8
6 6

4 4

2 2

2 6 10 14 18 22 26 30 34 38 42 46
Distancia Geofonos (m)
 FORMULAS DE REFRACCIÓN HORIZONTAL

X(i)
™ PARA LA ONDA DIRECTA: tD (i) = V
1

™ PARA LA ONDA REFRACTADA:

Λ λ(i) X(i) 2 · Z · V 22 − V 12 X(i)

tR (i) = + = ... = + = + t0
V1 V 2 V2 V 2 − V1 V2
™ LA DISTANCIA A PARTIR DE LA CUÁL LLEGAN ANTES LOS RAYOS
QUE SE PROPAGAN POR EL SEGUNDO MEDIO (DISTANCIA CRÍTICA)
SERÁ:
2 · Z · V 2 + V1
XC =
V 2 − V1
 CALCULOS PREVIOS

™ En una zona en cuyo subsuelo se encuentra una capa

superficial de litología sedimentaria sobre un medio volcánico,
se pretende determinar el espesor y la velocidad de dicho
recubrimiento mediante sísmica de refracción. Suponiendo
que el espesor de la capa superficial es mayor de 4 metros,
determinar la distancia máxima entre geófonos para que la
onda directa llegue antes que la onda refractada crítica en los
3 primeros geófonos.

™ Considerando la profundidad del primer medio en el subsuelo

de 2 metros, determinar el valor de la distancia crítica para los
casos de V2=1.5·V1, V2=2.0·V1 y V2=3.0·V1. Analizar los
resultados y recomendar un valor de la longitud del perfil.
 ESTRATIFICACIÓN PLANA NO HORIZONTAL

™ VELOCIDAD APARENTE del refractor (inverso de la pendiente de la

recta de llegadas de refracción) - depende de la velocidad real y del
buzamiento de la capa.
sen α12
V2D = V2
sen (β2 + α12)
= 2 · V 2A · V2B · cos β
V2
V 2A + V2B 2

™ METODOLOGÍA: Realización de registros situando la fuente en los

extremos del conjunto de geófonos y determinación de los tiempos
de llegada recíprocos.
™ La velocidad aparente de un refractor que muestra la dromocrona
para cada punto de tiro es igual al cociente entre la velocidad
superficial y el seno del ángulo que forma la recta de ajuste de dicho
refractor (denominado ángulo de emergencia).
∆x ∆s
V1AP = = = V0 > V0
∆t sen θ0 · ∆ t sen θ0
 ESTRATIFICACIÓN PLANA NO HORIZONTAL

DROMOCRONA

26 PTC 26
tAB = tBA
24 24
22 22
20 20
tAC = tCA
18 tCB = tBC 18

Tiempo (ms)
Tiempo (ms)

16 16
14 14
12 R1A 12
10 R2C R1C 10
R2B
8 8
6 6
tINT = tINT
4 4
2 2

2 6 10 14 18 22 26 30 34 38 42 46
PTA PTB
Distancia Geofonos (m)
 PRINCIPIOS DE LAS DROMOCRONAS

™ LEY DE LAS VELOCIDADES APARENTES

Las velocidades de propagación aparentes obtenidas para los tiros a
ambos lados del dispositivo son distintas, una menor y otra mayor
que la velocidad real del refractor.
™ PRINCIPIO DE RECIPROCIDAD
El tiempo de propagación de la onda sísmica desde un punto de tiro A
hasta un geófono B es el mismo que el que se obtiene situando el
punto de tiro en B y registrando con un geófono en A.
™ PRINCIPIO DE LOS TIEMPOS INTERCEPTADOS
Si en un mismo punto de la superficie se registran las dromocronas
en dispositivos extendidos es sentidos opuestos, los tiempos
interceptados por las rectas de ajuste, correspondientes a un mismo
refractor, son iguales.
™ PRINCIPIO DE PARALELISMO
Para refractores planos, las rectas correspondientes a un mismo
refractor en dromocronas obtenidas en el mismo perfil desde
distintos puntos de tiro tienen la misma pendiente (en la figura
R1A¦R1C y R2B¦R2C).
 RESOLUCIÓN Y LIMITACIONES

CAPAS OCULTAS - Espesor reducido CAPAS CIEGAS - Baja velocidad

INFLUENCIA TOPOGRÁFICA
 APLICACIONES

RIPABILIDAD DE DISTINTOS MEDIOS EN FUNCION DE LA VELOCIDAD

Excavadora tipo D-8-R (Caterpilar)
VELOCIDAD (x1000 m/s)
MEDIO LITOLOGICO
0.2 0.6 1.0 1.4 1.8 2.2 2.6 3.0 3.4 3.8 4.2 4.6
SUELOS
ARCILLAS
MORRENAS
ROCAS SEDIMENTARIAS
PIZARRA
ARENISCA
LIMOLITA
ARGILITA
CONGLOMERADO
BRECHA
CALICHE
CALIZA
ROCAS METAMORFICAS
ESQUISTOS
FILITAS
ROCAS IGNEAS
GRANITO
BASALTO
DOBRITA
MINERALES Y MENAS
CARBON
PIRITA
RIPABLE MARGINAL NO RIPABLE
 METODOLOGÍA DE CAMPO

™ ELECCIÓN DEL DISPOSITIVO

48: T ↓ ↓ ↓ ↓ ↓ ↓ T ↓ ↓ ↓ ↓ ↓ ↓ T ↓ ↓ ↓ ↓ ↓ ↓ T ↓ ↓ ↓ ↓ ↓ ↓ T
1 2 2 2 2 2 1 1 2 2 2 2 2 1 1 2 2 2 2 2 1 1 2 2 2 2 2 1

™ EXTENDIMIENTO DEL PERFIL

™ PROGRAMACIÓN DEL SISMÓGRAFO

™ REGISTRO DEL PERFIL

™ APROBACIÓN DE LAS TRAZAS

 SISMÓGRAFOS

Equipos electrónicos que realizan la detección del tiempo cero, la

recepción de señales eléctricas desde los geófonos, conversión A/D,
amplificación y almacenamiento en memoria.

™ CONVERSIÓN A/D
– Transformación de señal analógica en digital (100 µs)
– Multiplexado y desmultiplexado

™ AMPLIFICACIÓN
– Ganancia diferenciada en funcion de las distancia fiente-geófono
– Ganancia variable en función del tiempo (AGC) para alcanzar
amplitudes cuasi-constantes
– Amplificación binaria en modo de punto flotante

™ REGISTRO
– Almacenamiento de las trazas con su grado de amplificación
 SISMÓGRAFOS

™ Abem Terraloc Pro-2

 MANEJO DE SISMÓGRAFOS

Atlas Copco geophysics

ABEM TERRALOC SEISMIC SYSTEM

ABEM ABEM
TERRALOC TERRALOC

RAM RAM
DISK DISK

- PRI 7 8 9

+ COP 4 5 6

F MEN 1 2 3

ARM SET ENT 0 CLR

 PROGRAMACIÓN DE SIMÓGRAFOS

™ CONFIGURACIÓN DE REGISTRO
¾ Traza (polaridad) y canal
¾ Apilamientos (automático, ratificación, rápido, único, …)
¾ Ventanas de ruido
¾ Filtros analógicos

™ VISUALIZACIÓN
¾ Tipos: Total, promedio, limitada, AGC, normalizada …
¾ Escala de representación
¾ Trazo y Sombreado
 VISUALIZACIÓN

TOTAL / PROM / LIMIT / AGC NORMALIZADA

 APROBACIÓN DE LAS TRAZAS

™ MODIFICAR MODO Y ESCALA DE REPRESENTACIÓN HASTA

MÁXIMA RESOLUCIÓN DE PRIMERAS LLEGADAS

™ TRAZAR RECTAS (DIRECTA Y REFRACTADAS) Y ANOTAR

VELOCIDADES

™ MEDIR TIEMPOS FINALES Y ESTIMAR TIEMPOS

INTERMEDIOS

™ REALIZAR TIROS INTERMEDIOS E INVERSOS

™ COMPROBAR LEYES DE REFRACCIÓN Y REPETIR EN SU

CASO
TRAZAS REGISTRADAS

TIRO CENTRAL
 TECNICAS SÍSMICAS EN SONDEOS

UP–HOLE:
™ Fuente en la superficie y distintas posiciones de geófono en el
interior del sondeo.
™ Estudios de ondas directa (p y s): Determinación de la distancia que
recorren dichas ondas - Proyección de cada tiempo sobre el eje del
sondeo:
x
tk (cor) = tk · cos(arctan PT )
zk
siendo xPT la distancia entre el PT y la boca del sondeo.
™ La proximidad entre las distintas posiciones de geófonos, la
necesidad de repetir tiros cuando se tiene un solo geófono y la
resolución de los equipos y del tiempo cero, provoca cierta
indeterminación en las medidas
™ Para interpretar up-hole se debe realizar una subdivisión en tramos
suficientemente rectos para obtener posteriormente, mediante
aproximaciones lineales, las pendientes cuyos inversos
proporcionan la velocidad en cada estrato.
 UP-HOLE & CROSS-HOLE

UP-HOLE CROSS-HOLE
 TOMOGRAFÍA CRUZADA

Geofonos de superficie Geofonos de superficie

Posiciones de geofonos

Posiciones de geofonos
Puntos de Impacto
Puntos de Impacto

Tiro en posicion 2 Tiro en posicion 7