Перейти к:
Способ оценки состояния кожухотрубчатых теплообменников с применением машинного обучения для повышения энергоэффективности
https://doi.org/10.23947/2687-1653-2026-26-2-2237
EDN: JTIHTA
Аннотация
Введение. Кожухотрубчатые теплообменники широко применяются на установках подготовки и переработки нефти, а также в системах теплоснабжения, где их надежная работа во многом определяется интенсивностью термического загрязнения поверхности труб. В связи с этим актуальной задачей является своевременное определение момента, когда теплообменному оборудованию требуется техническое обслуживание, поскольку накопление отложений приводит к снижению эффективности теплопередачи и росту эксплуатационных затрат. Анализ литературных источников показывает, что проблема прогнозирования загрязнения теплообменных аппаратов достаточно активно исследуется. Однако существующие подходы, ориентированные на учет режимных параметров и физико-химических свойств теплоносителей, не всегда применимы к специфическим условиям нефтепереработки, где состав сырья, режимы эксплуатации и характер отложений существенно отличаются от типовых промышленных процессов. Таким образом, сохраняется пробел в научных знаниях, связанный с недостаточной разработанностью моделей, позволяющих надежно описывать изменение термического сопротивления загрязнений и коэффициента теплопередачи именно для кожухотрубчатых теплообменников нефтеперерабатывающих установок. Целью настоящего исследования является определение зависимости коэффициентов термического сопротивления загрязнений и теплопередачи от режимных параметров — в частности, от выходных температур теплоносителей и плотности нефти. Для достижения поставленной цели предполагается применение методов машинного обучения, позволяющих построить более адекватные прогностические модели и тем самым повысить обоснованность решений по техническому обслуживанию теплообменного оборудования.
Материалы и методы. В качестве объекта исследования рассмотрен кожухотрубчатый теплообменник установки переработки нефти малой мощности, в котором циркуляционное орошение дизельного топлива (210–235 °С) выступает горячим теплоносителем в межтрубном пространстве, а нефть — холодным теплоносителем в трубном пучке. Для расчета определялись площади сечений трубного и межтрубного пространств, линейные скорости потоков, критерии Рейнольдса, Прандтля и Нуссельта, а также коэффициенты теплоотдачи со стороны горячего и холодного теплоносителей. Коэффициент теплопередачи рассчитывался с учетом термического сопротивления стенки и загрязнений, теплообмен — по уравнениям теплового баланса и теплопередачи. Для оценки термического сопротивления загрязнений и коэффициента теплопередачи использовались методы машинного обучения и символьной регрессии, в том числе PySR, SISSO, FROLS и CatBoostRegressor, на основе данных о плотности нефти и температуре теплоносителей. Поверочный расчет теплообменника выполнен методом последовательных приближений с решением системы нелинейных уравнений в Matlab.
Результаты исследования. Поверочный расчёт кожухотрубчатого теплообменника для нефти трёх месторождений показал, что при росте термического сопротивления загрязнений от 0 до 0,002 (м²·°С)/Вт коэффициент теплопередачи снижается примерно с 93–95 до 81–83 Вт/(м²·°С). Установлено, что изменение выходных температур теплоносителей служит информационной основой для идентификации загрязнений. Методами символьной регрессии SISSO и PySR получены аналитические зависимости термического сопротивления от плотности нефти и выходных температур с RMSE 1,25⋅10⁻⁸ и 2,34⋅10⁻⁸ (м²·°С)/Вт соответственно. Для прогнозирования динамики теплопередачи построены NARX-модели на базе FROLS и CatBoostRegressor, валидированные на промышленных данных. Ошибка постпрогноза для алгоритма CatBoostRegressor RMSE = 0,03573 Вт/(м²∙°С), для алгоритма FROLS RMSE = 0,01296 Вт/(м²∙°С).
Обсуждение. Снижение коэффициента теплопередачи на 13 % при росте термического сопротивления загрязнений до 0,002 (м²·°С)/Вт согласуется с теоретическими моделями и экспериментальными данными других исследователей. Применение методов SISSO и PySR обеспечило существенно более высокую точность по сравнению с нейросетевыми подходами при малых обучающих выборках. Алгоритм FROLS превзошёл CatBoostRegressor по точности постпрогноза динамики коэффициента теплопередачи, что объясняется компактностью полиномиальной модели и гладким характером исследуемого процесса. Установленный порог снижения коэффициента теплопередачи на 25 % соответствует отраслевым нормам и обеспечивает переход к обслуживанию оборудования по фактическому состоянию.
Заключение. Разработана методика определения термического сопротивления загрязнений теплообменника по выходной температуре теплоносителей и плотности нефти. Методами символьной регрессии SISSO и PySR получены аналитические зависимости с RMSE 1,25⋅10⁻⁸ (м²·°С)/Вт и 2,34⋅10⁻⁸ (м²·°С)/Вт. Построены NARX-модели прогнозирования коэффициента теплопередачи (RMSE 0,01296 и 0,03573 Вт/(м²·°С), валидированные на промышленных данных. Снижение коэффициента теплопередачи на 25,2 % принято критерием необходимости очистки аппарата, что обеспечивает переход к обслуживанию по фактическому состоянию и повышение энергоэффективности нефтепереработки.
Ключевые слова
Для цитирования:
Тугашова Л.Г., Затонский А.В. Способ оценки состояния кожухотрубчатых теплообменников с применением машинного обучения для повышения энергоэффективности. Advanced Engineering Research (Rostov-on-Don). 2026;26(2):2237. https://doi.org/10.23947/2687-1653-2026-26-2-2237. EDN: JTIHTA
For citation:
Tugashova L.G., Zatonskiy A.V. Machine Learning-Based Condition Assessment Method for Shell-and-Tube Heat Exchangers to Improve Energy Efficiency. Advanced Engineering Research (Rostov-on-Don). 2026;26(2):2237. https://doi.org/10.23947/2687-1653-2026-26-2-2237. EDN: JTIHTA
Введение. Повышение энергоэффективности процессов подготовки и переработки нефти является актуальной и практически значимой проблемой, особенно в условиях ужесточения требований к энергоэффективности и к снижению себестоимости различного рода производств при ограничениях на сбыт нефти и нефтепродуктов за рубежом, введенных в последнее время. Производители нефти не могут влиять на эти ограничения, но могут усовершенствовать внутренние процессы на предприятиях с целью снижения себестоимости производства. В качестве конкретных возможностей очевидны снижение теплопотерь, расхода тепловой энергии, повышение коэффициента полезного действия установок. Последнее возможно как за счет обеспечения регламентируемой температуры процесса, так и за счет диагностики оборудования для снижения вероятностей внеплановых простоев. Эти подходы тесно связаны с процессами теплообмена, имеющими большое значение в нефтехимии. Слой отложений загрязнений на поверхностях труб в 1÷3 мм может существенно снизить эффективность теплопередачи. При этом ориентировочные потери тепла составляют 20–40 %. Для компенсации снижения теплопередачи и поддержания необходимых параметров процесса следует увеличивать расход теплоносителя, что ведет к росту потребления топлива и электроэнергии. По мнению экспертов, загрязнение теплообменника в 1 мм может привести к перерасходу топлива на 10–15 %.
Однако непосредственное определение состояния кожухотрубчатых теплообменников нефтепродуктов затруднено их рабочими температурами и невозможностью наблюдать внутренние отложения, коррозию или деформацию непосредственно. Поэтому ключевыми параметрами контроля загрязнения стенок труб являются коэффициент теплопередачи, перепад давления, температура теплоносителей на выходе.
Причиной возникновения отложений на поверхностях стенок аппаратов является присутствие в сырье асфальтосмолистых веществ, парафинов, воды, солей. Для избавления от отложений и накипи в полостях нефтепроводов и теплообменников применяются мехaнические, тепловые и химические способы очистки. Очистка теплообменника происходит с помощью промывки его растворителем парафинов. На объектах подготовки и переработки нефти применяется ультразвуковая очистка от парафиновых отложений.
Отложения на стенках труб приводят к ухудшению процесса теплообмена и снижению коэффициента теплопередачи. Существуют различные способы диагностики и оценки слоя загрязнений на поверхности труб теплообменных аппаратов. Например, экспериментально с помощью систем мониторинга определяют характеристики чистого и загрязненного теплообменника. Другой подход связан с построением математических моделей. Примером является модель Эберта-Панчала и ее усовершенствованные варианты.
В статье [1] приведен обзор моделей и примеры количественной оценки толщины отложений накипи в пластинчатых теплообменниках систем теплоснабжения. Авторами [2] выполнен аналитический учет влияния загрязнений на технико-экономические показатели установок промысловой подготовки нефти. В работе экспериментальным методом получены коэффициенты теплопередачи, теплоотдачи, термического сопротивления отложений в теплообменниках. В работах [1][2] предложенные модели касаются конкретных типов теплообменников и отложений, они не прогнозируют динамику загрязнений.
Для решения задач диагностики теплообменного оборудования широко применяются программные средства, в частности программы технологического моделирования. Например, в работе [3] авторы рассматривают проблему загрязнения теплообменных аппаратов установок гидроочистки дизельного топлива. Задача решается с использованием моделирующей программы Hysys. В исследовании [4] численное решение осуществлено путем построения тепловой модели на основе метода конечных элементов в программе Ansys. В статье [5] авторами предложены упрощенные формулы для расчета показателя эффективной теплопроводности пористых структур с применением Mathcad, Ansys. На сегодня использование вышеприведенных программ не представляется возможным в рамках импортозамещения.
В статье [6] использован программный комплекс «МиР ПиА» для моделирования показателей качества подготовки нефти к транспортировке. Для работы необходимо приобретение лицензии.
Ряд исследований направлен на улучшение конструктивных параметров применяемых в отрасли аппаратов. Улучшение конструкции, совершенствование аэродинамики трубных пучков должны обеспечивать при высокой надежности интенсификацию в них теплообмена. В работе [7] это связывается с перепадами давления в теплообменнике, которые могут привести к сбоям в работе оборудования. Также предложено внедрение алгоритмов автоматической корректировки режимов работы. Авторами статьи [8] предложена нейросетевая модель прогнозирования теплогидравлической эффективности поверхностных интенсификаторов теплообмена. Для нейросетевых моделей необходимы большие объемы обучающих данных по загрязнениям теплообменников, которые сложно собрать. Кроме того, подобные решения, требующие изменения конструкции аппарата, могут быть неприменимы для эксплуатируемого оборудования.
Авторами статьи [9] с помощью методов имитационного моделирования получена модель кожухотрубчатого теплообменника и предложены способы усовершенствования системы управления теплообменным аппаратом. Математическая модель кожухотрубчатого испарителя представлена как объект с сосредоточенными параметрами. Составлен тепловой баланс в дифференциальной форме с учетом транспортного запаздывания.
В [10] представлена математическая модель кожухотрубчатого теплообменника в виде ARMAX-модели для участия в системе управления. В этих моделях не учитываются нелинейности.
В функционировании теплообменного оборудования для повышения эффективности тепло- и массообмена большую роль играет коллоидное состояние сырья и продуктов, особенно тяжелых фракций и остатков [11]. Такое решение актуально для тяжелых фракций, а также требует установки дополнительного оборудования.
В результате анализа изученной литературы установлено, что все еще остаются малоисследованными вопросы прогнозирования коэффициентов термического загрязнения и теплопередачи теплообменников с учетом физико-химических свойств сырья, а также применения методов, работающих с малым объемом данных и учитывающих нелинейности. Для аппроксимации и прогнозирования представляется интересным направление комбинирования подходов, например, сочетание физико-химических (балансовых) методов и машинного обучения.
Целью данного исследования является разработка методики определения коэффициента термического сопротивления загрязнений теплообменного аппарата и коэффициента теплопередачи в зависимости от выходных температур нефти и дизельного топлива, от плотности нефти, а также прогнозирование коэффициента теплопередачи.
Для достижения поставленной цели необходимо решить следующие задачи:
- выполнить поверочный тепловой расчет кожухотрубчатого теплообменника установки переработки нефти малой мощности, т. е. получить математическую модель теплообменника в балансовой форме;
- с помощью полученной выше модели кожухотрубчатого теплообменника исследовать влияние физико-химических свойств теплоносителей и коэффициента термического загрязнения на коэффициент теплопередачи и выходные параметры (на примере образцов нефти с различных месторождений);
- получить зависимость коэффициента термического сопротивления загрязнений и коэффициента теплопередачи от выходных температур нефти и дизельного топлива в теплообменниках, а также от плотности нефти, т.е. построить аппроксимирующие зависимости с помощью символьной регрессии;
- используя полученные аппроксимирующие зависимости, экспериментальные данные действующей установки, а также описание объекта в виде ранее полученных обыкновенных дифференциальных уравнений, построить нелинейную динамическую модель, пригодную для прогнозирования коэффициента теплопередачи.
Материалы и методы. Для определения коэффициента термического сопротивления загрязнений и коэффициента теплопередачи в кожухотрубчатых теплообменниках предлагается применить методы машинного обучения, использующие малые объемы выборки, например, символьную регрессию, так как в реальных условиях работы установки сложно получить большие объемы информации по загрязнению теплообменников. Для прогнозирования коэффициента теплопередачи предлагается NARX-модель.
Объектом исследования является кожухотрубчатый теплообменник установки переработки малой мощности. При расчете кожухотрубчатого теплообменника применялись исходные данные по конструктивным параметрам теплообменного аппарата: поверхность теплообмена — 68 м², длина корпуса теплообменника — 7,2 м, длина трубок — 6 м.
Суммарное термическое загрязнение поверхности труб принималось в диапазоне от 0 (чистый теплообменник) до 0,002 (м²∙°С)/Вт (загрязненный теплообменник).
Для подобных установок циркуляционное орошение (ЦО) организовано следующим образом. Дизельное топливо с выхода основной колонны разделяется на два потока: продуктовый и циркуляционный. ЦО, проходя через кожухотрубчатый теплообменник, возвращается в основную колонну. Температура ЦО составляет 210–235 °С. Дизельное топливо (межтрубное пространство) является горячим теплоносителем, нефть (трубный пучок) — холодным теплоносителем.
Первый шаг расчета — это вычисление площади сечения трубного и межтрубного пространств. Второй шаг — определение линейных скоростей движения нефти и дизельного топлива.
Режимы движения потоков теплоносителей и интенсивности теплообмена определяются значениями критериев Рейнольдса, Прандтля, Нуссельта. Следующий шаг — вычисление коэффициента теплоотдачи от горячего потока (дизельного топлива) к стенке трубы a1 и коэффициента теплоотдачи от поверхности труб к нагреваемому сырью (нефти) a2. Расчет теплообменного аппарата выполняется по стандартным методикам. В связи с этим общеизвестные формулы и подробности расчета теплообменного аппарата не приводятся.
Коэффициент теплопередачи определяется по формуле:
(1)
где dt — толщина стенки теплообменных труб, м; lt — коэффициент теплопроводности, Вт/(м∙°С); a1, a2 — коэффициенты теплоотдачи, Вт/(м²∙°С).
Коэффициент теплопередачи зависит от коэффициентов теплоотдачи в теплообменных трубах, термического сопротивления теплопроводности стенки трубки.
При наличии загрязнений внутри и снаружи теплообменных труб коэффициент теплопередачи определяется следующей формулой:
(2)
где Sr — суммарное термическое сопротивление слоя загрязнений, (м²∙°С)/Вт.
Уравнение теплопередачи выглядит следующим образом:
(3)
где K — коэффициент теплопередачи, Вт/(м²∙°С); F — поверхность теплообмена, м²; DТ — средний температурный напор, °С.
Уравнение теплового баланса кожухотрубчатого теплообменника:
(4)
где cLC — теплоемкость дизельного топлива ЦО, кДж/(кг∙°С); cn — теплоемкость нефти, кДж/(кг∙°С); GLC — расход дизельного топлива ЦО; Gn — расход нефти, кг/с; T13 — температура на тарелке № 13 выхода дизельного топлива, °С; TLC — температура возврата ЦО, °С; Tn1 — температура нефти на входе в теплообменник, °С; Tn2 — температура нефти на выходе из теплообменника, °С.
По исходным конструктивным данным выполнен поверочный расчет кожухотрубчатого теплообменника. Определены коэффициенты теплоотдачи, коэффициент теплопередачи, решена система уравнений (3) и (4).
На рис. 1 приведена блок-схема алгоритма решения уравнений (3) и (4) кожухотрубчатого теплообменника.

Рис. 1. Блок-схема алгоритма решения системы уравнений теплообменника
Температура дизельного топлива и нефти на выходе определяется с применением метода последовательных приближений. Для этого после задания исходных данных в блоке 1 необходимо выбрать первое приближение выходной температуры горячего теплоносителя (дизельного топлива).
В блоке 2 вычисляются теплоемкость нефти и дизельного топлива с учетом фракционного состава нефти. В блоке 3 выполняется вычисление выходной температуры дизельного топлива посредством решения уравнения теплового баланса теплообменного аппарата (4). В блоке 4 рассчитываются следующие параметры: средний температурный напор, коэффициент теплопередачи чистого теплообменника (1) и коэффициент теплопередачи с учетом коэффициента термического сопротивления слоя загрязнений (2).
В следующем блоке, 5, определяется расчетная площадь теплообмена из уравнения теплопередачи (3). В блоке 6 осуществляется верификация полученного значения расчетной площади. Существенное расхождение между полученным и заданным значениями свидетельствует о некорректном выборе температуры нефти и дизельного топлива. Итерационный расчет продолжается до достижения расхождения между расчетной и заданной площадями теплообмена не более 5 %. Решение системы нелинейных уравнений (3) и (4) реализовано с применением численных методов в среде Matlab.
Результаты расчета определяются при заданных значениях входных температур нефти и дизельного топлива и расходов нефти и дизельного топлива. Задаются свойства нефти (плотность, вязкость), ее фракционный состав по истинной температуре кипения (ИТК), конструктивные параметры теплообменника, тепловая нагрузка, термическое сопротивление слоя загрязнений поверхности труб.
На эффективность процесса теплообмена влияет термическое сопротивление загрязнений поверхности труб. Кроме этого, на выходные параметры, в частности на температуру теплоносителей на выходе, оказывают влияние состав и физико-химические свойства сырья. В данном исследовании использовались физико-химические свойства нефти трех месторождений (таблица 1).
Таблица 1
Физико-химические свойства нефти трех месторождений [12]
|
Нефть |
Плотность ρ, кг/м³ |
Выход фракций до 350 °С, % |
Вязкость η, мПа∙с |
Содержание асфальтенов, % |
|
Арланская |
891 |
42,8 |
35,37 |
5,80 |
|
Туймазинская |
856 |
53,4 |
10,18 |
3,40 |
|
Кулешовская |
819 |
63,6 |
3,29 |
0,70 |
В таблице 1 приведены физико-химические свойства нефти различных месторождений. Плотность и вязкость связаны с содержанием асфальтенов. При увеличении молекулярной массы, температурного предела выкипания фракций плотность и вязкость возрастают. При этом в нефти с меньшей плотностью выход светлых фракций до 350 °С будет больше. Плотность и вязкость нефти из приведенных в таблице 1 ее физико-химических свойств могут измеряться на потоке и в лаборатории. Для установления связи между термическим коэффициентом сопротивления загрязнений, коэффициентом теплопередачи, а также плотностью и вязкостью нефти также использовались данные таблицы 1.
Для определения термического сопротивления загрязнений предлагается следующая методика. С использованием математической модели теплообменного аппарата, приведенной в работе [13], и формул (1–4) определяются входные и выходные параметры объекта. Они формируют обучающую выборку. Для определения термического сопротивления загрязнений и коэффициента теплопередачи в кожухотрубчатых теплообменниках предлагается применять искусственные нейронные сети [14] и методы машинного обучения, связанные с символьной регрессией и идентификацией систем. Примерами методов символьной регрессии являются SISSO [15], TorchSISSO [16], PySR [17], SyMANTIC, RGP. Для получения зависимостей между параметрами теплообменного аппарата выбраны методы PySR (High-Performance Symbolic Regression) и SISSO (Sure Independence Screening and Sparsifying Operator). Они позволяют находить математические выражения, связывающие исходные данные. Приведенные методы выбраны в связи с тем, что они подходят для работы с малыми наборами данных. Ключевые гиперпараметры метода PySR: binary_operators — список бинарных операторов ['+', '-', '*', '/', '^'], unary_operators — список унарных операторов ['sin', 'cos', 'exp', 'log', 'sqrt' и др.], maxsize — максимальный размер выражения, niterations —– число итераций, loss_function — функция потерь, population_size — размер популяции. Ключевые гиперпараметры метода SISSO: n_expansion — число расширений признаков, n_term — число слагаемых в итоговой формуле, k — число признаков для L0-регуляризации, initial_screening — начальный отбор признаков, operators — набор математических операций для генерации признаков.
Для задач идентификации и прогнозирования за определенный период выбрана модель NARX (Nonlinear Autoregressive with Exogenous Inputs) — нелинейная авторегрессионная модель с внешними входами. Для определения ее структуры применялись алгоритмы FROLS (Forward Regression Orthogonal Least Squares) и CatBoostRegressor.
В представленной работе в качестве факторов принимались плотность нефти, температура горячего и холодного теплоносителей на выходе. Выходным параметром зависимости является термическое сопротивление загрязнений, или коэффициент теплопередачи. По математическому выражению, полученному с помощью методов PySR или SISSO, CatBoostRegressor или FROLS, и измеренным значениям технологических параметров вычисляется целевой параметр, позволяющий определять коэффициент теплопередачи и сроки технического обслуживания теплообменников.
Результаты исследования. При различных значениях термического загрязнения, плотности и вязкости нефти из разных источников (данные из таблицы 1) в таблице 2 представлены значения коэффициентов теплопередачи.
Таблица 2
Значения коэффициентов теплопередачи
|
Термическое сопротивление, (м²∙°С)/Вт |
Коэффициент теплопередачи чистого и загрязненного теплообменников, Вт/(м² °С) |
||
|
Арланская нефть |
Туймазинская нефть |
Кулешовская нефть |
|
|
0,00000 |
95,27486 |
94,24999 |
93,19691 |
|
0,00018 |
93,97257 |
92,97539 |
91,95044 |
|
0,00036 |
92,70541 |
91,73479 |
90,73687 |
|
0,00055 |
91,47196 |
90,52687 |
89,55492 |
|
0,00073 |
90,27091 |
89,35035 |
88,40336 |
|
0,00091 |
89,10098 |
88,20401 |
87,28104 |
|
0,00109 |
87,96099 |
87,08672 |
86,18686 |
|
0,00127 |
86,84981 |
85,99737 |
85,11978 |
|
0,00145 |
85,76635 |
84,93495 |
84,07879 |
|
0,00164 |
84,70959 |
83,89845 |
83,06296 |
|
0,00182 |
83,67855 |
82,88695 |
82,07138 |
|
0,00200 |
82,67231 |
81,89954 |
81,1032 |
На следующем шаге вычислены выходные температуры теплоносителей (нефти и дизельного топлива) в зависимости от термического сопротивления слоя загрязнений и физико-химических свойств нефти (плотности и вязкости). На рис. 2 и 3 приведены полученные результаты.

Рис. 2. Температура нефти на выходе

Рис. 3. Температура дизельного топлива на выходе
На графиках, приведенных на рис. 2 и 3, видно, что при наличии термического загрязнения поверхности труб эффективность теплообмена снижается. Температура нагреваемой нефти на выходе не достигает требуемого значения, соответствующего чистому теплообменнику. Для поддержания теплового потока на заданном уровне необходимо увеличить расход или температуру горячего теплоносителя.
Можно заметить, что при увеличении коэффициента термического загрязнения увеличиваются разность и соотношение между выходными температурами нефти и дизельного топлива. Эта закономерность учитывается при выводе зависимости коэффициента термического загрязнения методом SISSO. Также меняется величина отношения выходных температур загрязненного теплообменника к температурам чистого теплообменника. Эта закономерность учитывалась при применении метода PySR.
Для обучения на вход PySR и SISSO подаются значения плотности и температуры нефти на выходе и температура дизельного топлива на выходе. Целевым параметром является термическое сопротивление загрязнений. Инструментом реализации выбран программный пакет Python.
Примеры математических зависимостей для определения термического сопротивления загрязнений, полученных с помощью SISSO:
(5)
где x1 — плотность нефти, кг/м³; x2 — выходная температура дизельного топлива, °C; x3 — выходная температура нефти, °C; y — термическое сопротивление, (м²∙°С)/Вт; k1–k4 — коэффициенты.
С помощью PySR получена математическая зависимость для определения термического сопротивления загрязнений:
(6)
где x0 — плотность нефти, кг/м³; x1 — отношение выходной температуры дизельного топлива загрязненного теплообменника к выходной температуре дизельного топлива чистого теплообменника; x2 — отношение выходной температуры нефти загрязненного теплообменника к выходной температуре нефти чистого теплообменника; y — термическое сопротивление, (м²∙°С)/Вт; k1–k3 — коэффициенты.
Для сравнения точности полученных зависимостей использована среднеквадратическая ошибка RMSE, определяющая среднюю разницу между прогнозируемыми и исходными данными. RMSE для выражения (5) составляет 1,250371·e⁻⁰⁸ (м²∙°С)/Вт. Для выражения (6) RMSE = 2,336612·e⁻⁰⁸ (м²∙°С)/Вт.
Математические выражения (5, 6) используются для нахождения значений термического сопротивления загрязнений по измеренным параметрам действующей установки.
Для практического применения также необходимо знать величины термического сопротивления загрязнений и коэффициента теплопередачи, определенные не только в данный момент, но и за некоторый период. Для получения последовательности значений термического сопротивления загрязнений и коэффициента теплопередачи необходимы архивные данные за прошлые периоды.
Обучающая и тестовая выборки формировались с помощью динамической модели теплообменного аппарата, приведенной в работе [18], с добавлением в модель зависимостей (5) или (6). После обучения и тестирования моделей существует возможность прогнозирования целевого параметра.
Для моделирования и прогнозирования нескольких значений целевого параметра за определенный период использовалась библиотека Python SysIdentPy. Модель NARX подходит для моделирования динамических систем, где используются как текущие, так и предыдущие значения входов и выходов.
Для выбора структуры модели применены алгоритмы FROLS и CatBoostRegressor. Подобраны лаги входных (xlag) и выходных (ylag) переменных. Значения лагов определены методом перебора.
Для алгоритма FROLS выбраны параметры:
метод оценки параметров модели — LeastSquares и TotalLeastSquares;
- order_selection=True — автоматический подбор порядка модели (количества значимых термов);
- n_terms=12— максимальное количество термов в итоговой модели;
- model_type="NARMAX" — тип модели;
- xlag=[ 1,1,1] — лаги для входных переменных;
- ylag=12 — лаги для целевой переменной;
- Polynomial — базовый класс для создания полиномиальных функций с выбором степени полинома (degree=2);
- информационный критерий — AIC(Akaike Information Criterion).
Для алгоритма CatBoostRegressor (регрессор на основе градиентного бустинга) определяются следующие гиперпараметры:
- iterations=500 — количество деревьев в ансамбле;
- learning_rate=0,1 — скорость обучения (шаг градиентного спуска);
- depth=6 — максимальная глубина деревьев;
- xlag=[ 1,1,1];
- ylag=12.
В данном случае решается задача прогнозирования параметра, вследствие чего получаются результаты, которых до настоящего момента не было. В связи с этим метод постпрогноза является наиболее подходящим. Обработанные данные разделены на две части — для построения модели и постпрогноза в соотношении 80/20.
После обучения модели выполнена оценка ее качества с помощью постпрогноза. Ошибка постпрогноза для алгоритма CatBoostRegressor RMSE = 0,03573 Вт/(м²∙°С), для алгоритма FROLS RMSE = 0,01296 Вт/(м²∙°С). При формировании обучающей и тестовой выборок использовался диапазон изменения термического сопротивления загрязнений от 0 до 0,0022 (м²∙°С)/Вт, плотности нефти — от 810 до 910 кг/м³. Полученные результаты для плотности нефти 900 кг/м³ приведены на рис. 4.
Для валидации моделей использованы промышленные данные за 12 отсчетов. Временные ряды измеренных температур нефти и дизельного топлива обработаны с применением медианного фильтра. В качестве дополнительной информации использовались экспертные оценки технологов (в части дат плановых очисток теплообменников).

Рис. 4. Реальное и прогнозируемое значения коэффициента теплопередачи
Для чистого теплообменника коэффициент теплопередачи составляет 95 Вт/(м²∙°С), а через 12 отсчетов — 77 Вт/(м²∙°С). За это время коэффициент теплопередачи снизился на 18,61 %. Через 18 отсчетов (на рис. 4 не показано) коэффициент теплопередачи снизился на 25,2 %, что свидетельствует о необходимости проведения очистки оборудования для поддержания эффективности работы. Подобным образом можно выполнить прогноз целевого параметра в случае изменения плотности нефти.
Сравнивая значения коэффициента теплопередачи для чистого и загрязненного теплообменников за определенный период, можно сделать вывод о снижении эффективности теплообмена. При снижении параметра на 25 – 30 % формируется информация о необходимости очистки теплообменного аппарата.
Обсуждение. Результаты поверочного теплового расчёта кожухотрубчатого теплообменника установки переработки нефти малой мощности подтверждают корректность построенной математической модели. При увеличении суммарного термического сопротивления загрязнений от 0 до 0,002 (м²·°С)/Вт коэффициент теплопередачи снижается приблизительно на 13 – 14 % для всех трёх исследований нефти из различных источников.
Установлена зависимость эффективности теплообмена от физико-химических свойств нефти. При увеличении плотности коэффициент теплопередачи несколько возрастает вследствие более высокой теплоёмкости и теплопроводности. Вместе с тем более тяжёлая нефть содержат больше асфальтосмолистых веществ, ускоряющих нарастание загрязнений. Таким образом, влияние плотности носит двойственный характер и требует дальнейшего изучения.
Методы символьной регрессии PySR и SISSO показали высокую точность восстановления аналитических зависимостей: RMSE составили 1,25⋅10⁻⁸ и 2,34⋅10⁻⁸ (м²·°С)/Вт соответственно. Выражение, полученное методом SISSO, представляется более практически применимым, поскольку не требует эталонных значений температур чистого теплообменника.
В задаче прогнозирования алгоритм FROLS в составе модели NARX превзошёл CatBoostRegressor по точности (RMSE = 0,013 против 0,036 Вт/(м²·°С), что объясняется лучшей приспособленностью ортогональной регрессии к малым выборкам с гладкими динамическими зависимостями. Установленный критерий очистки (снижение коэффициента теплопередачи на 25–30 %) позволяет перейти от реактивной к проактивной стратегии технического обслуживания.
Заключение. В работе представлена методика определения коэффициента термического сопротивления загрязнений и коэффициента теплопередачи кожухотрубчатого теплообменника установки переработки нефти малой мощности на основе методов машинного обучения. В ходе исследования выполнен поверочный тепловой расчёт кожухотрубчатого теплообменника с использованием итерационного алгоритма решения системы уравнений теплопередачи и теплового баланса. Показано, что с ростом коэффициента термического загрязнения поверхности труб коэффициент теплопередачи снижается, а выходные температуры теплоносителей отклоняются от значений, соответствующих чистому аппарату. Установлено также влияние физико-химических свойств нефти — плотности и вязкости — на выходные параметры теплообменника.
С применением методов символьной регрессии PySR и SISSO получены аналитические зависимости коэффициента термического сопротивления загрязнений от выходных температур нефти и дизельного топлива, а также от плотности нефти. Среднеквадратическая ошибка полученных выражений не превышает 1,25⋅10⁻⁸ и 2,34⋅10⁻⁸ (м²·°С)/Вт, что свидетельствует об их высокой точности.
Для прогнозирования коэффициента теплопередачи во времени применены алгоритмы FROLS и CatBoostRegressor в рамках модели NARX. По результатам постпрогноза, ошибка для алгоритма FROLS составила RMSE = 0,013 Вт/(м²·°С), для CatBoostRegressor — RMSE = 0,036 Вт/(м²·°С). На основе динамики прогнозируемого коэффициента теплопередачи показано, что снижение данного параметра на 25–30 % относительно значения чистого теплообменника является обоснованным критерием для назначения технического обслуживания.
Валидация разработанных моделей выполнена с использованием промышленных данных и экспертных оценок технологов, что подтверждает применимость предложенной методики в реальных условиях эксплуатации.
Практическая значимость работы состоит в том, что разработанная методика позволяет по оперативно измеряемым технологическим параметрам — выходным температурам теплоносителей и плотности нефти — определять степень загрязнения теплообменного аппарата и прогнозировать сроки его очистки, что способствует повышению энергоэффективности и сокращению внеплановых простоев установок переработки нефти.
В перспективе предполагается расширение применения методики для других типов теплообменного оборудования, а также её интеграция в автоматизированные системы управления технологическими процессами.
Список литературы
1. Elistratova Y, Seminenko A, Minko V. Relevance of Contamination Models for Diagnostics of Plate Heat Exchangers. Bulletin of Belgorod State Technological University named after V.G. Shukhov. 2020;5(10):33–40. https://doi.org/10.34031/2071-7318-2020-5-10-33-40
2. Антипов А.И., Голубев Л.Г., Мухтаров Я.С. Влияние отложений в теплообменной аппаратуре на термические сопротивления и энергетические затраты на объектах промысловой подготовки нефти. Известия высших учебных заведений. Проблемы энергетики. 2003;(5–6):28–35.
3. Демин А.М., Науменко А.П. Диагностирование теплообменного оборудования на основе режимных параметров установки гидроочистки дизельных топлив. Омский научный вестник. 2019;166(4):84–88. https://doi.org/10.25206/1813-8225-2019-166-84-88
4. Позевалкин В.В., Поляков А.Н. Реализация цифровой модели тепловых характеристик на основе температурного поля. Advanced Engineering Research (Rostov-on-Don). 2024;24(2):178–189. https://doi.org/10.23947/2687-1653-2024-24-2-178-189
5. Губарева К.В., Еремин А.В. Численное решение задачи теплопроводности в пористой пластине с топологией трижды периодических минимальных поверхностей. Advanced Engineering Research (Rostov-on-Don). 2025;25(1):23–31. https://doi.org/10.23947/2687-1653-2025-25-1-23-31
6. Konygin SB, Konovalenko DV, Kryuchkov DA. Estimating the Effect of Instrument Accuracy on the Accuracy of Mathematical Modeling of Quality Indicators of Oil Treatment for Transportation. IOP Conference Series: Materials Science and Engineering. 2020;862(3):032068. https://doi.org/10.1088/1757-899X/862/3/032068
7. Martins T, Spengler AW, Oliveira JLG, de Paiva KV, Seman LO. Active Control System to Prevent Malfunctioning Caused by the Pressure Difference in Gasket Plate Heat Exchangers Applied in the Oil and Gas Industry. Sensors. 2022;22(12):4422. https://doi.org/10.3390/s22124422
8. Gil’fanov KK, Shakirov RA. Neural Network Modeling of Thermal Hydraulic Efficiency of Promising Surface Heat Transfer Intensifiers. Russian Aeronautics. 2021;64(1):61–70. https://doi.org/10.3103/S1068799821010086
9. Tugashova LG, Zatonskiy AV. Comparison of Methods for Heat Exchanger Control. Theoretical Foundations of Chemical Engineering. 2021;55(1):53–61. https://doi.org/10.1134/S0040579520060226
10. Sahoo A, Radhakrishnan TK, Sankar Rao C. Modeling and Control of a Real Time Shell and Tube Heat Exchanger. Resource-Efficient Technologies. 2017;3(1):124–132. https://doi.org/10.1016/j.reffit.2016.12.001
11. Глаголева О.Ф., Капустин В.М., Пискунов И.В., Усманов МР. Регулирование агрегативной устойчивости сырьевых смесей и товарных нефтепродуктов (обзор). Нефтехимия. 2020;60(5):577–585.
12. Капустин В.М., Рудин М.Г., Кукес С.Г. Справочник нефтепереработчика. Москва: Химия; 2018. 416 с.
13. Тугашова Л.Г. Повышение энергоэффективности и производительности на объектах топливно-энергетического комплекса с применением методов моделирования и оптимизации. Монография. Альметьевск: АГТУ ВШН; 2025. 100 с.
14. Дерябин И.А., Королев В.В., Сорокин Г.С. Использование рекуррентной нейронной сети для решения обратной задачи теплопроводности применительно к расчету температуры оборудования с ВВЭР. Известия вузов. Ядерная энергетика. 2024;(4):144–154. https://doi.org/10.26583/npe.2024.4.12
15. Xiaolin Jiang, Guanqi Liu, Jiaying Xie, Zhenpeng Hu. Boosting SISSO Performance on Small Sample Datasets by Using Random Forests Prescreening for Complex Feature Selection. arXiv:2409.19209
16. Muthyala M, Sorourifar F, Paulson JA. TorchSISSO: A PyTorch-Based Implementation of the Sure Independence Screening and Sparsifying Operator for Efficient and Interpretable Model Discovery. Digital Chemical Engineering. 2024;13:100198. https://doi.org/10.1016/j.dche.2024.100198
17. Tonda A. Review of PySR: High-Performance Symbolic Regression in Python and Julia. Genetic Programming and Evolvable Machines. 2025;26:7. https://doi.org/10.1007/s10710-024-09503-4
18. Затонский А.В., Тугашова Л.Г. Идентификация параметров динамической модели теплообменного аппарата. Вестник Пермского национального исследовательского политехнического университета. Электротехника, информационные технологии, системы управления. 2020;33:99–114. https://doi.org/10.15593/2224-9397/2020.1.06
Об авторах
Л. Г. ТугашоваРоссия
Лариса Геннадьевна Тугашова, кандидат технических наук, доцент кафедры «Автоматизация и информационные технологии»
423462, г. Альметьевск, ул. Советская, 186 а
ResearcherID: AAG-3459-2019
Scopus Author ID: 57200701887
SPIN-код: 4382-0940
А. В. Затонский
Россия
Андрей Владимирович Затонский, доктор технических наук, профессор, заведующий кафедрой «Автоматизация и информационные технологии» Березниковского филиала
614990, г. Пермь, пр. Комсомольский, 29
Scopus Author ID: 56069394200
SPIN-код: 3197-0372
Разработана методика определения загрязнения теплообменников по температурным параметрам теплоносителей и плотности нефти. Применены методы машинного обучения и символьной регрессии для получения высокоточных аналитических зависимостей термического сопротивления. Построены динамические модели прогнозирования коэффициента теплопередачи, валидированные на производственных данных нефтепереработки. Установлено, что критерием необходимости очистки оборудования служит снижение коэффициента теплопередачи на двадцать пять процентов по сравнению с начальными показателями. Методика обеспечивает переход к обслуживанию оборудования по его фактическому техническому состоянию, повышая энергоэффективность и экономическую целесообразность операций в нефтепереработке и системах теплоснабжения.
Рецензия
Для цитирования:
Тугашова Л.Г., Затонский А.В. Способ оценки состояния кожухотрубчатых теплообменников с применением машинного обучения для повышения энергоэффективности. Advanced Engineering Research (Rostov-on-Don). 2026;26(2):2237. https://doi.org/10.23947/2687-1653-2026-26-2-2237. EDN: JTIHTA
For citation:
Tugashova L.G., Zatonskiy A.V. Machine Learning-Based Condition Assessment Method for Shell-and-Tube Heat Exchangers to Improve Energy Efficiency. Advanced Engineering Research (Rostov-on-Don). 2026;26(2):2237. https://doi.org/10.23947/2687-1653-2026-26-2-2237. EDN: JTIHTA
JATS XML






































