Transient response of a circular cavity in a poroelastic medium

This paper considers the transient response of a pressurized long cylindrical cavity in an infinite poroelastic medium. To obtain transient solutions, Biot's equations for poroelastodynamics are specialized for this problem. A set of exact general solutions for radial displacement, stresses, pore pressure and discharge are derived in the Laplace transform space by using analytical techniques. Solutions are presented for three different types of prescribed transient radial pressures acting on the surface of a permeable as well as an impermeable cavity surface. Time domain solutions are obtained by inverting Laplace domain solutions using a reliable numerical scheme. A detailed parametric study is presented to illustrate the influence of poroelastic material parameters and hydraulic boundary conditions on the response of the medium. Comparisons are also presented with the corresponding ideal elastic solutions to portray the poroelastic effects. It is noted that the maximum radial displacement and hoop stress at the cavity surface are substantially higher than the classical static solutions and differ considerably from the transient elastic solutions. Time histories and radial variations of displacement, hoop stress, pore pressure and fluid discharge corresponding to a cavity in two representative poroelastic materials are also presented.     

Shakedown in layered pavements under moving surface loads

Elastic–plastic deformations in pavements consisting of layers of different frictional materials are investigated. The upper bound, kinematic shakedown theorem is used to obtain estimates of the critical shakedown loads. Fully general three-dimensional deformations are considered. The influence of the loading distribution, interactions between loads, and the effect of varying the thickness, stiffness and strength of the layers are explored. Consequences of this investigation for particular existing designs of flexible pavements are investigated. In particular, it is found that the strength of the subgrade has no effect on the magnitude of the critical shakedown load.     

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 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.     

Optimization of grid-drilling using computer simulation

A computer simulation method has been developed to find efficient drilling grids for mineral deposits. A well-known ore deposit is used as a model to develop an efficient pattern for undiscovered ore bodies in the same area or in other prospects where similar geometry is suspected. The model for this study is the Austinville, Virginia deposit, a Mississippi Valley-type deposit composed of 17 ore bodies totaling 34 million short tons (30 million metric tons). The method employs a computer program that simulates drilling the model deposit with different patterns, including various levels of follow-up drilling. Follow-up holes are drilled in fences at one half the original spacing around holes in the grid that show ore-grade mineralization. Each pattern is drilled 100 times from random starting locations to provide a range of outcomes of drilling, including the best, worst, and most likely. For this study, patterns of 100 drill holes were composed of 10 fences spaced 1000–5000 feet (305–1524 m) apart, each with 10 holes spaced 200–1000 feet (61–305 m) apart. In all, 25 grids were used with zero to three levels of follow-up drilling. The 600/2000 grid, with drill holes spaced 600 feet (183 m) apart in fences spaced 2000 feet (610 m) apart, was compared with the 200/5000 grid because they represented contrasting outcomes. The 600/2000 grid penetrated many ore bodies consistently but with few multiple hits to individual ore bodies; whereas the 200/5000 grid inconsistently penetrated few ore bodies with many multiple hits. The 600/2000 grid was more efficient than the 200/5000 grid at hitting large ore bodies of 1,000,000 short tons or greater (900,000 metric tons or greater) and was made more effective by adding one cycle of follow-up drilling. The 600/2000 grid had a 97% chance of hitting one or more large ore bodies with at least one drill hole per ore body, and the 200/5000 grid had a 64% chance. Once hit, there was an 82% chance that the largest ore body would be penetrated by three or more holes when using the 600/2000 grid and an 88% chance using the 200/5000 grid.     

Recharge and salinization in the Madrid Basin

The purpose of this work was to reinvestigate the existing hydrogeological conceptual model of the basin of Madrid, Spain. A cumulative chemical isotopic diagram which enabled the distinction between different groups of water as well as calculation of the mode of their blending was used for this investigation. It was found that the groups of discharge were lighter in their isotopic composition than that of recharge. The previous explanation of this fact, backed by carbon-14 dating, was the long residence time due to flow lines going down to depths of more than 1000 m. This flow model assumes homogenous conditions to these depths. This assumption can not be supported by evidence from deep wells. Thus a modified model is suggested which maintains homogenous conditions only to about 300 m and a deep confined aquifer below containing paleowater. The higher degree of depletion of this water has been explained by a colder climate on top of an altitude effect. Another interesting observation was the correlation between the isotopic composition of the rains, the month of the rain event and the composition of the recharge group groundwater. It could be seen that the winter rains resemble the groundwater composition, which shows that practically all the spring and summer rains were evapotranspirated.     

THE EUROPEAN MACROSEISMIC SCALE (MSK-92)

A new updated version of the MSK macroseismic intensity scale has been prepared by a Working Group of the European Seismological Commission and has been published in April 1993 (European Macroseismic Scale 1992: updated MSK scale, 1993, ed. by G. Grünthal, Cahiers du Centre Européen de Geodynamique et de Séismologie, no. 7).