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.     

Readout Techniques for Drift and Low Frequency Noise Rejection in Infrared Arrays

Three different methods are presented to subtract thermal drifts and low-frequency noise from the signal of infrared array. The first is dead pixels with open Indium bumps, the second is reference output as implemented on the Hawaii2 multiplexer, and the third is dark pixels to emulate reference cells having a capacity connected to the gate of the unit cell field-effect transistor (FET). The third method is the most effective and yields a reduction in readout noise from15.4–9.4 erms. A novel method will be described to extend this readout technique to the Aladdin 1 K × 1 K InSb array. This revised version was published online in July 2006 with corrections to the Cover Date.     

The height,spectrum and mechanism of type-B red aurora and its bearing on the excitation of O(1S) in aurora

The height of the lower red border of type-B aurora has been determined by triangulation using TV cameras at two ground stations. A mean height of 91.4 ± 1.1 km was determined from a set of 12 measurements made under ideal conditions. A TV spectrograph was used simultaneously to seek possible spectral changes between 6400 and 6900 Å which would be indicative of changes in the vibrational distribution in the N₂ First Positive bands. No significant difference was found in this distribution between the spectra from 93 and 122 km. The height distribution of contributions to the OI 5577 Å emission relative to the N⁺₂ First Negative emission was modelled from 80 to 160 km. Contributions from electron impact on atomic O, O⁺₂ dissociative recombination and N₂(A)O energy transfer were included. Account was taken of recent laboratory data on O(¹S) quenching. It was concluded that these processes could explain the excitation of O(¹S) in normal aurora and the height distribution of OI 5577 Å in type-B red aurora. It was confirmed that the lifetime ofO(¹S) in type-B red auroral rapid time variations is about 0.5 s and it was found from the model that the observed time variation can be reproduced by the mechanisms considered, provided the concentration of NO in the auroral atmosphere is about 1 × 10⁹ at 95 km. Before reasonable certainty can be attained in the correctness of the interpretation it will however be necessary to have reliable simultaneous observations of neutral atmospheric composition particularly for O and NO as well as unchallengeable measurements of the yields of O(¹S) for the processes considered and for several other processes which have been suggested recently.     

The Expansion of the Hamiltonian of the Planetary Problem into the Poisson Series in All Keplerian Elements (Theory)

The study of the evolution of planetary systems, primarily of the Solar System, is one of the basic problems of celestial mechanics. The stability of motion of giant planets on cosmogonic time scales was established by numerical and analytical methods, but the question about the evolution of orbits of terrestrial planets and arbitrary solar-type planetary systems remained open. This work initiates a series of papers allowing one to advance in solving the problem of the evolution of the solar-type planetary systems on cosmogonic time scales by using powerful analytical tools. In the first paper of this series, we choose the optimum reference system and obtain the Poisson series expansion of the Hamiltonian of the problem in all Keplerian elements. We propose to use the integral representation of the corresponding coefficients or the Poisson processor means instead of conventionally addressing any possible special functions. This approach extremely simplifies the algorithm. The next paper of this series deals with the calculation of the expansion coefficients.     

Pulsar reflections within the Crab Nebula

Between 1997 August and October, the radio pulses from the Crab pulsar were followed by discrete moving echoes, which appear to be reflections from part of an ionized shell in the outer part of the Crab Nebula, crossing the line of sight to pulsar. Similar events have now been recognized in recordings from the past 30 yr, and it seems that the Nebula must contain a large number of ionized shell-like surfaces on a much finer scale than recognized hitherto.     

The interannual leading modes of the extratropical variability in the Southern Hemisphere simulated by the ECHAM-4 atmospheric model

 An ensemble of twenty-three 14-year experiments conducted with the ECHAM-4 GCM has been examined to test the model's capability to simulate the principal modes of interannual variability. The integrations were performed under specified monthly SST between 1979–1993. The analysis was focused on the Southern Hemisphere (SH) extratropics. Empirical orthogonal functions analysis (EOF) using seasonal anomaly fields has been performed to isolate the principal modes that dominate the southern extratropical variability at the interannual time scale. Leading patterns of 500 hPa geopotential height (z500) have been compared with those estimated from the ECMWF re-analysis dataset. The model is able to adequately reproduce the spatial pattern of the annular mode, but it represents the temporal variations of the oscillation less satisfactorily. The model simulation of the Pacific South American (PSA) pattern is better, both in the shape of the pattern and in the temporal evolution. To verify if the capability of the model to adequately simulate the temporal oscillation of the propagating patterns is related to the increased influence of the tropical external forcing, covarying SST-atmospheric modes have been identified by singular value decomposition (SVD). In winter (July-August-September, JAS) the tropical SST variability is highly correlated with the ENSO mode. In summer (January-February-March, JFM) the strength of the teleconnections is related to strong westerly anomalies, disrupted by a meridional out of phase relation near to South America. The large size of the ensemble was exploited by comparing the time-varying model spread and degrees of freedom of the simulated extratropical circulation. Results show that when the extratropical circulation has a few degrees of freedom, the reproducibility is relatively low and the ensemble is governed by a fairly robust zonally symmetric structure of dispersion. Received: 9 May 2000 / Accepted: 30 January 2001     

Long-term seismicity of Istanbul and of the Marmara Sea region

We examine the long-term seismicity of the Marmara Sea region in Turkey over a period of twenty centuries to attempt to answer two questions: (1) how well recorded in history is the earthquake activity of this region? and (2) does seismicity over this long period of time differ from that over the present century? We study this densely-populated and fast-developing region because it is one of the most tectonically active regions on the continents that appears today to be relatively quiescent.     

A viable mechanism for the production of energy in active galactic nuclei

It is suggested that a collapsing supermassive object, which acts as an ultra-high energy particle accelerator, is the precursor of an active galactic nucleus and that the gravitational energy released during the collapse of the object is locked in the quark-gluon plasma permeated by leptons into which the entire matter in the core of the object is converted as a result of the collapse. It is also pointed out that the collapse of the object to a space-time singularity is inhibited by Pauli's exclusion principle as well as by Heisenberg's uncertainty principle and that the object explodes, before it could collapse to a singularity, thereby releasing the enormous amount of energy locked in the quark-gluon plasma.