STREAMLINE

ABSTRACTO
La simulación aerodinámica se convierte en un método preciso y
rápido para la simulación de yacimientos hoy en día. Este método se
basa en desemparejar la ecuación de la presión de la ecuación de la
saturación y transformar la ecuación 3D a las ecuaciones 1D por lo
tanto introduce un factor de la velocidad para arriba en lo referente a la
simulación convencional.
el rastreo racionalizado en el reservorio es un paso principal en la
simulación optimizada. a medida que se traza una línea aerodinámica,
pasa a través de muchos bloques de cuadrícula. el enfoque numérico
de la ecuación de saturación de identificación requiere la captación de
la saturación en cada bloque de cuadrícula cuando se traza una línea
aerodinámica. De esta manera, la saturación versus el tiempo de la
información de vuelo registrado para cada línea aerodinámica trazada
del inyector al productor. Esta información se define en un tiempo
irregular de espacio de vuelo.
los esfuerzos pasados para mover la saturación hacia adelante se
basaban en transformar el tiempo irregular del espacio de vuelo en un
espacio regular. Transformar en una cuadrícula 1D normal da como
resultados promedios de saturaciones, y representa una fuente de
difusión numérica en el método aerodinámico.

En este trabajo presentamos un nuevo método para mover la

saturación hacia adelante basándonos en el tiempo irregular de vuelo
del espacio y comparamos los resultados con el método
común. El resultado de la saturación con este método será más
preciso y el promedio de las saturaciones y la difusión numérica
se reducirá.

INTRODUCCION
Racionalizar la simulación como alternativa a las simulaciones
convencionales se ha hecho más popular debido a varios factores.
Primero, la alta velocidad computacional y menos uso de la memoria
son necesarias. Esto se debe a las ecuaciones de presión de
desacoplamiento de las ecuaciones de saturación, por lo tanto el
número de veces que se debe calcular el campo de presión se reduce
y la simulación de campo de gran escala sería posible (Batycky 1997,
Peddibhotla 1997). En segundo lugar, la solución de ecuaciones
unidimensionales a lo largo de las aerodinámicas alivia el límite de
estabilidad de la cuadrícula subyacente (Batycky et al. 1997, Thiele et
al. 1996, King et al. 1998, Wang et al. 2000, Jang et al. 2002). Tercero,
el acercamiento aerodinámico preserva el frente agudo de la
saturación reduciendo la dispersión numérica,
también los efectos de la orientación de la rejilla reducidos al mínimo
en comparados a los métodos convencionales de la finito-diferencia
(Peddibhotla 1997).
en este trabajo, se estudiaron los diferentes esquemas numéricos que
se pueden aplicar a la ecuación de la saturación en el método
aerodinámico. la aproximación finito-diferente del primer derivado en la
ecuación diferencial es factor dominante para hacer un esquema de la
finito-diferencia más exacto, estable, y computacionalmente eficientes
con menos dispersión numérica. peddibhotla (1997) utilizó un
esquema de disminución total de la variación (TVD) para salve las
ecuaciones de flujo multifásicos para minimizar la dispersión numérica,
y prevenir cualquier ondas no físico.
Acentuamos en las rejillas irregulares para solucionar la ecuación de
la saturación a lo largo de las aerodinámicas. hemos demostrado en
este estudio que la transformación de las rejillas irregulares a las
rejillas regulares y viceversa puede crear efectos de esquemas
numéricos adicionales. también, la exactitud del solucionador
numérico para la ecuación de saturación con menos dispersión
numérica.
Este artículo consta de cuatro secciones. la primera sección establece
las ecuaciones de flujo de dos fases utilizadas en este trabajo. la
siguiente sección proporciona una descripción breve del método y los
procedimientos racionalizados. en la tercera sección se considera el
solucionador de solucionadores numéricos comunes para las rejillas
irregulares. Finalmente, llegamos a los resultados, las discusiones y
las conclusiones
Ecuaciones de Gobierno
La simulación aerodinámica se basa en dos pasos: principalmente, la
ecuación de la presión se resuelve, y entonces la saturación será
obtenida. Estas ecuaciones se derivan de la ecuación de balance de
masa para cada componente definido por el lago (1989) con los
siguientes supuestos:
--fluidos inmiscibles
–flujo incompresible
–la gravedad y las fuerzas capilares se descuidan
La ecuación de la presión elíptica es

y la ecuación de saturación hiperbólica,

Donde

ս representa la velocidad total (de Darcy), la fuente de q o los

términos del fregadero debido a los pozos, φ la porosidad, S la
saturación, K la permeabilidad. Finalmente. para dos fases (aceite y
agua)

Donde la movilidad de cada fase, se da como,

y el flujo fraccional de cada fase,

Método aerodinámico
La ecuación de saturación (EQ. 2) el método de optimización se
desempareja información de un conjunto de ecuaciones de 1D a lo
largo de las líneas aerodinámicas y este método se basa en un
procedimiento secuencial, mientras que en un método de diferencia
finita convencional esto, en lugar de la discreción y la solución
ecuación de saturación en una cuadrícula, las ecuaciones multipleone-
dimensionales se resuelven a lo largo de racionalizaciones
introduciendo el concepto de tiempo de vuelo. el tiempo del vuelo
definido como el requerido a una partícula viaja una distancia a lo
largo de la aerodinámica
en forma diferencial

Los métodos aerodinámicos se basan en una transformación de

coordenadas desde el espacio físico para agilizar las coordenadas de
vuelo para los cálculos de saturación. Esta transformación se basa en
el concepto de BI-streamfunctions y tiempo de vuelo.
el BI-streamfunctions ψ y χ introduced por Bear (1972) para expresar
el campo de la velocidad como,

Está claro que para fluidos incompresibles,

Por lo tanto

Una línea aerodinámica definida por intersección de un valor

constante para ψ con un valor constante de χ. Para aplicaciones
bidimensionales, una línea aerodinámica generada por

Esto se puede utilizar para encontrar el campo de velocidad como,

Como se mencionó, una transformación de coordenadas es requerida

en el método aerodinámico. Esta transformación se ve facilitada por el
hecho de que el jacobiano de la transformación de coordenadas toma
una forma simple extraordinaria utilizando EQ. 8 y EQ. 9

Se obtiene la siguiente relación entre el espacio físico y la hora de las

coordenadas de vuelo,
Una línea aerodinámica definida por intersección de un valor
constante para ψ con un valor constante de (t,ψ,χ). Para aplicaciones
bidimensionales, una línea aerodinámica generada por

Porque u es ortogonal tanto a ψ como a χ,

Ahora, EQ. 2 puede ser transformado en coordenada t como,

La ecuación de la conservación para la fase del agua en el flujo

incompresible de dos fases, lejos de fuentes y fregaderos, estaría
como:

Esta ecuación es pseudo ecuación de balance de masa 1D que se

resuelve a lo largo de las líneas aerodinámicas.
Procedimiento de simulación
Discutimos la idea de desemparejar del flujo flúido multi-dimensional
en ecuaciones unidimensionales múltiples. el procedimiento de la
simulación aerodinámica se resume s un diagrama de flujo en la Fig 1
En primer lugar, las movilidades, λ, se calculan utilizando la
distribución de saturación inicial conocida. Entonces, la ecuación de
presión (EQ. 1) se resuelve para dar la velocidad total, u, y la
distribución de la presión, p.
A continuación, el dominio de flujo se descompone en un número de
optimizaciones mediante el seguimiento de trayectorias aerodinámicas
desde el campo de velocidad calculado. El algoritmo de rastreo
utilizado aquí, es un método desarrollado por Pollock (1988).
este método implica el uso de una aproximación lineal de la velocidad
de la pieza con respecto a cada rejilla-bloque seguida por la
integración exacta para la posición. en este momento, todas las
posiciones a lo largo de una línea aerodinámica se pueden expresar
con TOF en vez de coordenadas cartesianos.
Ahora, la velocidad total se mantiene fija mientras que la saturación se
moverá hacia adelante por un paso de tiempo dado. La saturación se
puede obtener tanto analítica como numéricamente. La solución
analítica a la saturación es válida sólo en condiciones iniciales
uniformes. Finalmente, las soluciones se reasignan al dominio de flujo
original, y la distribución final de saturación se adquiere en función de
la ubicación y el tiempo.
Pseudo solución de saturación 1D
Una vez que las aerolíneas se trazan en el dominio, el siguiente paso
es resolver la ecuación multifase (EQ. 20) a lo largo de los
aerodinámicos. Esta ecuación se puede resolver analítica y
numéricamente. La solución analítica (ecuación de Buckley-Leverett) a
esta ecuación para el flujo incompresible de dos fases es dada cerca,
Por lo tanto, sólo es necesario definir el tiempo de vuelo a lo largo de
una línea aerodinámica para determinar la saturación de agua y los
flujos fraccionarios. La desventaja de este método es que se requiere
una condición inicial uniforme a lo largo de los aerodinámicos.

Aplicar la ley de Darcy para

Calcule el campo de la presión encontrar el campo de velocidad
usando EQ. 1 así que el campo de
la saturación está en 𝒕𝒏
Traza aerodinámica de
inyectores a productores

Recoger la información de
saturación actual de cada
gridblock que pasa la línea
aerodinámica
Promedio de todas las propiedades
aerodinámicas dentro de cada gridblock
para determinar la distribución de
saturación en promedio todas las Pasar la saturación en un
propiedades aerodinámicas dentro de solucionador numérico 1D y
cada gridblock para determinar la mover las saturaciones hacia
adelante por el tiempo.
distribución de saturación en 𝒕𝒏+𝟏

Fig. 1-procedimiento de simulación optimizada Vuelva a asignar el nuevo

perfil de saturación a la
cuadrícula de fondo.
Para superar esta hipótesis, el EQ. 20 debe ser obtenido
numéricamente.
solución numérica 1D
Los métodos numéricos comunes para el cálculo de la saturación en la
simulación aerodinámica se pueden dividir a cuatro pasos.

El primer paso es definir la distribución de saturación presente a lo

largo de una línea aerodinámica que se traza desde el inyector hasta
el productor. La saturación versus la información de tiempo de vuelo
registrada para cada línea aerodinámica se define en una cuadrícula
de t irregular. Esto se debe a las altas velocidades de flujo cerca de
los inyectores y las velocidades de flujo bajo en el centro del dominio.
Fig. 2 muestra la distribución de agua para el patrón de punto de
trimestre de cinco heterogéneos en un momento específico, la
saturación versus el tiempo de vuelo para una línea aerodinámica se
muestra en la Fig. 3.

En el segundo paso, la saturación versus el tiempo de vuelo se

transformó en una cuadrícula t regularmente espaciada. Esto
simplifica el cálculo de los flujos de entrenudo dentro del solucionador
numérico 1D, mientras que esta simplificación da como resultado un
promedio de saturaciones e introduce una fuente de difusión
numérica. El proceso de transformación que se basa en la ley de
conservación de masas es el siguiente (Batycky 1997),

El tercer paso es mover la saturación hacia adelante por Δt. La forma

discreta de EQ. 20 es la siguiente (Peddibhotla 1997),

donde 𝑭𝒊+𝟏/𝟐 el flujo en el límite entre los nodos i y i+1. La Fig. 4

muestra los nodos y flujos en una línea aerodinámica. Es evidente que
cualquier solución débil de EQ. 20 tiene una variación no creciente (le
Veque 1990). Esto es suficiente para obtener algunos resultados de
convergencia y además no garantiza ninguna generación de oscilación
espuria. Los métodos con esta propiedad se denominan métodos de
disminución total de la variación (TVD). Harten (1983) introdujo la
noción del esquema TVD y señaló un simple criterio algebraico para el
esquema. Cuando se busca una solución numérica, la variación de la
solución discreta debe disminuir. La variación total, TV ( 𝑭 𝒏+𝟏 ) de la
solución se define por:
y tan para el esquema de TVD, la condición siguiente debe ser
satisfecha:

Tenemos un esquema explícito general en forma conservadora

usando una función ϕ(r) de limitación TVD), (Sweby 1990)

Cuando las funciones de limitación se expresan en el cuadro 1. La

magnitud de la limitación depende de la suavidad de los datos y se
mide por el cociente del gradiente consecutivo r,

Sweby (1990) deriva las condiciones algebraicas en la función de

limitador que garantizan la propiedad TVD que es,

La ecuación anterior conduce a los limitadores TVD para las

formulaciones matemáticas de orden superior dadas en la tabla 2.
Estos esquemas de alta resolución (esquemas TVD) son
conservadores, generalmente (en la mayoría de los lugares)
precisión de segundo orden, menos dispersión numérica, no oscilante
en la naturaleza y capaz de resolver la discontinuidad en la solución.
El método de un solo punto por aguas arriba resulta en una solución
estable pero dispersa.
La saturación es determinada por,
Fig. 2-distribución de la saturación de agua a lo largo de un depósito
heterogéneo

Fig. 3-saturación de agua versus tiempo de vuelo para

aerodinamización trazada en la Fig. 2

Fig. 4-nodos y fundentes a lo largo de una línea aerodinámica

Tabla 1-funciones de limitación
Un solo punto aguas arriba Φ(r) =0
Ponderación de punto medio Φ(r) =1
Dos puntos aguas arriba Φ(r) =r
Leonard Scheme Φ(r) =r/3+2/3

Para el tiempo explícito, el tamaño máximo estable del paso del

tiempo (Δt) permitido para resolver EQ. 29 es determinado por el
criterio de la estabilidad de Courant-Fredrich-Lewy (Batycky 1997).
Peddibhotla (1997) usó el tercer plan de Leonard para su estudio y
mostró un buen acuerdo entre la solución TVD y la solución analítica,
mientras que el método de punto único de subida tiende a difundir los
frentes de saturación.
Tabla 2-TVD limitadores para las formulaciones matemáticas de
mayor orden
Un solo punto aguas arriba ------
Ponderación de punto medio Φ(r) =MAX(0,MIN(r,2))
Dos puntos aguas arriba Φ(r) =MAX(0,MIN(r,2))
Leonard Scheme Φ(r) = MAX(0,MIN(2.2r, 𝜱𝐿𝐸 ))

En el cuarto paso, después de que una solución aerodinámica se ha

movido adelante por Δt en los nodos t spaced regularmente, se
transfiere de nuevo a la rejilla irregular original. Una vez más, esta
transferencia es una fuente de difusión numérica Esta transformación
tiene un mayor efecto de promediación en el cambio del perfil de
saturación que cualquier difusión numérica asociada con el
solucionador 1D.

Una vez que las variables de saturación se transforman a la cuadrícula

t irregular, se pueden mapear a una nueva cuadrícula de saturación.
El procedimiento para calcular la saturación media de bloques de
cuadrícula es discutido por Batycky (1997).

1D solución numérica basada en los TOF irregulares

Como se ha comentado anteriormente, cuando se traza una línea
aerodinámica, la saturación versus el tiempo de grabación de la
información de vuelo se define en una cuadrícula de t irregular. La
transferencia de saturación versus tiempo de vuelo a los nodos de
espacio t regulares y viceversa con el fin de mover la saturación hacia
adelante conduce a la dispersión numérica. Además de la dispersión
numérica debido a la transformación (regular/irregular) hay efecto
de la difusión dentro del solucionador numérico 1D. Para reducir al
mínimo la última dispersión numérica, uno puede utilizar un método
más alto de la orden como el trabajo de Peddibehot (1997).

El uso de nodos espaciados irregularmente t es un método para lidiar

con este problema. Esto requiere un esquema de diferencia finita que
sea adecuado para nodos regulares e irregulares. Los métodos de
orden inferior pueden crear efectos de frotis numéricos adicionales de
a medida que la irregularidad del nodo se vuelve importante. Métodos
más altos de la orden, tales como tercer esquema de Leonard de la
orden (Leonard 1979) sea bueno en sus exactitudes del alto-orden y
en la reducción de efectos que manchan numéricos. Estos métodos,
sin embargo, están plagados por su generación de oscilaciones
espurias o retrasos en el vecindario de discontinuidades causadas por
limitaciones de monotonicidad adicionales (Patel et al. 1985).
El orden de precisión de los diferentes esquemas se indica en el
cuadro 3. Con el uso de nodos espaciados t irregulares, el orden de la
exactitud para el esquema de Leonard sería el primer y este esquema
produce los efectos que manchan numéricos causados por los nodos
irregular.

Tabla 3-orden de la exactitud para diversos métodos

EL PLAN Orden de la exactitud

Solo punto upstream Nodos regulares Primer orden
Nodos irregulares Primer orden

Ponderación de punto medio Nodos regulares Segundo orden

Nodos irregulares Primer orden

Dos puntos arriba Nodos regulares Segundo orden

Nodos irregulares Primer orden

Leonard Scheme Nodos regulares Tercera orden

Nodos irregulares Primer orden

Ahora, para mantener el orden de la exactitud para el solucionador

numérico, una aproximación de la tercera orden del primer derivado
espacial se puede obtener usando la extensión de la serie de Taylor
(Liu et al. 1995) de cuatro nodos. (la Fig. 4 muestra los nodos a lo
largo de una línea aerodinámica) Esta aproximación para el nodo τ yo
sería:

Los coeficientes, son 𝜶𝒊 , 𝜷𝒊 𝒚 𝜸𝒊 son:

donde Δt, Δ𝒕𝒊−𝟏 and Δ𝒕𝒊−𝟐 se muestran en la figura 4. Para los nodos
espaciados t regulares, estos coeficientes son los siguientes,

Se puede ver que la expansión cambió a Leonard Scheme.

Restricción de disminución de variación total
La aplicación de EQ. 30 para aproximar la primera derivada espacial
en EQ. 20 conduce a la siguiente forma explícita:

Para minimizar los artefactos tales como la dispersión numérica y

prevenir las oscilaciones no físicas, la variación total que disminuye
debe ser satisfecha; Así pues, el esquema explícito general en forma
conservadora sería (Liu et al. 1995):

donde está la función de limitación,

La medida de la suavidad de los datos, r, se determina como sigue,

La condición suficiente para que EQ. 33 sea la variación total

disminuyendo será como

donde es el paso de tiempo para la solución numérica de una

dimensión a lo largo de la línea aerodinámica
Resultados y discusiones
En primer lugar, se muestran los resultados del método de cálculo de
saturación común (basado en nodos t regulares). Este método según
lo discutido, transforma los espacios t irregulares de los espacios t
regulares. La Fig. 5 muestra una transformación de nodos irregulares
a uno normal para una línea aerodinámica.
Las diferencias en los perfiles de saturación entre los nodos regulares
e irregulares muestran el efecto del método de transformación donde
los valores locales de saturación de
alta y baja se promedian en ciertas ubicaciones a lo largo de la
coordenada t.

Fig. 5-saturación de agua versus tiempo de vuelo, efecto de la

transformación
En la Fig. 6 se muestra el efecto del número de nodos en la exactitud
de la transformación. A medida que aumenta el número de nodos, el
resultado de los nodos regulares casi coincide con el irregular. El
aumento del número de nodos, sin embargo, da como resultado un
aumento significativo del tiempo de CPU. Y los artefactos tales como
la dispersión numérica y la oscilación no física dentro de los
solucionadores numéricos 1D serán aumentados.
La Fig. 7 muestra el impacto de un gran número de nodos regulares
dentro del solucionador numérico 1D.
Fig. 6-efecto del número de nodos en la precisión de la transformación

Fig. 7-impacto de un gran número de nodos dentro del solucionador

numérico 1D
Como se mencionó anteriormente, la transformación a nodos
regulares simplifica el cálculo de flujos de entrenudo dentro del
solucionador numérico 1D, pero resulta en algunos problemas
numéricos.
para hacer frente a estos problemas, se puede utilizar los nodos
irregulares. Los solucionadores numéricos 1D dados en la tabla 1 se
convierten en precisión de bajo orden con nodos irregulares en
comparación con los nodos regulares. Esto puede ser
confirmado por la Fig. 8 y 9.
.

Fig. 8-precisión de dos métodos con nodos regulares

La Fig. 8 ilustra el porcentaje de error de cada nodo a lo largo de una
línea aerodinámica con espacios τ regulares. Como era de esperar, el
esquema de Leonard da menos error en este caso que un solo punto
río arriba método
Fig. 9-precisión de dos métodos con nodos irregulares
En la Fig. 9, se puede ver que Leonard Scheme ya no tiene una
precisión de tercer orden cuando se aplica a los nodos irregulares y
oscilaciones falsas aparecieron en los resultados. Este resultado
también fue confirmado por Patel et al. (1985). En este caso, el
método upstream de un solo punto sólo tiene un error mayor en
comparación con los nodos regulares.
Ahora, comparamos los resultados del solucionador numérico dado 1D
con el esquema de Leonard y el método ascendente del solo punto. La
Fig. 10 muestra el porcentaje de error entre el solucionador numérico
1D dado y la solución analítica para los nodos regulares a lo largo de
una línea aerodinámica. El rendimiento de este método es tan bueno
como el de Leonard Scheme con nodos regulares.

Fig. 10-exactitud del método dado y del esquema de Leonard con los
nodos regulares
La Fig. 11 ilustra el porcentaje de error de los diferentes métodos con
nodos irregulares. El solucionador numérico dado 1D mantiene su
exactitud con los nodos irregulares

Fig. 11-precisión de diferentes métodos con nodos irregulares

El esquema de Leonard y el método ascendente del solo punto
demuestran el error grande para los nodos irregulares, mientras que el
método dado no, consecuentemente, para reducir la dispersión
numérica, el método dado es mejor que método ascendente del solo
punto y esquema de Leonard.

Waterflood sintético
En esta sección presentamos resultados numéricos para los métodos
(basados en nodos regulares o irregulares) discutidos en secciones
anteriores. Los parámetros utilizados para evaluar la solución
numérica se dan en la tabla
4. En la Fig. 12 se muestra el campo de la permeabilidad sintética.
Fig. 13 muestra la distribución de la saturación de agua después de
0,4 PVI. La trama muestra los contornos de saturación para el
desplazamiento inmiscible. Consideramos que este resultado
representa una aproximación cercana a la solución analítica
mediante nodos irregulares. Es aparente la reducción del frotis de
dispersión numérica frontal y baja en este método en comparación con
los obtenidos con el método basado en nodos regulares.
Tabla 4-parámetros QFS para waterflood

Parámetros Valores
Dimensión en la I-dirección, ft 600
Dimensión en la J-dirección, ft 600
Espesor ft 10
Rejillas del reservorio 30x30x1
Presión del agujero inferior del 4000
inyector, PSI
Presión inferior del agujero del 1000
productor, PSI
Porosidad, fracción 0.2
Viscosidad del aceite, CP 3
Viscosidad del agua, CP 1
Saturación inicial de aceite, fracción 1
Saturación de aceite residual, fracción 0
Número de optimizaciones 200

Fig. 12-distribución de permeabilidad para QFS heterogéneos de

waterflood
El impacto de la dispersión numérica en la saturación del agua
especialmente en los dedos del agua es claro.
Fig. 13-distribución de saturación de agua para diferentes métodos
La Fig. 14 es una parcela de tarifa de aceite versus tiempo para tres
soluciones diferentes a lo largo de las aerodinámicas. El tiempo
temprano de la brecha del agua para el método común es debido a la
dispersión numérica. En comparación con la solución analítica, el
método dado conduce a una menor dispersión numérica y también la
distribución de la saturación de agua en este caso es más parecida a
la solución analítica.
Fig. 14-tasa de aceite versus tiempo para tres soluciones diferentes a
lo largo de aerodinamiza
Por los resultados numéricos considerados en esta sección, el método
se realiza bien y reduce drásticamente el frotis frontal.
Conclusión
El frotamiento frontal y la predicción inexacta del tiempo de avance es
uno de los principales problemas en la simulación de yacimientos.
Este problema es
resultado por la dispersión numérica asociada a los esquemas
numéricos.
Aquí, comparamos los algunos esquemas numéricos y enfatizamos en
el método para resolver la ecuación de saturación a lo largo de la línea
aerodinámica. Los resultados indicaron que los esquemas numéricos
comunes cuando se aplican a los nodos irregulares conducen a
errores más grandes con la dispersión numérica adicional. Se
comparó cada esquema numérico con la solución exacta (analítica)
para obtener el error asociado a cada esquema.
Los métodos de alto orden se convierten en bajo orden con nodos
irregulares mientras que el método dado mantiene su exactitud con las
rejillas irregulares.
también, la sensibilidad de este método a la irregularidad es baja y se
puede aplicar a cualquier nodo espaciado t irregular.
El Reconocimiento
A los autores les gustaría apreciar la dirección y el personal de la
Universidad de tecnología petrolera, el Departamento de ingeniería
química y petrolera de la Universidad de Calgary, y el Instituto de
investigación de la industria petrolera de Irán.