INTRODUCCIÓN
El procesamiento de datos sísmicos es claramente una parte vital en la
exploración sismológica, ya que recae en él tres etapas a realizar para cumplir el
propósito deseado.
Hablamos de la adquisición, el procesamiento y la interpretación. Donde,
conceptualizando a cada uno de ellos tenemos que el primero consiste en generar
una fuente de energía la cual se propagará en el subsuelo, y así, se obtiene
información que se registra en equipos especiales dependiendo la zona de estudio
donde se graba una señal sísmica, donde al final, obtendremos un sismograma
para después analizar las trazas sísmicas generadas. Describiendo ahora el
procesamiento, en el método sísmico de reflexión es una etapa compleja donde se
hacen diversos procesos basados en la naturaleza que los rige; en esta etapa se
recibe el producto del proceso anterior, o sea, los sismogramas, información
topográfica, etc. En el levantamiento de reflexión sísmica, las trazas sísmicas se
registran y el propósito del procesamiento sísmico puede verse como un intento de
reconstruir las diversas columnas de la Figura 4.6, moviéndose de derecha a
izquierda.
Eliminar el ruido.
Determinar el pulso de entrada y eliminarlo para dar el coeficiente de reflexión.
Determinar la función de velocidad para permitir la conversión de tiempo al eje
de profundidad
Determinación de las impedancias acústicas (o propiedades relacionadas) de
las formaciones.
Figura 1. El modelo convolucional de
la traza sísmica de reflexión, que
muestra la traza como la salida
convolucionada de una función de
reflectividad a las propiedades físicas
de las capas geológicas.
Al final del procesamiento, el resultado obtenido serán las secciones sísmicas a
partir de los sismogramas obtenidos y así generar una imagen que pueda ser
interpretada después. Donde, al interpretar estas secciones sísmicas se puede
extrapolar información geológica correspondiente (litología, estructuras,
formaciones, etc.).
1
Nuestro propósito en general al enfocarnos en este tipo de método de exploración
es obtener los datos necesarios para procesarlos de acuerdo con lo que nos
incumbe, para así generar la imagen que será posteriormente analizada y
finalmente interpretada, basándose naturalmente en la teoría respectiva existente.
En esta ocasión, se procesarán datos adquiridos de la línea TRV-434 en la costa
de Nueva Zelanda por NZOG, que datan del año 1986, y, que es una de las
provincias petroleras más abundantes de Nueva Zelanda. Esta línea tiene una
longitud de 22.7 kilómetros.
Se estará trabajando en el tiro de control 300.
TIPO DE FUENTE Pistola de Aire
PROFUNDIDAD DE ARRASTRE DE LA FUENTE 6 m
INTERVALO DE TIROS 25m
TIPO DE RECEPTORES Cable STREAMER
INTERVALO DE GRUPO 120m
PROFUNDIDAD DE ARRASTRE DE LOS RECEPTORES 13m
OFFSET CERCANO 258m
OFFSET LEJANO 3233m
SISTEMA DE GRABACIÓN DFS-5
FORMATO DE DATOS SEGD Cartuchos 3480
INTERVALO DE MUESTREO 2ms
LONGITUD DE GRABACIÓN 6000 ms
Figura 2. Tabla de información general de la adquisición de datos sísmicos
GEOMETRÍA
La geometría es una etapa verdaderamente importante y casi estrictamente
necesario para comenzar el procesamiento, donde, dependerá de las necesidades
de los datos. La geometría otorga la posición real a cada una de las trazas, dicha
ubicación no estará establecida en la fuente o en algún receptor, sino en el punto
medio común (CMP. Cualquier visualización de una colección de uno o más trazas
sísmicas se denomina sismograma. Una colección de tales trazas que representan
las respuestas de una serie de detectores a la energía de un disparo se denomina
colección de disparos. Una colección de trazas relacionadas con la respuesta
sísmica en un punto medio de la superficie se denomina colección de punto medio
2
común (CMP).) La recolección de las trazas sísmicas para cada CMP y su
transformación a un componente de la imagen presentada como una sección
sísmica es la tarea principal del procesamiento de reflexión sísmica.
Estas coordenadas (la ubicación del CMP) se asignarán a cada traza. Es
importante establecer la geometría en dentro de las primeras etapas del
procesamiento, ya que de esa manera identificaremos a las trazas por sus
coordenadas, y, como resultado de esto, se pueden agrupar los datos para el
proceso en general.
Figura 3. Tiro 300 sin agrupar
3
Figura 4. Tiro 300 agrupado por receptor común
GANANCIA
La ganancia es una función de tiempo escalada, en la cual, el escalamiento de la
función está basado en un criterio deseado. La señal al propagarse por el interior
de la tierra sufre una atenuación de las altas frecuencias como se ha mencionado.
Existen diferentes tipos de ganancia que se pueden aplicar. Una función de
ganancia g (t) consiste en una serie de escalares que se multiplican con las
amplitudes y se pueden derivar de los datos y son Para ello existen tres tipos de
ganancia: control de ganancia programada (programed gain control, PGC), el
control de ganáncia automática (AGC, automatic gain control) con amplitud
cuadrática media (root mean square RMS). Esta última es el tipo de ganancia más
comúnmente usado. Las AGC, PGC, etc., utilizan una medida estadística de la
amplitud dentro de una ventana. Esta medida estadística puede ser, el RMS, la
media, la mediana, etc. El cálculo de la medida estadística se realiza dentro de la
ventana, la ventana recorre toda la traza a intervalos de una muestra, generando
una serie numérica de la misma longitud que la traza. Esta serie se multiplica por
la serie de amplitudes de la traza. Si el tamaño de las ventanas es muy pequeño
puede causar una pérdida significante del carácter de la señal. En algunos casos,
se tiene, que las características de la amplitud varían de traza en traza. Esto
puede ser debido a diferencias de amplificación de los sensores, diferencias de
acoplamiento entre el sensor (geófono) y el suelo, etc. Para compensar, corregir
estas diferencias se aplica un balanceado entre trazas. Esto significa que todas las
trazas deben tener un nivel de amplitud RMS similar.
El tamaño de las ventanas es sumamente importante para poder visualizar de
mejor manera los eventos que hay en cada tiro, ya que las amplitudes que hay en
cada sismograma no es el mismo, las ventas adecuadas para información sísmica
van de 250- 500.
Figura 5. Tiro 300 sin ganancia
4
Figura 5. Tiro 300, aplicando la ganancia AGC, con la ventana más alta 250ms. Se observa que quita intensidad y
la energía se empieza a ver constante.
Figura 6. Tiro 300, aplicando la ganancia AGC, en la ventana más pequeña de 500ms. Vemos un gran cambio en
la visualización de los eventos a comparación de la ventana más grande , en esta observamos la ganancia se
define un poco mejor.
5
CORRECCIÓN POR DIVERGENCIA ESFÉRICA
El objetivo de esta corrección es reconstruir las amplitudes debido a la absorción
de los materiales y al decaimiento del frente de ondas con la distancia.
Para compensar esta pérdida de energía los datos se deben de corregir por este
efecto. La corrección es proporcional al tiempo (distancia recorrida). Esta
corrección se denomina divergencia esférica. Simplificando la fuente como un
punto generador de ondas esféricas, en un medio homogéneo la densidad de
energía decae de la forma 1/r2, donde r es el radio del frente de ondas. La
amplitud de la onda es proporcional la energía, la cual decae 1/r. En la práctica,
por lo general la velocidad se incrementa con la profundidad.
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 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. Aplicando la corrección por divergencia podemos observar
que las formas de onda que están a mayor profundidad se pueden distinguir de
una mejor manera, se corrigen las pérdidas de la energía, con respecto al tiempo
de propagación de la energía, en donde los eventos lineales a mayor profundidad
son ruidos. El objetivo de esta corrección es reconstruir las amplitudes debido a la
absorción de los materiales y al decaimiento del frente de ondas con la distancia.
Figura 7. Análisis de amplitud sin corrección por divergencia esférica en el tiro 300
6
Figura 8. Análisis de amplitud aplicada la corrección por divergencia esférica en el tiro 300,en donde se conserva
constante la energía en un rango de un solo color, a diferencia de la que no tiene convergencia esférica.
FILTRO PB (Pasa Bandas /Band Pass)
Un tipo relativamente simple de filtrado de frecuencia es el filtrado de paso
de banda. Se pasan todas las señales por encima de una frecuencia especificada,
o se pasan todas las frecuencias por debajo de otra frecuencia especificada, o se
producen ambas. Dichos filtros se designan generalmente como de paso alto,
paso bajo o paso de banda. Al diseñar filtros digitales que operan en el dominio
del tiempo, expresamos la característica de filtro deseada en función de la
frecuencia y luego calculamos la transformada de Fourier de esta función para
configurar un operador de filtro en función del tiempo.
El filtrado de paso de banda está diseñado para pasar la señal y rechazar el ruido.
Los escaneos de filtros se generan a partir de los datos, en los que se aplican
diferentes filtros estrechos de paso de banda y se grafican los resultados. El
geofísico diseña el filtro final para pasar las frecuencias que contienen energía
coherente (reflejos) y rechazar aquellas frecuencias que contienen principalmente
ruido y sin reflejos aparentes.
Los filtros que varían en el tiempo a menudo se aplican a los datos para pasar
frecuencias más altas a poca profundidad en la sección y frecuencias más bajas a
mayores profundidades.
7
Figura 9. Ejemplo de escaneos de filtro (derecha) en datos de campo
sin procesar (izquierda). El paso de banda del filtro utilizado para cada
panel está anotado en la cabecera de cada panel. Las señales de baja
frecuencia son de al menos 10 Hz (y probablemente más bajas). En este
caso, el corte alto del filtro debería variar en el tiempo: un valor superior
a 60-70 Hz podría cortarse después de 1,5 s; más alto de 50-60 Hz
podría cortarse después de 1,5 s; más alto de 50-60 Hz podría cortarse
Filtro pasa bajas: Se eliminan altas frecuencias, es decir, dejará pasar todas
las frecuencias que sean menores a Fc (frecuencia de corte) las cuales se
multiplicarán por 1, y las que son mayores que Fc se multiplicarán por 0 por
lo que se anularán.
Filtro pasa altas: Se eliminan bajas frecuencias, es decir que dejará pasar
todas las frecuencias que sean mayores a Fc las cuales se multiplicarán por
1, y las que son menores que Fc se multiplicarán por 0, lo que hará que
éstas se anulen. Contrario a filtro pasa bajas.
Filtro pasa banda: Eliminan altas y bajas frecuencias, es decir, al definirse
un intervalo de frecuencias, una mínima Fm y una máxima FM, las
frecuencias que se encuentran dentro de éste intervalo [Fm,FM] se
multiplicarán por 1, y las frecuencias que sean menores a Fm y mayores a
FM se multiplicarán por cero, por lo que se anularán. Éste es el filtro más
utilizado en procesamiento. Generalmente las frecuencias útiles se
encuentran en el intervalo de 10 Hz a 60 Hz
Filtro de rechazo: Contrario a filtro pasa bandas. Dado un intervalo de
frecuencias Fm y FM, éste filtro multiplicará por cero todas las frecuencias
contenidas en este intervalo [Fm,FM], y las restantes serán multiplicadas
por 1.
Figura10. Vista gráfica de filtro PB
8
Todo levantamiento sísmico se verá afectado por el ruido sísmico el cual puede se
causado por fenómenos del propio método (fuentes internas) como son las
refracciones de la onda, difracciones, ruido de pozo, etc,, así como por
condiciones ambientales o fuentes externas como son la marea, el viento, la
actividad humana, irregularidades en el subsuelo, etc. En consecuencia, se
obtendrá información que no será útil, además de afectar, para la interpretación de
los datos.
Por lo anterior, es importante aplicarle un filtro PB a los datos, generalmente se
utiliza el filtro pasabandas con un rango de frecuencias de 10 Hz a 60 Hz, para
suprimir el ruido contenido en los datos, pues si la relación señal ruido (S/R) es
pequeña entonces se tendrán registros pobres. El objetivo en obtener datos más
limpios.
Figura 11. Se utilizo el flujo 03_swell, el
cual es utilizado porque al ir moviéndose
el barco, el receptor unido al streamer va
rebotando y ésta funge como una fuente,
la cual generará energía que será
registrada y es parte del ruido que se tiene
que eliminar
Figura 12. Espectro de frecuencias
después de haber aplicado el filtro PB
9
Observando las imágenes, para el tiro de control 300:
Antes de aplicar el filtro PB, observamos toda la gama de frecuencias
obtenidas en campo. Se puede ver que los datos son muy ruidosos.
Después de aplicar el filtro PB, que en este caso fue pasabandas
[5,10]Hz, se eliminaron frecuencias por debajo de 5 Hz y las que
estaban por arriba de 10 Hz, a partir del cual obtuvimos datos más
limpios, sin tanto ruido.
DECONVOLUCIÓN
La deconvolución es un proceso que mejora la resolución temporal de los
datos sísmicos por medio de comprimir la ondícula de la fuente. Consiste en la
comprensió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 del medio sin alteración generada por el
impulso de entrada, mejorando la resolución temporal, atenuando o removiendo el
ruido coherente y reverberaciones de la fuente generada.
En este proceso la onda sísmica es comprimida a un pulso (pico) de retraso cero
la buena ejecución de este filtro depende de la longitud y de que la onda de
entrada sea de fase mínima
El coeficiente de reflexión está dado por:
𝜌2𝑣2 − 𝜌1𝑣1
𝐶𝑅 =
Donde: 𝜌2𝑣2 + 𝜌1𝑣1
ρ es la densidad de la capa
v es la velocidad de la capa
Es un algoritmo que se utiliza con el fin de aumentar la resolución temporal de las
reflexiones. Esto se consigue invirtiendo la película básica y convolución con la
con cada traza, el resultado es una compresión de la señal. La deconvolución
puede ser aplicada en diferentes etapas del procedimiento esto es de acuerdo con
las necesidades de los que se procesan la señal.
Selecciona una muestra (Tiro) de los datos de la zona de interés se hacen
pruebas con diferentes valores de longitudes de operador y gap, analizan los
resultados finales y la auto correlación de estos resultados el procesador decide la
longitud del operador y el valor gap a usar basado en las conclusiones de los
análisis.
Longitud de operador: Tamaño de la longitud que se tenga la perturbación.
Gap: separación que habrá entre las diferentes separaciones que pueda tener en
esta misma entrada.
10
Figura 13. Especificaciones de la deconvolución
Figura 14. Traza sísmica con
divergencia esférica, geometría sin
deconvolución
Figura 15. Traza sísmica con
11 divergencia esférica, geometría y
deconvolución.
FILTRO FK
Este filtro opera en el dominio de la frecuencia y el número de onda, es un filtro
espacial que ayuda a eliminar el ruido coherente que se presenta de manera lineal
y con velocidades bajas, como, por ejemplo, la onda generada por el aire y el
Ground roll, en la cual la densidad de energía dentro de un determinado intervalo
de tiempo es contorneado basado a una frecuencia versus número de onda. Se
usa para examinar la dirección y velocidad aparente de las ondas sísmicas y en el
diseño de un filtro de velocidad. Modifica la forma en que yo vea mis datos para
poder analizarlo en otros dominios, se aplica la transformada de Fourier en 2
dimensiones, modificando de un dominio de t (tiempo) a uno de f (frecuencia),
representando información en cuanto al espacio (distancia a número de onda)
analizo mis datos y los aplico, diseño el filtro, aplico el filtro, hago la transformada
inversa de mis dos dimensiones, cuando tengo mis datos en el dominio del tiempo
voy a tener información que no tiene cierta cantidad de información o señal que no
será útil, antes de aplicar un filtro f-k, es necesario tener un buen control de las
amplitudes.
Figura 16. Traza sísmica sin aplicación del filtro FK
12
Figura 17. Traza sísmica con aplicación del filtro FK
EDICIÓN
Esta parte del procesado es muy importante pues los resultados finales van a
depender de las modificaciones que se le hagan a la información sísmica. El
objetivo es obtener datos limpios, teniendo la precaución de no eliminar
información útil. La edición se refiere al control que debemos darle al contenido de
amplitudes.
Función Mute
Línea que se trazará en datos sísmicos paralimitar todo lo que esté fuera o dentro
de ésta.
13
enmudecido:
14
Figura 21. Visualizando, utilizando el archivo 01_mutes, y comparando el tiro de control original y el
enmudecido.
Generalmente, aún aplicando la función MUTE, encontramos datos muy
ruidosos, debido a diversos factores como el ruido ambiental. Por ello, es
necesario suprimir todas las trazas del tiro y esto se hace con el borrado
directo.
Flujo 5 Areal 2-D
En este podemos ver una cuadrícula que nos señala todos los tiros con los
que estamos trabajando, el punto rojo que aparece en el centro indica una
amplitud mayor, de acuerdo con la escala de colores.
15
Figura 22. Se puede observar una amplitud anómala en el tiro 500
Del mismo modo que en la función MUTE, el programa multiplicará por cero
esta traza señalada. Al realizar la operación, resulta lo siguiente:
Figura 23. Después de quitar la amplitud anómala.
Es de suma importancia realizar un buen tratamiento a las amplitudes pues
de lo contario nublará la demás información sísmica y dañará el
sismograma. Además, si no se lleva a cabo de forma correcta esta
aplicación, al aplicar el filtro FK se verá totalmente afectado.
Nota: Se permite quitar del 3% al 5% del volumen total de trazas que se
tienen.
16
DISEÑO DE ONDÍCULA
Se sabe que una ondícula es el pulso unidimensional que considerada una
respuesta básica de un reflector. Se origina como si fuera un paquete de energía
desde la fuente, teniendo un origen específico en tiempo, y retorna a los
receptores como una serie de eventos distribuidos en tiempo y energía. Es
importante considerar que la ondícula se encuentra en función del tiempo, junto a
sus dos espectros, el de amplitud y de fase
Figura 24. Forma de onda con la se perturbo el medio junto el espectro de amplitud y de fase
Para obtener una mejor visualización se retiran los primeros 14 ms, recorriendo la
forma de onda de manera que ahora inicia en su primer arribo. El espectro de
amplitud igual se recorre observe en (Fig 24)
Figura 25. Ondícula, espectro de amplitud y de fase mínima
Obteniendo estas dos ondículas se combinarán es decir aplicaremos un
“MarchFilt” que combinara la ondícula de fase original con la ondícula de fase
mínima, Esta se realizara usando las primeras 100 muestras. (Fig 25).
17
Figura 26. Ondícula original a fase mínima
Para finalizar con el diseño de la ondícula aplicaremos una convolución para
obtener una ondícula de fase mínima, que tiene las características adecuadas
de la señal. (Fig. n)
Figura 27. Ondícula de fase mínima con las características adecuadas de la señal para poder trabajar en ella
Flujo 10
En la primera imagen obtenemos nuestra traza sísmica sin datos es decir sin
ninguna ondícula
Figura 28. Traza sísmica donde aún no se encuentra aplicada la convolución
18
Al aplicar una autocorrelación para comparar la traza sísmica consigo misma así
mismo se nos facilitara detectar los patrones (forma de onda con la cual se
perturbo el medio). Notemos que en la autocorrelación se muestra una sección
oscura, conocida como “máximo primario”, que es la más parecida a la forma de
onda, los que se encuentran en la sección inferior son conocidos como los
máximos secundarios.
Esta forma de onda se va a encontrar modificada, en esta práctica se opta por
encontrar la más larga, que va a ser la longitud de nuestro operador, se nombran
2do máximo secundario tercer máximo secundario y así sucesivamente, esta
forma de onda tan larga que en nuestro caso tenemos una de 500 ms que no es
posible usarla, ya que se pueden producir alteraciones en las condiciones
elásticas, por lo mismo se debe buscar una más pequeña.
Figura 29. Traza sísmica donde se encuentra aplicada la convolución
19
Figura 30. Autocorrelación en la traza original
Flujo 11
Se modificarán los parámetros de la deconvolución donde se colocaron los
parámetros obtenidos de longitud de operador y del gap. Esto se realizó a la
opción DECONV dentro del flujo 11, para con ello colocar las siguientes longitudes
de operador con los valores de 50, 100, 150, 200, 300 y 600 (ms). El gap se deja
con un valor de 24ms por el momento para cada una de las seis líneas.
Recordemos que el objetivo principal de la deconvolución es hacer que las formas
de onda se compriman es decir que se vean mucho mejor definido.
20
Figura 31. Flujo 11 sin aplicar deconvolución con un filtro de 600 ms
Es importante mencionar que al final de este sismograma notamos una onda de
autocorrelación con baja nitidez es importante identificar esto ya que nos indicara
el máximo primario junto al 2do máximo secundario.
Figura 32. Flujo 11 con la autocorrelación
21
Finalmente se decidió en escoger que el mejor filtro es el de 200 ms y un gap de
48 ms eran lo más ideales para visualizar los tiros y sus componentes ya que se
puede observar de una mejor manera los eventos de reflexión si ser muy nítidos y
ni muy marcados.
Figura 33. Flujo 11 con un filtro de 200 ms y un gap 48 ms
APILADO
Al abrir el flujo 12 notamos que ya tiene una corrección NMO este es aplicada a un
conjunto de datos de desplazamiento largo generalmente crea un efecto palo de
hockey dando tiempos de viaje demasiado pequeños asociados a
desplazamientos grandes y causando ensanchamiento de la forma de
22
onda y pérdida de resolución. Es importante considerar que filtro 12 contenía y
previamente una deconvolución de 500 ms, corrección MNO, apilamiento y
autocorrección. De los cuales se modificaron parámetros y otros quedaron sin
modificar
Figura 34. Apilamiento del flujo 12
CONCLUSIÓN
El procesamiento de datos sísmico llevado a cabo en el presente trabajo surge de
una metodología para llegar al resultado " apilamiento" realizando un estudio de
los datos, cambiando los parámetros para una mejor visualización dando
resultados acertados. Cada uno de los procesos aplicados son muy importantes
para seguir avanzando al siguiente paso, desde el principio empezando por aplicar
la geometría para tener un orden congruente, la convolución donde conocemos
propiedades de la señal. Es importante considerar el preprocesado de datos
además de considerar el uso correcto del SOFTWARE CLARITAS para el
procesado de datos debido a que el software prepara los sismogramas de una
manera más visual, optima y automatizada para después llegar al pos-procesado
de datos.
23
Es importante mencionar que trabajar en modalidad en línea esta materia ha sido
un reto debido al trabajo en equipo, pero nos pudimos adaptar de buena manera
gracias a los videos de las clases pregrabadas donde podíamos consultar los
pasos a seguir de una manera autodidacta se complementó con videollamadas y
secciones para trabajar en conjunto de una manera más eficiente. Respecto a la
materia, nos pareció una materia muy interesante, debido a la complejidad de
cosas que se ven y aprenden debido a que estos temas son muy importantes para
un buen desarrollo profesional, siendo está una materia base para la carrera,
agradecerle por su enseñanza y empatía. Cabe mencionar que se pudo entender
con claridad cada uno de los conceptos, pasos y prácticas y así el equipo pudo
desarrollar las habilidades y cumplir con los objetivos para el curso de
preprocesado sísmico.
24
25