The Korean Society For New And Renewable Energy

[ Article ]

New & Renewable Energy - Vol. 12, No. s2, pp.36-43

ISSN: 1738-3935 (Print)

Article No. [2016-10-PV-005]

Print publication date Oct 2016

Received 31 May 2016 Revised 03 Aug 2016 Accepted 03 Aug 2016

DOI: https://doi.org/10.7849/ksnre.2016.10.12.s2.36

표준기상년 데이터의 법선면 직달일사량 불확도에 대한 고찰

김신영¹⁾ ; 이현진²⁾ ; 김현구¹⁾^{, *} ; 장길수³⁾ ; 윤창열¹⁾ ; 강용혁¹⁾ ; 강주식⁴⁾ ; 최종오⁴⁾

A Study on Uncertainty to Direct Normal Irradiance of Typical Meteorological Year Data

Shinyoung Kim¹⁾ ; Hyunjin Lee²⁾ ; Hyungoo Kim¹⁾^{, *} ; Gilsoo Jang³⁾ ; Changyeol Yun¹⁾ ; Yongheack Kang¹⁾ ; Chushik Kang⁴⁾ ; Jongoh Choi⁴⁾

1)Korea Institute of Energy Research
2)Kookmin University
3)Korea University
4)Korea Research Institute of Standards and Science

Correspondence to: Korea Institute of Energy Research E-mail: hyungoo@kier.re.kr Tel: +82-42-860-3376 Fax: +82-42-860-3543

Copyright © 2016 by the New & Renewable Energy
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Typical meteorological year (TMY) data consist of hourly values of solar radiation and meteorological elements over a 1 year period. This reflects the long term meteorological characteristics because it is derived from a long term meteorological database. The data were produced by applying the Sandia method and weighting factor of National Renewable Energy Laboratory (NREL) by new and renewable energy resource center in Korea Institute of Energy Research (KIER), and it was registered as the qualified reference standards in 2015. The reliability of the solar radiation is important as weighting factor of the solar radiation in the TMY data is 50%. The hourly global horizontal irradiance (GHI) among the solar radiation was provided by Korea Meteorological Administration (KMA), and direct normal irradiance (DNI) was predicted by applying Reindl-2 model.^[3] The uncertainty of the DNI is higher than the experimental value as it is a predicted value. The GHI, extra-terrestrial radiation, and zenith angle were surmised as the principal factors of uncertainty. The uncertainty on the zenith angle was highest among the 3 factors. The procedure of the uncertainty evaluation obeyed the guide to the expression of uncertainty in the measurement (GUM).^[8] The uncertainty on the DNI of TMY data in Daejeon was calculated and analyzed in this paper.

Keywords:

Typical Meteorological Year Data, Direct Normal Irradiance, Zenith Angle, Uncertainty, Reference Standards

키워드:

표준기상년 데이터, 법선면 직달일사량, 천정각, 불확도, 참조표준

1. 서 론

신재생에너지를 활용한 시스템의 발전량 및 효율에 대한 사전 평가를 하기 위해서는 기본적으로 신뢰성 있는 표준기상년 데이터가 필요하다. 표준기상년 데이터는 장기간 축적된 기상데이터베이스를 바탕으로 각 월마다 기상의 장기 특성을 반영하는 대표월을 산정하여 제작한 시간 간격의 연간 기상데이터이다. 기본적으로 일사량과 건구온도, 습구온도, 풍속 등 다양한 기상요소를 포함한다. 본 논문에서 쓰인 표준기상년 데이터는 한국에너지기술연구원의 데이터^[1]로, 2015년 유효참조표준데이터로 등록된 데이터이다. 1991년 1월 1일서부터 2010년 12월 31일까지 20년간의 장기 기상 데이터베이스를 바탕으로, Sandia method (1978)^[2]를 이용하고, NREL의 가중치^[2]를 적용하여 제작되었다. Table 1의 기상 요소별 가중치를 보면 일사량의 가중치가 50% 로 다른 요소들에 비해 비중이 높음을 알 수 있다. 일사량은 기상요소들 중에서 차지하는 비중이 가장 크므로 표준기상년 데이터가 신재생에너지 중 일사량의 영향을 많이 받는 태양광, 태양열 등 태양에너지 분석에 적합한 데이터임을 유추할 수 있다. 수평면 전일사량 및 건구온도, 습구온도 등의 데이터는 기상청으로부터 제공받은 측정 데이터를 이용하였으며, 법선면 직달일사량의 경우, Reindl-2 model(1990)^[3]를 이용하여 예측하였다. Liu & Jordan (1960)^[4]을 비롯한 다양한 성분별 관계식이 있지만^[5] 한국에너지기술연구원에서 2007~2009년 측정한 대전의 수평면 전일사량 데이터를 이용한 검증 결과 Reindl-2 model^[3]의 RMSE가 70.29W/m²로 가장 정확한 것으로 밝혀졌다.^[6] Reindl-2 model^[3]에서는 k_t의 구간에 따른 k_t와 k_d의 관계식을 통해 I_d를 구하고, 이를 이용하여 법선면 직달일사량을 구하기 때문에 k_t와 k_d의 상관관계가 중요하다. ^[3]에 의하면 미국과 유럽 내 총 6개 지점의 각각 1년 이상의 측정 기간 동안 수평면 전일사량, 산란 일사량 혹은 법선면 직달일사량, 주위온도, 습구온도 혹은 이슬점 온도의 실측 데이터를 바탕으로 k_d에 영향을 미치는 주요 변수를 선별하기 위해 단계적 회귀분석을 했다. 이를 통해 회귀식의 변곡점이 생기는 지점을 기준으로 k_t의 구간을 세 구간으로 나누었으며, k_d에 영향을 미치는 주요 변수가 k_t와 태양 고도임을 밝혀냈다. 또한 본 모델식은 지역적, 계절적인 영향을 받을 수 있지만 이러한 요인들은 k_d와 k_t, 그리고 태양 고도의 상관성에 매우 낮은 영향을 미칠 것이라고 밝히고 있다. 이에 본 논문의 대상 지점인 우리나라 대전 지점에도 동일한 세 구간으로 나누어 Reindl-2 model을 적용했으며, 법선면 직달일사량 실측 데이터와 예측 데이터의 비교를 통해 회귀식의 설명력 R²가 0.937로 그 상관성이 매우 높음을 확인했다.

Table 1.

Weighting values for FS statistics[2]

법선면 직달일사량은 예측한 값이기에 실측값보다 더 많은 불확실성이 존재한다. 특히 태양열을 활용한 건물 냉난방 설비의 성능을 평가하기 위해 표준기상년 데이터를 활용할 때 태양열 시스템에 영향을 미치는 주요소는 법선면 직달일사량이기에 이에 대해 신뢰할 수 있는 범위를 예측하는 것은 필수적인 일이다. 하지만 국내에서 배포하고 있는 표준기상년 데이터 중 대표적인 태양에너지학회의 대한민국표준기상데이터^[7]를 보면 데이터를 신뢰할 수 있는 범위가 기재되어 있지 않다. 또한, 한국에너지기술연구원의 표준기상년 데이터^[1]의 법선면 직달일사량과 비교를 해보면 태양에너지학회의 법선면 직달일사량은 한국에너지기술연구원의 법선면 직달일사량에 비해 전반적으로 낮은 경향을 보이며, RMSE 196.76W/m²의 차이가 있다. 이에 본 논문에서는 한국에너지기술연구원 표준기상년 데이터^[1]의 법선면 직달일사량에 대한 신뢰성을 검증하고 신뢰할 수 있는 범위를 알고자 불확도의 요인을 분석, 산출하고 그 결과를 분석했다. 불확도 평가 절차는 GUM^[8]을 준수하였으며, 국내에서 일사량이 비교적 높은 대전 지점(위도: 36°37′, 경도: 127°37′)을^[9] 대상 지점으로 선정하였다. 불확도 평가 결과를 다양한 관점에서 해석해보기 위해, 법선면 직달일사량의 증가에 따른 불확도와 불확도 비의 변화, 그리고 k_t의 구간에 따른 k_d와 법선면 직달일사량, 불확도의 변화 양상을 분석했다. 또한 이러한 변화 양상에 대하여 계절적 차이를 알아보기 위해 계절에 따른 결과를 비교, 분석했다.

2. 방법론

2.1 표준기상년 데이터

2.1.1 Sandia Method^[2]

1) Sandia method 제작 단계

① FS 통계(1971)^[10]를 이용하여 각 기상요소 별로 장, 단기 누적분포함수를 계산한다. 가중치 w_i를 고려한 10개 기상요소의 FS 합 기준으로 가장 낮은 값을 갖는 5개 후보월을 선정한다.

F S = 1 n F S ∑ i = 1 n F S δ i

(1)

W S = ∑ w i F S i

(2)

② 장기 평균값과 중앙값에 대한 근접성을 기준으로 5개 후보월에 등급을 부여한다.

③ 평균 건구온도와 수평면 전일사량의 장기 데이터를 바탕으로 지속성이 없는 후보월은 제외하고 2단계의 높은 등급을 가진 월을 대표월로 선정한다.

④ 선택한 12개 대표 월을 연결하여 표준기상년 데이터를 제작한다. 연결 시에, 월간 접점에서 발생하는 불연속성을 제거하기 위해 6시간 데이터는 선형 보간 하여 사용한다.

2.1.2 Reindl-2 model^[3,5,11]

시간별 법선면 직달일사량 I_n을 구하기 위해 먼저 실측한 시간별 수평면 전일사량 I_t를 이용하여 k_t를 구한다.^[5]

k t = I t I o

(3)

여기서 수평면 대기권 밖 일사량 I_o는 다음과 같다.^[11]

I o = I on cos θ z = I sc 1.00011 + 0.034221 cos B + 0.00128 sin B + 0.000719 cos 2 B + 0.000077 sin 2 B cos θ z

(4)

cos θ z = cos ϕ cos δ cos ω + sin ϕ sin δ

(5)

B = n - 1 360365

(6)

k_t의 구간에 따른 k_d는 다음과 같다.^[3,5]

k d = I d I t

(7)

I d I t = 1.02 - 0.254 k t + 0.0123 cos θ z k t ≤ 0.3

(8)

$I d I t = 1.4 - 1.749 k t + 0.177 cos θ z 0.3 < k t < 0.78 = 0.486 k t - 0.182 cos θ z k t ≥ 0.78$

최종적으로 시간별 법선면 직달일사량 I_n은 다음과 같다.

I n = I t 1 - k d / cos θ z

(9)

2.2 불확도 평가 절차^[3,5,8]

법선면 직달일사량을 구하는 수식 (9)에서 측정값은 시간별 수평면 전일사량 I_t이며, 계산값은 직달성분의 대기권 밖 일사량 I_on, 코사인 천정각 cosθ_z이다. 이 세 가지 요소는 법선면 직달일사량의 주요한 불확도 요인으로 추정된다. 이에 시간별 법선면 직달일사량 I_n에 대한 불확도 평가 절차는 다음과 같다.

2.2.1 모델식 수립

시간별 법선면 직달일사량의 모델식은 다음과 같이 k_t의 범위에 따른 세 구간의 모델식으로 나뉘어진다.

1) k_t≤0.3 인 경우

y = I n = 0.254 I on cos 2 θ z I t 2 - 0.0123 + 0.02 cos θ z I t

(10)

2) 0.3<k_t<0.78 인 경우

y = I n = - 1.749 I on cos 2 θ z I t 2 - 0.177 + 0.4 cos θ z I t

(11)

3) k_t≥0.78 인 경우

y = I n = - 0.486 I on cos 2 θ z I t 2 + 0.182 + 1 cos θ z I t

(12)

2.2.2 표준불확도

법선면 직달일사량 모델식의 변수이며 주요 불확도 요인인 수평면 전일사량, 직달성분의 대기권 밖 일사량, 그리고 코사인 천정각에 대한 표준불확도를 구한다. 시간별 수평면 전일사량은 기상청에서 1분마다 수신된 누적 일사량이며, 측정값의 통계적인 분석 즉 우연효과에 의한 불확도와 일사계의 교정불확도가 존재한다. 직달성분의 대기권 밖 일사량은 수식 자체의 불확도^[11]를 반영한다. 태양의 위치를 예측할 때 천정각을 구하게 되는데, 예측 방법은 Astronomical Almanac algorithm(1988)^[12]을 이용했으며, 코사인 천정각에 의한 표준불확도는 이 algorithm의 정확도를 반영한다.

1) I_t에 대한 표준불확도

u 2 I t = u 2 I t, 우 + u 2 I t, 교 = s n 2 + U I t 2 2 = s n 2 + 0.03 ⋅ I t 2 2

(13)

2) I_on에 대한 표준불확도

u I on = U I on 3 = 0.0001 ⋅ I on 3

(14)

3) cosθ_z에 대한 불확도 기여량 : [13]에 의해 θ_z=87°일 때, cosθ_z로 인해 발생하는 I_n의 불확도가 0.35%라고 가정하였다.

u cos θ z = U cos θ z 3 = 0.0035 ⋅ I n 3

(15)

2.2.3 합성표준불확도^[3,5,8]

합성표준불확도는 앞서 구한 표준불확도에 각 불확도 요인에 대한 감도계수 c_i를 곱하여 구한다.

u c 2 I n = c 12 u 2 I t + c 22 u 2 I on + u 2 cos θ z = c 12 u 2 I t, 우 + u 2 I t, 교 + c 22 u 2 I on + u cos θ z 2

(16)

1) k_t≤0.3 인 경우

c 1 = ∂ f ∂ I t = 0.508 I on cos 2 θ z I t - 0.0123 + 0.02 cos θ z

(17)

c 2 = ∂ f ∂ I on = 0.254 I t 2 cos 2 θ z I on 2

(18)

2) 0.3<k_t<0.78 인 경우

c 1 = ∂ f ∂ I t = 3.498 I on cos 2 θ z I t - 0.177 + 0.4 cos θ z

(19)

c 2 = ∂ f ∂ I on = - 1.749 I t 2 cos 2 θ z I on 2

(20)

3) k_t≥0.78 인 경우

c 1 = ∂ f ∂ I t = 0.972 I on cos 2 θ z I t + 0.182 + 1 cos θ z

(21)

c 2 = ∂ f ∂ I on = 0.486 I t 2 cos 2 θ z I on 2

(22)

따라서 합성표준불확도는 다음과 같다.

u c I n = c 12 u 2 I t, 우 + u 2 I t, 교 + c 22 u 2 I on + u cos θ z 2

(23)

2.2.4 확장불확도^[8]

확장불확도 U는 합성표준불확도 u_c(I_n)에 포함인자 k를 곱하여 구한다. 포함인자 k의 값은 y-U에서 y+U의 구간에서 요구되는 신뢰의 수준에 따라 결정되는데, 본 논문에서는 신뢰수준이 약 95%일 때를 가정하여 k=2의 값을 선택했다.

U = 2 ⋅ u c

(24)

따라서 측정결과는 Y=y±U로 표현된다.

3. 연구 결과

표준기상년 데이터의 시간별 법선면 직달일사량에 대한 불확도 평가 결과와 k_t 구간에 따른 k_d, 법선면 직달일사량, 불확도의 변화 그리고 계절별로 비교한 결과는 다음과 같다.

3.1 불확도 평가 결과

표준기상년 데이터의 시간별 법선면 직달일사량에 대한 불확도 평가 결과는 Figure 1과 같다. 시간별 불확도의 평균은 11.15W/m², 표준오차는 0.2W/m²이며, 시간별 법선면 직달일사량에 대한 불확도의 비는 평균 0.07, 이 비의 표준오차는 0.01이다. 시간별 법선면 직달일사량의 평균이 127.17W/m²이고, 표준오차가 2.44W/m²인 것에 비하면 불확도는 법선면 직달일사량에 비해 매우 낮음을 알 수 있다.

Fig. 1.

Hourly DNI of TMY data and uncertainty

Figure 2에서는 법선면 직달일사량에 따른 불확도 비의 변화를 보여준다. 법선면 직달일사량이 0W/m²에 가까울수록 법선면 직달일사량에 대한 불확도의 비가 커지는 경향을 보인다. 이는 낮은 일사량에 대한 불확실성이 크며, 그만큼 신뢰할 수 있는 수준이 낮음을 보여준다.

Fig. 2.

Relation between DNI and ratio uncertainty to DNI

법선면 직달일사량에 따른 불확도의 변화 양상은 Figure 3과 같다. 대체적으로 법선면 직달일사량이 증가할수록 불확도도 증가하는 경향을 보인다. 하지만 법선면 직달일사량이 800W/m² 내외에서 불확도가 10W/m²인 경우들도 있다.

Fig. 3.

Relation between DNI and uncertainty

이 원인을 찾아보기 위해 k_t의 구간에 따른 k_d, 법선면 직달일사량 그리고 불확도의 관계를 보니 Figure 4와 같았다. 법선면 직달일사량의 주요 변수 중 하나인 k_d는 k_t가 증가함에 따라 감소하는 경향을 보인다. 이는 하늘이 청명할수록 산란되는 일사량의 양이 줄어듬을 의미한다. k_t가 증가함에 따라 법선면 직달일사량은 증가하는 경향을 보인다. 하지만 k_t가 0.78 이상일 때부터 k_d는 불규칙하게 흩어지는데, 이에 따라 법선면 직달일사량도 산만하게 흩어짐을 볼 수 있다. 불확도도 k_t=0.78을 기준으로 끊어지는 현상을 볼 수 있다.

Fig. 4.

Relation of kt to kd and DNI and to uncertainty

특히 하늘이 청명한 조건에서 태양 고도가 낮아짐에 따라 태양복사에너지가 대기권을 투과하는 거리가 길어지고 이에 따라 산란되는 일사량이 많아지기 때문에 k_d의 값이 커진다. 태양 고도가 낮아질 때 천정각은 커지게 되는데 천정각의 경우 Figure 5와 같이 1년 동안 10°<θ_z<170°범위 내에서 연속적으로 변화한다. 하지만 불확도 평가 시 범위를 0°<θ_z<85°로 제한함으로써 불확도의 연속성이 끊어진 것으로 보인다. 천정각이 90°에 가까운 경우는 해가 뜨거나 지는 시점으로 이 때 태양빛은 낮보다 상대적으로 두꺼운 대기를 통과하게 되어 빛이 굴절되는 정도가 커지게 된다. 이에 불확도 또한 비선형적으로 커지는 현상이 발생하여 범위를 제한했다. k_t가 0.78 이상인 비율은 전체의 0.57%로 흩어지는 구간이 전체에서 차지하는 비중은 매우 낮다.

Fig. 5.

Variation of θz according to hours

3.2 계절별 변화 비교

계절별 k_t의 변화에 따른 k_d, 법선면 직달일사량, 그리고 불확도의 변화는 Figure 6~Figure 8과 같다. 전체적으로 서로 비슷한 변화 추이를 보인다. 공통적으로 k_t≥0.78에 흩어지는 경향을 보이는데, 특히 겨울에 그러한 경향이 가장 많이 보인다. 이는 사계절 중 겨울에 태양 고도가 가장 낮고, 천정각이 가장 크므로 생기는 결과로 보인다. Figure 9에서 보는 바와 같이 겨울은 천정각의 범위가 50°<θ_z<165°로, 불확도를 예측할 수 있는 범위가 적다. 반대로 여름은 사계절 중 태양 고도가 가장 높고, 천정각은 가장 낮음으로 k_t≥0.78 구간에서 비선형적인 경향이 가장 적게 나타난다.

Fig. 6.

Relation between kt and kd in four seasons

Fig. 7.

Relation between kt and DNI in four seasons

Fig. 8.

Relation between kt and uncertainty in four seasons

Fig. 9.

Variation of θz according to hours in four seasons

4. 결 론

표준기상년 데이터의 법선면 직달일사량에 대한 불확도는 평균 11.15W/m², 표준오차 0.2W/m²로 평균적으로 낮은 수치였다. 1년 중 k_t가 0.3 이하인 날은 66.95%, 0.3<k_t<0.78인 경우는 32.48%로 이 두 구간에서 k_t가 증가함에 따라 k_d는 감소했고, 법선면 직달일사량은 증가하는 경향을 보였다. 0.3<k_t<0.78 일 때가 k_t≤0.3 일 때보다 기울기의 절대값은 더 컸다. 불확도는 법선면 직달일사량이 증가할수록 증가하는 경향을 보였다. k_t가 0.78 이상인 구간은 전체의 0.57%에 불과했지만 이 구간에서는 k_d와 법선면 직달일사량, 그리고 불확도는 산만하게 흩어지는 경향을 보였다. 이는 천정각이 90°에 가까울수록 태양빛의 굴절률이 커지므로 천정각의 범위를 85°이하로 제한함으로써 얻어진 결과들로 유추된다.

계절별로 k_t의 변화에 따른 k_d, 법선면 직달일사량, 그리고 불확도의 변화를 보면 전체적인 양상은 서로 비슷한 경향을 보인다. 특히 겨울에는 태양의 고도가 사계절 중 가장 낮음으로 인해, 천정각이 가장 크므로 불확실성이 존재하는 범위가 넓어짐에 따라 k_t가 0.78 이상인 구간에서 k_d와 법선면 직달일사량, 그리고 불확도가 사계절 중 가장 많이 흩어지는 양상을 보였다.

따라서 천정각이 0°에 가까울 때, 즉 태양이 지표면 부근에 근접해있을 때, 법선면 직달일사량의 불확도에 빛의 굴절률이 미치는 영향과 더불어 본 논문에서는 법선면 직달일사량의 변수에 의한 불확도 요인만 고려했지만 앞으로 Reindl-2 model의 불확도를 반영한 불확도에 대한 연구가 추후 더 필요할 것으로 사료된다.

Nomenclature


c_i :	correlation coefficient
FS_i :	finkelstein-schafer statistic for index
I_d :	total diffuse irradiance on a horizontal surface, W/m²
I_n :	direct normal irradiance, W/m²
I_o :	extraterrestrial radiation, W/m²
I_sc :	solar constant, 1367 W/m²
I_t :	global solar irradiance, W/m²
k_d :	diffuse fraction, dimensionless
k_t :	clearness index, dimensionless
n_FS :	the number of daily readings in a month
n :	day of the year
U :	uncertainty
w_i :	weighting for index
x_i :	ith ordered sample value
δ_i :	absolute difference between the long-term CDF and the candidate month CDF at x_i
δ :	declination, degrees
φ :	latitude, degrees
θ_z :	zenith angle, degrees
ω :	hour angle, degrees
sp :	spring
su :	summer
f :	fall
w :	winter

subscript


CDF :	cumulative distribution function
DNI :	direct normal irradiance
FS :	finkelstein-schafer
GHI :	global horizontal irradiance
GUM :	guide to the expression of uncertainty in measurement
KIER :	korea institute of energy research
KMA :	korea meteorological administration
NREL :	national renewable energy laboratory
TMY :	typical meteorological year
WS :	weighted sum

Acknowledgments

본 연구는 한국에너지기술연구원의 주요사업(B6-2427) 및 산업통상자원부의 표준화기술개발사업(B6-7601)을 재원으로 수행한 연구과제의 결과입니다.

References

Korea Institute of Energy Research, (2015), “Typical meteorological data”, http://www.srd.re.kr/db/srdDataInfoView.do?scid=S06001003&sdid=3.
S. Wilcox, and W. Marion, (2008), “Users manual for TMY3 data sets”, NREL, The United States of America. [https://doi.org/10.2172/928611]
D. T. Reindl, W. A. Beckman, and J. A. Duffie, (1990), “Diffuse fraction correlations”, Solar Energy, 45(1), p1-7. [https://doi.org/10.1016/0038-092X(90)90060-P]
Benjamin Y. H. Liu, Richard C. Jordan, (1960), “The interrelationship and characteristic distribution of direct, diffuse and total solar radiation”, Solar Energy, 4(3), p1-19. [https://doi.org/10.1016/0038-092X(60)90062-1]
L.T. Wong, W.K. Chow, (2001), “Solar radiation model”, Applied Energy, 69, p191-224. [https://doi.org/10.1016/S0306-2619(01)00012-5]
Hyun Jin Lee, Shin Young Kim, Yong Heack Kang, Hyun Goo Kim, and Chang Yeol Yun, (2015), “Comparison of global horizontal and direct normal insolations in Daejeon, Korea”, ISES Solar World Congress 2015, P2-5-2.
The Korean solar energy society, Korean institute of architectural sustainable environment and bulding systems, (2015), “Typical whether data of the Republic of Korea”, http://kses.re.kr/data_06/list_hi.php.
Uncertainty committee of Korea research institute of standards and science, (2010), “ISO/IEC Guide 98-3: 2008, Guide to the expression of uncertainty in measurement (GUM: 1995)”, Korea research institute of standards and science, Republic of Korea.
Yong Heack Kang, Hyun Goo Kim, Deok Gi Cho, Chang Yeol Yun, Wan Soon Park, Chul Hyung Lee, Joon Pyo Lee, Hyoung Chan Kim, Euy Jun Lee, Eun Cheol Kang, (2015), “Korea renewable energy resource atlas”, Ministry of Science, ICT and Future Planning, Korea Institute of Energy Research, Republic of Korea.
J. M. Finkelstein, R. E. Schafer, (1971), “Improved goodnessof-fit tests”, Biometrika, 58(3), p641-645. [https://doi.org/10.1093/biomet/58.3.641]
John A. Duffie, William A. Beckman, (2013), “Solar engineering of thermal processes”, 4th ed., Wiley, The United States of America. [https://doi.org/10.1002/9781118671603]
Joseph J. Michalsky, (1988), “The Astronomical Almanac’s algorithm for approximate solar position (1950-2050)”, Solar Energy, 40(3), p227-235. [https://doi.org/10.1016/0038-092X(88)90045-X]
Ibrahim Reda, and Afshin Andreas, (2008), “Solar position algorithm for solar radiation applications”, National renewable energy laboratory(NREL), The United States of America.
Joon Bum Jee, Seung Woo Lee, Young Jean Choi, and Kyu Tae Lee, (2012), “The Generation of Typical Meteorological Year for Research of the Solar Energy on the Korean Peninsula”, New & Renewable Energy, 8(2), p14-23. [https://doi.org/10.7849/ksnre.2012.8.2.014]
Shin Young Kim, Hyun Jin Lee, Chang Yeol Yun, Jong Kyu Kim, Yong Heack Kang, and Dok ki Jo, (2014), “Generation of Typical Solar Radiation Year for Daejeon”, Journal of the Korean Society for New and Renewable Energy, p129.
Hyun Goo Kim, Ji Hwee Choi, Moon Seok Jang, and Wan Ho Jeon, (2010), “Uncertainty Analysis on Vertical Wind Profile Measurement of LIDAR for Wind Resource Assessment”, Journal of the Korean Society for New and Renewable Energy, p185.
Hyun Goo Kim, Yu Mi Kim, Shin Young Kim, Yong Heack Kang, and Chang Yeol Yun, (2015), “Calculation of Wind Shear Exponents in Jejudo and Uncertainty Evaluation”, Journal of the Korean Society for New and Renewable Energy, p82.

기상 요소	Sandia Method	NREL TMY
수평면 전일사량	12/24	5/20
법선면 직달일사량	0	5/20
건구온도 최대값	1/24	1/20
건구온도 최소값	1/24	1/20
건구온도 평균값	2/24	2/20
습구온도 최대값	1/24	1/20
습구온도 최소값	1/24	1/20
습구온도 평균값	2/24	2/20
풍속 최대값	2/24	1/20
풍속 평균값	2/24	1/20