CC BY
43
Вестник Дагестанского государственного технического университета. Технические науки. №2 (29), 2013
УДК 556.536:631.4
Магомедова А.В., Магомедова М.Р.
О ФАКТОРАХ, ОБУСЛОВЛИВАЮЩИХ ПРОЦЕСС ТРАНСПОРТА РУСЛОФОРМИРУЮЩИХ НАНОСОВ
Magomedova A.V., Magomedova M.R.
ABOUT FACTOR, CONDITIONNING PROCESS OF THE TRANSPORT CHANNEL-FORMING SEDIMENTS
В данной статье на основе анализа литературных источников и обобщения экспериментальных и натурных данных уточнены основные факторы, влияющие на процесс транспорта руслоформирующих наносов. Процесс транспорта руслоформирующих наносов является одним из сложных физических процессов, обусловленных множеством факторов: гидромеханических характеристик водного потока и физико-механических характеристик наносов, что является причиной отсутствия единой теории этого процесса, сложности многочисленных расчетных методов определения расхода наносов и большого расхождения в результатах расчетов по ним.
Ключевые слова: транспорт наносов,донные наносы, взвешенные наносы, диаметр частиц наносов, гранулометрический состав, турбулентный поток.
In given article on base of the analysis of the literary sources and generalization experimental and натурных given are elaborated main factors, influencing upon process of the transport channel-forming sediments. The Process of the transport channel-forming sediments is one of the complex physical processes, conditioned by ensemble factor: an hydrоmechanical of the features of the water flow and physico-mechanical features alluvium that is a reason of the absence to united theory of this process, difficulties of the multiple accounting methods of the determination of the consuptionsediment and big divergence in result calculation on him.
Key words: sediment transport, bed load, suspended load, diameter of the particles sediment, granulometric composition, turbulent flow.
Процесс транспорта руслоформирующих наносов в турбулентном потоке является одним из сложных физических процессов, обусловленных множеством факторов: гидромеханических характеристик водного потока и физико-механических характеристик наносов; геометрический размер частиц наносов, форма частиц, плотность частиц грунта, плотность русловых отложений, гидравлическая крупность, средняя скорость потока, допускаемая неразмывающая донная скорость потока, фактическая донная скорость потока, относительная шероховатость русла, коэффициент подвижности наносов, вес частиц, гранулометрический состав грунта, глубина потока, уклон, расход воды, грядовая форма дна, параметр турбулентности зерен наносов и т. д.
Теоретические и экспериментальные исследования движения наносов в турбулентном потоке (Г.А.Эйнштейн, М.А. Великанов, А. Шильдс, Е. Мейер-Петер, Г.К. Джильберт, И.В. Егиазаров, И.И. Леви, Ц.Е. Мирцхулава, А.В. Караушев, Н.А. Михайлова, В.Н. Гончаров, К.В. Гришанин, К.И. Россинский, А.В. Магомедова, Н.Н. Гришин, В.С. Вербицкий, Н.И. Алексеевский и др.) с самого начала велись в двух направлениях: изучение донного влечения наносов и изучение взвешивания частиц и их переноса. По мнению ряда исследователей, такое разделение не совсем правильно, так как вследствие переноса твердых частиц турбулентными (вихревыми) возмущениями в потоке постоянно происходит обмен частицами руслоформирующих наносов между дном и
58
Вестник Дагестанского государственного технического университета. Технические науки. №2 (29), 2013 основной толщей потока, за исключением частиц неруслоформирующих фракций; при этом основной формой перемещения твердых частиц считают скачкообразное движение наносов (сальтацию), а донное влечение и взвешивание - предельными случаями сальтации при бесконечно малой и бесконечно большой длинах скачков (М.А. Великанов). При неоднородном по крупности гранулометрическом составе наносов в потоке одновременно имеют место все три вида движения [1, 5, 8].
Характер движения твердых частиц после срыва со дна - влечение, качение, сальтация вблизи дна или движение в толще потока в качестве взвешенных наносов, -зависит от запаса кинетической энергии, сообщенной частицам при срыве, и от величины сил, действующих на них во время движения. Различие между движением влекомых (перекатывающихся с остановками по дну), полувзвешенных (сальтирующих) и взвешенных наносов (поднимающихся в толщу потока и периодически вновь опускающихся ко дну) заключается в высоте и длине совершаемых твердыми частицами скачков [2].
Гидродинамическое воздействие потока на неподвижные частицы грунта на дне, вызванные несимметричным обтеканием частиц жидкостью, значительным градиентом скоростей в придонной области, наличием вихреобразования и другими факторами, приводит к появлению гидродинамических сил и сил сопротивляемости отрыву. Если частица лежит вровень с другими, то, как полагают ряд исследователей [1, 2, 5, 10], основную роль в силовом воздействии потока на эту частицу играет импульс подъемной силы, который вызывает срыв частицы и подъем ее на небольшую высоту над дном. Под действием силы тяжести и сопротивления воды, после прекращения действия подъемной силы, как только частица поднимется над дном, ускоренное движение сменяется замедленным, разность давлений исчезает и частица снова опускается на дно, т.е. падение частиц совершается под действием силы тяжести, а их подъем - против силы тяжести (Великанов, 1955). Поэтому влекомые наносы перемещаются короткими скачками. Вращательное движение частиц может сопровождаться действием дополнительной подъемной силы, называемой силой Магнуса [4].
Показатели, характеризующие процесс транспорта наносов, объединяют, в основном, в две группы факторов: характеристики турбулентности потока и физико-механические свойства грунтов. Исследования существующих факторов велись на основе теоретических, экспериментальных и чисто эмпирических методов.
Основными характеристиками физико-механических свойств русловых грунтов и наносов являются: средний диаметр частиц, гидравлическая крупность, плотность материала частиц и русловых отложений, гранулометрический состав грунта [8].
Важнейшей характеристикой транспорта наносов является гидравлическая крупность - скорость равномерного падения зерен наносов в спокойной воде, зависящая от диаметра частиц, удельной плотности материала и формы зерен [7]. Сложный характер взаимодействия твердых частиц с водным потоком привел к многообразию формул для расчета гидравлической крупности наносов, начиная с 1850 г. Исследованием и созданием шкал гидравлической крупности занимались многие отечественные и зарубежные ученые: Б.В. Архангельский, Г. Н. Лапшин, В.Н. Гончаров и др. В.Н. Гончаров, детально рассмотрев вопросы о режимах осаждения, построил шкалу гидравлической крупности наносов в воде в зависимости от диаметра частиц наносов при температурах воды 5, 10, 15 и 200С. На основе этой шкалы для расчета на ЭВМ гидравлической крупности на языке FortranPowerStation создан модуль Fall_velocity.
Среди существующих формул наиболее приемлемым для практического использования является способ выражения гидравлической крупности через стандартную гидравлическую крупность, предложенный В.Н. Гончаровым:
1
W = Р\
2g(Рг ~P)d
, (1)
1,75р
где р - параметр турбулентности, показывающий во сколько раз меньше действительная скорость падения в сравнении с той, которой обладали бы наносы при турбулентном режиме их падения. Расчет параметра турбулентности р в зависимости от температуры (5, 10, 15 и 20 градусов по Цельсию) принят по шкале, предложенной В.Н. Гончаровым. Для численной реализации на ЭВМ, на языке FortranPowerStation создан модуль Fall parameter. При d >1,5 мм, р = 1; в диапазоне 0,15 - 1,5 мм параметр турбулентности определяется по формуле [3]:
1
р = — 3
р Р V
р 1 2g (2)
рх -р V1,75рс1 '
для частиц С <0,15 мм р определяется по формуле:
33,8
Р= I _ . (3)
д/1,75р(р1 - р)С3
С уменьшением крупности зерен значение р возрастает, доходя при С <0,01 мм до р = 300.
Плотность частиц грунта, определяющая содержание массы частиц в единице объема, играет важную роль в процессе транспорта руслоформирующих наносов. Плотность частиц грунта зависит от их минералогического состава, для глинистых частиц и песков изменяется соответственно в пределах 2000 -4500 кг/м3. Для речных песков в среднем принимают 2650 кг/м3. Плотность отложений наносов зависит от их фракционного состава: для заиленных песков - 1100.. .1500 кг/м3, для песчаных грунтов -1500.. .1800 кг/м3, для песчано-гравелистых грунтов - 1800.2100 кг/м3, для гравийно -галечниковых грунтов - 2000.2400 кг/м3 [8].
Речные наносы представляют собой частицы различных размеров и форм. Форма частиц зависит от диаметра наносов, от степени механической обработки в процессе движения и играет во взаимодействии частиц с потоком значительную роль. Существуют различные морфологические критерии, характеризующие форму частиц: удлиненность, уплощенность, окатанность, сферичность, шарообразность и др. [8].
Диаметр частиц наносов является важнейшей характеристикой наносов. Определяющим геометрическим размером частиц является средний диаметр, соответствующий размеру круглых или квадратных сит, а не самих частиц [7]. В.Н. Гончаровым предложена следующая классификация наносов по размерам частиц (табл. 1)
[3].
В реках всегда транспортируются наносы различного размера, поэтому определение лишь одного среднего диаметра частиц представляется недостаточным.
Одним из важных факторов, влияющих на процесс транспорта наносов, является их гранулометрический состав. Достоверную информацию о гранулометрическом составе грунта дают плотность и функция распределения размеров частиц [8]. Гранулометрический состав русловых грунтов и наносов характеризуется средневзвешенным диаметром частиц ёср, позиционными диаметрами интегральной гранулометрической кривой ё95, d99, которым соответствуют ординаты гранулометрической кривой Р = 95% и 99%. Коэффициент однородности грунта ко = dср / d5, где d5 - средний диаметр крупных фракций, содержание которых в грунте составляет 5%, он определяет шероховатость русла.
При к0> 0,8 грунт считается однородным, при к0< 0,8 - неоднородным. Диаметры частиц наносов в поверхностном слое, их форма и гранулометрический состав грунта обусловливают значение шероховатости русла. При зернистой шероховатости, образованной потоком в русле из однородного несвязного грунта, среднюю высоту выступов шероховатости определяют по зависимости [3]: А = 0,7 dср; для неоднородных по крупности несвязных грунтов А = 0,7 d5.
Таблица 1 - Классификация наносов по диаметрам
№ Наименование фракций Крупность, мм
№
п/п
1 Глина 0,0015 * 0,005
2 Ил мелкий 0,005 * 0,015
3 Ил крупный 0,015 * 0,05
4 Песок мелкий 0,05 * 0,15
5 Песок средний 0,15 * 0,50
6 Песок крупный 0,50 * 1,5
7 Гравий мелкий 1,5 * 5,0
8 Гравий крупный 5,0 * 15
9 Галька мелкая 15 * 50
10 Галька 50 * 150
11 Валуны мелкие 150 * 500
12 Валуны средние 500 * 1000
13 Валуны крупные >1000
На основе экспериментальных и натурных исследований предложено большое количество формул для определения неразмывающих скоростей потока. Впервые определение этих скоростей (непередвигающей и срывающей) дано В.Н. Гончаровым [3]: непередвигающая скорость - наибольшая скорость потока, при которой зерна на дне еще не перемещаются; срывающая скорость - та наименьшая скорость потока, при которой уже происходит непрерывный срыв отдельных зерен на дне.
Нормативными документами рекомендована для использования формула Ц.Е. Мирцхулава для определения допускаемой неразмывающей донной скорости потока для несвязных грунтов:
у^ = 1,25
А н.доп 7
2Щ У (А - + 2С ],
0,4Арпр ун '
(4)
где d - диаметр частиц грунта, м; рг и р - плотность материала частиц грунта и воды, кг/м3, Снун - нормативная усталостная прочность несвязного грунта на разрыв, Па; m\, щ,
Ы - соответственно коэффициенты условий работы русла, перегрузки и однородности сил сцепления. Коэффициент перегрузки для частиц d < 1мм определяем по эмпирической формуле [5]:
Пр — 1 + т ■ Г..2 // >3ч-|1/3 , (5)
для частиц и > 1 мм.р
Вестник Дагестанского государственного технического университета. Технические науки. №2 (29), 2013
1
0,3 + [у2/( 3)]
и > 1 мм пр — 4.
Результаты многочисленных исследований показали, что подвижность наносов хорошо характеризуется отношением средней скорости течения к гидравлической крупности наносов .
Режим движения наносов, их перемещение в потоках, взаимодействие текущей воды и русла в весьма значительной степени зависят от гидравлических элементов потока: от расхода потока, удельного расхода, ширины потока, глубины потока, средней скорости потока и распределения скоростей по вертикали, от характеристик турбулентности потока.
Начало и интенсивность движения придонных наносов зависят от соотношения допускаемой неразмывающей донной скорости УАндоп и осредненной во времени
фактической донной скорости ид . Последняя определяется в зависимости от закона распределения местных скоростей по глубине потока. При логарифмическом законе распределения местных скоростей [3]:
1^16,7 у / А +1)
V ~ 1§(6,15А / А) . (6)
Донная скорость потока ид на высоте выступов шероховатости А = 0,7и при однородном составе русловых отложений определяется по формуле:
и А —
1^8,8 А / и) . (7)
Согласно экспериментальным данным, мгновенные скорости течения подчинятся нормальному закону распределения[8]:
\ 1 Г (и - и)2
(и) = -рт~ехР [ —; (8)
/у л 1 Г (и)\
f (о) = еХр[-^2], (9)
где f (и) и f (о) - одномерные плотности распределения горизонтальной и(?) и
вертикальной о(?) составляющих мгновенной скорости; Си — \(и ) и (Уи — \(о ) -
средние квадратические отклонения этих составляющих; и — и — и, и = и -пульсационные скорости.
Для расчета транспорта руслоформирующих наносов необходимо знать структуру, размер, частоту и длительность крупномасштабных турбулентных вихревых образований в потоке. Статистические характеристики турбулентных образований получаем на основе анализа корреляционных и спектральных функций Эйлеровых и Лагранжевых скоростей течения.
По результатам многолетних исследований процессов осаждения твердых частиц в турбулентных потоках, проведенных в МГУ, как в натурных, так и лабораторных условиях, получено, что одним из важных параметров, характеризующих взвешенное состояние частиц является вертикальная скорость осаждения твердых частиц.
Вестник Дагестанского государственного технического университета. Технические науки. №2 (29), 2013
Средняя частота пульсаций Эйлеровой скорости уи зависит от средней скорости и глубины потока и определяется через опытное значение числа Струхаля Sh = 0,73:
уи = 0,73V / К.
Распределение максимальных значений мгновенной скорости течения описывается законом Релея [8]:
и — и (и — и )2 _
/т (ит) =еХР[--т 2 ] , ит > и , (10)
а 2а
где и - осредненная во времени местная скорость, ит - максимальное значение
продольной составляющей скорости.
С учетом закона Релея, выражение для математического ожидания максимальных значений продольной составляющей донной скорости течения, определяемое по формуле [8]:
и
т
= | ит ■ /т (ит №т , (11)
0
принимает вид:
йАт = и + 7^/2сиА; (12)
Аналогично, выражение математического ожидания максимальных значений вертикальной составляющей скорости в придонной области имеет вид:
»т =У[Л/2С0 , (13)
где сиА и с» - соответственно средние квадратические отклонения продольной и
вертикальной составляющих донной скорости, согласно экспериментальным данным равные:
= 0,33ид , с»= 0,15ид (14)
В результате натурных исследований на реках бассейнов Сев. Двины, Волги, Оби, Лены и Яны [9] построены графики зависимости числа Лохтина (Л) и морфометрического параметра Глушкова (Г) от удельной величины расхода влекомых наносов. Увеличение устойчивости русла соответствует снижению подвижности наносов, уменьшению расхода влекомых наносов и их среднегодовой концентрации, в тоже время величина параметра Глушкова, отражающая соотношение ширины ^р и глубиныА русла прямо пропорциональна удельной концентрации влекомых наносов в потоке.
Сложность механизма взаимодействия турбулентного потока с несвязным грунтом увеличивается в связи с тем, что транспорт наносов почти всегда сопровождается образованием на дне подвижных форм - рифелей и гряд, изменяющих кинематическую структуру потока, и следовательно, влияющих на характер и количественные закономерности движения твердых частиц. Экспериментально установлено, что расход влекомых наносов существенно меняется в зависимости от формы гряд (рифели, дюны, заструги и т.д.).
Формирование гряд, их тип, а также статистические характеристики их размеров определяются гидравлическим режимом, кинематической структурой потока и физико-механическими свойствами донных отложений. Отсутствие единой точки зрения в вопросе механизма образования и развития различных типов гряд, а также несовершенство теоретических моделей этого процесса обусловили широкое использование в русловых процессах эмпирических критериев и формул для определения типа и параметров гряд, полученных на основе опытных и натурных исследований.
В 2008 г. ГГИ[6] были выполнены натурные исследования на различных участках нижней части реки Кубань, в ходе которых измерялись геометрические характеристики гряд, глубины и скорости потока над грядами и гранулометрический состав наносов, слагающих эти гряды, в результате чего построен график зависимости и получена формула полного расхода наносов от расхода воды:
Qj = 0.0000002i?3 - 0.0002(?s + 0.53Q - 39 (15)
Учитывать все многообразие влияющих на транспорт наносов факторов практически невозможно, поэтому в существующих методах расчета ограничиваются основными параметрами.
Выводы
Обобщение экспериментальных данных ряда исследователей, анализ литературных источников привели к выводу, что процесс транспорта руслоформирующих наносов является одним из сложных физических процессов, обусловленных множеством факторов: гидромеханических характеристик водного потока и физико-механических характеристик наносов, что является причиной отсутствия единой теории этого процесса, сложности многочисленных расчетных методов определения расхода наносов и большого расхождения в результатах расчетов по ним.
На современном этапе развития гидравлики и инженерной гидрологии одним из оптимальных путей получения представлений о качественном и количественных закономерностях этого процесса является использование теории вероятностей, теории случайных процессов и методов многомерного статистического анализа.
Библиографический список:
1. Великанов М.А. Русловой процесс. М.: 1958. - 395 с.
2. Великанов М.А. Движение наносов. М.: МРФ СССР, 1948.
3. Гончаров В.Н. Движение наносов. - Л.: Гидрометеоиздат, 1938.
4. Гришин Н.Н. Механика придонных наносов. М.: Наука, 1982. - 160 с.
5. Караушев А.В. Теория и методы расчета речных наносов. Л.: Гидрометеоиздат, 1977. - 272 с.
6. Копалиани З.Д., Жук М.М. К проблеме оценки стока руслового материала в реках // Журнал университета водных коммуникаций. -2010.-С.74-89.
7. Лапшенков В.С. Русловая гидротехника. - Новочеркасск: Изд-во НГМА, 1999. - 408 с.
8. Магомедова А.В. Эрозионные процессы в руслах рек и каналов. - М.: ВЗПИ, 1990. -203 с.
9. Чалов Р.С. Сток наносов, транспортирующая способность потоков и их роль в формировании речных русел. // География и природные ресурсы.-2011. - №3.-С.20-27.
10. EinsteinH.A. Thebed-loadfunctionforsedimenttransportationinopenchannelflows. - U.S. Dept. AgricultureSoil Cons.Serv. Techn. Bull. - 1950.
