EVALUACIÓN DE DIFERENTES ALGORITMOS DE SUSTITUCIONES SUCESIVAS PARA EL CÁLCULO DE EQUILIBRIO LÍQUIDO-VAPOR CON ECUACIONES DE ESTADO / EVALUATION OF SEVERAL SUCCESSIVE SUBSTITUTION ALGORITHM FOR GASLIQUID EQUILIBRIUM WITH EQUATIONS OF STATE

Los calculos de equilibrio liquido-vapor con ecuaciones de estado son ampliamente utilizados en la industria petroleray consumen gran cantidad de tiempo en las simulaciones numericas. Por esta razon, en el presente trabajo se estudianlos principales modelos de calculo flash basados en sustituciones sucesivas y las diferentes partes que conforman estosalgoritmos. Los analisis son realizados sobre los diferentes metodos para determinar la fraccion de vapor, correlacionespara estimar las constantes de equilibrio y las diferentes propuestas para actualizar las relaciones de equilibrio en elproceso iterativo. Los resultados de la evaluacion muestran que el metodo de Leibovici- Neoschil con la rutina de Newton-Raphson, es la opcion mas eficiente para determinar la fraccion de vapor, ya que requiere un menor numero de iteraciones y,por ende, un menor tiempo de computo. Para las sensibilidades realizadas con respecto al tipo de correlacion utilizada paraestimar las relaciones de equilibrio iniciales, se encontro que la propuesta de Wilson es la mas efectiva, ya que se obtienela convergencia de un mayor numero de casos. Del analisis de los metodos acelerados de Risnes-Jessen, Mehra y Crowe-Nishio para la actualizacion de las relaciones de equilibrio se pudo concluir que el algoritmo de Mehra tipo I es el masrapido de todas las opciones, sin llegar a converger en ningun caso a la solucion trivial. Por ultimo, se propone un nuevoalgoritmo que combina los metodos de Crowe-Nishio, Mehra tipo I y Leibovici, y con el cual se obtiene convergencia deun 100 % de los casos y el menor tiempo de computo con respecto a las metodologias bajo estudio. ABSTRACT The calculations of vapor-liquid equilibrium with equations of state are widely used in the oil industry, and they consume alot time in numerical simulations. Therefore, in this paper several flash calculation models based on successive substitutionsand the different parts that defines these algorithms are studied. Basically, the evaluations are performed with differentmethods to determine the fraction of vapor, correlations to estimate the equilibrium constants and the various proposalsfor upgrading equilibrium ratios in the iterative process. The results of the evaluation show that the method of Leibovici-Neoschil with Newton-Raphson routine is the most efficient way to determine the fraction of vapor since it requires feweriterations and time. Sensitivity analyses are performed considering the correlation used to estimate initial equilibriumrelationships. It has been found that Wilson’s proposal is the most effective since the convergence has been achieved ina greater number of cases. The analysis of the accelerated methods of Risnes-Jessen, Mehra et al. and Crowe Nishio ofupdating the equilibrium relationships could conclude that the type I algorithm of Mehra is the fastest of all the options,without having to converge in any case to the trivial solution. Finally, a new algorithm which combines the methods fromCrowe Nishio, Mehra type I and Leibovici has been proposed. With this model the convergence is achieved in 100% of thecases. Additionally, the computation time is lower than any of the other methods under study. Keywords: Flash calculation, Liquid-Vapor equilibrium, Successive Substitution, Leibovici, Rachford-Rice, Equilibriumconstant, Equation of state.