ЭНЕРГОСБЕРЕГАЮЩИЕ ТЕХНОЛОГИИ ПОРТОВЫХ СООРУЖЕНИЙ НА ОСНОВЕ ПРИМЕНЕНИЯ ГЕОТЕРМАЛЬНЫХ ТЕПЛОВЫХ НАСОСОВ
Аннотация и ключевые слова
Аннотация (русский):
Энергосберегающие технологии портовых сооружений и судостроительных предприятий на основе геотермальных тепловых насосов являются весьма перспективным направлением. Особенностями прибрежных объектов является их нахождение в зоне влияния фильтрационных потоков грунтовых вод, значительно осложняющих внедрение энергосберегающих геотермальных систем. Представлены результаты математического моделирования температурного поля грунтового массива, формируемого при эксплуатации вертикального грунтового теплообменника геотермального теплового насоса на протяжении 5-ти лет. Разработанная методика базируется на уравнении нестационарной теплопроводности и установке граничных условий для различных режимов работы теплового насоса. Приведены результаты изменения температурного поля по трем режимам работы геотермальной установки: без изменения направления теплового потока, с изменением направления теплового потока (регенерация) и в условиях, осложненных влиянием фильтрационного потока грунтовых вод. Исследования проведены с целью построения расчётных моделей, адекватных реальным физическим процессам. Для обобщения полученных результатов использовался метод теории подобия.

Ключевые слова:
энергосберегающие технологии, прибрежные объекты, тепловой насос, геотермальная скважина, температурное поле грунта, математическая модель
Текст
Введение Применение энергосберегающих технологий на основе возобновляемых геотермальных энергоресурсов прибрежных объектов (таких как порты и судостроительные предприятия) осложняются влиянием водной среды в непосредственной близости от объекта и фильтрацией вод в грунт. Все это определяет, с одной стороны, низкий потенциал источника тепла, с другой - значительно осложняет проектирование геотермальных систем. Прогноз долгосрочной энергетической эффективности геотермального теплового насоса, характеристики которого ухудшаются вследствие изменения температурного поля пласта при длительной эксплуатации скважины [1-4], актуален вследствие активного внедрения энергосберегающих технологий. Моделирование процессов изменения температуры грунтового массива в условиях знакопеременного нестационарного теплового потока представляет собой чрезвычайно сложную задачу, требует физико-математического описания сложного температурного поля, сформировавшегося под действием разнообразных нестационарных механизмов, определяющих процесс теплосбора (теплоотвода), включая внешние климатические условия, параметры теплового насоса, изменение характеристик грунта. Разработка долгосрочного прогноза изменения параметров геотермальной теплонасосной установки, учитывающего эти влияющие факторы, являются целью настоящей работы. Особенности разработки математической модели геотермальной скважины Основной целью исследования является разработка методики прогнозирования температурного поля в условиях работы различных технологических режимов независимо от климатического региона. Математическое моделирование температурного поля грунта при использовании геотермальной энергии сводится к решению задачи нестационарной теплопроводности [5]. В формировании температурного поля участвует большое количество процессов, среди них: - фоновый тепловой поток Земли, величина которого в каждом районе определяется целым рядом геофизических особенностей; - солнечное излучение и температура наружного воздуха, от которых зависит температура верхних слоев грунта; - теплофизические свойства грунта, различные для каждого региона. Значительная масса грунтового массива и малое значение коэффициента температуропроводности значительно сокращает зону влияния скважины и снижает темп охлаждения или нагрева пласта; - нестационарность режима работы скважины, обусловленная работой теплового насоса на конкретный объект, потребность которого в тепле или холоде зависит от множества факторов технологического или климатического происхождения; - возможность переменного по направлению теплового потока, например, чередующаяся работа теплового насоса на нагрев и охлаждение; - искажение температурного поля в связи с фильтрационным потоком грунтовых вод. В цилиндрической системе координат, характерной для вертикальной скважины, уравнение нестационарной теплопроводности [5], с учетом естественного фонового потока Земли в виде компенсирующей поправки - источника и стока тепла, - примет вид , где q - источники и стоки тепла, обусловленные тепловыми потоками земли и тепловыделения через поверхность, Вт/м2; r - расстояние от оси скважины, м; t - температура, К; a - коэффициент температуропроводности, м2/с; t - время, с. В настоящее время известно несколько технологических схем работы теплового насоса: - в одном режиме (только на тепло/охлаждение); - в двух режимах (сезонное: зимой на тепло, летом на охлаждение); - смешанное (одновременно на тепло и на охлаждение; встречается крайне редко и в условиях производственных процессов). В рамках представленной статьи данный режим не рассматривается. При разработке математической модели нестационарного режима работы теплового насоса выделим следующие особенности [6-9]. Во-первых, тепловой поток постоянно меняется по случайному закону, определяемому внешними факторами, включая климатические, технологические, особенности работы оборудования и пр. Во-вторых, грунт имеет значительную тепловую инерцию, сглаживая, как минимум, суточные колебания температуры. График действительного режима работы теплового насоса для каждого объекта и региона может различаться, поэтому заменим в модели действительную нагрузочную кривую на дискретные участки с постоянной, средней за расчетный период тепловой нагрузкой. На рис. 1 представлено графическое изображение принятых расчетных условий. Q, Вт t, мес. Рис. 1. Зависимость тепловых нагрузок от времени в условиях постановки решаемой задачи При условии работы в одном режиме - только на тепло (охлаждение) - краевые условия постановки задачи будут выглядеть следующим образом. Начальные условия для температуры пласта разделены на 3 этапа. Первый этап - до начала эксплуатации установки. Начальное условие - во всех точках температура равна фоновой температуре пласта. Начало координат находится также в центре скважины: ; где t0 - температура грунта (фоновая), °С; начальное положение температурного поля грунта, м; t - время, с; t - температура грунта, °С; Когда тепловой насос включается циклически, чередуются периоды подвода (отвода) тепла с периодами остановок. Для второго этапа - этапа остывания - начальными условиями будет сложное температурное поле, сформировавшееся при стационарной работе теплового насоса в течение расчетного периода 0 < t < t1: где τ1 - время остывания, с. Третий этап - включение теплового насоса после перерыва. Начальным условием является сформировавшееся ранее температурное поле за период остывания пласта t1 < t < t2: где τ2 - время включения после перерыва, с. Предельные состояния температуры пласта достигаются к моменту завершения цикла, например, к окончанию отопительного сезона, в связи с чем расчетные значения, используемые для определения эксплуатационных характеристик, фиксировались на этот момент времени. Обработка результатов Расчеты проводись на базе прикладной программы MatLab с чередованием периодов включения теплонасосной установки (ТНУ) и ее остановки, с цикличностью, определяемой временем года, в первую очередь при работе в одном технологическом режиме, во вторую очередь - для двух. Время достижения квазистационарного режима определялось расчетом и находится в пределах 5-ти лет. Диапазон плотностей тепловых потоков также определялся условиями эксплуатации действующих скважин и ограничивался 100-500 Вт/м2. На рис. 2-4 представлены результаты расчета температурного поля 1-го, 3-го и 5-го года эксплуатации скважины в одном технологическом режиме, т. е. без изменения направления теплового потока при его плотности на стенке обсадной колонны скважины, равной 100 Вт/м2. а б Рис. 2. Изменение температурного поля за 6 месяцев эксплуатации а б Рис. 3. Изменение температурного поля за 3 года эксплуатации а б) б Рис. 4. Изменение температурного поля за 5 лет эксплуатации Представленные результаты расчета отражают рост температуры грунтового пласта. Например, за 2 года работы теплового насоса в одном технологическом режиме исходная температура грунта в среднем увеличивается на 2 °С (диапазон зависит от величины активной тепловой нагрузки). В отключенном режиме ТНУ диапазон температур в среднем составляет 1,5-2 °С. Изображение полей 5-го года эксплуатации (рис. 4) демонстрирует увеличение границ температуры до 10 °С при работающем режиме (зависит от нагружаемой величины), 0,5-2 °С - в режиме «простоя» (когда тепловой насос не совершает никакой работы). Окончательно температура стабилизируется на 3-м году эксплуатации. Результаты прогноза температурного поля в условиях работы с изменением направления теплового потока (т. е. когда тепловой насос работает круглогодично) представлены на рис. 5-6. Длительное воздействие привело к росту температур в районе скважины на 2-2,5 °С. а б Рис. 5. Изменение температурного поля после 1-го года эксплуатации а б Рис. 6. Изменение температурного поля после 3-го года эксплуатации Результаты исследования позволяют утверждать, что стабилизация температурного поля при переменной активной нагрузке в течение круглогодичного эксплуатационного периода наступает уже на 2-й год. Температура пласта колеблется в диапазоне 10-15 °С и зависит от длительности эксплуатации ТНУ и от значений нагрузки. Температурное поле может искажаться в связи с фильтрационным потоком грунтовых вод. На рис. 7 приведен график-прогноз температурного поля в условиях влияния фильтрационного потока грунтовых вод. Смещение поля T оси T пик а б Рис. 7. Изменение температурного поля под действием влияния фильтрационных потоков грунтовых вод Смещение экстремума температурного поля от оси скважины наблюдается в диапазоне 0,1-0,4 м (рис. 7, б). На основании результатов исследований установлено, что стабилизация температурного поля грунта при условии влияния грунтовых вод наступает уже через 3,5 года. При длительности работы теплового насоса и малых значениях скорости фильтрации смещение по оси Х не наблюдается. Ввиду разнообразия параметров грунта результаты исследований необходимо представлять в обобщенном виде как зависимость безразмерной температуры от числа Фурье, что позволит применять разработанную методику к любым влияющим параметрам грунта в исследуемом диапазоне. Полученные результаты представлены в виде графиков (рис. 8). Рис. 8. Зависимость безразмерной температуры θ от критерия Фурье Fo при значениях безразмерного теплового потока Q Как видно из графиков, тепловая нагрузка хорошо аппроксимируется полиномом второго порядка. Для того чтобы прогнозирование изменения температуры грунта применить независимо от климатических особенностей региона и времени эксплуатации, на основании обработки полученных результатов получено следующее критериальное уравнение для работы теплового насоса без изменения направления теплового потока: (1) где - критерий Фурье; q - безразмерная температура; - безразмерный тепловой поток [10]. Стоит отметить, что уравнение (1) справедливо в диапазоне Fo (4672; 23360). Для работы теплового насоса с измененным направлением теплового потока будет поправка на коэффициент регенерации, и уравнение выглядит следующим образом: (2) где kp - коэффициент регенерации, а с поправками на фильтрационный поток: (3) где С - критерий коэффициента кратности водообмена (безразмерная величина, введенная в рамках исследовательской работы). Предложенные модифицированные критериальные уравнения (1)-(3), могут быть применимы для любого климатического региона и позволяют рассчитать изменение температурного поля во времени. Выводы Применение геотермальных тепловых насосов для энергоснабжения прибрежных объектов является перспективным направлением. Длительная эксплуатация геотермальных скважин портовых сооружений и судостроительных предприятий в зоне влияния поверхностных и грунтовых вод имеет свои особенности: - во-первых, в условиях работы теплового насоса без изменения направления теплового потока (без регенерации) установлено, что за первые 3 года температура грунтового массива повышается, а начиная с 3-го года стабилизируется; - во-вторых, в условиях работы теплового насоса с изменением направления теплового потока (с регенерацией) температурное поле грунта стабилизируется уже на 2-м году эксплуатации; - в-третьих, при малых значениях скорости фильтрации смещение присутствует, но на формировании поля отрицательно не сказывается. Стабилизация температурного поля грунта при условии влияния грунтовых вод наступает уже через 3,5 года; - в-четвертых, обработка полученных результатов в виде критериальных уравнений выполняется для возможности распространения результатов экспериментальных и теоретических исследований на широкий круг объектов проектирования.
Список литературы

1. Шишкин Н. Д., Просвирина И. С. Оценка эффективности применения теплонасосных установок в системах теплоснабжения Астраханской области // Известия АЖКХ. 2000. № 4. С. 7.

2. Васильев Г. П. Теплохладоснабжение зданий и сооружений с использованием низкопотенциальной тепловой энергии поверхностных слоев земли: дис. … д-ра техн. наук. М., 2006. 423 с.

3. Руденко Н. Н., Фурсова И. Н. Моделирование температурного поля в грунте. URL: http://ivdon.ru/magazine/archive/n4y2013/2137).

4. Крылов В. А., Черноозерский В. А., Никитин А. А., Баранов И. В. Учет неравномерности температурного поля в геотермальной скважине теплового насоса // Вестник MAX. 2015. № 1. С. 75-80.

5. Лыков А. В. Теория теплопроводности: учебн. пособ. М.: Высшая школа, 1967. 600 с.

6. Крылов Д. А. Математическое моделирование температурных полей с учетом фазовых переходов в криолитозоне // Наука и образование: электронное научно-техническое издание. 2012. Апрель. № 4. URL: http://technomag.bmstu.ru/doc/354740.html.

7. Федякин В. Я., Карпов М. К. Использование грунтовых теплообменников в системах теплоснабжения // Ползуновский вестник. 2006. № 4. С. 98-103.

8. Малых В. В., Удалов С. Н., Захаров А. А. Методика расчета грунтового аккумулятора // Материалы науч.-практ. конф. «Энерго- и ресурсоэффективность малоэтажных зданий» (Новосибирск, 19-20 марта 2013 г.). Новосибирск: Ин-т теплофизики им. С. С. Кутателадзе СО РАН. 2013. С. 317-318.

9. Костиков А. О., Харлампиди Д. Х. Влияние теплового состояния грунта на эффективность теплонаносной установки с грунтовым теплообменником // Энергетика: экономика, технология, экология. 2009. № 1. С. 32-40.

10. Сапрыкина Н. Ю., Яковлев П. В. Постановка задачи определения температурного поля геотермального теплового насоса на объектах нефтяной промышленности // Новейшие технологии освоения месторождений углеводородного сырья и обеспечения безопасности экосистем Каспийского шельфа: материалы VI Междунар. науч.-практ. конф. Астрахань: Изд-во Астрахан. гос. техн. ун-та, 2015. С. 117-125.