• 제목/요약/키워드: 연속강우

검색결과 263건 처리시간 0.029초

가뭄 발생에 따른 논에서의 수문학적 거동 변화 분석 (Analysis of the Hydrological Responses Changes of Paddy Fields in Drought Season)

  • 박상현;김현준;장철희
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2020년도 학술발표회
    • /
    • pp.254-254
    • /
    • 2020
  • 우리나라의 주요 식량은 쌀이기 때문에 논 경작이 매우 활발한 국가이며 국토전체면적의 약 8.27%가 논으로 구성되어 있다. 그렇기 때문에 우리나라 유역을 대상으로 수문 모델을 모의함에 있어서 유역 내 논에서의 수문학적 반응을 살펴보는 것은 정확한 유출분석을 위해 반드시 필요하다. 본 연구에서는 물리적 매개변수 기반의 물 순환 해석 모형인 CAT(Catchment Hydrologic cycle Assessment Tool)을 이용하여 충청남도 보령시에 위치한 보령댐 유역을 대상으로 일 자료기반의 유출모의를 수행하였다. CAT은 논에서의 담수심과 물꼬높이 등을 고려한 모듈을 포함하고 있기 때문에 우리나라 논에서의 유출특성을 반영한 유역 유출 모의가 가능하다. 연구기간은 2000년부터 2017년으로 총 18년의 강우, 유출 및 기상자료를 수집하여 전체 기간에 대해 모의하였다. 전체 기간 모의 결과에 대한 유출률과 토양수분을 분석하여 연속된 평년기간(2003년 - 2005년)과 연속된 가뭄년 기간(2015년 - 2017년)으로 시나리오를 구분하였으며, 각 3년 기간으로 구분된 두 가지 시나리오에 대한 논에서의 수문학적 거동 변화를 토양수분과 지하수수위를 중심으로 비교 및 분석하였다.

  • PDF

강우-유출 예측모형 개발을 위한 자기조직화 이론의 적용 (Application of Self-Organizing Map Theory for the Development of Rainfall-Runoff Prediction Model)

  • 박성천;진영훈;김용구
    • 대한토목학회논문집
    • /
    • 제26권4B호
    • /
    • pp.389-398
    • /
    • 2006
  • 본 연구에서는 강우의 시 공간적 분포의 불규칙한 변동성을 고려한 강우-유출예측모형을 위해 인공신경망(Artificial Neural Networks: ANNs)의 기법의 일종인 자기조직화(Self Organizing Map: SOM) 이론과 역전파 학습 알고리즘(Back Propagation Algorithm: BPA을 복합적으로 이용하였다. 기존의 인공신경망 연구에서 야기된 저 갈수기의 유출량에 대한 과대평가, 홍수기의 유출량에 대한 과소평가, 예측값이 연속적으로 선행 유출량을 나타내는 Persistence 현상을 해결하기 위하여 패턴분류 성능을 지닌 SOM 이론을 예측모형의 전처리 과정으로 이용하였다. 먼저, 본 연구에서 제안한 방법은 SOM에 의해 강우-유출 관계를 분류하고, SOM에 의한 분류에 따라 각각의 모형을 구성한다. 개별적으로 구축된 모형은 유출량의 예측을 위해 각각의 양상에 따라 분류된 자료를 이용한다. 결과적으로 본 연구에서 제안한 방법은 과거의 인공신경망의 일반적인 적용에 의한 결과보다 더 나은 예측능력을 보여주었으며, 더불어 유출량의 과소 및 과대추정과 Persistence 현상과 같은 문제점이 나타나지 않았다.

정량강수모의를 이용한 실시간 유출예측 (Realtime Streamflow Prediction using Quantitative Precipitation Model Output)

  • 강부식;문수진
    • 대한토목학회논문집
    • /
    • 제30권6B호
    • /
    • pp.579-587
    • /
    • 2010
  • 기상청에서 제공하는 강우수치예보정보를 활용하여 10일이내의 중기유량예측을 수행하였다. 기상청의 원시예보자료로는 2일예보를 위한 RDAPS와 10일예측을 위한 GDAPS예측자료를 활용하였다. 수치예보의 정확도를 제고하기 위하여 강우상세 정보를 생산할 수 있는 강수진단모형(QPM)과 QPM모의결과에 내재된 계통적 편이를 제거하기 위하여 분위사상과정 (Quantile Mapping)을 적용하였다. QPM모의결과를 유출모형의 입력정보로 활용하기 위하여 일관적인 체계를 갖춘 유역강수 정보로 변환하여, 장기연속유출모형인 SSARR모형을 이용하여 금강유역내 주요지점에서의 유량예측을 수행하여 유량예측에 대한 검증을 수행하였다. 2006년 1월 1일부터 6월 20일까지 강수예측을 수행한 결과 2일예측인 RQPM의 경우 기간 총강수량을 기준으로 실적강우대비 89.7%의 강수모의값을 보임으로서 양호한 예측성능을 확인할 수 있었다. 유량예측모의에 있어서는 2일예측의 경우 일부 강우사상에서 예측누락과 예측오류가 발생하였지만 전반적으로 유량예측이 양호한 수준이었다. 다만, 하류지점의 경우 조절유량에 의한 유출모형보정의 어려움과 수위-유량관계곡선의 신뢰도저하등의 이유로 예측성능이 떨어지는 경우도 있었다. GQPM에 대한 10일강우예측은 첨두강수와 강수총량에 있어서 다소 과소한 모의값을 보이고 있으며, 강수보정효과도 RDAPS에 비하여 저조한 수준이었다. 이 부분은 강수예측의 사후보정으로는 한계가 있는 것으로 보여지며 원시예측모형의 안정화를 통하여 개선할 수 있는 부분으로 판단된다.

2023년 장마기간 동안 산악지역의 강우 특성 (Precipitation Characteristics in Mountainous Regions During Changma Period in 2023)

  • 김인혜;장근창;유병오
    • 한국농림기상학회지
    • /
    • 제26권3호
    • /
    • pp.161-173
    • /
    • 2024
  • 우리나라는 여름 장마철에 연 강우량의 50~65%가 집중되는 특징이 있다. 2023년 장마기간에는 전국적으로 660.2mm의 강우가 발생하였으며, 남부지방에서는 712.3mm가 기록되어 관측 이래로 가장 많은 비가 내렸다. 산림지역에서의 집중호우는 산사태를 유발하는 주요 원인이다. 따라서 장마철 산사태 재난 예방을 위해서는 산림지역에서 발생하는 강우를 정확하게 모니터링하고, 그 특성을 파악하는 것이 중요하다. 따라서 본 연구에서는 2023년 장마철 산악지역의 강우 특성을 살펴보기 위해 산림청 산악기상관측소(AMOS, Automatic Mountain Meteorology Observation Stations)의 자료를 활용하였으며, 비산악지역과의 강우량 비교를 위해 기상청의 강우량 관측 자료를 함께 활용하였다. 또한, 산사태 피해가 발생했던 경상북도 지역의 강우 분포를 확인하기 위해 MK-PRISM (Modified Korean-Parameter elevation Regressions an Independent Slopes Model)을 사용하여 강우량 분포도를 생산하였다. 강우량 관측 자료를 분석한 결과, 장마기간 중 강우량은 대체로 고도가 높은 산악지역에서 더 많은 것으로 나타났으며, 산사태가 발생했던 경북의 문경시, 봉화군, 영주시에서는 산악지역과 비산악지역 간의 강우량 차이가 최대 270mm였다. 특히, 이 지역에서는 산사태가 발생하기까지 50시간 이상의 연속 강우와 시간 당 30mm/hr 이상의 집중호우, 300mm 이상의 누적 강우가 발생했다. 강우분포도를 이용하여 산사태 발생 지역의 강우량을 비교한 결과, 장마기간 동안의 누적강우량은 평균 기후값 대비 2배 이상 많은 것으로 나타났다. 향후 효과적인 산사태 재난 대비를 위해서는 산림지역에서의 정확한 강우 관측과 함께 최근의 강우 특성을 고려한 기준 마련이 필요할 것으로 사료된다.

하천수의 생태독성을 파악하기 위한 황산화미생물의 이용가능성 평가 (Assessment of Biological Toxicity Monitoring in Water Using Sulfur Oxidizing Bacteria)

  • 강우창;오상은
    • 한국환경농학회지
    • /
    • 제31권2호
    • /
    • pp.170-174
    • /
    • 2012
  • 일반적으로 실제 하천의 EC 값은 일정하게 유지가 된다. 연속식, 반연속식 운전으로 실험한 결과 정상상태(steady state)를 보였다. 연속운전의 경우 유출수의 EC 값은 1.7 dS/m로 일정하게 유지되었으며(유입 EC = 0.1 dS/m) 반연속식운전의 경우 새로운 유입수가 유입됨으로 인해 초기 EC는 약 0.1 dS/m이고 미생물이 황을 산화함에 따라 황산염 이온이 생성되고 EC가 증가하여 30분 후 EC는 약 0.6 dS/m로 측정되었고 장시간 운전하여도 일정한 트랜드의 안정적인 값이 측정되었다. ${NO_2}^{-}-N$가 주입되면서 미생물이 저해를 받아 EC가 감소하는 경향을 보였으며 이를 통해 독성물질 주입여부를 판별할 수 있었다. 독성이 없는 유입수를 다시 유입시켰을 경우 약 4~5시간 후 반응조 내의 미생물이 원상태를 회복하여 유입 전 EC의 값이 다시 측정되었다. 황을 이용한 독성모니터링방법은 다른 여러 조건의 실험을 통해 문제점 파악 및 보완이 필요하지만 연속식 운전, 반연속식 운전 2가지 방법을 통해 독성물질을 신속하게 탐지할 수 있을 것으로 판단된다.

한강수계 하천에서의 시공간적 수질변화 특성 및 연속적 인공댐호의 경험적 모델 (Spatio-temporal Water Quality Variations at Various Streams of Han-River Watershed and Empirical Models of Serial Impoundment Reservoirs)

  • 전혜원;최지웅;안광국
    • 생태와환경
    • /
    • 제45권4호
    • /
    • pp.378-391
    • /
    • 2012
  • 본 연구에서는 2009~2010년 동안 한강 수계 인공댐 호소수 8개 지점 및 한강 본류 수계의 하천수 10개 지점의 총 18개 지점에서 인공댐 및 하류역 하천들에 대한 수질변이특성을 평가하고자 하였다. 이미 건설된 연속적 인공댐들이 상류로부터 하류까지 질소 (N), 인(P), 엽록소-a(CHL-a) 등의 화학적 수질 구배특성(Chemical gradient)에 미치는 영향을 분석하고, 수질 변수간의 경험적 모델(Empirical model)을 적용하여 하절기 몬순강우 특성 및 연별 집중강우 영향을 분석하였다. 한강수계는 상류에서 하류로 갈수록 수질이 악화되는 경향을 보였고 특히, 점오염원에 의한 오염물질의 유입은 수질변이에 가장 크게 작용하였으며, 특히 하류부의 중랑천 유입수는 수질을 급격하게 악화시키는 것으로 나타났다. 또한 하절기의 몬순강우는 하류역에서 총인(TP), 총질소(TN) 및 전기전도도(EC)값을 크게 낮추어 수질향상에 기여하였다. 5개의 인공댐들에서 총질소(TN), 총인(TP) 및 N:P 무게비에 대한 엽록소-a (CHL-a)의 경험적 모델식 평가에 따르면, 로그전환된 총인(TP)은 엽록소-a (CHL-a)농도 변이를 33.8% ($R^2$=0.338, p<0.001, 회귀식 기울기=0.710) 설명하였으나, 총질소 (TN) 변이는 21.4% 설명에 그쳤다 ($R^2$=0.214, p<0.001). 또한 N:P 무게비의 산정에 따르면, 본 인공호들에서는 모두 N:P 비가 29 이상으로 인 제한효과로 나타났으며, 총질소(TN)의 농도는 모든 수체에서 이미 $1000{\mu}g\;L^{-1}$을 상회하여 조류성장에 풍부한 것으로 나타나 엽록소-a (CHL-a)의 증감은 인(P)에 의해 조절되는 인 제한요인(P-limitation)으로 나타났다. 그러나 엽록소-a (CHLa)의 증감은 또한 개별적 인공댐에서 보여주는 바와 같이 부유물 증가에 의한 광투과도 저하, 수체류시간 감소에 의한 세척효과(Washing-out) 등도 부가적으로 영향을 미치는 것으로 나타났다. 결론적으로, 한강수계에서 하천 및 인공댐의 연속선 상에서 하천에 대한 댐의 이화학적 수질 특성에 미치는 영향은 미미한 것으로 나타났다. 반면, 유입 지천과 하천 유역의 점오염원으로부터의 방류수 유입은 이화학적 수질 변이에 크게 영향을 미치는 것으로 나타났다. 향후 연속적 인공댐의 수체에서 대형 무척추동물 및 어류와 같은 생물학적 지표를 이용한 영향평가도 수질과 같은 측면에서 연구될 필요가 잇는 것으로 사료되었다.

연속유출분석모형을 이용한 농업용 저수지 물 수지 분석 (Water Budget Analysis of Agricultural Reservoirs Using Continous Rainfall-runoff Model)

  • 이명진;정우창
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2011년도 학술발표회
    • /
    • pp.293-293
    • /
    • 2011
  • 고성군의 상수도 보급률은 2008년 기준으로 61%이며, 인근지역인 통영시의 92%에 비하여 매우 낮은 실정이다. 고성군 용수공급 전망을 살펴보면 2014년 일 최대 $25,830m^3$/일의 수요가 예상되나 현재 광역상수도를 전량 수수하는 고성군의 배분량은 $20,500m^3$/일로 $5,330m^3$/일이 부족하여 용수 확보가 필요한 실정이라고 볼 수 있다. 이러한 배경을 바탕으로 본 연구에서는 고성군의 218개소 저수지중 주변의 오염요소가 적은 42개소 저수지 유역에 대해서 강우-유출모형을 모의 수행하였으며, 토양과 관련된 각종 자료와 이를 GIS기법을 이용하여 매개변수를 산정하여, 연속유출분석모형인(NWS-PC)를 통하여 42개 유역의 유입량을 모의 수행하였다. 이러한 산정된 매개변수를 통해 2006년~2009년 동안의 강우량을 토대로 유출모의를 수행하였으며, 2008년 강우량을 토대로 극한 가뭄시의 경우도 모의 수행하였다. 최악의 가뭄시 최저 강우량을 2008년으로 선정하였기에 고성군의 평균유입량은 2,866.8천$m^3$/년으로 나타났다. 1990년~2005년도 평균 농업용수 이용량은 5,879.4천$m^3$/년이며, 2006년~2007년의 평균 농업용수 이용량은 2,186.2천$m^3$/년으로 2005년 이전 농경지 면적이 50%이상 줄어듦에 따라 이용량도 줄었다. 고성군의 평균 식수량은 406.3천$m^3$/년으로 산정되었으며, 고성읍, 하이면, 마암면, 거류면에서 유입량과 농업용수 및 식수량 차이에서 부족한 현상이 나타났으나 고성군 전체 평균으로 봤을 때는 274.3천$m^3$/년의 양이 남는 것으로 나타났다. 본 연구를 통하여 고성군 지역은 대부분의 농촌지역에서 생활용수를 마을상수도나 지하수에 의존하며 생활한다. 고성군 관내 농업용 저수지 중 식수겸용 저수지로 활용함으로써 인근 마을의 생활용수를 공급하게 되면 안정적인 마을상수원확보와 이상기후에 대비한 물 부족 대책을 수립할 수 있으며, 상습가뭄지역의 물 부족 문제가 해결하는데 기초자료로 활용할 수 있다고 판단된다.

  • PDF

하굿둑 실증시험 시 염지하수의 연직 프로파일링 관측 (During Demonstration Test of the Estuary Bank, Vertical Profiling Measurement in a saline groundwater)

  • 김병우;최두형
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2022년도 학술발표회
    • /
    • pp.447-447
    • /
    • 2022
  • 하굿둑 실증시험 및 개방 시 하굿둑으로부터 해안지역 대수층의 지하수 내 염분침투를 예방하고 지속가능한 지하수자원 이용과 관리를 위해서 담수(freshwater)와 염수(saltwater)의 경계면 변동 및 분포를 추적하기 위해 지화학적·물리적 접근 연구는 매우 중요하다. 해안 및 하굿둑 주변지역의 담·염수 경계면 특성에 관한 조사 및 분석은 직·간접적 관측방법으로 이뤄진다. 직접 관측방법에는 심도별 전기전도도(Electrical condutivity, EC) 측정 방법으로 1개 정밀센서(CTD-Diver)를 와이어에 연결하여 수동으로 측정하는 방식과 자동 관측센서(CTD-Diver)를 공내 여러 관심심도에 이격거리에 따른 관측센서를 여러 개 설치하여 특정 관심심도별 연속적인 상태를 파악하는 방법이다. 간접 관측방법에는 관심심도별 또는 특정 심도의 지하수 수질을 채취하여 지화학 분석하는 방법이 있다. 직·간접적 관측방법에는 관심심도별 구간의 이격거리와의 공백과 조사시기에 대수층의 해수침투를 파악할 수 있으나 연속적 관측과 예측은 매우 어렵다. 직접적 관측방법 중 1개 정밀센서와 와이어를 이용한 특정 시·공간에 대한 연직 프로파일링 관측은 가능하나, 연속적 관측과 예측 또한 매우 어렵다. 따라서 본 연구는 물리적 방법을 기반으로, 하굿둑 인근의 담수(freshwater)와 염수(saltwater)와의 경계면 분포특성에 관한 연구는 하굿둑 실증시험 시 지하수 내 염분침투에 따른 지하수 환경 등에 대한 전반적 특성 등을 분석하였다. 하구언 일대의 지하수 내 염분침투 및 지하수 오염관리를 위해 실증실험 등 인위적인 요소와 강우, 태풍, 가뭄 등 자연적인 요소에 의해 여러 형태의 분포특성을 보이는 담수와 염수의 분포 및 혼합대의 특성에 따른 담염수의 경계면 형태별 분석을 하는데 목적을 두었다. 해수침투의 예방과 지하수 오염 방지를 위해 심도별 자동 수질측정 장치를 개발하고, 이를 이용한 염지하수의 프로파일링 감시 시스템을 제안하였다.

  • PDF

Eco-Concrete를 이용한 SRT System의 사면수처리부와 평면수처리부의 정화효율분석 (Purification Efficiency of Slop & Plane Water Treatment Part of SRT System Using Eco-Concrete)

  • 장원근;박재영;최이송;장준영;오종민
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2006년도 학술발표회 논문집
    • /
    • pp.1860-1864
    • /
    • 2006
  • 본 연구는 강우시 발생되는 강우유출수와 합류식하수관거월류수에 의해 하천으로 유입되는 오염부하를 저감시키기 위한 공법으로, 고수부지 및 제방사면부와 둔치부를 형상화하여 pilot를 제작하였고, 연속적으로 시운전을 한 SRTS(Stormwater Runoff Treatment System)에 관한 것이다. SRT system 내부의 사면수처리부와 평면수처리부에는 다공성 콘크리트를 충진하였다. system 상부에는 식생을 조성하여 뿌리가 수면에 닿아 영양물질을 흡수하는 목적으로 사면수처리부와 평면수처리부에 각각 정육각형과 직사각형인 식생포트를 탈.부착이 가능하도록 고안하였다. 내부에서는 토양과 수처리조 사이에 연결관을 부착하였고, 모세관현상에 의해 토양이 수분을 흡수하도록 구성하였다. pilot plant는 유입부, 사면 수처리부, 평면 수처리부, 유출부로 나누었다. 유입부는 유입펌프와 V-notch로 구성하였고, 유입펌프는 2대를 설치하여 1시간 간격으로 연속적 유입으로 유량조절이 가능하도록 상호교대 운전을 하였다. 평면 수처리부$(W(1.0m){\times}(L(2.4m){\times}H(0.6m))$는 장방형의 접촉산화조로서 하부에 슬러지 침전 및 저류를 위한 hopper를 설치하여 슬러지의 원활한 수집 및 인발이 가능하도록 하였다. 유출부는 사각weir를 설치하였다. 강우유출수의 pH는 $7.27{\sim}7.92$이고, DO농도는 $7.12{\sim}7.88mg/l$로 관측되었다. 2차처리수의 pH는 평균7.4이고 DO농도는 최저 4.5 mg/l에서 최고 8.9 mg/l로 평균 6.8 mg/l로 관측되었다. 또한 강우유출수의 유입수의 T-N, T-P 농도는 각각 $17.5{\sim}22.5mg/l,\;8.9{\sim}11.4mg/l$의 범위이고, 2차 처리수의 유입수의 T-N, T-P 농도와 유사하였다.적인 방법론을 제시할 수 있을 것으로 사료된다.첨두홍수량을 저류하기 위해서 상대적으로 넓은 저류면적이 필요한 것으로 나타난다. 대등한 수위감소값의 홍수저감효과를 발휘하기 위해서 본 연구에서는 On-Line 저류지 면적은 Off-Line 저류지에 비 두배 이상이 필요한 것으로 보여졌다.들에 관한 정보는 종종 현장관측에서 조차 무시되는 경우가 많다. 이에 본 연구에서는 수질모형의 매개변수 중 특히 수리특성에 관련된 매개변수들이 수질에 미치는 영향을 파악하는 것을 목적으로 하고 있다. 이를 위해 적용된 수질모형은 QualKo를 사용하였으며, 대상 하천은 낙동강 본류 경남구간 시점 부근인 회천 합류 전부터 낙동강 본류 경남구간 종점 부근인 밀양강 합류 전까지의 경남 오염총량관리 기본계획 시 구축된 모형 매개변수를 바탕으로 분석을 수행하였다. 일차오차분석을 이용하여 수리매개변수와 수질매개변수의 수질항목별 상대적 기여도를 파악해 본 결과, 수리매개변수는 DO, BOD, 유기질소, 유기인 모든 항목에 일정 정도의 상대적 기여도를 가지고 있는 것을 알 수 있었다. 이로부터 수질 모형의 적용 시 수리 매개변수 또한 수질 매개변수의 추정 시와 같이 보다 세심한 주의를 기울여 추정할 필요가 있을 것으로 판단된다.변화와 기흉 발생과의 인과관계를 확인하고 좀 더 구체화하기 위한 연구가 필요할 것이다.게 이루어질 수 있을 것으로 기대된다.는 초과수익률이 상승하지만, 이후로는 감소하므로, 반전거래전략을 활용하는 경우 주식투자기간은 24개월이하의 중단기가 적합함을 발견하였다. 이상의 행태적 측면과 투자성과측면의 실증결과를 통하여 한국주식시장에 있어서 시장수익률을 평균적으로 초과할 수 있는 거래전략은 존재하므로 이러한 전략을 개발 및 활용할 수 있으며, 특히, 한국주식시장에 적합한 거래전략은 반전거래전략이고, 이 전략의 유용성은 투자자가 설정한 투자기간보다 더욱 긴 분석기간의 주식가격정보에 의하여 최대한 발휘될 수 있음을 확인하

  • PDF

건강과 자연농업-제230호

  • 정진영
    • 건강과 자연농업
    • /
    • 230호
    • /
    • pp.1-12
    • /
    • 2006
  • 2006년도 지회장단 연찬회 성료/미국산 유전자조작 쌀 세계로 퍼저나가/우리의 산야초/연속 강우에 따른 과일나무 잎의 갈변에 대하여/벼 멀칭재배의 이론과 기술/06'대구.경북도지부 지회장단 연찬회/생활습관으로 노화를 최소화하라/철원지회,친환경농산물 가공센터 준공/고령.우곡지회,그린수박 선별장 준공/아이들건강과 먹을거리 패러다임의 전환/미국에서 친환경 상품 소비 증가세/미국 LA에 한국인 대상 유기농 마켓 등장/아마존 닷컴 유기농식품 판매에 손길/두가지 비료로 유기농작물의 고품질 다수확 생산을 체계적으로 보다 쉽게 하는 방법/경북,10년간 1조7천억 투입해 친환경농업 육성/벼희잎마름병의 방제는 이렇게/빛을 이용한 해충의 방제/유기농업연수원 매각 한 후 본회 사옥 신축/화성지역 인증자 후기관리교육 실시/국내외 친환경농축산물의 생산 및 인증실태/친환경 유기낙농의 정착 방안/8월중 새식구 명단/2006 지회장단 연찬회 이모저모/유기 및 일반 농산물의 품질 비교에 관한 연구/SK건설, 유기농 주말농장 무상대여 서비스/항공사들, 유기농식단 등으로/중국에도 웰빙바람, 유기농이뛴다/유기질비료의 이해

  • PDF