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

Доклады БГУИР

2017, № 7

УДК 004.94: 004.9.032.26

Doklady BGUIR

2017, No. 7

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ РАСПРОСТРАНЕНИЯ ЗАГРЯЗНЯЮЩИХ ВЕЩЕСТВ ОТ ЭКОЛОГИЧЕСКИ ОПАСНЫХ ОБЪЕКТОВ С УЧЕТОМ ПОВЕРХНОСТНОГО СТОКА И ИНФИЛЬТРАЦИИ

П К. ШАЛЬКЕВИЧ1, С П. КУНДАС2, И.А. ГИШКЕЛЮК1

Международный государственный экологический институт имени А. Д. Сахарова Белорусского государственного университета, Республика Беларусь

2Белорусский национальный технический университет, Республика Беларусь.

Поступила в редакцию 19 сентября 2017

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

Ключевые слова: математическое моделирование, численные методы, поверхностный сток, инфильтрация.

Abstract. A mathematical model of pollutants' migration from environmentally hazardous facilities considering surface water flow and infiltration is developed. The methods for the numerical solution of differential equations describing the transport of pollutants and their entry into the groundwater from the surface of the soil due to the infiltration of precipitation are introduced. The developed mathematical model was implemented with the use of MATLAB and COMSOL software.

Keywords: mathematical modeling, numerical methods, surface water flow, infiltration.

Doklady BGUIR. 2017, Vol. 109, No. 7, pp. 45-52 Mathematical modeling of pollutants' migration

from environmentally hazardous facilities considering surface water flow and infiltration P.K. Shalkevich, S.P. Kundas, I.A. Gishkeluk

Введение

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

Ыатематическая модель распространения загрязняющих веществ от экологически опасного объекта с учетом поверхностного стока и инфильтрации

Разработанная математическая модель базируется на следующих уравнениях [1-4]. 1. Уравнение распространения загрязняющих веществ с поверхностным стоком:

СС)-АГ Е СС )-«/

& X ' У ' дх | х дх) ду

Е дСл

Еу ду

+ ^С = 0,

где С - концентрация загрязняющего вещества; и и V - скорость течения воды в направлении х и у соответственно; Ех и Еу - коэффициенты дисперсии в направлении х и у соответственно; ^ - скорость переноса массы вещества от источника.

2. Система уравнений, описывающая поверхностный сток:

дп д[(Л + л)и] д[(Л + л)и]

дt

дх

- + -

дх

= 0,

ди ттди тгди дп е др h + п PakW2cosю —+и—+V— = /V - е—1 ——^—з-

дt дх ду дх р дх 2 h + п

дV rтдV тrдV дп е др h + п sгnю —+и—+V— = /и - е—1 ——^—3--

дt дх ду дх р ду 2 h + п

иУи2 + V2 г_ ^и2 + V2

(д 2V

(л + п) с2

8

+ —

р

дх2

д V

+ -

а 2и ^

V

Г Я 2

ду2

2

- + -

д 2и

дх2 ду2

,

где х и у - горизонтальные координаты; t - время, и и V - скорость течения воды в направлении х и у соответственно; Л - глубина потока; п - высота поверхности воды; /- сила Кориолиса; е - ускорение свободного падения; р - плотность воды; ра - плотность воздуха; k - эмпирический коэффициент ветровой нагрузки; Жв - скорость ветра; ф - направление ветра; С - коэффициент Шези; е - коэффициент турбулентной вязкости.

3. Уравнение движения влаги в зоне аэрации с учетом инфильтрации:

дw дТ дрщ = у

дТ дt

дР1щ &

К^ СР, УТ * % дТ .

+

(

+У

р щ

К0 УР +р К0 Ку дР

У ГПс^ ГУ

К0 К1га

\

рЩ

Пщ П Р Пна

где К0 - коэффициент проницаемости пористой среды; Кщ (@гс) и К- (©^ - коэффициенты относительной фазовой проницаемости жидкости и водяного пара; ЦцС1 и ^ - коэффициенты динамической вязкости жидкости и водяного пара; е - ускорение свободного падения; г -координата в вертикальном направлении.

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

е,. аС+сс0»с

1гс

дt да

дt

+рь да = V(еftсDVC - иС)+Е ,

рь = рьР(КС - а) =

где 0йс - объемное влагосодержание грунта; С и а - концентрация загрязняющего вещества в растворе и адсорбированном состоянии; t - время; рь - плотность твердой фазы грунта; u - вектор скорости движения влаги в грунте; Е - функция, характеризующая интенсивность изменения концентрации загрязняющего вещества за счет процессов химического и радиоактивного превращения, поглощения вещества растениями и т. п.; К^ - коэффициент распределения, показывающий отношение концентрации адсорбированного вещества к концентрации вещества в растворе.

5. Уравнение движения грунтовых вод:

дН д К

ц— = —I Кх

дt дх \

дН

дх

д + — ду

(

К,

дН

ду

\

д

дН

+сг I К ^ Iа (х, у,t),

где ц - коэффициент водоотдачи (недостаток насыщения); Кх, Ку и К2 - коэффициенты фильтрации в направлении осей х, у и z соответственно; Н(х_у,2,0 - гидравлический напор; G (х,у,2,Л) - известная интенсивность пополнения (отъема) грунтовых вод.

6. Уравнение распространения загрязняющих веществ в грунтовых водах:

8КС )

8

8л

8х

8С \ 8

8С

= —\ еипу.— \ +—\ еиВу—

,гд у

8х) 8у\ 8

1ц г

8у

8

82

8С

+—\ еиПА.— \-

1ц у

82

К К )-8у К К у )-82 К К )+Q (х. у, 2,л).

где С - искомая концентрация загрязняющего вещества; Б, - коэффициент гидродинамической дисперсии; их, иу, - скорость движения грунтовых вод в направлении осей х, у и z соответственно; Q (х,у,2,Л) - известная интенсивность источников загрязнений.

7. Уравнения Сен-Венана, описывающие течение воды по поверхности грунта [5-8]:

8h 8uh 8vh

— +-+-= 0,

8л 8х 8у

8Ы 8 -+ —

8л 8х

+1 gh2\ + — ^иу) = hfr - gh—, 2 ) 8у х 8х

8Ну 8 8 ( 2 1

-+ — (Нич) + — \ Иу +— £

8л 8х 8у ^ 2

8C

= ¥у - gh 8ь,

8у

8л

:у(е11ч DVC )-v(uC).

8. Уравнение влагопереноса:

8w 8T 8w 8Р.

---^ = V

(

8Т 8л (

-V

8РЩ 8л

Ко Ку 8Р

VT

\

Ри

К0 Кщ_ ^^ Ко К

Чу

8Р

Чг

гд

Чу 8р

8Т

^ргд -РН,

Ид

Ко Кщ

Чйд

Рп^Б

9. Уравнение конвективной диффузии загрязняющих веществ с учетом кинетики сорбции:

е,. 8С+с ^

Ид

8л

8а

8л

+Рь § ^(е^с иС),

Рь — = РРь (КаС - а).

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

Численное решение дифференциальных уравнений, описывающих перенос

загрязняющих веществ

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

1. Уравнение конвективной диффузии [9-12]. Процедура решения с помощью метода конечных элементов уравнений конвективной диффузии заключается в следующем. Концентрация загрязняющих веществ (С) аппроксимируется полиномом

2

Р

V

С « С = £ С^, .

.=1

Применяя метод взвешенных невязок, уравнение (11) с граничными условиями записывается в виде

|N дС ап-| му^ - аЖ;) С;. дп+1 N..

п д п г

Интегрируя по частям член, характеризующий конвективный перенос и дисперсию, и принимая N. = - N., выражение (15) перепишется следующим образом:

| Nl дС an-JvNI

п д п

С учетом аппроксимаций

дС.....

т

п+у С"+1 - С"

К ,

Сп+у =(1 -у) С" +уСп+1 уравнение (16) запишется в виде

-- [(1 - у)е", + уе".+1 ]С;+1дп - [(1 - у)е", + уе"..+1 ]С"дп + (1 - у) -

^п п —> п

-(1 - У)-7 JVN. (е%DVNj -uNj )Спдп -уJvNI.БVNj -uNj .П = 0.

п п

Подставляя в систему уравнений (18) выражения (14) и применяя квадратурную формулу Гаусса-Лобатто для вычисления подынтегральных выражений, получим систему линейных алгебраических уравнений относительно Сп+1. Решение этой системы с использованием метода UMFPACK [1] дает значения концентрации загрязняющих веществ в момент времени п+1.

2. Уравнение Сен-Венана. Уровень воды (И) и компоненты вектора скорости течения (и и V) аппроксимируются полиномами:

МММ

И «И = £; и «и = £; V«V? = £vJNJ .

.=1 .=1 .=1

Применяя метод взвешенных невязок, уравнения (10) записываются в виде:

, дИ г дN дN. |^дN дN.

— дп + дп + дп = 0;

п дt п дх дх п дх дх

fдNI. Л 2 1 ,2^ ^ fдNI. дNi дЬ

глт дИ? fдNI дNi fдNI Г, 2 1 ,2 ^ ^ , ^ , дЬ

Ш..-дп +1—^—^иИЩ +1—^I h;v2 + -ghj I дп = И^ -^И,—.

п . дt п дх дх 1 3 3 п дУ I . . 2 . | у 3 ду

С учетом аппроксимаций (21) и численного интегрирования по квадратурной формуле Гаусса-Лобатто уравнения (20) сводятся к системе линейных алгебраических уравнений, решение которой дает значения уровня воды и компонент вектора скорости течения (и и V) в момент времени п+1. Трехмерное решение задачи достигается за счет специального преобразования членов уравнения N [9] к соответствующему виду и решения соответствующей для трехмерной задачи матрицы.

., Ип+у =(1 -у)Ип +уИп+1,

, ип+у=

Vв+у = (1 -у) Vй + ^п+\

дИ п+у И+1 -

дt AtИ

ди п+у и -1

57 К

п+у п+1 V - V

дt -п

п

Численное решение дифференциальных уравнений, описывающих поступления загрязняющих веществ в грунтовые воды с поверхности грунта за счет инфильтрации

атмосферных осадков

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

аппроксимируются полиномом:

м ~ м

= У w.N., С «С = У С.Ы.

1 1' 1 1

w « м>

а « а

= У а.Ы..

1 1

с г

1

,_1 ,_1 ,_1

Применяя метод взвешенных невязок в постановке Галеркина, уравнения (12) и (13) раничными условиями записываются в виде:

?+N г Р

+

I

+

| [ N ук^м^, + N.. ук^м^тр^ ] ап

п

(, + Чщ )аГ = 0;

аС

N0^ ал ап -Г

0 N1Pь аа J п

а?

(0и,,DVN, - иУК,) 0 0 0

С,

ап +

00

-щрькN ^^к^,.

Здесь

С,

ап = 0.

С _ащ. С . к ккр. к _

^Ьу ~ ^т ' ^ ^ ' Ьу ~ Г V _ , Л-,,, —

ат

ар.,

к0 ^ ^

-; к,.,

р к^ки+р к0 к ар

кг., Ну

Чв,

Чу ар

г.,

' ЛЧ Чу ат Ч 1щ

Искомые величины (давление, концентрация загрязняющего вещества в растворенном и адсорбированном виде) и их производные по времени в момент времени Лп + у аппроксимируются уравнениями:

ди '

а7

ип+1 - ип

д?п

_(1 -у) ип +уип

Здесь переменная и - искомая величина (давление, концентрация загрязняющего вещества в растворенном и адсорбированном виде). Тогда уравнение (23) запишется в виде

N. [(1 -у)0П, +У0Пц ]

д? •*

+(1 -у)дт I

пп

-(1 -у)дт I

N

( 0'+1 - 0П ^ 1щ 1щ

д?п

N1

Сп

п+1

ап --

д- I

д?п I

N [(1 -у)0п, +у0 0

VN (01,DVN, - ) 0 00

п

ап + у

N

( 0п+1 - 0п ^

1щ 1щ

дп

N,

N¡Pь N1Pь

0

ап +

Сп

п+1

ап -

Сп

ап - у Г

а". Л п

(01щ DVN, - ) 0 00

Сп

п+1

ап _ 0.

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

а

а

п

а

0

0

п

а

пп

ф

можно получить следующее приближение к искомому решению задачи:

[ рГ ] *+1 =[ РП+1 ] *+[>< ]* •

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

Ф( рС ) *

др.;1

дри

I I

W P+1 )|

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

С каждой k-й итерацией метода Ньютона - Рафсона искомые величины на временном шаге n+1 приближаются к значениям, которые являются решением системы уравнений (25).

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

Отметим, что в систему уравнений (25) входят подынтегральные выражения. Для их вычисления можно воспользоваться численным интегрированием, при котором определяемый интеграл заменяется квадратурной суммой. Использование квадратурных формул позволит, при программной реализации численного метода, оперативно изменять тип конечного элемента и степень многочленов, используемых в качестве базисных функций. Как показали предварительные расчеты, при применении квадратуры Гаусса-Лобатто получаемое решение задачи поступления загрязняющих веществ в грунтовые воды с поверхности грунта за счет инфильтрации атмосферных осадков обладает большей устойчивостью и точностью. Как и в предыдущем разделе, трехмерное решение задачи достигается за счет специального преобразования членов уравнения N [9] к соответствующему виду и решения соответствующей для трехмерной задачи матрицы.

Компьютерная реализация разработанной математической модели

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

Для проведения расчетов в системе MATLAB разработан механизм экспорта расчетной сетки из программы COMSOL. В программе COMSOL в Model Navigator выбирается вид компьютерной модели PDE, Coefficient Form. Перейдя в программе COMSOL в режим создания геометрии Draw Mode, с использованием инструментов программы для трехмерных геометрических построений создается геометрия расчетной области (рис. 1, а).

'38?-

л т »»> »* 7*«» ю» t

5!

а б

Рис. 1. Создание трехмерной геометрии расчетной области (а) и поверхностей, на которых задаются

граничные условия

В режиме работы с граничными условиями Boundary Settings, в программе COMSOL определяются поверхности, на которых будут заданы граничные условия (рис. 1, б). Далее

созданная трехмерная геометрия разбивается на конечно-элементную расчетную сетку в режиме программы Mesh Mode (рис. 2, а). Полученная расчетная сетка экспортируется в файл с использованием команды программы COMSOL File->Export->Mesh To File.

Рис. 2. Создание тетраидальной расчетной сетки (а) и визуализация концентрации загрязняющих

веществ в сечениях расчетной области

Затем в системе COMSOL в меню программы Postprocessing с использованием различных инструментов (сечения, построение графиков и др.) осуществляется представление и анализ результатов расчетов (рис. 2, б).

Верификация разработанных моделей с использованием экспериментальных данных, приведенных в [2], показала, что погрешность моделирования не превышает 5 %.

Заключение

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

Разработаны численные методы решения дифференциальных уравнений, описывающих поступления загрязняющих веществ в грунтовые воды с поверхности грунта за счет инфильтрации осадков и поверхностного стока. Особенностью численного метода является использование квадратурной формулы Гаусса-Лобатто для вычисления интегральных выражений, получаемых при применении метода взвешенных невязок для аппроксимации искомых величин. С использованием программных комплексов MATLAB и COMSOL осуществлена компьютерная реализация разработанной математической модели. Верификация разработанных моделей с использованием экспериментальных данных, приведенных в [2], показала, что погрешность моделирования не превышает 5 %.

Список литературы

1. Gishkeluk I.A., Kundas S.P., Hilko O.S. Numerical simulation of nonisothermal water and contaminant transport in natural disperse media // Protection and Restoration of the Environment XI: Proceedings of International Conference. Thessaloniki, 3-6 July, 2012. P. 721-729.

2. Компьютерное моделирование миграции загрязняющих веществ в природных дисперсных средах / С.П. Кундас [и др.]. Минск: МГЭУ им. А.Д. Сахарова, 2011. 212 с.

3. Гринчик Н.Н., Гишкелюк И.А., Кундас С.П. Моделирование тепломассопереноса и поверхностных явлений в капиллярно-пористых средах на основе уравнений двухфазной фильтрации и изотерм сорбции // Сб. научных статей «Современная наука 2010». 2011. № 1. С. 56-60.

4. Гишкелюк И.А., Кундас С.П. Численное моделирование неизотермического петреноса влаги и растворимых веществ в природных дисперсных средах // Экологический вестник. 2011. № 1 (15). С. 97-104.

5. Кивва С.Л., Железняк М.И. Численное моделирование двумерного открытого потока с подвижными границами: расчеты стока на водосборе и наката волн цунами на берег // Труды Международной конференции RDAMM-2001. 2001. Спец. выпуск. Т. 6. Ч. 2. C. 343-350.

6. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. М.: Физматлит, 2001. 732 с.

7. Randall D.A. The shallow water equations // Water. 2006. Vol 41. Is. 3. P. 695-719.

8. Куликовский А.Г., Погорелов Н.В., Семенов А.Ю. Математические вопросы численного решения гиперболических систем уравнений. М. : Физматлит, 2001. 608 с.

9. Зенкевич О., Морган К. Конечные элементы и аппроксимация. М., 1986. 318 с.

10. Zienkievich O.C., Taylor R.L. The Finite Element Method: The Basis. Oxford., 2000. 708 p.

11. Самарский А.А. Теория разностных схем. М.: Наука, 1983. 616 с.

12. Сегерлинд Л. Применение метода конечных элементов. М.: Мир, 1979. 392 с.

References

1. Gishkeluk I.A., Kundas S.P., Hilko O.S. Numerical simulation of nonisothermal water and contaminant transport in natural disperse media // Protection and Restoration of the Environment XI: Proceedings of International Conference. Thessaloniki, 3-6 July, 2012. P. 721-729.

2. Komp'juternoe modelirovanie migracii zagrjaznjajushhih veshhestv v prirodnyh dispersnyh sredah / S.P. Kundas [i dr.]. Minsk: MGJeU im. A.D. Saharova, 2011. 212 s.

3. Grinchik N.N., Gishkeljuk I.A., Kundas S.P. Modelirovanie teplomassoperenosa i poverhnostnyh javlenij v kapilljarno-poristyh sredah na osnove uravnenij dvuhfaznoj fil'tracii i izoterm sorbcii // Sb. nauchnyh statej «Sovremennaja nauka 2010». 2011. № 1. S. 56-60.

4. Gishkeljuk I.A., Kundas S.P. Chislennoe modelirovanie neizotermicheskogo petrenosa vlagi i rastvorimyh veshhestv v prirodnyh dispersnyh sredah // Jekologicheskij vestnik. 2011. № 1

5. Kivva S.L., Zheleznjak M.I. Chislennoe modelirovanie dvumernogo otkrytogo potoka s podvizhnymi granicami: raschety stoka na vodosbore i nakata voln cunami na bereg // Trudy Mezhdunarodnoj konferencii RDAMM-2001. 2001. Spec. vypusk. T. 6. Ch. 2. C. 343-350.

6. Landau L.D., Lifshic E.M. Teoreticheskaja fizika. M.: Fizmatlit, 2001. 732 s.

7. Randall D.A. The shallow water equations // Water. 2006. Vol. 41. Is. 3. P. 695-719.

8. Kulikovskij A.G., Pogorelov N.V., Semenov A.Ju. Matematicheskie voprosy chislennogo reshenija giperbolicheskih sistem uravnenij. M. : Fizmatlit, 2001. 608 s.

9. Zenkevich O., Morgan K. Konechnye jelementy i approksimacija. M., 1986. 318 s.

10. Zienkievich O.C., Taylor R.L. The Finite Element Method: The Basis. Oxford., 2000. 708 p.

11. Samarskij A.A. Teorija raznostnyh shem. M.: Nauka, 1983, 616 s.

12. Segerlind L. Primenenie metoda konechnyh jelementov. M.: Mir, 1979. 392 s.

Сведения об авторах

Шалькевич П.К., старший преподаватель Международного государственного экологического института имени А.Д. Сахарова Белорусского государственного университета.

Information about the authors

Shalkevich P.K., senior lecturer of International environmental institute named after A.D. Sakharov of Belarusian state university.

Кундас С.П., д.т.н., профессор, профессор кафедры «Теплогазоснабжение и вентиляция» Белорусского национального технического университета.

Гишкелюк И.А., к.т.н., доцент, доцент кафедры экологических информационных систем Международного государственного экологического института имени А.Д. Сахарова Белорусского государственного университета.

Kundas S.P., Dr.Sci, professor, professor of the department of heat and gas supply and air-conditioning of Belarusian national technical university.

Gishkeluk I.A., PhD, associate professor, associate professor of the department of ecological information systems of International environmental institute named after A.D. Sakharov of Belarusian state university.

Адрес для корреспонденции

220070, Республика Беларусь, г. Минск, ул. Долгобродская, д. 23, Международный государственный экологический институт имени А.Д. Сахарова Белорусского государственного университета тел. +375-17-230-68-97; e-mail: pavel.shalkevich@gmail.com Шалькевич Павел Константинович

Address for correspondence

220070, Republic of Belarus, Minsk, Dolgobrodskaya, 26 International environmental institute named after A.D. Sakharov of Belarusian state university tel. +375-17-230-68-97; e-mail: pavel.shalkevich@gmail.com Shalkevich Pavel Konstantsinovich