Научная статья на тему 'Модель и метод расчета системы обслуживания m/ h 2/ n- h 2 с ограничением времени пребывания'

Модель и метод расчета системы обслуживания m/ h 2/ n- h 2 с ограничением времени пребывания Текст научной статьи по специальности «Математика»

CC BY
414
82
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
СИСТЕМЫ МАССОВОГО ОБСЛУЖИВАНИЯ / QUEUING SYSTEMS / НЕТЕРПЕЛИВЫЕ ЗАЯВКИ / IMPATIENT CUSTOMERS / ОГРАНИЧЕННОЕ ПРЕБЫВАНИЕ / ГИПЕРЭКСПОНЕНЦИАЛЬНОЕ РАСПРЕДЕЛЕНИЕ / HYPEREXPONENTIAL DISTRIBUTION / ЧИСЛЕННЫЕ МЕТОДЫ / NUMERICAL METHODS / RESTRICTED SOJOURN TIME

Аннотация научной статьи по математике, автор научной работы — Уланов Александр Викторович

Предложена модель системы массового обслуживания с ограничением времени пребывания заявок M / Н 2/ n Н 2 в виде марковской цепи с непрерывным временем и метод расчета стационарного распределения и вероятности потери заявки в данной системе. Результаты верифицированы с помощью имитационной модели.

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

Похожие темы научных работ по математике , автор научной работы — Уланов Александр Викторович

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

Текст научной работы на тему «Модель и метод расчета системы обслуживания m/ h 2/ n- h 2 с ограничением времени пребывания»

УДК 519.872

А.В. Уланов

Научный руководитель: заслуженный деятель науки РФ, доктор технических наук, профессор Ю.И. Рыжиков

МОДЕЛЬ И МЕТОД РАСЧЕТА СИСТЕМЫ ОБСЛУЖИВАНИЯ MIH2ln-H2 С ОГРАНИЧЕНИЕМ ВРЕМЕНИ ПРЕБЫВАНИЯ

Предложена модель системы массового обслуживания с ограничением времени пребывания заявок М/Н2/п-Н2 в виде марковской цепи с непрерывным временем и метод расчета стационарного распределения и вероятности потери заявки в данной системе. Результаты верифицированы с помощью имитационной модели.

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

The paper presents the method of calculating the multi-channel queuing systems with impatient customers with Poisson arrival flow and second order hyperexponential distributions both service time and restricted sojourn time М/Н2/п-Н2. The results are verified using the simulation model.

Queuing systems, impatient customers, restricted sojourn time, hyperexponential distribution, numerical methods.

Введение. Во многих практических задачах теории очередей возникают ситуации, когда время пребывания заявки (время «терпения») в системе массового обслуживания (СМО) ограничено некоторой (в общем случае случайной) величиной. Примеры подобных ситуаций:

- в здравоохранении - при организации медицинской помощи на каком-либо ее этапе лечение больного может закончиться летальным исходом;

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

- в логистике - срок годности товаров при перемещении по цепям поставок ограничен;

- в управлении - информация со временем теряет актуальность и т.д.

Обстоятельный обзор современного состояния вопроса приводится, например, в [7]. Приходится констатировать, однако, что подавляющее большинство известных методов относится либо к чисто марковским системам М1М1п-М [2] (в дополнении к нотации Кендалла через дефис указан тип распределения времени терпения), либо ограничены весьма жесткими условиями - одноканальное обслуживание и большая загрузка [7]. Из работ по анализу многоканальных немарковских систем стоит выделить статью [8], в которой рассматривается диаграмма переходов системы с марковским обслуживанием и Н2-распределением терпения М/М/п-Н2, а также работу [3], в которой предложен метод расчета системы с М4Р-входящим потоком, матрично-фазовым (РИ) распределением обслуживания и марковским терпением. При этом не рассматривались системы с ограничением на полное (ожидание плюс обслуживание) время пребывания заявки в системе.

С точки зрения практической значимости наибольший интерес представляют многоканальные системы с простейшим входящим потоком и произвольными (немарковскими) распределениями времен

обслуживания и терпения. В данной статье представлен метод расчета многоканальной системы с простейшим входящим потоком и Н2-распреде-лениями обслуживания и предельного пребывания М/Н2/п-Н2. Применение Н2-распределения позволяет выровнять три начальных момента любого исходного (кроме Эрланга второго порядка), что обеспечивает приемлемую точность при расчете СМО [6].

Расчет системы М/Н2/п-Н2 проходит следующие этапы:

- построение диаграммы переходов;

- формирование в соответствии с диаграммой матриц интенсивностей переходов;

- составление уравнений баланса переходов, их решение и расчет стационарных характеристик СМО.

Остановимся подробнее на каждом этапе.

Диаграмма переходов. Рассмотрим фазовое представление состояний системы М/Н2/п-Н2 и переходов между ними. Пусть Н2-распределение обслуживания имеет параметры {у,, ц}, где {у,}, , = 1,2 - вероятности выбора фаз, {ц}, - = 1,2 - интенсивности показательного распределения времени обслуживания в них. Аналогично зададим параметры Н2-распределение терпения - {и,, у,}. Заметим, что параметры Н2-распределения могут принимать комплексные (попарно-сопряженные) значения, что подчеркивает фиктивный процесс расщепления процесса на фазы. Это допущение позволяет аппроксимировать распределения с коэффициентами вариации, меньшими единицы.

Для каждой вновь прибывшей в систему заявки необходимо зафиксировать ее тип по времени обслуживания (с вероятностью у: она будет первого типа и с вероятностью у2 - второго) и терпения (аналогично с вероятностями и1, и2). Расстановки заявок по типам обслуживания (,, -) и терпения (I, т) образуют микросостояния системы, которые характеризуется четырехразрядными «ключами» микросостояний (,, I, т), , + - = тт(к, п), I + т = к, где к - ко-

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

0-й ярус

1 -и ярус '

о-

00

о-

10

о-

01

10

2-й ярус

о-

20

-о-

-о-02

о-

20

о-

-о-

-о-

-о-02

-о-

20

11

02

'20

И

02

Рис. 1. Многослойное представление состояний системы М/Н2/п-Н2

Номер яруса соответствует количеству заявок в системе, двузначные коды по горизонтали - расстановку по типам обслуживания, справа по вертикали - по типам терпения. Начиная с яруса ] > п, прибывшие заявки становятся в очередь и не влияют на обслуживание; соответственно, диаграмма перестает расти «вширь», но увеличивается количество слоев каждого нижележащего яруса.

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

Изображение всех возможных переходов между микросостояниями на данной диаграмме приводит к полной потере наглядности. Поэтому ограничимся переходами из микросостояния с ключом (/, ], I, т) в смежные с ним. На рис. 2 показаны упомянутые переходы в «недогруженной» системе, т.е. когда имеются свободные каналы обслуживания. Интенсивность перехода из одного микросостояния в другое находится как произведение коэффициентов при соединяющих их стрелках.

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

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

Рис. 2. Переходы в недогруженной системе М/Н2/п-Н2

Уг

1/к /-1 /-1 /+1 т Уг

т/к

. /-1

— т

Уг

/+1 }~1

/ + 1 т- 1

1/к / /

/ т-1

О!

т к

/-1 т

Ик

/

-1

т к

У1

У1

У1

У1

]/и г/и

1-и/к и к

г>1

/ / 1

т /

туг

1-и/к

г/и

]/п

и/к

У1

Уг

«1

иг

. /+1

/ т+1

/

т+1

Рис. 3. Переходы в загруженной системеМ/Н2/п-Н2

При уходе заявки по завершению обслуживания следующее микросостояние выбирается исходя из того, какой она имела тип терпения (вероятности пропорциональны числу заявок соответствующего типа), и какого типа обслуживания требует головная заявка очереди (с вероятностью у1 первого типа, с вероятностью у2 - второго). При уходе по исчерпанию терпения для начала устанавливается, из канала он произошел или из очереди (вероятность каждого из этих событий пропорциональна числу каналов и длине очереди соответственно). Если «нетерпение» случилось в очереди, то уменьшается количество заявок соответствующего типа терпения; если в канале - выбирается какого типа по обслуживанию была ушедшая заявка, затем - какого типа обслуживания требует головная заявка очереди.

Матрицы переходов и уравнения баланса. На основе диаграмм (рис. 2, 3) могут быть получены диаграммы интенсивностей переходов в микросостояния (/, /, /, т), необходимые для записи уравнений баланса вероятностей состояний.

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

Обозначим через / множество всех возможных микросостояний обслуживания на /-м слое /-го яруса

системы. Номер слоя совпадает с количеством заявок второго типа терпения, / = 0, / . Через с/ обозначим количество микросостояний в В случае Н2-распределения с/ = шт(п, /) + 1. По прибытию заявок каждый «крайний» слой-приемник /-го яруса (нулевой и /-й) связан переходами лишь с одним крайним слоем-источником - нулевым и (/-1)-м соответственно - (/-1)-го яруса. Остальные его слои-

приемники / = 1, / — 1 связаны с двумя слоями-источниками вышележащего яруса. По уходам заявок каждый слой-приемник/-го яруса связан с двумя слоями-источниками (/+1)-го. В соответствии с диаграммой переходов и этим правилом формируются следующие матрицы интенсивностей инфинитези-мальных переходов:

А/ / х ст/+1], т=0, 1 - в Б/+1,+т (прибытие заявки), где т=0 соответствует первому типу терпения прибывшей заявки, т=1 - второму;

б™.. [ст/ хст/—- в —..—т (уход заявки по обслуживанию или исчерпанию терпения), причем для

слоев , = 1, - -1 возможен переход на два слоя вышележащего яруса, поэтому т = 0, 1, а для крайних слоев при , = 0 т = 0, при , = - т = 1;

в,, [ст х ст ] - ухода из состояний ,-го слоя яруса

- (в квадратных скобках здесь и далее указывается размер матриц).

Верхний индекс у первых двух матриц обозначает, какого типа терпения уходит (или приходит) заявка. Введем векторы-строки я-1 ={я-п, яj■ 2,..., я-, ст }

вероятностей нахождения СМО в состояниях ,-го слоя --го яруса. Теперь можно записать векторно-матричные уравнения баланса переходов между состояниями слоев для ярусов - = 0, Я -1:

Ко,0А.,0 =Ещ,тВтт , - = 0;

т=0 1

К-,г = X (К]-1,г-тА]-1,г-т + К]+1,г+тВ-+1,г+т ),

т=0

, = 0, -, - = 1, Я -1.

(1)

Для удобства записи уравнений будем считать, что у «крайних» слоев-приемников матрицы интен-сивностей несуществующих переходов и векторы-строки вероятностей состояний несуществующих слоев-источников нулевые.

Для последнего Я-го яруса, соответствующего предельной учитываемой длине очереди,

кЯ,Бя, , = X (к--1,,-тАт-г,,-т), , = 0, Я. (2)

Систему (1-2), дополненную условием нормировки, практически приходится расписывать покомпонентно. Даже для моделей с ограниченной очередью она характеризуется чрезвычайно высокой размерностью, и стандартные методы решения систем линейных алгебраических уравнений применительно к ней оказываются малоэффективными. Мы развили применительно к ситуации с нетерпеливыми заявками итерационный метод [5] решения системы (1-2), восходящий к Такахаси и Таками [9]. Из последней работы здесь используется только одна ключевая идея: расчет в терминах условных вероятностей микросостояний.

Положим 1 = я- ¡1 pj, где р- - суммарная вероятность наличия в системе ровно - заявок, и обозначим

= р-+7р-, = р--7р-.

(3)

Тогда систему (1-2) можно переписать относительно векторов условных вероятностей {- ,}, нормированных к единице в пределах яруса:

'0,0Б0,0 =Х Х011тВтт , - = 0;

т=0 1

, г = X (2--1,,-тАт-1,/-т + Х]']+1,/+тВт+1,г+т ),

1 , = 0, -, - = 1, Я -1;

'я,Бя, ,. = 2я XX ('--1,,.-тА1,.-т ), , = 0~Я, - = Я.

т=0

(4)

С помощью векторов-столбцов 1- = {1, 1, ... , 1}т размера с- для всех - могут быть записаны дополняющие систему (4) уравнения нормировки

X'/,,1 - =1 (5)

г = 0

и баланса суммарных интенсивностей переходов между (;-1)-м и--м ярусом

X-1,, X а- ,1 =1Х, X Я1- (6)

, = 0 т =0 V'=0 т=0 /

Перепишем уравнения системы (4) для - = 1, Я -1 в виде

'(к) Б = 2(к)'(к) Ат + Х(к)'(к-1) Вт )

'л,г - V-1,г--1,,-tj+1,г+т2";+1,г+т/'

т=0

г = 07,

где верхний индекс показывает номер итерации. Теперь ясно, что

= 2->Р],, + Х-^,,

(7)

где

Р' = 4"' '(к-1) Ат Б-1-

т=0 1

Р" = X '(к-1) Вт Б-1

Рг ¿^ у+1,г+^ у+1, г+^ у, г'

т=0 _

г = 0, -, - = 1, Я -1.

(8)

Осталось указать способ расчета {2( }} и {х Перепишем (6) с учетом (7):

j-l 1

^5 X Ь-1, г X А/Т-1. , 1 - =

г =0 т=0

= I 2(к> ^ Р7,г X В^г + Х-к> ^ Р"],г ^ В^г, -^^^

г=0 т=0

I=0 т =0

Отсюда следует пропорциональность

2 -) = с/ ),

(9)

т=0

т=0

Х

с коэффициентом

, X В™,, 11,-

1 =0 т=0

1-1 1

X ¿1-1,, X ^^ 1,, 11 -[Х Р1,, X1111-1

1=0 т=0

1=0 т=0

(10)

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

Подстановка (7) и (9) в (5) дает

1 )Х

1=0

1; = 1.

Итак,

1

=1 ХК,,+Р; ,, )11.

(11)

Удобным критерием прекращения итераций является условие

тах } - х(к-1) <е.

После прекращения итераций можно переходить к нахождению абсолютных значений вероятностей. Из определения чисел {х^} следуют равенства

Р< = Х1 -1 Р<-1, 1 = 1, Я.

(12)

Далее принимается р0= 1, а последующие вероятности считаются согласно (12) с одновременным накоплением суммы. Затем все вычисленные вероятности делятся на упомянутую сумму.

При известных л 11 = {л 1 п, лу ,■ 2,..., л 11 п } вероятность ухода заявки по нетерпению

Ршр = Г' X X [О'-1) +'У2 К, 11. (13)

1=1 1=0

Программная реализация. Описанный метод предполагает работу с переменным числом матриц переменной размерности (в зависимости от числа каналов и максимальной длины очереди). Поэтому для хранения матриц А, В и Б при программировании на Фортране 90 выделялась динамическая (а11о-са1аЪ1е) память.

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

источнику соответствует два слоя-приемника. Поскольку 1-й ярус включает в себя 1 + 1 слой, общее количество упомянутых матриц

Жв = 2Х (1 +1) = (Я + 3) Я.

1=1

Предельный размер матриц определяется по матрицам последнего (Я-го) яруса. Количество микросостояний его слоев-источников (строк матрицы переходов) составляет п+1, количество микросостояний слоев-приемников вышележащего (Я - 1)-го яруса (столбцов матрицы переходов) - п + 1. Итоговая характеристика измерений набора матриц В -[(Я + 3)Я, п + 1, п + 1]. По аналогичной методике рассчитывались размеры всех прочих наборов. Наборы векторов условных вероятностей t и диагоналей матриц Б упаковывались в двумерные массивы.

Поскольку параметры Н2-распределения могут принимать комплексные значения, тип данных для матриц А, В, Б, t задавался сотр1ех*16.

Для связи номера матрицы в наборе с номерами ярусов и слоев диаграммы были созданы массивы адресов - трехмерные для А и В и двумерные - для Б и В трехмерных массивах первое измерение соответствовало номеру яруса, с которого уходит заявка, второе и третье - номеру слоя-источника и слоя-приемника соответственно; в двумерных первое измерение - номеру яруса, второе - номеру слоя.

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

Для верификации полученной программы на Фортране 90 была разработана имитационная модель. Особенности ее программирования докладывались на VI Всероссийской конференции по имитационному моделированию [4].

Численные результаты. Было рассчитано стационарное распределение заявок в системе М/Н2/п-Н2 и вероятность потери заявки по нетерпению -формула (13), при следующих исходных данных:

- число каналов п = 3;

- среднее время обслуживания Ь1 = 1;

- среднее время терпения £1 = 2,0;

- коэффициент загрузки р=0,5;

- предельное число заявок в системе Я = 11.

Необходимая для расчета интенсивность входящего потока находилась по формуле 1 = рп/Ь1. Заметим, что коэффициент загрузки не учитывает уход нетерпеливых заявок, поэтому он теряет определяющую роль и условие р < 1 для существования стационарного режима не предъявляется.

По трем начальным моментам рассчитывались параметры гиперэкспоненты для аппроксимации детерминированного Б (коэффициент вариации V = 0) терпения и Эрланга третьего порядка Е3 (V = 0,577) распределения обслуживания.

В табл. 1 приведены результаты расчета стационарного распределения {р1}, 1 = 0,11 числа заявок в системе М/Б/3-Е3 предложенным в статье численным методом (Числ.) и с помощью имитационной модели (ИМ) для 1 млн испытаний. Выбор именно Б и Е3-распределений обусловлен тем, что при их аппрок-

1 =0

симации гиперэкспонентой параметры последней принимают комплексные значения. Преследуется цель показать, что такие парадоксальные значения не влияют на конечный результат.

Отклонение полученных численным методом вероятностей до порядка 10-5 от имитационной модели укладывается в 15 %, что принято считать вполне приемлемым [1]. Это говорит о правильности идеи и программной реализации предложенного расчетного метода. Заметим, что из-за несовершенства датчиков равномерно распределенных случайных чисел, на основе которых моделируются все прочие распределения, имитационная модель не может считаться идеальным эталоном, особенно при вычислении вероятностей редких событий. Дополнительную погрешность в численный метод вносит учет только трех начальных моментов исходных распределений. Также видно, что при заданных исходных данных вероятность длинной очереди имеет порядок 10-610-12, что говорит о сильном влиянии ухода нетерпеливых заявок на загрузку системы.

Рассмотрим вероятности потери «нетерпеливой» заявки (формула 13), которая является интегральной характеристикой СМО, при других распределениях обслуживания и терпения и различных коэффициен-

тах загрузки р (табл. 2). В скобках указаны коэффициенты вариации распределений.

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

Заключение. Перечислим основные результаты данной статьи.

1. Предложен метод расчета систем с ограничением времени пребывания заявок М/О/п-О с помощью Н2-аппроксимации времен обслуживания и терпения.

2. Установлено, что метод позволяет получить корректный результат даже в области комплексных параметров Н2-распределения.

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

Таблица 1

Вероятности р] системы ЫЮ/3-Е3

- ИМ Числ. - ИМ Числ. - ИМ Числ.

0 6,240е-1 6,288е-1 4 1,296е-3 1,173е-3 8 8,802е-8 3,436е-8

1 2,943е-1 2,912е-1 5 1,179е-4 1,042е-4 9 0,000е+0 1,970е-9

2 6,934е-2 6,792е-2 6 9,135е-6 8,018е-6 10 0,000е+0 1,046е-10

3 1,095е-2 1,078е-2 7 6,382е-7 5,507е-7 11 0,000е+0 5,128е-12

Влияние коэффициентов загрузки и вариации распределении на вероятность ухода заявки по «нетерпению»

Распределение терпения

и1 н Р Е3 (0.577) М (1) Н2 (2)

рц ° ИМ Числ. ИМ Числ. ИМ Числ.

Б 0,7 4,95е-2 5,18е-2 9,05е-2 9,04е-2 2,23е-1 2,22е-1

(0) 1,4 3,59е-1 3,72е-1 4,09е-1 4,10е-1 4,58е-1 4,56е-1

М 0,7 7,90е-2 8,18е-2 2,02-1 2,02е-1 2,39-1 2,38е-1

(1) 1,4 3,69е-1 3,77е-1 4,20е-1 4,20е-1 4,74е-1 4,71е-1

Н2 0,7 9,91е-2 1,15е-1 2,16-1 2,16е-1 2,48е-1 2,47е-1

(1,5) 1,4 3,79е-1 3,80е-1 4,27е-1 4,25е-1 4,81е-1 4,78е-1

Таблица 2

В настоящее время ведется работа по разработке метода расчета временных характеристик исследуемой СМО. Далее станет возможен расчет сети массового обслуживания с узлами типа М/О/п-О и ограничением полного времени пребывания заявки в ней.

Автор выражает искреннюю благодарность научному руководителю профессору Юрию Ивановичу Рыжикову за ряд ценных рекомендаций при написании данной статьи.

Литература

1. Башарин, Г.П. Анализ очередей в вычислительных системах / Г.П. Башарин, П.П. Бочаров, Я.А. Коган. - М., 1989.

2. Бочаров, П.П. Теория массового обслуживания: Учебник / П.П. Бочаров, А.В. Печинкин. - М., 1995.

3. Дудин, С.А. Модель функционирования колл-центра как система МАР/РН/М/Я-М с нетерпеливыми запросами / С.А. Дудин, О.С. Дудина // Проблемы передачи информации. - 2011. - № 47. - С. 68 - 83. - аог 10.1134/ 80032946011040053.

4. Рыжиков, Ю.И. Имитационное моделирование систем с «нетерпеливыми» заявками / Ю.И. Рыжиков, А.В. Уланов // Имитационное моделирование. Теория и практи-

ка: тр. VI Всероссийской конференции. - Казань, 2013. -С. 339 - 342.

5. Рыжиков, Ю.И. Итеративный метод расчета многоканальных систем с произвольным распределением времени обслуживания / Ю.И. Рыжиков, А. Д. Хомоненко // Проблемы управления и теории информации. - 1980. -№ 3. - С. 203 - 213.

6. Рыжиков, Ю.И. Опыт расчета сложных систем массового обслуживания / Ю.И. Рыжиков, А.В. Уланов // Информационно-управляющие системы. - 2009. - № 2. -С. 56 - 62.

7. Hoshi, K. TrafficPerfomance for a Time-Out Scheme Communication System / K. Hoshi, S. Iijima, Y. Takahashi, N. Komatsu // Proc. International Conference on Ultra Modern Communications "ICUMT'2009". - St. Petersburg. - 2009. -P. 1 - 6. - doi: 10.1109/ICUMT.2009.5345457.

8. Roubos, A. Call Centers with Hyperexponential Patience Modeling / A. Roubos, O. Jouini // International Journal of Production Economics. - 2013. - V. 141. - P. 307 - 315. -doi: 10.1016/j.ijpe.2012.08.011.

9. Takahashi, Y. Numerical Method for the Steady-State Probabilities of a GI/G/c Queuing System in General Class / Y. Takahashi, Y.A. Takami // Journal of the Operations Research Society of Japan. - 1976. - V. 19. - № 2. - P. 147 -155.

УДК 621.778.04

Р.А. Юдин, Н.И. Шестаков, И.Р. Юдин, Н.А. Тувалин

ОСОБЕННОСТИ РАСЧЕТА И РЕГУЛИРОВАНИЯ ПРОЦЕССОВ ДВУХСТАДИЙНОГО СЖИГАНИЯ ПРИРОДНОГО ГАЗА

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

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

The article presents the results of the algorithm for calculating the incomplete combustion of natural gas, which allows calculating the chemical composition of the products of incomplete combustion and temperatures developed by combustion, at various regime parameters. The algorithm used the dependence of thermodynamic equilibrium constant of the reaction of water vapor from the stoichiometric amounts of products of complete combustion. The article also provides opportunities to improve the automatic temperature control in furnace at two-stage burning of natural gas.

Oxide and carbon dioxide, hydrogen and water vapor flow rate of primary air, constant of thermodynamic equilibrium, inertia, automatic regulation.

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

ванно обеспечивает малоокислительный нагрев металлопродукции.

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

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