• 제목/요약/키워드: 유출계수

검색결과 779건 처리시간 0.038초

SWAT-K 모형을 이용한 연천지역의 지하수 개발가능량 추정 (Estimation of Groundwater Availability by Using the SWAT-K Model in Yeoncheon District, South Korea)

  • 이정은;김민규;정일문
    • 지질공학
    • /
    • 제34권2호
    • /
    • pp.263-277
    • /
    • 2024
  • 본 연구에서는 연천지역의 개발 가능한 지하수의 양을 추정하기 위해 분포형 수문모형 SWAT-K를 이용하여 토지이용, 토양분포 특성에 따른 함양량을 산정하였다. 모형의 검보정 결과 관측유량과 모의유량의 결정계수가 0.75~0.97로 양호하게 나타났으며 계산된 지하수 함양량은 겨울과 봄에 걸쳐 낮은 함양률을, 여름철에는 높은 함양률을 나타내는 시공간적 특성을 보였다. 계산된 연평균 함양률을 자연유량의 기저유출률과 비교하여 추정된 함양량의 타당성을 검토하였다. 지하수 개발가능량은 통계학적 빈도해석 기법을 이용한 10년 빈도 갈수시의 함양량으로 산정하여 연강수량의 11.5%인 것을 확인할 수 있었다.

Fe-loaded zeolite를 이용한 칼럼 실험에서의 Cd & Cr(VI) 동시제거 반응성 평가 (Simultaneous Removal of Cd & Cr(VI) by Fe-loaded Zeolite in Column System)

  • 이아라;이승학;박준범
    • 한국지하수토양환경학회지:지하수토양환경
    • /
    • 제11권1호
    • /
    • pp.14-22
    • /
    • 2006
  • Cd과 Cr(VI)으로 동시에 오염된 지하수 정화에 반응벽체 공법을 적용하기 위하여 환원능과 흡착능을 동시에 가지는 새로운 반응 물질인 Fe-loaded zeolite의 반응성을 실내 주상실험을 통해 평가하였다. Cd과 Cr(VI)을 동시에 포함한 오염용액을 반응물질이 충진된 칼럼에 주입하고, 이 칼럼을 통과한 유출수에서의 Cd과 Cr(VI) 농도를 측정하여 파과곡선을 얻는 방법으로 본 주상실험을 수행하였다. 유출수의 농도 분석 결과 Cd의 파과는 Cr(VI)의 파과곡선이 완료되는 시점까지 전혀 관찰되지 않아, 반응물질의 전반적인 효율은 Cr(VI)에 대한 반응성에 의해 결정됨을 알 수 있었다. Cr(VI)의 파과곡선은 파과 시작점이 다소 지체되고 파과의 진행 속도도 느려지며 파괴곡선의 상대농도 값이 1까지 증가하는 양상을 나타내었다. 이러한 Cr(VI) 파과곡선의 양상을 정량적으로 묘사하기 위하여 형상계수라고 하는 새로운 계수 ${\alpha}$${\beta}$를 정의하고, 이를 오염물질의 이동방정식에 적용함으로써 파과곡선의 지체 정도와 기울기 정도를 나타내었다. 초기 주입 용액의 농도가 증가할수록 ${\alpha}$값은 감소하였고, 이로부터 파과곡선의 시작 시점이 빨라 지는 것을 알 수 있었다. 또한 초기 주입 용액의 속도가 증가할수록 ${\beta}$ 값이 감소하였으며, 이로부터 파과곡선의 기울기가 점차 가팔라지는 것을 알 수 있었다. 이를 통해 Fe-loaded zeolite를 반응벽체 내 반응물질로 적용할 경우, Fe-loaded zeolite는 지하수 상에 존재하는 서로 다른 이온 형태의 중금속인 Cd과 Cr(VI)을 효과적으로 동시에 제거할 수 있음을 획인하였으며 이때, Fe-loaded zeolite를 통한 Cr(VI)의 파과는 지하수 흐름이 빠를수록 파과의 진행이 더욱 급격히 나타날 것이며 오염농도가 높은 지하수일수록 파과의 시점이 빠르게 나타나 반응물질은 그 사용 한계 에 다다르게 될 것임을 예측할 수 있었다.

제주도 지하수자원의 최적 개발가능량 선정에 관한 수리지질학적 연구

  • 한정상;김창길;김남종;한규상
    • 한국지하수토양환경학회:학술대회논문집
    • /
    • 한국지하수토양환경학회 1994년도 제주도 지하수자원
    • /
    • pp.184-215
    • /
    • 1994
  • 제주도 지하수자원의 산출특성을 규명키 위하여 총 455개 공의 자료를 전산처리하여 지역별 대수성 수리특성을 규명하였다. 제주도는 주로 현무암 내에 협재된 화산쇄설층, crinker층과 현무암의 1 및 2차 유효공극이 주 대수대의 역할을 하며 이들은 기저, 준기저및 상위대수층으로 구성되어 있다. 본도 대수충의 평균 투수량계수는 29,300m$^2$/일 이며 평균 저유계수는 0.12로써 자유면 대수층을 이루고 있다. 종합적인 물수지 분석을 실시한 바 본도에 부존된 지하수 부존량은 약 44억m$^3$이고, 년평균 강수량은 33.9억m$^3$으로써 이중 하천유출량은 6.38억m$^3$/년 이며, 증발산량은 12.56억m$^3$/년(37%)이고 지하수함양량은 년평균 강수량의 44.1%에 해당하는 14.94억m$^3$이다. 본도에 부존된 지하수의 최적 개발가능량(sustainable yield)을 각 지역별로 정량적으로 계산한 결과 그 양은 함양량의 41%에 해당하는 6.2억m$^3$/년(1,689,000 m$^3$/일)정도였으며 잔여 8.74억m$^3$/년(2,404,000m$^3$/일)은 해안이나 해저용천으로 유출된다. 특히 최근 심부 시추조사 자료에 의하면 EL-120$\pm$68m부근에 저투수성 해성 퇴적층(일명 세화리층)이 분포되어 있는 것으로 판명되었으며 과거 서귀포층군으로 알려진 저투수성 퇴적층이 북서부와 서부 일원에서 EL-70m 부근에 널리 분포되어 있어 서귀포층군자 세화리층의 명확한 구분이 필요하다. 만일 이러한 저투수성 퇴적층이 제주도의 기저층을 이루는 경우 제주도 내에 부존된 지하수는 주로 준기저 지하수일 것이며 이는 제주도 지하수의 산출특성에 결정적인 영향을 미칠 요인이다.rative processing at the best platform. Furthermore, from among the five structures utilized in Client/server architecture for distribution and cooperative processing of application between server and client this study presents two different data management methods under the Client/server environment; one is "Remote Data Management Method" which uses file server or database server and. the other is "Distributed Data Management Method" using distributed database management system. The result of this study leads to the conclusion that in the client/server environment although distributed application is assumed, the data could become centralized (in the case of file server or database server) or decentralized (in the case of distributed database system) and the data management method through a distributed database system where complete responsibility and powers with respect to control of data used by the user are given not only is it more adaptable to modern f

  • PDF

전북지역 농업용 하천유역의 수질과 부하량 특성 (Assessment of Water Quality and Pollutant Loads on Agricultural Watershed in Jeonbuk Province)

  • 엄미정;문영훈;안병구;신용규
    • 한국환경농학회지
    • /
    • 제27권2호
    • /
    • pp.111-119
    • /
    • 2008
  • 전북도내에서 농업용수원으로 사용되고 있는 소규모 하천들의 수질 및 오염부하량 조사를 통해 효율적인 수질개선을 위한 기초자료를 얻고자 본 시험을 수행하였다. 24개 하천의 수질은 EC가 $0.07{\sim}0.52$ dS/m, BOD가 $0.1{\sim}5.0$ mg/L, $COD_{Cr}$$0.6{\sim}17.7$ mg/L 등 이었으며, 범위는 다양하였으나 연도에 따른 평균값은 큰 차이를 보이지 않았다. 또한 ${SO_4}^{2-}$$5.7{\sim}52.7$ mg/L, $Cl^-$함량은 $0.9{\sim}57.5$ mg/L, $Na^+$함량은 $2.5{\sim}48.9$ mg/L의 범위를 보여 관개용수 수질로는 양호한 상태였다. 조사하천 대부분의T-N 함량은 5.0 mg/L 이하, T-P함량은 0.10 mg/L이하였으나 일부 하천에서 5.0 mg/L 이상의 질소를 함유하고 있는 경우가 있어 관개수질을 고려한 합리적인 시비가 필요한 것으로 나타났다. 강우기의 EC, $Ca^{2+},\;Mg^{2+},\;K^+$$Na^+$ 등의 함량은 감소하였지만, BOD, $COD_{Cr},\;COD_{Mn}$, T-N 및 T-P의 농도는 뚜렷한 차이를 보이지 않았으며, 높은 농도의 SS가 하천 유량의 증가와 함께 유출되고 있어 강우기의 지표면 유출이 수계에 영향을 주는 것으로 나타났다. 조사하천의 오염부하량은 BOD가 $9.6{\sim}757.9$ kg/day, T-N이 $51.2{\sim}1418.5$ kg/day, T-P가 $0.3{\sim}44.7$ kg/day로 유역 및 성분에 따라 많은 차이를 보이고 있었다. 갈수기에 비하여 강우기에 BOD는 $3.0{\sim}4.4$배, T-N는 $3.0{\sim}4.0$배, T-P는 $4.1{\sim}5.6$배의 부하량 증가를 보여 강우기에 주변 유역에서의 유기물을 비롯한 각종 영양염류의 유출이 많아지고 있었다. 논 면적비율이 30%이상인 하천유역(A그룹)이 30%이하인 하천(B그룹)유역이 비하여 영농기였던 7월의 강우기에 유량과 영양물질의 유달부하량이 상대적으로 적게 나타났으며, 갈수기 대비 강우기의 부하량 증가도 B그룹에서 유량이 5.7배 증가한 것에 비교하여 A그룹의 유량 증가는 2.3배로 적게 나타났다. BOD, T-N, T-P의 부하량과 유량과의 회귀식 결정계수($R^2$)는 BOD가 $0.8025^{**}$, T-N이 $0.9229^{**}$, T-P가 $0.7612^{**}$으로 유량과 T-N 부하량과의 상관이 가장 높게 나타났다.

충주댐 유역의 유출량에 대한 SWAT 모형의 예측 불확실성 분석 기법 비교 (Comparing Prediction Uncertainty Analysis Techniques of SWAT Simulated Streamflow Applied to Chungju Dam Watershed)

  • 조형경;박종윤;장철희;김성준
    • 한국수자원학회논문집
    • /
    • 제45권9호
    • /
    • pp.861-874
    • /
    • 2012
  • SWAT(Soil and Water Assessment Tool) 모형의 적용성 검증을 위해서는 매개변수 민감도 분석 및 검 보정, 예측 불확실성 분석을 필요로 한다. 최근 SWAT 모형의 불확실성을 분석하기 위한 다양한 기법들이 개발되었는데, 본 연구는 충주댐 유역(6,581.1 $km^2$)을대상으로유역출구점의 실측 일 유출량자료(1998~2003)를 바탕으로 SWAT 모형의 유출관련 매개변수에 대한 불확실성 분석을 실시하였다. 이때 사용된 분석 기법으로는 SUFI2 (Sequential Uncertainty FItting algorithm ver.2), GLUE(Generalized Likelihood Uncertainty Estimation), ParaSol (Parameter Solution) 등을 적용하였다. 이러한 기법은 모두 SWAT-CUP (SWAT-Calibration Uncertainty Program; Abbaspour et al., 2007) 모형에 탑재되어있으며, 모형의 결과로써 검 보정, 매개변수의 민감도 분석, 각종 목적 함수 및 불확실성의 범위 등이 자동으로 산출되므로 모형의 사용자가 불확실성 평가 기법의 분석 및 비교를 손쉽게 할 수 있다. 그 결과 대표적인 목적 함수인 결정 계수($R^2$; Legates and McCabe, 1999)와 NS (Nash and Sutcliffe, 1970) 모형효율은 모두 0.67에서 0.92 사이의 값을 나타내어 대체적으로 모의가 잘 이루어졌음을 알 수 있었다. 그러나 불확실성의 범위를 나타내는 지표인 p-factor 및 r-factor에서는 평가 기법 별로 그 차이가 확연하게 드러났다. 여기서 p-factor는 불확실성 범위에 실측치가 포함되는 비율이며, r-factor는 불확실성의 상대적인 범위로 각각 1과 0에 가까울수록 모의기법의 성능이 우수함을 의미한다. 세가지 알고리듬 중에서 SUFI2의 p-factor가 약 0.79로 가장 높게 나타났으며, ParaSol의 r-factor가 0.03으로 가장 작게 나타났다. 본 연구의 결과는 SWAT 모형을 이용한 수문 모의에서 수문분석에 따른 예측결과의 불확실성을 정량적으로 평가함으로서, 모형의 적용성 평가 및 모의결과의 신뢰성 확보에 근거자료로 활용이 가능할 것으로 판단된다.

선체 부식 및 손상 부재를 고려한 침몰선박의 구조 안전성 평가에 관한 연구 - 제7 해성호를 중심으로 - (Structural Safety Assessment of a Sunken Ship Considering Hull Corrosion and Damaged Members - Focus on the Sunken Ship 'No. 7 HaeSung' -)

  • 이승현;김원돈;서재준
    • 한국해양환경ㆍ에너지학회지
    • /
    • 제19권4호
    • /
    • pp.332-340
    • /
    • 2016
  • 침몰선박은 선체 내에 적재되어 있는 잔존연료 및 유해 화물 등의 유출로 인하여 주변 환경과 선박의 안전 항해에 피해 요인으로 작용할 수 있다. 대부분의 침몰선박은 해저에 침수된 상태로 비교적 안정적인 상태로 놓여있지만 선체 자중, 화물창의 잔존화물에 의한 하중과 해저면과의 접촉상태, 조류 및 해류 등의 영향으로 침몰선박 잔존연료 유출에 영향을 줄 수 있다. 침몰선박은 침몰 이후 부식이 진행되어 선체 구조부재의 두께가 지속적으로 감소하고 있는 상황이다. 따라서 침몰선박의 붕괴 가능성에 대한 구조안전성 평가가 필요하며 그 결과는 침몰선박의 위해도 평가요소 중 '유출가능성' 평가 항목에 반영되어야 한다. 본 연구에서는 침몰선박의 구조 안전성 평가 방법에 대한 단계별 절차를 제시하고, 집중관리대상 침몰선박으로 분류된 '제7해성호'를 대상으로 선체구조의 붕괴 가능성을 추정하기 위한 구조 안전성 평가 방안을 제시한다. 침몰 당시 선체 손상 상태와 침몰 이후 선체구조 부재의 부식량 추정을 바탕으로 대상 침몰선박의 선체구조 붕괴 가능성을 산정한 결과, 현 상태에서는 전체적인 선체구조 붕괴가 일어날 가능성이 낮은 편으로 추정되지만, 향후 시간에 따른 부식량 증가로 인한 선체구조 붕괴가 일어날 가능성은 다소 높은 것으로 판단된다. 따라서 부식 및 해양환경 변화가 침몰선체 구조 안전성에 미치는 영향 등의 연구는 지속적으로 수행 되어져야 할 것으로 사료된다.

유류 오염지역의 수산물 중 다환방향족탄화수소류 (PAHs) 분석 및 위해평가 (Analysis and risk assessment of polycyclic aromatic hydrocarbons (PAHs) in seafood from oil contaminated bay)

  • 정지윤;최찬웅;염태경;조경희;박세령;신호상;이광호;이효민
    • 분석과학
    • /
    • 제23권2호
    • /
    • pp.187-195
    • /
    • 2010
  • 허베이 스피리트호 유류유출 사고 후 지역주민들의 수산물 섭취로 인한 건강 위해도를 평가하기 위하여, 가스크로마토그래피-질량분석기(GC-MSD)를 이용하여 유류유출 지역 주민들이 직접 채취 및 구입한 수산물 중 다환방향족탄화수소류(PAHs) 16종을 분석하였다. 시료는 KOH로 가수분해하여 Methylene chloride로 추출하였다. 추출액 중의 16개의 PAHs는 실리카/플로리실 정제컬럼을 이용하여 Methylene chloride : n-헥산(1:9) 혼합용액으로 용출하였고 GC-MSD의 SIM을 이용하여 분석하였다. 16 개의 PAHs의 평균 회수율은 79~85%이었으며, 126개의 시료에서 $0.17\sim6.04\;{\mu}g/kg$ 수준으로 검출되었고, Benzo(a)pyrene의 독성등가계수(TEF)를 적용한 Benzo(a)pyrene 독성등가량($TEQ_{BaP}$)은 $0\sim0.91\;{\mu}gTEQ$/kg 수준이었다. 또한 유류유출 지역 주민들의 벤조피렌의 일일 평균 노출량은 $5.5{\times}10^{-8}\;mg/kg$ bw/day 이었고, PAHs 일일 평균 만성노출량은 $1.3{\times}10^{-5}\;mgTEQ$/kg bw/day이었다. 그리고 안전마진(MOE)과 초과발암위해도는 각각 $1.8{\times}10^6$, $9.8{\times}10^{-8}$으로 사람에게 위해영향발생 우려가 낮은 수준이었다.

Zone특성 분할을 통한 유형별 통행발생 모형개발 (Development of Trip Generation Type Models toward Traffic Zone Characteristics)

  • 김태호;노정현;김영일;오영택
    • 한국도로학회논문집
    • /
    • 제12권4호
    • /
    • pp.93-100
    • /
    • 2010
  • 통행발생은 4단계 모형의 처음 단계로 전체수요예측에 상당한 영향을 미치게 되므로 정확성이 무엇보다 필요한 단계라 할 수 있다. 현재 통행발생모형으로 도시교통 및 SOC시설 등의 계획에 널리 사용되고 있는 것은 선형회귀모형이며, 각종 사회경제지표와 통행발생량의 관계가 선형임을 전제로 한다. 하지만 급격한 도시개발이나 도시계획구조가 변경되었을 때 통행량을 추정하기 위한 사회경제지표 자료가 부족하여 추정된 통행량의 오차가 많을 수 있다. 이에 본 연구는 일반적으로 널리 사용되는 사회경제지표를 선형이란 가정을 하지 않고, 다양한 존의 특성을 반영할 수 있는 변수에 대한 시장분할을 토대로 새로운 유형별 통행발생모형을 개발하고자 한다. 본 연구에서는 교통수요예측의 처음 단계인 통행발생 모형의 예측력을 개선하기 위하여 존의 다양한 특성(토지이용, 사회경제적 등)을 고려하였다. 예측력 개선을 위한 시장분할 방법론으로는 통행 발생률을 기반으로 한 Data Mining(CART)방법과 회귀분석을 이용하였다. 연구의 결과를 살펴보면, 첫째, CART분석을 활용한 존 특성 분석결과, 유출통행은 사회경제적 요인(남녀상대비중, 연령대(22~29세))에 영향을 받고 있으며, 유입통행은 토지이용 요인(업무시설상대비중), 사회경제적 요인(3차 종사자상대비중)으로 나타났다. 둘째, 유형별 모형개발 결과 통행발생 계수 값은 유출의 경우 0.977~0.987(통행/인)이며, 유입의 경우 0.692~3.256(통행/인)로 나타나 유형구분이 필요한 것으로 나타났다. 셋째, 실측검증을 수행하였으며, 유출 및 유입의 경우 기존 모형보다 적합도가 높아진 것을 알 수 있다. 따라서 본 연구에서 개발한 유형별 통행발생모형이 기존 연구보다 우수한 것을 알 수 있었다.

수질정화용 인공습지 개방수역에서 유입수질에 따른 식물플랑크톤의 광합성특성 및 유기물생산력 (Photosynthetic Characteristics and Primary Production by Phytoplankton with Different Water Quality of Influent in Open Waters of Constructed Wetlands for Water Treatment)

  • 최광순;황길순;김동섭;김세원;김호준;조성주;박제철
    • 생태와환경
    • /
    • 제40권1호
    • /
    • pp.61-71
    • /
    • 2007
  • 본 연구는 국내 최대의 수질정화용 인공습지인 시화호 인공습지의 개방수역(open water)에서 식물플랑크톤에 의한 광합성특성 및 유기물생산력을 조사 평가하여 효율적인 습지관리를 위해 open water관리방안의 기초자료를 제공하기 위함이다. 식물플랑크톤의 광합성특성 및 1차 생산력은 유입수의 수질특성과 체류시간이 다른 반월천습지와 동화천습지의 open water에서 조사되었다. 조사기간 동안 시화호 인공습지내 5개 open water에서 식물플랑크톤에 의한 1차 생산력은 $481{\sim}11,275mgC\;m^{-2}\;day^{-1}$의 범위로 연중 부영양화 수준을 보였고 계절과 지점에 따라 변동이 컸다. 시화호 인공습지에서 식물플랑크톤의 P-I curve 모델계수를 통한 광합성특성은 두 습지로 유입되는 유입수의 상반되는 수질특성(인과 질소 농도)에 크게 영향을 받는 것으로 나타났다. 최대광합성속도 (Pmax)는 $42{\times}1,014 mgC\;m^{-3}\;hr^{-1}$로 chi. ${\alpha}$ 농도와 양의 상관(R=0.47)을 보였다. 식물플랑크톤의 광합성효율을 나타내는 동화계수(AN)두 습지 모두 고습지와 저습지의 AN값의 차이를 보이지 않았지만 동화천습지의 AN값이 평균 8.5gC $gChl^{-1}\;hr^{-1}$로 반월천습지의 평균 5.8gC $gChl^{-1}\;hr^{-1}$보다 높은 값을 보여 인 농도가 높은 동화천습지에서 식물플랑크톤의 광합성효율이 큰 것으로 나타났다. 두 하천의 수질특성으로부터 식물플랑크톤의 광합성효율은 질소보다 인 농도 그리고 작은 TN/TP와 관련이 있는 것으로 사료된다. 낮은 광도에 대한 광합성능력을 나타내는 초기기울기와 동화계수 사이에는 높은 양의 상관(R=0.81)을 보인 것으로 보아, 낮은 광조건에서 광합성능력이 큰 식물플랑크톤이 광합성효율도 큰 것으로 나타났다. 본 연구에서 시화호 인공습지의 open water에서 식물플랑크톤에 의한 유기물생산이 매우 높은 것으로 보아 식물플랑크톤이 습지의 수질 및 수처리효율에 영향을 미칠 것을 판단된다. 습지의 수처리효율을 향상시키기 위해서는 체류시간을 짧게 해주거나 open water의 면적을 줄여 식물플랑크톤의 증식을 억제하거나 또는 증식한 식물플랑크톤이 습지외부로 유출되지 않도록 하는 관리방안이 필요할 것으로 본다.

하구언 수문작동으로 인한 금강 하구역의 물리적 환경변화: II. 염분구조와 하구유형 (Physical Environment Changes in the Keum River Estuary by the Dyke Gate Operation: II. Salinity Structure and Estuary Type)

  • 이상호;권효근;최현용;양재삼;최진용
    • 한국해양학회지:바다
    • /
    • 제4권4호
    • /
    • pp.255-265
    • /
    • 1999
  • 금강하구언 수문작동이후 물성구조와 하구유형을 조사하기 위하여 1997년 6월, 7월 및 10월, 1998년 5월과 7윌에 CTD와 조류관측을 하였다. 조석과 조류의 관계는 고조전에 조류가 전류되어 낙조가 창조보다 1.5시간 정도 길다. 풍수기에는(5월부터 7월) 낙조시에 하구언에서 담수를 간헐적으로 대량 방류하여 하구 상류쪽의 수직적 염분차이와 염분의 시간변화가 커지고 하계 홍수시 2‰ 미만의 염분이 나타난다. 유출담수에 의한 염분약층의 발달로 2층구조를 보이며, 약층의 세기와 깊이는 저조시에 증가한다. 조석위상과 표층염분 변화의 관계는 히구 상류쪽에서 약 2.5 시간의 위상차이를 보이나 하구 하류쪽에서는 정상파 관계를 보인다. 하구 상류쪽에서의 위상차이는 낙조시 2~3시간 정도 방류된 담수에 의해 생성된 저염수가 조류에 의해 이동된 것을 의미한다. 갈수기에는 수직염분차이가 현저히 줄어든다. 관측된 표층유속과 염분구조 및 담수방류량을 사용하여 성층계수, 순환계수를 산출하였다. 금강하구는 부분혼합형 하구이고 풍수기와 갈수기에 따라 성층계수가 달라진다. 관측된 조류의 평균적 흐름은 상층에서 바다쪽으로, 저층에서 육지쪽으로 향하여 부분혼합형 하구의 전형적 순환형태와 일치하였다. 이러한 하구언 수문작동 후의 금강하구 유형과 순환은 하구역에서 저층부유사가 상류쪽으로 이동하며 퇴적되는 환경임을 제시한다.

  • PDF