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

Ключевые слова:
модель, анализ, холодильная установка, идентификация эксплуатационных параметров
Текст
Введение Системному анализу технологических процессов и идентификации эксплуатационных и аварийных ситуаций с использованием новых компьютерных технологий и методов искусственного интеллекта посвящены многие работы в связи с актуальностью данного вопроса. Однако они не учитывают особенности технологического процесса, связанные с влиянием параметров окружающей среды. Так, А. В. Филоненко [1] предлагает модель идентификации эксплуатационных ситуаций установки получения серы методом Клауса, в которой оценка ситуации производится на основании значений параметров технологического процесса. Для идентификации причин эксплуатационных ситуаций на объекте О. М. Проталинский [2] рассматривает лингвистическую переменную «Причина ситуации». Число ситуационных математических моделей в этом случае соответствует числу причин, которые могут привести к аварийной ситуации. В названных работах предлагается для всех входных и выходных координат использовать лингвистические переменные, представленные в виде нечетких множеств. Обработка информации производится с помощью математического аппарата нечетких множеств. А. Ф. Егоров [3] разработал модель, ориентированную на своевременное устранение неисправностей. Разработке системы управления аммиачной холодильной установки, учитывающей изменения параметров окружающей среды и обеспечивающей тем самым безопасность и экономичность эксплуатации, посвящена данная работа. Разработка системы управления Для отработки принципов идентификации эксплуатационных параметров с последующей подачей управляющего воздействия выбрана небольшая (сравнительно) тепловая машина, эксплуатационные параметры которой в значительной мере зависят от окружающей среды – аммиачная холодильная установка. На основе исследования литературы составлена блок-схема управления (рис. 1), режим управления реализуется простыми методами на основании комплексов «Да – Нет», «Больше – Меньше», «Если – то». Рис. 1. Блок-схема управления холодильной машиной Холодильная установка – это сложный комплекс оборудования, работающий в изменяющемся режиме, и для выбора определенного управляющего воздействия необходимо провести ряд технических расчетов, учитывающих изменение наружных условий. Разработанная блок-схема управления холодильной установкой представлена на рис. 2. Кроме того, за время эксплуатации технической системы меняется цена энергоносителей, и, таким образом, проектные значения разности температуры становятся нерациональными. Поэтому необходимо периодически проводить расчет рациональной разности значений температуры в теплообменных аппаратах с введением в исходные данные новых значений цен, а также остаточной стоимости оборудования. Эксплуатация холодильной системы в течение года должна вестись в рациональном режиме, а для этого необходимо анализировать годовой график эксплуатации. Другими словами, для выдачи управляющего воздействия по измеренным параметрам работы достаточно коротких расчетов, блок-схема которых приведена на рис. 2. После изменения расценок энергоносителей достаточно провести один расчет (например, в год) по определению рациональной разности значений температуры и внести их в базу рекомендуемых значений расчета. Термоэкономический анализ является наиболее совершенным по сравнению с другими методами определения эффективности энергопреобразующих систем. Термоэкономику рассматривают как метод, способный указать пути сокращения стоимости эксплуатации системы. При этом, как показывает практика проведения термоэкономического анализа [4], получить подобную информацию другими методами анализа невозможно. Рис. 2. Управление в режиме изменения входных параметров С целью создания экспертной системы управления холодильной установкой разработан алгоритм термоэкономического анализа одноступенчатой одноцелевой аммиачной системы, принципиальная схема которой приведена на рис. 3. Рис. 3. Принципиальная схема холодильной установки: а – окр. ср – окружающая среда; кр – компрессор; гр – градирня; к – конденсатор; Нвод – насос воды; Нтн – насос теплоносителя; и – испаритель; р.в – регулирующий вентиль; вент – вентилятор; во – воздухоохладитель; б – 1, 2, 3, 4 – характерные точки цикла Исходными данными являются: 1. Проектные и эксплуатационные характеристики: Vh – объем, описываемый поршнями компрессора, м3/с; Fк – площадь теплообменной поверхности конденсатора, м2; kк – коэффициент теплопередачи конденсатора, Вт/(м2·К); Fи – площадь теплообменной поверхности испарителя, м2; kи – коэффициент теплопередачи испарителя, м2; Сэ – стоимость электроэнергии, руб./(кВтч); Св – стоимость воды, руб. м3; Сэкс – стоимость эксплуатации и обслуживания, руб. ч; Сд – коэффициент дисконтирования и другие отчисления и коэффициенты, включая банковский процент и др.; Соб – стоимость здания и оборудования для производства холода, руб.·год. 2. Параметры, измеряемые в соответствии с требованиями ПБ 09-595-03 [5]: Ро – давление всасывания (кипения), МПа; Рк – давление конденсации, МПа; tвс – температура пара, всасываемого компрессором, оС; tн – температура пара на выходе из компрессора, оС; tж – температура жидкого хладагента на выходе из конденсатора, оС; tw1 – температура воды на входе в конденсатор, оС; tw2 – температура воды на выходе из конденсатора, оС; tокр. ср – температура наружного воздуха, оС; tкам – средняя температура по объему камеры, оС; – влажность воздуха, %; Nкр – мощность, затрачиваемая приводом компрессора, кВт·ч; Nвент.во – мощность, затрачиваемая приводом вентиляторов воздухоохладителей, кВт·ч; Nн. тн – мощность, затраченная приводом насоса теплоносителя, кВтч; Nн. воды – мощность, затраченная насосом воды, кВтч; N вент. гр – мощность, затраченная вентилятором градирни, кВтч. 3. Параметры (и характеристики), получаемые на основе измеренных величин: То – температура кипения, получаемая по давлению кипения (всасывания) по программе расчета свойств хладагента, К; Рк – давление конденсации, МПа, получаемое по программе расчета свойств хладагента – давление насыщения при температуре Тк = tw2 + (2÷4) + 273, К; Тк – температура конденсации (насыщения), получаемая по давлению конденсации (Рк) по программе расчета свойств хладагента, К; Токр.ср = tокр.ср + 273 – температура наружного воздуха, К; Ткам = tкам + 273 – средняя температура по объему камеры, ºС; tвс = (tвс + 273) – То – перегрев пара на всасывании компрессора, К; Тн = tн + 273 – температура нагнетания по шкале Кельвина, К; Т3 = tж + 273 – температура жидкого хладагента в т. 3; i1 – энтальпия пара холодильного агента в т. 1 цикла, кДж/кг, вычисляется по программе расчета свойств хладагента на основании ввода давления Ро и температуры Твс; S1 – энтропия пара хладагента в т. 1 цикла, кДж/(кгК), вычисляется по программе расчета свойств хладагента на основании ввода давления Ро и температуры Твс; i1’ – энтальпия пара на выходе из испарителя, кДж/кг, вычисляется по программе расчета свойств хладагента на основании ввода давления Ро и температуры Т1’ = То + + tвс /2; S1’ – энтропия пара хладагента в т. 3’ цикла, кДж/(кгК), вычисляется по программе расчета свойств хладагента на основании ввода давления Ро и температуры Т1 = То + tвс/2; i2 – энтальпия пара хладагента в конце сжатия (т. 2 цикла), вычисляется по программе расчета свойств хладагента на основании ввода давления Рк и температуры Тн, кДж/кг; S2 – энтропия пара хладагента в т. 2 цикла, кДж/(кгК), вычисляется по программе расчета свойств хладагента на основании ввода давления Р2 и температуры Т2; i3’ – энтальпия жидкого хладагента на линии насыщения (т. 3’ цикла), вычисляется по программе расчета свойств хладагента на основании ввода давления Рк кДж/кг; S3’ – энтропия жидкого хладагента на линии насыщения в т. 3’ цикла, кДж/(кгК), вычисляется по программе расчета свойств хладагента на основании ввода давления Рк; i3 – энтальпия жидкого переохлажденного хладагента (т. 3 цикла), вычисляется по программе расчета свойств хладагента на основании ввода давления Рк и температуры на выходе из конденсатора Т3, кДж/кг; S3’ – энтропия жидкого переохлажденного хладагента в т. 3 цикла, кДж/(кгК), вычисляется по программе расчета свойств хладагента на основании ввода давления Рк и температуры Т3, кДж/(кгК); i4 = i3 – энтальпия жидкого хладагента (т. 4 цикла), принимается равной i3, кДж/кг; S4 – энтропия жидкого хладагента после дросселирования (т. 4 цикла), кДж/(кгК), вычисляется по программе расчета свойств хладагента на основании ввода давления Ро и энтальпии i3, кДж/(кгК); 1 – плотность пара хладагента на всасывании в компрессор (т. 1 цикла), кг/м3, вычисляется по программе расчета свойств хладагента на основании ввода давления Ро и температуры Т1; кр – коэффициент подачи компрессора, вычисляемый по уравнению в зависимости от Рн и Ро; Vд = = Vhкр. – действительная объемная производительность компрессора, м3/с. Последовательность эксергетического расчета 1. Работоспособность (максимальная полезная работа, эксергия) системы, состоящей из источника работы и окружающей среды, определяется выражением lmax = e = (i1 – i0) – Tокр.ср(s1 – s0). Здесь индексы 1 и 0 относятся соответственно к начальному (неравновесному) и конечному (равновесному) состояниям этой системы. 2. Техническое совершенство машины, работающей по обратному циклу, определяется ее КПД: = евых / евх , где евых – эксергия вырабатываемого холода или теплоты; евх – полная энергия, затраченная в установке. 3. Изменение эксергии холодильного агента в испарителе [6]: еи = е1 – е4 = i1 – i4 – Токр.ср (S1 – S4). Это изменение эксергии складывается из двух составляющих: часть эксергии (етн) отводится теплоносителем (полезный расход эксергии): еиполезн = qо(1 – Токр.ср/Ттн). Другая часть эксергии теряется из-за необратимости процесса теплообмена в испарителе: dи = (е1 – е4) – еиполезн . 4. Изменение эксергии потока в трубопроводе (всасывающем): dтр = е1’ – е1 . 5. Подводимая в компрессоре эксергия вычисляется по уравнению: екр = Nе / Gха . Внутренние и внешние потери эксергии в компрессоре составят: dкр екр – ℓкртеор. 6. Изменение эксергии холодильного агента в конденсаторе: dк = е2 – е3 . 7. Потери эксергии из-за необратимости при дросселировании: dдрос = е3 – е4 . 8. Экономические составляющие. Для оценки эффективности эксплуатации холодильных установок необходимо использовать величину стоимости выработанной или потерянной эксергии в виде функции [7]: Ср = Цт Т + Цz Z , где Ср – стоимость выработанного холода; Цт – цена топлива (электроэнергии), затраченного на выработку холода; Цz, Z – стоимость и количество дисконтированных затрат, банковского процента и т. п. Т. В. Морозюк отмечает [4], что эта методика неоднократно была апробирована как для стран со стабильной экономической ситуацией, так и для стран переходного нестабильного экономического периода и дала положительные результаты. При этом рекомендуется проводить расчеты с учетом региональной политики ценообразования. 9. Проверка работоспособности программы термоэкономического расчета. Выполнена для примера, приведенного в [8] – расчет одноступенчатого цикла холодильной машины. Результаты расчетов имеют сходимость 2–3 %. Аналогичная проверка (с аналогичными результатами) выполнена для одноступенчатой холодильной машины, работающей в режиме кондиционирования воздуха, эксергетический расчет которой приведен в [8]. Разработанная программа термоэкономического расчета аммиачной холодильной установки дала удовлетворительную сходимость результатов анализа конденсатора аммиачной абсорбционной холодильной установки [9]. Указанные алгоритмы являются основополагающими для разрабатываемой автоматизированной системы управления холодильной установкой. Программа, заложенная в микроконтроллер в режиме реального времени, отслеживает эксплуатационные параметры, регулирует холодильный процесс не только безопасно, но и, благодаря проведенным расчетам, экономически целесообразно. Выводы 1. Алгоритм идентификации эксплуатационных параметров основывается на исследовании аммиачной холодильной установки, т. к. данная тепловая машина характеризуется четко выраженной зависимостью эксплуатационных параметров от характеристик окружающей среды. 2. На основании литературных данных выбран метод и принцип идентификации эксплуатационных параметров и разработана блок-схема управления. 3. Подготовлена программа определения эффективности холодильной установки в длительном периоде эксплуатации. 4. Проведен эксергетический расчет аммиачной холодильной установки. 5. Итоговая разработанная система управления учитывает изменения наружных условий, что обеспечивает безопасность и экономичность эксплуатации, в отличие от существующих систем управления аммиачными холодильными установками.
Список литературы

1. Филоненко А. В. Идентификация предаварийных ситуаций установки получения серы методом Клауса: автореф. дис. … канд. техн. наук. – Астрахань, 2005. – 20 с.

2. Проталинский О. М., Савицкая Т. В. Применение методов искусственного интеллекта при автоматизации технологических процессов. – Астрахань: Изд-во АГТУ, 2004. – 184 с.

3. Егоров А. Ф. Управление безопасностью химических производств на основе информационных технологий. – М.: Колос, 2004. – 416 с.

4. Морозюк Т. В. Теория холодильных машин и тепловых насосов. – Одесса: Негоциант, 2006. – 712 с.

5. ПБ 09-595-03. Правила безопасности аммиачных холодильных установок.

6. Бродянский В. М. Эксергетический метод термодинамического анализа. – М.: Энергия, 1973. – 324 с.

7. Оносовский В. В. Моделирование и оптимизация холодильных установок. – Л.: Изд-во Ленинград. ун-та, 1990. – 208 с.

8. Пахомов М. М. Эксергетический анализ парокомпрессионной холодильной установки. – СПб.: Изд-во Санкт-Петербург. гос. ун-та, 2009. – 22 с.

9. Гуиди Т. К. Анализ эффективности промышленных аммиачных холодильных систем на основе термодинамического метода: дис. … канд. техн. наук. – Астрахань, 2010. – 148 с.