ВЕСТН. САМАР. ГОС. ТЕХН. УН-ТА. СЕР. ТЕХНИЧЕСКИЕ НАУКИ. 2018. № 2 (58)
УДК 620.193
МОДЕЛИРОВАНИЕ РЕГУЛИРУЮЩИХ КЛАПАНОВ В ПРОГРАММНОМ ПРОДУКТЕ «МиР ПиА»
С.Б. Коныгин
Самарский государственный технический университет Россия, 443100, г. Самара, ул. Молодогвардейская, 244
Аннотация. Рассмотрены вопросы интеграции регулирующих клапанов в модели технологических установок нефтегазовых производств, разрабатываемые в программном продукте «МиР ПиА». В качестве примера показано создание модели фрагмента установки, состоящего из трубопровода с регулирующим клапаном. Представлена настройка последовательности расчета технологической схемы. Приведены результаты гидравлического расчета для двух случаев. Показано, что в первом случае установка может работать в нормальном режиме, а во втором не справляется с заданной производительностью. Внедрение таких моделей позволит производственному персоналу оперативно решать ряд проблем, таких как оценка возможности работы установки при использовании другого сырья или другой производительности, выявление причин несоответствия продукции нормативным требованиям и т. д.
Ключевые слова: моделирование технологических процессов, регулирующий клапан, гидравлический расчет.
В настоящее время актуальным направлением является создание математических моделей действующих технологических установок нефтегазовой и химической промышленности [1]. Их внедрение позволит производственному персоналу оперативно решать ряд проблем, таких как оценка возможности работы установки при использовании другого сырья или другой производительности, выявление причин несоответствия продукции нормативным требованиям и т. д. Для этого возможно использование современных программных продуктов, к числу которых может быть отнесена программная платформа «МиР ПиА» [2, 3].
Одной из проблем при разработке такого рода моделей является интеграция в них элементов автоматизации технологического процесса [4, 5]. В рамках настоящей статьи продемонстрирована интеграция регулирующего клапана в математическую модель, созданную с помощью платформы «МиР ПиА» [2]. Для простоты рассмотрения выбран простой фрагмент установки, приведенный на рис. 1. Здесь поток жидкости движется по трубопроводу из емкости 81 с высоким давлением в емкость 82 с низким давлением, проходя при этом аппарат воздушного охлаждения Т1. Расход жидкости регулируется клапаном Б1.
Включение в модель установки элементов автоматизации стало возможным после существенной переработки модуля решения технологической схемы в целом. Центральным моментом данной модернизации стала возможность включения аппаратов в циклы с прямым (значения давления вычисляются от начала к концу цикла) и обратным (значения давления вычисляются от конца к началу цикла) направлениями гидравлического расчета.
Коныгин Сергей Борисович (д.т.н., доцент), заведующий кафедрой «Машины и оборудование нефтегазовых и химических производств». 166
Рис. 1. Моделируемый фрагмент технологической установки Созданная математическая модель данного участка представлена на рис. 2.
Файл Лид Объекты Настройка Еасчет Отладка Справка
ОС^Н ж -4 >+ Л ¿Ъ И I чЩ^ Щ ® ® В I ?
! Технологическая схема | Контроль параметров 1 Расчеты оборудования 1
-- о! ш 0 ш ®; : и т и шш □1 11 В0
Обозначение объекта на схеме П1 П2 пз П4 П5 П6 П7 П8 П9 П10 П11 П12
Общий массовый расход кг/ч 30000.000 30000.000 30000.000 30000.000 зоооо.ооо ЗОООО.ООО ЗОООО.ООО 30000.000 46114.895 46114.895
Общий мольный расход км оль/ч 1665.279 0.000 1665.279 1665.279 1665.279 1665.279 1665.279 1665.279 0.000 1665.279 1592.365 1592.365
Общий объемный расход мЗ/ч 31.075 31.075 31.075 31.078 31.078 30.678 30.678 30.678 39459.046 46890.777
Температура с 90.000 90.000 90.001 90.033 90.034 69.788 69.789 69.789 19.239 74.173
Давление кгс/см2 3.000 3.000 2.958 1.208 1.166 1.043 1.000 1.00С 1.00С 1.00С
Мольный расход пара кмопь/ч 0.000 0.000 0.000 0 000 0.000 0.000 0.000 0.000 0.000 0.000 1592.365 1592.365
Мольный расход жидкости кмопь/ч 1665.279 0.000 1665.279 1665.279 1665.279 1665.279 1665.279 1665.279 0.000 1665.279
Общая мольная энтальпия кДж/моль 6.795 6.795 6.795 6.795 6.795 5.264 5.264 5.264 0.559 2.160
И«™,,™.™, п.р. кДж/моль 0.559 2.160
кДж/кг 377.163 377.164 377.164 377.165 377.165 292.203 292.205 292.205 19.310 74.581
кДж/кг кДж/кг 19.310 74.581
Общая молекулярная масса 18.015 18.015 18.015 18.015 18.015 18.015 18.015 18.015 28.960 28.960
Молекулярная масса пара г/моль 28.960 28.960
Молекулярная масса жидкости г/моль 18.015 18.015 18.015 18.015 18.015 18.015 18.015 18.015
Общая плотность Плотность жидкости кг/МЗ 965.407 965.407 965.404 965.304 965.302 977.901 977.898 977.898 1.169 0.983
кг/мЗ 1.169 0.983
кг/МЗ 965.407 965.407 965.404 965.304 965.302 977.901 977.898 977.898
Мольная теплоемкость жидкости Массовая теплоемкость пара кДж/моль. К 0.029 0.029
кДжУмоль.К 0.076 0.076 0.076 0.076 0.076 0.075 0.075 0.075
кДж/кг.К 1.005 1.008
кассовая теплоемкость жидкости кДж/кг.К 4.205 4.205 4.205 4.205 4.205 4.188 4.188 4.188
Теплопроводность пара Вт/м.К 0.026 0.030
Теплопроводность жидкости Вт/м.К мПа.с мПа.с 0.673 0.314 0.673 0.673 0.673 0.673 0.660 0.660 0.660
Динамическая вязкость пара 0.01820 0.02078
Динамическая вязкость жидкости 0.314 0.314 0.314 0.314 0.405 0.405 0.405
Кинематическая вязкость жидкости мм2/С 0.325 0.325 0.325 0.325 0.325 0.414 0.414 0.414
Поверхностное натяжение жидкости мНУм 60.816 60.816 60.816 60.810 60.810 64.519 64.518 64.518
Общая тепловая мощность МВт 3.143 0.000 3.143 3.143 3.143 3.143 2.435 2.435 0.000 2.435 0.247 0.955
Схема решена Масштаб: 100%
Рис. 2. Модель фрагмента технологической установки, созданная в программном
продукте «МиР ПиА»
Настройка направления гидравлического расчета производится в диалоговом окне (рис. 3). Из рисунка видно, что емкость 81 и трубопровод Р1 рассчитываются в прямом, а остальные объекты - в обратном направлении. Причем важно, чтобы в цикле обратного гидравлического расчета первым стоял регулирующий элемент Б1, а последним - элемент с фиксированным давлением Б2. При такой организации процесса моделирования на каждой итерации замыкающий последовательность гидравлического расчета клапан пытается найти такое положение затвора, которое создает требуемый перепад давления (рис. 4).
В качестве примера на рис. 5 представлены результаты расчета фрагмента установки при расходе протекающей по нему жидкости 30 т/ч. Перепады давления на отметках 50 и 100 м представляют собой гидравлические сопротивления регулирующего клапана и аппарата воздушного охлаждения соответственно. По-
ложение рабочей точки клапана (см. рис. 4) показывает, что работа при заданной производительности возможна.
Последовательность и точность расчета
Циклы и аппараты Количество итераций Погрешность расхода Погрешность температуры Погрешность давления
И Начало цикла [>>) 50 0.0001 0.01 0.01
2 31 50 1 1 1
3 Р1 50 1 1 1
4 Начало цикла (<<) 50 0.0001 0.01 0.01
5 П 50 1 1 1
6 Р2 50 1 1 1
7 Т1 50 1 1 1
3 РЗ 50 1 1 1
9 зг 50 1 1 1
10 Конец цикла
11 Конец цикла
Примечание: погрешности расхода и температуры для циклов измеряются в долям расхода и Кельвинам соответственног а погрешности для аппаратов измеряются в долях от соответствующих погрешностей циклов.
ОК
Отмена
Рис. 3. Диалоговое окно для настройки последовательности и точности расчета
Результаты расчета регулирующего клапана И
Сводные данные
Характеристика клапана
щя
0.8
0.6'
0.4'
0.2-
— Расходная характеристика о Рабочая точка
и
0.2
0.4
0.6
о.э
Положение затвора клапана Н / Н1СС
Рис. 4. Диалоговое окно с результатами расчета регулирующего клапана
р, кгс/см2
Ь, м
Рис. 5. Профиль давления при расходе жидкости 30 т/ч
В качестве другого примера рассмотрим работу данной линии при расходе жидкости 100 т/ч. Полученные результаты моделирования приведены на рис. 6. Видно, что гидравлические сопротивления элементов схемы существенно увеличись, что привело к отрицательному расчетному перепаду давления на клапане. Из этого следует вывод о невозможности работы линии при данном расходе жидкости.
р, кгс/см2
Ь, м
Рис. 6. Профиль давления при расходе жидкости 100 т/ч
Все сказанное позволяет сделать вывод о возможности построения математических моделей действующих технологических установок, имеющих в своем составе элементы автоматизации, в программном продукте «МиР ПиА».
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Иваняков С.В., Крючков Д.А. Применение программного продукта «МиР ПиА» для компьютерного моделирования систем сепарации нефти // Вестник Самарского государственного технического университета. Сер. Технические науки. - 2018. - № 1 (57). - С. 168-172.
2. Коныгин С.Б., Крючков Д.А. Моделирование и расчет процессов и аппаратов (МиР ПиА). Свидетельство о государственной регистрации программы для ЭВМ №2015613176.
3. Ghasem N. Computer Methods in Chemical Engineering // Taylor & Francis Group, LLC. 2012. -492 p.
4. Маршалов Е.Д., Нечаева О.А. Имитационное моделирование гидравлических систем с регулирующими органами // Вестник ИГЭУ, 2007. - № 4. - С. 84-87.
5. Капранова А.Б., Лебедев А.Е., Мельцер А.М., Солопов С.А., Неклюдов С.В. О методах расчета гидравлического сопротивления регулирующих органов при транспортировании однокомпо-нентных сред // Фундаментальные исследования. - 2016. - № 4 (ч. 1). - С. 52-60.
Статья поступила в редакцию 2 февраля 2018 года
CONTROL VALVES MODELLING USING THE MIR PIA SOFTWARE S.B. Konygin
Samara State Technical University
244, Molodogvardeyskaya st., Samara, 443100, Russian Federation
Abstract. This paper is about the integration of the control valves in the gas and oil plant models that can be created by MiR PiA software. As example, a model of the simple part of the plant is shown. This model contains the pipeline with the control valve. A configuration of the calculation sequence is shown. A hydraulics calculation results are described for two cases. In the first case hydraulic line can work in the normal mode. In the second case the pressure drop is too big. An adoption of these models can help staff solve several problems, such as different raw material using, productivity increase, assurance of product quality problems, etc.
Keywords: modeling of technological processes, control valve, hydraulic calculation.
Sergey B. Konygin (Dr. Sci. (Techn.)), Professor.