• Journal of the Korean Institute of Gas
• /
• v.16 no.3
• /
• pp.35-41
• /
• 2012

Manifold researches for carbon capture and storage (CCS) have been developed and large scale-carbon capture system can be performed recently. Hence, the technologies for $CO_2$ sequestration or storage become necessary to handle the captured $CO_2$. Among them, enhanced oil recovery using $CO_2$ can be a solution since it guarantees both oil recovery and $CO_2$ sequestration. In this study, the miscible flow of oil and $CO_2$ in porous media is modeled to analyze the effect of enhanced oil recovery and $CO_2$ sequestration. Based on Darcy-Muskat law, the equation is modified to consider miscibility of oil and $CO_2$ and the change of viscosity. Finite volume method is used for numerical modeling. As results, the pressure and oil saturation changes with time can be predicted when oil, water, and $CO_2$ are injected, respectively, and $CO_2$ injection is more efficient than water injection for oil recovery.

• Journal of Ocean Engineering and Technology
• /
• v.17 no.1
• /
• pp.8-15
• /
• 2003

본 논문에서는 투과성 유연막이 수면밑 일정한 깊이에 수평으로 잠겨있을 때 투과성 유연막에 의한 파랑제어성능을 살펴보았다. 해석 방법으로는 유체문제는 고유함수전개법 (Eigenfunction expansion method)을 사용하였고, 유연막과 파랑의 상호작용문제는 Newmann 이 제시한 유탄성 이론 (hydro-elastic theory)을 채택하였다. 막의 투과성 효과를 고려하기 위하여 수평막에서의 수직속도는 수평막 상하의 압력차에 선형적으로 비례하며 그들 사이에는 위상차가 없다고 가정한 Darcy 법칙을 사용하였다. 투과성 수평막의 설계변수 (초기장력, 길이, 잠긴 깊이, 공극율)와 입사파의 주파수를 바꿔가면서 반사율과 투과율 그리고 에너지 손실율을 살펴보았다.

• The Journal of Engineering Geology
• /
• v.27 no.4
• /
• pp.377-382
• /
• 2017

Field experiments were conducted to estimate streambed hydraulic conductivity at 15 sites in the Anseongcheon watershed, Korea. Seepage meters and piezometers were installed in the streambed at each site to measure the amount of stream water-groundwater exchange and the hydraulic gradient. The vertical hydraulic conductivity was then calculated using Darcy's formula. The measured stream water-groundwater exchange rates were $4.08{\times}10^{-6}$ to $1.49{\times}10^{-5}m/s$, and the vertical hydraulic gradients were 0.005 to 0.145. The data suggest the streambed hydraulic conductivity to be $7.80{\times}10^{-5}$ to $1.58{\times}10^{-3}m/s$. The results show significant differences in connectivity between stream and aquifer. Quantification of the hydraulic interconnection between stream and aquifer, and evaluation of the effects of groundwater development and utilization on the streamflow require hydrogeological investigations of the connection between stream and aquifer, including the hydraulic conductivity of the streambed. Various field testing and analysis methods for hydrogeological assessment also require further improvement.

• Journal of Korea Water Resources Association
• /
• v.45 no.1
• /
• pp.39-52
• /
• 2012

Surface-subsurface interactions are an intrinsic component of the hydrologic response within a watershed. In general, these interactions are considered to be one of the most difficult areas of the discipline, particularly for the modeler who intends simulate the dynamic relations between these two major domains of the hydrological cycle. In essence, one major complexity is the spatial and temporal variations in the dynamically interacting system behavior. The proper simulation of these variations requires the need for providing an appropriate coupling mechanism between the surface and subsurface components of the system. In this study, an approach for modelling surface-subsurface flow and transport in a fully intergrated way is presented. The model uses the 2-dimensional diffusion wave equation for sheet surface water flow, and the Boussinesq equation with the Darcy's law and Dupuit-Forchheimer's assumption for variably saturated subsurface water flow. The coupled system of equations governing surface and subsurface flows is discretized using the finite volume method with central differencing in space and the Crank-Nicolson method in time. The interactions between surface and subsurface flows are considered mass balance based on the continuity conditions of pressure head and exchange flux. The major module consists of four sub-module (SUBFA, SFA, IA and NS module) is developed.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2010.05a
• /
• pp.223-227
• /
• 2010

하천의 수위와 유량에 대한 정확한 정보는 이수, 치수와 같은 수자원 관리에 있어서 가장 기본적인 기초자료가 된다. 하천 수위와 유량 자료를 확보하기 위한 지속적인 직접 계측 방법은 많은 시간과 비용이 소모되며, 홍수 시에는 위험성이 존재하여 자료 확보가 불가능하다. 따라서 수치모형을 이용한 수위-유량 자료의 예측과 활용이 필요하며, 2차원 수치모형의 경계조건을 설정할 경우에도 활용할 필요성이 있다. 본 연구에서는 복단면 개수로 및 불규칙한 하상을 보이는 횡단면 상에서의 단위유량 예측을 위한 유한요소모형을 개발한다. 지배방정식은 Wark 등 (1990)이 제시한 운동량방정식을 이용하며, 단면형상과 Manning 조도계수, 그리고 수위를 알면, 결과적으로 흐름방향 단위유량의 횡방향 분포를 얻을 수 있다. 개발된 모형의 검증을 위해 실측 자료와 비교하며, 또한 Darcy-Weisbach 마찰계수를 이용하는 Darby와 Thorne (1996)의 모형 결과와도 비교한다. 검증된 모형의 알고리즘을 2차원 모형의 상류단 경계조건 설정에 활용하여, 절점별 유입유량을 차등분배 시켰을 때와 그렇지 못했을 때의 결과를 비교한다. 또한 수위를 경계조건으로 입력하였을 때 개발 모형을 활용하여 유입유량을 예측하는 활용방안을 제시한다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.386-386
• /
• 2017

지하수 해석을 위해 3차원 수치해석 모형이 적용되고 있으나, 다양한 매개변수를 명확히 적용하기에는 한계가 있다. 특히, 수리전도도는 지층의 투수성정도를 나타내는 계수로 지하수 분석에 매우 중요한 매개변수이다. 신뢰성의 확보를 위해 양수시험을 통해 산출된 결과를 이용하고 있으나 소수관정의 시험결과를 유역의 대푯값으로 적용하기에 불확실도가 매우 높고, 수위변화가 독특한 지역의 정확한 수리특성을 적용하기에 한계가 있다. 이 연구에서는 수리전도도 특성을 해석하기 위하여 3차원 수치해석모형을 적용하였으며, 모델보정 방법 중 Regularization(정규화)라고 불리는 Pilot point 기법을 사용하였다. 제주도와 같이 수리전도도 값이 지역별로 차이가 크고 동일 유역 내에도 다른 지질구조를 보이는 등 동일 매질에서 동일 투수성을 보이지 않는 다양함으로 실측값들을 적용하기에 어려운 곳에서 정규화라는 보정방법은 최적화된 방법이다. 지하수위는 중제주유역 내 위치한 12개소 수위관측정의 2016년 연평균수위를 적용하였다. 미계측지역은 제주도 등수위선자료를 이용하여 DEM을 구축하였으며, 임의지점 17개소를 선정하여 대표수위로 적용하였다. 중제주유역의 평균 수리전도도는 82.90 m/day로 분석되었으며, 유역의 동측 하류부는 최대 1,745 m/day로 비교적 큰 결과가 산출되었다. 유역의 중앙에 위치하는 OR관측정을 기준으로 동 서지역의 지하수위를 검토한 결과 서측은 지형구배에 따라 지하수위가 형성하고 있으나, 동측의 경우 상 하류의 표고차가 30m이상 발생되지만 지하수위는 유사한 형태를 보이고 있다. 지하수 흐름에 해석되는 Darcy의 방정식은 수리전도도와 동수경사는 반비례관계를 나타내며, 이 이론에 의해 상 하류 지하수위가 유사하게 형성되는 동측지역은 국부적으로 수리전도도가 높게 형성되는 것으로 확인되었다. 실무에서는 유역경계에 따라 평균화된 매개변수가 적용되므로 명확한 지하수 해석이 어렵고, 수리전도도와 같이 지역적 편차가 심한 매개변수는 다양한 연구를 통해 적용성검토가 수행되어야 한다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.3 no.4
• /
• pp.59-70
• /
• 1983

In this paper, it deals with runoff analysis in a mountainous watershed where the concentration time of the river is much Shorter than that of slopes. The momentum and continuity equation are used as the basic expressions and they are solved by the kinematic wave theory. By the introducing the distribution of concentration time, representing the characteristics of topographical and hydraulic factors, into a solution. A response function is derived. Using the optimization technique (simplex Method), the response functions derived in this paper are tested by comparing the observed and the estimated values, and shows promising.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.82-82
• /
• 2015

해안 대수층은 해수와 담수가 공존하는 지역으로 상대적으로 밀도가 큰 해수가 대수층의 담수 아래에 쐐기형태로 존재하게 된다. 이러한 쐐기형태의 해수와 담수의 경계면은 압력경도의 평형에 의해 경계면이 유지되며, 해수면 또는 지하수위가 변동할 경우 해수-담수 경계면의 균형이 무너짐과 더불어 압력경도의 평행이 이루어질 때 까지 해수-담수 경계면의 이동이 계속 진행된다. 수위 변화의 주요 원인으로는 지구온난화 및 기후변화로 인한 지속적인 해수면 상승과 도서지역의 인구증가 및 산업화로 인한 무분별한 지하수의 사용 등에 의한 지하수위 저하 등을 꼽을 수 있다. 이와 같은 원인으로 해안 및 도서지역에서는 해안 대수층의 해수침투거리가 증가하여 지하수 이용에 큰 어려움을 겪고 있다. 이에 해안 대수층의 해수침투 범위 및 거리를 추정하기 위한 많은 연구들이 다양한 분야에서 지속해서 이루어지고 있지만, 서로 밀도가 다른 해수와 담수가 공존하는 해안 대수층 내의 수리특성을 명확히 파악하기에는 아직까지 미흡한 점들이 많다. 과거에는 Darcy의 법칙 및 Ghyben-Herzberg 식에 근거한 이론적인 연구들이 주로 이루어졌고, 근래에 현장관측이나 수리모형실험이 국내 외적으로 수행되고 있으나, 모든 영역의 지하수의 특성을 조사하는 것이 사실상 불가능하다. 이에 최근에는 컴퓨터 성능의 비약적인 발전과 더불어 다양한 수치해석방법에 의한 수치모델들이 개발되어 시뮬레이션에 적용되고 있다. 하지만 거의 대부분의 수치모델은 해안 대수층 수리특성을 투수계수에 의존하고 있을 뿐, 대수층 내부의 해수-담수에 의한 밀도류의 유동특성을 전혀 고려하지 못한 채 정수압에 근거한 해수-담수 경계면에 대해 모의하고 있는 정도이다. 따라서 본 연구에서는 해안 대수층 내부의 유동현상을 투수계수에 의존하는 방법에서 탈피하여 대수층 매체의 입경, 공극, 형상 등을 고려할 수 있을 뿐만 아니라, 염분 및 온도차에 의한 밀도류를 해석할 수 있는 강비선형 수치모델을 개발하여 해수침투 현상을 직접 모의한다. 나아가 대부분의 이전 연구들에서 간과하고 있는 해안지역의 대표적 물리력인 파랑과 조석의 영향이 해안 대수층의 해수침투에 미치는 영향, 해안 대수층의 지하수위 및 해수면의 수위차에 의한 해수침투 특성 그리고 이를 제어 할 수 있는 새로운 대응기술을 제안하는 것을 목적으로 한다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.981-985
• /
• 2005

상수도 관망은 상수도 운용시스템에서 가장 중요한 부분이다. 그러므로 항상 수요자의 필요에 부응하는 충분한 수량과 좋은 수질로 공급할 수 있도록 적정한 압력 유지와 수량 확보가 되도록 유지 관리 되어야 한다. 이를 위해서는 상수도 관망의 합리적인 해석과 최적의 설계를 통해 수요자가 필요로 하는 양질의 물을 원활히 공급하고 재원의 낭비를 줄여야 한다. 그러나 우리나라의 많은 대도시들의 경우처럼 특히 상수도 시설이 거대해 지고 급수구역이 광범위하며 복잡한 배수시설을 가지고 있는 지역에 있어서는 모든 현상을 정확하게 파악하고 유지관리하기가 매우 어렵게 되므로 운영관리상 많은 문제점이 발생하게 된다. 이러한 문제점에 보다 효율적으로 대처하기 위해서는 가장 기본적으로 현재의 상수도 시설에 대한 수리적 상태를 파악하여야 한다. 따라서 본 논문에서는 대도시 상수도 시설의 효율적 운영 및 확장에 기본이 되는 현재의 상수도 시설에 대한 시설현황, 운영, 수리상태를 파악 하고, 밸브의 개폐에 따른 수리학적 특성 변화를 제시하고자 하였다. 관망 해석에 필요한 입력 자료는 실제 운영되고 있는 광주광역시의 상수도 급배수 관망 시설을 대상으로 산출 하였고 관망해석 프로그램인 Pipenet 98을 이용하여 시뮬레이션을 실시하였다. 본 논문에서 대상 구역에 대해 수리적 특성 변화를 모의한 결과 상수 관로내의 밸브 개폐도를 $100\%$ 유지시켰을 때 평균 수압은 $4.4kg/cm^2$으로 나타났으며 밸브 개폐도 $80\%$일 때의 평균 수압 $4.7kg/cm^2$, 밸브 개폐도 $60\%$일 때의 평균 수압은 $4.69kg/cm^2$으로 나타났다. 밸브 개폐도가 $100\%$일 때가 밸브를 $60\%$와 $80\%$ 개폐시켰을 때보다 $0.3kg/cm^2,\;0.29kg/cm^2$ 낮게 나타나 밸브를 전체 개방 했을 때 관로내의 수압이 상수설계기준에 적합한 수압을 유지함을 알 수 있다. 상수관로 설계 기준에서는 관로내 수압을 $1.5\~4.0kg/cm^2$으로 나타내고 있는데 $6kg/cm^2$보다 과수압을 나타내는 경우가 $100\%$로 밸브를 개방하였을 때보다 $60\%,\;80\%$ 개방하였을 때가 더 빈번히 발생하고 있으므로 대상지역의 밸브 개폐는 $100\%$ 개방하는 것이 선계기준에 적합한 것으로 나타났다. 밸브 개폐에 따른 수압 변화를 모의한 결과 밸브 개폐도를 적절히 유지하여 필요수량의 확보 및 누수방지대책에 활용할 수 있을 것으로 판단된다.

• Journal of Korea Water Resources Association
• /
• v.46 no.12
• /
• pp.1221-1234
• /
• 2013

The purpose of this study is to suggest the methodology for the computation of uplift pressure and discharge of the seepage flow under gravity dam. A 3-dimensional FDM model is developed for this purpose and this model can simulate the saturated Darcian flow in heterogeneous media. For the verification of the numeric model, test simulation has been executed and the mass balance has been checked. The error does not exceed 3%. Using the developed model, The uplift pressure and seepage flow discharge under gravity dam has been calculated. The uplift pressure shows the similar pattern, comparing with the result of flow-net method. As the length of grout curtain increases, the uplift pressure decreases linearly, but the seepage flow discharge shows the non-linear decreasing pattern. The coefficients of the formulas in the dam-design criteria have been analysed, and ${\alpha}=1/3$ corresponds to the value when the length of curtain grout is 70% of the aquifer height. The uplift pressure near the pressure relief drain has the big curvature vertically and horizontally. The developed model in this study can be used for the evaluation of the effects of seepage flow under gravity dam.