Razón de Gas Disuelto-Petróleo - Rs Ejemplo

 

Ahora un pequeño ejemplo de la entrada anterior, este es un ejercicio propuesto, el cual tiene como datos una presión de burbuja de 2100 psi, una temperatura de yacimiento de 200 °F, con una densidad API de 36 y una densidad del gas de 0.95.

 import numpy as np

 Pb = 2100

T = 200

API = 36

Yg = 0.95

 #Correlación Standing, M.B. para p = 4000 psi (p>pb), epor lo que Rs = Rsb

 

Rs1  = Yg*(((Pb/18.2)+1.4)*10**(0.0125*API-0.000917*T))**1.2048

 

#Correlación de Lasater, J.A.

 

Mo = 630 - 10*API

 

Yo = 141.5/(API+131.5)

 

G = (Pb*Yg)/(T+460)

 

if G < 3.29:

 

                Yg1  = 0.359*np.log(((1.473*Yg*Pb)/(T+460))-0.476)

 

elif G >= 3.29:

 

                Yg1 = (((0.121*Pb*Yg)/(T+460))-0.236)**0.281

 

else:

 

                print("Error")

 

Rs2 = (132755*Yo*Yg1)/(Mo*(1-Yg1))

 

#Correlación de Vázquez, M.E. y Beggs, H.D.

 

C1 = 0.0178; C2 = 1.187; C3 = 23.931

 

Rs3  = (C1*Yg*Pb**C2)*np.exp((C3*API)/(T+460))

 

#Correlación de Glaso, O.

 

F = (10)**((2.8869-(14.1811-3.3093*np.log10(Pb))**0.5))

 

Rs4 = Yg*(F*((API**0.989)/(T**0.172)))**1.2255

 

#Correlación de la TOTAL, C.F.P

 

C1 = 15.0057; C2 = 0.0152 ; C3 = 4.484*10**-4 ; C4 = 1.095

 

Rs5  = Yg*((Pb/C1)*10**((C2*API-C3*T)))**C4 

 

#Correlación de Al-Marhoun, M.A.

 

Rs6  = (185.84321*Pb*Yg**1.87784*Yo**(-3.1437)*(T+460)**(-1.32657))**1.3984

 

#Correlación de Dokla, M.E y Osman , M.E.

 

Rs7  = (0.11956*10**(-3)*Pb*Yg**1.01049*Yo**-0.107991*(T+460)**0.952584 )**1.3811

 

#Correlación de Petrosky, G.E., Jr. y Farshad, F.F.

 

Rs8  = (Yg**0.8439*((Pb/112.727)+12.34)*10**((7.916*10**-4*API**1.5410-4.561*10**-5*T**1.3911)))**1.73184

 

#Correlación de Kartoatmodjo, T. y Schmidt, Z.

 

C1 = 0.0315; C2 = 0.7587; C3 = 11.2895; C4 = 0.9143

 

Rs9 = C1*(Yg**C2)*(Pb**(1/C4))*(10**((C3*API)/(T+460)))

 

 

print("""El resumen de resultados es :

 

                Standing, M.B ---------------------------- """+str(round(Rs1,4))+""" ft^3/Bl

                Lasater, J.A. ---------------------------- """+str(round(Rs2,4))+""" ft^3/Bl

                Vásquez, M.E. y Beggs, H.D. -------------- """+str(round(Rs3,4))+""" ft^3/Bl

                Glaso, O. -------------------------------- """+str(round(Rs4,4))+"""  ft^3/Bl

                TOTAL, C.F.P. ---------------------------- """+str(round(Rs5,4))+""" ft^3/Bl

                Al-Marhoun, M.A. ------------------------- """+str(round(Rs6,4))+""" ft^3/Bl

                Dokla, M.E. y Osman, M.E. ---------------- """+str(round(Rs7,4))+""" ft^3/Bl

                Petrosky, G.E., Jr. y Farshad, F.F. ------ """+str(round(Rs8,4))+""" ft^3/Bl

                Kartoatmojo, T. y Schmidt, Z. ------------ """+str(round(Rs9,4))+""" ft^3/Bl""")

 

El programa imprime lo siguiente:

 

El resumen de resultados es:

 

                Standing, M.B ---------------------------- 616.1827 ft^3/Bl

                Lasater, J.A. ---------------------------- 408.0557 ft^3/Bl

                Vásquez, M.E. y Beggs, H.D. -------------- 547.6642 ft^3/Bl

                Glaso, O. -------------------------------- 535.9464  ft^3/Bl

                TOTAL, C.F.P. ---------------------------- 673.8101 ft^3/Bl

                Al-Marhoun, M.A. ------------------------- 711.3068 ft^3/Bl

                Dokla, M.E. y Osman, M.E. ---------------- 725.0022 ft^3/Bl

                Petrosky, G.E., Jr. y Farshad, F.F. ------ 584.8041 ft^3/Bl

                Kartoatmojo, T. y Schmidt, Z. ------------ 538.0479 ft^3/Bl