Zero-Emission Geothermal Power Generation: Experimental Study on Carbonate Mineralization through CO2-Andesitic Pyroclastic Rock Interaction at Oku-Aizu Geothermal Plant
Ryosuke Oochi 1, Kengo Aoyama 2, Akira Ueda 3, *
1. Life Cycle Valuation Headquarters, Shimizu Corporation, 2-16-1 Kyobashi, Chuoku, Tokyo 104-8370, Japan
2. Oku-Aizu Geothermal Co. Ltd., 1034-1, Sunagohara, Yanaizu, Kawanuma-gun, Fukushima 969-7321, Japan
3. Graduate School of Science and Engineering for Research, University of Toyama, 3490 Gofuku, Toyama 980-8555, Japan
Academic Editor: Andres Navarro Flores
Special Issue: Geothermal Energy Exploration and Production
Received: November 27, 2019 | Accepted: February 03, 2020 | Published: February 12, 2020
Journal of Energy and Power Technology 2020, Volume 2, Issue 1, doi:10.21926/jept.2001002
Recommended citation: Oochi R, Aoyama K, Ueda A. Zero-Emission Geothermal Power Generation: Experimental Study on Carbonate Mineralization through CO 2 -Andesitic Pyroclastic Rock Interaction at Oku-Aizu Geothermal Plant. Journal of Energy and Power Technology 2020; 2(1): 002; doi:10.21926/jept.2001002.
© 2020 by the authors. This is an open access article distributed under the conditions of the Creative Commons by Attribution License, which permits unrestricted use, distribution, and reproduction in any medium or format, provided the original work is correctly cited.
In order to achieve advancement in zero-emission geothermal power generation for fixing the CO2 in the geothermal fluids to the carbonate minerals after binary power generation and reinjection of the fluid into the reservoir, batch-type experiments and geochemical simulations of CO2-andesitic pyroclastic rock interaction (tuff from Oku-Aizu geothermal area) were performed at 150 °C with and without CO2 gas for 15 days. Despite the low Ca content (0.9% wt.%) in the rock, the Ca2+ concentrations in the solutions that were reacted with CO2 were observed to initially increase rapidly and then decreased after two days. In contrast, the Ca2+ concentrations in the solutions reacted with Ar gas remained almost constant. The solutions reacted with CO2 gas became oversaturated with respect to the carbonate minerals according to saturation degree with respect to calcite. These findings along with the other previously available data demonstrated that Ca2+ in calcite and glass in the rocks were dissolved into the fluid with CO2 gas and were fixed as carbonate minerals within a few days.
CO2; mineralization; zero emission; experiment; geochemical simulation
According to a report provided by the Japanese government, by 2050, Japan’s current geothermal power generation of 530 MW from 17 localities will increase tenfold. However, the current CO2 emissions from the geothermal power plants in Japan amount to ca. 500,000 tons/y . The global, CO2 emissions from geothermal plants have been discussed worldwide by considering the Life Cycle Assessment [2,3,4,5,6,7]. The current installed capacity worldwide is approximately 11 GW as of 2010, and the CO2 emissions from all the geothermal fields amount to 120 kg/MWh (11.5 Mtons/y)  and to 48 kg/MWh (4.6 Mtons/y)  as reported by LCA for the average scenario. In general, CO2 emissions from the geothermal plants in Japan are low; however, a few Japanese geothermal power plants emit levels of CO2 equivalent to the thermal power plants . In such geothermal power plants, CO2 emissions could be reduced to zero by reinjecting all the generated steam including CO2 and H2S as well as the brine used for power generation into the reservoir (the process is referred to as zero-emission geothermal power generation). Binary power generation involves generating electricity using a heat medium such as ammonia or pentane. However, no previous study has reported reinjecting large amounts of steam containing CO2 and H2S along with brine (~1,000 tons/h) into the reservoir. Therefore, it is necessary to assess what type of reaction would occur when the geothermal fluids containing large amounts of CO2 and H2S are reduced underground at the temperatures in the range of 150–250 °C. In geothermal power plants, brine is injected naturally into the reinjection wells because these wells are generally drilled at a lower elevation compared to the production wells. This enables brine injection without the requirement of a pump. In the zero-emission system, the pressure of the brine and steam is maintained at high levels, so that it may be returned directly to the injection wells.
Xu and Pruess (2008)  examined the behavior of CO2 fluid in an enhanced geothermal system using theoretical simulation with TOUGHREACT  and demonstrated carbonate precipitation at high temperature. In Iceland, sub-surface mineral storage of carbon and sulfur at the CarbFix site at 40–80 °C was reported to have been achieved. The CarbFix approach involves accelerating the mineralization of the acid gases (CO2 and H2S) injected into the basaltic layers [10,11,12,13,14,15]. Basaltic rocks prove to be advantageous for the mineralization of CO2 and other acid gases as they are relatively reactive compared to most of the other rock types. Their dissolution releases the divalent metal cations, such as Ca2+, Mg2+, and Fe2+, as well as assists in the neutralization of the acidic waters as the released cations provide something for the dissolved carbonate to react with [14,15]. Natural investigations on basalt–CO2 interaction have been reported in the oceanic flood basalts  and in Idaho Springs in Idaho, USA . Laboratory-based experimental approaches have also been performed in other parts of the world [18,19,20,21].
In Japan, the injection of CO2 into high-temperature geothermal areas composed of granitic rocks for the immobilization of minerals has been studied [22,23]. This technology, referred to as “Georeactor”, injects the CO2 dissolved in the groundwater to allow it to chemically react with the granitic rocks in order to fix the carbonate minerals below the ground. As with all the chemical reactions, this reaction also occurs faster at higher temperatures, and the solubility of the resultant carbonate minerals is sufficiently low for their formation at higher temperatures. In the Sumikawa geothermal field, Akita Prefecture, Japan, the Ca2+, Mg2+, and CO2 dissolved in the water react with andesitic to dacitic lava and the tuff reservoir rocks, resulting in the formation of a natural carbonate layer, suggesting the possibility of the immobilization of the carbonate minerals in the geothermal area . Yoo et al. (2013)  evaluated the change in permeability of the reservoir layer during the carbonate deposition through stepwise numerical calculations and provided predictions regarding when and where the carbonate deposits might be found in the reservoirs. The authors also demonstrated that the permeability around the injection point tends to significantly reduce upon carbonate accumulation, which might result in the insufficient injection of CO2.
The Yanaizu–Nishiyama Geothermal Power Plant in the Oku–Aizu area of Japan had an electricity generation capacity of 65 MW when it began operating in the year 1995, and the capacity was increased by 30 MW from August 2017 onward. The steam and brine from the production well are separated using a separator, and the power is generated by the steam. The brine is reinjected to the underground, and at this point, most of the CO2 gas in the steam is released into the atmosphere. In zero-emission geothermal power generation, a two-phase fluid of steam and brine is subjected to heat-exchange with a heat medium such as ammonia, and the heat medium is then sent to a turbine for power generation. The heat-exchanged two-phase fluid is reinjected directly to the underground without any component being released into the atmosphere. The geothermal reservoir contains open fractures associated with steeply dipping faults. Initially, the reservoir used to be water dominated, although, after the commencement of operations, it was changed to two-phase fluid (steam and brine) with a few steam dominated areas. The maximum temperature inside the wells is 340 °C . The estimated concentration of the total dissolved components in the reservoir fluid is ~2% (wt.%), while the CO2 concentration is 1% (wt.%), which is a high salt and gas concentration compared to the other Japanese geothermal fields; however, the H2S concentration is considerably low [26,27,28]. The amount of CO2 (450–550 g/kWh) discharged from this power plant is higher than that discharged by the other geothermal power plants (15–150 g/kWh), and corresponds to the amount of CO2 discharged by the LNG (liquefied natural gas) plants (~600 g/kWh) . Oochi et al. (2014)  analyzed the carbon and oxygen isotope composition of calcite in the wells and suggested that the source of these high levels of CO2 and salt in the geothermal fluids was of volcanogenic origin.
The present study carried out the experimental simulation of the behavior of CO2 during the reinjection of brine containing a large amount of CO2 into the geothermal reservoir of the Yanaizu–Nishiyama Power Plant for zero-emission operation. The study focused on understanding the interaction of CO2 with the andesitic pyroclastic rock (tuff) at 150 °C (the temperature corresponding to the reservoir temperature around the reinjection wells) and predicting the involved chemical reaction using geochemical simulations.
2. Outline of the Geology
The Oku–Aizu geothermal area is located in the mountainous region of the south-western margin of Aizu Basin in the central part of the Fukushima Prefecture, Japan (Figure 1). The basement of this area consists of early Triassic sedimentary and granitic rocks. Mid-Miocene marine (and partly non-marine) pyroclastic and volcanic rocks overlay this basement [30,31,32]. Three major north-wests to south-east trending faults occur in this geothermal field: the Chinoikezawa, Sarukurazawa, and the Oizawa faults (Figure 1); the production wells were drilled into the former two faults. In the vicinity of the Sunagohara area, which occurs north of the production area, a Na+–Cl−-type hot spring (Nishiyama Hot Spring) originates from the intersection of two faults, the Takiyagawa and the Oizawa faults .
Figure 1 Geological map and cross-section of the Oku–Aizu geothermal field .
Two stages of hydrothermal alteration have been recognized in the host rocks of this region. The first one is the mid-Miocene submarine volcanic related and the other is currently ongoing alteration from the hydrothermal activity. Large quantities of chlorite and sericite are commonly observed in the reservoir rocks (rhyolitic pyroclastic rocks present in the lower Miocene Takiya River Formation), which are thought to have been formed as a result of the submarine volcanic activity , along with shallow clay mineralization (smectite, sericite, kaolinite, and zeolite). Anhydrite is present in the high-temperature zone (100-200 °C) of the reservoir, and small amounts of calcite are present at both shallow and deep zones [26,27,28,32,33].
3. Experimental Procedure
A 100-mL Teflon beaker, which was used as the reaction vessel, was placed in a stainless-steel vessel (SUS; HUT–100; San-Ai Science Co., Ltd.) (Figure 2a). These vessels were connected to a CO2 and/or vacuum line via a high-pressure valve. An andesitic pyroclastic rock sample was obtained from a 1,691 m deep core sample collected from well 84N–3t located in the Oku–Aizu geothermal area. The microscopic observations and the chemical composition of the rock sample are presented in Figure 2b and Table 1, respectively. The microscopic observation of the rock sample was performed at the Mitsubishi Materials Techno Co., Ltd, which revealed a volcanic clastic texture and a composition consisting of quartz, plagioclase, potassium feldspar, iron hydroxide, zircon, apatite, volcanic glass, pumice, andesite, and sandstone, along with alteration minerals (chlorite, illite, calcite, and iron hydroxide). The CaO content in the rock sample was 0.9% (wt.%) (Table 1), which is lower than that of the other Japanese andesitic rocks such as a standard rock JA–1 (5.7% wt.%) .
Table 1 Components of the volcanic tuff used in the experiments.
Figure 2 (a) Experimental device used for the laboratory experiments and (b) microscopic observations of the rock sample; (upper) open and (lower) crossed Nicolls.
According to the analysis of the carbon and oxygen isotopic composition of the calcites in the core samples, the geothermal fluids in the Oku–Aizu geothermal plant are believed to have originated from meteoric water with oxygen shift and/or fossil seawater . This CO2–rock interaction experiment was the first test to examine the behavior of CO2 in the zero-emission geothermal system at Oku-Aizu. Therefore, tap water that originated from river water at the University of Toyama was used as a representative fluid for meteoric water. The composition of the tap water used has been presented in Table 2. It was observed that the tap water contained small amounts of Ca2+ and other minor cations.
Table 2 Chemical composition of the tap water collected at the University of Toyama.
The rock sample (7 g) and water (70 mL) were placed in the Teflon beaker in the SUS reaction vessel. A valve was attached to the lid of the vessel in order to either degas it or inject the gas in the vessel. After degassing the vessel for 1 min using a vacuum pump at <1 KPa, CO2 gas was injected into the vessel using a syringe pump. At 150 °C, the CO2 pressure was controlled and maintained at 10 MPa, at which 3% (wt.%) CO2 was dissolved in the water . Ar gas (without CO2) was also injected into the vessel for comparative experiments. The reaction vessel was subsequently placed in a high-temperature reaction furnace and the reaction (1 rpm) was allowed to occur at 150 °C for a predetermined period. Subsequent to the rapid cooling of the vessel, when the temperature reached below 100 °C, the water and rock samples were collected for analysis. An aliquot of the water sample was filtered through a 0.22-µm filter for chemical analyses. The rock sample was dried in an oven at 110 °C for two days, following which it was subjected to SEM (scanning electron microscope) and XRD (X-ray diffraction) analyses.
4. Analytical Procedures
The water samples were analyzed for Na+, K+, Ca2+, and Mg2+, using an atomic absorption spectrometer (AA–6650; Shimadzu Corporation), and also for Cl and SO4 using ion chromatography (761 Compact IC; Metrohm). The contents of Mg2+, Mn2+, Al3+, and total-Fe were measured using ICP-MS (7500 series) at the Okayama University. SiO2 content was measured using the molybdenum yellow method and an absorptiometer (UV-VIS Spectrophotometer; Shimadzu Corporation). The chemical compositions of the rock samples prior to and after conducting the experiments were observed under a fluorescence X-ray analyzer (Supermini 200; Rigaku). Rock surface observations were performed using a scanning electron microscope (SEM; MiniscopeTM 3030; Hitachi).
The reaction temperature for this experiment was maintained at 150 °C, assuming the depth for the CO2-rich fluid-rock interaction to be 1,000 m. The internal pressure of the vessel at this temperature was maintained at 10 MPa. The amount of CO2 that could be dissolved in the water samples at 150 °C was calculated using the geochemical calculation software PhreeqC . These calculations were then used to calculate the weight of CO2 to be placed inside the reaction vessel in the laboratory environment, and consequently, the initial filling pressure of CO2 was set at 4 MPa. Similarly, the initial filling pressure of Ar gas was set at 7.3 MPa. The reaction periods used were 1, 2, 5, 10, and 15 days.
5. Results and Discussion
5.1 Chemical Compositions of Solutions after CO₂–Rock Reaction
Table 3 and Figure 3 present the concentrations of chemical components in the solutions after the CO2–rock reaction experiments. The concentrations of Na+, Ca2+, and SiO2 in the solutions reacted with CO2 gas were observed to increase with the increasing reaction time, whereas, the K+ and Mg2+ concentrations remained constant (Figure 3a). In contrast, the chemical components in the solutions reacted without CO2 (i.e., solutions reacted with Ar gas) exhibited an increase within the first two days and remained constant thereafter (Figure 3b). Figure 4 presents the concentrations of the chemical components in the solutions with and without CO2. It is clearly observable that the solutions reacted with CO2 contained higher amounts of Ca2+ compared to the solutions reacted with Ar (Figure 4a). Mg2+ content was observed to have increased in the solutions reacted with CO2. In addition, in one case among the solutions reacted without CO2, Mg2+ content was observed to increase within the first two days and decrease thereafter (Figure 4b). In one case, Na+ and K+ in the solutions reacted both with and without CO2 were observed to increase within the first two days and remain almost constant thereafter.
Table 3 Chemical compositions of the solutions after the reactions.
Figure 3 Chemical compositions of the solutions during the reaction with (a) CO2 and (b) Ar.
The silica concentration in the solutions reacted with CO2 increased from 96 to 272 mg/L when the duration increased from 1 to 15 days and was higher than that in the solutions reacted without CO2 (57–134 mg/L) for the same period (Figure 4e). The solubility of quartz and amorphous silica at 150 °C was determined to be 125 mg/L and 619 mg/L, respectively , which implied that silica concentration in the solution reacted without CO2 was saturated with respect to quartz. In contrast, silica concentration in the solution reacted with CO2 continued to increase the solubility level of amorphous silica. Kuroda et al. (2009)  determined the chemical compositions of solutions after reactions with plagioclase and granodiorite at 150 °C for 15 days and reported that the silica concentrations in the solutions reacted with CO2 and N2 were 263 mg/L and 132 mg/L, respectively, which are similar to the values obtained in the present study.
Figure 4 Chemical compositions of the solutions during the reaction showing the evolution of the concentrations of (a) Ca2+, (b) Na+, (c) Mg2+, (d) K+, (e) SiO2, (f) SO42-, and (g) Cl-.
The rock sample used in the present study contained anhydrite and silicate minerals (i.e., volcanic glass and plagioclase), which implied that Ca2+ could have been released from both anhydrite and silicate minerals. The Ca2+ concentrations in the solutions were observed to remain almost constant after two days of reaction (Figure 3a). The Ca2+/SO42− mole ratios in the solutions after reaction with CO2 were determined to be in the range of 2.7–3.3, higher than those in the solutions reacted without CO2 gas (0.7–1.0 ratios; Table 3). These results indicated that Ca2+ in the solution after reaction with CO2 gas could be released because of the dissolution of anhydrite as well as because of the reaction with silicate minerals. Ueda et al. (2005)  reported similar results, i.e., the Ca2+/SO42− mole ratios in the solutions reacted with CO2 gas after reaction with granitic rocks and distilled water at 150 °C were higher compared to those in the solutions reacted with Ar gas.
5.2 Comparison of Ca2+ Concentration in Solutions after Andesitic Tuff–CO2 Interaction with other Rock Types
Several studies have been conducted on rock–CO2 interaction for carbon sequestration, although the majority of these studies have investigated the temperature conditions below 100 °C. Figure 5 compares the results of the reaction test with granodiorite–CO2–water at the temperatures of 150 °C and 200 °C [22,38]. In the reaction with granodiorite at 150 °C, the behavior of Ca2+ concentration exhibited the same trend as that observed for the andesitic tuff in the present study. The Ca2+ concentration increased to 200 mg/L within five days and decreased thereafter. On the other hand, Ca2+ concentration maintained a low value (100 mg/L; Figure 5) in the test conducted at 200 °C. At both the reaction temperatures, the solutions reacted with CO2 exhibited higher Ca2+ concentrations compared to those reacted in the absence of CO2.
Luhmann et al. (2017)  conducted flow-through experiments with basalt and CO2-rich brine at 150 °C and observed that Ca2+ concentration increased at the beginning of the reaction up to ca. 80 mg/L and then decreased rapidly to <1 mg/L with time (within 30 days). The authors noted that after the reaction, the solutions were enriched in Fe (up to 86 mg/L), undersaturated with respect to calcite and dolomite, and oversaturated with respect to siderite. In the experiments conducted on the granodiorite and CO2 reaction at 200 °C, total Fe concentration in the solutions after the reaction was observed to reach a maximum value of 0.7 mg/L, while the concentration was less than 0.1 mg/L in most of the cases . In these experiments, calcite was observed to be oversaturated in the solutions. These results suggested that carbonate precipitation during the interaction between CO2-rich brine and the rocks is controlled by siderite for the basaltic rocks and by calcite for the granitic rocks.
Griffioen and Appelo (1993)  demonstrated that for the solutions containing Fe2+, and for those containing Ca2+ with a high concentration of CO2, siderite is precipitated first, and is followed by calcite. Using the oxygen isotope compositions of the carbonate minerals and the XRD analytical data of the rock samples, Ueda et al. (2001)  demonstrated that during the infiltration of CO2-rich groundwater, siderite was precipitated at a shallower depth in the Sumikawa geothermal area, while calcite and dolomite were observed at deeper zones. The authors also demonstrated that the carbonate minerals precipitated with smectite during the interaction and acted as a cap rock for the reservoir. These results indicated carbonate formation during the interaction with CO2-rich water.
Figure 5 Comparison of Ca concentrations in solutions reacted with and without CO2 for Oku–Aizu tuff (this study) and granodiorite [22,38] at 150 °C.
5.3 Saturation of Minerals in Solution after CO₂–Rock Reaction
In order to investigate the possibility of precipitation of minerals including the carbonate minerals, the saturation index (SI) for several minerals was calculated with the geochemical software PhreeqC , using the obtained chemical composition values for the solution after the reaction. PhreeqC is a geochemical code that utilizes solution equilibrium models to calculate aqueous speciation and the saturation state of mineral phases and has been applied previously to CO2–rock reactions [40,41,42,43,44].
The SI has been represented by log (Q/K), where Q represents the activity product (Ion Activity Product), i.e., the activity product of the actual ions in the mineral precipitation/dissolution reaction equation, and K represents the solubility product. If the SI value is positive, the mineral tends to be oversaturated, i.e., it tends to precipitate, and if the SI value is negative, it indicates that the mineral is in a dissolved state.
The results of the calculations using the analyses of the solutions produced by the reactions between the Oku–Aizu tuff and the CO2-rich water are presented in Table 4 and Figure 6. The carbonate minerals, such as calcite, aragonite, and dolomite, were observed to be over-saturated and undersaturated in the solutions reacted with and without CO2, respectively. Anhydrite was always undersaturated in the solutions reacted with and without CO2. Ueda et al. (2005)  indicated that when no CO2 is present, anhydrite dissolves initially, adding Ca2+ to the solution, thereby causing dolomite, calcite, and anorthite to become oversaturated. However, when CO2 is added, anorthite becomes substantially undersaturated and dissolves, adding more Ca2+ to the solution and resulting in an even more positive saturation index (SI) for calcite and dolomite.
Figure 6 depicts that after the reaction occurred for 1 to 14 days with and without CO2, the solutions were always saturated or oversaturated with calcite and dolomite. In particular, the SI for dolomite in the solutions reacted with CO2 was observed to increase from approximately 0 to 2 with the increasing reaction time. In contrast, the SI of dolomite in the solutions reacted without CO2 remained less than zero. These results demonstrated that plagioclase and volcanic glass were dissolved by the reaction with the CO2-saturated water along with the formation of calcite and dolomite.
Table 4 Saturation indices of mineral phases calculated from chemical compositions of solutions after the reaction.
Figure 6 Saturation indexes for the minerals in the solution after reaction with (a) CO₂ and (b) Ar. The saturation index (SI) was calculated with the geochemical software PhreeqC .
5.4 SEM Analyses of Rock Surface after the Reaction
In the present study, after undergoing the rock–CO2/Ar reaction test, the rock surface was observed under an SEM (Figure 7a and 7b). Rhomboid crystals were observed on the surface of the rock samples reacted for 15 days (Figure 7a), indicating a high possibility of the presence of aragonite, compared to the rock sample reacted with Ar (Figure 7b). No clear carbonate mineral precipitates were observed in the samples on the other reaction days. Kuroda et al. (2009)  observed calcite crystals on the surface after reacting plagioclase and CO2 at 150 °C, although such crystals were not observed in the granodiorite samples subjected to the same experimental condition. The plagioclase in that study exhibited 21.6% (wt.%) Ca2+ content, which was higher than that in the granodiorite (1.2% wt.%) and the andesitic tuff (0.9% wt.%) examined in the present study (Table 1). When the saturation index for the solutions after the rock–CO2 reaction was considered, it was inferred that the solution was oversaturated with respect to calcite. These results suggested that in the present study, aragonite and calcite must have been precipitated after the reaction and observed only for the sample reacted with CO2 for 15 days. In the case of the other samples, the precipitated amount might have been too small to be observed in the SEM analyses.
Figure 7 SEM images of the surface of rock samples after reaction with (a) CO2 and (b) Ar.
5.5 Simulation of the Behaviors of Chemical Components in Solution during the Rock-CO2 Interaction
In order to theoretically predict the chemical reactions occurring among the rock, water, and CO2, PhreeqC code , using the experimental conditions (such as temperature, pressure, and the composition of the solution) and a thermodynamic database provided by the Lawrence Livermore National Laboratory, National Institute of America were used. Using this database, calculations could be applied to high temperatures of up to 300 °C. Besides the reaction temperature used in the experiments of the present study (150 °C), the reaction temperatures of 50 °C, 100 °C, 200 °C, and 250 °C were also used in the simulation. Calculated changes in the concentration of each element under the conditions identical to the experimental conditions were analyzed. The simulation model was created with an emphasis on Na+, Ca2+, and SiO2. In the model, when the rock sample was assumed to be a sphere of size 1 mm, the reactive surface area was calculated to be 0.0023 m2/g. However, using this value, the eluted amounts determined for the main components (Na+, Ca2+, and SiO2) were extremely small. Therefore, in the model, the surface areas of the initial minerals (quartz, albite, anorthite, etc.) were set to values 100–1,000 times greater than the prescribed values. Figure 8 illustrates that the experimental observation data were the most approximated among the simulation results, although the initial increase in Ca2+ estimated by the simulation was not observed. In the present study, the Ca2+ concentration in the solutions after the reaction was observed to gradually increase, similar to the previously reported data for labradorite at 120 °C . The present simulation did not agree sufficiently with the experimental results, although it provided the approximate change in chemical compositions (Figure 8).
Figure 8 Simulation results for the solutions after reaction with (a) CO2 and (b) Ar.
5.6 Investigation of CO₂ Mineralization
In the present study, the experimental temperature used for the Oku–Aizu tuff and water interaction was 150 °C. When conducting the zero-emission geothermal power generation, the temperature for reinjection of brine with high CO2 concentration into the underground is not necessarily 150 °C. Therefore, the temperature was changed using the parameters obtained from the simulation results in order to predict the Ca2+ concentration change in the injected brine during the interaction with the andesitic tuff. In order to predict the optimum temperature for CO2 mineralization, several other reaction temperatures (50, 100, 150, 200, and 250 °C) were analyzed (Figure 9), and it was observed that the lower the reaction temperature, the more Ca was eluted. This may have occurred because the solubility of calcium carbonate increased as the temperature decreased, implying that the elution is simply a dissolution of the precipitated carbonates.
Figure 9 Simulation results for the Ca2+ concentrations in the solutions of Oku–Aizu tuff-water interaction with CO2 at different temperatures.
Using both field and batch experimental studies, Mito et al. (2008)  analyzed the carbonate mineral fixation of CO2 at ~50 °C through geochemical reactions in a sandstone reservoir at Nagaoka, Japan (the place where the first Japanese pilot project of geological CO2 storage was conducted). Seawater was reacted with powdered rock samples (rock/water wt. ratio was 100) with and without CO2 at 10 MPa. The Ca2+ concentration was observed to increase from 406 mg/L to 529 mg/L within two days and remained constant thereafter. The SI values for calcite in the solutions after the reactions with and without CO2 were determined to be −2.7 and −0.4, respectively. Mito et al. (2008)  demonstrated using field observation and theoretical calculations that the dissolution of plagioclase and chlorite has great potential for enhancing the neutralization of the acidified formation water after CO2 injection and that the bivalent cations Ca2+, Mg, and Fe provided by these minerals are expected to increase the net carbonate mineral precipitation. The authors estimated that ~5% and half of the injected CO2 could be fixed as carbonate minerals and aqueous carbonate ions, respectively, within ten years. The authors also demonstrated that mineral trapping could occur from the early stage of CO2 storage. In the present study, Ca2+ concentration in the solution after rock–water interaction at 50 °C was estimated to be ca. 500 mg/L (Figure 9), similar to the results observed in the previously stated study, although the solutions used in the experiments were different for both the studies (tap water in the present study and seawater in the other study) .
5.7 New Scaling Problem of Carbonate during Zero-Emission Geothermal Operation
Ground facilities for zero-emission geothermal power generation are expected to be identical to the existing binary power generation facilities. The existing binaries generate electricity through heat exchange with a part of brine or between a two-phase fluid of steam and brine and an ammonia medium, and so on. In the zero-emission geothermal power generation facilities, all the fluid is reinjected into the underground through heat exchange.
In the existing geothermal plants, silica scaling is one of the major problems encountered during the geothermal operation worldwide. Silica in the brines becomes oversaturated with respect to amorphous silica after the vapor–brine separation at the surface, in which scenario, the silica in the brine precipitates in the brine pipes, reinjection wells, and the layer surrounding the reinjection wells.
As for the carbonate scales, it is well known that calcite and/or aragonite precipitate close to the flashing (vapor–brine separation) point in the production wells due to the oversaturation of the Ca2+–HCO3- concentration. Several scale inhibitors such as polyacrylic acid soda may be used to prevent the precipitation. In contrast, the carbonate scales do not precipitate in the reinjection wells due to the emission of most of the CO2 into the atmosphere along with the vapor at the geothermal plants. When injecting a large amount of brine with high CO2 concentration through the reinjection wells for zero-emission power generation at a geothermal power plant, the carbonate minerals such as calcite and aragonite must be precipitated in the reservoir layer close to the reduction wells within a few days. At lower temperatures, the reinjected brine might be more enriched in Ca2+ because of interaction with the reservoir rocks and may be precipitated as carbonate minerals during neutralization. In such a scenario, the permeability of the formation is lowered through the carbonate precipitation [25,46]. Therefore, it is necessary to control the temperature conditions and the CO2 concentration in the fluid for precipitation in the reservoir without precipitating the carbonate minerals in the formation around the reinjection well. In the present study, it was observed that Ca2+ could have been dissolved from the carbonates in the rocks and also from the glass at 150 °C. When the CO2-rich acidic fluid flows into the geothermal reservoir, it is expected that Ca2+ would be leached from glass and plagioclase as a result of the reaction with the surrounding rocks. Since the solubility of calcium carbonate decreases with increasing temperature, carbonate minerals are assumed to undergo precipitation.
Despite the low Ca2+ content (0.9% wt.%) in the sample rock used in the experiments, the Ca2+ concentration in the solutions at 150 °C increased sharply over the several days of observation, and the solutions became oversaturated with respect to carbonate minerals. This suggested that CO2 may be stored as carbonate minerals such as aragonite and calcite. The simulation results for the rock–brine interaction revealed that sufficient Ca2+ elution is expected even in a low-temperature region at 50 °C. The results demonstrated that CO2 emissions may be reduced to zero by reinjecting all the steam including CO2 and brine used for power generation into the reservoir and that the zero-emission geothermal power generation would serve as one of the most useful techniques for environmental assessments.
The experiments on rock-CO2 fluid interaction in this study were performed at the University of Toyama, with part of the simulation work carried out at the Shimizu Corporation. The authors would like to thank members of Oku–Aizu Geothermal Co., Ltd. and Mitsui Mining and Smelting Co. Ltd. for various kinds of convenience and cooperation during the rock sampling and data analyses. K. Yamaguchi of Mitsubishi Materials Corporation is acknowledged for his advice and help on the simulation of the rock-water reaction. We would like to express our gratitude to A. Ozawa of Mitsubishi Materials Techno Co., Ltd. for the microscopic observation of rock samples used in this study. We express our deep appreciation to J. Zhang, K. Marumo, M. Kusakabe, and K. Horikawa from the University of Toyama for their useful advice and guidance. Also, we would like to thank members of the Geochemical Laboratory of the University of Toyama, especially Y. Yamada, R. Isaji, E. Hanajima, and R. Ikeda.
We would like to thank J. Jin, an academic editor and two anonymous reviewers for their critical, kind and constructive comments and polishing English which improved the manuscript.
The authors would like to thank Enago (www.enago.jp) for the English language review for our first draft.
These authors contributed equally to this work.
The authors have declared that no competing interests exist.
- Muraoka H, Uchida T, Noda T. CO2 emissions and future directions of geothermal power plants in Japan. Japan Geotherm; 2011; Kyoto, Japan.
- Aksoy N, Gok OS, Mutlu H, Kılınc G. CO2 emission from geothermal power plants in Turkey. Proceedings World Geothermal Congress; 2015; Melbourne, Australia.
- Armannsson H, Fridriksson T, Kristjansson B. CO2 emissions from geothermal power plants and natural geothermal activity in Iceland. Geothermics. 2005; 34: 286-296. [CrossRef]
- Fridleifsson IB, Bertani R, Huenges E, Lund JW, Rybach L. The possible role and contribution of geothermal energy to the mitigation of climate change. IPCC Scoping Meeting on Renewable Energy Sources; 2008 January 20-25; Luebeck, Germany. 59-80.
- Bayer P, Rybach L, Blum P, Brauchler R. Review on life cycle environmental effects of geothermal power generation. Renw Sust Energ Rev. 2013; 26: 446-463. [CrossRef]
- Tomasini-Montenegro C, Santoyo-Castelazo E, Gujba H, Romero RJ, Santoyo E. Life cycle assessment of geothermal power generation technologies: An updated review. Appl Therm Eng. 2017; 114: 1119-1136. [CrossRef]
- Parisi ML, Ferrara N, Torsello L, Basosi B. Life cycle assessment of atmospheric emission profiles of the Italian geothermal power plants. J Clean Prod. 2019; 234: 881-894. [CrossRef]
- Xu T, Pruess K, Apps J. Numerical studies of fluid-rock interactions in enhanced geothermal systems (EGS) with CO2 as working fluid. In Proceedings, Thirty-third Workshop on Geothermal Reservoir Engineering Stanford University; 2008 January 28-30; California, USA.
- Xu T, Sonnenthal EL, Spycher N, Pruess K. TOUGHREACT— A simulation program for non-isothermal multiphase reactive geochemical transport in variably saturated geologic media: Applications to geothermal injectivity and CO2 geological sequestration. Comput Geosci. 2006; 32: 145-165. [CrossRef]
- Alfredsson HA, Hardarson BS, Franzson H, Gislason SR. CO2 sequestration in basaltic rock at the Hellisheidi site in SW Iceland: Stratigraphy and chemical composition of the rocks at the injection site. Minerl Mag. 2008; 72: 1-5. [CrossRef]
- Alfredsson HA, Oelkers EH, Hardarsson BS, Franzson H, Gunnlaugsson E, Gislason SR. The geology and water chemistry of the Hellisheidi, SW-Iceland carbon storage site. Int J Green Gas Con. 2013; 12: 399-418. [CrossRef]
- Gislason SR, Oelkers EH. Carbon storage in basalt. Science. 2014; 344: 373-374. [CrossRef]
- Gunnarson I, Aradottir ES, Oelkers EH, Clark DE, PorArnarson M, Sigfusson B, et al. The rapid and cost-effective capture and subsurface mineral storage of carbon and sulfur at the CarbFix2 site. Int J Greenh Gas Con. 2018; 79: 117-126. [CrossRef]
- Snæbjörnsdottir SO, Oelkers EH, Mesfin K, Aradottir ES, Dideriksen K, Gunnarsson I, et al. The chemistry and saturation states of subsurface fluids during the in situ mineralisation of CO2 and H2S at the CarbFix site in SW-Iceland. Int J Green Gas Con. 2017; 58: 87-102. [CrossRef]
- Gislason, SR, Wolff-Boenisch D, Stefansson A, Oelkers EH, Gunnlaugsson E, Sigurdardottir H, et al. Mineral sequestration of carbon dioxide in basalt: A pre-injection overview of the CarbFix project. Int J Greenh Gas Con. 2010; 4: 537-545. [CrossRef]
- McGrail PB, Schaef HT, Ho AM, Chien YJ, Dooley JJ, Davidson CL. Potential for carbon dioxide sequestration in flood basalts. J Geophys Res-Sol Ea. 2006; 111: B12. [CrossRef]
- Maskell A, Kampman N, Chapman H, Condon DJ, Bickla M. Kinetics of CO2–fluid–rock reactions in a basalt aquifer, Soda Springs, Idaho. Appl Geochem. 2015; 61: 272-283. [CrossRef]
- Wolff-Boenisch D, Gislason SR, Oelkers EH, Putnis CV. The dissolution of natural glasses at pH 4 and 10.6 and temperatures from 25 to 74 °C. Geochim Cosmochim Ac. 2004; 68: 4843-4858. [CrossRef]
- Wolff-Boenisch D, Gislason SR, Oelkers EH. The effect of crystallinity on dissolution rates and CO2 consumption capacity of silicates. Geochim Cosmochim Ac. 2006; 70: 858-870. [CrossRef]
- Kanakiya S, Adam L, Esteban L, Rowe MC, Shane P. Dissolution and secondary mineral precipitation in basalts due to reactions with carbonic acid. J Geophys Res-Sol Ea. 2017; 122: 4312-4327. [CrossRef]
- Luhman AJ, Tutolo BM, Tan C, Moskowitz BM, Saar MO, Seyfried Jr WE. Whole rock basalt alteration from CO2-rich brine during flow-through experiments at 150 °C and 150 bar. Chem Geol. 2017; 453: 92-110. [CrossRef]
- Ueda A, Kato K, Ohsumi T, Yajima T, Ito H, Kaieda H, et al. Experimental studies of CO2-rock interaction at elevated temperatures under hydrothermal conditions. J Geochem. 2005; 39: 417-425. [CrossRef]
- Xu T, Kaieda H, Ueda A, Sugiyama K, Ozawa A, Wan Y, et al. Fluid-rock interactions in enhanced geothermal systems with CO2 as working fluid: Modelling of geochemical changes induced by CO2 Injection into the Ogachi (Japan) EGS Site. GRC Trans. 2010; 34: 497-502.
- Ueda A, Ajima S, Yamamoto M. Isotopic study of carbonate minerals from the Sumikawa geothermal area and its application to water movement. GRSJ. 2001; 23: 181-196.
- Yoo SY, Kuroda Y, Mito Y, Matsuoka T, Nakagawa M, Ozawa A, et al. A geochemical clogging model with carbonate precipitation rates under hydrothermal conditions. Appl Geochem. 2013; 30: 67-74. [CrossRef]
- Seki Y. The physical and chemical structure of the Oku-Aizu geothermal system, Japan. Geochem J. 1991; 25: 245-265. [CrossRef]
- Nitta T, Suga S, Tsukagosi S, Adachi M. Geothermal resources in the Okuaizu, Tohoku district, Japan. Chinetsu. 1987; 24: 26-56.
- Seki Y, Adachi M. Stratigraphy and hydrothermal alteration based on well data from the Oku-aizu geothermal system,Japan. Bull Geol Surv. 1997; 48: 365-412.
- Oochi R, Ueda A, Kume H. A preliminary study of CO2 source from isotopic composition of carbonate minerals in Oku-aizu geothermal area. Japan Geothermal; 2014; Horosaki, Japan.
- Kitamura S, Sugawara Y, Suzuki Y, Fujii K, Itoh O, Takahashi K. Fukushima prefecture geological survey report, geology of the miyashita region. Report of Fukushima Prefecture; 1968; Japan.
- Ministry of International Trade and Industry. Wide area survey report of 1964, Nishi-Aizu area. 1970; Japan.
- Suzuki K, Shimazu M, Shimada I, Manabe K. In “Route No. 28, Geological data of Cenozoic northeast Honshu arc”. Kitamura, N. ed. 1986; Japan.
- Hayakawa N, Suzuki S, Oda S. Diagenesis and hydrothermal alteration of the Neogene rocks in the Aizu district. Mining Geol. 1977; 27: 367-378.
- Imai N, Terashima S, Itoh S, Ando A. 1994 compilation of analytical data for minor and trace elements in seventeen GSJ geochemical reference samples, “Igneous rock series”. Geostandard Newslett. 1995; 19: 135-213. [CrossRef]
- Duan Z, Sun R. An improved model calculating CO2 solubility in pure water and aqueous NaCl solutions from 273 to 533 K and from 0 to 2000 bar. Chem Geol. 2003; 193: 257-271. [CrossRef]
- Parkhurst DL, Appelo CAJ. Users guide to PHREEQC (Version2)–A computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations. Water-Resources Investigations Report. 1999; 99: 312.
- Fournier RO. Chemical geothermometers and mixing model for geothermal systems. Geothermics. 1977; 5: 41-50. [CrossRef]
- Kuroda Y, Yamada Y, Ueda A, Matsuoka T, Yamada N. Experimental research of plagioclase(rock)-gas-water interaction at hydrothermal conditions for CO2 mineralization. J Mner Petrol Sci. 2009; 38: 111-121. [CrossRef]
- Griffioen J, Appelo A. Nature and extent of carbonate precipitation during aquifer thermal energy storage. Appl Geochem. 1993; 8: 161-176. [CrossRef]
- Kampman N, Bertier P, Busch A, Snippe J, Harrington J, Pipich V, et al. Validating reactive transport models of CO2-brine-rock reactions in caprocks using observations from a natural CO2 reservoir. Energy Procedia. 2017; 114: 4902-4916. [CrossRef]
- Krupka KM, Cantrell KJ, McGrail BP. Thermodynamic data for geochemical modelling of carbonate reactions associated with CO2 sequestration-literature review. Department of Energy under Contract DE-AC05-76RL01830, PNNL-19766; 2010; U.S. [CrossRef]
- Labus K, Tarkowski R, Wdowin M. Modelling gas–rock–water interactions in carbon dioxide storage capacity assessment: a case study of Jurassic sandstones in Poland. Int J Environ Sci. Technol. 2015; 12: 2493-2502. [CrossRef]
- Pham THV, Aagaard P, Hellevang H. On the potential for CO2 mineral storage in continental flood basalts – PHREEQC batchand 1D diffusion–reaction simulations. Geochem Trans. 2012; 13: 1-5. [CrossRef]
- Wang X, Wang X, Hu W, Wan Y, Cao J, Lv C, et al. Supercritical CO2-involved water–rock interactions at 85°C and partial pressures of 10–20MPa: Sequestration and enhanced oil recovery. Energy Explor Exploit. 2017; 35: 237-258. [CrossRef]
- Mito S, Xue Z, Ohsumi T. Case study of geochemical reactions at the Nagaoka CO2 injection site, Japan. Int J Green Gas Control. 2008; 2: 309-318. [CrossRef]
- Noiriel C, Steefel CI, Yang L, Bernard D. Effects of pore-scale precipitation on permeability and flow. Adv Water Resour. 2016; 95: 125-137. [CrossRef]