해안지역에서 지하수 순환 특성 규명은 지하수 자원의 효율적 관리 측면에서 매우 중요하다. 본 연구에서는 물수지 분석법 을 이용하여 지하수 침누수량을 산정하였다. 증발산량은 Thornthwaite 방법으로 계산하였으며, 지표 직접유출량은 SCS-CN 방법으로 산정하였다. 지하수 저류량 변화는 229 mm/년로 30년 평균 강수량 1286 mm/년의 17.8%에 해당하며, 증발산량 은 693 mm/년 (53.9%), 지표 직접유출량은 124 mm/년 (9.6%)로 산정되었다. 강우량과 지하수 저류량 변화 사이에는 강우 량과 증발산량, 강우량과 지표 직접유출량에 비해 상관성이 높게 나타난다. 이는 증발산량 및 지표 직접유출량 보다 지하수 저류량 변화가 강우량에 대해서 더 민감하게 반응한다는 것을 반영한다.

본 연구에서는 해수-담수-해안대수층의 비선형 상호작용을 직접 해석할 있는 PBM(Porous Body Model) 기반의 3차원 N-S Solver인 LES-WASS-3D ver 2.0을 적용하며, 해안대수층의 해수침투모의를 수행하였다. 이와 같은 N-S Solver를 적용한 해안대수층의 해수침투모의는 국내 최초 수행되는 것일 뿐만 아니라, 국외적으로도 찾아보기 어려운 새로운 수치해석방법이라고 할 수 있다. 먼저 적용하는 수치모델을 검증하기 위하여 해안대수층의 해수-담수 경계면에 관한 수리모형실험결과와 비교·검토 하여 수치모델의 타당성 및 유효성을 확인하였다. 그리고 해수위 및 지하수위 변화를 고려한 해안대수층 내의 해수침투모의를 수행하여 해수위-지하수위 차와 해수위의 비(△h/h)의 증가에 따른 해안대수층 내의 유동장 그리고 해수-담수 경계면 분포 특성에 관하여 논의하였다. 또한 기존의 비확산 수치모델에서 도출할 수 없었던, △h/h에 따른 해안대수층 내의 연직 염분농도로 부터 해수침투 특성을 파악하였으며, 최종적으로 지표화 할 수 있는 △h/h가 해안대수층 내의 해수침투거리에 미치는 영향에 대하여 분석하였다. 이 결과로부터 △h/h가 작을수록 해안대수층 내의 해수침투가 약해지는 메커니즘을 이해할 수 있었다.

해안지역에서 지하수 해수침투 저감을 위하여 해수 양수 관정을 가동할 때 해수정으로 많은 양의 담수가 낭비되는 사례가 있다. 해수 침투를 저감시키고 낭비되는 담수를 유용하게 활용하기 위하여 이중 양수정이 제안된 바 있다. 이중 양수정은 해수 양수를 위한 해안 쪽의 해수정과 담수 양수를 위한 내륙 쪽의 담수정으로 구성된다. 선행 연구에서는 담수정의 해수 비율을 최소화시키기 위한 해수정의 양수량을 산정하였으나 담수정의 해수비율이 1~15% 범위로 나타나 담수정에서 담수를 확보하기 어렵다고 판단하였다. 본 연구에서는 이송이 지배적인 상황에서 이중 양수정의 효과를 분석하였다. 연구 결과 담수정에서 저염도 지하수를 확보할 수 있음을 밝혔다. 또한 해안 유역 차원의 적정 지하수관리를 위해서는 담수정의 해수비율을 최소화하는 것 외에도 해수정으로 낭비되는 담수량을 최소화하는 것도 함께 고려해야한다.

Due to tidal force, it is very difficult to estimate the hydraulic parameters of high permeable aquifer near coastal area in Jeju Island. Therefore, to eliminate the impact of tidal force from groundwater level and estimate the hydraulic properties, tidal response technique has been mainly studied. In this study we have extracted 38 tidal constituents from groundwater level and harmonic constants including frequency, amplitude, and phase of each constituent using T_TIDE subroutine which is used to estimate oceanic tidal constituents, and then we have estimated hydraulic diffusivity associated with amplitude attenuation factor(that is the ratio of groundwater level amplitude to sea level amplitude for each tidal constituent) and phase lag(that is phase difference between groundwater level and sea level for each constituent). Also using harmonic constants for each constituent, we made the sinusoidal wave and then we constructed the synthesized wave which linearly combined sinusoidal wave. Finally, we could get residuals(net groundwater level) which was excluded most of tidal influences by eliminating synthesized wave from raw groundwater level. As a result of comparing statistics for synthesized level and net groundwater level, we found that the statistics for net groundwater level was more insignificant than those of synthesized wave. Moreover, in case of coastal aquifer which the impact of tidal force is even more than those of other environmental factors such as rainfall and groundwater yield, it is possible to predict groundwater level using synthesized wave and regression analysis of residuals.

해안대수층의 염수침투에 의한 담수-염수 경계면 특성을, 가로 140 cm, 세로 70 cm, 두께 10 cm 의 직사각형 모래수조 모형을 이용하여 모의하였다. 대수층 모형의 특성과 염수의 염도에 따른 총 31가지의 실험을 수행하였으며, 위치흐름 해석법을 유도한 Glover (1959)의 이론식과 호도그래프를 이용한 Henry (1959) 해석식 및 위치흐름 해석법을 이용한 Strack (1976)의 해석해 결과와 비교하였다. 대수층의 수리전도도, 해수의