• Geophysics and Geophysical Exploration
• /
• v.20 no.3
• /
• pp.163-175
• /
• 2017

In land seismic exploration, irregular surface topography and weathering layer in near surface distorts the reflected signals of data. Therefore, typical land seismic data should be compensated for this distortion by static correction. To perform the static correction, near-surface velocity is required, which can be obtained by seismic refraction survey. However, land seismic data is often acquired in a limited form of geometry depending on the equipment availability, accessibility condition, and permission for the survey site. In this situation, refraction analysis should be performed using reflection data because it is impossible to acquire refraction-oriented data due to limited source and receiver geometry. In this study, we aimed to analyze the reliability of the results obtained by refraction traveltime tomography when using reflection data with a limited number of sources and receivers from irregular surface topography. By comparing the inversion result from irregular topography with that from flat surface, we found that the surface topography affects the reliability of the inversion results to some degree. We also found that the number of sources has little effect on the inversion results unless the number of sources are very small. On the other hand, we observed that velocity distortion occurred in the overlapped part of receiver arrays when using a limited number of receivers, and therefore suggested the size of the least overlapping ratio to avoid the velocity distortion. Finally, we performed numerical tests for the model which simulates the surface topography and acquisition geometry of the survey region and verified the reliability analysis of inversion results. We identified reliable areas and suspicious area of the inverted velocity model by applying the analysis results to field data.

• Geophysics and Geophysical Exploration
• /
• v.9 no.4
• /
• pp.261-270
• /
• 2006

Two-dimensional (2-D) inversion of magnetotelluric (MT) data for two survey lines having south-north direction from Jeju Island has been carried out. Broad band MT sounding curves with good quality could be gathered by performing audio-frequency magnetotelluric (AMT) survey during the MT survey and by operating the remote reference in Kyushu Island, Japan. Comparison of the 2-D inversion model using MT band only and that using both AMT and MT bands for the field data as well as for the data from numerical 2-D modeling said that high frequency information from AMT survey can be useful for interpreting not only the shallow part but also the deep structures, especially when the formation is resistive. The 2-D inversion models of field data show a thick layer having around 10 ohm-m in the depth of a few hundred meters throughout the survey area, which can be considered as the unconsolidated sedimentary layer. And they also show a conductive anomaly at the central part of each survey lines. It can be either the effect of the surrounding sea water, or the structures due to ancient volcanic events. But unfortunately by now, we do not have any further information about the anomaly.

• 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.26 no.4
• /
• pp.268-274
• /
• 2023

Finite source inversion is performed with a Green's function matrix and geodetic coseismic displacement. Conventionally, the Green's function matrix is constructed using the Okada model (Okada, 1985). However, for more realistic earthquake simulations, recent research has widely adopted the physics-based model, which can consider various material properties such as elasticity, viscoelasticity, and elastoplasticity. We used the physics-based software PyLith, which is suitable for earthquake modeling. However, the PyLith does not provide a mesh generator, which makes it difficult to perform finite source inversions that require numerous subfaults and observation points within the model. Therefore, in this study, we developed CPInterface (COMSOL-PyLith Interface) to improve the convenience of finite source inversion by combining the processes of creating a numerical model including sub-faults and observation points, simulating earthquake modeling, and constructing a Green's function matrix. CPInterface combines the grid generator of COMSOL with PyLith to generate the Green's function matrix automatically. CPInterface controls model and fault information with simple parameters. In addition, elastic subsurface anomalies and GPS observations can be placed flexibly in the model. CPInterface is expected to enhance the accessibility of physics-based finite source inversions by automatically generating the Green's function matrix.

• Journal of the Korean GEO-environmental Society
• /
• v.25 no.2
• /
• pp.5-14
• /
• 2024

Reinforced concrete bridge decks are the first to be damaged by vehicle loads and rain infiltration. Concrete deterioration primarily occurs owing to the corrosion of rebars and other metal components by chlorides used for snow and ice melting. The structural condition and concrete deterioration of the bridge decks within the pavement were evaluated using ground-penetrating radar (GPR) survey data. To evaluate concrete deterioration in bridges, it is necessary to develop GPR data analysis techniques to accurately identify deteriorated locations and rebar positions. GPR exploration involves the acquisition of reflection and diffraction wave signals due to differences in radar wave propagation velocity in geotechnical media. Therefore, a full-waveform inversion (FWI) method was developed to evaluate the deterioration of reinforced concrete bridge decks by estimating the radar wave propagation velocity in geotechnical media using GPR data. Numerical experiments using a GPR velocity model confirmed the deterioration phenomena of bridge decks, such as concrete delamination and rebar corrosion, verifying the applicability of the developed technology. Moreover, using the synthetic GPR data, FWI facilitates the determination of rebar positions and concrete deterioration locations using inverted velocity images.

• International Union of Geodesy and Geophysics Korean Journal of Geophysical Research
• /
• v.25 no.1
• /
• pp.57-81
• /
• 1997

Numerical prediction of nocturnal thermal high in summer of the 1995 near Taegu city located in a basin has been carried out by a non-hydrostatic numerical model over complex terrain through one-way double nesting technique in the Z following coordinate system. Under the prevailing westerly winds, vertical turbulent fluxes of momentum and heat over mountains for daytime hours are quite strong with a large magnitude of more than $120W/\textrm{m}^2$, but a small one of $5W/\textrm{m}^2$ at the surface of the basin. Convective boundary layer (CBL) is developed with a thickness of about 600m over the ground in the lee side of Mt. Hyungje, and extends to the edge of inland at the interface of land sea in the east. Sensible heat flux near the surface of the top of the mountain is $50W/\textrm{m}^2$, but its flux in the basin is almost zero. Convergence of sensible heat flux occurs from the ground surface toward the atmosphere in the lower layer, causing the layer over the mountain to be warmed up, but no convergance of the flux over the basin results from the significant mixing of air within the CBL. As horizontal transport of sensible heat flux from the top of the mountain toward over the basin results in the continuous accumulation of heat with time, enhancing air temperature at the surface of the basin, especially Taegu city to be higher than $39.3^{\circ}C$. Since latent heat fluxes are $270W/\textrm{m}^2$ near the top of the mountain and $300W/\textrm{m}^2$ along the slope of the mountain and the basin, evaporation of water vapor from the surface of the basin is much higher than one from the mountain and then, horizontal transport of latent heat flux is from the basin toward the mountain, showing relative humidity of 65 to 75% over the mountain to be much greater than 50% to 55% in the basin. At night, sensible heat fluxes have negative values of $-120W/\textrm{m}^2$ along the slope near the top of the mountain and $-50W/\textrm{m}^2$ at the surface of the basin, which indicate gain of heat from the lower atmosphere. Nighttime radiative cooling produces a shallow nocturnal surface inversion layer with a thickness of about 100m, which is much lower than common surface inversion layer, and lifts extremely heated air masses for daytime hours, namely, a warm pool of $34^{\circ}C$ to be isolated over the ground surface in the basin. As heat transfer from the warm pool in the lower atmosphere toward the ground of the basin occurs, the air near the surface of the basin does not much cool down, resulting in the persistence of high temperature at night, called nocturnal thermal high or tropical night. High relative humidity of 75% is found at the surface of the basin under the moderate wind, while slightly low relative humidity of 60% is along the eastern slope of the high mountain, due to adiabatic heating by the srong downslope wind. Air temperature near the surface of the basin with high moisture in the evening does not get lower than that during the day and the high temperature produces nocturnal warming situation.

• Journal of the Korean Society for Precision Engineering
• /
• v.30 no.5
• /
• pp.507-512
• /
• 2013

This paper proposes a new model to compensate for errors of a five-axis machine tool. A matrix with error components, that is, an error matrix, is separated from the error synthesis model of a five-axis machine tool. Based on the kinematics and inversion of the error matrix which can be obtained not by using a numerical method, an error compensation model is established and used to calculate compensation values of joint variables. The proposed compensation model does not need numerical methods to find the compensation values from the error compensation model, which includes nonlinear equations. An experiment using a double ball-bar is implemented to verify the proposed model.

• Computational Structural Engineering
• /
• v.8 no.4
• /
• pp.129-136
• /
• 1995

This paper presents a boundary element method for obtaining approximate solutions of 2-dimensional consolidation problems based on the Biot's linear theory. Laplace transform is applied to differential equation system in order to eliminate the time dependency. The boundary integral equations in transformed space are formulated and the fundamental solutions are shown in a closed form. In order to convert the transformed solutions to the ones in real space, the Hosono's numerical Laplace transform inversion method is applied. As a numerical example, a half-space consolidation problem subjected to a strip local load is selected and the applicability of the method is demonstrated through the comparison with the exact solutions.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.311-315
• /
• 2000

A numerical investigation is made of time-dependent buoyant convection in a square of a non-Boussinesq fluid. The density-temperature$({\rho}-T)$ relation is modeled by a quadratic function, with the maximum density ${\rho}_M$ at temperature $T_M$. The horizontal walls of the square are insulated, and a pulsating temperature $T_H=T_M+{\Delta}T'\;sin({\omega}{\tau})$ is imposed on the hot vertical sidewall. The temperature at the cold wall $T_c$ is constant. Extensive numerical solutions to the governing Navier-Stokes equations are portrayed. Resonance is identified by monitoring the amplitude of the mid-plane Nusselt number, $A(Nu^*)$. The primary resonance frequency is found by matching ${\omega}$ to the nondimensional basic mode $N_1$ of internal gravity oscillations. Due to the quadratic$({\rho}-T)$ relationship, the effective pulsation frequency for density, $2{\omega}$, is meaningful, which brings forth the secondary resonance frequency, i.e., $2{\omega}=N_1$

• Nuclear Engineering and Technology
• /
• v.24 no.3
• /
• pp.228-235
• /
• 1992

Some soluble species may not be solubility-limited or congruent-released with the matrix species. For example, during the operation of the nuclear reactor, the fission products can be accumulated in the fuel-cladding gap, voids, and grain boundaries of the fuel rods. In the waste package for spent-fuel placed in a geologic repository, the high solubility species of these fission products accumulated in the“gap”, e.g. cesium or iodine are expected to dissolve rapidly when ground water penetrates fuel rods. The time and space dependent mass transport for high solubility nuclides in the gap is analyzed, and its numerical illustrations are demonstrated. The approximate solution that is valid for all times is developed, and validated by comparison with an asymptotic solution and the solution obtained by the numerical inversion of Laplace transform covering the entire time span.