Stochastic modelling of groundwater flow and solute transport in aquifers

This paper presents an introductory overview of recently developed stochastic theories for tackling spatial variability problems in predicting groundwater flow and solute transport. Advantages and limitations of the theories are discussed. Lastly, strategies based on the stochastic approaches to predict solute transport in aquifers are recommended.     

Isolation of resonance in acoustic backscatter from elastic targetsusing adaptive estimation schemes

It has been shown that at certain frequencies the acoustic backscatter from elastic targets exhibits certain resonance behavior which closely relates to the physical properties of the target such as dimension, thickness, and composition. The purpose of this paper is to develop an automated approach for identifying the presence of resonance in the acoustic backscatter from an unknown underwater target by isolating the resonance part from the specular contribution. An adaptive transversal filter structure is used to estimate the specular part of the backscatter and consequently the error signal would provide an estimate of the resonance part. An important aspect of this scheme lies in the fact that it does not require an underlying model for the elastic return. The adaptation rule is based upon fast recursive least squares (RLS) learning. The approach taken in this paper is general in the sense that it can be applied to targets of unknown geometry and thickness and, further, does not require any a priori information about the target and/or the environment. Test results on acoustic data are presented which indicate the effectiveness of the proposed approach     

Application of micropolar plasticity to post failure analysis in geomechanics

A micropolar elastoplastic model for soils is formulated and a series of finite element analyses are employed to demonstrate the use of a micropolar continuum in overcoming the numerical difficulties encountered in application of finite element method in standard Cauchy–Boltzmann continuum. Three examples of failure analysis involving a deep excavation, shallow foundation, and a retaining wall are presented. In all these cases, it is observed that the length scale introduced in the polar continuum regularizes the incremental boundary value problem and allows the numerical simulation to be continued until a clear collapse mechanism is achieved. The issue of grain size effect is also discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

A continuum‐discrete hydromechanical analysis of granular deposit liquefaction

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.     

Smeared crack approach: back to the original track

This paper briefly reviews the formulations used over the last 40 years for the solution of problems involving tensile cracking, with both the discrete and the smeared crack approaches. The paper focuses on the smeared approach, identifying as its main drawbacks the observed mesh‐size and mesh‐bias spurious dependence when the method is applied ‘straightly’. A simple isotropic local damage constitutive model is considered, and the (exponential) softening modulus is regularized according to the material fracture energy and the element size. The continuum and discrete mechanical problems corresponding to both the weak discontinuity (smeared cracks) and the strong discontinuity (discrete cracks) approaches are analysed and the question of propagation of the strain localization band (crack) is identified as the main difficulty to be overcome in the numerical procedure. A tracking technique is used to ensure stability of the solution, attaining the necessary convergence properties of the corresponding discrete finite element formulation. Numerical examples show that the formulation derived is stable and remarkably robust. As a consequence, the results obtained do not suffer from spurious mesh‐size or mesh‐bias dependence, comparing very favourably with those obtained with other fracture and continuum mechanics approaches. Copyright © 2006 John Wiley & Sons, Ltd.     

Isotope hydrological study of mean transit time in the granitic Strengbach catchment (Vosges massif,France): application of the FlowPC model with modified input function

Measurements of ¹⁸O concentrations in precipitation, soil solution, spring and runoff are used to determine water transit time in the small granitic Strengbach catchment (0·8 km²; 883–1146 m above sea level) located in the Vosges Mountains of northeastern France. Water transit times were calculated by applying the exponential, exponential piston and dispersion models of the FlowPC program to isotopic input (rainfall) and output (spring and stream water) data sets during the period 1989–95. The input function of the model was modified compared with the former version of the model and estimated by a deterministic approach based on a simplified hydrological balance. The fit between observed and calculated output data showed marked improvements compared with results obtained using the initial version of the model. An exponential piston version of the model applied to spring water indicates a 38·5 month mean transit time, which suggests that the volume in the aquifer, expressed in water depth, is 2·4 m. A considerable thickness (>45 m) of fractured bedrock may be involved for such a volume of water to be stored in the aquifer. Copyright © 2005 John Wiley & Sons, Ltd.     

Reef exploration in the East Sengkang Basin,Sulawesi, Indonesia

Reconnaissance seismic shot in 1971/72 showed a number of well defined seismic anomalies within the East Sengkang Basin which were interpreted as buried reefs. Subsequent fieldwork revealed that Upper Miocene reefs outcropped along the southern margin of the basin. A drilling programme in 1975 and 1976 proved the presence of shallow, gas-bearing, Upper Miocene reefs in the northern part of the basin. Seismic acquisition and drilling during 1981 confirmed the economic significance of these discoveries, with four separate accumulations containing about 750 × 10⁹ cubic feet of dry gas in place at an average depth of 700 m. Kampung Baru is the largest field and contains over half the total, both reservoir quality and gas deliverability are excellent. Deposition in the East Sengkang Basin probably started during the Early Miocene. A sequence of Lower Miocene mudstones and limestones unconformably overlies acoustic basement which consists of Eocene volcanics. During the tectonically active Middle Miocene, deposition was interrupted by two periods of deformation and erosion. Carbonate deposition became established in the Late Miocene with widespread development of platform limestones throughout the East Sengkang Basin. Thick pinnacle reef complexes developed in the areas where reef growth could keep pace with the relative rise in sea level. Most reef growth ceased at the end of the Miocene and subsequent renewed clastic sedimentation covered the irregular limestone surface. Late Pliocene regression culminated in the Holocene with erosion. The Walanae fault zone, part of a major regional sinistral strike-slip system, separates the East and West Sengkang Basins. Both normal and reverse faulting are inferred from seismic data and post Late Pliocene reverse faulting is seen in outcrop.     

Microstructures of co-axially folded vein segments and crenulation cleavage： Evidence for dissolution phenomena in the Chamba Thrust Sheet, western Himalayas

The metasediments in the Chamba region experienced three phases of deformation: DF1, DF2 and DF3.Folded quartz veins are co-folded with the F2 crenulation folds. Their geometric and tectonic significance is studied by microstructures and shortening adjacent to the discrete crenulation cleavage, S2. Microstructures of folded vein segments, their geometric changes and truncation to cleavage (S2) are mainly due to pressure-solution phenomena and the estimated volume loss from reconstructed vein segments range from 16 to 25.5%,which is closely related to volume decrease (26%) calculated from the polydeformed slates of North Wales areas.