Obtención de imágenes RTM (Reverse Time Migration) en zonas estructuralmente complejas

Abstract

Reverse Time Migration es un método de migración en profundidad que implementa una solución de la ecuación de onda completa. En este trabajo se desarrolla una metodología para la aplicación de RTM pre-apilado a datos sintéticos y reales. Se presentan y analizan los componentes de RTM: Modelamiento de la onda acústica, condiciones de contorno y condición de imagen. Para el modelamiento de la onda se implementó un método de Diferencias Finitas (DF) de alto orden en el dominio espacio-tiempo, basado en un esquema de grilla alternada; en el proceso de extrapolación se aplicó un operador espacial adaptativo de longitud variable, en el cálculo de las derivadas espaciales para reducir el costo computacional en forma efectiva, con poca reducción de la exactitud de la solución numérica. Se estableció una condición de contorno ABC híbrido que combina dos de las soluciones más comunes al problema de las condiciones de contorno: el método predictivo y el método de atenuación. Las imágenes RTM se obtuvieron a partir de una Condición de Imagen por Correlación Cruzada con Separación de Campos. Se aplicó el algoritmo a datos sintéticos obtenidos del modelamiento de una estructura geológica compleja, teniendo control del modelo de velocidades. Sobre datos reales, se obtuvo una imagen de una estructura compleja que permitió conocer la ubicación real de los eventos de reflexión. La imagen sísmica obtenida, presenta mejores eventos sísmicos en comparación con los métodos convencionales de migración.

Abstract RTM is a depth migration method that implements a solution of the full wave equation. In this work, a methodology for the application of RTM to synthetic and real data is developed. We discuss the components of RTM: Acoustic wave modeling, boundary conditions and imaging condition. A Finite Difference Method (FDM) based on a staggered grid scheme was implemented for the wave modeling. A Hibryd ABC method that combines two of the most common solutions to the problem of boundary conditions: predictive boundary condition and absorbing boundary condition. RTM images were obtained using a Cross-correlation Imaging Condition with Wavefield Decomposition. The algorithm was applied to synthetic data obtained by modeling a complex geological structure, with a known velocity model. On real data, an image of the complex structure was obtained, which allowed us to know the real location of the reflection events. The obtained seismic image better resolved the seismic events compared with conventional methods of migration.