посмотрите пример расчёта от котельных работающих на природном газе. всё зависит от состава газа.
Программа ЭкоРасчета (версия 4.12 НПП Логус Москва
КОТЕЛЬНЫЕ
=========
Расчетные алгоритмы модуля основаны на нормативных материалах,
заложенных в "Методике определения выбросов загрязняющих веществ"
в атмосферу при сжигании топлива в котлах производительностью
менее 30 тонн пара в час или менее 20 ГКал в час", Москва, 1999г.
Расчетные формулы
-----------------
Сжигание природного газа
а)Оксиды азота NOx
M(NOx) = Bp (Bp) * Qr * Kr * bk * bt * ba * (1-bг) * (1-bд) * kп, г/с (т/год)
где:
М(NOx) - выброс оксидов азота NOx
Bp - расчетный расход топлива (куб.нм/с) для расчета выбросов г/с
Bp - расчетный расход топлива (тыс.куб.нм/год) для расчета выбросов т/год
Bp (Bp) = B (B) * (1-q4/100),
где B - фактический расход топлива, куб.нм/с
B - фактический расход топлива, тыс.куб.нм/год
При работе котла в соответствии с режимной картой Bp (Bp) = B (B)
Qr - низшая теплота сгорания топлива, МДж/куб.нм
Kr - удельный выброс оксидов азота, г/МДж
Kr = 0.01 * D^0.5 + 0.03 - для паровых котлов
где D - паропроизводительность котла (т/ч):
- средняя - для расчета выбросов т/год
- максимальная - для расчета выбросов г/с
Kr = 0.0113 * (Bp *Qr)^0.5 + 0.03 - для водогрейных котлов
bk - коэффициент, учитывающий конструкцию горелки
bk = 1 - для дутьевых горелок напорного типа
bk = 1.6 - для горелок инжекционного типа
bk = 0.7 - для горелок двухступенчатого сжигания
bt - коэффициент, учитывающий температуру подаваемого для горения воздуха
bt = 1 + 0.002 * (tгв - 30)
tгв - температура подаваемого для горения воздуха
ba - коэффициент, учитывающий влияние избытка воздуха
ba = 1 - при работе котла в соответствии с режимной картой
ba = 1 - 0.1 * (О2 - 5/Q_)^2 - 0.3 * (О2 - 5/Q_) - для дутьевых
горелок напорного типа и горелок двухступенчатого сжигания
при наличии данных о концентрации кислорода в дымовых газах
О2 - концентрация кислорода в дымовых газах за котлом, %
Q_ = Qф/Qн или Q_ = Dф/Dн
Qф, Dф, Qн, Dн - соответственно фактические и номинальные
тепловая нагрузка и паропроизводительность котла, МВт (т/ч)
ba = 0.577 * Sт^0.5 - для котлов с инжекционными горелками
Sт - разрежение в топке, кгс/кв.м (мм вод.ст.)
ba = 1.225 - в общем случае
bг - коэффициент, учитывающий влияние рециркуляции дымовых газов
через горелки
bг = 0.16 * r^0.5
r - степень рециркуляции дымовых газов, %
bд - коэффициент, учитывающий ступенчатый ввод воздуха в топочную
камеру
bд = 0.022 * д
д - доля воздуха, подаваемого в промежуточную зону факела, (% от
общего количества организованного воздуха)
kп - коэффициент пересчета
kп = 1 - при расчете выбросов г/с
kп = 0.001 - при расчете выбросов т/год
Примечание.
Нормирование выбросов оксидов азота с учетом их трансформации
в атмосферном воздухе в оксид и диоксид азота производится с
использованием экспериментально определенных коэффициентов
трансформации, а в случае отсутствия экспериментальных данных -
в соответствии с действующими нормативными документами.
б)Оксид углерода
M(CO) = B * Cco * (1 - q4/100), г/с
M(CO) = 0.001 * B * Cco * (1 - q4/100), т/год
где:
М(CO) - выброс оксида углерода
B - фактический расход топлива, куб.нм/с
B - фактический расход топлива, тыс.куб.нм/год
Cco - выход оксида углерода при сжигании топлива, г/куб.нм
q4 - потери тепла вследствие механической неполноты сгорания топлива, %
Cco = q3 * R * Qr, где
q3 - потери тепла вследствие химической неполноты сгорания
топлива, %
R = 0.5 - коэффициент, учитывающий долю потери тепла вследствие
химической неполноты сгорания топлива
в)Оксиды серы
M(SO2) = 0.02 * B (B) * pг * Sr * (1-n so2) * (1-n so2), г/с (т/год)
где:
М(SO2) - выброс оксидов серы
B - фактический расход топлива (куб.нм/с)*1000 для расчета выбросов г/с
B - фактический расход топлива (тыс.куб.нм/год) для расчета выбросов т/год
pг - плотность газа, кг/куб.нм
Sr = srab (srab) + 0.94 * H2S + Cs
Содержание серы в топливе на рабочую массу, %
srab - максимальное за год для расчета выбросов г/сек
srab - среднее за год для расчета выбросов т/год
H2S - содержание сероводорода в топливе на рабочую массу, %
Cs - содержание меркаптановой серы в топливе на рабочую массу, %
n so2 - доля оксидов серы, связываемых летучей золой в котле
(для газа = 0)
n so2 - доля оксидов серы, улавливаемых в мокром золоуловителе
- средняя за год для расчета выбросов т/год
- максимальная за год для расчета выбросов г/сек
г)Бенз(а)пирен
M(бенз(а)пирен) = с * Vсг * Bp (Вр) * kп, г/с (т/год)
где:
М(бенз(а)пирен) - выброс бенз(а)пирена
Bp - расчетный расход топлива (тыс.куб.нм/час) для расчета выбросов г/с
Bp - расчетный расход топлива (тыс.куб.нм/год) для расчета выбросов т/год
Bp = B * (1-q4/100) * 3.6
Bp = B * (1-q4/100)
где B - фактический расход топлива, куб.нм/с
B - фактический расход топлива, тыс.куб.нм/год
с = с*alfa/1.4
с - концентрация бенз(а)пирена в сухих продуктах сгорания природного
газа, мг/куб.нм
Для промтеплоэнергетических котлов малой мощности:
- при alfa = 1.08 - 1.25
c = (0.059 + 0.079*0.001*qv)*Kд*Kр*Kст / [1000*е^(3.8*(alfa - 1))]
- при alfa > 1.25
c = (0.032 + 0.043*0.001*qv)*Kд*Kр*Kст / [1000*е^(1.14*(alfa - 1))]
Для водогрейных котлов:
- при alfa = 1.05 - 1.25 и qv = 250 - 500 кВт/куб.м
c = (0.11*qv - 7)*Kд*Kр*Kст / [1000000 * е^(3.5*(alfa - 1))]
- при alfa > 1.25 и qv = 250 - 500 кВт/куб.м
c = (0.13*qv - 5)*Kд*Kр*Kст / [1000000*1.3*е^(3.5*(alfa - 1))]
где alfa - коэффициент избытка воздуха в продуктах сгорания на
выходе из топки
qv - теплонапряжение топочного объема (кВт/куб.м), при
сжигании непроектного топлива qv = Br*Qr/Vт, где
Br - расчетный расход топлива на номинальной нагрузке, куб.нм/с
Br = Br * (1 - q4/100),
Br - фактический расход топлива на номинальной нагрузке
котла, куб.нм/с
Qr - низшая теплота сгорания топлива, кДж/куб.нм
Vт - объем топочной камеры, куб.м
Kд - коэффициент, учитывающий влияние нагрузки котла на
концентрацию бенз(а)пирена в продуктах сгорания
Kр - коэффициент, учитывающий влияние рециркуляции дымовых
газов на концентрацию бенз(а)пирена в продуктах сгорания
Kст - коэффициент, учитывающий влияние ступенчатого сжигания
на концентрацию бенз(а)пирена в продуктах сгорания
Vсг - объем сухих дымовых газов, образующихся при полном сгорании
1 куб.нм топлива
Vсг = Voг + (alfa - 1) * Vo - Voh2o,
где:
Vo - объем воздуха при стехиометрическом сжигании 1 куб.нм
топлива, куб.нм/куб.нм
Voг - объем дымовых газов при стехиометрическом сжигании
1 куб.нм топлива, куб.нм/куб.нм
Voh2o - объем водяных паров при стехиометрическом сжигании
1 куб.нм топлива, куб.нм/куб.нм
kп = 10e-6 - для определения выбросов т/год
kп = 0.278*10e-3 - для определения выбросов г/c
д) Определение объемного расхода дымовых газов
Vрг = B * [k1 + k2*Qr + (alfa-1)*(k3 + k4*Qr)] * (273 + tp) / 273, куб.м/с
где:
B - секундный расход топлива, куб.нм/с
Qr - низшая теплота сгорания топлива, МДж/куб.нм
alfa - коэффициент избытка воздуха
tр - температура дымовых газов, град С
ki - численные коэффициенты, подобранные методом наименьших квадратов
k1 = 0.739
k2 = 0.278
k3 = 0.0864
k4 = 0.267
-----------------------------------------------------------------------------
Наименование источника выделения : Котёл
Номер источника выделения : 101
Характеристики топлива
----------------------
Топливо : природный газ
Газопровод : Игрим-Пунга-Серов-Нижний Тагил
Плотность сухого газа, кг/нм3: 0.746
Содержание сероводорода в топливе (%) : 0.0000
Содержание в топливе меркаптановой серы (% мас) 0.0000
Сернистость (%)
- среднегодовая : 0.02
- максимальная : 0.02
Топливо : проектное
Расход топлива :
- при максимальной нагрузке котла, куб. нм/с : 0.05000
- годовой, тыс. куб. нм/год : 5000.000
Эксплуатационные характеристики котлоагрегата
---------------------------------------------
Тип котла : паровой
Работа котла в соответствии с режимной картой
Производительность котлоагрегата (т/ч)
- средняя фактическая : 1.00000
- максимальная фактическая : 1.00000
- номинальная: 1.00000
Тип топки котла : КАМЕРНАЯ ТОПКА
Тип горелок котла :
Теплонапряжение топочного объема, кВт/куб.м : 0.0
Степень рециркуляции дымовых газов, % : 0.00
Температура подаваемого для горения воздуха, °С : 20
Доля воздуха, подаваемая в промежуточную зону факела, % : 0.00
Наличие золоуловителя : Отсутствует
Коэф.избытка воздуха в продуктах сгорания на выходе из топки : 1.11
Температура дымовых газов, °С : 670
Коэффициенты трансформации оксидов азота :
- в диоксид азота :
- для расчета выбросов т/год 0.800
- для расчета выбросов г/сек 0.800
- в оксид азота :
- для расчета выбросов т/год 0.130
- для расчета выбросов г/сек 0.130
----------------------------------------------------------------------
Потери тепла от механической неполноты сгорания
q4=0.1 %
Низшая теплота сгорания топлива
Qr=36.47 МДж/куб.нм
Удельный выброс оксидов азота для расчета выбросов т/год
Kr=0.01*1^0.5+0.03=0.04 г/МДж
Удельный выброс оксидов азота для расчета выбросов г/сек
Kr=0.01*1^0.5+0.03=0.04 г/МДж
Коэффициент, учитывающий конструкцию горелки
bk=0.7
Коэффициент, учитывающий температуру подаваемого для горения воздуха
bt=1
Коэффициент, учитывающий влияние избытка воздуха
ba=1
Коэффициент влияния рециркуляции дымовых газов через горелки
bг=0.16*0^0.5=0
Коэффициент, учитывающий ступенчатый ввод воздуха в топочную камеру
bд=0.022*0=0
Содержание серы в топливе на рабочую массу
среднегодовое: Sr=0.02+0.94*0+0=0.02 %
максимальное: Sr =0.02+0.94*0+0=0.02 %
Доля оксидов серы, связываемых летучей золой в котле
n so2=0
Коэффициент, учитывающий долю потери тепла вследствие химической
неполноты сгорания топлива
R=0.5
Выход оксида углерода при сжигании топлива
Cco=0.2*0.5*36.47=3.647 г/куб.нм
Объем дымовых газов при стехиометрическом сжигании 1 куб.нм топлива
Voг=10.86 куб.нм/куб.нм
Объем воздуха при стехиометрическом сжигании 1 куб.нм топлива
Vo=9.68 куб.нм/куб.нм
Объем водяных паров при стехиометрическом сжигании 1 куб.нм топлива
Voh2o=2.17 куб.нм/куб.нм
Коэффициент избытка воздуха в продуктах сгорания на выходе из топки
alfa=1.11
Объем сухих дымовых газов, образующихся при полном сгорании
1 куб.нм топлива
Vсг=10.86+(1.11-1)*9.68-2.17=9.7548001 куб.нм/куб.нм
Коэффициент, учитывающий влияние рециркуляции дымовых газов на
концентрацию бенз(а)пирена в продуктах сгорания
Kр=1
Коэффициент, учитывающий влияние ступенчатого сжигания на концентрацию
бенз(а)пирена в продуктах сгорания
Kст=1
Коэффициент, учитывающий влияние нагрузки котла на концентрацию
бенз(а)пирена в продуктах сгорания
Kд=1
Концентрация бенз(а)пирена в сухих продуктах сгорания природного
газа
c =1.11/1.4*0.001*(0.059+0.000079*0)/Exp(3.8*(1.11-1))*1*1*1=0.0000308 мг/куб.нм
Определение объемного расхода дымовых газов
Vрг = 0.05*[0.739+0.278*36.47+(1.11-1)*(0.0864+0.267*36.47)]*(273+670)/273=2.06532335 куб.м/с
----------------------------------------------------------------------
Расчет выбросов ЗВ
Азота диоксид
M(NO2)=0.8*5000*(1-0/100)*36.47*0.04*0.7*1*1*(1-0)*(1-0)*0.001=4.0846400 т/год
G(NO2)=0.8*0.05*(1-0/100)*36.47*0.04*0.7*1*1*(1-0)*(1-0)=0.0408464 г/сек
Азота оксид
M(NO)=0.13*5000*(1-0/100)*36.47*0.04*0.7*1*1*(1-0)*(1-0)*0.001=0.6637540 т/год
G(NO)=0.13*0.05*(1-0/100)*36.47*0.04*0.7*1*1*(1-0)*(1-0)=0.0066375 г/сек
Оксиды серы (в пересчете на SO2)
M(SO2)=0.02*5000*0.746*(1-0)*(1-0)*0.02=1.4920000 т/год
G(SO2)=0.02*0.05*0.746*(1-0)*(1-0)*0.02*1000=0.0149200 г/сек
Оксид углерода (CO)
M(CO)=0.001*5000*3.647*(1-0.1/100)=18.2167650 т/год
G(CO)=0.05*3.647*(1-0.1/100)=0.1821677 г/сек
Бенз(а)пирен
М(бенз(а)пирен)=0.000001*5000*(1-0.1/100)*0.0000308*9.7548001=0.0000015 т/год
G(бенз(а)пирен)=0.000278*0.05*3.6*(1-0.1/100)*0.0000308*9.7548001=1.501803163e-08 г/сек
----------------------------------------------------------------------
Shinkarenko писал(а)
От котельной на природном газе. Раньше всегда считали 4 вещества от сжигания природного газа (бензапирен, диоксид азота, оксид азота,углерод оксид), а сейчас стали просить обязательно считать еще и диоксид серы. Я вообще расчеты делаю в ЭРЕ, но там такого вещества не считает, звонила в техподдержку разработчикам программы, они говорят, что от природного газа не выделяется диоксид серы, нет таких методик. Если знаете как, то досчитывайте вручную это вещество. Я нашла что-то похожее только в НИИ Атмосфера 2012, и то там все не понятно.