Введение Получение дополнительной энергии за счет утилизации тепловых потерь двигателя внутреннего сгорания (ДВС) позволяет экономить топливо, расходуемое на работу судовой энергетической установки (СЭУ). Это соответственно ведет к уменьшению выбросов вредных веществ в атмосферу, способствует удовлетворению более жестких норм Международной морской организации (ИМО – International Maritime Organization, IMO), регламентирующих пределы этих выбросов. Именно поэтому важно обеспечить максимально возможную утилизацию тепловых потерь главного двигателя (ГД). В связи с существенным снижением температуры выхлопных газов современных ДВС и резким увеличением скорости низкотемпературной сернокислотной коррозии (до 1,5 мм/год при температуре стенки ниже 130 ºС при сжигании товарных сернистых топлив с водосодержанием Wr = 2 %) имеет место существенное снижение паропроизводительности утилизационных котлов (УК) и глубины утилизации теплоты газов. Уровень использования теплоты сжигания топлива может быть повышен путем глубокого охлаждения продуктов сгорания ниже точки росы, т. е. за счет максимального использования не только физической теплоты газов, но и скрытой теплоты содержащихся в них паров серной кислоты, конденсация которой и ограничивает температуру выхлопных газов при значениях температуры стенки tст ниже температуры точки росы паров серной кислоты (ТТР H2SO4). Именно поэтому снижение интенсивности низкотемпературной сернистой коррозии (НТК) при значениях температуры поверхности ниже ТТР H2SO4 является практически единственной возможностью снижения температуры уходящих из УК газов и повышения экологичности и экономичности СЭУ. Анализ последних исследований и публикаций В [1] представлены результаты экспериментальных исследований интенсивности коррозии и загрязнения поверхностей нагрева УК марки КУП-1300, установленного после газотурбинного двигателя ГТД М-25, при сжигании топлива ДС с содержанием серы Sr = 0,1 % при избытке воздуха a = 6,5 и 8,4 (при частичной мощности ГТД) при tст ниже ТТР H2SO4. Продолжительность коррозионных испытаний составляла 15…19 ч, скорость НТК находилась в области «кислотного пика» в пределах 0,24…0,32 г/(м2·ч), а скорость образования отложений – 0,27…0,74 г/(м2·ч), ТТР составляла 115 ºС. В [2] опубликованы данные по скорости НТК и роста отложений, полученные при проведении исследований на вспомогательном котле при a, соответствующих избытку воздуха в ДВС (2,3…2,5) при сжигании мазута с содержанием серы на уровне 1,6…1,7 %. Продолжительность испытаний составляла 6,5 ч, потеря массы образцов вследствие коррозии – 3,75 г/м2, а масса отложений коксосажистых частиц − 25 г/м2 в районе кислотного пика при tст = 110 ºС. Авторами [3] представлены результаты эксплуатационных испытаний котла КВГ-34К при сжигании товарного топлива М100В и водотопливной эмульсии (ВТЭ) на его основе с водосодержанием Wr = 30 % при a = 1,3…1,4. При сжигании ВТЭ наблюдалось снижение скорости отложений и изменение их структуры на поверхностях нагрева котла. Изменился характер коррозии, наблюдалось отсутствие язвенной и наличие только равномерной поверхностной коррозии. Но в [3] отсутствуют количественные данные по интенсивности развития этих процессов в зависимости от a, Wr и тем более tст, т. к. для этого необходимо проведение специальных исследований. По данным [4], при сжигании ВТЭ отложения становятся рыхлыми, легко удаляются водой или вообще отсутствуют вследствие уменьшения генерирования сажи и кокса, т. к. при горении капель ВТЭ наблюдаются микровзрывы капель, которые интенсифицируют процесс горения, что наблюдается и в камерах сгорания ДВС [5]. Но и в этих литературных источниках нет количественных данных по интенсивности НТК. В [6] представлены результаты исследований процесса коррозии, полученные на экспериментальной установке при сжигании ВТЭ на основе мазута М40 (Sr = 1,8 %, Wr = 30 %) при коэффициенте избытка воздуха a = 1,45 и при сжигании того же мазута при водосодержании Wr = 2,0 % и a, равном 1,05…1,15, т. е. при a, соответствующем режимам сжигания топлива в водогрейном котле. Продолжительность опытов составляла от 2 до 100 часов. В [6] представлены количественные зависимости скорости коррозии К в области кислотного пика при tст = 110 °С: сравнительный уровень коррозии при сжигании мазута с Wr = 2 % К = 0,63 г/(м2·ч) при a = 1,05, при сжигании ВМЭ К = 0,17 г/(м2·ч) при a = 1,45. В [7] даны количественные зависимости интенсивности загрязнения Кз в том же диапазоне a: при сжигании мазута с Wr = 2 % Кз = 30 г/(м2·ч), при сжигании ВМЭ Кз = 25 г/(м2·ч) при времени испытаний 8 часов. В опубликованных работах по эксплуатации ДВС на ВТЭ [5] не рассматриваются вопросы интенсивности коррозионных процессов в УК, определяющих надежность работы и тепловосприятие УК. Интенсивность процессов коррозии зависит в основном от содержания серы в топливе, коэффициента избытка воздуха, при котором идет горение, содержания воды в топливе (или ВТЭ), значений температуры поверхности нагрева. С учетом изложенного были поставлены следующие цели исследований: 1. Получение зависимостей удельных показателей интенсивности НТК ΔGк (по убыли металла) и загрязнений ΔGз (по приросту их массы) при влиянии: коэффициента избытка воздуха a (до 2,5…2,9), соответствующего сжиганию топлива в ДВС, содержания серы Sr (0,98…1,8 %) и водосодержания Wr (2…30 %). 2. Получение зависимостей скорости низкотемпературной коррозии К от tст, которые необходимы для проектирования поверхностей УК. 3. Оценка корректности переноса данных, полученных на экспериментальной установке, на рассмотрение развития этих процессов в реальных УК, установленных после ДВС. При планировании кратковременных экспериментальных исследований необходимо учесть механизм влияния a и Sr на интенсивность коррозионных процессов в рассматриваемом диапазоне их изменения. Интенсивность коррозионных процессов определяется влиянием a и Sr на образование SO3 и количество паров H2SO4, при конденсации которых начинается развитие НТК и осаждение твердых частиц. В пределах роста a от 1,01 до 1,1 наблюдается стремительный рост содержания SO3 (рис. 1, а), а значит, и рост количества паров H2SO4. Затем, при более высоких a, имеет место стабилизация содержания SO3. В пределах изменения содержания серы от 0,5 до 1,0 % имеет место значительный рост ТТР tр (рис. 1, б), значение которой также характеризует содержание паров H2SO4 в выхлопных газах. При Sr больше 1 % значения температуры tр стабилизируются. Обработка данных, представленных на рис. 1, указывает на то, что развитие коррозионных процессов в зависимости от a и Sr достоверно описывается степенными зависимостями, которые и были использованы при обработке данных исследований коррозионных процессов. Следовательно, вид зависимости и уровень развития интенсивности НТК и загрязнения низкотемпературных поверхностей нагрева (НТПН) УК возможно определить по данным на начальном участке исследуемых зависимостей от a и Sr. Поэтому следует считать корректным получение аппроксимационных уравнений DGк = f (a) и DGз = f (a) по данным исследований на экспериментальной установке при значениях a от 1,05 до 1,5, а затем оценить применимость этих уравнений к процессам в УК по значениям DGк и DGз при a = 2,9, полученным при проведении экспериментов на реальном УК на судне. Диапазон изменения Wr необходимо выбирать в пределах от 2 до 40 %, т. к. при увеличении Wr, по данным [6], имеет место существенное снижение интенсивности НТК и загрязнения [7] НТПН при tст ниже ТТР Н2SO4. Учитывая разные подходы к оценке водосодержания, необходимо отметить, что при приготовлении ВТЭ с водосодержанием 30 % смешивались вода и топливо в соотношении 300 г воды Gв и 700 г топлива GТ, что означает 30 % воды в 1 кг ВТЭ (обозначается Wr = 30 %). Этому значению Wr соответствует водотопливное отношение Gв/GТ = (300/700) 100 % = 42,86 %. а б Рис. 1. Влияние коэффициента избытка воздуха и содержания серы на механизм НТК: а – на концентрацию SO3; б – на ТТР (при a = 1,25): 1 – вид кривой по уравнению SO3 = 4,5487β0,2418 tр = 120,306Sr 0,143 R2 = 0,6566 R2 = 0,9755 На основании результатов кратковременных экспериментальных исследований были получены регрессионные уравнения динамики развития коррозии и загрязнения при влиянии a, Sr и Wr тоже в виде степенных зависимостей DGк = f (τ) и DGз = f (τ). Для построения искомых зависимостей DGк = f (a) и DGз = f (a) были выбраны результаты 8-часовых исследований, т. к. в течение этого времени воздействия потока газов наступает динамическое равновесие развития процессов. Как показывают представленные на рис. 2, а зависимости DGк = f (a), полученные на экспериментальной установке при одинаковом содержании серы (Sr = 1,5 %) при Wr = 2 % в топливе и Wr в эмульсии 15 и 30 %, с увеличением a наблюдается рост DGк при любых значениях Wr в эмульсии, хотя абсолютные значения DGк увеличиваются с ростом значения a. Если при сжигании стандартного топлива (Wr = 2 %) при увеличении a с 1,05 до 2,9 DGк возросло в 2,2 раза, то при сжигании ВТЭ с Wr = 30 % – в 1,9 раза, но на существенно более низком уровне интенсивности НТК. При сжигании ВТЭ с Wr = 30 % значение DGк в 2,8 раза ниже, чем при сжигании стандартного топлива при a = 2,9. Достоверность результатов, рассчитанных по уравнениям DGк = f (a), была оценена путем сравнения расчетного значения DGк при a = 2,9 с данными экспериментальных исследований при сжигании смеси топлив ДТ + ДЛ (Wr = 2 % и Sr = 1,5 %) a = 2,9 на судне во время ходовых испытаний на УК типа КУП 1000Р, установленного после ГД марки 8ДКРН 60/195-10. Расхождение составило 1,9 % (рис. 2, а). Следовательно, уравнения, полученные по данным результатов исследований НТК при a в диапазоне 1,05…1,5 (см. подрисуночные подписи), достоверно описывают ход процесса и значения скорости НТК и в УК после ДВС при больших значениях a. Это подтверждается одинаковыми значениями соотношений между величинами DGк при Wr, равных 2, 15 и 30 %, составляющих соответственно (3; 2; 1) при a = 1,05 и при a = 2,9, хотя абсолютные значения DGк при a = 2,9 естественно выше. Рис. 2. Зависимость удельной убыли металла ΔGк вследствие НТК при τ = 8 ч: а – от избытка воздуха α при разном водосодержании Wr ВТЭ при Sr = 1,5 %: Wr = 2 % – ΔGк = 13,9516 a 0,3841, R2 = 0,9999; Wr = 15 % – ΔGк = 8,9149 a 0,5134, R2 = 1; Wr = 30 % – ΔGк = 4,2105 a 0,568, R2 = 0,8852; б – от содержания серы Sr при разном водосодержании Wr ВТЭ при α = 2,9: Wr = 2 % – ΔGк = 9,8129 S 1,9308, R2 = 0,9989; Wr = 15 % – ΔGк = 7,4156 S 1,6821, R2 = 0,9895; Wr = 30 % – ΔGк = 5,2975 S 1,2173, R2 = 0,8956; в – от водосодержания Wr ВТЭ при разном содержании серы Sr при α = 2,9: Sr = 1,8 % – ΔGк = 33,7985е – 0,0341, R2 = 0,9975; Sr = 1,5 % – ΔGк = 22,783е – 0,0363, R2 = 0,9941; Sr = 0,98 % – ΔGк = 9,8871е – 0,0202, R2 = 0,9999 Следовательно, принятая методика определения зависимостей убыли метала DGк путем построения аппроксимационных уравнений, полученных при меньших значениях a, по которым представлено достаточное количество публикаций [1, 6], корректна и обеспечивает получение достоверных данных и при больших a. Поэтому считаем, что экспериментальные данные развития процессов коррозии и загрязнения в остальных режимах, полученные на экспериментальной установке при сжигании ВТЭ, при рассмотрении коррозионных процессов в УК также достоверны. Представленные на рис 2, б зависимости (удельной убыли массы металла вследствие коррозии) DGк = f (Sr) при одинаковом значении a = 2,9 при содержаниях Wr в топливе 2 % и эмульсии 15 и 30 % показали, что при сжигании топлива (Wr = 2 %) при увеличении Sr с 0,98 до 2 % DGк увеличилось в 3,9 раза, при Wr = 15 % – в 3,2 раза при тех же условиях. При сжигании ВТЭ с Wr = 30 % наблюдается незначительный рост DGк: при увеличении Sr с 0,98 до 2 % DGк увеличилось в 2,3 раза. При этом также необходимо отметить, что значение убыли массы металла DGк при сжигании ВМЭ (Wr = 30 %) в 1,7 раза ниже, чем при сжигании топлива (Wr = 2 %) при Sr = 0,98 % и в 3 раза ниже, чем при Sr = 2 %. На основании данных DGк = f (a) и DGк = f (Sr) (рис. 2, а, б) построены зависимости DGк = f (Wr) при Sr, равных 0,98, 1,5 и 1,8 % при a = 2,9 (рис. 2, в). Полученные зависимости показали, что с увеличением Wr в эмульсии DGк снижается. Чем выше содержание серы в топливе, тем круче снижается кривая зависимости DGк = f (Wr). При Sr = 0,98 % при увеличении содержания воды в топливе и ВТЭ от 2 до 30 % DGк снижается в 1,7 раза, при Sr = 1,8 % в тех же пределах изменения Wr DGк снижается в 2,6 раза. То, что значения DGк при разном Sr существенно сближаются при Wr = 30 % (рис. 2, в) объясняется тем, что при этом Wr в эмульсии создаются условия для пассивации поверхности металла и резкого снижения скорости НТК [8]. Результаты исследований и литературные данные [7, 9] показывают, что переход на сжигание ВТЭ с Wr = 30 % приводит к уменьшению количества отложений твердых золовых и сажистых частиц из потока газов в 5 раз по сравнению с Wr = 2 % вследствие улучшения качества горения топлива из-за микровзрывов капель эмульсии ВТЭ. По данным [7] на рис. 3 представлены графики изменения относительного уменьшения интенсивности отложений δотл, сульфатов δс, остаточной кислоты δк.ост и их суммы, т. е. загрязнений δз, в зависимости от Wr, которые позволяют оценить изменение их количества при сжигании ВТЭ с Wr от 2 до 30 % путем следующих построений. По данным исследований интенсивности загрязнения ΔGз при Wr = 2 %, полученным на экспериментальной установке при α в диапазоне 1,05...1,5, и при α = 2,9 в УК во время ходовых испытаний на судне, построена зависимость DGз(Wr=2 %)ДТ+ДЛ = f (α) при Wr = 2 % (рис. 4). Уменьшение значений ординат зависимости DGз(Wr=2 %)ДТ+ДЛ = f (α) в δз раз (по данным рис. 3 при Wr = 30 %) при каждом α (2,9; 2,8 ...) позволяет получить зависимость DGз(Wr=30 %)ДТ+ДЛ = f1 (α) при Wr = 30 %. С целью оценки влияния большего количества золы в топливе на загрязнения НТПН по экспериментальным данным построена зависимость DGз(Wr=2 %)М40 = φ (α) при сжигании М40. По принятой методике построена зависимость загрязнения при сжигании ВТЭ на основе М40 DGз(Wr=30 %)М40 = φ1 (α). Рис. 3. Зависимость удельного увеличения массы отложений от водосодержания ВТЭ Рис. 4. Зависимость удельного увеличения массы загрязнения ΔGз от α при Sr = 1,5 % при τ = 8 ч Сравнение значений интенсивности загрязнения для каждой марки топлива показывает, что при сжигании стандартных топлив с Wr = 2 % увеличение массы загрязнения в 1,2 раза выше, чем при сжигании ВТЭ с водосодержанием Wr = 30 % на основе этих же топлив. При сжигании М40 с большим содержанием золы и ВТЭ на основе этого же топлива DGз в 3 раза выше, чем при сжигании (ДТ + ДЛ). По данным [10], по мере нарастания отложений на НТПН с tст ниже ТТР H2SO4, приблизительно через 85 часов толщина слоя достигает 4 мм и стабилизируется, т. к. на поверхности слоя отложений температура становится равной и выше ТТР H2SO4. В соответствии с представленным механизмом образования отложений большое значение имеет размер частиц: частицы с радиусом менее 1,2·10-5 м при своем движении следуют линиям тока газов (и поэтому практически не осаждаются), а более крупные частицы как бы «пробивают пограничный слой» и осаждаются. При сжигании ВТЭ, благодаря микровзрывам капель ВТЭ, размеры частиц в основном меньшего размера и находятся на уровне 1,2·10-6 м, поэтому загрязнение НТПН резко снижается. Чтобы оценить совместное влияние 3-х рассмотренных факторов (a, Sr, Wr) на скорость коррозии в виде DGк, определить значение констант и весовых коэффициентов в уравнении многофакторной регрессии, определиться с факторами существенного влияния на DGк для статистической обработки экспериментальных данных, была использована система STATGRAPHICS Plus for Windows, которая обеспечивает доступ к полному набору статистических методов и предоставляет возможность провести расширенный регрессионный анализ. Модуль системы STATGRAPHICS Plus for Windows, посвященный планированию экспериментов, дает возможность учесть взаимодействие анализируемых факторов на величину DGк. Статистическая обработка данных позволила получить уравнение, учитывающее влияние на величину ΔGк коэффициента избытка воздуха a, содержания серы Sr в топливе и воды Wr в ВТЭ: По полученной диаграмме Парето (рис. 5) видно, что содержание серы Sr (фактор С), содержание Wr в эмульсии (фактор А) и a (фактор В) имеют статистически значимые эффекты. На это указывает то, что соответствующие столбцы пересекают вертикальную линию, которая соответствует 95 % уровню вероятности для определения значимости влияния факторов. Таким образом, подводя итог анализу результатов экспериментов и их статистической обработки, можно утверждать, что на величину DGк влияет прежде всего содержание Sr в топливе, а затем воды в ВТЭ Wr и избытка воздуха a. Рис. 5. Карта Парето Достоверность полученных зависимостей DGк и DGз при τ = 8 часов можно оценить только путем их сравнения с результатами длительных опытов. Исследования длительностью 100 часов, в течение которых стабилизируются все процессы, предоставляют возможность достоверно провести эти сравнения и оценить точность прогноза развития коррозионных процессов на длительное время при сжигании ВТЭ, что нас интересует прежде всего. Сравнительная оценка интенсивности НТК проведена по зависимостям скорости коррозии от температуры стенки. Результаты экспериментальных и расчетных исследований скорости коррозии за 100 часов в виде зависимостей K = f (tст) представлены на рис. 6 с указанием допустимого уровня НТК (0,2 мм/год) в течение этого периода времени, в который укладывается реальные скорости развития НТК при сжигании только ВМЭ с Wr = 30 %. Полученные зависимости показывают, что при сжигании топлива с Wr = 2 % увеличение a с 1,05 до 2,9 приводит к увеличению скорости коррозии в 2,2 раза (что соответствует данным рис. 2, а). При сжигании ВТЭ с Wr = 30 % при увеличении a с 1,45 до 2,9 скорость коррозии возрастает в 1,4 раза. Сравнение скоростей НТК при a = 2,9 показало, что при сжигании ВТЭ с Wr = 30 % при tст = 110 °С (область «кислотного пика») скорость НТК в 6,8 раза ниже, чем при сжигании стандартного топлива с Wr = 2 %. Рис. 6. Зависимость K = f (tст) при сжигании различных топлив: 1 – М100 (Sr = 2,0 %; Wr = 2,0 %) a = 1,05, τ = 100 ч [6]; 2 – ДТ + ДС (Sr = 1,5 %; Wr = 2,0 %) a = 2,9, τ = 100 ч; 3 – ВМЭ М40 (Sr = 1,8 %; Wr = 30,0 %) a = 1,45, τ = 88 ч [6]; 4 – ВМЭ ДТ + ДС (Sr = 1,5 %; Wr = 30,0 %) a = 2,9, τ = 88 ч На рис. 7 представлены результаты расчетных исследований скорости НТК, прогнозируемые на t = 1 000 ч по полученным регрессионным уравнениям при обеспечении постоянных значений избытка воздуха a и значений температуры стенки (которые необходимо использовать при расчетах поверхностей нагрева УК). В течение этого периода времени в допустимом уровне интенсивности НТК находится вся конденсационная поверхность нагрева в пределах tст от 140 до 70 °С при Wr = 30 %, что позволяет снизить температуру уходящих газов из УК до 80…90 °С. При Wr = 2 % допустимый уровень интенсивности НТК может быть получен при a = 1,05, который может быть достигнут только в энергетических котлах. Рис. 7. Зависимость K = f(tст) при сжигании различных топлив (τ = 1 000 ч): 1 – М100 (Sr = 2,0 %; Wr = 2,0 %) a = 1,05 [6]; 2 – ДТ + ДС (Sr = 1,5 %; Wr = 2,0 %) a = 2,9; 3 – ВМЭ М40 (Sr = 1,8 %; Wr = 30,0 %) a = 1,45 [6]; 4 – ВМЭ ДТ + ДС (Sr = 1,5 %; Wr = 30,0 %) a = 2,9 Выводы Увеличение водосодержания Wr от 2 до 30 % при значениях температуры стенки в области «кислотного пика» (tст = 100…110 °С) приводит к снижению скорости НТК в 1,7…2,6 раза при содержании серы в топливе 0,98…1,8 % при a = 2,9 (благодаря процессу пассивации вследствие абсорбции оксидов азота). Интенсивность загрязнения при сжигании ВТЭ с Wr = 30 % на основе стандартного топлива в 1,2 раза ниже, чем при сжигании топлива (Wr = 2 %) этой же марки. Интенсивность загрязнения при сжигании М40 с большим содержанием золы и ВТЭ на основе этого же топлива в 3 раза выше, чем при сжигании (ДТ + ДЛ). Статистическая обработка результатов исследований показывает, что на интенсивность НТК преимущественное влияние оказывает содержание серы и воды в ВТЭ. При сжигании ВТЭ с Wr = 30 % скорость НТК в 6,8 раза ниже, чем при сжигании стандартного топлива с Wr = 2 %. За 100 часов в допустимый уровень НТК (0,2 мм/год) укладывается скорость НТК только при сжигании ВТЭ с Wr = 30 %. За 1 000 часов в допустимом уровне НТК находится вся конденсационная поверхность нагрева в пределах tст от 140 до 70 °С при Wr = 30 %. Принятая методика определения зависимостей убыли металла DGк путем построения аппроксимационных уравнений, полученных при меньших значениях a, корректна и обеспечивает получение достоверных данных и при больших a.