• Tunnel and Underground Space
• /
• v.30 no.3
• /
• pp.256-269
• /
• 2020

The objective of this study is to evaluate the stress thresholds in crack development and the corrected fracture toughness of KURT granite under dry and saturated conditions. The stress thresholds were identified by calculation of inelastic volumetric strain from an uniaxial compression test. The corrected fracture toughness was estimated by using the Level II method (Chevron Bend specimen), suggested by ISRM (1988), in which non-linear behaviors of rock was taken into account. Average crack initiation stress(σ_ci) and crack damage stress(σ_cd) under a dry condition were 91.1 MPa and 128.7 MPa. While, average crack initiation stress(σ_ci) and crack damage stress(σ_cd) under a saturated condition were 58.2 MPa and 68.2 MPa. The crack initiation stress and crack damage stress of saturated ones decreased 36% and 47% respectively compared to those of dry specimens. A decrease in crack damage stress is relatively larger than that of crack initiation stress under a saturated condition. This indicates that the unstable crack growth can be more easily generated because of the saturation effect of water compared to the dry condition. The average corrected fracture toughness of KURT granite was 0.811 MPa·m^0.5. While, the fracture toughness of saturated KURT granite(K_CB) was 0.620 MPa·m^0.5. The corrected fracture toughness of rock in saturated condition decreases by 23.5% compared to that in dry condition. It is found that the resistance to crack propagation decreases under the saturated geological condition.

• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.382-393
• /
• 2020

The engineered barrier system of high-level radioactive waste disposal must maintain its performance in the long term, because it must play a role in slowing the rate of leakage to the surrounding rock mass even if a radionuclide leak occurs from the canister. In particular, it is very important to clarify gas dilation flow phenomenon clearly, that occurs only in a medium containing a large amount of clay material such as a bentonite buffer, which can affect the long-term performance of the bentonite buffer. Accordingly, DECOVALEX-2019 Task A was conducted to identify the hydraulic-mechanical mechanism for the dilation flow, and to develop and verify a new numerical analysis technique for quantitative evaluation of gas migration phenomena. In this study, based on the conventional two-phase flow and mechanical behavior with effective stresses in the porous medium, the hydraulic-mechanical model was developed considering the concept of damage to simulate the formation of micro-cracks and expansion of the medium and the corresponding change in the hydraulic properties. Model verification and validation were conducted through comparison with the results of 1D and 3D gas injection tests. As a result of the numerical analysis, it was possible to model the sudden increase in pore water pressure, stress, gas inflow and outflow rate due to the dilation flow induced by gas pressure, however, the influence of the hydraulic-mechanical interaction was underestimated. Nevertheless, this study can provide a preliminary model for the dilation flow and a basis for developing an advanced model. It is believed that it can be used not only for analyzing data from laboratory and field tests, but also for long-term performance evaluation of the high-level radioactive waste disposal system.

• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.320-334
• /
• 2020

This study presents an application of hydro-mechanical coupled Particle Flow Code 3D (PFC3D) to simulation of fluid injection induced fault slip experiment conducted in Mont Terri Switzerland as a part of a task in an international research project DECOVALEX-2019. We also aimed as identifying the current limitations of the modelling method and issues for further development. A fluid flow algorithm was developed and implemented in a 3D pore-pipe network model in a 3D bonded particle assembly using PFC3D v5, and was applied to Mont Terri Step 2 minor fault activation experiment. The simulated results showed that the injected fluid migrates through the permeable fault zone and induces fault deformation, demonstrating a full hydro-mechanical coupled behavior. The simulated results were, however, partially matching with the field measurement. The simulated pressure build-up at the monitoring location showed linear and progressive increase, whereas the field measurement showed an abrupt increase associated with the fault slip We conclude that such difference between the modelling and the field test is due to the structure of the fault in the model which was represented as a combination of damage zone and core fractures. The modelled fault is likely larger in size than the real fault in Mont Terri site. Therefore, the modelled fault allows several path ways of fluid flow from the injection location to the pressure monitoring location, leading to smooth pressure build-up at the monitoring location while the injection pressure increases, and an early start of pressure decay even before the injection pressure reaches the maximum. We also conclude that the clay filling in the real fault could have acted as a fluid barrier which may have resulted in formation of fluid over-pressurization locally in the fault. Unlike the pressure result, the simulated fault deformations were matching with the field measurements. A better way of modelling a heterogeneous clay-filled fault structure with a narrow zone should be studied further to improve the applicability of the modelling method to fluid injection induced fault activation.

• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.359-381
• /
• 2020

Within the framework of DECOVALEX-2019 Task D, full-scale engineered barriers experiment (FEBEX) at Grimsel Test Site was numerically simulated to investigate an applicability of implemented Barcelona basic model (BBM) into TOUGH2-MP/FLAC3D simulator, which was developed for the prediction of the coupled thermo-hydro-mechanical behavior of bentonite buffer. And the calculated heater power, temperature, relative humidity, total stress, saturation, water content and dry density were compared with in situ data monitored in the various sections. In general, the calculated heater power and temperature provided a fairly good agreement with experimental observations, however, the difference between power of heater #1 and that of heater #2 could not captured in the numerical analysis. It is necessary to consider lamprophyre with low thermal conductivity around heater #1 and non-simplified installation progresses of bentonite blocks in the tunnel for better modeling results. The evolutions and distributions of relative humidity were well reproduced, but hydraulic model needs to be modified because the re-saturation process was relatively fast near the heaters. In case of stress evolutions due to the thermal and hydraulic expansions, the computed stress was in good agreement with the data. But, the stress is slightly higher than the measured in situ data at the early stage of the operation, because gap between rock mass and bentonite blocks have not been considered in the numerical simulations. The calculated distribution of saturation, water content, and dry density along the radial distance showed good agreement with the observations after the first and final dismantling. The calculated dry density near the center of the FEBEX tunnel and heaters were overestimated compared with the observations. As a result, the saturation and water content were underestimated with the measurements. Therefore, numerical model of permeability is needed to modify for the production of better numerical results. It will be possible to produce the better analysis results and more realistically predict the coupled THM behavior in the bentonite blocks by performing the additional studies and modifying the numerical model based on the results of this study.

• Economic and Environmental Geology
• /
• v.51 no.1
• /
• pp.67-76
• /
• 2018

This study reviews three eco-friendly energy towns with hybrid thermal energy supply systems and borehole thermal energy storage (BTES) in Canada and Denmark. The district heating and cooling systems were designed by using multi-source energy for the higher efficiency and reliability as well as environment. ADEU (Alexandra District Energy Utility) located at the developing area in the city of Richmond, Canada was designed to supply district energy with the installation of 726 borehole heat exchangers (BHEs) and a backup boiler using natural gas. DLSC (Drake Landing Solar Community) located in the town of Okotoks, Canada is a district system to store solar thermal energy underground during the summer season by seasonal BTES with 144 BHEs. Brædstrup Solpark district heating system located in Denmark has been conducted energy supply from multiple energy sources of solar thermal, heat pump, boiler plants and seasonal BTES with 48 BHEs. These systems are designed based on social and economic benefits as well as nature-friendly living space according to the city based energy perspective. Each system has the energy center which distribute the stored thermal energy to each house for heating during the winter season. The BHE depth and ground thermal storage volume are designed by the heating and cooling load as well as the condition of ground water flow and thermophysical properties of the ground. These systems have been proved the reliance and economic benefits by providing consistent energy supply with competitive energy price for many years. In addition, the several expansions of the service area in ADEU and Brædstrup Solpark have been processed based on energy supply master plan. In order to implement this kind of project in our country, the regulation and policy support of government or related federal organization are required. As well as the government have to make a energy management agency associated with long-term supply energy plan.

• Tunnel and Underground Space
• /
• v.22 no.4
• /
• pp.266-275
• /
• 2012

In this study, a method is devised to implement a supercritical $CO_2$ ($scCO_2$) injection environment on a laboratory scale and to investigate the effects of $scCO_2$ on the properties of rock specimens. Specimens of shale and sandstone normally constituting the cap rock and reservoir rock, respectively, were kept in a laboratory reactor chamber with $scCO_2$ for two weeks. From this stage, a chemical reaction between rock surface and the $scCO_2$ was induced. The effect of saline water was also investigated by comparing three conditions ($scCO_2$-rock, $scCO_2-H_2O$-rock and $scCO_2$-brine(1M)-rock). Finally, we checked the changes in the properties before and after the reaction by destructive and nondestructive testing procedures. The swelling of shale was a main concern in this case. The experimental results suggested that $scCO_2$ has a greater effect on the swelling of the shale than pure water and brine. It was also observed that the largest swelling displacement of shale occurred after a reaction with the $H_2O-scCO_2$ solution. The results of a series of the destructive and nondestructive tests indicate that although each of the property changes of the rock differed depending on the reaction conditions, the $H_2O-scCO_2$ solution had the greatest effect. In this study, shale was highly sensitive to the reaction conditions. These results provide fundamental information pertaining to the stability of $CO_2$ storage sites due to physical and chemical reactions between the rocks in these sites and $scCO_2$.

• Tunnel and Underground Space
• /
• v.21 no.1
• /
• pp.33-40
• /
• 2011

This study focuses on having a temperature and water pressure effects on the change of material properties of rocks. Granite and limestone specimens from Gagok Mine were thermally treated with predetermined temperatures of 200, 300, 400, 500, 600 and $700^{\circ}C$ (excepting $700^{\circ}C$ for limestone) to estimate the reduction of material properties of rocks caused by heat. Specific gravity, effective porosity, elastic wave velocity, uniaxial compressive strength, Young's modulus and Poisson's ratio for pre-heated specimens were measured. With increasing temperature, material properties of both rock specimens change sequentially. Significant changes of specific gravity, effective porosity and elastic wave porosity occur above $400^{\circ}C$ for granite and $300^{\circ}C$ for limestone. Changes of uniaxial compressive strength, Young's modulus and Poisson's ratio seem to be similar to those of physical properties. GSI of 500, 600 and $700^{\circ}C$ specimens inferred by using uniaxial compressive strength and Young's modulus of preheated granite specimens is found to be 81, 66 and 58 each. In case of pre-heated limestone specimens of 400, 500 and $600^{\circ}C$, the corresponding GSI is 76, 71 and 65 each. 500, 600 and $700^{\circ}C$ granite specimens and 400, 500 and $600^{\circ}C$ limestone specimens were pressurized to 7.5 MPa and their effective porosity, elastic wave velocity, uniaxial compressive strength and Young's modulus were measured. The average value of material properties (mentioned above) of 500, 600 and $700^{\circ}C$ granite specimens under water pressure compared with material properties of non-pressurized pre-heated specimens exhibits the reduction of 7.6, 11.3 and 14.9%, respectively. In case of 400, 500 and $600^{\circ}C$ limestone specimens under water pressure, the average value of material properties decreases by 8.2, 13.8 and 21.9%, respectively.

• Korean Journal of Heritage: History & Science
• /
• v.47 no.1
• /
• pp.102-115
• /
• 2014

The stamped earth method is a typical ancient engineering technique which consists of in-filling wooden frame with layers of stamped earth or sand. This method has been universally used to construct earthen walls and buildings, etc. The purpose of this article is to understand the construction method and principles of the stamped earthen wall through analysis of various construction techniques of Pungnaptoseong Fortress(Earthen Fortification in Pungnap-dong). First of all, the ground was leveled and the foundations for the construction of the earthen wall were laid. The underground foundation of the earthen walls was usually constructed by digging into the ground and then in-filling this space with layers of mud clay. Occasionally wooden posts or paving stones which may have been used to reinforce the soft ground were driven in. The method of adding layers of stamped earth at an oblique angle to either side of a central wall is the most characteristic feature of Pungnaptoseong Fortress. Even though the traces of fixing posts, boards, and the hardening of earth - all signatures of the stamped earth technique - have not been identified, evidence of a wooden frame has been found. It has also been observed that this section was constructed by including layers of mud clay and organic remains such as leaves and twigs in order to strengthen the adhesiveness of the structures. The outer part of the central wall was constructed by the anti-slope stamped earth technique to protect central wall. In addition a final layer of paved stones was added to the upper part of the wall. These stone layers and the stone wall were constructed in order to prevent the loss of the earthen wall and to discharge and drain water. Meanwhile, the technique of cementing with fire was used to control damp and remove water in stamped earth. It can not be said at present that the stamped earth method has been confirmed as the typical construction method of Korean ancient earthen walls. If we make a comparative study of the evidence of the stamped earth technique at Pungnaptoseong Fortress with other archeological sites, progress will be made in the investigation of the construction method and principles of stamped earthen wall.

• Tunnel and Underground Space
• /
• v.28 no.6
• /
• pp.670-691
• /
• 2018

This study presents the research results of the BMT(Benchmark Model Test) simulations of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to predict fault reactivation and the coupled hydro-mechanical behavior of fault. BMT scenario simulations of Task B were conducted to improve each numerical model of participating group by demonstrating the feasibility of reproducing the fault behavior induced by water injection. The BMT simulations consist of seven different conditions depending on injection pressure, fault properties and the hydro-mechanical coupling relations. TOUGH-FLAC simulator was used to reproduce the coupled hydro-mechanical process of fault slip. A coupling module to update the changes in hydrological properties and geometric features of the numerical mesh in the present study. We made modifications to the numerical model developed in Task B Step 1 to consider the changes in compressibility, Permeability and geometric features with hydraulic aperture of fault due to mechanical deformation. The effects of the storativity and transmissivity of the fault on the hydro-mechanical behavior such as the pressure distribution, injection rate, displacement and stress of the fault were examined, and the results of the previous step 1 simulation were updated using the modified numerical model. The simulation results indicate that the developed model can provide a reasonable prediction of the hydro-mechanical behavior related to fault reactivation. The numerical model will be enhanced by continuing interaction and collaboration with other research teams of DECOVALEX-2019 Task B and validated using the field experiment data in a further study.

• Tunnel and Underground Space
• /
• v.31 no.4
• /
• pp.289-308
• /
• 2021

This study conducts coupled thermo-hydro-mechanical numerical modeling to investigate the maximum temperature and conditions for securing mechanical stability of the high-level radioactive waste repository when temperature criteria of bentonite buffer are 100℃ and 125℃, respectively. In case of temperature criterion of buffer as 100℃, the maximum temperatures at the interface between canister and buffer are calculated to be 99.4℃ and 99.8℃, respectively for a case with disposal tunnel spacing of 40 m and deposition hole spacing of 5.5 m and for the other case with disposal tunnel spacing of 30 m and deposition hole spacing of 6.5 m. In case of temperature criterion of buffer as 125℃, spacings of disposal tunnel and deposition hole could be decreased to 30 m and 4.5 m, respectively, which reduces the disposal area up to 55% compared to the disposal area of KRS⁺. According to analysis of mechanical stability for various disposal spacings, RMR of rock mass for KRS⁺ should be larger than 72.4 which belongs to good rock in RMR classification to prevent failure of rock mass. As disposal spacing is decreased, required RMR of rock mass is increased. In order to prevent failure of rock mass for a case with disposal tunnel spacing of 30 m and deposition hole spacing of 4.5 m, RMR larger than 87.3 is needed. However, mechanical stability of the repository is secured for all cases with RMR over 75 considering the enhancement of rock strength due to confining stress induced by swelling of the bentonite buffer and backfill.