Apuntes de Crudos Pesados UDO PDF

Exploración de Crudos Pesados PREPARADOR: JUAN ZAMORA
Guía de práctica de yacimientos de

crudos pesados
UDO-ANZOATEGUI Página 1
Inyección cíclica de vapor
La inyección cíclica de vapor, básicamente consiste en inyectar vapor a un pozo de petróleo durante un determinado
tiempo, generalmente de una a tres semanas; cerrar el pozo por un corto período de tiempo (3 a 5 días), y luego
ponerlo en producción. La Figura 7.1., es una representación esquemática de un proceso de inyección cíclica de vapor.
Una vez que el pozo es puesto en producción, este producirá a una tasa aumentada durante un cierto periodo de
tiempo, que en general, puede ser del orden de 4 a 6 meses, y luego declinará a la tasa de producción original. Un
segundo ciclo de inyección puede emplearse, y de nuevo la tasa de producción aumentará y luego declinará. Ciclos
adicionales pueden realizarse de una manera similar, sin embargo, el petróleo recuperado durante tales ciclos será cada
vez menor.
Aunque existen variaciones del proceso de inyección cíclica descrito, es evidente que se trata básicamente de un
proceso de estimulación, usualmente utilizado para petróleos pesados (8-15 ºAPI), puede utilizarse también para
yacimientos de cualquier tipo de crudo.
 Mecanismos en la recuperación cíclica de vapor
Los mecanismos involucrados en la producción de petróleo durante la inyección cíclica de vapor son diversos y
complejos. Sin duda. La reducción de la viscosidad del crudo en la zona calentada cercana al pozo, afecta
ampliamente el comportamiento de la producción.
Después de inyectar el volumen deseado de vapor (expresado como barriles equivalentes de agua), el pozo se cierra
aproximadamente por dos semanas. Esto es llamado el periodo de remojo (“soak period”), y el propósito es promover
una condensación parcial de todo el vapor inyectado para calentar la roca y los fluidos, así como también, permitir la
distribución uniforme del calor.
Durante la inyección de vapor y los periodos de remojo, la viscosidad del petróleo es disminuida dentro de la zona de
vapor, ocurriendo expansión térmica del petróleo y del agua.
Debido a la presurización de la arena, cualquier gas libre es forzado en solución, así, el gas disuelto en el petróleo
durante el periodo de producción juega un papel importante en la expulsión del petróleo. Esta fuerza expulsiva debida
al gas en solución, será relativamente efectiva solo en el primer ciclo, ya que para el segundo ciclo, mucho de este gas
puede haber sido producido.
Otros factores que contribuyen a la recuperación del petróleo en inyección cíclica son: el efecto de la temperatura
sobre las permeabilidades relativas de los fluidos que fluyen, el efecto del calentamiento mas allá de la zona
contactada por el vapor, la imbibición del agua caliente en estratos de baja permeabilidad, resultando flujo de
petróleo a los estratos permeables y finalmente al pozo, y la compactación de la roca-yacimiento en caso de existir.
Inyección continua de vapor
La inyección continua de vapor consiste en inyectar vapor en forma continua a través de algunos pozos y producir el
petróleo por otros. Los pozos de inyección y producción se perforan en arreglos, tal como en la inyección de agua. En
la actualidad se conocen varios proyectos exitosos de inyección continua de vapor en el mundo, muchos de los cuales
fueron inicialmente proyectos de inyección cíclica, que luego se convirtieron a inyección continua en vista de las
mejoras perspectivas de recuperación: 6-15% para cíclica vs. 40-50% para continua.
La inyección continua de vapor difiere apreciablemente en su comportamiento de la inyección de agua caliente, siendo
esta diferencia producto únicamente de la presencia y efecto de la condensación del vapor de agua. La presencia de la
fase gaseosa provoca que las fracciones livianas del crudo se destilen y sean transportados como componentes
hidrocarburos en la fase gaseosa. Donde el vapor se condensa, los hidrocarburos condensables también lo hacen,
reduciendo la viscosidad del crudo en el frente de condensación. Además, la condensación del vapor induce un
proceso de desplazamiento más eficiente y mejora la eficiencia del barrido.
 Mecanismos en la recuperación continua de vapor
Cuando se inyecta vapor en forma continua en una formación petrolífera, el petróleo es producido por causa de tres
mecanismos básicos: destilación por vapor, reducción de la viscosidad y expansión térmica, siendo la destilación
por vapor el más importante. Otros fenómenos que contribuyen a la recuperación de petróleo son la extracción con
solventes, empuje por gas en solución y desplazamientos miscibles por efectos de la destilación por vapor.
El petróleo en la vecindad del extremo de inyección es vaporizado y desplazado hacia delante. Una cierta fracción del
petróleo no vaporizado es dejada atrás. El vapor que avanza se va condensando gradualmente, debido a las pérdidas de
calor hacia las formaciones adyacentes, generando así una zona o banco de agua caliente, el cual va desplazando
petróleo y enfriándose a medida que avanza, hasta finalmente alcanzar la temperatura original del yacimiento (ver
figura 8.1). Desde este punto en adelante el proceso de desplazamiento prosigue tal como en la inyección de agua fría.
Así, se puede observar que se distinguen tres zonas diferentes: la zona de vapor, la zona de agua caliente y la zona de
agua fría.
Aún queda por evaluarse la formación de CO2 ( y de otros gases en menores cantidades)14 resultante de las
reacciones entre el vapor y el crudo (o de cualquier otra fuente), proceso conocido como acuatermólisis, el cual
también puede actuar como mecanismo de desplazamiento.
Combustión in Situ
La combustión in situ implica la inyección de aire al yacimiento, el cual mediante ignición espontánea o inducida,
origina un frente de combustión que propaga calor dentro del mismo. La energía térmica generada por éste método da
lugar a una serie de reacciones químicas tales como oxidación, desintegración catalítica, destilación y polimerización,
que contribuyen simultáneamente con otros mecanismos tales como empuje por vapor y vaporización, a mover el
petróleo desde la zona de combustión hacia los pozos de producción.
Se conocen dos modalidades para llevar a cabo la combustión in situ en un yacimiento, denominadas: combustión
convencional hacia adelante (forward combustión) debido a que la zona de combustión avanza en la misma dirección
del flujo de fluidos; y combustión en reverso o contracorriente (reverse combustion) debido a que la zona de
combustión se mueve en dirección opuesta a la del flujo de fluidos. En la primera de ellas, se puede añadir la variante
de inyectar agua alternada o simultáneamente con el aire, originándose la denominada combustión húmeda, la cual a
su vez puede subdividirse dependiendo de la relación agua/aire inyectado, en combustión húmeda normal, incompleta
y superhúmeda. Las mismas persiguen lograr una mejor utilización del calor generado por la combustión dentro del
yacimiento, reduciendo así los costos del proceso.
 Combustión Convencional
En este proceso los fluidos inyectados y el frente de combustión se mueven en el mismo sentido, es decir, del pozo
inyector hacia los pozos productores. Durante este proceso se forman dentro del yacimiento varias zonas
perfectamente diferenciables1, las cuales se indican en detalle en la figura 9.1. Estas zonas se originan por las altas
temperaturas generadas dentro del medio poroso, el cual se encuentra saturado inicialmente con agua, petróleo y gas.
En la combustión convencional, la ignición se induce en el pozo inyector, y una vez lograda, la temperatura aumenta
gradualmente hasta que se alcanza el punto de vaporización del agua. El vapor de agua generado se mezcla con la
corriente de gases, y fluye a través del yacimiento a la misma tasa con la que se transfiere calor desde el frente de
combustión. A esta temperatura ocurre el desplazamiento por destilación, de parte del petróleo.
Una vez que toda el agua se ha vaporizado, la temperatura en este punto aumenta progresivamente y la viscosidad del
crudo in situ disminuye, así mismo los volúmenes de petróleo y gas aumentan por expansión térmica. Este mecanismo
resulta en un incremento del flujo de petróleo dentro de la corriente líquida. Puede ocurrir también vaporización del
petróleo cuando la presión de vapor de sus componentes exceda la presión del sistema. Cuando la temperatura supera
los 350 °C (límite que depende del tipo de crudo y de las características del yacimiento) ya la mayor parte del petróleo
ha sido desplazado de allí hacia las zonas menos calientes del yacimiento y en el material pesado depositado comienza
a operarse la desintegración catalítica, de la que se origina un volumen adicional de hidrocarburos livianos.
Finalmente, la parte más pesada del crudo (coque) se consume como combustible, alcanzándose la máxima
temperatura de combustión.
Las dimensiones de las zonas de altas temperaturas (combustión, deposición de coque, desintegración catalítica y
evaporación) son en realidad pequeñas2,3 y su movimiento dentro del yacimiento obedece principalmente a dos
mecanismos de transferencia de calor: la conducción a través de la matriz sólida del yacimiento y la convección por
los gases que fluyen desde la zona caliente hacia la zona inalterada. A continuación de esta zona se encuentra la zona
de condensación, la cual también se conoce como meseta de vapor pues en ella la temperatura es más o menos
constante y cercana en magnitud al punto de ebullición del agua, a la presión parcial del agua en la fase de vapor.
Los mecanismos de producción que actúan durante este proceso son muy variados, destacándose el empuje por gas,
los desplazamientos miscibles, la vaporización y la condensación. Estos mecanismos son auxiliados por importantes
reacciones, tales como la oxidación, destilación, desintegración catalítica y polimerización, las cuales ocurren
simultáneamente en las zonas de combustión, coque y desintegración catalítica indicadas en la figura 9.1. En estas
zonas ocurre también un incremento en la presión debido principalmente al aumento en el volumen de los fluidos por
expansión térmica, lo cual produce un aumento de la tasa de flujo hacia los pozos productores.
La combustión convencional se recomienda para yacimientos poco profundos, entre 200 y 5000 pies (limitación
impuesta principalmente por los costos de compresión e inyección de aire), para crudos cuya gravedad oscile entre 8°
y 26° API, pues ello garantiza suficiente deposición de coque para mantener activo el frente de combustión6. Sin
embargo, se ha estudiado la posibilidad de aplicar este proceso a yacimientos de crudos más livianos (hasta 40° API),
siempre y cuando se trate de crudos de base nafténica o parafínica.
 Combustión en Reverso
En la combustión en reverso, el frente de combustión se mueve en dirección opuesta al flujo de aire. La combustión se
inicia en el pozo productor y el frente de combustión se mueve contra el flujo de aire. La Figura 9.5., presenta un
esquema simplificado de este proceso, indicándose las zonas formadas dentro del yacimiento.
El movimiento del frente de combustión es hacia las zonas de mayor concentración de oxígeno, y los fluidos
atraviesan dicho frente de combustión como parte de la corriente de gas, siendo transportados a través de las zonas
calientes hacia los pozos de producción por drenaje por gravedad y empuje por gas. El comportamiento de este
proceso es muy diferente al convencional, pues la zona de combustión no consume todo el combustible depositado
delante de ella; no obstante, parte de los componente livianos y medianos del petróleo in situ son utilizados como tal.
Casi no existe producción de monóxido o bióxido de carbono y las principales reacciones ocurridas durante la
oxidación del crudo originan compuestos oxigenados tales como aldehídos, ácidos, peróxidos, y otros.
El petróleo producido tiene características diferentes al crudo in situ, pues es más liviano y de menor viscosidad20.
Por esta razón, este proceso ha sido propuesto para aplicarlo en áreas bituminosas y en crudo extrapesados (de
gravedad menor de 8 °API), los cuales son prácticamente inmóviles a condiciones de yacimiento21. En estos
yacimientos la combustión convencional fracasaría, pues los bancos de líquido formados delante del frente de
combustión originarían un bloqueo de la permeabilidad específica al gas, impidiendo su circulación.
En comparación con la combustión convencional, la combustión en reverso es menos atractiva debido a
recuperaciones más bajas, mayor consumo de combustible, mayores pérdidas de calor, problemas de equipo más
severos y problemas asociados a la ignición espontánea en la cercanía del pozo inyector.
 Combustión Húmeda
Esta variante de la combustión convencional se genera al inyectar agua alternada o simultáneamente con el aire, una
vez que se ha logrado la ignición del crudo in situ. La combustión húmeda se clasifica en combustión húmeda normal,
incompleta y superhúmeda. Esta clasificación ha sido propuesta en base a las características de los perfiles de
temperatura y saturación originados durante el proceso
En el caso de combustión húmeda, normal o incompleta el agua inyectada al ponerse en contacto con la zona quemada
se evapora y fluye a través del frente de combustión como parte de la fase gaseosa, puesto que la máxima temperatura
del frente de combustión es, en este caso, mayor que la temperatura de evaporación del agua a la presión del sistema.
El proceso se denomina húmeda normal cuando el coque depositado se consume completamente. Por el contrario,
cuando el agua inyectada hace que el combustible depositado no se queme por completo, entonces se trata de una
combustión húmeda incompleta.
La combustión superhúmeda se logra cuando la cantidad de calor disponible en la zona quemada, no es suficiente para
vaporizar toda el agua inyectada al sistema. En este proceso, la máxima temperatura de combustión desaparece, y la
zona de vaporización-condensación se esparce por todo el medio poroso. Este tipo también ha sido denominada
combustión parcialmente apagada
Las ventajas que ofrece la combustión húmeda en cualquiera de sus variantes, es que reduce en cierto grado los
problemas operacionales propios del proceso, entre los cuales se cuentan: la producción de arena por las altas tasas de
gases, la corrosión por ácidos en los pozos y en las facilidades de producción, la oxidación en las tuberías por la
presencia de oxígeno en la corriente de gases inyectados y producidos, la formación de emulsiones, el taponamiento
de yacimiento por el coque depositado y otros. Estos problemas disminuyen debido a las menores cantidades de gases
producidos, por la disolución de los ácidos producidos en el gran volumen de agua existente (inyectada, connata y
formada), por el consumo completo de oxígeno y por la reducción de las temperaturas generadas
.
 Criterios para la selección del yacimiento en un proceso de combustión in situ
 Contenido de petróleo: dado que el frente de combustión puede consumir alrededor de 300 Bl/acres.ft del
petróleo inicial, al menos 600 Bl/acres.ft de petróleo deben estar presentes en el yacimiento. Esto implica una
porosidad del orden del 20% y una saturación porcentual del petróleo del 40%.
 Espesor: el espesor de arena neta no debe exceder los 50 pies. Espesores mayores de 50 pie requerirán
suficiente inyección de aire para mantener el frente de combustión moviéndose al menos a una velocidad de
0,25 pies/días, lo que sería excesivo con respecto a las limitaciones prácticas impuestas por el equipo de
compresión.
 Profundidad: la profundidad del pozo debe ser mayor de 200 pies. En general profundidades menores de 200
pies, podrían limitar severamente la presión a la cual el aire puede ser inyectado. Operaciones en yacimientos
profundos resultan en pozos altamente costosos, como también en gastos sustanciales en la compresión del
aire, por lo que las condiciones económicas pueden imponer profundidades prácticas del orden de 2.500 a
4.500 pies.
 Gravedad y viscosidad del petróleo: en general, petróleos de gravedades mayores de 40 °API no depositan
suficiente coque (combustible) para mantener un frente de combustión. Por otro lado, petróleos de gravedades
menores de 8 °API son generalmente muy viscosos para fluir delante del frente de combustión cuando la
temperatura del yacimiento prevalece sobre la temperatura de combustión. La recuperación de petróleo de
gravedades extremadamente bajas pueden ser posibles por medio de la combustión en reverso, donde el
petróleo producido fluye a través de la zona calentada y su composición es estructuralmente alterada.
 Permeabilidad: cuando la viscosidad del petróleo es alta (un yacimiento conteniendo un petróleo de 10
°API), una permeabilidad mayor de 100 md podría ser necesaria, especialmente si el yacimiento es somero y
la presión de inyección es limitada. Un crudo de gravedad entre 30 y 35 °API a una profundidad de 2.500 pies,
puede responder a un proceso de combustión in situ, aún con permeabilidades tan bajas como de 25 a 50 md.
 Tamaño del yacimiento: el yacimiento debe ser lo suficientemente grande, ya que si una prueba piloto a
pequeña escala tiene éxito, un éxito económico a gran escala puede ser esperado. Dependiendo del espesor de
la arena, el tamaño del yacimiento, podría ser aproximadamente de 100 acres.
 Ignición
Se refiere al inicio de la combustión del petróleo (combustible) en el yacimiento. La ignición puede ser de dos formas:
espontánea y artificial.
 Ignición Espontanea
Tal como su nombre lo indica, ocurre naturalmente cuando al aumentar la temperatura por efectos de la presión de
inyección de aire, se inicia la combustión. El que ocurra ó no ignición espontánea depende principalmente del tipo de
crudo, o sea, de su composición. Un parámetro importante referente a la ignición espontánea, es el tiempo requerido
para que esta ocurra.
 Ignición Artificial
Consiste en utilizar un calentador a gas o eléctrico, o productos químicos para lograr la ignición. Comparada con la
ignición espontánea, tiene la desventaja del uso de calentador y de que solo parte del espesor total de la formación es
puesto en ignición
CAPITULO II. PROPIEDADES TÉRMICAS DE ROCAS Y FLUIDOS
Viscosidad del Petróleo

o Ecuación de Andrade
𝑏 𝑏
𝜇 = 𝑎𝑒 𝑇 ∗ 𝑜 𝐿𝑛 𝜇 = 𝐿𝑛 𝑎 +
𝑇∗
Donde:
µ: Viscosidad del petróleo, cp
T*: Temperatura en grados absolutos, R= (°F + 460) K= (°C +273,1)
a y b son constantes
o Carta ASTM de viscosidad – temperatura

𝑇∗
𝐿𝑜𝑔 Log 𝜈1 + 0.8 = −𝑛𝐿𝑜𝑔 + Log 𝜈1 + 0.8
𝑇1 ∗
Donde:
𝜇 10 −𝑛
𝜈= = 𝜈1 + 0.8 − 0.8
𝜌𝑜
𝑇
𝐿𝑜𝑔
𝑇1
𝑛= 𝐿𝑜𝑔 Log 𝜈1 + 0.8 − 𝐿𝑜𝑔 Log 𝜈2 + 0.8
𝑇
𝐿𝑜𝑔 2
𝑇1
𝜌𝑜𝑠𝑐
𝜌𝑜 =
𝑇 − 68
1+
1885
ν: Viscosidad cinemática a la temperatura T*, centistokes
𝜌𝑜 : Densidad a la temperatura deseada pero en T(°F), gr/cc
𝜌𝑜𝑠𝑐 : Densidad a condiciones normales, gr/cc
n: Constante
Si se conoce solamente un punto de la viscosidad del crudo pesado medida, es factible hacer la predicción, usando la
técnica de Walter, si se dispone de la pendiente de la línea de viscosidad.
𝐿𝑜𝑔 Log 𝜈 + 0.7 = 𝑛𝐿𝑜𝑔 𝑇 + 𝑏
𝑚 = 0.3249 − 0.4106 𝑏
𝑏 +𝑚𝑙𝑜𝑔 (𝑇)
𝜈 = 1010 − 0.7
Viscosidad de las mezclas liquidas
o Método de Cragoe
2995.73
𝐿=
𝐿𝑛 𝜇 + 7.6009
𝐿𝑚 = 𝑓1 𝐿1 + 𝑓2 𝐿2
Donde:
L: Licuicidad
𝑓: Fracción en peso
µ: Viscosidad, cps
Viscosidad del Agua y del Vapor
o Ecuación de Gottfried
1776 − 𝑇𝑠
𝜇𝑤 =
26.5 ∗ 𝑇𝑠 − 89
Donde:
µw: Viscosidad del agua, cp
Ts: Temperatura de Saturación, °F
Para vapor seco y saturado: Para vapor húmedo:
𝜇𝑠 = 0.2 ∗ 𝑇𝑠 + 81.97 𝑋10−4 𝜇𝑤𝑠 = 𝑋𝑆𝑡 ∗ 𝜇𝑠 + 1 − 𝑋𝑆𝑡 ∗ 𝜇𝑤
𝜇𝑤𝑠 = 𝜇𝑤 − 𝑋𝑆𝑡 ∗ 𝜇𝑤 − 𝜇𝑠
Donde:
µs: Viscosidad del vapor seco saturado, cp
µws: Viscosidad del vapor húmedo y saturado, cp
Xst: Calidad del vapor en fracción
Efecto del gas en solución sobre la viscosidad del petróleo
Otras Unidades de Viscosidad

o Segundos Saybolt Universal, t:
𝐵
𝜈 = 𝐴𝑡 −
𝑡
A B Rango de t
0.226 195 32-100
0.22 135 >100
o Segundos Redwood N°1,t*:

𝐵
𝜈 = 𝐴𝑡 ∗ −
𝑡∗
A B Rango de t
0.26 179 34-100
0.247 50 >100
o Grados Engler, E:
1
1− 3
𝜈 = 𝐸(7.6) 𝐸
Densidad del petróleo

141.5 𝜌𝑜𝑠𝑐
𝛾𝑜 = 𝜌𝑜 =
𝐴𝑃𝐼 + 131.5 𝑇 − 68
1+
1885
Donde:
𝛾𝑜 : Gravedad especifica del petróleo adimensional
𝜌𝑜 : Densidad a la temperatura deseada pero en T(°F), gr/cc
𝜌𝑜𝑠𝑐 : Densidad a condiciones normales, gr/cc
Densidad del Agua:

1
𝜌𝑤 =
0.01602 + 0.000023 ∗ 𝐺
𝐺 = −6.6 + 0.0325 ∗ 𝑇 + 0.000657 ∗ 𝑇 2
𝜌𝑤 : Densidad del agua a T(°F), lb/pie3
Calor Especifico:
Hidrocarburos líquidos Hidrocarburos gaseosos

0.388 + 0.00045 ∗ 𝑇 𝑐𝑔 = 4 + 1.3 ∗ 𝑛 + 0.012 ∗ 𝑛 ∗ 𝑇
𝑐𝑜 =
𝛾𝑜
n: numero de átomos carbono/mol
𝛾𝑜 : Gravedad especifica del petróleo adimensional 𝐵𝑇𝑈
𝑐𝑔 : Calor específico del gas T(°K),
𝐿𝑏 ∗°𝐹
𝐵𝑇𝑈
𝑐𝑜 : Calor específico del petróleo T(°F),
𝐿𝑏 ∗°𝐹
Agua saturada Rocas
𝑐𝑤 = 1.0504 − 6.04 ∗ 10−4 ∗ 𝑇 + 1.79 ∗ 10−6 ∗ 𝑇 2 𝑐𝑟 = 0.18 + 0.00006 ∗ 𝑇
𝐵𝑇𝑈
𝑐𝑤 : Calor específico del agua T(°F), 𝑐𝑟 : Calor específico de la roca T(°F),
𝐵𝑇𝑈
𝐿𝑏 ∗°𝐹 𝐿𝑏 ∗°𝐹
Capacidad calorífica de las rocas
𝑀 = 1 − ∅ 𝜌𝑟 ∗ 𝑐𝑟 + ∅ ∗ 𝜌𝑜 𝑆𝑜 𝑐𝑜 + 𝜌𝑤 𝑆𝑤 𝑐𝑤 + 𝜌𝑔 𝑆𝑔 𝑐𝑔
Donde:
𝐵𝑇𝑈
𝑀: Capacidad calorífica de la roca
𝑝𝑖𝑒 3 ∗°𝐹
𝜌𝑛 : Densidad a la temperatura deseada pero en T(°F), Lb/pie3
Conductividad Térmica de líquidos y gases
Agua Líquidos
𝑇
𝐾𝑤 = 0.004 ∗ 𝑇 + 0.228 𝐾𝑙 = 0.5778 ∗ 0.0984 + 0.109 1 −
𝑇𝑏
𝐵𝑇𝑈 𝐵𝑇𝑈
𝐾𝑤 : Conductividad térmica del agua T (°K), 𝐾𝑙 : Conductividad térmica del liquido T (°K),
𝐿𝑏 ∗°𝐹∗𝑝𝑖𝑒 𝐿𝑏 ∗°𝐹∗𝑝𝑖𝑒
Tb: Temperatura de ebullición en °K
Fracciones de Petróleo y mezclas de Gases

hidrocarburos en general o Ecuación Gambill
2.48
o 𝐾𝑔 = 𝜇 ∗ 𝑐𝑝 +
Ecuación Cragoe 𝑀
𝐵𝑇𝑈
0.0677 ∗ 1 − 0.0003 ∗ 𝑇 − 32 𝐾𝑔 : Conductividad térmica del gas,
𝐾𝑜 = 𝐿𝑏 ∗°𝐹∗𝑝𝑖𝑒
𝛾𝑜 𝐿𝑏
𝜇𝑔 : Viscosidad del gas
𝑝𝑖𝑒 ∗𝑕𝑟
𝐵𝑇𝑈 𝐿𝑏
𝐾𝑜 : Conductividad térmica del petróleo T (°F) 𝑀: Peso Molecular del gas
𝐿𝑏 ∗°𝐹∗𝑝𝑖𝑒 𝐿𝑏 .𝑚𝑜𝑙
𝐵𝑇𝑈
𝑐𝑝 : Calor especifico del gas,
𝐿𝑏 ∗°𝐹
Vapor a altas Temperaturas
𝐾𝑣 = 0.5778 ∗ 10−4 ∗ 176 + 0.587 ∗ 𝑇 + 1.04 ∗ 10−3 ∗ 𝑇 2 − 4.51 ∗ 10−7 ∗ 𝑇 3

𝐵𝑇𝑈
𝐾𝑣 : Conductividad térmica del vapor T (°F),
𝐿𝑏 ∗°𝐹∗𝑝𝑖𝑒
Conductividad Térmica de las Rocas
Considerando el efecto de la densidad, Considerando el efecto de la porosidad,

saturación de fluidos y Temperatura saturación de fluidos y Temperatura
o Ecuación de Tikhomirov o Ecuación de Tikhomirov
6.36 6.36 ∗ 𝑒 0.6∗ 2.65∗ 1−∅ +𝑆𝑙

𝐾𝑕 = ∗ 𝑒 0.6∗ 𝜌𝑟 +𝑆𝑤
𝐾𝑕 =
𝑇 0.55 0.556 ∗ 𝑇 + 255.3 0.55
𝐵𝑇𝑈
𝐾𝑕 : Conductividad térmica de la roca T(°F) 𝐾𝑕 : Conductividad térmica de la roca T(°F)
𝐵𝑇𝑈
𝐻𝑟 ∗°𝐹∗𝑝𝑖𝑒
𝐻𝑟 ∗°𝐹∗𝑝𝑖𝑒
𝜌𝑟 : Densidad de la matriz o de la roca, Lb/pie3 𝜌𝑟 : Densidad de la matriz o de la roca, Lb/pie3

𝑆𝑤 : Saturación de Agua 𝑆𝑙 : Saturación de Líquido
Difusividad Térmica
𝐾𝑕
𝛼=
𝑀
Donde:
𝑃𝑖𝑒 2
𝛼: difusividad Térmica
𝑕𝑟
𝐵𝑇𝑈
𝑀: Capacidad calorífica de la roca
𝑝𝑖𝑒 3 ∗°𝐹
𝐵𝑇𝑈
𝐾𝑕 : Conductividad Térmica de la roca,
𝐻𝑟∗°𝐹∗𝑝𝑖𝑒
Saturación Residual remanente de Petróleo en la zona de vapor

1
𝑆𝑜𝑟𝑠𝑡 = ∗ 54 + 0.01004787 ∗ 𝐿𝑛𝜇𝑜 − 0.14881
𝑇
T: Temperatura en °F
EJERCICIOS RESUELTOS
2.1. La viscosidad de un petróleo es de 1700 cps a 60 °F y de 180 cps a 200 °F. Estimar la viscosidad a 400 °F y a
120 °F. La gravedad API del petróleo es de 25 °API a 60 °F. Utilice las ecuaciones de Andrade y de Walter y
la correlación generalizada de Viscosidad – Temperatura.
a) Ecuación de Andrade:
b
  ae t (1) 
 ln   ln(a)  b (2)
t
Las constantes a y b son obtenidas mediante un sistema de ecuaciones:
b b
ln(  )   ln(  )  
 b  5504,3758
T (°R) μ (cps) T T
520 1700 b b
ln(1700)   ln(180)  
 b  5504,3758
660 180 520 660
b 5504.3758
ln( a)  ln(  )  
 ln( a)  ln(1700) 
T 520
a= 0,0430
Al calcular viscosidades con datos obtenidos en ecuación 1:

5504.3758
  (0,0430)e 860
 25,891cps
5504.3758
  (0,0430)e 580
 568,925cps
b) Ecuación de Walter:
     (1  0,8)
10 n
 0,8
Donde :
 : Vis cos idad Cinemática (cSt )
 : Vis cos idad Absoluta (cps )
log( T )
T1
n: log[log( 1  0,8)]  log[log(  2  0,8)]
T2
log( )
T1
Se determina la densidad del petróleo

𝑔𝑟𝑎𝑚𝑜
𝜌𝑜𝑠𝑐 = 0.9042 ∗ 1 = 0.9042 𝑔𝑟/𝑐𝑐
𝑐𝑐
0.9042
𝜌𝑜@60 = = 0.9042
60 − 68
1+
1885
Se aplica el método de Walter
𝜇@60 1700
𝜈1 = = = 1880 𝑐𝑠𝑡
𝜌𝑜@60 0.9042
𝜇@200 180
𝜈2 = = = 213 𝑐𝑠𝑡
𝜌𝑜@200 0.8450
580
𝐿𝑜𝑔
𝑛1 = 520 𝐿𝑜𝑔 Log 1880 + 0.8 − 𝐿𝑜𝑔 Log 213 + 0.8 = 0.06768
660
𝐿𝑜𝑔
520
10 −0.06768
𝜈@120 = 1880 + 0.8 − 0.8 = 632.8436 𝑐𝑠𝑡
𝜇@120 = 𝜈@120 ∗ 𝜌𝑜@120 = 632.8436 ∗ 0.8799 = 556.9 𝑐𝑝𝑠
EC. DE WHALTER
API (°) 25 γo 0,9042 ρ (lb/pie3) 0,9042
Temperatura Temperatura Viscosidad Densidad Viscosidad Viscosidad Viscosidad

n
(°F) (°R) (Cps) (lb/pie3) (CSt) (CSt) (Cps)
60 520 1700 0,9042 1880 1700
120 580 0,8799 0,06768066 632,898562 556,9
200 660 180 0,8450 213 180
400 860 0,7688 0,31181871 38,7390011 29,8
2.2. ¿Cuál será la viscosidad de una mezcla de 30% en peso de bitumen de 1.000.000 cps y un solvente de 0.6 cps?
o Método de Cragoe
2995.73
𝐿1 = = 139.88
𝐿𝑛 1000000 + 7.6009
2995.73
𝐿2 = = 422.52
𝐿𝑛 0.6 + 7.6009
𝐿𝑚 = 𝑓1 𝐿1 + 𝑓2 𝐿2 = 0.3 ∗ 139.88 + 0.7 ∗ 422.52 = 337.73
2995.73
𝐿𝑛 𝜇𝑚 = − 7.6009 = 3.558 𝑐𝑝𝑠
𝐿𝑚
2.3. Obtenga una expresión de difusividad térmica de rocas saturadas, en función de temperatura, saturación de
agua y densidad de la roca seca. Considere una roca a 100 °F saturada 50 % con agua y 50 % con aire, la
porosidad 20 % y densidad de matriz de 2,65 gr/cc. Calcule los valores de α, K y M a temperatura de 100,
200, 300, 400 y 500 °F. Estime conveniente cualquier otro dato necesario. Tenga presente que al aumentar la
temperatura, la saturación de agua aumenta debido a la expansión térmica.
k 
 Conductividad ( BTU / h  pie   F )
 ( pie 2 / hr )  h
M
 Capacidad Calorífica ( BTU / pie 3  F )
Se empleará ecuación de Tikhomirov para hallar conductividad térmica, debido a que ofrece la expresión de
difusividad considerando Temp, Sw y ρr.
26,31 0, 6(  r  S w )
Kh  e
T 0,55
M  (1   )  r C r   (S w  w C r )
Nota: se asumirán saturaciones de agua de variando en 0,05 por cada 100 °F.
Temperatura (°F)
Temperatura (K)
Cw (BTU/PC-F)
Cr (BTU/PC-F)
(milicalorias/s-
M (BTU/pc-F)
(BTU/h-pie-F)
Conductividad
Conductividad
saturación
ρw (lb/PC)
α (pie2/hr)
cm-k)
100 310,9278 0,5 7,4126703 1,7920 0,186 1,0078 62,1347 30,8785 0,0580
200 366,4833 0,55 6,9780878 1,6870 0,192 1,0010 60,1607 37,4542 0,0450
300 422,0389 0,60 6,6535164 1,6085 0,198 1,0300 57,2986 38,0076 0,0423
400 477,5944 0,65 6,4053273 1,5485 0,204 1,0948 53,8069 38,7797 0,0399
500 533,15 0,70 6,2127743 1,5019 0,21 1,1954 49,9508 39,7345 0,0378
2.4. Una roca contiene una saturación de 0.5 de petróleo y 0.5 de agua, la gravedad °API del crudo es 25. La
viscosidad del petróleo a la temperatura original del yacimiento es 18000 cp y la porosidad de la roca es de
25%. Calcular la capacidad calorífica, la conductividad térmica, la difusividad térmica, el calor específico y
las densidades de los fluidos, a 300°F. suponer que la densidad de la roca ρr, es igual a 165 Lb/pie3, y l calor
especifico de la roca es 0.23 BTU/lb-°F
o Calculo de la Capacidad Calorífica M

141.5
𝛾𝑜 = = 0.90415
25 + 131.5
0.388 + 0.00045 ∗ 300 𝐵𝑇𝑈
𝑐𝑜 = = 0.55002
0.90415 𝐿𝑏 − °𝐹
𝐵𝑇𝑈
𝑐𝑤 = 1.0504 − 6.04 ∗ 10−4 ∗ 300 + 1.79 ∗ 10−6 ∗ 3002 = 1.03
𝐿𝑏 − °𝐹
𝑙𝑖𝑏𝑟𝑎 56.419
𝜌𝑜𝑠𝑐 = 0.90415 ∗ 62.4 = 56.419𝐿𝑏/𝑝𝑖𝑒 3 𝜌𝑜@300 = = 50.24 𝐿𝑏/𝑝𝑖𝑒 3
𝑝𝑖𝑒3 300 − 68
1+
1885
1
𝜌𝑤@300 = = 57.30𝐿𝑏/𝑝𝑖𝑒3
0.01602 + 0.000023 ∗ 62.28 𝐺 = −6.6 + 0.0325 ∗ 300 + 0.000657 ∗ 3002 = 62.28
𝐿𝑏
𝑀 = 1 − 0.25 165 ∗ 0.23 + 0.25 ∗ 50.24 ∗ 0.5 ∗ 0.55 + 57.30 ∗ 0.5 ∗ 1.03 = 39.29
𝑝𝑖𝑒 3− °𝐹
o Calculo de la Conductividad Térmica

- Ecuación de Tikhomirov
6.36 ∗ 𝑒 0.6∗ 2.65∗ 1−0.25 +1 𝐵𝑇𝑈

𝐾𝑕 = = 1.37
0.556 ∗ 300 + 255.3 0.55 𝐻𝑟 ∗ °𝐹 ∗ 𝑝𝑖𝑒
o Calculo de la Difusividad Térmica
1.37 𝑝𝑖𝑒 2
𝛼= = 0.035
39.29 𝑕𝑟
o Calculo de la Saturación Residual de Petróleo
1
𝑆𝑜𝑟𝑠𝑡 = ∗ 54 + 0.01004787 ∗ 𝐿𝑛18000 − 0.14881 = 0.13 (13%)
300
CAPITULO III. PROPIEDADES TÉRMICAS DEL AGUA Y DEL VAPOR
Temperatura de Saturación
𝑇𝑠 = 115.1 ∗ 𝑃𝑠0.225
Ts: Temperatura de saturación en °F
Calor Sensible o Calor Latente Calor Total o Entalpia

Entalpia del agua del Vapor Seco y
Saturada 𝐿𝑣 = 1318 ∗ 𝑃𝑠 −0.08774 Saturado
𝐻𝑠 = 𝐻𝑤 + 𝐿𝑣
𝐻𝑤 = 𝐶𝑤 ∗ 𝑇𝑠 − 32 𝐵𝑇𝑈
LV: Calor latente
𝐿𝑏 𝐵𝑇𝑈
𝐵𝑇𝑈 Ps: Presión de Saturación en Lpca HS: Entalpia del Vapor sat.
𝐿𝑏
Hw: Entalpia del agua sat.
𝐿𝑏
Ts: Temperatura de saturación °F
Entalpia de la Mezcla Entalpia Disponible Entalpia del Vapor

Sobrecalentado
𝐻𝑤𝑠 = 𝐻𝑤 + 𝑋 ∗ 𝐿𝑣 𝐻𝑟 = 𝐻𝑠 − 𝐶𝑤 ∗ 𝑇𝑟 − 32 𝐻𝑠𝑕 = 𝐻𝑠 + 𝐶𝑠 ∗ 𝑇𝑠𝑕 − 𝑇𝑠
HWS: Entalpia de la mezcla HWS: Entalpia de la mezcla
𝐿𝑏 𝐿𝑏
Tr: Temperatura del Yacimiento Hsh: Entalpia del Vapor
X: Calidad de la Mezcla 𝐵𝑇𝑈
en °F Sobrecalentado
𝐿𝑏
Tsh: Temperatura de
Sobrecalentamiento en °F
Volumen Especifico
𝜈𝑤𝑠 = 𝜈𝑆 ∗ 𝑋 + 1 − 𝑋 ∗ 𝜈𝑤
𝜈𝑤 = 0.01602 + 0.000023 ∗ −6.6 + 3.74075 ∗ 𝑃𝑠0.225 + 8.70394257 ∗ 𝑃𝑠0.45

𝜈𝑠 = 363.9 ∗ 𝑃𝑠 −0.9588
𝑝𝑖𝑒 3
𝜈𝑤𝑠 : Volumen especifico de la muestra
𝐿𝑏
𝑝𝑖𝑒 3
𝜈𝑤 : Volumen especifico del agua Saturada
𝐿𝑏
𝑝𝑖𝑒 3
𝜈𝑠 : Volumen especifico del Vapor Saturado
𝐿𝑏
Densidad del Vapor Calidad del Vapor
𝜌𝑠 = 0.0000440189 ∗ 𝑃𝑠0.9588 % 𝐶𝑙𝑜𝑟𝑢𝑟𝑜 𝑎𝑔𝑢𝑎 𝑑𝑒 𝑎𝑙𝑖𝑚𝑒𝑛𝑡𝑎𝑐𝑖𝑜𝑛

𝑋 = 100 −
% 𝐶𝑙𝑜𝑟𝑢𝑟𝑜 𝑝𝑎𝑟𝑡𝑒 𝐿𝑖𝑞𝑢𝑖𝑑𝑎
𝜌𝑠 : Densidad del Vapor, gr/cc
3.1. Utilizando las formulas dadas, calcular temperatura de saturación, entalpía del agua saturada, calor latente;
volumen específico y calor total para vapor húmedo con 85% de calidad y a presiones de 300, 1200 y 2500
lpca, respectivamente. Comparar las respuestas con los valores tabulados en las tablas y calcular el porcentaje
de error.
Ecuaciones de Farouq Ali
o Cálculo de la temperatura de saturación
𝑇𝑠 = 115,1 × 𝑃𝑠0,225
- A 300 Lpca:
0,225
𝑇𝑠 = 115,1 × 300 𝑙𝑝𝑐𝑎
𝑇𝑠 = 415,363 ℉
- Se repite el procedimiento anterior para las presiones de 1200 y 2500 lpca y se obtienen los
valores reflejados en la siguiente tabla:
Presión (Lpca) Temperatura de Saturación (ºF)
300 415,363
1200 567,402
1500 669,286
o Cálculo de la entalpía del agua saturada
𝑕𝑤 = 𝐶𝑤 𝑇𝑠 − 32
- Determinación del Calor Específico para agua saturada a 300lpca
𝐶𝑤 = 1,0504 − 6,05 × 10−4 . 𝑇 + 1,79 × 10−6 . 𝑇 2
𝐶𝑤 = 1,0504 − 6,05 × 10−4 415,363℉ + 1,79 × 10−6 415,363℉ 2
𝐶𝑤 = 1,108 𝐵𝑇𝑈 𝑙𝑏. ℉
Entonces se procede a determinar la entalpía del agua saturada a 300 Lpca:

𝑕𝑤 = 1,108 𝐵𝑇𝑈 𝑙𝑏. ℉ 415,363℉ − 32
𝑕𝑤 = 424,766 𝐵𝑇𝑈 𝑙𝑏𝑠
- Se repite el mismo procedimiento anterior para el resto de las presiones, obteniéndose los valores
reflejados en la siguiente tabla:
Presión (Lpca) Calor Específico (BTU/lb.ºF) Entalpía del agua saturada

(BTU/lbs)
300 1,108 424,766
1200 1,283 686,921
2500 1,447 922,153
o Cálculo de Calor Latente

- A 300 Lpca:
𝐿𝑣 = 1318 × 300 𝑙𝑝𝑐𝑎 −0,08774
𝐿𝑣 = 799,048 𝐵𝑇𝑈 𝑙𝑏𝑠
- Se aplica la misma ecuación anterior para el resto de las presiones y se obtiene:

Presión (Lpca) Calor Latente (BTU/lbs)
300 799,048
1200 707,535
2500 663,407
o Cálculo de volumen específico del vapor húmedo

- Entonces para 300 Lpca se tiene:
𝑉𝑠 = 363,9 × 300 𝑙𝑝𝑐𝑎 −0,9588
𝑉𝑠 = 1,534 𝑝𝑖𝑒 3 𝑙𝑏𝑠
𝑉𝑤 = 0,01602 + 0,000023 −6,6 + 3,74075 300 0,225 + 8,70394257 300 0,45
𝑉𝑤 = 0,019 𝑝𝑖𝑒 3 𝑙𝑏𝑠

𝑉𝑤𝑠 = 1,534 𝑝𝑖𝑒 3 𝑙𝑏𝑠 . 0,85 + 1 − 0,85 . 0,019 𝑝𝑖𝑒 3 𝑙𝑏𝑠
𝑉𝑤𝑠 = 1,307 𝑝𝑖𝑒 3 𝑙𝑏𝑠
- Se repiten los pasos anteriores para las presiones de 1200 y 2500 lpca y se obtiene lo siguiente:
Presión (Lpca) Vs (pie3/lbs) Vw (pie3/lbs) Vws (pie3/lbs)

300 1,534 0,019 1,307
1200 0,406 0,021 0,348
2500 0,201 0,023 0,174
o Cálculo de calor total

- A 300 Lpca:
𝑕𝑠 = 424,766 + 799,048
𝑕𝑠 = 1223,814 𝐵𝑇𝑈 𝑙𝑏𝑠
Se repite el procedimiento para el resto de las presiones y se obtiene:
Presión (Lpca) Calor Total (BTU/lbs)

300 1223,814
1200 1394,456
2500 1585,56
o Cálculo del porcentaje de error
Partiendo de la siguiente tabla se interpolan las presiones de 300, 1200 y 2500 lpca para obtener las
propiedades térmicas del agua a dichas presiones:
Hs Lv
Presión (lpca) Ts (°F) Hw (BTU/lb) Vs (pc/lb)
(BTU/lb) (BTU/lb)
276,75 410 385,83 1202,1 816,3 1,67
300 417 393,76 1202,82 809,05 1,5468
308,83 420 396,77 1203,1 806,3 1,5
1133,1 560 562,2 1186,4 624,2 0,3868
1200 567 581,55 1179,80 598,25 0,3396
1325,8 580 588,9 1177,3 588,4 0,3217
2365,4 660 714,2 1104,4 390,2 0,1442
2500 668 745,44 1077,44 332,00 0,1205
2708,1 680 757,3 1067,2 309,9 0,1115
Valores Interpolados
A continuación se presenta la tabla con todos los valores interpolados correspondientes a las
presiones de 300, 1200 y 2500 lpca:
TABLA DE VAPORES
Hs Lv
(BTU/lb) (BTU/lb)
300 417 394 1203 809 2
1200 566,943 581,551 1179,805 598,254 0,340
2500 667,855 745,437 1077,439 332,003 0,121
- Porcentaje de error de la temperatura de saturación a 300 lpca
𝑉𝑎𝑙𝑜𝑟 𝑡𝑎𝑏𝑢𝑙𝑎𝑑𝑜 − 𝑉𝑎𝑙𝑜𝑟 𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑑𝑜

%𝐸𝑟𝑟𝑜𝑟 = × 100
𝑉𝑎𝑙𝑜𝑟 𝑡𝑎𝑏𝑢𝑙𝑎𝑑𝑜
417 − 415,363
%𝐸𝑟𝑟𝑜𝑟 = × 100
417
%𝐸𝑟𝑟𝑜𝑟 = 0,452%
Se aplica la misma ecuación anterior para el resto de las propiedades térmicas a las diferentes
presiones y se obtiene la siguiente tabla:
ERRORES
Presión
Ts (°F) Hw (BTU/lb) Hs (BTU/lb) Lv (BTU/lb) Vs (pc/lb)
(lpca)
300 0,45% 32% 10% 1% 0,81%
1200 0,08% 50% 34% 18% 20%
2500 0,21% 58% 71% 100% 67%
Ecuaciones de Ejiogu y Fiori
o Cálculo de la Temperatura de Saturación
𝑇𝑠 = 116,79 × 𝑃0,2229
- A 1200 Lpca:
0,2229
𝑇𝑠 = 116,79 × 1200 𝑙𝑝𝑐𝑎 = 567,225 ℉
Se aplica la misma ecuación para la presión de 2500 lpca, pero para 300 lpca no aplica ya que el
rango de presiones para esta ecuación va desde 500 a 2500 lpca, y se obtiene:
Presión (lpca) Ts (°F)

300 (-)
1200 567,225
2500 668,046
o Cálculo de la entalpía del agua saturada
- Para presiones de 500 a 1500 lpca, se tiene la siguiente ecuación:
𝐻𝑤 = 77,036 × 𝑃0,28302
Entonces para 1200 lpca, se obtiene:
0,28302
𝐻𝑤 = 77,036 × 1200 𝑙𝑝𝑐𝑎 = 573,012 𝐵𝑇𝑈 𝑙𝑏𝑠
- Para presiones de 1500 a 2500 lpca, se tiene la siguiente ecuación:

𝐻𝑤 = 0,12038 × 𝑃 + 430,984
Entonces para 2500 lpca, se obtiene:
𝐻𝑤 = 0,12038 × 2500 𝑙𝑝𝑐𝑎 + 430,984 = 731,934 𝐵𝑇𝑈 𝑙𝑏𝑠
Presión (lpca) Hw (BTU/lb)

300 (-)
1200 573,012
2500 731,934
o Cálculo de calor total

1,73808
𝐻𝑠 = 1204,8 − 0,000197697 𝑃 − 453,23
- Para 1200 lpca:
𝐻𝑠 = 1204,8 − 0,000197697 1200 𝑙𝑝𝑐𝑎 − 453,23 1,73808
𝐻𝑠 = 1185,309 𝐵𝑇𝑈 𝑙𝑏𝑚
Se realiza el mismo procedimiento para la presión de 2500 lpca y tomando en cuenta que 300 lpca no se realiza por no
encontrarse dentro del rango de presiones para las cuales trabaja la ecuación, entonces se tiene:
Presión (lpca) Hs (BTU/lb)

300 (-)
1200 1185,309
2500 1092,365
o Cálculo de Calor Latente
- Para presiones de 500 a 2500 lpca:
𝐿𝑣 = 𝐻𝑠 − 𝐻𝑤
- A 1200 lpca:
𝐿𝑣 = 1185,309 − 573,012 = 612,297 𝐵𝑇𝑈 𝑙𝑏
Se aplica la misma ecuación anterior para la presión de 2500 lpca y se obtiene:
Presión (lpca) Lv (BTU/lb)

300 (-)
1200 612,297
2500 360,431
o Cálculo de volumen específico del vapor seco
490,386
𝑉𝑠 = − 0,04703
𝑃
Entonces para 1200 lpca, se tiene:
490,386
𝑉𝑠 = − 0,04703 = 0,362 𝑝𝑖𝑒 3 𝑙𝑏
1200 𝑙𝑝𝑐𝑎
551,74
𝑉𝑠 = − 0,0887
𝑃
Entonces para 2500 lpca, se tiene:
551,74
𝑉𝑠 = − 0,0887 = 0,132 𝑝𝑖𝑒 3 𝑙𝑏
2500 𝑙𝑝𝑐𝑎
Presión (lpca) Vs (pc/lb)

300 (-)
1200 0,362
2500 0,132
o Cálculo del porcentaje de error
Partiendo de la siguiente tabla de valores interpolados se puede determinar el porcentaje de error:
TABLA DE VAPORES
Hs Lv
(BTU/lb) (BTU/lb)
300 417 394 1203 809 2
1200 566,943 581,551 1179,805 598,254 0,340
2500 667,855 745,437 1077,439 332,003 0,121
El porcentaje de error viene dado por la siguiente ecuación:
𝑉𝑎𝑙𝑜𝑟 𝑡𝑎𝑏𝑢𝑙𝑎𝑑𝑜 − 𝑉𝑎𝑙𝑜𝑟 𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑑𝑜

%𝐸𝑟𝑟𝑜𝑟 = × 100
𝑉𝑎𝑙𝑜𝑟 𝑡𝑎𝑏𝑢𝑙𝑎𝑑𝑜
- Para la temperatura de saturación a 1200 lpca:
566,943 − 567,225
%𝐸𝑟𝑟𝑜𝑟 = × 100 = 0,050%
566,943
Se repite el mismo procedimiento anterior para el resto de las propiedades térmicas del agua a las diferentes presiones,
obteniéndose los siguientes valores:
ERRORES
Presión
Ts (°F) Hw (BTU/lb) Hs (BTU/lb) Lv (BTU/lb) Vs (pc/lb)
(lpca)
300 (-) (-) (-) (-) (-)
1200 0,050% 1,468% 0,467% 2% 6%
2500 0,029% 1,811% 1,385% 9% 10%
3.3. Se está inyectando vapor de 80% de calidad a una presión de 750 lpca en cierto yacimiento, a través de la
tubería de revestimiento. Si la presión del vapor en el fondo del pozo es de 600 lpca y las pérdidas se estiman
en un 46% del calor inyectado, calcular la calidad del vapor frente a la arena productora.
 Condiciones de Superficie
 Calor Sensible
𝑕𝑤 = 91 ∗ 𝑃0,2574 = 91 ∗ 750 𝑙𝑝𝑐𝑎 0,2574
𝑕𝑤 = 500,129 𝐵𝑇𝑈 𝑙𝑏
 Calor Latente
𝐿𝑣 = 1318 ∗ 𝑃 −0,08774 = 1318 ∗ 750 𝑙𝑝𝑐𝑎 −0,08774
𝐿𝑣 = 737,323 𝐵𝑇𝑈 𝑙𝑏
 Entalpia del vapor Húmedo
𝑕𝑤𝑠 = 𝑕𝑤 + ƒ𝑠1 ∗ 𝐿𝑣 = 500,129 𝐵𝑇𝑈 𝑙𝑏 + 0,8 ∗ 737,323 𝐵𝑇𝑈 𝑙𝑏

𝑕𝑤𝑠 = 1089,987 𝐵𝑇𝑈 𝑙𝑏
 Condiciones de Fondo
 Calor Sensible
𝑕𝑤 = 91 ∗ 𝑃0,2574 = 91 ∗ 600 𝑙𝑝𝑐𝑎 0,2574
𝑕𝑤 = 472,213 𝐵𝑇𝑈 𝑙𝑏
 Calor Latente
𝐿𝑣 = 1318 ∗ 𝑃 −0,08774 = 1318 ∗ 600 𝑙𝑝𝑐𝑎 −0,08774
𝐿𝑣 = 750,901 𝐵𝑇𝑈 𝑙𝑏
Al tener una pérdida de calor de 46%, al fondo solo llegara un 54% del calor inyectado por lo que la entalpia del vapor
húmedo será el 54% de la entalpia en condiciones de superficie:
𝑕𝑤𝑠 = 0,54 1089,987 𝐵𝑇𝑈 𝑙𝑏 → 𝑕𝑤𝑠 = 588,593 𝐵𝑇𝑈 𝑙𝑏
Partiendo entonces de la ecuación general de la entalpia del vapor húmedo se procede a calcular la calidad del vapor
frente a la arena productora con las condiciones de fondo ya calculadas:
𝑕𝑤𝑠 − 𝑕𝑤
𝑕𝑤𝑠 = 𝑕𝑤 + ƒ𝑠2 ∗ 𝐿𝑣 → ƒ𝑠2 =
𝐿𝑣
588,593 𝐵𝑇𝑈 𝑙𝑏 − 472,213 𝐵𝑇𝑈 𝑙𝑏

ƒ𝑠2 = = 0,155 (15,5%)
750,901 𝐵𝑇𝑈 𝑙𝑏
3.4. En un proyecto térmico se desea inyectar un fluido caliente a una temperatura no mayor de 600°F (basado en
condiciones de resistencia del equipo del pozo) en una formación a una temperatura de 110 °F y se requiere
suministra a la formación 5x106 BTU/h. las posibles escogencias son: aire caliente ( ce = 0.249 BTU/Lb-°F),
agua caliente (ce=1.195 BTU/Lb-°F), vapor saturado de 80% de calidad, vapor seco y saturado y vapor súper
calentado con 50°F de supercalentamiento. ¿Qué cantidad de fluido en Lb/h tendrá que ser inyectado en cada
caso?. Obtener las cantidades relativas con referencia al vapor húmedo.
 Aire Caliente (Ce= 0,249 BTU/Lb-°F)
𝑇𝑠 = 600°𝐹
𝐻𝑎𝑐 = 𝐶𝑎𝑐 ∗ 𝑇𝑠 − 32
𝐵𝑇𝑈
𝐻𝑎𝑐 = 0.249 ∗ 600 − 32 = 141.432
𝐿𝑏
𝐵𝑇𝑈
5𝑥106
𝑚= 𝑕𝑟 = 35352.67 𝐿𝑏
𝐵𝑇𝑈 𝑕𝑟
141.432
𝑙𝑏
 Agua Caliente (Ce= 1.195 0,249 BTU/Lb-°F)
𝑇𝑠 = 600°𝐹
𝐻𝑤 = 𝑤 ∗ 𝑇𝑠 − 32
𝐵𝑇𝑈
𝐻𝑤 = 1.195 ∗ 600 − 32 = 678.76
𝐿𝑏
𝐵𝑇𝑈
5𝑥106
𝑚= 𝑕𝑟 = 7366.3739 𝐿𝑏
678.76
𝑙𝑏
 Vapor saturado de 80%
Con la tabla 3.3 del libro de recuperación térmica de petróleo
Hw Lv
T (°F) BTU/Lb (BTU/Lb)
600 617 548.5
𝐻𝑤𝑠 = 𝐻𝑤 + 𝑋𝐿𝑣
𝐻𝑤𝑠 = 617 + 0.80 ∗ 548.5
𝐵𝑇𝑈
𝐻𝑤𝑠 = 1055.8
𝐿𝑏
𝐵𝑇𝑈
5𝑥10 6 𝑕 𝑟 𝐿𝑏
𝑚= 𝐵𝑇𝑈 = 4735.74
1055.8 𝑙𝑏 𝑕𝑟
 Vapor seco y saturado
Con la tabla 3.3 del libro de recuperación térmica de petróleo
Hw Lv
T (°F) (BTU/Lb) (BTU/Lb)
600 617 548.5
𝐻𝑠 = 𝐻𝑤 + 𝐿𝑣
𝐻𝑠 = 617 + 548.5
𝐵𝑇𝑈
𝐻𝑤𝑠 = 1165.5
𝐿𝑏
𝐵𝑇𝑈
5𝑥106
𝑚= 𝑕𝑟 = 4290 𝐿𝑏
1165.5
𝑙𝑏
3.5. Vapor a una presión de 200 lpca está siendo expandido a entalpía constante a través de una válvula, a la
presión atmosférica. El vapor resultante tiene una temperatura de 219 °F. Calcular la calidad del vapor
original. (El calor específico vapor a 219 °F es de 0,50 BTU/lb-ªF). Este proceso ilustra una técnica para
determinar la calidad del vapor.
A presión atmosférica (14,7 lpca), la Temperatura de saturación es 212 °F, lo que permite saber que el vapor
se encuentra sobrecalentado en 7 °F. Permitiendo hallar su entalpía con la siguiente ecuación:
𝐻𝑠𝑕 = 𝐻𝑠 + 𝐶𝑠 𝑇𝑠𝑕 − 𝑇𝑠
Para aplicar la formula anterior se debe buscar las propiedades del vapor a la condicion de sobrecalentamiento es decir
a 219°F
𝑇𝑠 = 115.1 ∗ 𝑃𝑠0.225
1
𝑇𝑠 0.225
𝑃𝑠 = = 17.44 𝐿𝑝𝑐𝑎
115.1
𝐿𝑣 = 1318 ∗ 17.44−0.08774 = 1025.59 𝐵𝑇𝑈/𝑙𝑏

𝐻𝑤 = 𝐶𝑤 ∗ 𝑇𝑠 − 32 = 1 ∗ 219 − 32 = 187 𝐵𝑇𝑈/𝑙𝑏
𝐻𝑠 = 187 + 1025.59 = 1212.59 𝐵𝑇𝑈/𝑙𝑏
𝐻𝑠𝑕 = 1212.59 + 0,5 219 + 212 = 1216.09 𝐵𝑇𝑈/𝑙𝑏
Ya que se trata de un proceso de entalpia constante entonces:
𝐻𝑠𝑕 = 𝐻𝑠𝑤 =1216.09 𝐵𝑇𝑈/𝑙𝑏
La ecuación de cálculo de entalpía para vapor húmedo es:
𝐻𝑠𝑤 = 𝐻𝑤 + 𝑋. 𝐿𝑣
Siendo X la calidad, valor que se desea despejar:
𝐻𝑤𝑠 − 𝐻𝑤
𝑋=
𝐿𝑣
A 200 Lpca, en tabla de vapores:
𝑇𝑠 = 115.1 ∗ 2000.225 = 379.15 °𝐹

𝐻𝑤 = 𝐶𝑤 ∗ 379.15 − 32 = 1 ∗ 379.15 − 32 = 347 𝐵𝑇𝑈/𝑙𝑏
𝐿𝑣 = 1318 ∗ 200−0.08774 = 827.986 𝐵𝑇𝑈/𝑙𝑏
1216.09 − 347
𝑋= =
827.986
Por lo tanto, se dice que el vapor inicial tiene una calidad de 107 %. Lo que indica que el vapor se encontraba
saturado o seco, cerca del sobrecalentamiento.
Capítulo IV. Perdidas de Calor durante la transmisión de fluidos calientes
SUPERFICIE
Coeficiente de radiación
𝑕𝑟 = 𝜍 ∗ 𝜀 ∗ 𝑇1 2 + 𝑇2 2 ∗ 𝑇1 + 𝑇2
𝜍 = 0.1713 𝑋10−8 𝐵𝑇𝑈/(𝐻𝑟 ∗ 𝑝𝑖𝑒 3 ∗ °𝐹)

𝐵𝑇𝑈
𝑕𝑟 : Coeficiente de radiación
𝑕𝑟∗𝑝𝑖𝑒 2∗°𝐹
T: Temperatura del acero °R
𝜀: emisibidad
Coeficiente de convección
0.25
3 𝑔 ∗ 𝛽𝑎 𝐶𝑝𝑎 ∗ 𝜇𝑎
0.53 ∗ 𝐾𝑕𝑎 ∗ 𝑑𝑒 ∗ 𝑇𝑠𝑢𝑟𝑓 − 𝑇𝑎 ∗
𝜈𝑎 2 𝐾𝑕𝑎
𝑕𝑐 =
𝑑𝑒
- Con viento
0.6
19.3 ∗ 𝐶𝑝𝑎 ∗ 𝜌𝑎 ∗ 𝑉𝑖𝑒𝑛𝑡𝑜
𝑕𝑐 = 10000 ≤ 8800 ∗ 𝑑𝑒 ∗ 𝑣𝑖𝑒𝑛𝑡𝑜 ≤ 50000
𝑑𝑒 0.4
0.805
0.0239 ∗ 𝐾𝑕𝑎 5280 ∗ 𝜌𝑎 ∗ 𝑉𝑖𝑒𝑛𝑡𝑜 ∗ 𝑑𝑒
𝑕𝑐 = ∗ 8800 ∗ 𝑑𝑒 ∗ 𝑣𝑖𝑒𝑛𝑡𝑜 ≥ 50000
𝑑𝑒 𝜇𝑎
Propiedades del Aire
𝐾𝑕𝑎 = 0.01328 + 2.471𝑥10−5 ∗ 𝑇 − 4.247𝑋10−9 ∗ 𝑇 2 𝜇𝑎 = 0.04 + 6.155𝑥10−5 ∗ 𝑇 − 1.22𝑋10−8 ∗ 𝑇 2

BTU 𝐿𝑏
𝐾𝑕𝑎: 𝜇𝑎 :
hr − pie − °F 𝑝𝑖𝑒 − 𝑕𝑟
𝑐𝑝𝑎 = 0.2382 + 1.39𝑥10−5 ∗ 𝑇 + 1.027𝑋10−8 ∗ 𝑇 2
𝐵𝑇𝑈
𝑐𝑝𝑎 :
𝐿𝑏 − °𝐹
𝜌𝑎 = 8.55865𝑥10−2 − 1.5531𝑥10−4 ∗ 𝑇 + 1.65602𝑋10−7 ∗ 𝑇 2 − 6.92225𝑥10−11 ∗ 𝑇 3
𝐿𝑏
𝜌𝑎 :
𝑝𝑖𝑒 3
𝛽𝑎 = 2.15844𝑥10−3 − 3.89367𝑥10−6 ∗ 𝑇 + 4.12773𝑋10−9 ∗ 𝑇 2 − 1.71867𝑥10−12 ∗ 𝑇 3
𝛽𝑎 : °𝐹−1
𝜇𝑎 𝑝𝑖𝑒 2
𝜈𝑎 =
𝜌𝑎 𝑕𝑟
Coeficiente de Transferencia de Calor

𝑟 𝑟 −1
𝑟𝑡𝑜 𝑟𝑡𝑜 ∗ 𝐿𝑛 𝑡𝑜 𝑟𝑡𝑜 ∗ 𝐿𝑛 𝑖𝑛𝑠 𝑟𝑡𝑜
𝑟𝑡𝑖 𝑟𝑡𝑜
𝑈𝑡𝑜 = + + +
𝑟𝑡𝑖 ∗ 𝑕𝑓 𝐾𝑕𝑠 𝐾𝑕𝑖𝑛𝑠 𝑟𝑖𝑛𝑠 ∗ 𝑕𝑟 + 𝑕𝑐
𝐵𝑇𝑈
𝑈𝑡𝑜 : Coeficiente de Transferencia de calor
𝐵𝑇𝑈 𝑟𝑖𝑛𝑠 : radio del aislante, pies

Khs: conductividad del acero
𝑕𝑟∗𝑝𝑖𝑒 ∗°𝐹
𝑟𝑡𝑜 : radio externo de la tubería, pies
𝐵𝑇𝑈
Khins: Conductividad del aislante
𝑕𝑟∗𝑝𝑖𝑒 ∗°𝐹 𝑟𝑖𝑛𝑠 = 𝑟𝑡𝑜 + ∆𝑟𝑖𝑛𝑠
Perdidas de Calor
𝑄 = 2𝜋𝑟𝑡𝑜 ∆𝐿 ∗ 𝑈𝑡𝑜 ∗ 𝑇𝑠𝑢𝑟𝑓 − 𝑇𝑎
𝐵𝑇𝑈
𝑈𝑡𝑜 : Coeficiente de Transferencia de calor
Q: Perdidas de calor (BTU / hr)
Tsurf: Temperatura en la superficie de la tubería (°F)
Coeficiente de Película
0.8 𝑐𝑝𝑓 ∗ 𝜇𝑓 0.4
𝑑𝑡𝑖 ∗ 𝐺
0.23 ∗ 𝐾𝑕𝑓 ∗ ∗
𝜇𝑓 𝑘𝑕𝑓
𝑕𝑓 =
𝑑𝑡𝑖
𝐵𝑇𝑈
𝑕𝑓 : Coeficiente de película
G: flujo másico (Lb / pie3)
Propiedades del agua en la tabla 4.6
Calidad (X)
𝑄
𝑓𝑠𝑡,𝑐𝑎𝑏𝑒𝑧𝑎𝑙 = 𝑓𝑠𝑡,𝑔𝑒𝑛 −
𝑊 ∗ 𝐿𝑣
4.1 Una tubería de 2000 pies de longitud y diámetro exterior igual a 2.25 pulg, transporta vapor a una tasa de 350
B/D (equivalentes de agua). La presión del vapor a la salida del generador es 1800 Lpca y la emisividad de la
superficie exterior de la tubería es igual a 1.0. Calcular las pérdidas de calor, considerando temperatura
ambiente igual a 0°F y velocidad del viento despreciable. Repetir para caso de tubería aislada con un aislante
de magnesio de 1 pulg de espesor y conductividad térmica de 0.04 BTU/h-pie-°f. Repetir para ambos casos
asumiendo un radio interno de la tubería, y posteriormente una Tsurf= 140.8 °F
Solución:
Tubería Desnuda
Calculo de la temperatura de saturación:
𝑇𝑠 = 115.1 ∗ 𝑃0.225
𝑇𝑠𝑢𝑟𝑓 = 𝑇𝑠 = 115.1 ∗ 18000.225 = 621.6°𝐹 = 1081.6°𝑅
Calculo de las propiedades físicas del aire a Tavg
621.6 + 0
𝑇𝑎𝑣𝑔 = = 310.8 °𝐹
2
𝐾𝑕𝑎 = 0.01328 + 2.471𝑥10−5 ∗ 𝑇 − 4.247𝑋10−9 ∗ 𝑇 2 (𝐸𝑐. 1)
BTU
𝐾𝑕𝑎 = 0.020549632
hr − pie − °F
𝜇𝑎 = 0.04 + 6.155𝑥10−5 ∗ 𝑇 − 1.22𝑋10−8 ∗ 𝑇 2 (𝐸𝑐. 2)
𝐿𝑏
𝜇𝑎 = 0.057951285
𝑝𝑖𝑒 − 𝑕𝑟
𝑐𝑝𝑎 = 0.2382 + 1.39𝑥10−5 ∗ 𝑇 + 1.027𝑋10−8 ∗ 𝑇 2 (𝐸𝑐. 3)
𝐵𝑇𝑈
𝑐𝑝𝑎 = 0.243512177
𝐿𝑏 − °𝐹
𝜌𝑎 = 8.55865𝑥10−2 − 1.5531𝑥10−4 ∗ 𝑇 + 1.65602𝑋10−7 ∗ 𝑇 2 − 6.92225𝑥10−11 ∗ 𝑇 3 (𝐸𝑐. 4)
𝐿𝑏
𝜌𝑎 = 0.051234502
𝑝𝑖𝑒 3
𝛽𝑎 = 2.15844𝑥10−3 − 3.89367𝑥10−6 ∗ 𝑇 + 4.12773𝑋10−9 ∗ 𝑇 2 − 1.71867𝑥10−12 ∗ 𝑇 3 (𝐸𝑐. 5)
𝛽𝑎 = 0.001295413 °𝐹 −1
𝜇𝑎
𝜈𝑎 = (𝐸𝑐. 6)
𝜌𝑎
𝑝𝑖𝑒 2
𝜈𝑎 = 1.131098836
𝑕𝑟
Tabla de Propiedades físicas del aire a Tavg=310.8°F
Propiedades del aire sin aislante

Kha (BTU/pie-hr-°F) 0.020549632
µa (Lb/pie-hr) 0.057951285
cpa (BTU/Lb-°F) 0.243512177
ρ (Lb/pie3) 0.051234502
βa (1/°F) 0.001295413
νa (pie2/hr) 1.131098836
Calculo de hc
0.25
0.53 ∗ 𝐾𝑕𝑎 ∗ 𝑑𝑒 ∗ 𝑇𝑠𝑢𝑟𝑓 − 𝑇𝑎 ∗
𝑕𝑐 = (𝐸𝑐. 7)
𝑑𝑒
Se sustituyen los valores sabiendo que:
g= 4.17x108 pie/h2
de=2.25 pulg= 0.1875 pies

0.25
4.17𝑥108 ∗ 0.001295 0.243512 ∗ 0.0579
0.53 ∗ 0.02055 ∗ 0.18753 ∗ 621.6 − 0 ∗
1.131092 0.02055
𝑕𝑐 =
0.1875
𝐵𝑇𝑈
𝑕𝑐 = 1.917733444
𝑕𝑟 − 𝑝𝑖𝑒 3 − °𝐹
Calculo de hr
𝑕𝑟 = 𝜍 ∗ 𝜀 ∗ 𝑇1 ∗2 + 𝑇2 ∗2 ∗ 𝑇1 ∗ + 𝑇2 ∗ (𝐸𝑐. 8)
Donde:
σ=0.1713x10-8 BTU/hr-pie2-°F4
Se sustituyen los valores:
𝑕𝑟 = 0.1713𝑥10−8 ∗ 1 ∗ 1081.62 + 4602 ∗ 1081.6∗ + 460∗
𝐵𝑇𝑈
𝑕𝑟 = 3.648108914
𝑕 − 𝑝𝑖𝑒 2 − °𝐹
Calculo de coeficiente de transferencia de calor total (U)
Si se toma en cuenta el área exterior de la tubería, se tiene:
𝑟 𝑟 −1
𝑟𝑡𝑜 𝑟𝑡𝑜 ∗ 𝐿𝑛 𝑡𝑜 𝑟𝑡𝑜 ∗ 𝐿𝑛 𝑖𝑛𝑠 𝑟𝑡𝑜
𝑈𝑡𝑜 = + + +
Considerando que no se conoce el radio interno de la tubería desnuda,

se puede despreciar la convección forzada en las paredes internas de
la tubería y además sabiendo que los coeficientes de película (hf) para
el agua caliente y para vapor son bastante altos, ofrecen muy poca
resistencia al flujo por lo que pueden despreciarse, pero en la parte 3
de este ejercicio se calculara para evaluar su efecto. A través de la
tubería existe conducción, pero ya que no se conoce el radio interno
de la tubería también se puede despreciar y por ultimo también se
desprecia la conducción a través del aislante ya que no se posee
aislante. Por lo que solo se tendrá un mecanismo de transferencia de
calor que es la convección y radiación entre la tubería y el medio
ambiente por lo que la ecuación anterior, resulta:
𝑟𝑡𝑜 −1
𝑈𝑡𝑜 =
𝑟𝑖𝑛𝑠 ∗ 𝑕𝑟 + 𝑕𝑐
Sabiendo
𝑟𝑖𝑛𝑠 = 𝑟𝑡𝑜 + ∆𝑟𝑖𝑛𝑠
Donde ∆rins es igual a 0 ya que no existe aislante por lo tanto rto=rins:

−1
1
𝑈𝑡𝑜 =
𝑕𝑟 + 𝑕𝑐
Se sustituye los valores:

−1
1
𝑈𝑡𝑜 =
3.648108914 + 1.917733444
𝐵𝑇𝑈
𝑈𝑡𝑜 = 5.565842358
𝑕 − 𝑝𝑖𝑒 2 − °𝐹
Calculo de las pérdidas de calor
𝑄 = 2𝜋𝑟𝑡𝑜 ∆𝐿 ∗ 𝑈𝑡𝑜 ∗ 𝑇𝑠𝑢𝑟𝑓 − 𝑇𝑎
𝑄 = 2𝜋 ∗ 0.1875 ∗ 2000 ∗ 5.565842358 ∗ 621.6 − 0
𝐵𝑇𝑈
𝑄 = 4075901.493 = 4.077𝑀𝑀𝐵𝑇𝑈/𝑕𝑟
𝑕𝑟
Tubería Aislada
Calculo de la temperatura en la superficie exterior del aislante

(Tsurf), la temperatura promedio para el cálculo de propiedades
del aire (Tavg) y diametro del aislante (rins)
621.6
𝑇𝑠𝑢𝑟𝑓 = = 310.8°𝐹 = 770°𝑅
2
310.8
𝑇𝑎𝑣𝑔 = = 155.4°𝐹
2
(4.25 𝑝𝑢𝑙𝑔)
𝑑𝑖𝑛𝑠 = = 0.354 𝑝𝑖𝑒
12
0.354
𝑟𝑖𝑛𝑠 = = 0.177 𝑝𝑖𝑒
2
Calculo de las propiedades físicas del aire a Tavg = 155.4 °F, aplicando las ecuaciones desde la 1 (Ec. 1) hasta
la ecuación 6 (Ec.6).
Propiedades del aire Con Aislante

Kha (BTU/pie-hr-°F) 0.017017378
µa (Lb/pie-hr) 0.049270263
cpa (BTU/Lb-°F) 0.240608076
ρ (Lb/pie3) 0.065190674
βa (1/°F) 0.001646594
νa (pie2/hr) 0.755786992
Calculo de hc con la ecuación 7, usando en este caso el diámetro externo del aislante (0.354 pie)
0.25
4.17𝑥108 ∗ 0.001646 0.240608 ∗ 0.04927
0.53 ∗ 0.01702 ∗ 0.3543 ∗ 310.8 − 0 ∗
0.7557872 0.01702
𝑕𝑐 =
0.354
𝐵𝑇𝑈
𝑕𝑐 = 1.484938145
𝑕𝑟 − 𝑝𝑖𝑒 3 − °𝐹
Calculo de hr con la ecuación 8
𝑕𝑟 = 0.1713𝑥10−8 ∗ 1 ∗ 770.82 + 4602 ∗ 770.8∗ + 460∗
𝐵𝑇𝑈
𝑕𝑟 = 1.698776884
𝑕 − 𝑝𝑖𝑒 2 − °𝐹
Si se toma en cuenta el área interior de la tubería, se tiene:
𝑟 𝑟 −1
𝑟𝑖𝑛𝑠 𝑟𝑖𝑛𝑠 ∗ 𝐿𝑛 𝑡𝑜 𝑟𝑖𝑛𝑠 ∗ 𝐿𝑛 𝑖𝑛𝑠 1
𝑈𝑡𝑖𝑛𝑠 = + + + (𝐸𝑐. 10)
𝑟𝑡𝑖 ∗ 𝑕𝑓 𝐾𝑕𝑠 𝐾𝑕𝑖𝑛𝑠 𝑕𝑟 + 𝑕𝑐
Considerando que no se conoce el radio interno de la tubería desnuda,

se puede despreciar la convección forzada en las paredes internas de la
tubería y el mecanismo de conducción a través de la tuberia. Pero ahora
si se posee un aislante como se observa en la figura entonces existe
conducción a través del aislante además existe convección y radiación
entre el aislante y el medio ambiente por lo que la ecuación 10, resulta:
𝑟𝑖𝑛𝑠 −1
𝑟𝑖𝑛𝑠 ∗ 𝐿𝑛 1
𝑟𝑡𝑜
𝑈𝑡𝑖𝑛𝑠 = +
𝐾𝑕𝑖𝑛𝑠 𝑕𝑟 + 𝑕𝑐
𝑟𝑖𝑛𝑠 −1
𝑟𝑖𝑛𝑠 ∗ 𝐿𝑛 1
𝑟𝑡𝑜
𝑈𝑡𝑖𝑛𝑠 = +
𝐾𝑕𝑖𝑛𝑠 𝑕𝑟 + 𝑕𝑐
−1
0.177
0.177 ∗ 𝐿𝑛 1
𝑈𝑡𝑖𝑛𝑠 = 0.09375 +
0.04 1.698776884 + 1.484938145
𝐵𝑇𝑈
𝑈𝑡𝑖𝑛𝑠 = 0.319452544
𝑕 − 𝑝𝑖𝑒 2 − °𝐹
𝑄 = 2𝜋𝑟𝑡𝑖 ∆𝐿 ∗ 𝑈𝑡𝑖𝑛𝑠 ∗ 𝑇𝑠 − 𝑇𝑎
𝑄 = 2𝜋 ∗ 0.1777 ∗ 2000 ∗ 0.318441842 ∗ 621.6 − 0
𝐵𝑇𝑈
𝑄 = 441922.8046 = 441.9𝑀𝐵𝑇𝑈/𝑕𝑟
𝑕𝑟
Calculo de Tsur* y U*
𝐾𝑕𝑖𝑛𝑠
𝑈∗ = 𝑟𝑖𝑛𝑠 𝐸𝑐. 11
𝑟𝑖𝑛𝑠 ∗ 𝐿𝑛
𝑟𝑡𝑜
𝑈 ∗ = 0.355081236
𝑄
𝑇𝑠𝑢𝑟𝑓 ∗ = 𝑇𝑠 −
𝐴𝑈 ∗
𝑇𝑠𝑢𝑟𝑓 ∗ = 62.371157 °𝐹
Dado que la diferencia entre Tsuf supuesto (310.8 °F) y Tsurf calculado (62.371 °F) difieren de una aproximación de
0.1°, los cálculos deben repetirse tomando en cuenta como nuevo Tsuf supuesto el calculado (62.371 °F), se deben
calcular todas las propiedades del aire nuevamente, los coeficiente de convección y radiación (hc y hr), el coeficiente
de calor transferido (U) para posteriormente calcular las pérdidas de calor (Q) y obtener un nuevo Tsurf calculado. Se
proceden las iteraciones hasta que la diferencia sea de 0.1°. A continuación se muestra las iteraciones
Tabla de Resultados:
Propiedades del Aire

Kha µa cpa hr hc U
Tavg (BTU/pie-hr- (Lb/pie- (BTU/Lb- ρ βa νa (BTU/pie2- (BTU/pie2- (BTU/pie2-
Iteraciones Tsur (°f) Tsur* (°F) (°F) °F) hr) °F) (Lb/pie3) (1/°F) (pie2/hr) hr-°F) hr-°F) hr-°F) Q (BTU/h)
1 310.800 62.371 155.400 0.01702 0.04927 0.24061 0.06519 0.00165 0.75579 1.69878 1.48494 0.31945 441922.80
2 62.371 99.349 31.186 0.01405 0.04191 0.23864 0.08090 0.00204 0.51800 0.81527 1.05129 0.29833 412701.25
3 99.349 90.338 49.675 0.01450 0.04303 0.23892 0.07827 0.00197 0.54972 0.91580 1.17237 0.30348 419822.56
4 90.338 92.257 45.169 0.01439 0.04276 0.23885 0.07890 0.00199 0.54187 0.89040 1.14694 0.30238 418305.65
5 92.257 91.837 46.129 0.01441 0.04281 0.23886 0.07877 0.00199 0.54354 0.89576 1.15254 0.30262 418638.03
6 91.837 91.928 45.918 0.01441 0.04280 0.23886 0.07880 0.00199 0.54317 0.89459 1.15132 0.30257 418565.65
7 91.928 91.908 45.964 0.01441 0.04280 0.23886 0.07879 0.00199 0.54325 0.89484 1.15159 0.30258 418581.44
Dif. Tsurf
(°F) 0.020
4.2 Una tubería de 3 pulg (diámetro exterior), transporta vapor a 600°F, si la emisividad de la superficie exterior
es igual a 0.8 y la temperatura ambiente 100°F calcular:
a) El coeficiente de transferencia de calor por radiación, hr

b) El coeficiente de transferencia de calor por convección, hc
c) La tasa de pérdidas de calor por pie de longitud, debida a convección y radiación
d) El efecto de la velocidad del viento de 30 millas/hr sobre las pérdidas de calor
Tubería Desnuda sin considerar la conductividad del acero
Calculo de las propiedades físicas del aire a T avg

600 + 100
𝑇𝑎𝑣𝑔 = = 350°𝐹
2
Tabla de Propiedades físicas del aire a Tavg=350°F
Propiedades del aire Khs =0

Kha (BTU/pie-hr-°F) 0.021408243
µa (Lb/pie-hr) 0.060048
cpa (BTU/Lb-°F) 0.244323075
ρ (Lb/pie3) 0.04854633
βa (1/°F) 0.001227614
νa (pie2/hr) 1.236921506
Calculo de hc
0.25
0.53 ∗ 𝐾𝑕𝑎 ∗ 𝑑𝑒 ∗ 𝑇𝑠𝑢𝑟𝑓 − 𝑇𝑎 ∗
𝑕𝑐 = (𝐸𝑐. 7)
𝑑𝑒
Se sustituyen los valores sabiendo que:
g= 4.17x108 pie/h2
de=3 pulg= 0.25 pies

0.25
4.17𝑥108 ∗ 0.001227 0.2443 ∗ 0.06005
0.53 ∗ 0.0214 ∗ 0.253 ∗ 600 − 100 ∗
1.23692 0.0214
𝑕𝑐 =
0.25
𝐵𝑇𝑈
𝑕𝑐 = 1.660402313
𝑕𝑟 − 𝑝𝑖𝑒 3 − °𝐹
Calculo de hr
𝑕𝑟 = 𝜍 ∗ 𝜀 ∗ 𝑇1 ∗2 + 𝑇2 ∗2 ∗ 𝑇1 ∗ + 𝑇2 ∗ (𝐸𝑐. 8)
Donde:
σ=0.1713x10-8 BTU/hr-pie2-°F4
𝑕𝑟 = 0.1713𝑥10−8 ∗ 0.8 ∗ 10602 + 5602 ∗ 1060∗ + 560∗
𝐵𝑇𝑈
𝑕𝑟 = 3.190652986
𝑕 − 𝑝𝑖𝑒 2 − °𝐹
Si se toma en cuenta el área exterior de la tubería, se tiene:
𝑟𝑡𝑜 𝑟𝑖𝑛𝑠 −1
𝑟𝑡𝑜 𝑟𝑡𝑜 ∗ 𝐿𝑛 𝑟𝑡𝑜 ∗ 𝐿𝑛 𝑟𝑡𝑜
𝑈𝑡𝑜 = + + + (𝐸𝑐. 9)
Considerando que no se tomara en cuenta la conductividad del acero ni la

convección forzada en el interior de la tubería y además sabiendo que los
coeficientes de película (hf) para el agua caliente y para vapor son bastante
altos, ofrecen muy poca resistencia al flujo por lo que pueden despreciarse.
Por último también se desprecia la conducción a través del aislante ya que
no se posee aislante. Por lo que solo se tendrá un mecanismo de
transferencia de calor que es la convección y radiación entre la tubería y el
medio ambiente por lo que la ecuación anterior, resulta:
𝑟𝑡𝑜 −1
𝑈𝑡𝑜 =
𝑟𝑖𝑛𝑠 ∗ 𝑕𝑟 + 𝑕𝑐
Sabiendo
𝑟𝑖𝑛𝑠 = 𝑟𝑡𝑜 + ∆𝑟𝑖𝑛𝑠
Donde ∆rins es igual a 0 ya que no existe aislante por lo tanto rto=rins:

−1
1
𝑈𝑡𝑜 =
𝑕𝑟 + 𝑕𝑐
Se sustituye los valores:

−1
1
𝑈𝑡𝑜 =
3.190652986 + 1.660402313
𝐵𝑇𝑈
𝑈𝑡𝑜 = 4.851055299
𝑕 − 𝑝𝑖𝑒 2 − °𝐹
𝑄 = 2𝜋𝑟𝑡𝑜 ∆𝐿 ∗ 𝑈𝑡𝑜 ∗ 𝑇𝑠𝑢𝑟𝑓 − 𝑇𝑎 𝐸𝑐. 10
𝑄
= 2𝜋 ∗ 0.125 ∗ 4.851055299 ∗ 600 − 100
∆𝐿
𝑄 𝐵𝑇𝑈
= 1905.004961
∆𝐿 𝑕𝑟−𝑝𝑖𝑒
Tubería Desnuda considerando la conductividad del acero
Se asume una temperatura en la superficie exterior del aislante (Tsurf), y se procede a calcular la temperatura
promedio para el cálculo de propiedades del aire (Tavg)
𝑇𝑠𝑢𝑟𝑓 = 550°𝐹 = 1010°𝑅
550 + 100
𝑇𝑎𝑣𝑔 = = 325°𝐹
2
Calculo de las propiedades físicas del aire a Tavg = 325 °F, aplicando las ecuaciones desde la 1 (Ec. 1) hasta
la ecuación 6 (Ec.6).
Propiedades del aire Khs ≠ 0

Kha (BTU/pie-hr-°F) 0.02141
µa (Lb/pie-hr) 0.05872
cpa (BTU/Lb-°F) 0.24380
ρ (Lb/pie3) 0.05023
βa (1/°F) 0.001227
νa (pie2/hr) 1.236921
Calculo de hc con la ecuación 7
𝐵𝑇𝑈
𝑕𝑐 = 1.63544
𝑕𝑟 − 𝑝𝑖𝑒 3 − °𝐹
Calculo de hr con la ecuación 8
𝐵𝑇𝑈
𝑕𝑟 = 2.86949
𝑕 − 𝑝𝑖𝑒 2 − °𝐹
𝑟𝑡𝑜 𝑟𝑖𝑛𝑠 −1
𝑟𝑡𝑜 𝑟𝑡𝑜 ∗ 𝐿𝑛 𝑟𝑡𝑜 ∗ 𝐿𝑛 𝑟𝑡𝑜
𝑈𝑡𝑜 = + + + (𝐸𝑐. 9)
Considerando que en este caso se tomara en cuenta la conductividad del acero pero sigue sin tomarse en cuenta la
convección forzada en el interior de la tubería ni la conducción a través del aislante ya que no se posee aislante. Por lo
que la ecuación 9, resulta:
𝑟𝑡𝑜 −1
𝑟𝑡𝑜 ∗ 𝐿𝑛 1
𝑟𝑡𝑖
𝑈𝑡𝑜 = +
𝐾𝑕𝑠 𝑕𝑟 + 𝑕𝑐
Sustituyendo los valores:
Donde:
𝐵𝑇𝑈
𝐾𝑕𝑠 = 26.2
𝑕 − 𝑝𝑖𝑒 − °𝐹
−1
0.125
0.125 ∗ 𝐿𝑛 1
𝑈𝑡𝑜 = 0.104 +
26.2 1.63544 + 2.86949
𝐵𝑇𝑈
𝑈𝑡𝑜 = 4.48719
𝑕 − 𝑝𝑖𝑒 2 − °𝐹
𝑄 = 2𝜋𝑟𝑡𝑜 ∆𝐿 ∗ 𝑈𝑡𝑜 ∗ 𝑇𝑠 − 𝑇𝑎
𝑄
= 2𝜋 ∗ 0.125 ∗ 6.083025224 ∗ 621.6 − 0
∆𝐿
𝐵𝑇𝑈
𝑄 = 1762.11720
𝑕𝑟 − 𝑝𝑖𝑒
Calculo de Tsur* y U*
𝐾𝑕𝑠
𝑈∗ = 𝑟𝑡𝑜 𝐸𝑐. 11
𝑟𝑡𝑜 ∗ 𝐿𝑛
𝑟𝑡𝑖
𝑈 ∗ = 1139.608338
𝑄
𝑇𝑠𝑢𝑟𝑓 ∗ = 𝑇𝑠 − (𝐸𝑐. 12)
𝐴𝑈 ∗
𝑇𝑠𝑢𝑟𝑓 ∗ = 598.031°𝐹
Dado que la diferencia entre Tsuf supuesto (550 °F) y Tsurf calculado (598.031 °F) difieren de una aproximación de
0.1°, los cálculos deben repetirse tomando en cuenta como nuevo Tsuf supuesto el calculado (598.031 °F), se deben
calcular todas las propiedades del aire nuevamente, los coeficiente de convección y radiación (hc y hr), el coeficiente
de calor transferido (U) para posteriormente calcular las pérdidas de calor (Q) y obtener un nuevo Tsurf calculado. Se
proceden las iteraciones hasta que la diferencia sea de 0.1°. A continuación se muestra las iteraciones:
Tabla Final
Kha cpa hr U
(BTU/pie- µa (Lb/pie- (BTU/Lb- νa (BTU/pie2- hc (BTU/pie2- (BTU/pie2-
Iteraciones Tsur °F Tsur* °F Tavg (°F) hr-°F) hr) °F) ρ (Lb/pie3) βa (1/°F) (pie2/hr) hr-°F) hr-°F) hr-°F) Q (BTU/h)
1 550.000 598.031 325.000 0.02086 0.05872 0.24380 0.05023 0.00127 1.16901 2.86949 1.63544 4.48719 1762.11720
2 598.031 597.887 349.016 0.02139 0.06000 0.24430 0.04861 0.00123 1.23422 3.17753 1.65948 4.81657 1891.46177
3 597.887 597.887 348.943 0.02139 0.05999 0.24430 0.04861 0.00123 1.23402 3.17657 1.65941 4.81555 1891.06090
Como se muestra en la tabla las pérdidas de calor no se ven afectadas por la conductividad del acero en gran medida, considerando este efecto las pérdidas de calor son
1891.06 BTU/pie.hora y sin considerar la conductividad del acero las pérdidas de calor son 1905 BTU/pie.hora, lo que genera un diferencia del 0.73 % por lo que el
mecanismo de convección por la conductividad del acero puede ser despreciable. Además esto puede ser fundamentando estudiándolo por la temperatura en la superficie de la
tuberia, ya que la temperatura real en la superficie de la tubería es de 597.887 °F una diferencia del 0.50% con respecto a la temperatura del vapor, por lo que se puede asumir
que Tsuf ≈ Ts
Considerando el efecto del Viento
Ya que se demostró en la parte anterior que las pérdidas de calor considerando o no considerando la conductividad del
hacer son relativamente iguales, en esta sección se estudiara sin considerar la conductividad del acero.
Calculo de hc
Va = 30 millas/hr
0.805
𝐾𝑕𝑎 5280 ∗ 𝜌𝑎 ∗ 𝑣𝑎 ∗ 𝑑𝑒
𝑕𝑐 = 0.0239 ∗ ∗ 8000 ∗ 𝑑𝑒 ∗ 𝑣𝑎 ≥ 50000 𝐸𝑐. 13
𝑑𝑒 𝜇𝑎
8000 ∗ 𝑑𝑒 ∗ 𝑣𝑎 = 60000
Propiedades del aire Tavg = 350°F
Propiedades del aire Khs =0

Kha (BTU/pie-hr-°F) 0.021408243
µa (Lb/pie-hr) 0.060048
cpa (BTU/Lb-°F) 0.244323075
ρ (Lb/pie3) 0.04854633
βa (1/°F) 0.001227614
νa (pie2/hr) 1.236921506
0.805
0.021408 5280 ∗ 0.04855 ∗ 30 ∗ 0.25
𝑕𝑐 = 0.0239 ∗ ∗
0.25 0.06005
𝐵𝑇𝑈
𝑕𝑐 = 8.6667
𝑕𝑟 − 𝑝𝑖𝑒 3 − °𝐹

−1
1
𝑈𝑡𝑜 =
8.6667 + 2.86949
𝐵𝑇𝑈
𝑈𝑡𝑜 = 11.53619
𝑕 − 𝑝𝑖𝑒 2 − °𝐹
𝑄
= 2𝜋 ∗ 0.125 ∗ 11.53619 ∗ 600 − 100
∆𝐿
𝐵𝑇𝑈
𝑄 = 4530.25
𝑕𝑟 − 𝑝𝑖𝑒
Con este resultado el mecanismo de convección aumenta notablemente con la velocidad del viento (hc= 8.667
BTU/hr-pie3-°F) y este tiene una alta incidencia en las perdidas de calor ya que aumentan las perdidas en un 237 %
en comparación con el cálculo hecho sin considerar el efecto del viento (1905 BTU/hr-pie) por lo que existe una
mayor transferencia de calor entre la superficie de la tubería y el aire de 30 millas/hora.
POZO
Emisibidad
−1
1 𝑟𝑡𝑜 1
𝜀= + ∗ −1
𝜀𝑡𝑜 𝑟𝑐𝑖 𝜀𝑐𝑖
o Tubería Aislada
−1
1 𝑟𝑖𝑛𝑠 1
𝜀= + ∗ −1
𝜀𝑖𝑛𝑠 𝑟𝑐𝑖 𝜀𝑐𝑖
𝜀𝑡𝑜 : emisibidad en la superfic
ie externa de la tubería (𝜀𝑡𝑜 ≈ 𝜀𝑐𝑖 )
Coeficiente de convección
0,049 ∗ 𝐾𝑕𝑎 ∗ (𝐺𝑟 ∗ 𝑃𝑟)0,333 ∗ (𝑃𝑟)0,074
𝑕𝑐 =
𝑟𝑐𝑖
𝑟𝑡𝑜 ∗ 𝐿𝑛 𝑟𝑡𝑜
(𝑟𝑐𝑖 − 𝑟𝑡𝑜)3 ∗ 𝑔 ∗ 𝜌𝑎2 ∗ 𝛽𝑎 ∗ (𝑇𝑡𝑜 − 𝑇𝑐𝑖) 𝐶𝑝𝑎 ∗ 𝜇𝑎

𝐺𝑟 = 𝑃𝑟 =
𝜇𝑎2 𝐾𝑕𝑎
o Tubería Aislada
0,049 ∗ 𝐾𝑕𝑎 ∗ (𝐺𝑟 ∗ 𝑃𝑟)0,333 ∗ (𝑃𝑟)0,074
𝑕𝑐 =
𝑟𝑐𝑖
𝑟𝑖𝑛𝑠 ∗ 𝐿𝑛
𝑟𝑖𝑛𝑠
(𝑟𝑐𝑖 − 𝑟𝑖𝑛𝑠)3 ∗ 𝑔 ∗ 𝜌𝑎2 ∗ 𝛽𝑎 ∗ (𝑇𝑖𝑛𝑠 − 𝑇𝑐𝑖) 𝐶𝑝𝑎 ∗ 𝜇𝑎
𝐺𝑟 = 𝑃𝑟 =
𝜇𝑎2 𝐾𝑕𝑎
𝑝𝑖𝑒 2
g: constante de la gravedad 4,17x108
𝑕𝑟
Tto: Temperatura en la superficie de la tubería, °F
Tci: Temperatura en la superficie interna del revestidor, °F(Tci = Tco)
Coeficiente de transferencia de calor
𝑟 𝑟 −1
𝑟𝑡𝑜 ∗ 𝐿𝑛 𝑡𝑜 𝑟𝑡𝑜 ∗ 𝐿𝑛 𝑖𝑛𝑠 𝑟𝑐𝑜 𝑟𝑕
𝑟𝑡𝑜 𝑟𝑡𝑖 𝑟𝑡𝑜 𝑟𝑡𝑜 𝑟𝑡𝑜 ∗ 𝐿𝑛 𝑟𝑡𝑜 ∗ 𝐿𝑛
𝑈𝑡𝑜 = + + + + 𝑟𝑐𝑖 + 𝑟𝑐𝑜
𝑟𝑡𝑖 ∗ 𝑕𝑓 𝐾𝑕𝑠 𝐾𝑕𝑖𝑛𝑠 𝑟𝑖𝑛𝑠 ∗ 𝑕𝑟 + 𝑕𝑐 𝐾𝑕𝑠 𝐾𝑕𝑐𝑒𝑚
𝐵𝑇𝑈
Khs: conductividad del acero
𝐵𝑇𝑈 𝑟𝑐𝑜 : radio externo del revestidor

Khins: Conductividad del aislante
𝑟𝑐𝑖 : radio interno del revestidor
𝐵𝑇𝑈
Khcem: Conductividad del cemento 𝑟𝑕 : radio del pozo
𝑕𝑟 ∗𝑝𝑖𝑒 ∗°𝐹
Temperatura del Pozo (Th)
𝐾𝑕𝑒 ∗ 𝑇𝑒
𝑇𝑠 ∗ 𝑓 𝑡 + 𝑟𝑡𝑜 ∗ 𝑈𝑡𝑜
𝑇𝑕 =
𝐾𝑕𝑒
𝑓 𝑡 + 𝑟𝑡𝑜 ∗ 𝑈𝑡𝑜
Te: Temperatura de la formación, °F

𝐵𝑇𝑈
Khe: Conductividad de la formación
Función Tiempo F(t)

𝑎∗𝑡
𝑓 𝑡 = 𝐿𝑛 − 0,29
𝑟𝑕
o Tubería Aislada (Khe = Kcem)
𝑎∗𝑡
𝑓 𝑡 = 𝐿𝑛 − 0,29
𝑟𝑐𝑜
Temperatura Externa del revestidor (Tco)
𝑟𝑕
𝑟𝑡𝑜 ∗ 𝑈𝑡𝑜 ∗ 𝐿𝑛 𝑟𝑐𝑜 ∗ (𝑇𝑠 − 𝑇𝑕)
𝑇𝑐𝑜 = 𝑇𝑕 +
𝐾𝑕𝑐𝑒𝑚
𝐾𝑕𝑖𝑛𝑠 ∗ 𝐹(𝑡) ∗ (𝑇𝑠 − 𝑇𝑖𝑛𝑠)
𝑇𝑐𝑜 = 𝑇𝑒 +
𝑟𝑖𝑛𝑠
𝐾𝑕𝑒 ∗ 𝐿𝑛 𝑟𝑐𝑜
o Tubería Aislada (Khe ≠ Kcem)
𝑟𝑕
Tco: Temperatura en la superficie externa del revestidor, °F(Tci = Tco)
Th: Temperatura del pozo, °F
Temperatura en la superficie del Aislante, Tins (Khe =Khcem)
o Tubería Aislada (Khe ≠ Kcem)

𝑟𝑖𝑛𝑠
𝑟𝑡𝑜 ∗ 𝑈𝑡𝑜 ∗ 𝐿𝑛( )
𝑇𝑖𝑛𝑠 = 𝑇𝑠 − 𝑟𝑡𝑜 𝑇𝑠 − 𝑇𝑕
Perdidas de calor (BTU / hr)
𝑄1 = 2 ∗ 𝜋 ∗ 𝑟𝑡𝑜 ∗ 𝑈𝑡𝑜 ∗ 𝑇𝑠 − 𝑇𝑕 ∗ ∆𝐿
2 ∗ 𝜋 ∗ 𝐾𝑕𝑒 ∗ 𝑇𝑕 − 𝑇𝑒 ∗ ∆𝐿 2 ∗ 𝜋 ∗ 𝐾𝑕𝑐𝑒𝑚 ∗ 𝑇𝑐𝑜 − 𝑇𝑕 ∗ ∆𝐿

𝑄2 = 𝑄3 =
𝑓(𝑡) 𝑟𝑕
𝐿𝑛
𝑟𝑐𝑜
𝑄1 = 2 ∗ 𝜋 ∗ 𝑟𝑡𝑜 ∗ 𝑈𝑡𝑜 ∗ 𝑇𝑠 − 𝑇𝑐𝑜 ∗ ∆𝐿
2 ∗ 𝜋 ∗ 𝐾𝑕𝑒 ∗ 𝑇𝑐𝑜 − 𝑇𝑒 ∗ ∆𝐿 𝑄3 =
2∗𝜋∗𝐾𝑕𝑐𝑒𝑚 ∗ 𝑇𝑡𝑜−𝑇𝑖𝑛𝑠 ∗∆𝐿
𝑄2 = 𝐿𝑛
𝑟𝑖𝑛𝑠
𝑓(𝑡) 𝑟𝑡𝑜
Perdidas de calor por RAMEY (BTU / hr)
2 ∗ 𝜋 ∗ 𝑟𝑡𝑜 ∗ 𝑈𝑡𝑜 ∗ 𝐾𝑕𝑒 𝐺𝑔 ∗ 𝑍 2

𝑄= ∗ 𝑇𝑠 − 𝑇𝑎 ∗ 𝑍 −
𝐾𝑕𝑒 + (𝑟𝑡𝑜∗ 𝑈𝑡𝑜 ∗ 𝑓(𝑡) 2
Ta: Temperatura ambiente, °F°
Z: profundidad, pies
Gg: Gradiente geotérmico
Calidad (X)
𝑄
𝑓𝑠𝑡,𝐹𝑜𝑛𝑑𝑜 𝑑𝑒𝑙 𝑝𝑜𝑧𝑜 = 𝑓𝑠𝑡,𝐶𝐴𝐵𝐸𝑍𝐴𝐿 −
𝑊 ∗ 𝐿𝑣
4.3 Vapor a 600°F es inyectado en un pozo a través de la tubería de inyección de 3 pulg. El pozo fue completado
con un revestidor de 9 pulg, 53/pie, N-60, en un hoyo de 12 pulg. La profundidad del pozo es 1000 pies y la
temperatura de la tierra es 100°F. Calcular las pérdidas de calor en el pozo después de 21 días de inyección
continua. Se dispone de la siguiente información adicional:
Rto=0.146 pies rci=0.55 pies
Rco=0.400 pies rh=0.5 pies
α=0.0286 pie2/h khe=1 BTU/h-°F-pie
εto=εci=0.9 Khcem=0.2 BTU/h-°F-pie
 Tubería desnuda
- Tiempo de inyección para t= 21 días= 504 horas
𝑎∗𝑡
𝑓 𝑡 = 𝐿𝑛 − 0,29 𝐸𝐶. 14
𝑟𝑕
2 ∗ 0,0286 ∗ 504
𝑓 𝑡 = 𝐿𝑛 − 0,29 = 2,43
0,500
Ts= Tsurf= 600 °F
Tci= 550 °F
−1
1 𝑟𝑡𝑜 1
𝜀= + ∗ −1 𝐸𝑐. 15
𝜀𝑡𝑜 𝑟𝑐𝑖 𝜀𝑐𝑖
1 0,146 1
𝜀=[ + ∗ − 1 ]−1
0,9 0,355 0,9
𝐸 = 0,864
Calculo de hr
𝑕𝑟 = 𝜍 ∗ 𝐸 ∗ 𝑇1∗2 + 𝑇2∗2 ∗ 𝑇1∗ + 𝑇2∗
𝑕𝑟 = 1,713𝑥10−9 ∗ 0,864 ∗ 1060∗2 + 1010∗2 ∗ 1060 + 1010
𝐵𝑇𝑈
𝑕𝑟 = 6,568
𝑕 𝑝𝑖𝑒 °𝐹
Calculo de las Propiedades del aire con las ecuaciones 1 hasta la 6
𝑇𝑡𝑜 + 𝑇𝑐𝑖
𝑇𝑎𝑣𝑔 =
2
600 + 550
𝑇𝑎𝑣𝑔 = = 575°𝐹
2
Propiedades del aire Sin Aislante

cpa (BTU/Lb-°F) 0.24959
ρ (Lb/pie3) 0.03788
βa (1/°F) 0.00096
νa (pie2/hr) 1.88400
Cálculo hc
𝐶𝑝𝑎 ∗ 𝜇𝑎
𝑃𝑟 = 𝐸𝑐. 16
𝐾𝑕𝑎
0,24959 ∗ 0,07135
𝑃𝑟 =
0,02608
𝑃𝑟 = 0,6828
(𝑟𝑐𝑖 − 𝑟𝑡𝑜)3 ∗ 𝑔 ∗ 𝜌𝑎2 ∗ 𝛽𝑎 ∗ (𝑇𝑡𝑜 − 𝑇𝑐𝑖)

𝐺𝑟 = 𝐸𝑐. 17
𝜇𝑎2
(𝑜, 355 − 0,146)3 ∗ 4,17𝑥108 ∗ 0,037872 ∗ 0,00096 ∗ (600 − 550)

𝐺𝑟 =
0,071352
𝐺𝑟 = 51477,661
0,049 ∗ 𝐾𝑕𝑎 ∗ (𝐺𝑟 ∗ 𝑃𝑟)0,333 ∗ (𝑃𝑟)0,074

𝑕𝑐 = 𝐸𝐶. 18
𝑟𝑐𝑖
0,049 ∗ 0,02608 ∗ (51477,661 ∗ 06828)0,333 ∗ (0,6828)0,074

𝑕𝑐 =
0,355
0,146 ∗ 𝐿𝑛 0,146
𝐵𝑇𝑈
𝑕𝑐 = 0,3126
𝑕 𝑝𝑖𝑒 °𝐹
- Cálculo de Uto
Para la cual se ha considerado despreciable la variación del coeficiente de transferencia de calor asociadas al material
de la tubería y del revestidor. En función a lo anterior, las siguientes consideraciones de temperaturas son necesarias:
Tciasumida=550 °F.
𝑟𝑕 −1
𝑟𝑡𝑜 ∗𝐿𝑛 1
𝑟𝑐𝑜
𝑈𝑡𝑜 = + 𝑕𝑐+𝑕𝑟
0,500
0,146 ∗ 𝐿𝑛 0,400 1
𝑈𝑡𝑜 = [ + ]−1
0,2 0,3126 + 6,568
𝐵𝑇𝑈
𝑈𝑡𝑜 = 3.23720
𝑕 − 𝑝𝑖𝑒2 − °𝐹
Cálculo de Th
𝑇 𝑠 ∗ 𝑓 𝑡 + 𝑟𝑡𝑜 ∗ 𝑈𝑡𝑜
𝑇𝑕 =
𝐾𝑕𝑒
𝑓 𝑡 + 𝑟𝑡𝑜 ∗ 𝑈𝑡𝑜
1 ∗ 100
600 ∗ 2,43 + 0,146 ∗ 3.23720
𝑇𝑕 =
1
2,43 + 0,146 ∗ 3.23720
𝑇𝑕 = 367.2788593 °𝐹
- Cálculo de Tco
𝑟𝑕
0,500
0,146 ∗ 3.23720 ∗ 𝐿𝑛 0,400 ∗ 600 − 367.2788593
𝑇𝑐𝑜 = 367.2788593 +
0,2
𝑇𝑐𝑜 = 489.998 °𝐹
Dado que la diferencia entre Tci supuesto (550 °F) y Tco calculado (489.998 °F) difieren de una aproximación de
0.1°, los cálculos deben repetirse tomando en cuenta como nuevo Tci supuesto el calculado (489.998 °F) ya que Tco ≈
Tci, se deben calcular todas las propiedades del aire nuevamente, los coeficiente de convección y radiación (hc y hr),
el coeficiente de calor transferido (U) y obtener un nuevo Tco calculado. Se proceden las iteraciones hasta que la
diferencia sea de 0.1°. A continuación se muestra las iteraciones
Tabla Final
Tuberia Desnuda
Kha µa cpa hr hc U
Tavg (BTU/pie- (Lb/pie- (BTU/Lb- ρ βa νa (BTU/pie2- (BTU/pie2- (BTU/pie2-
Iteraciones Th Tci Tco (°F) hr-°F) hr) °F) (Lb/pie3) (1/°F) (pie2/hr) Gr Pr hr-°F) hr-°F) hr-°F)
1 367.278859 550.000 489.998 575.000 0.02608 0.07136 0.24959 0.03788 0.00096 1.88400 51351.75319 0.68279 6.56758 0.31241 3.23720
2 363.345986 489.998 484.259 544.999 0.02549 0.06992 0.24883 0.03892 0.00098 1.79631 127733.18303 0.68267 6.02736 0.41342 3.13656
3 362.951439 484.259 483.684 542.130 0.02543 0.06978 0.24875 0.03903 0.00099 1.78792 136028.52378 0.68266 5.97790 0.42122 3.12665
4 362.911732 483.684 483.626 541.842 0.02542 0.06977 0.24875 0.03904 0.00099 1.78708 136871.18252 0.68266 5.97296 0.42199 3.12565
Cálculo de Q
𝑄1 = 2 ∗ 𝜋 ∗ 𝑟𝑡𝑜 ∗ 𝑈𝑡𝑜 ∗ 𝑇𝑠 − 𝑇𝑕 ∗ ∆𝐿
𝑄1 = 2 ∗ 𝜋 ∗ 0,146 ∗ 3,12565 ∗ 600 − 362,911732 ∗ 1000
𝐵𝑇𝑈
𝑄1 = 679803.7587
𝑕𝑟
2 ∗ 𝜋 ∗ 𝐾𝑕𝑒 ∗ 𝑇𝑕 − 𝑇𝑎 ∗ ∆𝐿 𝑄3 =
2∗𝜋∗𝐾𝑕𝑐𝑒𝑚 ∗ 𝑇𝑐𝑜 −𝑇𝑕 ∗∆𝐿
𝑄2 = 𝑟𝑕
𝐿𝑛 𝑟𝑐𝑜
𝑓(𝑡)
2𝜋 ∗ 0,2 483.626 − 362.91 ∗ 1000
2 ∗ 𝜋 ∗ 362.91 − 100 ∗ 1000 𝑄3 =
𝑄2 = 0,5
2,43 𝐿𝑛
0,4
𝑄2 = 679803.7587 𝑄3 = 679803.7587
𝑕𝑟 𝑕𝑟
Ya que el valor obtenido de Q de las tres ecuaciones anteriores es relativamente similar demuestra la
veracidad de los resultados.
Tubería aislada (sin considerar el cemento)
En este caso además de suponer despreciable el efecto de los términos que contienen hf y Khs, se considera que el
cemento tiene iguales propiedades térmicas que la tierra (Khcem = Khe)
Aislante de conductividad térmica igual a 0,04 BTU/h-pie-°F y espesor de 1,0 pulg (ξ=0,9)
Con f(t) y ε calculada mediante con la siguiente ecuación:
2 ∗ 0,0286 ∗ 504
𝑓 𝑡 = 𝐿𝑛 − 0,29 = 2.654
0,400
1 0,229 1
𝜀=[ + ∗ − 1 ]−1 = 0,846
0,9 0,355 0,9
El procedimiento pide suponer un valor de Tins (Temperatura al radio externo del aislante) y calcular un Tco a partir
de ahí, mediante la ecuación:
Asumiendo un Tins = 400
0,04 ∗ 2.654 ∗ 𝐿𝑛 ∗ 600 − 400

𝑇𝑐𝑜 = 100 + = 147.01
0,229
1 ∗ 𝐿𝑛 0,146
Calculo de las Propiedades del aire
Propiedades del aire Sin Aislante

cpa (BTU/Lb-°F) 0.24277
ρ (Lb/pie3) 0.05408
βa (1/°F) 0.00137
νa (pie2/hr) 1.03397
Calculo de hc y hr
𝑕𝑟 = 2.35376 𝑕𝑐 = 0.53172
𝑕 𝑝𝑖𝑒 °𝐹 𝑕 𝑝𝑖𝑒 °𝐹

−1
𝑟𝑖𝑛𝑠
𝑟𝑡𝑜 +
𝑈𝑡𝑜 =
𝐾𝑖𝑛𝑠 𝑟𝑖𝑛𝑠(𝑕𝑐 + 𝑕𝑟)
𝐵𝑇𝑈
𝑈𝑡𝑜 = 0.53489
𝑕 − 𝑝𝑖𝑒2 − °𝐹
Calculo de la Temperatura en la superficie del aislante
𝑇𝑖𝑛𝑠 = 234.77 °𝐹
Dado que la diferencia entre Tins supuesto (400 °F) y Tins calculado (234.77 °F) difieren de una aproximación de
0.1°, los cálculos deben repetirse tomando en cuenta como nuevo Tci supuesto el calculado (234.77 °F). A
continuación se muestra las iteraciones:
Tabla de Resultados:
Tuberia Aislada (Sin considerar el cemento)

hr hc U
(BTU/pie2- (BTU/pie2- (BTU/pie2-
Iteraciones Th °F Tins °F Tci °F Tco °F Gr Pr hr-°F) hr-°F) hr-°F)
1 185.7329 400 146.947325 146.95 267801.852 0.68842 2.35376 0.53172 0.53489
2 182.5829 234.7699 185.73287 185.73 79098.701 0.69222 1.74649 0.32925 0.51135
3 182.9556 248.1889 182.582944 182.58 102141.686 0.69186 1.78876 0.36074 0.51412
4 182.9141 246.6013 182.955628 182.96 99504.731 0.69190 1.78370 0.35735 0.51381
Si se Compara con el ejercicio anterior, se observa una gran diferencia en la distribución del perfil de temperaturas a manera que se aleja del pozo. Mientras que en el anterior,
la temperatura externa del revestidor era de 484 °F, acá es de 183 °F debido al efecto de aislante alrededor de la tubería. A continuación son calculadas las pérdidas de calor
asociadas a este ejemplo:
Ya que el valor obtenido de Q de las tres ecuaciones anteriores es relativamente similar demuestra la
veracidad de los .
Resumen de resultados
Sin Aislante Con Aislante Diferencia (%)

Q (BTU/Hr) 679803.7587 196599.9092 345.7803014
Tco (°F) 483.684 182.9556283 264.3721788
Como se esperaba, las pérdidas de calor son ampliamente reducidas por la colocación de un aislante al pozo. Para este
caso la disminución representa el 346 % del valor original
Tubería aislada considerando el cemento
 Tiempo de inyección
El tiempo de inyección se calculó para el caso de la tubería desnuda para la cual se obtuvo:
𝑓 𝑡 = 2,43
Para iniciar los cálculos se asume que la temperatura de la superficie del aislante es de 500°F y que la temperatura
interna del revestidor es de 450°F.
Tins= 500°F Tci=450°F
 Coeficiente de transferencia de calor por radiación

−1 −1
1 𝑟𝑖𝑛𝑠 1 1 0,229 1
𝜀= + ∗ −1 = + ∗ −1
𝜀𝑡𝑜 𝑟𝑐𝑖 𝜀𝑐𝑖 0,9 0,355 0,9
𝜀 = 0.8454
𝑕𝑟 = 𝜍 ∗ 𝜀 ∗ 𝑇1∗2 + 𝑇2∗2 ∗ 𝑇1∗ + 𝑇2∗
𝑕𝑟 = 1,713𝑥10−9 ∗ 0,8454 ∗ 9602 + 9102 ∗ 960 + 910
𝑕𝑟 = 4.73825 𝐵𝑇𝑈 𝑝𝑖𝑒 2 − 𝑕𝑟 − °𝐹
 Coeficiente de transferencia de calor por convección natural
Para este cálculo se utilizan las siguientes ecuaciones:
0,049 ∗ 𝐾𝑕𝑎 ∗ (𝐺𝑟 ∗ 𝑃𝑟)0,333 ∗ (𝑃𝑟)0,074

𝑕𝑐 =
𝑟𝑐𝑖
𝑟𝑡𝑜 ∗ 𝐿𝑛 𝑟𝑖𝑛𝑠
(𝑟𝑐𝑖 − 𝑟𝑖𝑛𝑠)3 ∗ 𝑔 ∗ 𝜌𝑎2 ∗ 𝛽𝑎 ∗ (𝑇𝑖𝑛𝑠 − 𝑇𝑐𝑖)

𝐺𝑟 =
𝜇𝑎2
𝐶𝑝𝑎 ∗ 𝜇𝑎
𝑃𝑟 =
𝐾𝑕𝑎
Para el cálculo del coeficiente de transferencia de calor por convección se necesitan las propiedades del aire.
Propiedades fisicas del Aire
Para el cálculo de las propiedades del aire se necesita una temperatura promedio que se calcula a partir de la
siguiente fórmula:
𝑇𝑖𝑛𝑠 + 𝑇𝑐𝑖 500 + 450

𝑇𝑎𝑣𝑔 = = → 𝑇𝑎𝑣𝑔 = 475°𝐹
2 2
Una vez calculadas las propiedades del aire se sustituye en las ecuaciones planteadas:
0.24712 ∗ 0.0664
𝑃𝑟 = = 0.68288
0.2406
(0.355 − 0.229)3 ∗ 4,17𝑥108 ∗ 0,041762 ∗ 0,00106 ∗ (500 − 450)

𝐺𝑟 =
0.066482
𝐺𝑟 = 17240.16022
0,049 ∗ 𝐾𝑕𝑎 ∗ (17240.16022 ∗ 0.68288)0,333 ∗ (0.68288)0,074

𝑕𝑐 =
0.355
0.229 ∗ 𝐿𝑛 0.229
𝑕𝑐 = 0.25938 𝐵𝑇𝑈 𝑝𝑖𝑒 2 − 𝑕𝑟 − °𝐹
 Calculo de coeficiente de transferencia de calor total

−1
𝑟𝑖𝑛𝑠 𝑟𝑕
𝑟𝑡𝑜 ∗ 𝐿𝑛 𝑟𝑡𝑜 𝑟𝑡𝑜 ∗ 𝐿𝑛
𝑈𝑡𝑜 = 𝑟𝑡𝑜 + + 𝑟𝑐𝑜
𝐾𝑕𝑖𝑛𝑠 𝑟𝑖𝑛𝑠 ∗ (𝑕𝑐 + 𝑕𝑟) 𝐾𝑕𝑐𝑒𝑚
−1
0.229 0.5
0.146 ∗ 𝐿𝑛 0.146 0.146 ∗ 𝐿𝑛
𝑈𝑡𝑜 = 0.146 + + 0.4
0.04 0.229 ∗ (𝑕𝑐 + 𝑕𝑟) 0.2
𝑈𝑡𝑜 = 0.51568 𝐵𝑇𝑈

𝑝𝑖𝑒 2 − 𝑕𝑟 − °𝐹
 Calculo de Th
𝑇𝑠 ∗ 𝑓 𝑡 +
𝑇𝑕 = 𝑟𝑡𝑜 ∗ 𝑈𝑡𝑜
𝐾𝑕𝑒
𝑓 𝑡 +
𝑟𝑡𝑜 ∗ 𝑈𝑡𝑜
1 ∗ 100
600 ∗ 2.43 +
𝑇𝑕 = 0.146 ∗ 0.51568 = 177,329°𝐹
1
𝑓 𝑡 +
0.146 ∗ 0.51568
Calculo de la temperatura externa del aislante (Tins)
𝑟𝑖𝑛𝑠
𝑟𝑡𝑜 ∗ 𝑈𝑡𝑜 ∗ 𝐿𝑛( )
𝑇𝑖𝑛𝑠 = 𝑇𝑠 − 𝑟𝑡𝑜 𝑇𝑠 − 𝑇𝑕
𝑇𝑖𝑛𝑠 = 240.745032 °𝐹
 Calculo de la temperatura externa del revestidor (Tco)
𝑟𝑕
𝑟𝑡𝑜 ∗ 𝑈𝑡𝑜 ∗ 𝐿𝑛 ∗ (𝑇𝑠 − 𝑇𝑕)
𝑇𝑐𝑜 = 𝑇𝑕 + 𝑟𝑐𝑜
0.5
0.146 ∗ 0.51568 ∗ 𝐿𝑛 ∗ (600 − 177.329)
𝑇𝑐𝑜 = 177.329 + 0.146
0.2
𝑇𝑐𝑜 = 212.8345°𝐹
El Tco calculado se encuentra muy alejado al Tci asumido, por lo que se procede a iterar para obtener el valor real de
Tci ya que Tco≅Tci. Tambin se observa que el Tins calculado se encuentra muy alejado del Tins asumido
Tuberia Aislada (Considerando al cemento)

hr hc U
Tavg (BTU/pie2- (BTU/pie2- (BTU/pie2-
Iteraciones Th Tins Tci Tco (°F) hr-°F) hr-°F) hr-°F)
1 177.329292 500 450 212.834504 475 4.73825 0.25938 0.51568
2 172.0125 240.745 212.834504 205.076543 226.79 1.87726 0.26806 0.47426
3 172.578802 265.4457 205.076543 205.902859 235.261 1.95045 0.34344 0.47862
4 172.527011 262.8148 205.902859 205.827288 234.359 1.94247 0.33709 0.47822
5 172.531807 263.0554 205.827288 205.834286 234.441 1.94320 0.33768 0.47826
6 172.531363 263.0331 205.834286 205.833639 234.434 1.94313 0.33763 0.47826
Se puede observar que en la iteración número 5, los valores de Tco y Tci son aproximadamente iguales por lo que se
toma este último como el valor real de Tci. De igual forma se toma el valor te Tins.
Se procede entonces al cálculo de las pérdidas de calor:
𝑄 = 2 ∗ 𝜋 ∗ 𝑟𝑡𝑜 ∗ 𝑈𝑡𝑜 ∗ 𝑇𝑠 − 𝑇𝑕 ∗ ∆𝐿
𝑄 = 2 ∗ 𝜋 ∗ 0.146 ∗ 0.47826 ∗ 600 − 172.53 ∗ 1000
𝑄 = 187529.2224 𝐵𝑇𝑈 𝑕𝑟
2 ∗ 𝜋 ∗ 𝐾𝑕𝑒 ∗ 𝑇𝑕 − 𝑇𝑎 ∗ ∆𝐿
𝑄=
𝑓(𝑡)
2 ∗ 𝜋 ∗ 1 ∗ 172.53 − 100 ∗ 1000

𝑄=
2.43
𝑄 = 187542.3852 𝐵𝑇𝑈 𝑕𝑟
2 ∗ 𝜋 ∗ 𝐾𝑕𝑐𝑒𝑚 ∗ 𝑇𝑐𝑜 − 𝑇𝑕 ∗ ∆𝐿
𝑄=
𝑟𝑕
𝐿𝑛 𝑟𝑐𝑜
2 ∗ 𝜋 ∗ 0.2 ∗ 205,83 − 172,53 ∗ 1000

𝑄=
0.5
𝐿𝑛 0.4
𝑄 = 187542.3852 𝐵𝑇𝑈 𝑕𝑟
Ya que el valor obtenido de Q de las tres ecuaciones anteriores es relativamente similar demuestra la veracidad de los
resultados.
En conclusión en esta parte se demuestra la influencia de despreciar el cemento a la hora de realizar los cálculos
generando los siguientes resultados
Sin Cemento Con Cemento Error (%)

Q (BTU/Hr) 196599.9092 187531.1318 4.835878383
Tco (°F) 182.956 205.833639 11.1148065
Tins (°F) 246.6012542 263.0330992 6.247063577
Th (°F) 182.9140565 172.5313633 6.017858475
En el cuadro se observa las diferencias entre considerar o no el

cemento que a pesar en las perdidas de calor no se observa una
diferencia considerable (error < 5%) si se observa en la caídas de
temperatura, pero se demuestra lo que establecen los fundamentos
teóricos que a través del cemento si existe una caída de
temperatura como se observa en la figura. Por lo que considerar el
efecto del cemento para el cálculo de las pérdidas de calor
generara valores más exactos y dentro de la realidad
Capítulo V, VII y VIII. Inyección cíclica y continua de vapor.
Capitulo V. Calentamiento de la formación por fluidos calientes
Modelo de Marx y Langenheim

𝑄 ∗ 𝑀𝑠 ∗ 𝑕𝑡 ∗ 𝐹1 𝑡𝑑
𝐴𝑠 = 𝐹1 =
4 ∗ 𝐾𝑕𝑜𝑏 ∗ 𝑀𝑜𝑏 ∗ 𝑇𝑠 − 𝑇𝑟 1 + 0.85 ∗ 𝑡𝑑
𝐴𝑠 4 ∗ 𝐾𝑕𝑜𝑏 ∗ 𝑀𝑜𝑏 ∗ 𝑡
𝑟𝑕 = 𝑡𝑑 =
𝜋 𝑀𝑠 2 𝑕2
As: Area calentada, piee ht: espesor total, pies
rh: radio calentado, pie Tr: temperatura de la formación, °F

𝐿𝑏 Td: tiempo adimensional
Ms: Capacidad calorífica de la formación
𝑝𝑖𝑒 −°𝐹
t: tiempo de inyección, horas
Mob: Capacidad calorífica de las capas supra y
𝐿𝑏
subyacentes Khob: Conductividad térmica de las capas sub y
𝑝𝑖𝑒 −°𝐹 𝐿𝑏
suprayacentes,
𝑕𝑟−𝑝𝑖𝑒 −°𝐹
F1: Función de tiempo adimensional
Modelo de Mandl y Volek
𝑄𝑖 ∗ 𝑀𝑠 ∗ 𝑕𝑡 2 ∗ 𝐹3 𝑇𝐷𝐶 ∗ 𝑀𝑠 2 ∗ 𝑕𝑡 2
𝑉𝑠 = 𝑇𝑐 =
4 ∗ 𝐾𝑕𝑜𝑏 ∗ 𝑀𝑜𝑏 ∗ (𝑇𝑠 − 𝑇𝑟) 4 ∗ 𝐾𝑕𝑜𝑏 ∗ 𝑀𝑜𝑏
1 𝑓𝑠 ∗ 𝐿𝑣
𝐹2𝑐 = 𝐵=
1+𝐵 𝐶𝑤 ∗ (𝑇𝑠 − 𝑇𝑟)
Vs: Volumen calentado, pie3 Tdc: tiempo adimensional critico (tabla 5.1)
F3: Función de tiempo adimensional Cw: conductividad térmica del agua𝐵𝑇𝑈 𝑙𝑏. ℉
Tasa de Inyección al tiempo (BTU/Hr)

350
𝑄𝑖 = ( )𝑖 ∗ 𝐶𝑤 ∗ 𝑇𝑠 − 𝑇𝑟 + 𝑓𝑠 ∗ 𝐿𝑣
24 𝑠𝑡
𝑖𝑠𝑡 : tasa de inyección B/d
Capítulo VII. Inyección continúa de vapor
Petróleo total recuperado de la zona de vapor, BN (Np)

∅ ∗ 𝑉𝑠 ∗ 𝑆𝑜𝑖 − 𝑆𝑜𝑟 𝑕𝑛
𝑁𝑝 = ∗ 𝐸𝑐 ∗
5,615 𝑕𝑡
Vs: Volumen calentado, piee Sor: Saturación de petróleo residual
Soi: Saturación de petróleo inicial Ec: Factor de Captura
Tasa de producción de petróleo (BN / día)
𝑄𝑖 ∗ ∅ ∗ (𝑆𝑜𝑖 − 𝑆𝑜𝑟) 𝑕𝑛 1 + 𝑂, 425 ∗ 𝑇𝐷

𝑄𝑜 = ∗ 𝐸𝑐 ∗ ∗ 𝐹2 𝐹2 =
𝑀𝑠 ∗ (𝑇𝑠 − 𝑇𝑟) 𝑕𝑡
(1 + 0,85 ∗ 𝑇𝐷 )2
Razón petróleo / Vapor acumulada

𝑁𝑝
𝑂𝑆𝑅 =
𝑖𝑠𝑡 ∗ 𝑡
Razón petróleo / Vapor instantánea
𝑄𝑜 ∗ 𝑡
𝑂𝑆𝑅 =
𝑖𝑠𝑡
Capítulo VIII. Inyección cíclica de vapor
Tasa de flujo de petróleo antes de la estimulación (BN / dia)
2 ∗ 𝜋 ∗ 1.127 ∗ 𝑕𝑘 𝑝𝑒 − 𝑝𝑤
𝑞𝑜𝑐 = 𝑟𝑒
𝜇𝑜𝑐 ∗ 𝐿𝑛
𝑟𝑤
Tasa de petróleo estimulado (BN / día)
2 ∗ 𝜋 ∗ 1.127 ∗ 𝑕𝑘 𝑝𝑒 − 𝑝𝑤
𝑞𝑜𝑕 =
𝑟𝑒 𝑟𝑕
𝜇𝑜𝑐 ∗ 𝐿𝑛 + 𝜇𝑜𝑕 ∗ 𝐿𝑛
𝑟𝑤 𝑟𝑤
Espeso aumentado de la formación
𝑕𝑡 ∗ 𝑡𝑑
𝑕=
𝐹1
𝐾𝑕𝑜𝑏 ∗ 𝑡
𝑥 = 𝑙𝑜𝑔
𝑀𝑜𝑏 ∗ 𝑟𝑕2
𝑣𝑟 = 0.180304 − 0.41269 ∗ 𝑥 + 0.18217 ∗ 𝑥 2 + 0.149516 ∗ 𝑥 3 + 0.024183 ∗ 𝑥 4
4 ∗ 𝐾𝑕𝑜𝑏 ∗ 𝑡
𝑦 = 𝑙𝑜𝑔
𝑀𝑜𝑏 ∗ 𝑕2
𝑣𝑧 = 0.474884 − 0.56832 ∗ 𝑦 − 0.239719 ∗ 𝑦 2 − 0.035737 ∗ 𝑦 3
Calor removido por el petróleo y el gas

𝐻𝑜, 𝑔 = 5.615 ∗ 𝜌𝑜 ∗ 𝑐𝑜 + 𝑅𝑔 ∗ 𝑐𝑔 ∗ 𝑇𝑎𝑣𝑔 − 𝑇𝑟
Calor removido por el agua y vapor
𝐻𝑤 ,𝑔 = 5.615 ∗ 𝑝𝑤 ∗ (𝑅𝑤 ∗ 𝐻𝑤 − 𝐻𝑟 + 𝑅𝑠 ∗ 𝐿𝑣)
1/0.225
337.10
𝑝𝑠 =
115.10
Agua producida en estado de vapor por barril de petróleo producido
Si Ps > Pw, entonces todo el agua es producida como vapor, si Ps< Pw, Rs debe ser calculado por la siguiente
ecuación:
1 𝜌𝑔 ∗ 𝑅𝑔 𝑃𝑠 𝑃𝑠
𝑅𝑠 = ∗ ∗ ≈ 1,6𝑥10−4
5.615 𝜌𝑤 𝑃𝑤 − 𝑃𝑠 𝑃𝑤 − 𝑃𝑠
Energía removida por medio de los fluidos producidos

1 𝑡 𝐻𝑓
𝜕= 2 𝑑𝑡
2 0 𝜋 ∗ 𝑟𝑕 ∗ 𝑕 ∗ 𝑀𝑠 ∗ (𝑇𝑠 − 𝑇𝑔)
Temperatura después de la estimulación
𝑇𝑎𝑣𝑔 = 𝑇𝑟 + 𝑇𝑠 − 𝑇𝑟 ∗ ((𝑉𝑟 − 𝑉𝑧) ∗ 1 − 𝜕 − 𝜕)
7.1. Vapor a 300 lpca y calidad igual a 80%, está siendo inyectado en una formación de espesor (total) igual a 35
pies y porosidad de 28%, a una tasa de 400 B/D (equivalente de agua). La capacidad calorífica de la arena y de
las capas supra y subyacentes es 33 BTU/pie3 F, la saturación inicial de petróleo es 65%, la saturación residual
de petróleo es 20%, la conductividad térmica de las formaciones adyacentes es 0,7 BTU/ h pie2 F y la
temperatura de la formación es de 85 F. Calcular el volumen de la zona de vapor, la recuperación total de
petróleo, la tasa de producción de petróleo y la razón petróleo/vapor al final de 100 días y 4.197,5 días luego de
iniciada la inyección. El calor específico del agua es de 1,02 BTU/ lb F. Considere que el espesor neto es igual
al espesor total y el factor de captura es igual al 100%.
100 DIAS
 Cálculo de la temperatura de saturación
𝑇𝑠 = 115,1 ∗ 𝑃𝑠 0,225
𝑇𝑠 = 115,1 ∗ 3000,225 = 415,3628 °𝐹
 Cálculo del calor latente de vaporización
𝐿𝑣 = 1318 ∗ 𝑃𝑠 −0,08774
𝐵𝑇𝑈
𝐿𝑣 = 1318 ∗ 300−0,08774 = 799,0479
𝑙𝑏𝑚
 Cálculo de B
𝑓𝑠 ∗ 𝐿𝑣
𝐵=
𝐶𝑤 ∗ (𝑇𝑠 − 𝑇𝑟)
0,8 ∗ 799,0479
𝐵= = 1,8970
1,02 ∗ (415,3628 − 85)
 Cálculo de F2c
1
𝐹2𝑐 =
1+𝐵
1
𝐹2𝑐 = = 0,3452
1 + 1,8970
Posteriormente con la tabla de funciones F1 y F2 de Marx y Langenheim se obtiene el valor de T DC interpolando con
F2, quedando así:
TDC F2
1,80 0,3496
1,8638 0,3452
1,90 0,3427
 Cálculo del tiempo
𝑇𝐷𝐶 ∗ 𝑀𝑠 2 ∗ 𝑕𝑡 2
𝑇𝑐 =
4 ∗ 𝐾𝑕𝑜𝑏 ∗ 𝑀𝑜𝑏
1,8638 ∗ 332 ∗ 352
𝑇𝑐 = = 26908,6125𝑕𝑟 = 1121,192 𝑑í𝑎𝑠
4 ∗ 0,7 ∗ 33
 Para 100 días= 2400 hr

T=100 días ≤ Tc=1121,192 𝑑í𝑎𝑠
Por lo tanto el método que se aplicará es el de Max y Langenheim
 Cálculo de la tasa de inyección de calor
𝑄𝑖 = 𝑖𝑠𝑡 ∗ 𝐶𝑤 ∗ 𝑇𝑠 − 𝑇𝑟 + 𝑓𝑠 ∗ 𝐿𝑣
𝐵 𝐵𝑇𝑈 𝐵𝑇𝑈 𝑙𝑏 1𝐷
𝑄𝑖 = 400 ∗ 1,02 ∗ 415,3628 − 85 𝐹 + 0,8 ∗ 799,0479 ∗ 350 ∗
𝐷 𝐿𝐵 𝐹 𝑙𝑏𝑚 𝐵𝑙𝑠 24𝑕
𝐵𝑇𝑈
𝑄𝑖 = 5,6945 𝑀𝑀
𝑕𝑟
 Cálculo de TD
4 ∗ 𝐾𝑕𝑜𝑏 ∗ 𝑀𝑜𝑏
𝑇𝐷 = ∗𝑡
𝑀𝑠 2 ∗ 𝑕𝑡 2
4 ∗ 0,7 ∗ 33
𝑇𝐷 = ∗ 2400 𝑕𝑟
332 ∗ 352
𝑇𝐷 = 0,1662
 Cálculo de F1
0,1662
𝐹1 =
1 + (0,85 ∗ 0,1662)
𝐹1 = 0,1234
 Cálculo del área de la zona de vapor calentada
𝑄𝑖 ∗ 𝑀𝑠 ∗ 𝑕𝑡
𝐴𝑠 =
4 ∗ 𝐾𝑕𝑜𝑏 ∗ 𝑀𝑜𝑏 ∗ (𝑇𝑠 − 𝑇𝑟)
5,6945𝑥106 ∗33∗35
𝐴𝑠 = 4∗0,7∗33∗(415,3628−85)
𝐴𝑠 = 26599,3904 𝑝𝑖𝑒 2
 Cálculo del volumen de la zona de vapor
𝑉𝑠 = 𝐴𝑠 ∗ 𝑕𝑡
𝑉𝑠 = 26599,3904 ∗ 35
𝑉𝑠 = 930978,662 𝑝𝑖𝑒 3
 Cálculo de la recuperación total de petróleo
∅ ∗ 𝑉𝑠 ∗ (𝑆𝑜𝑖 − 𝑆𝑜𝑟) 𝑕𝑛
𝑁𝑝 = ∗ 𝐸𝑐 ∗
5,615 𝑕𝑡
0,28 ∗ 930978,662 ∗ (0,65 − 0,20) 35

𝑁𝑝 = ∗1∗
5,615 35
𝑁𝑝 = 20891,0617 𝐵𝑁
 Cálculo tasa de producción de petróleo
𝑄𝑖 ∗ ∅ ∗ (𝑆𝑜𝑖 − 𝑆𝑜𝑟) 𝑕𝑛
𝑄𝑜 = ∗ 𝐸𝑐 ∗ ∗ 𝐹2
𝑀𝑠 ∗ (𝑇𝑠 − 𝑇𝑟) 𝑕𝑡
- Buscando F2 = ∂F1/∂TD
1 + 𝑂, 425 ∗ 𝑇𝐷
𝐹2 =
(1 + 0,85 ∗ 𝑇𝐷 )2
1 + 𝑂, 425 ∗ 0,1662
𝐹2 =
(1 + 0,85 ∗ 0,1662)2
𝐹2 = 0,6471
5,6945𝑥106 ∗ 0,28 ∗ (0,65 − 0,20) 35

𝑄𝑜 = ∗1∗ ∗ 0,6471
33 ∗ (415,3628 − 85) 35
𝑄𝑜 = 182,0273 𝐵𝑁𝐷
 Cálculo de la razón de petróleo-vapor acumulado
𝑁𝑝
𝑂𝑆𝑅 =
𝑖𝑠𝑡 ∗ 𝑡
20891,0617
𝑂𝑆𝑅 = = 0,5223
400 ∗ 100
Con este valor de la razón de petróleo-vapor acumulado se demuestra que el proyecto es exitoso.
 Cálculo del volumen de la zona de vapor, la recuperación total de petróleo, la tasa de producción de petróleo y
la razón petróleo/vapor a diferentes tiempos.
Qo, (calc)
T (días) Vs (pie3) Np (BND) OSR
(BND)
100 930978,6623 20891,0617 182,0273 0,5223
250 2024630,841 45432,4997 149,5649 0,4543
400 2961676,852 66459,7121 132,1406 0,4154
550 3803548,293 85351,2173 120,3954 0,3880

700 4577848,037 102726,4208 111,6543 0,3669
850 5300350,544 118939,2998 104,7657 0,3498
1000 5981235,789 134218,2920 99,1281 0,3355
160000
140000
120000
100000
Np (BN)
80000
Np (BN)
60000
40000
20000
0
0 200 400 600 800 1000 1200
tiempo (dias)
200
180
160
140
120
qo (B/D)
100
80 Qo (BND)
60
40
20
0
0 200 400 600 800 1000 1200
tiempo (DIAS)
Como se observa en las figuras, a medida que aumenta el tiempo de inyección la producción de crudo
aumenta debido a que el volumen de la zona de vapor es mayor, y por lo tanto se desplaza mayor cantidad de
crudo, pero caso distinto ocurre con la tasa, que a medida que pasa el tiempo esta disminuye debido a que
inicialmente la tasa es mayor debido a que se produce el crudo en las cercanías del pozo y poco a poco a
medida que pasa el tiempo la tasa decrece hasta un límite económico
4.197,5 DÍAS
Debido a que el tiempo de estudio es mayor al tiempo crítico, el volumen de la zona de vapor se calcula por el método
de Mandl y Volek
 Calculo del volumen de la zona de vapor
Para calcular el volumen de la zona de vapor se necesita del área que viene dada por la siguiente formula:
𝑄𝑖 ∗ 𝑀𝑠 ∗ 𝑕𝑡
𝐴𝑠 = ∗ 𝐹3
4 ∗ 𝐾𝑕𝑜𝑏 ∗ 𝑀𝑜𝑏 ∗ (𝑇𝑠 − 𝑇𝑟)
Para el cálculo de F3 se requiere de Td:
4 ∗ 𝐾𝑕𝑜𝑏 ∗ 𝑀𝑜𝑏 4 ∗ 0,7 ∗ 33

𝑇𝑑 = 2 2
∗𝑡 = ∗ 4197,5 → 𝑇𝑑 = 6,9777
𝑀𝑠 ∗ 𝑕𝑡 332 ∗ 352
Con el valor de Td se entra en la gráfica 5.5 para obtener el valor de F3.
𝐹3 = 2,05
15.459.665,13 ∗ 33 ∗ 35
𝐴𝑠 = ∗ 2,05 → 𝐴𝑠 = 441.704,7179 𝑝𝑖𝑒 2
4 ∗ 0,7 ∗ 33 ∗ (415,3628 − 100)
Se calcula finalmente el volumen de la zona de vapor
𝑉𝑠 = 𝐴𝑠 ∗ 𝑕𝑡 = 441.704,7179 ∗ 35 → 𝑉𝑠 = 15.459.665,13 𝑝𝑖𝑒 3
 Calculo de la recuperación acumulada de petróleo
La recuperación acumulada de petróleo se obtiene a partir de la siguiente ecuación:
Ø ∗ 𝑉𝑠 ∗ (𝑆𝑜𝑖 − 𝑆𝑜𝑟) 𝑕𝑛
𝑁𝑝 = ∗ 𝐸𝑐 ∗
5,615 𝑕𝑡
0,28 ∗ 15.459.665,13 ∗ (0,65 − 0,2) 35

𝑁𝑝 = ∗1∗ → 𝑁𝑝 = 346,913,2335 𝐵𝑁
5,615 35
 Calculo de la tasa de producción de petróleo
Debido a que la fórmula planteada para el cálculo de tasa no aplica para el método de Mandl y
Volek, se requiere calcular el Np a diferentes tiempos.
T (días) Td F3 (Grafica) As (pie2) Vs (pie3) Np (BN)

1500 2,49 1,05 226239,002 7918365,06 177687,266
2000 3,32 1,3 280105,431 9803690,08 219993,758
2500 4,16 1,5 323198,574 11311950,1 253838,951
3000 4,99 1,7 366291,717 12820210,1 287684,145
3500 5,82 1,85 398611,575 13951405,1 313068,04
4000 6,65 2 430931,432 15082600,1 338451,935
4197,5 6,98 2,05 441704,718 15459665,1 346913,233
4500 7,48 2,1 452478,004 15836730,1 355374,532
5000 8,31 2,3 495571,147 17344990,1 389219,725
5500 9,14 2,35 506344,433 17722055,1 397681,024
6000 9,97 2,5 538664,29 18853250,2 423064,919
450000
400000
350000
300000
250000
Np (BN)
200000
150000
100000
50000
0
0 1000 2000 3000 4000 5000 6000 7000
Tiempo (dias)
De la gráfica anterior se procede a obtener la pendiente de cada valor de Np en los diferentes tiempos, estos valore de
pendiente será la tasa.
T (días) Qo (BND)
1000
1500 94,01
2000 103,89
2500 66,1
100
3000 56,41
Log (q)
3500 60,44
4000 45,33
10
4197,5 34,77
4500 63,46
5000 15,11
1
5500 52,88
0 1000 2000 3000 4000 5000 6000
6000 49,36
tiempo de inyeccion (dias)
Debido a que la curva tiene una tendencia lineal con pendiente negativa, el grafico resultante aplica a una declinación
exponencial
8.1. En un proyecto de inyección cíclica de vapor se inyectaron 1.370 B/D de vapor (equivalente en agua) a una
presión de 1600 lpca y 80% de calidad en la cara de la arena productora. El espesor de la arena es de 55 pies,
la permeabilidad es de 2 darcy y el tiempo de inyección fué de 35 días. Calcular la tasa de producción
estimulada a los 10 días de terminada la inyección. La viscosidad del petróleo es de 156.210,4 cp a 60F y
220,13 cp a 200 F. La densidad del petróleo es 60 lbs/pie3, la densidad del agua es 62,4 lbs/pie3, el calor
especifico del petróleo es 0,5 BTU / Lb-°F, el calor especifico del gas es 0,01 BTU/PCN-°F, el calor
especifico del agua 1,0 BTU/ lb F, la relación gas/petróleo es 0,0 BN/PCN, y la relación agua/petróleo es 1,0
BN/BN. La presión en el fondo del pozo es 87 lpca, la presión en el límite exterior es 400 lpca y la
temperatura original del yacimiento es 70 F. La conductividad térmica y la difusividad térmica de las capas
supra y subyacentes son 1,4 BTU/h- pie- F y 0,04 h/pie2, respectivamente, y la capacidad calorífica de la
arena es 35 BTU/pie- F. El radio de drenaje es igual a 70 pies y el radio del pozo es de 0,333 pies. Considere
presión constante en el límite exterior del área de drenaje.
a) Calculo de Ts y Lv a 1600 lpca

𝑇𝑠 = 115,1 × 𝑃𝑠0,225
𝑇𝑠 = 605 °𝐹
𝐿𝑣 = 1318 × 𝑃𝑠−0,08774
𝐵𝑇𝑈
𝐿𝑣 = 690
𝐿𝑏
b) Calculo de Qi
350 ∗ 𝑖𝑠𝑡
𝑄𝑖 = ∗ 𝑐𝑤 ∗ 𝑇𝑠 − 𝑇𝑟 + 𝑋 ∗ 𝐿𝑣
24
350 ∗ 1370
𝑄𝑖 = ∗ 1 ∗ 605 − 70 + 0.8 ∗ 690
24
𝑄𝑖 = 21.71735 𝑀𝑀𝐵𝑇𝑈/𝑕
c) Calculo de td
4 ∗ 𝐾𝑕𝑜𝑏 ∗ 𝑀𝑜𝑏 ∗ 𝑡
𝑡𝑑 =
𝑀𝑠 2 𝑕2
4 ∗ 1.4 ∗ 35 ∗ 35 ∗ 24
𝑡𝑑 =
352 552
a) Calculo de F1
𝑡𝑑
𝐹1 =
1 + 0.85 ∗ 𝑡𝑑
𝐹1 = 0.0377
b) Calculo de As y rh
𝑄 ∗ 𝑀𝑠 ∗ 𝑕𝑡 ∗ 𝐹1
𝐴𝑠 =
4 ∗ 𝐾𝑕𝑜𝑏 ∗ 𝑀𝑜𝑏 ∗ 𝑇𝑠 − 𝑇𝑟
𝐴𝑠 = 15030 𝑝𝑖𝑒 2
𝐴𝑠
𝑟𝑕 =
𝜋
𝑟𝑕 = 69.168 𝑝𝑖𝑒
c) Suponiendo que la temperatura del yacimiento es igual a Tavg para iniciar los cálculos:
605 + 70
𝑇𝑎𝑣𝑔 = = 337.5 °𝐹
2
d) Calculo de las constante a y b, Ecuación de Andrade
𝑏
𝐿𝑛 156210.4 = 𝐿𝑛 𝑎 +
60 + 460
𝑏
𝐿𝑛 220.13 = 𝐿𝑛 𝑎 +
200 + 460
𝑎 = 5.66417𝑥10−9 𝑦 𝑏 = 1.60929𝑥104
e) Calculo de µoh a Tavg y µoc a Tr
−9 1.60929𝑥104 /797.1
𝜇𝑜𝑕 = 5.66417𝑥10 ∗𝑒 = 3.32 𝑐𝑝
−9 1.60929𝑥104 /530
𝜇𝑜𝑐 = 5.66417𝑥10 ∗𝑒 = 87.104 𝑐𝑝
f) Calculo de qoc
2 ∗ 𝜋 ∗ 1.127 ∗ 𝑕𝑘 𝑝𝑒 − 𝑝𝑤
𝑞𝑜𝑐 = 𝑟𝑒
𝜇𝑜𝑐 ∗ 𝐿𝑛
𝑟𝑤
𝐵
𝑞𝑜𝑐 = 0.654
𝐷
𝑡𝑑 = 0.0444
a) Calculo de qoh
2 ∗ 𝜋 ∗ 1.127 ∗ 𝑕𝑘 𝑝𝑒 − 𝑝𝑤
𝑞𝑜𝑕 =
𝑟𝑒 𝑟𝑕
𝜇𝑜𝑐 ∗ 𝐿𝑛 + 𝜇𝑜𝑕 ∗ 𝐿𝑛
𝑟𝑤 𝑟𝑤
𝑞𝑜𝑕 = 49.76 𝐵/𝐷
b) Calculo de vr y vz
𝑕𝑡 ∗ 𝑡𝑑
𝑕= = 64.8 𝑝𝑖𝑒𝑠
𝐹1
4 ∗ 𝐾𝑕𝑜𝑏 ∗ 𝑡
𝑦 = 𝑙𝑜𝑔 = −2.04
𝑀𝑜𝑏 ∗ 𝑕2
𝑣𝑧 = 0.474884 − 0.56832 ∗ 𝑦 − 0.239719 ∗ 𝑦 2 − 0.035737 ∗ 𝑦 3 = 0.94
𝐾𝑕𝑜𝑏 ∗ 𝑡
𝑥 = 𝑙𝑜𝑔 = −2.63
𝑀𝑜𝑏 ∗ 𝑟𝑕2
𝑣𝑟 = 0.180304 − 0.41269 ∗ 𝑥 + 0.18217 ∗ 𝑥 2 + 0.149516 ∗ 𝑥 3 + 0.024183 ∗ 𝑥 4 = 0.94
c) Calculo de ho,g
𝐻𝑜, 𝑔 = 5.615 ∗ 𝜌𝑜 ∗ 𝑐𝑜 + 𝑅𝑔 ∗ 𝑐𝑔 ∗ 𝑇𝑎𝑣𝑔 − 𝑇𝑟
𝐻𝑜, 𝑔 = 44.993 𝐵𝑇𝑈/𝐵𝑁
d) Calculo de Hws
1/0.225
337.10
𝑝𝑠 = = 118.62 𝑙𝑝𝑐𝑎
115.10
Como: ps>pw, entonces Ra=RW=1.0 BY/BN
e) Cálculo de HWS
𝐵𝑇𝑈
𝐻𝑊 = 1 ∗ 337.10 − 32 = 305.1
𝑙𝑏
𝐵𝑇𝑈
𝐻𝑟 = 1 ∗ 70 − 32 = 38
𝑙𝑏
𝐵𝑇𝑈
𝐿𝑣 = 1318 ∗ (118.6)−0.08774 = 866.8
𝑙𝑏
𝐻𝑤𝑎 = 5.615 ∗ 𝑝𝑤 ∗ (𝑅𝑤 ∗ 𝐻𝑤 − 𝐻𝑟 + 𝑅𝑠 ∗ 𝐿𝑣)
𝐻𝑤𝑎 = 5.615 ∗ 62.4 ∗ (1 ∗ 305.1 − 38 + 1 ∗ 866.8)

𝐵𝑇𝑈
𝐻𝑤𝑎 = 397.298
𝐵𝑁
d) Cálculo de σ
𝐵𝑇𝑈
𝐻𝑓 = 49.7 ∗ 44.993 + 397.298 = 22088.418
𝑙𝑏
𝑡
1 𝐻𝑓
𝜕= 𝑑𝑡
2 0 𝜋∗ 𝑟𝑕2 ∗ 𝑕 ∗ 𝑀𝑠 ∗ (𝑇𝑠 − 𝑇𝑔)
1 22088.418 ∗ 10
𝜕= = 6.79063𝑥10−3
2 2 ∗ 𝜋 ∗ (65.6)2 ∗ 64.8 ∗ 35 ∗ (604.19 − 70)
e) Cálculo de Tavg
𝑇𝑎𝑣𝑔 = 𝑇𝑟 + 𝑇𝑠 − 𝑇𝑟 ∗ ((𝑉𝑟 − 𝑉𝑧) ∗ 1 − 𝜕 − 𝜕)
𝑇𝑎𝑣𝑔 = 70 + 604 − 70 ∗ ((0.963 − 0.94) ∗ 1 − 6.79063𝑥10−3 − 6.79063𝑥10−3 )
𝑇𝑎𝑣𝑔 = 546.55 °𝐹
Dado que la diferencia entre Tavg supuesto (337.1 °F) y Tavg calculado (546.55 °F) difieren de una aproximación de
0.1°, los cálculos deben repetirse tomando en cuenta como nuevo Tavg supuesto el calculado (546.55 °F) Se proceden
las iteraciones hasta que la diferencia sea de 0.1°. A continuación se muestra las iteraciones:
Ho,g Hw,s
Tavg °F µoh (cp) qoh (B/d) (BTU/BN) Hw (BTU/Lb) (BTU/BN) Hf (BTU/d) ϭ Tavg* (°F)
337.1 3.3221176 49.7370051 44992.995 305.1 397291.3464 21997898.6 0.00671514 546.555584
546.555584 0.0497562 49.8775983 80275.78816 514.5555842 470679.5562 27480329.4 0.00838872 544.852904
544.852904 0.05112258 49.8775395 79988.97161 512.8529036 470082.9777 27436235.4 0.00837526 544.866598
IX. Combustión in Situ
Razón de las moléculas de CO2 a moléculas de CO (m)

%𝐶𝑂2
𝑚=
%𝐶𝑂
Razón de átomos de hidrogeno a átomos de carbono (n)
𝑂2𝑖
4∗ ∗ 𝑁2 − 𝐶𝑂2 + 0.5 ∗ 𝐶𝑂 + 𝑂2
1 − 𝑂2𝑖
𝐹𝐻𝐶 =
𝐶𝑂2 + 𝐶𝑂
Eficiencia de Utilización de Oxigeno
𝑂2𝑖 𝑂2
𝐸𝑜2 = 1 − ∗
1 − 𝑂2𝑖 𝑁2
Calor generado en la zona de combustión ∆Hc
1800
∆𝐻𝑐 = ∗ 94 − 67.9𝑚∗ + 31.2𝐹𝐻𝐶
12 + 𝐹𝐻𝐶
Requerimiento de aire para quemar el combustible depositado (α)
379 2𝑚 + 1 𝐹𝐻𝐶 𝐹𝑐
𝛼= ∗ + ∗
𝐸𝑜2 ∗ 𝑂2𝑖 2𝑚 + 2 4 12 + 𝐹𝐻𝐶
Flujo local de aire (Ua) Velocidad del frente de combustión (Vf)
𝑄𝑖 𝑄𝑖 𝑈𝑎
𝑈𝑎 = = 𝑉𝑓 =
2 ∗ 𝜋 ∗ 𝑟𝑓 ∗ 𝑕 𝐴 𝛼
Contenido de combustible
𝑚𝑓
𝐶𝑚 = 4 ∗ 𝑉𝑔 21 5
𝐿 𝐶𝑚 = ∗ ∗ %𝑁2 − %𝑂2 + 2 ∗ %𝐶𝑂2 + ∗ %𝐶𝑂
𝐴∗ 379 ∗ 100 ∗ 𝑉𝑏 79 2
30.48
Masa de combustible quemado Moles producidos

12 + 𝑛 ∗ 𝑚 + 1 ∗ 𝐶𝑂 𝑄𝑝 ∗ 𝑟𝑢𝑛 𝑡𝑖𝑚𝑒
𝑚𝑓 = ∗ 𝑚𝑜𝑙 𝑝𝑟𝑜𝑑 𝑚𝑜𝑙. 𝑝𝑟𝑜𝑑. =
100 379
Tasa de inyección y producción

𝑙𝑖𝑡𝑟𝑜 1000
𝑄𝑖 = 𝑄 ∗ ∗ 60
min⁡ 30,483
Solución de PAPER del profesor José Ramón Rodríguez
Analytical Model of the Combustion Zone in Oxygen-Enriched In Situ Combustion Tube

Experiments
Tabla de resultados obtenidos en el trabajo de investigación
Diámetro: 3pulg Qinj: 3 litro/ min
Esp. De tub: 1/16 pulg Run time: 6horas
Longitud: 91,51 cm Qp: 2,6 litro / min
Replicando los resultados antes mostrados
 Tasa de inyección y producción

𝑄𝑖 = 3 ∗ ∗ 60 = 6.356 𝑝𝑖𝑒 3 /𝑕𝑜𝑟𝑎
min⁡ 30,483
𝑄𝑝 = 2,773 ∗ ∗ 60 = 5.8756 𝑝𝑖𝑒 3 /𝑕𝑜𝑟𝑎
min⁡ 30,483
 Área y radio de la celda
2
3 2 ∗ 1/16
𝐴=𝜋∗ − = 0.0412 𝑝𝑖𝑒 2
24 12
3 2 ∗ 1/16
𝑟= − = 0.1145 𝑝𝑖𝑒
24 12
 Calculo de la razón de las moléculas de CO2 a moléculas de CO (m)

Para la corrida 1 (run 1)
%𝐶𝑂2 12.41
𝑚= = → 𝑚 = 3,0268
%𝐶𝑂 4.1
Se repite este paso con el resto de corridas
 Calculo de la razón de átomos de hidrogeno a átomos de carbono (n)
𝑂2𝑖
4∗ ∗ 𝑁2 − 𝐶𝑂2 + 0.5 ∗ 𝐶𝑂 + 𝑂2
1 − 𝑂2𝑖
𝐹𝐻𝐶 =
𝐶𝑂2 + 𝐶𝑂
0.21
4∗ ∗ 0.8122 − 0.1241 + 0.5 ∗ 0.041 + 0.0229
𝐹𝐻𝐶 = 1 − 0.21
0.1241 + 0.041
𝐹𝐻𝐶 = 1.173
 Calculo de la Eficiencia de Utilización de Oxigeno
𝑂2𝑖 𝑂2
𝐸𝑜2 = 1 − ∗
1 − 𝑂2𝑖 𝑁2
0.21 0.0229
𝐸𝑜2 = 1 − ∗
1 − 0.21 0.8122
𝐸𝑜2 = 0.894
 Calculo de los Moles producidos
5.8756 ∗ 6
𝑚𝑜𝑙. 𝑝𝑟𝑜𝑑. = = 0.093 𝑙𝑏𝑚 /𝑚𝑜𝑙
379
 Calculo de masa de combustible quemado
12 + 1.173 ∗ 3.0268 + 1 ∗ 4.1 𝑙𝑏𝑚
𝑚𝑓 = ∗ 0.093 = 0.24796
100 𝑐𝑜𝑚𝑏
 Calculo del Contenido de combustible o Fc
0.24796
𝐶𝑚 = = 2.04 𝑙𝑏𝑚/𝑝𝑖𝑒 2
89.92
0.0412 ∗
30.48
 Calculo del flujo de aire
𝑄𝑖 6.356 𝑝𝑖𝑒 3 /𝑕𝑜𝑟𝑎 𝑃𝐶𝑁

𝑈𝑎 = = = 154.296
𝐴 0.0412 𝑝𝑖𝑒 2 𝐻𝑜𝑟𝑎 ∗ 𝑝𝑖𝑒 2
 Calculo del requerimiento de aire para quemar el combustible depositado (α)
379 2𝑚 + 1 𝐹𝐻𝐶 𝐹𝑐
𝛼= ∗ + ∗
𝐸𝑜2 ∗ 𝑂2𝑖 2𝑚 + 2 4 12 + 𝐹𝐻𝐶
379 2 ∗ 3.027 + 1 1.173 2.04

𝛼= ∗ + ∗
0.894 ∗ 0.21 2 ∗ 3.027 + 2 4 12 + 1.173
𝑃𝐶𝑁
𝛼 = 365.489
𝑓𝑡3 𝑑𝑒 𝑟𝑜𝑐𝑎
Se repite este paso para el resto de corridas
 Calculo de la velocidad del frente de combustión (Vf)
𝑈𝑎
𝑉𝑓 =
𝛼
154.296 𝑝𝑖𝑒
𝑉𝑓 = = 0.422 = 12,86 𝑐𝑚/𝑕𝑟
365.489 𝑕𝑜𝑟𝑎
Se repite este paso para el resto de corridas
Replicando los resultados antes mostrados
 Tasa de inyección y producción

𝑄𝑖 = 3 ∗ ∗ 60 = 6.356 𝑝𝑖𝑒 3 /𝑕𝑜𝑟𝑎
min⁡ 30,483
𝑄𝑝 = 2,6 ∗ ∗ 60 = 5.509 𝑝𝑖𝑒 3 /𝑕𝑜𝑟𝑎
min⁡ 30,483

Apuntes de Crudos Pesados UDO PDF

Cargado por

Copyright:

Formatos disponibles

Apuntes de Crudos Pesados UDO PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Apuntes de Crudos Pesados UDO PDF

Cargado por

Copyright:

Formatos disponibles

Exploración de Crudos Pesados PREPARADOR: JUAN ZAMORA

Guía de práctica de yacimientos de

Inyección cíclica de vapor

 Mecanismos en la recuperación cíclica de vapor

Inyección continua de vapor

 Mecanismos en la recuperación continua de vapor

CAPITULO II. PROPIEDADES TÉRMICAS DE ROCAS Y FLUIDOS

Viscosidad del Petróleo

µ: Viscosidad del petróleo, cp

T*: Temperatura en grados absolutos, R= (°F + 460) K= (°C +273,1)

o Carta ASTM de viscosidad – temperatura

𝜌𝑜 : Densidad a la temperatura deseada pero en T(°F), gr/cc

𝜌𝑜𝑠𝑐 : Densidad a condiciones normales, gr/cc

𝐿𝑜𝑔 Log 𝜈 + 0.7 = 𝑛𝐿𝑜𝑔 𝑇 + 𝑏

Viscosidad de las mezclas liquidas

Viscosidad del Agua y del Vapor

µw: Viscosidad del agua, cp

Ts: Temperatura de Saturación, °F

Para vapor seco y saturado: Para vapor húmedo:

𝜇𝑠 = 0.2 ∗ 𝑇𝑠 + 81.97 𝑋10−4 𝜇𝑤𝑠 = 𝑋𝑆𝑡 ∗ 𝜇𝑠 + 1 − 𝑋𝑆𝑡 ∗ 𝜇𝑤

µs: Viscosidad del vapor seco saturado, cp

µws: Viscosidad del vapor húmedo y saturado, cp

Xst: Calidad del vapor en fracción

Efecto del gas en solución sobre la viscosidad del petróleo

Otras Unidades de Viscosidad

o Segundos Redwood N°1,t*:

Densidad del petróleo

𝛾𝑜 : Gravedad especifica del petróleo adimensional

𝜌𝑜 : Densidad a la temperatura deseada pero en T(°F), gr/cc

𝜌𝑜𝑠𝑐 : Densidad a condiciones normales, gr/cc

Densidad del Agua:

𝐺 = −6.6 + 0.0325 ∗ 𝑇 + 0.000657 ∗ 𝑇 2

𝜌𝑤 : Densidad del agua a T(°F), lb/pie3

Hidrocarburos líquidos Hidrocarburos gaseosos

Capacidad calorífica de las rocas

𝜌𝑛 : Densidad a la temperatura deseada pero en T(°F), Lb/pie3

Conductividad Térmica de líquidos y gases

Fracciones de Petróleo y mezclas de Gases

Vapor a altas Temperaturas

𝐾𝑣 = 0.5778 ∗ 10−4 ∗ 176 + 0.587 ∗ 𝑇 + 1.04 ∗ 10−3 ∗ 𝑇 2 − 4.51 ∗ 10−7 ∗ 𝑇 3

Conductividad Térmica de las Rocas

Considerando el efecto de la densidad, Considerando el efecto de la porosidad,

6.36 6.36 ∗ 𝑒 0.6∗ 2.65∗ 1−∅ +𝑆𝑙

𝜌𝑟 : Densidad de la matriz o de la roca, Lb/pie3 𝜌𝑟 : Densidad de la matriz o de la roca, Lb/pie3

Saturación Residual remanente de Petróleo en la zona de vapor

Al calcular viscosidades con datos obtenidos en ecuación 1:

 : Vis cos idad Absoluta (cps )

Se determina la densidad del petróleo

𝜇@120 = 𝜈@120 ∗ 𝜌𝑜@120 = 632.8436 ∗ 0.8799 = 556.9 𝑐𝑝𝑠

Temperatura Temperatura Viscosidad Densidad Viscosidad Viscosidad Viscosidad

60 520 1700 0,9042 1880 1700

120 580 0,8799 0,06768066 632,898562 556,9

200 660 180 0,8450 213 180

400 860 0,7688 0,31181871 38,7390011 29,8

o Calculo de la Capacidad Calorífica M

o Calculo de la Conductividad Térmica

6.36 ∗ 𝑒 0.6∗ 2.65∗ 1−0.25 +1 𝐵𝑇𝑈

o Calculo de la Difusividad Térmica

CAPITULO III. PROPIEDADES TÉRMICAS DEL AGUA Y DEL VAPOR