A new computational strategy for solving two‐phase flow in strongly heterogeneous poroelastic media of evolving scales

We develop a new computational methodology for solving two‐phase flow in highly heterogeneous porous media incorporating geomechanical coupling subject to uncertainty in the poromechanical parameters. Within the framework of a staggered‐in‐time coupling algorithm, the numerical method proposed herein relies on a Petrov–Galerkin postprocessing approach projected on the Raviart–Thomas space to compute the Darcy velocity of the mixture in conjunction with a locally conservative higher order finite volume discretization of the nonlinear transport equation for the saturation and an operator splitting procedure based on the difference in the time‐scales of transport and geomechanics to compute the effects of transient porosity upon saturation. Notable features of the numerical modeling proposed herein are the local conservation properties inherited by the discrete fluxes that are crucial to correctly capture the fingering patterns arising from the interaction between heterogeneity and nonlinear viscous coupling. Water flooding in a poroelastic formation subject to an overburden is simulated with the geology characterized by multiscale self‐similar permeability and Young modulus random fields with power‐law covariance structure. Statistical moments of the poromechanical unknowns are computed within the framework of a high‐resolution Monte Carlo method. Numerical results illustrate the necessity of adopting locally conservative schemes to obtain reliable predictions of secondary recovery and finger growth in strongly heterogeneous deformable reservoirs. Copyright © 2011 John Wiley & Sons, Ltd.     

Magnetic resonance sounding applied to aquifer characterization

Magnetic resonance sounding (MRS) is distinguished from other geophysical tools used for ground water investigation by the fact that it measures a magnetic resonance signal generated directly from subsurface water molecules. An alternating current pulse energizes a wire loop on the ground surface and the MRS signal is generated; subsurface water is indicated, with a high degree of reliability, by nonzero amplitude readings. Measurements with varied pulse magnitudes then reveal the depth and thickness of water saturated layers. The hydraulic conductivity of aquifers can also be estimated using boreholes for calibration. MRS can be used for both predicting the yield of water supply wells and for interpolation between boreholes, thereby reducing the number of holes required for hydrogeological modeling. An example of the practical application of MRS combined with two-dimensional electrical imaging, in the Kerbernez and Kerien catchments area of France, demonstrates the efficiency of the technique.     

Estimates of coda Q attenuation in the João Câmara area (Northeastern Brazil)

S coda wave of seventy-four local earthquakes recorded in a network of ten seismic stations were used to calculate coda Q attenuation (Q_c) in the João Câmara area (northeastern Brazil). The estimates show Q_c as a strong function of frequency in the range from 6.0 to 20.0 Hz. We found out that Q_c in João Câmara has a functional form given by Q_c= Q₀ f, where Q₀= 151 ± 99 and = 0.98 ± 0.05. If the standard deviations are taken into account,we conclude that there are no relevant changes in both Q₀ and values from one station to another. The estimated Q₀ values at the different stations suggest that the Samambaia fault is a boundary between two different seismic attenuation zones. In one side of the fault (left), where stations were installed in Pre-Cambrian terrain and thick sedimentary layer, the seismic attenuation is stronger than in the other side (stations installed in thin sedimentary layer and limestone outcrop).The anomalous Q₀ values in the left side of the Samambaia fault can be explained due to the presence of a shallow conductive layer in the upper crust( 10 km), such as proposed by Padilha et al. (1992). According to our results, if there is a conductive layer in the area, it probably spreads over João Câmara city and surrounding regions.However, more detailed investigation either with seismic methods (seismic attenuation,3D tomography with P and/or S wave velocities) or with other geophysical methods is needed to interpret the observed differences in Q₀ values between the two sides of the Samambaia fault.     

Rainfall interception by an isolated evergreen oak tree in a Mediterranean savannah

Redistribution of ground‐level rainfall and interception loss by an isolated Quercus ilex tree were measured over 2 years in a Mediterranean oak savannah. Stemflow, meteorological variables and sap flow were also monitored. Rainfall at ground level was measured by a set of rain‐gauges located in a radial layout centred on the tree trunk and extending beyond the crown limits. Interception loss was computed as the difference between the volume of rainwater that would reach the ground in the absence of the tree and the volume of water that actually fell on the ground sampling area (stemflow included). This procedure provided correct interception loss estimates, irrespective of rainfall inclination. Results have shown a clear non‐random spatial distribution of ground‐level rainfall, with rainwater concentrations upwind beneath the crown and rain‐shadows downwind. Interception loss amounted to 22% of gross rainfall, per unit of crown‐projected area. Stand interception loss, per unit of ground area, was only 8% of gross rainfall and 28% of tree evapotranspiration. These values reflect the low crown cover fraction of the stand (0·39) and the specific features of the Mediterranean rainfall regime (predominantly with few large storms). Nevertheless, it still is an important component of the water balance of these Mediterranean ecosystems. Copyright © 2005 John Wiley & Sons, Ltd.     

A parallelizable method for two‐phase flows in naturally‐fractured reservoirs

A parallelizable, semi‐implicit numerical method is proposed for the study of naturally‐fractured reservoir systems. It has proved to be computationally efficient in producing accurate numerical solutions for the dual‐porosity model for immiscible, two‐phase flow in such reservoirs. The method combines hybridized mixed finite elements, a new version of the modified method of characteristics, a sophisticated operator‐splitting procedure for separating the pressure calculation in the fractures from that of the saturation, another operator splitting to handle the interaction of the matrix blocks and the fractures, and domain decomposition iterative procedures for both the pressure and the saturation. It permits moderately long time steps for the pressure and the saturation in the fractures and matrix blocks by using short, inexpensive microsteps to treat the transport portion of the saturation equation in the fractures. This paper is devoted to the formulation of the method and a discussion of numerical results for five‐spot and vertical cross‐section examples.     

The transition between blueschists and lawsonite-bearing eclogites based on observations from Corsican metabasalts

Syntectonic eclogites, associated with blueschist parageneses, have recrystallized in metabasalts from eastern Corsica under very low temperatures (420°C). The evolution of these eclogites is recorded by the order of development of metamorphic minerals, as demonstrated by helicitic inclusions of lawsonite and actinolite in Fe-rich garnets, and the occurrence of fibrous Na-pyroxenes in the pressure shadows and cracks of boudinaged garnets, within a foliated matrix composed of glaucophane, Ca---Fe garnets and lawsonite.

A Schreinemakers analysis has been completed on the lawsonite-ferroglaucophane-actinolite-almandine-grossular-ferro-omphacite-chlorite system in a (T,P,μ_H2O)-space. The resulting three-dimensional relationships are applied to the studied samples. They are consistent with a prograde increase in temperature at low μ_H2O, producing eclogites from blueschists, and with a subsequent retrogressive decrease in pressure.     

Petrology and U/Pb geochronology of the Telohat migmatites,Aleksod, Central Hoggar,Algeria

The Aleksod region is composed of metasedimentary rocks and large areas of biotite and hornblende-bearing migmatites. Anatexis associated with the main deformation stages, occurred under high pressure and temperature conditions estimated at 13±2 Kbar and 750±50° C. The bulk mineralogical composition of the Telohat migmatites shows that their protolith was granodioritic. Internal structures of zircons and U-Pb data suggest a polyphased evolution, with a 2131±12 Ma age for the protolith and a 609±17 Ma age for the Pan-African tectono-metamorphic evolution, thus precluding any Kibaran event in the Aleksod area. Leucosomes are richer in Sr and display lower Rb, Zr, Nb, Y, Th, U and REE contents than melanosomes wherein accessory phases are stored. Eu contents are also lower in the leucosomes but in lesser proportion than the other rare earth's, leading to a significant positive anomaly. Petrogenetic modelling accounting for accessory mineral phases clearly shows that the trace element contents of leucosomes and melanosomes follow a distribution law consistent neither with equilibrium nor fractional melting. Their trace element patterns are best explained by the model of disequilibrium melting, with mixing of a few residual phases. The present results and previous Sr isotopic data as well raise the question of disequilibrium melting in anatexis of crustal material CRPG Contribution no 782     

Enviromental genotoxicity evaluation: Bayesian approach for a mixture statistical model

 The data analyzed in this paper are part of the results described in Bueno et al. (2000). Three cytogenetics endpoints were analyzed in three populations of a species of wild rodent – Akodon montensis – living in an industrial, an agricultural, and a preservation area at the Itajaí Valley, State of Santa Catarina, Brazil. The polychromatic/normochromatic ratio, the mitotic index, and the frequency of micronucleated polychromatic erythrocites were used in an attempt to establish a genotoxic profile of each area. It was assumed that the three populations were in the same conditions with respect to the influence of confounding factors such as animal age, health, nutrition status, presence of pathogens, and intra- and inter-populational genetic variability. Therefore, any differences found in the endpoints analyzed could be attributed to the external agents present in each area. The statistical models used in this paper are mixtures of negative-binomials and Poisson variables. The Poisson variables are used as approximations of binomials for rare events. The mixing distributions are beta densities. The statistical analyzes are under the bayesian perspective, as opposed to the frequentist ones often considered in the literature, as for instance in Bueno et al. (2000).     

Le passage Ordovicien–Silurien et la partie inférieure du Silurien (Sud-Est du Massif armoricain, France)The Ordovician–Silurian transition and the lower part of the Silurian (southeastern Armorican Massif, France)

A new section in the Silurian graptolitic ‘phtanites’ (black cherts) of Les Fresnaies at Chalonnes-sur-Loire (SE Armorican Massif) shows for the first time that these rocks (1) succeed conformably to Uppermost Ordovician (Hirnantian) glaciomarine deposits and (2) contain successive graptolite assemblages that characterise the base of the Silurian, the whole Rhuddanian and Aeronian stages and the lower part of the Telychian. To cite this article: J.M. Piçarra et al., C. R. Geoscience 334 (2002) 1177–1183.