Assessment of alluvial aquifer heterogeneity and development of stochastic hydrofacies models for the Hat Yai Basin in Southern Thailand

In previous studies, the groundwater flow models formulated for the Hat Yai Basin were conventional and deterministic because the geologic heterogeneity of the alluvial aquifer system in the basin had not yet been assessed. This paper describes an effort to develop hydrofacies models, such that the spatial variability of the aquifer system can be represented in a systematic way. Variogram parameters that characterize the alluvial aquifer heterogeneity were determined. Based on these variogram parameters, an indicator-based geostatistical approach was used to develop hydrofacies models using sequential indicator simulation. The hydrofacies models indicate three distinct aquifer units, namely Hat Yai, Khu Tao, and Kho Hong aquifers, which is in good agreement with a conceptual model, and incorporates spatial variability as observed in field data from borehole logs. The hydrofacies models can be used in groundwater modeling and simulations.     

Field characterization and data integration to define the hydraulic heterogeneity of a shallow granular aquifer at a sub-watershed scale

Providing a sound basis for aquifer management or remediation requires that hydrogeological investigations carried out to understand groundwater flow and contaminant transport be based on representative data that capture the heterogeneous spatial distribution of aquifer hydraulic properties. This paper describes a general workflow allowing the characterization of the heterogeneity of the hydraulic properties of granular aquifers at an intermediate scale of a few km². The workflow involves characterization and data integration steps that were applied on a 12-km² study area encompassing a decommissioned landfill emitting a leachate plume and its main surface water receptors. The sediments composing the aquifer were deposited in a littoral–sublittoral environment and show evidence of small-scale transitional heterogeneities. Cone penetrometer tests (CPT) combined with soil moisture and electrical resistivity (SMR) measurements were thus used to identify and characterize spatial heterogeneities in hydraulic properties over the study area. Site-specific statistical relationships were needed to infer hydrofacies units and to estimate hydraulic properties from high-resolution CPT/SMR soundings distributed all over the study area. A learning machine approach was used due to the complex statistical relationships between colocated hydraulic and CPT/SMR data covering the full range of aquifer materials. Application of this workflow allowed the identification of hydrofacies units and the estimation of horizontal hydraulic conductivity, vertical hydraulic conductivity and porosity over the study area. The paper describes and discusses data acquisition and integration methodologies that can be adapted to different field situations, while making the aquifer characterization process more time-efficient and less labor-intensive.     

Upscaling equivalent hydrofacies simulation based on borehole data generalization and transition probability geostatistics

The heterogeneity of hydrofacies is represented as spatial variability on different scales, and it has a significant impact on the behavior of groundwater flow and pollutant transport. However, effectively characterizing hydrofacies heterogeneity on different scales remains one of the most challenging problems in hydrogeology. In this study, an upscaling hydrofacies simulation (UHS) framework is proposed by integrating the upscaling borehole generalization (UBG) approach and transition probability geostatistics (TPG). A new UBG approach for generating virtual boreholes with equivalent hydrofacies information based on relatively high-density borehole lithological data is proposed, and the TPG is used to delineate the multiscale facies distribution. The results show that the UBG approach can significantly reduce borehole data volume while retaining the key equivalent hydrofacies information on a coarser scale. The UHS method can well characterize the overall distribution of equivalent hydrofacies on coarser scales, with the minor-component hydrofacies underestimated and the major-component hydrofacies overestimated to a lesser extent, and more equivalent facies appearing in strong heterogeneous areas. These results demonstrate that the UHS method can provide valuable capacity insights and advantages in characterizing hydrofacies heterogeneity on different scales using such high-density borehole lithological data.

     

Hierarchical dynamic stratigraphy in various Quaternary gravel deposits,Rhine glacier area (SW Germany): implications for hydrostratigraphy

The approach of dynamic stratigraphy aims to understand genetic processes that form stratigraphic units in a hierarchy of spatial and temporal scales. This approach was used to investigate Quaternary gravel deposits in terms of their sedimentology and in order to characterize the various sedimentary units in terms of their hydrogeological properties. Facies analysis within 62 gravel pits, laboratory permeability measurements of field samples and geophysical surveys (3-D georadar, 2-D seismic reflection) led to the detection and classification of sedimentary heterogeneity according to the following six scales whereby each scale can be translated into defined hydrostratigraphic units. (1) Particles and pores (micro scale) that reflect depositional and diagenetic fluid dynamics as well as source material behaviour (e.g. grain-size, roundness, lithological composition). This was found to be important for the hydrogeochemistry of groundwater in gravel aquifers (e.g. higher sorption capacity of carbon-rich limestone particles for organic pollutants). (2) Strata (meso scale) contain the recognition of sorting, fabric, texture and stratinomic features, which can give an indication of transport and depositional dynamics. Five major lithofacies groups, for example, were distinguished within fluvial gravel-bed deposits. Their variable hydraulic properties led to their subdivision into 12 hydrofacies types. They form the smallest mappable hydrostratigraphic units, which may result in either preferred pathways for fluid flow or flow barriers. (3) Depositional elements (macro scale) enable reconstruction of sedimentary/geomorphic elements and their dynamics within a depositional system (e.g. gravel-bed braided river systems are dominated by gravel sheet, gravel dunes and scour pool depositional elements). Hydrostratigraphically, the architecture of depositional elements influences the hydraulic connectivity and local permeability structure/distribution within an aquifer body. Five types of depositional elements in fluvial gravel-bed deposits were distinguished and their geometries/dimensions quantified. (4) Facies bodies (mega scale) composed of a stack of depositional elements and strata recording distinct environmental systems and their dynamics (e.g. a coarse-grained prograding delta system). Hydrostratigraphically, facies bodies represent major compartments of an aquifer. Six major types of meltwater-controlled facies bodies were identified in the study area. (5) Genetic sequences (mega scale) reflect the shifts of depositional environments caused by allocyclic changes (e.g. glacial advance recorded by a coarsening upward sequence) or autocyclic changes of landscape shaping events. These sequences may form separate hydrostratigraphic units or aquifer storeys. (6) Basin fill (giga scale) comprising the lateral and vertical stacking of facies bodies and genetic sequences controlled by either long-term glacier dynamics or short term flood events. The regional distribution of permeable gravel units and, for example, less permeable diamicts builds the larger scale hydrostratigraphy.     

Comparison of three geostatistical methods for hydrofacies simulation: a test on alluvial sediments

The hydrodispersive properties of porous sediments are strongly influenced by the heterogeneity at fine scales, which can be modeled by geostatistical simulations. In order to improve the assessment of the properties of three different geostatistical simulation methods (Sequential indicator simulation, SISIM; Transition probability geostatistical simulation, T-PROGS; Multiple point simulation, MPS) a comparison test at different scales was performed for a well-exposed aquifer analogue. In the analysed volume (approximately 30,000?m³) four operative hydrofacies have been recognised: very fine sand and silt, sand, gravelly sand and open framework gravel. Several equiprobable realizations were computed with SISIM, MPS and T-PROGS for a test volume of approximately 400?m³ and for the entire volume, and the different outcomes were compared with visual inspection and connectivity analysis of the very or poorly permeable structures. The comparison of the different simulations shows that the geological model is best reproduced when the simulations are realised separately for each highest rank depositional element and subsequently merged. Moreover, the three methods yield different images of the volume; in particular MPS is efficient in mapping the geometries of the most represented hydrofacies, whereas SISIM and T-PROGS can account for the distribution of the less represented facies.     

Assessing aquifer vulnerability to pollutants by electrical resistivity tomography (ERT) at a nitrate vulnerable zone in NE Spain

There is an increasing demand for groundwater vulnerability maps which illustrate the exposure of aquifers against pollution. These maps show areas of greatest potential for groundwater contamination on the basis of local subsurface conditions. Parameters affecting vulnerability are mainly permeability and thickness of each protective layer. For unconsolidated sediments, the permeability is strongly related to the clay content, which can be deduced from indirect resistivity methods, like electrical-imaging. Such geophysical methods can be of great help in groundwater vulnerability studies because they disturb neither the structure nor the dynamics of the soil. Sensibility analysis was performed of the electrical resistivity tomography method for accurately mapping soil media. Managers and public administrators may effectively use this method for assessing the potential risk of groundwater contamination. In the studied zone, electrical resistivity exhibits a wide range of variability that can be easily correlated to soil parameters, such as clay content and hydraulic conductivity. A numerical index of protection has been assessed from the geophysical information derived from 2D electrical resistivity tomography. This work represents a preliminary approach on the natural vulnerability evaluation of shallow aquifers at the Empordà basin (NE Spain) that is highly affected by diffuse pollution by nitrates.     

Stochastic groundwater flow analysis in the presence of trends in heterogeneous hydraulic conductivity fields

Due to changes in lithostatic pressure, differential fracturing across bedding planes and irregularities in depositional environments, hydraulic conductivity exhibits heterogeneities and trends at various spatial scales. Using spectral theory, we have examined the effect of trends in hydraulic conductivity on (1) the solution of the mean equation for hydraulic head, (2) the covariance of hydraulic head, (3) the cross-covariances of hydraulic head and log-hydraulic conductivity perturbations and their gradients, and (4) the effective hydraulic conductivity. It is shown that the field of hydraulic head is sensitive to the presence of trends in ways that cannot be predicted by the classical analysis based on stationary hydraulic conductivity fields. The controlling variables for the second moments of hydraulic head are the mean hydraulic gradient, the correlation scale of log-hydraulic conductivity and its variance, and the slope of the trend in log-hydraulic conductivity. The mean hydraulic gradient introduces complications in the analysis since it is, in general, spatially variable. In this respect, our results are approximate, yet indicative of the true role of spatially variable patterns of log-hydraulic conductivity on groundwater flow systems.     

Large to intermediate-scale aquifer heterogeneity in fine-grain dominated alluvial fans (Cenozoic As Pontes Basin,northwestern Spain): insight based on three-dimensional geostatistical reconstruction

Facies reconstructions are used in hydrogeology to improve the interpretation of aquifer permeability distribution. In the absence of sufficient data to define the heterogeneity due to geological processes, uncertainties in the distribution of aquifer hydrofacies and characteristics may appear. Geometric and geostatistical methods are used to understand and model aquifer hydrofacies distribution, providing models to improve comprehension and development of aquifers. However, these models require some input statistical parameters that can be difficult to infer from the study site. A three-dimensional reconstruction of a kilometer scale fine-grain dominated Cenozoic alluvial fan derived from more than 200 continuously cored, closely spaced, and regularly distributed wells is presented. The facies distributions were reconstructed using a genetic stratigraphic subdivision and a deterministic geostatistical algorithm. The reconstruction is only slightly affected by variations in the geostatistical input parameters because of the high-density data set. Analysis of the reconstruction allowed identification in the proximal to medial alluvial fan zones of several laterally extensive sand bodies with relatively higher permeability; these sand bodies were quantified in terms of volume, mean thickness, maximum area, and maximum equivalent diameter. These quantifications provide trends and geological scenarios for input statistical parameters to model aquifer systems in similar alluvial fan depositional settings.     

Estimation and scaling of hydrostratigraphic units: application of unsupervised machine learning and multivariate statistical techniques to hydrogeophysical data

Numerical models provide a way to evaluate groundwater systems, but determining the hydrostratigraphic units (HSUs) used in constructing these models remains subjective, nonunique, and uncertain. A three-step machine-learning approach is proposed in which fusion, estimation, and clustering operations are performed on different data sets to arrive at HSUs at different scales. In step one, data fusion is performed by training a self-organizing map (SOM) with sparse borehole hydrogeologic (lithology, hydraulic conductivity, aqueous field parameters, dissolved constituents) and geophysical (gamma, spontaneous potential, and resistivity) measurements. Estimation is handled by iterative least-squares minimization of the SOM quantization and topographical errors. Application of the Davies-Bouldin criteria to k-means clustering of SOM nodes is used to determine the number and location of discontinuous borehole HSUs with low lateral density (based on borehole spacing at 100 s m) and high vertical density (based on cm-scale logging). In step two, a scaling network is trained using the estimated borehole HSUs, airborne electromagnetic measurements, and numerically inverted resistivity profiles. In step three, independent airborne electromagnetic measurements are applied to the scaling network, and the estimation performed to arrive at a set of continuous HSUs with high lateral density (based on sounding locations at meter (m) spacing) and medium vertical density (based on m-layer modeled structure). Performance metrics are used to evaluate each step of the approach. Efficacy of the proposed approach is demonstrated to map local-to-regional scale HSUs using hydrogeophysical data collected at a heterogeneous surficial aquifer in northwestern Nebraska, USA.     

A multidisciplinary-based conceptual model of a fractured sedimentary bedrock aquitard: improved prediction of aquitard integrity

A hydrogeologic conceptual model that improves understanding of variability in aquitard integrity is presented for a fractured sedimentary bedrock unit in the Cambrian-Ordovician aquifer system of midcontinent North America. The model is derived from multiple studies on the siliciclastic St. Lawrence Formation and adjacent strata across a range of scales and geologic conditions. These studies employed multidisciplinary techniques including borehole flowmeter logging, high-resolution depth-discrete multilevel well monitoring, fracture stratigraphy, fluorescent dye tracing, and three-dimensional (3D) distribution of anthropogenic tracers regionally. The paper documents a bulk aquitard that is highly anisotropic because of poor connectivity of vertical fractures across matrix with low permeability, but with ubiquitous bed parallel partings. The partings provide high bulk horizontal hydraulic conductivity, analogous to aquifers in the system, while multiple preferential termination horizons of vertical fractures serve as discrete low vertical hydraulic conductivity intervals inhibiting vertical flow. The aquitard has substantial variability in its ability to protect underlying groundwater from contamination. Across widespread areas where the aquitard is deeply buried by younger bedrock, preferential termination horizons provide for high aquitard integrity (i.e. protection). Protection is diminished close to incised valleys where stress release and weathering has enhanced secondary pore development, including better connection of fractures across these horizons. These conditions, along with higher hydraulic head gradients in the same areas and more complex 3D flow where the aquitard is variably incised, allow for more substantial transport to deeper aquifers. The conceptual model likely applies to other fractured sedimentary bedrock aquitards within and outside of this region.     

Detection of the near surface structure through a multidisciplinary geophysical approach

The present survey aimed to image the subsurface structure, including karstic voids, and to evaluate the extent of the heterogeneities that can result in potentially dangerous collapse of road segments overlying these features. A multidisciplinary geophysical approach (seismic refraction, frequency domain electromagnetic and ground penetrating radar) in combination with a detailed geological survey indicated the presence of tectonic faults as well as velocity and conductivity anomalies along an old road within the area of Akrotiri at Chania (Crete). Due to the presence of subsurface fuel pipes, perpendicular to the direction of the road, 2D resistivity imaging was excluded from the applied geophysical methods. Interpretation of the geophysical data revealed that the section of the road investigated overlies prominent voids attributed mostly to karst features. The conductivity and velocity anomalies are interpreted to indicate an area where the host limestone rock has been downthrown by faulting and associated karstification. The continuation of this fault zone was observed in the slope of the road during later excavations. Interpretation, using geographic information systems (GIS) to integrate data, allowed these controls and relationships to be understood and monitored. The above methodology was proved successful for areas where the application of resistivity method is not possible.     

Digital simulation of the scale effect in hydraulic conductivity

 Measured hydraulic conductivity increases with the scale of testing, but the reason for this increase is not clear. Possibly, high-conductivity heterogeneities are more effective in raising hydraulic conductivity over the regional scale than at the local scale. Alternatively, borehole skin and storage effects, among others, can systematically bias the results of small-scale tests; thus, the increase may simply be an artifact of the test method. Radial-flow tests were simulated at various scales in digital models of fractured double-porosity media. The mean hydraulic conductivity increases until the radius of influence in the test exceeds the fracture spacing. Therefore, under radial flow, hydraulic conductivity is indeed dependent on measurement scale. The increase in hydraulic conductivity with scale is a natural consequence of heterogeneity. Over short distances, water converging toward a borehole must generally flow across heterogeneities. Therefore, small-scale tests tend to measure a weighted harmonic mean of the hydraulic-conductivity field. Over a large area, however, flow is primarily along high-conductivity heterogeneities. Therefore, large-scale tests approach a weighted arithmetic mean where high-conductivity heterogeneities have a greater influence. Received, April 1997 Revised, January 1998 Accepted, December 1997     

超高压变质岩物理性质的相关性对建立结晶岩区地球物理解释标尺的意义

大陆科学钻探为认识深部地壳的结构、组成、力学性质和物理状态提供了重要的基础数据。岩石物理性质的测量对于原位测井资料的校正和地表地球物理测量的解释具有重要意义,岩石物理性质和岩石学研究相结合,还能为岩石的变质作用以及变质过程中化学成分的迁移提供必要的约束。本文主要对CCSD100-3100m的综合岩石物理资料进行了相关性调查和聚类分析,并得出如下结论:(1)岩石物理性质中的地震波速度、密度、热导率之间具有强相关性,他们都受岩石中主要矿物的组成和含量所控制;(2)岩石的电阻率和磁化率受金属氧化物含量的控制,与岩性有相关性;孔隙度、渗透率与岩性没有相关性;(3)利用岩石的物理性质可以反演不同的岩性,其中地震波速度、密度和热导率对榴辉岩和片麻岩大类具有很好的分类效果,而结合磁化率和电阻率则能更好的区分出超基性岩、正片麻岩和副片麻岩。上述结论对综合地球物理解释中物理参数的选取和结晶岩区地球物理解释标尺的建立具有重要的意义。     

Coastal aquifer assessment based on geological and geophysical survey,northwestern Crete,Greece

Groundwater preservation comprises a major problem in water policy. The comprehension of the groundwater/hydraulic systems can provide the means to approach this problem. Generally, drilling is expensive and time-consuming. On the other hand, new techniques have been applied during the last few decades that provide useful information on the depth and quality of aquifers. Among them, transient electromagnetic method (TEM) is an appealing method that provides fast results with minimum field crew and solves several hydrogeological problems. Many portable systems for single-site measurements are commercially available. The TEM-Fast 48HPC was used for acquiring 106 soundings in the northwestern Crete in Greece for defining the hydrogeological characteristics of the study area, since there were no available data from boreholes. Detailed geological, hydrolithological and tectonic survey was applied prior to the geophysical measurements. All the data were integrated to produce a secure and reliable hydrogeological model for the study area prior to any future hydrowell. Specifically, geometrical and hydraulic data of the study area groundwater were acquired. Two unconnected aquifers were detected and their possible contamination due to saltwater intrusion was analyzed and eliminated. Moreover, a location for borehole construction and groundwater pumping based on the potential of the aquifer system was proposed. Finally, the contribution of TEM (and electrical resistivity tomography) geophysical methods in studying complex coastal aquifers is shown by this work.     

A combined study of radar facies,lithofacies and three‐dimensional architecture of an alpine alluvial fan (Illgraben fan,Switzerland)

Alluvial fans serve as useful archives that record the history of depositional and erosional processes in mountainous regions and thus can reveal the environmental controls that influenced their development. Economically, they play an important role as groundwater reservoirs as well as host rocks for hydrocarbons in deeply buried systems. The interpretation of these archives and the evaluation of their reservoir architecture, however, are problematic because marked heterogeneity in the distribution of sedimentary facies makes correlation difficult. This problem is compounded because the accumulated sedimentary deposits of modern unconsolidated fan systems tend to be poorly exposed and few such systems have been the focus of investigation using high‐resolution subsurface analytical techniques. To overcome this limitation of standard outcrop–analogue studies, a geophysical survey of an alpine alluvial fan was performed using ground‐penetrating radar to devise a scaled three‐dimensional subsurface model. Radar facies were classified and calibrated to lithofacies within a fan system that provided outcropping walls and these were used to derive a three‐dimensional model of the sedimentary architecture and identify evolutionary fan stages. The Illgraben fan in the Swiss Alps was selected as a case study and a network of ca 60 km sections of ground‐penetrating radar was surveyed. Seven radar facies types could be distinguished, which were grouped into debris‐flow deposits and stream‐flow deposits. Assemblages of these radar facies types show three depositional units, which are separated by continuous, fan‐wide reflectors; they were interpreted as palaeo‐surfaces corresponding to episodes of sediment starvation that affected the entire fan. An overall upward decline in the proportion of debris‐flow deposits from ca 50% to 15% and a corresponding increase in stream‐flow deposits were identified. The uppermost depositional unit is bounded at its base by a significant incision surface up to 700 m wide, which was subsequently filled up mostly by stream‐flow deposits. The pronounced palaeo‐surfaces and depositional trends suggest that allocyclic controls governed the evolution of the Illgraben fan, making this fan a valuable archive from which to reconstruct past sediment fluxes and environmental change in the Alps. The results of the integrated outcrop–geophysical approach encourage similar future studies on fans to retrieve their depositional history as well as their potential reservoir properties.     

Efficient upscaling of hydraulic conductivity in heterogeneous alluvial aquifers

An efficient method to upscale hydraulic conductivity (K) from detailed three-dimensional geostatistical models of hydrofacies heterogeneity to a coarser model grid is presented. Geologic heterogeneity of an alluvial fan system was characterized using transition-probability-based geostatistical simulations of hydrofacies distributions. For comparison of different hydrofacies architecture, two alternative models with different hydrofacies structures and geometries and a multi-Gaussian model, all with the same mean and variance in K, were created. Upscaling was performed on five realizations of each of the geostatistical models using the arithmetic and harmonic means of the K-values within vertical grid columns. The effects of upscaling on model domain equivalent K were investigated by means of steady-state flow simulations. A logarithmic increase in model domain equivalent K with increasing upscaling, was found for all fields. The shape of that upscaling function depended on the structure and geometry of the hydrofacies bodies. For different realizations of one geostatistical model, however, the upscaling function was the same. From the upscaling function a factor could be calculated to correct the upscaled K-fields for the local effects of upscaling.     

从一个布井失败案例探讨多种地球物理方法综合研究

为探讨多种地球物理方法结合的重要性,以一个应用地球物理方法布井失败案例为研究对象,通过比较高密度电阻率法、磁共振探测以及两种方法结合的探测结果与实际钻探资料的差异,分析了每种方法的优缺点,并借助经济成本计算,讨论了多种方法结合的可行性。结果表明:高密度电阻率法与磁共振探测各有优缺点。高密度电阻率法可很好地判断地层结构的整体变化,但对岩性的推断具有较强的不确定性。磁共振探测推断的地层岩性与钻孔资料基本吻合,但无法判断是否有断层或裂隙带存在。高密度电阻率法与磁共振探测结合推断的剖面尺度上地层岩性和厚度变化与钻孔资料基本一致,是揭示地层岩性及含水层位置的有效方法,说明应重视多种方法结合的研究及应用。相比不采用任何地球物理方法或单一方法,两种方法结合的费用更低,在实际应用中具有可行性。     

三角洲在坳陷盆地沉积中所占比例研究

在对全球现代三角洲体系统计分析的基础上,将坳陷盆地沉积体系简化为河流-三角洲-湖泊,建立沉积体系模型。后期沉积体系叠置在前期沉积体系之上,同时也对前期沉积体系进行改造。每期沉积体系体系均发育有河流-三角洲-湖泊,将其作为整体进行研究。假设盆地规模不变,沉积基准面下降-进积过程表现为河流沉积在盆地中所占比例增大,沉积基准面上升-退积过程表现为湖泊沉积在盆地中所占比例增大。变化过程中,三角洲在盆地中所占的比例变化较小;相对于沉积盆地,三角洲体系不占主体地位,而是处于一定的比例范围(三角洲总面积/盆地面积在0~10%之间)。因此三角洲的规模,取决于沉积盆地的规模,开展沉积相研究必须考虑宏观盆地尺度背景。     

Assessment of integrated electrical resistivity data on complex aquifer structures in NE Nuba Mountains – Sudan

Two geophysical techniques were integrated to map the groundwater aquifers on complex geological settings, in the crystalline basement terrain in northeast Nuba Mountains. The water flow is structurally controlled by the northwest–southeast extensional faults as one of several in-situ deformational patterns that are attributed to the collision of the Pan-African oceanic assemblage of the Nubian shield against the pre-Pan African continental crust to the west. The structural lineaments and drainage systems have been enhanced by the remote sensing technique. The geophysical techniques used are: vertical electrical soundings (VES) and electrical resistivity tomography (ERT), in addition to hydraulic conductivity measurements. These measurements were designed to overlap in order to improve the producibility of the geophysical data and to provide a better interpretation of the hydrogeological setting of the aquifer complex structure. Smooth and Block inversion schemes were attempted for the observed ERT data to study their reliability in mapping the different geometries in the complex subsurface. The VES data was conducted where ERT survey was not accessible, and inverted smoothly and merged with the ERT in the 3D resistivity grid. The hydraulic conductivity was measured for 42 water samples collected from the distributed dug wells in the study area; where extremely high saline zones were recorded and have been compared to the resistivity values in the 3D model.