CC BY
30
Вычислительные технологии
Том 17, № 1, 2012
Информационно-вычислительная технология для сценарных оценок динамики и качества атмосферы*
Э.А. Пьянова1, Л. М. Фалейчик2 1 Институт, вычислительной математики и математической геофизики СО РАН, Новосибирск, Россия
2Институт природных ресурсов, экологии и криологии СО РАН, Чита, Россия e-mail: [email protected], [email protected]
Представлено описание вычислительных технологий, развиваемых в ИВМиМГ СО РАН и ИПРЭК СО РАН для решения природоохранных задач.
Ключевые слова: численное моделирование, ГНС-технологии, динамика атмосферы, перенос примеси, монотонные схемы.
Введение
При хозяйственной деятельности человека с интенсивным техногенным вмешательством в природную среду возникает целый ряд природоохранных задач. Один из важных вопросов, который необходимо изучить перед началом промышленного освоения, — влияние проектируемых пародохозяйствеппых объектов на качество и режим циркуляции атмосферы.
Оценить воздействие на атмосферу формирующихся при разработке полезных ископаемых глубоких карьеров и гор отвалов, выбросов промышленных предприятий, проектируемых ГЭС и водохранилищ, АЭС и ТЭЦ с не замерзающими зимой водоемами-охладителями можно с помощью математического моделирования на основе сценарного подхода. Этот способ предполагает проведение серии численных экспериментов с моделированием набора сценариев, описывающих микроклиматические режимы, характерные для исследуемых территорий.
В работе представлены математическая модель для описания динамики атмосферы над неоднородной орографией и алгоритм получения на ее основе численных оценок вероятного влияния проектируемых антропогенных объектов на качество и режимы циркуляции воздуха. Предлагаемый информационно-вычислительный подход к оценке динамики и качества атмосферы в областях со сложным рельефом включает построение алгоритмов реализации и комплекса программ на базе вариационного принципа с использованием дискретно-аналитических схем для конвективно-диффузионных операторов, обладающих свойствами монотонности, транспортивности и консервативности. Новым здесь является также применение геоинформационных технологий в составе разработанного комплекса моделей и программ. Это существенно упростило процессы по-
* Работа выполнена при финансовой поддержке РФФИ (гранты № 11-01-00187 и 10-06-00060а), Программ фундаментальных исследований № 4 Президиума РАН и № 3 Отделения математических наук РАН.
лучения входной информации о подстилающей поверхности исследуемых территорий, а также расширило возможности анализа результатов и подготовки иллюстративного материала.
1. Модель динамики атмосферы и переноса примеси
Для моделирования микроклимата использовалась мезомаештабная негидростатическая модель динамики атмосферы [1]:
дри др' дру др'
—--Ь шургш = — —--Ь 1ру + /лии, —--Ь сиур^и = — —--1ри +
дЬ дх дЬ ду
дрт др' ы дрд' ы ы
—--\- шургои = — ——Ь Арг/ + —--\- сиург/ и = — Ьрги + р—Ф + А^хг ,
дЬ дг дЬ ср
дрд' , дд , дри дрь дри)
— + Лурд и = -ри,— - рФ + А,д , ^ + ^ + = (1)
Здесь Ь — время; х, у, г — декартовы координаты, и = (п,у,'ш), и, V, — компоненты вектора скорости ветра в направлениях х, у, г соответственно; §',д',р' — отклонения потенциальной температуры, массовой доли водяного пара и атмосферного давления от их фоновых значений; р(г) — заданная функция плотности воздуха; Q — фоновое значение массовой доли водяного пара; Ьш — скрытая теплота конденсации; Ф — скорость образования жидкой фазы; ср — удельная теплоемкость воздуха при постоянном давлении; I — параметр Кориолиса; Б — параметр стратификации; А — параметр плавучести. Фоновые поля, описывающие крупномасштабные атмосферные движения, предполагаются известными,
В начальный момент времени значения полей метеоэлементов приняты равными своим фоновым значениям, т, е, отклонения равны нулю. На боковых границах ставились однородные условия Неймана, На верхней границе задавались условия выхода метеоэлементов на их фоновые значения. Влияние орографических и термических неоднородностей подстилающей поверхности учитывалось в краевых условиях третьего рода на нижнем расчетном уровне модели [2], который совпадает с верхней границей приземного слоя,
В модели также рассчитывается температура подстилающей поверхности из уравнения баланса тепла на границе раздела атмосфера — почва. При этом учитываются суточный ход и неравномерное поступление прямой солнечной радиации на склоны различных экспозиции и уклона. Все расчеты проводились в предположении, что характеристики подстилающей поверхности являются кусочно-постоянными функциями в ячейках расчетной сетки, построенной на исследуемой области для численной реализации модели. Температура поверхности в ячейках, покрытых водой, в течение всего модельного времени (порядка двух суток) полагалась постоянной и задавалась в соответствии со средними значениями для рассматриваемых сезонов (использовались данные из справочной литературы). Наличие снежного покрова на исследуемых территориях учитывалось при задании альбедо подстилающей поверхности. Параметризация потоков тепла, влаги и импульса в приземном слое производится отдельно для каждой ячейки расчетной области с учетом характеристик поверхности и стратификации над ней.
Для моделирования пространственно-временной изменчивости пассивной примеси использовалось уравнение переноса с турбулентной диффузией
дс
— + иётс1с = Асс + Цх,г) (2)
с начальными и граничными условиями
дс дс
с(М) = с°(х), ¿ = 0, — = 0, х = 0,Х, — = 0, у = 0,У,
дх ду
дс дс
— = 0, х = Н, ис— = (Д. - и)с)с + Мх,у,г), г = 6(х,у),
где с — концентрация примеси; с0 — поле концентрации примеси в начальный момент времени; х = (х,у,г); и>с — скорость оседания прпмесн; /Зс — величина, характеризующая взаимодействие примеси с подстилающей поверхностью; f, ^ — функции источников,
гл л/ п\ . д д д д д д
Оператор Д„ (а = и,ь,т,{},д,с) имеет вид Да = —/¿а— + +
дх дх ду ду дх дх
где — коэффициенты турбулентной диффузии в горизонтальном и вертикальном
направлениях соответственно.
Система уравнений (1)-(2) с начальными и граничными условиями интегрировалась в области Д = Б х [0,Т], где Б ={0 < х < X, 0 < у < У, 8(х,у) + к < х < Н} — область изменения пространственных переменных, 6(х,у) — функция рельефа, к — высота приземного слоя, 0 < £ < Т — интервал времени.
2. Численная реализация модели
В данном разделе, следуя идеям и методике работ [3-5], приводится краткое описание методологии реализации модели гидротермодинамики атмосферы и переноса примеси над сложным рельефом. Численная реализация дифференциальных уравнений (1)-(2) с начальными и граничными условиями осуществлялась таким образом, чтобы разностные уравнения удовлетворяли балансным соотношениям, аналогичным тем, которые имеют место для исходной задачи в дифференциальной постановке, В качестве критерия для получения подходящих схем аппроксимации были выбраны соотношения баланса энергии и массы. Построение энергетически сбалансированных конечно-раз-ноетных уравнений для задачи (I) (2) проводилось на базе интегрального тождества [3]. Для этого системе уравнений гидротермодинамики атмосферы и переноса примеси сопоставлялся вариационный функционал (интегральное тождество) специальной структуры, который получается в результате домножения уравнений системы на весовые функции п*,у*,'ш*,'д*,д*,р*,с* и последующего интегрирования по частям по всей области Бг с учетом начальных и краевых условий.
Затем для полученного интегрального тождества I(ф, ф*) = 0, вде ф = (и, V, и>, д', р',с), ф* = (и*, V*, V*, §*, д*,р*, с*), строились дискретные аппроксимации, и из условия стационарности сумматорного функционала Iт(фт, ф*т) = 0 относительно произвольных и независимых вариаций сеточных функций ф*т получалась система конечно-разностных уравнений д1т (фт, ф*т)/дф*т = 0 (индекс т означает дискретизацию по пространству и времени).
При специальном выборе аппроксимаций вариационного функционала по времени [3] были получены схемы расщепления. Совместное использование интегрального
тождества и метода расщепления позволило упростить реализацию модели за счет сведения ее к последовательности более простых задач и обеспечить взаимное согласование всех этапов, на которые разбивалась исходная задача гидротермодинамики атмосферы и переноса примеси. Для аппроксимации подынтегральных выражений по времени использовались неявные схемы,
В соответствии с принципами расщепления по физическим процессам [4] решение системы уравнений (1)-(2) с начальными и краевыми условиями было разбито на решения трех задач: конвекции-диффузии, согласования метеополей и учета фазовых переходов влаги.
Поскольку концентрация примесей по физическому смыслу не может быть отрицательной, большое значение в реализации численных моделей, основанных на уравнении переноса с турбулентной диффузией, имеют конечно-разностные схемы, обладающие свойствами монотонности и транспортивности. Построение разностных схем для конвективно-диффузионных операторов первого этапа расщепления с указанными свойствами осуществлялось по методу, предложенному В,В, Пененко в [5]. При этом уравнения дополнительно расщеплялись по координатным направлениям, что позволило в дальнейшем решить одномерные сеточные уравнения экономичным методом прогонки.
Таким образом, для конвективно-диффузионных уравнений по каждому пространственному направлению были получены разностные схемы следующего вида:
(Д7+1 + - (и8А8В8)у+1(<р7+1 - <ру) + [и3А3)1_1^1 - = О,
(и3А3В3) > 0, (usAs) > 0, Я > 0, Us = u, v, т, (и3А3В3)1+1 = - (¡13['П(1 - сЬК ??)]/Дз)7+1 , (и3А3)1_1 = (¡13['П(1 + сЬК ??)]Дз)7_1 , Д7±1 = (Дв[1 Т ('ПсИг'п - 1)/??]/2)7±1 , = (иаАз/(2^))7±1,
Ав7 = (з7+1 — в1-1)/2, в = х, у, г, 7 = I, т, к. (3)
Положительность всех коэффициентов и диагональное преобладание обеспечивают монотонность и устойчивость разностных схем (3), Полученные численные схемы также являются балансными и транепортивными.
На этапе согласования полей метеоэлементов компоненты вектора скорости ветра выражались через отклонения давления и функции с предыдущего шага расщепления. После подстановки этих выражений в разностное уравнение неразрывности получался дискретный аналог уравнения Пуассона для нахождения отклонений давления, который реализовывалея итерационным методом сопряженных градиентов. Через найденные значения отклонений давления вычислялись компоненты скорости ветра, а затем и остальные искомые метеовеличины.
На третьем этапе расщепления в каждом узле расчетной сетки производился учет фазовых переходов влаги.
Построенные разностные схемы для уравнений динамики атмосферы и переноса примеси (1)-(2) имеют первый порядок аппроксимации по времени и второй по пространственным переменным.
Учет орографии в задаче осуществлялся с использованием идей метода фиктивных областей [1, 6, 7],
3. Подготовка данных о подстилающей поверхности
В численной модели рельеф местности и другие характеристики подстилающей поверхности исследуемой области задаются с помощью числовых значений в узлах регулярной прямоугольной расчетной сетки, В настоящее время информацию о высоте, крутизне и ориентации склонов в каждой точке модельной территории можно получить из разных источников. Возможности существующих ГИС-технологий позволяют использовать для этого и различные цифровые модели рельефа (ЦМР),
В качестве источника информации о рельефе были использовали данные SRTM (ftp://e0srp01u.ecs.nasa.gov), имеющие удовлетворительное для применяемой численной модели гидродинамики атмосферы разрешение (размер пиксела 30 х 30 м, погрешность измерения высоты не более 20 м).
Для автоматизации процесса извлечения из ЦМР информации о рельефе создан программный модуль, реализованный в ArcGIS 91 (Esri, Inc., USA) в виде расширения [8, 9], С помощью инструментов этого модуля строится регулярная расчетная сетка, для каждого узла которой из ЦМР извлекается и присваивается значение высоты соответствующей точки местности. Построенная расчетная сетка имеет географическую привязку, что дает возможность в дальнейшем результаты моделирования отобразить на топографическую основу и проанализировать в среде ГИС [8, 9]. Кроме того, из дополнительных тематических слоев, описывающих объекты гидрографии, транспортной сети, а также населенные пункты, технические сооружения, охраняемые территории и т.д., с использованием функциональности ГИС выделяются соответствующие им узлы сетки. Подготовленная таким образом информация о подстилающей поверхности экспортируется в текстовые файлы, применяемые затем в численной модели.
Таким образом, разработанный программный модуль позволяет извлекать из растровых слоев информацию о рельефе и преобразовывать ее в удобный для использования в численной модели вид и, тем самым, облегчает решение проблемы получения и подготовки информации о подстилающей поверхности при математическом моделировании природных процессов,
4. Организация численных экспериментов
Оценка воздействия, оказываемого промышленными объектами на атмосферу, выполнялась на основе метода прямого моделирования путем проведения двух серий сценарных расчетов. Первая серия описывает микроклимат исследуемой территории до антропогенного вмешательства в природную среду, вторая — отражает ситуацию после вмешательства. Различия в распределении полей метеоэлементов в соответствующих модельных сценариях каждой серии расчетов позволяют оценить зону и степень возможного влияния проектируемых объектов.
Режимы локальных циркуляции атмосферы, характерные для каждой исследуемой области, были получены с помощью мезомасштабной модели динамики атмосферы (1), Рассчитывался суточный ход в распределении полей основных метеоэлементов для различных сезонов года, т, е, сначала на входе в модель задавались распределения фоновых полей метеоэлементов, типичные для соответствующего сезона и рассматриваемой территории, затем на фоне этих климатических данных моделировались суточные колеба-
1 Программный продукт предоставлен Esri по гранту Conservation Program.
ния в распределении метеоэлементов вследствие термической неоднородности подстилающей поверхности. Неравномерный прогрев поверхности обусловливался наличием водных объектов на исследуемых территориях и неодинаковой инсоляцией склонов различных ориентации и наклона. Полученные в ходе моделирования квазистационарные поля метеоэлементов для ряда типичных ситуаций по каждому сезону года описывают микроклимат рассматриваемых территорий. Если в ходе исследования изучался вопрос о распространении примесей, то вместе с моделью гидротермодинамики использовалась также модель переноса примеси (2),
Рассчитанные таким образом поля метеоэлементов для ситуации без учета антропогенного объекта и для варианта, когда во входных данных модели учтены изменения термодинамических свойств подстилающей поверхности, вызванные антропогенным вмешательством, использовались для получения оценок влияния промышленных объектов на циркуляцию и качество атмосферы прилегающих территорий,
5. Использование ГИС-технологий для анализа результатов расчетов
Анализ результатов численных экспериментов и оценка размеров зоны влияния при-родно-технических объектов на окружающую территорию проводится в среде ГИС, Отображение результатов моделирования на топографической карте территории позволяет локализовать зоны существенного влияния промышленных объектов, а также установить причины некоторых особенностей распределения метеополей.
Результаты моделирования, сформированные в специальным образом организованный текстовый файл, импортируются в среду ГИС, где осуществляются их пространственная привязка и дальнейшее преобразование в тематические векторные и растровые слои, С использованием методов интерполяции строятся поверхности распределений метеоэлементов для каждого сценария и поверхности разностей соответствующих метеополей из параллельных сценарных расчетов [8]. Разности полей вычисляются для упрощения процесса локализации областей изменения метеоэлементов.
Кроме выявления зон значимых изменений в распределении каждого метеоэлемента в отдельности с помощью Калькулятора растров (Raster Calculator), Алгебры карт (Map Algebra) или оверлейных процедур из созданных ранее поверхностей разностей метеополей строятся различные интегральные показатели, учитывающие деформации нескольких метеополей одновременно. На результирующих поверхностях определяются области, в которых будет наблюдаться существенное изменение метеоэлементов в комплексе.
Так как единицы измерения и масштабы значений у всех метеоэлементов разные, то их необходимо нормировать с использованием одной из существующих метрик. Результатом нормировки являются новые растровые поверхности, все ячейки которых получают значения, определяемые, например, формулой
_ i fa ~ А«1 i та-\
где разным значениям индекса а соответствуют разные метеоэлементы Да1 — число, отклонения от которого в плане решения конкретных задач представляют наибольший интерес, Да2 _ величина разброса значений метеоэлемента <а. Можно применить
нормировку по дисперсиям (Д^) и математическим ожиданиям (Д^) либо по наилучшим или наихудшим значениям (Да1), Да2 = 'amax — 'aminj либо любую другую из существующих [10].
В нашем случае нормировка осуществлялась по формуле
ma = | 'a/'a max|, (4)
' a ' a max
симальпые значения, а итоговый суммарный показатель изменений комплекса метеоэлементов
n
ma
m = —--, (5)
n
где n — количество суммируемых полей.
Построенные в ГНС тематические векторные и растровые слои с их легендами, представляющие полученные результаты в виде поверхностей распределения основных моделируемых метеоэлементов на карте территории, организовываются и хранятся в геоинформационном каталоге проекта. Все слои имеют общую систему координат. Кроме того, в каталоге хранится используемая в численной модели информация о моделируемой территории, включая исходные данные о подстилающей поверхности, результаты модельных расчетов, а также слои топооеновы и другую информацию по данной территории, играющую в анализе вспомогательную роль, В этом же каталоге хранятся все построенные карты. Организация такого каталога дает возможность быстрого просмотра всех имеющихся данных по проекту, способствует оперативной подготовке любой уточняющей информации по нему, например, по отдельному участку модельной территории, дает возможность изменить способ отображения результатов, по-другому скомпоновать слои, участвующие в анализе, ввести дополнительные сведения о территории и т.д.
6. Примеры расчетов по предложенной методике
Рассмотрим примеры использования методики для оценки возможного изменения температуры, влажности и скорости ветра в зимнее и осеннее время в районе проектируемой Мотыгинекой ГЭС на Ангаре, Для проведения исследования была выбрана область моделирования размером 190 х 180 х 2.7 км, охватывающая территорию вокруг предполагаемой зоны затопления, а также часть территории со стороны нижнего бьефа плотины ГЭС,
Шаги горизонтальной регулярной сетки Дх = Ду = 1000 м, вертикальной — первые 15 шагов по 60 м, далее по 100, 150 и 200 м. Всего по вертикальной координате 29 узлов. Базовый шаг по времени ДЬ = 60 с. Для зимнего сценария расчетов стратификация фоновой атмосферы предполагалась сильно устойчивой. Фоновый ветер равен нулю. Температура воды в полынье 274 К, температура воздуха на поверхности 253 К, В численных экспериментах длина полыньи принималась равной 20 км, В расчетах осеннего сценария стратификация фоновой атмосферы полагалась устойчивой. Фоновый ветер равен нулю. Температура воды 280 К,
На рис, 1, о представлен фрагмент исследуемой территории с нанесенными изолиниями разности температур двух зимних сценарных расчетов (при наличии и отсутствии
Рис. 1. Изолинии разности температур, град (а), и массовой доли водяного пара, г/кг (б), на высоте 60 м над поверхностью; зимний сценарий
полыньи), иллюстрирующими зону отепляющего влияния полыньи. Так как температура воды в полынье выше температуры окружающей поверхности, то в прилегающей к ней зоне формируется локальная циркуляция по типу "острова тепла". Такая циркуляция приводит к тому, что на высоте приземного слоя над центром полыньи наблюдается повышение температуры атмосферы, достигающее в численных экспериментах 2-3 град. Повышение температуры воздуха на 1-2 град по результатам расчетов наблюдается на расстоянии порядка 1 км от полыньи.
В проведенных численных экспериментах для зимнего периода с учетом влияния незамерзающей полыньи Мотыгинской ГЭС с предполагаемой длиной 20 км повышение массовой доли водяного пара на 0.1 г/кг распространяется на 8-10 км от границ полыньи, при этом максимальное увеличение данного показателя на высоте приземного слоя наблюдается непосредственно над полыньей и составляет приблизительно 0.2 г/кг
100 Kilometers
Суммарный показатель Буферные зоны вокруг
---0.07 зоны затопления
-0.2 I I границы зон
— 0.5 ЕЖ зона затопления
Рис. 2. Суммарный показатель изменений метеорологических полей (температуры, влажности и скорости ветра) на высоте 60 м над поверхностью; осенний сценарий
(рис. 1,6). В целом над полыньей и в прилегающей к ней зоне рост относительной влажности с учетом низких зимних температур может достигать 15-17%.
На рис. 2 выделены зоны прогнозируемого комплексного изменения полей температуры, влажности и скорости ветра на высоте приземного слоя, полученное по расчетам для осеннего модельного сценария. Суммарный показатель изменений метеорологических полей вычислен по формулам (4), (5). Чем больше значение показателя, варьируемое от 0 до 1, тем существеннее ожидаемое комплексное изменение состояния атмосферы.
Анализ результатов численных экспериментов по оценке влияния на гидротермодинамический режим атмосферы прилегающих территорий водохранилища проектируемой Мотыгинской ГЭС позволил дать следующий прогноз: для территорий, удаленных от береговой линии на расстояние 10 км и более, возможны незначительные изменения микроклимата; в ближней зоне влияние водохранилища может быть весьма существенным. В зимний период это прежде всего касается районов, которые непосредственно прилегают к нижнему бьефу плотины и к незамерзающей полынье.
Заключение
Представленный программный комплекс и разработанная на его основе информационно-вычислительная технология для оценки состояния воздушного бассейна позволяют выполнять работы как по прогнозу и оценке возможных микроклиматических изме-
нений и атмосферных загрязнений территории в рамках экологических экспертиз технических проектов (водохранилища ГЭС, ТЭК, ТЭЦ, ГОК и др.), так и по оценке влияния уже существующих промышленных комплексов на качество геоэкологической среды прилегающей к объекту территории. Использование предложенной математической модели может оказать существенную помощь в изучении механизмов распространения загрязняющих веществ в сложных горных условиях и при решении других природоохранных задачах.
Предложенный подход был апробирован на конкретных объектах (Богучанекая и Мотыгинекая ГЭС на Ангаре) при подготовке информации для экспертных оценок [9, 11-13], результаты которых были использованы также для составления оценки воздействия на окружающую среду (ОВОС),
Авторы выражают благодарность В,В, Пененко и A.A. Фалейчику за полезные советы и обсуждения в ходе работы.
Список литературы
[1] Пененко В.В., Алоян А.Е. Модели и методы для задач охраны окружающей среды. Новосибирск: Наука, 1985. 256 с.
[2] Казаков А.Л., Лазриев Г.Л. О параметризации приземного слоя атмосферы и деятельного слоя почвы // Физика атмосферы и океана. 1978. Т. 4, № 3. С. 257-265.
[3] Пененко В.В. Методы численного моделирования атмосферных процессов. Л.: Гидро-метеоиздат, 1981. 352 с.
[4] Марчук Г.И. Методы вычислительной математики. М.: Наука, 1977. 456 с.
[5] Пененко В.В. Численные схемы для адвективно-диффузионных уравнений с использованием локальных сопряженных задач. Новосибирск, 1993 (Препр. ВЦ СО РАН; № 984).
[6] Алоян А.Е., Фалейчик A.A., Фалейчик Л.М. Алгоритм численного решения ме-зометеорологических задач в случае криволинейной области // Математические модели рационального природопользования. Новосибирск: Наука, 1989. С. 14-35.
[7] boyarshinova Е.А. Model of dynamics of atmosphere with monotone numerical schemes // Bull. Nov. Comput. Center. Num. Model, in Atmosph., etc., 2000. Iss. 6. P. 1-8.
[8] Фалейчик Л.М., Фалейчик A.A., Пъянова Э.А., Малюженко A.A. Использование ГНС-технологий для оценки зоны влияния искусственного водоема на микроклимат региона // VII Всерос. научно-практ. конф. "Кулагинские чтения". Чита: ЧитГУ, 2007. Ч. V. С. 48-53.
[9] Фалейчик Л.М., Пъянова Э.А. Использование ГНС-технологий и математического моделирования для оценки изменений состояния природной среды под влиянием хозяйственной деятельности человека // Вестник ЧитГУ. 2008. № 5(50). С. 117-127.
[10] Геоинформатика / Под ред. B.C. Тикунова. М.: Академия, 2005. 480 с.
[11] Фалейчик A.A., Фалейчик Л.М., Пъянова Э.А. Исследование влияния искусственного водоема на микроклимат прилегающих территорий с помощью математического моделирования // VII Всерос. научно-практ. конф. "Кулагинские чтения". Чита: ЧитГУ, 2007. Ч. V. С. 133-137.
[12] Пьянова Э.А., Фалейчик Л.М. Исследование влияния искусственного водоема на температурные режимы прибрежных территорий // Современные проблемы математического моделирования. Математическое моделирование и современные информационные технологии. Вып. 5. Ростов-на-Дону: Изд-во Южного федерального ун-та, 2007. С. 239-247.
[13] Pyanova Е.A., Faleychik L.M. Probable temperature regime changes of the atmosphere in Lower Angara region due to anticipated reservoirs // Intern. Conf. on Comput. Informat. Technologies for Environmental Sci. CITES 2009. Program^Abstracts. Russia, Krasnoyarsk, 2009. P. 35.
Поступила в редакцию 6 июня 2011 г., с доработки — 4 августа 2011 г.
