Método de procesamiento de datos sísmicos que proporciona una captación de imágenes comunes de desplazamiento de superficie.

Un método de procesamiento de datos sísmicos, que comprende:

- introducir unos datos sísmicos

(DS,G[t]) que incluyen, para cada una de las mismas una pluralidad de tomas en las localizaciones de origen respectivas, unas trazas sísmicas registradas en una pluralidad de localizaciones de receptor;

- aplicar un proceso de migración en profundidad a los datos sísmicos para obtener un primer conjunto de datos migrados que incluyen, para cada toma, los primeros valores migrados (MS[x, y, z]), asociados respectivamente con una pluralidad de posiciones del subsuelo;

- calcular los datos modulados desplazados (D'S,G[t]) multiplicando los datos sísmicos en cada traza sísmica por un desplazamiento horizontal entre las localizaciones de origen y de receptor para dicha traza sísmica;

- aplicar el proceso de migración en profundidad a los datos modulados desplazados para obtener un segundo conjunto de datos migrados que incluyen, para cada toma, unos segundos valores migrados (M'S[x, y, z]) asociados respectivamente con la pluralidad de posiciones de subsuelo;

- estimar, para cada toma, los valores de desplazamiento (hS[x, y, z]) asociados respectivamente con al menos alguna de las posiciones de subsuelo, por un proceso de división aplicado a los conjuntos de datos migrados primero y segundo; y

- estimar una captación de imágenes comunes en una posición horizontal, que comprende los valores migrados respectivos (Rx,y[z,h]) para los pares de parámetros incluyendo cada uno de los mismos un parámetro de profundidad y un parámetro de desplazamiento,

en el que el valor migrado para un par de parámetros en la captación de imágenes comunes en dicha posición horizontal es un primer valor migrado del primer conjunto de datos migrados asociados con una posición de subsuelo determinada por dicha posición horizontal y el parámetro de profundidad de dicho par de parámetros para una toma de tal manera que el valor de desplazamiento estimado asociado con dicha posición de subsuelo es el parámetro de desplazamiento de dicho par de parámetros.

Tipo: Patente Internacional (Tratado de Cooperación de Patentes). Resumen de patente/invención. Número de Solicitud: PCT/IB2011/001585.

Solicitante: Total SA.

Nacionalidad solicitante: Francia.

Dirección: 2 place Jean Millier, La Défense 6 92400 Courbevoie FRANCIA.

Inventor/es: BAINA,REDA, NICOLETIS,LAURENCE, GIBOLI,MATTEO.

Fecha de Publicación: .

Clasificación Internacional de Patentes:

  • SECCION G — FISICA > METROLOGIA; ENSAYOS > GEOFISICA; MEDIDA DE LA GRAVITACION; DETECCION DE... > Sismología; Prospección o detección sísmica o... > G01V1/28 (Tratamiento de los datos sísmicos, p. ej. para análisis, para interpretación, para corrección (G01V 1/48 tiene prioridad))
google+ twitter facebookPin it
Ilustración 1 de Método de procesamiento de datos sísmicos que proporciona una captación de imágenes comunes de desplazamiento de superficie.
Ilustración 2 de Método de procesamiento de datos sísmicos que proporciona una captación de imágenes comunes de desplazamiento de superficie.
Ilustración 3 de Método de procesamiento de datos sísmicos que proporciona una captación de imágenes comunes de desplazamiento de superficie.
Ilustración 4 de Método de procesamiento de datos sísmicos que proporciona una captación de imágenes comunes de desplazamiento de superficie.
Ilustración 5 de Método de procesamiento de datos sísmicos que proporciona una captación de imágenes comunes de desplazamiento de superficie.
Ilustración 6 de Método de procesamiento de datos sísmicos que proporciona una captación de imágenes comunes de desplazamiento de superficie.
Método de procesamiento de datos sísmicos que proporciona una captación de imágenes comunes de desplazamiento de superficie.

Texto extraído del PDF original:

DESCRIPCIÓN

Método de procesamiento de datos sísmicos que proporciona una captación de imágenes comunes de desplazamiento de superficie

Sector de la técnica La presente invención se refiere a unos métodos de procesamiento de datos sísmicos. Es aplicable en el campo de la prospección geofísica cuando se necesitan imágenes del subsuelo.

Estado de la técnica Para obtener imágenes del subsuelo, a menudo se usa un método sísmico, que consiste en crear y enviar las ondas sísmicas en el suelo usando fuentes tales como explosivos o camiones vibradores en tierra, o armas de aire comprimido en alta mar. Las ondas sísmicas atraviesan el suelo y rebotan o se reflejan en las discontinuidades geológicas en el subsuelo. Como resultado de ello, vuelven a la superficie, donde se registran usando matrices de tres geófonos (en tierra), o hidrófonos (en alta mar) de componente que se distribuyen de manera regular para cubrir zonas de varios kilómetros cuadrados.

La figura 1 ilustra en forma de diagrama un estudio de datos sísmicos con una fuente S de ondas sísmicas y una matriz de receptores G. Se muestra también un punto B del subsuelo que se supone que contribuye a la señal detectada por uno de los receptores G. Las coordenadas horizontales del punto B del subsuelo se indican por x, y (o una sola coordenada espacial si se consideran imágenes en 2D en lugar de imágenes en 3D), mientras que su profundidad se indica por z. La figura 1 proporciona también una representación simplificada (líneas discontinuas) de la propagación de las ondas sísmicas desde la fuente S al punto B y desde el punto B al receptor G. Las ondas se refractan en las discontinuidades de las capas geológicas en las que la impedancia acústica cambia y se refleja o se difracta en diferentes posiciones, incluyendo la del punto B. Los datos registrados en un estudio sísmico incluyen, para cada toma de una fuente S y para cada receptor G, una traza sísmica que es una serie temporal de la señal detectada por el receptor G. La traza para un número de tomas debe transformarse para proporcionar una imagen del subsuelo que será el resultado de apilar o integrar una gran cantidad de información. Una etapa importante de la transformación es la migración que consiste en la reordenación de los datos con respecto a un modelo, de tal manera que el apilamiento puede realizarse de manera coherente. El modelo es, en general, un mapa de la velocidad de propagación de las ondas acústicas en el subsuelo. No se conoce a priori y es un reto principal de todas las tecnologías de formación de imágenes sísmicas determinar un modelo que explicará correctamente los datos de campo después del apilamiento. En los métodos de migración de profundidad de pre-apilamiento (PSDM), se calculan los datos migrados para cada toma usando el modelo de velocidad y se disponen en un cubo de salida que contiene los valores migrados asociados con las posiciones en el subsuelo: Los cubos obtenidos para las diferentes tomas se analizan a continuación para comprobar la consistencia del modelo. El modelo puede corregirse y el proceso se repite hasta que se obtiene una imagen satisfactoria. Las captaciones de imágenes comunes (CIG) son herramientas populares para evaluar el campo de velocidades de migración y para la mejora de imágenes. Se hacen a partir de los datos extraídos de los cubos de salida, ordenados de una forma conveniente para el análisis con el fin de comprobar el modelo de velocidad. Un CIG es una estructura de datos bi-dimensional definida para una posición horizontal dada x, y, con un primer eje que representa la profundidad z y un segundo eje que representa un parámetro de dominio A contemplado para ordenar los datos de los cubos de salida. Contiene valores de reflectividad obtenidos a partir de los cubos de salida resultantes de la migración, que forman una imagen que puede analizarse para comprobar y/o corregir el modelo de velocidad. En esta imagen, un valor de píxel en un punto (z, A) representa un valor migrado derivado como una contribución de la posición de subsuelo x, y, z para una traza sísmica asociada con el parámetro de dominio A. Ejemplos de parámetros de dominio A usados incluyen el desplazamiento de superficie, es decir, la distancia entre la localización de origen para una toma y la localización del receptor que proporciona la traza pertinente para esa toma, o el ángulo de dispersión en la posición de subsuelo x, y, z. El cálculo de las captaciones de imágenes comunes no es sencillo en todos los métodos de extrapolación de campo de ondas. De hecho, hasta la fecha, los cubos de salida se ordenan de diferentes maneras en función del método de propagación empleado para realizar la migración. Incluso si no hay ninguna razón teórica para esto, la manera más práctica para producir las captaciones de imágenes comunes se selecciona de acuerdo con los diferentes tipos de migración. Por ejemplo, las CIG de ángulo de dispersión o las CIG de desplazamiento de superficie son de uso generalizado en las técnicas de tomografía basadas en rayos, para buscar el modelo de velocidad que proporcione las ‘captaciones planas’. Las figuras 2A-C ilustran el principio en el caso de las CIG de desplazamiento de superficie. Los inventores consideran un reflector plano en una posición x0, y0, z0 del subsuelo. Si el modelo de velocidad se estima de manera correcta, el reflector proporcionará un pico en la señal migrada en la profundidad z0 para todos los valores del desplazamiento h, proporcionando de este modo un evento plano en la CIG como se muestra en la figura 2A. Sin embargo, si se sobreestima la velocidad de migración en el modelo, el mismo pico observado al mismo tiempo en una traza sísmica se asociará con un mayor valor de profundidad z1, como se muestra en la figura 2B, y el evento correspondiente en la CIG tendrá una forma cóncava ya que la discrepancia en el tiempo de viaje de la onda aumenta con el desplazamiento h. Asimismo, si se subestima la velocidad de migración en el modelo, el pico se asociará con un valor de menos profundidad z2, como se muestra en la figura 2C, y el evento correspondiente en la CIG tendrá una forma convexa. Si el reflector tiene un ángulo de inclinación distinto de cero, su imagen en la CIG se desplazará de manera horizontal además de verticalmente. Se han desarrollado diferentes herramientas para analizar las CIG de desplazamiento de superficie con el fin de corregir el modelo de velocidad. Sin embargo, estas herramientas se han usado principalmente en métodos de migración basándose en la estimación de los tiempos de viaje entre los reflectores y la superficie. Se han desarrollado los métodos de migración más sofisticados para construir imágenes PSDM resolviendo la ecuación de onda con el fin de obtener amplitudes de reflector más precisas y un posicionamiento estructural. Por ejemplo, la migración en tiempo inverso (RTM) es una solución de migración de dos vías que puede describir de manera precisa la propagación de ondas en un medio complejo. Se usa cada vez más en la exploración sísmica en virtud de los avances en la potencia de los ordenadores y la programación.

Las herramientas de análisis mencionadas anteriormente no se usan con métodos PSDM de ecuación de onda, incluyendo el RTM, porque no se conoce cómo calcular las CIG de desplazamiento de superficie. En “Offset and angle-domain common image-point gathers for shot-profile migration”, Geofísica, vol. 67, Nº 3, 2002, pág. 883-889, J. Rickett y P. Sava establecieron la noción de las CIG de desplazamiento de subsuelo que requieren la extensión de la condición de formación de imágenes a través del cálculo de la función de correlación a lo largo de la dimensión horizontal espacial. Este tipo de captaciones es la forma más común de imágenes de migración basadas en el campo de onda de salida. Es el más adecuado para las técnicas basadas en el análisis de enfoque, que buscan la más alta correlación en el retraso de tiempo cero y/o de desplazamiento cero. En “Angle-domain common image gathers by wavefield continuation methods”, Geofísica, vol. 68, Nº 3, 2003, pág. 1065- 1074, P. Sava y S. Fomel proponen un método para derivar las CIG de ángulo de dispersión a partir de las CIG de desplazamiento de subsuelo. Más recientemente, los mismos autores, combinando tanto el tiempo como la profundidad, han introducido el concepto de condición de imagen ampliada en “Time-shift imaging condition in seismic migration”, Geofísica, vol. 71, Nº 6, pág. 209-217.

A pesar de sus beneficios, algunos problemas impiden el uso generalizado del desplazamiento de subsuelo y las captaciones de ángulo. En primer lugar, ya que el RTM es un proceso intensivo de ordenador, el cálculo de las correlaciones cruzadas para todas las localizaciones de CIG seleccionadas añade un considerable coste adicional. En segundo lugar, las captaciones de desplazamiento de subsuelo no pueden emplearse para la tomografía clásica. Los analistas no pueden beneficiarse de todo el arsenal de herramientas desarrolladas para las CIG de superficie o de ángulo clásicas (como RMO, Mute, el análisis AVO/AVA, etc.). Las CIG de ángulo de dispersión de subsuelo serían adecuadas para este objetivo. Sin embargo, estas captaciones implican un coste adicional para una transformada de Radón de dos parámetros, y no muestran el movimiento de salida cinemático correcto. ADLER F ET AL: “Nonlinear 3d tomographic least-squares inversion of residual moveout in Kirchhoff prestackdepth- migration common-image gathers”, GEOFÍSICA, SOCIEDAD DE GEOFÍSICOS DE EXPLORACIÓN, de Estados Unidos, vol. 73, Nº 5, sup, 1 de Septiembre de 2008 (2008-09-01), páginas VE13-VE23 describen un método de procesamiento de datos sísmicos, que comprende:  introducir los datos sísmicos que incluyen, unas trazas sísmicas registradas en una pluralidad de localizaciones de receptor;  aplicar un proceso de migración en profundidad a los datos sísmicos para obtener un primer conjunto de datos migrados que incluyen, los primeros valores migrados asociados respectivamente con una pluralidad de posiciones de subsuelo;  calcular un conjunto de datos migrados de profundidad modelados para obtener un segundo conjunto de datos migrados realizando una migración finita de mapa de desplazamiento;  equiparar el segundo conjunto de datos al primer conjunto de datos; y  estimar una captación de imágenes comunes. Sería deseable obtener captaciones de desplazamiento de superficie con diferentes tipos de métodos de migración, en particular con los métodos de campo de ondas que incluyen el RTM, con el fin de mantener las ventajas de los métodos de campo de ondas y, al mismo tiempo, enfrentarse a la limitación de la suposición asintótica de los métodos basados en rayos, mientras que ordenan los cubos migrados de la misma forma que las captaciones de desplazamiento de superficie clásicas.

Objeto de la invención

Se propone un método de procesamiento de datos sísmicos para proporcionar unas captaciones de imágenes comunes en el dominio de desplazamiento de superficie. El método comprende: - introducir unos datos sísmicos que incluyen, para cada una de una pluralidad de tomas en las localizaciones de origen respectivas, unas trazas sísmicas registradas en una pluralidad de localizaciones de receptor; - aplicar un proceso de migración en profundidad a los datos sísmicos para obtener un primer conjunto de datos migrados que incluyen, para cada toma, los primeros valores migrados asociados respectivamente con una pluralidad de posiciones del subsuelo; - calcular los datos modulados desplazados multiplicando los datos sísmicos en cada traza sísmica por un desplazamiento horizontal entre las localizaciones de origen y de receptor para dicha traza sísmica; - aplicar el proceso de migración en profundidad a los datos modulados desplazados para obtener un segundo conjunto de datos migrados que incluyen, para cada toma, unos segundos valores migrados asociados respectivamente con la pluralidad de posiciones de subsuelo; - estimar, para cada toma, los valores de desplazamiento asociados respectivamente con al menos alguna de las posiciones de subsuelo, por un proceso de división aplicado a los conjuntos de datos migrados primero y segundo; y - estimar una captación de imágenes comunes en una posición horizontal, que comprende los valores migrados respectivos para los pares de parámetros incluyendo cada uno de los mismos un parámetro de profundidad y un parámetro de desplazamiento.

El valor migrado para un par de parámetros en la captación de imágenes comunes en una posición horizontal es un primer valor migrado del primer conjunto de datos migrados asociados con una posición de subsuelo determinada por dicha posición horizontal y el parámetro de profundidad de dicho par de parámetros para un toma, de tal manera que el valor de desplazamiento estimado asociado con dicha posición de subsuelo es el parámetro de desplazamiento de dicho par de parámetros. Las captaciones de desplazamiento de superficie pueden obtenerse sin hacer referencia al paradigma de Kirchhoff tradicional que se basa en una aproximación de frecuencia infinita de la ecuación de onda y el trazado de rayos. El método puede aplicarse de este modo a diversos tipos de procesos de migración en profundidad, incluyendo los procesos de migración en profundidad de pre-apilamiento de ecuación de onda (WEPSDM). En una realización preferida, el proceso de migración en profundidad es un proceso de migración de tiempo inverso (RTM). El RTM proporciona un mejor propagador que el trazado de rayos. Por lo tanto, permite la producción de las CIG en áreas complejas como de sub-sal, de sobrecarga compleja, etc., algo más cerca de las captaciones de desplazamiento “verdaderas”. Además, ya que la ordenación de la salida es la misma, las captaciones de desplazamiento de Kirchhoff convencionales pueden reemplazarse por las CIG obtenidas por medio del método como la entrada de la mayoría de las herramientas de post-procesamiento convencionales. Debido a años de conocimiento en la migración de tipo Kirchhoff, este método permite una interpretación más fácil de los resultados migrados. En comparación con las CIG de subsuelo, requiere también mucho menos esfuerzo computacional para calcular las captaciones. Además, estas captaciones son potencialmente adecuadas para la tomografía de trayectoria de onda. Para lograr una estabilidad mejorada, el proceso de división usado para estimar los valores de desplazamiento asociados con una posición de subsuelo puede comprender, minimizar una función de coste definida por una variable de desplazamiento y los valores locales de los valores migrados primero y segundo en una vecindad de la posición de subsuelo. Otro aspecto de la invención se refiere a un sistema para procesar los datos sísmicos, que comprende un recurso informático configurado para realizar un método de análisis de datos sísmicos como se ha definido anteriormente. Sin embargo, otro aspecto de la invención se refiere a un producto de programa informático de un sistema para procesar unos datos sísmicos, que comprende unas instrucciones para realizar un método como se ha definido anteriormente cuando dicho producto de programa se ejecuta en una unidad de procesamiento informático del sistema. Otras características y ventajas del método y el aparato divulgado en el presente documento serán evidentes a partir de la siguiente descripción de unas realizaciones no limitantes, con referencia a los dibujos adjuntos.

Descripción de las figuras

La figura 1 es un diagrama esquemático que ilustra la adquisición de datos sísmicos. Las figuras 2A-C muestran unas CIG simplificadas calculadas usando diferentes valores de velocidad de migración.

La figura 3 es un diagrama de flujo de un método de un procesamiento de datos sísmicos de acuerdo con una realización de la invención. La figura 4 es un diagrama que ilustra la derivación de las CIG de desplazamiento de superficie de acuerdo con el método. Las figuras 5A-5B muestran los datos de entrada para el proceso de migración en un ejemplo sencillo simulado por medio del método propuesto.

Las figuras 6A-6B muestran los resultados de la migración en este ejemplo, y la figura 6C muestra la estimación de atributo resultante de la división de los resultados. La figura 7 es un diagrama que muestra parte de una CIG de desplazamiento de superficie obtenida aplicando el método en el ejemplo anteriormente mencionado.

Descripción detallada de la invención Una manera de obtener las CIG de desplazamiento de superficie para los métodos de migración de WEPSDM o de RTM sería calcular un migración por toma y por traza, necesitando un número de migraciones igual al número total de tomas medidas del número medio de receptores por toma. Esto es claramente poco práctico por el momento, especialmente para los casos de 3D. En su lugar, se propone usar una solución más factible, concretamente, la migración de atributo, también llamada migración doble. En el método de migración doble como se introdujo por N. Bleistein (“On the imaging of reflectors in the earth”, Geofísica, Vol. 52, Nº 7, julio 1987, pág. 931-942), se calculan dos migraciones con los mismos datos, la segunda implicando un operador de migración multiplicado por el ángulo de reflexión especular. A continuación, la división de las dos imágenes migradas proporciona el ángulo especular a lo largo de los reflectores. Un método similar puede usarse con el desplazamiento de superficie en lugar del ángulo de reflexión especular como el atributo migrado. La migración puede realizarse usando diversos métodos de PSDM incluyendo el RTM de toma de registro convencional. La metodología comprende entonces: - realizar una primera migración de RTM convencional (etapa 20 en la figura 3); - realizar una segunda migración de RTM de atributo de desplazamiento, en la que los datos se multiplican por los valores de desplazamiento antes de la migración (etapas 30 y 40 en la figura 3); - realizar una división de los dos datos migrados obtenidos, en un sentido de mínimos cuadrados, para obtener el desplazamiento (etapa 50 en la figura 3); - añadir la reflectividad obtenida a partir de los primeros datos migrados al subsuelo correspondiente y la localización de captación de desplazamiento proporcionada por el mapa de atributos obtenido para reconstruir las CIG de desplazamiento de superficie (etapa 60 en la figura 3).

De esta manera, los inventores son capaces de obtener las CIG de desplazamiento de superficie clásicas usando un método de migración de WEPDSM tal como de RTM, que es el mejor método de extrapolación disponible hoy en día para la migración sísmica. En la figura 3, se observa la entrada de trazas sísmicas en la etapa 10 a partir de las mediciones de campo DS,G[t], en la que S indica una localización de origen, G indica una localización de receptor y t es el tiempo. Cada traza se modula en la etapa 30 multiplicándola por el valor de desplazamiento correspondiente, concretamente el desplazamiento horizontal HS,G entre la localización de origen S y la localización de receptor G. Las trazas moduladas de desplazamiento son D’S,G[t] = HS,G.DS,G[t].

Los datos sísmicos DS,G[t] y los datos modulados de desplazamiento D’S,G[t] se migran, respectivamente, en las etapas 20 y 40 para proporcionar los datos de PSDM MS[x, y, z] y M'S[x, y, z]. El primer conjunto de datos migrados producidos en la etapa 20 incluye, para cada toma en una localización de origen S, un cubo de los valores migrados MS[x, y, z] asociado a las posiciones de subsuelo x, y, z. Del mismo modo, el segundo conjunto de datos migrados obtenidos en la etapa 40 usando el mismo proceso de migración en profundidad incluye otro cubo de valores migrados M'S[x, y, z] para cada toma.

Con el fin de estimar un valor de desplazamiento ĥS[x, y, z] para un valor migrado MS[x, y, z], es decir, un valor para una toma S y una posición de subsuelo x, y, z, se realiza un proceso de división en la etapa 50 para evaluar MS[x, y, z]/ M’S[x, y, z].

Una división sin procesar de los dos números puede dar lugar a problemas de estabilidad. En cambio, puede ser mejor lanzar la división como un conjunto de problemas de mínimos cuadrados locales. A continuación, el valor de desplazamiento ĥS[x, y, z] se encuentra minimizando una función de coste JS, x, y, z(h) definida en una vecindad  (x, y, z) centrada en la localización x, y, z. Una posible expresión de la función de coste JS, x, y, z(h) es: en la que es la señal analítica de la reflectividad, indicando H la transformada de Hilbert, y El tamaño de la vecindad  (x, y, z) es variable y puede estar en función de la aplicación. Se selecciona de tal manera que el valor del atributo migrado h puede suponerse razonablemente que sea constante a lo largo de  (x, y, z) para un toma dada.

Los valores de desplazamiento ĥS[x, y, z] pueden entonces estar contenidos como: en la que el superíndice (.)* indica la hermitiana. Téngase en cuenta que la parte superior del lado a mano derecha de la ecuación (2) es la correlación cruzada de las dos imágenes, la imagen migrada convencional y la imagen migrada de atributo. La parte inferior es la auto-correlación de la imagen migrada convencional, o la envoltura. De esta manera, se aumenta la estabilidad de la división.

Los valores de desplazamiento ĥS[x, y, z] obtenidos de este modo se usan para mapear los valores de reflectividad MS[x, y, z] a los desplazamientos correspondientes h, cuyos valores de reflectividad pueden, a continuación, disponerse como los CIG de desplazamiento de superficie en las posiciones horizontales x, y: Este proceso de calcular los CIG de desplazamiento de superficie se ilustra en el diagrama de la figura 4. El atributo de desplazamiento se obtiene por la división mapeada de la “migración doble”. Los dos paneles superiores del diagrama representan dos imágenes migradas de toma común que muestran un reflector horizontal. La primera imagen se obtiene por medio de una migración de tiempo inverso de los datos sísmicos planos, mientras que la segunda es la salida de la migración en la que los datos se han multiplicado por el desplazamiento h. En este panel, las amplitudes varían de manera lateral a lo largo del reflector y son proporcionales a la posición de receptor y, por lo tanto, al desplazamiento. En particular, comenzando desde el lado izquierdo, la amplitud es negativa, llega a cero exactamente por debajo de la posición de origen y posteriormente se convierte en positiva.

Para una posición específica (x0, z0) en la imagen migrada de toma (en este caso, la posición horizontal x0 puede ser de 2D, con los componentes x e y), se tiene un valor específico de la reflectividad R. En la misma posición en la imagen migrada de atributo, el valor de la reflectividad es R multiplicada por el desplazamiento h, que es el desplazamiento de la pareja origen/receptor que ha generado la reflexión especular. El desplazamiento es simplemente la división de las dos cantidades. Con estos cuatro valores (x0, z0, h, R), se puede construir ahora el cubo de desplazamiento migrado. El panel inferior del diagrama de la figura 4 representa una sección del punto medio común en la localización x0. La contribución de la reflectividad R se añade a la posición cuyas coordenadas están dadas por la pareja profundidad /desplazamiento (z0, h).

Un ejemplo numérico puede presentarse para ilustrar el método descrito anteriormente. En este ejemplo, el modelo es una velocidad constante de 2D de 5 km de profundidad y 8 km en un rectángulo de desplazamiento horizontal. La velocidad se supone que es de 2000 m/s y los datos de reflexión-P se generan por medio del modelado de Born acústico que introduce ocho discontinuidades cada 500 m a partir de 1 kilómetro. Se simula una anomalía de amplitud frente a un desplazamiento (AVO), en la que la reflexión cambia su polaridad cada 2000 m en un desplazamiento. Aprovechando la geometría 1 D de facto, puede representarse un estudio sísmico completo usando un solo origen. Ese origen se localiza a 3000 m, a una profundidad de 5 m por debajo de la superficie. Los receptores se localizan a la misma profundidad, en el lado derecho del origen cada 25 m, empezando a partir de la posición del origen. La figura 5A muestra la captación de toma común de entrada DS,G[t], mientras que la figura 5B muestra los datos multiplicados de atributo D’S,G[t]. Obsérvese cómo las amplitudes se enfatizan por la modulación de los desplazamientos más grandes. Para mayor simplicidad, se emplea un propagador de dominio de frecuencia para la migración en la simulación. Se toma un modelo de velocidad constante para la propagación, en el que la velocidad es mayor que la velocidad verdadera e igual a 2200 m/s. Los resultados obtenidos se comparan con el movimiento de salida teórico que puede calcularse de manera analítica para los modelos simples. Además, puede comprobarse la detección correcta de la anomalía AVO. Para cada conjunto de datos y para cada toma, se obtienen dos imágenes. La figura 6A-B muestra los resultados migrados de toma para los datos proporcionados en los que el origen se localiza a 3000 m. Los datos migrados planos MS[x, y, z] se muestran en la figura 6A, mientras que los datos migrados de atributo M'S[x, y, z] se muestran en la figura 6B. Como era de esperar, las amplitudes se enfatizan por la modulación que va de la izquierda a la derecha de la adquisición. Este comportamiento puede verse en la figura 6B que representa el atributo obtenido después de la división de las dos primeras imágenes. La estimación del atributo se produce en los puntos sobre los reflectores para los cuales hay un par especular de rayos procedentes del origen y del receptor hasta el punto de la superficie en el que se toman las imágenes. En las localizaciones en las que esto no se cumple, es decir, en las que no puede identificarse una reflectividad diferente de cero, los valores del atributo no son significativos. Por medio del mapa de atributos obtenido, se añade a continuación la contribución de la reflectividad en la posición de desplazamiento correspondiente. Para obtener la imagen migrada de desplazamiento final, se suman las contribuciones a lo largo de todas las tomas.

La figura 7 muestra un zoom de la captación de desplazamiento de superficie obtenida a partir de una migración en profundidad d pre-apilamiento de ecuación de onda de geófono de toma con los mismos parámetros de modelo que en la sección anterior. En particular, representa dos reflectores localizados a 1500 m y 2000 m de profundidad. Debido al error en la velocidad de migración, se descolocan y exhiben un movimiento de salida no plano. Se muestran también en la figura 7, mediante líneas continuas, las curvas teóricamente predichas z(h). Como la figura demuestra, la equiparación es bastante buena y la anomalía AVO se detecta correctamente. El método propuesto para obtener las captaciones de desplazamiento de superficie clásicas para diversas técnicas de PSDM, incluida la migración inversa en tiempo, ha demostrado poderse aplicar de manera exitosa en el contexto de la propagación de banda limitada. Los métodos propuestos no están en función de una implementación específica del método de extrapolación de campo de onda, ya que puede realizarse después de la migración de registro de toma en el dominio de Fourier, en el dominio de tiempo, etc. Una ventaja de este método es que permite el uso de mejores propagadores que los rayos para propagar el campo de ondas, haciendo posible adoptar la totalidad del arsenal de herramientas convencionales para el post- procesamiento desarrollado para los esquemas de migración asintótica. En comparación, el desplazamiento de subsuelo y el desplazamiento de subsuelo en ángulo de dispersión, con el desplazamiento de superficie para las CIG de RTM son una manera menos costosa de calcular, y preservan el movimiento de salida cinemático. Las realizaciones del método descrito en el presente documento pueden implementarse en cualquier tipo de ordenador u ordenadores y los componentes pueden implementarse como aplicaciones dedicadas o en arquitecturas cliente-servidor, incluyendo una arquitectura basada en la web, y pueden incluir programas funcionales, códigos, y segmentos de código. Cualquiera de los ordenadores puede comprender un procesador, una memoria para almacenar los datos de programa y ejecutarlo, un almacenamiento permanente tal como una unidad de disco, un puerto de comunicaciones para manejar las comunicaciones con dispositivos externos y dispositivos de interfaz de usuario, incluyendo una pantalla, un teclado, un ratón, etc. Normalmente, el procedimiento se realiza usando módulos de software que pueden almacenarse como instrucciones de programa o códigos legibles por ordenador ejecutables en el procesador en un medio legible por ordenador tales como una memoria de solo lectura (ROM), una memoria de acceso aleatorio (RAM), unos CD-ROM, unas cintas magnéticas, unos disquetes, y unos dispositivos ópticos de almacenamiento de datos. El medio de registro legible por ordenador también puede distribuirse a través de una red acoplada a los sistemas informáticos de manera que el código legible por ordenador se almacena y se ejecuta de manera distribuida. Este medio es legible por el ordenador, se almacena en la memoria, y se ejecuta por el procesador.

REIVINDICACIONES

1. Un método de procesamiento de datos sísmicos, que comprende: - introducir unos datos sísmicos (DS,G[t]) que incluyen, para cada una de las mismas una pluralidad de tomas en las localizaciones de origen respectivas, unas trazas sísmicas registradas en una pluralidad de localizaciones de receptor; - aplicar un proceso de migración en profundidad a los datos sísmicos para obtener un primer conjunto de datos migrados que incluyen, para cada toma, los primeros valores migrados (MS[x, y, z]), asociados respectivamente con una pluralidad de posiciones del subsuelo; - calcular los datos modulados desplazados (D’S,G[t]) multiplicando los datos sísmicos en cada traza sísmica por un desplazamiento horizontal entre las localizaciones de origen y de receptor para dicha traza sísmica; - aplicar el proceso de migración en profundidad a los datos modulados desplazados para obtener un segundo conjunto de datos migrados que incluyen, para cada toma, unos segundos valores migrados (M'S[x, y, z]) asociados respectivamente con la pluralidad de posiciones de subsuelo; - estimar, para cada toma, los valores de desplazamiento (ĥS[x, y, z]) asociados respectivamente con al menos alguna de las posiciones de subsuelo, por un proceso de división aplicado a los conjuntos de datos migrados primero y segundo; y - estimar una captación de imágenes comunes en una posición horizontal, que comprende los valores migrados respectivos (Rx,y[z,h]) para los pares de parámetros incluyendo cada uno de los mismos un parámetro de profundidad y un parámetro de desplazamiento, en el que el valor migrado para un par de parámetros en la captación de imágenes comunes en dicha posición horizontal es un primer valor migrado del primer conjunto de datos migrados asociados con una posición de subsuelo determinada por dicha posición horizontal y el parámetro de profundidad de dicho par de parámetros para una toma de tal manera que el valor de desplazamiento estimado asociado con dicha posición de subsuelo es el parámetro de desplazamiento de dicho par de parámetros. 2. El método de acuerdo con la reivindicación 1, en el que el proceso de división usado para estimar los valores de desplazamiento (ĥS[x, y, z]) asociados con una posición de subsuelo comprende minimizar una función de coste definida por una variable de desplazamiento y los valores locales de los valores migrados primero y segundo (MS[x, y, z]), (M'S[x, y, z]) en una vecindad de dicha posición de subsuelo.

3. El método de acuerdo con la reivindicación 1 o 2, en el que el proceso de migración en profundidad es un proceso de migración en tiempo inverso (RTM). 4. El método de acuerdo con la reivindicación 1 o 2, en el que el proceso de migración en profundidad es un proceso de migración en profundidad de pre-apilamiento de ecuación de onda.

5. Un sistema para procesar los datos sísmicos, que comprende un recurso informático configurado para realizar un método de acuerdo con una cualquiera de las reivindicaciones anteriores. 6. Un producto de programa informático de un sistema para procesar unos datos sísmicos, que comprende instrucciones para realizar un método de acuerdo con una cualquiera de las reivindicaciones 1-4, cuando dicho producto de programa se ejecuta en una unidad de procesamiento informático del sistema para procesar los datos sísmicos.