Upscaling of elastic properties for large scale geomechanical simulations

Large scale geomechanical simulations are being increasingly used to model the compaction of stress dependent reservoirs, predict the long term integrity of under‐ground radioactive waste disposals, and analyse the viability of hot‐dry rock geothermal sites. These large scale simulations require the definition of homogenous mechanical properties for each geomechanical cell whereas the rock properties are expected to vary at a smaller scale. Therefore, this paper proposes a new methodology that makes possible to define the equivalent mechanical properties of the geomechanical cells using the fine scale information given in the geological model. This methodology is implemented on a synthetic reservoir case and two upscaling procedures providing the effective elastic properties of the Hooke's law are tested. The first upscaling procedure is an analytical method for perfectly stratified rock mass, whereas the second procedure computes lower and upper bounds of the equivalent properties with no assumption on the small scale heterogeneity distribution. Both procedures are applied to one geomechanical cell extracted from the reservoir structure. The results show that the analytical and numerical upscaling procedures provide accurate estimations of the effective parameters. Furthermore, a large scale simulation using the homogenized properties of each geomechanical cell calculated with the analytical method demonstrates that the overall behaviour of the reservoir structure is well reproduced for two different loading cases. Copyright © 2004 John Wiley & Sons, Ltd.     

Geometry and kinematics of recent deformation in the Mondy–Tunka area (south-westernmost Baikal rift zone, Mongolia–Siberia)

We used satellite imagery and field data to investigate the south‐westernmost Baikal rift zone. We focus our study in the Mondy and Ikhe Ukhgun valleys, site of an Mw = 6.9 seismic event in 1950. Surface deformations are observed along the E–W‐trending Mondy strike‐slip fault and along the Ikhe Ukhgun thrust. The Mondy fault system is 80 km long and is composed of four segments 10–15 km long. These segments are characterized by subvertical planes with left‐lateral movements. The Ikhe Ukhgun thrust is 20 km long, dips 40° to the south and shows reverse movement with a left‐lateral component. These observations are consistent with the present‐day regional NNE–SSW compression and with the focal mechanism of the 1950 Mondy earthquake that was recently re‐evaluated. These features, like those observed in the Tunka basin, demonstrate a recent change of regional strain regime from transtension to transpression that we place before the Late Pleistocene.     

Preliminary organic petrography of Palaeozoic rocks from the Grand Banks,Newfoundland

Samples of lower Palaeozoic bedrock from the Grand Banks of Newfoundland were examined using reflected light microscopy and Rock-Eval pyrolysis. These samples contained organic material which included bitumen (structureless organic material), acritarchs, chitinozoa, scolecodonts and graptolites. The reflectance of the organic material increased from acritarchs towards graptolites with bitumen showing a wide range of reflectance. Multiple phases of oil migration through one sample were inferred by examination of the bitumens: the lowest reflecting being more recent than the highest reflecting bitumen. Combined reflectance and Rock-Eval analysis indicate that these samples are mature to overmature. The two oldest samples (Arenig-Llanvirn, Ordovician) are assessed as mature and as having potential for generation of liquid hydrocarbons, but must be considered as mainly gas prone.     

The stable isotope signature of kilometre-scale fracturedominated metamorphic fluid pathways, Mary Kathleen, Australia

Abstract Large calcite veins and pods in the Proterozoic Corella Formation of the Mount Isa Inlier provide evidence for kilometre-scale fluid transport during amphibolite facies metamorphism. These 10- to 100-m-scale podiform veins and their surrounding alteration zones have similar oxygen and carbon isotopic ratios throughout the 200 × 10-km Mary Kathleen Fold Belt, despite the isotopic heterogeneity of the surrounding wallrocks. The fluids that formed the pods and veins were not in isotopic equilibrium with the immediately adjacent rocks. The pods have δ¹³C_calcite values of –2 to –7% and δ¹⁸O_calcite values of 10.5 to 12.5%. Away from the pods, metadolerite wallrocks have δ¹⁸O_whole-rock values of 3.5 to 7%. and unaltered banded calc-silicate and marble wallrocks have δ¹³C_calcite of –1.6 to –0.6%, and δ¹⁸O_calcite of 18 to 21%. In the alteration zones adjacent to the pods, the δ¹⁸O values of both metadolerite and calc-silicate rocks approach those of the pods. Large calcite pods hosted entirely in calc-silicates show little difference in isotopic composition from pods hosted entirely in metadolerite. Thus, 100- to 500-m-scale isotopic exchange with the surrounding metadolerites and calc-silicates does not explain the observation that the δ¹⁸O values of the pods are intermediate between these two rock types. Pods hosted in felsic metavolcanics and metasiltstones are also isotopically indistinguishable from those hosted in the dominant metadolerites and calc-silicates. These data suggest the veins are the product of infiltration of isotopically homogeneous fluids that were not derived from within the Corella Formation at the presently exposed crustal level, although some of the spread in the data may be due to a relatively small contribution from devolatilization reactions in the calc-silicates, or thermal fluctuations attending deformation and metamorphism. The overall L-shaped trend of the data on plots of δ¹³C vs. δ¹⁸O is most consistent with mixing of large volumes of externally derived fluids with small volumes of locally derived fluid produced by devolatilization of calc-silicate rocks. Localization of the vein systems in dilatant sites around metadolerite/calc-silicate boundaries indicates a strong structural control on fluid flow, and the stable isotope data suggest fluid migration must have occurred at scales greater than at least 1 km. The ultimate source for the external fluid is uncertain, but is probably fluid released from crystallizing melts derived from the lower crust or upper mantle. Intrusion of magmas below the exposed crustal level would also explain the high geothermal gradient calculated for the regional metamorphism.     

The middle pleistocene deposits of Waverley Wood Pit,Warwickshire, England

A complex of channels underlying the Baginton-Lillington Gravel (Baginton Formation) at Waverley Wood Quarry, Warwickshire is described. Fossil pollen and plant macrofossils, Coleoptera, Ostracoda, Mollusca and Mammalia are described from the channel-fill deposits. Consideration of all the evidence allows the identification of four separate stages of channel fill which largely occurred under a cool temperate climate. At the top of Channel 2 evidence for a cold, continental climatic episode can be recognised, suggesting that the whole complex was deposited under a fluctuating climate at the end of a temperate stage. At two levels in the channels human artefacts were recovered confirming the presence of Palaeolithic people in Warwickshire during the deposition of the sediments. Amino-acid geochronology suggests an age within the ‘Cromerian Complex’ Stage for the channels. The small vertebrate and molluscan faunas indicate that the deposits are no older than the latter part of the ‘Cromerian Complex’ Stage of East Anglia. The regional stratigraphic significance of the Waverley Wood succession is outlined.     

A high‐resolution pollen and geochemical analysis of late Holocene human impact and vegetation history in southern Cumbria,England

The historic era, which in Cumbria begins with the Roman invasion of AD 71, is a frequently neglected period in palaeoecological research, but its study can bring benefits in improving knowledge of landscape history and in understanding the significance and limitations of palaeoecological records. Pollen and geochemical data are presented for late Holocene records from Deer Dyke and Hulleter Mosses in southern Cumbria. The records show initially low levels of anthropogenic impact, followed by a phase of forest clearance and mixed agriculture from the 7th to 11th centuries AD. The timing of these clearances suggests that they were initially Anglo‐Saxon in origin, rather than Norse. Further clearances in the 16th century AD are interpreted as a response to monastic dissolution and late Tudor population pressures; the landscapes reached their contemporary form following extensive clearances in the 17th century AD. Silicon and titanium concentrations at Deer Dyke Moss were used to reconstruct past levels of atmospheric dust loading, which is broadly related to soil erosion. Geochemical influx was found to peak during periods of landscape transition rather than from established land use. This relationship with pollen data is thought to reflect the predominantly low levels of anthropogenic impact in the region, which changes as substantial woodland clearances during the 16th century AD and continuous land use pressure since then have greatly increased the supply of airborne dust. Copyright © 2008 John Wiley & Sons, Ltd.     

VLA and MERLIN observations of RZ Cassiopeiae

We present radio interferometric observations of the Algol-type binary system RZ Cassiopeiae made with the VLA and MERLIN arrays at 6 cm over an incomplete orbital cycle of the system (1.195 d). We detected RZ Cas with both instruments. The images were unresolved in both cases, with angular extents comparable to the synthesized beams. The peak flux density in the VLA image was 1.14 mJy beam⁻¹ and in the MERLIN image it was 0.93 mJy beam⁻¹. The derived brightness temperatures are  4.02 × 10⁸  and  4.35 × 10⁸ K  and the effective electron energies are 0.347 and 0.346 MeV for the MERLIN and VLA data respectively. The radio light curve shows an interesting modulation centred close to the primary eclipse which seems to correlate with ASCA SIS observations of the system. The results can be interpreted as an emitting region on the outer hemisphere of the cool component aligned along the centroid axis of the binary system.     

Contribution of storms to groundwater recharge in the semi‐arid region of Karnataka,India

This paper describes the application of environmental isotopes and injected tracer techniques in estimating the contribution of storms as well as annual precipitation to groundwater recharge and its circulation, in the semi‐arid region of Bagepalli, Kolar district, Karnataka. Environmental isotopes ²H, ¹⁸O and ³H were used to study the effect of storms on the hydrological system, and an isotope balance was used to compute the contribution of a storm component to the groundwater. Some of the groundwater samples collected during the post‐storm periods were highly depleted in stable isotope content with higher deuterium excess relative to groundwater from the pre‐storm periods. Significant variation in deuterium excess in groundwater from the same area, collected in two different periods, indicates the different origin of air masses. The estimated recharge component of a storm event of 600 mm to the groundwater was found to be in the range of 117–165 mm. There was no significant variation in environmental tritium content of post‐storm and pre‐storm groundwater, indicating the fast circulation of groundwater in the system. After completion of the environmental isotope work, an injected radiotracer ³H technique was applied to estimate the direct recharge of total precipitation to the groundwater. The estimated recharge to the groundwater is 33 mm of the 550 mm annual precipitation during 1992. Copyright © 2003 John Wiley & Sons, Ltd.     

Magnetic Field Geometry in `Red' and `Blue' BL Lacs

We compare the systematics of the magnetic field geometry in the `red' low-energy peaked BL Lacs (LBLs) and `blue' high-energy peaked BL Lacs (HBLs) using VLBI polarimetric images. The LBLs are primarily `radio-selected' BL Lacs and the HBLs are primarily `X-ray selected'. In contrast to the LBLs, which show predominantly transverse jet magnetic fields, the HBLs show predominantly longitudinal fields. Thus, while the SED peaks of core-dominated quasars, LBLs and HBLs form a sequence of increasing frequency, the magnetic field geometry does not follow an analogous sequence. We briefly investigate possible connections between the observed parsec-scale magnetic field structures and circular polarization measurements in the literature on various spatial scales.     

Himalayan magmatism and porphyry copper–molybdenum mineralization in the Yulong ore belt, East Tibet

Summary ?The NW–SE-trending Yulong porphyry Cu–Mo ore belt, situated in the Sanjiang0 area of eastern Tibet, is approximately 400 km long and 35 to 70 km wide. Complex tectonic and magmatic processes during the Himalayan epoch have given rise to favorable conditions for porphyry-type Cu–Mo mineralization. Porphyry masses of the Himalayan epoch in the Yulong ore belt are distributed in groups along regional NW–SE striking tectonic lineaments. They were emplaced mainly into Triassic and Lower Permian sedimentary-volcanic rocks. K–Ar und U–Pb isotopic datings give an intrusion age range of 57–26 Ma. The porphyries are mainly of biotite monzogranitic and biotite syenogranitic compositions. Geological and geochemical data indicate that the various porphyritic intrusions in the belt had a common or similar magma source, are metaluminous to peraluminous, Nb–Y–Ba-depleted, I-type granitoids, and belong to the high-K calc-alkaline series. Within the Yulong subvolcanic belt a number of porphyry stocks bear typical porphyry type Cu–Mo alteration and mineralization. The most prominent porphyry Co–Mo deposits include Yulong, Malasongduo, Duoxiasongduo, Mangzong and Zhanaga, of which Yulong is one of the largest porphyry Cu (Mo) deposits in China with approximately 8 × 10⁶ tons of contained Cu metal. Hydrothermal alteration at Yulong developed around a biotite–monzogranitic porphyry stock that was emplaced within Upper Triassic limestone, siltstone and mudstone. The earliest alteration was due to the effects of contact metamorphism of the country rocks and alkali metasomatism (potassic alteration) within and around the porphyry body. The alteration of this stage was accompanied by a small amount of disseminated and veinlet Cu–Mo sulfide mineralization. Later alteration–mineralization zones form more or less concentric shells around the potassic zone, around which are distributed a phyllic or quartz–sericite–pyrite zone, a silicification and argillic zone, and a propylitic zone. Fluid inclusion data indicate that three types of fluids were involved in the alteration–mineralization processes: (1) early high temperature (660–420 °C) and high salinity (30–51 wt% NaCl equiv) fluids responsible for the potassic alteration and the earliest disseminated and/or veinlet Cu–Mo sulfide mineralization; (2) intermediate unmixed fluids corresponding to phyllic alteration and most Cu–Mo sulfide mineralization, with salinities of 30–50 wt% NaCl equiv and homogenization temperatures of 460–280 °C; and (3) late low to moderate temperature (300–160 °C) and low salinity (6–13 wt% NaCl equiv) fluids responsible for argillic and propylitic alteration. Hydrogen and oxygen isotopic studies show that the early hydrothermal fluids are of magmatic origin and were succeeded by increasing amounts of meteoric-derived convective waters. Sulfur isotopes also indicate a magmatic source for the sulfur in the early sulfide mineralization, with the increasing addition of sedimentary sulfur outward from the porphyry stock. Received August 29, 2001; revised version accepted May 1, 2002 Published online: November 29, 2002