• Economic and Environmental Geology
• /
• v.31 no.3
• /
• pp.235-247
• /
• 1998

A new technique of simultaneons inversion for 3-D seismic velocity structure by using direct, reflected, and refracted waves is applied to the center of the Korean Peninsula including Pyongnam Basin, Kyonggi Massif, Okchon Fold Zone, Taebaeksan Fold Zone, Ryongnam Massif and Kyongsang Basin. Pg, Sg, PmP, SmS, Pn, and Sn arrival times of 32 events with 404 seismic rays are inverted for locations and crustal structure. 5 ($1^{\circ}$ along the latitude)${\times}6$ ($0.5^{\circ}$ along the longitude) ${\times}8$ block (4 km each layer) model was inverted. 3-D seismic crustal velocity tomography including eight sections from the surface to the Moho, eight profiles along latitude and longitude and the Moho depth distribution was determined. The results are as follows: (1) the average velocity and thickness of sediment are 5.15 km/sec and 3-4 km, and the velocity of basement is 6.12 km/sec. (2) the velocities fluctuate strongly in the upper crust, and the velocity distribution of the lower crust under Conrad appears basically horizontal. (3) the average depth of Moho is 29.8 km and velocity is 7.97 km/sec. (4) from the sedimentary depth and velocity, basement thickness and velocity, form of the upper crust, the Moho depth and form of the remarkable crustal velocity differences among Pyongnam Basin, Kyonggi Massif, Okchon Zone, Ryongnam Massif and Kyongsang Basin can be found. (5) The different crustal features of ocean and continent crust are obvious. (6) Some deep index of the Chugaryong Rift Zone can be located from the cross section profiles. (7) We note that there are big anisotropy bodies near north of Seoul and Hongsung in the upper crust, implying that they may be related to the Chugaryong Rift Zone and deep fault systems.

• Journal of the Korean earth science society
• /
• v.30 no.2
• /
• pp.153-164
• /
• 2009

We examined the regional 1-D deep resistivity structure of the Korean Peninsula using MT data acquired at seven sites located in the Kyongsang Basin and Kyonggi Massif. At the sites located in the Kyongsang Basin, surrounding sea distorts observed MT response and hence this distortion, so called "sea effect", is corrected using an iterative tensor stripping method. The 1-D layered inversion results for the seven MT sites reveal 4 layered structure, which is composed of 1) near surface layer, 2) upper crust, 3) lower crust and upper mantle, and 4) asthenosphere from the surface downward. Conrad interface, which is a boundary between upper and lower crust, is distinctly identified beneath all the MT sites. Conrad interface depth is estimated to about be 17km in the Kyongsang Basin and about 12km in the Kyonggi Massif, while the upper crust of the Kyongsang Basin is about 5 times more resistive than that of the Kyonggi Massif. Finally, asthenosphere is inferred to exist below a depth of approximately 100km with a resistivity of 200-300 ohm-m.

• 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.14 no.4
• /
• pp.316-323
• /
• 2011

To test the applicability of resistivity survey methods for the archaeological prospection of a large-scale tumulus, a three-dimensional resistivity survey was conducted at the $3^{rd}$ tumulus at Bokam-ri, in Naju city, South Korea. Since accurate topographic relief of the tumulus and electrode locations are required to obtain a high resolution image of the subsurface, electrodes were installed after making grids by threads, which is commonly used in the archaeological investigation. In the data acquisition, data were measured using a 2 m electrode spacing with the line spacing of 1 m and each survey line was shifted 1 m to form an effective grid of 1 m ${\times}$ 1 m. Though the 3-D inversion of data, we could obtain the 3-D image of the tumulus, where we could identify the brilliant signature of buried tombs made of stones. The results were compared with the previous excavation results and we could convince that a 3-D resistivity imaging method is very useful to investigate a large-scale tumulus.

• Journal of the Korean Geotechnical Society
• /
• v.21 no.4
• /
• pp.71-81
• /
• 2005

Surface wave techniques were initially based on 2-D plane waves and were later improved to the techniques based the 3-D based cylindrical waves. However, body-wave interference, near-field effect and limited technology in surface wave measurements restricted the use of 3-D cylindrical waves to the 1-D evaluation of subgrade stiffness. In this study, by the numerical simulation of SASW measurements, the dispersion properties of surface waves including vertical, horizontal Rayleigh waves and Love waves were thoroughly investigated in the 3-D domain, and a new filter criteria to minimize the near-field effect was established, which led to CAP (common-array-profiling)-SASW technique. The CAP-SASW technique enabled the evaluation of subgrade stiffness fur a specific subgrade segment, not for a whole section of measurement array. Therefore, a contour plot of subgrade stiffness with a ground-truth quality can be obtained by the CAP-SASW technique. The procedure proposed in this study was verified by comparing the shear-wave velocity profiles with the shear-wave velocity profiles of downhole testing at two geotechnical sites.

• Geophysics and Geophysical Exploration
• /
• v.4 no.3
• /
• pp.70-79
• /
• 2001

The integral equation method is a powerful tool for numerical electromagnetic modeling. But the difficulty of this technique is the size of the linear equations, which demands excessive memory and calculation time to invert. This limitation of the integral equation method becomes critical in inverse problem. The conventional Born approximation, where the electric field in the anomalous body is approximated by the background field, is very rapid and easy to compute. However, the technique is inaccurate when the conductivity contrast between the body and the background medium is large. Quasi-linear, quasi-analytical and extended Born approximations are novel approaches to 3-D EM modeling based on the linearization of the integral equations for scattered EM field. These approximation methods are much less time consuming than full integral equation method and more accurate than conventional Born approximation. They we, however, still approximate methods for 3-D EM modeling. Iterative series methods such as modified Born, quasi-linear and quasi-analytical can be used to increase the accuracy of various approximation methods. Comparisons of numerical performance against a full integral equation and various approximation codes show that the iterative series methods are very accurate and almost always converge. Furthermore, they are very fast and easy to implement on a computer. In this study, extended Born series method is developed and it shows more accurate result than that of other series methods. Therefore, Iterative series methods, including extended Born series, open principally new possibilities for fast and accurate 3-D EM modeling and inversion.

• Journal of the Korean earth science society
• /
• v.31 no.3
• /
• pp.214-224
• /
• 2010

Gravity and Magnetic survey data were analyzed to investigate the geophysical characteristics and regional geological structures of the southwestern Yellow Sea. The set of data about the southwestern part of the Yellow Sea in Korea was one collected by the Korea Ocean Research and Development Institute (KORDI) in 2003, 2004, and 2005. The Yellow Sea has a few basins and the study area also includes parts of the Heuksan Basin and the East China Sea Basin. The bathymetry of the study area ranges from about ?40 m southwestward near China to about 150 m northeastward near Korea. The bathymetry has the gentle rise and fall and the smooth slope. The gravity anomalies, from sea surface gravity and satellite gravity data, reflect the basement rocks rather than the smooth bathymetry. The gravity anomalies are higher on Northeastern part of the study area and lower over the South of the Heuksan Basin. The analytic signal from the Bouguer anomaly shows higher anomalous zones near the boundaries of the basins. The magnetic anomalies and the analytic signal, from the magnetic data, suggest that the complex anomalies on the Northern part are attributed to the volcanic intrusions and that the smooth patterns in the Southern part are based on the lack of the intrusions. The power spectrum analysis of the Bouguer anomalies and the magnetic anomalies indicate that the depth to the Moho discontinuity varies from about 30.2 to 28.3 km and that the depths of the basement rocks and the Eocene discontinuity range from about 8.4 to 8 km and from about 1.5 to 1.7 km, respectively. The inversion of the Bouguer anomaly shows that the Moho depth to the Western part of the study area near China is slightly deeper than the Eastern part near Korea. The result of 2-D gravity modeling has a good coherence with the results of the analytic signal, the power spectrum analysis, and the inversion.

• Journal of The Korean Astronomical Society
• /
• v.36 no.spc1
• /
• pp.37-44
• /
• 2003

In this paper, we review recent studies on the magnetic helicity changes of solar active regions by photospheric horizontal motions. Recently, Chae(200l) developed a methodology to determine the magnetic helicity change rate via photospheric horizontal motions. We have applied this methodology to four cases: (1) NOAA AR 8100 which has a series of homologous X-ray flares, (2) three active regions which have four eruptive major X-ray flares, (3) NOAA AR 9236 which has three eruptive X-class flares, and (4) NOAA AR 8668 in which a large filament was under formation. As a result, we have found several interesting results. First, the rate of magnetic helicity injection strongly depends on an active region and its evolution. Its mean rate ranges from 4 to $17 {\times} 10^{40}\;Mx^2\;h^{-1}$. Especially when the homologous flares occurred and when the filament was formed, significant rates of magnetic helicity were continuously deposited in the corona via photospheric shear flows. Second, there is a strong positive correlation between the magnetic helicity accumulated during the flaring time interval of the homologous flares in AR 8100 and the GOES X-ray flux integrated over the flaring time. This indicates that the occurrence of a series of homologous flares is physically related to the accumulation of magnetic helicity in the corona by photospheric shearing motions. Third, impulsive helicity variations took place near the flaring times of some strong flares. These impulsive variations whose time scales are less than one hour are attributed to localized velocity kernels around the polarity inversion line. Fourth, considering the filament eruption associated with an X1.8 flare started about 10 minutes before the impulsive variation of the helicity change rate, we suggest that the impulsive helicity variation is not a cause of the eruptive solar flare but its result. Finally, we discuss the physical implications on these results and our future plans.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2008.10a
• /
• pp.119-124
• /
• 2008

Due to the recent development of instruments which enabled the measurements of subtle IP effect in the ground and analysis algorithms, complex resistivity (CR) method was expanding its application to various field. In this study, we applied the CR method to the test site where the ground reinforcement had been done by injecting the cement mortar for investigating the effect of ground reinforcement. For this site, resistivity monitoring and tomography survey was carried out while the ground reinforcement had been made by the grouting. From the result, the anomalous region that was shown on the result of resistivity 4D monitoring was coincident with those of phase section in the CR method, because the cement grouting material had the strong IP effects. It might be expected that the CR method should be very powerful surveying tool for the similar purpose.

• Geophysics and Geophysical Exploration
• /
• v.4 no.3
• /
• pp.59-69
• /
• 2001

Five kinds of modified electrode arrays were proposed to overcome the weak points of the commonly used arrays using dipole and/or pole in two-dimensional resistivity surveys. The modified pole-pole array was suggested to overcome the inefficiency caused by distant earthing in pole-pole array. Four kinds of modified arrays using dipole were designed to enhance the signal-to-noise ratio of the conventional dipole-dipole and pole-dipole arrays through boosting up the measured potential difference. In the numerical experiments using the two-dimensional modeling and inversion, the effects of the ambient electrical noise and the resolving power were examined and the results showed the validity of the modified arrays proposed in this study.