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.     

Hydrodynamic control of phytoplankton in low salinity waters of the James River estuary,Virginia, U.S.A.

Autotrophic biomass and productivity as well as nutrient distributions and phytoplankton cell populations in the James River estuary, Virginia, were quantified both spatially and temporally over a 17-month period. Emphasis was placed on the very low salinity region of the estuary in order to gain information on the fate of freshwater phytoplankters. Differing amounts of freshwater plant biomass are advected into the estuary as living material, DOC or POC and the demonstrated variability of this input must play an important role in marine biogeochemical cycling.Late summer and fall maxima in both chlorophyll a and the photosynthetic production of particulate organic carbon in very low salinity regions were inversely correlated with river discharge.During periods of low river discharge greater than 50% of the chlorophyll a biomass measured at 0‰ disappeared within a narrow range of salinity (0–2‰). Cell enumeration data suggest that species introduced from the freshwater end-member tend to comprise the bulk of the biomass removed. Confounding factors, which may contribute to the regulation of both the abundance and species of phytoplankters mid-river, include the flocculation of colloidal material with phytoplankton cells, the presence of the turbidity maximum and the growth of endemic phytoplankton populations.An inverse relationship exists between the phytoplankton abundance in very low salinity waters and the abundance of biomass measured in the lower portion of the river (estuary). Thus, autotrophic production in the fresh and very low salinity areas may indirectly regulate the onset on the spring bloom in the estuary by controlling the amount of nutrients available.     

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.     

Regolith Surface Reflectance: A New Attempt to Model

The reflectance coefficient of the regolith layer of celestial bodies has been studied in relation to the physical properties of regolith particles (size, refractive index, and packing density) on the basis of an accurate numerical radiative-transfer algorithm for a semi-infinite flat layer. Using the geometric-optics approximation, we have found that a shape mixture of randomly oriented spheroids can successfully model the single-scattering phase function of independent soil grains. In order to take into account the effect of packing density in a regolith layer, the concept of the so-called static structure factor was used. The main effect of increasing packing density is to suppress the forward-scattering peak of the phase function and to increase the albedo of the reflecting surface. We also investigated the influence of fine dust on the reflected light. An addition of small particles not only increases the surface albedo, but also changes the brightness profile and enhances the backscattering. Although the problem of unique solution, which is inherent in the retrieval of the properties of a medium from the measurements of the intensity of light scattered by this media, cannot be removed in the proposed model, the procedure used here, in contrast to widely used approximations, allows us to fit observational data with a set of real characteristics of the regolith. Semiempirical approaches are able to fit the measurements well with a small number of free parameters, but they do not explicitly contain crucial physical characteristics of the regolith such as grain sizes or the refractive index. We compared the numerical solution of the radiative-transfer equation with the Hapke approximation, which is most often used by investigators. The errors introduced by the Hapke model are small only for near-isotropic scattering by isolated particles. However, independent regolith grains are known to scatter light mainly in the forward direction.     

A unified model of dwarf nova outbursts based on the disk instability

A unified model for outbursts of dwarf novae is proposed based on the disk instability model in cataclysmic variable stars. In this model, two different intrinsic instabilities (i.e., the thermal instability and the tidal instability) within accretion disks are considered in non-magnetic cataclysmic variable stars. It is suggested that all of three sub-classes of dwarf novae (i.e., U Gem-type, Z Cam-type and SU UMa-type dwarf novae) may be explained in terms of two model parameters of the orbital period of the binary and of the mass transfer rate within the framework of the disk instability model.     

The Effect of Atmospheric Drag on Satellite Orbits During the Bastille Day Event

Geomagnetic storms driven by solar eruptions are known to have significant effects on the total density of the upper atmosphere in the altitude range 250–1000 km. This in turn causes a measurable effect on the orbits of resident space objects in this altitude range. We analyzed a sample of these orbits, both from sensor data and from orbital element sets, during the period surrounding the 14 July 2000 solar activity. We present information concerning the effects of this event on the orbits of resident space objects and how well accepted atmospheric models were able to represent it. As part of this analysis, we describe a technique for extracting atmospheric density information from orbital element sets. On daily time scales, the effect of geomagnetic activity appears to be more important than that of prompt radiation. However, the limitations in time and amplitude quantization of the accepted solar indices are evident. A limited comparison is also made with previous solar storm events.     

Introduction

A historical sketch of the nonlinear theory of variable stars is outlined briefly. The main break-through came from the hydrodynamic study of stellar pulsation. From a theoretical point of view, coupling needs to be discussed more carefully. The impact of new opacities on the astrophysical problems is also discussed.