DEVELOPMENT OF THE REFRIGERATION UNIT MANAGEMENT SYSTEM, CONSIDERING THE CHANGES OF ENVIRONMENT PARAMETERS
Journal: VESTNIK OF ASTRAKHAN STATE TECHNICAL UNIVERSITY. SERIES: MANAGEMENT, COMPUTER SCIENCE AND INFORMATICS
Volume 2013 № 1
, 2013
Rubrics: MANAGEMENT AND MODELING
Authors:
1.
Astrakhan State Technical University
Type:
Article
Pages:
from 59
to 64
Status:
Published
Received:
05.01.2013
Accepted:
25.01.2013
Published:
25.01.2013
Subject area:
VAC 05.12.2013 Системы, сети и устройства телекоммуникаций
VAC 05.13.01 Системный анализ, управление и обработка информации (по отраслям)
VAC 05.13.06 Автоматизация и управление технологическими процессами и производствами (по отраслям)
VAC 05.13.10 Управление в социальных и экономических системах
VAC 05.13.18 Математическое моделирование, численные методы и комплексы программ
VAC 05.13.19 Методы и системы защиты информации, информационная безопасность
UDK 621.565.2.003.13:[681.51.012:004.8]
GRNTI 20.01 Общие вопросы информатики
GRNTI 28.01 Общие вопросы кибернетики
GRNTI 49.01 Общие вопросы связи
GRNTI 50.01 Общие вопросы автоматики и вычислительной техники
GRNTI 82.01 Общие вопросы организации и управления
BBK 31.392.2-64:32.965,5
VAC 05.13.01 Системный анализ, управление и обработка информации (по отраслям)
VAC 05.13.06 Автоматизация и управление технологическими процессами и производствами (по отраслям)
VAC 05.13.10 Управление в социальных и экономических системах
VAC 05.13.18 Математическое моделирование, численные методы и комплексы программ
VAC 05.13.19 Методы и системы защиты информации, информационная безопасность
UDK 621.565.2.003.13:[681.51.012:004.8]
GRNTI 20.01 Общие вопросы информатики
GRNTI 28.01 Общие вопросы кибернетики
GRNTI 49.01 Общие вопросы связи
GRNTI 50.01 Общие вопросы автоматики и вычислительной техники
GRNTI 82.01 Общие вопросы организации и управления
BBK 31.392.2-64:32.965,5
Language:
Russian
Keywords:
model, analysis, refrigerating machine, identification of operational parameters
Abstract (English):
The developed methods of identification of operational parameters of technological installations allow development of the control systems ensuring their safety. However, these systems do not take into account the changes in external operational conditions. To process the algorithm of identification of operational parameters a heat machine – an ammoniac refrigerating machine is chosen, for this machine the dependence of operational parameters on the characteristics of the environment is clearly defined. On the basis of literary data the principle of identification of operational parameters is chosen and the management flowchart is developed. The model is chosen and the program of determination of efficiency of the refrigeration unit in the long period of operation is prepared.
The developed methods of identification of operational parameters of technological installations allow development of the control systems ensuring their safety. However, these systems do not take into account the changes in external operational conditions. To process the algorithm of identification of operational parameters a heat machine – an ammoniac refrigerating machine is chosen, for this machine the dependence of operational parameters on the characteristics of the environment is clearly defined. On the basis of literary data the principle of identification of operational parameters is chosen and the management flowchart is developed. The model is chosen and the program of determination of efficiency of the refrigeration unit in the long period of operation is prepared.
Keywords:
model, analysis, refrigerating machine, identification of operational parameters
model, analysis, refrigerating machine, identification of operational parameters
Введение Системному анализу технологических процессов и идентификации эксплуатационных и аварийных ситуаций с использованием новых компьютерных технологий и методов искусственного интеллекта посвящены многие работы в связи с актуальностью данного вопроса. Однако они не учитывают особенности технологического процесса, связанные с влиянием параметров окружающей среды. Так, А. В. Филоненко [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. Filonenko A. V. Identifikaciya predavariynyh situaciy ustanovki polucheniya sery metodom Klausa: avtoref. dis. … kand. tehn. nauk. - Astrahan', 2005. - 20 s.
2. Protalinskiy O. M., Savickaya T. V. Primenenie metodov iskusstvennogo intellekta pri avtomatizacii tehnologicheskih processov. - Astrahan': Izd-vo AGTU, 2004. - 184 s.
3. Egorov A. F. Upravlenie bezopasnost'yu himicheskih proizvodstv na osnove informacionnyh tehnologiy. - M.: Kolos, 2004. - 416 s.
4. Morozyuk T. V. Teoriya holodil'nyh mashin i teplovyh nasosov. - Odessa: Negociant, 2006. - 712 s.
5. PB 09-595-03. Pravila bezopasnosti ammiachnyh holodil'nyh ustanovok.
6. Brodyanskiy V. M. Eksergeticheskiy metod termodinamicheskogo analiza. - M.: Energiya, 1973. - 324 s.
7. Onosovskiy V. V. Modelirovanie i optimizaciya holodil'nyh ustanovok. - L.: Izd-vo Leningrad. un-ta, 1990. - 208 s.
8. Pahomov M. M. Eksergeticheskiy analiz parokompressionnoy holodil'noy ustanovki. - SPb.: Izd-vo Sankt-Peterburg. gos. un-ta, 2009. - 22 s.
9. Guidi T. K. Analiz effektivnosti promyshlennyh ammiachnyh holodil'nyh sistem na osnove termodinamicheskogo metoda: dis. … kand. tehn. nauk. - Astrahan', 2010. - 148 s.
