автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.01, диссертация на тему:Метод оптимального проектирования систем простых ректификационных колонн с заданной топологией

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

00501°«-'*'

л*'

БОГУЛА НЭЛЛИ ЮРЬЕВНА

МЕТОД ОПТИМАЛЬНОГО ПРОЕКТИРОВАНИЯ СИСТЕМ ПРОСТЫХ РЕКТИФИКАЦИОННЫХ КОЛОНН С ЗАДАННОЙ ТОПОЛОГИЕЙ

05.13.01 - Системный анализ, управление и обработка информации (химическая технология)

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

з май жа

Казань-2012

005016643

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

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

Официальные оппоненты

Ведущая организация

доктор технических наук, профессор, Зиятдинов Надир Низамович

доктор технических наук, профессор, заведующий кафедрой ФГБОУ ВПО «Российский химико-технологический

университет им. Д.И. Менделеева» Комиссаров Юрий Алексеевич

доктор технических наук, профессор, заведующий кафедрой ФГБОУ ВПО «Казанский государственный энергетический университет»

Лаптев Анатолий Григорьевич

ФГБОУ ВПО «Санкт-Петербургский государственный технологический институт (технический университет)»

Защита состоится «24» мая 2012 г. в «14» часов на заседании диссертационного совета Д 212.080.13 при Казанском национальном исследовательском технологическом университете по адресу: 420015, г. Казань, ул. К. Маркса, 68, зал заседаний Ученого совета (А-303).

Отзывы на автореферат (в двух экземплярах), заверенные гербовой печатью учреждения, просим отправлять по адресу: 420015, г. Казань, ул. К.Маркса, д. 68, ФГБОУ ВПО «Казанский национальный исследовательский технологический университет», секретарю диссертационного совета Д 212.080.13.

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

Автореферат разослан «24» апреля 2012 г.

Ученый секретарь диссертационного совета Д 212.080.13 д.т.н., профессор

А.В. Клипов

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

Задачи исследования

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

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

• Разработка алгоритма проектирования оптимальной системы простых ректификационных колонн с заданной топологией на основе предложенного подхода;

• Исследование эффективности -предложенного.метода на примере оптимального проектирования отдельных простых ректификационных колонн и их систем.

Научная новизна работы

• Предложен"подхбд к опредёлению нижних оценок критерия оптимальности для отстемы Простых ректификационных колонн, позволяющий перейти от дискретных поисковых переменных — числа тареЛок в колоннах, к непрерывным структурным параметрам. В основе подхода лежит модификация уравнения; связывающего рабочую и равновесные концентраций компонента :чёрез локальный эффективный коэффициент полезного действия тарелки и структурные параметры;

• Предложен подход для получения верхней оценки критерия оптимальности, в': котором для расчета числа тарелок в колоннах используются округленные значения суммы структурных параметров;

• Предложен подход для реализации процедуре ветвления;

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

Практическая значимость

• На основе предложенного метода разработан алгоритм оптимального проектирования систем простых ректификационных колонн с заданной топологией;

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

• Результаты работы переданы для внедрения в исследовательскую и проектную деятельность в инжиниринговую фирму "Инжехим";

• Разработанный алгоритм внедрен в учебный процесс кафедры математического моделирования и оптимизации химико-технологических процессов Санкт-Петербургского технологического института (технического университета) для использования в лабораторном практикуме дисциплины «Системный анализ и математическое моделирование химико-технологических процессов»..

Апробация работы :

Основные положения, и .результаты диссертационной работы докладывались и обсуждались на 21 ¡ 22, 23, 24 Международных научных конференциях «Математические методы ,в технике и технологиях» (Саратов, 2008; Псков, 2009; Саратов, 2010; Киев, 2011), Второй Между народной научно-практической конференции «Компьютерное моделирование в химической

технологии и устойчивое развитие» (Киев, 2010), Научных сессиях КНИТУ (Казань, 2009-2012г:д\). . ;

Публикации

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

Структура и объем работы - Диссертация общим объемом 150 страниц, состоит из введения, 4 глав основного текста, выводов, списка "использованной литературы из 129 наименований и приложений. Работа содержит 8 рисунков и 40 таблиц.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

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

Приведена потарелочная математическая модель многокомпонентной ректификации и обзор методов её расчета. Отмечено, что наиболее эффективным методом расчета многокомпонентной ректификации является метод «Inside-out».

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

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

укрепляющей части (s = 1) и ш/ тарелок исчерпывающей части (s = 2), где j - номер колонны. Тарелки нумеруются сверху вниз отдельно для каждой части колонны. Задача выбора оптимального числа тарелок, т.е. определения

значений inj, s = l, .v = 2 , в секциях каждой ректификационной колонны с одновременной оптимизацией по непрерывным переменным в системе сводится к решению следующей задами:

/= min "','»/)> О

X1,1t1 ,Hl' y-j

fl/(xVX) = 0. ^ = 1,2, j = 1,...,N , 1 <mJs<mi-"m, (2)

y/J(x' ,u')<0, j = \,...,N, (3)

FJ• =0, F'"-W8 =0, j,g,r = \,..,N, (4)

где

<PJ,o(xJ, и j),i = l, <pJsl(xJ,FJ),s = 2,

(fll(x',uJ, m

<p'sk (xJ), 1 <k< mJs, где s = 1,1 < к < m's, где s = 2, ^

<р'. (х1 > " 11 ,)• = 2,

j,r,g - номера ректификационных колонн; х1,и1 - переменные состояния и управляющие переменные в у-ой ректификационной колонне; номер укрепляющей (5 = I) и исчерпывающей (.у = 2) секций колонн; /я/ - число тарелок в укрепляющей или исчерпывающей секциях У-ой ректификационной колонны, которые в общем случае могут принимать любые целые значения в пределах от I до ш/1"™ соответственно; ;п/тах - максимальное число тарелок в секциях колонн, это значение задается на начальной итерации пользователем; /'(х1 ,и',/»/) - суммарные капитальные и эксплуатационные затраты у-ой ректификационной колонны; уравнения (2, 5) - математические модели укрепляющей и исчерпывающей секций 7-ой ректификационной колонны;

<р{0(х'\ и ¿) - математическая модель дефлегматора в у -ой колонне; <р'1к{х') - математическая модель к -ой тарелки в у -ой колонне; ) - матема-

тическая модель тарелки питания в у-ой колонне; Дх', и матема-

тическая модель кипятильника в у -ой колонне; ограничения типа неравенств (3) - проектные ограничения; уравнения (4) - соотношения, характеризующие структуру системы ректификационных колонн, которые означают, что р'< = 0 - потоком питания у-ой колонны является дистиллят я -ой колонны; р'' - IV =0 - потоком питания г-ой колонны является кубовый продукт g -ой колонны; N - число колонн в рассматриваемой системе.

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

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

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

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

В процессе оптимального проектирования СРК для расчета материально-теплового баланса системы целесообразно использование универсальных моделирующих программ. Наиболее подходящей универсальной моделирующей программой для решения задачи оптимального проектирования СРК является программа ишзнп.

В главе сформулированы цель и задачи исследования.

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

Рассмотрим математическую модель к -ой тарелки ректификационной колонны (рисунок 1).

" ▼ ц

Рисунок 1 - Тарелка ректификационной колонны

На данной блок-схеме и при дальнейшем описании модели тарелки ректификационной колонны индекс к ( к = \,..,ЫТ ) будет обозначать номер тарелки в колонне (нумерация сверху вниз), индекс /' - номер компонента смеси (г = ],..,N0), где ИТ - число тарелок в колонне, N0 - число компонентов в смеси. Массовые доли / -го компонента на к -ой тарелке в жидкой и газовой фазе обозначим через х-к и у1к соответственно. Тарелки в колонне связаны между собой потоками пара Ук и жидкости Ьк. Также имеется возможность подавать на тарелку поток питания состава . Энтальпия каждого потока (Ну - парового, Н,~ жидкостного) рассчитывается на основе температуры на тарелках Тк. Давление Рк определяется на каждой тарелке колонны. •

Математическая модель к -ой тарелки включает следующие уравнения:

Уравнение покомпонентного материального баланса:

Р^л + + УмУц-< - АЛ - УкУ* = 0' (б)

Уравнения фазового равновесия:

У'к = К*хн. (7)

Кл=К(Т0Р,,ха). (8)

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

Ул = У,л« + ЕАУл -Уим)-

Уравнения нормировки составов жидкой и паровой фаз:

= = (ю)

Уравнения теплового баланса: ¡РН'п + - ЦНЛ - УкНл = 0, к = , (11)

Нш=Н{Тк,Рк,хк), Нук = Н(Тк,Рк,ук), Нп=Н(ТгРг,2,к).

В уравнении (9) Е-к является оценкой разделительной способности контактного устройства и характеризует движущую силу процесса, определяемую кинетикой массопередачи и гидродинамической структурой взаимодействующих потоков пара и жидкости.

Рассмотрим уравнение (9). Введем переменные а'к, где каждая а'к соответствует к -ой тарелке, -ой части колонны. С помощью новых переменных модифицируем уравнение (9):

Уь ~ У*/' + (Уь ~ Уь*') • (12)

Отметим, при а\к = 0 мы получаем, что у^ = у./1, .х* = х*~', то есть при а'л =0 к -я тарелка отсутствует. И, наоборот, при а[к =1 тарелка присутствует. Используя новые переменные, мы можем сформулировать задачу оптимального проектирования как задачу поиска оптимальных значений переменных а1к и режимных переменных. Запишем задачу поиска оптимального числа тарелок в следующем виде:

/= чип/(*',«',«/) (13)

.г' .»'.а/

<р\х\и>,а1) = 0, 5 = 1,2, у = 1,..., N (14)

(И(лг'У)20, _/ = !,..„/V (15)

«/=[0,1],

где от/ - вектор, компонентами которого являются а'л , к = 1,...,т"ш.

Преобразование задачи (1) в задачу (13) состоит в том, что новые переменные могут принимать непрерывные значения, что позволит получать нижние оценки. В этом состоит преимущество перехода от переменных тг к переменным а'к.

Введем множества М'/ номеров тарелок в укрепляющей и исчерпывающей секциях колонны (я = 1,2), где / - номер шага метода ветвей и границ. Эти множества определяются следующим образом: М/'' содержит номера тарелок, для которых а'л принимают фиксированные значения 0 или 1. Набор =[0,1] при кеМ/'' обозначим как а'к. Эти значения определены на предыдущих итерациях метода ветвей и границ.

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

/' = пип £/'(*>'.<*/).кем', (16)

л- = 1,2, У = к = 1...,т>/ш* (17)

= (18) , [фиксированные (известные) значения О или 1,если к е Ми, [переменные значения О или 1,если к <£ М{~'.

Нижняя оценка оптимального значения критерия задачи (16) будет определяться решением следующей задачи:

/''= т (19)

<р> (х>, и' ,а'л) = 0,) = 1,..., N (20)

у\х\и>,а1к)< 0,у = 1.....Л^ (21)

- 1 ■ Д™ всех к £ М1/, принимают фиксированные (известные) значения 0 или 1 для всех к е А//-'.

Задача (19) имеет только непрерывные поисковые переменные и для её решения может быть использован один из методов нелинейного программирования.

Верхнюю оценку будем вычислять на основе значений параметров «*/, полученных при решении задачи (19) и а'л , полученных на предыду-

щей итерации. Введем величину р'

I _

2>;/ кчи;'

+ а'л , где р1 - сово-

купность ближайшего целого суммы структурных параметров, полученных в результате решения задачи (19), и суммы структурных параметров, найденных на предыдущей итерации. Положим для к - 1,...р/ = 1, для всех остальных а$к = 0. Верхняя оценка получается решением следующей задачи: >]' = гпш/' {х> ,и\а]к) (22)

X ,11

= 0, j = \,...,N, а*А =1, А = 1,...,р/,5 = 1,2 (23)

0. (24)

Рассмотрим принцип ветвления некоторой вершины. На каждом шаге проводится разбиение одной секции колонн, причем 5-ая секция колонны делится на 2 множества. В первом множестве «-ой секции колонны параметры а'1к в интервале ([т/-'/2] + 1 > '»!') варьируются, в интервале (1,

[»(/ '/2]) а'л равны I; где /и/ ' - число тарелок, для которых к z М'/. Во втором множестве - параметры а'к в интервале (1, J m'5J¡2^) варьируются, в

интервале (["^'/2] +1, ) равны 0.

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

Шаг 1. Начальный шаг / = 1. Задание начальных приближений поисковых переменных х, и, а'л . Формирование множеств с номерами висячих вершин Q = 0 . Задание точности решения задачи е. Вычисление верхней

//д и нижней оценок для нулевой вершины А°, <7 = 1. Выбор наименьшей верхней оценки Р° =if0. Построение потомков А],А12 вершины . Задание границ изменения параметров a'sl.

Шаг 2. Вычисление нижних //', fi[ и верхних ^оценок для вершин-потомков а1,А[ вершины Л'г ' решением задач (19) и (22), соответственно. Внесение их номеров в множество Q .

Шаг 3. Выбор наименьшей верхней оценки Р1 = min{P''\q\,q[}.

Шаг 4. Отброс неперспективных вершин, удовлетворяющих правилу //;>/>', {q,p}zQ.

Исключение номеров неперспективных вершин из множества Q .

Шаг 5. Поиск перспективной вершины, удовлетворяющей правилу //,f.' = min //''.

Шаг 6. Проверка окончания решения для перспективной вершины Afi

Если условие выполнено, то решение найдено. Перейти на шаг 8.

Если условие не выполнено, то q = q' , перейти на шаг 7.

Шаг 7. Ветвление перспективной вершины A1 = 1 + 1. Переход на

шаг 2.

Шаг 8. Решение получено. Значение критерия оптимальности , режимы работы колонн и число тарелок в секциях колонн берется из задачи вычисления верхней оценки для вершины А^.

В завершение главы приводится структура многоуровневого алгоритма оптимального проектирования ректификационных колонн. Дается описание каждого уровня, включающее описание методов последовательного квадратичного программирования и Inside-Out (рис. 2).

(

Рисунок 2 - Структура многоуровневого алгоритма В третьей главе дается обзор основных подходов к расчету критерия оптимальности в задаче оптимального проектирования СРК. На их основе сформулирован критерий оптимизации суммарных приведенных капитальных и эксплуатационных затрат:

1 N ( г »•! ^ ы "

1ХТ «Ь+МЬ+МЬ

где а'1к - структурный параметр; ртс1 - цена металла, руб/кг; / - срок окупаемости проектируемой системы, ч; М'кЬ - масса метгшла (кг), расходуемого на изготовление кипятильника у-ой колонны; М'пт1 - масса металла (кг), расходуемого на изготовление конденсатора у -ой колонны; М[ - масса металла (кг), расходуемого на изготовление одной тарелки и секции кожуха у-ой колонны высотой, равной межтарельчатому расстоянию, р„ - цена хлада-

гента, руб/м3; 6',{ - расход хладагента в ./-ой колонне, м3/ч; р, - цена теплоносителя, руб/Гкал; (2^ - тепловая нагрузка на куб у-ой колонны, Гкал/ч. Величина М, рассчитывалась по формуле:

м, = (*,„, +лл21ту /4) ■ 2мРпк,,

где И - межтарельчатое расстояние, м; <11тг -диаметртарелки, м; г11к1 толщина листа металла, м; ртс1 - плотность металла, кг/м3; п = 3.14. Расход хладагента в у -ой колонне рассчитывался по формуле:

Q!

■coud

-Рп-

где <211и1 - тепловая нагрузка на конденсатор, Ккал/ч; Д/ - заданная разность температур на концах холодного теплоносителя, °С; с!>я, - теплоемкость хладагента, Ккал/кг-°С; р№ - плотность хладагента, кг/м3. '

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

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

Колонна дебутапизации (рисунок 3)

Исходные данные: расход сырья - 10000 кг/ч; состав сырья (массовые доли): н-бутан - 0.35, н-пентан - 0.3, н-гексан - 0.35; температура сырья-80°С; давление сырья - 5 бар;

Рисунок 3 - Колонна дебутапизации трехкомпонентной смеси Давление верха колонны было выбрано с учетом возможности конденсации верхних продуктов доступным хладагентом и принято равным 4 бар. Эффективность тарелок для ректификационной колонны была принята: равной 0.8. Проектные требования наложены на содержание н-бутана в бута но-

вой

фракции: л£„ > 0.98 массовых долей. Точность решения е = |(?7'была принята равной 0.05.

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

Ход решения задачи приведен в таблице 1 и показан на рисунке 4.

Таблица 1 - Ход решения задачи

№ вершины Нижняя оценка Верхняя оценка точность

т\ ,п\ Флегмовое число Критерий

01 1485.7 8 7 4 1875.6 0.208

11 1111. 11 7 3,16 1214 0,008

12 2254 4 7 6 2617 0.138

21 979.7 10 10 1.95 1010.7 0.03

12 ■ 1916.2 8 2 4.5 1954.5 0.019

31 959.3 13 9 2.1 1078.9 0.11

32 1198.2 8 8 2.5 1663.6 0.28

41 920.2 12 12 1.84 961.7 0.043

42 962.2 12 9 2 1015.4 0.052

51 959.7 14 12 1.8 960.8 0.001

52 934.6 12 12 1.8 944 0.01

61 931.8 12 14 1.7 932.6 0.0008

62 934.3 12 12 1.8 934.6 0.0003

71 938.8 13 14 1.8 978 0.04

72. 933.1 12 14 1.8 955.9 0.0238

На рисунке 4 приведен дерево-граф решения задачи. Оптимальной является вершина 61. Найденное флегмовое число - 1.7, температура куба 100.5 °С, общее число тарелок в колонне 26, тарелка питания - 13. Материальный баланс колонны приведен в таблице 2.

Параметр потока Сырье с4 С5, С6

Температура. °С 80 42.6 100.5

Давление, бар 5 4 4.5

Массовый расход, кг/ч 10000 3365 6635

Масс, доля н-бутан 0.35 0.9929 0.0239

Масс, доля-н-пентан 0.3 0.0071 0.4485

Масс, доля н-гексан 0.35 0.0000 0.5275

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

938.6 978

Рисунок 4 - Дерево-граф решения задачи (в прямоугольниках верхнее число -значение нижней оценки, нижнее число - верхней оценки)

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

Qcoml öro/ir/

Q^zQZ ■ (26)

на качество целевого продукта

c,nh

> 0.98,

(27)

(28)

на коэффициент захлебывания в колонне к,w^O.82.

В ограничениях (25), (26) предельные тепловые нагрузки

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

В число неопределенных параметров Qi вошли массовые доли н-бутаиа (/ = 1), н-пентана (;' = 2), н-гексана ( /=3). Эти параметры определили размерность п0 - 3 области неопределенности

Т(5) = Щ : 0,"в." <в? +<f J = 1,....яв}, где = 0.35, в? = 0.3, 0," =0.35. Для оценки размера области работоспособности был вычислен индекс гибкости F(d) спроектированной колонны

l-'(d) = m^(ö:+S!'), (29)

max min max ш. (d, и, в) < 0,

eeT(t) и ' '

где S = ,i = l,...,3} — массив параметров, задающий размер области не-

определенности Т; у = ..,4 - номера ограничений ^(с/,м,0), представленных неравенствами (25) - (28); и - режим работы колонны, определяется флегмовым числом и температурой в кипятильнике колонны.

Решение задачи (29) S' = {¿¡"',i = 1,...,3} даёт максимальную область неопределенности Т, обладающую следующим свойством: при любом в е Т система колонн гарантирует выполнение всех требований по качеству разделения и условиям физической реализуемости.

В процессе вычисления индекса гибкости F(d) колонны учитывалось

требование = 1, отвечающее условиям нормировки состава сырья.

В результате вычисления индекса гибкости спроектированной колонны были найдены значения допустимых границ диапазонов изменения массовых долей для каждого из компонентов сырья, которые приведены в таблице 3. Область работоспособности формируется из любых значений неопределенных параметров из найденных диапазонов, при условии, что соблюдаются требования +02 £ 0.35 и на нормировку состава сырья 0t=\.

Таблица 3 - Допустимые границы диапазонов изменения массовых долей

' компонентов сырья

Неопределенный параметр e,N-s <9/v +S?

Массовая доля н-бутана 0.11 0.87

Массовая доля н-пентана 0.05 0.75

Массовая доля н-гексана 0.05 0.65

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

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

Параметры сырья: расход1 сырья, поступающего на разделение -100000 кг/ч; состав сырья (массовые доли): пропан - 0.2632 массовых долей, н-бутан - 0.3158 массовых долей, н-пентан - 0.2105 массовых долей, н-гексан - 0.2105 массовых долей; температура потока сырья- 105.3°С; давление потока сырья- 16 бар.

.Сырье

СЗ.

К1

£

К2

С4

С5_

КЗ

С6

Рисунок 5 - Система колонн для разделения четырехкомпонентной смеси

Требования на качество продуктов: содержание пропана в пропановой фракции хсл >0.98 массовых долей, содержание бутана в бутановой фракции г,- „ > 0.98 массовых долей, содержание пентана в пентановой фракции

хг,, >0.98 массовых долей.

Давление верха колонн было выбрано с учетом возможности конденсации верхних продуктов доступным хладагентом. Исходя из этого, оно было принято равным: для колонны К1 - 14 бар, для колонны К2 - 4 бар, для колонны КЗ - 1.5 бар.

Точность решения е была принята равной 0.05. Начальное приближение числа тарелок для каждой колонны было принято равным 20. В качестве поисковых переменных были приняты расходы флегмы в колонны, темпера-; туры кубов колонн, переменные а[к. Таким образом, общая размерность задачи вычисления нижней оценки была равна 64.

• Ход решения задачи приведен в таблице 4. Решение было получено'за 5 шагов. Значение относительной невязки между нижней и верхней оценкой критерия в этой вершине составляет 0.02, что удовлетворяет заданной точности 0.05. представленные следующие за вершиной 3-1 результаты ветвления подтверждают, что найденная вершина действительно является оптимальной.

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

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

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

СОСТОЯЩеЙ ИЗ 5-ТИ КОЛОНН. .■:..; ,

Таблица 4 - Ход решения задачи проектирования системы из трех _ректификационных колонн __

Верхняя оценка

№ Нижняя Колонна К1 Колонна К2 Колонна КЗ

вершины оценка т\ '"Г т\ 3 '»1 3 Критерий Точность

1-1 230685 7 5 6 5 5 5 249889 0.077

1-2 185170 3 5 5 5 5 5 255734 0.275

2-1 253540 7 5 7 5 5 5 300858 0.157

2-2 232573 7 5 3 5 5 5 273220 0.148

3-1 236429 7 5 3 5 7 5 241273 0.02

3-2 213925 7 5 3 5 3 5 252074 0.151

4-1 270606 7 7 3 5 7 5 281551 0.038

4-2 225381 7 3 3 5 7 5 313015 0.279

Таблица 5 - Значения приведенных капитальных и эксплуатационных затрат в вершинах графа ветвления

№ вершины Верхняя оценка критерия Приведенные капитальные затраты Эксплуатационные затраты

1-1 249889 142984 106905

1-2 255734 120209 135525

2-1 300858 172243 128435

2-2 273220 144442 128778

3-1 241273 134136 107137

3-2 232074 108333 123741

4-1 281551 168525 113026

4-2 313015 155880 157135

Таблица 6 - Оптимальные значения поисковых переменных

Параметр Колонна К1 Коломна К2 Колонна КЗ

Число тарелок в укрепляющей части 7 3 7

Число тарелок в исчерпывающей части 5 5 5

Флегмовое число 17.5 9.6 9.4

Температура куба, °С 125.4 105.8 99.5

Таблица 7 - Материальный баланс спроектированной системы

Параметры потока / поток Сырье с,н8 1 С4Н10 2 С,Н,2 С6НИ

Температура, °С 105.3 42.2 125.4 42.9 105.8 47.7 99.5

Давление, бар 16 14 15 4 5 1.5 2.5

Расход, кг/ч 100000 26780 73220 31130 42090 20530 21560

Состав, масс, доли

Пропан 0.2632 0.981 0.001 0.0016 0.0000 0.0000 0.0000

Бутан 0.3158 0.0192 0.4242 0.9859 0.0088 0.0181 0.0000

Пентан 0.2105 0.0000 0.2875 0.0124 0.4910 0.9817 0.0238

Гексан 0.2105 0.0000 0.2875 0.0000 0.5002 0.0003 0.9762

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

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

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

3. Разработаны подход и алгоритм вычисления нижней оценки критерия оптимальности для решения задачи ДННП, основанный па модификации уравнения', связывающего рабочую и равновесные концентрации компонента через локальный эффективный коэффициент полезного действия, с помощью введения структурных параметров. Такой подход позволяет переходить от дискретных структурных параметров к непрерывным, в результате чего задача йи>кней оценки критерия оптимальности сводится к задаче нелинейного программирования.

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

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

6. На основе предложенных подходов разработан метод оптимального проектирования СРК с заданной топологией.

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

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

Публикации в ведущих научных рецензируемых журналах, рекомендованных ВАК:

1. Островский r.ht. Оптимальное проектирование системы ректификационных колонн. / Г.М. Островский, H.H. Зиятдинов, Т.В. Лаптева, H.IO. Вогула//Доклады Академии Наук, 2010. Т. 431. № 6. С. 768-771.

2. Островский Г. М. Оптимальное проектирование системы ректификационных колони с заданной топологией. / Г.М. Островский, Н. Н. Зиятдинов, Т. В. Лаптева, Н. Ю. Богула // Теоретические основы химической технологии, 2011. Т. 45. №1. С. 88-97.

3. Зиятдииов Н.Н. Метод оптимального проектирования ректификационной колонны. / Н.Н. Зиятдииов. Н.Ю. Богула, Т.В. Лаптева, Г.М. Островский // Вестник Казанского технологического университета, 2011. №5, С. 118-123.

Прочие публикации по теме научного исследования:

4. Зиятдииов Н.Н. О подходе к решению задачи оптимального проектирования системы ректификационных колонн методом ветвей и границ. / Н.Н. Зиятдииов. Н.Ю. Богула, Г.М. Островский // Восточно-Европейский журнал передовых технологий. - Киев, 2010. №2/10 (44). С. 13-16.

5. Богула Н.Ю. Метод INS1DE-OUT для расчета однократного испарения / Н.Ю. Богула, Н.Н. Зиятдииов, Т.В. Лаптева, Г.М. Островский, Д.А. Рыжов // Математические методы в технике и технологиях. Сб. трудов 21-й Международ. науч. конф. - Саратов, 2008. - Т. 5. - С. 187-188.

6. Рыжов Д. А. Моделирование и оптимизация режимов работы узла разделения изоамилен-изопреновой фракции производства изопрена / Д.А. Рыжов, Н.Н. Зиятдинов, Т.В. Лаптева, В.А. Курбатов, Н.Ю. Богула // Математические методы в технике и технологиях. Сб. трудов 22-й Международ, науч. конф. -Псков, 2009. - Т. 10. - С. 79-81.

7. Богула Н.Ю. О подходе к решению задачи оптимального проектирования системы ректификационных колонн методом ветвей и границ / Н.Ю. Богула, Н.Н. Зиятдинов // Компьютерное моделирование в химической технологии и устойчивое развитие. Тезисы докладов второй межд. научно-практич. конф. -Киев: НТУУ «КПП», 2010.-С. 81-82.

8. Богула Н.Ю. Оптимальное проектирование ректификационной установки / Н.Ю. Богула, Н.Н. Зиятдинов // Математические методы в технике и технологиях. Сб. трудов 23-й Международ, науч. конф. - Саратов, 2010. - Т.2. -С.75-78.

9. Богула Н.Ю. Программный комплекс оптимального проектирования системы ректификационных колонн / Н.Ю. Богула // Сборник докладов участников программы УМНИК. - Белгород, 2010. - С. 14-17.

10. Богула Н.Ю. Метод оптимального проектирования систем разделения / Н.Ю. Богула, Н.Н. Зиятдинов, Д.А. Рыжов, Г.М. Островский // Математические методы в технике и технологиях. Сб. трудов 24-й Международ, науч. конф. - Саратов, 2011. - Т.2 - С. 36-39.

Заказ Ж,____.____Тираж/<7(7 экз.

Офсетная, лаборатория КНИТУ, 420015, Казань, К.Маркса, 63

61 12-5/3617 - /

КАЗАНСКИЙ НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ ТЕХНОЛОГИЧЕСКИЙ УНИВЕРСИТЕТ

МЕТОД ОПТИМАЛЬНОГО ПРОЕКТИРОВАНИЯ СИСТЕМ ПРОСТЫХ РЕКТИФИКАЦИОННЫХ КОЛОНН С ЗАДАННОЙ ТОПОЛОГИЕЙ

05.13.01 - Системный анализ, управление и обработка информации

(в химической технологии)

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

БОГУЛА НЭЛЛИ ЮРЬЕВНА

Диссертация на соискание ученой степени кандидата технических наук

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

Казань - 2012 г.

СОДЕРЖАНИЕ

Введение...........................................................................................5

Глава 1. Обзор методов математического моделирования и оптимального проектирования ректификационных колонн и систем...........................13

1.1. У прощенные методы проектирования ректификационных колонн...........15

1.2.Строгие методы оптимального проектирования ректификационных установок.................................................................................17

1.2.1. Математическая модель ректификационной колонны.....................17

1.2.2. Методы расчета фазового равновесия........................................21

1.2.3. Методы расчета эффективности тарелок.....................................25

1.2.4. Методы расчета ректификационных колонн.................................27

1.3. Формализованная постановка задачи оптимального проектирования систем ректификационных колонн.................................................30

1.4.Подходы к решению задачи оптимального проектирования ректификационных колонн и систем...............................................33

1.4.1. Поиск оптимальной тарелки питания при постоянном числе тарелок........33

1.4.2. Поиск оптимального числа тарелок, оптимальной тарелки питания и режимов работы колонны........................................................35

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

1.5.1. Методы нелинейного программирования...................................37

1.5.2.Методы дискретно-непрерывного программирования.....................44

1.6. Решение задач исследования гибкости проектируемого технологического объекта..........................................................................................47

1.7. Возможности универсальных моделирующих программ для оптимального проектирования системы ректификационных колонн.................................50

Выводы........................................................................................55

Постановка задачи исследования........................................................56

Глава 2. Разработка многоуровневого алгоритма оптимального проектирования систем ректификационных колонн............................57

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

2.1.1 Подход к вычислению нижней оценки критерия оптимальности...........................................................................57

2.1.2 Подход к вычислению верхней оценки критерия оптимальности.....60

2.1.3 Процедура ветвления............................................................60

2.1.4 Алгоритм метода ветвей и границ оптимального проектирования систем ректификационных колонн..................................................61

2.2 Метод последовательного квадратичного программирования ..............67

2.3 Двухциклический метод расчета колонн Inside-Out.............................68

2.4 Многоуровневый алгоритм оптимального проектирования СРК с

использованием УМП Unisim............................................................70

Выводы.......................................................................................72

Глава 3. Решение задач оптимального проектирования простых ректификационных колонн..............................................................74

3.1 Построение критерия оптимальности проектируемой системы ректификационных колонн................................................................74

3.2 Оптимальное проектирование колонны дебутанизации........................77

3.2.1. Постановка задачи оптимального проектирования колонны дебутанизации...........................................................................77

3.2.2. Ход решения задачи.............................................................79

3.2.3 Результаты решения для различных вариантов состава сырья.......................................................................................81

3.2.4 Оптимальное проектирование колонны дебутанизации для неравновесных тарелок.................................................................82

3.3 Анализ гибкости оптимального проекта колонны дебутанизации...........84

3.3.1 Алгоритм вычисления индекса гибкости....................................84

3.3.2 Результаты решения задачи анализа гибкости.............................92

3.4. Задача оптимального проектирования колонны разделения пропан-

пропил еновой фракции..................................................................93

Выводы.......................................................................................97

Глава 4. Решение задач оптимального проектирования систем ректификационных колонн..............................................................98

4.1. Оптимальное проектирование системы из двух колонн для выделения пропановой и бутановой фракций из четырехкомпонентной смеси...........................................................................................98

4.2. Оптимальное проектирование системы из трех ректификационных колонн для разделения четырехкомпонентной смеси........................................102

4.3. Оптимальное проектирование газофракционирующей установки ГФУ-300..107

4.3.1. Постановка и решение задачи...............................................107

4.3.2. Оценка достоверности полученного решения............................117

Выводы......................................................................................120

Основные результаты и выводы......................................................121

Список использованных источников................................................123

Приложение 1. Вычислительный эксперимент. Решение задачи 3.2 для

разных составов сырья....................................................................136

Приложение 2. Вычислительный эксперимент. Решение задач

оптимального проектирования ректификационных колонн...................143

Справки и акты.............................................................................151

ВВЕДЕНИЕ

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

Особое внимание при этом уделяется вопросам оптимального проектирования процессов ректификации [1], поскольку системы ректификационных колонн для разделения многокомпонентных смесей являются высоко металло- и энергоемкими установками [1-4].

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

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

Цели работы

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

Задачи исследования

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

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

• Разработка алгоритма проектирования оптимальной системы простых ректификационных колонн с заданной топологией на основе предложенных подходов;

• Исследование эффективности предложенного метода на примере оптимального проектирования отдельных простых ректификационных колонн и их систем.

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

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

■ Предложен подход для получения верхних оценок критерия оптимальности, в котором для расчета числа тарелок в колоннах используются округленные значения суммы структурных параметров;

■ Предложен подход для реализации процедуры ветвления;

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

Практическая значимость

■ На основе предложенного метода разработан алгоритм оптимального проектирования систем простых ректификационных колонн с заданной топологией;

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

■ Результаты работы переданы для внедрения в исследовательскую и проектную деятельность в ООО «Инженерно-Внедренческий Центр «Инжехим»;

■ Разработанный алгоритм внедрен в учебный процесс кафедры математического моделирования и оптимизации химико-технологических процессов Санкт-Петербургского технологического института (технического университета) для использования в лабораторном практикуме дисциплины «Системный анализ и математическое моделирование химико-технологических процессов».

Апробация работы

Основные положения и результаты диссертационной работы докладывались и обсуждались на Международных научных конференциях «Математические методы в технике и технологиях»: ММТТ-21 (Саратов, 2008 г.), ММТТ-22 (Псков, 2009 г.), ММТТ-23 (Саратов, 2010 г.); ММТТ-23 Киев, 2011; Международной научно-практической конференции «Компьютерное моделирование в химической технологии и устойчивое развитие» (Киев, 2010).

Публикации

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

Структура и объем работы

Диссертация общим объемом 150 страниц, состоит из введения, 4 глав основного текста, выводов, списка использованной литературы из 129 наименований и приложений. Работа содержит 21 рисунок и 62 таблицы.

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

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

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

Приведена потарелочная математическая модель многокомпонентной ректификации и методы её расчета. Отмечено, что в настоящее время наиболее эффективным методом расчета многокомпонентной ректификации является метод «inside-out».

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

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

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

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

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

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

В главе сформулированы цель и задачи исследования.

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

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

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

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

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

1) колонны дебутанизации трехкомпонентной смеси: бутан, пентан, гексан. Проведена оценка области ее работоспособности при изменении условий эксплуатации;

2) пропан-пропиленовой колонны;

3) колонны для выделения воды из смеси моно- и диэтано л аминов;

4) колонны дебутанизации пятикомпонентной смеси.

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

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

2) система из трех простых ректификационных колонн для разделения четырехкомпонентной смеси;

3) системы г�