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.     

济南市地下水资源现状分析及开发利用建议

本文采用较新的统计资料,对济南市地下水资源的开发利用现状及开采潜力,以行政区划为单元进行了分析论证.针对地下水资源开发利用中存在的环境水文地质问题,提出了开发利用规划和建议.     

Late Quaternary landscapes in Central Australia: sedimentary history and palaeoecology of Puritjarra rock shelter

Puritjarra rock shelter provides a long record of late Quaternary vegetation in the Australian arid zone. Analysis of the sedimentary history of this rock shelter is combined with reanalysis of charcoal and phytolith records to provide a first‐order picture of changing landscapes in western Central Australia. These show a landscape responding to increasing aridity from 45 ka with deflation of clay‐rich red palaeosols (<45 ka) and sharp declines in grassland and other vegetation at 40–36 ka, and at the beginning of the Last Glacial Maximum (LGM) (24 ka). Vegetation in the catchment of the rock shelter recovered after 15 ka with expansion of both acacia woodland and spinifex grasslands, registering stronger summer rainfall in the interior of the continent. By 8.3 ka re‐vegetation of local palaeosols and dunes had choked off sediment supply to the rock shelter and the character of the sediments changed abruptly. Poaceae values peaked at 5.8 ka, suggesting the early–mid Holocene climatic optimum in Central Australia is bracketed between 8.3 and 5.8 ka. Local vegetation was disrupted in the late Holocene with a sharp decline in Poaceae at 3.8 ka, coinciding with an abrupt intensification of ENSO. Local grasslands recovered over the next two millennia and by 1.5 ka the modern vegetation appears to have become established. Copyright © 2009 John Wiley & Sons, Ltd.     

The second order characteristic condition for plastic flow

In this paper the second order characteristic (discontinuous bifurcation) condition is derived for the granular flow (fully plastic) equations. This second order bifurcation equation is shown to be formally identical to the first order localization requirement during steady elastoplastic deformation provided the elastic compliance tensor is substituted for the product of the plastic multiplier with the flow Hessian. For isotropic yield and flow functions the invariant form of the characteristic condition is given in detail, as well as an alternative expression in adapted co‐ordinates. The characteristic condition can be regarded as defining a hardening function which is maximized to identify the critical angles. When the method is applied to 3D Coulomb flow, Mohr's 3D fracture plane conditions are obtained uniquely. Copyright © 2003 John Wiley & Sons, Ltd.     

Three-dimensional calculations of high- and low-mass planets embedded in protoplanetary discs

We analyse the non-linear, three-dimensional response of a gaseous, viscous protoplanetary disc to the presence of a planet of mass ranging from 1 Earth mass (1 M_⊕) to 1 Jupiter mass (1 M_J) by using the zeus hydrodynamics code. We determine the gas flow pattern, and the accretion and migration rates of the planet. The planet is assumed to be in a fixed circular orbit about the central star. It is also assumed to be able to accrete gas without expansion on the scale of its Roche radius. Only planets with masses   M _p≳ 0.1 M_J  produce significant perturbations in the surface density of the disc. The flow within the Roche lobe of the planet is fully three-dimensional. Gas streams generally enter the Roche lobe close to the disc mid-plane, but produce much weaker shocks than the streams in two-dimensional models. The streams supply material to a circumplanetary disc that rotates in the same sense as the orbit of the planet. Much of the mass supply to the circumplanetary disc comes from non-coplanar flow. The accretion rate peaks with a planet mass of approximately 0.1 M_J and is highly efficient, occurring at the local viscous rate. The migration time-scales for planets of mass less than 0.1 M_J, based on torques from disc material outside the Roche lobes of the planets, are in excellent agreement with the linear theory of type I (non-gap) migration for three-dimensional discs. The transition from type I to type II (gap) migration is smooth, with changes in migration times of about a factor of 2. Starting with a core which can undergo runaway growth, a planet can gain up to a few M_J with little migration. Planets with final masses of the order of 10 M_J would undergo large migration, which makes formation and survival difficult.