FACULTAD DE CIENCIA Y TECNOLOGÍA

CARRERA DE INGENIERIA DE PETROLEO Y GAS NATURAL

ESTUDIO DE LA ESTIMACION DE HIDROCARBUROS APLICANDO EL

METODO DINAMICO
ELL
MATERIA: SIMULACION DE RESERVORIOS (PGP-310)
GRUPO: # 1
UNIVERSITARIOS:
- Cali Fernández Giovanna Natividad
- Fernández Garcia Orlando
- Icallu Nuñez Jhonny Javier
- Palacios Nuñez Adalberto Jhonatan
- Palacios Rodriguez Juan Fernando
- Paniagua Herrera Alan
- Peñaranda Castro Melisa Sidney
- Rodriguez Apaza José Manuel
- Rodriguez Escudero Estrella
- Tito Patiño Alfredo
DOCENTE: ING. JHON ALEX LEON SENO
Sucre - Bolivia

i
 INDICE
INTRODUCCIÓN ...........................................................................................................................1
ANTECEDENTES ...........................................................................................................................2
JUSTIFICACIÓN ............................................................................................................................2
OBJETO DE ESTUDIO ..................................................................................................................2
OBJETIVOS .....................................................................................................................................3
OBJETIVO GENERAL ..............................................................................................................3
OBJETIVOS ESPECÍFICOS......................................................................................................3
MARCO TEORICO Y CONCEPTUAL ........................................................................................3
1. Ecuación de Balance de Materia .............................................................................................3
1.1. Consideraciones y aplicaciones .........................................................................................4
1.2. Limitaciones de la ecuación de balance de materia ........................................................4
1.3. Aplicación de la ecuación de balance de materia ............................................................6
2. Derivación de la Ecuación de Balance de Materia .................................................................7
2.1. Aplicación de la ecuación de balance de materia: ...........................................................8
3. Ecuación de balance de materia para yacimientos de gas .....................................................8
4. Ecuación general de balance de materia para yacimientos de petróleo .............................11
5. Otros usos de la ecuación de balance de materia .................................................................14
6. Lineación de la ecuación de Balance de Materia mediante el método de Havlena y Odeh
.....................................................................................................................................................15
6.1. Predicción del comportamiento y recobro final de un yacimiento...............................15
6.1.1. Intrusión de agua ..........................................................................................................15
6.1.2. Modelo de estado estable (Schilthuis) .........................................................................15
6.1.3. Modelo de estado inestable (Hurst and Van Everdingen) .........................................18
7. Método de Fetkovich para Acuíferos Finitos .......................................................................22
CONCLUSIONES ..........................................................................................................................26
RECOMENDACIONES ................................................................................................................26
BIBLIOGRAFIA ............................................................................................................................27
ANEXOS .........................................................................................................................................28
EJERCICIO #1 (APLICANDO PARA UN RESERVORIO DE GAS)......................................28
EJERCICIO #2 (APLICANDO PARA UN RESERVORIO DE PETROLEO) ........................31

ii
INTRODUCCIÓN
El estudio dinámico se encarga de caracterizar apropiadamente un yacimiento, puede mejorar
drásticamente su desarrollo y su productividad. La metodología de caracterización se desarrolla
en dos etapas; una etapa de caracterización estática y otra de caracterización dinámica. En la
primera se define las características físicas del volumen de roca a condiciones estáticas, mientras
que en la segunda se describe la interacción de los fluidos dentro del volumen de roca a
condiciones Dinámicas. Modelo estático; en esta etapa se define la geometría del yacimiento y
se describen los parámetros petrofísicos, modelo Dinámico; esta etapa analiza la interacción
dinámica roca- fluido del yacimiento; el propósito fundamental es desarrollar metodologías que
permitan comprender de una manera integral como se desplazan los fluidos en el sistema poroso
de la roca. La elaboración del modelo dinámico se inicia durante el desarrollo de la fase II y la
fase III de un estudio integrado. La fase II corresponde a la caracterización física y energética
del yacimiento, fundamentándose en la aplicación de ingeniería convencional de yacimiento, la
cual permitirá realizar el análisis de los datos históricos de producción y presión del yacimiento,
análisis de fluidos y de núcleos disponibles, para la generación de los siguientes modelos:
presiones, termodinámico e interacción roca-fluido, entre otros. así como también, involucra un
estudio de factibilidad, en cual permitirá al equipo multidisciplinario, luego de un análisis
técnico y financiero, iniciar la realización del modelaje numérico del yacimiento. La fase III
involucra la etapa de simulación numérica del yacimiento, la cual se realiza mediante la
construcción de un modelo numérico a partir de la integración de los modelos estáticos e
ingeniería convencional de yacimientos, realizados durante la fase de caracterización fase II
para un ciclo o estudio específico, la generación del modelo dinámico, se inicia con la
recolección y análisis de los datos referentes al medio poroso del yacimiento continua con la
construcción del modelo de simulación numérica cotejado y culmina con la evaluación de
escenarios que involucran el esquema de recuperación base y esquemas de explotación que
determinan la potencialidad de producción del yacimiento. Posteriormente los resultados de las
evaluaciones son entregados a la gerencia responsable de elaborar el modelo.

1
ANTECEDENTES
La ecuación de balance de materia (EBM), resulta de un balance volumétrico que iguala la
producción acumulada de fluidos, expresada como una extracción, y la expansión de los fluidos
como resultado de una caída de presión en el yacimiento. Se basa en la ley de conservación de
la materia, que establece que la masa de un sistema cerrado permanece siempre constante (a
menos que se produzcan cambios en el sistema). La forma general de la ecuación de balance de
materia (EGBM) fue desarrollada inicialmente por Schilthuis en 1936. La EBM establece que
la diferencia entre la cantidad de fluidos iniciales en el yacimiento y la cantidad de fluidos
remanentes en el yacimiento es igual a la cantidad de fluidos producidos. Una vez obtenida la
información necesaria de este proyecto estaremos capacitados para poder resolver ejercicios que
conciernen al tema tanto de la forma general como también linealmente.

JUSTIFICACIÓN
Este trabajo investigativo presenta el desarrollo del modelo dinámico en un yacimiento, donde
se aplicó nuestros conocimientos adquiridos a lo largo de la carrera, así mismo se especifican
las ecuaciones empleadas en la caracterización de los fluidos presentes, (Ecuación de Balance
de Materia EBM) para evaluar el volumen original de hidrocarburos y conocer los
mecanismos de empuje asociados a los yacimientos para replantear estrategias de producción,
buscando siempre maximizar el factor de recuperación de hidrocarburos.

OBJETO DE ESTUDIO
Los modelos de un yacimiento como tal, se basan en la recolección de información y la
revisión geológica los cuales son datos importantes y estratégicos para conocer los volúmenes
originales de hidrocarburos. donde la ecuación de balance de materia es una herramienta
fundamental para el ajuste del comportamiento de un yacimiento

2
OBJETIVOS
OBJETIVO GENERAL
Determinar el volumen de reservorio aplicando el método de balance de materia (E.B.M.)

OBJETIVOS ESPECÍFICOS
 Describir los factores involucrados en la Ecuación de Balance de Materia para yacimientos
 Demostrar la derivación del Balance de Materia aplicando la ecuación de Schilthuis
 Realizar un ejemplo aplicando la ecuación de balance de materia.

MARCO TEORICO Y CONCEPTUAL

El estudio dinámico se encarga de caracterizar apropiadamente un yacimiento, puede mejorar
dramáticamente su desarrollo y su productividad. La metodología de caracterización se
desarrolla en dos etapas; una etapa de caracterización estática y otra dinámica. En la primera se
define las características del volumen de roca a condiciones estáticas, mientras que en la segunda
se describe la integración de los fluidos dentro del volumen de roca a condiciones dinámicas.

1. Ecuación de Balance de Materia

La EBM se reconoce desde hace mucho tiempo como una de las herramientas básicas para
interpretar y predecir el comportamiento de los yacimientos, utilizando para ello un modelo tipo
tanque que los describe basándose en la relación que debe existir siempre en un yacimiento de
petróleo y gas que produce a condiciones de equilibrio. Es decir, las variaciones de presión y
saturación se consideran semejantes en tiempo, pero no en posición. Aunque en algunos casos
es posible resolver simultáneamente para encontrar el volumen inicial de hidrocarburos y el
influjo de agua, generalmente, alguno de estos datos debe ser conocido a través de la data o
métodos que no dependen de los cálculos del balance de materiales. Uno de los usos más
importantes de las ecuaciones es predecir el efecto de la tasa de producción y/o la tasa de
inyección (de petróleo o agua) en la presión del yacimiento; por lo tanto, es deseable conocer
en adelanto el volumen de petróleo inicial y la razón m mediante datos precisos de núcleos y
perfiles. La presencia de un acuífero se indica usualmente mediante evidencia geológica, sin
embargo, el balance de materiales puede ser usado para detectar la existencia de acuíferos
mediante el cálculo del volumen inicial de hidrocarburos en períodos de producción sucesivos,

3
asumiendo cero influjo de agua. A menos que otros factores de complicación estén presentes, la
constancia en el valor calculado de N y/o G indica la presencia de un yacimiento volumétrico,
y los valores continuamente cambiantes de N y G indican la presencia de un acuífero.

1.1. Consideraciones y aplicaciones

 Volumen poroso constante. No existe compactación ni subsidencia.
 El PVT es representativo del yacimiento
 Proceso isotérmico
 cw y cf son despreciables
 Se considera equilibrio termodinámico entre el gas y el petróleo a presión y temperatura
de yacimiento
 Dimensión cero

Las consideraciones mencionadas según bibliografía utilizada coinciden en la mayoría de

autores reconocidos en el campo de Hidrocarburos, pero para profundizar a continuación se da
una breve explicación de ciertas limitaciones para el uso de la EBM, gracias a las mismas se
puede realizar una aplicación correcta con resultados exactos y precisos, tales resultados llevan
a uno de nuestros objetivos de toda Ingeniería, en especial para la honorable carrera de Ingeniera
en Petroleó y Gas Natural. Ya que una variación de cifras puede causar terribles e irreversibles
catástrofes.

1.2. Limitaciones de la ecuación de balance de materia

Algunas de las limitaciones y suposiciones más importantes con respecto al desarrollo y
aplicación de la EBM son las siguientes: 1. El EBM considera que el yacimiento tiene un
volumen poroso constante ocupado por diferentes fluidos. Así, en un yacimiento de gas seco, el
espacio poroso está ocupado por diferentes fluidos, por gas y agua connata; en uno de petróleo
no saturado o en el punto de saturación, por petróleo y agua connata; y, si existe una capa de
gas, por petróleo, agua connata y gas libre. Con posterioridad al comienzo de la producción del
yacimiento, parte del volumen poroso, considerado constante, puede estar ocupado por agua de
invasión en el caso de que exista un acuífero activo. 2. Supone constante la temperatura del
yacimiento, esto es, que en el proceso de producción no ocurre un cambio considerable de

4
temperatura cuando se extrae petróleo y gas. 3. Considera que existe condiciones de equilibrio
en el yacimiento en cualquier tiempo. Esto significa que los efectos de caída de presión alrededor
del pozo no se toman en cuenta y que la saturación de líquido es uniforme a través de la zona
de petróleo. En otras palabras, a un tiempo en particular, la razón de permeabilidades relativas
gas-petróleo (krg/kro) es constante en toda la zona de petróleo, lo cual incluye la suposición de
que no existe segregación por gravedad. En el caso de yacimientos que tienen una capa de gas
inicial, se considera que no ocurre bonificación en los pozos y que, por tanto, los volúmenes de
la capa de gas y de la zona de petróleo no cambian en el tiempo. Se supone, también, que
cualquier gas que sale de la capa, debido a su expansión, se distribuye uniformemente a través
de la zona de petróleo. 4. Las propiedades PVT disponibles o estimadas deben ser
representativas del yacimiento y relacionar la producción con los datos de vaciamiento. Esto
quiere decir que el mecanismo de liberación del gas que se aplica a las muestras de fluidos para
determinar los datos PVT, debe ser similar al de yacimiento. Usualmente, el mecanismo más
representativo es la vaporización diferencial pues toma en cuenta los cambios que ocurren en la
composición de los fluidos. Ahora bien, en el caso de yacimientos de petróleo volátiles, que
dependen fundamentalmente de la presión, tales cambios se consideran despreciables. 5. La
recuperación es independiente de la tasa. 6. Supone que la producción es totalmente una
consecuencia de liberación del gas en solución y de la expansión del gas liberado del petróleo y
de una capa de gas inicial, cuando disminuye la presión del yacimiento. Esto incluye la hipótesis
de que no hay inyección de fluidos, que el agua es inmóvil, que no hay producción ni entrada
de agua, y que el agua del yacimiento y la compresibilidad de la roca son muy pequeñas y se
pueden despreciar. Nótese que esta suposición solo es válida por debajo del punto de burbujeo
y que estos efectos deben ser considerados por encima del punto de burbujeo. 7. No considera
el factor geométrico del yacimiento, ya que resulta casi imposible determinar la distribución de
los fluidos en la estructura o en los pozos. 8. Requiere cierto grado de explotación del yacimiento
a fin de disponer de suficientes datos de producción y de presión. 9. Generalmente, tanto para
yacimientos de gas seco como para yacimientos de petróleo, el factor volumétrico del agua en
la formación y la solubilidad (razón gas disuelto-agua) se consideran iguales a la unidad y a
cero, respectivamente La exactitud de los valores calculados depende de la precisión de los datos
disponibles para sustituir en la ecuación y en las diversas hipótesis subyacentes en las

5
ecuaciones. Una de dichas hipótesis es el logro del equilibrio termodinámico en el yacimiento,
principalmente entre el petróleo y su gas en solución. También se asume implícitamente que la
data de PVT usada en el balance de materiales se obtiene mediante procesos de liberación de
gas que casi duplican los procesos de liberación de gas en el yacimiento, en el pozo y en los
separadores en superficie, por ende, se producen errores considerables en los resultados del
balance de materiales.

1.3. Aplicación de la ecuación de balance de materia

 Determinar hidrocarburos iniciales en el yacimiento
 Estimar el petróleo (POES) y el gas (GOES) en sitio.
 Estimar el tamaño de la capa de gas.
 Estimar la presencia, tipo y tamaño de acuífero.
 Predecir el comportamiento de presión conociendo la historia de producción del
yacimiento. Estimar las profundidades de los contactos agua-petróleo, gas-petróleo y
agua-gas.
 Predecir el comportamiento futuro de los yacimientos.
 Evaluar We conociendo N o G.
 Predecir el comportamiento y recobro final de yacimientos especialmente los que producen
por gas en solución o depleción.
 Evaluar factores de recobro.

En forma práctica, es normal que en los yacimientos de aceite el efecto de la disminución del
volumen poroso disponible a hidrocarburos por compactación se desprecie cuando se está por
debajo del punto de burbuja. Por su simplicidad, es normal cometer abusos al momento de usar
los modelos de balance de materia. Algunos abusos típicamente cometidos al utilizar estos
métodos se indican a continuación:

 Usar los modelos de balance de materia con poca información.

 Usar información inconsistente con las consideraciones.
 No complementar los estudios con otras herramientas pertinentes, como estudios
volumétricos y de simulación de yacimientos.

6
2. Derivación de la Ecuación de Balance de Materia
La literatura sobre ingeniería de yacimientos contiene numerosos estudios sobre la forma
general de la EBM, pero en general todos están basados en la presentación original desarrollada
por Schilthuis. En su forma más simple la EBM está estructurada para mantener un balance de
todos los materiales que entran, salen y se acumulan en volumen de roca que contiene petróleo,
agua y gas.

Debido a que este se considera constante, la suma algebraica de los cambios de volumen
(incluida la producción e inyección) del petróleo, el gas libre y el agua debe ser igual a cero. En
otras palabras, la ecuación se deriva considerando un balance volumétrico del yacimiento que
iguala lo que sale (producción acumulada) con la expansión que experimentan los fluidos como
consecuencia de una caída de presión finita.

Las siguientes suposiciones son necesarias para la obtención de la EGBM:

El yacimiento puede tener una fase inicial de vapor y una fase inicial de líquido.
El gas se puede disolver en la fase líquida.
El aceite puede ser volátil en la fase de vapor.

7
 El agua puede invadir al yacimiento desde el acuífero durante la producción.
El agua y la roca son compresibles. Si bien cualquier conjunto consistente de unidades es
satisfactorio, los términos utilizados en estas ecuaciones se expresan en unidades de campo
Cambio neto en el volumen de producción de petróleo y agua= al cambio de volumen
de la capa de gas + Gas en solución liberado + Cambio de volumen en el petróleo +
Expansión del agua connata + Expansión de la roca

2.1. Aplicación de la ecuación de balance de materia:

1. Determinar hidrocarburos iniciales en el yacimiento
2. Evaluar We conociendo N o G
3. Predecir el comportamiento y recobro final de yacimientos especialmente los que
producen por gas en solución o depleción
4. Evaluar factores de recobro

3. Ecuación de balance de materia para yacimientos de gas

El balance de materia expresa que la masa de hidrocarburos producidos es igual a los moles
iniciales menos los remanentes, matemáticamente;

𝑥𝑝 = 𝑥𝑖 − 𝑥𝑟

El volumen inicial de hidrocarburos es igual al volumen remanente (la cantidad de masa si

varia), de modo que;

𝑉𝑖 = 𝑉𝑟

El volumen remanente considera el volumen inicial menos el agua que intruye más el agua que
se produce, es decir:

𝑉𝑟 = 𝑉𝑖 − (𝑊𝑒 − 𝛽𝑤 𝑊𝑃 )

Luego, la ecuación de balance de materia, para yacimientos de gas, con el factor volumétrico
dado en 𝑝𝑖𝑒 3 /𝑝𝑐𝑛.

𝑉𝑖 = 𝑉𝑟; 𝐺𝛽𝑔𝑖 (𝐺 − 𝐺𝑃 )𝛽𝑔

8
Simplificando;

𝐺𝛽𝑔𝑖 = 𝐺𝛽𝑔 − 𝐺𝑃 𝛽𝑔

𝐺(𝛽𝑔𝑖 − 𝛽𝑔 ) = −𝐺𝑃 𝛽𝑔

𝛽𝑔 − 𝛽𝑔𝑖 𝐺𝑃
=
𝛽𝑔 𝐺

Si el factor volumétrico se expresa en 𝑝𝑐𝑛/𝑝𝑖𝑒 3

𝐺 𝐺 − 𝐺𝑃
=
𝛽𝑔𝑖 𝛽𝑔

Re arreglando;

𝐺𝛽𝑔 = 𝐺𝛽𝑔𝑖 − 𝐺𝑃 𝛽𝑔𝑖

Dividiendo por G;

𝐺𝑃
𝛽𝑔 = 𝛽𝑔𝑖 − 𝛽
𝐺 𝑔𝑖

Puesto que el factor volumétrico es definido por;

𝑝𝑇
𝛽𝑔 = 𝑐𝑡𝑒
𝑧

Dicha constante incluye la temperatura del yacimiento que para un ejemplo dado se mantiene
invariable luego, la Ec. (5.2) se convierte en:

𝑝 𝑃𝑖 𝑃𝑖 𝐺𝑃
𝑐𝑡𝑒 = 𝑐𝑡𝑒 − 𝑐𝑡𝑒
𝑧 𝑧𝑖 𝑧𝑖 𝐺

9
 Figura N°1: P/Z

𝑷 𝒑𝒊 𝒑𝒊
= 𝑮𝑷 +
𝒛 𝒛𝒊 𝑮 𝒛𝒊

Los resultados al graficar dichos datos se ajustarán a una línea recta si el yacimiento es
volumétrico, sin embargo, dos casos pueden ocurrir, como se muestra en la Fig.

CASO 1:

5. Datos erróneos

6. Condensación retrógrada

7. Drenaje o escape de fluidos del yacimiento

8. Se está drenando más de la cuenta

9. Puede existir un yacimiento con presiones anormales (sobrepresionado)

10. Existe caída de presión a medida que hay reducción de permeabilidad y se pierde
comunicación

CASO 2:

11. Datos erróneos

12. Subsidencia

13. Hay comunicaciones o escapes hacia el yacimiento a lo largo de fallas o problemas

operacionales

10
Si se extrapola la recta en la figura anterior a P/Z = 0 (ó 14.7) se obtiene G. Observe que la
pendiente es;

𝒑𝒊
𝒎=
𝒛𝒊 𝑮

4. Ecuación general de balance de materia para yacimientos de petróleo

Normalmente está escrita en bases volumétricas: “los hidrocarburos iniciales = hidrocarburos
remanentes + agua que ha intruido”

• El yacimiento es considerado un tanque con cero dimensión

• Herramienta simple - si se necesitan más respuestas se recurre a un método más sofisticado.

La Fig. 2 muestra un tanque (yacimiento) con capa de gas a condiciones iniciales. Allí, m es el
tamaño de la capa de gas, ésta se obtiene a partir de registros, datos de corazones, datos de
completamiento, presiones de fondo y mapas estructurales. También se ha presentado una
disminución en el volumen poroso disponible a hidrocarburos causados por expansión de roca
y fluido.

Fig. 2. Esquema del balance de materia en yacimientos de petróleo

11
La ecuación de balance de materia puede escribirse como:

El volumen original de hidrocarburos + volumen remanente de hidrocarburos + el volumen de

agua intruida desde un acuífero + la reducción del volumen poroso debido a la expansión de
roca y fluido

La ecuación de balance de materia se expresa en unidades de yacimiento. Los términos que en

ella intervienen se definen como sigue:

Np=Petróleo producido, BF
N = Petróleo original in-situ, BF
G = Gas inicial en el yacimiento
m = Tamaño inicial de la capa de gas o volumen inicial de la capa de gas/volumen de la zona
de petróleo (N)
Np = Petróleo producido acumulado, BF
Gp = Gas producido acumulado, pcn
Wp = Gas producido acumulado, BF
Rp = Relación gas-petróleo acumulada, Gp/Np, pcn /BF
Rs = Relación gas-petróleo, pcn/BF
ßo, ßw = Factor volumétrico de formación del petróleo y del agua, bbl/BF
ßg = Factor volumétrico de formación del gas, bbl/pcn
We = Intrusión acumulada de agua, BF
Sw = Saturación de agua, fracción
cw, co, cg = Compresibilidad del agua, del petróleo y de gas, 1/psi
cf = Compresibilidad del volumen poroso, 1/psi
P = Presión estática del yacimiento, psia
∆P = Pi - P
i = inicial

Volumen de hidrocarburos iniciales = Volumen original de aceite + volumen original de

gas (en la capa), bbl

12
Volumen original de aceite = N ßoi, bbl
Volumen original de gas = mN ßoi, bbl

De modo que:

Volumen de hidrocarburos iniciales = N ßoi (1 + m)

Volumen de hidrocarburos remanentes = capa de gas original expandida + petróleo

remanente + gas neto liberado, bbl

Asumiendo que no se produce gas de la capa de gas. Existen problemas cuando el yacimiento
es delgado y cuando existe conificación.

β𝑔
Capa de gas original expandida = mNβ 𝑜𝑖 , bbl
β 𝑔𝑖

Petróleo remanente = (N − N 𝑃 )β 𝑜, 𝑏𝑏𝑙

Gas neto liberado = Gas liberado - gas libre producido

gas liberado = 𝑁(R 𝑠𝑖 − R 𝑠 )β 𝑔, 𝑏𝑏𝑙

gas libre producido = N 𝑃 (R 𝑃 − R 𝑠 )β 𝑔, 𝑏𝑏𝑙

Gas en solución producido = N 𝑃 R 𝑠, 𝑏𝑏𝑙

Gas neto liberado = 𝑁(R 𝑠𝑖 − R 𝑠 )β 𝑔, − N 𝑃 (R 𝑃 − R 𝑠 )β 𝑔, 𝑏𝑏𝑙

𝑔 β
𝐕𝐨𝐥𝐮𝐦𝐞𝐧 𝐡𝐢𝐝𝐫𝐨𝐜𝐚𝐫𝐛𝐮𝐫𝐨𝐬 𝐫𝐞𝐦𝐚𝐧𝐞𝐧𝐭𝐞𝐬 = (N − N 𝑃 )β 𝑜 + mNβ 𝑜𝑖 β + 𝑁 (R 𝑠𝑖 −
𝑔𝑖

R 𝑠 )β 𝑔 − N 𝑃 (R 𝑃 − R 𝑠 )β 𝑔, 𝑏𝑏𝑙

𝐈𝐧𝐭𝐫𝐮𝐬𝐢ó𝐧 𝐧𝐞𝐭𝐚 𝐝𝐞 𝐚𝐠𝐮𝐚 = (𝐖 𝒆 − 𝐖 𝑷 )𝛃 𝒘

Dake define el decremento del volumen poroso disponible a hidrocarburos debido a la expansión
de roca y fluido que ocurre por el decremento de presión a partir del concepto de compresibilidad
como:

13
 1 dV c𝑤 S𝑤 + c𝑓
C𝑒 =− − 𝑑(𝐻𝐶𝑉𝑃) = (1 + 𝑚)𝑁β 𝑜𝑖 ( ) 𝛥𝑃
V dP 1 − S𝑤

Involucrando todos los términos en mención, la ecuación de balance de materia queda:

β
(1 + 𝑚)𝑁β 𝑜𝑖 = (N − N 𝑃 )β 𝑜 + mNβ 𝑜𝑖 ( 𝑔 ) + 𝑁(R 𝑠𝑖 − R 𝑠 )β 𝑔 − N 𝑃 (R 𝑃 − R 𝑠 ) +
β 𝑔𝑖

c𝑤 S𝑤 +c𝑓
(W 𝑒 − W 𝑃 )β 𝑤 + (1 + 𝑚)𝑁β 𝑜𝑖 ( ) 𝛥𝑃
1−S𝑤

Agrupando los términos con N al lado izquierdo y con Np al lado derecho, se tiene:

β𝑔 c𝑤 S𝑤 +c𝑓
N [β 𝑜 − β 𝑜𝑖 + 𝑚β 𝑜𝑖 ( − 1) + β 𝑔 (R 𝑠𝑖 − R 𝑠 ) + β 𝑜𝑖 (1 + 𝑚) ( ) 𝛥𝑃] = N 𝑃 [β 𝑜 +
β 𝑔𝑖 1−S𝑤

β 𝑔 (R 𝑃 − R 𝑠 )] − (W 𝑒 − W 𝑃 )β 𝑤

Despejando N;

N 𝑃 [β 𝑜 + β 𝑔 (R 𝑃 − R 𝑠 )] − (W 𝑒 − W 𝑃 )β 𝑤
N =
β𝑔 c𝑤 S𝑤 + c𝑓
β 𝑜 − β 𝑜𝑖 + 𝑚β 𝑜𝑖 (β − 1) + β 𝑔 (R 𝑠𝑖 − R 𝑠 ) + β 𝑜𝑖 (1 + 𝑚) ( 1 − S ) 𝛥𝑃
𝑔𝑖 𝑤

5. Otros usos de la ecuación de balance de materia

1) Desarrollar una relación entre los fluidos producidos con la presión para predecir el
comportamiento del yacimiento

2) Verificar la existencia de la capa de gas

3) Comparar con el petróleo calculado volumétricamente. Esto no significa que los valores sean
cercanos. Balance de materia siente la presencia de una falla (2 yacimientos) o puede haber una
extensión en el yacimiento que el método volumétrico no detecta.

4) No es bueno tratar de hallar N con la Ecuación de balance de materia en yacimientos con

alta intrusión de agua porque Pi - P es pequeño.

La ecuación de balance de materia normalmente se parte en dos. Puede ser significativo.

El efecto de la disminución del volumen poroso disponible a hidrocarburos es despreciable

cuando se está por debajo del punto de burbuja.

14
6. Lineación de la ecuación de Balance de Materia mediante el método de Havlena y
Odeh
6.1. Predicción del comportamiento y recobro final de un yacimiento
6.1.1. Intrusión de agua
Algunas indicaciones de la existencia del empuje o intrusión de agua son:

Existe una zona subyacente de agua

Existe suficiente permeabilidad para soportar el movimiento de agua, usualmente mayor
a 50 md.
Aumento de la producción de agua a medida que transcurre el tiempo
Balance de materia es el mejor indicador para detectar la intrusión de agua

Existen tres métodos básicos para predecir y/o estimar la cantidad de agua que intruye al
yacimiento. Ellos son: El método de Schilthuis, el método de Hurst y Van Everdingen y el
método de Fetckovick.

6.1.2. Modelo de estado estable (Schilthuis)

Es el método más simple que se utiliza para fines de estimación de la intrusión de agua. Se
prefiere usar al principio. Se asume acuífero gigante y altamente permeable (Mayor de 50
md) para que la presión nunca caiga. Este método integra la ley de Darcy, la cual es estado
estable. Schilthuis arranca con la ley de Darcy:

Q W=CS (Pi−P)

Figura N°3: Esquema de flujo estable

15
La intrusión de agua resulta como la sumatoria de:

𝑡
𝑊𝑒 = 𝐶𝑠 ∫ (𝑃𝑖 − 𝑃)𝑑𝑡
0

𝑛
1
(𝑊𝑒 )𝑛 = 𝐶𝑠 ∑ [𝑃𝑖 − (𝑃𝑗−1 + 𝑃𝑗 )] ∆𝑡𝑗
2
𝑗=1

Siendo;

Qw = flujo de agua a través del WOC

Cs = Constante del acuífero. Contiene las constantes de Darcy: k, µ, A y sus unidades son
bbl/tiempo/psia.

Pi = Presión del acuífero

P = Presión estática del yacimiento.

Pj = presión estática a tj

Pj-1 = Presión estática entre tj-1 y tj.

La ecuación de balance de materia es:

𝑁𝑃 [𝐵° + 𝐵𝑔 (𝑅𝑃 − 𝑅𝑆 )] + 𝑊𝑃 ∗ 𝐵𝑊 − 𝑊𝑒
𝑁=
𝐵𝑔 (𝑅𝑆𝑖 − 𝑅𝑆 ) − (𝐵°𝑖 − 𝐵° )

El término del denominador se le conoce como expansibidad, D. Para estimar si existe o no

intrusión de agua se toma la ecuación con We = 0 y se calcula una Na (petróleo in-situ
aparente) para cada periodo de tiempo (puede ser 4 o 5 años). Este valor de Na es cada vez
más creíble puesto que se dispone de más datos.

𝑁𝑃 [𝐵° + 𝐵𝑔 (𝑅𝑃 − 𝑅𝑆 )] + 𝑊𝑃 ∗ 𝐵𝑊
𝑁=
𝐷

16
 Figura N°4: Método de Schilthuis Figura N°5: Determinación de intrusión de agua

Si no hay intrusión de agua el Na permanece constante para todos los periodos. Pero si existe
intrusión de agua el Na se incrementa en cada periodo. El método de Schilthuis se usa con la
ecuación de balance de materia para:

Calcular el petróleo original, N

Hallar la constante del acuífero, Cs, la cual relaciona la rata de intrusión de agua por
psi de caída de presión en el contacto agua-petróleo.

𝑊𝑒
Si 𝑁𝑎 = 𝑁 + , entonces:
𝐷

∆𝑝∆𝑡
𝑁𝑎 = 𝑁 + 𝐶𝑠 (∑ )
𝐷

Si se asume que N y Cs son constantes, la ecuación se comporta linealmente. Esto es, si se

grafica Na contra Σ(∆P∆t/D) resultará el comportamiento mostrado en la Fig. 5.12. El mismo
tratamiento es hecho para gas, pero D’ es diferente:

𝐷 " = 𝐵𝑔 − 𝑅𝑆𝑖

Figura N°6: Gráfico de Na vs. Σ(∆P∆t/D)

17
 Si no da línea recta es porque el modelo no es el apropiado. Si es aproximadamente recta el
modelo sirve.

𝐺𝑝 ∗ 𝐵𝑔 + 𝑊𝑃 + 𝐵𝑤
𝐺𝑎 =
𝐵𝑔 − 𝐵𝑔𝑖

Si el yacimiento es de gas o petróleo hay dos formas para obtener la “mejor” línea
recta: Al mejor criterio del ingeniero
Mínimos cuadrados

6.1.3. Modelo de estado inestable (Hurst and Van Everdingen)

Se asume:

Existe flujo radial de agua a través del acuífero

La caída de presión es constante a través del acuífero durante todo el tiempo
Las propiedades del acuífero son constantes y uniformes

Figura N°7: Esquema de estado inestable

Basado en estas consideraciones la ecuación es:

2 [
𝑊𝑒 = 2𝜋𝛼ℎ∅𝑐𝑒 𝑟𝑦𝑡𝑜 ∆𝑃 𝑄𝑡𝐷 ]

α = 0 − 1. Representa la extensión a la cual el acuífero rodea el yacimiento

We = Cantidad de agua intruida en cm3

h = Espesor del estrato en cm

ce = Compresibilidad efectiva del acuífero, 1/atm.

18
ryto = Radio del yacimiento de petróleo o gas, cm

∆P = Caída de presión constante a través del acuífero, atm

QtD = Función acumulativa de entrada de agua

tD = Tiempo adimensional La ecuación para tiempo adimensional es:

La ecuación de intrusión de agua en unidades de campo es:

2 [
𝑊𝑒 = 1.119𝛼ℎ∅𝑐𝑒 𝑟𝑦𝑡𝑜 ∆𝑃 𝑄𝑡𝐷 ]

Este método se basa en la solución de la ecuación de difusividad. Aplica el principio de

superposición que indica que las soluciones son aditivas. El principio de superposición al
contacto agua petróleo original es mostrado en la siguiente figura. Para usar superposición,
la curva se aproxima a una serie de pequeños incrementos de presión de modo que los pasos.

∆𝑃1 = 𝑃° − 0.5(𝑃° + 𝑃1 ) = (𝑃° + 𝑃1 )

∆𝑃2 = 0.5(𝑃° + 𝑃1 ) − (𝑃1 + 𝑃2 ) = 0.5 (𝑃° + 𝑃1 )

∆𝑃3 = 0.5(𝑃1 + 𝑃2 ) − 0.5 (𝑃2 + 𝑃3 ) = 0.5 (𝑃1 + 𝑃3 )

∆𝑃4 = 0.5(𝑃2 + 𝑃3 ) − 0.5 (𝑃3 + 𝑃4 ) = 0.5 (𝑃2 + 𝑃4 )

Generalizando;

1
∆𝑃𝑗 = (´𝑃 − 𝑃𝑗 )
2 𝑗−2

Para j = 1

1
∆𝑃1 = (𝑃 − 𝑃𝑗 )
2 °

19
Cuando la ecuación de Hurst and Van Everdingen se arregla para permitir variaciones de
presión queda:

Figura N° 9: Método de Hurst y van Everdingen

𝑊𝑒 = 𝐶𝑉 ∑ ∆𝑃𝑗 𝑄[𝐴(𝑡𝑛 − 𝑡𝑗−1 )]

𝐽=1

0.00633𝑘
𝐴= 2
∅𝑢𝑤 𝑐𝑒 𝑟𝑦𝑡𝑜

Si los intervalos de tiempo son igualmente espaciados:

𝑛

𝑊𝑒 = 𝐶𝑉 ∑ ∆𝑃𝑗 𝑄[𝐴(𝑛 − 𝑗 + 1)∆𝑡]

𝐽=1

Ahora, puesto que Na = N + We / D

𝐶𝑉 ∑𝑛𝐽=1 ∆𝑃𝑗 𝑄[𝐴(𝑛 − 𝑗 + 1)∆𝑡]

𝑵𝒂 = +
𝑫

Esta ecuación tiene tres incógnitas a saber: N, Cv y A. Estas dos últimas son las constantes
de Hurst and Van Everdingen. Los valores de la entrada de agua adimensional se encuentran
tabulados.

20
 Figura N°9: Determinación de A

Cálculo de N, Cv y A

1) Si no se especifican las dimensiones del yacimiento y del acuífero, asuma comportamiento

infinito. Lea los 𝑄𝑡𝐷

2) Calcule un valor aproximado de A. Use datos similares a los del yacimiento.

3) Escoja un incremento del tiempo, ∆t, normalmente 365 días.

4) Evalúe j A ∆t. Siendo j un entero sucesivo.

5) Evalúe 𝑄[(𝑗)(𝐴)(∆𝑡)] para cada j A ∆t

6) Evalúe ∆𝑃𝑗 para cada cambio de presión

7) Evalúe ∆𝑃𝑗 𝑄[(𝑗)(𝐴)(∆𝑡)]

8) Divida el valor obtenido en el paso 7 por D

9) Calcule 𝑁𝑎 𝑜 𝐺𝑎

10) Grafique 𝑁𝑎 vs. El resultado del paso 8. Ver esquema en Fig. 5.14.

11) Si el valor de A asumido fue correcto entonces se obtendrá una línea recta. En caso contrario,
estime un nuevo valor de A y repita los pasos 3 a 10. Cuando se obtenga una línea cercanamente

21
recta entonces A es cercanamente correcta. N es el intercepto y Cv es la pendiente. Como se
muestra en la siguiente figura. 5.10.3.

7. Método de Fetkovich para Acuíferos Finitos

Los resultados de este modelo se aproximan bastante a los del método de Hurst y Van
Everdingen para acuíferos finitos. Sin embargo, la teoría de Fetkovich es más simple y su
aplicación más fácil. Se utiliza un procedimiento de ensayo y error para evaluar las constantes
de intrusión de agua con datos de producción. El índice de productividad es adecuado para
describir la intrusión de agua desde un acuífero finito a un yacimiento de hidrocarburos; la rata
de intrusión de agua es directamente proporcional a la caída de presión entre la presión promedia
del acuífero y la presión en el contacto agua-aceite:

𝛿𝑊
𝑞𝑤 = 𝐽(𝑃̅𝑎 − 𝑃𝑓 ) = 𝛿𝑡𝑒 (1)

Donde:

𝑞𝑤 , 𝐵𝑃𝐷 Rata de intrusión de agua

𝐵𝑃𝐷
𝐽, 𝐼𝑛𝑑𝑖𝑐𝑒 𝑑𝑒 𝑝𝑟𝑜𝑑𝑢𝑐𝑡𝑖𝑣𝑖𝑑𝑎𝑑 𝑑𝑒𝑙 𝑎𝑐𝑢𝑖𝑓𝑒𝑟𝑜
𝑝𝑠𝑖

𝑃̅𝑎 , 𝑝𝑠𝑖𝑎 𝑃𝑟𝑒𝑠𝑖𝑜𝑛 𝑝𝑟𝑜𝑚𝑒𝑑𝑖𝑎 𝑑𝑒𝑙 𝑎𝑐𝑢𝑖𝑓𝑒𝑟𝑜

𝑃𝑓, 𝑝𝑠𝑖𝑎 𝑃𝑟𝑒𝑠𝑖𝑜𝑛 𝑒𝑛 𝑒𝑙 𝑐𝑜𝑛𝑡𝑎𝑐𝑡𝑜 𝑎𝑔𝑢𝑎 − 𝑎𝑐𝑒𝑖𝑡𝑒

𝑊𝑒, 𝑏𝑏𝑙 𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑎𝑐𝑢𝑚𝑢𝑙𝑎𝑑𝑜 𝑑𝑒 𝑖𝑛𝑡𝑟𝑢𝑠𝑖𝑜𝑛 𝑑𝑒 𝑎𝑔𝑢𝑎

𝑡, 𝑑𝑖𝑎𝑠 𝑇𝑖𝑒𝑚𝑝𝑜

Usando consideraciones de compresibilidad, la intrusión acumulada de agua es proporcional a

la caída de presión total en el acuífero:

𝑊𝑒 = 𝑉𝑎 𝑐𝑒 (𝑃𝑖 − 𝑃̅𝑎 ) (2)

𝑉𝑎 , 𝐵𝑃𝐷 𝑉𝑜𝑙𝑢𝑚𝑒𝑛 𝑖𝑛𝑖𝑐𝑖𝑎𝑙 𝑑𝑒 𝑎𝑔𝑢𝑎 𝑒𝑛 𝑒𝑙 𝑎𝑐𝑢𝑖𝑓𝑒𝑟𝑜

22
 1
𝑐𝑒 , 𝐶𝑜𝑚𝑝𝑟𝑒𝑠𝑖𝑏𝑖𝑙𝑖𝑑𝑎𝑑 𝑒𝑓𝑒𝑐𝑡𝑖𝑣𝑎 𝑑𝑒𝑙 𝑎𝑐𝑢𝑖𝑓𝑒𝑟𝑜, 𝑐𝑤 + 𝑐𝑓
𝑝𝑠𝑖

𝑃𝑖 , 𝑝𝑠𝑖𝑎 𝑃𝑟𝑒𝑠𝑖𝑜𝑛 𝑖𝑛𝑖𝑐𝑖𝑎𝑙 𝑑𝑒𝑙 𝑎𝑐𝑢𝑖𝑓𝑒𝑟𝑜

Derivando la Ec. (2) y combinando el resultado con la Ec. (1), resulta:

𝜕𝑃̅𝑎
−𝑉𝑎 𝑐𝑒 = 𝐽(𝑃̅𝑎 − 𝑃𝑓 )
𝜕𝑡

Si se asume que Pf permanece constante a lo largo del intervalo de interés, separando variables
e integrando, se tiene:

𝐽𝑡
−( )
𝑃̅𝑎 = 𝑃𝑓 + (𝑃𝑖 − 𝑃𝑓 )𝑒 𝑉𝑎 𝑐𝑒

Separando variables e integrando:

𝐽𝑡
−( )
𝑊𝑒 = (𝑃𝑖 − 𝑃𝑓 )𝑉𝑎 𝑐𝑒 [1 − 𝑒 𝑉𝑎 𝑐𝑒 ]

Como se asumió 𝑃𝑓 constante sobre el periodo de interés, el principio de superposición debe

usarse para manejar cambios de presión en el contacto agua-petróleo. Una forma incremental
de la ecuación anterior puede usarse sin superposición. Esta modificación es:

𝐽𝛥𝑡𝑛
−( )
𝛥𝑊𝑒,𝑛 = (𝑃̅𝑎,𝑛−1 − 𝑃̅𝑓,𝑛 )𝑉𝑎𝐶 𝑒 [1 − 𝑒 𝑉 𝐶
𝑎 𝑒 ] (5.10)

donde: Intrusión de agua incremental al paso de tiempo n ( 𝑡𝑛 − 𝑡𝑛−1 ) Presión promedia

del acuífero al final del paso de tiempo n - 1 el cual también es el arranque del
𝛥𝑊𝑒,𝑛, 𝑏𝑏𝑙
paso de tiempo n Presión promedia del contacto agua-aceite durante el paso de
𝑃̅𝑎,𝑛−1 , 𝑝𝑠𝑖𝑎 tiempo n𝑃̅𝑓,𝑛 = (𝑃𝑓,𝑛−1 − 𝑃𝑓,𝑛 )/2

𝑃̅𝑓,𝑛 , 𝑝𝑠𝑖𝑎 Tamaño del paso de tiempo (n)

𝛥𝑡, 𝑑𝑖𝑎𝑠

Si 𝛥𝑡𝑛 = 𝛥𝑡 para todos los pasos de tiempo, entonces la expresión:

23
 𝐽𝛥𝑡 𝐽𝛥𝑡𝑛
−( ) −( )
𝑉𝑎𝐶 𝑒 𝑉𝑎𝐶 𝑒
𝑉𝑎𝐶 𝑒 [1 − 𝑒 ] = 𝑉𝑎𝐶 𝑒 [1 − 𝑒 ] = 𝑐𝑡𝑒

Entonces la ec. 5.10 se transforma en:

𝛥𝑊𝑒,𝑛 = 𝐴(𝑃̅𝑎,𝑛−1 − 𝑃̅𝑓,𝑛 ) (5.11)

siendo A:

𝐽𝛥𝑡
−( )
𝑉𝑎𝐶 𝑒
𝐴 = 𝑉𝑎𝐶 𝑒 [1 − 𝑒 ]

La Ec. 5.8 puede rescribirse como:

∑𝑛−1
𝑗=1 𝛥𝑊𝑒,𝑗
𝑃̅𝑎,𝑛−1 = 𝑃𝑖 [1 − ] (5.12)
𝑉𝑎𝐶 𝑒 𝑃𝑖

La presión promedia en el contacto agua-petróleo durante el paso de tiempo n es:

𝑝̅𝑓,𝑛 = (𝑝𝑓,𝑛−1 + 𝑝𝑓,𝑛 )/2

La Ec. 5.11 proporciona la presión promedia del acuífero al final del paso de tiempo n-1.
Escribiendo una ecuación similar para la presión promedia del acuífero al final del paso de
tiempo n:

∑𝑛
𝑗=𝑖 ∆𝑊𝑖,𝑗
𝑝̅𝑎,𝑛 = 𝑝𝑖 [1 − ] (5.13)
𝑣𝑎 𝑐𝑒𝑝𝑖

Sustituyendo la Ec. 5.11 en 5.13, se tiene:

𝑛
1
𝑝̅𝑎,𝑛 = 𝑝𝑖 − ∑ ∆𝑊𝑒,𝑗
𝑣𝑎 𝑐𝑒
𝑗=𝑖

𝑛
1
𝑝̅𝑎,𝑛 = 𝑝𝑖 − ∑ 𝐴(𝑝̅𝑎,𝑗−1 − 𝑝̅𝑓,𝑖 )
𝑣𝑎 𝑐𝑒
𝑗=𝑖

24
 𝑛 𝑛

𝑝̅𝑎,𝑛 = 𝑝𝑖 − 𝐵 (∑ 𝑝̅𝑎,𝑗−1 − ∑ 𝑝̅𝑓,𝑖 )

𝑗=𝑖 𝑗=𝑖

donde:

𝐴
𝐵=
𝑣𝑎 𝑐𝑒

𝐽∆
−( 𝑡 )
𝐵 =1−𝑒 𝑣𝑎 𝑐𝑒

Rescribiendo la Ec. 5.8 al nivel de tiempo n, se tiene:

𝑊𝑒,𝑛 = 𝑣𝑎 𝑐𝑒 (𝑝𝑖 − 𝑝̅𝑎,𝑛 )

Luego, la intrusión de agua incremental se puede hallar de:

∆𝑊𝑒,𝑛 = 𝑊𝑒,𝑛 − 𝑊𝑒,𝑛−1

Tipo de acuífero Flujo radial Flujo lineal

Finito – no flujo 𝜃 0.003381𝑘𝑤ℎ
0.00708𝑘ℎ (360) 𝐽=
𝐽= 𝜇𝐿
𝑟
𝜇[ln (𝑟 𝑎𝑐 ) − 0.75]
𝑦𝑡𝑜

Finito – Presión constante 𝜃 0.001127𝑘𝑤ℎ

0.00708𝑘ℎ (360) 𝐽=
𝐽= 𝑟 𝜇𝐿
𝜇ln(𝑟 𝑎𝑐 )
𝑦𝑡𝑜

K en md, L=Longitud del acuífero lineal

Tabla N°1: índices de productividad

Como puede apreciarse en la tabla 1 el índice de productividad requiere condiciones de estado

pseudoestable. Luego, este método ignora los efectos del periodo transitorio en los cálculos de
intrusión de agua lo que obviamente introduce errores en los cálculos. Sin embargo, como se
hace por pasos, el método da resultados cercanos al de Hurst y van Everdingen.

25
CONCLUSIONES
Con base en los resultados de investigación recopilado de este trabajo, se pueden presentar las
conclusiones siguientes:

El objetivo principal de este trabajo ha sido presentar la aplicación de la ecuación de Balance de

Materia es una alternativa muy confiable para la ecuación del volumen original de un yacimiento.
También se llegó a la conclusión que este método normalmente es aplicable en etapas avanzadas
de explotación, donde el yacimiento presenta caídas de presión (ΔP) originadas por la extracción
de fluidos del yacimiento (aceite, gas y agua).
Se consiguió realizar dos ejercicios aplicando las ecuaciones que estiman el volumen de los fluidos
de reservorio con el método de Schilthuis normalmente conocida como método del balance de
materia.

RECOMENDACIONES
Se recomienda hacer un control de calidad en la información antes de proceder con los cálculos
de la EBM.
Para plantear estas ecuaciones de Balance de Materia se requiere que la información sea de
calidad y que haya poca dispersión en los datos, ya que la falta de esta cualidad en la
información puede llegar a ofrecer resultados indeseados, ya que la correlación es iterativa.
Para una mejor explotación del yacimiento y una buena administración del mismo, se requiere
conocer el comportamiento y la interacción con su entorno del mismo, para prevenir errores en
su explotación y así determinar cuál es el mejor método de recuperación para ese yacimiento.

26
BIBLIOGRAFIA
Bastidas C, F. (08 de Noviembre de 2013). Modelo de simulacion dinamico del yacimiento G-8 GF-
5, campo Guafita, Area Sur, Cuenca Barinas- Apure. Recuperado el 28 de Junio de 2022,
de http://saber.ucv.ve/bitstream/10872/5788/1/T.E.G.pdf
Escobar Macualo, F. H. (2002). FUNDAMENTOSDE INGENIERIADE YACIMIENTOS. Neiva -
Huila - Colombia: Editorial Universisad Surcolombiana.
BaGarcía., M. G. (2 de Octubre de 2012). Caracterización Dinámica de Yacimientos. Obtenido de
https://www.academia.edu/6799478/Caracterizaci%C3%B3n_Din%C3%A1mica_de_Ya
cimientos_de_Hidrocarburos

27
ANEXOS
EJERCICIO #1 (APLICANDO PARA UN RESERVORIO DE GAS)
Dados:

P Gp*MMM pcn Componente X

4000 0 C1 86 %
3825 34 C2 7%
3545 75 C3 3%
3388 88 C4 2%
3192 120 C5 2%
2997 145
2738 192 T=618 °R

Estime el gas inicial.

Yacimiento de gas no volumétrico

We es diferente de cero. Existen producción a la expansión y/o empuje de gas.

𝑮𝑷 𝜷𝒈𝒊 = 𝑮[𝜷𝒈𝒊 − 𝜷𝒈 ] + 𝑾𝒆 − 𝑾𝑷 𝜷𝒘

Las incógnitas son G y We. Asuma que:

F: extracciones

E: expansiones

Entonces;

𝑭 = 𝑮𝑷 𝜷𝒈𝒊 + 𝑾𝑷 𝜷𝒘

𝑬 = 𝜷𝒈𝒊 − 𝜷𝒈 , 𝒍𝒖𝒆𝒈𝒐;

𝑭 = 𝑮𝑬𝒈 + 𝑾𝒆

De donde;

28
 𝑬 𝑾𝒆
= 𝑮+ ↔ 𝒚 = 𝒎𝒙 + 𝒃
𝑬𝒈 𝑬𝒈

Para solucionar y conocer las incógnitas se utiliza el método de Hurst & Van Everdingen quienes
resolvieron la ecuación de difusividad para dos casos:

1) Presión constante (cálculos de We)

2) Caudal constante (análisis de presiones de fondo)

Con el método de Hurst & Van Everdingen, se tiene una intrusión de agua adimensional dada
por;

𝑾𝒆𝑫 = 𝜷 ∑ ∆𝒑𝒋 ∗ 𝑸𝒕𝑫

Siendo QtD es el caudal adimensional obtenido por Hurst & Van Everdingen y ß es la constante
de intrusión de agua. Utilizando este nuevo término, se tiene:

𝑭 𝑾𝒆𝑫
= 𝑮+𝜷
𝑬𝒈 𝑬𝒈

Una vez conocido We se pueda estimar el valor de G.

Ecuación de Havlena y Odeh para un yacimiento no volumétrico

r = 1500 ft, φ = 25 %, ct = 5x10-5 /pisa, Pi = 3274 psia, viscosidad = 30 cp y permeabilidad

0.01 md. Los resultados de entrada y salida son:

29
 t, días P, psia ∆𝑷𝒊 ∆𝑷𝒋 tp 𝑸𝒕𝑫 ∑ ∆𝑷𝒋 𝑸𝒕𝑫

14 3094 180 90 2.36 2.719 244.71

28 2984 110 145 4.72 4.358 876.62
42 2889 95 102.5 7.08 5.789 1469.99
53.62 2819 70 82.5 9.04 7.107 2056.3

Donde;

∆Pi + ∆Pi−1
∆Pj =
2

0.006323kt
tD = 2
∅μct rw

Estando en unidades de campo. La constante de intrusión de agua, se calcula mediante:

2
𝜃
𝛽 = 1.119∅𝑐𝑒 𝑟𝑦𝑡𝑜 ℎ
360

Siendo θ el ángulo subtendido por la circunferencia del yacimiento semicircular, los demás
parámetros se dan en unidades de campo.

co So + cw Sw + cg Sg + cf
ce
1 − Sw

Puesto que se trata del acuífero; la ecuación para la compresibilidad efectiva queda:

𝒄𝒆 = 𝒄𝒘 + 𝒄𝒇

We también puede ser estimada por flujo fraccional, así:

𝒒𝒘 𝟏
𝒇𝒘 = =
𝒒𝒘 + 𝒒𝒐 𝒌𝒈 𝝁𝒘 𝜷𝒘
𝟏+
𝒌𝒘 𝝁𝒐 𝜷𝒐

30
EJERCICIO #2 (APLICANDO PARA UN RESERVORIO DE PETROLEO)
Calcule el petróleo in-situ para un yacimiento de petróleo cuya presión incial es de 5000psia y
su presión de burbuja es de 2750 psia. Otros datos son:
𝛽𝑜𝑖 = 1.305 𝑏𝑏𝑙/𝐵𝐹
Ø = 10%
𝐶𝑟 = 3𝐸 − 7 𝑝𝑠𝑖𝑎
𝛽𝑜𝑏 = 1.35 𝑏𝑏𝑙/𝐵𝐹
𝐶𝑜 = 15.3𝐸 − 6 𝑝𝑠𝑖𝑎
𝑇 = 240°𝐹
𝑆𝑤 = 21.6 %
𝑐𝑤 = 3.5𝐸 − 6 𝑝𝑠𝑖𝑎
𝑅𝑠𝑖 = 500 𝑝𝑐𝑛/𝐵𝐹
A una presión de 3350 psia se han producido 1510000 BF de petróleo y el factor volumétrico
es de 1.338 bbl/BF. A una presión de 1500 psia se han producido 6436000 BF de petróleo y
3732000000 de pcn de gas. El factor volumétrico es de 1.250 bbl/BF, el gas en solución es 375
bbl/BF y el factor de desviación del gas es de 0.9. Se considera que la intrusión de agua es cero.
Halle el petróleo original a 3350 psia y a 1500 psia. Realice los cálculos considerando cw y
expansión de roca. Efectúe los cálculos considerando la compresibilidad del petróleo. Para el
segundo caso parta la ecuación de balance de materia y no. Desprecie los efectos de cw y cf .
Solución: Primero que todo es necesario calcular cf , ya que el dato dado fue cr, la
compresibilidad de la matriz.
Ø = 1 − 𝑉𝑟 (5.3)
1 𝑑𝑉𝑟
𝐶𝑟 = − 𝑉 (5.4)
𝑟 𝑑𝑃

1 𝑑Ø
𝐶𝑓 = − Ø 𝑑𝑃 (5.5)

Derivando la Ec. (5.3) con respecto a P y sustituyendo este resultado en la Ec. (5.5)
1 𝑑𝑉
𝐶𝑓 = − Ø (− 𝑑𝑃𝑟 ) (5.6)

Despejando dVr/dP de la Ec. (5.4) y reemplazando ese resultado en la ecuación (5.6) se tiene:
1
𝐶𝑓 = − (𝐶𝑟 ∗ 𝑉𝑟 )
Ø
De la Ec. (5.3) Vr = 1 - φ, luego:

31
 1
𝐶𝑓 = − 𝐶𝑟 (1 − Ø) = 2.7𝐸 − 6 𝑝𝑠𝑖𝑎
Ø
Cuando P > Pb, se tiene m = 0, Rp = Rs, Gp = Rs*Np, Rs = Rsi y We = 0. La ecuación de
balance de materia queda entonces:
𝑁𝑝 𝛽𝑜
𝑁=
𝑐𝑤 𝑆𝑤 + 𝑆𝑓
𝛽𝑜 − 𝛽𝑜𝑖 + 𝛽𝑜𝑖 ( 1 − 𝑆 )∆𝑃
𝑤

a) Puesto que para el primer caso se desprecia cw y cf:

𝑁𝑝 𝛽𝑜 1510000(1.338)
𝑁= = = 61.2𝑀𝑀𝐵𝐹
𝛽𝑜 − 𝛽𝑜𝑖 1.338 − 1.305

b) Considerando cw y cf

𝑁𝑝 𝛽𝑜 1510000(1.338)
𝑁= = = 47.5𝑀𝑀𝐵𝐹
𝛽𝑜 − 𝛽𝑜𝑖 3.5𝐸 − 6(0.216) + 2.7𝐸 − 6
1.338 − 1.305 ( ) 1650
1 − 0.216
∆= 29%

c) Involucrando la compresibilidad del petróleo. Puesto que ßoi-ßo está en la ecuación de

balance de materia y el factor volumétrico de formación del petróleo puede determinarse
por:
𝛽𝑜 = 𝛽𝑜𝑖 [1 + 𝐶𝑜 (𝑃𝑖 − 𝑃)]

Haciendo este reemplazo, la ecuación de balance de materia queda:

𝑁𝑝 𝛽𝑜
𝑁= =
𝛽𝑜𝑖 𝐶𝑜 ∆𝑃 + 𝛽𝑜𝑖 𝐶𝑒 ∆𝑃

Definiendo:
𝐶𝑒′ = 𝐶𝑜 + 𝐶𝑒

Luego;

𝑁𝑝 𝛽𝑜
𝑁= = 47.6 𝑀𝑀𝐵𝐹
𝛽𝑜𝑖 𝐶𝑒′ ∆𝑃

d) Sin romper la ecuación de balance de materia y considerando los efectos de la expansión

de agua connata y roca:

32
 𝑇𝑧 (460 + 20)0.9
𝛽𝑔 = 0.00503 = (0.00503) = 0.00211 𝑏𝑏𝑙
𝑝 1500

𝑅𝑝 = 𝐺𝑝 + 𝑁𝑝 = 3732/6.436 = 580𝑝𝑐𝑛/𝐵𝐹. 𝐸𝑛𝑡𝑜𝑛𝑐𝑒𝑠;

6.436[1.25 + 0.0021(580 − 37)]

𝑁=
3.5 ∗ 0.216 + 2.7
1.25 − 1.305 + 0.00211(500 − 375) + 1.305 ( ) 𝑥10−6 ∗ 3500
1 − 0.216
= 47.3 𝑀𝑀𝑝𝑐𝑛

e) Partiendo la Ecuación de Balance de Materia

Por encima del punto de burbuja;
𝑁𝑝 𝛽𝑜
𝑁= =
𝛽𝑜𝑖 𝛽𝑜𝑒 ∆𝑃

De donde;

𝑁 ∗ 𝛽𝑜 ∗ 𝐶𝑜𝑒 ∆𝑃 (47.6𝑀𝑀)(1.305)(19.71𝐸 − 6)(5000 − 2750)

𝑁𝑝 = = = 2.04 𝑀𝑀𝐵𝐹
𝛽𝑜 1.35

𝑁𝑝 = 𝑁𝑝 ∗ 𝑅𝑠 = (2.04𝑀𝑀)(500) = 1.022𝑀𝑀𝑝𝑐𝑛

e) El punto de burbuja es considerado ahora el punto de arranque. El petróleo producido

desde el Punto de burbuja hasta la presión de 1500 psia es (6.436 - 2.04) MMBF = 4.396
MMBF.

El gas producido sería (3732 - 1020) MMpcn = 2712 MMpcn. Rp será = 2712/4.396 =
617 pcn/BF.
4.396𝑀𝑀[1.25 + 0.0021(617 − 375)]
𝑁@𝑃𝑏 = = 45.21 𝑀𝑀𝐵𝐹
3.5 ∗ 0.26 + 0.27
1.25 − 1.35 + 0.00211(500 − 375) + 1.35 ( ) 𝑥10−6 ∗ 1250
1 − 0.216

𝑁 = 𝑁𝑠𝑜𝑏𝑟𝑒 + 𝑁𝑏𝑎𝑗𝑜 𝑃𝑏 = 47.25 𝑀𝑀𝐵𝐹

f) Despreciando las compresibilidades del agua y volumen poroso.

4.396 𝑀𝑀 [1.25 + 0.00211(617 − 375)]

𝑁= = 47.26𝑀𝑀𝐵𝐹
1.25 − 1.35 + 0,00211(500 − 375)

Luego N=(2.04+47.26)MM=49.3, habiendo una diferencia del 4%

33