2020
Химическая технология и биотехнология
№ 4
АВТОМАТИЗАЦИЯ ПРОЦЕССОВ ХИМИЧЕСКОЙ ТЕХНОЛОГИИ
Б01: 10.15593/2224-9400/2020.4.13 УДК 004.896
Е.А. Муравьева, М.Д. Квиринг
Уфимский государственный нефтяной технический университет, филиал в г. Стерлитамаке, Россия
РАЗРАБОТКА СИСТЕМЫ АДАПТИВНОГО НЕЧЕТКОГО УПРАВЛЕНИЯ ПЛАСТИНЧАТЫМИ ТЕПЛООБМЕННИКАМИ И ЕЕ МОДЕЛИРОВАНИЕ В СРЕДЕ МДТЬДБ БИМШМК
На сегодняшний день на многих предприятиях одной из стадий производства является нагрев жидкости. Для достижения необходимой температуры жидкости чаще всего используются пластинчатые теплообменники, благодаря своей высокой теплопередаче, надежности, меньшему объему по сравнению с другими аппаратами. Однако вследствие воздействия на объект управления внешних возмущающих воздействий происходит отклонение температуры жидкости, получаемой на выходе из теплобменника, от целевых значений. Для поддержания необходимой температуры жидкости без ошибок в системе теплообменников предлагается вместо РЮ-регуляторов использовать нечеткие регуляторы. Такие регуляторы способны в большей степени компенсировать возникающие возмущения и дать более точный результат конечного значения температуры. Данные нечеткие регуляторы способны повысить качество процесса регулирования, сократить время нагрева жидкости, учитывая внешние воздействия и взаимные влияния и адаптируясь к ним, что приводит к стабильному качеству производимой предприятием продукции.
В данной статье описана разработка адаптивного управления системой взаимосвязанных пластинчатых теплообменников на базе нечетких регуляторов. Представлена технологическая схема работы теплообменных аппаратов и описана их взаимосвязь. Также реализован способ расчета температуры жидкости на выходе из теплообменников в зависимости от заданных входных параметров. Проведено моделирование разработанной системы адаптивного нечеткого управления и сравнение полученных результатов с расчетными значениями.
Представлено описание процесса моделирования разработанной системы в среде МЛТЬЛБ. Система управления содержит четыре нечетких регулятора. Для работы каждого из них составлены правила нечеткого управления, в процессе которого решается задача нагрева жидкости до необходимой температуры.
Ключевые слова: теплообменник, управление, нечеткий регулятор, клапан, параметр, температура, давление.
E.A. Muravyova, M.D. Kviring
Ufa State Petroleum Technological University Branch of the University in the City of Sterlitamak, Russian Federation
DEVELOPMENT OF AN ADAPTIVE FUZZY CONTROL SYSTEM FOR PLATE HEAT EXCHANGERS AND ITS MODELING IN MATLAB SIMULINK
Today, in many enterprises, one of the stages ofproduction is heating the liquid. To achieve the required liquid temperature, plate heat exchangers are most often used, due to their high heat transfer, reliability, and smaller volume compared to other devices. However, due to the influence of external disturbances on the control object, the temperature of the liquid obtained at the outlet of the heat exchanger deviates from the target values. To maintain the required liquid temperature without errors in the heat exchanger system, it is proposed to use fuzzy controllers instead of PID controllers. Such regulators are able to compensate for the resulting disturbances to a greater extent, and get a more accurate result of the final temperature value. These fuzzy controllers are able to improve the quality of the control process, reduce the heating time of the liquid, taking into account external influences and mutual influences, and adapting to them, which leads to a stable quality of the products produced by the enterprise.
This article describes the development of adaptive control of a system of interconnected plate heat exchangers based on fuzzy controllers. The technological scheme of operation of heat exchangers is presented and their interrelation is described. A method for calculating the temperature of the liquid at the outlet of heat exchangers depending on the specified input parameters is also implemented. The developed adaptive fuzzy control system is modeled and the results are compared with the calculated values.
The process of modeling the developed system in the Matlab environment is described. The control system contains four fuzzy controls. For the operation of each of them, rules offuzzy control are drawn up in the process of which the problem of heating the liquid to the required temperature is solved.
Keywords: heat exchanger, fuzzy controller, parameter, temperature, valve, pressure, control.
Введение. В настоящее время промышленное производство обладает высокой сложностью. Применение большого количества объектов управления сопровождается различными внешними воздействиями, которые вызывают отклонение технологических параметров от заданных значений.
Большинство отраслей промышленности, таких как нефтеперерабатывающая, химическая, нефтехимическая, энергетическая, используют в своих производствах теплообменные аппараты, так как нагревание различных жидкостей до необходимой температуры - одна из
важнейших стадий на производстве. Самыми распространенными теп-лообменными аппаратами являются пластинчатые теплообменники, которые получили широкое применение благодаря своей высокой теплопередаче, надежности, меньшим по сравнению с другими аппаратами габаритам.
В большинстве случаев технологический процесс нагрева жидкости в теплообменниках осуществляется под управлением классических РГО-регуляторов, которые не всегда способны в полной мере отработать внешние воздействия на объект управления и взаимное влияние технологических параметров. Чаще всего это происходит из-за того, что РГО-регуляторы не всегда способны оптимальным образом компенсировать возникающие возмущения, следствием чего является нестабильность работы объекта управления (теплообменных аппаратов) и отклонение температуры жидкости от заданного значения на выходе из теплообменников [1, 2]. Таким образом, жидкость будет недостаточно нагретой или, напротив, перегретой, что может привести к возникновению брака на производстве.
Описание технологического объекта управления. Объектом управления является система из трех пластинчатых теплообменников. Разработанная технологическая схема взаимодействия теплообменных аппаратов представлена на рис. 1. Система теплообменных аппаратов содержит три пластинчатых теплообменника ТП1, ТП2 и ТП3. В качестве нагреваемой жидкости принята вода, в качестве теплоносителя -греющий пар.
Рис. 1. Технологическая схема работы системы пластинчатых теплообменников
Из коллектора пара греющий пар подается на каждый теплообменник. Холодная вода из соответствующего коллектора подается в каждый теплообменник насосами. Нагретая в теплообменниках вода поступает в общий коллектор и подается в следующий цех на технологические нужды.
В качестве регулируемой величины выбрана температура нагретой воды на выходе из каждого теплообменника, которые измеряются соответствующими датчиками температуры ТЕ1, ТЕ2 и ТЕ3.
На температуру нагреваемой воды влияют такие параметры, как давление греющего пара, расход холодной воды на входе в теплообменник и расход нагретой воды на выходе из теплообменника. Значения технологических параметров, идентичных для всех трех теплообменников, должны находиться в следующих диапазонах:
температура нагреваемой воды, °С.............................................70±5
давление греющего пара на входе
в теплообменник Рп, кПа...............................................от 200 до 300
расход холодной воды на входе
в теплообменник т/ч.....................................................от 4 до 7,5
расход нагретой воды на выходе
из теплообменника т/ч..................................................от 4 до 7,5
Для регулирования указанных выше параметров используются следующие исполнительные устройства: насосы с частотными преобразователями Чх, служащие для подачи холодной воды в теплообмен-ные аппараты, регулирующие клапаны в линиях подачи греющего пара ¥п, регулирующие клапаны в линиях отвода нагретой воды ¥т.
Система теплообменников ТП1-ТП3 функционирует следующим образом. Вода нагревается первым теплообменником ТП1. В ситуации, когда один теплообменник не способен нагреть воду до необходимой температуры, вода поступает на несколько теплообменников сразу, благодаря чему возможно достижение требуемой тепловой мощности за меньшее время. Это позволяет сократить потребление энергоресурсов за счет управления каждым теплообменником в составе общей системы [3-5].
Для осуществления процесса регулирования температуры в трех теплообменниках разработана система адаптивного управления с четырьмя нечеткими регуляторами НР1, НР2, НР, НР4 (рис. 2).
Задача регулятора НР1 заключается в оценке технологической ситуации, которую он осуществляет на основании значений давления пара и расходов холодной и нагретой воды. Результатом такой оценки является вывод о том, какой из нечетких регуляторов НР2-НР4 необходимо включить для эффективного управления пластинчатыми теплообменниками, а также сокращения энергоресурсов в процессе нагрева жидкости.
Рис. 2. Блок-схема нечеткого адаптивного управления
Регуляторы НР2-НР4 служат для регулирования входных параметров.
На вход нечеткого регулятора НР2 будут подаваться 4 параметра: Рп - давление пара в теплообменнике ТП1; - расход горячей воды из теплообменника ТП1; - расход холодной воды в теплообменник ТП1; Т - температура внутри теплообменника ТП1.
На вход нечеткого регулятора НР3 будут подаваться 8 параметров: Рп - давление пара в теплообменнике ТП1; - расход горячей воды из теплообменника ТП1; - расход холодной воды в теплообменник ТП1; Т - температура внутри теплообменника ТП1; Рп2 - давление пара в теплообменнике ТП2; Рг2 - расход горячей воды из теплообменника ТП2; Рх2 - расход холодной воды в теплообменник ТП2; Т2 - температура внутри теплообменника ТП2.
На вход нечеткого регулятора НР4 будут подаваться 12 параметров: Рп - давление пара в теплообменнике ТП1; - расход горячей
воды из теплообменника ТП1; Fx - расход холодной воды в теплообменник ТП1; T - температура внутри теплообменника ТП1; Рп2 - давление пара в теплообменнике ТП2; Fr2 - расход горячей воды из теплообменника ТП2; Fx2 - расход холодной воды в теплообменник ТП2; T2 - температура внутри теплообменника ТП2; Рп3 - давление пара в теплообменнике ТП2; Fr3 - расход горячей воды из теплообменника ТП2; Fk3 - расход холодной воды в теплообменник ТП2; T3 - температура внутри теплообменника ТП2.
Для реализации системы адаптивного управления разработана блок-схема, описывающая взаимодействие четырех нечетких регуляторов (см. рис. 2). В блок-схеме представлено взаимодействие 4 нечетких регуляторов. НР1 служит для анализа текущей ситуации в системе теплообменников, в зависимости от подачи давления пара, холодной воды и выхода нагретой жидкости. Для повышения эффективности производства в системе присутствуют три идентичных теплообменных аппарата. Исходя из поступивших на НР1 параметров, последний определяет, сколько теплообменников из трех требуется включить в работу.
Адаптивное управление призвано сократить количество правил в базе регуляторов НР2-НР4, повысив тем самым эффективность и точность управления системой из трех теплообменных аппаратов.
После блока «Параметры ТП1» на выходе имеет значение температуры в теплообменнике 1, полученной при регулировании параметров посредством клапанов ¥п, Уг и частотного преобразователя Чх.
Значение температуры после системы теплообменников определяется по формуле
V + V + Ч
^з—, (1)
где Vu - степени открытия клапанов пара V^V^ соответственно после теплообменников НР2-НР4; Vi- - степени открытия клапанов нагретой воды V^V^ соответственно после теплообменников НР2-НР4; Чх - это диапазон регулирования скорости частотными преобразователями Чх-Чх3 соответственно после теплообменников НР2-НР4.
Для управления системой теплообменников сигналы от датчиков поступают на НР1, который обрабатывает их и после анализа решает, какой НР нужно активизировать, чтобы поддержать заданные значения параметров.
Таким образом, сигнал поступает на один из трех НР. Если сигнал от НР1 поступает на НР2, который управляет параметрами тепло-
обменника ТП1, то в работу запускается только теплообменник ТП1. Если сигнал от НР1 поступает на НР3, который управляет параметрами теплообменников ТП1 и ТП2, то в работу запускаются теплообменники ТП1 и ТП2. Если сигнал от НР1 поступает на НР4, который управляет параметрами теплообменников ТП1, ТП2 и ТП3, то в работу запускаются все теплообменники: ТП1, ТП2 и ТП3.
Для достижения оптимальной температуры необходимо регулировать степень открытия клапанов пара Уп1-Упз, клапанов выхода нагретой воды Уг1-Угз, частоту преобразователей Чх1-Чх3. В итоге система адаптивного управления теплообменниками поддерживает значение температуры нагретой жидкости на уровне70 ± 5 °С.
Таким образом, совместная работа регуляторов НР1-НР4 позволяет поддерживать заданную температуру нагреваемой воды.
Разработка базы правил управления. Для реализации нечетких регуляторов НР1-НР4 сформулированы эвристические правила, описывающие работу системы управления. В качестве примера приводятся по одному правилу для каждого регулятора. Остальные правила составляются аналогичным способом.
Правило управления для НР1. Если давление подаваемого в теплообменник ТП1 пара высокое и предполагается повышение интенсивности подачи пара в диапазоне от 260 до 300 кПа и расход поступающей в теплообменник ТП1 холодной воды высокий и предполагается повышение интенсивности подачи холодной воды в диапазоне от 6 до 7,5 т/ч и расход нагретой воды на выходе из теплообменника ТП1 низкий и предполагается понижение интенсивности выхода горячей воды в диапазоне от 4 до 5,5 т/ч и температура выходящей из теплообменника жидкости высокая, то включен только НР2, который управляет теплообменником ТП1.
Правило управления для НР2. Если давление подаваемого в теплообменник ТП1 пара среднее в диапазоне от 230 до 270 кПа и расход поступающей жидкости средний в диапазоне от 5 до 6,5 т/ч и расход нагретой жидкости средний в диапазоне от 5 до 6,5 т/ч и температура выходящей из теплообменника жидкости средняя, то клапан пара Уп полуоткрыт в диапазоне от 50 до 85 % и клапан нагретой жидкости Уг полуоткрыт в диапазоне от 50 до 85 %, и скорость частотного преобразователя для двигателя насоса на подачу холодной воды теплообменника ТП1 работает на средней мощности Чх в диапазоне от 30 до 70 %.
Правило управления для НР3. Если давление пара Рп среднее в диапазоне от 230 до 270 кПа и расход поступающей жидкости ¥х низкий в диапазоне от 4 до 5,5 т/ч и расход нагретой жидкости ¥г высокий в диапазоне от 6 до 7,5 т/ч и температура выходящей из теплообменника жидкости низкая у первого теплообменника, и давление пара Рп2 среднее в диапазоне от 230 до 270 кПа и расход поступающей жидкости ¥х2 низкий в диапазоне от 4 до 5,5 т/ч и расход нагретой жидкости Бг2 высокий в диапазоне от 6 до 7,5 т/ч и температура выходящей из теплообменника жидкости низкая у второго теплообменника, то клапан пара Уп полуоткрыт в диапазоне от 50 до 85 % и клапан нагретой жидкости Уг открыт в диапазоне от 65 до 100 % скорость частотного преобразователя для двигателя насоса на подачу холодной воды теплообменника ТП1 работает на средней мощности Чх в диапазоне от 30 до 70 % и клапан пара Уп2 открыт в диапазоне от 65 до 100 %, и клапан нагретой жидкости Уг2 открыт в диапазоне от 65 до 100 % и скорость частотного преобразователя для двигателя насоса на подачу холодной воды теплообменника ТП2 работает на средней мощности Чх2 в диапазоне от 30 до 70 %.
Правило управления для НР4. Если давление пара Рп низкое в диапазоне от 200 до 240 кПа и расход поступающей жидкости ¥х высокий в диапазоне от 65 до 100 % и расход нагретой жидкости ¥г высокий в диапазоне от 65 до 100 % и температура выходящей из теплообменника жидкости низкая у первого теплообменника и давление пара Рп2 низкое в диапазоне от 200 до 240 кПа и расход поступающей жидкости ¥х2 высокий в диапазоне от 65 до 100 % и расход нагретой жидкости ¥г2 высокий в диапазоне от 65 до 100 % и температура выходящей из теплообменника жидкости низкая у второго теплообменника, и давление пара Рп3 низкое в диапазоне от 200 до 24, кПа и расход поступающей жидкости ¥х3 высокий в диапазоне от 65 до 100 % и расход нагретой жидкости ¥г3 высокий в диапазоне от 65 до 100 % и температура выходящей из теплообменника жидкости низкая у третьего теплообменника, то клапан пара Уп открыт в диапазоне от 65 до 100 % и клапан нагретой жидкости Уг полузакрыт в диапазоне от 25 до 60 % и скорость частотного преобразователя для двигателя насоса на подачу холодной воды теплообменника ТП1 работает на средней мощности Чх в диапазоне от 30 до 70 % и также клапан пара Уп2 открыт в диапазоне от 65 до 100 % и клапан нагретой жидкости Уг2 полузакрыт в диапазоне от 25 до 60 % и скорость частотного преобразователя для двигателя на-
соса на подачу холодной воды теплообменника ТП2 работает на средней мощности Чх2 в диапазоне от 30 до 70 % и также клапан пара Уп3 открыт в диапазоне от 65 до 100 %, и клапан нагретой жидкости Уг3 полузакрыт в диапазоне от 25 до 60 %, скорость частотного преобразователя для двигателя насоса на подачу холодной воды теплообменника ТП3 работает на средней мощности Чх3 в диапазоне от 30 до 70 %.
Для формирования базы правил нечеткого вывода определены следующие входные лингвистические переменные:
Рп - давление греющего пара на входе в теплообменник в диапазоне от 200 до 300 кПа;
иРп - интенсивность подачи пара в теплообменник в диапазоне от 108 000 до 180 000 м/ч;
- расход холодной воды, поступающей в теплообменник, в диапазоне от 20 до 60 об/с;
- интенсивность подачи холодной воды, поступающей в теплообменник в диапазоне от 3170 до 5952 м/ч;
- расход горячей воды, выходящей из теплообменника, в диапазоне от 20 до 60 об/с;
Црг - интенсивность горячей воды, выходящей из теплообменника, в диапазоне от 3170 до 5952 м/ч;
Т - температура жидкости, поступающей в теплообменник, в диапазоне от 60 до 80 °С.
В качестве выходных лингвистических переменных определены следующие:
НР2 - включен или выключен (1/0);
НР3 - включен или выключен (1/0);
НР4 - включен или выключен (1/0).
Чх - диапазон регулирования скорости частотным преобразователем для двигателя насоса на подачу холодной воды теплообменника в диапазоне от 20 до 60 об/с;
Уг - степень открытия клапана горячей воды теплообменника в диапазоне от 0 до 100 %;
Уп - степень открытия клапана пара теплообменника в диапазоне от 0 до 100 %;
Определены термы-множества для каждой лингвистической переменной:
Рп = {Рп1, Рп2, Рп3},
где Рп1 - «Давление пара = Низкое», где 200 < Рп1 < 240 кПа;
Рп2 - «Давление пара = Среднее», где 230 < Рп2< 270 кПа; Рп3 - «Давление пара = Высокое», где 260 < Рп3< 300 кПа.
¥х = {¥х1, ¥х2, ¥Х3},
где ¥х1 - «Расход жидкости = Низкий», где 4 < ¥х1 < 5,5 т/ч; ¥х2 - «Расход жидкости = Средний», где 5 < ¥х2 < 6,5 т/ч; ¥х3 - «Расход жидкости = Высокий», где 6 < ¥х3 < 7,5 т/ч.
¥г = {¥г1, ¥г2, ¥р3},
где ¥г1 - «Расход жидкости = Низкий», где 4 < ¥г1 < 5,5 т/ч; ¥г2 - «Расход жидкости = Средний», где 5 < ¥г2 < 6,5 т/ч; ¥г3 - «Расход жидкости = Высокий», где 6 < ¥г3 < 7,5 т/ч. Т = {71, 72, Т3},
где 71 - «Температура жидкости на входе в ТП1 = Низкая», где 60 < Т1 < 68 °С;
72 - «Температура жидкости на входе в ТП1 = Средняя», где 66 < Т2 < 74 °С;
73 - «Температура жидкости на входе в ТП1 = Средняя», где 72 < Т3 < 80 °С.
иРп = {ЦРп1, иРп2, иРп3},
где иРп1 - «Интенсивность подачи пара = Низкий», где 108 000 <
< иРп1 < 135000 т/ч;
иРп2 - «Интенсивность подачи пара = Средний», где 131 000 <
< иРп2 < 158 000 т/ч;
иРп3 - «Интенсивность подачи пара = Высокий», где 153 000<
< иРп3 < 180 000 т/ч.
и¥х = {и¥х1, и¥х2, и¥х3},
где и¥х1 - «Интенсивность подачи пара = Низкий», где 3170 <
< и¥х1 < 4359 т/ч;
и¥х2 - «Интенсивность подачи пара = Средний», где 3961 < и¥х2 <
< 5149 т/ч;
и¥х3 - «Интенсивность подачи пара = Высокий», где 4762 < и¥х3 <
< 5952 т/ч.
и¥г = {и¥Г1, и¥г2, и¥г3},
где и¥г1 - «Интенсивность подачи пара = Низкий», где 3170 <
< и¥г1 < 4359 т/ч;
и¥г2 - «Интенсивность подачи пара = Средний», где 3961 < и¥г2 <
< 5149 т/ч;
и¥г3 - «Интенсивность подачи пара = Высокий», где 4762 < и¥г3 <
< 5952 т/ч.
Уп = {Уп1, Уп2, Уп3, Уп4},
где Уп1 - «Клапан пара = Закрыт», где 0 < Уп1 < 35 %;
Уп2 - «Клапан пара = Полузакрыт», где 25 < Уп2 < 60 %;
Уп3 - «Клапан пара = Полуоткрыт», где 50 < Уп3 < 85 %;
Уп4 - «Клапан пара = Открыт», где 65 < Уп4 < 100 %.
Уг = {Уг1, Уг2, Уг3, Уг4},
где Уг1 - «Клапан пара = Закрыт», где 0 < Уг1 < 35 %;
Уг2 - «Клапан пара = Полузакрыт», где 25 < Уг2 < 60 %;
Уг3 - «Клапан пара = Полуоткрыт», где 50 < Уг3 < 85 %;
Уг4 - «Клапан пара = Открыт», где 65 < Уг4 < 100 %.
Чх = {ЧХ1, Чх2, Чх3},
где Чх1 - «Скорость частотного преобразователя для двигателя насоса = Низкая», где 0 < Чх1 < 40 %;
Чх2 - «Скорость частотного преобразователя для двигателя насоса = = Среднее», где 30 < Чх2 < 70 %;
Чх3 - «Скорость частотного преобразователя для двигателя насоса = Высокое», где 60 < Чх3 < 100 %.
Таким образом, рассмотренные выше в качестве примера эвристические правила будут иметь следующий вид.
Правило управления для НР1:
«ЕСЛИ Рп = Рп3 И иРп = иРп3 И Рх = Рх3 И ЦРх = ЦРхэ И ^ = Рг1 И и^г = и^г1 И Т=Т3, ТО НР = НР2».
Правило управления для НР2:
«ЕСЛИ Рп = Рп2 И Рх = Рх2 И Рг = Рг2 И Т=Т2, ТО Уп = Уп3 И Чх =
= Чх2 И Уг = Уг3».
Правило управления для НР3:
«ЕСЛИ Рп = Рп2 И Рх = Рх1 И Рг = Рг3 И Т=Т1 И Рп2 = Рп2 И Рх2 = = Рх1 И Рг2 = Рг3 И Т=Т1, ТО Уп = Уп2 И Чх = Чх2 И Уг = Уг2 И Уп2 = Уп2 И Чх2= Чх2 И Уг2 = Уг2».
Правило управления для НР4:
«ЕСЛИ Рп = Рп1 И Рх = Рх3 И Рг = Рг3 И Т = Т1 И Рп2 = Рп1 И Рх2 = = Рх3 И Рг2 = Рг3 И Т2=Т1 И Рп3 = Рп1 И Рх3 = Рх3 И Рг3 = Рг3 И Т?=Т1, ТО Уп = Уп4 И Чх = Чх2 И Уг = Уг2 И Уп2 = Уп4 И Чх2 = Чх2 И Уг2 = Уг2 И Уп3 = = Уп4 И Чх3 = Чх2 И Уг3 = Уг2».
Моделирование нечетких регуляторов в среде МАТЬАБ. Следующим этапом после составления базы правил является их реализация в составе нечетких регуляторов.
Для моделирования регулятора с использованием нечеткой логики использована среда МЛТЬЛВ. В качестве примера рассмотрено мо-
делирование нечеткого регулятора НР2. Число входных параметров регулятора НР2 - четыре, выходных - три. Для каждого входного и выходного параметра задан диапазон измерения. В соответствии с диапазонами заданы термы для входных параметров [6-9]. Задание термов для выходных параметров составлено аналогичным образом.
Также созданы модели нечетких регуляторов НР1, НР3 и НР4.
Разработка адаптивной нечеткой системы управления в среде MATLAB Simulink. Для реализации адаптивной системы нечеткого управления использована среда MATLAB Simulink с использованием спроектированного НР fuzzy-logic controller (нечеткий логический регулятор).
На примере приведенных выше правил управления для регуляторов НР1 и НР2 рассмотрена работа системы адаптивного нечеткого управления.
На первом этапе заданы исходные параметры: Рп = Рп3, F = ^п3, F = FX2, F = Fr2, T=T3.
Структурная схема модели формирования входных параметров регуляторов представлена на рис. 3. В модели использованы следующие типовые блоки: Slider Gain (ползунковый регулятор), Const (задание постоянных величин), Constant (задаются площади поперечных сечений), Mux (соединение сигналов), Display (отображение значения сигнала в виде числа), Subsyste (создание подсистем).
Рис. 3. Структура модели формирования исходных параметров
Сформированные моделью сигналы объединяются в блоке Mux и поступают на НР1, выходным сигналом является число теплообменников, которые необходимо включить. На его вход подаются следующие параметры: = Fti3 = 5000, Fx = Fx2 = 6.06, Fr = Fr2 = 6.08, Т = Т3 = = 78.89. На дисплее ТП1 отображена «1». Это означает, что теплообменник ТП1 включен, остальные теплообменники выключены. Данная работа регулятора соответствует правилу «ЕСЛИ Рп = Рп3 И иРп = = иРпз И Fx = FX2 И UFx = UFX2 И Fr = Fr2 И UFT = UF^ И T=T3, ТО НР = НР2». Результат такой работы регулятора НР1, полученный на модели, представлен на рис. 4.
Рис. 4. Структура модели блока вывода НР1
Следующим этапом моделирования работы системы является подача входных параметров на блок НР2, которые совпадают с входными параметрами для НР1. Нечеткий регулятор определяет, каким диапазонам и каким правилам соответствуют входные значения, и выдает значения степени открытия клапанов Уп, Уг, Ух и частоту преобразования для работы насосов Чх [10-14]. Значения, полученные в процессе моделирования работы регулятора НП2, представлены на рис. 4. Полученные значения управляющих сигналов означают следующее: Уп = = 42,5 - полузакрыт, Чх = 50 - средняя частота, Уг = 42,5 - полузакрыт. Полученные значения управляющих параметров находятся в диапазонах, обозначенных в следующем правиле для НР2: «ЕСЛИ Рп = Рп3 И Рх = Рх2 И Рг = Рг2 И Т=Т3, ТО Уп = Уп1 И Чх = Чх2 И Уг = Уг2» (рис. 5).
ж
Fuzzy Lsgic Ourtrul 1вг with RtilEviewEfl
HP2
1 bi. I
г"
-Jx
1
Рис. 5. Структура модели блока вывода НР2
По значениям входных и выходных параметров для НР2 выполняется расчет температуры жидкости на выходе из теплообменников. Алгоритм расчета содержит следующие шаги:
1. Расчет давления пара на выходе из теплообменников:
Р = (РЙ-РИМ!+Р, = (300 - 200>'425 + 200 = 242,5, (2) пк 100 п1 100
где Рпк - конечное давление пара; Рп2, Рп1 - конечное и начальное значения диапазона изменения давления, подаваемого в систему теплообменников пара, соответственно; Уп - степень открытия клапана в линии подвода пара.
2. Расчет отношения конечного давления пара к заданному:
Р 242 5
АР = =-— = 0,846, (3)
Рп 286,6
где ДР - отношение конечного давления пара к заданному; Рпк - конечное давление пара; Рп - начальное давление пара.
По аналогичным формулам рассчитываются значения отношений конечного расхода пара к заданному значению Д¥х и Д¥г.
3. Расчет расхода холодной жидкости на выходе из теплообменников:
¥„ = +¥х, = (7=5-42:50+4=
100 100
где ¥хк - конечное значение расхода холодной жидкости; ¥х2, ¥х1 - конечное и начальное значения диапазона изменения расхода холодной жидкости, подаваемого в систему теплообменников, соответственно; Чх - скорость частотного преобразователя на линии подачи холодной жидкости.
4. Расчет отношения конечного расхода холодной воды к заданному:
¥ 5 75
А¥ = -— = 0,9488, (5)
х ¥х 6,06
где Д¥х - отношение конечного значения расхода холодной жидкости к заданному; ¥хк - конечный расход холодной жидкости; ¥х - начальный расход холодной жидкости.
5. Расчет расхода горячей жидкости на выходе из теплообменников:
(F^ - Fri) • Ч:
(7.5 - 4) • 42,5 100
+ 4 = 5,4875, (6)
гк
100
где Ргк - конечное значение расхода горячей жидкости; Рг2, Рг1 - конечное и начальное значения диапазона изменения горячей жидкости, выходящей из теплообменника, соответственно; Уг - степень открытия клапана выходящей из теплообменника горячей жидкости.
6. Расчет отношения конечного расхода холодной жидкости к заданному:
где АРг - отношение конечного значения расхода горячей воды к заданному; Ргк - конечный расход горячей воды; Рг - начальный расход горячей воды.
Таким образом, получены следующие значения АР = 0,846, АРх = = 0,9488, АРг = 0,902.
7. Определение среднего арифметического найденных значений:
где п - число параметров, значения которых определяются в процессе работы системы адаптивного нечеткого управления.
Расчет конечной температуры горячей жидкости на выходе из теплообменников осуществляется согласно выражению
ТВЬКод = Ср.арифм(ДР, AF, AFx) • Т = 0,899 • 78,89 = 70,9, (9)
где Т - это заданная температура холодной жидкости на входе в теплообменник.
Полученное значение температуры 70,9 °С входит в необходимый диапазон технологических параметров, представленный в начале статьи.
На рис. 6 приведен фрагмент схемы модели нахождения значения конечной температуры горячей жидкости на выходе из теплообменника ТП3, реализующий формулу (9). Данная схема реализуется в блоке Subsystem, результат работы которого представлен на рис. 7.
AF =
F* = 5,4875 F 6,08
= 0,902,
(7)
Ср. арифм(ДР, AFT, ДРх) =
Рис. 6. Фрагмент схемы модели для расчета конечной температуры горячей жидкости на выходе из теплообменника ТП3
-to -гтп- 1п2 1пЗ In40ut1 1пБ
-► | 71.Э3531663715|
1пв ЛчЯ- Тв ыход
!—to
Subsystem
Рис. 7. Вывод конечной температуры горячей жидкости на выходе из теплообменника ТП3 на дисплее
Результаты проведения исследования для НР1 и НР2. После
сбора следующих данных: входных значений (Рп, иРп, ¥х, и¥х, ¥г, и¥г, 7) и выходных значений (Уп, Уг, Чх, Гвыход> в МЛТЬЛВ 81шиНпк, получены графики термов в правилах нечеткого вывода для каждого НР [15-20]. В качестве примера на рис. 8 представлены графики термов для регулятора НР1.
р _ 287 п.=6 06 Рд=6 08 Т = 78.9 Рп = 1.59е ♦005 ирх = 4.8 е+003 рг _
/ \ А N 1
/ 1 - А 1
/ 1 Л А А 1
Л Л 1А 1 1
|/\ А 1
|/\ 1А А 1
Л С2 |/\ Л, - 1
Л |/\ Л, А 1 1
Л |/\ Л, 1А А 1 1
Л |/\ Л, 1 1
Л С2 |Л 1 1 1 _1 1
— - 1 « 1 " '
-
Рис. 8. График термов для НР1
Обработка результатов исследования. Для определения зависимости выходного значения температуры на выходе из теплообменников Твь1ход от входной на входе в первый теплообменник ТП1 Т произведено девять измерений. В ходе эксперимента давление греющего пара, расход холодной воды и расход горячей воды оставались неименными; изменялись только значения входной температуры холодной воды, поступающей в первый теплообменник ТП1, в пределах от 72 до 80 °С. Результаты измерений представлены в таблице. Зависимость выходной температуры от входной представлена на рис. 9.
Результаты измерения температуры для НР2
Номер измерения Н Р2
Т Т 1 выход
1 72,8 65,46
2 73,46 66,05
3 74,38 66,87
4 75,94 68,29
5 76,87 69,1
6 77,79 69,94
7 78,34 70,44
8 79,08 71,1
9 79,9 71,84
На рис. 9 видно, что полученная зависимость линейная. Таким образом, при увеличении входной температуры пропорционально увеличивается температура на выходе, при этом значения выходной температуры находятся в заданном диапазоне. Зависимости выходных
температур от входных, получаемые при использовании остальных правил нечеткого вывода для НР2, НР3, НР4, аналогичны и также имеют линейный характер.
73 727!
70
^ВЫХОД 69 " 68676665-1-(-.-I-(-(-(-(-*
72 73 74 75 76 77 78 79 80 81 Т
Рис. 9. Зависимости выходной температуры от входной
Заключение. В ходе исследования была разработана система адаптивного управления теплообменными аппаратами на базе нечетких регуляторов, с помощью которой производилось регулирование температуры нагреваемой жидкости - воды.
Созданы базы правил нечеткого управления для четырех регуляторов. Проведено моделирование системы адаптивного управления нечеткими регуляторами температуры в среде МЛТЬЛВ и получены графики термов входных параметров для НР1: давления пара Рп, расхода холодной воды ¥х, расхода горячей воды ¥г, температуры на входе 7, интенсивности подачи пара иРп, интенсивности подачи холодной воды и¥х, интенсивности подачи горячей воды и¥г для каждого из трех теплообменников.
В результате моделирования системы адаптивного нечеткого управления теплообменными аппаратами в среде МЛТЬЛВ 81шиНпк установлено, что полученные в ходе работы системы значения температуры удовлетворяют заданным. Это свидетельствует о корректной работе разработанной системы.
Таким образом, можно сделать вывод, что разработанная адаптивная нечеткая система управления системой из трех теплообменных аппаратов является эффективной для производств, использующих теплообменники, так как повышается эффективность использования энергоресурсов за счет управления каждым теплообменником в составе общей системы.
Список литературы
1. Соловьев К. А., Муравьева Е.А., Султанов Р.Г. Адаптация нечеткого регулятора // Нефтегазовое дело. - 2014. - № 3. - С. 123-128.
2. Прядко П.В., Муравьева Е.А. Система управления сушильной печью для цеолитсодержащих катализаторов на основе 5-мерного нечеткого регулятора с интервальной неопределенностью // 63-я научно-техническая конференция студентов, аспирантов и молодых ученых УГНТУ. - Уфа, 2012. -С. 358.
3. Габитов Р.Ф., Каяшев А.И., Муравьева Е.А. Управление печью прокалки цеолитсодержащих катализаторов крекинга нефтепродуктов методом плавающего горизонта с использованием нейросетевой модели // Автоматизация, телемеханизация и связь в нефтяной промышленности. - 2010. -№ 12.- С. 12-19.
4. Разработка метода адаптации нечеткого регулятора, предназначенного для регулирования технологических параметров / Е.А. Муравьева, Д. С. Ра-дакина, М.И. Шарипов, Э.Р. Еникеева, Э.М. Артыкаева // Вестник Чувашского университета. - 2018. - № 1. - С. 137-145.
5. Jain M., Madeira A., Martins Manuel A. A fuzzy modal logic for fuzzy transition systems // Electronic Notes in Theoretical Computer Science. - 2020. -Vol. 3481. - P. 85-103.
6. Fuzzy controller design of micro-unmanned helicopter relying on improved genetic optimization algorithm / Y. Hu, Y. Yang, S. Li, Y. Zhou // Aerospace Science and Technology. - 2020. - Vol. 98. - № 105685.
7. J An improved genetic algorithm optimization fuzzy controller applied to the wellhead back pressure control system / H. Liang, J. Zou, K. Zuo, M. Khan // Mechanical Systems and Signal Processing. - 2020. - Vol. 142. - № 106708.
8. Talbi N. Design of Fuzzy Controller rule base using Bat Algorithm // Energy Procedia. - 2019. - Vol. 162. - P. 241-250.
9. Lamamra K., Batat F., Mokhtari F. A new technique with improved control quality of nonlinear systems using an optimized fuzzy logic controller // Expert Systems with Applications. - 2020. - Vol. 145, 1. - № 113148.
10. Alsafadi L.A., Chulin N.A., Mironova I.V. Synthesis of fuzzy controller based on simple PID controller // Procedia Computer Science. - 2019. -Vol. 150. - P. 28-38.
11. Sreedharan D., Paul V., Thottungal R. Mathematical modelling of polymer electrolyte membrane fuel cell and fuzzy-based intelligent controllers for performance enhancement // Computers & Electrical Engineering. - 2019. - Vol. 77. -P. 354-365.
12. Farajdadian S.M. Hassan Hosseini Design of an optimal fuzzy controller to obtain maximum power in solar power generation system // Solar Energy. -2019. - Vol. 182. - P. 161-178.
13. Назаров И.С. Математическое моделирование нечеткого регулятора // Информационные технологии в науке, управлении, социальной сфере и медицине: сб. науч. тр. V Междунар. науч. конф. - Томск, 2018. - С. 113-118.
14. Мефедова Ю.А., Моисеева Ю.Р. Применение нечеткой логики в системах управления // Актуальные проблемы и пути развития энергетики, техники и технологий: сб. тр. IV Междунар. науч.-практ. конф. - СПб., 2018. -С.94-97.
15. Бузикаева А.В., Черный С.П. Особенности моделирования нечетких регуляторов для систем электропривода переменного тока // Электропривод на транспорте и в промышленности: тр. II Всерос. науч.-практ. конф. -Хабаровск, 2018. - С. 271-275.
16. Каримов И.И. Концептуальная модель нечеткого регулятора для определения возмущений // Современные технологии в нефтегазовом деле -2019: сб. тр. междунар. науч.-техн. конф.: в 2 т. - Октябрьский, 2019. -С.300-302.
17. Алимова М.М. Нечеткая модель управления технологическими параметрами теплогенерирующих установок // Интернаука. - 2019. - № 8 (90). -С.34-35.
18. Карабцов Р. Д. Оптимизация нечетких систем управления с использованием генетического алгоритма // Прикладная математика и фундаментальная информатика. - 2018. - Т. 5, № 2. - С. 43-50.
19. Шилкина С.В., Фокина Е.Н. Контроллер нечеткой логики в управлении технологическими процессами // Вестник Сибирского государственного автомобильно-дорожного университета. - 2018. - Т. 15, № 1 (59). - С. 106-114.
20. Потапов Н.С., Маслов А.А. Оценка адекватности модели адаптивной системы управления с применением средств математического моделирования // Наука - производству: материалы междунар. науч.-практ. конф. / Мурманский гос. техн. ун-т. - Мурманск, 2018. - С. 77-84.
References
1. Solov'ev K.A., Murav'eva E.A., Sultanov R.G. Adaptacija nechetkogo reguljatora [Adaptation of a fuzzy controller]. Oil and Gas business, 2014, no 3, pp.123-128.
2. Prjadko P.V., Murav'eva E.A. Sistema upravlenija sushil'noj pech'ju dlja ceolitsoderzhashhih katalizatorov na osnove 5-mernogo nechetkogo reguljatora s interval'noj neopredelennost'ju [Control system of a drying oven for zeolite-containing catalysts based on a 5-dimensional fuzzy controller with interval uncertainty]. 63rd scientific and technical conference of students, postgraduates and young scientists of USNTU, 2012, pp. 358.
3. Gabitov R.F., Kajashev A.I., Murav'eva E.A. Upravlenie pech'ju prokalki ceolitsoderzhashhih katalizatorov krekinga nefteproduktov metodom plavaju-shhego gorizonta s ispol'zovaniem nejrosetevoj modeli [Control of the furnace for calcining zeolite-containing catalysts for cracking petroleum products by the floating horizon method using a neural network model]. Automation, telemechanization and communication in the oil industry, 2010, no 12, pp. 12-19.
4. Murav'eva E.A., Radakina D.S., Sharipov M.I., Enikeeva Je.R., Artykae-va Je.M. Razrabotka metoda adaptacii nechetkogo reguljatora, prednaznachennogo dlja regulirovanija tehnologicheskih parametrov [Development of a method for adapting a fuzzy controller designed for regulating technological parameters]. Bulletin of the Chuvash University, 2018, no 1, pp. 137-145.
5. Jain M., Madeira A., Martins Manuel A. A Fuzzy Modal Logic for Fuzzy Transition Systems. Electronic Notes in Theoretical Computer Science, 2020, vol. 3481, pp.85-103.
6. Hu Y., Yang Y., Li S., Zhou Y. Fuzzy controller design of micro-unmanned helicopter relying on improved genetic optimization algorithm. Erospace Science and Technology, 2020, vol. 98, no 105685.
7. Liang H., Zou J., Zuo K., Khan M. J An improved genetic algorithm optimization fuzzy controller applied to the wellhead back pressure control system. Mechanical Systems and Signal Processing, 2020, vol. 142, no 106708.
8. Talbi N. Design of Fuzzy Controller rule base using Bat Algorithm. Energy Procedia, 2019, vol. 162, pp. 241-250.
9. Lamamra K., Batat F., Mokhtari F. A new technique with improved control quality of nonlinear systems using an optimized fuzzy logic controller. Expert Systems with Applications, 2020, vol. 145, 1. no 113148.
10. Alsafadi L. A., Chulin N. A., Mironova I. V. Synthesis of Fuzzy Controller Based on Simple PID Controller. Procedia Computer Science, 2019, vol. 150, pp. 28-38.
11. Sreedharan D., Paul V., Thottungal R. Mathematical modelling of polymer electrolyte membrane fuel cell and fuzzy-based intelligent controllers for performance enhancement. Computers & Electrical Engineering, 2019, vol. 77, pp. 354-365.
12. Farajdadian, S. M. Hassan Hosseini Design of an optimal fuzzy controller to obtain maximum power in solar power generation system. Solar Energy, 2019, vol. 182, pp. 161-178.
13. Nazarov I.S. Matematicheskoe modelirovanie nechetkogo reguljatora [Mathematical modeling of a fuzzy controller]. In the collection of scientific papers of the V International scientific conference: Information technologies in science, management, social sphere and medicine, 2018, pp. 113-118.
14. Mefedova Ju.A., Moiseeva Ju.R. Primenenie nechetkoj logiki v sistemah upravlenija [Application of fuzzy logic in control systems]. Proceedings of the IV International scientific and practical conference: Actual problems and ways of development of energy, engineering and technologies, 2018, pp. 94-97.
15. Buzikaeva A.V., Chernyj S.P. Osobennosti modelirovanija nechetkih reguljatorov dlja sistem jelektroprivoda peremennogo toka [Features of modeling fuzzy controllers for alternating current electric drive systems]. SB.: electric Drive in transport and industry proceedings of the II all-Russian scientific and practical conference, 2018, pp. 271-275.
16. Karimov I.I. Konceptual'naja model' nechetkogo reguljatora dlja opredelenija vozmushhenij [Conceptual model of a fuzzy controller for determining perturbations]. In the collection: Modern technologies in oil and gas business-
2019 proceedings of the international scientific and technical conference: in 2 volumes, 2019, pp. 300-302.
17. Alimova M.M. Nechetkaja model' upravlenija tehnologicheskimi parametrami teplogenerirujushhih ustanovok [Fuzzy model of control of technological parameters of heat generating plants]. Internauka, 2019, no 8 (90), pp. 34-35.
18. Karabcov R.D. Optimizacija nechjotkih sistem upravlenija s ispol'zo-vaniem geneticheskogo algoritma [Optimization of fuzzy control systems using a genetic algorithm]. Applied mathematics and fundamental computer science, 2018, vol. 5, no 2, pp. 43-50.
19. Shilkina S.V., Fokina E.N. Kontroller nechjotkoj logiki v upravlenii tehnologicheskimi processami [Controller of fuzzy logic in process control]. Bulletin of the Siberian state automobile and road University, 2018, vol. 15. no 1 (59), pp. 106-114.
20. Potapov N.S., Maslov A.A. Ocenka adekvatnosti modeli adaptivnoj sistemy upravlenija s primeneniem sredstv matematicheskogo modelirovanija [Assessment of the adequacy of the adaptive control system model using mathematical modeling tools], Nauka - proizvodstvu. Materialy mezhdunarodnoi konferentsii. Murmansk, 2018, pp. 77-84.
Получено 05.09.2020
Об авторах
Муравьева Елена Александровна (Стерлитамак, Россия) - доктор технических наук, профессор, заведующий кафедрой «Автоматизированные технологические и информационные системы» Уфимского государственного нефтяного технического университета, филиал в г. Стерлитамаке (453118, г. Стерлитамак, пр. Октября, 2, e-mail: [email protected]).
Квиринг Маргарита Дмитриевна (Стерлитамак, Россия) - магистрант кафедры «Автоматизированные технологические и информационные системы» Уфимского государственного нефтяного технического университета, филиал в г. Стерлитамаке (453118, г. Стерлитамак, пр. Октября, 2, e-mail: [email protected]).
About the authors
Еlena A. Muravyova (Sterlitamak, Russian Federation) - Doctor of Technical Sciences, Professor, Head of the Department "Automated technological and information systems" of Ufa State oil Technical University, Branch in Sterlitamak (2, October av., Sterlitamak, 453118, e-mail: [email protected]).
Margarita D. Kviring (Sterlitamak, Russian Federation) - Undergraduate Student of the Department "Automated technological and information systems" of the Ufa State oil Technical University, Branch in Sterlitamak (2, October av., Sterlitamak, 453118, e-mail: [email protected]).