В статье

Цель работы

Метод определения плотности солнечного излучения на горизонтальной поверхности при чистом небе на территории Вьетнама

Результаты расчёта по методике

Расчёт плотности солнечного излучения, поступающего на ориентированные поверхности

Расчёт энергетического потенциала

Заключение

В настоящее время традиционные ископаемые энергоресурсы — уголь, нефть и природный газ — остаются основным источником производства электрической энергии во многих странах мира, включая Республику Вьетнам [2]. Однако их использование сопровождается значительными выбросами парниковых газов, что является одной из главных причин изменения климата и загрязнения атмосферы. Кроме того, запасы ископаемых ресурсов постепенно истощаются, тогда как во Вьетнаме уже наблюдаются серьёзные трудности в обеспечении растущих энергетических потребностей.

Согласно прогнозам правительства Вьетнама, в 2025 году ожидается рост спроса на электроэнергию примерно на 10 % по сравнению с предыдущим годом, что обусловлено продолжением экономического роста и индустриализации страны [3]. В таких условиях развитие и внедрение возобновляемых источников энергии, особенно ветровой и солнечной, становится приоритетным направлением энергетической политики и научных исследований Вьетнама.

Республика Вьетнам расположена между 8° и 23° северной широты [4], находится вблизи экватора и получает значительное количество солнечного излучения — в среднем от 4 до 5,5 кВт·ч/м2 в день в различных регионах со средним количеством солнечных часов от 2000 до 2500 часов в год [5, 6] (рис. 1). Особенно благоприятные условия наблюдаются в центральной и южной частях страны.

Карта Вьетнама с цветовой заливкой от зелёного на севере до тёмно-красного на юге и шкалой сумм радиации
Рис. 1. Распределение интенсивности солнечной радиации на территории Вьетнама [6]

Это делает Вьетнам одним из государств в Юго-Восточной Азии с наибольшим потенциалом для развития возобновляемой энергетики, в частности, солнечной [5]. К концу 2022 года установленная мощность солнечных электростанций во Вьетнаме превысила 16 500 МВт, что составляет около 25 % общей энергетической мощности страны [1].

Провинция Нинь Тхуан, расположенная в центрально-южной части Вьетнама [7], является одним из регионов с наибольшим уровнем солнечной активности, что делает её приоритетной зоной для развития фотоэлектрических (PV) систем. Поэтому оценка потенциала солнечной энергии, как на национальном, так и на региональном уровне (на примере провинции Нинь Тхуан), является крайне важной задачей. Она позволяет определить участки с наилучшими условиями солнечной радиации, что способствует рациональному планированию, выбору оптимальных площадок для размещения солнечных электростанций, повышению инвестиционной эффективности, снижению эксплуатационных затрат и оптимизации генерации электроэнергии. Кроме того, результаты таких исследований могут служить научной основой для разработки стратегий развития возобновляемой энергетики, обеспечивающих баланс между экономическим ростом и охраной окружающей среды, а также способствовать привлечению внутренних и иностранных инвестиций, особенно в сфере «зелёных» технологий.

Одним из ключевых этапов оценки эффективности солнечных энергетических систем является определение потенциального количества солнечной энергии, доступной на единицу площади в течение суток на горизонтальных и наклонных поверхностях.

В данном исследовании используется метод Берда (Bird), применяемый для расчёта солнечной радиации на горизонтальной поверхности, как в условиях ясного неба, так и при реальных атмосферных условиях, а также методика, предложенная Лю и Джорданом, для расчёта радиации на наклонных поверхностях с целью определения оптимального угла наклона солнечной батареи для территории Вьетнама [8, 9].

Цель работы

Цель исследования заключается в разработке и адаптации методики расчёта прихода солнечной энергии на горизонтальные и наклонные поверхности для различных временных промежутков (суточного, месячного и среднегодового), а также в создании атласов среднегодового прихода солнечной радиации на горизонтальные поверхности территории Вьетнама с использованием математических моделей прихода солнечной радиации, метеорологических данных базы NASA [10].

Метод определения плотности солнечного излучения на горизонтальной поверхности при чистом небе на территории Вьетнама

В этой работе использована модифицированная методика Берда. Эта методика применяется для оценки поступления солнечной радиации на горизонтальную поверхность в условиях идеально прозрачной атмосферы, без учёта облачности и атмосферных загрязнений. Методика Берда характеризуется высокой точностью при расчётах внеземного и прямого солнечного излучения и позволяет учитывать астрономические параметры (угол склонения Солнца, широта местности, часовой угол) в течение светового дня (рис. 2).

Схема небесной сферы с солнечными лучами, экватором и обозначенными углами широты места, часового угла и склонения Солнца
Рис. 2. Схема для определения углов: широты места, часового угла и склонения Солнца

Именно благодаря своей простоте реализации, стабильности результатов и высокой корректности при отсутствии атмосферы методика широко используется как базовая при теоретических оценках энергетического потенциала региона.

Для расчёта используются данные:

  • географическая долгота, широта местности λм = 108,8° в. д. и φм = 11,4° с. ш.;
  • атмосферное давление pвозд = 1010 гПа;
  • наблюдаемые среднемесячные дневные значения сумм прямого Eгорiпрч, диффузного Eгорiдифч СИ, поступающих на горизонтальную поверхность, кВт/м2.

Склонение и часовой угол определяются по формуле:

δn = 23,5 sin[360(n − 81)/365]; (1)
ωc = kωc (tcист − 12), (2)

где n — порядковый номер дня года, отсчитываемого от 1 января; kωc — коэффициент для перевода значения часового угла Солнца в градусы, kωc = 15°/ч; tcист — истинное солнечное время (в часах).

Косинус угла падения солнечных лучей на горизонтальную поверхность определяется:

cosθz = cosδncosφмcosωc + sinφмsinδn. (3)

Время восхода и захода Солнца и продолжительность светового дня.

В моменты восхода и захода Солнца угол падения лучей на горизонтальную поверхность равен 90°, из чего следует, что cosθz = 0. Подставляя указанные значения в выражение (3), можно получить:

tcзах,восх = ±ωc/kωc + 12, (4)

откуда длительность светового дня:

Tc = tcзах − tcвосх = (2/kωc)·arccos[−tgφм·tgδn]. (5)

Значение атмосферной массы, скорректированное по высоте местности над уровнем моря, определяется по формуле:

M*i = p / {1013[cosθzi + 0,15(93,885 − θzi)−1,25]}. (6)

Среднемесячная плотность потока прямого СИ при чистом небе для i-го дня определяется:

Eгорiпрч = Ec⊥cosθziτriτO3iτгазiτH2Oi, (7)

где τri — коэффициент пропускания, учитывающий релеевское рассеяние; τO3i, τгазi и τH2Oi — коэффициенты пропускания, учитывающие поглощение солнечного излучения озоном, газовой смесью (O2, N2, CO2) и водяным паром.

Пять колоколообразных кривых плотности солнечного излучения в течение суток для января, марта, июля, октября и декабря
Рис. 3. Плотности солнечного излучения в моменты восхода, полудня и захода Солнца на 15-е число характерных месяцев года (январь, март, июль, октябрь, декабрь)

Коэффициенты пропускания атмосферы определяются по выражениям:

Система из четырёх формул коэффициентов пропускания атмосферы для релеевского рассеяния, озона, газовой смеси и водяного пара

где dO3 — толщина озонового слоя при нормальных температуре и давлении, dO3 = 0,32 см для средних широт; dH2Oi — показатель содержания пара в вертикальном столбе атмосферы, для зимы dH2O = 1,3 г/см2, для лета dH2O = 3,0 г/см2, для погоды в Нинь Тхуан принято постоянное значение dH2O = 3,0 г/см2.

Результаты расчёта по методике

С учётом географического положения Вьетнама, расположенного между 8° и 23° с. ш. и 102° и 110° в. д., то есть вблизи экватора, страна обладает значительным преимуществом для развития солнечной энергетики. Для детального анализа суточной динамики инсоляции в условиях «ясного неба» проведены теоретические расчёты плотности солнечного излучения.

На рис. 3 представлены значения плотности солнечного излучения E [Вт/м2], рассчитанные для трёх ключевых моментов суток (восход, полдень и заход Солнца) на 15-е число характерных месяцев (январь, март, июль, октябрь и декабрь). Такая визуализация позволяет оценить сезонные колебания солнечной активности и служит основой для последующего анализа энергетического потенциала в различные периоды года.

Для верификации теоретической модели и учёта влияния реальных климатических факторов в статье проведено сопоставление расчётных значений с фактическими данными многолетнего спутникового мониторинга базы NASA Prediction of Worldwide Energy Resources (POWER). Данные извлекались по географическим координатам (широте и долготе) исследуемых точек, что позволило получить точные значения суточной инсоляции за расчётный период [10].

На рис. 4 представлена среднемесячная суммарная плотность солнечного излучения [кВт·ч/м2], падающего на горизонтальную поверхность в провинции Нинь Тхуан для каждого месяца года. Сплошная линия соответствует значениям, рассчитанным по методике Берда с использованием основных коэффициентов атмосферного пропускания: молекулярного рассеяния, озонового слоя, водяного пара, аэрозолей и др. Пунктирная линия отображает фактические данные из базы NASA POWER за 2023 год, которые учитывают влияние облачности и метеорологических условий.

График с двумя кривыми — расчётной и фактической — суммарного прихода солнечного излучения по месяцам с максимумом в марте
Рис. 4. Суммарный приход солнечного излучения на горизонтальную поверхность в провинции Нинь Тхуан для каждого месяца года

Из рис. 4 видно, что максимальные значения как по расчётной, так и по фактической кривой наблюдаются в период с апреля по август. Минимальное значение зафиксировано в декабре, что соответствует наиболее короткому световому дню и меньшей солнечной активности.

Видно, что значения, полученные по модели, превышают реальные данные на протяжении всего года. Это объясняется тем, что методика Берда ориентирована на расчёты в условиях ясного неба, без учёта облачности, а данные NASA отражают реальное климатическое воздействие на поступление солнечной радиации.

Тем не менее, качественное совпадение формы кривых подтверждает адекватность выбранной модели, а расхождения могут использоваться для оценки влияния погодных условий при проектировании солнечных установок.

Суточная плотность энергии солнечного излучения [Вт·ч/м2 в день] в течение года на горизонтальную поверхность в провинции Нинь Тхуан показана на рис. 5.

График суточной плотности энергии излучения за год с сильными колебаниями фактических данных вокруг расчётной кривой
Рис. 5. Суммарная плотность энергии солнечного излучения на горизонтальную поверхность в провинции Нинь Тхуан в течение года

Для учёта влияния атмосферных факторов на поступление солнечного излучения была применена схема поправочного пересчёта, основанная на использовании коэффициента K, определяемого как отношение суммарной фактической радиации за месяц (по данным NASA POWER) к расчётному значению по модели Берда.

Расчёт коэффициента K проводился в два этапа. На первом этапе определялись суточные значения поправочного коэффициента Kгорiпр для каждого дня в месяце:

Kгорiпр = EгорiпрР / Eгорiпрч. (9)

На втором этапе, для получения обобщающих показателей по характерным месяцам, вычислялось средневзвешенное значение коэффициента за месяц:

Kмес = (1/n)∑i=1nKгорiпр dt, (10)

где Kгорiпр — поправочный коэффициент для i-го дня; n — количество дней в месяце.

После определения коэффициентов для каждого месяца, была проведена обратная процедура: теоретические почасовые значения умножались на соответствующий Kгорiпр для получения приближённых значений фактической радиации.

Табл. 1. Фактические и теоретические значения суммарного излучения по месяцам и определение поправочного коэффициента

Месяц Wтеор, кВт·ч/м2 Кмес Wфакт, кВт·ч/м2
1 189,734 0,70 132,440
2 189,222 0,84 158,210
3 231,073 0,89 204,750
4 229,603 0,77 176,930
5 247,88 0,73 181,950
6 238,920 0,68 162,750
7 247,532 0,61 150,590
8 239,520 0,78 187,520
9 229,795 0,67 153,160
10 217,238 0,70 152,290
11 188,502 0,66 123,530
12 183,140 0,79 144,450

Совокупное значение по всем часам месяца при этом точно совпало с исходной суммой по данным NASA (табл. 1), что и ожидаемо с математической точки зрения. Это подтверждает корректность подхода и позволяет использовать такой способ пересчёта в случаях отсутствия подробных почасовых измерений, для приближенной оценки выработки фотоэлектрических установок.

Из табл. 1 видно, что максимальное значение солнечной радиации Wфакт приходится на март, равно 204,75 кВт·ч/м2. Вьетнам расположен между 8° и 23° с. ш., то есть близко к экватору. В период весеннего равноденствия (около 21 марта) солнечные лучи падают почти перпендикулярно к поверхности Земли, особенно в южных регионах страны, таких как Нинь Тхуан (около 11,4° с. ш.). Это обеспечивает максимальный угол возвышения Солнца и, как следствие, наибольшую плотность солнечной радиации.

Расчёт плотности солнечного излучения, поступающего на ориентированные поверхности

Для расчёта плотности прямого солнечного излучения, поступающего на поверхность фотоэлектрической установки в условиях ясного неба, используется следующее выражение:

Eследпр = Ec⊥τriτO3iτгазiτH2OiKследпр. (11)

Среднемесячные дневные суммы прямого солнечного излучения на следящую поверхность при ясном небе вычисляются интегрированием по времени в течение светового дня:

Wследjпрч = ∑t=tмвосхtмзахEследjпрч(t) dt. (12)

Для перехода от расчётов при ясном небе к реальным атмосферным условиям вводится коэффициент коррекции, рассчитываемый для каждого месяца:

Kследjпр = Wследjпрн / Wследjпрч. (13)

Для неподвижной ориентированной поверхности фотоэлектрической установки плотность прямого солнечного излучения в реальных условиях рассчитывается по выражению:

EнаклiпрР = EследiпрчKследiпрcosθi, (14)

где cosθi — косинус угла падения лучей на ориентированную поверхность.

Для рассматриваемого случая ориентации фотоэлектрических модулей на юг (γ = 0) косинус угла падения солнечных лучей на наклонную поверхность определяется выражением:

cosθнаклi = sinδnsinφмcosβ − sinδncosφмsinβ + cosδncosφмcosβ cosωci + cosδnsinφмsinβ cosωci. (15)

Для оценки влияния угла наклона солнечной панели на величину годового прихода солнечной энергии были проведены расчёты суммарной радиации Wгод на наклонную поверхность при изменении угла β от 0° до 20° с шагом 1°.

Расчёты выполнены для трёх характерных регионов Вьетнама [северного (φ = 20° с. ш.), центрального (φ = 11,4° с. ш.) и южного (φ = 9° с. ш.)], что позволяет оценить изменение оптимального угла наклона вдоль широтного градиента страны.

Для количественной оценки влияния отклонения угла наклона от оптимального значения проведён расчёт относительных потерь солнечной энергии ΔW %, определяемых как процентное отклонение годового прихода солнечной радиации Wβ от максимального значения Wmax:

ΔW % = (Wmax − Wβi)/Wmax·100 %. (16)

Расчёт выполнен для диапазона углов β = 0–20° и для трёх регионов Вьетнама (северного, центрального и южного).

Результаты приведены в табл. 2. Как видно, оптимальный угол наклона солнечной панели для территории Вьетнама находится в диапазоне от 0° до 2°.

Табл. 2. Годовые значения суммарной солнечной радиации Wгод на наклонной поверхности и относительные потери ΔW при различных углах наклона β для регионов Вьетнама

β Wгод, кВт·ч/м2 (на севере) ΔW, % (на севере) Wгод, кВт·ч/м2 (на юге) ΔW, % (на юге) Wгод, кВт·ч/м2 (в центре) ΔW, % (в центре)
0 1143,983 0,048 1234,746 0 1230,035 0
1 1144,459 0,006 1233,861 0,072 1229,149 0,072
2 1144,534 0 1232,586 0,175 1227,879 0,175
3 1144,366 0,014 1230,962 0,306 1226,226 0,310
4 1143,796 0,064 1228,950 0,469 1224,191 0,475
5 1142,879 0,145 1226,564 0,663 1221,773 0,672
6 1141,613 0,255 1223,804 0,886 1218,974 0,899
7 1139,999 0,396 1220,671 1,140 1215,795 1,158
8 1138,038 0,568 1217,167 1,424 1212,237 1,447
9 1135,730 0,769 1213,292 1,738 1208,301 1,767
10 1133,076 1,001 1209,047 2,081 1203,989 2,117
11 1130,078 1,263 1204,434 2,455 1199,303 2,499
12 1126,734 1,555 1199,454 2,858 1194,243 2,910
13 1123,048 1,877 1194,108 3,291 1188,811 3,351
14 1119,020 2,229 1188,399 3,754 1183,011 3,823
15 1114,651 2,610 1182,328 4,245 1176,843 4,324
16 1109,942 3,022 1175,897 4,766 1170,310 4,856
17 1104,895 3,463 1169,108 5,316 1163,414 5,416
18 1099,512 3,934 1161,962 5,895 1156,159 6,006
19 1093,793 4,433 1154,462 6,502 1148,545 6,625
20 1087,742 4,962 1146,611 7,139 1140,577 7,273

Для северного региона (φ = 20° с. ш.) максимальный приход солнечной радиации достигается при угле наклона около 2°, тогда как для центрального (φ = 11,4° с. ш.) и южного (φ = 9° с. ш.) регионов оптимальным является горизонтальное расположение поверхности (β ≈ 0°).

Отметим, что для большинства крышных фотоэлектрических установок, эксплуатируемых во Вьетнаме, угол наклона панелей обычно составляет от 10° до 15°. Это связано с конструктивными особенностями кровельных систем и необходимостью установки модулей с небольшим наклоном, обеспечивающим естественный сток дождевой воды и облегчённое обслуживание, включая очистку и техническое обслуживание панелей.

Таким образом, полученные расчётные значения оптимальных углов (0–2°) и определённые потери прихода при изменении угла наклона хорошо согласуются с реальными условиями монтажа солнечных батарей на крышах зданий во Вьетнаме.

Расчёт энергетического потенциала

Для оценки солнечных ресурсов территория Вьетнама была разделена на три характерные зоны с использованием геоинформационной системы (ГИС-технологии, QGIS). Районирование проведено на основе данных о среднегодовой плотности солнечного излучения, полученных из ресурсов Global Solar Atlas.

Данный подход позволил визуализировать пространственное распределение энергии и выделить регионы с наибольшим потенциалом для размещения фотоэлектрических систем.

На первом этапе с использованием программного комплекса QGIS была разработана карта районирования территории страны (рис. 6). В соответствии с распределением плотности солнечного излучения, территория Вьетнама была разделена на три основные зоны: северную, центральную и южную.

Карта Вьетнама с тремя цветовыми зонами — жёлтой, оранжевой и красной — с легендой площади и диапазонов радиации зон
Рис. 6. Распределение солнечного энергетического ресурса по трём зонам Вьетнама

Для количественной оценки ресурсов в рамках выделенных зон были рассчитаны два вида потенциала: валовый и технический.

Валовый потенциал WΣвэс характеризует суммарную солнечную энергию, поступающую на всю горизонтальную поверхность региона при идеальных условиях. Результаты расчёта валового потенциала для каждой зоны представлены в табл. 3. При этом в качестве расчётного значения Wi использовался удельный годовой приход суммарной солнечной радиации на оптимально наклонённую поверхность, что позволяет оценить максимально достижимый энергетический ресурс в каждой из выделенных зон.

Табл. 3. Валовый потенциал солнечной энергии каждой зоны территории Вьетнама

Зоны Площадь зоны Si, км2 Wгодi, ГВт·ч/км2 WiSi, ГВт·ч в год
1 (северный) 163 013 1114,534 181 683 530,94
2 (центральный) 103 402 1230,035 127 188 079,07
3 (южный) 62 698 1234,746 77 416 104,71
Итого     386 287 714,72

Технический потенциал Wтi определяет долю энергии, которая может быть фактически преобразована в электрическую с учётом ограничений по площади использования земель и эффективности современных фотоэлектрических систем.

В табл. 4 приведены значения технического потенциала в зависимости от коэффициента использования территории h (от 0,5 до 1,5 %).

Общий валовый потенциал солнечной энергии на оптимальную поверхность всей территории Вьетнама определяется как:

WΣвэс = ∑WiSi, (17)

то есть WΣвэс = 386 287,715 ТВт·ч/год.

Для перехода от валового потенциала к техническому потенциалу Wiт необходимо учитывать ограничения по площади, доступной для строительства фотоэлектрических систем (принятой как доля h от общей площади территории), а также общий коэффициент полезного действия (КПД) фотоэлектрических систем ηФЭС. В данном исследовании принят коэффициент использования площади и суммарный КПД системы ηФЭС = 18 %:

Wтi = WiSihηФЭС. (18)

Учитывая неопределённость в распределении земельного фонда, в данной работе проведён анализ технического потенциала для диапазона коэффициента использования площади h от 0,5 до 1,5 %. Данный интервал соответствует наиболее вероятным сценариям развития солнечной энергетики без ущерба для сельскохозяйственных и лесных угодий Вьетнама.

Следовательно, общий технический потенциал солнечной энергии на оптимальную поверхность территории Вьетнама определяется по формуле:

WΣт = ∑Wтi. (19)

Результаты расчётов при различных значениях h представлены в табл. 4. Это позволяет оценить вариативность энергетического выхода в зависимости от интенсивности использования доступных территорий.

Табл. 4. Потенциал солнечной энергии, приходящейся на территорию Вьетнама

Площадь h / потенциал солнечной энергии Wтi 0,005 0,007 0,009 0,01 0,011 0,012 0,015
Зона 1 (север) 163 515 228 921 294 327 327 030 359 733 392 436 490 546
Зона 2 (центр) 114 469 160 257 206 045 228 939 251 832 274 726 343 408
Зона 3 (юг) 69 674 97 544 125 414 139 349 153 284 167 219 209 023
WΣт 347 659 486 723 625 786 695 318 764 850 834 381 1 042 977

Как показывают расчёты, даже при минимально принятом значении h = 0,5 % (0,005), суммарный технический потенциал исследуемых регионов составляет примерно 348 ТВт·ч/год. Данный показатель уже на 60 % превышает общий объём годовой выработки электроэнергии во Вьетнаме, который в 2023 году достиг 217 ТВт·ч [11]. При увеличении коэффициента h до базового уровня 1 % (0,01), потенциал возрастает до 695 ТВт·ч/год, что более чем в три раза покрывает энергетические потребности страны. Это открывает широкие возможности не только для полной декарбонизации энергетического сектора, но и для экспорта электроэнергии в соседние страны региона.

Заключение

  1. Разработана и апробирована методика расчёта прихода солнечной радиации на горизонтальные и наклонные поверхности территории Вьетнама с использованием модели Берда для условий ясного неба и изотропной модели Лю и Джордана для наклонных площадок.
  2. Проведённые расчёты показали, что оптимальный угол наклона солнечных панелей на территории Вьетнама изменяется в пределах 0–2°: около 2° на севере и близкий к 0° в центральной и южной частях страны. При этом относительные потери годового прихода солнечной энергии при изменении угла наклона β в диапазоне 0–20° не превышают 8 %.
  3. Показано, что даже при фиксированном положении панелей без сезонного регулирования угла наклона энергетические потери незначительны, что делает такие конструкции рациональными для климатических условий Вьетнама.
  4. Разработанная методика и полученные результаты могут быть использованы для оценки солнечного потенциала и выбора оптимальных параметров установки фотоэлектрических систем, а также при создании атласов солнечной радиации для других регионов Юго-Восточной Азии.
  5. На основе ГИС-технологий и данных Global Solar Atlas проведено районирование территории Вьетнама, позволившее оценить валовый потенциал в 386 287 ТВт·ч/год. Установлено, что при использовании лишь 1 % площади территории технический потенциал составляет 695 ТВт·ч/год, что в 2,5–3 раза превышает современное электропотребление страны. Это подтверждает стратегическую возможность полного обеспечения энергетических нужд Вьетнама за счёт солнечной генерации.