Inversión del campo de ondas sísmicas PP y tiempos de llegada de ondas convertidas PS para la estimación de un modelo de velocidades P y S en el subsuelo: Un enfoque Level Set
No Thumbnail Available
Date
2024-03-06
Authors
Journal Title
Journal ISSN
Volume Title
Publisher
Universidad Industrial de Santander
Abstract
La velocidad de onda cortante (Vs) es una propiedad fundamental de los medios elásticos cuya estimación a partir de ondas convertidas PS es desafiante y requiere modelar la interfaz donde ocurre la conversión de P a S. Este artículo presenta una tomografía PS donde los puntos de conversión/reflexión de ondas sísmicas corresponden al reflector geológico modelado con la función Level set en su nivel cero (ϕ(x,z)=0). El método propuesto pretende una inversión estable de Vs en un entorno de adquisición sísmica utilizando receptores multicomponente. Se utilizan modelos sintéticos que simulan verdaderos Vs, Vp y la ubicación del reflector geológico.
Los tiempos de viaje de las ondas PS convertidas y las ondas PP reflejadas, tanto para datos observados como modelados (problema directo), se calculan utilizando la metodología propuesta por Rawlinson y Sambridge. Este método utiliza los tiempos de llegada de las ondas P desde la fuente sísmica hasta cada punto en el reflector como tiempos iniciales, que originan los tiempos Tps y Tpp. Estos tiempos se determinan como soluciones a la ecuación eikonal mediante el método Fast Marching.
La metodología plantea un funcional definido a partir de los residuales de los tiempos de viaje (modelados-observados), en función de las variables ϕ, Vs y Vp. Se propone minimizar el funcional utilizando el método de variaciones sobre el funcional aumentado, basado en el Lagrangiano asociado. La inversión toma los multiplicadores de Lagrange como variables adjuntas que actualizan las variables Vs, Vp, y ϕ(x,z), haciendo que los residuales de los tiempos se minimizen en cada iteración. El rendimiento del algoritmo se evalúa para modelos con geometrías de reflectores sinclinales, sinusoidales y monoclinas. La tomografía propuesta estima Vs y las posiciones de los reflectores, lo que puede contribuir a corregir las estáticas y mejorar la caracterización litológica de la superficie cercana.
Description
Keywords
Tomografía Sísmica, Ecuación Eikonal, Level Set, Onda Convertida PS, Estado Adjunto, Fast Marching, Función Green Empírica, Teoría Inversa, Modelos de Velocidad de Onda Cortante Vs