Diseño de un filtro inverso para la atenuación del ruido en imágenes sísmicas

Abstract

A partir de los modelos de velocidad de onda P y de densidad, se calcularon datos sintéticos representativos de un área con tectonismo salino intenso en aguas profundas del Golfo de México Norte. Para este fin, se simuló la propagación de ondas acústicas a través de un domo de sal aislado, perteneciente a este medio geológico. El algoritmo de simulación numérica de la propagación de ondas acústicas, utiliza el método de diferencias finitas de segundo orden. Además, se agregaron fronteras absorbentes tipo CPML, para eliminar las difracciones espurias, no físicas, provenientes de las fronteras del dominio computacional. De esta forma, se obtuvo el registro sísmico del medio geológico previamente descrito y se identificaron los eventos múltiples y primarios que generan las interfaces del modelo, correspondientes a las unidades litológicas principales. Posteriormente, se aplicó la transformada de Hough, para detectar automáticamente los múltiples mediante el reconocimiento de patrones parabólicos que producen el mejor ajuste en la zona de mayor amplitud. Empleando lógica difusa, se implementaron filtros en las regiones delimitadas por las parábolas identificadas. Finalmente, se empleó un artificio numérico, consistente en emplear una frontera absorbente superior, cercana a la superficie libre, para suprimir por completo los múltiples asociados con dicha superficie y conseguir una imagen libre de estos efectos. En la etapa final, para la remoción de los múltiples se aplica un cociente que atenúa la amplitud de los eventos en las regiones previamente delimitadas empleando lógica difusa y funciones de membresía, que acentúan los efectos del filtrado en función de la amplitud. Por otro lado, al incorporar una frontera superior tipo CPML en la simulación numérica de la propagación de ondas acústicas, se consiguió eliminar los múltiples asociados con la superficie libre.