Учебное пособие: Тепловой расчет котлов ДКВР
Название: Тепловой расчет котлов ДКВР Раздел: Промышленность, производство Тип: учебное пособие | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
ВведениеВ руководстве изложена методика теплового расчета котлов типа ДКВР, не содержащих пароперегревателя. Некоторые величины, входящие в расчет (например, коэффициенты ослабления излучения дымовыми газами, коэффициенты теплоотдачи при обтекании конвективных пучков и др.), вычисляются с помощью программы «BoilerTools». Задание по курсовому проекту должно содержать параметры пара на выходе из агрегата и его паропроизводительность, температуру питательной воды и состав топлива. Для выполнения теплового расчета понадобятся чертежи котлоагрегата и значения ряда параметров, характеризующих процессы в котельной установке. Эти материалы опубликованы на сайте кафедры. Рекомендуется так же использовать нормативный метод [1]. Объем и состав дымовых газовДля твердого и жидкого топлива объемы воздуха и продуктов сгорания, отнесенные к одному килограмму топлива, вычисляются по приведенным ниже формулам. По составу топлива определяется минимально необходимое количество воздуха V вз о для полного сжигания топлива V вз о = 0,0889 (C p + 0,375S p )+ 0,265H p − 0,0333O pВ приведенной формуле H p ,C p ,O p и S p - состав рабочей массы топлива, выраженный в процентах. В результате полного сжигания топлива в этом количестве воздуха образуются продукты сгорания, объем которых, отнесенный к килограмму топлива, равен: азот N p o V N 2 = 0,79V вз + 0,8 100 трехатомные газы C p + 0,375S p V RO 2 = 1,866 100 водяные пары V H O o 2 = 0,111H p + 0,0124W p + 0,0161V вз oЗдесь W p - содержание влаги в рабочей массе топлива, выраженное в процентах. Для газообразного топлива аналогичные величины определяются по следующим формулам: минимально необходимое количество воздуха для полного сжигания газа V вз o = 0,0476 0,⎡ ⎢ 5CO + 0,5H 2 +1,5H S 2 + ∑ ⎛ ⎜⎝m + n 4⎞ ⎟⎠C H m n −O 2 ⎦⎥⎤⎣ содержание азота в продуктах сгорания o N 2 V N 2 = 0,79V вз + 100 содержание трехатомных газов в продуктах сгорания VRO 2 = 0,01[CO2 +CO + H S2 +∑mC Hm n ]содержание водяных паров в продуктах сгорания V H Oo 2 = 0,01⎡⎢⎣H S H 2 + 2 +∑n 2 C H m n + 0,124d T ⎥⎦⎤+ 0,0161V взoЗдесь C H m n ,H 2 ,CO и др. – состав газообразного топлива, отнесенный к кубическому метру топлива, выраженный в процентах, d T - влагосодержание газообразного топлива, отнесенное к кубическому метру топлива. При избытке воздуха α>1объем водяных паров в продуктах сгорания равен V H O 2 =V H Oo 2 + 0,0161(α−1)V взo Полный объем дымовых газов равен V Г =V RO 2 +V N 2 +V H O 2 + −(α 1)V взo It –диаграмма продуктов сгоранияЗадаются коэффициент избытка воздуха в топкеα Т и всеми присосами Δα по ходу дымовых газов. Эти параметры приведены в таблицах ниже. Коэффициент избытка воздуха
Присосы воздуха в газоходах котла
Вычисляется действительное количество воздуха V д , подаваемого в топку V д =αТ V взо Выполняется расчет энтальпий продуктов сгорания для температур в диапазоне 100 – 2000 С. Для этого удобно вначале вычислить энтальпию продуктов сгорания при α=1. I Гo = (ct )N 2V N 2 + (ct )RO 2V RO 2 + (ct )H O H O 2 V 2 Здесь (ct ) - произведение средней теплоемкости газа в интервале температур 0 - t и температуры газа t . Затем вычисляется энтальпия продуктов сгоранияI Г при заданной температуре и коэффициенте избытка воздуха α по формуле I Г = I Гo + (α−1)(ct )вз вз V о Коэффициент избытка воздуха возрастает по мере движения продуктов сгорания по тракту котлоагрегата вследствие присосов. Для коэффициента избытка воздуха в топке α Т энтальпия продуктов сгорания строится в диапазоне температур 1500 – 2000 С. На выходе из топки коэффициент избытка увеличивается на величину присосов Δα Т . Следует построить кривую для значения α=α Т + Δα Т в диапазоне температур 700 – 1200 С. На выходе из конвективного пучка коэффициент избытка воздуха α=α Т + Δα Т + Δα К , где Δα К - присосы в конвективном пучке. Для этого значения αследует построить кривую энтальпии продуктов сгорания для температур 300 – 700 С. Коэффициент избытка воздуха уходящих газов равен αух =αТ + ΔαТ + ΔαК + ΔαЭгде Δα Э - присосы в экономайзере. Кривая энтальпии для этого значения коэффициента избытка воздуха строится в температурном диапазоне 100 – 400 С. Так выполняется построение It -диаграмма дымовых газов. Расход топливаВычисляется располагаемое тепло, отнесенное к единице топлива (килограмму или кубометру) Q РР = Q НР +Q ф Здесь Q Н Р - теплота сгорания топлива, Q ф = с Т Т t - физическое тепло топлива, равное произведению теплоемкости топлива и его температуры. Теплоемкости топлива приведены в таблице Теплоемкость топлив
Задается температура уходящих газов t ух . Если не предполагается устанавливать экономайзер, то температура дымовых газов на выходе из котла принимается на 30 – 40 С выше температуры кипения воды в барабане котла, При использовании экономайзера - согласно таблице, приведенной ниже. Для твердых топлив температура уходящих газов выбирается по приведенной влажности топлива W П . Она равна W П = 4,19 10⋅ 3 W Р Q Н Ргде Q Н Р - теплота сгорания топлива, кДж/кг (куб.м), W Р - влажность рабочей массы топлива. Температура уходящих газов
Затем вычисляются потери тепла с уходящими газами Q 2 = −I ух αух вз I о где I вз о - энтальпия минимально необходимого количества воздуха для полного сгорания топлива (α= 1) , отнесенная к температуре в котельной (температуру воздуха в котельной принять 30 С, при этой температуре (ct ) вз = 34,0 кДж/м3 ). На It –диаграмме продуктов сгорания для определения энтальпии уходящих газов используется кривая, соответствующая α ух . Далее вычисляется значение q 2 , равное Q 2 q 2 = P ⋅100% . Q PЗатем задаются величины химического недожога топлива q 3 и механического уноса топлива q 4 . Соответствующая информация приведена в таблице. Потери с механическим уносом и химическим недожогом
Потери тепла в окружающую среду q 5 находятся с помощью графика, представленного на рис.1. По оси абсцисс указана паропроизводительность котла. Кривая 1 относится к котлу с хвостовыми поверхностями, кривая 2 – собственно к котлу (без хвостовых поверхностей). Рис. 1 Потери тепла в окружающую среду Для твердых топлив дополнительно оценивается величина потерь тепла в результате удаления шлаков q 6 . Она равна (1− a ун )(сt )зл A p q 6 = P % Q P Здесь a ун - доля золы топлива в уносе ( принять 0,1 ), (ct )зл - энтальпия золы ( принять 550 кДж/кг золы ), Ap - зольность рабочей массы твердого топлива (в процентах). Теперь можно оценить к.п.д. котлоагрегата ηК =100% − q 2 − q 3 − q 4 − q 5 − q 6Расход топлива котлоагрегатом B равен D i ( ПП − i ПВ ) + G пр (i ′ − i ПВ ) B = Р 100ηК Q Р Здесь D и G пр - соответственно паропроизводительность котлоагрегата и расход продувки, i ПП - энтальпия перегретого пара (если пароперегреватель отсутствует, принимается i ПП = i ′′), i ′′и i ′- соответственно энтальпия сухого насыщенного пара и кипящей жидкости (параметры на линии насыщения при давлении в барабане котла), i ПВ - энтальпия питательной воды на входе в котел. Температуру питательной воды принять равной 100 С. Адиабатная температура сгоранияТепло, вносимое в топку Q Т , равно Р ⎛100%− q 3 − q 4 ⎞⎟⎟+αТ вз I оQ Т = Q Р ⎜⎜ 100%− q 4 ⎠ ⎝ где I вз о - энтальпия теоретически необходимого количества воздуха при температуре воздуха в котельной. С помощью It –диаграммы продуктов сгорания по кривой, соответствующей α Т , определяется адиабатная температура сгорания так, как это показано на рисунке. Температура дымовых газов на выходе из топкиПо чертежам котлоагрегата вычисляется суммарная поверхность стен топки F ст . Затем вычисляются площади поверхностей топки, занятых экранами Экранов в топке может быть несколько. Например, фронтальный, задний, боковые и др. Площадь поверхности i -ого настенного экрана обозначим. F i Э . Она равна F i Э = [(т −1)s + d H ] Э Здесь m - количество экранных труб на экранируемой поверхности, s - шаг экранных труб, d и H Э - соответственно наружный диаметр и высота экранной трубы. Последняя величина определяется по чертежам котла. Средняя тепловая эффективность экранов равна ∑ψi F iЭ ψср = i . F ст где ψ i = x i ξ i , x i - угловой коэффициент i -ого экрана (определяется по номограмме), ξ i - коэффициент загрязнения этой экранной поверхности. Угловые коэффициенты однорядных гладкотрубных экранов определяются по графику, приведенному на рис.2 .На этом графике кривые 1 – 4 учитывают излучение обмуровки топки, кривая 5 - нет. При сжигании газа принятьξ= 0,65 , при сжигании мазута ξ= 0,55 и при сжигании углей в слое ξ= 0,60. По составу дымовых газов находятся объемные доли трехатомных газов V RO V H O 2 r RO 2 =r H O 2 = ГГ Здесь V RO 2 и V H O 2 - объемы трехатомных газов в продуктах сгорания, V Г - объем дымовых газов при значении коэффициента избытка воздуха на выходе из топки α=α Т + Δα Т . Далее вычисляется эффективная степень черноты факела a ф . При сжигании газообразных и жидких топлив указанная степень черноты вычисляется по формуле aф = maсв + (1− m a) гПри сжигании газа m = 0,1 , при сжигании мазута m = 0,55. В приведенной формуле a св - степень черноты светящегося пламени, a г - степень черноты несветящегося пламени. Они определяются по формулам a св = 1− exp( (− k r г г + k c )ps ) a г = 1− exp(−k r ps г г ) Здесь k г - коэффициент ослабления лучей трехатомными газами, k c - коэффициент ослабления лучей сажистыми частицами, r г = r RO 2 + r H O 2 - объемная доля трехатомных газов в продуктах сгорания, p - давление в топке ( p = 1 атм), s - эффективная толщина излучающего слоя в топке. Она равна V Т s = 3.6F ст Здесь V T - объем топки в кубических метрах. Рис.2 Угловые коэффициенты экранов кривая 1 – e > 1,4d кривая 2 - e = 0,8d кривая 3 - e = 0,5d кривая 4 - e = 0 кривая 5 – e > 0,5d Коэффициенты ослабления лучей трехатомными газами и сажистыми частицами можно вычислить с помощью программы. Коэффициент k c , в частности, зависит от соотношения содержания в топливе углерода и водорода, определяемого как C Р m Р ∑ m n = 0,12 C H H nЛевое отношение действительно для мазутов, правое – для газов. При сжигании твердых топлив эффективная степень черноты факела рассчитывается по формуле aф =1−exp(−kps)где k - коэффициент ослабления лучей топочной средой, p - давление в топке ( p = 1 атм), s - эффективная толщина излучающего слоя в топке (см. выше). Коэффициент ослабления лучей равен k = k r г г + k зл μзл +κ1κ2 Здесь k зл - коэффициент ослабления лучей зольными частицами, μ зл - безразмерная концентрация золы в дымовых газах, κ 1 и κ 2 - коэффициенты, зависящие от рода топлива и способа его сжигания. Так при сжигании бурых и каменных углей κ 1 = 0,5 . При слоевом сжигании твердого топливаκ 2 = 0,03. Коэффициент k г определяется так же, как при сжигании газа или мазута (см. выше). Для вычисления k зл с помощью программы нужно, в частности, знать содержание золы в топливе AP , средний диаметр частиц золы и величину ее уноса. При слоевом сжигании этот диаметр равен 20 мкм, а унос составляет 10%. Объем воздуха, подаваемого в топку, принять равным V д (см. выше). После определения a ф следует вычислить степень черноты топки a Т по формуле, приведенной ниже a ф a Т = 1− (1− a ф )(1−ψср )(1−ρ)где ρ= R - отношение площадей зеркала горения слоя топлива R к полной поверхности стен топки F ст Fст . Расчетный расход топлива равен q 4 B P = B (1− ) 100 Далее вычисляется критерий Больцмана ϕB P (Vc )срBo = 3 σψср F ст T а –8 2 4 q 5 Здесь σ= 5,73 10 Вт/м К , ϕ= (1− ) - коэффициент сохранения тепла, (Vc )ср - средняя ηK + q 5 теплоемкость продуктов сгорания в интервале температур от t ′′до t a . Последняя величина оценивается как Q T − I ′′ t a − t ′′(Vc )ср = где t a и t ′′ - соответственно адиабатная температура и температура дымовых газов на выходе из топки (в градусах Цельсия), I ′′ - энтальпия продуктов сгорания на выходе из топки. Последняя величина является целью расчета настоящего раздела. Поэтому необходимо предварительно задать величину t ′′ ( диапазон ожидаемых значений этой температуры 800 – 1100 С) и далее действовать согласно рисунку Кривая 1 на рисунке соответствует α Т , кривая 2 – значению α=α Т + Δα Т . Теперь вычисляется отношение температур T ′′и T a (температуры в градусах Кельвина) по следующей формуле T ′′ Bo 0,6 = 0,6 0,6T a M a T + Bo Здесь коэффициент M при сжигании газа и мазута равен M = 0,54 − 0,2 x T а при сжигании твердых топлив M = 0,59 − 0,5x TВ обеих формулах x T есть отношение высоты расположения горелки от пода топки к высоте топки. При сжигании твердого топлива в тонком слое x T = 0, а при сжигании в толстом слое x T = 0,14. С помощью It –диаграммы дымовых газов (кривая α=α Т + Δα Т ) по величине T ′′или t ′′ = T ′′ - 273 определяется значение I ′′ (см. рисунок выше). Если значение t ′′будет сильно отличаться от величины этой температуры, принятой при расчете средней теплоемкости дымовых газов (Vc ) ср (см. выше в этом разделе), то выполняют коррекцию принятого значения, и расчет повторяют. Количество тепла, воспринятое в топке за счет излучения, отнесенное к единице топлива, равно Q л =ϕ(Q T − I ′′) Расчет пучка кипятильных трубЗадается температура продуктов сгорания на выходе из пучка кипятильных труб t K ′′ . С помощью этой температуры и кривой на It –диаграмме, соответствующей значению α=α Т + Δα Т + Δα К , находится энтальпия продуктов сгорания за конвективным пучком I K ′′ . И далее определяется тепловосприятие пучка кипятильных труб Q K ′ по формуле теплового баланса Q K ′ = I K ′ − I K ′′ + ΔαK I взоЗдесь I K ′ - энтальпия продуктов сгорания на выходе из топки ( I K ′ = I ′′). Исходя из конструктивных размеров котлоагрегата, необходимо определить минимальное живое сечение для прохода дымовых газов в пучке кипятильных труб. С помощью программы определяется коэффициент теплоотдачи α K со стороны продуктов сгорания для коридорного пучка труб при поперечном обтекании его газами с температурой t ср = 0,5(t K ′ + t K ′′ ), где t K ′ температура дымовых газов на выходе из топки. Коэффициент теплопередачи от дымовых газов к пароводяной смеси в кипятильных трубах равен k K =ψ(αK +αл ) Здесь ψ- коэффициент тепловой эффективности (ψ= 0,60 – 0,65). Величину коэффициента теплоотдачи излучением α л так же следует вычислить с помощью программы. Температура стенок труб пучка принимается равной t w = t нс + Δt , где t нс - температура нагреваемой среды. Для котельного пучка она равна температуре кипения воды в барабане котла. При сжигании твердого или жидкого топлива Δ =t 60C , при сжигании газа Δ =t 25C . Далее с помощью уравнения теплопередачи оценивается тепловосприятие пучка кипятильных труб Q K ′′k F K K Δt лог Q K ′′ =B P Здесь F K - площадь поверхности теплообмена кипятильного пучка (принимается согласно конструктивным размерам котлоагрегата), Δt лог - температурный напор. Он вычисляется по следующей формуле. Δt ′ − Δt ′′ Δt лог = Δt ′ln Δt ′′ где Δt ′ = t K ′ − t S , Δt ′′ = t K ′′ − t S , t S - температура кипения воды в барабане котла. Если величины Q K ′ и Q K ′′ сильно различаются (более, чем на 5%), то необходимо скорректировать энтальпию продуктов сгорания I K ′′ . Так, если Q K ′ > Q K ′′ , то I K ′′ или t K ′′ следует увеличить, и наоборот. Тепло, расходуемое котлом на генерацию пара, определяется по формуле D i ( ПП − i В ) + G пр (i ′ − i ПВ ) Q KA = B P Здесь i ПB - энтальпия питательной воды на входе в котел. Если пароперегреватель отсутствует, то на выходе из котла пар – насыщенный и его энтальпия принимается равной i ПП = i ′′. Принимая далее QK ≅ Q K ′ ≅ Q K ′′ , определяем тепло, расходуемое на подогрев питательной воды в экономайзере Q Э = Q КА − Q л −Q K Если величина Q Э близка к нулю, то экономайзер не нужен, в противном случае вычисляется энтальпия дымовых газов за экономайзером по формуле I Э ′′ = I Э ′ − Q Э + ΔαЭ вз I о где I Э ′ и I Э ′′ - энтальпия дымовых газов соответственно на входе и выходе из экономайзера, Δα Э - присосы в экономайзере. Очевидно, что I Э ′ = I K ′′ . Если I Э ′′ и энтальпия дымовых газов на выходе из котла I ух , принятая ранее (см. раздел «Расход топлива»), примерно одинаковы (с точностью до 5%), то приступают к расчету экономайзера. В противном случае корректируется значение температуры уходящих газовt ух , и расчет повторяется, начиная с пункта «Расход топлива». Расчет чугунного экономайзера ВТИ Энтальпия воды i В на выходе из экономайзера составляет B Q P Эi B = i ПВ + D + G пр По величине i В и давлению питательной воды, равному давлению в барабане котла, определяется температура воды на выходе из экономайзера t B . Температуры дымовых газов на входе t Э ′ и выходе из экономайзера t Э ′′ определяются с помощью It – диаграммы продуктов сгорания (см. рисунок ниже). На рисунке кривая 1 соответствуетα=αТ + ΔαТ + ΔαК , кривая 2 - αух =αТ + ΔαТ + ΔαК + ΔαЭ . Средняя температура продуктов сгорания в экономайзере равна t ср = 0,5(t Э ′ + t Э ′′) Далее определяется объем дымовых газов V Э , поступающий в экономайзер 273+ t 2730 ср V Э = 0,5(α α+ ух )V B вз P Среднюю скорость продуктов сгорания в экономайзере следует принять равной w Э = 7 – 8 м/сек. Тогда площадь живого сечения для прохода газов равна V Э S Э = w Э Выбирается трубы экономайзера (см. таблицу ниже), и вычисляется количество труб в одном ряду поперек потока газов. S Э z Э = s тр где s тр - живое сечение для прохода газов одной трубы. Если z Э не целое число, то выполняется округление до целого значения, и корректируются величины S Э и wЭ . V Э S Э = z s Э тр w Э = S Э Размеры трубы чугунного экономайзера
С помощью программы вычисляется коэффициент теплопередачи k Э экономайзера. Величина поверхности теплообмена экономайзера равна Q B Э PF Э = k Э Ψ Δt лог В этом выражении логарифмический напор вычисляется по формуле Δt 1 − Δt 2Δ t лог = Δt ln Δt 2где Δt 1 = t Э ′ − t B и Δt 2 = t Э ′′ − t ПВ . Коэффициент пересчета ψот противоточной схемы к более сложной определяется по номограмме, приведенной на рис.3. Эта номограмма применима для теплообменников, схема течения теплоносителей в которых организована так, как это указано на рис. 4. Предварительно вычисляются две разности температур: для дымовых газов υ′−υ′′ и для воды t ′′− t ′. Большая из этих разностей – это τ б , меньшая обозначена как τ м . Далее вычисляются величины P и R , равные соответственно τм τбP = R =υ′ − t ′ τм Если число ходов больше четырех принимается ψ = 1,0 . Количество рядов труб по ходу потока дымовых газов равно F Э z = f z Э ЭЛитература 1. Тепловой расчет котельных агрегатов (нормативный метод), «Энергия», Москва, 1973 г., 300 стр. Рис. 3 Коэффициент ψ для перекрестного тока. Рис. 4 Схемы течения теплоносителей |