• Clean Technology
• /
• v.25 no.3
• /
• pp.206-222
• /
• 2019

Recently, various multiphase flows have been developed, and among them some models have been commercialized. However, most of them have been developed based on a pressure-based approach; therefore, various numerical difficulties were involved inherently. Accordingly, in order to overcome these numerical difficulties, a multiphase flow model, MultiPhaSe flow (MPS), following a fractional-flow based approach was developed. In this study, by combining a contaminant transport module describing an enhanced dissolution effect of a surfactant with MPS, a MultiPhaSe flow and TranSport (MPSTS) model was developed. The developed model was verified using the analytical solution of Clement. The MPSTS model can simulate the process of surfactant enhanced aquifer remediation including interphase mass transfer and contaminant transport in multiphase flow by using the coupled particle tracking method and Lagrangian-Eulerian method. In this study, a surfactant was used in a non aqueous phase liquid (NAPL) contaminated area, and the effect of hydro-geological heterogeneity in the layered media on remediation efficiency was studied using the developed model. According to the numerical simulation, when hydraulic conductivity in a lower layer is 10 times, 20 times, and 50 times larger than that in an upper layer, the concentration of dissolved diesel in the lower layer is much higher than that in the upper layer because the surfactant moves faster along the lower layer owing to preferential flow; thus, the surfactant enhances dissolution of residual non aqueous phase liquid in the lower layer.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.3
• /
• pp.97-106
• /
• 1989

This paper is on the modeling of two-dimensional groundwater flow, which is the first step of the development of Dynamic System Model for groundwater flow and pollutant transport in subsurface porous media. The particular features of the model are its versatility and flexibility to deal with as many real-world problems as possible. Points as well as distributed sources/sinks are included to represent recharges/pumping and rainfall infiltrations. All sources/sinks can be transient or steady state. Prescribed hydraulic head on the Dirichlet boundaries and fluxes on Neumann or Cauchy boundaries can be time-dependent or constant. Sources/sinks strength over each element and node, hydraulic head at each Dirichlet boundary node and flux at each boundary segment can vary independently of each other. Either completely confined or completely unconfined aquifers, or partially confined and partially unconfined aquifers can be dealt with effectively. Discretization of a compound region with very irregular curved boundaries is made easy by including both quadrilateral and triangular elements in the formulation. Large-field problems can be solved efficiently by including a pointwise iterative solution strategy as an optional alternative to the direct elimination solution methed for the matrix equation approximating the partial differential equation of groundwater flow. The model also includes transient flow through confining leaky aquifers lying above and/or below the aquifer of interest. The model is verified against three simple cases to which analytical solutions are available. The groundwater flow model shall be combined with the model of pollutant transport in subsurface porous media. Then the combined model, with the applications of the Eigenvalue technique and the Dynamic system theory, shall be improved to the Dynamic System Model which can simulate the real groundwater flow and the pollutant transport accurately and effectively for the analyses and predictions.

• 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.

• Journal of Korea Water Resources Association
• /
• v.40 no.6 s.179
• /
• pp.431-446
• /
• 2007

Groundwater flow in a basin is greatly affected by many hydrogeological and hydrological characteristics of the basin. A groundwater flow model for the Kap-cheon basin ($area=648.3km^2$) in the Geum river basin was established using MODFLOW by fully considering major features obtained from observed data of 438 wells and 24 streams. Furthermore, spatial groundwater recharge distribution was estimated employing accurately calibrated watershed model developed using SWAT, a physically semi-distributed hydrological model. Model calibration using observed groundwater head data at 86 observation wells yielded the deterministic coefficient of 0.99 and the water budget discrepancy of 0.57%, indicating that the model well represented the regional groundwater flow in the Kap-cheon basin. Model simulation results showed that groundwater flow in the basin was strongly influenced by such factors as topological features, aquifer characteristics and streams. The streams in mountainous areas were found to alternate gaining and losing steams, while the streams in the vicinity of the mid-stream and down-stream, especially near the junction of Kap-cheon and Yudeong-cheon, areas were mostly appeared as gaining streams. Analysis of water budget showed that streams in mountainous areas except for the mid-stream and up-stream of Yudeong-cheon were mostly fed by groundwater recharge while the streams in the mid and down-stream areas were supplied from groundwater inflows from adjacent sub-basins. Hence, it was concluded that the interactions between surface water-groundwater in the Kap-cheon basin would be strongly inter-connected with not only streams but also groundwater flow system itself.

• Economic and Environmental Geology
• /
• v.38 no.2 s.171
• /
• pp.177-185
• /
• 2005

Data from 122 pumping tests were obtained from 100 boreholes in granites, volcanic rocks, metamorphic rocks, and Cretaceous and Tertiary sedimentary rocks, and then were analyzed using AQTESOLV. Results from 86 of the 122 tests ($71\%$) have an analytical solution corresponding to Theis (1935), Cooper-Jacob (1946), Papadopulos-Cooper (1967), Hantush (1962), Moench (1985), or Hantush-Jacob (1955), whereas the remaining 36 results ($29.5\%$) do not correspond to any of the analytical methods. Of the 86 results, only 17 match the Theis and Cooper-Jacob methods, indicating that the basic methods fer pumping test analysis are useful far only $14\%$ of the total data. This suggests that analytical solutions derived using leaky boundary conditions are appropriate for the analysis of pumping test data in fractured aquifers in this study. Furthermore, the results show the importance of carefully selecting an appropriate model for the analysis of pumping test data. Results from the 122 pumping tests were also analyzed using the GRF model. Using the Barker method, the results show that 77 of the 122 tests ($63\%$) have dimensions ranging between 1.1-2.9. Of these 77 solutions, ($39(44.2{\%})$) have a fractional dimension of 1.1-1.9, ($26(6.5{\%})$) show 2-dimensional radial flow also applicable to the Theis method, and ($38(49.3{\%})$) have dimensions of 2.1-2.9. The results show that groundwater flows according to a fractional flow dimension in fractured aquifers.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2011.05a
• /
• pp.8-9
• /
• 2011

On 4 September 2010 an earthquake of magnitude 7.1 on the Richter scale occurred on the Canterbury Plains in the South Island of New Zealand. The Canterbury Plains are an area of extensive groundwater and spring fed surface water systems. Since the September earthquake there have been several thousand aftershocks (Fig. 1), the largest being a 6.3 magnitude quake which occurred close to the centre of Christchurch on 22February 2011. This second quake caused extensive damage to the city of Christchurch including the deaths of 189 people. Both of these quakes had marked hydrological impacts. Water is a vital natural resource for Canterburywith groundwater being extracted for potable supply and both ground and surface water being used extensively for agricultural and horticultural irrigation.The groundwater is of very high quality so that the city of Christchurch (population approx. 400,000) supplies untreated artesian water to the majority of households and businesses. Both earthquakes caused immediate hydrological effects, the most dramatic of which was the liquefaction of sediments and the release of shallow groundwater containing a fine grey silt-sand material. The liquefaction that occurred fitted within the empirical relationship between distance from epicentre and magnitude of quake described by Montgomery et al. (2003). . It appears that liquefaction resulted in development of discontinuities in confining layers. In some cases these appear to have been maintained by artesian pressure and continuing flow, and the springs are continuing to flow even now. In spring-fed streams there was an increase in flow that lasted for several days and in some cases flows remained high for several months afterwards although this could be linked to a very wet winter prior to the September earthquake. Analysis of the slope of baseflow recession for a spring-fed stream before and after the September earthquake shows no change, indicating no substantial change in the aquifer structure that feeds this stream.A complicating factor for consideration of river flows was that in some places the liquefaction of shallow sediments led to lateral spreading of river banks. The lateral spread lessened the channel cross section so water levels rose although the flow might not have risen accordingly. Groundwater level peaks moved both up and down, depending on the location of wells. Groundwater level changes for the two earthquakes were strongly related to the proximity to the epicentre. The February 2011 earthquake resulted in significantly larger groundwater level changes in eastern Christchurch than occurred in September 2010. In a well of similar distance from both epicentres the two events resulted in a similar sized increase in water level but the slightly slower rate of increase and the markedly slower recession recorded in the February event suggests that the well may have been partially blocked by sediment flowing into the well at depth. The effects of the February earthquake were more localised and in the area to the west of Christchurch it was the earlier earthquake that had greater impact. Many of the recorded responses have been compromised, or complicated, by damage or clogging and further inspections will need to be carried out to allow a more definitive interpretation. Nevertheless, it is reasonable to provisionally conclude that there is no clear evidence of significant change in aquifer pressures or properties. The different response of groundwater to earthquakes across the Canterbury Plains is the subject of a new research project about to start that uses the information to improve groundwater characterisation for the region. Montgomery D.R., Greenberg H.M., Smith D.T. (2003) Stream flow response to the Nisqually earthquake. Earth & Planetary Science Letters 209 19-28.

• Economic and Environmental Geology
• /
• v.40 no.4
• /
• pp.403-418
• /
• 2007

The objective of this research was to characterize the fate and transport of Cr(VI) contaminated groundwater in the Daejeon industrial area. Five subsidiary monitoring wells were newly installed and two existing wells were utilized for the investigation and the reduction process of Cr(VI) contaminated groundwater of the Daejeon(Mun-pyeong) national groundwater monitoring station. The Cr(VI) concentrations at the shallow aquifer well of the station were in the range of 3.2-4.5 mg/L indicating continuous contamination. However, Cr was not detected at the deep bedrock well and the other monitoring wells except MPH-1 and 3. The Cr(VI) concentrations of MPH-1 and MPH-3 were below the drinking water guideline value (0.05 mg/L). Therefore, the plume of the Cr(VI) contaminated groundwater was predicted to be confined within the narrow boundary around the station. The soluble/exchangeable Cr(VI) concentrations were below the detection limit in all core and slime samples taken from the five newly installed wells. Although the exact source of contamination was not directly detected in the study area, the spatial Cr(VI) distribution in groundwater and characteristics of the core samples indicated that the source and the dispersion range were confined within the 100 m area from the monitoring station. The contamination might be induced from the unlined landfill of industrial wastes which was observed during the installation of an subsidiary monitoring well. For the evaluation of the natural attenuation of Cr(VI), available reduction capacities of Cr(VI) with an initial concentration of 5 mg/L were measured in soil and aquifer materials. Dark-gray clay layer samples have high capacities of Cr(VI) reduction ranging from 58 to 64%, which is obviously related to organic carbon contents of the samples. The analysis of reduction capacities implied that the soil and aquifer materials controlled the dispersion of Cr(VI) contamination in this area. However, some possibilities of dispersion by the preferential flow cannot be excluded due to the limited numbers of monitoring wells. We suggest the removal of Cr(VI) contaminated groundwater by periodical pumping, and the continuous groundwater quality monitoring for evaluation of the Cr(VI) dispersion should be followed in the study area.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.8 no.4
• /
• pp.23-32
• /
• 1988

The rates of infiltration contributed to the flow fo water in an unconfined aquifer under the partially penetrated stream at an ungaged station and the corresponding base flow in channel are coupled by using the hydraulic and/or hydrologic characteristics obtained from the geomorphologic and soil maps. For the determination of groundwater flow, the linearized model which is originally Boussinesq's nonlinear equation is applied in this study. Also, a stream flow routing model for base flow in channel is based on a simplification of the Saint-venant. The distributed runoff model with piecewise spatial uniformity is presented for obtaining its solution based on a finite difference technique of the kinematic wave equations. The method developed in this study was tested to the Bocheong watershed(area : $475.5km^2$) of the natural stream basin which is one of tributaries in Geum River basin in Korea. As a result, it is suggested that the rationality of hydro-graph separation according to a wide variability in hydrogeologic properties be worked out as developing the physically based subsurface model. The results of the present model are shown to be possible to simulate a base flow due to an arbitrary rate of infiltration for ungaged basins.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2002.04a
• /
• pp.167-172
• /
• 2002

Seoul subway system has been constructed to solve traffic difficulties of Seoul metropolitan, and now is the major public transportation. However, the more line has added in the system the deeper the bottom of the tunnel base. And a huge amount of groundwater along the line has seeping into the tunnel. Several subway stations has pumping system to extract the groundwater to the outside and consequently, groundwater table along the line has declined gradually. Groundwater table has dropped about 40 meters at some areas, There was some study for the proper usage of the abstracted groundwater and the project to use the groundwater has launched already by the local government. However. more serious problem is expected on quality degradation of soil and groundwater as the decline of groundwater table along the subway line. This study suggests that the detailed groundwater environmental study should be made as soon as possible for this. If there is any pollution leaking at the surface area of the groundwater depression, the pollution will be seep into the subway tunnel in some day even though the time will be different with the soil material and hydraulic characteristics of the aquifer. And the polluted area of the soil and groundwater would be enlarged along the pathway The study on possibility of the soil subsidence and reducing surface water flow in small creek were also needed. This study suggest one of the counter measurement that restoring the declined groundwater table after groundwater environmental study

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2004.04a
• /
• pp.114-118
• /
• 2004

울산지역의 대수층 분포특성과 지하수 유동특성에 관해 연구하기 위해 크리깅 기법을 이용하였다. 울산지역의 시추자료 1,783개 지점과 지하수위자료 총 1,171개 지점에 대하여 2년간(2002년～2003) 자료를 획득하여 분석하였다. 본 연구지역의 표고 충적층 하부경계 심도, 충적층, 하부경계 심도, 충적층 층후, 풍화대 층후, 충적층-풍화대 층후, 대수층 단면, 지하수위 그리고 지하수유동 분석을 하기 위해 베리오그램 분석한 결과 풍화대 하부경계 심도에 지수형모델(exponential model), 나머지 지하수위를 포함한 성분들은 모두 구상형모델(shperical model)이 가장 적합하게 나타났다. 울산지역 대수층의 경우 산악지역은 얕은 충적층과 풍화대 분포가 나타나는 반면, 남구는 충적층, 북구는 풍화대 발달이 우세하게 나타났다. 그리고 울산지역의 충적층과 풍화대 층후 분포특성이 울산단층을 따라 층후가 두껍게 나타났다. 지하수 유동은 고지형의 북구와 동구 산악지역에서 바다에 인접한 내륙지역으로 지하수 유동특성을 보였으며, 울산의 4개구는 산악지역, 태화강, 동천강, 울산단층 그리고 충적층과 풍화대가 잘 발달된 대수층을 따라 지하수유동 특성이 잘 나타났다.