CC BY
60
УДК 663.2+665.2/002.3:681.513.6(043.3)
АЛГОРИТМ ОПЕРАТИВНОГО УПРАВЛЕНИЯ ПЕРЕРАБОТКОЙ ПАРТИЙ СЫРЬЯ НА МАСЛОЭКСТРАКЦИОННОМ ЗАВОДЕ
процессов
Кубанский государственный технологический университет, Краснодар, Россия
В статье предлагается алгоритм оперативного управления переработкой партий сырья на маслоэкстракционном заводе
Ключевые слова: СЫРЬЁ, УПРАВЛЕНИЕ, ПРОИЗВОДСТВО, МОДЕЛЬ, ПЕРЕРАБОТКА, ОПТИМИЗАЦИЯ, МАСЛОЭКСТРАКЦИОННОМ ЗАВОД
UDC 663.2+665.2/002.3:681.513.6(043.3)
ALGORITHM OPERATIONAL CONTROL OF PROCESSING OF LOT OF RAW MATERIAL AT OIL EXTRACTION PLANTS
Asmaev Michael Petrovich
professor of the Department of Automation of
production processes
Kuban State Technological University,
Krasnodar, Russia
This article offers the logical probabilistic approach of administrative decisions acceptance in electronic libraries creation
Keywords: RAW, MANAGEMENT, PRODUCTION, MODEL, PROCESSING, OPTIMIZATION, AND OIL EXTRACTION PLANT
Асмаев Михаил Петрович
профессор кафедры автоматизации производственных
В настоящее время особую актуальность приобретают задачи оптимизации сельскохозяйственного производства на различных этапах его функционирования [1,2]. Широким классом указанных задач являются многовариантные задачи управления технологическими процессами, характерной особенностью которых является необходимость перебора большого числа вариантов для определения оптимального графика технологического процесса. Поэтому процедура оптимального решения крайне трудоемка и практически не реализуема на современных ЭВМ для большинства многовариантных задач [1]. Это объясняется тем, что время реализации алгоритма растет экспоненциально с ростом размерности задачи.
Пищевая промышленность страны представляет собой одну из крупных и сложных отраслей в структуре экономики. Значительное место в общем объеме производства пищевой промышленности занимает переработка скоропортящегося сельскохозяйственного сырья.
Анализ большого числа производств, перерабатывающих скоропортящегося сельскохозяйственную продукцию, позволил
представить их обобщенную структуру в виде двухэлементной схемы: сырьевое отделение - цех переработки.
Объект хранения в сырьевом отделении - скоропортящееся сельскохозяйственное сырье - является сложной биологической системой, подверженной влиянию различных случайных факторов.
Исходя из сути протекающих при хранении превращений и методов организации контроля, процесс хранения скоропортящегося сырья рассматривается как последовательные переходы из начального состояния (момент поступления на предприятие) в конечное (момент подачи на промышленную переработку ) под действием совокупности внешних и внутренних факторов. Сам процесс перехода носит устойчивый характер, протекает сравнительно медленно, в связи с чем может характеризоваться как квазистатический; параметры, определяющие процесс перехода носят случайный характер.
На основе общности закономерностей квазистатических процессов перехода для скоропортящегося сырья предлагается характеризовать их условными кинетическими характеристиками - изменением функции потерь во времени. Эти характеристики можно считать первыми приближенными аналогами процесса хранения при анализе и классификации их. Нестационарный случайный характер процессов хранения, отсутствие достаточно оперативной информации об изменении качественных и количественных характеристик сырья под влиянием как внешних, так и внутренних воздействий предопределяет использование адаптивных методов идентификации с прогнозированием. Для идентификации предлагается использовать два вида моделей с прогнозированием стохастический и детерминированный варианты.
В качестве алгоритмической основы прогнозирования используются реккурентные соотношения, позволяющие проводить адаптацию в темпе с технологическим процессом.
Формализованное описание задачи идентификации представляется как движение точки (параметры сырья) внутри полосы, границами которой служат кинетические характеристики (предельные значения функции потерь). Положение точек внутри полосы ограничений определяется параметрами модели возмущений. Общая траектория такого перехода определяется условными кинетическими характеристиками, а вероятностные состояния объекта в каждый момент времени зависят от влияния неконтролируемых стохастических воздействий. Проведенные исследования показывают, что действующие на объект возмущения относительно малы и, следовательно, невелики отклонения траектории реального случайного процесса от эталонной траектории, которая является решением детерминированной модели. На основе этого случайный процесс рассматривается как линейная функция от соответствующих значений эталонной детерминированной траектории
П (дг) = Л(г) • хт + %(г), (1)
где П(А1) - потери сырья и ценного компонента в нем;
А(11) - матрица параметров модели;
^текущее время;
X - момент времени начала хранения;
^ - значение возмущений в виде шума, некоррелированного с параметрами модели.
В свою очередь:
Л(г) = [ао(г), ^(г), ^)], М = г - гх, ,
хт = (1, Мг, Мг2),
где П(А1) - потери сырья и ценного компонента в нем; а;(1=0,1,2) - параметры модели.
Задача идентификации параметров функции потерь в этом случае формулируется следующим образом [2]. В начальный момент t0 известны значения элементов матрицыA(t0).По данным лабораторного анализа поступающим в момент времени ti(i=1,2,...,n),n -число шагов идентификации в виде значений ni=n(ti-tix), можно спрогнозировать значение ni+1 на момент времени ti+1 посредством идентификации элементов матрицы A(ti) так, чтобы погрешность между прогнозом ni+1 была минимальной в смысле критерия J:
n
J = I (П +1 -П1) ® min, (2)
i +1
где A(ti) - матрица параметров идентифицированных в момент времени ti.
В свою очередь:
Пг +1 = A(h) ,
Dt = ti +1 - ^, ,
XT = (1, Dt, Dt2).
Для определения элементов матрицы Ai=A(ti) был использован метод обобщенных фиктивных помех для линейных моделей с независимым шумом.
Дискретные значения моментов времени ti определялись из соотношения:
тг = t0 - (i - 1)h,
t0 = 0, , (3)
i = 1,..., n,
где h-шаг идентификации; n - номер шага.
Начало хранения сырья принимается равным времени начала процесса идентификации, т.е. 1;х=1;0=0.
Анализ процесса хранения скоропортящегося
сельскохозяйственного сырья позволил принять гипотезу, что шум ц-случайная нормально-распределенная величина. В исследованиях ц присваивали случайные значения из интервала (а,Ь), причем полагали
а = 0,
Ь = П ■ х п х 0,1,,
I ' "
цооЫ (0,1).
Проведенный анализ показывает, что для синтеза систем управления подобного класса объектов целесообразно использовать принцип детерминированной эквивалентности, основным преимуществом которого является возможность декомпозиции сложной стохастической задачи на две более простые - стохастического оценивания (прогнозирования) и детерминированного управления. Целевая функция задачи оптимального управления ТК, организованного по схеме однопродуктового потока имеет вид:
п
(4)
П Х (АТХ) = а01 + аи + а2, (АТЦ )2, (5)
Ь- • м..
П (АТЦЦ) = Ьи + ехр(— ), (6)
и и аТП
ц
Д грп _ грП грХ {Н\
[Ц = ТИ -ТИ , (7)
где ATnij,TXij - соответственно моменты времени начала переработки и хранения ij -й партии сырья;
ij - номер партии перерабатываемой в i -ю очередь;
n - общее число партий сырья в цехе хранения;
Mij - масса ij -й партии сырья;
a0i,a1ia2i вш, в^- параметры модели;
Qj - стоимость сырья и ценного компонента в нем.
Ограничения задачи вытекают из физической сути протекающих технологических процессов.
Задача управления ТК формулируется следующим образом:
* . * * . * ^
определить последовательность p j=j bJ 2,...,j n переработки партий сырья и соответствующую этой последовательности оптимальные моменты времениТ^*,Т^2*,... ,TPjn*, чтобы
n
&(Tß*,Tj2*,..., j *) = mm = CJj [ П x (DTx) + Пn (ATn)], <П j = Ji, J2'...' Jn >'
ТТПЛ ,Tn,...,Tn >, (8)
4 i1' i2' ' in
tmin £ Tn £ tmax
liJ £ TiJ £ liJ '
где Pj=j1,j2,.,jn- некоторая произвольная последовательность подачи сырья в промышленную переработку (если 1,2,...,n -номера партий, то j1,j2,. .,jn некоторая перестановка из элементов 1,n).
Большая размерность задачи не позволяет использовать классические оптимизационные методы. В связи с этим разработан приближенный метод решения при условии, что параметры модели -детерминированные величины.
На первом шаге методом случайной направленной выборки генерируется последовательность Pj переработки партий сырья, для
формирования которой использовалась идея «цепного» метода. Определяется расстояние между перестановками как минимальное число транспозиций необходимое для перехода от одной перестановки к другой. К - окрестностью перестановки называется множество перестановок, для которых расстояние до заданной перестановки не больше К. Суть метода в следующем. Выбрав в качестве начальной последовательности некоторую перестановку р0 в первой окрестности р0, ищем перестановку р1 такую, чтобы Ф(р1)<Ф(р0). Если такая перестановка отсутствует, то поиск ведем во второй окрестности р0. После определения последовательности р1 осуществляем поиск перестановки р2 в первой окрестности р0. Продолжая поиск нужное число раз, получим требуемую совокупность последовательностей {р}. Начальная последовательность р0 формируется посредством ранжирования номеров партий в порядке убывания скорости порчи сырья.
Для оценки эффективности алгоритма предложено использовать коэффициент направленности
* = Мф КА =-
М
ф (9)
где Мф, Мф - математическое ожидание соответственно величин Ф(р.), Ф(ж.), ]=1,к.
При КА^0 алгоритм случайной выборки назовем алгоритмом случайной направленной выборки.
На втором шаге приближенным методом динамического программирования определяются оптимальные для данной
последовательности р моменты времени Ц-, Ц-,..., -,
минимизирующие целевую функцию (4) при заданной последовательности Прямой ход метода состоит в получении рекуррентных зависимостей ф;:
тПо=ф (ТП+1)' (10)
где i=1,n-1,T0ij - условно- оптимальные управления такие что
Ф(ТЯ ТП20'..,ТП +1) = тп Ф(ТД0' ТП20'...' т"ю-Т], +1>. (11)
Неоднократная последовательность реализации шагов 1,2 позволяет
из совокупности (как правило достаточно нескольких десятков)
*
сгенерированных последовательностей выбрать лучшую л^^рь
j*2, — ,j*n- с меньшим значением целевой функции (4), При этом оценивается вероятность ¥ того, что найденное решение Ф(я^) окажется в V - окрестности оптимального Фопт, т.е.
Ф(р*) - Фопт
-1 * °Пт <п. (12)
Ф(р0)
Значение V задается из практических соображений, а значение ¥ определяется из соотношения
(1 ^ 2 ¥ = 1--,
12 )
(13)
где ъ - число значений целевой функции Ф^), попавших в V - окрестность минимально из найденных Ф(я^).
Соотношение (13) следует из возможности аппроксимации закона распределения Ф(^) в области малых вероятностей равномерным. Справедливость такого решения подтверждена экспериментом.
Приближенные методы оптимизации позволили разработать алгоритм двухуровневой оптимизации процесса переработки и хранения сырья. На первом внешнем уровне алгоритм случайной направленной выборки определяет последовательность поступления партий сырья на промышленную переработку, на втором внутреннем - приближенным методом динамического программирования определяются моменты времени начала переработки каждой партии сырья.
Адаптивный алгоритм оперативного управления рассматриваемым классом объектов сводится к детерминированному поиску значений управления в квазистатических состояниях объекта. Однако, в отличие от известных задач вариационного исчисления, здесь в качестве искомых параметров, характеризующих управляющие воздействия, выбираются не моменты переключения, а абсолютные значения управлений на заданных интервалах. При этом расчет осуществляется после вычисления соответствующих траекторий предполагаемого движения объекта. Здесь предполагается использовать метод прогнозирующего управления как вычислительную процедуру для реализации адаптивных алгоритмов, состоящих из двух этапов.
На первом этапе осуществляется поиск детерминированного или стохастического оценивания (в зависимости от принятой модели процесса хранения) значения потерь сырья при хранении на каждый последующий шаг ti+1 по данным текущих измерений. На втором этапе определяются абсолютные значения оптимальных управлений для рассчитанных значений прогнозируемых потерь Zn(ti+1)®min в соответствии с описанными ранее алгоритмическими процедурами. Указанная вычислительная процедура повторяется всего интервала [0,T] переработки очередной партии сырья. Для остальных партий, находящихся в очереди на переработку, с помощью адаптивных алгоритмов уточняются значения текущих потерь. В момент окончания переработки возникает
необходимость решения задачи по выбору очередной партии сырья. К этому времени за промежуток [0,Т] условия и ход процесса изменились. Управляющее устройство, имея информацию о величине потерь в каждой партии, решает эту задачу. Рассмотренная последовательность операций продолжается до окончания переработки всей массы сырья.
Литература
1 Опыт создания и эксплуатации АСУТП маслодобывающего предприятия на базе Днепропетровского маслоэкстракционного завода / Я.М. Лесов [и др.]. ВНИИТЭИАгропром. Серия Масложировая промышленность. Вып. 5. М., 1987. 24 с.
2 Асмаев М.П., Пиотровский Д.Л. Автоматизированные информационно-управляющие системы. Краснодар, ГОУ ВПО КубГТУ, 2009. 192 с.
