• Geophysics and Geophysical Exploration
• /
• v.22 no.2
• /
• pp.80-87
• /
• 2019

In this paper, we propose a method to apply the polynomial fitting for regional-residual separation of microgravity data based on the characteristics of gravity anomaly without a prior information. Since the microgravity survey is usually carried out in small regions, it is common to approximate regional anomaly by the first-order polynomial plane. However, if the regional anomaly patterns are difficult to be approximated to a first-order plane, the complete gravity anomaly is divided into small zones enough to approximate first-order plane by means of Parasnis density estimation method. The regional-residual separation is then applied on the splitted zones individually. When the gravity anomalies can be splitted spatially, we showed that the residual anomalies can be more effectively extracted based on the regional geological structures by regional anomaly separation from each of the divided regions, rather than applying the entire data set at one time.

• Geophysics and Geophysical Exploration
• /
• v.22 no.2
• /
• pp.49-55
• /
• 2019

In this study, we presented a reference-shot subset method for stable convergence of full waveform inversion using a cyclic-shot subsampling technique. Full waveform inversion needs repetitive modeling of wave propagation and thus its calculation time increases as the number of sources increases. In order to reduce the computation time, we can use a cyclic-shot subsampling method; however, it makes the cost function oscillate in the early stage of the inversion and causes a problem in applying the convergence criteria. We introduced a method in which the cost function is calculated using a fixed reference-shot subset while updating the model parameters using the cyclic-shot subsampling method. Through the examples of full waveform inversion using the Marmousi velocity model, we confirmed that the convergence of cost function becomes stable even under the cyclic-shot subsampling method if using a reference-shot subset.

• Geophysics and Geophysical Exploration
• /
• v.22 no.1
• /
• pp.29-36
• /
• 2019

Regional airborne magnetic survey is very cost-effective mapping tool. Magnetic anomaly maps have abundant information, which are an important tool for understanding the geological evolution and mineral exploration. For this reason, the governments of many countries have made significant investment in the acquisition of airborne geophysical data over many decades. KIGAM (Korea Institute of Geoscience and Mineral Resources) began nationwide airborne magnetic mapping programme in 1982, and completed in 2017. The obtained magnetic data was reprocessed and magnetic database was built in 2018. In addition, the magnetic anomaly map of Korea with a scale of 1:1,000,000 was published. In this paper, we introduced a new magnetic anomaly map of Korea through describing the changing survey parameters during data acquisitions and history of data processing.

• Geophysics and Geophysical Exploration
• /
• v.22 no.1
• /
• pp.12-20
• /
• 2019

The first-arrival traveltime calculation method based on the damped wave equation overcomes the shortcomings of ray-tracing methods. Since this algorithm needs to solve the damped wave equation, numerical modeling is essential. However, it is not desirable to use the finite-difference method (FDM), which has good computational efficiency, for simulating the land seismic data because of irregular topography. Thus, the finite-element method (FEM) which requires higher computational cost than FDM has been used to correctly describe the irregular topography. In this study, we computed first-arrival traveltimes in an irregular topographic model using FDM incorporating embedded boundary method (EBM) to overcome this problem. To verify the accuracy and efficiency of the proposed algorithm, we compared our results with those of FEM. As a result, the proposed method using EBM not only provided the same accuracy as the FEM but also showed the improved computational efficiency.

• Geophysics and Geophysical Exploration
• /
• v.22 no.2
• /
• pp.56-61
• /
• 2019

In this paper, we studied efficient modeling strategies when we simulate the 3D time-domain acoustic wave propagation using a cell-based finite difference method which can handle the variations of both P-wave velocity and density. The standard finite difference method assigns physical properties such as velocities of elastic waves and density to grid points; on the other hand, the cell-based finite difference method assigns physical properties to cells between grid points. The cell-based finite difference method uses average physical properties of adjacent cells to calculate the finite difference equation centered at a grid point. This feature increases the computational cost of the cell-based finite difference method compared to the standard finite different method. In this study, we used additional memory to mitigate the computational overburden and thus reduced the calculation time by more than 30 %. Furthermore, we were able to enhance the performance of the modeling on several media with limited density variations by using the cell-based and standard finite difference methods together.

• Journal of The Geomorphological Association of Korea
• /
• v.25 no.4
• /
• pp.63-75
• /
• 2018

The lower marine terrace 1 and 2 surfaces distributed between Ulsan and Pohang coast in the southeastern coast of the Korean penninsula have been correlated with MIS 5e and 5a (or 5c) by amino acid dates, 14C dates, wide-spread tephra correlation and pollen analysis respectively. In this study, to test the reliability of the OSL method for the estimation of the numerical burial age of marine sediment deposits, we analyzed the samples from the marine terraces which have been known as typical marine terraces formed during MIS 5e and MIS 5a in the above-mentioned coast. The burial ages of the marine deposit of the lower marine terrace 1 and 2, with paleoshoreline altitudes of 18m to 19m and 10m to 11m respectively, both showed about the same age of 60 ka BP. The lower marine terraces 1 and 2, however, were divided into two terrace surfaces by a clear terrace cliff. Besides, the OSL dates of the lower and upper parts of the lower marine terrace 2 of the Bonggil coast showed the reversed burial ages. In the lower marine terrace 1 of the Sanhari coast, almost the same burial ages were derived from both the lower part (marine rounded gravel layer) and the upper part (terrestrial angular gravel layer) of the terrace deposit. Therefore, at the present time, judging from only the OSL dates measured in this study, it could be argued that the OSL method is not the best for the estimation of forming periods of the lower marine terraces 1 and 2 and their classification.

• Geomechanics and Engineering
• /
• v.30 no.1
• /
• pp.11-26
• /
• 2022

One of the most important issues in tunneling, is the squeezing phenomenon. Squeezing can occur during excavation or after the construction of tunnels, which in both cases could lead to significant damages. Therefore, it is important to predict the squeezing and consider it in the early design stage of tunnel construction. Different empirical, semi-empirical and theoretical-analytical methods have been presented to determine the squeezing. Therefore, it is necessary to examine the ability of each of these methods and identify the best method among them. In this study, squeezing in a part of the Alborz service tunnel in Iran was estimated through a number of empirical, semi- empirical and theoretical-analytical methods. Among these methods, the most robust model was used to obtain a database including 300 data for training and 33 data for testing in order to develop a machine learning (ML) method. To this end, three ML models of Gaussian process regression (GPR), artificial neural network (ANN) and support vector regression (SVR) were trained and tested to propose a robust model to predict the squeezing phenomenon. A comparative analysis between the conventional and the ML methods utilized in this study showed that, the GPR model is the most robust model in the prediction of squeezing phenomenon. The sensitivity analysis of the input parameters using the mutual information test (MIT) method showed that, the most sensitive parameter on the squeezing phenomenon is the tangential strain (ε_θ^α) parameter with a sensitivity score of 2.18. Finally, the GPR model was recommended to predict the squeezing phenomenon in tunneling projects. This work's significance is that it can provide a good estimation of the squeezing phenomenon in tunneling projects, based on which geotechnical engineers can take the necessary actions to deal with it in the pre-construction designs.

• The Journal of Engineering Geology
• /
• v.32 no.3
• /
• pp.363-376
• /
• 2022

An area under construction for a living facility collapsed around 12:48 KST on 13 January 2021 in Sa-dong, Ansan-si, Gyeonggi-do. There were no casualties due to the rapid evacuation measure, but part of the temporary retaining facility collapsed, and several cracks occurred in the adjacent road on the south side. This study used the potential of synthetic aperture radar (SAR) satellite for surface property changes that lies in backscattering characteristic to map the collapsed structure. The interferometric SAR technique can make a direct measurement of the decorrelation among different acquisition dates by integrating both amplitude and phase information. The damage proxy map (DPM) technique has been employed using four high-resolution Constellation of Small Satellites for Mediterranean basin Observation (COSMO-SkyMed) data spanning from 2020 to 2021 during ascending observation to analyze the collapse of the construction. DPM relies on the difference of pre- and co-event interferometric coherences to depict anomalous changes that indicate collapsed structure in the study area. The DPMs were displayed in a color scale that indicates an increasingly more significant ground surface change in the area covered by the pixels, depicting the collapsed structure. Therefore, the DPM technique with SAR data can be used for damage assessment with accurate and comprehensive detection after an event. In addition, we classify the amplitude information using support vector machine (SVM) and maximum likelihood classification algorithms. An investigation committee was formed to determine the cause of the collapse of the retaining wall and to suggest technical and institutional measures and alternatives to prevent similar incidents from reoccurring. The report from the committee revealed that the incident was caused by a combination of factors that were not carried out properly.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.24 no.3
• /
• pp.279-292
• /
• 2022

Unexpected anomalies in the geotechnical investigation at design stage may cause problems during tunnel excavation. Therefore, it is important to predict the ground condition ahead of a tunnel face during tunnel excavation in order to prevent tunnel collapse. Despite the fulfillment of an electrical resistivity survey at the tunnel face, the existing electrical resistivity survey program can produce distorted results by the limitation of tunnel modelling. In this background, this study develops a modelling program for an electrical resistivity survey considering the tunnel shape. Numerical simulation and inverse calculation were performed for the electrical resistivity survey in the tunnel using the developed program. As a result, it was proved that the developed program could predict accurately the anomalous object's location and condition ahead of the tunnel face.

• The Journal of Engineering Geology
• /
• v.32 no.2
• /
• pp.271-280
• /
• 2022

Landfast sea ice forms near coastlines in polar regions. Continuous monitoring of this sea ice is important, as it plays a key role in the marine ecosystem and affects the operation of nearby research stations. This study detected landfast sea ice around Jang Bogo research station in East Antarctica by stacking interferometric coherence images of Sentinel-1 synthetic aperture radar (SAR) data with 6-, 12- and 18-day temporal baselines. A total of 50 landfast sea ice maps were generated covering July 2017 to June 2018. The time series revealed regional differences in the timing of the maximum extent as well as growth rate of landfast sea ice. Overall, detecting landfast sea ice using interferometric SAR coherence seems promisingly feasible; however, limitations remain owing to low backscattering coefficients from new and smooth sea ice surfaces and subtle movements of sea ice in contact with the Campbell Glacier Tongue.