담수 주입에 따른 담-염수 혼합대에서의 수질 특성 변화

1. 서 언

해안지역은 해수와 담수가 공존하는 지역으로서 밀도가 큰 해수층이 담수 대수층 하부로 침투하여 나타나는 수리적 평형관계를 갖게 된다. 이와 같은 해수 침투는 두 매질의 밀도 차이에 의하여 자연적으로 나타나거나 지하수 양수나 주입 등 인위적 요인에 의하여 조절되는 특성을 갖고 있다. 해수와 담수 경계면은 기본적으로 Ghyben-Herzberg의 이론에 의하여 나타나게 되는데, 지하수의 과다한 양수 등은 이 경계면을 보다 내륙으로 이동시키게 되어 지하수내 염분을 증가시키게 된다.

전세계적으로 해수침투 및 대책에 대한 연구 및 사업이 다양하게 수행되었는데, 미국 서부 캘리포니아 지역의 경우, 1951년에 맨하탄 비치 해안의 방치공을 활용한 주입 시험을 실시한 바 있으며 이를 토대로 1953년부터 9개의 함양정과 54개의 지하수 관측정을 활용한 주입을 통하여 해수침투를 방지할 수 있는 해수방벽 프로젝트(Seawater barrier project)를 시작하였으며, 이후 미국 태평양 해안 및 플로리다 지역 등에 담수 주입을 통한 해수침투 차단 사업이 시행되어 왔다(Lipshie and Larson, 1995). 주입정을 통한 담수의 함양은 해수 침입에 대한 벽체를 형성하는 주요 방법으로서, 피압이나 자유면 또는 층상의 대수층에서 모두 적용 가능한 방법으로 알려져 있다(Mahesha, 2001). 레바논의 Jieh 및 Damour 지역에서는 주입정 함양 방식을 활용한 3열의 차단벽 설치 효과를 모델을 통하여 평가하고 해수침투를 차단하였다(Masciopinto, 2013). 나이지리아의 Lagos 섬 지역에서는 해안지역에서의 지구물리탐사를 통하여 해수 침투를 방지하기 위한 적정 지하수 개발 심도를 평가한 바 있으며 아울러 136 m 깊이의 차수벽 설치를 제안한 바 있다(Akinbinu, 2015). 해안지역 뿐 아니라 내륙으로 해수가 유입될 수 있는 경로는 운하, 배수로, 태풍 등이 해당된다. 운하의 경우에는 배의 이동을 위하여 깊이 굴착하게 됨으로써 해수가 유입되는 경우가 있으며, 멕시코의 Sabine 호에서는 호수와 연결된 대규모의 운하로 인하여 호수내로 해수가 유입되기도 하였다(Barlow, 2008). 아메리카 대륙의 대표 운하인 파나마 운하의 갑문에는 담수를 해수가 대체하는 방식으로 물을 이송시켜 배를 이동시키는 시스템을 가동하고 있으며, 갑문의 운행시 운하 내에는 염수가 유입되게 되는데 이를 방지하기 위하여 갑문내에 공기 순환 셀 시스템에 의한 염분 분리 기술을 도입하고 있다(Jones and Dawson, 1973). 최근에는 해안 대수층에서의 지하수 주입에 의한 해수유입 차단이 대수층의 수질 확보와 해수침투로 인한 농업 활동에 미치는 영향 등을 경제적 관점에서 평가하여 해수침투 대책의 필요성을 강조한 바 있다(Siarkos et al., 2017).

국내의 경우에도 다양한 해수침투 관련 연구가 수행된 바 있다. Lee (1992)는 제주도내 해수와 담수간의 유체동력학적 평형에 대한 연구를 수행하였으며, Hwang et al. (1999)은 임해지역의 지하 비축기지 주변에서의 지하수 유동에 해수침투가 미치는 영향을 2차원 수치모델을 활용하여 평가한 바 있다. Lee et al. (2001)은 부산 동래지역 암반 대수층내에 나타나는 해수침투 특성을 지하수의 이화학 특성 자료를 활용하여 평가하였으며, Park et al. (2002)는 서부 해안 지역에서 채취된 229개의 시료를 활용하여 해안으로부터의 거리와 심도에 따른 광역적인 Seawater mixing index (SMI)를 도출한 바 있다. Park et al. (2008)은 해안 대수층내의 단층의 존재가 해수 침투에 어떤 영향을 미치는지 모델을 통하여 분석하였는데, 단층면에 평행한 방향이 지하수 및 염분 이동에 통로 역할을 할 수 있음을 밝힌 바 있다. Lee and Moon (2008)은 서천-군산 지역의 해수침투 특성을 수질데이터의 시계열 및 요인분석을 통하여 평가하였는데, 이들 지역은 충적/풍화대 대수층에 비하여 암반 대수층을 통한 해수침투가 활발하게 이루어지고 있음을 규명하였다. 또한, 해수 침투의 범위를 파악하는 수단으로서 전기탐사나 전기비저항탐사와 같은 지구물리탐사 방법을 활용하여 담수와 해수의 경계를 결정하였으며, 이화학 분석자료와 전기비저항탐사를 융합한 연구를 수행하기도 하였다(Hagemeyer and Stewart, 1990; Hwang et al., 2000; Choudhury and Saha, 2004; Sherif et al., 2006).

이상과 같은 해수침투 특성에 대한 연구뿐 아니라 해수침투 방지를 위한 연구도 수행되었는데, Kim et al. (2003)는 화성시 해안지역에서의 양수정내 염수가 침투하지 않도록 조석 효율과 Ghyben-Hertzberg 이론을 이용한 적정 양수량을 제시한 바 있으며, Kang et al. (2005)는 연안 대수층에서의 해수침투를 효과적으로 방지하기 위하여 인공 하수처리수의 대수층 함양수로서의 활용성을 Sand box 모델을 활용하여 평가한 바 있다. Jung et al. (2014)는 지하수를 양수하여 해수 침투를 방지하는 대책을 시행하는 경우에, 2중 양수정 방식으로서 해안측의 해수정과 내륙측의 담수정을 복합 설치 운영하는 방식을 설치한다면 내륙측에서 보다 저염소 지하수를 취수함으로써 물 활용성을 제고할 수 있다는 연구를 수행한 바 있다. Lee and Kim (2015)에서는 부산광역시 수영구 지역의 해안 대수층내의 해수침투 저감을 위하여 기존의 방치공을 통한 담수 주입시의 효과를 정량적으로 평가하기 위하여 다차원 수리동역학적 분산 수치모델링을 수행한 바 있다.

본 연구에서는 해수가 포함된 수로 주변에서 염수(이하 수로내의 물이므로 해수 보다는 염수로 기술)유입을 차단하기 위하여, 우물 방식의 지하수 주입을 실시할 경우에 나타나는 주입 전후의 담-염수 혼합층에서의 지하수 수질 특성 변화를 분석하는데 목적을 두었다.

2. 연구지역

본 연구지역은 인천광역시 계양구 상야동 일대로서, GI-Waterway와 굴포천 연결수로의 합류지점에 해당하며, 서측에 구릉성 산지가 존재하나 대부분 평야 지대로서 충적층이 잘 발달된 지역이다(그림 1). GI-Waterway는 한강 하류에서 서해로 연결되는 국내 최초의 인공뱃길로써 주운수로의 길이는 18 km, 폭 80 m, 수심 6.3 m이며, 2011년 9월 말부터 수질관리를 위해 서해안의 해수와 한강 물을 혼합하여 시범 운영하였으며, 2012년 하반기부터 본격 운영을 개시하였다.

Fig. 1.

Location of the study area and injection/observation wells.

인근 지역의 토지는 대부분 벼농사 또는 비닐하우스 재배로 활용되고 있으며, 주 농업용수 공급원은 인공 수로를 통하여 유입되는 한강 또는 굴포천의 물이며 일부 지역은 지하수 관정이 사용되고 있다. 본 지역의 기반암은 선캠브리아기의 호상흑운모편마암이며 상부에 제4기의 충적층이 잘 발달하고 있다. 본 지역의 충적층은 현장 시추자료에 의하면 26~33 m의 두께로 분포하며 상부로부터 일반적으로 표토, 점토, 실트, 모래 및 풍화토의 순으로 구성되어 있다(그림 2). 연구지역에는 GI-Waterway 운영시 설치된 M-19 지하수 관측정이 설치 운영중에 있었으며, 본 연구를 위하여 3개의 관측정(OB-1, OB-2, OB-3)과 1개의 주입정(IW)을 추가로 설치하였다.

Fig. 2.

Vertical distribution of sediments at the injection and observation wells.

3. 대수층의 수리 및 수질 특성

GI-Waterway는 평수위 2.7 El.m 및 2.23 ~ 3.07 El.m의 범위 내에서 운영되고 있으며 갑문의 운영 조건에 따라 7 ~ 8일 주기로 상승 하강을 하고 있다. 연구지역 인근의 지하수위는 GI-Waterway의 물을 채우기 이전에는 수로의 수위 및 주변 관정에서의 양수 영향 등으로 -2.14 ~ 2.50 El.m의 범위를 나타냈으나, 이후에는 1.02 ~ 3.74 El.m 정도로 평균 2.65 El.m를 보이고 있다. 본 지역은 평야지대로서 지형 경사에 의한 지하수 흐름보다는 주로 굴포천과 GI-Waterway의 수위 및 관정의 양수 영향을 받는 지역에 해당한다.

대수층의 수리특성을 파악하기 위하여, 관측정 3개소(OB-1, OB-2, OB-3)에 대하여 순간수위변화시험을 실시하였으며, 평균 수리전도도는 OB-1이 1.44 × 10^-4 cm/sec, OB-2가 1.32 × 10^-4 cm/sec, OB-3이 1.13 × 10^-4 cm/sec로 분석되었다.

시추공내 지하수의 흐름 방향을 파악하기 위하여 주 대수층으로 파악되는 지표하 심도 15.0~16.0 m 구간의 모래층을 대상으로 유향유속검층을 실시한 결과, GI-Waterway에 가까운 OB-1와 OB-2 및 굴포천 연결수로에 가까운 M-19는 유속이 느리면서 SSE 방향으로 주 흐름이 형성되어 있으며, 상대적으로 양 하천으로부터 먼 곳에 위치한 OB-3은 유속이 빠르면서 NNE 방향의 주흐름이 형성되고 있다(표 1). 이와 같은 결과는 굴포천 연결수로에 의하여 영향을 받는 NNE 방향의 빠른 흐름이 연구지역을 지배하고 있으나, GIWaterway와 굴포천 연결수로가 만나는 구간으로서 하천수의 흐름 방향이 충돌하기 때문에 다양한 방향을 보이는 것으로 파악된다.

Table 1.

Result of flow logging at groundwater wells.

지난 6년간의 M-19 관측정의 전기전도도의 변화를 보면, 지표하 11 m (EL. -4 m) 지점은 약 400 μS/cm 정도에서 큰 변화가 없이 유지되고 있으나, 지표하 28 m (EL. -21 m)지점의 경우에는 약 5,000 μS/cm까지 증가하는 특성을 보였으며, 이 기간 동안 GI-Waterway는 약 25,000 μS/cm까지 증가하였다(표 2; 그림 3). 2012년 이후 지표수의 전기전도도가 높아지기 시작한 점을 고려할 때, M-19의 깊은 심도에서만 전기전도도가 증가한 점은 대수층 하부로 염수의 유입이 발생하였음을 보여준다.

Table 2.

Result of hydrogeochemical anaysis for water samples before an injection test.

Fig. 3.

Long-term changes in electric conductivity at M-19 and GI-Waterway.

시간 경과에 따른 연구지역내에 설치된 지하수 관측정 M-19의 지하수의 이화학 특성을 보면, 지표하 11 m 지점에서 채취된 시료는 대부분 Ca(Mg)HCO₃ 유형의 범주에 해당되나 시간이 경과하면서 Ca(Mg)Cl(SO₄) 유형으로 변해가는 것을 볼 수 있다. 지표하 28 m 지점의 시료는 초기를 제외하고 후반기에는 Cl_SO₄ 유형으로 변화되어 GI-Waterway의 시료와 근접하고 있음을 알 수 있다(그림 4).

Fig. 4.

Piper diagrams for water samples at M-19.

한편, 각 관측정 및 주입정내에서의 전기전도도의 수직적인 분포를 파악하기 위하여 CTD diver를 활용한 검층을 실시하였는데, M-19 관측정의 2015년 4월에 실시한 결과에서는 15 m 깊이에서부터 전기전도도 농도가 급격히 증가하기 시작하였음을 보여주며, 2015년 9월에는 약 5,000 μS/cm를 초과하는 것으로 나타났다(그림 5(a)~(d)). 반면에, OB-1의 경우에는 2015년 11월에 심도 15~23 m 구간에서 약 17,000 μS/cm의 높은 값을 보이다가 23 m 이상 깊이에서는 다시 줄어드는 특성을 보이고 있다(그림 5(e)). OB-2의 경우에도 OB-1과 유사한 패턴을 보이는데 약 10,000 μS/cm까지 증가된 전기전도도가 약 28 m 하부에서는 7,000 μS/cm 정도로 감소되는 특성을 보인다(그림 5(f)). OB-3에서는 나머지 2개 관측정에 비하여 전기전도도의 증가폭이 크지 않은 것으로 나타났으며, 주입정 IW에서는 OB-1과 유사한 패턴을 보이는 것으로 나타났다(그림 5(g)). 고 전기전도도의 출현 심도는 OB-3이 OB-1 보다 약간 깊으며 그 값도 작은데 이는 시추주상도에서 보듯이 OB-3의 모래층의 출현심도 및 두께와 관련이 있는 것으로 추정되며, 아울러, 수로에서부터의 거리도 OB-3이 상대적으로 먼 것이 영향을 미쳤을 수 있다. 그러나, 이들 중 어느 것이 주 요인인지에 대해서는 명확하지 않아 보인다.

Fig. 5.

Result of electric conductivity logging at each well before an injection test.

이상에서 보는 바와 같이, 주입 시험 이전까지 연구지역의 지하수는 당초에는 전형적인 지하수-지표수 혼합 수질 특성을 갖고 있었으나 GI-Waterway의 물을 채우면서 염수 유입 영향을 받는 과정이 지난 수년간 지속되어 왔음을 알 수 있다. 이는 본 지역이 GI-Waterway와 굴포천 연결수로가 만나는 돌출 구간이면서 충적층으로 이루어진 지역으로서 GI-Waterway의 타 구역보다 지표수의 대수층내 유입이 용이한 지역이라는 점과 대수층의 수직 분포 및 거리 등이 염수 유입 및 분포에 영향을 미친 것으로 보인다.

4. 주입정의 설계 및 현장 실험

현장 주입 실험은 주입정(IW)을 활용하여 수행되었으며, 주변 관측정(OB-1, OB-2, OB-3 및 M-19)에서 지하수위 및 수질 변화를 관측하였다. 주입정은 굴착경 300 mm, 내경 150 mm, 심도 26 m로서, 지표하 14~19 m, 21~24 m 구간에 유공관(슬롯 크기: 2.0 mm)을 설치하고, 외부는 2.0~5.0 mm의 여재를 부설하였다. 주입수는 생활용 수질기준을 만족하는 인근의 농업용수로내의 지표수를 사용하였으며, 중력에 의한 비압력식 주입을 실시하였다.

주입시험은 예비시험, 단기시험 및 장기시험의 3단계로 진행되었으며, 장기시험 이후에 전기전도도 검층 및 간이수질 분석을 실시하였다(표 3). 장기 시험시 주입량은 초기에 약 74.2 m³/d를 보이다가 5일을 경과하면서 약 34 m³/d의 일정한 값으로 유지되었다(그림 6).

Table 3.

Details on the injection tests.

Fig. 6.

Changes of injection rate for a long-term test.

5. 주입에 따른 수질 변화 특성

각 관측정에서 주입전(6월 27일), 주입 종료 직후(7월 9일) 및 주입 종료 후 7일 경과(7월 16일) 시의 전기전도도 검층을 수행하여 주입 효과의 변화를 파악해 보았다(그림 7). 상부와 하부 구간의 수질 변화 특성이 차이를 보이므로, OB-1의 경우는 15 m, OB-2는 15 m, OB-3은 17 m를 기준으로 상하 구간의 변화를 비교해 보았다.

Fig. 7.

Vertical distribution of electric conductivity at three monitoring wells.

상부 구간의 경우, GI-Waterway에 인접한 OB-1은 5,500 ~ 16,300 μS/cm의 전기전도도를 보이다가 주입 종료 직후에는 3,800 ~ 7,100 μS/cm까지 감소하였으며, 주입 종료 후 7일 경과 후에도 2,200 ~ 8,700 μS/cm 정도로 유지되고 있다. OB-2의 경우에는 800 ~ 9,200 μS/cm의 분포를 보였는데, 주입 종료 직후에는 600 ~ 800 μS/cm까지 하강하였으며, 주입 종료 후 7일 경과 시에도 500 ~ 600 μS/cm 정도로 유지되고 있다. 주입정에 가장 가까이 위치한 OB-3의 경우에는, 2,000 ~ 5,300 μS/cm의 분포를 보였는데, 주입 종료 직후에는 600 ~ 1,800 μS/cm까지 하강하였으며, 주입 종료 후 7일 경과 시에는 800 ~ 1,200 μS/cm로 전체적으로 재상승되는 현상이 나타났다.

하부 구간의 경우, OB-1에서는 15,000 ~ 17,500 μS/cm의 전기전도도를 보이다가 주입 종료 직후에는 4,500 ~ 8,900 μS/cm까지 감소하였으며, 주입 종료 후 7일 경과 시에는 재 상승하여 11,000 ~ 15,600 μS/cm을 보이고 있다. OB-2의 경우에도 하부에 8,500 ~ 10,300 μS/cm의 분포를 보였는데, 주입 종료 직후에는 4,300 ~ 5,200 μS/cm까지 하강하였으며, 주입 종료 후 7일 경과 시에 4,700 ~ 5,700 μS/cm로 약간 상승하는 것으로 나타났다. 주입정에 가장 가까이 위치한 OB-3의 경우에는, 하부에 6,100 ~ 7,400 μS/cm의 분포를 보였는데, 주입 종료 직후에는 2,800 ~ 4,100 μS/cm까지 하강하였으며, 주입 종료 후 7일 경과시에도 2,100 ~ 3,900 μS/cm로 비교적 안정적으로 유지되는 것으로 나타났다.

이와 같이, 주입전에 비하여 주입 종료 직후의 전기전도도의 수직 분포를 보면, OB-1은 약 15 m 깊이, OB-2는 약 15 m 깊이, OB-3은 약 17 m 깊이에서 계단상으로 급격한 변화를 보이는 것으로 나타났다. 이는 주입전에 전기전도도의 수직 분포가 OB-1 및 OB-3에서 완만하게 변하는 점과 차이가 있다. 이와 같은 주입 종료 직후의 급격한 전기전도도의 변화가 발생한 것은 시추주상도에서 보듯이 이들 심도에서 주 대수층인 모래층이 발달하기 때문으로서 모래층을 통한 확산이 용이했음을 의미한다. 주입 종료 후 7일 경과 시에는 GI-Waterway에 가장 가깝고 주입정에서 가장 먼 OB-1공에서 염수의 재유입이 빠르게 나타났으며, OB-2 및 OB-3에서는 주입 효과가 상대적으로 지속되었음을 보여준다. 또한, 동일 모래층이라고 하더라도 부분적인 투수성의 차이 및 경계부 등에서 주입수의 확산 및 염수의 재유입 정도에 차이가 나타나기도 한다.

상부 심도와 하부 심도로 구분하여 구간별 전기전도도의 평균값을 토대로 주입 전후의 값을 비교하여 주입 효과를 평가해 보았다(표 4). 상부 심도의 경우에는 OB-3에서 75.0%, OB-2에서 58.8%, OB-1에서 55.8%, M-19에서 16.7%의 순으로 효과가 나타났는데, 주입정으로부터의 거리에 따라 주입 효과가 반비례하는 것으로 평가되었다. 반면에, 하부 심도의 경우에는 OB-1에서 63.7%, OB-3에서 52.9%, OB-2에서 46.1%, M-19에서 30.3%의 순으로 효과가 나타났으며, 이는 주입정과의 거리와 무관하게 나타났다. 이와 같은 특징은 하부 모래층(그림 2의 sand 층)의 지하수 흐름이 상부보다 이방성이 크기 때문으로 보인다.

Table 4.

Effect of freshwater injection according to electric conductivity logging at four monitoring wells.

담수와 염수의 혼합대의 분포를 파악하기 위하여 쌍극자 배열에 의한 전기비저항탐사를 2측선씩 주입전(6월 27일), 주입중(7월 2일) 및 주입 종료 후(7월 9일)로 구분하여 수행하였다(그림 1, 8). 측선의 전극 간격은 5 m, 조사 심도는 25 m로 설정하였다.

Fig. 8.

Results of electric resistivity survey along L-1 and L-2.

측선 L-1의 경우, 비저항의 분포범위는 주입전에 13.5 ~ 88.1 Ω․m, 주입중에 3.3 ~ 105.0 Ω․m, 주입후는 16.6 ~ 61.9 Ω․m로 나타났다. OB-1 주변의 비저항의 변화 양상을 보면, 주입전에는 심도 5.0 ~ 17.0 m 정도에서 중간 정도의 비저항대가 존재하나, 주입중에는 측점 No. 10과 No. 11 사이에서 담수 주입 효과가 잘 나타나는 반면에 OB-1 주변에서는 하부로부터 저비저항대의 영역이 넓어지면서 저비저항대로 변화되고 있다. 이와 같이 초기 주입 4일 후인 7월 2일의 저비저항대가 넓어지는 것은 담수 주입시 예상되는 현상과 차이를 보이는데, 이는 주입으로 인한 지하수 흐름 변화로 하부 구간에서 염수화된 물의 이동이 발생하였거나 탐사 결과의 공간적 오차로 인한 것으로 추정된다.

측선 L-2의 경우에는, 비저항의 분포범위는 주입전에 9.39 ~ 120.0 Ω․m, 주입중에 13.3 ~ 83.7 Ω․m, 주입 종료 후는 5.7 ~ 313.0 Ω․m로 나타났다. OB-3 및 IW 인근의 변화 양상을 보면, 주입전에는 심도 7 m 부근에서 중간 내지 높은 비저항대를 보이고 심도 10.0 m 하부에서는 낮은 비저항대가 존재하나, 주입중에는 OB-3의 10 m 하부에서 약간의 상승이 나타나고 있어 주입 효과가 반영된 것으로 보인다. 또한, 주입 종료 후에는 높은 비저항이 유지되어 주입 효과가 지속되고 있음을 알 수 있다.

이상 두 측선에 대한 비저항 분포의 주입 시간에 따른 변화를 보면, 주입 초기에는 담수의 물리적 이동이 우세하기 보다는 주입으로 인한 기존 염수 성분의 이동이 발생하는 것으로 보이며, 측선 L-2에서 보듯이 주입 후 지체 시간을 보이면서 담수 효과가 발생하는 것으로 분석된다. 다만, 주입에 따른 수질 변화와 전기비저항탐사 결과가 명확히 일치하지 않는 것은 탐사의 한계점으로 보인다.

주입 전(6월 21일), 중(7월 1일), 종료 후(7월 9일 및 16일)의 지하수 시료에 대한 이화학 분석 결과를 Piper diagram에 도시해 보았다(그림 9). 상부 구간시료는 지표하 10 m, 하부 구간 시료는 지표하 22 m에서 채취하였다. 주입 전(6월 21일)의 수질을 보면, 상부 구간의 경우에는, OB-1과 OB-3에서 염수의 혼합을 보이는 영역에 도시되었는데 이는 그림 7에서도 유사하게 나타나고 있다. 반면에 OB-2와 M-19는 지하수의 영역에 도시되어 있음을 알 수 있다. 하부 구간의 경우에는 모든 관측정이 염수의 혼합 특성을 보이고 있으며 전기전도도 검층 결과에서도 알 수 있다.

Fig. 9.

Piper diagram for water quality before, during, and after an injection.

주입 중 및 후의 수질 변화를 보면, 상부 구간의 경우 대부분의 관측정(OB-1, OB-2, M-19)에서 주입 중인 7월 1일의 수질 특성이 주입 전보다 다이아그램의 우측에 위치하여 염수 성분이 확대된 것처럼 보이며, 주입 종료 후에는 주입의 효과로 인하여 좌측으로 다시 이동하여 주입 효과가 발생함을 알 수 있다. 하부 구간의 경우에는, OB-1은 주입 중(7월 1일)에는 담수화 현상이 보이지 않으나 주입 종료 직후인 7월 9일에는 주입 효과가 잘 나타났으며, 7월 16일에는 다시 염수화 됨을 알 수 있다. 주입수의 Na⁺ 농도는 32.57 mg/L, Cl^- 농도는 49.2 mg/L, Ca²⁺는 22.09 mg/L 등으로서 염수화된 지하수의 농도에 비하여 매우 낮아 주입으로 인한 희석 또는 배수에 의하여 수질 변화가 발생된 것으로 보인다. OB-2, OB-3 및 M-19의 경우에는 주입에 따른 수질 변화 효과가 잘 나타나지 않는데, 이는 그림 7에서 보듯이 주입에 따른 수질변화의 폭이 상대적으로 적기 때문으로 보인다. 즉, 상부 구간이 하부 구간에 비하여 주입에 따른 이화학적 특성 변화 및 효과가 뚜렷한데, 이는 모래층으로 구성된 하부 구간(그림 2의 sand 층) 특성이 반영된 것이며 아울러 하부가 상부 구간에 비하여 염수 침입이 우세하였기 때문으로 보인다.

이상 주입 전, 중 및 후의 전기전도도 검층 및 이화학 분석 결과에서 보면, 주입 초기에는 인근 관측정까지 담수의 확산이 이루어지지 못하나 일정 시간이 지나면서 주입 효과가 반영되는 것을 알 수 있으며, 중단 시에도 일정 기간 동안 주입 영향이 지속됨을 알 수 있다.

6. 결 언

본 연구는 GI-Waterway와 굴포천 연결수로가 만나는 지점에서 해수가 혼합된 염수의 주변 대수층으로의 유입 특성을 파악하고, 관정을 통한 담수의 주입시 나타나는 지하수 수질 변화 및 염수 차단 효과를 평가하였으며 다음과 같은 결론을 도출하였다.

첫째, 본 지역은 GI-Waterway와 굴포천 연결수로가 만나는 꼭지점 형태의 충적층 지역으로서 지표수의 영향을 쉽게 받는 지역이며, 해수가 혼합된 지표수의 영향으로 담수 이후 대수층내 지하수 수질 변화가 지속되었던 지역으로서 지하수의 이화학특성은 지하수~지표수 성분에서 담-염수 혼합대 지하수 성분으로 변화하고 있다.

둘째, 담수 주입시 각 관측정에서의 주입 효과는 투수성이 양호한 하부의 모래층(그림 2의 sand 층)을 통하여 잘 나타나고 있다.

셋째, 담수 주입 전후의 전기비저항탐사를 활용한 염수 분포 범위를 평가한 결과, 비저항값은 전기전도도와 반비례 관계를 보이므로 추가적인 연구를 통하여 염수 침입 및 차단 효과에 대한 확인 수단으로 활용 가능할 것이다.

넷째, 지하수 관측정에서의 전기전도도 검층과 이화학 특성 변화를 보면, 주입 효과는 일정 시간이 경과한 이후에 나타나기 시작하며 주입 중단 이후에도 주입 효과가 지속되는 것을 알 수 있었다. 이와 같은 특성을 고려할 때, 인공함양 방식의 주입에 의한 염수 침입 방지의 설계를 위해서는 장기간의 주입 시험과 주입량 및 간격 등을 고려한 다양한 형태의 현장 시험을 통하여 주입 주기나 기간을 설정할 필요가 있다.

References

• Akinbinu, V.A., (2015), Delineation of saline water intrusion to safe-guide inland groundwater resources, Ocean & Coastal Management, 116, p162-168. [https://doi.org/10.1016/j.ocecoaman.2015.07.005]
• Barlow, P.M., (2008), Preliminary Investigation: Saltwater Barrier-Lower Sabine River, Sabine River Authority of Texas, Retrieved 2012-12-09.
• Choudhury, K., and Saha, D.K., (2004), Integrated geophysical and chemical study of saline water intrusion, Ground Water, 42(5), p671-677. [https://doi.org/10.1111/j.1745-6584.2004.tb02721.x]
• Hagemeyer, R.T., and Stewart, M., (1990), Resistivity investigation of saltwater intrusion near a major sea-level canal, in S. H. Ward ed., Geotechnical and Environmental Geophysics, vol. II: SEG, p67-77.
• Hwang, H.S., Lee, S.K., Ko, D.C., Kim, Y.S., and Park, I.J., (2000), Detection of sea-water intrusion caused by tidal action using DC resistivity monitoring, Geophysics and Geophyiscal Exploration, 3(1), p1-6, (in Korean with English abstract).
• Hwang, Y.S., Bae, H.S., Seo, D.I., Kim, K.S., and Kim, C.S., (1999), Groundwater flow characteristics affected by the seawater intrusion near simulated underground storage caverns in the coastal area, Journal of Engineering Geology, 9(1), p17-29, (in Korean with English abstract).
• Jones, M.L., and Dawson, C.E., (1973), Salinity-temperature profiles in the Panama canal locks, Marine Biology, 21, p86-90. [https://doi.org/10.1007/BF00354602]
• Jung, E.T., Lee, S.J., Lee, M.J., and Park, N.S., (2014), Effectiveness of double negative barriers for mitigation of seawater intrusion in coastal aquifer: Sharp-interface modeling investigation, Journal of Korea Water Resources Association, 47(11), p1087-1094. [https://doi.org/10.3741/JKWRA.2014.47.11.1087]
• Kang, J.K., Lee, S.J., and Kim, C.G., (2005), Investigation on hydraulic properties according to artificial recharge and extraction, Journal of Korean Society of Environmental Engineers, 27(9), p995-1005, (in Korean with English abstract).
• Kim, Y.Y., Lee, K.C., and Jung, H.J., (2003), Hydrogeological properties of coastal aquifer in Hwaseong area: implication for tidal effects, Journal of the Geological Society of Korea, 39(3), p385-390, (in Korean with English abstract).
• Lee, B.J., Cho, B.W., Moon, S.H., and Lim, M.T., (2001), Seawater enchroachment into bedrock aquifer in dongnae, Pusan, Journal of the Geological Society of Korea, 37(3), p407-420, (in Korean with English abstract).
• Lee, B.J., and Moon, S.H., (2008), Integrated approach for evaluating the characteristics of seawater intrusion using factor analysis and time series analysis: Seocheon-Gunsan area, Journal of the Geological Society of Korea, 44(2), p219-232, (in Korean with English abstract).
• Lee, J.M., and Kim, J.M., (2015), Three-dimensional numerical modeling of impacts of fresh water injection on groundwater flow and salt transport in a coastal aquifer due to urbanization of Suyeong-gu, Busan, Korea, Journal of the Geological Society of Korea, 51(2), p203-219, (in Korean with English abstract). [https://doi.org/10.14770/jgsk.2015.51.2.203]
• Lee, K., (1992), On the hydrodynamic balance of fresh and salt waters in the Cheju Island, Journal of the Geological Society of Korea, 28(6), p649-664, (in Korean with English abstract).
• Lipshie, S.R., and Larson, R.A., (1995), The west coast basin, Dominguez Gap, and Aalamitos seawater-intrusion barrier system, Los Angeles and Orange counties, California, AEG News, 38(4), p25-2.
• Mahesha, A., (2001), An overview of control of salt water intrusion in coastal aquifers, ISH Journal of Hydraulic Engineering, 7(1), p58-71. [https://doi.org/10.1080/09715010.2001.10514690]
• Masciopinto, C., (2013), Management of aquifer recharge in Lebanon by removing seawater intrusion from coastal aquifers, Journal of Environmental Management, 130, p306-312. [https://doi.org/10.1016/j.jenvman.2013.08.021]
• Park, J.H., Kihm, J.H., Kim, H.T., and Kim, J.M., (2008), Three-dimensional numerical simulation of impacts of fault existence on groundwater flow and salt transport in a coastal aquifer, Buan, Korea, Journal of Soil and Groundwater Environment, 13(5), p33-46, (in Korean with English abstract).
• Park, S.C., Yun, S.T., Chae, G.T., and Lee, S.K., (2002), Hydrogeochemistry of shallow groundwaters in western coastal area of Korea: A study on seawater mixing in coastal aquifer, Journal of Soil and Groundwater Environment, 7(1), p63-77.
• Siarkos, I., Latinopoulos, D., Mallios, Z., and Latinopoulos, P., (2017), A methodological framework to assess the environmental and economic effects of injection barriers against seawater intrusion, Journal of Environmental Management, 193, p532-540. [https://doi.org/10.1016/j.jenvman.2017.02.051]
• Sherif, M., Mahmoudi, A.E., Garamoon, H., Kacimov, A., Akram, S., Ebraheem, A., and Shetty, A., (2006), Geoelectrical and hydrogeochemical studies for delineating seawater intrusion in the outlet of Wadi Ham, UAE, Environmental Geology, 49(4), p536-551. [https://doi.org/10.1007/s00254-005-0081-4]