• Proceedings of the Korea Water Resources Association Conference
• /
• 2022.05a
• /
• pp.299-299
• /
• 2022

수문학적 가뭄 발생의 직접적 영향은 강수부족량이나, 다양한 사회경제적 인자들은 수문학적 가뭄에 간접적으로 영향을 미치고 있다. 물관리 선진기관에서는 인간의 활동 및 물관리 방식에 따라 수문학적 가뭄을 심화시키거나 완화시킬 수 있음을 인지하고, 인간의 물사용이 가뭄에 미치는 영향을 평가하기 위한 다양한 연구가 이루어지고 있다. 본 연구에서는 강수량 및 미래의 인구변화에 따른 수문학적 가뭄의 영향의 정도를 판단함으로써, 인간의 활동이 가뭄에 미치는 영향을 정량적으로 제시하고자 한다. 충정북도 시군지역을 대상지역으로 선정하였으며, 시군 장래인구 추정값을 미래 인구자료로, 미래 유출량이 산정되어 제공되는 RCP 4.5와 RCP 8.5시나리오를 활용하여 미래 가뭄상황 예측하였다. 강수량 및 인구변화가 수문학적 가뭄에 미치는 영향 평가를 위하여 코플라함수 기반의 베이지안 네트워크 모형이 활용하였다. 베이지안 네트워크는 강수량, 인구밀도, 수문학적 가뭄사이의 관계 도출을 위하여 활용되었으며, 베이지안 네트워크 내의 결합확률의 산정을 위하여 코플라 함수가 활용되었다. 미래의 강수량 및 인구밀도의 변화에 따른 수문학적 가뭄의 영향 관계를 분석한 결과는 다음과 같다. 강수량이 인구밀도보다 수문학적 가뭄의 발생에 영향을 미치며, 약 0.2~0.3 정도 발생확률이 크게 산정되었다. 두 인자를 동시에 고려할 경우, 강수량이 적고, 인구밀도가 높아지는 조건(F(강수량)=0.1, F(인구밀도)=0.9)에서는 조건부 CDF 변화율이 크게 나타나, 곧 수문학적 가뭄의 위험성이 높음을 확인할 수 있었다. 인구밀도는 수문학적 가뭄의 발생 위험성을 높이 알려져 있으나, 정량적으로 그 값을 제시한 연구 사례는 찾기 어렵다. 이에 따라 본 연구에서는 가뭄의 영향정도를 정량적으로 표현하였으며, 한 인자만의 영향이 아닌 두 개 이상의 인자들의 복합적인 영향 정도를 제시함으로써 수치적인 비교가 가능하게 하였다. 미래 추정 인자가 인구자료가 한정적이라 인구 자료만을 활용하여 수문학적 가뭄에 미치는 영향을 분석하였으나, 다른 사회경제적 지표를 활용하여 미래 변화에 따른 미래 수문학적 가뭄의 영향 정도의 비교 및 분석 결과를 바탕으로 가뭄 대응 우선순위 선정을 위한 연구자료로 활용 가능할 것으로 사료된다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2022.05a
• /
• pp.6-6
• /
• 2022

최근 수재해에 대응하기 위한 물관리 환경은 기후변화에 따른 홍수 피해 심화와 댐과 하천 시설의 노후화 점증, 하천관리일원화 등 정책적 변화, 그리고 포스트코로나 디지털 혁신 등 복합적 대전환 시대 진입에 따라 복잡다단한 양상을 보이고 있다. 디지털 트윈은 디지털 대전환(digital transformation) 시대 다양한 산업 영역에서 지능화와 생산성 향상을 목적으로 도입되고 있다. 본 국가 시범사업에서는 170 km에 달하는 섬진강 유역 전체를 대상으로 홍수에 대응하기 위한 디지털 트윈 플랫폼(K-Twin SJ)을 구축하고 있다. 본 플랫폼은 국가 인프라 지능정보화 사업의 일환으로 시작되었으며, 공간정보와 시설물 모델링, 홍수 분석 등 수재해에 대응하기 위한 수자원 분야의 다학제적인 강소기업들과 K-water에서 컨소시엄을 구성하여 추진하고 있다. 본 사업의 내용은 섬진강 댐-하천 유역에 대하여 고정밀도 3D 공간정보화, 실시간 물관리 데이터 연계, 홍수 분석 시뮬레이션, AI 댐 운영 최적화, AI 사면 정보 생성, 하천 제방 안전성 평가, AI 지능형 CCTV 영상분석, 간이 침수피해 예측, 드론 제약사항 조사 체계 개발을 포함하고 있다. 물관리 데이터와 하천 시설정보를 트윈 플랫폼 상에서 위치기반으로 시각화 표출하기 위해서는 유역의 공간정보를 3차원으로 구축하는 과정이 필수적이다. 따라서 GIS 기반의 섬진강 하천 중심 공간정보 구축을 위해 유역의 국가 정사영상과 5m 수치표고모형(DEM)은 최신성과를 협조 받아 적용하였으며, 홍수 분석을 위한 하천 중심 공간정보는 신규 헬기에 LiDAR 매핑을 수행하여 0.5m 급 DEM을 신규 구축하였다. 또한 하천 시설물 중 섬진강댐과 79개 주요 하천 횡단 교량과 3개 보 시설을 지상기준점 측량과 드론 매핑, 패턴 방식의 경량화 작업을 통해 트윈에 탑재할 수 있는 시설물 3D 객체 모델을 제작하였다. 홍수 분석을 위해서는 섬진강 유역에 대해 K-Drum, K-River, K-Flood 모델을 구축하였으며, AI 하천 수위 예측 학습 모델을 개발하였다. 섬진강 디지털 트윈 유역 물관리 플랫폼을 통해 데이터 기반의 똑똑한 물관리를 구현하고자 한다.

• Journal of the Korean Geotechnical Society
• /
• v.40 no.2
• /
• pp.81-93
• /
• 2024

In Korea, numerous large-scale infrastructure construction projects and housing site developments are being undertaken. However, due to limited land availability, sourcing high-quality backfill materials that meet the standards for railroad and road embankment compaction and mechanically stabilized earth (MSE) retaining wall construction poses significant challenges. Concurrently, there has been an increase in structural failures of many MSE retaining walls, attributed primarily to reduced bearing capacity and impaired drainage performance, resulting from inadequate backfill compaction. This study aimed to analyze the structural performance and safety of an MSE retaining wall using recycled soil as backfill. We conducted small-scale model tests utilizing 3D printing technology combined with two-dimensional numerical analysis. The study quantitatively evaluated the MSE retaining wall's performance concerning the recycled soil mixing ratio and reinforcement installation methods. Furthermore, the utility of 3D printing was confirmed through the production of an experimental wall designed to facilitate easy reinforcement attachment, mirroring the conditions of actual MSE retaining wall construction.

• Journal of Korea Water Resources Association
• /
• v.50 no.7
• /
• pp.455-465
• /
• 2017

Understanding of the two-dimensional velocity field is crucial in terms of analyzing various hydrodynamic and fluvial processes in the riverine environments. Until recently, many numerical models have played major roles of providing such velocity field instead of in-situ flow measurements, because there were limitations in instruments and methodologies suitable for efficiently measuring in the broad range of river reaches. In the last decades, however, the advent of modernized instrumentations started to revolutionize the flow measurements. Among others, acoustic Doppler current profilers (ADCPs) became very promising especially for accurately assessing streamflow discharge, and they are also able to provide the detailed velocity field very efficiently. Thus it became possible to capture the velocity field only with field observations. Since most of ADCPs measurements have been mostly conducted in the cross-sectional lines despite their capabilities, it is still required to apply appropriate interpolation methods to obtain dense velocity field as likely as results from numerical simulations. However, anisotropic nature of the meandering river channel could have brought in the difficulties for applying simple spatial interpolation methods for handling dynamic flow velocity vector, since the flow direction continuously changes over the curvature of the channel shape. Without considering anisotropic characteristics in terms of the meandering, therefore, conventional interpolation methods such as IDW and Kriging possibly lead to erroneous results, when they dealt with velocity vectors in the meandering channel. Based on the consecutive ADCP cross-sectional measurements in the meandering river channel. For this purpose, the geographic coordinate with the measured ADCP velocity was converted from the conventional Cartesian coordinate (x, y) to a curvilinear coordinate (s, n). The results from application of A-VIM showed significant improvement in accuracy as much as 41.5% in RMSE.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.30-36
• /
• 2006

Japan's first pilot-scale $CO_2$ sequestration experiment has been conducted in Nagaoka, where 10400 t of $CO_2$ have been injected in an onshore aquifer at a depth of about 1100 m. Among various measurements conducted at the site for monitoring the injected $CO_2$, we conducted time-lapse crosswell seismic tomography between two observation wells to determine the distribution of $CO_2$ in the aquifer by the change of P-wave velocities. This paper reports the results of the crosswell seismic tomography conducted at the site. The crosswell seismic tomography measurements were carried out three times; once before the injection as a baseline survey, and twice during the injection as monitoring surveys. The velocity tomograms resulting from the monitoring surveys were compared to the baseline survey tomogram, and velocity difference tomograms were generated. The velocity difference tomograms showed that velocity had decreased in a part of the aquifer around the injection well, where the injected $CO_2$ was supposed to be distributed. We also found that the area in which velocity had decreased was expanding in the formation up-dip direction, as increasing amounts of $CO_2$ were injected. The maximum velocity reductions observed were 3.0% after 3200 t of $CO_2$ had been injected, and 3.5% after injection of 6200 t of $CO_2$. Although seismic tomography could map the area of velocity decrease due to $CO_2$ injection, we observed some contradictions with the results of time-lapse sonic logging, and with the geological condition of the cap rock. To investigate these contradictions, we conducted numerical experiments simulating the test site. As a result, we found that part of the velocity distribution displayed in the tomograms was affected by artefacts or ghosts caused by the source-receiver geometry for the crosswell tomography in this particular site. The maximum velocity decrease obtained by tomography (3.5%) was much smaller than that observed by sonic logging (more than 20%). The numerical experiment results showed that only 5.5% velocity reduction might be observed, although the model was given a 20% velocity reduction zone. Judging from this result, the actual velocity reduction can be more than 3.5%, the value we obtained from the field data reconstruction. Further studies are needed to obtain more accurate velocity values that are comparable to those obtained by sonic logging.

• Journal of the Korean Geotechnical Society
• /
• v.32 no.10
• /
• pp.41-53
• /
• 2016

Stability of the braced earth wall in the composite ground, which is composed of the jointed base rocks and the soil strata depends on the earth pressure acting on it. In most cases, the earth pressure is calculated by the empirical method, in which base rocks are considered as a soil strata with the shear strength parameters of base rocks. In this case the effect of the joint dips of the jointed base rocks is ignored. Therefore, the calculated earth pressure is smaller than the actual earth pressure. In this study, the magnitude and the distribution of the earth pressure acting on the braced wall in the composite ground depending on the joint dips of the base rocks and the ratio of soil strata and base rocks were experimentally studied. Two dimensional large-scale model tests were conducted in a large scale test facility (height 3.0 m, length 3.0 m and width 0.5 m) by installing 10 supports in a scale of 1/14.5. The test ground was presumed with the base rock ratio of the composite ground of 65%:35% and 50%:50% and with the joint dips for each base rock layer, $0^{\circ}$, $30^{\circ}$, $45^{\circ}$ and $60^{\circ}$, respectively. And then finite element analyses were performed in the same condition. As results, the earth pressure on the braced wall increased as the base rock layer's joint dips became larger. And earth pressure at the rock layer increased as the rock rate became larger. The largest earth pressure was measured when the base rock rate was 50% (R50) and the rock layer's joint dips was $60^{\circ}$. Based on these results, a formular for the calculation of the earth pressure in the composite ground could be suggested. Distribution of earth pressure was idealized in a quadrangular form, in which the magnitude and the position of peak earth pressure depended on the rock ratio and the joint dips.

• Geophysics and Geophysical Exploration
• /
• v.5 no.4
• /
• pp.280-290
• /
• 2002

Three-dimensional (3-D) electromagnetic (EM) modeling algorithm has been developed using finite element method (FEM) to acquire more efficient interpretation techniques of EM data. When FEM based on nodal elements is applied to EM problem, spurious solutions, so called 'vector parasite', are occurred due to the discontinuity of normal electric fields and may lead the completely erroneous results. Among the methods curing the spurious problem, this study adopts vector element of which basis function has the amplitude and direction. To reduce computational cost and required core memory, complex bi-conjugate gradient (CBCG) method is applied to solving complex symmetric matrix of FEM and point Jacobi method is used to accelerate convergence rate. To verify the developed 3-D EM modeling algorithm, its electric and magnetic field for a layered-earth model are compared with those of layered-earth solution. As we expected, the vector based FEM developed in this study does not cause ny vector parasite problem, while conventional nodal based FEM causes lots of errors due to the discontinuity of field variables. For testing the applicability to high frequencies 100 MHz is used as an operating frequency for the layer structure. Modeled fields calculated from developed code are also well matched with the layered-earth ones for a model with dielectric anomaly as well as conductive anomaly. In a vertical electric dipole source case, however, the discontinuity of field variables causes the conventional nodal based FEM to include a lot of errors due to the vector parasite. Even for the case, the vector based FEM gave almost the same results as the layered-earth solution. The magnetic fields induced by a dielectric anomaly at high frequencies show unique behaviors different from those by a conductive anomaly. Since our 3-D EM modeling code can reflect the effect from a dielectric anomaly as well as a conductive anomaly, it may be a groundwork not only to apply high frequency EM method to the field survey but also to analyze the fold data obtained by high frequency EM method.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.26 no.1
• /
• pp.59-78
• /
• 2024

In tunnel construction, the stability is evaluated by the settlement of adjacent structures and ground, but the shear strain of the ground is the main factor that determines the failure mechanism of the ground due to the tunnel excavation and the change of the operating load, and can be used to review the stability of the tunnel excavation and to calculate the reinforcement area. In this study, a twin tunnel excavation was simulated on a soft ground in an urban area through a laboratory model test to analyze the behavior of the twin tunnel excavation on the adjacent pile grouped foundation and adjacent ground. Both the displacement and the shear strain of ground were obtained using a close-range photogrammetry during laboratory model test. In addition, two-dimensional finite element numerical analysis was performed based on the model test. The results of a back-analysis showed that the maximum shear strain rate tends to decrease as the horizontal distance between the pillars of the twin tunnel and the vertical distance between the toe of the pile group and the crown of the tunnel were decreased. The impact of the second tunnel on the first tunnel and pile group was decreased as the horizontal distance between the pillars of the twin tunnel was increased. In addition, the vertical distance between the toe of the pile group and the crown of the tunnel had a relatively greater impact on the shear strain results than the horizontal distance of the pillars between the twin tunnels. According to the results of the close-range photogrammetry and numerical analysis, the settlement of adjacent pile group and adjacent ground was measured within the design criteria, but the shear strain of the ground was judged to be outside the range of small strain in all cases and required reinforcement.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.25 no.2
• /
• pp.63-75
• /
• 2013

The Inner Port Phase 2 area of the Pyeongtaek-Dangjin Port is enclosed by a total of three permeable sea-walls, and the disposal site to the east of the Inner Port Phase 2 is also enclosed by two permeable sea-walls. The maximum tidal range measured in the Inner Port Phase 2 and in the disposal site in May 2010 is 4.70 and 2.32 m, respectively. It reaches up to 54 and 27%, respectively of 8.74 m measured simultaneously in the exterior. Regression formulas between the difference of hydraulic head and the rate of interior water volume change, are induced. A three-dimensional numerical hydrodynamic model for the Asan Bay is constructed incorporating a module to compute water discharge through the permeable sea-walls at each computation time step by employing the formulas. Hydrodynamics for the period from 13th to 27th May, 2010 is simulated by driving forces of real-time reconstructed tide with major five constituents($M_2$, $S_2$, $K_1$, $O_1$ and $N_2$) and freshwater discharges from Asan, Sapkyo, Namyang and Seokmoon Sea dikes. The skill scores of modeled mean high waters, mean sea levels and mean low waters are excellent to be 96 to 100% in the interior of permeable sea-walls. Compared with the results of simulation to obstruct the flow through the permeable sea-walls, the maximum current speed increases by 0.05 to 0.10 m/s along the main channel and by 0.1 to 0.2 m/s locally in the exterior of the Outer Sea-wall of Inner Port. The maximum bottom shear stress is also intensified by 0.1 to 0.4 $N/m^2$ in the main channel and by more than 0.4 $N/m^2$ locally around the arched Outer Sea-wall. The module developed to compute the flow through impermeable seawalls can be practically applied to simulate and predict the advection and dispersion of materials, the erosion or deposion of sediments, and the local scouring around coastal structures where large-scale permeable sea-walls are maintained.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.17 no.4
• /
• pp.384-398
• /
• 2015

In this paper, the high-resolution Weather Research and Forecasting/Noah-MultiParameterization (WRF/Noah-MP) modeling system is configured for the Cheongmicheon Farmland site in Korea (CFK), and its performance in land and atmospheric simulation is evaluated using the observed data at CFK during the 2014 special observation period (21 August-10 September). In order to explore the usefulness of turning on Noah-MP dynamic vegetation in midterm simulations of surface and atmospheric variables, two numerical experiments are conducted without dynamic vegetation and with dynamic vegetation (referred to as CTL and DVG experiments, respectively). The main results are as following. 1) CTL showed a tendency of overestimating daytime net shortwave radiation, thereby surface heat fluxes and Bowen ratio. The CTL experiment showed reasonable magnitudes and timing of air temperature at 2 m and 10 m; especially the small error in simulating minimum air temperature showed high potential for predicting frost and leaf wetness duration. The CTL experiment overestimated 10-m wind and precipitation, but the beginning and ending time of precipitation were well captured. 2) When the dynamic vegetation was turned on, the WRF/Noah-MP system showed more realistic values of leaf area index (LAI), net shortwave radiation, surface heat fluxes, Bowen ratio, air temperature, wind and precipitation. The DVG experiment, where LAI is a prognostic variable, produced larger LAI than CTL, and the larger LAI showed better agreement with the observed. The simulated Bowen ratio got closer to the observed ratio, indicating reasonable surface energy partition. The DVG experiment showed patterns similar to CTL, with differences for maximum air temperature. Both experiments showed faster rising of 10-m air temperature during the morning growth hours, presumably due to the rapid growth of daytime mixed layers in the Yonsei University (YSU) boundary layer scheme. The DVG experiment decreased errors in simulating 10-m wind and precipitation. 3) As horizontal resolution increases, the models did not show practical improvement in simulation performance for surface fluxes, air temperature, wind and precipitation, and required three-dimensional observation for more agricultural land spots as well as consistency in model topography and land cover data.