ES2223156T3 - Metodo de determinacion de cambios espaciales en la estratigrafia de la estructura del subsuelo, litologia y contenido de fluidos y de reduccion del ruido sismico. - Google Patents
Metodo de determinacion de cambios espaciales en la estratigrafia de la estructura del subsuelo, litologia y contenido de fluidos y de reduccion del ruido sismico.Info
- Publication number
- ES2223156T3 ES2223156T3 ES99201782T ES99201782T ES2223156T3 ES 2223156 T3 ES2223156 T3 ES 2223156T3 ES 99201782 T ES99201782 T ES 99201782T ES 99201782 T ES99201782 T ES 99201782T ES 2223156 T3 ES2223156 T3 ES 2223156T3
- Authority
- ES
- Spain
- Prior art keywords
- seismic
- data
- horizons
- grid
- amplitudes
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
- 238000000034 method Methods 0.000 title claims abstract description 99
- 239000012530 fluid Substances 0.000 title claims description 40
- 230000009467 reduction Effects 0.000 title description 5
- 230000008859 change Effects 0.000 claims abstract description 43
- 238000004458 analytical method Methods 0.000 claims abstract description 21
- 239000011435 rock Substances 0.000 claims description 64
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000001914 filtration Methods 0.000 claims description 13
- 238000009826 distribution Methods 0.000 claims description 12
- 238000005457 optimization Methods 0.000 claims description 4
- 238000004590 computer program Methods 0.000 claims description 2
- 230000009189 diving Effects 0.000 claims description 2
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims 1
- 230000000875 corresponding effect Effects 0.000 description 8
- 229930195733 hydrocarbon Natural products 0.000 description 7
- 150000002430 hydrocarbons Chemical class 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 239000004215 Carbon black (E152) Substances 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 238000006073 displacement reaction Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000008021 deposition Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000003628 erosive effect Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008520 organization Effects 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000013479 data entry Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000002513 implantation Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000003860 storage Methods 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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Método para determinar y analizar cambios espaciales en el subsuelo de la tierra y para reducir ruido en datos sísmicos, que comprende las etapas de: a) obtener datos sísmicos; b) obtener un modelo terrestre del subsuelo con capas representativas de la estructura del subsuelo; c) definir una serie de microcapas definidas por una serie de horizontes de microcapas entre los límites de las capas del modelo terrestre; d) determinar puntos de la cuadrícula de los horizontes de microcapas a lo largo de un horizonte de microcapas; e) determinar las amplitudes de los datos sísmicos en los puntos de la cuadrícula de los horizontes de microcapas; f) determinar en cada punto de la cuadrícula de los horizontes de microcapas una o más determinaciones representativas de la velocidad de cambio de las amplitudes de los datos sísmicos y/o determinaciones de la dirección en que cambian las amplitudes de los datos sísmicos y/o determinar en cada punto de la cuadrícula de los horizontes de microcapas una o mássalidas del filtro; g) almacenar como salida dicha amplitud de los datos sísmicos, dichas determinaciones y dichas salidas del filtro en cada punto de la cuadrícula de los horizontes de microcapas para análisis adicional del subsuelo;
Description
Método de determinación de cambios espaciales en
la estratigrafía de la estructura del subsuelo, litología y
contenido de fluidos y de reducción del ruido sísmico.
Esta invención se refiere a la materia general de
reducir ruido en datos de propiedades de las rocas, sísmicos y
deducidos de sísmicos; deduciendo diversos datos nuevos a partir de
datos de propiedades de las rocas, sísmicos y deducidos de
sísmicos, de entrada, que destacan cambios espaciales en
estructura, estratigrafía, litología y fluidos de las rocas, del
subsuelo; y al análisis e interpretación de tales datos.
Los datos sísmicos se adquieren para proporcionar
información acerca de la estructura, estratigrafía, litología y
fluidos contenidos en las rocas, del subsuelo. Los registros de
datos sísmicos adquiridos son la respuesta de una fuente de energía
sísmica después de pasar a través, y que se están reflejando, por
las rocas en el subsuelo. Se pueden adquirir datos sísmicos en, o
cerca de, la superficie terrestre o se pueden adquirir a lo largo
de pozos de sondeo. Después de la adquisición, los datos sísmicos
se procesan típicamente para una serie de trazas sísmicas, en que
cada traza representa la respuesta sísmica en una cierta posición
x, y de la superficie. La traza misma consta de una serie de
muestras de la respuesta sísmica, normalmente ordenadas para
corresponder con tiempo de desplazamiento sísmico creciente o,
después de conversión a profundidad, profundidad creciente.
Dependiendo de la geometría de adquisición, las trazas sísmicas
normalmente se procesan y se organizan para formar líneas con
trazas espaciadas regularmente a lo largo de la superficie. Los
datos sísmicos a lo largo de tales líneas, se pueden considerar
como secciones a través de la tierra. Los datos sísmicos se
refieren como datos sísmicos en 2D cuando las líneas están en
direcciones diferentes o están separadas distantes en relación con
el espaciamiento de las trazas. Los datos sísmicos se refieren como
datos sísmicos en 3D cuando la adquisición es de manera que el
procesamiento da como resultado una serie de líneas sísmicas que
están organizadas secuencialmente y en que las posiciones de las
trazas x, y forman una cuadrícula regular y de manera que el
espaciamiento de las líneas sísmicas, en general, está dentro del
mismo orden de magnitud que el espaciamiento de las trazas dentro
de las líneas. En la práctica, las líneas por las que se adquieren
los datos se denominan líneas alineadas y las líneas ortogonales a
las líneas alineadas se refieren como líneas cruzadas. La Fig. 1
muestra una sección sísmica con una serie de trazas de datos
sísmicos tomados a partir del cubo de los datos sísmicos en 3D, de
la que la cuadrícula x,y se muestra en la figura 2. Las series de
datos sísmicos en 2D y 3D se analizan y se interpretan con
posterioridad, en general, en estaciones de trabajo de ordenador
con programas informáticos especializados, para revelar la
estructura, estratigrafía, litología y fluidos, del subsuelo y para
predecir así la posición, estructura, estratigrafía, litología y
distribución de fluidos de yacimientos de hidrocarburos, acuíferos
asociados y otras características del subsuelo de interés. La Fig.
3 muestra una interpretación estructural de los datos sísmicos de
la Fig. 1. Esta interpretación delimita la zona del yacimiento
total, dentro de la que se correlacionan amplitudes sísmicas altas
con arenas petrolíferas. La interpretación también muestra
relaciones estructurales y estratigráficas. Las relaciones
estructurales típicamente tienen que ver con fallamiento, por
ejemplo, en la Fig. 3 la interpretación muestra cómo se dislocan
las capas definidas por los horizontes, por las fallas. Las
relaciones estratigráficas típicamente tienen que ver con
deposición y erosión. Por ejemplo, una interpretación puede mostrar
cómo una superficie de erosión trunca capas que se encuentran
más
abajo.
abajo.
Las amplitudes de los datos sísmicos se
determinan principalmente por la intensidad de la reflexión de las
ondas sísmicas en los límites de la capa. La intensidad de la
reflexión, a su vez, se determina por cambios en ciertos parámetros
físicos de las rocas cuando van de una capa a la siguiente. Estos
parámetros físicos se determinan por las propiedades físicas de la
matriz de la roca, es decir, la roca con los poros de la roca
vacíos y los fluidos contenidos en los poros, referidos
conjuntamente como "datos de propiedades de las rocas". Los
cambios en la matriz de la roca pueden estar causados por cambios en
la litología (composición mineral de la roca y dislocación). Los
cambios en los fluidos surgen de cambios en el tipo de fluido:
agua, petróleo y gas; o cambios en las propiedades de los tipos de
fluidos. Usando algoritmos de ordenador modernos, se pueden estimar
datos de propiedades de las rocas a partir de las amplitudes de los
datos sísmicos. Los datos de propiedades de las rocas que se pueden
estimar directamente a partir de datos sísmicos incluyen: impedancia
acústica, impedancia elástica, velocidad de ondas de presión,
velocidad de ondas de cizallamiento y densidad. También se pueden
deducir datos de propiedades de las rocas, adicionales, directa o
indirectamente, usando relaciones funcionales, estadísticas u otras,
entre las diferentes propiedades de las rocas. Se pueden usar
directamente datos de propiedades de las rocas deducidos de
sísmicos, para analizar cambios en la litología y los fluidos en
las capas. También, se mantiene información acerca de la estructura
y la estratigrafía y, con frecuencia, incluso se exalta en relación
con los datos sísmicos. El uso de datos de propiedades de las rocas
deducidos de sísmicos, en análisis e interpretación del subsuelo
con frecuencia se prefiere, por lo tanto, sobre el uso de datos de
reflexión sísmica o se hace junto con análisis e interpretación del
subsuelo de datos sísmicos. Por la misma razón, el método de la
materia se aplica preferiblemente a datos de propiedades de las
rocas deducidos de sísmicos. La Fig. 4 muestra la misma sección que
la Fig. 1, pero no muestra una sección por un cubo de la propiedad
de las rocas de impedancia acústica procedente de los datos de
reflexión sísmica. Los cambios en impedancia acústica dan como
resultado cambios en los coeficientes de reflexión sísmica. En
otras palabras, la impedancia acústica es una propiedad de la capa
mientras los coeficientes de reflexión sísmica tienen que ver con
las interfases de las capas. El análisis de las diferencias entre
los datos sísmicos y los datos de impedancia acústica revela que
las arenas que soportan petróleo y sus límites se pueden
interpretar con más precisión a partir de los datos de impedancia
acústica que a partir de los datos sísmicos mismos.
Para caracterizar un horizonte o plano de falla
interpretado, se pueden calcular el buzamiento y el azimut en cada
punto del horizonte. Como se ilustra en la Fig. 5, el buzamiento en
un punto del horizonte, es el ángulo desde la vertical al vector
gradiente de un plano tangente a la superficie del horizonte, en el
punto del horizonte. El azimut es el ángulo de la proyección del
vector gradiente sobre un plano horizontal calculado en el sentido
de las agujas del reloj, en general, en relación con el Norte.
Un aspecto clave de los datos de propiedades de
las rocas, sísmicos y deducidos de sísmicos, es que, en general,
estos datos no contienen suficiente información para determinar, en
cada muestra, todos los datos requeridos acerca de la estructura,
estratigrafía, litología y fluido en esa muestra. Se proporciona
información adicional para analizar e interpretar variaciones
espaciales en los datos de propiedades de las rocas, sísmicos y
deducidos de sísmicos. Por ejemplo, se puede determinar a partir de
la naturaleza de un cambio espacial si el cambio es debido a un
cambio en la estructura, por ejemplo, una falla, o debido a un
cambio en la litología o en los fluidos. El problema es que la
información acerca de las variaciones espaciales con frecuencia no
se aprecia fácilmente o se cuantifica en seguida a partir de los
datos de propiedades de las rocas, sísmicos o deducidos de
sísmicos. Esto motiva la necesidad de métodos que generen datos que
destaquen cambios espaciales en la estructura, estratigrafía,
litología y fluidos, del subsuelo y de métodos para analizar e
interpretar tales datos.
Se han descrito métodos que se centran en
calcular ciertas determinaciones de la discontinuidad espacial
usando sólo datos sísmicos. Estos métodos no utilizan la información
captada en una interpretación de los datos sísmicos. El método
propuesto se aparta de los métodos existentes por el empleo de un
modelo del subsuelo, basado en una interpretación disponible, para
conducir el cálculo de nuevos tipos de determinaciones de cambios
espaciales en estructura, estratigrafía, litología y fluidos, del
subsuelo. Estas determinaciones proceden de cambios en las
amplitudes de datos sísmicos o datos de propiedades de las rocas
deducidos de sísmicos, a lo largo de horizontes. Una de tales
determinaciones es el gradiente de las amplitudes, por distinción
con el gradiente de una superficie geométrica, referido como el
gradiente de la propiedad. Este gradiente de la propiedad se
determina en cada punto del horizonte, ajustando una superficie a
amplitudes en el punto del horizonte y puntos del horizonte
circundantes y calculando el gradiente de esta superficie. Los
gradientes grandes corresponden a cambios laterales rápidos. Un
método alternativo para caracterizar los cambios de amplitud es
suavizando los datos de amplitud a lo largo del horizonte por
filtración, y después tomando la diferencia de los datos filtrados
y de entrada como determinación de la velocidad de cambio de las
amplitudes. Tales operaciones de filtración también reducen ruido y
como tales proporcionan una nueva manera de reducir ruido en los
datos de propiedades de las rocas sísmicos y deducidos de sísmicos,
de entrada.
La presente invención proporciona un
procedimiento nuevo y mejorado para reducir ruido en datos de
propiedades de las rocas, sísmicos y deducidos de sísmicos; para
generar datos que revelen información acerca de la variación
espacial de la estructura, estratigrafía, litología y contenido en
fluido, del subsuelo, a partir de datos de propiedades de las rocas
sísmicos y deducidos de sísmicos; y para analizar e interpretar
esos datos. El método usa como entrada datos de propiedades de las
rocas sísmicos o deducidos de sísmicos y una interpretación de
estos datos. A partir de la interpretación se construye un modelo
del subsuelo, en lo siguiente referido como modelo terrestre. Los
cálculos del procedimiento se conducen por este modelo terrestre. En
resumen, las etapas principales para la versión en 3D del
procedimiento son:
Obtener datos de propiedades de las rocas
sísmicos o deducidos de sísmicos; interpretar estos datos para
definir los horizontes y las fallas que determinan límites de capas
de interés; a partir de los horizontes y las fallas y las
relaciones estratigráficas y estructurales entre los horizontes y
las fallas construir un modelo terrestre en que los horizontes y
las fallas de entrada formen los límites de las capas principales
del modelo terrestre; guiado por las relaciones estratigráficas y
estructurales, subdividir cada capa principal del modelo terrestre
en microcapas unidas por horizontes de microcapas definidos en
puntos de la cuadrícula x, y de propiedades de las rocas sísmicas o
deducidas de sísmicas, de manera que los horizontes de microcapas
definan la estructura interna de cada una de las capas principales
del modelo terrestre; para cada punto de la cuadrícula de cada
horizonte de microcapas encontrar las coordenadas espaciales de una
serie de puntos de la cuadrícula, circundantes, sobre el horizonte
de microcapas, que definan juntos un segmento de la superficie del
horizonte de microcapas; rotar el segmento de la superficie del
horizonte de microcapas en la dirección de la línea alineada y la
línea transversal por un intervalo de ángulos definidos por el
usuario, con un tamaño de etapa definido por el usuario alrededor
del punto de definición actual para definir las coordenadas
espaciales del segmento de la superficie del horizonte de
microcapas, rotado; para este segmento de la superficie del
horizonte de microcapas, rotado, extraer las correspondientes
amplitudes de los datos de propiedades de las rocas, sísmicos o
deducidos de sísmicos; a partir de estas amplitudes extraídas
calcular una o múltiples determinaciones para la velocidad de
cambio de estas amplitudes a lo largo del segmento de horizonte de
microcapas rotado, actual, calculando el tamaño del gradiente o
usando diversos filtros, la dirección del gradiente y las salidas
del filtro; repetir para todos los ángulos tanto en la dirección de
la línea alineada como en la de la línea transversal; para cada
determinación de velocidad de cambio determinar los ángulos en que
es mínima la determinación de velocidad de cambio y para esos
ángulos almacenar como datos de salida el buzamiento y el azimut del
ángulo de rotación, la determinación de velocidad de cambio, la
dirección del gradiente y las salidas del filtro; repetir este
procedimiento para todos los puntos de la cuadrícula sobre cada
horizonte de microcapas y para todos los horizontes de la
microcapa; imprimir los resultados en forma de una serie de
horizontes de microcapas en que cada punto del horizonte de
microcapas contiene los correspondientes datos de salida del
procedimiento; generar salida adicional por interpolación de los
horizontes de microcapas para las cuadrículas de propiedades de las
rocas sísmicas y deducidas de sísmicas; analizar e interpretar los
datos de salida que contienen la información sobre cambios
espaciales para predecir variaciones laterales en la estructura,
estratigrafía, litología y distribución de fluidos, del subsuelo. El
análisis e interpretación de los datos de salida interpolados, se
puede hacer sobre estaciones de trabajo sísmicas clásicas, usando
herramientas de visualización e interpretación de sección, mapa y
3D. Los propios horizontes de microcapas se pueden analizar e
interpretar de una nueva manera, según la cual los horizontes de
microcapas se visualizan en vista de mapa (véanse las Fig. 8, 9 y
10) o en 3D, y en que el usuario puede someter a un ciclo de
operaciones por la pila de horizontes de microcapas para volver a
visualizar los cambios a lo largo de los horizontes de
microcapas.
Los datos de salida generados destacan
información acerca de variaciones laterales en estructura,
estratigrafía, litología y distribución de fluido, de subsuelo, no
directamente evidentes en los datos de propiedades de las rocas,
sísmicos o deducidos de sísmicos, de entrada. Adicionalmente, los
datos de salida generados contienen versiones filtradas de los datos
de propiedades de las rocas, sísmicos o deducidos de sísmicos, de
entrada, en que se reduce ruido y, por lo tanto, se pueden usar
para mejorar el análisis e interpretación de datos de propiedades
de las rocas, sísmicos o deducidos de sísmicos, clásicos.
El procedimiento se conduce por un modelo
geológico, que en general no tenderá a captar cambios detallados en
la estructura, estratigrafía, litología y distribución de fluido,
del subsuelo. El procedimiento compensa las faltas de precisión
geométricas en el modelo geológico por el procedimiento de
perturbación del ángulo especificado. De hecho, se pueden generar
datos útiles con el procedimiento de la materia, en base a un
modelo muy simple. En su forma más simple tal modelo constará de
una capa unida por dos horizontes paralelos. El procedimiento
también se puede llevar a cabo sin perturbación del ángulo
(intervalo de ángulos = 0). En tal caso, se puede usar la salida
para evaluar los cambios espaciales en la estructura,
estratigrafía, litología y distribución de fluidos en relación con
el modelo geométrico mismo. En la aplicación práctica, el
procedimiento se puede aplicar veces consecutivas en que la salida
se usa para mejorar el modelo geométrico, que a su vez se usa para
generar datos mejorados acerca de cambios espaciales en la
estructura, estratigrafía, litología y distribución de fluidos, del
subsuelo.
La invención es aplicable en particular a la
prospección, desarrollo y producción de hidrocarburos, para
determinar la estructura, estratigrafía, litología y distribución
de fluidos en yacimientos de hidrocarburos y acuíferos asociados y
para determinar el movimiento de fluido a partir de estudios
sísmicos repetidos en el tiempo sobre un yacimiento a medida que se
agota. Los datos generados por el procedimiento revelan cómo cambia
espacialmente la estructura, estratigrafía, litología y
distribución de fluidos, y cómo son de rápidos tales cambios. Los
cambios rápidos pueden indicar límites de yacimientos totales o
límites entre diferentes unidades de yacimiento y contactos de
fluido. Los cambios más suaves pueden indicar, por ejemplo,
porosidad creciente o decreciente y el % de roca que soporta
hidrocarburo en relación con roca que no soporta hidrocarburo. Los
datos generados por el procedimiento pueden revelar detalles acerca
de cambios espaciales en la estructura, estratigrafía, litología y
distribución de fluidos que no se detectan fácilmente cuando se
trabaja con los datos de propiedades de las rocas, sísmicos y
deducidos de sísmicos. De manera similar, cuando se aplica el
procedimiento para repetir estudios sísmicos, se puede conseguir
detección exaltada del movimiento en el tiempo de los límites de
fluido.
El procedimiento no se limita a la aplicación en
prospección, producción y desarrollo de hidrocarburos. Cualquier
análisis e interpretación de datos de propiedades de las rocas,
sísmicos o deducidos de sísmicos, con el propósito de determinar la
estructura, estratigrafía, litología y distribución de fluidos, del
subsuelo, se puede beneficiar del procedimiento.
La Fig. 1 es un ejemplo de una sección sísmica a
partir de un cubo sísmico en 3D. Se muestran trazas con posiciones x
que oscilan de 32 a 110 y una posición y 251. La respuesta sísmica
en cada traza se muestra como función del tiempo t de
desplazamiento. Las amplitudes sísmicas se rellenan con negro para
una desviación positiva y no se rellenan para una desviación
negativa. El tamaño de las desviaciones indica la intensidad de las
reflexiones subyacentes desde la tierra.
La Fig. 2 muestra la cuadrícula x,y de las trazas
que comprenden el cubo en 3D a partir del cual se origina la sección
en la Fig. 1. La línea en la posición y 251 y las posiciones x que
oscilan de 32 a 110 se corresponden con la sección sísmica de la
Fig. 1.
La Fig. 3 muestra la sección sísmica de la Fig. 1
con superposición de una interpretación de algunas de las
principales características estructurales y estratigráficas del
subsuelo. Tales interpretaciones se obtienen normalmente usando
estaciones de trabajo de gráficos de ordenador con programas
informáticos de interpretación sísmica especiales. Los horizontes y
las fallas se digitalizan en la pantalla del ordenador guiados por
los datos sísmicos visualizados. La interpretación sobre la Fig. 3
muestra el buzamiento de las unidades estratigráficas dentro de las
que se encuentran las arenas del yacimiento y muestra algunas de
las fallas que están dislocando el yacimiento en diferentes
unidades.
La Fig. 4 muestra una sección de datos de
propiedades de las rocas, deducidos de sísmicos, en este caso la
impedancia acústica, a lo largo de la misma línea como en la
Fig.1.
La Fig. 5 ilustra la definición del buzamiento y
el azimut para caracterizar la geometría de una superficie.
La Fig. 6 es un organigrama que muestra las
etapas del procedimiento en una realización del nuevo método.
La Fig. 7 es un ejemplo del modelo terrestre
construido a partir de la interpretación del de la Fig. 3, que
muestra cómo se conducen los horizontes de microcapas dentro de las
capas definidas por la interpretación de entrada.
La Fig. 8 muestra el tamaño absoluto del
gradiente de la propiedad sobre los puntos de la cuadrícula de un
horizonte de microcapas. Las amplitudes superiores indicadas por las
flechas muestran claramente una falla.
La Fig. 9 muestra el azimut del gradiente de la
propiedad sobre los puntos de la cuadrícula de un horizonte de
microcapas diferente. Las amplitudes superiores indicadas por las
flechas muestran claramente fallas múltiples.
La Fig. 10 muestra el tamaño absoluto del
gradiente de la propiedad sobre los puntos de la cuadrícula de aún
otro horizonte de microcapas que está pasando por algunas arenas del
yacimiento. La magnitud del gradiente de la propiedad indica
claramente diferencias en la continuidad lateral de las arenas del
yacimiento. En el área A la magnitud del gradiente de la propiedad
delimita una arena del yacimiento con buena continuidad lateral,
mientras en el área B la rápida variación lateral de la magnitud
del gradiente de la propiedad indica que las arenas del yacimiento
presentan continuidad lateral deficiente.
Esta invención se puede expresar en muchas formas
diferentes. La exposición y descripción de la invención en los
dibujos y en esta descripción son ilustrativas y aclaratorias de la
misma, y se pueden hacer diversos cambios en la secuencia de las
etapas del procedimiento, de los parámetros en el procedimiento y
de los detalles del procedimiento, sin apartarse del alcance de la
invención.
La Fig. 6 ilustra las etapas en el nuevo
procedimiento en la forma de un organigrama. Hay tres tipos de datos
de entrada para el procedimiento: datos sísmicos, datos de
propiedades de las rocas deducidos de sísmicos y un denominado
"modelo terrestre". La entrada obligatoria consta de un modelo
terrestre y al menos una serie de datos de propiedades de las rocas,
sísmicos o deducidos de sísmicos.
La etapa 100 muestra la entrada de los datos
sísmicos. Si bien esto no es una condición necesaria, se asume que
los datos de entrada se han procesado a la forma de una serie de
datos en 2D o 3D para interpretación de datos sísmicos, es decir,
las series de datos constan de series de trazas sísmicas organizadas
en líneas en que cada traza tiene una coordenada x,y y cada punto
de datos de la traza corresponde a un cierto tiempo de
desplazamiento sísmico o profundidad (t o z).
Tradicionalmente se obtiene una serie de datos
del procesamiento para análisis e interpretación adicionales. Se
aplican, cada vez más, procedimientos especiales para exaltar
ciertas características del subsuelo en relación con una serie de
datos procesados de manera rutinaria. Por ejemplo, se pueden generar
series de datos sísmicos que contengan información acerca de ondas
sísmicas que se hayan reflejado por un cierto intervalo de ángulos,
ya que éstos pueden exaltar la posibilidad de detectar límites de
fluido. Por lo tanto, las múltiples series de datos sísmicos que
cubren la misma zona del subsuelo están disponibles cada vez más
para análisis e interpretación. Otro ejemplo son datos sísmicos de
tiempo continuo, en que se adquieren datos sísmicos múltiples veces
sobre un yacimiento en producción. En este caso las diferencias
entre las series de datos sísmicos pueden revelar información
acerca de cambios en límites de fluidos a medida que se producen
hidrocarburos.
La etapa 110 muestra la entrada de datos de
propiedades de las rocas deducidos de sísmicos. Los datos de
propiedades de las rocas deducidos de sísmicos, se pueden usar
directamente para analizar cambios en la litología y los fluidos en
capas. También, se mantiene información acerca de la estructura y
la estratigrafía y, con frecuencia, incluso se exalta en relación
con datos sísmicos. El uso de datos de propiedades de las rocas
deducidos de sísmicos, en análisis e interpretación de subsuelo se
prefiere, por lo tanto, con frecuencia, sobre el uso de datos de
reflexión sísmica o se hace junto con análisis e interpretación del
subsuelo, de datos sísmicos. Por la misma razón, el método de la
materia se aplica preferiblemente a datos de propiedades de las
rocas deducidos de sísmicos.
Se puede utilizar cualquiera de los tipos de
series de datos de propiedades de las rocas, sísmicos y deducidos de
sísmicos, descritos y otras series de datos similares, en el
procedimiento descrito en detalle a continuación. En la práctica,
el procedimiento se aplica a diversas series de datos, ya que cada
serie de datos puede destacar diferentes aspectos de los cambios
espaciales en la estructura, estratigrafía, litología y fluidos, del
subsuelo.
La etapa 120 es para dar entrada al modelo
terrestre para el área del subsuelo de interés. La Fig. 7 muestra un
ejemplo. El modelo terrestre se determina a partir de una
interpretación de horizontes y fallas del subsuelo a partir de los
datos sísmicos, una parte de los cuales se muestra en la figura 1
y/o datos de propiedades de las rocas, deducidos de sísmicos,
correspondientes, una parte de los cuales se muestra en la figura 4
y/o datos de salida a partir del procedimiento de la materia y de
las relaciones geométricas entre estos horizontes y fallas. A
partir de esta información se puede generar un modelo terrestre que
consta de una serie de capas en que se determinan los límites, para
cada capa, a partir de los horizontes y las fallas de entrada y, en
que, en el interior de cada capa, se genera una serie de horizontes
que definen las denominadas microcapas para esa capa. La geometría
de estos horizontes de microcapas se puede modelar de diferentes
maneras para captar el modelo deposicional en el interior de cada
capa del modelo terrestre. En general, se usan los horizontes de
entrada del modelo terrestre, horizontes de referencia definidos
por el usuario o una combinación de estos horizontes, para guiar la
generación de los horizontes de microcapas. En el caso de
deposición compleja, se pueden aplicar métodos más sofisticados para
generar los horizontes de microcapas. La resolución vertical de
estos horizontes de microcapas puede variar, pero requiere que esté
bajo para rodear el intervalo de la muestra de tiempo o
profundidad, vertical, de los datos de propiedades de las rocas,
sísmicos o deducidos de sísmicos, de entrada, para maximizar la
información extraída por el procedimiento de la materia. Los
horizontes de microcapas se pueden definir sobre la misma cuadrícula
x,y que las cuadrículas de propiedades de las rocas, sísmicas o
deducidas de sísmicas, de entrada, si bien esto no es una condición
necesaria. El modelo terrestre descrito puede oscilar de simple, en
el caso de que el objetivo sea captar los cambios espaciales
principales en la estructura, estratigrafía, litología y fluidos,
del subsuelo, usando el procedimiento de la materia o puede ser muy
detallado en el caso de que el objetivo sea captar cambios
espaciales más suaves en la estructura, estratigrafía, litología y
fluidos, del subsuelo. En su manera más simple, el modelo terrestre
estará basado en dos horizontes paralelos con las microcapas en
medio también pasando paralelas a los horizontes de entrada. Se
deducen modelos terrestres muy detallados por interpretación
detallada de los datos sísmicos, datos de propiedades de las rocas
deducidos de sísmicos y datos generados por el procedimiento de la
materia, disponibles. A medida que los modelos terrestres llegan a
ser más detallados, mejora la precisión con que la
microestratificación del modelo terrestre capta la estructura y la
estratigrafía del subsuelo.
La siguiente etapa es generar los datos de salida
a lo largo de los horizontes de microcapas y horizontes de
microcapas rotados, que se rotan localmente para reducir óptimamente
ruido y para captar variaciones espaciales en estructura,
estratigrafía, litología y fluidos, del subsuelo. La salida consta
de diversas determinaciones de la velocidad de cambio de los datos
sísmicos o datos de propiedades de las rocas deducidos de sísmicos,
en relación con los horizontes de microcapas u horizontes de
microcapas rotados localmente, la dirección en que la velocidad de
cambio es la mayor para cada una de estas determinaciones, el
buzamiento y el azimut de rotación correspondientes y la salida del
filtro de filtros aplicados. El algoritmo se presenta en la presente
memoria en la forma más adecuada para exposición, pero se puede
modificar fácilmente en otras formas sin apartarse del alcance de
la invención. Las etapas de cálculo clave para la versión en 3D del
procedimiento, como se aplican a un cubo de propiedades de las
rocas, sísmico o deducido de sísmico, de entrada, son:
Para cada horizonte de microcapas
Definir el horizonte de microcapas seleccionado
como el "horizonte actual" (etapa 200).
Para cada punto de la cuadrícula en el horizonte
actual:
- Definir el punto de la cuadrícula seleccionado como el "punto de la cuadrícula actual", encontrar sus coordenadas (x, y, t o z) asociadas y por interpolación de los datos de entrada de propiedades de las rocas, sísmicos o deducidos de sísmicos, calcular la amplitud en el punto de la cuadrícula actual (etapa 210).
- Encontrar, en base a criterios de distancia especificados por el usuario, típicamente entre 25 y 500 metros, una serie de puntos de la cuadrícula sobre el horizonte actual alrededor del punto de la cuadrícula actual y determinar sus coordenadas (x, y, t o z). Esta serie de puntos, que incluye el punto de la cuadrícula actual, se define como el "segmento del horizonte de microcapas actual" (etapa 220).
- Para cada ángulo en la dirección de la línea alineada, calculado a partir de un intervalo e incremento de ángulos de rotación del segmento del horizonte de microcapas de la línea alineada especificada por el usuario:
- Para cada ángulo en la dirección de la línea transversal, calculado a partir de un intervalo e incremento de ángulos de rotación del segmento del horizonte de microcapas de la línea transversal especificado por el usuario:
- Rotar el segmento del horizonte de microcapas, actual, alrededor del punto de la cuadrícula, actual, de acuerdo con el ángulo de la línea alineada y la línea transversal, actual, para definir el "segmento del horizonte de microcapas, actual, rotado" y las coordenadas (x', y', t' o z') de los puntos de la cuadrícula del segmento del horizonte de microcapas, actual, rotado (etapa 230).
- Por interpolación de los datos de entrada de propiedades de las rocas, sísmicos o deducidos de sísmicos, calcular la amplitud en cada punto del segmento del horizonte de microcapas, actual, rotado, para definir la serie de datos que constan de (a x', y', t' o z') en cada punto, en que a es la amplitud para ese punto calculada por la interpolación (etapa 240).
- A partir del segmento del horizonte de microcapas, actual, rotado, y los datos (a, x', y', t' o z') de los puntos del segmento, calcular el gradiente de las amplitudes (a) o, si el usuario ha definido ponderaciones para los puntos del segmento del horizonte de microcapas, actual, rotado, el gradiente se calcula teniendo en cuenta la ponderación (por ejemplo, un punto que tiene una ponderación w en relación con el punto actual, se cuenta con una ponderación w en el cálculo del gradiente). El gradiente se refiere como "gradiente de la propiedad" para distinguirlo del gradiente geométrico de la superficie definida por los puntos (x', y', t' o z'). El gradiente de la propiedad calculado se refiere como gradiente de la propiedad para el caso tanto no ponderado como el ponderado (etapa 250).
- Calcular el tamaño absoluto del gradiente de la propiedad. El tamaño absoluto del gradiente de la propiedad se almacena como una determinación de la velocidad de cambio (etapa 260).
- Calcular la dirección del gradiente de la propiedad por su proyección sobre el plano x, y y tomando después el ángulo en relación con el norte. El ángulo se almacena (etapa 270).
- Filtrar las amplitudes del segmento del horizonte de microcapas, actual, rotado, por una serie de filtros definidos por el usuario. Los filtros adecuados incluyen un filtro promedio y cualquiera de diversos filtros no lineales tales como un filtro de valor mediano. Las ponderaciones definidas por el usuario se pueden asignar en la operación del filtro similar a la manera en que se usan en el cálculo del gradiente de la propiedad ponderado (etapa 280).
- Después de la filtración, para cada tipo de filtro aplicado, almacenar el valor de salida del filtro. Además, para cada tipo de filtro, calcular la diferencia entre la salida del filtro y el valor (a) de la amplitud en el punto de la cuadrícula actual, y almacenar estos valores como determinaciones adicionales de la velocidad de cambio (etapa 290).
- Almacenar los ángulos de rotación de la línea alineada y la línea transversal con las correspondientes determinaciones de velocidad de cambio, dirección del gradiente de la propiedad y salidas del filtro (etapa 300).
- Terminar el ciclo sobre todos los ángulos de la línea transversal.
- Terminar el ciclo sobre todos los ángulos de la línea alineada.
Para cada determinación de velocidad de
cambio:
Definir la determinación de velocidad de cambio
como la "determinación de la velocidad de cambio actual"
(etapa 400).
(etapa 400).
Recuperar del almacenamiento y recopilar en una
lista, para todas las combinaciones de ángulos de la línea alineada
y la línea transversal, los siguientes datos almacenados para el
punto de la cuadrícula actual: ángulo de la línea alineada, ángulo
de la línea transversal, valor de la determinación de la velocidad
de cambio actual, la dirección del gradiente de la propiedad y, en
caso de filtración, el valor de salida del filtro y la diferencia
entre el valor de salida del filtro y el valor de la amplitud en el
punto de la cuadrícula actual (etapa 410).
A partir de la lista, encontrar la entrada para
la que es mínimo el valor absoluto de la determinación de la
velocidad de cambio actual (etapa 420). A partir del ángulo de la
línea alineada y el ángulo de la línea transversal, que
corresponden con esta entrada, calcular el buzamiento y el azimut.
Este buzamiento y azimut determinan la perturbación del segmento
del horizonte de microcapas, actual, rotado, que corresponde a la
determinación mínima de velocidad de cambio, en relación con el
segmento del horizonte de microcapas actual. En esta entrada,
extraer también el valor de la determinación de la velocidad de
cambio actual, la dirección del gradiente de la propiedad, la
amplitud del punto de la cuadrícula actual y, en caso de
filtración, el valor de salida del filtro y la diferencia entre el
valor de salida del filtro y el valor de la amplitud en el punto de
la cuadrícula actual. Todos estos valores calculados y extraídos se
refieren conjuntamente como los "atributos de salida" (etapa
430).
Almacenar los atributos de salida en el punto de
la cuadrícula actual (etapa 440).
Terminar el ciclo sobre todas las determinaciones
de velocidad de cambio.
Terminar el ciclo sobre todos los puntos de la
cuadrícula del horizonte de microcapas actual.
Terminar el ciclo sobre todos los horizontes de
microcapas.
Para cada determinación de velocidad de
cambio:
Sacar la información almacenada a un fichero
soportado sobre disco del ordenador o soportado en memoria del
ordenador, organizada de manera que a partir de ella se extraigan
los horizontes fácilmente, en que cada horizonte corresponda a un
horizonte de microcapas y tenga asociado con él, en cada punto de
definición del horizonte, los correspondientes atributos de salida.
Esta salida se refiere además como "salida de horizontes de
microcapas" (etapa 600).
Terminar el ciclo sobre todas las determinaciones
de velocidad de cambio.
Para cada determinación de velocidad de
cambio:
Interpolar espacialmente cada atributo de salida
para los puntos de la cuadrícula o subcuadrícula (x, y, t o z) de
los datos de propiedades de las rocas sísmicos o deducidos de
sísmicos, de entrada, usados en los cálculos del procedimiento.
Alternativamente, se usa una cuadrícula definida por el usuario.
Esta salida se refiere además como "salida cuadriculada" (etapa
610).
Sacar la salida cuadriculada a un archivo
soportado sobre disco de ordenador o soportado en la memoria del
ordenador (etapa 620).
Terminar el ciclo sobre todas las determinaciones
de velocidad de cambio.
Este algoritmo se aplica a todos los cubos de las
propiedades de las rocas, sísmicos y deducidos de sísmicos, de
entrada, disponibles.
Las Fig. 8, Fig. 9 y Fig. 10 muestran ejemplos de
salidas del nuevo método. La Fig. 8 muestra el tamaño absoluto del
gradiente de la propiedad sobre los puntos de la cuadrícula de un
horizonte de microcapas. Las amplitudes superiores indicadas por
las flechas muestran claramente una falla. La Fig. 9 muestra el
azimut del gradiente de la propiedad sobre los puntos de la
cuadrícula de un horizonte de microcapas diferente. Las amplitudes
superiores indicadas por las flechas muestran claramente fallas
múltiples. La Fig. 10 muestra el tamaño absoluto del gradiente de
la propiedad sobre los puntos de la cuadrícula de aún otro
horizonte de microcapas que está pasando por unas arenas del
yacimiento. La magnitud del gradiente de la propiedad indica
claramente diferencias en la continuidad lateral de las arenas del
yacimiento. En el área A, la magnitud del gradiente de la propiedad
delimita una arena de yacimiento con buena continuidad lateral,
mientras en el área B, la rápida variación lateral de la magnitud
del gradiente de la propiedad indica que las arenas del yacimiento
presentan continuidad lateral deficiente.
El gradiente de la propiedad especificado
anteriormente se calcula fijando una superficie a través de las
amplitudes de los puntos del segmento del horizonte de microcapas,
actual, rotado, en que se usan las ponderaciones suministradas para
dar énfasis o quitar énfasis a la contribución de ciertos puntos en
el segmento. Están disponibles muchos algoritmos diferentes para
fijación de superficie. La superficie que se tiene que fijar también
puede variar y puede estar en un plano, ranuras bicúbicas u otra
superficie paramétrica. En la realización preferida se ofrece al
usuario una elección de opciones para el algoritmo fijado y el tipo
de superficie. El gradiente de la propiedad se calcula a partir de
la superficie en el punto de definición actual y sale en la forma
del tamaño absoluto del gradiente y el correspondiente azimut.
En la etapa final, los datos de salida
cuadriculados se cargan en estaciones de trabajo de gráficos, con
programas informáticos de interpretación de datos sísmicos
disponibles, clásicos, para análisis e interpretación adicionales
de la información revelada en los datos de salida sobre la
estructura, estratigrafía, litología y fluidos, del subsuelo y los
cambios espaciales en los mismos.
La serie de horizontes de microcapas permite un
método alternativo y novedoso de análisis e interpretación de datos
del subsuelo, en base a la visualización de cualquiera de los
valores de los atributos de cualquiera de los horizontes de
microcapas en una vista de mapa, como se ilustra en la Fig.8, Fig.
9 y Fig. 10 y en la visualización animada generada por la ejecución
secuencialmente a través de una pila de horizontes de microcapas
ordenados en tiempo o profundidad. Si la geometría de las
microcapas se ajusta a la estratigrafía, y asumiendo que se conoce
la posición de los límites estructurales, las variaciones laterales
que están conduciendo, entonces, el procedimiento de la materia sin
rotación angular, destacarán específicamente variaciones en la
litología y el fluido. Para muchos entornos geológicos no es
prácticamente factible o exige mucho tiempo construir un modelo con
tal precisión que las microcapas se ajusten con precisión a la
estratigrafía del subsuelo. El objeto deseable de encontrar
variaciones a lo largo de la estratigrafía se consigue entonces con
el procedimiento de rotación angular descrito. La rotación de las
microcapas compensa las faltas de precisión del modelo de entrada,
a fin de que se pueda calcular la variación lateral real de datos
de propiedades de las rocas, sísmicos y deducidos de sísmicos, a lo
largo de la estratigrafía, con modelos terrestres prácticamente
disponibles. Cuando se visualizan los atributos a lo largo de los
horizontes de microcapas o cuando se ejecuta una pila de horizontes
de microcapas secuencialmente en el tiempo o profundidad y los
atributos se visualizan para cada horizonte de microcapas, se
obtienen nuevas percepciones en procedimientos deposicionales
geológicos y en cambios en la estratigrafía y litología que no se
revelan en el análisis e interpretación de la salida de atributos
cuadriculada usando métodos de análisis e interpretación
sísmicos,
clásicos.
clásicos.
Es fácilmente obvio cómo se puede aplicar el
algoritmo anterior a datos de propiedades de las rocas, sísmicos y
deducidos de sísmicos, en 2D. Esto se consigue por eliminación del
ciclo de perturbación angular de la línea transversal. También es
fácilmente obvio que se puedan generar datos adicionales valiosos
aplicando las operaciones de filtración descritas cuando la
determinación de la velocidad de cambio actual es el gradiente de la
propiedad. Los filtros se aplican entonces al segmento del
horizonte de microcapas, actual, rotado, para el que es mínimo el
valor absoluto del gradiente de la propiedad y las salidas del
filtro se almacenan como adiciones adicionales de los atributos de
salida. Además, lo anterior describe las determinaciones de
velocidad de cambio procedentes del gradiente de la propiedad, con
o sin ponderaciones, y operaciones de filtración, con o sin
ponderaciones. Se pueden deducir otras diversas determinaciones de
cambio de velocidad como una función de las amplitudes y
ponderaciones del segmento del horizonte de microcapas, actual, y
el uso de cualquier determinación de velocidad de cambio procedente
de tales funciones se considera parte del método.
El procedimiento hacer salir a lo largo de cada
horizonte de microcapas los valores de los datos de propiedades de
las rocas, sísmicos o deducidos de sísmicos, de entrada, obtenidos
por interpolación. Esto permite un método de análisis e
interpretación en que el usuario ejecuta consecutivamente la serie
de horizontes de microcapas y en que en el mapa o modo de vista en
3D, se contornea el horizonte o se da color codificado con estos
valores de amplitud o los valores de cualquiera de los atributos de
salida. En el caso de que los horizontes de microcapas se pudieran
ajustar exactamente a un horizonte dado, los resultados serían los
mismos que los que se obtendrían con una técnica aplicada
comúnmente en que las amplitudes se extraen justo a lo largo de un
horizonte que se copia hacia arriba y/o hacia abajo en un intervalo
especificado por el usuario. Con este nuevo análisis y método de
interpretación se obtienen resultados superiores cuando la
estructura y la estratigrafía del subsuelo son de manera que los
horizontes de microcapas no se pueden ajustar más tiempo a uno de
los horizontes de entrada.
Si se aplica un filtro, el procedimiento también
produce la salida del filtro como datos para uso adicional en
análisis e interpretación. Los datos de propiedades de las rocas,
sísmicos o deducidos de sísmicos, siempre contienen ruido. El
procedimiento de filtración por horizontes de microcapas proporciona
un nuevo método de reducción de ruido. La base para esto es que los
puntos de los datos a lo largo del horizonte de microcapas tenderán
a ser más parecidos que si se pasa un filtro en una dirección que
no se ajusta a la estratigrafía. Como resultado, el
"procedimiento de filtración estratigráfica" puede revelar
características escondidas en el ruido que puede que no se recuperen
con filtros aplicados clásicos, que no tienen en cuenta la
geometría local. El procedimiento de reducción de ruido se exalta
además por el uso del procedimiento de rotación del horizonte de
microcapas. En el procedimiento de rotación cada segmento del
horizonte de microcapas (el segmento del horizonte de microcapas
actual) usado en la operación de filtración se rota para encontrar
los puntos que, de acuerdo con la determinación aplicada, sean tan
parecidos como sea posible. Aplicar un filtro a los puntos del
segmento del horizonte de microcapas, rotado, que corresponden al
valor óptimo de la determinación aplicada, exalta además la
potencia de reducción de ruido del filtro en relación con el uso del
mismo filtro a lo largo de horizontes de microcapas, ya que estos
pueden contener errores en su orientación con respecto a la
estratigrafía.
Dos variaciones importantes de los métodos
anteriores son como sigue. La primera variación es que se implante
el procedimiento de perturbación angular óptimo con un procedimiento
de optimización en lugar del procedimiento de barrido descrito. La
función objetivo es el valor absoluto de la determinación de
velocidad de cambio lateral, actual, entonces y el procedimiento de
optimización pretenderá encontrar los ángulos de la línea alineada y
de la línea transversal bajo los cuales se minimice la función
objetivo. Tales procedimientos de optimización se describen bien y
están fácilmente disponibles para su implantación.
La segunda variación es en el caso de que se
reemplace el procedimiento de interpolación para calcular amplitudes
a partir de las cuadrículas de propiedades de las rocas, sísmicas o
deducidas de sísmicas, de entrada, por un procedimiento en que las
amplitudes se calculan por promedio alrededor de la microcapa
actual. Una realización es tomar como intervalo promedio el
intervalo vertical desde la mitad entre el horizonte de microcapas
actual y el horizonte de microcapas siguiente por encima a la mitad
de la microcapa actual y la microcapa siguiente por debajo. Otra
realización es tomar como intervalo promedio el intervalo vertical
a bien, el horizonte de microcapas que se encuentra siguiente por
encima o que se encuentra siguiente por debajo, para crear
promedios sobre el espesor de las microcapas. Las operaciones
promediadas proporcionan robustez frente al ruido, exaltan
características con pequeña variación vertical en relación con las
de con variación vertical mayor y permiten trabajar con valores
promedio de propiedades de las rocas, sísmicos o deducidos de
sísmicos, en base al espesor de microcapa.
La presente invención no se limita a la
realización descrita y a las variaciones descritas anteriormente;
los derechos solicitados se determinan por las siguientes
reivindicaciones, dentro del alcance de las cuales se pueden prever
muchas modificaciones. Se desea cubrir por las reivindicaciones
adjuntas todas dichas modificaciones implicadas dentro del alcance
de las reivindicaciones.
Claims (24)
1. Método para determinar y analizar cambios
espaciales en el subsuelo de la tierra y para reducir ruido en datos
sísmicos, que comprende las etapas de:
- a)
- obtener datos sísmicos;
- b)
- obtener un modelo terrestre del subsuelo con capas representativas de la estructura del subsuelo;
- c)
- definir una serie de microcapas definidas por una serie de horizontes de microcapas entre los límites de las capas del modelo terrestre;
- d)
- determinar puntos de la cuadrícula de los horizontes de microcapas a lo largo de un horizonte de microcapas;
- e)
- determinar las amplitudes de los datos sísmicos en los puntos de la cuadrícula de los horizontes de microcapas;
- f)
- determinar en cada punto de la cuadrícula de los horizontes de microcapas una o más determinaciones representativas de la velocidad de cambio de las amplitudes de los datos sísmicos y/o determinaciones de la dirección en que cambian las amplitudes de los datos sísmicos y/o determinar en cada punto de la cuadrícula de los horizontes de microcapas una o más salidas del filtro;
- g)
- almacenar como salida dicha amplitud de los datos sísmicos, dichas determinaciones y dichas salidas del filtro en cada punto de la cuadrícula de los horizontes de microcapas para análisis adicional del subsuelo;
2. Método de acuerdo con la reivindicación 1, en
el que se repiten las etapas d-g para todos los
horizontes de microcapas en el modelo terrestre.
3. Método de acuerdo con la reivindicación 1 ó 2,
en que la etapa f) comprende:
- f1)
- definir para un punto de la cuadrícula de los horizontes de microcapas un segmento de la superficie de microcapas que consta de una serie de puntos de la cuadrícula de los horizontes de microcapas sobre el horizonte de microcapas que rodea y puede incluir dicho punto de la cuadrícula de los horizontes de microcapas;
- f2)
- determinar las amplitudes de los datos sísmicos, en los puntos del segmento;
- f3)
- determinar a partir de las amplitudes de los datos sísmicos, en los puntos del segmento, una o más determinaciones representativas de la velocidad de cambio de las amplitudes de los datos sísmicos, sobre ese punto de la cuadrícula de los horizontes de microcapas y/o una determinación representativa de la dirección en que cambian las amplitudes de los datos sísmicos.
4. Método de acuerdo con la reivindicación 3, que
comprende, después de la etapa f3, la etapa de filtración de las
amplitudes de los datos sísmicos en los puntos del segmento con uno
o más filtros para generar uno o más valores de amplitud filtrada
en el punto de la cuadrícula de los horizontes de microcapas.
5. Método de acuerdo con la reivindicación 3 ó 4,
que comprende además:
- f4)
- rotar la serie de puntos del segmento alrededor del punto de la cuadrícula de los horizontes de microcapas a una serie rotada de puntos del segmento;
- f5)
- determinar las amplitudes de los datos sísmicos en la serie rotada de puntos del segmento;
- f6)
- determinar a partir de las amplitudes de los datos sísmicos, en los puntos del segmento rotado, una o más determinaciones representativas de la velocidad de cambio de las amplitudes de los datos sísmicos en ese punto de la cuadrícula de los horizontes de microcapas y una determinación representativa de la dirección en que cambian las amplitudes de los datos sísmicos y determinar a partir de la rotación determinaciones adicionales representativas de la dirección en que cambian las amplitudes de los datos sísmicos;
- f7)
- filtrar las amplitudes de los datos sísmicos en los puntos del segmento rotado con uno o más filtros para generar uno o más valores de amplitud filtrada en el punto de la cuadrícula de los horizontes de microcapas.
- f8)
- determinar para cada rotación el valor absoluto de una determinación representativa seleccionada de la velocidad de cambio de las amplitudes de los datos sísmicos;
- f9)
- determinar la rotación para la que es mínimo el valor absoluto de dicha determinación;
- f10)
- seleccionar para salida en el punto de la cuadrícula de los horizontes de microcapas la amplitud de los datos sísmicos y las determinaciones y salidas del filtro calculadas para dicha rotación, para la cual es mínimo el valor absoluto de dicha determinación.
6. Método de acuerdo con cualquiera de las
reivindicaciones 1-5, en el que se determina la
amplitud de los datos sísmicos en el punto de la cuadrícula de los
horizontes de microcapas y los puntos del segmento, por
interpolación de amplitudes de los datos sísmicos de puntos de la
cuadrícula vecinos de la cuadrícula sísmica.
7. Método de acuerdo con cualquiera de las
reivindicaciones 1-5, en el que se determina la
amplitud de los datos sísmicos en el punto de la cuadrícula de los
horizontes de microcapas y los puntos del segmento, por cálculo del
promedio de las amplitudes de los datos sísmicos de puntos de la
cuadrícula vecinos de la cuadrícula sísmica.
8. Método de acuerdo con la reivindicación 7, en
el que el cálculo del promedio se lleva a cabo tomando el intervalo
vertical desde la mitad entre el horizonte de microcapas actual y
el horizonte de microcapas siguiente, por encima, a la mitad de la
microcapa actual y la microcapa siguiente, por debajo.
9. Método de acuerdo con la reivindicación 7, en
el que el cálculo del promedio se lleva a cabo tomando como
intervalo el intervalo vertical para el horizonte de microcapas,
bien siguiente que se encuentra por encima o siguiente que se
encuentra por debajo, para crear promedios sobre el espesor de las
microcapas.
10. Método de acuerdo para cualquiera de las
reivindicaciones 1-5, en el que una de las
determinaciones representativas de la velocidad de cambio de las
amplitudes sísmicas es la magnitud del gradiente determinado a
partir de las amplitudes de los datos sísmicos en los puntos del
segmento y en el que se determina una de las determinaciones
representativas de la dirección en que cambian las amplitudes
sísmicas a partir de la proyección del gradiente sobre un plano
horizontal y tomando el ángulo en relación con una dirección del
compás predeterminada.
11. Método de acuerdo con la reivindicación 10,
en el que se ponderan el gradiente y la dirección de los cálculos
del gradiente, de acuerdo con las ponderaciones asignadas a los
puntos del segmento sobre el que se basa el cálculo del
gradiente.
12. Método de acuerdo con cualquiera de las
reivindicaciones 1-4, en el que se obtienen
determinaciones representativas adicionales de la velocidad de
cambio de las amplitudes sísmicas, filtrando primero los datos de la
amplitud sísmica de los puntos del segmento con un filtro lineal o
no lineal, para determinar el valor filtrado en el punto de
definición, seguido por su sustracción de la amplitud de los datos
sísmicos en el punto de definición.
13. Método de acuerdo con la reivindicación 12,
en que se ponderan los cálculos del filtro de acuerdo con
ponderaciones asignadas a los puntos del segmento usados en la
filtración.
14. Método de acuerdo con cualquiera de las
reivindicaciones 1-13, que comprende: determinar
determinaciones adicionales representativas de la dirección en que
cambian las amplitudes de los datos sísmicos a partir de la
rotación por cálculo del ángulo de buzamiento de la dirección de
rotación en relación con la vertical y por el cálculo del ángulo
azimutal de la dirección de rotación por proyección de la dirección
de rotación sobre un plano horizontal y tomar el ángulo en relación
con una dirección del compás predeterminada.
15. Método de acuerdo con cualquiera de las
reivindicaciones 1-14, en el que las rotaciones se
llevan a cabo por la rotación paso a paso de la serie de puntos del
segmento sobre una pluralidad de ángulos de rotación en la dirección
de la línea para datos sísmicos en 2D o la dirección de la línea
alineada y la línea transversal u otra serie de direcciones no
paralelas para datos sísmicos en 3D, en que para cada dirección se
determinan los ángulos de rotación a partir de un valor de partida
del ángulo de rotación, valor final del ángulo de rotación e
incremento del ángulo, definido por el usuario.
16. Método de acuerdo con cualquiera de las
reivindicaciones 1-15, en el que se usa un
procedimiento de optimización para encontrar una serie de
rotaciones de manera que se encuentre eficazmente el mínimo del
valor absoluto de dicha determinación.
17. Método de acuerdo con cualquiera de las
reivindicaciones 1-16, en el que se almacena cada
elemento de los datos de salida en su correspondiente horizonte de
microcapas, en su correspondiente punto de la cuadrícula de los
horizontes de microcapas, de manera que para cada horizonte de
microcapas seleccionado los elementos de los datos de salida
correspondientes estén disponibles independientemente.
18. Método de acuerdo con cualquiera de las
reivindicaciones precedentes, en el que para cualquier horizonte de
microcapas, se recupera cualquier elemento de datos de salida
almacenado en ese horizonte de microcapas y se representa
visualmente para una parte o todo del horizonte de microcapas, para
analizar e interpretar la variación espacial de la estructura,
estratigrafía, litología y contenido de fluido, del subsuelo.
19. Método de acuerdo con cualquiera de las
reivindicaciones precedentes, en el que para cualquier pila de
horizontes de microcapas organizados secuencialmente en tiempo o
profundidad, se recupera cualquier elemento de datos de salida
almacenado en dicha pila de horizontes de microcapas y se representa
visualmente, secuencialmente, para una parte o todo de los
horizontes de microcapas, horizonte de microcapas por horizonte de
microcapas, para analizar e interpretar la variación espacial de la
estructura, estratigrafía, litología y distribución de fluidos, del
subsuelo.
20. Método de acuerdo con cualquiera de las
reivindicaciones precedentes, en el que para cualquier pila de
horizontes de microcapas se recupera cualquier elemento de datos de
salida almacenado en dicha pila de horizontes de microcapas y se
vuelve a cuadricular en un cuadrícula regular que puede ser la misma
que la cuadrícula de los datos sísmicos de entrada.
21. Método de acuerdo con cualquiera de las
reivindicaciones precedentes, en el que se usa en lugar de datos
sísmicos, datos de propiedades de las rocas deducidos de
sísmicos.
22. Método de acuerdo con cualquiera de las
reivindicaciones precedentes, en el que el procedimiento se repite
para las diferentes determinaciones alternativas para la velocidad
de cambio, para diferentes definiciones de los puntos del segmento
que rodean a cada punto de definición de los horizontes de
microcapas, y para diferentes ponderaciones de los puntos del
segmento para generar una pluralidad de series de datos de salida
que, cada uno de una manera diferente, destaque cambios espaciales
en la geometría, estratigrafía, litología y distribución de fluido,
del subsuelo.
23. Método de acuerdo con cualquiera de las
reivindicaciones precedentes, en el que el procedimiento se aplica
iterativamente, en que en cada etapa se redefine el modelo
terrestre, y en que tales redefiniciones del modelo terrestre están
basadas en la interpretación y análisis de datos de salida del
procedimiento generados sobre la base de modelos terrestres
precedentes.
24. Método de acuerdo con cualquiera de las
reivindicaciones precedentes, en el que el método se expresa como un
programa de ordenador.
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EP99201782A EP1058131B1 (en) | 1999-06-03 | 1999-06-03 | Method of determining spatial changes in subsurface structure stratigraphy, lithology and fluid content and of reducing seismic noise |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| ES2223156T3 true ES2223156T3 (es) | 2005-02-16 |
Family
ID=8240275
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| ES99201782T Expired - Lifetime ES2223156T3 (es) | 1999-06-03 | 1999-06-03 | Metodo de determinacion de cambios espaciales en la estratigrafia de la estructura del subsuelo, litologia y contenido de fluidos y de reduccion del ruido sismico. |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US6401042B1 (es) |
| EP (1) | EP1058131B1 (es) |
| DE (1) | DE69919848T2 (es) |
| DK (1) | DK1058131T3 (es) |
| ES (1) | ES2223156T3 (es) |
Families Citing this family (27)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2772136B1 (fr) * | 1997-12-08 | 2000-01-07 | Elf Exploration Prod | Methode de detection d'une discontinuite geologique presente dans un milieu par utilisation du flux optique |
| GB2372567B (en) * | 2001-02-22 | 2003-04-09 | Schlumberger Holdings | Estimating subsurface subsidence and compaction |
| US7248259B2 (en) * | 2001-12-12 | 2007-07-24 | Technoguide As | Three dimensional geological model construction |
| US6885942B2 (en) * | 2003-01-09 | 2005-04-26 | Schlumberger Technology Corporation | Method to detect and visualize changes in formation parameters and borehole condition |
| US7206782B1 (en) * | 2003-01-29 | 2007-04-17 | Michael John Padgett | Method for deriving a GrAZ seismic attribute file |
| US7013218B2 (en) * | 2003-07-16 | 2006-03-14 | Siesmic Micro-Technology, Inc. | System and method for interpreting repeated surfaces |
| US7463552B1 (en) * | 2003-11-25 | 2008-12-09 | Michael John Padgett | Method for deriving 3D output volumes using filters derived from flat spot direction vectors |
| CA2572449C (en) | 2004-07-01 | 2014-05-06 | Exxonmobil Upstream Research Company | Hydrodynamics-based gridding geologic modeling (hydro-gridding) |
| US7492968B2 (en) * | 2004-09-07 | 2009-02-17 | Siemens Medical Solutions Usa, Inc. | System and method for segmenting a structure of interest using an interpolation of a separating surface in an area of attachment to a structure having similar properties |
| US7373251B2 (en) * | 2004-12-22 | 2008-05-13 | Marathon Oil Company | Method for predicting quantitative values of a rock or fluid property in a reservoir using seismic data |
| EP1707993B1 (fr) * | 2005-03-29 | 2009-08-19 | Total S.A. | Procédé et programme de recherche de discontinuites géologiques |
| RU2009136585A (ru) * | 2007-03-05 | 2011-04-10 | Пэредайм Джиофизикал (Люксембург) С.А.Р.Л. (Lu) | Основанная на модели томография с сохранением времени |
| US8310481B2 (en) * | 2007-10-12 | 2012-11-13 | Edward Ernest Bailey | Computer aided design method for enhancement of local refinement through T-splines |
| GB2468088B (en) * | 2007-11-27 | 2012-08-15 | Exxonmobil Upstream Res Co | Method for determining the properties of hydrocarbon reservoirs from geophysical data |
| CA2705277C (en) * | 2007-12-18 | 2017-01-17 | Exxonmobil Upstream Research Company | Determining connectivity architecture in 2-d and 3-d heterogeneous data |
| AU2008340399B2 (en) | 2007-12-21 | 2013-09-26 | Exxonmobil Upstream Research Company | Method and apparatus for analyzing three-dimensional data |
| EP2252903A4 (en) | 2008-03-10 | 2018-01-03 | Exxonmobil Upstream Research Company | Method for determing distinct alternative paths between two object sets in 2-d and 3-d heterogeneous data |
| AU2009244721B2 (en) | 2008-05-05 | 2013-09-26 | Exxonmobile Upstream Research Company | Systems and methods for connectivity analysis using functional obejects |
| US8094515B2 (en) * | 2009-01-07 | 2012-01-10 | Westerngeco L.L.C. | Seismic data visualizations |
| CN101852866B (zh) * | 2010-05-13 | 2012-01-04 | 中国石油天然气股份有限公司 | 一种叠后地震数据滤波方法 |
| RU2573166C2 (ru) | 2010-05-28 | 2016-01-20 | Эксонмобил Апстрим Рисерч Компани | Способ сейсмического анализа углеводородных систем |
| US20130138350A1 (en) * | 2011-11-24 | 2013-05-30 | Manoj Vallikkat Thachaparambil | Enhanced termination identification function based on dip field generated from surface data |
| US10108762B2 (en) * | 2014-10-03 | 2018-10-23 | International Business Machines Corporation | Tunable miniaturized physical subsurface model for simulation and inversion |
| CN107728204B (zh) * | 2016-08-11 | 2019-03-29 | 中国石油化工股份有限公司 | 基于叠前纵波各向异性的裂缝预测方法及系统 |
| CN108038881B (zh) * | 2018-01-11 | 2018-12-18 | 山东科技大学 | 一种利用平行光评价岩石结构面起伏程度的计算方法 |
| US11156744B2 (en) | 2019-01-10 | 2021-10-26 | Emerson Paradigm Holding Llc | Imaging a subsurface geological model at a past intermediate restoration time |
| US10520644B1 (en) | 2019-01-10 | 2019-12-31 | Emerson Paradigm Holding Llc | Imaging a subsurface geological model at a past intermediate restoration time |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4849887A (en) * | 1987-08-28 | 1989-07-18 | Amoco Corporation | Horizon velocity analysis |
| US5056066A (en) * | 1990-06-25 | 1991-10-08 | Landmark Graphics Corporation | Method for attribute tracking in seismic data |
| GB2278920B (en) * | 1993-06-07 | 1996-10-30 | Geco As | Method of determining earth elastic parameters in anistropic media |
| GB2322704B (en) * | 1994-07-07 | 1998-12-09 | Geco As | Method of Processing seismic data |
| US5930730A (en) * | 1994-12-12 | 1999-07-27 | Amoco Corporation | Method and apparatus for seismic signal processing and exploration |
| CN1186647C (zh) * | 1996-04-12 | 2005-01-26 | 环球核心实验室有限公司 | 处理和探测地震信号的方法及装置 |
| FR2765692B1 (fr) * | 1997-07-04 | 1999-09-10 | Inst Francais Du Petrole | Methode pour modeliser en 3d l'impedance d'un milieu heterogene |
-
1999
- 1999-06-03 ES ES99201782T patent/ES2223156T3/es not_active Expired - Lifetime
- 1999-06-03 EP EP99201782A patent/EP1058131B1/en not_active Expired - Lifetime
- 1999-06-03 DK DK99201782T patent/DK1058131T3/da active
- 1999-06-03 DE DE69919848T patent/DE69919848T2/de not_active Expired - Lifetime
-
2000
- 2000-06-02 US US09/586,047 patent/US6401042B1/en not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| EP1058131A1 (en) | 2000-12-06 |
| DK1058131T3 (da) | 2004-09-20 |
| DE69919848D1 (de) | 2004-10-07 |
| US6401042B1 (en) | 2002-06-04 |
| DE69919848T2 (de) | 2005-09-15 |
| EP1058131B1 (en) | 2004-09-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| ES2223156T3 (es) | Metodo de determinacion de cambios espaciales en la estratigrafia de la estructura del subsuelo, litologia y contenido de fluidos y de reduccion del ruido sismico. | |
| US8392163B2 (en) | Virtual petroleum system with salt restoration functionality | |
| US10013800B1 (en) | Systems and methods for coordinated editing of seismic data in dual model | |
| CA2279266C (en) | Method for determining barriers to reservoir flow | |
| US8793113B2 (en) | Method and apparatus for near well structural modeling based on borehole dips | |
| EP2778725B1 (en) | Systems and methods to build sedimentary attributes | |
| US20160124116A1 (en) | Generation of structural elements for subsurface formation using stratigraphic implicit function | |
| US20090295792A1 (en) | Virtual petroleum system | |
| CN105388526B (zh) | 一种陆相湖盆的古地质结构恢复方法 | |
| US20090299709A1 (en) | Virtual petroleum system | |
| EP3293552B1 (en) | System and method for editing geological models by switching between volume-based models and surface-based structural models augmented with stratigraphic fiber bundles | |
| US20140222403A1 (en) | Geologic model via implicit function | |
| US20160131800A1 (en) | Modeling fluid-conducting fractures in reservoir simulation grids | |
| BR112012009045B1 (pt) | métodos para gerar uma grade de simulação tridimensional para um modelo de reservatório, aparelho de gradeamento de simulação e simulador de reservatório | |
| CN111624651B (zh) | 基于古地貌约束的储层预测方法及装置 | |
| NO335056B1 (no) | Metode for beregning av 3D gittermodeller av et reservoar | |
| EP2778724A2 (en) | Systems and methods for building axes, co-axes and paleo-geographic coordinates related to a stratified geological volume | |
| Shaban et al. | Comparison between curvature and 3D strain analysis methods for fracture predicting in the Gachsaran oil field (Iran) | |
| Sala et al. | The Alpstein in three dimensions: fold-and-thrust belt visualization in the Helvetic zone, eastern Switzerland | |
| CN114814944A (zh) | 基于散度和旋度的弹性纵横波场分离和逆时偏移成像方法 | |
| US10209380B2 (en) | Methods and systems for juxtaposition across geological discontinuities | |
| US20090299703A1 (en) | Virtual petroleum system | |
| CN104111476A (zh) | 构建地层速度场的方法及装置 | |
| CN109324354A (zh) | 一种确定砂体连通性的方法、装置及系统 |