'CORRELACIONES PARA ESTIMAR z'

'METODO DE PAPAY'Function Z_papay(SGg, P_psi, T_R) P_pc = 756.8 - (131.8 * SGg) - (3.6 * SGg ^ 2) Psr = P_psi / P_pc T_pc = 169.2 + (349.5 * SGg) - (74# * SGg ^ 2) Tsr = T_R / T_pc Z_papay = 1 - ((3.52 * Psr) / (10 ^ (0.9813 * Tsr))) + ((0.274 * (Psr ^ 2)) / (10 ^ (0.8157 * Tsr)))End Function

Function T_R(T_F) T_R = T_F + 459.67End Function

'METODO DE DRANCHUCK-ABOU-KASSEM'Function Z_Dranchuk(Gsg, P_psi, T_R)Static T2, T3, T4, T5, rho, Zl, F, FPConst A1 = 0.3265, A2 = -1.07, A3 = -0.5339, A4 = 0.01569Const A5 = -0.05165, A6 = 0.5475, A7 = -0.7361, A8 = 0.1844Const A9 = 0.1056, A10 = 0.6134, A11 = 0.721 P_pc = 756.8 - (131.8 * Gsg) - (3.6 * Gsg ^ 2) Ppr = P_psi / P_pc T_pc = 169.2 + (349.5 * Gsg) - (74# * Gsg ^ 2) Tpr = T_R / T_pcZ_Dranchuk = 1If Ppr > 0 Then T2 = A1 * Tpr + A2 + A3 / (Tpr ^ 2) + A4 / (Tpr ^ 3) _ + A5 / (Tpr ^ 4) T3 = A6 * Tpr + A7 + A8 / Tpr T4 = -A9 * (A7 + A8 / Tpr) T5 = A10 / (Tpr ^ 2) rho = 0.27 * Ppr / Tpr Do Zl = Z_Dranchuk F = rho * (Tpr + T2 * rho + T3 * rho ^ 2 + T4 * rho ^ 5 _ + T5 * rho ^ 2 * (1# + A11 * rho ^ 2) _ * Exp(-A11 * rho ^ 2)) - 0.27 * Ppr FP = Tpr + 2 * T2 * rho + 3 * T3 * rho ^ 2 _ + 6 * T4 * rho ^ 5 + T5 * rho ^ 2 _ * Exp(-A11 * rho ^ 2) * (3 + A11 * rho ^ 2 _ * (3 - 2 * A11 * rho ^ 2)) rho = rho - F / FP Z_Dranchuk = 0.27 * Ppr / Tpr / rho Loop While Abs(Z_Dranchuk - Zl) >= 0.00005End IfEnd Function

'METODO DE BRILL AND BEGGS'Function Z_BRILL(SGg, P_psi, T_R) P_pc = 756.8 - (131.8 * SGg) - (3.6 * SGg ^ 2) Psr = P_psi / P_pc T_pc = 169.2 + (349.5 * SGg) - (74# * SGg ^ 2) Tsr = T_R / T_pc A = (1.39 * (Tsr - 0.92) ^ 0.5) - 0.36 * Tsr - 0.1 B = ((0.62 - 0.23 * Tsr) * Psr) + (((0.066 / (Tsr - 0.86)) - 0.037) * Psr ^ 2) + ((0.32 / (10 ^ (9 * (Tsr - 1)))) * Psr ^ 6) C = 0.132 - (0.32 * (Log(Tsr) / Log(10))) D = 10 ^ (0.3106 - (0.49 * Tsr) + (0.1824 * (Tsr ^ 2))) Z_BRILL = A + ((1 - A) / Exp(B)) + (C * (Psr ^ D))End Function

'Rs CORRELACION DE STANDING'Function Rs_standing(SGg, P_psia, T_F, API, R_sb)F = ((R_sb / SGg) ^ 0.83) * (10 ^ (0.00091 * T_F - 0.0125 * API))Pb = 18.2 * (F - 1.4)Select Case P_psia Case Is < Pb Rs_standing = SGg * ((P_psia / 18.2 + 1.4) * 10 ^ (0.0125 * API - 0.00091 * T_F)) ^ 1.2048 Case Is >= Pb Rs_standing = SGg * ((Pb / 18.2 + 1.4) * 10 ^ (0.0125 * API - 0.00091 * T_F)) ^ 1.2048End SelectEnd Function

'Bo CORRELACION DE SATANDING'Function Bo_Standing(R_sb, Ge_gas, P_psi, T_F, API, Co) n = ((R_sb / Ge_gas) ^ 0.83) * (10 ^ (0.00091 * T_F - 0.0125 * API)) Pb = 18.2 * (n - 1.4) Ge_oil = 141.5 / (API + 131.5)Select Case P_psi Case Is Pb F = (R_sb * (Ge_gas / Ge_oil) ^ (1 / 2)) + 1.25 * T_F Bob = 0.9759 + (12 * (10 ^ -5) * (F ^ 1.2)) Bo_Standing = Bob * Exp(Co * (Pb - P_psi))End SelectEnd Function

'uo Correlacin de Vazquez, M.E. y Beggs, H.D.'Function u_o_VazquezyBeggs(T_F, API, R_sb, SGg, P_psi) F = ((R_sb / SGg) ^ 0.83) * (10 ^ (0.00091 * T_F - 0.0125 * API)) Pb = 18.2 * (F - 1.4) Select Case P_psi Case Is = 15 Z = 3.0324 - 0.02023 * API y = 10 ^ Z X = y * T_F ^ -1.163 u_od = 10 ^ X - 1 u_o_VazquezyBeggs = u_od Case Is Pb Z = 3.0324 - 0.02023 * API y = 10 ^ Z X = y * T_F ^ -1.163 u_od = 10 ^ X - 1 A = 10.715 * (R_sb + 100) ^ -0.515 B = 5.44 * (R_sb + 150) ^ -0.338 u_ob = A * (u_od) ^ B M = 2.6 * P_psi ^ 1.187 * Exp(-11.513 - 8.98 * 10 ^ -5 * P_psi) u_o_VazquezyBeggs = u_ob * (P_psi / Pb) ^ M End Select'funcion valida para: Presin, lpca. 141 - 9515'Razn gas disuelto-petrleo, PCN/BN 9.3 - 2199'Gravedad del petrleo, "API 15.37 - 59.5'Gravedad especifica del gas, (aire = 1) 0.511,- 1.351'Viscosidad, cp. 0.1178- 148End Function

'uo Correlacin de Beggs-Robinson'Function u_o_Begg_Robinson(T_F, API, R_sb, SGg, P_psi) F = ((R_sb / SGg) ^ 0.83) * (10 ^ (0.00091 * T_F - 0.0125 * API)) Pb = 18.2 * (F - 1.4) Select Case P_psi Case Is = 15 Z = 3.0324 - 0.02023 * API y = 10 ^ Z X = y * T_F ^ -1.163 u_od = 10 ^ X - 1 u_o_Begg_Robinson = u_od Case Is