서 론

급격한 도시화와 산업화로 인해 환경에 대한 인식이 낮아져 하천은 점차 치수와 이수의 목적으로 획일화된 형태로 변형되었다. 특히 도시하천은 도심의 하수도 역할을 하는 단순한 배수로로 전략하게 되어 비점오염원에서의 오염물질 유입에 따른 수질오염 등으로 인해 생태기능 상실 및 생활환경에 직접적으로 영향을 미치게 되었다. 이에 따라 하천환경 개선의 필요성이 인식되어 각 지방자치단체는 매년 생태복원하천과 같은 사업을 진행하여 콘크리트를 이용한 인위적인 직선형 하천에서 각 지형에 적합한 형태 즉, 자연적인 굴곡형 물길, 수초대 및 정화 습지 등과 같은 자연정화 기능을 갖춘 하천복원사업을 활발히 진행 중이다.

하천아래 하상퇴적층의 적당한 깊이에 수평으로 한 개 혹은 여러 개의 집수관을 설치하고 이들을 집수정에 연결하여 이로부터 여과수를 확보하는 하상여과공법은 하천복원사업 중 수질개선, 취수량 증대 및 하천의 건천화 방지 효과를 얻기 위해 국내 지방자치단체(대전 갑천, 서울 홍제천 및 울산 태화강 등)에서 이용되고 있으며, 이를 통해 긍정적인 효과를 얻고 있다(Kim, 2007). 하상여과시설 시공을 위한 설계는 하천현황, 지반조사 등을 통해서 기본적인 현장자료를 획득하게 되고, 이를 바탕으로 수리학적 수치해석모델링을 통하여 목표취수량을 산정하게 된다. 하상여과시설의 취수량은 집수정의 개수, 길이, 방향 및 수평 집수관의 위치 등과 같은 설계인자에 따른 영향과 미립토사의 유입, 유기물의 축적에 따른 투수성 저하 및 주변지하수의 수위 등과 같은 시공지역의 환경에 따른 영향을 받게 된다(Sung et al., 2003). 이와 같은 하상여과의 설계인자를 파악하기 위해 Sohn et al. (2011)은 모델링 및 수치해석을 이용하여 하상여과용 인공하천 설계연구를 통해 인공하천의 길이가 길어질수록 폐색방지에 유리하나, 빠른 하천유속을 발생시켜 하상의 폐색을 방지하는 방법은 효율이 떨어진다고 보고하였고, Ahn et al. (2005)은 모형 집수관 실험을 통해 집수관의 길이와 직경이 증가하면 산출유량은 증가하지만, 한계생산성은 감소한다고 보고하였다. 이외에도 하천과 지하수를 대상으로 수질특성 및 지구화학적 해석을 이용하여 주요 성분의 지화학적인 거동(Kim and Nam, 2009), 오염원의 기원해석(Yun et al., 2003), 오염원 및 유해성 평가(Na and Son, 2005) 등과 같은 연구가 다양하게 진행 중이나, 하상여과 시설에 미치는 설계인자 해석은 부족한 실정이다. 이에 따라서, 본 연구는 현재 하상여과공법을 이용하여 시공 중인 충남 아산시 곡교천을 연구대상지역으로 선정하여 하천과 주변지하수의 이화학원소 함량과 안정동위원소 함량비 등과 같은 현장자료를 통해서 지구화학적 해석을 이용한 하상여과 설계인자에 미치는 영향인자를 파악하고, 기초자료 및 가능성을 제시하고자 한다.

연구방법

시료채취 및 분석방법

연구대상지역에 대하여 하상여과에 미치는 설계인자를 파악하고자 하천수와 주변지하수에 대하여 시료를 채취하여 분석을 실시하였다. 하천수는 하천흐름방향에 따라 상류, 중류, 하류에서 각각 채취하였으며, 주변 지하수는 하천으로부터 거리에 따른 영향을 고려하여 9개의 지하수 관정에서 멸균시료용기를 이용하여 2 L 이상 채수하였다(Fig. 1). 주요 양이온, 음이온 및 안정동위원소시료에 대하여 0.45 μm 여과지를 이용하여 부유물질을 제거하였다. 양이온시료는 시료용기 벽면에 양이온이 흡착되는 것을 방지하기 위해 농질산을 첨가하여 pH를 2 이하로 조절하였으며, 시료채수 이후 분석까지 모든 운반 및 저장은 아이스박스와 냉장고를 이용하여 4℃ 이하로 유지하였다.

Fig. 1.River and groundwater sampling sites.

온도, pH, Eh, DO, EC, TDS와 같은 물리적 특성은 현장에서 채수 즉시 다항목 수질측정기(U-5000, HORIBA)를 사용하여 현장에서 측정하였고, 중탄산 이온함량은 적정법을 이용하여 현장에서 pH meter(D-51, HORIBA)를 이용하여 직접 측정하였다. 양이온(Ca, Fe, K, Mg, Na, Si 등)의 함량은 유도결합플라즈마(ICP1000-III, SHIMADZ), 음이온(F, Cl, NO2, Br, NO3, PO4, SO4)의 함량은 이온크로마토그래피(883 Basic IC Plus, metrohm)를 이용하여 분석하였다. 안정동위원소인 질소안정동위원소(δ15N)는 가스-질량분석기(FIAstar 5000, Foss Tecator)를 이용하여 분석하였다. 채취된 하천수 및 주변 지하수 시료내 미생물의 존재 유·무에 따른 탈질산화반응의 영향을 조사하기 위해 이들 시료는 0.2 μm 여과지를 이용하여, 글루타르알데히드(glutaraldehyde)와 에틸 알콜(ethyl alcohol, 75, 85, 95 및 100%)을 각각 이용하여 고정 및 단계별 탈수를 실시한 후, SEM(Hitachi, S4800)을 이용하여 관찰하였다.

결과 및 토의

연구지역의 지질 및 지역적 현황

연구지역인 충남 아산시 곡교천의 지질은 선캠브리아기의 호상편마암이 기반암으로 분포하고 있고, 중생대의 쥬라기 화강암류인 각섬석화강암, 흑운모화강암, 반상화강암과 백악기의 흑운모화강암이 관입하였으며, 이들 모든 지층을 관입한 백악기의 산성, 염기성 맥암류 및 제4기 충적층으로 구성되어 나타난다(Kim et al., 1979). Kim et al. (1979)에 의하면 호상편마암은 아산시 전역에 걸쳐 넓은 분포를 보이며, 크게 북측과 남측으로 나누어 암상의 차이를 보인다. 북측은 온양편마암인 호상편마암이 광역변성작용으로 인해 전반적으로 화강암질화되어 산출되며, 남측은 호상편마암이 넓은 분포를 보이고 있다.

본 연구지역은 곡교천 하류지역을 중심으로 주변 논, 밭, 과수원 등의 농업활동이 활발하며, 축산농가 등의 비점오염원이 산재되어 있다. 또한 대상 지역의 남쪽으로는 아산시청이 위치해 있으며, 이로 인한 아파트 등의 주거 밀집지역이 존재한다. 대상하천 남쪽으로 약 500 m 거리 내에 대형 제지 공장과 수 개의 주유소 등이 위치하고 있다(Fig. 1). Chon and Ahn (1998)은 대상 연구지역에 대하여 GIS를 이용한 지하수 오염 평가에 관한 연구를 수행한 결과, 대상 지역에 대하여 기존 공장에서 사용되고 있는 지하수는 주변 민가에서 사용하고 있는 개발심도보다 깊어, 지하수오염이 진행될 경우 주민들의 건강에 심각한 문제를 야기할 수 있음을 제시하였다.

일반적으로 지하수는 주변의 암석과 계속 반응하므로 대수층의 지질에 따라 수리지구화학적 특성에 차이가 발생된다. 이러한 현상은 대수층의 지하수뿐만 아니라 지표수도 기반암 종류에 따라 지구화학적 특성이 다르게 나타난다. 연구지역에서 시료채취가 이루어진 지하수공에 대한 지질을 지질도를 통해 확인한 결과 하천에 의해 형성된 충적층으로 나타났다. 국가지하수정보센터에서 아산시 지하수 정보신청을 통해 심도 자료를 확인한 결과 시료를 채취한 관정의 심도는 10~100 m 범위로 나타났다. 이들 지역은 분급이 양호하여 투수성이 높아 지하수 산출성이 좋으나, 미고결 퇴적물로 이루어진 충적층에서의 지하수 수위는 지표면에 가깝게 상승되어 있어 지표수가 지하로 침투되어 자체정화 시간이 짧아 지하수의 수질은 지표수 오염원의 영향을 받기 쉽고, 하상여과 설계인자에 영향을 미칠 것으로 판단된다.

연구지역의 지화학적 특성

하천수와 하천 인근지역의 관정에서 채취된 지하수 시료에 대한 물리적 특성은 Table 1에 나타냈다. 수온의 경우 하천수 17.96~18.42℃, 지하수 14.42~21℃ 범위를 나타내고 있다. pH는 하천수 7.11~11.53, 지하수 6.54~8.43 범위로 나타난다. 하천 중류(OR-2)의 경우 pH가 11.53으로 높게 나타난다. 하천수의 높은 pH는 오염물질의 유입에 의해서 나타날 수 있으나, 대부분의 지표수는 대기와 충분히 접하고 있으므로 CO2의 용해에 의한 탄산의 형성으로 실제 그와 같이 높은 pH가 나타나기 어렵지만(Kim and Nam, 2009), NO3함량 및 영양염류가 과다하게 녹아있을 경우 하천에 존재하는 조류의 광합성 작용에 의해 CO2가 소비되어 pH가 증가될 수 있다(Kim et al., 2002). 본 연구지역 하천 중류의 pH는 암상의 변화보다는 시료 채취 시 녹조 등의 조류가 하천 표면에 형성된 것으로 보아 생물학적 요인이 더 크게 작용한 것으로 판단된다. 지하수 pH는 평균 7.41로 중성 및 약산성의 특징을 보인다. 하천수와 지하수의 Eh는 각각 평균 160.00, 175.22 mV로 전이환경으로 나타난다. EC는 하천수 40.5~81.8 mS/m (평균: 56.5 mS/cm), 지하수 0.1~78.7mS/m (평균: 43.5 mS/m) 범위로 나타난다. Lee and Woo (2002)에 의하면 오염되지 않은 일반수(원수)의 전기전도도는 5~50 mS/m 범위로 나타나며 국지적인 매질의 차이, 매립물질의 차이 및 지표환경의 차이에 의해 수질변화가 일어난다고 보고하였다. 본 연구지역 중 하천 하류(OR-1)는 하상여과시공으로 인해 곡교천 하류지역을 차단하여 국지적으로 물이 순환되지 못하여 하천 상류, 중류보다 높게 나타나며, 지하수의 경우 주유소, 제지공장 등과 같은 지표환경의 영향으로 인해 전기전도도가 높게 나타난 것으로 판단된다.

Table 1.Summary of the physical properties of the study area.

하천수와 지하수시료에 대한 양·음이온의 주성분 원소들의 함량을 Table 2, 3에 나타냈다. 이 중, Freeze and Cherry (1979)는 NO3를 생활하수, 축산폐수 등과 같은 오염물질의 유입으로 지하수를 부화시키는 전형적인 오염물질로 규정하였다. 이러한 NO3를 분석한 결과, 하천수와 지하수함량은 각각 9.03~15.59 mg/L, 0.38~173.95 mg/L로 나타났다. 하천수의 경우 주변의 특정 오염원의 존재보다는 하상여과 시공으로 인하여 하천 바닥에 축적되었던 유기물들이 하천으로 유입된 결과로 사료된다. 지하수는 큰 분산범위를 보이며, 하천주변(OG-1, OG-2, OG-3, OG-4)과 도심지역(OG-5, OG-6, OG-7, OG-8, OG-9)에 따라서 NO3의 함량이 크게 차이가 나타나는데 도심지역 지하수의 경우 주변의 산업시설, 농업활동 및 거주시설에 따른 영향으로 판단되며, 특히, OG-7의 경우 농업활동이 진행 중인 과수원 하부 지하수공으로 비료사용으로 인해 높게 나타난 것으로 해석된다. SO4분석결과, 하천수와 지하수는 각각 76.97~126.05mg/L (평균: 108.97mg/L), 5.37~123.43mg/L (평균: 54.33 mg/L) 범위로 높은 함량을 보인다. 일반적으로 SO4 함량이 높은 지역은 광산지역의 황철석과 같은 황화광물의 산화작용(Park and Cho, 2010), 증발암 용해, 강수, 섬유산업과 펄프제조 공정에서의 공장폐수 및 비료에서 공급 된다(Cortecci et al., 2002). 연구지역 지하수의 SO4 함량은 하천인근 지하수(평균: 75.25 mg/L)와 도심지 지하수(평균: 37.58 mg/L)간에 차이가 나타난다. 하천인근 지하수의 경우 물-암석반응보다는 하천수의 혼합, 인간 활동에 의한 부화 및 지구화학적인 변화보다는 하상여과공사에 따른 오염물질의 유입과 하상퇴적물의 용존 이온흡착, 포화에 의한 용출작용(Kim and Nam, 2009)으로 생각된다. 하천수와 주변지하수의 물리·화학적 특성을 이용한 하상여과에 미치는 설계인자(산화-환원 전위의 차이, EC의 변화 및 상대적 변화량 등)를 파악하기 위해서는 상대적 거리에 따른 시료채취와 강수기록 등 집중적이고 체계적인 자료가 요구될 것으로 판단된다.

Table 2.Summary of water chemistry (cations) in the study area.

Table 3.N.D: not detected

수리 화학적 특성

연구지역의 하천수와 주변 지하수에 대한 수질특성을 파악하기 위해 분석된 수질자료를 이용하여 Piper diagram에 도시한 결과, 양이온은 Ca 우세형, 음이온은 Cl 우세형을 보인다. 양이온은 밀집형태를 음이온은 분산형태로 나타난다. 수질의 유형비는 Ca-Na-HCO3, Ca-Cl, Ca-HCO3-Cl 유형을 보이고 있다(Fig. 2). 지하수의 Ca, Cl, HCO3 및 Na 농도분포를 비교한 결과, Cl을 제외한 나머지 원소는 크게 차이가 나타나지 않는다. Cl의 경우, 하천 인근지하수는 도심지역(평균, Cl < 63 mg/L)보다 비교적 높게 분포되며, 하천 인근지하수는 120 mg/L 이상인 이상분포대가 형성된다.

Fig. 2.Piper diagram showing the distribution of major cations and anions in the study area.

배경수질의 차이를 구분하는 주된 원소를 확인하기위해 Scholler diagram을 이용하여 주성분 이온함량을 도시해 본 결과(Fig. 3), 하천수와 지하수의 차이는 K, HCO3, 지하수(하천 인근, 도심지역)에서는 NO3가 구별되는 것으로 나타난다. 하천수의 HCO3는 지하수에 비해 암석, 토양 및 유기물기원과의 반응에 따른 결과로 생각되며, K는 하천수의 경우 지표의 오염물질, 지하수는 자연적인 지질매체와의 반응으로 판단된다.

Fig. 3.Schoeller diagram showing the distribution of major cations and anions in the study area, for: (a) river water, (b) Riverside groundwater, and (c) urban groundwater.

Yun et al. (2003)은 Na와 Cl의 상대적인 농도는 이들 성분의 기원을 규명하는데 유용한 매개변수로 이용된다고 보고하였다. 대수층의 Na-장석류(NaAlSi3O8)가 풍화작용에 의해 Na가 공급되는 경우, Cl/Na의 몰농도비는 약 0.5로 감소하게 되며, 해무기원인 연무질 암염이 용해되는 경우는 Cl/Na의 몰농도 비가 1에 근접한 값을 나타낸다. 연구 지역 하천수와 지하수의 Na와 Cl 몰농도 비를 도시한 결과(Fig. 4) Cl의 성분이 Na에 비해 높은 몰농도 비를 나타내며, 도심 지하수의 경우 해수의 혼합에 따른 영향보다는 주변 산업시설 및 거주시설에 따른 인위적 오염영향으로 판단되며, 하천수와 하천 인근지하수는 염분이 해무기원인 연무질 암염의 용해와 수암 반응에 의해 복합적으로 공급된 것으로 사료된다.

Fig. 4.Relative micromol ratios of sodium and chloride in the river water, groundwater, and urban groundwater samples. The fitted lines show the expected relative concentrations of species assuming aerosol or host rock.

대수층의 수문구조에 따라 수질특성에 차이가 나타나며 이를 기반으로 하천수의 지하유입 혹은 지하수의 하천유입을 판단한다(Simpson and Herczeg, 1991; Lamontagne et al., 2002). 본 연구지역의 수리 화학적 특성 중 K, HCO3, NO3 및 Cl 원소가 하상여과 설계 인자에 영향을 줄 것으로 판단된다.

질소안정동위원소

하천수와 지하수에 대하여 δ15N (‰)를 분석한 결과 +2.94 ~ +14.87‰ 범위를 보이며, 하천수 평균 11.88‰, 지하수 평균 9.06‰로 나타난다(Table 4). 지하수는 하천 주변지하수 +3.95 ~ +12.35‰ (평균: 8.04‰), 도심 지하수 +2.94 ~ +14.87‰ (평균: 9.30‰) 범위를 보인다. 질소안정동위원소는 자연 존재비(δ15N) 값을 이용한 유래추정, δ15N값을 이용한 질산성질소의 기원별 구성비의 추정 등과 같은 방법을 이용하여 질산염의 기원을 추적하는데 많이 이용된다(David, 2001; Mayer et al., 2002) 또한, 지하수내 δ15N값은 지하수가 이동하는 과정에서 암모니아 휘발반응, 질산화작용, 탈질산화작용, 이온교환반응 및 식물신진대사 등의 여러 가지 반응에 의해 변화하기 때문에 질소동위원소를 이용하여 정량적인 오염원을 파악하는데 이용된다(Hauck et al., 1972).

Table 4.Nitrogen isotope data.

연구지역의 NO3의 농도에 대한 δ15N값을 도시하여 보았을 때 NO3의 농도, 시료채취거리, 주변 지류의 영향 및 δ15N값의 차이에 따라 A, B, C그룹으로 구분된다(Fig. 5). B그룹은 특징적으로 다른 두 그룹보다 δ15N-NO3의 값이 낮게 도시된다. B그룹의 경우 지류주변의 지하수관정으로 지류의 유동방향에 따라 NO3의 농도가 감소하고, 하천 및 지류의 유입에 따른 희석작용을 통해서 δ15N값이 증가된 것으로 판단되며, 이는 탈질산화작용에 의한 질소동위원소의 분별작용으로 해석된다. C그룹은 농업활동 및 거주시설 증가로 인해 NO3의 증가 및 탈질산화작용에 따른 δ15N의 값이 증가되는 것으로 판단된다. C그룹은 인위적인 오염에 의해 NO3가 높게 나타나며, δ15N값은 8.69, 10.7‰로 자연토양, 생활하수 유래 등과 같은 유기질 질소의 영향으로 추정된다. A그룹은 하천수와 하천주변지하수로 δ15N값은 높고 NO3 값은 낮게 나타난다. Mariotti et al. (1988)에 의하면, 질산염의 농도가 감소하고 반대로 δ15N값이 증가하는 요인은 탈질산화작용, 지하수의 희석작용, 식물의 동화작용 및 미생물의 고정화에 의한 것으로 해석하였다. 연구지역 A그룹의 하천수의 경우 SEM을 통해서 하천수내 미생물의 존재를 파악한 결과 Fig. 6과 같이 하천수 상류(OR-3), 중류(OR-2), 하류(OR-1)에서 각각 미생물이 관찰되었다. 이들 미생물들은 탈질산화작용을 일으키는 미생물로 판단되며, 반응성이 큰 NO2를 분석한 결과 평균 1.26 mg/L로 나타난 것으로 보아, 탈질산화반응에 의한 δ15N값의 증가로 판단된다.

Fig. 5.Relationship between NO3 (mg/L) and δ15N (‰) at each sampling site.

Fig. 6.SEM images of bacteria in the river at the following sites: (a) upstream (OR-3), (b, c) midstream (OR-2) and (d) downstream (OR-1).

서로 다른 기원의 질산염들은 특징적인 질소동위원소조성을 보이기 때문에 동위원소는 여러 경우에서 질산염 오염원의 직접적인 파악의 수단을 제공하며, 이러한 특징을 이용하여 하상여과에 미치는 영향인자로 사용하였으나, 질소변환반응에는 암모니아화작용, 질산화작용, 탈질산화작용 및 동화작용 등 다양한 반응에 의해 변하기 때문에 정량적인 해석을 위해 추가적으로 산소, 탄소 동위원소 자료 및 주변 환경조사가 동시에 이루어진다면, 하상여과에 미치는 설계인자 해석에 정확성을 높일 수 있을 것으로 판단된다.

결 론

연구지역 지질은 충적층으로 이루어져 분급이 양호하며 투수성이 높아 지하수 산출성이 좋으나, 미고결 퇴적물로 이루어진 충적층에서의 지하수 수위는 지표면에 가깝게 상승되어 있어 지표수 오염원의 영향을 받기 쉽고, 하상여과 취수량에 영향을 미칠 것으로 판단된다. 수리화학적 특성은 배경수질의 차이를 구분한 결과 하천수와 지하수의 차이는 K, HCO3가 지하수에서는 Cl, NO3로 하천인근 지하수와 도심지 지하수로 구분되며, 하천 인근지하수 Cl 원소가 이상분포대를 형성되었다. 연구지역은 NO3의 농도에 대한 δ15N값에 따라 A, B, C그룹으로 구분되며, A그룹은 하천수의 기원 및 탈질산화반응, B그룹은 지류의 유동방향에 따른 탈질산화작용에 의한 분별작용, C그룹은 인위적인 오염과 자연토양, 생활하수 유래의 질산성질소로 판단된다. 기존의 하상여과에 미치는 설계인자를 파악하기 위해 대부분 수치해석 및 모델링을 통한 실내조사가 진행되었으나, 하상여과 시공지역의 하천수와 주변 지하수를 대상으로 지구화학적 해석 및 주변 환경조사 등과 같은 현장조사를 통해서도 하상여과에 미치는 설계인자 파악이 가능한 것으로 판단된다.