Factor de volumétrico del agua, Bw Ejemplo

 

Un rápido ejemplo en el cual se utiliza una presión de 5000 psi, presión en el punto de burbuja de 2500 psi, a una temperatura de 200 °F, con una salinidad de 20000 ppm y la compresibilidad del agua de 3.43*10^-6 1/psi.

 

import numpy as np

 

p = 5000 # psi

pb = 2500 # psi

T = 200 # °F

S = 20000/10000 # ppm

Cw = 3.43*10**-6

 

# Las siguientes correlaciones funcionan por debajo de la presión de burbuja (p<pb)

 

# Correlación de McCain, W.D, Jr.

 

VwT= -1.0001*(10**-2)+1.33391*(10**-4)*T+5.50654*(10**-7)*(T**2)

Vwp= (-1.95301*(10**-9)*p*T)-(1.72834*(10**-13)*(p**2)*T)-(3.58922*(10**-7)*p)-(2.25341*(10**-10)*(p**2))

Bwb1 = (1+Vwp)*(1+VwT)

 

# Correlación de McCoy, R.L.

 

# Considerando agua pura libre de gas

 

A = 0.9947+5.8*(10**-6)*T+1.02*(10**-6)*(T**2)

B = -4.228*(10**-6)+1.8376*(10**-8)*T-6.77*(10**-11)*(T**2)

C = 1.3*(10**-10)-1.3855*(10**-12)*T+4.285*(10**-15)*(T**2)

 

Bwp= A+B*p+C*(p**2)

 

Bwb = 1+S*(5.51*(10**-8)*p+(5.47*(10**-6)-1.95*(10**-10)*p)*(T-60)-(3.23*(10**-8)-8.5*(10**-13)*p)*(T-60)**2)*Bwp

 

Bwb2 = Bwp*Bwb

 

# Considerando agua pura saturada con gas

 

A =0.9911+6.35*(10**-5)*T+8.5*(10**-7)*(T**2)

B = -1.093*(10**-6)-3.497*(10**-9)*T+4.57*(10**-12)*(T**2) 

C = -5.0*(10**-11)+6.429*(10**-13)*T-1.43*(10**-15)*(T**2)

 

Bwp= A+B*p+C*(p**2)

 

Bw = 1+S*(5.51*(10**-8)*p+(5.47*(10**-6)-1.95*(10**-10)*p)*(T-60)-(3.23*(10**-8)-8.5*(10**-13)*p)*(T-60)**2)*Bwp

 

Bwb3 = Bwp*Bwb

 

# Para obtener el factor volumétrico por encima de la presión de burbuja se utiliza la siguiente expresión

 

Bw = Bwb3*np.exp(Cw*(pb-p))

 

# Resumen de resultados

 

print("""

Resumen de Resultados

----------------------------------------------------------------

Por debajo del punto de burbuja

 

Correlación de McCain, W.D, Jr. ................ """+str(round(Bwb1,5))+""" B@Y/B@S

 

Correlación de McCoy, R.L.                                                                   

           

            Considerando agua pura libre de gas ........ """+str(round(Bwb2,5))+""" B@Y/B@S

            Considerando agua pura saturada con gas .... """+str(round(Bwb3,5))+""" B@Y/B@S

 

Por encima del punto de burbuja ................ """+str(round(Bw,5))+""" B@Y/B@S""")

 

El código devuelve lo siguiente:

Resumen de Resultados

----------------------------------------------------------------

Por debajo del punto de burbuja

Correlación de McCain, W.D, Jr. ................ 1.02806 B@Y/B@S

Correlación de McCoy, R.L.

Considerando agua pura libre de gas ........ 1.0217 B@Y/B@S

Considerando agua pura saturada con gas .... 1.03103 B@Y/B@S

Por encima del punto de burbuja ................ 1.02223 B@Y/B@S

[Finished in 690ms]