Contenido
1. Introducción ..................................................................................................... 2
2. Corrección estática .......................................................................................... 4
2.1 Corrección estática por datum ................................................................. 4
2.2 Corrección estática por refracción ............................................................ 6
3. Divergencia esférica ........................................................................................ 8
4. Filtro f-k ......................................................................................................... 10
5. Deconvolución ................................................................................................ 12
6. Filtro pasa-banda .......................................................................................... 15
7. Corrección dinámica ...................................................................................... 17
8. Análisis de velocidades.................................................................................. 19
8.1 Primera iteración. ................................................................................... 21
8.2 Segunda iteración. .................................................................................. 23
9. Conclusión ..................................................................................................... 25
10. Bibliografía ................................................................................................. 26
1
1. INTRODUCCIÓN
La Tierra no es un cuerpo estático, sino que ha estado y está sujeta a una
continua actividad, es decir, la Tierra es un sistema dinámico. Esto se debe a la
interacción de las fuerzas geológicas del interior de la Tierra. La corteza terrestre
está constituida por las llamadas placas tectónicas, las cuales, en su interacción,
producen los terremotos y volcanes. Como producto de esta interacción entre
placas, se generan los terremotos, que consiste en el movimiento de la superficie
terrestre debido a la propagación de las ondas sísmicas. La presencia de estas
ondas es detectada por los instrumentos llamados sismógrafos, que básicamente
consisten en un sistema mecánico: masa, resorte y amortiguador (un sistema de
segundo orden) y un sistema de grabación de datos, ya sea en un medio analógico
(aguja sobre tambor con papel ahumado) o mediante un medio digital (PC con
conversor analógico digital).
Con base en los principios de la sismología, la exploración sísmica emplea las
ondas elásticas que se propagan a través del terreno y que han sido generadas
artificialmente. Su objetivo es el estudio del subsuelo en general, lo cual permite
obtener información geológica de los materiales que lo conforman. La prospección
sísmica es una herramienta de investigación poderosa, ya que con ella se puede
inspeccionar con buena resolución desde los primeros metros del terreno (sísmica
de alta resolución o sísmica superficial; shallow seismic) hasta varios kilómetros
de profundidad (sísmica profunda; deep seismic). Así, para la sísmica profunda
se utilizan fuentes de energía muy potentes (explosivos o camiones vibradores)
capaces de generar ondas elásticas que llegan a las capas profundas del subsuelo,
mientras que para la sísmica superficial se utilizan martillos de impacto, rifles
sísmicos y explosivos de baja energía. De manera que el diseño de una campaña
sísmica (equipo y material a utilizar) está en función del objetivo del estudio.
Según esto, la sísmica profunda se emplea en la detección de reservorios
petrolíferos (ya sea terrestre o marítima), grandes estructuras geológicas
(plegamientos montañosos, zonas de subducción, etc.), yacimientos minerales,
domos salinos, etc. Mientras que la sísmica superficial tiene mucha aplicación
en la obra pública y la ingeniería civil. La prospección sísmica se basa en el
mismo principio que la sismología, consiste en generar ondas sísmicas mediante
una fuente emisora y registrarlas en una serie de estaciones sensoras (geófonos)
distribuidas sobre el terreno. A partir del estudio de las distintas formas de onda
y sus tiempos de trayecto, se consiguen obtener imágenes del subsuelo que luego
se relacionan con las capas geológicas (secciones sísmicas, campos de
velocidades, etc.). El desarrollo de la teoría sísmica se remonta a 1678 cuando se
enuncia la Ley de la Elasticidad de Hooke, mucho antes de la existencia de
2
instrumentos capaces de realizar medidas significativas. Sin embargo, no es sino
hasta 1845 cuando, Robert Mallet, realiza los primeros intentos de medición de
las velocidades sísmicas a través de “terremotos artificiales”, usando pólvora
negra como fuente de energía y recipientes de mercurio como receptores. En 1899
Knott desarrolla la teoría sísmica de la reflexión y la refracción. Pero, es en 1910
cuando las diferencias entre las ondas S y P se da a conocer por A. Mohorovicic,
quien las identifica y las relaciona con la base de la corteza, el Moho.
La sísmica de reflexión nace gracias a los primeros trabajos realizados por
Reginald Fesseden, en 1913, con el fin de detectar icebergs. Pero no fue sino
hasta 1927 cuando el método de reflexión se convierte en una técnica comercial
de exploración geofísica. En 1919, Ludger Mindtrop aplicó para una patente
sobre el método de refracción y ya hacia 1930 todos los domos salinos
superficiales habían sido descubiertos mediante esta técnica de exploración.
Rieber (1939) introduce la idea del procesado de datos sísmicos usando una
grabación de densidad variable y foto celdas para la reproducción de las trazas
sísmicas. Sin embargo, es en 1953, cuando las cintas magnéticas se hicieron
comercialmente disponibles, que se dio el paso al inicio del procesamiento de
datos; difundiéndose rápidamente en los años siguientes. Hasta este momento
no se empleaba la geometría CMP, la cual es usada por primera vez en 1956. A
finales de los 70, coincidiendo con el auge informático y el desarrollo tecnológico,
los nuevos soportes digitales y la nueva instrumentación representaron otro
cambio significativo en el campo de la sísmica. Desde entonces no se ha dejado
de trabajar en la continua mejora de las técnicas de adquisición y procesamiento
de datos. En la actualidad toda la adquisición se realiza en formato digital y los
datos son procesados antes de su interpretación lo cual hace al procesado de
datos algo fundamental para la localización de estructuras y así obtener
resultados certeros en los objetivos deseados.
3
2. CORRECCIÓN ESTÁTICA
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ático.
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 CORRECCIÓN ESTÁTICA POR DATUM
Se conoce como Corrección Estática a las correcciones verticales invariante en
tiempo, estas aplicadas para colocar fuentes y detectores sobre un mismo plano
horizontal imaginario (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
4
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
𝐸𝑙𝑒𝑣.𝑓𝑢𝑒𝑛𝑡𝑒 − 𝐸𝑙𝑒𝑣.𝑁𝑅 𝐸𝑙𝑒𝑣.𝑟𝑒𝑐𝑒𝑝𝑡𝑜𝑟 − 𝐸𝑙𝑒𝑣.𝑁𝑅
𝐶𝑒𝑠𝑡 = +
𝑉 𝑉
𝐸𝑙𝑒𝑣.𝑓𝑢𝑒𝑛𝑡𝑒 − 𝐸𝑠𝑝.𝑤 𝐸𝑙𝑒𝑣.𝑤 − 𝐸𝑙𝑒𝑣.𝑁𝑅
𝐶𝑒𝑠𝑡 = ( + )
𝑉𝑤 𝑉𝑤
𝐸𝑙𝑒𝑣.𝑟𝑒𝑐𝑒𝑝𝑡𝑜𝑟 − 𝐸𝑠𝑝.𝑤 𝐸𝑙𝑒𝑣.𝑤 − 𝐸𝑙𝑒𝑣.𝑁𝑅
+( + )
𝑉𝑤 𝑉1
Donde:
𝐶𝑒𝑠𝑡 = 𝐶𝑜𝑟𝑟𝑒𝑐𝑐𝑖ó𝑛 𝐸𝑠𝑡á𝑡𝑖𝑐𝑎
𝐸𝑙𝑒𝑣.𝑓𝑢𝑒𝑛𝑡𝑒 = 𝐸𝑙𝑒𝑣𝑎𝑐𝑖ó𝑛 𝑑𝑒 𝑓𝑢𝑒𝑛𝑡𝑒
𝐸𝑙𝑒𝑣.𝑁𝑅 = 𝐸𝑙𝑒𝑣𝑎𝑐𝑖ó𝑛 𝑑𝑒𝑙 𝑁𝑅 𝑜 𝐷𝑎𝑡𝑢𝑚
𝐸𝑠𝑝.𝑤 = 𝐸𝑠𝑝𝑒𝑠𝑜𝑟 𝑑𝑒 𝑙𝑎 𝑐𝑎𝑝𝑎 𝑑𝑒 𝑖𝑛𝑡𝑒𝑚𝑝𝑒𝑟𝑖𝑠𝑚𝑜
𝐸𝑙𝑒𝑣.𝑤 = 𝐸𝑙𝑒𝑣𝑎𝑐𝑖ó𝑛 𝑑𝑒 𝑙𝑎 𝑐𝑎𝑝𝑎 𝑑𝑒 𝑖𝑛𝑡𝑒𝑚𝑝𝑒𝑟𝑖𝑠𝑚𝑜
𝑉𝑤 = 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑑𝑒 𝑙𝑎 𝑐𝑎𝑝𝑎 𝑑𝑒 𝑖𝑛𝑡𝑒𝑚𝑝𝑒𝑟𝑖𝑠𝑚𝑜
𝑉1 = 𝑉𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑑𝑒 𝑙𝑎 𝑐𝑎𝑝𝑎 1
5
2.2 CORRECCIÓN ESTÁTICA POR REFRACCIÓN
La corrección estática por refracción consiste en calcular el valor estática con
base a los primeros quiebres de tiempo de arribo de los registros de campo.
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 funte-
receptor no muy grande, para tener la seguridad que se está en una zona donde
se tiene 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 puedes
variar significativamente de una interfaz 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 reflector, en la cual, para las
correcciones estáticas, está definida como una interface entre la capa de
intemperismo hasta la cima del basamento.
6
Ilustración 2.1.- Tiro de control 40, con corrección estática por refracción aplicada.
En la ilustración 1.1 se puede observar el tiro de control 40 con la corrección
estática por refracción aplicada, para llegar a este resultado se tienen que
graficar las dromocronas y un modelo de capas planas (ilustración 1.2) y se
reduce los puntos donde la pendiente cause picos pronunciados. Una vez
eliminado estos picos modifican las capas a una más realistas (ilustración 1.3).
Ilustración 2.2.- Dromocronas y modelo de capas planas.
7
Ilustración 2.3.- Dromocronas y modelo de capas con corrección por refracción.
3. DIVERGENCIA ESFÉRICA
Podemos considerar a la energía sísmica como un frente de onda (esfera de
energía expandiéndose desde el punto de tiro en todas direcciones), o una serie
de rayos (líneas dibujadas perpendicularmente al frente de onda que muestran
la dirección de la propagación de la energía) emanados de un tiro o disparo, la
perturbación en estos puntos tiene la misma fase, cuando el movimiento
ondulatorio alcanza los puntos que componen a un frente de onda, cada partícula
del frente se convierte en una fuente secundaria de ondas, este proceso se repite,
lo que hace posible la propagación de la energía a través de un medio.
La corrección por divergencia esférica se realiza debido a que la energía se
distribuye radialmente, mientras los frentes de onda se alejan de la fuente (Fig.
4.) esta disipación de la energía es más fuerte cuando el área de dispersión
aumenta, generando que las amplitudes de la señal disminuyan al alejarse de la
fuente, esto es, a mayor offset la energía decrece y por consecuencia las
amplitudes registradas también. La amplitud es inversamente proporcional al
offset cuando se asume un modelo de velocidad constante (Divergencia esférica=
1/distancia) en un medio homogéneo e isotrópico, pero los contrastes litológicos
son variados en la realidad y las capas más profundas poseen una mayor
compactación (lo que genera una mayor velocidad) debido a la presión litostática,
por lo cual la recuperación de las amplitudes para las líneas de este trabajo se
realizó asumiendo un aumento en la velocidad de propagación de la energía con
la profundidad.
8
Figura 4. Propagación de los frentes de onda
(Tomada de Nordmann, 2007).
En la (ilustración 2.1) se muestra el tiro de control 40 donde se observa que a
partir de los 3000 milisegundos se pierde la mayoría se información, debido a la
distancia que ha viajado la onda generando la perdida de energía.
Ilustración 3.1.- Tiro de control 40.
9
Después de haber hecho una corrección por divergencia esférica (ilustración 2.2),
podemos observar datos después de los 3000 milisegundos, llegando a tener
información más allá de los 5000 milisegundos.
Ilustración 3.2.- Tiro 40 de control con corrección por divergencia esférica.
4. FILTRO F-K
Este tipo de filtro es útil para eliminar el ruido coherente que presenta una
tendencia lineal. El filtro se basa en la transformada de Fourier 2D de los datos
sísmicos. Los datos pasan del espacio (t-x) al espacio (f-k). Los eventos que poseen
la misma inclinación en el espacio (t-x), sin importar su localización, serán
mapeados en una línea recta en dirección radial en el espacio (f-k).
A este filtro también se le conoce también como filtro de velocidad ya que en el
espacio en que opera (frecuencia, número de onda), se discriminan los distintos
eventos por estar alineados según rectas cuyas pendientes definen las distintas
velocidades. De esta forma los eventos lineales de baja velocidad en el espacio (t-
x), son mapeados como rectas con ángulos menores, mientras que los eventos
lineales de alta velocidad se encuentran en sectores angulares mayores.
Normalmente los eventos lineales en el espacio (f-k) son discriminados haciendo
10
uso de polígonos, los cuales determinan los eventos a eliminar. La ilustración 3.1
muestra el tiro de control 40, y cómo se atenúan los ruidos lineales en el espacio
f-k, resaltando de esta forma las reflexiones.
a)
b)
Ilustración 4.1.- a) Tiro de control 40 sin filtro f-k
b) Tiro de control 40 con filtro f-k
11
5. DECONVOLUCIÓN
Para entender el concepto de deconvolución se debe tener claro que es el modelo
convolucional de la tierra (Fig. 5), el cual radica en lo siguiente: las señales
sísmicas, al introducirse en el subsuelo se ven afectadas por reflexiones,
refracciones, reverberaciones, difracciones, múltiples, reflexiones en interfaces
de estratos delgados, absorción no lineal y ruido adicional, parte de energía que
viaja hacia el subsuelo, es reflejada al encontrar cambios en la densidad y la
velocidad (contrastes de IA) entre capas de diferente composición, a esta energía
que regresa se le conoce como coeficiente de reflexión, la energía que continua su
camino es el coeficiente de transmisión.
De tal manera que, finalmente la señal registrada en superficie es una
convolución del impulso producido por la fuente con los factores mencionados, lo
que produce cambios a la señal de entrada: 𝑆(𝑡) = 𝑘δ(t) a(t)z(t) + u(t)
Donde:
𝑆(𝑡) = 𝑂𝑛𝑑í𝑐𝑢𝑙𝑎
𝑘δ(t) = 𝐼𝑚𝑝𝑢𝑙𝑠𝑜 − δ
a(t) = Absorciónrápidadealtasfrecuencias
z(t) = Otros factores
u(t) = 𝑅𝑢𝑖𝑑𝑜
El coeficiente de reflexión está dado por la siguiente fórmula:
ρ2v2 − ρ1v1
𝐶𝑅 =
ρ2v2 + ρ1v1
Donde:
ρ = Es la densidad de la capa
v = es la velocidad de la capa
La cantidad de energía reflejante es el coeficiente entre la diferencia de las
impedancias acústicas y su suma.
12
La deconvolución consiste en la compresión de la ondícula básica, con el propósito
de eliminar el efecto generado por la fuente de energía. Este proceso limpia la
señal para que se pueda reducir a un spike o pico, el cual es la respuesta
reflectiva del medio sin la alteración generada por el impulso de entrada,
mejorando la resolución temporal, atenuando o removiendo ruido coherente, y
reverberaciones de la fuente generadas por la capa de agua.
Figura 5. Modelo de convolución tomando en cuenta solo la respuesta
reflectiva de las capas.
13
a)
b)
Ilustración 5.1.- a) Tiro de control 40 antes del proceso de deconvolución
b) Tiro de control 40 después del proceso de deconvolución.
Los óvalos muestran las zonas donde es más apreciable los cambios efectuados
por el proceso de deconvolución.
14
6. FILTRO PASA-BANDA
Es uno de los filtros más empleados, tiene como finalidad dejar pasar la señal en
una banda limitada de frecuencias de manera que se aceptan las frecuencias que
contienen energía de reflexión coherente y se rechaza aquellas frecuencias
asociadas al ruido sísmico (ondas superficiales, aérea, ruido, ambiental, etc.).
Figura 6. Configuración de filtro pasa-banda
15
a)
b)
Ilustración 6.1.- a) Tiro de control 40, sin filtro pasa-banda
b) Tiro de control 40, después de filtro pasa-banda.
Los óvalos muestran las zonas donde es más apreciable los cambios efectuados
por el proceso del filtro pasa-banda, donde se eliminan las frecuencias fuera de
rango, dando esa sensación de tener más tendencia.
16
7. CORRECCIÓN DINÁMICA
Antes de que las trazas sean apiladas (sumadas), la reflexión necesita ser
ajustada a su verdadero tiempo de viaje de ida y vuelta, en otras palabras, una
corrección de tiempo debe ser efectuada para compensar el tiempo de viaje de la
componente horizontal de velocidad de la onda sísmica grabada, lejano a la
estación de disparo (Manzoni, 2001).
Entonces, este tiempo hay que corregirlo a través de la ecuación:
∆𝑡𝑛𝑚𝑜 = 𝑡(𝑥) − 𝑡0
Donde 𝑡(𝑥) es el tiempo registrado, ∆𝑡𝑛𝑚𝑜 es la corrección y 𝑡0 el tiempo vertical
e la reflexión (de ida y vuelta). Recordemos que el tiempo registrado o tiempo de
viaje viene expresada por:
𝑥2
𝑡 2 (𝑥) = 𝑡02 +
𝑉2
Donde x es la distancia entre las posiciones de la fuente y el receptor, v es la
velocidad del medio que se encuentran por sobre la interface reflectora y 𝑡0 es
dos veces el tiempo de viaje, para esto es útil la función de velocidad, ya que ella
nos da información de la velocidad del medio, la cual no es conocida a priori, pero
a través del análisis de velocidad como se verá más adelante, se puede
determinar la velocidad que logra corregir cada evento geológico primario.
La corrección NMO de la diferencia entre t(x) y t0 quedaría entonces como:
2 𝑥2
∆𝑡𝑛𝑚𝑜 = √(𝑡02 + ) − 𝑡0
𝑉2
Si no se aplica dicha corrección, las reflexiones se cancelarán cuando el CDP, sea
sumado o apilado, debido a que un CDP está constituido por trazas de diferentes
aperturas o distancias entre puntos fuentes y receptores (offset), generando
interferencia destructiva, con lo cual no se apreciará una amplitud resultante,
(Manzoni, 2001).
En las trazas de amplia cobertura, la señal sísmica, será reforzada en el
momento de su apilamiento y los ruidos tenderán a sumarse destructivamente,
pero la corrección dinámica, es también un proceso que logra atenuar los
múltiples y reverberaciones.
La traza sísmica que resulta de la aplicación de la corrección dinámica, (NMO)
después del apilamiento, se coloca en el punto medio entre las fuentes y los
17
receptores. Esto funciona, cuando las interfaces de reflexión son horizontales,
por lo tanto, las velocidades de apilamiento son iguales a las velocidades
cuadráticas medias (VRMS) de cada capa geológica. Cuando la capa es inclinada,
la traza es posicionada también en el punto medio entre la fuente y receptor,
pero no pasa por el punto de reflexión (figura 7), (Pereira, J., 1987) Esta situación
es una de las razones por la que las secciones sísmicas que se obtienen en áreas
muy plegadas geológicamente estén deformadas y requieran del proceso de
migración.
Figura 7. Corrección por NMO de las trayectorias fuente-receptor para una capa buzante
(NMO), (Modificado de Pereira, 1987)
La figura 7 muestra una capa buzante, la geometría de este tipo, en este caso
según Levin (1971), el tiempo de viaje de la fuente C al reflector en el punto O,
hasta llegar al receptor 32 para una capa que buza con un ángulo φ viene
expresada como:
𝑥 2 𝑐𝑜𝑠 2 𝜑
𝑡 2 (𝑥)
= 𝑡02 +
𝑉2
18
Esta ecuación corresponde ciertamente a una hipérbola, sin embargo, la
velocidad para la corrección NMO viene dado por la velocidad del medio dividido
por el coseno del ángulo buzante (para un plano en 2D):
𝑉
𝑉𝑛𝑚𝑜 =
𝑐𝑜𝑠𝜑
Considerando estos aspectos, la corrección NMO de la diferencia entre t(x) y t0
quedaría entonces como:
2 𝑥 2 𝑐𝑜𝑠 2 𝜑
∆𝑡𝑛𝑚𝑜 = √(𝑡02 + ) − 𝑡0
𝑉𝑛𝑚𝑜 2
8. ANÁLISIS DE VELOCIDADES
Antes de que las trazas sean sumadas o apiladas, la reflexión de la onda sísmica
necesita ser ajustada a su verdadero tiempo de viaje, y esto se hace haciendo
correcciones dinámicas (o de NMO). Para ello se necesita de la información de la
velocidad con la que viajan las ondas sísmicas, y la forma de obtener esta
información es a través del proceso de estimación de velocidades o análisis de
velocidad. Mediante el análisis de velocidad, se puede determinar la velocidad
que logra corregir cada evento geológico primario, (Manzoni, 2001), de tal forma
que se pueda llevar este, a una posición donde la reflexión obtenida fuese
equivalente a la de tener a la fuente en el mismo sitio que el receptor (cero offset),
(Pereira, J., 1987). A esta velocidad se le denomina velocidad de apilamiento.
El análisis de velocidad tiene por objeto de seleccionar las funciones de velocidad,
como función del tiempo, las cuales se necesitarán para hacer las correcciones
dinámicas para poder realizar un apilado preliminar. Por otro lado, las trazas de
una sección sísmica apilada son equivalentes a la utilización de distancias
fuentereceptor, igual a cero, (cero offset).
Para el análisis de velocidad se necesitan los datos más reales posibles, de esta
forma, lograremos comparar estos resultados y ver en qué medida se han ido
mejorando las imágenes. Para ello se necesitan los datos previamente
deconvolucionados y corregidos estáticamente, y con una pequeña recuperación
de ganancia (Automatic Gain Control, AGC), podremos obtener una imagen
aceptable para realizar este análisis de velocidad.
19
La mejor forma de estimar las velocidades es utilizando una medida numérica
de la coherencia de cada una de las trazas que conforman la traza apilada. Una
medida común de coherencia es la de semblanza. La semblanza no es más que la
relación que existe entre la energía de las trazas apiladas con la suma de las
energías de las trazas individuales dentro de una ventana de tiempo.
Ilustración 8.1. Primer modelo de velocidades.
20
8.1 PRIMERA ITERACIÓN.
Para llevar a cabo este análisis, los grupos de puntos comunes en profundidad
(CDPs) en que se deben agrupar para la obtención de una información
suficientemente representativa, depende de análisis que se desea hacer y de las
reflexiones geológicas que se desean considerar en el análisis, por lo general se
debería hacer junto con el intérprete geológico, de tal forma de escoger las
reflexiones más importantes de la zona de estudio.
Es importante que cuando se realicen los piques o el análisis de velocidad y se
detecten la presencia de múltiples, se debe seleccionar un punto o una función
de velocidad que no corrija dinámicamente esos múltiples, los cuales al no ser
apilados se atenuaran mediante el proceso de cancelación, en la suma de las
trazas de un punto común en profundidad. Los múltiples se pueden reconocer
mediante el análisis de coherencia, ya que como se dijo anteriormente,
usualmente los múltiples, presentan coherencia a una velocidad inferior a la de
los eventos primarios, o presentan igual velocidad que el evento primario a un
tiempo mayor de reflexión.
21
Ilustración 8.2. Segundo modelo de velocidades, después de una primera iteración.
Dentro del ovalo azul se puede apreciar cambios que comienzan a ser
considerables dentro de la tendencia.
22
8.2 SEGUNDA ITERACIÓN.
En su segunda iteración es necesario realizar un segundo pase de análisis de
velocidad, sobre los datos que previamente han sido corregidos, esto se hace con
el fin de obtener una función de velocidades que corresponda a los datos reales
de la geología de la zona. En consecuencia, con este segundo paso de análisis de
velocidades, será mucho más fácil seleccionar y diseñar una nueva función de
velocidades, que corresponde a la verdadera información de la geología, puesto
que, al incrementar la coherencia lateral de una reflexión sísmica, también
aumentará el nivel de energía convergente de las amplitudes en esas reflexiones
sísmicas.
La selección de esta segunda función de velocidad, por punto común en
profundidad, se hará de manera similar a como se realizó para el primer pase de
análisis de velocidad, para lo cual se necesitará los datos previamente corregidos
por estáticas de campo y posteriormente corregidos por las primeras estáticas
residuales. Luego que se haga la selección de la segunda función de velocidad, se
deberá implementar análogamente a como se realizó para el primer paso de
análisis de velocidades, un control de calidad de estas, mediante la
representación de la función en paneles de isovelocidades.
Aquí, se podrá establecer una comparación entre los paneles de isovelocidades
producto del primer análisis de velocidades y los paneles del segundo análisis de
velocidad, en ellos se podrá hacer una evaluación de los procesos que hasta aquí
hemos realizados, 120 además de visualizar las variaciones de ambas funciones.
De esta forma podemos asegurar la buena selección de las funciones de velocidad.
23
Ilustración 8.3. Tercer modelo de velocidades, a partir de una segunda iteración.
Dentro del ovalo azul se puede apreciar como ha habido una mejora contundente
comparado con los dos modelos anteriores (ilustración 9.1 e ilustración 9.2).
24
9. CONCLUSIÓN
Es importante que cada uno de los procesos aplicados durante el procesado de
datos sean los correctos en función del objetivo a estudiar, existen diferentes
procesos que pueden realizarse a nuestros sismogramas, sin embargo, los
mostrados en este trabajo son una base para obtener un modelo de isovelocidades
que nos proporciona información sobre las características de nuestra de zona de
interés.
Para verificar los piques de velocidad realizados en el último proceso, el
programa permite mostrar los valores numéricos de los piques de velocidad en
un diagrama de contorno, llamado también gráfico de isovelocidades (ilustración
9.1).
Estos mapas de contorno representan la tendencia y las variaciones de las
velocidades de las capas de rocas en el subsuelo, en ellos se pueden apreciar la
coherencia de las funciones de velocidad seleccionadas por el analista geofísico y
determinar la existencia de posibles valores anómalos, que pudieran dar lugar a
un apilado erróneo. Abruptos cambios en la tendencia de estas líneas de
isovelocidades, son indicativos de variaciones laterales de velocidad, producto de
la actividad geológica que hubo en el área.
En la ilustración 9.2 podemos ver los datos luego de aplicada la primera
iteración, en ella se puede observar que los eventos presentan una mejora en
cuanto a la continuidad de los eventos, lo cual quiere decir que el proceso es
adecuado. De esta manera podemos actualizar el último modelo.
Luego de este segundo modelo, en la ilustración 9.3 podemos ver los datos
después de un segundo proceso de iteración, en el, de nuevo se puede observar
una mejora en cuanto a la continuidad de los eventos, sobre todo en las partes
más profundas, lo cual quiere decir que tanto los procesos llevados a cabo como
el modelo de isovelocidades respectivo a una segunda iteración han sido
satisfactorios.
25
10. BIBLIOGRAFÍA
- Gallardo, A. 2003. Procesamiento de datos sísmicos sobre topografía irregular al noreste
del estado de Anzoátegui (Tesis para el grado de Ingeniero). Universidad Central de
Venezuela, Venezuela.
- Jiménez, P. Optimización en la secuencia de procesamiento de datos sísmicos marinos
(Tesis para el grado de Ingeniero). Facultad de Ciencias, Universidad Nacional de
Ingeniería, Perú.
- Manzoni G., S (2001). Secuencia básica de procesamiento de Data Sísmica para obtener
datos finales orientada a estudiantes de Ingeniería Geofísica (Trabajo especial de Grado).
Universidad Central de Venezuela, Venezuela.
- Pereira, J. (1994). Algunas Ideas para Mejorar la Calidad de los Datos Sísmicos. Memoria
VII. Congreso Venezolano de Geofísica, Venezuela.
- Teixidó-iUllod, Ma. T. 2000. Caracterización del subsuelo mediante sísmica de reflexión
de alta resolución (Tesis Doctoral). Departamento de Geodinámica y Geofísica,
Universidad de Barcelona, España.
- Vargas, Y. 2011. Estructura en el sector noroeste del golfo de california a partir de
sísmica de reflexión (Tesis para el grado de Maestro). Facultad ciencias de la tierra,
Universidad Autónoma de Nuevo León, México.
- Zarate, A. 2011. Análisis de la corrección estática en el procesado sísmico de reflexión
(Tesis para el grado de Ingeniero). Escuela Superior de Ingeniería y Arquitectura Unidad
Ticomán, Instituto Politécnico Nacional, México.
26
