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
Application number
MX2011007084A
Other languages
English (en)
Inventor
Shengwen Jin
Shiyong Xu
Original Assignee
Landmark Graphics Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Landmark Graphics Corp filed Critical Landmark Graphics Corp
Publication of MX2011007084A publication Critical patent/MX2011007084A/es

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/003Seismic data acquisition in general, e.g. survey design
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/34Displaying seismic recordings or visualisation of seismic data or attributes
    • G01V1/345Visualisation of seismic data or attributes, e.g. in 3D cubes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V20/00Geomodelling in general
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/67Wave propagation modeling
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/67Wave propagation modeling
    • G01V2210/673Finite-element; Finite-difference
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/67Wave propagation modeling
    • G01V2210/675Wave equation; Green's functions
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/67Wave propagation modeling
    • G01V2210/679Reverse-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)

  1. 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.
MX2011007084A 2009-01-19 2009-01-19 Adquisicion de datos y migracion antes del apilamiento con base en analisis de visibilidad sismica. MX2011007084A (es)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

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