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]