Numerical considerations on the modeling of source and boundary conditions for the frequency domain visco-acoustic wave equation solution

Abstract

ABSTRACT: Seismic modeling is an important step in the process used for imaging Earth subsurface. Current applications require accurate models associated with solutions of the wave propagation equation in real media. Unfortunately, it is common not to find in the technical literature deep discussions on the impact of specific details associated with the physical modeling of some crucial ingredients of the process, such as seismic source term and boundary conditions. In this paper, we discuss some issues related to the modeling of wave propagation in visco-acoustic media using finite differences. We focus our attention on two major elements of the modeling problem that are associated to the source term and the boundary conditions. We show that the source term can be modeled using a scale parameter that controls the spread of energy and shows that this parameter is a function of frequency and position of the source. As to boundary conditions, we show that Perfectly Matched Layer (PML) parameters are also frequency dependent. For both cases, seismic source scale parameter and PML model parameters we provide values and functions that optimize the performance of the approach for problems where visco-acoustic wave propagation is required. Frequency domain Full Waveform Inversion (FWI), or Reverse Time Migration (RTM) processes that depend fundamentally on the appropriate modeling of the wave-field are potential fields of application of these results.

RESUMEN: El modelamiento sísmico es un ingrediente importante en el proceso de construcción de imágenes del subsuelo. Aplicaciones actuales requieren soluciones precisas asociadas con la solución de la propagación de ondas en medios realistas. Desafortunadamente noes común encontrar en la literatura técnica, discusiones y análisis suficientemente profundas asociadas al impacto detalles específicos asociados con el modelamiento físico de ingredientes claves para este proceso como lo son el término de la fuente y las condiciones de frontera. En este trabajo se discuten algunos de esos detalles relacionados con el modelamiento de la propagación de ondas en un medio visco-acústico usando diferencias finitas. Nuestra atención se enfoca en dos aspectos importantes del modelamiento asociados con el término de la fuente y las condiciones de frontera. Se muestra que el término de la fuente puede ser modelado a través del uso de parámetros de escala que controlan la distribución de energía y se muestra que este parámetro es una función de la frecuencia y de la posición de la fuente. Para las condiciones de frontera se muestra que los parámetros del modelo de fronteras de tipo PML (Perfectly Matched Layer) también son dependientes de la frecuencia. Tanto para el parámetro de escala de la fuente como los parámetros que modelan la PML se provee valores y funciones que optimizan el rendimiento de dichas aproximaciones para problemas donde se requiere el modelamiento de la propagación de ondas en medios visco-acústicos. Inversión de onda completa (FWI) o Migración Reversa en Tiempo (RTM) son escenarios que dependen fuertemente del modelamiento apropiado del campo de onda y son campos potenciales de aplicación para los resultados que se presentan.