MX2011007084A - Adquisicion de datos y migracion antes del apilamiento con base en analisis de visibilidad sismica. - Google Patents
Adquisicion de datos y migracion antes del apilamiento con base en analisis de visibilidad sismica.Info
- Publication number
- MX2011007084A MX2011007084A MX2011007084A MX2011007084A MX2011007084A MX 2011007084 A MX2011007084 A MX 2011007084A MX 2011007084 A MX2011007084 A MX 2011007084A MX 2011007084 A MX2011007084 A MX 2011007084A MX 2011007084 A MX2011007084 A MX 2011007084A
- Authority
- MX
- Mexico
- Prior art keywords
- visibility
- seismic
- receiver
- source
- data
- Prior art date
Links
- 238000013508 migration Methods 0.000 title claims abstract description 42
- 230000005012 migration Effects 0.000 title claims abstract description 40
- 238000004458 analytical method Methods 0.000 title claims abstract description 16
- 238000000034 method Methods 0.000 claims abstract description 56
- 230000006870 function Effects 0.000 claims description 15
- 230000015654 memory Effects 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 abstract 1
- 238000012545 processing Methods 0.000 description 11
- 238000006073 displacement reaction Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000010304 firing Methods 0.000 description 3
- 238000005755 formation reaction Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000004880 explosion Methods 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/30—Analysis
-
- 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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/34—Displaying seismic recordings or visualisation of seismic data or attributes
- G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- 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/67—Wave propagation modeling
-
- 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/67—Wave propagation modeling
- G01V2210/673—Finite-element; Finite-difference
-
- 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/67—Wave propagation modeling
- G01V2210/675—Wave equation; Green's functions
-
- 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/67—Wave propagation modeling
- G01V2210/679—Reverse-time modeling or coalescence modelling, i.e. starting from receivers
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Un análisis de visibilidad sísmica de unas estructuras del subsuelo seleccionadas se emplea para determinar lugares de superficie que ofrecen una alta visibilidad de eventos objetivo. Estos lugares se pueden utilizar como base para la adquisición de datos de estudios sísmicos adicionales y/o seleccionar indicios existentes para volver a migrar con métodos de migración más sofisticados. Ya sea con el uso, se espera que los datos migrados nuevamente ofrezcan imágenes mejoradas del evento objetivo. En algunas modalidades, la determinación de la visibilidad incluye el uso de un propagador basado en una ecuación de onda para encontrar, para cada uno de los múltiples disparos simulados, un campo de ondas de reflexión del evento objetivo en un modelo sísmico; y calcular, para cada una de las múltiples posiciones del receptor, una señal de contribución de cada campo de ondas de reflexión. La determinación de la visibilidad además incluye convertir cada señal de contribución en un valor de visibilidad de la fuente-receptor. Debido a que adquirir y/o volver a migrar los datos se limita a la región seleccionada, el esfuerzo de representación gráfica para el evento objetivo se reduce significativamente.
Description
ADQUISICIÓN DE DATOS Y MIGRACIÓN ANTES DEL APILAMIENTO CON BASE EN ANÁLISIS DE VISIBILIDAD SÍSMICA
ANTECEDENTES DE LA INVENCIÓN
Los científicos e ingenieros a menudo emplean estudios sísmicos para la exploración, investigación geofísica y proyectos de ingeniería. Los estudios sísmicos pueden proporcionar información acerca de las estructuras subterráneas, incluyendo los límites de formación, tipos de roca, y la presencia o ausencia de depósitos de líquido. Esta información ayuda mucho a la búsqueda de agua, depósitos geotérmicos y depósitos de minerales tales como hidrocarburos y minerales. Particularmente, las compañías petroleras a menudo invierten en extensos estudios sísmicos para seleccionar los lugares para los pozos petroleros de exploración .
Los estudios sísmicos convencionales emplean fuentes de energía sísmica artificiales, como cargas de disparo, armas de aire comprimido o fuentes de vibración para generar ondas sísmicas. Cuando se disparan las fuentes se crea un "evento" sísmico, es decir, un pulso de energía sísmica que se propaga
como las ondas sísmicas desde la fuente hacia abajo dentro de la tierra. Las fallas y los límites entre las diferentes formaciones crean diferencias en la impedancia acústica que provocan reflexiones parciales de las ondas sísmicas. Un arreglo de sensores sísmicos detecta y registra estas reflexiones para su posterior análisis. Entonces, se aplican técnicas de procesamiento sofisticadas a las señales registradas para extraer una imagen u otra representación de la estructura del subsuelo.
Desafortunadamente, los analistas sísmicos a menudo encuentran que ciertas características del subsuelo son pobremente retratadas o inadecuadamente distinguibles. En tales circunstancias, las únicas soluciones son llevar a cabo una técnica de procesamiento más sofisticada o insistir en la adquisición de datos adicionales en el área previamente investigada. Cada una de estas soluciones puede ser muy cara en términos de tiempo y dinero.
SUMARIO DE LA INVENCIÓN
De conformidad, en la presente se divulgan sistemas y métodos para realizar el análisis de visibilidad sísmica de las estructuras del subsuelo seleccionadas. Estos sistemas y
métodos identifican las posiciones de la fuente sísmica y receptor que mejor puedan revelar los detalles de la estructura del subsuelo. Estas posiciones se pueden utilizar como base para la adquisición de nuevos datos sísmicos y/o someter a un subconjunto de datos existentes a un procesamiento de datos más sofisticado. Debido a que la región de adquisición de datos y/o procesamiento es muy reducida, los gastos asociados se reducen al mínimo.
Algunas modalidades ilustrativas del método incluyen un método de estudio sísmico que incluye: determinar la visibilidad de un evento objetivo en función de las posiciones de la fuente sísmica y del receptor; y adquirir datos sísmicos en una región seleccionada, al menos en parte, para incluir las posiciones que tienen valores de visibilidad por encima de un umbral. El evento objetivo puede ser ilustrado con base en los datos sísmicos recientemente adquiridos. Las modalidades del método ilustrado también incluyen un método de migración sísmica que incluye: determinar la visibilidad de un evento objetivo en las posiciones de la fuente y del receptor de los indicios en un estudio sísmico existente; y volver a migrar los indicios que tienen valores de visibilidad por encima de un umbral para ilustrar el evento objetivo. En ambos casos, la determinación
de la visibilidad puede incluir el uso de un propagador con base en ecuación de onda para encontrar, para cada uno de los múltiples disparos simulados, un campo de ondas de reflexión del evento objetivo en un modelo sísmico; y calcular, para cada una de las múltiples posiciones de receptor, una señal de contribución de cada campo de ondas de reflexión. La determinación de la visibilidad puede incluir además la convertir cada señal de contribución en un valor de visibilidad fuente-receptor.
BREVE DESCRIPCIÓN DE LOS DIBUJOS
Se podrá obtener una mejor comprensión de las diversas modalidades divulgadas cuando la descripción detallada sea considera en conjunto con los dibujo que la acompañan, en los que :
La figura 1 muestra un entorno de estudio sísmico ilustrativo;
La figura 2 muestra una forma geométrica de una fuente sísmica y receptor ilustrativos;
La figura 3 muestra indicios sísmicos ilustrativos;
La figura 4 muestra una imagen de fondo migrada para un modelo sísmico ilustrativo;
La figura 5a ilustra la energía sísmica de propagación de una fuente a un objetivo;
La figura 5b muestra un gráfico ilustrativo de la visibilidad del receptor;
Las figuras 6a a 6c muestran gráficos ilustrativos de la visibilidad de fuente bajo criterios diferentes;
La figura 7 muestra una función de visibilidad fuente-receptor ilustrativa;
La figura 8 muestra una imagen de fondo migrada ilustrativa para un subconjunto de datos sísmicos migrados;
La figura 9 muestra un diagrama de flujo de un método de análisis de visibilidad ilustrativo; y
La figura 10 muestra un sistema análisis de visibilidad ilustrativo.
Mientras que la invención es susceptible de varias modificaciones y formas alternativas, las modalidades específicas de la misma se muestran a modo de ejemplo en los dibujos y serán descritas en detalle en la presente. Sin embargo, se deberá entender que los dibujos y la descripción detallada no están destinados a limitar la invención a la forma particular divulgada, sino por el contrario, la intención es cubrir todas las modificaciones, equivalentes y
alternativas que entran dentro del alcance de las reivindicaciones que se acompañan.
DESCRIPCIÓN DETALLADA DE LA INVENCIÓN
Esta divulgación proporciona diversos métodos y sistemas de análisis de visibilidad que identifican los lugares de la fuente sísmica y/o receptor que mejor puedan medir los particulares de una o más características del subsuelo seleccionado. Los analistas puede centrar sus esfuerzos de adquisición y procesamiento en estas regiones para mejorar el detalle de imagen para las características seleccionadas. Los sistemas y los métodos descritos se entienden mejor cuando se describe en un contexto de uso ilustrativo.
En consecuencia, la figura 1 muestra un entorno de estudio sísmico ilustrativo en el que se colocan un arreglo de receptores sísmicos (102) en un arreglo espaciado-aparte en la superficie de la tierra (104) para detectar las ondas sísmicas. Los receptores (102) se acoplan de forma inalámbrica o mediante cable a una unidad de adquisición de datos (106) que recibe, procesa y almacena los datos de la
señal sísmica recogidos por los receptores. Una fuente de energía sísmica (108) (por ejemplo, un camión vibrador) se dispara en múltiples posiciones para generar ondas de energía sísmica que se propagan a través de la tierra (110) y que se reflejan de las discontinuidades de impedancia acústica para alcanzar a los receptores (102) . Estas discontinuidades pueden ser creadas por las fallas, los límites entre las camas de formación y los límites entre los fluidos de la formación. Las discontinuidades aparecerán como puntos brillantes en la representación de la estructura del subsuelo que se deriva a partir de los datos de señal sísmica.
La figura 1 muestra además un modelo del subsuelo ilustrativo que se utiliza como ejemplo en esta divulgación. En este modelo, la tierra tiene cuatro capas de formación relativamente planas con un límite fuertemente curvado entre las capas tercera y cuarta. La velocidad del sonido en cada una de las capas desde arriba hasta abajo es de 2000 m/s, 3500 m/s, 2800 m/s y 4000 m/s, respectivamente. Se deberá tener en cuenta que para propósitos de ilustración, la figura 1 no está dibujada a escala. La extensión horizontal del arreglo de receptor se espera que esté limitada a un par de kilómetros o menos, mientras que la dimensión horizontal del
modelo del subsuelo es de dieciséis kilómetros y la dimensión vertical es de cinco kilómetros.
La figura 2 muestra una geometría ilustrativa para las posiciones de la fuente (202) y las posiciones del receptor (204), como se verían desde una vista superior. Geometrías viables de estudio sísmico son infinitas en variedad y pueden recogerse en una base sector-por-sector, una base de geometría móvil, una base de arreglo en movimiento y varias combinaciones de las mismas. El mensaje principal aquí es que el número de señales del receptor adquirido en respuesta a cada disparo de la fuente sísmica sea bastante grande y cuando el número de disparos de la furnte se tiene en cuenta, el número resultante de los indicios puede alcanzar fácilmente a millones.
La figura 3 muestra algunos indicios de señal sísmica
(S1-S3) recibidos ilustrativos. Los indicios pueden representar desplazamiento, velocidad, aceleración, presión o alguna otra medida de la intensidad de la energía sísmica en función del tiempo. La señal recibida por cada receptor suele ser muestreada y digitalizada a entre 8 y 32 bits de resolución a una velocidad de cerca de 500 muestras por segundo para una duración de unos 30 segundos después de cada disparo. En algunos casos, los receptores captan los datos de
múltiples componentes, aumentando aún más la cantidad de datos muestra para cada indicio. Los datos de indicio pueden ser filtrados y comprimidos antes de su almacenamiento. Los datos del estudio sísmico almacenados se transportan o comunican de otra manera a un centro de procesamiento de datos .
Una red de computadoras en el sitio de procesamiento de datos procesa los datos para estimar la distribución volumétrica de las velocidades del sonido mediante técnicas conocidas. Véase, por ejemplo, Jon F. Claerbout, Fundamentos de Procesamiento de Datos Geofísicos, p. 246-56, que se incorpora a la presente como referencia. Por otra parte, la distribución de velocidad puede estar disponible independiente de otras fuentes, por ejemplo, registros de pozo. Con la distribución de velocidad estimada, el sitio de procesamiento de datos "migra" los indicios sísmicos, traduciéndolos de funciones de tiempo a funciones de profundidad .
Existen varias técnicas de migración, incluyendo las que están basadas en la migración de la ecuación de onda de un solo sentido (one-way EM, por sus siglas en inglés) y la migración inversa en tiempo basada en la ecuación de onda de sentido completo (RTM, por sus siglas en inglés) . La
migración one-way WEM es una técnica muy popular y ampliamente utilizada debido a que es eficaz en muchos casos y es relativamente barata en términos de complejidad computacional . Sin embargo, en áreas con estructuras complejas, especialmente aquellas que generan fuertes ondas anuladas (por ejemplo, las ondas de prisma) y reflexiones múltiples (que pueden crear ondas dúplex) , la migración one-way WEM simplemente falla al representar gráficamente las estructuras complejas. Esta situación puede agravarse en los estudios sísmicos marinos ya que el uso de receptores con cobertura azimutal estrecha limita la cantidad de compensación de la línea de cruce.
La migración RTM es capaz de hacer frente a esos problemas de representación gráfica. (Véase, por ejemplo, E. Baysal, D.D. Kosloff y J.W.C. Sherwood, "Migración inversa en tiempo", Geofísica, 48, 1514
[1983]; G.A. McMechan, "Migración por extrapolación de los valores límite dependientes del tiempo", Prospección Geofísica, 31, 413 -420
[1983]; y N.D. Whitmore, "Representación gráfica iterativa de profundidad por la propagación hacia atrás del tiempo", SEG Resúmenes Expandidos, 2, 382-385
[1983]) . En años recientes, la RTM se ha vuelto más atractiva para el procesamiento de imágenes de profundidad previo apilamiento en medios
complejos y estructuras subsal . Sin embargo, en comparación con el método one-way EM, la RTM es computacionalmente costosa y requiere de un sitio de procesamiento de datos para tener computadoras con memorias y discos de gran capacidad. Por otra parte, la RTM se hace aún más difícil cuando la migración de componentes de alta frecuencia del campo de onda debido a la dispersión numérica del esquema de diferencias finitas .
La figura 4 muestra una imagen de fondo migrada (402) para el modelo sísmico que se ilustra en la figura 1. Un centenar de disparos se simularon mediante modelaje de frente de diferencias finitas. La posición de disparo fue de 6 km hasta 14 km con un intervalo de 80 metros de tiro. La apertura del receptor para cada tiro cubrió todo el modelo. El intervalo del receptor es de 10 metros. La imagen de profundidad de la figura 4 se obtuvo mediante la migración de todos los tiros con la abertura total. En general, el evento curveado está bien representado gráficamente, pero la amplitud del evento de pendiente empinada (404) es débil. Para obtener una mejor imagen de este evento objetivo, no es necesario obtener grandes cantidades de datos de estudio nuevo o volver a procesar todos los datos existentes con la RTM. En su lugar, se puede centrar solamente en los indicios
sísmicos que tienen contribuciones significativas a la proyección gráfica del evento objetivo. Los métodos y sistemas de análisis de visibilidad sísmica descritos a continuación proporcionarán una identificación cuantitativa de los indicios que proporcionan esas contribuciones.
La figura 5 corresponde a una ilustración esquemática de trazado de rayos de la energía sísmica de propagación sísmica. El rayo (502) muestra la propagación de energía de una fuente a un extremo del evento objetivo y de regreso a la superficie, mientras que el rayo (504) muestra la propagación de energía de la fuente hasta el otro extremo del evento objetivo y de regreso a la superficie. En la práctica, la migración de ecuación de onda se emplea para propagar la energía sísmica hacia abajo así como de regreso hacia la superficie. Más específicamente, el análisis de visibilidad se pone en ejecución preferentemente utilizando un propagador basado en la ecuación de onda en lugar de un enfoque de alta frecuencia asintótica basado en rayos. El método de la ecuación de onda es más preciso y proporciona un resultado más confiable.
El análisis de visibilidad se lleva a cabo en dos fases. En primer lugar el campo de onda de un tiro simulado se propaga hacia abajo y el software mide el campo de ondas de
reflexión en el evento objetivo. En la segunda fase, el campo de ondas de reflexión se propaga hacia atrás y el software mide la contribución del objetivo en las señales registradas por cada receptor. El receptor-fuente de visibilidad V(s,r) del evento objetivo se mide al integrar el cuadrado de la señal de contribución medida Csr(t) (similar a elevar al cuadrado y sumar cada uno de los valores muestra para un rastro en la figura 3) :
donde r es la posición del receptor, s es la posición de la fuente y Csr(t) es la señal de contribución medida en función del tiempo entre el tiempo de disparo t = 0 y el final del intervalo de registro t = T. Las posiciones del receptor y del disparo simulado pueden ser uniformemente espaciadas a lo largo del área del modelo o pueden ser personalizadas con el entorno del estudio contemplado (por ejemplo, una geometría de buque sísmico marino) .
La visibilidad del receptor VR(r) se define como la visibilidad de la fuente-receptor V(s,r) para una posición determinada de la fuente s = S:
VR (r) = V{S, r) .
La figura 5b muestra la visibilidad del receptor (506) del evento objetivo (404) para la posición de disparo de la fuente que se muestra en la figura 5a. Se puede observar que para esta posición de la fuente, la visibilidad del receptor se limita principalmente a una zona cercana bien definida sobre el evento objetivo.
El término de visibilidad de la fuente Vs(s) se define en esta descripción como la suma de la visibilidad de la fuente-receptor V(s,r) en todas las posiciones del receptor
{R} :
re{R}
La figura 6 muestra la visibilidad de la fuente de evento objetivo (404) bajo tres supuestos diferentes. La figura 6a muestra la visibilidad de la fuente con apertura completa del receptor (es decir, todos los receptores pueden responder a todos los disparos de la fuente) . La figura 6b muestra la visibilidad de la fuente asumiendo una abertura de ±2 km (es decir, sólo los receptores dentro de los 2 kilómetros de la posición de la fuente puede responder a los disparos de la fuente) . La figura 6c muestra visibilidad de la fuente con desplazamiento a cero (es decir, el receptor
está co-ubicado con la fuente)1. En cada caso, la visibilidad
varia con la posición de la fuente, y los disparos en la zona
cercana de arriba y a la derecha del evento objetivo que más
contribuyen a la visibilidad del evento objetivo. Una
comparación de estas funciones de visibilidad de la fuente
permite cuantificar fácilmente los efectos de la apertura del
receptor. Una gran cantidad de visibilidad se pierde si sólo
se considera el caso con desplazamiento a cero.
La figura 7 muestra un mapa de la función de visibilidad
de la fuente-receptor V(s,r) para el evento objetivo (404) .
El eje horizontal indica la ubicación del receptor, mientras
que el eje vertical indica la ubicación de la fuente. Los
indicios de las posiciones de la fuente-receptor que se
encuentran dentro de la región (702) contribuyen más con la
visibilidad del evento (404) . La función de visibilidad a lo
largo de la linea (704) es la función de visibilidad del
receptor VR(r) para la posición de la fuente que se muestra
en la figura 5a. La función de visibilidad de la fuente con
desplazamiento a cero que se muestra en la figura 6c se
encuentra a lo largo de la linea a 45° (708) . La visibilidad
1 Por otra parte, se observa que los cálculos de la visibilidad con desplazamiento cero se pueden simplificar utilizando el supuesto del "reflector de explosión", en el que los eventos objetivo son tratados como una serie de puntos de explosión distribuidos cada uno con una fuerza igual a la reflectividad del evento objetivo.
de la fuente que se muestra en la figura 6c se encuentra mediante la integración horizontal a través de toda la figura. La función de visibilidad de la fuente que se muestra en la figura 6b se encuentra mediante la integración horizontal entre las lineas (706) .
La apertura del receptor máxima (que corresponde a la longitud del cable en los estudios sísmicos) se puede seleccionar mediante el ajuste de la separación entre las líneas (706) para capturar la mayor parte de la zona no cero bajo la función de visibilidad. Las posiciones de la fuente pueden seleccionarse para capturar la mayor parte de la zona no cero bajo la función de visibilidad de la fuente. Con esta estrategia para seleccionar los indicios (y, en su caso, la adquisición de datos) para la migración de profundidad previo apilamiento reduce en gran medida la cantidad de esfuerzo necesario para mejorar la representación gráfica del evento objetivo. La figura 8 muestra una imagen de la migración de profundidad utilizando sólo los datos seleccionados, lo que demuestra que el comportamiento de la amplitud de la imagen controlada de visibilidad es más equilibrada.
La figura 9 muestra un método de análisis de visibilidad ilustrativo que puede ser llevado a cabo por un sistema informático de forma automática o bajo la dirección de un
usuario, tal como un analista de depósito. Empezando en el bloque (902), el sistema obtiene los datos que representan una imagen de migración de profundidad. En muchos casos, estos datos han sido generados por el propio sistema basado en un estudio sísmico anterior. En el bloque (904) el sistema identifica uno o más eventos objetivo, por ejemplo, características que inadecuadamente se hayan representado de manera gráfica. En algunas implementaciones , el sistema identifica los eventos objetivo mediante la visualización de la representación gráfica de migración de profundidad a un usuario y solicitar la entrada del usuario sobre las áreas que parecen estar definidas adecuada o inadecuadamente.
En el bloque (906), el sistema selecciona un método de migración que es más sofisticado que el utilizado para generar la migración de datos original. Por ejemplo, la migración original podría haber empleado la migración one-way WEM, pero el sistema puede ser capaz de aplicar la migración RTM de onda completa. El usuario puede seleccionar el método de migración deseado en donde múltiples mejoras estén disponibles.
En el bloque (908), el sistema determina la visibilidad de la fuente-receptor V(s,r), utilizando el método de migración seleccionada para simular disparos en la estructura
tentativa identificada en los datos de profundidad migrados originales. Como se mencionó anteriormente, la visibilidad está determinada por el cálculo de los campos de onda de reflexión de los eventos objetivo para cada una de las múltiples posiciones de la fuente y luego medir las contribuciones de señal de estos campos de onda de reflexión para las señales medidas en cada una de las múltiples posiciones del receptor.
En el bloque (910), el sistema identifica los indicios existentes cuyas posiciones de la fuente-receptor tienen visibilidades de los eventos objetivo por encima de un umbral determinado. El umbral puede ser programado, basado en un valor de visibilidad máxima o seleccionado para capturar una fracción predeterminada (por ejemplo, el 90%) del área bajo la superficie de la visibilidad multidimensional . En el bloque (912), el sistema aplica el método de migración seleccionada para identificar los indicios de alta visibilidad. Debido a que se espera que los indicios identificados representen un pequeño subconjunto de los datos disponibles, el uso de un método de migración más sofisticado pueden ser eminentemente factible.
En el bloque (914), el sistema determina si los eventos objetivo se han representado gráficamente de manera adecuada
y, de ser así, el método se desplazará hasta el bloque (922). En algunas implementaciones , el sistema hace esta determinación al desplegar los datos de profundidad migrados a un usuario y solicitar retroalimentación del usuario. Si el evento objetivo sigue siendo representado gráficamente de manera insuficiente, se espera necesaria la adquisición de datos adicionales. En consecuencia, en el bloque (916), el sistema identifica una región del estudio y otros parámetros del estudio basado por lo menos en parte en los cálculos de visibilidad de la fuente-receptor. En algunas implementaciones, el rango conveniente de posiciones de la fuente y del receptor puede ser determinado al dibujar un rectángulo (para los estudios de tierra) o un paralelogramo (para estudios marinos en donde la posición de receptor varía con la posición de la fuente) que encierre sustancialmente la mayor parte de la región de valores de alta visibilidad.
En el bloque (918), el sistema obtiene los datos de los indicios del nuevo estudio y en el bloque (920) el método de migración seleccionada se aplica a generar una nueva imagen de los datos de profundidad migrados de la región que contiene los eventos objetivo. En el bloque (922), una imagen combinada se sintetiza y se despliega. La imagen combinada incluye la estructura completa identificada a partir de los
datos migrados originales, pero también incluye las imágenes de los eventos objetivo en los datos que acaban de migrar. Un ingeniero de yacimientos puede evaluar el potencial de producción con las estructuras de interés adecuadamente definidas para análisis.
La figura 10 muestra un sistema de análisis de visibilidad ilustrativo en la forma de una computadora (50) que tiene uno o más de dispositivos de entrada (54) y uno o más dispositivos de salida (56) a través de los cuales puede interactuar con un usuario. El software (ilustrado como discos portátiles de almacenamiento de información (52)) configura el (los) procesador (es ) de la(s) computadora ( s ) para recibir comandos del usuario y recuperar datos de manera responsable de red o de almacenamiento interno, procesa los datos y genera imágenes para el usuario para ver y analizar. Al aplicar los métodos descritos, el software normalmente puede distribuir la carga de procesamiento de los datos entre varios equipos interconectados por una red.
La descripción anterior se basa en un modelo sísmico 2D con fines explicativos. En la práctica, es de esperar que un volumen 3D sea representado gráficamente y que cada una de las posiciones de la fuente y el receptor sean especificadas en términos de al menos dos coordenadas espaciales. Como una
consecuencia, se espera que el mapa de visibilidad de la fuente-receptor (véase la figura 7) tenga al menos cuatro dimensiones espaciales. Sin embargo, los principios subyacentes son los mismos.
En resumen, unos sistemas y métodos de análisis de visibilidad sísmica han sido descritos. Estos sistemas y métodos cuantitativamente identifican las posiciones de la fuente y del receptor deseables en la superficie de un evento objetivo en medios complejos. La fuerza de la visibilidad para una geometría fuente-receptor determinada indica si un evento objetivo es visible o invisible con esa geometría. Tal conocimiento es aplicado al diseño de estudios de adquisición y a la migración de profundidad antes del apilamiento. Los experimentos de visibilidad proporcionan los siguientes puntos :
Para un evento objetivo determinado, los -datos de estudio fuera del área de alta visibilidad no son necesarios.
• Para un evento objetivo determinado, volver a migrar los indicios fuera del área de alta visibilidad no son necesarios .
• Si un evento objetivo determinado no tiene un área de alta visibilidad, no puede ser reconstruido con el método de migración seleccionada.
Numerosas variaciones y modificaciones serán evidentes para los expertos en la materia una vez que la descripción anterior sea completamente apreciada. Se pretende que las siguientes reivindicaciones deban interpretarse para comprender todas las variaciones y modificaciones.
Claims (21)
- REIVINDICACIONES Un método de estudio sísmico que comprende: determinar la visibilidad de un evento objetivo en función de las posiciones de la fuente y del receptor sísmicos; y adquirir datos sísmicos en una región seleccionada, al menos en parte para incluir las posiciones que tienen valores de visibilidad por encima de un umbral. El método de la reivindicación 1, que comprende: representar gráficamente el evento objetivo basado en los datos sísmicos adquiridos. El método de la reivindicación 2, en donde determinar la visibilidad de un evento objetivo incluye: encontrar, para cada uno de múltiples disparos simulados, un campo de ondas de reflexión del evento objetivo en un modelo sísmico; calcular, para cada una de las múltiples posiciones del receptor, una señal de contribución de cada campo de ondas de reflexión; y convertir cada señal de contribución en un valor de visibilidad de la fuente-receptor. El método de la reivindicación 3, en donde el modelo sísmico está determinado con base en al menos en parte sobre los datos migrados de un estudio sísmico. El método de la reivindicación 3, en donde el evento objetivo es una característica del subsuelo que es seleccionado por un usuario como inadecuadamente representado de manera gráfica. El método de la reivindicación 3, en donde el valor de visibilidad es un cuadrado integrado de la señal de contribución . El método de la reivindicación 1, en donde el umbral se selecciona para capturar de una fracción predeterminada del área bajo la función de visibilidad multidimensional . Un método de migración sísmica que comprende: determinar la visibilidad de un evento objetivo en las posiciones de la fuente y el receptor de indicios en un estudio sísmico existente; y volver a migrar los indicios que tienen valores de visibilidad por encima de un umbral para representar gráficamente el evento objetivo. El método de la reivindicación 8, en donde determinar la visibilidad de un evento objetivo incluye: encontrar, para cada uno de múltiples disparos simulados, un campo de ondas de reflexión del evento objetivo en un modelo sísmico; calcular, para cada una de las múltiples posiciones del receptor, una señal de contribución de cada campo de ondas de reflexión; y convertir cada señal de contribución en un valor de visibilidad de la fuente-receptor. El método de la reivindicación 9, en donde el modelo sísmico está determinado con base en al menos en parte sobre los datos migrados de un estudio sísmico existente, en donde los datos migrados se obtienen mediante un primer método de migración. El método de la reivindicación 10, en donde convertir incluye determinar un cuadrado integrado de la señal de contribución. El método de la reivindicación 10, en donde la visibilidad se determina utilizando un segundo método de migración diferente al primer método de migración. El método de la reivindicación 10, en donde volver a migrar se realiza mediante un segundo método de migración diferente al primer método de migración. El método de la reivindicación 8, en donde el evento objetivo es una característica del subsuelo que es seleccionado por un usuario como inadecuadamente representado de manera gráfica. 15. El método de la reivindicación 8, en donde el umbral está seleccionado para capturar una fracción predeterminada del área bajo la función de visibilidad. 16. Un sistema de análisis de visibilidad que comprende : una memoria que almacena un software de análisis de visibilidad; y por lo menos un procesador acoplado a la memoria para ejecutar el software de análisis de visibilidad, en donde el software configura el al menos un procesador para : determinar, para cada uno de múltiples disparos, un campo de ondas de reflexión de al menos un evento objetivo en un modelo sísmico; calcular, para cada una de las múltiples posiciones del receptor, una señal de contribución de cada campo de ondas de reflexión; y convertir cada señal de contribución a un valor de visibilidad de la fuente-receptor. 17. El sistema de la reivindicación 16, en donde además el software configura al procesador para visualmente mapear la visibilidad en función de al menos una posición de la fuente y posición del receptor. 18. El sistema de la reivindicación 16, en donde el modelo sísmico está determinado con base en al menos en parte sobre los datos migrados de un estudio sísmico, en donde los datos migrados se obtienen mediante un primer método de migración. 19. El sistema de la reivindicación 18, en donde la visibilidad se determina utilizando un segundo método de migración diferente al primer método de migración. 20. El sistema de la reivindicación 16, en donde además el software configura al procesador para identificar indicios de un estudio sísmico existente para volver a migrar con un método de migración mejorado, en donde los indicios son identificados con base en al menos en parte sus valores de visibilidad. 21. El sistema de la reivindicación 16, en donde además el software configura al procesador para identificar posiciones de la fuente y del receptor para la adquisición de datos sísmicos adicionales, en donde las posiciones son identificadas con base en al menos en parte los valores de visibilidad de la fuente-receptor.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/US2009/031400 WO2010082938A1 (en) | 2009-01-19 | 2009-01-19 | Data acquisition and prestack migration based on seismic visibility analysis |
Publications (1)
Publication Number | Publication Date |
---|---|
MX2011007084A true MX2011007084A (es) | 2011-11-04 |
Family
ID=42340028
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
MX2011007084A MX2011007084A (es) | 2009-01-19 | 2009-01-19 | Adquisicion de datos y migracion antes del apilamiento con base en analisis de visibilidad sismica. |
Country Status (8)
Country | Link |
---|---|
US (1) | US9329288B2 (es) |
EP (1) | EP2376947B1 (es) |
CN (1) | CN102282481B (es) |
AU (1) | AU2009337134B2 (es) |
BR (1) | BRPI0924056A2 (es) |
CA (1) | CA2747146C (es) |
MX (1) | MX2011007084A (es) |
WO (1) | WO2010082938A1 (es) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2443436B (en) * | 2006-11-04 | 2010-07-21 | Westerngeco Seismic Holdings | Method of determing properties of the earth using turning wave data |
US8116168B1 (en) | 2008-06-18 | 2012-02-14 | Halliburton Energy Services, Inc. | Hybrid one-way and full-way wave equation migration |
CN102282481B (zh) | 2009-01-19 | 2014-11-12 | 兰德马克图形公司 | 基于地震能见度分析的数据采集和叠前偏移 |
WO2010120301A1 (en) | 2009-04-16 | 2010-10-21 | Landmark Graphics Corporation | Seismic imaging systems and methods employing a fast target-oriented illumination calculation |
US8830788B2 (en) | 2011-02-24 | 2014-09-09 | Landmark Graphics Corporation | Sensitivity kernal-based migration velocity analysis in 3D anisotropic media |
US9164184B2 (en) | 2011-05-27 | 2015-10-20 | Conocophillips Company | Reciprocal method two-way wave equation targeted data selection for seismic acquisition of complex geologic structures |
WO2012166733A1 (en) | 2011-05-27 | 2012-12-06 | Conocophillips Company | Two-way wave equation targeted data selection for improved imaging of prospects among complex geologic structures |
EP2715404A4 (en) * | 2011-05-27 | 2015-01-28 | Conocophillips Co | ABOUT A BIDIRECTIONAL WAVE EQUALIZATION DATA SELECTION FOR THE SEISMIC COLLECTION OF COMPLEX GEOLOGICAL STRUCTURES |
WO2013048585A1 (en) * | 2011-09-28 | 2013-04-04 | Conocophillips Company | Reciprocal method two-way wave equation targeted data selection for improved imaging of complex geologic structures |
US10061046B2 (en) * | 2012-02-06 | 2018-08-28 | Ion Geophysical Corporation | Integrated passive and active seismic surveying using multiple arrays |
GB2510873A (en) | 2013-02-15 | 2014-08-20 | Total Sa | Method of modelling a subsurface volume |
GB2510872A (en) * | 2013-02-15 | 2014-08-20 | Total Sa | Method of modelling a subsurface volume |
US20150226868A1 (en) * | 2014-02-10 | 2015-08-13 | Microseismic, Inc. | Method for microseismic event moment magnitude estimation |
CN104570055A (zh) * | 2015-02-04 | 2015-04-29 | 中国海洋石油总公司 | 一种***震源控制*** |
EP3371749A1 (en) | 2015-11-06 | 2018-09-12 | Google LLC | Regularizing machine learning models |
FR3059700A3 (fr) * | 2016-12-02 | 2018-06-08 | Al Khatib Habib Agissant Au Nom Et Pour Le Compte De La Soc Spotlight En Cours De Formation | Procede d’amelioration des systemes actifs ultralegers de detection sismique |
US11137509B2 (en) * | 2018-02-02 | 2021-10-05 | Fairfield Industries, Inc. | Seismic imaging with a temporal decomposition imaging condition |
CN112782754A (zh) * | 2019-11-06 | 2021-05-11 | 中国石油天然气集团有限公司 | 多源同步采集数据模拟方法及装置 |
US11474267B2 (en) | 2020-06-11 | 2022-10-18 | China Petroleum & Chemical Corporation | Computer-implemented method and system employing compress-sensing model for migrating seismic-over-land cross-spreads |
Family Cites Families (42)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5014230A (en) | 1986-06-05 | 1991-05-07 | Xox Corporation | Solid-modeling system using topology directed subdivision for determination of surface intersections |
US4890242A (en) | 1986-06-05 | 1989-12-26 | Xox Corporation | Solid-modeling system using topology directed subdivision for determination of surface intersections |
US4953142A (en) * | 1989-01-06 | 1990-08-28 | Marathon Oil Company | Model-based depth processing of seismic data |
US5138584A (en) | 1989-09-06 | 1992-08-11 | Chevron Research & Technology Company | Migration of seismic turning waves |
US5274605A (en) | 1992-06-26 | 1993-12-28 | Chevron Research And Technology Company | Depth migration method using Gaussian beams |
US5490120A (en) | 1993-04-30 | 1996-02-06 | Union Oil Company Of California | Method of imaging overturned waves |
US5530679A (en) | 1993-05-10 | 1996-06-25 | Western Atlas International, Inc. | Method for migrating seismic data |
US5544126A (en) | 1993-11-04 | 1996-08-06 | Berryhill; John R. | System and method of seismic shot-record migration |
US6002642A (en) | 1994-10-19 | 1999-12-14 | Exxon Production Research Company | Seismic migration using offset checkshot data |
US5784334A (en) * | 1996-03-13 | 1998-07-21 | Atlantic Richfield Company | Method and system for detecting hydrocarbon reservoirs using amplitude versus offset analysis of seismic signals |
GB9609612D0 (en) * | 1996-05-08 | 1996-07-10 | Geco Prakla Uk Ltd | A method of and apparatus for determining the quality of seismic data |
US6446007B1 (en) | 1998-11-05 | 2002-09-03 | Exxonmobil Upstream Research Company | Method for controlled-amplitude prestack time migration of seismic data |
US6021094A (en) | 1998-12-03 | 2000-02-01 | Sandia Corporation | Method of migrating seismic records |
US6687618B2 (en) | 2000-08-07 | 2004-02-03 | 3D Geo Development, Inc. | Typing picks to horizons in migration velocity analysis |
WO2002023222A1 (en) | 2000-09-15 | 2002-03-21 | Nutec Sciences, Inc. | Illumination weighted imaging condition for migrated seismic data |
US20020033832A1 (en) * | 2000-09-18 | 2002-03-21 | Rafail Glatman | Method for computer modeling of visual images and wave propagation |
US6611761B2 (en) | 2000-12-19 | 2003-08-26 | Schlumberger Technology Corporation | Sonic well logging for radial profiling |
US6864890B2 (en) | 2002-08-27 | 2005-03-08 | Comoco Phillips Company | Method of building and updating an anisotropic velocity model for depth imaging of seismic data |
US6778909B1 (en) | 2002-10-05 | 2004-08-17 | 3Dgeo Development, Inc. | Seismic data processing systems and methods |
US6904368B2 (en) * | 2002-11-12 | 2005-06-07 | Landmark Graphics Corporation | Seismic analysis using post-imaging seismic anisotropy corrections |
US6819628B2 (en) | 2003-04-07 | 2004-11-16 | Paradigm Geophysical (Luxembourg) S.A.R.L. | Wave migration by a krylov space expansion of the square root exponent operator, for use in seismic imaging |
FR2858063B1 (fr) * | 2003-07-21 | 2005-09-16 | Geophysique Cie Gle | Procede d'estimation du taux de couverture d'illumination dans le domaine migre |
US7065004B2 (en) | 2004-04-22 | 2006-06-20 | Pgs Americas, Inc. | Horizon-based residual depth migration velocity analysis |
US7602937B2 (en) | 2004-06-08 | 2009-10-13 | International Electronic Machines Corporation | Image-based visibility measurement |
US7376517B2 (en) | 2005-05-13 | 2008-05-20 | Chevron U.S.A. Inc. | Method for estimation of interval seismic quality factor |
US20060270537A1 (en) * | 2005-05-27 | 2006-11-30 | Karol Fisher | Exercising system with trellis |
US20070162249A1 (en) | 2006-01-06 | 2007-07-12 | Min Lou | Traveltime calculation in three dimensional transversely isotropic (3D TTI) media by the fast marching method |
US7196969B1 (en) | 2006-02-09 | 2007-03-27 | Pgs Geophysical As | Three-dimensional two-way acoustic wave equation pre-stack imaging systems and methods |
US7952960B2 (en) | 2006-10-03 | 2011-05-31 | Bp Corporation North America Inc. | Seismic imaging with natural Green's functions derived from VSP data |
US8120991B2 (en) | 2006-11-03 | 2012-02-21 | Paradigm Geophysical (Luxembourg) S.A.R.L. | System and method for full azimuth angle domain imaging in reduced dimensional coordinate systems |
US7400553B1 (en) * | 2006-11-30 | 2008-07-15 | Shengwen Jin | One-return wave equation migration |
US20080137480A1 (en) | 2006-12-07 | 2008-06-12 | Woodside Energy Limited | Method of Building a Subsurface Velocity Model |
US8194498B2 (en) | 2008-01-18 | 2012-06-05 | Westerngeco L.L.C. | Using a wave propagator for transversely isotropic media |
US7663972B2 (en) | 2008-02-22 | 2010-02-16 | Pgs Geophysical As | Method for three dimensional seismic travel time tomography in transversely isotropic media |
US8116168B1 (en) | 2008-06-18 | 2012-02-14 | Halliburton Energy Services, Inc. | Hybrid one-way and full-way wave equation migration |
US8296069B2 (en) | 2008-10-06 | 2012-10-23 | Bp Corporation North America Inc. | Pseudo-analytical method for the solution of wave equations |
CN102282481B (zh) | 2009-01-19 | 2014-11-12 | 兰德马克图形公司 | 基于地震能见度分析的数据采集和叠前偏移 |
WO2010120301A1 (en) | 2009-04-16 | 2010-10-21 | Landmark Graphics Corporation | Seismic imaging systems and methods employing a fast target-oriented illumination calculation |
US8332156B2 (en) | 2009-07-10 | 2012-12-11 | Chevron U.S.A. Inc. | Method for propagating pseudo acoustic quasi-P waves in anisotropic media |
US8406081B2 (en) | 2009-09-25 | 2013-03-26 | Landmark Graphics Corporation | Seismic imaging systems and methods employing tomographic migration-velocity analysis using common angle image gathers |
AU2009354764B2 (en) | 2009-11-02 | 2014-03-20 | Landmark Graphics Corporation | Seismic imaging systems and methods employing a 3D reverse time migration with tilted transverse isotropy |
US8830788B2 (en) | 2011-02-24 | 2014-09-09 | Landmark Graphics Corporation | Sensitivity kernal-based migration velocity analysis in 3D anisotropic media |
-
2009
- 2009-01-19 CN CN200980154752.0A patent/CN102282481B/zh not_active Expired - Fee Related
- 2009-01-19 EP EP09838520.6A patent/EP2376947B1/en active Active
- 2009-01-19 AU AU2009337134A patent/AU2009337134B2/en not_active Ceased
- 2009-01-19 MX MX2011007084A patent/MX2011007084A/es active IP Right Grant
- 2009-01-19 US US13/142,478 patent/US9329288B2/en active Active
- 2009-01-19 BR BRPI0924056-0A patent/BRPI0924056A2/pt not_active IP Right Cessation
- 2009-01-19 CA CA2747146A patent/CA2747146C/en active Active
- 2009-01-19 WO PCT/US2009/031400 patent/WO2010082938A1/en active Application Filing
Also Published As
Publication number | Publication date |
---|---|
EP2376947A1 (en) | 2011-10-19 |
CN102282481B (zh) | 2014-11-12 |
EP2376947A4 (en) | 2015-09-02 |
US20110273959A1 (en) | 2011-11-10 |
BRPI0924056A2 (pt) | 2015-07-14 |
AU2009337134A1 (en) | 2011-07-07 |
CA2747146A1 (en) | 2010-07-22 |
US9329288B2 (en) | 2016-05-03 |
EP2376947B1 (en) | 2021-08-25 |
WO2010082938A1 (en) | 2010-07-22 |
CA2747146C (en) | 2016-10-04 |
AU2009337134B2 (en) | 2013-11-21 |
CN102282481A (zh) | 2011-12-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2009337134B2 (en) | Data acquisition and prestack migration based on seismic visibility analysis | |
Groos et al. | Application of a complete workflow for 2D elastic full-waveform inversion to recorded shallow-seismic Rayleigh waves | |
US10989825B2 (en) | Method and system for determining source signatures after source ghost removal | |
AU2008282547B2 (en) | Method for determining seismic data quality | |
EP3710867B1 (en) | Noise attenuation of multiple source seismic data | |
US11609349B2 (en) | Determining properties of a subterranean formation using an acoustic wave equation with a reflectivity parameterization | |
Neducza | Stacking of surface waves | |
WO2012141805A2 (en) | Stable shot illumination compensation | |
Xiang et al. | Improving the resolution of impedance inversion in karst systems by incorporating diffraction information: A case study of Tarim Basin, China | |
Jardin et al. | Depth seismic imaging using reflection and first arrival traveltime tomography: Application to a deep profile across the Northern Emirates Foothills | |
US20230305176A1 (en) | Determining properties of a subterranean formation using an acoustic wave equation with a reflectivity parameterization | |
Song | Locating Petroleum Sources Using Dsp Techniques | |
O'Neill et al. | Shallow surface wave inversion? research and applications in Western Australia |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FG | Grant or registration |