• 제목/요약/키워드: 지점평균확률강우량

검색결과 45건 처리시간 0.034초

유전자 알고리즘을 이용한 강우강도식의 매개변수 추정에 관한 연구 (Parameter Estimation of Intensity-Duration-Frequency Curve using Genetic Algorithm)

  • 신주영;김수영;김태순;허준행
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2007년도 학술발표회 논문집
    • /
    • pp.142-146
    • /
    • 2007
  • 본 연구에서는 강우강도식의 매개변수를 보다 효율적으로 산정하기 위해서 유전자알고리즘을 적용한 매개변수 산정법을 제시하였으며, 지속기간의 장, 단기간에 따른 매개변수의 변화를 고려하기 위하여 다목적 유전자알고리즘을 적용하여 매개변수를 추정한 후 기존의 강우강도식에 의한 결과와 비교해 보았다. 매개변수는 지점빈도해석을 통해 산정된 확률강우량을 사용하여 추정하였고, 유전자 알고리즘의 목적함수로는 Nash & Sutcliffe Index, Root Mean Square Error(RMSE), Relative Root Mean Square Error(RRMSE), 결정계수, 평균들을 사용하여 가장 효율적인 형태의 목적함수를 구성하였다. 그 결과 기존의 매개변수 추정 방법들에 비해 유전자알고리즘을 이용한 경우에 더 정확한 강우량값을 산정할 수 있었고, 특히 다목적 유전자 알고리즘을 사용할 경우 장기간과 단기간에 걸쳐서 동시에 정확도를 향상시킬 수 있는 매개변수를 추정할 수 있는 것으로 나타났다.

  • PDF

BLRPM을 이용한 서울지점 시단위 강우의 월별 경년변화 연구 (A Study on The Monthly Trend of Seoul Hourly Rainfall Using BLRPM)

  • 황석환;김중훈;유철상;정성원;유도근
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2009년도 학술발표회 초록집
    • /
    • pp.1162-1166
    • /
    • 2009
  • 본 논문에서는 공간상의 한 점에 대한 연속시간 강우모델로 시간단위의 통계특성을 재현하는데 유용하다고 알려진 BLRPM(Bartlett-Lewis Rectangular Pulse Model)을 이용하여, 측우기 자료를 포함하는 서울지점의 일강우 자료를 분해하고 지속시간별로 집성하여 모의된 시단위 단기 강우의 정량적 장기변화 특성을 분석해 보았다. 모의자료의 전이확률과 발생빈도 분석결과, 측우기 관측계열(CWK)에 비해 근대 우량계 관측계열(MRG)의 무강우 지속기간은 길어지고 강우지속기간은 짧아졌음을 알 수 있다. 이는 과거에 비해 근래의 강우지속기간이 짧아졌음을 의미하고, 근래 강우량의 양적증가와 같이 고려해 보면 근래의 강우강도가 과거보다 매우 높아졌음을 증명해 준다. 1960년대 전후로 구분한 자료계열(BCC와 ACC)에 대하여 지속시간별 모의 강우량의 평균과 분산의 비(증감률)를 월별로 비교한 결과, 전반적으로 9월이 증가의 정도가 가장 컸고 8월도 비교적 큰 증가 경향을 보였다. 그러나 6월은 오히려 약간 감소한 경향을 나타냈다.

  • PDF

강우의 통계학적 다운스케일링을 위한 공간특성 분석 (Spatial Characteristics for Statistical Downscaling of Rainfall Data)

  • 이정은;이정우;김철겸;김남원
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2018년도 학술발표회
    • /
    • pp.166-166
    • /
    • 2018
  • 수자원 분야의 기후변화 연구에서 유출분석을 위한 장기유출모형의 입력자료로 일단위의 기상자료가 요구된다. 이러한 일자료의 생성을 위해 통계학적 다운스케일링 기법 중 추계학적 기상모의모형이 가장 널리 적용되고 있다. 또한, 유역단위의 합리적 유출분석을 위해서는 기상모의모형을 이용한 일자료 발생시 기상관측지점 간의 공간상관성 확보가 선행되어야 한다. 이러한 문제점을 극복하기 위한 다지점 추계학적 기상모형의 개발 및 적용에 앞서 기존모형의 강우 발생과 크기와 관련된 주요요소들의 공간적인 특성을 분석하고자 하였다. 따라서, 본 연구에서는 국내 기상청 지점의 관측자료를 중심으로 모형의 강우발생과 관련된 강우/무강우 발생확률, 강우크기와 관련된 월강우량의 평균값, 월평균 강우량의 표준편차, 왜곡도를 산정하였다. 이를 중심으로 전국에 걸친 공간특성 분석을 통하여 다지점 추계학적 기상모의모형의 개발 및 적용시 고려해야 될 사항을 도출하고자 하였다.

  • PDF

기후변화시나리오를 반영한 남북공유하천유역의 미래 수문특성 변화 전망 (Prospects of future changes of hydrological characteristics in South-North Korea river basin according to climate change scenarios)

  • 염웅선;박동혁;안재현
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2020년도 학술발표회
    • /
    • pp.266-266
    • /
    • 2020
  • 본 연구에서는 기후변화로 인한 남북공유하천유역의 미래 수문특성 변화를 전망하기 위해 ArcGIS 프로그램을 통해 산정된 격자형 수문특성 매개변수를 분포형 모형인 GRM에 적용하여 임진강유역의 미래 유출수문특성 변화를 분석하였다. 분포형 모형에 사용되는 강우량 자료는 기상관측소 단위로 상세화된 13개 전지구 기후 모델 중 RCP4.5, 8.5 시나리오의 공유하천유역 인접 11개 관측소별 빈도해석 결과를 시·공간적으로 분포하여 사용하였다. 또한 미래기간별 유출특성 변화추이를 분석하기 위하여 참조기간(1981-2005), 21세기 전반기(F1, 2011-2040), 중반기(F2, 2041-2070), 후반기(F3, 2071-2100)로 구분하여 분석을 실시하였다. 분석 결과 본 연구의 대상지점인 임진강유역은 기후변화로 인해 확률강우량이 증가하여 유역의 유출수문특성에 직접적인 영향이 있을 것으로 예측되었다. RCP 4.5 시나리오에서는 21세기 후반기인 F3에 확률강우량 및 유출량의 증가추세가 줄어들 것으로 전망되나, 참조기간 대비 F1에서 20.4%, F2에서 35.7%, F3에서 34.6%의 평균 유출량 증가율을 보였으며, RCP 8.5 시나리오에서는 F1에서 19.9%, F2에서 38.3%, F3에서 67.8%로 지속적인 증가가 전망되었다. 또한 첨두홍수량 발생시각은 참조기간 대비 약 4.6~13.3% 감소가 예상되었다. 기후변화로 인한 홍수량의 변화는 재해위험을 증가시킬 수 있으며, 이러한 상황에서 남한과 북한의 협력을 통한 유역통합관리의 필요성은 점차 커질 것으로 보인다. 이를 위해서는 정확한 수문학적 분석을 선행하여야 하며, 본 연구가 남북공유하천유역의 재해위험을 평가하는 기초자료로 활용될 수 있을 것으로 판단된다.

  • PDF

군집분석과 이변량 확률분포를 이용한 가뭄빈도해석 (Drought Frequency Analysis Using Cluster Analysis and Bivariate Probability Distribution)

  • 유지영;김태웅;김상단
    • 대한토목학회논문집
    • /
    • 제30권6B호
    • /
    • pp.599-606
    • /
    • 2010
  • 우리나라의 경우 짧은 강우 관측 기간으로 인해 지점빈도해석에 의한 가뭄분석은 불확실성이 크다는 문제가 있어 가뭄에 대한 지역빈도해석의 적용이 필요하다. 본 연구에서는 강우의 경향이나 크기를 고려하기 위해 단기가뭄을 나타내어 줄 수 있는 3개월 이동평균강우량과 장기가뭄을 나타내어 줄 수 있는 12개월 이동평균강우량을 산정한 후, 가뭄분석을 수행하기 위한 가뭄특성변수를 추출하였다. 가뭄특성변수를 이용하여 주성분분석과 군집분석을 수행하여 가뭄의 동질성을 갖는 관측지점들을 구분하였다. 또한, 본 연구에서는 지역별 가뭄빈도해석을 위해 이변량 확률분포함수를 적용하였으며, 가뭄 특성(가뭄 지속기간과 심도)의 상호 관계를 고려하여 지역적 가뭄특성을 종합적으로 판단하였다. 또한 이변량 핵밀도함수의 적용을 통해 가뭄 발생의 분포 및 경향성을 가장 근접하게 나타내어 줄 수 있는 결합 확률밀도함수를 추정하고, 군집지역별 2개월, 5개월, 10개월, 20개월의 가뭄지속기간과 5년, 10년, 20년, 50년, 100년의 재현기간에 따른 지역적 가뭄특성을 분석하였다. 그 결과 금강하류, 영산강의 일부 권역 및 남해안 일대에서 상대적으로 큰 가뭄심도가 발생하는 것으로 나타났다.

국내 강우의 이상치검정 방법의 비교 연구 (A Comparative Study on Methods for Outlier Test of Rainfall in Korea)

  • 이정식;신창동
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2018년도 학술발표회
    • /
    • pp.359-359
    • /
    • 2018
  • 이상치는 표본자료에서 크게 어긋나 다른 자료들로부터 떨어져 표시되는 자료로써, 실제로 발생할 확률이 매우 낮은 자료로 정의되고 있다. 설계홍수량을 산정하기 위하여 적용하고 있는 극치계열의 연최대치 강우자료에는 기계오작동 및 엔지니어의 표독오류가 발생하고 있으며, 기후변화에 따른 거대태풍 및 국지적인 집중호우 발생 등으로 인한 극치값 등에서 이상치가 관측되고 있다. 통상 이상치들은 통계분석시 자료 본연의 특성을 왜곡시켜 편향된 결과를 산정할 수 있으므로 빈도해석시 이상치해석 절차를 수행하여 자료의 적정성을 확인하여야 한다. 현재 실무에서는 설계홍수량 산정요령과 하천설계기준 해설 등에서 관련 내용을 기술하고 있지만, 국내 강우자료의 기록연수의 부족으로 인하여 빈도해석시 이상치 해석이 미수행되고 있어 이상치에 따른 자료편의가 발생하면 결과물인 확률강우량이 왜곡되게 산정될 수 있다. 따라서, 본 연구에서는 국내 주요 도시의 강우자료를 이용하여 이상치검정을 수행하였다. 대상지점으로는 서울, 부산, 대전, 대구, 인천, 광주, 울산 등의 비교적 긴 관측년수를 보유하고 있는 광역시를 선정하였으며, 지속기간은 10분, 1~24시간의 25개 강우자료를 적용하였다. 이상치검정 방법으로는 타 방법에 비하여 이상치 검정력이 뛰어난 것으로 알려진 2가지 방법을 채택하였으며, 표본자료의 평균과 표준편차로 표준화된 z값을 이용하여 상 하 한계선를 초과하는 값을 확인하는 z-Score 방법중 향상된 중위수 절대편차(MAD)에 의한 수정 z-Score 방법(Hoaglin, 1993)과 Box-Plot 방법(Tukey, 1969)을 적용하였다. Box-Plot 방법(Tukey, 1969)은 전체 자료를 25%씩 사분위로 구분하는 방법으로 정렬된 자료계열을 중앙값, 박스, 수염(whiskers), 이상치로 구분한다. 정렬된 25~75% 값들을 박스로 포함하여 외곽의 수염값들을 이상치로 분류하며, 특히 사분위수의 도식화로 데이터의 분포를 파악하기 좋으며, 이상치들의 위치와 자료의 비대칭 여부를 쉽게 파악할 수 있다. 본 연구의 수행으로 수정 z-Score 방법의 경우에는 서울과 대구지점에는 이상치가 없으며, 부산지점에는 13개, 대전지점 7개, 인천지점 5개, 광주지점 32개, 울산지점 26개가 나타났다. Box-Plot 방법으로는 서울지점 35개, 부산지점 39개, 대전지점 32개, 대구지점 38개, 인천지점 51개, 광주지점 61개, 울산지점 65개의 이상치가 분석되었다. 연구를 수행한 결과, 수정 z-Score 방법에 비하여 Box-Plot 방법에 의한 이상치가 더 많이 발생하였으며, 각각의 방법으로 지속기간 및 연도별 이상치 발생자료를 확인하였다. 방법별 이상치 발생현황 등을 분석하여 지점별 발생횟수를 분석하였으며, 추후 지점 및 자료의 보완이 수행되면 활용성을 증대시킬 수 있을 것으로 판단된다.

  • PDF

기후변화를 고려한 방재성능목표 강우량 개선 방향 (Improvement of Disaster Prevention Performance Target Rainfall Considering Climate Change)

  • 이정훈;김경민;김상단
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2018년도 학술발표회
    • /
    • pp.175-175
    • /
    • 2018
  • 우리나라에서 발생하는 대규모 자연재해의 상당부분은 강우에 의한 홍수피해이다. 최근 이러한 홍수피해는 기후변화와 더불어 극한강우 현상의 빈발에 의한 새로운 재해양상으로 전개되고 있으며, 이에 따라 정부에서도 재해발생시 원상복구의 개념이 아닌 항구복구의 개념으로 복구사업을 수행하고 있다. 그러나 설계에 기후변화에 대한 영향을 반영하고 있지 못하기 때문에 기후변화에 의하여 미래에 발생할 극한강우로 반복적인 피해가 예상되고 있으므로 기존의 방재성능목표 강우량의 설정 방법에 대한 개선이 필요하다. 전 세계적으로 이러한 기후변화에 의한 현상을 모의하기 위한 연구로 전지구기후모델(Global Climate Model, 이하 GCM)과 지역기후모델(Reginal Climate Model, 이하 RCM)을 사용하고 있다.우리나라 기상청에서도 CMIP5 국제사업의 표준 실험체계를 통해 전지구 기후변화 시나리오 산출을 위해서 영국 기상청 해들리센터의 GCM인 HadGEM2-AO를 도입하였다. 또한 한반도 기후변화 시나리오를 산출하기 위해 HadGEM3-RA 모형을 이용하여 전지구 기후변화 시나리오를 역학적으로 상세화하고 이를 한반도에 대해 12.5km 공간 해상도로 일 자료를 제공하고 있다. 하지만 유역규모 혹은 지점규모에서 사용하기 위해서는 이러한 일자료의 시 공간적인 상세화기법이 요구된다. 본 연구에서는 기후변화를 고려한 방재성능목표 강우량 개선 방향을 제안하기 위해 다양한 연구단에서 도출된 상세화 결과를 수집하고 비교분석을 통해 기후변화를 고려하고자 하였다. 다양한 연구기관에서 생산된 미래 확률 전망을 살펴본 결과, 동일한 GCM자료를 사용하더라도 상세화 방법론에 따라 서로 다른 결과가 도출되는 것을 확인하였다. 미래 예측의 불확실성을 고려하면 특정한 방법론이 우수하다고 평가하기는 어려움에 따라 앙상블 평균을 활용한 개선방향을 제안한다. 본 연구의 결과는 전국 지자체의 강우특성만을 고려한 것으로, 연안지역의 경우 해수면 상승을 고려하여 추가적인 대책이 필요할 것으로 판단된다.

  • PDF

청계천 유역의 옥상녹화 적용에 따른 유출저감 분석 (Analysis of Runoff Reduction applying Green Roof in the CheonggyeCheon Watershed)

  • 박구영;양정석;이재범;방용선
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2015년도 학술발표회
    • /
    • pp.400-400
    • /
    • 2015
  • 기후변화로 인해 예측하기 어려운 강우패턴의 변화와 도시화로 인한 불투수면의 증가로 도심지는 매년 침수위험에 노출되어 있다. 배수시설의 개선은 인구밀도가 높고 유동인구가 많은 도심지에서 이루어지기 쉽지 않은 상황이며, 강우패턴의 변화는 예측과 대처가 어려운 상황이다. 이러한 이유로, 침투증진을 통해 직접유출수를 줄이는 연구 중 LID(Low Impact Development)에 대한 관심이 높아지고 있다. LID기법은 도시화로 증가된 불투수면을 투수면으로 대체하여 저류, 침투, 여과, 증발산을 유도하여 물 순환 회복에 기여할 수 있는 방법으로 옥상녹화, 투수성포장, 침투 트렌치 등의 기술요소가 있다. LID 기술요소에 대해 국외의 경우 설치방안과 평균적인 저감효과를 메뉴얼로 제시하고 있지만, 강우 및 토지의 지역적 편차가 큰 국내에 적용하기에는 어려운 상황이다. 또한, LID 모델링의 경우 유역 내 적용 면적에 따른 유출저감효과를 제시하는 연구결과는 다수 제시되고 있지만, 적용 면적과 지점에 따라 효율성을 제시하는 점에서는 다소 미흡한 상황이다. 따라서 LID 기술요소 별 설계 사례와 GIS를 바탕으로 유역 내 적용 지점과 면적을 산출하고, 적용 지점 및 면적에 따른 저감 효율을 분석하는 것이 필요하다. 본 연구는 SWMM모형을 이용하여 LID 기술요소 중 옥상녹화를 적용하여, 강우 강도와 적용지점 및 면적에 따른 유출저감 효율을 분석하고자 한다. 연구지역은 청계천 유역으로 총 면적의 75% 이상이 주거지, 상업지, 교통시설 등을 포함한 불투수면으로 이루어져있는 도시 지역이다. 강우자료는 10년, 30년, 80년 빈도 확률강우량을 적용하였으며, 옥상녹화를 적용하기 위한 지점과 면적은 GIS를 기반으로 산정하였다.

  • PDF

금강유역(錦江流域) 주요지점(主要地点)의 제(諸) 수문량(水文量)에 관(關)한 연구(硏究)(I) (Studies on the Some Hydrological Quantities of Principal Locations in the Basin of Geum River(I))

  • 안병기;조성섭
    • 농업과학연구
    • /
    • 제2권1호
    • /
    • pp.281-300
    • /
    • 1975
  • 1. 금강유역(錦江流域)의 평균년강수량(平均年降水量)은 1203.0mm이고, 면적강수량(面積降水量)을 산출(算出)할 때 많이 사용(使用)하는 Thiessen법(法)의 Areal weight는 표(表)-과 같다. 2. 최대연강수량(最大年降水量)은 1501~2000mm인 곳이 24개관측소중(個觀測所中)며 17개소(個所)나 되어 가장 많은 71%이고, 나머지는 1500mm이하(以下)인 곳이 3개소(個所)이며, 2001mm이상(以上)인 곳은 4개소(個所)이다. 3. 최대일강우량(最大日降雨量)은 201~300mm인 곳이 15개소(個所)로 가장 많은 63%의 분포(分布)이다. 4. 최대(最大)2일간연속강우량(日間連續降雨量)은 300mm를 중심(中心)으로 전후(前後) 각각 50%씩 분포(分布)된다. 5. 최대(最大)3일간연속강우량(日間連續降雨量)은 301~400mm인 곳이 15개소(個所)로 63%의 분포(分布)이며 이는 최대일강우량(最大日降雨量)보다 100mm가 더 많은 강우량(降雨量)이다. 6. 최대연속강우량(最大連續降雨量)은 401~600mm인 곳이 14개소(個所)로 가장 많은 58%의 분포(分布)이며 이는 최대일강우량(最大日降雨量)보다 200mm가 더 많은 강우량(降雨量)이다. 7. 확률일강우량(確率日降雨量)은 표(表)-5와 같고, 무주(茂朱)를 중심(中心)으로한 대청(大淸)댐유역(流域)의 확률일강우량(確率日降雨量)이 가장 많으며, 다음으로는 강경(江景) 공주(公州) 및 부여(扶餘)를 연결(連結)하는 지방(地方)이 다우(多雨)인 편(便)이다. 8. 6월(月)1일(日)부터 9월(月)10일(日)까지 100일간(日間)의 논벼 관개기간중(灌漑期間中)의 강우량(降雨量)은 601~800mm 인 곳이 16개소(開所)로 76%의 분포(分布)이며 이는 표(表)-6과 같다. 9. 유효우량(有效雨量)은 501~600mm인 곳이 15개소(個所) 71%의 분포(分布)이며, 유효율(有效率)은 66~75%인 곳이 13개소(個所)로 62%의 비율(比率)이고, 다음으로는 76~85%인 곳이 7개소(個所)로 33%의 비율(比率)이다. 10. 확률유효우량(確率有效雨量)은 각 지방(地方)의 연평균유효우량(年平均有效雨量)은 보다 100mm정도(程度)씩 적고, 301~500mm인 곳이 14개소(個所)로 67%의 분포(分布)이다. 11. 갈수년도(渴水年度)는 유역내(流域內)를 총괄적(總括的) 동일(同一)하게 지정(指定)할수 없다고 여겨진다. 12. 금강유역(錦江流域)의 s/S는 0.53~0.74로 되어 수자원면(水資源面)에서 비교적(比較的) 안정(安定) 편(便)이다.

  • PDF

극치강수량의 시공간적 특성을 이용한 지역빈도분석 (Regionalization of Extreme Rainfall with Spatio-Temporal Pattern)

  • 이정주;권현한;김병식;윤석영
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2010년도 학술발표회
    • /
    • pp.1429-1433
    • /
    • 2010
  • 수공구조물의 설계, 수자원 관리계획의 수립, 재해영향 검토 등을 수행할 때, 재현기간에 따른 확률개념의 강우량, 홍수량, 저수량 등을 산정하여 사용하게 되며, 보통 대상지역의 장기 수문관측 자료를 이용하여 수문사상의 확률분포를 산정한 후 재현기간을 연장하여 원하는 설계빈도에 해당하는 양을 추정하게 된다. 미계측지역 또는 관측자료의 보유기간이 짧은 지역의 경우는 지역빈도 분석 결과를 이용하게 된다. 지역빈도해석을 위해서는 강우자료들의 동질성을 파악하는 것이 가장 기본적인 과정이 되며 이를 위해 통계학적인 범주화분석이 선행되어야 한다. 지점 빈도분석의 수문학적 동질성 판별을 위해 L-moment 방법, K-means 방법에 의한 군집분석 등이 주로 사용되며 관측소 위치좌표를 이용한 공간보간법을 적용하여 시각화하고 있다. 강수량은 시공간적으로 변하는 수문변량으로서 강수량의 시간적인 특성 또한 강수량의 특성을 정의하는데 매우 중요한 요소이다. 이러한 점에서 본 연구를 통해 강수지점의 공간적인 좌표 및 강수량의 양적인 범주화에 초점을 맞춘 기존 지역빈도분석의 범주화 과정에 덧붙여 시간적인 영향을 고려할 수 있는 요소들을 결정하고 이를 활용할 수 있는 범주화 과정을 제시하고자 한다. 즉, 극치강수량의 발생 시기에 대한 정량적인 분석이 가능한 순환통계기법을 이용하여 관측 지점별 시간 통계량을 산정하고, 이를 극치강수량과 결합하여 시 공간적인 특성자료를 생성한 후 이를 이용한 군집화 해석 모형을 개발하는데 연구의 목적이 있다. 분석 과정에 있어서 시간속성의 정량화 및 일반화는 순환통계기법을 사용하였으며, 극치강수량과 발생시점의 속성자료는 각각의 평균과 표준편차를 이용하였다. K-means 알고리즘을 이용해 결합자료를 군집화 하고, L-moment 방법으로 지역화 결과에 대한 검증을 수행하였다. 속성 결합 자료의 군집화 효과는 모의데이터 실험을 통해 확인하였으며, 우리 나라의 58개 기상관측소 자료를 이용하여 분석을 수행하였다. 예비해석 단계에서 100회의 군집분석을 통해 평균적인 centroid를 산정하고, 해당 값을 본 해석의 초기 centroid로 지정하여, 변동적인 클러스터링 경향을 안정화시켜 해석이 반복됨에 따라 군집화 결과가 달라지는 오류를 방지하였다. 또한 K-means 방법으로 계산된 군집별 공간거리 합의 크기에 따라 군집번호를 부여함으로써 군집의 번호순서대로 물리적인 연관성이 인접하도록 설정하였으며, 군집간의 경계선을 추출할 때 발생할 수 있는 오류를 방지하였다. 지역빈도분석 결과는 3차원 Spline 기법으로 도시하였다.

  • PDF