The temperature in the Valley of Geysers (Kamchatka) geothermal reservoir calculated using the feldspar Na-K-geothermometer has been steadily increasing over the past 10 years on average from 165 to 235 °C, which is close to the temperature values of a hydrothermal explosion of the steam and water mixture. For the analysis of chemical geothermometers, TOUGHREACT-simulation was used, with the help of which the previously known Na-K feldspar geothermometer was reproduced on a single-element model and new formulas were obtained for three Na-K geothermometers: zeolite, smectite, and based on volcanic glass. Data of chemical analysis for the period 1968-2018, in which the chloride ion is considered as an inert tracer of geofiltration processes, indicates that after 2007 a significant inflow of infiltration water (its mass fraction is estimated from 5 to 15 %) into the Geyser reservoir. It is assumed that the Na-K increased values of the feldspar geothermometer are not the result of the temperature increase in the Geyser reservoir, but the effect of smectite water dilution.
Numerical 3D model of Mutnovsky geothermal field (Dachny springs), which consist of 517 elements and partially takes into account double porosity, was developed in 1992-1993 using computer program TOUGH2. Calibration of the model was based on data from test yield of the wells and initial distribution of temperature and pressure in the reservoir. This model was used for techno-economic justification of power plant construction (Mutnovskaya GeoES, 2002). The model was recreated in the program PetraSim v.5.2, the calibration was carried out using additional data on production history before year 2006 and inversion iTOUGH2-EOS1 modeling. Comparison of reservoir parameters, estimated using inversion modeling, with previous parameter estimations (given in brackets) showed the following: upflow rate of heat-transfer agent in natural conditions 80.5 (54.1) kg/s, heat flux enthalpy 1430 (1390) kJ/kg, reservoir permeability 27∙10 –15 -616∙10 –15 (3∙10 –15 -90∙10 –15 ) m 2 . Inversion modeling was also used to estimate reinjection rates, inflow of meteoric water in the central part of geothermal field and compressibility of reservoir rocks.