В настоящее время ведется работа по разработке метода расчета временных характеристик исследуемой СМО. Далее станет возможен расчет сети массового обслуживания с узлами типа MIGIn-G и ограничением полного времени пребывания заявки в ней.
Автор выражает искреннюю благодарность научному руководителю профессору Юрию Ивановичу Рыжикову за ряд ценных рекомендаций при написании данной статьи.
Литература
1. Башарин, Г.П. Анализ очередей в вычислительных системах I Г.П. Башарин, П.П. Бочаров, Я.А. Коган. - М., 1989.
2. Бочаров, П.П. Теория массового обслуживания: Учебник I П.П. Бочаров, А.В. Печинкин. - М., 1995.
3. Дудин, С.А. Модель функционирования колл-центра как система MAPIPHINIR-N с нетерпеливыми запросами I С.А. Дудин, О.С. Дудина II Проблемы передачи информации. - 2011. - № 47. - С. 68 - 83. - аог 10.1134! 80032946011040053.
4. Рыжиков, Ю.И. Имитационное моделирование систем с «нетерпеливыми» заявками I Ю.И. Рыжиков, А.В. Уланов II Имитационное моделирование. Теория и практи-
ка: тр. 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.
В черной металлургии и машиностроении при термохимической обработке стали широко используют открытый малоокислительный нагрев, применяя двухстадийное сжигание природного газа в печах с последовательными и параллельными зонами сжигания и дожигания. Окончательную термообработку проводят в атмосфере продуктов неполного сгорания природного газа, а предварительную - в атмосфере продуктов их дожигания, что гарантиро-
ванно обеспечивает малоокислительный нагрев металлопродукции.
Расчет дожигания продуктов неполного сгорания, при котором обеспечивается полнота сгорания топлива при минимальных избытках воздуха, не представляет трудностей. Для получения итоговых химических реакций полного сгорания любого вида топлива достаточно ограничиться прямым уравниванием начальных и конечных элементов горения в
левой и правой частях уравнений этих реакций. При неполном сгорании топлива такой подход неприемлем, так как химический состав продуктов неполного сгорания во многом зависит от термодинамических особенностей процессов горения, поэтому его необходимо рассчитывать с их учетом.
Так, до последнего времени при расчетах химического состава продуктов неполного сгорания топлива использовали методику В.Ф. Копытова [2], основанную на принципе термодинамического равновесия компонентов продуктов неполного сгорания по обратимой физико-химической реакции, иногда именуемой реакцией водяного пара:
СО2 + Н2 О СО + Н2О.
(1)
Константа термодинамического равновесия этой реакции является табулированной величиной, значения которой в зависимости от температуры определяют по формуле:
и тт
к = —= / Ф,
Усо ' 2Н,О
(2)
где означенные неизвестные являются объемами указанных компонентов продуктов неполного сгорания топлива.
По методике В. Ф. Копытова для определения химического состава продуктов неполного сгорания составляют материальный баланс химических элементов, участвующих в реакции (1). Первоначально составляют систему, состоящую из трех линейных алгебраических уравнений, соответственно для углерода, водорода и кислорода. В левой части каждого из этих уравнений записывают количество данного элемента в продуктах неполного сгорания, а в правой - в исходном топливе и воздухе, подаваемом на горение, при различных коэффициентах расхода первичного воздуха а!.
Четвертым уравнением этой системы является формула (2), при подстановке в которую трех неизвестных образуется квадратное уравнение, которое следует решать относительно одного из компонентов реакции (1).
В монографиях [1], [3] показано, что рассматриваемую систему уравнений можно решить относительно одного из неизвестных компонентов, используя в правой части уравнений вместо элементов исходного топлива, стехиометрические параметры его продуктов сгорания. Полученное в [3] решение для газа произвольного химического состава имеет следующий вид:
к =-
(СО2Шах + СО2 - Усо )
Здесь СО2шах и Н2Ошах - максимальные объемы диоксида углерода и водяных паров в продуктах сгорания, соответствующие стехиометрическим условиям горения; СО, СО2 - объемы оксида и диоксида углерода в топливе; 8О2 - объем сернистого газа в продуктах сгорания, равный объему сероводорода (Н28) в топливе; а! -коэффициент расхода первичного воздуха. Для топлива, содержащего только углеводородные горючие газы и негорючие компоненты, формула (3) имеет более компактный вид:
к =
(СО2шах + СО2 - Усо) •[(!-а.) • (2СО2шах + Н^) - уго] УсО[а1(2СО2шах + Н.Ошах) - 2СО2шах + уго]
(3.1)
Для метана СО2шах = 1, Н2Ошах = 2, поэтому:
(1 - Усо) • [4(1 -а1) -Усо]
к =-
Усо(4а1 - 2 + Усо)
(3.2)
к,усо2 - [3 - 2к1 (2а, -1)]усо + 4(1 - а,). (4)
Квадратное уравнение (3.2) имеет следующее решение:
Усо =
3 - 2к1 (2а1 -1) - ^[3 - 2к1 (2а1 -1)]2 - 16к1 (1 - а1)
2к1
(5)
где к1 = 1 - к. Ее зависимость от температуры приведена в таблице.
Таблица
Зависимость от температуры к1 = (1 - к) = Ж (1)
Температура, °С 1200 1300 1400 1500 1600 1700 1800
Конста нта равнове сия, к 0,607 0,655 0,703 0,74 0,758 0,776 0,789
Температура, °С 1900 2000 2100 2200 2300 2400
Конста нта равнове сия, к! 0,801 0,813 0,825 0,837 0,841 0,845
Усо [а1(2со2шах + н2ош
х) - 2СО2шЯХ + а1(28О2 - СО) - 2О2 + Усо ]
[(1 -а1)(2СО2шах + Н2Ош
х)-а1(28О2 - СО) + 2О2 - Усо ].
(3)
Выявленные закономерности позволили разработать алгоритм расчета химического состава продуктов неполного сгорания топлива и температур горения при различных коэффициентах расхода первичного воздуха а1 и
температурах его нагрева. Структура алгоритма представлена на рисунке. Блок ввода исходных данных 1 алгоритма содержит три секции. В первой секции заданы табулированные зависимости константы термодинамического равновесия от температуры [3, с. 387], во второй - химический состав компонентов газообразного топлива и низшая теплота сгорания, а в третьей - табулированные зависимости теплоемкости компонентов продуктов неполного сгорания и воздуха от их температуры [3, с. 385].
В блоке 2 алгоритма по формулам, полученным в монографии [3, с. 78] и работе [5], рассчитывают стехиометрические объемы продуктов полного сгорания топлива:
СО2ШЯХ = СО + П(СпН2П+2 + СПН2П + СПНП ),
И2Ошах = Н2 + + (П + 1)СПН2П+2 + ПВДП + 0,5пСпип ,
(6)
где СпН2п+2, СпН2п, СпНп - содержание алканов, алкенов и алкинов в исходном топливе в м^м3. Необходимый объем воздуха определяют по формуле:
V0 = 4,76 (СО2шах + 0,5 Н2Ошах + 1,5 ЪБ - 0,5 СО).
(7)
В блоке 3 проводят расчет фактических объемов азота и воздуха, предварительно задавая коэффициент расхода первичного воздуха а1. Кроме того, в блоке (3) по формуле (3) и ее частным случаям составляют квадратное уравнение, подобное (4) для метана, корнем которого является уСО.
Для определения его численного значения предварительно задают, искомое значение температуры горения топлива, по которому в секции 1 блока 1 определяют табулированное значение константы равновесия к1. Результаты расчетов блока 3 направляют в блок 4, где рассчитывают объемы трех остальных компонентов продуктов неполного сгорания (СО2, Н2, Н2О) и направляют результаты расчетов в блок 5. По рассчитанным объемам оксида углерода и водорода в блоке 5 определяют количество тепла, выделяемого при экзотермических реакциях неполного сгорания топлива:
Qэ = 0нР - 12645 уСО - 10786 ин, кД^м3, (8)
где 12645 кДжЛм3 - теплота сгорания оксида углерода, 10786 кДжЛм3 - теплота сгорания водорода.
Полученное количество тепла экзотермической реакции неполного сгорания топлива - Qэ в блоке 5 передают в блоки 6 и 7, где соответственно рассчитывают температуры для холодных газо-воздушных смесей и калориметрические температуры, развиваемые при неполном сжигании топлива при нагреве воздуха. В эти блоки из блоков 2 - 4 также передают полные данные об объемах компонентов продуктов сгорания. Для расчета калориметрической температуры в блок 7 из блока 3 дополнительно передают
результаты расчета объема воздуха, соответствующее заданному коэффициенту его расхода аь
В блоки 6, 7 из третьей секции блока исходных данных 1 также передают табулированные значения теплоемкости компонентов продуктов неполного сгорания, соответствующие заданным температурам.
Кроме того, в блок 7 из третьей секции блока 1 направляют табулированное значение теплоемкости воздуха при заданной температуре его нагрева. Калориметрическую температуру рассчитывают по формуле:
к =_QL+Q±_
кал С • V + С • V + С • V + С • V + С •V '
(9)
где Qв - тепло, вносимое с нагретым воздухом. В зависимости от полученных температур задаются новыми значениями /х и /кал и, начиная с блока 3, используя метод итераций, повторяют расчет, обеспечивая заданную сходимость результатов. Для уменьшения числа итераций целесообразно использовать данные, полученные в результате расчетов по предлагаемому алгоритму для метана в монографии [3, с. 344 - 369] и работах [4] - [6].
Наибольший практический интерес представляют результаты расчета сжигания природного газа, приведенные в монографии [3, с. 362 - 363], при коэффициенте расхода первичного воздуха а1 = 0,60 и работе с «холодным» воздухом. Это обусловлено тем, что указанные условия весьма близки к верхней концентрационной границе воспламенения метана -15 %.
Так, физико-химическая реакция при а! = 0,60 соответствует 14,9 % содержания метана в горючей смеси:
СН4 + 1,2 О2 + 4,512 N = 0,2585 СО2 + 1,1415 Н2О +
+ 0,7415 СО + 0,8585 Н2 + 4,512 М2 + 17183,9 кДж.
(10)
Расчеты показали, что этой реакции соответствует температура 1494 °С. Кроме того, расчеты в монографии [3, с. 368 -369], показали, что при нагреве воздуха до 400 °С коэффициент расхода первичного воздуха можно снизить до а1 = 0,50, так как этим условиям отвечает калориметрическая температура 1494 °С. Это убедительно свидетельствует о возможности реализации физико-химической реакции:
СН4 + О2 + 3,76 1Ч2 = 0,1598 СО2 + 0,8402 Н2О +
+ 0,8402 СО + 1,1598 Н2 + 3,76 М2 + 15218,8 кДж.
(11)
Из (11) следует, что подогрев воздуха до 400 °С позволяет повысить верхний предел воспламенения метана с 15 до 17,4 %. При этом согласно работам [1] - [2] подавляется образование сажистого углеро-
да. Кроме того, согласно реакциям (10) - (11) указанное снижение коэффициента а1 позволяет увеличить содержание оксида углерода в сухих продуктах сгорания с 11,5 до 16,2 %, а водорода - с 13,0 до 19,6 % и тем самым существенно повысить восстановительные свойства продуктов неполного сгорания и, как следствие, качество поверхности металлопродукции.
Следует отметить, что коэффициенты перед всеми компонентами продуктов неполного сгорания физико-химических уравнений (10) - (11) приведены с точностью до четвертого знака после запятой. Это убедительно свидетельствует о том, что они получены достаточно сложным путем с использованием численных методов расчета.
1. Исходные данные
Константа термодинамического равновесия: к =А(г) Химический состав и теплота сгорания топлива (9нР), кДж / м3 Теплоемкость продуктов сгорания и воздуха: С1, в = Аг), кДж/м3 К
2. Стехиометрические объемы компонентов продуктов сгорания:
СО2max, Н2Ошax, N2 ,
воздуха горения: ¥в0, м3
3. РасчетуСО = /(к, а!), м , (1 - к)у2 - Ву + С = 0, ^ В,С = / (СОгшах, Н2Ош3х,а1)
4. Расчет х, и, г =А (у), м3:
ХСО2 = СО2т + СО2шах - УC0,
N2 = а1N20 + N2, , ¥в = а! Ув°, м3
5. Расчет тепла экзогазовых реакций горения, кДж/м3 :
0Э = бнР -12645усо - 10800ин
6. Жаропроизводительность и температуры горения, °С:
Гмакс, X = Qэ /
7. Калориметрическая температура, °С: Гкал = (9з + Св^вГв) /
8. Вывод результатов
СО2шах, Н2Ошах, ¥N2°, Ув0, Усо2, Усо, Ун2о, Ун2, VN2, У802, Qэ, 1макс, 1кал
Рисунок. Алгоритм расчета неполного сгорания топлива
В настоящее время печи открытого малоокислительного нагрева с последовательными зонами сжигания и дожигания в зависимости от сортамента металлопродукции работают с коэффициентами расхода первичного воздуха а1 = 0,50 - 0,60. При а1 = 0,50 дожигание продуктов неполного сгорания проводят с коэффициентами расхода вторичного воздуха а2 = 0,55 - 0,60, а при а1 = 0,60 - с а2 = 0,45 - 0,50. В обоих вариантах суммарный коэффициент расхода воздуха составляет аЕ = 1,05 - 1,10.
В ходе эксплуатации протяжных печей метизного производства с последовательными зонами сжигания и дожигания при а1 = 0,50 - 0,55 была выявлена высокая инерционность изменения температуры зон окончательного нагрева проволоки, несмотря на существенное изменение расходов газа и воздуха при заданном соотношении. При этом зонах окончательного нагрева, в которых осуществляется дожигание продуктов неполного сгорания и через сопла подают только вторичный воздух, при изменении всех режимных параметров печи температура изменялась практически мгновенно.
В связи с этой особенностью был сделан вывод о целесообразности организации автоматической связи блока регулирования соотношения природного газа и первичного воздуха с контурами регулирования температуры в обеих зонах печи, а не только с зоной окончательного нагрева, как было общепринято в существующих системах автоматического регулиро-
вания печей открытого малоокислительного нагрева металлопродукции с двухстадийным сжиганием природного газа. Эксплуатация протяжных печей метизного производства показала, что указанная дополнительная связь позволяет исключить негативное влияние инерционности изменения температуры зоны окончательного нагрева на качество работы автоматики.
Литература
1. Ващенко, А.И. Окисление и обезуглероживание стали / [А.И. Ващенко и др.]. - М., 1972.
2. Копытов, В.Ф. Сжигание газа с недостатком воздуха / В.Ф. Копытов // Сталь. - 1954. - № 3. - С. 57 - 62.
3. Юдин, Р.А. Научно-технические основы процессов горения / Р. А. Юдин. - Череповец, 2013.
4. Юдин, Р.А. Основные закономерности неполного сгорания природного газа / [Р.А. Юдин и др.] // Вестник Череповецкого государственного университета. - 2013. -№ 3. - Т. 1. - С. 26 - 29.
5. Юдин, Р.А. Особенности горения топлива произвольного химического состава / Р.А. Юдин, И.Р. Юдин // Вестник Череповецкого государственного университета. -2012. - № 1. - С. 15 - 18.
6. Юдин, Р.А. Разработка алгоритма расчета неполного сгорания природного газа / [Р.А. Юдин и др.] // Вестник Череповецкого государственного университета. - 2013. -№ 3. - Т. 2. - С. 31 - 35.