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.     

A shipborne- type wave-height meter for oceangoing vessels, using microwave Doppler radar

A wave-height meter using a simple microwave Doppler radar,simeq10mW in power and 10.525 GHz in frequency, is proposed so that we can measure oceanic waves effectively while the ship is steaming. It was first applied to the measurement of the variation of water level generated in a wave tank, which suggested that it is adequately applicable to the measurement of oceanic waves. A field test was carried out off the cape of Nojimazaki by installing the Doppler radar 5 m above the sea level at the bow of the ship. The result agreed reasonably well with that measured simultaneously by the ultrasonic wave-height meter installed at the same position. Another test is running successfully on a larger ship with the wave-height meter installed at 9 m above the sea level. The significant wave height measured by the present meter is being compared with that observed visually by the navigation officers.     

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.     

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.     

Cokriging nonstationary data

Universal cokriging is used to obtain predictions when dealing with multivariate random functions. An important type of nonstationarity is defined in terms of multivariate random functions with increments which are stationary of orderk. The covariance between increments of different variables is modeled by means of the pseudo-cross-covariance function. Criteria are formulated to which the parameters of pseudo-cross-covariance functions must comply so as to ensure positive-definiteness. Cokriging equations and the induced cokriging equations are given. The study is illustrated by an example from soil science.     

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.     

Seismotectonic significance of the 29 January 1989 Etne earthquake, SW Norway

The integrated analysis of geological, seismological and field observations with lineament data derived from satellite images allows the identification of a possible seismogenic fault zone for an earthquake which occurred near Etne in southwestern Norway, on 29 February 1989. The hypocentre of the earthquake was located at the mid-crust at a depth of 13.8±0.9 km which is typical of small intraplate earthquakes. The Etne earthquake occurred as a result of normal faulting with a dextral strike-slip component on a NW–SE trending fault. Available geological and lineament data indicate correlation of the inferred seismogenic fault with the NW–SE trending Etne fault zone. An aeromagnetic anomaly related to the Etne fault zone forms a regional feature intersecting both Precambrian basement and allochthonous Caledonian rocks. Based on these associations the occurrence of the Etne event is ascribed to the reactivation of a zone of weakness along the Etne fault zone. Slope-instabilities developed in the superficial deposits during the Etne event demonstrate the existence of potentially hazardous secondary-effects of such earthquakes even in low seismicity areas such as southwestern Norway.     

Thaw modification of frost-fissure wedges,Richards Island,Pleistocene Mackenzie Delta,Western Arctic Canada

Thaw modification is the general process whereby frost-fissure wedges are modified during thaw, and by which frost-fissure pseudomorphs may develop. Specific processes of thaw modification are inferred from ice-wedge pseudomorphs, composite-wedge pseudomorphs and deformed sand wedges in the Pleistocene Mackenzie Delta: i.e. thermal erosion, collapse, subsidence, refreezing, loading, buoyancy, spreading, folding and shearing. Thaw modification is believed to result in selective preservation of pseudomorphs and wedges. Sand wedges are more likely to be preserved than are ice-wedge pseudomorphs or compositewedge pseudomorphs, because only those sand wedges that penetrate massive ice or icy sediments are prone to thaw modification. Furthermore, whereas ice wedges preferentially develop in ice-rich, fine-grained sediments (thaw-sensitive), their pseudomorphs appear to be selectively preserved in ice-poor, coarse-grained sediments (thaw-stable).