Journal Search Engine
Search Advanced Search Adode Reader(link)
Download PDF Export Citaion korean bibliography PMC previewer
ISSN : 2288-0992(Print)
ISSN : 2288-100X(Online)
Protected Horticulture and Plant Factory Vol.22 No.4 pp.328-334
DOI : https://doi.org/10.12791/KSBEC.2013.22.4.328

온실 내 잉여 태양열을 이용한 공기열원 히트펌프 성능향상

권진경*,강금춘,문종필,강연구,김충길,이수장
농촌진흥청 국립농업과학원 에너지환경공학과
본 연구에서는 온실 내부의 태양 잉여열과 외부의 공기열을 선택적으로 열원으로 이용함으로써 히트펌프의 성능을 향상시키고, 온실의 환기 지연을 통해 이산화탄소 시용비용을 절감할 수 있는 온실 공조시스템을 개발하고자 하였다. 본 시스템의 축열 과정은 태양 잉여열을 이용하는 내부순환모드와 외기열을 이용하는 외부순환모드가 온실 내부온도에 따라 자동으로 절환되도록 구성하였으며, 히트펌프 가동, 축열모드 절환, 난방 가동을 위한 6개의 온도값을 입력함으로써 축열과 난방이 자동으로 수행되도록 설계하였다. 단동온실을 대상으로 무환기조건에서 기초시험을 수행한 결과, 태양 잉여열을 이용한 축열은 약 11시부터 시작되어 평균 3시간 30분 정도 유지되었으며, 주간의 온실 내부온도는 환기를 수행하지 않음에도 대부분 약 20~28℃ 범위를 유지하였다. 주간 내부순환모드에서 시스템의 난방성능계수는 약 3.35로 야간 외부순환모드의 2.46 및 주간 외부순환모드의 2.67에 비해 각각 36% 및 25% 향상됨을 확인하였다. 본 시스템의 개선사항으로 태양 잉여열의 효율적 이용을 위해 축열조 관리온도를 상승시킬 수 있는 고효율 히트펌프의 적용이 필요하며, 온실의 무환기 운용에 따른 과습환경의 조성을 방지하고 태양 잉여열 수준이 높은 시기에 온실의 온도상승을 방지하기 위해 강제환기를 운전모드에 추가할 필요가 있는 것으로 판단되었다.

Performance Improvement of an Air Source Heat Pump by Storage of Surplus Solar Energy in Greenhouse

Jin Kyung Kwon*, Geum Chun Kang, Jong Pil Moon, Youn Ku Kang, Chung Kil Kim, Su Jang Lee
Energy & Environmental Engineering Division, National Academy of Agricultural Science, RDA, Suwon 441-707, Korea
Received September 25, 2013; Revised October 11, 2013; Accepted October 15, 2013

Abstract

A greenhouse heating system to improve heat pump performance using inside and outside air of greenhouseas a heat source selectively and cut CO2 enrichment costs by delay of greenhouse ventilation was developed.In this system, thermal storage modes divided into inside circulation mode using surplus solar energy and outside circulationmode using outside air heat. The thermal storage modes were designed to be switched mutually according toinside greenhouse temperature and six temperature values were input to control the heat pump operating, thermalstorage mode switching and greenhouse heating automatically. Operating characteristics of this system were tested ina plastic greenhouse of non-ventilation condition. The results of test showed that the inside circulation mode began atabout 11:00 and lasted for about 210 minutes and inside greenhouse temperature was maintained between 20~28℃in spite of non-ventilation. System heating COP of the inside circulation mode in the daytime was 3.35, which was36% and 25% higher than that of the outside circulation modes in the nighttime and daytime respectively.

 

0023-01-0022-0004-7.pdf3.30MB

서 론

 2011년 기준 국내 시설원예 전체면적 52,393ha 중 가온면적은 16,263ha로, 2008년과 비교할 때 각각 53,408ha, 13,329ha로써 전체면적의 감소에도 불구하고 가온면적의 비중은 꾸준히 증가하고 있다(MAFRA, 2012a, b). 국내 가온온실의 난방연료는 약 90%를 유류에 의존하고 있어 국제유가의 변동에 취약한 구조를 가지고 있으며, 생산비 중 광열비의 비중이 높아 시설농가의 경영압박 요인으로 작용하고 있다. 이에 농림축산식품부는 시설원예농가의 난방비 부담 경감을 위하여 2010년부터 지열원 히트펌프를, 2012년부터 공기열원 히트펌프를 각각 보급대상기종에 포함하여 설치 지원사업을 추진하고 있다(MAFRA, 2012c). 히트펌프는 원예시설의 난방, 야간냉방, 제습 등의 다양한 환경관리에 적용 가능한 다목적공조기기로서(Gracia 등, 1998; Willits와 Gurjer, 2004), 국내의 경우 원예시설용 히트펌프와 관련한 연구는 대부분 지열원 히트펌프를 대상으로 수행되고 있다(Kang 등, 2007; Ryou 등, 2008). 지열원 히트펌프는 안정적인 공조 성능을 가지나 지중 열교환기 매설비용 등 초기 설비투자비가 높은 단점이 있다. 반면 공기를 열원으로 하는 공기열원 히트펌프는 설치비가 상대적으로 낮고 소형화가 가능하여 설치장소의 제약을 적게 받는 장점이 있으나, 외기온에 따라 성능변화가 커지는 단점이 있다(Marsh와 Singh, 1994). 원예시설용 공기열원 히트펌프는 일본을 중심으로 활발한 연구가 수행되었으며 초기설비비와 운용비를 절감하기 위해 히트펌프와 보조 난방기를 병용 운전하는 하이브리드 운용기술이 개발되어 다양한 시설재배에 적용되고 있다(Kawashima 등, 2008; Kawashima 등, 2011; Tong 등, 2011).

 한편 온실은 동절기에도 태양복사로 인해 내부온도가 대부분 작물의 생육적온 범위를 초과하게 되므로 환기를 실시하여 온실 내부의 열을 외부로 배출시킨다. 이와 같이 환기에 의해 버려지는 태양 잉여열을 회수, 축열하여 야간 난방에 활용함으로써 무가온 온실의 작기 연장이나 가온온실의 난방부하 경감에 활용하고자 하는 연구가 수행되었다(Beshada 등, 2006; Lee 등, 2011). 온실에서 발생하는 태양 잉여열을 수치모델로 분석한 연구결과에 의하면 지역별로 연간 태양 잉여열은 연간 난방부하의 약 50~200%에 달해 난방 열원으로 충분한 활용가치가 있는 것으로 평가되었다(Suh 등, 2011). 또한 작물재배 시 증수를 위해 이산화탄소를 시용하는 온실의 경우 고온피해를 방지하기 위해 실시하는 환기가 이산화탄소의 손실로 이어져 시용효과의 경감과 시용비용이 증가하게 되므로 효율적인 태양 잉여열의 축열은 난방에너지 절감뿐 아니라 환기 지연을 통한 이산화탄소 시용비용의 절감효과도 기대할 수 있다.

 본 연구에서는 공기열원 히트펌프의 성능 향상과 환기지연이 가능한 온실 공조시스템을 개발하기 위한 기초연구로서 온실내부의 태양 잉여열과 외부 공기열을 선택적으로 집열 및 축열하여 온실 난방에 사용할 수 있는 공기열원 온실 공조 시스템을 설계하였으며, 온실에 대한 적용시험을 통해 히트펌프의 성능개선 효과 및 시스템가동 시의 온실 내부 환경변화를 분석하였다.

재료 및 방법

 태양 잉여열 및 외기열 축열식 히트펌프 공조시스템의 설계 및 성능시험 대상 공시온실을 Fig. 1에 점선으로 나타내었다. 공시온실은 수원시 입북동에 소재한 07-단동-3형의 작물을 재배하지 않는 비닐하우스(폭 7.0m, 측고 1.4m, 동고 3.3m, 길이 50m)로 폴리에틸렌 필름으로 2중 피복되어 있으며 보온커튼은 설치되어 있지 않다. 공조시스템의 용량을 결정하기 위해 식(1)을 이용하여 온실의 최대 난방부하를 계산하였다. 설계외기온은 수원의 최근 5년 최저기온 평균값인 −15.9℃(KMA, 2013)를, 난방설정온도는 15℃를 기준으로 하였으며, 폴리에틸렌 필름 2중 피복의 열관류율은 3.6W/m2  · ℃, 환기전열계수는 0.35W/m2  · ℃, 지면전열계수는 0.28W/m2  · ℃, 설계지중온도는 -4.3℃, 보온피복 열절감률은 0, 풍속보정계수는 1.0을 사용하였다. 계산한 결과 해당 온실의 최대난방부하는 약 73kW였으며 이를 히트펌프의 냉동톤으로 환산하면 약 20RT가 된다.

 

Fig. 1. Experimental greenhouse in Suwon.

 여기서, Qg: 온실 최대난방부하(W), Ag: 온실 피복면적(m2), As: 온실 상면적(m2), ht: 열관류율(W/m2 · ℃), hv: 환기전열계수(W/m2 · ℃), hs: 지면전열계수(W/m2 · ℃), Ts: 난방설정온도(℃), Td: 설계외기온(℃), Tg: 설계지중온도(℃), fr: 보온피복 열절감율, fw: 풍속보정계수

 본 연구에서 주간에는 온실내부 태양 잉여열과 외기열을 주로 이용하는 전축열 운전을, 야간에 축열조의 열을 이용한 온실 난방과 외기열에 의한 축열을 동시에 수행하는 부분축열 운전을 모두 이용하여 24시간 가동하는 온실 난방시스템을 고려하였으므로 히트펌프는 계산된 최대 난방부하의 50%인 35kW(10RT)급 공기 대 물 히트펌프(삼영종합기기, SY-HS010AC)로 선정하였다.

 본 연구에서 개발하고자 하는 온실 난방시스템의 개략도를 Fig. 2에 나타내었다. 시스템의 축열운전은 온실내부의 태양 잉여열을 히트펌프의 열원으로 사용하기 위해 온실 상부의 공기를 흡기하여 히트펌프 증발기를 통과시킨 후 다시 온실 내부로 배기시키는 내부순환모드와 온실 외부의 공기열을 열원으로 사용하기 위해 온실 외부공기를 흡기하여 다시 외부로 배기하는 외부순환모드로 대별된다.

Fig. 2. Schematic of greenhouse heating system using surplus solar heat and outside air heat.

 각각의 운전모드에서 공기유동을 분리하기 위해 히트펌프 상부에 내부순환 및 외부순환용 덕트를 설치하였다. 덕트는 600 × 600mm 단면의 아연도금강판 재질로 각 덕트에는 순환팬과 모드 절환용 댐퍼를 설치하여 각각의 운전모드에서 해당 덕트의 순환팬 가동과 댐퍼 개방이 이루어진다. 내부순환모드에서 공기는 온실 마구리벽 쪽의 흡기구에서 흡기되고 길이 15m, 직경 500mm의 비닐덕트를 통해 온실 내부로 배기된다. 히트펌프를 이용하여 태양 잉여열 및 외기열로부터 흡수한 열량은 축열펌프로 축열조에 저장한 후 야간에 방열펌프와 팬코일유닛을 이용하여 온실내부로 공급한다. 축열조는 디퓨져가 설치되지 않은 혼합 개방형으로 용량은 15t이며, 온실내부의 팬코일유닛은 9.3kW 용량으로 온실 길이방향으로 8대를 설치하였다.

 축열운전을 수행하는 동안 외부순환모드와 내부순환모드의 상호절환은 온실 내부에서 실시간으로 측정한 온도를 기준으로 제어기가 해당 모드의 댐퍼를 개방하고 순환팬을 가동시킴으로써 수행된다. Table 1에는 본 시스템의 구성 요소들에 대한 사양을 나타내었으며, Fig. 3에는 상기의 단동 비닐하우스에 실제 설치한 외부의 히트펌프, 순환시스템, 축열조 등과 내부의 팬코일유닛 등을 나타내었다.

Table 1. Components of present system.

Fig. 3. Experimental facilities.

 본 시스템의 축열 성능을 분석하기 위해 히트펌프의 응축기 입, 출구 물배관과 축열조 내부 중앙 1.2m 높이에 온도센서(PT100, Ahlborn, Germany)를 설치하였으며, 출구 물배관에는 유량계(FVA915, Ahlborn, Germany)를 설치하였다. 내부순환 및 외부순환용 덕트의 입, 출구와 온실 내 중앙부의 0.5m, 2.5m 높이에 각각 온습도센서(FHA646, Ahlborn, Germany)를 설치하여 히트펌프 증발기 통과 전후의 공기 온습도와 축열 및 난방 시의 온실내부 온습도 변화를 측정하였다. 태양 잉여열의 분석을 위하여 온실 내, 외부에 각각 일사량계(LP02, Hukseflux, Netherlands)를 설치하였으며, 전체 시스템의 난방성능계수를 측정하기 위하여 시스템의 주전원부에 전력량계(CW240, Yokogawa, Japan)를 설치하였다. 상기의 모든 데이터는 1분 간격으로 측정하였으며, 모든 측정값을 데이터수집기(Almemo, Ahlborn, Germany)로 수집하여 분석하였다.

 본 시스템의 난방성능계수는 전체 시스템의 소비전력과 응축기 입, 출구의 물온도차 및 유량을 이용하여 아래의 식(2)에 의해 계산하였다.

 

 여기서 COPh: 시스템 난방성능계수, PHP: 시스템 소비전력(kW), ρw: 열전달매체 밀도(kg/m3), Qw: 열전달매체유량(m3/s), cw: 열전달매체의 비열(kJ/kg · ℃), Tw,i, Tw,o: 각각 응축기 입, 출구 열전달매체 온도(℃)

 본 시스템의 축열 및 난방을 위한 운전 알고리즘을 Fig. 4와 같이 설계하였다. 히트펌프와 축열펌프는 축열조의 상한 설정온도에서 정지하고, 하한 설정온도에서 가동하도록 하였으며, 상, 하한 설정온도는 히트펌프의 토출가스 고온 조건을 고려하여 각각 45℃, 40℃로 설정하였다. 태양 잉여열을 이용하는 내부순환모드 축열과외기열을 이용하는 외부순환모드 축열은 온실 내부의 온도를 기준으로 자동절환 되도록 하였다. 시설원예에서 가온재배를 하는 채소류 및 화훼류는 대부분 20~30℃의 생육적온 범위를 가지며 최고한계온도는 약 35℃ 정도이므로 모드 전환 시 온도변화는 지연시간을 가지고 후 행하는 조건을 고려하여 일출 후 온실 내부온도가 27℃가 되기 전까지는 외부순환모드로 축열하며, 27℃에 도달하면 다시 22℃로 하강할 때 까지 내부순환모드로 자동절환 되도록 하였다. 이와 같이 온실 내부온도에 따른 축열모드 절환을 통해 상대적으로 고온인 태양 잉여열을 열원으로 이용함으로써 히트펌프의 성능을 높이고, 무환기 조건에서도 온실 내부 온도를 작물의 생육적온 범위내에서 유지할 수 있도록 하였다. 가온재배 시설작목의 야간 설정온도 범위인 약 8~18℃의 중간정도 온도를 고려하여 난방 설정온도는 14℃를 기준으로 온실 내부 온도가 13℃에 도달하면 팬코일유닛과 방열펌프가 작동하고 15℃에서는 정지하도록 하였다.

Fig. 4. Schematic of operating algorism of experimental system.

 본 시스템은 히트펌프 가동, 축열모드 절환, 난방 가동을 제어하기 위해 각각 2개씩, 총 6개의 온도값을 제어기에 입력하면 시스템의 축열 및 난방 운전이 자동으로 수행되도록 하였으며, 본 성능시험에 적용한 제어용 온도값들을 Table 2에 나타내었다.

Table 2. Temperature setup for operation of experimental system.

결과 및 고찰

 공시온실에 대해 Table 2의 제어 온도값을 적용하여 축열 및 난방시험을 수행한 결과를 Fig. 5에 나타내었다. 시험은 2013년 2월 20일부터 26일까지 7일간 수행하였으며 시험 기간 중 온실은 주야간 모두 환기를 하지 않는 밀폐형 온실로 운용하였다. 측정 데이터는 축열조 온도가 안정화되고 청명일이 유지된 2월 22일 18시부터 25일 16시까지 약 3일간의 결과를 나타내었다.

Fig. 5. Analysis of system performance and measurements of greenhouse environment by experiments.

 Fig. 5(a)에는 온실 내, 외부온도 및 히트펌프 응축기 입출구 온도차를 나타내었으며, 위쪽에는 축열운전의 내, 외부순환모드 구간을 표시하였다. 시험기간 3일간 평균외기온은 −0.4℃(최저 −7.4℃, 최고 9.4℃), 평균 외부상대습도는 63%였다. 팬코일유닛과 방열펌프가 가동되어 야간에 난방을 수행한 시간은 3일간 하루 평균 15시간 15분(17시 45분~익일 9시)이었다. 시스템의 축열운전은 야간 및 오전의 외부순환모드 과정을 거쳐 2월 23일 10시 57분에 온실 내부온도가 27℃에 도달함으로서 내부순환모드로 절환되었으며, 11시 05분에 온실 내부온도가 22℃로 하강하며 다시 외부순환모드로 절환되었다. 이 후 축열모드 간 절환을 5회 반복한 후 일사량이 최대에 근접하는 12시부터 3시간 32분 동안 내부순환모드로 축열을 수행하였으며, 약 28분간의 외기순환모드 축열과정을 거친 후 16시 1분에 축열조 온도가 상한온도인 45℃에 도달하여 히트펌프는 37분간 가동을 중단하였다. 이 후 축열조 온도가 하강함에 따라 다시 외기순환모드 축열과정을 다시 시작하여 익일 온실 내부온도가 27℃에 도달할 때까지 외부순환모드를 유지하였다. 이와유사한 축열과정이 24일과 25일에도 반복되었으며 24일에는 최초 10시 48분부터 3시간 25분간, 25일에는 최초11시 27분부터 3시간 4분간 내부순환모드로 태양 잉여열을 이용한 축열을 수행하였다. 이 때 주간의 온실은 환기를 실시하지 않음에도 불구하고 내부순환모드를 통해 흡기온도보다 약 9℃ 정도 낮아진 저온공기가 온실내부로 배기되어 내부온도는 작물생육적온 범위 이내인 20~28℃를 유지하였으나, 25일 주간에는 온실 내부의 잉여열이 존재함에도 불구하고 축열조 온도가 상한온도에 일찍 도달하여 히트펌프는 가동을 중지하였다. 이로 인해 온실 내부온도는 일반적인 작물의 최고한계온도인 35℃까지 상승하였으며 이와 같은 현상은 태양 잉여열수준이 상대적으로 높은 시기에 더욱 심화될 것으로 예상되었다. 태양 잉여열의 효율적 이용과 온실 온도관리를 위해서는 축열조의 용량을 가변적으로 운용하거나 축열조 관리온도의 범위를 넓게 운용할 필요가 있으며, 현실적으로는 중온수 제조용 고효율 히트펌프를 채용하여 축열조 관리온도를 상승시키는 방법으로 시스템을 개선하는 것이 타당할 것으로 판단되었다. 응축기 입출구 온도차는 내부순환모드의 경우 평균 약 6.7℃로 외부순환모드의 약 4.5℃보다 49% 높아 응축기 토출열량이 상승하였으며, 외부순환모드의 경우 약 40분간의 축열운전과 약 6분간의 제상운전이 반복됨을 알 수 있다.

 Fig. 5(b)에는 시스템의 전력소비량과 응축기 토출열량의 측정치를 나타내었다. 상대적으로 고온의 공기를 히트펌프 열원으로 사용하는 내부순환모드가 외부순환모드에 비해 응축기 토출열량 뿐 아니라 전력소비량도 높게 나타난 것은 증발기의 열교환 증가가 압축기의 압축부하증가를 수반하기 때문이다. 그러나 내부순환모드의 응축기 토출열량이 외부순환모드보다 약 24~49% 상승한 반면 전력소비량의 증가는 약 8~14%에 그쳐 내부순환모드에서 시스템의 성능이 향상되는 것으로 나타났다.

 Fig. 5(c)에는 식(2)를 통해 계산한 시스템의 난방성능계수를 나타내었다. 난방성능계수는 야간의 외부순환모드에서 약 2.46, 일출 후 주간의 외부순환모드에서 약 2.67 정도의 평균값을 나타내었으며, 주간의 내부순환모드에서는 평균 약 3.35로 야간 및 주간의 외부순환모드에 비해 각각 36% 및 25% 증가함을 확인하였다. 내부순환모드에서 난방성능계수는 응축기 입구온도가 증가할수록 감소하였으며 이는 응축기의 열교환 감소로 인한 압축기의 압축부하 증가가 Fig. 5(b)에 나타낸 바와 같이 전력소비량을 증가시켰기 때문이다.

 Fig. 5(d)에는 온실 내부의 상대습도와 온실 내, 외부의 일사량의 변화를 나타내었다. 3일간 온실 내부의 상대습도는 평균 76%였으며, 일몰 약 1시간 전부터는 90%를 상회하고 야간 난방 시간대에도 85%에 근접하는 과습 환경이 조성되는 것으로 나타났다. 이는 축열기간 동안 온실을 무환기로 운용함으로서 온실 내 토양에서 증발한 수분이 내부에 축적된 결과이며, 실제 작물재배시 작물의 증산량이 더해지면 온실 내부는 더욱 과습해질 것으로 예상된다. 온실 내, 외부의 일사량 측정치로 계산한 3일간의 평균 광투과율은 62%였으며, 일출 후 온실 내부의 과습 환경에 의해 피복재 표면에 발생한 응축물방울이 광투과율과 온실의 온도상승을 저해하는 요인으로 판단되었다. 따라서 본 시스템의 작물재배 적용 시, 온실 내부의 상대습도가 급증하고 태양 잉여열 수준이 낮아지는 오후 일몰 전 시간대에 내부순환용 덕트의 흡기구와 외부순환용 덕트의 배기구를 직결하여 온실 내부의 습공기를 외부로 강제 배출하는 배습모드를 추가할 필요가 있는 것으로 판단되었다. 이는 태양 잉여열 수준이 히트펌프의 성능을 초과하여 내부순환모드의 가동에도 불구하고 온실 내부온도가 작물 생육적온범위를 초과할 수 있는 시기에 강제환기에도 활용할 수 있을 것이다. 일반적으로 이산화탄소는 작물의 광합성속도가 높은 오전 시간대를 대상으로 일출 후 약 4시간 정도를 시용하므로 본 시스템의 내부순환모드 축열을 통해 보통의 온실에서 오전 중에 개시되는 환기를 지연시켜 오후 시간대에 배습 환기를 수행함으로써 시용된 이산화탄소의 손실과 온실의 과습을 방지할 수 있을 것으로 판단되었다.

사 사

 본 연구는 농촌진흥청 국립농업과학원 농업과학기술연구개발사업(과제번호: PJ00853001)의 지원에 의해 이루어진 것임.

Reference

1.Beshada, E., Q. Zhang, and R. Boris. 2006. Winter performance of a solar energy greenhouse in southern Manitoba. Canadian Biosystems Engineering 48:5.1-5.8.
2.Gracia, J.L., De la Plaza, L.M. Narvas, R.M. Benavente, and L. Luna. 1998. Evaluation of the feasibility of alternative energy sources for greenhouse heating. J. Agric. Engng Res. 69:107-114.
3.Kang, Y.K., Y.S. Ryou, G.C. Kang, Y. Paek, and Y.J. Kim. 2007. Heating performance of horizontal geothermal heat pump system for protected horticulture. J. of Biosystems Eng. 32(1):30-36 (in Korean).
4.Kawashima, H., M. Takaichi, M. BaBa, K. Yasui, and Y. Nakano. 2008. Effects of energy saving and the reduction of carbon dioxide emissions with a hybrid-heating system using an air-to-air heat pump for greenhouse heating. Bulletin of the National Institute of Vegetable and Tea Science 7:27-36 (in Japanese).
5.Kawashima, H., M. Takaichi, and K. Yasuba. 2011. Performance of air-to-air heat pump for tomato greenhouse cooling and reduction of night-cooling load by mulch to control soil heat flux. Bulletin of the National Institute of Vegetable and Tea Science 10:95-104 (in Japanese).
6.Korea Meteorological Administration (KMA), 2013. Climate data-Annual data extreme value. http://www.kma.go.kr/weather/climate/extreme_yearly.jsp (in Korean).
7.Lee, S.H., Y.S. Ryou, J.P. Moon, N.K. Yun, J.K. Kwon, S.J. Lee, and K.W. Kim. 2011. Solar energy storage effectiveness on double layered single span plastic greenhouse. J. of Biosystems Eng. 36(3):217-222 (in Korean).
8.Marsh, L.S. and S. Singh. 1994. Economics of greenhouse heating with a mine air-assisted heat pump. Trans. of the ASAE 37(6):1959-1963.
9.Ministry of Agriculture, Food and Rural Affairs (MAFRA). 2012a. Greenhouse status for the vegetable grown in facilities and the vegetable productions in 2011. ed. Gwacheon, Korea (in Korean).
10.Ministry of Agriculture, Food and Rural Affairs (MAFRA). 2012b. Cultivation status of floricultural crop in 2011. ed. Gwacheon, Korea (in Korean).
11.Ministry of Agriculture, Food and Rural Affairs (MAFRA). 2012c, Agriculture, forestry and fisheries enterprise enforcement manual. http://manual.agrix.go.kr/home/index.php (in Korean).
12.Ryou, Y.S., Y.K. Kang, G.C. Kang, Y.J. Kim, and Y. Paek. 2008. Cooling performance of horizontal type geothermal heat pump system for protected hoticulture. Journal of Bio- Environment Control 17(2):90-95 (in Korean).
13.Suh, W.M., Y.H. Bae, Y.S. Ryou, S.H. Lee, H.T. Kim, Y.J. Kim, and Y.G. Yoon. 2011. Estimation of surplus solar energy in greenhouse (II). Journal of Bio-Environment Control 20(2):83-92 (in Korean).
14.Tong, Y., T. Kozai, N. Nishioka, and K. Ohyama. 2011. Greenhouse heating using heat pumps with a high coefficient performance (COP). Biosystems Engineering 106:405- 411.
15.Willits, D.H. and Y.R. Gurjer. 2004. Heat pumps for the heating and night-cooling of greenhouse crops: a simulation study. Trans. of the ASAE 47(2):575-584.