Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Типовые задачи оперативного управления непрерывным производством

..pdf
Скачиваний:
6
Добавлен:
12.11.2023
Размер:
19.07 Mб
Скачать

деление плановых заданий по интервалам планиро­

вания.

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

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

1)проверяется, имеет ли решение задача технико­ экономического планирования с заданными плановыми показателями и неограниченными запасами компонен­ тов. Отсутствие решения говорит о необходимости пере­ смотра общего планового задания;

2)проверяется, имеет ли решение задача технико­ экономического планирования с заданными плановыми заданиями и агрегированными поставками сырьевых компонентов. Отсутствие решения означает необходи­ мость корректировки величин поставок сырья и плано­ вых заданий с учетом решения задачи 1.

Кроме того, решение задачи 2 позволяет оценить не­ равномерность заданного графика поставки сырьевых компонентов и, если возможно, скорректировать этот график;

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

мости изменения агрегированных плановых заданий на весь планируемый период. Если решение существует, оно дает распределение плановых заданий по интерва­ лам и возможность скорректировать заданные плановые задания;

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

Более сложные процедуры анализа и устранения не­ совместности (см. § 6-3) не вводились в постоянное алго­ ритмическое обеспечение системы.

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

 

 

 

 

 

 

 

 

 

 

 

Таблица 6-3

 

Показатель качества

 

 

 

 

 

Марка бензина

 

 

 

 

 

 

 

1

 

2

 

3

 

 

 

 

 

 

 

 

 

 

 

Октановое число

%

 

 

 

 

76

 

72

58

 

 

Содержание серы,

 

 

 

 

 

0,14

 

0,14

0,30

Цена, руб/т

 

 

 

 

 

 

33

 

29

23

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица

6-4

Показатель качества

 

 

 

 

 

Компонент

 

 

 

 

I

 

2

3

 

4

5

6

7

1 8

 

 

 

 

 

 

Октановое число

%

49

70

83

 

75

58

62

72

 

45

Содержание

серы,

0,05

1,00

0,40

0,05 0,05

0,02

0,02

 

0,02

Себестоимость, руб/т

15

28

29

 

30

20

23

27

 

16

 

 

 

 

 

 

 

 

 

 

 

Таблица

6-5

Интервалы

 

 

 

 

Компонент

 

 

 

 

 

1

2

 

3

4

 

5

 

6

7

 

8

 

 

 

 

 

 

 

 

I

1410

2200

 

860

910

 

 

500

2020

 

36

 

и

 

920

1120

 

800

900

 

550

2060

 

 

ш

1240

1140

 

910

 

 

 

480

2580

 

42

 

 

 

430

 

 

 

IV

1450

970

 

840

 

 

520

3350

 

37

Объем емко­ 1500

2000

1000

1500 1000

1000

2000

 

100С

стей, т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица

 

6-6

Марка

 

 

 

 

Компонент

 

 

 

 

Общее ко­

бензина

1

2

 

3

4

 

5

 

6

 

7

личество,

 

 

 

 

 

8

 

т

1

 

2172

1810

 

 

 

3348,5

 

 

7330

2

813,6 471,4

254,2

 

6661,5

 

8200

?

5020

2730

766,6

175,8

 

2050

1150

11 890

Значения октановых чисел, содержание серы в исходных ком­ понентах и их себестоимость представлены в табл. 6-4. Предпо­ лагается, что эти значения являются усредненными и не изме­ няются на всем периоде планирования.

на

Период

планирования

разбит

 

I

К

Ш

 

W

четыре

временных

интервала.

 

 

 

 

 

 

 

В табл. 6-5 заданы

 

графики

 

 

 

 

 

 

поступления

исходных

 

компонен­

 

 

 

 

 

 

тов во время каждого из четырех

 

 

 

 

 

 

временных

интервалов.

 

 

началь­

 

 

 

 

 

 

ные

Предполагается,

что

 

 

 

 

 

 

запасы компонентов

равны

 

 

 

 

 

 

нулю.

ограничения на объ­

 

 

 

 

 

 

 

Верхние

 

 

 

 

 

 

емы емкостного парка для каж­

 

 

 

 

 

 

дого из компонентов

приведены

 

 

 

 

 

 

в последней строке табл. 6-5.

 

 

 

 

 

 

 

 

Решение задачи

текущего пла­

Рис.

6-7.

Решение

детализиро­

нирования

для всего

 

планируе­

мого периода по критерию макси­

ванной

задачи

планирования

мизации

прибыли

 

приведено

выпуска

трех

марок

бензина.

в табл. 6-6.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Согласно этом> решению необходимо выпустить бензина марки

1—7330 т, 2—8200 т и 3—М 890

т.

Прибыль

при этом

достигает

104 430 руб.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 6-7

 

 

Марка бензина

 

 

Марка бензина

 

 

 

 

 

 

 

 

 

Остаток

 

 

 

 

" Остаток

Компонент

1

 

2

 

 

3

в емко­

1

2

 

3

в емко­

 

 

 

стях, т

 

стях, т

 

 

 

 

I интервал

 

 

 

II интервал

 

 

1

 

 

149,8

 

1410.0

0

 

 

218,8

920.0

0

 

2

498.3

 

 

751.4

1299,0

 

643,7

1556,5

 

3

 

27,2

 

 

334,5

 

176.8

 

 

957,7

 

4

415.3

 

 

 

 

 

494,7

 

 

 

1394,7

 

5

 

 

 

 

 

500.0

0

 

150.8

 

 

0

 

6

768.3

1132,2

 

0

 

 

399,2 ,

 

7

 

119.5

0

 

 

1800,0

259.0

$

 

8

 

 

 

 

 

 

0

 

 

 

360.0

0

Общее количе­

1681,9

1309,2

 

2780,9

 

 

2346,4

2582,8

 

ство

 

 

 

 

 

 

 

 

 

 

 

 

 

Компонент

 

 

III

 

интервал

 

 

 

IV интервал

 

1

 

 

1240,0

0

о

, -,

820,7

1875,8

3

1073,4

194,3

0

4

1394,7

—-

0

5

 

430,0

0

и

.

.—,

480,0

0

7

2580,0

.—1

0

8

 

 

420,0

0

 

268,7

1181,3

0

308,0

648,0

1889,8

617,9

222,0

0

— .

— .

0

520,0

0

0

3349,8 0

370,0 0

 

 

 

 

"

Общее количе­ 5648,1

.—

3585,0

4544,4 2941.3

ство

 

 

 

 

Для детализации полученного плана по четырем временным интервалам была применена описанная выше двухинтервальная про­ цедура. В результате рассмотрения трех двухинтервальных задач получено решение детализированной задачи, представленное в табл. 6-7. Как видно из рис. 6-7, полученное оптимальное де­ тализированное решение существенно отличается от равномерного по интервалам IIV.

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

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

6-7. Имитационное моделирование календарного планирования

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

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

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

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

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

Под имитационным моделированием понимают про­ ведение на ЭВМ численных экспериментов с моделями сложных систем [100] с целью предсказания поведения системы на основе поведения модели.

Для моделирования необходимо:

1)построить модель объекта, выделить управляемые

инеуправляемые переменные, указать параметры мо­ дели в некоторый момент времени;

2)описать способ изменения этих параметров во вре­ мени под влиянием неуправляемых переменных, изме­ няющихся случайно или детерминированно. Задать про­

грамму изменения детерминированных неуправляемых переменных;

3) задать алгоритм (или несколько алгоритмов) вы­ бора управляемых переменных (алгоритмы управления).

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

Эти задачи решаются путем моделирования (провер­ ки) последствий решения при некоторой реализации зна­ чений случайных параметров модели. Затем процедура повторяется для других реализаций значений случайных параметров модели. Таким образом, при каждой реали­ зации решается детерминированная задача.

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

В качестве примера алгоритма решения детермини­ рованной задачи можно указать на алгоритм скользя­ щего метода планирования (см. § 6-5), когда задача календарного планирования при случайных параметрах качества сырья, поставок и сбыта решается последова­ тельно с помощью двухиитервальной процедуры, начи­ ная с первого интервала. Имитационное моделирование

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

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

Рис. 6-8. Пример схемы промыш­

ленного

предприятия, состоящего

из трех

технологических

комплек­

сов.

 

 

А — запасы продуктов; /. 2,

3 — ком­

плексы.

 

 

Для календарного планирования моделируемые

алгоритмы должны быть достаточно простыми,

чтобы

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

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

Решим для данного производства задачу технико­ экономического планирования (см. § 6-1). В этой задаче используется модель, которая не учитывает запасов про­ дуктов между комплексами. Требуется найти оптималь­ ный текущий план при средних значениях материальных потоков на входе и выходе отдельных комплексов и сред­ них значениях других параметров модели (производи­ тельности агрегатов, качества продуктов и т. п.). Ука­ занная задача представляет собой задачу линейного программирования. Переменными являются материаль­ ные потоки Xitu—/г-й компонент (или поток), направ­

ляемый от /-го комплекса к s-му. Решение — оптималь­ ный технико-экономический план Х*ма—определяет все материальные потоки.

Теперь рассмотрим задачу календарного планирова­ ния. С целью упрощения модели, а следовательно, и алгоритма, предположим, что в течение всех интервалов времени для каждого производственного комплекса со­ храняются соотношения между всеми материальными потоками, полученными в результате решения задачи технико-экономического планирования. Тогда в качест­ ве переменных модели календарного планирования бу­ дем рассматривать только интенсивности использования Яг/, i= l, ..., N, t= 1, ..., Г; Я// обозначает интенсивность использования i-ro производства в t-й временной интервал.

Переменных Я« значительно меньше, чем переменных Xiks для всех интервалов. Значения всех потоков в f-м интервале могут быть представлены в виде

Критерий задачи

T N I N к,

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

Ограничения задачи включают запасы в каждом интервале zPt, полученные с учетом запасов предыду­ щего интервала zvt-1

O^Zp/^ZpMaKc, t= 1, ...,

T,

p= 1, ..., P\

Zpt----Zpt-i~\~ 2

— 2

 

n't"

n ~

T

QI

й нагрузки произ-

внешних потоков ^

ВОДСТВ

где 2Рмакс— максимальные запасы р-го продукта; р+— множество входных потоков для продукта р; р~— мно­ жество выходных потоков для продукта р; Я/мин — мини­

мальная

нагрузка производства; Q* — внешние потоки

(сырье и готовые продукты).

Такое

упрощение задачи позволяет исключить все

ограничения, связанные со структурой производства. Эффективность полученного плана и его реализуе­

мость можно оценить путем имитационного моделиро­ вания работы производства.

В том случае, когда полученная модель, содержащая NXT переменных А,«, имеет все же большую размер­ ность, то возможно дальнейшее упрощение путем отка­ за от решения задачи для всего горизонта планирования и решение подзадач последовательно для каждого t-ro интервала. Для каждого интервала задача имеет неболь­ шую размерность, число переменных равно числу ком­ плексов. Например, для схемы на рис. 6-8 число пере­ менных равно трем. При моделировании такой задачи на каждом интервале можно учитывать влияние случай­ ных переменных: изменение производительности обору­ дования, поставок сырья, отгрузки готовой продук­ ции и т. п.

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

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

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

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

В результате моделирования этого алгоритма на ЭВМ с использованием множества реализаций находит­ ся приемлемый календарный план.

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

вкоторых значительная часть оборудования находится

вплановом ремонте.

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

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

Пограммирование имитационных моделей на ЭВМ осуществляется с помощью универсальных языков типа ФОРТРАНа, АЛГОЛа, PL/1, Ассемблера или специали­ зированных языков моделирования. При применении универсальных языков для выполнения различных стан­ дартных операций, например вычисления математиче­ ских ожиданий и дисперсий переменных модели, получе­ ния чисел с равномерным, нормальным или другим за­ коном распределения, используются подпрограммы. Обычно имеется библиотека стандартных подпрограмм, написанных на ФОРТРАНе. В программе моделирова­ ния описываются алгоритмы функционирования системы и .имеются обращения к соответствующим подпрограм­ мам. Программы имитационного моделирования доста­ точно сложные. В различных программах имитационно­ го моделирования часто выполняются одинаковые алгоритмы. В связи с этим были разработаны специали­ зированные языки для программирования имитацион­ ных моделей: СИМСКРИПТ, DYNAMO, CSL и др.

Программа моделирования на таком языке короче, чем программа на ФОРТРАНе, и может быть быстрее отлажена. Большие и средние ЭВМ третьего поколения обычно имеют в составе программного обеспечения язык

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

6-8. Календарное планирование работы периодических производств

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

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

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

Т[1* (время начала его первого рабочего цикла).

Введем'следующие обозначения: T{kl)— время начала k-ro рабочего цикла t-ro а гр е га т а ;^ —продолжитель-

Соседние файлы в папке книги