POSTSTACK

11. STACKING
1. Introducción:
- Estáticas
- Estáticas Residuales en combinación con análisis de velocidad
- Eliminación de ruido
 Sección aproximada zero offset

2. Estáticas
I) Correcciones estáticas (static corrections): Son correcciones en tiempo debido a la topografía (existencia
de desniveles entre fuente de energía y receptores) y a las variaciones de velocidad y espesor de la capa
meteorizada (de baja velocidad o weathering). Se establece un plano horizontal de referencia (PR o
datum) por debajo de ésta, que es el nuevo cero de los tiempos(figura de más abajo) y si el PR varía de
cota se lo llama flotante. Al tiempo total se le restan los tiempos que tarda el frente de onda en recorrer
la capa meteorizada hacia abajo y hacia arriba. Se hace una corrección por emergencia y otra por
incidencia, única para cada estaca o estación, de modo tal que se aplica la estática calculada a cada una
de las trazas de un registro individual. Otras veces se emplea un plano en el aire, por encima de la
topografía, y matemáticamente se rellena ese espacio con una velocidad de reemplazo semejante a la
que se estimó para la primera capa de alta velocidad Vc (la que infrayace a la meteorizada: consolidada o
subweathering). La velocidad de la capa meteorizada Vm (del orden de 400 a 1200 m/s) suele variar
lateralmente, mientras que la capa infrayacente (de unos 1700 a 2000 m/s) generalmente varía menos.
Los problemas que ocasiona la capa somera se resumen en la figura inmediata. La caída de alta
frecuencia es inatajable. Su espesor se calcula por sísmica de refracción, con dromocronas horizontales
(suele aplicarse la técnica de sumación, expuesta en el Tema 12) y/o dromocronas verticales (en pozos,
up o downholes). Estos últimos son más costosos y por lo tanto generalmente se hacen unos pocos en
puntos estratégicos para el control fino de las dromocronas horizontales. Si los registros de sísmica de
reflexión son lo suficientemente detallados para niveles someros, pueden usarse para el cálculo de
dromocronas horizontales, evitándose así la necesidad de hacer sísmica de refracción, aunque esto
sucede raramente debido a su muy distinta profundidad de interés y consecuentemente su escasa
confiabilidad. En el mar o lagos sólo debe considerarse su variación de profundidad, más el efecto de
marea. No existe bajo el agua una capa de baja velocidad como en tierra firme. La velocidad en el agua
promedia los 1500 m/s y la capa de sedimentos del fondo, embebida en agua, tiene una velocidad algo
mayor.
II) Correcciones estáticas residuales: Recurso de ajuste fino de las correcciones estáticas, si éstas no
hubieran sido muy buenas por las carencias en los datos de las estáticas de campo. Es un proceso que se
hace por comparación de trazas contiguas, donde el algoritmo busca eventuales saltos de tiempo
constante a lo largo de toda una traza respecto a
sus vecinas y, de haberlo, los corrige. En su
aplicación debe cuidarse de mejorar la
consistencia en el colgado de los tiempos de las
trazas, pero sin atenuar los reales rasgos
estratigráficos o estructurales del área. Como con
toda cosmética, tampoco emocionarse y que
termine pareciendo una mascarita de carnaval. A
la izquierda un ejemplo de este proceso opcional.
3. Método para realizar el Stack
. Stack medio

. Stack ponderado

Stack diverso/ Min-max excluidos

Método para realizar el Stack

 11. FILTROS ESPECIALES
11a. Entendiendo las múltiples:
FLUJO DE PROCESAMIENTO DE DATOS

OBJETIVOS:
- Entender el efecto de la velocidad en el filtrado del ruido
- Analizar la transformación de los datos desde el dominio t-x al dominio f-k
- Aplicación del proceso fk filtering
- Conocer la secuencia básica de procesamiento
- Procesar un dato símico real
INTRODUCCIÓN
- Las reflexiones múltiples se producen cada vez que la energía del shot se refleja más de una vez en
cualquier capa
- Las capas con fuerte contraste de impedancia acústica pueden causar reflexiones múltiples
- Las capas de agua también generan reflexiones múltiples. Es decir, en marina, la interfaz aire-agua y
la interfaz de agua-fondo permite causar una secuencia de fuertes reverberaciones.
- Estas (múltiples) aparecen con regularidad sucesivamente en los tiempo posteriores

CARACTERÍSTICAS DE LAS REFLEXIONES MÚLTIPLES

 Las reflexiones múltiples coninúan y se repiten con el mismo intervalo de tiempo y disminuyendo
gradualmente en amplitud.
  Las reflexiones múltiples tiene velocidades más lentas que las reflexiones primarias y por lo tanto,
toman más tiempo para alcanzar largos offset y tendrá grandes valores de move-out que aparecen
con más curvatura.
PROBLEMAS CAUSADOS POR LAS REFLEXIONES MÚLTIPLES
 Las múltiples de corto período modificarán el espectro de frecuencias de los datos, agregando
información no deseada al detalle fino en cada reflexión primaria
 La Deconvolución es utilizado para eliminar las reflexiones múltiples de corto período
 Las múltiples de largo período generan copias de los eventos primarios en tiempos posteriores en el
registro sísmico
 La discriminación de velocidad se usa para eliminar múltiples de largo período de datos sísmicos
porque su trayectoria de tiempo de viaje es mucho más larga que las reflexiones primarias.

ATENUACIÓN DE MÚLTIPLES: CORRECCIÓN NMO

 A) Antes de la corrección NMO, ambas las reflexiones primarias y múltiples tienen curvas de move-
out hiperbólico..
 B) Después de la corrección NMO con la velocidad de las reflexiones primarias, las múltiples son
todas sub-corregidas mientras los eventos pirmarios se aplanan
 Las reflexiones múltiples siempre tienen más move-out, es decir, mayor curvatura que la reflexiones
primarias. Es este move-out diferencial entre las reflexiones primarias y los múltiplos que ayuda para
eliminar múltiplos
TÉCNICAS DE MÚLTIPLES: CORRECCIÓN NMO
 Cuando las trazas en el gathers son stacked, las
reflexiones primarias se sumarán de manera coherente.
Debido a que las reflexiones múltiples aún tienen move-
out, no se sumaran de forma coherente
 Por lo tanto, el éxito del stacking de CMP en atenuación
de múltiples depende del grado de diferencia entre las
velocidades de los primarias y los múltiples.
TÉCNICAS DE MÚLTIPLES: FILTRO FK
 El gather puede ser corregido mediante el uso de una
función de velocidad que se encuentre entre las
funciones de velocidad primaria y múltiple.
 Esto tiene el efecto de sobre-corregir las primarias y subcorregir las múltiples
 Los eventos que se inclinan hacia arriba y hacia la derecha, después de aplicar la corrección de NMO
son las reflexiones primarias y las que se inclinan hacia abajo son las reflexiones múltiples.
 La diferencia en dip se usa para separar reflexiones primarias de reflexiones múltiples. Se usa un
filtro FK para eliminar los dispositivos en el dominio FK, y estas son las reflexiones múltiples.
TÉCNICAS DE-MÚLTIPLES: TRANSFORMADA PARABÓLICA
 La transformada parabólica de-múltiple se puede usar
para atenuar múltiples de período relativamente corto
 Esto permite la discriminación basada en la curvatura en
lugar del dip. Transformada parabólica proporciona una
mejor separación de reflexiones primarias y múltiples.
  La transformada suma los datos a lo largo de parábolas constantes. Las reflexiones primarias y
múltiples están bien separadas en el dominio parabólico

1)NMO correction en gather

2) Transformada Parabólica del registro
3)Sustracción de múltiples y transformación Inversa

11.b Filtros Especiales FK (Filtración de Velocidad)

INTRODUCCIÓN:

 Propósito del filtro fk: Atenual la energía no deseada como ,

- Ruido coherente lineal (ground roll).
- Ondas guiadas.
- Energía dispersa lateralmente.

 Estructuras inclinadas en el dominio (t-x) puede ser separadas en el dominio (f-k) por sus dips.

 El filtro fk se aplica a las trazas ordenadas por offsets

- Shots gathers.
- CDP gathers
- Stacks data

 Trabaja como un filtro 2D: proceso multicanal

FILTROS ESPECIALES
Transformación de Fourier Temporal (tf) Aplicaciones:

f-kx trasnformación
Transformación de Fourier Espacial (xk)

Aplicaciones:

Transformación Radón (τ-px ) Transformada de radón lineal y

tranformada de radón parabólica

DOMINIOS DE TRANSFORMACIÓN

DUALIDAD: Temporal vs Espacial

 Transformada de Fourier vs Espacial
- Intervalo de muestreo, Δt
- Taza de muestreo (Frecuencia de muestreo) f=1/ Δt
 - Muestreo que va preservar todas las frecuencias hasta la frecuencia de Nyquist
 Tranformada de Fourier Espacial
- Intervalo de muestreo, Δx
- Taza de muestreo (frecuencia de muestreo k=1/ Δx
- Muestreo que va a preservar todas las frecuencias hasta la frecuencia de Nyquist
DEFINICIÓN DE PARÁMETROS

GRÁFICO EN EL PALNO F-K

Se usará la siguiente convención :
- Dominio t-x=Eventos con downdip hacia la derecha se les asigna dip positivo, mientras que los
eventos con updip hacia la derecha se asignan dip negativos
- Dominio f-k= Los dips positivos se asignan al cuadrante de la derecha, que corresponde a números
de onda positivos, mientras que los dips negativos se asignan al cuadrante izquierdo, lo que
corresponde a números de onda negativos
- Zero dip es equivalente a cero números de onda
- La magnitud de los spickes corresponde al pico de amplitud de las sinusoides que componen las
trazas en la sección
- El plano f-k representa el espectro amplitud 2-d de la sección en el dominio t-x
- Este proceso se describe matemáticamente con la transformada de Fourier 2-d
RELACIONES PRACTICAS
CRITERIOS SOBRE EL MUESTREO ESPACIAL

- Desde un punto de vista práctico, las mediciones deben llevarse a cabo de tal manera que los
eventos en trazas separadas Δx puedan correlacionarse como si fueran del mismo horizonte o Punto
de reflexión en el subsuelo
- Para un componente de frecuencia dada, el tiempo de demora (Δt) entre dos mediciones
consecutivas debe ser menor que la mitad del período (T/2) de esa componente de frecuencia para
permitir una correlación de dos reflexiones como procedentes del mismo horizonte.

Delay tiene máximo

ALIASING ESPACIAL: INFLUENCIA DE LA FRECUENCIA

- Cada sección tiene una frecuencia única asociada a un n° de onda.
- A estas secciones zero-offset se las pueden considerar como representaciones de ondas planas, que
se propagan en un ángulo único desde la vertical y lleva una señal monocromática.
- El frente de onda se define como la línea de fase constante, mientras que la dirección de
propagación es perpendicular al frente de onda.
- Dado que un campo de ondas sísmicas (t-x) es una superposición de muchos dips(Δt/Δx) y
frecuencias (f), es equivalente a la síntesis de muchas componentes de onda planas.
- En este sentido, el significado físico de la transformada de Fourier 2D es importante, ya que es un
análisis (descomposición) de un campo de ondas en sus componentes de ondas planas.
- Para una dip dado, todos los componentes de frecuencia se grafican en el
- plano f-k a lo largo de una línea recta que pasa por el origen
- Cuanto más alto es el dip, (en el plano f-k) la línea recta radial se acerca al eje del wavenumber (k).
- Para los dips de 9, 12 y 15 ms/trace, las frecuencias con alias espaciales se ubican a lo largo del
segmentos lineales que rodean los cuadrantes opuestos del plano f-k.
- Cuanto más empinado es el dip, más baja es la frecuencia a la que se produce el aliasing espacial.

INFLUENCIA DEL DIP Y ΔX SOBRE EL ALIASING

 11c. FILTROS

EVENTOS NO LINEALES

ALIASING ESPACIAL
PROCESO DE FILTRADO DE VELOCIDADES
- Separación de la señal con el ruido en base a los
diferentes DIPS o VELOCIDADES APARENTES
- Definir las líneas de corte de velocidades
positivas y negativas. Todo lo que esté entre
estas líneas queda, y lo demas se eliminan
- Multiplicar los datos en la zona de paso por 1 y
en la zona de rechazo por 0
- Aumente linealmente el multiplicador de 0 o 1 en
zonas de rampas estrechas justo fuera de las 2
líneas de corte

- Tranformar el dominio F-K en x-t nuevamente

CUESTIONES PRÁCTICAS

La Transformada de Fourier produce aliasing (ruido envolvente)

 Extender los datos rellenando con trazas ceros.

 El gather de entrada se incrementara en por 4 trazas.

 Esto aumenta el costo, pero elimina el aliasing.

El ancho del abanico o polígono no debe ser demasiado estrecho.

 Si el ancho de banda de la zona de deyección es estrecha la respuesta (t, x) del dip filter tendrá muchos valores distintos de
cero.

 El espectro de amplitud del fk filter no debe tener limites abruptos.

La zona de reyeccion puede cubrir componentes de frecuencias aliadas

 Esto influye en el rendimiento del fk filte

VER EJERCICIOS

11.C Filtros especiales -TAU-PI (1)

LA TRANSFORMADA SLANT-STACK

- La transformada de Fourier 2-D es una forma de descomponer un campo de ondas en sus componentes de
ondas planas, cada uno con una frecuencia única y viajando en un ángulo único con respecto a la dirección
vertical.
- Vamos a presentar otra forma de descomponer el campo de ondas en sus componentes de ondas de planas,
usando el dominio de parámetro rayo p.
- Para un shot gather común, la descomposición en sus componentes de onda , se puede lograr mediante la
aplicacón lineal moveout (diferentes p) y sumar las amplitudes para todos los offsets.
- Este procedimiento se llama slant stacking (apilamiento inclinado)
- En el apilamiento inclinado se reemplaza el eje de offsets x por el parámetro de rayo p.
- El parámetro de rayo ρ es el inverso de la velocidad de fase horizontal
- El grupo de trazas con un rango de valores ρ es llamado gather slant stack

ASPECTOS FÍSICOS PARA SLANT STACK:ONDAS PLANAS

Sen ϴ=v.Δt/Δx  Δt= Sen ϴ* Δx /v

 Δt= p* Δx
Donde Δt es el tiempo de retraso
Fig. 1 geometría del raypath-frente de ondas línea de receptores

- Si el parámetro pi se concoe, entonces el rayo se puede trazar en un modelo en capas horizontales, como una
función de la velocidad conocida. V(z)
- Pi=0 onda que viaja verticalmente
- En general los receptroes para todos los offsets registran ondas planas de muchos valores de Pi
- Para descomponer un shot gather en sus componentes de ondas planas, se deben sumar a lo largo de cada
trayectoria inclinada las amplitudes de las trazas en el shot gather, cada una con un tiempo de demora único
definido por la ec.  Δt= p* Δx

TRANSFORMADA SLANT STACK

se usan dos paso para sintetizar las ondas planas, sumando amplitudes a lo largo
de path inclinados para todos los offsets x.
1° La CORRECCIÓN MOVEOUT LINEAL (LMO) se aplica a los datos a través
de una transformada de coordenadas definida por
τ = t – px donde p es el parámetro del rayo, x es el offset, t es el tiempo de viaje doble, y τ es el
tiempo de intercepción en p = 0

2° los datos se SUMAN según el eje del offset x

Al repetir la corrección de LMO para un rango de valores de p y realizar la suma de la ecuación, se constituye un
completo slant stack gather.

GATHER SLANT STACK

Un SLANT-STACK gather se conoce como τ-pi Gather, que consiste de TODAS LAS COMPONENTES DE DIPS,
dentro de un rango específicado de valores de pi, en el datos de offset originales
GATHER SLANT STACK: TRANSFORMADA LINEAL TAU-PI

- La transformada tau-pi es un caso especial de la transformada donde los datos se descomponen en una serie de
líneas rectas que se grafican como puntos en el dominio TAU-P y los eventos hiperbólicos se grafican como
curvas elípticas.
- Este proceso es usado para referirse como slant stacking que produce TAU-P domain, donde los datos de
entrada son sumado a lo largo de una serie de lineas rectas
- El proceso TAU.P es común de aplicar antes de la deconvolucón predictiva para la supresión de múltiples.
- La Tau P, transform puede ser usada también para aislar y filtras ondas guiadas, refraciones y otros tipos
deinterferencia.
- El filtrado en el dominio TAU-P es generalmente más caro que en el dominio FK y puede generar mejores
resultados.
- La hipérbola en el dominio T-X se tranforma en una serie de valores tangenciales tau-pi
- El resultado es una elipse en el dominio tau-pi
-

-
P1 Y P2 Son
eventos primarios
W es el fondo del mar con sus múltiples

ATENUACIÓN DE MÚLTIPLES OBJETIVOS

- Identificar diferentes tipos de múltiples, explicar sus características y las técnicas para atenuarlas
- Aplicar el método de las trans tau-p PARA REMOVER LAS MÚLTIPLES
- Tareas en los diferentes dominios para remover las multiples

TÉCNICAS DE ATENUACIÓN DE MÚLTIPLES

1. STACK
2. STACK ponderado (weighted stack) / Stack de eliminación de ruido de trazas internas
3. Deconvolución predictiva
 4. Predicción de la ecuación de onda
5. τ-p lineal con deconvolucón
6. τ-p parabólico (transformada de RADON
7. Discriminación de Amplitud y Frecuencia
8. Modelado relacionado con la superficie y sustracción
TIPOS DE MULTIPLES

CARACTERÍSTICA DE LAS MÚLTIPLES

 Considerar discriminación por velocidad

 Aplicar corrección NMO para remover múltiples
 Naturaleza repetitiva
 Con características estadísticas
TRAYECTORIA DE LOS RAYOS

USO DE LAS TRANSFORMADA DE RADON

- Filtro de Velocidad
- Supresión de múltiples
- Interpolación de trazas
- Análisis de Ondas Guiadas
TRANSFORMADA LINEAL τ-P
- Cada uno de los stack a lo largo de la pendiente pi, nos da una Amplitud de la posición
correspoindiente en el domino τ-p
- Luego repetir los mismos valores de pi, en el siguiente tau

FILTRO TAU-P
1. Registros shots transformados al domino tau-pi-
2. Aplicamos deconvolución predictiva en el
dominio tau-pi.
3. Inversa de la transformada lineal

TAU-P PARABOLIC TRANSFORM

  Es una versión de alta resolución de la transformada de Radon convencional
 Permite preservar las amplitudes de la primaria mejor que la transformada convencional y también
nos da una más completa remoción de las múltiples.
 Resiste el aliasing espacial de las múltiples que ocurren con algunas geometrías de adquisición
 Reduce la necesidad de la interpolación entre trazas antes de remover las múltiples.
RESUMEN:
 La señal múltiple se registró, debido a la energía atrapada en la capa de agua o en las capas de rocas.
 Los eventos de reflexiones múltiples obscurecen o enmascaran a los eventos primarios desde los
reflectores más profundos.
 Las múltiples tienen mayor normal moveout que la reflexiones primarias porque una gran parte de su
trayectoria esta cerca de la superficie con baja velocidad de la roca.
 Cuando el dato es stacked, las reflexiones múltiples tienen dips que son múltiplo del dip de las
reflexiones primarias
 El cambio de periodicidad de la múltiple con el tiempo y el offset x, dificulta la habilidad del proceso
de deconvolución para removerla.

12. MIGRACIÓN
1. Introducción

 La migración es el proceso de reconstrucción de una sección sísmica para que los eventos de
reflexión se vuelvan a colocar bajo su ubicación de superficie correcta y en un tiempo de reflexión
vertical corregido.
Objetivos :
- Colapso de la difracción
- Corregir el dip y la posición de las capas buzantes
- Mejorar resolución (Focalización de la energía)
2. Distorsión geométrica
1. Puntos difractor, usando el principio de Huygens: todos los objetos finitos pueden ser considerados
como una suma de puntos dispersos.
2. Dipping layer
3. Valley

2. Diagrama de Frente de Onda.

. Uso del diagrama de frente de onda y de la curva de máxima convexidad. Para migrar un evento reflexión
específico.
. El evento es tangencial a la curva de máxima convexidad y la posición migrada del evento esta en la
intersección entre el frente de onda y el ápice de la curva
3. Capa Buzante
La pendiente en la sección de stack siempre es menos comparado
con el dip real de una interfase

Sinclinal

4. Migración en tiempo/Profundidad

MIGRACIÓN EN TIEMPO

- No hay cambios de velocidad laterales

- Diagrama frente de onda (wavefront-charting )
- Migración por difracción (migración por kirchoff)
- Diferencia finita. Ecuación de onda (migración de diferencia finita)
- Migración fk
- Cambio de fase (gazdag )
- Migración stolt
MIGRACIÓN EN PROFUNDIDAD
Cambios laterales de velocidad
MIGRACIÓN DE KIRCHHOFF (Migración por
difracción)
MIGRACIÓN FK
EL HORIZONTE KX ES INVARIABLE EN LA MIGRACIÓN
. La migración en el dominio Fk involucra el mapeo de las líneas de frecuencia
constante AB en el plano (w-kx) a los círculos AB´en el
plano kz-kx
. una hipérbola de difracción está representada por un
área triangular invertida en el domino de frecuencia. La
migración convierte el área triangular en una forma
circular

. El Análisis FK de la hipérbola de difracción, se BASA EN LA PRESENTACIÓN

DE LA HIPERBOLA como una SERIE DE DIPPING DISCRETOS
. La asíntota se mapea a lo largo de la línea radial que representa EL LIMITE ENTRE LA PROPAGACIÓN
Y LA REGIÓN EVANESCENTES

5. PRE/POST STACK MIGRATION

Con variación lateral de velocidad Capas planas
PRESTACK migración en profundidad POST STACK migración en tiempo

BARATO  POST TIME

STACK MIGRACIÓN
CARO  PRE DEPTH
La migración POST STACK es aceptable cuando la sección stackeada es equivalente a una sección zero
offset
- PRESTACK FULL TIME MIGRATION
- PRESTACK DEPTH MIGRATION
Migración 3D
POST-STACK
DEPTH MIGRATION

PRE STACK
DEPTH MIGRATION
 12.B ENTENDIENDO LA MIGRACIÓN
 Cuando el reflector es horizontal, el punto de reflexión se encuentra en el punto medio entre el offset
fuente-receptor

T=2Z/V

El tiempo de viaje, T (two-way), para el evento reflejado es el doble de la distancia recorrida (dos
veces la profundidad, Z) dividido por la velocidad de la capa.

 Pero si el reflector se inclina, entonces el punto de reflexión no está directamente debajo del punto
medio entre el offset fuente-receptor, sino que se encuentra pendiente arriba del offset fuente-
receptor

T=2Z cos ϴ/V

El tiempo de viaje, Tw, (two,way) del evento reflejado para el reflector inclinado es igual al doble de
la profundidad, Z, multiplicado por el coseno del ángulo de inclinación dividido por la velocidad de
la capa.
2. Punto de reflexión aparente:
 En el caso de un reflector inclinado, si la energía sísmica viaja 0.3 s hacia el verdadero punto de reflexión
y 0.3s hacia el verdadero punto de reflexión y 0.3s hasta el receptor-el total es de 0.6 s
 Cuando estos datos se muestran en la sección sísmica, la reflexión se trazará 0,6s por debajo del
punto medio fuente-receptor (punto de reflexión aparente, A) pero no en el punto de reflexión real
( es decir, pendiente arriba del offset fuente-receptor)
 Cuando se aplican las correcciones NMO y DMO para ubicar todas las fuentes y receptores en la
posición del punto medio, se conoce el tiempo de reflexión aparente, pero no sus verdaderos puntos
de reflexión.
 Entonces, las reflexiones de los eventos dipping aparecen ubicadas con error en la sección sísmica no
migrada. Deben moverse”pediente arriba” a lo largo de una curva hiperbólica para colocarlos en su
punto de reflexión correcto y la forma de esta hipérbola depende del campo de velocidad.

3. Desplazamiento Espacial y de tiempo

El proceso de migración se usa para reposicionar los datos desde su punto de reflexión aparente, hasta su
verdaero punto de reflexión, tanto en el espacio y el tiempo. El desplazamiento espacial usado para
reposicionar los datos es calculado usando:
 ΔX se puede calcular picando los dips en la sección sísmica no migrada y
utilizando las velocidades de stacking

 El tiempo también se calcula picando los dips en la sección sísmica no migrada y

utilizando las velocidades de stacking derivadas de CMP/CDP gather con el
análisis de velocidad.

 A partir de estas ecuaciones, tenga en cuenta que cualquier evento se

moverá hacia arriba y hacia abajo en el tiempo.

 Si tomamos las dos ecuaciones de migración, podemos ver que el

término de dip no aparece en la segunda ecuación.

 Entonces podemos derivar todas las posiciones posibles para el evento

migrado, conociendo el tiempo del evento en la sección no inmigrada y
la velocidad.

4. Velocidad de migración
 La velocidad de migración se define como la velocidad que optimiza el reposicionamiento de la
energía reflejada a las ubicaciones correctas.
 Las velocidad de stacking se usan normalmente como base para determinar un “campo de velocidad”
para la migración.
 Cuando se lleva a cabo un análisis de velocidad de un survey 3D, el volumen 3D de las velocidades
determinadas se conoce como campo de velocidad.
Tenga en cuenta que:
 Si los datos se migran con una velocidad demasiado baja, los datos estarán sub-migrados
 Si los datos se migran con una velocidad que es demasiado alta, los datos estarán sobre-migrado
 Sección no migrada  pobre resolución

5. Estructuras geológicas
ANTICLINALES
 LOS ANTICLINALES se muestran más anchos en las secciones sísmicas no migrads más allá de su verdadera
psoición y se cruzan por los puntos de reflexión. Cuanto mayor sean los dips, más amplia aparecerá la
estructura en la sección no inmigrada.
 En la sección sísmica migrada, aparece como un
anticlinal comprimido. Si las velocidades utilizadas para
migrar son incorrectas, la estructura final migrada
puede ser más estrecha o más ancha que la verdadera
estructura.
 Además, velocidades de migración incorrectas puede
conducri a incertidumbres laterales y verticales cuando
se utilizan los datos sísmicos para hacer well-to seísmic
tie durante la interpretación sísmica

 LOS SINCLINALES, Tienen una forma cóncava. En la sección sísmica no migrada, un sinclinal aparece más
estrecho de lo que realmente es y también todo el sinclinal se mostrará expandido en las sección sísmica
migrada.

Las velocidades incorrectas conducirán a una forma incorrecta:

- Si las velocidades son demasiadas ALTAS, entonces el sinclinal será más amplio de lo que debería ser
despues de la migración
- Si las velocidades que son demasiado BAJAS sobre-migran la estructura, dejándola estrecha.

 BURIED-FOCUS ANTICLINAL(foco anticlinal enterrado) es un sinclinal apretado que en una sección sísmica
tiene la forma sísmica de un anticlinal y esto puede confundirse con un potencial trampa de petroleo y o de
gas en una sección sísmica.
Tenga en cuenta que las trayectorias de rayos en el centro del
sinclinal formarán un “bow-tie” porque las reflexiones desde
cualquier lado de la estructura se cruzan. Esto resulta en el foco
anticlinal enterrado y será corregido por el proceso de
migración.

 FALLAS afectas la respuesta sísmica de dos maneras

- 1° Lugar, debido al dip del plano de falla, las reflexiones desde el plano de falla no están ubicadas en el punto
medio del offset fuente-receptor, sino en un punto pendiente arriba ene l plano de falla. Por lo tanto, la
respuesta sísmica para un plano de falla inclinado es lo mismo que para una interfaz inclinada.
- 2° Lugar, la falla puede causar que la formación termine abruptamente. La energía sísmica directa en el
punto de determinación de la formación, se refleja en todas las direcciones.

Como la distancia desde el punto aumenta, el tiempo de viaje aumenta y la reflexión resultante produce una
curva de difración la curva de difracción aparece como un anticlinal (fig. 12.13)

- El ápice de la curva de difracción se encuentra en el punto real de reflexión. Este se llama punto de
difracción. Todas las reflexiones registradas en la curva de difracción se pueden atribuir a este punto.
- Uno de los objetivos del proceso de migración sísmica es colapsar todos los eventos reflejados en la curva de
difracción, llevarlos de nuevo al verdadero punto de reflexión y restaurar la brusquedad de la reflexión de
falla y la falla se muestra con más apariencia en la sección migrada.
- La migración también reposiciona un reflector inclinado, ya sea una interfaz estratigráfica o un plano falla a
su ubicación espacial correcta.
- En la figura 12.14 se puede ver las curvas de difracción en la sección no migrada. A la derecha está la sección
migrada y observe que el proceso de migración ha colapsado
las vurvas de difracciones y mueve los eventos inclinados
lateralmente en la dirección de pendiente arriba y tiempos
hacia arriba o antes en tiempo. Es decir, la imagen migrada
tiene una mayor inclinación que imagen no migrada. La
cantidad de movimiento depende de la inclinación (dip), la
velocidad y el tiempo del evento. Si el campo de velocidad es
incorrecto, el plano de falla no se moverá al lugar correcto
- Si los datos se migran con una velocidad más rápida que la velocidad correcta, veremos difracciones en la
sección migrada convertida en formas sinclinales llamadas sonrisa. Si las velocidades son demasiada lentas,
la sub-migración dejará un ceño residual en los datos. Migrando con la correcta velocidad se colapsa todas
las difracciones en un punto
- Tenga en cuenta también que si hay un único punto de evento reflejado en la sección no migrada que no
corresponde a un evento real, por ejemplo, un pico, este migrará a una sonrisa parabólica en la sección
migrada.

6. Migración en tiempo
La migración en tiempo es el reposicionamiento de los datos sísmicos reflejados desde sus puntos de
reflexión aparente, desde la mitad del camino entre la fuente y el receptor hasta sus verdaderos puntos de
reflexión en el espacio y tiempo
 La migración en tiempo utiliza tres principios:
- Método de migración de Kirchoff
- Método migración de diferencias finitas
- Método de migración de transformada de Fourier.
Pero solo discutiremos dos de los principios de migración en tiempo
MIGRACIÓN DE KIRCHOFF
Todos los métodos de migración derivan de la solución a la ecuación de onda. La ecuación de onda es
expresada a continuación.
Donde:
u: es el campo de ondas sísmicas,
(x,y,z): son las tres coordenadas de un sistema tridimensional.
t: es el tiempo;
V: es la velocidad de propagación o la relación de
velocidad a la que la onda sísmica se propaga.
 La migración de Kirchhoff usa la integración
para resolver la ecuación de onda.
 El método de migración de Kirchhoff utiliza la
geometría y el principio de Huggen para
colapsar la difracción y reubicar los datos
registrados. Considera el aplace de la curva de
difracción como la ubicación del verdadero
punto de reflexión (figura 12.16)
  La migración de Kirchhoff colapsa la difracción al sumar las amplitudes a lo largo de la curva de
difracción y luego colocando la suma, en el punto de refleción verdadero (el apice de la curva de
difracción). Esta pondrá la energía que pertenece al difractor en su posición correcta en la traza
migrada. Este proceso se realiza para cada traza migrada, y la entrada consiste en la seción no
inmigrada.
 Hay dos parametros que son importantes para el método de migración de Kirchhoff
- El primer parametro esta relacionado con la curva de difracción en sí misma, se llama parámetro
ancho de apertura. El éxito de la migración de Kirchhoff requiere una buena definición de la
difracción. Es importante tener en cuenta que la forma de la curva de difracción varía con el tiempo y
la distancia.
- Se debe incluir un número suficiente de trazas en la ventan a de procesamiento para definir
completamente la curva de difracción. Si la curva de difracción está parcialmente registrada o si el
punto difractor está tan cerca de la superficie, es decir solo se registra una curva pequeña, entonces la
curva de difracción puede colapsar en el punto incorrecto y los datos pueden ser incorrectamente
migrados.
- El modelo de velocidad es el segundo parametro critico cuando se utiliza una migración de Kirchhoff
- La migración de Kirchhoff es extremadamente sensible a pequeños cambios en el modelo de
velocidad. Por ejemplo, usando una velocidad 5% más rápida que la velocidad verdadera, dará como
resultado la sobre-migración de los datos y también una velocidad 5% más lenta que la verdadera
velocidad, resultará en la sub-migración de los datos.
 Tenga en cuenta también que el método de migración de Kirchhoff no puede manejar datos sísmicos
con baja relación señal-ruido y datos con variación de velocidad lateral.
 La migración de Kirchhoff es la más comúnmente utilizada para el procesamiento de profundidad.
 Es particularmente atractiva debido a su velocidad y capacidad orientada a objetivos, lo que permite
una construcción del modelo de velocidad eficiente y de actualización de modelos.
 Además, la migración de Kirchhoff ofrece adaptabilidad para geometrías de adquisición irregulares,
flexibilidad en el manejo de anisotropías, modelos de velocidad de ondas convertidas y topografía de
superficie.

MIGRACIÓN POR DIFERENCIAS-FINITAS

 Este segundo método de migración utiliza la diferenciación para obtener una solución a la ecuación
de onda y generalmente se conoce como métodos de diferencia finita.
 Diferencia finita también conocida como la migración de
continuación descendente que mapea la curva de difracción (figura
12.18) registrada en una profundidad a otra, como si esta se registrara en
la nueva profundidad (extrapolada) y predice el cambio en la forma de
onda generada, con el cambio de profundidad.
 En otras palabras, la migración de diferencia finita extrapola
datos de entrada (secciones CMP stack o pre stack)usando incrementos
finitos de profundidad y predice cómo se verían los datos en la nueva
profundidad.
 Migración como un proceso: Downward continuation

 La entrada al proceso de migración, es la

sección sísmica stack de zero offset de acuerdo a
como fue registrada en la superficie, conocido el
modelo de velocidad y la ecuación de onda que
 describe el proceso de propagación, el campo de ondas es reconstruido a un nivel de profundidad
más bajo (Δz debajo de la superficie).

 Este proceso es llamado downdward continuation que significa que el campo de ondas es
recalculado, para simular una situación donde los receptores están debajo de la superficie.

En este nivel de profundidad, se selecciona la

parte del campo de ondas que corresponde a
tiempo cero. De acuerdo a la definición del
proceso de migración esta parte del campo de
ondas estaría migrado.

 El proceso downward continuation, sigue

con la propagación del campo de ondas desde el
nivel de profundidad Δz hasta el nivel de
profundidad 2 Δz y asi siguiendo hasta obtener la
máxima profundidad requerida.
 Hay dos parámetros que son importantes para el métodos de migración con diferencias finitas:
 El parámetro más importante es el tamaño Δz de los intervalos en el proceso de continuación
descendente, porque si los intervalos son demasiado pequeños se producirá una sub-migración. En
este caso, la difracción está incompletamente colapsada y el evento sísmico no se reposiciona lo
suficiente up-dip. Además, el uso de un intervalo pequeño no mejora significativamente los
resultados de la migración, pero requiere más tiempo y costo de la computadora.
 El modelo de velocidad es el segundo parámetro importante. Usando una velocidad demasiado
rápida causa la sobre-migración, mientras que la velocidad que es demasiado lenta causa la sub-
migración.
 Tenga en cuenta que la migración de diferencias finitas no es sencible a las variaciones de velocidad
como Kirchhoff migration methods.
 La migración de diferencia finita puede acomodar variaciones menores de velocidad lateral.
 El método de migración de diferencias finitas puede manejar datos con pequeñas variaciones lateral
de velocidad y también puede manejar datos ocn baja relación señal/ruido.
 Tenga en cuenta también que la migración de diferencias finitas toma mucho tiempo de computación
para procesar.

7. Migración en Tiempo
 En la migración en profundidad, los datos se migran en el dominio de la profundidad y el dato de
entrada para el proceso son datos sísmicos no migrados en el dominio del tiempo. El resutlado es el
dato migrado en el dominio de la profundidad.
 Tenga en cuenta que la conversión de tiempo a profundidad requiere de la velocidad promedio. La
velocidad promedio puede ser derivada de la velocidad de stacking calculando primero la velocidad
interválica.
 La migración de profundidad es extremadamente sensible al modelo de velocidad, por lo que el
modelo detallado de velocidad del subsuelo es necesario. El modelo de velocidad incorrecto produce
inexactitud en los resultados de la migración. Sin embargo, si el modelo de velocidad es preciso,
entonces los resultados finales del proceso de migración en profundidad son una mejora sobre los
resultados de los métodos de migración en tiempo..
 Tenga en cuenta que la migración en profundidad puede manejar variaciones de velocidad tanto
vertical como lateral, así como variación de velocidad moderada
 Notese también que la migración en profundidad requiere un tiempo de procesamiento más
prolongado, que conducen a un mayor costo de procesamiento de datos.
  La migración en profundidad requiere de datos muy precisos de velocidad de intervalo. La figura
12.20 muestra un modelo de velocidad de intervalo utilizado para la migración de profundidad.
 Hay tres principios de migración en profundidad:
-migración de profundidad de diferencia finita (finite-difference depth migration)
-migración de profundidad teórica de rayos (ray-theoretical depth migration)
-migración de profundidad de imágenes por trazado de rayos (Image-ray tracing depth migration)
DIFERENCIA ENTRE TIME Y DEPTH MIGRATION
 La migración de profundidad tiene en cuenta la curvatura de rayos mientras que la migración de
tiempo no hace eso.
 Las migraciones de profundidades se pueden aplicar a los datos pre-stack y post-stack
 Las tecnicas de migración en profundidad requieren más tiempo de computación para ejecutarse, que
la migración en el tiempo. Tenga en cuenta que la conversión de la sección de tiempo a una sección
de profundidad mejora la imagen produciendo un espesor más realista entre los reflectores-alta
velocidad produce un espesor pequeño aparente en la sección de tiempo pero la conversión de
profundidad se reserva esto.
 Tenga en cuenta también que, en aras de la consistencia, la sección de profundidad final se convierte
de nuevo en tiempo ya que esa
es la escala habitual para una
sección sísmica.
 En las secciones sísmicas que
se muestran en la figura 12.21,
capas correspondientes
velocidades de intervalo más
bajas, aparecen más juntas en
la sección de profundidad que
en la sección de tiempo. Capas
correspondiente a velocidades
de intervalo más altas aparecen
más separadas en la sección de
profundidad que en la sección
de tiempo.

[Link]-stack depth migration

 Algunas estructuras son imposibles de stack correctamente, sin primero migrar los datos antes stack.
 Las migraciones de pre-stack time y pre-stack depth usan todos los datos generalmente en planos de
offset común, como entrada al proceso de migración.
 Tenga en cuenta que los gathers CDP generados a partir de los datos prestack se utilizan para el análisis
AVO con la preservación de la amplitud como una función de offset.
PSTM (PRE->STACK TIME MIGRATION)
 Si el campo de velocidades varía lateralmente ligeramente, la migración
en tiempo pre-apilamiento (PSTM) dará una buena imagen de las
estructuras geológicas en el subsuelo (figura 12.22).
PSDM (PRE-STACK DEPTH MIGRATION)
si la velocidad varía rápidamente lateralmente, se requiere la migración pre-
stack depth para lograr una clara imagen de las estructuras geológicas en el
subsuelo

9. Post-stack Migration
 Si las capas son más o menos planas con pequeñas variaciones en la velocidad, una migración post-
stack proporcionará una buena imagen de la estructura geológica (fig.
12.24)

 Tenga en cuenta que la migración prestack mejora mucho los

resultados, si la estructura dentro del área es altamente compleja.

 Al igual que la migración en profundidad, la migración pre-stack

requiere más tiempo de computadorea y por lo tanto cuestan más que la
migración post-stack

COMPARACIÓN POST-STACK MIGRATION – PRESTACK MIGRATION

 La figura 12.25 compara las secciones
sísmicas 3D, post-stack migradas (izquierda) y PSTM.
En el basamento fracturado con geología compleja, el
conjunto de datos procesados con PSTM proporciona
mucho mejor imagen.

10. Filtros finales

Los filtros se aplican al final de la secuencia de procesamiento de datos para realzar el stacking final o
sección migrada. Recuerde que el término filtrado generalmente significka iltrado de frecuencia. El filtro se
usa para reducir o atenuar los componentes de frecuencia donde el ruido domina sobre los datos sísmicos
BAND-PASS FILTER
 Los pasos de procesamiento de datos, como la deconvolución, han intentado mantener el ancho de
banda de los datos sísmicos tan aplios como sea posible. Al hacerlo, los procesos de deconvolución
pueden sobre-compensar el espectro en lo stimepo iniciales y sub-compesar para los tiempo finales
y por lo tanto puede haber enfatizado demasiado el ruido más allá del espectro de la señal
 Además la claridad de la reflexión de los datos sísmicos puede cambiar de traza a traza por un
número de razones tales como ruidos diferentes, amplitud reducida del stack debido al dip y
diferencias de traza a traza en la de-convolución. Por estos motivos, el stack final o la sección
migrada necesitan ser filtradas.
 Filtrando las componentes de frecuencias más dominados por el ruido, hara que los datos se vean
más limipios y mejores. El filtro final aplicado a los datos se llama bandpass filter.
  Un filtro de pasa banda atenúa las frecuencias por debajo de su frecuencia de corte baja y por encima
de su corte de alta frecuencia. La diferencia entre estas dos frecuencias de corte es la banda de paso.
TYPES OF BAND PASS FILTERS
 Al hacer una elección de ancho de banda del filtro para aplicar
a los datos, es importante tener en cuenta que las frecuencias
bajas son importantes para comprender el carácter de
reflexión y las altas frecuencias son importantes para mejorar
la resolución de los datos
 Observe el fondo más limpio y las amplitudes enfocadas de las
reflexiones sísmicas después de filtras, asi como los bordes de falla conservados.
 Tenga en cuenta que el filtrado ha eliminado los componente de frecua dominados por el ruido en
datos sísmicos sin filtro y hace que la sección filtrada ( a la derecha) se más limpia y mejor.
 HAY DOS TIPOS DE FILTROS pasa banda que se aplican comúnmente alos datos sísmicos. Ellos
son los filtros de fase mínima y de fase cero.
 Cuando un filtro se aplica para mejor el stack final o la sección símica migrada, el filtro de fase cero
es aplicado
 Cuando un filtro se aplica a los datos símicos en los paso de rpocesamiento intermedio, se aplica un
filtro de fase mínima, ya que los pasos de desconvolución pueden depender de que la wavelet sea de
fase mínima
 La wavelet de fase mínima no tiene amplitud en el tiempo cero (figura 12.27ª), mientras que la
wavelet de fase cero tiene amplitud en el tiempo cero (Figura 12.27b)
 La wavelet de fase cero es preferida para la interpretación debido a que el pico central fuerte en el
tiempo cero es fácil de relacionar con el reflector de interés.
 El tiempo de reflexión para una wavelet de fase cero se encuentra en el pico o valle central (-) por lo
que los tiempos de reflexión son mucho más fáciles y precisos para wavelets de fase cero. Esta ses
una de las razones porque las wavelets de fase cero generalmente se generan en el procesamiento de
datos sísmicos.
APPLICACIONES DE BAND PASS FILTERING
Hay tres aplicaciones comunes de filtrado pasa banda:
-Band pass filtering simple. Se aplica un solo filtro a todas las trazas y en todo tiempo.
-time-variant Filtering: este permite a los geofísicos especificar un filtro variable en el tiempo que es
luego aplicado a todas las trazas.
Time and spatially variant filterin: este aplica una serie de filtros variables en tiempo y que varían a lo
largo de la línea sísmica.
TIME -VARIANT FILTER
- El contenido de frencuencia de la señal sísmica disminuye a medida que aumenta la profundidad
porque las componentes de frecuencia más altas de la señal sísmica son absorbidas a medida que la
señal se propaga más profundamente
 Debido a esto, el ruido en el stack o en la sección migrada puede ser de banda ancha. Por esta razón
los geofísicos prueban el filtro óptimo varias veces en la sección que proporcione una buena relación
señal/ruido.
- Por lo tanto, los geofísicos necesitan aplicar filtros variantes en tiempo
 Tambien por el espectro de la señal y la señal-ruido cambian con el timepo y aveces con el espacio, y
es más probable que un filtro no funcione para toda la sección.
- Entonces, el filtro final probablemente sea variante en tiempo.
 13. POST PROCESO
OBJETIVOS:
 Eliminación de efectos por :
- FILTRADO
- STACKING
- MIGRACIÓN
 Ajustar las amplitudes para el ploteo
 Mejorar la coherencia de las reflexiones
OPERACIONES USADAS FRECUENTEMENTE
 Filtro de frecuencia dependiente en el tiempo
 Filtro de coherencia
 AGC
FILTROS DE FRECUENCIA DEPENDIENTE DEL TIEMPO
RAZONES DEL FILTRADO DESPUÉS DEL STACKING:
 - Debido al uso de diferentes filtros se genera ruido adicional
- El contenido de frecuencia de los datos cambia en los diferentes pasos del procesamiento (NMO
streching)
- Gradualmente el filtro de frecuencia variables con el tiempo no se usa antes del stacking
FILTROS DE COHERENCIA
Objetivos: supresión de ruidos, eliminación de dips extremos, mejorar la amplitud y continuidad de las
reflexiones
Filtros de coherencias: FK Filters, fx deconvolucion (Wiener-Levingson)
AJUSTES DE AMPLITUD
AMPLITUD DE LAS REFLEXIONES
-Dependiendo de los objetivos de la interpretación
-estruras claras
-preservar las amplitudes y continuidad de las reflexiones
ADAPTAR LAS AMPLITUDES PARA PRINTING
-Agc
ALMACENAMIENTO DE LOS DATOS
Almacenamiento Digital : Formtao SEGY
Archivo ASCII
Archivo LAS
PLOTEO DE LOS DATOS
DATOS
INFO HEADER (Función de velocidad, Fold, Pto de intersección con otros perfiles, localización de pozos
de logueo , topografia, numero de shot

TEMA 5
CORRELACIÓN DE DATOS SÍSMICOS Y DE POZO

INTRODUCCIÓN: RESUMEN/OBJETIVOS/RESULTADOS
VINCULACIÓN A TRAVÉS DE:
 WELL SHOOT
 SONICO INTEGRADO
 SISMOGRAMA SINTÉTICO
 VSP

ELEMENTO USADOS EN LA CORRELACIÓN

1) Herramientas de Logueo: dipmeter (caliper)
Sónico
Densidad
 Microresistividad
2) Curvas y sismogramas
1. Curva T-Z derivada de las check shots
2. Simograma sintético (de well logs)
3. Sismograma: derivado del VSP por un corridor stack
4. Sección sísmica, derivada del offset VSP

1. CURVA T-Z DERIVADA DE LOS CHECK SHOTS

2. SISMOGRAMA SINTÉTICO: DERIVADO DE LOS WELL LOGS

DATOS DERIVADOS DE LOS

WELLS LOGS:

 Chequear las anomalías de

lecturas
 Los logs deben ser
adecuadamente editados
- Washouts: derrumbes
- Cicle Skips: Saltos de ciclos
- Zonas Invadidas
DATOS DERIVADOS DE
LOS CHECK SHOTS
- Chequear las anomalías de
Velocidades interválicas.
- Verificar la calidad de los datos originales
- Verificar la calidad de los picks de primeros arribos
3. SISMOGRAMA: DERIVADO DEL VSP

DISEÑOS GENERALES:
o Checkshots espaciado de geófono aproximadamente 100m
 fuente próxima
al pozo

o Zero offset VSP: Pozo vertical  espaciado

de geófono 20m
 fuente
próxima al pozo

o Walkabove VSP, Pozo horizontal

o Offset VSP, Imaging:

Pozo vertical
 o Walkaway VSP:Avo,
Pozo vertical

o 3D VSP, Pozo vertical

VSP: CONCEPTOS BÁSICOS

o Calculo del tiempo de
transito TT
 o One-way vs two way time

o VSP meno (-) TT

o VSP mas (+) TT

PROCESAMIENTO BASICO VSP

 Secuencia de procesamiento
- Preparación de los datos
- Edición de los datos
- Técnicas de stacking: median, mean

- Normalización, Tar, Filter

 Separación de los campos de ondas

- Estimación de las ondas directas: a) sustracción del tiempo de tránsito para el
alineamiento vertical de todas las energías downgoing
b) aplicacíon del filtro mediana para realzar la energía downgoing en fase y suprimir
todas las energías fuera de fase
c) cambiar cada traza hacia su one way time original
- Sustracción de las ondas directas
- Realce de las ondas reflejadas
- FK filters: Filtro de velocidad
 Deconvolución
- Cálculo del operados sobre la onda directa
- Aplicación del operador sobre la onda directa
- Aplicación del operador sobre la onda reflejada
 Corridor stack
- Ejemplo