A continuum‐discrete hydromechanical analysis of granular deposit liquefaction

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.     

Weathering process at the natural fission reactor of Bangombé

 The uranium deposits in the basin of Franceville (Gabon) host the only natural fission reactors known in the world. Unique geological conditions favoured a natural fission reaction 2 Ga ago. This was detected by anomalous isotopic compositions of uranium and rare earth elements (REE), which are produced by the fission reaction. In total, 16 reactor zones were found. Most of them are mined out. The reactor zone of Bangombé, is only 10–11 m below the surface. This site has been influenced by surface weathering processes. Six drill cores have been sampled at the site of the reactor zone of Bangombé during the course of the study and only one drill core (BAX 08) hit the core of the reactor. From these data and previous drilling campaigns, the reactor size is estimated to be 10 cm thick, 2–3 m wide and 4–6 m long. The migration of fission products can be traced by the anomalous isotope ratios of REE because of the fission process. The ¹⁴⁹Sm/¹⁴⁷Sm ratio close to the reactor zone is only 0.28 (normal: 0.92) because of the intense neutron capture of ¹⁴⁹Sm and subsequent transmutation, whereas ¹⁴⁷Sm is enriched by the fission reaction. Similar changes in isotopic patterns are detectable on other REE. The isotope ratios of Sm and Nd of whole rock and fracture samples surrounding the reactor indicate that fission-genic REE migrated only a few decimetres above and mainly below the reactor zone. Organic matter (bitumen) seems to act as a trap for fission-genic REE. Additional REE-patterns show less intense weathering with increasing depth in the log profile and support a simple weathering model. Received: 26 November 1999 · Accepted: 2 May 2000     

Effect of element size on the static finite element analysis of steep slopes

The accuracy of the computed stress distribution near the free surface of vertical slopes was evaluated in this study as a function of the element size, including aspect ratio. To accomplish this objective, a parametric study was carried out comparing stresses computed using the finite element method (FEM) to those obtained from a physical model composed of photoelastic material. The results of the study indicate a reasonable agreement between a gelatin model and the FEM model for shear stresses, and an overall good agreement between the two models for the principal stresses. For stresses along the top of the slope, the height of the element tends to be more important than width or aspect ratio, at least for aspect ratios up to 4. In all cases, the greatest difference between the two models occurs in the vicinity of the slope. Specifically, if H is defined as the slope height, an element height of H/10 appears to be adequate for the study of stresses deep within the slope, such as for typical embankment analyses. However, for cases where tensile stresses in the vicinity of the slope face which are critical, such as for the stability analysis of steep slopes, element heights as small as H/32, or higher‐order elements, are necessary. Copyright © 2001 John Wiley & Sons, Ltd.     

A micromechanical finite element model for linear and damage‐coupled viscoelastic behaviour of asphalt mixture

This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd.     

Binogeochemical process of major elements in sinking particulate of Nansha coral reef lagoons,South

Ｂｉｎｏｇｅｏｃｈｅｍｉｃａｌ　ｐｒｏｃｅｓｓ　ｏｆ　ｍａｊｏｒ　ｅｌｅｍｅｎｔｓ　ｉｎ　ｓｉｎｉｎｇ　ｐａｒｔｉｃｕｌａｔｅ　ｏｆ　Ｎａｎｓｈａ　ｃｏｒａｌ　ｒｅｅｆ　ｌａｇｏｏｎｓ，Ｓｏｕｔｈ　TXＢｉｎｏｇｅｏｃｈｅｍｉｃａｌｐｒｏｃｅｓｓｏｆｍ...     

雪宝顶碱性花岗岩岩石地球化学与成矿控制

对雪宝顶碱性花岗岩具有高丰度的W、Sn、Be含量和W、Sn矿化花岗岩的微量元素组合特征进行了分析，结果显示：稀土元素具有∑REE含量与LREE／HREE比值偏低和负Eu异常的特点；元素比值K／Rb，U／Th，Nb／Ta的规律性变化以及矿脉产状特征等表明W-(Sn)-Be脉状矿床形成于岩浆演化晚期，与岩浆期后的热液活动有关。     

Growth Aspects and Magnesium Carbonate Concentrations of the Reef Building Octocoral Heliopora coerulea (Pallas) after Transplantation

Abstract. An experiment to test C have's hypothesis on temperature dependent MgCO₃ concentrations in marine calcareous skeletons was undertaken by transplanting colonies of Heliopora coerulea. The corals were transferred from a shallow lagoon with high average temperatures to the outer bottom of a fringing reef, 16 m depth, with low average temperatures. Against expectation, a statistically significant increase of average mol % concentrations from 0.48 to 0.54 % for upper growth regions and from 0.49 to 0.62% for lower growth regions was obtained. Temperature records during the experimentation period and the above data indicate that the mol% MgCO₃ concentration in the skeleton of living Heliopora colonies increases when transplanted to an environment with a lower temperature average than that of their original growth site, and with seasonal temperature fluctuations which are smaller in the experimental site than in the control site. Temperature fluctuations at the control site were up to 2.5 times greater than at the experimental site for most of the observation months. This supports the conclusion that the parameters a) minimum temperature and b) temperature stability of the environment are essential for the incorporation of Mg into the carbonate skeleton matrix of marine organisms. Within this study the results of a quantitative chemical analysis of the carbonate skeleton of H. coerulea are presented and growth data of transplanted colonies compared with those from the control site.     

Simulations of 3D DC borehole resistivity measurements with a goal-oriented <Emphasis Type="Italic">hp</Emphasis> finite-element method. Part II: through-casing resistivity instruments

We simulate direct current (DC) borehole resistivity measurements acquired in steel-cased deviated wells for the assessment of rock formation properties. The assumed data acquisition configuration considers one current (emitter) and three voltage (collector) electrodes that are utilized to measure the second difference of the electric potential along the well trajectory. We assume a homogeneous, 1.27-cm-thick steel casing with resistivity equal to 10^− 5 Ω· m. Simulations are performed with two different numerical methodologies. The first one is based on transferring two-dimensional (2D) axisymmetric optimal grids to a three-dimensional (3D) simulation software. The second one automatically produces optimal 3D grids yielded by a 3D self-adaptive goal-oriented algorithm. Both methodologies utilize high-order finite elements (FE) that are specially well-suited for problems with high-contrast coefficients and rapid spatial variations of the electric field, as it occurs in simulations that involve steel-cased wells. The method based on transferring 2D-optimal grids is efficient in terms of CPU time (few seconds per logging position). Unfortunately, it may produce inaccurate 3D simulations in deviated wells, even though the error remains below 1% for the axisymmetric (vertical) well. The method based on optimal 3D grids, although less efficient in terms of CPU time (few hours per logging position), produces more accurate results that are validated by a built-in a posteriori error estimator. This paper provides the first existing simulations of through-casing resistivity measurements in deviated wells. Simulated resistivity measurements indicate that, for a 30° deviated well, measurements in conductive layers 0.01 Ω· m) are similar to those obtained in vertical wells. However, in resistive layers (10,000 Ω· m), we observe 100% larger readings in the 30° deviated well. This difference becomes 3,000% for the case of a 60° deviated well. For this highly-deviated well, readings corresponding to the conductive formation layer are about 30% smaller in magnitude than those in a vertical well. Shoulder effects significantly vary in deviated wells.     

胶莱盆地含金建造地球化学特征

胶莱盆地蓬家夼金矿床及宋家沟金矿床均赋存于下白垩统莱阳群一段含金建造中,该含金建造由砾岩、含砾砂岩、碳质泥岩及碳酸盐岩等组成,Au100.2～241.0ng/g,Ag1.50～3.73μg/g,Zn560～1021μg/g,Pb356～364μg/g,Cu60～394μg/g,As11.0～25.4μg/g,同时富硫及有机碳。含金建造K-Rb、K-Ba呈正相关,K-Sr、K-Ca不相关,Sr/Ba<1,B、V关系均显示了陆相淡水沉积特征。La/Y-REE组成显示了沉积特征。围岩、矿石及矿石中黄铁矿及磁铁矿的稀土元素分布模式相似,该含金建造是矿源层也是容矿层,是区内重要的成矿地质条件及找矿标志。     

Post-mid-Cretaceous eastern North Sea evolution inferred from 3D thermo-mechanical modelling

The Late Cretaceous–Cenozoic evolution of the eastern North Sea region is investigated by 3D thermo-mechanical modelling. The model quantifies the integrated effects on basin evolution of large-scale lithospheric processes, rheology, strength heterogeneities, tectonics, eustasy, sedimentation and erosion.

The evolution of the area is influenced by a number of factors: (1) thermal subsidence centred in the central North Sea providing accommodation space for thick sediment deposits; (2) 250-m eustatic fall from the Late Cretaceous to present, which causes exhumation of the North Sea Basin margins; (3) varying sediment supply; (4) isostatic adjustments following erosion and sedimentation; (5) Late Cretaceous–early Cenozoic Alpine compressional phases causing tectonic inversion of the Sorgenfrei–Tornquist Zone (STZ) and other weak zones.

The stress field and the lateral variations in lithospheric strength control lithospheric deformation under compression. The lithosphere is relatively weak in areas where Moho is deep and the upper mantle warm and weak. In these areas the lithosphere is thickened during compression producing surface uplift and erosion (e.g., at the Ringkøbing–Fyn High and in the southern part of Sweden). Observed late Cretaceous–early Cenozoic shallow water depths at the Ringkøbing–Fyn High as well as Cenozoic surface uplift in southern Sweden (the South Swedish Dome (SSD)) are explained by this mechanism.

The STZ is a prominent crustal structural weakness zone. Under compression, this zone is inverted and its surface uplifted and eroded. Contemporaneously, marginal depositional troughs develop. Post-compressional relaxation causes a regional uplift of this zone.

The model predicts sediment distributions and paleo-water depths in accordance with observations. Sediment truncation and exhumation at the North Sea Basin margins are explained by fall in global sea level, isostatic adjustments to exhumation, and uplift of the inverted STZ. This underlines the importance of the mechanisms dealt with in this paper for the evolution of intra-cratonic sedimentary basins.