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

Для цитирования: Григорьян О.П., Магомедова М.Р. ИМИТАЦИОННОЕ МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ РАСПРОСТРАНЕНИЯ ВЗВЕШЕННЫХ ЧАСТИЦ ДЛЯ ОПРЕДЕЛЕНИЯ ЗОН ПОВЫШЕННОЙ МУТНОСТИ ПРИ ОЦЕНКЕ УЩЕРБА ВОДНЫМ БИОЛОГИЧЕСКИМ РЕСУРСАМ РЕКИ ЧЕРЕК. Вестник Дагестанского государственного технического университета. Технические науки. 2016; 42 (3):92-100. DOI:10.21822/2073-6185-2016-42-3-92-100.

For citation: Grigoryan O.P., Magomedova M.R. IMITATION MATHEMATICAL MODELING OF SUSPENDED PARTICLES EXPANSION TO DETERMINE THE AREAS OF INCREASED TURBIDITY IN THE ASSESSMENT OF DAMAGE TO WATER BIOLOGICAL RESOURCES OF THE CHEREK RIVER. Herald of Dagestan State Technical University. Technical Sciences. 2016;42 (3):92-100. (In Russ.) DOI:10.21822/2073-6185-2016-42-3-92-100

ТЕХНИЧЕСКИЕ НАУКИ

ИНФОРМАТИКА, ВЫЧИСЛИТЕЛЬНАЯ ТЕХНИКА И УПРАВЛЕНИЕ

УДК 556.536:631.4

DOI: 10.21822/2073-6185-2016-42-3-92-100

1 2 Григорьян О.П. , Магомедова М.Р.

Дагестанский филиал Каспийского научно-исследовательского института рыбного хозяйства, 367022,г. Махачкала, ул. Абубакарова, д.104 e-mail:dokaspiy@mail.ru, 2Дагестанский государственный технический университет, 367015 г. Махачкала, пр. И. Шамиля, 70, n-

guseinova@mail. ru ИМИТАЦИОННОЕ МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ РАСПРОСТРАНЕНИЯ ВЗВЕШЕННЫХ ЧАСТИЦ ДЛЯ ОПРЕДЕЛЕНИЯ ЗОН ПОВЫШЕННОЙ МУТНОСТИ ПРИ ОЦЕНКЕ УЩЕРБА ВОДНЫМ БИОЛОГИЧЕСКИМ РЕСУРСАМ РЕКИ ЧЕРЕК Аннотация. Цель. Оценить ущерб водным биоресурсам от эксплуатации Нижне-Черекских ГЭС(Кашхатау, Аушигерской). Методы. Приведены результаты апробации имитационной математической модели распространения взвешенных частиц для определения зон повышенной мутности при оценке ущерба водным биологическим ресурсам реки Черек от эксплуатации каскада Нижне-Черекских ГЭС. Для установления основных характеристик годового и сезонного стока в створах гидротехнических сооружений на реке Черек использованы сведения о стоке в опорных гидрологических створах Гидрометеослужбы, обобщённые в справочниках «Ресурсы поверхностных вод». Результат. Математическая модель, разработанная на основе теории вероятностей и теории выбросов случайных процессов, с учетом нормального закона распределения горизонтальной и вертикальной составляющих мгновенных скоростей течения и закона Релея для распределения их максимумовиспользована для расчета «шлейфа мутности» при промыве отстойника головного узла и при глубинном промыве головного узла Кашхатау ГЭС. Вывод. В результате проведенных исследований получено, что негативное влияние «шлейфа мутности» в реке Черек, будет прослеживаться на протяжении приблизительно 3 км.

Ключевые слова: транспорт наносов, взвешенные наносы, гранулометрический состав,

скорость течения, расход наносов, турбулентный поток, средняя глубина

* 1 * 2 Olga P. Grigoryan , Milada R. Magomedova

1Federal State Budgetary Scientific Institution,

Caspian Fisheries Research Institute, Daghestan branch,

104, Abubakarov Ave, Makhachkala,367022,

2Daghestan State Technical University,

701. Shamil Ave, Makhachkala, 367015

IMITATION MATHEMATICAL MODELING OF SUSPENDED PARTICLES EXPANSION TO DETERMINE THE AREAS OF INCREASED TURBIDITY IN THE ASSESSMENT OF DAMAGE TO WATER BIOLOGICAL RESOURCES OF THE CHEREK RIVER

Abstract. Aim. Assessment of damage to aquatic bioresources. Methods. This article presents the results of testing of the simulation mathematical model of propagation of suspended particles to determine the areas of increased turbidity in the assessment of damage to Cherek River water biological resources because of Nizhne-Cherek HPS Cascade exploitation. In order to establish the basic characteristics of annual and seasonal runoff in the hydraulic structures dam locations on the river Cherek are used data on runoff in support hydrological dam locations of the Hydrometeorological Service, generalized in the reference book "Surface water resources." Results. A mathematical model developed on the basis of probability theory and stochastic processes emissions theory, based on the normal law of distribution of horizontal and vertical components of the instantaneous flow velocity and Rayleigh law for the distribution of their peaks is used to calculate the " turbidity plume " in the washing of the head unit settler and deep washing of Kashkhatau HPS head unit. Conclusion. The studies found that the negative impact of "turbidity plume " in the Cherek river will be traced for approximately 3 km.

Key words: drift transport, suspended drift, granulometric composition, current speed, drift discharge, turbulent flow, average depth

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

Река Черек - правобережный приток р. Баксан образуется от слияния рек Черек Безен-гийский и Балкарский. Истоки притоков находятся в ледниках центральной части Большого Кавказа. Соединившись в один поток, Черек Балкарский и Безенгийский текут одним руслом, лишь в отдельных местах разбиваясь на ряд рукавов. Длина реки составляет 79 км, водосборная площадь равна 3070 км2. Бассейн расположен в пределах высот от 4000 м в верховьях до 350 м - в устьевой части. Площадь водосбора в створе плотины Советской (Кашхатау) ГЭС составляет 1350 км2.

Материалом для проведения расчета ущерба водным биологическим ресурсам от эксплуатации Кашхатау (Советской) и Аушигерской ГЭС послужили комплексные гидробиологические исследования за ряд лет, проводившиеся ДФ ФГБНУ «КаспНИРХ», в том числе результаты гидробиологических съемок, осуществленных в районе влияния гидроэлектростанций в летне-осенний период 2015 г.

При возведении каскада Нижне-Черекских ГЭС в основном русле реки Черек построены плотины Кашхатау ГЭС (головная) и Аушигерской ГЭС (резервная).

Примерное расстояние от плотины Кашхатау ГЭС до плотины Аушигерской ГЭС (по руслу реки) составляет 8 км. Расстояние от плотины Аушигерской ГЭС до точки сброса отработанной Аушигерской ГЭС воды составляет также 8 км.

Из-за забора большей части расхода реки Черек в деривацию на участке реки от водозаборной плотины Кашхатау ГЭС до точки сброса отработанной Аушигерской ГЭС воды происходит снижение площади русла реки на 68%.

Протяжённость находящегося в зоне воздействия каскада Нижне-Черекских ГЭС участка реки Черек составляет 16 км. При средней ширине русла реки 30 м, площадь русла реки Черек, подвергающегося воздействию каскада ГЭС составляет порядка 480 000 м2 (или 48 га).

Потеря площади русла реки от забора большей части расхода воды в деривацию каскада Нижне-Черекских ГЭС (68%) составит 32,64 га или 326 400 м .

Для определения параметров и расчётных величин годового стока р. Черек был использован 62-х летний ряд наблюдений за стоком в опорном гидрологическом створе -Кашкатау. Параметры кривой обеспеченности и средние годовые расходы воды различной расчетной обеспеченности реки Черек в створе плотины Кашхатау ГЭС приведены в таблице 1.

Таблица 1-Параметры кривой обеспеченности и средние годовые расходы воды

Створ л 2 А, км Мо, л/скм 2 Параметры кривой обеспеченности Расходы воды, м/с различной обеспеченности, Р, %

Q, м /с Cv Cs

10 50 75 95

Советская ГЭС 1350 30,2 40,8 0,13 2Cv 47,7 40,6 37,1 32,6

Внутригодовое распределение стока р.Черек в створе плотины головного узла Советской (Кашхатау) ГЭС за характерные по водности годы приведено в таблице 2.

IV V VI VII VIII IX X XI XII I II III Год

Средние месячные расходы воды, м3/с

Средний по водности год

15,0 29,9 72,5 114 114 58,6 23,9 15,9 12,5 10,8 10,3 9,85 40,8

Маловодный год Р~75%

18,8 34,2 76,6 96,5 92,9 48,6 21,8 14,0 10,9 10,8 9,18 8,29 37,1

Маловодный год Р~95%

16,7 30,3 65,9 84,2 85,3 43,7 19,2 11,9 9,22 8,83 7,68 6,53 32,6

Параметры кривой распределения вероятностей превышения максимальных расходов и расчётные расходы воды приведены в таблице 3. _Таблица 3-Параметры кривой распределения вероятности превышения_

Створ Параметры кривой распределения вероятности превышения Расходы воды ( м/с) различной вероятности превышения, Р%

Qср.макс., м3/с а Cs 0,1 0,5 1 3 5 10

Кашхатау (Советская) ГЭС 202 0,26 4Cv 459 392 364 319 309 271

При сбросе в русло реки Черек отработанной гидроэлектростанциями каскада Нижне-Черекских ГЭС воды не образуется «шлейф мутности», так как в гидроагрегаты поступает отстоянная вода, которая сбрасывается в реку с малой скоростью (кинетическую энергию воды забирают гидроагрегаты).

Постановка задачи. «Шлейф мутности» образуется при промыве отстойника головного узла Кашхатау ГЭС и при глубинном промыве головного узла Кашхатау ГЭС. Мутность смываемой с отстойника воды составляет 41 кг/м3.

Основная часть промывов отстойника проводится в летние 3 месяца. В этот период в средний по водности год средний сток реки Черек составляет 104 м3/с. Из них в деривацию за-

3 3

бирается 65 м /с, на промыв отстойника расходуется 15 м/с воды, расход оставшейся в реке воды составляет 24 м3/с.

Использованная для промыва отстойника вода с повышенной мутностью перемешивается с оставшейся в реке водой и образуется общий сток 39 м3/с. При этом мутность воды снижается и составит в месте смешивания 15,8 кг/м .

С учётом мутности самой воды в реке (примерно 0,2 кг/м ) мутность воды в нижнем бьефе при промывке отстойника составляет 16 кг/м .

Среднегодовая продолжительность промыва отстойника составляет 210 часов (35 промывов по 6 часов каждая). Треть (33%) создающих мутность воды наносов имеют линейные размеры не более 0,2 мм, а фракции с размерами не более 2 мм составляют 95% от общей массы смываемых наносов.

Один раз в год в паводковый период проводится глубокий промыв наносов головного узла Кашхатау ГЭС средней продолжительностью 5 суток. Половина (50%) создающих мут-

ность воды наносов имеют линейные размеры не более 0,2 мм, а фракции с размерами не более 7 мм составляют 90% от общей массы смываемых наносов.

Во время проведения глубинного промыва наносов головного узла Кашхатау ГЭС гидростанции каскада не работают. Одновременно с глубинным промывом наносов головного узла Кашхатау ГЭС открываются затворы на плотине Аушигерской ГЭС и мутная вода «транзитом» проходит через резервный головной узел.

В связи с тем, что резервный головной узел заполнен наносами полностью, с него практически наносы не удаляются при глубоком промыве глубинным промывом наносов головного узла Кашхатау ГЭС.

Структура взвешенных наносов реки Черек представлен следующим гранулометрическим составом (табл.4):

d, мм >2 2-1 1-0,5 0,5-0,2 0,2-0,1 0,1-0,05 0,05-0,01 <0.01

% 4,4 6,0 3,8 18,9 17,1 8,3 21,6 19,9

Исходя из стохастической природы турбулентности водных потоков и случайного характера процессов распространения взвешенных наносов, при разработке математической модели транспорта взвешенных наносов турбулентным потоком при промыве отстойника головного узла и при глубинном промыве головного узла Кашхатау ГЭС, был принят стохастический подход с использованием теории вероятностей и теории выбросов случайных процессов с учетом нормального закона распределения горизонтальной и вертикальной составляющих мгновенных скоростей течения и закона Релея для распределения их максимумов [5, 8, 9, 12].

Методы исследования. Математическая модель распространения взвешенных наносов в реке Черек и алгоритм расчета расхода наносов в естественном деформируемом речном русле основаны на представлении о едином механизме движения придонных и взвешенных наносов как стохастическом процессе взаимодействия турбулентных возмущений водного потока с частицами грунта [1, 2, 3, 4, 11,15].

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

(у — а)в =

f

1 - Ф

— а

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

ст.,

У в =

1 — Ф

У в =

Ф

*3' ,

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

уВ =42^уВ f (хз), 1 — Ф(хз ) '

Тв = --

yl2^yBf (хз ).

Средняя частота пульсаций лагранжевой вертикальной скорости жидкости равна:

-л V Ув =-

23И ;

Среднее расстояние вдоль потока между смежными взвешивающими возмущениями в придонной области:

L» = —f

V .

Средняя концентрация указанных возмущений над поверхностью русла:

1

Т

V

'В — а ;

Гидравлическую крупность частицы определяем по формуле В.Н. Гончарова[3]:

1

о =

и

2g(рг - p)d

1.75р

где d - диаметр частиц грунта, м, Р - плотность воды, Р = 1000 кг/м3, Рг - плотность матеР 3 2 0 риала частиц грунта, г = 1800 кг/м , g - ускорение свободного падения, g = 9,81 м/с , ^ -

параметр турбулентности, зависит от крупности частиц грунта.

Дальность переноса донных и взвешенных наносов определяем по зависимостям:

/ = Vt -{Р'\ Р/2)0 ln

gpr - Р) (Рг + Р/ 2)о

^ дон = —Ас tдон ( \

g

ln

1 + g(рг -р) Vje (рГ + р / 2 )о2 ю в

. + g(Рг -Р) — t

1 + -Г—-—Ас t

(рг + р/ 2)о

дон

Для практических расчетов взвешенных частиц формулу (10) можно использовать в виде:

^ аз = ^аз^аз

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

N

где

gt = 2 qt к bk

k=i

Чт , К = Чдонк + Чв:

Чт , К , Чдон, к h

юн,к'ччз,к- соответственно удельные расходы руслоформирующих, донных и

К

взвешенных наносов, к - ширина к-й струи.

Длину участка, на котором выпадут наносы, т.е. путь на протяжении которого выпадут наносы рассматриваемой крупности, определяем по формуле:

L = HK

v

ср

где н - средняя глубина на участке длиной L, Ув - гидравлическая крупность наносов,

V

ср - средняя скорость течения воды.

Определяя предельную концентрацию взвешенных наносов над единицей площади дна

в наносонесущем слое потока, как Pnb nb lb , а скорость переноса их, как Vß3 ^вз / tß:

с учетом

параметров Раз и ^аз, формулу для предельного расхода взвешенных наносов используем в виде:

qS3 = qdanVereVe

P d i

вз вз вз

P ddi дон

Обсуждение результатов. Разработка и выполнение программ на ЭВМ осуществляется в студии разработчика Microsoft Developer Studio (MDS) на алгоритмическом языке FortranPowerStation (стандарт языка Фортран-90) [7, 8, 10].

Разработанная математическая модель распространения взвешенных наносов в русле реки Черек реализована в программном комплексе Sediment Transport Cherek.

Массив охватывает следующий диапазон изменения гидравлических характеристик потока и русла (табл.5-8): [13-14]

1. Характеристики «шлейфа мутности» при промыве отстойника головного узла Кашхатау ГЭС: расход воды Q = 39 м3/с, объем наносов V = 7 400 м3 , глубина потока h = 0,7 м, ширина потока b = 30 м, средняя скорость потока V = 1,9 м/с, диаметр частиц наносов d = 0,002 м, мутность воды ц=16 кг/м.3

В результате объемы загрязняемой воды, протекающей в шлейфах взвеси с концентрацией частиц наносов выше заданной, — «проточные объёмы» (ПО), м3получаем в виде:

Таблица 5 - Объемы загрязняемой воды

Концентрация взвеси в воде, кг/м3

> 0,100 > 0,050 > 0,020 > 0,010 > 0,005 >0,001

518400 259200 103680 51840 25920 5184

Расстояние от створа гидроузла до места выпадения осадков с концентрацией частиц наносов выше заданной, представлено в таблице 6:

Таблица 6 - Расстояние от створа гидроузла до места выпадения осадков

Концентрация взвеси в воде, кг/м3

Операция > 0,100 > 0,050 > 0,020 > 0,010 > 0,005 >0,001

L, м 52,1 104,2 260 521 1041,7 5208,3

Таблица 7 - Время существования шлейфов взвеси с концентрацией _выше заданной, час.

Концентрация взвеси в воде, кг/м3

> 0,100 > 0,050 > 0,020 > 0,010 > 0,005 >0,001

0,076 0,15 0,38 0,76 1,15 1,76

Таблица 8 - Площади дна (м2), покрываемые слоем отложений взвешенных

веществ

Толщина слоя осадков, м

> 0,100 > 0,050 > 0,020 > 0,010 > 0,005 >0,001

1563 3126 7800 15630 31251 156249

2. Характеристики «шлейфа мутности» при глубинном промыве представлены в таблицах 9-12: расход воды Q = 100 м3/с, объем наносов V = 7 400 м3 , глубина потока h = 0,9 м, ширина потока Ь = 50 м, средняя скорость потока V = 2,2 м/с, диаметр частиц наносов d = 0,007 м, мутность воды ц= 10,6 кг/м.

В результате объемы загрязняемой воды, протекающей в шлейфах взвеси с концентрацией частиц наносов выше заданной, — «проточные объёмы» (ПО), м3получаем в виде показателей, представленных в таблицах 9-12.

Таблица 9 - Объемы загрязняемой воды, протекающей в шлейфах взвеси с _концентрацией частиц наносов выше заданной_

Концентрация взвеси в воде, кг/м3

> 0,100 > 0,050 > 0,020 > 0,010 > 0,005 >0,001

381600 190800 76320 38160 19080 3816

Таблица 10 - Расстояние от створа гидроузла до места выпадения осадков _с концентрацией частиц наносов выше заданной _

Концентрация взвеси в воде, кг/м3

Операция > 0,100 > 0,050 > 0,020 > 0,010 > 0,005 >0,001

L, м 28 56 139 279 557 2786

Таблица 11 - Время существования шлейфов взвеси с концентрацией

выше заданной

Концентрация взвеси в воде, кг/м3

> 0,100 > 0,050 > 0,020 > 0,010 > 0,005 >0,001

0,03 0,07 0,2 0,35 0,7 1,35

2

Таблица 12 - Площади дна (м ), покрываемые слоем отложений

взвешенных веществ

Толщина слоя осадков, м

> 0,100 > 0,050 > 0,020 > 0,010 > 0,005 >0,001

1400 2800 6950 13950 27850 139300

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

Произведена оценка ущерба водным биоресурсам от эксплуатации Нижне-Черекских ГЭС (Кашхатау, Аушигерской). Результаты проведенных исследований показали, что негативное влияние «шлейфа мутности» при однократном смыве селевых наносов [16-17] на одной электростанции, с учетом критических показателей концентрации взвешенных веществ воде, будет прослеживаться на протяжении почти 3 км в проточном объеме около 710 тыс. м3.

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

Библиографический список:

1. Алексеевский Н.И. Формирование и движение речных наносов. М.: Изд-во МГУ. 1998. - С.95-114.

2. Барышников Н.Б., Пагин А.О. Гидравлическое сопротивление речных русел// Журнал университета водных коммуникаций. - СПб, т.2 , 2010. - С.90-93.

3. Гришин Н.Н. Механика придонных наносов. М.: Наука, 1982. - С.102-119.

4. Копалиани З.Д., Костюченко А.А. Расчеты расхода донных наносов в реках: сборник работ по гидрологии. - СПб.: Гидрометеоиздат, 2004. - № 27. - С. 25-40.

5. Клавен А.Б., Копалиани З.Д. Экспериментальные исследования и гидравлическое моделирование речных потоков и руслового процесса. - СПб.: Нестор-История, 2011. - С. 103-107.

6. Магомедова М.Р. Математическое моделирование движения придонных наносов в открытых руслах: Махачкала: Алеф, Овчинников М.А., 2014. - С.53-68.

7. Магомедова А.В., Гусейнова М.Р., Насибова Н.В. Разработка программного комплекса для расчета транспорта руслоформирующих наносов в речных руслах с использованием ГИС-технологий //Вестник Дагестанского государственного технического университета. Технические науки - 2009. - Т.15 - №4- С.147-157.

8. Магомедова А.В., Магомедова М.Р. О факторах, обусловливающих процесс транспорта руслоформирующих наносов //Вестник Дагестанского государственного технического университета. Технические науки - 2013. - Т.29. - №2- С.58-64.

9. Магомедова М.Р. Практическое применение авторской модели транспорта минеральных частиц//Вестник Дагестанского государственного технического университета. Технические науки - 2015. -Т.37, - №2 - С.84-91.

10. Магомедова А.В. Эрозионные процессы в руслах рек и каналов. - М.: ВЗПИ, 1990. -С.98-108.

11. Магомедова М.Р. Разработка программного комплекса для обработки многомерных данных движения руслоформирующих наносов // Научно-методический электронный журнал «Концепт». - 2016. - Т. 15. - С. 1951-1955. - URL: http://e-koncept.ru/2016/96309.htm.

12. Магомедова М.Р. Программный комплекс для расчета расхода придонных наносов в речных руслах. А.с. 2013610555 Рос. Федерация, № 2013616724; заявл. 10.01.13; опубл. 17.06.13.

13. Ржаницын Н.А. Руслоформирующие процессы рек. - Л.: Гидрометеоиздат, 1985. -С.127-130.

14. Brian W.D., Peter F.F. Grain Size, Sediment Transport Regime, and Channel Slope in Alluvial Rivers. The Journal of Geology, 2011, vol. 106, no. 6, pp.662-673.

15. Benoit C., Magnus L.A General formula for non-cohesive bed load sediment transport. Estua-rine, Coastal and Shelf Science. 2005, pp.251-258.

16. Kopaliani Z.D. Problem of bed load discharge assessment in rivers. Proc. 10-th Inter Symp. on River Sedimenteion. Moscow, 2007, vol. 3, pp. 175-181.

17. Wilcock P.R., Crowe J.C. Surface-based transport model for mixed-size sediment. Journal of Hydraulic Engineering. 2003, vol.129, no.2, pp.120-128.

References:

1. Alekseevskiy N.I. The formation and movement of river sediments. Moscow: Moscow State University. 1998, pp.95-114.

2. Baryshnikov N.B., Pagin A.O. Hydraulic resistance of river channels. Zhurnal universiteta vodnyh kommunikacij. St. Petersburg [Journal of University of water communications] 2010, vol. 2, pp. 90-93.

3. Grishin N.N. The mechanics of the bottom sediment. Moscow: Nauka, 1982, pp.102-119.

4. Kopaliani Z.D., Kostyuchenko A.A. Calculations of the flow bed load in rivers: a collection of papers on hydrology. St. Petersburg: Gidrometeoizdat, 2004, no. 27, pp.25-40.

5. Klaven A.B., Kopaliani Z.D. Experimental studies and hydraulic modeling of river flows and channel process. St. Petersburg: Nestor-Istoriya, 2011, pp.103-107.

6. Magomedova M.R. Matematicheskoe modelirovanie dvizhenija pridonnyh nanosov v otkrytyh ruslah. [Mathematical modeling of the bottom sediment motion in open channels]. Makhachkala: Aleph, M.A Ovchinnikov, 2014, pp.53-68.

7. Magomedova A.V., Huseynova M.R, Nasibova N.V. Development of software for transport calculation river bed forming sediments in river beds using GIS technologies. Vestnik Dage-stanskogo gosudarstvennogo tekhnicheskogo universiteta. Tekhnicheskie nauki. [Herald of

Daghestan State Technical University. Technical science.] 2009, vol.15, no.4, pp.147-157.

8. Magomedova A.V., Magomedova M.R. Factors that determine the process of river bed forming sediment transport. Vestnik Dagestanskogo gosudarstvennogo tekhnicheskogo universi-teta. Tekhnicheskie nauki. [Herald of Daghestan State Technical University. Technical science], 2013, vol.29, no.2, pp.58-64.

9. Magomedova M.R. The practical application of the author's model of mineral particles transport. Vestnik Dagestanskogo gosudarstvennogo tekhnicheskogo universiteta. Tekhniches-kie nauki. [Herald of Daghestan State Technical University. Technical science], 2015, vol.37, no.2, pp.84-91.

10. Magomedova A.V. Jerozionnye processy v ruslah rek i kanalov. [Erosion processes in rivers and canals]. Moscow: VZPI, 1990, pp.98-108.

11. Magomedova M.R. Development of software for processing of multidimensional data of river bed forming sediments movement. Scientific and methodical electronic journal"Concept", 2016, vol.15, pp.1951-1955. URL:http://ekoncept.ru/2016/96309.htm.

12. Magomedova M.R. Software for calculating the flow rate of bottom sediment in the river beds. AS 2013610555 R.Federation, number 2013616724; appl. 01/10/13; publ. 06.17.13.

13. Rzhanitsyn N.A. Rusloformirujushhie processy rek. [River bed forming processes]. Leningrad: Gidrometeoizdat, 1985, pp.127-130

14. Brian W.D., Peter F.F. Grain Size, Sediment Transport Regime, and Channel Slope in Alluvial Rivers. The Journal of Geology, 2011, vol. 106, no. 6, pp.662-673.

15. Benoit C., Magnus L. A general formula for non-cohesive bed load sediment transport. Estua-rine, Coastal and Shelf Science. 2005, pp.251-258.

16. Kopaliani Z.D. Problem of bed load discharge assessment in rivers. Proc. 10-th Inter Symp. on River Sedimenteion. Moscow, 2007, vol. 3, pp. 175-181.

17. Wilcock P.R., Crowe J.C. Surface-based transport model for mixed-size sediment. Journal of Hydraulic Engineering. 2003, vol.129, no.2, pp.120-128.

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

Григорьян Ольга Петровна - заместитель директора по научной работе.

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

Authors information.

Olga P. Grigoryan - the Deputy Director on scientific work.

Milada R. Magomedova-candidate of technical Sciences, Dean of the faculty of oil, gas and environmental engineering, doctoral student in building structures and hydraulic structures of the Faculty of Architecture and Construction.

Конфликт интересов Conflict of interest

Авторы заявляют об отсутствии конфликта интересов. The authors declare no conflict of interest.

Поступила в редакцию 01.06.2016. Received 01.06.2016.

Принята в печать 10.08.2016. Accepted for publication 10.08.2016.