• Tunnel and Underground Space
• /
• v.28 no.5
• /
• pp.400-425
• /
• 2018

This study presents the research results and current status of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to simulate the coupled hydro-mechanical behavior of fault, including slip or reactivation, induced by water injection. The first research step of Task B is a benchmark simulation which is designed for the modelling teams to familiarize themselves with the problem and to set up their own codes to reproduce the hydro-mechanical coupling between the fault hydraulic transmissivity and the mechanically-induced displacement. We reproduced the coupled hydro-mechanical process of fault slip using TOUGH-FLAC simulator. The fluid flow along a fault was modelled with solid elements and governed by Darcy's law with the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by creating interface elements between two separating rock blocks in FLAC3D. A methodology to formulate the hydro-mechanical coupling relations of two different hydraulic aperture models and link the solid element of TOUGH2 and the interface element of FLAC3D was suggested. In addition, we developed a coupling module to update the changes in geometric features (mesh) and hydrological properties of fault caused by water injection at every calculation step for TOUGH-FLAC simulator. Then, the transient responses of the fault, including elastic deformation, reactivation, progressive evolutions of pathway, pressure distribution and water injection rate, to stepwise pressurization were examined during the simulations. The results of the simulations suggest 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 collaboration and interaction with other research teams of DECOLVAEX-2019 Task B and validated using the field data from fault activation experiments in a further study.

• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.430-440
• /
• 2000

A series of studies such as geological logging data analysis, detailed geological survey, rock mass evaluation, in-situ and laboratory tests, rock strength and mechanical properties of the rock were concerned. The stability of the slope were carried out inorder to design the pit slope and individual benches using the stereographic projection analysis and numerical methods in Roto Pit of Pasir coal field. The bedding plane was one of the major discontinuities in the Roto Pit and the dip of which is about 60$^{\circ}$ in the northern part and 83$^{\circ}$ in the southern part. The dip of bedding becomes steeper from north to south. The plane and toppling failures are presented in many slopes. In laboratory test the average uniaxial compressive strength of mudstone was 9MPa and that of weak sandstone was 10MPa. In-situ test showed that the rocks of Roto north mining area are mostly weak enough to be classified in grade from R2(weak) to R3(medium strong weak) and the coal is classified in grades from R1(Very weak) to R2(Weak). The detailed stability analysis were carried out on 4 areas of Roto north (east, west, south and north), and 2 areas of Roto south(east and west). In this paper, the minimum factor of safety was set to 1.2 which is a general criterion for open pit mines. Using the stereographic projection analysis and the limit equilibrium method, slope angles were calculated as 30∼36$^{\circ}$ for a factor of safety greater than 1.2. Then these results were re-evaluated by numerical analysis using FLAC. The final slope angles were determined by rational described above. A final slope of 34 degrees can guarantee the stability for the eastern part of the Roto north area, 33 degrees for the western part, 35 degrees for the northern part and 35 degrees for the southern part. For the Roto south area, 36 degrees was suggested for both sides of the pit. Once the pit slope is designed based on the stability analysis and the safety measures, the stability of slope should be checked periodically during the mining operations. Because the slope face will be exposed long time to the rain fall, a study such aspreventive measures against weathering and erosion is highly recommended to be implemented.

• Tunnel and Underground Space
• /
• v.10 no.3
• /
• pp.447-456
• /
• 2000

Rock bridge in rock masses can be considered as one of several types of opening-mode fractures, and also it has been known to have a great influence on the stability of structures in rock mass. In the beginning of researching a rock bridge it used to be studied only in characteristics of its behavior, as considering resistance of material itself. However the distribution pattern of rock bridges, which can affect the stability of rock structures, is currently researched with a fracture mechanical approach in numerical studies. For investigating the effect of rock bridges on the development pattern of hydraulic fractures, the author analyzed numerically the stress state transition in rock bridges and their phenomena with a different pattern of the rock bridge distributions. From the numerical studies, a two-crack configuration could be defined to be representative of the most critical conditions for rock bridges, only when cracks are systematic and same in their length and angle. Moreover, coalescence stresses and onset of propagation stresses could be known to increase with decreasing s/L ratio or increasing d/L ratio. The effect of pre-existing crack on hydraulic fracturing was studied also in numerical models. Different to the simple hydraulic fracturing modeling in which the fractures propagated exactly parallel to the maximum remote stress, the hydraulic fractures with pre-existing cracks did not propagate parallel to the maximum remote stress direction. These are representative of the tendency to change the hydraulic fractures direction because of the existence of pre-existing crack. Therefore s/L, d/L ratios will be identical as a function effective on hydraulic fractures propagation, that is, the K$_1$ value increase with decreasing s/L ratio or increasing d/L ratio and its magnification from onset to propagation increases with decreasing s/L ratio. The scanline is a commonly used method to estimate the fracture distribution on outcrops. The data obtained from the scanline method can be applied to the evaluation of stress field in rock mass.

• Tunnel and Underground Space
• /
• v.21 no.6
• /
• pp.471-479
• /
• 2011

The velocity horizontal response spectra using the observed ground motions from the recent 5 macro earthquakes, equal to or larger than 4.8 in magnitude, around Korean Peninsula were analysed and then were compared to the acceleration horizontal response spectra, seismic design response spectra (Reg Guide 1.60), applied to the domestic nuclear power plants, and finally the Korean Standard Design Response Spectrum for general structures and buildings. 102 velocity horizontal ground motions, including NS and EW components, were used for velocity horizontal response spectra and then normalized with respect to the peak velocity value of each ground motion. First, the results showed that velocity horizontal response spectra have larger values at the range of medium natural period, but acceleration horizontal response spectra have larger values at the range of short natural periods. Secondly, the results also showed that velocity horizontal response spectra exceed Reg. Guide 1.60 for longer natural periods bands less than 6-7 Hz. Finally, the results were also compared to the Korean Standard Response Spectrum for the 3 different soil types(SC, SD, and SE soil type) and showed that velocity horizontal response spectra revealed much higher values for the frequency bands below 1.5(SC), 2.0(SD), and 3.0(SE) seconds, respectively, than the Korean Standard Response Spectrum. The results suggest that the fact that acceleration, velocity, and displacement horizontal response spectra have larger values at the range of short, medium, and long natural periods, respectively, can be applied consistently to those form domestic ground motion, especially, the velocity ground motion. Information on response spectrum at such medium range periods can be very important since the domestic design of buildings and structures emphasizes recently medium and long natural periods than short one due to increased super high-rise buildings.

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

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.37 no.1
• /
• pp.9-17
• /
• 2017

Characteristics of aerodynamic damping ratios of a helical $180^{\circ}$ model which shows better aerodynamic behavior in both along-wind and across-wind responses on a super tall building was investigated by an aeroelastic model test. The aerodynamic damping ratio was evaluated from the wind-induced responses of the model by using Random Decrement (RD) technique. Further, various triggering levels in evaluation of aerodynamic damping ratios using RD technique were also examined. As a result, it was found that when at least 2000 segments were used for evaluating aerodynamic damping ratio for ensemble averaging, the aerodynamic damping ratio can be obtained more consistently with lower irregular fluctuations. This is good agreement with those of previous studies. Another notable observation was that for square and helical $180^{\circ}$ models, the aerodynamic damping ratios in along-wind direction showed similar linear trends with reduced wind speeds regarding of building shapes. On the other hand, for the helical $180^{\circ}$ model, the aerodynamic damping ratio in across-wind direction showed quite different trends with those of the square model. In addition, the aerodynamic damping ratios of the helical $180^{\circ}$ model showed very similar trends with respect to the change of wind direction, and showed gradually increasing trends having small fluctuations with reduced wind speeds. Another observation was that in definition of triggering levels in RD technique on aerodynamic damping ratios, it may be possible to adopt the triggering levels of "standard deviation" or "${\sqrt{2}}$ times of the standard deviation" of the response time history if RD functions have a large number of triggering points. Further, these triggering levels may result in similar values and distributions with reduced wind speeds and either may be acceptable.

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

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