GIS-based hydrogeological databases and groundwater modelling

Reliability and validity of groundwater analysis strongly depend on the availability of large volumes of high-quality data. Putting all data into a coherent and logical structure supported by a computing environment helps ensure validity and availability and provides a powerful tool for hydrogeological studies. A hydrogeological geographic information system (GIS) database that offers facilities for groundwater-vulnerability analysis and hydrogeological modelling has been designed in Belgium for the Walloon region. Data from five river basins, chosen for their contrasting hydrogeological characteristics, have been included in the database, and a set of applications that have been developed now allow further advances. Interest is growing in the potential for integrating GIS technology and groundwater simulation models. A "loose-coupling" tool was created between the spatial-database scheme and the groundwater numerical model interface GMS (Groundwater Modelling System). Following time and spatial queries, the hydrogeological data stored in the database can be easily used within different groundwater numerical models. Electronic Publication     

Comparison of aquifer vulnerability assessment techniques. Application to the Néblon river basin (Belgium)

Five different methods for assessing intrinsic aquifer vulnerability were tested in a case study and their results compared. The test area was a slightly karstified district in the Condroz region of Belgium. The basin covers about 65 km² and the karst aquifer provides a water-supply of about 28,000 m³d^-1. The methods tested were: EPIK (Doerfliger et al. 1999), DRASTIC (Aller et al. 1987), 'German method' (von Hoyer and Söfner 1998), GOD (Foster 1987) and ISIS (Civita and De Regibus 1995). The results are compared and critically examined. From the analysis, it seems that reducing the number of parameters is unsatisfactory, due to the variety of geological conditions. The various methods produce very different results at any given site. As only physically-based methods can be checked for their reliability, it is clear that future vulnerability mapping techniques must incorporate such methods.     

Amplitude Preserving Kirchhoff Migration: A Traveltime Based Strategy

Amplitude versus offset information is a key feature to seismic reservoir characterization. Therefore amplitude preserving migration was developed to obtain this information from seismic reflection data. For complex 3-D media, however, this process is computationally expensive. In this paper we present an efficient traveltime based strategy for amplitude preserving migration of the Kirchhoff type. Its foundations are the generation of traveltime tables using a wavefront-oriented ray-tracing technique, and a generalized moveout relation for 3-D heterogeneous media. All required quantities for the amplitude preserving migration are computed from coarsely gridded traveltime tables. The migration includes the interpolation from the coarsely gridded input traveltimes onto the fine migration grid, the computation of amplitude preserving weight functions, and, optionally, the evaluation of an optimized migration aperture. Since ray tracing is employed for the traveltime computation the input velocity model needs to be smooth, i.e. velocity variations of spatial dimensions below the wavelength of the considered reflection signals are removed. Numerical examples on simple generic models validate the technique and an application to the Marmousi model demonstrates its potential to complex media. The major advantage of the traveltime based strategy consists of its computational efficiency by maintaining sufficient accuracy. Considerable savings in storage space (10⁵ and more for 3-D data with respect to no interpolation at all) can be achieved. The computational time for the stack can be substantially reduced (up to 90% in 3-D) with the optimized migration aperture since only those traces are stacked which really contribute to the image point under consideration.     

Anharmonicity of inner-OH stretching modes in hydrous phyllosilicates: assessment from first-principles frozen-phonon calculations

The theoretical anharmonicity of the inner-OH stretching modes of kaolinite and lizardite is determined using a frozen-phonon approach. The anharmonicity parameter, defined as the difference between half the first overtone frequency and the frequency of the fundamental transition, is −95 and −92 cm⁻¹ for kaolinite and lizardite, respectively. These values are consistent with the experimental measurements. The very good agreement usually observed between harmonic calculations of vibrational frequencies at the density functional theory (DFT) level, using the generalized gradient approximation (GGA), and experimental positions of vibrational bands results from the fortuitous, but almost exact, cancelation of two types of errors. One is related to the GGA approximation and the other to the neglect of anharmonicity.     

Sensitivity analysis for the EPIK method of vulnerability assessment in a small karstic aquifer, southern Belgium

Applying the EPIK parametric method, a vulnerability assessment has been made for a small karstic groundwater system in southern Belgium. The aquifer is a karstified limestone of Devonian age. A map of intrinsic vulnerability of the aquifer and of the local water-supply system shows three vulnerability areas. A parameter-balance study and a sensitivity analysis were performed to evaluate the influence of single parameters on aquifer-vulnerability assessment using the EPIK method. This approach provides a methodology for the evaluation of vulnerability mapping and for more reliable interpretation of vulnerability indices for karst groundwater resources. Received, March 1999/Revised, December 1999, February 2000/Accepted, February 2000     

Surface modes in the infrared spectrum of hydrous minerals: the OH stretching modes of bayerite

The theoretical infrared (IR) and Raman spectra of bayerite (β-Al(OH)₃) are computed in the density functional theory framework, using the linear response theory. The results are consistent with the occurrence of six non-equivalent OH groups in a bayerite structure with space group P2₁/n. Similar to gibbsite, the transmission powder IR spectrum of bayerite in the region of the OH stretching bands is found to depend on the shape of particles. In particular, the broadening of the strong band observed at about 3,460 cm⁻¹ in the spectrum of Al hydroxides is related to the electrostatic charges occurring at the surface of the polarized dielectric particles. The experimental correlation observed between the shape of this band and morphological parameters has therefore a physical, instead of chemical, origin.     

Effects of fluctuating daily surface fluxes on the time-mean oceanic circulation

The effect of fluctuating daily surface fluxes on the time-mean oceanic circulation is studied using an empirical flux model. The model produces fluctuating fluxes resulting from atmospheric variability and includes oceanic feedbacks on the fluxes. Numerical experiments were carried out by driving an ocean general circulation model with three different versions of the empirical model. It is found that fluctuating daily fluxes lead to an increase in the meridional overturning circulation (MOC) of the Atlantic of about 1 Sv and a decrease in the Antarctic circumpolar current (ACC) of about 32 Sv. The changes are approximately 7% of the MOC and 16% of the ACC obtained without fluctuating daily fluxes. The fluctuating fluxes change the intensity and the depth of vertical mixing. This, in turn, changes the density field and thus the circulation. Fluctuating buoyancy fluxes change the vertical mixing in a non-linear way: they tend to increase the convective mixing in mostly stable regions and to decrease the convective mixing in mostly unstable regions. The ACC changes are related to the enhanced mixing in the subtropical and the mid-latitude Southern Ocean and reduced mixing in the high-latitude Southern Ocean. The enhanced mixing is related to an increase in the frequency and the depth of convective events. As these events bring more dense water downward, the mixing changes lead to a reduction in meridional gradient of the depth-integrated density in the Southern Ocean and hence the strength of the ACC. The MOC changes are related to more subtle density changes. It is found that the vertical mixing in a latitudinal strip in the northern North Atlantic is more strongly enhanced due to fluctuating fluxes than the mixing in a latitudinal strip in the South Atlantic. This leads to an increase in the density difference between the two strips, which can be responsible for the increase in the Atlantic MOC.     

Assessment of shrinkage–swelling influences in coal seams using rank-dependent physical coal properties

Characterization of coal reservoirs and determination of in-situ physical coal properties related to transport mechanism are complicated due to having lack of standard procedures in the literature. By considering these difficulties, a new approach has been developed proposing the usage of relationships between coal rank and physical coal properties. In this study, effects of shrinkage and swelling (SS) on total methane recovery at CO₂ breakthrough (TMRB), which includes ten-year primary methane recovery and succeeding enhanced coalbed methane (ECBM) recovery up to CO₂ breakthrough, and CO₂ sequestration have been investigated by using rank-dependent coal properties. In addition to coal rank, different coal reservoir types, molar compositions of injected fluid, and parameters within the extended Palmer & Mansoori (P&M) permeability model were considered. As a result of this study, shrinkage and swelling lead to an increase in TMRB. Moreover, swelling increased CO₂ breakthrough time and decreased displacement ratio and CO₂ storage for all ranks of coal. Low-rank coals are affected more negatively than high-rank coals by swelling. Furthermore, it was realized that dry coal reservoirs are more influenced by swelling than others and saturated wet coals are more suitable for eliminating the negative effects of CO₂ injection. In addition, it was understood that it is possible to reduce swelling effect of CO₂ on cleat permeability by mixing it with N₂ before injection. However, an economical optimization is required for the selection of proper gas mixture. Finally, it is concluded from sensitivity analysis that elastic modulus is the most important parameter, except the initial cleat porosity, controlling SS in the extended P&M model by highly affecting TMRB.