• Tunnel and Underground Space
• /
• v.23 no.5
• /
• pp.383-391
• /
• 2013

A number of true triaxial tests on rock samples have been conducted since the late 1960 and their results strongly suggest that the intermediate principal stress has a considerable effect on rock strength. Based on these experimental evidence, various 3-D rock failure criteria accounting for the effect of the intermediate principal stress have been proposed. Most of the 3-D failure criteria, however, are focused on the phenomenological description of the rock strength from the true triaxial tests, so that the associated strength parameters have little physical meaning. In order to confirm the likelihood that the intermediate principal stress dependency of rock strength is related to the presence of weak planes and their distribution to the preferred orientation, true triaxial tests are simulated with the transversely isotropic rock model. The conventional Mohr-Coulomb criterion is extended to its anisotropic version by incorporating the concept of microstructure tensor. With the anisotropic Mohr-Coulomb criterion, the critical plane approach is applied to calculate the strength of the transversely isotropic rock model and the orientation of the fracture plane. This investigation hints that the spatial distribution of microstructural planes with respect to the principal stress triad is closely related to the intermediate principal stress dependency of rock strength.

• Tunnel and Underground Space
• /
• v.26 no.2
• /
• pp.59-67
• /
• 2016

In this paper, we perform a numerical analysis to evaluate the potential of fault reactivation during gas production from hydrate bearing sediments and the moment magnitude of induced seismicity. For the numerical analysis, sequential coupling of TOUGH+Hydrate and FLAC3D was used and the change in effective stress and consequent geomechanical deformation including fault reactivation was simulated by assuming that Mohr-Coulomb shear resistance criterion is valid. From the test production simulation of 30 days, we showed that pore pressure reduction as well as effective stress change hardly induces the fault reactivation in the vicinity of a production well. We also investigated the influence of stress state conditions to a fault reactivation, and showed that normal fault stress regime, where vertical stress is relatively greater than horizontal, may have the largest potential for the reactivation. We tested one simulation that earthquake can be induced during gas production and calculated the moment magnitude of the seismicity. Our calculation presented that all the magnitudes from the calculation were negative values, which indicates that induced earthquakes can be grouped into micro-seismic and as small as hardly perceived by human beings. However, it should be noted that the current simulation was carried out using the highly simplified geometric model and assumptions such that the further simulations for a scheduled test production and commercial scale production considering complex geometric conditions may produce different results.

• Tunnel and Underground Space
• /
• v.31 no.4
• /
• pp.221-242
• /
• 2021

A numerical study of the performance assesment of coupled thermo-hydro-mechanical (THM) processes in improved Korean reference disposal system (KRS⁺) for high-level radioactive waste is conducted using TOUGH2-MP/FLAC3D simulator. Decay heat from high-level radioactive waste increases the temperature of the repository, and it decreases as decay heat is reduced. The maximum temperature of the repository is below a maximum temperature criterion of 100℃. Saturation of bentonite buffer adjacent to the canister is initially reduced due to pore water evaporation induced by temperature increase. Bentonite buffer is saturated 250 years after the disposal of high-level radioactive waste by inflow of groundwater from the surrounding rock mass. Initial saturation of rock mass decreases as groundwater in rock mass is moved to bentnonite buffer by suction, but rock mass is saturated after inflow of groundwater from the far-field area. Stress changes at rock mass are compared to the Mohr-Coulomb failure criterion and the spalling strength in order to investigate the potential rock failure by thermal stress and swelling pressure. Additional simulations are conducted with the reduced spacing of deposition holes. The maximum temperature of bentonite buffer exceeds 100℃ as deposition hole spacing is smaller than 5.5 m. However, temperature of about 56.1% volume of bentonite buffer is below 90℃. The methodology of numerical modeling used in this study can be applied to the performance assessment of coupled THM processes for high-level radioactive waste repositories with various input parameters and geological conditions such as site-specific stress models and geothermal gradients.

• Tunnel and Underground Space
• /
• v.20 no.4
• /
• pp.267-276
• /
• 2010

In modern rock engineering practice, fully grouted rock bolting is actively employed as a major supporting system, so that understanding the behavior of fully grouted rock bolts is essential for the precise design of rock bolting. Despite its importance, the supporting mechanism of rock bolts has not been fully understood yet. Since most of existing analytical models for rock bolts were developed by drastically simplifying their boundary conditions, they are not suitable for the bolts of in-situ condition. In this study, 3-D elastic FE analysis of fully grouted rock bolts has been conducted to provide insight into the supporting mechanism of the bolt. The distribution of shear and axial stresses along the bolt are investigated with the consideration of different boundary conditions including three different displacement boundary conditions at the bolt head, the presence of intersecting rock joints, and the variation of elastic modulus of adjacent rock. The numerical result reveals that installation of the faceplate at the bolt head plays an important role in mobilizing the supporting action and enhancing the supporting capabilities of the fully grouted rock bolts.

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