Научный журнал
Успехи современного естествознания
ISSN 1681-7494
"Перечень" ВАК
ИФ РИНЦ = 0,823

ЭКОЛОГО-ГИДРОФИЗИЧЕСКОЕ ОБОСНОВАНИЕ ВЛИЯНИЯ КОЭФФИЦИЕНТА ВЕРТИКАЛЬНОГО ТУРБУЛЕНТНОГО ОБМЕНА НА СОДЕРЖАНИЕ РАСТВОРЕННОГО КИСЛОРОДА В ПРИДОННОМ СЛОЕ МЕЛКОВОДНОГО ВОДОЕМА

Никитина А.В. 1 Сухинова Т.Г. 2 Проценко С.В. 3 Семенякина А.А. 1 Бедная Т.А. 3
1 Южный федеральный университет
2 Южный научный центр Российской академии наук
3 Донской государственный технический университет
В статье представлены результаты гидрохимических наблюдений, полученных в ходе экспедиции на исследовательском судне «Денеб» в восточной части Азовского моря и Таганрогской бухты в июле 2017 г. В ходе экологических и гидрофизических исследований были обнаружены аномально низкие уровни растворенного кислорода в нижнем слое Азовского моря. Загрязняющие вещества скапливаются в донных отложениях водоемов. В теплое время года, при благоприятном температурном режиме донные осадки прогреваются, иловые микроорганизмы освобождают вещества, вызывающие выделение сероводорода. Результатом распада существенного количества органических веществ за относительно небольшой период времени является значимое снижение концентрации кислорода, растворенного в воде, что может служить причинами отрицательных изменений водной экосистемы: замора рыбы, ухудшения качества воды. Экспериментальным путем получены распределения коэффициента вертикального турбулентного обмена. Проблемы обеспечения стабильного функционирования водных экосистем, предотвращения уменьшения биологического разнообразия, предупреждения негативных последствий антропогенного воздействия, обеспечения безопасного водопользования в процессе социально-экономического развития территории приобретают в настоящее время особую актуальность. Проведен численный эксперимент по определению коэффициента вертикального турбулентного обмена для некоторых точек Азовского моря. Получено соответствие результатов расчетов на основе статистических данных о поле скоростей водного потока и подсеточной модели Смагоринского с данными натурных измерений. Изучено влияние изменения значений коэффициента вертикального турбулентного обмена на содержание растворенного кислорода в придонном слое мелководного водоема. Выявлена прямая зависимость между значениями коэффициента вертикального турбулентного обмена и концентрацией растворенного кислорода в нижнем слое.
коэффициент вертикального турбулентного обмена
турбулентный поток
концентрация растворенного кислорода
1. Алексеенко Е.В. Сравнительный анализ классических и неклассических моделей гидродинамики водоемов с турбулентным обменом / Е.В. Алексеенко, Б.В. Сидоренко, О.В. Колгунова, А.Е. Чистяков // Известия ЮФУ. Технические науки. – 2009. – № 8 (97). – С. 6–18.
2. Sukhinov A.I. Mathematical modeling of sediment transport in the coastal zone of shallow reservoirs / A.I. Sukhinov, A.E. Chistyakov, E.A. Protsenko // Mathematical Models and Computer Simulations. – 2014. – Т. 6, № 4. – С. 351–363.
3. Дегтярева Е.Е. Программная реализация трехмерной математической модели транспорта взвеси в мелководных акваториях / Е.Е. Дегтярева, Е.А. Проценко, А.Е. Чистяков // Инженерный вестник Дона. – 2012. – Т. 23, № 4–2 (23). – С. 30.
4. Никитина А.В. Дифференциально-игровая модель предотвращения заморов в мелководных водоемах / А.В. Никитина, М.В. Пучкин, И.С. Семенов, А.И. Сухинов, Г.А. Угольницкий, А.Б. Усов, А.Е. Чистяков // Управление большими системами: сборник трудов. – 2015. – № 55. – С. 343–361.
5. Проценко Е.А. Двумерная конечно-разностная модель формирования наносов в прибрежной зоне водоема и ее программная реализация / Е.А. Проценко // Инженерный вестник Дона. – 2010. – Т. 13, № 3. – С. 23–31.
6. Проценко Е.А. Модель и алгоритмы решения задачи о транспорте наносов/ Е.А. Проценко // Известия ЮФУ. Технические науки. – 2009. – № 8 (97). – С. 71–75.
7. Protsenko S. Mathematical modeling of wave processes and transport of bottom materials in coastal water areas taking into account coastal structures / S. Protsenko, T. Sukhinova // MATEC Web Conf. Volume 132, 2017 XIII International Scientific-Technical Conference «Dynamic of Technical Systems» (DTS-2017) DOI. URL: https://doi.org/10.1051/matecconf/201713204002.
8. Белоцерковский О.М. Турбулентность: новые подходы / О.М. Белоцерковский, А.М. Опарин, В.М. Чечеткин. – М.: Наука. – 2003. – 286 с.
9. Smagorinsky J. General circulation experiments with the primitive equations: 1. The basic experiment / J. Smagorinsky // Monthly Weather Review. – 1963. – Vol. 91, № 3. – P. 99–164.

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

Происходящие в настоящее время увеличение объёмов промышленного производства, сброс производственных отходов и загрязняющих веществ в водоемы резко увеличивают нагрузку на естественные экосистемы. Следствием растущего антропогенного воздействия на природу, аномальных климатических изменений, естественного старения водного объекта являются риски нарушения естественной жизнедеятельности в экосистеме водного объекта, возникновения необратимых отрицательных изменений водной экосистемы, процессов эвтрофикации водоемов – перенасыщения водоема биогенами, что влечет за собой быстрый рост фитопланктона. Экосистема мелководного водоема имеет сложную структуру, в рамках которой происходят взаимодействия органических и неорганических веществ с живыми организмами. Среди всех гидрохимических параметров, влияющих на биологическую продуктивность Азовского моря, особое место отводят растворенному кислороду – источнику дыхания гидробионтов и фактору, определяющему полноту и скорость минерализации органических веществ в водоеме. Особенностью кислородного режима Азовского моря является его нестабильность, которая определяется рядом постоянных и сезонных факторов: поглощением атмосферного кислорода поверхностным слоем воды, неравномерным притоком речных и черноморских вод, развитием и затуханием жизнедеятельности организмов, циркуляционными процессами, определяющими возможности проникновения кислорода в придонные слои, и другими процессами [1–2].

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

Летом в условиях плотностной стратификации возникает усиленное снижение концентрации кислорода в придонных горизонтах Азовского моря (в отдельных районах до 60 %). Ведущими факторами формирования кислородного режима в поверхностных водах являются главным образом продукционные процессы, в придонных – биохимическое разложение донных отложений [4–5]. Принимая во внимание процессы вертикальной турбулентной диффузии в стратифицированных водоемах, кислород в процессе инвазии попадает в нижние слои. Цель настоящей работы состоит в проверке влияния гидродинамических условий на изменение концентрации растворенного кислорода в придонном слое мелководного водоема.

Описание измерительной базы. Маршрут экспедиции

Экспедиция с целью проведения научно-исследовательских работ в акватории Азовского моря «Денеб» (далее НИС «Денеб») Южного научного центра РАН (далее ЮНЦ РАН). Исследования проводились на 17 станциях в восточной части Азовского моря (рис. 1).

nik1a.tif nik1b.tif

а) б)

Рис. 1. Научно-исследовательские работы в Азовском море (июль 2017 г.): а) маршрут экспедиции; б) НИС «Денеб»

Оборудование НИС «Денеб»: гидрологический CDT-зонд SEACAT SBE19; регистратор течений RCM 9LW; узколучевой параметрический профилограф SES-2000 light (Innomar Technologie GmbH); гравитационная грунтовая трубка 13,540 B с возможностью установления поршневой системы (Piston/Gravity corer Model 13,540 B); автоматический проточный анализатор San++ с сэмплером SA1100, вмещающим 2×50 позиций для образцов; глубоководный пробоотборный комплекс карусельного типа; оборудование НИС «Денеб» для гидробиологических и литологических исследований: батометры Молчанова и Нискина для отбора проб воды; дночерпатели Петерсена и Ван-Виина для отбора донных отложений; бентосная драга; планктонные сети Апштейна и Джеди для отбора проб планктона; икорная сеть для отбора проб ихтиопланктона; набор сетей и драг для проведения ихтиологических исследований; грунтовая прямоточная трубка с возможностью взятия колонки осадков 2–2,5 м.

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

Оборудование НИС «Денеб» для ихтиологических исследований: донный бимтрал; пелагический распорный трал (28 м по горизонтали, 8 м по вертикали); гидрофизический зонд SBE19 plus; датчик растворенного кислорода SBE43; ПК (персональный компьютер) с программным обеспечением для подключения гидрофизических зондов; батометры 3 л и 5 л; комплектное оборудование для измерения концентрации кислорода методом Винклера; профилограф (ADCP) WHS 600.

Содержание выполненных экспедиционных работ

Экспедиционные работы проводились с целью получения систематической информации влияния химических и физических свойств воды на экологическое состояние акватории Азовского моря. Программа наблюдений за качеством морских вод состояла в следующем. Была найдена концентрация нефтяных углеводородов, мг/дм3 (мг/л); растворенного кислорода, мг/дм3 (мг/л, %); водородного показателя рН; хлорированных углеводородов, в том числе пестицидов, мкг/дм3 (мкг/л); тяжелых металлов, мкг/дм3 (мкг/л); фенолов, мкг/дм3 (мкг/л); других веществ, мкг/дм3 (мкг/л). Были выявлены показатели содержания биогенных веществ, свойственных для района проведения экспедиции: нитритного азота, мкг/дм3 (мкг/л); кремния, мкг/дм3 (мкг/л); солености воды, %; температуры воды и воздуха, °С; скорости и направления ветра, м/с; прозрачности воды (единицы цветности); волнения моря (баллы); проведение визуальных наблюдений за состоянием поверхности морского водного объекта.

Измерение пульсаций вектора скоростей водного потока

При моделировании полей течений водного потока мелководных водоемов возникает проблема параметризации коэффициента турбулентного обмена. От турбулентного перемешивания во многом зависит способность водоема к самоочищению, турбулентный обмен влияет на транспорт наносов и взвешенных веществ и на интенсивность деформации донной поверхности [6–7].

В ходе экспедиции были получены данные о пульсациях скоростей в некоторых точках мелководных систем на основе зонда ADCP (Acoustic Doppler Current Profiler) WHS600 Sentinel. ADCP-профилограф был зафиксирован на свободной поверхности водоема, измерения проводились от свободной поверхности до дна. Работа профилографа основана на доплеровском эффекте. Из значений лучевых скоростей вычисляют пульсации компонент скорости nik02.wmf.

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

Концентрация растворенного кислорода в Азовском море

№ станции

Координаты станции

Измеренные значения растворенного кислорода

1(49)

46 ° 39, 569' N

37 ° 20, 743' E

глубина, м

0

4

11,5

О2, мг/л

9,02

8,53

5,32

2(48’)

46 °29,725' N 36 °59,786' E

глубина, м

0

6

12

О2, мг/л

12,69

13,03

0,66

3(44’)

46 ° 19, 693' N 36 ° 59, 683' E

глубина, м

0

6

12

О2, мг/л

5,75

4,94

0,61

12(42)

45 ° 59, 906' N 37 ° 20, 128' E

глубина, м

0

5,5

11

О2, мг/л

9,21

9,87

1,41

Коэффициент вертикальной турбулентной вязкости

Проведенные натурные эксперименты, в том числе измерение пульсации скоростей водного потока, позволяют при построении математических моделей учитывать коэффициент вертикального турбулентного обмена, неоднородный по глубине. Результаты, наиболее близкие к экспериментальным натурным данным, получены с помощью параметризации О.М. Белоцерковского [8].

Если вертикальное разрешение расчетной сетки позволяет воспроизвести все механизмы моделируемого процесса до масштабов вязкой диссипации очень мелких вихрей, то для расчета турбулентности используют уравнения Навье – Стокса и находят путем численного моделирования [3, 8]. Скорость диссипации может быть представлена как

nik03.wmf

где nik04.wmf – осредненный тензор скоростей деформации: nik05.wmf

Выражение для турбулентной вязкости имеет вид

nik06.wmf

где Cs – константа, значение которой зависит от однородной изотропной турбулентности. Значение константы Cs должно гарантировать соответствие результатов моделирования с экспериментальными измерениями. Значения константы Cs в ряде случаев изменяются от 0,1 до 0,2, при течении на канале до случая изотропной турбулентности соответственно.

Если размер сетки достаточно мал, то nik07.wmf равна средней скорости деформации в параметризации Прандтля D, что позволяет перейти к аналогу Прандтля в виде

nik08.wmf

пренебрегая градиентами nik09.wmf, nik10.wmf, nik11.wmf, имеем равенство

nik12.wmf.

Численный расчет коэффициента вертикальной турбулентной вязкости

Были выполнены численные эксперименты на основе описанных выше подходов для всех точек, nik13.wmf в которых проводились натурные измерения в ходе экспедиционных работ (рис. 2).

nik2.tif

Рис. 2. Коэффициент вертикального турбулентного обмена в точках 1, 2, 3, 12, по горизонтали – значения в м2/с, по вертикали – уровень воды от свободной поверхности ко дну

Было изучено влияние изменения значений коэффициента вертикального турбулентного обмена на содержание растворенного кислорода в придонном слое мелководного водоема [1, 3]. Анализ результатов моделирования процессов гидрофизики в Азовском море показал, что наибольшему значению коэффициента вертикального турбулентного обмена соответствует наибольшая концентрация растворенного кислорода в придонном слое Азовского моря (точка 12), наименьшему значению концентрации растворенного кислорода (точка 3) соответствует наименьшее значение. На глубине около и более трех метров значения данного коэффициента близки к 0, что означает пониженный турбулентный обмен по вертикали в указанной области и объясняет отсутствие кислорода в придонном слое центрально-восточной части Азовского моря.

Выводы

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

Получено соответствие результатов расчетов на основе статистических данных о поле скоростей водного потока и подсеточной модели Смагоринского [9] с данными натурных измерений. Изучено влияние изменения значений коэффициента вертикального турбулентного обмена на содержание растворенного кислорода в придонном слое мелководного водоема.

Работа выполнена при финансовой поддержке Программы фундаментальных исследований Президиума РАН I.33П, проект № 0256-2015-0078.


Библиографическая ссылка

Никитина А.В., Сухинова Т.Г., Проценко С.В., Семенякина А.А., Бедная Т.А. ЭКОЛОГО-ГИДРОФИЗИЧЕСКОЕ ОБОСНОВАНИЕ ВЛИЯНИЯ КОЭФФИЦИЕНТА ВЕРТИКАЛЬНОГО ТУРБУЛЕНТНОГО ОБМЕНА НА СОДЕРЖАНИЕ РАСТВОРЕННОГО КИСЛОРОДА В ПРИДОННОМ СЛОЕ МЕЛКОВОДНОГО ВОДОЕМА // Успехи современного естествознания. – 2018. – № 1. – С. 115-119;
URL: http://natural-sciences.ru/ru/article/view?id=36660 (дата обращения: 16.06.2021).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074