• Tunnel and Underground Space
• /
• v.31 no.5
• /
• pp.309-332
• /
• 2021

Natural barrier systems surrounding the geological repository for the high-level radioactive waste should guarantee the hydraulic performance for preventing or delaying the leakage of radionuclide. In the case of the behavior of a crystalline rock, the hydraulic performance tends to be decided by the existence of discontinuities, so the coupled hydro-mechanical(HM) processes on the discontinuities should be characterized. The discontinuum modelling can describe the complicated behavior of discontinuities including creation, propagation, deformation and slip, so it is appropriate to model the behavior of a crystalline rock. This paper investigated the coupled HM processes in discontinuum modelling such as UDEC, 3DEC, PFC, DDA, FRACOD and TOUGH-UDEC. Block-based discontinuum methods tend to describe the HM processes based on the fluid flow through the discontinuities, and some methods are combined with another numerical tool specialized in hydraulic analysis. Particle-based discontinuum modelling describes the overall HM processes based on the fluid flow among the particles. The discontinuum methods that are currently available have limitations: exclusive simulations for two-dimension, low hydraulic simulation efficiency, fracture-dominated fluid flow and simplified hydraulic analysis, so it could be improper to the modelling the geological repository. Based on the concepts of various discontinuum modelling compiled in this paper, the advanced numerical tools for describing the accurate coupled HM processes of the deep geological repository should be developed.

• Tunnel and Underground Space
• /
• v.17 no.5
• /
• pp.411-427
• /
• 2007

Reliability in tunnel analysis is necessary to accomplish technically sound design and economical construction. For this, a thorough understanding of the construction procedure including the ground-support interaction has to be obtained. This paper describes a proper modelling technique to simulate the behavior of the steel fiber reinforced shotcrete (SFRS) which maintain the supporting capability in post-failure regime. The additional supporting effect of the steel support was also verified by 3-D analyses and a new load distribution factor were proposed. The use of the plastic moment limit (PML) alone can eliminate the occurrence of the awkwardly high tensile stress in the shotcrete and can successfully model the post-peak ductile behavior of the SFRS. But with this method, moment is limited whenever the stress caused by moment reaches tensile strength of the shotcrete irrespective of the stress by axial force. Therefore, it was necessary to find a more comprehensive method which can reflect the influence of the moment and axial force. This can be accomplished by the proper use of "liner element" which is the built-in model in FLAC. In this model, the peak and residual strength as well as the uniaxial compressive strength of the SFRS can be specified. Analyses were conducted with these two models on the 2-lane road tunnels excavated in class IV and V rock mass and results were compared with the conventional elastic beam model. Results showed that both models can reflect the fracture toughness of the SFRS which could not be accomplished by the elastic beam model.

• Tunnel and Underground Space
• /
• v.18 no.5
• /
• pp.343-354
• /
• 2008

A new method is suggested herein to measure the virgin earth stresses by means of a borehole. This novel concept is basically a combination of borehole stress relieving and borehole fracturing techniques. The destressing of the borehole is achieved by means of inducing thermal tensile stresses at the borehole periphery by using a cryogenic fluid such as Liquid Nitrogen($LN_2$). The borehole wall eventually develops fractures when the induced thermal stresses exceed the existing compressive stresses at the borehole periphery in addition to the tensile strength of the rock. The above concept is theoretically analyzed for its potential applicability to interpret in situ stress levels from the tensile fracture stresses and the corresponding borehole wall temperatures. Coupled thermo-mechanical numerical simulations are also conducted using FLAC3D, with thermal option, to check the validity of the proposed techniques. From the preliminary theoretical and numerical analysis, the method suggested for the measurement of in situ stresses appears to be capable of accurate estimation of the virgin stresses by monitoring tensile crack formation at a borehole wall and recording the wall temperatures at the time of crack initiation.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.19 no.4 s.74
• /
• pp.429-440
• /
• 2006

This paper comprises two specific objectives. The first is to examine the applicability of ordinary kriging interpolation(OK) to the p-adaptivity of the finite element method that is based on variogram modeling. The second objective Is to present the adaptive procedure by the hierarchical p-refinement in conjunction with a posteriori error estimator using the modified S.P.R. (superconvergent patch recovery) method. The ordinary kriging method that is one of weighted interpolation techniques is applied to obtain the estimated exact solution from the stress data at the Gauss points. The weight factor is determined by experimental and theoretical variograms for interpolation of stress data apart from the conventional interpolation methods that use an equal weight factor. In the p-refinement, the analytical domain has to be refined automatically to obtain an acceptable level of accuracy by increasing the p-level non-uniformly or selectively. To verify the performance of the modified S.P.R. method, the new error estimator based on limit value has been proposed. The validity of the proposed approach has been tested with the help of some benchmark problems of linear elastic fracture mechanics such as a centrally cracked panel, a single edged crack, and a double edged crack.

• Tunnel and Underground Space
• /
• v.34 no.4
• /
• pp.330-354
• /
• 2024

In site characterization of high-level radioactive waste, discontinuity(joint) distribution and rock mass classification, which are key evaluation parameters in the rock engineering field, were evaluated using deep boreholes in the Wonju granite and Chuncheon granite, which belong to Mesozoic Jurassic era. To evaluate joint distribution characteristics, fracture zones and joint surfaces extracted from ATV data were used, and major joint sets were evaluated along with joint frequency according to depth, dip direction, and dip. Both the Wonju and Chuncheon granites that were studied showed a tendency for the frequency of joints to increase linearly with depth, and joints with high angles were relatively widely distributed. In addition, relatively large amounts of weathering tended to occur even in deep depth due to groundwater inflow through high-angle joints. RQD values remained consistently low even at considerable depth. Meanwhile, joint groups with low angles showed different joint characteristics from joint sets with high angles. Rock mass classification was performed based on RMR system, and along with rock mass classification for 50 m intervals where uniaxial compressive strength was performed, continuous rock mass classification according to depth was performed using velocity log data and geostatistical techniques. The Wonju granite exhibited a superior rock mass class compared to the Chuncheon granite. In the 50 m interval and continuous rock mass classification, the shallow part of the Wonju granite showed a higher class than the deep part, and the deep part of the Chuncheon granite showed a higher class than the shallow part.

• Tunnel and Underground Space
• /
• v.34 no.4
• /
• pp.374-392
• /
• 2024

In order to successfully select a site for deep geological disposal of high-level radioactive waste, it is important to perform the stepwise approach along with the systematic selection and survey of evaluation parameters of geological environmental characteristics suitable for the domestic geological environment. In this study, we evaluated the characteristics of hydraulic conductivity, which is considered the most important evaluation parameter in the field of hydrogeology, targeting a deep-depth rock aquifer where actual disposal facilities are expected to be located. In particular, for the first time in Korea, we obtained in-situ pressure-flow data by directly conducting hydraulic tests in boreholes at depths ranging from 500 m to 750 m in various rock types distributed in Korea (granite/volcanic rock/gneiss/mudstone). And we derived hydraulic conductivity values by rock types and depth using verified analytical methods. For this purpose, precision hydraulic testing equipment developed in-house through this study was used, and detailed investigation procedures based on standard test methods were applied to field tests. As a result of the analysis, the average hydraulic conductivity value was found to be in the range of 10^-9 m/s in all granite/volcanic rock/gneiss areas. In the mudstone area, an average hydraulic conductivity value of 10^-11 m/s was derived, which was about 100 times (2 orders of magnitude) lower than that of the fractured rock aquifers. Moreover, permeability tended to slightly decrease with depth in fractured rock aquifers (granite and volcanic rock areas) containing many rock fractures. The gneiss area tended to have large local differences in permeability according to the composition of the stratum and the development of fracture zones rather than depth. In mudstone areas with weak fracture development, there was no significant variation in rock permeability according to depth. The hydraulic conductivity results by various rock types and depth presented in this study are expected to be utilized in building a foundational database for the site selection, design, and construction of disposal facilities in Korea.

• Tunnel and Underground Space
• /
• v.13 no.1
• /
• pp.52-63
• /
• 2003

Discontinuities such as faults, fractures and joints in rock mass play the dominant role in the mechanical and hydraulic properties of the rock mass. The key factors that influence on the flow of groundwater are hydraulic and geometric characteristics of discontinuities and their connectivity. In this study, a program that analyzes groundwater flow in the 3D discontinuity network was developed on the assumption that the discontinuity characteristics such as density, trace length, orientation and aperture have particular distribution functions. This program generates discontinuities in a three-dimensional space and analyzes their connectivity and groundwater flow. Due to the limited computing capacity In this study, REV was not exactly determined, but it was inferred to be greater than 25$\times$25$\times$25 ㎥. By calculating the extent of aperture that influences on the groundwater flow, it was found that the discontinuities with the aperture smaller than 30% of the mean aperture had little influence on the groundwater flow. In addition, there was little difference in the equivalent hydraulic conductivity for the the two cases when considering and not considering the boundary effect. It was because the groundwater flow was mostly influenced by the discontinuities with large aperture. Among the parameters considered in this study, the length, aperture, and orientation of discontinuities had the greatest influence on the equivalent hydraulic conductivity of rock mass in their order. In case of existence of a fault in rock mass, elements of the equivalent hydraulic conductivity tensor parallel to the fault fairly increased in their magnitude but those perpendicular to the fault were increased in a very small amount at the first stage and then converged.