Распределение нагрузок на тэц с поперечными связями с учетом потокораспределения воды
На правах рукописи
РОМАШОВА ОЛЬГА ЮРЬЕВНА РАСПРЕДЕЛЕНИЕ НАГРУЗОК НА ТЭЦ С ПОПЕРЕЧНЫМИ СВЯЗЯМИ С УЧЕТОМ ПОТОКОРАСПРЕДЕЛЕНИЯ ВОДЫ 05.14.14 – Тепловые электрические станции, их энергетические системы и агрегаты
Автореферат диссертации на соискание ученой степени кандидата технических наук
Томск – 2007
Работа выполнена в Томском политехническом университете на кафедре атомных и тепловых электрических станций
Научный консультант: кандидат технических наук, доцент Беляев Л.А.
Официальные оппоненты:
доктор технических наук, профессор Андрюшин А. В.
доктор физико-математических наук, профессор Логинов В. С.
Ведущая организация: ОАО «КУЗБАССЭНЕРГО» (г. Кемерово)
Защита состоится 30 мая 2007 года в 12 часов на заседании диссертационного со вета К 212.269.04 при Томском политехническом университете по адресу:
634050, г. Томск, пр. Ленина, 30, корпус 4, ауд. 406.
С диссертацией можно ознакомиться в научно-технической библиотеке Томско го политехнического университета Автореферат разослан «29» апреля 2007 г.
Ученый секретарь диссертационного совета А.С. Заворин
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность темы. Одно из главных направлений повышения эффективности топливоиспользования на ТЭЦ – внутристанционная оптимизация режимов, кото рая дает значительную экономию топлива без дополнительных капитальных вло жений.
На современном этапе развития энергетики вопрос распределения нагрузок особенно актуален по причине возросшей конкуренции на рынке производства электрической и тепловой энергии, а также в связи с тем, что существенно изме нились объемы отпуска теплоты от теплоэлектроцентралей по сравнению с теми, которые были заложены при проектировании станций.
Задача распределения нагрузок на промышленно-отопительных ТЭЦ, имеющих в своем составе разнотипное турбинное оборудование и отпускающих тепло в виде пара разного потенциала и с горячей водой на несколько тепломаги стралей, относится к многомерным, поэтому очень сложна и окончательно не ре шена.
Кроме того, для ТЭЦ с поперечными магистралями проблема усугубляется еще и тем, что на характеристики турбин существенное влияние оказывают связи турбоустановок с общестанционными коллекторами пара и воды и гидравличе ские характеристики последних.
С учетом этого работа является весьма актуальной и своевременной.
Работа выполнялась в соответствие с основными направлениями научной деятельности Томского политехнического университета «Разработка методов и средств повышения надежности и эффективности эксплуатации энергетических объектов», в русле критических технологий Российской федерации «Технологии производства топлива и энергии из органического сырья».
Цель работы Целью работы является решение задачи распределения нагрузок и выбора состава работающего оборудования на ТЭЦ с поперечными связями.
Решаемые вопросы для достижения поставленной цели 1) Моделирование гидравлического тракта питательной воды и расчет потокораспределения.
2) Определение влияния потокораспределения питательной воды на энергетические характеристики групп турбоустановок.
3) Разработка математической модели многомерного распределения на грузок между теплофикационными турбинами.
4) Исследование эффективности отпуска тепла с сетевой и подпиточной водой.
Научная новизна полученных результатов состоит в следующем 1. Доказано влияние потокораспределения питательной воды на энергетические характеристики групп турбинного оборудования и распределение нагрузок.
2. Разработан алгоритм многомерного распределения нагрузок на ТЭЦ на основе метода динамического программирования с целью выбора оптимального состава работающих турбоустановок. Впервые используется аппарат двумерного динами ческого распределения тепловых нагрузок в прямой постановке с видоизменен ным критерием Беллмана.
3. Впервые предложена схема решения задачи потокораспределения для выбора состава работающих насосов.
4. Впервые аналитически получена формула для расчета оптимального распреде ления подогрева сетевой воды в двухступенчатой сетевой установке турбин типа ПТР- при использовании пара нижнего отбора на общестанционные нужды.
5. Обоснован способ распределения тепловой нагрузки между параллельно рабо тающими теплофикационными турбинами с двухступенчатым подогревом сете вой воды. Дано теоретическое обоснование полученного эффекта и выявлены границы эффективности применения данного способа.
Практическая значимость 1. Разработанные модели потокораспределения питательной воды могут быть ис пользованы при решении задач расчета нормативных удельных расходов топ лива (НУРТ) для ТЭЦ с поперечными связями.
2. Разработанные программные комплексы используется на Ново-Кемеровской ТЭЦ для выбора оптимального состава работающего оборудования и распре деления нагрузок.
3. Результаты исследований используются в учебном процессе в Томском поли техническом университете в дисциплинах «Режимы работы и эксплуатации ТЭС» и «Методы оптимизации и расчеты на ЭВМ технико-экономических за дач» для студентов специальности «Тепловые электрические станции».
4. Результаты исследований режимов турбоустановки ПТР-80-130/13 могут быть использованы для выбора оптимальной загрузки отопительных отборов в соот ветствии с изменением параметров теплосети.
Достоверность результатов обеспечивается :
- применением апробированных математических моделей и надежных методов вычислений;
- сравнением результатов, полученных различными методами расчета;
- хорошей сходимостью результатов вычислений с экспериментальными данны ми, полученными другими исследователями.
На защиту выносится 1). Научно-методические основы решения многомерной задачи распределения на грузок на ТЭЦ.
2). Результаты численных экспериментов по влиянию потокораспределения пита тельной воды на энергетические характеристики групп турбинного оборудования.
3). Аналитическое определение оптимального расчетного давления в нижнем ото пительном отборе при проектировании турбоустановок на совместный подогрев сетевой и подпиточной воды.
4). Способ повышения эффективности ступенчатого подогрева сетевой воды.
Апробация работы. Основные положения диссертационной работы лично представлялись и докладывались автором на IX,X,XI и XII Всероссийских науч но-технических конференциях «Энергетика: экология, надежность, безопасность» (Томск 2003, 2004, 2005, 2006 г.г.);
II и IV семинарах вузов Сибири и Дальнего Востока по теплофизике и теплоэнергетике (Томск, 2002;
Владивосток, 2005);
международной научно-технической конференции «Электроэнергия и будущее цивилизации», (Томск: 2004).
Публикации. Основные положения и результаты исследований, представ ленных в диссертации, опубликованы в 17-ти работах, среди которых 2 статьи в рецензируемых изданиях, а также материалы докладов на конференциях разного уровня.
Объем и структура диссертации. Диссертация состоит из введения, пяти глав, заключения, списка литературы (110 наименований) и приложений. Работа содержит 129 страниц с приложениями, 17 таблиц и 30 рисунков.
Личное участие автора. Автором выполнены полностью самостоятельно работы по разработке математических моделей решения задач оптимального рас пределения нагрузки между турбоагрегатами ТЭЦ сложной технологической структуры и потокораспределеня питательной воды в тракте подогревателей вы сокого давления, их алгоритмизации и созданию пакетов прикладных программ, разработаны методики проведения вычислений и обработки их результатов, про ведены все вычислительные эксперименты и анализ полученных данных. Поста новка задач исследований, обсуждение методики вычислительных экспериментов и осмысление полученных результатов выполнены с участием научного руково дителя к.т.н. Беляева Л.А.
ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ
Во введении обоснована актуальность темы исследования, сформулированы основные цели и задачи и намечены пути их решения.
В первой главе дается анализ современного состояния рассматриваемого вопроса.
Установлено, что промышленно-отопительные ТЭЦ с поперечными связями характеризуются сложной тепловой схемой и разнообразием основного и вспомо гательного оборудования. Оптимизация режимов работы таких ТЭЦ возможна с использованием математических моделей, адекватно отражающих работу элемен тов тепловой схемы, а также учитывающих тепловые и гидравлические процессы, характерные для станции с поперечными связями.
Наибольшую трудность при математическом описании станции представля ют энергетические характеристики теплофикационных турбин, экономичность которых на переменных режимах зависит от большого числа факторов.
Доказано, что задача распределения нагрузок на промышленно-отопительных ТЭЦ с поперечными магистралями имеет размерность больше трех, поэтому в общей постановке окончательно не решена.
Проанализированы способы распределения нагрузок между теплофикацион ными турбоагрегатами. Обзором установлено, что универсальным методом опти мизации, пригодным для любого вида ЭХ турбин и позволяющим учитывать лю бые ограничения, накладываемые на оптимизируемые переменные, является ме тод динамического программирования. Он успешно применяется для одномерно го распределения электрических нагрузок, однако использование его для ТЭЦ на талкивается на значительные трудности вычислительного плана.
Проанализированы эффективность ступенчатого подогрева сетевой воды и режимы загрузки противодавленческих турбин на действующих ТЭЦ в течение года.
На основе проведенного анализа поставлены задачи исследования.
Во второй главе показано влияние потокораспределения питательной воды на энергетические характеристики теплофикационных турбоустановок и распре деление нагрузок между ними.
На станциях с поперечными связями расход питательной воды Gпв через группы подогревателей высокого давления (ПВД) отдельных установок функцио нально не связан с расходом острого пара D0 на них.
В таблицае 1 приведены результаты испытаний для одной из ТЭЦ с началь ным давлением P0 =12,8 МПа, подтверждающие сказанное.
Таблица Сводная таблица значений расхода острого пара и питательной воды по отдельным тур боустановкам одной из ТЭЦ ОАО «Кузбассэнерго» P0 =12,8 МПа Тип Режим работы Ст. турбо-агрегата 1 2 3 № D0, D0, D0, D0, Gпв, Gпв, Gпв, Gпв, т/час т/час т/час т/час т/час т/час т/час т/час 1 2 3 4 5 6 7 8 9 7 ПТР-80-130/13 385 340 314 234 325 277 - 9 Р-50-130/7 - - - - - - - 10 Р-50-130/13 - - 204 290 156 393 187 11 ПТ-50-130/7 200 - 242 362 223 350 183 12 ПТ-50-130/7 220 320 230 327 - - 230 13 Р-50-130/18 290 390 290 367 - - - 14 ПТ-150-130/18 - - 638 603 605 612 296 Не смотря на разнообразие способов представления ЭХ, все они получены при равенстве расхода питательной воды расходу острого пара, а несоответствие значений Gпв и D0 учитывается введением соответствующих поправок.
Однако область использования таких поправок слишком мала ( ± 10 %) и не соот ветствует фактическому отклонению параметров в реальных режимах.
При значительном отклонении Gпв от D0 ЭХ отдельных турбин, получен ные при равенстве Gпв = G0, теряют физический смысл. Для ТЭЦ с поперечными связями они должны быть «разомкнутыми» по рабочему телу:
Qтуi = f ( N Эi, DПi, DТi, PПi, P, Gпвi,...) (1) Тi hпвi = f ( N Эi, DПi, DТi, Gпвi, PПi, PТi, Gовi,...), (2) и служить основой энергетической характеристики станции или отдельных групп оборудования QТУi = DТУi h0i Gпвi hпвi. (3) Для построения энергетической характеристики группы турбин необходимо совместное моделирование параллельно работающих агрегатов с выбором в каче стве основного показателя удельного расхода тепла на группу оборудования Qтуi q=. (4) N Эi Физический смысл влияния Gпв на показатели работы отдельных турбин от ражен в поправках к ЭХ. Добавим, что это влияние неоднозначно. При увеличе нии Gпв растет регенеративная выработка отборов на ПВД из-за увеличения рас хода греющего пара, однако энергетический эффект сдерживается ростом недог ревов в подогревателях. Так как давление пара в регенеративных отборах опреде ляется основным потоком, работающим в турбине, его изменение незначительно.
Учитывая практически прямо пропорциональную зависимость температуры на сыщения от давления и экспоненциальный характер зависимости недогрева от расхода воды, изменение температуры воды за ПВД отдельных турбин от расхода носит экстремальный характер. Расчеты подтверждают его наличие при снижении Gпв до (30-50) % от D0. Кроме того, в большинстве режимов на влияние потоко распределения накладывается влияние КПД проточной части вследствие измене ния объемных пропусков пара.
Влияние потокораспределения питательной воды на распределение нагру зок анализировалось на примере совместной работы двух однотипных турбин - Т 110-130 и ПТР-80-130/13 в характерных режимах работы. Критерием оптимально го распределения нагрузок во всех режимах принимался удельный расход тепла на группу q. Потокораспределение питательной воды задавалось относительным расходом воды на первую турбоустановку g пв1 = Gпв1 / Gпв при изменении Gпв1 от нуля до суммарного Gпв = Gпв1 + Gпв 2 при соблюдении условия равенства сум марных расходов питательной воды и острого пара на турбины Gпв = D0.
Оптимальному распределению нагрузок при заданном потокораспределении питательной воды g пв1 соответствует минимальное значение удельного расхода тепла на группу q опт. Изменение тепловой экономичности определялось относи опт тельно режима Gпв = G0 с показателем q (q опт q0 ) опт опт q = 100. (5) опт q На рис. 1 приведены выборочные результаты оптимального распределения электрической мощности при работе турбоустановок Т-110 в конденсационном режиме, которые позволяют установить, что до суммарных нагрузок, составляю щих примерно (60-70) % от суммы номинальных мощностей, неравномерная за грузка по питательной воде разных турбин оказывается более выгодной по срав нению с равномерной. При суммарной мощности выше 70 % оптимальное распре 1, N Э=120 мВт qОПТ, % N Э1ОПТ, % NЭ=120 мВт 90 1, 0, 60 50 0, 0, 0, gПВ 10 g ПВ 0 0 0,2 0,4 0,6 0,8 -0, 0 0,2 0,4 0,6 0,8 -0, Рис.1. Зависимость оптиальной загрузки Рис. 2. Перерасход тепла в оптимальных кон первой турбины от потокораспределе- денсационных режимах работы двух турбин ния питательной воды при совместной Т-110-130 в связи с изменением работе двух установок Т-110-130 потокораспределения в конденсационном режиме деление между турбинами – равномерное и не зависит от потокораспределения.
Расчет и анализ совместной работы турбоустановок в разных режимах по зволяет обобщить полученные результаты - преимущественно загружать по элек трической и тепловой мощности следует ту турбину, которая в соответствии с гидравлическим сопротивлением питательного тракта имеет большее значение Gпв. Объясняется это тем, что в большинстве случаев рост регенеративной выра ботки отборов на ПВД приводит к снижению относительных приростов тепла на установку.
Кроме того, по результатам исследований можно сделать вывод, что эконо мичность совместной работы турбин меняется при изменении g пв1 (рис. 2). В це лом отметим, что оптимальная загрузка агрегатов при несоответствии расходов Gпв и D0 при незначительной экономии в отдельных режимах, в большинстве случаев приводит к увеличению расхода тепла в среднем на (0,5-1,5) % с ростом до (2,5-3) % в отдельных режимах по сравнению с оптимальными режимами при Gпв = G0. Неучет потокораспределения снижает эффективность оптимизации в среднем на (0,3-1,0) %. При сочетании неблагоприятных факторов рост тепла в свежем паре может достигать (3-5) %.
В третьей главе дано математическое описание гидравлической модели питательного тракта ТЭЦ с поперечными связями, отмечены методические осо бенности потокораспределения питательной воды и представлен алгоритм расчета гидравлического контура параллельно работающих питательных насосов.
Математическая модель расчета гидравлической схемы ПТ может быть включена в расчет ЭХ турбин и котлов при их компьютерном моделировании ли бо использоваться при расчете поправок к ЭХ.
Принципиальная схема питательного тракта ТЭС с поперечными связями в общем виде приведена на рис. 3. Математическая модель распределения пита тельной воды по группам ПВД отдельных турбин строится на базе теории гидрав лических цепей (ГЦ). Топология схемы задается количеством ветвей n, узлов m, контуров c и матрицами: A - размера ( n n ) соединений ветвей в узлах;
и B ( c n ) – контуров. Движение рабочего тела в гидравлической цепи в установив шемся режиме происходит в соответствии с двумя законами Кирхгофа. Парамет рами, характеризующими режим, являются давления в узлах pi и расходы (пото ки ) gi на участках системы. Математической моделью питательного тракта ТЭЦ являются системы смешанных (линейных и нелинейных) уравнений, которые в матричной форме имеют следующий вид AG = Q;
(6) B Y = 0 ;
(7) Y + H = S Gm G (8) относительно неизвестных векторов G = ( g1, g 2,..., g n ) и Y = ( y1, y2,..., yn ), где G – вектор размера n расходов на участках цепи;
Q – вектор размера ( m 1 ) притоков (уравнительный коллектор деаэраторов) и оттоков (котлоагрегаты) сре ды в узлах;
Y - вектор размера n потерь давления на ветвях;
H - вектор размера n напоров питательных насосов;
S - диагональная матрица ( n n ) сопротивлений ветвей;
Gm - диагональная матрица ( n n ) модулей расходов рабочего тела. Для решения системы уравнений (6)-(8) используются методы контурных расходов и их «увязочные» варианты. Исходными данными для расчета являются давление в одном из узлов Pm, компоненты векторов Q и H, матрица сопротивлений S.
Отличие гидравлической модели ПТ от аналогичных, применяемых для расчета теплоснабжающих систем, в наличии большого числа (до (7-15)) парал лельно работающих мощных питательных насосов, в связи с чем гидравлическое описание питательного тракта дополняется напорными характеристиками каждо го агрегата H = H (G ) (9) N = N (G ).
Моделирование ГЦ с учетом уравнения (9) выполняется на основе последо вательного расчета контуров ПВД и ПЭН итерационным путем, в процессе кото рого выходные параметры одного контура являются входными для расчета друго го.
Этапы моделирования представлены на рисунках 4-6.
1). Суммарное сопротивление напорного коллектора принимается равным нулю, определяется состав параллельно работающих насосов и их напор H пар.
2). Составляется матричная модель «упрощенного» контура ПВД (напорный коллектор сосредоточен в одном узле) и для него решается задача потокораспре деления. Вектор притоков (оттоков) среды Q в узлах этой схемы формируется с учетом заданных значений расхода воды на котлоагрегаты DКАj (оттоки) и сум GПВД GПВД = GПВ g СН, где марного расхода через группы ПВД j GКА-4 GКА-3 G GКА-7 GКА-8 GКА- GКА-4 GКА-3 GКА- GКА- GКА-6 GКА-8 GКА- GКА-5 КА- 68 5 9 4 10 3 2 ТГ-2 ТГ- ТГ-3 ТГ-4 ТГ-1 ТГ-3 ТГ- ТГ- 1 2 5 3 ПН-1 ПН- PКОЛ=PД+HПАР GСН GПВ GПВД Рис.3. Принципиальная схема питатель- Рис.4. Расчетная гидравлическая схема ного тракта ТЭЦ с поперечными связями упрощенного контура ПВД gПВД1 GКА-6 GКА-7 GКА-8 GКА- gПВД2 GКА-4 GКА-3 GКА- gПВД3 gПВД gСН P5 15 8 13 4 3 5 P1 P2 P3 P4 P6 P7 P 7 6 ТГ-2 ТГ-3 ТГ- gСН1 ТГ- ПН- 2 1 2 3 2 1 5 6 12 10 9 7 6 P 9 15 10 11 12 PД GПН- GПН- GПВ Рис.5. Расчетная гидравлическая схема Рис.6. Расчетная гидравлическая схема контура ПЭН контура ПВД g СН - суммарный расход питательной воды на собственные нужды, обес j печивающий непрерывную продувку котлов и впрыск в РОУ. Таким образом оп ределяются расходы питательной воды по всем ветвям контура, включая ветви с ПВД g ПВД.
j 3). Составляется модель группы параллельно работающих насосов. Она описывает суммарную характеристику и решает задачу приведения характеристик отдельных агрегатов группы к некоторому (произвольно выбранному) узлу на порного коллектора. Оттоки воды из коллектора определяются с учетом найден ных выше значений g ПВД и известных g CН. На этом этапе определяются значе j j ния давления во всех узлах коллектора, а также подачи и напоры всех работаю щих насосов.
4). Выполняется расчет потокораспределения для контура ПВД, включающе го напорный коллектор. За Pm принимается давление в одном из узлов коллекто ра. В результате расчета уточняются значения g ПВДi.
Расчетный цикл повторяется, начиная с п.3 с уточненными значениями рас хода питательной воды через ПВД.
Моделирование приведенной напорной характеристики заключается в ите рационном расчете последовательности простых контуров, на которые разбивает ся исходная схема, для каждого из которых помимо выполнения двух законов Кирхгофа однозначно решается задача соответствия подач и напоров насосов.
Простой контур включает две ветки с насосами и участок коллектора между ни ми. Два соседних простых контура имеют общую ветку с насосом, что позволяет принимать выходные параметры одного контура за входные – соседнего. Гидрав лический расчет простого контура выполняется на основе совместного решения уравнений системы (6)-(9) и предназначается для расчета его гидравлики при лю бом значении подачи одного из насосов g1. Моделирование внешнего контура увязывает общий заданный расход на насосную группу с суммой подач всех включенных в работу насосов.
На базе приведенной методики реализованы в виде компьютерной про граммы матричные модели питательного тракта действующих ТЭЦ с начальным давлением P0 =8,8 МПа и P0 =12,8 МПа, имеющих в своем составе 5-8 групп ПВД и до десяти параллельно работающих питательных насосов. Погрешность расчета расходов составляет не более (3-5) %.
В четвертой главе выполнена постановка задачи распределения нагрузок на ТЭЦ с поперечными связями, обосновано использование метода динамическо го программирования, разработан алгоритм последовательного распределения те пловых и электрических нагрузок на базе декомпозиции исходной схемы ТЭЦ с целью снижения размерности оптимизационной задачи.
В общем случае на промышленно-отопительных ТЭЦ распределению под лежат электрическая мощность станции N, а также производственные и отопи тельные нагрузки разных параметров. Размерность оптимизационной задачи рав на (Z+Y+1), где Z- число групп турбин с одинаковым давлением пара в производ ственном отборе, Y- количество присоединенных к ТЭЦ тепломагистралей, объе диненных параметрами температурного графика и пропускной способностью гидравлической системы.
Распределение нагрузок сводится к минимизации функции Z Y S Z J Y K F = Qs ( N s, DПs, DТs ) + Q j ( N j, DПj ) + Qk ( N k, DТk y ) + z y z z =1 y =1 s =1 z =1 j =1 y =1 k = Z R Z Y G L + Qr ( DПr z ) + Qg ( DПg z, DТg y ) + Ql ( Nl ) z =1 r =1 z =1 y =1 g =1 l = (10) при следующих ограничениях I + J + K + R +G N i = N, N imin N i N imax ;
(11) i = m D Пi Z + DРОУ Z +1 = Dпр Z + DПБj + DРОУ Z 1 + Dс.н.Z ;
(12) i = DПi Z DПi Z DПi Z, z = 1 Z ;
min max (13) S + K +G Y DТ i + DРОУ = DТ y 1, + Dс.н. ;
(14) i =1 y = H DТi y = DТ y, 0 DТi y DТi y, y = 1 Y.
max (15) i = min max DТi DТi DТi (16) DОБj + DПБj = D БУj ;
(17) l = D БУj ;
Y (18) DТ j = = DТ Y + Dс.н.1,2, …, (19) DТ где n – общее число турбогрегатов станции, N – электрическая мощность отдель ных турбин;
DП – тепловая нагрузка производственного отбора;
DТ – тепловая нагрузка отопительного отбора;
S – число конденсационных турбин с двумя ре гулируемыми (производственным и отопительным) отборами;
J - число турбин с одним производственным отбором;
K - число турбин с одним отопительным от бором;
R, G – число противодавленческих турбин типа Р- и ПТР- соответственно;
L - число конденсационных турбин, m – количество турбин, имеющих производ ственный отбор с заданным давлением PZ ;
Н – количество турбин, отпускающих тепло из отопительных отборов на бойлерные установки тепломагистрали с пара метрами Y, l –число бойлерных установок, отнесенных к тепломагистрали Y;
Z DП Z - тепловые нагрузки турбоагрегатов с регулируемым отбором PП z ;
Dпр суммарная производственная нагрузка в паре потенциала PZ ;
DОБ, DПБ - значе ния расхода пара на основные и пиковые бойлера;
D БУ суммарный расход пара на бойлерную;
DТ Y - cуммарный отпуск пара из отопительных отборов турбин, питаемых бойлерные установки, отнесенные к тепломагистрали Y;
DРОУ Z +1, DРОУ Z 1 - расходы пара от РОУ c параметрами PZ +1 / PZ, подключенной к паровой магистрали более высокой ступени давления и с параметрами PZ / PZ 1 в ниже стоящий коллектор соответственно;
Dс.н.1,2 - расход пара на общестанционные те плообменники низкого потенциала (подогревателей сырой и добавочной воды, атмосферные деаэраторы и т.д.).
Универсальным методом оптимизации, пригодным для любого вида ЭХ турбин, является метод динамического программирования. Однако использование его для многомерных задач наталкивается на значительные трудности вычисли тельного плана.
Анализ состава и параметров различных ТЭЦ позволил предложить схему решения, в основе которой лежит известный в математическом моделировании прием – декомпозиция сложной структуры, позволяющий разбить исходную за дачу на составляющие меньшей сложности и требующие более простых алгорит мов. При этом для поиска глобального оптимума строится итерационная схема последовательного решения частных задач.
Такой подход не противоречит установленным на основе обзора литературы закономерностям последовательно распределения нагрузок на ТЭЦ с поперечны ми связями, однако имеет ряд отличий.
1). Впервые используется аппарат двумерного динамического программирования в прямой постановке для распределения тепловых нагрузок.
2). Достаточное условие сходимости алгоритма к глобальному оптимуму прове ряется в итерационном расчете.
3). За критерий оптимизации при работе по тепловому графику принимается ни максимум выработки на тепловом потреблении, а удельный расход тепла на груп Qi min. Эта величина, во-первых, учитывает, помимо пу оборудования q = Ni максимума выработки мощности на тепловом потреблении, экономичность кон денсационного потока, во-вторых, является обобщающим критерием при распре делении нагрузок по тепловому и электрическому графикам, и, в третьих, механизмом, включающим в рабо Z1 Y1 Z2 Y Ввод : N, D, D, D,D,...
П Т П Т ту итерационный цикл последова NК *= тельного распределения.
NТ 4). В отличии от предыдущих ра Распределение NК *=NК N NТ тепловых нагрузок бот предлагается ни комбинация Да различных методов оптимизации, N, Q, Теплофикационная NК =N -NТ группа 1 Тi i D,D Z1 Y Z1 Y а единый метод – динамического D,D Пi Тi П Т Q / N min Конденсационная программирования. группа i Тi Распределение N, Q, Использование аппарата ди Теплофикационная Тi i конденсационной группа D,D Z2 Y мощности NК намического программирования Пi Тi Z2 Y D,D П Т q min K дает возможность автоматически Q / N min i i Тi... определять оптимальный состав N N, Q, Эi Теплофикационная Тi i NК *=NК группаn Zn Yn турбинного цеха. В соответствии с D,D Zn Yn D,D Пi Тi Нет П Т Да разработанным алгоритмом (рис.7) Вывод: D, D, N, N,...
Q / N min Пi Тi Тi Эi i Тi исходная тепловая схема разбива Рис.7. Алгоритм распределения электрических и ется на пересекающиеся подмно тепловых нагрузок между турбоугрегатами ТЭЦ жества (группы оборудования).
в итерационном цикле Принципиально различают две группы, объединенные между собой конденсаци онными турбинами с отпуском тепла. Первая группа включает теплофикационные турбины различных типов (Т-, П-, ПТ-, Р-, ПТР-). Вторая группа - конденсацион ные турбины (К-, Т-, П-, ПТ-). Первая группа разбивается на подгруппы из агрега тов, объединенных отпуском производственного пара одного давления и (или) отопительной нагрузки на независимую тепломагистраль. Распределение тепло вых нагрузок осуществляется методом поступенчатого нагружения: последова тельно распределяются тепловые нагрузки D П Z 1, D П Z 2, …, при этом PП 1 PП 2..., избыток пара более низкого потенциала передается в группы более Z Z высокого потенциала вплоть до использования для покрытия заданных нагрузок РОУ свежего пара.
Распределение электрической мощности осуществляется в итерационном цикле последовательно в названных группах, в то время как тепловая нагрузка параллельно в подгруппах «теплофикационной» группы. Распределение тепловой нагрузки дает «привязанную» к ней электрическую мощность станции на тепло вом потреблении, значение которой определяет конденсационную мощность, рас пределение которой осуществляется между турбинами конденсационной группы.
Декомпозиция исходной структуры позволяет использовать аппарат динамиче ского программирования для малого числа турбоагрегатов при распределении не более двух видов нагрузки. Распределение электрической нагрузки производится на основе одномерного динамического программирования.
Для распределения тепловых нагрузок применяется метод двумерного дина мического программирования, который заключается в пошаговом построении функций согласно рекуррентным соотношениям Беллмана. Используя принцип оптимальности, распределение нагрузок сводят к решению n двумерных задач.
На первом шаге многостадийного процесса строятся функции f Q1 ( DП, DТ ) = Q1 ( DП1, DТ 1 ) ;
(20) f N1 ( DП, DТ ) = NТ 1 ( DП1, DТ 1 ) (21) DП1min DП1 DП1max ;
DТ 1min DТ 1 DТ 1max.
при На каждом последующем шаге в соответствии с рекуррентными соотноше ниями Беллмана определяются зависимости f QK ( DП, DТ ) = QК ( DПK, DТK ) + fQK 1 [( DП DП ( K 1) ),( DТ DТ ( K 1) )] ;
(22) f NK ( DП, DТ ) = NТK ( DПK, DТK ) + f N K 1 [( DП DПK ),( DТ DТK )] (23) DПK min DПK DПK max ;
DТK min DТK DТK max, с использованием которых решаются оптимизационные задачи Qk ( DПk, DТk ) + fQk 1 (( DП DПk ),( DТ DТk )) k ( DП, DТ ) = min FK = min NТk ( DПk, DТk ) + f Nk (( DП DПk ),( DТ DТk )) при изменении DПk min DПk DПk max ;
DТk min DТk DТk max, в результате чего определяются дискретные функции оптимальных значений DПk опт ( DП, DТ ) и DТk опт ( DП, DТ ).
Рис.8. Годовая выработ ка электроэнергии тур 1 2 боагрегатами ТЭЦ P0 = 12,8 МПа 6 1 – ПТР-80-130/13;
2 - Р 4 50-130/7;
3 - Р-50 5 130/13;
4 - ПТ-50-130/7;
5- ПТ-50-130/7;
6 - ПТ 50-130/18.
А) Оптимальное распределение Б) Фактическая загрузка % 1 2 3 4 5 6 7 8 9 10 11 Рис.9. Экономия тепла на ТЭЦ Рис.10. Оптимальная загрузка отопительных отбо P0 = 12,8 МПа по месяцам ров турбин ТЭЦ P0 = 8,8 МПа в зависимости от тем отчетного периода пературы наружного воздуха:
1 – Т-25-90, 2 – ПТ-25-90, 6 – Т-50-90, 7 – ПТ-25-90, 11,12,14 – Т-100-90.
Как видно из математического описания, особенность критерия оптимизации видоизменило рекуррентную формулу Беллмана по сравнению с ее классическим представлением.
На рисунках 8-10 представлены результаты выбора оптимального состава турбинного цеха двух ТЭЦ сложной технологической структуры с начальным давлением P0 = 12,8 МПа и P0 = 8,8 МПа, имеющих в своем составе разнообразное турбинное оборудование.
В пятой главе исследуется эффективность отпуска тепла с сетевой и подпи точной водой на ТЭЦ с поперечными связями Выбор оптимальных режимов работы турбоустановки типа ПТР-80- К покрытию отопительной нагрузки на промышленно-отопительных ТЭЦ широко привлекаются противодавленческие турбины. Значительный энергетиче ский эффект при их эксплуатации может быть получен лишь при эффективной за грузке отопительных отборов.
В работе исследовались режимы работы турбоустановки ПТР-80-130 (рис.
11), отработавший пар которой используется для общестанционных целей. В та ких условиях на распределение нагрузки между нижним и верхним отопительны ми отборами, помимо параметров теплосети, влияют также расходы пара в обще станционный коллектор и в пусковой конденсатор для нагрева подпиточной воды теплосети Gх.
Приведены результаты расчета схемы турбоустановки ПТР-80-130 при за данном расходе сетевой воды Gсв = 1200 т/ч и изменении расхода G х от нуля до максимальной пропускной способности конденсатора 1000 т/ч для трех режимов работы теплосети: 90/40;
100/50;
120/70. Давление в производственном отборе принято номинальным, расход - 100 т/час. Температура химически очищенной воды на входе в пусковой конденсатор – 30°С. На рис.12 представлен график сни жения доли подогрева сетевой воды в нижнем подогревателе при росте G х для разных параметров теплосети. Рост экономичности с увеличением расхода подпи точной воды в большинстве режимов двухступенчатого подогрева подтверждает ся увеличением выработки мощности на тепловом потреблении (рис.13) и сниже нием удельного расхода тепла на выработку электроэнергии q (рис.14). Однако при высоких температурах сетевой воды проявляется неоднозначный характер влияния G х на экономичность турбины. С одной стороны, увеличение расхода подпитки через пусковой конденсатор, при прочих равных условиях, увеличивает нагрузку ПК и пропуск пара через промежуточный отсек, что снижает давление Pт1 и увеличивает выработку мощности на тепловом потреблении. С другой сто роны, увеличение тепловой нагрузки конденсатора при увеличении G х повышает температурный напор заданной теплообменной поверхности, что ведет к росту Pт1 и обратному влиянию на экономичность. До определенных значений темпера туры теплосети наблюдается прирост мощности на базе теплового потребления, однако при высоких значениях P 2 рост температурного напора конденсатора на чинает компенсировать выигрыш от снижения Pт1. Экстремальная зависимость q от расхода G х в режиме температур теплосети 120/70 приведена на рис.14. Опти мальный расход подпитки через конденсатор в этом режиме составляет чуть бо лее половины максимальной пропускной способности конденсатора. При одно ступенчатом подогреве сетевой воды, когда нет эффекта от понижения давления отработавшего пара при росте G х, увеличение нагрузки ПК приводит к росту удельного расхода тепла q из-за увеличения температурного напора конденсато ра.
С целью более полного использования отопительных отборов турбины ПТР 80 в течение года оправданным является перевод турбины на подогрев подпиточ ной воды теплосети, при этом нагрев подпитки целесообразно осуществлять в пределах одной турбоустановки от параметров ХВО до температуры насыщения атмосферных деаэраторов подпитки. С целью повышения экономичности требу ется реконструкция схемы для организации двухступенчатого подогрева подпи точной воды в сетевых подогревателях с подключение ДПТС к верхнему отопи тельному отбору. Действующая схема, предусматривающая питание атмосферных деаэраторов из противодавления не обеспечивает надежности работы турбины P0 t0 DКОЛ PКОЛ ЧВД ЧНД Gx1 Gx,tx Рис.11.
ПК PП Принципиальная PТ2 PТ PКОЛ схема отпуска ПХВ тепла от ТУ ПТР-80-130/ DТ ВСП DТ2 НСП GСВ ДПТС 0,85 0, е 0,8 0, 100/50°С 0,75 0, 120/70°С 120/70°С 0,7 0, 0, 0, 100/50°С 0, 0, 0, 0, Gx, т/ч Gx, т/ч 0, 0, 0 100 200 300 400 500 600 700 800 900 0 100 200 300 400 500 600 700 800 900 Рис.12. Отношение доли подогрева в ниж- Рис.13. Доля выработки электроэнергии ней ступени к общему подогреву сетевой на тепловом потреблении паром сетевой воды в зависимости от расхода установки в зависимости от расхода под подпиточной воды питочной воды qе, % qе, % 2 0, 90/40°C Gxов, т/ч 900 1000 1100 1200 1300 1400 Gx, т/ч -0, 0 100 200 300 400 500 600 700 800 900 - - 100/50°C - -1, 120/70°C - Рис.14. Относительное изменение удельного Рис.15. Относительное изменение удель расхода тепла на выработку электроэнергии ного расхода тепла на выработку электро В режиме совместного подогрева сетевой и энергии при переводе турбоустановки подпиточной воды теплосети ПТР-80 на подогрев подпитки теплосети по причине разогрева и вибрации хвостовой части. Подключение деаэратора и верхнего сетевого подогревателя к одному отбору позволит сократить потери ра бочего тела до минимума, а увеличение площади поверхности нижней ступени сетевой установки приведет к снижению температурного напора ПК и пониже нию противодавления. Рассматривалась работа турбины при изменении расхода G х при нагреве в пределах турбоустановки ПТР-80 от температуры t x =30°С до температуры насыщения ДПТС 104°С. Наибольшее значение G х соответствует максимальной нагрузке отопительных отборов турбины ПТР-80-130. Расчетными исследованиями получено, что во всем диапазоне нагрузок основной нагрев под питочной воды осуществляется в нижней ступени. В диапазоне расходов G =(800-1400) т/ч тепловая экономичность – наивысшая. Оптимальное значение G составляет (1100-1200) т/ч, что соответствует (70-80) % отопительной на грузки турбины.
Факторы, определяющие эффективность ступенчатого подогрева сете вой воды при совместном отпуске тепла с сетевой и подпиточной водой Известно, что максимальная эффективность ступенчатого подогрева достига ется при равном подогреве по ступеням. Этот вывод нельзя использовать для тур боустановок, тепловой схемой которых предусмотрен отпуск пара из нижнего от бора на общестанционные нужды – для подогрева подпиточной воды теплосети или котлов (рис.10). При проектировании таких турбин расчетные параметры нижнего отопительного отбора следует определять с учетом максимальной выра ботки электроэнергии в течение года на базе совместного подогрева сетевой, до бавочной и подпиточной воды теплосети.
С использованием известной методики получена формула для расчета опти мальной температуры сетевой воды за нижним подогревателем турбоустановки GСВ GX tC 2 + t2 + (tос t1 ) + (t X 1 + t K 2 t1 ) GСВ + G X GСВ + G X opt tC1 =.
t + GСВ tос + GX t X opt tC1 = C В случае t1 = t2 ;
t2 + t K = 2 t1.
Способ повышения эффективности ступенчатого подогрева сетевой во ды на ТЭЦ при совместной работе турбоустановок Подогрев сетевой воды в нижнем и верхнем подогревателях двухступенча той сетевой установки в течение отопительного сезона значительно отличается от равномерного. При снижении температуры наружного воздуха тепловая нагрузка нижнего теплофикационного отбора и доля подогрева воды в НСП растет (до 70 80 % от номинальной тепловой нагрузки сетевой установки). При уменьшении расхода сетевой воды при прочих равных условиях неравномерность подогрева увеличивается. При совместной работе однотипных теплофикационных турбин по тепловому графику для обеспечения равномерного подогрева сетевой воды пред лагается сетевые установки отдельных турбин включать по воде последовательно.
Для этого каждая турбина переводится на одноступенчатый подогрев с регулиро ванием давления PТ 1 в нижнем отборе. Нижний и верхний сетевые подогреватели одной турбоустановки включаются по сетевой воде и пару параллельно.
Совместная работа двух турбоустановок Т-110-130 по тепловому графику по предложенной схеме дает повышение электрической мощности на (3-4) %, вы работки электроэнергии на тепловом потреблении - на (5-6) % при общем сниже нии удельного расхода тела на турбоустановку. Мощность теплофикационных отборов увеличивается на (4-9) %. При снижении нагрузки отопительных отборов эффект от равномерного распределения при последовательном соединении сете вых установок возрастает.
Выводы 1. На основе обобщенного анализа загрузки ПВД действующих промыш ленно-отопительных ТЭЦ с поперечными магистралями в различных режимах ра боты установлено, что отклонение расхода питательной воды от расхода острого пара для отдельных турбоустановок может достигать ± (50-100) % и тем самым влияет на энергетические характеристики групп турбинного оборудования и рас пределение нагрузки между ними.
2. Разработаны алгоритмы многомерного распределения нагрузок на ТЭЦ, на базе которых созданы программные комплексы с использованием метода ди намического программирования. Предложена схема двумерного распределения тепловых нагрузок в прямой постановке с видоизмененным критерием Беллмана.
3. Впервые сформулированы научно-методологические основы гидравличе ского расчета питательного тракта ТЭЦ с поперечными связями, которые дают возможность определять расходы воды в трактах ПВД и выбирать состав рабо тающих питательных насосов.
4. На основании расчетных исследований с использованием программных комплексов на примере нескольких ТЭЦ сложной структуры и разнотипными турбинными установками показано, что оптимальное распределение нагрузки да ет экономию топлива на уровне (4-12) %, в основном, за счет выбора состава ра ботающего оборудования.
5. Аналитически обоснованное распределение подогрева сетевой воды в двухступенчатой установке турбин, отпускающих тепло с сетевой и подпиточной водой может быть использовано на стадии технико-экономического обоснования расчетных параметров тепловых схем.
Основное содержание работы
отражено в следующих публикациях 1. Л.А.Беляев, О.Ю.Ромашова. Использование тепла промперегрева для увеличения отпуска тепла от турбоустановки // Известия Томского политехнического университета, 2002. – Т.
305. - № 2. - С.114-120.
2. Л.А.Беляев, О.Ю.Ромашова. Потокораспределение в тракте питательной воды ТЭЦ с попе речными связями // Известия Томского политехнического университета, 2002. – Т. 305. - № 2. - С.187-190.
3. Беляев Л.А., Ромашова О.Ю. Методические особенности расчета потокораспределения пи тательной воды на ТЭЦ с поперечными связями //Энергетика: Экология, надежность, безо пасность. Материалы докладов девятой Всероссийской научно-технической конференции Томск, ТПУ, 3-5.12.2003. - Томск: Изд. ТПУ, 2003. - с. 160-162.
4. Беляев Л.А., Ромашова О.Ю. Численное моделирование характеристики группы питатель ных насосов ТЭЦ //Энергетика: Экология, надежность, безопасность. Материалы докладов девятой Всероссийской научно-технической конференции - Томск, ТПУ, 3-5.12.2003. Томск: Изд. ТПУ, 2003. - с. 163-165.
5. Ромашова О.Ю., Габидуллин О.Р. Перераспределение поверхностей нагрева между нижней и верхней ступенями сетевой установки теплофикационных турбин //Энергетика: Экология, надежность, безопасность. Материалы докладов девятой Всероссийской научно технической конференции - Томск, ТПУ, 3-5.12.2003. - Томск: Изд. ТПУ, 2003. - с. 239-242.
6. Беляев Л.А., Ромашова О.Ю. Приложение динамического программирования к оптимизации режимов ТЭЦ на основе декомпозиции исходной структуры станции //Энергетика: Эколо гия, надежность, безопасность: Материалы докладов десятой Всероссийской научно технической конференции - Томск, 8-10 дек. 2004. - Томск: Изд. ТПУ, 2004. - с. 189-191.
7. Ромашова О.Ю., Габидуллин О.Р. Способ повышения эффективности ступенчатого подог рева сетевой воды на ТЭЦ при совместной работе турбоустановок //Энергетика: Экология, надежность, безопасность: Материалы докладов десятой Всероссийской научно технической конференции - Томск, 8-10 дек. 2004 г. - Томск: Изд. ТПУ, 2004. - с. 222-225.
8. Ромашова О.Ю., Жмакин Д.А. Исследование режимов работы турбоустановки ПТР-80 130/13 в составе ТЭЦ с поперечными связями на математической модели //Энергетика: Эко логия, надежность, безопасность: Материалы докладов десятой Всероссийской научно технической конференции - Томск, 8-10 дек. 2004. - Томск: Изд. ТПУ, 2004. - с.
9. Ромашова О.Ю., Мухаммадеев К.М. Выбор расчетных параметров отопительных отборов противодавленческих турбин //Энергетика: Экология, надежность, безопасность: Материа лы докладов десятой Всероссийской научно-технической конференции - Томск, 8-10 дек.
2004 г. - Томск: Изд. ТПУ, 2004. - с. 225-228.
10. Ромашова О.Ю. Способ получения пиковой мощности на промышленно-отопительных ТЭЦ //Электроэнергия и будущее цивилизации: Материалы Международной научно-технической конференции - Томск,. - Томск: Изд. ТПУ, 2004. - с. 438-440.
11. Ромашова О.Ю., Мухаммадеев К.М. Перевод турбины ПТР-80-130/13 на подогрев подпи точной воды теплосети //Энергетика: экология, надежность, безопасность: Материалы док ладов XI Всероссийской научно-технической конференции - Томск, 7-9 дек. 2005. - Томск:
Изд. ТПУ, 2005. - с. 242-245.
12. Ромашова О.Ю., Волков В.А. Распределение отопительной нагрузки на ТЭЦ //Энергетика:
экология, надежность, безопасность: Материалы докладов XI Всероссийской научно технической конференции - Томск, 7-9 дек. 2005. - Томск: Изд. ТПУ, 2005. - с. 245-247.
13. Мухаммадеев К.М., Ромашова О.Ю. Выбор схемы и параметров теплофикационных турбин при проектировании на подогрев подпитки //Современные техника и технологии: Труды XII Международной научно-практической конференции студентов, аспирантов и молодых уче ных в 2-х т. - Томск, 27-31.03.2006. - Томск: Изд. ТПУ, 2006. - с. 389-392.
14. Беляев Л.А., Ромашова О.Ю., Волков В.Н. Распределение тепловых нагрузок между турби нами ТЭЦ СХК методом динамического программирования //Энергетика: экология, надеж ность, безопасность: Материалы докладов двенадцатой Всероссийской научно-технической конференции - Томск, ТПУ, 6-8 дек. 2006. - Томск: Изд. ТПУ, 2006. - с. 220-222.
15. Л.А. Беляев, А.А. Боярский, О.Ю. Ромашова. Выбор оптимального состава турбоагрегатов ТЭЦ на период. // Материалы двенадцатой Всероссийской научно-технической конферен ции «Энергетика: экология, надежность, безопасность». – Томск: Изд-во ТПУ, -2006. -465 с., с. 216-220.
16. Боярский А.А., Беляев Л.А., Ромашова О.Ю., Мухаммадеев К.М. Выбор оптимальных ре жимов работы турбоустановки ПТР-80-130/13 //Повышение эффективности производства и использования энергии на Дальнем Востоке: Материалы IV семинара Вузов Сибири и Даль него Востока по теплофизике и теплоэнергетике - Владивосток, ДВГТУ, 6-9 сент. 2005. Владивосток: Изд-во ДВГТУ, 2006. - с. 183-190.
17. Боярский А.А., Беляев Л.А., Пушкин С.В., Ромашова О.Ю. Особенности распределения на грузок на ТЭЦ с поперечными связями //Повышение эффективности производства и исполь зования энергии на Дальнем Востоке: Материалы IV семинара Вузов Сибири и Дальнего Востока по теплофизике и теплоэнергетике - Владивосток, ДВГТУ, 6-9 сентября 2005. Владивосток: Изд-во ДВГТУ, 2006. - с. 190-201.
Подписано к печати 27.04.2007. Формат 6084/16. Бумага «Классика».
Печать RISO. Усл. печ. л. 1,16. Уч.- изд. л. 1,05.
Заказ. Тираж 100 экз.
ИЗДАТЕЛЬСТВО ТПУ. 634050, г. Томск, пр. Ленина, 30.