Имрв «Поток» 1. 0 - имитационная математической модель для расчёта распространения и седиментации технологических наносов в водотоках при определении вреда водным биологическим ресурсам

Денисенко Олег Сергеевич; Живчиков Виталий Георгиевич

УДК 004.94

03.00.00 Биологические науки

ИМРВ «ПОТОК» 1.0 - ИМИТАЦИОННАЯ МАТЕМАТИЧЕСКОЙ МОДЕЛЬ ДЛЯ РАСЧЁТА РАСПРОСТРАНЕНИЯ И СЕДИМЕНТАЦИИ ТЕХНОЛОГИЧЕСКИХ НАНОСОВ В ВОДОТОКАХ ПРИ ОПРЕДЕЛЕНИИ ВРЕДА ВОДНЫМ БИОЛОГИЧЕСКИМ РЕСУРСАМ

Денисенко Олег Сергеевич к.б.н, БРШ-код: 3403-3102 ФГБНУ «Азовский научно исследовательский институт рыбного хозяйства», Краснодарское отделение, г. Краснодар, Россия denisenko_o_s@azniirkh. т

Живчиков Виталий Георгиевич к.с.-х.н., БРШ-код: 6132-6032 ООО «Азово-Черноморский научный центр рыбохозяйственных исследований», г. Краснодар, Россия, [email protected]

В статье приведено описание и алгоритм работы имитационной математической модели ИМРВ «Поток» 1.0, предназначенной для расчёта параметров распространения и седиментации технологических наносов в водотоках. Актуальность имитационной модели обусловлена необходимостью определения интегральных показателей (объём воды, протёкшей через области шлейфа с заданной концентрацией взвеси; средний слой наилка; удельная плотность заиления дна в расчётной зоне) при расчете вреда водным биологическим ресурсам от осуществления различных видов хозяйственной деятельности в соответствии с методикой исчисления размера вреда, причинённого водным биологическим ресурсам, утвержденной приказом Федерального агентства по рыболовству №1166 от 25.11.2011 г. Использование имитационной математической модели ИМРВ «Поток» 1.0 позволит повысить достоверность расчетов при определении вреда водным биологическим ресурсам, а также привести деятельность расчёта параметров распространения и седиментации технологических наносов в водотоках в соответствие с действующими нормативными правовыми актами Федерального агентства по рыболовству

Ключевые слова: ИМИТАЦИОННОЕ МОДЕЛИРОВАНИЕ, ИМРВ «ПОТОК» 1.0, РАСЧЕТ ВРЕДА, ВОДНЫЕ БИОРЕСУРСЫ, ИХТИОФАУНА, ШЛЕЙФ ДОПОЛНИТЕЛЬНОЙ МУТНОСТИ, ИНТЕГРАЛЬНЫЕ ПОКАЗАТЕЛИ

UDC 004.94 Biological sciences

SMDS "POTOK" 1.0 - A SIMULATION MATHEMATICAL MODEL FOR CALCULATION OF DISTRIBUTION AND SEDIMENTATION OF TECHNOLOGICAL ALLUVIUMS IN WATERCOURSES FOR DETERMINATION OF HARM TO WATER BIOLOGICAL RESOURCES

Denisenko Oleg Sergeyevich Сand. Biol. Sci., SPIN-code: 3403-3102 FSBSI Azov Sea research fisheries Institute, Krasnodar branch, Krasnodar, Russia denisenko_o_s@azniirkh. ru

Zhivchikov Vitaly Georgiyevich Cand. Agr. Sci., SPIN-code: 6210-0979 JSC Azovo-Chernomorsky Scientific Center of Fishery Researches, Krasnodar, Russia, zvg-vint@mail. ru

The article gives the description and the algorithm of the simulation model SMDS "Potok" 1.0, designed to calculate the parameters of distribution and sedimentation of technological alluviums in watercourses. The relevance of the simulation model is due to the need to determine the integral indicators (the volume of water flowing through the areas of the plume with a predetermined suspension concentration, the average layer of the silt, the specific density of the bottom silt in the calculated zone) when calculating the harm to water biological resources from the implementation of various economic activities in accordance with the methodology of calculation of damage caused to water biological resources approved by the order of the Federal Agency for Fisheries No. 1166 of 25.11.2011. Usage of the simulation model of SMDS "Potok" 1.0 will increase the reliability of calculations in determining the damage to water biological resources, as well as bring the calculations of parameters of distribution and sedimentation of technological alluviums in watercourses in accordance with the current regulatory legal acts of the Federal Agency for Fisheries

Keywords: SIMULATION MODELING, SMDS "POTOK" 1.0, CALCULATION OF HARM, WATER BIOLOGICAL RESOURCES, ICHTHYOFAUNA, PLUME OF ADDITIONAL TURBIDITY, INTEGRAL INDICATORS

Бок 10.21515/1990-4665-129-046

Введение

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

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

В соответствии с Методикой исчисления размера вреда, причиненного водным биологическим ресурсам, утвержденной Федеральным агентством по рыболовству приказом № 1166 от 25.11.2011 г. [12], сведения о количественной оценке интенсивности факторов воздействия, объёмов, площадей и зон их влияния выполняются посредством имитационного (математического) моделирования с

использованием специальных компьютерных программ.

Учитывая вышеизложенное, была разработана имитационная математическая модель ИМРВ «Поток» 1.0, предназначенная для расчёта параметров распространения и седиментации технологических наносов в водотоках. Федеральной службой по интеллектуальной собственности (РОСПАТЕНТ) имитационная математическая модель ИМРВ «Поток» 1.0 внесена в реестр программ для ЭВМ (свидетельство о государственной регистрации программы для ЭВМ № 2017613750, правообладатель ООО «Азово-Черноморский научный центр рыбохозяйственных исследований»). Сертификатом соответствия РОССЯи.04ЖИГ0.00081 (изготовитель ООО «Азово-Черноморский научный центр рыбохозяйственных исследований», срок действия сертификата - по 21.06.2019 г.) подтверждено соответствие требованиям следующих нормативных документов: ГОСТ 19179-73, ГОСТ 17.1.1.02-77, ГОСТ 826793, ГОСТ 12536-2014, ГОСТ Р ИСО/МЭК 9126-93, ГОСТ 17.1.3.07-82, ГОСТ Р ИСО/МЭК 12119-2000, ГОСТ 27065-86, ГОСТ Р ИСО/МЭК 159102002, ГОСТ 17.1.3.13-86.

Материал и методы

В качестве методической основы при разработке имитационной математической модели ИМРВ «Поток» 1.0 положены нормативные документы и литературные данные [1-15].

В модели ИМРВ «Поток» 1.0 указанные методики являются взаимодополняющими друг друга, при этом моделирование производится с расчётом всего необходимого набора площадей Р{ и поперечных размеров облака дополнительной взвеси В21.

Расчёт ведётся для величины превышения мутности над фоновыми значениями в предположении, что фоновая мутность соответствует гидравлическим характеристикам водотока и не приводит к заилению русла. При этом дополнительная (технологическая) мутность

рассматривается как загрязняющее вещество.

Модельный расчёт распространения и седиментации технологических наносов в водотоке посредством ИМРВ «Поток» 1.0, ведётся последовательно по этапам (рис. 1).

I ЭТАП: расчёт стартовой мутности в створе производства работ и расчёт времени воздействия пятна мутности на участок водотока в пределах его распространения

П ЭТАП: расчёт зон распространения частиц грунта разной крупности вниз по течению от створа производства земляных работ

Ш ЭТАП: расчёт массы грунта, оседающего по участкам реки ниже створа производства земляных работ, и расчёт изменения мутности вниз по водотоку

ИСХОДНЫЕ ДАННЫЕ: объём перемещаемого (извлекаемого) грунта; производительность землеройной техники; характеристика механического состава грунта; процентуноса грунта; морфометрические характеристики реки (ширина и глубина); скорость течения

ИСХОДНЫЕ ДАННЫЕ: гранулометрический состав грунта (по фракциям, в % массы); стандартная гидравлическая крупность частицрасчётного диаметра; температура воды; температурный поправочный коэффициент к значению гидравлической крупности; мор фо метр ические характеристики реки (ширина и глубина); скорость

ИСХОДНЫЕ ДАННЫЕ: результаты расчёта I и П этапов

IV ЭТАП: расчёт толщины слоя шщка. образовавшегося в результате производства земляных работ, и массы грунта, отложившегося на единицу площади дна водотока (плотности заиления)

V ЭТАП: расчёт площадей зон распространения частиц грунта разной крупности вниз по течению от створа производства земляных работ ирасчёт объёмов воды, загрязнённой взвесью частиц разной крупности

ВЫВОД РЕЗУЛЬТАТОВ

ИСХОДНЫЕ ДАННЫЕ: результаты расчёта Ш этапа; объём перемещаемого (извлекаемого) грунта;

производительность землеройной техники; процентуноса грунта; морфометрические характеристики реки (ширина и глубина); скорость течения

Рисунок 1 - Схема последовательности обработки данных с помощью

модели ИМРВ «Поток» 1.0

Результаты и обсуждение

Все естественные водотоки несут в своей толще и передвигают по дну твёрдые частицы грунта (или наносы).

Наносы включают в себя взвешенные наносы (наносы, переносимые водным потоком во взвешенном состоянии) и влекомые наносы (наносы, перемещаемые водным потоком в придонном слое и движущиеся путем скольжения, перекатывания или сальтации) [3].

Причиной взвешивания и переноса наносов во взвешенном состоянии является наличие в потоке вертикальных составляющих мгновенных или пульсационных скоростей, величина которых превосходит гидравлическую крупность «U» взвешиваемых наносов или, иными словами: превышает скорость выпадения частиц этих наносов в спокойной воде.

Передвижение наносов во влекомом состоянии (донных наносов) связано с наличием значительных придонных скоростей потока и вихревых образований воды, отделяющихся от дна.

Согласно теории потока, если в прямолинейном открытом русле, сложенном из однородных частиц, повышать величину расхода Q, то при определённой скорости частицы начнут выноситься потоком. Далее при постоянном расходе воды будет происходить седиментация частиц, и дно русла потока станет стабильным, сложенным при этом из частиц, находящихся в предельном состоянии равновесия, а течение потока станет спокойным, и в нём не будут содержаться наносы.

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

По истечении некоторого времени, расходы наносов в створах I-I и II-II станут одинаковыми и отметка дна русла в исследуемой зоне «а» станет стабильной.

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

Рисунок 2 - Схема распределения частиц наносов в русле водотока

Содержание (концентрация) взвешенных веществ - наносов в единице объёма смеси воды с наносами, выражаемое в весовых единицах -мг/л, г/м , называется мутностью воды [3, 6].

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

Большую часть расхода наносов Qн - около 90% составляют взвешенные наносы и только 10% составляют донные наносы, наиболее крупные, образующие гряды [15].

Таким образом, вымытые водотоком частицы грунта будут не только осаждаться под действием силы тяжести, но и одновременно смешиваться с водой водотока, а поэтому на некотором расстоянии от места сброса

концентрация взвешенных частиц (мутность), будет снижаться.

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

Допустимое содержание взвешенных частиц в водных объектах для разных типов водопользования представлено в таблице 1 [1, 2].

Таблица 1 - Требования к содержанию взвешенных веществ в воде водотоков и водоёмов различного назначения

КАТЕГОРИЯ ВОДОПОЛЬЗОВАНИЯ

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

Содержание взвешенных веществ (мг/л) не должно увеличиваться больше чем на:

0,25 0,75 0,25 0,75

Оценку воздействия проведения работ в водотоке на концентрацию технологических наносов (и мутность) следует выполнять с учётом типа водного объекта и технологии работ [11, 13, 15]. Тип водного объекта

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

- гидрометеорологические данные, включая расчётные расход и уровень воды, морфометрические характеристики русла, скорость течения и фоновые значения мутности;

- геологическая информация (осреднённый гранулометрический состав разрабатываемых грунтов и донных отложений);

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

Такие исходные данные полностью согласуются с требованиями с п. 34 гл. III «Методики исчисления размера вреда, причинённого водным биологическим ресурсам» [12].

Общетеоретические основы применяемой математической модели

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

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

Далее для расчётов используется двумерная (усреднённая по глубине) модель, полученная из трёхмерного уравнения конвекции-диффузии плотности взвеси ф.

В основе модели лежит выражение единичного адекватного и турбулентного переноса вещества через контрольную поверхность в водной среде г3 [г/(м •с)]:

Г = ^ - / 9п (1)

где: уп - компонента осредненной скорости течения, направленная по нормали п к контрольной поверхности; я - концентрация вещества в точке пересечения нормали с

контрольной поверхностью; В - коэффициент турбулентной диффузии в той точке.

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

В пределах каждого элемента концентрация взвеси я неизменна в пространстве и может меняться только лишь во времени [13].

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

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

Каждую клетку расчётной сетки обозначим индексами г и у, показывающими её положение соответственно по координатным направлениям х и г (рис. 3). На этом же рисунке показаны обозначения граней расчётных элементов а и Ь - тыльная и фронтальная грани, считая по оси х; с и С - соответственно тыльная и фронтальная по направлению оси г.

В случае стационарного режима течений уравнение неразрывности течений контрольного элемента у будет иметь вид:

Чгуа - ЧгуЬ + Чгуо - ЧгуС = 0 (2)

где: qija - расход воды через грань а для контрольного элемента у; дф - расход воды через грань Ь для контрольного элемента у; qiJ■с - расход воды через грань с для контрольного элемента у;

qijd - расход воды через грань d для контрольного элемента ij.

Рисунок 3 - Схема расположения расчетного элемента в двухмерной

модели

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

Ча = ^ ' На , (3)

где: ча - расход воды через грань а;

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

Соответственно этому условию должно удовлетворять с определенной степенью точности поле скоростей, используемое для расчета концентрации.

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

Эя = к ■ (4)

где: я - концентрация взвешенного вещества; ? - время;

к - коэффициент уравнения. Для дальнейшего упрощения математических выкладок введём функцию ф, полученную на основании уравнения (4):

$ =к■ 4, (5)

/

где: 1 - концентрация взвешенного вещества в пределах контрольного элемента 1 в интервал времени с порядковым номером к; At - интервал времени.

Результирующий вертикальный секундный перенос чяп твёрдых частиц, относящихся к фракции п (по значению гидравлической крупности ип) через единицу граничной поверхности дна и водных масс может быть выражен, согласно Методике (1987), следующей формулой:

Чяп = (ип + Еп ) ■ Яп - Еп ■ Эазм.п , (6)

где: яп - средняя по вертикали частная мутность (отвечающая концентрации частиц крупностью ип; $взм.п - частная мутность взмыва для той же фракции, определяемая по формуле (7):

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

а ■ 8 а ■ Ь ■ N ■ ¥т

_ взм.п взм „ _ взм.п

взм.п 1 взм.п 1 , ч

100 или 100 , (7)

где: 8взм - общая мутность взмыва;

авзм.п - процентное содержание п-й фракции в составе взвешенных

частиц; Ь - коэффициент уравнения;

N - безразмерное характеристическое число, зависимое от

коэффициента Шези (С) и её производной (М); Гг - число Фруда.

Параметр Еп вычисляется по формуле (8):

и ■ Г Е =-п-

п 1 - Г , (8)

где: Г - гидромеханический параметр (находится по справочной таблице), зависящий от коэффициента Шези (С) и безразмерной величины:

а = и,

у , (9)

где: V- средняя скорость.

Расчёты ведутся по средней гидравлической крупности частиц (ип). Формула (6) может быть разделена на две части, одна из которых определяет взмыв частиц (индекс «в»), а другая - их осаждение на дно (индекс «о»). При этом согласно Методике (1987), из формулы (6) непосредственно получаем следующие функции модели:

Р

тув тт

нф

к (ип + Еп ).. &

(10)

. (11) Общее расчётное уравнение модели расчета изменения

^+1

концентрации взвешенных частиц ( ф ) приобретает вид:

к+1 „ к . . . . ЬV = Ьк + Рг]а - Рф + Рус - Рг]Л + Рут + Рг]в ~ Р.о

(12)

/

где: . - концентрация взвешенного вещества в расчётном элементе . в

интервал времени порядковым номером к;

к к к к Рф, рФ, р]с, и Рф - дифференциальные функции концентраций

взвешенных частиц на соответствующих гранях а, Ь, с и d

расчётного элемента ф в интервал времени порядковым номером

к;

к

Рг]т - дифференциальная функция концентрации взвешенных частиц дополнительно поступающих в расчётный элемент ф в интервал

времени порядковым номером к.

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

Модельный расчёт концентрации взвеси в пределах шлейфа дополнительной мутности

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

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

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

Б = Б + (Б Б ) е[~вя (и )Х/вД]

= + (Бяг - ) ■е , (13)

где: Бп - концентрация частиц расчётной крупности, соответствующая транспортирующей способности потока на рассматриваемом участке водотока;

8Н - концентрация взвешенных загрязняющих частиц, относящихся к ¿-ой фракции, в начальном поперечнике загрязнённой струи при X = 0;

Вд - действующий фронт работ;

и - средняя гидравлическая крупность взвешенных загрязняющих частиц (расчётная крупность);

Е1 - коэффициент, определяемый в зависимости от гидравлической

крупности частиц и гидромеханического параметра (Г) для ¿-ой фракции;

X - длина участка, в конце которой вычисляется Бх1;

Qд - действующий расход воды в пределах фронта работ.

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

Определив по формуле (14) максимальный диаметр частиц грунта, которые будут уноситься потоком, получим размер фракций грунта, по которым в дальнейшем будем проводить дальнейшие расчёты:

С - ^ 1000

V 3,6 у

Н , (14)

где: V- средняя скорость водотока;

И - средняя глубина водотока.

Для выделенных фракций вычисляется средняя гидравлическая крупность и.

Далее расчёт переноса наносов выполняется по формуле (15), позволяющей получить распределение концентрации наносов вдоль оси и вместе с тем охарактеризовать процессы взмыва и осаждения:

С С _1_/"С о ч А-Вд ■ (и +Е)■ ЛХ / Qд ]

8 К = БТ + (БН - Бт )■ е , (15)

где: - общая мутность (суммарная концентрация всех взвешенных в воде фракций) в конце расчётного участка длиной ЛХ;

Бт - мутность, соответствующая транспортирующей способности потока на рассматриваемом участке водотока;

Бн - мутность в начальном створе расчетного участка (в створе для которого определены Qд и ВД);

Bд - действующий фронт работ;

u - средняя гидравлическая крупность взвешенных загрязняющих частиц (расчётная крупность);

Е - коэффициент, определяемый в зависимости от гидравлической крупности частиц и гидромеханического параметра (Г);

Х - длина участка, в конце которой вычисляется SК;

Qд - действующий расход воды в пределах фронта работ.

Учитывая тот факт, что для случаев ведения производственных работ в водотоках формируется шлейф дополнительной мутности, то уравнение (15) трансформируется в следующую зависимость:

= (Sн - Sт ) ■ ^^/ *] . (16)

Расчёт распределения концентрации мутности по оси Х ведётся по формуле (16) от участка к участку. В пределах каждого участка гидравлические и морфометрические характеристики принимаются однородными. Полученное по расчёту значение мутности в конце первого расчётного участка является начальным для второго и т. д. Расчёт выполняется до получения значений мутности, обусловленных поставленной задачей (достижения допустимых значений).

Особенности расчёта концентраций взвесей при работах по разработке подводных каналов

Модель «Поток» 1.0 позволяет учитывать вид производимых в водотоке работ, при этом рассматриваются два основных вида земляных работ:

1) разработка подводных каналов;

2) отсыпка дамб.

Для каждого из них предлагается свой способ расчёта. В основу расчётов положено уравнение баланса наносов транзитной струи.

Первый способ расчётов позволяет оценить последствия разработки подводных каналов, осуществляемой способом землечерпания. На участке

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

В расчётах по первому способу важное значение имеет определение действующей ширины Вд. Эта величина определяется как проекция фронта проводимых работ на направление перпендикулярное оси Х, и соответствует ширине потока в начальном створе. Под фронтом работ понимается участок производства работ протяженностью ВФ в плоскости зеркала водотока и глубиной ЬФ, измеряемой по перпендикуляру к ВФ. Для вычисления Вд предлагается использовать формулу (17):

ВД _ Вф sin b + Ьф eos b (17)

где: Вф - фактическая ширина фронта работ;

Ьф - фактическая глубина разработки грунта;

в - угол между линией фронта проводимых работ и направлением

течения (изменяется от 0 до 90°).

Значение мутности взмыва находится по формуле (18):

„ _ 0,15-N-v2 ""'_ h , (18)

где: h - средняя глубина водотока;

v - средняя скорость течения в водотоке;

N - безразмерное характеристическое число, зависимое от

коэффициента Шези (С) и её производной (М). Мутность в начальном створе расчётного участка рассчитывается по формуле (19):

_ P-P0-m-(0,01-^Рг) Л Н _ О

Уд , (19)

где: P - производительность одного ковша машины;

р0 - плотность смеси воды и грунта в ковше;

т - понижающий коэффициент к производительности землеройных машин при разработке грунта из под воды;

Qд- действующий расход воды в пределах фронта работ;

- суммарный процент содержания фракций транспортируемых во взвешенном состоянии.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

При выполнении расчётов учитывается, что в большинстве случаев

р0<р грунта.

Особенности расчёта концентраций взвесей при работах по отсыпке дамб

Для учёта особенностей распространения и седиментации дополнительной мутности при работах по отсыпке дамб, согласно Методике (1987), в формулы (17)-(19) вносятся следующие изменения. В формулу мутности взмыва (18) вводится дополнительный множитель (1 - вту), учитывающий угол откоса у, формирующегося при отсыпке грунта в голове дамбы:

„ _ 0,3 ■ 2

вЗМ. 7/1 • \

к ■(1 - вт Г) . (20)

Поскольку поступление наносов в районе откоса дамбы происходит только за счёт естественного размыва под влиянием течения водотока, начальная мутность в этом случае будет равна мутности взмыва на естественном откосе с углом у, то есть:

SН _ $взм (21)

Под фронтом работ в данном случае понимается ширина зоны, расположенной в торцевой части дамбы, где производится отсыпка грунта. Принимается, что направление течения перпендикулярно оси дамбы, т. е. направлено параллельно её торцевой части. При этом ширина области взмучивания, обозначается также Вд, и соответствует размерам конуса отсыпки грунта. Значение Вд определяется по формуле (22):

Я = h

вд ~-

tgg , (22)

где: у - угол откоса насыпи дамбы.

Модельный расчёт площади и ширины шлейфа дополнительной мутности

Для проведения дальнейших расчётов и определения интегральных параметров проводится расчёт необходимого набора площадей Pi и поперечных размеров облака взвеси B2i. Согласно экспериментам угол расхождения потока взвеси составляет 11-14°, обычно его принимают равным 13°.

Для расчёта площади занятой взвесью, и поперечных размеров облака взвеси, необходимо задать следующие параметры (рис. 4):

bcp - средняя ширина живого сечения русла;

- ширина облака взвеси в створе проведения работ;

- расстояние от левого берега до облака взвеси;

о: - угол расхождения потока.

В расчётах ширина B1 принимается равной фактической ширине фронта работ Вф. В случае малых водотоков величину В1 можно считать равной ширине реки.

Для проведения расчётов русло считается прямолинейным. В случае расстояния до створа 150 м и менее можно принять гипотезу о том, что интенсивность диффузионного перемешивания зависит линейно от интенсивности конвективного переноса и составляет приблизительно 23,1% (что соответствует тангенсу угла в 13°). При расстояниях до створа более 150 м преобладание конвективного переноса над диффузионными процессами усиливается.

Рисунок 4 - Схема расположения основных параметнов распространения

взвеси в водотоке

Площадь облака Р(Ь) рассчитывается по формуле (23), при этом применяемые величины соответствуют обозначениям рисунка 4:

= хЬ + Д((£т2 1дау/2,0 < Ь <141е>(^уЦ1)/2 +

-141 ),141 <

(23)

=

¿2 =

Ъср - В -Ду

Ширина облака ^(Ь) находится по формуле (24):

В 2( Ь) = вх +

\Ь • tga, 0 < Ь < Ь I Ду, Ь < Ь

Ь • (а, 0 < Ь < Ь1 ] -В1 -Ду, Ь2 <Ь|

^ ^ ^ . (24)

Границами расчётных зон распространения мутности принимаются створы, до которых происходит полное оседание на дно выделенных фракций грунта.

Значения толщины слоя наилка и удельной плотности заиления рассчитываются как средние в зоне, поэтому эти значения относятся к створу реки в середине зоны.

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

Модельный расчёт интегральных показателей Определение интегральных показателей (объём воды, протёкшей через области шлейфа с заданной концентрацией взвеси; средний слой наилка; удельная плотность заиления дна в расчётной зоне) проводится по Методу ГГИ [11], для этого рассчитываются границы зон полного осаждения фракций заданного размера (Ь):

НУ

Ьг =

и , (25)

где: Н - средняя глубина водотока;

у - средняя скорость течения в водотоке;

ы{ - средняя гидравлическая крупность взвешенных загрязняющих частиц (расчётная крупность). Далее каждой зоны определяется масса грунта каждой фракции

■ осевшая в ¿-ой зоне:

^ч = - У ]'") /о

Ь , (26)

где: - масса грунта, поступившая в поток ¿-ой фракции;

Ь - длина расчётной зоны;

Ь - расстояние от створа работ до створа полного оседания расчётной фракции;

у ^

¿^ ¿(]-п) - масса грунта данной фракции, осевшая в предыдущих расчётных зонах выше по течению.

Суммированием ■ получаем массу грунта, осевшую в каждой зоне.

В результате формируется интегральная сумма массы осевшего грунта в русле от створа проведения работ до каждого расчётного створа с учётом массы транзитного грунта.

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

Работы могут производиться последовательно или одновременно на нескольких участках водотока. Если работы производятся одновременно в разных створах речного русла, облака взвеси могут пересекаться; в этих областях плотности будут суммарными. Отложения на дне, образующиеся от производства работ на различных участках также будут суммироваться. Для упрощения расчётов будем считать, что протяженность участка, на котором проводятся работы, незначительна по сравнению с предполагаемой длиной облака мутности, и, можно полагать, что взвесь поступает в воду в одном створе. Для более сложных случаев общий результат можно получить, как суперпозицию результатов расчёта отдельных простых конфигураций. При этом необходимо сохранять промежуточные результаты для получения в итоге решения общей задачи.

Исходные данные для проведения моделирования вносятся в соответствующие графы программы (рис. 5).

Рисунок 5 - Пример формирования исходных данных для моделирования распространения облака взвеси, образующейся в ходе работ в русле

водотока

Получаемые в результате моделирования промежуточные результаты имеют следующий вид (рис. 6).

Максимальный диаметр частиц, уносимых потоком с1т, мм = 1.19Е-06

Средняя гидра отческая крупность транспортируемых частиц ис, м/сек. = 3,96875Е-05

Концентрация взвеси соответствующая транспортной способности водотока Э, г/л = 0,013268

Время воздействия допогнительной мутности на водоток!, сек. = 295,24 Расстояние до створа достаточного осаждения и разбавления взвеси Хт, м = 23.97066 Расстояние до створа осаждения и разбавлении взвеси до концентрации 100 мг/л, XI00, м = 15,52551 Расстояние до створа осаждения и разбавления взвеси до концентрации 20 мг/л, Х20, м = 18,30349 Толдина слоя приращения осажденных взвесей до концентрации ПДКИ пдк, мм = 32.25597 Толдина слоя осажденных взвеоей до створа с концентрацией 100 мг/л, И1 00, мм = 32,252 Толдина слоя осажденных взвесей до створа с концентрацией 20 мг/п. И20. мм = 32.2552

Рисунок 6 - Пример промежуточных результатов моделирования распространения облака взвеси, образующейся в ходе работ в русле

водотока

Изменение концентрации взвеси в воде в зависимости от расстояния от места производства работ отображается также в графическом виде (рис.

7).

Рисунок 7 - Пример графического изображения данных об изменении концентрации взвеси в воде в зависимости от расстояния от фронта работ

Полученные в ходе имитационного моделирования итоговые данные вносятся в таблицу (табл. 2).

Таблица 2 - Пример формирования итоговой таблицы рассчитанных параметров распространения и интегральных показателей распространения шлейфа мутности

Средние расстояния распространения взвеси от источника, м