secuencia de procesamiento convencional hasta migración pre

Secuencia de Procesamiento convencional hasta Migración Pre – Apilamiento 

en tiempo del Levantamiento Sísmico Arecuna Oeste 04G – 3D. 

Yilmar Solano Vesga. 

UNIVERSIDAD SIMÓN BOLÍVAR 

Coordinación de Ingeniería Geofísica 

Departamento de Ciencias de la Tierra. 

Por: 

Proyecto de Grado 

Presentado ante la Ilustre Universidad Simón Bolívar 

Como requisito Parcial para optar al título de 

Ingeniero Geofísico 

Sartenejas, Septiembre de 2006. 

– .




UNIVERSIDAD SIMÓN BOLÍVAR 

Coordinación de Ingeniería Geofísica 

Departamento de Ciencias de la Tierra. 

Por: 

Proyecto de Grado 

Presentado ante la Ilustre Universidad Simón Bolívar 

Como requisito Parcial para optar al título de 

Ingeniero Geofísico 

Tutor Académico: Ph. D. Milagrosa Aldana 

Tutor Industrial: Ph. D. Horacio Galan. 

Sartenejas, Septiembre de 2006. 

– .



Este trabajo ha sido aprobado en nombre de la Universidad Simón Bolívar por el 

siguiente jurado calificador: 


__________________________________ 

Presidente: Prof. Mario Caicedo 

__________________________________ 

Tutor Académico: Prof. Milagrosa Aldana. 

__________________________________ 

Tutor Industrial: Dr. Horacio Galan.



!" ! # ! $ ! % ! & $ !' ( $ 

$ ( # % $ $ ) ! !" $ * 

Por 

Yilmar Solano Vesga 

Los datos del levantamiento sísmico Arecuna Oeste 04G – 3D fueron procesados con la 

secuencia convencional hasta migración pre – apilamiento en tiempo con el algoritmo de 

Kirchhoff, generando un volumen sísmico apilado con una apropiada resolución de los eventos, 

permitiendo futuras interpretaciones estructurales en área. 

El acondicionamiento previo de los datos se fundamentó en el mejoramiento de la 

relación señal – ruido, atenuando el efecto del ground roll en los registros con el empleo de 

distintas funciones en el dominio f,t. Luego, sobre estos datos se aplicaron las funciones básicas 

(deconvolución, balanceo de amplitud, ganancia por divergencia esférica, etc.) para la 

generación de los datos a migrar. Por otro lado, para el cálculo de las correcciones de estáticas de 

refracción se dispuso de los primeros quiebres seleccionados con anterioridad al inicio de este 

proyecto, cada diez (10) disparos, y los demás seleccionados con la metodología de redes 

neuronales. Se realizaron tres (3) análisis de velocidades con sus respectivas correcciones de 

estáticas residuales y se aplicaron funciones que ecualizaran el espectro de amplitud para 

obtener gathers que generaran una sección migrada con una apropiada resolución. Posterior al 

tercer análisis de velocidades, se realizaron tres iteraciones de migración pre – apilamiento en 

tiempo para seleccionar los residuos de velocidad y refinar el cubo de velocidades, obteniéndose 

a partir de éste una sección migrada con reflectores continuos y coherentes fáciles de interpretar. 

Se comparó la sección migrada post – apilamiento con la pre –apilamiento y se propuso 

seleccionar los primeros quiebres en todos los disparos para mejorar las correcciones de estáticas 

de refracción y obtener una sección con mayor continuidad en los eventos. 

Con las nuevas correcciones estáticas se volvió a realizar los análisis de velocidades y las 

secciones apiladas no mostraron gran diferencia en comparación con las anteriores. Esto sugirió 

que fue apropiada la solución de las primeras correcciones de estáticas de refracción, de la 

selección de los primeros quiebres dada por la metodología de redes neuronales, que se aplicó al 

cubo migrado obtenido. 

Yilmar Solano Vesga Vesga. 

iv



A Dios, por la salud que me brinda y por guiarme de forma grata, paciente, esmerada, 

conveniente y sabia en el camino que me trazo. El cual ha sido el mejor que pudo escoger, 

logrando ofrecerme oportunidades que sabré aprovechar. 

A mis padres, por su incondicional apoyo, fortaleza, cariño, confianza, dedicación, 

consejos y esmero que me han brindado. Por toda la ayuda que me han dado en los momentos 

más difíciles y su enorme colaboración para satisfacer todas mis necesidades; desde buscarme en 

la universidad, acompañarme a donde fuera para cumplir una meta, atender cariñosamente a 

mis amigos en casa, hasta trasnocharse cuando terminaba tarde un trabajo. 

A mi hermana, por soportarme cuando más agobiada y cansada estaba, por ser esa 

persona en quien siempre confiare y contare de forma incondicional, por hacerme reír con sus 

cuentos y anécdotas locas. 

A Teddy, Daniel Rosales, por acompañarme siempre, aguantando las risas en los mejores 

momentos y las lagrimas y quejas cuando estaba preocupada o desmoralizada. Por ayudarme en 

todo lo que le pedía, logrando que mi computadora sobreviviera a cualquier problema. Por 

hacerme sentir mejor cuando más cansada estaba, por sus palabras de tranquilidad cuando más 

las necesitaba, por todo el cariño y dedicación que me ha brindado. Y a su familia por ofrecerme 

ayuda tecnológica (internet) de forma ilimitada. 

A mi familia, por su colaboración incondicional. A mi abuela, por acompañarnos cuando 

más la necesitamos en casa, a mi tío Pedro por su apoyo y consejos ante cualquier inconveniente, 

a Mary por consentirme con sus tortas, a Alex por contar con el cual en lo que fuera, a Frank, a 

mi tío Edgar por estar disponible cuando más lo necesitábamos. 

A mis amigos de Charallave, por hacer más gratos los viajes de regreso a casa. Al COPA, 

el cual espero que prevalezca siempre. Al ventubus por ser fiel y siempre dispuesta a llevarme y 

traerme. 


v



A mis amigos del tercer piso, Geofísica, por los gratos momentos que compartimos 

(fiestas, pijamadas, campo, viajes) y por los más agobiantes (parciales, exposiciones, defensas, 

entrega de trabajos,…) que supimos superar, por su ayuda, chistes, cuentos, moralejas, consejos 

y su cariño: Cindy, Debora, Alejandra, Carolina, Cristina, Luis, Josex, Jostin, Adelo, Victor, 

Pedro, Walter, Andrés. 

A los del Centro de Procesamiento Geofísico en PDVSA Guaraguao, por su ayuda 

incondicional y las respuestas que me dieron cuando más las necesitaba, por hacerme sentir 

como en casa, a Johndelcy Guzmán, Anny Lara, la Sra. Xiomara Barrios, Ernesto Oropeza, 

Francisco Quintero, el Sr. Manuel Mendoza, Daniel Salas, el Sr. Horacio Galan, el Sr. Miguel Gil 

y especialmente a Luis Figuera por todo su apoyo, consejos, tiempo y dedicación permitiendo la 

culminación de este proyecto. A los asesores de Paradigm, Hector Pastini y al Sr. Jaime Yañes, 

por su conveniente orientación metodológica y la enseñanza que me brindaron. A los de soporte 

de Paradigm, Zahary, Yadelin y Eduardo, por toda su colaboración en los momentos más 

difíciles para entenderse con el software. A todos aquellos del edificio y Puerto la Cruz que me 

brindaron ayuda gustosamente cuando se las pedí. 

Y a todos aquellos que confiaron en mi, gracias. 


vi



Resumen. iv 

Agradecimientos. v 

Índice. vii 

Índice de Figuras. x 

I. Introducción. 1 

II. Fundamentos Teóricos. 3 

II.1 Propagación de Ondas. 3 

II.2 Prospección por Reflexión. 5 

II.3 Prospección por Refracción. 6 

II.3.1 Principio de Huygens. 6 

II.3.2 Ley de Snell. 7 

II.3.3 Principio de Fermat. 8 

II.4 Teoría del Modelo Convolucional. 9 

II.5 Cálculo de las Correcciones de Estáticas de Refracción. 10 

II.5.1 Selección de las Primeras Llegadas. 12 

II.5.2 Definición de Refractores. 13 

II.6 Recuperación de Amplitudes. 16 

II.6.1 Corrección por Divergencia Esférica. 16 

II.6.2 Corrección por Ganancia Exponencial. 18 

II.7 Deconvolución. 18 

II.7.1 Deconvolución Tipo Spike. 19 

II.7.2 Deconvolución Predictiva. 19 

II.8 Filtros de Frecuencias. 20 

II.8.1 Filtro Paso – Banda. 21 

II.8.2 Filtro de Corte Bajo y Alto. 21 

II.8.3 Filtro Tipo Knoch. 22 

II.8.4 Filtro Anti – Alliasing. 22 

II.9 Análisis de Velocidades. 22 

II.10 Migración. 25 

II.10.1 Principios de la Migración. 26 

II.10.2 Algoritmo Kirchhoff de Migración en Tiempo. 30 

II.10.2.1 Migración por Suma de Difracciones. 33 


vii



II.10.2.2 Migración por Superposición de Semicírculos. 34 

II.10.3 Parámetros del Algoritmo de Kirchhoff. 34 

III. Características Geomorfológicas y Geológicas del Levantamiento Sísmico Arecuna 

Oeste 04G – 3D. 

III.1 Localización del Levantamiento. 36 

III.2 Marco Geológico del Área. 37 

III.2.1 Formaciones Geológicas presentes en la zona. 39 

III.2.1.1 Formación Freites. 39 

III.2.1.2 Formación Las Piedras. 40 

III.2.1.3 Formación Mesa. 40 

III.3 Geomorfología del Área de Estudio. 41 

III.3.1 Mesa Conservada. 41 

III.3.2 Mesa Disectada. 41 

III.4 Características Hidrológicas. 42 

III.5 Importancia Económica del Levantamiento. 42 

IV. Descripción de los Datos Disponibles del Levantamiento Sísmico Arecuna Oeste - 

. 

IV.1 Parámetros de Adquisición del Proyecto. 44 

IV.2 Geometría del Proyecto. 46 

IV.3 Modelo de la Capa Meteorizada según Refracciones de Campo. 47 

V. Metodología y Análisis de los Resultados. 52 

V.1 Acondicionamiento de los Datos. 52 

V.1.1 Cálculo de las Estáticas de Refracción. 52 

V.1.2 Incremento Relación Señal – Ruido. 57 

V.1.3 Secuencia de Procesamiento. 65 

V.1.3.1 Análisis de Velocidades y Generación de Correcciones de Estáticas 

Residuales. 

V.1.3.1.1 Primer Análisis de Velocidades. 75 

V.1.3.1.2 Primer Cálculo de Correcciones de Estáticas Residuales. 77 

V.1.3.1.3 Segundo Análisis de Velocidades. 81 

V.1.3.1.4 Segundo Cálculo de Correcciones de Estáticas Residuales. 82 

V.1.3.1.5 Tercer Análisis de Velocidades. 85 

V.2 Migración Pre – Apilamiento en Tiempo con el Algoritmo de Kirchhoff. 91 

V.2.1 Cálculo de los Parámetros para la Migración. 91 


viii 

36 

44 

73



V.2.1.1 Análisis del Rango de Offset. 91 

V.2.1.2 Definición del Datum Flotante. 94 

V.2.1.3 Generación del Cubo de Velocidades. 96 

V.2.2 Primera Iteración. 98 

V.2.3. Determinación de los Residuos de Velocidad. 102 

V.2.4. Segunda Iteración. 104 

V.2.5. Tercera Iteración. 105 

V.2.6. Migración Post – Apilamiento en Tiempo en el Algoritmo de Kirchhoff. 107 

V.3 Generación de los Datos con Nuevas Correcciones de Estáticas, según la selección 

de las Primeras llegadas en todos los Disparos del Levantamiento. 

V.3.1 Cálculo de las Estáticas de Refracción. 112 

V.3.2 Análisis de Velocidades y Generación de Correcciones de Estáticas Residuales. 115 

V.3.2.1 Primer Análisis de Velocidades. 115 

V.3.2.2 Primer Cálculo de Correcciones de Estáticas Residuales. 117 

V.3.2.3 Segundo Análisis de Velocidades. 120 

V.3.2.4 Segundo Cálculo de Correcciones de Estáticas Residuales. 121 

V.3.2.5 Tercer Análisis de Velocidades. 124 

VI. Conclusiones y Recomendaciones. 130 

VII Referencias Bibliográficas. 133 

Apéndice A. 136 

Apéndice B. 141 

Apéndice C. 142 

Apéndice D. 144 


ix 

112



Fig II.1. Principio de Huygens. Tomado de http://es.wikipedia.org 6 

Fig II.2. Repartición de Energía en una interfase. 7 

Fig II.3. Modelo de las superficies de referencia de los datos sísmicos. Tomado de Estándares 

de presentación del Dato Sísmico en Venezuela. Ministerio de Energía y Minas. 

Fig II.4. Simulación recorrido onda refractada. 12 

Fig II.5. Representación gráfica (domocrónicas) para la estimación del número de reflectores. 13 

Fig II.6. Representación gráfica del recorrido del frente de ondas en el fenómeno de refracción. 13 

Fig II.7. Propagación del Frente de Onda en el Subsuelo (Divergencia Esférica). 17 

Fig II.8. Diseño de Filtro Pasa – Banda, en la izquierda. Diseño de Filtro Pasa – Banda con 

pendientes verticales, causando el fenómeno de Gibbs, en la derecha. 

Fig II.9. Geometría Convencional de los Levantamientos Sísmicos. 23 

Fig II.10. Registro Sintético en el que se observa (a) la reflexión de una interfase horizontal, y 

(b) la reflexión de la misma interfase pero corregida por nmo. 

Fig II.11. Principios geométricos de la Migración. 27 

Fig II.12. (a) Difractores de una estructura en profundidad. (b) Difracciones de una sección en 

tiempo. 

Fig II.13. (a) Trayectoria de los rayos desde el punto difractor hasta la superficie. (b) Sección en 

tiempo de una difracción. 

Fig II.14. (a) Impulso en una sección en tiempo. (b) Respuesta Impulso en el dominio migrado. 

(c) Curva de Difracción en una sección en tiempo. (d) Respuesta de migración a una curva de 

difracción en tiempo. 

Fig III.1. Mapa de Ubicación del Levantamiento Arecuna Oeste 04G – 3D. 36 

Fig III.2. Perfil Norte – Sur de la Cuenca Oriental de Venezuela. Tomado de WEC, 1997. 38 

Fig III.3. Perfil del Sistema Petrolero de la Faja del Orinoco. Tomado de Audemar et Al, 1985. 38 

Fig III.4 Mapa Estructural de las áreas de interés del Cuadrángulo Arecuna. Considerando que 

el levantamiento Arecuna Oeste 04G – 3D se ubica al Noreste del mapa, enmarcando las zonas 

1 y 2. 

Fig IV.1. Mapa de la dirección de las líneas de tiro y las receptoras. Los puntos rojos 

representan la ubicación de los disparos y las líneas negras en dirección Norte – Sur 

corresponden con los receptores. 

Fig IV.2. Mapa de Cobertura del Proyecto 46 

Fig IV.3. Gráfico tiempo vs offset de las refracciones superficiales. Proporcionando velocidad y 


x 

+ 

11 

21 

24 

31 

31 

33 

43 

45



espesor de la capa somera y la velocidad del primer refractor. 48 

Fig IV.4. Mapa de Contornos de las velocidades de la capa superficial (en la izquierda) y del 

primer refractor (en la derecha). 

Fig IV.5. Mapa de Contornos de las velocidades del segundo refractor. 49 

Fig IV.6. Mapa de Contornos de Espesores de la capa superficial (en la izquierda) y del primer 

refractor (en la derecha). 

Fig IV.7. Mapa Contornos de elevación para fuentes (en la izquierda) y receptores (en la 

derecha). 

Fig IV.8. Superficie de la base de la capa superficial. 51 

Fig IV.9. Superficie de la base del Primer Refractor. 51 

Fig IV.10. Superficie Topográfica. 51 

Fig IV.11. Superposición de las superficies de las figuras IV.8, IV.9 y IV.10. 51 

Fig. V.1. Tiempos de las primeras llegadas (tiempo vs offset). En la figura se observa la línea de 

tendencia por mínimos cuadrados, con una tolerancia máxima de factor 4. 

Fig. V.2. Selección primeras llegadas para los disparos 1163 y 1164. 53 

Fig. V.3. Curva tiempo de retardo para el disparo 1163 (en la derecha) y 1164 (en la izquierda). 54 

Fig. V.4. Curva tiempo de retardo para el disparo 1163 y 1164. 54 

Fig V.5. Correcciones de Estáticas totales, en la izquierda en el dominio de disparos y en la 

derecha de receptores. 

Fig V.6. Correcciones de Estáticas de refracción en el domino de cdp. 56 

Fig V.7. Efecto del Ground Roll en los registros del levantamiento. 57 

Fig V.8. Efecto de la aplicación del Low Frequency Array Filtering. 59 

Fig V.9. Comparación entre la aplicación de las funciones Low Frequency Array Filtering 

(disparo central) vs Time-Frequency Noise Supresión (disparo lado derecho). 

Fig V.10. Secuencia en la que fueron aplicadas las funciones para atenuar el ruido. 62 

Fig V.11. (a) Comparación disparo 9 con geometría (en la izquierda) vs disparo con atenuación 

del Ground Roll (en la derecha). (b) Comparación disparo 22 con geometría (en la izquierda) vs 

disparo con atenuación del Ground Roll (en la derecha). 

Fig V.12. Disparo con geometría, en la izquierda, y su correspondiente espectro de amplitud, en 

la derecha. 

Fig V.13. Disparo con atenuación del Ground Roll, en la izquierda, y su correspondiente 

espectro de amplitud, en la derecha. 

Fig V.14. Flujo de Trabajo para aplicar la deconvolución. 65 

Fig V.15. Comparación disparo 3250 sin deconvolución vs disparo 3250 con Deconvolución 


xi 

49 

50 

50 

53 

56 

61 

63 

63 

64



consistente con superficie aplicada. 66 

Fig V.16. Flujo de trabajo para balancear el espectro de amplitud de los datos. 66 

Fig V.17. Comparación shot_gather vs disparo con balanceo de amplitud consistente con 

superficie. 

Fig V.18. Espectro de amplitud de shot_gather (esquina superior izquierda), disparo con 

balanceo de amplitud consistente con superficie 

Fig V.19. Comparación de los shots_gather (esquina superior izquierda) y los disparos con 

balanceo de amplitud consistente considerando diferentes componentes con común fuente, 

receptor, offset o cdp. 

Fig V.20. Espectro de amplitud de la figura V.19; se distingue el mismo contenido de frecuencia 

sin (a) o con (b, c y d) la aplicación del balanceo de amplitud consistente con superficie. 

Fig V.21. Comparación disparo con balanceo de amplitud consistente con superficie vs disparo 

con balanceo espectral. 

Fig V.22. Comparación disparo con balanceo de espectro vs disparo con remoción de ruido con 

la función Time-Frequency Noise Supresión. 

Fig V.23. Comparación disparos con remoción de ruido con la función Time-Frequency Noise 

Supresión vs disparo con balanceo de amplitud sin consistencia con superficie. 

Fig V.24. Ecuaciones para determinar las estáticas según la elevación. 73 

Fig. V.25. Flujo de trabajo para la generación de un cubo aplicado corregido con estáticas de 

refracción y por elevación. 

Fig V.26. Sección apilada con correcciones de estáticas por refracción (lado izquierdo) y con 

correcciones de estáticas por elevación (lado derecho). 

Fig V.27. Flujo de trabajo para realizar el primer análisis de velocidades y generar el volumen 

corregido por nmo con éste análisis. 

Fig V.28. Selección de las primeras velocidades (vel_1A), en la izquierda. En la derecha, 

primera sección apilada. 

Fig V.29. Flujo de trabajo para generar un volumen piloto. 77 

Fig V.30. Flujo de trabajo para determinar las primeras correcciones de estáticas residuales. 78 

Fig V.31. Comparación sección apilada corregida con el primer análisis de velocidades sin 

estáticas residuales (stack_vel1A), en la izquierda, y con correcciones residuales 1, en la 

derecha. 

Fig V.32. Mapa de residuos A, segunda corrección. Lado izquierdo para receptores, lado 

derecho para fuentes. 

Fig V.33. Flujo de trabajo para realizar el segundo análisis de velocidades y generar el segundo 


xii 

68 

68 

69 

70 

70 

71 

72 

74 

74 

76 

76 

79 

80



cubo piloto. 81 

Fig. V.34. Flujo de trabajo para determinar las segundas correcciones de estáticas residuales. 82 

Fig V.35. Flujo de trabajo para aplicar las segundas correciones de estáticas residuales y generar 

un cubo apilado. 

Fig V.36. Comparación sección apilada, en la izquierda, con estáticas residuales 1 y primer 

análisis de velocidades y en la derecha con correcciones residuales 2 y segundo análisis de 

velocidades. 

Fig V.37. Mapa de residuos, segunda corrección. Lado izquierdo para estaciones, lado derecho 

para fuentes. 

Fig V.38. Flujo de trabajo para determinar el tercer análisis de velocidades y generar el volumen 

apilado corregido por nmo con éste análisis de velocidades. 

Fig V.39. Comparación sección apilada con segundo análisis de velocidades y segundas 

estáticas residuales (en la izquierda) y sección apilada con tercer análisis de velocidades y con 

segundas estáticas residuales aplicadas, stack_vel3A (en la derecha). 

Fig V.40. Flujo de trabajo para determinar los gathers con correcciones de estáticas según una 

traza piloto y para generar el volumen apilado aplicando ésta corrección. 

Fig V.41. Comparación sección apilada con tercer análisis de velocidades y segundas estáticas 

residuales, stack_vel3A (en la izquierda) y sección apilada con tercer análisis de velocidades, 

segundas estáticas residuales aplicadas y corrección de estáticas traza a traza, stack_cdptrim A 

(en la derecha). 

Fig V.42. Flujo de trabajo para aplicar un balanceo de amplitudes a los datos y removerles la 

corrección por nmo; generando los gathers de entrada para la migración pre – apilamiento. 

Fig V.43. Cubo de velocidades de apilamiento, primer, segundo y tercer análisis, 

respectivamente. 

Fig V.44. Selección del rango de offset con mayor cobertura. 92 

Fig V.45. Mapa de cobertura del rango de offset seleccionado. 93 

Fig V.46. Selección Pesos de las trazas. 93 

Fig V.47. Parámetros para el cálculo del Datum Flotante. 94 

Fig V.48. Superficie del Datum Flotante. 95 

Fig V.49. Parámetros para la generación de los time slice de velocidades rms (recuadro de la 

izquierda) y para los del Cubo de Velocidades rms (recuadro de la derecha). 

Fig V.50. Cubo de velocidades rms para la migración pre – apilamiento. 97 

Fig V.51. Parámetros de la migración pre – apilamiento en tiempo con el algoritmo de 

Kirchhoff. 


xiii 

83 

83 

84 

85 

86 

87 

88 

89 

89 

96 

100



Fig V.52. Comparación entre la sección apilada stack_cdptrimA, en la izquierda, y la sección 

apilada primera iteración de la migración pre – apilamiento en tiempo, en la derecha. 

Fig V.53. Selección de residuos de velocidad para la primera iteración. 102 

Fig V.54. Parámetros para la generación del cubo de residuos de velocidades. En la izquierda se 

muestra el menú para crear los slice de los residuos de velocidades y en la derecha el menú para 

la generación del cubo. 

Fig V.55. Cubo de velocidades rms, determinado con los segundos residuos de velocidades. 104 

Fig V.56. Comparación sección apilada stack_cdptrimA y sección apilada de la tercera iteración 

de la migración, inline 121, respectivamente. 

Fig V.57. Comparación sección apilada stack_cdptrimA y sección apilada de la tercera iteración 

de la migración, inline 81, respectivamente. 

Fig V.58. Flujo de trabajo para generar el cubo apilado para la migración post – apilamiento. 107 

Fig V.59. Menú para el cálculo de la migración post – apilamiento en tiempo con el algoritmo 

de Kirchhoff. 

Fig V.60. Comparación stack_cdptrimA, stack_cdptrimA con post – proceso 

(stack_post_migración) y migración post – apilamiento en tiempo, de izquierda a derecha 

respectivamente. 

Fig V.61. Comparación sección apilada de migración pre – apilamiento en tiempo y migración 

post – apilamiento en tiempo. 

Fig V.62. Comparación sección apilada, sección migrada post apilamiento y sección migrada 

pre apilamiento, en el dominio de tiempo, de izquierda a derecha respectivamente. 

Fig. V.63. Selección primeras llegadas para los disparos 4535 y 4536. 113 

Fig V.64. Correcciones de Estáticas totales, en la izquierda en el dominio de disparos y en la 

derecha de receptores. 

Fig V.65. Correcciones de Estáticas de refracción B en el domino de cdp. 114 

Fig V.66. Flujo de trabajo para realizar el primer análisis de velocidades y el primero cubo 

piloto. 

Fig V.67. Sección apilada con las correcciones de estáticas manuales (stack_vel1B), lado 

izquierdo. Sección apilada con las correcciones de estáticas automáticas (stack_vel1A), lado 

derecho. 

Fig V.68. Flujo de trabajo para determinar las primeras correcciones de estáticas residuales. 117 

Fig V.69. Comparación sección apilada corregida con vel_1A y aplicadas las correcciones 

residuales res1A, lado izquierdo; con la sección apilada corregida con vel_1B y aplicadas las 

correcciones residuales res1B, lado derecho. 


xiv 

101 

103 

105 

106 

108 

108 

110 

111 

114 

116 

116 

118



Fig V.70. Mapa de residuos, primera corrección. Lado izquierdo para estaciones, lado derecho 

para fuentes. 


cubo piloto. 

Fig V.72. Flujo de trabajo para determinar las segundas correcciones de estáticas residuales. 121 

Fig V.73. Comparación entre sección apilada corregida con segundo análisis de velocidades y 

aplicadas segundas correcciones residuales, pero corregidas con estáticas de refracción 

diferentes. 

Fig V.74. Mapa de residuos, segunda corrección. Lado izquierdo para receptores, lado derecho 

para fuentes. 


cubo piloto. 

Fig V.76. Comparación entre las secciones apiladas corregidas con el tercer análisis de 

velocidades pero con diferentes correcciones de estáticas de refracción aplicadas. 

Fig V.77. Flujo de trabajo para generar los datos corregidos con estáticas según una traza piloto 

y para generar un volumen apilado aplicando éstas correcciones. 

Fig V.78. Comparación sección apilada con tercer análisis de velocidades y segundas estáticas 

residuales, stack_vel3B (en la izquierda) y sección apilada con tercer análisis de velocidades, 

segundas estáticas residuales aplicadas y corrección de estáticas traza a traza, stack_cdptrim B 

(en la derecha) 

Fig V.79. Comparación entre la sección apilada corregida con segundas estáticas residuales y 

tercer análisis de velocidad. En la izquierda con estáticas de refracción según primeros quiebres 

seleccionados manualmente en todos los disparos, y en la derecha con estáticas de refracción 

determinadas con primeros quiebres seleccionados cada diez (10) disparos. 


xv 

119 

120 

122 

123 

124 

125 

126 

127 

128



Las áreas de interés exploratorio requieren datos geofísicos y geológicos para la 

interpretación del subsuelo, debido a que generan perfiles, volúmenes o secciones que permiten 

caracterizar y delimitar cualquier yacimiento. En la mayoría de los estudios los datos que 

facilitan el análisis de cualquier zona prospectiva son los levantamientos sísmicos, éstos generan 

imágenes en tiempo o profundidad que representan las estructuras del subsuelo, además, de 

proporcionar las características geométricas y las propiedades físicas de los estratos que son de 

gran utilidad en los estudios petrofísicos. Una sección o un cubo sísmico de buena calidad son 

fundamentales en la toma de decisiones para la ejecución de un pozo, siendo éste el que permite 

corroborar la información dada por la sísmica. Ante ello, es indispensable la realización de un 

procesamiento sísmico elaborado y detallado que genere secciones confiables para futuras 

interpretaciones. 

En este proyecto se pretende aplicar la secuencia convencional de procesamiento sísmico, 

hasta la migración pre – apilamiento en tiempo con el algoritmo de Kirchhoff, al Levantamiento 

sísmico Arecuna Oeste 04G – 3D, perteneciente a la estatal PDVSA (Petróleos de Venezuela 

S.A.). El cubo consta de un área de 70 Km 2 y la zona se encuentra enmarcada en el Municipio 

Francisco de Miranda al Sur Oeste del Edo. Anzoátegui, aproximadamente a 45 Km al Sur de la 

población el Tigre. 

El procesamiento sísmico será realizado con el software Focus 5.3 y la migración pre – 

apilamiento con el software Geodepth de Paradigm. La metodología a emplear se fundamenta en 


1



el acondicionamiento de los datos con el objetivo de mejorar la relación señal – ruido, debido a 

que los registros muestran gran influencia del ground roll. Por otro lado, por ser un 

levantamiento tan próximo al área de la Faja del Orinoco, los problemas de estáticas son difíciles 

de solventar. Seguidamente, se pretende aplicar las funciones convencionales para la generación 

de datos que proporcionen secciones con buena resolución de las reflexiones (ganancia por 

divergencia esférica, deconvolución, filtros de frecuencias, balanceo de amplitudes, ecualización 

de espectro, etc.). Además, se realizarán los respectivos análisis de velocidades, así como las 

correcciones estáticas residuales con el propósito de incrementar la coherencia y continuidad en 

los eventos de interés. 

Finalmente se generará el cubo de velocidades con el que se ejecutará la migración pre – 

apilamiento en tiempo con el algoritmo de Kirchhoff. Y éste cubo de velocidades será refinado 

con la selección de los residuos de velocidades dados por las diferentes iteraciones de la 

migración. El propósito final es obtener un cubo sísmico migrado en tiempo de buena resolución 

en el que las futuras interpretaciones estructurales sean factibles de realizar, proporcionando 

información confiable para posteriores perforaciones. 

El proyecto de pasantía larga fue realizado en las instalaciones del Centro de Procesamiento 

Geofísico de la Gerencia de Exploración y Producción, Departamento de Geofísica de PDVSA, 

división de Oriente, Puerto la Cruz, Estado Anzoátegui, bajo la tutela industrial del Dr. Horacio 

Galan y el Ing. Luis Figuera y la tutela académica de la Prof. Milagrosa Aldana. 


2



Los métodos de prospección símica se fundamentan en la propagación de un frente de 

ondas en el subsuelo, ésta propagación dará una respuesta interpretable gracias a las diferentes 

propiedades físicas de los estratos; tales como, las velocidades de propagación de las ondas en 

las capas rocosas, el contenido de fluidos, el volumen de espacio vacío (porosidad), la 

disposición estructural de las capas, etc. 

Estas propiedades permiten a los analistas generar modelos certeros, en los que se 

describen las condiciones geológicas en las que ocurrió la génesis de algún área de interés 

exploratorio, el arreglo estructural en el que se encuentra depositado, la descripción litológica, 

las características físicas de los fluidos que saturan dichas formaciones, y definir los limites del 

yacimiento analizado. 

Algunos de los principios básicos de la teoría de reflexión sísmica, la cual permite 

generar los datos empleados para las caracterizaciones y delineamientos de zonas en el subsuelo 

con interés económico se presentan a continuación. 

* ( !' # # * 

La teoría del movimiento ondulatorio, considerada como la propagación de energía de 

un lugar a otro sin transferencia de materia, mediante ondas mecánicas (propagándose en un 

medio material) o electromagnéticas (no requieren propagarse en algún medio material), es el 

principio teórico de los estudios de prospección sísmica. 

Onda Elástica: deformación que viaja a través de un medio elástico. 

Onda Sísmica: deformación que viaja a través del subsuelo. 


3



Frente de Onda: conjunto de puntos en el espacio que son alcanzados simultáneamente 

por una onda. Si se trazan líneas perpendiculares a los frentes de onda, se observa que dichas 

líneas indican la dirección en la que viajan las ondas. Estas líneas son denominadas rayos y son 

útiles para describir las trayectorias de la energía sísmica. 

La Propagación de Ondas se fundamenta en la solución del problema dinámico, es decir, 

la ecuación de onda: 

2 

2 

y( x, 

t) 

1 ∂ y( 

x, 

t) 

∂x 


4 

∂ 

2 

= 

V 

2 

∂t 

2 

……1 

Una onda al encontrar un cambio en las propiedades elásticas del medio en que viaja, 

reparte su energía según los principios de reflexión y refracción, éste fenómeno de partición de 

energía es empleado para generar los modelos geofísicos a partir de las distintas interfases que 

reporten la propagación de las ondas sísmicas en el subsuelo. 

La partición de la energía en el contacto entre dos medios elásticos, solo es resulta 

considerando varias condiciones de borde, tales como la continuidad de los desplazamientos 

normales y tangenciales, y la continuidad de los esfuerzos normales y tangenciales en la 

interfase [16], las ecuaciones que resuelven este fenómeno son las de Zoeppritz, las cuales 

dependen de parámetros como velocidades sísmicas, densidades de las rocas adyacentes a la 

interfase y las variaciones del módulo de Poisson; considerando que dichos factores están 

asociados a la litología y los fluidos que contienen las formaciones.



Las ecuaciones de Zoeppritz se aplican a una reflexión entre dos semi-espacios y no 

incluyen los efectos de interferencia de la ondícula causada por la estratificación; el sistema de 

ecuaciones se fundamenta en suposiciones como: frentes de ondas planos, interfase plana y 

uniforme, semi-medios homogéneos, isotrópicos y perfectamente elásticos. 

sinθ 

1 

− cosθ 

sin 2θ 

1 

1 

cos 2θ 

1 

cosφ 

sinφ1 

α1 

cos 2φ1 

β1 

β1 

− sin 2φ1 

α 

1 

− sinθ 

− cosθ 

2 

2 

ρ2β 

2 α1 

sin 2θ 

2 

2 

ρ1β1 

α 2 

ρ2α 

2 − cos 2θ 

2 

ρ α 


5 

1 

1 

1 

velocidadonda 

P. 

ρ = densidad. 

2 

cosφ 

− sinφ2 

ρ2β 

2α1 

− cos 2φ2 

ρ1β1 

ρ2β 

2 − sin 2φ2 

ρ α 

A = coeficiente 

de reflexión de la onda P 

B = coeficiente 

de reflexión de la onda Sv 

C = coeficiente 

de transmisión 

de la onda P 

D = coeficiente 

de transmisión 

de la onda Sv. 

α = 

β = velocidadonda 

S. 

1 

1 

2 

A 

B 

C 

D 

= 

− sinθ 

1 

− cosθ 

1 

− sin 2θ 

1 

− cos 2φ 

1 

……2 

Pero, por la complejidad del sistema de ecuaciones de Zoeppritz, se han realizado 

diferentes aproximaciones, como las de Aki y Richards (1980), Shuey (1985), Pan y Gadner 

(1987), Hilterman (1987). 

* ( !!' ( , -' * 

Técnica generada según las reflexiones que ocurren en las interfases de los estratos, 

debido a los cambios en las propiedades físicas de los mismos, permitiendo determinar la 

ubicación y los espesores de las estructuras geológicas que contengan yacimientos de interés 

exploratorio. El principio básico es la medición del tiempo que demora la onda en viajar desde el



punto donde fue generada (coordenada de la fuente – shot – ) hasta el lugar donde se encuentra 

el receptor (coordenada del geofono). [16] 

Mediante la reflexión sísmica se obtienen excelentes resultados en zonas donde se 

presenta sal, hielo, rocas ígneas y metamórficas; mientras que en zonas con sedimentos arcillos, 

los resultados disminuyen en resolución. La adquisición, el procesamiento y la interpretación de 

los datos empleando esta técnica son complejos y costosos. Sin embargo, esta técnica es 

considerada como la más eficiente para la exploración geofísica del subsuelo. [6] 

* ( !!' ( , !!' * 

Técnica fundamentada en las refracciones críticas de las ondas símicas ocurridas en la 

interfase entre dos medios elásticos diferentes. Estas son las primeras ondas en ser registradas 

por los receptores después de ser generado el pulso de energía por la fuente; debido a que solo 

se consideran aquellas ondas que siguen la trayectoria de tiempo mínimo. 

siguientes: 

Las leyes que rigen la propagación y trayectoria de éstas ondas sísmicas son las 

* 

“Todo punto alcanzado por un frente de onda actúa como fuente de nuevas ondas” 

(Huygens). Éste postulado proporciona un método geométrico para determinar, a partir de una 

frente de onda de forma conocida, la forma que adoptará éste frente de onda en un instante 

posterior. 

Fig II.1. Principio de Huygens. Tomado de http://es.wikipedia.org 


6



* 

“Independientemente de cuál sea el radio de la circunferencia considerada para un rayo 

de luz dado, la razón de los desplazamientos en uno y otro medio es la misma y depende sólo 

de la naturaleza de los medios elegidos y del orden en que se los hayan dispuesto” (W. Snell). 

n sinθ = n sinθ 

……3 


7 

1 

1 

2 

Donde θ son los ángulos de incidencia de los rayo y n es el índice de refracción de los 

medio y puede ser expresado como, 

2 

c 

n = ……4 

V 

Donde c es la velocidad de propagación de la luz en el vacío y V la velocidad de 

propagación de ondas en el medio; en el caso de análisis geofísicos, ésta velocidad corresponde 

con la propagación de ondas de cuerpos en los estratos del subsuelo. Entonces la ecuación 3, en 

función de la ecuación 4, queda expresada según: 

c 

V 

1 

c 

sinθ θ 

V sinθ = V sinθ 

……6 

1 = sin 2 ……5 2 1 1 2 

V2 

n 1 

n 2 

Fig II.2. Repartición de Energía en una interfase.



La ecuación 6 establece que si un rayo pasa de un medio de menor velocidad a otro de 

mayor velocidad el ángulo de incidencia aumentará. Análogamente, si en un rayo pasa de un 

medio de mayor velocidad a uno de menor velocidad, el ángulo de incidencia disminuirá. 

* 

“La trayectoria de un rayo de luz que va de un punto a otro es la correspondiente al 

tiempo mínimo” (P. Fermat). Éste principio establece que una onda para dirigirse de un punto a 

otro recorre la trayectoria en la que se emplee el tiempo mínimo de recorrido. 

Finalmente, debido a las propiedades de propagación de las ondas refractadas, la 

prospección de refracción es empleada para el estudio de zonas someras (entre 0 y 100 m de 

profundidad). Ésta técnica es eficiente para delimitar la interfase entre dos medios con fuerte 

contraste de velocidad de propagación de ondas sísmicas, por ello es utilizada para definir el 

modelo de la capa meteoriza respecto al primer refractor de interés. [6] 

Empleando los tiempos de la selección de los primeros quiebres dados por las 

refracciones de las ondas en la interfase entre la base de la capa meteorizada y el primer estrato 

rocoso de alta velocidad, se grafican las domocrónicas, las cuales permiten determinar los 

velocidades y espesores de las capas someras. Con tales parámetros se genera el modelo de la 

capa meteorizada y se determinan las correcciones de estáticas que deben aplicarse a los datos 

sísmicos para remover el efecto que causa dicho estrato, debido a que como es una capa de baja 

velocidad provoca distorsión en la continuidad de la hipérbola que representa las reflexiones de 

lo eventos de interés. 


8



* ) # # % "! * 

La traza sísmica real es producto de la convolución entre la serie de reflectividad y una 

función que representa las distorsiones e interferencias de los efectos causados por la 

propagación del impulso artificial en el subsuelo, además se le debe agregar un ruido aleatorio. 

modelo. 

Por tanto, la expresión matemática que representa el sismograma real viene dada por: 

S = w r + n 


9 

( ) ( ) * [16] 

t t ( t ) ( t) 

……6 

Donde w(t) = resultado neto de la distorsión convolucional de todos los efectos del 

r(t) = serie de reflectividad del subsuelo 

n(t) = ruido aditivo. 

Este modelo fue construido según las siguientes asunciones [18]: 

1a. “La Tierra esta constituida por capas horizontales de velocidad constante”. 

1b. “La fuente genera una onda plana compresional que incide perpendicularmente en las capas. 

Bajo tales circunstancias, las ondas generadas no son truncadas”. 

2. “La onda generada por la fuente es de tipo estacionaria, no cambia de forma durante la 

propagación”. 

3. “La componente del ruido es cero, n(t) (0)”. 

4. “La forma de la onda generada por la fuente es conocida”. 

5. “La serie de reflectividad es un proceso aleatorio. Esto implica que el sismograma presenta 

características de ondículas sísmicas en las que su autocorrelación y espectro de amplitud son 

similares”. 

6. “La ondícula sísmica es de fase mínima. Por consiguiente, tiene una fase mínima inversa”. 

Muchas de las asunciones anteriores no son cumplidas. Sin embargo, el principal 

objetivo del procesamiento sísmico es eliminar o atenuar el efecto causado por w(t) y determinar 

correctamente la serie de reflectividad.



La función w(t) esta compuesta por diferentes elementos, 

( wF 

( t) 

* wR( 

t ) * wP( 

t) 

) . A( 

t) 

w = [16] ……7 


10 

( t) 

Donde wF(t) = Ondícula de Campo, considerando fuente, fantasmas (focos 

considerados como fuentes y receptores secundarios), respuesta de los instrumentos, respuestas 

de los arreglos de fuentes y receptores. 

wR(t) = reverberaciones causadas por el entrampamiento de la señal en 

capas de baja velocidad (múltiples de largo periodo). 

wP(t) = efecto de propagación (absorción, múltiples de período corto). 

A(t) = factor de amplitud (atenuación de la amplitud sísmica debido a la 

propagación del frente de onda en el subsuelo). 

* .!" # !! # $.$! # , !!' * 

La determinación de este parámetro es fundamental para la remoción del efecto que 

causa la capa meteorizada en la sección apilada. Esta capa somera de baja velocidad distorsiona 

y degrada los eventos, dificultando el análisis de los intérpretes. 

Las correcciones de estáticas son aplicadas para remover los tiempos de retardo que 

causan la capa meteorizada, debido a las variaciones de espesor, velocidad y elevaciones de 

fuentes y receptores. Por tanto, permiten referenciar los datos sísmicos a una misma superficie 

de adquisición (datum plano o datum flotante); ver figura II.3.



Fig II.3. Modelo de las superficies de referencia de los datos sísmicos. Tomado de Estándares de presentación del 

Dato Sísmico en Venezuela. Ministerio de Energía y Minas. 

El cálculo de las correcciones estáticas por refracción viene dado por la determinación de 

los espesores y velocidades de la capa meteorizada, definiendo el modelo geológico de la capa 

somera considerando los primeros arribos de las ondas compresivas refractadas en la base de la 

capa meteorizada y las distancias fuente – receptor. Los pulsos refractados, simplificando el 

modelo con una sola capa, viajan desde la fuente a la interfase con la velocidad de la capa 

meteorizada, luego como ondas refractadas, viajan en la interfase con la velocidad de la capa 

refractante y finalmente, viajan desde la interfase con la velocidad de la capa meteorizada al 

receptor (Cerveny y Ravindra, 1971). La simulación del recorrido realizado por una onda 

refractada es mostrada en la figura II.4. 


11



Fig II.4. Simulación recorrido onda refractada. 

Según la ley de Snell para ondas críticamente refractadas, se tiene 

sin 

( θ ) sin( 

90° 

) 

V1 

= ……8 sin( 

θ ) = 1 

V w 


12 

V 

1 

V 

w 

……9 

Considerando la geometría de la figura II.4, los primeros arribos, la definición del 

refractor y el cálculo de la velocidad del mismo; todos estos parámetros obtenidos mediante la 

graficación de las domocrónicas, se puede determinar el espesor y velocidad de la capa 

meteorizada. 

* 

Fuente 

Etapa en la cual se asigna el tiempo de llegada del primer frente de ondas, según la 

misma fase del refractor. La interpretación de una sección apilada con correcciones de estáticas 

proveniente de una excelente selección de primeros quiebres será más sencilla si se compara con 

una sección con selección de primeros quiebres deficientes. Por ello, es recomendable dedicarle 

tiempo y esfuerzo a esta fase del procesamiento. 

θ θ 

V w 

V 1 

Receptor 

Z



* 

La determinación del número de refractores se obtiene mediante el análisis gráfico de la 

densidad, comportamiento y tendencia de los valores estimados con los tiempos de las primeras 

llegadas vs los offsets (ver figura II.5). La dependencia al offset permite excluir las reflexiones de 

las capas someras, las llegadas directas y refracciones de la primera capa. En cuanto a la 

delimitación de los offsets lejanos, se debe considerar un umbral donde la pendiente de la línea 

de tendencia varíe; con ello se define trabajar con un modelo de una sola capa. 

Primeras Llegadas (ms) 

Tiempo de las Primeras Llegadas vs Offset 

Offset (m) 


13 

dos capas 

una capa 

Fig II.5. Representación gráfica (domocrónicas) para la estimación del número de reflectores. 

Fig II.6. Representación gráfica del recorrido del frente de ondas en el fenómeno de refracción. 1 

1 Tomada de Alan Lembang, Practical Aspect of Seismic Data Adquisition. 2003



De la gráfica de la figura II.5 se puede demostrar que la velocidad del medio en que viaja 

la onda viene dada por el inverso de las pendientes de las rectas en la gráfica; y el modelo de la 

capa meteorizada puede ser completado considerando la imagen de la figura II.6, de ella se 

puede demostrar que el espesor de la capa meteorizada puede ser determinado según las 

siguientes ecuaciones: 

t = t 

t = 

V 

según la ecuación 9, sin( 

) 

t = 

t = 

+ t 

+ t 


14 

AB 

w 

BC 

2Z 

x − 2Z 

tanα 

……10 

c 

+ 

cosα 

V 

V w 

1 

c 

CD 

α c = , tenemos 

V 

2Z 

t = 

+ 

V sinα 

cosα 

1 

x 2Z 

+ 

V V cosα 

1 

x 

V 

1 

1 

1 

c 

2Z 

cosα 

c 

+ = 

V sinα 

c 

c 

c 

x 

V 

1 

x 

V 

c 

1 

2Z 

sinα 

c 

− 

V cosα 

1 

− sinα 

c 

sinα 

1 

1 

w 

c 

2Z 

cosα 

c 

+ 

V 

……11 

entonces, los valores de las primeras llegadas corresponden con la siguiente ecuación: 

t 

2Z 

cosα 

c 

i = ……12 

Vw 

Por tanto, el espesor de la capa meteorizada es calculado según la siguiente expresión: 

Z 

i 

Vwt 

= 

2 

i 

1− 

1 

V 

V 

w 

1 

2 

[2]……13 

donde Vw es la velocidad de la capa meteorizada y V1 la velocidad del primer refractor. 

Es de considerar que la aplicación de las correcciones de estáticas debe ir paralelamente 

con la aplicación de estáticas por elevación, puesto que variaciones abruptas en los espesores de



la capa meteorizada pueden provocar distorsiones apreciables en los tiempos de llegada de las 

reflexiones [12]. Se sabe que si se tienen elevaciones abruptas lo más probable es tener espesores 

también abruptos de la capa meteorizada; de igual forma en el caso análogo. Por tanto, en zonas 

con elevaciones sin grandes variaciones la aplicación de correcciones de estáticas por elevación 

sería suficiente, mientras que en levantamientos con elevaciones moderadas la calidad de los 

datos sísmicos aumentará si se aplican correcciones de estáticas por refracción. Sin embargo, 

debe considerarse que en zonas con intensa actividad geológica (efectos compresionales, 

causando sobrecorriemento de capas, o afloramientos de material muy consolidado, producto 

de erosión) ninguno de los dos métodos de correcciones de estáticas, refracción o elevación, 

proporcionarán una sección apilada fácil de interpretar [12]. 

Por otro lado, la presencia de la capa meteorizada puede no afectar en gran magnitud la 

sección apilada, lo que causará el mayor inconveniente será la variación en el espesor y la 

velocidad de dicha capa. Es por ello que son calculadas las estáticas de corto y largo periodo. 

Las estáticas de largo periodo son causadas por cambios laterales graduales, las cuales 

introducen estructuras aparentes en la sección sísmica en tiempo. Mientras que las estáticas de 

corto periodo son ocasionadas por cambios laterales rápidos, los cuales introducen variaciones 

de tiempo dentro de las trazas individuales de un punto común en profundidad [12]. Estos dos 

tipos de estáticas corrigen cada traza sísmica, posicionándola en tiempo sobre el modelo 

geológico diseñado (superficie teórica idealizada). 


15



* !"( !' # ( $"# 

Proceso necesario para compensar las pérdidas por atenuación de amplitud debido a la 

propagación del frente de ondas en el subsuelo. Es decir, cada vez que la onda se propaga va 

incidiendo en alguna interfase de contraste acústico y parte del ella es reflejada y parte es 

transmitida; por tanto, la onda que continúa descendiendo pierde energía, por el contenido de 

fluidos en el espacio poroso y la distancia recorrida por dicho frente de onda. [18] 

Las funciones de ganancia permiten compensar las pérdidas de amplitud que sufren las 

reflexiones por dispersión esférica, por la reflexión misma y por pérdidas térmicas. Dichas 

funciones pueden ser la corrección por divergencia esférica o por ganancia exponencial. 

* 

La energía producida por un disparo se va atenuando a medida que se propaga en el 

subsuelo debido a la ampliación del frente de onda; de tal manera que a mayor profundidad el 

radio de dicha esfera será mayor, mientras que la amplitud disminuye como lo muestra la figura 

II.7. La divergencia esférica que se produce desde el lugar de disparo, provoca que la energía 

disminuya inversamente al cuadrado de la distancia que el frente de onda ha viajado. [16] 


16



Fig. II.7. Propagación del Frente de Onda en el Subsuelo (Divergencia Esférica). 

La divergencia esférica varía según sea el medio en que se propaguen los impulsos: 

1.a) Si la propagación del frente de ondas ocurre en un medio homogéneo, el 

decaimiento de la amplitud viene dado por: 

L Z 

( ) = o 

Z 


17 

1 

L 1 

( t ) V. 

t 

= ……14 

1.b) Si el medio es heterogéneo, el decaimiento de las amplitudes viene dado por: 

L 

V 

1 

( t) 

= ……15, 

2* 

V( 

t) 

. t 

V 

2* 

( t) 

= 

Divergencia 

Esférica 

N 

2 

Vi 

t 

t 

i= 

1 0 

i 

……16 

donde V1 es la velocidad rms del primer reflector al que se le aplicará la corrección y t es 

V1.t 

el tiempo doble de viaje de la onda hasta la profundidad dada Z. 

Z 

Disparo 

V2.t 

= 1/Z 2



N 1 

2 

Finalmente, el factor de corrección por amplitud será: = G ( 1− 

R ) 

* ! 


18 

∏ − 

A( t) 

( t) 

i 

i= 

1 

Función matemática aplicada a los valores de amplitud para compensar la atenuación de 

las mismas. Es de considerar que la amplitud relativa será alterada con la aplicación de dicha 

función, con el objetivo de mejorar la relación señal – ruido. 

La ecuación que rige la recuperación de las amplitudes mediante la ganancia exponencial 

en el dominio de decibles (Anthony, 1983), es: 

Ganancia ( db) 

= a. 

t + 20. 

b. 

log( t) 

+ c ……17 

Donde t es el tiempo, a el factor de atenuación, b el factor de divergencia esférica y c el 

factor de ganancia constante. 

La corrección se fundamenta en aplicar la curva inversa dada por la ecuación anterior; de 

tal forma que varían ligeramente las amplitudes, así como las amplitudes verdaderas relativas, 

disminuyendo la atenuación de las mismas. 

* ! % "!' 

Su principal objetivo es comprimir las ondículas sísmicas, logrando mejorar la resolución 

de las reflexiones primarias de manera que la sección apilada sea una representación confiable 

de la función de reflectividad del subsuelo, reduciendo las reverberaciones de la señal sísmica y 

el ancho de la ondícula [18]. Por otro lado, nivela el espectro de amplitud, tanto de las 

componentes coherentes de la señal, como del ruido coherente o no coherente. Este último hecho 

puede ser corregido con la aplicación posterior de algún filtro. 

V 

V 

1 

2* 

t 

[16]



Matemáticamente, actúa como un filtro inverso que intenta obtener de la señal de salida 

dada por el subsuelo, la señal de entrada (el impulso generado en superficie) y poder deducir la 

serie de reflectividad que representa las interfases en el subsuelo. 

* " # 

Se fundamenta en convertir cada reflexión primaria en un impulso con una longitud 

predictiva de una muestra de duración en tiempo y asume que la señal reflejada de cada 

interfase, es de fase mínima [15]. La técnica consiste en emplear la teoría de filtro predictivo y 

generar el filtro más apropiado de tal manera que el resultado de la convolución, entre las trazas 

originales y dicho filtro, sean trazas similares a la tipo spike [12]. La habilidad del analista estará 

en definir adecuadamente los parámetros del filtro predictivo para que la señal de salida sea la 

mejor representación la serie de reflectividad. Sin embargo, en muchas ocasiones deben 

elaborarse gran cantidad de pruebas y hallar los parámetros de deconvolución que mejor se 

adapten a los datos analizados. 

* 

Este tipo de deconvolución permite controlar la resolución de la ondícula deseada, 

debido a que puede atenuar eventos periódicos empleando una distancia predictiva. Usa 

información de la primera parte de la traza para predecir y remover el efecto de la ondícula de la 

parte final de la traza sísmica [12]. De tal manera, que la distancia predictiva discrimina parte de 

la ondícula que se desea preservar y la que se desea remover. 


19



La atenuación de efectos no deseados mediante la aplicación de ésta deconvolución se 

fundamenta en la teoría del filtro de Wiener – Levinson. La longitud del operador del filtro es 

determinado con la solución del sistema matricial de Robinson y Treitel, 1980. 

r 

r 

r 

r 

0 

1 

2 

. 

. 

. 

n−1 

n−2 

n−3 


20 

r 

r 

r 

r 

1 

0 

1 

. 

. 

. 

r 

r 

r 

r 

2 

1 

0 

. 

. 

. 

. 

. 

. 

. 

. 

. 

. 

. 

. 

. 

. 

. 

. 

. 

. 

. 

. 

. 

. 

. 

. 

r 

r 

n−1 

n−2 

. 

. 

. 

. 

r 

0 

. 

f 

f 

f 

f 

. 

. 

. 

0 

1 

2 

n−1 

= 

r 

r 

α + 1 

r 

r 

α 

α + 2 

. 

. 

. 

α + n−1 

……18 

donde r(t) es la autocorrelación de la señal de entrada, f(t) es el filtro de Wiener – Levinson 

y r(t+α) la croscorrelación (ver apéndice A) entre la salida deseada desplaza en una fase α y la 

salida deseada. Como el operador predictivo es un filtro que actuará en la señal de entrada en 

un tiempo posterior a t, determinando la traza en un tiempo t + α . Entonces, el inicio de la 

aplicación de la deconvolución predictiva debe ser diseñado con el objetivo de minimizar los 

efectos caudados por los ruidos periódicos no deseados. 

* + $ # + !" ! 

Funciones, de fase cero y simulando una respuesta impulsiva, aplicadas a los datos 

sísmicos con el propósito de remover componentes de frecuencias desligadas de la señal de las 

reflexiones primarias. En el dominio f,t la discriminación entre las frecuencias asociadas al ruido 

de las asociadas a las reflexiones sísmicas es más sencilla, por ello el análisis de estas funciones 

se aplica en dicho dominio, y los datos son transformados mediante la Transformada Rápida de 

Fourier.



El principio básico es eliminar las frecuencias correspondientes al ruido no deseado de 

los datos a procesar, en la mayoría de los algoritmos se diseñan funciones boxcar con las cuales 

serán convolucionados los datos para identificar de forma concreta la señal que se desea 

remover. Considerando que el rango de frecuencia no deseado es comúnmente conocido, tal 

como el del ground roll. Sin embargo, no deben eliminarse bandas de frecuencias sin un análisis 

detallado y minucioso que permita corroborar que las componentes removidas están 

verdaderamente ligadas al ruido. 

* $ % 

Es diseñado para eliminar bandas de frecuencias, definiendo un máximo y un mínimo, 

pero sin incluir valores con pendientes verticales en el dominio f,t (como en la figura II.8b); 

debido a que se originaría la deformación de la señal, según el fenómeno de Gibbs (ver apéndice 

A.1.7). Por ello, es necesario diseñar con cuatro (4) componentes la banda de frecuencias a 

remover (como en la figura II.8a). 

F1 

F2 

Fig II.8. (a) Diseño de Filtro Pasa – Banda. 

(b) Diseño de Filtro Pasa – Banda con pendientes verticales, causando el fenómeno de Gibbs. 

* % & ' 

Tipo de filtro fundamentado en eliminar umbrales de frecuencias definidos por el 

analista, ya sean de valores bajos o altos, respectivamente. 


21 

a b 

F3 

F4 

F1 

F2



* ( ) 

Se emplea para eliminar algún valor determinado de frecuencia, comúnmente la 

componente asociada (60 Hz) a la presencia de cables de alta tensión o cualquier estructura 

eléctrica que sea registrada durante la grabación del levantamiento. 

* ' *' 

Aplicado a los datos sísmicos para preservar el periodo de Nyquist y cumplir con la 

condición de muestreo. La cual propone que para registrar adecuadamente las frecuencias de 

una señal de banda limitada se debe cumplir con la siguiente ecuación: 

f N 

* . # !# # 

1 

≤ ……19 

2∆t 

Procedimiento en el procesamiento convencional basado en corregir en tiempo la curva 

hiperbólica que representa las reflexiones en las interfases. Esta corrección, conocida como 

normal moveout (NMO), requiere de la selección de velocidades que mejor horizontalicen los 

eventos representados hiperbólicamente por el incremento entre la distancia fuente – receptor 

(offset). De tal forma, que los datos serán transformados a cero offset mediante la disminución del 

tiempo de viaje de un rayo reflejado en algún punto del reflector en profundidad (Castellanos, 

2004). Por tanto esta corrección se fundamenta en simular incidencia normal, generando datos 

que puedan ser apilados. 


22 

Donde fN = frecuencia de Nyquist 

∆t = periodo de muestreo



Se sabe que por la geometría de los levantamientos sísmicos, el tiempo de viaje de la 

reflexión de una interfase horizontal que viaja desde la fuente hasta el receptor viene dado por la 

ecuación 20 y puede ser demostrada según la figura II.9. 

4Z 

x 

t + 

V V 

2 2 

2 

= 2 2 ……20 [16] 

Fig II.9. Geometría Convencional de los Levantamientos Sísmicos. 

La representación de la ecuación anterior en un gráfico offset vs tiempo será una 

hipérbola, como muestra la figura II.10a. Y los análisis de velocidades serán realizados para 

obtener trazas corregidas por normal moveout, como muestra la figura II.10b. Esta corrección 

viene dada según las siguientes ecuaciones: 

t t z t − = ∆ ……21 

t 

( 0) 

= ……22 

( 


23 

) 

0 

Z 

2V 

Con las ecuaciones 20, 21 y 22, tenemos: 

Donde Z = profundidad hasta la interfase 

x = distancia fuente – receptor 

V = velocidad de propagación del medio 

Donde ∆t = corrección por normal moveout. 

t(Z) = tiempo doble de viaje de la onda 

desde la fuente, a la interfase, hasta el receptor 

t(0) = tiempo doble de viaje de la onda si 

incidiera perpendicularmente en los estratos



t 

2 

( Z ) 

= t 

2 

( 0) 

x 

+ 

V 

2 

x 

( t( 

Z ) + t0 

)( t( 

Z ) − t0 

) = 2 


24 

t 

2 

2 

V 

x 

2 

∆t = 

……23 2 

V ( t( 

Z ) + t0 

) 

[16] 

Fig II.10. Registro Sintético en el que se observa (a) la reflexión de una interfase horizontal, y 

(b) la reflexión de la misma interfase pero corregida por nmo. 

Las ecuaciones anteriores asumen que la trayectoria de las ondas sísmicas ocurre en 

capas horizontales, la onda se propaga con trayectoria recta y la velocidad de propagación es 

constante desde el plano de referencia, hasta la interfase reflectora. [15] 

a 

b 

La aplicación de la corrección nmo requiere conocer la velocidad del medio en el que se 

propagan las ondas. Mediante los análisis de velocidades este parámetro será estimado, pero 

bajo la figura de velocidad de apilamiento. La cual es la velocidad que corrige cada evento, 

llevándolo a la posición donde la reflexión obtenida estaría si los datos estuviesen adquiridos a 

cero (0) offset, [15]. La selección de estas velocidades será guiada según la máxima semblanza, la 

cual permite estimar los valores que mejor horizontalizan los eventos en los registros; 

considerando a que a mayor profundidad el espectro de semblanza resulta cada vez menos 

x 

t(0) 

∆t 

t(Z)



confiable, pues se hace más difícil seguir los máximos de coherencia para determinar la 

velocidad correcta. Este espectro de semblanza consiste en generar varias hipérbolas de la 

ecuación 20 y corregir por nmo los datos, apilarlos y desplegar los máximos picos de amplitud o 

también llamados máximos de coherencia. De manera que el máximo de coherencia en la 

semblanza, indicará el valor de la velocidad de apilamiento que mejor horizontalizará los 

eventos de los datos analizados. 

En el caso que los cdp muestreen distintos puntos de una misma interfase horizontal, la 

velocidad de apilamiento será igual a la velocidad cuadrática media (rms), ver apéndice B. Pero 

si la capa buza, la traza será posicionada en el punto medio fuente – receptor, y el cdp no estará 

posicionado en el punto de reflexión, [15]. Por ello, en zonas con elevada complejidad estructural 

los resultados no presentan buena resolución, requiriendo el proceso de migración. 

!' 

Los datos sísmicos facilitan el análisis de cualquier zona de interés, permitiendo la 

interpretación de la misma con el uso de perfiles, volúmenes o secciones sísmicas. Estas 

imágenes sísmicas se generan mediante el despliegue de las reflexiones en su posición correcta, 

lo cual se fundamenta en dos (2) procesos: apilamiento y migración. 

La migración se obtiene con las diversas soluciones de la ecuación de onda que describe 

la propagación de ondas elásticas a través de las rocas. Esta emplea un modelo de velocidad 

para redistribuir la energía sísmica reflejada, desde la posición supuesta en el punto medio a su 

verdadera posición; no obstante, se remueve el efecto de propagación de ondas llevando los 

reflectores con buzamiento a su posición verdadera en espacio y tiempo y se colapsa la energía 


25



de difracción a su punto de origen, cuya difracción ocurre por discontinuidades repentinas en 

los distintos tipos de roca. Por tanto, el objetivo de la migración sísmica es proporcionar una 

imagen con mayor resolución espacial de los eventos. [17] 


26 

+ 

Sea un reflector con buzamiento real en profundidad como el mostrado en la geometría 

de la figura II.11a, la sección donde los receptores y fuentes están ubicados en la misma posición 

(sección a cero offset) es mostrada en la figura II.11b. Ésta figura corresponde con la sección en 

tiempo del reflector real, pero muestra un buzamiento aparente menor. Esta imagen es 

construida si se grafica la primera y última reflexión en el reflector buzante. 

La primera llegada de incidencia normal grabada en A es mostrada por el punto C’. Con 

el transcurso del tiempo se siguen grabando reflexiones normales del reflector CD y la última 

llegada será graba en B y corresponde con el punto D’. El reflector C’D’ corresponde con la 

representación del reflector real CD en una sección en tiempo, en la imagen de la figura II.11b se 

observa que el buzamiento del reflector es más pequeño y está desplazado tanto vertical como 

lateralmente. 

El cambio de posición será realizado por el proceso de migración (figura II.11c), 

provocando un aumento en la inclinación de los reflectores, una disminución en la longitud del 

reflector y se mueve el evento en dirección buzamiento arriba. [18] 

Las variaciones causadas por la migración en la representación de los reflectores son 

fundamentales para representar, lo mejor posible, a los eventos geológicos reales. Estas



imágenes serán obtenidas asumiendo secciones a cero offset y los rayos deben seguir patrones 

de incidencia normal a las interfases. Además, el buen resultado del proceso de migración 

sísmica estará condicionado a la buena aproximación entre la sección en tiempo y la sección a 

cero offset, la relación señal – ruido y las velocidades empleadas para corregir el normal moveout. 

a 

b 

c 

Fig II.11. Principios geométricos de la Migración. 


27



Los fundamentos matemáticos del proceso de migración sísmica son basados en la teoría 

de ecuación de onda acústica en dos (2) dimensiones: 

∂ 

∂x 

1 

− 

V 


28 

2 

2 

∂ 

+ 

∂z 

2 

2 

2 

∂ 

∂t 

2 

2 

P 

( x, 

z, 

t ) = 0 ……24 

Donde x es la variable espacial horizontal, z la variable espacial vertical (profundidad), t 

la variable del tiempo y P(x,z,t) el campo de onda que se propaga en un medio con densidad 

constante y velocidad de onda compresional V(x,z). 

El principal objetivo del algoritmo de migración es determinar la reflectividad P(x,z,0) a 

partir del campo de onda P(x,0,t) grabado en superficie. [18] 

Aplicando la transformada de Fourier al campo de onda en el plano x, asumiendo que no 

existe variación lateral de velocidad, se tiene la siguiente ecuación: 

P 

= P , 

( k , z, 

ω ) ( x, 

z t ) 

Y la transformada inversa, será: 

x 

( x, 

z, 

t ) ( k , z ω ) 

e 

( ik xx 

−iωt 

) 

dxdt ……25 

( −ik 

x+ 

iωt 

) 

x 

P = P , e dkxdω 

x 

……26. 

Si el operador diferencial de la ecuación 24 es aplicado en la ecuación 26, se tiene: 

∂ 

∂z 

2 

2 

P 

2 

ω 

V 

2 

( z, 

ω ) + − k x P( 

k z, 

ω ) = 0 ……27 

kx 2 

x 

La solución de la ecuación diferencial 24, para el campo de onda ascendente, será: 

2 

ω 

2 

= k . z 

2 

V 

P( k z ) P( 

k ) − i − 

x , ω 

exp 

x 0 , ω 

x 

……28



Por otro lado, se define el número de onda vertical kz según la ecuación 29. 

2 

k z = − k 2 

V 

ω 


29 

2 

x 

……29 

Sustituyendo la ecuación 29 en la ecuación 28, se tiene: 

P 

( k z , ω ) = P( 

k 0 , ω ) 

x 

x 

e 

( − ik z ) 

z 

……30 

El principio físico representado por la ecuación 30 es que el campo de onda que se 

propaga en profundidad (P(x,z,ω))puede ser determinado a partir del campo de onda grabado en 

superficie (P(x,0,t)). Además, éste campo de onda emergente grabado en superficie “puede ser 

descompuesto en un plano de ondas monocromático” [18], donde Kz es dependiente de la 

profundidad, entonces: 

2 

k z( 

z) 

= − k 2 

V 

ω 

2 

x 

……31 

Sustituyendo la ecuación 31 en la ecuación 30, se tiene: 

P 

= 

( k z , ω ) ( k 0 , ω ) 

x 

P 

x 

e 

− i 

z 

0 

k 

z 

( z ) 

dz 

……32 

Con esta expresión se representa que el campo de onda puede continuar en profundidad 

en un medio donde la velocidad varía verticalmente con la profundidad. Al mismo tiempo, el 

proceso de migración incluye la condición de imagen, ésta condición establece que se deben 

sumar de todas las frecuencias cada vez que se realiza la extrapolación del campo de onda en 

profundidad. [18]



El tipo de migración a elegir, tanto el algoritmo como el dominio de aplicación, depende 

de la complejidad del objetivo y las estructuras de sobrecarga de interés. Dicha selección se 

fundamentará en los cambios laterales de velocidad, debido a que en zonas donde los estratos 

no presentan leves buzamientos y los cambios laterales de velocidad son moderados, se aplica la 

migración en tiempo. Mientras que, en zonas de elevada complejidad estructural con 

variaciones laterales significativas de velocidad, se emplea la migración en profundidad. Las 

técnicas de migración más usadas son por suma de difracciones, por diferencias finitas en el 

dominio de Fourier, etc. 

' ( )) + " 

El proceso de migración según el algoritmo de Kirchhoff es generado mediante el 

método de difracción o el de sumatorias hiperbólicas. La geometría del método puede ser 

ilustrada con el modelo de reflectores explotantes dado por Loewenthal en 1976. Éste modelo 

asume que los reflectores emiten rayos y cada punto que constituye un reflector del subsuelo 

puede visualizarse como una fuente de energía que genera un tren de onda ascendente; los 

patrones de rayo son trazados a incrementos constantes de tiempo; y los frentes de onda son 

grabados en superficie. Entonces, un reflector sísmico puede ser observado como si estuviese 

compuesto de puntos difractores (figura II.12a) y la migración de esa sección sísmica será 

realizada colapsando las hipérbolas de difracción a su ápice (figura II.12b). 


30



Fig II.12. (a) Difractores de una estructura en profundidad. (b) Difracciones de una sección en tiempo. 

Finalmente, cada punto de la sección migrada es generado al sumar todas las amplitudes 

a lo largo de la curva de difracción, cuya curvatura estará gobernada por una función de 

velocidad, y asignando el resultado al ápice de la hipérbola [13]. 

Por otro lado, la curva de difracción se obtiene a partir del trazado de rayos desde el 

punto de difracción en el subsuelo hacia todos los puntos en superficie (según la figura II.13a) y 

la ecuación que describirá la curva de difracción (según la figura II.13b) viene dada por: 

Donde t 

0 

a b 

( x − x ) 


31 

2 

0 

2 

2 

t = t0 

+ ……33 

V 

2z 

= 0 y V la velocidad del medio hasta el punto difractor. 

V 

Z 

Z0 

a 

Fig II.13. (a) Trayectoria de los rayos desde el punto difractor hasta la superficie. 

(b) Sección en tiempo de una difracción. 

t 

T0 

b



En caso de ser una sección a cero offset para datos antes de apilar, la ecuación que define 

la curva de difracción viene dada por la siguiente ecuación, conocida como la ecuación de raíz 

cuadrara doble: 


32 

2 2 

( x − x ) 2 ( x − x ) 

r 

2 

RMS 

2 

s 

2 

RMS 

2 

t = t0 

+ + t0 

+ ……34 

ν 

ν 

Donde xr y xs son las posiciones de los receptores y las fuentes, respectivamente. Ésta 

ecuación define el tiempo de viaje de la trayectoria de un rayo que viaja desde la fuente, pasa 

por el punto difractor y llega hasta el receptor. Con el valor de velocidad rms y la ecuación de 

raíz cuadrada doble se determina la superficie de difracción. Es de considerar, que el algoritmo 

de migración antes de apilar presenta fuerte dependencia a la precisión de la solución dada por 

los análisis de velocidad. [4] 

Por otro lado, si el problema de la figura II.13 es visto de forma inversa, según la figura 

II.14. Es decir, se selecciona un punto fijo (x,t) en el espacio de entrada a la migración, entonces, 

la respuesta en el espacio de salida (sección migrada) de ese punto será una elipse, según la 

ecuación 33. Por ello, un punto o una fuente en una sección en tiempo constituyen un 

semicírculo en espacio migrado. Éste semicírculo es denominado respuesta impulso de la migración 

(ver figura II.14a e II.14b). En un esquema alternativo, una curva de difracción en la sección en 

tiempo, será un punto en el dominio migración (ver figura II.14c e II.14d); éstos modelos 

originan dos métodos de migración [18]:



a b c d 

x 

Fig II.14. (a) Impulso en una sección en tiempo. (b) Respuesta Impulso en el dominio migrado. (c) Curva de 

Difracción en una sección en tiempo. (d) Respuesta de migración a una curva de difracción en tiempo. 

+ ,fue el primer algoritmo de migración 

implantado en la era de las computadoras. Se basa en la suma de las amplitudes a lo largo de las 

trayectorias hiperbólicas que definen a las curvas de difracción. La ecuación que describe esas 

trayectorias, asumiendo un modelo de velocidades en capas horizontales y usando la velocidad 

rms en el ápice de la hipérbola, viene dada por la siguiente expresión: 

2 

4x 

Vrms 

2 2 

( ) = t0 

+ ……35 

2 


33 

t x 

El comportamiento del frente de onda a lo largo de la hipérbola de difracción está afectado 

por diferentes factores, los cuales se describen a continuación [18]: 

a) Factor de Oblicuidad: asociado a la dependencia de las amplitudes con el ángulo en que 

cada evento llega a cada receptor. Antes de la suma, las amplitudes deben ser 

reescaladas por el coseno del ángulo θ, según la figura II.13d. 

b) Factor de Divergencia Esférica: asociada al decaimiento de las amplitudes, debido a que la 

energía decae en un factor de 1/r 2, por tanto las amplitudes deberán ser reescaladas por 

un factor de 1/r antes de la suma, para compensar el efecto de propagación de ondas en 

tres (3) dimensiones. 

θ 

t0



c) Factor de la Forma de la Onda: la onda resultante de la sumatoria debe tener un espectro de 

fase y amplitud proporcional a la frecuencia a 90° de fase constante (para levantamientos 

3D). 

+ - ,implantada antes de la era 

de las computadoras, se fundamenta en la repartición de amplitudes a lo largo de elipses en el 

espacio migrado de salida, superponiendo las respuestas impulsos. La imagen es generada con 

la construcción de la envolvente de todas las elipses por interferencia constructiva. Y los bordes 

de las curvas son cancelados por interferencia destructiva. 

parámetros [18]: 

. ' ( )) 

En la migración con el algoritmo de Kirchhoff, deben considerarse los siguientes 

a. Apertura de Migración: asociado a la extensión espacial que tendrá la hipérbola de 

difracción, viene medida en términos de número de trazas del trayecto que abarcará la 

hipérbola. El ancho de éste parámetro, incrementa con la profundidad, lo que implica 

que eventos más profundos migrarán más que los someros; es decir, se tendrá mejor 

corrección para eventos profundos en comparación con los someros. 

b. Máximos de buzamientos de migración: corresponde al rango de buzamientos que el 

algoritmo considerará, la magnitud puede ser estimada según conocimientos geológicos. 

Dentro de las ventajas de este algoritmo es que maneja adecuadamente grandes 

buzamientos y obtiene, en el desarrollo de las iteraciones, líneas migradas y sus 

correspondientes secciones de cdp (ver apéndice C) que permiten el proceso iterativo en la 


34



construcción del modelo de velocidades; además el tiempo de respuesta (generación de la 

imagen migrada) es muy rápido en comparación con otros algoritmos. Por otro lado, las 

desventajas son que para el cálculo de los tiempos de viaje considera solo las primeras llegadas y 

no maneja adecuadamente variaciones de velocidades laterales. 


35



! ! !" !" #$ #$ !# 

!# 

El proyecto de prospección sísmica Arecuna Oeste 04G – 3D se encuentra 

aproximadamente a 45 Km al Sur de la población El Tigre, ubicándose político-territorialmente 

en el Municipio Francisco de Miranda al Sur Oeste del Estado Anzoátegui. 

Se encuentra delimitado geográficamente por las coordenadas de los vértices del 

recuadro fucsia mostrado en la figura III.1. 

12° 

8° 

72° 

60° 

Fig III.1. Mapa de Ubicación del Levantamiento Arecuna Oeste 04G – 3D. 


36



% % ! ! & & !'%! 

!'%! 

El levantamiento sísmico esta localizado en la Cuenca Oriental de Venezuela, la cual es 

una depresión geográfica y estructural situada en el Centro-Este del país, limitada al Norte por 

la Serranía del Interior; al sur con el Escudo de Guayana; al Este con el Golfo de Paria y el 

Océano Atlántico, y al Oeste con el arco de El Baúl. El inicio de la cuenca se sitúa en el Devono – 

Carbonífero, para dicho período se reconocen tres (3) periodos sedimentarios. 

Se desarrollo como estructura de geosinclinal a partir del Cretácico con la invasión de los 

mares sobre el basamento peneplanado. Durante su formación, el eje del geosinclinal comenzó a 

desplazarse desde el Norte hacia el Sur hasta llegar a la parte central de Anzoátegui y Monagas, 

en el Plioceno. Durante éste periodo en gran parte de la cuenca los ambientes que se 

establecieron fueron los fluvio-deltáicos y neríticos costeros. La posición del eje del geosinclinal 

controló la deposición y el carácter de los sedimentos, los cuales se vuelven generalmente más 

marinos de Sur a Norte y de Oeste a Este. La culminación del proceso sedimentario en la Cuenca 

Oriental esta representado por la Formación Mesa, de ambiente continental de sedimentos 

lutíticos y limolíticos. La figura III.2 muestra un perfil en dirección Norte – Sur de la cuenca, 

distinguiendo algunas Formaciones que la conforman. 

Específicamente, el proyecto se encuentra enmarcado dentro de la Faja del Orinoco, 

ubicada al Sur de la Cuenca Oriental, y es considerada como la acumulación de crudos pesados 

y extrapesados más grande del mundo (ver perfil del sistema petrolero en la figura III.3). Los 

yacimientos de petróleo de estas áreas se encuentran en arenas no consolidadas, de origen 

deltáico, con porosidades promedio de 30% [3]. 

Los sedimentos presentes en la zona de estudio son de origen terciario, cuaternario y 

cretácico, según el reporte de campo. La litología característica esta conformada por areniscas, 

arcillas, margas, lutitas, calizas, limolitas, gravas y cuarcitas. [5]. 


37



2 Audemar et Al., 1985. 

Fig III.2. Perfil Norte – Sur de la Cuenca Oriental de Venezuela. Tomado de WEC, 1997. 

Fig III.3. Perfil del Sistema Petrolero de la Faja del Orinoco. 2 


38



+ ! ' ! ( $ / * 


39 

/012 

Edad: Terciario (Mioceno Medio a Mioceno Tardío Basal) 

Descripción litológica: En la localidad tipo, Hedberg et al. describen lutitas 

físiles con areniscas en el tope y la base, que permiten la subdivisión de la 

unidad en tres intervalos. Jam y Santos (1987) observaron que la formación se vuelve más 

arenosa en cuanto más se adelgaza en el área Cerro Negro de la Faja Petrolífera del Orinoco. Las 

variaciones verticales en el contenido de arcilla, lutita, limo y siderita permiten una correlación 

extraordinaria, sobre la mayor parte del flanco sur de la cuenca, por variaciones en la curva de 

resistividad amplificada de los perfiles eléctricos. 

Contactos: suprayace concordantemente a la Formación Oficina en casi toda su 

extensión, excepto en el área de Anaco. Hacia el tope, la Formación es concordante con la 

Formación Las Piedras suprayacente. El contacto superior se identifica con el tope del mas alto 

nivel marino de Freites. Funkhouser et al., (1948) observaron que Las Piedras (Sacacual) 

sobrelapa discordantemente a Freites en el área de Anaco. 

Paleoambiente: En la mayor parte de la cuenca se presenta un ambiente marino somero 

en su proporción inferior, pasando a ambientes de aguas algo más profundas en la parte media. 

La parte superior corresponde a ambientes de aguas llanas. Hedberg y Sass, 1936 (fide De Sisto, 

1961, p. 101), postularon un ambiente estancado y ligeramente ácido por debajo del nivel de 

oleaje. Audemard et al. (1985) identifican que el carácter de la formación es nerítico. Hacia el 

margen sur de la cuenca, en el sector Cerro Negro de la faja petrolífera del Orinoco, la 

Formación disminuye de espesor al paso que se hace mucho más arenosa. En dicho sector, el 

ambiente de sedimentación es de tipo litoral, según la fauna de moluscos (Jam y Santos, 1987). 

Importancia económica: Las lutitas y arcillas de la formación constituyen el sello 

principal de los reservorios petrolíferos de la Formación Oficina en el área Mayor de Oficina, y 

contienen, en la parte inferior, arenas productoras de hidrocarburos en algunos campos del 

parte media y sur de esa área.




40 

/012 

Edad: Mioceno Tardío – Plioceno. 

Descripción litológica: En la localidad tipo, la Formación consiste en 

areniscas micáceas, friables, de grano fino, interlaminada con lutitas, 

arcilitas sideríticas, lutitas ligníticas y lignitos (Hedberg, 1950). También se encuentran algunas 

calizas arenosas. Constituida por sedimentos de origen marino de textura fina y mal 

consolidados. A pesar de tener una composición litológica consistente y uniforme, el 

predominio de uno u otro material varía tanto vertical como horizontalmente. 

Contactos: En la parte norte de la cuenca de Maturín, la formación es concordante sobre 

la Formación La Pica. En el flanco sur de la cuenca, yace concordantemente sobre la Formación 

Freites. En el extremo sur de la cuenca, cerca del río Orinoco, la formación yace posiblemente 

sobre el basamento ígneo metamórfico, al desaparecer la Formación Freites por acuñamiento 

(Jam y Santos, 1987). 

Paleoambientes: En el área del campo Pedernales, la formación fue depositada en un 

ambiente deltáico a marino somero (Barnola, 1960, p. 561). 

Importancia económica: Las arenas de la Formación son productoras de petróleo pesado 

en los campos de Orocual, Manresa y Pirital, en el norte del estado Monagas. 

+ 

/012 

Edad: Cuaternario (Pleistoceno) 

Descripción Litológica: conformada por arenas de grano grueso y gravas 

con mucho cemento ferruginoso, generalmente con alto grado de 

cementación y endurecimiento, dando lugar a conglomerados muy duros, arenas brillantes con 

fuerte estratificación cruzada, lentes discontinuos de arcilla finamente arenosa y lentes de 

limolita. 

Contactos: suprayace a la Formación Las Piedras en forma concordante y transicional.



Paleoambientes: Según González de Juana (1946), la Formación es producto de una 

sedimentación fluvio-deltáica y paludal, resultado de un extenso delta que avanzaba hacia el 

este en la misma forma que avanza hoy el delta del Río Orinoco. El mayor relieve de las 

cordilleras septentrionales desarrolló abanicos aluviales que aportaban a la sedimentación 

clásticos de grano más grueso, mientras que desde el sur el aporte principal era de arenas. 

Coplanarh (1974) considera que los sedimentos de la formación representan depósitos 

torrenciales y aluviales, contemporáneos con un levantamiento de la Serranía del Interior. 

, ) # 0 # $"# 

El levantamiento se ubica en la Provincia Fisiográfica de los Llanos, Subregión Llanos 

Orientales. En el sector domina el relieve de Altiplanicie de Mesa, con pendientes planas (1 a 

2%) y onduladas (3 a 5%). Esta Altiplanicie esta conformada por dos tipos de relieve [5]: 

+ : formada a partir de areniscas y arcillas, mostrando los 

afloramientos de la Formación Las Piedras. Constituye un relieve casi plano a moderadamente 

ondulado, ligeramente afectado por una erosión laminar generalizada que lo convierte en un 

medio morfodinámicamente inactivo. Existen áreas más erosionadas hacia el Este de la zona 

cubierta por el levantamiento sísmico. 

+ : con coluviones retomados de la Formación Mesa. El paisaje esta 

formado por la alternancia de los topes de mesa con sus farallones. Se localiza en el centro Norte 

del levantamiento, su litología está conformada por areniscas, material del cuaternario que ha 

sido afectado por la erosión regresiva ocurrida al retirarse las aguas, originando un relieve 

moderadamente disectado y un medio morfodinámicamente activo con problemas graves de 

erosión laminar, lo que ocasiona la presencia de cárcavas. 


41



!$ )$! # ' ! 

El área de estudio se localiza en la cuenca del río Orinoco, en la subcuenca del Río Pao. El 

principal cuerpo de agua presente en la zona es el Río Hamaca. El cual presenta escorrentía en 

dirección Norte – Sur, régimen dendrítico, poca velocidad y pendiente. Este río presenta una 

marcada influencia por el régimen estacional de las precipitaciones de la zona; y es factor 

determinante en la presencia de la unidad de vegetación Morichal. [5] 

( $ ! ! ' ! # % $ $ 

El estudio geológico realizado en el levantamiento adyacente por el Este (Cuadrángulo 

Arecuna) al proyecto Arecuna Oeste 04G – 3D, sugirió la exploración del sector analizado en este 

informe. El reporte indica la presencia de dos (2) zonas de interés exploratorio ubicadas en el 

área enmarcada por el levantamiento Arecuna Oeste. 

Según la figura III.4, el área denotada por el número uno (1) permitió redefinir los límites 

de ocho (8) yacimientos, correspondientes a las arenas H3L,4; I3; I4,5; I6; J3; N2; U1 y U3 (ver 

apéndice D). Dichos límites indicaron la posibilidad de perforar nuevos pozos hacia el sector 

Oeste para drenar eficientemente los yacimientos [8]. En cuanto al área denotada por el número 

dos (2), la situación no es tan clara puesto que solo se presenta un pozo (MFA–92) y los límites 

de los yacimientos se dificulta. Sin embrago, nuevos estudios deben ser dirigidos a la generación 

de pozos en dicha zona, debido a que se presentan ocho (8) yacimientos de interés reportados 

por el pozo MFA–92 , dados por las arenas L3,4; R1; R2; S1,2; S3,4; S5 y U3, pero sus límites no 

son conocidos con precisión. [8] 


42



2 

1 

Fig III.4 Mapa Estructural de las áreas de interés del Cuadrángulo Arecuna. Considerando que el levantamiento 

Arecuna Oeste 04G – 3D se ubica al Noreste del mapa, enmarcando las zonas 1 y 2. [8] 


43



. $ # #1" !' # 2 !$ 

El proyecto Arecuna Oeste 04G – 3D fue definido como un levantamiento sísmico 

tridimensional para la adquisición de los datos, realizada por la compañía BGP Internacional of 

Venezuela, S.A. 

En el levantamiento las líneas de tiro se desarrollaron con un acimut de 90° y en las 

líneas receptoras el acimut fue de 180°, el estacado se incremento de Oeste a Este y de Norte a 

Sur, respectivamente. El área del levantamiento sísmico recubre 70 Km 2 con un diseño 

geométrico cuadrado, ver figuras III.5 e III.6, según los siguientes parámetros: 


44 

( 

Líneas de Tiro Líneas Receptoras 

Número de líneas 33 36 

Separación entre líneas 250 m 250 m 

Separación entre estacas 50 m 50 m 

800 canales 100 canales con 8 líneas receptoras 

Los patrones de perforación para la adquisición corresponden con dos (2) modelos: uno 

con un pozo de 10 m de profundidad, cargado con 800 gr y un fulminate; y el otro un pozo de 

10 m de profundidad cargado con 1 Kg y un fulminate. Los explosivos empleados para la 

detonación fueron: 

)



- Booster de Pentolita Biodegradable GEOPRIME TM de 1 Kg y 800 gr. Producto 

autodestructible con bacterias consumidoras de pentolita en caso de falla de tiro y 

biodegradable. Velocidad de detonación 8.230 m/s, densidad 1,6 g/cc y alta resistencia al 

agua. [5] 

- Detonadores Eléctricos Sismográficos de 3 m y 17 m de longitud. Detonadores de 

sensibilidad normal con filamento AWG-24 (∅ 0,5mm), resistencia a la presión hidrostática de 

10 Kg/cm2. [5] 

Fig IV.1. Mapa de la dirección de las líneas de tiro y las receptoras. Los puntos rojos representan la ubicación de los 

disparos y las líneas negras en dirección Norte – Sur corresponden con los receptores. 


45



* $) # 2 !$ * 

Fig IV.2. Mapa de Cobertura del Proyecto. 

El levantamiento dispone de: 

Número de Inline 321 

Dirección de las Inline Este – Oeste 

Incremento de las Inline Norte – Sur 

Número de Crossline 351 

Dirección de las Crossline Norte – Sur 

Incremento de las Crossline Oeste – Este 

Intervalo de Muestreo 2 ms 

Longitud del Registro 3 s 

Offset 50 – 2800 m 

Tamaño del Bin 25 x 25 m 

Número de Fuentes 5509 

Número de Estaciones 5760 

Número de Trazas ~ 3.648.000 


46



* # # ( $ / # 3 , !! # ( * 

El Levantamiento Sísmico dispone de datos de refracción somera, estos datos permiten 

generar un modelo inicial de las capas superficiales. El tendido fue realizado con 500 gr de carga 

a 1 m de profundidad, generando 70 refracciones. Los parámetros del tendido son los 

siguientes [5]: 

Parámetros del tendido: 999 m 

Números de canales: 48 

0ffset: 0.5 m 

Profundidad: 1 m 

Carga: 500 gr. 

Muestreo: 0.25 milisegundos 

Ganancia de grabación: 36 dB 

Ganancia de ploteo: ACG o fix gain 

Características del instrumento: 

Cámara: Stratavisor Geometrics 

Geófono: SM_24 


47



Fig IV.3. Gráfico tiempo vs offset de las refracciones superficiales. Proporcionando velocidad y espesor de la capa 

somera y la velocidad del primer refractor. 

Con las domocrónicas de gráficos como el mostrado en la figura IV.3 se genera el mapa 

de contornos de las velocidades (ver figura IV.4 y IV.5) y espesores de la capa superficial y el de 

velocidades del primer y segundo refractor (ver figura IV.6). 


48

958000 

957000 

956000 

955000 

954000 

953000 

952000 

951000 

950000 



949000 

336000 337000 338000 339000 340000 341000 342000 343000 344000 345000 

660 956000 

640 

620 

600 

955000 

580 

540 

954000 

520 

480 

953000 

460 

360 

340 951000 


49 

560 

500 

440 

420 

400 

380 

320 

300 

958000 

957000 

952000 

950000 

949000 

336000 337000 338000 339000 340000 341000 342000 343000 344000 345000 

Fig IV.4. Mapa de Contornos de las velocidades de la capa superficial (en la izquierda) 

y del primer refractor (en la derecha). 

958000 

957000 

956000 

955000 

954000 

953000 

952000 

951000 

950000 

949000 

336000 337000 338000 339000 340000 341000 342000 343000 344000 345000 

Fig IV.5. Mapa de Contornos de las velocidades del segundo refractor 

Los valores de velocidades mostrados en la figura IV.4 permiten distinguir que el primer 

refractor presenta propiedades físicas muy similares a la capa superficial. Mientras que las 

velocidades observadas en la figura IV.5 son de mayor magnitud que las de los estratos 

suprayacentes, lo que sugiere que el modelo de la capa meteorizada esta constituido por una 

sola capa. 

1805 

1800 

1795 

1790 

1785 

1780 

1775 

1770 

1765 

1760 

1755 

1750 

1745 

1740 

1735 

1730 

1725 

1720 

1715 

1710 

1705 

1700 

940 

920 

900 

880 

860 

840 

820 

800 

780 

760 

740 

720 

700 

680 

660 

640 

620 

600

958000 

957000 

956000 

955000 

954000 

953000 

952000 

951000 

950000 



949000 

336000 337000 338000 339000 340000 341000 342000 343000 344000 345000 


50 

-4 

-5 

-6 

-7 

-8 

-9 

-10 

-11 

-12 

-13 

-14 

-15 

-16 

-17 

-18 

-19 

-20 

-21 

-22 

-23 

-24 

-25 

958000 

957000 

956000 

955000 

954000 

953000 

952000 

951000 

950000 

949000 

336000 337000 338000 339000 340000 341000 342000 343000 344000 345000 

Fig IV.6. Mapa de Contornos de Espesores de la capa superficial (en la izquierda) y de la capa meteorizada (en la 

derecha). 

Con los espesores de las capas someras (según la figura IV.6) y la superficie topográfica 

(ver figura IV.7) se puede generar el modelo del primer refractor (fig IV.8 y IV.9), permitiendo 

comparar los valores de las velocidades y las características geométricas del modelo generado 

por las refracciones según los primeros quiebres y el modelo dado por las refracciones someras; 

por tanto, se verificará el apropiado resultado de las correcciones de estáticas dadas por el 

modelo que sugerirá la selección de las primeras llegadas. 

958000 

957000 

956000 

955000 

954000 

953000 

952000 

951000 

950000 

949000 

336000 337000 338000 339000 340000 341000 342000 343000 344000 345000 

250 

245 

240 

235 

230 

225 

220 

215 

210 

205 

200 

195 

190 

185 

180 

175 

170 

165 

160 

155 

150 

145 

140 

135 

949000 

336000 337000 338000 339000 340000 341000 342000 343000 344000 345000 

Fig IV.7. Mapa de Contornos de elevación para fuentes (en la izquierda) y receptores (en la derecha). 

958000 

957000 

956000 

955000 

954000 

953000 

952000 

951000 

950000 

0 

-5 

-10 

-15 

-20 

-25 

-30 

-35 

-40 

-45 

-50 

-55 

-60 

-65 

-70 

250 

245 

240 

235 

230 

225 

220 

215 

210 

205 

200 

195 

190 

185 

180 

175 

170 

165 

160 

155 

150 

145 

140 

135 

130



Fig IV.8. Superficie de la base de la capa superficial. 

Fig IV.9. Superficie del tope del Primer Refractor. 

Fig IV.10. Superficie Topográfica. 

Fig IV.11. Superposición de las superficies de las figuras IV.8, IV.9 y IV.10. 


51 

220 

215 

210 

205 

200 

195 

190 

185 

180 

175 

170 

165 

160 

155 

150 

145 

140 

135 

205 

200 

195 

190 

185 

180 

175 

170 

165 

160 

155 

150 

145 

140 

135 

130 

230 

225 

220 

215 

210 

205 

200 

195 

190 

185 

180 

175 

170 

165 

160 

155 

150 

145



! # ! $ # # $ * 


52 

' 

.!" # $.$! # , !!' * 

Las correcciones de estáticas por refracción son calculadas considerando la selección de 

los primeros quiebres en los disparos del levantamiento, debido a que en la selección de las 

primeras llegadas viene inherente la información de las refracciones, el efecto de la capa 

meteoriza y la topografía. La selección de los primeros arribos empleando el software Focus 5.3 

puede ser automatizada eligiendo un rango de análisis (cada diez disparos, para este proyecto) y 

el programa con la metodología de redes neuronales estimará la solución para los disparos no 

seleccionados. Se debe considerar que se dispuso de la selección de las primeras llegadas 

realizada previamente a la iniciación de este proyecto (ver fig. V.1). 

La solución dada por el programa serán las estáticas totales para los disparos y las 

estaciones, por tanto es necesario separarlas en el dominio de las de periodo corto y largo, y 

estas últimas permitirán remover el efecto de la capa meteorizada. 

Con la aplicación de la función Generalized 3D Refraction Statics Generation del software 

Focus 5.3 se obtiene la solución del modelo de la capa meteorizada mediante una aproximación 

por mínimos cuadrados, considerando el cálculo de los tiempos de retardo, consistencia con 

superficie y topografía [11]. Dicho modelo fue generado considerando un solo refractor con 

velocidad próxima a 1800 m y con contribución de offset entre 500 y 2800 m. Este modelo 

s 

permite obtener las componentes de periodo largo asociadas a la geometría y velocidad del 

primer refractor y las estáticas por refracción de corto periodo, asociadas a la capa más 

superficial y los efectos causados por la topografía [1].



No obstante, con estas dos últimas componentes se evitará el efecto de la capa 

meteorizada mediante la determinación de las estáticas para referenciar la sísmica al datum plano 

o al datum flotante. 

En las figuras V.1, V.3 y V.4 se muestra el gráfico de las domocrónicas generadas por el 

software a partir de la selección de los primeros quiebres y en la figura V.2 se observa la interfaz 

para seleccionar las primeras llegadas. 

Fig. V.1. Tiempos de las primeras llegadas (tiempo vs offset). En la figura se observa 

la línea de tendencia por mínimos cuadrados, con una tolerancia máxima de factor 4. 

0.15 s 

Fig. V.2. Selección primeras llegadas para los disparos 1163 y 1164. 


53 

0.7 s



Fig. V.3. Curva tiempo de retardo para el disparo 1163 (en la derecha) y 1164 (en la izquierda). 

Fig. V.4. Curva tiempo de retardo para el disparo 1163 y 1164. 

En la figura V.4 se observan dos tendencias en las primeras llegadas enmarcadas por el 

cuadro rojo, pudiendo interpretarse como la errónea selección en la fase del primer quiebre. Sin 

embargo, la figura V.2 demuestra que la selección del primer arribo es correcta. En ésta imagen 

se observa la continuidad en la selección de la primera llegada para dos disparos diferentes, en 

uno el primer arribo se registra a los 0.15 s y en el otro a los 0.7 s. Además, los datos se 

aproximan adecuadamente a una la línea de tendencia (según la fig V.3). Por tanto, la diferencia 

en las líneas de tendencia de las primeras llegadas de los disparos analizados (fig V.4) es debida 


54



a los distintos tiempos de arribo, posiblemente asociados a diferencias litológicas en las capas 

superficiales, zonas con distinto coeficiente de compactación, etc. 

Comparando el modelo de la capa meteorizada dado por las refracciones de campo con 

el generado por las refracciones según los primeros quiebres, se distingue que este último 

modelo no proporcionó buena información para offsets cercanos, por tanto las primeras capas 

con velocidades entre 400 y 800 m no son reconocidas; mientras que en las refracciones de 

s 

campo sí fueron analizadas (según la figura IV.4 del capítulo IV, sección IV.3). 

Por otro lado, el rango de velocidades del segundo refractor (según modelo de campo, 

ver fig IV.5 del capítulo IV, sección IV.3) es muy similar al obtenido con el modelo de las 

refracciones según los primeros quiebres, permitiendo garantizar el apropiado cálculo de las 

correcciones de estática de refracción dadas por la función aplicada. 

Posterior a la selección de las primeras llegadas y la generación del modelo de la capa 

meteorizada, se emplea la función CDP Mean Statics Computation, la cual permite estimar las 

correcciones de estáticas para referenciar los datos a una misma superficie. Dicha función 

necesita los datos de entrada en el dominio de cdp (ver apéndice C) y los primeros arribos 

editados (sin puntos divergentes a la línea de tendencia). La solución será dos correcciones de 

estáticas, una de ellas de menor magnitud permitiendo referenciar los datos al datum flotante, y 

disminuyendo la distorsión en la hipérbola que representa la reflexión de los eventos, sin 

olvidar que esta corrección debe ser aplicada después de corregir por normal moveout para 

referenciar los datos al datum plano. Y la otra corrección de estática, de mayor componente en el 

dominio de fuente y receptor, que permite referenciar los datos al datum final. En las figuras V.5 

y V.6 se observan las estáticas obtenidas en el dominio de estaciones, receptores y cdp. 


55



Fig V.5. Correcciones de Estáticas totales, en la izquierda en el dominio de disparos, en un rango de -39 a 57 ms; 

y en la derecha de receptores, en un rango de -52 a 44 ms. 

Fig V.6. Correcciones de Estáticas de refracción en el domino de cdp, al removerlas permiten referenciar 

los datos al datum flotante y al aplicarlas al datum plano. Oscilan en el rango de -86 a 96 ms. 

Las correcciones de estáticas de refracción en el dominio de receptores, disparos y cdp 

(figuras V.5, V.6) presentan un comportamiento muy similar entre sí; además, se distingue que 

tienen una tendencia análoga al de la superficie topográfica (ver figura IV.7 del capítulo IV, 

sección IV.3) en el carácter de correcciones de estáticas. Es decir, para elevaciones topográficas 


56



altas presentes en la zona central del levantamiento, la corrección de estática es con signo 

negativo; mientras que para menores elevaciones de superficie, observadas en los laterales de la 

sección, la corrección de estática es con signo positivo. 

! $ !' 4 " # * 

Disponiendo de los shots con geometría se propone acondicionar los datos de entrada 

para la migración pre – apilamiento en tiempo. Dicha preparación estará fundamentada en la 

eliminación completa o parcial del ground roll y/o todos aquellos efectos asociados a anomalías 

en el momento de grabación (ruido ambiental, artificial, etc), con el objetivo de obtener datos 

que permitan definir con mejor resolución los eventos en el volumen sísmico. 

Empleando diferentes funciones, basadas en la eliminación de ciertas frecuencias, se 

logró obtener registros en los cuales el cono de ruido es atenuado. La figura V.7 permite 

caracterizar el efecto del ground roll en los registros del levantamiento Arecuna Oeste 04G - 3D, 

discriminando que la mayor frecuencia de esta señal es próxima a 18 Hz. 

Fig V.7. Efecto del Ground Roll en los registros del levantamiento. El recuadro amarillo indica 

la zona del análisis espectral mostrada en la ventana del lado derecho. 


57



Las funciones empleadas para la atenuación del cono de ruido son las siguientes: 

1) Low Frequency Array Filtering [11]: se fundamenta en atenuar el ruido de onda superficial, 

cambiando los datos del dominio x,t al x,f mediante la Transformada Rápida de Fourier. Las 

componentes de las frecuencias asociadas a la onda superficial, caracterizadas por presentar 

elevadas amplitudes para frecuencias bajas, son convueltas con una función tipo caja (boxcar); el 

analista determinará la longitud del arreglo para ésta función según la siguiente ecuación, para 

atenuar las componentes de las frecuencias contenidas en la banda que se desea filtrar: 

LongArreglo 


58 

V 

= 

∆x. 

f 

donde V es la velocidad de la capa superficial, ∆x el bin del proyecto (en este caso 25x25) 

y f la frecuencia. Este último parámetro corresponde a dos (2) valores de frecuencias, los cuales 

serán el umbral, mínimo y máximo, que definirán la banda de frecuencia que se desea filtrar. 

Los parámetros definidos en el proyecto, con mejores resultados después de diferentes 

pruebas, fue una velocidad de 1800 m/s y el rango de frecuencia en el que actuó la función fue 

entre 2 y 18 Hz.



Fig V.8. Efecto de la aplicación del Low Frequency Array Filtering. En el recuadro de la izquierda se distingue la 

variación en la señal asociada al ground roll (disparo con geometría en la izquierda y disparo con el filtro aplicado en 

la derecha). En el recuadro de la derecha se muestra el espectro de amplitud de los disparo analizados, observando la 

atenuación de las bajas frecuencias (0 – 18 Hz), rango en el que fue aplicado el filtro. 

El registro de la figura V.8 muestra el efecto, sobre el cono de ruido, causado por la 

aplicación de la función descrita con anterioridad. Debido a que el comportamiento de las 

reflexiones asociadas al ground roll fue modificado, como lo muestra el área enmarcada por los 

cuadros rojos, se podría atenuar éste ruido con la aplicación de futuros apilamientos. En el 

espectro de amplitud de la misma figura, se distingue la desproporción de amplitudes entre las 

frecuencias superiores a 20 Hz y las bajas frecuencias, estas últimas características del ruido en 

los datos (según la figura V.7); ante ello es necesaria la aplicación de otras funciones que resalten 

las frecuencias de las reflexiones. 

Es importante considerar que en algunos casos usar sólo este filtro no es suficiente, por la 

no uniformidad en el muestreo del ground roll, por tanto el ruido no es eliminado absolutamente. 


59



2) Time-Frequency Noise Supresión [11]: su eficiencia se basa en la eliminación del ruido sin afectar a 

las muestras y/o trazas cercanas. Dicha supresión es aplicada muestra por muestra en distintas 

componentes de frecuencias obteniendo un balanceo adecuado del espectro, considerando que 

esta función trabaja en el dominio f-t. 

La función se basa en descomponer las muestras de la traza sísmica en diferentes 

componentes de frecuencia dentro de una ventana en tiempo; emplea la Transformada Rápida de 

Fourier dentro de esta ventana y construye el espectro de amplitud de cada traza de entrada. 

El usuario define el número de muestras que contendrá la ventana de análisis y las 

frecuencias en que será aplicada la función. Esta ventana, en el dominio f,t, estará centrada en la 

muestra de la salida deseada y se aplica una función Gaussiana a todas las muestras restantes de 

la ventana. La muestra de salida resultante en el dominio f,t tendrá asociado un número 

complejo que corresponderá con diferentes sub-bandas de frecuencia. 

La supresión del ruido estará condicionada a la comparación entre las medianas de las 

amplitudes de cada sub-banda y la mediana del ruido. Ésta última corresponde con la mediana 

de todas las medianas de las amplitudes en las sub-bandas de frecuencias. Entonces, si la 

mediana de alguna sub-banda no sobrepasa el umbral determinado es inalterada; pero si el 

valor comparado es mayor al umbral, la amplitud de ésta banda será reescalada. Este proceso 

debería causar un balanceo apropiado en el espectro de amplitud. 

Los parámetros elegidos para aplicar ésta función a los datos fueron 32 trazas de análisis, 

centrando la ventana en la traza número 15; el tiempo de aplicación correspondió con dos 

ventanas (de 100 a 700 ms y de 700 a 2500 ms), el rango de frecuencias en el que actuó fue de 12 a 

100 Hz. 


60



Fig V.9. Comparación entre la aplicación de las funciones Low Frequency Array Filtering (disparo central) vs Time- 

Frequency Noise Supresión (disparo lado derecho). El espectro de amplitud, recuadro derecho, muestra la atenuación 

de las frecuencias por debajo de 20 Hz causada por la aplicación de la función Time-Frequency Noise Supresión 

En la figura V.9 se distingue el efecto causado por la aplicación de la función Time- 

Frequency Noise Supresión, donde la escala del espectro de amplitud del lado derecho de la figura 

indica la leve atenuación de las frecuencias por debajo de 20 Hz. Sin embargo, en el dominio del 

tiempo (en los registros) la atenuación de éstas frecuencias no está evidenciada con rasgos 

resaltantes. 

Con las dos posibles metodologías, descritas con anterioridad, para el acondicionamiento 

de los datos, se realizaron diferentes pruebas para determinar en que orden debían aplicarse. A 

dicho procedimiento se agregó un balanceo de amplitud consistente con superficie. 

Este balanceo es realizado a través de tres paneles distintos (Surface-Consistent Amplitude 

Balancing Analysis + Surface-Consistent Amplitude Balancing Solution + Surface-Consistent Amplitude 


61



Balancing Application). El primer módulo calcula, en ventanas de tiempo, el promedio de las 

amplitudes de cada traza, en valores absolutos, promedios pesados o rms. El siguiente módulo 

analiza la información generada por el primero y reduce la información de la amplitud 

consistente con superficie en dos fases, primero reconoce las trazas no contribuyentes y luego 

reduce la información excluyendo dichas trazas, grabando en la base de datos un factor escalar 

que describe el modelo consistente con superficie considerando la contribución de trazas con 

común receptor, fuente, cdps y/o offset. Finalmente, el tercer módulo aplica el modelo diseñado 

por el panel anterior empleando el factor escalar a los datos sísmicos. [11] 

la figura V.10. 

No obstante, el orden en que se ejecutaron las metodologías ya descritas corresponde con 

Balanceo de amplitud 

Low Frequency Array 

Filtering. 

Time Frequency 

Noise Supression 

Shot con geometría 



Fig V.10. Secuencia en la que fueron aplicadas las funciones para atenuar el ruido. 


62 


Filtering. 

1 2 3 4 5 6 


Filtering. 

Balanceo de amplitud 

Shot_Gather 



Balanceo de amplitud



El flujo con mejor resultado después de diferentes pruebas, según la figura V.10, fue el 

designado con el número 2. En las figuras V.11, V12 y V.13 se observa la atenuación del cono de 

ruido en los registros del levantamiento. 

(a) (b) 

Fig V.11. (a) Comparación disparo 9 con geometría (en la izquierda) vs disparo con atenuación del Ground Roll (en 

la derecha). (b) Comparación disparo 22 con geometría (en la izquierda) vs disparo con atenuación del Ground Roll 

(en la derecha). 

Fig V.12. Disparo con geometría, en la izquierda, y su correspondiente espectro de amplitud, en la derecha. 


63



Fig V.13. Disparo con atenuación del Ground Roll, en la izquierda, 

y su correspondiente espectro de amplitud, en la derecha. 

Sí se compara la figura V.12 con la figura V.13 se observa la apropiada atenuación del 

ground roll en los registros del levantamiento. En estas imágenes la relación entre las amplitudes 

de las frecuencias por encima a 20 Hz y las asociadas al ground roll presentan una mejor 

proporción si se compara el disparo original (fig V.12) con el disparo después de aplicar las 

funciones basadas en la atenuación del ruido (fig V.13). Sin embargo, como el espectro de 

amplitud de la figura V.13 muestra desproporción en las amplitudes de baja frecuencia respecto 

a las amplitudes de frecuencias mayores a 35 Hz y en los registros se distingue que el cono de 

ruido no fue suprimido absolutamente, en posteriores procesos se aplicarán funciones que 

permitan balancear aún más el espectro de amplitud de los registros para mejorar los resultados 

obtenidos. Por otro lado, es de considerar que los datos del levantamiento sísmico Arecuna 

Oeste 04G – 3D fueron adquiridos cercanos a la faja del Orinoco, y dicha zona se caracteriza por 

presentar baja relación señal – ruido, así como problemas de estáticas, según informes de 

procesamiento y reprocesamiento en levantamientos aledaños al área de este proyecto 3. Por ello 

se aplicaron distintas funciones para mejorar la relación señal – ruido. 

3 Informe de Procesamiento Arecuna 3D Faja Petrolífera del Orinoco, Tensor Geofísica C.A. 1997. Informe de 

Reprocesamiento Arecuna 3D. PGS Tensor Geofisica C.A. 1998. 


64



!" ! # ! $ * 

Los shot_gather, con el mejoramiento de la relación señal – ruido, obtenidos según la 

sección anterior serán procesados con la secuencia convencional para obtener los datos de 

entrada de la migración pre – apilamiento en tiempo. 

figura V.14. 

1 

2 

3 

La primera función aplicada será la Deconvolución. El flujo de trabajo se muestra en la 

Shot_gather 

Ganancia por Divergencia Esférica (1) 

Deconvolución (2) 

Edición de trazas (3) 

Shot_Deconvolución 

Fig V.14. Flujo de Trabajo para aplicar la deconvolución. 

Corrección por Divergencia esférica, para compensar la pérdida de energía por la 

propagación del frente de onda en el subsuelo. 

Aplicación de Deconvolución consistente son superficie*, fundamentada en remover el 

efecto de la ondícula del pulso, generando mejor definición de los eventos. La empleada fue 

tipo spike, debido a que proporcionó mejores resultados después de realizar varias pruebas, 

con una longitud de operador de 150 ms. Éste parámetro fue determinado con diferentes 

pruebas, considerando que valores muy elevados podían remover reflexiones de interés, pero 

valores muy pequeños causarían impulsos ruidosos. [18] 

Edición automática de trazas para eliminar aquellas trazas con amplitudes igual a cero (0). 

* Para la aplicación de la Deconvolución consistente con superficie el software emplea el análisis previo de la 

autocorrelación de los registro considerando las componentes de los receptores, fuentes y offsets. Por ello debe ser 

aplicada la función Mean Shot, Receiver Station and Offset Autocorrelation Generation previamente a la 

deconvolución. Dicha función calcula la autocorrelación consistente con superficie (considerando información de los 

disparos, receptores y offsets) y los escalares consistentes con superficie en ventanas en tiempo diseñadas por el 

usuario. En este proyecto los parámetros de cálculo fueron offset dentro de los rango de 50 a 400 m y 500 a 2800 m y 

tiempo dentro de los rangos de 100 a 1500 ms y 1000 a 3000 ms. 


65



Fig V.15. Comparación disparo 3250 sin deconvolución vs disparo 3250 con Deconvolución consistente con 

superficie aplicada. En el recuadro de la derecha se observa el espectro de amplitud, distinguiendo que el 

correspondiente al disparo con deconvolución presenta un espectro más balanceado. 

En la figura V.15 se observa el efecto causado por aplicación de la deconvolución, 

distinguiendo un registro en el que se visualizan más eventos si es comparado con el disparo del 

lado izquierdo de la misma figura, según el área encerrada por los cuadros rojos. Además, el 

espectro de amplitud de esta imagen indica el incremento en las distintas componentes de 

frecuencias del registro con deconvolución. 

Posterior a la aplicación de la Deconvolución, se continuó con la secuencia de 

procesamiento. La figura V.16 muestra las funciones aplicadas a los datos. 

Shot_Deconvolución 

Balanceo de Amplitudes Consistente con Superficie (a) 

Balanceo Espectral (b) 

Supresión de ruido (Time-Frequency Noise Supresión) (c) 

Balanceo de Amplitud en ventanas en tiempo (d) 

shot_balanceados 

Fig V.16. Flujo de trabajo para balancear el espectro de amplitud de los datos. 


66

a 

b 

c 



El balanceo de amplitud consistente con superficie es calculado mediante tres funciones 

diferentes. La primera determinará el promedio de las amplitudes en una ventana en tiempo 

diseñada por el usuario. La segunda función calculará la solución del modelo consistente con 

superficie considerando diferentes parámetros, como común offset, cdp, receptor, y fuente. Y 

el tercer módulo aplicará a los datos la solución determinada. La solución con mejor resultado 

en el proyecto corresponde a la que considera las componentes con común receptor, fuente, 

cdp y offset. Es de considerar que la solución fue calculada para tres casos diferentes, con 

común recetor y fuente; receptor, fuente y offset; y receptor, fuente, offset y cdp. Ver figuras 

V.17 a V.20. 

La función Balanceo del Espectro se fundamenta en remover el efecto del ruido generado 

por la fuente y solventar el decaimiento de las amplitudes para bandas de frecuentas de mayor 

magnitud. La función balancea todas las componentes del espectro a un mismo nivel, 

provocando que la energía total del ruido decrezca cuando las frecuencias asociadas al ruido 

son de mayor magnitud que las frecuencias de la señal sísmica de interés [11] . La función 

transforma los datos al dominio f,t mediante la Transformada Rápida de Fourier, separa la 

información en bandas de frecuencia y asigna un porcentaje a las amplitudes de dichas 

bandas, éstos parámetros son definidos por el analista. Calcula un escalar para cada banda 

según el punto medio de las mismas y de vuelta en el dominio x,t aplica la corrección 

causando la ecualización del espectro. Los convenientes resultados proporcionados por la 

aplicación de ésta función son debidos a que su efecto varía continuamente en el dominio del 

tiempo, lo cual se asocia adecuadamente con la característica particular del comportamiento 

del ruido (variar con el tiempo). Ver figura V.21. 

Se removió el ruido mediante la función Time-Frequency Noise Supresión, la cual analiza 

sub-bandas de frecuencias, calcula la mediana de las amplitudes en dicha ventana en el 

dominio f,x y si algún valor sobrepasa el umbral defino es reescalada. La aplicación de esta 

función es sugerida con el objetivo de obtener registros con mejor relación señal – ruido. Ver 

figura V.22 

d La aplicación del Balanceo de Amplitudes considerando los valores en ventanas de 

tiempo, diseñadas por el usuario, son estimadas sin analizar la consistencia con superficie. Por 

tanto, será un reescalamiento de las amplitudes considerando el máximo valor de la ventana 

analizada (Time-Variant Equalization). La función balanceará las amplitudes de la traza según 

el promedio calculado en la ventana de tiempo, éste valor será asignado a la muestra central e 

interpolado en las demás muestras de la ventana. El objetivo de aplicar esta función es 

justificado debido a que no serán datos para análisis AVO. Ver figura V.23. 


67



Fig V.17. Comparación shot_gather vs disparo con balanceo de amplitud consistente con superficie 

considerando las componentes con común fuente, receptor; fuente, receptor y offset; 

y fuente, receptor, offset y cdp, de izquierda a derecha respectivamente. 

Fig V.18. Espectro de amplitud de shot_gather (esquina superior izquierda), disparo con balanceo de amplitud 

consistente con superficie considerando común fuente y receptor (esquina superior derecha); fuente, receptor y offset 

(esquina inferior izquierda); y fuente, receptor, offset y cdp (esquina inferior derecha). 


68



El espectro de amplitud de los datos después de aplicar el balanceo de amplitud 

consistente con superficie muestra que esta función no modificó el contenido de frecuencia de 

los disparos (según la figura V.18), sólo influyó en la magnitud de las amplitudes (ver figura 

V.17). Esta apropiada influencia se distingue si se analizan los registros con más detalle, ver 

figuras V.19 y V.20. En estas imágenes se distingue la atenuación de las amplitudes del cono de 

ruido respecto a las de la señal sísmica. 

Fig V.19. Comparación de los shots_gather (esquina superior izquierda) y los disparos con balanceo de amplitud 

consistente con superficie considerando las componentes con común fuente y receptor (esquina superior derecha); 

fuente, receptor y offset (esquina inferior izquierda); y fuente, receptor, offset y cdp (esquina inferior derecha). 

Discriminando que la solución más apropiada (buena relación señal – ruido) corresponde con el balanceo consistente 

con común fuente, receptor, offset y cdp 


69



(a) (b) (c) (d) 

Fig V.20. Espectro de amplitud de la figura V.19; se distingue el mismo contenido de frecuencia 

sin (a) o con (b, c y d) la aplicación del balanceo de amplitud consistente con superficie. 

Por otro lado, el efecto causado por la aplicación del balanceo espectral es mostrado en la 

figura V.21. En ésta imagen se distingue que las amplitudes asociadas al cono de ruido son 

ecualizadas respecto a las asociadas a la señal sísmica, según el área encerrada por los recuadros 

rojos. Además, el espectro de amplitud de los registro presenta una proporción equivalente de 

amplitudes para frecuencias comprendidas en el rango de 10 a 80 Hz. 

Fig V.21. Comparación disparo con balanceo de amplitud consistente con superficie vs disparo con balanceo 

espectral. En el recuadro de la derecha se observa una adecuada nivelación de todas las frecuencias en el espectro de 

amplitud. 


70



Según el espectro de amplitud de la figura V.21, se consideró apropiada la aplicación de 

la función Time-Frequency Noise Supresión, con el objetivo de ecualizar aún más éste espectro. El 

efecto causado por la aplicación de esta función se muestra en la figura V.22; en esta imagen se 

observan amplitudes ecualizadas apropiadamente para frecuencias comprendidas en el rango 

de 15 a 80 Hz, aproximadamente. Y es de considerar que las frecuencias en las que fue aplicada 

ésta función fue de 12 a 100 Hz. Como el espectro de amplitud mostrado en la figura V.22 

presenta un adecuado balanceo de amplitudes en todo el rango de frecuencias de interés, no 

fueron aplicadas más funciones para atenuar el efecto del ground roll. 

Fig V.22. Comparación disparo con balanceo de espectro vs disparo con remoción de ruido con la función Time- 

Frequency Noise Supresión. El recuadro de la derecha muestra un espectro balanceado para todas las frecuencias. 

Finalmente, se aplicó un balanceo de amplitudes sin consistencia con superficie para 

balancear las amplitudes de la señal presente en los registros respecto a las del cono de ruido. 

Los resultados obtenidos son mostrados en la figura V.23. Esta figura permite distinguir la 

apropiada atenuación del ground roll en los registros del levantamiento, así como el espectro de 

amplitud con todas las componentes ecualizadas dentro del rango de frecuencia de interés (de 

20 a 80 Hz). 


71



Fig V.23. Comparación disparo con remoción de ruido con la función Time-Frequency Noise Supresión vs disparo 

con balanceo de amplitud sin consistencia con superficie. En el recuadro de la derecha se distingue que el espectro de 

amplitud no es afectado por la aplicación de esta función, como es de esperar. 

En la figura V.23 se observa el resultado final de la aplicación de las distintas funciones 

para atenuar el ground roll (mas no suprimir de forma absoluta), mejorando la relación señal – 

ruido. 


72



' . 3 . 

Disponiendo de los disparos con un espectro de amplitud balanceado, se procede a 

realizar el primer análisis de velocidad. Es de considerar que en esta etapa del procesamiento los 

datos son reordenados en el dominio de cdp y son aplicadas las correcciones de estáticas. 

Antes de realizar los análisis de velocidades se debe comprobar el correcto cálculo de las 

estáticas de refracción. Para ello, se comparará el apilado bruto aplicando las estáticas por 

elevación (las cuales sólo consideran el efecto de la topografía) y las estáticas de refracción. 

figura V.24. 

El cálculo de las estáticas por elevación corresponde con las siguientes fórmulas en la 

ED 

− ETp 

Para fuentes: E T = + Tuphole 

V 

Para receptores: 


73 

E 

T 

= 

r 

E 

D 

− E 

V 

r 

Tp 

Fig V.24 Ecuaciones para determinar las estáticas según la elevación. 

Posterior al cálculo de las estáticas por elevación, se realiza un análisis de velocidades 

preliminar para generar la sección apilada. Este análisis preliminar definirá el comportamiento 

regional de la velocidad en el área de estudio; por tanto, se seleccionarán velocidades en 

distancias elevadas. En el proyecto se eligieron 100 inline por 100 crossline, conformando un 

radio de 5 Km. El flujo de trabajo para realizar la selección de las velocidades y generar los 

volúmenes apilados corresponde con la figura V.25. 

Donde: ET = estática total. 

ED = elevación del datum (250 m). 

ETp = elevación topográfica. 

Vr = velocidad de reemplazo (1900m/s) 

Tuphole = tiempo de uphole.



cdp_balanceados 

Aplicación Estáticas Totales por Refracción Aplicación Estáticas Totales por Elevación 

Remoción Estáticas de LP por Refracción Remoción Estáticas de LP por Elevación 

Selección de velocidades (v_refracción) Selección de velocidades (v_elevación) 

Corrección por NMO 

Apilamiento 

Aplicación Estáticas de LP por Refracción Aplicación Estáticas de LP por Elevación 

Sección apilada con Est. por Refracción Sección apilada con Est. por Elevación 

Fig. V.25. Flujo de trabajo para la generación de un cubo aplicado corregido 

con estáticas de refracción y por elevación. *LP: largo periodo. 

Fig V.26. Sección apilada con correcciones de estáticas por refracción (lado izquierdo) y con correcciones de 

estáticas por elevación (lado derecho). La imagen permite distinguir la adecuada resolución de los eventos 

proporcionada por las correcciones de estáticas por refracción. 


74



Las secciones mostradas en la figura V.26 permiten distinguir la apropiada solución dada 

por las correcciones de estáticas de refracción en comparación con las dadas por las correcciones 

por elevación. La sección de la izquierda de ésta imagen, a la que se le aplicó las correcciones de 

refracción, presenta mejor resolución en los eventos observados en las zonas laterales de la 

sección, enmarcadas por los recuadros rojos en la figura. 

' . 3 . 

El campo de velocidades será determinado con los cdp_balanceados, seleccionando las 

velocidades que mejor horizontalizan los eventos según la máxima coherencia en la semblanza. 

La selección de estas velocidades debe ser minuciosa y elaborada en rangos de inline y crossline 

pequeños; sin embargo, no debe exagerarse en el número de cdp agrupados, debido a que se 

pueden presentar variaciones laterales en las propiedades físicas de las rocas. Por tanto, la 

definición del espaciado entre las inline, las crosssline y el número de cdp a agrupar en los cuales 

se definirán las velocidades debe ser cuidadoso y considerando los cambios estructurales y 

litológicos en el área de estudio. 

En el proyecto, la grilla que se definió para la selección de velocidades fue de cada diez 

(10) inline por cada diez (10) crossline; es decir, en un radio de 500 m. El número de cdp 

agrupados fue de nueve (9) debido a que la zona no presenta cambios estructurales de mayor 

relevancia. 

El flujo de trabajo con el que es estimado el primer análisis de velocidad, el cual 

permitirá generar el primer volumen apilado, se muestra en la figura V.27. 


75

1 

2 



Cdp_balanceados 

Filtro (4-10-40-60) (1) Aplicación Estáticas Totales por Refracción 

AGC (500) (2) Remoción Estáticas por Refracción de Largo Periodo 

Selección de Velocidades (vel_1A) Corrección por NMO (vel_1A) 


76 

Apilamiento 

Aplicación Estáticas de Largo Periodo 

Stack_vel1A 

Fig V.27. Flujo de trabajo para realizar el primer análisis de velocidades 

y generar el volumen corregido por nmo con éste análisis. 

Filtro de frecuencias por bandas, iniciando en 4 Hz y eliminando frecuencias mayores a 60 Hz. 

Corrección escalar para incrementar la amplitud de la señal de interés. La aplicación de esta 

función (AGC), ecualiza los promedios de las amplitudes a lo largo de la traza sísmica, 

reescalando las muestras a lo largo de toda la traza, según los parámetros definidos. 

Fig V.28. Selección de las primeras velocidades (vel_1A), en la izquierda. En la derecha, primera sección apilada.



La sección apilada mostrada en la figura V.28, permite distinguir la continuidad en los 

eventos proporcionada por el primer análisis de velocidades. Se espera que con los posteriores 

análisis se incremente la resolución de la imagen, diferenciando eventos con mayor continuidad. 


77 

. . . 

Con el primer volumen apilado (stack_vel1A) se determinarán las primeras estáticas 

residuales. El cálculo de estas correcciones necesita la elaboración de un cubo, denominado 

piloto, que presente buena coherencia en los eventos. 

Son varias las funciones que pueden generar el cubo piloto. La usada en el proyecto fue: 

Nombre Fundamento teórico 

3D FXY Spatial Prediction Filtering 

Filtro predictivo en 3 dimensiones, fundamentado en la 

suposición de que localmente los eventos son planos y que las 

características de la ondícula exhiben una variación espacial 

gradual. El tamaño de la ventana de aplicación y el incremento en 

el desplazamiento son designados por el usuario según ventanas 

en tiempo [11] . 

El volumen apilado calculado previamente (stack_vel1A) será el dato de entrada para la 

elaboración de las trazas piloto, con las cuales se correlacionarán las trazas originales y 

permitirán la determinación de las primeras estáticas residuales: 

Stack_vel1A 

AGC (1000) 

3D FXY Spatial Prediction Filtering 

Piloto1A 

Fig V.29. Flujo de trabajo para generar un volumen piloto.



La función proporcionada por el software Focus 5.3 para calcular las estáticas residuales 

es 3D Surface-Consistent Statics Computation. El cálculo se fundamenta en estimar la diferencia en 

tiempo consistente con superficie entre las muestras y la traza piloto, 

∆ t = s + r + S + R ( 4) 

El algoritmo para determinar esta corrección se muestra a continuación, 


Aplicación Estáticas Totales por Refracción 

Remoción Estáticas por refracción de Largo Periodo 

Corrección por NMO (vel_1A) 

Mute * 


AGC (800) 

Filtro (4-8-40-45) 

Análisis Automático de Estáticas Residuales Consistentes con Superficie 

Estáticas Residuales (según piloto1A, res1A) 

* Mute: enmudecimiento o eliminación de muestras en las trazas 

para tiempos menores al de corrección por normal moveout. 

Fig V.30. Flujo de trabajo para determinar las primeras correcciones de estáticas residuales. 

4 Relación descrita por Ronen and Claerbout (1985) and Yilmaz (1987). 


78 

∆t = variación en tiempo de las trazas, 

s = componente estática de la fuente, 

r = componente estática del receptor, 

S = componente estática estructural, 

R = componente estática residual.



Fig V.31. Comparación sección apilada corregida con el primer análisis de velocidades sin estáticas residuales 

(stack_vel1A), en la izquierda, y con correcciones residuales 1, en la derecha. La sección con las estáticas residuales 

aplicadas muestra mayor coherencia en la continuidad de los eventos. 

En la figura V.31 se muestra la comparación entre la sección con el primer análisis de 

velocidades, pero sin y con la aplicación de las primeras correcciones de estáticas residuales. En 

las zonas enmarcadas por los cuadros rojos se observa el incremento en la continuidad de los 

eventos, proporcionada por el efecto de las estáticas residuales. 

La figura V.32 muestra las correcciones de estáticas residuales calculadas en el dominio 

de disparos y receptores. La imagen permite distinguir que los rangos de valores en tiempo son 

pequeños, y las zonas con mayor variación en la corrección corresponden con los laterales del 

levantamiento 


79



Fig V.32. Mapa de residuos A, segunda corrección. Lado izquierdo para receptores, lado derecho para fuentes. 


80



' . 3 . 

La selección del segundo análisis de velocidades será realizada aplicando las primeras 

estáticas residuales. El flujo de trabajo para obtener la sección apilada se muestra en el diagrama 

de la figura V.33, paralelamente se muestra el algoritmo para la selección del segundo análisis 

de velocidades: 



Remoción Estáticas por Refracción de Largo Periodo 

Aplicación Primeras Estáticas Residuales (res1A) 

Filtro (4-10-50-60) Corrección por NMO (vel_2A) 

AGC (500) Mute 

Selección Velocidades (vel_2A) AGC (800) 


81 

Apilamiento 


Stack_vel2A 

AGC (800) 

FXY 

Piloto2A 

Fig V.33. Flujo de trabajo para realizar el segundo análisis de velocidades y generar el segundo cubo piloto. 

La verificación de la adecuada selección del segundo análisis de velocidades fue 

realizada observando el volumen stack_vel2A, en éste volumen los eventos se presentan con 

mayor resolución si se compara con el volumen stack_vel1A. Además, el cubo piloto2A será el



volumen con el que se realizará la correlación para el cálculo de las segundas estáticas 

residuales. 


82 

. . . 

Las segundas estáticas residuales son determinadas con el algoritmo de la figura V.34, las 

cuales serán aplicadas para generar el volumen apilado con el tercer análisis de velocidades. 




Aplicación Estáticas Residuales 1 (res1A) 


Mute 

Aplicación Estáticas por Refracción de Largo Periodo 

AGC (500) 

Filtro (4-8-40-50) 

Análisis Automático de Estáticas Residuales Consistentes con Superficie (según piloto2A) 

Estáticas Residuales 2 (res2A) 

Fig. V.34. Flujo de trabajo para determinar las segundas correcciones de táticas residuales. 

Seguidamente, las estáticas residuales calculadas según la figura V.34 serán aplicadas a 

los datos, generando un volumen apilado corregido con las segundas estáticas residuales. En la 

figura V.35 se muestra el algoritmo para generar este volumen.






Aplicación Estáticas Residuales 1 y 2 (res1A; res2A) 


Mute 

AGC (800) 

Apilamiento 


stack_vel2_res2A 

Fig V.35. Flujo de trabajo para aplicar las segundas correciones de estáticas residuales y generar un cubo apilado. 

Fig V.36. Comparación sección apilada, en la izquierda, con estáticas residuales 1 y primer análisis de velocidades y 

en la derecha con correcciones residuales 2 y segundo análisis de velocidades. La sección con las segundas estáticas 

residuales aplicadas muestra mayor continuidad en algunos eventos. 


83



En la figura V.36 se compara la sección apilada corregida con el primer análisis de 

velocidades y aplicadas las primeras correcciones residuales con la sección apilada corregida 

con el segundo análisis de velocidades y las segundas correcciones residuales aplicadas. Las 

zonas enmarcadas por los recuadros rojos permiten distinguir el apropiado efecto causado por 

el segundo análisis de velocidades y el incremento en la continuidad en los eventos dado por las 

segundas estáticas residuales. Sin embargo, para eventos más profundos, los cuadros azules 

delimitan una zona en la que el análisis de velocidad no resolvió adecuadamente ese evento. 

Las segundas correcciones de estáticas residuales en el dominio de fuentes y receptores 

son mostradas en la figura V.37. Esta imagen muestra valores pequeños de corrección y se 

distingue que se presentan de forma uniforme en todo la sección del levantamento; además las 

pequeñas magnitudes de éstas correcciones sugieren la culminación del cálculo de correcciones 

residuales. 

Fig V.37. Mapa de residuos, segunda corrección. Lado izquierdo para estaciones, lado derecho para fuentes. 


84



" ' . 3 . 

El tercer análisis de velocidades es deteminado según el diagrama en la figura V.38; 

paralelamente se muestra el algoritmo para generar el volumen apilado corregido con estas 

nuevas velocidades seleccionadas. 




Aplicación Primeras Estáticas Residuales (res1A) 

Aplicación Segundas Estáticas Residuales (res2A) 

Filtro (4-10-40-50) Corrección por NMO (vel_3A) 


Selección Velocidades (vel_3A) AGC (800) 


85 

Apilamiento 


Stack_vel3A 

Fig V.38. Flujo de trabajo para determinar el tercer análisis de velocidades 

y genera la sección apilada corregida por nmo con éste análisis de velocidades.



Fig V.39. Comparación sección apilada con segundo análisis de velocidades y segundas estáticas residuales (en la 

izquierda) y sección apilada con tercer análisis de velocidades y con segundas estáticas residuales aplicadas, 

stack_vel3A (en la derecha). La sección de la derecha muestra mejor continuidad en los eventos, pero sin mayor 

resolución en las estructuras, indicando la posible finalización del ciclo de cálculo de correcciones residuales. 

Las secciones apiladas de la figura V.39 muestran la comparación entre el efecto causado 

por el segundo y tercer análisis de velocidades, con segundas correcciones residuales aplicadas. 

La sección de la derecha de la figura anterior presenta eventos continuos, mostrando el 

apropiado efecto causado por las segundas correcciones de estáticas residuales, según los 

eventos encerrados por los recuadros rojos en las secciones. Sin embargo, estructuralmente no es 

mayor la diferencia entre las dos imágenes, lo que sugirió la posible finalización del ciclo de 

cálculo de correcciones residuales y análisis de velocidades. 

No obstante, con los resultados obtenidos se puede continuar con el procesamiento de 

los datos para la migración pre – apilamiento, debido a que el cubo de velocidades seleccionado 

(vel_3A) genera un volumen apilado con reflectores continuos y coherentes. 


86



Seguidamente, se aplicó la función 3D Non – Surface – Consistent Statics Computation para 

determinar pequeñas variaciones en tiempo entre una traza piloto y las trazas originales, pero 

sin considerar la consistencia con superficie. La aplicación de esta función debe ser con los datos 

referenciados al datum plano (ver diagrama en la figura V.40). 




Aplicación Primeras y Segundas Estáticas Residuales (res1A y res2A) 



3D Non – Surface – Consistent Statics Computation (estática) 

Aplicación de Estáticas no consistentes con superficie (estática) 

Remoción de Estáticas por refracción de Largo Periodo 

Cdp_trim A Apilamiento 


87 


Stack_cdptrim A 

Fig V.40. Flujo de trabajo para determinar los gathers con correcciones de estáticas según una traza piloto 

y para genera la sección apilada aplicando ésta corrección. 

El volumen stack_cdptrimA permitirá visualizar el efecto causado por la corrección de 

estáticas no consistentes con superficie, sugiriendo los apropiados resultados obtenidos si se 

aplica esta función (ver figura V.41).



Fig V.41. Comparación sección apilada con tercer análisis de velocidades y segundas estáticas residuales, 

stack_vel3A (en la izquierda) y sección apilada con tercer análisis de velocidades, segundas estáticas residuales 

aplicadas y corrección de estáticas traza a traza, stack_cdptrim A (en la derecha) 

La comparación entre la sección apilada sin y con la aplicación de las estáticas no 

consistentes con superficie mostradas en la figura V.41, permite distinguir el apropiado efecto 

causado por las estáticas determinadas a partir de una traza piloto, según lo muestran las zonas 

enmarcadas por los recuadros rojos en las secciones, en las que se observa la mayor continuidad 

en los reflectores que proporcionan las correcciones aplicadas. 

Posteriormente, se aplica una función que ecualizará todas las amplitudes (Ensemble 

Equalization) para compensar anomalías en los registros, posiblemente causadas por variaciones 

en el arreglo de la fuente, en el tamaño de la carga o profundidad [11]. La función consiste en 

designar una ventana en tiempo donde las amplitudes sean elevadas, calcula el promedio de las 


88



amplitudes dentro de la ventana y determina un escalar; este parámetro es aplicado a la muestra 

central y multiplicará a todas las muestras dentro de la ventana, por tanto se conservan las 

amplitudes relativas; es de considerar, que esta función debe ser aplicada con los datos 

referenciados al datum plano. (Ver algoritmo de trabajo en la figura V.42). 

Cdp_trim 


Ensemble Equalization 


Remoción corrección por NMO (vel_3A) 

Gathers a migrar A 

Fig V.42. Flujo de trabajo para aplicar un balanceo de amplitudes a los datos y removerles la corrección por nmo; 

generando los gathers de entrada para la migración pre – apilamiento. 

migrar A. 

Seguidamente se ejecutará la migración pre – apilamiento en tiempo con los gathers a 

Fig V.43. Cubo de velocidades de apilamiento, primer, segundo y tercer análisis, respectivamente. 


89



En la figura V.43 se muestran los cubos de velocidades generados en los tres análisis 

realizados, observando que los valores no presentan elevados gradientes, siendo una zona con 

variaciones leves, tanto laterales como verticales, de velocidades. Además, se distingue que en 

los diferentes análisis los rangos de velocidad no variaron drásticamente, debido a que para el 

mismo punto (inline 181, crossline 181, tiempo ∼ 1000 ms) la velocidad varía según la tabla V.1. 

Tabla V.1. Comparación entre los valores de Velocidad de los distintos análisis de realizados. 

Velocidad de Apilamiento(m/s) 

Primer Análisis 1995 

Segundo Análisis 1980 

Tercer Análisis 1944 


90



!' ( $ ( ! $ 5 !&& ,,* 

.!" . $ ( !' * 

Disponiendo de los cdps en datum flotante y sin corrección por normal moveout (gathers a 

migrar A), es necesario establecer algunos parámetros antes de realizar la migración pre – 

apilamiento en tiempo. Estos parámetros se describen a continuación. 

' . : se elige el rango de offset que contribuya 

apropiadamente en la migración; la determinación de éste rango estará condicionada al mapa de 

cobertura. 

La elección de éste parámetro dependerá del mínimo y máximo valor de offset que 

presente el mapa (ver figura V.44). En ésta imagen se distingue la selección del rango de offset 

con mayor contribución en el levantamiento, y la figura V.45, permite visualizar la cobertura del 

rango de offset elegido. El rango seleccionado corresponde con: 

Rango de offset seleccionado 130 – 2680 (m) 

En caso de necesitar estimar el peso de las trazas para compensar la geometría de 

adquisición, primero se debe definir el rango de azimut, el cual corresponderá con el rango de 

offset ya determinado. Posterior a ésta definición, se selecciona la opción del cálculo de pesos y 

se escoge el rango con mayor contenido de datos; dicha selección es guardada en la base de 

datos del proyecto y si se desea usar en alguna migración sólo basta seleccionar dicha opción. La 

figura V.46 muestra la selección de los pesos de las trazas; en ella se observa que este parámetro 

es elegido por rango de offset y se seleccionan los niveles con mayor contribución según el 

histograma de los pesos. 


91



Se debe considerar que al emplear el peso de las trazas como un parámetro en la 

migración los resultados pueden no ser los esperados; en algunos casos el resultado no es el más 

conveniente debido a que se incrementa el valor de la trazas con menor cobertura respecto a las 

que presentan cobertura eficiente, y en ocasiones al incrementar estos valores se agregan efectos 

no deseados en la sección migrada. Pero en otros casos, aplicar este parámetro aporta 

información importante contenida en las trazas con baja cobertura, que al ser resaltada con los 

pesos contribuye apropiadamente a los resultados finales (a nivel de apilado). 

Fig V.44. Selección del rango de offset con mayor cobertura. 


92



Fig V.45. Mapa de cobertura del rango de offset seleccionado. 

Fig V.46. Selección Pesos de las trazas. 


93




94 

* 

Este parámetro constituye la superficie a la que los datos sísmicos están referenciados, en 

caso de no estar al datum plano. Este dato es determinado mediante la remoción de las estáticas 

por refracción de largo período, referenciando la sísmica a una superficie con menor gradiente 

que la superficie topográfica en la que fueron adquiridos los datos. 

En el software Geodepth de Paradigm se establece la superficie del datum flotante; el 

software realiza el cálculo analizando las estáticas de largo período guardadas en la base de 

datos y la velocidad de reemplazo, así como la elevación de la superficie de adquisición 

(SRD = -250m). En la figura V.47 se observa el menú para calcular del datum flotante. En él se 

distingue que las estáticas de largo periodo son cargadas de la base de datos, denominadas 

cdpmean en este proyecto, y el valor de velocidad de reemplazo empleado fue 1900 (m/s). 

Es de considerar que la superficie del datum flotante puede ser determinada por el 

software de dos formas distintas, usando la velocidad de reemplazo y las estáticas de largo 

periodo (metodología empleada en este proyecto) o cargando dicha superficie definida en la 

base de datos de la sísmica, la cual podría ser definida con otro software. 

Fig V.47. Parámetros para el cálculo del Datum Flotante.



En la figura V.48 se muestra la representación de la superficie del datum flotante. Esta 

superficie es una aproximación suavizada en tiempo de la superficie de adquisición. La imagen 

permite distinguir que el comportamiento de los gradientes es muy similar al de la topografía 

(según figura IV.7 del capítulo IV, sección IV.3), pero de menor magnitud. También se distingue 

que los valores negativos corresponden con la zona central, en dirección Norte – Sur, área en la 

que se presentan las mayores elevaciones, por lo que la corrección de estáticas debe ser negativa. 

De manera análoga, los valores positivos de estáticas se asocian a los laterales del levantamiento, 

correspondientes a las zonas con menor altitud. 

Fig V.48. Superficie del Datum Flotante. 


95



4 3 

Se creó el cubo de velocidades tipo rms con los parámetros a conveniencia. Con una 

interpolación por mínima curvatura según la Ecuación de Dix (ver menú del lado izquierdo en 

la figura V.49) se generan time slice cada 120 ms, parámetro definido para este proyecto. A los 

time slice determinados se les aplica un filtro para suavizar el gradiente. El radio de filtrado en 

este proyecto para los time slice fue de 1 Km. Posteriormente, para la generación del cubo se 

presenta la opción de suavizar el valor de velocidades en dirección tanto lateral (inline y 

crossline) como vertical (por número de muestras); el radio escogido fue de 5 inline, 5 crossline 

por cada 5 muestras. El menú del lado derecho de la figura V.49 muestra las opciones 

seleccionadas para el cálculo del volumen de velocidades rms. 

Fig V.49. Parámetros para la generación de los time slice de velocidades rms (recuadro de la izquierda) 

y para los del Cubo de Velocidades rms (recuadro de la derecha). 


96



Por otro lado, considerando la transformación de velocidades de apilamiento a tipo rms, 

como la zona de estudio presenta mayormente capas horizontales la diferencia entre estos dos 

valores no debe ser grande. En la figura V.50 se muestra el cubo de velocidades rms con el que se 

ejecutará la migración pre – apilamiento en tiempo y se distingue que para el mismo punto 

(inline 181, crossline 181, tiempo ∼ 1000 ms) la diferencia no es notable entre el valor de velocidad 

de apilamiento (según la figura V.43) y tipo rms (ver tabla V.2). 

Tabla V.2. Comparación entre los valores de Velocidad de Apilamiento y Velocidad RMS. 

Velocidad de Apilamiento: Tercer Análisis 1944 (m/s) 

Velocidad tipo rms 1946 (m/s) 

Fig V.50. Cubo de velocidades rms para la migración pre – apilamiento. 

El cubo de velocidades mostrado en la figura V.50, presenta gradientes suaves tanto en la 

dirección lateral como en la vertical, sugiriendo una zona con baja complejidad estructural. 


97



$ !' * 

Según el algoritmo a utilizar se establecerán los parámetros para el cálculo de la sección 

migrada. Usando el algoritmo de Kirchhoff se definen los siguientes datos mostrados en la 

figura V.51. 

En ésta imagen se observa que se definió: 

a. La aproximación del tiempo de viaje tipo hiperbólica; 

b. La magnitud del filtro stretch media; éste parámetro esta asociado al valor del 

enmudecimiento que se empleará en el operador de migración. El valor puede variar entre 

un rango de cero (no se aplica mute) y siete (máximo valor de mute). Es de considerar que si 

se elige el mayor valor para el filtro stretch, el tiempo de migración disminuye. 

c. Un filtro anti-alias según bandas de frecuencias, con magnitud media; éste parámetro 

es empleado para evitar el fenómeno de alliasing ocurrido durante la migración, el cual es 

generado cuando el operador de migración es mayor que el espaciamiento entre las trazas de 

entrada y el contenido de frecuencia, es decir, cuando los ángulos con los que arriban las 

ondas difractadas son mayores respecto al intervalo de espaciamiento y a la máxima 

frecuencia del operador. Es de considerar que el fenómeno de alliasing en la migración 

puede introducir artefactos no deseados en la sección, principalmente, en las zonas 

someras [11]. El filtro anti-alliasing es usado para evitar la degradación en la imagen, causada 

por el efecto del filtro corte alto que se le aplica a las bajas frecuencia, lo que origina el 

incremento en los ángulos de la curva del operador dada por las trazas de entrada [11]. 

Aplicar el apropiado filtro anti-alliasing asegura convenientes resultados en la resolución de 

las imágenes, debido a que se puede atenuar el ruido provocado por el alliasing al definir 

aperturas muy elevadas para el operador de migración. 


98



d. Se migraron desde la inline 1 hasta la 321 con un paso de 10 inline, para todos las 

crossline; 

e. El valor de apertura fue de 3000 m, este dato fue el más apropiado según las 

diferentes pruebas realizadas (con valores desde 2000 m hasta 6000 m, con incrementos de 

1000 m); 

f. Los offset de entrada corresponden al rango de 1400 a 2400 m, debido a que es el que 

presenta mejor cobertura en el levantamiento; 

g. Los offset de salida corresponden con el rango de 130 a 2680 m, con un incremento de 

85 m; las magnitudes seleccionadas para éste parámetro fueron escogidas con el objetivo de 

obtener buena distribución de offset según los rangos que presenta el levantamiento. 


99



Fig V.51. Parámetros de la migración pre – apilamiento en tiempo con el algoritmo de Kirchhoff. 


100



Fig V.52. Comparación entre la sección apilada stack_cdptrimA, en la izquierda, y la sección apilada primera 

iteración de la migración pre – apilamiento en tiempo, en la derecha. La sección correspondiente a la primera 

iteración de la migración presenta eventos a tiempo 1,25 s que no son resueltos en la sección sin migrar; y la 

continuidad de los eventos puede ser mejorada con la selección de residuos de velocidades en las siguientes 

iteraciones. 

La figura V.52 muestra la comparación entre la sección apilada y la sección apilada 

migrada, respectivamente. En la imagen se distingue que el proceso de migración mejora la 

resolución de los eventos a 1.25 s, debido que en ésta zona (enmarcada por los cuadros rojos en 

la sección) se presentan leves buzamientos y están a mayor profundidad, y éstas son las 

principales característica estructurales que corrige este proceso. Con próximas iteraciones que 

permitan la selección de los residuos de velocidad se espera incrementar la continuidad en los 

eventos que el proceso de migración no muestra adecuadamente, según las zonas encerradas 

por los cuadros azules de la figura descrita. 


101



$ !' # #" # !# #. 

Desplegando los datos migrados de la primera iteración se seleccionarán los residuos de 

velocidad según la máxima semblanza, de tal forma que el moveout de las reflexiones sea 

corregido. La selección de los residuos debe realizarse en las mismas inline y crossline en las que 

fueron seleccionados los análisis de velocidades. En el segundo panel de izquierda a derecha de 

la figura V.53 se observa la desviación de la máxima semblanza respecto a la línea de tendencia 

que representa a la apropiada corrección del moveout; seleccionando estos residuos de velocidad 

se podrá corregir la continuidad de los eventos. 

Posterior a esta selección, se crea el cubo de velocidades residuales para obtener el nuevo 

volumen apilado corregido con los residuos seleccionados. Este volumen deberá presentar 

eventos con mejor definición. 

Fig V.53. Selección de residuos de velocidad para la primera iteración. 


102



Una vez seleccionados los residuos se procede a generar el cubo de velocidades. Primero, 

se obtendrán slices de las velocidades iniciales modificadas según los residuos. El menú para 

generar estas secciones se muestra en el recuadro de la izquierda de la figura V.54; luego se 

interpola por mínima curvatura para generar el volumen de velocidades completo, según el 

menú mostrado en el recuadro de la derecha de la figura V.54. 

Fig V.54. Parámetros para la generación del cubo de residuos de velocidades. En la izquierda se muestra el menú 

para crear los slice de los residuos de velocidades y en la derecha el menú para la generación del cubo. 


103



" # $ !' * 

Con la selección de los residuos y el cubo de velocidad generado a partir de los mismos, 

se ejecuta una segunda iteración sin variar los parámetros para el cálculo de la migración. 

Luego, los datos de la segunda iteración, de igual forma como la descrita en la sección V.2.2, son 

analizados. Este análisis sugirió la selección de nuevos residuos, debido a que se observaron 

desviaciones de la máxima semblanzan respecto a la línea de tendencia. Con los segundos 

residuos seleccionados se generó el cubo de velocidades para una nueva iteración, éste es 

mostrado en la figura V.55. En ella se distingue el bajo gradiente, tanto lateral como vertical, que 

presentan las velocidades en la zona de estudio. 

Fig V.55. Cubo de velocidades rms, determinado con los segundos residuos de velocidades. 


104



! $ !' * 

Con el cubo de velocidades generado con las segundas correcciones de residuos se 

realiza la tercera iteración de la migración, sin variar los parámetros descritos en la sección V.2.2. 

En la imagen de la figura V.56 se muestra la sección apilada de la tercera iteración de la 

migración pre – apilamiento en tiempo. La culminación de la realización de posteriores 

iteraciones es sugerida con el despliegue de los gathers migrados por la tercera iteración y la 

observación de residuos de velocidades próximos a la línea de tendencia, por tanto no se 

consideró necesaria la selección de nuevos residuos de corrección. 

Fig V.56. Comparación sección apilada stack_cdptrimA y sección apilada de la tercera iteración de la migración, 

inline 121, respectivamente. En la sección de la derecha se distingue mejor resolución de los eventos profundos, 

posiblemente por ser la zona con mayor buzamiento de la sección y es en dicha área donde la migración arroja 

mejores resultados. 


105



Fig V.57. Comparación sección apilada stack_cdptrimA y sección apilada de la tercera iteración de la migración, 

inline 81, respectivamente. En la sección de la derecha se distingue mejor resolución de los eventos profundos de 

igual forma que en la figura V.56. 

En las figuras V.56 y V.57 se muestra la comparación entre las secciones apiladas y las 

secciones apiladas migradas con tres (3) iteraciones obtenidas a partir de los datos gathers a 

migrar A. En las zonas enmarcadas por los cuadros rojos se distingue la mayor continuidad de 

los eventos que proporciona el proceso de migración, y en las zonas delimitadas por los cuadros 

azules, se observa que la resolución estructural que ofrece este proceso es apropiada. Sin 

embargo, en las zonas encerradas por los recuadros fucsia se visualiza que la migración no 

arrojó resultados apropiados, debido a que en esas secciones se perdió la continuidad en los 

eventos lograda en las secciones apiladas; este hecho quizás pueda ser atribuido a la errónea 

selección de velocidades en los análisis de residuos, lo que sugiere que se debe ser más 

cuidadoso y detallado en esta fase del proyecto, ya que la migración pre – apilamiento es muy 

sensible a la selección de velocidades. 


106



!' $6 ( $ ( ! $ # 5 !&& ,,* 

La migración post – apilamiento en tiempo con el algoritmo de Kirchhoff fue ejecutada 

con el objetivo de comparar la resolución dada por la migración pre – apilamiento en tiempo. El 

cubo de entrada para esta migración fue la sección apilada stack_cdptrimA con un post – 

proceso, según la figura V.58, para incrementar la continuidad y coherencia de los eventos en el 

volumen. 

Stack_cdptrimA 

AGC (800) 

FXY 

Stack_post_migración 

Fig V58. Flujo de trabajo para generar el cubo apilado para la migración post – apilamiento. 

La figura V.59 muestra el menú para el cálculo del cubo migrado post – apilamiento; en 

ella se observa que se migraron todas las inline, con todas las correspondientes crossline, al igual 

que en la migración pre – apilamiento se aplicó una magnitud media de filtro strech, un filtro 

anti – alias según bandas de frecuencias y la misma apertura de 3000 m. Es de considerar, que 

los parámetros definidos en esta migración corresponden con los mismos fundamentos teóricos 

expuestos en la sección V.2.2. 


107



Fig V.59. Menú para el cálculo de la migración post – apilamiento en tiempo con el algoritmo de Kirchhoff. 

Fig V.60. Comparación stack_cdptrimA, stack_cdptrimA con post – proceso (stack_post_migración) 

y migración post – apilamiento en tiempo, de izquierda a derecha respectivamente. 


108



En las imágenes de la figura V.60 se observa que en las zonas enmarcadas por los 

recuadros rojos a medida que son aplicados los diferentes procesos (post – procesos y 

migración), la continuidad y coherencia de los eventos incrementan apropiadamente, generando 

volúmenes que permiten interpretaciones estructurales sencillas de analizar. 

Sin embargo, si se compara la sección migrada post – apilamiento con una sección 

migrada pre – apilamiento (ver figura V.61), se observa una sección con mayor continuidad y 

coherencia en los reflectores para la migración post – apilamiento. Ante este hecho, se sugirió 

verificar la metodología con que fueron generados los datos de entrada de la migración pre – 

apilamiento en tiempo. Procesamientos previos en zonas cercanas al levantamiento sísmico 

Arecuna Oeste 04G – 3D sugieren que las correcciones de estáticas son un factor determinante 

en la calidad de las secciones apiladas en el área de estudio. Se propuso seleccionar los primeros 

arribos pero sin aproximación automática, sino seleccionando las primeras llegadas en todos los 

disparos del levantamiento. 


109



Fig V.61. Comparación sección apilada de migración pre – apilamiento en tiempo 

y migración post – apilamiento en tiempo. 

Por otro lado, se podría sugerir que la sección obtenida con la migración post 

apilamiento en tiempo puede ser considerada como un complemento de la sección migrada pre 

– apilamiento en tiempo según figura V.62. En ésta figura, se distingue que es mayor la 

continuidad en los eventos de la sección central de la imagen comparada con la sección de la 

derecha. Sin embargo, los eventos con buzamientos notables, presentes por debajo de los 1,2 s en 

las secciones (áreas delimitadas por los recuadros fucsia), son desplazados en tiempo y espacio, 

debido al efecto causado por el proceso de migración sísmica. En la sección migrada pre – 

apilamiento éstos eventos se distinguen con una definición más apropiada que no es observada 

en la sección migrada post – apilamiento, permitiendo generar una imagen con reflectores que 

posiblemente intentan ser la mejor representación de los eventos geológicos en el subsuelo, lo 


110



que podría hacer la diferencia para obtener mayor aproximación y certeza en próximas 

perforaciones que sean ejecutadas según lo sugieran interpretaciones hechas con este volumen 

migrado pre – apilamiento en tiempo. 

Fig V.62. Comparación sección apilada, sección migrada post apilamiento y sección migrada pre apilamiento, en el 

dominio de tiempo, de izquierda a derecha respectivamente. 


111



!' # # $ ! " % ! !! # $.$! 7 3 !!' 

# ( # $ # # ( # % $ $ * 

.!" # $.$! # , !!' * 

Las correcciones estáticas por refracción son estimadas con la selección de los primeros 

arribos en todos los shots del levantamiento, para un total de 5509 shots. Como fue explicado en 

la sección V.1.1, el programa proporciona las estáticas por refracción totales para fuentes y 

receptores, permitiendo separarlas en el dominio de las de periodo corto y largo. Estas 

correcciones estáticas permitirán referenciar los datos sísmicos a una superficie de menor 

gradiente que la superficie de adquisición, además de remover el efecto de la capa meteorizada, 

debido a que se distorsiona la continuidad y coherencia en la hipérbola que representa las 

reflexiones de interés. 

Luego de seleccionar los primeros quiebres en todos los disparos del levantamiento 

(5509), de igual forma que en la sección V.1.1, si se agrupan varios disparos en una sección 

tiempo vs offset se observan dos tendencias distintas en las primeras llegadas; pero no son 

asociadas a la mala selección de la fase del arribo, como lo demuestra la figura V.63. En ésta 

imagen se observa la correcta selección del primer quiebre, los valores de las primeras llegadas 

en los disparos analizados, para el shot 4535 corresponden aproximadamente con 0.10 ms; 

mientras que en el shot 4536 corresponde con 0.13 ms. Lo que sugiere que la diferencia en las 

líneas de tendencia de las domocrónicas está posiblemente asociada a variaciones litológicas y/o 

propiedades físicas intrínsecas de la zona analizada. 


112



0.10 s 

Fig. V.63. Selección primeras llegadas para los disparos 4535 y 4536. 

Después de la selección de las primeras llegadas se emplea la función Generalized 3D 

Refraction Statics Generation del software Focus 5.3 y se obtiene la solución del modelo de la capa 

meteorizada. Dicho modelo fue generado considerando un solo refractor con velocidad entre 

1700 y 1900 m con contribución de offset entre 350 y 2800 m. 

s 

Posteriormente, se aplicó la función CDP Mean Statics Computation, la cual determina las 

correcciones estáticas para referenciar los datos a una misma superficie. Los datos de entrada 

para aplicar esta función deben estar en el dominio de cdp. La solución será dos correcciones 

estáticas, una de menor magnitud que referencia los datos al datum flotante, provocando la 

disminución en la distorsión de hipérbola que representa las reflexiones de los eventos, y la otra 

corrección estática, de mayor componente en el dominio de fuente y receptor, referencia los 

datos al datum final. Las siguientes figuras (V.64 y V.65) muestran las estáticas totales en el 

dominio de receptores, fuentes y cdp. 


113 

0.13 s



Fig V.64. Correcciones de estáticas totales, en la izquierda en el dominio de disparos, en un rango de -40 a 58 ms; 

y en la derecha de receptores, en un rango de -53 a 43 ms. 

Fig V.65. Correcciones de estáticas de refracción B en el domino de cdp. Oscilan en el rango de -91 a 95,5 ms. 


114



De igual forma que en la sección V.1.1, la tendencia de las correcciones de estáticas de 

refracción en el dominio de receptores, fuentes y cdp corresponde adecuadamente, según su 

carácter de correcciones de estáticas, con el comportamiento de la topografía. Las figuras V.64 y 

V.65 permiten distinguir que para la zona central, con mayor elevación topográfica, las 

correcciones de estáticas serán de menor valor (negativas); mientras que en los laterales del 

levantamiento, con menor elevación, las correcciones serán de mayor magnitud (positivas). 

Comparando las nuevas correcciones de estáticas de refracción con las calculas en la 

sección V.1.1 se observa la similitud en los rangos de magnitud (ver tabla V.3). 

#" * 

Tabla V.3. Comparación de Estáticas de Refracción según la selección de primeros quiebres 

en cada diez disparos y en todos los disparos del Levantamiento. 

Estáticas por Refracción Selección cada 10 shots Selección todos los shots 

Totales por Fuente -39,1 a 57,4 ms -40,0 a 58,0 ms 

Totales por Receptor -51,7 a 44,4 ms -53,0 a 43,0 ms 

Totales por cdp -86,0 a 96,0 ms -91,0 a 95,5 ms 

. # !# # 2 !' # !! $.$! 

' . 3 

La selección de velocidades será generada con los mismos datos (shots_balanceados) a los 

que se les aplicó la secuencia de acondicionamiento en la sección V.1.2. Es de considerar que con 

las nuevas correcciones de estáticas calculadas no se realizó el apilado bruto para comparar 

dichas correcciones con la solución de las estáticas por elevación, debido a que las magnitudes 

de estas últimas correcciones son muy similares a las correcciones de estáticas por refracción 

calculadas con anterioridad (según la tabla V.3) 


115



El algoritmo para seleccionar las primeras velocidades se muestra en la figura V.66: 


Filtro (4-10-40-60) Aplicación Estáticas Totales por Refracción 

AGC (500) Remoción Estáticas de Largo Periodo 

Selección de Velocidades (vel_1B) Corrección por NMO (vel_1B) 


116 

Apilamiento 


Stack_vel1B 

AGC (800) 

FXY 

Piloto1B 

Fig V.66. Flujo de trabajo para realizar el primer análisis de velocidades y el primero cubo piloto. 

Fig V.67. Sección apilada con las correcciones de estáticas manuales (stack_vel1B), lado izquierdo. Sección apilada 

con las correcciones de estáticas automáticas (stack_vel1A), lado derecho. La imagen no muestra mucha diferencia 

en la resolución de los eventos proporcionada por la solución manual en comparación con la solución automática.



La imagen de la figura V.67 muestra la comparación entre la sección apilada con las 

nuevas correcciones de estáticas aplicadas y la sección apilada con las correcciones de estáticas 

calculadas en la sección V.1.1. En rasgos generales, las dos secciones presentan la misma 

continuidad en los eventos, sin mostrar mayor diferencia en la resolución dada por las secciones. 

. . 

Para el cálculo de las primeras correcciones residuales se emplea el Pilto1B como cubo 

modelo para la autocorrelación que permite estimar dichas correcciones. La figura V.68 muestra 

el flujo de trabajo para determinar éstas correcciones. Y en la figura V.69 se muestra el resultado 

de aplicar las correcciones residuales calculadas. 



Remoción Estáticas de Largo Periodo 

Corrección por NMO (vel_1B) 

Mute 


AGC (800) 

Filtro (4-8-40-50) 

Análisis Automático de Estáticas Residuales Consistentes con Superficie 

Estáticas Residuales (según piloto1B, res1B) 

Fig V.68. Flujo de trabajo para determinar las primeras correcciones de estáticas residuales. 


117



Fig V.69. Comparación sección apilada corregida con vel_1A y aplicadas las correcciones residuales res1A, lado 

izquierdo; con la sección apilada corregida con vel_1B y aplicadas las correcciones residuales res1B, lado derecho. 

En la imagen no se distingue mayor diferencia entre la sección de la derecha y la de la izquierda, no se observan 

otras estructuras o eventos distintos en la sección de la derecha como se hubiese querido, dando el primer indicio de 

que las primeras correcciones estáticas dan buena resolución a los datos del levantamiento procesado. 

La figura V.69 permite observar la comparación entre las secciones apiladas corregidas 

con el primer análisis de velocidades, aplicadas las primeras correcciones residuales y 

corregidas con las estáticas de refracción determinadas en la sección V.1.1 y en la sección V.3.1, 

respectivamente. Las zonas limitadas por los cuadros rojos en la figura muestran mayor 

continuidad en los eventos de la sección corregida con las últimas estáticas de refracción. Sin 

embargo, la resolución de la sección en el lado derecho de la figura V.69 no es tan apropiada si 

es comparada con la sección corregida con las estáticas calculas en la sección V.1.1, debido a que 

no se presentan nuevas estructuras o eventos distintos como se hubiese esperado. 

En la figura V.70 se muestran las primeras correcciones de estáticas de refracción en el 

dominio de receptores y fuentes. La imagen permite distinguir los pequeños valores dados por 

estas correcciones, sin presentar un comportamiento característico en el área de estudio. 


118



Fig V.70. Mapa de residuos, primera corrección. Lado izquierdo para estaciones, lado derecho para fuentes. 


119



' . 3 

Disponiendo de las primeras correcciones residuales se puede estimar el segundo 

análisis de velocidades, para determinar una sección con mayor resolución. La figura V.71 

muestra el algoritmo de trabajo para determinar el segundo análisis de velocidades y, 

paralelamente, muestra el flujo para generar el segundo volumen piloto. 




Aplicación Primeras Estáticas Residuales (res1B) 

Filtro (4-10-50-60) Corrección por NMO (vel_2B) 


Selección Velocidades (vel_2B) AGC (800) 


120 

Apilamiento 


Stack_vel2B 

AGC (800) 

FXY 

Piloto2B 

Fig V.71. Flujo de trabajo para realizar el segundo análisis de velocidades y generar el segundo cubo piloto.

figura V.72. 



. . 

Seguidamente, se calculan las segundas estáticas residuales, según el algoritmo de la 




Aplicación Estáticas Residuales 1 (res1B) 


Mute 


AGC (500) 

Filtro (4-8-40-50) 

Análisis Automático de Estáticas Residuales Consistentes con Superficie (según piloto2B) 

Estáticas Residuales 2 (res2B) 

Fig V.72. Flujo de trabajo para determinar las segundas correcciones de estáticas residuales. 

Las segundas estáticas residuales (res2B) serán aplicadas (ver figura V.72) a los datos 

para aumentar la coherencia y continuidad en los eventos de la sección. 


121



Fig V.73. Comparación entre sección apilada corregida con segundo análisis de velocidades y aplicadas segundas 

correcciones residuales, pero corregidas con estáticas de refracción diferentes. La sección de la izquierda esta 

corregida con las estáticas de refracción calculadas en la sección V.3.1, con las segundas correcciones residuales 

(res2B) y por nmo con vel_2B. Mientras que, la sección de la derecha esta corregida con las estáticas de refracción 

estimadas automáticamente, con segundas correcciones residuales (res2A) y por nmo con vel_2A. 

En la figura V.73 se observa la diferencia entre las secciones apiladas corregidas con el 

segundo análisis de velocidades, segundas correcciones residuales y las estáticas de refracción 

calculas en la sección V.1.1 y en la V.3.1. Las áreas enmarcadas por los recuadros rojos muestran 

la apropiada resolución dada por las nuevas correcciones de estáticas de refracción. Es de 

considerar, que aunque las magnitudes de las correcciones de estáticas calculadas en la sección 

V.1.1 y V.3.1 presentan sólo leves variaciones entre sí (ver tabla V.3), la diferencia en la 

continuidad de los eventos observada en la figura V.73 sugiere la severa influencia de las 

correcciones de estáticas en el área de estudio. 


122



En la figura V.74 se muestran las segundas correcciones de estáticas residuales en el 

dominio de receptores y disparos. En la imagen se observa que el rango de estas correcciones es 

de menor magnitud que el de las primeras correcciones de estáticas residuales; también se 

distingue que no se presentan con un comportamiento uniforme en toda la zona, sino como 

pequeñas correcciones locales. 

Fig V.74. Mapa de residuos, segunda corrección. Lado izquierdo para receptores, lado derecho para fuentes. 

Por otro lado, comparando el rango en que oscilan las correcciones residuales dadas con 

las dos metodologías (determinando estáticas de refracción según primeros quiebres 

seleccionados cada diez shots o en todos los shots) se observa que los valores pertenecen a rangos 

muy similares (ver tabla V.4). 

Primera Corrección 

Segunda Corrección 

Tabla V.4. Comparación entre los rangos de las Correcciones de Estáticas Residuales 

Estáticas Residuales Selección cada 10 shots Selección todos los shots 


123 


Totales por Receptor -5,2 a 4,8 ms -5,1 a 3,5 ms 


Totales por Receptor -6,3 a 3,6 ms -7,3 a 3,8 ms



La tabla V.4 indica la pequeña variación en las correcciones residuales dadas con la 

nueva selección de los primeros quiebres. 

" ' . 3 

Finalmente, para ser equivalente la comparación con las secciones apiladas de la sección 

V.1.3.1.5, se realiza el tercer análisis de velocidades y la sección apilada corregida con dichas 

velocidades es obtenida con el diagrama mostrado en la figura V.75. 




Aplicación Primeras Estáticas Residuales (res1B) 

Aplicación Segundas Estáticas Residuales (res2B) 

Filtro (4-10-40-50) Corrección por NMO (vel_3B) 


Selección Velocidades (vel_3B) AGC (800) 


124 

Apilamiento 


Stack_vel3B 

Fig V.75. Flujo de trabajo para realizar el segundo análisis de velocidades y generar el segundo cubo piloto.



Fig V.76. Comparación entre las secciones apiladas corregidas con el tercer análisis de velocidades pero con 

diferentes correcciones de estáticas de refracción aplicadas. La sección de la izquierda esta corregida con las estáticas 

por refracción según todos los shots, con las estáticas residuales res2B y por nmo con vel_3B. La sección de la 

derecha esta corregida con segundas correcciones residuales (res2A) y por nmo con vel_3A. Las imágenes permiten 

distinguir que la diferencia entre las dos secciones es leve. 

En la figura V.76 la sección apilada del lado izquierdo no muestra gran diferencia en la 

continuidad de los eventos al ser corregidos con el tercer análisis de velocidades y las segundas 

correcciones residuales, como muestra la zona limitada por los recuadros rojos, permitiendo 

finalizar con el ciclo del cálculo de estáticas residuales y los análisis de velocidades respectivos. 

A nivel de secciones apiladas los resultados con la aplicación de estas nuevas 

correcciones son similares, debido a que no se observan estructuras resaltantes como era de 

esperar, lo que podría sugerir que las secciones migradas con estas últimas correcciones de 

estáticas presentarán una apariencia similar que las obtenidas en la sección V.1.3.1.5. 


125



Posteriormente, se aplicó la función 3D Non – Surface – Consistent Statics Computation, la 

cual determina variaciones en tiempo de una traza piloto con las trazas originales, pero sin 

considerar la consistencia con superficie, generando correcciones de estáticas para cada traza 

que serán aplicadas con el objetivo de incrementar la coherencia en los eventos. El algoritmo 

para determinara éstas correcciones se muestra en la figura V.77. 




Aplicación Primeras y Segundas Estáticas Residuales (res1B y res2B) 



3D Non – Surface – Consistent Statics Computation (estáticaB) 

Aplicación de Estáticas no consistentes con superficie (estáticaB) 

Remoción de Estáticas de Largo Periodo 

Cdp_trim B Apilamiento 


126 


Stack_cdptrim B 

Fig V.77. Flujo de trabajo para generar los datos corregidos con estáticas según una traza piloto 

y para generar un volumen apilado aplicando éstas correcciones. 

El volumen stack_cdptrimB permite diferenciar el efecto de la corrección de estáticas no 

consistentes con superficie, sugiriendo el beneficio de aplicar esta función (ver la figura V.78). En 

esta imagen, como lo muestran las áreas enmarcadas por los recuadros rojos, la continuidad de 

los eventos es incrementada con la aplicación de la función descrita con anterioridad.



Fig V.78. Comparación sección apilada con tercer análisis de velocidades y segundas estáticas residuales, 

stack_vel3B (en la izquierda) y sección apilada con tercer análisis de velocidades, segundas estáticas residuales 

aplicadas y corrección de estáticas traza a traza, stack_cdptrim B (en la derecha) 

Finalmente, si la sección apilada corregida con las estáticas de refracción determinadas 

con la selección de los primeros quiebres para todos los disparos (manualmente) es comparada 

con la sección apilada corregida con las estáticas de refracción obtenida con la selección de los 

primeros quiebres en cada diez (10) disparos, se distingue que la diferencia en la resolución y 

continuidad de los eventos no es de gran relevancia (ver figura V.79). En ésta figura se observa 

que la diferencia en la resolución de los eventos en las dos secciones no es de gran magnitud; 

además, con la aplicación de las nuevas estáticas de refracción no se generaron nuevas 

estructuras como era de esperar. Por ello, no se consideró necesaria la realización de la 

migración pre – apilamiento en tiempo con los nuevos gathers obtenidos (cdp_trimB). 


127



Fig V.79. Comparación entre la sección apilada corregida con segundas estáticas residuales y tercer análisis de 

velocidad. En la izquierda con estáticas de refracción según primeros quiebres seleccionados manualmente en todos 

los disparos, y en la derecha con estáticas de refracción determinadas con primeros quiebres seleccionados cada diez 

disparos y el resto de forma automática. En la imagen no se observa gran diferencia en la resolución de los eventos. 

En las figuras V.76, V.78 y V.79 se observa que la aplicación de las nuevas correcciones 

de estáticas de refracción, calculadas según la sección V.3.1, no otorgan mayor resolución en la 

líneas apiladas si se compara con las líneas apiladas corregidas con las estáticas de refracción 

calculadas en la sección V.1.1, debido a que no se distinguen nuevos eventos o reflectores con 

mayor continuidad y coherencia y/o diferentes estructuras geológicas como se hubiese 

esperado, lo que permite sugerir la apropiada resolución dada por las correcciones de estáticas 

de refracción calculadas a partir de la selección de los primeros quiebres en cada diez (10) 

disparos e interpolando, por la metodología de redes neuronales, en los disparos no 

seleccionados. Además, se generó un cubo sísmico migrado pre – apilamiento en el dominio de 

tiempo que permite una interpretación estructural del área, y ésta podría ser complementada si 


128



se incluye la información dada por el cubo migrado post – apilamiento en tiempo. Así, esta 

interpretación podría ser la primera proposición para promover futuras perforaciones en la zona 

delimitada por el levantamiento sísmico Arecuna Oeste 04G – 3D. 


129



1.) La atenuación del ground roll en los registros del levantamiento necesitó la aplicación de 

una secuencia de funciones que permitieron filtrar el efecto del cono de ruido en los 

datos procesados. 

2.) La ecualización del espectro de los registros fue aplicada con el objetivo de compensar 

todas las componentes de las frecuencias de la señal para disminuir el efecto causado por 

frecuencias anómalas, tales como las generadas por los ruidos de los disparos. 

3.) Las amplitudes de los datos fueron balanceadas considerando la consistencia con 

superficie, tomando en cuenta la influencia de las componentes con común receptor, 

fuente, offset y cdp; y considerando el comportamiento de las reflexiones según ventanas 

en tiempo y rangos de offset. La aplicación de estas funciones no influye en el propósito 

de los próximos análisis que se realizarán con éstos datos, debido a que los mismos no 

serán empleados para estudios AVO, sino para interpretaciones estructurales. 

4.) La corrección de estáticas de refracción realizada según la selección de los primeros 

arribos en todos los disparos y en cada diez (10) disparos, interpolando por redes 

neuronales las primeras llegadas en los disparos no seleccionados, generaron secciones 

apiladas muy similares. Esto sugiere que para los datos del levantamiento procesado la 

metodología de interpolación por redes neuronales proporcionó una solución confiable, 

aún cuando se presumía que la zona ameritaba la aplicación de una minuciosa corrección 

de estáticas por refracción. 

5.) Los análisis de velocidad no fueron difíciles de realizar, debido a que la zona estudiada 

no es un área con elevada complejidad estructural, por tanto las velocidades no 


130



presentaron elevados gradientes. Es de considerar que en algunas áreas fue confusa la 

selección de la velocidad que permitía resolver el basamento, debido a que en dichas 

zonas éste presentaba ondulaciones causadas por paleoerosiones. 

6.) Las imágenes obtenidas con la migración post – apilamiento en tiempo con el algoritmo 

de Kirchhoff presentaron eventos con mejor continuidad que en las secciones migradas 

pre – apilamiento en tiempo. Esta última metodología no proporcionó una imagen 

apropiada del subsuelo como era de esperar, lo que podría asociarse posiblemente a la 

sensibilidad de la migración pre – apilamiento a las pequeñas variaciones en las 

velocidades que no fueron detectadas en los tres (3) análisis de velocidades y/o a la baja 

relación señal – ruido que presenta la zona estudiada. Sin embargo, el volumen generado 

con la migración pre – apilamiento en tiempo presenta eventos buzantes profundos 

(aproximadamente a los 1,2 s en el volumen) con buena definición que pueden ser 

analizados en futuras interpretaciones estructurales para promover perforaciones 

productivas. 

7.) Por la gran influencia de las correcciones de estáticas en la resolución de los eventos en 

una imagen sísmica, se sugiere que previo al cálculo de las correcciones de estáticas de 

refracción, se analice detalladamente los datos de refracción de campo y generar un 

modelo de la capa meteorizada para comparar y corroborar el adecuado cálculo de las 

correcciones de estáticas calculadas a partir de la selección de las primeras llegadas. 

8.) Para el cálculo de la migración pre – apilamiento en tiempo se puede considerar la 

opción de usar el peso de las trazas, lo cual permite compensar la geometría de 

adquisición del levantamiento. En algunos casos, como en este proyecto, al considerar 

esta opción se agregaron efectos no deseados en la sección apilada. Sin embargo, se 

recomienda realizar siempre la prueba de comparación de las secciones migradas 


131



considerando o no los pesos de las trazas, debido a que se pueden presentar zonas en las 

que las trazas con baja cobertura quizás contener información de gran interés, y al ser 

incrementados sus pesos agregarán dicha información en la sección final. 

9.) Interpretaciones previas en el levantamiento Arecuna, el cual colinda por el Este con el 

levantamiento Arecuna Oeste 04G – 3D, sugirieron el estudio del levantamiento 

procesado en este proyecto. Debido a que se identificaron zonas de interés exploratorio 

que se extienden hacia ésta área, se recomienda la realización de la migración en 

profundidad del levantamiento, calibrar ese volumen con futuros registros de pozos y 

realizar la interpretación petrofísica que permita discriminar y caracterizar posibles 

estructuras con potencial económico de interés presentes en la zona. 


132




133 

(( ' 

[1] Adan, E. y Mikereit, B. A fast method for accurately determining 3-D refraction static 

corrections. University of Toronto. Toronto, USA. 

[2] Alan Lembang. Practical Aspect of Seiemic Data Adquisition. Bucharest. Marzo, 2003. 

[3] Ancheta M., Arquímedes J. Evaluación de las Estrategias de Explotación del Yacimiento 

U1/HAM-3, Campo Arecuna, Área Hamaca, Faja del Orinoco, Durante la Fase de Agotamiento 

Natural. Tesis de Grado. Universidad de Oriente, marzo de 2005. 

[4] Biondo, L., Biondi. 3-D Seismic Imaging. Stanford University. Septiembre, 2004. 

[5] BGP Internacional of Venezuela, S.A.. Volumen I: Operaciones Levantamiento Sísmico 

Arecuna Oeste 04G – 3D. Adquisición de Datos Sísmicos. Reporte Final. Estado Anzoátegui, 

Venezuela. Marzo – Mayo de 2005. 

[6] Cavada, José. Guía de Prospección Sísmica por Refracción. Facultad de Ingeniería de la 

Universidad Central de Venezuela. Versión 1.4. Marzo de 2000. 

[7] Claerbout, J.F., Imaging the Earth’s Interior. Geophysics Department Stanford University. 

California, USA.1984.



[8] Geo Exploración, S.A. Revisión Geológica del Área Central y Nor – Central del 

Cuadrángulo Arecuna. Julio, 1998. 

[9] Gresillaud,A., Suaudeau, E., Granger, P.. C-wave imaging : Kirchhoff PSTM versus DMO. 

CGG, R&D Department, Francia. 

[10] Kim, Y.C., Gonzalez, R., Migration Velocity Analysis whit The Kirchhoff Integral. 

Geophysics, 56, pp. 365 -370. 1991. 

[11] Manual de Ayuda del software Focus 5.3 de Paradigm 

[12] Manzoni, Salvador. Secuencia Básica de Procesamiento de Data Sísmica para Obtener 

Datos Finales, Orientada a Estudiantes de Ingeniería Geofísica. Tesis de Grado. Universidad 

Central de Venezuela, julio de 2001. 

[13] Pacheco B, Carlos F. Evaluación de la Sensibilidad a la Estimación de Velocidades de las 

Técnicas de Migración en Zonas Complejas. Tesis de Grado. Universidad Simón Bolívar. 

Noviembre de 1997. 

[14] PDVSA Intevep, PDVSA Exploración y Producción. Léxico Estratigráfico Electrónico de 

Venezuela. Última Actualización en marzo de 1999. www.pdv.com/lexico 

[15] Pereira, J., del Pino, E. y otros. Curso Básico de Procesamiento Sísmico de Reflexión. 

Publicado por Lagoven, filial de Petróleos de Venezuela. 1990. 


134



[16] Regueiro, José. Sísmica para Ingenieros y Geólogos de Yacimiento. CIED. 1997. 

[17] Schneider, W., Integral Formulation for Migration in two and three Dimensions. 

Geophysics, 43, pp 49 - 76. 1978. 

[18] Yilmaz, O., Seismic Data Processing. Society of Explorations Geophysicists, Tulsa, USA. 

Segunda Edición. 

[19] Yilmaz, O., Chambres, R.. Migration Velocity Analysis by Wavefiekd Extrapolation. 

Geophysics, 49, pp. 1664 . 1674. 1984. 


135



* # + " 


136 

8 

Constituye el fundamento teórico para el análisis de señales, continuas o discretas. “En el 

se establece que cualquier función periódica de cuadrado integrable puede ser expresada como 

la suma pesada (combinación lineal o superposición) de un número, en principio infinito, de 

funciones trigonométricas cuyos períodos dependen en una forma precisa del período de la 

función original”. 5 

. +" !' # " # # $ 9 * 

2 

f ∈ L ) en 

Sea una función de variable real y valores reales, es de cuadrado integrable ( [ a,b] 

el intervalo [ a, b] 

, si f satisface la condición: 

. # + " * 

Toda función periódica ( t T ) ( t ) 

b 

a 

f 

dt 

( f = f ) 

2 

( t) 

< ∞ (1)……1 

2 T T 

+ y de cuadrado integrable ( ∈ L − , ) 

aproximarse en forma única por una serie trigonométrica con la siguiente ecuación, 

donde ω 

2π 

T 

a 

∞ 

∞ 

f 

2 n= 

1 

n= 

1 

n 

0 

( t ) + an 

cos( 

ωnt ) + bn 

sin( 

ω t) 

0 = y n 0 

= ……2 

f puede 

ω = nω 

se denomina frecuencia fundamental y subarmónicos de f. 1 

Expresada en forma exponencial, será: 

f 

( t ) 

∞ 

= 

n= 

−∞ 

C 

n 

e 

iw t 

n 

C 

1 

= 

T 

T 

, donde: n 

( t ) 

2 

−T 

2 

f 

e 

−iw 

t 

n 

dt 

2 2



. , # # + " 

Sea una función periódica (fT(t)), de cuadrado integrable en [ , ] T 


137 

2 2 

T − , que pueda ser 

expresada en términos de la Serie de Fourier y una función no periódica (f(t)), de cuadrado 

integrable, tal que: 

f T = f en [ , ] T − 

( t ) ( t ) 

En términos de serie de Fourier, tenemos 

f 

f 

( t ) 

( t') 

∞ 

= 

n= 

−∞ 

∞ 

= 

n= 

−∞ 

C 

C 

Si T → ∞ ; entonces fT ( t ) → f ( t' 

) 

Por otro lado, en el intervalo [ , ] T 

2π 

Además, como ωn+ − ωn 

= = ∆ω 

T 

También, 

Entonces: 

∞ 

n= 

−∞ 

→ 

n 

n 

e 

e 

2 2 

T 

1 T 

T 

iωnt 

, donde: Cn 

= f ( t ) 

1 T 

2 

−T 

2 

iωnt 

' 

2 

, donde: Cn 

= f ( t ') 

T −T 

2 

f 

= 

1 

T 

e 

e 

T 

2 

−iω 

t 

n 

−iω 

t' 

2 2 

T − : ( t ) ( t ') 

1 , y si T ∞ 

∞ 

−∞ 

∞ 

n= 

−∞ −T 

2 

f 

n 

dt 

dt' 

e 

→ ; entonces ∆ ω = dω 

−iω 

nt ' iωnt 

dt ' e 

5 Mario Caicedo y Milagrosa Aldana. Aplicaciones Geofísicas de los Filtros Digitales. Versión 1.1, Sartenejas. 

Julio de 2002



f 

= 

∆ω 

2π 

( t ) ( t') 

lim 

T →∞ 

f 

∞ 

n= 

−∞ −T 

2 

1 

= 

2π 

∞ 

( t ) ( t') 

( t') 

−∞ 

T 

2 

dt'. 


138 

f 

∞ 

−∞ 

e 

−iω 

t' 

n 

dω. 

f 

i 

dt' 

. e 

e 

−iωt 

' 

ω t 

e 

n 

iωt 

= 

1 

2π 

∞ 

−∞ 

dω. 

e 

iωt 

1 

2π 

∞ 

−∞ 

f 

e 

−iωt 

' 

La expresión anterior permite definir la ecuación de la Transformada de Fourier: 

f 

f 

∞ 

1 −iωt 

' 

ω = f e dt' 

2π 

( t )( ) ( t ') 

. , # # % # + " 

Esta definida con la siguiente ecuación: 

f 

1 

= 

2π 

∞ 

−∞ 

iωt 

( t ) f ( ω ) e dω 

. , # # + " ( +" ! ! $ /02 * 

−∞ 

Sea una serie de tiempo con N muestras (X(kT), k = 0, 1, 2,…,N-1) tomadas con un 

intervalo de muestreo T, la Transformada de Fourier Discreta será una señal discreta de N 

muestras en el dominio de la frecuencia ( X (xf), k = 0, 1, …, N-1) registrada con un intervalo de 

muestreo 

1 

f = y viene dada por la siguiente ecuación: 

NT 

2πi 

− 

N 

Si se define W = e 

Entonces, 

X 

( xf ) 

N −1 

= 

N −1 

n= 

0 

e 

kn 

−2πi 

X ( nT ) 

N 

nk 

X ( xf ) = W x( 

k ) , con n= 0, 1, … N-1 

n= 

0 

dt'



. , # # + " % ( +" ! ! $ /02 * 

Se define con la siguiente ecuación, 

( nT ) 

−1 

n= 

0 


139 

X 

. + ' # 99 * 

= 

1 N 

N 

e 

np 

2πi 

N 

Comportamiento que representa la serie de Fourier de una función de cuadrado 

integrable, pero con saltos finitos en el límite de la discontinuidad y consiste en: 

figura, 

1. Se observa un salto finito de la serie aproximadamente en el límite de la discontinuidad 

de f(t). 

2. El aumento del número de términos de la serie no disminuye la intensidad o valor del 

salto, sino que lo corre al límite de la misma. 

Matemáticamente, sea la extensión periódica de la función (1) 

F 

( t ) 

cuya serie de Fourier viene dada por: 

S 

4 

π 

∞ 

= 

−1, 

+ 1, 

X 

d 

p 

NT 

− π < t < 0 

……3 

0 < t < π 

( t ) = 

sin[ 

( 2m 

+ 1) 

. t] 

m= 

0 

1 

2m 

+ 1 

……4 

Si se suman algunos términos de la serie se puede generar la gráfica en la siguiente 

Gráfico de la suma de los seis (6) primeros términos de la serie de la ecuación 4. Se observa que cerca de la 

dicontinuidad se presenta la anomalía asociada al fenómeno de Gibbs. (1)



. . 

Superposición que se produce entre las altas y bajas frecuencias de la Transformada de 

Fourier de una señal y que obedece a un muestreo incorrecto de la misma. 

. , !' $ .$! # % "!' : 

Sean dos (2) señales causales f y g de variable real y valores reales, la convolución (y) 

entre estas dos señales viene dada por la siguiente ecuación, 

y s 

( t ) ( f g) 

( t ) = f ( t− 

s ) g ( ) ds 


140 

∞ 

= * ; donde s es un desplazamiento en tiempo. 

. , !' $ .$! # !' : 

−∞ 

Sean dos funciones f y g de variable real y valores complejos, la correlación entre ellas 

será (h), definida por la ecuación: 

h u 

( ) ( ) ( ) ( ) du g f g f corr 

t 

= t+ 

u 

= , 

∞ 

−∞



;* , !' # !# # ) ! /052 

;* !# ## ( $ * 


141 

8 ; 

Velocidad asociada a la relación tiempo – offset en el proceso de reflexión, asociada a la 

corrección por normal moveout. Viene dada por la siguiente ecuación. 

V 

nmo 

= 

t 

2 

x 

x 

− t 

donde x es la distancia fuente – receptor, t0 tiempo de reflexión para incidencia normal 

(x=0) y tx tiempo de reflexión para todo x. 

;* !# # $) ! * 

Velocidad relacionada a las propiedades físicas de cualquier material. 

;* !# # $ %. ! * 

Velocidad de propagación en cualquier formación rocosa, en muchas ocasiones es 

asociada con la litología. 

;* !# # * 

Velocidad interválica a través de una sección geológica integrada por más de una 

formación o litología. Su magnitud viene dada por el promedio cuadrático pesado de las 

velocidades interválicas: 

V 

i= 

1 

i 

2 

0 

n 

rms = 

2 

Vi 

ti 

i= 

1 

n 

t




142 

8 

* # $ # $ ) ! # * 

El reordenamiento de los datos en el dominio de common midpoint, pero incorrectamente 

llamados common depth point, es necesario debido a que los programas de procesamiento realizan 

los análisis y cálculos de distintos parámetros en este dominio. 

Los datos de campo vienen referenciados según las coordenadas de fuentes y receptores 

y estas coordenadas deben ser transformadas a coordenadas de procesamiento. Según la 

geometría en la siguiente figura, se muestra el camino recorrido por un rayo reflejado en algún 

punto de un reflector plano en profundidad. 

Fig C.1. Representación de los trayectos recorridos por los rayos generados de la fuente, reflejados en el reflector en 

profundidad y recibidos por los receptores. Se muestra la relación de las coordenadas referenciadas al punto 0, donde 

s es la coordenada de la fuente; y es la coordenada del punto medio; g es la coordenada del receptor y h la distancia 

fuente – receptor. [18] 

Entonces, cada traza grabada en un receptor será asignada al punto medio entre la 

distancia fuente – receptor. Y todas las trazas pertenecientes a la misma ubicación conformarán 

un cmp gather. [18]



Como en los levantamientos sísmicos, la ubicación de fuentes y receptores varían 

espacialmente, el número de pares fuente – receptor que generen trazas que se ubiquen en un 

mismo punto medio definirá la cobertura de ese punto (fold). 

En el caso en que los reflectores sean horizontales y no se presente variación horizontal 

de velocidades, el common mid point (cmp) es equivalente al common depth point (cdp). Por ello, 

como una de las suposiciones en la teoría del modelo convolucional es que la propagación de los 

frentes de ondas registrados en los levantamientos símicos ocurre en capas planas, homogéneas 

e isótropas; los cálculos de los distintos parámetros, realizados por las funciones aplicadas 

durante el procesamiento, deben ejecutarse con los datos en el dominio de cdp. 

cmp 

cdp 

Fig. C.2. Representación del cmp y cdp. [18] 


143




144 

8 

Breve descripción geológica de las Arenas descritas en el levantamiento vecino, por el 

Este, del proyecto Arecuna Oeste 04G – 3D. Información tomada de [8]. 

* < 7=: cuerpo arenoso tipo canal, con espesores de hasta 30’ en algunas 

zonas. Hacia el Este y Oeste el cuerpo arenoso no se reporta la presencia de arenas con espesores 

significativos, lo que sugiere la presencia de una llanura de inundación en las direcciones 

señaladas. 

*>



* : intervalo representado fundamentalmente por un cuerpo arenoso en el 

sector central del área. Se asocia a un ambiente fluvial de ríos entrelazados. 

* B: cuerpo arenoso del que derivan dos (2) canales, permitiendo interpretar 

una llanura de inundación de ríos entrelazados. 

* >: cuerpo representado por cinco (5) canales de arena, ubicados en la parte 

central y Oeste del levantamiento. Hacia el Este en lo pozos no se muestran cuerpos arenosos, 

indicando una llanura de inundación en un sistema de ríos entrelazados. 

* B7>: al Sur se presentan dos (2) canales. Los cuales convergen hacia la 

parte central y se encuentran flanqueados por áreas pobres en arenas, lo que sugiere una llanura 

de inundación en un sistema de ríos entrelazados. 

* B

secuencia de procesamiento convencional hasta migración pre

Create successful ePaper yourself

Delete template?

Save as template?