Закономерности изменения гидравлических сопротивлений земляных русел каналов при эксплуатации

УДК 620.178.620.179

ЗАКОНОМЕРНОСТИ ИЗМЕНЕНИЯ ГИДРАВЛИЧЕСКИХ СОПРОТИВЛЕНИЙ ЗЕМЛЯНЫХ РУСЕЛ КАНАЛОВ ПРИ ЭКСПЛУАТАЦИИ

© 2011 г. Ю.М. Косиченко

Российский научно-исследовательский институт проблем мелиорации

The Russian Scientific Research Institute of Land Improvement Problems

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

Ключевые слова: гидравлические сопротивления; земляные русла; каналы; эмпирические зависимости.

Analysis for changes in hydraulic resistance of ground channels in operation at full-scale observational data channels for southern Russia was carried out. Typical regularities for its changes are established. Empirical dependences are received. The empirical relationships are associated with well-known formulas of other authors and field data.

Keywords: hydraulic resistance; ground channels; channels; empirical relationships.

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

Гидравлические сопротивления русел каналов непосредственно входят в формулу потерь напора на трение.

Изучению гидравлических сопротивлений земляных русел каналов посвящено наибольшее количество работ: В. С. Алтунина [1], Г. В. Железнякова [2], И. Ф. Карасева [3], В. Т. Чоу [4] и других [5 - 9]. Такой интерес к земляным руслам каналов обусловлен сложностью процессов, происходящих при движении потока: различные деформации русел, русловые явления, заиление, подмывы берегов, зарастание водной растительностью, обрушение откосов.

Однако имеющиеся данные исследований земляных русел в основном касаются крупных каналов средней Азии с характерными руслоформирующими-ся процессами [1, 3, 5]. Для каналов Юга России при их эксплуатации часто превалирующими факторами, влияющими на пропускную способность, а следовательно, и сопротивление, являются местные деформации откосов и реже дна русла, значительное заиление отдельных участков и интенсивное зарастание русел растительностью в летний период [3, 9, 10].

Для анализа закономерностей изменения гидравлических сопротивлений используем данные натурных исследований по 15 каналам в основном Северного Кавказа, за исключением некоторых данных по каналам Украины, полученным эксплуатационными организациями и рядом исследователей (И. Ф. Кара-сев, И. А. Долгушев, Ю. М. Косиченко, Ю. И. Иовчу, Н. Ю. Полякова, К. Г. Гурин). При этом отдельные каналы относились к категории крупных (Большой Ставропольский канал I очередь - БСК-1, Донской

магистральный канал, Невинномысский и Терско-Кумский каналы) с расходами от 75 до 200 м3/с, а остальные - к категории средних и малых (магистральные и распределительные) с расходами от 0,1 до 25,0 м3/с. Для подавляющего большинства каналов их русла сложены из связных грунтов (суглинки, супеси). Срок эксплуатации каналов составляет от 37 до 60 лет. Вычисления гидравлических характеристик потока в руслах каналов на основании натурных данных (Q, ш, и, I) производились по следующим формулам: - коэффициент шероховатости русла у 1/ J/2

и

- коэффициент Шези по Н. Н. Павловскому

С =1 Яу ; п

- коэффициент гидравлического трения

Х =

U 2

- число Рейнольдса

Re =

4UR

где Q - расход канала; ш - площадь живого сечения; и - средняя скорость потока; i - уклон дна русла; R -гидравлический радиус; J - уклон свободной поверхности; V - кинематическая вязкость воды (1-10"6 м2/с).

Результаты собранных натурных данных по каналам в земляном русле Юга России и вычислений их гидравлических характеристик приведены в обобщенном виде на рис. 1 и 2 в качестве графиков зависимостей п = п(0), Х = Х(0), п = п(сД/Я) и Х = Х(Re). При этом отметим, что последние две зависимости

n =

V

представлены с безразмерными аргументами, что принято в современной гидравлике и механике жидкости [7].

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

n 0,06

0.05 0,04 0,03

0,02 0,01

0,32 0.18

0,14 0.10

0,06

\

Л А \ о А Т -1 +-11 -2 *-12 -3 0-13 -4 »-14

\ • и  С * -5 Ф-1 -6 -7 -я 5

т Л • Л -9 -10

♦ - т т  — <v> --5 '

30

60

90

а

120

150

180

Q, м7с

+ V

\ о-1 + -2 • А-3 с? -11 -12 -13

т-4 т -5 ф  -6 «,-7 -14 -15

• • \ ♦-8 • -9 А-10

/ f м* ш — CJS- — ТО* -

30 60

90

120 150 180

Q, м7с

б

Рис. 1. Графики зависимостей п = п^) (а) и X = (б) для каналов в земляном русле: 1 - БСК-1; 2 - Донской МК; 3 -Азовский МК; 4 - Нижне-Донской МК; 5 - Пролетарский; 6 - Невинномысский; 7 - Терско-Кумский; 8 - Солдатский I; 9 - Сев. Донец-Донбасс; 10 - Теплушка; 11 - Солдатский II; 12 - Распределитель № 67; 13 - Баксан-Малка; 14 - Бг-Р-6;

15 - Днепр-Донбасс

Для земляных русел каналов (Азовский МК и Солдатовский I), где в процессе исследований было отмечено сильное зарастание водной растительностью (камышом) [9, 10], значения коэффициентов п достигают максимума 0,055 - 0,060, которые превышают для соответствующих расходов при нормальном состоянии земляных русел (Нижне-Донской МК, Бг-Р-6 и Солдатовский II) в 2,2 - 2,4 раза. Однако такое большое влияние зарастания каналов на коэффициен-

ты шероховатости характерно только для малых и средних каналов с расходами от 1,0 до 75,0 м3/с. Для крупных каналов при расходах более 75,0 м3/с влияние зарастания существенно уменьшается. Так, при расходах от 75,0 до 90,0 м3/с коэффициенты п максимально повышаются в 1,2 - 1,45 раза, а при Q от 90,0 до 200,0 м3/с - в 1,20 раз.

Кроме того, при одном и том же состоянии русел, в том числе нормальном, коэффициенты п с увеличением расходов снижаются. Это также подтверждается исследованиями на 1-й очереди Большого Ставропольского канала [10].

Аналогичные закономерности распределения также наблюдаются для гидравлических сопротивлений А в зависимости от расхода (рис. 1, б), числа Рей-нольдса (рис. 2, б) и для коэффициентов шероховатости п в зависимости от безразмерного параметра С/Л (рис. 2, а).

0,05 0,04 0,03 0.02 0,01

8

X 0,18

0,15

0,12

0,09

0,06

0,03

\ х о-1 +-11 е-2 4-12 а-3 0-13

^ \ ч т-4 т-14 0-5  -6

S ^ ♦-8 »-9 4-10

—— с _ э-о-Т,

10

12 14

а

16

1В

с! Г*

м

V 0—1 е-2 1-3 ч— * 9" 11 12 3

\ f • \ 4-5  6 е-7

• 1 + •-о • -9 а-10

ч ЧА

~ -D-Ch

-W3&—

Rex106 б

Рис. 2. График зависимостей п = п

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

1. Зона изменчивости коэффициентов шероховатости и гидравлических сопротивлений каналов в земляных руслах при расходах от 1,0 до 75,0 м3/с имеет расширенный диапазон их значений. Далее она существенно сужается и при расходах более 75,0 -120,0 м3/с практически остается постоянной.

2. Нижняя граница зоны шероховатости и гидравлических сопротивлений представляет собой плавную

n

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

3. Влияние зарастания русел каналов особенно сильно проявляется вначале зон п и X при относительно небольших расходах от 1,0 до 20,0 - 30,0 м3/с. Здесь значения коэффициентов шероховатости и гидравлических сопротивлений при зарастании превышают их по сравнению с нормальным состоянием русел каналов более чем в 2 раза.

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

Согласно принятым в гидравлике критериям, область сопротивления для открытых потоков является квадратичной, если удовлетворяется условие Прандт-ля - Никурадзе [8, 11]

Igl^l >1,83,

где и* =^JgRJ - средняя динамическая скорость потока, м/с; А = Кэ = (25п)6 - высота эквивалентной шероховатости, м; V - кинематическая вязкость воды (при t = 20 °С, V = 1 • 10-6 м2/с); Я - гидравлический радиус, м; J - уклон свободной поверхности.

Согласно проведенным расчетам, для крупного канала в земляном русле БСК-1 (при Q = 182,6 м3/с, и = 1,17 м/с, h = 5,5 м, J = i = 0,00015, п = 0,0235, Я = 3,43 м) получим и* = 0,0710 м/с, А = 0,044 м, ^ (и*А / V) = 3,50, для малого канала Бг-Р-6 (при Q = 1,93 м3/с, и = 0,39 м/с h = 1,39 м, J = i = 0,00014, п = 0,0259) имеем и* = 0,0327 м/с, А = 0,073 м, Я = 0,78 м ^ (и*А / V) = 3,38.

Таким образом, для указанных каналов по условию (1) должна наблюдаться квадратичная область сопротивления. Аналогичная зона сопротивлений получается и для других каналов в земляном русле, данные по которым представлены на графиках п и X (см. рис. 1 и 2).

Однако, как показано в работе О. М. Айвазяна [8], для течений в земляных каналах в связных грунтах и в песчаном ложе свойственна не квадратичная, как до сих пор считалось, а смешанная (переходная) зона сопротивления. В этой зоне кривые сопротивлений в логарифмических координатах характеризуются восходящими или нисходящими кривыми связи с числами Рейнольдса ^е). В квадратичной зоне кривые сопротивлений приобретают вид прямых горизонтальных линий относительно оси абсцисс, не зависящих от числа Re.

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

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

стоянии (рис. 3). Анализ этой графической зависимости показывает, что в начальной части на значительном протяжении наблюдается осредненная нисходящая кривая связи, которая при некотором предельном значении ^ Re = 7,0 практически становится горизонтальной. Подобные кривые сопротивлений также были получены И. Ф. Карасевым [3] для ряда каналов в земляных руслах и руслах рек.

1,85 1.75

1.65 1,55 1,45 1,35

л ч V о- -1

r ч Ф- -2 5 "7

о т- -14

Cl г» С Ô о

о -t- I 9 0

6,0 6,1 6,2 6,3 6,4 6,5 6,6 6,7 6,8 6,9 7,0

lgRe

Рис. 3. График зависимости коэффициентов сопротивления X от Re для земляных русел каналов при нормальном их состоянии

Считая предельным значением между переходной и квадратичной областью сопротивления lg Re = 7,0, где использованы данные канала БСК-1 (n = 0,0226, R = 3,43, А = 0,0329 м, i = 0,00015), получим U = = 0,0710 м/с, а критерий (1) составит

-g ( ^ 'g ( ^^Г^ ] = ^

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

-g ( >3,35.

Используя ранее представленные натурные данные по каналам в земляных руслах (см. рис. 1 и 2), найдем эмпирические зависимости при нормальном их состоянии (без зарастания) и при зарастании. Обработка данных и получение эмпирических зависимостей проводились с использованием пакета программ Microsoft Excel в виде степенной и полиноминальной связи.

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

4);

-7 ГУ2 1 1 Л-4 r\ I л Г\Л 1 Л i D2

n = 0,027g-°,°4

n = 2 • 10-7 Q2 -1 • 10-4Q + 0,0434 (R2 = 0,765) ;

)036-^ + 0,0651 (R2 = 0,827) ; (3) Vg

n = 0,00006

2

1 C Л

.fg.

X = 0,063Q-0,15 (R2 = 0,759) ; X = 2,707 Re-0,28 (R2 = 0,846),

v

где п - коэффициент шероховатости русла; Q - расход канала м3/с; С - коэффициент Шези, м0'5/с; X - коэффициент гидравлических сопротивлений (коэффициент Дарси); Re - число Рейнольдса.

Для практического использования могут быть рекомендованы лишь более точные формулы (2), (3) и (4), имеющие достаточно высокую достоверность аппроксимации > 0,80).

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

п = 0,422

Г л-1,121

1 С л

.я.

Х = 495,2 • Re"

^2 = 0,831); ,45 (R2 = 0,895) .

С целью сопоставления результатов определения коэффициентов гидравлического сопротивления каналов в земляном русле при нормальном их состоянии и при зарастании водной растительностью по рекомендуемым формулам (4) и (5) и известным формулам других авторов (И.Ф. Карасева, О.М. Айвазяна, Ю.А. Ибад-Заде, Марки, Штриклера, В. С. Боровкова и др.) проведем соответствующие расчеты и сравним с натурными данными.

Результаты расчетов коэффициентов сопротивлений X по формулам различных авторов и их отклонений от натурных данных 5Х, %, приведены в таблице.

Полученные результаты расчетов и их сопоставление с натурными данными для незаросших русел позволяет заключить, что наиболее близкие значения коэффициента X получены по формулам И.Ф. Карасева, Ю.А. Ибад-Заде и (4) (отклонение от натурных данных по ним не превышает ±11 %). По другим формулам (О.М. Айвазяна, Марки и Штриклера) отклонение от натурных данных составляет до 40 - 90 % и более. Существенное отклонение X по формуле О.М. Айвазяна от результатов натурных исследований объясняется тем, что в этой зависимости учтены данные только по каналам Средней Азии. Большое отклонение по формулам зарубежных авторов (Марки и Штриклера), на наш взгляд, объясняется недостаточным учетом в последних гидравлических условий и характеристик.

Анализ данных расчетов коэффициента X для заросших русел земляных каналов свидетельствет о большом расхождении результатов с натурными данными, за исключением формулы (5) (- 6,6 %) и И.Ф. Карасева (20,9 %). Для остальных зависимостей расхождение составляет до 85 - 450 %, что не позволяет их рекомендовать к применению для частично заросших русел каналов в условиях эксплуатации каналов Северного Кавказа.

Сопоставление расчетных значений коэффициентов гидравлических сопротивлений X по формулам различных авторов и натурных данных для каналов в земляных руслах

Канал Расчетные данные X / 5*,, % по формулам Натурные данные X

И.Ф. Карасева О.М. Айвазяна Марки Штриклера Ю.А. Ибад-Заде

Для незаросших русел

Большой Ставропольский канал (БСК-1) 0,0278* 10,2 0,0203 -19,3 0,0350 39,0 0,0970 285,0 0,0224 -11,3 0,0263 4,3 0,0252

Распределительный канал Бг-Р-6 (ПК-89) 0,0596 11,6 0,0342 -93,6 0,0727 36,1 0,076 43,4 0,0596 10,9 0,0536 0,4 0,0534

Для заросших русел

Канал И.Ф. Карасева А.Д. Асановой В.С. Боровкова Н.И. Турян-ской А. Найта

Распределительный канал Бг-Р-6 (ПК-82) 0,0473 20,9 0,327 446,8 0,110 83,9 0,13 122,4 0,314 425,0 0,058 -6,6 0,0597

Расчеты выполним для следующих каналов и исходных данных: канал БСК-1 в нормальном состоянии, Q = 186,15 м3/с, i = 0,00015, Re = 15,43-106, R = = hср = 3,2 м, п = 0,0226, X = 0,0252; распределительный канал Багаевско-Садковской ОС Бг-Р-6 (ПК-89) без зарастания, Q = 1,94 м3/с, i = 0,000145, Re = = 1,21 -106, п = 0,0249, X = 0,05, R = 0,75 м; для участка Бг-Р-6 (ПК-82) при зарастании Q = 2,34 м3/с, Re = = 1,32-106, п = 0,0268, X = 0,0598, i = 0,000145, R = = 0,83 м.

Выводы

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

такое влияние незначительно и в основном обусловлено деформациями русла.

2. Проведенный анализ натурных данных коэффициентов шероховатости и гидравлических сопротивлений земляных русел каналов Юга России позволил установить, что коэффициенты n земляных русел при нормальном их состоянии изменяются от 0,0210 до 0,0260, а при зарастании - от 0,0250 до 0,060, соответственно гидравлические коэффициенты X - от 0,0250 до 0,050 и от 0,0350 до 0,230.

3. Установлено, что для подавляющего количества земляных русел каналов Юга России наблюдается переходная (доквадратичная) область сопротивления, для которой получено уточненное значение критерия перехода к квадратичной области по условию Пранд-таля - Никурадзе. Очевидно, что для русел с переходной областью сопротивления не могут быть применимы широко известные в гидравлике формулы квадратичного сопротивления Н.Н. Павловского, И.И. Агро-скина, Маннинга и других авторов.

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

пользования для гидравлических расчетов каналов, аналогичных условиям Юга России.

Литература

1. Алтунин В. С. Мелиоративные каналы в земляных руслах. М., 1979. 255 с.

2. Железняков Г.В. Пропускная способность русел каналов и

рек. Л., 1981. 311 с.

3. Карасев И.Ф. Русловые процессы при переброске стока. Л., 1975. 288 с.

4. Чоу В.Т. Гидравлика открытых каналов : перевод с англ., М., 1969. 464 с.

5. Абальянц С.Х. Устойчивые и переходные режимы в искусственных руслах. Л., 1981. 240 с.

6. Ибад-Заде Ю.А. Водопроводные каналы. М., 1975. 192 с.

7. Айвазян О.М. Основы гидравлики равномерных течений / НИЦ «Регулярная и хаотическая динамика»; Институт компьютерных исследований. М.; Ижевск, 2006. 152 с.

8. Айвазян О.М. Зона гидравлического сопротивления земляных каналов // Гидротехническое строительство. 1987. № 11. С. 54 - 58.

9. Долгушев И.А. Повышение эксплуатационной надежности

оросительных каналов. М., 1975. 136 с.

10. Косиченко Ю.М. Каналы переброски стока России. / НГМА. Новочеркасск, 2004. 470 с.

11. Зегжда А.П. Гидравлические потери на трение в каналах и трубопроводах. Л.; М., 1975. 278 с.

Поступила в редакцию 29 марта 2011 г.

Косиченко Юрий Михайлович - д-р техн. наук, профессор, заместитель директора по науке ФГНУ «РосНИ-ИПМ». Тел. (8635)26-51-11.

Kosichenko Yuri Mihailovich - Doctor of Technical Sciences, professor, Deputy Director for Science FSSE «RSRILIP». Ph. (8635)26-51-11._