WO2013100345A1 - 냉난방부하 추정을 위한 시간별 기상데이터 예측방법 - Google Patents

냉난방부하 추정을 위한 시간별 기상데이터 예측방법 Download PDF

Info

Publication number
WO2013100345A1
WO2013100345A1 PCT/KR2012/009001 KR2012009001W WO2013100345A1 WO 2013100345 A1 WO2013100345 A1 WO 2013100345A1 KR 2012009001 W KR2012009001 W KR 2012009001W WO 2013100345 A1 WO2013100345 A1 WO 2013100345A1
Authority
WO
WIPO (PCT)
Prior art keywords
dimensionless
hourly
relative humidity
weather data
temperature
Prior art date
Application number
PCT/KR2012/009001
Other languages
English (en)
French (fr)
Inventor
유성연
윤홍익
한규현
김태호
Original Assignee
(주)가교테크
충남대학교산학협력단
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by (주)가교테크, 충남대학교산학협력단 filed Critical (주)가교테크
Priority to US14/369,639 priority Critical patent/US20140365128A1/en
Publication of WO2013100345A1 publication Critical patent/WO2013100345A1/ko

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01WMETEOROLOGY
    • G01W1/00Meteorology
    • G01W1/10Devices for predicting weather conditions
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K17/00Measuring quantity of heat
    • G01K17/06Measuring quantity of heat conveyed by flowing media, e.g. in heating systems e.g. the quantity of heat in a transporting medium, delivered to or consumed in an expenditure device
    • G01K17/08Measuring quantity of heat conveyed by flowing media, e.g. in heating systems e.g. the quantity of heat in a transporting medium, delivered to or consumed in an expenditure device based upon measurement of temperature difference or of a temperature
    • G01K17/20Measuring quantity of heat conveyed by flowing media, e.g. in heating systems e.g. the quantity of heat in a transporting medium, delivered to or consumed in an expenditure device based upon measurement of temperature difference or of a temperature across a radiating surface, combined with ascertainment of the heat transmission coefficient
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Ecology (AREA)
  • Atmospheric Sciences (AREA)
  • Environmental Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Combustion & Propulsion (AREA)
  • Chemical & Material Sciences (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Air Conditioning Control Device (AREA)

Abstract

본 발명은 냉난방부하 추정을 위한 시간별 기상데이터 예측방법에 관한 것으로, 기상청으로부터 과거의 기상데이터를 확보하는 기상데이터 획득단계(S100)와; 상기 기상데이터 획득단계(S100)로부터 획득된 기상데이터를 분석하여 필요한 데이터를 추출하는 기상데이터 분석 및 추출단계(S200)와; 상기 기상데이터 분석 및 추출단계(S200)에서 추출된 기상데이터를 무차원화시켜 무차원값을 산출하는 무차원값 산출단계(S300)와; 상기 무차원값 산출단계(S300)에 의해 산출된 무차원값으로부터 상관관계를 나타내는 상관식 결정단계(S400) 및; 상기 상관식 결정단계(S400)에서 구한 시간별 무차원값으로부터 익일 시간별 기상데이터를 예측하는 익일 시간별 기상데이터 예측단계(S500)로 이루어지되, 상기 익일 시간별 기상데이터 예측단계(S500)에서 사용되는 최고, 최저 상대습도 및 일사량은 퍼지 알고리듬에 의해 추정하는 것을 특징으로 하는 난방부하 추정을 위한 시간별 기상데이터 예측방법으로 이루어지는 것을 특징으로 한다. 본 발명은 상기와 같이 구성에 의해 실측값을 사용하지 않고도 기상청에서 제공하는 최고 온도와 최저 온도만으로 시간별 외기온도, 상대습도 및 일사량을 정확하게 예측할 수 있다.

Description

냉난방부하 추정을 위한 시간별 기상데이터 예측방법
본 발명은 냉난방부하 추정을 위해 사용되는 시간별 기상데이터 예측방법에 관한 것으로, 더욱 상세하게는 건물의 냉난방부하를 산출할 때 사용되는 시간별 외기온도, 상대습도 및 일사량 등을 기상청에서 제공하는 최고온도와 최저온도에만 근거하여 정확하게 예측함으로써 냉난방 시스템을 효과적이며 경제적으로 운전할 수 있도록 하는 시간별 기상데이터 예측방법에 관한 것이다.
건물의 설계 단계에서 냉난방 시스템의 용량을 결정하기 위해서는 건물의 냉난방 부하를 필히 계산하여야 하며, 또한 건물의 공조설비 운전 단계에서도 운전에 필요한 에너지 소비량과 비용 등에 기초한 효율적이고 경제적인 제어를 위해서는 냉난방 부하에 대한 예측이 중요하다.
건물의 열부하는 실내 공간을 목표로 하는 환경으로 유지하기 위해 투입 또는 제거하여야 하는 열량으로서 이때 이러한 열부하를 계산하거나 예측하는 데에 가장 큰 영향을 주는 요소 중의 하나는 대상 건물이 위치하는 지역의 기상상태라고 할 수 있다.
따라서 건물의 열부하를 정확히 예측하기 위해서는 장기적이고 신뢰성이 있는 기상데이터를 확보하는 것이 무엇보다 중요하며, 이러한 기상 데이터가 부정확한 경우 열부하 계산 또는 예측 결과에 대한 신뢰도를 저하시킬 수 있다.
건물의 초기 설계단계에서 냉난방 시스템의 용량을 결정하기 위한 목적으로 냉난방부하를 계산할 때 필요한 기상데이터는 일반적으로 장기간의 측정데이터에 근거하여 작성되어 그 지역의 기후를 대표할 수 있는 표준기상데이터(typical climatic data)이다.
표준기상데이터는 통계처리방법 및 데이터의 구조에 따라 TRY, TMY, WYEC, TMY2, WYEC2 등과 같이 매우 다양한 형식으로 개발 및 제시되고 있으며 실제로 건물의 초기 설계단계에서 건물의 에너지 시뮬레이션에 사용되고 있는데, 이러한 표준기상데이터의 공통된 특징은 주요 기상요소들의 분포가 장기적인 분포와 가까운 일관성 있는 분포를 보인다는 것이다.
반면, 건물의 초기 설계단계와는 달리 실제 건물의 공조설비 운전단계에서 익일의 냉난방부하를 예측하여 냉난방 시스템의 운전에 적용하려면 다양한 외부환경에 대한 불확실성 및 운전 시나리오를 고려하여야 하며, 따라서 이러한 특성이 반영된 부하예측이 가능하도록 하기 위해서는 표준기상데이터와는 다른 기상데이터 형태가 필요하다. 즉 표준기상데이터 같은 일관성 및 주기성을 유지하는 동시에 실시간 외부환경의 변화를 고려한 시간별 기상데이터를 예측하는 것이 신뢰성 있는 익일 냉난방부하 예측의 필수조건이다.
이러한 이유로 익일 냉난방부하를 예측하기 위한 방법들이 다수 제안되고 있는데, 제안된 방법들은 대부분 기상예보를 통해 익일의 외기온도값을 얻은 다음, 이 값과 실제의 외기온도를 측정하여 그 차이를 반영한 온도 보정계수를 통해 업데이트시키는 방법을 채택하고 있으며, 일사량에 대해서는 고려하지 않고 있다는 단점이 있다.
이에 따라 본 발명자 등은 외기 온도를 실제로 측정하지 않고도 기상청에서 발표한 5년간의 기상데이터를 이용하여 익일 시간별 온도와 비습도를 예측하는 방법을 제안하여 특허등록(특허 제949044호, 발명의 명칭: 냉방시스템의 최적 운전방법)을 받은바 있는데, 이 방법은 종래에 비해 익일의 외기온도를 정확하게 예측할 수 있다는 장점은 있으나 상대습도 예측의 정확성이 다소 떨어지며 냉난방부하 계산에 필요한 일사량이 반영되지 않고 있다.
본 발명은 상기와 같은 종래의 기술이 가지는 문제점을 개선하기 위해 안출된 것으로, 본 발명은 외기온도 등의 기상상태 실제로 측정하지 않고도 기상청에서 제공하는 최고 온도와 최저 온도만으로 시간별 외기온도, 상대습도, 일사량(이하 이들을 "기상데이터"라 통칭한다)을 정확하게 예측할 수 있도록 한 것으로, 이를 통하여 익일의 냉난방 부하를 더욱 정확하게 예측할 수 있도록 하는 시간별 기상데이터 예측방법을 제공하는 것을 그 목적으로 한다.
상기와 같은 본 발명의 목적을 달성하기 위한 본 발명의 시간별 기상데이터 예측방법은 기상청으로부터 과거의 기상데이터를 확보하는 기상데이터 획득단계와; 상기 기상데이터 획득단계로부터 획득된 기상데이터를 분석하여 필요한 데이터를 추출하는 기상데이터 분석 및 추출단계와; 상기 기상데이터 분석 및 추출단계에서 추출된 기상데이터를 무차원화시켜 무차원값을 산출하는 무차원값 산출단계와; 상기 무차원값 산출단계에 의해 산출된 무차원값으로부터 상관관계를 나타내는 상관식 결정단계 및; 상기 상관식 결정단계에서 구한 시간별 무차원값으로부터 익일 시간별 기상데이터를 예측하는 익일 시간별 기상데이터 예측단계로 이루어지되, 상기 익일 시간별 기상데이터 예측단계에서 사용되는 최고, 최저 상대습도 및 일사량은 퍼지 알고리듬에 의해 추정하는 것을 그 기술적 특징으로 한다.
또한 본 발명은 외기온도(
Figure PCTKR2012009001-appb-I000001
), 상대습도(
Figure PCTKR2012009001-appb-I000002
), 일사량(
Figure PCTKR2012009001-appb-I000003
)을 무차원화시키는 것을 또 다른 기술적 특징으로 한다.
이에 더하여 본 발명은 무차원 외기온도, 무차원 상대습도, 무차원 일사량을 각각 시간과의 상관식에 의해 산출하는 것을 또 다른 기술적 특징으로 한다.
본 발명은 외기온도 등의 기상상태를 실제로 측정하지 않고도 기상청에서 제공하는 최고 온도와 최저 온도만으로도 시간별 외기온도, 상대습도, 일사량 등의 기상데이터를 예측하여 익일의 냉난방 부하를 정확하게 산출할 수 있는 효과가 있다.
또한 본 발명은 퍼지 알고리듬에 의해 최고/최저 상대습도 및 일사량을 추정함으로써 더욱 정확하고 합리적으로 시간별 기상데이터를 예측할 수 있다.
이에 더하여 본 발명은 외기온도, 상대습도, 일사량에 대해 무차원화(normalize)시켜 사용하기 때문에 주어진 기상데이터의 범위에 영향을 받지 않으며, 따라서 그 적용대상과 장소에 영향을 받지 않는다.
도 1은 본 발명에 따른 시간별 기상데이터 예측방법을 나타낸 순서도,
도 2는 외기온도 분포를 나타낸 그래프,
도 3은 상대습도 분포를 나타낸 그래프,
도 4는 일사량 분포를 나타낸 그래프,
도 5는 5년 동안의 시간별 무차원 외기온도 변화를 나타낸 그래프,
도 6은 5년 동안의 시간별 무차원 상대습도 변화를 나타낸 그래프,
도 7은 5년 동안의 시간별 무차원 일사량 변화를 나타낸 그래프,
도 8은 예측된 외기온도와 실제로 측정된 외기온도를 비교한 그래프,
도 9는 예측된 상대습도와 실제로 측정된 상대습도를 비교한 그래프,
도 10은 예측된 일사량과 실제로 측정된 일사량을 비교한 그래프,
도 11 내지 13은 각각 2008년의 7월과 8월에 대해 대전지역의 예측된 외기온도, 상대습도 및 일사량과 기상청에서 실제로 측정한 데이터를 비교한 그래프이다.
이하에서는 첨부된 도면을 참조하여 본 발명의 구성과 바람직한 실시예를 더욱 상세히 설명한다.
본 발명은 실측값을 사용하지 않고도 기상청에서 제공하는 기상데이터만으로 익일의 냉난방 부하를 정확하게 예측할 수 있도록 하는 시간별 기상데이터 예측방법에 관한 것으로 이를 위해 본 발명은 도 1에 도시된 바와 같이 크게 기상데이터 획득단계(S100), 기상데이터 분석 및 추출단계(S200), 무차원값 산출단계(S300), 상관식 결정단계(S400) 및 익일 시간별 기상데이터 예측단계(S500)를 포함한다.
(1) 기상데이터 획득단계(S100)
이 단계는 기상청으로부터 냉난방 부하계산에 필요한 기초적인 기상데이터를 획득하는 단계이다.
기상청에서는 전국의 기상 관측소에서 외기 온도, 상대습도, 일사량, 풍향, 풍속 등을 측정하여 기록하고 있는데, 이 중 온도와 상대습도는 1시간마다 순간치를 측정하여 온도는 ˚C, 상대습도는 %로 나타내고 있으며, 일사량은 일출시간으로부터 일몰시간까지 매 분별로 측정하여 1시간 적산한 것을
Figure PCTKR2012009001-appb-I000004
로 나타내고 있다.
본 발명에서는 냉방부하가 발생하는 6월, 7월, 8월, 9월 대전지역의 외기온도와 상대습도, 일사량을 분석하기 위해 기상청에서 측정했던 2003년부터 2007년까지의 5년간의 기상데이터를 활용하였다.
(2) 기상데이터 분석 및 추출단계(S200)
이 단계는 상기 기상데이터 획득단계(S100)에서 획득한 기상데이터를 분석하여 필요한 기상데이터만 추출하는 단계이다.
도 2 내지 도 4는 각각 대전지역의 2007년 7월과 8월의 2개월간의 외기온도, 상대습도 및 일사량을 나타낸 그래프를 나타낸 것이다.
도 2 내지 도 4를 살펴보면 외기온도의 경우 최고 온도는 13시부터 15시 사이에 나타나며, 새벽 4시부터 6시 사이에 최저 온도가 나타남을 알 수 있다. 또한 최고 온도와 최저 온도 사이에서 단조증가와 단조감소를 하면서 매일 일정한 양상으로 변화하고 있음을 알 수 있다.
상대습도의 변화를 온도와 관련시켜 살펴보면 외기온도가 높을 때에는 상대습도가 낮고, 반대로 외기온도가 낮을 때는 상대습도가 높은 것이 일반적인 경향임을 알 수 있다.
일사량의 경우 일출시간인 5시와 일몰시간인 19시를 경계로 11시와 12시 사이에 최고값을 가진다.
결과적으로 외기온도와 상대습도, 일사량은 모두 최대값과 최소값 사이에 단조 증가 및 단조 감소를 하는 일정한 변화 양상을 갖고 있음을 알 수 있는데, 이에 따라 본 발명에서는 이러한 특성을 이용하여 외기온도와 상대습도 및 일사량의 기상데이터를 예측한다.
(3) 무차원값 산출단계(S300)
이 단계는 외기온도, 상대습도 및 일사량을 예측하기 위해 외기온도, 상대습도 및 일사량의 무차원값을 산출하는 단계로서 이를 위해 본 발명에서는 하루 동안의 시간별 외기온도, 상대습도 및 일사량을 무차원화시킨다.
① 무차원 외기온도(
Figure PCTKR2012009001-appb-I000005
)
먼저 무차원 외기온도를 산출하는 방법에 대해 설명하면, 각 월에 대해 하루 동안의 시간별 외기온도를 아래의 수학식 1에 대입하여 외기온도를 무차원화시키며, 이에 의해 도 5와 같은 그래프를 얻을 수 있고, 이때 최고온도는 1, 최저온도는 -1이 되며 그 이외의 외기온도는 +1∼-1의 범위 내의 무차원값을 가지게 된다.
도 5는 2003년부터 2007년까지의 대전 지역에 있어서 6월부터 9월까지의 무차원 외기온도(
Figure PCTKR2012009001-appb-I000006
)를 수학식 1에 의해 산출한 결과를 나타낸 그래프이다.
[수학식 1]
Figure PCTKR2012009001-appb-I000007
여기서
Figure PCTKR2012009001-appb-I000008
는 무차원 외기온도이며,
Figure PCTKR2012009001-appb-I000009
는 시간별 외기온도이고,
Figure PCTKR2012009001-appb-I000010
는 하루 중 최고온도이며,
Figure PCTKR2012009001-appb-I000011
는 최고온도와 최저온도의 산술 평균온도이다.
② 무차원 상대습도(
Figure PCTKR2012009001-appb-I000012
)
온도 예측과 마찬가지로 상대습도를 예측하기 위해 하루 동안 시간별 상대습도를 아래의 수학식 2에 의해 무차원화시켜 무차원 상대습도(
Figure PCTKR2012009001-appb-I000013
)를 구한다.
도 6은 2003년부터 2007년까지의 대전 지역에 있어서 6월부터 9월까지의 무차원 상대습도(
Figure PCTKR2012009001-appb-I000014
)를 수학식 2에 의해 산출한 결과를 나타낸 그래프이다.
[수학식 2]
Figure PCTKR2012009001-appb-I000015
여기서
Figure PCTKR2012009001-appb-I000016
는 무차원 상대습도,
Figure PCTKR2012009001-appb-I000017
는 시간별 상대습도,
Figure PCTKR2012009001-appb-I000018
는 하루 중 최고 상대습도,
Figure PCTKR2012009001-appb-I000019
는 최고 상대습도와 최저 상대습도의 산술 평균 상대습도이다.
③ 무차원 일사량(
Figure PCTKR2012009001-appb-I000020
)
같은 방법으로 일사량을 예측하기 위해 하루 동안 시간별 누적 일사량을 아래의 수학식 3에 의해 무차원화시켜 무차원 일사량(
Figure PCTKR2012009001-appb-I000021
)을 구한다.
도 7은 2003년부터 2007년까지 5년간의 대전 지역에 있어서 6월부터 9월까지의 무차원 일사량(
Figure PCTKR2012009001-appb-I000022
)을 수학식 3에 의해 산출한 결과를 나타낸 그래프이다.
[수학식 3]
Figure PCTKR2012009001-appb-I000023
여기서
Figure PCTKR2012009001-appb-I000024
는 무차원 일사량,
Figure PCTKR2012009001-appb-I000025
는 시간별 누적 일사량,
Figure PCTKR2012009001-appb-I000026
는 하루 중 최대 시간별 누적 일사량이다.
(4) 상관식 결정단계(S400)
이 단계는 위의 무차원값 산출단계(S300)에서 구한 시간별 무차원값들로부터 시간과 무차원값 사이의 상관식을 구함으로써 이들 간의 상관관계를 결정하는 단계이다.
도 5 내지 도 7에서 확인할 수 있는 바와 같이 외기온도와 상대습도에 대한 월별 무차원값 곡선은 6월, 7월, 8월, 9월 모두 시간에 따라 일정한 경향을 보이고, 일사량의 경우에도 모든 달의 일출 시간은 5시로 같게 나타나면서 동시에 시간에 따라 일정한 경향을 보이며, 다만 9월의 일사량 분포는 다른 달에 비해 11시 이후의 값들이 2시간 앞으로 치우쳐 있다.
이에 따라 본 발명에서는 무차원 외기온도, 상대습도 및 일사량의 시간과의 상관관계를 이용하는데 이를 위해 무차원 외기온도(
Figure PCTKR2012009001-appb-I000027
), 무차원 상대습도(
Figure PCTKR2012009001-appb-I000028
), 무차원 일사량(
Figure PCTKR2012009001-appb-I000029
)과 시간에 대한 상관관계를 각각 아래의 수학식 4 내지 6에 의해 구한다.
[수학식 4]
Figure PCTKR2012009001-appb-I000030
여기서
Figure PCTKR2012009001-appb-I000031
는 무차원 외기온도,
Figure PCTKR2012009001-appb-I000032
는 상관계수,
Figure PCTKR2012009001-appb-I000033
는 시간이다.
[수학식 5]
Figure PCTKR2012009001-appb-I000034
여기서
Figure PCTKR2012009001-appb-I000035
는 무차원 상대습도,
Figure PCTKR2012009001-appb-I000036
는 상관계수,
Figure PCTKR2012009001-appb-I000037
는 시간이다.
[수학식 6]
Figure PCTKR2012009001-appb-I000038
여기서
Figure PCTKR2012009001-appb-I000039
는 무차원 일사량,
Figure PCTKR2012009001-appb-I000040
는 상관계수,
Figure PCTKR2012009001-appb-I000041
는 시간이다.
상기의 수학식 4 내지 6에 의해 생성된 각각의 상관식은 정확성을 높일 수 있도록 6차 다항식 형태로 계산되며, 아래의 표 1 내지 3에 나타난 상관계수(
Figure PCTKR2012009001-appb-I000042
)는 2003년부터 2007년까지 5년간의 대전 지역에 있어서 6월부터 9월까지의 각각의 상관식에 대한 상관계수를 구한 것으로서, 이와 동일한 방법에 의해 1월부터 12월까지의 각각의 상관식에 대한 상관계수도 구할 수 있다.
아래의 표 1은 외기온도에 대한 상관계수, 표 2는 상대습도에 대한 상관계수, 표 3은 일사량에 대한 상관계수이다.
표 1
6월 7월 8월 9월
Figure PCTKR2012009001-appb-I000043
-0.62 -0.66 -0.62 -0.65
Figure PCTKR2012009001-appb-I000044
-0.00582 0.1018 0.11 0.18
Figure PCTKR2012009001-appb-I000045
-0.096 -0.14 -0.15 -0.2
Figure PCTKR2012009001-appb-I000046
0.027 0.033 0.035 0.045
Figure PCTKR2012009001-appb-I000047
-0.0023 -0.002 -0.003 -0.0039
Figure PCTKR2012009001-appb-I000048
8.29E-5 1.0E-4 1.03E-4 1.4E-4
Figure PCTKR2012009001-appb-I000049
1.6E-7 -1.29E-6 -1.33E-6 -1.86E-6
표 2
6월 7월 8월 9월
Figure PCTKR2012009001-appb-I000050
0.6045 0.74028 0.735 0.76831
Figure PCTKR2012009001-appb-I000051
0.0548 -0.0073 -0.08855 -0.21334
Figure PCTKR2012009001-appb-I000052
0.0813 0.09098 0.13419 0.19735
Figure PCTKR2012009001-appb-I000053
-0.026 -0.02656 -0.03419 -0.04399
Figure PCTKR2012009001-appb-I000054
0.0024 0.00239 0.00299 0.00366
Figure PCTKR2012009001-appb-I000055
-9.27E-5 -8.71E-5 -1.09E-4 -1.30E-4
Figure PCTKR2012009001-appb-I000056
1.252E-6 1.131E-6 1.437E-6 1.682E-6
표 3
6월 7월 8월 9월
Figure PCTKR2012009001-appb-I000057
-0.02673 -0.03755 -0.03447 -0.00429
Figure PCTKR2012009001-appb-I000058
0.18488 0.2269 0.22167 0.12743
Figure PCTKR2012009001-appb-I000059
-0.13492 -0.15505 -0.15539 -0.11915
Figure PCTKR2012009001-appb-I000060
0.03178 0.03475 0.03542 0.03229
Figure PCTKR2012009001-appb-I000061
-0.00289 -0.00307 -0.00317 -0.00322
Figure PCTKR2012009001-appb-I000062
1.12E-4 1.17E-4 1.2E-4 1.35E-4
Figure PCTKR2012009001-appb-I000063
-1.57E-6 -1.63E-6 -1.72E-6 -2.02E-6
(5) 익일 시간별 기상데이터 예측단계(S500)
이 단계는 위의 상관식 결정단계(S400)에서 구한 시간별 무차원값을 아래의 수학식 7 내지 9에 각각 대입하여 익일 하루 동안의 시간별 외기온도, 상대습도 및 일사량의 변화를 각각 변화를 예측하는 단계이다.
[수학식 7]
Figure PCTKR2012009001-appb-I000064
여기서
Figure PCTKR2012009001-appb-I000065
는 익일 예측 외기온도,
Figure PCTKR2012009001-appb-I000066
는 상관식(수학식 4)으로 구한 무차원 외기온도,
Figure PCTKR2012009001-appb-I000067
,
Figure PCTKR2012009001-appb-I000068
는 각각 익일 최대온도와 익일 평균온도이다.
[수학식 8]
Figure PCTKR2012009001-appb-I000069
여기서
Figure PCTKR2012009001-appb-I000070
는 익일 예측 상대습도,
Figure PCTKR2012009001-appb-I000071
는 상관식(수학식 5)으로 구한 무차원 상대습도,
Figure PCTKR2012009001-appb-I000072
,
Figure PCTKR2012009001-appb-I000073
는 각각 익일 최대 상대습도와 익일 평균 상대습도이다.
[수학식 9]
Figure PCTKR2012009001-appb-I000074
여기서
Figure PCTKR2012009001-appb-I000075
는 익일 예측 일사량,
Figure PCTKR2012009001-appb-I000076
는 상관식(수학식 6)으로 구한 무차원 일사량,
Figure PCTKR2012009001-appb-I000077
,
Figure PCTKR2012009001-appb-I000078
는 각각 익일 최대 일사량과 평균 일사량이다.
(6) 최대,최저값 추정단계
이 단계는 위의 익일 시간별 기상데이터 예측단계(S500)에서 수학식 7 내지 9의 예측식을 이용하는 데에 필요한 최고 상대습도, 최저 상대습도, 최고 일사량 및 최저 일사량을 추정하는 단계이다.
위 수학식 7 내지 9를 이용하여 익일 시간별 외기온도, 상대습도, 일사량을 예측하기 위해서는 입력데이터인
Figure PCTKR2012009001-appb-I000079
,
Figure PCTKR2012009001-appb-I000080
,
Figure PCTKR2012009001-appb-I000081
,
Figure PCTKR2012009001-appb-I000082
,
Figure PCTKR2012009001-appb-I000083
,
Figure PCTKR2012009001-appb-I000084
등을 알아야 하는데, 이들 값 중 최고 외기온도와 최저 외기온도는 기상청의 일기예보로부터 쉽게 알 수 있는 반면, 최고 상대습도와 최저 상대습도, 최고 일사량과 최저 일사량 등은 미리 알 수 없기 때문에 추정하여야 한다.
이를 위해 본 발명에서는 "온도가 높고, 운량(雲量)이 많아지면 상대습도가 낮아진다", "온도가 높고, 운량이 적어지면 일사량이 높아진다"와 같은 온도, 운량 및 상대습도, 온도, 운량 및 일사량간의 애매한 특성을 퍼지 알고리즘(Fuzzy Algorithm)에 적용함으로써 상대습도와 일사량의 최고,최저값을 추정한다.
최고,최저 상대습도 및 일사량을 추정하기 위해 사용되는 퍼지 알고리즘의 입력변수로는 최고온도와 최저온도, 운량을 사용하는데, 여기서 운량은 기상청의 날씨 예보를 통해 퍼지(fuzzy) 정량화할 수 있으며, 본 발명에서는 날씨 예보에서 말하는 '맑음'은 구름의 양이 0∼2.5, '구름 조금'은 2.5∼5, '구름 많음'은 5∼7.5, '흐림'은 7.5∼10 범위의 값으로 정량화한다.
퍼지 알고리듬을 적용하기 위해 출력변수에 대한 멤버십은 표 4 및 5에 보인 값을 사용하고, 추론방법으로서는 Min-Max 방법을 사용하며, 역퍼지화 방법으로는 무게중심법을 사용하는데, 이들 Min-Max 방법 및 무게중심법은 이미 잘 알려져 있으므로 이에 대한 상세한 설명은 생략한다.
아래의 표 4는 상대습도를 구하기 위한 멤버십이고, 표 5는 일사량을 구하기 위한 멤버십이다.
표 4
Figure PCTKR2012009001-appb-T000001
표 5
Figure PCTKR2012009001-appb-T000002
위에서 설명한 본 발명의 익일 외기온도, 상대습도, 일사량의 예측방법의 정확성을 검증하기 위해 본 발명자들은 본 발명에 따라 계산된 익일 시간별 값들과 실제 기상청에서 측정한 값들을 비교해 보았다.
도 8은 2008년 7월 30일 대전 지역에 있어서 수학식 7을 통해 계산된 예측 외기온도와 기상청의 측정 외기온도를 시간별로 비교한 결과를 나타낸 그래프로서 본 발명의 예측 외기온도와 기상청의 실측 외기온도가 거의 일치함을 알 수 있다.
도 9는 2008년 7월 30일 대전 지역에 있어서 수학식 8을 통해 계산된 상대습도 예측값과 기상청의 실측 상대습도를 시간별로 비교한 결과를 나타낸 그래프이다. 이 그래프로부터 예측값의 시간별 분포와 실측값의 시간별 분포경향은 매우 잘 일치하고 있음을 알 수 있다. 다만 예측 상대습도가 실측 상대습도보다 전반적으로 약 5% 정도 낮은 분포를 나타내었는데, 이는 예측을 위한 최고,최저 상대습도를 퍼지 알고리즘에 의해 추정하는 과정에서 추정된 값이 실제 측정값보다 낮게 추정되었기 때문에 전체적인 예측 결과도 실측 결과보다 낮은 것으로 판단되며, 이러한 오차는 퍼지 알고리듬의 멤버 함수의 조정 등을 통해 해소할 수 있다.
도 10 역시 2008년 7월 30일 대전 지역에 있어서 수학식 9를 통해 계산된 예측값과 기상청에서 실측한 시간별 누적 일사량을 비교한 결과를 나타낸 그래프로서 예측값과 실측값이 상당한 정도로 일치되고 있음을 알 수 있다.
또한 도 11 내지 13은 각각 2008년의 7월과 8월의 2달 동안에 대해 대전지역의 예측된 외기온도, 상대습도 및 일사량과 기상청에서 실제로 측정한 데이터와 비교한 것으로 이들 그래프로부터 예측값과 실측값이 비교적 일치하고 있음을 알 수 있다.
이상 설명한 바와 같이 본 발명은 실측값을 사용하지 않고도 기상청에서 제공하는 최고 온도와 최저 온도만으로도 시간별 외기온도, 상대습도 및 일사량을 정확하게 예측할 수 있으며, 이를 통하여 익일의 냉난방 부하를 더욱 정확하게 예측할 수 있다.

Claims (3)

  1. 기상청으로부터 과거의 기상데이터를 확보하는 기상데이터 획득단계(S100)와;
    상기 기상데이터 획득단계(S100)로부터 획득된 기상데이터를 분석하여 필요한 데이터를 추출하는 기상데이터 분석 및 추출단계(S200)와;
    상기 기상데이터 분석 및 추출단계(S200)에서 추출된 기상데이터를 무차원화시켜 무차원값을 산출하는 무차원값 산출단계(S300)와;
    상기 무차원값 산출단계(S300)에 의해 산출된 무차원값으로부터 상관관계를 나타내는 상관식 결정단계(S400) 및;
    상기 상관식 결정단계(S400)에서 구한 시간별 무차원값으로부터 익일 시간별 기상데이터를 예측하는 익일 시간별 기상데이터 예측단계(S500)로 이루어지되,
    상기 익일 시간별 기상데이터 예측단계(S500)에서 사용되는 최고, 최저 상대습도 및 일사량은 퍼지 알고리듬에 의해 추정하는 것을 특징으로 하는 냉난방부하 추정을 위한 시간별 기상데이터 예측방법.
  2. 청구항 1에 있어서,
    상기 무차원값 산출단계(S300)에서 산출되는 무차원값은 무차원 외기온도(
    Figure PCTKR2012009001-appb-I000085
    ), 무차원 상대습도(
    Figure PCTKR2012009001-appb-I000086
    ), 무차원 일사량(
    Figure PCTKR2012009001-appb-I000087
    )로서 아래의 수학식 1 내지 3에 의해 각각 구하는 것을 특징으로 하는 냉난방부하 추정을 위한 시간별 기상데이터 예측방법.
    [수학식 1]
    Figure PCTKR2012009001-appb-I000088
    여기서
    Figure PCTKR2012009001-appb-I000089
    는 무차원 외기온도,
    Figure PCTKR2012009001-appb-I000090
    는 시간별 외기온도,
    Figure PCTKR2012009001-appb-I000091
    는 하루 중 최고온도,
    Figure PCTKR2012009001-appb-I000092
    는 최고온도와 최저온도의 산술 평균온도이다.
    [수학식 2]
    Figure PCTKR2012009001-appb-I000093
    여기서
    Figure PCTKR2012009001-appb-I000094
    는 무차원 상대습도,
    Figure PCTKR2012009001-appb-I000095
    는 시간별 상대습도,
    Figure PCTKR2012009001-appb-I000096
    는 하루 중 최고 상대습도,
    Figure PCTKR2012009001-appb-I000097
    는 최고 상대습도와 최저 상대습도의 산술 평균 상대습도이다.
    [수학식 3]
    Figure PCTKR2012009001-appb-I000098
    여기서
    Figure PCTKR2012009001-appb-I000099
    는 무차원 일사량,
    Figure PCTKR2012009001-appb-I000100
    는 시간별 누적 일사량,
    Figure PCTKR2012009001-appb-I000101
    는 하루 중 최대 시간별 누적 일사량이다.
  3. 청구항 2에 있어서,
    상기 상관식 결정단계(S400)에서 사용되는 무차원 외기온도(
    Figure PCTKR2012009001-appb-I000102
    ), 무차원 상대습도(
    Figure PCTKR2012009001-appb-I000103
    ), 무차원 일사량(
    Figure PCTKR2012009001-appb-I000104
    )에 대한 상관식은 각각 아래의 수학식 4 내지 6이 사용되는 것을 특징으로 하는 냉난방부하 추정을 위한 시간별 기상데이터 예측방법.
    [수학식 4]
    Figure PCTKR2012009001-appb-I000105
    여기서
    Figure PCTKR2012009001-appb-I000106
    는 무차원 외기온도,
    Figure PCTKR2012009001-appb-I000107
    는 상관계수,
    Figure PCTKR2012009001-appb-I000108
    는 시간이다.
    [수학식 5]
    Figure PCTKR2012009001-appb-I000109
    여기서
    Figure PCTKR2012009001-appb-I000110
    는 무차원 상대습도,
    Figure PCTKR2012009001-appb-I000111
    는 상관계수,
    Figure PCTKR2012009001-appb-I000112
    는 시간이다.
    [수학식 6]
    Figure PCTKR2012009001-appb-I000113
    여기서
    Figure PCTKR2012009001-appb-I000114
    는 무차원 일사량,
    Figure PCTKR2012009001-appb-I000115
    는 상관계수,
    Figure PCTKR2012009001-appb-I000116
    는 시간이다.
PCT/KR2012/009001 2011-12-29 2012-10-30 냉난방부하 추정을 위한 시간별 기상데이터 예측방법 WO2013100345A1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US14/369,639 US20140365128A1 (en) 2011-12-29 2012-10-30 Method for predicting hourly climatic data to estimate cooling/heating load

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
KR1020110146439A KR101141027B1 (ko) 2011-12-29 2011-12-29 냉난방부하 추정을 위한 시간별 기상데이터 예측방법
KR10-2011-0146439 2011-12-29

Publications (1)

Publication Number Publication Date
WO2013100345A1 true WO2013100345A1 (ko) 2013-07-04

Family

ID=46271278

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/KR2012/009001 WO2013100345A1 (ko) 2011-12-29 2012-10-30 냉난방부하 추정을 위한 시간별 기상데이터 예측방법

Country Status (3)

Country Link
US (1) US20140365128A1 (ko)
KR (1) KR101141027B1 (ko)
WO (1) WO2013100345A1 (ko)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150073737A1 (en) * 2012-04-20 2015-03-12 Tatsuki Inuzuka Power monitoring apparatus and power monitoring method
CN104236620B (zh) * 2013-06-14 2018-06-15 富泰华工业(深圳)有限公司 具有温湿度传感器的电子装置及湿度显示方法
KR101515003B1 (ko) * 2014-12-31 2015-04-24 (주)가교테크 일사량 예측방법
JP6467953B2 (ja) * 2015-01-30 2019-02-13 中国電力株式会社 気温予測システム、気温予測方法およびプログラム
WO2018111071A1 (es) * 2016-12-16 2018-06-21 Aguilar Fraga Omar Algoritmo regional para la automatización del cálculo de las temperaturas mínimas en zonas agrícolas utilizando factores difusos
KR101928516B1 (ko) 2017-11-13 2018-12-12 현대건설 주식회사 Rts-sarek 프로그램을 이용한 시각별 냉난방 부하 도출 방법
CN110738380B (zh) * 2018-07-18 2023-11-07 浙江盾安节能科技有限公司 热负荷控制方法、装置及系统
KR102035818B1 (ko) 2019-01-07 2019-10-24 주식회사 나라컨트롤 클라우드 기반 건물 냉난방 제어로직 자동 생성방법
CN110390470A (zh) * 2019-07-01 2019-10-29 北京工业大学 建筑气候分区方法及设备
CN113033016B (zh) * 2021-04-12 2021-10-22 北京信息科技大学 基于混合式的中央空调负荷预测方法和装置
CN113551375B (zh) * 2021-07-23 2022-09-06 贵州汇通华城股份有限公司 一种用于电子工厂的冷负荷预测方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100327995B1 (ko) * 1999-09-17 2002-03-09 양해원 심야 전력을 이용하는 축열식 난방 기기의 난방 부하량 예측 제어기 및 그 방법
KR100753141B1 (ko) * 2007-01-22 2007-08-30 충남대학교산학협력단 냉난방부하 추정을 위한 온도와 습도 예측방법

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6307500B1 (en) * 1999-08-13 2001-10-23 University Corporation For Atmospheric Research Method of moment estimation and feature extraction for devices which measure spectra as a function of range or time
US6748327B1 (en) * 2001-10-11 2004-06-08 Larry Watson Method for spatial estimation of solar radiation flux density and other solar correlated meteorological variables
IE20070331A1 (en) * 2006-05-03 2007-11-14 Lightwave Technologies Ltd A method of optimising energy consumption
KR100830095B1 (ko) * 2007-11-12 2008-05-20 충남대학교산학협력단 냉방부하 예측방법
EP2356496A4 (en) * 2008-11-13 2015-05-27 Univ Saint Louis APPARATUS AND METHOD FOR COMMUNICATING ENVIRONMENTAL PREDICTION INDICATORS TO EMERGENCY RESPONSE MANAGERS

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100327995B1 (ko) * 1999-09-17 2002-03-09 양해원 심야 전력을 이용하는 축열식 난방 기기의 난방 부하량 예측 제어기 및 그 방법
KR100753141B1 (ko) * 2007-01-22 2007-08-30 충남대학교산학협력단 냉난방부하 추정을 위한 온도와 습도 예측방법

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DO YOUNG, HAN ET AL.: "Real - Time Building Load Prediction by the On - Line Weighted Recursive Least Square Method", KOREAN JOURNAL OF AIR-CONDITIONING AND REFRIGERATION ENGINEERING, vol. 12, no. 6, 2000, pages 609 - 615 *
TAEK-BEOM, KOH: "The Study on Cooling Load Forecast of Ice-Storage System using Neural Network", JOURNAL OF KOREA INDUSTRIAL INFORMATION SYSTEMS SOCIETY, 2006, pages 115 - 118 *

Also Published As

Publication number Publication date
US20140365128A1 (en) 2014-12-11
KR101141027B1 (ko) 2012-05-03

Similar Documents

Publication Publication Date Title
WO2013100345A1 (ko) 냉난방부하 추정을 위한 시간별 기상데이터 예측방법
WO2016108420A1 (ko) 일사량 예측방법
WO2016114477A1 (ko) 예측 일사량을 이용한 냉난방부하 예측방법
WO2017010778A1 (ko) 대형 운용 시스템의 상태 기반 예방정비 장치 및 방법
WO2020096141A1 (ko) 초단기예보모델 기반의 레이더 융해층 고도 자료 동화 방법, 이를 수행하기 위한 기록 매체 및 장치
WO2014104574A1 (ko) 선형배열 영상 센서와 자세제어 센서 간의 절대 오정렬 보정방법
WO2010041836A2 (en) Method of detecting skin-colored area using variable skin color model
WO2014073884A1 (ko) 전산유체역학을 이용한 대기유동장 수치모의 시스템 및 방법
WO2021210838A1 (ko) 다양한 오믹스 데이터 분석에 기초한 생체나이 예측방법 및 시스템
WO2014119817A1 (ko) 계장화 압입 시험법을 이용한 잔류응력 평가 방법과 이를 포함하는 컴퓨터 프로그램이 기록된 기록매체 및 이 기록매체가 실행되어 계장화 압입 시험을 수행하는 압입 시험장치
WO2019066139A1 (ko) 지상 기반 라이다, 라이다 측정오차 보정 장치 및 방법
Srivastava et al. Estimation of global solar radiation in Uttar Pradesh (India) and comparison of some existing correlations
Chowdhury et al. Prediction and comparison of monthly indoor heat stress (WBGT and PHS) for RMG production spaces in Dhaka, Bangladesh
WO2017095101A1 (ko) 건물 에너지 성능 측정 평가 시스템 및 구동 방법
WO2013089535A1 (ko) 다중선형 회귀분석을 이용한 지상라이다 자료 필터링 방법
CN110631710A (zh) 一种无挡片红外测温方法及装置
CN110781611B (zh) 基于Lambert W函数的光伏组件温度实时计算方法
Piliougine et al. Series resistance temperature sensitivity in degraded mono–crystalline silicon modules
WO2018199656A1 (ko) 변전소의 자산 관리 방법
WO2017002999A1 (ko) 복수의 온도센서가 구비된 라디오존데와 이를 이용한 온도 측정 방법 및 그 보정 시스템과 보정 방법
CN102621180B (zh) 一种节能门窗性能测试方法
WO2014035001A1 (ko) 전력 관리 방법 및 시스템
CN113865750A (zh) 非接触式设备的测温校准方法及非接触式设备
WO2009154586A1 (en) Process for prediction of parameter sensor response
Steeley et al. Ambient temperature corrected dynamic transmission line ratings at two PG&E locations

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 12863839

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 14369639

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 12863839

Country of ref document: EP

Kind code of ref document: A1