• 한국지구물리탐사학회:학술대회논문집
• /
• 2006.06a
• /
• pp.161-166
• /
• 2006

Vertical gravity gradient measurement offers greater structural resolution and detectability than gravity alone. Practical difficulties of field measurement of vertical gravity gradient have raised questions of its accuracy and utility. But, modern automated gravimeter of $1\;{\mu}Gal$ sensitivity makes it easier to measure vertical gradient with required accuracy. It is particularly effective to engineering and environmental problems which target shallow subsurface structure. This paper attempts to apply the vertical gravity gradient technique to high resolution density mapping. The method was generally reviewed and numerical inverse modeling was executed for comparing with conventional gravity. And actual vertical gravity gradient data surveyed overt karstic cavity area at Muan was analysed and interpreted.

• Journal of the Korean Geophysical Society
• /
• v.8 no.4
• /
• pp.211-214
• /
• 2005

Low resistivity zone is observed at the lower part of a CFRD (Concrete Face Rockfill Dam). Generally, CFRD tends not to have any saturated zone within the body, but the result of resistivity survey shows that it is possible for the dam to be saturated under 20m depth with water. The level of reservoir was under 10 m from the crest. We suspect that this result may come from the wrong 2D inversion process ignoring the 3D geometry of dams. For the analysis of possibility of distortion by different geometry, we perform the 3D forward modeling for the dam and apply the 2D inversion process. And then we check the point of traditional interpretation of resistivity data. By the analysis, it is found that the result of 2D inversion process of 3D geometry of dams, seems to have deep relation with the reservoir level, and the complex 3D structure hide some internal electrical anomaly of dams from resistivity information.

• Geophysics and Geophysical Exploration
• /
• v.10 no.2
• /
• pp.111-123
• /
• 2007

Unlike the conventional reverse time migration method which is implemented by simply extrapolating wavefield in reverse time, this paper presents a derivation of another reverse time migration operator as the exact adjoint of the presumed forward wavefield extrapolation operator. The adjoint operator is obtained by formulating the forward time extrapolation operator in an explicit matrix equation form and then taking the adjoint to this matrix equation followed by determining the corresponding operator. The reverse time migration operator as the exact adjoint to the implied forward operator can be used not only as a migration algorithm but also as an adjoint operator which is required in the imaging through an inversion such as least-squares migration.

• 한국지구물리탐사학회:학술대회논문집
• /
• 1999.08a
• /
• pp.21-36
• /
• 1999

Recently resistivity survey is widely used for site investigations in the field of civil engineering. Since such application area requires accurate interpretation tools especially in the area of complicated geology and rough terrain topography, we developed a three-dimensional (3-D) resistivity inversion code, which can reconstruct real earth structures. Furthermore, the inversion code gives resolution-enhanced images by applying the ACB(Active Constraint Balancing) method. With the help of this inversion code, 3-D resistivity survey is now used as new techniques for site investigations in civil engineering problem. By imaging the 3-D resistivity distribution, we could get useful informations such as depth distribution of basement rock, distribution of weak zone, fractures and cavities which is crucial to civil engineers.

• Geophysics and Geophysical Exploration
• /
• v.6 no.1
• /
• pp.13-22
• /
• 2003

TSP (Tunnel Seismic Profiling) survey is a technique for imaging and characterizing geological structures ahead of a tunnel face. The seismic modeling algorithm and the synthetic data could be helpful for TSP surveys. However, there is few algorithm to describe the propagation of the elastic waves around the tunnel. In this study, existing 2-dimensional seismic modeling algorithm using finite element method was modified to make a suitable algorithm for TSP modeling. Using this algorithm, TSP modeling was practiced in some models. And the synthetic data was analyzed to examine the propagation characteristics of the elastic waves. First of all, the modeling for the homogeneous tunnel model was practiced to examine the propagation characteristics of the direct waves in the vicinity of the tunnel. And the algorithm was applied to some models having reflector which is perpendicular or parallel to the excavation direction. From these, the propagation characteristics of the reflected waves were examined. Furthermore, two source-receiver arrays were used in respective models to investigate the properties of the two arrays. These modeling algorithm and synthetic data could be helpful in interpreting TSP survey data, developing inversion algorithm and designing new source-receiver arrays.

• Geophysics and Geophysical Exploration
• /
• v.23 no.1
• /
• pp.22-37
• /
• 2020

We analyze the comprehensive three-dimensional (3D) magnetic structure characteristics from the seafloor to the deep layer of the Tofua Arc (TA) 12 seamount in the Tonga Arc, Southwestern Pacific, using bathymetric and geomagnetic data, and magnetization vector inversion (MVI) results. The seafloor features surrounding TA 12 highlight a NW-SE-oriented elliptical caldera at the summit of the seamount, two small cones in the depressed area. A large-scale sea valley is present on the western flank of the seamount, extending from these caldera cones to the southwestern base of the seamount. TA 12 seamount exhibits a low magnetic anomaly in the caldera depression, whereas a high magnetic anomaly is observed surrounding the low magnetic anomaly across the caldera summit. It is therefore presumed that there may be a strong magnetic material distribution or magma intrusion in the caldera. The 3D MVI results show that the high anomaly zones are mainly present in the surrounding slopes of the seamount from the seafloor to the -3,000 m (below the seafloor) level, with the magnetic susceptibility intensity increasing as the seafloor level increases at the caldera depression. However, small high anomaly zones are present across the study area near the seafloor level. Therefore, we expect that the magma ascent in TA 12 seamount migrated from the flanks to the depression area. Furthermore, we assume that the complex magnetic distribution near the seafloor is due to the remnant magnetization.

• Tunnel and Underground Space
• /
• v.13 no.6
• /
• pp.465-475
• /
• 2003

Groundwater inflow into the caverns constructed in fractured rock mass was simulated by numerical modeling, NAPSAC (DFN, discrete fracture network model) and NAMMU (CPM, continuous porous media model), a finite-element software package for groundwater flow in 3D fractured media developed by AEA Technology, UK. The input parameters for modeling were determined on surface fracture survey, core logging and single hole hydraulic test data. In order to predict the groundwater inflow more accurately, the anisotropic hydraulic conductivity was considered. The anisotropic hydraulic conductivities were calculated from the fracture network properties. With a minor adjustment during model calibration, the numerical modeling is able to reproduce reasonably groundwater inflows into cavern and the travel length and times to the ground surface along the flow paths in the normal, dry and rainy seasons.

• The Journal of Engineering Geology
• /
• v.31 no.4
• /
• pp.533-540
• /
• 2021

Gravity field analysis and density modeling were performed to evaluate the internal state of the rock mass, which is the cause of cut slope collapse. The shape of the weathered zone and the depth of basement could be confirmed from the complete Bouguer anomaly and density model. The basement depth at the center of the cut slope calculated using the Euler deconvolution inverse method is 30 m, which is about 10 m deeper than the surrounding area. In addition, the depth of basement and the thickness of the weathered zone are similar to the boundary between low resistivity and high resistivity in dipole-dipole survey. From the study results, gravity field analysis and density modeling recognizes the internal state of the rock slope and can be used for slope safety analysis, and is particularly suitable as a method to determine the shape of weathered zones in interpreting the safety of cut slopes

• Geophysics and Geophysical Exploration
• /
• v.14 no.1
• /
• pp.7-17
• /
• 2011

We have developed an inversion code for three-dimensional (3D) resistivity tomography including the anisotropy effect. The algorithm is based on the finite element approximations for the forward modelling and Active Constraint Balancing method is adopted to enhance the resolving power of the smoothness constraint least-squares inversion. Using numerical experiments, we have shown that anisotropic inversion is viable to get an accurate image of the subsurface when the subsurface shows strong electrical anisotropy. Moreover, anisotropy can be used as additional information in the interpretation of subsurface. This algorithm was also applied to the field dataset acquired in the abandoned old mine area, where a high-rise apartment block has been built up over a mining tunnel. The main purpose of the investigation was to evaluate the safety analysis of the building due to old mining activities. Strong electrical anisotropy has been observed and it was proven to be caused by geological setting of the site. To handle the anisotropy problem, field data were inverted by a 3D anisotropic tomography algorithm and we could obtain 3D subsurface images, which matches well with geology mapping observations. The inversion results have been used to provide the subsurface model for the safety analysis in rock engineering and we could assure the residents that the apartment has no problem in its safety after the completion of investigation works.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2005.09a
• /
• pp.93-116
• /
• 2005

SPT-Uphole tomography method was introducedand verified in this paper. In SPT-Uphole method, SPT (Standard Penetration Test) which is common in site investigation, was used as a source and several surface geophones in line were used as receivers. Shear wave velocity (Vs) distribution map which has triangular shape around the boring point can be obtained by tomography inversion. The factors for obtaining reliable result of SPT-Uphole tomography are exact travel time information and accurate inversion method. To establish of the SPT-Uphole tomography procedure, the most reliable method for obtaining exact travel time information and verification of tomography inversion method were studied by using theoretical travel time information and finite element method (FEM) analysis. finally, SPT-Uphole tomography method was performed at the weathered soil site in Kimje. By comparing with several boring data including SPT-N value, feasibility of this method was verified in the field.