Geothermal investigations in an area of induced seismic activity, northern São Paulo State, Brazil

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Applications of finite element computer modelling and thermal infrared remote sensing to the study of geothermal anomalies

Buried heat sources can be investigated by examining thermal infrared images and comparing these with the results of theoretical models which predict the thermal anomaly a given heat source may generate. Key factors influencing surface temperature include the geometry and temperature of the heat source, the surface meteorological environment, and the thermal conductivity and anisotropy of the rock. In general, a geothermal heat flux of greater than 2% of solar insolation is required to produce a detectable thermal anomaly in a thermal infrared image. A heat source of, for example, 2-300K greater than the average surface temperature must be a t depth shallower than 50m for the detection of the anomaly in a thermal infrared image, for typical terrestrial conditions. Atmospheric factors are of critical importance. While the mean atmospheric temperature has little significance, the convection is a dominant factor, and can act to swamp the thermal signature entirely. Given a steady state heat source that produces a detectable thermal anomaly, it is possible to loosely constrain the physical properties of the heat source and surrounding rock, using the surface thermal anomaly as a basis. The success of this technique is highly dependent on the degree to which the physical properties of the host rock are known. Important parameters include the surface thermal properties and thermal conductivity of the rock. Modelling of transient thermal situations was carried out, to assess the effect of time dependant thermal fluxes. One-dimensional finite element models can be readily and accurately applied to the investigation of diurnal heat flow, as with thermal inertia models. Diurnal thermal models of environments on Earth, the Moon and Mars were carried out using finite elements and found to be consistent with published measurements. The heat flow from an injection of hot lava into a near surface lava tube was considered. While this approach was useful for study, and long term monitoring in inhospitable areas, it was found to have little hazard warning utility, as the time taken for the thermal energy to propagate to the surface in dry rock (several months) in very long. The resolution of the thermal infrared imaging system is an important factor. Presently available satellite based systems such as Landsat (resolution of 120m) are inadequate for detailed study of geothermal anomalies. Airborne systems, such as TIMS (variable resolution of 3-6m) are much more useful for discriminating small buried heat sources. Planned improvements in the resolution of satellite based systems will broaden the potential for application of the techniques developed in this thesis. It is important to note, however, that adequate spatial resolution is a necessary but not sufficient condition for successful application of these techniques.

Difficulties in distinguishing coherent from fragmental kimberlite: A case study of the Muskox pipe (Northern Slave Province, Nunavut, Canada)

Kimberlite drill core from the Muskox pipe (Northern Slave Province, Nunavut, Canada) highlights the difficulties in distinguishing coherent from fragmental kimberlite and assessing the volcanological implications of the apparent gradational contact between the two facies. Using field log data, petrography, and several methods to quantify crystal and xenolith sizes and abundances, the pipe is divided into two main facies, dark-coloured massive kimberlite (DMK) and light-coloured fragmental kimberlite (LFK). DMK is massive and homogeneous, containing country-rock lithic clasts (~ 10%) and olivine macrocrysts (~ 15%) set in a dark, typically well crystallised, interstitial medium containing abundant microphenocrysts of olivine (~ 15%), opaques and locally monticellite, all of which are enclosed by mostly serpentine. In general, LFK is also massive and structureless, containing ~ 20% country-rock lithic clasts and ~ 12% olivine macrocrysts. These framework components are supported in a matrix of serpentinized olivine microphenocrysts (10%), microlites of clinopyroxene, and phlogopite, all of which are enclosed by serpentine. The contact between DMK and LFK facies is rarely sharp, and more commonly is gradational (from 5 cm to ~ 10 m). The contact divides the pipe roughly in half and is sub-vertical with an irregular shape, locally placing DMK facies both above and below the fragmental rocks. Most features of DMK are consistent with a fragmental origin, particularly the crystal- and xenolith-rich nature (~ 55-65%), but there are some similarities with rocks described as coherent kimberlite in the literature. We discuss possible origins of gradational contacts and consider the significance for understanding the origin of the DMK facies, with an emphasis on the complications of alteration overprinting of primary textures.

Analytical solutions for movement of cold water thermal front in a heterogeneous geothermal reservoir

In the present study an analytical model has been presented to describe the transient temperature distribution and advancement of the thermal front generated due to the reinjection of heat depleted water in a heterogeneous geothermal reservoir. One dimensional heat transport equation in porous media with advection and longitudinal heat conduction has been solved analytically using Laplace transform technique in a semi infinite medium. The heterogeneity of the porous medium is expressed by the spatial variation of the flow velocity and the longitudinal effective thermal conductivity of the medium. A simpler solution is also derived afterwards neglecting the longitudinal conduction depending on the situation where the contribution to the transient heat transport phenomenon in the porous media is negligible. Solution for a homogeneous aquifer with constant values of the rock and fluid parameters is also derived with an aim to compare the results with that of the heterogeneous one. The effect of some of the parameters involved, on the transient heat transport phenomenon is assessed by observing the variation of the results with different magnitudes of those parameters. Results prove the heterogeneity of the medium, the flow velocity and the longitudinal conductivity to have great influence and porosity to have negligible effect on the transient temperature distribution. (C) 2013 Elsevier Inc. All rights reserved.

A field and reactive transport model study of arsenic in a basaltic rock aquifer

The use of geothermal energy as a source for electricity and district heating has increased over recent decades. Dissolved As can be an important constituent of the geothermal fluids brought to the Earth's surface. Here the field application of laboratory measured adsorption coefficients of aqueous As species on basaltic glass surfaces is discussed. The mobility of As species in the basaltic aquifer in the Nesjavellir geothermal system, Iceland was modelled by the one-dimensional (1D) reactive transport model PHREEQC ver. 2, constrained by a long time series of field measurements with the chemical composition of geothermal effluent fluids, pH, Eh and, occasionally, Fe- and As-dissolved species measurements. Di-, tri- and tetrathioarsenic species (As(OH)S22-, AsS₃H^2-, AsS33- and As(SH)4-) were the dominant form of dissolved As in geothermal waters exiting the power plant (2.556μM total As) but converted to some extent to arsenite (H₃AsO₃) and arsenate HAsO42- oxyanions coinciding with rapid oxidation of S2- to S2O32- and finally to SO42- during surface runoff before feeding into a basaltic lava field with a total As concentration of 0.882μM following dilution with other surface waters. A continuous 25-a data set monitoring groundwater chemistry along a cross section of warm springs on the Lake Thingvallavatn shoreline allowed calibration of the 1D model. Furthermore, a series of ground water wells located in the basaltic lava field, provided access along the line of flow of the geothermal effluent waters towards the lake. The conservative ion Cl^- moved through the basaltic lava field (4100m) in less than10a but As was retarded considerably due to surface reactions and has entered a groundwater well 850m down the flow path as arsenate in accordance to the prediction of the 1D model. The 1D model predicted a complete breakthrough of arsenate in the year 2100. In a reduced system arsenite should be retained for about 1ka. © 2011 Elsevier Ltd.

Feasibility Study of Solar Driven Underground Cooling System

In the United States the peak electrical use occurs during the summer. In addition, the building sector consumes a major portion of the annual electrical energy consumption. One of the main energy consuming components in the building sector is the Heating, Ventilation, and Air-Conditioning (HVAC) systems. This research studies the feasibility of implementing a solar driven underground cooling system that could contribute to reducing building cooling loads. The developed system consists of an Earth-to-Air Heat Exchanger (EAHE) coupled with a solar chimney that provides a natural cool draft to the test facility building at the Solar Energy Research Test Facility in Omaha, Nebraska. Two sets of tests have been conducted: a natural passively driven airflow test and a forced fan assisted airflow test. The resulting data of the tests has been analyzed to study the thermal performance of the implemented system. Results show that: The underground soil proved to be a good heat sink at a depth of 9.5ft, where its temperature fluctuates yearly in the range of (46.5°F-58.2°F). Furthermore, the coupled system during the natural airflow modes can provide good thermal comfort conditions that comply with ASHRAE standard 55-2004. It provided 0.63 tons of cooling, which almost covered the building design cooling load (0.8 tons, extreme condition). On the other hand, although the coupled system during the forced airflow mode could not comply with ASHRAE standard 55-2004, it provided 1.27 tons of cooling which is even more than the building load requirements. Moreover, the underground soil experienced thermal saturation during the forced airflow mode due to the oversized fan, which extracted much more airflow than the EAHE ability for heat dissipation and the underground soil for heat absorption. In conclusion, the coupled system proved to be a feasible cooling system, which could be further improved with a few design recommendations.

Continuous geochemical monitoring by mass-spectometer in the Campi Flegrei geothermal area. An application at Pisciarelli-Solfatara (diffuse and fumarolic gases) and at the mud gases during drilling of the CFDDP pilot hole

During this work has been developed an innovative methodology for continuous and in situ gas monitoring (24/24 h) of fumarolic and soil diffusive emissions applied to the geothermal and volcanic area of Pisciarelli near Agnano inside the Campi Flegrei caldera (CFc). In literature there are only scattered and in discrete data of the geochemical gas composition of fumarole at Campi Flegrei; it is only since the early ’80 that exist a systematic record of fumaroles with discrete sampling at Solfatara (Bocca Grande and Bocca Nuova fumaroles) and since 1999, even at the degassing areas of Pisciarelli. This type of sampling has resulted in a time series of geochemical analysis with discontinuous periods of time set (in average 2-3 measurements per month) completely inadequate for the purposes of Civil Defence in such high volcanic risk and densely populated areas. For this purpose, and to remedy this lack of data, during this study was introduced a new methodology of continuous and in situ sampling able to continuously detect data related and from its soil diffusive degassing. Due to its high sampling density (about one measurement per minute therefore producing 1440 data daily) and numerous species detected (CO2, Ar, 36Ar, CH4, He, H2S, N2, O2) allowing a good statistic record and the reconstruction of the gas composition evolution of the investigated area. This methodology is based on continuous sampling of fumaroles gases and soil degassing using an extraction line, which after undergoing a series of condensation processes of the water vapour content - better described hereinafter - is analyzed through using a quadrupole mass spectrometer

Secondary Microseisms Characterization and Green’s Function Extraction at the Larderello-Travale Geothermal Field (Italy)

Over the past ten years, the cross-correlation of long-time series of ambient seismic noise (ASN) has been widely adopted to extract the surface-wave part of the Green’s Functions (GF). This stochastic procedure relies on the assumption that ASN wave-field is diffuse and stationary. At frequencies <1Hz, the ASN is mainly composed by surface-waves, whose origin is attributed to the sea-wave climate. Consequently, marked directional properties may be observed, which call for accurate investigation about location and temporal evolution of the ASN-sources before attempting any GF retrieval. Within this general context, this thesis is aimed at a thorough investigation about feasibility and robustness of the noise-based methods toward the imaging of complex geological structures at the local (∼10-50km) scale. The study focused on the analysis of an extended (11 months) seismological data set collected at the Larderello-Travale geothermal field (Italy), an area for which the underground geological structures are well-constrained thanks to decades of geothermal exploration. Focusing on the secondary microseism band (SM;f>0.1Hz), I first investigate the spectral features and the kinematic properties of the noise wavefield using beamforming analysis, highlighting a marked variability with time and frequency. For the 0.1-0.3Hz frequency band and during Spring- Summer-time, the SMs waves propagate with high apparent velocities and from well-defined directions, likely associated with ocean-storms in the south- ern hemisphere. Conversely, at frequencies >0.3Hz the distribution of back- azimuths is more scattered, thus indicating that this frequency-band is the most appropriate for the application of stochastic techniques. For this latter frequency interval, I tested two correlation-based methods, acting in the time (NCF) and frequency (modified-SPAC) domains, respectively yielding esti- mates of the group- and phase-velocity dispersions. Velocity data provided by the two methods are markedly discordant; comparison with independent geological and geophysical constraints suggests that NCF results are more robust and reliable.

(Table 1) Results of geothermal measurements in the Black Sea 2000 and 2001 expeditions

According to the World Ocean Program in the northeastern part of the continental slope of the Black Sea geothermal, seismologic and seismic studies were carried out. An analysis of heat flow distribution allowed to distinguish a negative geothermal anomaly near the Dzhubga area, where the Russia-Turkey pipeline was being constructed. During seismological observations (August-September 1999, September 2001) more than 1200 seismic events were recorded. They proved high tectonic activity of the region under study, which stimulates gravitational sediment transport on the continental slope. The seismo-acoustic survey carried out in the area of the geothermal anomaly revealed no reflecting horizons within the sedimentary cover. This may be related to turbidite-landsliding processes. Results of modeling of the heat flow anomaly showed that it had originated approximately 1000 years ago due to a powerful landslide. This also suggests a possibility of an avalanche displacement of sedimentary masses in the area of the pipeline at present.

An isotopic analysis of geothermal brine and calcite scaling from the Blue Mountain geothermal field, Winnemucca, Nevada

Understanding, and controlling, the conditions under which calcite precipitates within geothermal energy production systems is a key step in maintaining production efficiency. In this study, I apply methods of bulk and clumped isotope thermometry to an operating geothermal energy facility in northern Nevada to see how those methods can better inform the facility owner, AltaRock Energy, Inc., about the occurrence of calcite scale in their power plant. I have taken water samples from five production wells, the combined generator effluent, shallow cold-water wells, monitoring wells, and surface water. I also collected calcite scale samples from within the production system. Water samples were analyzed for stable oxygen isotope composition (d18O). Calcite samples were analyzed for stable oxygen and carbon (d13C) composition, and clumped isotope composition (D47). With two exceptions, the water compositions are very similar, likely indicating common origin and a well-mixed hydrothermal system. The calcite samples are likewise similar to one another. Apparent temperatures calculated from d18O values of water and calcite are lower than those recorded for the system. Apparent temperatures calculated from D47 are several degrees higher than the recorded well temperatures. The lower temperatures from the bulk isotope data are consistent with temperatures that could be expected during a de-pressurization of the production system, which would cause boiling in the pipes, a reduction in system temperature, and rapid precipitation of calcite scale. However, the high apparent temperature indicated by the D47 data suggests that the calcite is depleted in clumped isotopes given the known temperature of the system, which is inconsistent with this hypothesis. This depletion could instead result from disequilibrium isotopic fractionation during the aforementioned boil events, which would make both the apparent d18O-based and D47-based temperatures unrepresentative of the actual water temperature. This research can help improve our understanding of how isotopic analyses can better inform us about the movement of water through geothermal systems of the past and how it now moves through modern systems. Increased understanding of water movement in these systems could potentially allow for more efficient utilization of geothermal energy as a renewable resource.

Self-powered desalination of geothermal saline groundwater:technical feasibility

This theoretical study shows the technical feasibility of self-powered geothermal desalination of groundwater sources at <100 °C. A general method and framework are developed and then applied to specific case studies. First, the analysis considers an ideal limit to performance based on exergy analysis using generalised idealised assumptions. This thermodynamic limit applies to any type of process technology. Then, the analysis focuses specifically on the Organic Rankine Cycle (ORC) driving Reverse Osmosis (RO), as these are among the most mature and efficient applicable technologies. Important dimensionless parameters are calculated for the ideal case of the self-powered arrangement and semi-ideal case where only essential losses dependent on the RO system configuration are considered. These parameters are used to compare the performance of desalination systems using ORC-RO under ideal, semi-ideal and real assumptions for four case studies relating to geothermal sources located in India, Saudi Arabia, Tunisia and Turkey. The overall system recovery ratio (the key performance measure for the self-powered process) depends strongly on the geothermal source temperature. It can be as high as 91.5% for a hot spring emerging at 96 °C with a salinity of 1830 mg/kg.

Exergoeconomic analysis of solar organic rankine cycle for geothermal air conditioned net zero energy buildings

This study is an attempt at achieving Net Zero Energy Building (NZEB) using a solar Organic Rankine Cycle (ORC) based on exergetic and economic measures. The working fluid, working conditions of the cycle, cycle configuration, and solar collector type are considered the optimization parameters for the solar ORC system. In the first section, a procedure is developed to compare ORC working fluids based on their molecular components, temperature-entropy diagram and fluid effects on the thermal efficiency, net power generated, vapor expansion ratio, and exergy efficiency of the Rankine cycle. Fluids with the best cycle performance are recognized in two different temperature levels within two different categories of fluids: refrigerants and non-refrigerants. Important factors that could lead to irreversibility reduction of the solar ORC are also investigated in this study. In the next section, the system requirements needed to maintain the electricity demand of a geothermal air-conditioned commercial building located in Pensacola of Florida is considered as the criteria to select the optimal components and optimal working condition of the system. The solar collector loop, building, and geothermal air conditioning system are modeled using TRNSYS. Available electricity bills of the building and the 3-week monitoring data on the performance of the geothermal system are employed to calibrate the simulation. The simulation is repeated for Miami and Houston in order to evaluate the effect of the different solar radiations on the system requirements. The final section discusses the exergoeconomic analysis of the ORC system with the optimum performance. Exergoeconomics rests on the philosophy that exergy is the only rational basis for assigning monetary costs to a system’s interactions with its surroundings and to the sources of thermodynamic inefficiencies within it. Exergoeconomic analysis of the optimal ORC system shows that the ratio Rex of the annual exergy loss to the capital cost can be considered a key parameter in optimizing a solar ORC system from the thermodynamic and economic point of view. It also shows that there is a systematic correlation between the exergy loss and capital cost for the investigated solar ORC system.