Введение. Формирование условий для социально-экономической стабильности в обществе предполагает необходимость формирования достаточных объемов и рациональной структуры продовольственных ресурсов. Важная роль в решении этой задачи принадлежит круглогодичному обеспечению населения высококачественной и разнообразной овощной продукцией [1–3].В настоящее время в сельском хозяйстве активно ведется строительство новых и эксплуатация существующих теплиц в модульном (индивидуальном) и блочном исполнениях [4]. Решением проблемы выращивания рассады и овощей в защищенном грунте является проектирование не только оптимальной конструкции зимней теплицы, но и системы отопления, которая бы обеспечивала комфортные условия для растений при резких колебаниях температуры наружного воздуха в ночное и дневное время холодного периода года [5].Современные тенденции тепличного строительства направлены на децентрализацию теплоснабжения зимних теплиц, т. е. на обогрев помещений без каких-либо промежуточных теплоносителей. Данное техническое решение позволяет увеличить коэффициент использования топлива примерно в 1,5÷2,0 раза в сравнении с традиционными системами отопления [6]. При этом использование солнечной энергии (гелиотеплицы) [7, 8], тепловых насосов [9, 10] и геотермальных источников энергии [11] для обогрева зимних теплиц по ряду причин, к которым можно отнести, в первую очередь, зависимость производства теплоты от географических и климатических условий местности, является во многих случаях нецелесообразным вариантом.Высокая стоимость энергоресурсов, которая порой приводит к нерентабельности всего тепличного производства, заставляет руководителей предприятий переходить на энергоэффективные технологии обогрева, например, с использованием в качестве источников теплоты длинноволновых инфракрасных излучателей, работающих в основном на газообразном топливе [12]. Это позволяет одновременно создать условия «солнечного» отопления в теплице при минимальных энергозатратрах и использовать продукты сгорания газообразного топлива (CO2) для питания растений. Однако, как показали проведенные исследования, применение исключительно инфракрасного обогрева почвы не позволяет, особенно в зимний период года, поддерживать в теплице оптимальную (или допустимую) температуру внутреннего воздуха, согласно данным табл. 1 [13].Таблица 1Температурно-влажностный режим теплицы (до плодоношения)КультураТемпература воздуха, ºСТемпература грунта, ºСОтносительная влажность воздуха, %деньночьсолнечнопасмурноОгурец (зимне-весенний оборот)22÷2420÷2217÷1820÷2470÷75Огурец (осенний оборот)25÷2622÷2319÷2022÷2470÷75Томат (зимне-весенний оборот)22÷2419÷2016÷1718÷2060÷65Томат (осенний оборот)24÷2618÷2016÷1818÷1960÷70Салат кочанный20÷2316÷181015÷1670÷80Редис20÷227÷95÷615÷1660÷70Укроп, шпинат17÷188÷125÷615÷1665÷80  Поэтому имеет смысл рассмотреть вариант комбинированного отопления зимней теплицы, когда тепловой режим почвы формируется за счет потолочных инфракрасных излучателей, а подогрев внутреннего воздуха до требуемых значений (табл. 1) осуществляется с помощью автономных конвективных обогревателей (например, электрокалориферных установок). Здесь стоит оговориться, что режим их работы неравномерный в течение отопительного сезона (зависит от температуры наружного воздуха) и возможно полное отключение конвективных обогревателей в осенний и весенний периоды года, когда система лучистого отопления справляется самостоятельно.Если коснуться существующих методик расчета лучистого отопления зимних теплиц, то они не лишены недостатков. Так, например, в работе [12] совершенно не рассматриваются процессы массопереноса, происходящие в теплице с инфракрасными излучателями, которые напрямую влияют на тепловлажностный режим культивационного сооружения. Кроме того, не учтена специфика самого инфракрасного излучения, связанная, например, с многократностью отражений лучистого потока теплоты от поверхностей теплообмена. Также не стоит забывать, что поддержание требуемых температурных условий почвы при инфракрасном обогреве вовсе не дает гарантий соблюдения климатических условий в самом помещении в соответствии с табл. 1, так как температура внутреннего воздуха является зависимой величиной. То же самое относится и к методике расчета, представленной в научной статье [3], где акцентировано внимание не на самом отоплении зимней теплицы, а на теплообмене излучением между источником теплоты и окружающими объектами (в первую очередь, полом). В публикации [14] уделено внимание только температуре воздуха в помещении при работе инфракрасных излучателей. Таким образом, остается актуальной проблема разработки метода расчета децентрализованной лучисто-конвективной системы отопления, обеспечивающей требуемые параметры микроклимата (табл. 1) в течение всего календарного года.Алгоритм расчета лучисто-конвективного (комбинированного) отопления зимней теплицы. Одним из вариантов решения обозначенной проблемы, в первую очередь, для эксплуатируемых зимних теплиц, является установка потолочных инфракрасных излучателей при существующей водяной системе отопления (для новых теплиц, очевидно, имеет смысл применение автономных конвективных обогревателей), что позволит уменьшить тепловую нагрузку традиционной системы теплообеспечения. Сущность метода расчета лучисто-конвективного (комбинированного) отопления заключается в определении всех тепловых и массовых потоков, поступающих и уходящих из теплицы в окружающую среду, и нахождении неизвестных параметров микроклимата, формирующих тепловлажностный режим культивационного сооружения (табл. 1).На рис. 1 представлена расчетная схема лучисто-конвективного отопления зимней теплицы.Тепловлажностный режим зимней теплицы (рис. 1), согласно табл. 1, определяется температурой tвн.в, °С , и относительной влажностью φвн.в , %  , внутреннего воздуха, а также температурой поверхности почвы 1, равной tпов, °С . Влагосодержание поверхности почвы 1 имеет постоянное значение Wпов , кг/кг. Температура и относительная влажность наружного воздуха также известны и соответственно составляют tвн.в, °С , и φвн.в , % . Температура внутренней поверхности ограждения 2 теплицы tогр , ºС, неизвестна и находится через решение системы уравнений теплового баланса относительно рассматриваемой поверхности. Рис. 1. Расчетная схема лучисто-конвективного отопления теплицы: 1 – почва; 2 – ограждение; 3 – источник инфракрасного излучения (излучатель); 4 – автономные конвективные обогреватели; 5 – вентиляционный проем (приточный, вытяжной)Количество теплоты, поступающей в единицу времени на внутреннюю поверхность ограждения 2 зимней теплицы, равно:Qогр = (1 – (A1/(1 – kотр)))Q изл+ Qлуч + Qконв2, Вт , (1)где kотр = (1 – А1)(1 – А2)φ21/(1 – φ22(1 – А2))  – коэффициент, учитывающий многократное отражение теплового излучения от поверхности почвы 1 и внутренней поверхности ограждения 2 теплицы; А1 и А2 – соответственно коэффициенты поглощения поверхности почвы 1 и внутренней поверхности ограждения 2 теплицы; φ21 – коэффициент облученности с внутренней поверхности ограждения 2 теплицы на поверхность почвы 1; φ22 – коэффициент самооблученности внутренней поверхности ограждения 2 теплицы; Qизл – лучистая составляющая тепловой мощности комбинированной системы отопления, Вт; Qлуч – результирующее тепловое излучение, возникающее между поверхностью почвы 1 и внутренней поверхности ограждения 2 теплицы, Вт; Qконв2 – конвективная составляющая теплообмена между внутренним воздухом и внутренней поверхностью ограждения 2 теплицы, Вт.Лучистая составляющая тепловой мощности комбинированной системы отопления, согласно действующему своду правил [15], не должна превышать максимальной величины (полагаем, что все инфракрасное излучение, испускаемое излучателем 3, падает на поверхность почвы 1 теплицы), равной:Qизл = qизлFпов , Вт,                    (2)где qизл  – максимальная интенсивность инфракрасного излучения, Вт/м2. В соответствии со сводом правил [15] величина qизл  не должна превышать 150 Вт/м2 на постоянных и 250 Вт/м2 на непостоянных рабочих местах; Fпов – площадь поверхности почвы 2 теплицы, м2.Тепловой поток излучением между поверхностью почвы 1 и внутренней поверхностью ограждения 2 теплицы при условии, что tпов > tогр  и φ12=1  (все тепловое излучение, идущее с поверхности почвы 1, полностью попадает на внутреннюю поверхность ограждения 2 теплицы), определяется по закону Стефана-Больцмана:Qлуч = с0εпр12Fпов((Tпов/100)4 – (Tогр/100)4) , Вт, (3)где с0 – коэффициент излучения абсолютно черного тела, 5,67 Вт/(м2·К4); εпр12  – приведенный относительный коэффициент теплового излучения поверхности почвы 1 и внутренней поверхности ограждения 2 теплицы; Tпов= tпов+273,15  и Tогр= tогр+ 273,15  – соответственно абсолютная температура поверхности почвы 1 и внутренней поверхности ограждения 2, К.Конвективная составляющая теплообмена между внутренним воздухом и внутренней поверхностью ограждения 2 теплицы рассчитывается по закону конвективного теплообмена Ньютона-Рихмана:Qконв2 = αвн(tвн.в – tогр)Fогр,  Вт,        (4)где αвн  – коэффициент теплоотдачи внутренней поверхности ограждения 2 теплицы, Вт/(м2·К); Fогр  – суммарная площадь ограждения 2 теплицы, м2.В то же время тепловой поток, уходящий через ограждение 2 теплицы, по уравнению теплопередачи будет равен:Qогр = ((tогр – tн.в)/Rt)Fогр,  Вт,            (5)где Rt = Rогр + 1/αн  – сопротивление теплопередаче ограждения 2 теплицы, м2·К/Вт; Rогр  – термическое сопротивление ограждения 2 теплицы, м2·К/Вт; αн  – коэффициент теплоотдачи наружной поверхности ограждения 2, Вт/(м2·К). По данным таблицы 6 свода правил [16] αн= 23 Вт/(м2·К).Конвективная составляющая тепловой мощности комбинированной системы отопления находится через уравнение теплового баланса зимней теплицы:Q изл+Q нагр=Qогр+Qвент+Qгр , Вт,   (6)где Q нагр – конвективный тепловой поток, идущий от автономных конвективных обогревателей 4, Вт; Qвент  – потери теплоты с вентиляционным воздухом, уходящим из теплицы через вытяжной вентиляционный проем 5, Вт; Qгр  – потери теплоты в глубинные горизонты грунта, Вт.Потери теплоты с уходящим вентиляционным воздухом численно равны расходу теплоты на нагрев воздуха, поступающего из окружающей среды в теплицу через приточный вентиляционный проем 5:Qвент = Gв(h вн.в– hн.в) , Вт,             (7)где  Gв – массовый расход сухой части влажного воздуха, участвующего в воздухообмене в теплице, кг/c; h вн.в и hн.в – соответственно удельная энтальпия внутреннего и наружного влажного воздуха, Дж/кг.Потери теплоты в грунт упрощенно вычисляются по уравнению:Qгр=(tвн.в –tн.в)Σi=1n(Fi/Ri) ,    Вт,         (8)где Fi  – расчетная площадь i-ой зоны грунта в теплице при их общем количестве n, м2; Ri  – сопротивление теплопередаче i-ой зоны грунта в теплице при их общем количестве n, м2·К/Вт. Для первой зоны сопротивление теплопередаче равно RI=2,1 м2·К/Вт; для второй – RII= 4,3  м2·К/Вт; для третьей – RIII=8,6 м2·К/Вт; для четвертой – RIV=14,2  м2·К/Вт [17].Уравнение теплового баланса поверхности почвы 1 связывает между собой тепловые и массообменные процессы, происходящие в теплице:А1Qизл1 – kотр=Qлуч+Qконв1+Qисп+Qгр , Вт,       (9)где Qконв1  – тепловой поток, вызванный конвективным теплообменом между поверхностью почвы 1 и приземным слоем воздуха в теплице (на рис. 1 обозначен как Qконв ), Вт; Qисп  – тепловой поток, расходуемый на испарение влаги с поверхности почвы 1, Вт.Конвективный теплообмен между поверхностью почвы 1 и приземным слоем воздуха в теплице также подчиняется закону Ньютона-Рихмана:Qконв1 = αк(tпов – tвн.в)Fпов , Вт,       (10)где αк  – коэффициент теплоотдачи поверхности почвы 1, Вт/(м2·К).Тепловой поток, связанный с процессом парообразования на поверхности почвы 1, рассчитывается по уравнению:Qисп = Gиспr , Вт,                    (11)где Gисп  – испарение влаги с поверхности почвы 1, кг/с; r – удельная теплота фазового превращения (перехода), Дж/кг.При стационарном тепловлажностном режиме в зимней теплице интенсивность испарения влаги с поверхности почвы 1 должна быть равна количеству воды, уходящему из теплицы вместе с вентиляционным воздухом, т. е.:G исп=Gвент=Gв(dвн.в – dн.в) , кг/с,      (12)где Gвент т – потери влаги с уходящим из теплицы воздухом, кг/с; dвн.в  и  dн.в  – соответственно влагосодержание внутреннего и наружного влажного воздуха, кг/кг.Если на данном этапе расчета не учитывать количество влаги, поглощаемое растениями в теплице, то расход воды на полив почвы будет равен:Gпол = Gисп , кг/с,                   (13)где Gпол  – расход воды на полив почвы, кг/ч.Стоит отметить, что при разработке метода расчета комбинированной системы отопления зимней теплицы, принято во внимание ряд допущений: во-первых, расчет теплообмена через ограждение культивационного сооружения по уравнению теплопередачи (5) выполняется как для плоской стенки, ввиду относительно малой толщины конструкции, обычно не превышающей 32 мм; во-вторых, потолочные инфракрасные излучатели рассматриваются как линейные источники теплоты, и поэтому затенения ими ограждающих конструкций теплицы не учитываются.Пример реализации способа. Рассмотрим на примере индивидуальной промышленной теплицы «Фермер 7.5» (рис. 2) способ расчета комбинированной системы отопления.Рис. 2. Промышленная теплица «Фермер 7.5» и схема разбивки грунта на зоныИсходные данные для выполнения расчета:1. Размеры (габариты) теплицы: ширина a = 7,45 м; длина b = 8,40 м; высота h = 3,80 м.2. Параметры поверхности почвы (томаты до плодоношения в зимне-весенний оборот, согласно данным табл. 1): температура tпов=18 ºС ; коэффициент поглощения А1 = 0,65 (коэффициент отражения R1 = 0,35); коэффициент теплового излучения (степень черноты) ε1 = А1 = 0,65.3. Параметры ограждения: материал – сотовый поликарбонат с термическим сопротивлением ограждения Rогр=0,25  м2·К/Вт; коэффициент поглощения А2 = 0,94 (коэффициент отражения R2 = 0,06); коэффициент теплового излучения (степень черноты) ε2 = А2 = 0,94.4. Расчетные параметры внутреннего воздуха (табл. 1): температура tвн.в=22℃ ; относительная влажность φвн.в=65% .5. Расчетные параметры наружного воздуха для условий г. Вологды [18]: температура tн.в=–32°С ; относительная влажность φн.в=85%. 6. Параметры воздухообмена: принудительная вентиляция; кратность воздухообмена для индивидуальной теплицы – nв = 1,0 м3/мин на 1 м2 площади помещения [19].7. Параметры орошения: коэффициент орошения почвы kор = 1 (орошается вся поверхность почвы в теплице).8. Условие пребывания людей в рабочей зоне: непостоянное пребывание рабочего персонала на своих местах.Расчет лучисто-конвективного отопления промышленной теплицы «Фермер 7.5» (рис. 2) выполнен в математическом редакторе “Mathcad”.По итогам программного расчета лучисто-конвективного отопления промышленной теплицы «Фермер 7.5» получены следующие результаты:1. Требуемая мощность инфракрасного излучения Qизл≈15,6  кВт и конвективной системы отопления Qнагр≈133,7  кВт.2. Необходимый расход воды на полив почвы Gпол≈58,5  кг/ч.3. Потери теплоты: через ограждение теплицы Qогр≈23,4  кВт; с вентиляционным воздухом, уходящим из теплицы в окружающую среду, Qвент≈124,2  кВт; в грунт Qгр≈1,82  кВт.4. Температура внутренней поверхности ограждения теплицы tогр≈8,8℃ .5. Потери теплоты от теплообмена излучением между поверхностью почвы и внутренней поверхностью ограждения теплицы Qлуч≈1,96  кВт.6. Потери теплоты от конвективного теплообмена между внутренним воздухом и поверхностью почвы в теплице Qконв1≈0,08  кВт (в данном случае tвн.в>tпов °С ); между внутренним воздухом и внутренней поверхностью ограждения теплицы Qконв1≈16,1  кВт.7. Потери теплоты, связанные с процессом испарения влаги с поверхности почвы в теплице, составили Qисп≈40,0  кВт.На рис. 3 показана зависимость расчетной тепловой мощности системы лучисто-конвективного отопления Qсум = Qизл + Qнагр , Вт, и температуры внутренней поверхности ограждения ,℃ , зимней теплицы (рис. 2) от температуры наружного воздуха tн.в,℃ , в течение отопительного сезона. Температура внутренней поверхности ограждения tогр  определяется исходя из решения системы взаимосвязанных уравнений тепломассопереноса (1), (5), (6), (9), (12) и (13).Рис. 3. Тепловая мощность системы отопления и температура ограждения теплицыСогласно рис. 3, с повышением температуры наружного воздуха tн.в,℃ , наблюдается уменьшение суммарной тепловой мощности комбинированной системы отопления Qсум , Вт, что является вполне закономерным. Отопление служит для компенсации тепловых потерь в помещении, которые снижаются при росте температуры окружающей среды [20–23]. Ввиду того, что ограждение теплицы не обладает высокими теплозащитными свойствами (данная проблема является одной из ключевых в тепличном строительстве) температура tогр , ºС, в определенных случаях может опускаться ниже 10 ºС. Поэтому стоит рассмотреть вариант предварительного нагрева наружного воздуха для нужд вентиляции теплицы в случае неблагоприятных климатических условий.Результаты исследования. Конвективная составляющая тепловой мощности комбинированной системы отопления теплицы оказалась гораздо больше лучистой. Это объясняется не только низкой температурой наружного воздуха tн.в,℃ , но и ограничениями по интенсивности инфракрасного излучения. С ростом температуры наружного воздуха доля лучистого отопления в тепловом балансе комбинированного отопления, очевидно, будет неуклонно увеличиваться. В любом случае рассмотренный способ обогрева позволяет поддерживать как требуемую температуру почвы, так и температуру внутреннего воздуха в зимней теплице без подключения к системе централизованного теплоснабжения (СЦТ). Как показали расчеты, расход теплоты на нужды воздухообмена оказался достаточно большим, поэтому имеет смысл рассмотреть также вариант предварительного подогрева приточного воздуха до расчетной температуры перед его подачей в теплицу.Заключение и выводы. Предложенный способ расчета может быть использован для проектирования комбинированных систем отопления культивационных сооружений, включающих как конвективные отопительные приборы, так и инфракрасно-лучистые обогреватели.Перспективой дальнейших исследований является разработка инженерной методики расчета системы лучисто-конвективного отопления теплицы культивационного сооружения, которая позволит проектировщику в максимально сжатые сроки и по предельно упрощенному алгоритму, без решения сложной системы взаимосвязанных уравнений теплового и материального балансов, определять искомые характеристики при заданных расчетных условиях.