An elastoplastic constitutive model based on the Modified Cam Clay(MCC)model is developed to describe the mechanical behaviour of soils cemented via microbially induced calcite precipitation(MICP).It considers the inc...An elastoplastic constitutive model based on the Modified Cam Clay(MCC)model is developed to describe the mechanical behaviour of soils cemented via microbially induced calcite precipitation(MICP).It considers the increase of the elastic stiffness,the change of the yield surface due to MICP cementation and the degradation of calcium carbonate bonds during shearing.Specifically,to capture the typical contraction-dilation transition in MICP soils,the original volumetric hardening rule in the MCC model is modified to a combined deviatoric and volumetric hardening rule.The model could reproduce a series of drained triaxial tests on MICP-treated soils with different calcium carbonate contents.Further,we carry out a parametric study and observe numerical instability in some cases.In combination with an analytical analysis,our numerical modelling has identified the benefits and limitations of using MCCbased models in the simulation of MICP-cemented soils,leading to suggestions for further model development.展开更多
The study presents a comprehensive coupled thermo-bio-chemo-hydraulic(T-BCH)modeling framework for stabilizing soils using microbially induced calcite precipitation(MICP).The numerical model considers relevant multiph...The study presents a comprehensive coupled thermo-bio-chemo-hydraulic(T-BCH)modeling framework for stabilizing soils using microbially induced calcite precipitation(MICP).The numerical model considers relevant multiphysics involved in MICP,such as bacterial ureolytic activities,biochemical reactions,multiphase and multicomponent transport,and alteration of the porosity and permeability.The model incorporates multiphysical coupling effects through well-established constitutive relations that connect parameters and variables from different physical fields.It was implemented in the open-source finite element code OpenGeoSys(OGS),and a semi-staggered solution strategy was designed to solve the couplings,allowing for flexible model settings.Therefore,the developed model can be easily adapted to simulate MICP applications in different scenarios.The numerical model was employed to analyze the effect of various factors,including temperature,injection strategies,and application scales.Besides,a TBCH modeling study was conducted on the laboratory-scale domain to analyze the effects of temperature on urease activity and precipitated calcium carbonate.To understand the scale dependency of MICP treatment,a large-scale heterogeneous domain was subjected to variable biochemical injection strategies.The simulations conducted at the field-scale guided the selection of an injection strategy to achieve the desired type and amount of precipitation.Additionally,the study emphasized the potential of numerical models as reliable tools for optimizing future developments in field-scale MICP treatment.The present study demonstrates the potential of this numerical framework for designing and optimizing the MICP applications in laboratory-,prototype-,and field-scale scenarios.展开更多
As is known, high-level radioactive waste （HLW） is commonly heat-emitting. Heat output from HLWwilldissipate through the surrounding rocks and induce complex thermo-hydro-mechanical-chemical（THMC） processes. In hi...As is known, high-level radioactive waste （HLW） is commonly heat-emitting. Heat output from HLWwilldissipate through the surrounding rocks and induce complex thermo-hydro-mechanical-chemical（THMC） processes. In highly consolidated clayey rocks, thermal effects are particularly significantbecause of their very low permeability and water-saturated state. Thermal impact on the integrity of thegeological barriers is of most importance with regard to the long-term safety of repositories. This studyfocuses on numerical analysis of thermal effects on hydro-mechanical properties of clayey rock using acoupled thermo-mechanical multiphase flow （TH2M） model which is implemented in the finite elementprogramme OpenGeoSys （OGS）. The material properties of the numerical model are characterised by atransversal isotropic elastic model based on Hooke＇s law, a non-isothermal multiphase flow model basedon van Genuchten function and Darcy＇s law, and a transversal isotropic heat transport model based onFourier＇s law. In the numerical approaches, special attention has been paid to the thermal expansion ofthree different phases： gas, fluid and solid, which could induce changes in pore pressure and porosity.Furthermore, the strong swelling and shrinkage behaviours of clayey material are also considered in thepresent model. The model has been applied to simulate a laboratory heating experiment on claystone.The numerical model gives a satisfactory representation of the observed material behaviour in thelaboratory experiment. The comparison of the calculated results with the laboratory findings verifies thatthe simulation with the present numerical model could provide a deeper understanding of the observedeffects. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
Dilatancy-controlled gas flow in preferential pathways plays a key role in the safety analysis of radioactive waste repositories.This is particularly the case for bentonite,an often-preferred barrier material.Gas flow...Dilatancy-controlled gas flow in preferential pathways plays a key role in the safety analysis of radioactive waste repositories.This is particularly the case for bentonite,an often-preferred barrier material.Gas flow in preferential pathways is characterized by localization and spontaneous behavior,which is challenging to simulate in numerical models due to strong hydro-mechanical coupling.To analyze a laboratory experiment in the framework of the DECOVALEX-2023 project,this study introduced a new approach of combining continuous modelling methods with spatial material properties derived from material heterogeneities and experimental observations.The proposed model utilized hydro-mechanical spatial distributions,namely Young’s modulus and gas entry pressure,and elastoplasticity combined with a linear swelling model.A conceptual strain-dependent permeability approach simulated dilatancycontrolled gas flow based on hydro-mechanical coupling.To test the effectiveness of the presented approach,a gas injection test in a compacted,saturated bentonite sample was simulated using the opensource code OpenGeoSys 5.8 and compared with experimental observations.The presented methodology is capable of simulating localized gas flow in preferential pathways.The spatial distributions of Young’s modulus and gas entry pressure affect the swelling pressure,relative permeability and,in combination with the strain-dependent permeability model,also the intrinsic permeability.展开更多
Simulation of subsurface energy system involves multi-physical processes such as thermal,hydraulical,and mechanical(THM)processes,and requires a so-called THM coupled modeling approach.THM coupled modeling is commonly...Simulation of subsurface energy system involves multi-physical processes such as thermal,hydraulical,and mechanical(THM)processes,and requires a so-called THM coupled modeling approach.THM coupled modeling is commonly performed in geothermal energy production.However,for hydrocarbon extraction,we need to consider multiphase flow additionally.In this paper,we describe a three-dimensional numerical model of non-isothermal two-phase flow in the deformable porous medium by integrating governing equations of two-phase mixture in the porous media flow in the reservoir.To account for inter-woven impacts in subsurface condi-tions,we introduced a temperature-dependent fluid viscosity and a fluid density along with a strain-dependent reservoir permeability.Subsequently,we performed numerical experiments of a ten-year water flooding pro-cess employing the open-source parallelized code,OpenGeoSys.We considered different well patterns with colder water injection in realistic scenarios.Our results demonstrate that our model can simulate complex interactions of temperature,pore pressure,subsurface stress and water saturation simultaneously to evaluate the recovery per-formance.High temperature can promote fluid flow while cold water injection under non-isothermal conditions causes the normal stress reduction by significant thermal stress.Under different well patterns the displacement efficiency will be changed by the relative location between injection and production wells.This finding has provided the important reference for fluid flow and induced stress evolution during hydrocarbon exploitation under the environment of large reservoir depth and high temperature.展开更多
Due to extensive water pollution in Chinese rivers and lakes,large efforts have to be made to improve the quality of drinking water and manage the sewage water treatment process.We propose a general workflow for integ...Due to extensive water pollution in Chinese rivers and lakes,large efforts have to be made to improve the quality of drinking water and manage the sewage water treatment process.We propose a general workflow for integrating a large number of heterogeneous data sets relating to various hydrological compartments into a Virtual Geographic Environment(VGE).This allows both researchers and stakeholders to easily access complex data collections in a unified context,find interrelations or inconsistencies between data sets and evaluate simulation results with respect to other observations or simulations in the same region.A prototype of such a VGE has been set up for the region around Chao Lake,containing more than 20 spatial data sets and collections as well as first simulation result.The prototype has been successfully presented to researchers and stakeholders from China and Germany.展开更多
基金funded by the German Research Foundation(DFG)(Grant No.NA 330/20e1).
文摘An elastoplastic constitutive model based on the Modified Cam Clay(MCC)model is developed to describe the mechanical behaviour of soils cemented via microbially induced calcite precipitation(MICP).It considers the increase of the elastic stiffness,the change of the yield surface due to MICP cementation and the degradation of calcium carbonate bonds during shearing.Specifically,to capture the typical contraction-dilation transition in MICP soils,the original volumetric hardening rule in the MCC model is modified to a combined deviatoric and volumetric hardening rule.The model could reproduce a series of drained triaxial tests on MICP-treated soils with different calcium carbonate contents.Further,we carry out a parametric study and observe numerical instability in some cases.In combination with an analytical analysis,our numerical modelling has identified the benefits and limitations of using MCCbased models in the simulation of MICP-cemented soils,leading to suggestions for further model development.
基金support from the OpenGeoSys communitypartially funded by the Prime Minister Research Fellowship,Ministry of Education,Government of India with the project number SB21221901CEPMRF008347.
文摘The study presents a comprehensive coupled thermo-bio-chemo-hydraulic(T-BCH)modeling framework for stabilizing soils using microbially induced calcite precipitation(MICP).The numerical model considers relevant multiphysics involved in MICP,such as bacterial ureolytic activities,biochemical reactions,multiphase and multicomponent transport,and alteration of the porosity and permeability.The model incorporates multiphysical coupling effects through well-established constitutive relations that connect parameters and variables from different physical fields.It was implemented in the open-source finite element code OpenGeoSys(OGS),and a semi-staggered solution strategy was designed to solve the couplings,allowing for flexible model settings.Therefore,the developed model can be easily adapted to simulate MICP applications in different scenarios.The numerical model was employed to analyze the effect of various factors,including temperature,injection strategies,and application scales.Besides,a TBCH modeling study was conducted on the laboratory-scale domain to analyze the effects of temperature on urease activity and precipitated calcium carbonate.To understand the scale dependency of MICP treatment,a large-scale heterogeneous domain was subjected to variable biochemical injection strategies.The simulations conducted at the field-scale guided the selection of an injection strategy to achieve the desired type and amount of precipitation.Additionally,the study emphasized the potential of numerical models as reliable tools for optimizing future developments in field-scale MICP treatment.The present study demonstrates the potential of this numerical framework for designing and optimizing the MICP applications in laboratory-,prototype-,and field-scale scenarios.
基金supported by BMWi (Bundesministerium für Wirtschaft und Energie,Berlin)
文摘As is known, high-level radioactive waste （HLW） is commonly heat-emitting. Heat output from HLWwilldissipate through the surrounding rocks and induce complex thermo-hydro-mechanical-chemical（THMC） processes. In highly consolidated clayey rocks, thermal effects are particularly significantbecause of their very low permeability and water-saturated state. Thermal impact on the integrity of thegeological barriers is of most importance with regard to the long-term safety of repositories. This studyfocuses on numerical analysis of thermal effects on hydro-mechanical properties of clayey rock using acoupled thermo-mechanical multiphase flow （TH2M） model which is implemented in the finite elementprogramme OpenGeoSys （OGS）. The material properties of the numerical model are characterised by atransversal isotropic elastic model based on Hooke＇s law, a non-isothermal multiphase flow model basedon van Genuchten function and Darcy＇s law, and a transversal isotropic heat transport model based onFourier＇s law. In the numerical approaches, special attention has been paid to the thermal expansion ofthree different phases： gas, fluid and solid, which could induce changes in pore pressure and porosity.Furthermore, the strong swelling and shrinkage behaviours of clayey material are also considered in thepresent model. The model has been applied to simulate a laboratory heating experiment on claystone.The numerical model gives a satisfactory representation of the observed material behaviour in thelaboratory experiment. The comparison of the calculated results with the laboratory findings verifies thatthe simulation with the present numerical model could provide a deeper understanding of the observedeffects. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
基金This research was conducted within the DECOVALEX-2023 projectDECOVALEX is an international research project comprising participants from industry,government,and academia,focusing on development of understanding,models and codes in complex coupled problems in sub-surface geological and engineering applications.DECOVALEX-2023 is the current phase of the project.The authors appreciate the DECOVALEX-2023 Funding Organisations Andra,BASE,BGE,BGR,CAS,CNSC,COVRA,US DOE,ENRESA,ENSI,JAEA,KAERI,NWMO,RWM,SÚRAO,SSM and Taipower for their financial and technical support of the work described in this paper.The statements made in the paper are,however,solely those of the authors and do not necessarily reflect those of the Funding Organisations.This work was further supported by the German Federal Ministry for Economic Affairs and Climate Action(BMWK).
文摘Dilatancy-controlled gas flow in preferential pathways plays a key role in the safety analysis of radioactive waste repositories.This is particularly the case for bentonite,an often-preferred barrier material.Gas flow in preferential pathways is characterized by localization and spontaneous behavior,which is challenging to simulate in numerical models due to strong hydro-mechanical coupling.To analyze a laboratory experiment in the framework of the DECOVALEX-2023 project,this study introduced a new approach of combining continuous modelling methods with spatial material properties derived from material heterogeneities and experimental observations.The proposed model utilized hydro-mechanical spatial distributions,namely Young’s modulus and gas entry pressure,and elastoplasticity combined with a linear swelling model.A conceptual strain-dependent permeability approach simulated dilatancycontrolled gas flow based on hydro-mechanical coupling.To test the effectiveness of the presented approach,a gas injection test in a compacted,saturated bentonite sample was simulated using the opensource code OpenGeoSys 5.8 and compared with experimental observations.The presented methodology is capable of simulating localized gas flow in preferential pathways.The spatial distributions of Young’s modulus and gas entry pressure affect the swelling pressure,relative permeability and,in combination with the strain-dependent permeability model,also the intrinsic permeability.
基金supported by the National Key Research and Development Project(No.2023YFE0110900)National Natural Science Foundation of China(No.42320104003,42077247).
文摘Simulation of subsurface energy system involves multi-physical processes such as thermal,hydraulical,and mechanical(THM)processes,and requires a so-called THM coupled modeling approach.THM coupled modeling is commonly performed in geothermal energy production.However,for hydrocarbon extraction,we need to consider multiphase flow additionally.In this paper,we describe a three-dimensional numerical model of non-isothermal two-phase flow in the deformable porous medium by integrating governing equations of two-phase mixture in the porous media flow in the reservoir.To account for inter-woven impacts in subsurface condi-tions,we introduced a temperature-dependent fluid viscosity and a fluid density along with a strain-dependent reservoir permeability.Subsequently,we performed numerical experiments of a ten-year water flooding pro-cess employing the open-source parallelized code,OpenGeoSys.We considered different well patterns with colder water injection in realistic scenarios.Our results demonstrate that our model can simulate complex interactions of temperature,pore pressure,subsurface stress and water saturation simultaneously to evaluate the recovery per-formance.High temperature can promote fluid flow while cold water injection under non-isothermal conditions causes the normal stress reduction by significant thermal stress.Under different well patterns the displacement efficiency will be changed by the relative location between injection and production wells.This finding has provided the important reference for fluid flow and induced stress evolution during hydrocarbon exploitation under the environment of large reservoir depth and high temperature.
基金provided by the German Federal Ministry of Education and Research(BMBF)CLIENT program‘International Partnerships for Sustainable Innovations’for the project‘Managing Water Resources for Urban Catchments’[grant number 02WCL1337A]the Helmholtz Research Network‘Research Centre for Environmental Information Science’[grant number HIRN0001].
文摘Due to extensive water pollution in Chinese rivers and lakes,large efforts have to be made to improve the quality of drinking water and manage the sewage water treatment process.We propose a general workflow for integrating a large number of heterogeneous data sets relating to various hydrological compartments into a Virtual Geographic Environment(VGE).This allows both researchers and stakeholders to easily access complex data collections in a unified context,find interrelations or inconsistencies between data sets and evaluate simulation results with respect to other observations or simulations in the same region.A prototype of such a VGE has been set up for the region around Chao Lake,containing more than 20 spatial data sets and collections as well as first simulation result.The prototype has been successfully presented to researchers and stakeholders from China and Germany.
