Universidad Simón Bolívar Departamento de Ciencias de la Tierra Procesamiento Sísmico MODELADO Y PROCESAMIENTO SÍSMICO DE UNA ESTRUCTURA DE CORRIMIENTO Y ESTRUCTURA TIPO FLOR POSITIVA UTILIZANDO EL PAQUETE SEISMIC UNIX Profesor: Realizado por: Ilich García Ada Zamora #06-40493 Carla Barbara #06-39202 Sartenejas, marzo de 2011
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Universidad Simón Bolívar
Departamento de Ciencias de la Tierra
Procesamiento Sísmico
MODELADO Y PROCESAMIENTO SÍSMICO DE UNA
ESTRUCTURA DE CORRIMIENTO Y ESTRUCTURA TIPO
FLOR POSITIVA UTILIZANDO EL
PAQUETE SEISMIC UNIX
Profesor: Realizado por:
Ilich García Ada Zamora #06-40493
Carla Barbara #06-39202
Sartenejas, marzo de 2011
RESUMEN
El siguiente trabajo trata sobre el modelado y procesamiento básico de una data sísmica no
real de una estructura compleja utilizando el programa Seismic Unix (SU), para evaluar la
capacidad y utilidad de dicho programa para modelar y aplicar una secuencia de procesamiento
básica a una zona geológica de interés.
El modelo estructural elegido fue una zona que presenta una estructura tipo flor positiva
junto a una estructura de corrimiento, dicho modelo se graficó en papel milimetrado para calcular
coordenadas ficticias para su modelado sísmico en el programa SU. Luego se procedió a realizar el
trazado de rayos, el diseño sintético de adquisición sísmica, el análisis de velocidades, la
corrección NMO, y la migración post-apilamiento, todo esto utilizando dicho programa.
En el informe se señalan varias definiciones para la mayor comprensión de proyecto y una
explicación detallada del programa a utilizar y los pasos realizados para la obtención del los
resultados.
TABLA DE CONTENIDO
- Introducción
- Marco Teórico
- Marco Metodológico
Generación del modelo geológico del subsuelo
Generación del trazado de rayos
Modelado de Propagación de ondas acústicas
Análisis de Velocidad, Corrección NMO y Sección Apilada
Mute y Migración de la Data Sísmica
- Resultados y Análisis
Modelo Geológico en 2D generado con Seismic Unix
Trazado de Rayos y Propagación del frente de ondas a través del modelo
Trazado de Rayos y Frentes de Ondas
Sismogramas Sintéticos
Propagación de ondas acústicas en diferencias finitas para el modelo generado
Shots Gathers
Ordenamiento por CMP
Análisis de Velocidad, Corrección NMO y Sección Apilada
Mute y Migración de la Data Sísmica
Conclusiones
Referencias
Apéndice 1: 2_perfil_tri.sh
Apéndice 2: 6_varios_shots_rayos.ah
Apéndice 3: 8_varios_shots_acusticas.sh
Apéndice 4: 9_nmo.sh
Apéndice 5: 10_miga.sh
TABLA DE IMÁGENES
• Figura 1: Pasos básicos del procesamiento sísmico (Ilich García, 2008)
Figura2: Estructuras tipo flor negativa (A) y flor positiva (B).
Figura3: Modelo geológico asignado
Figura4: Merfil de Lentitud Cuadrática
Figura5: Perfil de Velocidad
Figura6: Modelo Triangular
Figura7: Trazado de Rayos en 2km para la interfaz 2
Figura8: Frente de Ondas en 2km para la interfaz 2
Figura9: Trazado de Rayos en 2km para la interfaz 3
Figura10: Frente de Ondas en 2km para la interfaz 3
Figura11: Trazado de Rayos en 2km para la interfaz 6
Figura12: Frente de Ondas en 2km para la interfaz 6
Figura13: Trazado de Rayos en 3km para la interfaz 3
Figura14: Frente de Ondas en 3km para la interfaz 3
Figura15: Trazado de Rayos en 3km para la interfaz 4
Figura16: Frente de Ondas en 3km para la interfaz 4
Figura17: Trazado de Rayos en 4km para la interfaz 8
Figura18: Frente de Ondas en 4km para la interfaz 8
Figura19: Trazado de Rayos en 2km para la interfaz 12
Figura20: Frente de Ondas en 2km para la interfaz 12
Figura21: Sismograma sintético en 2km para la interfaz 2
Figura22: Sismograma sintético en 2km para la interfaz 3
Figura23: Sismograma sintético en 2km para la interfaz 6
Figura24: Sismograma sintético en 3km para la interfaz 3
Figura25: Sismograma sintético en 3km para la interfaz 4
Figura26: Sismograma sintético en 4km para la interfaz 8
Figura27: Sismograma sintético en 2km para la interfaz 12
Figura28: Secuencia de imágenes de la película generada
Figura29: Shots Gathers correspondientes a un disparo en posición 2000m
Figura30: Shots Gathers correspondientes a un disparo en posición 3200m
Figura31: Shots Gathers correspondientes a un disparo en posición 4000m
Figura32: Shots Gathers correspondientes a un disparo en posición 5200m
Figura33: CMP Gathers
Figura35: CMP gather antes de la corrección (izquierda) y CMP gather después de la
corrección NMO
Figura36: Mapa de Semblanza para CMP=4000m
Figura37: CMP gather antes de la corrección (izquierda) y CMP gather después de la
corrección NMO
Figura38: Mapa de Semblanza para CMP=5000m
Figura39: CMP gather antes de la corrección (izquierda) y CMP gather después de la
corrección NMO
Figura 40: Sección apilada
Figura 41: Sección apilada sin migrar (izquierda) y sección apilada migrada (derecha)
Figura 41: Sección apilada sin migrar (izquierda) y sección apilada migrada (derecha)
INTRODUCCIÓN
Como estudiantes de las Geociencias, sabemos que lo que se encuentra en el subsuelo es un
misterio que continuamente tratamos de revelar. Por esta razón, se han ideado una gran cantidad
de métodos para tratar de obtener información del subsuelo.
El paquete Seismic Unix, constituye una herramienta valiosa para poder realizar simulaciones
de modelos estructurales, así como de los resultados obtenidos en caso de hacer una adquisición
sísmica. Y además de los resultados, ofrece la opción de poder procesarlos, todo ello para tener
interpretaciones más acertadas de nuestra eterna incógnita.
En este trabajo, se trabajará con SU con el propósito de analizar los resultados obtenidos y
ampliar nuestros conocimientos sobre el procesamiento sísmico.
MARCO TEÓRICO
El procesamiento sísmico tiene como objetivo principal transformar las observaciones sísmicas
adquiridas en campo en secciones o volúmenes de datos para poder realizar interpretaciones que
conlleven a la elaboración de un modelo del subsuelo y ver si es rentable la exploración y/o
producción de hidrocarburos. Durante el procesamiento de la data se obtienen parámetros
importantes como velocidades, amplitudes y fase.
Figura 1: Pasos básicos del procesamiento sísmico (Ilich García, 2008)
Para la realización de este trabajo se utilizó el programa Seismic Unix, el cual se operó en el
sistema operativo Ubuntu
Seismic Unix (SU), es un paquete de software libre y de fuente abierta, que realiza modelado y
procesamiento de una data sísmica. Fue desarrollado en el Centro de onda Fenómenos (CWP) en
la Escuela de Minas de Colorado (CSM) en 1987 por los científicos Jack K. Cohen y Shuki Ronen,
entre varios más. Seismic Unix fue programado en lenguaje C y trabaja con base en sistemas
operativo UNIX.
SU contiene herramientas para el modelado, procesamiento, transformadas de Fourier,
filtrado de señales, transformadas de ondículas, entre otras funciones.
Al utilizar la herramienta de modelado, se pueden crear perfiles de propiedades del subsuelo,
formando: modelos uniformemente muestreados y modelos triangulares. También, SU genera la
data sísmica mediante el trazado de rayos y la propagación de ondas acústicas.
Con la herramienta de procesamiento sísmico podemos visualizar y ajustar el header, hacer
ordenamiento por CMP, análisis de velocidad, corrección NMO, apilamiento y migración.
Para poder utilizar este programa y correr los códigos para modelar, es necesario tener claro
ciertos conceptos enunciados a continuación:
- Rayo: Cada una de las líneas, generalmente rectas, que parten del punto en que se
produce una determinada forma de energía y señalan la dirección en que esta se propaga.
Línea de luz que procede de un cuerpo luminoso
- Frente de Ondas: Se denomina frente de onda al lugar geométrico en que los puntos del
medio son alcanzados en un mismo instante por una determinada onda. Dada una onda
propagándose en el espacio o sobre una superficie, los frentes de onda pueden
visualizarse como superficies o líneas que se desplazan a lo largo del tiempo alejándose de
la fuente sin tocarse.
- Shot Gather: Conjunto de trazas que representan la respuesta de una serie de detectores
desde un mismo disparo.
- Common Midpoint (CMP): Es el punto situado en el centro de los diferentes pares fuente-
receptor que corresponden a reflexiones de un mismo punto del subsuelo
- Offset: Es la distancia entre la fuente y el receptor
- Velocidad NMO: Es la velocidad empleada para la corrección de “Normal Move Out” cuando la distancia fuente receptor es pequeña.
- Migración: Es el proceso de reubicación de lo estratos a su posición original.
Modelo geológico generado: Estructuras en flor
La presencia de inflexiones y saltos a lo largo de una falla promueve la formación de
estructuras denominadas duplexes). Las rocas que se encuentran dentro de la zona de influencia
de una inflexión o salto se pueden fallar progresivamente, generando sistemas imbricados en
relación a la falla principal.
El análisis del proceso de formación de duplexes puede realizarse a través de dos
mecanismos: Deformación planar y Deformación no planar.
La formación de duplexes por deformación no planar se encuentra en la compensación
volumétrica por hundimientos y levantamientos. Así en duplexes extensionales se formaran
hundimientos, y en los compresionales, levantamientos. Unas de las estructuras más notables
derivadas de la deformación no planar son las denominadas flores positivas (positive flowers) y
flores negativas (negative flowers) (Woodcock y Fischer, 1986), que a su vez se corresponden con
los casos de duplex compresional y extensional, respectivamente. La geometría interior de las
fallas secundarias en un duplex muestra que estas convergen en profundidad hacia una zona de
falla única. En el caso de una estructura en flor negativa (fig.A), las fallas secundarias dentro del
duplex tendrán un comportamiento mixto, transcurrente-normal, por lo cual podemos decir que la
zona interior está en extensión y sometida a hundimiento. Por el contrario, en una estructura tipo
flor positiva (fig.B), las fallas interiores se comportarán como transcurrentes-inversas, por lo cual
podemos decir que la zona está en compresión y por lo tanto sometida a levantamiento. Esta
estructura tiende a crecer con el aumento del desplazamiento experimentado por las fallas que la
confeccionan.
Figura 2: Estructuras tipo flor negativa (A) y flor positiva (B). Note la subsidencia y alzamiento respectivo en A y B, así mismo como el carácter mixto (transcurrente-normal y transcurrente-inverso) del fallamiento asociado.
MARCO METODOLOGICO
Generación del modelo geológico del subsuelo
El modelo asignado para realizarle el procesamiento fue el de una estructura debido a inversión tectónica positiva, generando una estructura en flor con dos fallas inversas en corrimiento.
Las dimensiones tomadas para el modelo son de 8km de longitud y 5 km de profundidad, tomando como referencia el tamaño real de la estructura. Para poder graficarlo, lo dibujamos sobre papel milimetrado y obtuvimos las coordenadas de cada segmento.
Figura3: modelo geológico asignado.
Ya que el modelo presenta 3 estratos, se le asignaron velocidades a cada uno. Esto se hizo
tomando como referencia velocidades características de litologías comunes. La primera capa tiene
una velocidad de 1500m/s, la segunda una de 2000m/s y la tercera 3000m/s.
En Seismic Unix, el comando Trimodel es utilizado para generar la representación del subsuelo
con las coordenadas y velocidades dadas.
Una vez obtenidas las coordenadas y velocidades del modelo, estas se insertaron en el código
22_perfil_tri.sh. Además, en el código se colocaron los valores de lentitud cuadrática, la cual viene
dada por la relación L=1/V2
Generación del trazado de rayos
Para la generación del trazado de rayos, simulando un arreglo Split spread para varios disparos, se utilizo el comando TRIRAY para el modelo construido con el comando TRIMODEL.
Luego de introducir los valores de nuestro modelo, se introdujeron los siguientes parámetros: interf= 13 Número de interfaces, incluyendo los bordes superior e inferior nangle= 20 Número de ángulos de despegue fangle= -60 Primer ángulo de despegue langle= 60 Último ángulo de despegue nxz= 500 Número de (x,z) en el archivo rayfile fpeak= 10 Frecuencia pico de la ondícula Ricker en Hz offset_max=1 Offset máximo (debe ser un entero) ID=1 Intervalo entre puntos de disparos (debe ser un entero) n1=1000 Número de muestras en z d1=0.005 Intervalo de muestreo en z
Modelado de Propagación de ondas acústicas
Para el modelado de propagación de ondas acústicas en diferencias finitas de nuestro modelo, se utilizó el comando SUFDMOD2, simulando varios disparos para un arreglo split spread.
En este código, los parámetros introducidos debían estar en metros, por lo que transformamos los valores del modelo y de las velocidades a esta unidad
Estos parámetros fueron:
offset_max=2000 Offset máximo
IR=20 Intervalo entre grupos de receptores
fmax= 20 Frecuencia máxima de la ondícula fuente
tmax= 5 Tiempo máximo de grabación
zs= 0 Coordenada z de la fuente
hsz= 0 Coordenada z de la linea horizontal del sismograma (profundidad
de los receptores)
ID=400 Intervalo entre puntos de disparos
Longitud=8000 Longitud del modelo en metros
Profundidad=5000 Profundidad del modelo en metros
fz= 0 Coordenada z de la primera muestra
min_shot=1500 Valor mínimo del shot que se desea graficar
max_shot=6500 Valor máximo del shot que se desea graficar
min_cmp=1500 Valor mínimo del cmp que se desea graficar
max_cmp=6500 Valor máximo del cmp que se desea graficar
clip= 0.49 Valor máximo de la amplitud de la traza en la gráfica de los shots
gathers formato wigle trace display
Análisis de Velocidad, Corrección NMO y Sección Apilada
Para realizar el análisis de velocidades y la corrección NMO se utilizó el archivo CMP arrojado
por el modelado de la propagación de ondas acústicas, para utilizar los CMP Gathers
correspondientes a nuestro modelo. Adicionalmente, se proporcionaron ciertos parámetros que
permitieran generar los mapas de semblanza en función de los cuales se seleccionan las
velocidades para la corrección NMO.
ND= 10 Número de disparos de la adquisición sísmica
cdpmin= 2000 Primer CMP
cdpmax= 6000 Último CMP
dcdp= 500 Intervalo entre CMP
nv= 500 Número de velocidades
dv= 10 Intervalo entre las velocidades
fv= 500 Primera velocidad
vf = fv + ((nv*dv) - dv) = 10990
Las velocidades de corrección para cada CMP Gather se seleccionaban a partir del mapa de
semblanza correspondiente, oprimiendo la tecla “s” con el cursor encima de los puntos de mayor
semblanza que estuvieran dentro del rango de tiempo abarcado por las hipérbolas de las
reflexiones.
Una vez seleccionadas todas las velocidades y horizontalizados los CMP Gathers, se genera
una sección apilada de los mismos.
Mute y Migración de la Data Sísmica
La migración y el mute se realizaron para la sección apilada generada previamente. Para realizar el mute, se selecciona el área escogida pulsando la tecla “s” y colocando el cursor
en dos puntos que unen la línea por encima de la cual se desea eliminar las amplitudes.
En el caso de la migración deben ser escogidos los parámetros siguientes: v = velocidad de migración (debe corresponder a un aproximado de la velocidad de la capa a migrar, aunque no necesariamente debe ser igual)
t = tiempo para el cual se realiza la migración.
RESULTADOS Y ANALISIS
Modelo Geológico en 2D generado con Seismic Unix
El modelo tectónico a generar corresponde a una estructura de flor de dos capas que muestra
inversión tectónica.
El modelado se hizo utilizando Seismic Unix (SU) a través del comando Trimodel. Luego de correr el código 2_perfil_tri.sh, se obtuvieron los perfiles de lentitud cuadrática y velocidad:
Figura4: Perfil de Lentitud Cuadrática
Figura5: Perfil de Velocidad
Se obtuvo también, el modelo triangular:
Figura6: Modelo Triangular
En este modelo, se pueden apreciar los triángulos generados para crear el perfil por medio del método de Triangulación de Delaunay, que corresponden a las líneas grises que atraviesan todo el modelo.
Se aprecia la diferencia en las velocidades de cada capa, con velocidades distintas en cada una.
Las diferentes velocidades se presentan con los colores blanco, gris y negro.
Trazado de Rayos y Propagación del frente de ondas a través del modelo
Una vez creado el modelo triangular, para generar el trazado de rayos es necesario convertirlo en un modelo uniforme. Esto se logra utilizando el código de 6_varios_shots_rayos.sh.
En el modelo generado no se presentan las líneas que definen el perfil, sólo se aprecian las
interfaces de los tres estratos que presentan velocidades distintas. Trazado de Rayos y Frentes de Ondas Para simular el trazado de rayos sobre el modelo creado anteriormente, se asignaron ciertos
parámetros de adquisición, considerando un arreglo Split Spread. Con el propósito de poder abarcar gran parte del perfil, se escogió el rango de ángulos para los
rayos entre -60° y 60° y un offset máximo de 1km. A continuación, se presentan algunas imágenes de los rayos y los frentes de ondas generados
en las interfaces del modelo.
Figura7: Trazado de Rayos en 2km para la interfaz 2
Figura8: Frente de Ondas en 2km para la interfaz 2
Figura9: Trazado de Rayos en 2km para la interfaz 3
Figura10: Frente de Ondas en 2km para la interfaz 3
Figura11: Trazado de Rayos en 2km para la interfaz 6
Figura12: Frente de Ondas en 2km para la interfaz 6
Figura13: Trazado de Rayos en 3km para la interfaz 3
Figura14: Frente de Ondas en 3km para la interfaz 3
Figura15: Trazado de Rayos en 3km para la interfaz 4
Figura16: Frente de Ondas en 3km para la interfaz 4
Figura17: Trazado de Rayos en 4km para la interfaz 8
Figura18: Frente de Ondas en 4km para la interfaz 8
Figura19: Trazado de Rayos en 2km para la interfaz 12
Figura20: Frente de Ondas en 2km para la interfaz 12
En las imágenes anteriores se puede observar la reflexión y la refracción del rayo al conseguirse
con las distintas interfaces, así como también los frentes de onda generados.
Sismogramas Sintéticos
En la práctica, las ondas reflejadas al llegar a la superficie son detectadas por los receptores y con
esta información se puede generar un sismograma sintético. Experimentalmente, para este modelo,
el sismograma sintético se puede obtener mediante las simulaciones en SU.
Los sismogramas generados a través de estas simulaciones representan sismogramas del tipo
Beam Gaussiano, por lo que es necesario definir una frecuencia pico (10Hz), así como el intervalo de
muestreo (0.005s) y el número de muestras (1000). Esto con el propósito de definir la longitud del eje
vertical del sismograma.
Figura21: Sismograma sintético en 2km para la interfaz 2
Figura22: Sismograma sintético en 2km para la interfaz 3
Figura23: Sismograma sintético en 2km para la interfaz 6
Figura24: Sismograma sintético en 3km para la interfaz 3
Figura25: Sismograma sintético en 3km para la interfaz 4
Figura26: Sismograma sintético en 4km para la interfaz 8
Figura27: Sismograma sintético en 2km para la interfaz 12
Podemos observar que para algunas interfaces, la información recibida en lo geófonos es casi nula,
como es el caso de la figura XX, esto se debe a que las ondas reflejadas no llegan hasta los geófonos
en la superficie, o perdieron mucha energía en el camino y su llegada es muy débil para ser
detectada.
Propagación de ondas acústicas en diferencias finitas para el modelo generado
Más allá del trazado de rayos y de la propagación de los frentes de onda, SU permite crear una
propagación de ondas acústicas y de esta manera generar tanto los shot gathers correspondientes,
como realizar el ordenamiento por CMP gathers y crear un película que muestra la trayectoria de
la onda desde su generación hasta sus reflexiones y refracciones. Para lograr esto, se establecieron
parámetros de adquisición específicos.
En figura siguiente podemos observar 8 snapshots de la película que simula la propagación de
ondas acústicas en el medio, para el disparo ubicado en 4000m. En ella se puede apreciar el punto
donde se genera la onda y como esta se va propagando y reflejando en cada una de las interfaces.
Figura28: Secuencia de imágenes de la película generada
Shots Gathers
Los Shot Gathers obtenidos para esta simulación fueron generados en formato Gray en en
formato Wigle, ello para ver de dos formas los resultados obtenidos. A continuación se presentan
los Shots Gathers para los disparos a 2000m, 3200m, 4000m y 5200m
En esta figura se pueden observar en los Shots gathers la onda directa, (aprox entre 0 y 1.2
segundos) y una primera reflexión al conseguir una interfaz (a los 2 segundos)
Figura29: Shots Gathers correspondientes a un disparo en posición 2000m
En este caso, se puede observar la onda directa (entre 0 y 1.5seg), una primera reflexión
aproximadamente a los 2seg y una tenue segunda reflexión a los 3 seg.
Figura30: Shots Gathers correspondientes a un disparo en posición 3200m
Aquí se puede observar nuevamente la onda directa (0 a 1.5seg), una primera reflexión (1.7 a 2.2
seg), una segunda reflexión (3seg.) y hasta una tercera reflexión (4 a 5seg).
Figura31: Shots Gathers correspondientes a un disparo en posición 4000m
En este caso se observa igualmente la onda directa (0 a 1.5seg), se aprecia una primera
reflexión aproximadamente entre 1.6 y 2 seg., una segunda reflexión a los 3 seg. y una tercera
reflexión a los 4 seg.
Figura32: Shots Gathers correspondientes a un disparo en posición 5200m
En todos los Shot Gathers mostrados, se observa que los que presentan mejor resolución
son los que se encuentran en formato Gray, por los que resulta más conveniente estudiar estos
shots gathers, ya que podemos obtener una interpretación más precisa.
Ordenamiento por CMP
Luego de generar los shots gathers, estos se ordenan en el dominio de los CMP gathers
para poder realizar el procesamiento de los mismos. En este ordenamiento de pueden apreciar
todas las trazas generadas de acuerdo a su CMP.
Figura33: CMP Gathers
Análisis de Velocidad, Corrección NMO y Sección Apliada
Luego de ordenar las trazas de acuerdo a los CMP, el siguiente paso en el procesamiento es el
análisis de velocidades y la corrección NMO. Con SU, es posible realizar estas secuencias de
procesamiento para la simulación realizada del modelo de inversión tectónica.
Para realizar el análisis de velocidades, podemos usar el espectro de velocidades. El atributo
utilizado para generar estos espectros en este caso, es la semblanza.
La semblanza se define como una medida de la coherencia de las velocidades. En este caso, se
puede representar con escala de colores. El color azul, representa semblanza 1, lo cual indica que
el valor de velocidad es muy coherente y por lo tanto confiable, a medida que la semblanza toma
colores más cálidos, va perdiendo coherencia.
A continuación se presentan algunos mapas de semblanza y las correcciones NMO hechas
Figura34: Mapa de Semblanza para CMP=3000m
Figura35: CMP gather antes de la corrección (izquierda) y CMP gather después de la corrección NMO
Figura36: Mapa de Semblanza para CMP=4000m
Figura37: CMP gather antes de la corrección (izquierda) y CMP gather después de la corrección NMO
Figura38: Mapa de Semblanza para CMP=5000m
Figura39: CMP gather antes de la corrección (izquierda) y CMP gather después de la corrección NMO
El hecho de que los eventos se hayan podido horizontalizar indica que la sección ha sido bien
corregida. En las figuras anteriores, se observa la diferencia entre las secciones corregidas y sin
corregir.
Esto indica que el método de corrección NMO basado en mapas de semblanza es realmente
útil a la hora de procesar datos sísmicos. Esto se evidencia con los ejemplos dados.
Al aplicar esta corrección a todos los CMP se genera una sección apilada de la data sísmica.
Esta sección es:
Figura 40: Sección apilada
Mute y Migración de la Data Sísmica
Estos representan los últimos pasos en el procesamiento sísmico. El mute se realiza con el
propósito de eliminar amplitudes no deseadas en la sección, aquellas que no aportan información
relevante, como es el caso se las primeras llegadas.
La migración se realiza para mover los eventos buzantes a sus posiciones correctas.
La migración se obtuvo tomando en cuenta la velocidad y el tiempo. Esta sección comprende
a los resultados obtenidos al realizar la migración con un tiempo de 1.5seg y una velocidad de
1000m/s.
Figura 41: Sección apilada sin migrar (izquierda) y sección apilada migrada (derecha)
Migración para un tiempo de 1.5 seg, y una velocidad de 2000m/s:
Figura 42: Sección apilada sin migrar (izquierda) y sección apilada migrada (derecha)
CONCLUSIONES
Mediante el uso de las herramientas de SU, se pueden generar modelos geológicos de
estructuras y realizar simulaciones con los mismos. Sin embargo, al momento de realizar el
modelo, se debe tener el mayor cuidado posible, suministrando al programa las
coordenadas a graficar, ya que un pequeño error en esta etapa, evita que podamos
obtener el modelo respectivo.
Al momento de interpretar los shots gahters generados, las reflexiones se pueden ver
mejor en el shot gather de formato Gray que en el formato Wigle.
Gracias a la utilización de gráfico de espectro de velocidad, se escogieron valores con
mayor coherencia para aplicar la corrección NMO a las trazas, pudiéndose así
horizontalizar los eventos para corregir por NMO.
Los resultados obtenidos fueron certeros. Esto demuestra el alto grado de confiabilidad
que posee esta simulación.
Los códigos utilizados (suministrados por la ayudante docente) para generar las
simulaciones, fueron de gran ayuda para obtener resultados adecuados.
REFERENCIAS
Regueiro, José. Sísmica para Ingenieros y Geólogos. El Método Sísmico de Reflexión.
Universidad Simón Bolívar. (2007).
Judson L. Ahern . Seismic Reflection Method (2007)
Excess Geophysics . Migration in practice (1999)
Pérez, Verónica. Utilización del paquete Seismic Unix para la enseñanza de la sísmica.
Universidad Simón Bolívar. Departamento de Ciencias de la Tierra. (2008).
Vieira, Adelo. Modelado básico de datos sísmicos preapilados y su correspondiente
secuencia de procesamiento utilizando los programas del paquere Seismic Unix.
Universidad Simón Bolívar. Departamento de Ciencias de la Tierra. (2007).
García, Illich. Procesamiento Sísmico. Fundamentos del procesamiento de senales (2009)
Urdaneta, Vanessa. Procesamiento Sísmico. Análisis de velocidad (2011)