• Geophysics and Geophysical Exploration
• /
• v.5 no.3
• /
• pp.178-184
• /
• 2002

For economy and convenience, seismic refraction survey is widely used in surveying for large civil engineering work. The purpose of this study is to obtain the numerical responses of various models using Kirchhoff migration, and to analyze its application to the real data processing. Synthetic traveltime curve was calculated by vidale's algorithm, and various models such as 2 or 3 layer model and irregular topography model are tested to simulate the response of real structure. In order to compare the effect of initial velocity model, true velocity models, inversion results by tomography, smooth velocity models are used as an initial guess. The responses of model data show that the algorithm of this study is more sensitive to initial velocity model than the reflection survey, so choosing a suitable initial velocity model will be the most important thing in real data processing.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2023.05a
• /
• pp.424-424
• /
• 2023

댐, 저수지 등 수자원 시스템분석 시 가장 기초가 되는 유입량 자료는 실측 수위(저수량)와 방류량을 역산하여 산정된다. 이 중 댐 수위는 수표면 진동으로 인해 변동이 크며, 특히, 급격한 수위 변화가 발생하는 홍수기에는 수위-저수량 변환 시 큰 오차가 발생하여 유입량 진동이 더욱 커지게 된다. 하지만 홍수기 저수지 운영 효과 분석 등 관련 연구를 위해서는 시간 간격이 짧은 10분 또는 1시간 단위의 유입량 자료가 필요함에 따라 관련 연구 수행 시 이동평균법(Moving Average) 등을 통해 실측 유입량 자료를 보정하여 사용하는 것이 일반적이다. 데이터 평활화를 위해 이동평균법을 적용하면 데이터의 변동을 효과적으로 줄일 수는 있지만 실측자료와 비교하였을 때 첨두 유입량이 큰 폭으로 감소하거나, 첨두 유입량 발생시간이 지체되는 문제가 발생한다. 본 연구에서는 저수지 유입량과 같이 변동이 큰 수문자료의 평활화를 위해 가우시안 가중 이동평균법(Gaussian-weighted moving average technique), 사비츠키-골레이 필터링기법(Savitzky-Golay filtering technique) 등 필터링 기법을 댐 유입량 보정에 적용하고, 이에 따른 효과를 분석하고자 하였다. 이를 위해 2020년 8월에 발생한 홍수사상을 대상으로 충주댐, 합천댐 등 다목적댐 유입량 자료를 수집하고, 보정을 수행하였다. 필터링 기법의 적용 효과 분석을 위해서는 실측자료와 이동평균법을 적용하여 보정한 결과와 비교하였고, 추가적으로 비교적 변동이 작은 일 단위 유입량 자료와의 양적 비교를 진행하였다. 그 결과 이동평균법을 적용하였을 때보다 필터링 기법을 적용하였을 때 실측자료와의 양적 차이가 작고, 첨두 유입량 및 첨두 유입 발생시간에서도 차이를 큰 폭으로 감소시킬 수 있는 것으로 확인되었다.

• Geophysics and Geophysical Exploration
• /
• v.16 no.1
• /
• pp.18-26
• /
• 2013

Shear-wave velocity ($v_s$) structures beneath two seismic stations, JJU and JJB on the flanks of the volcano Halla on Jeju island, Korea, were estimated by receiver-function inversion and H-${\kappa}$ stacking applied to 150 teleseismic events ($M_W{\geq}5.5$) recorded since 2007. $P_S$ waves converted at the Moho discontinuity does not appear clearly for northwesterly back-azimuths ($207{\sim}409^{\circ}$, average $308^{\circ}$) at station JJU and southeasterly back-azimuths ($119{\sim}207^{\circ}C$, average $163^{\circ}$) at station JJB. This may be due to a gradual velocity increase at Moho or heterogeneity within the crust. The $v_s$ models derived by inversion of receiver functions indicate a distinct low velocity layer ($v_s{\leq}3.5km/s$; LVL) within the crust and a gradual increase in $v_s$ in the depth interval of 30 to 40 km. Within the radius of 18 km beneath station JJB, the LVL occurs at depths of 14 ~ 26 km and the 'Moho' ($v_s{\geq}4.3km/s$) is at 34 km depth. Ten kilometers to the west, within the radius of 16 km beneath station JJU, both the LVL and the Moho are significantly shallower, at depths of 14 to 24 km and 30 km, respectively. H-${\kappa}$ analyses for stations JJU and JJB yield estimated crustal thickness of 29 and 33 km and $v_p/v_s$ ratios of 1.64 and 1.75, respectively. The lesser $v_p/v_s$ ratio was derived for rocks nearest to th peak of the volcano.

• Journal of Energy Engineering
• /
• v.19 no.3
• /
• pp.143-150
• /
• 2010

Line source method of borehole system assumes the entire surrounding medium is uniform. However, thermal properties of grouting region are considerably different from those of surrounding soil. In this study we investigate the effect of grouting materials on the solution of line source method with the aid of numerical analysis. This numerical model generates the temperature of borehole fluid with which line source solution can be obtained. Then this solution can be compared with input condition of numerical model. The results of this comparison show that thermal conductivity and borehole thermal resistance of line source solution are approximately 86% and 91% of the input condition of numerical model. Chart method is developed in this study to find the numerical input conditions (thermal conductivity and borehole thermal resistance) from the line source solution. Thermal response test of test borehole is conducted, the results of which are approximately consistent with the Chart method. Thermal property changes of grouting materials on the line source solution are also examined.

• Geophysics and Geophysical Exploration
• /
• v.13 no.4
• /
• pp.397-406
• /
• 2010

Resistivity method is a practical and effective geophysical technique to detect leakage zones in embankment dams. Generally, resistivity survey conducted along the crest assumes that the embankment dam has a 2D structure. However, the 3D topography of the embankment distorts significantly resistivity data measured on anywhere of the dam. This study evaluates the influence from 3D effects created by specific dam geometry and effects of water level fluctuations through the 3D finite element modeling technique. Also, a comparison between different locations of survey line are carried out, and topographic correction technique is developed for the resistivity data obtained along the embankment dam. Furthermore, using synthetic resistivity data for an embankment dam model with leakage zone, detectability of leakage zones is estimated through 2.5D inversion.

• Journal of the Korean Geophysical Society
• /
• v.5 no.1
• /
• pp.63-71
• /
• 2002

The integral equation method is a powerful tool for electromagnetic numerical modeling. But the difficulty of this technique is the size of their linear equations, which demands excessive memory and calculation time to invert. This limitation of the integral equation method becomes critical in inverse problem. To overcome this limitation, a lot of approximation and series methods, such as conventional Born, modifed Born and extended Born, were developed. But all the methods need volume integration of Green tensor, which is very time consuming. In electromagnetic theory, Green tensor rapidly decreases as the distance between source and field cell increases. Therefore, the source cell which are far away from the field cell does not make an effect on the electric field of the field cell. Consequently, by ignoring the effect of Green tensor due to far away source cells, computing time for electromagnetic numerical modeling can be reduced dramatically. Comparisons of this new method against a full integral equation, extended Born approximation and series code show that the method is accurate enough much less time consuming.

• Geophysics and Geophysical Exploration
• /
• v.4 no.4
• /
• pp.124-132
• /
• 2001

Electrical resistivity survey has been conducted for delineating geological fault structure in Kaekok-ri near Kyungju. In general, electrical resistivity survey has an advantage of searching buried faults and its traces compared with other geophysical survey methods. Distribution of electrical conductivity in the ground is influenced by the ratio of pores, groundwater and clay minerals. These properties are evidenced indirectly to explain for weathering condition, faults and fracture Bones. Thus the electrical resistivity survey can be an effective method to find buried faults. We have carried out two dimensional (2-D) interpretation by means of smoothness-constrained least-squares and finite element method. Field data used in this paper was acquired at Kaekok-ri, Wuedong-eup, Kyungju-si, where is Ulsan Fault and is close to the region in which debatable quaternary fault traces were found recently. The dipole-dipole array resistivity survey which could show the 2-D subsurface electrical resistivity structure, was carried out in the area with three lines. The results showed good property of fault, fracture zone and fault traces which we estimated were congruous with the results. Through this study, 2-D electrical resistivity survey interpretation for fault is useful to apply.

• Geophysics and Geophysical Exploration
• /
• v.15 no.1
• /
• pp.1-7
• /
• 2012

To evaluate methods of determining near-surface shear-wave velocities (${\nu}_s$), we derived dispersion curves of Rayleigh waves generated by both passive and active sources in Chuncheon, Korea. Microtremors were recorded for 5 minutes in each of four triangular arrays with radii of 5 ~ 40 m. Those data were analyzed using the Spatial Autocorrelation method. Rayleigh waves were also generated by a hammer source and recorded in the same area for 2 s using 24 4.5-Hz geophones. Multichannel Analysis of Surface Waves was applied to those data. Velocity spectra were derived with relatively high signal-to-noise ratios in the frequency ranges of 7 ~ 19 and 11 ~ 50 Hz for the microtremors and synthetically generated Rayleigh waves, respectively. The resultant dispersion curves were combined as one and then input to inversion to derive shear wave velocities that were compared with a lithology log from a nearby well. Shearwave velocities in the top soil and soft-rock layers are almost constant with values of 221 and 846 m/s, respectively; while the inverse-modeled ${\nu}_s$ increases linearly in the gravelly sand, cobbles, and weathered-rock layers. If rock type is classified based on shear-wave velocity, the inversion-derived boundary between weathered-rock and soft rock may be about 5 m deeper than in the well log.

• Resources Recycling
• /
• v.29 no.3
• /
• pp.11-23
• /
• 2020

This study aims to investigate breakage and liberation characteristics of iron ore from Shinyemi mine, Jeongseon by ball mill. Parameters of breakage functions for three grade samples of iron ore were obtained using single-sized-feed breakage test and back-calculation based on nonlinear programming. The results showed that with the increase in the grade of iron ore, the breakage rate factor decrease whereas the particle size sensitivity decreases. This results from retardation of microcrack-propagation by magnetite grain in the ore. Breakage distribution analysis showed that the breakage mechanism appear to be impact fracture dominant with the increase of grade owing to the stress distribution effect by magnetite grain. Degree of liberation (DOL) increased with the increase in grade and decrease in particle size, respectively. Using the breakage function and size-DOL relationship, a model that can predict time-dependent-DOL is established. When scale-up factors from operating condition are available, the model is expected to be capable of predicting size and DOL with time in actual mining process.

• Geophysics and Geophysical Exploration
• /
• v.11 no.1
• /
• pp.60-67
• /
• 2008

In this paper we consider an application of the method of electromagnetic (EM) migration to the interpretation of a typical marine controlled-source (MCSEM) survey consisting of a set of sea-bottom receivers and a moving electrical bipole transmitter. Three-dimensional interpretation of MCSEM data is a very challenging problem because of the enormous number of computations required in the case of the multi-transmitter and multi-receiver data acquisition systems used in these surveys. At the same time, we demonstrate that the MCSEM surveys with their dense system of transmitters and receivers are extremely well suited for application of the migration method. In order to speed up the computation of the migration field, we apply a fast form of integral equation (IE) solution based on the multigrid quasi-linear (MGQL) approximation which we have developed. The principles of migration imaging formulated in this paper are tested on a typical model of a sea-bottom petroleum reservoir.