In this work, the authors thermodynamic equilibrium calculations developed a computer program in MATLAB of vapor-liquid multicomponent phases below the language for aqueous speciation and critical temperature of wat...In this work, the authors thermodynamic equilibrium calculations developed a computer program in MATLAB of vapor-liquid multicomponent phases below the language for aqueous speciation and critical temperature of water (647 K) and pressures up to 2 kilobar with two computational routines. The first routine calculates the aqueous speciation from the equations of mass balance, charge balance and mass action, forming a system of N equations and N unknowns describing the aqueous phase with N unknown chemical species, which is solved with Newton-Raphson method using the Debye-Hiickel and Bdotactivity coefficients that solved simultaneously with Henry's Law, the Soave-Redlich-Kwong equation of state and a Poynting correction term, calculate chemical potentials and composition of the phases of the system, which are used at the Ratchford-Rice flash routine that iteratively recalculates the distribution of gases in liquid, until convergence. The chemical composition data of gas and brine phases of production wells at standard conditions for 1998 and 2007 of the Los Humeros geothermal field are modeled in the described manner, deducing that fluids inside the pipes of the reservoir are corrosive and cause the precipitation of newly formed minerals such as iron and calcium sulfates and silica, at reservoir conditions.展开更多
文摘In this work, the authors thermodynamic equilibrium calculations developed a computer program in MATLAB of vapor-liquid multicomponent phases below the language for aqueous speciation and critical temperature of water (647 K) and pressures up to 2 kilobar with two computational routines. The first routine calculates the aqueous speciation from the equations of mass balance, charge balance and mass action, forming a system of N equations and N unknowns describing the aqueous phase with N unknown chemical species, which is solved with Newton-Raphson method using the Debye-Hiickel and Bdotactivity coefficients that solved simultaneously with Henry's Law, the Soave-Redlich-Kwong equation of state and a Poynting correction term, calculate chemical potentials and composition of the phases of the system, which are used at the Ratchford-Rice flash routine that iteratively recalculates the distribution of gases in liquid, until convergence. The chemical composition data of gas and brine phases of production wells at standard conditions for 1998 and 2007 of the Los Humeros geothermal field are modeled in the described manner, deducing that fluids inside the pipes of the reservoir are corrosive and cause the precipitation of newly formed minerals such as iron and calcium sulfates and silica, at reservoir conditions.