• Geophysics and Geophysical Exploration
• /
• v.17 no.2
• /
• pp.49-57
• /
• 2014

3D magnetic inversion, based on the assemblage of 2D forward modeling and inversion as a practical technique to reflect the a priori information, was conducted to investigate the spatial distribution features of black-shale related and pyrometamorphic uranium deposit and several lithological units of Ogcheon Super Group in an area of Geumsan. By using the 3D visualization technique with suitable susceptibility interval and horizontal slice map, the spatial distribution of magnetic susceptibility corresponded to the black shale related uranium bearing lithological units, Black Slate member was well coincided with a information of uranium deposit. Also, even though it is indirect indicator for the detetction of uranium deposits interbedded in Gray Hornfels member, spatial susceptibility distribution which shows the south-east magnetic linearment corresponding to the Majeon-ri formation and Dark Gray Slate were matched well. From this investigation, we inferred that maximum depth extension which Black Slate member can be separately recognized with respect to adjacent Dark Gray Slate with strong magnetic susceptibility anomaly is about 150m with reference elevation level of 306m. In addition, Majeon-ri formation located south of Black Slate member revels relatively high magnetic susceptibility range but shows high spatial susceptibility fluctation. And, as an intrusive rocks, Jurassic Biotite Granite shows relatively low magnetic susceptibility characteristics. On the contrary, Cretaceous granite distributed in soutthern part of the study area shows the relatively high susceptibility distribution.

• Geophysics and Geophysical Exploration
• /
• v.15 no.3
• /
• pp.136-143
• /
• 2012

In this case study, we present the various and consistent processing techniques for the reasonable interpretation of aeromagnetic data. In the processing stage, we especially focused on the three major respects. First, in the low latitude area, severe artifacts are occurred as a result of reduction to the pole technique. To overcome this problem, variable alternative methods were investigated. From the comparison of each technique, we concluded that energy balancing method gives more fruitful result. Second, because of limited a priori information, it is nearly impossible to employ detailed geological survey due to wide and thick spreading of soils in the survey area. So we especially investigated the new techniques such as extracting slope, curvature and aspect information mainly used in GIS field as well as conventional methods. Finally, by using the Euler deconvolution, we extracted the depth information on the magnetic anomalous body. From the synthetic analysis between depth information and previous discussed results, the detailed future survey area was proposed. We think that a series of processing techniques discussed in this study may perform an important role in the domestic and abroad resource development project as a useful guideline.

• Geophysics and Geophysical Exploration
• /
• v.15 no.3
• /
• pp.129-135
• /
• 2012

The simultaneous-source full waveform inversion improves the applicability of full waveform inversion by reducing the computational cost. Since this technique adopts simultaneous multi-source for forward modeling, unwanted events remain in the residual seismograms when the receiver geometry of field acquisition is different from that of numerical modeling. As a result, these events impede the convergence of the full waveform inversion. In particular, the streamer-type data with limited offsets is the most difficult data to apply the simultaneous-source technique. To overcome this problem, the global-correlation-based objective function was suggested and it was successfully applied to the simultaneous-source full waveform inversion in time domain. However, this method distorts residual wavefields due to the modified objective function and has a negative influence on the inversion result. In addition, this method has not been applied to the frequency-domain simultaneous-source full waveform inversion. In this paper, we apply a timedamping function to the observed and modeled data, which are used to compute global correlation, to minimize the distortion of residual wavefields. Since the damped wavefields optimize the performance of the global correlation, it mitigates the distortion of the residual wavefields and improves the inversion result. Our algorithm incorporates the globalcorrelation-based full waveform inversion into the frequency domain by back-propagating the time-domain residual wavefields in the frequency domain. Through the numerical examples using the streamer-type data, we show that our inversion algorithm better describes the velocity structure than the conventional global correlation approach does.

• Geophysics and Geophysical Exploration
• /
• v.22 no.1
• /
• pp.21-28
• /
• 2019

We perform a 3D tomographic inversion using surface wave dispersion curves to obtain S-velocity model and radial anisotropy beneath Saudi Arabia. The Arabian Peninsula is geologically and topographically divided into a shield and a platform. We used event data with magnitudes larger than 5.5 and epicentral distances shorter than $40^{\circ}$ during 2008 ~ 2014 from the Saudi Geological Survey. We obtained dispersion curves by using the multiple filtering technique after preprocessing the event data. We constructed SH- and SV-velocity models and consequently radial anisotropy model at 10 ~ 60 km depths by inverting Love and Rayleigh group velocity dispersion curves with period ranges of 5 ~ 140 s, respectively. We observe high-velocity anomalies beneath the Arabian shield at 10 ~ 30 km depths and low-velocity anomalies beneath the Arabian platform at 10 km depth in the SV-velocity model. This discrepancy may be caused by the difference between the Arabian shield and the Arabian platform, that is, the Arabian shield was formed in Proterozoic thereby old and cold, while the Arabian platform is covered by predominant Paleozoic, Mesozoic, and Cenozoic sedimentary layers. Also we obtained radial anisotropy by estimating the differences between SH- and SV-velocity models. Positive anisotropy is observed, which may be generated by lateral tension due to the slab pull of subducting slabs along the Zagros belt.

• Geophysics and Geophysical Exploration
• /
• v.24 no.4
• /
• pp.180-193
• /
• 2021

The complex resistivity method is an exploration technique that can obtain various characteristic information of underground media by measuring resistivity and phase in the frequency domain, and its utilization has recently increased. In this paper, a three-dimensional inversion algorithm for the CR data was developed to increase the utilization of this method. The Poisson equation, which can be applied when the electromagnetic coupling effect is ignored, was applied to the modeling, and the inversion algorithm was developed by modifying the existing algorithm by adopting comlex variables. In order to increase the stability of the inversion, a technique was introduced to automatically adjust the Lagrangian multiplier according to the ratio of the error vector and the model update vector. Furthermore, to compensate for the loss of data due to noisy phase data, a two-step inversion method that conducts inversion iterations using only resistivity data in the beginning and both of resistivity and phase data in the second half was developed. As a result of the experiment for the synthetic data, stable inversion results were obtained, and the validity to real data was also confirmed by applying the developed 3D inversion algorithm to the analysis of field data acquired near a hydrothermal mine.

• Geophysics and Geophysical Exploration
• /
• v.25 no.3
• /
• pp.109-119
• /
• 2022

The computational efficiency of reverse time migration (RTM) based on numerical modeling is not secured due to the high-frequency band of several hundred Hz or higher for data acquired through a three-dimensional (3D) ultra-high-resolution (UHR) seismic survey. Therefore, this study develops an RTM program to derive high-quality 3D geological structures using UHR seismic data. In the traditional 3D RTM program, an excitation amplitude technique that stores only the maximum amplitude of the source wavefield and a domain-limiting technique that minimizes the modeling area where the source and receivers are located were used to significantly reduce memory usage and calculation time. The program developed through this study successfully derived a 3D migration image with a horizontal grid size of 1 m for the 3D UHR seismic survey data obtained from the Korea Institute of Geoscience and Mineral Resources in 2019, and geological analysis was conducted.

• Geophysics and Geophysical Exploration
• /
• v.26 no.3
• /
• pp.114-125
• /
• 2023

The physics-informed neural network (PINN) has been proposed to overcome the limitations of various numerical methods used to solve partial differential equations (PDEs) and the drawbacks of purely data-driven machine learning. The PINN directly applies PDEs to the construction of the loss function, introducing physical constraints to machine learning training. This technique can also be applied to wave equation modeling. However, to solve the wave equation using the PINN, second-order differentiations with respect to input data must be performed during neural network training, and the resulting wavefields contain complex dynamical phenomena, requiring careful strategies. This tutorial elucidates the fundamental concepts of the PINN and discusses considerations for wave equation modeling using the PINN approach. These considerations include spatial coordinate normalization, the selection of activation functions, and strategies for incorporating physics loss. Our experimental results demonstrated that normalizing the spatial coordinates of the training data leads to a more accurate reflection of initial conditions in neural network training for wave equation modeling. Furthermore, the characteristics of various functions were compared to select an appropriate activation function for wavefield prediction using neural networks. These comparisons focused on their differentiation with respect to input data and their convergence properties. Finally, the results of two scenarios for incorporating physics loss into the loss function during neural network training were compared. Through numerical experiments, a curriculum-based learning strategy, applying physics loss after the initial training steps, was more effective than utilizing physics loss from the early training steps. In addition, the effectiveness of the PINN technique was confirmed by comparing these results with those of training without any use of physics loss.

• Journal of Soil and Groundwater Environment
• /
• v.7 no.3
• /
• pp.33-44
• /
• 2002

A field technique for assessing the transmissive fractures in an aquifer was applied to a fractured rock formation in Youngchun area Korea. Geological mapping and detailed acoustic borehole teleview(BHTV) logging were performed to obtain information about the fractures. The study area consists predominantly of two types of fractures. The fracture sets of low angle partings such as bedding and sheeting plains have strike N70-80$^{\circ}$W, 25$^{\circ}$-30$^{\circ}$SW and N3S$^{\circ}$W, 12$^{\circ}$NE, respectively. In areas of high fractures, on the other hand, the major fracture sets show strike N80$^{\circ}$W and dip 70$^{\circ}$-85$^{\circ}$SW, N10$^{\circ}$E.85$^{\circ}$SE in sedimentry rocks, N40-50$^{\circ}$E.85$^{\circ}$SE/85$^{\circ}$NE, N70$^{\circ}$E.80$^{\circ}$SE, and N7$^{\circ}$-75$^{\circ}$W.80$^{\circ}$SW in granites and volcanic rocks. Injection tests have been performed to identify discrete production zones and quantify the vertical distribution of hydraulic conductivity. The calculated hydraulic conductivities range from 3.363E-10 to 2.731E-6, showing that the difference between maximum and minimum value is four order of magnitude. Dominant section in hydraulic conductivity is extensively fractured. Geophysical logging was carried out to clarify characterization of the distribution of fracture zones. Transmissive fractures were evaluated through the comparison of the results obtained by each method. The temperature logs appeared to be a good indicator that can distinguish a high transmissive fractures from a common fractures in hydraulic conductivity. In numerous cases, evidence of fluid movement was amplified in the temperature gradient log. The fracture sets of N70-80$^{\circ}$W.60-85$^{\circ}$NE/SW N75-80$^{\circ}$W.25-30$^{\circ}$SW, N50-64$^{\circ}$W.60-85$^{\circ}$NE, N35-45$^{\circ}$E.65-75$^{\circ}$SE, and N65-72$^{\circ}$E.80$^{\circ}$SE/60$^{\circ}$NW were idenfied as a distinct transmissive fractures through the results of each tests.

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

This study is the case of groundwater development based on the geometrical fracture model of target area established only through geological fracture mapping technique. A fracture mapping of $9km^2$, eastern Munsan, has been conducted to determine geological and hydrological factors for new water well placement in the Gyeonggi gneiss complex. Geophysical exploration was not applicable because of small restricted area and dense underground utilities at the site. Form line mapping on the basis of foliation orientation and rock type revealed a synform of NS fold axis bearing to the south. An EW geological cross-section passed through the site area shows a F2 synform as a double-wall ice cream spoon shape. Three regional faults of $N20^{\circ}E,\;N30^{\circ}W$, and NS have been dragged into the site to help understand extensional fault paths. The $N20^{\circ}E$ fault with dextral sense is geometrically interpreted as a western fault of two flexural conjugate type-P shear faults in the F2 synformal fold. The NE cross-section reveals that a possible groundwater belt in the western limb of super-posed fold area is formed as a trigonal prism within 100 m depth of the intersectional space between the $N20^{\circ}E$ fault plane and the weakly sheared plane of transposed foliation. Another possible fault for water resource strikes $N40^{\circ}E$. Recommended sites for new water well placement are along the $N20^{\circ}E\;and\;N40^{\circ}E$ faults. As a result of fracture mapping, 145 ton/day of water can be produced at one well along the $N20^{\circ}E$ fault line. Exploration of groundwater in the area is succeeded only using with geological fracture mapping and interpretation of geological cross-section, without any geophysical survey. Intersection of fault generated with the F2 synformal fold and foliation supply space of groundwater reserver.

• Geophysics and Geophysical Exploration
• /
• v.16 no.3
• /
• pp.154-162
• /
• 2013

Vertical electrical sounding and 2D electrical resistivity survey were applied for evaluating the characteristics of alluvial layers at a groundwater artificial recharge site. The fine particles in alluvial layer, main target layer of groundwater artificial recharge, may cause clogging phenomena. In this case, electrical resistivity method is an effective technique to verify the spatial distribution of low-resistivity layers, such as saturated silts and clays. On the other hand, much attention should be paid to interpret the resistivity data in unconsolidated layers, because thick clayey overburden sometimes produces a masking effect on underlying interbedded resistive sands and gravels. Considering these points, we designed 35 points arranged in a grid form for vertical electrical sounding and 10 lines for 2D electrical resistivity survey, and concentrated our effort on enhancing the vertical and horizontal resolution of resistivity images. According to the results, 15 meters thick layers consisting of sands and gravels are located in 30 meters below ground. And the spatial distribution of silts and clays are mapped, which may cause clogging. Consequently, this approach can contribute to design and determine the location and depth of injection and observation wells for groundwater artificial recharge.