QrtzGeotherm: A revised algorithm for quartz solubility geothermometry to estimate geothermal reservoir temperature and vapor fraction with multivariate analytical uncertainty propagation ?
The quartz solubility geothermometry to calculate geothermal reservoir temperature and vapor fraction with multivariate analytical uncertainty propagation is programmed as two classes, SiO2TD and QrtzGeotherm in Visual Basic in Visual Studio 2010 (VB.NET). The class, SiO2TD calculates the total discharge concentration, SiO2TD and its uncertainty, SiO2TDErr from the analytical concentration of silica, SiO2msd and uncertainty, SiO2msdErr of separated water, sampled after N-separations of vapor and liquid. The class, QrtzGeotherm uses the following properties as input parameters: (i) HResreservoir enthalpy (kJ/kg), (ii) HResErruncertainty in the reservoir enthalpy (kJ/kg), (iii) SiO2TDtotal discharge silica concentration (ppm), (iv) SiO2TDErruncertainty in the total discharge silica concentration (ppm) (v) GeoEqnumber of quartz solubility regression equation, (vi) TempGuessa guess value of the reservoir temperature (C). The properties corresponding to the output parameters are (i) TempResreservoir temperature (K), (ii) TempResErruncertainty in the reservoir temperature (K), (iii) VaporResreservoir vapor fraction and (iv) VaporResErruncertainty in the reservoir vapor fraction. Similarly, it has a method, SiO2Eqn(EqNo, Temp) to provide the silica solubility as function of temperature corresponding to the regression equation. Four quartz solubility equations along the liquidvapor saturation curve: (i) a quadratic equation of 1/T and pressure, (ii) a linear equation relating log SiO2to the inverse of absolute temperature (T), (iii) a polynomial of T including logarithmic terms and (iv) temperature as a polynomial of SiO2including logarithmic terms are programmed. A demonstration program, QGeotherm is written VB.NET. Similarly, the applicability of classes SiO2TD and QrtzGeotherm in MS-Excel is illustrated considering Los Azufres geothermal field as an example.
qrtzgeotherm vb.net solubility regression equations quartz geothermometry steamtables