• 제목/요약/키워드: 유량의 관리범위

검색결과 141건 처리시간 0.035초

강우 및 기온 변동성에 대한 집중형 강우-유출 모형 민감도 평가 (Runoff projections under precipitation and temperature variability)

  • 우동국;조지현;이가림;이송희;노성진
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2022년도 학술발표회
    • /
    • pp.349-349
    • /
    • 2022
  • 기후변화가 고착화되면서 강우와 기온 변동으로 인한 가뭄 및 홍수 발생이 점차 증가하고 있다. 유역 단위의 유출량 예측은 기후변화로 인한 자연재해에 대비하기 위한 수자원 관리의 시작이라 할 수 있다. 하지만, 기후변화와 유출모형의 불확실성은 정확한 유출 분석을 어렵게 한다. 본 연구에서는 위에 제시된 불확실성을 완화하기 위하여 기후 스트레스 시나리오에 따른 두 개의 집중형 수문모형, 즉 airGR과 IHACRES를 이용하여 강우 및 온도 변화에 따른 유출량 변화를 비교, 분석하였다. 연구 대상 지역은 합천댐과 섬진강댐 유역이며, 각 모형을 NSE(Nash-Sutcliffe Efficiency) 및 KGE(Kling Gupta Effieicncy)를 목적함수로 하여 매개변수를 최적화를 하였다. 모형의 보정과 검정은 20년(1995년~2014년)의 유출 자료를 활용하였으며, 보정 및 검정 기간은 각각 6:4 비율로 설정하였다. 두 모형 모두 보정과 검정 기간에 비교적 높은 신뢰도(NSE>0.7, KGE>0.8)를 보여, 모형이 과거 사상을 재현하기에 적합하고, 모의 결과가 비교적 유사함을 확인하였다. 다음으로, 기후 스트레스 시나리오를 구축하기 위해 위 20년 입력 자료를 바탕으로, 강우는 -50%에서 +50%의 범위를 1%씩 구분하였으며, 기온은 0℃에서 8℃까지 0.1℃ 범위로 하여 총 8,181개의 시나리오를 구축하였다. 이후, 기후 스트레스 시나리오에 따른 두 모형의 풍수량, 최대 유량, 평수량을 비교, 분석하였다. 기후 스트레스 영향을 반영한 풍수량과 연최대유량의 경우, 강우 증가에 따른 유출 증가 등의 패턴은 두 모형에서 비슷하였으나, 강우와 기온의 변화가 커질수록 더욱 상이한 결과를 얻었다. 이와 반대로, 평수량의 경우 강우와 온도의 변화가 증가함에 따라 더욱 유사한 결과를 얻었다. 즉, 유역의 탄력적 기후변화 대응을 위해서는 모형의 불확실성에 대한 정량적 평가가 필요하다는 것을 시사한다.

  • PDF

산림 지역의 비점오염물질 유출 특성 (Characteristics of Non-Point Pollutants from Forest Landuse)

  • 김지연;김지현;정민경;지용대;황재엽;박수영;유재정;김태승
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2011년도 학술발표회
    • /
    • pp.287-287
    • /
    • 2011
  • 모니터링 자료의 부족으로 인하여 다양한 토지이용에서 발생하는 비점오염물질의 관리에 어려움을 겪고 있는 실정이다. 현재 환경부에서는 현행 토지계의 원단위를 세부적으로 분류하여 재산정하기 위하여 지목별로 장기 모니터링이 수행되고 있다. 특히, 산림 지역의 경우 도시 및 축산지역에 비하여 강우유출수의 농도는 낮더라도 유량적인 측면에 보았을 때 전체 수계에 대한 부하량 기여도는 매우 높다고 볼 수 있다. 따라서 본 연구는 장기모니터링의 일환으로 산림지역에 대한 비점오염물질 유출 특성을 파악하기 위하여 모니터링 및 분석을 실시하였으며, 이러한 결과는 향후 비점오염원 평가기반을 마련하고자 한다. 본 연구는 활엽수지역을 대상으로 2010년 4월부터 10월까지 총 16회에 걸쳐 모니터링이 수행되었으며, 시료의 성분 변화를 막기 위해 냉장기능이 있는 자동채수기를 이용하여 시료를 채취하였다. 수질분석항목은 BOD, COD, DOC, SS, T-N, $NO_3$-N, $NH_3$-N, T-P, $PO_4$-P로 총 9가지 항목을 분석하였다. 강우사상에 대한 모니터링 결과, 총강우량은 7.0~76.5mm, 강우지속시간은 1~30hr, 평균 강우강도는 0.88~18.50mm/hr의 범위를 보이고 있으며, EMC(Event Mean Concentration, 유량가중평균농도)결과 BOD는 0.4~2.4mg/L, T-N은 1.156~14.777mg/L, T-P는 0.009~0.562mg/L인 것으로 나타났으며, SS는 1.8~71.9mg/L 로 비교적 높은 값을 나타내는 것으로 분석되었다. 농도 변화 및 유출경향의 패턴을 볼 때, 유량이 증가함에 따라 농도도 점점 증가하여 첨두유량이 발생된 후 감소하는 경향을 나타내는 것으로 분석되었다. 또한 우리나라의 경우, 시험유역을 대개 산지 소유역에 설치하는 경우가 많아서 일반적으로 지연시간이 짧은 경우가 많기 때문에 이 지역 역시 강우가 내린 후 계류유출량의 증가에 영향을 주는 강우의 유출속도는 비교적 빠른 것으로 나타났다. 그리고 단기 수문곡선상에서 강우량이 많을 시 유출이 빠르게 일어나 첨두 유량에 도달하는 시간이 짧고, 강우량이 적을 시에는 첨두 유량의 출현시간이 늦어지는 것을 볼 수 있었다.

  • PDF

농업용 저수지에서의 질화세균의 계절적인 변화 (Seasonal Variations of Nitrifying Bacteria in Agricultural Reservoir)

  • 이희순;이영옥
    • 생태와환경
    • /
    • 제35권3호통권99호
    • /
    • pp.152-159
    • /
    • 2002
  • 농업용 저수지로 축조된 문천지에서 질소순환에 관여하는 질화세균의 연중변화를 현장그대로의 세균상을 반영하는 FISH법으로 정량분석하였다. 아울러 질화세균군의 계절적인 변화와 이화학적인 환경요인과의 상관성을 규명해보기 위해 수온,pH,용존 유기물(DOC),chi-a등이 측정되었다. 총세균수는 정점 1에서 3월에 급격히 상승한 경우를 제외하고는 정점간의 큰 차이가 없이 $0.8{\sim}1.5{\times}10^6\;cells/ml$ 범위에서 변화하였고 eubacteria의 총세균수에 대한 비율는 $44.9{\sim}79.5%$, 질화세균이 eubacteria에서 차지하는 비율은 $1.0{\sim}7.4%$에 그쳤다. 암모니아산화세균은 큰 계절적인 변동이 없이 $1.1{\sim}3.0{\times}10^4\;cell/ml$범위에서 변화한 반면에 아질산산화세균은 정점에 관계없이 겨울철에 증가하는 양상을 나타냈다($1.3{\sim}5.7{\times}10^4\;cells/ml$).또한 암모니아산화세균보다 아질산산화세균이 대체적으로 더 많은 개체수를 나타냈다. 아울러 이화학적인 환경요인 중, 수온은 총세균수와 양의 상관성 (r=0.355, p<0.05), 그리고 아질산산화세균과는 음의 상관성 (r= -0.416, p<0.05)을 보였다. 마찬가지로 아질산산화세균은 pH와도 음의 상관성 (r= -0.568,p<0.01)을 나타냈다. 아울러 DOC는 총세균수(r = 0.586,p<0.000),  eubacteria (r = 0.448, p<0.01)와 긴밀한 상관관계에 있는 것으로 나타났다. 따라서 문천지에서는 수온, DOC 그리고 pH가 수계 세균군집구조를 변화시키는 중요한 요인으로 작용함을 알 수 있었다. 또한 FISH법은 느리게 성장하는 질화세균군을 정량 분석할 수 있는 유용한 기법으로 평가된다.기간동안 강우의 빈도 및 강도는 하천의 육수학적 현상 변화와 패턴에 중요한 요소로서 작용하며, 이는 동아시아몬순기후대에 속하는 하천들에서 유사하게 발생하는 하천의 중요한 특성으로 사료된다.다. 정화효율과 수리학적조건간의 상관계수($R^2$)는 수리학적 체류시간과 0.016-0.731,일처리유량과는 0.015-0.868을 나타내었으며, 시간당 정화량과 수리학적 조건간의 상관계수($R^2$)는 수리학적 체류시간과는 0.173-0.763,일처리유량과는 0.209-0.770의 범위를 나타내었다. 정화효율과 수리학적 부하조건간의 상관계수($R^2$)Tt 0.5 이상을 나타내는 각 수생식물 습지별 수질항목은 체류시간과 일처리유량에 대해각각 20%,정화속도와 수리학적 조건간의 상관계수는 체류시간에 대해 53%, 일처리유량에 대해73%가 0.5이상을 보이고 있어 시간당 정화량과 수리학적 조건간의 상관관계가 정화효율과의 상관관계보다 좀더 유의성 있게 나타났다. 이것은 높은 수리학적 부하조건이 영양염류 등의 정화효율에는 크게 영향을 미치지 않음을 보여주고 있으며, 따라서 비교적 낮은 농도의 영양염류를 가지고 있고, 많은 처리수량을 요구하는 부영양화된 저수지의 수질개선을 위해서는 높은 수리학적 부하조건에서 시간당 정화량을 늘리는 관리방법이 경제적이며, 이에 초점을 맞추어 나가야 할 것으로 사료된다.도로 검출되어 이를 고려한정수공정의 관리가 필요하리라 생각된다.$yr^{-1}$)의 71%가 감소되어야 할 필요성이 제기되었다. 또한 여름철 심층 산소 고갈이 야기되었고, 이 시기에 퇴적물로 용출된 인이 식물플랑크톤 성장에 이용될 수 있기때문에 퇴적물에 대한 관리도 수행될 필요가 있다.rsity)지수는 2000년 7월에 2.22로 가장 높았으며, 생물량도 조사기간중 36,640 cells ${\cdot}\;mL^{-1}$로 가장 많았다.

배수시설을 고려한 토석류 피해지 분석 (Analysis of Debris Flow Damaged Areas Considering Drainage Facilities)

  • 강배동;전계원;장창덕;조항일
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2023년도 학술발표회
    • /
    • pp.140-140
    • /
    • 2023
  • 전 세계적인 기후변화로 인해 태풍과 집중호우의 발생빈도와 규모가 지속적으로 증가하고 있다. 이러한 태풍과 집중호우는 국토의 64%가 산지로 이루어져 있는 우리나라에서 토석류 재해를 발생시켜 하류부에 많은 피해를 입힌다. 토석류는 계곡을 따라 흐르는 특성을 가지고 있어 콘크리트 사방댐, 슬릿트형 사방댐과 같은 횡단구조물과 기슭막이와 같은 종단구조물의 설치는 토석류피해 저감에 매우 유효하게 작용하고 있다. 따라서 토석류 재해의 발생 위험이 높은 유역에서의 토사유출량의 예측과 재해저감시설의 규모와 형상, 위치 설정은 토석류 피해를 예방하기 위한 중요한 요소이다. 본 연구에서는 2019년 10월 토석류 재해가 발생한 지역의 피해를 분석하기 위해 2차원 수치모형인 Hyper KANAKO 모형을 이용하였으며, 해당 지역의 수치지형도를 이용하여 수치표고모델(DEM)과 토석류 피해를 저감하기 위한 배수시설을 고려한 지형자료를 제작하였다. 그리고 모형의 주요 입력변수인 첨두유량은 토석류 피해가 발생하였을 때의 강우 자료와 현장조사 자료를 이용하여 산정하였다. Hyper KANAKO 모형의 결과로 나타나는 확산범위, 퇴적량 등을 분석하여 배수시설의 유무에 따른 토석류 피해 저감효과를 분석하였다.

  • PDF

국내 고속도로 강우 유출수의 EMCs 및 유출 부하량 산정 (Estimation of Pollutant EMCs and Loadings in Highway Runoff)

  • 김이형;고석오;이병식;김성길
    • 대한토목학회논문집
    • /
    • 제26권2B호
    • /
    • pp.225-231
    • /
    • 2006
  • 최근 들어 오염총량관리제의 시행과 함께 비점오염원 관리의 중요성이 커지고 있다. 오염총량이란 하천수질의 개선을 위하여 도입하였으며 유역이란 개념을 도입하였다. 하천수계로 유입되는 오염물질 부하량은 점오염원 부하량, 비점오염 부하량 및 안전율의 합으로 산정이 되고 있다. 그러나 일반적으로 유역 및 강우 특성에 기인하여 유출되는 비점오염원의 특수성 때문에 안전율을 보통 10% 정도로 두고 있다. 이러한 불확실성은 다양한 요인에 기인하는데, 모니터링을 통한 기초자료의 부족이 가장 원인으로 보고되고 있다. 따라서 본 연구는 다양한 비점오염원 중에서 포장율이 높아 오염물질의 유출이 높은 고속도로 지역에서의 오염물질 유출 부하 원단위 산정을 위하여 가장 필수적으로 요구되는 오염물질의 EMC 및 유출부하량 산정을 위하여 수행되었다. 고속도로의 경우, 높은 포장율과 많은 자동차의 운행으로 인하여 건조기간 동안 오염물질들의 축적이 높은 지역이며, 강우시 집중적으로 인근 수계로 유입되기에 비점오염원 관리에서 매우 중요한 토지이용으로 고려되고 있다. 그러나 현재 국내에서는 비점오염원에 관한 기초 연구가 미흡하여 고속도로 등을 포함한 각종 토지이용에서의 비점오염원에 관한 기초 자료의 부족이 매우 심각한 상황이다. 본 연구를 수행하기 위하여 국내 5개의 고속도로 지점에서 모니터링이 수행되었으며, 각 지점별 자동 유량 및 강우량 측정계를 설치하여 자료를 수집하였다. 모니터링 자료는 고속도로에서의 비점오염물질 EMC 및 부하량 산정을 위하여 사용되었다. EMC에 관한 95% 신뢰구간을 살펴보면, TSS의 경우 45.52-125.76 mg/L, COD는 52.04-95.48 mg/L, TN은 1.77-4.48 mg/L 그리고 TP의 경우 0.29-0.54 mg/L의 범위로 나타났다. 부하량에 관한 95% 확신범위는 TSS가 $712.7-2,418.4mg/m^2$의 범위를 보였으며, COD의 경우 $684.1-1,779.6mg/m^2$의 범위로 나타났다.

국내 LID 기법 적용확대의 문제점 및 개선방안 (Issues and Improvements in Extended Application of LID technologies)

  • 최종수;현경학;이정민;강명수;정승권
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2012년도 학술발표회
    • /
    • pp.242-246
    • /
    • 2012
  • 우리나라의 경우 비점오염원에 의한 오염부하 기여도가 점차 증가하면서 '06년 관련 법규 개정을 통해 일정 규모 이상의 개발사업에 대해서는 비점오염저감시설 설치를 의무화하였으며, '11년 기준 국내에는 약 6,000여개 이상의 저감시설이 설치되었다. 해당 시설에 대한 모니터링 결과 장치형에 비해 자연형이 처리효율 및 유지관리 측면에서 유리한 것으로 보고되면서 환경부는 개발사업에 대해 저영향개발기법(LID)의 적용을 의무화하는 방안을 추진 중이다. 개발사업에 있어 LID는 토지이용기법을 통해 우수 유출과 오염원 유출을 저감하고자 하는 기법으로, 비용 및 유지관리 측면에서는 유리하나 기존 단위시설 적용기법에 비해 계획적인 요소가 강해 적용 확대 및 의무화가 용이하지 않은 한계가 있다. 따라서, 본 연구에서는 LID 적용 확대와 관련된 문제점을 분석하고 적용 의무화를 위해 필요한 대안을 모색해 보고자 하였다. 이를 위해 해당 기법 자체가 가지는 한계와 제도적인 문제점을 분석하였으며, LID 기법의 확대 적용을 위해 필요한 대안을 제안하였다. 분석 결과, LID 기법은 기존 장치형 시설에 비해 계획기법의 측면이 커 기술개발이 어려울 뿐만 아니라, 저감시설의 적용여부에 대한 검증도 곤란한 것으로 분석되었다. 또한, 처리효율 계량화가 어려워 오염총량제와 연계된 오염부하 삭감효과 인정이 곤란해 확대 적용에 한계가 있는 것으로 나타났다. 상기의 분석결과에 근거할 때, LID 기법의 확대 적용은 국내의 제도적 여건 및 LID 기법 자체의 한계를 감안하여 단계적으로 추진되어야 할 것으로 판단되며, 조속한 확대 시행을 위해서는 다음의 방안 마련을 제안한다. 첫째, 적용 여부에 대한 판단이 모호한 단점을 보완하기 위해 적정 적용여부를 심의할 수 있는 인증제도 도입을 제안한다. 이 제도는 친환경건축물인증제도 등 기존 유사 제도에 항목을 추가하는 방법과 별도의 제도를 신설하는 방법이 가능할 것이다. 둘째, 처리효율 정량화를 위해 기존 농도 위주의 처리효율 개념을 벗어나 유량 개념의 처리효율 산정이 필요하다. 즉, 처리용량 범위내에서 외부에 유출없이 우수저류 및 침투가 가능한 경우 전량 처리되는 것으로 간주하는 방식이다. 연구결과로 제안한 LID 기법 적용에 대한 인증제도 마련과 유량개념의 처리효율 산정방안이 도입된다면, 보다 조속한 시일내 LID 기법의 확대 적용이 가능할 것이다.

  • PDF

필지내 벼 수량의 공간변이 해석 (Spatial Variability Analysis of Paddy Rice Yield in Field)

  • 이충근;우메다미키오;정인규;성제훈;김상철;박우풍;이용범
    • Journal of Biosystems Engineering
    • /
    • 제29권3호
    • /
    • pp.267-274
    • /
    • 2004
  • 논 포장에 있어서 벼 수량의 연차별 공간변이를 살펴보기 위해 포장면적, 조사 년도, 조사방법 및 재배품종 등이 각각 다른 안성포장, 수원포장, 쿄토포장에서 수집된 수량정보를 이용하여 기술통계와 함께 공간통계학적 개념을 도입하여 해석한 결과 아래와 같은 결과를 얻었다. (1) 포장면적, 조사 년도, 조사방법 그리고 품종에 따라 한 필지의 포장 내 수량분포는 100-946 kg/10a의 큰 변이를 보였고, 최대수량과 최소수량의 차이는 272-653 kg/10a를 보였으며, 변이계수는 5.9-22.4%의 큰 공간변이를 보였다. (2) 수량빈도 분포로부터 보면 약 90 % 이상의 수량이 10 a당 350-850kg의 수량범위를 나타내었다. 이 결과에 따르면 4조식 콤바인의 평균 작업속도를 0.8 m/s로 가정하고 평균수량이 350-850kg/10a이라고 한다면, 수확량 모니터링 시스템 개발에 있어서 초당 곡립유량은 0.34-0.82kg으로 산정 되므로 곡물유량센서를 개발함에 있어서 이 수량범위를 작동범위에 포함시켜야 할 것으로 판단되었다. (3) 수량정보의 공간의존성 해석결과에 의하면 수량의 공간구조를 나타내는 Q값은 0.20-0.97이며, 그 의존거리를 표시하는 범위는 6.9-53.3m로 나타났다. 따라서, 이 시험결과에 의하면 수량정보의 이론적인 샘플링 간격은 6.9 m 이내가 적정하다고 판단되었다. (4) 논 포장에 있어서 유사한 포장관리를 함에도 불구하고 수량의 연차별 차이가 발생하는 또 다른 주원인은 연차별 기상조건과 병$.$해충 발생 정도로 생각되며, 정밀농업을 논 포장에 적용함에 있어서 포장정보 외의 이러한 외적인 정보들도 데이터베이스화할 필요성이 있었다.

SWAT 모형의 수문.수질 모듈 개선에 의한 BOD 모의 (A BOD Simulation by Improved Hydrology and Water Quality Module in SWAT)

  • 김남원;신아현
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2009년도 학술발표회 초록집
    • /
    • pp.2068-2072
    • /
    • 2009
  • 본 연구에서는 국내에서 널리 사용되는 대표적 유역모형인 SWAT 모형의 한계를 인식하고 유출해석모듈과 하도수질모듈을 개선한 SWAT-K(Korea) 모형이 현재 오염총량관리 대상물질로 지정되어 있는 BOD 모의에 미치는 영향을 분석하기 위하여 충주댐유역을 대상으로 각 개선모듈별 BOD 모의 결과를 분석하였다. 유출해석모듈의 개선결과 하도의 흐름을 따라 이동하는 BOD 항목의 모의 특성 상 유출 구조의 개선에 의하여 BOD 부하량의 지체가 감소되며 첨두 부하량이 증가한 것을 확인할 수 있었다. 하도수질모듈의 개선결과 저유량 시 매개변수 산정에도 불구하고 증가하지 않던 값들이 모형의 개선에 의하여 정상적인 범위로 상승한 것을 확인할 수 있었다. 전체적인 개선에 의하여 각 모듈별 개선효과와 함께 충주댐 지점의 개선 전 후 일대일 비교를 수행한 결과 결정 계수($R^2$)가 0.54에서 0.88로 상승하며 개선 후 모의값이 실측값에 근사하게 모의되는 것을 확인할 수 있었다.

  • PDF

과류차단형 LPG용기용 밸브의 차단성능에 관한 실험적 연구 (A Experiment Study of Cut-off Capability for L.P Gas Cylinder Cut-off Type Over-flow Valve)

  • 김진구;노국장;허문희;이성일;임종국;이인찬
    • 대한안전경영과학회:학술대회논문집
    • /
    • 대한안전경영과학회 2004년도 추계학술대회
    • /
    • pp.53-60
    • /
    • 2004
  • 본 연구에서는 현재 일반가정이나 영업소에서 사용중인 액화석유가스 용기용밸브의 문제점을 제시하고, L.P 가스사고를 예방할 수 있는 LPG용 과류밸브의 과류차단장치를 개발${\cdot}$제작함에 있어 고압가스 안전관리기준 통합고시 ${\lceil}$제12장 4절 과류차단형 액화석유가스 용기용밸브${\rfloor}$에서 규정한 과류차단 성능시험을 실시하는데 필요한 여러 가지 변수를 선정하여 차단장치의 형상에 따른 최적의 차단 유량 범위를 찾고자 실험을 실시하였으며, 보다 궁극적으로는 일반 용기용밸브가 가진 문제점을 조속히 해결하고 진정한 의미의 안전밸브를 개발하는데 목적이 있다.

  • PDF

계절 ARIMA모형을 이용한 과거 유입량 분석기간 적용성 연구 (The past Inflow data Period Validit Analysis Using Seasonal ARIMA Model)

  • 김건순;이충대
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2010년도 학술발표회
    • /
    • pp.1410-1414
    • /
    • 2010
  • 최근 들어 가뭄과 국지성 호우 등의 기상이변이 지속적으로 발생하고 있으며, 이는 국민 삶의 발전과 향상에 밀접한 관계가 있는 것으로 전세계적으로 이에 대한 관심이 증가하고 있는 추세이다. 특히 댐의 효율적 관리와 안정적인 운영은 홍수피해 방지, 안정적인 용수공급과 같은 국민 생활과 밀접한 관계를 가지고 있어 수자원의 효율적인 운영과 이용은 장기적인 관점을 통하여 수립해야 한다. 이와 같이 댐 유입량의 예측은 유출모형의 목적 중 중요한 부분으로 확정론적 모형이 시 혹은 일유량과 같은 매우 짧은 시간의 유출을 예측하는데 주로 사용되지만 이는 매개변수의 추정이 불가능하거나 실제유역에서의 측정이 불가능 할 경우에는 모형적용에 한계가 있다. 이에 반해 추계학적 모형에 의한 유출예측은 장기간의 유출을 과거자료의 통계학적 특성변수를 매개변수로 하여 예측하는 방법으로 모형의 적용에 필요한 매개변수가 적어 그 적용성이 간편한 장점이 있다. 본 연구에서는 계절형 ARIMA모형을 적용하여 과거자료의 적용범위, 매개변수의 산정, 적합성 판정에 대하여 판단하고, 이 모형이 월유입량의 예측에 적합한지를 검토하였다.

  • PDF