Printer Friendly

Groundwater Chemistry and Overpressure Evidences in Cerro Prieto Geothermal Field.

Dedicated in memory of Dr. Ramiro Rodriguez C.

1. Introduction

Different authors have observed that the origin of some brines around the world can be caused by diagenesis evolution and porewater trapped during burial [1-3]. Commonly, porewater shows characteristics according to the depositional environment [4]. The different stages of diagenesis are composed of burial and compaction of material in sedimentary basins; the sandstones and shales undergo changes in their physical, chemical, textural, and mineralogical properties which are reflected in the sediment density, the compaction of the granular package, and loss of porosity [5]; additionally, a chemical alteration occurs in sandstones generating cementation and lithification as a product of chemical precipitation affecting detrital grains, dissolution, recrystallization, or mineralogical alteration [6]. Silica is the most abundant cementing agent in sandstones, compared to calcium carbonate, either calcite or aragonite, since the latter dissolves more easily in contact with groundwater, mainly the second [7]. In the diagenetic stages, interstitial fluids are constantly lost in most shales [8]; this phenomenon can occur at low depth and can last for thousands or hundreds of thousands of years or have a longer duration and greater depth; but when it suddenly occurs and the permeability of the medium is so low that it does not allow the interstitial fluid to leak, an increasing degree of stress can be generated [9]; this geological event is known as overpressure [10]. Overpressure is an abnormally high pressure in the subsoil that exceeds the hydrostatic pressure at a certain depth; this pressure is carried out in the pores where the pressure of the interstitial fluids increases as the overcoat increases [11, 12]. Overpressure indicates that the high pressures developed during compaction do not dissipate efficiently [10] and can generate hydraulic fracturing in the system (stress applied to compressible rock and fluid expansion). These processes generate overpressured fluids related to porosity reduction, changes on porewater flow, and diagenetic reactions due to compaction and disequilibrium in a sedimentary basin [12-14]. Also overpressure can be related to chemical compaction due to changes in mineralogy (e.g., ion exchange, dissolution/precipitation) or to diagenetic processes and fluid expansion by thermal water in pore space.

Cerro Prieto Geothermal Field (CP), located in northwestern Mexico (32[degrees]24'43"N, 115[degrees]14'41"W), is a brine with high-temperature geothermal system characteristics. Several studies about the origin and behavior of CP groundwater have been reported [15-18]. According to geological evidences, a large accumulation of sedimentary material from a continental and marine origin, overlying the depositional basin, is related to the origin of brine [17]. The sedimentary material shows diagenetic evolution evidences and porewater trapped between sediment grains during burial processes. The porewater is saline with high [Cl.sup.-], [Na.sup.+], [Ca.sup.+], and [K.sup.+] concentration; in geothermal brines this characteristic is common [1, 2, 10]; the composition depends mainly on the primary origin, mineralogical composition of the sediments and their modification due to diagenetic processes (e.g., facies distribution), and hydrothermal characteristics [15, 19]. Among the most evident diagenetic processes in CP are cementation, mineral replacement, recrystallization, authigenesis, and growth of concretions and nodules [16]. On the other hand, in deep sedimentary basins as CP, mechanical processes of deformation related to burial mechanisms are common; likewise, hydrostatic and lithostatic conditions increase with depth due to an increase of the superposed fluids hydraulically connected, through the pore and the pressure exerted by sediments overload [10]. If pore pressure in deep aquifers like CP is higher than expected from hydrostatic conditions, anomalous pressure (overpressure) can be generated; overpressure is common mainly within 2-4.5 km depth [10-12].

According to geochemical evidences, the origin of geothermal brine at CP could be governed by mixing processes related to a hydrothermal environment and the sedimentary material located at depth which shows burial diagenesis evolution with hydrogeochemical evidences of an overpressurized environment. The aim of the present study was to evaluate hydrogeochemical behavior of geothermal groundwater and its relation with diagenetic processes, including overpressure.

2. Localization

Cerro Prieto Geothermal Field (CP), located in Mexicali Valley, SE of Mexicali City, in Baja California State, Mexico (Figure 1), is a Basin of Salton Sea [17]. Climate is arid with temperatures up to 40[degrees]C in July and to 4[degrees]C in winter. The average annual precipitation is 55 mm/year and the average annual evaporation is 2200 mm/year [23]. Groundwater at CP is extracted from geothermal wells that are in constant exploitation to generate electricity. The Comision Federal de Electricidad of Mexico (CFE) operates and manages the Geothermal Field. CP power production is up to 720 MW and is composed of five individual units: CP1, CP2, CP3, CP4, and CP5; each unit has a total capacity of production with a specific number of production wells. Sampled wells are indicated in Figure 1. Their localization in the five individual CP units is included in Table 1. All the wells are located in zone "beta," 1500 to 3100 m depth [17].

2.1. CP Geology. The lithology around CP is composed of gneiss (quartz-feldspars), shale (quartz-mica), marble, amphibolite, and quartzite from Permic to Jurassic [24] and metamorphic, granitic, and granodiorite rocks which are intruded by batholitic rocks [16], together with dacite and andesite from Miocene and rhyodacite from Quaternary [25]. The tectonic basin was filled by sedimentary material that, due to burial, compaction, and diagenesis processes, evolved to gray shales from Late Miocene (shales and silt shales that vary from light gray to black); this unit overlies the granitic basement and the mafic intrusive and is interlayered by permeable sandstones (composed of quartz and feldspars; arkoses type). The thickness is near 3000 m [21, 26]. Immediately above a layer of brown shale (shale and silt shale), the gray shale covers interlayered permeable sandstones and sands cemented by carbonates, about 500 m thickness [21]. In these zones, a rapid distribution of geothermal fluids enhances the recharge. Erratic mudstone layer and unconsolidated clastic sediments (clays, silts, sand, and little gravel) overlie all the previous units. The thickness of these units is between 400 and 2500 m [26]. Sedimentary material shows evidence of the diagenetic evolution and recrystallization processes due to an incipient low grade metamorphism.

The geological evolution of CP is a complex blend of rifting, rapid deltaic sedimentation, and large scale strike-slip faulting located within the Salton Basin [16, 19, 21, 27]. The Geothermal Field is placed in a shear zone where NW-SE and NE-SW fault systems intersect. The more important faults are Cucapa, Imperial, Cerro Prieto, and Michoacan [21]. This fault system is part of a major regional lineament that penetrates deep into the crustal and basement rocks and serves as conduit for geothermal flow. The system originates in a tectonic basin of 5200 m depth, filled by alluvial and deltaic sediments from Tertiary to Quaternary [21].

Vonder Haar and Howard [27] observed that, in sandstone and shale units, a mineral dissolution/precipitation took place along microfractures, originating secondary porosity and newly precipitated hydrothermal minerals, causing a reduction of permeability. Likewise, Elders et al. [20] observed cementation, mineral replacement, recrystallization, authigenesis, and growth of concretions and nodules; these processes are related to diagenesis.

2.2. CP Hydrogeology. Some authors consider that CP brine may have been formed from marine evaporite dissolved and partly by evaporated Colorado River water [16, 19, 21, 28-31]. However, according to geological evidences, a large accumulation of sedimentary material, overlying the depositional basin, from a continental and marine origin and mixing with meteoric water [30-33], was related to the origin of brine [17]. The sedimentary material shows diagenetic evolution evidences and porewater trapped between grains during burial processes. Isotopic evaluations ([sup.18]O, [sup.2]H) and chemical analysis ([Cl.sup.-] and Br) elaborated by Coplen [34] and Birkle et al. [35] suggest that Salton Sea was the probable predecessor of high chlorinated groundwater of CP.

3. Methodology

A hydrogeological and hydrogeochemical study was carried out in geothermal groundwater samples from CP following standard methods of APHA-AWWA [38]. Water samples were collected from geothermal wells that are in constant exploitation. Temperature, pH, and conductivity were measured in the field during the summer of 2010 and calibrated to the water temperature at each site. The chemical analyses for the major elements, B and Si[O.sub.2], were performed at the Analytical Chemistry Laboratory of the Geophysics Institute, UNAM, Mexico (the laboratory participates in international calibration exercises of chemical analyzes of geothermal waters). Boron was determined by colorimetry through reactions with carminic acid (Method 4500-B C) APHAAWWA [38]. Si[O.sub.2] was determined by atomic absorption spectrophotometry with flame and UV-visible spectroscopy (molybdosilicic acid method). Major ions were analyzed following standard methods [38]. HC[O.sub.3.sup.-] and C[O.sub.3.sup.2-] were determined by volumetry (titrating with HCl), [Ca.sup.+] and [Mg.sup.2+] were determined by volumetry (titrating with EDTA), [Cl.sup.-] was determined by potentiometry with selective electrodes (4500-[Cl.sup.-]) [38], [Na.sup.+] and [K.sup.+] were determined by atomic emission spectrophotometry (3500-[Na.sup.+] and [K.sup.+]), and S[O.sub.4.sup.2-] was determined by turbidimetry (4500-S[O.sub.4.sup.2-]). Analytical quality was assessed through ionic balance (less than 10%) and the use of certified (NIST) reference solutions.

In order to evaluate brine evolution of CP, different techniques were used: (a) Carpenter [39] evaluated the behavior of major elements using a plot with concentrations as a function of dissolved chloride concentration, considering the composition of seawater during evaporation and diagenesis; using chemical results of CP groundwater a similar evaluation was elaborated; (b) Davisson and Criss [40] devised a diagram to determine the geochemical evolution of mineralogy in brines applying an evaluation of [Na.sub.(deficit)] and [Ca.sub.(excess)] in water samples; hydrogeochemical results of samples from CP were evaluated with this diagram; (c) Boschetti [3] considers B-Cl concentrations to explain the evolutionary process in the groundwater; B-Cl diagram was used to determine the dominant geological environment in CP. A geochemical simulation with the measured water concentrations was carried out using the Phreeqc[c] program to determine saturation indexes.

In order to determine pressure conditions, a calculation was performed using hydrostatic and lithostatic properties from CP, considering geological characteristics of the study area and theoretical information about some basin environments. To estimate pressure conditions ([P.sub.l](h)) exerted by a geological column at a depth (h) offshore in a geologic formation, (1) was used [12]:

[mathematical expression not reproducible] (1)

It is necessary to take into account pure water column weight at sea level h = 0, depth of seawater column [h.sub.w], seawater density [[rho].sub.sea] = 1100 kg/[m.sup.3], gravity constant g = 9.78 (m/[s.sup.2]), and material density (rock or sediments) [[rho].sub.b]. The integral about overburden weight of sediments can be replaced by a sum of the individual weights of layers [12]:

[mathematical expression not reproducible] (2)

where it is important to consider the number of layers (ni) with thickness [d.sub.i] (i is the layer number measured in m), rocks density [[rho].sub.ri] (km/[m.sup.3]), porosity [[phi].sub.i], and water density [[rho].sub.w] (which can change with salinity variation, while temperature and pressure dependence is relatively small or negligible).

To estimate pressure in a reservoir unit from CP it was considered that the study area is located onshore, few meters above seawater level; therefore the integral including the weight of seawater column is zero; the distinctive strata units, their thickness, and depth of each sampled well were considered; the porosity values reported by Olson [41] and Hiriart Le Bert [42] in the geologic material from CP were used (between 0.15 and 0.25), for lithological units without reported porosity; theoretical values proposed by distinct authors studied unlike geological material and depth were used (e.g., mudstone and shales, slate, quartz, feldspars, uncemented sandstones, or sandstones as reservoirs) [43-46].

4. Results and Discussion

4.1. Major Elements in CP Groundwater. Table 1 shows the major solutes concentration in CP groundwater.

In Table 1 chemical results from the studied wells are shown. Total dissolved solids were calculated from conductivity measures showing values between 12945.5 and 49251.2 mg/L. These values are comprised between saline water (>1,000 mg/L) and brine (>35,000). All the sampled wells are located at beta reservoir [17, 47]. To explain geochemical variations in CP brine it is necessary to evaluate the composition and removal of solutes by salt precipitation according to diagrams proposed by Carpenter [39] where the circle represents the solute-chloride composition of normal seawater; the line represents the limit between evaporation-dilution curves of seawater and freshwater (Figure 2).

Figure 2 shows that groundwater in CP is dominated by high concentrations of [Na.sup.+], [K.sup.+], [Ca.sup.+], and [Cl.sup.-]. When [Cl.sup.-] concentration increases [Na.sup.+], [K.sup.+], and [Ca.sup.+] concentrations also increase in ratios 1:1, 1:1, and 2: 1, respectively (Figures 2(a), 2(b), and 2(c)). Bicarbonate ion, [Mg.sup.2+], and S[O.sub.4.sup.2-] have a low concentration (Figures 2(d), 2(e), and 2(f)), [Na.sup.+] values are parallel and near coincident with their respective evaporation-dilution curve (Figure 2(a)), [Ca.sup.+] values cross the evaporation-dilution curve (Figure 2(b)), and [K.sup.+] also increases with [Cl.sup.-] but the values are enlarged relative to the seawater evaporation-dilution curve (Figure 2(c)). In CP brine, [Mg.sup.2+] concentrations lie well below the seawater evaporation trajectory indicating significant depletion of the element. According to Hanor [48] and Kharaka and Hanor [1] [Mg.sup.2+] concentrations in brines decrease when temperature increases in the subsurface, and when alkalinity decreases [49]. Evaporation of continental waters has a more variable concentration range in water samples from CP due to evaporating mixtures of continental and marine water with meteoric water [50].

In CP brine depleted and enriched concentrations of some major elements are a consequence of reactions linked with the hydrothermal processes and water-rock interactions. Major elements concentrations are controlled by the alteration and formation of minerals like feldspar-K, plagioclases, quartz, biotite, amphibole, chlorite, pyrite, wairakite, prehnite, muscovite, epidote, and talc, as reported in CP by Elders et al. [20] and Izquierdo et al. [16] and shown by calculated saturation index values (Figures 3 and 4).

The [K.sup.+] origin is restrained by alteration of feldspar-K, illite, and biotite and by muscovite formation. Very low S[O.sub.4.sup.2-] and HC[O.sub.3.sup.-] concentrations in CP could be inhibited by the water interactions with anhydrite, dolomite, talc, pyrite, or calcite (Table 1; Figures 2 and 3); Elders et al. [20] report pyrite formation at depth. Low concentration of [Mg.sup.2+] and high concentration of [Ca.sup.+] could be related to dolomitization of limestone as major source of [Ca.sup.+]; and low [Mg.sup.2+] contents are associated with the evolution of chlorites and micas when temperature and depth increase according to mineralogy reported in CP (Figure 3); similar behavior has been reported previously in hydrothermal brines with diagenetic evolution evidences [51] (Figure 3). Albite reactions at depth at high temperature can be linked with the slight [Na.sup.+] decrease [16, 30,31].

4.2. Saturation Index. Results of saturation index calculations are shown in Figure 4. From these results amorphous silica (Si[O.sub.2] am), albite, k-feldspar, and in some cases k-mica show a behavior close to equilibrium with the fluid. Besides, quartz, chalcedony, talc, and crysocole are oversaturated. Dolomite, calcite, and aragonite are undersaturated at some sites and oversaturated at others, in agreement with Figures 2 and 3.

4.3. [Na.sub.(deficit)]-[Ca.sub.(excess)] Plot. Figure 5 was used to explain the initial composition of brines and the nature of fluid-rock interactions. In the diagram, Basinal Fluid Line (BFL) is a straight line with a unit slope that indicates a 2Na-1Ca exchange relationship [40]; BFL represents the effect of plagioclase albitization on water composition. Seawater evaporation trajectory is a representation of the natural trends for seawater evaporation which is formed by large positive [Na.sub.(deficits)] and small negative [Ca.sub.(excess)]; reactions involving seawater evaporation follow a vertical descent and afterwards produce large deficits along a horizontal trend. Halite dissolution in seawater or freshwater can produce negative values along a slope of 1: 4.

To determine the origin and geochemical evolution of CP brine, an evaluation of Na(deficit) and [Ca.sub.(excess)] was applied to explain the initial composition and nature of fluid-rock interaction (Figure 5). All the analyzed water samples from CP were located right and over the seawater evaporation trajectory (SET) in the [Na.sub.(deficits)]-[Ca.sub.(excess)] diagram (Figure 5) indicating that the fluids are a product of brine that passed the point of halite precipitation evaporated (along of [Na.sub.(deficits)] axis). The horizontal line of CP has a large positive [Na.sub.(deficits)] and a small negative [Ca.sub.(excess)], but the fluid is more enriched in [Ca.sub.(excess)] than expected from seawater evaporation.

Dolomitization produces elevated Ca contents, increasing [Ca.sub.(excess)] without changing the [Na.sub.(deficit)], and can explain the phenomenon observed in CP. Other dissolved minerals could interact in the evolutionary processes to brine aside from the above-mentioned ones, such as calcite, anhydrite, quartz, halite, and illitization to smectite transformation. Often, these geochemical processes are treated separately but in some situations the origin is linked to a mixing of processes and is not mutually exclusive [1]. Illite formation involves reactions relevant for diagenesis: (a) release of water during transformation of feldspar to kaolinite and smectite; (b) potassium and silica budget [2, 52]. Some of these processes can occur in CP.

Geology of CP shows equilibrium with clays (Figure 4); groundwater shows high concentrations of [K.sup.+]. Also reactions between calcite, illite, and [K.sup.+] to form K-feldspar could be involved in brine evolution. According to the results (Figures 2 and 5) the origin of CP groundwater is a consequence of an evolution of dissolved evaporative products (e.g., halite), residual water remaining during the dissolution precipitation of seawater evaporites, and different water-rock interactions (e.g., clays, siltstones, and shales) but could have a slight contribution to exchange reactions 1Ca for 2Na in the aquifer by albitization of plagioclases, which could change the ionic composition of fluids (Figure 5), according to the following reaction proposed by Carpenter [39], Hanor [48], and Demir and Seyler [53]:

Ca[Al.sub.2][Si.sub.2][O.sub.8(an)] + 4Si[O.sub.2(aq)] + 2[Na.sup.+] = 2NaAl[Si.sub.3][O.sub.8](ab) + [Ca.sup.2+].

4.4. Overpressure

4.4.1. B-Cl Plot. Some hydrogeochemical evidences (Figures 2 and 5) shown in a B-Cl plot (Figure 6) indicate that overpressurized fluids participate in the evolutionary process of the CP geothermal brine. Boron behavior helps to define these geologic environments because B is adsorbed by clay minerals and is released into the fluid in a deep environment, mainly when tectonic stresses by vertical and/or lateral compaction are high and temperature increases with depth, and stronger geologic deformations are generated [10, 12,54].

4.4.2. Pressure Conditions Estimated due to Overburden. Some authors [55, 56] estimated pressure at depth at CP considering the host rock density at a specific depth. The obtained values were in the range between 0.5 and 42 MPa. Those calculations considered only geothermal water conditions from the extraction zone and no lithologic information was used from the stratigraphic columns of each site. To confirm hydrogeochemical evidence of overpressure, a calculation to determine pressure conditions was performed using hydrostatic and lithostatic properties from CP and theoretical information about some basin environments (Figure 7). It is necessary to consider that normal pressure increases with depth according to the hydrostatic pressure gradient (10 MPa/Km), higher or lower values of this gradient, and their associated depths are considered abnormal pressures (e.g., overpressure). Lithostatic pressure is equivalent to the total charge of the overlaying sediments in a geological formation and increases according to the lithostatic pressure gradient (23 MPa/Km) [10-12, 57].

In general, overpressured systems can take place when porewater is not expelled from rock at a proper rate, remaining under hydrostatic pressure. Overpressured volume rock must be trapped by low permeability layers where fluids move slowly, even when there is high pressure in the environment. The overpressure affects the effective stress that acts between the grains within the rock and generates a change in the compaction. In many areas of active sedimentation rate around the world, porewater pressure in deep groundwater (>1km) is higher than would be expected from hydrostatic circumstances [10, 57-59]. It is necessary to consider thermal expansion in the pore space (increasing volume) and increment of the system temperature by thermal conditions and by fluid movement and mineral phases transition [11]. According to Swarbricket al. [11] and Kauerauf and Hantschel [12] secondary overpressure by chemical cementation may occur at large depths when porosity is reduced by dissolution-diffusive transport-precipitation of silica cement (temperature affects diffusion-precipitation rate) or by fluid expansion processes when gas or thermal solutions are originated in highly permeable facies interconnected locally at certain depth levels generating compaction, rearrangement of grains, and reduction in the pore space. Some of these conditions could occur in CP aquifer [22].

The results obtained confirm that in CP overpressure is present (Figures 7 and 8). Positive anomalies with an increase of depth are shown between 0.6 and 3.1 km; this phenomenon can be caused by a hydraulic seal. Overpressure coincides with lithology variation in CP [22], when sedimentary material composed of sand and clay changes to shale and sandstone due to clay and sand diagenesis and shale compaction (Figure 7). In CP, sandstone is hosted by a low permeability layer of shale and siltstones creating adequate conditions for this process. Overpressure in CP brine could occur as a consequence of a rapid sedimentation, deposition, and accumulation rates of fine-grained material along time and due to the stress increase in sediments, compressibility, or expansion of fluids by hydrothermalism.

Geological evidence in CP indicates that mineralogical changes occur at depth, mainly by diagenetic processes (Figure 5). According to the geochemical behavior observed in CP, it is possible that mineral dissolution precipitation processes generated this phenomenon; osmosis, buoyancy, and tectonic or magmatic process can generate changes in mineralogy (e.g., feldspar to illite or smectite to illite conversion), diagenesis, and carbonate or silicate cementation. In CP, sandstones and shales units reveal clogging, mineral dissolution, and mineral precipitation along microfractures as indicated by Vonder Haar and Howard [27]. In the study area the calcite dissolution and/or cementation more likely develops where cold water interacts with hotter rocks or precipitation of quartz and k-feldspar occurs when hot waters interact with colder rocks [16, 20]. Exchange of [Na.sup.+] by [K.sup.+] with a higher ion radius is carried out in mineral transformations and increases the volume of the solid matrix [12]; this process is controlled by temperature and [K.sup.+] availability in minerals.

The permeability of sandstones in CP facilitates reactions between rocks and hydrothermal fluids (e.g., dissolution of some minerals); these reactions can reduce or increase the porosity and also generate secondary fracturing or microfracturing and modify physical properties [16, 21], which could cause overpressure (Figures 7 and 8). Sandy shale and siltstone facies in CP are most amenable to increased microfracturing. In sandstones (where high temperature dominates) and shales from CP, mineral dissolution precipitation takes place along microfractures, which originates secondary porosity and newly precipitated hydrothermal minerals, causing a reduction of permeability.

5. Geothermometry

Alkali feldspar geothermometers are the most used tool to determine chemical equilibrium in fluids at depth in a geothermal system [60]. Na/K and Na-K-Ca geothermometers were developed to evaluate the temperature in high enthalpy geothermal systems [61]; these geothermometers are less affected by chemical reequilibration in mixing zones, but the calculated temperature may be affected by mixing with cold water or by deposition of aluminum-rich minerals or clays [62]. Na/K geothermometers are adapted for temperatures between 180 and 350[degrees]C. Na/K geothermometers have been generally used to estimate temperatures at CP [36, 37, 62]. We applied (4) and (5) [61, 62] to evaluate the temperature of CP reservoir:

T = (1178/(l.47 + log (Na/K))) - 273.15 (4)

T = (1289/(1.615 + log (Na/K))) - 273.15. (5)

The temperature range of CP reservoir, calculated with Na/K geothermometer, varies between 236 and 306[degrees]C for (4) and between 251 and 318[degrees]C for (5) (Table 2).

6. Conclusions

CP shows intermediate brine characteristics (Na-Ca-Cl) with high [K.sup.+], [Ca.sup.+], and salinity contents; the relation [Na.sup.+]/[Cl.sup.-] is less than 1.

The porewater composition in CP evolved from its primary origin and was modified by the interaction with minerals of the sedimentary material. Brine characteristics were acquired by deep-burial diagenesis processes and low grade metamorphism at high temperatures. Results show geochemical evidence of overpressured fluids due to compaction.

Groundwater samples from CP show a mixing of marine and continental water; this situation is partially related to a continental and evaporative precursor. The hydrogeochemical evidence indicates that the sedimentary material has porewater between grains, which was trapped during burial processes. The diagenetic processes could have generated high concentrations of [Cl.sup.-], [Na.sup.+], [K.sup.+], and Ca. Calcium enrichment, [Na.sup.+] depletion, and [K.sup.+] release could have a relation with a contribution of exchange reactions 1Ca for 2Na in the aquifer by albitization of plagioclases or by illitization processes, respectively, and precipitation of secondary minerals. High [K.sup.+] and low [Mg.sup.2+] contents are related to alteration of feldspar-K, illite, biotite, and muscovite formation. Ca-Na exchange with plagioclases could be a geochemical control on the fluids of CP and may directly explain slight [Ca.sub.(excess)] and [Na.sub.(deficit)] in the brine.

Overpressure in CP is related to burial mechanisms; secondary overpressure is related to chemical pressure by mineralogical changes and by fluid expansion, which increases with depth. The magnitude of overpressure may be produced by some characteristics of deposit formation (burial), permeability evolution of sedimentary material, and the compressibility of rock and fluid. Secondary overpressure in the system related to chemical pressure and porosity changes due to mineral dissolution can be generated at large depths. Fluid expansion takes place in the reservoir which generates compaction, rearrangement of grains, and reduction of pore space. Boron release at overpressure conditions can be related to high contents of [K.sup.+] in water.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


The authors thank Aguayo A., Ceniceros N., and Cruz O. for chemical determinations. The authors acknowledge Comision Federal de Electricidad for support on sampling within the Cerro Prieto Geothermal Field.


[1] Y. K. Kharaka and J. S. Hanor, "Deep Fluids in the Continents: I. Sedimentary Basins," Treatise on Geochemistry, vol. 5-9, pp. 1-48, 2003.

[2] T. Boschetti, L. Toscani, O. Shouakar-Stash et al., "Salt Waters of the Northern Apennine Foredeep Basin (Italy): Origin and Evolution," Aquatic Geochemistry, vol. 17, no. 1, pp. 71-108, 2011.

[3] T. Boschetti, "Application of brine differentiation and Langelier-Ludwig plots to fresh-to-brine waters from sedimentary basins: Diagnostic potentials and limits," Journal of Geochemical Exploration, vol. 108, no. 2, pp. 126-130, 2011.

[4] A. Arche, "Sedimentologia del proceso fisico a la cuenca sedimentaria," in Consejo Superior de Investigaciones Cientficas, 978-84-00-09145-3, pp. 1-1273, Madrid, Spain.

[5] S. J. Blott and K. Pye, "Particle shape: A review and new methods of characterization and classification," Sedimentology, vol. 55, no. 1, pp. 31-63, 2008.

[6] A. Ceriani, A. Di Giulio, R. H. Goldstein, and C. Rossi, "Diagenesis associated with cooling during burial: An example from Lower Cretaceous reservoir sandstones (Sirt basin, Libya)," AAPG Bulletin, vol. 86, no. 9, pp. 1573-1591, 2002.

[7] F. W. Witkowski, D. J. Blundell, P Gutteridge, A. D. Horbury, N. H. Oxtoby, and H. Qing, "Video cathodoluminescence microscopy of diagenetic cements and its applications," Marine and Petroleum Geology, vol. 17, no. 10, pp. 1085-1093, 2000.

[8] C. H. Moore, "Carbonate reservoirs. Porosity evolution and diagenesis in a sequence stratigraphic framework," Developments in Sedimentology, vol. 55, 2001.

[9] M. S. Fantle, K. M. Maher, and D. J. Depaolo, "Isotopic approaches for quantifying the rates of marine burial diagenesis," Reviews of Geophysics, vol. 48, no. 3, Article ID RG3002, 2010.

[10] K. M. Hiscock and V F. Bense, Hydrogeology Principles And Practice, Wiley Blackwell, 2nd edition, 2014.

[11] R. E. Swarbrick, M. J. Osborne, and G. S. Yardley, "Comparison of overpressure magnitude resulting from the main generating mechanisms," in Pressure Regimes in Sedimentary Basins and Their Prediction, A. R. Huffman and G. L. Bowers, Eds., vol. 76, pp. 1-12, American Association of Petroleum Geologists Memoir, 2002.

[12] A. I. Kauerauf and T. Hantschel, Fundamentals of Basin and Petroleum Systems Modeling, Springer Science & Business Media, 2009.

[13] O. Walderhaug, P. A. Bjorkum, P. H. Nadeau, and O. Langnes, "Quantitative modelling of basin subsidence caused by temperature-driven silica dissolution and reprecipitation," Petroleum Geoscience, vol. 7, no. 2, pp. 107-113, 2001.

[14] A. Zanella, P R. Cobbold, and C. Le Carlier de Veslud, "Physical modelling of chemical compaction, overpressure development, hydraulic fracturing and thrust detachments in organic-rich source rock," Marine and Petroleum Geology, vol. 55, pp. 262-274, 2014.

[15] E. Portugal and M. P Verma, "Hidroqulmica de la laguna de evaporacion en Cerro Prieto," in Ingenieria hidraulica en Mexico, vol. 16, pp. 153-174, Baja California, Mexico, 2 edition, 2001.

[16] G. Izquierdo, A. Aragon, E. Portugal, V. M. Arellano, J. de Leon, and J. Alvarez, "Mineralogia de la zona mineralizada de silice-epidota (ZMSE) del yacimiento geotermico de Cerro Prieto B.C Mexico," Geotermia, vol. 19, no. 2, pp. 2-12, 2006.

[17] V. M. Arellano, R. M. Barragan, A. Aragon, M. H. Rodriguez, and A. Perez, "The Cerro Prieto IV (Mexico) geothermal reservoir: Pre-exploitation thermodynamic conditions and main processes related to exploitation (2000-2005)," Geothermics, vol. 40, no. 3, pp. 190-198, 2011.

[18] M. A. Armienta, R. Rodriguez, N. Ceniceros et al., "Groundwater quality and geothermal energy. The case of Cerro Prieto Geothermal Field, Mexico," Journal of Renewable Energy, vol. 63, pp. 236-254, 2014.

[19] J. M. Camacho-Hernandez, "Zonas de alteracion hidrotermal y condiciones actuales del yacimiento: un enfoque para determinar zonas productoras al oriente del campo geotermico de Cerro Prieto, BC," Geotermia. Revista Mexicana de Geoenergia, vol. 22, no. 2, pp. 35-45, 2009.

[20] W. A. Elders, J. R. Hoagland, S. D. McDowell, and J. M. Cobo, "Hydrothermal mineral zones in the geothermal reservoir of Cerro Prieto," Geothermics, vol. 8, no. 3-4, pp. 201-209, 1979.

[21] H. H. Lira, "Actualization del modelo geologico conceptual del yacimiento geotermico de Cerro Prieto, BC," Geotermia, vol. 18, no. 1, pp. 37-46, 2005.

[22] A. L. Pena, C. I. Puente, and C. E. Diaz, "Modelo geologico del campo geotermico de Cerro Prieto. Geothermal-Energy," Comision Federal de Electricidad, pp. 29-52, 1979, .pdf.

[23] E. Portugal, J. Alvarez, and B. I. Romero, "Hydrochemical and isotopical tracers in the lacustrine aquifer of the Cerro Prieto area, Baja California, Mexico," Journal of Geochemical Exploration, vol. 88, no. 1-3, pp. 139-143, 2006.

[24] M. Siem, The structure and petrology of Sierra El Mayor [Master, thesis], University of San Diego State, Northeastern Baja Calfornia, Mexico, 1992.

[25] A. L. Quintanilla-Montoya and F. Suaorez-Vidal, "Cerro Prieto and its relation to the Gulf of California spreading centers," Ciencias Marinas, vol. 22, no. 1, pp. 91-110, 1996.

[26] R. J. M. Cobo, Configuration de los cuerpos litologicos de lodolita, lutita cafe, lutita gris, zonas de silice y epidota, y sus relaciones con la tectonica del campo geotermico de Cerro Prieto. Proceedings of the Third Symposium on the Cerro Prieto Geothermal Field, Mexico, 1981.

[27] S. Vonder Haar and J. H. Howard, "Intersecting faults and sandstone stratigraphy at the cerro prieto geothermal field," Geothermics, vol. 10, no. 3-4, pp. 145-167, 1981.

[28] A. Manon, E. Mazor, M. Jimenez, A. Sanchez, J. Fausto, and C. Zenizo, "Extensive geochemical studies in the geothermal field of Cerro Prieto, Mexico," Tech. Rep. LBL-7019, 1977

[29] E. Mazor and A. Manon M., "Geochemical tracing in producing geothermal fields: A case study at Cerro Prieto," Geothermics, vol. 8, no. 3-4, pp. 231-240, 1979.

[30] A. H. Truesdell, R. O. Rye, F. J. Pearson Jr. et al., "Preliminary isotopic studies of fluids from the Cerro Prieto geothermal field," Geothermics, vol. 8, no. 3-4, pp. 223-229, 1979.

[31] A. H. Truesdell, J. M. Thompson, T. B. Coplen, N. L. Nehring, and C. J. Janik, "The origin of the Cerro Prieto geothermal brine," Geothermics, vol. 10, no. 3-4, pp. 225-238, 1981.

[32] H. A. Truesdell, M. J. Lippmann, and H. Gutierrez-Puente, "Evolution of the Cerro Prieto Reservoirs under exploitation," in Proceedings of the Anual Meeting of the Geothermal Resources Council, pp. 1-7, Burlingame, Calif, USA, 1997

[33] M. J. Lippmann, A. H. Truesdell, and K. Pruess, "The control of fault H on the hydrology of the Cerro Prieto III Area," in Proceedings of the Twenty-fifth Workshop on Geothermal Reservoir Engineering Standford University, Standford, Calif, USA, 2000.

[34] T. B. Coplen, "Origin of geothermal waters in the Imperial Valley of southern California. Cooperative Investigation of Geothermal Resources in the Imperial Valley and their Potential Value for Desaltine of Water and other purposes," R. W. Rex, Ed., Rwerslde Report IGPP-UCR-72-33, pp. E1-E31, University of California, 1972.

[35] P Birkle, E. P Marin, D. L. Pinti, and M. C. Castro, "Origin and evolution of geothermal fluids from Las Tres Virgenes and Cerro Prieto fields, Mexico - Co-genetic volcanic activity and paleoclimatic constraints," Applied Geochemistry, vol. 65, pp. 36-53, 2016.

[36] D. Nieva and R. Nieva, "Developments in geothermal energy in Mexico-part twelve. A cationic geothermometer for prospecting of geothermal resources," Heat Recovery Systems and CHP, vol. 7, no. 3, pp. 243-258, 1987.

[37] S. P Verma and E. Santoyo, "New improved equations for Na/K, Na/Li and SiO2 geothermometers by outlier detection and rejection," Journal of Volcanology and Geothermal Research, vol. 79, no. 1-2, pp. 9-23, 1997.

[38] AWWA, APHA, and WWF, Standard methods for the Examination of Water and Wastewater, American Public health Association, The American Water Works Association, Association Water Environment Federation, Washington, D.C., USA, 2005.

[39] A. B. Carpenter, "Origin and chemical evolution of brines in sedimentary basins Oklahoma," Geological Survey Circular, vol. 79, pp. 78-88, 1978.

[40] M. L. Davisson and R. E. Criss, "Na-Ca-Cl relations in basinal fluids," Geochimica et Cosmochimica Acta, vol. 60, no. 15, pp. 2743-2752, 1996.

[41] E. R. Olson, "Oxygen and Carbon isotopes studies of calcite from the Cerro Prieto Geothermal Field," in Proceedings of the First Symposium on the Cerro Prieto Geothermal Field, Baja California, Mexico, 1978.

[42] Hiriart Le Bert, Evaluacion de la Energia Geotdrmica en Meoxico Informe para el Banco Interamericano de Desarrollo y la Comision Reguladora de Energia, 2011.

[43] W. J. Harrison and L. L. Summa, "Paleohydrology of the Gulf of Mexico basin," American Journal of Science, vol. 291, no. 2, pp. 109-176, 1991.

[44] E. D. Pittman and R. E. Larese, "Compaction of lithic sands: experimental results and applications," The American Association of Petroleum Geologists Bulletin, vol. 75, pp. 1279-1299, 1991.

[45] J. Gluyas and C. A. Cade, "Prediction of porosity compacted sands," in Reservoir Quality Prediction in Sandstones and Carbonates, A. Kupecz, Gluyas J., and S. Bloch, Eds., vol. 69, pp. 19-28, AAPG, Memoir, 1997.

[46] S. N. Ehrenberg and P. H. Nadeau, "Sandstone vs. carbonate petroleum reservoirs: A global perspective on porosity-depth and porosity-permeability relationships," AAPG Bulletin, vol. 89, no. 4, pp. 435-445, 2005.

[47] R. M. Prol-Ledesma, C. Arango-Galvan, and M.-A. Torres-Vera, "Rigorous analysis of available data from cerro prieto and las tres virgenes geothermal fields with calculations for expanded electricity generation," Natural Resources Research, vol. 25, no. 4, pp. 445-458, 2016.

[48] J. S. Hanor, "Origins of saline fluids in sedimentary basins," in Geofluids: Origin, Migration and Evolution of Fluids in Sedimentary Basins, J. Parnell, Ed., pp. 151-174, Geological Society of London, 1994.

[49] A. W. Hounslow, Water Quality Data, Analysis and Interpretation, Taylor & Francis Group, 1995.

[50] J. N. Valette-Silver, J. M. Thompson, and J. W. Ball, "Relationship between water chemistry and sediment mineralogy in the Cerro Prieto Geothermal Field: A preliminary report," Geothermal Energy, pp. 263-273, 1981.

[51] K. H. Wolf and G. V. Chilingarian, "Chapter 1 Introduction," Developments in Sedimentology, vol. 51, no. C, pp. 1-17, 1994.

[52] P. Kaur, N. Chaudhri, A. W. Hofmann et al., "Two-stage, extreme albitization of A-type granites from Rajasthan, NW India," Journal of Petrology, vol. 53, article egs003, no. 5, pp. 919-948, 2012.

[53] I. Demir and B. Seyler, "Chemical composition and geologic history of saline waters in Aux Vases and Cypress Formations, Illinois Basin," Aquatic Geochemistry, vol. 5, no. 3, pp. 281-311, 1999.

[54] G. V Chilingarian, T. F. Donaldson, and T. F. Yen, "Subsidence due to fluid withdrawal," in Developments in Petroleum Science, vol. 519, 41 Elsevier Sciences, 1995.

[55] A. Garcia, F. Ascencio, G. Espinosa, E. Santoyo, H. Gutierrez, and V. Arellano, "Numerical modeling of high temperatura deel Wells in the Cerro Prieto geotermal fiels, Mexico," Geofisica International, vol. 38, pp. 251-260, 1999.

[56] M. J. Lippmann, A. H. Truesdell, and H. Gutierrez-Puente, "What will a 6 km deep well at Cerro Prieto find," in Proceedings of the twenty-first Workshop on Geothermal Reservoir Engineering, Stanford University, 1997, https://pangea.stanford .edu/ERE/pdf/IGAstandard/SGW/1997/Lippmann.pdf.

[57] M. J. Osborne and R. E. Swarbrick, "Mechanisms for generating overpressure in sedimentary basins: A reevaluation," AAPG Bulletin, vol. 81, no. 6, pp. 1023-1041, 1997

[58] A. M. Stueber and L. M. Walter, "Origin and chemical evolution of formation waters from Silurian-Devonian strata in the Illinois basin, USA," Geochimica et Cosmochimica Acta, vol. 55, no. 1, pp. 309-325, 1991.

[59] N. H. Mondol, K. Bjorlykke, J. Jahren, and K. Hoeg, "Experimental mechanical compaction of clay mineral aggregates Changes in physical properties of mudstones during burial," Marine and Petroleum Geology, vol. 24, no. 5, pp. 289-311, 2007

[60] R. O. Fournier and J. J. Rowe, "Estimation of underground temperatures from the silica content of water from hot springs and wet-steam wells," American Journal of Science, vol. 264, no. 9, pp. 685-697, 1966.

[61] R. Sonney, Groundwater flow, heat and mass transport in geothermal systems of a Central Alpine Massif. The cases of Lavey-les Bains, Saint-Gervais-les-Bains and Val d'Illiez. Geochemistry. Universite de Neuchatel, 2010.

[62] F. DAmore and S. Arnorsson, "Isotopic and chemical techniques in geothermal exploration, development and use. Sampling methods, data handling, interpretation," in Geothermometry, S. Arnorsson, Ed., pp. 152-199, International Atomic Energy Agency, Vienna, Austria, 2000.

Ivan Morales-Arredondo, Maria Aurora Armienta, and Nuria Segovia

Instituto de Geofisica, Universidad Nacional Autonoma de Mexico, Cd. Universitaria, 04510 Mexico City, Mexico

Correspondence should be addressed to Ivan Morales-Arredondo;

Received 19 May 2017; Revised 30 October 2017; Accepted 19 November 2017; Published 18 December 2017

Academic Editor: Ian Clark

Caption: Figure 1: Localization of Cerro Prieto Geothermal Field and sampled wells.

Caption: Figure 2: Variations in dissolved (a) [Na.sup.+], (b) [Ca.sup.+], (c) [K.sup.+], (d) [Mg.sup.2+], (e) HC[O.sub.3.sup.-], and (f) S[O.sub.4.sup.2-] concentrations as a function of dissolved chloride concentration. The solid line on each plot is the seawater evaporation-dilution curve for each cation; circle is the ionic composition of seawater.

Caption: Figure 3: Mineralogy and paragenesis reported in lithology from CP, depth of sampled geothermal wells (modified by [20]), and estimated porosity considering geological characteristics.

Caption: Figure 4: Saturation index diagrams.

Caption: Figure 5: Diagram [Na.sub.(deficit)] and [Ca.sub.(excess)]. Diamonds correspond to the CP sampled wells.

Caption: Figure 6: CP overpressure: chloride versus boron concentration (mg/L; fields and paths from [3]). Dashed lines represent binary mixings where specific geological environments dominate. Diamonds correspond to the CP sampled wells.

Caption: Figure 7: CP geology cross section. Stratigraphic column was modified from Lira [21] and Pena et al. [22].

Caption: Figure 8: Plots represent pressure estimated conditions at depth: pressure conditions were calculated with the following values: porewater density [[rho].sub.w] = 1040 kg/[m.sup.3], shale and mudstone density [[rho]] = 2700 kg/m3, sandstone density [[rho]] = 2720 kg/[m.sup.3], slate density [[rho]] = 2750 kg/[m.sup.3], and sedimentary material density [[rho]] = 1650 kg/[m.sup.3] [10, 12].
Table 1: Major elements concentration values at CP
groundwater samples measured in summer 2010.

            (individual     Depth          T                   CE
Well           units)        (Km)    ([degrees]C)     pht    mS/cm

M 104            CPS         1.73        88.4        7.73    20.45
M 117A           CPS         2.89        91.2         7.7      35
M 119A           CP2         3.1         87.9        6.79      35
M 127            CP2         2.5         76.8        7.61     24.8
M 133A           CPS         2.49        80.8        7.49    35.45
M 148A           CP2         2.99        85.8        7.43     26.8
M 155            CPS         2.9         80.5        7.44    36.65
M 198            CP4         2.8         91.8        7.52    21.85
M 200            CPS         2.84        82.8        7.42    25.85
112              CPI         2.65         86         7.25    33.05
233              CP2         2.64        88.9        7.53     4.34
311              CPS         2.94        85.1        7.56      43
323              CPS         2.9         94.6        7.67    26.75
343              CPS          2          95.8        7.48     28.5
403              CP4         2.89        82.8        7.32     16.1
407              CP4         2.99        88.4        7.91    18.95
611              CPS         2.52        88.4        7.87    52.85
222 D           CP 2         3.14        92.2        7.74    29.75
E-23             CP2         2.95        71.1                 58.5
E-29             CP2         2.69        83.3        7.12     79.4
E-47A            CPS         2.89        91.7        7.47      61
T 350A           CP2         2.9          82         6.88     82.2
T 395            CP2         2.65        80.1        7.92     4.88
T 400            CPI         2.1         89.6         7.9    19.73
611*A            CPS         2.52        87.6        7.85     52.4

                TSD          Ca        Mg         Na         K
Well         Calculated    (mg/L)    (mg/L)     (mg/L)    (mg/L)

M 104         20613.2       261.6     50.9       6250     1187.5
M 117A        17315.5       140.2     34.7       5575      1185
M 119A        34506.9       379.1     38.3       9850     2682.5
M 127         21457.8       288.3     31.1      6087.5    1352.5
M 133A        36545.4       819.4     53.9      10750      1840
M 148A        27327.2       379.1     28.7       8300     1652.5
M 155         38863.8       735.5     44.9      11325      2050
M 198         25601.2       187.6     56.9      7587.5    1577.5
M 200         22151.3       183.6     25.1      6312.5     1690
112           31216.4       533.1     47.9       9125     1992.5
233           23971.6       189.5     29.9       7275      1490
311           23548.8       288.3     23.9      7137.5    1432.5
323           27018.9       286.3     53.9      7437.5    1482.5
343           34309.7       789.8     59.9      9712.5    1447.5
403            15503        167.8     59.9       4540      1030
407           20042.6       108.6     35.9       5815      1185
611           33262.9       730.5     53.9       9875      1585
222 D         28719.9       296.2     27.5      8062.5     1907
E-23          39344.1       580.5     21.6      11175     2717.5
E-29          44491.3       785.8     14.4      13200     2962.5
E-47A         35074.6       647.6     16.8      10650     2247.5
T 350A        49251.2      1032.6     22.8     14312.5    2782.5
T 395         23924.2       337.6     27.5      7037.5     1420
T 400         12945.5       122.4     33.5       3910      560.5
611*A         33217.5       704.9     37.1     10087.5     1595

             HC[O.sub.3]     [Cl.sup.-]     S[O.sub.4]
Well           (mg/L)          (mg/L)         (mg/L)       % IB

M 104           41.5          11962.5          13.5        -2.8
M 117A          18.4            9375           0.09        3.3
M 119A          10.7           20150           0.09        -4.5
M 127           28.4          12812.5          5.9         -6.7
M 133A           19            22200           14.2        -5.6
M 148A          21.3           15875           7.1         -2.7
M 155           27.2           23800           0.09        -6.8
M 198           16.6           15250           0.09        -5.6
M 200           28.4           12950           0.09        -5.2
112             14.8           18250           0.09        -3.6
233             23.7           13700           3.6         -2.6
311             11.8           13550           5.4         -2.5
323             37.9           16750           0.09       -10.8
343             14.3           21587            29         -9.4
403             42.6            8950           4.4         -3.2
407              19           11887.5          8.5          -7
611              19            20100           22.4        -5.2
222 D            9.5           17300           0.09        -7.9
E-23            28.4           23750           0.09        -6.6
E-29            21.3           26600           0.09        -4.1
E-47A           16.6           20500           0.09        -2.1
T 350A          15.4           29700           0.09        -5.7
T 395            4.7          14062.5          6.3         -4.6
T 400           20.7            7570           33.6        -5.2
611*A           20.1           19875           22.5         -4

               B      Si[O.sub.2]
Well        (mg/L)       (mg/L)

M 104        13.26       851.7
M 117A       14.65       980.1
M 119A       17.42       1401.7
M 127        18.4        847.4
M 133A       19.86       853.9
M 148A       18.45       1074.3
M 155        27.92       890.2
M 198        28.8        915.9
M 200        22.35       973.7
112          15.8        1260.5
233          16.6        1262.6
311          21.9        1091.4
323           22         973.7
343          17.56       667.7
403          28.85       727.6
407          16.7        969.4
611          21.61       877.4
222 D        11.25       1114.9
E-23         24.2        1080.7
E-29         26.55        918
E-47A        21.97       995.1
T 350A       28.55       1393.1
T 395        10.25       1016.5
T 400        9.11        680.5
611*A         22         877.4

Table 2: Results of Na/K geothermometers.

            T [degrees]C     T [degrees]C
                 Na/K             Na/K
Well             [36]             [37]

M 104            264              279
M 117A           277              290
M 119A           306              318
M 127            282              295
M 133A           254              268
M 148A           269              283
M 155            259              274
M 198            274              288
M 200            304              316
112              280              293
233              273              286
311              270              284
323              270              284
343              240              255
403              284              297
407              272              286
611              247              262
222 D            289              302
E-23             292              305
E-29             283              296
E-47A            276              290
T 350A           267              281
T 395            271              285
T 400            236              251
611*A            246              260
COPYRIGHT 2018 Hindawi Limited
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2018 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research Article
Author:Morales-Arredondo, Ivan; Armienta, Maria Aurora; Segovia, Nuria
Geographic Code:0GULF
Date:Jan 1, 2018
Previous Article:The Empirical and Theoretical Miscible Characterization Method in Gas-Enhanced Oil Recovery.
Next Article:Near-Surface Geochemical Anomalies Integrated with Seismic and Well Data over the Contact of the Outer Carpathians and the Carpathian Foredeep (SE...

Terms of use | Privacy policy | Copyright © 2022 Farlex, Inc. | Feedback | For webmasters |