Damage caused by frost heave leads to costly maintenance in cold regions, like Hokkaido, Japan. Therefore, thestudy of the frost mechanism with experimental and numerical methods has been of great interest. Numerousmo...Damage caused by frost heave leads to costly maintenance in cold regions, like Hokkaido, Japan. Therefore, thestudy of the frost mechanism with experimental and numerical methods has been of great interest. Numerousmodels have been developed to describe the freezing process of saturated soil, which differs from the partiallysaturated conditions in the field. In fact, most subsurface soils are unsaturated. The freezing process of partiallysaturated soils is more complex than saturated soils, as the governing equations show strongly nonlinear characteristics. This study proposes a thermo-hydro-mechanical coupled model considering the heat transfer, waterinfiltration, and deformation of partially saturated soil to reproduce the freezing process of partially saturatedfrost susceptible soils distributed in Hokkaido. This model better considers the water-ice phase change and soilfreezing characteristic curve (SFCC) during freezing under field conditions. The results from the multiphysicssimulations agree well with the frost heave and water migration data from frost heave tests of Touryo soil andFujinomori soil. In addition, this study discussed the influence of the various factors on frost heave amount,including temperature gradients, overburden pressures, water supply conditions, cooling rates, and initial saturation. The simulation results indicate that the frost heave ratio is proportional to the initial degree of saturationand is inversely proportional to the cooling rate and overburden pressure.Moreover, simulation under the open system generates much more frost heave than under the closed system.Finally, the main features of the proposed model are revealed by simulating a closed-system frost heave test. Thesimulation results indicate that the proposed model adequately captures the coupling characteristics of water andice redistribution, temperature development, hydraulic conductivity, and suction in the freezing process. Togetherwith the decreased hydraulic conductivity, the increased suction controls the water flow in the freezing zone. Theinflow water driven by cryogenic suction gradient feeds the ice formation, leads to a rapid increase in total watercontent, expanding the voids that exceed the initial porosity and contributing to the frost heave.展开更多
Based on fluid mechanics, thermodynamics and damage mechanics, thermal-hydro-mechanical(THM) coupling damage model of brittle rock is established by analyzing THM coupling mechanism, where THM coupling damage variable...Based on fluid mechanics, thermodynamics and damage mechanics, thermal-hydro-mechanical(THM) coupling damage model of brittle rock is established by analyzing THM coupling mechanism, where THM coupling damage variable DTHM is dominated by TH coupling damage variable DTH, TM coupling damage variable DTM and HM coupling damage variable DHM, and DTH is firstly expressed in term of dimensionless total thermal conductivity of the water Nu. Permeability test, uni-axial compression test and THM coupling test are conducted to measure the permeability, elastic modulus and THM coupling stress-strain curves of brittle rock. The tested values of THM coupling elastic modulus THME? are in good agreement with the predicted values of THM coupling elastic modulus ETHM, which can verify the newly established THM coupling damage model.展开更多
One of the most suitable ways under study for the disposal of high-level radioactive waste (HLW) is isolation in deep geological repositories. It is very important to research the thermo-hydro- mechanical (THM) coupl...One of the most suitable ways under study for the disposal of high-level radioactive waste (HLW) is isolation in deep geological repositories. It is very important to research the thermo-hydro- mechanical (THM) coupled processes associated with an HLW disposal repository. Non-linear coupled equations, which are used to describe the THM coupled process and are suited to saturated-unsaturated porous media, are presented in this paper. A numerical method to solve these equations is put forward, and a finite element code is developed. This code is suited to the plane strain or axis-symmetry problem. Then this code is used to simulate the THM coupled process in the near field of an ideal disposal repository. The temperature vs. time, hydraulic head vs. time and stress vs. time results show that, in this assumed condition, the impact of temperature is very long (over 10 000 a) and the impact of the water head is short (about 90 d). Since the stress is induced by temperature and hydraulic head in this condition, the impact time of stress is the same as that of temperature. The results show that THM coupled processes are very important in the safety analysis of an HLW deep geological disposal repository.展开更多
In examining potential host rocks for such purposes as the disposal of high-level radioactive wastes,it is important to understand the coupled thermo-hydro-mechanical（THM） behavior of a porous medium.A rigorous and ...In examining potential host rocks for such purposes as the disposal of high-level radioactive wastes,it is important to understand the coupled thermo-hydro-mechanical（THM） behavior of a porous medium.A rigorous and fully unified coupled thermo-hydro-mechanical model for unsaturated porous media is required to simulate the complex coupling mechanisms involved.Based on modified Darcy＇s and Fourier＇s laws,equations of mechanical equilibrium,mass conservation and energy conservation are derived by introducing void ratio and volumetric liquid water content into the model.The newly derived model takes into account the effects of temperature on the dynamic viscosity of liquid water and void ratio,the influence of liquid flow on temperature gradient（thermo-osmosis）,the influence on mass and heat conservation equations,and the influence of heat flow on water pressure gradient and thermal convection.The new coupled THM constitutive model is constructed by a finite element program and is used to simulate the coupled behavior of a tunnel during excavation,ventilation and concrete lining stages.Oil and gas engineering,underground disposal of nuclear waste and tunnel engineering may be benefited from the development of the new model.展开更多
As one of the most important ways to reduce the greenhouse gas emission,carbon dioxide（CO2）enhanced gas recovery（CO2-EGR） is attractive since the gas recovery can be enhanced simultaneously with CO2sequestration.B...As one of the most important ways to reduce the greenhouse gas emission,carbon dioxide（CO2）enhanced gas recovery（CO2-EGR） is attractive since the gas recovery can be enhanced simultaneously with CO2sequestration.Based on the existing equation of state（EOS） module of TOUGH2 MP,extEOS7C is developed to calculate the phase partition of H2O-CO2-CH4-NaCl mixtures accurately with consideration of dissolved NaCI and brine properties at high pressure and temperature conditions.Verifications show that it can be applied up to the pressure of 100 MPa and temperature of 150℃.The module was implemented in the linked simulator TOUGH2MP-FLAC3 D for the coupled hydro-mechanical simulations.A simplified three-dimensional（3D）1/4 model（2.2 km×1 km×1 km） which consists of the whole reservoir,caprock and baserock was generated based on the geological conditions of a gas field in the North German Basin.The simulation results show that,under an injection rate of 200,000 t/yr and production rate of 200,000 sm3/d,CO2breakthrough occurred in the case with the initial reservoir pressure of 5 MPa but did not occur in the case of 42 MPa.Under low pressure conditions,the pressure driven horizontal transport is the dominant process;while under high pressure conditions,the density driven vertical flow is dominant.Under the considered conditions,the CO2-EGR caused only small pressure changes.The largest pore pressure increase（2 MPa） and uplift（7 mm） occurred at the caprock bottom induced by only CO2injection.The caprock had still the primary stress state and its integrity was not affected.The formation water salinity and temperature variations of ±20℃ had small influences on the CO2-EGR process.In order to slow down the breakthrough,it is suggested that CO2-EGR should be carried out before the reservoir pressure drops below the critical pressure of CO2.展开更多
Geo-energy and geo-engineering applications,such as improved oil recovery(IOR),geologic carbon storage,and enhanced geothermal systems(EGSs),involve coupled thermo-hydro-mechanical(THM)processes that result from fluid...Geo-energy and geo-engineering applications,such as improved oil recovery(IOR),geologic carbon storage,and enhanced geothermal systems(EGSs),involve coupled thermo-hydro-mechanical(THM)processes that result from fluid injection and production.In some cases,reservoirs are highly fractured and the geomechanical response is controlled by fractures.Therefore,fractures should explicitly be included into numerical models to realistically simulate the THM responses of the subsurface.In this study,we perform coupled THM numerical simulations of water injection into naturally fractured reservoirs(NFRs)using CODE_BRIGHT and TOUGH-UDEC codes.CODE_BRIGHT is a finite element method(FEM)code that performs fully coupled THM analysis in geological media and TOUGH-UDEC sequentially solves coupled THM processes by combining a finite volume method(FVM)code that solves nonisothermal multiphase flow(TOUGH2)with a distinct element method(DEM)code that solves the mechanical problem(UDEC).First,we validate the two codes against a semi-analytical solution for water injection into a single deformable fracture considering variable permeability based on the cubic law.Then,we compare simulation results of the two codes in an idealized conceptual model that includes one horizontal fracture and in a more realistic model with multiple fractures.Each code models fractures differently.UDEC calculates fracture deformation from the fracture normal and shear stiffnesses,while CODE_BRIGHT treats fractures as equivalent porous media and uses the equivalent Young’s modulus and Poisson’s ratio of the fracture.Finally,we obtain comparable results of pressure,temperature,stress and displacement distributions and evolutions for the single horizontal fracture model.Despite some similarities,the two codes provide increasingly different results as model complexity increases.These differences highlight the challenging task of accurately modeling coupled THM processes in fractured media given their high nonlinearity.展开更多
The main objective of this paper is to investigate and analyse the thermo-hydro-mechanical（THM） coupling phenomena and their influences on the repository safety.In this paper,the high-level waste（HLW） disposal con...The main objective of this paper is to investigate and analyse the thermo-hydro-mechanical（THM） coupling phenomena and their influences on the repository safety.In this paper,the high-level waste（HLW） disposal concept in drifts in clay formation with backfilled bentonite buffer is represented numerically using the CODE BRIGHT developed by the Technical University of Catalonia in Barcelona.The parameters of clay and bentonite used in the simulation are determined by laboratory and in situ experiments.The calculation results are presented to show the hydro-mechanical（HM） processes during the operation phase and the THM processes in the after-closure phase.According to the simulation results,the most probable critical processes for the disposal project have been represented and analyzed.The work also provides an input for additional development regarding the design,assessment and validation of the HLW disposal concept.展开更多
基金This research was supported in part by Grant-in-Aids for Scientific Research(A,16H02360)and(B,17H03307)from the Japan Society for the Promotion of Science(JSPS)KAKENHI.
文摘Damage caused by frost heave leads to costly maintenance in cold regions, like Hokkaido, Japan. Therefore, thestudy of the frost mechanism with experimental and numerical methods has been of great interest. Numerousmodels have been developed to describe the freezing process of saturated soil, which differs from the partiallysaturated conditions in the field. In fact, most subsurface soils are unsaturated. The freezing process of partiallysaturated soils is more complex than saturated soils, as the governing equations show strongly nonlinear characteristics. This study proposes a thermo-hydro-mechanical coupled model considering the heat transfer, waterinfiltration, and deformation of partially saturated soil to reproduce the freezing process of partially saturatedfrost susceptible soils distributed in Hokkaido. This model better considers the water-ice phase change and soilfreezing characteristic curve (SFCC) during freezing under field conditions. The results from the multiphysicssimulations agree well with the frost heave and water migration data from frost heave tests of Touryo soil andFujinomori soil. In addition, this study discussed the influence of the various factors on frost heave amount,including temperature gradients, overburden pressures, water supply conditions, cooling rates, and initial saturation. The simulation results indicate that the frost heave ratio is proportional to the initial degree of saturationand is inversely proportional to the cooling rate and overburden pressure.Moreover, simulation under the open system generates much more frost heave than under the closed system.Finally, the main features of the proposed model are revealed by simulating a closed-system frost heave test. Thesimulation results indicate that the proposed model adequately captures the coupling characteristics of water andice redistribution, temperature development, hydraulic conductivity, and suction in the freezing process. Togetherwith the decreased hydraulic conductivity, the increased suction controls the water flow in the freezing zone. Theinflow water driven by cryogenic suction gradient feeds the ice formation, leads to a rapid increase in total watercontent, expanding the voids that exceed the initial porosity and contributing to the frost heave.
基金Project(11072269) supported by the National Natural Science Foundation of ChinaProject(20090162110066) supported by the Research Fund for the Doctoral Program of Higher Education of China
文摘Based on fluid mechanics, thermodynamics and damage mechanics, thermal-hydro-mechanical(THM) coupling damage model of brittle rock is established by analyzing THM coupling mechanism, where THM coupling damage variable DTHM is dominated by TH coupling damage variable DTH, TM coupling damage variable DTM and HM coupling damage variable DHM, and DTH is firstly expressed in term of dimensionless total thermal conductivity of the water Nu. Permeability test, uni-axial compression test and THM coupling test are conducted to measure the permeability, elastic modulus and THM coupling stress-strain curves of brittle rock. The tested values of THM coupling elastic modulus THME? are in good agreement with the predicted values of THM coupling elastic modulus ETHM, which can verify the newly established THM coupling damage model.
文摘One of the most suitable ways under study for the disposal of high-level radioactive waste (HLW) is isolation in deep geological repositories. It is very important to research the thermo-hydro- mechanical (THM) coupled processes associated with an HLW disposal repository. Non-linear coupled equations, which are used to describe the THM coupled process and are suited to saturated-unsaturated porous media, are presented in this paper. A numerical method to solve these equations is put forward, and a finite element code is developed. This code is suited to the plane strain or axis-symmetry problem. Then this code is used to simulate the THM coupled process in the near field of an ideal disposal repository. The temperature vs. time, hydraulic head vs. time and stress vs. time results show that, in this assumed condition, the impact of temperature is very long (over 10 000 a) and the impact of the water head is short (about 90 d). Since the stress is induced by temperature and hydraulic head in this condition, the impact time of stress is the same as that of temperature. The results show that THM coupled processes are very important in the safety analysis of an HLW deep geological disposal repository.
基金Supported by the National Natural Science Foundation of China (50579087,50720135906, 50539050)CAS/SAFEA International Partnership Program for Creative Research Teams
文摘In examining potential host rocks for such purposes as the disposal of high-level radioactive wastes,it is important to understand the coupled thermo-hydro-mechanical（THM） behavior of a porous medium.A rigorous and fully unified coupled thermo-hydro-mechanical model for unsaturated porous media is required to simulate the complex coupling mechanisms involved.Based on modified Darcy＇s and Fourier＇s laws,equations of mechanical equilibrium,mass conservation and energy conservation are derived by introducing void ratio and volumetric liquid water content into the model.The newly derived model takes into account the effects of temperature on the dynamic viscosity of liquid water and void ratio,the influence of liquid flow on temperature gradient（thermo-osmosis）,the influence on mass and heat conservation equations,and the influence of heat flow on water pressure gradient and thermal convection.The new coupled THM constitutive model is constructed by a finite element program and is used to simulate the coupled behavior of a tunnel during excavation,ventilation and concrete lining stages.Oil and gas engineering,underground disposal of nuclear waste and tunnel engineering may be benefited from the development of the new model.
基金funded by the National Natural Science Foundation of China(Grant No.NSFC51374147)the German Society for Petroleum and Coal Science and Technology(Grant No.DGMK680-4)
文摘As one of the most important ways to reduce the greenhouse gas emission,carbon dioxide（CO2）enhanced gas recovery（CO2-EGR） is attractive since the gas recovery can be enhanced simultaneously with CO2sequestration.Based on the existing equation of state（EOS） module of TOUGH2 MP,extEOS7C is developed to calculate the phase partition of H2O-CO2-CH4-NaCl mixtures accurately with consideration of dissolved NaCI and brine properties at high pressure and temperature conditions.Verifications show that it can be applied up to the pressure of 100 MPa and temperature of 150℃.The module was implemented in the linked simulator TOUGH2MP-FLAC3 D for the coupled hydro-mechanical simulations.A simplified three-dimensional（3D）1/4 model（2.2 km×1 km×1 km） which consists of the whole reservoir,caprock and baserock was generated based on the geological conditions of a gas field in the North German Basin.The simulation results show that,under an injection rate of 200,000 t/yr and production rate of 200,000 sm3/d,CO2breakthrough occurred in the case with the initial reservoir pressure of 5 MPa but did not occur in the case of 42 MPa.Under low pressure conditions,the pressure driven horizontal transport is the dominant process;while under high pressure conditions,the density driven vertical flow is dominant.Under the considered conditions,the CO2-EGR caused only small pressure changes.The largest pore pressure increase（2 MPa） and uplift（7 mm） occurred at the caprock bottom induced by only CO2injection.The caprock had still the primary stress state and its integrity was not affected.The formation water salinity and temperature variations of ±20℃ had small influences on the CO2-EGR process.In order to slow down the breakthrough,it is suggested that CO2-EGR should be carried out before the reservoir pressure drops below the critical pressure of CO2.
基金financial support received from the“Iran’s Ministry of Science Research and Technology”(PhD students’sabbatical grants)funding from the European Research Council under the European Union’s Horizon 2020 Research and Innovation Program through the Starting Grant GEoREST(www.georest.eu)(Grant Agreement No.801809)+1 种基金support by the Korea-EU Joint Research Program of the National Research Foundation of Korea through Grant No.NRF2015K1A3A7A03074226funded by the Korean Government’s Ministry of Science and Information and Communication Technology(ICT)in the framework of the European Union’s Horizon 2020 Research and Innovation Program(Grant No.691728)。
文摘Geo-energy and geo-engineering applications,such as improved oil recovery(IOR),geologic carbon storage,and enhanced geothermal systems(EGSs),involve coupled thermo-hydro-mechanical(THM)processes that result from fluid injection and production.In some cases,reservoirs are highly fractured and the geomechanical response is controlled by fractures.Therefore,fractures should explicitly be included into numerical models to realistically simulate the THM responses of the subsurface.In this study,we perform coupled THM numerical simulations of water injection into naturally fractured reservoirs(NFRs)using CODE_BRIGHT and TOUGH-UDEC codes.CODE_BRIGHT is a finite element method(FEM)code that performs fully coupled THM analysis in geological media and TOUGH-UDEC sequentially solves coupled THM processes by combining a finite volume method(FVM)code that solves nonisothermal multiphase flow(TOUGH2)with a distinct element method(DEM)code that solves the mechanical problem(UDEC).First,we validate the two codes against a semi-analytical solution for water injection into a single deformable fracture considering variable permeability based on the cubic law.Then,we compare simulation results of the two codes in an idealized conceptual model that includes one horizontal fracture and in a more realistic model with multiple fractures.Each code models fractures differently.UDEC calculates fracture deformation from the fracture normal and shear stiffnesses,while CODE_BRIGHT treats fractures as equivalent porous media and uses the equivalent Young’s modulus and Poisson’s ratio of the fracture.Finally,we obtain comparable results of pressure,temperature,stress and displacement distributions and evolutions for the single horizontal fracture model.Despite some similarities,the two codes provide increasingly different results as model complexity increases.These differences highlight the challenging task of accurately modeling coupled THM processes in fractured media given their high nonlinearity.
基金financed and supported by the German research institute "Gesellschaft für Anlagen-und Reaktorsicherheit (GRS) mbH"
文摘The main objective of this paper is to investigate and analyse the thermo-hydro-mechanical（THM） coupling phenomena and their influences on the repository safety.In this paper,the high-level waste（HLW） disposal concept in drifts in clay formation with backfilled bentonite buffer is represented numerically using the CODE BRIGHT developed by the Technical University of Catalonia in Barcelona.The parameters of clay and bentonite used in the simulation are determined by laboratory and in situ experiments.The calculation results are presented to show the hydro-mechanical（HM） processes during the operation phase and the THM processes in the after-closure phase.According to the simulation results,the most probable critical processes for the disposal project have been represented and analyzed.The work also provides an input for additional development regarding the design,assessment and validation of the HLW disposal concept.
