Velocities of compressional and shear waves in limestones

Carbonate rocks are important hydrocarbon reservoir rocks with complex textures and petrophysical properties (porosity and permeability) mainly resulting from various diagenetic processes (compaction, dissolution, precipitation, cementation, etc.). These complexities make prediction of reservoir characteristics (e.g. porosity and permeability) from their seismic properties very difficult. To explore the relationship between the seismic, petrophysical and geological properties, ultrasonic compressional‐ and shear‐wave velocity measurements were made under a simulated in situ condition of pressure (50 MPa hydrostatic effective pressure) at frequencies of approximately 0.85 MHz and 0.7 MHz, respectively, using a pulse‐echo method. The measurements were made both in vacuum‐dry and fully saturated conditions in oolitic limestones of the Great Oolite Formation of southern England. Some of the rocks were fully saturated with oil. The acoustic measurements were supplemented by porosity and permeability measurements, petrological and pore geometry studies of resin‐impregnated polished thin sections, X‐ray diffraction analyses and scanning electron microscope studies to investigate submicroscopic textures and micropores. It is shown that the compressional‐ and shear‐wave velocities (V_p and V_s, respectively) decrease with increasing porosity and that V_p decreases approximately twice as fast as V_s. The systematic differences in pore structures (e.g. the aspect ratio) of the limestones produce large residuals in the velocity versus porosity relationship. It is demonstrated that the velocity versus porosity relationship can be improved by removing the pore‐structure‐dependent variations from the residuals. The introduction of water into the pore space decreases the shear moduli of the rocks by about 2 GPa, suggesting that there exists a fluid/matrix interaction at grain contacts, which reduces the rigidity. The predicted Biot–Gassmann velocity values are greater than the measured velocity values due to the rock–fluid interaction. This is not accounted for in the Biot–Gassmann velocity models and velocity dispersion due to a local flow mechanism. The velocities predicted by the Raymer and time‐average relationships overestimated the measured velocities even more than the Biot model.     

Dissolution-filling mechanism of atmospheric precipitation controlled by both thermodynamics and kinetics

Affected by structural uplift,the Ordovician carbonate rockbed in the Tarim Basin,China,was exposed to dissolution and reformation of atmospheric precipitation many times,and formed a large quantity of karst caves serving as hydrocarbon reservoir.However,drilling in Tahe area showed that many large karst caves,small pores and fractures are filled by calcite,resulting in decrease in their reservoir ability.Calcite filled in the karst caves has very light oxygen isotopic composition and87Sr/86Sr ratio.Its 18OPDB ranges from 21.2‰to 13.3‰with the average of 16.3‰and its87Sr/86Sr ratio ranges from0.709561 to 0.710070 with the average of 0.709843.The isotope composition showed that calcite is related to atmospheric precipitation.Theoretic analyses indicated that the dissolving and filling actions of the precipitation on carbonate rocks are controlled by both thermodynamic and kinetic mechanisms.Among them,the thermodynamic factor determines that the precipitation during its flow from the earth surface downward plays important roles on carbonate rocks from dissolution to saturation,further sedimentation,and finally filling.In other words,the depth of the karstification development is not unrestricted,but limited by the precipitation beneath the earth surface.On the other hand,the kinetic factor controls the intensity,depth,and breadth of the karstification development,that is,the karstification is also affected by topographic,geomorphologic,climatic factors,the degree of fracture or fault,etc.Therefore,subject to their joint effects,the karstification of the precipitation on the Ordovician carbonate rocks occurs only within a certain depth(most about 200 m)under the unconformity surface,deeper than which carbonate minerals begin to sedimentate and fill the karst caves that were formed previously.     

CO2-induced dissolution of low permeability carbonates. Part II: Numerical modeling of experiments

We used the 3D continuum-scale reactive transport models to simulate eight core flood experiments for two different carbonate rocks. In these experiments the core samples were reacted with brines equilibrated with pCO₂ = 3, 2, 1, 0.5 MPa (Smith et al., 2013 [27]). The carbonate rocks were from specific Marly dolostone and Vuggy limestone flow units at the IEAGHG Weyburn-Midale CO₂ Monitoring and Storage Project in south-eastern Saskatchewan, Canada. Initial model porosity, permeability, mineral, and surface area distributions were constructed from micro tomography and microscopy characterization data. We constrained model reaction kinetics and porosity–permeability equations with the experimental data. The experimental data included time-dependent solution chemistry and differential pressure measured across the core, and the initial and final pore space and mineral distribution. Calibration of the model with the experimental data allowed investigation of effects of carbonate reactivity, flow velocity, effective permeability, and time on the development and consequences of stable and unstable dissolution fronts.The continuum scale model captured the evolution of distinct dissolution fronts that developed as a consequence of carbonate mineral dissolution and pore scale transport properties. The results show that initial heterogeneity and porosity contrast control the development of the dissolution fronts in these highly reactive systems. This finding is consistent with linear stability analysis and the known positive feedback between mineral dissolution and fluid flow in carbonate formations. Differences in the carbonate kinetic drivers resulting from the range of pCO₂ used in the experiments and the different proportions of more reactive calcite and less reactive dolomite contributed to the development of new pore space, but not to the type of dissolution fronts observed for the two different rock types. The development of the dissolution front was much more dependent on the physical heterogeneity of the carbonate rock. The observed stable dissolution fronts with small but visible dissolution fingers were a consequence of the clustering of a small percentage of larger pores in an otherwise homogeneous Marly dolostone. The observed wormholes in the heterogeneous Vuggy limestone initiated and developed in areas of greater porosity and permeability contrast, following pre-existing preferential flow paths.Model calibration of core flood experiments is one way to specifically constrain parameter input used for specific sites for larger scale simulations. Calibration of the governing rate equations and constants for Vuggy limestones showed that dissolution rate constants reasonably agree with published values. However the calcite dissolution rate constants fitted to the Marly dolostone experiments are much lower than those suggested by literature. The differences in fitted calcite rate constants between the two rock types reflect uncertainty associated with measured reactive surface area and appropriately scaling heterogeneous distribution of less abundant reactive minerals. Calibration of the power-law based porosity–permeability equations was sensitive to the overall heterogeneity of the cores. Stable dissolution fronts of the more homogeneous Marly dolostone could be fit with the exponent n = 3 consistent with the traditional Kozeny–Carman equation developed for porous sandstones. More impermeable and heterogeneous cores required larger n values (n = 6–8).     

Permeability Evolution During Non-linear Viscous Creep of Calcite Rocks

Permeability, storage capacity and volumetric strain were measured in situ during deformation of hot-pressed calcite aggregates containing 10, 20, and 30 wt% quartz. Both isostatic and conventional triaxial loading conditions were used. The tests were performed at confining pressure of 300 MPa, pore pressures between 50 to 290 MPa, temperatures from 673 to 873 K and strain rates of 3 × 10⁻⁵ s⁻¹. Argon gas was used as the pore fluid. The initial porosities of the starting samples varied from 5% to 9%, with higher porosity correlated to higher quartz content. Microstructural observations after the experiment indicate two kinds of pores are present: 1) Angular, crack-like pores along boundaries between quartz grains or between quartz and calcite grains and 2) equant and tubular voids within the calcite matrix. Under isostatic loading conditions, the compaction rate covaries with porosity and increases with increasing effective pressure. Most of the permeability reduction induced during compaction is irreversible and probably owes to plastic processes. As has been found in previous studies on hot-pressed calcite aggregates, permeability, k, is nonlinearly related to porosity, ϕ. Over small changes in porosity, the two parameters are approximately related as k ∝ ϕⁿ. The exponent n strongly increases as porosity decreases to a finite value (from about 4 to 6% depending on quartz content), suggesting a porosity percolation threshold. When subjected to triaxial deformation, the calcite-quartz aggregates exhibit shear-enhanced compaction, but permeability does not decrease as rapidly as it does under isostatic conditions. During triaxial compaction the exponent n only varies between 2 and 3. Non-isostatic deformation seems to reduce the percolation threshold, and, in fact, enhances the permeability relative to that at the same porosity during isostatic compaction. Our data provide constraints on the governing parameters of the compaction theory which describes fluid flow through a viscous matrix, and may have important implications for expulsion of sedimentary fluids, for fluid flow during deformation and metamorphism, and melt extraction from partially molten rocks.     

Effect of supercritical CO2 on carbonates: Savonnières sample case study

CO₂ geosequestration is an efficient way to reduce greenhouse gas emissions into the atmosphere. Carbonate rock formations are one of the possible targets for CO₂ sequestration due to their relative abundance and ability to serve as a natural trapping reservoir. The injected supercritical CO₂ can change properties of the reservoir rocks such as porosity, permeability, tortuosity, and specific surface area due to dissolution and precipitation processes. This, in turn, affects the reservoir characteristics, i.e., their elastic properties, storage capacity, stability, etc. The tremendous progresses made recently in both microcomputed X‐ray tomography and high‐performance computing make numerical simulation of physical processes on actual rock microstructures feasible. However, carbonate rocks with their extremely complex microstructure and the presence of microporosity that is below the resolution of microcomputed X‐ray tomography scanners require novel, quite specific image processing and numerical simulation approaches. In the current work, we studied the effects of supercritical CO₂ injection on microstructure and elastic properties of a Savonnières limestone. We used microtomographic images of two Savonnières samples, i.e., one in its natural state and one after injection and residence of supercritical CO₂. A statistical analysis of the microtomographic images showed that the injection of supercritical CO₂ led to an increase in porosity and changes of the microstructure, i.e., increase of the average volume of individual pores and decrease in the total number of pores. The CO₂ injection/residence also led to an increase in the mean radii of pore throats, an increase in the length of pore network segments, and made the orientation distribution of mesopores more isotropic. Numerical simulations showed that elastic moduli for the sample subjected to supercritical CO₂ injection/residence are lower than those for the intact sample.     

Reconstruction of the diagenesis of the fluvial-lacustrinedeltaic sandstones and its influence on the reservoir quality evolution

The reservoir quality of Jurassic and Triassic fluvial and lacustrine-deltaic sandstones of the Yanchang Oil Field in the Ordos Basin is strongly influenced by the burial history and facies-related diagenetic events. The fluvial sandstones have a higher average porosity (14.8%) and a higher permeability (12.7 × 10⁻³ μm²) than those of the deltaic sandstones (9.8% and 5.8 ×10⁻³) μm², respectively). The burial compaction, which resulted in 15% and 20% porosity loss for Jurassic and Triassic sandstones, respectively, is the main factor causing the loss of porosity both for the Jurassic and Triassic sandstones. Among the cements, carbonate is the main one that reduced the reservoir quality of the sandstones. The organic acidic fluid derived from organic matter in the source rocks, the inorganic fluid from rock-water reaction during the late diagenesis, and meteoric waters during the epidiagenesis resulted in the formation of dissolution porosity, which is the main reason for the enhancement of reservoir-quality.     

Impact of Poroelastic Response of Sandstones on Geothermal Power Production

During geothermal power production using a borehole doublet consisting of a production and injection well, the reservoir conditions such as permeability k, porosity φ and Skempton coefficient B at the geothermal research site Gross Schoenebeck/Germany will change. Besides a temperature decrease at the injection well and a change of the chemical equilibrium, also the pore pressure p _p will vary in a range of approximately 44 MPa ± 10 MPa in our reservoir at ?3850 to ?4258 m depth. This leads to a poroelastic response of the reservoir rocks depending on effective pressure p _eff (difference between mean stress and pore pressure), resulting in a change in permeability k, porosity φ and the poroelastic parameter Skempton coefficient B. Hence, we investigated the effective pressure dependency of Flechtinger sandstone, an outcropping equivalent of the reservoir rock via laboratory experiments. The permeability decreased by 21% at an effective pressure range from 3 to 30 MPa, the porosity decreased by 11% (p _eff = 6 to 65 MPa) and the Skempton coefficient decreased by 24% (p _eff = 4 to 25 MPa). We will show which mechanisms lead to the change of the mentioned hydraulic and poroelastic parameters and the influence of these changes on the productivity of the reservoir. The most significant changes occur at low effective pressures until 15 to 20 MPa. For our in situ reservoir conditions p _eff = 43 MPa a change of 10 MPa effective pressure will result in a change in matrix permeability of less than 4% and in matrix porosity of less than 2%. Besides natural fracture systems, fault zones and induced hydraulic fractures, the rock matrix its only one part of geothermal systems. All components can be influenced by pressure, temperature and chemical reactions. Therefore, the determined small poroelastic response of rock matrix does not significantly influence the sustainability of the geothermal reservoir.     

Effect of mineral reactions on the hydraulic properties of unsaturated soils: Model development and application

The selective radius shift model was used to relate changes in mineral volume due to precipitation/dissolution reactions to changes in hydraulic properties affecting flow in porous media. The model accounts for (i) precipitation/dissolution taking place only in the water-filled part of the pore space and further that (ii) the amount of mineral precipitation/dissolution within a pore depends on the local pore volume. The pore bundle concept was used to connect pore-scale changes to macroscopic soil hydraulic properties. Precipitation/dissolution induces changes in the pore radii of water-filled pores and, consequently, in the effective porosity. In a time step of the numerical model, mineral reactions lead to a discontinuous pore-size distribution because only the water-filled pores are affected. The pore-size distribution is converted back to a soil moisture characteristic function to which a new water retention curve is fitted under physically plausible constraints. The model equations were derived for the commonly used van Genuchten/Mualem hydraulic properties. Together with a mixed-form solution of Richards’ equation for aqueous phase flow, the model was implemented into the geochemical modelling framework PHREEQC, thereby making available PHREEQC’s comprehensive geochemical reactions. Example applications include kinetic halite dissolution and calcite precipitation as a consequence of cation exchange. These applications showed marked changes in the soil’s hydraulic properties due to mineral precipitation/dissolution and the dependency of these changes on water contents. The simulations also revealed the strong influence of the degree of saturation on the development of the saturated hydraulic conductivity through its quadratic dependency on the van Genuchten parameter α. Furthermore, it was shown that the unsaturated hydraulic conductivity at fixed reduced water content can even increase during precipitation due to changes in the pore-size distribution.     

Effect of pore geometry on Gassmann fluid substitution

Although there is no assumption of pore geometry in derivation of Gassmann's equation, the pore geometry is in close relation with hygroscopic water content and pore fluid communication between the micropores and the macropores. The hygroscopic water content in common reservoir rocks is small, and its effect on elastic properties is ignored in the Gassmann theory. However, the volume of hygroscopic water can be significant in shaly rocks or rocks made of fine particles; therefore, its effect on the elastic properties may be important. If the pore fluids in microspores cannot reach pressure equilibrium with the macropore system, assumption of the Gassmann theory is violated. Therefore, due to pore structure complexity, there may be a significant part of the pore fluids that do not satisfy the assumption of the Gassmann theory. We recommend that this part of pore fluids be accounted for within the solid rock frame and effective porosity be used in Gassmann's equation for fluid substitution. Integrated study of ultrasonic laboratory measurement data, petrographic data, mercury injection capillary pressure data, and nuclear magnetic resonance T₂ data confirms rationality of using effective porosity for Gassmann fluid substitution. The effective porosity for Gassmann's equation should be frequency dependent. Knowing the pore geometry, if an empirical correlation between frequency and the threshold pore‐throat radius or nuclear magnetic resonance T₂ could be set up, Gassmann's equation can be applicable to data measured at different frequencies. Without information of the pore geometry, the irreducible water saturation can be used to estimate the effective porosity.     

The role of gypsum and/or dolomite dissolution in tufa precipitation: lessons from the hydrochemistry of a carbonate–sulphate karst system

The precipitation of freshwater carbonates (tufa) along karstic rivers is enhanced by degassing of carbon dioxide (CO₂) downstream of karstic springs. However, in most karstic springs CO₂ degassing is not enough to force the precipitation of tufa sediments. Little is known about the role of dissolution of gypsum or dolomite in the hydrochemistry of these systems and how this affects the formation of tufa deposits. Here we present a monitoring study conducted over a year in Trabaque River (Spain). The river has typical karst hydrological dynamics with water sinking upstream and re‐emerging downstream of the canyon. Mixing of calcium–magnesium bicarbonate and calcium sulphate waters downstream of the sink enhances the dissolution of carbonates and potentially plays a positive role in the formation of tufa sediments. However, due to the common‐ion effect, dissolution of dolomite and/or gypsum causes precipitation of underground calcite cements as part of the incongruent dissolution of dolomite/dedolomitization process, which limits the precipitation of tufa sediments. Current precipitation of tufa is scant compared to previous Holocene tufa deposits, which likely precipitated from solutions with higher saturation indexes of calcite (SIcc values) than nowadays. Limited incongruent dissolution of dolomite/dedolomitization favours higher SIcc values. This circumstance occurs when waters with relatively high supersaturation of dolomite and low SO₄^2? composition sink in the upper sector of the canyon. In such a scenario, the process of mixing waters enhances the exclusive dissolution of limestones, preventing the precipitation of calcite within the aquifer and favouring the increase of SIcc values downstream of the springs. Such conditions were recorded during periods of high water level of the aquifers and during floods. This research shows that the common‐ion effect caused by the dissolution of gypsum and/or dolomite rocks can limit [or favour] the precipitation of tufa sediments depending on the occurrence [or not] of incongruent dissolution of dolomite/dedolomitization. Copyright © 2016 John Wiley & Sons, Ltd.     

The porosity trend and pore sizes of the rocks in the continental crust of the earth: Evidence from experimental data on permeability

The depth trends of permeability are constructed from the measurements of the tight rocks typical of the basement of the continental crust at temperatures up to 600°C and pressures up to 200 MPa. It is established that the permeability decreases with depth. The statistical processing of the experimental data yielded the generalized dependence logk = ?12.6-3.23H ^0.223. The method is suggested and, based on the experimental data on permeability, the estimates are obtained for the effective diameters of the pore channels and effective porosity at the PT parameters corresponding to the in situ deep zones of the continental crust. It is found that porosity decreases with depth, while distinct depth dependence of the pore sizes is not observed. The dependence of porosity on the depth is approximated by the relationship logφ = ?0.65–0.1H + 0.0019H ². The porosity is estimated at a few percent for a depth of 10 km with a decline to 0.01–0.1% at 35 km. The estimates of porosity retrieved from the experimental data agree with the theoretical calculations based on the present-day ideas of the structure of the discrete media and with the results of magnetotelluric sounding. Thus, according to three independent estimates, the porosity of the rocks of the continental crust decreases with depth. At the same time, in both the intermediate and lower crust there are intervals where the porosity values suggest the presence of fluid-saturated horizons at these depths.     

Geochemistry and isotopes of fluids from sulphur springs, Valles Caldera, New Mexico

Detailed geochemistry supported by geologic mapping has been used to investigate Sulphur Springs, an acid-sulfate hot spring system that issues from the western flank of the resurgent dome inside Valles Caldera. The most intense activity occurs at the intersection of faults offsetting caldera-fill deposits and post-caldera rhyolites. Three geothermal wells in the area have encountered pressures <1 MPa and temperatures of 200°C at depths of 600 to 1000 m. Hot spring and fumarole fluids may discharge at boiling temperatures with pH 1.0 and SO₄ 8000 mg/l. These conditions cause argillic alterations throughout a large area.Non-condensible gases consist of roughly 99% CO₂ with minor amounts of H₂S, H₂, and CH₄. Empirical gas geothermometry suggests a deep reservoir temperature of 215 to 280°C. Comparison of ¹³C and ¹⁸O between CaCO₃ from well cuttings and CO₂ from fumarole steam indicates a fractionation temperature between 200 and 300°C by decarbonation of hydrothermally altered Paleozoic limestone and vein calcite in the reservoir rocks. Tritium concentrations obtained from steam condensed in a mudpot and deep reservoir fluids (Baca #13, 278°C) are 2.1 and 1.0 T.U. respectively, suggesting the steam originates from a reservoir whose water is mostly >50 yrs old. Deuterium contents of fumarole steam, deep reservoir fluid, and local meteoric water are practically identical even though ¹⁸O contents range through 4‰, thus, precipitation on the resurgent dome of the caldera could recharge the hydrothermal system by slow percolation. From analysis of D and ¹⁸O values between fumarol steam and deep reservoir fluid, steam reaches the surface either (1) by vaporizing relatively shallow groundwater at 200°C or (2) by means of a two-stage boiling process through an intermediate level reservoir at roughly 200°C.Although many characteristics of known vapor-dominated geothermal systems are found at Sulphur Springs, fundamental differences exist in temperature and pressure of our postulated vapor-zone. We propose that the reservoir beneath Sulphur Springs is too small or too poorly confined to sustain a “true” vapor-dominated system and that the Sulphur Springs system may be a “dying” vapor-dominated system that has practically boiled itself dry.     

Aragonite dissolution on the Bermuda pedestal: Its depth and geochemical significance

The depth distribution of pteropod and planktonic foram tests, and fine-grained (<62 μm) aragonite, high-Mg calcite (12 mode mol.% MgCO₃), and low-Mg calcite has been determined for surface sediments of an area of the eastern slope of the Bermuda pedestal. Over the range 1800–3000 m, fine-grained aragonite and fine-grained high-Mg calcite gradually disappear relative to fine-grained low-Mg calcite, and pteropods gradually disappear relative to planktonic forams. This is interpreted as preferential dissolution of aragonite (and high-Mg calcite) relative to low-Mg calcite over this depth range. Coarse aragonitic debris derived from shallow-water organisms living on the Bermuda platform does not show consistent disappearance over the same depth zone. Chemical analyses of bottom water samples taken at the same time as some of the sediment samples indicate that the degree of saturation with respect to aragonite Ω_A over the zone of aragonite disappearance ranges from 0.55 to 0.85; i.e. major dissolution occurs only at Ω_A values distinctly less than one. These results lend credence to the hypothesis that CaCO₃ dissolution in the oceans, both as aragonite and as calcite, takes place mainly as a response to complex chemical kinetic phenomena and not as a result of the simple attainment of undersaturation (thermodynamic hypothesis) or the resuspension of bottom sediment (hydrodynamic hypothesis).     

Geochemical evolution of groundwater in the carbonate aquifers of Sierra de Segura (Betic Cordillera, southern Spain)

Sierra de Segura (Betic Cordillera), with a total area of over 3000 km², is the source of the two principal rivers in southern Spain, the Guadalquivir and the Segura. Due to the orographic effect of these mountains, precipitations are considerably more abundant than in nearby lowland areas, where the climate is semi-arid. Sierra de Segura is constituted of Mesozoic and Cenozoic sedimentary rocks, among which there are thick limestone–dolomitic formations which have given rise to extensive outcrops of permeable materials. In geomorphological terms, there is a large plateau intensively karstified that constitutes the main recharge area. Discharge takes place via a large number of springs, of which the 50 most important add up to a mean spring flow of about 13,500 l/s. The active geochemical processes in aquifers of Sierra de Segura, with their corresponding time sequence, are: dissolution of CO₂, dissolution of calcite, incongruent dissolution of dolomite, dedolomitization, exsolution of CO₂, and precipitation of calcite. More evolved water has higher temperature, magnesium content and Mg/Ca ratio; therefore, these parameters can be utilised as indicators of the degree of hydrochemical evolution. In addition, a good correlation between water temperature and magnesium concentration (or Mg/Ca ratio) indicates that an increase in temperature accelerates the kinetics of the dissolution of dolomite. Finally, the distribution of the temperatures in the vadose zone, determined by atmospheric thermal gradient, implies an apparent stratification of the predominant hydrochemical processes and of the groundwater physical and chemical characteristics.     

δCa evolution in a carbonate aquifer and its bearing on the equilibrium isotope fractionation factor for calcite

To improve understanding of Ca isotope transport during water-rock interaction on the continents, we measured dissolved δ⁴⁴Ca values along a 236 km flow path in the Madison aquifer, South Dakota, where fluids have chemically evolved according to dolomite and anhydrite dissolution, calcite precipitation, and Ca-for-Na ion-exchange over a timescale spanning ~ 15 kyr. We used a reactive transport model employing rate data constrained from major ion mass-balances to evaluate the extent to which calcite precipitation and ion-exchange fractionate Ca isotopes. Elevated δ⁴⁴Ca values during the initial and final stages of water transport possibly result from calcite precipitation under supersaturated conditions and Ca-for-Na ion-exchange, respectively. However, for the bulk of the flow path, δ⁴⁴Ca values evolve by mixing between anhydrite and dolomite dissolution, with no fractionation during calcite precipitation under saturated conditions. We attribute the absence of Ca isotope fractionation to the long timescale of water-rock interaction and slow rate of calcite precipitation, which have enabled fluids to chemically and isotopically equilibrate with calcite. We therefore conclude that the equilibrium Ca isotope fractionation factor between calcite and water (Δ_cal–w) is very close to zero. To the extent that the Madison aquifer typifies other groundwater systems where calcite slowly precipitates from solutions at or near chemical equilibrium, this study suggests that groundwater contributions to δ⁴⁴Ca variability on the continents can be modeled according to simple mixing theory without invoking isotope discrimination.     

NMR Measurements in Carbonate Rocks: Problems and an Approach to a Solution

Carbonate rocks are well known for their complex petrophysical behavior where, in contrast to siliciclastic rocks, different parameters, including porosity and permeability, usually are not directly related. This behavior is the result of thorough reorganization of porosity during diagenesis, and it turns prediction of reservoir quality of carbonate rocks into a challenge. The study presented here deals with the problem of utilizing NMR techniques in prediction of petrophysical properties in carbonates.We employ a visual porosity classification as a priori knowledge for better interpreting NMR data for prediction purposes. This allows for choice of suitable T₂ cutoff values to differentiate movable from bound fluids adapted for the specific carbonate rock, thus resulting in better interpretation of NMR data. The approach of using a genetic pore type classification for adapting the conventional method for T₂ cutoff determination, which originally was developed for siliciclastic rocks, is promising. Similarly, for permeability determination on the basis of NMR measurements, the classification of carbonate rocks based on porosity types also shows potential. The approach implemented here has the promise to provide a basis of standardized interpretation of NMR data from carbonate rocks.Acknowledgment We are grateful to Baker Hughes INTEQ, Celle, for permission to publish the results of this work. This study was part of I.S.s Masters Thesis at Hannover University that was carried out in collaboration with Baker Hughes INTEQ, Celle, Germany. We appreciate comments of an anonymous referee and PAGEOPH editor Brian J. Mitchell.     

Compressional- and shear-wave velocities as a function of confining stress in dry sandstones

A laboratory study was carried out to investigate the influence of confining stress on compressional- and shear-wave velocities for a set of rock samples from gas-producing sandstone reservoirs in the Cooper Basin, South Australia. The suite of samples consists of 22 consolidated sublitharenites with helium porosity ranging from 2.6% to 16.6%. We used a pulse-echo technique to measure compressional- and shear-wave velocities on dry samples (cylindrical 4.6 × 2 cm) at room temperature and at elevated confining stress (≤ 60 MPa). Compressional- and shear-wave velocities in samples increase non-linearly with confining stress. A regression equation of the form V = A ? Be^?DP gives a good fit to the measured velocities with improved prediction of velocities at high confining stresses compared with equations suggested by other studies. The predicted microcrack-closure stresses of the samples show values ranging from 70 MPa to 95 MPa and insignificant correlation with porosity, permeability or clay content. There is a positive correlation between change in velocity with core porosity and permeability, but this association is weak and diminishes with increasing confining stress. Experimental results show that pore geometry, grain-contact type, and distribution and location of clay particles may be more significant than total porosity and clay content in describing the stress sensitivity of sandstones at in situ reservoir effective stress. The stress dependence of Cooper Basin sandstones is very large compared with data from other studies. The implication of our study for hydrocarbon exploration is that where the in situ reservoir effective stress is much less than the microcrack-closure stress of the reservoir rocks, the variation of reservoir effective stress could cause significant changes in velocity of the reservoir rocks. The velocity changes induced by effective stress in highly stress-sensitive rocks can be detected at sonic-log and probably surface-seismic frequencies.     

Hydrogeochemistry of seepage water collected within the Youngcheon diversion tunnel,Korea: source and evolution of SO4‐rich groundwater in sedimentary terrain

In the Youngcheon Diversion Tunnel area, South Korea, 46 samples of tunnel seepage water (TSW) and borehole groundwater were collected from areas with sedimentary rocks (mainly sandstone and shale) and were examined for hydrogeochemical characteristics. The measured SO₄ concentrations range widely from 7·7 to 942·0 mg/l, and exceed the Korean Drinking Water Standard (200 mg/l) in about half the samples. The TDS (total dissolved solid) content generally is high (171–1461 mg/l) from more shale‐rich formations and also reflects varying degrees of water–rock interaction. The water is classified into three groups: Ca? SO₄ type (61% of the samples collected), Ca? SO₄? HCO₃ type (15%) and Ca? HCO₃ type (24%). The Ca? HCO₃ type water (mean concentrations=369 mg/l Ca, 148 mg/l HCO₃ and 23 mg/l SO₄) reflected the simple reaction between CO₂‐recharged water and calcite, whereas the more SO₄‐rich nature of Ca? SO₄ type water (mean concentrations=153 mg/l Ca, 66 mg/l HCO₃ and 416 mg/l SO₄) reflected the oxidation of pyrite in sedimentary rocks and fracture zones. Pyrite oxidation resulted in precipitation of amorphous iron hydroxide locally within the tunnel as well as in high concentrations of Ca (mean 153 mg/l) and Na (mean 49 mg/l) for TSW, and is associated with calcite dissolution resulting in pH buffering. The pyrite oxidation required for the formation of Ca? SO₄ type water was enhanced by the diffusion of oxygenated air through the fractures related to the tunnel's construction. The subsequent outgassing of CO₂ into the tunnel resulted in precipitation of iron‐bearing carbonate. Copyright © 2001 John Wiley & Sons, Ltd.     

Exploration drilling and reservoir model of the Platanares geothermal system, Honduras, Central America

Results of drilling, logging, and testing of three exploration core holes, combined with results of geologic and hydrogeochemical investigations, have been used to present a reservoir model of the Platanares geothermal system, Honduras. Geothermal fluids circulate at depths ≥ 1.5 km in a region of active tectonism devoid of Quaternary volcanism. Large, artesian water entries of 160 to 165°C geothermal fluid in two core holes at 625 to 644 m and 460 to 635 m depth have maximum flow rates of roughly 355 and 560 l/min, respectively, which are equivalent to power outputs of about 3.1 and 5.1 MW(thermal). Dilute, alkali-chloride reservoir fluids (TDS ≤ 1200 mg/kg) are produced from fractured Miocene andesite and Cretaceous to Eocene redbeds that are hydrothermally altered. Fracture permeabillity in producing horizons is locally greater than 1500 and bulk porosity is ≤ 6%. A simple, fracture-dominated, volume-impedance model assuming turbulent flow indicates that the calculated reservoir storage capacity of each flowing hole is approximately 9.7 × 10⁶ l/(kg cm⁻²), Tritium data indicate a mean residence time of 450 yr for water in the reservoir. Multiplying the natural fluid discharge rate by the mean residence time gives an estimated water volume of the Platanares system of ≥ 0.78 km³. Downward continuation of a 139°C/km “conductive” gradient at a depth of 400 m in a third core hole implies that the depth to a 225°C source reservoir (predicted from chemical geothermometers) is at least 1.5 km. Uranium-thorium disequilibrium ages on calcite veins at the surface and in the core holes indicate that the present Platanares hydrothermal system has been active for the last 0.25 m.y.