Phase variance as a seismic quality-control attribute

Este artículo introduce la varianza de fase como un atributo de control de calidad sísmico basado en estadística circular que cuantifica automáticamente la dispersión de fase local en datos terrestres pre-apilados, permitiendo una evaluación objetiva y dependiente de la frecuencia de la fiabilidad de la señal para optimizar flujos de trabajo como la inversión de onda completa y la AVO.

Lo esencial

Imagina que estás intentando escuchar una conversación clara en medio de una fiesta muy ruidosa, pero con un giro extraño: el ruido no solo es volumen, sino que distorsiona el tono de las voces de manera impredecible.

Esta es la historia que cuentan Akshika Rohatgi, Andrey Bakulin y Sergey Fomel en su artículo sobre la "Varianza de Fase".

Aquí tienes la explicación, traducida a un lenguaje sencillo y con analogías creativas:

1. El Problema: El "Efecto Espejo Roto"

Cuando los geofísicos exploran el subsuelo (para buscar petróleo, gas o entender la tierra), disparan ondas sonoras y escuchan los ecos. El problema es que la capa superior de la tierra (el "near-surface") es como un suelo lleno de baches, rocas y arena suelta.

• La analogía: Imagina que lanzas una pelota de tenis contra una pared lisa; rebota perfectamente. Ahora, imagina lanzarla contra una pared llena de agujeros, clavos y espejos rotos. La pelota rebota, pero en direcciones locas y con tiempos extraños.
• La realidad: En la tierra, estas irregularidades hacen que las ondas sísmicas lleguen a los sensores con un "desfase". No es solo que lleguen tarde (como un retraso en una llamada de Zoom); es que la forma de la onda se ha deformado.

2. El Viejo Método: "El Promedio Mágico"

Durante décadas, los geofísicos han intentado arreglar esto usando un método llamado "deconvolución consistente con la superficie".

• La analogía: Es como si todos los invitados a la fiesta tuvieran un micrófono defectuoso. El método antiguo asumía que todos los micrófonos de la misma mesa tenían el mismo defecto. Así que, calculaban el promedio del error de la mesa y lo corregían para todos.
• El fallo: Esto funciona bien para errores grandes y predecibles, pero falla estrepitosamente cuando el problema es local y caótico (como un micrófono que se rompe solo en un segundo específico). Además, nadie tenía una regla para medir cuánto estaba "locura" la fase. Solo miraban las ondas y decían: "Parece mejor" o "Parece peor", basándose en la intuición humana.

3. La Nueva Solución: "La Varianza de Fase" (El Nuevo Termómetro)

Los autores proponen una nueva forma de ver el problema. En lugar de tratar cada traza sísmica como una línea individual, la tratan como un grupo de personas (un conjunto) y preguntan: "¿Qué tan de acuerdo están todos en cuanto al momento exacto en que ocurre el sonido?".

Para entender esto, usen la analogía de las agujas de un reloj:

• Señal perfecta: Imagina 100 personas, cada una sosteniendo una aguja de reloj. Si todas apuntan exactamente a las 12:00, están coherentes. La "varianza" es cero.
• Ruido total: Si el ruido es muy fuerte, las agujas apuntan a direcciones totalmente aleatorias (una a las 3, otra a las 9, otra a las 11). Están descoherentes. La "varianza" es máxima (1.0).
• La Varianza de Fase: Es simplemente un número que mide qué tan dispersas están las agujas.
  • Si el número es bajo (cerca de 0): ¡Excelente! Las ondas están sincronizadas y la imagen será clara.
  • Si el número es alto (cerca de 1): ¡Peligro! Las ondas están desordenadas. No importa cuánto amplifiques el volumen (la amplitud), la información es basura.

4. La Gran Revelación: "El Truco de la Amplitud"

El descubrimiento más interesante del papel es que el volumen no es lo mismo que la claridad.

• La analogía: Imagina que tienes una radio con mala señal. Puedes subir el volumen al máximo (aumentar la amplitud) y escucharás un estruendo fuerte, pero seguirás sin entender la canción.
• Lo que descubrieron: Los métodos tradicionales de procesamiento sísmico a menudo "suben el volumen" de las frecuencias altas, haciendo que los datos parezcan más nítidos visualmente. Pero, al medir la Varianza de Fase, descubrieron que, aunque el volumen subió, las "agujas de los relojes" seguían apuntando a todas partes.
  • Conclusión: Puedes tener una señal muy fuerte y muy ruidosa al mismo tiempo. La nueva métrica nos dice: "Oye, puedes usar estas frecuencias bajas porque la fase es estable, pero ignora esas frecuencias altas porque, aunque suenan fuertes, son solo ruido".

5. ¿Por qué es importante esto?

Antes, los geofísicos tenían que adivinar si sus datos eran buenos o malos. Ahora tienen un termómetro objetivo:

1. Diagnóstico automático: Pueden escanear todo el volumen de datos y decir exactamente dónde la información es confiable y dónde es basura, sin depender de la vista humana.
2. Ahorro de tiempo y dinero: Evitan intentar hacer análisis complejos (como inversión de ondas completas) en zonas donde la fase es tan desordenada que los resultados serían falsos.
3. Mejor imagen: Al saber exactamente dónde está el "ruido de fase", pueden aplicar filtros más inteligentes para limpiar solo lo que realmente necesita limpieza.

En resumen

Este artículo nos dice que para ver el subsuelo con claridad, no basta con hacer el sonido más fuerte. Necesitamos asegurarnos de que el ritmo (la fase) de las ondas esté sincronizado.

La Varianza de Fase es la herramienta que nos permite medir ese ritmo. Es como pasar de mirar una foto borrosa y decir "se ve un poco mejor" a tener un informe científico que dice: "El 80% de la imagen es nítida, pero el 20% superior es un caos total, así que no intentes leer lo que hay ahí".

Es un cambio de paradigma: dejar de mirar solo la intensidad del sonido y empezar a medir la confianza en el momento exacto en que ocurre.

Resumen técnico

Resumen Técnico: Varianza de Fase como Atributo de Control de Calidad Sísmico

1. Planteamiento del Problema

Las ondas sísmicas registradas en tierra sufren distorsiones severas debido a la heterogeneidad de la superficie cercana (near-surface). Estas heterogeneidades introducen perturbaciones de fase específicas por traza y dependientes de la frecuencia que persisten incluso después de aplicar procesos avanzados de tiempo.

• Limitaciones de los métodos convencionales: El procesamiento sísmico tradicional se basa en la deconvolución consistente con la superficie y correcciones estáticas residuales. Estos métodos asumen que las distorsiones son consistentes con la superficie (repetibles entre fuentes y receptores) y tratan la variabilidad de fase como un problema de desplazamiento temporal (cinemático). Sin embargo, la heterogeneidad a pequeña escala induce dispersión hacia adelante e interferencias que generan perturbaciones de fase no consistentes con la superficie y dependientes de la frecuencia. Estas no pueden corregirse solo con desplazamientos de tiempo.
• Falta de métricas objetivas: Actualmente, no existe una medida directa y cuantitativa de la fiabilidad de la fase. La calidad de la fase se evalúa indirectamente mediante el comportamiento de la amplitud o la inspección visual, lo que deja sin diagnosticar el desorden de fase residual, especialmente en entornos de alta dispersión (como desiertos o rocas cristalinas).
• El problema de la "fase envuelta": La fase sísmica es una variable circular (rango $[-\pi, \pi]$). Los métodos lineales tradicionales fallan al tratar con ruido o incoherencia, y el "desenrollado" (unwrapping) de fase es inestable en regiones ruidosas.

2. Metodología Propuesta

Los autores proponen un nuevo marco basado en estadística circular para tratar la fase sísmica como una variable aleatoria circular dentro de conjuntos de trazas locales (ensembles).

• Enfoque Estocástico: En lugar de estimar la fase traza por traza, se modela la fase en cada frecuencia como una muestra de una distribución subyacente dentro de un conjunto de trazas vecinas.
• Definición de la Varianza Circular ($V$):
  • Se normalizan los espectros complejos de las trazas a magnitud unitaria para aislar la información de fase.
  • Se calcula la longitud del vector resultante medio ($\bar{R}$), que mide la alineación de los vectores de fase unitarios.
  • La Varianza Circular se define como: $V(\omega) = 1 - \bar{R}(\omega)$.
  • Interpretación:
    • $V \approx 0$: Coherencia de fase perfecta (todas las trazas alineadas).
    • $V \approx 1$: Incoherencia total (ruido puro, distribución uniforme en el círculo).
• Ventajas del método:
  • No requiere desenrollado de fase (funciona directamente con fases envueltas).
  • Es independiente de modelos de onda o suposiciones globales.
  • Proporciona una clasificación traza a traza y frecuencia a frecuencia.
  • Utiliza la distribución de von Mises como modelo conceptual (análogo circular a la distribución normal) para interpretar la dispersión, aunque el cálculo de $V$ es no paramétrico.

3. Contribuciones Clave

1. Nuevo Atributo de Control de Calidad (QC): Introducción de la "varianza de fase" como una métrica cuantitativa, automática y objetiva para evaluar la fiabilidad de la fase en datos pre-apilados (prestack).
2. Cambio de Paradigma: Pasar de ver la variabilidad de fase como un artefacto de traza individual a una propiedad medible del conjunto (ensemble) que refleja la física de la dispersión por heterogeneidades.
3. Definición de Banda de Frecuencia Efectiva: Capacidad de distinguir entre el aumento de amplitud (bandwidth) y la coherencia real de la fase. Permite definir qué frecuencias son útiles para flujos de trabajo sensibles a la fase (como Inversión de Forma de Onda Completa - FWI, o AVO).
4. Validación Sintética y de Campo: Demostración de que la varianza circular recupera fielmente perturbaciones de fase impuestas y revela limitaciones ocultas en el procesamiento convencional.

4. Resultados

• Pruebas Sintéticas: En experimentos controlados donde se impusieron perturbaciones de fase conocidas (variando la varianza circular), el atributo calculado coincidió casi perfectamente con los valores impuestos, validando su precisión para cuantificar el ruido de fase dependiente de la frecuencia.
• Aplicación a Datos de Campo (Tierra):
  • Análisis Pre y Post-Procesamiento: Al comparar datos prestack antes y después del procesamiento convencional (deconvolución, estáticas, atenuación de ruido), se observó que:
    • El procesamiento reduce la varianza de fase principalmente en frecuencias bajas a intermedias y en offsets medios a lejanos.
    • En frecuencias altas, la varianza de fase permanece alta (cerca de 1.0), indicando que la fase sigue dominada por ruido, a pesar de que los espectros de amplitud muestren una "extensión" de la banda de frecuencias.
    • En los offsets cercanos (dentro del cono de ruido), la varianza sigue siendo muy alta incluso después del procesamiento.
  • Comparación con SNR: Mientras que la relación señal-ruido (SNR) basada en apilamiento muestra una mejora general, la varianza de fase revela que esta mejora no es uniforme en el espectro de frecuencias. Las métricas basadas solo en amplitud pueden sobreestimar la banda de frecuencias utilizable.

5. Significado e Impacto

El trabajo establece las estadísticas circulares como una base fundamental para el análisis de fase sísmica moderna, especialmente en la era de la adquisición con receptores puntuales (point-receivers) y muestreo denso, donde la variabilidad de fase es más evidente.

• Mejora en la Toma de Decisiones: Permite a los geofísicos definir objetivamente la "banda de frecuencia efectiva" basada en la coherencia de fase, no solo en la energía. Esto es crucial para evitar inyectar ruido de fase en procesos sensibles como la Inversión de Forma de Onda Completa (FWI) o el análisis AVO.
• Guía para el Procesamiento: Ofrece una herramienta para monitorear paso a paso el flujo de trabajo de procesamiento, identificando qué pasos mejoran o degradan la integridad de la fase.
• Habilitación de Nuevos Métodos: Facilita el desarrollo de técnicas de enmascaramiento y condicionamiento de fase basadas en objetivos explícitos (ej. alcanzar un umbral de varianza específico) en lugar de criterios visuales subjetivos.

En conclusión, la varianza de fase transforma la evaluación de la calidad de los datos sísmicos de un arte subjetivo a una disciplina estadística rigurosa, revelando limitaciones físicas que los métodos tradicionales de amplitud y SNR pasan por alto.