Fusion technology is a key to maximize innovative potential of geoscience for many challenging issues today that require integrated multi-disciplinary approach. Successful fusion technological advance can be achieved when interdisciplinary cooperation is firmly established. In order to establish firm the context of inter-disciplinarity that is still feeble, it is urgent to continuously develop geoscientific models and systematic infra for interdisciplinary cooperation such as well-prepared geo-spatial database and knowledge base network that can support multi-lateral cooperation between multiple disciplines and multi-phase international cooperation.

The visualization of three dimensional geophysical data is forcing a revolution in the way of working, and allowing the discovery and production of hydrocarbons at much lower costs than previously thought possible. There are many aspects of this revolution that are behind the scenes, such as the database structure, the storage and retrieval of data, and the exchange of data among programs. Also the user had changes where the interpreter (or manager, or processor) actually looks at and somehow interacts with the data. The use of opacity in volume rendering, and how its judicious application can assist in imaging geologic features in three dimensional seismic data. This revolutionary development of new technology is based on the philosophy of synergy of inter-disciplines of the oil industry. Group interaction fostered by large room visualization environments enables the integration of disciplines we strive for, by putting the petrophysicist, geologist, geophysicist, and reservoir engineer in one place, looking at one image together, without jargon or geography separating them. All these tools developed in the oil industry can be applied into the civil engineering industry also such as the prior geological and geophysical survey of the constructions. Many examples will show how three dimensional geophysical technology might make a revolution in the oil business industry now and in future. This change can be considered as a fusion process at data, information, and knowledge levels.

Fire accidents in tunnels and underground installations have resulted in significant loss in terms of life and costs. Several large research projects are launched in Europe to find cost effective ways to upgrade and make such installation more robust against accidents from fires and explosions. For one single tunnel, severe fire accidents are rare and operators and users will have very limit experience, when they do occur. There is a trend to solve the problem with high tech instrumentation and monitoring techniques and to rely on the fire brigade to play a major role in the rescue operation. These precautionary measures are very difficult to validate for the severe fires. For protection against the most critical fires, the safety has to be considered in the design of the tunnel and the technical installation validated so that it operates as indented.

Rock mass characterization is an integral part of rock engineering design. Much of the information for rock mass characterization comes from field fracture mapping and data collecting. This paper describes two technologies that can be used to assist with the field mapping and data collecting activities associated with rock mass characterization: digital image processing and 3D laserscanning. The basis for these techniques is described, as well as the results of field case studies and an analysis of the error in estimating fracture orientation.

The global demand for underground facilities has increased substantially in the past decades, and a substantial number of underground projects have had to deal with challenging ground conditions in urban environments. Particularly challenging are weak and unstable water bearing soils. Advancements in shielded TBM tech-nology have led to significant improvements regarding the ability to control ground deformations in soft ground. Nonetheless, ground collapse may occur even when the most advanced TBM designs are employed if unexpected adverse ground conditions are encountered or if insufficient stabilizing pressure is transferred to the tunnel face. This paper reviews common approaches for face stability and face pressure transmission calculations, and provides an overview of some of the latest technological developments and considerations for soft ground TBM applica-tions.

It is only possible through the image analysis of borehole wall and the core recovered from borehole constructed in rock mass that the real information about geologic characteristics in rock mass is directly obtained in primary research. Monitoring apparatus with multi-functional utility has implemented and applied in-situ condition for finding the geologic condition of target area. But, this apparatus is very expensive to be applied at the risk of loss during monitoring and cause hard work for moving them to the determined position. This paper shows the underground imaging from the borehole information obtained by a borehole camera with the simple utility and low cost enough to investigate the characteristics of borehole wall. Monitoring for this has been done in open-pit mine located at the northeastern part of Fukuoka Prefecture in Japan, and finally the three dimensional imaging of geological discontinuity was discussed relative to the field condition.

A new algorithm was developed to interpret joint orientations from a pair of images of the rock slope to overcome the limitation of photographing direction as in the parallel stereophotogrammetric system and to maximize the range of image measurement. This algorithm can be regarded as a modified multistage convergent photographing system. To determine camera parameters in the perspective projection equation that are the major elements in the photogrammetric technique, a new concept was developed by using three ground control points and single ground guide point. This method could be considered to be very simple when compared with other existing methods which use a number of ground control points and complicated analysis processes.

Profiles of naturally fractured surfaces of three sedimentary rock samples were plotted from the measured data using a mechanical profilometer. Fractal dimension of these profiles were computed and statistical F-test indicates that fractal dimension (FD) values can be used as a parameter for distinguishing the rock types. The comparison between FD values and a commonly used profile-roughness parameter called the Mayer's $Z_2$ parameter shows the superiority of the FD values as roughness estimator. Two-dimensional fractal roughness parameters of the same naturally fractured rock surfaces were also studied from their scanning electron microscopic (SEM) images at various magnification levels. The most suitable level of magnification of the SEM images for the study of the 2-D fractal roughness parameter was identified. The values of 2-D fractal roughness parameter for three different rocks were also compared using different methods of fractal dimensioning.

The study tasked by Ryan Lode Mines, Inc. involved analysis of pit slope stability for two future pits, Ryan and Curlew. A geological discontinuity survey and groundwater information were obtained through a field program. Twenty core logs containing geotechnical information were used for rock mass classification. The kinematic analysis was performed based on a friction angle (${\Phi}=35^{\circ}$), the distribution of geological structures, and a dry slope condition. Factors of safety of pit slopes in two future mines were determined using the limit equilibrium method. The mine slopes and benches designed by Mine Development Associates (MDA) were analyzed. The analysis indicated that both pits should have an overall safety factor above 1.0, provided the slopes are kept dry. However, slopes in both pits exceeding 91.4 m (300 ft) high will become critical, when water fills the cracks and discontinuities.

In this paper, simulated annealing technique was used to estimate the rock joint characteristics, RMR(rock mass rating) values, to overcome the defects of ordinary kriging. Ordinary kriging reduced the variance of data, so lost the characteristics of distribution. Simulated annealing technique could reflect the distribution feature and the spatial correlation of the original data. Through the comparisons between three times simulations, the uncertainty of the simulation could be quantified, and sufficient results were obtained.

Experimental blast studies were carried out in a limestone quarry to study the effects of blasting on structures. To have an in-depth understanding of the possible relation between parameters like vibrations, frequency and scale distance and ten trial blast were conducted. 29 monitoring stations were located in such a pattern to give a true representation of blast induced vibrations for the entire mining in this area. The vibrations were monitored in the vicinity of structures surrounding the quarry in the direction towards the village, road, railway line, office building, etc. Scale distances were determined to identify the maximum charge permissible to cause damage to structures.

In recent, as Republic of Korea has been interested in environmental problem increasingly and became a member of many organizations or institutions related to environmental preservation such as a Ramsar convention, fundamental and completed methods to prevent ground water's drying up and leakage in tunnel excavation are requested. In this paper, we have studied the anticipated problems by tunnel excavation under the wetland and described the effective designed method to maintain the wetland's ecosystem environment. To accomplish this purpose, firstly, we investigated the wetland's ecosystem, ground's hydraulic properties and analysed the foreign similar case for tunnel excavation near the wetland. And by numerical analysis, we analyzed the runoff and infiltration quantity of water and hydraulic behaviour properties by saturation and unsaturation concept in rock mass and wetland. Finally, we established the effective countermeasure to minimize the ecosystem's bad influence by tunnel excavation.

In the construction of mountain tunnels, reaction forces of the legs of steel arch supports against the ground are often expected to support the ground being excavated. In these cases, a stress concentration occurs in the ground directly under the support legs. If the bearing capacity of the ground is insufficient or displacement is not effectively constrained, the local failure of the ground under the support legs or settlement of the tunnel supports due to large deformation could result. It is therefore necessary to reinforce the support legs to reduce settlement. As a means of reducing settlement, wing-ribbed steel arch supports are well used. In this study, with the aim of finding a way to quickly reduce the settlement of steel arch support legs, effectiveness of a new type of wing ribs to reinforce steel arch supports was investigated through laboratory testing.

It is the common practice to reinforce excessively the secondary tunnel lining due to the lack of rational insights into the ground loosening loads. The main load of the secondary lining for drained-type tunnels is the ground loosening. The main cause of the load for secondary tunnel lining is the deterioration of the primary support members such as shotcrete, steel ribs, and rockbolts. Accordingly, the development of the analysis model to consider the ground-primary supports-secondary lining interaction is very important for the rational design of the secondary tunnel lining. In this paper, the interaction is conceptually described by the simple mass-spring model and the load transfer from the primary supports to the ground and the secondary lining is showed by the characteristic curves including the secondary lining reaction curve for the theoretical solution of a circular tunnel. And also, the application of this model to numerical analysis is verified in order to review the potential tool for practical tunnel problems with the complex conditions like non-circular shaped tunnels, multi-layered ground, sequential excavation and so on.

Discrete joint network approach has widely been used to investigate the hydraulic behavior of jointed rock masses. In general, joints will undergo deformation due to stress redistribution induced by construction of underground openings, hence joint aperture is often assumed to have a probability distribution rather than to be a constant value. In real situations, however, it is more reasonable to take into account the effect of stress change on aperture values by calculating joint deformation. In this report, a mechanical process has been developed to determine the joint opening or closure based on a statistically generated joint network model. By performing numerical analyses, some significant results on the hydro-mechanical behavior of jointed rock masses have been summarized.

Rock bolts are widely used as a supplementary roof support system in hard-rock mining since a long time. Since the performance of fully grouted passive bolts depends on bond strength, in the present investigation extensive laboratory pull-out as well as push-out tests were conducted varying the bolt diameter, length and cement-water mixing ratios of grout. The load-displacement curves were developed and were verified with the numerical results obtained from finite element analysis using ALGOR software. Numerical models were validated for push-out tests and a detailed analysis was carried out to know the displacement, stress, strain distribution along the bolt.

The precise prediction of ground displacement plays an important role in planning and constructing tunnels. In this study, an equation for predicting the surface and crown settlement is suggested by examining the theories of ground movement caused by tunnel excavation. From the 3D numerical modeling, the reinforcement effect of UAM (Umbrella Arch Method) is quantitatively analyzed with respect to deformation modulus and overburden. By using a regression technique for the numerical results, an equation for predicting the settlement is suggested.

Blind hydraulic backfilling is a commonly used technique for subsidence control of the strata over unapproachable waterlogged underground excavations. In this investigation model studies on all the three variants of this technique, namely, hydro-pneumatic or air-assisted gravity backfilling, pumped-slurry backfilling and simple gravity backfilling, have been carried out in fully transparent models of the underground excavations. On examination of the filling process, it was revealed that in all the three cases, the basic process of filling occurs by sand transport along one or more meandering channels. The relative influence of sand, water and air flow rates on the area of filling from a single inlet point and the hydraulic pressure loss per unit length were studied in details. In hydro-pneumatic backfilling process, the air bubbles while moving upward through the meandering channels provide an additional buoyant force over and above the available hydraulic head. In this way the area of filling from a single borehole may be quite large even at small flow rates of water. During actual field implementation the injected air, if not released completely from the rise side holes, may cause troubles by way of creating potholes on the surface. The pumped-slurry technique has shown its capability of filling a relatively larger area at faster rate, especially when high-volume, low-pressure method was selected. But simple gravity filling was also found to be equally effective method as slurry pumping, especially when flow rates were high. In the second and third method discussed above, examination of variations of injection pressure was also done and its relation with physical phenomenon was also attempted. Some empirical relationships were also developed using multivariate regression with a view to help the practicing engineers.

Fractures in the form of micro cracks are commonly found in natural rocks. A rock behaves in a complex way due to fracture; in particular, the anisotropic strength of a rock material is significantly influenced by the presence of these fractures. Therefore, it is essential to understand the failure mechanism of a fractured rock. In this study, a fractured rock is formulated in terms of fabric tensor based on geometric and mechanical simplifications. In this way, position, density and shape of fractures can be determined by the fabric tensor so that rocks containing multi-fractures can successfully be modeled. Also an index to evaluate the degree of anisotropy of a fractured rock is proposed. Hence, anisotropic strength of a rock containing fractures under uniaxial compression condition is estimated through a series of numerical analyses for the multi-fractured model. Numerical investigations are carried out by varying the fracture angle from $0^{\circ}\;to\;90^{\circ}$ and relationship between uniaxial compression strength and the degree of anisotropy is investigated. By comparing anisotropic strength of numerical analysis with analytic solution, this study attempts to understand the failure mechanism of rock containing fractures.

Mechanical behavior of rock joints usually can be characterized by small-scale laboratory shear tests due to economical and technical limitations, but their applicability to the behaviour of rock mass has been always questioned by a number of researchers because of scale effect. Though there have been several researches regarding the scale effect, it has been a controversial problem how to apply the result of small-scale laboratory shear test directly to field design from different conclusions among researchers. In order to grasp the trend of scale effect of shear behavior, a series of direct shear tests on replicas of natural rock joint surfaces made of gypsum cement with different size and roughness were conducted and analyzed. Result showed that as the size of the specimen increased, average peak shear displacement increased, but average shear stiffness and average peak dilation angle decreased. As for the dependency of scale on shear strength, the degree of scale effect was dependent on normal stress and roughness of rock joint. For the condition of low normal stress and high roughness, decrease of average peak shear strength with increasing size of joint was evident.

Excavation by TBM can be characterized by a rock-machine interaction during the cutting process on a small scale, but on a large scale the interaction between the rock mass and TBM becomes very significant. For the planning and evaluation of TBM tunnelling it needs to understand rock fracture mechanism by a cutter or cutters on a small scale, and to estimate penetration rate, advance rate and utilization on a large scale. In this study rock chipping mechanism due to cutter-penetration is analysed by numerical simulation, showing that rock chipping is mainly occurred by tensile failure. Also, through the analysis of factors that affect on TBM procedures in various assessment systems, it is determined that the key elements that should be considered in the planning and evaluation of TBM tunnelling are classified into rock properties, the geological structures and properties of rock mass, and the structural and functional specifications of the machine. The user-friendly assessment tool is developed, so that penetration rate, advance rate and TBM utilization are evaluated from various input data. The tool developed in this study can be applied to a practical TBM tunnelling by understanding TBM tunnelling procedures.

Ground and rock mass considered in tunnelling have characteristics such as uncertainty, heterogeneity and structural complexity because they have been formed undergoing various geological events for a long period. So, it is difficult for engineers to predict behaviors of rock mass in tunneling. In the paper the authors describe the development of an integrated expert system prototype for site investigation, design and construction in tunnelling and introduce the case applying this system to the tunnel construction site under construction. Geostructure Research Group in Korea Institute of Construction Technology (KICT) has developed the system during the past 4 years. The system mainly consists of several modules which is related to the design, construction and management of tunnelling. The test site, Neung-dong tunnel is located in Ulsan, Korea. The geology map shows it may confront big fault zone whose width is over kilometres. With the networking system of ITIS, various information of face mapping, monitoring and other construction task can be transmitted into the database and GIS Server at real time. And necessary analyses can be carried out with the modules equipped in the system.

This paper presents a neural network based approach to disseminating information relating to experimental and field observations in engineering. Although the methodology is generic and can be applied to many areas of engineering science, attention is focussed here solely on geotechnical engineering applications. Field data relating to the settlement of foundations presented by Burland and Burbidge (1985) which led to their well known equation for calculation of settlement, now included in most text books, is re-visited. A part of the data, chosen randomly, is used to train an Artificial Neural Network (ANN), which relates foundation settlement to various causes as identified by the authors. Predictions are made for situations for which data were not used in training. These indicate sufficient accuracy when compared to the original field data. Accuracy of predictions is further improved when all the data are included in the training set. The finally trained ANN is shown to represent these data more accurately than the Burland and Burbidge equation. Based on the above heuristic example, an ANN is presented as an alternative to developing equations and design rules in geotechnical engineering practice. Significant advantages are shown to arise by using this methodology. Ease of updating the ANN, as and when additional data becomes available, being the most important one. Loss of transparency, however, seems to be the main disadvantage.

In this paper, a mathematical framework based on a homogenisation technique to simulate 'umbrella arch reinforcement system' (UARS) and its implementation into a 3D Finite Element program that can consider stage construction situations are presented. The constitutive model developed allows considering the main design parameters of the problem and only requires geometrical and mechanical properties of the constituents. Additionally, the use of a homogenisation approach implies that the generation of the Finite Element mesh can be easily produced and that re-meshing is not required as basic geometrical parameters such as the orientation of the pipes are changed. The model developed is used to simulate tunnelling with the UARS. From the analyses, the effects of the main design parameters on the elastic and the elastoplastic analyses are considered.

It can be said that rock mass properties are characterized not by a mean value but by values with variation due to its characteristic uncertainty. This characteristic is one of the most important parts for the design of underground structures, but yet to be fully examined. Stochastic finite element method (SFEM) has been developed in order to take the randomness of structural systems into account. Using SFEM, the response variability of structural system can be obtained and it leads probabilistic stability of structure to be analyzed. In this study, displacements response variability of circular opening with hydrostatic stress field are analyzed in terms of rock mass properties having a certain mean and a standard deviation using the SFEM. The analyzed response variability shows that the necessity of probabilistic stability analysis of underground structures using reliable mean value and standard deviation of deformation modulus.

We considered the variation of elastic wave velocity due to the anisotropy of rock materials and stress level for acoustic emission (AE) source location in cylindrical rock specimens. Elastic wave velocity and AE were measured for Keochang granite and Yeosan marble under various axial stresses and confining pressures. Partition approximation method was suggested and it was compared with the difference approximation method and the least square method.

Laboratory testing of representative rock specimens is of fundamental necessity for the successful design and/or assessment of facilities associated with many kinds of underground exploitation, including the geological disposal of radioactive nuclear waste. As a fundamental and systematic study, a series of measurements of the physical, mechanical and hydraulic properties of Inada granite and Shirahama sandstone, two rock types that are widely available in Japan, have been performed. This paper presents the results of a study of the effective porosity, density, compressive and shear wave velocity, unconfined compressive strength and permeability of the two rocks. The anisotropy and the effects of confining pressure on the permeability of the rocks, as well as the relationships among the physical, mechanical and hydraulic properties, are also investigated and discussed.

Pipe jacking is a name for a method to excavate a tunnel by pushing pipe into the ground from an especial pit. Size of tunnels in this method is different from under 900mm (microtunneling) to more than 3,000mm. Method of excavation is also different from hand digging to use of any kind of tunnel boring machines such as slurry and earth pressure balance (EPB) machines. Slurry pipe jacking was firmly established as a special method for the nondisruptive construction of the underground tunnels in urban area. During the pipe jacking and microtunneling process, the jacking load is an important parameter, controlling the pipe wall thickness, need to and location of intermediate jacking station, selection of jacking frame and lubrication requirements. The main component of the jacking load is due to frictional resistance. In this paper the skin friction between pipe surface and surrounding condition also lubricant quality based on a few fundamental tests, were considered. During this study unconfined compressive strength test, dynamic friction measurement test and direct shear box test were raised for one of the largest diameter slurry pipe jacking project in Fujisawa city in Japan. It could be concluded that in slurry pipe jacking, prediction of frictional forces are mainly dependent on successful lubrication, its quality and lubricant strength parameters. Conclusions from this study can be used for the same experiences.

The stress corrosion index of Kumamoto andesite are evaluated by two types of testing method. One is the uniaxial compression test under various water vapor pressures, and the other is the double torsion (DT) test under a constant water vapor pressure. For the uniaxial compression tests, the uniaxial compressive strength increases linearly with decreasing water vapor pressure on the double logarithmic coordinates. As the results, the stress corrosion index obtained is estimated 44. On the other hand, in the DT test, the relaxation (RLX) test and the constant displacement rate (CDR) test were conducted. For the CDR test, as the displacement rate of loading point increases, the crack velocity increases. However, the fracture toughness is constant regardless of the change in displacement rate and the average fracture toughness is evaluated $2.07MN/m^{3/2}$. For the RLX test, the crack velocity-stress intensity factor curves are smooth and linear. The stress corrosion index estimated from the curves is 37. Comparing stress corrosion indexes in the uniaxial compression test and the DT test, there is no significant difference in these values, and they are considered to be in coincident each other regardless of testing methods. Therefore, it is concluded that stress corrosion is one of material constants of rock.

The stereographic projection method accesses the stability of slope roughly, so it is difficult to apply its result in the technical designing of rock slopes. To solve this difficulty, using the limit analysis method, we quantified the stereographic projection with the safety level graded in the daylight envelope, which represents the failure possible area in the net.

For Taejon LNG Pilot Cavern being constructed to verify the technical aspects for storing LNG in lined rock cavern, various numerical studies were carried out to estimate the temperature profile and to understand thermo-mechanical behaviour in the rock around the cavern. With the help of Claesson's analytical solution and numerical models, the extent of zero degree isotherm and possible boil-off rate of gas to be stored were estimated. Even though the tensile stress by cooling down is very large compared to the tensile strength of the rock, it has been shown that possible rock yielding might bring about the dramatic reduction of the stress.

As there are many advantages on underground caverns, such as safety and operation, they can also be used for gas storage purpose. When liquefied gas is stored underground, the cryogenic temperature of the gas will affect the stability of the storage cavern. In order to store the liquefied gas successfully, it is essential to estimate the exact temperature distribution of the rock mass around the cavern. In this study, an analytic solution and a conceptual model that can estimate three-dimensional temperature distribution around the storage cavern are suggested. When calculating the heat transfer within a solid, it is likely to consider the solid as the intersection of two or more infinite or semi-infinite geometries. Therefore heat transfer solution for the solid is expressed by the product of the dimensionless temperatures of the geometries, which are used to form the combined solid. Based on the multi-dimensional transient heat transfer theory, the analytic solution is successfully derived by assuming the cavern shape to be of simplified geometry. Also, a conceptual model is developed by using the analytic solution of this study. By performing numerical experiments of this multi-dimensional model, the temperature distribution of the analytic solution is compared with that of numerical analysis and theoretical solutions.

A two-dimensional BEM code, $FRACOD^{2D}$, was applied to simulate fracture initiation and propagation processes in a rock pillar during an in situ heater test of a rock pillar planned at the $\"{A}sp\"{o}$ Underground Rock laboratory of SKB, in Southern Sweden. To take the advantage of conventional BEM for simulating fracturing processes, but without efforts for domain integral transformation, a hybrid approach is developed to simulate the fracturing processes in rock pillar under coupled thermo-mechanical loading. The code FRACOD was used for simulating the fracture initiation and propagation processes with its boundary tractions reflecting the effects of the initial and redistributed thermomechanical stresses in the domain of interest at multiple excavation and heating steps were produced by a special algorithm of stress inversion, based on resultant thermo-mechanical stress fields at each excavation and heat loading step by a FEM code without considering fracturing processes. This hybrid approach can take the advantages of both types of numerical methods and avoids their shortcomings for fracturing process simulation and domain effects, respectively. In this paper, we present the hybrid approach for the stress, displacements, and fracturing processes at sequential excavation and heating steps of the in situ heater test as a predictive modelling, the formulation of the fracturing models and the predictive results. Two sections of borehole depth, 0.5 m and 1.5 m below the tunnel floor are considered. The pillar area is modelled with the FRACOD and the stress field produced by excavation and heating is transferred with corresponding boundary stresses. From the modelling results, the degree of fracturing and damage are evaluated for 120 days of heating. Dominated shear fracturing in the vicinity of the central pillar was observed from the models at both sections, but spalled area appears to be limited. Based on the modelling results, a sensitivity study for the effect of pre-existing fractures in the vicinity of the holes is also conducted, and the initiation and evolution of EDZ around the deposition holes are investigated using this particular numerical technique.

In this study, it was purposed to develop the new method for the prediction of pollutant concentration in road tunnels. The new method was the use of artificial neural network with the back-propagation algorithm which can model the non-linear system of tunnel environment. This network system was separated into two parts as the visibility and the CO concentration. For this study, data was collected from two highway road tunnels on Yeongdong Expressway. The tunnels have two lanes with one-way direction and adopt the longitudinal ventilation system. The actually measured data from the tunnels was used to develop the neural network system for the prediction of pollutant concentration. The output results from the newly developed neural network system were analysed and compared with the calculated values by PIARC method. Results showed that the prediction accuracy by the neural network system was approximately five times better than the one by PIARC method. ill addition, the system predicted much more accurately at the situation where the drivers have to be stayed for a while in tunnels caused by the low velocity of vehicles.

Tunnelling in poor and heterogeneous ground is a difficult task. Even with a good geological investigation, uncertainties with respect to the local rock mass structure will remain. Especially for such conditions, a reliable short-term prediction of the conditions ahead and outside the tunnel profile are of paramount importance for the choice of appropriate excavation and support methods. The information contained in the absolute displacement monitoring data allows a comprehensive evaluation of the displacements and the determination of the behaviour and influence of an anisotropic rock mass. Case histories and with numerical simulations show, that changes in the displacement vector orientation can indicate changing rock mass conditions ahead of the tunnel face (Schubert & Budil 1995, Steindorfer & Schubert 1997). Further research has been conducted to quantify the influence of weak zones on stresses and displacements (Grossauer 2001). Sellner (2000) developed software, which allows predicting displacements (GeoFit$\circledR$). The function parameters describe the time and advance dependent deformation of a tunnel. Routinely applying this method at each measuring section allows determining trends of those parameters. It shows, that the trends of parameter sets indicate changes in the stiffness of the rock mass outside the tunnel in a similar way, as the displacement vector orientation does. Three-dimensional Finite Element simulations of different weakness zone properties, thicknesses, and orientations relative to the tunnel axis were carried out and the function parameters evaluated from the results. The results are compared to monitoring results from alpine tunnels in heterogeneous rock. The good qualitative correlation between trends observed on site and numerical results gives hope that by a routine determination of the function parameters during excavation the prediction of rock mass conditions ahead of the tunnel face can be improved. Implementing the rules developed from experience and simulations into the monitoring data evaluation program allows to automatically issuing information on the expected rock mass quality ahead of the tunnel.

Specification of reinforcement method was suggested according to the ground condition and tunnelling environment such as adjacent building and surface settlement. Tunnel database consists of 8 different groups of data according to the tunnel construction situations and major problems of ground. A tunnel countermeasure expert system based on client/server system was developed with on-line. The expert system provides proper solution to the each construction sites backing up the information of the tunnelling and ground information through Internet. The effective factors of tunnel construction were shown by the analyzing relationship and partial relationship between face stability and RMR factors. This study will be very helpful to make the most of in-situ data and suggest proper applicability of tunnel reinforcement system escaping from the dependence of some experienced experts for the absent of guide.

A definite shut-in pressure in hydraulic fracturing techniques is needed for obtaining the correct information on the in-situ stress regimes in rock masses. The relation between the behaviour of hydraulically induced fractures and the condition of remote stress is considered to be major reasons of an ambiguous shut-in pressure in hydraulic fracturing pressure-time history curves. This paper describes the results of a series of numerical analyses carried out using UDEC(Universal Distinct Element Code, Itasca), which is based on the discrete element method, to compare several methods for determining the shut-in pressure during hydraulic fracturing. The fully coupling of hydraulic and mechanical analysis was applied, and the effects of four different discontinuity geometries in numerical modelling have been investigated for this purpose. The effects of different remote stress regimes and different physical properties on hydraulic fracture propagation have been also analyzed. Several methods for obtaining shut-in pressure from the ambiguous shut-in curves have been applied to all the numerical models. The graphical intersection methods, such as (P vs. t) method, (P vs. log(t)) method, (log(P) vs. log(t)) method, give smaller values of the shut-in pressure than the statistical method, (dP/dt vs. P). Care should be taken in selecting a method for shut-in pressure, because there can be existed a stress anomaly around the wellbore and fracturing from the wellbore by a constant flow rate may have a more complicate mechanism.

A Natural cavities were found at shallow depth during construction of a huge bridge in Moon-Kyung, Korea. The distribution patterns of cavities in the Moon-Kyung limestone were investigated carefully with a supplementary field job such as a structural geological survey, a geophysical survey, and a rock mechanical test in laboratory or field. A structural geological mapping produced a detail geological map on this area. It suggested that there were three faults in this area, and these faults had an influence on the mechanism of natural cavities. Among many kinds of geophysical surveys, an electrical resistivity prospecting was applied firstly on the specific area that was selected by results from the geological survey. Many evidences for cavities were disclosed from this geophysical data. Therefore, a seismic tomography was tested on the target area, which was focused by results from the electrical resistivity prospecting and was believed to have several large cavities. A distinct element numerical simulation using the UDEC was followed on the target area after completing all of field surveys. Data from field tests were directly dumped or extrapolated to numerical simulations as input data. It was verified from numerical analysis that several natural cavities underneath the foundation of the bridge should be reinforced. Based on the project result, finally, most of foundations for the bridge were re-examined and the cement grouting reinforcement was constructed on several foundations among them.

The relationship between the reflection number obtained from seismic reflective survey and the rock strength value obtained from TBM excavation is examined, and the procedure of the conversion from the reflection number to the rock strength value is proposed. Subsequently, geostatistical method is employed to evaluate the rock properties ahead of the tunnel face and around the tunnel with good precision, using both the seismic reflective survey data and the TBM driving data for the purpose of the tunnel driving and enlargement. The applicability of this evaluation method is examined at the actual tunnel site.

Borehole radar is one kind of GPR, but it can be used in deep boreholes, and it has many advantages compared with low frequency borehole EM tools, and surface GPR. We have developed various techniques on borehole radar for environment study. The hardware development includes broadband radar system with the functions of polarimetry and inteferometry. By using these systems, we tested the measurements to applications such as subsurface fracture characterization, subsurface cavity detections. In this paper, we will describe the advantages of the advanced radar technology for environment studies, and show some experiment results.

Remote sensing can provide invisible information in addition to acquire wide-view image data from space. Synthetic Aperture Radar (SAR) transmits microwave to the earth from a satellite and collects the reflected echo from the surface. Interferometric processing of SAR data can detect the subtle land deformation. The information of the surface movement by SAR is useful to monitor the volcanic activity, extended subsidence of urbanized area and the prediction of the earthquake caused by crustal deformation, and it complements the conventional levelling and GPS technique. PSInSAR (Permanent Scatterers Interferometric SAR) is one of interferometric techniques to be applied to practical projects in Japan. In this paper, the projects of land deformation monitoring are shown after the explanations of the PSInSAR principle. Tokai earthquake risk assessment is the first example. PSInSAR detects the subduction of crustal deformation of the adjacent area of new assumed epicenter region of the Tokai Earthquake. The extended subsidence of the urbanized area was implemented by using Japanese satellite data i.e. JERS that has so much data the surrounding of Japan as the archive. We examine the relationship between the geological structure and settlement at Nohbi basin including Nagoya city.

Synthetic seismograms recorded with a logging-while-drilling (LWD) tool in the presence of slow formation are computed by the discrete wavenumber method. Monopole, dipole, and quadrupole logging tools are simulated with a source center frequency of 4 kHz. The modes in the responses are identified and characterized with time and frequency semblance plots. Numerical results show that, to obtain the formation shear velocity, we need to correct the peak velocities of the multipole modes in the semblance plots by using analytical dispersion curves.

[ $CO_2$ ] sequestration in oil reservoirs can be one of the most effective strategies for long-term removal of greenhouse gas from atmosphere. This paper presents an advantage of the localized nonlinear approximation of integral equation solutions for inverting crosswell electromagnetic data, which are observed as a part of pilot project of $CO_2$ flooding at the Lost Hills oil field in central California, U.S.A. To monitor the migration of $CO_2$, we have used 2-D cylindrically symmetric and 2.5-D tomographic inversion methods. These two schemes produce nearly the same images if the borehole separation is large compared with the skin depth. However, since the borehole separation is much less than five skin depths in this $CO_2$ injection experiment, the 2.5-D model seems to be more reliable than the 2-D model. In fact, the pre-injection 2.5-D image is more successfully compared with induction logs observed in the two wells than the 2-D model. From the time-lapse crosswell imaging, we can confirm the replacement of brine with $CO_2$ makes a decrease of conductivity.

Evaluating stability of dam foundations is one of the prime areas of any rock engineering investigations. Despite best engineering efforts in the design and construction of dam foundations, the foundation regime of a constructed dam suffers deterioration due to continuous erosion from backwater current of dam discharge and dynamic effects of loading and unloading process. Even during construction, development of frequent cracks due to sudden thermal cooling of concrete blocks is not uncommon. This paper presents two case studies from India and Bhutan. In the first case, the back current of water discharge from the Srisailam dam in India had continuously eroded the apron and has eaten into the dam foundation. In the second case with dam construction at Tala Hydroelectric Project in Bhutan, sudden overflow of river during the construction stage of dam had led to development of three major cracks across the dam blocks. This was ascribed to adiabatic cooling effect of concrete blocks overlain by chilled water flow. Non-destructive evaluation of rock mass condition in the defect regime by the borehole GPR survey helped in arriving at the crux so as to formulate appropriate restoration plan.

The estimation of seawater intrusion into deep aquifers has been becoming an important subject in terms of site characterization for geological disposal of radioactive waste. Conventional direct-current resistivity methods have been used for ground water explorations and recently have been applied to environmental problems. However, electromagnetic methods are more practical and useful for such a deep investigation. We consider audio-frequency magnetotelluric (AMT) and surface-to-borehole electromagnetic (EM) tomography methods as promising tools for the investigation of deep aquifer. These methods were tested in the Hasunuma area, Chiba Prefecture, Japan. Although the study area is in an urban area, high-quality AMT data were acquired, which was mainly accomplished by night-time data recording and remote-reference data processing. One-dimensional inversion results of the AMT data revealed two extremely conductive zones, which is consistent with the electrical conductivity profile of pore water in core samples. It can be interpreted as the seawater intrusions into both zones. However, the chemical analysis of the groundwater sampled in the deep zone suggests that this groundwater must be fossil seawater that had been confined during sedimentation processes. In addition, the permeability coefficient of the deep layer is very low. Thus the deep conductive zone corresponds to the fossil seawater regarded as being difficult to flow.

This study presents a practical procedure for the Bayesian inversion of geophysical data by Markov chain Monte Carlo (MCMC) sampling and geostatistics. We have applied geostatistical techniques for the acquisition of prior model information, and then the MCMC method was adopted to infer the characteristics of the marginal distributions of model parameters. For the Bayesian inversion of dipole-dipole array resistivity data, we have used the indicator kriging and simulation techniques to generate cumulative density functions from Schlumberger array resistivity data and well logging data, and obtained prior information by cokriging and simulations from covariogram models. The indicator approach makes it possible to incorporate non-parametric information into the probabilistic density function. We have also adopted the MCMC approach, based on Gibbs sampling, to examine the characteristics of a posteriori probability density function and the marginal distribution of each parameter. This approach provides an effective way to treat Bayesian inversion of geophysical data and reduce the non-uniqueness by incorporating various prior information.

Hydrogeological survey and geochemical analysis were carried out in Phra Yun area, Northeast Thailand, which is a typical salt-affected area for an understanding of hydrogeological groundwater behaviours. Geological survey reveals the presence of G1 and F1 faults. Electromagnetic ground conductivity prospecting shows that the high conductivity zones of 15 mS/cm or more are distributed at underground of the G1 and F1 faults where saline groundwater is discharged. The distribution patterns of tritium concentration show that high tritium concentration zones of groundwater were recharged from pond and river. On the assumption that the annual average tritium concentration of precipitation in Northeast Thailand is same as tritium concentration of precipitation in Tokyo and groundwater flows as piston flow, the age of recharging precipitation of groundwater with 15 TU in 1997 could be estimated at 1967-1970 years. The velocity of groundwater flow was calculated to be $5.3{\times}10^{-7}\;m/s\;and\;2.1{\times}x10^{-6}\;m/s$ respectively from a duration time of 30 years and distance of groundwater flow 500m -2000m from the pond and river to the investigation wells. Because the estimated values of velocity of groundwater flow are compatible with the hydraulic conductivities, it is considered that 30 years is a reasonable period for recharging groundwater.

We have evaluated the extent of saltwater intrusion from electrical resistivity distribution in a coastal aquifer system in the southeastern part of Busan, Korea. This aquifer system is divided into four layers according to the hydrogeologic characteristics and the horizontal extent of intruded saltwater is determined at each layer through the geostatistical interpretation of electrical resistivity data. In order to define the statistical structure of electrical resistivity data, variogram analysis is carried out to obtain best generalized covariance models. IRF-k (intrinsic random function of order k) kriging is performed with covariance models to produce the plane of spatial mean resistivities. The kriged estimates are evaluated by cross validation to show a good agreement with the true values and the statistics of cross validation represented low errors for the estimates. In the resistivity contour maps more than 5 m below the surface, we can see a dominant direction of saltwater intrusion beginning from the east side. The area of saltwater intrusion increases with depth. The northeast side has low resistivities less than 5 ohm-m due to the presence of saline water in the depth range of 20 m through 70 m. These results show that the application of geostatistical technique to electrical resistivity data is useful for assessing saltwater intrusion in a coastal aquifer system.

A infiltration experiment of river water has been conducted to evaluate the applicability of electrical resistivity monitoring methods in an area containing gravelly deposits in Nagaoka, Japan. Apparent resistivity data, which are inverted to obtain the resistivity distribution, are measured with a newly developed system. This system can collect 490 data in an hour and be controlled with PC to store the data. Subsurface resistivity sections, which are obtained from two-dimensional nonlinear inversion of time-lapse apparent resistivity data, enable us to estimate the direction of the flow and the rate of infiltration. The infiltration rate is estimated to be $4.4{\times}10^4m/s$ in the early stage of the experiment when the infiltration process is dominant.

Continuous seismic refraction, reflection and echo-sounder surveys conducted at Koyna Project site provided geotechnical information which helped in choosing the alignment for Head race tunnel and in designing and choosing the site for Lake Tap. Seismic refraction survey both on land and in shallow water determined depths to bedrock and helped in inferring the bedrock quality. Seismic reflection survey mapped the subsurface stratigraphy with high resolution. Reservoir-bed and bedrock contours drawn from the results of the survey helped in choosing the tunnel alignment and the lake tap position cost effectively. It was inferred from the results of the survey that the geology and the quality of rock do not change unexpectedly around the site for extension of Head race tunnel and the lake tapping. The bedrock levels evaluated by seismic survey agreed remarkably well with those inferred in boreholes having Rock Quality Designation 90 percent or more.

In gravity data correction process, mass effect of the upper part of base level is removed with Bouguer density. Usually, Bouguer density is estimated as a mean density in the field area. But, this may causes a serious problem when ore body is in the area. To overcome this problem, we tried to apply a new method mixing up mass corrections and inversion (3DGTI). 3-D Gravity Terrain Inversion (3DGTI) includes information of topography and distribution of Bouguer density. For this method does not remove the mass effect above base level, it is no longer useless to use Bouguer density. Numerical model tests have shown that the 3DGIT successfully retrieves the anomalous subsurface density distribution of both surface and deeper layers. Model tests shows that this method shows better results than those of conventional one, especially when main target is ore body. The inversion result well delineates the three-dimensional shape of the intruded granite body and basement.

Clay minerals show a distinct induced-polarization phenomenon, which is one of the most important factors for predicting groundwater flow and contaminant transport. This paper presents a step-by-step process to estimate Cole-Cole parameters from spectral induced-polarization (IP) data measured on the surface of three-dimensional earth. First, the inversion of low-frequency resistivity survey data is made to identify the dc resistivity ${\rho}_dc$ of a volume having IP effects. The other parameters, chargeability m, time constant $\tau$, and frequency dependence c, are sought for the polarizable volume. Next, using multi-frequency data, c can be obtained as high or low asymptotes of the slope of log phase vs. log frequency. Further, for low m, intrinsic $\tau$ is approximated by apparent one, ${\tau}_a$, which is derived from the relation ${{\omega}{\tau}}_a$=1 at an angular frequency $\omega$, where the imaginary component of spectral IP data has an extreme value. Finally, to obtain intrinsic m a two-step linearized procedure has been derived. For a body of given $\tau$ and c, forward modeling with a progression of m values yields a plot of observed vs. intrinsic imaginary components for a frequency. Since this plot is essentially linear, to extract the intrinsic imaginary component is quite simple with an observed value. Using the plot of intrinsic imaginary component vs. m, intrinsic m is determined. We present a synthetic example to illustrate that the Cole-Cole parameters can be recovered from spectral IP data.

A three-dimensional (3D) inversion technique has been developed for interpretation of magnetotelluric (MT) data. The inversion method is based on the linearized least-squares (Gauss-Newton) method with smoothness regularization. In addition to the underground 3D resistivity distribution, static shifts are also treated as unknown parameters in the inversion. The forward modeling is by the staggered-grid finite difference method. A Bayesian criterion ABle is applied to search the optimum trade-off among the minimization of the data misfit, model roughness and static shifts. The method has been applied to several MT datasets obtained at geothermal fields in Japan and other Asian countries. In this paper, two examples will be discussed: one is the data at the Ogiri geothermal area, southwestern Japan, and the other is at the Pohang low-enthalpy geothermal field, southeastern Korea. The inversion of the Ogiri data has been performed stably, resulting in a good fitting between the observed and computed apparent resistivities and phases. The recovered 3D resistivity structure is generally similar to the two-dimensional (2D) inversion models, although the deeper portion of the 3D model seems to be more realistic than that of the 2D model. The 3D model is also in a good agreement with the geological model of the geothermal reservoirs. 3D interpretation of the Pohang MT data is still preliminary. Although the fitting to the observed data is very good, the preliminary 3D model is not reliable enough because the station coverage is not sufficient for a 3D inversion.

This paper deals with the recent achievements in AMT data acquisition and processing, which have improved the quality of the collected data,. Comparing AMT with VES, IP and logging data, it is concluded through some case histories that qualified results can be achieved in probing the karst water with AMT technique.

Inversion of traveltime requires an efficient algorithm for computing the traveltime as well as its $Frech\hat{e}t$ derivative. We compute the traveltime of the head waves using the damped wave solution in the Laplace domain and then present a new algorithm for calculating the $Frech\hat{e}t$ derivative of the head wave traveltimes by exploiting the numerical structure of the finite element method, the modem sparse matrix technology, and SWEET algorithm developed recently. Then, we use a properly regularized steepest descent method to invert the traveltime of the Marmousi-2 model. Through our numerical tests, we will demonstrate that the refraction tomography with large aperture data can be used to construct the initial velocity model for the prestack depth migration.

Ground subsidence has occurred in the downtown of Muan-eup in Korea. Integrated geophysical survey, including two-dimensional resistivity, CSMT(Controlled source magnetotelluric), magnetic, borehole logging, GPR and resistivity tomography, has been conducted to investigate the cause of subsidence and ground conditions. Since the target area is in the city downtown, there were no spaces for surface geophysical methods. To get regional geology and to facilitate the detailed geophysical interpretation in the survey area, two-dimensional resistivity, CSMT and magnetic surveys have been applied in the outer region of the downtown. From these results, we could accurately define the Gwangju fault system and estimate the geologic conditions in the downtown. For the detailed survey of the downtown area, resistivity tomography and borehole logging data have been acquired using a few tens of densely located boreholes. Among these survey results, borehole logging data provided the guide to classification of the rock type and we could define the geologic boundary of granite and limestone formations. From the resistivity tomograms of 42 sections, which are densely located enough to be interpreted in a three-dimensional manner, we could delineate the possible weak zones or cavities in the limestone formations. In particular, resistivity tomograms in the subsided area showed the real image of ground subsidence. The map of hazardous zone has been derived from the joint interpretation of these survey results and we could provide the possible reinforcement strategy in this area.

To develop an easy, low-cost method for evaluating the degree of weathering by radioactive prospecting, radioactive prospecting and the investigation of the degree of weathering were carried out in the southern Kitakami massif of Iwate Prefecture, Japan, in which weathering granitic rocks was distributed. Fifty outcrops in the study area were selected, and strength of the gamma-rays emitted from the weathering bedrock of $^{40}K,\;^{214}Bi,\;and\;^{208}Tl$ was measured for 15 minutes at each point. At the same points, soil hardness was measured on the surface of the outcrop with a Yamanaka soil penetration tester. In addition, 100cc samples of each outcrop were taken with the sampler. The samples were analyzed by XRD, and the kind of the rock-forming minerals containing K was identified. We then compared the degree of weathering and the radioactive prospecting results by using K as an indicator. The relation between $^{40}K/^{208}Tl$ gamma rays counting rate by the radioactive prospecting and the hardness index showed a positive correlation as a result of the investigation, and the correlation coefficient ($R^2$) was 0.67.Moreover, when $^{40}K/^{208}Tl$ gamma rays counting rate emitted from the bedrock was low, the number of rock-forming mineral species containing K was also low. Thus, it was found that $^{40}K/^{208}Tl$ gamma rays counting rate measured by the radioactive prospecting could be used as an indicator of the degree of weathering.

Helicopter-borne electromagnetic (HEM) systems were originally developed for the exploration of mineral deposits. The frequency range of a conventional HEM system for mineral exploration, however, is relatively low and so not invariably suitable for its application to the fields of civil engineering because of its poor resolution in the shallower part of the earth. A DIGHEM HEM system was acquired by Nippon Engineering, with the frequencies chosen by the senior author. The five frequencies range from 220 Hz (the lowest) to 137,500 Hz (the highest). These frequencies improve the resolution of materials in the shallower part while maintaining a depth of investigation of greater than 100 m. This paper describes six case histories of geological and geotechnical surveys for civil engineering using HEM. These case histories include HEM surveys for investigating landslide, an alluvial area, root selection of road construction, areas related to dam and tunnel construction, and the simultaneous joint inversion of HEM and CSAMT data for a deep tunnel. These survey results show that HEM has sufficient resolution in both horizontal and vertical directins to contribute significantly to outlining the regional geology and its engineering problems.

Geophysical well logging techniques which are useful for delineating permeability of geological formation have been reviewed. A new technique for obtaining permeability using conductivity log technique has been discussed. This conductivity logging technique has been tested by monitoring the conductivity change within the model hole using borehole environment water and incoming-outgoing water of different salinity with constant flow rate by maintaining balance between inflow and outflow. Conductivity variation features depended mainly on flow rate, density contrasts due to salinity and temperature contrasts between fluid within the hole and incoming-outgoing fluid. The results of the experiment show uniform change of fluid conductivity within bore hole with time, and a fairly good correlation between the flow rate and the conductivity change rate. This conductivity logging technique is expected to be an efficient tool for determining permeability.

As we know, roof is composed of heterogeneous rock. When roof fractures, a large amount of energy would be released in the form of seismic wave. How to identify the abnormal signal of seismic wave is a much difficult problem, there are many methods used usually, such as Fourier Transformation, filter technique etc., but abnormal signal can't be recognized accurately. In this paper, multi-resolution wavelet technique is used to identify the first and second variation point, based on the Lipschitz $\alpha$. A living example analysis shows, multi-resolution wavelet technique can identify the abnormal signal of seismic wave effectively in different scale, and the omen of roof fall can be grasped in order to forecast the roof fall accurately. It provides a new idea for the predication of catastrophe on rock mechanics and engineering.

Refraction tomography requires an algorithm for efficiently computing the traveltimes and their $Fr\'{e}chet$ derivatives. We have attempted to solve the damped wave field using the frequency domain finite element model ing and then invoked the reciprocity theorem to calculate the $Fr\'{e}chet$ derivative of the traveltime with respect to the subsurface parameter. Then, we used a damped least square method to invert the traveltimes of the Marmousi 2 model. Numerical tests demonstrate that the refraction tomography with large aperture data can be used to estimate the smooth velocity model for the prestack depth migration.

The accurate mapping of the basement is one of the most crucial factors in construction of harbour facilities and bridges in the coastal areas. In shallow waters, the seismic reflection method often fails to image the basement geometry beneath the sediment cover in many cases. We present the result of a shallow marine seismic refraction survey using two ships, l2-channel hydrophone arrays deployed on the bottom and a borehole sparker or percussion powder as sources. Velocity structure could be computed by tomography algorithm since more than 6 different source points had been applied for one spread. The comparison of the results of the refraction survey with drilling logs demonstrates remarkable consistency in basement geometry. It thus appears that the refraction method in this study is an efficient and cost-effective way to investigate the basement structure in coastal area, river, and lake.

A combination of drilling, hydrogeochemical survey, geophysical survey and the numerical modelling for the flow and transport of groundwater was performed to evaluate the seawater intrusion in Baeksu-eup, Yeonggwang-gun, Korea. The survey area extends to over 24 $km^2$. Twelve wells were also drilled for the collection of geologic, geochemical, hydrologic, and geophysical logging data to delineate the degree and vertical extent of seawater intrusion. To evaluate and map the salinity in a coastal aquifer, geophysical data and hydrogeochemical results were used. Layer parameters derived from VES data, various in situ physical properties from geophysical well loggings, and the estimated equivalent NaCl concentration were used as the useful input parameters for the numerical simulation with density-dependent flow. Our multidisciplinary approach for evaluating the seawater intrusion can be considered as a valuable attempt to enhancing the utilization of various data and the reliability of numerical ground modelling.

Waveform inversion requires extracting a reliable low frequency content of seismic data for estimating of the low wave number velocity model. The low frequency content of the seismic data is usually discarded or neglected because of the band-limited response of the source and the receivers. In this study, however small the spectral of the low frequency seismic data is, we assume that it is possible to extract a reliable phase information of the low frequency from the seismic data and use it in waveform inversion. To this end, we exploit the frequency domain finite element modeling and source-receiver reciprocity to calculate the $Frech\`{e}t$ derivative of the phase of the seismic data with respect to the earth model parameter such as velocity, and then apply a damped least squares method to invert the phase of the seismic data. Through numerical example, we will attempt to demonstrate the feasibility of our method in estimating the correct velocity model for prestack depth migration.

The traditional and still most widely used, test methods for concrete structures are destructive method, such as coring, drilling or otherwise removing part of the structure to permit visual inspection of the interior. While these methods are highly reliable, they are also time consuming and expensive, and the defects they leave behind often become focal point for deterioration. In this study, tomography by theoretical inversion method in case of elastic wave using impact-echo method among concrete non-destruction test method was made. Taken model experiments are theoretical inversion method and time domain and frequency domain test on pier test model at laboratory level. Also experiment concerning frequency domain on 3 kinds of tunnel model with I-dimension form was carried out.

KIGAM (Korea Institute of Geoscience and Mineral Resources) launched a new project to develop the low-enthalpy geothermal water in the area showing high geothermal anomaly, north of Pohang city, for large-scale space heating from KORP (Korea Research Council of Public Science & Technology) funding. Surface geologic and geophysical surveys including Landsat TM image analysis, gravity, magnetic, Magnetotelluric (MT) and controlled-source audio-frequency MT (CSAMT) and self-potential (SP) methods have been conducted and the possible fracture zone was found that would serve as deeply connected geothermal water conduit. By the end of 2003, two test wells of 1 km depth will be drilled and various kinds of borehole survey along with additional MT measurements and sample analysis will follow and then the detailed subsurface condition is to be characterized. Next step would be drilling the production well of 2 km depth and all further steps remain to be determined depending upon the results of the test well studies.

Wheat, barley, rye, and maize were grown in field and pot experiments at various non-contaminated soils in order to establish uptake rates for added selenate, and to find baseline concentrations for various soil types. Edible parts (grains) and stalks of the crops were analyzed separately for Se, as well as for Ca, Cu, Fe, Mn, P, S, and Zn. The addition of Na-selenate in admixture with the NPK 20:8:8 fertilizer had no influence on the composition of the other elements investigated. The proportions of added nitrate: selenate, and sulfate:selenate were kept constant. The Se- uptake rate differed among the cereals tested, it was highest for winter wheat. Utilization of added Se in the field ranged from $0,4-4,7\%$, and and in the pots from $3,3-5,4\%$, it was markedly lower in clay soil. Whereas P and Zn were preferably found in the grains, Ca-Fe-Mn-S got enriched in the stalks. For the fields, the location had some influence upon Fe, Mn, and Zn, whereas it was not important for P, S, Cu, and strikingly, Ca. Pot and field experiments on similar soils led to different results, except for P and S. Maize (whole grains) was significantly lower in Ca, Cu, and Mn, and might even cause trace element deficiencies, if exclusively fed. Few correlations between the trace elements investigated led to the conclusion that most element contents were governed by plant metabolism. Variations of mobile Fe in the soils were balanced by uptake into the stalks. The data are compared with data from other presumably non-contaminated sites.

The effect of reductive treatment with elemental iron on the extent of mineralization by Fenton oxidation was studied for the explosive 2,4,6-trinitrotoluene (TNT) and hexahydro-l,3,5-trinitro-l,3,5-triazine (RDX) using a completely-stirred tank reactor (CSTR). The results support the hypothesis that TNT and RDX are reduced with elemental iron to products that are oxidized more rapidly and completely by Fenton's reagent. Iron pretreatment enhanced the extent of TOC removal by approximately $20\%\;and\;60\%$ for TNT and RDX, respectively. Complete TOC removal was achieved for TNT and RDX solutions with iron pretreatment under optimal conditions. On the other hand, without iron pretreatment, complete mineralization of TNT and RDX solutions were not achieved even with much higher $H_2O_2$ and $Fe^{2+}$ concentrations. The bench-scale iron treatment-Fenton oxidation integrated system showed more than $95\%$ TOC removal for TNT and RDX solutions under optimal conditions. The proposed zero-valent iron-Fenton process was evaluated with pink water from the Iowa Army ammunition plant. Results from batch and column experiments show that TNT, RDX, and octahydro-l,3,5,7-tetranitro-l,3,5,7-tetrazocine (HMX) were completely removed from the pink water and that triaminotoluene (TAT) and ${NH_4}^+$ were recovered as products in reduction with zero-valent iron. By using an integrated system, $83.3\pm4.2\%$ of TOC was removed in a CSTR with 10 mM of $Fe^{2+}$ and 50 mM of $H_2O_2$. These results suggest that the reduction products of TNT and RDX are more rapidly and completely mineralized by Fenton oxidation and that a sequential iron treatment-Fenton oxidation process may be a viable technology for pink water treatment.

Soil-gas radon level and other atmospheric factors have been measured at residual soil profiles that overlie granite bedrock which consists of major geology in Korea for 6 months from November, 2000 to April, 2001. Seasonal variations of soil-gas radon concentration are generally of greater magnitude than day-to-day fluctuations. The highest radon concentrations of 5,131 pCi/L measured during winter season and the lowest radon concentrations of 107 pCi/L during spring season. Two study areas, Bongcheon-dong(granite bedrock) and Seongnam-Yongin(gneiss bedrock) were investigated to assess the radon potential according to their field survey and emanation tests. The mean values of radon decrease in sequentially from Suji-A(813 pCi/L)>Suji-B(757 pCi/L)>Bundang-B(691 pCi/L)>Bundang-A(643 pCi/L)>Bongcheon-dong(513 pCi/L). Estimated soil-gas radon potential using maximum radon emanation ratios of each study area decreases in the order of Bongcheondong(950 pCi/L)>Suji-B(524 pCi/L)>Bundang-A(437 pCi/)>Bundang-B(259 pCi/L)>Suji-A(230 pCi/L) areas. The values of indoor radon and its daughter product concentrations in Bongcheon-dong area show that indoor basement rooms in poor ventilation condition could be classified as extremely high radon risk location of more than 4 pCi/L Rn and 0.02 WL.

This study has been carried out to evaluate the passive treatment systems for acid mine drainage in Korea and to suggest, if possible, the method for the improvement. 35 passive treatment systems in 27 mines have been constructed since 1996. SAPS, being the main process, was combined with more than one of processes such as anaerobic wetland, aerobic wetland, and oxidation pond for the construction of passive treatment system. Problems observed during the operation include the poor sulfate removal ratio, overflow, leakage, unusabless of the whole system, and inefficiency. The reasons of the poor sulfate removal ratio are believed that the low temperature during the winter prohibits the SRB activity and HRT for bacterial sulfate reduction is insufficient. An alternative method In Adit Sulfate Reducing System which enables to keep the temperature constant at about $15^{\circ}C$ was suggested. IASRS is the methods of placing the SAPS inside the adit, which enables the temperature around the system constant can be maintained. The experiments using the laboratory scaled model systems made up of four sections showed high efficiencies in pH control and metal removal ratios, but showed still low sulfate removal ratio of about $23\%$ also with high COD at the beginning of the operation.

The objective of this study is to investigate the vertical variation of heavy metals in tailings from the Sangdong W mine. Tailings samples were taken at 6 drilling sites with 50cm intervals up to 21 meters in depth and dried at room temperature. The pH value, loss-on-ignition and water contents were measured. In addition, chemical compositions of the samples were determined by AAS after 0.1 N HCl leaching and ICP-AES after aqua regia leaching. The pH values were in the range of 7.4 to 9.5 due to chemical reactions of carbonate minerals. The ranges of heavy metals (mg/kg) extracted by 0.1N HCl were from 0.17 to 0.93 for Cd, 0.04 to 4.39 for Cu, 0.03 to 10.9 for Pb and 0.06 to 14.1 for Zn and those extracted by aqua regia were $3.10\~10.5,\;23.61\~251,\;63.7\~337\;and\;42.6\~134$ for Cd, Cu, Pb and Zn, respectively. Generally, the metal concentrations in tailings extracted by 0.1N HCl decreased with depth, whilst those extracted by aqua regia have a tendency to increase with depth in some case. Those trends might be due to the change of oxidation-reduction condition of the tailings.

In order to assess the risk of adverse health effects on human exposure to arsenic and heavy metals influenced by past mining activities, environmental geochemical surveys were undertaken in the abandoned metal mine areas (Dongil Au-Ag-Cu-Zn, Okdong Cu-Pb-Zn, Songcheon Au-Ag, Dongjung Au-Ag-Pb-Zn, Dokok Au-Ag-Cu and Hwacheon Au-Ag-Pb-Zn mines). Arsenic and other heavy metals were highly elevated in the tailings from the Dongil, the Songcheon and the Dongjung mines. High concentrations of heavy metals except As were also found in tailings from the Okdong, the Dokok and the Hwacheon mines. These significant concentrations can impact on soils and waters around the tailing dumps. Risk compounds deriving from mine sites either constitute a toxic risk or a carcinogenic risk. The hazard index (H.I.) of As in the Dongil, the Okdong, the Songcheon and the Hwacheon mine areas was higher value more than 1.0. In the Okdong and the Songcheon mine areas, H.I. value of Cd exceeded 1.0. These values of As and Cd were the highest in the Songcheon mine area. Therefore, toxic risks for As and Cd exist via exposure (ingestion) of contaminated soil, groundwater and rice grain in these mine areas. The cancer risk for As in stream or ground water used for drinking water from the Songcheon, the Dongil, the Okdong, the Dongjung and the Hwacheon mine areas was 3E-3, 8E-4, 7E-4, 2E-4 and 1E-4, respectively.

LCA (Life Cycle Assessment), that unifies the scale of various environmental impacts, and simulated annealing are applied to optimizing electrolysis of wastewater from PCB (Printed Circuit Board) production. The changes of environmental impact can be quantified with LCA and the total changes of environmental impacts can be expressed as a function of power consumed, Cu recycled, $Cl_2$, NOx and SOx discharged through restriction of feasible reactions. In a single-variate condition, the environmental optimum can be easily obtained through plotting and comparing each environmental impact value. In 8V potentiostatic electrolysis, the lowest environmental impact can be achieved after 90min. To optimize a multi-variate conditional system, simulated annealing can be applied and this can give the quick and near optimum in complex systems, where many input and output materials are involved, through experimentally measured values without a theoretical modeling.

Mine waters, surface waters and groundwaters were sampled around seven Au-Ag mine areas (Dongil, Okdong, Dongjung, Songcheon, Ssangjeon, Dogok and Gubong Au-Ag mines). The main contamination sources of As in these abandoned Au-Ag mines can be suggested as mine tailings and waste rocks including the sulfide gangue minerals (arsenopyrite). The relatively high concentration of As in mine waters was shown in the Dongil (524 ${\mu}g/L$) and the Dogok (56 ${\mu}g/L$) mine areas. Arsenic concentrations in stream waters from the Dongil ($0.9\~118{\mu}g/L$), the Songchon ($0.8\~63{\mu}g/L$), the Ssangjeon ($1.6\~109{\mu}g/L$) and the Gubong ($3.6\~63{\mu}g/L$) mine areas exceeded the permissible level for stream water in Korea. Groundwaters collected from the Dongil ($0.9\~64{\mu}g/L$ ), the Okdong ($0.2\~69{\mu}g/L$) and the Gubong ($0.5\~101{\mu}g/L$) mine areas contained high As concentration to cause the arsenicosis in these areas. In As speciation, the concentration ratios of As(III) to As(total) present up to $75\%$ and $100\%$ in stream waters from the Okdong and the Songcheon mines, and $70\%$ in groundwaters from the Okdong and the Dongjung mines. Arsenic concentration decreases downstream from the tailing dump correlatively with pH and Fe concentration. Highly elevated As concentrations are found in the dry season (such as April and March) than in the wet season (September) due to the dilution effect by heavy rain during summer in stream waters from the Dongil and the Songcheon mine areas.

There are some reports which are about measuring the discontinuity from 3D model of rock surface. To generate 3D model is mainly based on remote sensing technique like laser scanning and digital stereo photogrammetry. It is obtained the DEM of the rock slope using above techniques in this study, and examined a suitability and improvement of the photogrammetry for the rock slope by overlap the DEM. It seems that accuracy of DEM generated by photogrammetry technique is influenced by the geometry of rock and type of joint.

In order to investigate the level of heavy metal contamination and the seasonal variation of metal concentrations in soils and sediments influenced by past mining activities, tailings, soil and sediment samples were collected from the Hwacheon mine in Korea. The main pollution sources in this mine site are suggested as tailings and mine waste rocks. Elevated levels of Cd, Pb and Zn were found in soils and sediments. In a study of seasonal variation on the heavy metals in soils and sediments, heavy metals were higher enriched collected from before rainy season ($2^{nd}$ sampling) than after rainy season ($1^{st}$ sampling). Also, in order to estimate the microbial effects on Cd speciation in sediments, bacteria which can adsorb Cd was isolated and Cd adsorption characteristics of isolated bacteria in Cd solution was evaluated. The Cd bioremoval efficiency in Cd solution (5 ppm) by bacteria was more than $90\%$. Bioremoval efficiency in single metal solution was higher than that in mixed metal solution of Pb and Zn.

Indigenous bacterial mediation of As in contaminated sediment after biostimulation with a variety of carbon sources was investigated under aerobic and anaerobic conditions. Under the aerobic condition with lactate supply, indigenous bacteria increased the amount of total As extracted from the sediment and most dissolved As existed as As(V). Under the anaerobic, glucose-supplied condition, dissolved As diminished with time likely due to production of As sulfide(s) and subsequent precipitation, which resulted from bacterial reduction of ${SO_4}^{2-}$. The results implied that bacterial natural attenuation of As in subsurface has a potential to be practically applied.

The Dongjung Au-Ag-Cu mine area was seriously contaminated with As and heavy metals-Cd, Cu, Hg, Mn, Pb and Zn etc. Those elements were highly accumulated in plants grown at farmland as well as farmland soil. Stream waters and groundwater which has been used as drinking water around the mine site contain high levels of heavy metals, especially As. As a result of human health risk assessment using EHS(Extraction of Heavy metals in Stomach and Small intestine) test for bioaccessible contents of heavy metals, there is a potential of cancer and adverse effects on human health for the residents of the mine area.

There are many stone monuments which are weathered by natural or artificial factors in Korea. The partly weathering in stone could accelerate the rate of weathering, so it is demanded to keep them from the further weathering. Consolidation is evaluated as one of the efficient treatments which have a good effectiveness in stone monuments. But following the former researches, the effectiveness of consolidation could be different, related to the kinds of stone or consolidants. Therefore, in this study, the change of properties was monitored in the long term for the exact evaluation of the stability of the consolidated stone. It is estimated that the pore structure of consolidated stone was filled by consolidant, according to the measurements of ultrasonic velocity, and the digital image analysis of the sample was shown that the color property of stone surface has varied during the curing.

The objective of this study is to investigate the contamination levels and dispersion patterns of arsenic and heavy metals and to estimate the bioaccessible fraction of the metals in soil and plant samples in the vicinity of the abandoned Songcheon Au-Ag mine. Tailings, soils, plants (Chinese cabbage, red pepper, soybean, radish, sesame leaves, green onion, lettuce, potato leaves, angelica and groundsel) and waters were collected around the mine site. After appropriate preparation, all samples were analyzed for As, Cd, Cu, Pb and Zn by ICP-AES and ICP-MS. Elevated levels of As and heavy metals were found in tailings. Mean concentrations of As in agricultural soils were higher than the permissible level. Especially, maximum level of As in farmland soil was 513 mg/kg. The highest concentrations of As and Zn were found in Chinese cabbage (6.7 mg/kg and 359 mg/kg, respectively). Concentrations of As, Cd, and Zn in most stream waters which are used for drinking water around this mine area were higher than the permissible levels regulated in Korea. Maximum levels of As, Cd and Zn in stream waters were 0.78 mg/L, 0.19 mg/L and 5.4 mg/L, respectively. These results indicate that mine tailings can be the main contamination sources of As and heavy metals in the soil-water system in the mine area. The average of estimated bioaccessible fraction of As in farmland soils were $3.7\%$ (in simulated stomach) and $10.8\%$ (in simulated small intestine). The highest value of bioaccessible fraction of metal in farmland soils was $46.5\%$ for Cd.

During the last few decades, the frequency and distribution of rockfalls have been increased in Korea due the development of rock slope in mountain areas. Although the scale of, falling block of rockfall was small, there were some casualties of lives and loss of properties. In order to reduce damage from rockfall, analysis on rock slope indanger of rockfall should be carried out. Thus, the simulation softwares for rockfall behavior analysis have been introduced. In this study, geotechnical investigation and input data for rockfall simulation are described with the comparison among four commercially available rockfall simulation softwares. Finally, rockfall simulation works are described by exmining a case of inaccessible rock-slope of potential rockfall in Korea.

The objective of this study is to implement geological functions on the Internet such as the visualization, search and analysis of geological data. To improve their performance and reduce the cost during development, special concerns should be given to the dynamic integration and interoperability of them. XML Web Service technologies can be adopted and OpenGIS specifications should be considered for the integration of distributed GIS components and their interoperability. Because standard Internet protocols and XML Messages are used in the Web Service, it makes the integration among different hardware and software systems possible. OpenGIS specifications also enable users to get benefits from geographic information and services across any network, platform and application. A prototype of efficient Internet GIS applications has been suggested in this study.

Some well-known artificial weathering tests such as freezing-thawing, acid immersion, and salt crystallization are adopted to examine the change of rock properties during the processes of artificial weathering. Granites and other rock types of limestone, marble and basalt collected from different quarries in south Korea were sampled for this study. All tests were performed up to 30 cycles and physical properties were measured after experiencing every ten cycles of artificial weathering tests. During the tests, the variation trends of rock properties were too variable to draw generalized variation patterns but it can be concluded that weathering agents have different effect on rock properties depending on weathering circumstance and time. Even in short terms of salt crystallization tests, some rocks were severely deformed and then burst, and in the early stages of salt weathering, recrystallized salts filling pores and cracks in rocks could be a important factor affecting rock properties.

Recently, the recycled amount of electric, mechanical parts, and appliances in artifacts has increased. These products use valuable rare metals such as platinum group metals and gold, which are included occasionally as additives. Rare metals are maldistributed in the world and most of them are produced in small quantities. A small amount of rare metals used in the appliances causes a large loss of rare metal resources because of the lack of an economically recycling method. The present recycling technologies including physical and chemical separation methods that are considered for recycling of electric, mechanical parts and appliances.

Titanium dioxide ($TiO_2$) nanoparticles were prepared by the oxidation of titanium tetrachloride ($TiCl_4$) in a diffusion flame reactor. The average diameter of particles was 15 to 30 nm and mass fraction of anatase ranged from $40\;to\;80\%$. Effects of particle size and phase composition of those $TiO_2$ nanoparticles on photocatalytic properties such as decomposition of methylene blue and bacteria gas were investigated. The degree of decomposition of methylene blue by the $TiO_2$ nanoparticles under the illumination of the black light was directly proportional to the anantase mass fraction, but inversely to the particle size. The decomposition of bacteria by the $TiO_2$ nanoparticles under the illumination of the fluorescent light showed the same trend as in the case of the methylene blue.

In Taiwan there are 21 Municipal Solid Waste Incinerators (MSWI) built to treat 80% of the MSW nationwide. Approximately 2,000 tons of incineration ashes of municipal waste contain reaction ash and fly ash (3:1 by weight)will be produced daily, and this may cause a serious waste problem. According to EPA regulations, reaction ash and fly ash produced after incineration should be properly treated. Landfill capacity barely meets the general demands. More efficient actions should be planned and taken. The study found 'reclamation' should be the optimal solution to this problem. Only limited research and previous successful experiences are available among other countries. An incinerator in Northern Taiwan is chosen for this study to make environmental bricks from the reaction ash and fly ash. From the previous tests, the results of strength test were measured. From the previous test results, the fly ash products have not reached the desired strength; hence, reaction ash is chosen for further pilot study. In the experiment, incineration ashes, cement and gravel are mixed in the ratio of 1:1:1(by weight), to ground concretization aggregate and pelletization aggregate, the concrete products made from the aggregates were of the strength of 108 $kgf/cm^2$ and 142 $kgf/cm^2$ individually. For the purpose of making nonstructural walls which met the State Building Standards. In the study, 50 tons of concrete products was yielded from aggregate and environmental bricks. Further observation and supervision are recommended to ascertain the resource recycling and reclamation. EPA has planned to build three 'Recycling Plants' in northern, middle and southern Taiwan to develop efficient techniques to produce concrete products, sub-base course, soundproofing wall, gravel, artificial fishing reefs, tiles, drainage, bricks and etc. This experiment of the demonstrative plant solves the problem of the incineration ashes and opens another opportunity to reclaim them.

The variable flotation response ores from different deposits results basically from mineralogical association and their differences. Development of new techniques for analyzing the metallurgical performance of flotation and other concentration processes is demanded even in the treatment of rather simple ores such as porphyry ores. Diagnostic metallurgical analysis can be used to quantify the most possible recovery processes. Several porphyry copper/gold ores around the world were used to examine the responses in flotation, gravity separation and cyanidation in order to define the linkage between the recovery processes for both copper and old values. Laboratory batch flotation, gravity separation and cyanidation tests were carried out on these samples. All results were used to correlate the relative recovery of copper and gold, and to predict the highest possible metal recovery in the system. The metallurgical predictions were made according to the flotation conditions used and gravity separation. The results of various concentration processes on each porphyry ore samples are presented and discussed. All seven samples have shown significantly different gold/copper metallurgy. The grade/recovery relationships of gold and copper in the laboratory batch tests for the best results and the plants are given in the Figures below. The results of laboratory tests show that the copper recoveries converged to about $90\%$, but the gold recoveries were spread over $55-80\%$, except the K S ore. Series of standard cyanidation tests on the flotation concentrate samples and gravity separation using Knelson Separator on heads ores were carried out to cross-link the metallurgy and mineralogy of gold in the porphyry ores.

Processing and smelting of copper containing sulphide concentrates result in the accumulation of impurities into various process streams. All primary copper smelters and refineries around the world produce significant amounts of slag, dust, sludge, residues and others, which contain copper and precious metals. The recovery of these valuable metals is essential to the overall economics of the smelting process. Physical, chemical and mineralogical characterization of particular slag and Cottrell dusts from primary smelters and $Dor\'{e}$ furnace (TBRC) slag and Pressure Leached Anode slimes from a copper refinery have been carried out to understand the basic behind the recovery processes. Various process options have been evaluated and adapted for the treatment of slag from different smelting furnaces and Cottrell dusts as well as the intermediate products from copper refineries. Besides the hydro- or pyro-metallurgical treatments, the above mentioned physical separation options such as magnetic, gravity separation, flotation and precipitation flotation processes have been successfully identified and adapted as the possible process options to produce a Cu-rich or precious metal-rich concentrates for in-house recycling and other valued by-product for further treatment. The results of laboratory, pilot plant and production operations are presented, and incorporation of several alternative flowsheet is discussed in this paper.

The general red-mud minerals consist of hematite, sodalite, anatase, quartz, gibbsite and miner impurities. This gives serious environmental damage for the ocean disposal. It mixed with chloride compound and the content of chlorine is about 2,000-3,000ppm. This paper can be suggested the chloride reduction technology that is applied basically mineral processing by physical separation. Then it can be possible to produce the totally 24wt. $\%$ useful red-mud which chloride content is less then 400ppm.

The influence of Poisson's ratio on the tensile strength of brittle materials is neglected in many studies. When brittle materials are loaded in compression or impact, substantial tensile stresses are induced within the materials. These tensile stresses are responsible for splitting failure of the materials. In this paper, the state of stress in a spherical particle due to two diametrically opposed forces is analyzed theoretically. A simple equation for the state of stress at the center of the particle is obtained. An analysis of the distribution of stresses along the z-axis due to distributed pressures and concentrated forces, and on diametrically horizontal plane due to concentrated forces, shows that it is reasonable to propose the tensile stress at the center of the particle at the point of failure as a tensile strength of the particle. Moreover, the tensile strength is a function of the Poisson's ratio of the material. As the state of stress along the z-axis in an irregular specimen tends to be similar to that in a spherical particle compressed diametrically with the same force, this tensile strength has some validity for irregular particles as well. Therefore, it can be proposed as the tensile strength for brittle materials generally. The effect of Poisson's ratio on the tensile strength is discussed.

This study was performed to develop the method for producing industrial coal sources by cleaning Korean anthracite. Laboratory hindered-settling separation column was set and three coal samples were used for tests. Tests were conducted to evaluate the effects of the major operating variables, teeter water flow rate and relative column pressure (set point). Additional tests were performed to elevate the yield and properties of the products using air bubble injecting process. In results, nice products were obtained with high teeter water flow rate and air bubble injection. Also, model of continuous hindered-settling separation process was established to assist the evaluation of the equipment and several operating variables, such as dispersion, teeter water flow rate, feeding rate, etc.

Solvent extraction experiments for separation of impurities from Ni-rich solution were carried out for manufacturing of high purity Ni compounds from acid leaching solution of spent Ni-Cd secondary battery. Artificial and leaching solutions were used as aqueous phases and PC88A saponified by sodium in kerosene were used as organic phase. The extraction order is Zn>Mn>Co>Ni and extraction percentage of metal ions was increased with increase of the concentration of extractant, initial pH of aqueous phase and ratio of O/A. The separation of cobalt, zinc and manganese from nickel was effectively accomplished at the condition of extraction stage=l, O/A=1 and initial pH 5.0 with 1.0 $mol/dm^3$ PC88A saponified to $50\%$ with NaOH.

Lanthanum oxyfluoride can be synthesized by mechanochemical (MC) reaction between lanthanum oxide ($La_2O_3$) and polytetrafluoroethylene (PTFE, ($({CF_2CF_2}_n)$) in air using a planetary mill. MC reaction between the two materials induced from intensive grinding operation. The MC reaction is almost finished by 240min, and the products ground for 240min or more are composed of LaOF, amorphous $La(CO_3)F$ and amorphous carbon (C). Heating this MC reaction products at $600^{\circ}C$ enables us to eliminate amorphous C and decompose $La(CO_3)F$ into LaOF, so that pure LaOF material can be obtained as the final product. The average particle size of the final product (purified LaOF) is around few ten nanometers.

The retardation of heavy metals in a mixture of fly ash and bentonite was studied as a potential barrier material for a landfill. Column tests were conducted using synthetic leachate having 100 mg/L and 50 mg/L of lead (Pb) and cadmium (Cd), respectively. Results indicated that the mixture had obvious retardation ability for heavy metals. To investigate the retardation factor caused by adsorption, batch adsorption tests were conducted at various concentrations. Test results were correlated with both Langmuir and Freundlich isotherms. The adsorption of the lead ion was applicable to the Langmuir isotherm and the adsorption of the cadmium ion was applicable to the Freundlich isotherm. In addition, based on experimental results, the migration characteristics of heavy metals through the bed of fly ash and bentonite mixture were investigated using the PHREEQC, a reactive transport model, under the real conditions of the landfill liner.

Abrasive powders were recovered from electrical industry sludge by simple physical separation for their recycling. The raw electrical industry sludge was filter pressed, dried, dispersed and then classified by air classifier at various conditions. The three kinds of particles with different particle size distribution were classified by controlling rotor speed and air volumes of the classifier. The recovered abrasive powders, which are classified at 5,000,9000 and 13,000 rpm of rotor speed, are almost same properties to raw pumice, garnet and rouge powders, respectively. The results of particle size analysis, X-ray diffraction and SEM observation show that the recovered powders can be reused as an abrasive powders.

Palladium is widely used for several applications and recovery of palladium from secondary sources becomes increasingly important since palladium is one of maldistributed platinum group metals. Electrochemical recovery of dense palladium metal sheet from Pd leaching solution is a simple and easily controlled method. The surface morphology of the recovered Pd metal was significantly affected by current density and temperature. Dense deposit morphology was in higher stress state regardless of preparation condition under $55^{\circ}C$. Rising temperature up to $70^{\circ}C$ had a stress releasing effect besides densification of Pd deposit.

Optimized recovery of heavy minerals from the near shore sands of Korean Yellow Sea was investigated using physical processing technologies such as gravity concentration and magnetic separation. The head samples were subjected to the three stages effective separation; Head sample was first treated by a spiral separator to recover rough heavy mineral concentrates, which are contained minerals like ilmenite, zircon and rare earth minerals. Much higher beneficiation processes were subsequently taken by wilfley table and magnetic separation according to their magnetic field responses. Heavy minerals were effectively recovered by wilfley table and subsequent recleaning of heavy minerals by magnetic separations was conducted. Qualitative and relative-quantitative analyses of their constituent elements were doing using XRD and XRF.

We obtained amorphous calcium carbonate through the carbonation reaction of $Ca(OH)_2$, and through this reaction, observed changes in particle shape and phase by electric conductivity, XRD and TEM analysis. According to the result of the analysis, in the first declining stage of electric conductivity, amorphous calcium carbonate that has formed is coated on the surface of $Ca(OH)_2$ and obstructs its dissolution, and in the first recovery stage of electric conductivity, amorphous calcium carbonate is dissolved and re-precipitated and forms chains of fine calcite particles linearly joined. In the second decline of conductivity, viscosity increases due to the growth of chains of calcite particles, and finally the calcite particles are dissolved and separated into colloidal crystalline calcite, thereby increasing electric conductivity again.

This paper presents Lattice-Boltzmann simulation techniques for single-phase and two-phase fluid flow in porous media. Numerical experiments were performed in a digital rock sample from X-ray microtomography. Computed results showed very good agreement with laboratory measurements of permeability and relative permeability. Two applications using these simulation techniques show the potential of the Lattice-Boltzmann flow simulation to solve many difficult problems coupled with fluid flow in porous media.

This paper examines structural breaks and asymmetries of prices of four domestic petroleum products, i.e., gasoline, kerosene, diesel, and bunker-C, following the changes in the pricing policies pertaining to petroleum products in Korea from the government-controlled pricing system to the market pricing system. We use the monthly wholesale market price data for the sample period between July 1988 and December 2001. Using the four methods: the Chow test, the CUSUM/CUSUMQ tests, the Bayesian approach and the Dufour test, the structural behaviors of the petroleum product prices are examined. We found that structural change occurred in all petroleum products, with the exception of Kerosene, at the point of pricing policy change from government-controlled to the spot-price related pricing system. We, also conducted asymmetric analysis using the Borenstein, Cameron, and Gilbert (1997)'s model and found evidences of price asymmetry for all four product types, but in different pattern for each product.

A laboratory experiment using artificial fracture rocks was used to understand the 3-dimensional dispersion of a tracer and the mixing process in a fractured network. In this experiment, 12cm polystyrene foam cubes with two electrodes for monitoring electric conductivity (EC) were used as artificial fractured rocks. Distilled water with 0.5mS/m was used as a tracer in water with 35mS/m and the difference of EC between the tracer and the water was monitored by a multipoint simultaneous measurement system of electrical resistance. The results showed that even if the fracture arrangement pattern was not straight in the direction of the flow, the tracer did not diffuse along individual fractures and an oval tracer plume, which was the distribution of tracer concentrations, tended to be form in the direction of the flow. The vertical cross section of the tracer distribution showed small diffusivity in the vertical direction. The calculated total tracer volume passing through each measurement point in the horizontal cross section showed while that the solute passed through measurement points near the direction of hydraulic gradient and in other directions, the passed tracer volumes were small. Using Peclet number as a criterion, it was found that the mass distribution at the fracture intersection was controlled in the stage of transition between the complete mixing model and the streamline routing model.

One field identified by an inverse method is one of multiple candidate solutions those are independently obtained through a specific estimation technique. While averaging of optimized fields can provide a better description of the spatial feature of an unknown field, it deteriorates the flow and transport characteristics of the optimized fields. As a result, the averaged field is not suited for modeling aquifer performances. Based on genetic algorithm, an optimal field synthesis technique is developed, which combines diversely optimized fields into a refined group of fields. Each field in the population is paired, and a sub-region of each field is exchanged by crossover operation to create a group of synthesized fields of enhanced modeling capability. The population of the fields is evolved till the synthesized fields become sufficiently similar. Applications of the optimal field synthesis to synthetic cases indicate that the objective functions of the fields assessing the modeling capabilities are further reduced after the optimal field synthesis. The identified fields from various inverse techniques may yield a range of modeling results under varied flow situations. The uncertainty is narrowed down through the optimal field synthesis and the associated modeling results converge on that of the reference field. The developed inverse modeling facilitates the construction of a reliable simulation model and hence trustworthy predictions of the future performances.

Groundwater preferentially flows through sediment layers with high permeability such as colluvium. Its flow paths are called groundwater vein streams. An underground temperature survey is a method to locate vein streams by underground temperature anomalies associated with flowing groundwater. A groundwater flow system near an irrigation reservoir located in the upper part of a landslide block was surveyed with this method. After a geomembrane lining was installed in the reservoir, the total cross-sectional area of the vein streams in the aquifer decreased to as little as 0.35 times that before installation of the liner. A change in groundwater quality also indicated that the mixing of groundwater with leaked water from the reservoir stopped after installation of the lining.

In this study, an experimental apparatus has been designed and set-up to analyse the dissociating phenomena of hydrate in porous rock using electric heating method supplied at downhole. The electric heat injecting experiments have been performed to investigate the heat transfer within the core, the dissociating phenomena of hydrate, and the productivities of dissociated gas and water. These experiments were under constant heat injecting method as well as preheating methods. From the experimental results, it is seen that the hydrates is dissociated along the phase equilibrium curve and dissociation of hydrate is accelerated with heat. The injected heat is consumed for the dissociation and also it is lost together with outflow of the dissociated gas and water. From the investigation of gas producing behavior for various heat injecting methods, as the injected heat is greater, dissociation is accelerated faster at outlet and hence the initial gas production becomes higher. Also, it is shown that the initial gas productivity under the constant heating method is better, however, the energy efficiency is low because of smaller amount of the produced gas comparing to the amount of heat injected. In the experiments of preheating method, it was seen that gas production only initial stage is different with the preheating time, but the producing behaviors of gas production are similar.

In the sub-surface environments, detection of the movement of contaminant substances and recharge of groundwater by rainfall are very important factors which contain porosity and effective porosity of porous media. In this paper, the applicability of permittivity methods and proposed dielectric mixing models (DDMs) are discussed. This study showed that the ratio of effective porosity to porosity of Toyoura and River sands were 0.856 and 0.843. From the relationships between the relative porosity and effective porosity, all measured values can be confirmed to outside the range to about 0.800 for Toyoura and River sands under all experiments by FDR and FDR-V systems. In the study, this permittivity equipment can be considered to be good enough to measure determining the physical parameters of saturated soils. Consequently, this permittivity method can be contributed to estimate a porosity and effective porosity of saturated porous media because it is easy and instantaneous than previous in-situ methods.

For Taejon LNG Pilot Cavern being constructed to verify the technical aspects for storing LNG in lined rock cavern, efficiency tests of groundwater drainage system composed of many pumps and boreholes were performed around the cavern before and after the construction of concrete lining. Through evaluation of water balance and monitoring of pressures and flowrates, even if the present drainage system is very good for reducing water entries into the cavern, non-negligible water is still flowing in the floor of the cavern concrete due to heavy rainfall. To improve the drainage efficiency, additional drainage holes and some grouting were planned.

For the successful characterization of the fractured reservoir, it is inevitable to describe fracture properties more precisely. Although the deterministic features, of mega-trend faults can be relatively definite, the background features of minor fractures are not easily analysed in spite of the various data of these features. In this study, two different methods, statistical and fractal methods, are used to construct the minor fracture system over the entire field. After completing these tasks, with the identified deterministic features and background features, total fracture system for fractured reservoir is constructed and updated with the aid of dynamic data obtained from well test.

This paper presents the development of reservoir characterization model equipped with parallelized genetic algorithm, and its application for a heterogeneous reservoir system with integration of the well data and multi-phase production data. A parallel processing method performed by PC-cluster was applied to the developed model in order to reduce time for an inverse calculation. By utilizing the developed model, we performed the inverse calculation with the production data obtained from three layered reservoir system to estimate porosity and permeability distribution. As a result, the pressures observed at well almost identical to those calculated by the developed model. Also, it was confirmed that parallel processing could be applied for reservoir characterization study efficiently.

In general, well log and core data have been utilized for reservoir characterization. These well data can provide valuable information on reservoir properties with high vertical resolution at well locations. While the seismic surveys cover large areas of field but give only indirect features about reservoir properties. Therefore it is possible to estimate the reservoir properties guided by seismic data on entire area if a relationship of seismic data and well data can be defined. Seismic attributes calculated from seismic surveys contain the particular reservoir features, so that they should be extracted and used properly according to the purpose of study. The method to select the suitable seismic attributes among enormous ones is needed. The stepwise regression and fuzzy curve analysis based on fuzzy logics are used for selecting the best attributes. The relationship can be utilized to estimate reservoir properties derived from seismic attributes. This methodology is applied to a synthetic seismogram and a sonic log acquired from velocity model. Seismic attributes calculated from the seismic data are reflection strength, instantaneous phase, instantaneous frequency and pseudo sonic logging data as well as seismic trace. The fuzzy curve analysis is used for choosing the best seismic attributes compared to sonic log as well data, so that seismic trace, reflection strength, instantaneous frequency, and pseudo sonic logging data are selected. The relationship between the seismic attribute and well data is found out by the statistical regression method and estimates the reliable well data at a specific field location derived from only seismic attributes. For a future work in this study, the methodology should be checked an applicability of the real fields with more complex and various reservoir features.