автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.06, диссертация на тему:Оптимальное управление режимами работы ТЭС Вьетнама с учетом экологических факторов

На правах рукописи

ЛЕКУАНГХОА

ОПТИМАЛЬНОЕ УПРАВЛЕНИЕ РЕЖИМАМИ РАБОТЫ ТЭС ВЬЕТНАМА С УЧЕТОМ ЭКОЛОГИЧЕСКИХ ФАКТОРОВ

Специальность 05.13.06 - Автоматизация и управление технологическими процессами и производствами (отрасль: энергетика)

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

Москва, 2003 г.

Работа выполнена на кафедре "Автоматизированные системы управления тепловыми процессами" Московского энергетического института (технического университета).

Научные руководители: доктор технических наук,

Ведущая организация: Ивановский государственный энергетический

Защита диссертации состоится "27" ноября 2003 г. в 15 час. 00 мин. в аудитории Б-205 на заседании диссертационного совета Д 212.157.14 при Московском энергетическом институте (техническом университете) по адресу: г. Москва, Красноказарменная ул., дом 17.

Отзывы в двух экземплярах, заверенные печатью организации, просим присылать по адресу: 111250, Москва, Красноказарменная ул., дом 14. Ученый Совет МЭИ.

С диссертацией можно ознакомиться в библиотеке МЭИ (ТУ).

Автореферат разослан "_"_ 2003 года

Ученый секретарь

диссертационного Совета Д 212.157 14

профессор Аракелян Эдик Койрунович; доктор технических наук, Нгуен Ван Мань.

Официальные оппоненты: доктор технических наук,

профессор Балакирев Валентин Сергеевич;

кандидат технических наук,

профессор Марченко Евгений Михайлович

университет.

к.т.н, доцент

Буров В.Д.

1 бою

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы:

Одной из актуальных задач оперативного управления в процессе эксплуатации тепловых электростанций является задача оптимального распределения нагрузок (ОРН) между параллельно работающими агрегатами.

Многочисленные исследования, направленные на решение задачи ОРН базируются на поиске минимальных затрат расхода топлива. Вместе с тем в последнее время широко обсуждаются учет других немаловажных факторов, таких как надёжность работы оборудования, экологические ограничения и учет реального технического состояния оборудования и т.д.

Современные условия работы ТЭС отличаются ужесточением требований к предельно допустимым выбросам в атмосферу, введением платы за выбросы. Учет экологических факторов особенно важен при неблагоприятных метеорологических условиях. Все эти факторы делают актуальной задачу регулирования вредных выбросов в темпе с изменением режимов по нагрузке и топливу.

Учет факторов экологии усложняет эту задачу в связи с большим числом переменных. При этом необходимо решить многокритериальную задачу оптимизации, т.е. критериями оптимизации являются минимум расхода топлива и минимум выбрасываемых вредных веществ в окружающую среду.

Во Вьетнаме до сих пор при эксплуатации электростанций не уделялось внимания оптимальному распределению нагрузки между энергоблоками. ТЭС Вьетнама работают с равным распределением нагрузки между энергоблоками, что приводит к низкой экономичности при эксплуатации электростанций. Кроме того, для ТЭС Вьетнама экологический фактор остается острой проблемой. Все угольные ТЭС Вьетнама расположены в регионах с высокой плотностью населения, заведомо напряженной экологической обстановюЛ1.ндцИОНдЛЬцд;51"

Указанное обстоятельство вызывает необходимость решения задачи оптимального распределения нагрузки между энергоблоками ТЭС Вьетнама с учетом экологических факторов. Цель работы:

Целью диссертационной работы является разработка методических положений оперативного управления режимами работы ТЭС с учетом экологических факторов, и их реализация применительно к ТЭС Вьетнама. Методы исследования:

При решении поставленных задач использованы основные положения теории вероятности, статистических методов построения функциональных зависимостей, методы вычислительной математики и обработки экспериментальных данных, современные математические методы оптимизации для решения многокритериальной задачи оптимизации, вычислительные программы на ЭВМ. Научная новизна:

- разработана методика определения энергетических и экологических характеристик энергоблоков на основе аппроксимации моделей кусочно-полиномиальными функциями;

- разработана методика решения многокритериальной задачи оптимизации по разным методам: главного фактора, весовых коэффициентов и минимаксного подхода;

- разработан метод решения многокритериальной задачи оптимального распределения нагрузки ТЭС на основе применения оврагоперешагового алгоритма аффинных проекций;

- разработана комплексная программа для идентификации энергетических характеристик энергоблоков и решения многокритериальной задачи оптимизации распределения нагрузки между энергоблоками станции;

- впервые проведена полная оценка влияния вредных выбросов ТЭС Вьетнама на биосферу и определены энергетические и экологические характеристики ТЭС Вьетнама на основе статистических фактических

"■V 4

данных и особенностей эксплуатации энергоблоков. На их основе реализовано решение многокритериальной задачи оптимизации распределения нагрузки между энергоблоками ТЭС Вьетнама с учетом экологических факторов.

Степень достоверности:

Предложенные методики и полученные результаты обоснованы теорией оптимизации и практической эксплуатацией энергоустановок.

Практическая ценность работы:

Предложенные методики позволяют решить задачу идентификации энергетических характеристик энергоблоков и многокритериальную задачу оптимизации режимов работы теплоэлектростанции.

Предложенные методы расчетов, разработанные программы и полученные результаты могут быть внедрены в практической эксплуатации ТЭС Вьетнама.

Апробация работы:

Основные положения и результаты диссертационной работы докладывались и обсуждались на заседании кафедры АСУ ТП МЭИ, на Международной научной конференции «Control - 2000» (Москва, 2000 г.), Международном конгрессе «Энергетика - 3000» (Обнинск, 2000 г.), Международном 6-ом научном симпозиуме Вьетнамской научно-технической ассоциации в РФ (Москва, 2002 г.)

Публикация:

По теме диссертации опубликованы 6 работ, в том числе статьи в журнале «Вестник» МЭИ.

Структура и обьем работы:

Диссертация состоит из введения, трех глав, заключения, списка литературы, состоящего из 92 наименований, и приложения.

Общий объем 156 страниц включает 126 страниц текста, 28 рисунков и 36 таблиц.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Во введении обосновывается актуальность темы диссертационной работы, формулируются ее цели и задачи, излагается краткое содержание глав диссертации.

В первой главе проведены анализ особенностей графика потребления энергии, в том числе графики нагрузки угольных ТЭС Фалай (Вьетнама) и анализ основных особенностей общей задачи оптимизации режимов работы теплоэнергетических установок в современных и реальных условиях.

Сформулирована в наиболее общей форме многокритериальная задача оптимального распределения отпускной электрической и тепловой энергии нагрузки между параллельно работающими блоками теплоэлектростанции. В общем случае эта задача выражается следующим образом:

Fi (x)=Zfu(Pj^j) -> min

J=1 х

П

Fi (х) — min

J-1 х

П

Fs (х) f3i(pJ,qj) -у min

J=1 x

(1.1)

при условиях:

I Pj=Pr, (1.2)

J-l г*

„ min „ ^ „ max min -- max e1™" ! — 10 „ il 1\

Pj < Pj<Pj , q} <qj<qj ,PjiFqj<Sj , j - 1,2,..., n, (1.3) где Fi(x), F2(x), F}(x) - критерии оптимальности по расходу топлива, надежности и экологии; f^ (pj,qj), f (pj.qj), f (pj.qj) - характеристики j-го

энергоблока по расходу топлива, надежности и экологии; рр qj -электрическая и тепловая мощность j-ro блока; Pi, - суммарные требуемые - электрическая и тепловая мощности, заданные для всей системы; р™т, дтах, q/nn, q™**, 5jraax - допустимые пределы нагрузки; х = IP.q} ~ {xi,...,x„y, р = {р!,...,р„}\q = {qi,...,q„}.

В этой главе проведен анализ математических методов оптимизации, их достоинств и недостатков при решении задачи оптимизации АСУ ТП на

практике. На основе анализа пригодности и эффективности существующих математических методов оптимизации предлагается применение оврагоперешагового алгоритма аффинных проекций (ОАГТ) автора Нгуен Ван Мань (Вьетнама), как эффективного метода для решения многокритериальной задачи оптимизации режимов работы ТЭС. Сущность алгоритма ОАП состоит в следующем: Оврагоперешаговый метод итерационной оптимизации строится на основе понятий направления изменения целевой функции и правила регулировки шага. Общая схема метода описывается следующим уравнением: ж*"1"1 = хк+ ак+18к, к = 0,1,..., (1.4)

где хк, х|с+1 - начальная и конечная точки к+1-й итерации; як - вектор направления изменения переменных в и-мерном пространстве Е"; осм -длина шага итерации.

Существенное отличие оврагоперешагового метода от известных методов безусловной оптимизации состоит в правиле регулировки шага, согласно которому начальная и конечная точки на каждой итерации находятся на двух сторонах точки минимума функции в поисковом направлении. Такое правило оптимизации называется "принципом перешагивания", а длина шага итерации, определяемая по принципу перешагивания, называется "оврагоперешагом".

Пусть имеется минимизируемая функция У(х). Обозначим функцию одной переменной относительно к+1-й итерации:

к(а) = Дхк + а як), (1.5)

где х* - начальная точка данной итерации; як - улучшающее направление движения поисковой точки; а - шаг итерации. Производная от (1.5) равна:

=7"/(хк = (у-/(хк),8к)<0. (1.6)

Неравенство (1.6) означает, что в окрестности точки «= 0 функция к(а) является убывающей по а, таким же образом, как целевая функция убывает по направлению вк. Если Дх) ограничена снизу, то также и И(а) ограничена снизу. При этом, существует по крайней мере одно такое значение а, что: а = аг§гшп/г(а), (1.7)

а>0

где а - точка некоторого локального минимума функции И(а).

Откуда, по принципу перешагивания (в направлении вк), можно записать условие определения оврагоперешага «сп:

А'(«)|^>0, ЦаП^ЦО). (1.8)

Процесс перешагивания на каждой к-й итерации происходит следующим образом (рис. 1.1).

МО)

Рис. 1.1. Определение оврагоперешага И(а°п) сГ; А(а)=-Дхк+а8к).

Ыа\ О

При соблюдении условия (1.8) поисковая точка перемещается от начальной точки хк к конечной х1с+1 по итерационному уравнению: хк+1 = хк+ ак+18к, где ак+[ = сГ.

Направление поиска во оврагоперешаговом алгоритме аффинных проекций (ОАП) определяется по перпендикуляру пирамиды, образуемой из г антиградиентов предыдущих итераций, следующем образом:

Пусть, вк является линейным сочетанием г антиградиентов, т.е.

8к=1Д,Ук-'+|, г<п, (1.9)

1=1

где V1"' = - У1(хк"1 ) - антиградиент функции на (к4+1) -й интеграции; г = к+1, если к <п, иначе г = п, Л;- весовые коэффициенты.

Итерационный процесс алгоритма ОАП поясняется следующим образом (рис. 1.2).

Рис. 1.2. Геометрическая интерпретация процесса минимизации по алгоритму

ОАИ. V1 = -У./(х'), 5° = У°.

XV V1 V1'

Начиная с точки х°, поисковая точка перемешается по антиградиенту

8°=У° со оврагоперешагом а\ до точки х1 с антиградиентом V1.

В точке х1 вычисляется антиградиент V1 и строится направление в1 как перпендикуляр треугольника с вершиной в точке х1 и основанием У'-У°. По полученному направлению в1, с точки х1 выполняется перешагивание до точки х2. Запоминается еще антиградиент V1 для построения пирамиды.

В точке х2, определяется антиградиент V2, который вместе с У',У° образует пирамиду вершины х2. Строят перпендикуляр этой пирамиды, получая вектор 82 с началом в ее вершине и концом на ее основании. По направлению в2 опять осуществляется поиск с перешагиванием до получения х3. На этом шаге дополняется еще антиградиент V2 в память для построения пирамиды в следующих итерациях.

Процесс выполняется таким образом до к = и-й итерации. Дальше, при к > и, на каждой итерации сохраняются только и последних антиградиентов, необходимых для построения пирамиды и ее перпендикуляра. Итерационный процесс продолжается до тех пор, пока не будут удовлетворены критерии остановки.

Во второй главе предложена методика идентификации энергетических и экологических характеристик энергоблоков теплоэлектростанции.

Энергетические и характеристики энергоблоков являются весьма сложными и нелинейными. В зависимости от конкретных технологических и эксплуатационных условий эти характеристики обычно имеют несколько диапазонов, каждый из которых характеризуется отдельной гладкой функцией. На практике статическую характеристику энергоблоков можно с

высокой точностью аппроксимировать последовательным соединением конечного числа кусочно-гладких функций:

ЯЛ А) = £/Ч(Р,А),

(2.1)

4=1

где Р - нагрузка энергоблока;/^,А) - характерная функция в диапазоне Бф А) = 0, если РйБч; - q-й диапазон изменения переменных; С? - число отличающихся диапазонов; А = {аиа2,...,ач} - вектор коэффициентов.

т,

Рис. 2.1. Характеристика энергоблока с кусочно-гладкой моделью. Р - переменная нагрузки

ф .

/¡т

01 О] Оз

Анализ показал вполне удовлетворительным кусочно-полиномиальную модель второго порядка с тремя диапазонами (рис.2.1):

'/¡(Р),еслиР</>, ЯЛА)= /2(ПеслиР,</><;Р2, (2.2)

/3(Р),еслкР>Р2,

где /¡(А^Р^&Р) - функции, характеризующиеся в 1-м, 2-м и 3-м диапазонах соответственно; Рщт, Ртах - нижние и верхние пределы общей области изменения нагрузки; Р\, Р2 - границы между диапазонами.

Вектор искомых коэффициентов А = {а], а2,..., ач} определяется минимизацией среднеквадратичного отклонения:

т

ДА^ВЯЛ'и)-/»]2

тт,

(2.3)

где У - статистическое значение характеристики при нагрузке ?(1) с ограничениями:

Для минимизации функции (2.3) особенно эффективен оврагоперешаговый алгоритм аффинных проекций (ОАП). На основе статистических данных, решена задача определения энергетических и экологических характеристик энергоблоков для ТЭС Фалай Вьетнам со следующими ограничениями:

число блоков: п = 4;

60 Мвт = Р(!)тш < Р(,) < Р("тах =110 Мвт, и получены следующие энергетические и экологические характеристики энергоблоков:

Расходно-топливные характеристики:

- Блок 1: В,(Р) = 19,144 - 0,609 Р + 22,624 Р2 ;

- Блок 2: В2(Р) = 39,913 - 46,376 Р + 53,723 Р2;

- Блок 3: В3(Р) = 8,412 + 20,082 Р + 11,339 Р2;

- Блок 4: В4(Р) = 46,123 - 57,952 Р + 59,830 Р2, Характеристики по выбрасьшаемой золе:

- Блок 1: Gi(P) = 0,135 - 0,129 Р + 0,104 Р2;

- Блок 2: G2(P) = 0,094 - 0,010 Р + 0,087 Р2;

- Блок 3: G3(P) = 0,129 - 0,199 Р + 0,184 Р2;

- Блок 4: G4(P) = 0,071 + 0,103 Р + 0,055 Р2. Графики этих характеристик показаны на рис. 2.2,2.3, 2.4 и 2.5.

ОТ Об ОЛ I

Рис.2.2. Расходная характеристика блоков 1 и 2 Рис.2.3. Расходная характеристика блоков 3 и 4 I А . статистические значения блоков 1 и 2 ■ ^ - статистические значения блоков 3 и 4

01 us ол

Рис.2.4. Характеристика по золе блоков 1 и 2 ■ А . статистические значения блоков 1 и 2

Рис.2.5. Характеристика по золе блоков 3 и 4 I А . статистические значения блоков 3 и 4

В этой главе также проведено исследование влияния вредных выбросов ТЭС Вьетнама на биосферу:

- проведен анализ экологических особенностей ТЭС Вьетнама, в том числе характеристик топлива;

- проведены расчеты количества выбрасываемых вредных веществ угольных ТЭС в атмосферу и распределения их концентраций в атмосфере в условиях климата Вьетнама. Распределение выбросов определяется на основе применении метода Паскуилла - Гиффорда:

А, = —£-_ехр [- - fR. (2-4)

лауа2и 2сг у 2сгг

Результаты расчета концентрации вредных веществ выбросов и графики их распределения для угольных ТЭС Фалай Вьетнама в приземном слое воздуха на оси струи приведены в табл. 2.1 и рис. 2.6.

Полученные результаты показывают, что концентрация кратковременных выбросов угольных ТЭС Вьетнама (в виде дымовых газов, содержащих золы, оксиды, серы, азота) в приземном слое воздуха в населенных пунктах на расстоянии 2 - 4 км от ТЭС в 1,5 - 8 раз больше нормы. Вследствие радиоактивных выбросов естественный радиоактивный фон в населенных пунктах на расстоянии 2 - 4 км от ТЭС в 20 - 40 раз выше, чем средний фон по стране.

Таблица 2.1

Расстояние, м Твердые частицы, г/м' Радионуклиды, Бк/м3 вСЬ, г/м3 СО, г/м3 N0;, г/м3

300 1,3.10'14 2,1.10 м 7,8.10"'4 1.1.10-'3 4,5.10""

500 5,3.10"ш 8,5.10 ю 1,4.10-® 1,9.10"® 7,8.10""

700 2.3.10-05 3,7.10"" 5,9.10-® 8,2. Ю"05 3,4.10"06

900 1,2.10-" 1,9. Ю-04 3,1. ю-* 4,3.10"04 1,8.10"05

1,0Е + 3 4,6.10^ 4,2.10"04 7,7.10"04 1,4. Ю"03 6,9.10 ю

2,0 Е + 3 6,2.10-" 5,1.10^ 10,9.10-" 8,6.10-® и.10-04

З.ОЕ + З 6,0.10 м 3,2.10-™ 8,3. Ю"04 7,4.10"03 1,5.10"°*

5,0 Е + 3 З.З.Ю"04 2,0 Ю"04 5,3. Ю-04 1,6. Ю^3 5Д10-05

7,0 Е + 3 7,9.10 м 1,3.10"04 2.1.10-04 ЗДЮ"01 1,2.10-®

9,0 Е + 3 5,9.10"05 9,5. Ю-05 1,6.10^ 2,2.10-04 9,2.10-®

1,0 Е +4 5,0. Ю-05 8,0.10"05 1,3.10"04 1ДЮ"04 7,7.10-*

2,0 Е + 4 1,9.10-05 ЗДЮ"05 5,2. Ю-05 7,3.10"™ ЗДЮ"06

3,0 Е1-4 1,2.10 м 1,9.10-® ЗДЮ"05 4,8. Ю"05 2,0.10"06

5,0 Е 4 4Д.10"06 7Д10"06 1,4.10^ 1,9.10-" 7,9.10-°'

7.0Е + 4 гл.ю-06 ЗДЮ"06 6,8.10-® 9,5. Ю"06 ЗДЮ"07

9,0 Е + 4 ьгло-04 2,0. Ю"06 4ДЮ-06 5,8.10"® гдт-07

1,0 Е + 5 5,7.10"07 1,4.10"ое ЗДЮ"06 4,2.10"06 1,7.10"07

Рис. 2.6. Распределение концентрации вредных выбросов ТЭС "Фялай" в приземном слое воздуха.

ЬовА»

I"1

5-.

со

да/

|/ " N02 ~"

Н твердые

Л/ частхцы -

Расстояниях^ Ьовх

В третьей главе приведена общая методика решения многокритериальной задачи оптимального распределения нагрузки с учетом экологических факторов между энергоблоками станции. Такая задача формулируется в следующем виде:

Г=1 1 .

Х={р1.....ра},

(3.1)

при условиях:

(3.2)

(3.3)

(*ГЛ>.И>2.....п),

>■1

„тт - ;_1 „

где п - число подсистем (энергоблоков относительно станции или станций относительно энергосистемы); ^(х), ^(х), - критерии по расходу топлива, экологии; /у(ру), /2\(Р]), - характеристики .¡-ой подсистемы по расходу топлива и выбросам; р} - электрическая нагрузка ;-ой подсистемы;

их допустимые пределы; - заданная суммарная

электрическая нагрузка системы в целом; х = {/зь..., рп) = {хи..., *„} -переменные нагрузки.

Для решения задачи (3.1) - (3.3), как правило, преобразуют векторный критерий (3.1) в единый скалярный вид, что приводит в результате к однокритериальной задаче минимизации с ограничениями. В работе использованы наиболее употребляемые в инженерной практике подходы скаляризации векторного критерия оптимальности. Такими являются

„тт ш»

Р} >Р\

подходы: "главного критерия", "весовых коэффициентов" и "минимаксный подход".

- По подходу "главного критерия", из всех критериев оптимальности выбирается один самый важный, а остальные включаются как ограничения. В результате, исходная задача преобразуется в форму:

№ = Ш —► min, (3.4)

X

при дополнительных ограничениях: f, (я) < ff, гфт, (3.5)

где х ={р\,...,ра}, fm (х) - главный критерий;^ (t£m) - допустимые пределы ухудшения неглавных критериев.

- По подходу "весовых коэффициентов" составляется компромиссный критерий, состоящий из суммы всех критериев, умноженных на свои соответствующие весовые коэффициенты. Форма задачи при этом приводится к виду:

m=tr,m —► min, (з.б)

<=1

где yt - положительные весовые коэффициенты.

- По минимаксному подходу максимальное из значений всех критериев принимается каждый раз за общее, относительно которого осуществляется минимизация. Такая задача оптимизации записывается в форме:

/(х)= max ft (х) —► min, (3.7)

Q

при ограничениях: Y.P}=ps> (*j = ft>j=1.2>-,n)> (3-8)

j-i

j = l,...,n, (3.9)

Для решения многокритериальной задачи оптимального распределения

нагрузки между энергоблоками ТЭС предлагается метод решения на

основе применения алгоритма ОАП. Метод состоит из следующих шагов:

1. Нормализация критериев оптимальности

Отдельные критерии имеют различные физические сущности с разными размерностями и отличаются по величине. Для решения задачи

сначала осуществляется нормализация критериев, т.е. приведение их к стандартному безразмерному виду, в пределе от 0 до 1, по формуле: F(x)-Fmm Fix) Fmm

\ =Fr-Fr, (3.10)

1 1 "l 1 где ^т,п,^гаах - нижний и верхний пределы i-ro критерия; //¡(х)его нормализованная форма.

Откуда вместо (3.1) получается эквивалентная задача:

Я, (х), Н2 (х) -> min (3.11)

2. Формирование единого критерия оптимальности

Формируют единый критерий по весовому подходу, суммируемый из частных критериев с весом:

/(х)= (х)+угН2 (х) min, (3.12)

X

где у\, уг~ весовые коэффициенты, определенные на основе экспертных оценок.

- По минимаксному подходу, единый критерий определяется по соотношению:

/ (х) = max Я, (х) -> min (3.13)

i

- По подходу "главного фактора" выбирается наиболее "важный", а для остальных устанавливают определенные пороги допустимого ухудшения:

/(х) = Ят (х) -*■ min, (3.14)

с включением дополнительных ограничений:

Я,(х) < Hf, 1 < i < 2, i * m, (3.15)

где Я|Д, - допустимый предел ухудшения i-ro критерия

3. Формирование эквивалентной безусловной целевой функции

/(х)min, х= {хих2,...,х2п}, (3.16)

х

при ограничениях:

Я,(х)<Я,д, l<i<2, i*m, (3.17)

(3.18)

Х™<Х1<Х™, ] = 1,2,...,п, (3.19)

где ^ = Рр х™'т = рГ", х™ = р™, j =1, 2,..., п. Ограничения (3.17) имеют место только в подходе главного критерия.

На основе штрафных функций для решения задачи (3.16) - (3.19) формируется эквивалентная целевая функция:

Ах) = Дх) + с^ (х) + сЛ (х)+сх (ж)], (3.20)

где Сх, сх, Ся - штрафные коэффициенты;

• штрафная функция по ограничению ухудшения неглавных критериев:

^я(х) = 1[|Я,(х)-ЯЛ+Я,(х)-Я,д]; (3.21)

1=1

• штрафная функция по дисбалансам между производством и потреблением:

.1=1

• штрафная функция по ограничению нагрузок подсистемы:

^ (х)=£ ки - *г+1 ^ - хг I ■+1 ^ - *г |}- (3-23)

Для минимизации функции (3.20) принимается оврагоперешаговый алгоритм ОАП. Решение задачи оптимального распределения нагрузки между, блоками ТЭС Вьетнама с учетом экологических факторов реализуется в двух вариантах: при режиме распределения нагрузки по всем существующим работающим энергоблокам и при отключении одного из блоков.

Результаты решения задачи при режиме распределения нагрузки по всем существующим работающим энергоблокам для ТЭС Фалай Вьетнама, состоящей из четырех энергоблоков с установленной мощностью 440 Мвт, приведены в табл. 3.1. Значения расхода топлива и выбросов при оптимальном и фактическом распределении нагрузки приведены в табл. 3.2.

Результаты решения задачи при отключении одного из энергоблоков с суммарной мощностью 256 Мвт приведены в табл. 3.3.

Таблица 3.1

Сумм, мощность, (Мвт) Весовой подход Г\=П =0,5 Минимаксный подход Подход "главного фактора" (расхода топлива)

Блок 1 Блок 2 БлокЗ Блок 4 Блок 1 Блок 2 БлокЗ Блок 4 Блок I Блок 2 Блок 3 Блок 4

248 68 60 60 60 68 60 60 60 62,34 60 65,66 60

256 71,32 64,46 60,22 60 74,55 60 61,45 60 66,59 60 69,41 60

272 81,59 63,22 66,05 61,14 83,51 60,89 67,6 60 67,51 71,03 73,46 60

292 86,97 69,05 71,9 64,08 83,4 70,63 69,53 68,44 77,14 75,08 62,69 77,09

312 93,32 72,92 77,74 68,02 82,41 77,31 73,23 79,05 82,40 77,3 73,2 79,1

336 100,98 77,49 84,77 72,76 88,73 79,98 85,81 81,48 88,76 79,97 85,81 81,46

348 104,82 79,77 88,29 75,12 91,9 81,28 92,14 82,68 91,91 81,29 92,13 82,67

364 109,92 82,68 92,94 78,46 96,1 83,07 100,53 84,3 96,13 83,07 100,53 84,27

380 110 87,3 99,83 82,87 100,32 84,86 108,96 85,86 100.34 84,85 108,96 85,85

400 110 92,92 108,41 88,67 110 89,75 110 90,25 110 89,75 110 90,25

Таблица 3.2

Суммарная мощность, Мвт Весовой подход У\ - /2 =0,5 Минимаксный подход Подход главного фактора по Фактические данные

Расход топлива В (т/ч) Выбросы G (т/ч) Расход топлива В (т/ч) Выбросы G(т/ч) расходу топлива выбросам Расход топлива В (т/ч) Выбросы G (т/ч)

В (т/ч) G (т/ч)

248 117,8082 0,44872 118.0281 0,44683 118,0184 0,44613 118,215 0,450

256 121,1900 0,44986 j 121,6215 0,44474 122,0488 0,44875 123,252 0,454

272 124,0512 0,45471 123,8071 0,45112 123,5113 0,45097 126,681 0,467

292 130,5237 0,47067 130,6786 0,47193 129,9828 0,46888 133,880 0,486

312 137,2235 0,49015 137,0309 0,49132 137,0746 0,48932 140,818 0,506

336 145,0738 0,51480 145,2282 0,52023 145,0262 0,51970 149,116 0,534

348 150,0254 0,53157 149,8635 0,53255 150,1153 0,52959 153,163 0,549

364 157,3454 0,55713 157,7725 0,55963 157,2682 0,55697 160,306 0,570

380 164,1920 0,57835 164,7822 0,58002 164,7777 0,57992 166,921 0,593

400 175,3679 0,61085 178,2176 0,61562 174,9976 0,61538 177.104 0,623

Таблица 3.3

Суммарная мощность, Мвт Блок 1 Л, Мвт Блок 2 p¡, Мет БлокЗ pi, Мвт Блок 4 р4, Мвт При оптимальной нагрузке Фактические данные

в (т/ч) | G(t/4) в (т/ч) | g (т/ч)

Весовой подход, yi = = 0,5

256 97,85 76,07 82,08 0 118,21 0,432 123,25 0,454

103,52 79,43 0 73,05 122,45 0,442

100,34 0 84,47 71,19 119,70 0,440

0 83,92 94,61 77,47 122,26 0,445

Минимаксный подход

256 89,18 80,15 86,67 0 118,62 0,438 123,25 0,454

91,98 81,31 0 82,71 121,96 0,444

88,73 0 85,79 81,48 119,16 0,439

0 81,29 92,08 82,63 122,06 0,447

Подход главного фактора (расхода топлива)

256 89,28 80,11 86,61 0 118,86 - 123,25 0,454

92,02 81,29 0 82,69 121,84 -

89,13 0 85,69 81,18 119,07 -

0 80,88 91.74 83,38 121,76 -

Подход главного фактора (выбросов)

256 НО 64,07 81,93 0 - 0,433 i 123,25 | 0,454

110 86 0 60 - 0,448

107,19 0 | 86,69 62,12 - 0,439

0 96,15 | 92,2 62,65 - 0,446

Из полученных результатов видно, что техническая эффективность метода оптимизации зависит не только от подходов решения задачи, но и от рабочих нагрузок ТЭС или от относительной мощности ТЭС: АР = (Р/Ртш)-Значения 0,6 < АР < 0,8 дают наилучшие результаты оптимизации по расходу топлива и выбросам. Для ТЭС Фалай Вьетнама, при оптимальном распределении нагрузки по всем работающим энергоблокам расход топлива уменьшается на 1% - 2,9% и количество выбросов сокращается на 1% -3,6%.

- Если нагрузка станции меньше 60% установленной мощности, следует выбрать вариант с отключением четвертого блока при решении задачи оптимального распределения нагрузки. Эти оптимальные решения приносят более 3% экономии расхода топлива и более 3,5% сокращения количества выбросов

В приложении описан пакет комплексной программы для идентификации энергетических характеристик энергоблоков и решения многокритериальной задачи оптимизации распределения нагрузки между энергоблоками станции.

ЗАКЛЮЧЕНИЕ

Основные результаты диссертационной работы заключаются в следующем:

1. Проведен анализ основных особенностей задачи оптимального управления режимами работы теплоэнергетических установок в современных и реальных условиях. Сформулирована в общей форме многокритериальная задача оптимального распределения отпускной электрической и тепловой энергии нагрузки между параллельно работающими блоками теплоэлектростанции.

2. Проведены анализ математических методов решения задачи оптимизации и сравнение их достоинств и недостатков при решении задачи оптимизации АСУ ТП на практике. На основе анализа пригодности и эффективности существующих математических методов

оптимизации предлагается применение оврагоперешагового алгоритма аффинных проекций (ОАП) как эффективного метода для решения многокритериальной задачи оптимизации режимов работы ТЭС.

3. На основе анализа математических моделей идентификации объектов и особенностей эксплуатации энергоблоков, предложен эффективный алгоритм определения энергетических характеристик энергоблоков ТЭС путем моделирования характеристик энергоблоков кусочно-гладкой функцией, на основе которого определены энергетические и экологические характеристики энергоблоков ТЭС Вьетнама.

4. Проведен анализ влияния вредных выбросов ТЭС Вьетнама на биосферу на основе расчетов количества выбрасываемых вредных веществ и распределения их концентраций в атмосфере.

5. Предложена общая методика решения многокритериальной задачи оптимизации. Построен принцип решения многокритериальной задачи оптимизации на основе трех способов скаляризации: "весового суммирования", "минимаксного определения" и "главного фактора".

6. Предложен эффективный метод решения многокритериальной задачи оптимального распределения нагрузок между параллельно работающими энергетическими установками на основе применения оврагоперешагового алгоритма.

7. Реализовано решение многокритериальной задачи оптимизации распределения нагрузки между энергоблоками ТЭС Вьетнама с учетом экологических факторов на основе фактических данных по расходу топлива, выбросам и особенностей эксплуатации энергоблоков. При оптимальном распределении нагрузок ТЭС, расход топлива уменьшается примерно на 2%, а количество выбросов - на 3% относительно режима равномерного распределения по всем энергоблокам.

8. Разработана комплексная программа для идентификации энергетических характеристик энергоблоков и решения многокритериальной задачи оптимизации распределения нагрузки между энергоблоками станции.

№16 0 10

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ *

1. Аракелян Э.К., Хоа JI.K. Мань Н.В. Оптимизация режима работы тепловой электростанции по экономическому и экологическому^^^-^ ^ критериям. // Вестник МЭИ - 2002. - № 4. - с. 25-30. | £о(С>

2. Jle Куанг Хоа, Нгуен Хай Нинь. Общая оценка вредных выбросов ТЭС Вьетнама (на английском языке). // Седьмая международная научно-техническая конференция студентов и аспирантов, «Радиоэлектроника, электротехника и энергетика», Москва, 27-28 февраля 2001 г.: Тезисы докладов. - Москва: МЭИ, 2001 - С. 219-221.

3. Ле Куанг Хоа, Нгуен Хонг Вьет. Исследование энергетических характеристик ТЭС Вьетнама. // Восьмая международная научно-техническая конференция студентов и аспирантов «Радиоэлектроника, электротехника и энергетика», Москва, 28 февраля - 1 марта 2002 г.: Тезисы докладов. - Москва: МЭИ, 2002 - С. 219-220.

4. Ле Куанг Хоа, Нгуен Хоиг Вьет. Распределение концентрации выбросов ТЭС Вьетнама в атмосфере (на вьетнамском языке). // Международный 6-ой научный симпозиум Вьетнамской научно-технической ассоциации в РФ, Москва, 24 апреля 2002 г.: Сборник докладов. - Москва: Изд. Творчество, 2002 - С. 278 - 285.

5. Ле Куанг Хоа. Сравнительная оценка радиоактивных выбросов АЭС и угольных ТЭС. // Международный конгресса «ЭНЕРГЕТИКА - 3000», Обнинск, 16-20 октября 2000 г.: Тезисы докладов. - Обнинск: ИАТЭ, 2000 -С. 86-87.

6. Хоа Л. К., Мань Н. В. Оценка радиоактивных выбросов угольных ТЭС во Вьетнаме. // Труды международной научной конференции «CONTROL -2000», Москва, 26-28 сентября 2000 г. - Москва: МЭИ, 2000 - С. 181-184.

Подписано к печати ЛЗ ■ №' л -Печ. Л i. Xb Тираж № Заказ

Типография МЭИ, Красноказарменная, 13. 20

Введение.

Глава 1. Общие положения оптимизации режимов работы теплоэнергетических установок и методы оптимизации.

1.1. Графики потребления энергии.

1.2. Общие положения оптимизации режимов работы теплоэнергетических установок.

1.3. Комплексная оптимизация режимов работы теплоэнергетических установок.

1.3.1. Иерархические свойства задачи оптимизации в энергетической системе.

1.3.2. Оптимизация режимов работы ТЭС.

1.4. Оптимальное распределение нагрузки между энергоблоками ТЭС.

1.5. Характеристики исходной информации. Расходные топливные характеристики электростанции.

1.6. Известные математические методы решения задачи оптимизации.

1.6.1. Общая постановка.

1.6.2. Метод решения задачи оптимизации.

1.6.3. Метод поиска оптимума путем решения системы нелинейных уравнений.

1.6.4. Метод множителей Лагранжа.

• 1.6.5. Метод модифицированной функции Лагранжа.

1.6.6. Метод динамического программирования.

1.6.7. Метод направленного поиска.

1.7. Овраперешаговый метод оптимизации.

1.7.1. Принцип перешагивания поиска минимума.

1.7.2. Общая схема овраперешагового метода.

1.7.3. Овраперешаговый алгоритм аффинных проекций.

1.8. Выводы.

Глава 2. Исследование энергетических характеристик ТЭС Вьетнама и влияния их вредных выбросов на биосферу.

2.1. Идентификация энергетических характеристик энергоблоков теплоэлектростанций.

2.1.1. Аппроксимация линейной комбинацией функций.

2.1.2. Аппроксимация кусочно-гладкой функцией.

2.2. Исследование энергетических характеристик энергоблоков угольных ТЭС Вьетнама.

2.2.1. Настоящее положение эксплуатации ТЭС Вьетнама.

Определение энергетических и экологических характеристик энергоблоков ТЭС.

2.2.2. Оценка ошибки моделей энергетических и экологических характеристик.

2.3. Исследование влияния вредных выбросов угольных ТЭС Вьетнама на биосферу.

2.3.1. Экологические особенности угольных ТЭС Вьетнама.

2.3.2. Определение количества выбросов в атмосферу угольных ТЭС Вьетнама.

2.4. Определение распределения вредных выбросов в атмосфере угольных ТЭС Вьетнама.

2.4.1. Факторы, влияющие на распределение выбросов в атмосфере.

• 2.4.2. Методы расчета распределения выбросов в атмосфере.

2.4.3. Расчет распределения вредных веществ в атмосфере от угольных ТЭС в условиях климата Вьетнама.

2.5. Выводы.

Глава 3. Оптимизация распределения нагрузки между энергоблоками угольных ТЭС Вьетнама с учетом экологических факторов. ф 3.1. Задача оптимизации распределения нагрузки между энергоблоками

ТЭС с учетом экологических факторов.

3.2. Общая методика решения многокритериальной задачи оптимизации.

3.3. Выбор целесообразного математического метода решения комплексной оптимизации.

3.4. Метод решения задачи оптимального распределения нагрузки между энергоблоками ТЭС на основе применения алгоритма ОАП.

3.4.1. Нормализация критериев оптимальности.

3.4.2. Формирование единого критерия оптимальности.

3.4.3. Формирование эквивалентной безусловной целевой функции.

3.5. Решение задачи оптимизации распределения нагрузки между энергоблоками ТЭС Вьетнама с учетом экологических факторов.

3.5.1. Особенности графика нагрузки ТЭС Вьетнама.

3.5.2. Алгоритм решениа задачи оптимизации распределения нагрузки между энергоблоками ТЭС Вьетнама с учетом факторов по выбросам.

3.5.3. Решение задачи оптимизации распределения нагрузки между энергоблоками ТЭС Вьетнама при режиме распределения нагрузки по всем существующим работающим энергоблокам.

3.5.4. Решение задачи оптимизации распределения нагрузки между энергоблоками ТЭС Вьетнама при отключении одного из блоков.

3.6. Сравнительная количественная оценка технической эффективности метода оптимизации режимов работы ТЭС.

Щ 3.7. Оценка значимости коэффициентов и адекватности полученных уравнений регрессии.

3.8. Выводы.

9 Главной и актуальной задачей оперативного управления в процессе эксплуатации тепловых электростанций и электроэнергетических систем является оптимизация режимов работы энергетических установок электростанций. В центре этого процесса стоит задача оптимального распределения нагрузок (ОРН) между параллельно работающими агрегатами (или подсистемами). В любой момент времени в энергосистеме может произойти включение или выключение того или иного потребителя, генератора или линий электропередачи, что приводит к дисбалансу между выработкой электроэнергии и потреблением. Для подержания баланса можно применить плановое изменение нагрузки энергетических единиц (энергоблоков или электростанций) по графику.

Многочисленные исследования, направленные на решение задачи ОРН базируются на поиске минимальных затрат условного топлива или стоимости топлива. Вместе с тем в последнее время широко обсуждаются учет других важных факторов, таких как надёжность работы оборудования, экологические ограничения и учет реального технического состояния оборудования и т.д.

Современные условия работы ТЭС отличаются - ужесточением требований к предельно допустимым выбросам в атмосферу, введением платы за выбросы, повышенной концентрации выбрасываемых вредных веществ в окружающую среду. Учет экологических факторов особенно важен при неблагоприятных метеорологических условиях, способствующих накоплению вредных веществ в приземном слое атмосферы. Все эти • факторы наряду с необходимостью активного участия в поддержании переменной составляющей суточного графика нагрузки делают актуальной задачу регулирования вредных выбросов в темпе с изменением режимов по нагрузке и топливу. Эта задача в настоящее время становится одной на главных для ТЭС, особенно для угольных ТЭС, выбрасывающих большое количество твёрдых и газообразных загрязнителей, в том числе Ф радиоактивных веществ.

Во Вьетнаме до сих пор при эксплуатации электростанций не уделялось внимания оптимальному распределению нагрузки между энергоблоками. ТЭС Вьетнама работают с равным распределением нагрузки между энергоблоками, что приводит к большему расходу топлива и низкой экономичности при равной суммарной мощности электростанции.

Кроме того, для ТЭС Вьетнама экологический фактор остается острой проблемой. Все угольные ТЭС Вьетнама расположены в регионах с высокой плотностью населения, заведомо напряженной экологической обстановкой и вблизи населенных пунктов. Причем ТЭС Вьетнама имеют низкий коэффициент золоуловителей и высокую зольность топлива, поэтому количество летучей золы, выбрасываемой в атмосферу, составляет значительную величину. Экспериментальные данные показывают, что концентрация вредных веществ, выбрасываемых от ТЭС, в населенных пунктах на расстоянии 2 - 3 км от ТЭС в 2 - 7 раз больше нормы.

Указанное обстоятельство вызывает необходимость постановки задачи оптимизации режимов работы ТЭС Вьетнама с учетом экологических факторов. Вследствие отсутствия возможностей режимных мероприятий, направленных на подавление выбросов, воздействие на перераспределение электрической мощности между блоками остается не только повышением экономичности но и единственным способом сдерживания роста выбросов.

При учете факторов экологии, задача оптимального распределения нагрузок (ОРН) между энергоблоками электростанций гласит следующим • образом: при заданной суммарной нагрузке системы требуется определить такие нагрузки энергоблоков, при которых текущие затраты на производство электроэнергии являются наименьшими с наименьшим количеством выбросов вредных веществ в окружающую среду.

Учет факторов экологии усложняет задачу оптимального распределения нагрузки между энергоблоками в связи с большим числом оптимизируемых ф переменных. При этом необходимо решить многоцелевую многокритериальную) задачу оптимизации, т.е. критериями оптимизации являются минимум расхода топлива и минимум выбрасываемых вредных веществ в окружающую среду.

Существующие известные методы решения задачи оптимального распределения нагрузки между энергоблоками в основном базируются на технико-экономических критериях. Применение этих методов для решения многокритериальной задачи оптимизации становится затруднительным.

Целью диссертационной работы является разработка методических положений оперативного управления режимами работы ТЭС с учетом экологических факторов, и их реализация применительно к ТЭС Вьетнама.

Диссертация состоит из введения, трех глав, заключения и приложения.

Основные результаты диссертационной работы заключаются в следующем:

1. Проведен анализ основных особенностей задачи оптимизации режимов работы теплоэнергетических установок в современных и реальных условиях.

2. Сформулирована в наиболее общей форме многокритериальная задача оптимального распределения отпускной электрической и тепловой энергии нагрузки между параллельно работающими блоками теплоэлектростанции.

3. Проведен анализ математических методов решения задачи оптимизации. Проведено сравнение достоинств и недостатков математических методов при решении задачи оптимизации АСУ ТП на практике. На основе анализа пригодности и эффективности существующих математических методов оптимизации предлагается применение «оврагоперешагового» алгоритма аффинных проекций (ОАП) как эффективного метода для решения многокритериальной задачи оптимизации режимов работы ТЭС.

4. На основе анализа математических моделей идентификации объектов и особенностей эксплуатации энергоблоков, предложен эффективный алгоритм определения энергетических характеристик энергоблоков ТЭС путем моделирования характеристик энергоблоков кусочно-гладкой функцией: энергетические характеристики энергоблоков аппроксимируются моделью кусочно-полиномиальной функции и выражаются полиномом не выше второго порядка; вектор искомых коэффициентов полинома определяется путем минимизации квадратических отклонений между моделью и статистическими данными при применении алгоритма ОАП.

5. Проведено определение энергетических и экологических характеристик ТЭС Вьетнама на основе статистических фактических данных и особенностей эксплуатации энергоблоков.

6. Проведена оценка влияния вредных выбросов ТЭС Вьетнама на биосферу: проведение анализа экологических особенностей ТЭС Вьетнама, в том числе характеристик топлива; проведение расчетов количества выбрасываемых вредных веществ в атмосферу, в том числе количества радиоактивных веществ в составе выбросов угольных ТЭС Вьетнама; определение распределения в атмосфере приземных концентраций выбросов вредных веществ угольных ТЭС в условиях климата Вьетнама.

7. Сформулирована общая задача оптимизации распределения нагрузки с учетом экологических факторов между подсистемами относительно энергосистемы или между энергоблоками относительно станции.

8. Предложена общая методика решения многокритериальной задачи оптимизации. Построен принцип решения многокритериальной задачи оптимизации на основе подходов скаляризации единого критерия оптимальности: "весового суммирования", "минимаксного определения" и "главного фактора": по подходу "весовых коэффициентов" (по принципу Парето) составляется компромиссный критерий, как линейная свертка, состоящая из суммы всех критериев, умноженных на свои соответствующие весовые коэффициенты; по минимаксному походу, максимальное из значений всех критериев принимается каждый раз за общее, относительно которого осуществляется минимизация; по подходу "главного критерия", из всех критериев оптимальности выбирают один самый важный, а остальные включаются как ограничения, т.е. для них определяются некоторые пороги и составляются ограничивающие неравенства.

9. Предложен эффективный метод решения многокритериальной задачи оптимального распределения нагрузок между параллельно работающими энергетическими установками на основе применения оврагоперешагового алгоритма. Метод состоит из следующих шагов:

- нормализация критериев оптимальности, то есть проведение критериев к стандартному безразмерному виду, в пределе от 0 до 1;

- формирование единого критерия оптимальности; формирование эквивалентной безусловной целевой функции на основе штрафных функций. Минимизация целевой функции осуществляется применением алгоритма ОАП.

10. Реализовано решение многокритериальной задачи оптимизации распределения нагрузки между энергоблоками ТЭС Вьетнама с учетом экологических факторов по разным подходам на основе фактических данных по расходу топлива, выбросам и особенностей эксплуатации энергоблоков. И. Разработана комплексная программа для идентификации энергетических характеристик энергоблоков и решения многокритериальной задачи оптимизации распределения нагрузки между энергоблоками станции.

Предложенная методика и разработанная комплексная программа могут быть применены для решения задач оптимизации режимов работы не только энергетических установок, но и энергосистем в целом с учетом различных факторов, например, фактора надежности установок и экологического фактора.

ЗАКЛЮЧЕНИЕ

Главным итогом настоящей диссертационной работы является разработка метода решения многокритериальной задачи оптимального распределения нагрузок между параллельно работающими блоками электростанции. Для ТЭС Вьетнама эта задача оптимизации решается с учетом экологических факторов.

1. Александров А.Г. Синтез регуляторов многомерных систем. М.: Машиностроение, 1986.

2. Александров Е.Е. К вопросу параметрической оптимизации линейных регулируемых систем//Изд. Вуз. Электромеханика, 1990. № 6. С. 84-87.

3. Аминов Р.З. Векторная оптимизация режимов работы электростанций. -М.: Энергоатомиздат, 1994

4. Анализ и оптимальный синтез на ЭВМ систем управления /Под ред. A.A. Воронова и И.А. Орурка. -М.: Наука, 1984.

5. Андрющенко А.И., Аминов Р.З. Оптимизация режимов работы и параметров тепловых электростанций. — М.: Высш. школа, 1983

6. Аоки М. Введение в методы оптимизации. — М.: Наука, 1977.

7. Аракелян Э.К., Кормилицын В.И., Самаренко В.Н. Оптимизация режимов оборудования ТЭЦ с учетом экологических ограничений //Теплоэнергетика, 1992. № 2. С. 29-34.

8. Аракелян Э.К., Мань Н.В., Хунг Н.Ч. Оптимальное распределение нагрузки между параллельно работающими энергетическими блоками с учетом фактора надежности //Вестник МЭИ. 1997. № 3. С. 15-20.

9. Аракелян Э.К., Старшинов В.А. Повышение экономичности и маневренности оборудования тепловых электростанций. — М.: Изд. МЭИ, 1993.

10. Бакластов A.A. и др. Алгоритм оптимального распределения нагрузок между параллельно работающими агрегатами //Сб. Науч. Тр. МЭИ, 1987. № 142. С. 61-64.

11. Бакушинский А.Б., Гончарский A.B. Некорректные задачи. Численные методы и приложения. М.: МГУ, 1989.

12. Балакирев В.С., Дудников Е.Г., Цирлин А.М. Экспериментальное определение динамических характеристик промышленных объектов управления. -М.: Энергия, 1967.

13. Банди Б. Методы оптимизации. Перевод с англ. М.: «Радио и связь», 1988

14. Бард И. Нелинейное оценивание параметров. Перевод с англ., под ред. Горского В.И. М.: Мир, 1979.

15. Бахвалов Н. С. Численные методы. М.: Наука, 1973.

16. Бейко И.В. и др. Методы и алгоритмы решения задач оптимизации. — Киев: Вища школа, 1983

17. Бородюк В.П., Вощинин А.П, Иванов А.З. Статистические методы в инженерных исследованиях (лабораторный практикум). Под редакцией Круга Г.К. М.: «Высшая школа», 1983.

18. Бызора НЛ., Гаргев Е.К., Иванов В.Н. Экспериментальное исследование атмосферной диффузии и расчеты рассеяния примеси. -JL: Гидрометеоиздат, 1991.

19. Васильев Ф.П. Численные методы решения экстремальных задач. М.: Наука, 1980.

20. Веников В.А., Журавлев В. Г., Филиппова Г.А. Оптимизация режимов электростанций и энергосистем. М.: Энергоатомиздат, 1990, с.347.

21. Внуков А.К. Защита атмосферы от выбросов энергообъектов. Справочник. -М.: Энергоатомиздат, 1992.

22. Волков Э.П., Лысков М.Г., Фетипова Е.И. Методы расчета приземных концентраций вредных веществ в атмосферном воздухе. М.: Изд. МЕИ, 1991

23. Гаряев А.Б. Распространение опасных веществ при промышленных авариях. М.: Изд. МЭИ, 2000.

24. Гилл Ф., Мюррей Ц., Райт М. Практическая оптимизация. М.: Мир, 1985

25. Гук Ю.Б. Основы надежности электроэнергетических установок. Л.: Энергоатомиздат, 1990

26. Дейч А.М. Методы идентификации динамических объектов. — М.: Энергия, 1979.

27. Демьянов В.Ф. Васильев Л.В. Недифференцируемая оптимизация. М.: Наука, 1981

28. Доклад о исследовании и мониторинге окружающей среды в регионах угольных ТЭС Вьетнама (на вьетнамском языке). Вьетнамский Центр анализа и окружающей среды. Ханой, 11/2000г.

29. Доклад об эксплуатации и эконо-технических параметрах ТЭС Фалай. Вьетнамское ведомство по энергетике (на вьетнамском языке). Ханой: ЕЩ 2/2001

30. Дорофеюк А.А. Алгоритмы автоматической классификации //Автоматика и телемеханика, 1971. № 12.

31. Евтушенко Ю.Г. Методы решения экстремальных задачи их приложения в системах оптимизации. М.: Мир, 1982

32. Зайтендейк Г. Методы возможных направлений. Изд. Иностранной литературы, 1963

33. Защита атмосферы от промышленных выбросов: Справочник./ Под ред. С. Калветрта и Г. Инглумда. Часть 2. М.: Металлургия, 1986

34. Иванов АЛ., Круг Г.К., Филаретов Г.Ф. Статистические методы в инженерных исследованиях. Планирование первого порядка для регрессионных экспериментов. М.: МЭИ, 1978.

35. Карманов В.Г. Математическое программирование. М.: Наука, 1986.

36. Козлов В.Ф. Справочник по радиационной безопасности. — М.: Энергоатомиздат, 1991.

37. Кюнци Г.П., Крелле В. Нелинейное программирование. Изд. «Советское радио», 1965

38. Лавритев О.И., Никифоров А.Д. Аналитический обзор процедур решения многокритериальных задач математического программирования //Экономика и математические методы. 1986. Т. XXII. Вып. 3. С. 508-523

39. Ле Куанг Хоа. Сравнительная оценка радиоактивных выбросов АЭС и угольных ТЭС. Доклад на международном конгрессе "ЭНЕРГЕТИКА -3000", г. Обнинск, Россия. 16-20 октября 2000г.

40. Левенталь Г.В., Попырин Л.С. Оптимизация теплоэнергетических установок. М.: Энергоатомиздат, 1972.

41. Льюнг Л. Идентификация систем: Теория для пользователя. М.: Мир, 1991.

42. Льюнг Леннарт. Идентификация систем: Теория для пользователя. Перевод с англ., под ред. Цыпкина Я.З. М.: Мир, 1991.

43. Мань Н.В. Оврагоперешаговый метод нелинейной оптимизации по биссектрисным направлениям // Тез. Док. 16-й Науч. Конф. Секция прикладной математики. Ханой : ХПИ, 1989

44. Мань Н.В. Оврагоперешаговый метод нелинейной оптимизации по перпендикулярным направлениям // Тез. Док. 16-й Науч. Конф. Секция прикладной математики. Ханой : ХПИ, 1989.

45. Мань Н.В. Поисковые методы оптимизации систем управления недетерминированными объектами (на примере теплоэнергетики). Авт. Реферат д.т.н. Дис. -М.: МЭИ, 1999.

46. Мань Н.В. Применение оврагоперешагового метода оптимизации для идентификации передаточной функции объекта управления // Теплоэнергетика, 1995, № 6, с. 71-77.

47. Мань Н.В. Принцип перешагивания в нелинейной оптимизации /Материалы Док. Науч. Семинара по Прикладной математике. Институт математики и ХПИ. Ханой, 1990

48. Мань Н.В., Аракелян Э.К., Хунг Н.Ч. Оптимизация фактического режима эксплуатации теплоэнергетических установок // Вестник МЭИ, 1988, №4, с. 56-61.

49. Мань Н.В., Ий Б.М. Метод биссектрисных направлений для решения задачи безусловной оптимизации // Журнал Математики, 1987, № 4, с. 113.

50. Метеорология и атомная энергия. — Л.: Гидрометеоиздат, 1971.

51. Методы оптимизации с приложениями к механике космического полета. Под ред. Дж. Лейтмана. М.: Изд. "Наука", 1965.

52. Нгуен Ван Мань. Поисковые методы оптимизации систем управления недетерминированными объектами. Автореф. докт. дисс. Москва, 1999.

53. Невчаев В.В. Совершенствование методики расчета показателей готовности оборудования электростанций // Электрические станции, 1976, №8, с: 4-7

54. ОНД — 86. Методика расчета концентрации в атмосферном воздухе вредных веществ, содержащихся в выбросах предприятий. — Л.: Гидрометеоиздат, 1986.

55. Опойцев В.И. Идентификация статических объектов кусочно-линейными функциями //Автоматика и телемеханика, 1970. № 5.

56. Перельман И.И., Поляков О.А. Идентификация объекта методами кусочной аппроксимации //Автоматика и телемеханика, 1968. № 10.

57. Пледнев Г.П. Автоматизированное управление объектами тепловых электростанций. -М.: Энергоатомиздат, 1981.

58. Пледнев Г.П. Автоматизированные системы управления объектами тепловых электростанций. М.: МЭИ, 1995.

59. Плетнев Г.П., Щедеркина Т.Е. Управление электрической нагрузкой энергоблоков ТЭС с учетом эксплуатационных ограничений //Изв. Вуз. Энергетика, 1983. № 5. С. 55-60

60. Плетнев Г.П., Щедеркина Т.Е., Горбачев A.C. Автоматическое управление вредными выбросами в переменных режимах ТЭС //Теплоэнергетика, 1995. № 4. С. 54-56.

61. Подиновский В.В., Ногин В .Д. Парето-оптимальные решения многокритериальных задач. М.: Наука, 1982.

62. Поляк Б.Т. Введение в оптимизацию. М.: Наука, 1983

63. Попырин Л.С. Математическое моделирование и оптимизация теплоэнергетических установок. М.: Энергия, 1978

64. Построение математических моделей химико-технологических объектов /Под ред. Е.Г. Дудникова. JL: Химия, 1970

65. Прокопенко А.Г. Стационарные, переменные и пусковые режимы энергоблоков ТЭЦ. М.: энергоатомиздат, 1990, с.315.

66. Райбман Н.С., Чадеев В.М. Построение моделей процессов производства. М.: Энергия, 1975.

67. Реклейтис Г.В., Рейвиндрэн А., Рэгсделл K.M. Оптимизация в технике. — М.: Мир, 1986.

68. Росляков П.В., Егорова JI.E., Ионкин И.Л. Методы расчета выбросов вредных веществ с дымовыми газами котлов М.: Изд. МЭИ, 2000.

69. Ротач В.Я. Теория автоматического управления теплоэнергетическими процессами. М.: Энергоатомиздат, 1985.

70. Рузинов Л.П. Статистические методы оптимизации химических процессов. М.: «Химия», 1972.

71. Сейдж Э.П., Мелса Д.Л. Идентификация систем управления. Перевод с англ., под ред. Райбмана Н.С. М.: Изд. "Наука", 1974.

72. Современные методы идентификации систем. Под ред. Эйкхоффа. Перевод с англ. под ред. Я.З. Цыпкина. М.: Мир, 1983.

73. Стерман Л.С. и др. Тепловые и атомные электрические станции. М.: Энергоатомиздат, 1995

74. Стратегия развития добычи угольного топлива в Вьетнаме до 2005 г. (на вьетнамском языке). Министерство промышленности Вьетнама, Ханой, 4/2000

75. Теверовский Е. Н., Артемова Н. Е., Бондарев А. А. и др. Допустимые выбросы радиоактивных и химических веществ в атмосферу. / Под ред. Теверовского Е. Н. и Терновского И. А — М.: Энергоатомиздат, 1985

76. Теория управления. Идентификация и оптимальное управление. Спиди К., и другие. Перевод с англ. М.: Изд. "Наука", 1973.

77. Тихонов А.Н. Решение некорректно поставленных задач и метод регуляризации//ДАН СССР, 1963. Т. 151. № 3. С. 184-187.

78. Фиакко A.B., Мак-Кормик Г.П. Нелинейное программирование. Методы последовательной безусловной оптимизации. М.: Мир, 1972

79. Химмельблау Д. Прикладное нелинейное программирование. Перевод с англ. под редакцией M.JI. Быховского. М.: Мир, 1975.

80. Хоа JI.K. и Мань Н.В. Оценка радиоактивных выбросов угольных ТЭС во Вьетнаме. Труды международной конференции "CONTROL 2000", Москва. 26-29 сентября 2000г.

81. Хуиг Н.Ч. Многоцелевая оптимизация режимов работы теплоэнергетических установок. Авт. Реферат. Канд. Дис. -М.:, 1997.

82. Черноруцкий И.Г. Оптимальный параметрический синтез. Электрические устройства. JL: Энергоатомиздат, 1987

83. Шитов Н.Ф., Гаряев А. Б. Определение поля концентрации примеси от промышленных источников при изменяющихся метеоусловиях. М.: Издательство МЭИ, 1999.

84. Эйкхофф П. Основы идентификации систем управления. — М.: Мир, 1975.

85. ATlas of Vietnam, Hanoi, February 2001.

86. Bellman RE. Dynamic programming. Princeton university press, 1957.

87. Master Plan of Power Development for 2000 2010. Ministry of Industry of Vietnam, Hanoi, 2/2000

88. Nguyen Van Manh, Bui Minh Tri. Method of "cleft-overstep" by perpendicular direction for solving the unconstrained non-linear optimization problem // ACTA Math.Vietnam, 1990, № 2, p. 73-83.

89. Nguyen Van Manh. The affine projection method for solving non-linear optimization problems //NCST Proceedings, 1992., № 2, P. 53-60.

90. Pasqill F. Atmospheric diffusion. Van Nostr. Co. Ltd. L., 1962.

91. Reklaitis G.V., Ravindran A., Ragsdell K.M. Engineering optimization. Method and application. New York: Wiley - Interscience, 1983.

92. The environmental issues of existing and planned power plants in Vietnam. VINATOM. Hanoi, July 1997.