CC BY
80
УДК 621.311.25:551.521.1
В. В. Чертищев, Г.Н. Филиппова Использование справочных метеорологических данных для расчета теплопроизводительности солнечных коллекторов
На протяжении многих столетий солнечная энергия использовалась для различньх целей. Горячее водоснабжение для бытовых нужд с помощью солнечных водонагревательных установок стало нормой жизни в ряде стран, в том числе в Австралии, Израиле, Японии. Отопление помещений для создания необходимого комфорта также может осуществляться с помощью системі, использующих солнечную энергию [1] .
В настоящее время из активных систем утилизации солнечной энергии для отопления и горячего водоснабжения наиболее доступными являются системы на основе солнечных коллекторов. Солнечный коллектор представляет собой теплообменное устройство, в котором происходит поглощение солнечной энергии с преобразованием ее в теплоту циркулирующего теплоносителя. Конструктивно солнечные коллекторы отличаются значительным разнообразием как формы (плоские, параболические), так и технологических параметров (с прозрачным одинарным и двойным покрытием!, селективные, вакуумированные) . Они работают с естественной и принудительной циркуляцией теплоносителя. Основной показатель эффективности солнечного коллектора - это его КПД [ 2 ] . Наибольшее влияние на КПД, как известно, оказывают:
1) метеорологические факторы - интенсивность солнечной радиации и температура наружного воздуха;
2) конструктивные характеристики - поглощательная способность абсорбера, пропускательная способность прозрачного покрытия, теплоизоляционные свойства материалов и конструкции;
3) монтажно-технологические параметры -ориентация, температура теплоносителя на входе в коллектор, расход теплоносители.
Прежде чем принять решение об использовании солнечной энергии для этих целей, необходимо оценить имеющееся количество солнечной энергии, зависящее, в частности, от угла наклона облучаемой солнцем поверхности.
В работе предлагается методика оценки располагаемого количества солнечной энергии, основанная одновременно:
а) на закономерностях взаимного движения Земли и Солнца, позволяющих рассчитывать ин-
тенсивность прямого солнечного излучения для каждого момента времени вне земной атмосферы;
б) на использовании региональных справочных метеоданных для оценки прозрачности атмосферы и потока, рассеянного атмосферой и поверхностью Земли излучения.
Данная методика лежит в основе разработанного нами программного модуля для расчета тепловой мощности солнечного коллектора. В программном модуле для каждого момента времени вычисляется результирующая (складывающаяся из потока прямой солнечной радиации, солнечной радиации, рассеянной в атмосфере, и солнечной радиации, отраженной от поверхности) энергетическая освещенность наклонной поверх-ности, и по заданным (или вычисленным по соответствующему алгоритму) температурам наружного воздуха и входящего в коллектор теплоносителя вычисляется КПД; затем по соотношению (2) рассчитывается тепловая мощность коллектора.
КПД коллектора находится из балансового соотношения для энергии [2] :
Л - Е. + рк • и • (Твх - Тнаруж ) = та - Е., (1)
где первое слагаемое в левой части - полезная часть энергии, которая передается потребителям; второе слагаемое - часть энергии, теряющаяся теплоотдачей в окружающую атмосферу, правая часть соотношения - та часть падающей на коллектор солнечной энергии, которая воспринимается коллектором. Здесь Твх и Тшруж - температуры входящего в коллектор теплоносителя и наружного воздуха; - и - произведение общего коэффициента теплопотерь коллектора и коэффициента эффективности поглощающей панели (то и другое на 1 м2 площади поверхности коллектора) ; - Ха - произведение оптического КПД и
коэффициента эффективности поглощающей панели коллектора. По КПД рассчитывается тепловая мощность солнечного коллектора
дк =л- е. . (2)
Суммарный поток солнечной радиации на данную поверхность складывается из потоков прямого излучении, излучения, рассеянного атмосферой, и излучения, отраженного от поверхности
Использование справочные метеорологическ^к данные.
Земли и различных близко расположенных предметов. Обычно (за исключением зимних месяцев и отдельных районов) наибольшую долю в суммарной радиации составляет прямая радиация [2-3] .
Результирующая энергетическая освещенность Es наклонной поверхности при безоблачном небе, если принять допущении, что небо может рассматриваться как изотропный источник радиации и что рассеяние света земной поверхностью близко к изотропному, может быть представлена в виде
Es = E + Edh • (l + cos b)/ 2 + + (Edh + !m >r -(l - cos b)/2, (3)
где первое слагаемое - энергетическая освещенность наклонной поверхности прямой солнечной радиацией; второе слагаемое - энергетическая освещенность данной поверхности рассеянной радиацией; третье слагаемое - энергетическая освещенность данной поверхности отраженной радиацией [3-5] . Первое слагаемое определяется
выражением E = Im • COS i .
Входящий в эту формулу i-угол между солнечными лучами и нормалью к облучаемой поверхности рассчитывается по соотношению [3, 6]
cosi = sin в • cos а •
• cos(a - an) + cos в • sin а, (4)
где произведение COS a • COS(a — an) определяется по выражению
cos acos(a — an) = cos an (sin j (cos W cos w+
+ cos 8 sin w) — cos j sin 8) — sinan (5)
(cos 8 cos w — cos W sin w) + sin j sin 8,
здесь; j - широта местоположения точки; j -дневной часовой угол, а угол склонения Солнца и угловая высота Солнца над горизонтом определяются по формулам
sin 8 = sin W • sin e, cos 8 = sin W • cos e (6)
(e - угол наклона земной оси) ;
sin a = cos j • (cos 8 • sin w + cos W • cos w). (7)
I - интенсивность солнечной радиации у земной поверхности;
Im = Io • es • Pm, m = 1/sin a (8)
оптическая масса атмосферы;
I0 = 1,353 кВт/м 2 - солнечная постоянная;
es = 1 — 0,034 • sin W (9)
- коэффициент эксцентричности орбиты, W - часовой угол.
Поток солнечного излучения, достигающий по-вepxнocти Земли, зависит не только от времени суток, времени года и месторасположения, но также и от погоды, и в облачную погоду не может быть рассчитан по формуле (3) . Такие параметры, как коэффициент прозрачности атмосферы, коэффициент отражения излучения земной повер-xнocтью, диффузная составляющая инсоляции, рассчитываются по данным соответствующе мeтeopoлoгичecкиx измерений. Они различны для pазныx регионов.
При разработке методики оценки пржода солнечной радиации на наклонную площадку было отдано предпочтение наиболее полному использованию основный справочный мeтeoдaнныx, приведении в «Справочнике по климату СССР» [ 7 ] для каждого месяца, для времени суток с тpexчaco-вым интервалом, а также данный по облачности и по альбедо земной пoвepxнocти. С учетом облачности расчетная формула приняла вид
E, = E• P„ + Edh '(l + cosb)/2 +
+ (Edh + ]m •P.. •(l + cos(a - an))/ 2^
• r •(l — cos b)/2, (10)
где a — an - разность азимутов солнца и нормали к облучаемой пoвepxнocти, b - угол между наклонной пoвepxнocтью и горизонтальной плоскостью), r - коэффициент отражения радиации земной пoвepxнocтью.
Коэффициент прозрачности атмосферы рассчитывается по метеоданным об интенсивности прямого солнечного излучения с использованием соотношения
P = exp(sin a • Ь (Im /(I0 • es ))). (11)
Энергетическая освещенность горизонтальной пoвepxнocти рассеянной в атмосфере солнечной радиацией может быть рассчитана по формуле
Edh = D • P.. + Im • Di • (l — P„). (12)
Здесь P.. - вероятность того, что солнце не
закрыто облаками, и диффузная составляющая энергетической освещенности при безоблачном небе берутся непосредственно из метеорологичес-:hmx данный за длительный срок наблюдений. Диффузная составляющая энергетической освещенности при облачной погоде представляется в виде произведения интенсивности света у по-
вepxнocти Земли Im и некоторой функции двуx переменный - времени суток и времени года, эта
функция рассчитывается по соотношению
О; = ^ - (Е + Б) - Рис)/(1т - (1 - Рис)), (13)
где р - суммарная солнечная радиация при средних условиях облачности, усредненная по метеорологическим данным за десятилетний период для каждого месяца, для времени суток с трехчасовым интервалом.
Приведенные формулы лежат в основе разработанного программного модуля расчета энергетической освещенности произвольно ориентированной площадки, в нашем случае - воспринимающей плоской поверхности солнечного коллектора. При
вычислении Р,Рис,ЕёЬ для всех промежуточных
моментов времени используется двумерная интерполяция.
Расчет по описанной методике позволяет найти теплопроизводительность солнечного коллектора для заданного угла наклона к горизонту и оценить располагаемое количество энергии, а для различных углов наклона солнечного коллектора позволяет выбрать оптимальный для предполагаемого сезона использования угол наклона коллектора в любой местности, для которой известны перечисленные выше метеоданные.
Литература
1. Даффи Дж. А., Бекман У.А. Тепловые процессы с использованием солнечной энергии: Пер. с англ. / Под ред. Ю.Н. Малевского. М., 1977.
2. Валов М.И., Казанджан Б.И. Использование солнечной энергии в системах теплоснабжения: Монография. М., 1991.
3 . Харченко Н. В. Индивидуальные солнечные установки!. М., 1991.
4. Кондратьев К.Я., Пивоварова З.И., Федорова М.П. Радиационный режим наклонных поверхностей: Монография/ Под ред. К.Я. Кондратьева. Л., 1978.
5 . Крейт Ф., Блэк У. Основы теплопередачи: Пер. с англ. М., 1983.
6. Бринкворт Б. Дж. Солнечная энергия для человека. М., 1976.
7. Справочник по климату СССР. Вып. 20: Томская, Новосибирская, Кемеровская области, Алтайский край и Горно-Алтайская автономная область. Метеорологические данные за отдельные годы. Ч. 1: Солнечная радиация, радиационный баланс и солнечное сияние . Новосибирск, 1977.
