Nuevos métodos y tecnologías en el tratamiento de rutina
El procesamiento de corrección dinámica y estática se ha desarrollado desde el procesamiento de puntos de muestra completa hasta métodos de corrección de alta precisión (o alta fidelidad). La idea básica de este método es calcular los valores de corrección sin redondear, calcular los valores de muestra en múltiplos no enteros de intervalos de muestreo mediante métodos de interpolación de alta precisión y luego corregir a múltiplos enteros de puntos de muestreo.
Para superar la distorsión de tensión después de la corrección dinámica, se propone un método de corrección dinámica de la onda de reflexión en movimiento general.
La corrección estática se ha desarrollado desde la corrección estática de campo hasta la corrección estática de refracción y la corrección estática cromatográfica. La corrección estática tomográfica es un método de corrección estática de rayos de flexión que tiene una gran adaptabilidad a la complejidad del medio de la superficie.
En el método de superposición convencional, se añaden la superposición que preserva la amplitud y la superposición ponderada adaptativa. Este método permite calcular y superponer automáticamente coeficientes de peso en función de la calidad del registro sísmico de cada canal y de cada momento. Se ha desarrollado un método que puede reemplazar la superposición sin superposición: el ajuste polinomial AVA (o AVO). El resultado adecuado de este método es directamente el coeficiente de reflexión con compensación cero, que tiene una alta resolución en relación con la superposición.
3.6.1.2 El filtrado digital mejora la relación señal-ruido
El filtrado digital tradicional se basa principalmente en el filtrado de frecuencia unidimensional y el filtrado F-k bidimensional. Los filtros digitales recientemente desarrollados incluyen: filtrado por división de frecuencia con transformada de ondas, filtrado por transformada K-L, filtrado por descomposición SVD, filtrado de ajuste polinomial, filtrado de dominio F-x, filtrado por transformada CT, filtrado de continuación de campo de ondas y filtrado de mediana. Algunos de estos filtros son lineales y otros no lineales. Con el desarrollo de la teoría y la tecnología del análisis de señales, surgirán nuevos métodos de filtrado.
3.6.1.3 Filtrado inverso para mejorar la resolución vertical
El filtrado inverso de datos sísmicos es un punto difícil en el procesamiento de datos sísmicos. La razón principal que afecta el efecto de filtrado inverso es que los supuestos teóricos son inconsistentes con la situación real. Si se pudieran reducir los requisitos de condición o se pudiera determinar una ondícula sísmica real, se mejoraría el efecto de deconvolución, por lo que se desarrollaron una serie de nuevos métodos de deconvolución.
Para cambiar los supuestos del filtrado inverso, se desarrollaron el filtrado inverso homomórfico y el filtrado inverso de mínima entropía. Todos abandonan la suposición de que las ondas sísmicas son ondas de fase mínima y las secuencias de coeficientes de reflexión son ruido blanco. El filtrado inverso homomórfico transforma los registros del dominio del tiempo en el dominio de la frecuencia, realiza operaciones logarítmicas y obtiene espectros logarítmicos y secuencias de espectros logarítmicos, convirtiendo la relación de convolución en una relación de suma simple. Cuando la secuencia espectral logarítmica de la wavelet y la secuencia espectral logarítmica del coeficiente de reflexión se distribuyen en diferentes ubicaciones (nueva hipótesis), se pueden separar fácilmente. La idea básica del filtrado inverso de entropía mínima es encontrar un operador de filtro inverso de modo que la salida de su convolución con el registro sísmico sea un pico grande con la forma más simple y signo y posición desconocidos (mejorando así la resolución), es decir , minimizando la señal de salida de "entropía" (la entropía es un término tomado de la termodinámica). En teoría de la información, la entropía mínima significa que la señal tiene el mejor orden, la regularidad más fuerte y la menor incertidumbre. La premisa de este método requiere que el espacio entre las capas reflectantes sea grande y desigual, y que la forma de la ondícula sísmica permanezca sin cambios.
A través del estudio del mecanismo de filtrado geodésico, nos dimos cuenta de que debido al efecto de transformación del filtrado geodésico sobre las ondas sísmicas, las ondas sísmicas varían en el tiempo (diferentes profundidades y capas poco profundas), por lo que se utilizó un método para la inversa. Se propuso un filtrado inverso que varía por partes en el tiempo. Reconociendo que el filtrado terrestre es causado principalmente por la absorción del medio, se propone el método de filtrado Q inverso para compensar la absorción. Este es en realidad un filtrado inverso que varía en el tiempo, pero es mucho más refinado que el inverso segmentado que varía en el tiempo. filtración.
Sabemos que mientras conozcamos la ondícula sísmica, es fácil encontrar el factor de filtrado inverso mediante cualquier método. Por lo tanto, otra dirección de desarrollo del filtrado inverso es la tecnología de extracción de ondas. En la exploración costa afuera, la instalación de un geófono cerca de la fuente del terremoto para registrar directamente las ondas y el uso de datos de registro de pozos y perfiles sísmicos verticales para obtener ondas son métodos para obtener ondas de fuentes de información distintas al registro convencional, mientras que los métodos de promedio estadístico multicanal, automáticos Los métodos de correlación y los métodos de búsqueda de raíces polinómicas son métodos para extraer directamente wavelets de registros convencionales. En términos generales, las ondas sísmicas obtenidas mediante el primer método son confiables, mientras que el segundo método a menudo necesita hacer algunas aproximaciones o introducir suposiciones.
Tanto la teoría como la práctica han demostrado que entre las wavelets con el mismo espectro de amplitud, las wavelets de fase cero tienen la resolución más alta. Por lo tanto, existe un método de filtrado inverso para mejorar la resolución vertical, en lugar de obtener la secuencia del coeficiente de reflexión, convierte la wavelet de fase distinta de cero en el registro en una wavelet de fase cero, lo que mejora la resolución sin reducir la señal. -Relación de ruido. Las ricas características dinámicas aún se conservan para la interpretación. Este método se llama procesamiento wavelet. Aunque existen muchos métodos, son esencialmente filtrado inverso.
3.6.1.4 Migración de imágenes para mejorar la resolución lateral
En los últimos años, con el desarrollo de una gran capacidad de almacenamiento y computadoras de alta velocidad, la tecnología de migración de imágenes ha logrado grandes avances. La dinámica principal es desde imágenes de migración bidimensionales a pseudo-3D (tridimensionales de dos pasos) a verdaderas (completas) tridimensionales, desde la migración posterior a la pila hasta la migración previa a la pila, y desde la migración temporal a la migración en profundidad; La migración en profundidad tridimensional previa al apilamiento se ha utilizado ampliamente en la producción. Dado que la migración de profundidad previa al apilamiento 3D es el retorno del campo de ondas de tiempo 3D completo que tiene en cuenta la refracción de rayos de la interfaz, las imágenes estructurales 3D son más realistas y confiables, la estructura es más clara y la precisión de la exploración sísmica mejora enormemente. Otro desarrollo de las imágenes de migración son las imágenes de migración de múltiples componentes.
La migración en profundidad previa al apilamiento requiere un modelo de velocidad confiable. Con el desarrollo simultáneo de la tecnología de migración, han surgido métodos para establecer modelos de velocidad de migración.
3.6.1.5 Extracción de parámetros de velocidad
En los últimos años, la tecnología de extracción de parámetros de velocidad también se ha desarrollado rápidamente.
En términos de análisis de velocidad, ha surgido un método de análisis de velocidad adecuado para interfaces inclinadas, llamado análisis de velocidad continuo de escaneo inclinado. Se ha desarrollado un método de análisis de velocidad adecuado para cualquier situación de interfaz compleja. La idea principal es que eventualmente se unificarán los cálculos iterativos directos e inversos y el ajuste de hipérbola. Además, también se ha estudiado en profundidad la interpretación automática del espectro de velocidad. Constantemente surgen métodos para analizar automáticamente los datos del espectro de velocidad utilizando varias herramientas matemáticas más recientes, como el reconocimiento de patrones, la teoría de grafos y las matemáticas difusas, y se han logrado ciertos resultados.
En el cálculo de la velocidad de la capa, se implementa el método de calcular la velocidad de la capa capa por capa utilizando la transformación τ-p. La velocidad de la capa también se puede calcular utilizando el método de continuación de la ecuación de onda, que es más preciso. Utiliza características dinámicas como la amplitud y la forma de onda de las ondas sísmicas. Particularmente dignos de mención son los diversos métodos de inversión de velocidad para ecuaciones ondulatorias y no ondulatorias desarrollados en los últimos años, como la inversión de Bourne, la inversión de Leitoff, la inversión lineal generalizada, la inversión tomográfica, la tecnología de reconstrucción algebraica, etc. Usando estos métodos, la estructura de velocidades finas del medio subterráneo, que está estrechamente relacionada con la estructura y la litología, se puede obtener directamente de los registros sísmicos.