Compresibilidad del Aceite por Debajo de la Presión de Burbuja Ejemplo

En esta entrada ejemplificaremos lo visto en la anterior entrada, recuerde que la entrada anterior es sobre la compresibilidad del aceite por debajo de la presión de burbuja. En esta ocasión solo se desarrolla el ejemplo de la compresibilidad, en entradas siguientes se subirá un ejemplo completo obteniendo todas las propiedades necesarias para llevar a cabo el ejemplo.

 

import numpy as np
 
P = 2000    #psi
Pb = 2500   #psi
T = 180     #°F
γAPI = 31   #°API
γg = 0.95
γo = 0.871
Psc = 628.10 # psi
Tsc = 434.44 #°R
Z = 0.762
Bg = 0.00123
Bo = 1.32170
 
Rs = 516 # PCN
 
# Correlación de McCain, W.D., Jr. y Cols
 
A = -7.573-1.450*np.log(P)-0.383*np.log(Pb)+1.402*np.log(T+460)+0.256*np.log(γAPI)+0.449*np.log(Rs)
 
Co1 = np.exp(A)
 
# Si no se cuenta con la presión de burbuja
 
B  = -7.663-1.497*np.log(P)+1.115*np.log(T+460)+0.533*np.log(γAPI)+0.184*np.log(Rs)
 
Co2 = np.exp(B)
 
# Si no se cuenta con la presión de burbuja ni con la relación de subilidad.
 
C = -7.114-1.394*np.log(P)+0.981*np.log(T+460)+0.770*np.log(γAPI)+0.446*np.log(γg)
 
Co3 = np.exp(C)
 
 
# Correlación Standing, M.B.
 
D = 14.4*(10**-5)*(γg/γo)**0.5*(Rs*(γg/γo)**0.5+(1.25*T))**0.2
 
D1 = Rs/(0.83*P+21.15)
 
Co4 = -(1/Bo)*(D*D1-Bg*D1)
 
 
# Correlación Vázquez, M.E. y Beggs, H.D.
 
C1 = 4.67*10**-4; C2 = 1.187; C3 = 1.337*10**-9
 
E = C1+C3*(T-60)*(γAPI/γg)
 
E1 = C2*(Rs/P)
 
Co5 = -(1/Bo)*(E*E1-Bg*E1)
 
 
# Correlación Glaso, O.
 
F = Rs*(γg/γo)**0.526+0.968*T
 
F = (2.91329-0.55366*np.log10(F))*10**(-6.58511+2.91329*np.log10(F)-0.27683*np.log10(F)**2)*(1/F)*(γg/γo)**0.526
 
F1 = 2.02777*(Rs/(P*(14.1811-3.3093*np.log10(P))**0.5 ))
 
Co6 = -(1/Bo)*(F*F1-Bg*F1)
 
 
# Correlación Total, C.F.P.
 
C4 = 1.095
 
G =4.857*10**-4+(17.569*10**-9)*(T-60)*(γAPI/γg)
 
G1  =C4*(Rs/P)
 
Co7 = -(1/Bo)*(F*F1-Bg*F1)
 
 
 
# Correlación Al-Marthoun, M.A.
 
F = (Rs**0.74239)*(γg**0.323294)*(γo**-1.20204)
 
H =(1.35556*(10**-3)*F+4.723*(10**-6)*F**2)/Rs
 
H1 =1.3984*(Rs/P)
 
Co8 = -(1/Bo)*(H*H1-Bg*H1)
 
 
# Dokla, M.E. y Osman, M.E.
 
F = (Rs**0.773572)*(γg**0.40402)*(γo**-0.882605)
 
I =(1.08126*(10**-3)*F+5.887*(10**-6)*F**2)/Rs
 
I1 =1.3811*(Rs/P)
 
Co9 = -(1/Bo)*(I*I1-Bg*I1)
 
 
# Petroskky, G.E., Jr y Farshad, F.F.
 
F = Rs**0.3738*((γg**0.2914)/(γo**0.6265))+0.24626*(T**0.5371)
 
J = 83.313*(10**-6)*(F**2.0936)*(Rs**-0.6262)*((γg**0.2914)/(γo**0.6265))
 
J1 =1.73184*(Rs/(P+139.05))
 
Co10 = -(1/Bo)*(J*J1-Bg*J1)
 
 
# Correlación artoatmodjo, T. Schmidt, Z.
 
C4 = 0.9143
 
F = (Rs**0.755)*(γg**0.25)*(γo**-1.5)+0.45*T
 
K = 1.1325*(10**-4)*(((F**1.5)-0.45*T*(F**0.5))/Rs)
 
K1 = (1/C4)*(Rs/P)
 
Co11 = -(1/Bo)*(K*K1-Bg*K1)
 
 
print("""
 
El resumen de resultados es:
 
Correlaciones                                                                                                                    Co - 1/psi
 
Correlación de McCain, W.D., Jr. y Cols ............ """+str(round(Co1,6))+"""
Correlación de McCain, W.D., Jr. y Cols ............ """+str(round(Co2,6))+"""
Correlación de McCain, W.D., Jr. y Cols ............ """+str(round(Co3,6))+"""
Correlación Standing, M.B. ......................... """+str(round(Co4,6))+"""
Correlación Vázquez, M.E. y Beggs, H.D. ............ """+str(round(Co5,6))+"""
Correlación Glaso, O. .............................. """+str(round(Co6,6))+"""
Correlación Total, C.F.P. .......................... """+str(round(Co7,6))+"""
Correlación Al-Marthoun, M.A. ...................... """+str(round(Co8,6))+"""
Dokla, M.E. y Osman, M.E. .......................... """+str(round(Co9,6))+"""
Petroskky, G.E., Jr y Farshad, F.F. ................ """+str(round(Co10,6))+"""
Correlación artoatmodjo, T. Schmidt, Z. ............ """+str(round(Co11,6)))
 
El código devuelve  lo siguiente:
 
El resumen de resultados es:
 
            Correlaciones                                                                                                                    Co - 1/psi
 
            Correlación de McCain, W.D., Jr. y Cols ............ 0.000144
            Correlación de McCain, W.D., Jr. y Cols ............ 0.000142
            Correlación de McCain, W.D., Jr. y Cols ............ 0.000159
            Correlación Standing, M.B. ......................... 0.000154
            Correlación Vázquez, M.E. y Beggs, H.D. ............ 0.000176
            Correlación Glaso, O. .............................. 0.000142
            Correlación Total, C.F.P. .......................... 0.000142
            Correlación Al-Marthoun, M.A. ...................... 0.000214
            Dokla, M.E. y Osman, M.E. .......................... 0.000194
            Petroskky, G.E., Jr y Farshad, F.F. ................ 0.000222
            Correlación artoatmodjo, T. Schmidt, Z. ............ 0.000169