МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
УДК 519.21
А.С.Сорокин
О ВЫБОРЕ ОПТИМАЛЬНЫХ ПАРАМЕТРОВ ТРАНСПОРТИРОВАНИЯ И ОБОРУДОВАНИЯ ПРИ ГИДРОТРАНСПОРТЕ
В данной работе предлагается общий подход к выбору оптимальных параметров оборудования, применяемого при дальнем гидротранспорте, и оптимальных параметров самого процесса транспортирования.
Основной задачей при проектировании является получение либо минимальных суммарных относительных удельных затрат с учетом амортизационных отчислений; связанных с абразивным износом трубопровода и оборудования [1], либо максимальной экономии металла при капитальных минимальных затратах на сооружение трубопроводов и минимальных издержках на их эксплуатацию [2], либо полной максимальной прибыли при работе системы в целом в течение, например, года с учетом надежности системы и затрат на ее эксплуатацию [3] .
В работах [4 , 5] для оперативного определения оптимальных режимов работы системы транспортирования в качестве критерия предлагается минимум удельной энергоемкости транспортирования (величины удельных затрат энергии на транспортирование одной тонны твердого материала на расстояние одного километра по трубопроводу заданного диаметра).
Нам представляется искусственным разделение выбора оптимальных параметров процесса транспортирования и выбора оптимальных параметров оборудования. Более естественным является построение совместной модели транспортирования, сущность которой составляют четыре взаимосвязанные между собой задачи нелинейного программирования.
Опишем конструктивное построение этой модели.
Входные параметры модели: необходимая производительность системы транспортирования Q , протяженность транспортирования L, топология трассы транспортирования (к, ^, ] = 1, Р), вектор экономической информации , вектор информации о надежности Ц2 , вектор технологической информации Цз , вектор гидравлической информации Ц4, вектор прочностной информации Ц5.
Выходные параметры модели: Р- оптимальное число перекачивающих станций;
(,] = 1, Р) — оптимальные напоры на каждой
станции; Qj - оптимальные производительности одного агрегата на каждой станции; пг- - оптимальное число рабочих агрегатов; тг- - оптимальное число резервных агрегатов; гг- - оптимальное (ожидаемое) число агрегатов, находящихся в ремонте, на каждой перекачивающей станции; V -оптимальная скорость транспортирования гидравлической смеси; й - оптимальный диаметр трубопровода; 8 - оптимальная толщина стенок труб; $ - объёмная оптимальная концентрация пульпы.
В результате проигрывания модели получаем: координаты размещения перекачивающих станций вдоль трассы транспортирования (к], 1\, i = 1,Р); коэффициенты готовности перекачивающих станций (КГ, ] = 1,Р); полную прибыль при работе системы транспортирования в целом с учетом надежности и затрат на её эксплуатацию Г ; приведенные затраты, учитывающие капитальные вложения, эксплуатационные расходы и коэффициент эффективности
С } +ЕК общую мощность, потребляемую системой транспортирования Ж.
Укрупненная блок-схема алгоритма поиска оптимальных параметров оборудования и процесса транспортирования при дальнем гидротранспорте представлена на рис. 1 (прямоугольниками изображены подсистемы, каждая из которых имеет внешние входы и потоки компонентов из всех остальных подсистем; из каждой подсистемы выходят компоненты, направляемые во все подсистемы и внешние выходы).
Приступим к детализации общей топологической структуры системы транспортирования, выделяя из общей схемы каждую из четырех подсистем.
Подсистема резервирования цепи транспортирования (Первая подсистема)
Математическая модель подсистемы резервирования цепи транспортирования содержит следующие функциональные ограничения:
ат.. + Ьг < С, ] = 1, Р,
Рис. 1. Информационная блок - схема алгоритма определения оптимальных параметров оборудования и процесса транспортирования
где а,, Ь], С] - компоненты Ц}.
Целевой функцией для подсистемы резервирования является коэффициент готовности .
-1—1
(24п)т+1 — 5т+1
КГ =
0,17 • 24т пт+1(552п + 29) — 5 • 5п
К Г (п,,т,,г, )=
1+-
1
гр2 Еак(т1,п)р1 к=0
КГ =
г((5г)т+1 — 24т+1) 5 ^ 5тгт+2 — 24т(4,83г + 92)'
причем
Р,
[х +Хп +^п у
Р2
п
М +»'„ (+х)
где Х,Х„ , М , Мп — компоненты вектора Ц2. Кроме того, ак(т,,г) равно Г при 0<к
г!
т—к -; -
г при т-г < к < т.
<т-г и равно (т—к)!'
В заключение отметим, что величины т] и г,> 0 имеют целочисленные значения.
На рис. 2-4 представлены результаты моделирования подсистемы резервирования цепи транспортирования - зависимости коэффициента готовности как от числа рабочих углесосов, так и от числа углесосов, находящихся в ремонте, представляемые формулами:
г(24п + 5г)
1ЦИ6НТ готовности
КГ =
92п2+24пг+5г2
г, углесосы в ревкцнте _
^^^^^^"раВочие углесосы
О о
Рис. 2. Зависимость коэффициента готовности, как от числа рабочих углесосов, так и от числа углесосов, находящихся в ремонте
Подсистема процесса транспортирования при дальнем гидротранспорте (Вторая подсистема )
Рассмотрим ограничения, входящие в состав математической модели подсистемы процесса транспортирования при дальнем гидротранспорте.
Рис.3. Зависимость коэффициента готовности, Рис.4. Зависимость коэффициента готовности,
как от числа рабочих углесосов, так и от числа как от числа углесосов, находящихся в ремонте,
углесосов, находящихся в резерве
так и от числа углесосов, находящихся в резерве
дельта, толщина трув^гмм
15 " .9
Э,консистенция
Рис.5. Зависимость приведенных затрат от консистенции пульпы и от толщины трубопровода
Рис. 6. Зависимость приведенных затрат от скорости течения пульпы и от толщины трубопровода
80000
затраты, тыс. руб.
.30000
а, диаметр трубы, мм
Рис. 7. Зависимость приведенных затрат от скорости течения пульпы и от диаметра трубопровода
Пропускная способность системы в целом по твердому топливу (выход готовой продукции)
П> 40Т,
где QТ - нормативно обоснованное (технологически возможное) задание по выпуску готовой про-
Рис.8. Зависимость приведенных затрат от консистенции пульпы и от диаметра трубопровода
дукции.
Пропускная способность системы по пульпе
П2¥8 > 4QТ,
причем Q -необходимая производительность системы транспортирования.
Ограничение общему износу трубопровода:
ёУп < /,
где /- компонента вектора Цз [6].
Предположение о безаварийном течении процесса транспортирования приводит нас к неравенству (исключение забучивания трубопровода)
V6 > Её2$,
причем Е -компонента вектора Ц4 .
Ограничения по прочности трубопровода:
ёНн < 1.635<тТ8, , = 1, Р,
где Н{ и 0,'т — компоненты векторов Ц4 и Ц5 .
Помимо функциональных ограничений имеются позиционные ограничения:
на объемную концентрацию пульпы:
0.12<8 <0.8, на скорость транспортирования пульпы;
V <5.
Величины ё и 8 задаются в виде соответствующих дискретных рядов.
Целевая функция для подсистемы процесса транспортирования представляет собой приведенные затраты, учитывающие капитальные вложения в систему и её эксплуатационные расходы:
Г($ V, ё ,8) = с; + ЕК =
В;ё V2 + В2ё V + В3ё V +
В4ё V (1 — $) + В^п + В6 Ь8(ё + 8) +
Ба2 Ь + КёЬ + КЬ +
До Ь
• +
11/ 3/
с + Б її
2 пг [ б;2у‘ 74 с/4 + б;3дс ^ +
1
б;4 зу ^ а34+б;5 я 2а54 V “14].
В этом соотношении Б /, Б/, / = 12,15 яв
У У ^
ляются компонентами векторов Ц1 , Цз .
Результаты моделирования процесса транспортирования
На рис. 5 представлена зависимость приведенных затрат от консистенции пульпы и от толщины трубопровода:
с;+ ЕК = 1,35 • 104^2 + 3375^ +1052 +15755 +1,64 • 105 .
На рис. 6 - зависимость приведенных затрат от скорости течения пульпы и от толщины трубопровода:
С1 + ЕК = 16^14 + 3303 V +10,752 +16205 +1,6 • 105 .
На рис. 7 - аналогичная зависимость от скорости течения пульпы и от диаметра трубопровода:
С + ЕК = -
37
+ С2(3,2У + 7) +
с + 20 + а (0,1У1,4 + 83,8) + 519.
На рис. 8 представлена зависимость приведенных затрат от консистенции пульпы и от диаметра трубопровода:
, 37
С + ЕК =----------+
1 с + 20
й2(0,6ї2 + 0,15& + 7,23) + 84,8$ + 519 .
На рис. 9 - зависимость приведенных затрат от консистенции пульпы и от скорости транспортирования пульпы:
С1 + ЕК = 270052у + 6755У + 16,5У1’4 + 1035У +1,7 • 105 .
На рис. 1о аналогичная зависимость от толщины и от диаметра трубопровода:
С'+ ЕК = ■
37
- + 7,7$2 +
с+20 с1(10,55 +11,3) +10,552 + 4,5.
300000
зтраты, тыс.руб.
100000
0 скорость пульпы, м/с 5 □
Рис. 9. Зависимость приведенных затрат от консистенции пульпы и от скорости транспортирования пульпы
Рис.10. Зависимость приведенных затрат от толщины и от диаметра трубопровода
Подсистема функционирования системы дальнего гидротранспорта (Третья подсистема)
Математическая модель функционирования системы дальнего гидротранспорта имеет следующие ограничения:
1 < п, < N0 — т, — г, ' = 1,п . Величина № - максимальное число устанавливаемых агрегатов на ] -й перекачивающей станции. Кроме того, величины принимают натуральные значения.
Для этой подсистемы целевой функцией будет полная прибыль при работе системы транспортирования в целом с учетом надежности и затрат на ее эксплуатацию:
= Г1Е а‘п‘Кг —
' = 1
В1Е ап,[ — К г )— Е ат>— Е ь/р
причем
КГ = П КГ (пг , тг , Гг )•
Кроме того, величины Г; и В; - компоненты векторов Ц1 и Ц3 .
Подсистема распределения напоров вдоль трассы транспортирования (Четвертая подсистема)
Эта подсистема представляет собой математическую модель распределения напоров вдоль трассы транспортирования с несколькими перекачивающими станциями.
Модель содержит систему функциональных ограничений в виде двойных неравенств:
4 <Е*1 <В,,' = 1Р .
1=1
где
А = НП+, - Н о + 2 Н] .
1=1
Б, = НН - Но +2 Н],
1 = 1
причем Н0, Нг, НН - компоненты вектора q4.
Кроме того,
Т У
НТ = ^ I И/ + ж1/)
где И], 11 - координаты размещения перекачивающих станций, а остальные величины, входящие в эту формулу, являются компонентами вектора q4.
Далее, может быть полезно соотношение :
2
КГ = 1 + [п 70 У-------.
Г 70
Целевая функция подсистемы представляет собой общие затраты мощности, потребляемой всей системой транспортирования. Эта функция имеет различные представления в зависимости от типа гидравлических машин, применяемых для транспортирования полезного ископаемого.
В случае применения лопастных машин целевая функция принимает вид :
Р у П . 2 п г
102 п . г = 1 'г
где Иг - коэффициент полезного действия гидрав-
лических машин, а величины а. и
Ь1 -
ком-
поненты вектора Ц4 .
В случае применения поршневых машин целевая функция имеет вид:
Р
у п2
ж =2
102
і = 1
X
Подсистема определения производительности одного агрегата, на каждой перекачивающей станции представляет собой соотношение:
Qi = , ' = Гр-.
' п,
Подсистема подготовки информации о размещении перекачивающих станций вдоль трассы транспортирования описывается с помощью метода Шухова [7-9] .
От моделей подсистем можно перейти к модели дальнего гидротранспорта в целом. Это означает, что моделью дальнего гидротранспорта
становится совокупность моделей подсистем. Кроме того, нужно еще к этой совокупности присоединить ограничения, описывающие систему в целом.
Рассмотрим эти ограничения. Отметим, что каждая из первых четырех подсистем представляет собой самостоятельную задачу математического (нелинейного) программирования с различными целевыми функциями. Каждая из этих задач может иметь и может не иметь решения. Разрешимость задачи в целом существует только при условии разрешимости всех частных задач нелинейного программирования.
Далее, выходы отдельных подсистем являются одновременно входами для других подсистем. Отождествление некоторых входов и выходов различных подсистем порождает ограничения материального баланса потока.
В качестве критерия оптимальности для выбора оптимального числа перекачивающих станций
І = 1
выбран: найти наименьшее натуральное Р, при
котором все частные задачи нелинейного программирования системы транспортирования разрешимы одновременно.
Математическая постановка задачи: необходимо минимизировать функцию затрат Р на технологическую схему.
т-2 р\+1 р\+2
р=2 2
к=0 г=1+р1 / =1 + р\+!
2 с1х1,
(1)
где с,! - затраты на 1-е звено технологии, Ху -
альтернативный коэффициент, принимающий значения 0 или 1.
Ограничения на надежность:
т — 2 Рк +1 рк+2
ЕЕ Е|пКуХу> 1пк, (2)
к = 0 , = к + рк 1 = к+ рк+1 где К - предельный уровень надежности, Ку -надежность у- го звена технологической схемы.
Ограничения на технологические факторы (например, водоснабжение): т— 2 р\+1 рк+2
ЕЕ Е1п О/ Х1 < 1п Q, (3)
к=0 г=;+р\ ]=;+рк+1 где Q - предельный уровень технологического фактора, Qij - технологический фактор у- го звена технологической схемы.
Единственность перехода (единственность выбора одного из рк+1 вариантов звеньев) от
одного типа технологических процессов к другому обеспечивается условием:
рк+1 рк+2
2 2 X/ =1 к = 0 т - 2.
г = 1+ рк 1=1 + рк+1 Неразрывность технологической цепи:
рк+1 рк+3
(4)
Е Ху = Е Ху, 1 =1+р;+^ р;+2;
, =; + рк 1=1 + рк+2
к = 0, т — 3. (5)
Принадлежность к заданному графу определяется условиями:
0 < Хг] < 1 для р, 1 )е Г, (6)
= 0 для І )є С5 Г.
(7)
В такой постановке задача отыскания оптимальной технологии сведена автором к классической задаче линейного целочисленного программирования.
Необходимо отметить, что ограничения (2) и
(3) являются характерными для определенного класса задач, они могут отсутствовать или по одному или оба одновременно.
Постановка задачи допускает введение дополнительных ограничений (характерных только для данной конкретной задачи); наличие же условий
(4) - (7) является обязательным при постановке
каждой задачи линейного целочисленного программирования.
В последнее время получил развитие метод отыскания решений задачи целочисленного программирования, обладающий по сравнению с вышеописанными методами, большей быстротой сходимости - метод вектора спада [10] (см. также [12]). Сущность этого метода состоит в следующем.
Рассмотрим множество Мц векторов вида
2 = (2;,2я ) где 2 р = 1, Ц, - целые числа.
Расстояние между произвольными элемента-
ми
2(1) = (2?)41)) и 2(2Ц^2^)
я ) V 1 ’ я
множества Мц определим таким образом:
г(22), ^(2)=Е 2р;)— 2р2). (8)
р=
Нетрудно доказать, что величина г (2(;), 2[2)) удовлетворяет аксиомам расстояния.
Таким образом, множество Мц можно рассматривать как метрическое пространство с метрикой, определяемой по формуле (8).
Окрестностью Ьр(2) точки г £Мц радиуса р называем совокупность точек х пространства Мц, удовлетворяющих неравенству (2, х )< р.
Расстояние, определяемое по формуле (8), позволяет сделать вывод о том, что окрестность Ьр(г) состоит из точки 2 и множества точек Мц, отстоящих от нее на расстоянии 1,2,...,[р]
([р] - целая часть числа р). Поэтому в дальней-
шем можно ограничиться рассмотрением окрестностей целочисленного радиуса.
Пусть множество ^ Мц состоит из точек
г = (г;, 22,..., 2Ч ), удовлетворяющих условиям (2)
- (7). Точку 2 е Яц назовем точкой локального минимума функции Е(г) определенной в области ^ Мц, если существует такое р>0 , что
Е (2) > Е (2) для всех 2 ^ Ьр П| и множество Ьр 1^1 содержит хотя бы одну точку, от-
личную от 2.
В качестве вектора спада функции Е(2), определенной в области Яд, в некоторой окрестности ьр(2) точки 2 е Яц рассмотрим вектор
Я(2) = (а(2, 2(;))| Ар, 2(я))),
где 22 \ , = 1, £, - точки множества Ьр 1^1 Я ,
отличные от г.
Л(ъ,ъ(>) = р(г(г >)- р(г) і = 1, g.
Если все
Л(,Ъ(^)> 0, і = 1
то
Р (г )= тіп Р (г(г))
гі Є Ьр (г )пй,
Если точка г Є Я, не является минимумом функции Р(г), то при помощи вектора спада Я(г) можно найти такую точку г ^ Ьр ^ Я, ,
что Р(г1 )< Р(г), т.е. л(г, г1 )< 0.
Алгоритм отыскания экстремального решения задачи по методу вектора спада запишется следующим образом.
Шаг 1. Выбираем начальное приближение
ъ(0ЧгСи”)...,г(')є Я,.
Шаг 2. На каждом очередном п+1 - м шаге алгоритма осуществляем следующие действия. Шаг 2.1. Выбираем радиус г.
Шаг 2.2. Рассмотрим окрестность
Ьр(г(г)^М, радиуса г с центром в точке г(п.
(п)
Точки множества Ьр(г1'"п) обозначим таким
образом:
ъ
(п,Р1,трРа,тра ) =
_ (пр1,тяра,тра ) ъ (п,р1тдратРа) ^1 ’
V
ъ (ътрРа ,тра ) = І
ъх ' + т р , если / = рі, р рі ’ ■> /^і’
>і ’
если І Ф Рі,
где ^ р, = 1, д; Р; <Р2 < ... <Ра , тр, е{0,+1,—1,+2,—2,...,+р,—р};
, = 1, а; а = тт(р, ц) (9)
Шаг 2.3. Выбираем одну из точек окрестности ьр(2(п)) , удовлетворяющую условиям (2)-(7)
Л(ъ^,ъ
п ) 1п,р1,тяра ,тРа
) < 0,
(10)
где
Л(ъ(п ) ъ(п,Р1, тР1 ра ,тРа )) =
а т—2 рк+1 рк+2
Е(Е Е Е с1 (х1—4)).
,1 =; к = 0 , =; + рк 1 =; + рк+1
Обозначив эту точку через 2(п+1 , переходим к шагу 2.2. Заметим, что эта точка принадлежит множеству Яц. Это приводит к уменьшению значения целевой функции Рф по
(п)
сравнению с ее значением в 2 .
Шаг 2.4. Если ни одна из точек окрестности
Ьр (2 (”)) не удовлетворяет условиям (2) - (7),
(10), то при п > 1 точка 2(п^ есть искомое решение задачи, а при п = 0 :
а) 2[п ^ - точка минимума функции Р (2), определенной на Яд , если существует хотя бы одна
точка окрестности
отличная от ъ
(п)
удовлетворяющая условиям (2) - (7) одновременно;
б) в противном случае переходим к шагу 1, выбирая другое начальное решение.
В заключение заметим, что оптимизация сложных систем, имеющих контуры с обратными связями, систем с противотоками, какой является предлагаемая совместная модель транспортирования, представляет довольно трудную самостоятельную задачу [11-12].
Перечисленные выше выходные параметры модели составляют так называемое оптимальное компромиссное решение задачи нелинейного программирования с несколькими целевыми функциями. Это решение является оптимальным компромиссом между оптимальными относительно различных целевых функций решениями задач.
Итак, предлагаемый подход позволяет найти не локальные, а глобальный оптимум для параметров оборудования и процесса транспортирования при дальнем гидротранспорте.
Итак, установлено, что возрастание экономической эффективности функционирования (прибыли) обслуживаемой системы возможно при введении резервных элементов в i-ю подсистему. Наряду с увеличением годовых эксплуатационных затрат увеличивается вероятность исправного состояния системы и уменьшается вероятность аварийного простоя.
Предложенный алгоритм был реализован на ЭВМ в математическом пакете БЕЫУЕ6.10.
Разработанный метод выбора оптимального резервирования может быть использован при расчетах на надежность любых систем, где необходимо учесть влияние надежности на экономичность.
Разработанная модель позволяет оптимизировать кратность резервирования оборудования технологических блоков водоснабжения и гидравлического подъёма.
Построена совместная модель транспортирования. В результате проигрывания её получаем: координаты размещения перекачивающих станций вдоль трассы транспортирования, коэффициенты готовности перекачивающих станций, полную прибыль при работе системы транспортирования в целом с учетом надежности и затрат на ее эксплуатацию, приведенные затраты, учитывающие капитальные вложения, эксплуатационные расходы, коэффициент эффективности и общую мощность, потребляемую системой транспорти-
и
и
рования.
На предварительном этапе математического моделирования дальнего гидротранспорта моделью является совокупность моделей подсистем, каждая из которых представляет собой самостоя-
тельную задачу математического (нелинейного) программирования с различными целевыми функциями. Разрешимость задачи в целом существует только при условии разрешимости всех частных задач нелинейного программирования.
СПИСОК ЛИТЕРАТУРЫ
1. Покровская В.Н., Бойцов В.С. Определение экономической целесообразности параметров гидравлического трубопроводного транспорта //Строительство трубопроводов. №2, 1972, с. 18 -19.
2. Кшондзер Э.Г. Определение оптимальных параметров магистральных трубопроводов для гидравлического транспорта сыпучих материалов //Строительство трубопроводов. №9, 1972, с. 15-16.
3. Сорокин А.С., Разгильдеев Г.И. Выбор оптимальных параметров процесса транспортирования и оборудования при дальнем гидротранспорте / Труды ВНИИГидроуголь, вып. 35, Новокузнецк, 1975. с. 86 - 94.
4. Спиваковский А.О., Смолдырев А.Е., Зубакин Ю.С. Автоматизация трубопроводного транспорта в горной промышленности. - М.: Недра, 1972. - 343 с.
5. Стахеев Н.А., Витебский В.Я. Анализ и оценка технологических систем на основе теории надежности // Изв. ВУЗов. Горный журнал, № 12, 1966.
6. Силин Н. А., Витошкин Ю. К. Гидротранспорт угля по трубам и методы его расчета. Киев. 1964.
7. Черникин В.И. Перекачка вязких и застывающих нефтей. М., 1958. - 234 с.
8. Яблонский В.С., Новоселов В.Ф., Галлеев В.Б., Закиров Г.З. Проектирование, эксплуатация и ремонт нефтепродуктопроводов. - М.: Недра, 1965. - 323 с.
9. Сорокин А.С. Применение методов теории вероятностей к исследованию некоторых процессов производства. Труды 4-ой междунар. конф. Кибернетика и технологии XXI века. Воронеж, 2003. с. 312323.
10. Архипова Т.Т., Рощин В.А., Сергиенко И.В. О решении одной задачи целочисленного программирования. Кибернетика, № 1, 1973.
11. Чуев Ю.В., Спехова Г.П. Технические задачи исследования операций. М., 1971. - 223 с.
12. Сорокин А.С. Алгоритм решения систем уравнений Колмогорова (Оценка качества системы) //Вторая Всероссийская научная конф. Проектирование инженерных и научных приложений в среде МЛТЬЛБ. - М., 2004. с. 389 - 397.
□ Автор статьи:
Сорокин Андрей Семенович
- канд. физ.-мат.наук, доц., ст.научн. сотр. ( филиал КузГТУ, г. Новокузнецк)
УДК 622:272
В.В. Крюкова
НОВЫЙ подход к моделированию процессов и подсистем
УГОЛЬНОЙ ШАХТЫ
В основе указанного подхода лежит методология имитационного моделирования сложной системы «шахта» на сетях Петри (СП). Под методологией здесь понимаем набор методов, применяемых для описания и исследования процессов и подсистем горного производства угольной шахты, объединенных единой концепцией. Основа концепции - объектноориентированный взгляд на процессы сложной производст-
венной системы «шахта» с позиции теории СП.
Метод - определенная процедура создания моделей с помощью четко специфицированных нотаций. Нотация - графический язык для описания и представления моделей процессов и подсистем, в ее основе модифицированные сети Петри.
Исследуемая сложная система «шахта», ее подсистемы и процессы в данной методологии рассматриваются с разных
точек зрения. Результат такого рассмотрения и есть модель подсистемы или процесса.
Формально методологию можно представить набором М={{А, N}, I}, где А - архитектура (строение) полуформальных моделей с графическими нотациями с учетом объектной ориентированности; N - нотации, используемые для описания конкретных моделей процессов и подсистем; I -инструментальные средства, -