У меня 24 значения GHI и температуры 2 м для данной фотоэлектрической станции, расположенной в
ilat = -31,389
Я использую пакет SolaR для получения выработки электроэнергии для подключенной сетис использованием фиксированных панелей и системы отслеживания одной оси.
1) Первое, что я применил, было преобразование моего набора данных с созданием следующего data.frame, называемого P1, где G0 - GHI, а Ta - температура на 2м.
time G0 Ta
1980-01-01 00:00:00 0,0 20,7
1980-01-01 01:00:00 0,0 19,8
1980-01-0102:00:00 0,0 19,0
1980-01-01 03:00:00 0,0 18,3
1980-01-01 04:00:00 0,0 17,4
1980-01-01 05:00:00 0,0 16,8
1980-01-01 06:00:00 0,0 16,0
1980-01-01 07:00:00 0,0 15,4
1980-01-01 08:00:00 0,0 14,9
1980-01-01 09:00:00 10,8 14,5
1980-01-01 10:00:00 166,915.0
1980-01-01 11:00:00 400,6 19,7
1980-01-01 12:00:00 632,8 24,2
1980-01-01 13:00:00 840,0 26,9
1980-01-01 14:00:00 1007,8 28,9
1980-01-01 15:00:00 1113,0 30,4
1980-01-01 16:00:00 1152,5 31,4
1980-01-01 17:00:00 1123,0 32,2
1980-01-01 18:00:00 1027,5 32,8
1980-01-01 19:00: 00 872,5 32,9
1980-01-01 20:00:00 670,5 32,5
1980-01-01 21:00:00 435,2 31,4
1980-01-01 22:00:00 199,1 29,9
1980-01-01 23:00:00 21,7 27,8
2) После этого я преобразую класс своего объекта в объект "Meteo", используяФункция dfI2Meteo пакета solaR выполняется следующим образом:
облучение <- dfI2Meteo (P1, time.col = "time", lat = ilat, source = "P1", format = '% Y /% m /% d% H:% M:% S ') </p>
3) Затем выходная мощность для фиксированной системы по умолчанию рассчитывается с использованием функции prodGCPV пакета solaR следующим образом:
p.fixed <- prodGCPV (ilat, modeRad = "bdI", dataRad = облучение, образец = "час") </p>
4) Когда я ищу почасовую выработку электроэнергии с использованием набора данных P1, я получаю следующееПоведение:
Tc Pac Pdc EffI
1980-01-01 00:00:00 20,70000 0,000 0,000 0,0000000
1980-01-01 01:00:00 19,80000 0,000 0,000 0,0000000
1980-01-01 02:00:00 19,00000 0,000 0,000 0,0000000
1980-01-01 03:00:00 18,30000 0,000 0,000 0,0000000
1980-01-01 04:00: 00 17,40000 0,000 0,000 0,0000000
1980-01-01 05:00:00 16,80000 0,000 0,000 0,0000000
1980-01-01 06:00:00 NaN NaN NaN NaN
1980-01-01 07:00:00 NaN NaN NaN NaN
1980-01-01 08:00:00 NaN NaN NaN NaN
1980-01-01 09:00: 00 14,82380 0,000 0,000 0,0000000
1980-01-01 10:00:00 20,00939 3277,467 3735,144 0,9043499
1980-01-01 11:00:00 31,77288 7664,552 8471,650 0,9324471
1980-01-01 12:00:00 43.50324 11526,131 12708,465 0,9347510
1980-01-01 13:00:00 53,19436 14825,572 16378,712 0,9329045
1980-01-01 14:00: 00 60.75024 17127.105 18966.287 0.9306945
1980-01-01 15:00:00 NA NA NA NA
1980-01-01 16:00:00 NA NA NA NA
1980-01-01 17:00:00 NA NA NA NA
1980-01-01 18:00:00 NA NA NA NA
1980-01-01 19:00:00 NA NA NA NA
1980-01-01 20:00:00 32,50000 0,000 0,0000,0000000
1980-01-01 21:00:00 31,40000 0,000 0,000 0,0000000
1980-01-01 22:00:00 29,90000 0,000 0,000 0,0000000
1980-01-01 23:00:00 27,80000 0,000 0,000 0,0000000
В заключение я получаю некоторые значения «NaN» в часы, предшествующие восходу солнца, и значения NA после нашего местного полудня (15 UTC), где мощность PV должна быть важной.
Может кто-нибудь помочьмне понять, что не так с процедурой?
Большое спасибо!
Кристиан