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

Размещено на http://www.allbest.ru//

Размещено на http://www.allbest.ru//

Технический университет «Gh. Asachi» Ясы

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

Андриан Прикоп

1. Введение

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

Существует несколько методов для автоматического контроля оросительных каналов. Преимущества и недостатки этих методов были обсуждены разными учёными /Buyalski в 1991 году; Goussard в 1993 году/.

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

2. Основные уравнения

2.1 Уравнения непостоянного движения в открытых руслах (уравнения Сен-Венана)

Уравнения Сен-Венана описывают движение воды в открытых руслах. Они состоят из уравнения неразрывности и уравнения движения (динамики). Попа в 1997 г. показал, что дифференциальные уравнения непостоянного движения в открытых руслах могут иметь вид:

. (1)

, (2)

где А - площадь поперечного сечения; Т - время; Q - расход; х - горизонтальная координата в направлении течения; g - гравитационное ускорение; h - глубина воды; Jf - продольный уклон дна канала; Sf - уклон трения; U - средняя скорость воды.

2.2 Уравнение течения из-под затвора

Истечение из-под вертикальных затворов при условии свободного и затопленного течения может быть описано следующими уравнениями /Киселев 1988/:

; (3)

, (4)

где ц - коэффициент который отражает потери напора (ц = f(Fr); Fr - число Фруда); е - коэффициент сжатия (е = f(a/Ham)); a - открытие затвора; Ham - глубина воды в верхнем бьефе; Hav - глубина воды в нижнем бьефе; b -пролет отверстия перекрываемого затвором; м - коэффициент расхода.

Метод интегрирования уравнений течения

Уравнения (1) и (2) не могут быть решены аналитическим методом, но они могут быть решены методом последовательного приближения. Расчетная схема показана на рис. 1: существует некоторая сеть в точке P, проектируемая в координатах x и t; x - пространственный шаг, t - шаг времени между точками сети.

Знакомые кондиции UL, hL и UR, hR на время t = t, использованы для решения неизвестных UP, hP, после шага времени Дt, то есть на время t = t + Дt.

Этот метод позволяет определить глубину воды hP в определённое время Дt с помощью уравнения, приведённого выше, после чего определяется скорость Up в точке P.

Рис. 1. Схема в конечных разностях использованная в расчетной схеме

Для стабильности алгоритма необходимо придерживаться критерия Куранта

, (5)

где с - скорость распространения.

Контрольный алгоритм

На рисунке 2 показана генеральная схема одного бьефированного канала с вертикальными затворами с электрическим приводом. Номер бьефов (m+1), а номер затворов (m).

Течение в канале описывается уравнением Сен-Венана (1, 2). Расход, вытекающий из-под вертикального затвора при условии свободного и затопленного течения, определяется уравнениями (3, 4). Коэффициент расхода для условия затопленного течения равен /Киселев 1988/

, (6)

где n - взаимоотношение между площадью отверстия и площадью поперечного сечения вверх по течению от него; m - взаимоотношение между площадью отверстия и площадью поперечного сечения вниз по течению от него; ж - коэффициент сопротивления (ж 0,06); - коэффициент сжатия.

Размещено на http://www.allbest.ru//

Размещено на http://www.allbest.ru//

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

затвор бьеф канал оросительный

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

К исходному моменту (t = 0) постоянное течение по каналу постепенно варьирует Q = Q(t = 0)).

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

3.Пример вычисления

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

Физические параметры канала: длина канала - 9000 м (бьеф1 - 2000 м; бьеф2 - 5000 м; бьеф3 - 2000 м); ширина по дну - 2 м; продольный уклон дна канала - 0,0005; коэффициент заложения откоса - 1,5; коэффициент шероховатости по Мэнинга - 0,018; ширина затвора - 2 м. В канале нет боковых входов и выходов.

Параметры контрольного алгоритма: уровень воды (заданный) вниз по течению от затвора 1: H1 = 0,8 м; уровень воды (заданный) вниз по течению от затвора 2: H2 = 0,5 м; временной шаг: t = 0,01 с; пространственный шаг: x = 10 м; число узлов = 903; скорость поднятия затворов = 0,05 t , м/с; ширина мёртвой полосы = 0,04 м ( 2 cм); начальные условия: Q (t = 0) = 1,2 м3/с; общее время = 7000 с.

Размещено на http://www.allbest.ru//

Размещено на http://www.allbest.ru//

Рис.3. Канал с автоматическим контролем вниз по течению

Рис. 4. Гидрограф притока в конце канала вверх по течению

Математическая модель была разработана с помощью программы Матлаб для анализа поведения канала под автоматическим контролем вниз по течению (рис. 5,6).

Выводы

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

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

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

Размещено на http://www.allbest.ru//

Размещено на http://www.allbest.ru//

Рис. 5. Начальные условия - движение постоянное постепенно варьирует линия свободной площади воды (t = 0); hn = 0,6 м; hcr = 0,32 м

Рис. 6. Линия свободной площади воды: Т = 7000 с

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

Вuyalski, C. P., Ehler, D. G., Falvey, H. T., Rogers , D. C. and Serfozo, E. A. (1991). Canal systems automation manual. Vol. I. U.S. Bureau of Reclamation, Denver, Colorado.

Goussard, J. (1993). Automation of canal irrigation systems. ICID, New Delhi, ISBN: 81-85068-44-5.

Graf, W., H., Altinakar, M., S. (1993) Hidraulique fluviale; Tome I. Press Polytechniques et Universitairea de Romandes, Laussane, ISBN: 2-88074-261-7.

Kiselev P. G. (1988) Оndreptar pentru calcule hidraulice. Editura Tehnicг, Bucuresti.

Marcoie, N., Pricop, A., Cismaru, C. (2002). Model єi program оn Matlab pentru analiza hidraulicг a canalelor de irigaюii echipate cu regulatoare de nivel. Lucrгrile Sesiunii єtiinюifice omagiale, 17-18 mai, U.S.A.M.V. Bucureєti, Editura BREN, ISBN: 973-648-020-8, p. 403-414.

Popa, R. (1997). Elemente de hidrodinamica rвurilor. EDP, Bucuresti, 1997, ISBN: 973-30-5705-3.

Размещено на Allbest.ru