MX2012002684A

MX2012002684A - Sistema y metodo de modelado de formacion de hidrocarbono.

Info

Publication number: MX2012002684A
Application number: MX2012002684A
Authority: MX
Inventors: John Killough
Original assignee: Landmark Graphics Corp
Priority date: 2009-09-02
Filing date: 2009-09-02
Publication date: 2012-06-13
Also published as: AU2009351922B2; CN102576374A; CA2770602A1; AU2009351922A1; US20120179439A1; US9390207B2; EP2457187B1; CA2770602C; WO2011028197A1; BR112012004477A2; EP2457187A4; EP2457187A1; CN102576374B

Abstract

Modelado de formación de hidrocarbono. Al menos algunas de las modalidades ilustrativas son métodos incluyendo la simulación de la reacción de la formación a la extracción de hidrocarbonos durante una pluralidad de pasos de tiempo. Con cada paso de tiempo, el método ilustrativo hace una primera aproximación de migración de saturaciones para al menos un bloque de retícula de la pluralidad de bloques de retícula (en donde la migración de la saturación dentro de al menos un paso de tiempo, abarca más de un bloque de retícula), y después corrige la primera aproximación de migración de saturaciones para un efecto no considerado en la primera aproximación.

Description

SISTEMA Y MÉTODO DE MODELADO DE FORMACIÓN DE HIDROCARBONO ANTECEDENTES Para · maximizar la producción de hidrocarbono a partir de los depósitos de hidrocarbono, las empresas petroleras y de gas estimulan las técnicas de extracción de depósitos usando modelos de depósito, y después implementan la extracción real con base en los resultados identificados. La complejidad y precisión del modelado de depósitos ha incrementado a medida que la tecnología de computadoras ha avanzado, y a medida que las técnicas de modelado de depósitos ha mejorado.

En el medio de modelado de depósitos, existen ventajas recíprocas entre la precisión de modelado de depósitos y velocidad de simulaciones de operación usando el modelo de depósito. Los modelos de depósitos más precisos son más complejos y toman más tiempo en producir resultados. Los modelos de depósitos menos complejos pueden producir resultados más rápidamente o con menos costo computacional, pero pueden no tomar en cuenta adecuadamente las acciones y reacciones geofísicas. Por lo tanto, cualquier técnica que desempeñe más precisamente y más rápidamente el modelado de depósitos es altamente codiciado.

BREVE DESCRIPCIÓN DE LOS DIBUJOS Para una descripción detallada de modalidades ejemplares, ahora se hará referencia a los dibujos que se acompañan en los que: La Figura 1 muestra un método de acuerdo con al menos algunas modalidades; y La Figura 2 muestra un sistema de cómputo de acuerdo con al menos algunas modalidades.

NOTACIÓN Y NOMENCLATURA Ciertos términos se usan a lo largo de la siguiente descripción y reivindicaciones para referirse a componentes de sistema en particular. Como apreciará un experto en la materia, las empresas petroleras y de gas pueden referirse a un componente mediante nombres diferentes. Este documento no pretende distinguir entre componentes que difieren en nombre pero no en su función.

En la siguiente discusión y en las reivindicaciones, los términos "incluyendo", "que incluye", "comprendiendo" y "que comprende" se usan en un modo abierto, y por lo tanto deberá ser interpretado para significar "incluyendo, pero no limitado a...". También, el término "acopla" o "que se acopla" pretende significar ya sea una conexión indirecta o directa. Por lo tanto, si un primer dispositivo se acopla a un segundo dispositivo, esa conexión puede ser a través de una conexión directa o a través de una conexión indirecta mediante otros dispositivos y conexiones.

"Saturaciones" querrá decir la proporción relativa de componentes modelados. Por lo tanto, las saturaciones pueden ser saturaciones de agua e hidrocarbono dentro del volumen modelado, o saturaciones de diferentes hidrocarbonos dentro del volumen modelado.

DESCRIPCIÓN DETALLADA La siguiente discusión está dirigida a varias modalidades de la invención. Aunque una o más de estas modalidades pueden ser preferidas, las modalidades divulgadas no deberán ser interpretadas, o usadas de otra manera, como limitantes del alcance de la divulgación, incluyendo las reivindicaciones. Además, un experto en la materia entenderá que la siguiente descripción tiene una amplia aplicación, y la discusión de cualquier modalidad solo pretende ser ejemplar de dicha modalidad, y no pretende intimar que el alcance de la divulgación, incluyendo las reivindicaciones, está limitada a dicha modalidad .

Las técnicas de la materia relacionadas para el modelaje de formaciones de hidrocarbonos tienen limitaciones inherentes. Por ejemplo, una técnica para modelaje de depósitos es conocida como la técnica de diferencia finita. Las técnicas de diferencias finitas modelan los depósitos como una pluralidad de bloques de retícula de tamaño particular. Las ecuaciones diferenciales que predicen la presión de hidrocarbonos y/o agua dentro de cada bloque de retícula son resueltas. Con base en las presiones calculadas, la velocidad de flujo de fluido en cada cara de cada bloque de retícula es calculada. Sin embargo, la técnica de diferencia finita está limitada en el sentido de que el modelo no puede contar fácilmente para un flujo de hidrocarbono y/o agua que atraviesa más de un bloque de retícula dentro de un periodo modelado (ej . , paso de tiempo) . Dependiendo del tamaño del bloque de retícula y la velocidad del movimiento del fluido, el paso de tiempo para la técnica de diferencia finita puede ser limitado a un tamaño extremadamente pequeño (e . , un día o menos). Al modelar la extracción de depósitos durante la vida de un depósito, que puede ser del orden de 20 años o más, los pasos de tiempo en el orden de un día o menos pueden ser excesivamente pequeños.

Otra técnica de modelaje de depósitos, que no tiene la limitación de paso de tiempo corto, es la técnica de flujo laminar (también conocida como la técnica de Euler-Lagrangian) . La técnica optimizada inicialmente usa una técnica de tipo de diferencia finita para determinar las presiones y el fluido de flujo (ej . , velocidades) en los limites de bloque de retícula (la porción de Euler) , pero después usa las velocidades para interpolar tridimensionalmente el flujo de fluido a través de muchos bloques de retículas (la porción de Lagrangian) . Por ejemplo, los flujos de fluido interpolados pueden fluir muchos bloques de retículas durante el periodo modelado, de ahí el término "flujo laminar". El paso de tiempo para la técnica de flujo laminar puede abarcar periodos de tiempo significativamente más largos (ej ., 60, 90, 180 días), y por lo tanto, puede modelar más rápidamente la reacción del depósito a técnicas de extracción particulares. Sin embargo, para poder usar la técnica de flujo laminar, se ignora algo de la física de flujo de fluido. Por ejemplo, la técnica de flujo laminar no cuenta fácilmente para; gravedad; cambios en la permeabilidad relativa como cambios de saturación de agua; cómo afecta la presión capilar en el flujo de fluido en el medio poroso; o el flujo del fluido transversal al flujo laminar (flujo transversal) .

Una técnica para tomar en cuenta al menos algo de la física del flujo del fluido junto con la técnica de flujo laminar es conocida como "división de operador". La división de operador puede ser conceptualizada como un proceso de dos pasos; primero la técnica tradicional de flujo laminar es aplicada y el fluido "movido" a lo largo del flujo laminar en el modelo; y después la física del flujo del fluido (tal como la gravedad) es aplicada al fluido estacionario en la nueva ubicación. En este respecto, la división de operador es referida como una técnica "explícita", significando que las soluciones a las ecuaciones relacionadas con la física del flujo del fluido se resuelven de manera secuencial, en lugar de simultáneamente, con las ecuaciones relacionadas con la técnica de flujo laminar. Mientras que la división de operador para considerar la física del flujo del fluido de otra manera ignorada por la técnica de flujo laminar mejora la exactitud del modelo, el incremento en la exactitud está limitado .

En un documento de la Society of Petroleum Engineers (Sociedad de Ingenieros Petroleros) titulado "Timestep Selection During Streamline Simulation Through Transverse Flux Correction", los autores, Osaka, Datta-Gupta y King, describen el desempeño de la técnica de flujo laminar que implícitamente considera el flujo transversal. Sin embargo, al añadir consideración del flujo transversal, la técnica de Osaka et al., se vuelve limitada en el tamaño del paso de tiempo. Para que el sistema Osaka et al. permanezca numéricamente estable, el paso de tiempo debe ser seleccionado de forma tal que "la onda más rápida no debe pasar a través de toda una celda durante un paso de tiempo". Por lo tanto, aunque Osaka et al. discute una técnica de flujo laminar "mejorada", un beneficio principal de la técnica de flujo laminar - la capacidad de usar grandes pasos de tiempo - se pierde.

Las varias modalidades están dirigidas a sistemas y métodos, junto con instrucciones de almacenamiento de medios de almacenamiento legibles por computadora, que desempeñan el modelado de depósitos con los beneficios tanto de tomar en cuenta implícitamente el fenómeno físico tal como la permeabilidad relativa y la presión capilar, como también la capacidad del uso de grandes pasos de tiempo. La descripción dará primero un vistazo en palabras, seguido por un tratamiento más matemático.

Las varias modalidades están dirigidas a dividir lógicamente la formación en una pluralidad de volúmenes, o bloques de retícula. En modalidades particulares, el número de bloques de retícula puede estar en el orden de millones de bloques de retícula, pero tales bloques de retícula pueden ser más o menos equivalentemente usados. En algunas modalidades, los bloques de retícula son del mismo volumen con base en la actividad del movimiento de hidrocarbonos y/o agua dentro del bloque de retícula. Por ejemplo, se pueden usar retículas de bloque más pequeñas en áreas "activas", en donde se pueden usar bloques de retícula más grandes en áreas con muy poco o nada de movimiento de fluidos .

Para cada retícula de bloque, y tomando en cuenta las presiones inherentes de formación así como la presión asociada con las fuentes (e . , pozos de inyección) y disipadores (ej., pozos de producción), la presión de los fluidos en cada límite de bloque de retícula es calculada. En modalidades particulares, la presión se calcula usando la técnica de diferencias finitas (ej . técnica Euleriana) . Con base en las presiones en cada límite de bloque de retícula, o más precisamente las diferencias en las presiones consideradas a través de los límites de bloque de retícula, se determinan las velocidades de flujo.

Con base en las velocidades de flujo, la progresión de las saturaciones de' fluido (o masa) se determina durante el paso de tiempo. Dicho de otra manera, las saturaciones (o masas) en cada retícula de bloque al final del paso de tiempo son determinadas. En modalidades particulares, la determinación del progreso de las saturaciones usa la técnica Lagrangiana, y por lo tanto, los fluidos modelados pueden "fluir" a través de múltiples bloques de retícula. Sin embargo, como se comentó anteriormente, la técnica Lagrangiana de este paso del proceso, no cuenta para muchas propiedades físicas del fluido que causan exactitud de la saturación de agua calculada en cada bloque de retícula. Por ejemplo, la saturación de agua calculada no toma en cuenta: gravedad; cambios en la permeabilidad relativa como cambios de saturación de agua; cómo afecta la presión capilar al flujo de fluido en el medio poroso; o el flujo del fluido transversal a la corriente de flujo capilar (flujo transversal) .

De acuerdo con las varias modalidades, el llevar a cabo los pasos iniciales de manera similar a la técnica de flujo laminar, representa un cálculo aproximado de la primera aproximación de la migración de la saturación (ej . , saturación de agua) en la formación modelada, y la primera aproximación es entonces modificada o corregida para tomar en cuenta algunos o todos los efectos físicos arriba notados. Sin embargo, la corrección de dichos efectos físicos no deberá afectar adversamente la longitud del paso de tiempo, como parece ser el caso en la técnica del documento de Osaka et al. arriba citado. En particular, con los resultados de la primera aproximación, las varias modalidades calculan un valor que es el cambio en la saturación dentro de cada bloque de retícula multiplicado por el volumen poroso de la celda dividido entre el tamaño del paso de tiempo. El valor es una indicación del flujo del fluido que ha ocurrido durante un periodo de tiempo. En seguida, y otra vez dentro de cada bloque de retícula, una velocidad total de los fluidos es determinada. Y finalmente, el método resulta resolver ecuaciones simultáneas de Buckley-Leverett modificadas para incluir al menos una, pero en particular abarca una pluralidad, de consideraciones, tales como permeabilidad relativa como entre los hidrocarbonos y el agua en el bloque de retícula, presión capilar, gravedad, o flujo transversal. El cálculo del flujo del fluido, la velocidad del fluido y la resolución de Buckley-Leverett, pueden llevarse a cabo múltiples veces hasta que el valor sea reducido (y en algunos casos minimizado) , pero en algunos casos una sola iteración es suficiente.

Al resolverse, las ecuaciones proporcionan correcciones a la determinación de saturación de agua. A diferencia de Osaka et al., las varias modalidades no resultan en inestabilidad numérica. Dicho de otra manera, las correcciones no imponen limitaciones de paso de tiempo debido a que las correcciones pueden "mover" las saturaciones a través de los límites de bloque de retícula.

Vayamos ahora a un tratamiento más matemático de la corrección de acuerdo con las varias modalidades. Con las saturaciones de agua de la primera aproximación completas, el método calcula, para cada bloque de retícula, un valor residual del flujo del fluido (ej . , agua) usando sustancialmente la siguiente ecuación: en donde Ru(i,j,k) es el valor residual para un bloque de retícula en particular, Sw(i,j,k)n+1 es la saturación calculada para el bloque de retícula en particular, Pv (i ,j ,k) es el volumen poroso dentro 'del bloque de retícula, ñt es el tamaño del paso de tiempo, y Bw es el factor de volumen de formación, y n es el paso de tiempo.

En seguida, la velocidad total de fluidos en las interfaces de cada bloque de retícula, se determinan usando sustancialmente la siguiente ecuación: = ??? , (2) en donde ut es la velocidad total, Ti son los tiempos de transmisibilidad de la movilidad de flujo laminar de la fase i, ?f es el gradiente potencial en la interfaz de cada bloque de retícula, y en donde la fase i es aceite (o) , agua (w) y/o gas (g) .

Finalmente, las ecuaciones de Buckley-Leverett para cada celda se resuelven usando el número de velocidad total calculado a partir de la ecuación (2), con soluciones iterativamente determinadas hasta que el error o valores residuales lleguen a un valor predeterminado, tal como un mínimo. La ecuación de Buckley-Leverett toma sustancialmente la siguiente forma: en donde Sw es saturación (ej . saturación de agua) , ut es la velocidad de fluido total, fw es el flujo de fluido fraccional, t es tiempo, y n es el paso de tiempo. Otra vez, la ecuación (3) sirve como ejemplo pero no limita la técnica a la solución de sólo saturaciones de agua. Otras saturaciones y/o composiciones también podrían resolverse.

De acuerdo con las varias modalidades, los efectos físicos adicionales a ser considerados, se incluyen en las ecuaciones para el flujo de fluido fraccional fu y/o término de velocidad de fluido ut. Por ejemplo, en un sistema bifásico (ej., aceite y agua) en donde la corrección aplicada es para enfocar la permeabilidad relativa entre agua y aceite dentro del bloque de retícula, y considerando una sola dimensión, el f„ podría tomar la siguiente forma: en donde kw es la permeabilidad relativa del agua (dada por la siguiente ecuación) , Kro es la permeabilidad relativa del aceite (dada por la siguiente ecuación) , µ„ es la viscosidad del agua, y µ? es viscosidad del aceite. Sin embargo, las permeabilidades relativas no son constantes en la Ecuación (4). La ecuación relacionada con la permeabilidad relativa del agua Krw Y la permeabilidad relativa del aceite Kro pueden tomar la forma: con los varios parámetros según se¦ definen con anterioridad .

Para no complicar indebidamente la descripción, la expansión particular del flujo de fluido fraccional en las ecuaciones (4), (5) y (6) toma en cuenta la permeabilidad relativa, una sola dimensión y sólo dos fases; sin embargo, un experto en la materia, ahora entendiendo la metodología, podría expandir las consideraciones a múltiples dimensiones y multi-fases/multi-componentes, así como podría tomar en cuenta otros efectos, tales como: gravedad, presión capilar y flujo transversal.

Usando el análisis de Von Neumann, la solución de las anteriores ecuaciones implícitas para saturaciones y/o composiciones, lleva a un método estable incondicional sin limitación alguna en cuanto a paso de tiempo.

La Figura 1 muestra un método de acuerdo con al menos algunas modalidades. En particular, el método inicia (bloque 100) y procede a la formulación de un modelo lógico de una formación subterránea de hidrocarbonos con base en los datos de una formación subterránea real de hidrocarbonos, el modelo comprendiendo una pluralidad de bloques de retícula (bloque 104). En seguida, el método ilustrativo simula la reacción de la formación a extracción de hidrocarbonos durante una pluralidad de pasos de tiempo (bloque 108): al hacer una primera aproximación de la migración a saturación (ej . , saturación de agua) para al menos un bloque de retícula de la pluralidad de bloques de retícula, en donde la migración de la saturación durante al menos un paso de tiempo abarca más de un bloque de retícula (bloque 112) ; y después corrigiendo la primera aproximación de migración de saturación para un efecto no considerado en la primera aproximación (bloque 116) . Finalmente, el método ilustrativo despliega una representación visual de una locación del límite de la saturación de agua (bloque 120) , y el método termina (bloque 124).

Muchas variaciones de las técnicas arriba descritas son posibles. Por ejemplo, en modalidades particulares, los tamaños de los bloques de retícula pueden ser reducidos en áreas activas para reducir los tiempos computacionales . Aún como otro ejemplo, en las áreas en donde hay muy poco o no hay movimiento de fluido como el que se identificó en la primera aproximación, los tamaños de bloque de retícula podrían ser agrandados en dichas áreas, y/o el sistema podría reprimirse de resolver las ecuaciones de Buckley-Leverett en las áreas identificadas como que tienen muy poco o no tienen movimiento de fluido, podría asumirse un sistema bicomponente (ej . , aceite y gas) , y en donde se espera el movimiento de fluido relevante, la complejidad del modelo podría ser incrementada para contar para componentes múltiples (ej . , metano, hexano, butano, etc.)- Más aún, para tipos de fluido particulares y en porosidades particulares de la arenisca, habrá muy poco o nada de movimiento de los fluidos. En las modalidades alternativas, las técnicas antes descritas podrían ser usadas para modelar inicialmente todas las áreas, y después no llevar a cabo ningún modelado en las áreas en donde no es probable ningún movimiento, debido a la viscosidad de fluido y la relación de porosidad (en lugar de una porosidad de recorte arbitrario de la materia relacionada) .

La Figura 2 ilustra a mayor detalle un sistema de cómputo (200), que es ilustrativo de un sistema de cómputo en el que las varias modalidades pueden ser practicadas. El sistema de cómputo (200) comprende un procesador (202), y el procesador se acopla a una memoria principal (204) a manera de un dispositivo de puente (208) . Más aún, el procesador (202) se puede acoplar a un dispositivo de almacenamiento a largo plazo (210) (ej., una unidad de disco duro, disco "floppy", tarjeta de memoria) mediante el dispositivo de puente (208). Los programas ejecutables por el procesador (202) pueden ser almacenados en el dispositivo de almacenamiento (210) y accedidos por el procesador (202) cuando sea necesario. El programa almacenado en el dispositivo de almacenamiento (210) puede comprimir programas para implementar las varias modalidades de la presente descripción, incluyendo programas para implementar la respuesta de la formación de modelaje a las técnicas de extracción. En algunos casos, los programas son copiados desde el dispositivo de almacenamiento (210) hacia la memoria principal (204), y los programas son ejecutados desde la memoria principal (204). Asi, tanto la memoria principal (204) como el dispositivo de almacenamiento (210) se consideran medios de almacenamiento legibles por computadora. Los resultados del modelado mediante el sistema de cómputo (200) pueden ser enviados a un dispositivo de proyección que puede hacer una representación para ser vistos por un ingeniero de depósitos u otra persona experta en la materia.

De la descripción aquí provista, aquellos expertos en la materia son fácilmente capaces de combinar el software creado como se describió con el hardware de cómputo apropiado (incluyendo sistemas de cómputo paralelos) para crear un sistema de cómputo de fines especiales y/o sub-componentes de cómputo con fines especiales, de acuerdo con las varias modalidades, para crear un sistema de cómputo con fines especiales y/o subcomponentes de cómputo para llevar a cabo los métodos de las varias modalidades y/o para crear un medio legible por computadora que almacena un programa de software para implementar los aspectos del método de las varias modalidades .

La anterior discusión pretende ser ilustrativa de los principios y varias modalidades de la presente invención. Numerosas variaciones y modificaciones serán aparentes para los expertos en la materia una vez que la anterior divulgación sea totalmente apreciada. Se pretende que las siguientes reivindicaciones sean interpretadas para abarcar dichas variaciones y modificaciones.

Claims

REIVI DICACIO ES

1. Un método que comprende: la formulación de un modelo lógico de una formación subterránea de hidrocarbonos basada en los datos de una formación de hidrocarbonos real, el modelo comprendiendo una pluralidad de bloques de retícula; la simulación de una reacción de la formación a la extracción de hidrocarbonos durante una pluralidad de pasos de tiempo, la simulación dentro de cada paso de tiempo: haciendo una primera aproximación de la migración de saturaciones para al menos un bloque de retícula de la pluralidad de bloques de retícula, en donde la migración de las saturaciones dentro de al menos un paso de tiempo, abarca más de un bloque de retícula; y después corrigiendo ¦ la primera aproximación de migración de saturaciones para un efecto no considerado en la primera aproximación; y desplegando una representación visual de una locación de saturaciones.

2. El método de la reivindicación 1 en donde la corrección comprende, además, la corrección de al menos uno seleccionado del grupo que consiste de: gravedad; permeabilidad relativa mientras cambia la saturación del agua dentro de un bloque de retícula; presión capilar; y flujo transversal.

3. El método de la reivindicación 1 en donde la corrección comprende, además, la resolución simultánea de una pluralidad de ecuaciones de Buckley-Leverett, las ecuaciones de Buckley-Leverett para al menos algunos de los bloques de retícula, en donde cada ecuación de Buckley-Leverett toma en cuenta un efecto no considerado en la primera aproximación.

4. El método de la reivindicación 3, en donde la resolución simultánea de la pluralidad de ecuaciones de Buckley-Leverett comprende, además, la resolución de las ecuaciones de Buckley-Leverett que corrigen al menos uno seleccionado del grupo de: gravedad; permeabilidad relativa mientras cambia la saturación del agua dentro de un bloque de retícula; presión capilar; y flujo transversal.

5. El método de la reivindicación 1, en donde el modelado comprende, además, el modelado de una reacción para al menos un punto de extracción de hidrocarbonos.

6. El método de la reivindicación 1, en donde la corrección de la primera aproximación comprende, además, hacer la primera aproximación de forma tal, que la migración de las saturaciones durante el periodo de tiempo abarque más de un bloque de retícula.

7. El método de la reivindicación 1, en donde la corrección de la primera aproximación comprende, además, la corrección de la primera aproximación de forma tal, que la migración de las saturaciones durante el periodo de tiempo abarque más de un bloque de retícula.

8. Un medio legible por computadora que almacena un programa que, al ser ejecutado por uno o más procesadores, hace que el (los) procesador (es) : simulen la reacción de un modelo de formación a la extracción de hidrocarbonos durante una pluralidad de pasos de tiempo, la simulación dentro de cada paso de tiempo hace que el procesador: haga una primera aproximación de la migración de saturación para al menos un bloque de retícula de una pluralidad de bloques de retícula, y en al menos una región del modelo de formación, la migración de la saturación abarca más de un bloque de retícula; y corrija la primera aproximación de la migración de saturación, para un efecto que no está considerado en la primera aproximación.

9. El medio legible por computadora de la reivindicación 8, en donde cuando el procesador corrige, el programa, además, hace que el procesador corrija al menos uno seleccionado del grupo que consiste de: gravedad; permeabilidad relativa mientras cambia la saturación del agua dentro de un bloque de retícula; presión capilar; y flujo transversal.

10. El medio legible 'por computadora de la reivindicación 8, en donde cuando el programa corrige, el programa, además, hace que el procesador resuelva simultáneamente una pluralidad de ecuaciones de Buckley-Leverett, las ecuaciones de Buckley-Leverett son cada una para al menos algunos de los bloques de retícula, cada ecuación de Buckley-Leverett resuelta toma en cuenta un efecto no considerado en la primera aproximación.

11. El medio legible por computadora de la reivindicación 10, en donde cuando el procesador resuelve simultáneamente, el programa, además, hace que el procesador resuelva simultáneamente la pluralidad de ecuaciones de Buckley-Levett que corrigen al menos uno seleccionado del grupo que consiste de: gravedad; permeabilidad relativa mientras cambia la saturación del agua dentro de un bloque de retícula; presión capilar; y flujo transversal.

12. El medio legible por computadora de la reivindicación 8, en donde cuando el procesador hace la primera aproximación, el programa hace, además, que el procesador modele la reacción a al menos algún punto de extracción de hidrocarbonos .

13. El medio legible por computadora de la reivindicación 8, en donde cuando el procesador corrige, el programa, además, hace que el procesador corrija la primera aproximación de forma tal que la migración de las saturaciones cambie impuesta por la corrección para al menos una región que abarca más de un bloque de retícula.

14. Un sistema de cómputo que comprende: un procesador; una memoria acoplada al procesador, la memoria almacena un programa que, al ejecutarse por el procesador, hace que el procesador: simule la reacción de un modelo de formación para extracción de hidrocarbonos durante una pluralidad de pasos de tiempo, la simulación dentro de cada paso de tiempo hace que el procesador: haga una primera aproximación de migración de saturaciones para al menos un bloque de retícula de una pluralidad de bloques de retícula, y en al menos una región del modelo de formación la migración de saturación abarca más de un bloque de retícula; y después corrija la primera aproximación de migración de saturaciones para un efecto no considerado en la primera aproximación .

15. El sistema de cómputo de la reivindicación 14, en donde cuando el procesador corrige, el programa, además, hace que el procesador corrija al menos uno seleccionado del grupo que consiste de: gravedad; permeabilidad relativa mientras cambia la saturación del agua dentro de un bloque de retícula; presión capilar; y flujo transversal.

16. El sistema de cómputo de la reivindicación 14, en donde cuando el programa corrige, el programa, además, hace que el procesador resuelva simultáneamente una pluralidad de ecuaciones de Buckley-Leverett , las ecuaciones de Buckley-Leverett son cada una para al menos algunos de los bloques de retícula, y cada ecuación de Buckley-Leverett resuelta toma en cuenta un efecto no considerado en la primera aproximación.

17. El sistema de cómputo de la reivindicación 16, en donde cuando el procesador resuelve simultáneamente, el programa, además, hace que el procesador resuelva simultáneamente la pluralidad de ecuaciones de Buckley-Levett que corrigen al menos uno seleccionado del grupo que consiste de: gravedad; permeabilidad relativa mientras cambia la saturación del agua dentro de un bloque de retícula; presión capilar; y flujo transversal.

18. El sistema de cómputo de la reivindicación 14, en donde cuando el procesador simula, el programa, además, hace que el procesador modele la- reacción a al menos un punto de no inyección de hidrocarbonos y al menos un punto de extracción de hidrocarbonos.

19. El sistema de cómputo de la reivindicación 14, en donde cuando el procesador corrige, el programa, además, hace que el procesador corrija la primera aproximación de forma tal que la migración del cambio de saturación impuesto por la corrección para al menos una región, abarca más de un bloque de retícula.