ENERGY SAVING TECHNOLOGIES OF PORT FACILITIES ON THE BASIS OF GEO-THERMAL HEAT PUMPS
Authors:
1.
Astrakhan Institute of Civil Engineering
2.
Saint-Petersburg Mining University
(Doctor of Technical Sciences, Professor; Professor of the Department of Thermal Engineering and Heat Power Engineering)
Type:
Article
Pages:
from 116
to 124
Status:
Published
Received:
20.02.2017
Accepted:
25.02.2017
Published:
25.02.2017
Subject area:
UDK 697.7
GRNTI 45.01 Общие вопросы электротехники
GRNTI 55.42 Двигателестроение
GRNTI 55.45 Судостроение
GRNTI 73.34 Водный транспорт
GRNTI 44.31 Теплоэнергетика. Теплотехника
GRNTI 45.01 Общие вопросы электротехники
GRNTI 55.42 Двигателестроение
GRNTI 55.45 Судостроение
GRNTI 73.34 Водный транспорт
GRNTI 44.31 Теплоэнергетика. Теплотехника
Language:
Russian
Keywords:
energy-saving technologies, coastal objects, heat pump, geothermal well, the tem-perature field of soil, mathematical model
Abstract (English):
Energy-saving technologies of port facilities and shipyards on the basis of geothermal heat pumps are a very promising direction. The features of offshore facilities are their presence in the zone of influence of seepage flow of ground water, significantly complicating the introduction of energysaving geothermal systems. The results of mathematical modeling of temperature field of soil mass, generated during the operation of the vertical ground heat exchanger of geothermal heat pump over 5 years are presented. The developed method is based on the equation of unsteady heat conduction and setting of boundary conditions for the different operating modes of the heat pump. The results of the temperature field changes were suggested in three modes of operation of geothermal plants: without changing the direction of heat flow, with a change in the direction of heat flow (regeneration) and in the condition of complicated influence of filtration flow of groundwater. The studies were conducted in order to build computational models that are adequate to real physical processes. The method of similarity theory is used to generalize the results.
Energy-saving technologies of port facilities and shipyards on the basis of geothermal heat pumps are a very promising direction. The features of offshore facilities are their presence in the zone of influence of seepage flow of ground water, significantly complicating the introduction of energysaving geothermal systems. The results of mathematical modeling of temperature field of soil mass, generated during the operation of the vertical ground heat exchanger of geothermal heat pump over 5 years are presented. The developed method is based on the equation of unsteady heat conduction and setting of boundary conditions for the different operating modes of the heat pump. The results of the temperature field changes were suggested in three modes of operation of geothermal plants: without changing the direction of heat flow, with a change in the direction of heat flow (regeneration) and in the condition of complicated influence of filtration flow of groundwater. The studies were conducted in order to build computational models that are adequate to real physical processes. The method of similarity theory is used to generalize the results.
Keywords:
energy-saving technologies, coastal objects, heat pump, geothermal well, the tem-perature field of soil, mathematical model
energy-saving technologies, coastal objects, heat pump, geothermal well, the tem-perature field of soil, mathematical model
Введение Применение энергосберегающих технологий на основе возобновляемых геотермальных энергоресурсов прибрежных объектов (таких как порты и судостроительные предприятия) осложняются влиянием водной среды в непосредственной близости от объекта и фильтрацией вод в грунт. Все это определяет, с одной стороны, низкий потенциал источника тепла, с другой - значительно осложняет проектирование геотермальных систем. Прогноз долгосрочной энергетической эффективности геотермального теплового насоса, характеристики которого ухудшаются вследствие изменения температурного поля пласта при длительной эксплуатации скважины [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. Shishkin N. D., Prosvirina I. S. Ocenka effektivnosti primeneniya teplonasosnyh ustanovok v sistemah teplosnabzheniya Astrahanskoy oblasti // Izvestiya AZhKH. 2000. № 4. S. 7.
2. Vasil'ev G. P. Teplohladosnabzhenie zdaniy i sooruzheniy s ispol'zovaniem nizkopotencial'noy teplovoy energii poverhnostnyh sloev zemli: dis. … d-ra tehn. nauk. M., 2006. 423 s.
3. Rudenko N. N., Fursova I. N. Modelirovanie temperaturnogo polya v grunte. URL: http://ivdon.ru/magazine/archive/n4y2013/2137).
4. Krylov V. A., Chernoozerskiy V. A., Nikitin A. A., Baranov I. V. Uchet neravnomernosti temperaturnogo polya v geotermal'noy skvazhine teplovogo nasosa // Vestnik MAX. 2015. № 1. S. 75-80.
5. Lykov A. V. Teoriya teploprovodnosti: uchebn. posob. M.: Vysshaya shkola, 1967. 600 s.
6. Krylov D. A. Matematicheskoe modelirovanie temperaturnyh poley s uchetom fazovyh perehodov v kriolitozone // Nauka i obrazovanie: elektronnoe nauchno-tehnicheskoe izdanie. 2012. Aprel'. № 4. URL: http://technomag.bmstu.ru/doc/354740.html.
7. Fedyakin V. Ya., Karpov M. K. Ispol'zovanie gruntovyh teploobmennikov v sistemah teplosnabzheniya // Polzunovskiy vestnik. 2006. № 4. S. 98-103.
8. Malyh V. V., Udalov S. N., Zaharov A. A. Metodika rascheta gruntovogo akkumulyatora // Materialy nauch.-prakt. konf. «Energo- i resursoeffektivnost' maloetazhnyh zdaniy» (Novosibirsk, 19-20 marta 2013 g.). Novosibirsk: In-t teplofiziki im. S. S. Kutateladze SO RAN. 2013. S. 317-318.
9. Kostikov A. O., Harlampidi D. H. Vliyanie teplovogo sostoyaniya grunta na effektivnost' teplonanosnoy ustanovki s gruntovym teploobmennikom // Energetika: ekonomika, tehnologiya, ekologiya. 2009. № 1. S. 32-40.
10. Saprykina N. Yu., Yakovlev P. V. Postanovka zadachi opredeleniya temperaturnogo polya geotermal'nogo teplovogo nasosa na ob'ektah neftyanoy promyshlennosti // Noveyshie tehnologii osvoeniya mestorozhdeniy uglevodorodnogo syr'ya i obespecheniya bezopasnosti ekosistem Kaspiyskogo shel'fa: materialy VI Mezhdunar. nauch.-prakt. konf. Astrahan': Izd-vo Astrahan. gos. tehn. un-ta, 2015. S. 117-125.
