MODELING OF FORMATION OF ACID WATER DISCHARGED FROM HIGH TEMPERATURE GEOTHERMAL RESERVOIR

Similar documents
Geochemical monitoring of the response ofgeothermal reservoirs to production load examples from Krafla, Iceland

Geochemical Modeling of Acidic Geothermal Fluids using SOLVEQ and CHIM-xpt

Geochemical Characteristics of Reservoir Fluid from NW-Sabalan Geothermal Field, Iran

predictive mineral discovery*cooperative Research Centre A legacy for mineral exploration science Mineral Systems Q3 Fluid reservoirs

FLUID-MINERAL EQUILIBRIA AND SUBSURFACE TEMPERATURE EVALUATION OF SELECTED HOT SPRINGS, JIANGXI PROVINCE, SE-CHINA

CHEMICAL GEOTHERMOMETERS FOR GEOTHERMAL EXPLORATION

Geochemical Modelling of Low-Temperature Geothermal Fields from Bihor County, Romania

Prediction of Calcite Scaling at the Oguni Geothermal Field, Japan: Chemical Modeling Approach

MULTICOMPONENT EQUILIBRIUM MODELS FOR TESTING GEOTHERMOMETRY APPROACHES

Graduate School of Science and Engineering, Kagoshima University

Fluid Geochemistry at the Nir Geothermal Field, Nw-Iran

Key words: seismic activity, permeability, heat source, Kakkonda Granite, Iwate Volcano, Japan 2. SEISMIC ACTIVITY IN HIGH PERMEABLE ZONE

Enhancement of Silica-Enthalpy Mixing Model to Predict Enthalpy of Geothermal Reservoir

REDUCING COST AND ENVIRONMENTAL IMPACT OF GEOTHERMAL POWER THROUGH MODELING OF CHEMICAL PROCESSES IN THE RESERVOIR

1) In-flow of very hot fluid of more than 300T at the

WAMUNYU EDWARD MUREITHI I13/2358/2007

APPLICATION OF GEOCHEMICAL METHODS IN GEOTHERMAL EXPLORATION

APPLICATION OF GEOCHEMICAL METHODS IN GEOTHERMAL EXPLORATION

THE IMPLICATION OF 34 S AND 18 O ISOTOPE SYSTEMATICS ON THE SULFUR CYCLE IN MAHANAGDONG HYDROTHERMAL SYSTEM (PHILIPPINESS)

Geochemical Monitoring of the Lakes District Area

Capabilities of TOUGH Codes for Modeling Geologic Sequestration and Leakage of CO 2

Effective geochemical methods for identifying geothermal systems in the western branch of the EARS

THE FLUID CHARACTERISTICS OF THREE EXPLORATION WELLS DRILLED AT OLKARIA-DOMES FIELD, KENYA

Geol. 656 Isotope Geochemistry

WATER-ROCK INTERACTION IN THE BAKKI LOW-TEMPERATURE GEOTHERMAL FIELD, SW-ICELAND

THE CHEMICAL CHARACTERISTICS OF GEOTHERMAL FLUIDS IN JIAODONG PENINSULA, SHANDONG, CHINA

APPLICATION OF GEOCHEMICAL METHODS IN GEOTHERMAL EXPLORATION. Halldór Ármannsson November 2007

GEOCHEMISTRY OF RWENZORI HOT SPRINGS. Vincent Kato Department of Geological Survey and Mines, Entebbe, Uganda

N. Todaka 1, C. Akasaka 1, T. Xu 2, and K. Pruess 2. Proceedings World Geothermal Congress 2005 Antalya, Turkey, April 2005

4 Carbonates and CO2

Takigami geothermal system, northeastern Kyushu, Japan

ALTERATION ZONATION OF SILICA MINERALS IN A GEOTHERMAL SYSTEM A NUMERICAL SIMULATION BASED ON REACTION-TRANSPORT MODEL

SUBSURFACE HYDROTHERMAL ALTERATION AT THE ULUBELU GEOTHERMAL FIELD, LAMPUNG, SOUTHERN SUMATRA, INDONESIA. Suharno 1, 2 and PRL Browne 2

PERMEABLE FRACTURES AND IN-SITU STRESS AT THE KAKKONDA GEOTHERMAL FIELD, JAPAN

Effect of chemical composition to large scale CO 2 Injection in Morrow Sandstone, Farnsworth Hydrocarbon Field, Texas, USA

DISSOLUTION, TRANSPORT AND PRECIPITATION OF SILICA IN GEOTHERMAL SYSTEM

Geochemical Modeling of Acid Fluids in Los Humeros Geothermal Field, Mexico

Model Study of Formation of the Cap Rocks for Geothermal System Using ChemTOUGH2

HYDROTHERMAL FLUID FLOW AND MINERAL ALTERATION IN A FRACTURED ROCK UNDER MULTIPHASE H 2O- CO 2 MIXTURE CONDITIONS

Luigi Marini and Claudio Pasqua. ELC-Electroconsult S.p.A. Via 1 Maggio 41, I Baranzate, Milano, Italy

TWO DIFFERENT ROLES OF FRACTURES IN GEOTHERMAL DEVELOPMENT

GEOCHEMICAL CHARACTERISTICS OF ACID FLUIDS IN MT PINATUBO, PHILIPPINES

Metals in geothermal fluids and metal precipitations in surface pipelines in Iceland can these been utilized?

Isotope and Gas Geochemistry of Dieng Geothermal Field, Indonesia

Geochemical and Hydrological Assessment of Thermal Features in the Vicinity of the Champagne Pool in the Waiotapu Geothermal System

Redox, ph, pe OUTLINE 9/12/17. Equilibrium? Finish last lecture Mineral stability Aquatic chemistry oxidation and reduction: redox

Create custom rock (Rock1) and fluid (Fluid1) compositions. 1. Copy the folder Module3 to your project folder located in Library\Gems3\projects.

Update on the Stable Isotope Geochemistry of the Leyte Geothermal Production Field

Lecture 11: Petrology of Mineralizing Aqueous Solutions GY303 Igneous & Metamorphic Petrology

Preliminary Assessment of the Ohinewai Low Enthalpy Geothermal System, New Zealand

GEOCHEMICAL INTERPRETATION OF THERMAL FLUID DISCHARGE FROM WELLS AND SPRINGS IN BERLÍN GEOTHERMAL FIELD, EL SALVADOR

Geochemistry & mineralogy of late-metamorphic shear zones:

The Initial-State Geochemistry as a Baseline for Geochemical Monitoring at Ulubelu Geothermal Field, Indonesia

Geothermometer, Geoindicator and Isotope Monitoring in Lahendong Wells during

Numerical Simulation Study of the Mori Geothermal Field, Japan

THREE DIMENSIONAL MODELING OF GEOELECTRICAL STRUCTURE BASED ON MT AND TDEM DATA IN MORI GEOTHERMAL FIELD, HOKKAIDO, JAPAN

HYDROTHERMAL ALTERATION IN THE SUNAGOHARA FORMATION, OKUAIZU GEOTHERMAL SYSTEM, JAPAN

Geochemistry: water-rock interaction, water and gas geochemistry, isotopes, geothermometry

Hydrothermal Chemistry/ Reverse Weathering. Marine Chemistry Seminar

THE MINERALOGIC CONSEQUENCES AND BEHAVIOR OF DESCENDING ACID- SULFATE WATERS: AN EXAMPLE FROM THE KARAHA-TELAGA BODAS GEOTHERMAL SYSTEM, INDONESIA

Geochemical and Isotopic Characteristics of Geothermal Water in the Gonghe Basin of the Northeast Tibetan Plateau, China

40 th Anniversary Chemical Characteristics of Geothermal Fluids in Jiaodong Peninsula, Shandong, China. Tingting Zheng Student UNU GTP

Insights to High Temperature Geothermal Systems (New Zealand) from Trace Metal Aqueous Chemistry

David I. Norman and Joseph N. Moore

GEOCHEMISTRY CHANGE DURING CIRCULATION TEST OF EGS SYSTEM

Kizito Maloba Opondo. Kenya Electricity Generating Company

JosC Tenorio Mejia', Franco D'Amore' and Jane Y. Gerardo3

Active and fossil geothermal systems in Iceland

Simulation of Mineral Precipitation and Dissolution in the 5-km Deep Enhanced Geothermal Reservoir at Soultz-sous-Forêts, France

organisms CaCO 3 + H 2 O + CO 2 shallow water

WATER-ROCK INTERACTION IN A LOW-ENTHALPY BACK-RIFT GEOTHERMAL SYSTEM, NEW ZEALAND

Lecture 5. Introduction to Stable Isotopes

Lawrence Berkeley National Laboratory Lawrence Berkeley National Laboratory

Heat (& Mass) Transfer. conceptual models of heat transfer. large scale controls on fluid movement. distribution of vapor-saturated conditions

Shuichi HATTORI Director, 1st Construction Division, Japan Railway Construction, Transport and Technology Agency

THE GEOTHERMAL POTENTIAL OF BILIRAN ISLAND, PHILIPPINES

HIGH TEMPERATURE HYDROTHERMAL ALTERATION IN ACTIVE GEOTHERMAL SYSTEMS A CASE STUDY OF OLKARIA DOMES

INTERPRETATION OF CHEMICAL COMPOSITION OF GEOTHERMAL FLUID FROM THE GEOTHERMAL FIELD OF BARANSKY VOLCANO, ITURUP ISLAND, RUSSIA

Geochemistry and Isotope Study of Discharged Geothermal Fluids, NW Sabalan Geothermal Field, NW Iran

CO 2 -water-rock reactivity at hydrothermal temperatures: The BigRig2 experiment

Carbonates and Silica Scaling Model Applied to Hot Water Pipes for Sanitary Use at Hotel Río Perlas, Costa Rica

GAS EQUILIBRIA CONTROLLING H 2 S IN DIFFERENT PHILIPPINE GEOTHERMAL FIELDS

Activity and Concentration

Japan Engineering Consultants, Inc., Energy and Industrial Technology Development Organization,Tokyo, Japan

Numerical Modelling in Predictive Mineral Discovery: Geochemical Models

GEOTHERMAL AND HOT SPRING WATER ORIGIN DETERMINATION USING OXYGEN AND HYDROGEN STABLE ISOTOPE IN THE TOYOHIRAKAWA CATCHMENT, HOKKAIDO, JAPAN

Various Geoscientific Investigations of Low-Enthalpy Geothermal Sites in the United Arab Emirates

(4) Give an example of important reactions that are responsible for the composition of river water.

A New Method of Evaluation of Chemical Geothermometers for Calculating Reservoir Temperatures from Thermal Springs in Nevada

THE FORMATION OF THE PANGUNA PORPHYRY COPPER DEPOSIT, with an appendix on the. Frieda porphyry copper prospect, New Guinea.

EQUILIBRIUM VERSUS KINETICS IN C02 DOMINATED HYDROTHERMAL ALTERATION F. MAY

Sian Nichols*, Honours New Celebration. Kate Moran, Honours

UNDERSTANDING THE RESISTIVITIES OBSERVED IN GEOTHERMAL SYSTEMS

FLUID STRATIGRAPHY OF THE COSO GEOTHERMAL RESERVOIR

Search and Discovery Article #50999 (2014)** Posted August 18, Abstract

2. California Energy Company, Inc. 3. Thermochem, Inc.

Types of Metamorphism!

INTERPRETATION OF INTERFERENCE EFFECTS IN THREE PRODUCTION WELLS IN THE KAWERAU GEOTHERMAL FIELD, NEW ZEALAND. Lynell Stevens and Kevin J Koorey

APPLICATION OF GEOCHEMICAL METHODS IN GEOTHERMAL EXPLORATION

GEOCHEMICAL INTERPRETATION OF THE MASAYA-GRANADA-NANDAIME CHEMICAL DATA, NICARAGUA

Transcription:

MODELING OF FORMATION OF ACID WATER DISCHARGED FROM HIGH TEMPERATURE GEOTHERMAL RESERVOIR Kohei Akaku 1, Kaichiro Kasai 2 and Katsuto Nakatsuka 3 1 JAPEX Research Center, 11 Hamada, Mihamaku, Chibashi 2610025, Japan 2 Japan Metals and Chemicals Co., Ltd. 72 Sasamori, Ukai, Takizawamura, Iwategun, Iwateken 0200172, Japan 3 Department of Geoscience and Technology, Graduate School of Engineering, Tohoku University Aoba, Aramaki, Aobaku, Sendaishi 9808579, Japan Key words: acid water, chemical equilibrium, waterrock reaction, boiling in high temperature geothermal reservoir, Fushime geothermal system, Kakkonda geothermal system ABSTRACT The Ferich acid geothermal waters (ph = 35) associated with large enthalpy fluids are discharged from the high temperature reservoirs (over 300 o C) in the Fushime (Yamagawa) and Kakkonda geothermal systems. In contrast, the moderate temperature reservoirs in these geothermal systems produce neutral to slightly alkaline waters with a little or no excess enthalpies. The geothermal water in Fushime is derived from seawater, while that in Kakkonda is believed to be of meteoric origin. In order to study the formation processes of these saline and dilute acid waters, numerical modeling based on chemical equilibrium for liquidsolid phases was carried out. The models involve boiling in the reservoir caused by productioninduced pressure decrease near the well and consequent reactions with wall rocks. The calculated fluid compositions through the models closely agree with the compositions of the fluids discharged from the high temperature reservoirs in Fushime and Kakkonda. The calculation results indicate that the H + supply from sphalerite precipitation during boiling is the most probable cause of the observed acidity in the saline Fushime waters. On the other hand, the H + supply from the reactions with pyrite and magnetite and from redox reaction possibly contributes to the observed acidity in the dilute Kakkonda waters. The models demonstrate that productioninduced boiling and consequent chemical reactions in the reservoirs are able to produce the observed acidity, without any contribution of acidic volatiles such as HCl and SO 2. 1. INTRODUCTION As the geothermal energy exploitation reached to the high temperature and deep reservoirs, acid waters of ph 3 to 5 have been found in the Fushime (Yamagawa) and Kakkonda geothermal field, Japan (Akaku et al., 1991; Yanagiya et al., 1996). The geothermal water in the Fushime system is derived from seawater whereas that in Kakkonda is believed to be of meteoric origin. The higher temperature reservoirs in these fields (over 300 o C) produce acid waters containing remarkable Fe2+. Local nearwell boiling caused by productioninduced pressure decrease is common in these reservoirs. On the other hand, the moderate temperature reservoirs in these fields produce neutral to slightly alkaline waters containing much lower Fe 2+. The single liquid phase feed or a little excess enthalpy discharge is usual. Since the acid waters sometimes introduce engineering problems related to corrosion and scale formation, it is important to know the source of acidity of these waters. This paper aims to report model calculations to investigate the effects of productioninduced boiling in reservoir and consequent reaction processes with wall rocks. 2. OUTLINE OF THE GEOTHERMAL SYSTEMS 2.1 The Fushime geothermal system The Fushime geothermal field is located close to the southeastern coast of the Satsuma Peninsula in Kyushu. The temperatures measured in some wells exceed 350 o C at depths below 1,500 m. The wells penetrate thick dacitic to andesitic tuffs and lavas. The most productive reservoir is situated around the intrusive dacite found at depths below 1,600 m (e.g., Okada et al., 2000). The ph of the waters discharged from wells is in the range of 3.8 8.0 at room temperature, after steam separation. The acid waters (ph < 5) are discharged from higher temperature wells above 300 o C, with large excess enthalpy. The temperatures at which the waters are in equilibrium with Nafeldspar and Kfeldspar (330 o 350 o C) are higher than the quartz geothermometer temperatures (270 o 325 o C). The waters are almost exactly saturated with anhydrite at the quartz geothermometer temperatures (Akaku et al., 1991). It is possible that the slow readjustment of Na + /K +, relative to reequilibration with quartz and anhydrite, in response to decreasing temperatures by productioninduced reservoir boiling is a cause of the discrepancy. In contrast, the neutral waters are discharged from the lower temperature wells (usually less than 300 o C), with small excess enthalpies. As the chloride concentrations, Br/Cl ratios and δd in the high temperature reservoir waters are close to those of seawater, it is believed that the geothermal fluid originates from seawater. Various combinations of boiling of the seawater reacted with rocks at 350 o C and dilution by groundwater are believed to occur in the geothermal system (Akaku et al., 1991, Okada et al., 2000). The chemical composition of the waters is characterized by depletion in Mg 2+ and, whereas enrichment in K +, Ca 2+, Fe 2+, Mn 2+, Zn 2+, Pb 2+, SiO 2, etc. over those of seawater. This characteristic is common in both acid and neutral waters. However, it is noticeable that the acid waters contain much larger Fe 2+ than the neutral waters. The principal alteration minerals observed in the production zone 895

are quartz, albite, Kfeldspar, sericite, chlorite, epidote and anhydrite. An abundance of Kfeldspar is the most distinctive feature around the feed points of the wells. There is no remarkable difference in alteration mineral assemblage between the high temperature (over 300 o C) and the moderate temperature production zones (Yagi and Kai, 1990). 2.2 The Kakkonda geothermal system The Kakkonda geothermal field is located inland in the northeastern part of Japan. The temperature profiles of the geothermal system have a steep increase at depths around 1,500 m, suggesting that two different convection cells exist. The upper and lower convection cells are called the shallow and deep reservoir respectively. The shallow reservoir in Miocene dacitic tuffs and shales, whose temperature is from 230 o to 260 o C, is very permeable. The deep reservoir in PreTertiary formations and a neogranitic pluton (the Kakkonda granite), whose temperature ranges from 300 o to 350 o C, is less permeable. Singlephase liquid feed is common in the shallow reservoir, while excess enthalpy conditions (twophase feed) commonly develop during production in the deep reservoir (Hanano and Takanohashi, 1993). Yanagiya et al. (1996) reported that the waters from the shallow reservoir are slightly alkaline (ph = 89, measured at room temperature, after steam separation), but those from the deep reservoir are acidic (ph = 35). The acid waters contain remarkable Fe 2+ than the neutral waters. The Na/K geothermometry (e.g., Giggenbach, 1988) suggests that the deep waters are in equilibrium with Na/Kfeldspars at about 300 o C, while most of deep waters are supersaturated with respect to quartz at 300 o C. The origin of the both waters is probably meteoric based on the isotopic compositions. The chloride concentrations in the shallow reservoir waters were about 500 mg/kg before they were disturbed by reinjection of waste water. The deep reservoir waters have about 1,000 mg/kg of chloride, although there are some uncertainties caused by excess enthalpy discharge conditions (Yanagiya et al., 1996). Yanagiya et al. (1996) interpreted that the acid waters from the deep reservoir have not fully attained to the equilibrium with the minerals in wall rocks. A possible sustained supply of magmatic volatiles such as HCl from a heat source is considered. The immature deep waters may be neutralized by the reactions with wall rocks during ascent. The deepest well in Kakkonda, WD1a, drilled by the New Energy and Industrial Technology Development Organization (NEDO) penetrated the Quaternary granite, a heat source of the geothermal system. The maximum temperature recorded in the well was over 500 o C (e.g., Doi et al., 1998). Emission of CO 2 and S was experienced during drilling. Although the hypersaline NaFeKMnCaCl type brine was obtained from the bottom of the well after the circulation of river water, its ph was 5.2. The sustained supply of magmatic volatiles from a heat source was not confirmed (e.g., Kasai et al., 1998). The secondary minerals in the Kakkonda system were formed through both of the contact metamorphism by the Quaternary granite and the hydrothermal alteration. Abundant and common hydrothermal minerals in the deep reservoir are quartz, albite, chlorite, muscovite, anhydrite and pyrite. Kfeldspar, magnetite, pyrrhotite, hematite and actinolitetremolite are also observed as hydrothermal alteration minerals in the deep reservoir (Doi et al., 1998). Acid alteration minerals such as alunite, sulfur, kaolinite and pyrophyllite are not observed in the deep reservoir and heat source, also arguing against the sustained supply of acidic volatiles. 3. MODEL CALCULATIONS 3.1 Calculation procedures Numerical modeling is performed by using CHILLER (Reed, 1982; Spycher and Reed, 1992), a multicomponent chemical equilibrium calculation program. For simulating effects of productioninduced boiling in reservoir, we envision the following processes. The liquid water fully equilibrated with the alteration minerals in wall rocks begins to boil owing to productioninduced pressure decrease near the well. The boiling in reservoir sometimes involves isothermal process due to heat transfer from wall rocks to the fluid. In response to the changes of physical conditions by boiling, the reactions of the fluid with the alteration minerals proceed to establish a new equilibrium until it enters the well. However, some reactions are too slow to attain equilibrium. In the well and surface equipment, boiling and consequent cooling are continued. They result in precipitation of scale. The chloride concentration of the primary reservoir water (singlephase liquid) for the Fushime case is set to 19,000 ppm, being equal to that in seawater. For the Kakkonda case, it is set to 1,000 ppm. The CO 2 concentrations of the waters are set to 1,230 ppm and 250 ppm, based on the measured total discharge concentrations for the Fushime well SKG16 and for the Kakkkonda well19 (Akaku et al., 1991, Yanagiya et al., 1996). The compositions of the full equilibrium waters are computed by adding the principal alteration and scale minerals observed in the fields into these waters at 350 o C, until the waters are saturated with them. The minerals used in the calculations are quartz, albite, Kfeldspar, muscovite, clinochlore (Mg endmember of chlorite), anhydrite, clinozoisite (Ca endmember of epidote), pyrite, magnetite, sphalerite and galena. The isoenthalpic reservoir boiling calculations from 350 o through 250 o C are conducted for the Fushime and Kakkonda case. Precipitation and dissolution of quartz, muscovite, clinochlore, anhydrite, pyrite, magnetite, sphalerite and galena are allowed, while albite, Kfeldspar and clinozoisite are suppressed in the calculations. The isothermal reservoir boiling calculation with a stepwise increment of enthalpy at 300 o C is tested only for the Kakkonda case. The liquid water fully equilibrated with the minerals mentioned above at 300 o C is used as a starting solution. The precipitation and dissolution of clinochlore, anhydrite, pyrite, magnetite, sphalerite and galena are only allowed in the isothermal calculation. The redox reaction 896

for is allowed in both of the isoenthalpic and isothermal models. The fluids computed at 300 o C in the isoenthalpic models boil to 100 o C. Precipitation of scale minerals (sphalerite and galena, commonly observed in these fields) is only allowed, while other minerals are suppressed in the calculations. The redox reaction for is also disallowed for kinetic reasons. The same type of the calculations for the fluids computed at the enthalpies of 1,680 kj/kg and 2,057 kj/kg in the isothermal model are also carried out. The ph at 20 o C is calculated by cooling of the boiled waters (separated from phase) at 100 o C. Precipitation of any minerals is disallowed in the cooling calculations. 3.2 Results for the Fushime system At the beginning of boiling, CO 2, S and strongly partition into phase (Fig. 1a). The S desing drives dissolution of pyrite from wall rock and precipitation of magnetite (Fig. 1b). An increase in aqueous Fe 2+ concentration is caused by the dissolution of pyrite (Fig. 1d), but most of Fe 2+ dissolved from pyrite enters precipitating magnetite (reaction 1). At 270 o C, pyrite turns to precipitate, whereas magnetite dissolves. An increase in concentration is caused by the combined effects of pyrite dissolution and redox reaction (reaction 2) at the early stage of boiling. However, dissolution of anhydrite supplies at the later stage (Fig. 1b and 1d). The continuous precipitation of sphalerite (reaction 3) and quartz is driven by temperature decrease. A small amount of galena also precipitates. The aqueous Mg 2+ concentration increases with the dissolution of clinochlore, while the Al 3+ concentration decreases with the precipitation of muscovite (Fig. 1b and 1d). The water equilibrated with the principal alteration and scale minerals has a ph of 5.6 at 350 o C (Fig. 1c). The ph decreases to 5.1 at 300 o C and to 5.3 at 250 o C, as boiling and reactions proceed ( in Fig. 1c). The ph decrease results from the excess supply of H + from the reactions 1, 2 and 3 over the H + consumption by CO 2 desing (reaction 4). Therefore, the separate calculation in which all these reactions are disallowed shows the ph increase with decreasing temperatures ( w/o mineral and redox reac. in Fig. 1c). If sphalerite and galena are suppressed, a slight decrease in ph at high temperatures and a steep increase at low temperatures are caused by the reactions with pyrite and magnetite ( w/o sph & gal in Fig. 1c). It is clearly demonstrated that the H + supply from sphalerite precipitation is the most effective cause of the ph decrease in the model. 3FeS 2 + 6 O Fe 3 O 4 + 0.5 + 5.5 S + H + (1). (pyrite) (aq) (magnetite) (aq) (aq) (aq) O + 0.25 S + 0.25 + 0.5H + (2). (aq) (aq) () (aq) (aq) ZnCl 2 + S ZnS + 2H + + 2Cl (3). (aq) (aq) (sphalerite) (aq) (aq) + H + CO 2 + O (4). (aq) (aq) () (aq) The calculated ph at 20 o C is 3.3 when precipitation of sphalerite and galena during boiling from 300 o C to 100 o C is allowed (Fig. 2). If their precipitation is disallowed, a higher ph of 4.3 is obtained. These values are fairly close to the measured ph for the waters discharged from the Fushime wells (ph = 3.84.7). The calculated fluid composition also closely agrees with the observed compositions, except for the minor discrepancies in Fe 2+ and Pb 2+ concentrations. 3.3 Results for the Kakkonda system The calculation results of isoenthalpic reservoir boiling for the Kakkonda case are shown in Fig. 3. The quantity of sphalerite precipitation from dilute Kakkonda type water is remarkably smaller than that from saline Fushime type water (Fig. 3b). Desing of S from aqueous phase drives dissolution of pyrite and precipitation of magnetite (reaction 1) from 350o to 290 o C and its reverse reaction proceeds at low temperatures. The ph change during boiling is also different from the Fushime case (Fig. 3c). The ph decreases very steeply at higher temperatures, due to the combined effects of reactions 1 and 2 (from 6.4 at 350 o C to 5.3 at 320 o C). The reverse progress of reaction 1 causes the ph increase at lower temperatures ( in Fig. 3c). Therefore, the separate calculation in which the reactions with pyrite and magnetite are disallowed causes a continuous ph decrease as boiling proceeds, owing to the H + supply from redox reaction ( w/o Fe minerals in Fig. 3c). It is obvious that the ph of the dilute Kakkonda type water is scarcely influenced by the precipitation of sphalerite and galena because the quantity of their precipitation is small. The calculated ph at 20 o C is on the range of the measured ph for the waters discharged from the Kakkonda wells (ph = 3. 4.5). However, the calculated Fe 2+ concentration is far apart from the observed concentrations (Fig. 4). The calculated concentrations of and Al 3+ are also lower than the observed concentrations. It is indicated that the isoenthalpic model cannot fully explain the processes occurred in the Kakkonda deep reservoir. One possibility to explain the observed high concentrations of Fe 2+ may be that isothermal boiling process in the reservoir leads to more vigorous desing of S from the liquid phase to cause the further dissolution of pyrite. Isothermal boiling results in continuous dissolution of pyrite and precipitation of magnetite from enthalpies of 1,345 kj/kg (the enthalpy of the liquid water at 300 o C) through 2,057kJ/kg (Fig. 5b). The concentrations of S and increase with increasing enthalpies owing to the progress of reaction 1. The ph decrease is mainly due to the reactions with pyrite and magnetite, though some contribution of redox reaction can be seen in Fig. 5c. The concentrations of Fe 2+, Mg 2+, Al 3+, Zn 2+ and Pb 2+ increase owing to the dissolution of pyrite, clinochlore, sphalerite and galena. A small amount of anhydrite precipitation scarcely influences on the Ca 2+ concentration (Fig. 5b and 5d). The calculated ph at 20 o C for the fluid computed at the enthalpy 897

of 1,680 kj/kg in the isothermal model is 3.4, and it is on the range of the measured ph (from 3.2 to 4.5). The calculated composition closely agrees with the observed compositions, but the calculated Fe 2+ concentration is still lower than the observed concentrations (Fig. 4). The fluid computed at the enthalpy of 2,057 kj/kg contains much higher Fe 2+, being similar to the observed concentrations. However, the calculated ph at 20 o C (2.9) is lower than the measured ph. In addition, the small disagreements with respect to Mg 2+ and Al 3+ concentrations can be seen in Fig. 4. 4. DISCUSSION The fluid composition calculated through the isoenthalpic model for the Fushime type water is very similar to those of the fluids discharged from the wells, except for the minor disagreements in Fe 2+ and Pb 2+. The reason for the disagreements is not clear, but it is possible that the isothermal boiling process reconciles the disagreement in Fe 2+ as we describe in the Kakkonda case. It is also possible that the equilibrium with magnetite and galena is not fully achieved in the real system. Another possibility is that the waters discharged from wells may be enriched in Fe 2+ by dissolution from steel equipment. The calculation results indicate that the H + supply from sphalerite precipitation is the most probable cause to form acid waters in the Fushime system, as already demonstrated by Akaku et al., (1991, 1997). The H + supply from dissolution of pyrite and precipitation of magnetite occurred at high temperatures is also able to contribute to the acidity, but it is less effective than the sphalerite precipitation in the isoenthalpic boiling model. The isoenthalpic model for the Kakkonda case forms acid water only at high temperatures, and fails to explain the observed high concentrations of Fe +. In contrast, the isothermal boiling model shows the continuous ph decrease and the Fe + increase with increasing enthalpies. The calculated composition through the isothermal model fit the observed compositions more closely. This may suggest that the isothermal model is more appropriate for the Kakkonda case. The both models indicate that dissolution of pyrite and precipitation of magnetite driven by S desing is effective to form acid water. In addition, the H + supply from redox reaction is able to cause the ph decrease. The possibility that acidic volatiles such as HCl and SO 2 contribute to form acid water discharged from the Fushime and Kakkonda wells was discussed in Reed (1991), Yanagiya et al. (1996) and Akaku et al., (1997). However, the contribution of acidic volatiles is not necessary to explain the observed acidity in our models. 5. CONCLUSIONS The model calculations indicate that temperature decrease and desing caused by productioninduced boiling in reservoir drive the reactions of fluid with minerals in wall rocks. The calculated compositions through the models closely agree with those of the fluids discharged from the high temperature reservoirs in the Fushime and Kakkonda geothermal systems. The models are able to explain the observed acidity, without taking account of acidic volatiles such as HCl and SO 2. The H + supply from sphalerite precipitation is the most probable cause of the observed acidity in the saline Fushime waters. On the other hand, the H + supply from dissolution of pyrite and precipitation of magnetite and from the redox reaction possibly contributes to the observed acidity in the dilute Kakkonda waters. ACKNOWLEDGMENTS The authors thank M. H. Reed and Y. Yamada for helpful discussions. REFERENCES Akaku, K., Reed, M. H., Yagi, M., Kai, K. and Yasuda, Y. (1991). Chemical and physical processes occurring in the Fushime geothermal system, Kyushu, Japan. Geochem. Jour. Vol.25 (4), pp.315333. Akaku, K., Kasai, K., Nakatsuka, K. and Uchida, T. (1997). The source of acidity in water discharged from high temperature geothermal reservoirs in Japan. Proc. 22 nd Workshop on Geothermal Reservoir Engineering. Stanford University, Stanford. pp.427434. Doi, N., Kato, O., Ikeuchi, K., Komatsu, R., Miyazaki, S., Akaku, K. and Uchida, T. (1998). Genesis of the plutonichydrothermal system around Quaternary granite in the Kakkonda geothermal system, Japan. Geothermics, Vol.27 (5/6), pp.663690. Giggenbach, W. F. (1988). Geothermal solute equilibria. Derivation of NaKMgCa geoindicators. Geochim. Cosmochim. Acta, Vol.52 (12), pp.27492765. Hanano, M. and Takanohashi, M. (1993). Review of recent development of the Kakkonda deep reservoir, Japan. Proc. 18 th Workshop on Geothermal Reservoir Engineering. Stanford Univ., Stanford. pp.2934. Kasai, K., Sakagawa, Y., Komatsu, R., Sasaki, M., Akaku, K. and Uchida, T. (1998). The origin of hypersaline liquid in the Quaternary Kakkonda granite, sampled from WD1a well, Kakkonda geothermal system, Japan. Geothermics, Vol.27 (5/6), pp.631645. Okada, H., Yasuda, Y., Yagi, M and Kai, K. (2000). Geology and fluid chemistry of the Fushime geothermal field, Kyushu, Japan. Geothermics, Vol.29 (2), pp.279311. Reed, M. H. (1982). Calculation of multicomponent chemical equilibria and reaction processes in systems involving minerals, es and aqueous phase. Geochim. Cosmochim. Acta, Vol.46 (4), pp.513528. Reed, M. H. (1991). Computer modeling of chemical processes 898

in geothermal systems: Examples of boiling, mixing and waterrock reaction. In: Application of geochemistry in geothermal reservoir development, F. D Amore (Coordinator), UNITAR/UNDP Centre on small Energy Resources, Rome, pp.275297. Spycher, N. F. and Reed, M. H. (1992). Users guide for CHILLER: A program for computing waterrock reactions, boiling, mixing and other reaction processes in aqueousmineral systems. University of Oregon, Eugene, Oregon. 64pp. Yagi, M. and Kai, K. (1990). Hydrothermal alteration in the Fushime geothermal field, Kyushu, Japan. Proc. 12th NZ Geothermal Workshop. Auckland, pp.3133. 1E+2 1E+1 2E3 0E+0 2E3 7 6 H O 2 CO 2 H S 2 gal spl Mgchl anhy py w/o mineral and H redox reac. 2 w/o mineral reactions w/o spl & gal qz mag musc (a) (b) (c) Yanagiya, S., Kasai, K., Brown, K. L., Giggenbach, W. F. (1996). Chemical characteristics of deep geothermal fluid in the Kakkonda geothermal system, Iwate Prefecture, Japan. Chinetsu, Vol. 33 (1), pp.118. 5 2E+0 Cl Na + K + Ca ++ (d) SiO 2 Mg ++ HS Zn ++ Fe ++ Pb ++ Al +++ 250 300 350 Temperature (deg. C) Figure 1. Numerical reservoir boiling of the Fushime type water. Isoenthalpic boiling proceeds from the right to left side of the graphs. The starting solution (1kg) is fully equilibrated with the alteration and scale minerals (quartz, albite, Kfeldspar, muscovite, clinochlore, anhydrite, clinozoisite, pyrite, magnetite, sphalerite and galena) at 350 o C. (a) Gas phase composition. (b) Mineral precipitates and dissolves. (c) ph. (d) Total aqueous molality of component species. Mineral abbreviations: qz, quartz; musc, muscovite; py, pyrite; Mgchl, clinochlore (Mg end member of chlorite); anhy, anhydrite; spl, sphalerite; gal, galena; mag, magnetite. 899

3E+0 FUSHIME 3E+0 isothermal boiling 2057kJ/kg KAKKONDA w/o scale ppt. measured H Na K Ca Mg Fe Al Zn Pb Cl SO4SiO2CO2H2S H2 isothermal boiling 1680kJ/kg isoenthalpic boiling measured H Na K Ca Mg Fe Al Zn Pb Cl SO4SiO2CO2H2S H2 Figure 2. Comparison between calculated and measured compositions for the Fushime case. Dashed line indicates the composition when scale precipitation (sphalerite and galena) is disallowed. Activity of H + (10 ph ) at 20 o C, aqueous total molality of Na +, K +, Ca 2+, Mg 2+, Fe 2+, Al 3+, Zn 2+, Pb 2+, Cl, and SiO 2 at 100 o C and mole fraction of CO 2, S and in phase at 100 o C. 1E+2 1E+1 (a) O CO 2 H S 2 (b) qz gal spl musc Mgchl mag 0E+0 7 6 5 Na + anhy w/o mineral and redox reac. Cl Al +++ HS Fe ++ Pb ++ Zn ++ SiO 2 K + Mg ++ w/o Fe minerals Ca ++ 250 300 350 Temperature (deg. C) py ph Figure 3 Isoenthalpic boiling of the Kakkonda type water. Boiling proceeds from the right to left side of the graphs. See caption to Fig. 1 for abbreviations. (c) (d) Figure 4. Comparison between calculated and measured compositions for the Kakkonda case. See caption to Fig. 2 for abbreviations. Dashed lines indicate the compositions when scale precipitation (sphalerite and galena) is disallowed. 1E+2 1E+1 2E3 0E+0 2E3 7 6 5 4 Cl Na + O CO 2 Ca ++ Al +++ Mg ++ Fe ++ Zn ++ py Pb ++ Mgchl S anhy 1300 1400 1500 1600 1700 1800 1900 2000 2100 Enthalpy (kj/kg) mag w/o mineral and redox reac. w/o Fe minerals (a) (b) (c) (d) SiO 2 K + HS Figure 5. Isothermal boiling of the Kakkonda type water. Boiling proceeds from the left to right side of the graphs. The starting solution (1kg) is fully equilibrated with the alteration and scale minerals at 300 o C. 900