• 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.

• Tunnel and Underground Space
• /
• v.32 no.1
• /
• pp.30-58
• /
• 2022

The deep geological repository for high-level radioactive waste disposal is a multi barrier system comprised of engineered barriers and a natural barrier. The long-term integrity of the deep geological repository is affected by the coupled interactions between the individual barrier components. Erosion and piping phenomena in the compacted bentonite buffer due to buffer-rock interactions results in the removal of bentonite particles via groundwater flow and can negatively impact the integrity and performance of the buffer. Rapid groundwater inflow at the early stages of disposal can lead to piping in the bentonite buffer due to the buildup of pore water pressure. The physiochemical processes between the bentonite buffer and groundwater lead to bentonite swelling and gelation, resulting in bentonite erosion from the buffer surface. Hence, the evaluation of erosion and piping occurrence and its effects on the integrity of the bentonite buffer is crucial in determining the long-term integrity of the deep geological repository. Previous studies on bentonite erosion and piping failed to consider the complex coupled thermo-hydro-mechanical-chemical behavior of bentonite-groundwater interactions and lacked a comprehensive model that can consider the complex phenomena observed from the experimental tests. In this technical note, previous studies on the mechanisms, lab-scale experiments and numerical modeling of bentonite buffer erosion and piping are introduced, and the future expected challenges in the investigation of bentonite buffer erosion and piping are summarized.

• Economic and Environmental Geology
• /
• v.55 no.1
• /
• pp.63-76
• /
• 2022

Dacitic tuffs, 97 to 118 m thick, were recovered from the lower part of the subsurface Seongdongri Formation, Janggi Basin, which was drilled to assess the potential for underground storage of carbon dioxide. The tuffs are divided into four depositional units(Unit 1 to 4) based on internal structures and particle componentry. Unit 1 and Units 3/4 are ignimbrites that accumulated in subaerial and subaqueous settings, respectively, whereas Unit 2 is braided-stream deposits that accumulated during a volcanic quiescence, and no dacitic tuff is observed. A series of analysis shows that mordenite and clinoptilolite mainly fill the vesicles of glass shards, suggesting their formation by replacement and dissolution of volcanic glass and precipitation from interstitial water during burial and diagenesis. Glass-replaced clinoptilolite has higher Si/Al ratios and Na contents than the vesicle-filling clinoptilolite in Units 3. However, the composition of clinoptilolite becomes identical in Unit 4, irrespective of the occurrence and location. This suggests that the Si/Al ratio and pH in the interstitial water increased with time because of the replacement and leaching of volcanic glass, and that the composition of interstitial water was different between the eastern and western parts of the basin during the formation of the clinoptilolite in Units 1 and 3. It is also inferred that the formation of the two zeolite minerals was sequential according to the depositional units, i.e., the clinoptilolite formed after the growth of mordenite. To summarize, during a volcanic quiescence after the deposition of Unit 1, pH was higher in the western part of the basin because of eastward tilting of the basin floor, and the zeolite ceased to grow because of the closure of the pore space as a result of the growth of smectite. On the other hand, clinoptilolite could grow in the eastern part of the basin in an open system affected by groundwater, where braided stream was developed. Afterwards, Units 3 and 4 were submerged under water because of the basin subsidence, and the alkali content of the interstitial water increased gradually, eventually becoming identical in the eastern and western parts of the basin. This study thus shows that volcanic deposits of similar composition can have variable distribution of zeolite mineral depending on the drainage and depositional environment of basins.