Научная статья на тему 'Альтернативная модель приоритетного обслуживания для оценки пропускной способности аэродрома'

Альтернативная модель приоритетного обслуживания для оценки пропускной способности аэродрома Текст научной статьи по специальности «Компьютерные и информационные науки»

CC BY
128
38
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ОРГАНИЗАЦИЯ ВОЗДУШНОГО ДВИЖЕНИЯ / АЭРОДРОМНЫЕ ОПЕРАЦИИ / ПРИОРИТЕТНОЕ ОБСЛУЖИВАНИЕ / ВРЕМЯ ОЖИДАНИЯ ОБСЛУЖИВАНИЯ / ВЕРОЯТНОСТЬ ОТКАЗА В ОБСЛУЖИВАНИИ

Аннотация научной статьи по компьютерным и информационным наукам, автор научной работы — Митрофанов Артем Константинович, Петров Максим Александрович, Рудельсон Лев Ефимович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по компьютерным и информационным наукам , автор научной работы — Митрофанов Артем Константинович, Петров Максим Александрович, Рудельсон Лев Ефимович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Model of priority service of aerodrome operations is analyzed. The assumption of the legality of replacing the random variable queue length at its most probable value is used. The model is characterized in that the take-off and landing are described not by means of only safety priorities, but by technical and flight characteristics as well. Evaluations of the service waiting time of service and the probability of failure (or re-service) were obtained.ICAO experts have put forward the concept of "System Wide Information Management" and "Collaborative Air Traffic Management", which require for their introduction of new instruments operational management decision making for air traffic services. It should provide participants with the joint organization of new software for rapid analysis of computer results. One of the possible directions of such support is discussed in this article.Analysis of the two priority model shows that the efficiency of priority service is reduced in the case of a large spread of values of average request processing time. As a result of the loss probability for both types of requests become comparable or even approximately equal.Experiment with the five priority flows confirms the order of request losses in accordance with an increase of flow index, identified by formulas for an arbitrary number of input streams. The tendency of the values for the partial converging efficiency criterion of service can only be observed in non-real case where the values of the average service time varies several hundred times. Moreover, this effect can easily be compensated by a slight increase in volume of storage buffer. It becomes obvious that the best conditions for the maintenance of the flow of requests are provided in the case of approximate equality of the average processing times; that is, when all the coefficients of the sequence tend to a single value.The obtained formulas have a sufficient degree of generality and can be used to analyze the priority system in a wide range of data acquisition and processing tasks.

Текст научной работы на тему «Альтернативная модель приоритетного обслуживания для оценки пропускной способности аэродрома»

Civil Avition High TECHNOLOGIES

Vol. 19, № 04, 2016

УДК 629.735.015:681.3

АЛЬТЕРНАТИВНАЯ МОДЕЛЬ ПРИОРИТЕТНОГО ОБСЛУЖИВАНИЯ ДЛЯ ОЦЕНКИ ПРОПУСКНОЙ СПОСОБНОСТИ АЭРОДРОМА

А.К. МИТРОФАНОВ, М.А. ПЕТРОВ, Л.Е. РУДЕЛЬСОН

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

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

ВВЕДЕНИЕ

Проблемными элементами авиатранспортной системы остаются аэродромы как конечные и отправные пункты всех полетов. Недостатки нередко связывают с качеством организации воздушного движения [1]. Район аэродрома - сложный элемент [2] структуры воздушного пространства (ВП), здесь одновременно набирают высоту вылетающие воздушные суда (ВС), снижаются приземляющиеся и горизонтально движутся транзитные рейсы. В [3] решена задача исследования аэродромных операций с учетом приоритетов рейсов. Получены выражения для оценки вероятностных характеристик работы системы. Уточнен термин «отказ в обслуживании» в приложении к аэродромным операциям. Это не прекращение деятельности, а постановка в очередь на дообслуживание, т. е. либо задержка вылета, либо кружение в зоне ожидания, либо уход на второй круг при посадке и т. д. В такой постановке отказ - это выброс, нарушение стационарного процесса. Каждый борт должен быть обслужен, но в реальной системе стопроцентная вероятность недостижима. Необслуженных рейсов практически не бывает, но если система организована плохо, то она сама себя перегружает потоком перенаправления ВС на запасной аэродром, повторных попыток приземления, ожидания вылета и посадки. Задача состоит в такой организации аэродромных операций, чтобы вероятности дообслуживания стремились к минимуму.

В данной статье в развитие упомянутых результатов к анализу аэродромных операций привлекается учет различий обслуживаемых ВС по их летно-техническим характеристикам (ЛТХ), таким как нормативное время выполнения операции, категория спутной турбулентности, скороподъемность и т. д. Соответственно, общая очередь ВС, рассмотренная в [3], представляется раздельными «секциями», что создает возможности формализации, не использованные ранее.

ПОСТАНОВКА И ФОРМАЛИЗАЦИЯ ЗАДАЧИ

Используем, следуя [3, 4], допущение о правомерности замены случайной величины длины очереди ее наиболее вероятным значением Ь и рассмотрим модель процесса появления и удаления клиентов (заявок) в системе. Специфика обслуживания состоит в том, что вероятность потери заявки произвольного к-го потока, к = 1, т, определяется не динамически изменяющимся остатком буферного накопителя (БН), свободным от заявок /-го приоритета, г = 1, к — 1,

Vol. 19, № 04, 2016

Civil Avition High TECHNOLOGIES

а фиксированным размером собственной секции (очереди). Тогда резкие изменения поведения (случайные выбросы [5]) системы массового обслуживания (СМО) определяются таким соотношением параметров, при которых за сеанс обслуживания одной заявки у'-го приоритета, у = к, т, в к-й секции разделенного БН образуется очередь, наиболее вероятная длина Ьк которой превышает количество гк имеющихся в ней мест для ожидания (рис. 1).

Пусть на вход СМО поступают два простейших входных потока клиентов (заявок): прилетающие и вылетающие ВС, образуя две раздельные очереди в зоне ожидания посадки и старта. Для каждого входящего потока выделены свои секции объемами г] для заявок высшего приоритета (ЗВП) и г2 для заявок низшего приоритета (ЗНП). Потоки - пуассоновские с интенсивностями Х] и Х2 соответственно. Времена обслуживания - экспоненциальные с параметрами ¡и] и /л2, загрузка = р] + р2 < 1. В силу последнего условия р] < 1 и, следовательно, вероятность освобождения системы от ЗВП, при котором на обслуживание назначается ЗНП, конечна. За время обработки, среднее значение Т2 составляет Т2 = 1/и2, в секции первого БН образуется очередь, наиболее вероятная длина которой Q¡ = (1 + }2) р]у, где } = Г /Т - коэффициент вариации (отношение среднеквадратического отклонения Г к среднему значению времени ожидания Т), у = /л] /¡л2. Введение обозначения Qi длины очереди (в отличие от символа Ьк) продиктовано стремлением формально подчеркнуть, что при переходе от потока к потоку она не суммируется, как это считалось в модели с общим БН [3], а определяется отдельно по каждому потоку. Если Q1 = (1 + })р1у< г1, т. е. накапливающаяся очередь может быть размещена

в соответствующей секции БН, то, согласно введенному выше допущению, на характеристики обслуживания ЗВП второй поток не оказывает существенного влияния. Следовательно, при выполнении неравенства Q1 < г1 вероятность щ потери заявки первого типа может быть приближенно оценена как

Рис. 1. СМО с относительными приоритетами и приемом заявок в раздельные секции очереди

П =р +1(Р+Р), если (1 + })ру<V.

Действие формулы справедливо слева от точки излома укр = г] / (1 + }2) р]. Для исследования поведения функции ж] = /(7) справа от укр рассчитаем по аналогии с [3] долю потерь ЗВП за время обслуживания одной ЗНП: <5 = 1 - г] / (1 + }2) р]у, а также вероятность потери ЗВП, пропорциональную загрузке р2 системы потоком ЗНП и вероятности Р2 обслуживания ЗНП в условиях приоритетности ЗВП, являющейся дополнением до единицы вероятности п2 потери ЗНП. На этом аналогия с моделью СМО с общим БН исчерпывается. Вероятность п2 потери ЗНП зависит от соотношения длины Q2 очереди, образующейся во второй секции, с количеством мест ожидания в ней, а не определяется отсутствием очереди ЗНП справа от укр. Накопление ЗНП происходит во время занятости СМО обслуживанием ЗВП. Рассчитаем наиболее вероятное значение длительности этого времени. Длину q¡ очереди ЗВП можно оценить как

1(1 + }2)рху< г1, если условие (0; < г1) выполняется, ^ = 1

[V в противном случае.

Civil Avition High TECHNOLOGIES

Vol. 19, № 04, 2016

Очередь q1 образуется в единичном сеансе обслуживания ЗНП. Наиболее вероятное значение времени освобождения системы от qi ЗВП можно представить как Тч = (1 + За этот период в секции БН, отведенной для первого потока, могут поступить и быть принятыми еще q1'заявок, причем q1' = (1 + $2)Х1ТЧ = (1 + 'д2^1Х1Т1 = (1 + $2) q1 р1. Заметим, что и за время обслуживания очереди длиной q1' в системе накапливается новая порция q1" ЗВП:

= (1 + = (1 + Я) Р12.

Вообще за период занятости системы ЗВП с учетом убывающего потока приращений q1', q1",..., qP будет обслужена очередь, наиболее вероятное значение q1* длины которой складывается из частичных сумм q1(i): q1* = (1 + Я2) (1 + р1 + р12 + ... + р/ + ... + р/) q1.

Последовательность |р/} представляет собой геометрическую прогрессию со знаменателем р1 < 1, а величина q1 в исследуемой области (справа от укр) равна г1. Следовательно,

. = (1 + Я2)(1 — А+1) 1 — А

Показатель степени при р11 определяется из тех соображений, что минимальная длина очереди (1 + Я2)р < 1; ] I [ = — 1п[(1 + г?2)г;]/1п р, способной загрузить систему, не превосходит единицы. Тогда наиболее вероятное значение длительности Т периода занятости системы обслуживанием очереди q1* ЗВП с учетом убывающих приращений ее длины составляет

Т = (1 + Я) д* Т = (1 + Я)(1 — Р1+1) г1Т1/(1 — р).

За это время в секции БН объемом г2 мест для ожидания, отведенной ЗНП, накопится очередь, наиболее вероятная длина Q2 которой составляет

Q2 = (1 + Я) Х2 Т = (1 + Я)(1 — Р1+1) гТ^К- А) = а + Я)(1 — Р1+1) Ру(1 — р,).

Очевидно, что если Q2 > г2, то за время освобождения СМО от q1* ЗВП вторая секция БН будет переполняться, что повлечет за собой дополнительные потери, доля которых равна

Р = 1 — г2 У(1 — Р1)

Ь2 1 , . _п2\/1 ,;+1> '

r (1+^2)(i-рГ1) р

2

При Q2 < г2, потери ЗНП определяются объемом г2 собственной секции БН и суммарной загрузкой системы р^, так как их обслуживание возможно после освобождения СМО от ЗВП

,2 = АК1"А), если(1+Я2)(1—а1+1)а <г2.

2 1 — РГ1 Г(1 — А) 2

При этом вероятность Р2 обслуживания ЗНП составляет

Р2 = 1 — = 1 — Р*(1 — А) = -1—А+,, если (1 + Я)(1 — А+')гР < г2.

2 2 1 — РГ1 1 — РГ1 7(1 — А) 2

В случае Q2 > г2, когда каждый сеанс обслуживания ЗНП приводит к переполнению обеих секций, формула для вероятности Р2 ее обслуживания учитывает вызываемые этим сеансом

Том 19, № 04, 2016_Научный Вестник МГТУ ГА

Vol. 19, № 04, 2016 Civil Avition High TECHNOLOGIES

дополнительные потери. Она рассчитывается как пересечение вероятностей нахождения в БН хотя бы одного места для ожидания и доли потерь ЗНП вследствие его переполнения при

1-Р2 Г Г(1-А) 1-Р2+1 г (1+^)0-/^)р Получив вспомогательные величины, определяющие характеристики обслуживания ЗВП, т. е. долю потерь вследствие переполнения первой секции в сеансе обслуживания одной ЗНП и вероятности событий, при которых могут происходить потери ЗВП, можно записать

1 — /Л, / /11 — // I

каждом обслуживании ЗНП Р2 = (1-п2)(1-£) = -—-^—-> Г.

п =

р+1(1 -р.)

1 -рГ

рр1+.(1 -р.)

если (1 + $2)р1у< г.;

1 -р

r+2

+ р2

1-

(1+ &2)РГ

1 -р?

1 -р?"

л „2 ч (1 +^2)(1 -р1+1) гхрг s если (1 + $1)р1у> r и--——у < г.

р+'(1 -р)

1 -р+

1-

Y(1 -рг)

1 -р^Т2 Г2

(1 + & )р1 у

Г(1 -р)

1 -р+1 Г1 (1 + &)(1 -р1+1) гр

если (1 + 1&)р1у> г и

(1 + &2)(1 -р1+1) гр Y(1 -р1)

> Л.

Отметим характерное отличие полученного выражения от формулы для дисциплины с общим БН. Вероятность потери ЗВП в зависимости от соотношения параметров теперь имеет либо одну, либо две точки излома: обязательную при укр = г1/(1+^2)р1, когда переполняется собственная зона и значение щ резко возрастает, и другую, когда переполняется секция ЗНП и нарастание щ замедляется - в точке укр2 = (1 + ^2)(1 - Р+1) г1 р2/г2 (1 - р).

Приближенное выражение для вероятности потери ЗНП фактически уже получено при анализе возможных аспектов обслуживания ЗВП:

п р^(1 -р^) + О 1 -р? П =-:--+ О-

1 -рГ

1 -рГ

1 - /2

Y(1 -р)

qx (1 + &2)(1 -р[+1) рг

где

8 = [0, если (1 + ^2)(1 - р[+1)трг/7(1 - р) < Г2, = |(1 + )рУ, если (1 + )рху< гх, [1 в противном случае, 1 в противном случае.

Приведенные выражения подчеркивают зависимость вероятностей жг потери заявки от соотношения у параметров обслуживания различных потоков в системах с относительными приоритетами. Обосновать ее можно достаточно простыми рассуждениями. С возрастанием у при неизменности других параметров снижается интенсивность потока ЗНП и, как следствие, вероятность образования очереди, превышающей размер второй секции БН обратно пропорциональна у). При этом период занятости системы ЗВП фиксируется сверху объемом г1 выделенной им секции БН. В этих условиях вероятность потери ЗНП уменьшается, и с ростом у повышается частота обслуживания таких заявок. В итоге увеличивается вероятность события, при котором теряются ЗВП, что и приводит к резкому возрастанию вероятности щ их потери. Скорость возрастания щ зависит от условий обслуживания ЗНП; если в процессе освобождения системы от ЗВП вторая секция БН не переполняется, то вероятность л2 потерь ЗНП незна-

Civil Avition High TECHNOLOGIES

Vol. 19, № 04, 2016

чительна и щ возрастает резко; в противном случае рост щ ограничивается фиксированным объемом г2 буферной зоны потока ЗНП. При этом начинает играть роль соотношение г2/г1 объемов секций БН, аналитически выведенное выше и установленное эмпирически в известных работах [6].

Расчеты по формулам подтверждаются результатами моделирования СМО с одинаковыми п, т и р/. На рис. 2 представлены графики для системы с параметрами г1 = 2; г2 = 8; на рис. 3 - с одинаковыми размерами секций БН г1 = г2 = 5; на рис. 4 - с г1 = 8 и г2 = 2. Нетрудно видеть, как с увеличением соотношения у интенсивностей обслуживания снижается эффективность приоритетных систем обсуждаемого класса. Вероятности потери заявок различных потоков не только выравниваются, как это было в модели с общим БН, но обнаруживается недопустимая тенденция: ЗВП теряются чаще ЗНП. Выведенная ниже формула для произвольного числа потоков показывает богатую палитру возможных деформаций отношений предпочтения. Так, для трехприоритетной системы в наихудших условиях как по критерию л, так и с точки зрения ожидания в очереди оказываются заявки второго потока.

ОБОБЩЕНИЕ МОДЕЛИ НА СЛУЧАЙ ПРОИЗВОЛЬНОГО КОЛИЧЕСТВА ПОТОКОВ

Распространим полученный результат на произвольное количество входящих потоков. Пусть имеем одноканальную СМО (п = 1, т. е. на аэродроме одна полоса), обслуживающую с относительным приоритетом т пуассоновских входных потоков с интенсивностями Х/ каждый (/ = 1, т ). Для наглядности сначала будем считать, что назначение приоритетов выполнено в соответствии с возрастанием величин у/. В дальнейшем это ограничение снимается за счет введения последовательности символов ¿¡к. Времена обслуживания заявок всех потоков независимы и распределены экспоненциально с параметрами причем у/ = / ¡и/. Каждому /-му приоритетному потоку выделена собственная зона БН объемом на г/ мест для ожидания. Рассчитаем вероятность лк потери заявки произвольного к-го потока в функции параметров системы.

5 10 20

Рис. 2. Параметры: п = 1; т = 2; р/ = 0,45 и г1 = 2 < г2 = 8

50 100

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

5 10 20

Рис. 3. Параметры: п = 1; т = 2; р/ = 0,45 и г = г2 = 5

50 100

5 10 20

50 100

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Рис. 4. Параметры: п = 1;

т = 2; р/ = 0,45 и г = 8 > г2 = 2

0.4

0.8

0.3

0.3

0.6

0.2

0.2

0.4

0.1

0.1

0.2

I 2

I 2

I 2

Vol. 19, № 04, 2016

Civil Avition High TECHNOLOGIES

Условие излома кривой лк, определяемого процессом заполнения к-й секции БН, принимает вид Qk = (1+}2) Хк Т, где Т - время освобождения СМО от заявок более высоких, чем к-й, приоритетов. По аналогии с двухприоритетной моделью нетрудно заключить, что общее количество ^ изломов кривой лк = /(у) может колебаться от одного до ^тах = т - к + 1, т. е. зависит не только от условий обслуживания заявок низших, чем к-й, приоритетов (с индексами у > к), но и от места к-го потока в шкале приоритетов.

Для определения величины времени ожидания Т рассмотрим процесс образования оче-

к-1

реди Ьк-1 заявок более высоких, чем к, приоритетов, исходя из соотношения Ьк-1 = ^ Qi, где

Q =

(1 + ^2)Гк

' р >

если Q < гг,

V Y J

ri в противном случае.

Наиболее вероятное значение TQi занятости системы обслуживанием очереди Qi заявок составит TQi = (1+}2) Qi Ти где Т. - среднее время обслуживания заявки /-го типа, и для очереди

к-1

к-1

Ьк-1 можно приближенно записать ТЬ = ^ Т& = (1 + }2) ук ^ (( pi / ).

г'=1

i=1

В течение этого времени в системе образуется приращение АЬ' очереди заявок потоков

т

высших приоритетов с индексами I < к. В силу стационарности и условия ^ р < 1 интенсивно-

i=\

сти Х. поступлений таких заявок меньше параметров /л. их обслуживания и суммарная очередь Ьк-1 вместе с приращениями АЬ', АЬ", ..., АЬ1' со временем укорачивается:

k-1 f Pi ^ v Y i у

i=1

k -1 f P. ^ 2

AL = (1 + Yk ^^ < 1. AL' = (1 + d2)Yk , ..., AL" = (1 + tf2)Yk £

V Y - у

k-1 ' P V v Yi j

Вообще, при расчете Ьк-1 можно рассмотреть ^ монотонно убывающих приращений суммарной очереди АЬ— заявок /-х уровней приоритетности, более высоких, чем к-й, (! < к). Соответственно, процедура вычисления Ьк-1 заканчивается на I-м шаге, когда (как и в модели с общим БН) получаем приращение очереди, не превышающее единицы. Тогда считаем, что система полностью освободилась от ЗВП, принадлежащих потокам с индексами I < к, и может приступить к обслуживанию заявок к-го типа. Так как поступление заявок разных типов пропорционально ин-тенсивностям входящих потоков, для выполнения неравенства Ьк-1 < 1 достаточно потребовать,

чтобы (1 + }2)Гк (р/К)) Я , откуда ](. [ = 1п Я - 1п £ Я^ 1п

(1 + ^2)Yk р Yi

Тогда наиболее вероятная длина Ь очереди всех заявок с индексами ! < к, образующейся в системе за время обслуживания одной заявки у'-го типа, у = к, т, с учетом I приращений ее

к-1

длины при освобождении от Ьк-1 ЗВП, составит Ь = (1 + }2) ук ^

Г Л1'+

P

V Yi J

« р ^ 1

Yi

, а пери-

i У

од Т занятости системы обслуживанием накопленной очереди Ь заявок -

T = (1 + tf2)Yk IT

1 -

/ nA +1

Yi

1-

Pi

(1)

1=1

=1

i=1

Y

i=1

Civil Avition High TECHNOLOGIES Vol. 19, № 04, 2016

За это время в k-й секции разделенного БН образуется очередь заявок к-го типа, наиболее

к-1

вероятная длина Qk которой равна Qk = (1 + tf2)рк ^ уг

( V'+1 Р

V Y У

« р ^ 1

Y

' У

Очередь Qk определяет условие обязательного излома (по переполнению собственной

к-1

секции БН) кривой щ = /ОХ Qk = (1 + tf2) рк £ у

'=1

sh +1

1 -

р. 1 -р

V Y у V Y у

При выполнении последнего неравенства потери заявок к-го типа определяются суммарной загрузкой системы первыми к потоками и объемами гк собственных секций БН:

\kf

\

\rk +

Пк = IX Р 1 - I Р 1 - IX Р

V '=1 у V ¿=1 Л V ¿=1 У

если Qk < Гк. Справа от точки излома потери возрастают за счет переполнения к-й секции БН в сеансе обслуживания одной заявки у-го типа, ] = к, т, причем доля £к потерь составляет

4 = = 1 - rj (1+ tf2) рк Ц y

1 -

Г Л1' +1 р

V Y У

р

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

\

1 -

V У' у

Вероятность Рк события, при котором за время Ту обработки одной заявки]-го типа (] > к) теряются & заявок к-го типа, пропорциональна загрузке Р] потоком ]-го приоритета и вероятности обслуживания заявки к-го типа при условии возможного прореживания]-го потока:

( k \ rk 1 - (z р)

Pk = (1 -п)(1 -4)=—'

k-1

1 - II р

V '=1

\ rk +1 Tk!

(1 + tf2) рк ^Y

Л +1

д 1 -Р

V Y у V Y у

Прореживание]-го потока имеет место при нарушении сходного с (2) условия Qj < г; если же оно выполняется, то Рк = 1 - пк. Окончательная форма записи приближенного выражения для оценки вероятности пк потери заявки произвольного к-го потока принимает вид составного выражения 2:

П

" к У'к ( к \ z р 1 - z р

ч '=1 У V ' =1 У ( к \rk +1 1 - (z р

V '=1 У

( к Yk f к \

Zр 1 -Zр| m 1 1 + Р

j=1

рк

п = __^L

П f к Vk +1

если для всех j = к +1, m справедливо (1 + tf ) Y ^ гк;

Yk

rk Yk

1 HI р

. '=1

к-1

(1+tf2) р Zy

1 -

f р

I Г,)

1' +1

(1 + tf2) Yj рк

1

i=1

=1

Vol. 19, № 04, 2016

Civil Avition High TECHNOLOGIES

где

S

jk

0, если (1 + } V — < г.,

П

1 в противном случае,

(2)

2 рк

если хотя бы для одного у = к + 1,т, (1 + } )у, — > гк и очередь

k-1

Qk = (1 + ^2) Pk Sy

i=1

1 -

f P

IY J

li +1'

' k Yk / k Л r k Yk

S P 1 - S A I и 1 - Is P

-+SS Pj v 1=1 J j=1

rk -Yk

(1 + ^2)• Yj •Pk

если хотя бы для одного у, ] = к +1, т, (1 + г2) ~ > гк и 0к >

rk •

Наиболее вероятное время Т ожидания аэродромной операции рассчитывается по формуле (1), утверждающей, что до взятия на обслуживание заявки к-го типа СМО должна освободиться от всех заявок более высокого приоритета, чем к-й. Обратим внимание на неизменность степени множителя (1 + т}2 ) на каждом шаге расчета приращений очередей ЗВП при ее выводе. Приращения АЬ', АЬ", ..., АЬ1' суммарной очереди Ьк-1, определяются как средние -

Y k S

k-1 f pi V

1=1

v Y i у

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

вероятной величине всякии раз происходит заново и не зависит от количества предшествовавших итераций.

Выражение (2) дает возможность приближенно рассчитывать вероятности потери заявок в системах с относительными приоритетами и приемом заявок в раздельные секции БН при известных значениях р, у, г' по каждому входящему потоку. Результаты вычислений по формулам подтверждаются статистическим моделированием. Тем не менее необходимо еще раз подчеркнуть, что допущение о правомерности замены случайной величины длины очереди ее наиболее вероятным значением, на котором основан их вывод, и, следовательно, методы расчета объемов БН, вносит в результаты дополнительную погрешность. Для СМО с параметрами п = 1; т = 2; г = 10; р1 = р2 = 0,45 ошибка относительно результатов моделирования достигает 15 % в окрестности точки излома у = г / (1 +}2) р1 и менее 10 % - в остальном диапазоне изменения у.

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

Civil Avition High TECHNOLOGIES

Vol. 19, № 04, 2016

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

ЗАКЛЮЧЕНИЕ

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

СПИСОК ЛИТЕРАТУРЫ

1. Этап определения проекта SESAR (Исследование организации воздушного движения в едином Европейском небе). Документ Евроконтроля DLM-0607-001-02-00, 2007.

2. Федеральные авиационные правила «Организация воздушного движения в РФ». Утв. Минтрансом РФ 25.11.2011, № 293. - Официальный текст. - Режим доступа: http://dvmtu-favt.ru/about/documents/18/200.

3. Митрофанов А.К., Петров М.А., Рудельсон Л.Е. Модель приоритетного обслуживания аэродромных операций. // Научный вестник МГТУ ГА. 2015. № 221.

4. Башарин Г.П. Об однолинейной системе с эрланговским распределением времени обслуживания для заявок различных видов // Известия АН СССР. Техническая кибернетика. 1966. № 3.

5. Кузьмин В.П. Метод оценки вероятностей больших выбросов случайных процессов // Ученые записки ЦАГИ. Вып. 3-4. Том XXXIV. 2003.

6. Таненбаум Э., Вудхалл А. Операционные системы. Разработка и реализация. СПб.: Питер, 2007.

7. Руководство по концепции общесистемного управления информацией. Документ ИКАО Doc 10039 AN/511. Монреаль, 2015.

8. Руководство по совместной организации потоков воздушного движения. Документ ИКАО Doc 9971 AN/485. 2-е изд. Монреаль, 2014.

ALTERNATIVE MODEL FOR EVALUATING AERODROME CAPACITY BASED ON PRIORITIZED QUEUE SERVICE

Mitrofanov A.K., Petrov M.A., Rudel'son L.E.

Model of priority service of aerodrome operations is analyzed. The assumption of the legality of replacing the random variable queue length at its most probable value is used. The model is characterized in that the take-off and landing are described not by means of only safety priorities, but by technical and flight characteristics as well. Evaluations of the service waiting time of service and the probability of failure (or re-service) were obtained.

Vol. 19, № 04, 2016

Civil Avition High TECHNOLOGIES

ICAO experts have put forward the concept of "System Wide Information Management" and "Collaborative Air Traffic Management", which require for their introduction of new instruments operational management decision making for air traffic services. It should provide participants with the joint organization of new software for rapid analysis of computer results. One of the possible directions of such support is discussed in this article.

Analysis of the two priority model shows that the efficiency of priority service is reduced in the case of a large spread of values of average request processing time. As a result of the loss probability for both types of requests become comparable or even approximately equal.

Experiment with the five priority flows confirms the order of request losses in accordance with an increase of flow index, identified by formulas for an arbitrary number of input streams. The tendency of the values for the partial converging efficiency criterion of service can only be observed in non-real case where the values of the average service time varies several hundred times. Moreover, this effect can easily be compensated by a slight increase in volume of storage buffer. It becomes obvious that the best conditions for the maintenance of the flow of requests are provided in the case of approximate equality of the average processing times; that is, when all the coefficients of the sequence tend to a single value.

The obtained formulas have a sufficient degree of generality and can be used to analyze the priority system in a wide range of data acquisition and processing tasks.

Key words: air traffic management, airdrome operations, priority service, the service waiting time of service, the probability of a denial of service.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

REFERENCES

1. The ATM. Target. Concept. D3. DLM-0607-001-02-00, 2007. SESAR Definition Phase -Deliverable 3. Issued by the EUROCONTROL SESAR Consortium or the SESAR Definition Phase Project - September 2007.

2. Federal Aviation Regulations "Air Traffic Management in the Russian Federation." Approved. RF Ministry of Transport 25.11.2011, № 293. The official text. http://dvmtu-favt.ru/about/documents/18/200 Access (in Russian).

3. Mitrofanov A.K., Petrov M.A., Rudel'son L.E. Model of priority service for airdrome's operations. The MSTUCA Scientific Bulletin. 2015. No. 221 (in Russian).

4. Basharin G.P. About certain single-line system with Erlang distribution for service time onto different types of requests. Proceedings of the Academy of Sciences of the USSR, Technical Cybernetics. 1966. № 3 (in Russian).

5. Kuz'min V.P. Method of estimation of the probability of large emissions of random processes. Scientific notes TsAGI. Vol. 3-4. That XXXIV. 2003. (in Russian).

6. Tanenbaum E, Woodhull A. Operating Systems. Development and implementation. SPb: Peter, 2007 (in Russian).

7. Manual on System Wide Information Management (SWIM) Concept. Doc. 10039 AN/511, ICAO, 2012.

8. Manual on Collaborative Air Traffic Flow Management, Doc. 9971, ICAO, Montreal, Canada, 2012.

СВЕДЕНИЯ ОБ АВТОРАХ

Митрофанов Артем Константинович, диспетчер-стажер Московского центра автоматизации управления воздушным движением, аспирант МГТУ ГА, электронный адрес: [email protected].

Петров Максим Александрович, аспирант МГТУ ГА, электронный адрес: [email protected].

Рудельсон Лев Ефимович, д.т.н., профессор МГТУ ГА, электронный адрес: [email protected].

i Надоели баннеры? Вы всегда можете отключить рекламу.