• 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.

• 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.

• Journal of the Korean earth science society
• /
• v.32 no.6
• /
• pp.548-559
• /
• 2011

3-D Multi-geophysical surveys were carried out around the Hwasan caldera at the Euisung Sub-basin. To overcome the limitations of resolutions in previous studies, dense gravity data and magnetotelluric (MT) data were obtained and analyzed. In this study, the independent inversion models from gravity and MT data were integrated using correlation and classification approaches for 3-D imaging of the geologic structures. A Structure Index (SI) method was proposed and applied to the integration and classification analyses. This method consists of Type Angle (TA) and Type Intensity (TI) values, which are estimated by the spatial correlation and abnormality of the physical properties. The SI method allowed the classification analysis to be effectively performed. Major findings are as follows: 1) pyroclastic rocks around the central area of the Hwasan caldera with lower density and resistivity than those of neighboring regions extended to a depth of around 1 km, 2) intrusive igneous rocks with high resistivity and density were imaged around the ring fault boundary, and 3) a basement structure with low resistivity and high density, at a depth of 3-5 km, was inferred by the SI analysis.

• Geophysics and Geophysical Exploration
• /
• v.12 no.4
• /
• pp.361-366
• /
• 2009

Layered-earth Green's functions in electormagnetic (EM) surveys play a key role in modeling the response of exploration targets. They are computed through the Hankel transforms of analytic kernels. Computational precision depends upon the choice of algebraically equivalent forms by which these kemels are expressed. Since three-dimensional (3D) modeling can require a huge number of Green's function evaluations, total computational time can be influenced by computational time for the Hankel transform evaluations. Linear digital filters have proven to be a fast and accurate method of computing these Hankel transforms. In EM modeling for 3D inversion, electric fields are generally evaluated by the secondary field formulation to avoid the singularity problem. In this study, three components of electric fields for five different sources on the surface of homogeneous half-space were derived as primary field solutions. Moreover, reflection coefficients in TE and TM modes were produced to calculate EM responses accurately for a two-layered model having a sea layer. Accurate primary fields should substantially improve accuracy and decrease computation times for Green's function-based problems like MT problems and marine EM surveys.