• Geophysics and Geophysical Exploration
• /
• v.1 no.1
• /
• pp.49-56
• /
• 1998

Seismic tomography has been widely used as high resolution subsurface imaging techniques in engineering applications. Although most of the techniques have been using travel time inversion, waveform method is being driven forward owing to the progress of computational environments. Although full-waveform inversion method has been known as the best method in terms of model resolving power without high-frequency restriction and weak scattering approximation, it has practical disadvantage that it is apt to get stuck in local minimum if the initial guess is far from the actual model and it consumes so much time to calculate. In this study, 2-D full-waveform inversion algorithm in acoustic medium is developed, which uses result of traveltime tomography as initial model. From the application on synthetic data, it is proved that this approach can efficiently reduce the problem of conventional approaches: our algorithm shows much faster convergence rate and improvement of model resolution. Result of application on physical modeling data also shows much improvement. It is expected that this algorithm can be applicable to real data.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.3C
• /
• pp.143-155
• /
• 2006

SPT-Uphole tomography method was introduced for the evaluation of near subsurface shear wave velocity (Vs) distribution map. In SPT-Uphole method, SPT (Standard Penetration Test) which is common in geotechnical site investigation was used as a source and several surface geophones in line were used as receivers. Vs distribution map which is the triangular shape around the boring point can be developed by tomography inversion. To obtain the exact travel time information of shear wave component, a procedure using the magnitude summation of vertical and horizontal components was used based on the evaluation of particle motion at the surface. It was verified that proposed method could give reliable Vs distribution map through the numerical study using the FEM (Finite Element Method) model. Finally, SPT-Uphole tomography method was performed at the weathered soil site where several boring data with SPT-N values are available, and the feasibility of proposed method was verified in the field.

• The Journal of the Acoustical Society of Korea
• /
• v.13 no.6
• /
• pp.75-82
• /
• 1994

The ray paths and travel times of sound wave in the ocean depend on the physical properties of the propagating media. Ocean Acoustic Tomography(OAT), which is inversely estimate the travel time variations between fixed sources and receivers the physical properties of the corresponding media can he understood. To apply ocean survey technology by using the OAT, the tomographic procedure requires forward problem that variation of the travel times be identified with the variability of the medium. Also, received signals must be satisfied the necessary conditions of ray path stability, identification and resolution in order for OAT to work. The canonical ocean has been determined based on the historical data and its travel time and ray path are used as reference values. The sound speed of canonical ocean in the East Sea is about 1523 m/s at the surface and 1458 m/s at the sound channel axis(400m). Sound speeds in the East Sea are perturbed by warm eddy whose horizontal extension is more than 100 km with deeper than 200 m in depth scale. In this study, an acoustic source and receiver are placed at the depth above the sound channel axis, 350 m, and are separated by 200 km range. Ray paths are identified by the ray theory methed in a range dependent medium whose sound speeds are functions of a range and depth. The eigenray information obtained from interpolation between the rays bracketing the receiver are used to simulate the received signal by convolution of source signal with the eigenray informations. The source signal is taken as a 400 Hz rectangular pulse signal, bandwidth is 16 Hz and pulse length is 64 ms. According to the analysis of the received signal and identified ray path by using numerical model of underwater sound propagation, simulated signals satisfy the necessary conditions of OAT, applied in the East Sea.

• Geophysics and Geophysical Exploration
• /
• v.10 no.4
• /
• pp.314-321
• /
• 2007

Recently, multi-dimensional joint inversion of geophysical data based on fundamentally different physical properties is being actively studied. Joint inversion can provide a way to obtaining much more accurate image of the subsurface structure. Through the joint inversion, furthermore, it is possible to directly estimate non-geophysical material properties from geophysical measurements. In this study, we developed a new algorithm for jointly inverting dc resistivity and seismic traveltime data based on the multiple constraints: (1) structural similarity based on cross-gradient, (2) correlation between two different material properties, and (3) a priori information on the material property distribution. Through the numerical experiments of surface dc resistivity and seismic refraction surveys, the performance of the proposed algorithm was demonstrated and the effects of different regularizations were analyzed. In particular, we showed that the hidden layer problem in the seismic refraction method due to an inter-bedded low velocity layer can be solved by the joint inversion when appropriate constraints are applied.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2008.03a
• /
• pp.142-152
• /
• 2008

SPT-Uphole method was introduced for the evaluation of near subsurface shear wave velocity (Vs) profile. In SPT-Uphole method, SPT (Standard Penetration Test) which is common in geotechnical site investigation was used as a source and several surface geophones in line were used as receivers. 1D shearwave velocity profile can be obtained in the manner of downhole method, Vs distribution map which is the triangular shape around the boring point can be developed by tomography inversion. To obtain the exact travel time information of shear wave component, a procedure using the magnitude summation of vertical and horizontal components was used based on the evaluation of particle motion at the surface. It was verified that proposed method could give reliable Vs distribution map through the numerical study using the FEM (Finite Element Method) model. Finally, SPT-Uphole method was performed and the feasibility of proposed method was verified in the field.

• Economic and Environmental Geology
• /
• v.57 no.1
• /
• pp.1-15
• /
• 2024

The Melanesian region in the western Pacific is dominated by complex plate tectonics, with the largest oceanic plateau, the OntongJava plateau, and a hotspot, the Caroline Islands. To better understand the complex geodynamics of the region, we estimate P- and S-velocity models and 𝛿 (V_P/V_S) model by using relative teleseismic travel times measured at seismometers on land and the seafloor. Our results show high-velocity anomalies in the subduction zones of the Melanesian region to a depth of about 400 km, which is thought to be subducting Solomon Sea, Bismarck, and Australian plates along plate boundaries. Along subduction zones, positive 𝛿 (V_P/V_S) anomalies are found, which may be caused by partial melting due to dehydration. A broad high-velocity anomaly is observed at 600 km depth below the Ontong-Java plateau, with a negative 𝛿 (V_P/V_S) anomaly. This is thought to be a viscous and dry remnant of the Pacific plate that subducted at 45-25 Ma, with a low volume of fluids due to dehydration for a long period in the mantle transition zone. Beneath the Caroline Islands, a strong low-velocity anomaly is obseved to a depth of 800 km and appears to be connected to the underside of the remnant Pacific plate in the mantle transition zone. This suggests that the mantle plume originating in the lower mantle has been redirected due to the interaction with the remnant Pacific plate and has reached its current location. The mantle plume also has a positive 𝛿 (V_P/V_S) anomaly, which is thought to be due to the influence of embedded fluids or partial melting. A high-velocity anomaly, interpreted as an effect of the thick lithosphere beneath the Ontong-Java plateau, is observed down to 300 km depth with a negative 𝛿 (V_P/V_S) anomaly, which likely indicate that little fluid remains in the melt residue accumulated in the lithosphere.

• Nuclear Engineering and Technology
• /
• v.55 no.2
• /
• pp.681-689
• /
• 2023

Purpose: Muons are characterized by a strong penetrating ability and can travel through thousands of meters of rock, making them ideal to image large volumes and substances typically impenetrable to, for example, electrons and photons. The feasibility of 3D image reconstruction and material identification based on a cosmic ray muons tomography (MT) system with triangular bar plastic scintillator detectors has been verified in this paper. Our prototype shows potential application value and the authors wish to apply this prototype system to 3D imaging. In addition, an MT experiment with the same detector system is also in progress. Methods: A simulation based on GEANT4 was developed to study cosmic ray muons' physical processes and motion trails. The yield and transportation of optical photons scintillated in each triangular bar of the detector system were reproduced. An image reconstruction algorithm and correction method based on muon scattering, which differs from the conventional PoCA algorithm, has been developed based on simulation data and verified by experimental data. Results: According to the simulation result, the detector system's position resolution is below 1 ~ mm in simulation and 2 mm in the experiment. A relatively legible 3D image of lead bricks in size of 20 cm × 5 cm × 10 cm used our inversion algorithm can be presented below 1× 10⁴ effective events, which takes 16 h of acquisition time experimentally. Conclusion: The proposed method is a potential candidate to monitor the cosmic ray MT accurately. Monte Carlo simulations have been performed to discuss the application of the detector and the simulation results have indicated that the detector can be used in cosmic ray MT. The cosmic ray MT experiment is currently underway. Furthermore, the proposal also has the potential to scan the earth, buildings, and other structures of interest including for instance computerized imaging in an archaeological framework.

• Geophysics and Geophysical Exploration
• /
• v.9 no.4
• /
• pp.279-290
• /
• 2006

The borehole radar methods used to tunnel detection are mainly classified into borehole radar reflection, directional antenna, crosshole scanning, and radar tomography methods. In this study, we have investigated the feasibility and limitation of each method to tunnel detection through case studies. In the borehole radar reflection data, there were much more clear diffraction signals of the upper wings than lower wings of the hyperbolas reflected from the tunnel, and their upper and lower wings were spreaded out to more than 10m higher and lower traces from the peaks of the hyperbolas. As the ratio of borehole diameter to antenna length increases, the ringing gets stronger on the data due to the increase in the impedance mismatching between antennas and water in the boreholes. It is also found that the reflection signals from the tunnel could be enhanced using the optimal offset distance between transmitter and receiver antennas. Nevertheless, the borehole radar reflection data could not provide directional information of the reflectors in the subsurface. Direction finding antenna system had a advantage to take a three dimensional location of a tunnel with only one borehole survey even though the cost is still very high and it required very high expertise. The data from crosshole scanning could be a good indicator for tunnel detection and it could give more reliable result when the borehole radar reflection survey is carried out together. The images of the subsurface also can be reconstructed using travel time tomography which could provide the physical property of the medium and would be effective for imaging the underground structure such as tunnels. Based on the results described above, we suggest a cost-effective field procedure for detection of a tunnel using borehole radar techniques; borehole radar reflection survey using dipole antenna can firstly be applied to pick up anomalous regions within the borehole, and crosshole scanning or reflection survey using directional antenna can then be applied only to the anomalous regions to detect the tunnel.

• Geophysics and Geophysical Exploration
• /
• v.16 no.3
• /
• pp.109-118
• /
• 2013

On and near the 23-m high earthen Cheongcheon dam in Boryeong City, Korea, short seismic refraction and surface-wave profiles were conducted using a 5-kg sledgehammer. From vertical and horizontal components of the seismic waves, near-surface P-wave velocities (${\nu}_p$) and S-wave velocities (${\nu}_s$) were derived by inverting first-arrival refraction times and dispersion curves of Rayleigh waves. Average ${\nu}_p$ and ${\nu}_s$ for the Jurassic sedimentary basement were determined to be 1650 and 950 m/s at a depth of 30 m directly beneath the dam and 1650 m/s and 940 m/s at a depth of 10 m at the toe of the dam, respectively. The dynamic Poisson's ratio for these strata were therefore in the range of 0.24 to 0.25, which is consistent with ratios for consolidated sedimentary strata. Near a 45-m borehole 152 m downstream from the dam crest, an SH tomogram indicates a refraction boundary with an average ${\nu}_s$ of 870 m/s at depths of 10 ~ 12 m. At this site, the overburden comprises the upper layer with relatively constant ${\nu}_p$ and ${\nu}_s$ around 500 and 200 m/s, respectively, and the lower layer in which both ${\nu}_p$ and ${\nu}_s$ increase with depth almost linearly. The dynamic Poisson's ratios for the overburden were in the range of 0.30 to 0.43.

• Economic and Environmental Geology
• /
• v.39 no.2 s.177
• /
• pp.111-128
• /
• 2006

Three dimensional crustal structure and source features of earthquake hypocenters on the Korean peninsula were investigated using P and S-wave travel time tomography. The main goal of this research was to find Vp/Vs anomalies at earthquake hypocenters as well as those of crustal structure of basins and deep tectonic settings. This allowed fer the extrapolation of more detailed seismotectonic force from the Korean peninsula. The earthquake hypocenters were found to have high Vp/Vs ratio discrepancies (VRD) at the vertical sections. High V/p/Vs ratios were also found in the sedimentary basins and beneath the Chugaryong Rift Zone (CRZ), which was due to mantle plume that subsequently solidified with many fractures and faults which were saturated with connate water. The hypocenters of most earthquakes were found in the upper crust for Youngwol (YE), Kyongju (KE), Hongsung (HE), Kaesong (KSE), Daekwan (DKE), and Daehung (DHE) earthquakes, but near the subcrust or the Moho Discontinuity for Mt. Songni (SE), Sariwon (SRE) and Mt. Jiri (JE) earthquakes. Especially, we found hot springs of the Daekwan, Daehung and Unsan regions coincide with high VRD. Also, this cannot rule out the possibility that there are some partial meltings in the subcrust of this region. High VRD might indicate that many faults and fractures with connate water were dehydrated when earthquakes took place, reducing shear modulus in the hypocenter areas. This is can be explained by due to the fact that a point source which is represented by the moment tensor that may involve changes in volume, shear fracture, and rigidity. High Vp/Vs ratio discrepancies (VRD) were also found beneath Mt. Backdu beneath 40 km, indicating that magma chamber existed beneath Mt. Backdu is reducing shear modulus of S-wave velocity.