CC BY
126
Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. Т. 62, № 4 (2019), с. 341-361 Energetika. Proc. CIS Higher Educ. Inst. and Power Eng. Assoc. V. 62, No 4 (2019), pp. 341-361 341
https://doi.org/10.21122/1029-7448-2019-62-4-341-361 УДК 662.765.1
Новый метод расчета состава древесного генераторного газа обращенного процесса газификации
Е. М. Кашин1*, В. Н. Диденко1*
'-Ижевский государственный технический университет имени М. Т. Калашникова (Ижевск, Российская Федерация)
© Белорусский национальный технический университет, 2019 Belarasian National Technical University, 2019
Реферат. В статье излагается новый метод определения состава древесного генераторного газа, получаемого в газогенераторах обращенного процесса газификации. Анализируются недостатки существующих расчетных методов, главным из которых является плохое согласование результатов расчетов с опытными данными. Обосновывается приоритет основных химических реакций, протекающих при газификации древесного топлива. Выделяются три активные зоны газификации: окислительно-восстановительная, восстановительная и зона взаимодействия продуктов газификации между собой и с углеродом топлива. Окислительно-восстановительная зона в общем случае состоит из двух подзон: в первой протекают реакции образования водяного газа, а вторая появляется при подаче в газогенератор избыточного воздуха. Предлагаемый метод расчета компонентов генераторного газа - это совокупность модифицированного балансового и добавленного методов расчета концентраций продуктов химических реакций по константам равновесия этих реакций в активных зонах газификации с различными температурами. В модифицированном балансовом методе рассматриваются первичные процессы превращения древесины и влажного воздуха в компоненты генераторного газа в первой подзоне окислительно-восстановительной зоны. Модифицированный балансовый метод основан на уравнениях материального баланса углерода, водорода, кислорода, влаги, азота и теплового баланса системы. В добавленном методе определяются концентрации компонентов генераторного газа во второй подзоне окислительно-восстановительной зоны, а также в восстановительной зоне и зоне взаимодействия продуктов газификации между собой и с углеродом топлива. Совокупность этих двух методов позволяет расчетным путем с большой точностью определять выход генераторного газа, концентрации его компонентов, расход топлива и воздуха, а также ряд других характеристик газогенератора.
Ключевые слова: концентрация компонента газа, уравнение баланса углерода, уравнение баланса водорода, уравнение баланса кислорода, уравнение баланса влаги, константа равновесия химической реакции, гомогенная реакция, гетерогенная реакция
Для цитирования: Кашин, Е. М. Новый метод расчета состава древесного генераторного газа обращенного процесса газификации / Е. М. Кашин, В. Н. Диденко // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. 2019. Т. 62, № 4. С. 341-361. https://doi. org/10.21122/1029-7448-2019-62-4-341-361
Адрес для переписки Address for correspondence
Кашин Евгений Михайлович Kashin Evgeny M.
Ижевский государственный технический Kalashnikov Izhevsk State Technical
университет имени М. Т. Калашникова University
ул. 7-я Подлесная, 32а, 32a 7 Podlesnaya str.
426069, г. Ижевск, Российская Федерация 426069, Izhevsk, Russian Federation
Тел.: +7 3412 77-60-55 доб. 8102 Tel.: +7 3412 77-60-55 / 8102
[email protected] [email protected]
The New Method for Determining the Composition of Wood Gas Produced in Gas Generators of the Inverted Gasification Process
E. M. Kashin1*, V. N. Didenko1*
'-Kalashnikov Izhevsk State Technical University (Izhevsk, Russian Federation)
Abstract. The article presents a new method for determining the composition of wood generator gas produced in gas generators of the inverted gasification process. The shortcomings of the existing calculation methods are analyzed, the main of which is the insufficient harmonization of the calculation results with the experimental data. The authors substantiate the priority of the main chemical reactions occurring during gasification of wood fuel. There are three active zones of gasification, viz.: a redox zone, a reduction zone and a zone of interaction of gasification products with each other and with the carbon of the fuel. In general, a redox zone consists of two subzones: in the first one reactions of water gas formation occur, and the second one appears when excess air is supplied to the gas generator. The proposed method for calculating the components of the generator gas is a set of a modified balance method and an added method for calculating the concentrations of chemical reaction products by the equilibrium constants of these reactions in the active gasification zones with different temperatures. The modified balance method considers the primary processes of wood and moist air transformation into components of the generator gas in the first subzone of the redox zone. The modified balance method is based on the equations of material balance of carbon, hydrogen, oxygen, moisture, nitrogen and thermal balance of the system. The added method determines the concentrations of the components of the generator gas in the second subzone of the redox zone, as well as in the reduction zone and the zone of interaction of the gasification products with each other and with the fuel carbon. The combination of these two methods makes it possible to calculate with greater accuracy the output of the generator gas, the concentration of its components, fuel and air consumption, as well as a number of other characteristics of the gas generator.
Keywords: concentration of gas component, equation of carbon balance, equation of hydrogen balance, equation of oxygen balance, equation of water balance, balance constant of chemical reaction, homogeneous reaction, heterogeneous reaction
For citation: Kashin E. M., Didenko V. N. (2019) The New Method for Determining the Composition of Wood Gas Produced in Gas Generators of the Inverted Gasification Process. Energetika. Ргос. CIS Higher Educ. Inst. and Power Eng. Assoc. 62 (4), 341-361. https://doi.org/10.21122/1029-7448-2019-62-4-341-361 (in Russian)
Введение
Газогенератор обращенного процесса газификации по отношению к газогенератору прямой схемы обладает рядом преимуществ, главными из которых являются: работа на различных видах топлива, низкое содержание смол в продуктах газификации, загрузка топлива без остановки процесса газификации, сонаправленность центробежных сил и потока генераторного газа (для газогенераторов роторного типа [1]). Считается, что в газогенераторе обращенной схемы абсолютно все элементы топлива, включая влагу, участвуют в процессе образования газа. Для расчета состава газа обращенного процесса газификации используется общепринятый балансовый метод, основанный на составлении материального и теплового балансов по заданному составу топлива [2-5]. Метод имеет ряд недостатков:
• результаты расчетов плохо согласуются с опытными данными (табл. 1);
• полагается, что влага диссоциирует на водород и кислород, хотя температура в древесных газогенераторах намного ниже 2000 °C;
• не учитывается химическое взаимодействие продуктов газификации друг с другом и с углеродом вновь поступающего нагретого топлива в других температурных областях реальной активной зоны газификации;
• в явном виде не учитывается коэффициент избытка воздуха;
• содержание метана в продуктах газификации не определяется и в окончательных расчетах задается величиной СН4 = 1,5 %;
• не учитывается влага воздуха, используемого в качестве газифицирующего агента;
• для широкого диапазона изменения температуры в зоне газификации константа равновесия химической реакции водяного газа принимается средним значением в интервале от 1,7 до 2,5.
Перечисленные недостатки балансового метода [2] делают актуальной разработку нового метода расчета состава древесного генераторного газа (табл. 1).
Таблица 1
Сухой состав генераторного газа (%) по расчетам методики Н. Г. Юдушкина [2], по опытным данным Д. Б. Гинзбурга [3] и В. А. Лямина [6]
Dry composition of the generator gas (%) calculated by Yudushkin's method [2], according to Ginzburg's experimental data [3] and Lyamin's data [6]
CO CO2 H2 O2 CH4 n2
Н. Г. Юдушкин [2]. Топливо - береза, влажность W = 28 %, T = 1500 K, K = 2,49
19,1 11,4 16,6 - 1,5 51,5
Д. Б. Гинзбург [3]. Топливо - береза, W = 28 %
27,5 7,6 11,4 0,4 1,9 50,8
В. А. Лямин [6]. Топливо - береза, W = 15-20 %
30,8 4,3 5,6 0,3 1,4 57,6
Основная часть
Разработан метод расчета состава генераторного газа, представляющий совокупность балансового метода [2], модифицированного авторами статьи, и добавленного метода расчета концентраций продуктов химических реакций по константам равновесия этих реакций в активных зонах газификации с различными температурами.
Для разработки нового метода расчета состава древесного генераторного газа потребовалось установить приоритеты и границы температурных зон протекания химических реакций в газогенераторе.
Способность веществ химически взаимодействовать друг с другом при T, P = const определяется величиной и знаком изменения свободной энтальпии (энергии Гиббса)
AG = AH - T AS,
где АН - изменение энтальпии (теплового эффекта) химической реакции; AS - изменение энтропии химической реакции.
Химическая реакция при условиях T, Р = const возможна, если AG < 0. При AG = 0 реакция находится в состоянии равновесия. При AG > 0 химическая реакция невозможна. Чем больше по модулю отрицательная величина AG, тем больше реакционная способность веществ, и доля реакции в системе химических превращений.
Доля 7-й реакции в системе химических превращений определяется отношением
В случаях, когда АО > 0 и химическая реакция невозможна, значение АО7 и доля реакции принимаются равными 0.
Величины АНр7, ASрi, АО, и для 7-х реакций, полученные по значениям изменений энтальпии и энтропии реагентов и продуктов реакций [7], приведены в табл. 2.
В предлагаемом методе, в отличие от существующего, выделяются не одна, а три активные зоны газификации (рис. 1, табл. 2):
1 - окислительно-восстановительная зона (ОВЗ), состоящая в общем случае из двух подзон (1.0 и 1.1). В подзоне 1.0 (БОВП) протекают реакции образования водяного газа. Температура в окислительно-восстановительной зоне в зависимости от теплоизоляции газогенератора составляет 60-85 % от адиабатной температуры горения топлива. При подаче в газогенератор избыточного воздуха (а > 1) появляется подзона 1.1 (ПДО), где идут дополнительные окислительные реакции;
2 - восстановительная зона (ВЗ), где отсутствует окислитель и протекают только восстановительные реакции. По величине АО этих химических реакций получено, что температура в восстановительной зоне выше 980 °С;
3 - зона взаимодействия продуктов газификации между собой и с углеродом топлива (ЗВП). По величине АО химических реакций, приведенных в табл. 2, получено, что температура в этой зоне ниже 980 °С.
Рис. 1. Структура активной зоны газификации в газогенераторе роторного типа при центральной подаче газифицирующего агента (обращенная схема газификации)
Fig. 1. The structure of the gasification zone in rotary gas generator at the central supply of the gasifying agent (reversed gasification scheme)
AG
Газ Зола
Таблица 2
Приоритет химических реакций в различных температурных зонах газогенератора
The priority of chemical reactions in different temperature zones of the gas generator
№ реакции в активной зоне Реакция кДж/моль ASpi, Дж/(моль-К) T, K AG,, Дж/моль ZAG, Дж/моль доли ед.
1. Окислительно-восстановительная зона
1.0 Балансовая окислительно-восстановительная подзона
0 CO2 + H2 о CO + H2O 41,17 41,94 1500 -21740 -21740 1
1.1. Подзона дополнительного окисления избыточным кислородом при а > 1
1 C + O2 о CO2 -393,51 2,83 1500 -397755 -1446410 0,27
2 2C + O2 о 2CO -221 178,29 1500 -488435 0,34
3 2CO + O2 о 2CO2 -172,51 24,77 1500 -209665 0,14
4 2H2 + O2 о 2H2O -483,68 -88,75 1500 -350555 0,24
2. Восстановительная зона
1 CO2 + H2 о CO + H2O 41,17 41,94 1200 -9158 -95810 0,10
2 С + CO2 о 2CO 172,51 175,46 1200 -38042 0,40
3 C + H2O о CO + H2 131,34 133,52 1200 -28884 0,30
4 C + 2H2O о CO2 + 2H2 90,17 91,58 1200 -19726 0,21
3. Зона взаимодействия продуктов между собой и с углеродом топлива
1 CO + H2O о CO2 + H2 -41,17 -41,94 950 -1327 -5267 0,25
2 C+2H2 о CH4 -74,85 -80,75 950 0 0,00
3 CO + 3H2 о CH4 + H2O -206,19 -214,27 950 -2633,5 0,50
4 CO2 + 4H2 о CH4 + 2H2O -165,02 -172,33 950 -1306,5 0,25
Предложенный метод позволяет вычислить выход, концентрации компонентов генераторного газа, определить расход топлива и воздуха, а также ряд других характеристик работы газогенератора.
Генераторный газ в газогенераторе считается идеальным, поэтому объемное содержание компонентов в газе соответствует их мольному содержанию. Состав генераторного газа зависит от активной зоны и имеет следующий наиболее общий вид:
СО + С02 + СН4 + Н2 + Н 2 О+N 2 + 0^ = 100%,
где СО, С02, СН4, Н2, Н20, N ,02з6 - процентное содержание соответствующих газов по объему.
В отличие от [2], состав газа включает 02з6 и СН4. Наличие в составе газа неизрасходованного кислорода 02з6 возможно только на выходе из подзоны 1.0 (БОВП) окислительно-восстановительной зоны 1 при а > 1. Метан СН4 появляется только на выходе из зоны 3.
Модифицированный балансовый метод
(материальный и тепловой балансы)
Исходными составляющими материального баланса являются влажные топливо и воздух, а конечными - влажный генераторный газ, недожог и зола. В отличие от [2], в качестве газифицирующего агента рассматривается не сухой, а влажный воздух. Тепловой баланс составляется по низшей теплоте сгорания топлива, определяемой расчетным путем.
Элементарный состав твердого органического топлива (по рабочей массе)
Ср + Нр + Ор + N + Sp + Wр + Ар = 100 %,
где Ср, Нр, Ор, Nр, 8р ,Wр, Ар - массовые процентные содержания элементов в исходном рабочем твердом топливе. Для древесного топлива 8р = 0.
Для удобства работы с малоразмерными коэффициентами уравнений расчет производится для 100 кмоль влажного генераторного газа.
Для определения выхода ивг из 1 кг топлива, содержания компонентов генераторного газа СО, СО2, Н2, Н2О, N и избыточного кислорода О2 в конце первой активной зоны требуется семь уравнений, составленных из материального и теплового балансов процесса газификации, учитывающих константу равновесия реакции образования водяного газа.
Переход элементов топлива (Ср, Нр, Ор, №, WрP) и воздуха (О2 ,N2, влага) в генераторный газ осуществляется по четырем основным уравнениям материального баланса углерода, водорода, кислорода и влаги.
Объем 100 кмоль генераторного газа при нормальных физических условиях (^ = 0 0С, Р = 760 мм рт. ст.) по закону Авогадро
Увх = 100• 22,4 м3.
Масса рабочего топлива, кг, необходимая для получения Увг газа:
„ УВг 100 • 22,4
От =^ =-,
и и
в.г в.г
где ивг - выход влажного газа из 1 кг топлива, м3/кг.
Уравнение баланса углерода
Баланс углерода заключается в том, что количество углерода до и после химических реакций не изменяется. Количество углерода (кмоль), содержащегося в топливе массой От, определяется по уравнению
= О Ср 100• 22,4Ср 1,867Ср
Ст = "
100ЦС ив.г 100 42 ив.г
где цС - молекулярная масса углерода.
В первой активной зоне углерод в количестве Сд переходит в газообразные компоненты СО и СО2
Сд = Ср ( 100 - 5 )/100,
где Срд - действительное процентное содержание углерода по рабочей массе, необходимое для образования 100 кмоль влажного генераторного газа; 5 - потери углерода, связанные с образованием золы и шлака, %. Тогда уравнение баланса углерода запишется в виде
1,867Сд
--д = СО + СО2. (1)
Уравнение баланса водорода
Весь водород Нр твердого топлива после химического взаимодействия превращается в водород генераторного газа (содержанием водорода в шлаке пренебрегается ввиду незначительности).
Количество водорода (кмоль) в массе От рабочего топлива определяется по уравнению
нр = От Нр = 100 • 22,4Нр = 11,2Нр
100ц„ ивг100•2 и
в.г
где цН - молекулярная масса водорода.
Тогда уравнение баланса водорода запишется в следующем виде
112Нр=Н,_. (2)
ив.г
Уравнение (2) отличается от [2] тем, что не включает водород от диссоциации паров воды на Н2 и О2.
Уравнение баланса кислорода
Количество кислорода (кмоль) в рабочей массе топлива
Ор = ОтОр = 100 • 22,4Ор = 0,7Ор
т = 100 ц02 = ив.г 100^32 =~UB7,
где Цо - молекулярная масса кислорода.
Количество воздуха (кмоль), необходимое для сжигания рабочей массы топлива От при коэффициенте избытка воздуха а:
ив.г
У= Mв_ От V Рв а 100 • 22,4 У0 рв а
Цв Цв ив.г Цв
где Мв - масса воздуха, необходимая для сжигания От (кг) топлива при коэффициенте избытка воздуха а, кг ; У0 - теоретический объем воздуха, необходимый для полного сжигания 1 кг рабочего топлива при а = 1 м3/кг; рв - плотность влажного воздуха при нормальных физических условиях, кг/м3; цв - кажущаяся молекулярная масса влажного воздуха, кг/моль;
рв =—---0,0013 фР^; Цв = 28,96-10,944 фр,
^ Т Т В
В - барометрическое давление, Па; ф - относительная влажность, доли ед.; Рн - давление насыщения водяного пара при температуре Т, Па; Лв -удельная газовая постоянная сухого воздуха, Дж/(кг-К).
Давление насыщения водяного пара Рн при температуре Т определяется по таблицам насыщенного водяного пара или по приближенной зависимости
Рн = 479 + (11,52 +1,62 (Т - 273))2.
Полагается, что объемная доля кислорода в воздухе равна 0,21, и уравнение баланса кислорода записывается в следующем виде:
0 70р
--+ 0,21у в = С02 + 0,5 (СО + Н20) + 0,2^, (3)
ив.г
где 0,21ув - количество кислорода в подаваемом воздухе; 0,21у^эб - количество кислорода избыточного воздуха в составе продуктов газификации в окислительно-восстановительной зоне 1;
уИзб = 100 • 22,4 У о Рв в р = | 0 пРи а^1;
в ив.г Цв , 1(а-1) пРи а>1.
В (3), в отличие от [2], в явном виде учтен коэффициент избытка воздуха а и отсутствует О2 от диссоциации паров воды.
Уравнение баланса влаги
Количество влаги (кмоль), содержащейся в топливе массой От: н 0Р = 100 • 22,4 Wр = 1,244Wр
2 т = ивТ 100 • 18 = ив, . Количество влаги (кмоль), содержащейся в подаваемом воздухе:
в MН0
ЦН20
где 0 - масса водяного пара, содержащегося во влажном воздухе, кг; цН 0 - молекулярная масса воды, кг/моль.
Масса водяного пара в воздухе определяется по известному влаго-содержанию воздуха d
d 100 • 22,4 V рв а
MH O =---———.
H2° 1 + d u„
Тогда количество влаги (кмоль), содержащейся в воздухе:
в = 100 • 22,4 V рв а
vh2o =, , .
2 1 +d ив.г M'H2O
Окончательно уравнение баланса влаги запишется следующим образом:
1 244Wp
1,244W + VH20 =H2O. (4)
ив.г 2
В (4), в отличие от [2], учтено влагосодержание подаваемого воздуха.
Так как диссоциация водяного пара на водород и кислород возможна только при температурах выше 2000 °С, недостижимых для древесных газогенераторов, то вся влага топлива и влажного воздуха полностью переходит в водяной пар генераторного газа.
Уравнение баланса азота
Замыкающим уравнением материального баланса является уравнение баланса азота
100 - (CO + C02 + H2 + H20 + °2з6)=N2.
Количество азота до и после процессов горения и газификации остается практически неизменным, так как азот не участвует в химических реакциях, поэтому количество азота в генераторном газе N2 принимается равным количеству азота в воздухе N2 [8, 9]
N2 ^N2 = 0,79vв = 0,79•100•22,4^.
ив.г М'в
Принимается, что объемная доля азота в воздухе составляет 0,79.
И окончательно уравнение баланса азота запишется в виде
100 - (CO + C02 + H2 + H2O + 0,21 vиз5) = 0,79 Vв. (5)
Уравнения теплового баланса
Содержание СО, Н2О, СО2 и Н2 в окислительно-восстановительной зоне 1 зависит от температуры и определяется константой равновесия химической реакции образования водяного газа:
С0+Н20 ^С02 + Н2;
к =С01Н20. (6)
' С02 • Н 1 '
По результатам исследований [4, 10-13] константу равновесия химической реакции можно найти по уравнению
2207 2
^ Ку =--^ + 0,9115^ T - 9,738 •Ю-4 T + 1,487 •Ю-7 T2 + 0,098.
Вторым уравнением теплового баланса является уравнение теплоты процесса газификации. Из определения условного КПД газогенератора следует
ив.г ОТ
Пг =-
О
где ивг - выход влажного газа из 1 кг топлива, м3/кг; Ор - низшая теплота сгорания рабочей массы газифицируемого твердого топлива; Онг - низшая теплота сгорания влажного газа, кДж/м3, определяемая по известной зависимости
Онвг = 127,7С0 + 107,6Н2 + 356,7СН4.
Поскольку метан в окислительно-восстановительной зоне 1 не образуется из-за высоких температур, уравнение теплового баланса запишется в виде
ОР П
= 127,7С0 + 107,6Н 2. (7)
Низшая теплота сгорания рабочей массы газифицируемого твердого топлива, как и в [2, 5, 8-10], определялась в кДж/кг по формуле Д. И. Менделеева
ОнР = 340Ср + 1035НР -109(0Р -Sp)-25WР.
Таким образом, выход влажного газа из 1 кг рабочего топлива и состав генераторного газа на выходе из окислительно-восстановительной зоны 1 определяются из совместного решения уравнений (1)-(7).
Метод расчета концентраций продуктов химических реакций
по константам равновесия этих реакций
в активных зонах газификации с различными температурами
Обозначения (в соответствии с табл. 2): /,] - номер реакции и активной зоны; ] = 1 соответствует окислительно-восстановительная зона 1 (ОВЗ); ] = 1.0 - балансовая окислительно-восстановительная подзона 1.0 (БОВП),
ив.г
рассматриваемая в модифицированном балансовом методе; ] = 1.1 - подзона дополнительного окисления 1.1 (ПДО); ] = 2 - восстановительная зона 2 (ВЗ); ] = 3 - зона взаимодействия продуктов между собой и углеродом топлива 3 (ЗВП); Х - общее обозначение компонента генераторного газа; [Х\у -конечная мольная доля компонента Х, в % на 100 кмоль генераторного газа (далее - концентрация), создаваемая в результате протекания 1-й реакции в ]-й зоне; №Х)ц - начальная мольная доля (далее - концентрация) компонента Х, в % на 100 кмоль генераторного газа, для 7-й реакции в ]-й зоне; - доля 7-й реакции из системы химических превращений в ]-й зоне, доли ед.; К7]- - константа равновесия 7-й реакции в ]-й зоне; [ХЪ,/ - мольная доля (далее - концентрация) компонента Х, в % на 100 кмоль генераторного газа, создаваемая в результате протекания всех реакций в]-й зоне.
1. Окислительно-восстановительная зона 1
ОВЗ в общем случае состоит из двух подзон: балансовой окислительно-восстановительной подзоны 1.0 (БОВП) и подзоны дополнительного окисления 1.1 (ПДО), существующей только при а > 1.
Изложенный выше модифицированный балансовый метод рассматривает процессы перехода элементов твердого топлива (древесины) и газифицирующего агента (влажного воздуха) в компоненты генераторного газа только в БОВП без учета дополнительного окисления углерода и образующихся горючих газов СО и Н2 в случае избыточного воздуха (а > 1).
При а > 1 в составе генераторного газа на выходе из подзоны 1.0 присутствует неизрасходованный кислород 02з6 . В данной работе полагается, что этот неизрасходованный кислород является причиной дополнительных реакций окисления, локализованных в ПДО. Тогда с учетом наличия подзоны 1. 1 конечная концентрация СО на выходе из ОВЗ в общем случае запишется как:
[coL =
[CO]0,1.0 при
(о)
[CO]E11 при а > 1.
Концентрации [С02 ]21, [Н2 ] [Н 20 ] определяются аналогично (8)
[ N 2^1 =100-( [С°],,1 + [С02 ] + [Н2 ] + [ Н20] ).
1.1. Подзона дополнительного окисления избыточным кислородом при а > 1
Реакция 1 ПДО: С + 02 ^ С02.
Выражение константы равновесия химической реакции 1 и ее численное значение [4, 10-14] запишутся в следующем виде:
К [С02 ]■•■ ■ ■ (9)
Ки1 ^^^ТГГ' (9)
20582 8 т
lg К11Л =-^— 0,302 lg Т + 0,143 • 10-3 Т - 0,24 •Ю-7 Т2 + 0,622;
[002 ]ил =(^002 )ил + а [02 ] ил =(^02 )
- а,
1,1.1
где а - количество О2, прореагировавшее в реакции 1 ПДО, в % от 100 кмоль газа[15].
Избыточный кислород в ПДО распределяется между четырьмя протекающими реакциями в соответствии с их приоритетом, т. е. пропорционально величине югдл. Тогда начальная концентрация О2 для реакции 1 в ПДО
(^02 ) 11Л = [02 ] 0
] 0,1.0 ®1,1.1.
Начальная концентрация СО2 для реакции 1 также берется с учетом доли этой реакции в системе химических превращений ПДО
(Я002 ) 1 1Л = [002 ] 0,1,
-2/ 1,1.1 ь '2.10,1.0 ю 1,1.1 -
Результатом решения (9) является искомая величина а. Реакция 2 ПДО: 20 + 02 ^ 200.
Выражение константы равновесия химической реакции 2 и ее численное значение [4, 10-14] запишутся в следующем виде:
[00]2
К -=кЬ (10)
^ К 2,1.1 = 116г35,1 + 2,1656lg Т - 0,94 •Ю-3 Т - 0,876 •Ю-6 Т2 + 3,394;
[°°] 2,1.1 =(Яс0 ) 2,1.1 + 2Ь; [02 ]2,1.1 =(Я02) 2ДЛ - Ь
где Ь - количество О2, прореагировавшее в реакции 2 ПДО, выраженное в % от 100 кмоль газа [15];
(^2 ) 2,1.1 =[02 ] 0,1.0 Ю 2,1.1.
Начальная концентрация СО для реакции 2 также берется с учетом доли реакции 2 в системе химических превращений ПДО
(^0 ) 2,1.1 =[СО] 0,1.0 ® 2,1.1. Результатом решения (10) является искомая величина Ь.
Реакция 3 ПДО: 2СО + 02 ^ 2С02.
Выражение константы равновесия химической реакции и ее численное значение [4, 10-14] запишутся в следующем виде:
К 3,1 = [С°2]- ; (1,)
' № 3,11 '[02 ] 3,1.1
29530 5
lgК31Л =-^--2,769lgТ +1,225 • 10-3Т - 0,1356 • 10-6Т2 - 2,15;
[С0]3,1.1 = (^С )„.! -[02 ]3,1.1 =К),,! - й;
[С02 ] 3,1.1 =( V ) 3,1Л + 2й,
где й - количество О2, прореагировавшее в реакции 3 ПДО, выраженное в % от 100 кмоль газа [15\;
(Ro2 ) 31Л = [O2 ] 0,1.0 ® 3,1.1.
Начальные концентрации всех участвующих в реакции 3 компонентов газа также берутся с учетом доли реакции 3 в системе химических превращений ПДО:
(^СС ) 3,1.1 =[С0] 0,1.0 ю 3,1.1; (^ ) 3,1Л = [С02 ] 0,1.0 « 3,1.1.
Результатом решения (11) является искомая величина й. Реакция 4 ПДО: 2Н2 + 02 -о- 2Н20.
Выражение константы равновесия химической реакции и ее численное значение [4, 10-14] запишутся в следующем виде:
К ^4,1-1 . (12)
К4,1.1 = 7-~2-----; (12)
[Н2 ] 4,1.1 ^ [02 ] 4,1.1
lgК41Л = 251^16,1 -0,9466 ^Т- 7,216•Ю-4Т +1,618•Ю-7 Т2 -1,714;
[Н 2 ] 4,1.1 =( ^ ) 4ДЛ - 2е; [02 ] 4,1.1 = ( ^ ) 4,1.1 - в; [Н20] 4,1.1 = Ко ) ^ +
где е - количество О2, прореагировавшее в реакции 4 ПДО, выраженное в % от 100 кмоль генераторного газа [15\;
(^ )4,1.1 = [02 ]0,1.0 Ю4; (^ )4,1.1 = [Н2 ]0,1.0 «4 ; (^О )4Д Л = [Ц*0]^ «4.
Результатом решения (12) является искомая величина е.
Вычисление итоговой конечной концентрации СО на выходе из ПДО производится по следующей зависимости:
[ССки = №.0-2 (^с )т +1 [Ш],,,
/¿,1.1
¿=1 ¿ =1 Аналогично вычисляются [СС2 ]Е 11, [Н2 ]Е 11 и [Н2С]Е 11
[N2 к,1 = 100 -([СС]х,1Л + [СС2 ]ЕД1 + [Н2 ]ЕДЛ + [Н2С]ЕДЛ ).
2. Восстановительная зона Реакция 1 ВЗ: СС2 + Н2 ^ СС + Н2С.
Выражение константы равновесия химической реакции и ее численное значение [4, 10-14] запишутся в следующем виде:
К1,2 = [СС2 ]и-[Нг ],,2 ' ( '
2207 2
lgК12 =--^ + 0,9115 ^Т-9,738-10-4Т +1,487-10-7Т2 + 0,098;
[СС2 1,2 =()12 - I; [Н2 ]1,2 = (^Н2 )12 -I;
[СС1,2 =(^СС )1,2 + I; [Н2С]и =(^2С )1,2 + I,
где I — количество Н2, прореагировавшее в реакции 1 ВЗ, выраженное в % от 100 кмоль генераторного газа [15].
Начальные концентрации всех компонентов газа, участвующих в ¿-й реакции ВЗ, берутся на выходе из ОВЗ с учетом доли данной реакции в системе химических превращений ВЗ, т. е. пропорционально величине «¿,2:
Кс2 )и = [СС2 ],, «1,2; К )и = [Н 2 ],д «1,2 ;
(^С )1,2 =[СС]Е,1 «1,2; (*Н2С )1,2 = [Н2О],д «1,2.
Результатом решения (13) является искомая величина I Реакция 2 ВЗ: С+СС2 ^ 2СС.
Выражение константы равновесия химической реакции и ее численное значение [4, 10-14] запишутся в следующем виде:
[СС]2 2
К22 = у-^; (14)
2 2 [СС2 ]2,2
8947 7
lgK2, 2 =--2,4673 lgT-1,0824•Ю-3T + 0,116•Ю-6T2 + 2,772;
[C02 ] 2,2 = (^C°2 ) 2 2 - g; [C0]2,2 = (*C° )2,2 + 2
где g - количество СО2, прореагировавшее в реакции 2 ВЗ, выраженное в % от 100 кмоль генераторного газа [15];
(Rc°2 )2,2 = [C02 L,1 «2,2; (*C° )2,2 = [C0L,1 ® 2,2.
Результатом совместного решения (14) является искомая величина g. Реакция 3 ВЗ: C + H20 ^ CO + H2.
Выражение константы равновесия химической реакции и ее численное значение [4, 10-14] запишутся в следующем виде:
K '"I IH2. ,|5)
32 = [Нг0]з,22 ' ('5)
lg K3 2 = - 6740,5 +1,5561 lg T -1,092 •Ю-4 T - 0,37 •Ю-7 T2 + 2,554;
Н[Н20]з,2 = (RH2° )з,2 - [C°]3,2 = (R° )з,2 + [H2 ]з,2 = К ^ + h
где h - количество Н2О, прореагировавшее в реакции 3 ВЗ, выраженное в % от 100 кмоль генераторного газа [15];
(RH2° )з,2 = [Н20],д «3,2; (RC0 )3,2 = [C0]Z,! «3,2; (RH2 )З,2 = [Н2 L,1 «3,2 •
Результатом решения (15) является искомая величина h. Реакция 4 ВЗ: C + 2Н 20 ^ C02 + 2Н 2.
Выражение константы равновесия химической реакции и ее численное значение [4, 10-14] запишутся в следующем виде:
K [С°2 ]4,2 •[Н2 ]4,2 . (16)
K4,2 =---12-; (16)
[Н 20]4,2
4533 3
lg K42 =--+ 0,6446lg T + 8,616 •Ю-4 T -1,858-10-7 T2 + 2,336;
[H2°]4,2 = (RH2° )4,2 - 2 j; [C02 ]4,2 = (V 1,2 + j; [H2 ^,2 = (+ 2 j
где 2j - количество Н2О, прореагировавшее в реакции 4 ВЗ, выраженное в % от 100 кмоль генераторного газа [15];
(Ян2о )4,2 = [Н20]ЕЛ ®4,2; ()4,2 = [С02 ^ ®4,2; (^ \2 = [Н2 ^ ®4,2.
Результатом решения (16) является искомая величина Для вычисления конечного значения компонента генераторного газа СО на выходе из восстановительной зоны использована следующая зависимость:
№2 = №: -I (^0 )■ 2 + 1 [С0],2.
7=1 ' 7 =1
Аналогично вычисляются [С02 ]Е 2, [Н2 ]Е 2, [Н20]Е 2
[ N2 ]х,2 = 100 "(№,2 + [С02 Ц + [ Н2 Ь,2 + [Н 2 0]2,2 ) •
3. Зона взаимодействия продуктов между собой и с углеродом топлива
Реакция 1 ЗВП: С0 + Н20 ^ С02 + Н2.
Выражение константы равновесия химической реакции и ее численное значение [4, 10-14] запишутся в следующем виде:
= [С0]1,3 .[Н20]1,3; (17)
2207 2
^К1з = —у— 0,9115lgТ + 9,738-10-4Т-1,487-10-7 Т2 - 0,098;
[С0]1,з =(^0 )1,з " [Н2°]и = (^0 )1,з" к; [С02 ]1,3 = ( V )1,з + к; [Н2 ]1,3 = (^ )1,з + к.
где к - количество Н2О, прореагировавшее в реакции 1 ЗВП, выраженное в процентах от 100 кмоль генераторного газа [15].
Начальные концентрации всех компонентов 7-й реакции ЗВП берутся на выходе из ВЗ с учетом доли данной реакции в системе химических превращений ЗВП, т. е. пропорционально величине ю7,з:
(^0 )1,з = [С0]Е,2 ®1,3; (*Н20 )и = [Н20]^,2 ®1,з;
(^ )и = [С02 ]Е,2 ®1,з; К )1,з = [Н2 ]Е,2 ®1,з.
Результатом решения (17) является искомая величина к.
Реакция 2 ЗВП: С + 2Н2 ^ СН4.
Выражение константы равновесия химической реакции и ее численное значение [4, 10-14] запишутся в следующем виде:
[СН 4 ] 23
к2,3 = -—"Т1; (18)
[Н2 ] 2,3
3348
^к2 3 = —--5,57 lgТ + 0,186-10-2 Т-1,095-10-7 Т2 +11,79;
[Н 2 ] 2,3 = К ) 23 - 2т'; [СН4 ]2,3 = Кн )2 3 + ^
где 2m - количество Н2, прореагировавшее в реакции 2 ЗВП, выраженное в % от 100 кмоль генераторного газа [15];
К )23 = [Н2 ]Е,2 «2^ (Rh< )23 = [СН4 ]
2,3 L ^1,2 ®2,3.
Поскольку метан СН4 отсутствует среди компонентов генераторного газа в ОВЗ и ВЗ, то () = 0.
Результатом решения (18) является искомая величина т. Реакция 3 ЗВП: СО+3Н2 ^ СН4 + Н20.
Выражение константы равновесия химической реакции и ее численное значение [4, 10-14] запишутся в следующем виде:
Кз,з=1СНЫНф; (19)
[co] 3,3 ■[ h2 ]
lgK33 = 7,14 lgT + 0,188■ÍO-2T-0,094■ÎO-6T1 + 8,64;
9874 2 6 2
--7,14 lgT + 0,188-10-2^ пппл
T
[H2 ]з,з = (RH )з,з - 3n; [C°]3,3 = ()з,з - n;
[CH4 ]3,3 = (rch4 )3 3 + n; [н2°]з,з = (RH2° )3 3 + n
где 3n - количество Н2, прореагировавшее в реакции 3 ЗВП, выраженное в % от 100 кмоль генераторного газа [15];
(RC° )3,3 = [C0]s,2 «3,3; (RH2 )зз =[ Н2 ]z,2 «3,3 ;
(RCH4 )33 = [СН4 ]Е,2 «3,3; (RH2O ) 33 = [Н20]
33 L 2 JZ,2 « 3,3.
Поскольку метан СН4 отсутствует среди компонентов генераторного газа в ОВЗ и ВЗ, то (ЛСН ) = 0.
Результатом решения (19) является искомая величина п.
Реакция 4 ЗВП: С02 + 4Н2 ^ СН4 + 2Н20.
Выражение константы равновесия химической реакции и ее численное значение [4, 10-14] запишутся в следующем виде:
п2
[СН4 ]4 3 -[Н 2 0]4з
К = у--мз; (20)
[СО]4,з -[Н2]4,з
7674 , 7 т
lg К43 = —--6,23 ^ Т + 0,906 -10-3 Т + 0,596 -10-7 Т2 + 8,79;
[Н2 ]4,3 = (^ ^ - 4Р; [СО2 ]4,3 = («СО2 ^ " Р;
[СН4 ],3 = (^ ^ + Р; [Н2О],3 = (V ^ + 2А
где 4р - количество Н2, прореагировавшее в реакции 4 ЗВП, выраженное в % от 100 кмоль генераторного газа [15];
(КСО2 )4,3 = [СО2 ]Е,2 «4,3 ; (^ )4,3 = [Н2 ]е,2 «4,3 ; (^ )4,3 = [СН4 ]е,2 «4,3; (^О ^ = [Н2О]Е,2 «4,3 .
Поскольку метан СН4 отсутствует среди компонентов генераторного газа в ОВЗ и ВЗ, то (ЛСН4) = 0.
Результатом решения (20) является искомая величина р. Для вычисления конечного значения компонента генераторного газа СО на выходе из зоны взаимодействия продуктов между собой и углеродом топлива использовано следующее выражение:
№,3 =№,2-I (^О ) г3 +1 [С°] г,3.
г=1 ' г=1
Аналогично вычисляются [СО2 ]Е 3, [Н2 ]Е 3, [Н2О]Е 3
[N2 ]е,3 = 100 - ([СО]^,3 + [СО2 ]е,3 + [Н2 ]Е,3 + [Н2О]Е,3 ) .
Результаты расчетов по разработанному методу сравнивались с данными газового анализа сухого состава генераторного газа и приведены в табл. 3. Сравнение расчетных и опытных данных осуществлялось при одинаковых значениях содержания азота К2, не участвующего ни в одной из химических реакций. Данные опытов Д. Б. Гинзбурга [3] по газификации березы влажностью 28 % (табл. 1) сопоставлялись с результатами расчетов по разработанному методу при коэффициенте избытка воздуха а = 0,33. Получено приемлемое совпадение расчетных и опытных значений состава генераторного газа по всем его компонентам.
Таблица 3
Сухой состав генераторного газа из березы влажностью 28 % по разработанному методу*
The dry composition of the generator gas produced from birch with moisture 28 %
calculated by this method
а СО, % СО2, % Н2, % О2, % СН4, % N2, % QH, кДж/м3
0,1 35,36 10,79 14,45 0,00 6,40 33,00 8353
0,2 30,88 9,18 12,44 0,00 5,61 41,89 7282
0,3 27,41 7,96 10,87 0,00 5,01 48,75 6456
0,31 27,10 7,86 10,73 0,00 4,96 49,36 6383
0,32 26,81 7,76 10,60 0,00 4,90 49,93 6313
0,33 26,52 7,66 10,47 0,00 4,85 50,50 6244
0,34 26,23 7,56 10,33 0,00 4,81 51,07 6176
0,35 25,95 7,47 10,21 0,00 4,76 51,61 6110
0,4 24,65 7,03 9,62 0,00 4,54 54,17 5801
0,5 22,40 6,28 8,60 0,00 4,16 58,57 5269
0,6 20,49 5,60 7,53 0,00 3,89 62,49 4813
0,7 18,91 5,10 6,81 0,00 3,63 65,55 4442
0,8 17,27 4,65 5,60 0,00 3,58 68,90 4083
0,9 16,15 4,29 5,16 0,00 3,36 71,04 3816
1,0 15,10 3,98 4,61 0,00 3,22 73,09 3573
1,2 15,66 5,03 4,40 0,00 2,53 72,39 3375
1,4 15,61 5,14 3,25 0,00 2,27 73,73 3154
1,6 15,98 5,20 2,62 0,00 2,11 74,09 3077
1,8 16,34 5,38 2,56 0,00 1,89 73,83 3036
2,0 16,62 5,46 2,32 0,00 1,74 73,86 2993
3,0 17,43 5,54 1,07 0,00 1,39 74,57 2836
"Температура воздуха 283 К, относительная влажность воздуха 70 %. Температуры по зонам [5]: окислительно-восстановительная 1500 К, восстановительная 1275 К, зона взаимодействия продуктов друг с другом и углеродом топлива 950 К.
ВЫВОДЫ
1. Предложен новый метод расчета состава древесного генераторного газа, представляющий совокупность модифицированного общепринятого балансового метода [2] и добавленного метода расчета концентраций продуктов химических реакций по константам равновесия этих реакций в активных зонах газификации с различными температурами. В данном методе различают не одну, а три активные зоны газификации: окислительно-восстановительную, восстановительную и зону взаимодействия продуктов газификации между собой и с углеродом вновь поступающего нагретого топлива. В материальном и тепловом балансе устранены основные недостатки общепринятого метода [2].
2. По величине изменения свободной энтальпии (энергии Гиббса) ряда химических реакций процесса газификации древесины проведена оценка
химического сродства реагентов, выявлены преобладающие (основные) химические реакции, установлен их приоритет в конкретных температурных условиях (активных зонах) и найдены доли этих реакций в системе химических превращений.
3. Метод расчета концентраций продуктов химических реакций по константам равновесия этих реакций позволяет определять концентрации продуктов основных химических реакций в каждой из активных зон по начальным концентрациям исходных веществ или их соединений (газов), температуре и константам равновесия этих реакций.
4. Разработанный метод расчета компонентов генераторного газа, состоящий из балансового метода (материальный и тепловой балансы) и метода расчета концентраций продуктов химических реакций по константам равновесия этих реакций, позволяет более достоверно моделировать процесс газификации древесного топлива в газогенераторе обращенного процесса и получать приемлемое согласование расчетных и опытных данных по составу генераторного газа.
ЛИТЕРАТУРА
1. Газогенератор твердого топлива: пат. 2579112 РФ, МПК F23G5/027, СШ3/00 / Е. М. Кашин, В. Н. Диденко; опубл. 27.03.2016.
2. Юдушкин, Н. Г. Газогенераторные тракторы. Теория, конструкция и расчет / Н. Г. Юдуш-кин. М.: Гос. науч.-техн. изд-во машиностроит. лит-ры, 1955. 244 с.
3. Гинзбург, Д. Б. Газификация низкосортного топлива / Д. Б. Гинзбург. М.: Гос. изд-во лит-ры по строит. матер., 1950. 173 с.
4. Федосеев, С. Д. Полукоксование и газификация твердого топлива / С. Д. Федосеев, А. Б. Чернышев. М.: Гос. науч.-техн. изд-во нефтян. горно-топлив. лит-ры, 1960. 327 с.
5. Коллеров, Л. К. Газомоторные установки / Л. К. Коллеров. М., СПб.: Гос. науч.-техн. изд-во машиностроит. лит-ры, 1951. 240 с.
6. Лямин, В. А. Газификация древесины / В. А. Лямин. М.: Лесная промышленность, 1967. 263 с.
7. Равдель, А. А. Краткий справочник физико-химических величин / А. А. Равдель, А. М. Пономарева. 8-е изд., перераб. Л.: Химия, 1983. 232 с.
8. Головков, С. И. Энергетическое использование древесных отходов / С. И. Головков, И. Ф. Коперин, В. И. Найденов. М.: Лесная промышленность, 1987. 221 с.
9. Токарев, Г. Г. Газогенераторные автомобили / Г. Г. Токарев. М.: МАШГИЗ, 1955. 204 с.
10. Дешалит, Г. И. Расчеты процессов газификации топлив / Г. И. Дешалит. Харьков: Изд-во Харьковского Ордена Трудового Красного Знамени гос. ун-та имени А. М. Горького, 1959. 167 с.
11. Богданов, Н. Н. Полукоксование и газификация торфа / Н. Н. Богданов. М.; Л.: Гос. энергетич. изд-во, 1947. 271 с.
12. Лавров, Н. В. Физико-химические основы горения и газификации топлива / Н. В. Лавров. М.: Гос. науч.-техн. изд-во по черной и цветной металлургии, 1957. 289 с.
13. Лавров, Н. В. Введение в теорию горения и газификации топлива / Н. В. Лавров, А. П. Шуригин. М.: Изд-во Академии наук СССР, 1962. 217 с.
14. Совершенствование техники и технологии процесса газификации высоковлажных древесных отходов / Н. Ф. Тимербаев [и др.]. Казань: Изд-во КНИТУ, 2013. 92 с.
15. Иноземцев, Н. В. Основы термодинамики и кинетики химических реакций / Н. В. Иноземцев. М.: Военная Академия механизации и моторизации Красной Армии имени Сталина, 1940. 257 с.
Поступила 18.05.2018 Подписана в печать 27.02.2019 Опубликована онлайн 31.07.2019
REFERENCES
1. Kashin E. M., Didenko V. N. (2016) The Gasifier of Solid Fuel. Patent No 2579112 of the Russian Federation (in Russian).
2. Yudushkin N. G. (1955) Gas Generator Tractors. Theory, Design, Calculation. Moscow, State Scientific-and-Technical Publishing House of Mechanical Engineering Literature. 244 (in Russian).
3. Ginzburg D. B. (1950) The Gasification of Low-Grade Fuel. Moscow, State Publishing House of Literature on Construction Materials. 173 (in Russian).
4. Fedoseev S. D., Chernyshev A. B. (1960) Semi-Coking and Gasification of Solid Fuel. Moscow, State Scientific-and-Technical Publishing House of Oil and Mine Fuel Literature. 327 (in Russian).
5. Kollerov L. K. (1951) Gas Engine Installations. Sankt-Petersburg, State Scientific-and-Technical Publishing House of Mechanical Engineering Literature. 240 (in Russian).
6. Lyamin V. A. (1967) Wood Gasification. Moscow, Lesnaya Promyshlennost' Publ. 263 (in Russian).
7. Ravdel' A. A., Ponomareva A. M. (1983) Concise Reference Book of Physical and Chemical Magnitudes. 8th ed. Leningrad, Khimiya Publ. 232 (in Russian).
8. Golovkov S. I., Koperin I. F., Naidenov V. I. (1987) Energy Use of Wood Waste. Moscow, Lesnaya Promyshlennost' Publ. 221 (in Russian).
9. Tokarev G. G. (1955) Gas Generator Automobiles. Moscow, MASHGIZ Publ. 204 (in Russian).
10. Deshalit G. I. (1959) The Calculations of the Fuel Gasification Processes. Kharkov, Publishing House of the Kharkiv Labour Red Banner Order State University named after A. M. Gor-ky. 167 (in Russian).
11. Bogdanov N. N. (1947) Semi-Coking and Gasification of Peat. Moscow, Leningrad, State Power Engineering Publ. 271 (in Russian).
12. Lavrov N. V. (1957) Physical and Chemical Bases of Combustion and Gasification of Fuel. Moscow, State Scientific-and-Technical Publishing House of Ferrous and Nonferrous Metallurgy. 289 (in Russian).
13. Lavrov N. V., Shurigin A. P. (1962) Introduction to the Theory of Combustion and Gasification of Fuel. Moscow, Publishing House of USSR Academy of Sciences. 217 (in Russian).
14. Timerbaev N. F., Safin R .G., Hisameeva A. R., Ishakov T. D. (2013) Improvement of the Equipment and Technology of Gasification Process of High-Moist Wood Waste. Kazan, KNITU Publ. 92 (in Russian).
15. Inozemtsev N. V. (1940) Basics of Thermodynamics and Kinetics of Chemical Reactions. Moscow, Military Red Army Mechanization and Motorization Academy named after Stalin. 257 (in Russian).
Received: 18 May 2018
Accepted: 27 February 2019
Published online: 31 July 2019
