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 PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 claims abstract description 77
- 238000005259 measurement Methods 0.000 claims abstract description 37
- 230000008569 process Effects 0.000 claims abstract description 15
- 230000004913 activation Effects 0.000 claims abstract 11
- 230000000704 physical effect Effects 0.000 claims abstract 8
- 230000006870 function Effects 0.000 claims description 29
- 238000010276 construction Methods 0.000 claims description 9
- 238000006073 displacement reaction Methods 0.000 claims description 5
- 238000005096 rolling process Methods 0.000 claims description 3
- 238000001994 activation Methods 0.000 claims 8
- 238000000926 separation method Methods 0.000 description 24
- 238000005474 detonation Methods 0.000 description 22
- 238000003384 imaging method Methods 0.000 description 10
- 238000012545 processing Methods 0.000 description 9
- 238000005094 computer simulation Methods 0.000 description 8
- 238000004088 simulation Methods 0.000 description 7
- 230000001934 delay Effects 0.000 description 6
- 230000004044 response Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000001914 filtration Methods 0.000 description 5
- 230000005012 migration Effects 0.000 description 5
- 238000013508 migration Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000010304 firing Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000002310 reflectometry Methods 0.000 description 3
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 238000005314 correlation function Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 206010010904 Convulsion Diseases 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013079 data visualisation Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 238000007654 immersion Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/003—Seismic data acquisition in general, e.g. survey design
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/003—Seismic data acquisition in general, e.g. survey design
- G01V1/005—Seismic 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/24—Recording seismic data
- G01V1/247—Digital recording of seismic data, e.g. in acquisition units or nodes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
- G01V1/368—Inverse filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
- G01V1/3808—Seismic data acquisition, e.g. survey design
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/12—Signal generation
- G01V2210/121—Active source
- G01V2210/1212—Shot
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/12—Signal generation
- G01V2210/127—Cooperating multiple sources
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/14—Signal detection
- G01V2210/142—Receiver location
- G01V2210/1423—Sea
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6161—Seismic or acoustic, e.g. land or sea measurements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface 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).
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)
- 5101520253035404550556065REIVINDICACIONES1. 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; yprocesar 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. 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. 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. 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. 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. 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. 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. 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. 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.510152025303540455055
- 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; ouna 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. 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. 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. 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. 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. 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; ouna 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.
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)
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)
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 |
-
2014
- 2014-07-11 US US14/329,431 patent/US9772413B2/en active Active
- 2014-07-11 EP EP14744730.4A patent/EP3036566B1/en not_active Not-in-force
- 2014-07-11 ES ES14744730.4T patent/ES2671789T3/es active Active
- 2014-07-11 MX MX2015016629A patent/MX359179B/es active IP Right Grant
- 2014-07-11 BR BR112015030104A patent/BR112015030104A2/pt not_active IP Right Cessation
- 2014-07-11 WO PCT/US2014/046361 patent/WO2015026451A2/en active Application Filing
- 2014-07-11 AU AU2014309376A patent/AU2014309376B2/en not_active Ceased
- 2014-07-11 EP EP18154578.1A patent/EP3351972A1/en not_active Withdrawn
- 2014-07-11 DK DK14744730.4T patent/DK3036566T3/en active
- 2014-07-11 CA CA2913496A patent/CA2913496C/en not_active Expired - Fee Related
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 |