ES2671789T3 - Búsqueda simultánea de fuentes durante la adquisición sísmica y la inversión sísmica - Google Patents

Búsqueda simultánea de fuentes durante la adquisición sísmica y la inversión sísmica Download PDF

Info

Publication number
ES2671789T3
ES2671789T3 ES14744730.4T ES14744730T ES2671789T3 ES 2671789 T3 ES2671789 T3 ES 2671789T3 ES 14744730 T ES14744730 T ES 14744730T ES 2671789 T3 ES2671789 T3 ES 2671789T3
Authority
ES
Spain
Prior art keywords
source
sources
record
simultaneous
measurement
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
ES14744730.4T
Other languages
English (en)
Inventor
Christine E. Krohn
Partha S. Routh
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
ExxonMobil Upstream Research Co
Original Assignee
ExxonMobil Upstream Research Co
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by ExxonMobil Upstream Research Co filed Critical ExxonMobil Upstream Research Co
Application granted granted Critical
Publication of ES2671789T3 publication Critical patent/ES2671789T3/es
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/282Application of seismic models, synthetic seismograms
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/003Seismic data acquisition in general, e.g. survey design
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/003Seismic data acquisition in general, e.g. survey design
    • G01V1/005Seismic data acquisition in general, e.g. survey design with exploration systems emitting special signals, e.g. frequency swept signals, pulse sequences or slip sweep arrangements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/24Recording seismic data
    • G01V1/247Digital recording of seismic data, e.g. in acquisition units or nodes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/303Analysis for determining velocity profiles or travel times
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362Effecting static or dynamic corrections; Stacking
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/364Seismic filtering
    • G01V1/368Inverse filtering
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
    • G01V1/3808Seismic data acquisition, e.g. survey design
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/10Aspects of acoustic signal generation or detection
    • G01V2210/12Signal generation
    • G01V2210/121Active source
    • G01V2210/1212Shot
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/10Aspects of acoustic signal generation or detection
    • G01V2210/12Signal generation
    • G01V2210/127Cooperating multiple sources
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/10Aspects of acoustic signal generation or detection
    • G01V2210/14Signal detection
    • G01V2210/142Receiver location
    • G01V2210/1423Sea
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/30Noise handling
    • G01V2210/32Noise reduction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616Data from specific type of measurement
    • G01V2210/6161Seismic or acoustic, e.g. land or sea measurements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/622Velocity, density or impedance
    • G01V2210/6222Velocity; travel time
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/66Subsurface modeling

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Acoustics & Sound (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Oceanography (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)

Abstract

Un procedimiento para deducir un modelo de subsuperficie de velocidad u otra propiedad física, que comprende: obtener (702) datos sísmicos adquiridos en una medición haciendo avanzar una fuente anterior y una fuente posterior de seguimiento a lo largo de una línea de fuente con un tendido no fijo de receptores, y activar cada fuente a intervalos seleccionados a lo largo de la línea de fuente, en el que cada ubicación de activación para la fuente posterior se sitúa sobre una ubicación de activación previa para la fuente anterior, y en el que se utiliza una codificación de campo seleccionada para todos los tiros de fuente siendo la codificación de campo para la fuente posterior siempre la misma que la aplicada anteriormente a la fuente anterior en la misma ubicación de activación, dando como resultado una serie de registros de tiro, registrando cada uno simultáneamente ambas fuentes pero con cada tiro registrado por menos de un tendido completo de receptores; y procesar los registros de tiro por etapas que comprenden: (i) construir (703) a partir de los registros de tiro una pluralidad de pseudorregistros de súpertiro, construidos de modo que cada uno tiene datos de un tendido completo de receptores, en el que cada pseudorregistro de súpertiro incluye una pluralidad de ubicaciones de fuente; (ii) codificar (705) cada pseudorregistro de súpertiro y acumularlo para formar un registro de fuentes simultáneas de datos medidos; (iii) utilizar un ordenador para simular (706) el registro de fuentes simultáneas de datos medidos, utilizando la misma codificación utilizada en (ii) y también la codificación de campo, y utilizando un modelo de subsuperficie supuesto de velocidad u otra propiedad física; y (iv) comparar (707) el registro de fuentes simultáneas simulado con el registro de fuentes simultáneas de datos medidos, y determinar a partir de aquí un ajuste al modelo de subsuperficie de velocidad u otra propiedad física.

Description

5
10
15
20
25
30
35
40
45
50
55
60
65
DESCRIPCION
Búsqueda simultánea de fuentes durante la adquisición sísmica y la inversión sísmica Campo de la invención
Esta divulgación se refiere en general al campo de la prospección geofísica y, más particularmente, al procesamiento de datos sísmicos tal como se utiliza en la exploración de hidrocarburos. Específicamente, la divulgación se refiere a un procedimiento para adquirir datos, a un costo de adquisición sísmica reducido, utilizando fuentes simultáneas en el campo, y a continuación construyendo pseudorregistros de fuente que cumplan mejor con los requisitos para utilizar una búsqueda simultánea adicional de fuentes para simulaciones por ordenador o modelado directo como parte de procedimientos de inversión iterativa que actualizan un modelo de subsuperficie, tales como FWI (Full Wavefield Inversión, inversión de forma de onda completa) o LSRTM (Least-Square Reverse Time Migration, migración de tiempo inverso de mínimos cuadrados), con una reducción adicional en los costos computacionales.
Antecedentes de la invención
Adquisición y luego separación
La búsqueda simultánea de fuentes, también denominada búsqueda de fuentes combinada, es un procedimiento de adquisición sísmica emergente para reducir los costos de adquisición y mejorar el muestreo espacial. Convencionalmente, se adquieren mediciones localizando una única fuente puntual o una disposición de fuentes puntuales en una única ubicación de fuente, detonando las fuentes al mismo tiempo y a continuación registrando la respuesta durante el tiempo necesario para que las fuentes dejen de detonar tras lo cual se produce un tiempo de escucha en el que se registra todo lo que devuelve el objetivo de la subsuperficie. Opcionalmente, es posible repetir la detonación de las fuentes y pueden registrarse múltiples registros en la misma ubicación. Entonces se mueve la disposición de fuente a otra ubicación, y se repite el proceso. El costo de la adquisición de datos sísmicos con este procedimiento secuencial está relacionado con el tiempo necesario para registrar cada ubicación de fuente individual y el número de tales ubicaciones, y a menudo este costo limita la capacidad para registrar datos con un muestreo fino. Mediante la detonación de una o varias fuentes puntuales en diferentes ubicaciones de fuente al mismo tiempo o casi al mismo tiempo dentro del mismo registro de datos, pueden reducirse el tiempo y costo de adquisición y aumentarse el muestreo. Esto puede denominarse adquisición simultánea. Originalmente, cuando se introdujo el procedimiento, las fuentes en cuestión se excitaban exactamente al mismo tiempo o simultáneamente. En la actualidad, se utiliza el mismo término también para una adquisición en la que las fuentes detonan dentro de la misma ventana de tiempo que otra fuente aunque la detonación de las fuentes no sea simultánea en el tiempo y difiera con cierto retardo de tiempo. En general, las fuentes que detonan casi al mismo tiempo dentro del mismo registro breve forman una disposición espacial o de área extendida, sin ninguna expectativa de que las posiciones de las fuentes puntuales individuales estén próximas entre sí. La desventaja de la adquisición simultánea es la necesidad de mitigar la energía solapante o interferencia entre las fuentes en diferentes ubicaciones mediante una combinación de codificación de fuente en el campo y mediante técnicas de filtrado y separación de fuentes en el procesamiento. El procesamiento convencional requiere registros individuales para cada ubicación de fuente y éstos deben extraerse o separarse de los registros de datos registrados.
La búsqueda simultánea de fuentes se utiliza más comúnmente para fuentes Vibroseis con funciones de barrido largo, que pueden codificarse fácilmente. Con el procedimiento Vibroseis, cada vibrador individual puede accionarse con un barrido que difiere en cierto modo de los barridos para otros vibradores dentro de la disposición, por ejemplo utilizando las diferencias en la fase de barrido, función pseudoaleatoria, velocidad de barrido, intervalo de frecuencia de barrido, tiempo de inicio, etc. Algunos procedimientos requieren múltiples barridos y múltiples registros por ubicación para la separación. En el caso especial de que el número de barridos sea mayor que o igual al número de vibradores, entonces los registros de fuentes individuales pueden extraerse casi perfectamente de los múltiples registros combinados aplicando un filtro inverso tal como se describe para el procedimiento HFVS en Sallas, et al. (patente estadounidense n.° 5.721.710). Con este procedimiento y procedimientos similares, es crítico que las fuentes y los receptores no se muevan durante los múltiples barridos. Este procedimiento proporciona registros separados de alta calidad, porque la separación está bien planteada; existen tantos registros de entrada o barridos como registros de salida o sismogramas separados. Pero como son necesarios múltiples barridos, el procedimiento no es eficaz y los costos son mucho mayores que en el caso de los procedimientos con un solo barrido. La desventaja de realizar un solo barrido es que la separación no está bien planteada, y que habrá cierto ruido de interferencia residual tras la extracción de los sismogramas de fuentes. El problema de las interferencias se agrava por el hecho de que la firma o salida de los vibradores se relaciona de manera imperfecta con la señal piloto deseada por distorsión y la adición de armónicos y de que la señal real es desconocida. El ruido de interferencia se mitiga normalmente con un procedimiento de inversión y separación de datos iterativo (Neelamani, et al., patente estadounidense n.° 8.248.886) o mediante filtrado (Huo et al., publicación de patente estadounidense n.° 2012/0290214).
5
10
15
20
25
30
35
40
45
50
55
60
65
El documento US2007165486 se refiere a procedimientos y sistemas para adquirir datos sísmicos marinos con cable remolcado. Un procedimiento y sistema comprende colocar una pluralidad de embarcaciones de remolque de sólo fuente y una o varias embarcaciones de remolque de cable de fuente para adquirir una medición sísmica acimutal amplia y/o completa sin que sea necesario que el tendido repita un trayecto una vez atravesado. Otro procedimiento y sistema permite medir una característica geológica submarina utilizando un tendido sísmico marino, realizando el tendido al menos una vuelta durante la medición y disparando y registrando durante la vuelta.
La búsqueda simultánea de fuentes también puede utilizarse para fuentes impulsivas aunque hay menos procedimientos y procedimientos menos potentes para codificar las fuentes impulsivas. Hay poco beneficio de ahorro de costos para el uso de búsqueda simultánea de fuentes para la adquisición en tierra con dinamita, aunque el uso de una búsqueda simultánea de fuentes para cañones de aire en la adquisición marina puede resultar beneficioso, especialmente para la adquisición acimutal amplia. El uso de tiempos de detonación aleatorios para fuentes marinas que detonan casi simultáneamente pero están ubicadas en diferentes embarcaciones se dio a conocer por Vaage (patente estadounidense n.° 6.906.981). Más recientemente, se ha propuesto la búsqueda simultánea de fuentes para el disparo de embarcación múltiple de mediciones marinas acimutales amplias (WAZ) (Beasley et al., “A 3D simultaneous source field test processed using alternating projections: a new active separation method”, Geophysical Prospecting 60, 591-601 (2012)). La búsqueda simultánea de fuentes es la única manera de poder adquirir puntos de fuente con una separación fina (por ejemplo 25-m) en una única pasada de los cables. Sin la búsqueda simultánea de fuentes, son necesarias múltiples pasadas y la medición tarda mucho más y los costos son significativamente mayores.
Se ilustra una configuración para una medición marina WAZ, en la figura 1 para mostrar el beneficio de la búsqueda simultánea de fuentes. La figura muestra la línea de fuente 123, que atraviesa una embarcación de fuente, y la línea de receptor 121, que atraviesa una embarcación que tira de múltiples cables de hidrófonos. Las dos embarcaciones se mueven en paralelo a la misma velocidad, normalmente como mínimo de 6 nudos. En la figura, la primera posición de las embarcaciones se muestra en negro y las posiciones futuras se muestran en gris. La embarcación de fuente se dispara en la posición 103 mientras que la embarcación de receptor está en la posición 101, y la respuesta se registra normalmente durante aproximadamente 10 s. Durante estos 10 s, las embarcaciones están en movimiento. Algunos segundos más tarde, la embarcación de fuente alcanza el punto de tiro siguiente en 113, normalmente a 20-40 m del tiro anterior, y la embarcación de receptor alcanza la posición 111. Para registrar acimuts más amplios, la embarcación de receptor puede realizar 4 pasadas de la línea de receptor 121 mientras que una embarcación de fuente atraviesa las líneas de fuente 123-126 en secuencia. Se trata de una opción cara aunque puede proporcionar un muestreo fino de fuentes para cada línea de fuente, por ejemplo, un intervalo de fuente de 25-m. Alternativamente, pueden utilizarse 4 embarcaciones de fuente, y las fuentes que se disparan cambian entre las líneas. Por ejemplo, puede dispararse una fuente en la posición 103, a continuación en la posición 114, 135, 136, y a continuación 143. Si una de las fuentes dispara en modo cambiante cada 25-m, entonces el intervalo de fuentes a lo largo de cada línea (de 103 a 143) es de 200 m, mucho más basto que en el caso de una medición de separación fina. No es posible disparar y registrar con una separación de tiro más fina, porque cuando se ha adquirido el registro completo las embarcaciones se han movido decenas de metros a lo largo de las líneas de navegación. Puede registrarse una medición de separación fina con un disparo simultáneo detonando las cuatro fuentes dentro del mismo registro de tiempo pero con un pequeño retardo aleatorio o fluctuación en la posición o el tiempo de disparo. Por ejemplo, las fuentes pueden detonarse en las posiciones 103, 104, 105 y 106 para formar un registro con energía de fuente solapante. Entonces para el registro siguiente se detonan las fuentes en 113, 114, 115 y 116, etc.
La fluctuación es una forma de codificación que permite eliminar parcialmente la interferencia mediante filtrado en el procesamiento. Como las embarcaciones están en movimiento, un retardo en el tiempo de detonación significa necesariamente un ligero desplazamiento en la posición de detonación alrededor del intervalo de búsqueda de fuentes nominal como se determina mediante la velocidad de la embarcación. En lugar de requerir una sincronización en el tiempo entre embarcaciones puede ser más sencillo desde el punto de vista funcional implementar retardos de tiempo aleatorios generando un “patrón previo” de posiciones de búsqueda de fuentes a lo largo de cada línea con variaciones de posición aleatorias alrededor del intervalo de fuentes nominal. Durante la adquisición, cada embarcación dispara independientemente de las demás embarcaciones en posiciones de búsqueda de fuentes predeterminadas. Con este procedimiento, se predetermina la posición de detonación exacta pero no el tiempo de detonación, aunque el resultado sigue siendo una aleatorización en el tiempo. En la presente invención, se entiende que la aleatorización de la posición o el tiempo de búsqueda de fuentes son equivalentes. En cualquier caso, es importante determinar la posición de detonación y el tiempo de detonación real y estos valores junto con otras características de búsqueda de fuentes comprenden la función de codificación.
El registro de datos combinado obtenido con la búsqueda simultánea de fuentes debe separarse en registros individuales para cada fuente para un procesamiento convencional. En la figura 2 se muestra un diagrama de flujo del proceso estándar, y el proceso se ilustra en la figura 3. En la etapa 201, se obtienen registros de fuentes de longitud Treguo > Tescucha con múltiples excitaciones de fuente durante el registro. Durante la búsqueda de fuentes se utiliza algún tipo de esquema de codificación de campo tal como posición o tiempos de inicio con fluctuación. Tescucha es el tiempo necesario para que la energía se desplace desde la fuente al objetivo y luego al receptor. En la figura 3, 304 es una simple ilustración de un registro de fuente para una única fuente. En estos diagramas y en diagramas
5
10
15
20
25
30
35
40
45
50
55
60
65
posteriores, la respuesta de una única fuente se ilustra con un evento lineal 302 y un evento hiperbólico 303. En el registro de fuentes simultáneas 309, las cuatro fuentes detonan con pequeños retardos de tiempo y un evento lineal y uno hiperbólico de cada una de las cuatro fuentes interfieren. Supóngase que estas 4 fuentes están a distancias transversales largas en las líneas de fuente 123-126 en la figura 1, pero sólo se muestran una embarcación y una fuente en la vista de sección transversal de 309. A continuación, en la etapa 202, se determinan las funciones de codificación, incluyendo las posiciones de fuente y los tiempos de inicio, para las fuentes que contribuyen a los registros. El tiempo y la ubicación de fuente proyectados sobre la ventana de registro se indica mediante el rayo de sol 301 para la única fuente. La variación de posición de fuente correspondiente es relativamente pequeña en comparación con la escala de la figura y no se ilustra en el diagrama. El registro de fuentes simultáneas 309 se genera por fuentes en posiciones proyectadas 305, 306, 307 y 308, teniendo cada una un desplazamiento de tiempo reducido entre sí. Entonces, en la etapa 203, se utiliza la función de codificación para extraer registros de fuentes individuales, uno para cada fuente que empieza en el tiempo de detonación de esa fuente y que continúa durante el tiempo de escucha apropiado Tescucha. Para datos Vibroseis, esta extracción puede incluir el proceso de correlación mediante el barrido particular utilizado para esa fuente. En este ejemplo marino, el único registro 309 se copia 4 veces y a continuación se desplaza en el tiempo de modo que el registro comienza (tiempo cero) en el tiempo de detonación para cada fuente respectiva. Por ejemplo, el registro 315 se obtiene copiando el registro 309 y realizando un desplazamiento en el tiempo hasta el tiempo de la fuente 305. El registro 316 corresponde a la fuente 306; el registro 317 corresponde a 307 y 318 corresponde a 308. Esta etapa se denomina a veces pseudoseparación o pseudodescomposición. En este ejemplo, nada del ruido de interferencia se ha eliminado en esta fase. A continuación, en la etapa 204 se utilizan procedimientos de procesamiento adicionales para filtrar la energía de interferencia no deseada en cada agrupamiento de fuentes, o para utilizar un esquema de inversión dispersa para mejorar la separación de los datos, dando como resultado un sismograma separado para cada fuente puesto que se ha registrado independientemente de las demás fuentes. Entonces, en la etapa 205, es posible invertir u obtener imágenes de manera convencional de los agrupamientos de fuentes separados.
Puede utilizarse el mismo procedimiento de procesamiento indicado en la figura 2 para un registro en tierra o en el fondo del mar. Con datos de tierra o del fondo del mar, ahora es posible registrar la respuesta de los receptores de manera continua. Los receptores inalámbricos contienen memoria y un reloj y pueden registrar sin detenerse durante semanas o meses. Durante este tiempo, las fuentes detonan y se registran sus tiempos de detonación. Pueden utilizarse múltiples fuentes y para reducir el costo y el tiempo de adquisición, estas fuente pueden detonar de modo que la respuesta se solape en el tiempo como se muestra en la figura 4A. La figura 4A muestra un único registro continuo largo 401 con múltiples excitaciones de fuente ilustradas por rayos de sol. A diferencia del caso del cable marino, las fuentes no se detonan a pequeños intervalos en comparación con la longitud del registro, y por tanto, no se cumplen las condiciones para el uso del término disposición de área. Pero como en el caso del cable marino, la etapa de pseudoseparación inicial 203 implica extraer ventanas con el tamaño de la longitud del registro deseada Tescucha que empieza en el tiempo de detonación de una de las fuentes tal como se muestra en la figura 4B. Por ejemplo se copia una ventana correspondiente a 412 que comienza en el tiempo de detonación de fuente 402 y se extrae para obtener un registro de fuente 422. Tiene energía de interferencia de otras fuentes en 433 y 424. A continuación se copia una ventana 413 y se extrae para la fuente 403 para obtener el registro de fuente 423, y se extrae una ventana 414 para la fuente 404 obteniendo el registro de fuente 424, etc.
La búsqueda simultánea de fuentes seguida de la separación de fuentes también puede utilizarse como ayuda con simulación de datos sísmicos más cara desde el punto de vista computacional o modelado directo como se describe en Neelamani et al. (patente estadounidense n.° 8.248.886). Tal modelado directo es un componente de formación de imágenes sísmicas o inversión sísmica siendo la salida una imagen de reflectividad o de propiedades de formación tales como la velocidad sísmica de la subsuperficie. El modelado directo utiliza un modelo de velocidad detallado y calcula los campos de onda complejos generados teóricamente por cada fuente. Puede ahorrarse un tiempo de cálculo considerable reduciendo el número de fuentes a modelar a la vez utilizando la búsqueda simultánea de fuentes con algún tipo de esquema de codificación, y a continuación separando los datos en los sismogramas de fuentes individuales. Este procedimiento es idéntico a la adquisición de campo, aunque hay más opciones de esquemas de codificación cuando se realizan en el ordenador, y la secuencia codificada específica para una fuente se conoce perfectamente. Un esquema de codificación común es utilizar un ajuste a escala aleatorio en el que la salida de cada fuente se multiplica aleatoriamente por más o menos uno. Este esquema no puede implementarse físicamente en el campo para fuentes de impulso tales como cañones de aire o explosivos.
Como se describió anteriormente, puede utilizarse la búsqueda simultánea de fuentes para reducir los costos para adquirir datos sísmicos en el campo o para simular datos sísmicos en el ordenador. Esto implica registrar uno o varios registros compuestos que contienen interferencia de múltiples fuentes. Puede ser un registro corto con fuentes excitadas muy próximas en el tiempo y que forman una disposición de fuentes espacial. También puede ser un registro largo continuo con fuentes individuales excitadas a intervalos aleatorios o fijos. Para una formación de imágenes e inversión convencional, el registro compuesto debe separarse en agrupamientos de fuentes individuales. Normalmente, esto implica una pseudoseparación extrayendo una ventana alrededor del tiempo de detonación de las fuentes y a continuación el uso de operaciones de filtrado o inversión para eliminar la interferencia o ruido de interferencia. En el caso especial, en el que el número de registros es el mismo o mayor que el número de fuentes individuales dentro de una disposición espacial, la separación es bastante buena, aunque la adquisición de múltiples
5
10
15
20
25
30
35
40
45
50
55
60
65
registros es cara. Con menos registros, existe el problema de que la separación es imperfecta quedando algo de ruido de interferencia o habiendo eliminado una señal importante por el filtrado o la inversión.
Inversión sin separación
La búsqueda simultánea de fuentes también se utiliza para ahorrar costos computacionales asociados con la formación de imágenes y la inversión de datos sísmicos. En estos procedimientos, los agolpamientos de fuentes sísmicas individuales adquiridos secuencialmente, es decir, un tiro de fuente o disposición de fuente a la vez, se codifican en el ordenador y se suman para formar un registro de fuentes simultáneas que a continuación se utiliza para formar una imagen de reflectividad sísmica o para determinar propiedades de subsuperficie. El uso de este procedimiento para aumentar la velocidad y reducir el costo de la migración convencional (no iterativa y no mejora un modelo de subsuperficie) se da a conocer por Ober et al. (patente estadounidense n.° 6.021.094) y el uso del procedimiento en inversión se da a conocer por Krebs, et al. (patente estadounidense n.° 8.121.823). La interferencia entre fuentes también es un problema para este uso de la búsqueda simultánea de fuentes y tal interferencia se manifiesta como ruido en las salidas de inversión y formación de imágenes. La interferencia puede minimizarse en cierto modo optimizando las funciones de codificación por ordenador, tal como utilizando ajuste a escala aleatorio en lugar de rotación de fase, aunque los resultados pueden no ser tan buenos como con el uso secuencial de las fuentes individuales que requiere más cálculo.
La búsqueda simultánea de fuentes es particularmente útil para la inversión, tal como inversión de forma de onda completa (FWI) y migración de tiempo inverso de mínimos cuadrados (LSRTM). Estos procedimientos, a diferencia de la formación de imágenes tradicional, funcionan para actualizar de manera iterativa un modelo experimental para minimizar una función de desajuste de datos. El modelo son las propiedades de subsuperficie tales como velocidad para FWI, o la reflectividad para LSRTM. Obsérvese que la función de desajuste se calcula sin separación de fuentes. Como tanto el modelado directo como el procedimiento de actualización de modelo requieren mucho cálculo, la búsqueda simultánea de fuentes resulta muy ventajosa. Normalmente todas las fuentes en la medición o todas las fuentes en un barrido o línea de navegación se codifican y suman para obtener una disposición de fuentes simultáneas muy grande. Para minimizar el ruido de interferencia y para mejorar los resultados, las fuentes pueden volver a codificarse y volver a sumarse en cada iteración y a continuación utilizarse para una actualización de modelo (Krebs, patente estadounidense n.° 8.121.823). Cada grupo de datos codificados y sumados puede denominarse implementación de los datos. Se consiguen los mejores resultados y una interferencia reducida cuando se utilizan múltiples implementaciones en el proceso iterativo.
En la figura 5 se muestra un proceso típico para el uso de la búsqueda simultánea de fuentes en inversión. En la etapa 501, se obtienen varios registros de campo, cada uno con la misma extensión de tendido Ltendido y duración de registro Tescucha. La duración de registro Tescucha deberá incluir el tiempo necesario para que las ondas sísmicas se desplacen desde la fuente al objetivo y a continuación a los receptores. Para cada registro puede utilizarse una disposición de área o fuente única. A continuación se codifican los registros por ordenador, preferiblemente con un esquema de codificación aleatorizado en la etapa 502. Por ejemplo, los registros pueden multiplicarse aleatoriamente por más o menos 1 o aplicarse una rotación de fase por un factor aleatorio. A continuación se suman o acumulan todos los registros en la línea de navegación o barrido o en toda la medición, formando un registro de fuentes simultáneas. Esto se denomina implementación de los datos. A continuación, en la etapa 503, se simula la respuesta sísmica en el ordenador para el tiempo de desplazamiento Tescucha para todas las fuentes a la vez utilizando el esquema de codificación por ordenador. Esta etapa utiliza un modelo inicial o actualizado. Los registros simulados y medidos se comparan en la etapa 504, y se utiliza la comparación o función de desajuste para actualizar el modelo de propiedades o imágenes de subsuperficie. Si son necesarias múltiples iteraciones (etapa 506), es preferible volver a la etapa 502 y volver a codificar los registros de campo, llevando a cabo una segunda implementación de los datos. Cambiando la codificación de cada iteración, se reducen los artefactos y el ruido residual.
El uso de la búsqueda simultánea de fuentes para la inversión iterativa supone que el tendido de receptores y la longitud del registro son fijos, es decir, que todos los receptores registran para todas las fuentes durante la misma longitud de tiempo de modo que los registros puedan sumarse entre sí. El ordenador se utiliza para aplicar un modelo directo a todas las fuentes para todos los receptores como si se iniciaran al mismo tiempo o casi al mismo tiempo. Si los datos de fuente puntual no se registran con un tendido fijo, por ejemplo si se utilizan diferentes ubicaciones de receptor para registrar diferentes tiros, entonces el caso de modelado directo no coincide con el caso de datos de campo. Esto puede crear problemas porque la función de desajuste, la diferencia entre los datos de campo y los datos del modelado directo, estará dominada por la energía que falta entre los datos de modelado directo y los datos medidos y no podrá utilizarse para actualizar el modelo experimental. Los datos de campo registrados por un cable marino son particularmente problemáticos, porque el cable de receptor se mueve con la embarcación y no es fijo. Con un registro en tierra o en el fondo del mar se alcanza más comúnmente un tendido fijo, aunque incluso en este caso puede obtenerse un tendido rodante en el que las líneas de receptor activas cambian con la posición de fuente y pueden no cumplirse los supuestos de un tendido fijo.
La figura 6 ilustra el problema con la adquisición de datos convencionalmente con un cable marino y a continuación el uso de una búsqueda simultánea de fuentes para reducir el esfuerzo de cálculo requerido en la inversión. En 61,
5
10
15
20
25
30
35
40
45
50
55
60
65
se detona una fuente en la posición 602 y se capta un registro 601. En este ejemplo, la embarcación avanza entonces hacia la posición 604 y capta el registro 603 y a continuación a la ubicación 606 para el registro 605. Los receptores están en movimiento de modo que los receptores reales están en diferentes ubicaciones a lo largo de la línea de fuente. Si todas las trazas se disponen por sus verdaderas posiciones a lo largo de las líneas de navegación, se codifican y suman, se obtiene un agolpamiento de fuentes simultáneas 610. A continuación, si las tres fuentes (622, 623, 624, correspondientes a 602, 604 y 606) se excitan simultáneamente en el ordenador, se obtiene el registro simulado 612. Existe un desajuste inmediato entre los datos medidos y los simulados. Los datos medidos no incluyen trazas para receptores a la izquierda de cada fuente, ni miden los desplazamientos más largos. Por tanto, los eventos mostrados para los datos simulados, por ejemplo 644 y 645 y 646, faltan en el registro de datos medidos 610. Se han propuesto varios procedimientos para realizar la inversión en el caso del tendido no fijo incluyendo Rickett, et al. (publicación de patente estadounidense n.° 2012/0215506), que propone separar los datos simulados antes de calcular el desajuste, y Routh et al. (publicación de patente estadounidense n.° 2012/0143506) que propone utilizar la función objetivo de correlación cruzada. Estas dos soluciones son compromisos que no solucionan totalmente el problema. Puede haber errores en la separación de los datos para el primer enfoque, y la función objetivo de correlación cruzada es menos sensible a la información de amplitud en los datos en comparación con la función objetivo de mínimos cuadrados convencional.
Otros intentos publicados para tratar el fallo del supuesto de receptor fijo incluyen (1). “Hybrid method for full waveform inversion using simultaneous and sequential source method”, de Routh et al., patente estadounidense n.° 8.437.998; (2) “Simultaneous source encoding and source separation as a practical solution for full wavefield inversion”, de Routh et al., publicación estadounidense n.° 2012/0073825; (3) “Orthogonal source and receiver encoding”, de Krebs, et al., publicación estadounidense n.° 2013/0238246; (4) Haber et al., “An effective method for parameter estimation with PDE constraints with multiple right hand sides”, publicación previa - UBC en la dirección de Internet
http://www.math.ubc.ca/~haber/pubs/PdeOptStochV5.pdf (2010).
En esta sección se ha comentado la generación de un agrupamiento de fuentes simultáneas en el ordenador a partir de datos registrados secuencialmente en el campo. Krebs, et al. (patente estadounidense n.° 8.121.823) enseñó que pueden utilizarse los registros codificados de campo adquiridos con una disposición de fuentes de área codificada en un registro corto en la inversión tal como se adquiere, sin la etapa de separación comentada en la sección de “Adquisición y después separación” de este documento. Al no separar los datos, no se incluyen los errores de los procesos de separación en las etapas de inversión o formación de imágenes. Tales errores pueden incluir una pérdida o eliminación de determinadas componentes de reflexión que son importantes, por ejemplo pueden eliminarse difracciones de inmersión pronunciadas por error y limitar la capacidad para obtener imágenes nítidas de las terminaciones de capa con pequeños fallos. Sin embargo, sigue existiendo el problema de que en el campo no pueden conseguirse determinados procedimientos de codificación potentes disponibles en el ordenador, tales como ajuste a escala aleatorio. Además, si todas las fuentes se adquieren simultáneamente en el campo con un conjunto de funciones de codificación, el patrón de codificación es fijo y no puede cambiarse cada iteración para obtener múltiples implementaciones de los datos. Finalmente, no siempre se alcanzan los requisitos para utilizar la búsqueda simultánea de fuentes para la inversión cuando se utiliza en el campo la búsqueda simultánea de fuentes. Como se comentó anteriormente es un requisito que los datos se registren con un tendido fijo, no móvil de los receptores para una longitud de tiempo corta fija. Los problemas de los tendidos móviles tal como se ilustró anteriormente para fuentes marinas son incluso peores cuando se realiza una búsqueda simultánea de fuentes en el campo. Además, no es práctico utilizar simulación por ordenador para simular exactamente los datos según se adquieren de manera continua mediante receptores inalámbricos en tierra durante semanas, como se ilustra en la figura 4A. Para ser prácticos, la longitud del registro simulado en el ordenador deberá ser corta, del orden de algunos segundos. Del registro 4A no puede extraerse un registro de este tipo con un pequeño grupo de fuentes registradas durante un periodo de tiempo al tiempo que están aisladas de otras fuentes.
La presente invención utiliza la búsqueda simultánea de fuentes en el campo superando los problemas de los tendidos no fijos y los tiempos de registro largos para obtener una pluralidad de pseudorregistros de súperfuente que pueden codificarse y acumularse por ordenador para obtener múltiples implementaciones de los datos que pueden cambiarse con cada iteración de la inversión.
Sumario de la invención
Esta invención es un procedimiento para adquirir datos, a un costo de adquisición sísmica reducido, utilizando fuentes simultáneas en el campo, y a continuación construir pseudorregistros de fuente que cumplen mejor con los requisitos para utilizar la búsqueda simultánea adicional de fuentes para simulaciones por ordenador o modelado directo como parte de una inversión iterativa, tal como FWI (Full Wavefield Inversion) o LSRTM (Least-Squares Reverse Time Migration), con una reducción adicional de los costos computacionales. Al cumplir mejor con los requisitos de la búsqueda simultánea de fuentes para FWI o RTM, se reducen los artefactos y la interferencia en la salida. El procedimiento puede utilizarse para la adquisición con cable marino y otras geometrías de tendido no fijas para adquirir desplazamientos positivos y negativos y para mitigar el problema de la “falta de datos” para la FWI de fuentes simultáneas. También puede utilizarse para datos en tierra para superar los problemas con los tendidos móviles y registros continuos largos, donde un registro continuo largo significa un registro de datos demasiado largo para una simulación eficaz por ordenador.
5
10
15
20
25
30
35
40
45
50
55
60
65
Una primera forma de realización de la invención es un procedimiento según la reivindicación 1. Una segunda forma de realización de la invención es un procedimiento según la reivindicación 11.
La primera forma de realización de la invención descrita anteriormente puede utilizarse sin las características adicionales de la segunda forma de realización, por ejemplo cuando se procesan datos en los que todas las fuentes iluminan un tendido completo de receptores. El modelo de velocidad actualizado o ajustado resultante del presente procedimiento inventivo puede utilizarse para migrar los datos sísmicos para generar una imagen de la subsuperficie, o para otros fines de procesamiento e interpretación de datos sísmicos relativos a la exploración de hidrocarburos.
Breve descripción de los dibujos
La presente invención y sus ventajas se entenderán mejor haciendo referencia a la siguiente descripción detallada y los dibujos adjuntos en los que:
la figura 1 es un diagrama esquemático de las posiciones de cuatro embarcaciones de fuente y una embarcación de cable de receptor a medida que se mueven hacia la derecha a lo largo de líneas predeterminadas para una medición de adquisición acimutal amplia convencional;
la figura 2 es un diagrama de flujo de las etapas convencionales utilizadas para procesar datos adquiridos con una búsqueda codificada simultánea de fuentes en el campo y a continuación separados en el procesamiento;
la figura 3 es una ilustración de registros sísmicos adquiridos para datos marinos con una búsqueda codificada simultánea de fuentes según la figura 1, seguida del procedimiento de separación de fuentes de la figura 2;
las figuras 4A-4B ilustran un registro sísmico adquirido con la búsqueda codificada simultánea de fuentes para tierra (4A) y a continuación separado en agrupamientos de fuentes (4B) según el procedimiento de la figura 2;
la figura 5 es un diagrama de flujo de las etapas para utilizar la búsqueda codificada simultánea de fuentes en el ordenador para la inversión iterativa de datos sísmicos adquiridos con un tiro cada vez como se describe en el documento US 8.121.823;
la figura 6 ilustra la dificultad de utilizar las etapas en la figura 5 para utilizar la búsqueda simultánea de fuentes para la inversión utilizando adquisición marina de una sola fuente convencional, en la que los cables en movimiento violan el supuesto de receptor fijo necesario para el procedimiento de la figura 5;
la figura 7 es un diagrama de flujo que muestra las etapas básicas en una forma de realización del presente procedimiento inventivo;
la figura 8 es un diagrama que muestra un ejemplo de adquisición marina y la posición de embarcación de fuente y receptor y los cables en una forma de realización del presente procedimiento inventivo;
la figura 9 es una ilustración de una forma de realización marina de la presente invención que muestra los agrupamientos de fuentes tal como se adquieren y la construcción de un pseudorregistro de súperfuente;
la figura 10 es una ilustración de varios pseudorregistros de súperfuente que podrían construirse utilizando la forma de realización marina de la presente invención;
la figura 11 es una ilustración que muestra la construcción de un registro medido con fuentes simultáneas y la simulación por ordenador simultánea de todas las fuentes en una línea de navegación;
la figura 12 muestra una ilustración de dos registros de campo registrados en tierra o en el fondo del mar con un tendido en movimiento;
la figura 13 ilustra la generación de un tendido fijo a partir de los registros de campo de la figura 12;
la figura 14 es una ilustración de dos pseudorregistros de súperfuente construidos a partir del registro de campo combinado de la figura 13;
la figura 15 ilustra la codificación y acumulación por ordenador de los pseudorregistros de súperfuente de la figura 14 y su simulación en el ordenador con todas las fuentes detonando simultáneamente; y
la figura 16 muestra los resultados de una prueba de modelo 2D de una forma de realización marina de la presente invención.
5
10
15
20
25
30
35
40
45
50
55
60
65
Por las restricciones de la ley de patentes respecto al uso del color, la figura 16 es una reproducción en blanco y negro de un original en color. La invención se describirá junto con formas de realización de ejemplo. Sin embargo, en la medida en que la siguiente descripción detallada sea específica de una forma de realización particular o un uso particular de la invención, será únicamente a modo ilustrativo y no deberá interpretarse como limitativo del alcance de la invención. Por el contrario, se pretende cubrir todas las alternativas, modificaciones y equivalentes que puedan incluirse en el alcance de la invención, tal como se define en las reivindicaciones adjuntas.
Descripción detallada de las formas de realización de ejemplo
La invención se describe en primer lugar en su forma básica, a continuación se describirán formas de realización específicas para datos marinos y en tierra. Esta invención utiliza la búsqueda simultánea de fuentes en el campo para mejorar la capacidad para utilizar adicionalmente la búsqueda simultánea de fuentes en la inversión iterativa reduciendo los efectos del ruido de interferencia y aproximándose mejor a la adquisición mediante un tendido fijo de receptores. La invención construye, a partir de registros de adquisición, lo que pueden denominarse pseudorregistros de súperfuente (o súpertiro), cada uno con la misma duración y extensión espacial. Cada pseudorregistro de súperfuente contiene la energía registrada de múltiples fuentes, estando energizada cada fuente con un esquema de codificación de campo (por ejemplo, desplazamientos en el tiempo aleatorios, posiciones de fuente aleatorias, rotaciones de fase, función de barrido, u otros procedimientos) y cada registro se construye mediante las operaciones de función ventana, desplazamiento en el tiempo, suma y adición de los registros de campo originales. La medición se adquiere de modo que permite construir estos pseudorregistros de súpertiro de modo que las fuentes puedan simularse de manera apropiada simultáneamente en un ordenador. En particular, la energía sísmica que se generaría por una fuente sintética y registrada dentro de una distancia predeterminada Dfuente dentro del tendido y la duración de tiempo Tfuente se representa dentro del pseudorregistro de súpertiro medido. Este requisito puede requerir la repetición de algunos grupos de puntos de fuente en diferentes tendidos de receptor con la misma codificación que la utilizada previamente. Los múltiples registros de súpertiro se codifican entonces por separado en el ordenador, preferiblemente con un ajuste a escala aleatorio tal como multiplicando por +1 o -1 seleccionado aleatoriamente, y a continuación se suman y utilizan para la inversión o formación de imágenes. Preferiblemente, el esquema de codificación por ordenador se cambia en iteraciones posteriores de la inversión de la inversión o formación de imágenes.
En el diagrama de flujo de la figura 7 se proporcionan las etapas básicas en una forma de realización del presente procedimiento inventivo. En la etapa 701, se planifica una medición sísmica que incluye elegir el tipo de fuente y un número y separación óptimos de las ubicaciones de búsqueda de fuentes deseadas y un esquema de codificación de campo. El esquema de codificación de campo puede incluir cualquier parámetro relacionado con la excitación de fuente incluyendo su ubicación, tiempo de detonación, componentes de frecuencia, fase de fuente, etc. Además, se seleccionan dos parámetros: Tfuente y Dfuente. Tfuente es un tiempo mayor que o igual a Tescucha y, además, es un tiempo lo suficientemente largo como para que la amplitud de esa fuente se reduzca hasta un nivel en el que ya no proporciona un ruido de interferencia sustancial. Normalmente, Tfuente puede seleccionarse para ser igual a Tescucha a menos que se haya adquirido un registro continuo largo tal como en la figura 12. Dfuente está relacionado porque es igual o mayor que la distancia de interés de desplazamiento máximo o igual a la distancia en la que se reduce la amplitud de fuente de modo que ya no proporciona un ruido de interferencia sustancial. Estos parámetros se utilizan para guiar la adquisición y posterior generación de pseudorregistros de modo que los pseudorregistros puedan simularse con precisión dentro de una distancia Dfuente y tiempo Tfuente de cualquier ubicación de fuente y tiempo de detonación. De manera ideal los parámetros son lo más pequeños posible al tiempo que son lo suficientemente grandes como para incluir retornos desde el objetivo. Parámetros más pequeños significan requisitos menos estrictos para la adquisición y construcción de pseudoagrupamientos.
En la etapa 702, se obtienen uno o varios registros de campo que se generan con la búsqueda “simultánea” de fuentes de modo que la energía de las diferentes fuentes se solape parcialmente en el tiempo. Dicho de otro modo, las fuentes no tienen que activarse exactamente de manera simultánea, y los desplazamientos de tiempo reducidos entre las mismas son una manera de realizar la codificación de campo a la que se hace referencia en la etapa 702. Un registro de campo es normalmente el conjunto, o un subconjunto, de los datos registrados por los receptores activos (en movimiento o estacionarios) en un periodo de tiempo, con un tiempo de inicio y un tiempo de parada y sin espacios. Los registros de campo pueden ser registros discretos de una duración de tiempo fija o pueden ser un único registro de tiempo continuo. Si se mueve el tendido de registro durante la adquisición, entonces preferiblemente se repiten algunos de los puntos de fuente dentro de la distancia Dfuente del límite del primer tendido en el segundo tendido con el mismo esquema de codificación utilizado previamente de modo que toda la energía dentro de la distancia Dfuente se registre en ambos conjuntos de posiciones de tendido de modo que puedan añadirse entre sí.
Entonces, en la etapa 703, se construyen una pluralidad de lo que pueden denominarse pseudorregistros de súperfuente de extensión y duración fija. Preferiblemente, la extensión del registro abarcaría la anchura de medición, como si la medición se hubiera registrado mediante un tendido fijo de receptores por la anchura de la medición, y la duración de registro sería al menos igual de larga que el tiempo para que la energía sísmica se propague desde la fuente al objetivo y a los receptores con el desplazamiento o la distancia útil máxima desde la fuente. El proceso de construcción puede incluir operaciones tales como extracciones de diversas ventanas de tiempo y regiones de traza
5
10
15
20
25
30
35
40
45
50
55
60
65
a partir de los registros de campo. Además, puede agregarse un relleno de trazas cero y puede añadirse un relleno de tiempo antes o después de las ventanas. Un objetivo de la construcción de un pseudorregistro de súperfuente es que cada ubicación de receptor dentro de una distancia de desplazamiento Dfuente con respecto a la ubicación de la fuente tenga datos apropiados, es decir, datos que se habrían registrado si hubiera habido un tendido de receptores fijo cuando se produjo el tiro de fuente. Normalmente, los datos de cada registro de campo aparecerán en al menos un pseudorregistro de súperfuente. Entonces las diversas ventanas pueden añadirse o sumarse entre sí para formar un pseudorregistro de súpertiro. Entonces, en la etapa 704, cada tiro que influye en o contribuye a la región de interés se identifica junto con su función de codificación de campo, y tiempo de inicio con respecto al tiempo cero del pseudorregistro de súpertiro. Puede suponerse que los tiros de contribución o influencia son aquéllos para los cuales la fuente se excita dentro de la distancia Dfuente y un tiempo Tfuente desde el límite de la región de interés. Esta información se combina con la función de codificación por ordenador y se utiliza para la etapa de simulación por ordenador 706.
A continuación, en la etapa 705, se codifican por ordenador los diferentes pseudorregistros de súpertiro, preferiblemente, aunque no necesariamente (funcionará cualquier esquema de codificación incoherente), mediante ajuste a escala aleatorio en el que se multiplican aleatoriamente por más o menos uno (+1). Entonces se suman entre sí todos los pseudorregistros para formar un registro de fuentes simultáneas. Entonces se utiliza el ordenador para calcular la simulación de modelado directo en una etapa para todas las fuentes dentro del registro simultáneo, que se identificaron en la etapa 704, como si todas las fuentes se hubieran detonado simultáneamente o casi simultáneamente (etapa 706). (Dicho de otro modo, se resuelve una ecuación de propagación de ondas con condiciones iniciales y límite apropiadas y suponiendo un modelo de velocidad de subsuperficie, utilizando procedimientos numéricos tales como diferencia finita iterativa). La simulación por ordenador se realiza utilizando los esquemas de codificación por ordenador y de campo combinado, es decir, lo que podría denominarse codificación doble. Cuando en la etapa 706 se simula un registro de fuentes simultáneas correspondiente a un registro medido con fuentes simultáneas de la etapa 705, se genera el registro de fuentes simultáneas simulado utilizando una combinación de la codificación por ordenador que se utilizó en la etapa 705 combinada con la codificación de campo de la etapa 704/702. En la etapa 707, se comparan los registros registrados de 703 y los registros simulados de 706 por una región de interés, y se utilizan los resultados para actualizar el modelo de subsuperficie. Si son necesarias más iteraciones de la formación de imágenes o inversión como se determina en la etapa 708, entonces preferiblemente se repite la etapa de codificación por ordenador 705 con una función de codificación diferente.
Forma de realización marina
En esta sección, se describe una forma de realización particular que supera el problema del tendido en movimiento para un cable marino, que se ilustró en la figura 6. Esto explicará adicionalmente con una ilustración lo que quiere decirse con planificar una adquisición para permitir la pseudoconstrucción de súperfuente en la etapa 701, y cómo se construyen los pseudorregistros de súperfuente en la etapa 703. En esta forma de realización, se utiliza un patrón de detonación y geometría de adquisición de modo que pueden construirse los agrupamientos de súpertiro y aproximarse mejor a la geometría de tendido fijo necesaria para la simulación por ordenador. En la etapa 701, se planifica la adquisición lo que modifica la geometría de adquisición mostrada en la figura 1, ubicando una fuente en la parte posterior del cable para cada fuente en la parte anterior del cable tal como se muestra en la figura 8. A medida que avanzan las embarcaciones de fuente y cable, las fuentes posteriores siguen la misma línea de fuente que las fuentes anteriores. Obsérvese que las fuentes 803 y 807 siguen la línea de fuente 801, y las fuentes 804 y 808 siguen la línea de fuente 802. Esto es diferente de las geometrías acimutales amplias convencionales que pueden utilizar una fuente posterior (Treadgold, et al., “Implementing A wide Azimuth Towed Streamer Field Trial, The What, Why, and Mostly How of WATS in Southern Green Canyon”, SEG Expanded Abstracts, 2901-2903 (2006)), pero ubica las fuentes posteriores en diferentes líneas de navegación. Preferiblemente las distancias de desplazamiento cercanas son las mismas para las fuentes anterior y posterior. Por ejemplo, la distancia desde la fuente 803 al receptor más cercano 805 es la misma que desde la fuente 807 a su receptor más cercano 806.
Todas las fuentes detonan en el mismo intervalo de fuentes pero con diferentes retardos de tiempo aleatorios o posiciones aleatorias alrededor de la ubicación de fuente nominal, y se registra un único registro de longitud fija como se ilustra como 901 en la figura 9. En este ejemplo, dos embarcaciones de fuente, 803 y 804 están delante del cable y dos embarcaciones de fuente 807 y 808 detonan en la parte posterior del cable. La posición del tiro proyectado sobre el registro de fuente desde dos fuentes anteriores se muestra mediante los rayos de sol en 903. Las embarcaciones avanzan, detonando las cuatro fuentes cada intervalo de tiro, por ejemplo cada 25-m, con retardos de tiempo aleatorios o una fluctuación de posición aleatoria para cada grupo de tiros. Este retardo de tiempo o esta fluctuación de posición aleatorios se consideran una función de codificación. Cuando las fuentes posteriores se aproximan a la posición ocupada previamente por las fuentes anteriores 903, entonces no detonan aleatoriamente, sino que detonan lo más cerca de la misma posición que sea posible funcionalmente y con los mismos retardos de tiempo con respecto a otras fuentes posteriores que tenían las fuentes anteriores entre sí en la misma ubicación. Por ejemplo, el registro de fuente 906 se registra con las fuentes posteriores detonando con el mismo patrón 904 que cuando esta posición estaba ocupada por las fuentes anteriores para el registro 903 en 901. Y de nuevo cuando las fuentes posteriores alcanzan la posición 907 detonan con la misma función de codificación y en la misma ubicación que cuando las fuentes anteriores estaban en 905; véase el registro 908. Este patrón continúa preferiblemente para toda la línea de navegación.
5
10
15
20
25
30
35
40
45
50
55
60
65
En la etapa 703, se construye un pseudorregistro de súperfuente. Cada registro que se registró con fuentes posteriores y anteriores idénticas en la misma posición tiene una alineación en el tiempo para coincidir con la sincronización de fuente y se añaden. Las trazas pueden rellenarse (es decir, pueden añadirse ceros) al final o al principio. Esto da lugar a un registro largo 912 que abarca toda la línea de navegación con fuentes separadas por aproximadamente la longitud del cable. Ahora, este pseudorregistro de súperfuente se aproxima mejor a un tendido fijo porque se registran los desplazamientos tanto positivos como negativos de cada posición de fuente hasta una distancia de Dfuente. Aquí Dfuente es naturalmente la longitud de cable. Ahora todas estas fuentes pueden simularse simultáneamente en el ordenador, por ejemplo poniendo grupos de fuentes en 922, 932, 942 y 952.
En la etapa 703, se construyen pseudorregistros de súperfuente adicionales, teniendo cada uno la misma extensión espacial y duración de tiempo que se ilustra en la figura 10. Si el primer tiro fue en la posición 0, entonces comienza el primer registro de súperfuente. En el siguiente intervalo, por ejemplo 25-m con respecto al primero, comienza el segundo registro de súperfuente, sus tiros están desplazados un intervalo de tiro con respecto al primero. Esto continúa hasta que se utiliza toda la longitud de cable y las fuentes posteriores alcanzan la primera ubicación de las fuentes anteriores. Por ejemplo, los registros de súperfuente 1001-1004 tienen posiciones de fuente que están desplazadas por un intervalo de fuentes. En general, el número de registros de súpertiro que pueden generarse está relacionado con la longitud del cable dividida por el intervalo de tiro. Puede ser necesario rellenar trazas, por ejemplo en 1006, para construir un tamaño de tendido fijo para cada súperregistro.
En la etapa 704, se determinan la información de codificación y ubicación de fuente, incluyendo desplazamientos de tiempo, para cada pseudorregistro de súperfuente con respecto a los límites del pseudorregistro. Por ejemplo, el tiempo de inicio de cada fuente se ajusta por los desplazamientos de tiempo utilizados para formar el pseudorregistro y ahora se relaciona con el tiempo cero del pseudorregistro. Esta información se utilizará en la etapa 706, combinada con la codificación por ordenador utilizada en la etapa 705, para simular simultáneamente los pseudorregistros codificados.
En la etapa 705, cada pseudorregistro de súperfuente que contiene muchos tiros se codifica en el ordenador. Preferiblemente esto se realiza multiplicando aleatoriamente por +1 o -1. Alternativamente, pueden utilizarse rotaciones de fase u otra forma de codificación. Entonces se acumulan o suman los pseudorregistros codificados, tal como se muestra en la ilustración de 1112 en la figura 11, para formar un registro de fuentes simultáneas para toda la línea de navegación.
En la etapa 706, todas las fuentes en la línea de navegación se simulan por ordenador a la vez utilizando una combinación de la codificación de campo determinada en la etapa 704/702 y la codificación por ordenador utilizada en la etapa 705. Pueden conseguirse ahorros adicionales en el costo computacional limitando la región del modelo utilizado en una simulación de línea de navegación única. Esta simulación se ilustra con los rayos de sol en 1116. Esto puede implicar extender o rellenar el espacio de modelado con una región adicional tal como se indica en 1113, lo que permite que el modelado directo genere todos los bits de energía registrados en la ventana de datos 1118. A continuación, en la etapa 707, se compara el registro de fuentes simultáneas medido por la región de interés (1118 en el ejemplo) con el registro de fuentes simultáneas simulado y se utiliza el resultado para actualizar el modelo de subsuperficie. Como la energía simulada de todas las posiciones de fuente que influyen en la región de interés está presente en los datos tanto medidos como simulados, al menos hasta una distancia y un tiempo de Dfuente y Tfuente desde cada posición de detonación de fuente, se evita el problema de que la función de desajuste o comparación se distorsione por artefactos por crear el registro de fuentes simultáneas. Construyendo los pseudorregistros, se cumple con los requisitos de un tendido fijo y una duración de traza reducida. Con iteraciones adicionales de la etapa de inversión o formación de imágenes, la codificación por ordenador se cambia preferiblemente repitiendo la etapa 705 con una codificación diferente y formando así múltiples implementaciones de los datos que reducen adicionalmente la interferencia y los artefactos.
Es sencillo modificar el ejemplo mostrado en este caso para diferentes requisitos de adquisición que implican más embarcaciones de cable y conjuntos de cables y más embarcaciones de fuente. La posición de fuente y el patrón de detonación se eligen de modo que puedan construirse pseudoagrupamientos de súperfuente aproximándose mejor al registro por un tendido fijo para la línea de navegación o para toda la medición. Por ejemplo, pueden utilizarse fuentes de remolque de embarcaciones adicionales delante de la embarcación de cable y más atrás del extremo del cable y detonarse simultáneamente con las demás fuentes para registrar datos de desplazamiento más largo para la línea de navegación. Además, las embarcaciones pueden ubicarse en ambos lados de las embarcaciones de cable para aproximarse mejor a un tendido fijo transversal permitiendo que todas las fuentes en la medición se simulen simultáneamente en un esfuerzo de modelado directo computacional.
Forma de realización en tierra
A continuación se considerará la adquisición de datos en tierra o en el fondo del mar. En este caso los sensores no se mueven como en un cable marino, sino que el grupo de receptores activos puede cambiar durante la medición. Esto se denomina a menudo rodamiento del tendido. En la figura 12, se muestra un registro en una posición de tendido como 1212 y el segundo registro correspondiente a un conjunto diferente de posiciones de traza se muestra
5
10
15
20
25
30
35
40
45
50
55
en 1213. A menudo los receptores pueden registrar realizando de manera continua un registro largo para cada posición de tendido. En la etapa 701, se planifica la medición y se identifica Tfuente y Dfuente. Puede activarse la búsqueda de fuentes aleatoriamente o en un patrón. Preferiblemente, algo de la energía generada de múltiples posiciones de fuente se solapa en el tiempo y el espacio para reducir el tiempo y costo de adquisición. Cuando se mueve el tendido, se repiten los puntos de fuente dentro de una distancia de Dfuente desde el límite para el nuevo tendido utilizando los mismos desplazamientos de tiempo relativos. Por ejemplo si las fuentes 1201, 1202 y 1203 están dentro de la distancia Dfuente desde el límite de modo que no se registre toda la energía dentro de la distancia Dfuente desde una fuente, se repiten estas fuentes en la nueva posición de tendido con el mismo desplazamiento de tiempo relativo y la misma función de codificación utilizada previamente en esa posición de fuente. Por ejemplo 1221, 1222 y 1223 son una repetición de 1201, 1202 y 1203. Esto asegura que el registro combinado ha captado toda la señal importante, tal como 1224, de cada fuente.
En la etapa 703, se construyen varios pseudorregistros de súperfuente de tamaño y duración fijos. Preferiblemente la duración es más larga que Tescucha pero lo suficientemente corta para simularse de manera eficaz en el ordenador. La construcción para este ejemplo puede ilustrarse en dos etapas. En la figura 13 se combinan los dos registros de la figura 12 para obtener un tendido fijo. Como se han repetido algunos tiros cerca del límite con la misma función de sincronización y codificación relativa, se ha garantizado que toda la energía dentro de una distancia de Dfuente y un tiempo de Tfuente se capte en el registro combinado. A continuación, en la figura 14 se extraen ventanas de tiempo aisladas preferiblemente de longitud Tfuente. No es necesario aislar la energía medida de ninguna fuente dentro de la ventana, ni es necesario que ningún registro comience en el tiempo de detonación de ninguna fuente como en el caso de la pseudodescomposición. Por ejemplo se extraen las ventanas 1404 y 1406 del registro en la figura 13. Se ilustran sólo dos registros, aunque preferiblemente se extrae un número mayor. Se rellenan los registros en tiempo negativo por otra duración de tiempo de Tfuente tal como se muestra como 1403 y 1405 para obtener pseudosúperregistros 1401 y 1402.
Entonces, en la etapa 704, se identifican las ubicaciones de fuente, la función de codificación y el tiempo de detonación con respecto al tiempo cero de los registros de súperfuente 1401 y 1402. Las fuentes deberán estar dentro de un tiempo de Tfuente desde la parte superior de la ventana o dentro de una distancia de Dfuente desde los límites del registro. En la etapa 705, los registros de súperfuente se codifican y suman por ordenador, obteniendo un registro de fuentes simultáneas medido como se ilustra en la figura 15 como 1501. Entonces, en la etapa 706, se utilizan las posiciones de búsqueda de fuentes y las funciones de codificación de campo y ordenador para generar un registro de fuentes simultáneas simulado como se ilustra en 1502. Esta simulación es eficaz porque todas las fuentes se simulan simultáneamente durante un breve periodo de tiempo. A continuación, en la etapa 707, se compara la región de interés medida 1504 con la región de interés simulada 1506, y se utiliza el resultado para actualizar un modelo de subsuperficie. Como toda la energía dentro de la distancia Dfuente y Tfuente de cada posición de fuente se representa en ambos registros, se evitan los artefactos de la creación del registro simultáneo. Con iteraciones posteriores de la formación de imágenes e inversión, preferiblemente los registros de súperfuente se combinarían con diferentes funciones de codificación formando diferentes implementaciones de todos los datos, reduciendo el ruido de interferencia.
Ejemplos
La figura 16 muestra un ejemplo de modelo de velocidad en 2D. Por las restricciones de la ley de patentes respecto al uso del color, la figura 16 es una reproducción en blanco y negro de una visualización de datos en la que la velocidad se representa cuantitativamente en una escala de colores. Los datos se han simulado por ordenador con una fuente anterior y posterior y codificación de campo con retardos de tiempo aleatorios de hasta 200 ms utilizando el modelo de superficie cercana mostrado en 1601. El uso del procedimiento convencional para la inversión de fuentes codificadas simultáneas con los datos de cable en movimiento no proporcionará los resultados correctos (no mostrado) debido al fallo del supuesto de receptor fijo. Entonces se construyen los agrupamientos de súperfuente como se describe en la forma de realización marina, y se realiza una inversión de baja frecuencia utilizando el presente procedimiento inventivo con los resultados mostrados en 1602. Puede observarse que el modelo deducido por la inversión de datos utilizando el presente procedimiento inventivo es comparable de manera muy favorable con el modelo “real” 1601.
La aplicación anterior se refiere a formas de realización particulares de la presente invención para su ilustración. Los expertos en la técnica reconocerán fácilmente que en formas de realización preferidas de la invención, al menos algunas de las etapas en el presente procedimiento inventivo se realizan en un ordenador, es decir, la invención se implementa por ordenador.

Claims (15)

  1. 5
    10
    15
    20
    25
    30
    35
    40
    45
    50
    55
    60
    65
    REIVINDICACIONES
    1. Un procedimiento para deducir un modelo de subsuperficie de velocidad u otra propiedad física, que comprende:
    obtener (702) datos sísmicos adquiridos en una medición haciendo avanzar una fuente anterior y una fuente posterior de seguimiento a lo largo de una línea de fuente con un tendido no fijo de receptores, y activar cada fuente a intervalos seleccionados a lo largo de la línea de fuente, en el que cada ubicación de activación para la fuente posterior se sitúa sobre una ubicación de activación previa para la fuente anterior, y en el que se utiliza una codificación de campo seleccionada para todos los tiros de fuente siendo la codificación de campo para la fuente posterior siempre la misma que la aplicada anteriormente a la fuente anterior en la misma ubicación de activación, dando como resultado una serie de registros de tiro, registrando cada uno simultáneamente ambas fuentes pero con cada tiro registrado por menos de un tendido completo de receptores; y
    procesar los registros de tiro por etapas que comprenden:
    (i) construir (703) a partir de los registros de tiro una pluralidad de pseudorregistros de súpertiro, construidos de modo que cada uno tiene datos de un tendido completo de receptores, en el que cada pseudorregistro de súpertiro incluye una pluralidad de ubicaciones de fuente;
    (ii) codificar (705) cada pseudorregistro de súpertiro y acumularlo para formar un registro de fuentes simultáneas de datos medidos;
    (iii) utilizar un ordenador para simular (706) el registro de fuentes simultáneas de datos medidos, utilizando la misma codificación utilizada en (ii) y también la codificación de campo, y utilizando un modelo de subsuperficie supuesto de velocidad u otra propiedad física; y
    (iv) comparar (707) el registro de fuentes simultáneas simulado con el registro de fuentes simultáneas de datos medidos, y determinar a partir de aquí un ajuste al modelo de subsuperficie de velocidad u otra propiedad física.
  2. 2. El procedimiento según la reivindicación 1, en el que los datos de medición adquiridos incluyen al menos una línea de fuente adicional, adquirida simultáneamente, paralela, en el que cada línea de fuente tiene una fuente anterior y una fuente posterior de seguimiento con ubicaciones de activación de fuente coincidentes para todas las líneas de fuente, en el que cada ubicación de activación para una fuente posterior en cualquier línea de fuente se sitúa sobre una ubicación de activación previa para la fuente anterior en esa línea de fuente, y en el que se utiliza una codificación de campo seleccionada para todos los tiros de fuente siendo la codificación de campo para cada fuente posterior siempre la misma que la aplicada anteriormente a su fuente anterior correspondiente en la misma ubicación de activación.
  3. 3. El procedimiento según la reivindicación 1, en el que la fuente anterior y la fuente posterior están separadas por una distancia d fija que subdivide en un número entero N de intervalos de tiro, siendo tales intervalos de tiro el intervalo de tiro nominal o promedio para la medición.
  4. 4. El procedimiento según la reivindicación 3, en el que se construyen N + 1 pseudorregistros de súpertiro en (i), cada uno de duración y extensión espacial fijas, siendo la duración mayor que o igual a un tiempo de escucha utilizado en la medición y, además, es un tiempo mínimo para que las vibraciones de tiro de fuente se reduzcan hasta un nivel predeterminado; y siendo la extensión espacial mayor que o igual a una distancia de desplazamiento de receptores máxima utilizada en la medición y, además, es una distancia mínima por la que se reduce la amplitud de fuente hasta un nivel predeterminado.
  5. 5. El procedimiento según la reivindicación 4, en el que se repiten (i)-(iii) una pluralidad de veces utilizando otros registros de tiro de otras ubicaciones de tiro a lo largo de la línea de fuente, y a continuación se realiza (iv) utilizando los resultados de todas las repeticiones de (i)-(iii).
  6. 6. El procedimiento según la reivindicación 5, en el que todos los datos adquiridos a lo largo de la línea de fuente se incluyen en al menos un pseudorregistro de súpertiro.
  7. 7. El procedimiento según la reivindicación 1, que comprende además tomar el modelo de velocidad ajustada de (iv) y repetir las etapas (i)-(iv) utilizando una codificación por ordenador diferente en (ii).
  8. 8. El procedimiento según la reivindicación 1, en el que la codificación de campo comprende uno o varios de desplazamientos en el tiempo aleatorios, posiciones de fuente aleatorias, rotaciones de fase y diferentes funciones de barrido.
  9. 9. El procedimiento según la reivindicación 1, en el que la construcción en (i) comprende utilizar uno o varios de función ventana, desplazamiento en el tiempo, relleno de tiempo o trazas y adición.
    5
    10
    15
    20
    25
    30
    35
    40
    45
    50
    55
  10. 10. El procedimiento según la reivindicación 1, en el que la medición sísmica comprende:
    una medición marina en la que una embarcación en movimiento remolca una o varias fuentes y uno o varios cables de receptores; o
    una medición en tierra o del fondo del mar que implica un tendido rodante de receptores a medida que avanza la medición a lo largo de una línea de medición.
  11. 11. Un procedimiento para llevar a cabo una medición sísmica con una pluralidad de fuentes simultáneas y un tendido no fijo de receptores, y deducir un modelo de subsuperficie de velocidad u otra propiedad física, que comprende:
    (a) activar dos o más fuentes con una codificación de campo seleccionada aplicada a cada fuente en una ubicación a lo largo de una línea de fuente con un primer tendido parcial de receptores;
    (b) activar las dos o más fuentes, o fuentes idénticas, duplicando la codificación de campo en (a), en la misma ubicación a lo largo de la línea de fuente con un segundo tendido parcial de receptores, en el que el primer tendido parcial y el segundo tendido parcial, cuando se combinan, forman un tendido completo de receptores para la ubicación a lo largo de la línea de fuente;
    (c) repetir (a) y (b) para una o varias ubicaciones adicionales a lo largo de la línea de fuente; a continuación procesar los registros de datos sísmicos registrados por los receptores de las activaciones en (a)-(c) con etapas que comprenden:
    (d) construir a partir de los registros de tiro una pluralidad de pseudorregistros de súpertiro, construidos para ser de la misma duración y extensión espacial, y de modo que cada uno tiene datos de un tendido completo de ubicaciones de receptor, en el que cada pseudorregistro de súpertiro incluye una pluralidad de ubicaciones de fuente;
    (e) codificar por ordenador cada pseudorregistro de súpertiro y sumarlos para formar un único registro de fuentes simultáneas de datos medidos;
    (f) utilizar un ordenador y un modelo de subsuperficie de velocidad u otra propiedad física, simulando con una operación de modelado directo un registro de fuentes simultáneas de datos sintéticos para predecir el registro de fuentes simultáneas de datos medidos; y
    (g) comparar el registro de fuentes simultáneas de datos medidos con el registro de fuentes simultáneas simulado, determinar un desajuste y utilizar el desajuste para ajustar el modelo de subsuperficie.
  12. 12. El procedimiento según la reivindicación 11, que comprende además tomar el modelo de subsuperficie ajustado de (g) y repetir las etapas (e)-(g) utilizando una codificación por ordenador diferente en (e).
  13. 13. El procedimiento según la reivindicación 11, en el que la codificación de campo comprende uno o varios de desplazamientos en el tiempo aleatorios, posiciones de fuente aleatorias, rotaciones de fase y diferentes funciones de barrido.
  14. 14. El procedimiento según la reivindicación 11, en el que la construcción en (d) comprende utilizar uno o varios de función ventana, desplazamiento en el tiempo, relleno de tiempo o trazas y adición.
  15. 15. El procedimiento según la reivindicación 11, en el que la medición sísmica comprende:
    una medición marina en la que una embarcación en movimiento remolca una o varias fuentes y uno o varios cables de receptores; o
    una medición en tierra o del fondo del mar que implica un tendido rodante de receptores a medida que avanza la medición a lo largo de una línea de medición.
ES14744730.4T 2013-08-23 2014-07-11 Búsqueda simultánea de fuentes durante la adquisición sísmica y la inversión sísmica Active ES2671789T3 (es)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201361869292P 2013-08-23 2013-08-23
US201361869292P 2013-08-23
PCT/US2014/046361 WO2015026451A2 (en) 2013-08-23 2014-07-11 Simultaneous sourcing during both seismic acquisition and seismic inversion

Publications (1)

Publication Number Publication Date
ES2671789T3 true ES2671789T3 (es) 2018-06-08

Family

ID=51230232

Family Applications (1)

Application Number Title Priority Date Filing Date
ES14744730.4T Active ES2671789T3 (es) 2013-08-23 2014-07-11 Búsqueda simultánea de fuentes durante la adquisición sísmica y la inversión sísmica

Country Status (9)

Country Link
US (1) US9772413B2 (es)
EP (2) EP3036566B1 (es)
AU (1) AU2014309376B2 (es)
BR (1) BR112015030104A2 (es)
CA (1) CA2913496C (es)
DK (1) DK3036566T3 (es)
ES (1) ES2671789T3 (es)
MX (1) MX359179B (es)
WO (1) WO2015026451A2 (es)

Families Citing this family (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9188693B2 (en) * 2012-03-26 2015-11-17 Apache Corporation Method for acquiring marine seismic data
US10422922B2 (en) * 2012-05-24 2019-09-24 Exxonmobil Upstream Research Company Method for predicting rock strength by inverting petrophysical properties
US20150063064A1 (en) * 2013-09-03 2015-03-05 Pgs Geophysical As Methods and systems for attenuating noise in seismic data
US10036818B2 (en) 2013-09-06 2018-07-31 Exxonmobil Upstream Research Company Accelerating full wavefield inversion with nonstationary point-spread functions
WO2015168130A1 (en) 2014-04-28 2015-11-05 Westerngeco Llc Wavefield reconstruction
DE102014108789A1 (de) * 2014-06-24 2016-01-07 Byk-Gardner Gmbh Mehrstufiges Verfahren zur Untersuchung von Oberflächen sowie entsprechende Vorrichtung
US10838092B2 (en) * 2014-07-24 2020-11-17 Exxonmobil Upstream Research Company Estimating multiple subsurface parameters by cascaded inversion of wavefield components
US10073183B2 (en) 2014-10-20 2018-09-11 Pgs Geophysical As Methods and systems that attenuate noise in seismic data
US9784867B2 (en) * 2015-04-01 2017-10-10 Schlumberger Technology Corporation Seismic data processing
MX2017014013A (es) * 2015-05-01 2018-03-23 Schlumberger Technology Bv Prospección con fuentes directivas de vibradores marinos.
EP3292428A4 (en) 2015-05-05 2019-06-12 Services Petroliers Schlumberger Removal of acquisition effects from marine seismic data
CN108431636B (zh) 2015-12-02 2021-02-05 斯伦贝谢技术有限公司 平均至少相隔二十米且成对排列的陆地地震传感器与相邻多分量地震传感器
GB2547940A (en) * 2016-03-04 2017-09-06 Robertsson Johan Source separation method
WO2017218723A1 (en) 2016-06-15 2017-12-21 Schlumberger Technology Corporation Systems and methods for acquiring seismic data with gradient data
GB2565744A (en) * 2016-08-19 2019-02-20 Halliburton Energy Services Inc Full waveform inversion of vertical seismic profile data for anisotropic velocities using pseudo-acoustic wave equations
GB2558630A (en) * 2017-01-12 2018-07-18 Seismic Apparition Gmbh Method for dealiasing data
US20190011583A1 (en) * 2017-07-06 2019-01-10 Chevron U.S.A. Inc. System and method for full waveform inversion of seismic data
WO2019017957A1 (en) * 2017-07-20 2019-01-24 Halliburton Energy Services, Inc. MODELING DIPLES FOR ELECTRIC AND / OR MAGNETIC FIELDS
US10788597B2 (en) * 2017-12-11 2020-09-29 Saudi Arabian Oil Company Generating a reflectivity model of subsurface structures
EP3746815A1 (en) * 2018-02-02 2020-12-09 Fairfield Industries, Inc. Seismic imaging with a temporal decomposition imaging condition
US11656377B2 (en) * 2018-03-30 2023-05-23 Cgg Services Sas Visco-acoustic full waveform inversion of velocity and Q
US11231516B2 (en) * 2018-05-15 2022-01-25 Exxonmobil Upstream Research Company Direct migration of simultaneous-source survey data
WO2019224623A1 (en) 2018-05-24 2019-11-28 King Abdullah University Of Science And Technology Method for partial differential equation inversion of data
US11372123B2 (en) 2019-10-07 2022-06-28 Exxonmobil Upstream Research Company Method for determining convergence in full wavefield inversion of 4D seismic data
US11914101B2 (en) 2020-01-31 2024-02-27 ExxonMobil Technology and Engineering Company Method for partitioning a search direction when using least squares reverse time migration
US11320557B2 (en) 2020-03-30 2022-05-03 Saudi Arabian Oil Company Post-stack time domain image with broadened spectrum
US11614555B2 (en) 2020-09-14 2023-03-28 China Petroleum & Chemical Corporation Method and system for connecting elements to sources and receivers during spectrum element method and finite element method seismic wave modeling
CN112415579A (zh) * 2020-11-03 2021-02-26 中国石油天然气集团有限公司 同时源随机激发方法、系统及装置
US20230194736A1 (en) * 2021-12-16 2023-06-22 Chevron U.S.A. Inc. System and method for robust seismic imaging
CN114296134B (zh) * 2021-12-24 2023-03-31 西安交通大学 一种深度卷积网络地震资料解混方法及系统
CN115060769B (zh) * 2022-06-07 2024-04-02 深圳大学 一种基于智能反演的隧道围岩裂隙及松动检测方法、系统

Family Cites Families (205)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3812457A (en) 1969-11-17 1974-05-21 Shell Oil Co Seismic exploration method
US3864667A (en) 1970-09-11 1975-02-04 Continental Oil Co Apparatus for surface wave parameter determination
US3984805A (en) 1973-10-18 1976-10-05 Daniel Silverman Parallel operation of seismic vibrators without phase control
US4168485A (en) 1974-08-12 1979-09-18 Continental Oil Company Simultaneous use of pseudo-random control signals in vibrational exploration methods
US4545039A (en) 1982-09-09 1985-10-01 Western Geophysical Co. Of America Methods for seismic exploration
US4675851A (en) 1982-09-09 1987-06-23 Western Geophysical Co. Method for seismic exploration
US4575830A (en) 1982-10-15 1986-03-11 Schlumberger Technology Corporation Indirect shearwave determination
US4594662A (en) 1982-11-12 1986-06-10 Schlumberger Technology Corporation Diffraction tomography systems and methods with fixed detector arrays
US4562540A (en) 1982-11-12 1985-12-31 Schlumberger Technology Corporation Diffraction tomography system and methods
FR2543306B1 (fr) 1983-03-23 1985-07-26 Elf Aquitaine Procede et dispositif pour l'optimisation des donnees sismiques
US4924390A (en) 1985-03-04 1990-05-08 Conoco, Inc. Method for determination of earth stratum elastic parameters using seismic energy
US4715020A (en) 1986-10-29 1987-12-22 Western Atlas International, Inc. Simultaneous performance of multiple seismic vibratory surveys
FR2589587B1 (fr) 1985-10-30 1988-02-05 Inst Francais Du Petrole Procede de prospection sismique marine utilisant un signal vibratoire code et dispositif pour sa mise en oeuvre
US4707812A (en) 1985-12-09 1987-11-17 Atlantic Richfield Company Method of suppressing vibration seismic signal correlation noise
US4823326A (en) 1986-07-21 1989-04-18 The Standard Oil Company Seismic data acquisition technique having superposed signals
US4686654A (en) 1986-07-31 1987-08-11 Western Geophysical Company Of America Method for generating orthogonal sweep signals
US4766574A (en) 1987-03-31 1988-08-23 Amoco Corporation Method for depth imaging multicomponent seismic data
US4953657A (en) 1987-11-30 1990-09-04 Halliburton Geophysical Services, Inc. Time delay source coding
US4969129A (en) 1989-09-20 1990-11-06 Texaco Inc. Coding seismic sources
US4982374A (en) 1989-10-23 1991-01-01 Halliburton Geophysical Services, Inc. Method of source coding and harmonic cancellation for vibrational geophysical survey sources
GB9011836D0 (en) 1990-05-25 1990-07-18 Mason Iain M Seismic surveying
US6005916A (en) 1992-10-14 1999-12-21 Techniscan, Inc. Apparatus and method for imaging with wavefields using inverse scattering techniques
US5469062A (en) 1994-03-11 1995-11-21 Baker Hughes, Inc. Multiple depths and frequencies for simultaneous inversion of electromagnetic borehole measurements
GB2322704B (en) 1994-07-07 1998-12-09 Geco As Method of Processing seismic data
US5583825A (en) 1994-09-02 1996-12-10 Exxon Production Research Company Method for deriving reservoir lithology and fluid content from pre-stack inversion of seismic data
US5924049A (en) 1995-04-18 1999-07-13 Western Atlas International, Inc. Methods for acquiring and processing seismic data
EP0766836B1 (en) 1995-04-18 2003-01-29 Western Atlas International, Inc. Uniform subsurface coverage at steep dips
US5719821A (en) 1995-09-29 1998-02-17 Atlantic Richfield Company Method and apparatus for source separation of seismic vibratory signals
US5721710A (en) 1995-09-29 1998-02-24 Atlantic Richfield Company High fidelity vibratory source seismic method with source separation
US5822269A (en) 1995-11-13 1998-10-13 Mobil Oil Corporation Method for separation of a plurality of vibratory seismic energy source signals
US5715213A (en) 1995-11-13 1998-02-03 Mobil Oil Corporation High fidelity vibratory source seismic method using a plurality of vibrator sources
US5790473A (en) 1995-11-13 1998-08-04 Mobil Oil Corporation High fidelity vibratory source seismic method for use in vertical seismic profile data gathering with a plurality of vibratory seismic energy sources
US5838634A (en) 1996-04-04 1998-11-17 Exxon Production Research Company Method of generating 3-D geologic models incorporating geologic and geophysical constraints
US5798982A (en) 1996-04-29 1998-08-25 The Trustees Of Columbia University In The City Of New York Method for inverting reflection trace data from 3-D and 4-D seismic surveys and identifying subsurface fluid and pathways in and among hydrocarbon reservoirs based on impedance models
GB9612471D0 (en) 1996-06-14 1996-08-14 Geco As Method and apparatus for multiple seismic vibratory surveys
US5878372A (en) 1997-03-04 1999-03-02 Western Atlas International, Inc. Method for simultaneous inversion processing of well log data using a plurality of earth models
US5999489A (en) 1997-03-21 1999-12-07 Tomoseis Inc. High vertical resolution crosswell seismic imaging
US6014342A (en) 1997-03-21 2000-01-11 Tomo Seis, Inc. Method of evaluating a subsurface region using gather sensitive data discrimination
US5920828A (en) 1997-06-02 1999-07-06 Baker Hughes Incorporated Quality control seismic data processing system
FR2765692B1 (fr) 1997-07-04 1999-09-10 Inst Francais Du Petrole Methode pour modeliser en 3d l'impedance d'un milieu heterogene
GB2329043B (en) 1997-09-05 2000-04-26 Geco As Method of determining the response caused by model alterations in seismic simulations
US5999488A (en) 1998-04-27 1999-12-07 Phillips Petroleum Company Method and apparatus for migration by finite differences
US6219621B1 (en) 1998-06-30 2001-04-17 Exxonmobil Upstream Research Co. Sparse hyperbolic inversion of seismic data
US6388947B1 (en) 1998-09-14 2002-05-14 Tomoseis, Inc. Multi-crosswell profile 3D imaging and method
FR2784195B1 (fr) 1998-10-01 2000-11-17 Inst Francais Du Petrole Methode pour realiser en 3d avant sommation, une migration de donnees sismiques
US6574564B2 (en) 1998-10-01 2003-06-03 Institut Francais Du Petrole 3D prestack seismic data migration method
US6225803B1 (en) 1998-10-29 2001-05-01 Baker Hughes Incorporated NMR log processing using wavelet filter and iterative inversion
US6021094A (en) 1998-12-03 2000-02-01 Sandia Corporation Method of migrating seismic records
US6754588B2 (en) 1999-01-29 2004-06-22 Platte River Associates, Inc. Method of predicting three-dimensional stratigraphy using inverse optimization techniques
WO2000048022A1 (en) 1999-02-12 2000-08-17 Schlumberger Limited Uncertainty constrained subsurface modeling
US6058073A (en) 1999-03-30 2000-05-02 Atlantic Richfield Company Elastic impedance estimation for inversion of far offset seismic sections
FR2792419B1 (fr) 1999-04-16 2001-09-07 Inst Francais Du Petrole Methode pour obtenir un modele optimal d'une caracteristique physique dans un milieu heterogene, tel que le sous-sol
GB9927395D0 (en) 1999-05-19 2000-01-19 Schlumberger Holdings Improved seismic data acquisition method
US6327537B1 (en) 1999-07-19 2001-12-04 Luc T. Ikelle Multi-shooting approach to seismic modeling and acquisition
FR2798197B1 (fr) 1999-09-02 2001-10-05 Inst Francais Du Petrole Methode pour former un modele d'une formation geologique, contraint par des donnees dynamiques et statiques
EP2296013B1 (en) 1999-10-22 2016-03-30 CGG Services (NL) B.V. Method of estimating elastic and compositional parameters from seismic and echo-acoustic data
US6480790B1 (en) 1999-10-29 2002-11-12 Exxonmobil Upstream Research Company Process for constructing three-dimensional geologic models having adjustable geologic interfaces
FR2800473B1 (fr) 1999-10-29 2001-11-30 Inst Francais Du Petrole Methode pour modeliser en 2d ou 3d un milieu heterogene tel que le sous-sol decrit par plusieurs parametres physiques
US6836448B2 (en) 2000-01-21 2004-12-28 Schlumberger Technology Corporation System and method for seismic wavefield separation
CA2397664C (en) 2000-01-21 2008-11-18 Schlumberger Canada Limited System and method for estimating seismic material properties
US6826486B1 (en) 2000-02-11 2004-11-30 Schlumberger Technology Corporation Methods and apparatus for predicting pore and fracture pressures of a subsurface formation
FR2805051B1 (fr) 2000-02-14 2002-12-06 Geophysique Cie Gle Methode de surveillance sismique d'une zone souterraine par utilisation simultanee de plusieurs sources vibrosismiques
GB2359363B (en) 2000-02-15 2002-04-03 Geco Prakla Processing simultaneous vibratory seismic data
US6687659B1 (en) 2000-03-24 2004-02-03 Conocophillips Company Method and apparatus for absorbing boundary conditions in numerical finite-difference acoustic applications
US6317695B1 (en) 2000-03-30 2001-11-13 Nutec Sciences, Inc. Seismic data processing method
AU1176802A (en) 2000-10-17 2002-04-29 Westerngeco Llc Method of using cascaded sweeps for source coding and harmonic cancellation
WO2002047011A1 (en) 2000-12-08 2002-06-13 Ortoleva Peter J Methods for modeling multi-dimensional domains using information theory to resolve gaps in data and in theories
FR2818753B1 (fr) 2000-12-21 2003-03-21 Inst Francais Du Petrole Methode et dispositif de prospection sismique par emission simultanee de signaux sismisques obtenus en codant un signal par des sequences pseudo aleatoires
FR2821677B1 (fr) 2001-03-05 2004-04-30 Geophysique Cie Gle Perfectionnements aux procedes d'inversion tomographique d'evenements pointes sur les donnees sismiques migrees
US6751558B2 (en) 2001-03-13 2004-06-15 Conoco Inc. Method and process for prediction of subsurface fluid and rock pressures in the earth
US6927698B2 (en) 2001-08-27 2005-08-09 Larry G. Stolarczyk Shuttle-in receiver for radio-imaging underground geologic structures
US6545944B2 (en) 2001-05-30 2003-04-08 Westerngeco L.L.C. Method for acquiring and processing of data from two or more simultaneously fired sources
US6882958B2 (en) 2001-06-28 2005-04-19 National Instruments Corporation System and method for curve fitting using randomized techniques
GB2379013B (en) 2001-08-07 2005-04-20 Abb Offshore Systems Ltd Microseismic signal processing
US6593746B2 (en) 2001-08-27 2003-07-15 Larry G. Stolarczyk Method and system for radio-imaging underground geologic structures
US7672824B2 (en) 2001-12-10 2010-03-02 Westerngeco L.L.C. Method for shallow water flow detection
US7069149B2 (en) 2001-12-14 2006-06-27 Chevron U.S.A. Inc. Process for interpreting faults from a fault-enhanced 3-dimensional seismic attribute volume
US7330799B2 (en) 2001-12-21 2008-02-12 Société de commercialisation des produits de la recherche appliquée-Socpra Sciences et Génie s.e.c. Method and algorithm for using surface waves
US6842701B2 (en) 2002-02-25 2005-01-11 Westerngeco L.L.C. Method of noise removal for cascaded sweep data
GB2387226C (en) 2002-04-06 2008-05-12 Westerngeco Ltd A method of seismic surveying
FR2839368B1 (fr) 2002-05-06 2004-10-01 Total Fina Elf S A Methode de decimation de traces sismiques pilotee par le trajet sismique
US6832159B2 (en) 2002-07-11 2004-12-14 Schlumberger Technology Corporation Intelligent diagnosis of environmental influence on well logs with model-based inversion
US6906981B2 (en) 2002-07-17 2005-06-14 Pgs Americas, Inc. Method and system for acquiring marine seismic data using multiple seismic sources
FR2843202B1 (fr) 2002-08-05 2004-09-10 Inst Francais Du Petrole Methode pour former un modele representatif de la distribution d'une grandeur physique dans une zone souterraine, affranchi de l'effet de bruits correles entachant des donnees d'exploration
WO2004034088A2 (en) 2002-10-04 2004-04-22 Paradigm Geophysical Corporation Method and system for limited frequency seismic imaging
GB2396448B (en) 2002-12-21 2005-03-02 Schlumberger Holdings System and method for representing and processing and modeling subterranean surfaces
US20040225483A1 (en) 2003-02-24 2004-11-11 Michal Okoniewski Fdtd hardware acceleration system
US6735527B1 (en) 2003-02-26 2004-05-11 Landmark Graphics Corporation 3-D prestack/poststack multiple prediction
US6999880B2 (en) 2003-03-18 2006-02-14 The Regents Of The University Of California Source-independent full waveform inversion of seismic data
US7184367B2 (en) 2003-03-27 2007-02-27 Exxonmobil Upstream Research Company Method to convert seismic traces into petrophysical property logs
ATE535826T1 (de) 2003-04-01 2011-12-15 Exxonmobil Upstream Res Co Geformte hochfrequenz-vibrationsquelle
US7072767B2 (en) 2003-04-01 2006-07-04 Conocophillips Company Simultaneous inversion for source wavelet and AVO parameters from prestack seismic data
NO322089B1 (no) 2003-04-09 2006-08-14 Norsar V Daglig Leder Fremgangsmate for simulering av lokale prestakk dypmigrerte seismiske bilder
GB2400438B (en) 2003-04-11 2005-06-01 Westerngeco Ltd Determination of waveguide parameters
US6970397B2 (en) 2003-07-09 2005-11-29 Gas Technology Institute Determination of fluid properties of earth formations using stochastic inversion
US6882938B2 (en) 2003-07-30 2005-04-19 Pgs Americas, Inc. Method for separating seismic signals from two or more distinct sources
GB2405473B (en) 2003-08-23 2005-10-05 Westerngeco Ltd Multiple attenuation method
US6944546B2 (en) 2003-10-01 2005-09-13 Halliburton Energy Services, Inc. Method and apparatus for inversion processing of well logging data in a selected pattern space
US6901333B2 (en) 2003-10-27 2005-05-31 Fugro N.V. Method and device for the generation and application of anisotropic elastic parameters
US7046581B2 (en) 2003-12-01 2006-05-16 Shell Oil Company Well-to-well tomography
US20050128874A1 (en) 2003-12-15 2005-06-16 Chevron U.S.A. Inc. Methods for acquiring and processing seismic data from quasi-simultaneously activated translating energy sources
US7359283B2 (en) 2004-03-03 2008-04-15 Pgs Americas, Inc. System for combining signals of pressure sensors and particle motion sensors in marine seismic streamers
US7791980B2 (en) 2004-05-21 2010-09-07 Westerngeco L.L.C. Interpolation and extrapolation method for seismic recordings
FR2872584B1 (fr) 2004-06-30 2006-08-11 Inst Francais Du Petrole Methode pour simuler le depot sedimentaire dans un bassin respectant les epaisseurs des sequences sedimentaires
EP1617309B1 (en) 2004-07-15 2011-01-12 Fujitsu Limited Simulation technique with local grid refinement
US7646924B2 (en) 2004-08-09 2010-01-12 David Leigh Donoho Method and apparatus for compressed sensing
US7480206B2 (en) 2004-09-13 2009-01-20 Chevron U.S.A. Inc. Methods for earth modeling and seismic imaging using interactive and selective updating
FR2876458B1 (fr) 2004-10-08 2007-01-19 Geophysique Cie Gle Perfectionnement aux traitements sismiques pour la suppression des reflexions multiples
GB2422433B (en) 2004-12-21 2008-03-19 Sondex Wireline Ltd Method and apparatus for determining the permeability of earth formations
US7373251B2 (en) 2004-12-22 2008-05-13 Marathon Oil Company Method for predicting quantitative values of a rock or fluid property in a reservoir using seismic data
US7230879B2 (en) 2005-02-12 2007-06-12 Chevron U.S.A. Inc. Method and apparatus for true relative amplitude correction of seismic data for normal moveout stretch effects
EP1859301B1 (en) 2005-02-22 2013-07-17 Paradigm Geophysical Ltd. Multiple suppression in angle domain time and depth migration
US7840625B2 (en) 2005-04-07 2010-11-23 California Institute Of Technology Methods for performing fast discrete curvelet transforms of data
US7271747B2 (en) 2005-05-10 2007-09-18 Rice University Method and apparatus for distributed compressed sensing
US7405997B2 (en) 2005-08-11 2008-07-29 Conocophillips Company Method of accounting for wavelet stretch in seismic data
EP1941386A4 (en) 2005-10-18 2010-03-17 Sinvent As IMAGING OF GEOLOGICAL RESPONSES DATA WITH FLOW PROCESSORS
US7373252B2 (en) 2005-11-04 2008-05-13 Western Geco L.L.C. 3D pre-stack full waveform inversion
FR2895091B1 (fr) 2005-12-21 2008-02-22 Inst Francais Du Petrole Methode pour mettre a jour un modele geologique par des donnees sismiques
US7400552B2 (en) 2006-01-19 2008-07-15 Westerngeco L.L.C. Methods and systems for efficiently acquiring towed streamer seismic surveys
GB2436626B (en) 2006-03-28 2008-08-06 Westerngeco Seismic Holdings Method of evaluating the interaction between a wavefield and a solid body
US7620534B2 (en) 2006-04-28 2009-11-17 Saudi Aramco Sound enabling computerized system for real time reservoir model calibration using field surveillance data
US20070274155A1 (en) 2006-05-25 2007-11-29 Ikelle Luc T Coding and Decoding: Seismic Data Modeling, Acquisition and Processing
US7725266B2 (en) 2006-05-31 2010-05-25 Bp Corporation North America Inc. System and method for 3D frequency domain waveform inversion based on 3D time-domain forward modeling
US7599798B2 (en) 2006-09-11 2009-10-06 Westerngeco L.L.C. Migrating composite seismic response data to produce a representation of a seismic volume
WO2008042081A1 (en) 2006-09-28 2008-04-10 Exxonmobil Upstream Research Company Iterative inversion of data from simultaneous geophysical sources
EP2089738B1 (en) 2006-12-07 2019-01-02 Council of Scientific & Industrial Research A method for computing an exact impulse response of a plane acoustic reflector at zero offset due to a point acoustic source
ATE543109T1 (de) 2007-01-20 2012-02-15 Spectraseis Ag Zeitumkehr-reservoir-lokalisierung
US8248886B2 (en) 2007-04-10 2012-08-21 Exxonmobil Upstream Research Company Separation and noise removal for multiple vibratory source seismic data
US7715986B2 (en) 2007-05-22 2010-05-11 Chevron U.S.A. Inc. Method for identifying and removing multiples for imaging with beams
US7974824B2 (en) 2007-06-29 2011-07-05 Westerngeco L. L. C. Seismic inversion of data containing surface-related multiples
JP2009063942A (ja) 2007-09-10 2009-03-26 Sumitomo Electric Ind Ltd 遠赤外線カメラ用レンズ、レンズユニット及び撮像装置
US20090070042A1 (en) 2007-09-11 2009-03-12 Richard Birchwood Joint inversion of borehole acoustic radial profiles for in situ stresses as well as third-order nonlinear dynamic moduli, linear dynamic elastic moduli, and static elastic moduli in an isotropically stressed reference state
US20090083006A1 (en) 2007-09-20 2009-03-26 Randall Mackie Methods and apparatus for three-dimensional inversion of electromagnetic data
WO2009067041A1 (en) 2007-11-19 2009-05-28 Steklov Mathematical Institute Ras Method and system for evaluating the characteristic properties of two contacting media and of the interface between them based on mixed surface waves propagating along the interface
US20090164186A1 (en) 2007-12-20 2009-06-25 Bhp Billiton Innovation Pty Ltd. Method for determining improved estimates of properties of a model
CA2707526C (en) 2008-01-08 2016-01-26 Exxonmobil Upstream Research Company Spectral shaping inversion and migration of seismic data
US8577660B2 (en) 2008-01-23 2013-11-05 Schlumberger Technology Corporation Three-dimensional mechanical earth modeling
AU2009226041B2 (en) 2008-03-21 2013-09-05 Exxonmobil Upstream Research Company An efficient method for inversion of geophysical data
EP2105765A1 (en) 2008-03-28 2009-09-30 Schlumberger Holdings Limited Simultaneous inversion of induction data for dielectric permittivity and electric conductivity
WO2009120401A1 (en) 2008-03-28 2009-10-01 Exxonmobil Upstream Research Company Characterizing spatial variablility of surface waves in seismic processing
US8275592B2 (en) 2008-04-07 2012-09-25 Westerngeco L.L.C. Joint inversion of time domain controlled source electromagnetic (TD-CSEM) data and further data
US8494777B2 (en) 2008-04-09 2013-07-23 Schlumberger Technology Corporation Continuous microseismic mapping for real-time 3D event detection and location
US8345510B2 (en) 2008-06-02 2013-01-01 Pgs Geophysical As Method for aquiring and processing marine seismic data to extract and constructively use the up-going and down-going wave-fields emitted by the source(s)
US20110182141A1 (en) 2008-08-14 2011-07-28 Schlumberger Technology Corporation Method and system for monitoring a logging tool position in a borehole
US8559270B2 (en) 2008-08-15 2013-10-15 Bp Corporation North America Inc. Method for separating independent simultaneous sources
US8295124B2 (en) 2008-08-15 2012-10-23 Bp Corporation North America Inc. Method for separating independent simultaneous sources
US20100054082A1 (en) 2008-08-29 2010-03-04 Acceleware Corp. Reverse-time depth migration with reduced memory requirements
US8296069B2 (en) 2008-10-06 2012-10-23 Bp Corporation North America Inc. Pseudo-analytical method for the solution of wave equations
US7616523B1 (en) 2008-10-22 2009-11-10 Pgs Geophysical As Method for combining pressure and motion seismic signals from streamers where sensors are not at a common depth
US9213119B2 (en) 2008-10-29 2015-12-15 Conocophillips Company Marine seismic acquisition
US20100118651A1 (en) 2008-11-10 2010-05-13 Chevron U.S.A. Inc. Method for generation of images related to a subsurface region of interest
US20100142316A1 (en) 2008-12-07 2010-06-10 Henk Keers Using waveform inversion to determine properties of a subsurface medium
US8095345B2 (en) 2009-01-20 2012-01-10 Chevron U.S.A. Inc Stochastic inversion of geophysical data for estimating earth model parameters
US9052410B2 (en) 2009-02-12 2015-06-09 Conocophillips Company Multiple seismic signal inversion
US20110299361A1 (en) 2009-02-17 2011-12-08 Changsoo Shin Apparatus and method for imaging subsurface structure
US8352190B2 (en) 2009-02-20 2013-01-08 Exxonmobil Upstream Research Company Method for analyzing multiple geophysical data sets
US9075163B2 (en) 2009-04-17 2015-07-07 Westerngeco L.L.C. Interferometric seismic data processing
US8176360B2 (en) 2009-08-11 2012-05-08 Texas Memory Systems, Inc. Method and apparatus for addressing actual or predicted failures in a FLASH-based storage system
US20110044127A1 (en) 2009-08-19 2011-02-24 Clement Kostov Removing free-surface effects from seismic data acquired in a towed survey
US8923093B2 (en) 2009-08-25 2014-12-30 Westerngeco L.L.C. Determining the quality of a seismic inversion
US20110131020A1 (en) 2009-09-09 2011-06-02 Conocophillips Company Dip guided full waveform inversion
GB2486121B (en) 2009-10-01 2014-08-13 Halliburton Energy Serv Inc Apparatus and methods of locating downhole anomalies
US9244181B2 (en) 2009-10-19 2016-01-26 Westerngeco L.L.C. Full-waveform inversion in the traveltime domain
GB2490051B (en) 2009-12-07 2015-04-01 Geco Technology Bv Simultaneous joint inversion of surface wave and refraction data
WO2011091216A2 (en) 2010-01-22 2011-07-28 Schlumberger Canada Limited Real-time formation anisotropy and dip evaluation using tri-axial induction measurements
EP2529254A4 (en) 2010-01-25 2017-08-02 CGG Veritas Services (U.S.) Inc. Methods and systems for estimating stress using seismic data
US20130098608A1 (en) 2010-01-29 2013-04-25 Robert Barnum Temporary field storage of gas to optimize field development
US8265875B2 (en) 2010-01-29 2012-09-11 Westerngeco L.L.C. Interpolation of periodic data
US8537638B2 (en) 2010-02-10 2013-09-17 Exxonmobil Upstream Research Company Methods for subsurface parameter estimation in full wavefield inversion and reverse-time migration
WO2011112932A1 (en) 2010-03-12 2011-09-15 Cggveritas Services (Us) Inc. Methods and systems for performing azimuthal simultaneous elatic inversion
US8680865B2 (en) 2010-03-19 2014-03-25 Schlumberger Technology Corporation Single well reservoir imaging apparatus and methods
US20110235464A1 (en) 2010-03-24 2011-09-29 John Brittan Method of imaging the earth's subsurface during marine seismic data acquisition
US8223587B2 (en) 2010-03-29 2012-07-17 Exxonmobil Upstream Research Company Full wavefield inversion using time varying filters
US9176244B2 (en) 2010-03-31 2015-11-03 Schlumberger Technology Corporation Data set inversion using source-receiver compression
KR101167715B1 (ko) 2010-04-30 2012-07-20 서울대학교산학협력단 복소 구배 최소자승법에의한 파형 역산을 이용한 지하 구조의 영상화 장치 및 방법
US8576663B2 (en) 2010-04-30 2013-11-05 Schlumberger Technology Corporation Multicomponent seismic inversion of VSP data
US8694299B2 (en) 2010-05-07 2014-04-08 Exxonmobil Upstream Research Company Artifact reduction in iterative inversion of geophysical data
US8756042B2 (en) 2010-05-19 2014-06-17 Exxonmobile Upstream Research Company Method and system for checkpointing during simulations
WO2011158652A1 (ja) 2010-06-15 2011-12-22 電気化学工業株式会社 透光性硬質基板積層体の製造方法
US20110320180A1 (en) 2010-06-29 2011-12-29 Al-Saleh Saleh M Migration Velocity Analysis of Seismic Data Using Common Image Cube and Green's Functions
US8612188B2 (en) 2010-07-12 2013-12-17 The University Of Manchester Wave modelling
EA026517B1 (ru) * 2010-08-02 2017-04-28 Бп Корпорейшн Норт Америка Инк. Способ сейсмической разведки
US9195783B2 (en) 2010-08-16 2015-11-24 Exxonmobil Upstream Research Company Reducing the dimensionality of the joint inversion problem
US20120051176A1 (en) 2010-08-31 2012-03-01 Chevron U.S.A. Inc. Reverse time migration back-scattering noise removal using decomposed wavefield directivity
US8437998B2 (en) 2010-09-27 2013-05-07 Exxonmobil Upstream Research Company Hybrid method for full waveform inversion using simultaneous and sequential source method
CN103119552B (zh) 2010-09-27 2016-06-08 埃克森美孚上游研究公司 同时源编码和源分离作为全波场反演的实际解决方案
BR112013002847A2 (pt) 2010-09-27 2016-06-07 Exxonmobil Upstream Res Co método híbrido para inversão de forma de onda completa que usa um método de fonte simultânea e sequencial
AU2011310635B2 (en) 2010-09-28 2014-09-18 Shell Internationale Research Maatschappij B.V. Earth model estimation through an acoustic Full Waveform Inversion of seismic data
KR101797451B1 (ko) 2010-12-01 2017-11-14 엑손모빌 업스트림 리서치 캄파니 상호상관 목적 함수를 통한 해양 스트리머 데이터에 대한 동시 소스 반전
US9134442B2 (en) 2010-12-16 2015-09-15 Bp Corporation North America Inc. Seismic acquisition using narrowband seismic sources
IT1404170B1 (it) * 2011-02-10 2013-11-15 Eni Spa Metodo di indagine sismica del sottosuolo
US9702994B2 (en) 2011-02-18 2017-07-11 Westerngeco L.L.C. Waveform inversion by multiple shot-encoding for non-fixed spread geometries
KR101931488B1 (ko) 2011-03-30 2018-12-24 엑손모빌 업스트림 리서치 캄파니 스펙트럼 성형을 이용하는 전 파동장 반전의 수렴 레이트
US20120275267A1 (en) 2011-04-26 2012-11-01 Ramesh Neelamani Seismic Data Processing
CA2833968C (en) 2011-05-13 2017-12-12 Saudi Arabian Oil Company Coupled time-distance dependent swept frequency source acquisition design and data de-noising
US20120316791A1 (en) 2011-06-08 2012-12-13 Chevron U.S.A. Inc. System and method for seismic data inversion by non-linear model update
US20120316844A1 (en) 2011-06-08 2012-12-13 Chevron U.S.A. Inc. System and method for data inversion with phase unwrapping
US20120316790A1 (en) 2011-06-08 2012-12-13 Chevron U.S.A. Inc. System and method for data inversion with phase extrapolation
US9075159B2 (en) 2011-06-08 2015-07-07 Chevron U.S.A., Inc. System and method for seismic data inversion
US9075162B2 (en) * 2011-11-10 2015-07-07 Pgs Geophysical As Method and system for separating seismic sources in marine simultaneous shooting acquisition
US9091788B2 (en) 2011-11-14 2015-07-28 Cggveritas Services Sa Device and method for de-blending simultaneous shooting data with apex shifted radon transform
US9103943B2 (en) 2011-11-28 2015-08-11 Fugro-Geoteam As Acquisition and processing of multi-source broadband marine seismic data
US9176930B2 (en) 2011-11-29 2015-11-03 Exxonmobil Upstream Research Company Methods for approximating hessian times vector operation in full wavefield inversion
AU2013230789B2 (en) 2012-03-08 2016-02-11 Exxonmobil Upstream Research Company Orthogonal source and receiver encoding
US20130311149A1 (en) 2012-05-17 2013-11-21 Yaxun Tang Tomographically Enhanced Full Wavefield Inversion

Also Published As

Publication number Publication date
MX359179B (es) 2018-09-18
EP3351972A1 (en) 2018-07-25
AU2014309376B2 (en) 2016-11-17
CA2913496A1 (en) 2015-02-26
US9772413B2 (en) 2017-09-26
EP3036566A2 (en) 2016-06-29
WO2015026451A3 (en) 2015-04-30
US20150057938A1 (en) 2015-02-26
WO2015026451A2 (en) 2015-02-26
DK3036566T3 (en) 2018-07-23
AU2014309376A1 (en) 2016-03-03
BR112015030104A2 (pt) 2017-07-25
EP3036566B1 (en) 2018-04-04
CA2913496C (en) 2018-08-14

Similar Documents

Publication Publication Date Title
ES2671789T3 (es) Búsqueda simultánea de fuentes durante la adquisición sísmica y la inversión sísmica
US11774618B2 (en) Marine vibrator directive source survey
US10520623B2 (en) Methods and systems for marine survey acquisition
EP3625593B1 (en) Non-uniform optimal survey design principles
EP2972512B1 (en) Seismic acquisition using phase-shifted sweeps
US9310503B2 (en) Methods to process seismic data contaminated by coherent energy radiated from more than one source
ES2968710T3 (es) Imágenes sísmicas pasivas
Peng et al. Deblending of simulated simultaneous sources using an iterative approach-an experiment with variable-depth streamer data
US9453928B2 (en) Methods and computing systems for processing data
RU2570827C2 (ru) Гибридный способ для полноволновой инверсии с использованием способа одновременных и последовательных источников
Ziramov et al. Kevitsa Ni-Cu-PGE deposit, North Finland? A seismic case study
Lim et al. High resolution anisotropic Earth model building on conventional seismic data using full waveform inversion: a case study offshore Australia
Tarasov et al. Modern techniques in crosshole seismic testing (Cst) fulfilment: Case study
Lim et al. High Resolution Anisotropic Earth Model Building on Conventional Seismic Data Using Full Waveform Inversion: A Case Study Offshore Australia
Plessix et al. A land example of full waveform inversion and distance separated simultaneous sweeping
Sava et al. Wavefield tomography without low frequency data