автореферат диссертации по строительству, 05.23.16, диссертация на тему:Совершенствование методов прогноза гидродинамического и термического режимов водоемов-охладителей

ИИ4609502

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

Исаенков Александр Юрьевич

СОВЕРШЕНСТВОВАНИЕ МЕТОДОВ ПРОГНОЗА ГИДРОДИНАМИЧЕСКОГО И ТЕРМИЧЕСКОГО РЕЖИМОВ ВОДОЕМОВ-ОХЛАДИТЕЛЕЙ

Специальность 05.23.16 - Гидравлика и инженерная гидрология

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

3 0 СЕН 2010

Москва, 2010

004609502

Работа выполнена в лаборатории динамики русловых потоков и ледотермики Института водных проблем РАН (ГУ ИВП РАН)

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

Ведущая организация:

ГНУ ВНИИГиМ Россельхозакадемии

Защита состоится «18» октября 2010 г. в 16 ч. 30 мин. на заседании Диссертационного совета Д 220.045.02 во ФГОУ ВПО Московском государственном университете природообустройства (МГУП) по адресу: 127550, Москва, ул. Прянишникова, д. 19, аудитория 201. Факс: +7 (495) 976-10-46 (в Диссертационный совет Д 220.045.02 для Евдокимовой И.М.).

С диссертацией можно ознакомиться в библиотеке ФГОУ ВПО Московского государственного университета природообустройства (МГУП) по адресу: 127550, Москва, ул. Прянишникова, д. 19.

ведущий научный сотрудник

Елена Ивановна Дебольская

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

ГУ «ГГИ», лаборатория речных наводнений Бузия Владимир Александрович

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

ФГОУ ВПО МГУП, кафедра вычислительной

техники и прикладной математики

Верхоглядов Андрей Александрович

Автореферат разослан

Ученый секретарь Диссертационного совета, кандидат технических наук, доцент

И.М. Евдокимова

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

Актуальность исследования

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

Развитие численных методов позволило добиться существенных достижений как в теоретическом, так и в практическом направлениях этих исследований. Интенсивное развитие математического моделирования в области океанологии и гидрологии, происходящее в последнее Бремя, предоставляет все новые возможности для исследования природных процессов, влияющих на состояние внутренних водоемов, антропогенного воздействия на них и оценки возможности их использования. В то же время при всем многообразии существующих моделей встает проблема выбора наиболее подходящей из них для адекватного отражения процессов, происходящих в конкретном природном объекте, всех значимых параметров и действующих сил, формирующих его гидро- и термодинамическую структуру. Вышеизложенными обстоятельствами и определяется актуальность данной работы.

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

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

термических параметров не является равномерной и может вносить существенные отклонения в результаты расчетов, проведенных с использованием двухмерных моделей. Цель работы .

Цель настоящего исследования - модификация и адаптация глобальной математической модели для расчетов характеристик гидродинамического и термического режимов водоемов-охладителей при изменяющихся внешних воздействиях и батиметрии. Задачи исследования

Для достижения поставленной цели потребовалось решение ряда задач:

- выполнение сравнительного анализа существующих методов, основанных на экспериментально полученных данных, и численных методов моделирования гидродинамического и термического режимов водоемов-охладителей;

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

- адаптация трехмерной математической модели для условий конкретного водоема-охладителя;

- верификация модели но данным натурных измерений;

- проведение модельных расчетов;

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

Научная новизна

- Адаптация глобальной трехмерной модели для исследования особенностей гидродинамического и термического режимов водоемов-охладителей, связанных с неравномерностью вертикального распределения их параметров.

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

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

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

- Получены новые данные натурных исследований гидродинамического и термического режимов мелководных водоемов-охладителей.

- Выявлены новые закономерности гидродинамического режима водоемов-охладителей.

Практическая значимость работы

Практическая значимость исследования определяется возможностью применения полученных результатов для оценки рисков возникновения теплового загрязнения вод и нарушения экологического равновесия в водоемах-охладителях. Даны рекомендации эксплуатирующим организациям по наиболее оптимальному использованию охлаждающей способности водоемов-охладителей Шатурской ГРЭС, сохранению их экологической стабильности; оценке целесообразности строительства повых струеналравляющих сооружений и проведения дноуглубительных работ с вводом нового энергоагрегата ПГУ-400 при обеспечении режима эксплуатации водоемов-охладителей в соответствии с нормативными требованиями по охране поверхностных вод.

Достоверность полученных результатов

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

Результаты, полученные в рамках диссертационной работы, докладывались и обсуждались на заседаниях семинаров и Ученого Совета ИБП РАН, на VII конференции «Динамика и термика рек, водохранилищ и прибрежной зоны морей».

Публикации

Результаты, полученные в диссертации, опубликованы в пяти печатных работах, в том числе четыре в изданиях, включенных в Перечень ВАК. Структура и объем диссертации

Работа состоит из введения, пяти глав, заключения, списха использованной литературы, включающего 92 наименования. Полный объем диссертации составляет 127 страниц, включая 30 иллюстраций и 21 таблицу.

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

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

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

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

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

Изложена сущность и механизм переноса тепла в водоемах в летний период. Рассмотрены условия пространственной направленности переноса.

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

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

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

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

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

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

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

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

В разделе 1.2 рассматриваются водоемы-охладители. Указываются их характеристики, особенности их термического режима, приводится классификация.

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

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

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

Гидротермику водоема-охладителя определяют тепло- и массообмен от водной поверхности в атмосферу и процессы переноса тепла непосредственно в самой водной массе водоема. Между зтими процессами существует обратная связь, так как температура на поверхности определяет интенсивность охлаждения, которая в свою очередь влияет на температурное поле в водохранилище, условия стратификации,' развитие плотностных течений и турбулентный обмен.

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

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

водоема-охладителя соответственно; f - коэффициент трения на поверхности раздела слоев; в - циркуляционный расход; ц - коэффициент разбавления на водовыпуске; Ъ -коэффициент температурного расширения воды; ДГ - температурный перепад водоема-охладителя, принимаемый равным температурному перепаду на конденсаторах турбин; g -ускорение свободного падения.

На неглубоких водоемах-охладителях сброс подогретых вод с тепловых или атомных электростанций вносит искажение в нормальное распределение температур воды и скоростей течения. Подогретые воды летом растекаются по поверхности водоема, образуя в его верхних слоях термоклин мощностью в 3-5 м с наибольшей температурой вблизи поверхности воды. Вглубь водоема температура воды уменьшается с градиентом 1,5-2,0°С/м. К границам распространения подогретых вод по площади водохранилища температурные градиенты убывают до нормальных. Влияние теплых вод сказывается в зависимости от размеров сброса на более или менее значительном расстоянии.

где - толщина верхнего слоя при штиле; I, В, Н- длина, ширина и средняя глубина

(2)

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

В разделе 2.1 анализируются метода гидротермических расчетов водохранилищ и водоемов-охладителей, основанные на экспериментально полученных данных. Анализируются допущения, используемые в этих методах.

Рассмотрен метод суперпозиции, предложенный А.И. Пеховичем для решения гидротермических задач при сложной картине распределения температуры воды в водоеме. Метод охватывает многообразные условия, и применим для расчетов водохранилищ различных типов, условий наполнения, видов перемешивания, условий водной поверхности.

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

В разделе рассматриваются методы расчета водоемов-охладителей, в основе которых лежит метод теплового баланса. Анализируется глубокая методика гидротермического расчета водоемов-охладителей, водозаборных и водовыпускных сооружений, предложенная ВНИИГ им. Б.Е. Веденеева. Рассматриваются задачи, для решения которых осуществляется расчет температурного режима водоемов-охладителей (определение температуры охлажденной циркуляционной воды, необходимой площади водоема-охладителя, предельной мощности ГРЭС). Обозначаются необходимые условия для проведения расчетов, требуемые исходные данные. Рассматривается необходимость учета ряда факторов при расчете температуры водоема-охладителя по уравнению теплового баланса, включения ряда составляющих в уравнение теплового баланса при существенном изменении объема воды в водоеме-охладителе и соответствующем изменении теплосодержания.

В разделе 2.2 цроводится анализ методов численного моделирования при расчете водоемов-охладителей и особенностей их применения. Рассматриваются работы ряда научных организаций, занятых исследуемой проблемой (Институт вычислительного моделирования СО РАН, Институт гидродинамики СО РАН, ИВП РАН, НИИЭС, ВНИИГ им. Б.Е. Веденеева, ИФТПЭ АН Литвы).

Рассматривается применение одномерной модели, построенной на использовании уравнения энергии для описания температурного поля водоемов-охладителей и

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

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

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

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

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

В разделе 2.3 проводится обзор глобальных математических моделей, тенденция к разработке которых сложилась в последние годы за рубежом. Разработкой таких моделей занимаются крупные научные коллективы; модели широко распространены и апробированы на многих природных объектах. Приводится обоснование необходимости применения трехмерной модели для решения поставленных в работе задач. Требованием к такой модели является возможность работы с мелководными объектами и выделения вертикальной

структуры водоема. Отмечается модель ОЕТМ (лулууу.ееЬд.еи). как наиболее перспективная модель с точки зрения поставленных задач. Подчеркивается наличие у нее широких возможностей для расчетов и ряд других преимуществ.

В третьей главе представлено описание натурных исследований на водоемах-охладителях Шатурской ГРЭС. Дано описание объектов исследований, приводится методика проведения и основные результаты исследований.

В разделе 3.1 приводится обзор объектов исследования - Шатурской ГРЭС и водоемов-охладителей электростанции. Рассматриваются физико-географическое положение, климатические особенности района исследований, геологические, геоморфологические, гидрогеологические особенности, гидрологические характеристики водоемов-охладителей Шатурской ГРЭС (рис. 1). Указывается характеристика каждого из озер (Святого, Муромского, Белого, Черного) и системы охлаждения в целом. Приводится историческая справка и технические характеристики электростанции.

Рис. 1. Система водоемов-охладителей Шатурской ГРЭС.

В разделе 3.2 представлены характеристики измерительной аппаратуры (навигадиониой, гидрометрической, метеорологической) и методика проведения исследований.

В разделе 3.3 приводятся результаты натурных исследований на водоемах-охладителях Шатурской ГРЭС.

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

Определены морфометрические характеристики озер (размеры, конфигурация, площади зеркала); произведена батиметрическая съемка, в результате которой установлен характер распределения глубин в озерах, особенности донного рельефа, осложняющие циркуляцию вод в озерах, рассчитаны объемы озер при разных отметках уровня воды. Данные измерения позволили оценить Шатурские озера как мелководные.

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

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

Рассчитанные параметры водообмена позволяют предположить, что воды озер не могут иметь существенных отличий по своему составу, что полностью подтверждают результаты гидрохимического анализа. Выявлено, что по термическому режиму Шатурские озера близки к природным.

Определена зависимость температурный режима озер от метеорологических условий и теплового потока, создаваемого выпусками воды ГРЭС. Замечено, что характерной особенностью термического режима озер является частое превышение температуры воды

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

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

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

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

Приведена батиметрия водоемов-охладителей Шатурской ГРЭС (рис. 2), горизонтальное и вертикальное разрешение модели.

конвекции постепенно вовлекается в процесс испарения. Определена роль каждого из озер в

г,ш

40С

201

60С

Рис. 2. Батиметрия водоемов-охладителей Шатурской ГРЭС.

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

В пятой главе изложены результаты применения модифицированной для водоемов-охладителей модели йЕТМ и даны рекомендации по объектам исследования.

Отметается предстоящее в связи с планируемым вводом в эксплуатацию на Шатурской ГРЭС новой парогазовой установки (ПГУ-400) увеличение расхода поступающей нагретой воды, сокращение времени достижения водозабора тепловым потоком, повышение температуры воды в водоемах-охладителях и на водозаборе и возможные негативные последствия для работы станции и экологического состояния водоемов-охладителей.

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

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

Построено распределение поверхностпых температур воды при различных сценариях расположения струеналравляющих дамб и выемок в различные моменты времени (через 2, 5 и 10 дней) после начала водовыпуска с расходом <3 = 54 м3/'с и температурой = 33°С и экстремальных метеоусловиях: начальная температура воды в озерах ТВОды = 25°С, температура воздуха независимо от времени суток Тр.01„ = 25°С, скорость ветра V = 0 м/с. Анализ данных моделирования позволяет сделать заключение, что во всех случаях температуры на водозаборе различаются незначительно. В то же время различие температур в оз. Святом существенно, и самые низкие значения соответствуют сценарию при наличии выемок в оз. Святом. Отмечается, что этот же сценарий наиболее эффективен с точки зрепия нагревания воды на водозаборе.

2000 4000 X, М

Рис. 3. Батиметрия Шатурских водоемов-охладителей с различным расположением стрзгенаправляюБЗЯх дамб и наличием выемки: а - современное состояние; б-с удлиненной дамбой в оз. Святом и прорезью в существующей поперечной дамбе; в - с удлиненной дамбой в оз. Святом без поперечной дамбы; г - с существующим расположением дамб и наличием выемки в оз. Святом.

Построено распределение осредненных по глубине температур воды в более экстремальных условиях через 10 дней после начала водовыпуска с расходом 60 м'/с и температурой 30°С, начальная температура воды в озерах Тюль, = 22°С, температура воздуха независимо ог времени суток Твозд = 35°С, скорость ветра V = 0 м/с. Анализ графических результатов моделирования позволяет сделать вывод, что удлинение продольной дамбь1 будет иметь неблагоприятное воздействие на температурный режим в оз. Черном и особенно на нагревание оз. Святого. Установлено, что самым благоприятным сценарием в этих условиях является сценарий при отсутствии дамб.

Результаты расчетов для нормальных метеоусловий (температура воздуха 22°С) при современных нагрузках на ГРЭС (расход 20 м^с) позволили установить, что при нормальных условиях удлинение продольной дамбы приводит к незначительному снижению температуры воды на водозаборе в оз. Черном, но при этом нагревание воды в оз. Святом происходит гораздо быстрее, и распределение температуры по вертикали становится более однородным.

По результатам наблюдений на оз. Святом построено распределение по времени и глубине температур в геометрическом центре озера и в среднем по озеру для случаев с существующей и удлиненной дамбой (рис. 4 и 5). Из рисунков видно, что удлинение дамбы приводит к ускорению установления гомотермид с более высокими температурами.

Отмечается существенная роль ветрового воздействия и влажности воздуха в тепловом режиме Шатурских водоемов-охладителей. Наиболее неблагоприятным для работы ГРЭС является ветер северного направления. При моделировании рассчитаны варианты с разными направлениями и скоростями ветра, разной температурой точки росы и облачности для учета испарения. Анализ численных экспериментов показал, что учет этих параметров дает результаты, соответствующие наблюдениям. Установлено, что наибольшее влияние имеют ветровое воздействие и влажность, меньшее влияние оказывает облачность.

Рис. 4, Распределение температур: а - в вертикальном столбе в центре оз. Святого; б -в среднем по оз. Святому в случае существующей дамбы (с учетом суточного хода температур воздуха).

Рис. 5. Распределение температур: а-в вертикальном столбе в центре оз. Святого; б -в среднем по оз. Святому в случае удлиненной дамбы.

Построено распределение поверхностной температуры воды при существующем расположении дамб в экстремальных условиях через 5 (рис. 6, а, б) и 10 дней (рис. 6, в, г) после начала водовыпуска в отсутствие ветра (рис. 6, а, в) и при ветре северного направления со скоростью 5 м/с (рис. б, б, г). Замечено, что ветер значительно повышает температуру воды на водозаборе и, соответственно, на водовьшуске, но понижает температуру воды в оз. Святом. Анализ результатов моделирования показывает, что влияние новой продольной

дамбы усиливает воздействие северного ветра в сторону уменьшения температуры воды на водозаборе на 1°С, но увеличивает перегрев оз. Святого по сравнению с ситуацией, соответствующей современному расположению дамб, на 9°С (рис. 7).

Рис. б. Распределение поверхностной температуры воды в Шатурских озерах при существующем расположении дамб в экстремальных условиях через 5 дней {а, 6) и 10 дней (в, г) после начала водовыпуска в отсутствие ветра (а, в) я при ветре северного направления со скоростью 5 м/с (б, г) № = 54 м3/с, Тд= 33°С, Т80да = 25°С, Твозд = 25°С).

О 2000 4000 ][)М

т. °с

25.2 25.0 24.8 24.6

\

а I

1 3 5 7 I

т °с

32.0

24.0

/

%1 /-

1 / ь / /

/ } !

/ ( / б

, суш.

1 3 5 7 суш.

Рис. 7. Графики изменения температуры воды во времени в оз. Святом для различных сценариев расположения струенаправдяющих дамб: а - современное состояние дамб, в отсутствие ветра; б - с удлиненной дамбой в оз. Святом без поперечной дамбы, при ветре северного направления со скоростью 5 м/с.

Заключение

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

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

струи, связывающей систему водоемов-охладителей. Эти особенности могут бшъ существенны при моделировании теплового потока.

Измерение глубин и картографическое определение площади озер позволили рассчитать объемы озер и оценить их как мелководные.

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

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

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

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

4. Для верификации модели бЕТМ в применении к системе водоемов-охладителей Шатурской .ГРЭС проводились численные эксперименты с заданием в качестве начальных условий реальных метеорологических условий и известных параметров выпусков нагретой вода с ГРЭС, при которых проводились натурные измерения вертикальных распределений скоростей течений и температуры воды по акватории водоемов. Сравнение рассчитанных и

измеренных параметров продемонстрировало их хорошее соответствие и существование неоднородной вертикальной термической структуры.

5. Трехмерная численная модель ОЕТМ позволила описать наиболее важные термодинамические и гидродинамические процессы в исследуемом объекте применительно к различным условиям. Модель позволяет прогнозировать уровни воды, скорости, температуру. Особая ценность модели заключается в возможности ее использования для расчетов в условиях мелководья, учета стратификации различного происхождения, сложной топографии дна и берегов, задании сложного метеорологического и антропогенного форсинга.

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

6. Численные эксперименты с использованием трехмерной модели показали не только существование вертикальной неоднородности распределения температур воды и скоростей течении, но и значительное ее влияние на общий температурный режим водоемов-охладителей, и опровергли результаты моделирования с применением двухмерной модели.

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

Результаты трехмерного математического моделирования показали не только неэффективность строительства дополнительных струенаправляющих сооружений, но и ж негативное влияние на экологическое состояние водоемов-охладителей.

Публикации по теме диссертации

1. Дебслъския Е.И. Использование программы ОЕТМ для математического моделировании термического режима Шатурских озер-охладителей [Текст] / Е.И. Дебольская, О.Я. Масликова, А.Ю. Исаенков // Труды VII конференции «Динамика и термита рек, водохранилищ и прибрежной зоны морей». - РФФИ. ИБП РАН. -2009.-С. 168-176.

2. Грицук И.И. Модели прогнозирования транспорта наносов водными потоками на урбанизированных территориях в зимшш период [Текст] / И.И. Грицук, В.К. Дебольский, А.Ю. Исаенков, Н.К. Пономарев // Гидротехническое строительство. X» 1.-2010,-С. 38-43.

3. Дебольская Е.И. Исследование гидро- и термодинамических процессов в системе Шатурских озер-охладителей [Текст] / Е.И. Дебольская, О.Я. Масликова, А.Ю. Исаенков, И.С. Кузнецов // Гидротехническое строительство. - № 3. - 2010. - С. 2-10.

4. Дебольская ЕЛ Математическое моделирование деформаций русла в нижних бьефах ГЭС, расположенных в криолнтозоне, при катастрофических наводнениях в условиях ледовых затруднений [Текст] / Е.И. Дебольская, О Л. Масликова, А.Ю. Исаенков // Природообустройство.3.-2010. - С. 51-58.

5. Дебольская Е.И. Математическое моделирование транспортирующей способности подледных потоков [Текст] / Е.И. Дебольская, А.Ю. Исаенков //' Водные ресурсы. - Т. 37, № 5. -2010. С. 558-567. '

Подписано в печать 13.09.2010. Формат 60x84/16. Бумага типографская № 2. Печать - ризография. Усл. печ. л. 1,4 Тираж 100 экз. Заказ Л О ВЭ.

Московский государственный институт электроники и математики 109028, Москва, Б.Трехсвятительский пер., 3.

Центр оперативной полиграфии (495) 916-88-04, 916-89-25

Введение.

Глава 1. Современное состояние исследований.

1.1. Гидродинамические и термические процессы во внутренних водоемах в летний и зимний периоды.

1.1.1. Динамика течений в водоемах со свободной поверхностью.

1.1.2. Перенос тепла в водоемах в летний период.

1.1.3. Термические колебания в водоемах в летний период.

1.1.4. Замерзание водоемов.

1.1.5. Динамика течений и перенос тепла в водоемах, покрытых льдом.

1.2. Водоемы-охладители.

1.2.1. Показатели эффективности схемы использования водоема-охладителя

1.2.2. Типы водоемов-охладителей.

1.2.3. Особенности термического режима водоемов-охладителей.

1.2.4. Термическая структура водоемов-охладителей.

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

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

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

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

В работе рассматривается задача моделирования гидродинамических и термических процессов в водоемах-охладителях на основе натурных наблюдений и математического моделирования. Объектом исследования выбрана система водоемов-охладителей Шатурской ГРЭС. Для расчетов используется модифицированная для водоемов-охладителей Глобальная Эстуарийная Турбулентная Модель ОЕТМ (www.getm.eu). В модели используются различные возможности для расчета турбулентности (все известные модели замыкания), что позволяет досконально исследовать именно вертикальную структуру. Модель позволяет варьировать основные параметры, отвечающие за изменение термодинамического режима водоема, такие, как рельеф дна, температура, расход воды, поступающей при сбросах с электростанции и метеоусловия.

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

Научная новизна полученных результатов заключается в следующем:

- Адаптация глобальной трехмерной модели для исследования особенностей гидродинамического и термического режимов водоемов-охладителей, связанных с неравномерностью вертикального распределения их параметров.

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

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

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

- Получены новые данные натурных исследований гидродинамического и термического режимов мелководных водоемов-охладителей.

- Выявлены новые закономерности гидродинамического режима водоемов-охладителей.

Практическая значимость исследования определяется возможностью применения полученных результатов для оценки рисков возникновения теплового загрязнения вод и нарушения экологического равновесия в водоемах-охладителях. Даны рекомендации эксплуатирующим организациям по наиболее оптимальному использованию охлаждающей способности водоемов-охладителей Шатурской ГРЭС, сохранению их экологической стабильности; оценке целесообразности строительства новых струенаправляющих сооружений и проведения дноуглубительных работ с вводом нового энергоагрегата ПГУ-400 при обеспечении режима эксплуатации водоемов-охладителей в соответствии с нормативными требованиями по охране поверхностных вод.

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

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

В третьей главе представлено описание натурных исследований на водоемах-охладителях Шатурской ГРЭС. Дано описание объектов исследований, приводится методика проведения и основные результаты исследований.

Четвертая глава посвящена рассмотрению трехмерной модели, применяемой в настоящей работе для расчета изучаемых водоемов-охладителей. Приведены краткая история и перспективы развития модели, общая характеристика модулей, основные уравнения гидродинамики и переноса. Описывается применение модели в рамках диссертации.

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

В работе использован и проработан обширный материал, подготовленный на базе многочисленных теоретических разработок, натурных и лабораторных исследований динамики течений, процессов переноса тепла в водохранилищах и водоемах-охладителях, а также процессов ледообразования и развития ледовых явлений, проводившихся за последнее столетие как отечественными, так и зарубежными учеными. По теме исследования проведен анализ публикаций: монографий и коллективных работ, сборников статей и материалов конференций, содержащих как обобщение достижений науки, так и материалы исследований авторов, а также первоисточников этих работ и официальных данных научно-исследовательских и проектных организаций в сети Интернет. Используется материал обширных натурных и лабораторных исследований на водоемах-охладителях Шатурской ГРЭС, проводившихся в 2008-2009 гг. в ИВП РАН с участием автора. С перечнем использованных при написании работы источников и именами авторов можно ознакомиться в списке использованной литературы.

5.4. Общие выводы по главе и рекомендации

Численные эксперименты с использованием трехмерной модели показали не только существование вертикальной неоднородности распределения температур воды и скоростей течений, но и значительное ее влияние на общий температурный режим водоемов-охладителей.

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

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

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

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

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

Заключение

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

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

Измерение глубин и картографическое определение площади озер позволили рассчитать объемы озер и оценить их как мелководные.

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

Исходя из порядка величин расходов воды на водовыпусках и расходов воды в местах сужений, определены параметры водообмена водоемов-охладителей. Время водообмена увеличивается с уменьшением отношения объема озера к площади его зеркала. Уменьшение степени охлаждения вод связано с увеличением расходов воды на водовыпуске. Шатурские озера, исходя из показателей водообмена, нельзя отнести к природным бессточным озерам или озерам' со стоком. Тем не менее, эти озера' по термическому режиму близки к природным.

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

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

4. Для верификации модели вЕТМ в применении к системе водоемов-охладителей Шатурской ГРЭС проводились численные эксперименты с заданием в качестве начальных условий реальных метеорологических условий и известных параметров выпусков нагретой воды с ГРЭС, при которых проводились натурные измерения вертикальных распределений скоростей течений и температуры воды по акватории водоемов. Сравнение рассчитанных и измеренных параметров продемонстрировало их хорошее соответствие и существование неоднородной вертикальной термической структуры.

5. Трехмерная численная модель ОЕТМ позволила описать наиболее важные термодинамические и гидродинамические процессы в исследуемом объекте применительно к различным условиям. Модель позволяет прогнозировать уровни воды, скорости, температуру. Особая ценность модели заключается в возможности ее использования для расчетов в условиях мелководья, учета стратификации различного происхождения, сложной топографии дна и берегов, задании сложного метеорологического и антропогенного форсинга.

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

6. Численные эксперименты с использованием трехмерной модели показали не только существование вертикальной неоднородности распределения температур воды и скоростей течений, но и значительное ее влияние на общий температурный режим водоемов-охладителей, и опровергли результаты моделирования с применением двухмерной модели.

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

Результаты трехмерного математического моделирования показали не только неэффективность строительства дополнительных струенаправляющих сооружений, но и их негативное влияние на экологическое состояние водоемов-охладителей.

1. Вернадский, Н.М. Опыт теории и практики расчета прудов-холодильников Текст. / Н.М. Вернадский, Б.В. Проскуряков // Материалы по гидрологии, гидрографии и водным силам СССР. Вып. V. - М.-Л.: Госэнергоиздат. - 1931.

2. Васильев, О.Ф. Гидротермический режим водоемов-охладителей тепловых и атомных электростанций Текст. / О.Ф. Васильев, В.И. Квон, И.И. Макаров // Изв. АН СССР, сер. энерг. и трансп. -№ 4. 1976. - С. 102-110.

3. ВНИИГ им. Б.Е. Веденеева Электронный ресурс. // ОАО «ВНИИГ им. Б.Е. Веденеева»: официальный сайт. Режим доступа: http://vniig.ru/

4. Гавриш, П.Д. Водохранилища и водооградительные сооружения ГАЭС, ТЭС и АЭС Текст. / П.Д. Гавриш и др.; под ред. Т.П. Доценко // М.: Энергоатомиздат. 1989. -192 с.

5. Гончаров, В.Н. Динамика русловых потоков Текст. / В.Н. Гончаров // Л: Гидрометеоиздат. 1962.

6. Готлиб, Я.Л. Термика водохранилищ и гидротермические расчеты Текст. / Я.Л. Готлиб // Сборник лекций. Гидрология озер и водохранилищ. Часть 2. Водохранилища. М.: МГУ. - 1975. - С. 142-149.

7. Гришанин, К.В. Динамика русловых потоков Текст. / К.В. Гришанин // Л.: Гидрометеоиздат. 1969. - 428 с.

8. Губин, М.Ф. Гидроаккумулирующие электростанции (основные вопросы проектирования, исследования и строительства) Текст. / М.Ф. Губин // Гидроэнергетика. Т. 2. - М.: ВИНИТИ. - 1975.

9. Дебольский, В.К. Статистическое описание турбулентного движения в реках Текст. / В.К. Дебольский, E.H. Долгополова, O.A. Замай, A.C. Орлов // Вод. ресурсы. №4.- 1986.

10. Донченко, Р.В. Исследования и расчеты интенсивности шугообразования и зажорных явлений на зарегулированных участках рек Текст. / Р.В. Донченко, А.М. Филиппов, М.И. Баюсова, Н.С. Чачина // Тр. IV Всесоюз. гидрол. съезда. — Т. 6. J1.- 1976.-С. 351-360.

11. Донченко, Р.В. Ледовый режим водохранилищ СССР Текст. / Р.В. Донченко // Тр. ГГИ.-Вып. 187,- 1971.-С. 3-108.

12. Донченко, Р.В. Ледовый режим рек СССР Текст. / Р.В. Донченко // Л.: Гидрометеоиздат. 1987.

13. Дрижюс, М.-Р. Гидротермический режим водохранилищ-охладителей Текст. / М.-Р. Дрижюс // Вильнюс: Москлас. 1985. - 166 с.

14. Дрижюс, М.-Р. Численное моделирование теплового сброса в прямоугольный бассейн Текст. / М.-Р. Дрижюс, О.П. Кайрис // Тр. АН ЛитССР. Сер. Б. Т. 6 (139).- 1983.-С. 62-68.

15. ИВМ СО РАН Электронный ресурс. // Институт вычислительного моделирования СО РАН: официальный сайт. Режим доступа: http://icm.krasn.ru

16. Караушев, A.B. Расчет ветрового нагона в водохранилищах Текст. / A.B. Караушев // Тр. ГГИ. Вып. 35 (89). - 1952. - С. 94-123.

17. Караушев, A.B. Речная гидравлика Текст. / A.B. Караушев // JL: Гидрометеоиздат. -1969.-417 с.

18. Карты Google Электронный ресурс, Карты. // Google maps. Режим доступа: http://maps.google.com

19. Квон, В.И. Температурно-стратификационное течение в проточном водоеме Текст. /

20. B.И. Квон // Метеорология и гидрология. -№7. 1979.

21. Козлов, Д.В. Лед пресноводных водоемов и водотоков Текст. / Д.В. Козлов // М.: МГУП. 2000. - 263 с.

22. Козлов, Д.В. Основы гидрофизики Текст. / Д.В. Козлов // М.: МГУП. 2004. - 246 с.

23. Колесников, А.Г. Метод расчета температуры воды для периода осеннего охлаждения водохранилищ Текст. / А.Г. Колесников, A.A. Пивоваров // Вестн. МГУ. Сер. 3.-№2.- 1960.

24. Колесников, А.Г. Ход температуры воды в водохранилище в зимний период Текст. / А.Г. Колесников // ДАН СССР. Т. 92, № 1. - 1953.

25. Корытникова, H.H. О прогнозе зимних температур в непроточных водоемах Текст. / H.H. Корытникова // Изв. АН СССР, сер. геогр. и геофиз. № 6. - 1940.

26. Кочергин, В.П. Теория и методы расчета океанических течений Текст. / В.П. Кочергин // М.: Наука. 1978.

27. Крицкий, С.Н. Гидрологические основы речной гидротехники Текст. / С.Н. Крицкий, М.Ф. Менкель // Изд. АН СССР. 1950.

28. Крицкий, С.Н. Зимний термический режим водохранилищ, рек и каналов Текст. /

29. C.Н. Крицкий, М.Ф. Менкель, К.И. Россинский // M.-JL: Госэнергоиздат. 1947.

30. Ляхтер, В.М. Гидравлические исследования численными методами Текст. // В.М. Ляхтер, А.Н. Милетеев // Водные ресурсы. -№ 3. 1981. - С. 60-79.

31. Маккавеев, В.М. Теория процессов перемешивания при турбулентном движении свободных потоков и вопросы зимнего режима рек Текст. / В.М. Маккавеев // Зап. ГГИ.-Т. 5.-Л.: Гидрометеоиздат. 1931.-С. 75-115.

32. Марчук, Г.И. Математическое моделирование поверхностной турбулентности в океане Текст. / Г.И. Марчук, В.П. Кочергин, В.И. Климок, В.А. Сухоруков // Изв. АН СССР, сер. физ. атм. и океана. Т. 12. - № 8. - 1976.

33. Методические указания по технологическим расчетам водоемов-охладителей: РД 153-34.2-21.144-2003 Текст. Санкт-Петербург: ОАО «ВНИИГ им. Б.Е. Веденеева». - 2004.

34. Милетеев, А.Н. О расчете гидротермического режима нестратифицированных водоемов-охладителей Текст. / А.Н. Милетеев // Энергетическое строительство. № 3.- 1979.-С. 47-51.

35. Монин, A.C. Статистическая гидромеханика Текст. / A.C. Монин, A.M. Яглом // М.: Наука.-Т.1.- 1965.

36. О дрова, Т.В. Гидрофизика водоемов суши Текст. / Т. В. Одрова // Л.: Гидрометеоиздат. 1970.

37. Пехович, А.И. Расчет неустановившегося температурного режима в пластине с помощью графиков Текст. / А.И. Пехович, В.М. Жидких // Л.: Энергия. 1963.

38. Пехович, А.И. Расчеты теплового режима твердых тел Текст. / А.И. Пехович, В.М. Жидких // Л.: Энергия. 1963.

39. Пехович, А.И. Тепловой расчет глубоких водохранилищ в безледоставный период Текст. / А.И. Пехович // М.: Гидротехническое строительство. -№ 11. 1959.

40. Пивоваров, A.A. Предвычисление средних по вертикали подледных температур воды проточных водоемов Текст. / A.A. Пивоваров // Вестн. МГУ. Сер. 3. -№ 4. — 1966.

41. Пивоваров, A.A. Термика замерзающих водоемов Текст. / A.A. Пивоваров // Изд. МГУ.-1972.

42. Пиотрович, В.В. Об уточнении определений некоторых осенних ледовых явлений Текст. / В.В. Пивоваров // Метеорология и гидрология. № 9. - 1953.

43. Пиотрович, В.В. Расчет сроков ледостава и очищения от льда озер и водохранилищ Текст. / В.В. Пиотрович // JL: Гидрометеоиздат. 1955.

44. Рекомендации по термическому расчету водохранилищ: П 78-79 / ВНИИГ Текст. -Л.: ВНИИГ им. Б.Е. Веденеева. 1979.

45. Россинский, К.И. Движение донных наносов Текст. / К.И. Россинский // Тр. ГГИ. -Вып. 160.- 1968.

46. Россинский, К.И. Коэффициент теплопереноса в малопроточных водоемах в период летнего нагрева Текст. / К.И. Россинский, М.П. Кравцова // Тр. ГГИ. Вып. 162. -Л.: Гидрометеоиздат. - 1969.

47. Россинский, К.И. Температурный режим глубоких водоемов Текст. / К.И. Россинский // Тр. III Всесоюзн. гидролог, съезда. Т. IV. - Л.: Гидрометеоиздат. -1958.

48. Россинский, К.И. Термический режим водохранилищ Текст. / К.И. Россинский // М.: Наука.- 1974.- 240с.

49. Саркисян, A.C. Основы теории и расчета океанических течений Текст. / A.C. Саркисян // Л.: Гидрометеоиздат. 1966.

50. Саркисян, A.C. Численный анализ и прогноз морских течений Текст. / A.C. Саркисян // Л.: Гидрометеоиздат. 1977.

51. Соколов, A.C. Оценка температурного режима водозабора тепловых и атомных электростанций по данным моделирования гидротермического режима водоемов-охладителей Текст. / Соколов A.C. // Известия ВНИИГ им. Б.Е.Веденеева. Т.240. -2002.-С. 185-190.

52. Сперанская, A.A. Исследование турбулентности малопроточных водохранилищ Текст. / A.A. Сперанская // Тр. III Всесоюзн. гидрол. съезда. Т. 4. - Л.: Гидрометеоиздат. - 1959.

53. Сперанская, A.A. О влиянии устойчивости на характеристики вертикального турбулентного обмена в пресном водоеме Текст. / A.A. Сперанская // Метеорология и гидрология. № 8. - 1966.

54. Указания по термическому расчету водохранилищ: ВСН 46-71 / Минэнерго СССР Текст. Л.: Энергия, Ленинградское отд. - 1972.

55. Фидман, Б.А. Потоки с неоднородной плотностью (обзор) Текст. / Б.А. Фидман // Тр. ГГИ. — Вып. 162. Л.: Гидрометеоиздат.- 1969.

56. Филатов, H.H. Некоторые особенности турбулентного обмена в озерах Текст. / H.H. Филатов // Изменчивость гидрофизических полей в озерах. Л.: Наука. - 1978. - С. 89-116.

57. Филатова, Т.Н. Исследование течений в озерах и водохранилищах Текст. / Т.Н. Филатова // JL: Гидрометеоиздат. 1972. - 319 с.

58. Хинце, И.О. Турбулентность: ее механизм и теория Текст. / И.О. Хинце; пер.с англ. О.В. Яковлевского; под ред. Г.Н.Абрамовича // М.: Физматгиз 1963. - 680 с.

59. Шатурская ГРЭС Электронный ресурс. // ОАО «ОГК-4»: официальный сайт. -Режим доступа: http://ogk4.ni/?obHd2300&id=1687

60. Шатурская ГРЭС Электронный ресурс. // Свободная энциклопедия ВикипедиЯ. -Режим доступа:http://ru.wikipedia.org/wiki/%D0%A8%D0%B0%D 1 %82%D 1 %83%D 1 %80%D 1 %81 %D 0%BA%D0%B0%D 1 %8F %D0%93%D0%A0%D0%AD%D0%A 1

61. Шатурские озера: озера Московской области Электронный ресурс. // Все водоемы России: по материалам сайта Fisr.ru. Режим доступа: http://fisr.ru/articleold/a-67.html

62. Штокман, В.Б. Вертикальное распространение тепловых волн в море и косвенные методы определения коэффициента теплопроводности Текст. / В.Б. Штокман // Тр. Балт. НИРО. Вып. 3. - Калининград. - 1959.

63. Шуляковский, Л.Г. О механизме вскрытия рек, текущих на север Текст. / Л.Г. Шуляковский II Тр. Центр, инст. Прогнозов. Вып. 6 (33). - М., Гидрометеоиздат. -1947.

64. Шуляковский, Л.Г. Появление льда и начало ледостава на реках, озерах и водохранилищах. Расчеты для целей прогнозов Текст. / Л.Г Шуляковский // Л.: Гидрометеоиздат. 1960. - 216 с.

65. Baumert, H. Hysteresis of turbulent kinetic energy in nonrotational tidal flows: A model study Text. / H. Baumert, G. Radach // J. Geophys. Res. Vol. 97. - 1992. - P. 36693677.

66. Blumberg, A.F. A description of a coastal ocean circulation model, in Three dimensional ocean models Text. / A.F. Blumberg, G.L. Mellor; ed. by N. S. Heaps // Washington, D.C.: American Geophysical Union. 1987. - P. 1-16.

67. Bryan, K. A numerical model for the study of the world ocean Text. / K. Bryan // J. Computat. Phys. Vol. 4. - 1969. - P. 347-376.

68. Burchard, H. Applied turbulence modelling in marine waters Text. / H. Burchard // Lecture Notes in Earth Sciences. Vol. 100. - Springer - Berlin - Heidelberg - New York. -2002.

69. Burchard, H. Hybridisation between o and z coordinates for improving the internal pressure gradient calculation in marine models with steep bottom slopes Text. / H. Burchard, O. Petersen // Int. J. Numer. Meth. Fluids. Vol. 25 - 1997. - 1003-1023.

70. Burchard, H. On the performance of a mixed-layer model based on the k-e turbulence closure Text. / H. Burchard, H. Baumert // J. Geophys. Res. Vol. 100. - 1995. - P. 85238540.

71. Burchard, H. Three-dimensional modelling of estuarine turbidity maxima in a tidal estuary Text. / H. Burchard, K. Bolding, M. R. Villarreal 11 Ocean Dynamics. Vol. 54. - 2004. -P. 250-265.

72. Burchard, H. Turbulenzmodellierung mit Anwendungen auf thermische Deckschichten im Meer und Strömungen in Wattengebieten: Ph.D. thesis Text. / H. Burchard // Institut für Meereskunde, Universität Hamburg. Report 95/E/30. - GKSS Research Centre. - 1995.

73. Cox, M.D. A primitive equation, 3-dimensional model for the ocean Text. / M.D. Cox // University of Princeton. Geophysical Fluid Dynamics Laboratory. — Tech. Rep. 1. -Princeton, N. J. 1984. - 75 p.

74. Devik, O. Thermische und dynamische Bedingungen der Eisbildung im Wasserlaufen Text. / O. Devik // Geophys. Puhl. 1932.

75. Duwe, K. Modellierung der Brackwasserdynamik eines Tideästuars am Beispiel der Unterelbe: Ph.D. thesis Text. / К. Duwe; Universität Hamburg // Hydromod Publ. No. 1. -Hamburg: Wedel. 1988.

76. Fischer, H.B. Mixing in Inland and Coastal Waters Text. / H.B. Fischer, etc. // New York London - Toronto - Sydney - San Francisco: Academic Press. - 1979. - 483 p.

77. Graf, W.H. Hydraulic of sediment transport Text. / W.H. Graf, // N.-Y.: McGraw-Hill. -1971.-513 p.

78. General Estuarine Transport Model (Глобальная Эстуарийная Модель Переноса) Электронный ресурс. // GETM: официальный сайт. Режим доступа: http ://•www. getm. eu

79. Haidvogel, D.B. Numerical Ocean Circulation Modelling Text. / D.B. Haidvogel, A. Beckmann // Series on Environmental Science and Management. Vol. 2. - London: Imperial College Press. - 1999.

80. Jerlov, N.G. Optical oceanography Text. / N.G. Jerlov // Elsevier. 1968.

81. Jirka, G.H. Mathematical predictive models for cooling ponds and lakes. Part A: Model development and design consideration Text. / G.H. Jirka, etc. // Ralph M. Parsons Lab. for Water Res. and Hydrodyn. Rep. N 238. Cambridge: MIT. - 1978. - 109 p.

82. Jirka, G.H. Steady-state estimation of cooling pond performance Text. / G.H. Jirka, M. Watanabe // J. Hydraul. Div. Vol. 106, N HY6. - 1980. - P. 1116-1123.

83. Kantha, L.H. Numerical models of oceans and oceanic processes Text. / L.H. Kantha, C.A. Clayson // International Geophysics Series. Vol. 66. - Academic Press. - 2000.

84. Kondo, J. Air-sea bulk transfer coefficients in diabatic conditions Text. / J. Kondo // Bound. Layer Meteor. Vol. 9. - 1975. - P. 91-112.

85. Kvon, V.l. Hydrothermical simulation of turbulent temperature stratified flows in water bodies Text. / V.l. Kvon // Proc. 20th Congr. IAHR. Vol. 4. - M. - 1983. - P. 10-19.

86. Mathieu, P.-P. The role of topography in small well-mixed bays, with application to the lagoon of Mururoa Text. / P.-P. Mathieu, etc. // Cont. Shelf Res. Vol. 22. - 2001. - P. 1379-1395.

87. Paulson, C.A. Irradiance measurements in the upper ocean Text. / C.A. Paulson, J.J. Simpson // J. Phys. Oceanogr. Vol. 7. - 1977. - P. 952-956.

88. Pedlosky, J. Geophysical fluid mechanics Text. / J. Pedlosky // 2. ed. New York: Springer. - 1987.

89. Taylor, G.I. The dispersion of matter in turbulent flow through a pipe Text. / G.I. Taylor // Proc. Roy. Soc. Ser. A. Vol. 223. - 1954. - P. 446-468.

90. Zanke, U.C.E. Grundlagen der Sedimentbewegung Text. / U.C.E. Zanke // Berlin: Springer-Verlag. 1982. - S. 402.

91. Zukauskas, A. Warm-water spreading in cooling ponds of steam power plants Text. / A. Zukauskas // Heat transfer and Turbulent Buoyant Convection. Vol. 1. - Washington: Hemisphere. - 1977. - P. 139-157.