MX2011008294A - Formacion de imagenes con mediciones vectoriales. - Google Patents
Formacion de imagenes con mediciones vectoriales.Info
- Publication number
- MX2011008294A MX2011008294A MX2011008294A MX2011008294A MX2011008294A MX 2011008294 A MX2011008294 A MX 2011008294A MX 2011008294 A MX2011008294 A MX 2011008294A MX 2011008294 A MX2011008294 A MX 2011008294A MX 2011008294 A MX2011008294 A MX 2011008294A
- Authority
- MX
- Mexico
- Prior art keywords
- pressure
- measurements
- image
- seismic
- pressure gradient
- Prior art date
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 29
- 238000003384 imaging method Methods 0.000 title description 5
- 238000000034 method Methods 0.000 claims abstract description 56
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 33
- 238000009530 blood pressure measurement Methods 0.000 claims abstract description 16
- 239000002245 particle Substances 0.000 claims description 16
- 238000012545 processing Methods 0.000 claims description 10
- 230000005012 migration Effects 0.000 claims description 8
- 238000013508 migration Methods 0.000 claims description 8
- 230000008030 elimination Effects 0.000 claims description 6
- 238000003379 elimination reaction Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 229910000831 Steel Inorganic materials 0.000 claims 1
- 239000010959 steel Substances 0.000 claims 1
- 230000006870 function Effects 0.000 description 36
- 238000005755 formation reaction Methods 0.000 description 24
- 238000010586 diagram Methods 0.000 description 4
- 238000004613 tight binding model Methods 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 229930195733 hydrocarbon Natural products 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000001364 causal effect Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 125000001183 hydrocarbyl group Chemical group 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 229910052704 radon Inorganic materials 0.000 description 2
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 2
- 230000001174 ascending effect Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000003643 water by type Substances 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/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Oceanography (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measuring Fluid Pressure (AREA)
Abstract
Una técnica incluye recibir datos sísmicos, los cuales son indicativos de mediciones de presión y mediciones de gradiente de presión adquiridas en un estudio sísmico de al menos una formación subterránea. La técnica incluye modelar una imagen de la(s) formación(es) subterránea(s) en función de las mediciones de presión y de las mediciones de gradiente de presión. La técnica incluye determinar la imagen con base en la modelación.
Description
FORMACIÓN DE IMÁGENES CON MEDICIONES VECTORIALES"
Campo de la Invención
La invención se refiere en términos generales a la formación de imágenes con mediciones vectoriales.
Antecedentes de la Invención
La exploración sísmica implica hacer estudios de formaciones geológicas subterráneas para yacimientos de hidrocarburos. Un estudio típicamente implica implementar fuentes sísmicas y sensores sísmicos en determinadas posiciones. Las fuentes generan ondas sísmicas, las cuales se propagan en las formaciones geológicas que crean cambios de presión y vibraciones en su trayecto. Los cambios en las propiedades elásticas de la formación geológica dispersan las ondas sísmicas, cambiando su dirección de propagación y otras propiedades. Parte de la energía emitida por las fuentes alcanza los sensores sísmicos. Algunos sensores sísmicos son sensibles a los cambios de presión (hidrófonos), otros al movimiento de partículas (por ejemplo, geófonos y/o acelerómetros), y estudios industriales pueden implementar únicamente un tipo de sensores o ambos. En respuesta a los eventos sísmicos detectados, los sensores generan señales eléctricas para producir datos sísmicos. El análisis de los datos sísmicos puede indicar después la presencia o ausencia de posiciones probables de yacimientos de hidrocarburos.
Algunos estudios son conocidos como estudios "marinos" porque
se realizan en ambientes marinos. Sin embargo, los estudios "marinos" pueden realizarse no solamente en ambientes de agua salada, sino también en aguas dulces y salobres. En un tipo de estudio marino, llamado estudio de "barrido reticular", un arreglo de cables marinos y fuentes que contienen sensores sísmicos es remolcado detrás de una nave de estudios.
Breve Descripción de la Invención
En una modalidad de la invención, una técnica incluye la recepción de datos sísmicos, los cuales son indicativos de mediciones de presión y mediciones de gradiente de presión adquiridas en un estudio sísmico de al menos una formación subterránea. La técnica incluye modelar una imagen de la(s) formación(es) subterránea(s) en función de las mediciones de presión y las mediciones de gradiente de presión. La técnica incluye determinar la imagen con base en el modelado.
En otra modalidad de la invención, un sistema incluye una interfaz y un procesador. La interfaz recibe datos sísmicos, los cuales son indicativos de mediciones de presión y mediciones de gradiente de presión adquiridas en un estudio sísmico de al menos una formación subterránea. El procesador procesa los datos sísmicos para determinar una imagen de la(s) formación(es) subterránea(s) con base en un modelo de la imagen en función de las mediciones de presión y las mediciones del gradiente de presión.
Las ventajas y otras características de la invención se volverán
aparentes a partir de los siguientes dibujos, descripción y reivindicaciones.
Breve Descripción de las Figuras
La Figura 1 es un diagrama esquemático de un sistema de adquisición de datos sísmicos marinos de acuerdo con una modalidad de la invención.
La Figura 2 es un diagrama de flujo que representa gráficamente una técnica para procesar datos sísmicos de acuerdo con una modalidad de la invención.
Las Figuras 3 y 4 son diagramas de flujo que representan gráficamente técnicas para generar una imagen de al menos una formación subterránea de acuerdo con modalidades de la invención.
La Figura 5 es un diagrama esquemático de un sistema de procesamiento de datos sísmicos de acuerdo con una modalidad de la invención.
Descripción Detallada de la Invención
La Figura 1 representa gráficamente una modalidad 10 de un sistema de adquisición de datos sísmicos marinos de acuerdo con algunas modalidades de la invención. En el sistema 10, una nave de estudios 20 remolca uno o más cables marinos sísmicos 30 (un cable marino 30 a manera de ejemplo se representa gráficamente en la Figura 1) detrás de la nave 20. Los cables marinos sísmicos 30 pueden tener varios miles de metros de largo y pueden contener diversos cables de
refuerzo (no se muestran), así como también cableado y/o circuitería (no se muestran) que pueden utilizarse para soportar la comunicación a lo largo de los cables marinos 30. En general, cada cable marino 30 incluye un cable primario en el cual se instalan los sensores sísmicos 58 que registran las señales sísmicas.
De acuerdo con modalidades de la invención, los sensores sísmicos 58 pueden ser únicamente sensores de presión o pueden ser sensores sísmicos de multí-componente. Para el caso de los sensores sísmicos de multi-componente, cada sensor es capaz de detectar un campo ondulatorio de presión y al menos un componente de un movimiento de partículas que se encuentra asociado con señales acústicas que están próximas al sensor sísmico de multi-componente. Los ejemplos de movimientos de partículas incluyen uno o más componentes de desplazamiento de partícula, uno o más componentes (componentes en línea (x), transversal (y) y vertical (z) (ver los ejes 59, por ejemplo)) de una velocidad de partícula y uno o más componentes de una aceleración de partícula.
Dependiendo de la modalidad particular de la invención, el sensor sísmico de multi-componente puede incluir uno o más hidrófonos, geófonos, sensores de desplazamiento de partículas, sensores de velocidad de partícula, acelerómetros, sensores de gradiente de presión, o combinaciones de los mismos.
Por ejemplo, de acuerdo con algunas modalidades de la invención, un sensor sísmico de multi-componente particular puede incluir un hidrófono para medir la presión y tres acelerómetros alineados
ortogonalmente para medir tres componentes ortogonales correspondientes de velocidad y/o aceleración de partícula cercanos al sensor sísmico. Se observa que el sensor sísmico de multi-componente puede implementarse como un solo dispositivo o puede implementarse como una pluralidad de dispositivos, dependiendo de la modalidad particular de la invención . Un sensor sísmico de multi-componente particular también puede incluir sensores de gradiente de presión, los cuales constituyen otro tipo de sensores de movimiento de partículas. Cada sensor de gradiente de presión mide el cambio en el campo ondulatorio de presión en un punto particular con respecto a una dirección particular. Por ejemplo, uno de los sensores de gradiente de presión puede adquirir datos sísmicos indicativos de, en un punto particular, la derivada parcial del campo ondulatorio de presión con respecto a la dirección transversal, y el otro sensor de gradiente de presión puede adquirir, un punto particular, datos sísmicos indicativos de los datos de presión con respecto a la dirección en línea.
El sistema de adquisición de datos sísmicos marinos 10 incluye una fuente sísmica 1 04 que puede formarse a partir de uno o más elementos de fuente sísmica, tales como martillos neumáticos, por ejemplo, los cuales se conectan a la nave - de estudios 20. Alternativamente, en otras modalidades de la invención, la fuente sísmica 1 04 puede operar independientemente de la nave de estudios 20, porque la fuente sísmica 104 ' puede acoplarse a otras naves o boyas, por mencionar unos ejemplos.
Debido a que los cables marinos sísmicos 30 son remolcados
detrás de la nave de estudios 20, las señales acústicas 42 (una señal acústica 42 a manera de ejemplo se representa gráficamente en la Figura 1), frecuentemente denominada "disparos", son producidas por las fuentes sísmicas 104 y se dirigen hacia abajo a través de una columna de agua 44 dentro de los estratos 62 y 68 debajo de la superficie del agua 24. Las señales acústicas 42 se reflejan desde las diversas formaciones geológicas subterráneas, tales como una formación 65 a manera de ejemplo que se representa gráficamente en la Figura 1.
Las señales acústicas incidentes 42 que son adquiridas por las fuentes 40 producen señales acústicas reflejadas correspondientes, u ondas de presión 60, las cuales son detectadas por los sensores sísmicos 58. Se observa que las ondas de presión que son recibidas y detectadas por los sensores sísmicos 58 incluyen ondas de presión "ascendente" que se propagan hacia los sensores 58 sin reflejos, así como también ondas de presión "descendente" que son producidas por reflejos de las ondas de presión 60 provenientes de una frontera aire-agua 31.
Los sensores sísmicos 58 generan señales (señales digitales, por ejemplo), llamadas "vestigios", las cuales indican las mediciones adquiridas del campo ondulatorio de presión y el movimiento de partículas (si los sensores son sensores de movimiento de partículas). Los vestigios son registrados y pueden procesarse al menos parcialmente por una unidad procesadora de señales 23 que se despliega en la nave de estudios 20, de acuerdo con algunas
modalidades de la invención. Por ejemplo, un sensor sísmico de multi-componente particular puede proporcionar un vestigio, el cual corresponde a una medición de un campo ondulatorio de presión por su hidrófono; y el sensor puede proporcionar uno o más vestigios que corresponden a uno o más componentes de movimiento de partículas, las cuales son medidas por sus acelerómetros.
El objetivo de la adquisición sísmica es conformar una imagen de un área de estudios para propósitos de identificación de formaciones geológicas subterráneas, tales como la formación geológica a manera de ejemplo 65. El análisis subsecuente de la representación puede revelar posiciones probables de yacimientos de hidrocarburos en formaciones geológicas subterráneas. Dependiendo de la modalidad particular de la invención, las porciones del análisis de la representación pueden realizarse en la nave de estudios sísmicos 20, tal como por la unidad procesadora de señales 23.
La formación de imágenes sísmica se utiliza para propósitos de producción de una imagen de la subsuperficie. Como se describe a continuación, los datos de medición de presión así como también los datos de gradiente de presión (de otro modo llamados datos vectoriales) pueden utilizarse para propósitos de generación de la imagen. Con objeto de preparar los datos sísmicos para la formación de imágenes, primeramente deben realizarse diversos pasos de procesamiento de datos. Por ejemplo, de acuerdo con algunas modalidades de la invención, los datos sísmicos adquiridos por los- sensores sísmicos pueden procesarse de acuerdo con una técnica 100 que se representa
gráficamente en la Figura 2.
Haciendo referencia a la Figura 2, la técnica 100 incluye procesar
(bloque 104) los datos sísmicos para eliminar, o atenuar, el ruido. A este respecto, la eliminación de ruido se refiere a la eliminación del ruido relacionado con la adquisición, tal como el ruido de vibración, ruido inducido por el flujo, ruido proveniente de naves sísmicas cercanas, ruido proveniente de plataformas petroleras cercanas, etc.
Después de que se atenúa el ruido, la técnica 100 incluye procesar (bloque 108) los datos sísmicos para eliminar las ondas directas y cualesquier múltiples ondas de reflejo. Más específicamente, el campo ondulatorio registrado puede descomponerse en ondas que se dispersan varias veces: la onda directa, la cual no está dispersa; ondas primarias, las cuales se dispersan una vez; ondas dispersas dos veces; etc. Cuando se genera una imagen de la subsuperficie, se supone que todas las ondas excepto las ondas primarias han sido eliminadas de los datos sísmicos. Se observa que el bloque 108 puede implementarse utilizando múltiples técnicas convencionales de eliminación.
Después de la eliminación múltiple que se realiza en el bloque
108, los datos sísmicos se procesan para determinar un modelo de velocidad de fondo, de conformidad con el bloque 112. Por ejemplo, el análisis de velocidad de migración puede realizarse para derivar el modelo de velocidad de fondo.
Después de la atenuación de ruido (bloque 104), la eliminación múltiple (bloque 108) y la derivación del modelo de velocidad de fondo (bloque 112), los datos sísmicos resultantes pueden procesarse después
(bloque 116) para determinar la imagen subsuperficial.
Como se describe en la presente, pueden utilizarse los datos de gradiente de presión por sí mismos o en conjunto con los datos de medición de presión para propósitos de generación de la imagen subsuperficial. Se supone que se han eliminado los múltiplos de los datos procesados, ya sea que estos datos incluyen o no los datos de medición de presión. Además, se supone en la presente que se ha derivado un modelo adecuado de velocidad de fondo utilizando, por ejemplo, análisis de velocidad de migración. La imagen subsuperficial puede derivarse de la siguiente manera. Las ondas de presión son modeladas por una aproximación de Born de primer orden, como se establece a continuación:
Ec.1 u((r,s,o)) = I ( 2g(r,x, a))c(x)cb3(x) g(x, s, ?)??,
donde "cb" representa el modelo de velocidad de fondo, el cual se supone que es conocido; "c" representa la perturbación de imagen, o función; "u ' representa la aproximación de Born de primer orden; - "g(r,s,w)" representa la función de Green, la cual corresponde al medio de fondo de las ondas excitadas en la fuente s y registradas en el receptor r.
La Ec. 1 supone que se ha aplicado la desconvolución de la fuente. Sin embargo, se observa que la desconvolución de la fuente no es esencial, ya que las técnicas que se describen en la presente pueden volver a aplicarse cuando no se ha aplicado la desconvolución de la fuente. Para el ejemplo específico representado gráficamente en la Ec.
1, una de las funciones de Green en la Ec. 1 se convoluciona con la forma de onda de la fuente.
El gradiente de presión puede derivarse al tomar la derivada de la presión u¡ con respecto a la posición del receptor r, comó se establece a continuación:
La Ec. 2 describe, al primer orden, la propagación de las ondas del gradiente disperso (es decir, el "vector") a través del medio.
Dadas las ondas de presión modelada y de gradiente de presión descritas con anterioridad, las cuales son funciones de la función de imagen c, la función de imagen c puede derivarse a través de un proceso de inversión que involucra minimizar una función de costos. En particular, la función de costos (llamada "F" en la presente) puede minimizarse con respecto a la función de imagen c, expuesta a continuación:
dF _
Ec.3 de
La función de costos F puede describirse de la siguiente manera:
donde "o", s °"ry" y "°"rz" son funciones de ponderación.
Dependiendo de la modalidad particular de la invención, una o más de las funciones de ponderación pueden establecerse igual a cero, aunque al menos una de las funciones de ponderación tiene un valor diferente a cero. Alternativamente, las funciones de ponderación pueden ser proporcionales a las matrices de covarianza, las cuales posiblemente pueden ser dependientes de la frecuencia. Alternativamente, pueden utilizarse otras técnicas para seleccionar los valores de ponderación.
En la Ecuación 4, las funciones w, wx, wy y wz son funciones de ponderación de formación de imágenes. Pueden establecerse iguales a uno, que es el caso de la mayoría de métodos de formación de imágenes que no se basan en métodos de rayos. Si las funciones de Green se calcula utilizando teoría de rayos, entonces las funciones de ponderación pueden determinarse utilizando el método desarrollado en Beylkin, G., 1985, Imaging of Discontínuities in The Inverse Scattering Problem by Inversión of a Causal Generalized Radon Transform (Formación de imágenes de discontinuidades en el problema de dispersión inversa por la inversión de una transformada de Radon generalizada causal), J. Math. Phys., 26, 99-108, en la cual se toman en cuenta los elementos fuera de la diagonal de la Hessiana. En la Ecuación 4, "u0(r,s,c«j)" representa la función de Green convolucionada con la forma de onda de la fuente.
A partir de la Ecuación 5, la imagen c puede derivarse utilizando, por ejemplo, las técnicas establecidas en Tarantola, a., 1984. Inversión of Seismic Reflection Data In The Acoustic Approximation (Inversión de
datos de reflejo sísmico en la aproximación acústica), Geophycis, 49, 1259-1266 y en Publicación PCT No. WO2008/081156, titulada "ACCURATE SEISMIC PROCESSING TECHNIQUES" ("TÉCNICAS DE PROCESAMIENTO SÍSMICO PRECISO"), la cual se publicó el 10 de julio de 2008. En particular, la diagonal de las ecuaciones normales puede tomarse, como lo hace Tarantola. Alternativamente, pueden utilizarse métodos asintóticos para tomar encuentra los elementos fuera de la diagonal eligiendo la función de ponderación de manera apropiada, como se establece en Beylkin.
Para propósitos de derivación de una expresión explícita de la función de imagen c, puede utilizarse el método de Tarantola y la función de ponderación puede seleccionarse igual a uno. En este caso, la función de imagen c(x) puede representarse de la siguiente manera:
donde "A" y "B" son de la siguiente manera:
B = o)2^(eír.(r, s, ?)— vr.(r, s, oi)gri(r, x, ?)? (?, s, (ú))drdsdu) Ec.7
Si únicamente se utiliza el gradiente vertical de los datos de presión, la función de imagen c(x) puede escribirse de la siguiente manera:
3 S ?2(?(t,5,? - Uoir.s^^g r.x.oú^gCx^.co^drdsdai Ec.8 c(x c0(x) J a)4,a2)g(r,x,(D)g(x,s,ü))drdsdcú\
Debe observarse que la función de imagen c(x) puede asumir muchas formas diferentes, dependiendo de la modalidad particular de la invención. Además de las técnicas expuestas con anterioridad, la formación de imágenes también puede aplicarse sobre puntos medios, o puntos medios y compensaciones (donde el punto medio se define como m = (s+r)/2 y la compensación se define como h= (s-r)/2), una técnica que se realiza típicamente en la migración de Kirchhoff. Alternativamente, la condición de formación de imágenes para cada fuente puede aplicarse separadamente, como se realiza en la migración de ecuación de ondas.
En otras modalidades de la invención, la función de Green puede calcularse utilizando una técnica de modelado de formas de onda, tal como la teoría de rayos, la formación de imágenes de Kirchhoff, técnicas de haces, las cuales originan las condiciones de formación de imágenes de haces, ecuaciones unidireccionales, las cuales genera la migración de ecuaciones de ondas y diferenciación finita, lo cual da como resultado una migración de tiempo inverso. Por ejemplo, para la siguiente función de Green:
g(r,s) = A(r,s)eicúT^ Ec.9
La formación de imágenes de Kirchhoff produce la siguiente función de imagen:
Éc.10
En el dominio del tiempo, la función de imagen puede expresarse
alternativamente de la siguiente manera:
/ [drz{r,s,t = T(r,x) + T(x,s)) - vrz(r,s,t = T(r,x) +T(x,s))^A(r,x)A(x,s)drds c(x) = c0 (x) 2 j ?2s?. \s(co A (r, x)A(x, s) | 2 drdsdco '
Ec.11
La migración de la ecuación de ondas puede representarse de la siguiente manera:
• j(drz(r,s^)-vrz(r,s^ ^gtz(r,x^)g}z(x,s^ drdsd(ti Ec-12 0 f (ú2a2\g?z(r,x,tú)gz(x,s,ü))\drdsdü)
donde "g4rz(r,x, )" representa la función de Green ascendente derivada del punto de dispersión hacia el receptor de la velocidad de partícula vertical y "g4rz(x,s, )" representa la función de Green correspondiente para formar el punto de dispersión a la fuente.
En resumen, la Figura 3 representa gráficamente una técnica 130 que puede aplicarse generalmente de acuerdo con las modalidades de la invención. La técnica 130 incluye la recepción (bloque 134) de datos sísmicos como los cuales son indicativos que mediciones de presión y gradiente de presión. Se modela una imagen de la subsuperficie (bloque 138) en función de las mediciones de presión y gradiente de presión. Se determina una imagen de la subsuperficie (bloque 142) con base en este modelo.
Más específicamente, de acuerdo con algunas modalidades de la invención, puede utilizarse una técnica 150 que se representa gráficamente en la Figura 4. De conformidad con la técnica 150, la primera aproximación de Born de una onda de presión se modela en función de la imagen, de conformidad con el bloque 154, y la onda de
gradiente de presión se modela como la derivada de la aproximación de Born de primer orden de la onda de presión con respecto a la posición del receptor en función de la imagen, de conformidad con el bloque 158. Después se determina una función de costos (bloque 162) con base en las mediciones de presión y gradiente de presión. Subsecuentemente se determina (bloque 166) un mínimo de la función de costos con respecto a la imagen para propósitos de determinación de la función de imagen c(x).
Haciendo referencia a la Figura 5, de acuerdo con algunas modalidades de la invención, un sistema de procesamiento de datos 320 puede implementar al menos parte de las técnicas que se describen en la presente, tales como técnicas relacionadas para recibir datos sísmicos indicativos de las mediciones de presión y gradiente de presión adquiridas en un estudio sísmico de al menos una formación subterránea; modelar una imagen de la(s) formación(es) subterránea(s) en función de las mediciones de presión y gradiente de presión; con base en el modelado, determinar una imagen de la(s) formación(es) subterránea(s); modelar una onda de presión de una aproximación de Born de primer orden; modelar una onda de gradiente de presión como una derivada de la aproximación de Born de primer orden de una onda de presión con respecto a la posición del receptor; ponderar selectivamente los componentes de una función de representación de imágenes asociada con diferentes componentes de movimiento de partículas; desplegar la imagen en una pantalla 374 del sistema 320; etc.
El sistema 320 puede ubicarse en uno de los cables marinos 30, en cada cable marino 30, distribuido entre los cables marinos 30, en la fuente sísmica 104, en la nave de estudios 30, en una instalación terrestre remota, etc. De acuerdo con algunas modalidades de la invención, el sistema 320 puede incluir un procesador 350, tal como uno o más microprocesadores y/o microcontroladores.
El procesador 350 puede acoplarse a una interfaz de comunicaciones 360 para propósitos de recepción de datos indicativos de mediciones sísmicas, parámetros de modelo, parámetros geofísicos, parámetros de estudios, etc. Los datos pertenecientes a las mediciones sísmicas pueden ser datos de presión, datos de multi-componente, etc.
Como ejemplo no limitante, la interfaz 360 puede ser una interfaz de bus serial USB, una interfaz de red, una interfaz de medio extraíbíe (tal como una tarjeta flash, CD-ROM, etc.) o una interfaz de almacenamiento magnético (interfaces IDE o SCSI, como ejemplos). Consecuentemente, la interfaz 360 puede asumir numerosas formas, dependiendo de la modalidad particular de la invención.
De acuerdo con algunas modalidades de la invención, la interfaz 360 puede acoplarse a una memoria 340 del sistema 320 y puede almacenar, por ejemplo, diversos conjuntos de datos 348 de entrada y/o salida involucrados con las técnicas que se describen en la presente. La memoria 340 puede almacenar instrucciones de programa 344, las cuales cuando son ejecutadas por el procesador 350, pueden ocasionar que el procesador 350 realice al menos parte de las técnicas que se describen en la presente y desplegar los resultados obtenidos mediante
la(s) técnica(s) en la batalla 374 del sistema 320, de acuerdo con algunas modalidades de la invención. Como se observa en la Figura 5, el sistema 320 puede incluir una interfaz de pantalla 370 que acopla el dispositivo de pantalla 374 al sistema 320.
Aunque se ha descrito la presente invención con respecto a un número limitado de modalidades, aquellos expertos en la materia, que tienen el beneficio de esta descripción, observarán numerosas modificaciones variaciones a partir de la misma. Se pretende que las reivindicaciones anexas cubran todas esas modificaciones y variaciones que se encuentran dentro del espíritu y alcance de la presente invención.
Claims (22)
1 . Un método que comprende: recibir datos sísmicos indicativos de las mediciones de presión y las mediciones de gradiente de presión adquiridas en un estudio sísmico de al menos una formación subterránea; modelar una imagen de al menos una formación subterránea en función de las mediciones de presión y las mediciones de gradiente de presión; y con base en el modelado, determinar la imagen.
2. El método según la reivindicación 1 , donde el modelado comprende: modelar una onda de presión en función de la imagen; modelar una onda de gradiente de presión en función de la imagen; determinar una función de costos con base en las mediciones de presión, las mediciones de gradiente de presión, la onda de presión modelada y la onda de gradiente de presión; y determinar un mínimo de la función de costos con respecto a la imagen.
3. El método según la reivindicación 2, donde el acto de modelar la onda de presión comprende: determinar una aproximación de Born de primer orden de la onda de presión.
4. El método según la reivindicación 2 , donde el acto de modelar la onda de gradiente de presión comprende: determinar una derivada de una aproximación de Born de primer orden de la onda de presión con respecto a una posición del receptor.
5. El método según la reivindicación 1 , donde la imagen se modela como una función que comprende los componentes asociados con diferentes ejes de gradiente, comprendiendo el método además: ponderar selectivamente los componentes.
6. El método según la reivindicación 5, donde el acto de ponderar selectivamente comprende: asignar una ponderación diferente a cero al menos a uno de los componentes.
7. El método según la reivindicación 5, donde el acto de ponderar selectivamente comprende: asignar una ponderación unitaria a cada uno de los componentes.
8. El método según la reivindicación 1 , donde los datos sísmicos comprenden datos adquiridos en un estudio sísmico remolcado.
9. El método según la reivindicación 1 , que comprende además: desplegar la imagen.
10. El método según la reivindicación 1 , donde las mediciones de gradiente de presión comprenden mediciones de velocidad de partícula.
1 1 . El método según la reivindicación 1 , que comprende además: procesar los datos sísmicos para eliminar el ruido proveniente de las mediciones de presión y de gradiente de presión.
12. El método según la reivindicación 1 , que comprende además: procesar los datos sísmicos para eliminar las múltiples ondas de reflejo provenientes de las mediciones de presión y de gradiente de presión.
1 3. El método según la reivindicación 1 , que comprende además: procesar los datos sísmicos para determinar un modelo de velocidad de fondo.
14. Un sistema, que comprende: una interfaz para recibir datos sísmicos indicativos de mediciones de presión y mediciones de gradiente de presión adquiridas en un estudio sísmico de al menos una formación subterránea; y un procesador para procesar los datos sísmicos a fin de determinar una imagen de al menos una formación subterránea con base en un modelo de la formación en función de las mediciones de presión y las mediciones de gradiente de presión.
1 5. El sistema según la reivindicación 14, donde la imagen se modela como una función que comprende componentes asociados con diferentes ejes de gradiente, y el procesador se adapta para ponderar selectivamente los componentes.
16. El sistema según la reivindicación 15, donde el procesador se adapta para asignar una ponderación diferente acero al menos a uno de los componentes.
17. El sistema según la reivindicación 15, donde el procesador se adapta para asignar una ponderación unitaria a cada uno de los componentes.
18. El sistema según la reivindicación 15, donde los datos sísmicos comprenden datos adquiridos en un estudio sísmico remolcado.
19. El sistema según la reivindicación 14, que comprende además: un cable marino para adquirir los datos sísmicos.
20. El sistema según la reivindicación 19, que comprende además: una nave para remolcar el cable marino.
21 . El sistema según la reivindicación 14, que comprende además: un dispositivo de pantalla, donde el procesador se adapta para desplegar la imagen en el dispositivo de pantalla.
22. El sistema según la reivindicación 14, donde el procesador se adapta además para procesar los datos sísmicos a fin de realizar al menos una de las siguientes: análisis de velocidad de migración, eliminación de múltiples reflejos y atenuación del ruido. RESUM EN U na técnica incluye recibir datos sísmicos, los cuales son indicativos de mediciones de presión y mediciones de gradiente de presión adquiridas en un estudio sísmico de al menos una formación subterránea. La técnica incluye modelar una imagen de la(s) formación(es) subterránea(s) en función de las mediciones de presión y de las mediciones de gradiente de presión. La técnica incluye determinar la imagen con base en el modelado.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US12/367,328 US8451687B2 (en) | 2009-02-06 | 2009-02-06 | Imaging with vector measurements |
| PCT/US2010/023095 WO2010091115A2 (en) | 2009-02-06 | 2010-02-03 | Imaging with vector measurements |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| MX2011008294A true MX2011008294A (es) | 2011-10-17 |
Family
ID=42540313
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| MX2011008294A MX2011008294A (es) | 2009-02-06 | 2010-02-03 | Formacion de imagenes con mediciones vectoriales. |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US8451687B2 (es) |
| EP (1) | EP2394187A2 (es) |
| AU (1) | AU2010210611A1 (es) |
| BR (1) | BRPI1008096A2 (es) |
| MX (1) | MX2011008294A (es) |
| WO (1) | WO2010091115A2 (es) |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9046626B2 (en) | 2011-01-10 | 2015-06-02 | Westerngeco L.L.C. | Performing reverse time imaging of multicomponent acoustic and seismic data |
| WO2012160430A2 (en) * | 2011-05-24 | 2012-11-29 | Geco Technology B.V. | Data acquisition |
| US9291734B2 (en) | 2011-06-10 | 2016-03-22 | International Business Machines Corporation | Full waveform inversion using combined shot data and no scratch disk |
| US9291735B2 (en) | 2011-06-10 | 2016-03-22 | Globalfoundries Inc. | Probablistic subsurface modeling for improved drill control and real-time correction |
| US9063248B2 (en) | 2011-06-10 | 2015-06-23 | International Business Machines Corporation | RTM seismic imaging using combined shot data |
| US8983779B2 (en) | 2011-06-10 | 2015-03-17 | International Business Machines Corporation | RTM seismic imaging using incremental resolution methods |
| US9176930B2 (en) * | 2011-11-29 | 2015-11-03 | Exxonmobil Upstream Research Company | Methods for approximating hessian times vector operation in full wavefield inversion |
| US8868347B2 (en) * | 2012-01-06 | 2014-10-21 | Baker Hughes Incorporated | Forward elastic scattering in borehole acoustics |
| CN102636811B (zh) * | 2012-04-10 | 2014-01-29 | 恒泰艾普石油天然气技术服务股份有限公司 | 一种海上二维地震资料中多次波的消除方法 |
| NO345402B1 (no) * | 2012-05-10 | 2021-01-18 | Cggveritas Services Sa | Fremgangsmåte, innretning og prosesseringsalgoritme for eliminasjon av multipler og støy fra marine seismiske data |
| US8971149B2 (en) * | 2012-05-10 | 2015-03-03 | Cggveritas Services Sa | Method, device and processing algorithm for multiple and noise elimination from marine seismic data |
| CN103105623B (zh) * | 2012-12-13 | 2013-08-21 | 东北石油大学 | 一种地震勘探中的数据波形处理方法 |
| US20140200820A1 (en) * | 2013-01-11 | 2014-07-17 | Westerngeco L.L.C. | Wavefield extrapolation and imaging using single- or multi-component seismic measurements |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4752916A (en) * | 1984-08-28 | 1988-06-21 | Dan Loewenthal | Method and system for removing the effect of the source wavelet from seismic data |
| US6021092A (en) * | 1997-06-09 | 2000-02-01 | Western Atlas International, Inc. | Method for deriving surface consistent reflectivity map from dual sensor seismic data |
| GB0106091D0 (en) * | 2001-03-13 | 2001-05-02 | Geco As | A method of determining the orientation of a seismic receiver, a seismic receiver, and a method of seismic surveying |
| 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 |
| 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 |
| 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 |
| GB2445162B (en) | 2006-12-28 | 2009-06-24 | Westerngeco Seismic Holdings | Geophysical signal processing |
| WO2008112036A1 (en) * | 2007-03-09 | 2008-09-18 | Ikelle Luc T | Imaging of multishot seismic data |
| US7426438B1 (en) * | 2007-03-16 | 2008-09-16 | Westerngeco L.L.C. | Technique to provide seismic data-based products |
| US9279899B2 (en) * | 2007-07-18 | 2016-03-08 | Westerngeco L.L.C. | System and technique to estimate physical propagation parameters associated with a seismic survey |
| 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 |
-
2009
- 2009-02-06 US US12/367,328 patent/US8451687B2/en not_active Expired - Fee Related
-
2010
- 2010-02-03 WO PCT/US2010/023095 patent/WO2010091115A2/en not_active Ceased
- 2010-02-03 EP EP10739081A patent/EP2394187A2/en not_active Withdrawn
- 2010-02-03 AU AU2010210611A patent/AU2010210611A1/en not_active Abandoned
- 2010-02-03 BR BRPI1008096A patent/BRPI1008096A2/pt not_active Application Discontinuation
- 2010-02-03 MX MX2011008294A patent/MX2011008294A/es not_active Application Discontinuation
Also Published As
| Publication number | Publication date |
|---|---|
| WO2010091115A2 (en) | 2010-08-12 |
| WO2010091115A3 (en) | 2010-11-18 |
| BRPI1008096A2 (pt) | 2016-03-08 |
| US20100202250A1 (en) | 2010-08-12 |
| US8451687B2 (en) | 2013-05-28 |
| AU2010210611A1 (en) | 2011-08-18 |
| EP2394187A2 (en) | 2011-12-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| MX2011008294A (es) | Formacion de imagenes con mediciones vectoriales. | |
| EP2283385B1 (en) | Jointly interpolating and deghosting seismic data | |
| AU2010213986B2 (en) | Using a rotation sensor measurement to attenuate noise acquired by a streamer-disposed sensor | |
| CN102016643B (zh) | 用于衰减双传感器地震拖缆中的低频噪声的方法 | |
| AU2010213830B2 (en) | Reconstructing seismic wavefields | |
| EP2191302B1 (en) | Removing vibration noise from multicomponent streamer measurements | |
| CN102062874B (zh) | 用于海洋地震拖缆数据的完全带宽消幻影的方法 | |
| AU2008254385B2 (en) | System and technique to remove perturbation noise from seismic sensor data | |
| MX2011006036A (es) | Uso de inversion de forma de onda para determinar las propiedades de un medio en el subsuelo. | |
| MX2012013721A (es) | Eliminar fantasmas de datos sísmicos. | |
| WO2005114258A1 (en) | Interpolation and extrapolation method for seismic recordings | |
| JPS58172571A (ja) | 海洋の地震探査方法および装置 | |
| WO2008008632A2 (en) | Time lapse marine seismic surveying employing interpolated multicomponent streamer pressure data | |
| EP2073043A2 (en) | Technique and system to cancel noise in measurements provided by sensors of a multi-component streamer | |
| US20100211322A1 (en) | Interpolating a pressure wavefield along an undersampled direction | |
| MX2011008541A (es) | Reconstruccion de un campo de onda sismica. | |
| EP2180345A1 (en) | Removing seismic interference using simultaneous or near simultaneous source separation | |
| AU2011205327A1 (en) | Deghosting and interpolating seismic data | |
| US20100172208A1 (en) | System and technique for local in-sea processing of particle motion data | |
| AU2015224508B2 (en) | Deghosting and interpolating seismic data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FA | Abandonment or withdrawal |