Viscosidad del Aceite Ejemplo

 En esta entrada, se pone en práctica lo aprendido en la entrada “Viscosidad del Aceite”.

 import numpy as np
 
# DATOS
 
YAPI = 31
T = 180   # °F
P = 4000  # PSI
Pb = 2500 # PSI
Rs = 675  # PCN/BN
 
# CRUDO MUERTO / SIN GAS
 
# Correlación de Beal, C.
 
a = 10**(0.43+(8.33/YAPI))
Mod1 = (0.32+((1.8*10**7)/(YAPI**4.53)))*(360/(T+200))**a
 
# Correlación de Beggs, H.D. y Robinson, J.R
 
z = 3.0324-0.02023*YAPI
y = 10**z
x = y*T**(-1.163)
Mod2 = 10**x-1
 
# Correlación de Glaso, O.
 
Mod3  = (3.141*10**10)*(T**-3.444)*(np.log10(YAPI))**(10.313*np.log10(T)-36.447)
 
# Correlación de Egbogad, E.O.
 
Mod = 1.8653-0.025086*YAPI-0.5644*np.log10(T)
 
Mod4 = 10**(10**Mod)-1
 
# Correlación de Kartoatmodjo, T. y Schimidt, Z.
 
Mod5 = (16.0*10**8)*(T**-2.8177)*np.log10(YAPI)**(5.7526*np.log10(T)-26.9718)
 
# CRUDO SATURADO / CON GAS
 
# Correlación de Chew, J.N. y Connally, C.A., Jr.
 
A = 10**(Rs*((2.2*10**-7)*Rs-7.4*(10**-4)))
b = (0.68/10**(8.62*10**(-5)*Rs))+(0.25/10**(1.1*10**(-3)*Rs))+(0.062/10**(3.74*10**(-3)*Rs))
 
Mob1 = A*Mod1**b
 
# Correlación de Beggs, H.D y Robinson, J.R.
 
a = 10.175*(Rs+100)**(-0.515)
b = 5.44*(Rs+150)**(-0.338)
 
Mob2 = a*Mod1**b
 
# Correlación de Kartoatmodjo, T. y Schmidt, Z.
 
b = 10**(-0.00081*Rs)
A = (0.2001+(0.8428*10**(-0.000845*Rs)))*(Mod1**(0.43+0.5165*b))
Mob3 = -0.06821+0.9824*A+40.34*(10**-5)*(A**2)
 
# CRUDO BAJOSATURADO
 
# Correlación de Beal, C.
 
Mo1 = (0.024*Mob1**1.6+0.038*Mob1**0.56)*(0.001)*(P-Pb)+Mob1
 
# Correlación de Vázquez, M.E. y Beggs, H.D.
 
m = 2.6*P**1.187*np.exp(-11.513-8.98*10**(-5)*P)
Mo2 = Mob1*(P/Pb )**m
 
# Correlación de Kartoatmodjo, T. y Schmidt, Z.
 
Mo3  = 1.00081*Mob1+1.127*(10**-3)*(P-Pb )*(-65.17*(10**-4)*(Mob1**1.8148)+0.038*(Mob1**1.59) )
 
print("""
---------------------------------------------------------------------
                                               Correlación                                          | Mod, cp| Mob, cp|  Mo,cp |   
---------------------------------------------------------------------
Beal, C. ............................... | """+str(round(Mod1,4))+""" | ------ | """                    +str(round(Mo1,4))+""" |
Beggs, H.D. y Robinson, J.R ............ | """+str(round(Mod2,4))+""" | """+str(round(Mob2,4))+""" | ------ |
Glaso, O. .............................. | """+str(round(Mod3,4))+""" | ------ | ------ |
Egbogad, E.O. .......................... | """+str(round(Mod4,4))+""" | ------ | ------ |
Kartoatmodjo, T. y Schimidt, Z. ........ | """+str(round(Mod5,4))+""" | """+str(round(Mob3,4))+""" | """+str(round(Mo3,4))+""" |
 Chew, J.N. y Connally, C.A., Jr. ...... |  ----  | """                    +str(round(Mob1,4))+"""  | ------ |
 Vázquez, M.E. y Beggs, H.D. ........... |  ----  | ------ | """                                        +str(round(Mo2,4))+""" |""")

# El código regresa lo siguiente

El resumen de resultados es:

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

Correlación | Mod, cp| Mob, cp|  Mo,cp |

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

Beal, C. ............................... | 2.6545 | ------ | 0.8158 |

Beggs, H.D. y Robinson, J.R ............ | 3.0354 | 0.5726 | ------ |

Glaso, O. .............................. | 2.7582 | ------ | ------ |

Egbogad, E.O. .......................... | 3.4954 | ------ | ------ |

Kartoatmodjo, T. y Schimidt, Z. ........ | 2.6279 | 0.6681 | 0.7794 |

 Chew, J.N. y Connally, C.A., Jr. ...... |  ----  | 0.745  | ------ |

 Vázquez, M.E. y Beggs, H.D. ........... |  ----  | ------ | 0.8751 |

[Finished in 1.9s]