автореферат диссертации по химической технологии, 05.17.08, диссертация на тему:Моделирование статики и динамики жидкофазной адсорбции на цеолитах

На правах рукописи

Заиченко Надежда Викторовна

□□3171572

МОДЕЛИРОВАНИЕ СТАТИКИ И ДИНАМИКИ ЖИДКОФАЗНОЙ АДСОРБЦИИ НА ЦЕОЛИТАХ

Специальность 05.17.08 Процессы и аппараты химических технологий

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

О 5 [{№ 2003

Уфа - 2008

003171572

Работа выполнена в Уфимском государственном нефтяном техническом университете

Научный руководитель доктор технических наук профессор

Самойлов Наум Александрович

Официальные оппоненты доктор технических наук профессор

Умергалин Талгат Галеевич, кандидат технических наук старший научный сотрудник Балаев Александр Всеволодович

Ведущая организация государственное унитарное предприятие

«Институт нефтехимпереработки» Республики Башкортостан

Защита состоится «25» июня 2008 года в Ю'ОО часов на заседании совета по защите докторских и кандидатских диссертаций Д 212 289 03 при Уфимском государ ственном нефтяном техническом университете по адресу ул Космонавтов, корп 1, г Уфа, Республика Башкортостан, 450062

С диссертацией можно ознакомиться в библиотеке Уфимского государственного нефтяного технического университета

Автореферат разослан «24» мая 2008 года

Ученый секретарь диссертационного совета

КГ Абдульминев

Общая характеристика работы

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

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

Диссертационная работа выполнена в соответствии с координационным пла ном Научного совета по адсорбции и хроматографии РАН на 2005-2008 гг

Цель работы Основными целями диссертации являются'

- поиск достаточно простых способов получения кинетических и статических параметров для разработки математической модели динамики жидкофазной адсорбции на цеолитах,

- оптимизация работы адсорбционной колонны

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

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

- оценить характер лимитирующего диффузионного сопротивления при жидкофазной адсорбции на цеолитах,

- получить кинетические параметры, необходимые для расчета динамики жидкофазной адсорбции,

- разработать метод расчета динамики адсорбции органических примесей цеолитами на основе минимального объема экспериментальных данных о кинетике адсорбции, доказать адекватность расчетов по разработанной модели данным контрольных экспериментов по динамике жидкофазной адсорбции, выполнить оценку параметрической чувствительности математической модели,

- разработать метод оптимизации адсорбционной установки, учитывающий затраты на адсорбционный процесс затраты на адсорбент, перекачку, металл конструкций и регенерацию сорбента

Научная новизна Показано, что общий объемный коэффициент массопере-дачи, рассчитанный как величина, обратная первому статистическому моменту кинетической кривой статики адсорбции, является параметром, независимым от концентрации сорбируемой примеси ароматического углеводорода в бинарной смеси с парафиновым растворителем

Общий объемный коэффициент массопередачи использован при разработке математической модели динамики жидкофазной сорбции ароматических углеводородов в гептане

Выполнена оценка параметрической чувствительности математической модели динамики жидкофазной адсорбции бинарной смеси углеводородов

Разработан метод оптимизации работы адсорбционной установки с учетом затрат на основное оборудование, адсорбент, стоимость перекачки сырья и регенерацию адсорбента.

Практическая пенност ь Предложен способ определения кинетических коэффициентов на основе минимального объема экспериментальных данных и способ расчета коэффициента массопередачи при жидкофазной адсорбции ароматических углеводородов в растворе н-гептана на цеолитах МаХ с разработкой методических материалов, которые переданы в ГУЛ «Институт нефтехимпереработки РБ» для использования в научных исследованиях, проводимых институтом Предложена модель оптимизации жидкофазной адсорбции углеводородов на цеолитах, учитывающая вклад каждого вида затрат (на основную-аппаратуру—адеор- -

бент, перекачку очищаемого сырья и топливо) в общие затраты на реализацию адсорбционного процесса

Апробация работы Основные результаты диссертационной работы докладывались на VI конгрессе нефтегазопромышленников "Нефтегазопереработка и нефтехимия" (Уфа, 2005), III Всероссийской научно-практической конференции" Нефтегазовые и химические технологии " (Самара, 2005), Всероссийском семинаре "Термодинамика поверхностных явлений и адсорбции" (Иваново, 2005), конференции "Перспективы развития химической переработки горючих ископаемых" (Санкт-Петербург, 2006), Международной научно-практической конференции "Нефтегазопереработка и нефтехимия" (Уфа, 2006), III Всероссийской научной конференции "Теория и практика массообменных процессов химической технологии (Марушкинские чтения)" (Уфа, 2006), Всероссийской научной конференции "Переработка углеводородного сырья. Комплексные решения. (Левиптерские чтения)" (Самара, 2006), XI Всероссийском симпозиуме с участием иностранных ученых "Актуальные проблемы теории адсорбции, пористости и адсорбционной селективности" (Москва-Клязьма, 2007), 56-й, 57-й и 58-й научно-технических конференциях студентов, аспирантов и молодых ученых (Уфа, 2005, 2006 ,2008)

Публикации По материалам диссертации опубликовано 2 статьи в журналах, 15 тезисов докладов

Объем и структура диссертации Работа изложена на 199 страницах, состоит из введения, четырех глав, выводов и трех приложений Диссертация содержит 34 рисунка, 32 таблицы и список использованных источников из 151 наименования

Основное содержание работы

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

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

Во второй главе представлены результаты исследования кинетики сорбции ароматических углеводородов (бензола, толуола, п-ксилола, изопропилбензола и нафталина) из раствора в н-гептане цеолитом №Х при температуре 20°С в диапазоне концентраций 1-80 % об В главе приводится методика проведения экспериментов, характеристика используемых реагентов и исследуемого цеолита, параметры активации адсорбента

По значениям адсорбционной способности цеолита, полученным экспериментальным путем, строились кинетические кривые с определением величины равновесной адсорбции для каждой концентрации при времени контакта т—>к> (т = 24 часа), по значениям которых строились изотермы избыточной (гиббсовской) адсорбции ароматических углеводородов цеолитом ЫаХ из расгвора в н-гептане (рисунок 1)

0,3

^ 0,25

£ 0 2 а '

0,15 0,1 0,05 0

0 ♦ бензол ■ толуол а п-ксилол 100 О ИНЬ -О нафталин

С0, % об

Рисунок 1 - Зависимость избыточной (гиббсовской) адсорбции (аг) от концентрации ароматических углеводородов (Со) в растворе н-гептана

Из рисунка 1 видно, что в области низких концентраций ароматического углеводорода величина избыточной адсорбции увеличивается, достигая максимума при С0 = 9-18% об Кроме того, замечено, что активность цеолита №Х по бензолу несколько ниже, чем по толуолу, п-ксилолу, изопропилбензолу и нафталину, позиция максимума у гомологов бензола правее и выше, чем у бензола Эти факты можно объяснить усилением адсорбционных сил гомологов бензола за счет специфического взаимодействия я-связеи с ионами кристаллической структуры цеолита из-за наличия у них дипольного момента

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

Так, для выявления вклада миграции молекул адсорбата в общее сопротивление диффузионному массопереносу кинетика адсорбция рассматривалась в координатах [г] - [т1/2], где у = а /~ величина относительной адсорбции (рисунок 2)

1 1

у* / / • / л 1

/ / \ / /■ / /А \/ а / / / У* / / X—/—. А,/ ж < -1

О 5 10 15 20 т0'5

Концентрация толуола, % об • ♦ 0,8 * 5,4 • 21,3 " 48,8 * 78,8

Рисунок 2 - Зависимость относительной адсорбции толуола из раствора в н-гептане (у) от времени (т0'5) адсорбции (шаг по оси абсцисс для каждой концентрации толуола в пределах 0,8-78,8 % об. равен 2°5 мин)

Как следует из рисунка 2, относительная адсорбция толуола из раствора в н-гептане (/) в диапазоне у= 0 - (0,7 - 0,8) пропорциональна т0,3 независимо от на-

чальной концентрации толуола в растворе Пропорциональное изменение у с ростом т0'5 свидетельствует о существенной роли миграции молекул адсорбата в первичной пористой структуре цеолита в общем сопротивлении диффузионному мас-сопереносу при адсорбции из жидкой фазы Аналогичные данные получены и для других исследуемых ароматических углеводородов, Из рисунка 2 видно, что с увеличением концентрации ароматического углеводорода (независимо от его строения) в смеси возрастает крутизна кинетических кривых, а это, в свою очередь, свидетельствует об активации переноса

Известно, что при сорбции углеводородов цеолитом ИаХ из жидкой фазы для исследования возможности лимитирования сорбции внутридиффузионным сопротивлением можно воспользоваться уравнением

У2 =3<р-3<рг (1)

где <р = г/г„ - относительное время

Для обработки экспериментальных данных по уравнению (1) целесообразно пользоваться его линейной формой

тогда при справедливости уравнения (1) экспериментальные точки в координатах [(1-УГ ] - [т ]Д0ЛЖНЫ располагаться на прямой линии (рисунок 3)

О

20

40

60

Концентрация толуола. % об

■0,8% об ■21,3% об

• 3,8% об ■ 36,8% об

■5,4% об -48,8% об

X, мин

■ 9,3% об • 78,8% об

Рисунок 3 - График в координатах [(1-у2)"3] - [т] при сорбции толуола в концентрациях от 0,8 до 78,8 % об в растворе н-гептана на цеолите ЫаХ

Анализ изменения адсорбции как функции времени в координатах [(1-у2)"3] -

[т] показывает (рисунок 3), что для рассмотренных процессов это условие выполняется, то есть внутридиффузное сопротивление являегся практически преобладающим в общем сопротивлении сорбции Аналогичные данные получены и для других исследуемых ароматических углеводородов

Анализ кинетики адсорбции из жидкой фазы обычно связан с необходимостью проведения большого объема длительных экспериментов с последующим построением на их базе серии кинетических кривых, обобщающей изотермы избыточной адсорбции Однако для формирования информационного массива по статике адсорбции с позиций дальнейшего расчета динамики сорбции достаточно построить кинетические кривые при двух начальных концентрациях адсорбтива Со при низкой Сон = (1-4 % об) и высокой Сов (50-80 % об), с определением соответствующих равновесных значений величин гиббсовской адсорбции аГн и агв По величине адсорбции аГв при высокой концентрации ароматического углеводо рода Сов (50-80 % об) определяется величина максимальной полной адсорбции Омах при С0 = 100 % об, учитывая градиент адсорбции Эа/сС0 и экстраполируя линию, проходящую через точки с координатами [агв, С0в и 0, 100] (см рисунок 1) Это позволило далее решением системы уравнений

по величинам адсорбции аГн и амах при концентрациях низкой Сон = (1—4) % и Со = 100 %, рассчитать коэффициенты А и В уравнения изотермы адсорбции (уравнение Ленгмюра).

где ад— ленгмюровская (полная) адсорбция, характеризующая величину равновесной адсорбции примеси адсорбентом

artH=ABCj{\ + BCj ашх =АВ 100/(1+5 100)

(3)

(4)

Экстраполяцией прямолинейного участка гиббсовской изотермы на ось ординат получены предельные величины адсорбции для исследованных систем (таблица 1)

Таблица 1 - Параметры изотермы адсорбции ароматических углеводородов

цеолитом ЫаХ из раствора в н-гептане

Адсорбтив Концентрация адсор-бтива, при максимальной аг, С0, % об. Максимальная адсорбционная способность, г г"1 Коэффициенты уравнения изотермы адсорбции (4)

избыточная (гиббсов-ская), аг полная (ленгмю-ровская), ал А, % масс. в, % об''

Бензол 17,3 0,140 0,170 17,5 0,364

Толуол 9,4 0,250 0,285 30,0 0,196

П-ксилол 10,6 0,196 0,225 24,0 0,150

Изопропилбензол 10,0 0,164 0,182 27,0 0,586

Нафталин 7,5 0,120 0,130 13,2 0,774

Как видно из таблицы 1, величина ал бензола занимает промежуточное положение среди исследованных гомологов С одной стороны, наличие алкильных цепей приводит к возрастанию молекулярной массы молекулы адсорбата и, соответственно, приводит к росту активности сорбента, а с другой стороны - к увеличению диффузионного сопротивления при прохождении сорбируемой молекулы в окна пористой структуры цеолита и, как следствие, к снижению его активности

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

г)-(К/5)2, (5)

где Ц; - эффективный коэффициент диффузии, см^/ с, у - величина относительной адсорбции к моменту времени т, т - время адсорбции, с, V - объем гранулы цеолита, см3, 8 - площадь наружной поверхности гранул, см2.

Кроме эффективного коэффициента диффузии Ое, характеризующего сорбцию в микропорах, также рассчитывался коэффициент диффузии бинарной смеси ароматический углеводород - //-гептан (Вр) в межгрануляционном пространстве и широких транспортных порах гранул по уравнению Вильке-Ченга в модификации Шайбеля

БР =8 2-10"8 • ¡[1 + (3 У2/71)3/3]г|/(и2У11/3), (6)

где ц - вязкость к-гептапа, сПз, Т - температура адсорбции, К, VI - молярный объем диффундирующего ароматического углеводорода, см3 моль'1, У2 - молярный объем //-гептана, см3 моль"1 (таблица 2)

Таблица 2 - Коэффициенты диффузии ароматических углеводородов

при их сорбции цеолитом №Х из раствора в н -гептане

Адсорбтив Коэффициент диффузии Б 108, см2 с"'

эффективный в порах Ое в растворе Е>р

максимальный минимальный

Бензол 281,5 31,2 4196,3

Толуол 228,9 31,2 3460,1

П-ксилол 224,5 28,4 2731,8

Изопропилбензол 217,6 31,2 2228,9

Нафталин 193,1 31,2 2714,5

Результаты расчета коэффициентов диффузии показали, что в широком интервале значений концентрации ароматических углеводородов в смеси с н-гептаном и продолжительности адсорбционного процесса зависимость эффективного коэффициента диффузии ароматических углеводородов в гранулах цеолита от времени адсорбции носит экстремальный характер (рисунок 4) При этом меньшие величины коэффициентов диффузии Ос и соответственно большее диффузи-

онное сопротивление в гранулах цеолита достигаются на заключительной стадии сорбционного процесса

Так, например, величина коэффициента диффузии толуола в растворе в н-гептане равна 3460,1 10"8 см2 с"1, т е в 15 раз выше максимально наблюдаемой величины Ое, равной в зоне максимума 228,9 10'8 см2 с'1, и в 110 раз выше асимптотических значений Бе, составляющих около 30 10'8 см2 с"1 Это позволяет считать вполне корректным принятое ранее допущение о том, что лимитирующее диффузионное сопротивление сосредоточено в микропористом пространстве гранул цеолита, а не в объеме потока смеси между гранулами, и пренебречь диффузионными явлениями в межгрануляционном пространстве слоя сорбента при разработки модели динамики процесса

Концентрация -»-0,8 -«-3,8 -*~5,4 -*~9,3 -»-21,3 -°-Зб,8 толуола, % об -Ж-48,8 -»-78,8

Рисунок 4 - График зависимости коэффициента диффузии (Бе) смеси толуол-гептан при концентрациях толуола (0,8 -78,8) % об. от времени

адсорбции (т)

Наличие максимума у зависимости Бе=((т) при сорбции ароматических углеводородов из растворов в н-гептане, по-видимому, можно объяснить наличием воздуха в порах цеолита на начальной стадии сорбции (10-30 минут) Препятствуя транспорту адсорбата через макропоры к микропорам, пузырьки воздуха вносят в процесс дополнительное диффузионное сопротивление По мере вытеснения воздуха величина дополнительного сопротивления уменьшается и Ос возрастает, дос-

тигая максимальной величины, после удаления большей части воздуха эффективный коэффициент диффузии начинает уменьшаться, так как при заполнении ад-сорбатом пористой структуры цеолита растет и диффузионное сопротивление массопереносу С увеличением концентрации адсорбтива позиция максимума у ароматических углеводородов сдвигается правее и выше, что также свидетельствует об активации переноса по мере возрастания количества молекул адсорбтива в исходной смеси

По мере возрастания молекулярного веса адсорбата в гомологическом ряду ароматических углеводородов и увеличения динамической вязкости коэффициент диффузии у гомологов бензола, которые можно рассматривать как бензол с сопряженными от 1 до 4 атомами углерода, уменьшается, то есть возрастает сопротивление массопереносу В области низких концентраций диффузионное сопротивление гомологов бензола при наличии коротких боковых цепей выше, чем у бензола, а в области высоких концентраций - диффузионные сопротивления ароматических углеводородов близки и асимптотически стремятся к некоторой постоянной величине Бс тш при времени сорбции более 100 мин (см рисунок 4)

Исследование кинетики жидкофазной адсорбции ароматических углеводородов цеолитами К'аХ показало, что кинетика сорбции, лимитируемая внутренней диффузией, описывается уравнением Глюкауфа

йа/йг = ра(ал-а), (7)

где а - текущее значение активности адсорбента к моменту времени т от начала сорбции, а«, - равновесное значение гиббсовской адсорбции |30 - коэффициент массопередачи, который можно интерпретировать как функцию отклика гидродинамической модели на импульсное возмущение

Коэффициент массопередачи может быть рассчитан как величина, обратнаг первому статистическому моменту кинетической кривой

Д, = }(1 - а/а. Уг/]г(1-а/<0^ (8)

о о

при следующих допущениях 1) гранула адсорбента является элементарной точкой системы "адсорбент-раствор"; 2) концентрация ароматического углеводорода в

гранулах цеолита в любой момент времени т есть величина постоянная, 3) концентрация ароматического углеводорода в растворе практически не изменяется, несмотря на процесс адсорбции.

Параметр ß0 можно рассматривать в качестве общего эффективного объемного коэффициента массопередачи, который учитывает внутридиффузионное сопротивление в макро- и микропорах адсорбента Определенная таким образом величина ßo практически не зависит от концентрации адсорбтива (рисунок 5), поэтому она может быть использована при моделировании динамики адсорбции как константная характеристика во всем диапазоне изменения концентраций адсорбтива в слое сорбента от сырьевой до проскоковой

ßo, мин

0,03 0,02 0,01 0

0 20 40 60 80 Ч°г ™ ♦ Бензол " Толуол 1 Изопропилбензол

Рисунок 5 - Зависимость общего объемного коэффициента массопередачи ß0 от концентрации извлекаемого компонента в сырье С0 при сорбции ароматических углеводородов цеолитом NaX из раствора в н-гептане

Отношение внутреннего сопротивления массопереносу в зерне адсорбента к внешнему сопротивлению из жидкости путем диффузии через пограничный вязкий слой к поверхности зерна определялось с использованием критерия Био (BiD)

BiD={t2-ß0)/D3, (9)

где I - характерный линейный размер зерна (эффективный радиус), D3- коэффициент диффузии молекул внутри зерна адсорбента

Результаты расчетов (таблица 3) еще раз подтверждают тот факт, что скорость адсорбции ароматических углеводородов в растворе /i-гептана на цеолите

№Х зависит практически только от диффузии молекул внутри адсорбента (В1 > 30), т с контролируется стадией внутридиффузионного массопереноса

Таблица 3 - Расчет критерия Био при адсорбции ароматических

углеводородов в н-гептане на цеолите NaX

Адсорбтив Начальная концентрация С,% об Критерий Био (Вгц)

при т = 2 мин при

максимальном минимальном

Бензол 1,4 1846,795 1266,792 5135,587

17,28 1972,130 918,871 6769,638

75,53 2206,697 953,293 8403,688

Толуол 3,8 5163,672 1429,390 6827,996

21,3 18830,487 1322,640 7878,457

78,8 3419,469 976,797 7878,457

п-Ксилол 2,96 1299,945 1299,945 5777,536

40,56 2095,823 997,025 6040,151

81,85 3698,498 1062,481 7878,458

Изопропилбензол 2,05 2453,993 1641,346 6565,381

39,88 5994,538 1091,938 7615,843

80,45 8205,745 890,379 6827,997

Нафталин 2,45 1575,338 1451,762 6565,381

8,3 4672,990 1196,285 8141,073

В третьей главе рассмотрена реализация алгоритма моделирования динамики жидкофазной адсорбции в неподвижном слое сорбента

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

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

Такое внутридиффузионное сопротивление характерно, что было подтверждено во второй главе диссертации, для процесса деароматизации растворов углеводородов цеолитами В этом случае внешнедиффузионным сопротивлением можно пренебречь, что позволяет при разработке модели динамики адсорбции пользоваться информацией, получаемой при исследовании кинетики и статики адсорбции

Исходными данными для расчета динамики сорбции при разделении бинарной смеси, например, углеводородов цеолитами в сорбционной колонне высотой К (см) и диаметром Б (см) с загрузкой адсорбента в (г) при расходе сырья V (см3/мин) с концентрацией адсорбтива С0 (% об), являются значения коэффициентов А (% масс) и В (% об"') уравнения Ленгмюра, коэффициента массопередачи р (мин"1), доля свободного пространства в слое сорбента с учетом порозности слоя сорбента и удельного объема транспортных пор и микропор в гранулах цеолита, а также скорость потока чу (см/мин) в расчете на полное сечение колонны

Полученный в процессе выполнения экспериментов по кинетике и статике адсорбции информационный базис позволяет моделировать динамику адсорбции в неподвижном слое сорбента на основе уравнений материального баланса, решаемых методом численного интегрирования В этом случае слой адсорбента разбивается по высоте на N элементарных ячеек толщиной ДЬ, в каждой из которых просчитывается отработка соответствующего участка слоя за исследуемую (за данную) продолжительность сорбционного процесса т (мин) в течение малых интервалов времени Ат = ДЪЛу В связи с дискретностью расчета процесса величина ДЬ составляет 1 см (несколько диаметров гранул адсорбента), величина Дт при этом обычно соответствует (0,5-2) мин, составляя 0,2-0,5 % от исследуемого временного диапазона и обеспечивая выполнение не менее 200-500 шагов интегрирования в ходе расчета каждой адсорбционной ячейки При этом движение очищаемого потока по колонне в целом характеризуется режимом идеального вытесне-

ния, а в пределах одной ячейки - режимом идеального смешения Расчеты г-той сорбции за время сорбционного процесса т позволяют сформировать массив конечных значений количества жидкости и концентрации в ней адсорбируемого вещества на выходе из ячейки с малыми интервалами времени Дт. Последующая переадресовка массива конечных значений параметров ¿-той ячейки (расход потока жидкости и концентрация в ней адсорбируемого компонента) в массив начальных значений соответствующих параметров (г+1)-вой ячейки, распределяемых во времени процесса т, позволяет выполнить расчет по выше рассмотренному алгоритму (;+1)-вой ячейки Аналогично выполняется расчет всех ячеек адсорбционной ко лонны

В ходе расчета динамики сорбции последовательно для каждой ячейки рассматриваются три задачи

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

2) выход первой порции очищаемого продукта из полностью заполненной ячейки (при этом количество выходящего продукта меньше количества поступающего в ячейку продукта),

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

Проверка разработанной модели осуществлялась путем расчета динамики адсорбции бензола из раствора в я-гептане цеолитами КаХ Расчетные параметры сравнивались с экспериментальными данными по адсорбции бензола из раствора в н-гептане цеолитами №Х в динамических условиях Результаты моделирования процесса адсорбции бензола из раствора в //-гептане показали, что разработанный метод расчета динамики сорбции позволяет удовлетворительно как качественно, так и количественно рассчитать позицию адсорбционного фронта в слое адсорбента и выходную кривую сорбции- изоплану адсорбции, по которой рассчитываются основные-характеристики адсорбционного процесса динамическая активность ад-

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

Объективную оценку адекватности решения задачи по разработанной модели выполняли по статистическому критерию Фишера. Из результатов оценки адекватности следует, что расчет по предлагаемой математической модели динамиче-скои активности адсорбента до условного проскока бензола в рафинат ав, полной динамической активности адсорбента аЕ и длины зоны массопередачи Ь0 является адекватным с уровнем вероятности ошибочности принятия гипотезы об адекватности модели 0,05 (таблица 4)

Успешное практическое использование разработанной модели динамики сорбции зависит от точности экспериментальных данных о кинетике адсорбции Поскольку используемые при моделировании параметры модели определяются по экспериментальным данным с неизбежной погрешностью, то в диссертации была выполнена оценка параметрической чувствительности математической модели для четырех, используемых в расчете, параметров коэффициента массопередачи Р, констант А и В уравнения изотермы адсорбции Ленгмюра и суммарной порозно-сти слоя адсорбента е. Варьирование параметров А, В, е и Р в пределах ±10 % относительно базовой величины выполнялось по очереди

Таблица 4 - Результаты оценки адекватности модели динамики адсорбции

Уровень вероятности ошибочности принятия гипотезы об адекватности модели Критерий Фишера

табличное значение расчетное значение для

проскоко-вой активности полной активности длины зоны массопередачи

0,05 234 166,6 4,67 13,8

Установлено, что при изменении кинетических параметров (коэффициента массопередачи, константы В уравнения изотермы адсорбции Ленгмюра и суммар ной порозности слоя адсорбента) на 10 % параметры динамики сорбционного процесса (проскоковая динамическая активность адсорбента, полная динамиче-

екая активность адсорбента, фактор симметричности и длина зоны массопередачи) изменяются не более чем на 7 % Исключением является константа А уравнения Ленгмюра, при изменении которой на 10%, динамические параметры процесса могут изменяться на 13,5 % относительных

Четвертая глава посвящена разработке математической модели оптимизации работы адсорбционной установки на основе предварительных экспериментов на лабораторной или пилотной аппаратуре при параметрах существенно отличающихся от промышленных

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

К = За+Зп+Зм+Зр, (10)

где За - затраты на адсорбент, руб /год, Зп - затраты на перекачку, руб /год, Зм - затраты на металл конструкций, руб /год, Зр - затраты на регенерацию, руб /год

На основе зависимости длины зоны массопередачи от размеров слоя сорбента и условий его регенерации разработана целевая функция, минимизация которой осуществлялась методом неопределенных множителей Лагранжа Дифференцирование функции Лагранжа по оптимизируемым параметрам (диаметру аппарата, его высоте, температуре регенерации адсорбента) и множителю Лагранжа дало ряд нелинейных уравнений, преобразование которых дало обобщенное уравнение, позволившее путем закрепления значения диаметра аппарата вычислить реальные значения остальных параметров высоты аппарата и температуры регенерации Расчеты производились для серии значений диаметра от 0,4 до 6 м при использовании двух видов топлива (мазута и метана) при регенерации адсорбента

Разработанная технико-экономическая модель была проверена путем сравнения расчетных данных, полученных на основе лабораторных опытов по осушке изопрена цеолитами №А, и материалов действующего объекта - адсорбционной установки ЦГФУ Нижнекамского нефтехимического комбината, состоящей из двух аппаратов производительностью 125 м3/ч и предназначенной для глубокой осушки углеводородов С2-С7 цеолитами №А в жидкой фазе Диаметр одного про-

мышленного аппарата 2,6 м, высота слоя -8м Продолжительность стадии адсорбции составляет 20 суток, остановка на плановый ремонт - 31 день

Из полученных результатов расчетов установлено, что оптимальными параметрами для этой установки являются диаметр аппарата - 3,4 м, высота слоя -6,22 м, температура регенерации - 339,74°С, при этом значение критерия оптимальности при использовании в качестве топлива при регенерации адсорбента метана равно 2086 тыс руб /год (рисунок 6), а также что при регенерации адсорбента экономически более выгодно использовать в качестве топлива газ внутризавод-

V

ского производства

Диаметр, м 1 - общие затраты, 2 - высота аппарата

Рисунок 6 - Зависимость общих затрат и высоты аппарата от его диаметра при использовании в качестве топлива метана

Выводы

1 Исследована кинетика сорбции ароматических углеводородов (бензола, толуола, п-ксилола, изопропилбензола и нафталина) из раствора в н-гептане цеолитом КтаХ при температуре 20°С в диапазоне концентраций 1-80 % об. и получена информационная база данных (статических и кинетических характеристик процес-

са адсорбции), необходимая для моделирования динамики сорбции ароматически: углеводородов цеолитами ЫаХ

2 Рядом методов расчетного анализа обосновано, что процесс жидкофазной сорбции ароматических углеводородов цеолитами №Х лимитируется внутридиф-фузионным сопротивлением в пористой структуре сорбента, которое на 1-2 порядка превышает сопротивление транспорту ароматических углеводородов в макропорах

3 Установлено, что общий объемный коэффициент массопередачи, рассчитанный как величина, обратная первому статистическому моменту кинетической кривой статики адсорбции, является параметром, независимым от концентрации сорбируемой примеси ароматического углеводорода в бинарной смеси с парафиновым растворителем

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

5 Выполнена оценка параметрической чувствительности математической модели динамики жидкофазной адсорбции бинарной смеси углеводородов, показавшая, что погрешность определения кинетических параметров (коэффициент массопередачи, константы А и В уравнения изотермы адсорбции Ленгмюра и суммарная порозность слоя адсорбента) оказывается выше погрешности расчета параметров сорбционного процесса за исключением константы А уравнения Ленгмюра

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

7 Определены оптимальные размеры адсорбера для процесса жидкофазной осушки углеводородов - сырья ЦГФУ Нижнекамского нефтехимического комбината

Содержание работы опубликовано в 13 научных трудах, из них 2 статьи опубликованы в перечне ведущих рецензируемых научных журналов и изданий в соответствии с требованиями ВАК Минобрнауки России

1 Тебенькова, ЕБ Исследование кинетики адсорбции бензола цеолитами NaX из раствора в н гептане / Е Б Тебенькова, Н В Заиченко, Н А. Самойлов // 56-я научно-техническая конференция студентов, аспирантов и молодых ученых сб тез докл / УГНТУ - Уфа, 2005 - Кн 2 - С 186

2 Самойлов, Н А. Параметры кинетики жидкофазной адсорбции ароматических углеводородов цеолитами NaX /НА Самойлов, Н В Заиченко, Е Б Тебенькова // Нефтегазопереработка и нефтехимия, тез. докл VI конгресса нефтегазо-промышленников/ ГУП ИНХП РБ - Уфа, 2005,- С. 400-401.

3 Самойлов, Н А Особенности диффузии ароматических углеводородов при их сорбции цеолитами их растворов /НА Самойлов, Н В Заиченко // Нефтегазовые и химические технологии тез докл III Всерос. науч -практ конф./ СГТУ -Самара, 2005 - С 124

4 Самойлов, IIА Влияние концентрации адсорбируемой примеси на коэффициент массопередачи в кинетике адсорбции углеводородов цеолитами из жидкой фазы /НА Самойлов, Н В Заиченко, Е Б Тебенькова // Термодинамика поверхностных явлений и адсорбции тез докл Всерос семинара / ИГХТУ - Иваново, 2005-С 16-17.

5 Самойлов, HAK аналитическому расчету изотерм избыточной и полной адсорбции /НА Самойлов, Н.В Заиченко, Е Б Тебенькова // Термодинамика поверхностных явлений и адсорбции тез, докл Всерос. семинара / ИГХТУ.- Иваново, 2005.-С 17-18

6 Тебенькова, Е Б Исследование кинетики адсорбции нафталина цеолитами NaX из раствора в н гептане / Е Б Тебенькова, Н В Заиченко,-Н А Самойлов // 57-

я научно-техническая конференция студентов, аспирантов и молодых ученых сб тез. докл / УГНТУ.- Уфа, 2006 - Кн 2 - С. 8

7 Тебенькова, Е Б Исследование кинетики адсорбции изопропилбензола цеолитами КтаХ из раствора в н гептане / Е Б Тебенькова, Н В Заиченко, Н А Самойлов / / 57-я научно-технической конференции студентов, аспирантов и молодых ученых сб тез докл / УГНТУ -Уфа, 2006 - Кн 2 - С 9.

8 Тебенькова, Е Б Исследование кинетики адсорбции толуола цеолитами ЫаХ из раствора в п гептане / Е Б Тебенькова, Н В Заиченко, Н А Самойлов // 57-я научно-техническая конференция сгуденгов, аспирантов и молодых ученых сб тез докл / УГНТУ - Уфа, 2006 - Кн 2 - С 10

9 Тебенькова, ЕБ Исследование кинетики адсорбции п-ксилола цеолитами №Х из раствора в н гептане / Е Б Тебенькова, Н В Заиченко, Н А Самойлов // 57-я научно-технической конференция студентов, аспирантов и молодых ученых- сб тез докл/УГНТУ-Уфа, 2006-Кн 2 - С 11

10 Заиченко, НВ Лимитирующее диффузионное сопротивление при жидко фазной адсорбции ароматических углеводородов цеолитами ИаХ / Н.В. Заиченко, Н А Самойлов // Перспективы развития химической переработки горючих ископаемых тез докл. конф / ХПГИ. - СПб, 2006 - С. 260.

11 Заиченко, Н В Оценка диффузионного сопротивления при жидкофазной адсорбции ароматических углеводородов цеолитами ЫаХ / НВ Заиченко, НА Самойлов // Нефтегазопереработка и нефтехимия тез докл междунар науч -практ конф / ГУЛ ИНХП РБ - Уфа, 2006 - С 272-273

12 Заиченко, НВ Параметрическая чувствительность модели динамики жидкофазной адсорбции / НВ, Заиченко, Н.А. Самойлов // Теория и практика массо-обменных процессов химической технологии (Марушкинские чтения) тез докл III Всерос науч конф / УГНТУ.-Уфа, 2006.-С. 16-18

13 Заиченко, Н В Математическая модель динамики адсорбции / Н В Заиченко, Н А. Самойлов / / Переработка углеводородного сырья Комплексные решения (Левинтерские чтения) тез докл. Всерос науч конф - Самара, 2006.-С. 143-144

14 Самойлов, H.A. Моделирование динамики адсорбции цеолитами из жидкой фазы на основе данных статики адсорбции /НА Самойлов, Н В Заиченко /' Журн. прикладной химии - 2006 - Т 79, № 7 - С 1137-1143

15 Самойлов, H.A. Моделирование оптимальных условий реализации адсорбционного процесса /НА Самойлов, Д Н. Митягин, Н В Заиченко // Актуальные проблемы теории адсорбции, пористости и адсорбционной селективности тез докл XI Всерос симпозиума с участием иностр ученых - Москва-Клязьма, 2007 -С 44

16 Митягин, Д Н Разработка метода оптимизации адсорбционной установки / Д Н Митягин, Н.В Заиченко, Н А Самойлов // 58-я научно-техническая конференция студентов, аспирантов и молодых ученых сб тез докл / УГНТУ - Уфа, 2008 - Кн 2 - С 28

17 Заиченко, Н В Кинетика адсорбции ароматических углеводородов из растворов в н гептане цеолитами NaX / Н В Заиченко, Н А Самойлов // Изв вузов Нефть и газ - 2008 - №3 - С 65-73

Подписано в печать ZO.OG.oi Бумага офсетная Формат 60x84 1/16 Гарнитура «Тайме» Печать трафаретная Уел печ л. 1 Тираж 9 0 Заказ 9$.

Типография Уфимского государственного нефтяного технического университета

Адрес типографии 450062, Республика Башкортостан, г Уфа, ул Космонавтов,!

Введение.

Глава 1. Литературный обзор.

1.1 Статика и кинетика адсорбции.

1.1.1 Адсорбционные силы. Типы адсорбции.

1.1.2 Кинетическая кривая адсорбции.

1.1.3 Изотермы адсорбции. Теории адсорбции.

1.1.4 Коэффициент диффузии.

1.1.5 Энергия активации диффузионного процесса. Теплота адсорбции

1.1.6 Адсорбция из жидких сред.

1.1.7 Коэффициент массопередачи.

1.2 Динамика адсорбции и десорбции. Моделирование кинетики и динамики процессов адсорбции. Оптимизация адсорбционных процессов.

1.2.1 Динамика адсорбции.

1.2.2 Моделирование кинетики и динамики процессов адсорбции.

1.2.3 Десорбция.

1.2.4 Оптимизация адсорбционных процессов.

1.3 Выводы по главе.

Глава 2. Определение характеристик кинетики адсорбции ароматических углеводородов из н-гептана цеолитами NaX.

2.1 Методика проведения экспериментов.

2.2 Экспериментальные данные и обсуждение результатов.

2.3 Выводы по главе.

Глава 3. Моделирование динамики адсорбции ароматических углеводородов из н-гептана цеолитами NaX.

3.1 Особенности разработки модели динамики сорбции.

3.2 Оценка параметрической чувствительности модели динамики адсорбции.

3.3 Выводы по главе.

Глава 4. Оптимизация процессов жидкофазной адсорбции на цеолитах

4.1 Разработка модели оптимизации адсорбционной колонны.

4.2 Выводы по главе.

Выводы.

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

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

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

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

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

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

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

И, наконец, практически отсутствуют публикации по оптимизации адсорбционных процессов.

Отсутствие целостной теории о процессе сорбции затрудняет поиск оптимального технологического режима:

На основании анализа публикаций последних 40 лет можно определить перспективные направления исследования, сформулировать цель работы:

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

- оптимизация работы адсорбционной колонны.

Для достижения поставленной цели необходимо решить следующие задачи:

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

- оценить характер лимитирующего диффузионного сопротивления при жид-кофазной адсорбции на цеолитах;

- получить кинетические параметры, необходимые для расчета динамики жидкофазной адсорбции;

- разработать метод расчета динамики адсорбции органических примесей цеолитами на основе минимального объема экспериментальных данных о кинетике адсорбции; доказать адекватность расчетов по разработанной модели данным контрольных экспериментов по динамике жидкофазной адсорбции; выполнить оценку параметрической чувствительности математической модели;

- разработать метод оптимизации адсорбционной установки, учитывающий затраты на адсорбционный аппарат: затраты на адсорбент, перекачку, металл конструкций и регенерацию.

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

1) показано, что общий объемный коэффициент массопередачи, рассчитанный как величина, обратная первому статистическому моменту кинетической кривой статики адсорбции, является параметром, независимым от концентрации сорбируемой примеси ароматического углеводорода в бинарной смеси с парафиновым растворителем;

2) общий объемный коэффициент массопередачи использован при разработке математической модели динамики жидкофазной сорбции ароматических углеводородов в гептане;

3) выполнена оценка параметрической чувствительности математической модели динамики жидкофазной адсорбции бинарной смеси углеводородов;

4) разработан метод оптимизации работы адсорбционной установки с учетом затрат на основное оборудование, стоимость перекачки сырья и регенерацию адсорбента.

выводы

1. Исследована кинетика сорбции ароматических углеводородов (бензола, толуола, п-ксилола, изопропилбензола и нафталина) из раствора в //-гептане цеолитом NaX при температуре 20°С в диапазоне концентраций (1—80) % об. и получена информационная» база данных (статических и кинетических характеристик процесса адсорбции), необходимая для моделирования-динамики сорбции ароматических углеводородов цеолитами NaX.

2. Рядом методов расчетного анализа обосновано, что процесс жидкофазной сорбции ароматических углеводородов цеолитами NaX лимитируется внутридиффу-зионным сопротивлением в пористой структуре сорбента, которое на 1-2 порядка превышает сопротивление транспорту ароматических углеводородов в макропорах.

3. Установлено, что общий объемный коэффициент массопередачи, рассчитанный как величина, обратнаяiпервому статистическому моменту кинетической кривой статики адсорбции, является^ параметром, независимым от концентрации сорбируемой примеси ароматического углеводорода в бинарной смеси с парафиновым растворителем.

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

5. Выполнена оценка параметрической чувствительности математической модели динамики жидкофазной-адсорбции бинарной смеси углеводородов, показавшая, что погрешность определения кинетических параметров (коэффициент массопередачи, константы ,4 и В уравнения изотермьь адсорбции Ленгмюра и суммарная порозность слоя адсорбента) оказывается выше погрешности расчета параметров сорбционного процесса за исключением константы А уравнения1. Ленгмюра. ' t

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

7. Определены оптимальные размеры адсорбера для процесса жидкофазной осушки углеводородов - сырья ЦГФУ Нижнекамского нефтехимического комбината.

1. Кельцев, Н.В. Основы адсорбционной техники / Н.В. Кельцев. — М.: Химия, 1984-592 с.

2. Тимофеев, Д. П. Кинетика адсорбции / Д.П. Тимофеев. — М.: Изд-во АН СССР, 1962-252 с.

3. Моделирование внутридиффузионной кинетики адсорбции на цеолитах / С.Г. Разборова, Е.А. Устинов, В.Ю. Челышев и др. // Журн. прикладной химии.-1990.-Т. LXIII,№3.- С. 588-591.

4. Беринг, Б.П. Особенности адсорбции смесей паров на микропористых адсорбентах / Б.П. Беринг, В.В. Серпинский, С.И. Суринова // Цеолиты, их синтез,свойства и применение: Мат. II Всесоюзного совещания по цеолитам.- М., 1965. —1. С. 30-36.

5. Равновесная адсорбция азота, кислорода и их смеси на цеолитах / А.К. Акулов, А.К. Агабалян, Е.А. Устинов, С.Ф. Гребенников // Журн. прикладной химии.- 1987.-Т. LX, №8.- С. 1707-1710.

6. Ritter, James F. Equilibrium theory for hysteresis-dependent fixed-bed;adsorptions / J. F. Ritter, R. T. Yang // Chemical Engineering Science.- 1991.-Vol. 46, No2.- P. 563-574.

7. Якимова, A.B. О расчете коэффициентов диффузии при жидкофазной адсорбции ш-парафинов цеолитами1 / А.В: Якимова, Астахов;, ВЗД: Лукин // Журн. прикладной химии.- 1983.-Т. LVI, №7.- С. 1488-1491.

8. Пономарев, А С. Исследование скорости стационарной диффузии органических павров в стержнях из спрессованного цеолита NaX / АС. Пономарев,

9. Д.П. Тимофеев // Цеолиты, их синтез, свойства и применение: Мат. II Всесоюзноiго совещания по цеолитам. М., 1965. — С. 338-340. .

10. Кожушко,чВ;В1, Метод расчета коэффициентов^ диффузииша основе кинетических кривых адсорбции / В.В: Кожушко, А.В. Новосельский,. В.Д. Лукин // Журн. прикладной,химии.- 1988.-гГ. LXI, №2.- С. 430-432. ■

11. Ерашко; И:Т. Исследование кинетики сорбции н-пентана и диэтилового эфира на цеолитах СаХ и NaX / И.Т. Ерашко // Цеолиты, их синтез, свойства и применение: Мат. II Всесоюзного совещания по цеолитам. Mi, 1965: - С. 340-■345.

12. Аксельруд, Е.А:.Исследование кинетикшадсорбции; бензойной; и серной? кислот из; водных растворов- / Е.А. Аксельруд, А.Е. Кащеев; MtB. Левицкая? // Журн. прикладной химии.- 1990:-T.,EXIII; №5:-С. 1092-1095:

13. Carleton; F.D. Adsorption in non-isobaric fixed beds/ F.D. Carleton, L,S. Ker-shenbaum, W.A. Wakeham // Chemical Engineering Science.- 199E-Vol.33; No 9.- P. 1239-1246.

14. An analysis of a nonisothermal one-component sorption in a single adsorbent particle / A. Brunovska, V. Hlavacek, J-. Ilavsky, J. Valtyni //Chemical Engineering Science.- 1978.-Vol. 33, No 10.-P. 1358-1391.

15. Массообмен при периодической адсорбции из растворов / Г.А. Аксель-руд, М. Б. Левицкая, В.Д. Степанов, А.Е. Кащеев// Журш прикладной химии,-1988.-Т. LXI, №3.- С. 658-661.

16. Кац, Б.М. Влияние предадсорбции паров органических соединений на поглотительную способность макропористых анионитов по диоксиду углерода / Б.М. Кац, Е.К. Малиновский, Г.А. Артюшин // Журн. прикладной химии.-' 1986.-Т. LIX, №1.- С. 120-123.

17. Жданов, А.А. К вопросу об адсорбции газов и паров / А.А. Жданов, Н.Т.

18. Селиванов//Журн. прикладной химии.- 1984.-Т. LVII, №8.- С. 1864-1867. 1

19. Кутлаева, Т.В. Влияние природы обменных катинов цеолитов на теплоту адсорбции диоксида серы и аммиака / Т.В. Кутлаева, С.А. Ануров, Д.А. Лобанов // Изв. вузов. Химия и хим. технология,- 1992.- Т. 35, № 6.- С. 115-116.

20. Устинов, Е.А. Закономерности динамики циклических адсорбционных процессов очистки при больших длинах слоя / Е.А. Устинов, А.К. Акулов // Журн. прикладной химии,- 1981.- Т. LIV, №6.- С. 1270-1276.

21. Адсорбция органических веществ из воды / A.M. Когановский, Н.А. Клименко, Т.М. Левченко, И.Г. Рода. Л.: Химия, 1990. — 256 с.

22. Когановский, A.M. Кинетика адсорбции органических веществ и йх смеIсей из водных растворов / A.M. Когановский, P.M. Марутовский // Адсорбция иадсорбенты: Тр.шестой Всесоюз. конф. по теоретическим вопросам адсорбции. i

23. М.: Наука, 1987.-С. 135-141.

24. Heink, W. Application of zeugmatography to study kinetics of physical adsorption / W. Heink, J. Karger, H. Pfeifer // Chemical Engineering Science.- 1978.-Vol. 33, No 8.- P. 1019-1023.

25. Buzanowski, Marek A. Approximations for intraparticle diffusion rates in cyclic adsorption and desorption / Marek A. Buzanowski, Ralph T. Yang // Chemical Engineering Science.- 1991.-Vol. 46, No 10.- P. 2589-2598.

26. Эльтекова, Н.А. Термодинамические исследования состояния адсорбированного вещества / Н.А. Эльтекова, Ю.А. Эльтеков. // Журн. физ. химии,- 1995.-Т. XXXIX, №6.- С. 22-33.

27. Эльтекова, Н.А. Особенности жидкостной хроматографии при отрицательных значениях/Н.А. Эльтекова, Ю.А. Эльтеков // Журн. физ. химии.- 1996.-Т. 70, №3.-С. 532-538.I

28. Акулов, А.К. Динамика короткоцикловой безнагревной адсорбционной осушки воздуха на цеолитах/ А.К. Акулов // Журн. прикладной химии.- 1989.-Т. LXII, №10.- С. 2258-2262.

29. Устинов, Е.А. Динамика адсорбции воды цеолитом в процессах разделения воздуха методом PSA/ Е.А. Устинов, Ю.Г. Карандашова // Журн. прикладной химии.- 1994.-Т. 67, №6.- С. 949-952.

30. Юсубов, Ф.В. Идентификация параметров внешней и внутренней диффузионной динамики адсорбции/ Ф.В. Юсубов, Р.И. Зейналов, Ч.Ш. Ибрагимов //I

31. Журн. прикладной химии.- 1994.-Т. 67, №7.- С. 1121-1123.

32. Карандашова, Ю.Г. Моделирование равновесной адсорбции газов цеолитом в присутствии воды/ Ю.Г. Карандашова, Е.А. Устинов, Н.Ф. Федоров // Журн. прикладной химии.- 1994.-Т. 67, № 6.- С. 1032-1034.

33. Титов, А.А. Исследование структуры барботажного слоя и массообмена в трехфазной системе при адсорбции на провальных щелевых тарелках/ А.А. Титов, Э.М. Ганбаров //Хим. промышленность.- 1995.- № 8,- С. 426-431.

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

35. Ауашриа, G. Расчет коэффициента массоотдачи в жидкой фазе на струй-но-направленных тарелках с продольными перегородками/ Ауашриа С. //Химия и тёхнология;тЬплив;№масел.-4993г.г-№-Г1.-;С.1:5-1'6:.

36. Dudukovic, М.Р. A zone model for reactions of solid particles with strongly adsorbing species / M.P. Dudukovic, M.S. Lamba //Chemical Engineering Science.-1978.-Vol. 33, No 4.- P. 471-478.

37. Определение истинных коэффициентов массоотдачи в насадочных адсорбционных аппаратах/ И.В. Макина, Н!ЩРельперин; ВСоколов; B.JI. Новобратский //Хим. промышленность.- 1986,- Лг« 11.- С. 679-681.

38. Норкина, Л.А. Закономерность изменения коэффициента массопередачи в зависимости от константы фазового равновесия в абсорбере Вентури/JI.A. Норкина, Б.П. Волгин, Ф'С. Югай //Журн. прикладной химии.- 1980.- т. ЫЩ№12.- С.2703-2706. 'j

39. Супрун, В.М. Механизм массоотдачи от поверхности к псевдоожижен-ному слою/ В.М. Супрун //Хим. промышленность.- 1982.г№ 6.- С. 368-369.

40. Акулов, А.К. Кинетика массопереноса при адсорбционном разделении газовых смесей/ А.К. Акулов, Е.А. Устинов//Журн. прикладной химии.- 1985.-Т. LVIII, №7.- С. 1517-1522.

41. Ikeda, Kenji. Performance of the nonisothermab fixed bed adsorption column with nonlinear isotherms^ / Kenji-1 Ikeda // Chemical Engineering Science.- 1979.-Vol. 34, No 7,-P. 941-949.

42. Do, D.D. Hierarchy of rate models; for adsorption and? desorption in bidisipersed structured sorbents / D:D. Do // Chemical Engineering Science.- 1990--Vol. 45, No 5.- P. 1373-1381. •

43. Hills, J:H; Non-isothermal adsorption in a pellet / J.H. Hills // Chemical -Engineering Science.- 1991.-Vol. 46, No 1.- P. 69-74.

44. Лукин, В.Д. Математическая модель процесса адсорбции в аппарате с центробежным разделением фаз/ В.Д. Лукин; H.IT. Чусов // Журн. прикладной химии.- 1981.- Т. LIV, №1.- С. 102-107.

45. Беспалов, А.П. Адсорбция аммиака на эластичной тарелке провального типа большого свободного сечения/ А.П. Беспалов, Е.Б. Зотов, B.C. Лебедев // Журн. прикладной химии.- 1991.- Т. 64, №10.- С. 2199-2200. '

46. Исследование непрерывных процессов адсорбции в аппаратах с движущимися слоями/ А.А. Себалло, Е.А. Устинов, Р.А'. Багиров и др. // Журн. прикладной химии.- 1984.-T.LVII, №5.- С. 1053-1057.

47. Pressureswing- adsorption for a system, with a Langmuir-Freundlich isotherm / X. Lu, X. Madey, D. Rothstein, M. Jaroniec // Chemical Engineering Science.- 1997.-Vol. 45, No 4.-P. 1097-1103.

48. Gladden, L.F. Influence of pellet structure on selectivity during pressure swingadsorption separations / L.F. Gladden // Chemical Engineering Science.- 1991.-Vol. 46,i1. No 10.-P. 2455-2464.

49. Farooq, S. Numerical simulation of a kinetically controlled pressure swing adsorption bulk separation process based on a diffusion model / S. Farooq, D.M. Ruthven // Chemical Engineering Science.- 1991.-Vol. 46, No 9.- P. 2213-2224.

50. Влияние внутренней диффузии на массопередачу в процессах адсорбции / Р.Я. Фишер, В.М. Кисаров, Р.В. Корчагина, А.И. Субботин // Журн. прикладной химии.- 1981.- Т. LIV, №12.- С. 2701-2704.

51. Математическая модель сорбции во'взвешенном слое в системе жидкость-твердое тело / В.П: Кононов, И.О. 11ротодьконов, Н.А. Марцулевич, П:Г. Романков //Журн. прикладной химии.- 1980.- Т. LIII, №11.- G. 2562-2565:

52. Moon, Нее: Adsorption of unknown substances from aqueus solutions / Нее Moon, Heung Chul Park, Ohi Tien // Chemical Engineering Science.- 1991.-Vol*. 46, No 1.- P. 23-31.

53. Gorius, A. Perturbative solutions for non-linear Fixed-bed adsorption-II. Approximate analytical expressions for the dynamics of micro-columns / A. Gorius, M. Bailly, D; Tondeurv// Chemical'Engineering Science.- 199L-VoL 46,,No*2.- P. 685-692.

54. A simplified competitive equilibrium, adsorption model' / F.A. Digiano, G.Baldauf,. B. Frick, H. Sontheimer // Chemical Engineering Science.- 1978.-Vol. 33, No 12.-P. 1667-1673.

55. Samoilov, N. A. Mathematical model of temperature's anomalite at coadsorpthtion of binary mixture / N.A. Samoilov // 13 International Zeolite Conference: Recent Research Reports: French Zeolite Group.- Montpellier, 2001.- 16-R-01

56. Localization of water molecules and sodium ions in Na-modernite by thermally stimulated current mtasurement / S .Devautour, G: Maurin, A. Abdoulaye, F.• th

57. Henn, J.C. Giuntini // 13 Internationale Zeolite Conference: Recent Research Reports: French Zeolite Group.- Montpellier, 2001.- 17-R-02

58. Давид, И. Уравнение диффузии в случае переменного коэффициента диффузии / И. Давид, 3. Гропшяну, М. Неагу // Журн. прикладной химии.- 1986.Т. LIX, JV<> 9.- С. 2044-2047.

59. Массообмен при адсорбции в жидкостно-газовых струйных аппаратах с удлиненной камерой смешения/ А.Н1 Блазнов^ В;А. Куничащ,Д-В;.Чащилов| Ю;Н; Денисов//Журн: прикладной химии;-2003.-Т. 76; №1.-С. 101-106:

60. Самойлов, Н.А. Влияние высоты слоя цеолита NaX на его адсорбционную емкость при лабораторных исследованиях / Н.А. Самойлов; Л:Ф. Фоминых //Технология нефти и газа: Вопросы,фракционирования/ УНИ;-М., 1967.- О. 79-83.

61. Захарова, 11.В. Разделение амиленов изо- и нормального строения на синтетических цеолитах / Н;В. Захарова, B.C. Виноградова, Л.С. Кофман // Цеолиты, их синтез, свойства и применение:,Мат. II Всесоюзного совещания по цеолитам. М., 1965-. - С. 274-276.

62. Процесс очистки изобутилена адсорбцией синтетическими цеолитами / B.C. Виноградова, JI.O: Дымент, В.Н. Комарова и др. // Цеолиты, их синтез, свойства и применение: Мат. II Всесоюз. совещания по цеолитам. М., 1965. — С. 277283

63. Фоминых, Л.Ф. Определение длины зоны массопередачи при адсо'рбции /I

64. Л.Ф: Фоминых, Н;А. Самойлов //Технология нефти и газа: Вопросы фракциони-рования/УНИ .- Уфа, 1975'.- Вып. 3.- С. 309-316:

65. Тонкая очистка алканов С4-С6 от сернистых соединений / В.Г. Урбан,

66. Золотарев, П.П. Проблемы теории динамики адсорбции в неподвижныхIслоях / П.П. Золотарев // Адсорбенты, их получение, свойства и применение: Тр. V Всесоюз. совещания по адсорбентам. М.,1985 - С. 85-91.

67. Filippov, L.K. Technological computation of frontal modes in adsorption separation of multicomponent mixtures / L.K. Filippov, I.V. Filippova // Chemical Engineering Science.- I990.-Vol. 45, No 4.- P. 927-933.

68. Gorius, A. Perturbative solutions for non-liner fixed-bed. adsorption-I. Approximate analytical solutions for asymptotic fronts / A. Gorius, M. Bailly, D. Tondeur //. Chemical Engineering Science.- 1991.-Vol. 46, No 2.- P. 677-684.t

69. Frey, Douglas D. The entropy cjndition* for the dynamics of nonlinear multi-component sorption in porous media / Douglas D. Frey // Chemical Engineering Science.- I990.-Vol. 45, No l.-P. 131-142.

70. Акулов, А.К. Организация циклического адсорбционного процесса очистки сорбирующегося газа от близкой по сорбируемости примеси / А.К. Акулов, Е.А. Устинов // Журн. прикладной химии.- 1986.-Т. LIX, №3.- С. 583-590.

71. Акулов, А.К. Кинетика и динамика адсорбционного разделения'бинарной смеси газов / А.К. Акулов// Журн. прикладной химии.- 1988.-T.LXI, №7.- С. Г525-1529.

72. Акулов, А.К. Динамика короткоцикловой безнагревной адсорбционной осушки воздуха на цеолитах / А.К. Акулов // Журн. прикладной химии.- 1989.-Т. LX, № 10.- С. 2258-2262.

73. Сольчани, П. Моделирование многокомпонентной адсорбции без теплообмена / П. Сольчани, Р. Мохилла // Журн. прикладной химии.- 1986.-Т. LIX, №10.-С. 2226-2230.I

74. Устинова, Л.И. Динамика адсорбции углекислого газа на цеолите / Л.И. Устинова, В.Н. Лепилин, Т.И. Аладьина // Журн. прикладной химии.- 1980.- Т. LIII, №7.- С. 1489-1494.

75. Бобков, С.П. Применение вероятностных подходов для моделирования технологических процессов / С.П. Бобков // Изв. вузов. Химия и хим. технология.-2005.- Т. 48, №7.- С. 105-112.

76. Таганов, И.Н. Моделирование процессов массо- и энергопереноса / И.Н. Таганов. Л.: Химия, 1979. - 204 с.

77. Ходоров, Е.И. Разработка метода расчета аппаратов непрерывного действия для адсорбционной очистки жидкостных потоков / Е.И. Ходоров, В.Л. Ци-перман //Хим. промышленность.- 1999.- № 1.- С 48-53.

78. Лукин, В.Д. Автоматизированная система обработки экспериментальных данных по адсорбционному равновесию / В.Д. Лукин, А.В. Новосельский // Журн. прикладной химии.- 1986.-Т. LIX, № 6,- С. 1215-1240.

79. Глинский, В.А. Математическая модель сорбционного процесса в псев-доожиженном слое при неидеальном перемешивании твердой фазы / В.А. Глинский, И.О. Протодьяконов, Ю.Г. Чесноков // Журн. прикладной химии.- 1981.-- Т. LIV, №1.- С. 107-111.I

80. Лукин, В.Д. Информационный метод описания изотерм адсорбции на макропористых адсорбентах / В.Д. Лукин, А.В. Новосельский, В.А. Астахов // Журн. прикладной химии.- 1984.-Т. LVII, №5.- С. 993-996.

81. Устинов, Е.А. Моделирование процессов сорбционногофазделения газовых смесей в движущемся слое адсорбента / Е.А. Устинов// Журн. прикладной, химии.- 1995.-Т. 68, №5.- С 774-780.

82. Акулов, А.К. Моделирование циклического адсорбционного процесса с вакуум-термической регенерации / А.К. Акулов, Е.А. Устинов // Журн. прикладной химии.- 1981.- Т. LIV, №11.- С. 2437-2442.

83. Глинский, В.А. Математическая модель процесса адсорбции в неоднородном псевдоожиженном слое / В.А. Глинский, И.О: Протодьяконов, Ю.Г. Чес-ноков // Журн. прикладной химии.- 1980,- Т. LIII, №10.- С. 2370-2373.

84. Лукин, В.Д. О расчете динамики адсорбции на основе вариационных .принципов / В.Д. Лукин, А.В. Новосельский // Журн. прикладной химии.- 1981.Т. LIV, №11.- С. 2422-2426.

85. Моделирование сорбционных процессов разделения смесей / А.А. Се-балло, В.А. Шибаев, О.Л. Иванова и др. // Журн. прикладной химии.- 1988.-Т. LXI, №10.-С. 2244-2252.

86. Аксельруд, Г.А. Математическая модель изотермической адсорбции в условиях постоянной внешней концентрации адсорбтива / Г.А-. Аксельруд, В.Д. Степанов // Журн. физ. химии,- 1976.- № 8,- С. 2027-2031.

87. Татаренко, А.А. Моделирование процесса одноцентровой адсорбции на энергетически неоднородной поверхности / А.А. Татаренко, В.И. Быков // Изв. вузов. Химия и хим. технология.- 2000.- Т.43, № 6.- С.54-55.

88. Кисаров, В.М. Десорбция «-ксилола в неизотермических условиях / В.М. Кисаров, Р.Я. Фишер, Р.В. Корчагина //Хим. промышленность.- 1981.- № 1.- С. 22-23.

89. Математическая модель процесса регенерации в аппарате с неподвижным слоем ионита/JI.A. Зенькевич, В.А. Константинов, А.И. Волжинский, Н.Н. Смирнов //Журн. прикладной химии.- 1984.-Т. LVII, №10.- С. 2377-2380.

90. Устинов, Е.А. Моделирование десорбции из неподвижного слоя адсорбента нагретым потоком инертного газа / Е.А. Устинов // Журш прикладной химии.- 1981 .- Т. LIV, №9.- С. 1997-2002.I

91. Богословская, Н.С. Динамика десорбции^воды из цеолита нагретым потоком инертного газа / Н.С. Богословская, Е.А. Устинов, А.А. Себалло // Журн. прикладной химии.- 1987.-Т. LX, №12.- С. 2044-2047.

92. Устинов, Е.А. Моделирование динамики десорбции в движущихся слоях адсорбентов/ Е.А. Устинов, Р.А. Багиров, А.А. Себалло // Журн. прикладной химии.- 1982'.- Т. LV, №8. С. 1774-1779.

93. Акынбеков, Е.К. Расчет непрерывного процесса десорбции в аппарате с псевдоожиженным слоем / Е.К. Акынбеков // Изв. вузов. Химия и хим. технология.- 2001.- Т.44; №4.- С 97-99:

94. Математическая модель процессов десорбции с программированным нагревом слоя адсорбента / А.А. Себалло, В.М. Кисаров, В.А. Шибаев и др. // Хим. промышленность.- 1982.- № 12.- С. 746-748.

95. Шелыгин, A.JI. Кинетика терморадиционной десорбции толуола из активного угля АР-3 / A.JI. Шелыгин, Н.М. Корольков // Журн. прикладной химии.-1982.- Т. LV, №4,- С. 967-970.

96. Bunke, G. Breakthrough curves in the cyclic steady state for adsorption cysteines with concavt isotherms/ G. Bunke, D. Gelbin // Chemical Engineering Science.t1978.-Vol. 33, No 1.-P. 101-108.

97. Реутский, В.А. Расчет процессов изотермической десорбции / В.А. Ре-утский, Б.С. Сажин // Журн. прикладной химии.- 1988.-Т. LXI, №11.- С. 25222526.

98. Моделирование процессов вытеснительной десорбции/ А.А. Себалло,

99. B.А. Шибаев, O.JI. Иванова и др. // Журн. прикладной химии.- 1988.-Т. LXI, №7.1. C. 1617-1621.

100. Акулов, А.К. Динамика вытеснения кислорода азотом из слоя цеолита СаА / А.К. Акулов, А.К. Агабалян // Журн. прикладной химии.- 1988.-Т. LXI, №7.-С. 1621-1624.

101. Устинов, Е.А. Моделирование десорбции органических веществ из неподвижного слоя адсорбента водяным паром/ Е.А. Устинов, А.А. Себалло, Т.Г. Плаченов // Журн. прикладной химии.- 1982. Т. LV, №1.- С. 115-122.

102. Закономерности адиабатической десорбции органических веществ водяным паром в процессах рекуперации/ Е.А. Устинов, А.А. Себалло, В.М. Кисаров, Т.Г Плаченов // Журн. прикладной химии.- 1985.-Т. LVIII, №4.- С. 803-808.

103. Сипаров, С.В. Кинетика десорбции в бипористом зерне сорбента для случая нелинейного термического уравнения-сорбции / С.В. Сипаров, С.В. Бул-дырев, И.О. Протодьяконов^// Журн. прикладной химии.- 1983.- T.LVI, №12- С. 2665-2669.

104. Сипаров, С.В. Десорбция в бипористом* зерне сорбента'для случая нелинейного термического уравнения сорбции / С.В. Сипаров, И.О.1 Протодьяконов // Журн. прикладной химии.- 1983.- Т. LVI, №9.- С. 2005-2009.

105. Николаев, В.В. Регенерация цеолитных адсорбентов органическими растворителями /В.В. Николаев, А.Н. Вшивцев, А.Д. Шахов //Химия и технология топлив и масел.- 1997.- №6.- С. 27-28.

106. Медведев, С.Ф. Оценка задержки сорбента в колонных адсорберах, работающих в разных гидродинамических режимах / С.Ф. Медведев, А.В. Арефьев,I

107. С.М. Балакирев //Хим. промышленность.- 1988.- № 2.- С. 109-111.

108. Ходоров, Е.И. Разработка модели транспорта сорбента в секционированных адсорберах циклического режима / Е.И. Ходоров, Е.В. Тарабуткин //Хим. промышленность.- 1994.- № 2.- С. 128-132.

109. Интенсификация процесса массообмена в адсорбционной установке с подвижным зернистым слоем/ М.И. Сатаев, Ф.Е. Алтынбеков, М.Ш. Оспанов, К.К. Наурызбаев //Изв. вузов. Химия и хим. технология.- 2002.- Т. 45, №1.- С. 150151.• • . 145

110. Романов, Н.В; Исследование процесса разделения двухкомпонентных углеводородных смесей; на адсорбционных установках квазинепрерывного'дейст-вия / Н.В. Романов, Н.В; Лисицын, В.В. Кашмет // Журн. прикладной химии.-1988.- Т. LXI, №10.- С. 2255-2261. ; .

111. Андриевский, Б Л. Оптимизация разделения i компонентов врежиме фронтальной«адсорбции;/ Б.П. Андриевский, В*Б: Варганов, Юй: Титаренко' //Хим. промышленность.- 1990.- № 4.- С. 238-241.

112. Самойлов, Н. А. Особенности диффузии ароматических углеводородов при их сорбциище'олитами их растворов / Н.А. Самойлов, ILB. Заиченко // Нефтегазовые ^химические технологии: тез. докл. III Всерос. науч.-практ. конф./ СЕТУ.- Сама-, ра, 2005.- С. 124.

113. Самойлов, Hi А. К аналитическому расчету изотёрмшзбыточной и полной; адсорбции / Н:А. Самойлов, Н;В: Заиченко; Е.Б. .Тебенькова // Термодинамика поверхностных явлений и адсорбции: тез-, докл. Всерос: семинара / ИГХТУ.- Иваново, 2005.-С. 17-18.1. I .

114. Тебенькова, Е.Б. Исследование кинетики адсорбции изопропилбензола цеолитами NaX из раствора в н.гептане / Е.Б. Тебенькова, Н:В. Заиченко, Н.А.

115. Самойлов // 57-я научно-техническая конференция студентов, аспирантов и молодых ученых: сб. тез. докл./ УГНТУ.- Уфа, 2006.-Кн.2. С. 8.

116. Заиченко, Н.В. Кинетика адсорбции ароматических углеводородов из растворов в н.гептане цеолитами NaX / Н.В. Заиченко, Н.А. Самойлов // Изв. вузов. Нефть и газ. 2008.- №3. - С. 64-72.

117. Рабинович, В.А. Краткий химический справочник / В.А. Рабинович, З.Я. Хавин. Л.: Химия, 1991 - 432 с.

118. Влияние природы связующего на кинетику адсорбции диоксида углерода на цеолите NaA / Б.М. Эстрин, А.В: Нестеров, Е.А. Васильева и др. // Журн. прикладной химии.- 1991.-Т. 64, № 8.- С. 1764-1766.

119. Рид, Р. Свойства газов и жидкостей: определение и корреляция / Р. Рид, Т. Шервуд. Л.: Химия, 1971 - 702 с.

120. Самойлов, Н.А. Расчет кинетических кривых сорбции методом статистических моментов / Н.А. Самойлов // Журн. физ. химии.- 1976.- Т. 50, № Ю,- С. 2640-2641.

121. Заиченко, Н.В. Математическая модель динамики адсорбции / Н.В. Заиченко, Н.А. Самойлов / /Переработка углеводородного сырья. Комплексные решения: (Левинтерские чтения): тез. докл. Всерос. науч. конф. Самара, 2006.- С. 143-144.

122. Самойлов, Н.А. Моделирование динамики адсорбции цеолитами из жидкой фазы на основе данных статики адсорбции / НА. Самойлов, Н.В. Заиченко // Прикладная химия.- 2006.- Т. 79, № 7.- С. 1137-1143.

123. Закгейм, А.Ю. Введение в моделирование химико-технологических процессов: Математическое описание процессов / А.Ю. Закгейм. -М., Химия, , 1973 223 с.

124. Золотарев, П.П. Основы теории кинетики адсорбции бипористыми адсорбентами / П.П. Золотарев / /Адсорбция в микропорах: Труды пятой конференции по теоретическим вопросам адсорбции.- М., 1983. — С. 124-128.

125. Самойлов, Н.А. Перемещение и характер фронта адсорбции бензола цеолитами из смеси бензол — н-гептан / Н.А. Самойлов, Л.Ф. Фоминых // Доклады нефтехимической секции /Баш. респ. правление ВХО им. Д.И. Менделеева:-1971,-Вып. 6.-С. 174-178.

126. Кафаров, В.В. Методы кибернетики в химии и химической технологии / В.В. Кафаров. -М.: Химия , 1971 495 с.

127. Батунер, JI.M. Математические методы в химической технике / JI.M. Баjтунер, М.Е. Позин. Д.: Химия, 1971 - 823 с.

128. Самойлов, Н.А. Моделирование оптимальных условий реализации адсорбционного процесса / Н.А. Самойлов, Д.Н. Митягин, Н.В. Заиченко // Актуальные проблемы теории адсорбции, пористости и адсорбционной селективности: тез. докл.

129. XI Всерос. симпозиума с участием иностр. ученых.- Москва-Клязьма, 2007.- С. 44.

130. Митягин, Д.Н. Разработка метода оптимизации адсорбционной установки / Д.Н. Митягин, Н.В. Заиченко, Н.А. Самойлов // 58-я научно-техническая конференция студентов, аспирантов и молодых ученых: сб; тез. докл./ УГНТУ.-Уфа, 2008.- Кн.2.- С. 28

131. Самойлов, Н.А. Некоторые направления интенсификации адсорбционных процессов / Н.А. Самойлов // Журн. прикладной химии.- 1998.- Т. 71; №1. — С. 99-104.

132. Расчёт и многокритериальная оптимизация статических режимов мас-сообменных процессов на примере абсорбции в производстве газоразделения /

133. A.Ю. Торгашов, В.П. Кривошеев, Ю.К. Машунин, Ч.Д. Холланд // Изв. вузов. Нефть и газ. 2001.- №3. - С. 82-86.

134. Самойлов, Н.А. Формирование оптимального рабочего пространства вIадсорберах с неподвижным и движущимся слоем адсорбента / Н.А. Самойлов // Журн. прикладной химии. 2001.- Т. 74, № 10. - С. 1655-1663.

135. Самойлов, Н.А. Определение оптимальных размеров вертикальных адсорбционных колонн со стационарным слоем адсорбента / Н.А. Самойлов // Хим. промышленность.--!988. № 2. - С. 40-42.

136. Лукин, В.Д. Циклические адсорбционные процессы: (теория и расчет) /

137. B.Д. Лукин; А.В. Новосельский.- Л!: Химия , 1989 255 с.

138. Касаткин, А.Г. Основные процессы и аппараты химической технологии / А.Г. Касаткин. М.: Альянс, 2004. - 751 с.I

139. Самойлов, Н.А. Исследование адсорбционной жидкофазной рчистки бензина-растворителя цеолитами: Дис. . канд. техн. наук: 05.17.08/Н. А. Самойлов. Куйбышев, 1969. — 250 с.

140. Прайс-лист ГУП «Ишимбайский специализированного химического завода катализаторов». — http://www.ceolite.ru.

141. Павлов, К.Ф. Примеры и задачи по курсу процессов и аппаратов химиIческой технологии / К.Ф. Павлов, П.Г. Романков, А.А. Носков. — М.: ООО "Рус-МедиаКонсалт", 2004.-576 с.

142. Прайс-лист ОАО «Стальторг». http://www.staltorg-74.ru/pricerussia.

143. Прайс-лист ОАО «Уфимский НПЗ». http://www.riccom.ru/demand-sale 1814 order 2.htm.