CC BY
31
УДК 502.3:504.61
ОЦЕНКА ЭКОЛОГИЧЕСКИХ РИСКОВ РЕИНДУСТРИАЛИЗАЦИИ НА ЮГО-ВОСТОКЕ ЗАБАЙКАЛЬСКОГО КРАЯ
REINDUSTRIALIZATION IN SOUTH-EASTERN TRANSBAIKALIE: ENVIRONMENTAL RISK ASSESSMENT
Э. А. Пьянова,
Институт вычислительной
математики и математической геофизики СО РАН, г. Новосибирск [email protected]
E. Pyanova,
Institute of Computational
Mathematics and Mathematical Geophysics, Russian Academy of Sciences, Siberian Branch (ICMMG SB RAS), Novosibirsk.
П
Инг
Л. М. Фалейчик,
Институт природных ресурсов,
экологии и криологии СО РАН, Забайкальский государственный университет, г. Чита [email protected]
L. Faleychik,
Institute of Natural Resources, Ecology and Cryology, Russian Academy of Sciences, Siberian Branch (INREC SB RAS), Transbaikal State University, Chita
А. А. Фалейчик,
Читинский институт Байкальского государственного университета, Забайкальский государственный университет, г. Чита [email protected]
A. Faleychik,
Chita Institute of the Baikal State University, Transbaikal State University, Chita
Перспективы реиндустриализации в Забайкальском крае в первую очередь связывают с инвестиционными проектами в горнодобывающей отрасли края. Их реализация добавляет новые и увеличивает существующие экологические риски. На Юго-Востоке края создается крупный горнопромышленный кластер. Необходимы исследования по оценке последствий воздействия его объектов на окружающую среду и расположенные здесь особо охраняемые природные территории.
На базе комплексного использования методов математического моделирования локальных атмосферных процессов и современных геоинформационных технологий, на основе сценарного подхода проведена серия численных экспериментов для прогнозной оценки возможного загрязнения атмосферы территории Юго-Востока выбросами котельной Быстринского ГОКа — первого из запланированных к вводу объектов кластера. Выполнен многовариантный анализ развития локальных атмосферных циркуляций, типичных для исследуемой территории в зимний период. Моделировалось распространение двух видов пассивной примеси: «легкой» (невесомой) и «тяжелой» (со средней скоростью оседания 0,4 м/с). Анализ результатов сценарных расчетов показал, что в зимний период вынос «тяжелых» загрязнений от котельной ГОКа за его пределы будет незначительным; северо-западный фоновый поток будет способствовать накоплению относительно небольших концентраций загрязнений у границ Борзинского заказника
Ключевые слова: реиндустриализация; Юго-Восток .Забайкальского края; экологические риски; особо охраняемые природные территории (ООПТ); оценка антропогенного воздействия; горнопромышленный комплекс (ГПК); загрязнение атмосферы; математическое моделирование локальных атмосферных процессов; географические информационные системы (ГИС); геоинформационные технологии
Reindustrialization perspectives in Transbaikal region are associated primarily with investment projects in the regional mining. Their implementation adds new environmental risks and enhances the existing ones. A major mining cluster is created in the South-East of the Transbaikal region. The research is needed to assess the effects of its facilities on the environment and nature protected areas located here.
A series of numerical experiments was conducted for prognostic assessment of air pollution by emissions of Bystrinsky mining factory (BMF) boiler in the South-East territory. The experiments were done on the basis of mathematical modeling of the local atmospheric processes and use of modern geo-informational technologies, also a scenario approach was applied. A multivariate analysis of the local atmospheric circulations typical for this area in winter was carried out. A distribution of two types of passive atmospheric pollutants was simulated: «light» (weightless) and «heavy» (with an average sedimentation rate of 0,4 m/s).
The results of the analysis scenario calculations have showed that «heavy» passive pollutant from the BMF boiler beyond it will be negligible in winter; but northwest basic state thread will contribute to the accumulation of relatively small concentrations of contaminants at the borders of Borzinsky reserve area
Key words: reindustrialization; South-East Transbaikalie; ecological risks; nature protected areas (NPA); estimation of human impact; mining complex; air pollution; mathematical modeling of local atmospheric processes; Geographic Information System (GIS); GIS technology
Работа выполняется при частичной поддержке РГНФ, проект № 16-02-00102а (экологические риски реиндустриализации для Юго-Востока Забайкальского края, интерпретация результатов сценарных расчетов), РФФИ, проект № 14-01-00125-a (численные алгоритмы), Программы фундаментальных исследований Президиума РАН 1.33П (математическая модель, авторская модификация и адаптация), проекта IX.88.1.6 Фундаментальных исследований СО РАН (геоинформационный инструментарий для получения и обработки данных моделирования). Расчеты выполнены с использованием ресурсов ЦКП ССКЦ СО РАН
Концепция реиндустриализации экономики России активно обсуждается не только в научных кругах, но и в российских средствах массовой информации. Несмотря на разницу подходов и убедительность используемых аргументов, все мнения сходятся к необходимости и безальтерна-тивности курса на реиндустриализацию (неоиндустриализацию) [5; 9]. Ядро этой концепции — идея резкого скачкообразного роста промышленного потенциала страны за счет глубокой структурной перестройки и модернизации экономики на базе качественного технологического обновления материального производства, преимущественного развития современных высоких технологий, интеграции научной, производственной и образовательных сфер [ 1].
Долгосрочные перспективы развития Забайкальского края связывают в основном с реализацией инвестиционных проектов в приоритетных отраслях экономики, определенных Стратегией социально-экономического развития Забайкальского края на период до 2030 г., в первую очередь в минерально-сырьевом комплексе и горно-
рудной промышленности. Таким образом, наиболее вероятным путем новой индустриализации применительно к Забайкальскому краю может стать обращение, следуя основным положениям Программы «Сибирь», на новом современном технологическом уровне к проектам освоения природных богатств территории. Такой сценарий развития реиндустриализации вполне возможен: на территории края находится целый ряд привлекательных для инвесторов месторождений полезных ископаемых [7; 10]. Ключевым приоритетом развития горнодобывающей отрасли объявлено создание горнорудных комплексов, которые должны повысить эффективность и комплексность использования природных ресурсов. Однако такой путь развития существенно увеличивает экологические риски. Поэтому требуется выявление и очень существенная их проработка на всех этапах освоения территории; необходимо достигать компромисса экономических и экологических интересов, чтобы исключить экологически неравноценный обмен или хотя бы снизить его риски [2; 3].
В настоящее время на территории края достаточно быстро формируются и развиваются два горнопромышленных кластера. Один из них — на севере края, в труднодоступном районе с очень ранимыми и слабо восстанавливающимися экосистемами. В развитии этого кластера существенную роль играет Байкало-Амурская магистраль. Другой горнопромышленный кластер — на более доступной территории юго-востока края. В первоначальном варианте здесь планировалось создание Забайкальского территориального горно-металлургического комплекса из пяти ГОКов. Одним из первых объектов горнопромышленного комплекса на Юго-Востоке будет строящийся Быстринский ГОК, срок ввода в эксплуатацию которого намечен на 2017 или 2018 г.
С предыдущих этапов освоения Юго-Востока Забайкальского края (ЮВЗ) уже накоплен целый комплекс экологических проблем: почти вся его территория в той или иной степени связана с функционированием предприятий горнопромышленного комплекса (ГПК). Кроме того, здесь находятся десятки заброшенных и законсервированных объектов: нерекультиви-рованные карьеры, отвалы и отстойники, представляющие собой опасные источники вредных элементов и деградации земель.
Вместе с тем Юго-Восточное Забайкалье выделяется особым своеобразием природных комплексов. Здесь расположены особо охраняемые природные территории (ООПТ) федерального и регионального значений [4], которые могут оказаться в зоне существенного экологического риска. Здесь значительно выше риск необратимых изменений природных систем и их экологических функций. Поэтому оценка воздействия как действующих, так и планируемых горнорудных предприятий на природные комплексы приобретает особый смысл.
На предыдущих этапах исследований оценены масштабы воздействия на природные комплексы объектов горнорудной промышленности Юго-Востока Забайкалья, оставшихся на территории со времен предыдущего освоения [7]. С использованием ГИС-технологий выявлены и оценены тер-
ритории, в той или иной степени уже подвергшиеся антропогенному влиянию. Однако проведенные исследования в меньшей степени касаются вопросов, связанных с загрязнением воздушного бассейна территории выбросами предприятий ГПК. Учет информации об атмосферном загрязнении, которое в значительной степени зависит от характера рельефа местности и проявлений особенностей мезоклимата, позволит уточнить границы областей влияния объектов ГПК на окружающую территорию.
В данной статье особое внимание уделено использованию методов математического моделирования, разработанных для решения задач охраны окружающей среды [6; 8], в комплексе с возможностями современных геоинформационных технологий не только для анализа существующего положения дел, но и для прогнозных оценок вариантов развития ситуаций. Центральным ядром используемого подхода является информационно-вычислительная технология получения, обработки и анализа многофакторной информации о территории, начиная с получения исходных данных для численного моделирования и заканчивая пространственным анализом и представлением результатов сценарных расчетов по математическим моделям в виде, удобном и понятном для конечных пользователей, принимающих решения.
Для математического описания гидротермодинамических режимов атмосферы и процессов переноса примеси над выбранной территорией использовалась разработанная в ИВМиМГ СО РАН мезомасштабная негидростатическая модель. Основными ее уравнениями являются три уравнения движения; уравнения для переноса тепла и влаги; уравнение неразрывности; уравнение состояния; уравнение конвекции-диффузии пассивной примеси. Взаимодействие атмосферы с подстилающей поверхностью учитывается с помощью схем параметризации в рамках модели приземного слоя. Весь комплекс моделей адаптирован авторами к сложным орографическим и климатическим условиям Юго-Востока Забайкальского края.
Задача атмосферной динамики и переноса пассивной примеси реализовывалась численно на основе вариационно-разностного подхода [8]. В результате построены энергетически сбалансированные численные схемы второго порядка точности по пространству. Для аппроксимации конвективно-диффузионных операторов задачи использовались монотонные дискретно-аналитические схемы [12]. Вариационный подход обеспечивает точный учет естественных краевых условий на границах области моделирования. Более детальное описание моделей и методов построения конечно-разностных уравнений можно найти в [11] и других работах.
Геоинформационная часть информационно-вычислительного комплекса для получения сценарных оценок включает в себя два блока. Первый блок предназначен для получения входной информации для численной модели — построение в модельной области координатно привязанной к географическому пространству двумерной расчетной сетки, извлечение из существующих цифровых моделей данных о релье-
фе, категориях землепользования и типов подстилающей поверхности области моделирования. Второй блок — для отображения и анализа результатов расчетов в среде ГИС. Экспорт в геоинформационную среду результатов сценариев численного моделирования и дальнейший их анализ позволяет получить карты пространственных распределений метеоэлементов и концентраций атмосферных примесей для каждого сценария; выявить зоны, в том числе и далеко не очевидные участки, потенциально подверженные загрязнению; определить места с наибольшим потенциалом изменений характеристик микроклимата и качества атмосферы.
В статье даны некоторые предварительные оценки возможного атмосферного загрязнения в результате строительства первого из запланированных объектов горнопромышленного кластера на Юго-Востоке — Быстринского ГОКа. Сценарные расчеты и оценки получены для территории (расчетной области) 100 х 100 км2, представленной на рис. 1.
Рис. 1. Территория Забайкальского края, объекты проектируемого горнопромышленного кластера
и область моделирования
Fig. 1. The territory of Transbaikal region, the objects of the projected mining cluster and the modeling area
По оценкам [7], в результате строительства ГОКа зона полного разрушения природных комплексов составит не менее 120 км2 и будет непосредственно примыкать к северной границе Борзинского заказника. Ее границы будут находиться в 25 км от юго-западной границы Урюмканского заказника и примерно в 35 км от восточной границы Туровского заказника.
В силу природно-климатических особенностей Забайкалья наиболее опасным с точки зрения возникновения высоких уровней загрязнения воздушного бассейна является зимний период. Зимой Забайкальский край и его юго-восток находятся под действием отрога Сибирского антициклона. Это обусловливает высокую повторяемость температурных инверсий и маловетреной погоды, что, в свою очередь, способствует застойным явлениям, особенно в котловинах. По этой причине при разработке сценариев численных экспериментов в первую очередь рассмотрен зимний период.
Для численного моделирования атмосферных циркуляций и переноса атмосферных примесей в области моделирования (рис. 1) использовалась расчетная сетка 251х251х50 узлов с шагами по горизонтали Ах = Ау = 400 м, по вертикали шаг Az от 30 м над рельефом до 200 м в верхней части расчетной области. Расчеты по вертикали велись до высоты порядка 4200 м над поверхностью.
При организации сценариев численных экспериментов рассмотрены типичные для юго-восточной части края варианты метеорологических ситуаций. Поскольку и зимой, и летом для этой территории характерна большая повторяемость штилей, преобладают западные, северо-западные и юго-западные ветры, то мы рассмотрели следующие зимние сценарии: (1) штиль при устойчивой стратификации фоновой атмосферы; (2) штиль при инверсионном распределении температуры в нижних слоях фоновой атмосферы; (3) юго-западный фоновый ветер 4 м/с при инверсионном распределении температуры в нижних слоях фоновой атмосферы; (4) северо-западный фоно-
вый ветер 4 м/с при устойчивой стратификации фоновой атмосферы.
Моделирование штилевых сценариев (1) и (2) подразумевает, что внешнего по отношению к рассматриваемой территории набегающего ветрового потока нет, и все атмосферные циркуляции в условиях сложной орографии изучаемой области формируются только за счет неравномерного прогрева подстилающей поверхности. В сценариях с фоновым потоком (3) и (4) задавались начальные распределения поля скорости ветра по высоте.
При реализации сценарных расчетов воспроизводился суточный ход метеорологических параметров и на их фоне рассчитывались процессы переноса пассивной примеси от точечного источника загрязнений. В качестве источника выбросов в наших сценариях рассматривалась труба котельной высотой 60 м. Выбросы при взрывных работах и пыление отвалов не учитывались. Во всех экспериментах моделировалось распространение двух видов пассивной примеси: «легкой» (невесомой) и «тяжелой» (со средней скоростью оседания 0,4 м/с). К первой можно отнести выбросы СО, вторая может имитировать неорганическую пыль, составляющую значительную долю всех выбросов котельной. Во всех экспериментах мощность источника задавалась как 1 условная единица выброса за 1 единицу времени (сек); шаг по времени А^ = 60 с. В начальный момент времени поле примеси по всей расчетной области полагалось нулевым.
Результаты расчетов полей вектора скорости ветра и концентраций пассивных примесей по всем четырем зимним сценариям представлены на рис. 2...6, изолинии на которых отражают падение концентраций примеси в 10 раз.
Как показали результаты численных экспериментов для штилевых сценариев, невесомая примесь имеет заметное вертикальное и горизонтальное рассеивание, в отличие от «тяжелой» (рис. 2), даже низкие концентрации которой локализуются в пределах долины, где располагается Быстрин-ский ГОК. Численные оценки масштабов распространения невесомой примеси от точ-
ки выброса показали, что изолиния со значением 1 достигает высот 200...400 м над поверхностью. В предположении, что суммарная мощность выбросов СО котельной составит примерно 40,4 г/с, то изолинии со значением 1 усл. ед. будет соответствовать концентрация СО порядка 0,51 мг/м3.
Сравнивая результаты расчетов по штилевым сценариям (1) и (2), можно
отметить, что слой рассеивания примеси при фоновом инверсионном распределении температуры имеет меньшую высоту, чем в сценарии без инверсии (рис. 3). Такое поведение поля примеси было ожидаемо, поскольку инверсия является тормозящим фактором при развитии вертикального перемешивания.
Рис. 2. Вертикальные разрезы полей вектора скорости ветра и концентраций: а) невесомой пассивной и б) «тяжелой» примесей через 9 ч с момента включения источника выбросов;
сценарий(1)
Fig. 2. Vertical sections of fields of wind velocity vector and concentration: a) weightless passive and b) «heavy» pollutants after 9 hours from the moment of the emission
source appearance; scenario (1)
в
N
2
35 — 3 iil/c
Y, km
2 -
1 - :
Jljll m j|
35 _ J m/c 40 45 50 55
Y, km
Рис. 3. Вертикальные разрезы полей вектора скорости ветра и концентраций невесомой пассивной примеси через 9 ч с момента включения источника выбросов: а) сценарий (1); б) сценарий (2)
Fig. 3. Vertical sections of fields of wind velocity vector and concentrations of passive tracer after 9 hours from the moment of the emission source appearance: a) scenario (1); b) scenario (2)
Кроме того, расчеты для штилевых сценариев показали, что «тяжелая» примесь практически не поднимается выше окружающих склонов (рис. 2, б) и ее большая часть оседает в радиусе 1...1,5 км от источника выбросов. Особенно это характерно для моделируемой ситуации с атмосферной инверсией (рис. 4, б). И даже фоновый ветровой поток, как показали расчеты по
сценарию с юго-западным ветром (3), не способен при таких инверсионных условиях значительно усилить вынос загрязнения из долины, в которой располагается котельная (рис. 5, б). Изолинии со значением 1 усл. ед. для концентрации «тяжелой» примеси (пыли) соответствуют 0,7 мг/м3, если предположить по оценочным данным, что мощность источника пыли 57 г/с.
Рис. 4. Поля «легкой» (а) и «тяжелой» (б) примесей на высоте 30 м над поверхностью через 32 ч с момента включения источника выбросов; сценарий (1)
Fig. 4. Fields of «light» (a) and «heavy» (b) pollutants at a height of 30 m above the surface after 32 hours from the moment of the emission source appearance; scenario (1)
Рис. 5. Поля «легкой» (а) и «тяжелой» (б) примесей на высоте 30 м над поверхностью через 32 ч с момента включения источника выбросов; сценарий (3)
Fig. 5. Fields of «light» (a) and «heavy» (b) pollutants at a height of 30 m above the surface after 32 hours from the moment of the emission source appearance; scenario (3)
В модельных расчетах по сценарию (3) наблюдается интересный эффект формирования локальных атмосферных цир-куляций. Несмотря на юго-западный фоновый ветер, ожидаемый перенос примесей в
северо-восточном направлении к Урюмкан-скому заказнику практически отсутствует (рис. 5). Юго-западный фоновый поток трансформируется за счет взаимодействия с поверхностью, в долине р. Ильдикан об-
разуются ветровые потоки, направление которых совпадает с направлением вытя-нутости межгорной долины с востока на запад. Этот локальный ветер и способствует переносу «тяжелой» и «легкой» примесей в западном направлении против основного фонового потока.
Численные эксперименты по сценарию с северо-западным фоновым ветром (4) показали, что основные концентрации как «тяжелой», так и «легкой» примесей тоже локализуются в радиусе 1.. .1,5 км от трубы котельной (рис. 6).
Рис. 6. Изолинии пассивных примесей на высоте приземного слоя через 16 ч с начала работы источника:
а) «легкая» и б) «тяжелая» примеси; сценарий (4) Fig. 6. Isolines of passive pollutants at the height of the surface layer after 16 hours since the beginning of
operation of the source: a) «easy» and b) «heavy» pollutants; scenario (4)
При этом северо-западный фоновый поток способствует переносу загрязнений к границам Борзинского заказника, и хотя в случае «тяжелой» примеси, как показали численные расчеты, границ этой ООПТ могут достигать только очень незначительные концентрации (рис. 6, б), за продолжительный период времени это может привести к уже заметному негативному антропогенному воздействию на природный комплекс заказника. Концентрации «легкой» примеси, в предположении нулевой начальной фоновой концентрации, у границ Борзинского заказника составили 0,051 мг/м3 (0,1 усл. ед.). Так как в представленных оценках учитывалась только работа котельной Быстринского ГОКа, пыление отвалов и проведение взрывных работ здесь не рассматривалось, то, учитывая и их в качестве источников загрязнения атмосферы, можно говорить о возможном увеличении значений концентраций примесей вблизи границ заказника. Поэтому эта метеоситуация требует более детального изучения и учета при планировании работ ГОКа.
Выводы. Таким образом, в представленной статье приведены результаты начального этапа исследований по прогнозной оценке состояния и качества атмосферы Юго-Востока Забайкальского края, где продолжают осваиваться рудные запасы его недр и планируются к вводу новые горнодобывающие и перерабатывающие комплексы. Первый из заявленных к вводу в эксплуатацию предприятий ГПК — Быстринский ГОК. Добыча руды открытым способом и переработка сырья будут, очевидно, сопровождаться загрязняющими окружающую среду выбросами, в том числе и в атмосферу.
Анализ результатов модельных расчетов переноса загрязняющих выбросов котельной ГОКа (другие источники загрязнений здесь не рассматривались) по перечисленным зимним сценариям показал, что при штилевых фоновых условиях в зимний период вынос «тяжелых» загрязнений (например, неорганической пыли) от котельной ГОКа за его пределы будет незначительным. Однако такие метеороло-
гические условия будут неблагоприятными для проветривания долины, в которой располагается ГОК. При юго-западном фоновом ветре выраженного переноса примесей в направлении Урюмканского заказника не наблюдалось: за счет рельефа формируется локальный ветер, способствующий переносу «тяжелой» и «легкой» примесей в западном направлении против основного фонового потока. Метеоситуция с северо-западным фоновым ветром — неблагоприятная для природных систем Борзинского заказника: фоновый поток способствует переносу загрязнений к их границам. Хотя основные концентрации как «тяжелой», так и «легкой» примесей локализуются в радиусе 1.1,5 км от трубы котельной, в незначительных концентрациях примеси могут достигать границ данной ООПТ, и за продолжительный период времени это может привести к накоплению негативного воздействия на ее природный комплекс.
Следует отметить, что рассмотренные сценарии далеко не исчерпывают все особенности формирования и развития метеорологических режимов над исследуемой областью. Предполагается продолжить
Список литературы_
исследования сценарными расчетами процессов формирования локальных атмосферных циркуляций и переноса примесей в различные сезоны, а также дальнейшее изучение процессов атмосферного переноса загрязняющих примесей и от других объектов ГОКа.
Широкие возможности комплексиро-вания математического и геоинформационного моделирований могут быть использованы и для решения многих других задач выявления и оценки экологических рисков. На очереди следующий проект в русле ре-индустриализации Юго-Востока Забайкальского края — освоение Бугдаинского месторождения рудного золота, молибдена, свинца и серебра, строительство на его базе Бугдаинского ГОКа. В планах дальнейших исследований авторов — прогнозная оценка загрязнения атмосферы окружающей территории выбросами этого объекта, а также оценка кумулятивного воздействия на качество атмосферы территории Юго-Востока этих и остальных планируемых к вводу ГОКов данного горнопромышленного кластера.
1. Бодрунов С. Д., Гринберг Р. С., Сорокин Д. Е. Реиндустриализация российской экономики: императивы, потенциал, риски // Экономическое возрождение России. 2013. № 1. С. 19—49.
2. Глазыгрина И. П. Принцип «загрязнитель платит» и экологически неравноценный обмен // Вестник ЗабГУ. 2016. Т. 22. № 3. С. 93-100.
3. Глазырина И. П., Забелина И. А., Клевакина Е. А. Экологическая составляющая экономического развития: приграничные регионы России и Китая // ЭКО. 2014. № 6. С. 5-24.
4. Кирилюк О. К. Эколого-географические основы развития и современное состояние сети особо охра-няемыгх природныгх территорий Восточного Забайкалья / / Вопросы современной науки и практики. 2009. № 8 (22). С. 144-151.
5. Мартынов А. В. Новая индустриализация: взаимодействие экономической и социальной политик // Проблемы теории и практики управления. 2014. № 2. С. 25-34.
6. Марчук Г. И. Математическое моделирование в задачах охраны окружающей среды. М.: Наука, 1982. 320 с.
7. Минерально-сыгрьевой сектор Азиатской России: как обеспечить социально-экономическую отдачу? / под ред. акад. РАН В. В. Кулешова. Новосибирск: ИЭОПП СО РАН, 2015. 352 с.
8. Пененко В. В., Алоян А. Е. Модели и методы для задач охраны окружающей среды. Новосибирск: Наука, 1985. 256 с.
9. Сухарев О. С. Реиндустриализация экономики России и технологическое развитие // Национальные интересы: приоритеты и безопасность. 2014. № 10. С. 2-16.
10. Glazyrina I. P., Kalgina I. S., Lavlinskii S. M. Problems in the Development of the Mineral and Raw Material Base of Russia's Far East and Prospects for the Modernization of the Region's Economy in the Framework of Russian-Chinese Cooperation // Regional Research of Russia. 2013. Vol. 3. No. 4. Р. 405-413.
11. Penenko V. V., Pyanova E. A., Faleychik L. M. Simulation of atmospheric dynamics and air quality in the Baikal region / / Proc. SPIE 9292, 20th International Symposium on Atmospheric and Ocean Optics: Atmospheric Physics, 929247 (November 25, 2014); doi:10.1117/12.2074998.
12. Penenko V., Tsvetova E. Discrete-analytical methods for the implementation of variational principles in environmental applications // J. of Computation and Applied Mathematics. 2009. Vol. 226, iss. 1. P. 319—330.
List of literature_
1. Bodrunov S. D., Grinberg R. S., Sorokin D. E. Ekonomicheskoe vozrozhdenie Rossii (Economic revival of Russia), 2013, no 1, pp. 19-49.
2. Glazyrina I. P. Vestn. Zab. Gos. Univ. (Transbaikal State University Journal), 2016, vol. 22, no 3, pp. 93-100.
3. Glazyrina I. P., Zabelina I. A., Klevakina E. A. EKO (ECO), 2014, no. 6, pp. 5-24.
4. Kirilyuk O. K. Voprosy sovremennoy nauki i praktiki. Universitet im. V.I. Vernadskogo (Problems of Contemporary science and practice. Vernadsky University), 2009, no. 8 (22), pp. 144-151.
5. Martynov A. V. Problemy teorii i praktiki upravleniya (Theoretical and practical aspects of management), 2014, no. 2, pp. 25-34.
6. Marchuk G. I. Matematicheskoe modelirovanie vzadachah ohrany okruzhayushhey sredy [Mathematical modeling in environmental problems]. Moscow: Nauka, 1982. 320 p.
7. Mineralno-syrevoy sektor Aziatskoy Rossii: kak obespechit sotsialno-ekonomicheskuyu otdachu? [The mineral-raw material sector of the Asian Russia: how to ensure the socio-economic return?]; ed. by V. V. Kule-shov. Novosibirsk: IEIE SB RAS, 2015, 352 p.
8. Penenko V. V., Aloyan A. E. Modeli i metody dlya zadach ohrany okruzhayuschey sredy [Models and methods for environment protection problems]. Novosibirsk: Nauka, 1985. 256 p.
9. Sukharev O. S. Natsionalnye interesy: prioritety i bezopasnost (National interests: priorities and security), 2014, no. 10, pp. 2-16.
10. Glazyrina I. P., Kalgina I. S., Lavlinskii S. M. Regional Research ofRussia (Regional Research of Russia), 2013, vol. 3, no. 4, pp. 405-413.
11. Penenko V. V., Pyanova E. A., Faleychik L. M. Simulation of atmospheric dynamics and air quality in the Baikal region (Simulation of atmospheric dynamics and air quality in the Baikal region): Proc. SPIE 9292, 20th International Symposium on Atmospheric and Ocean Optics: Atmospheric Physics, 929247 (November 25, 2014); doi:10.1117/12.2074998.
12. Penenko V., Tsvetova E. J. of Computation and Applied Mathematics (J. of Computation and Applied Mathematics), 2009, vol. 226, iss. 1, pp. 319-330.
Коротко об авторах_
Пьянова Эльза Андреевна, канд. физ.-мат. наук, научный сотрудник, Институт вычислительной математики и математической геофизики СО РАН (ИВМиМГ СО РАН), г. Новосибирск, Россия. Область научных интересов: охрана окружающей среды, математическое моделирование, численные методы и эксперименты, природные ресурсы, географические информационные системы pianova@ngs. ги
Фалейчнк Лариса Михайловна, канд. техн. наук, доцент, старший научный сотрудник, Институт природных ресурсов, экологии и криологии СО РАН (ИПРЭК СО РАН), доцент, Забайкальский государственный университет, г. Чита, Россия. Область научных интересов: геоинформационные системы и технологии (ГИС), геоэкология, математическое моделирование [email protected]
Фалейчнк Андрей Анатольевич, канд. физ.-мат. наук, доцент, Читинский институт Байкальского государственного университета, доцент, Забайкальский государственный университет, г. Чита, Россия. Область научных интересов: математическое моделирование, численные методы и эксперименты, геоэкология [email protected]
Briefly about the authors_
Elza Pyanova, candidate of physical and mathematical sciences, research scientist, Institute of Computational Mathematics and Mathematical Geophysics, Russian Academy of Sciences, Siberian Branch (ICMMG SB RAS), Novosibirsk, Russia. Sphere of scientific interests: environment protection, mathematical modeling, calculus of approximations and experiments, natural resources, GIS
Larisa Faleychik, associate professor, candidate of technical sciences; senior research scientist, Institute of Natural Resources, Ecology and Cryology, Russian Academy of Sciences, Siberian Branch (INREC SB RAS); associate professor, Transbaikal State University, Chita, Russia. Sphere of scientific interests: Geographic Information System (GIS), GIS technologies, geoecology, mathematical modeling
Andrey Faleychik, associate professor, candidate of physical and mathematical sciences, associate professor, Chita Institute of the Baikal State University, associate professor, Transbaikal State University, Chita, Russia. Sphere of scientific interests: mathematical modeling, calculus of approximation and experiments, geoecology
Образец цитирования _
Пьянова Э.А., Фалейчик Л.М., Фалейчик A.A. Оценка экологических рисков реиндустриализации на Юго-Востоке Забайкальского края // Вестн. Забайкал. гос. ун-та. 2016. Т. 22. № 9. С. 22—32.
