FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
UNIDAD I
Capítulo V: MÉTODOS DE INTERPRETACIÓN DE REFRACCIÓN SÍSMICA
Práctica 9. Problemas sobre correcciones en refracción sísmica
1. INTRODUCCIÓN
El procesamiento sísmico consiste en la elección de los parámetros y algoritmos adecuados para
su aplicación a datos sísmicos adquiridos en el campo, con el fin de obtener secciones sísmicas
sin errores por efectos topográficos. Para obtener un procesamiento de datos convencional se
necesita los procesos de: pre-apilamiento, apilamiento y post-apilamiento.
Dentro de la etapa de pre-apilamiento es necesario resaltar la corrección estática. En dicha
corrección, normalmente se asignan valores de velocidad para la capa de intemperismo y la capa
firme, de acuerdo con los datos obtenidos en las pruebas de campo. Sin embargo, estos valores
pueden cambiar lateralmente, ya sea por una topografía abrupta o por cambio de los materiales
que conforman dichas capas.
1.1 OBJETIVO
• Determinar la velocidad de las ondas de las capas superficiales meteorizadas y
promediarlas
• Calcular los espesores de las capas meteorizadas y promediarlas para sus efectos de
cálculo
• Realizar las correcciones estáticas en Excel
2. METODOLOGÍA
Un factor fundamental y necesario para el proceso de datos sísmicos es la exactitud de la
medida de tiempo, considerando que las posiciones de los puntos de tiro y receptores varían a
lo largo de la línea debido a la topografía del terreno; esto afecta a todos los reflectores
encontrados a diferentes profundidades que están representando imágenes de estructuras
influenciadas por topografía, por lo tanto, es necesario corregir estas falsas representaciones en
tiempo. Para lograr esto se corrigen todos los datos a un plano denominado nivel de referencia
de cierta elevación, con lo que, todos los puntos de tiro y localización de receptores deben tener
un tiempo de corrección asociado, este tiempo de corrección es llamado ESTÀTICA.
O también se puede considerar las refracciones estáticas ya que permiten determinar espesores
y las velocidades cercanas a la superficie al efectuar el análisis de los primeros arribos. Con
esto se puede definir dos tipos de corrección estática: corrección estática por DATUM y
corrección estática por refracción.
2.1 CORRECCION ESTATICA POR DATUM.
Se conocen como Corrección Estática a las correcciones verticales invariantes en tiempo, estas
aplicadas para colocar fuentes y detectores sobre un mismo plano horizontal imaginario
FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
(Datum). Tiene como objetivo cancelar los efectos sísmicamente indeseables de la porción
superficial de la Tierra, que afectan a todas las ondas que viajan a través de los estratos
superficiales. Estos efectos son debidos a variaciones de elevaciones, variaciones erráticas del
estrato superficial (capa de intemperismo), sean de espesor o de cambios laterales de velocidad.
FIGURA 1. Muestra gráficamente la corrección estática por Datum
Figura 2. a) Modelo antes de la corrección estática por Datum, b) modelo después de la corrección
estática por Datum.
2.2 CORRECCIÓN ESTÁTICA POR REFRACCIÓN.
La corrección estática por refracción consiste en calcular el valor de estática en base a los
primeros quiebres de tiempo de arribo de los registros de campo.
FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
Para efectuar este cálculo se cuenta con la posición y elevación de los puntos de tiro y
receptores, con esto se puede deducir la velocidad y profundidad de las capas cercanas a la
superficie, para un mejor resultado se toman distancia fuente-detector no muy grande, para tener
la seguridad que se está en una zona donde se tienen refracciones, ya que a distancias mayores
llegan tanto ondas reflejadas y refractadas.
Figura 3 Comportamientos de la onda conforme a la corrección.
En la práctica una capa simple del modelo superficial frecuentemente es suficiente para resolver
las anomalías estáticas. Sin embargo, las capas pueden variar significativamente de una
interface plana a una interface de forma irregular. La suposición de la capa simple para la
superficie cercana tampoco funciona cuando hay un cambio lateral en la composición de la roca
asociado con afloramientos. Las complejidades de una capa simple del modelo superficial
pueden deberse a los siguientes casos:
a. Variaciones rápidas en la elevación de fuentes y receptores.
b. Variaciones laterales en la velocidad de la capa de intemperismo.
c. Variaciones laterales en la geometría del refractor, en la cual, para las correcciones estáticas,
está definida como una interface entre la capa de intemperismo hasta la cima del basamento.
2.2.1 Primeros quiebres (first breaks).
La energía refractada asociada con la base de la capa de intemperismo frecuentemente
constituye los primeros arribos en el registro de tiro común. Los primeros quiebres (FB) ocurren
con distintos grados de variabilidad, dependiendo del tipo de fuente y de las condiciones de la
superficie. Son distorsionados en la presencia de irregularidades en la base de la capa de
intemperismo y esto se vuelve en una onda dividida cuando no hay un fuerte contraste de
velocidad entre la capa de intemperismo y el subestrato
La desviación de la trayectoria lineal del tiempo de arribo de los primeros quiebres puede
atribuirse a las diferencias de las elevaciones a lo largo del perfil (línea). Así que de los primeros
quiebres pueden inferirse irregularidades en la superficie, así como un refractor de forma
variable o variaciones de velocidad en la capa superficial.
FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
Estos pueden seleccionarse en forma automática, interactiva o manual, o como una
combinación de ambos. Para hacer los FB confiables, primero es necesario aplicar la corrección
lineal (LMO) a los datos. Una vez que se seleccionaron los FB hay que guardar la corrección
LMO. Note que la efectividad de ambos reflexión y refracción basada en el método de corrección
estática depende de la confiabilidad del proceso de selección de FBs.
La energía refractada está asociada con los “los frentes de onda”. Si los arribos refractados son
observados en un gather de tiro común, la geometría de la capa superficial se supone simple,
sin embargo, la teoría de rayo no puede afirmarse para variaciones de la longitud de onda corta
en base a la capa de intemperismo que es mucho más pequeña que la longitud del cable. Estas
variaciones son manipuladas en una corrección subsecuente (llamada estáticas residuales
usando la distorsión de tiempo de viaje asociado con reflexiones en CMP gather corregido por
NMO).
2.2.2 Comportamiento de los primeros arribos
Curvas de tiempo de tránsito, La sísmica de refracción utiliza como materia prima los llamados
primeros arribos. El objetivo de la sismología de refracción dentro de un estudio de reflexión es
determinar las variaciones de velocidad en función de la profundidad para las capas someras.
Figura 4 Muestra la relación entre arribos de refracción y reflexión con el mismo tiempo de viaje.
FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
Debido al aumento en la velocidad de las capas más profundas, los eventos de refracción que
viajan a mayores profundidades pueden llegar a confundirse con los eventos de reflexión de
trayectorias cortas, por lo tanto, para el estudio de refracción se consideraran únicamente los
llamados “primeros arribos” los cuales tienen trayectorias de viaje cortos.
Todos los eventos de sísmica que se describirán son Ondas-P, también llamadas ondas de
compresión u ondas de cuerpo. Los primeros eventos en llegar a los receptores siguen una
trayectoria horizontal (onda directa), los cuales viajan con la velocidad de la primera capa. (1).
En la gráfica dichos eventos definen un segmento de recta cuya pendiente es 1/v1.
El frente de ondas generado de la fuente “S” llega a la interface entre las capas 1 y 2. La
trayectoria del rayo que incide en el medio 2 con el ángulo 𝑖1,2= 𝑠𝑖𝑛−1𝑣1𝑣2⁄ por la ley de Snell,
se refracta en el segundo medio de manera que el ángulo de refracción de este medio es 90°.
Este ángulo es llamado “ángulo crítico”, 𝑖1,2=𝑖𝑐𝑟𝑖𝑡=𝑠𝑒𝑛−1(𝑣1/𝑣2). Los primeros arribos (2) son
generados en la parte superior de la segunda capa los cuales viajaran a lo largo de la parte
superior de la capa 2 con una velocidad 𝑣2 (𝑣2>𝑣1).
En la gráfica los primeros arribos se ven representados como una línea recta de pendiente 1/𝑣2.
Como ya se mencionó 𝑣2>𝑣1, así que 1/𝑣2<1/𝑣1, es decir, la pendiente de esta onda será menor
que la de la onda directa.
Las ondas que inciden con un ángulo crítico y que regresaran inmediatamente a la superficie
como una reflexión crítica (2𝑐), forman una distancia ente a y 2𝑐, que es conocida como
“distancia crítica”.
Los arribos que se encuentran entre la distancia crítica (𝑋𝑐𝑟𝑖𝑡𝑖𝑐𝑎) y la distancia de cruce (𝑋𝑐𝑟𝑜𝑠𝑠)
o distancia crossover y viajan horizontalmente en la primera capa son registrados como
segundos arribos, como se muestra en la figura 4.
2.3 RETARDOS DE TIEMPO
Los tiempos conocidos como “delay time” son los tiempos que tarda en llegar la energía a cada
horizonte refractor dependen de las velocidades de cada capa; estos se pueden utilizar para
obtener la corrección estática por refracción.
Asumiendo una vertical ficticia y un frente de onda que va de la fuente A al receptor D con
dirección horizontal (como se muestra en la figura 5). Y si se debe corregir la distancia X que
recorre con una velocidad 𝑣2 entonces:
La figura 5 muestra una corrección que no es más que la diferencia del tiempo de una trayectoria
inclinada con velocidad 𝑣1 y el tiempo a lo largo de la proyección de una trayectoria inclinada
con dirección de la onda que tiene una velocidad 𝑣2.
FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
En la figura de abajo se ilustra el concepto de retardo de tiempo. Si se asume un frente de ondas
que viaja de “a” a “d” con velocidad 𝑣2 (𝑇2= 𝑋𝑣2⁄) y que solo la trayectoria que va de “b” a “c” es
recorrida con una velocidad alta 𝑣2, entonces se debe de hacer la corrección aumentado más
tiempo. El tiempo requerido es la diferencia entre el tiempo que pasa a lo largo de la inclinación
“ab” con una velocidad de 𝑣1 y la trayectoria del componente horizontal (“eb”) que viaja a la
velocidad 𝑣2. Esta diferencia es llamada”Retardo de tiempo”. Para el refractor horizontal, el
retardo total de tiempo es el retardo de tiempo asociado con el segmento descendente del rayo
(“ab”), más el retardo de tiempo asociado con el segmento del rayo ascendente, (“cd”). Para el
segmento “bc”, la expresión matemática del tiempo de intercepción para cuando X= 0 es el
tiempo que pasa a lo largo “ab” con velocidad 𝑣2 más el tiempo que transcurre a lo largo de “ca”
con velocidad 𝑣1, (para el caso del segmento “abca” no hay ningún evento).
Figura 5. Concepto de retado de tiempo
Esta diferencia, o “corrección” es llamada retardo de tiempo (D), y para las trayectorias
descendentes es 𝐷𝑑.
De la figura 5 la ley de Snell
FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
Cuando la trayectoria de rayos es simétrica el valor del segmento es doble, y su expresión es:
La expresión de tiempo de intercepción es idéntica, por lo tanto, el retardo en tiempo es igual al
tiempo de intercepción, entonces 𝑇2 se puede expresar como:
Si se considera la figura 5 y la ecuación 2.1 el significado matemático del tiempo de intercepción
cuando X=0 es: el tiempo a lo largo de “ab” con velocidad 𝑣1 menos el tiempo a lo largo de “bc”
con velocidad 𝑣2 más el tiempo a lo largo de “ca” con velocidad 𝑣1. Esta definición también se
usa para el retardo de tiempo y se puede aplicar a varios segmentos de trayectoria de rayos en
un modelo de varias capas.
Si se desea expresar el comportamiento del retardo en tiempo de una tercera capa. Entonces
se debe de asumir que el primer arribo de la onda es generada en la parte superior de la capa
n y el retado de tiempo 𝐷3,𝑛 para este segmento está definido como:
FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
Por lo tanto, la expresión de retardo en tiempo para el segmento descendente de la interface 3
es:
La doble notación se puede utilizar, si se considera que el retardo de tiempo está asociado con
un segmento de rayo descendente de la capa profunda (figura 6), de la tercera capa o en un
modelo con n interfaces. En donde el primero número corresponde al número de interface que
se encuentran durante la trayectoria de rayos, y el segundo número a las capas (n).
FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
FIGURA 6 El cálculo del Retardo de tiempo para un segmento en la tercera capa se muestra
en esta figura. El retardo total de tiempo a lo largo de toda la trayectoria de rayos es la suma
de los retardos de cada segmento de la trayectoria de rayo. Nota el ángulo ∝3,𝑛 no es el
ángulo crítico.
3. METODOLOGIA
PASO 1, First Breaks & Coordinates, se tiene la siguiente tabla
PASO 2: Charge Depth Correction, se corrige la profundidad de carga.
SP-RCV Dist.: Separación de estación a estación
First Breack: primeros tiempos de llegada, Considere di < dc
FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
Si di > dc; ti+(Pc/Varrive )x1000
FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
PASO 3: Raw Time Breaks Graph, grafique sus datos
FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
THICKNESS CALCULATION
Relative Coordinate of Y axis for for Intercept time calculation: 0m
Direct Profile Inverse Profile
Break Points Break Points
x0= 0 y0= 0 x0= 302 y0= 0
x1= 4 y1= 6 x1= 296 y1= 8
x2= 127 y2= 88 x2= 206 y2= 69
x3= 302 y3= 165 x3= 0 y3= 165
Direct Profile Inverse Profile
Intercept Times (msec) Intercept Times (msec)
i1 3.3 i1 3.9
i2 32.1 i2 24.3
Direct Profile Inverse Profile
Velocities (m/sec) Velocities (m/sec)
V1 667 V1 750
V2 1500 V2 1475
V3 2273 V3 2146
Direct Profile Inverse Profile
Thickness (m) Thickness (m)
Z1 1.2 Z1 1.7
Z2 28.5 Z2 20.3
Average Velocities (m/sec) Average Thicknesses (m)
V1 708 Z1 1
V2 1488 Z2 24
V3 2209
Promedio de espesores de la primera capa (Zw)= Z1
Promedio de las velocidades primera capa Vw = V1
Promedio de las velocidades tercera capa Vsw = V3
PASO 3: Static Correction, realizar la corrección estática en su tabla de Excel, siguiendo el orden
señalado
FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
Volver a graficar los perfiles de topografía y de tiempo:
FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
3. MATERIALES
• Software de Excel o cualquier software de hoja calculo libre
• Plantilla de Excel para el desarrollo
• Tabla de velocidades con elevación de cada receptor.
4. PRÁCTICA
La figura de abajo se muestra las primeras llegadas a estaciones de geófonos
separadas 100 m de fuentes a 25 m de profundidad en cada extremo de la dispersión.
De fuentes a 25 m de profundidad en cada extremo de la extensión. El grupo de
FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
geófonos de cada extremo no se registra debido al ruido del pozo. El tiempo del pozo
está en la tercera de la derecha. Las elevaciones de cada grupo se indican en la parte
superior.es de 500 m/s.
1. Calcular su distancia critica
2. Realizar las correcciones estáticas de tiempos
3. Realizar las correcciones estáticas por elevación
4. Calcular las velocidades y esposares de cada capa, con los nuevos datos
corregidos
5. Determinar las capas meteorizadas
XA(m) XB(m) E(m) tA(ms) tB(ms)
0 1000 1125 — 433
100 900 1130 66 398
200 800 1127 106 357
300 700 1128 151 317
400 600 1125 184 270
500 500 1120 218 220
600 400 1120 257 172
700 300 1125 313 143
800 200 1133 367 118
900 100 1138 418 83
1000 0 1140 466 —
FACULTAD DE GEOLOGIA, GEOFÍSICA Y MINAS
ESCUELA PROFESIONAL DE INGENIERIA GEOLÓGICA
PROSPECCIÓN GEOFÍSICA I
5. BIBLIOGRAFÍA DE REFERENCIA
- Cantos Figuerola, J., 1972. Tratado de geofísica aplicada (p 204-234).
- Burger, H., 1992. Exploration geophysics of the shallow subsurface (p 57-133).
Prentice Hall.
- Dobrin, M., 1976. Introduction to geophysical prospecting (p 292-338). McGraw-Hill.
- Dong-Joo, M. y Changsoo, S., 2006. Refraction tomography using a waveform-
inversion back-propagation technique. Geopyhsics, V.71 Nº3 (p 21-30).
- Sheriff, R., 1991. Encyclopedic dictionary of exploration geophysics (376 p). Society
of Exploration Geophysicists.
- Sheriff, R. y Geldart, L., 1995. Exploration seismology (p 95-106, 425-450). Cambridge
University Press.
- Smith, P., 1975. Temas de geofísica (p 26-53). Editorial Reverté.
- Telford, W., Geldart, L., Sheriff, R. y Keys, D., 1976. Applied geophysics (p 363-369 y
415-425). Cambridge University Press.
