• Title/Summary/Keyword: 수문특성분석

Search Result 1,284, Processing Time 0.038 seconds

Baseflow Separation from Complex Hydrograph using Eckhart Filter (Eckhart filter를 이용한 복합수문곡선에서의 기저유출분리)

  • Kang, Boo-Sik;Moon, Su-Jin
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2009.05a
    • /
    • pp.1292-1295
    • /
    • 2009
  • 일반적으로 유역의 기저유출 산정에는 적용의 간편성을 이유로 수평직선분리법, N-day법등을 주로 사용해 왔으며 이를 단기호우사상에 대한 모의에 적용해 왔다. 그러나 수평직선분리법이나 N-day에의한 기저유량의 산출은 연구자의 주관성이 반영될 수 있는 가능성이 다분하며 총 유출에 대한 기저 유출의 기여가 상대적으로 크게 되는 장기유출모의에 이용하기에는 효율적이지 못할 뿐 아니라 분석을 수행하는 사람에 따라서 그 결과에 많은 차이가 있을 수 있다. 그래서 최근에는 이러한 문제를 해결하기 위해서 다양한 'Digital Filtering' 방법이 수문곡선 분리에 많이 사용되고 있다. Eckhardt 필터는 출력의 신호대 잡음비가 최소화를 목적함수로 사용하는 최적선형필터의 형태를 갖고 있으며, BFImax 변수 값을 이용하여 수문분석 시 대수층별 특성을 반영할 수 있는 장점을 지니고 있다. 따라서 BFImax 변수 값에 따라서 분리된 직접유출과 기저유출 값에는 상당한 차이가 발생할 수 있는데, 정확한 수문분석을 위해서는 연구대상 유역 내 대수층의 특성에 가장 부합되는 BFImax 변수 값을 구하여야 하는 어려움이 있다. 이러한 문제를 해결하기 위하여 감수곡선 분석방법과 최적화 기법을 이용하여 BFImax 변수 값을 결정해 주는 모듈을 사용하고 있다. 금강유역의 수문곡선분리를 위해 먼저 금강 유역 내 관측소의 강우 및 수위 일 자료를 수집하고 수위자료는 동향, 옥천, 천천, 청성, 호탄지점의 자료를 이용하였다. 관측단위는 일 단위를 사용하였다. 분리된 수문곡선에 대한 검증은 지하수관측이나 중간유출에 대한 관측에 한계가 있기 때문에 정량적인 검증은 쉽지 않은 상황이고 수문곡선에 대한 시각적 판단에 의지하고 있다. 하지만 효과적인 기저유출분리가 가능해진다면 하천 건천화에 대한 정확한 진단과 함께 기후변화에 따른 지하수 및 지표수 영향에 대한 보다 신뢰도 있는 전망이 가능해지므로 매우 중요한 연구 분야라 판단된다.

  • PDF

An Analysis of the Temporal Pattern according to Hydrologic Characteristics of Short-Duration Rainfall (단시간강우의 수문학적 특성에 따른 시간분포 분석)

  • Lee, Jung-Sik;Shin, Chang-Dong;Chang, Jin-Wook
    • Journal of the Korean Society of Hazard Mitigation
    • /
    • v.6 no.3 s.22
    • /
    • pp.57-68
    • /
    • 2006
  • The objective of this study is to analyze the temporal pattern characteristic of short-duration rainfall defined as a rainfall durations of 6 hours or less by the Huff's 4th quartile distribution. To analyze the temporal pattern characteristic of short-duration rainfall, the rainfall data are classified by rainfall duration and rainfall type(Changma, Typhoon, Severe rain storm, Frontal storm) and change of rainfall segment. Also, the results of this study compared with result of research work of Korea Institute of Construction Technology(1989) and Ministry of Construction & Transportation(2000). The conclusions of this study are as follows; (1) Short-duration rainfall with duration of 6 hours or less is found to be most prevalent frist-quartile storms. (2) In the case of rainfall type, Changma and Severe rain storms and Frontal storm is found second-quartile storms, and Typhoon is found third-quartile storms. (3) In the result by change of sixth segment storms, the type of temporal pattern of rainfall is found to be most prevalent two sixth parts, (4) Comparative analysis of the results shows that shapes of the dimensionless cumulative curves and values are different from those of existing researches.

An operation of Pangyo experimental catchment to assess the variation of stream hydrology and water quality by urbanization (도시화에 따른 하천 수문.수질 변동 영향 평가를 위한 판교 시험유역의 운영)

  • Jang, Cheol-Hee;Kim, Hyeon-Jun;Kim, Chul-Gyum;Noh, Seong-Jin
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2006.05a
    • /
    • pp.1865-1869
    • /
    • 2006
  • 도시화는 수문학적으로 산림이나 농경지와 같은 투수지역을 건물, 도로 등의 불투수지역으로 변화시키는 것이며, 이로 인하여 홍수파의 도달시간이 줄어들고 첨두유량이 증가하는 등의 수문변화를 수반하게 된다. 도로나 건물 등이 대부분을 차지하고 있는 도시지역에서는 지표면이나 식생으로부터 대기중으로 방출되는 증발산량이 농촌이나 산림지역보다 상대적으로 적으며, 강우시 토양중의 침투량과 지표면의 저류량도 도시지역에서는 매우 적게 나타난다. 본 연구에서는 도시화로 인하여 대규모의 토지이용 변화가 예상되는 신도시 개발예정지구인 판교 운중천 유역을 시험유역으로 선정하여 장기적인 수문/수질 모니터링을 수행하고 있으며, 향후 판교 시험유역에서 계측된 자료와 유역수문모형을 활용하여 개발 전후의 홍수 및 장기유출 특성을 분석하여 도시화가 하천의 수문 및 수질에 미치는 영향을 정량적으로 비교 분석하고자 한다. 이를 위해 2004년과 2005년에 시험유역내의 운중천과 금토천에 매송2교, 삼평교, 판교교, 운중저수지, 저수지 용수로, 내동교 등 총 6개의 수위관측소와, 매송2교와 내동교의 2개 우량관측소를 설치하여 실시간 계측을 수행하고 있으며, 인근의 기존에 운영되고 있던 탄천 본류의 성남과 궁내 수위관측소와 운중천의 낙생 우량관측소 자료를 활용하여 유역 수문자료를 수집하고 있다. 수집된 자료를 이용하여 장단기 유출분석을 수행한 결과, 택지개발공사가 시작되기 전인 현재로서는 택지개발 대상지구인 운중천 유역이나 비대상지구인 금토천 유역, 그리고 하류의 운중천 본류 유역 모두 강우시나 무강우시 비슷한 유출 특성을 가지고 있는 것으로 나타났다. 향후 도시화가 점차적으로 진행됨에 따라서 이러한 유출 특성의 변화 양상을 파악해 나간다면, 도시화에 따른 유출 특성의 변화는 물론 수질 변화에 대해서도 정성적, 정량적인 규명이 가능할 것으로 생각된다.TEX>의 범위이고, 2차 처리수의 유입수의 T-N, T-P 농도와 유사하였다.적인 방법론을 제시할 수 있을 것으로 사료된다.첨두홍수량을 저류하기 위해서 상대적으로 넓은 저류면적이 필요한 것으로 나타난다. 대등한 수위감소값의 홍수저감효과를 발휘하기 위해서 본 연구에서는 On-Line 저류지 면적은 Off-Line 저류지에 비 두배 이상이 필요한 것으로 보여졌다.들에 관한 정보는 종종 현장관측에서 조차 무시되는 경우가 많다. 이에 본 연구에서는 수질모형의 매개변수 중 특히 수리특성에 관련된 매개변수들이 수질에 미치는 영향을 파악하는 것을 목적으로 하고 있다. 이를 위해 적용된 수질모형은 QualKo를 사용하였으며, 대상 하천은 낙동강 본류 경남구간 시점 부근인 회천 합류 전부터 낙동강 본류 경남구간 종점 부근인 밀양강 합류 전까지의 경남 오염총량관리 기본계획 시 구축된 모형 매개변수를 바탕으로 분석을 수행하였다. 일차오차분석을 이용하여 수리매개변수와 수질매개변수의 수질항목별 상대적 기여도를 파악해 본 결과, 수리매개변수는 DO, BOD, 유기질소, 유기인 모든 항목에 일정 정도의 상대적 기여도를 가지고 있는 것을 알 수 있었다. 이로부터 수질 모형의 적용 시 수리 매개변수 또한 수질 매개변수의 추정 시와 같이 보다 세심한 주의를 기울여 추정할 필요가 있을 것으로 판단된다.변화와 기흉 발생과의 인과관계를 확인하고 좀 더 구체화하기 위한 연구가 필요할 것이다.게 이루어질 수 있을 것으로 기대된다.는 초과수익률이 상승하지만, 이후로는 감소하므로, 반전거래전략을 활용하는 경우 주식투자기간은 24개월이하의 중단기가 적합함을 발견하였다. 이상의 행태적 측면과 투자성과측면의 실증결과를 통하여 한국주식시장에 있어서 시장수익률을 평균적으로 초과할 수 있는 거래전략은 존재하므로 이러한 전략을 개발 및 활용할 수 있으며, 특히, 한국주식시장에

  • PDF

The Analysis of Hydrological Response Structure Based on Spatial Correlation of Extracted Geomorphic Variables by Using DEM (DEM에 의해 추출된 지형인자의 공간상관성을 기반으로 한 수문학적 응답구조의 해석)

  • Kim, Joo-Cheol;Choi, Yong-Joon;Kim, Jae-Han
    • Journal of Korean Society for Geospatial Information Science
    • /
    • v.18 no.2
    • /
    • pp.69-78
    • /
    • 2010
  • The hydrological-geomorphic character is closely related with runoff in basin. A development of GIS greatly helps investigating about mechanism between theirs. We analyze local slope and hillslope length which are related with hydrological response. But variation of these geomorphic variables has very wide range at each pixel. So there's a limit as to use directly. Therefore we investigate a relation between hydrological response and distributed geomorphic variables according to statistical character of distributed map considering spatial correlation. As a result, the local slope affects peak discharge, and the hillslope length affects peak discharge and time, mean and variance of hydrological response. Henceforth these hydrological-geomorphic analyze methods can be improved that hydrology response is directly analogized with DEM data.

Analysis evapotranspiration characteristics of mountain and farmland (산지와 농경지의 증발산량 특성 분석)

  • Lee, JungHoon;Kim, Kiyoung;Lee, YeonKil;Jung, SungWon
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2017.05a
    • /
    • pp.235-235
    • /
    • 2017
  • 수문 순환과 물 수지에 관한 연구는 강수량, 지표유출량, 지하수, 토양수분 및 증발산량 등에 대한 관측이 이루어질 때 실제로 규명될 수 있다. 하지만, 수문 순환과 물수지 평가에 중요한 부분을 차지하는 증발산량의 경우 관측값보다 단순한 가정이나 경험식에 의한 추정값을 사용하고 있다. 따라서, 수문 순환과 물수지의 정량적인 분석을 위해서는 수문 순환 과정에서 상당부분을 차지하는 증발산량의 측정(실측)과 자료의 축척이 필요한 실정이다. 본 연구는 국토교통부의 기초수문자료 구축사업의 일환으로 수행되었으며, 수문자료의 다양화 목적을 가지고 에디공분산 기술을 사용하여 증발산량을 직접 관측하고 있다. 관측지점은 한반도의 약 70%를 차지하는 산림지 중 대표적 식생 기능 형태인 혼효림으로 구성된 지점(설마천 관측소, 2007년 8월부터)과, 인위적인 관개가 이루어지는 농경지(청미천 관측소, 2008년 8월부터)에서의 증발산량 측정을 수행하였다. 관측된 자료를 활용하여 관측소별로 연도별 증발산량 특성을 분석한 결과는 다음과 같다. 설마천 관측소(주변식생 : 혼효림)에서 산정된 연증발산량은 2008년 471.7mm, 2009년 408.4mm, 2010년 489.4mm, 2011년 387.0mm, 2012년 323.3mm, 2013년 293.3mm, 2014년 360.9mm, 2015년 419.6mm, 2016년 566.9mm이고, 발생한 강수량 대비 증발산 비율은 18.9%~56.2%범위로 산정되었다. 청미천 관측소(주변식생 : 농경지)에서 산정된 연증발산량은 2009년 571.8mm, 2010년 650.6mm, 2011년 523.9mm, 2012년 509.8mm, 2013년 467.9mm 2014년 533.9mm, 2015년 600.5mm, 2016년 588.0mm이고, 발생한 강수량 대비 증발산량의 비율은 25.6%~71.4%범위로 산정되었다. 강수량 대비 증발산량 비율의 최대값은 설마천 관측소는 2014년, 청미천 관측소는 2015년에 발생하였다. 평균 증발산량 비율은 산림지인 설마천 관측소보다 논경지인 청미천 관측소가 평균 13.3%정도 높은 특성을 보였다.

  • PDF

Evaluation for the Effectiveness of Distributed Hydrologic Model in Urban Sub-drainages (도시하천 소배수구역에서의 분포형 수문모형의 효율성 평가)

  • Kim, Chung-Soo;Kim, Hak-Su
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2007.05a
    • /
    • pp.1829-1833
    • /
    • 2007
  • 최근 강우-유출 모형은 유역의 지형적, 물리적 상황을 최대한 반영한 입력자료를 이용할 수 있고, 최종 유출지점에서의 해석 결과뿐만 아니라 해석하고자 하는 유역내의 시간적, 공간적 수문 요소 분포 특성을 이해할 수 있도록 물리적 이론에 기초한 분포형 수문모형을 개발하고 지속 발전시키고 있다. 분포형 수문모형은 집중형 수문모형에 비해 매개변수 갯수를 현저히 줄일 뿐만 아니라 실지 지형자료에 입각한 매개변수 선정이 가능하다. 분포형 지형자료와 강우자료로부터 추정된 초기 매개변수의 값에 의한 모의가 매개변수가 개념적, 경험적 의미가 큰 집중형 모형에 비해 정확하기 때문에 미세한 매개변수의 조정만으로도 유역의 유출량을 모의할 수 있다. 본 연구에서는 이러한 여러 가지 장점을 지닌 분포형 수문모형이 유역크기가 작고, 시간 공간적으로 지형입력자료의 변동성이 작은 도시하천 소배수구역에서의 모의 효율성 및 적용성을 평가하였다. 레이더강우를 이용한 분포형 수문모형으로서 미국 University of Oklahoma에서 개발한 $Vflo^{TM}$을 사용하여 중랑천 유역의 3개 소배수구역(월계1 배수구역, 군자 배수구역, 어린이대공원 배수구역)에서 유출 모의를 수행하였다. 각 소배수구역 특성별 모의 결과 비교, 분석을 수행하였다. 이를 통해 도시하천 소배수구역에서의 유출 모의 및 분포형 수문모형 선정에 적절한 근거로 활용될 수 있을 것으로 판단된다.

  • PDF

Runoff Characteristics Analysis in the ImJin-River Basin (2005년 임진강 유역의 유출특성 분석)

  • Kim, Dong-Phil;Jung, Sung-Won;Lee, Seok-Ho;Lee, Sang-Cheol
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2006.05a
    • /
    • pp.1756-1761
    • /
    • 2006
  • 본 연구의 목적은 임진강 유역을 대상으로 2005년에 발생한 유출특성을 분석하는데 있다. 임진강 유역은 유역전체 면적중 2/3이상이 북한지역에 위치하고 있어서 전체적인 수문관측 자료를 수집하기 어려운 상황에 있으나, 가용한 수문관측 자료 및 유량측정성과를 이용하여 5개 지점에 대한 유출특성을 분석하였다. 임진강 유역에서는 30분 단위의 연속적인 우량 및 수위관측(건설교통부)을 하고 있으며, 기상청에서는 1시간 단위의 기상관측을 하고 있다. 우량 및 수위자료의 수집과 검토를 통하여 유역면적우량과 5개 수위관측소의 수위를 최종 확정하였다. 그리고 2005년에 측정한 유량측정성과를 통하여 유량측정성과에 대한 불확실도, 기본 수리특성 분석, 수위-유량관계곡선식을 개발하여 유량을 산정하였다. 산정된 유량은 월별, 주요 호우사상별 유출률 분석을 통하여 합리적인 유출률 범위내에서 재 조정절차를 거쳐 최종 유량을 확정하였다. 산정된 유출률은 57.0%(영중)${\sim}$72.1%(적성)의 유출률을 보였으며, 2004년 임진강 유량측정용역 보고서의 유출률과 비교한 결과 일부지점을 제외하고는 비교적 일치하는 것으로 나타났다. 임진강 유역의 일부인 설마천 시험유역에서의 2005년 연간 유출률은 60.9%로 상기의 값은 매우 타당하다고 할 수 있다. 그리고, 남북한 수문자료의 공유 및 기술적 교류, 유역 전체의 정밀한 수문관측 및 신뢰성 있는 수문자료 생성이 이루어진다면 보다 정확한 유출특성 거동을 분석할 수 있을 것이며, 효과적인 치수 및 이수계획의 수립 등 수자원 개발을 위한 기반을 확보할 수 있을 것이다.수리모형실험과 연계하여 댐 설계에 효율적으로 사용될 수 있을 것으로 전망된다.지는 것으로 평가되었다. 그러나 본 연구결과를 통하여 투수성 포장과 지하수에 관련된 매개변수의 집적과 분석결과는 현장기술 적용 시 매개변수의 유용한 선택과 도시유역의 물 순환 건전화 대안기술 적용에 효과적인 방법론을 제시할 수 있을 것으로 사료된다.첨두홍수량을 저류하기 위해서 상대적으로 넓은 저류면적이 필요한 것으로 나타난다. 대등한 수위감소값의 홍수저감효과를 발휘하기 위해서 본 연구에서는 On-Line 저류지 면적은 Off-Line 저류지에 비 두배 이상이 필요한 것으로 보여졌다.들에 관한 정보는 종종 현장관측에서 조차 무시되는 경우가 많다. 이에 본 연구에서는 수질모형의 매개변수 중 특히 수리특성에 관련된 매개변수들이 수질에 미치는 영향을 파악하는 것을 목적으로 하고 있다. 이를 위해 적용된 수질모형은 QualKo를 사용하였으며, 대상 하천은 낙동강 본류 경남구간 시점 부근인 회천 합류 전부터 낙동강 본류 경남구간 종점 부근인 밀양강 합류 전까지의 경남 오염총량관리 기본계획 시 구축된 모형 매개변수를 바탕으로 분석을 수행하였다. 일차오차분석을 이용하여 수리매개변수와 수질매개변수의 수질항목별 상대적 기여도를 파악해 본 결과, 수리매개변수는 DO, BOD, 유기질소, 유기인 모든 항목에 일정 정도의 상대적 기여도를 가지고 있는 것을 알 수 있었다. 이로부터 수질 모형의 적용 시 수리 매개변수 또한 수질 매개변수의 추정 시와 같이 보다 세심한 주의를 기울여 추정할 필요가 있을 것으로 판단된다.변화와 기흉 발생과의 인과관계를 확인하고 좀 더 구체화하기

  • PDF

Drought propagation assessment with WRF-Hydro model : from meteorological drought to hydrological drought (WRF-Hydro 모형을 활용한 가뭄전이 분석)

  • Lee, Jaehyeong;Kim, Yeonjoo
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2019.05a
    • /
    • pp.51-51
    • /
    • 2019
  • 기후변화 따라 과거에 경험하지 못했던 이상 수문기상 및 기상재해가 유발되며, 그 피해의 규모는 매년 증가한다. 그 중 가뭄은 미국 해양기상청(NOAA)이 선정한 20세기 최대 자연재해 중 상위 5위 안에 랭크되었으며 가뭄의 피해와 영향력은 막대하다고 언급하였다. 가뭄은 발생과정과 피해 영향에 따라 기상학적, 농업적, 수문학적, 사회경제학적 가뭄으로 구분할 수 있으며 직 간접적으로 영향을 미치고 기상학적 가뭄에서 수문학적 가뭄 등으로 가뭄의 종류가 변화되며 이를 가뭄전이라고 부른다,. 본 연구에서는 기상학적 가뭄에서 수문학적 가뭄으로 전이되는 과정을 분석하기 위하여 Weather Research and Forecasting and Model Hydrological modeling extension package (WRF-Hydro) 모형을 한반도 대상으로 구축하였다. 모형의 정확성을 향상시키기 위해 충주댐, 소양강댐, 용담댐, 남강댐의 유입량과 모형 유출량을 비교 분석하였으며, 유출에 영향을 미치는 지면 유출, 표면 거칠기와 같은 파라미터를 보정하여 주었다. 위와 같이 구축, 보정된 모형을 활용하여 모의된 유출량을 이용하여 수문학적 가뭄지수 Standardized Streamflow Drought Index(SSFI)를 도출하여 기상학적 가뭄지수 Standardized Precipitation Index(SPI)와 비교하여 기상학적 가뭄이 수문학적 가뭄으로 전이되는 과정을 가뭄의 빈도, 강도, 특성 등에 초점을 맞추어 분석하였다.

  • PDF

Seasonality Analysis of Soil Moisture using Cyclostationary Empirical Orthogonal Function (CSEOF 분석을 이용한 토양수분의 계절성 분석)

  • Cho, Eunsaem;Lee, Hyoungtaek;Lee, Myungseob;Lee, Youngju;Yoo, Chulsang
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2016.05a
    • /
    • pp.282-282
    • /
    • 2016
  • 지표수문해석모형이란 전 지구를 대상으로 수문해석 및 예측이 가능한 분포형 수문모형이다. 본 연구에서는 CSEOF(Cyclostationary Empirical Orthogonal Functions) 분석 방법을 이용하여 지표수문해석 모형 중 하나인 VIC(Variable Infiltration Capacity)모형의 토양수분 모의 성능을 평가해보고자 한다. 이를 위하여 먼저 남한에 대한 VIC 모형으로 모의한 토양수분 예측 결과와 관측자료를 수집하였다. 모의 성능 평가 기간은 1976년부터 2006년까지이다. 이후 본 연구에서는 수집된 VIC 모형의 예측 결과와 관측 자료에 대한 CSEOF 분석을 수행하여 각 자료의 월별 주된 변동 특성을 추출하였다. VIC 모형의 예측 결과와 관측자료의 상관관계는 CSEOF 분석 결과에 대한 Pattern Correlation으로 정량화되었다. 이와 더불어 본 연구에서는 모형의 모의 성능 평가에 주로 사용되는 NRMSE(Nomalized Root Mean Square Error)를 산정하여 예측 결과의 오차를 평가하였다. Pattern Correlation과 NRMSE를 모두 고려하여 VIC 모형의 성능을 평가해본 결과, 건기에 해당하는 기간과 우기에 해당하는 기간의 모의 성능이 다르게 나타났다. 본 연구의 결과는 추후에 지표수문해석 모형의 예측 결과를 이용하는 기후변화 관련 연구에 활용될 수 있을 것으로 판단된다.

  • PDF

Spatial Characteristics of Meteorological Data in Jeju Island (제주지역 기상자료의 공간특성 분석)

  • Lee, Jeong Eun;Lee, Jeongwoo;Kim, Chul Gyum;Chung, Il Moon
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2021.06a
    • /
    • pp.375-375
    • /
    • 2021
  • 제주지역의 경우에는 지형학적으로 한라산을 중심으로 한 산악지형을 이루고 있으며, 이러한 산악지형은 수문기상요소에 고도에 따른 영향을 미치게 된다. 따라서, 제주지역의 수문기상요소에 대한 산악효과를 정량적으로 파악하기 위해서는 해안 저지대에서부터 중산간, 고지대에 이르기까지 관측지점이 고르게 분포되어 있어야 한다. 제주 전역의 과거 운영되었거나, 현재 운영중인 기상관련 관측지점의 수는 125개에 이른다. 강수의 경우, 기상청(ASOS, AWS), 제주도 재난본부에서 관측하고 있으며, 강수 이외의 기상요소는 기상청에서 관측되어지고 있다. 본 연구에서는 제주전역의 기상관측지점에 대하여 강수량, 온도, 습도를 중심으로 제주지역의 산악효과 및 공간적인 특성을 분석하고자 하였다. 고도에 따른 기상요소의 정량적인 상관관계를 제시하고자 하였으며, 제주 전역에 걸친 수문기상요소들의 공간상관성 분석을 토대로 공간적인 분포특성을 도출하고자 하였다.

  • PDF