A rapid method for hydraulic profiling in unconsolidated formations

Information on vertical variations in hydraulic conductivity ( K ) can often shed much light on how a contaminant will move in the subsurface. The direct-push injection logger has been developed to rapidly obtain such information in shallow unconsolidated settings. This small-diameter tool consists of a short screen located just behind a drive point. The tool is advanced into the subsurface while water is injected through the screen to keep it clear. Upon reaching a depth at which information about K is desired, advancement ceases and the injection rate and pressure are measured on the land surface. The rate and pressure values are used in a ratio that serves as a proxy for K . A vertical profile of this ratio can be transformed into a K profile through regressions with K estimates determined using other techniques. The viability of the approach was assessed at an extensively studied field site in eastern Germany. The assessment demonstrated that this tool can rapidly identify zones that may serve as conduits for or barriers to contaminant movement.     

Analysis of slug tests in formations of high hydraulic conductivity

A new procedure is presented for the analysis of slug tests performed in partially penetrating wells in formations of high hydraulic conductivity. This approach is a simple, spreadsheet-based implementation of existing models that can be used for analysis of tests from confined or unconfined aquifers. Field examples of tests exhibiting oscillatory and nonoscillatory behavior are used to illustrate the procedure and to compare results with estimates obtained using alternative approaches. The procedure is considerably simpler than recently proposed methods for this hydrogeologic setting. Although the simplifications required by the approach can introduce error into hydraulic-conductivity estimates, this additional error becomes negligible when appropriate measures are taken in the field. These measures are summarized in a set of practical field guidelines for slug tests in highly permeable aquifers.     

沉积盆地构造热演化研究进展：回顾与展望

构造热演化模拟是研究沉积盆地的重要手段之一，其模型依赖于沉积盆地的成因机制.裂谷盆地构造热演化的定量模型在描述盆地沉降和热流演化方面取得了极大的成功，实现了构造和热的完美结合.而前陆盆地的定量模型更多关注的是构造沉降，在构造与热的结合方面尚不够完善.关于克拉通盆地目前还没有很成熟的定量模型，构造热演化研究程度远远低于裂谷盆地和前陆盆地.随着我国陆域海相沉积盆地油气勘探的突破，对海相沉积盆地热体制的研究迫在眉睫.而我国陆域海相沉积盆地，如塔里木和四川盆地，演化历史长且复杂，是古生代海相克拉通与中、新生代前陆盆地组成的叠合盆地.现有的关于沉积盆地构造热演化的单一模式难以适应复杂的构造—热历史.对我国陆域海相大型沉积盆地进行深入全面的动力学分析，发展叠合盆地的构造—热演化模型，建立相应的构造热演化模式及模拟方法技术，将是一项具有开拓意义并极具挑战性的工作.     

Evaluation of streambed hydraulic conductivity heterogeneity in an urban watershed

It has long been understood that streambed hydraulic conductivity plays an important role in surface-subsurface solute exchange. Using a portable falling head permeameter in situ, we estimated the horizontal hydraulic conductivity, K, of the near-surface streambed sediments at a total of 85 locations encompassing two depth intervals: 7.5–10 and 10–12.5 cm. The measurements were conducted in an 80 m reach of Indian Creek, a small urban stream in Philadelphia, PA, USA. We found that the ln K data within each sediment layer were Gaussian, but the combined data set was not. The results indicated that while the mean hydraulic conductivity decreased with depth, horizontal heterogeneity (e.g. the variance) increased with depth. This strong contrast between layers suggests that they should be treated as separated entities in modeling studies. Variogram analyses across the stream suggested symmetry with respect to the thalweg in the upper layer and fractality in the lower layer. The variograms along the streams suggested that the K data are random.     

Probabilistic modeling of aquifer heterogeneity using reliability methods

A probabilistic model of groundwater contaminant transport is presented. The model is based on coupling first- and second-order reliability methods (FORM and SORM) with a two-dimensional finite element solution of groundwater transport equations. Uncertainty in aquifer media is considered by modeling hydraulic conductivity as a spatial random field with a prescribed correlation structure. FORM and SORM provide the probability that a contaminant exceeds a target level at a well, termed the probability of failure. Sensitivity of the probability of failure to basic variabilities in grid block conductivity is also obtained, at no additional computational effort. The effect of the choice of the predetermined target level at the observation well is provided, along with its effect on the relevant sensitivity information. Considerable saving in computational time was achieved by superimposing a coarse random variables mesh on a finer numerical mesh. The presence of regions of lower conductivity on the probabilistic event is analyzed, and the regions in which conductivity most affects the results are identified.     

Influence of heterogeneity on unsaturated hydraulic properties: (1) local heterogeneity and scale effect

Spatial heterogeneity is ubiquitous in nature, which may significantly affect the soil hydraulic property curves. The models of a closed‐form functional relationship of soil hydraulic property curves (e.g. VG model or exponential model) are valid at point or local scale based on a point‐scale hydrological process, but how do scale effects of heterogeneity have an influence on the parameters of these models when the models are used in a larger scale process? This paper uses a two‐dimensional variably saturated flow and solute transport finite element model (VSAFT2) to simulate variations of pressure and moisture content in the soil flume under a constant head boundary condition. By changing different numerical simulation block sizes, a quantitative evaluation of parameter variations in the VG model, resulting from the scale effects, is presented. Results show that the parameters of soil hydraulic properties are independent of scale in homogeneous media. Parameters of α and n in homogeneous media, which are estimated by using the unsaturated hydraulic conductivity curve (UHC) or the soil water retention curve (WRC), are identical. Variations of local heterogeneities strongly affect the soil hydraulic properties, and the scale affects the results of the parameter estimations when numerical experiments are conducted. Furthermore, the discrepancy of each curve becomes considerable when moisture content becomes closer to a dry situation. Parameters estimated by UHC are totally different from the ones estimated by WRC. Copyright © 2012 John Wiley & Sons, Ltd.     

Effects of sedimentary layer on earthquake source modeling from geodetic inversion

Ground deformation as observed with GPS or InSAR has been broadly inverted in constraining source parameter of earthquakes. However, for earthquakes occurring beneath sedimentary basins, the very slow sub-surface shear velocity (v S , down to 200 m/s) may cause substantial bias to earthquake source inversion if simple crustal models are used. For Bohai basin, Sichuan basin and rock-sites, we test effects of sub-surface shear velocity structure on ground deformation, and find that up to a factor of 2 overest...     

Influence of heterogeneity on unsaturated hydraulic properties (2) – percentage and shape of heterogeneity

In subsurface porous media, the soil water retention curve (WRC) and unsaturated hydraulic conductivity curve (UHC) are two important soil hydraulic property curves. Spatial heterogeneity is ubiquitous in nature, which may significantly affect soil hydraulic property curves. The main theme of this paper is to investigate how spatial heterogeneities, including their arrangements and amounts in soil flumes, affect soil hydraulic property curves. This paper uses a two‐dimensional variably saturated flow and solute transport finite element model to simulate variations of pressure and moisture content in soil flumes under a constant head boundary condition. To investigate the behavior of soil hydraulic property curves owing to variations of heterogeneities and their arrangements as well, cases with different proportions of heterogeneities are carried out. A quantitative evaluation of parameter variations in the van Genuchten model (VG model) resulting from heterogeneity is presented. Results show that the soil hydraulic properties are strongly affected by variations of heterogeneities and their arrangements. If the pressure head remains at a specific value, the soil moisture increases when heterogeneities increase in the soil flumes. On the other hand, the unsaturated hydraulic conductivity decreases when heterogeneities increase in the soil flumes under a constant pressure head. Moreover, results reveal that parameters estimated from both WRC and UHC also are affected by shapes of heterogeneity; this indicates that the parameters obtained from the WRC are not suitable for predicting the UHC of different shapes in heterogeneous media. Copyright © 2011 John Wiley & Sons, Ltd.     

Quantifying geological uncertainty for flow and transport modeling in multi-modal heterogeneous formations

In this work, we address the problem of characterizing the heterogeneity and uncertainty of hydraulic properties for complex geological settings. Hereby, we distinguish between two scales of heterogeneity, namely the hydrofacies structure and the intrafacies variability of the hydraulic properties. We employ multiple-point geostatistics to characterize the hydrofacies architecture. The multiple-point statistics are borrowed from a training image that is designed to reflect the prior geological conceptualization. The intrafacies variability of the hydraulic properties is represented using conventional two-point correlation methods, more precisely, spatial covariance models under a multi-Gaussian spatial law. We address the different levels and sources of uncertainty in characterizing the subsurface heterogeneity, and explore their effect on groundwater flow and transport predictions. Typically, uncertainty is assessed by way of many images, termed realizations, of a fixed statistical model. However, in many cases, sampling from a fixed stochastic model does not adequately represent the space of uncertainty. It neglects the uncertainty related to the selection of the stochastic model and the estimation of its input parameters. We acknowledge the uncertainty inherent in the definition of the prior conceptual model of aquifer architecture and in the estimation of global statistics, anisotropy, and correlation scales. Spatial bootstrap is used to assess the uncertainty of the unknown statistical parameters. As an illustrative example, we employ a synthetic field that represents a fluvial setting consisting of an interconnected network of channel sands embedded within finer-grained floodplain material. For this highly non-stationary setting we quantify the groundwater flow and transport model prediction uncertainty for various levels of hydrogeological uncertainty. Results indicate the importance of accurately describing the facies geometry, especially for transport predictions.     

Seismic failure modeling of concrete dams considering heterogeneity of concrete

Study on the failure process of high concrete dams subjected to strong earthquakes is crucial to reasonable evaluation of their seismic safety. Numerical simulation in this aspect involves dynamic failure analysis of big bulk concrete dam subjected to cyclic loading. The Rock Failure Process Analysis (RFPA) proposed by C.A. Tang, with successful applications to failure modeling of rock and concrete specimens mainly subjected to static loading, is extended for this purpose. For using the proposed model, no knowledge on the cracking route needs to be known beforehand, and no remeshing is required. Simulation of the whole process of elastic deformation, initiation and propagation of microcracks, severe damage and ultimate failure of concrete dams in earthquakes with a unified model is enabled. The model is verified through a shaking table test of an arch dam. Finally a practical gravity dam is employed as a numerical example. Considering the uncertainty in ground motion input and concrete material, typical failure process and failure modes of gravity dam are presented. Several small cracks may occur due to tension particularly at dam neck, dam faces and dam heel, and a few of them evolve into dominant ones. Relatively smaller earthquake may cause damage to the dam neck while a bigger one may bring on cracks at lower parts of the dams. Cracking at the dam bottom may incline to a direction almost perpendicular to the downstream face after propagating horizontally for a certain distance when the shaking is strong enough.     

A probabilistic framework for floodplain mapping using hydrological modeling and unsteady hydraulic modeling

ABSTRACT

Prediction of design hydrographs is key in floodplain mapping using hydraulic models, which are either steady state or unsteady. The former, which require only an input peak, substantially overestimate the volume of water entering the floodplain compared to the more realistic dynamic case simulated by the unsteady models that require the full hydrograph. Past efforts to account for the uncertainty of boundary conditions using unsteady hydraulic modeling have been based largely on a joint flood frequency–shape analysis, with only a very limited number of studies using hydrological modeling to produce the design hydrographs. This study therefore presents a generic probabilistic framework that couples a hydrological model with an unsteady hydraulic model to estimate the uncertainty of flood characteristics. The framework is demonstrated on the Swannanoa River watershed in North Carolina, USA. Given its flexibility, the framework can be applied to study other sources of uncertainty in other hydrological models and watersheds.     

Seismic failure modeling of concrete dams considering heterogeneity of concrete

Study on the failure process of high concrete dams subjected to strong earthquakes is crucial to reasonable evaluation of their seismic safety. Numerical simulation in this aspect involves dynamic failure analysis of big bulk concrete dam subjected to cyclic loading. The Rock Failure Process Analysis (RFPA) proposed by C.A. Tang, with successful applications to failure modeling of rock and concrete specimens mainly subjected to static loading, is extended for this purpose. For using the proposed model, no knowledge on the cracking route needs to be known beforehand, and no remeshing is required. Simulation of the whole process of elastic deformation, initiation and propagation of microcracks, severe damage and ultimate failure of concrete dams in earthquakes with a unified model is enabled. The model is verified through a shaking table test of an arch dam. Finally a practical gravity dam is employed as a numerical example. Considering the uncertainty in ground motion input and concrete material, typical failure process and failure modes of gravity dam are presented. Several small cracks may occur due to tension particularly at dam neck, dam faces and dam heel, and a few of them evolve into dominant ones. Relatively smaller earthquake may cause damage to the dam neck while a bigger one may bring on cracks at lower parts of the dams. Cracking at the dam bottom may incline to a direction almost perpendicular to the downstream face after propagating horizontally for a certain distance when the shaking is strong enough.     

用矢量有限元方法模拟随钻测井仪在倾斜各向异性地层中的电磁响应

用矢量有限元法模拟和分析了电各向异性倾斜沉积岩层的随钻测井电磁响应,为很好地模拟钻铤、井眼和激励线圈等,采用了基于圆柱坐标下的六面体剖分;同时采用了高阶叠层矢量基函数来描述每个剖分单元内的场,有利于减少剖分单元数和未知量.三维有限元程序采用数值模式匹配(当处于简单地层时)以及时域有限差分法计算的结果进行了验证,显示吻合非常好.最后仿真了不同各向异性系数和不同倾角的各向异性地层的随钻电阻率测井响应,有助于对各向异性及地层倾角的校正,以及对储层作出正确评价.     

Influence of inclination and permeability of solitary woody riparian plants on local hydraulic and sedimentary processes

The role of solitary woody riparian plants with respect to local erosion and deposition of sediments is investigated. A focus is laid on the characteristics ‘inclination’ and ‘permeability’ of the plant's projected frontal area. Therefore, two experimental studies using cylindrical obstacles were carried out in a laboratory flume, one aiming at inclination, the other at permeability. The first series revealed that the total amount of mobilized sediment around the cylinder on average decreased by 8–10% per 5° increasing inclination in streamwise direction. Locations of maximum scour depth simultaneously shifted downstream. A horseshoe vortex system, causing the frontal and lateral scouring, ceased to exist below inclinations of 25–30°. The second series revealed that with increasing permeability, frontal scour incision is delayed, and the eroded sediment volume is significantly reduced. With permeable obstacles, two system states were observed: first, frontal scouring with leeside deposition at higher flow velocities and, second, moderate leeside scouring at lower flow velocities. For up‐scaling and comparison, a field study focussing on fluvial obstacle marks at poplars and willows in secondary channels of the River Loire was additionally conducted. A modified analytical model enabled us to quantify the amount of deposited sediments leeside of the plants. Leeside sediment ridges are significantly stabilized and have a higher preservation potential when covered by pioneer vegetation. Under such conditions, they may indeed induce the development of stable islands. Eventually, ‘sediment ridge width’ turned out to be a suitable indicator for leeside deposited sediment volume, irrespective of spatial scale. Copyright © 2012 John Wiley & Sons, Ltd.     

Effect of input data in hydraulic modeling for flood warning systems

The main objective of this study was to quantify the error associated with input data, including various resolutions of elevation datasets and Manning’s roughness for travel time computation and floodplain mapping. This was accomplished on the test bed, the Grand River (Ohio, USA) using the HEC-RAS model. LiDAR data integrated with survey data provided conservative predictions, whereas coarser elevation datasets provided a positive difference in the travel time (11.03–15.01%) and inundation area (32.56–44.52%). The minimum differences in travel time and inundation area were 0.50–4.33% and 3.55–7.16%, respectively, when the result from LiDAR integrated with survey data was compared with a 10-m DEM integrated with survey data. The results suggest that a 10-m DEM in the channel and LiDAR data in the floodplain combined with survey data would be appropriate for a flood warning system. Additionally, Manning’s roughness of the channel section was found to be more sensitive than that of the floodplain. The decrease in inundation area was highest (8.97%) for the lower value of Manning’s roughness.     

Inverse modeling of coastal aquifers using tidal response and hydraulic tests

Remediation of contaminated aquifers demands a reliable characterization of hydraulic connectivity patterns. Hydraulic diffusivity is possibly the best indicator of connectivity. It can be derived using the tidal response method (TRM), which is based on fitting observations to a closed-form solution. Unfortunately, the conventional TRM assumes homogeneity. The objective of this study was to overcome this limitation and use tidal response to identify preferential flowpaths. Additionally, the procedure requires joint inversion with hydraulic test data. These provide further information on connectivity and are needed to resolve diffusivity into transmissivity and storage coefficient. Spatial variability is characterized using the regularized pilot points method. Actual application may be complicated by the need to filter tidal effects from the response to pumping and by the need to deal with different types of data, which we have addressed using maximum likelihood methods. Application to a contaminated artificial coastal fill leads to flowpaths that are consistent with the materials used during construction and to solute transport predictions that compare well with observations. We conclude that tidal response can be used to identify connectivity patterns. As such, it should be useful when designing measures to control sea water intrusion.     

Hydrochemical characteristics and groundwater evolution modeling in sedimentary rocks of the Tono mine,Japan

Based on geochemical data collected by Japan Nuclear Cycle Development Institute (JNC) in the Tono uranium mine, a conceptual groundwater evolution model developed by JNC is tested to evaluate whether equilibrium-based concepts of water–rock interaction are consistent with observed variations in the mineralogy and hydrochemistry of the Tono mine area. The chemical evolution of the groundwaters is modeled assuming local equilibrium for selected mineral–fluid reactions, taking into account the rainwater origin of these solutions. Results suggest that it is possible to interpret approximately the actual groundwater chemistry (i.e., pH, Eh, total dissolved concentrations of Si, Na, Ca, K, Al, carbonate and sulfate) if the following assumptions are adopted (a) CO₂ concentration in the gas phase contacting pore solutions in the overlying soil zone=10⁻¹ atm, and (b) minerals in the rock zone that control the solubility of respective elements in the groundwater include: chalcedony (Si), albite (Na), kaolinite (Al), calcite (Ca and carbonate), muscovite (K) and pyrite (Eh and sulfate). This result helps to build confidence in the use of simplified geochemical modeling techniques to develop an understanding of dominant geochemical reactions controlling groundwater chemistry in rocks similar to those that could be used for the geological disposal of radioactive wastes. It is noted, however, that the available field data may not be sufficient to adequately constrain parameters in the groundwater evolution model. In particular, more detailed information characterizing certain site properties are needed to improve the model. For this reason, a model that accounts for ion-exchange reactions among clay minerals, and which is based on the results of laboratory experiments, has also been evaluated in the present study. Further improvement of model considering ion-exchange reactions are needed in future, however.     

Inverse modeling of a radial multistep outflow experiment for determining unsaturated hydraulic properties

Modeling flow and solute transport in the unsaturated zone on the basis of the Richards equation requires specifying values for unsaturated hydraulic conductivity and water potential as a function of saturation. The objectives of the paper are to evaluate the design of a transient, radial, multi-step outflow experiment, and to determine unsaturated hydraulic parameters using inverse modeling. We conducted numerical simulations, sensitivity analyses, and synthetic data inversions to assess the suitability of the proposed experiment for concurrently estimating the parameters of interest. We calibrated different conceptual models against transient flow and pressure data from a multi-step, radial desaturation experiment to obtain estimates of absolute permeability, as well as the parameters of the relative permeability and capillary pressure functions. We discuss the differences in the estimated parameter values and illustrate the impact of the underlying model on the estimates. We demonstrate that a small error in absolute permeability, if determined in an independent experiment, leads to biased estimates of unsaturated hydraulic properties. Therefore, we perform a joint inversion of pressure and flow rate data for the simultaneous determination of permeability and retention parameters, and analyze the correlations between these parameters. We conclude that the proposed combination of a radial desaturation experiment and inverse modeling is suitable for simultaneously determining the unsaturated hydraulic properties of a single soil sample, and that the inverse modeling technique provides the opportunity to analyze data from nonstandard experimental designs.     

Geochemical modeling of evaporation processes on Mars: Insight from the sedimentary record at Meridiani Planum

New data returned from the Mars Exploration Rover (MER) mission have revealed abundant evaporites in the sedimentary record at Meridiani Planum. A working hypothesis for Meridiani evaporite formation involves the evaporation of fluids derived from the weathering of martian basalt and subsequent diagenesis. On Earth, evaporite formation in exclusively basaltic settings is rare. However, models of the evaporation of fluids derived from experimentally weathering synthetic martian basalt provide insight into possible formation mechanisms. The thermodynamic database assembled for this investigation includes both Fe²⁺ and Fe³⁺ in Pitzer's ion interaction equations to evaluate Fe redox disequilibrium at Meridiani Planum. Modeling results suggest that evaporation of acidic fluids derived from weathering olivine-bearing basalt should produce Mg, Ca, and Fe-sulfates such as jarosite and melanterite. Calculations that model diagenesis by fluid recharge predict the eventual breakdown of jarosite to goethite as well as the preservation of much of the initial soluble evaporite component at modeled porosity values appropriate for relevant depositional environments (< 0.30). While only one of several possible formation scenarios, this simple model is consistent with much of the chemical and mineralogical data obtained on Meridiani Planum outcrop.     

A semi-analytical model of expected areal-average infiltration under spatial heterogeneity of rainfall and soil saturated hydraulic conductivity

A semi-analytical model for the estimate of expected areal-average infiltration rate at hillslope scale is presented. It accounts for spatial heterogeneity of the saturated hydraulic conductivity, K_s, and rainfall rate, r. The K_s field is characterized by a lognormal probability density function while the rainfall rate r is represented by a uniform distribution between two extreme values. The model formulation relies upon the use of cumulative infiltration as the independent variable which is then expressed as a function of an expected time for use in practical applications. The solution is applicable for those ranges of r and K_s that allow for neglecting the infiltration of surface water running downslope into pervious soils (run-on process). The model was tested by comparisons with Monte Carlo simulations carried out for a variety of coefficients of variation of r and K_s over a clay loam soil and a sandy loam soil. The model was found to be very reliable both with coupled spatial variability of r and K_s and when only one variable is characterized by spatial heterogeneity while the other is uniform.