• Title/Summary/Keyword: 강우 침투 및 유출

Search Result 229, Processing Time 0.032 seconds

Watershed priority evaluation for water circulation management (물순환 관리를 위한 소권역 우선순위평가)

  • Kim, Seokhyeon;Kim, Sinae;Gwak, Jihye;Lee, Hyunji;Kim, Hakkwan;Kang, Moon-Seong
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2022.05a
    • /
    • pp.499-499
    • /
    • 2022
  • 급격한 도시화와 이로인한 불투수면의 증가는 유역의 물순환을 왜곡시키고 있다. 직접유출의 증가와 침투량 감소로 이어지며 이는 지하수 함유량의 감소와 하천건천화를 유발한다. 환경부는 더 이상의 물순환 왜곡을 막기위해 물환경보전법상의 물순환율을 정의하고 물순환 관리목표를 설정하였으며, 제3차 강우유출 비점오염원관리 종합대책에서 물순환 사업 추진을 위한 소권역 별 우선순위산정을 과제로써 제시하고 있다. 대표적인 물순환 사업은 저영향개발기법 설치이다. 저영향 개발기법은 토지 및 공간을 이용해야하기 때문에 설치비용이 높고 위치선정에 제약이 많으며 유지관리도 어렵다. 이에 물순환 사업에는 단순 불투수면 뿐만아니라 사회, 경제적 요소까지 고려하여야한다. 본 연구에서는 물순환 우선순위를 산정하기위해 불투수면뿐만 아니라 사회, 경제적 요소까지 고려한 종합적인 우선순위를 산정하였다. 유역 물순환 평가를 위해 PSR framework을 이용하였다. PSR framework는 OECD가 개발한 지속가능성 평가 개념이며, Pressure, State, Response 세 가지 요소로 구분해 평가하게된다. PSR framework의 기본 개념은 인간의 활동들이 환경에 압력 (P)를 주고, 이로 인해 자연의 질과 영향 (S)을 미치며, 이에대한 회복을 위해 인식과 행동을 통해 정책과 제도 등을 통해 반응 (R)한다는 것이다. 본 연구에서는 PSR framework의 평가요소 Pressure는 불투수면 및 강우, State는 물순환, 하수관거, 수질, 수생태계, Response는 협력, 인식, 재정, 토지로 구분하여 평가하였다. 최종적으로 불투수면적 감축 대상 소권역('30년 무대책 불투수면적률 25% 이상)소권역에 대하여 최종우선순위를 산정하였다.

  • PDF

Determination on the component arrangement of a hybrid rain garden system for effective stormwater runoff treatment (강우유출수 처리를 위한 하이브리드 빗물정원 시스템의 구성요소 배열 연구)

  • Flores, Precious Eureka D.;Geronimo, Franz Kevin F.;Alihan, Jawara Christian P.;Kim, Lee-Hyung
    • Journal of Wetlands Research
    • /
    • v.19 no.3
    • /
    • pp.271-278
    • /
    • 2017
  • Low impact development (LID) technology has been recently applied for the treatment of nonpoint source pollutants. Rain garden is one of the widely used LIDs since it utilizes various mechanisms such as biological and physico-chemical treatment to reduce pollutants. However, problem such as clogging has been one of the issues encountered by the rain garden that do not undergo constant maintenance. Therefore, this research was conducted to develop and determine the component arrangement of a rain garden system for a more efficient volume and pollutant reduction. Two hybrid rain garden systems having different characteristics were developed and evaluated to determine the optimum design and arrangement of the system. The results showed that the components arranged in a series manner showed a volume reduction of 93% and a pollutant reduction efficiency of approximately 99%, 93% and 95% was observed for particulates, nutrients and heavy metals, respectively. While when the system is connected in a combined series-parallel, the volume and average pollutant reduction efficiency for the TSS, nutrients and heavy metals are 65%, 94%, 80% and 85%, respectively. Moreover, the component arrangement in the order of sedimentation tank, infiltration tank and plant bed exhibited a high pollutant reduction efficiency compared when the infiltration tank and plant bed were interchanged. The findings of this research will help in the further development and optimization of rain garden systems.

Development of immersive flood routing model using three-dimensional virtual reality visualization (3차원 가상현실 시각화를 활용한 몰입형 홍수추적 모형 개발)

  • Son, Sangyoung;Hwang, Sooncheol
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2022.05a
    • /
    • pp.96-96
    • /
    • 2022
  • 가상현실 (Virtual Reality, VR) 기술은 3차원 가상공간 내에서의 높은 몰입감에 기반한 체험을 바탕으로 다양한 분야에서 활용되고 있다. 소방훈련이나 태풍, 지진 등 재해 대응훈련과 같이 인명피해의 위험이 있는 재해에 대한 VR 기술을 활용한 방재교육은 위험성을 동반하지 않으면서도 현장감에 기반한 높은 교육적 효과를 창출할 수 있다. 한국전자통신연구원에서는 VR 기술을 이용하여 소방훈련을 위한 실감 소방훈련 시뮬레이터를 개발한 바 있으며 목동재난체험관에서는 홍수, 태풍, 지진 등 다양한 재해에 대한 안전교육을 위한 자연재해 가상현실체험을 운영하고 있다. 이외에도 전국 지자체 및 교육청에서는 방재교육을 목적으로 VR 기술을 활용하고 있다. 그러나 기존의 VR을 활용한 수재해 방재교육은 범람의 수리학적 특성과 함께 수해지의 지형적 특성을 적절히 반영하지 못하는 단점을 가지고 있다. 이는 방재교육이나 경각심을 부각하는 데엔 효과적이나 실질적인 방재 가이드라인을 제시하는 데엔 한계가 있다. 본 연구는 몰입형 파랑해석모형인 Celeris Base를 기반으로 3차원 가상현실 시각화를 활용한 수리학적 홍수추적 모형을 개발하였다. 3차원 가상현실 시각화는 Unity3D를 이용하여 모의환경 내에 구현되었다. 강우-유출 과정의 수리학적 해석을 위해 동수역학 수치모형의 연속방정식 내에 강우와 침투에 대한 항을 추가하였다. 침투모형으로는 Horton 모형, Green-Ampt 모형과 함께 사면의 기울기를 고려한 Green-Ampt 모형을 적용하였다. 실제 유역에서의 홍수추적 모의결과는 관측값과 비교적 잘 일치함을 확인하였다. 개발된 모형은 VR 방재교육을 통해 일반인의 수재해 대응능력 향상에 기여함과 동시에 정확성 높은 홍수추적 모의결과에 기반한 홍수대책 마련에도 활용 가능할 것으로 기대된다.

  • PDF

Application of the Water Cycle Analysis Model for The Urban Chtahment Using WEP ( : Foucus in Multifuntional Administrative City of Before Developed) (WEP 모형을 이용한 도시 물순환 해석 (개발 전 행정중심 복합도시를 대상으로))

  • Lee, Yong-Jun;Jang, Cheol-Hee;Noh, Seon-Gjin;Kim, Hyeon-Jun
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2008.05a
    • /
    • pp.1264-1269
    • /
    • 2008
  • 최근에 이르러 농촌 인구의 도시유입으로 도시의 주거공간이 절대적으로 부족해지고 있으며, 이는 대도시의 아파트 가격 상승을 유발하는 등 부수적인 사회/경제적 문제를 일으키고 있다. 정부나 지방자치단체에서는 이러한 문제를 해결하기 위하여 대규모의 신도시 개발계획을 추진하고 있으며, 이러한 신도시 개발은 도시화로 인하여 환경적, 수문학적으로 영향을 미치는 것으로 조사되고 있으며 도시물순환계는 강우가 지표에 도달 한 후 바로 하천으로 유출되는 것과 침투하여 서서히 유출하는 자연경로 그리고 상수도 및 하수도를 통해 하천으로 배수되는 인공적 경로에 의해서 형성된다. 일반적으로 신도시 개발 전후의 수문순환을 평가하는 방법 중의 하나는 개발예정지에 대한 장단기의 수문 관측을 통하여 개발 전과 개발 후의 유출특성과 수질부하를 정량적으로 비교하는 것이다 (한국건설기술연구원, 2004). 따라서 본 연구에서는 행정중심 복합도시 물순환 해석을 하기 위해 대상유역을 100m 크기의 정방형 격자로 구분하고 기상조건, 지표면 조건, 하천, 토양, 지하 대수층, 농업용수 등의 광범위한 입력 자료를 구축한 후 이를 물리적인 기반의 공간분포모형인 WEP 모형에 적용하여 개발 전 홍수 및 유출특성을 규명하고자 한다. 이를 위해 대상지역에 대한 장기적인 수문모니터링 자료를 바탕으로 1996년에서 2007년 까지 모의 하여 실측자료가 있는 $2006{\sim}2007$년을 대상으로 비교한 결과 일부 갈수기 자료를 제외 하고, 전반적인 높은 상관성을 나타냈다. 하지만 일부 홍수위에서 금강으로 부터의 배수위 영향으로 인한 차이가 발생해 추가적인 검토가 필요하며, 본 연구에서 도출된 결과들은 향후 도시유역에서의 수자원 분배와 물자원 관리 등의 적응전략을 수립하는데 있어 본 연구가 도움이 될 것이라고 판단한다.

  • PDF

An Analysis of Soil Moisture Using Satellite Image and Neuro-Fuzzy Model (위성영상과 퍼지-신경회로망 모형을 이용한 토양수분 분석)

  • Yu, Myung-Su;Choi, Chang-Won;Yi, Jae-Eung
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2012.05a
    • /
    • pp.154-154
    • /
    • 2012
  • 지표에서의 토양수분은 작은 구성비를 가짐에도 불구하고 여러 수문 현상을 연계하는 매우 중요한 인자로써 최근 관련 연구가 활발하게 진행되고 있다. 토양수분은 침투나 침루를 통하여 강우와 지하수를 연결하는 기능을 함과 동시에 강우사상에 따른 유출특성에 직접적인 영향을 미치며 증발산을 통하여 에너지 순환을 연결하는 중요한 기능을 한다. 토양수분을 측정하는 방법에는 세타 탐침(Theta Probe), 장력계, TDR(Time Domain Reflectrometry) 등이 이용되고 있으며, 광역 토양수분자료의 보다 정확한 공간 변동성의 관측을 위하여 항공원격탐사와 인공위성 원격탐사기술이 개발되어 적용되고 있다. 인공위성 영상은 자료의 분석이 간편하며, 공간자료이므로 공간 변화를 분석하는 데 있어 매우 편리하다. 그 중 MODIS(Moderate Resolution Imaging Spectroradiometer) 위성영상은 저해상도 영상으로 극궤도 위성인 Terra와 Aqua 위성에 장착되어 있으며, NASA에서 필요한 정보를 받아 사용할 수 있다. 본 연구에서는 유역의 물리적 지형자료와 같은 방대한 양의 자료 수집 없이도, 모형이 구축되면 인공위성자료와 강우자료만으로도 신뢰성 높은 결과를 단시간 내에 효율적으로 산정할 수 있는 자료 지향형 모형인 ANFIS(Adaptive Neuro-Fuzzy Inference System)를 사용하였다. 사용된 퍼지변수로는 시험유역의 토양수분 관측자료와 강수량 및 인공위성 자료인 MODIS NDVI(Normalize Difference Vegetation Index), MODIS LST(Land-Surface Temperature) 영상을 이용하였다. MODIS NDVI는 시간 해상도 8일, 공간해상도 250 인 Level 3 영상이며, MODIS LST는 시간 해상도 1일, 공간해상도 1 km인 Level 3 영상을 사용하였다. 위성자료를 사용하기 위해 Korea TM 좌표체계로 변환한 뒤, 토양수분 관측지점이 속한 각 셀의 속성값을 추출하였다. 위성자료와 수집된 자료 및 토양수분자료와의 관계를 분석하기 위하여 입력자료를 다양한 방법으로 구성하여 입력 변수를 생성하였다. 생성된 입력 변수와 ANFIS 모형을 연계하여 각각의 토양수분 산정모형을 구축하고 대상지점에 대한 토양수분을 산정 및 비교 분석하였다.

  • PDF

Stability Analysis of Geocell Reinforced Slope During Rainfall (강우 시 지오셀 보강 사면의 안정성 평가에 관한 연구)

  • Shin, Eun-Chul;Kim, Jang-Ill
    • Journal of the Korean Geosynthetics Society
    • /
    • v.16 no.4
    • /
    • pp.33-41
    • /
    • 2017
  • In this study, the increment effect of safety factor according to increasing of horizontal permeability coefficient is analyzed when geocell is installed on the slope for protection. To evaluate the horizontal permeability and reinforcement effect, the laboratory tests such horizontal permeability test were conducted. According to the laboratory test results, as the porosity rate of geocell increases, the coefficient of horizontal permeability is also increased. And also, regardless of the different types of filled materials, the coefficient of horizontal permeability is improved in a geocell reinforced ground compare with the non-reinforced ground. Laboratory test results and the rainfall intensity were applied to the numerical modeling of slope for seepage analysis and stability analysis of slope by using Soilworks, numerical analysis program. As a result of the slope stability analysis, it is confirmed that the installed geocell on the slope facilitates the drainage of water on the surface of slope. Hence, the ground water elevation is suppressed. Therefore, the safety factor of the slope is increased by the increasing of the internal friction angle, apparent cohesion, and coefficient of horizontal permeability by reinforcing the slope with geocell.

Evalution of Input and Output Amount of Silica in Sandy Paddy Soil during Growing Periods of Rice Plant (벼 재배과정중 사질답에서 규소의 유, 출입량 평가)

  • Seo, Young-Jin;Kim, Jong-Su;Kim, Chang-Bae;Park, Man;Lee, Dong-Hoon;Choi, Choong-Lyeal;Choi, Jung
    • Korean Journal of Soil Science and Fertilizer
    • /
    • v.35 no.3
    • /
    • pp.162-168
    • /
    • 2002
  • The silica uptake of rice plant(Oryza sativa L.) and a potential capacity of paddy soil as a source of silica supply for rice plant was studied under field experimental condition. Amount of Si from irrigation water, rainfall, which was a source of silica supply, and percolating rate of water through paddy soil profile as a source of silica loss from paddy soil, was investigated. Rice plants were grown in three different conditions, which included Control as non-fertilized, RDA's recommandation(NPK, Urea-Fused phophate-Potassium chloride=$110-45-57kg\;ha^{-1}$) and RDA's recommandation with supplement of silicate fertilizer(NPK+Si, Urea-Fused phosphate-Potassium chloride+Si=$110\;N-45\;P_2O_5-57\;K_2O+267.2\;Si\;kg\;ha^{-1}$). An amount of silica supply from rainfall was $0.5kg\;ha^{-1}$ and average amount of silica supplied from irrigation water $42.5kg\;ha^{-1}$, ranging from $28.1kg\;ha^{-1}$ to $58.8kg\;ha^{-1}$. Silica amount percolated through the soil profile have uniform trend comparatively showing $62.9kg\;ha^{-1}$ in Control, $64.8kg\;ha^{-1}$ in NPK treatment and $62.9kg\;ha^{-1}$ in NPK+Si treatment. Silica uptaked by Rice plant was $335.6kg\;ha^{-1}$ in Control, $406.6kg\;ha^{-1}$ in NPK+Si treatment and $471.1kg\;ha^{-1}$ in NPK+Si treatment. The difference between an amount of Si input(from rainfall, irrigation water and silicate fertilizer) and an amount of Si output(percolated Si in soil, uptaked Si by rice plant) was $357.4kg\;ha^{-1}$ in control, $412.1kg\;ha^{-1}$ in NPK treatment and $238.2kg\;ha^{-1}$ in NPK+Si treatment. Results of our study imply that paddy soil is a potential pool as a source of Si supply during growing periods of rice plant.

The Water Circulation Improvement of Apartment Complex by applying LID Technologies - Focused on the Application of Infiltration Facilities - (LID 기술 적용을 통한 공동주택단지 물순환 개선 연구 - 침투시설 적용을 중심으로 -)

  • Suh, Joo-Hwan;Lee, In-Kyu
    • Journal of the Korean Institute of Landscape Architecture
    • /
    • v.41 no.5
    • /
    • pp.68-77
    • /
    • 2013
  • Change in the Global Climate causes flood, drought, heavy snow, and rainfall patterns in the Korean Peninsula. A variety of alternatives related to climate change are considered. The foreign researchers are interested in Low Impact Development(LID); the utilization of water resources and eco friendly development, over 10 years ago. The research and development of related technology has been advanced to apply LID techniques in order to develop several projects in the country. However, sharing of technology or system that can be used easily in the private sector is insufficient. The performance of the elements of LID Technology has not been fully agreed. LID elements of this technology are easy to apply to Apartment complex. The elements are classified technology. The infiltration of elements performs the functions of apartment complex landscaping space technology applied to the target. The water cycle simulation(SWMM 5.0) and technology the implementation of the effectiveness is also verified. For this purpose, three different places in apartment complex to target by SWMM5.0 U.S. EPA conducted utilizing simulated rainfall and applying LID techniques before and after the simulated water cycle (infiltration, surface evaporation, and surface runoff) were conducted. The importance of green space in the LID techniques of quantitative and qualitative storm water control as well as the role of Apartment Housing is to promote Amenity. Remember that the physical limitations of apartment complex and smooth water circulation system for the application of LID integrated management techniques should be applied. To this end, landscapes, architecture, civil engineering, environmental experts for technical consilience between the Low Impact Development efforts are required.

Application and Evaluation of improving techniques for watershed water cycle using downscaled climate prediction (상세화 기후전망자료를 활용한 유역 물순환 개선 기술 적용 및 평가)

  • Jang, Cheol Hee;Kim, Hyeon Jun;Cho, Jae Pil
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2019.05a
    • /
    • pp.334-334
    • /
    • 2019
  • 기후변화에 능동적으로 대처하기 위해서는 기후변화에 따른 수자원가용량의 변화를 정량적으로 평가할 수 있어야 한다. 평가결과의 신뢰도를 높이기 위해서 기후변화 시나리오는 지역기후 및 유역특성에 적합한 결과를 포함하여야 한다. 또한, 기후변화가 유역의 물순환계에 미치는 영향이 있다면, 물순환 개선 기술을 통해 지속가능한 유역 물환경을 구축하는 것이 필요하다. 유역 물순환 개선 기술은 기후변화가 진행 중에 있거나 예상되는 지역에 대하여 강우로부터 발생되는 유출을 지연, 저류, 침투시켜 지속가능한 물순환 체계를 유지하고 회복하도록 하는 기법이라 할 수 있다. 한국건설기술연구원에서는 기후변화에 따른 영향을 평가하고 적응 대책을 수립하기 위한 실무적인 유역 물순환 개선 및 평가 모형인 CAT3(Catchment hydrologic cycle Assessment Tool 3)을 개발하였으며 본 모형은 침투시설, 저류시설, 습지, 빗물저장시설과 같은 물순환 개선시설에 대한 효과를 정량적으로 평가할 수 있다. 본 연구에서는 팔당댐 상류의 경안천 유역을 대상으로 APCC 기후변화 시나리오 통계적 상세화 자료를 활용하여 물순환 개선 기술의 적용성을 평가하였다. 통계적 상세화 자료는 APCC에서 개발된 AIMS(APCC Integrated Modeling Solution) 플랫폼을 이용하였다. AIMS는 다양한 기후정보를 기반으로 사용자 관점에서 상세화를 수행할 수 있는 장점이 있다. 상세화 기법은 SDQDM(Spatial Disaggregation Quantile Delta Mapping) 방법을 이용하였다. 상세화된 기후자료는 과거자료의 재현성 및 미래 기간에 대한 왜곡도를 평가하기 위해 극한기후지수(Climate Index)를 이용하는데 본 연구에서는 장기간에 걸친 수자원가용량의 평가 및 예측을 위해 연강수량(PRCPTOT)을 사용하였으며 증발산량의 평가 및 예측에 영향을 미치는 온도 관련 극한기후지수는 평균기온 개념의 DTR(TMAX&TMIN)을 이용하였다. 통계적 상세화 과정을 통해 최종적으로 HadGEM2-CC, INMCM4, CanESM2 시나리오를 선택하였으며 각 시나리오별 물순환 개선 기술을 적용한 후 미래의 수문학적 변동성을 평가하였다.

  • PDF

Comparison of Discharge Characteristics of NPS Pollutant Loads from Urban, Agricultural and Forestry Watersheds (유역으로부터 배출되는 비점원 오염부하의 토지이용에 따른 특성 비교)

  • Yur, Joong-Hyun;Kim, Geon-Ha;Kim, Jeong-Kon
    • Proceedings of the Korea Water Resources Association Conference
    • /
    • 2005.05b
    • /
    • pp.1214-1218
    • /
    • 2005
  • 지금까지 국내에서의 수질개선을 위한 노력은 점오염원에 대한 저감만을 중심으로 진행되어 왔기 때문에, 보다 획기적인 수질개선을 위해서는 비점오염에 대한 연구와 관리가 진행되어야 한다. 그러므로 본 연구에서는 도시지역, 농촌지역, 임야지역의 토지이용특성이 다른 3유역을 대상으로 현장연구를 실시하여 비점오염물질의 발생특성 및 유역별 비교를 실시하였다. 측정은 도시지역의 가장 큰 오염원인 CSOs(Combined Sewer Overflows)에 대하여 측정을 실시하였고, 농촌지역 및 임야지역의 경우 각 유역의 출구 지점에서 측정을 실시하였다. 강우특성이 다른 15개의 강우사상을 대상으로 유량 및 SS, TCOD, TN, TP의 항목에 대하여 한 강우사상당 $15\~20$회의 측정을 실시하였다. 각각의 강우사상에 대하여 EMCs(Event Mean Concentrations)를 산출하여 도시지역, 농촌지역, 임야지역의 각 유역에 대한 확률별 EMCs를 산정한 결과 3개 유역의 EMCs는 도시지역>농촌지역>임야지역 순으로 나타났다. $EMC_{TCOD}$는 도시지역과 농촌임야지역간에 가장 큰 차이를 보이는 것으로 나타났으며, $EMC_{TN}$에서 가장 작은 차이가 나타났다. 각 유역별 EMCs의 로그-정규 확률그래프의 분산계수를 비교한 결과 농촌임야지역은 도시지역에 비하여 오염물질의 농도 변화가 강우특성에 따라 보다 큰 변화를 보이는 것으로 나타났다. 연구유역에서 $EMC_{TN}$의 발생확률 $50\% 값이 도시지역은 17.0mg/L, 농촌임야지역은 4.5mg/L로 나타났으며, 이는 유사한 유역특성을 나타내는 타 지역에서의 연구자료를 분석 값과 매우 유사한 크기를 가지는 것으로 나타났다.를 분석하였다. 실험을 수행하여 보다 정밀한 공식으로 개선할 수 있었다.$10,924m^3/s$ 및 $10,075m^3/s$로서 실험 I의 $2,757m^3/s$에 비해 통수능이 많이 개선되었음을 알 수 있다.함을 알 수 있다. 상수관로 설계 기준에서는 관로내 수압을 $1.5\~4.0kg/cm^2$으로 나타내고 있는데 $6kg/cm^2$보다 과수압을 나타내는 경우가 $100\%$로 밸브를 개방하였을 때보다 $60\%,\;80\%$ 개방하였을 때가 더 빈번히 발생하고 있으므로 대상지역의 밸브 개폐는 $100\%$ 개방하는 것이 선계기준에 적합한 것으로 나타났다. 밸브 개폐에 따른 수압 변화를 모의한 결과 밸브 개폐도를 적절히 유지하여 필요수량의 확보 및 누수방지대책에 활용할 수 있을 것으로 판단된다.8R(mm)(r^2=0.84)$로 지수적으로 증가하는 경향을 나타내었다. 유거수량은 토성별로 양토를 1.0으로 기준할 때 사양토가 0.86으로 가장 작았고, 식양토 1.09, 식토 1.15로 평가되어 침투수에 비해 토성별 차이가 크게 나타났다. 이는 토성이 세립질일 수록 유거수의 저항이 작기 때문으로 생각된다. 경사에 따라서는 경사도가 증가할수록 증가하였으며 $10\% 경사일 때를 기준으로 $Ro(mm)=Ro_{10}{\times}0.797{\times}e^{-0.021s(\%)}$로 나타났다.천성 승모판 폐쇄 부전등을 초래하는 심각한 선천성 심질환이다. 그러나 진단 즉시 직접 좌관상동맥-대동맥 이식술로 수술적 교정을 해줌으로써 좋은 성적을 기대할 수 있음을 보여주었다.특히 교사들이 중요하게 인식하는 해

  • PDF