Numerical estimation of REV and permeability tensor for fractured rock masses by composite element method

The Monte Carlo method is used to generate parent stochastic discrete fracture network, from which a series of fractured rock samples of different sizes and orientations are extracted. The fracture network combined with a regular grid forms composite element mesh of the fractured rock sample, in which each composite element is composed of sub‐elements incised by fracture segments. The composite element method (CEM) for the seepage is implemented to obtain the nodal hydraulic potential as well as the seepage flow rates through the fractured rock samples. The application of CEM enables a large quantity of stochastic tests for the fractured rock samples because the pre‐process is facilitated greatly. By changing the sizes and orientations of the samples, the analysis of the seepage characteristics is realized to evaluate the variation of the permeability components, the existence of the permeability tensor and the representative element volume. The feasibility and effectiveness are illustrated in a numerical example. Copyright © 2008 John Wiley & Sons, Ltd.     

Biaxial test simulations using a packing of polygonal particles

The mechanical response of cohesionless granular materials under monotonic loading is studied by performing molecular dynamic simulations. The diversity of shapes of soil grains is modelled by using randomly generated convex polygons as granular particles. Results of the biaxial test obtained for dense and loose media show that samples achieve the same void ratio at large strains independent of their initial density state. This limit state resembles the so‐called critical state of soil mechanics, except for some stress fluctuations, which remain for large deformations. These fluctuations are studied at the micro‐mechanical level, by following the evolution of the co‐ordination number, force chains and the fraction of the sliding contacts of the sample. Copyright © 2007 John Wiley & Sons, Ltd.     

Geostatistical estimation of chemical contamination in stream sediments: The case study of Vale das Gatas mine (northern Portugal)

Based on an environmental geochemistry case study carried out in the neighbourhood of a W–Sn abandoned mine, the pollution in stream sediments was modelled through a Global Contamination Index. Such an index permits one to summarize the combination of deleterious elements in a single variable, obtained by the projection of samples onto the first axis of a PCASD (Principal Components Analysis of Standardized Data) applied to the entire n × p matrix containing the available concentrations of p = 16 elements in the set of n = 220 collected samples.In order to provide a sound basis for a coherent planning of the remediation process which will be put in operation in the affected area, it is necessary to balance the costs of reclaiming with the probabilities of exceeding the upper limits accepted for concentrations of environmentally harmful elements in sediments. Given these limits, they are back-transformed in the index values, providing a practical threshold between ‘clean’ and ‘contaminated’ samples. On the other hand, the minimum dimension of the cell to be reclaimed is restrained by the selected remediation process to be applied in the affected area. Hence, to meet the constraints of such a remediation process, it is required to estimate the probabilities of exceeding the index threshold in technologically meaningful sub-areas. For this end, the Indicator Block Kriging technique was applied, producing a series of maps where sub-areas to be reclaimed can be spotted for different probability levels. These maps, on which the decision making remediation agency can rely for its cost-benefit analysis, take into account both the spatial structure of ‘clean’ vs. ‘contaminated’ samples and the constraints of the reclaiming process.     

Methodology of hydrogeological characterization of deep carbonate aquifers as potential reservoirs of groundwater. Case of study: the Jurassic aquifer of <Emphasis Type="Italic">El Maestrazgo</Emphasis> (Castellón,Spain)

A methodology for the characterization of deep carbonate aquifers has been developed and applied to El Maestrazgo Jurassic aquifer in Castellón, Spain. Characterization of these aquifer formations, located at more than 300 m deep, consisted of a previous phase of compilation, analysis and synthesis of the existing information about the area, followed by a coordinated combination of different speciality studies: geology, stratigraphy, structural analysis, hydrogeology, hydrochemistry, geophysics and remote sensing. Geological studies included geological mapping, definition of stratigraphical units and facies and structural analysis. The aim of the hydrogeology study was to define aquifer formations, recharge area, aquifer points inventory and groundwater flow directions for the establishment of piezometric and water quality observation nets. Special techniques were applied, like thermal infrared aerial images and the evaluation of submarine groundwater discharge by means of natural radium isotopes. Hydrochemical techniques, including majority elements characterization and stable isotopes (¹⁸O, ²H and ³H) determination, allowed classifying hydrochemical facies and establishing a renewal pattern for water within the system. Geophysics was useful in determining the aquifer geometry, the features of the basement and the petrophysical characteristics of the geological formations. Preliminary results show an important tectonic complexity and the possibilities for groundwater uses in the area of study.     

Flexural isostatic response of the Alps to increased Quaternary erosion recorded by foreland basin remnants,SE France

We test the hypothesis that flexural isostatic compensation of the mass removed by enhanced Quaternary erosion is responsible for uplift of the Western European Alps and their forelands. We use two well‐preserved and well‐dated (1.8 Ma) abandonment surfaces of foreland basin remnants in SE France (the Chambaran and Valensole plateaux) as passive benchmarks for tilting of the foreland. Estimating their initial slope from morphometric scaling relationships, we determine bulk post‐depositional tilting of 0.5–0.8% for these surfaces. The calculated isostatic response of the Alpine lithosphere to erosional unloading, using the method recently proposed by Champagnac et al. [Geology 35 (2007) 195–198] , yields a predicted tilting of 0.3–0.4% in the considered areas, explaining approximately half of the determined post‐depositional tilting. Such long‐term deformation being insensitive to cyclic loading/unloading because of glaciations, we suspect the other half to be related to as yet undetermined long‐wavelength and long‐lived tectonic process(es).     

On the basic principles of geophysical radioactive measurements

Summary The process of radon diffusion in space and time through rocks is treated in detail from the point of view of radiometric prospecting for oil. Taking into account the depth and diffusivity conditions encountered in nature it is stated that vertical radon diffusion reaching from the oil deposits up to the surface cannot be expected. Time intervals as required for setting in of a stationary concentration state and decay conditions of the radon lead to the conclusion that the penetration of emanometric measurements seldom exceeds the depth of 5–10 m.     

Zur Ermittlung der bestanschmiegenden Ebene zu einer Anzahl von,aus Erdbeben bestimmten Richtungen

Zusammenfassung Es werden mit Hilfe der sphärischen Trigonometrie und der Ausgleichsrechnung Methoden angegeben, die die Berechnung der bestanschmiegenden Ebene an eine Anzahl von Geraden mit einem gemeinsamen Punkt ermöglichen.

---

Summary With the aid of the spherical trigonometry and the adjustment, methods are given to determin the best-fitting plane for a given set of directions through a common point.