Лазерный оптико-акустический октанометр в системах оптимизации процесса компаундирования товарных бензинов Текст научной статьи по специальности «Математика»

УДК 665.7

ЛАЗЕРНЫЙ ОПТИКО-АКУСТИЧЕСКИЙ ОКТАНОМЕТР В СИСТЕМАХ ОПТИМИЗАЦИИ ПРОЦЕССА КОМПАУНДИРОВАНИЯ ТОВАРНЫХ БЕНЗИНОВ

© 2009 В. Н. Астапов Самарский государственный аэрокосмический университет

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

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

1. Постановка задачи оптимизации

Оптимальное решение задачи компаундирования бензинов является целью получения товарного бензина требуемого качества при минимизации его стоимости. Стоимость

товарного бензина Съ определяется стоимостью сг и количеством каждого из входящих в него компонентов: Съ сгиг . Поэтому цель задачи управления процессом смешения математически можно представить в следующем виде [1]:

СТи ® min, (1)

(2)

которым поступают компоненты смеси (предполагается, что все эти трубопроводы одинаковы).

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

| иі (к +1) - иі (к)|< Ь, і = 1, п,

(3)

где С = |с1,..., еп\ - вектор стоимости входящих в товарный бензин компонентов; и =| их,.., ип |Т - вектор величин потоков смешиваемых компонентов; п - число смешива-

где к - моменты дискретного времени к = 0,1, ...; запись и (к) обозначает зависимость величины и от конкретного момента времени; Ь - известная величина, задаваемая оператором-технологом.

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

емых компонентов п = 2,16.

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

На управляющие воздействия накладываются естественные ограничения, связанные с физическим смыслом этих величин:

qTU (к) = G,

(4)

где г = 1, ..., п, - компоненты вектора и;

d - известная величина, равная максимальной пропускной способности трубопроводов, по

где q - вектор размерности п, все компоненты которого равны единице; G - заранее заданная величина, которая устанавливается оператором-технологом.

Качество товарного бензина нормируется государственными стандартами или стандартами предприятия. Среди нормируемых показателей основным является октановое число QS товарного бензина. Требования,

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

(5) | 8 (к) |< і, "к > 0.

(9)

где - соответственно верхняя и ниж-

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

Известно, что чем ниже октановое число нефтепродуктов, тем меньше его стоимость. Поэтому желательно, чтобы октановое число получаемого товарного бензина соответствовало нижней допустимой границе интервала (5):

& = & .

(6)

Математическая модель технологического процесса смешивания бензиновых фракций представляет собой нелинейное уравнение [2]:

N 1 N -1 N /(

е=Ъе>,+-— її,. 1 а

М "=;+1

е 0 - б,0

■ї,,

-Х].ХГ,

I=1

(7)

где а - октановое число товарного бензина,

- октановые числа чистых компонентов,

Х], х, - мольная доля] -го (,-го) компонента в растворе.

Данное уравнение может быть описано следующим векторным уравнением:

( (к) = и (к )е*^ти (к )]-1 + б (к):

(8)

где ( =|(*,..,(П |т - вектор октановых чисел компонентов товарного бензина в смешиваемых потоках; q =|1,...,1|т - вектор размерности п, все компоненты которого равны единице; 8(к) = 8(и(к), Є) - величина, равная второму слагаемому правой части соотношения (7) и характеризующая степень нелинейности процесса; к = 0, 1, ... - моменты дискретного времени.

Практика компаундирования топлив свидетельствует, что величина 8, характеризующая нелинейность процесса, не превышает некоторой постоянной величины і :

Решение задачи (1) при ограничениях (2) - (6), (9) может быть получено путем вычисления вектора оптимальных управляющих воздействий методами линейного программирования, если известен вектор октановых чисел Qi компонентов смеси модели

(8) .

Однако октановые числа Qi смешиваемых компонентов точно неизвестны, а известны лишь их интервальные оценки следующего вида:

(10)

для всех г = 1,2, ... , п.

Поэтому задача оптимизации процесса компаундирования сводится к поочередному решению следующих задач:

- идентификации неизвестных параметров Q модели (8), например по [1];

- решению задачи линейного программирования (1) при ограничениях (2) - (6) и

(9).

2. Идентификация параметров математической модели

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

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

октановым числом с первого цикла смешения.

Используя генетический алгоритм (ГА), решим систему уравнений ЛОА-октаномет-ра в общем виде (в виде системы нелинейных интегральных уравнений).

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

В этом случае для определения концентраций компонент в К-компонентной бензиновой смеси необходимо решить систему из 2К линейных уравнений (11). Здесь неизвестными являются значения концентраций компонент бензиновой смеси и значения показателей неселективного поглощения.

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

пары каналов показатели Д. можно принять равными константе. В матричной форме система линейных алгебраических уравнений лазерного ОА-октанометра имеет вид [4]:

KN +Ь = у,

(11)

к1(їм ) к2(їм )••• кк (їм )

N - М/2-мерный вектор долей (концентраций) компонент бензиновой смеси (М/2=К):

N =

N,

N

м/2,

b - M-мерный вектор показателей фонового неселективного поглощения:

b

b (4)

Ь &) = b (Л)

b (їм-l) b (їм ) = b (їм-l)

y - М-мерный вектор правых частей системы уравнений (11):

ХЛ) " >1 "

У (ї2) У 2

У (їм-1) yM-1

АК) _ _ Ум _

где К - матрица коэффициентов поглощения компонент исследуемой бензиновой смеси размерностью М х М/2:

ВД) к2(ї1)... кк (Л)

|к1(ї2) к2 (ї2 ) •• • кк (ї2)

K =

Перейдем от системы (11) к системе разностных уравнений (вычитая попарно уравнения, для которых можно положить показатели b j = const). В этом случае вместо (11) имеем

(12)

где у - М/2-мерный вектор правых частей системы уравнений (12):

Уі 1 1у - 2у і

к11 к12 ••• к1к У 2 4 У - 3

к21 2 к 2 2 к у =

У(M-1)/2 ) _ Ум/2 _ ум-3 Ум-2

км 1 к к Лм2 ••• мк _ 5 _ ум-1 - ум _

К - матрица разностей коэффициентов поглощения компонент исследуемой бензиновой смеси размерностью М/2 х М/2:

К =

к11 к12 - к1к

к21 2 2 к2к =

кМ 1 к ЛМ 2 ... к Мк _

- к21 2 2 к - 2 • к1к к 2 к -

- к41 2 - 2 к . к3к - к4к

к - к

Л( М -1)1 ЛМ1

к - к

(М -1)2 ЛМ 2

к - к

Л( М -1)к Мк

1

2К

2 К

X(у-- у? )2>

где ? - индивид, для которого вычисляется целевая функция; у. - значение измеренного сигнала ву-м спектральном канале; у- - рассчитанное значение сигнала в у-м спектральном канале, полученное путем подстановки значения гена индивида ? в левую часть системы уравнений (11):

у?=Х +Ь

(15)

Для уравнения (12) квазирешение находится из условия

(13)

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

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

Рассмотрим особенности построения ГА для решения задачи (11) методом подбора квазирешения. Каждый индивид, кандидат на решение ? {т, g} (т - номер индивида в популяции, g - номер поколения) включает в себя все неизвестные параметры задачи - гены, разделенные на две хромосомы длиной К (рис. 1), одна из которых - вектор неизвестных концентраций компонент, а вторая - вектор неселективного поглощения. Хромосомы кодировались с плавающей точкой - это наиболее подходящее кодирование для данной задачи.

За целевую функцию (функцию пригодности индивидов) возьмем функционал невязки вида

Ы? - ген 1-й неизвестной концентрации индивида ?; Ь? ] - ген у-го неизвестного показателя фонового поглощения индивида ?.

Необходимо решить задачу минимизации функции F(S) вида (14).

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

Введем физические ясные ограничения:

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

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

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

0 < N < 1,1 = 1,2,..., К;

.=1

(16)

ь, > о, у = 1, 2,..., 2К.

(14)

Ареал в некоторых случаях можно сократить дополнительно. Например, значения концентраций (долей) компонент N. товарно-

Индивид S{m,g}

Рис.1. Неизвестные параметры задачи лазерного оптико-акустического октанометра

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

В начале работы ГА сформируем стартовую популяцию. Затем на каждой итерации к индивидам популяции будем применять генетические операторы отбора, репродукции, мутации и селекции.

Рассмотрим конкретную реализацию используемых в данном алгоритме генетических операторов.

Задание стартовой популяции. На начальном этапе работы алгоритма необходимо сформировать стартовую популяцию. Для этого случайным образом формируется М индивидов ?{т,0}, т.е. задается М вариантов векторов начальных концентраций компонент бензиновой смеси и векторов неселективного поглощения:

?{т,0} = (N1 т’0),

Д( т,0),

N (т ,0)

1У 2 ,

Р 2т,0),.

N (т ,0) г ’

р( т,0)

(т,0). К ;

Р Кт,0)),

(17)

где т = 1, ... , М.

Для полученных индивидов ?{т, 0} стартовой популяции вычисляется целевая функция (14). Стартовая популяция сортируется по возрастанию целевой функции F(S{m,0}) (наименьший индекс т соответствует наименьшему значению целевой функции) и используется как популяция предков на первой итерации работы алгоритма.

Отбор. Отбор индивидов из популяции предков для участия в репродукции осуществляется с использованием ранжирования по нормальному закону. Расчет вероятности выбора Р(?{т^}) индивида осуществляется по формуле

ЗХ, 1 < т < М -1

2 3 г ( -9х к М Л ехр I 2М

2

1х, т = М.

(18)

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

Репродукция. К каждой паре индивидов ?{т1,}, ?{т2^} популяции предков, прошедших оператор отбора, применяется оператор репродукции - модифицированный вариант равномерного кроссинговера и образовывается один новый индивид ? ’{I, g} популяции потомков. Заполнение хромосом нового индивида осуществляется следующим образом:

N{m(,^, апёе р1 = 0;

апёе р1

N1'" 2, g}, апёе р1 = 1;

Рг

}, апёе р 2 = 0;

р{т 2^ }(1 + w), апёе р2 = 1,

(19)

где 1=1, ., Мс - номер индивида популяции потомков; N{l•g} - значения генов концентраций индивида, подвергаемого мутации; N{tl'g} тШ - значения генов концентраций индивида после мутации; ^} - значения ге-

нов показателей неселективного поглощения

индивида, подвергаемого мутации; Д{

){1 ^} г тШ

значения генов показателей неселективного поглощения индивида после мутации; Р3,р4 - случайные величины, принимающие значения 0 и 1; н - случайная величина, распределенная по нормальному закону с математическим ожиданием, равным нулю, и дис-

~ “2 персией, равной а , причем значение дисперсии случайной величины н меняется (адаптируется) в процессе эволюции.

Для индивидов популяции потомков, прошедших оператор мутации, вычисляются значения целевой функции (14).

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

Таким образом, на каждой итерации работы генетического алгоритма сначала к популяции предков из М индивидов ?{т^} применяются операторы отбора и репродукции, в результате чего образуется популяция потомков из Мс индивидов ?’{1^}. Популяция потомков подвергается мутации. Далее в результате применения оператора селекции М лучших индивидов образуют популяцию предков следующего поколения ?{т^=1}.

Адаптация дисперсии мутации. Один раз в 0М поколений осуществляется измене-

2

ние значения дисперсии а случайной величины Ж в зависимости от относительной частоты / (/=Сйе/0МР где Оскс - число поколений, в которых произошло уменьшение целевой функции) уменьшения целевой функ-

ции F(S{0'g}) лучшего индивида ?{0^} поко-ления, посчитанной за эти Ом итераций. Причем, если частота уменьшения целевой фун-

2

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

а2 а, апёе / > /р; а2/ а, апёе / < /,

где а - новое значение дисперсии мутации;

2

а - предыдущее значение дисперсии мутации; а > 1,0 < /р < 1 - некоторые числа.

Типовые значения параметров а и /р, применимые для широкого класса задач, могут быть взяты из [4]:

/р @ 0,2; а @ 1,22.

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

F(?{0, ^) < О,.

(20)

Здесь ?{0, g} - лучший индивид поколения g; - порог остановки итераций алгорит-

ма.

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

Для устранения возможного зацикливания алгоритма или чрезмерно долгой его работы условие (20) может быть скомбинировано с условием окончания работы по достижению некоторого числа поколений:

О = О ,

таху

где Отах - поколение, в котором должна быть остановлена работа алгоритма.

По окончании итераций хромосомы лучшего индивида £{0, g} (с точки зрения наименьшего значения целевой функции) выбираются за компоненты концентраций (долей) бензиновых компонент N и показатели

неселективного поглощения искомого квазирешения.

3. Этапы работы ЛОА-октанометра в режиме определения долей компонентов в бензиновой смеси

Работа ЛОА-октанометра состоит из следующих этапов.

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

2. Измерение величины поглощения лазерного излучения анализируемой смесью.

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

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

На втором этапе анализа осуществляется измерение поглощения анализируемой бензиновой смеси. Программное обеспечение, управляющее процессом измерения, со-

стоит из двух частей. Одна часть входит в состав программного комплекса, который реализует алгоритм определения НСКИ и на основе выбранного НСКИ выдает команды на измерение для второй части программного обеспечения.

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

На третьем этапе обработки результатов измерений осуществляется определение долей компонентов анализируемой бензиновой смеси.

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

4. Калибровка ЛОА-октанометра

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

При калибровке ЛОА-октанометра определяется постоянная оптико-акустической ячейки (ОА-ячейки) К:

иоп

п?п

иII

иоїї и о

N* (к™ - к£ У

(21)

где Nзо - концентрация (доля) изооктана в калибровочной смеси (#шо=0,95); к(0’\, к/

- коэффициенты поглощения изооктана в оп-, /-спектральных каналах соответствен-

но; иоп, ио// - амплитуды акустического сигнала, измеренные в оп-, о//-спектральных

каналах соответственно; и™ , и°{ - амплитуды опорного сигнала, измеренные в оп-, о//-спектральных каналах соответственно.

Для калибровки ЛОА-октанометра ОА-ячейка заполняется смесью известного состава. Коэффициенты поглощения компонент, входящих в состав этой смеси, известны с высокой точностью. В качестве калибровочной смеси принято использовать 2-х компонентную смесь изооктана (95%) и н-гептана

(5%).

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

Поэтому к калибровке ЛОА-октаномет-ра предъявляются высокие требования, зак-

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

Библиографический список

1. Астапов В. Н., Бакан Г. М., Сальников Н. Н. Оценивание с помощью эллипсоидов параметров линейной регрессии при линейных ограничениях на вектор входных переменных //Автоматика. - 1993. - №1. - С.28-34.

2. Астапов В. Н., Бакан Г. М., Коцюба А. Т., Одинцова Е. А. Математическое моделирование технологического процесса смешивания бензиновых фракций //Автоматика.

- 1992. - №5. - С. 31 - 37.

3. Астапов В. Н. Физические основы оптико-акустического метода контроля качества углеводородных топлив // Научные труды VII международной научно-практической конференции «Фундаментальные и прикладные проблемы приборостроения, информатики и экономики». Кн. «Приборостроение». Московская государственная академия приборостроения и информатики. - М.: МГАПИ, 2005. - С. 19-23.

4. Schoeneburg E., Heinmann F., Feddersen S. Genetische Algorithmen und Evolutionsstrategien: Eine Einfuerung in Theorie und Praxis der simulierten Evolution. - Bonn; Paris; Readin; Mass. [u.a.]: Addison-Wesley, 1994.

References

1. Astapov V. N, Bakan G. M, Salnikov N. N. Estimation of linear regression parameters with the aid of ellipsoids with linear constraints on the vector of input variables. // Automatics. - 1993. - No. 1. pp. 28-34.

2. Astapov V. N, Bakan G. M, Kotsyu-ba A. T., Odintsova Ye. A. Mathematical modeling of petroleum fraction mixing process // Automatics. - 1992. - No. 5. - pp. 31-37.

3. Astapov V. N. Physical foundations of the optico-acoustical method of hydrocarbon fuel quality control // Transactions of the VIIth

international practical science conference «Fundamental and applied problems of instrument engineering, information science and economics». “Priborostroyeniye”. Moscow State Academy of Instrument Engineering and information science. Мoscow: МSАIE, 2005, pp. 19-23.

4. Schoeneburg E., Heinmann F., Feddersen S. Genetische Algorithmen und Evolutionsstrategien: Eine Einfuerung in Theorie und Praxis der simulierten Evolution. - Bonn; Paris; Readin; Mass. [u.a.]: Addison-Wesley, 1994.

LASER OPTICO-ACOU STICAL OCTANE METER IN SYSTEMS OF OPTIMIZING THE PROCESS OF COMMERCIAL PETROL COMPOUND-FILLING

© 2009 V. N. Astapov Samara State Aerospace University

The paper presents a methodology of constructing systems of optimal commercial petrol compound-filling. A mathematical model of the “octane number” and a method of identification of the octane numbers of the components being mixed using a laser optico-acoustical octane meter are given. For this purpose the inverse problem of petroleum mixture analysis using a genetic algorithm is solved.

Laser optico-acoustical octane meter, optimization, compound-filling, commercial petroleum, identification, mathematical model, spectral channels of measurement, calibration.

Информация об авторе Астапов Владислав Николаевич, кандидат технических наук, доцент кафедры электронных систем и устройств, Самарский государственный аэрокосмический университет, email: asta_09@mail.ru. Область научных интересов: устройства контроля качества и автоматизированные системы управления в нефтеперерабатывающей промышленности.

Astapov Vladislav Nicolayevitch, candidate of technical science, associate professor of the department «Electronic systems and devices», Samara State Aerospace University, e-mail: asta_09@mail.ru. Area of research: quality control devices and automatic control systems in petroleum industry.