ИСТИНА

В настоящее время в низовьях р.Терек происходит интенсивная перестройка формы продольного профиля русла за счет сокращения длины русла спрямляющей прорезью через Аграханский полуостров и соответственного увеличения уклонов водной поверхности. Кроме того, проектируемое сооружение выше Каргалинского гидроузла руслового- водохранилища обусловит аккумуляцию наносов в верхнем створе и приведет к вертикальным переформированиям русла (эрозии) в нижнем бьефе» Для расчета вертикальных переформирований применяется обычная система одномерных уравнений неравномерного движения и деформации. В уравнении движения основное внимание обращалось на определение коэффициента Шези. На основании специальных четырехлетних наблюдений за гидравликой потока на эталонных участках реки была получена двучленная формула для определения коэффициента Шези: Эта формула учитывает как Зернистую, так и грядовую шероховатость, при этом под грядами понимается весь комплекс русловых образований от рифелей до побочней. В уравнении деформации используются также специально разработанные для условий нижнего Терека формула для" расчета стока рус лообразующих наносов. Для распределения руслообразующих наносов в потоке по вертикали показана целесообразность использования формулы Рауза-Великанова с учетом зависимости постоянной" Кармана от гидравлической крупности. Для вычисления придонной мутности для каждой фракции наносов с учетом их содержания в донном грунте также получена региональная формула. Для ее построения было использовано более 200 измерений придонной мутности вакуумным батометром (на горизонте 5-7 см от дна в диапазоне глубин 1,0-8,0 м и скоростей 0,4-2,3м/с). Показано увеличение мутности с ростом числа Рейнольдса, а не числа Фруда, как в большинстве формул. Решение уравнения движения при заданных по сетке отметках дна, ширинах русла, гранулометрическом составе донных наносов, начальной и конечной глубине и расходе воды ведется методом Рунге-Кутта. Последовательность операций следующая: а) рассчитываются отметки свободной поверхности снизу вверх; в) рассчитанная конечная глубина сравнивается с заданной; в) вносятся поправки в коэффициенты. Далее цикл повторяется до достижения необходимой точности определения конечной глубины. Для ускорения расчетов необходимо выделить по длине реки несколько отрезков, в начале и конце которых можно задать глубины при данном расходе воды (по эмпирическим связям или по формуле Шези), В результате расчета получаются гидравлические характеристики, -необходимые для решения уравнения деформации. Уравнение деформации решается в конечных разностях за дискретный отрезок времени (декаду), в результате чего вносятся поправки в начальные отметки дна и в гранулометрический состав активного слоя грунта. Для нового расхода воды последовательность операций повторяется. Изложенный алгоритм реализован на ЭВМ БЭСМ-6. Рассчитаны изменения продольного профиля реки Терек и отметки уровня паводковых вод для разных вариантов регулирования стока водохранилищами и для периодов разной водности реки.