This study was performed to suggest to suggest suitable design conditions of water curtain system through analysis on pressure down in boreholes by hydraulic tests carried out I construction fields for underground oil storages. The influence by hydraulic conductivities of rock mass around boreholes on pressure down in boreholes was analysed. The relationship between array of boreholes and their pressure down was also analysed. Groundwater flow analysis on crude oil and LPG storages was carried out to evaluate results of field tests and to investigate distribution of hydraulic gradient in rock mass around cavern using finite difference method. As the results, hydraulic tests showed that pressure down in boreholes was inverse proportional to the hydraulic conductivity of surrounding rock mass. The rate of pressure down of boreholes was not influenced by water curtain system more than 20m over cavern and was proportional to installation interval of boreholes. The hydraulic gradient in rock mass around cavern was proportional to distance and interval of boreholes and its value was not satisfactory to oil tightness condition in case of no water curtain system.

Field instrumentation and numerical analysis by the finite difference method were applied to estimate the relaxed zone in a subway tunnel of shallow depth in soft rock, excavated by NATM. The convergence and ground displacement can be used to estimate the deformation behavior and the relaxed zone. Parameters for the several models previously suggested were measured using regression analysis techniques adopting a function of time and the face advance. The estimated relaxed zone by the MPBX and FDM analysis were 1.5~3.0 m and 1.5~2.0 m, respectively. It was concluded that the visco-elastic model and the time-dependent elasto-plastic model correlate very well ($r^2$>0.9) with results of the numerical analyses.

Rock contains discontinuities at all scales. These discontinuities make rock behave in a complex way. This paper discusses a new approach to underground design based on the theory of rock fracture mechanics. The mechanism of deformation and failure of coal was studied by observing the distributions of length, orientation and spacing of the pre-existing as well as stress-induced cracks. Different types of crack information. The crack information is dependent on the scale used. The cracks propagate along the intersections of the pre-existing cracks, and both extensile and shear crack growth occur depending on the direction of the load relative to the bedding planes. An analytical model that takes into account both shear and extensile crack growth was developed to predict the nonlinear stress-strain behavior of coal including strain-hardening and strain-softening.

In the case of a fire disaster in a uni-directional road tunnel, it is important to determine the critical ventilation velocity to prevent the backflow travelling toward the tunnel exit where vehicles are stopped. The critical ventilation velocity is horizontal velocity to prevent hot smoke from moving toward the tunnel exit. According to Froude modelling, the model tunnel whcih was 300mm in diameter and 21 m in length was made of acryl tubes. Inner section of acryl tubes was clothed with polycarbonate. 1/20 scaled model vehicles were installed to simulate the situation that vehicles are stopped in the tunnel exit. Methanol in a pool type burner was burned in the middle of tunnel to simulate a fire hazard. In this study, the basis of determining the critical ventilation velocity is the ventilation flow rate that is able to maintain the allowable CO concentration in the tunnel section. We assumed that the allowable CO concentration was backflow dispersion index. Futhermore, We intended to find out CO distribution and temperature distribution according as we changed ventilation velocity. The results of this study were that no backflow happened when ventilation velocity was 0.52 m/s in the case of 5.75 kW. If we adapt these results of a fire disaster releasing 10MW heat capacity in real tunnel which is 400m in length, no backflow happens when ventilation velocity is 2.31m/s. After we figured out dimensionless heat release rate and dimensionless ventilation velocity of model test and those of real test to verify experimental correctness, we tried to find out correlation between experimental results of model tunnel and those of real tunnel.

A three dimensional boundary element method-BEAP3D was applied to the stability analysis of the mine openings not only to improve the stability during mining operations but also to serve the evaluation of the mine openings for further utilization. Stability analysis on the stability of the room-and-pillar stopes underneath of the old mine openings and the openings to be created by the newly proposed sublevel stoping method at the Nowhado Pyrophyllite Mine, showed that rock mass around the old and new stopes would be stable. Six stopes of a sublevel stoping designed for the Choongmu Limestone Quarry would be stable, too. A sublevel stoping method consisting of six stopes was similarly suggested for the Keumpyung Quartzite Mine. The stability can be guaranteed through out six stopes. Since mining starts from the bottom 1st sublevel to the uppermost sublevel, the safety of the stopes will improve together with the mining process. It would highly be recommended to investigate in-situ rock properties and the rock stresses for future studies. Even though the rock around the uppermost part and bottom of all the stopes have a very high factor of safety, spot reinforcements such as rock bolting would be recommended to mitigate the intermediate and minor principal stresses acting in a tensile mode.

Comprehensive monitoring system for the safe and economical excavation of underground opening has been established by employing the 3 independent models each of which can i) predict the ultimate convergence, ii) assess the in-situ stresses and the elastic modulus of excavating rock, iii) calculate the time-dependent opening behavior with respect to the face advance rate and support pressure at the equilibrium state. Accuracy of each model has been verified through illustrative examples. The step-by-step procedures of comprehensive monitoring system for analyzing the rock behavior and the optimum support installation has been explained. The capability and applicability of this system to the practical excavation also has been discussed.

The reliable estimation of the system parameters and the accurate prediction of the system behavior are important to design underground structures safely and economically. Especially, the elastic modulus and the in-situ stresses are very important parameters in predicting the behavior of the underground structure. Therefore, the back analysis using the field measurement data is developed to determine accurately the elastic modulus and the in-situ stresses of the underground structural system in this study. A back analysis using the combined finite and boundary element is developed. It can consider the far field boundary condition and is efficient in computation. In this study, a back analysis is performed to predict behaviors of underground structures for the real construction site. The comparison between the results of the back analysis with field measurement data and the obtained material properties from the field test shows good agreement for the real construction site.

Nemerical algorithms were developed to analyze the behavior of the double lining as well as ground mass separately or simultaneously. A lining interface element was especially developed, verified and applied to the study on the coupled interaction of shotcrete and the concrete lining. It could be known fro parameter studys on double lining support systems that as the contact surface between shotcrete and concrete lining was rougher, the more decreased bearing capacity against the cracking of the system. If the thickness of the shotcrete increased, the bearing capacity of the double lining also increased linearly with the thickness. If the thickness of the concrete lining increased, the bearing capacity of the double lining had the relationship of the characteristic S-shape of a sigmoid function with the thickness. When the thickness increased over a given value, it was not useful to increase more the thickness because bearing capacity had no remarkable change. It could be concluded that the behavior of the shotcrete and concrete lining was generally reversed before and after the ratio of horizontal to vertical earth preassure of 1.0 and 0.5 respectively. Therefore, we could guess that the movement which two shotcrete and concrete lining deflect toward each other around the crown caused a friction between two linings and thus this disadvantageous effect could contribute to reducing the bearing capacity against the cracking.