CC BY
28
УДК 502.3:504.61
1)01: 10.21209/2227-9245-2017-23-12-41-50
ИСПОЛЬЗОВАНИЕ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ ДЛЯ ОЦЕНКИ ЭКОЛОГИЧЕСКИХ ПОСЛЕДСТВИЙ РЕИНДУСТРИАЛИЗАЦИИ ЭКОНОМИКИ РЕГИОНА
USE OF MATHEMATICAL MODELING FOR ESTIMATION OF ENVIRONMENTAL CONSEQUENCES OF THE REGIONAL ECONOMY RE INDUSTRIALIZATION
Э. А. Цьянова, Институт вычислительной
математики и математической геофизики Сибирского отделения Российской академии наук, г. Новосибирск pUmova@t igs. тг i
Е. Pyanova, Institute of Computational Mathematics and Mathematical Geophysics, Russian Academy of Sciences, Siberian Branch, Novosibirsk
A. /1. Фалейчик, 3 a 6 a it к ал ьский государственный университет, г. Чита faa55@bk. ru
A. Faleychik,
Transbaikal State, i University, Chita
JI. M. Фалейчик, I Iitcinumym. природных ресурсов, экологии и криологии Сибирского отделения Российской академии
itat/к, г. Чита; Забайкальский государственный университет, г. Чита lfaleychik@bk, ru
L. Fateychik, Institute of Natural Resources, Ecology and Cryology; Russian Academy of Sciences, Siberian Branch, Chita: Transbaikal State University, Chita
Проведено исследование формирования атмосферных циркуляции и процессов распространения загрязняющих примесей к Читино-Иигодииской котловине в условиях сложной орографии. Рассмотрены особенности формирования гидрометеорологического режима атмосферы, важные сточки зрения опасности возникновения высоких уровней загрязнения воздушного бассейна города Читы. Разработаны сценарии возможного развития локальных атмосферных циркуляций в зимний и летний периоды. В зимнем сценарии моделируются процессы распределения загрязняющих выбросов от трех высотных источников в отсутствие фонового ветра. Летние сценарии моделируют процессы распространения дымовых шлейфов от одиночного очага лесного пожара в окрестностях города Читы, 11редставлены результаты сценарных расчетов, вы полненных на основе мезомасштабной негидростатической модели динамики атмосферы н переноса примеси в областях со сложным рельефом. Воспроизведены особенности развития атмосферных циркуляции над территорией в разных условиях. Показано, что климатические особенности и характеристики подстилающей поверхности существенно влияют на характер локальных циркуляции и на перераспределение атмосферных загрязнений над территорией города. Отмечено, что особую актуальность исследования в этой области приобретают в условиях происходящих процессов реин дуетри а л изации экономики региона, 1! частности при переходе предприятий на новые технологии и новые виды топлива
Ключевые слова: экономическое развитие региона: оценка антропогенного воздействия; загрязнение атмосферы; математическое моделирование; геогрифические ин^юрмационные системы; геоипформациошше технологии; реиндустриализация
© Э. А. Пьянова,Л. А. Фалейчик, Л. М. Фалейчик, 2017
The atmospheric circulation generation ami the pollutants transport processes in the Chitino-Ingodinsky basin under complex orography conditions arc studied. The orographic and climate features of the hydro meteorological regime formation of Ihe atmosphere air; important factors to creation of the pollution high levels iri the Chi la city air basin. The authors have created scenarios of possible development of local atmospheric circulations during Ihe winter and summer periods. The simulation on the winter scenarios was supposed that Ihe background wind was absent. Summer scenarios simulate the atmospheric transport of smoke from a forest fire in the Chita neighborhood. The results of scenario calculations which have been performed on the basis of a numerical model of atmospheric dynamics and impurity transport are presented. It was noted that the research in this field have a special significance in the context of ongoing processes of re industrialization of the region's economy, in particular, when businesses switch to new technologies and new fuel
Key words: economic development of region; anthropogenic impact assessment; air pollution; mull) modeling; geographical information systems; geoinformation technologies; reinduslrializatton
ведение. Процессы ренндустриализа-ции экономики России носят крайне неоднородный характер. Забайкалье не относится к числу лидеров экономического развития среди сибирских регионов. Тем не менее, данные процессы неизбежно набирают обороты, к связи с чем проблема экологической оценки результатов хозяйственной деятельности человека для всех природных сред, и частности для атмосферного воздуха и крупных городах, по-прежнему является актуальной. Рациональнее заранее тщательно проанализировать возможные неблагоприятные последствия тех или иных вариантов принимаемых решений, нежели впоследствии тратить ресурсы на исправление ошибок. Зачастую результаты неправильных решений имеют необратимый характер. 13 крупных городах, особенно Сибирского региона, таких примеров достаточно много.
Наиболее распространенными причинами низкого качества городского воздуха являются большие объемы загрязняющих выбросов от промышленных предприятий, автотранспорта и других источников-за-грязиителей. Важно отметить, что качество воздуха также существенно зависит от погодных условий. I Гримерами неблагоприятных метеорологических факторов являются безветренная погода, инверсионные условия, неблагоприятные направления ветра, способствующие переносу загрязнителей от источников к селитебным территориям, и др. Безусловно, эти и другие (рак-торы учитываются и ходе планирования
развития городов, однако существующие методики оценки их влияния не всегда точны. Например, в основе используемых методик расчета уровней загрязнения атмосферы городов лежат данные розы ветров, получаемые в результате статистической обработки наблюдений стационарных метеостанций. И то же время сложная орография территории, климатические и другие особенности региона могут существенно трансформировать ветровые потоки. Следовательно. такая локальная трансформация ветровых потоков приведет к более сложному распределению загрязняющих выбросов внутри городской агломерации. Роза ветров, получаемая на основе информации с фоновой метеорологической станции, обычно располагающейся на окраине города, может не совпадать с розой ветров для отдельных районов городской застройки [в: 16]. Это может быть обусловлено не только влиянием сложной орографии территории, но и фактически имеющейся городской застройкой.
Для изучения взаимосвязей между погодными процессами, антропогенным вмешательством и качеством атмосферы I! городских агломерациях широко применяются методы численного моделирования [2; 5; 7—!>; 14; 19], позволяющие на основе численного эксперимента провести анализ различных сценариев развития метеоси-туаций, спрогнозировать и оценить па его основе риски превышений предельно допустимых концентраций загрязнений при различных метеорологических условиях.
Численное моделирование может выступать инструментом для получения оценок вклада различных источников в загрязнение городской атмосферы, при этом источники могут быть как природными, так и антропогенными.
I) настоящей работе представлена технология проведения сценарных оценок качества атмосферы в областях со сложным рельефом, которая использовалась при исследованиях процессов переноса примесей в атмосфере г. Чита.
Методология проведения сценарных оценок и используемые данные. При изучении процессов формирования неблагоприятных метеорологических ситуаций I; районе г. Чита мы использовали сценарный подход, суть которого заключается в том, чтобы описать мезоклимат исследуемой территории набором сценариев. Результаты расчетов по этим сценариям будут отражать особенности развития метеорологических ситуаций, типичных для исследуемых регионов. Па основе сценарного моделирования можно проанализировать и варианты состояния атмосферы, приводящие к экстремально высоким уровням ее загрязнения.
11рп моделировании мы выделяем крупномасштабные атмосферные процессы, в большей мере связанные с годовым вращением Земли вокруг Солнца, и циркуляции локального характера, связанные с суточным вращением вокруг своей оси. Крупномасштабные поля метеоэлементов мы считаем известными и задаем их распределение на основе данных климатических справочников и результатов расчетов с помощью крупномасштабных моделей. Суточные изменения атмосферных характеристик мы рассматриваем как отклонения значений метеополей от их фоновых крупномасштабных значений. Отклонения отражают эволюцию метеорологических полей, вызванную факторами локального характера, прежде всего, рельефом местности и неод!юродностям и тем I юратуры подстилающей поверхности.
Для моделирования развития атмосферных циркуляции и процессов распространения примесей в Читино-Иигоди некой
впадине, в которой расположен г. Чита, мы использовали мезомасштабную негидростатическую модель динамики атмосферы и переноса примеси ft областях со сложным рельефом, разрабатываемую в ИВМиМГ СО РАН с непосредственным участием авторов. Основными уравнениями модели являются три уравнения движения для нахождения компонент вектора скорости ветра, уравнения переноса тепла и влаги, уравнения неразрывности. IÎ модели также рассчитывается тем i юратура i юдстилаю-щей поверхности. Для описания переноса пассивной примеси используется конвективно-диффузионное уравнение с начальными и граничными условиями.
Численные схемы для реализации модели построены на базе вариационных принципов с использованием сопряженных задач ¡8]. Вариационный подход обеспечивает энергетическую сбалансированность разностных схем, а также точный учет естественных краевых условий на границах области моделирования. Для учета орографии исследуемой территории используются идеи метода «фиктивных» областей [ I ]. Более детальное описание модели динамики атмосферы и переноса примеси и используемых методов численной реализации приведено в ряде работ [10: 1 I ; 17; 18 ].
Для проведения численных исследовании модель адаптировалась к орографии и другим характеристикам подстилающей поверхности выбранной территории. И нашей работе данные о высотах рельефа в узлах расчетной сетки получены с использовали ем ГИС-технологий 112] из находящихся в свободном доступе данных интерфероме-трической съемки земного шара — Shuttle Radar Topography Mission (SRTM) [15].
Геоинформационная часть информационно-вычислительного комплекса сценарного моделирования состоит из двух блоков. Один из них предназначен для получения входной информации для численной модели. Здесь строится координат-но привязанная к общегеографическому пространству двумерная расчетная сетка для области моделирования, из доступных цифровых моделей извлекаются данные о
рельефе, категориях землепользования и типов подстилающей поверхности области моделирования, вся полученная информация привязывается к узлам расчетной сетки н экспортируется 15 форматы, используемые в численном моделировании. Другой блок предназначен для отображения, обработки и анализа в среде ГИС результатов численных экспериментов. Здесь с использованием как системного, так и специально созданного авторского инструментария результаты моделирования привязываются к географическому пространству, строятся поля распределений исследуемых характеристик атмосферы. Геоинформационная обработка результатов сценариев численного моделирования и дальнейший их анализ позволяют построить карты пространственных распределений метеоэлементов и концентраций атмосферных примесей для каждого сценария . Это особенно важно для выявления далеко не очевидных участков, потенциально подверженных загрязнению, что позволяет определить места с наибольшим потенциалом изменений характеристик микроклимата и качества атмосферы.
Городская застройка, замерзшие или открытые водоемы и другие особенности территории учитываются в модели динамики атмосферы при задании соответствующих различий в значениях параметров, описывающих характеристики подстилающей поверхности: значения альбедо и параметра шероховатости поверхностей, различия начального распределения температуры н др. Эти различия имеют как пространственное, так и сезонное распределение.
Таким образом, проводя серии сценарных расчетов для каждого времени года, мы получаем п накапливаем данные о характере распределения полей метеоэлемеп-тов, описывающих особенности развития атмосферных циркуляции для региона. И дальнейшем эти данные можно использовать в качестве начальных условий при реализации более сложных сценариев развития метеоситуа1 цш.
Используемая модель и сценарный подход позволяют изучать закономерности формирования нолей концентрации при-
месей как от действующих, так и от потенциально возможных источников в регионе. Использование условных единиц в расчетах концентраций загрязнений позволяет, зная мощности выбросов каждого вида примесей для каждого из источников, достаточно просто получить оценочные значения концентраций и в натуральных единицах.
Для настройки модели к климатическим условиям региона необходимо задать распределение крупномасштабных полей метеоэлементов. В настоящей работе для расчета сценариев задавались различные варианты распределения метеополей, опирающиеся на значения климатических норм для соответствующих полей метеоэлемен-тов в каждом из исследуемых периодов года. Уже с таким набором параметров и распределений фоновых полей решалась задача установления квазипериодического режима, отражающего суточное поведение полей метеоэлеметов над выбранной территорией. Практически проверено, что для выхода на квазипериодический режим достаточно 2...3 сут модельного времени.
Результаты исследования и их обсуждение. Проведенные численные эксперименты можно условно разбить на две части. Первая из них касается сценариев, связанных с условиями, характерными для зимнего периода, коша территория региона подвержена воздействию Азиатского антициклона. И летний период для атмосферы г. Чита большую угрозу представляют лесные пожары в окрестностях города. Второй блок расчетных сценариев ориентирован па моделирование этих ситуаций.
Зимний сценарный расчет. Важная черта мезоклимата многих городов сибирского региона — формирование температурных инверсий (приземных и/или приподнятых) в зимний период [4|. Инверсионное распределение температуры затрудняет вертикальное перемешивание слоев воздуха, что приводит к застойным явлениям в атмосфере. Рельеф местности способствует формированию инверсионных условий. Исследованию влияния инверсионных условий на качество атмосферы в городах России посвящены работы [3: 13].
Для моделирования атмосферных циркуляции в окрестностях г. Чита выбрана облает!» 100*100 км- {рис. I). По вертикали расчеты велись до высоты 4,2 км. Па рис. 1 выделены городская застройка, водоемы и три рассматриваемых источника загрязнения атмосферы.
ж
<0 .УГ
Чита lu
рц I
, ГА- ъ
ш
В? - Щ ™
о
* у <*
! т * т
Рис. 1. Область моделирования / Fig. 1. Area of modelling
Численные расчеты проводились на сетке 200x200x50 узлов с шагами 500 м по горизонтали, по вертикали шаг был неравномерным: от 30 м над поверхностью, постепенно увеличивающийся к верхней границе области до 200 м.
Зимний сценарий рассмотрен со штилевыми ветровыми фоновыми условиями. Температура поверхности в начальный момент времени полагалась равном —25 "С. Стратификация атмосферы задавалась устойчивой.
В эксперименте предполагалось, что развитие циркуляции происходит только из-за температурных неоднородностей подстилающей поверхности в условиях сложного рельефа, поскольку фоновый набегающий поток отсутствует.
В рамках описанного сценария на фоне моделируемой метеорологической ситуации мы рассмотрели процесс переноса пассивных невесомых (не имеющих соб-
ственной скорости осаждения) примесей от трех точечных источников. Местоположение модельных источников на расчетной сетке задано максимально близко к локализации промышленных труб двух Т011 и котельной г. Чита. Источник № I имитирует трубу высотой 120 м над поверхностью земли, источники № 2 и № 3 соответствуют трубам высотой 90 м. В численном эксперименте мощность каждого источника задавалась равной I условной единице в секунду. Предполагалось, что в начальный момент во всей области концентрации примесей я вляются нулевым и.
Анализ результатов моделирования показал, что I; рамках рассматриваемого зимнего сценария I! Читино-Ингодинской котловине формируется инверсионное распределение температуры: подстилающая поверхность и прилегающий к ней воздух холоднее, чем слой воздуха, располагают щиеся выше. При таких условиях вертикальное перемешивание в нижних слоях затруднено, что способствует накоплению загрязняющих примесей от источников выбросов, располагающихся внутри инверсионного слоя. В нашем сценарии источники выбросов располагаются достаточно высоко над поверхностью, благодаря этому и очень слабому вертикальному рассеиванию загрязняющие выбросы остаются либо в верхней части инверсионного слоя, либо рассеиваются выше него.
На рис. 2 представлены изолинии концентраций пассивных примесей па высоте 30 м, 90 м и I 20 м над поверхностью в дневное время, когда температурная приземная инверсия наиболее ослаблена за счет нагрева подстилающей поверхности. Па рисунках видно, что концентрации на уровне 30 м значительно ниже, чем па высоте 90 и 120 м. Таким образом, сформировавшийся инверсионный слой может оказаться благоприятным метеорологическим фактором при наличии только высотных источников загрязнения, однако при отсутствии низких.
Также можно отметить, что от источника № 1 загрязнение преимущественно переносится на юг и юго-запад по долине, 1; противоположную сторону от основной го-
родской застройки, а от двух других источников примеси рассеиваются в юго-восточном направлении. Соответственно, можно сделать вывод, что в условиях рассмотренного зимнего сценария менее всего на ка-
Летние сценарные расчеты. Предложенная технология использовалась авторами и для изучения процессов распространения дымовых шлейфов от пожара в окрестностях г. Чита в летнее время. Область моделирования и параметры расчетное! сетки совпадают с областью и сеткой, использующимися в зимнем сценарии. Источник возгорания задавался в районе с. Сивяково, к юго-западу от г. Чита.
Здесь рассмотрены два сценарных расчета: сценарий с южным фоновым ветровым потоком 7 м/с и сценарий с юго-западным фоновым ветром 5 м/с.
В расчетах с фоновым ветром формирование атмосферных циркуляцпй происходит за счет неравномерного прогрева орографически сложной подстилающей поверхности и взаимодействия возникающих воздушных течений с фоновым набегающим потоком. В численных экспериментах значение скорости фонового ветра задавалось и фиксировалось на верхней границе расчетной области. Значения крупномасштабных фоновых полей метеоэлементов задавались в соответствии с данными о погодных условиях I! регионе 10—12 июля
чество атмосферы г. Чита влияют выбросы от источника № I. в то время как выбросы от источников № 2 и № 3 рассеиваются в атмосфере центральной н юго-восточной части города.
2015 г. В это время, по сообщениям прес-службы «Забмедиа» со ссылкой па Гослес-службу, в Читинском районе, недалеко от с. Сивяково уже несколько дней действовал пожар. Первый очаг возгорания обнаружен 7 июля 2015 г. К 10 июля столб дыма стал хорошо виден в г. Чита п. начиная сотого времени, дымовой шлейф начал распространяться над котловиной, в которой расположен город.
В численных экспериментах мы рассмотрели одиночный очаг возгорания. Начальная локализация пожара I! модели задавалась в соответствии с данными о пожарной ситуации на 12 июля 2015 г. В сценарных расчетах при параметризации разрастающегося очага пожара, как источника выброса моделируемой невесомой пассивной примеси, опирались на оценочные данные о скорости распространения огня. Так как динамика и тип пожара детально не моделировались, то для имитации выбросов от него в экспериментах полагалось, что непосредственно над очагом па высоте 30 м примесь имела начальную скорость вертикального распространения 1 м/с. Для простоты положено, что в каждой ячейке сеточной
Рис. 2. Изолинии концентраций пассивных примесей (усл. ед.) в 15 ч местного времени: а) на высоте 30 м над поверхностью; Ь) 90 м; с) 120 м / Fig 2. Isolines of passive pollutant concentration (conventional units) at 3 p.m. of local time: a) at the height of 30 m above the surface; b) 90 m; c) 120 m
области, находящейся в зоне возгорания, происходит выброс мощностью I условная единица в секунду.
Горизонтальные сечения полей примеси на высоте 30 м над подстилающей по-
верхностью, демонстрирующие результаты численных экспериментов по этим сценариям, представлены на рис. 3.
Рис. 3. Изолинии поля примеси на высоте 30 м над поверхностью: а) сценарий с южным фоновым ветром 7 м/с, Ь) сценарий с юго-западным фоновым ветром 5 м/с / Fig. 3. Isolines of passive pollutant concentration (conventional units) at the height of 30 m above the surface: a) south basic state wind 7 m/s;
b) south-western basic state wind 5 m/s
Анализ результатов расчетов для модельного сценария с южным фоновым ветровым потоком показал, что дымовой шлейф от выбранного нами случая пожара |; результате взаимодействия с локальными атмосферными циркуляциями, формирующимися в течение суток над неравномерно прогревающейся поверхностью, будет распространяться в основном на северо-запад поперек Читиио-Ингодинской впадины и перенесется через Яблоновый хребет. Перенос дымовой взвеси в северном направлении, захватывающий атмосферу над окраиной города, I; эксперименте наблюдался только в ранние утренние часы.
На рис. 3, I) представлены результаты расчетов для сценария с юго-западным внешним набегающим фоновым потоком. I Госкольку направление фонового ветра совпадает с основным направлением Чити-но-Иигоди некой впадины, то этот ветровой поток претерпевает суточные трансформации I! значительно меньшей степени, чем южный фоновый. В результате дымовой шлейф в продолжение всего модельного
времени переносился с различной степенью интенсивности, обусловленной ослаблением юго-западного ветра в дневные часы восходящими потоками над склонами долины, северо-восточном направлении на город.
Таким образом, результаты моделирования и основном соответствуют реальным картинам задымления городской атмосферы, которые регулярно наблюдают жители Читы летом.
Выводы. Анализ результатов численных экспериментов позволяет сделать вывод, что используемая модель в целом адекватно описывает как процессы формирования приземных н приподнятых инвер-сий в Читиио-Ингодинской впадине, так и процессы распространения загрязняющих примесей 15 ее атмосфере.
Эксперименты показали, что даже 15 отсутствие фонового набегающего потока изолинии концентраций загрязнений имеют достаточно сложный характер, обусловленный влиянием рельефа п температурных иеоднородностей подстилающей поверхности. Таким образом, можно констатиро-
вать, что местные локальные атмосферные циркуляции оказывают существенное влияние на перераспределение примесей над территорией города.
Используемая модель и сценарный подход позволяют детально изучать закономерности формирования полей концентрации примесей как от действующих, так и от потенциально возможных источников в регионе. Предположение об отсутствии собственной скорости осаждения у примесей
(невесомая примесь) позволяет изучать ветровые движения н их суточную трансформацию по деформации изолиний примесей па картах.
Использование методов математического моделирования и современных геоинформационных технологий особенно актуально при изучении изменений качества воздушной среды территории 15 результате перехода предприятий на новые технологии и новые виды топлива.
Список литературы
1. Ллоян Л. Е., Фалсйчик Л- А., Фалейчик Л. М. Алгоритм численного решения метеорологических задач в случае криволинейной области // Математические модели рационального природопользования. 11о-восибирск: Наука. Сиб. отд-е, 1989. С. 14—35.
2. Барт А. Д., Старченко А. В., Фаалиев А. 3. Информационно-вычислительная система для краткосрочного прогноза качества воздуха над территорией г. Томск // Оптика атмосферы и океана. 2012, Т. 25. № 7. С. 594-601.
3. Битехтина М, А., Михайлюта С, I!.. Леженин А. А., Тасейко О. Ii, Эволюция пограничного слоя и особенности загрязнения атмосферы города // Вестник Кемеровского государственного университета. 2012. Т. 4. J4b2. С. 143-148.
4. Климатические характеристики условий распространения примесей в атмосфере / под ред. I). Ю, Безуглой, VI. Ii. Берлин да. Ленинград: Гидрометеои.чдат. I !)83. 328 е.
5. Марчук Г, И. Математическое моделирование в проблеме окружающей среды, М,: Наука, М)82. 320 с.
6. Михайлюта С. В., Леженин А. Л., Тасейко О. В.. Битехтина М. А. Экологическая индустрия: ветровые потоки в городской застройке Красноярска // Инженерная экология. 2012. .№ 3. С. 26—37.
7. 11ененко Ii, Ii. Развитие концепции природоохранного прогнозирования для оперативного анализа и прогноза критических ситуаций в городских агломерациях // Интерэкспо Гео-Сибирь. 2017. Т. 4. № I, С. 70-73.
8. 11ененко В. В., Алоян А. Е. Модели и методы для задач охраны окружающей среды. 11овоснбирск: Наука, 1985. 256с.
9. Пененко В. В., Цветова Е. А. Математическое моделирование климато-экологическнх процессов урбанизированных территорий // Оптика атмосферы и океана. 2017. Т. 30, N> 06, С. 509—514,
10. Пьянова U. Д. Исследование трансформации воздушного потока над термически и орографически неоднородной подстилающей поверхностью // Вычислительные технологии. 2005, Т, 10. .№ ,'i. С. 106— III,
I I. Пьянова 0. А., Фалейчик .'I. М. Информационно-вычислительная технология для сценарных оценок динамики и качества атмосферы // Вычислительные технологии. 2012. Т. 17. .№ I. С. 109— 11 it.
12. Фалейчик Л. М. Геоинформационное обеспечение численного моделирования локальных атмосферных процессов // Вестник Новосибирского государственного университета. Сер. Информационные технологии. 2012. Т. 10. J& 2. С. 14-24.
13. Ячмёнева II. В., Гольвей А. 10, 11овторяемость инверсий и их влияние на уровень загрязнения атмосферного воздуха в г. Челябинске // Вестник Челябинского государственного униве|)ситета. 2011. .№ 5 (22(1), С. 84-89. '
14. Integrated Systems of Meso-Meteorological and Chemical Transport Models / ed. A. Baklanov, A. Makhura, Ii. Sokhi. Springer-Verlag Berlin I leidelberg, 2011.
15. Jarvis Д., Reuter II. I., Nelson A., Guevara Ii. 1 lole-filled SRTM for the globe. Version 4. 2008, Avail-able al: http://www,4nusrtm.csi.cgiar.org {Date of access: 12,09.2017).
16. Mikhailuta S. V., Lczbenin A.A., Pittd Д.. TaseäkoO. V, Urban wind Heids: Phenomena in transformation //Urban Climatc. 2017. Vol. 19. P. 122-140.
17. Penenko V.. Tsvetova Ii. IMscrete-analytical metbods for the Implementation of variational principles in environmental applications // J. of Сотри tation and Applied Mathematics. 2009. Vol. 226. Issue 2. P. 319—330.
18. Pyanova E, A., Penenko V. V.. Faleychik L. M. Scenario studies of local atmospheric circulations in t ho Krasnoyarsk region // Proc, SPIE 10035, 22nd International Symposium on Atmospheric and Ocean Optics: Atmospheric Physics. 100356D (November 29, 2016).
10. Pyanova E. A., Pcnenko V. V., Faleychik !.. M. Simulation of winter mesoelimates in Krasnoyarsk urban agglomeration // Proc. SPIE 10466, 23rd International Symposium on Atmospheric and Ocean Optics: Atmospheric Physics. 104666J (30 November 2017).
References_
1. Alovan A. E., Faleychik A. A., FaleychikL. M. Matematicheskiemodeliratsionalnogoprirodopolzovaniya (Mathematical models of rational environmental use), Novosibirsk: Nauka, I !)8!t. pp. 14—35.
2. Hart A. A., Starehen ко A. V., Fazliev A. Z. Optika atmosfery i океана {Atmospheric and occanic optics), 2012, no. 25 (7), pp. 504-601.
3. BitekhtinaM. A., MikhailutaS. V.. LezheninA. A.,Taseiko О. V. VeslnikKernerovskogogosudurslvennogo universiteta (Bulletin of the Kemerovo State University), 2012, no. 4(2), pp. 143—148.
4. Klimaticheskie khamkteristiki usloviy rasprostraneniya primesey v atmosfere (Climatic characteristics of the conditions for the propagation of impurities in the atmosphere) / ed, E. Yu. Bezuglay a, M. E. Berlyand). Leningrad, Gidrometeoizdat, 1983, 328 p.
Marchuk I. Matemalicheskoe modetirovanie v probleme okruzhayuschey sredy (Mathematical Modeling in the Problem of Environment). Moscow: Nauka, 1982. 320 p.
6. Mikhailuta S. V., Lezlienin A. A., Taseiko О. V., Bitekhtina M, A. Inzhenernaya ekologiya (Engineering ecology), 2012, no 3, pp. 26—37.
7. Pcnenko V. V, Interekspo Geo-Sibir (Interexpo Geo-Siberia), 2017, no. 4(1), pp. 70—73.
8. Pen en ко V. V., Aloyan A. E. Modell i melody dlya zadach ohrany okruzhayuschey sredy (Models and methods for environment protection problems). Novosibirsk: Nauka, 1985. 256 p.
9. Penenko V. V., Tsvelova E. A. Optika atmosfery i океана (Atmospheric and oceanic optics), 2017, no. 30 (6), pp. 500-514.
10. PyanovaE. A. Vychislitelnyetehnologii. (ComputationalTechnologies), 2005, vol. 10. no. 3, pp. 106— III.
11. Pyanova E. A., Faleychik L. M. Vychislitelnye tehnologii (Computational Technologies). 2012, no. 17 (1), pp. 109-119.
12. Faleychik L.M. VeslnikNovosibirskogogosudarstvennago universiteta. Ser. Informacionnye tehnologii (Novosibirsk State University Journal. Sor. Information Technologies), 2012, no. 10 (2), pp. 14—24.
13. Yachmeneva N. V., Golvey A. Yu. Vest nil; Chelyabinskogo gosudarstvermogo universiteta (Bulletin of the Chelyabinsk State University), 201 I, no. 5 (220), pp. 84-89.
14. Integrated Systems of .11eso-Meteorological and Chemical Transport Models (Integrated Systems of Me so-Meteorological and Chemical Transport Models) / cd. A. Baklanov, A. Makhura, Ii. Sokhi. Springer-Verlag Berlin Heidelberg, 2011.
15. JarvisA., Reuter II. I., Nelson A., Guevara E. Hole filled SRTM for the globe (Hole filled SRTM for the globe). Version 4. 2008. Available al: http://www.nnusrlm.esi.cjiiar.org (Date of access: 12.09.2017).
Hi. MikhailutaS. V., LezheninA. A., Pittd A., Taseiko О. V. Urban Climate (Urban Climate), 2017, vol. 19, pp. 122-140.
17. Penenko V,, Tsvelova E,./. of Computation and Applied Mathematics (J, of Computation and Applied Mathematics), 2009, vol. 22«, issue 2, pp. 319-330.
18. Pyanova E. A., Penenko V. V., Falcyehik L. M. Proc. SPIE 10035, 22nd International Symposium on Atmospheric and Ocean Optics: Atmospheric Physics {Proc. SPIE 10035, 22nd International Symposium on Atmospheric and Ocean Optics: Atmospheric Physics), 1003561) (November 29, 2016).
19. Pyanova E. A.. Penenko V, v., Falevehik L. M. Proc, SPIE 10406. 23rcl International Symposium on Atmospheric and Ocean Optics: Atmospheric Physics (Proc. SPIE 10466, 23i4l International Symposium on Atmospheric and Ocean Optics: Atmospheric Physics). 104666J (November 30, 2017).
Коротко об авторах_
Нышшш D.iii.tii Андреевна, канд. физ.-мат. наук, науч. сотрудник, Институт пычислител1.ИОЙ математики н математической геофизики Сибирского отделения Российской академии наук, г. Новосибирск, Россия. Область научных интересов: охрана окружающей среды, математическое моделирование, численные методы н эксперименты, природные ресурсы, географические информационные системы pianova@Bgs.ru
■ [■и ленчи к Андреи Анатольевич, канд. фш-мат. наук, доцент, Забайкальский государственный университет, г. Чита, России. Область научных интересов: математическое моделирование, численные методы и эксперименты, геоэкология faa55@bk.ru
([■нленчнБ Лариса Михайловна, канд. техн. наук, доцент, ст. науч. сотрудник, Институт природных ресурсов, экологии и криологии Сибирского отделения Российской академии наук, г. Чита, Россия: доцент, Забайкальский государственный университет, г. Чита, Россия. Область научных интересов: геоинформационные системы и технологии, геоэкология, математическое моделирование lfaleychik@bk.ru
Briefly about the authors___
lil/.a Рун nova, candidate of physical and mathematical sciences, research scientist, Institute of Computational Mathematics and Mathematical Geophysics, Russian Academy of Sciences, Siberian Branch ( I CM M G .Sli RAS), Novosibirsk, Russia. Sphere of scientific interests: environment protection, mathemalical modeling, calculus of approximations and experiments, iialural resources, geographic information system
Andrey l itlevehik, candidate of physical and mathematical sciences, associate professor, Transbaikal State University, Chita, Russia. Sphere of scientific interests: mathematical modeling, calculus of approximation and experiments, geoecology
I ¿irisa Fhleychik, candidate of technics I sciences, associate professor; senior research scientist, Institute of Natural Resources, Ecology and Cryology, Russian Academy of Sciences, Siberian Branch (INRECSB RAS), Cliita, Russia; associate professor, Transbaikal State University, Chita, Russia. Sphere of scientific interests: geographic information system and technologies, geoecology, mathematical modeling
Образец цитирования _
Пьянова Э. А., Фалейник А. А., ФалеймикЛ. М. Использование математического моделирования для оценки экологических последствий реиндустриализации экономики региона // Нести. Забайкал. гос. ун та. 2017. Г. 23. № 12. С. 41-50. DOI: 10.21209/2227 9245 2017-23-12-41 50.
Pyanova 11., Faleychïk A., Faleychik L. Use of mathematical modeling for estimation of environmental consequences of the regional economy ^industrialization // Thinsbaikal Stale University Journal. 2017, vol. 23, no. 12, pp. 41-50. DOI: 10.21209/2227-9245-2017-23-12-41-50.
Работа выполнена при частичной финансовой поддержке Программы фундаментальных исследований 11резидиума РАИ 1.33П, РФФИ, проект № 17-01-00137. Адаптация авторской модели для условий Забайкальского края выполнена в рамках проекта РГНФ, .№ 16-02-00102а. Расчеты выполнены с использованием ресурсов Ц1\П ССКЦИВМнМГСОРЛН.
Дата поступлении статьи: 08.12.2017 г. Дата опубликовании статьи: 25.12.2017 г.
