MXPA06001757A - Metodo para atenuar los multiples de capas de agua. - Google Patents

Abstract

Un metodo para atenuar los multiples de capa de agua de un grupo de rastros de datos sismicos. El metodo incluye pronosticar los multiples de capa de agua del lado del receptor en el grupo de rastros de datos sismicos con el uso de un operador convolutivo derivado de un modelo de capa de agua, extraer, en forma adaptable, los multiples de capa de agua del lado del receptor del grupo de rastros de datos sismicos, pronosticar una pluralidad de multiples de capa de agua del lado de fuente con el uso del operador convolutivo derivado de un modelo de capa de agua y extraer, en forma adaptable, los multiples de capa de agua del lado del receptor y los multiples de capa de agua del lado de la fuente del grupo de rastros de datos sismicos para generar una pluralidad de primarias en el grupo de rastros de datos sismicos.

Description

MÉTODO PARA ATENUAR LOS MÚLTIPLES DE CAPAS DE AGUA CAMPO DE LA INVENCIÓN Las modalidades de la presente invención se relacionan en general con el campo de datos sísmicos y más en particular a la predicción y retiro de múltiples en la inspección de reflexión sísmica.

ANTECEDENTES DE LA INVENCIÓN La inspección sísmica es un método para determinar la estructura de formaciones subterráneas en la tierra. La inspección sísmica típicamente utiliza las fuentes de energía sísmica que generan ondas sísmicas y receptores sísmicos que están colocados estratégicamente para detectar las ondas sísmicas. Las ondas sísmicas se propagan en formaciones de tierra, en donde una porción de las ondas se refleja de las interfaces entre las formaciones subterráneas. La amplitud y polaridad de las ondas reflejadas se determinan por las diferencias en la ¡mpedancia acústica entre las capas de roca que comprenden las formaciones subterráneas. La impedancia acústica de una capa de rocas es el producto de la velocidad de propagación acústica dentro de la capa y la densidad de la capa. Las ondas sísmicas reflejadas se detectan por los receptores sísmicos, que convierten las ondas reflejadas en señales eléctricas representativas. Típicamente, las señales se transmiten por un medio eléctrico, óptico, de radio a dispositivos que graban las señales. A través del análisis de las señales grabadas, la forma, posición y composición se pueden determinar las formaciones subterráneas. La inspección sísmica terrestre es un método para determinar la estructura de las formaciones subterráneas debajo de la superficie terrestre. Típicamente, la fuente de energía sísmica utilizada en la inspección sísmica terrestre es un aparato con la capacidad de entregar una serie de impactos o vibraciones mecánicas a la superficie terrestre o la detonación de una carga explosiva cerca de la superficie terrestre, mientras el receptor sísmico utilizado es un sensor de movimiento, tal como un geófono o un acelerómetro. Las fuentes sísmicas y los receptores sísmicos típicamente se colocan en la superficie terrestre, aunque la fuente o el receptor se puede colocar en un orificio de pozo para determinar los perfiles verticales. Tanto las fuentes sísmicas, como los receptores sísmicos típicamente se vuelven a colocar para abarcar el área de inspección.

La inspección sísmica marina es un método para determinar la estructura de las formaciones subterráneas que subyacen a los cuerpos de agua. La inspección sísmica marina típicamente utiliza fuentes de energía sísmica y receptores sísmicos ubicados en el agua y se jalan detrás de una barca o se colocan en el fondo del agua desde la barca. La fuente de energía, típicamente, es un dispositivo explosivo o un sistema de aire comprimido que genera energía sísmica, los cuales están configurados para propagar como ondas sísmicas a través del cuerpo de agua y dentro de las formaciones de tierra por debajo del fondo del agua. Conforme las ondas sísmicas chocan con las interfaces entre las formaciones subterráneas, una porción de las ondas sísmicas se refleja de regreso a través de la tierra y el agua hacia los receptores sísmicos, para ser detectadas, transmitidas y grabadas. Los receptores sísmicos utilizados típicamente en la inspección sísmica marina son sensores de presión, como hidrófonos. Tanto las fuentes como los receptores se pueden recolocar para abarcar el área de inspección. Las ondas sísmicas, sin embargo, no se reflejan solamente de las interfaces entre las formaciones subterráneas, como sería deseable. Las ondas sísmicas también se reflejan del fondo del agua y de la superficie del agua, y las ondas reflejadas resultantes continúan reflejándose por sí mismas. Las ondas que se reflejan muchas veces se llaman "múltiples". Las ondas que se reflejan múltiples veces en la capa de agua entre la superficie del agua sobre y por debajo del fondo del agua se llaman "múltiples de fondo del agua" o "múltiples de capas de agua". Los múltiples de capas de agua se reconocen como un problema en el procesamiento e interpretación sísmica marina, y en consecuencia, se han desarrollado muchos métodos de atenuación de múltiples para manejar los múltiples de capas de agua. En general, los métodos actuales de atenuación de múltiples pueden dividirse en dos categorías: métodos de filtración y métodos de predicción de ecuación de onda. Los métodos de filtración se basan en la periodicidad de los múltiples o en las diferencias importantes de velocidad entre las reflexiones primarias ("primarias") y los múltiples. Las reflexiones primarias son aquellas ondas sísmicas que solamente se reflejan una vez, desde el fondo del agua o de una interfaz entre las formaciones subterráneas, antes de ser detectadas por el receptor sísmico. La desconvolución predictiva, es un método de filtración que supone que los múltiples son periódicos, mientras las primarias no. Esta suposición usualmente se cumple para los datos desde profundidades menores a 500 msec (aproximadamente 360 m) y aproximadamente una geología de subsuperficie en capas. En áreas en donde las profundidades de agua son mayores que 500 msec, en donde la diferencia de velocidad entre las primarias y los múltiples es importante, se pueden utilizar los métodos de filtración de velocidad tales como el tau-p y filtración f-k, en donde la variable f representa la frecuencia, k representa el número de onda, p representa el parámetro de rayo y tau representa el tiempo de intercepción del desplazamiento cero. Sin embargo, los métodos de filtración por lo general, requieren de una determinación o por lo menos un pronóstico educado, de las velocidades de propagación de onda en el medio sub-superficial a través del cual pasan las ondas sísmicas reflejadas en su trayecto desde la fuente sísmica a un receptor. Estas velocidades pueden diferir mucho de un tipo de medio a otro. Además, la desconvolución predictiva con frecuencia llega a un retiro accidental de primarias en el proceso de retirar los múltiples. Además, la desconvolución predictiva con frecuencia, falla en tomar en cuenta el factor no lineal en la reflexión, que por lo general es provocada por múltiples de jalón. Los métodos de ecuación de onda, por lo general, utilizan el fenómeno de propagación de onda física para pronosticar y extraer los múltiples de ios datos. Los métodos de ecuación de onda pueden ser muy exactos, pero también son costosos y consumidores de tiempo comparados con los métodos de filtración. Los métodos de ecuación de onda utilizan el hecho de que las primarias y los múltiples están relacionados físicamente. Los métodos de ecuación de onda, por lo general, pueden manejar geometrías complejas y no necesitan o necesitan poca información acerca de las propiedades de la sub-superficie. Sin embargo, algunos métodos de ecuación de onda, requieren información estructural, es decir, información acerca de la estructura de la sub-superficie, cuya determinación es la razón principal para hacer la exploración sísmica. Otros métodos de ecuación de onda requieren la forma de la forma de onda que no será una función pura de delta debido a las reverberaciones y la limitación de ancho de banda de frecuencia. Algunos métodos de ecuación de onda requieren la información de forma de onda fuente y la estructural. Cuando la forma de la forma onda fuente no se toma en cuenta, es decir, se ignora, los múltiples pronosticados no coincidirán con los múltiples reales y los intentos para retirar los múltiples serán fallidos. Por lo tanto, existe la necesidad en la técnica de un método para pronosticar con exactitud y remover los múltiples en un grupo de rastros de datos sísmicos.

BREVE DESCRIPCIÓN DE LA INVENCIÓN Las modalidades de la presente invención por lo general, están dirigidas a un método para atenuar los múltiples de capas de agua desde un grupo de rastros de datos sísmicos. El método incluye pronosticar una pluralidad de múltiples de capa de agua en el lado receptor en el grupo de rastros de datos sísmicos, con el uso de un operador convolutivo derivado del modelo de capa de agua, extraer, en forma adaptable los múltiples de capa de agua del lado del receptor del grupo de rastros de datos sísmicos, pronosticar una pluralidad de múltiples de capa de agua del lado de la fuente con el uso del operador convolutivo derivado del modelo de capa de agua y, extraer, en forma adaptable, los múltiples de capa de agua del lado del receptor y los múltiples de capa de agua del lado de la fuente del grupo de rastros de datos sísmicos para generar una pluralidad de primarias en el grupo de rastros de datos sísmicos. En una modalidad, el método incluye formar el grupo de rastros de datos sísmicos en un dominio t-x, transformar el grupo de rastros de datos sísmicos del dominio t-x en un dominio tau-p, convolutar el grupo de rastros de datos sísmicos con un operador convolutivo para pronosticar un primer grupo de múltiples de capa de agua en el grupo de rastros de datos sísmicos, extraer, en forma adaptable, el primer grupo de múltiples de capa de agua del grupo de rastros de datos sísmicos, remover la primaria de fondo de agua del grupo de rastros de datos sísmicos, convolutar el operador convolutivo con el grupo de rastros de datos sísmicos después de que se ha extraído en forma adaptable el primer grupo de múltiples de capa de agua de los rastros de datos sísmicos y después de que la primaria de fondo de agua ha sido retirada del grupo de rastros de datos sísmicos para pronosticar un segundo grupo de múltiples de capas de agua en el grupo de rastros de datos sísmicos, agregar el primer grupo de múltiples de capa de agua al segundo grupo de múltiples de capa de agua, transformar la suma del primer grupo de múltiples de capa de agua y el segundo grupo de múltiples de capa de agua del dominio tau-p al dominio t-x, y extraer, en forma adaptable, la suma transformada del primer grupo de múltiples de capa de agua y el segundo grupo de múltiples de capa de agua del grupo de rastros de datos sísmicos en el dominio t-x para generar una pluralidad de primarias en el grupo de rastros de datos sísmicos. Las modalidades de la presente invención también están dirigidas a un método para generar un operador convolutivo configurado para ser aplicado a un grupo de rastros de datos sísmicos. El método incluye generar un valor estimado de un tiempo de trayecto de dos vías de un desplazamiento cero en una capa de agua y un valor estimado de la reflexión del fondo del agua desde un modelo de capa de agua, y computar el operador convolutivo con el uso de los valores estimados del tiempo de trayecto de dos vías de desplazamiento cero y la reflexión del fondo del agua.

BREVE DESCRIPCIÓN DE LOS DIBUJOS La siguiente descripción detallada hace referencia a los dibujos acompañantes, los cuales se describen brevemente. La Figura 1 ilustra un diagrama esquemático de una inspección sísmica marina. La Figura 2 ilustra un método para atenuar los múltiples de capa de agua de un grupo de rastros de datos sísmicos de conformidad con una modalidad de la invención. La Figura 3 ¡lustra una formación de los rastros de datos sísmicos en el dominio t-x de conformidad con una modalidad de la invención. La Figura 4 ilustra una formación de rastros de datos sísmicos en un dominio tau-p de conformidad con una modalidad de la invención. La Figura 5 ilustra un método para derivar al operador convolutivo o de reflexión de conformidad con una modalidad de la invención. Mientras esta invención se describe como ejemplificativa para varias modalidades y dibujos ilustrativos, las personas experimentadas en la técnica reconocerán que la invención no está limitada a las modalidades o dibujos descritos. Se debe entender que los dibujos y la descripción detallada no tienen la intención de limitar la invención a la forma particular descrita, por el contrario, la intención es abarcar todas las modificaciones, equivalentes y alternativas que caigan dentro del espíritu y alcance de la presente invención, como se define por las reivindicaciones anexas. Los subtítulos utilizados tienen el propósito de organizar y no tienen la intención de limitar el alcance de la descripción o de las reivindicaciones. A través de esta especificación, la palabra "puede" se utiliza en un sentido permisible (es decir, que tiene el potencial para.\.), mejor que en un sentido obligatorio (es decir, con significado debe...). De forma similar, las palabras "incluir", Incluyen" e "incluye" significan incluyendo, pero no limitando.

DESCRIPCIÓN DETALLADA DE LA INVENCIÓN La Figura 1 muestra un diagrama esquemático de una inspección 100 sísmica marina. Se exploran las formaciones subterráneas, como la 102 y 104, por debajo de un cuerpo de agua 106. Las fuentes 108 de energía sísmica y los receptores 110 sísmicos están colocados en el cuerpo de agua 106, típicamente por una o más recipientes sísmicos (no mostrados). Una fuente 108 sísmica, tal como una pistola de aire, crea ondas sísmicas en el cuerpo de agua 106 y una porción de las ondas sísmicas viaja hacia abajo a través del cuerpo de agua 106 hacia las formaciones subterráneas 102 y 104 debajo del cuerpo de agua 106. Cuando las ondas sísmicas alcanzan un reflector sísmico, una porción de las ondas sísmicas se refleja hacia arriba y una porción de las ondas sísmicas continúa hacia abajo. El reflector sísmico puede estar en el fondo 112 de agua o en una de las interfaces entre la formación subterránea, como la interfaz 114 entre las formaciones 102 y 104. Cuando las ondas reflejadas que viajan hacia arriba alcanzan la interfaz de agua/aire en la superficie 116 del agua, una porción mayor de las ondas se reflejan hacia abajo. Para continuar en esta manera, las ondas sísmicas pueden reflejar múltiples veces entre el fondo 112 de agua o las interfaces de formación por debajo y la superficie 116 del agua. Cada vez que las ondas reflejadas se propagan más allá de la posición del receptor 110 sísmico, el receptor 110 detecta las ondas reflejadas y genera señales representativas. Las reflexiones primarias son aquellas ondas sísmicas que se han reflejado una vez desde el fondo 112 del agua o de una ¡nterfaz entre las formaciones subterráneas, antes de ser detectadas por el receptor 110 sísmico. Las reflexiones primarias contienen la información deseada acerca de las formaciones subterráneas, la cual busca la inspección sísmica marina. Un ejemplo de la reflexión primaria se muestra en la Figura 1 por trayectos de rayos 120 y 122. Las ondas sísmicas se generan por una fuente 108 sísmica y una porción de las ondas viaja hacia abajo a través del cuerpo de agua 106 y dentro de la formación 102 subterránea a lo largo del trayecto de rayo 120. Una porción de las ondas sísmicas se refleja desde la interfaz 114 entre formaciones 102 y 104. Las ondas reflejadas viajan hacia arriba a través de la formación 102 y de regreso al cuerpo de agua 106 a lo largo del trayecto de rayo 122. Las ondas reflejadas viajan más allá del receptor 110, el cual detecta las ondas y genera una señal representativa. Los múltiples de capa de agua son aquellas ondas que han sido reflejadas múltiples veces entre la superficie 116 del agua y el fondo 112 del agua antes de ser detectadas por un receptor 110. Un ejemplo de los múltiples de capa de agua se muestra al inicio del trayecto de rayo 130. Las ondas sísmicas se generan por una fuente 108 sísmica y una porción de las ondas viaja hacia abajo a través del cuerpo de agua 106 a lo largo del trayecto de rayo 130. Una porción de las ondas sísmicas se refleja desde el fondo del agua 112 y viaja de regreso hacia arriba a través del cuerpo de agua 106 a lo largo del trayecto de rayo 132. Una porción de las ondas reflejadas se refleja desde la superficie 116 del agua y viaja de regreso hacia abajo a través del cuerpo de agua 106 a lo largo del trayecto de rayo 134. Una porción de las ondas reflejadas dos veces se refleja otra vez desde el fondo 112 del agua y viaja de regreso hacia arriba a través del cuerpo de agua 106 a lo largo del trayecto de rayo 136. Una porción de las ondas reflejadas tres veces se refleja otra vez desde la superficie 116 del agua y viaja de regreso hacia abajo a través del cuerpo de agua 106 a lo largo del trayecto de rayo 138. Una porción de las ondas reflejadas cuatro veces se refleja otra vez desde el fondo 112 del agua y viaja de regreso hacia arriba a través del cuerpo de agua 106 a lo largo del trayecto de rayo 140. Una porción de las ondas reflejadas cinco veces se refleja otra vez desde la superficie 116 del agua y viaja de regreso hacia abajo a través del cuerpo de agua 106 a lo largo del trayecto de rayo 142. Una porción de las ondas reflejadas seis veces se refleja otra vez desde el fondo 112 del agua y viaja de regreso hacia arriba a través del cuerpo de agua 106 a lo largo del trayecto de rayo 144. Las ondas reflejadas siete veces viajan más allá del receptor 110, el cual detecta las ondas y genera una señal representativa. Los múltiples de capa de agua que inician en el trayecto de rayo 130 es un múltiple de capa de agua sencillo ya que las ondas sísmicas nunca van por debajo del fondo 112 del agua, es decir, se pueden quedar dentro del cuerpo de agua o en el cuerpo de agua 106. El múltiple de capa de agua que empieza en el trayecto de rayo 130 es un múltiple del orden tres, ya que los múltiples contienen tres reflexiones desde la fuente 116 del agua. Los múltiples de capa de agua que van por debajo del fondo 112 del agua con frecuencia son llamados como múltiples de capa de agua de jalón. La Figura 1 ilustra un ejemplo de un múltiples de capa de agua de jalón que empieza en el trayecto de rayo 150. Las ondas sísmicas se generan por una fuente 108 sísmica y una porción de las ondas viaja hacia abajo a través del cuerpo de agua 106 y dentro de la formación 102 subterránea a lo largo del trayecto de rayo 150. Una porción de las ondas sísmicas se refleja desde el fondo 112 del agua y viaja de regreso hacia arriba a través del cuerpo de agua 106 a lo largo del trayecto de rayo 152. Una porción de las ondas reflejadas viaja desde la superficie 116 del agua y viaja de regreso hacia abajo a través del cuerpo de agua 106 y dentro de la formación 102 subterránea a lo largo del trayecto de rayo 154.

Una porción de las ondas reflejadas dos veces viaja desde la interfaz 114 y viaja de regreso hacia arriba a través de la formación 102 y el cuerpo de agua 106 a lo largo del trayecto de rayo 156. Una porción de las ondas reflejadas tres veces se refleja otra vez desde la superficie 116 del agua y viaja de regreso hacia abajo a través del cuerpo de agua 106 a lo largo del trayecto de rayo 158. Una porción de las ondas reflejadas cuatro veces se refleja otra vez desde el fondo 112 del agua y viaja de regreso hacia arriba a través del cuerpo de agua 106 a lo largo del trayecto de rayo 160. Una porción de las ondas reflejadas cinco veces se refleja otra vez desde la superficie 116 del agua y viaja de regreso hacia abajo a través del cuerpo de agua 106 a lo largo del trayecto de rayo 162. Una porción de las ondas reflejadas seis veces se refleja otra vez desde el fondo 112 del agua y viaja de regreso hacia arriba a través del cuerpo de agua 106 a lo largo del trayecto de rayo 164. Las ondas reflejadas siete veces entonces viajan más allá del receptor 110, el cual detecta las ondas y genera una señal representativa. Los múltiples de capa de agua de jalón que inician en el trayecto de rayo 150 es un múltiple de orden tres, ya que el múltiple contiene tres reflexiones desde la superficie 116 del agua. Una de las tres reflexiones desde la superficie 116 del agua se ubican en la fuente 108 lateral, mientras las otras dos están ubicadas en el receptor 110 lateral. La Figura 2 ilustra un método 200 para atenuar los múltiples de capa de agua desde un grupo de rastros D de datos sísmicos de conformidad con una modalidad de la invención. Los rastros D de datos pueden ser clasificados en grupos, como grupos de disparo, grupos de punto medio común, o grupos receptores. Un grupo de disparo se define como un sub-grupo de rastros de datos que tienen una ubicación de disparo (fuente) común. El grupo de rastros D de datos sísmicos contiene una reflexión primaria del fondo del agua, un grupo completo de reflexiones primarias de la sub-superficie y todas las reflexiones de los múltiples de capa de agua (y otras) asociadas. Como tal, el grupo de rastros D de datos se puede representar como: D = W/(1+R) + P/(1+R)2 (1) en donde P representa las reflexiones primarias (excepto la primaria del fondo del agua) y todos los múltiples de capa de no agua en el grupo de rastros D de datos sísmicos, R representa un operador convolutivo o de reflexión, el cual consiste de un factor de desplazamiento de tiempo que es igual al tiempo de trayecto de dos vías en el cuerpo de agua o capa 106 de agua, y un factor de amplitud, que es igual a la reflexión del fondo del agua; W representa la reflexión primaria del fondo del agua en el grupo de rastros D de datos sísmicos ; P/(1 + R)2 representa las reflexiones primarias (excepto la primaria del fondo del agua) y todos los múltiples asociados, incluyendo los múltiples de capa de agua del lado del receptor y del lado de la fuente; (1+R)2 por lo general, se reconoce como el operador Backus y W/(1+R) representa la reflexión primaria del fondo del agua y todos los múltiples de capa de agua asociados. La ecuación (1) se puede aplicar en el dominio t-x en el desplazamiento cero, o en el dominio tau-p para cada valor del parámetro p de rayo. Cuando se aplica en el dominio tau-p, sin embargo, el factor de desplazamiento de tiempo asociado con R es el tau de tiempo de intercepción asociado con la reflexión W primaria del fondo del agua. La ecuación (1) puede re-escribirse como: W + P = D - M1 - M2 (2) en donde M1 = RD y M2 = R (D - M1- W). M1 representa un sub-grupo de múltiples de capa de agua equivalentes a múltiples del lado del receptor, tal como (D-M1-W) representa los rastros de datos grabados con los múltiples M1 del lado del receptor y la reflexión W primaria del fondo del agua retirada. M2 representa los múltiples restantes, es decir los múltiples del lado de fuente. W + P representa los datos primarios, es decir, eventos que no son múltiples de capa de agua, que es el objetivo de la presente invención. Un evento, por lo general, se define como fluctuaciones de presión grabadas que están asociadas con una ondas sísmica que pasa el receptor 110. Típicamente, estas fluctuaciones de presión son de corta duración (décimas de milísegundos) comparadas con el tiempo que se toma para que las ondas sísmicas viajen desde la fuente 108 al receptor 110 (centésimas de milísegundos para varios segundos). De este modo, las fluctuaciones de presión por lo general, aparecen como grabaciones distintas en los rastros D de datos sísmicos. El tiempo del evento normalmente será del tiempo de la generación de las fluctuaciones de presión. Un evento puede ser una primaria o un múltiple.

Con referencia otra vez a la Figura 2, en el paso 210, el grupo de rastros D de datos sísmicos formado en el dominio t-x, en donde t representa la cantidad del tiempo de trayecto de las ondas sísmicas desde la fuente 108 al receptor 110 y x representa la distancia de desplazamiento, es decir, la distancia horizontal, entre la fuente 108 y el receptor 110. En la Figura 3, se ilustra una formación 300 de los rastros D de datos sísmicos en el dominio t-x. El eje vertical representa la cantidad de tiempo de trayecto (t) como se describe antes, con la cantidad del tiempo de trayecto que aumenta en la dirección hacia abajo. El eje horizontal representa la distancia (x) de desplazamiento entre la fuente 108 y el receptor 110. Las curvas en la formación 300 indican los tiempos a los cuales se graban las reflexiones para cada evento como una función de x. La curva 310 representa la reflexión primaria del fondo del agua. Las curvas 311 y 312 representan los múltiples de capa de agua sencillos de primer y segundo orden. La curva 320 representa los múltiples de capa de agua sencillos de tercer orden mostrados en la Figura 1 como trayecto de rayos 130-144. La curva 315 representa la reflexión primaria de la ¡nterfaz 114, mostrada en la Figura 1 como trayecto de rayos 120 y 122. La curva 330 representa el múltiples de capa de agua de jalón mostrado en la Figura 1 como trayecto de rayos 150-164. La curva 321 representa el múltiple de capa de agua de jalón de primer orden y la curva 322 representa el múltiples de capa de agua de jalón de segundo orden.

Con referencia otra vez a la Figura 2, en el paso 220, los rastros D de datos sísmicos se transforman en un dominio tau-p por la transformación lineal Radon. De esta manera, el análisis con el uso de la ecuación (1) se puede aplicar en los datos de desplazamiento cero. Una formación 400 de los rastros D de datos en el dominio tau-p se ilustra en la Figura 4. El eje vertical representa tau, que es el tiempo en el cual la línea recta definida por tau y p intercepta el eje de desplazamiento cero en el dominio t-x y es el mismo que el tiempo de desplazamiento cero del evento cuando p =0. El eje horizontal representa el parámetro p de rayo, que se relaciona con un ángulo entre el trayecto de rayo y la dirección vertical. Las curvas 410, 411, 412, 415, 420, 421, 422 y 430 en el dominio tau-p corresponden a las curvas 310, 311, 312, 315, 320, 321, 322 y 330 en el dominio t-x, respectivamente. En el paso 230, el operador convolutivo R o de reflexión se aplica en el grupo de rastros D de datos sísmicos para generar un valor de M1 estimado, el cual representa los múltiples de capa de agua del lado del receptor en el grupo de rastros D de datos sísmicos. Esto es, el grupo de rastros D de datos sísmicos se convoiutan con el operador convolutivo R para pronosticar un valor de los múltiples M1 de capa de agua del lado del receptor. De esta forma, el valor de los múltiples M1 de capa de agua del lado del receptor se determinan en forma pronosticada. El operador convolutivo R puede derivarse de un método 500 descrito en la Figura 5.

El siguiente párrafo define el método 500 para derivar el operador convolutivo R o de reflexión de conformidad con una modalidad de la presente invención. En el paso 510, el tiempo de trayecto de dos vías del desplazamiento cero en la capa de agua y la reflexión en el fondo 112 del agua se calculan a partir de un modelo de capa de agua. El tiempo de trayecto de desplazamiento cero en la capa de agua también puede capturarse de las auto-correlaciones de los rastros de datos de desplazamiento cercado, corregidos de movimiento. La reflexión en el fondo 112 del agua puede ser escalable, que representará el fondo rígido del agua, definido como una transición suave del agua a un material razonablemente homogéneo. La reflexión en el fondo 112 del agua también puede estar con base en una serie de tiempos. Una serie de tiempos por lo general, se utiliza para la reflexión del fondo del agua cuando el fondo del agua es de transición, es decir, que consiste de muchas capas de sedimentos. La reflexión en el fondo del agua puede ser dependiente de p. Para cada p, se calculan la reflexión del fondo del agua y el tiempo de trayecto de desplazamiento cero en la capa de agua. También, en el paso 520, el operador convolutivo R o de reflexión se computa con base en los valores estimados del tiempo de trayecto de desplazamiento cero en la capa del agua y la reflexión en el fondo 112 del agua. El operador convolutivo R por lo general, se computa en el dominio tau-p. En una modalidad, el operador convolutivo R es la reflexión del fondo del agua desplazada en tiempo por el tiempo asociado con la capa de agua. En el dominio tau-p, el tiempo asociado con la capa de agua depende en el parámetro p de rayo y es equivalente al tiempo de evento W de la Figura 4. Este tiempo se puede computar como una función elíptica de] tiempo de trayecto de desplazamiento cero y el parámetro p de rayo, determinado como un cálculo de la velocidad del agua. Por ejemplo, cuando la reflexión del fondo del agua es 0.5 y el tiempo de trayecto asociado con la capa de agua para una p determinada es 200 ms, entonces el operador convolutivo R consistirá de una cresta de amplitud 0.5 en el tiempo 200 ms. Ya que el tiempo del agua y los cálculos de reflexión se derivan de un modelo, opuesto a la derivación estadística de los rastros D de datos sísmicos, por sí mismos, la predicción de los múltiples con base en el operador convolutivo R se puede describir como determinativa. Con referencia otra vez a la Figura 2, en el paso 240, el valor pronosticado de los múltiples M1 de capa de agua del lado del receptor se extraen en forma adaptable de los rastros D de datos sísmicos. Los múltiples M1 de capa de agua del lado del receptor pueden extraerse en forma adaptable de los rastros D de datos sísmicos con el uso de medias cuadradas o de un algoritmo de filtración adaptable de mínimas cuadradas. La extracción adaptable se puede configurar para permitir errores en los valores calculados del tiempo de trayecto de desplazamiento cero y la reflexión en el fondo 112 del agua. En una modalidad, ya que los múltiples M1 de capa de agua en el lado del receptor y los múltiples M2 de capa de agua del lado de la fuente por lo general, contienen los mismos eventos pero con diferentes amplitudes, tanto los múltiples M1 de capa de agua del lado del receptor y los múltiples M2 de capa de agua del lado de la fuente se pueden extraer en forma adaptable de los rastros D de datos con el uso de un algoritmo de filtro adaptable más agresivo. En tal caso, el procesamiento continúa en el paso 280 con los múltiples del lado del receptor pronosticados que reemplazan a los múltiples sumados. Los siguientes pasos se dirigen a la dependencia no lineal de los múltiples en la reflexión. En el paso 250, la reflexión W primaria del fondo del agua se retira o se elimina del grupo de rastros D de datos sísmicos. Al remover la reflexión W primaria del fondo del agua, cada amplitud asociada con la primaria del fondo del agua se reemplaza con cero. En una modalidad, este paso se puede alcanzar al eliminar una zona de rastros de datos sísmicos en el dominio tau-p que contienen la reflexión W primaria del fondo del agua. La reflexión W primaria del fondo del agua puede también removerse con otros métodos, como crear un modelo de la reflexión W primaria del fondo del agua y extraerla del grupo de rastros D de datos sísmicos, u otros métodos de filtración conocidos por las personas experimentadas en la técnica. En el paso 260, el operador convolutivo R o de reflexión (el cual se deriva del método 500 antes descrito) se aplica o se convoluta con los rastros de datos sísmicos después de que los múltiples M1 de capa de agua del lado del receptor han sido extraídos en forma adaptable de los rastros D de datos sísmicos y la primaria del fondo del agua ha sido removida de los rastros D de datos sísmicos. Esta convolución genera un valor calculado de M2, el cual representa los múltiples de capa de agua del lado de la fuente. De esta forma, el valor de los múltiples M2 de capa de agua del lado de la fuente se pronostican en forma determinativa. En el paso 270, los múltiples M1 de capa de agua del lado del receptor se agregan a los múltiples M2 de capa de agua del lado de la fuente para generar todos los múltiples M de capa de agua en los rastros D de datos grabados. En el paso 280, todos los múltiples M de capa de agua se transforman en el dominio tau-p para el dominio t-x por una transformación Radon inversa. En el paso 290, todos los múltiples M de capa de agua se extraen en forma adaptable de los rastros D de datos sísmicos para generar las reflexiones P primarias en los rastros de datos sísmicos. En otra modalidad, todos los múltiples M de capa de agua pueden permanecer en el dominio tau-p y extraerse en forma adaptable de los rastros D de datos sísmicos, seguidos por la transformación del resultado de la extracción adaptable al dominio t-x. Mientras lo anterior está dirigido a las modalidades de la presente invención, estas y otras modalidades de la invención pueden ser contempladas sin apartarse del alcance básico de la misma, el alcance es determinado por las reivindicaciones que siguen.

Claims (20)

REIVINDICACIONES

1. Un método para atenuar múltiples de capa de agua de un grupo de rastros de datos sísmicos, caracterizado porque comprende: pronosticar una pluralidad de múltiples de capa de agua del lado del receptor en el grupo de rastros de datos sísmicos con el uso de un operador convolutivo derivado de un modelo de capa de agua; extraer, en forma adaptable, los múltiples de capa de agua del lado del receptor del grupo de rastros de datos sísmicos; pronosticar una pluralidad de múltiples de capa de agua del lado de fuente con el uso de un operador convolutivo derivado del modelo de capa de agua; y extraer, en forma adaptable, los múltiples de capa de agua del lado del receptor y los múltiples de capa de agua del lado de la fuente del grupo de rastros de datos sísmicos para generar una pluralidad de primarias en el grupo de rastros de datos sísmicos.

2. El método de conformidad con la reivindicación 1, caracterizado porque pronosticar la pluralidad de múltiples de capa de agua del lado del receptor comprende: formar el grupo de datos sísmicos en el dominio t-x; transformar el grupo de rastros de datos sísmicos del dominio t-x al dominio tau-p; y convolutar el grupo de rastros de datos sísmicos con un operador convolutivo para pronosticar los múltiples de capa de agua del lado del receptor.

3. El método de conformidad con la reivindicación 2, caracterizado porque pronosticar los múltiples de capa de agua del lado de la fuente comprende: remover la primaria del fondo del agua del grupo de rastros de datos sísmicos; y convolutar el operador convolutivo con el grupo de rastros de datos sísmicos después de que los múltiples de capa de agua del lado del receptor han sido extraídos en forma adaptable del grupo de rastros de datos sísmicos y después de que la primaria del fondo del agua ha sido removida del grupo de rastros de datos sísmicos para pronosticar los múltiples de capa de agua del lado de la fuente.

4. El método de conformidad con la reivindicación 2, caracterizado porque el grupo de rastros de datos sísmicos se transforma al dominio tau-p con el uso de una transformación Radon lineal.

5. El método de conformidad con la reivindicación 3, caracterizado porque además comprende: agregar los múltiples de capa de agua del lado del receptor a los múltiples de capa de agua del lado de la fuente; transformar la suma de los múltiples de capa de agua del lado del receptor y los múltiples de capa de agua del lado de la fuente del dominio tau-p al dominio t-x.

6. El método de conformidad con la reivindicación 5, caracterizado porque la extracción adaptable de los múltiples de capa de agua del lado del receptor y los múltiples de capa de agua del lado de la fuente del grupo de rastros de datos sísmicos comprende extraer en forma adaptable, la suma de los múltiples de capa de agua del lado del receptor y los múltiples de capa de agua del lado de la fuente en el dominio t-x para generar las primarias en el grupo de rastros de datos sísmicos en el dominio t-x.

7. El método de conformidad con la reivindicación 3, caracterizado porque remover la primaria del fondo del agua comprende reemplazar cada amplitud asociada con la primaria del fondo del agua con cero.

8. El método de conformidad con la reivindicación 3, caracterizado porque la suma de los múltiples de capa de agua del lado del receptor y los múltiples de capa de agua del lado de la fuente se transforma en el dominio t-x con el uso de una transformación lineal inversa Radon.

9. El método de conformidad con la reivindicación 1, caracterizado porque el operador convolutivo se computa con el uso de un tiempo de trayecto de dos vías de desplazamiento cero en una capa de agua y una reflexión en el fondo de agua calculado del modelo de capa de agua.

10. El método de conformidad con la reivindicación 9, caracterizado porque el operador convolutivo es el valor estimado de la reflexión del fondo del agua desplazada en tiempo por el valor estimado del tiempo de trayecto en la capa de agua en el dominio tau-p.

11. Un método para atenuar los múltiples de capa de agua de un grupo de rastros de datos sísmicos, caracterizado porque comprende: formar el grupo de rastros de datos sísmicos en el dominio t-x; transformar el grupo de rastros de datos sísmicos del dominio t-x al dominio tau-p; convolutar el grupo de rastros de datos sísmicos con un operador convolutivo para pronosticar un primer grupo de múltiples de capa de agua en el grupo de rastros de datos sísmicos; extraer, en forma adaptable, el primer grupo de múltiples de capa de agua del grupo de rastros de datos sísmicos; remover la primaria de fondo de agua del grupo de rastros de datos sísmicos; convolutar el operador convolutivo con el grupo de rastros de datos sísmicos después de que el primer grupo de múltiples de capa de agua ha sido extraído en forma adaptable de los rastros de datos sísmicos y después de que la primaria del fondo del agua ha sido removida del grupo de rastros de datos sísmicos para pronosticar un segundo grupo de múltiples de capa de agua en el grupo de rastros de datos sísmicos; agregar el primer grupo de múltiples de capa de agua al segundo grupo de múltiples de capa de agua; transformar la suma del primer grupo de múltiples de capa de agua y el segundo grupo de múltiples de capa de agua del dominio tau-p al dominio t-x; y extraer, en forma adaptable, la suma transformada del primer grupo de múltiples de capa de agua y el segundo grupo de múltiples de capa de agua del grupo de rastros de datos sísmicos en el dominio t-x para generar una pluralidad de primarias en el grupo de rastros de datos sísmicos.

12. El método de conformidad con la reivindicación 1, caracterizado porque el primer grupo de múltiples de capa de agua comprende uno o más múltiples de capa de agua del lado del receptor.

13. El método de conformidad con la reivindicación 1, caracterizado porque el segundo grupo de múltiples de capa de agua comprende uno o más múltiples de capa de agua del lado de la fuente.

14. El método de conformidad con la reivindicación 1, caracterizado porque remover la primaria del fondo del agua comprende reemplazar cada amplitud asociada con la primaria del fondo del agua. con cero.

15. El método de conformidad con la reivindicación 1, caracterizado porque la pluralidad de rastros de datos sísmicos se transforma al dominio tau-p con el uso de una transformación lineal Radon.

16. El método de conformidad con la reivindicación 1, caracterizado porque el operador convolutivo se deriva de un modelo de capa de agua.

17. El método de conformidad con la reivindicación 1, caracterizado porque el operador convolutivo se deriva de un tiempo de trayecto de dos vías de desplazamiento cero en la capa de agua y la reflexión en el fondo del agua calculados de un modelo de capa de agua.

18. Un método para generar un operador convolutivo configurado para ser aplicado a un grupo de rastros de datos sísmicos; caracterizado porque comprende: generar un valor estimado de un tiempo de trayecto de dos vías de desplazamiento cero en una capa de agua y un valor estimado de una reflexión del fondo del agua de un modelo de capa de agua; y computar un operador convolutivo con el uso de los valores estimados de tiempo de trayecto de dos vías de desplazamiento cero y la reflexión del fondo del agua.

19. El método de conformidad con la reivindicación 18, caracterizado porque el operador convolutivo se configura para convolutar con el grupo de rastros de datos sísmicos para pronosticar un grupo de múltiples en el grupo de rastros de datos sísmicos.

20. El método de conformidad con la reivindicación 18, caracterizado porque computar el operador convolutivo comprende desplazar el valor estimado de la reflexión del fondo del agua en tiempo por el valor estimado del tiempo de trayecto en la capa de agua en el dominio tau-p.