Modeling the heat development in hydrating CPB structures

Cemented paste backfill (CPB) is a mixture of dewatered tailings, hydraulic binders and water. In addition to contributing to the stability of mine workplaces, CPB greatly benefits the environment by minimizing surface tailings disposal. Hence, it has become one of the most commonly used ways in mine backfilling around the world. Temperature can significantly affect the mechanical properties of cemented backfill. A source of heat in CPB is produced by binder hydration. Hence, a FLAC based numerical model is developed to predict and analyse the heat developed by hydrating CPB structures. To validate the model, results of the developed model are compared with three case studies (mathematical, laboratory, and field investigations). The validation results show a good agreement between the developed model and these cases. The effects of stope geometry, thermal properties of both rock and CPB, filling rate, binder content and initial boundary conditions are also investigated.     

A coupled chemo‐thermo‐hygro‐mechanical model of concrete at high temperature and failure analysis

A hierarchical mathematical model for analyses of coupled chemo‐thermo‐hygro‐mechanical behaviour in concretes at high temperature is presented. The concretes are modelled as unsaturated deforming reactive porous media filled with two immiscible pore fluids, i.e. the gas mixture and the liquid mixture, in immiscible–miscible levels. The thermo‐induced desalination process is particularly integrated into the model. The chemical effects of both the desalination and the dehydration processes on the material damage and the degradation of the material strength are taken into account. The mathematical model consists of a set of coupled, partial differential equations governing the mass balance of the dry air, the mass balance of the water species, the mass balance of the matrix components dissolved in the liquid phases, the enthalpy (energy) balance and momentum balance of the whole medium mixture. The governing equations, the state equations for the model and the constitutive laws used in the model are given. A mixed weak form for the finite element solution procedure is formulated for the numerical simulation of chemo‐thermo‐hygro‐mechanical behaviours. Special considerations are given to spatial discretization of hyperbolic equation with non‐self‐adjoint operator nature. Numerical results demonstrate the performance and the effectiveness of the proposed model and its numerical procedure in reproducing coupled chemo‐thermo‐hygro‐mechanical behaviour in concretes subjected to fire and thermal radiation. Copyright © 2005 John Wiley & Sons, Ltd.     

Multiphysics modeling of arching effects in fill mass

A numerical modeling study is conducted to assess and gain a better understanding of the arching effects of field cemented tailings backfill (CTB). An integrated multiphysics model is developed that can illustrate and capture the changes in the material properties of CTB, consolidation behavior of CTB mass, and the shear behavior at the CTB/Rockwall interface. The predictive capability of the model has been successfully verified with comparisons of the predicted results with monitoring data taken from a series of field studies. The model is then used to simulate a series of applications that are relevant to CTB in practice.     

Unsaturated hydraulic properties of cemented tailings backfill that contains sodium silicate

Recycling the mine waste (tailings) into cemented tailings backfill has economical and environmental advantages for the mining industry. One of the most recent types of cemented tailings backfill is gelfill (GF), a backfill that contains sodium silicate as chemical additive. GF is typically made of tailings, water, binder and chemical additives (sodium silicate gel). It is a promising mine tailings backfill technology. From a design point of view, the environmental performance or durability of GF structures is considered as a key factor. Due to the fact that GF structures are cementitious tailings, their durability and environmental performance depend on their ability to resist the flow of aggressive elements (water and oxygen). Thus, understanding the unsaturated hydraulic properties of GF is essential for a cost-effective, environmentally friendly and durable design of GF structures. However, there is a lack of information with regards to unsaturated hydraulic properties of GF, the factors that affect them and their evolution with time. Hence, the unsaturated hydraulic properties (water retention curve (WRC) or water characteristic curve, air entry value (AEV), residual water content, unsaturated hydraulic conductivity) of GF are investigated in this paper. GF samples of various compositions and cured in room temperature for different times (3, 7, 28, and 90 days) are considered. Saturated hydraulic conductivity and microstructural tests have been conducted; WRCs are measured by using a WP4-T dewpoint potentiameter and the saline solution method. Unsaturated hydraulic conductivity is predicted using the van Genuchten (1980) equation. The water retention curve (WRC) is determined as the relationship between volumetric water content and suction for each GF mix and curing time. The van Genuchten (1980) equation is used to simulate the WRC to best-fit the experimental data. AEV and residual water content are also computed for each mix and curing time. Furthermore, functions are developed to predict the evolution of AEVs, residual water content and fitting parameters of the van Genuchten model with degree of hydration. Important outcomes have been achieved with regards to unsaturated hydraulic properties. The unsaturated hydraulic conductivity of GF was calculated to decrease when the suction, binder content, and degree of hydration increase. The effects of binder content and degree of hydration are more obvious at low suction ranges. The obtained results would contribute to a better design and assessment of the durability and environmental performance of GF structures.     

Coupled Modeling of Temperature Distribution and Evolution in Cemented Tailings Backfill Structures that Contain Mineral Admixtures

Cemented paste backfill (CPB, a mixture of tailings, water and binder) is widely utilized to fill underground mine voids. To achieve a good, economical performance, one approach is to proportionally use mineral admixtures such as fly ash and slag as partial substitutes for Portland cement. Binder hydration is one of the most significant factors that can generate heat within hydrating CPB structures, which in turn, influences the mechanical and hydraulic properties of CPB, as well as the pore structure within CPB. However, the temperature evolution due to the hydration of Portland cement that contains fly ash or slag is different from that of hydration with solely Portland cement. Hence, in consideration of the heat generated by both binder hydration and transferred between CPB and its surrounding media, a numerical model is developed to predict and determine the temperature development within CPB that contains mineral admixtures. After that, data from field and laboratory studies are employed to validate the developed model. The validation results demonstrate a good consistency between the model and the field and laboratory studies. Consequently, the proposed model is applied to simulate and determine the temperature evolution with time via mineral admixtures, binder content, initial rock and CPB temperatures, stope geometry, backfilling rate, curing time and backfilling strategy. The obtained results will contribute to better designs and preparation of CPB mixtures, as well as predict the temperature distribution within CPB structures.     

Consolidation and incipient oxidation of alkaline arsenopyrite-bearing mine tailings,Macraes Mine,New Zealand

Fine grained (ca. 15 μm), arsenopyrite-bearing mine tailings have been exposed to drying and oxidation for 4 a pending relocation. The tailings are still partly covered by a pond of decanted pore waters. The water table in drying tailings has lowered by 1–3 m and desiccation cracks up to 2 cm wide have formed on the 1 m scale, extending through the unsaturated zone. Tailings in the unsaturated zone have similar pore water contents to saturated tailings: typically 16–32 wt% water. Saturated tailings retain alkaline pH (ca. 10) from the mine cyanidation plant, but pH lowers progressively towards ca. 7 near the surface, or near desiccation cracks, in the unsaturated zone. The redox state of the tailings changes in parallel with pH, with an empirical relationship: Eh(mV)=−55 pH+290. Water in the remnant decant pond reflects this relationship also. Unsaturated tailings have variable but low permeabilities, typically 10⁻³ to 10⁻⁴ m/day, and more permeable horizons have allowed incursion of oxygenated air and/or rain water from desiccation cracks. Sulphide grains in all tailings examined are unaltered. Sulphides and solutions in the tailings are out of thermodynamic equilibrium predicted from the redox–pH conditions, due to kinetic constraints. Incursion of rain water locally facilitates deposition from pore waters of insoluble Fe oxide and arsenate minerals, thus fixing As in the dry unsaturated tailings.     

A coupled poroplastic damage model accounting for cracking effects on both hydraulic and mechanical properties of unsaturated media

Damage induced by microcracking affects not only the mechanical behaviour of geomaterials but also their hydraulic properties. Evaluating these impacts is important for many engineering applications, such as the safety assessment of radioactive waste disposal facilities. This paper presents a new constitutive model accounting simultaneously for the impact of damage on hydraulic and mechanical properties of unsaturated poroplastic geomaterials. The hydro‐mechanical coupling is formulated by means of the thermodynamic framework for partially saturated media, extended by taking into account isotropic damage and plasticity. State and complementary laws are governed by the so‐called plastic effective stress and equivalent pore pressure. Assuming a bimodal pore size distribution for cracked porous media, the hydraulic part (water retention curve and hydraulic conductivity) is modelled using phenomenological functions of damage variable. The participation of damage on both mechanical and hydraulic part enables this model to describe bilateral couplings between them. This coupled model is then validated against a number of experimental data obtained from Callovo‐Oxfordian argillite, which is the possible host rock for a radioactive waste disposal in France. Parametric studies are also carried out to check the consistency and to better demonstrate the bilateral couplings in the model. Copyright © 2015 John Wiley & Sons, Ltd.     

尾矿连续沉积过程中超孔隙水压力演化的温度-渗流-力学-化学耦合模型

尾矿胶结充填技术可避免矿渣的地表堆存、改善采场围岩稳定性、提高矿石回采率,因此在国内外的矿山开采中得到了广泛应用。充填系统的稳定性取决于尾矿连续充填过程中多场耦合作用造成的孔隙水压力演化。本研究基于经典Biot孔隙弹性理论建立尾矿的温度-渗流-力学-化学耦合模型框架,进而提出连续沉积过程中尾矿热-化学固结的一维超压模型,并推导化学反应造成水压变化的临界温度闭合公式。通过分析不同沉积速率、不同初始和边界温度条件下的充填过程,揭示了多场耦合作用对尾矿超孔隙水压力的影响机理。     

泥岩核废料处置库温度-渗流-应力耦合参数敏感性分析

针对比利时HADES地下实验室PRACLAY现场加热试验,应用温度-渗流-应力耦合弹塑性模型,模拟现场加热过程中泥岩核废料处置库的水力学响应特征。采用单因素分析法,就泥岩热、水、力学参数对核废料处置库围岩孔压、温度、有效应力的影响进行了三维有限元分析。并基于参数敏感性分析结果,就温度、渗流、应力三场两两耦合作用对处置库围岩水力学响应的影响程度进行了系统分析。研究结果表明：泥岩热、水、力学参数中,渗透系数、弹性模量以及导热系数对加温所导致的超孔压的值影响较大;凝聚力、内摩擦角以及热膨胀系数对孔压的影响较小,但会显著影响围岩的有效应力;导热系数对围岩温度场的分布有决定性影响,温度传递的差异会显著影响围岩的孔压和有效应力;不同的热、水、力学参数对孔压、温度以及有效应力的影响机制是不同的,温度、渗流、应力三场两两耦合作用对围岩水力学响应的影响程度也存在显著的差异性。温度场对应力场、温度场对渗流场的耦合效应十分显著,加热后,围岩超孔压的产生以及热膨胀导致的有效应力变化会显著影响处置库的稳定。该研究结果在一定程度上可以为核废料处置库泥岩的热、水、力学参数的确定及耦合机制分析提供科学依据。     

Consolidation in spatially random unsaturated soils based on coupled flow‐deformation simulation

This paper integrates random field simulation of soil spatial variability with numerical modeling of coupled flow and deformation to investigate consolidation in spatially random unsaturated soil. The spatial variability of soil properties is simulated using the covariance matrix decomposition method. The random soil properties are imported into an interactive multiphysics software COMSOL to solve the governing partial differential equations. The effects of the spatial variability of Young's modulus and saturated permeability together with unsaturated hydraulic parameters on the dissipation of excess pore water pressure and settlement are investigated using an example of consolidation in a saturated‐unsaturated soil column because of loading. It is found that the surface settlement and the pore water pressure profile during the process of consolidation are significantly affected by the spatially varying Young's modulus. The mean value of the settlement of the spatially random soil is more than 100% greater than that of the deterministic case, and the surface settlement is subject to large uncertainty, which implies that consolidation settlement is difficult to predict accurately based on the conventional deterministic approach. The uncertainty of the settlement increases with the scale of fluctuation because of the averaging effect of spatial variability. The effects of spatial variability of saturated permeability k_sat and air entry parameters are much less significant than that of elastic modulus. The spatial variability of air entry value parameters affects the uncertainties of settlement and excess pore pressure mostly in the unsaturated zone. Copyright © 2016 John Wiley & Sons, Ltd.     

Experiment and validation of numerical simulation of coupled thermal,hydraulic and mechanical behaviour in the engineered buffer materials

Evaluation of the coupled heat transfer, water flow and stress changes in the engineered clay barrier is an important issue in the performance assessment of the high‐level radioactive waste disposal. To demonstrate the function of the engineered barrier system, the large‐scale experiment is conducted, which is called Big Bentonite facility (BIG‐BEN). The facility consists of an electric heater surrounded by glass beads, carbon steel overpack, buffer material and man‐made rock. The buffer is a mixture of bentonite and sand. The heater is operated at 0·8 kW. Water is injected from the interface between the buffer and the man‐made rocks at the pressure of 0·05 MPa. The duration of the experiment is 20 months. The change in temperature and swelling pressure are continuously monitored and gravimetric water content is measured by sampling. The coupled thermal, hydraulic and mechanical processes are simulated with a finite element code THAMES, which can simulate the fully coupled phenomena in the saturated and unsaturated clay under anisothermal condition. To examine the validity of the code, all the parameters used in the model are evaluated from the other laboratory tests. The simulated results are compared with the measured ones without calibration of the parameter values using the results from the BIG‐BEN experiment. It can be concluded that the changes in temperature and gravimetric water content within the buffer can be simulated reasonably well and that the mechanical effect such as swelling pressure is difficult to realize. Copyright © 2000 John Wiley & Sons, Ltd.     

Temperature dependence of the reactivity of cemented paste backfill

The environmental performance of cemented paste backfill (CPB; a mixture of tailings, water and binder), which contains sulphide mineral-bearing tailings, is strongly influenced by its reactivity. However, our understanding of the reactivity of CPB under various thermal loading conditions as well as its evolution with time is limited. Hence, a laboratory investigation is conducted to study the effects of curing and ambient (atmospheric) temperatures on the reactivity of CPB. Oxygen consumption (OC) tests are conducted on CPB specimens cured at different temperatures to study their reactivity. Furthermore, microstructural analyses (e.g., x-ray diffraction (XRD), mercury intrusion porosimetry, and thermogravimetry/derivative thermogravimetry) are performed to assess the microstructural characteristics of the tested CPBs. The results show that the reactivity of CPB is temperature-dependent. As the curing temperature increases, the reactivity generally decreases. The reactivity is also affected by the ambient temperature. The reactivity increases as the atmospheric temperature increases. However, the extent of the effect of the temperature depends on the curing time and is generally more pronounced at the early ages. Furthermore, the presence of sulphate in the pore water of CPB can significantly affect the reactivity of CPB cured at high temperatures (50 °C). The findings of this study will therefore help to better assess and predict the environmental behavior of CPB under various field thermal conditions.     

Experimental soft‐sediment deformation caused by fluidization and intrusive ice melt in sand

Identifying the driving mechanisms of soft‐sediment deformation in the geological record is the subject of debate. Thawing of ice‐rich clayey silt above permafrost was proved experimentally to be among the processes capable of triggering deformation. However, previous work has failed so far to reproduce similar structures in sand. This study investigates fluidization and intrusive ice formation from soil models in the laboratory. Experimental conditions reproduce the growth of ice‐cored mounds caused by pore water pressure increase during freeze‐back of sand in a permafrost context. Excess pore water pressure causes hydraulic fracturing and the development of water lenses beneath the freezing front. Later freezing of the water lenses generates intrusive ice. The main structures consist of sand dykes and sills formed when the increase in pore water pressure exceeds a critical threshold, and soft‐sediment deformations induced by subsidence during ice melt. The combination of processes has resulted in diapir‐like structures. The experimental structures are similar to those described in Pleistocene sites from France. These processes constitute a credible alternative to the seismic hypothesis evoked to explain soft‐sediment deformation structures in other European regions subjected to Pleistocene cold climates.     

尾矿渗透破坏及其导波监测试验研究

采用自主设计的尾矿渗透破坏仪,模拟上游水位不断升高,水力梯度持续增加的情况下尾矿发生渗透破坏的全过程。总结了尾矿渗透破坏现象和演化过程,分析了尾矿中的孔隙水压力和弯曲导波参数随时间的变化规律,研究了尾矿渗透破坏过程中弯曲导波信号b值和分形维数的演化规律。根据孔隙水压力和弯曲导波特征参数的变化特征,提出尾矿渗透破坏等级的评价标准。结果表明：（1）在尾矿渗透破坏过程中,孔隙水压力随着水力梯度的增加而增加,但当水力梯度增加到一定程度时,孔隙水压力会出现突降现象。（2）尾矿渗透破坏是一个循序渐进的过程,在这一过程中的弯曲导波参数表现为基本无信号、微弱信号、强信号和剧烈信号,分别对应尾矿正常运行阶段、发展阶段、破坏前期和最终渗透破坏阶段。（3）导波的b值和分形维数在尾矿最终渗透破坏前一直处于低水平状态,在临近渗透破坏时剧烈震荡。（4）通过对弯曲导波参数分析,提出了尾矿渗透破坏程度的量化指标,将尾矿渗透破坏分为4个预警级别,并用蓝、黄、橙、红4种警报颜色对应相应级别的预警。     

A three‐phase thermo‐hydro‐mechanical finite element model for freezing soils

Artificial ground freezing (AGF) is a commonly used technique in geotechnical engineering for ground improvement such as ground water control and temporary excavation support during tunnel construction in soft soils. The main potential problem connected with this technique is that it may produce heave and settlement at the ground surface, which may cause damage to the surface infrastructure. Additionally, the freezing process and the energy needed to obtain a stable frozen ground may be significantly influenced by seepage flow. Evidently, safe design and execution of AGF require a reliable prediction of the coupled thermo‐hydro‐mechanical behavior of freezing soils. With the theory of poromechanics, a three‐phase finite element soil model is proposed, considering solid particles, liquid water, and crystal ice as separate phases and mixture temperature, liquid pressure, and solid displacement as the primary field variables. In addition to the volume expansion of water transforming into ice, the contribution of the micro‐cryo‐suction mechanism to the frost heave phenomenon is described in the model using the theory of premelting dynamics. Through fundamental physical laws and corresponding state relations, the model captures various couplings among the phase transition, the liquid transport within the pore space, and the accompanying mechanical deformation. The verification and validation of the model are accomplished by means of selected analyses. An application example is related to AGF during tunnel excavation, investigating the influence of seepage flow on the freezing process and the time required to establish a closed supporting frozen arch. Copyright © 2013 John Wiley & Sons, Ltd.     

Desiccation shrinkage of non‐clayey soils: multiphysics mechanisms and a microstructural model

Analysis of macroscopic desiccation shrinkage experiments indicates that most, but not all of the shrinkage during drying occurs while soil is still saturated. Shrinkage practically ceases and air starts to penetrate the soil, when the water content is still quite high, for example, above 20% for the tested soils. The remaining, unsaturated drying process occurs with a much‐reduced shrinkage rate. In this context, we examine data of the pore system evolution as represented by the mercury porosimetry experimental results. The process is then modeled as a two‐stage process of deformation and evacuation of a two‐tube vessel system driven by the external evaporation flux. In the first stage, Poiseuille flow occurs through the vessels. The amount of water evaporated in this stage equals to the reduction of volume of the vessel through the deformation of its walls. This stage ends when a negative water pressure (suction) required to further deform the vessel reaches a critical value at which air enters the pore space. Two physical interpretation of such threshold are discussed. In the subsequent stage, evaporation proceeds with a receding liquid/vapor interface starting from the open end, incrementally emptying the vessel but with a marginal water flow and vessel deformation. The leading variables of the process are identified, and a quantifiable multiphysics meso‐scale scenario of models is established. Copyright © 2012 John Wiley & Sons, Ltd.     

Semi‐analytical solution to one‐dimensional consolidation for unsaturated soils with semi‐permeable drainage boundary under time‐dependent loading

This paper presents semi‐analytical solutions to Fredlund and Hasan's one‐dimensional consolidation of unsaturated soils with semi‐permeable drainage boundary under time‐dependent loadings. Two variables are introduced to transform two coupled governing equations of pore‐water and pore‐air pressures into an equivalent set of partial differential equations, which are easily solved by the Laplace transform. The pore‐water pressure, pore‐air pressure and settlement are obtained in the Laplace domain. Crump's method is adopted to perform the inverse Laplace transform in order to obtain semi‐analytical solutions in time domain. It is shown that the present solutions are more general and have a good agreement with the existing solutions from literatures. Furthermore, the current solutions can also be degenerated into conventional solutions to one‐dimensional consolidation of unsaturated soils with homogeneous boundaries. Finally, several numerical examples are provided to illustrate consolidation behavior of unsaturated soils under four types of time‐dependent loadings, including instantaneous loading, ramp loading, exponential loading and sinusoidal loading. Parametric studies are illustrated by variations of pore‐air pressure, pore‐water pressure and settlement at different values of the ratio of air–water permeability coefficient, depth and loading parameters. Copyright © 2017 John Wiley & Sons, Ltd.     

Some numerical issues using element‐free Galerkin mesh‐less method for coupled hydro‐mechanical problems

A new formulation of the element‐free Galerkin (EFG) method is developed for solving coupled hydro‐mechanical problems. The numerical approach is based on solving the two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Spatial variables in the weak form, i.e. displacement increment and pore water pressure increment, are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on a penalty method. Numerical stability of the developed formulation is examined in order to achieve appropriate accuracy of the EFG solution for coupled hydro‐mechanical problems. Examples are studied and compared with closed‐form or finite element method solutions to demonstrate the validity of the developed model and its capabilities. The results indicate that the EFG method is capable of handling coupled problems in saturated porous media and can predict well both the soil deformation and variation of pore water pressure over time. Some guidelines are proposed to guarantee the accuracy of the EFG solution for coupled hydro‐mechanical problems. Copyright © 2008 John Wiley & Sons, Ltd.     

Influence of gas generation in electro‐osmosis consolidation

This paper presents a numerical model for the elasto‐plastic electro‐osmosis consolidation of unsaturated clays experiencing large strains, by considering electro‐osmosis and hydro‐mechanical flows in a deformable multiphase porous medium. The coupled governing equations involving the pore water flow, pore gas flow, electric flow and mechanical deformation in unsaturated clays are derived within the framework of averaging theory and solved numerically using finite elements. The displacements of the solid phase, the pressure of the water phase, the pressure of the gas phase and the electric potential are taken as the primary unknowns in the proposed model. The nonlinear variation of transport parameters during electro‐osmosis consolidation are incorporated into the model using empirical expressions that strongly depend on the degree of water saturation, whereas the Barcelona Basic Model is employed to simulate the elasto‐plastic mechanical behaviour of unsaturated clays. The accuracy of the proposed model is evaluated by validating it against two well‐known numerical examples, involving electro‐osmosis and unsaturated soil behaviour respectively. Two further examples are then investigated to study the capability of the computational algorithm in modelling multiphase flow in electro‐osmosis consolidation. Finally, the effects of gas generation at the anode, the deformation characteristics, the degree of saturation and the time dependent evolution of the excess pore pressure are discussed. Copyright © 2016 John Wiley & Sons, Ltd.