Моделирование и оптимизация установки выделения товарного пропилена

или тройной комплекс ТФМ-красителя взаимодействуют между собой, образуя слабые комплексы, которые по величине константы устойчивости можно охарактеризовать как комплексы с перено-

сом заряда. Данные комплексы, вероятно, могут быть центрами фотогенерации носителей заряда при спектральной сенсибилизации поли^-эпо-ксипропилкарбазола ТФМ-красителями.

СПИСОК ЛИТЕРАТУРЫ

1. Давиденко Н.А., Деревянко Н.А., Ищенко А.А., Кувшин-ский Н.Г., Павлов В.А., Чуприна Н.Г. Ионные и мероцианино-вые красители в качестве сенсибилизаторов. Голографические регистрирующие среды на основе пленок поли^-эпоксипро-пилкарбазола // Журн. науч. и прикл. фотографии. - 2002. -Т. 47. - № 6. - С. 29-37.

2. Шафферт Р. Электрофотография. - М.: Мир, 1968. - 448 с.

3. Ермолаев В.Л., Бодунов В.Н., Свешников Е.Б., Шахвердов Т.А. Безызлучательный перенос энергии электронного возбуждения. - Л.: Наука, 1977. - 310 с.

4. Демидов К.Б., Акимов И.А., Газиев З.А. Исследование сенсибилизированной фотографической чувствительности в поли-N-винилкарбазоле с помощью Штарк-спектроскопии // Журн. науч. и прикл. фото- и кинематографии. - 1985. - Т. 30.

- № 1. - С. 56-59.

5. Williams D.J. Photoconductivity in polymers. An Inderdisciplinary approach / Ed. by A.V. Patsis, D.A. Seanor. - USA: Technomic, 1976. - 274 p.

6. Орлов И.Г., Кошелев К.К., Матасова Г.И., Яковлев В.Б. Физико-химические свойства полимерных светочувствительных систем // Журн. науч. и прикл. фото- и кинематографии. - 1982.

- Т. 27. - № 1. - С. 34-37.

7. Румянцев Б.М., Балобанов Е.И., Букин Ю.И., Семенова Л.В., Мельничук Л.А., Юдина Г.И. О механизме сенсибилизации красителем внутреннего фотоэффекта в полимерных полупроводниках // X координационное совещание по органическим полупроводникам. - Косов, 1979. - С. 66-67.

8. Коботаева Н.С., Сироткина Е.Е., Микубаева Е.В. Спектральная сенсибилизация фотопроводимости поли-^эпоксипро-пилкарбазола и дифенилгидразонов бензальдегида тройными комплексами трифенилметановых красителей // Химия высоких энергий. - 2005. - Т. 39. - № 5. - С. 362-366.

9. Пилипенко А.Р., Тананайко М.Ш. Разнолигандные и разноме-тальные комплексы и их применение в аналитической химии. - М.: Химия, 1973. - 305 с.

10. Блюм И.А. Экстракционно-фотометрические методы анализа с применением основных красителей. - М.: Наука, 1970. - 219 с.

11. А.с. 503200 СССР. МКИ G03G 5/06. Электрофотографический материал / А.И. Ундзенас, В.И. Гайдялис, И.Б. Сидаравичюс, Р.И. Каволюнас, И.И. Зданавичус, Н.К. Дуобинис. Заявлено 08.02.1972; Опубл. 1976, Бюл. № 6. - 110 с.

12. Акимов И.А. Функции локальных центров фотопроводника в акте спектральной сенсибилизации // Журн. науч. и прикл. фото- и кинематографии. - 1980. - Т. 25. - № 3. - С. 228-250.

13. Гурьянова Е.Н., Гольдштейн И.П., Ромм И.П. Донорно-акцеп-торная связь. - М.: Химия, 1973. - 400 с.

УДК 541.122.4:66.048

МОДЕЛИРОВАНИЕ И ОПТИМИЗАЦИЯ УСТАНОВКИ ВЫДЕЛЕНИЯ ТОВАРНОГО ПРОПИЛЕНА

А.А. Ляпков, Ю.В. Шефер

Томский политехнический университет E-mail: alex@chtd.tpu.ru

Проведено моделирование процесса непрерывной многокомпонентной ректификации в тарельчатой колонне и оптимизация параметров ее работы, обеспечивающая выделение товарного пропилена требуемого качества. Проанализированы основные факторы, влияющие на качество товарного пропилена. Сделан вывод о том, что для достижения нужного качества пропилена необходимо либо снижать расход по питанию, либо повышать давление в колонне, не меняя конструкцию и конфигурацию оборудования. Заменив существующие тарелки на более эффективные (например, решетчатого типа) и изменив обвязку колонн таким образом, чтобы сделать из них две самостоятельные, можно существенно улучшить работу установки.

Разделение жидких и газообразных смесей на индивидуальные компоненты имеет большое значение в химической и нефтехимической промышленности. Особенно важно выделение компонентов в чистом виде в производстве полимеров с повышенными физико-химическими свойствами. Для разделения смесей на индивидуальные компоненты среди других методов широкое применение нашел метод ректификации [1-6]. Достижения в развитии термодинамики необратимых процессов дали возможность выразить и количественно оценить довольно сложную по своей природе движу-

щую силу процесса массопередачи в условиях ректификации [7-12].

В течение последних 30 лет в сырьевой базе отечественной и мировой нефтехимии ведущая роль принадлежит низшим олефинам - этилену и пропилену. Основным источником их производства служит процесс термического пиролиза углеводородов с водяным паром. Именно на установках пиролиза получают сегодня первичные продукты, обеспечивающие сырьем производство пластических масс, синтетических смол, каучуков и волокон.

Процесс пиролиза непрерывно совершенствуется, и в настоящее время эксплуатируются высокоэффективные этилен-пропиленовые установки производительностью 300...500 тыс. т в год и более. При выборе сырья пиролиза в нашей стране был взят курс преимущественно на использование пря-могонного бензина. На современной отечественной крупнотоннажной этилен-пропиленовой установке типа ЭП-300, расположенной на территории северного промузла г. Томска, мощностью по этилену 300 тыс. т в год, работающей на прямогонном бензине, получают одновременно 139 тыс. т в год пропилена. Продукты, получаемые на установке ЭП-300, являются сырьем для крупных производств пластмасс - заводов «Полиэтилен» и «Полипропилен».

В последние годы в связи с обостряющимся в стране дефицитом бензина все более актуальным направлением становится использование широкой фракции легких углеводородов на установках пиролиза высокой мощности, запроектированных на бензиновое сырье. Наиболее перспективным на данный момент считается пиролиз прямогонного бензина с добавкой широкой фракции легких углеводородов. Перспективность использования такого вида сырья объясняется:

• низкой себестоимостью по сравнению с прямо-

гонным бензином;

• большей селективностью по целевым продуктам;

• имеющейся сырьевой базой в пределах области.

При использовании широкой фракции легких углеводородов в качестве добавки к сырью - бензину в количестве 20...25 об. % снижается себестоимость пропилена.

По проекту в составе сырья, идущего на газоразделение в пропиленовую колонну, содержание пропилена должно быть не менее 92 % об., а по результатам анализа работы установки за 2004 г. содержание пропилена колебалось в пределах 62...87 об. %. Снижение содержания пропилена привело к увеличению расхода питания пропиле-новой колонны с 21,8 до 23,1 т/ч.

В связи с изменением состава сырья, назрела необходимость поиска оптимальных параметров проведении процесса ректификации и возможной реконструкции оборудования. Настроенное на разделение продуктов пиролиза прямогонного бензина, оно перестало обеспечивать необходимое качество выпускаемой продукции. На старой установке невозможно получение товарного пропилена с содержанием последнего более 99,8 об. %.

Для обеспечения выпуска пропилена заданного качества при сохранении или даже некотором увеличении объема его выпуска с помощью созданной математической модели были проанализированы основные факторы, влияющие на работу установки выделения товарного пропилена.

Рис. 1. Схема колонны и потоков жидкости и пара в ней: 1) тарелка питания; 2) контактное устройство; 3) дефлегматор; 4) кипятильник. Расходы: F, FV, FL - питания полный, в паровой и жидкой фазах, соответственно; V,, Vе, V" - пара в ¡-ом сечении колонны, в верхней и нижней частях колонны, соответственно; Ц, и,Ц" - жидкости в ¡-ом сечении колонны, в верхней и нижней частях колонны, соответственно; D -дистиллята; - кубового остатка

В качестве основных параметров, по которым проводилась оптимизация работы установки, были выбраны:

• расход и состав питания;

• давление и температура в колонне;

• расход флегмы;

• тип контактных устройств.

Оптимизация производилась с целью выявления параметров, наиболее влияющих на процесс ректификации и на качество получаемого пропилена. Также в ходе выполнения работы были рассмотрены следующие вопросы:

• определение максимального расхода и минимального состава питания, при котором возможно ещё получение товарного пропилена высшего сорта на данной установке без замены оборудования;

• подбор типа и количества тарелок для возможной их замены, чтобы справиться с новым сырьём и ещё большей нагрузкой на будущее. На основании теоретических [1-6, 13-20] и

экспериментальных данных была построена мо-

дель колонны выделения товарного пропилена, которая состоит из модели тарелки питания, модели контактных устройств, кипятильника и дефлегматора (см. схему на рис. 1).

При построении модели были сделаны следующие допущения [21-25]:

• давление в пределах данного контактного устройства постоянно;

• жидкость в колонне находится при температуре кипения, а пар - при температуре точки росы;

• разделяемые смеси близки к идеальным;

• физико-химические свойства компонентов усреднены в диапазоне изменения концентраций на текущем контактном устройстве;

• потоки жидкости и пара постоянны;

• для описания гидродинамики жидкой фазы принята ячеистая модель, а для паровой фазы -модель идеального вытеснения;

• равновесная кривая линейна в пределах изменения концентрации на контактном устройстве;

• локальный КПД постоянен в пределах контактного устройства;

• массопередача по фазам независимая, диффузионные сопротивления аддитивны, т. е. общий коэффициент массопередачи может быть разложен на частные;

• коэффициенты массоотдачи постоянны в пределах контактного устройства;

• унос жидкости постоянен;

• теплоты смешения потоков жидкости пренебрежимо малы;

• режим работы контактных устройств - адиабатический.

Математическая модель учитывает такие факторы как тепло- и массоперенос в виде дифференциальных уравнений и гидродинамику в виде критериальных уравнений.

Уравнения, описывающие контактное устройство, аналогичны приведенным в [23].

Величина уноса зависит от конструктивных особенностей контактного устройства и физико-химических свойств компонентов.

Расчет основных физических параметров для пара и жидкости вели по общепринятым уравнениям, построенным на основе уравнений состояния идеального газа и критериальных параметров.

На рис. 2 представлена блок-схема программы оптимизации параметров работы колонны многокомпонентной ректификации.

Программа функционально состоит из следующих расчетных блоков:

• расчета физических свойств веществ (вязкость, плотность и др.), описанных в форме полиномиальных аппроксимаций;

• описания математической модели тарелки питания;

Рис. 2. Блок-схема программы оптимизации установки выделения товарного пропилена

• расчета контактного устройства с подбором температуры, при которой жидкость данного соста-

ва кипит на данном контактном устройстве, эффективности контактного устройства, расчета теплот фазовых переходов при испарении жидкости и теплоёмкости жидкости, стекающей с предыдущего контактного устройства; •подбора расхода пара и флегмы; •вывода данных.

В созданной программе рассчитывается реальное количество контактных устройств, необходимых для получения конечного продукта с содержанием пропилена 99,8 об. %.

приблизительно равно уменьшению питания на 0,9 т/ч (23,1 до 22,2 т/ч). При этом нужно учитывать экономические и производственные факторы. Колонна должна перерабатывать всё сырье, поступающее на разделение, а это ведёт к повышению расхода греющего пара и соответственно росту переменных издержек производства, увеличению себестоимости продукции. Уменьшение подачи бензина в печи пиролиза снизит нагрузку по питанию на всех стадиях производства, понизит давление в ректификационных колоннах и уменьшит расход греющего пара в теплообменниках.

Рис. 3. Зависимость количества тарелок от давления и расхода флегмы

Анализ полученных результатов позволил выявить следующие закономерности. Увеличение давления с 1,4 до 2,6 МПа, рис. 3, ведёт к уменьшению количества тарелок, необходимых для достижения заданной степени разделения. В свою очередь, увеличение расхода флегмы с 0,4 до 0,5 т/ч ведёт лишь к минимальному изменению состава дистиллята. Скорее всего, это связано с большой разницей по расходу пара и жидкости в колонне. Если рассмотреть состояние системы на верхней тарелке, то можно заметить, что через её отверстия проходит приблизительно 14,3 т/ч пара и только 0,5 т/ч жидкости на орошение колонны. Своим малым объёмом флегма не может значительно обогатить жидкость на тарелке. Вследствие этого невозможно добиться улучшения содержания пропилена в дистилляте путём увеличения расхода флегмы; ее нужно подавать ровно столько, сколько необходимо для проведения процесса ректификации на самой верхней тарелке.

Изменяя давление и расход питания, рис. 4, в качестве отклика получаем, что на эффективность работы тарелок сильное влияние оказывает увеличение давления. Однако ещё большее влияние на КПД тарелок оказывает снижение расхода питания. Так, увеличение давления с 2,1 до 2,2 МПа

Рис. 4. Зависимость эффективности контактного устройства от давления и расхода питания колонны

Уменьшение подачи прямогонного бензина в печь пиролиза приведёт к уменьшению себестоимости конечных продуктов на всех стадиях производства и к улучшению их качества. Применение подобных моделей для отыскания оптимальных режимов работы всех колонн ЭП-300 позволит удерживать качество целевых продуктов на стабильно высоком уровне.

Обеднение сырья диктует необходимость увеличения давления (рис. 5) таким образом, чтобы количество тарелок для достижения требуемой чистоты товарного пропилена осталось неизменным. Поскольку увеличение давления ведёт к удорожанию товарного пропилена, то с уменьшением содержания в сырье пропилена остаётся только увеличить количество тарелок для получения товарного пропилена нужного качества.

Учитывая, что сырьё, поступающее на разделение, имеет колебания по содержанию пропилена от партии к партии поступающего на пиролиз прямо-гонного бензина, необходимо обеспечить качественно стабильный выход продукта даже при самом низком его содержании в сырье. Поэтому было интересно посмотреть, как на этих параметрах отразится изменение конструкции тарелок.

Рис. 5. Зависимость количества тарелок от содержания пропилена в сырье и давления

Экономические расчеты пяти различных типов тарелок показали, что наиболее эффективными являются решетчатые тарелки. Для достижения требуемой степени разделения их нужно почти в два раза меньше тех, что применяются сейчас.

При реконструкции колонны выделения товарного пропилена планировалось заменить 201 тарелку ситчатого типа на 250 решетчатых провальных тарелок типа ECMD. Проектная фирма гарантировала стабильное получение пропилена с содержанием в дистилляте - 99,8 об. % при его содержании в сырье минимум до 61,1 об. %.

Согласно нашим расчетам при замене имеющихся тарелок на решетчатые, их будет необходимо всего 123 для выделения пропилена высшего сорта. Вследствие этого, можно рекомендовать вместо одной разрезной колонны, изменив обвязку, сделать две, снизив при этом расход по питанию, например, до 18,0 т/ч в каждой из колонн. При такой реконструкции две новые колонны будут способны переработать на 56 % больше сырья состава 61,1 об. % и смогут разделять ещё более обедненное сырьё (таблица).

Сравнительные расчеты колонны, оборудованной пятью различными типами тарелок, связывающие их условную стоимость и количество необходимое для достижения нужной степени разделения

СПИСОК ЛИТЕРАТУРЫ

1. Александров И. А. Ректификационные и абсорбционные аппараты. Методы расчета и основы конструирования. - М.: Химия, 1971. - 296 с.

2. Баготуров С.А. Теория и расчет перегонки и ректификации. -М.: Гостоптехиздат, 1961. - 436 с.

показали, что решетчатые тарелки почти в два раза более эффективны, чем применяемые сейчас.

Таблица. Сравнение параметров работы колонн

Показатель Рабочая колонна Колонна по плану реконструкции Две колонны

Тарелок выше питающей, шт. 131 162 81+81

Тарелок ниже питающей, шт. 70 88 44+44

Расход питания, т/ч 23,1 23,1 18,0+18,0

Расход дистиллята,т/ч 13,8 13,8 10,0+10,0

Содержание пропилена в дистилляте, об. % 99,70 99,9976 99,8964

Выводы

Наибольшее влияние на качество получаемого пропилена оказывает расход питания, причем его снижение с 23,1 до 22,2 т/ч оказывает на эффективность контактного устройства такое же влияние, как и увеличение давления с 2,1 до 2,2 МПа.

При увеличении давления с 1,4 до 2,6 МПа происходит увеличение содержания пропилена в дистилляте, но это ведет к повышенному расходу греющего пара.

При обеднении пропаном сырья, идущего на разделение, необходимо большее количество тарелок при прочих равных условиях для получения содержание пропилена в дистилляте не менее 99,8 об. %.

Для достижения нужного качества пропилена можно изменять два параметра: если возможно -снизить расход по питанию до 22...21 т/ч, если нет - то повышать давление в колонне до 2,2...2,3 МПа, не меняя конструкцию и конфигурацию оборудования.

Замена существующих тарелок на более эффективные, например на решетчатые, позволит получать продукт существенно лучшего качества.

Применив более эффективный тип тарелки и изменив обвязку колонн таким образом, чтобы сделать из них две самостоятельные, можно существенно улучшить работу установки. Например, снизив расход по питанию до 17...18 т/ч, можно получать пропилен высшего сорта из более обедненных продуктов пиролиза, а при сохранении содержания пропилена в сырье на существующем уровне общий расход питания может быть увеличен по сравнению с действующим в 1,5 раза.

3. Chan Willie K., Boston Joseph F., Evans Lawrence B. Select the right software for modeling separation processes // Chem. Eng. Progr. - 1991. - Т. 87. - № 9. - Р. 63-69.

4. Щеголев В. В. Классификация колонных аппаратов. Интенсификация процессов тепло- и массообмена при проведении процессов ректификации // Химическая промышленность. -2005. - № 8. - С. 389-396.

5. Iglesias M., Tojo J., Dominguez A. Rectification de sistemas multi-componentes. 1. Simulacion de una columna empleando contribucion de grupos // Afinidad. - 1995. - V. 52. - № 455. - Р. 12-18.

6. Леонтьев В.С. Компьютерное моделирование процессов ректификации // Химическая промышленность. - 2005. - № 7. -С. 334-346.

7. Карапетьянц М.Х. Химическая термодинамика. - М.: Химия, 1975. - 582 с.

8. Рид Р., Праусниц Дж., Шервуд Т. Свойства газов и жидкостей: Справочное пособие. - Л.: Химия, 1982. - 592 с.

9. Рабинович Г.Г., Рябых П.М., Хохряков П.А. и др. Расчеты основных процессов и аппаратов нефтепереработки: Справочник / Под ред. Е.Н. Судакова. - М.: Химия, 1979. - 568 с.

10. Холланд Ч.Д. Многокомпонентная ректификация. - М.: Химия, 1969. - 348 с.

11. Хиао-Мин Ху, Фу Шен. Применение оценки скорости переноса для расчета бинарных коэффициентов массопередачи и моделирование многокомпонентной ректификации // Теоретические основы химической технологии. - 1992. - Т. 26. - № 4. - С. 486-493.

12. Платонов В.М., Берго Б.Г. Разделение многокомпонентных смесей. Расчет и исследование ректификации на вычислительных машинах. - М.: Химия, 1965. - 368 с.

13. Taylor R., Achuthan K., Lucia A. Complex-domain distillation calculations // Comput. and Chem. Eng. - 1998. - V. 22. - № 12. -Р. 1731-1732.

14. Петлюк Ф.Б., Серафимов Л.А. Многокомпонентная ректификация. Теория и расчет. - М.: Химия, 1983. - 304 с.

15. Diwekar U.M., Malik R.K., Madhavan K.P. Optimal reflux rate policy determination for multicomponent batch distillation columns // Comput. and Chem. Eng. - 1987. - V. 11. - № 6. - Р. 629-637.

16. Ульянов Б.А., Щелкунов Б.И. Процессы и аппараты химической технологии: Гидравлика контактных тарелок. - Иркутск: Изд-во Иркут. ун-та, 1996. - 160 с.

17. Павлечко В.Н. Комплексная модель эффективности ректификационных тарелок. 5. Перекрестное движение фаз при перемешивании жидкости // Инженерно-физический журнал. -2001. - Т. 74. - № 3. - С. 177-180.

18. Павлечко В.Н. Комплексная модель эффективности ректификационных тарелок. 7. Взаимосвязь отдельных параметров // Инженерно-физический журнал. - 2002. - T. 75. - № 1. - С. 112-116.

19. Pamfil V. Modele matematice ale corelatiei Gilliland // Rev. Chim. - 1983. - V. 34. - № 11. - Р. 1022-1024.

20. Rooney J.M. Simulating batch distillation // Chem. Eng. - 1984. -Т. 91. - № 10. - P. 61-64.

21. Alopaeus V., Aittamaa J. Appropriate simplifications in calculation of mass transfer in a multicomponent rate-based distillation tray model // Ind. and Eng. Chem. Res. - 2000. - V. 39. - № 11. -Р. 4336-4345.

22. Sharif M., Shah N., Pantelides C.C. On the design of multicom-ponent batch distillation columns // Comput. and Chem. Eng. -1998. - V. 22. - Appl. - Р. 69-76.

23. Анисимов И.В., Бодров В.И., Покровский В.Б. Математическое моделирование и оптимизация ректификационных установок. - М.: Химия, 1975. - 214 с.

24. Дьяконов С.Г., Лаптев А.Г., Данилов В.А. Определение объемных коэффициентов массоотдачи с помощью математической модели при расчете тарелок с прямоточными клапанами // Химическая промышленность. - 1991. - № 8. - С. 499-501.

25. Jimenez L., Basualdo M. S., Gomez J.C., Toselli L., Rosa M. Nonlinear dynamic modeling of multicomponent batch distillations: a case study // Braz. J. Chem. Eng. - 2002. - V. 19. - № 3. - Р. 307-317.

УДК 536.46

МОДЕЛИРОВАНИЕ АКУСТИЧЕСКОЙ НЕУСТОЙЧИВОСТИ В КАМЕРЕ РАКЕТНОГО ДВИГАТЕЛЯ

ТВЕРДОГО ТОПЛИВА

Р. Мырзакулов*, М.Ж. Козыбаков**, К.О. Сабденов

*Евразийский национальный университет, г. Астана **Шымкентский социально-педагогический университет, г. Шымкент Томский политехнический университет E-mail: sabdenovko@mail.ru

Изучается возникновение, развитие и последствия акустической неустойчивости в камере ракетного двигателя с твердотопливным зарядом. Рассматриваются относительно низкочастотные колебания с периодом, намного превышающим период собственных колебаний камеры. Но частота изменения термодинамических параметров находится в широких пределах и может быть сравнима с собственной частотой зоны горения. Неустойчивость может приводить к автоколебательному горению, или к хаотическому режиму, или к погасанию горения.

Введение

Еще в конце 30-х гг. прошлого века разработчики системы залпового огня «Катюша» столкнулись со странным явлением. Как оказалось, ряд оригинальных технических решений, который позволил бы достичь больших скоростей реактивных снарядов, не осуществим из-за особенностей горения топлива. Оно, превосходно сгорая на открытом воз-

духе, при определенных условиях «не желало» гореть в камере, или же горело настолько нестабильно, что совершенно приводило в негодность двигатель. Как оказалось, наличие окислителя и горючего в составе топлива не гарантирует равномерное протекание химической реакции, необходимо еще, чтобы внутренние закономерности механизма горения оптимально сочетались с конструкционными особенностями самого двигателя ракеты. Даль-