Analysis of the hydration of a bentonite seal in a deep radioactive waste repository

A deep geological repository for nuclear waste requires the backfilling and sealing of shafts and galleries to block any preferential path for radioactive contaminants. The paper presents the coupled hydromechanical analyses of an in situ test carried out in the HADES underground laboratory in Mol, Belgium. The test examines the effectiveness of an expansive clay seal in a horizontal borehole specifically drilled for this purpose. The analysis covers the phase of seal hydration up to saturation and subsequent pore pressure equilibration. Hydraulic and mechanical constitutive laws suited to expansive clay materials have been chosen for the analyses with all the parameters determined independently of the in situ test. A quite good agreement has been found between test observations and computed results suggesting that the numerical formulation employed is able to reproduce the main features of a real sealing system.     

Interaction between the Nazca and South American plates and formation of the Altiplano–Puna plateau: Review of a flexural analysis along the Andean margin (15°–34°S)

The results of a two-dimensional flexural analysis applied to the Andean margin, which is based on the correlation between topography and Bouguer anomaly, are here reviewed in order to characterize rigidity variations across and along the forearc–arc transition of the Central Andes and to understand the role of the forearc in the formation of the Altiplano Plateau. The forearc has maximum rigidities between 15° and 23°S. Forearc rigidity decreases gradually southward and sharply toward the plateau. The main orogen (elevations higher than 3000 m) is very weak along the entire Central Andes. A semi-quantitative interpretation of these trends, based on the relationship between flexural rigidity and the thermo-mechanically- and compositionally-controlled strength of the lithosphere, allows the following conclusions to be made: (1) across-strike rigidity variations are dominated by the thermal structure derived from the subduction process; (2) the forearc constitutes a strong, cold and rigid geotectonic element; (3) southward weakening of the forearc is directly related to the decreasing thermal age of the subducted slab; (4) very low rigidities along the main orogen are caused by the existence of a thick, quartz-rich crust with a low strain rate-to-heat flow ratio; (5) the strength of the plateau lithosphere is localized in an upper-crustal layer whose base at 15 km could be correlated with a P-to-S seismic wave converter (TRAC1 of Yuan et al., 2000 [Yuan, X., Sobolev, S., Kind, R., Oncken, O. et al. 2000. Subduction and collision processes in the Central Andes constrained by converted seismic phases. Nature, V 408, 21/28 Diciembre, p. 958–961]); (6) the forearc–plateau rigidity boundary corresponds to a zone of changing thermal conditions, eastward-increasing crustal thickness and felsic component in the crust, and low strain-rate deformation, which correlates with a west-verging structural system at the surface. These conclusions suggest that the rigid forearc acts as a pseudo-indenter against the weak plateau and allows the accumulation of ductile crustal material that moves westward from the eastern foreland. This pseudo-indenter is geometrically represented by a crustal-scale triangular zone rooted at TRAC1. This model allows the integration of existing contradictory ideas on the dynamics of forearc–plateau interaction that are related to the relative importance of upper-crustal compressive structures and lower crustal accumulation below the forearc.     

Modelling a geo-composite cell using discrete analysis

Reinforced soil structures can be used to protect infrastructures against rockfalls. An innovative way to implement such protection is to construct the front face using a cellular assembly. This paper focuses on the cellular scale, with the cell a complex heterogeneous material, composed of a wire netting box filled with rocky particles. The behaviour of a single cell was modelled using the discrete element method, thus accounting for the interaction between the rocky particles and the interaction between the box and the rocky particles. A constitutive model was developed and calibrated along confined compression loading paths, then preliminary elements of validation were obtained from the simulation of unconfined compression tests.     

Can tensile stress develop in fractured multilayers under compressive strain conditions?

To contribute to the understanding of how opening-mode fractures (joints) form and open or close at depth in layered rocks, we present a 2D numerical study aiming to determine whether tensile stress can develop in pre-fractured elastic multilayers submitted to biaxial compressive strain conditions.First, we investigate the role of the elastic and geometrical properties of the layers on the development of tensile stress in models with five bonded layers and containing one open fracture in the central layer. Our results indicate that, in absence of elastic contrast (in Young's modulus) between the layers, no tensile stress develops in the models. However, when the fractured layer is stiffer than the two adjacent layers directly above and below, a lobe of horizontal tensile stress develops centered on the pre-existing fracture. The creation of this tensile stress is contingent upon the partial closing of the fracture. The levels of tensile stress and the thickness of the lobe of tensile stress increase logarithmically with an increase in the elastic contrast and are systematically larger for a larger Soft/Stiff ratio (ratio of the total thickness of the soft layers with the total thickness of the stiff layers).Second, we investigate the role of fracture interaction in the development of tensile stress in models containing a pair of open fractures. We observe that the levels of tensile stress in the region between the fractures are systematically higher than those observed in identical models containing a single fracture. This increase in tensile stress is very large for small elastic contrasts between the layers but diminishes when the elastic contrast increases. Furthermore, the spacing between the pre-existing fractures plays an important role in the stress distribution in the region between them. When the fracture spacing is equal to or lower than 1.15 times the height of the fractured layer for the experimental conditions chosen, the lobes of tensile stress centered on the fractures coalesce. This results in the formation of vast areas of tensile stress in models under remote compressive loading conditions. Such tensile areas are likely to allow the initiation and propagation of subsequent opening-mode fractures.The results obtained provide new insights into the formation of joints in layered rocks in compressive environments, with important consequences on fluid flow.     

Salt and ice crystallisation in porous sandstones

Salt and ice crystallisation in the pore spaces causes major physical damage to natural building stones. The damaging effect of these processes can be traced back to physically induced stress inside of the rock while crystallizing. The increasing scientific research done during the past century has shown that there are numerous parameters that have an influence on the weathering resulting from these processes. However, the working mechanisms of the stress development within the rock and its material dependency are still subject to discussion. This article gives an overview of salt and ice weathering. Additionally, laboratory results of various sandstones examined are presented. Salt crystallisation tests and freeze/thaw tests were done to obtain information about how crystallisation weathering depends on material characteristics such as pore space, water transportation, and mechanical features. Simultaneous measuring of the length alternating during the salt and ice crystallisation has revealed detailed information on the development of crystal in the pore spaces as well as the development of stress. These findings can help to understand the damaging mechanisms.     

Irreversible thermodynamic and viscoelastic model for power-law relaxation and attenuation of rocks

A constitutive law for viscoelastic behaviour of rocks is derived from irreversible thermodynamics. To this model, two specific parameters are introduced; one is an internal state variable which is a variable concerning the microstructures such as defects in crystals or microcracks, and the other is a temperature reduced time obtained by normalizing the various temperature behaviours. A large number of internal state variables have the respective relaxation times and show the respective time evolutions, while a set of the time evolutions generates temporal power-law behaviour of rocks. The time evolutions of internal states are regarded as dynamics of elements of the generalized Maxwell model, and the stress–strain relation is represented by a response function following a temporal power-law in terms of linear system theory. This relation is inversely formulated to investigate the source field from output data. This model enables us to explain experimentally-based constitutive laws for transient and steady-state behaviour of rocks (e.g., lherzolite) following a temporal power-law and for attenuation behaviour of polycrystals (e.g., olivine) represented by a relation between the quality factor and frequency. Both laws show power-laws on deformation time or frequency depending on the fractal structure in polycrystals or rocks, and the experimental high-temperature behaviours can be extrapolated to long deformation time or high frequency behaviour.     

Assessing the EU ETS with a bottom-up,multi-sector model

The European Emissions Trading System (EU ETS) is the central pillar of the EU response against climate change. This trading mechanism is considered, from the theoretical point of view, as the most cost-effective method to reduce GHG. However, previous studies show that the agents who participate in these markets may behave in a way that may lead to inefficient CO₂ prices, creating doubts about the static and dynamic efficiency of the system. This article analyses these possible anomalies by first trying to model the ETS in a more realistic way, addressing some of the limitations of previous models, and second, by comparing the results with real market transactions. For this, a bottom-up, multi-sector model has been built, which represents the EU ETS in an integrated, cross-sectoral way, paying particular attention to the interactions among the most emissions intensive industries. The results show the benefits of this modelling approach and how it better reflects real market conditions. Some preliminary conclusions regarding the behaviour of the agents in the ETS market are also presented.

POLICY RELEVANCE

Low allowance prices in the EU ETS have put into question the dynamic efficiency of the EU ETS system, prompting various ideas for structural reform. However, determining the right reform also requires estimating correctly how agents will respond to it. This article proposes a tool to realistically simulate the EU ETS under the assumption of rational agents, and compare it to real market outcomes, in order to understand better the behaviour of agents in this carbon market, and therefore how to design better policies.     

原位低温拉曼光谱技术在人工合成CaCl2-H2O和MgCl2-H2O体系流体包裹体分析中的应用Ⅱ:低温下流体包裹体相变行为的研究

本文利用低温原位拉曼技术,对CaCl2-H2O体系和MgCl2-H20体系人工合成流体包裹体进行了研究.结果表明:对于盐浓度不同的溶液而言,可采用不同的冷冻方式有效采集低温拉曼光谱;通过系统采集不同温度下的拉曼光谱,可以直接准确地测定包裹体中流体的成分和低温相变过程.人工合成包裹体原位低温拉曼光谱的研究,为将该技术应用于天然包裹体分研究奠定了理论基础,可以预见,该技术必将在流体包裹体研究领域发挥其它方法不可替代的重要作用.     

考虑围岩软化特性和应力释放的圆形隧道黏弹塑性解

将围岩的塑性应变软化特性引入到考虑应力释放的圆形隧道黏弹塑性解中,并且在围岩的软化和残余强度阶段考虑围岩的塑性体积膨胀特性,提出了考虑塑性软化以及塑性体积膨胀和围岩应力释放的圆形隧道弹塑性解。当软化系数k = ∞、膨胀系数h = s时,该解转化为黏弹-脆塑性解;当k = 0、h = s时,则转化为黏弹-理想塑性解,进一步令h = s = 1,则转化为不考虑塑性体积膨胀的黏弹-理想塑性解。通过具体实例计算,分析了掌子面与研究断面间距x、围岩的软化系数k、膨胀系数h和s、支护结构等对围岩塑性区、破碎区半径和变形的影响。当开挖面与研究断面间距x在(0~4)D（D为隧道直径）范围内,随着时间增加塑性圈和破碎区迅速增大;超过4D,塑性区和破碎区半径增量逐渐变小,趋于稳定值;围岩中包含塑性区和破碎区时,二者半径的比值只取决于围岩的性质,与支护结构无关,但支护结构可以限制塑性区及破碎区的范围;考虑应变软化和塑性体积膨胀时,围岩径向位移和塑性区及破碎区半径均大于不考虑应变软化和塑性体积膨胀时的结果;软化系数k增大,围岩位移、塑性区和破碎区半径增加、塑性区半径和破碎区半径之间的比值变小。得到的结果对于隧道工程设计和施工具有一定的指导性和参考价值。