УДК 626.824-52
05.00.00 Технические науки
ГИДРАВЛИЧЕСКИЙ РАСЧЕТ ЛИНЕЙНО РАСПРЕДЕЛЕННЫХ ПОПУСКОВ СБРОСНЫХ РАСХОДОВ ВОДЫ ИЗ ВОДОХРАНИЛИЩА
Иваненко Юрий Георгиевич д-р техн. наук, профессор кафедры «Водоснабжение и использование водных ресурсов»
SPIN-код автора: 2715-4838 E-mail: [email protected] Новочеркасский инженерно-мелиоративный институт имени А.К. Кортунова ФГБОУ ВО Донской ГА У, Новочеркасск, Россия
Ткачёв Александр Александрович д-р техн. наук, профессор кафедры «Гидротехническое строительство» SPIN-код автора: 4732-0239 E-mail: [email protected] Новочеркасский инженерно-мелиоративный институт имени А.К. Кортунова ФГБОУ ВО Донской ГА У, Новочеркасск, Россия
Гурин Константин Георгиевич к.т.н., профессор кафедры «Водоснабжение и использование водных ресурсов» SPIN-код автора: 9483-0262 E-mail: [email protected] Новочеркасский инженерно-мелиоративный институт имени А.К. Кортунова ФГБОУ ВО Донской ГА У, Новочеркасск, Россия
Иваненко Дмитрий Юрьевич аспирант
E-mail: [email protected] Новочеркасский инженерно-мелиоративный институт имени А.К. Кортунова ФГБОУ ВО Донской ГА У, Новочеркасск, Россия
Цель настоящей работы заключается в разработке методов гидравлического расчета экстремальных расходов и глубин воды в критических створах водотоков, отводящих воду от водосбросных сооружений, для принятого закона регулирования линейно распределенных попусков сбросных расходов воды из водохранилищ. Задача регулирования линейно распределенных попусков сбросных расходов воды из водохранилищ основана на гидравлическом расчете процесса распространения и трансформации длинных волн, описываемых уравнениями Сен-Венана. Эти уравнения нелинейны и в общем случае не имеют точного решения. Для получения приближенных решений их линеаризуют. Установлено, что задача интегрирования этих уравнений предполагает
UDC 626.824-52 Technical Sciences
HYDRAULIC DESIGN OF LINEARLY DISTRIBUTED RELEASES OF STORED WATER DISCHARGES
Ivanenko Yury Georgievich
Dr.Sci.Tech., Professor of the chair "Water supply and water resources use" SPIN-code: 2715-4838 E-mail: [email protected]
Novocherkassk Engineering Meliorative Institute named after A.K.Kortunov, FSBEIHEDonskoi SAU, Novocherkassk, Russia
Tkachev Aleksandr Aleksandrovich Dr.Sci.Tech., Professor of the chair "Hydraulic Engineering construction" SPIN-code: 4732-0239 E-mail: [email protected]
Novocherkassk Engineering Meliorative Institute named after A.K.Kortunov, FSBEI HE Donskoi SAU, Novocherkassk, Russia
Gurin Konstantin Georgievich
Candidate of Tech.Sci, Professor of the chair "Water
supply and water resources use"
SPIN-code: 9483-0262
E-mail: gurin-knstantin@rambler. ru
Novocherkassk Engineering Meliorative Institute named
after A.K.Kortunov, FSBEI HE Donskoi SAU,
Novocherkassk, Russia
Ivanenko Dmitry Yurevich postgraduate
E-mail: [email protected]
Novocherkassk Engineering Meliorative Institute named after A.K.Kortunov, FSBEI HE Donskoi SAU, Novocherkassk, Russia
The purpose of the article is to develop methods for hydraulic design of extreme discharges and water depths in critical gage lines of waterways diverting water from spillways for the accepted law to regulate linearly distributed releases of stored water discharges. The task of control for linearly distributed releases of stored water discharges is based on hydraulic design for the process of propagation and transformation of long waves described with Saint Venant equations. These equations are nonlinear and have no accurate solution in general case. To get approximate solutions they are linearized. It is established that the task of integration for these equations supposes to get their preliminary solution in the form of complete integral. Applying the complete integral theory we get analytical solution that describe the process of transformation for one direction waves moving in
предварительное получение их общего решения в виде полного интеграла. С применением теории полного интеграла получены аналитические решения, которые описывают процесс трансформации волн одного направления, движущихся в бесконечно длинном призматическом канале полуограниченной протяженности с начальным равномерным режимом течения воды. Приведен пример расчета изменений расхода воды и глубины потока в двух фиксированных створах нижнего бьефа гидроузла на расстоянии 5000 и 10000 м от начального створа. Гидравлический расчет осуществлен по выведенным в работе аналитическим формулам и принятому в качестве образца классическому методу характеристик, который рассматривается как аналоговый. Сравнение результатов расчетов по двум методам позволило определить максимальную относительную погрешность, которая для расчетных расходов воды и глубин не превышает 3,5 %. Внедрение разработанного метода гидравлического расчета линейно распределенных попусков сбросных расходов на водотоках, отводящих воду от водосбросных сооружений, позволит оптимизировать холостые и нетехнологические сбросы воды из водохранилищ
Ключевые слова: ДИФФЕРЕНЦИАЛЬНОЕ УРАВНЕНИЕ, ЧАСТНЫЕ ПРОИЗВОДНЫЕ, ИНТЕГРИРОВАНИЕ, ГИДРАВЛИЧЕСКИЙ РАСЧЕТ, СБРОСНЫЕ РАСХОДЫ, ВОДОХРАНИЛИЩЕ, ГИДРОУЗЕЛ
РО!: 10.21515/1990-4665-125-013
Введение.
Гидравлический расчет расходов и глубин воды водотоков в нижних бьефах гидроузлов, осуществляющих попуски сбросных расходов воды из водохранилищ, представляет собой сложный процесс, основанный на применении закона регулирования режима водоподачи. Одновременное прямое измерение гидравлических характеристик в контролируемых створах водотоков затруднено даже при установившемся режиме течения воды. Из-за отсутствия проверенных и практически осуществимых методов измерения расходов и глубин воды в расчетных створах водотоков в настоящее время ощущается недостаток данных для уверенного прогнозирования экстремальных значений расходов и глубин воды в нижних бьефах гидроузлов [1-4].
infinitely long prism canal of semi-bounded extent with initial uniform regime of water flow. An example of changes in both water discharge and flow depths in two fixed downstream gage lines of hydraulic works at the distance of 5000 and 10000 m from the initial gage line is presented. Hydraulic design was carried out according to the derived analytical formulas accepted as an example classical method of characteristics that is considered as analogous. Comparison of the design results on two methods allowed to determine the maximum relative error that does not exceeds 3,5 % for the design discharges and water depths. Introduction of the developed method for hydraulic design of linearly distributed releases of discharges in waterways diverting water from spillways will make it possible to optimize idle and non-technological stored water discharges
Keywords: DIFFERENTIAL EQUATION, INDIVIDUAL DERIVATIVES, INTEGRATION, HYDRAULIC DESIGN, DISCHARGES, RESERVOIR, HYDRAULIC WORKS
Цель настоящей работы заключается в разработке методов гидравлического расчета экстремальных расходов и глубин воды в критических створах водотоков, отводящих воду от водосбросных сооружений, для принятого закона регулирования линейно распределенных попусков сбросных расходов воды из водохранилищ.
Постановка и решение задачи
Задача регулирования линейно распределенных попусков сбросных расходов воды из водохранилищ основана на гидравлическом расчете процесса распространения и трансформации длинных волн, описываемых уравнениями Сен-Венана [4]:
Э£+_ (и 2 _с 2)ВЩ_ _ ?Ю(,0 _и1_}=0, (1)
дt дх дх Сш Я
^+В Ш.=0, (2)
д х д?
где Q - расход воды, м/с; ? - время, с; и - средняя скорость течения воды в сечении, м/с; х - пространственная координата, м; С - скорость распространения начального возмущения, м/с; В - ширина водотока по урезу воды, м; Н - глубина потока, м; g - ускорение силы тяжести, м/с ;
ю - площадь поперечного сечения, м ; г0 - уклон дна водотока; Сш -коэффициент Шези, м0,5/с; Я - гидравлический радиус, м.
В качестве граничных условий на внешних границах рассматриваемого водотока применяются функции изменения расходов или глубин во времени в виде Q(t) и Н (?) или зависимости вида Q=/(Н). Во внутренних узлах рассматриваются условия равенства глубин, баланса расходов и др.
В качестве начального условия принимается состояние гидравлического режима, близкое к установившемуся.
Дифференциальные уравнения (1) и (2) - нелинейные и в общем
случае не имеют точного решения. Для получения приближенных решений эти уравнения линеаризуют [5-7].
Разложим значения расхода и глубины на составляющие в виде:
H = И + ДИ, Q=Qo + ДQ, (3)
где ДH - малое возмущение глубины потока, м; ДQ - малое возмущение
расхода воды, м3/с.
Уравнения возмущенного течения запишем в виде:
^+2и 0 М - (и 02 - Со2) ВО + 'Щ0 ДQ - (2 00 - '»С^ФВ ДH = 0, (4)
Э? Эх Эх U 0
В^М = 0, (5)
Э? Эх
где параметры с индексом «0» отвечают режиму невозмущенного течения воды. Их наименования и единицы измерения приведены выше;
Ф - параметр, определяемый элементами невозмущенного течения и зависящий от формы поперечного сечения русла.
Можно определить, в частности, для трапецеидального сечения [1]:
2лД+т2 Ь+2тИ.
Ф=(1+2 у)
0
(6)
Ь+2И 0л/1+т2 (Ь+тИ 0) И 0 где у - показатель степени в формуле Шези, зависящий от величины
коэффициента шероховатости и гидравлического радиуса; т - заложение откоса; Ь - ширина русла по дну.
В уравнениях (4), (5) параметры ДИ и ДQ носят название малых возмущений глубины и расхода и отвечают режиму возмущенного течения. Они характеризуют изменение глубины и расхода относительно начального положения.
Преобразуем дифференциальное уравнение (4) к виду:
^+2и -(и2 -С2)В0ЭДИ=-П^-(и ИДЯ], (7) Э? Эх Эх
где П - параметр, определяемый по зависимости:
(Р-^ _ У)
П= 3АН_; (8)
[ ВЩ - ^ 0М С°)]
Р=^ У=(2 g<o -С Ф), С0 =.
g— - параметры формы русла;
В0
тт ' I \ <^0 0 0 /■> - 0 -1
и 0 )
АQ и АН - осредненное значение возмущений по расходу и глубине
потока на участке дифференцирования.
Произведем подстановку вида:
К Р
АQ=(и0 М С0)В0АН_Пе_П? _Пе (и0±С0), (9)
где К, Р - постоянные;
е - основание натурального логарифма, е = 2,718; Система уравнений (5), (7) приводится к одному дифференциальному уравнению в частных производных первого порядка относительно функции DQ:
^+(и 0 м С0) М=Ке _П ?. (10)
д? дх
Общим решением дифференциального уравнения (10) является соотношение в виде полного интеграла:
АQ=Ф[ х _ (и0 + С0)? ] _ П е _П? + С, (11)
где функция Ф и постоянный параметр С определяются из граничных условий [8].
Применив теорию полного интеграла, будем искать решение, которое описывает процесс трансформации волны одного направления, движущейся в бесконечно длинном призматическом канале с начальным равномерным режимом течения воды [9, 10].
Пусть граничное условие в створе возмущения задается в виде:
Д0(0, 0=£г, (12)
где £ - постоянный параметр.
Подставив выражение (12) в уравнение полного интеграла (11), получим:
К -
П г
£ г=-Ф[(и0 М С)] -Пе + С.
(13)
В уравнениях (9)-(13) все обозначения и единицы измерения приведены выше.
Уравнение (13) продифференцируем по независимой переменной г:
-Ш
£ = -Ф (и 0 + С 0) + Ке" Из (14) найдем выражение для г :
1
г=-=1п П
£+Ф(и 0 М С0)
К
(14)
(15)
Исключив из двух соотношений: (13) и (15) переменную г, получим:
-11п
П
£ + Ф (и 0 М С0)
К
Ф (и 0 М С 0)
П
[£+ф (и+С0)]
П
1п
£ + Ф(и0 МС0)
К
(16)
+С.
Из соотношения (16) определяется значение постоянного параметра С. Подставив соотношение (16) в уравнение (11), найдем:
ДQ=Ф[х - (и0 + С0)г ] - К е-П г - £ 1п
П
П
£+Ф(и0 М С0)
К
Ф(и 0 М С0)
П
1п
£+Ф(и0 М С0) + [£+Ф(и0 М С0)]
К П
(17)
Продифференцируем соотношение (17) по параметру Ф и полученное выражение приравняем к нулю:
^ = [х_(и0 +С0)?]__ ^0 +С0)
д^ 1 0 ^ П[|+Ф(и0 +С0)]
_ (и 0 + С0) 1п П
| + Ф(и0 + С0)
К
Ф(и0 + С0)2 + (и 0 + С0) = 0
(18)
Ф=_К_еП ^^ J__I_. (19)
П[|+Ф(и 0 + С0)] П
Из выражения (18) определим параметр Ф:
^ х_(и0 +Ср) ? ] (и0 +С0)
и + С0) (и, + С0)
Подставив соотношение (19) в выражение (17), получим искомое решение в виде:
_[ х_(и +С0)? ] _
АQ=К еП^0+С0Г _ К е _П ? _х _(и0 + ^] (20)
П П (и +С0) ' '
Глубина воды определяется из соотношения (9) следующим образом:
В0 АН =_ К е _ х _ (и 0 + С0)? ] +
0 П(и 0 +С0) (иа +С0)2
П (21)
+ _ Р еи ±0>У
П(и 0 + С0) '
Связь между постоянными параметрами К и р находится из граничных условий:
АQ=0, АН=0 при ?=0, х=0. (22)
Подставляя выражение (22) в соотношения (20) и (21), найдем:
Р=_ К.
Соотношения (20) и (21) можно преобразовать к виду:
,[ х_(и 0 +Ср)? ]
ю=КеП (и0 +С) _КеП? _х_(и0 + Ср)?] (23)
П П (и, +С0) ' ' 7
-[ х-Ц+С0)1 ]
В ДЯ = _ К_ е - А X - (Цо м Ср)г ]
П(Ц м Со) (Ц м Со)
- П х
К_е[(Цо ±Со)]
(24)
П(Цо м Со)
Постоянный параметр К определяется из граничных условий:
ДО = о, ДЯ=о, I=Т, х=(Цо ±Со)Т, (25)
где Т - время распространения начального возмущения до расчетного створа х, вычисляется по формуле:
х
Т=---. (26)
(Ц ± Со)
Подставляем выражение (25) в соотношение (23) или (24) и находим:
Р=±-2ХПСоТ-. (27)
± 2ПСоТ _
(ио мСо)(^(ЦомСо) -е-ПТ)
Введем соотношение (27) в выражения (23) и (24). После чего можно определить:
П х-(Ц о+Со) г ] _
Х (е (иомСо) - е-П1 )
ДО =----[х - (и + С )1] ± 2ХС Т——=--——}-, (28)
(ТТ ИГ \ I 2ПСоТ _ ' к '
(ио + Со) (е±(Цо^ - е-ПТ)
П [ х-(ио мСо)г ] -П х X (е (ЦомСо) - е(Цо±Со))
Во ДЯ = - (и * ), [х - (и о М Со)1 ] ± 2^СТ(е 2ПСТ " е -1. (29)
(ио М Со) (е^- е-ПТ)
Соотношения (28) и (29) являются точными решениями уравнений (4), (7), удовлетворяющими также граничному условию (12). Сетка координат для полученных решений является прямоугольной. В соотношениях (28), (29) берется верхний знак, если направление движения возмущения совпадает с направлением течения воды в водотоке. В противном случае берется знак нижний. Параметр П определяется из уравнения (8), рассматриваемого в виде:
(Р^О-^ - т)
п = _в°о(н-Но)-, (30)
[ (д- бо) - (ио + Со)] Во( Н - Но) У о
где Q и Н - средние значения расхода и глубины на расчетном участке, которые рассчитываются следующим образом:
Q=-
п^к + Qо
2
Н=Нк + Н о
2
где Qk - расход воды, определяемый из краевого условия (13) при значении t = Ц, м/с; Нк - глубина воды, соответствующая расходу, определяемому из краевого условия (13) при значении t = tk, м; Q0 -начальное значение расхода на расчетном участке, м /с; Н0 - начальное значение глубины на расчетном участке, м.
Внедрение и оценка эффективности
Аналитические решения (28) и (29), полученные из предположения изменения расхода в створе возмущения по линейному закону (12), могут быть обобщены для реальных гидрографов попусков. В этом случае кривую гидрографа с любой степенью точности заменяют линейно аппроксимирующей ломаной линией. Зная решение нестационарной задачи для одной ступени, можно получить суммарное решение для всех ступеней.
Разработанный метод гидравлического расчета экстремальных расходов и глубин воды в критических створах водотоков, отводящих воду от водосбросных сооружений, для принятого закона регулирования линейно распределенных попусков сбросных расходов воды из водохранилищ обеспечивает:
- универсальность, заключающуюся в возможности определения
расходов и глубин воды в фиксированных створах при любых режимах неустановившегося, медленно изменяющегося режима течения воды;
- возможность повышения точности определения расходов и глубин воды. Это достигается использованием в качестве расчетных формул соотношений, являющихся точными аналитическими решениями уравнений неустановившегося, плавно изменяющегося режима течения воды.
Внедрение разработанного метода на водотоках с применением аналитических решений дифференциальных уравнений Сен-Венана позволит оптимизировать процессы водоизмерения и минимизировать холостые и нетехнологические сбросы воды.
Рассмотрим пример расчета по разработанному методу. При сбросе воды из водохранилища изменение расхода в начальном створе (х = 0) нижнего бьефа гидроузла подчиняется закону (12):
Q=520+0,05^ 0£t<3600 с.
Рассчитаем изменения расхода воды и глубины потока в фиксированных створах х = 5000 м их = 10000 м, если в начальный момент в призматическом русле реки на расчетном участке наблюдался режим течения, близкий равномерному. Морфометрические характеристики и гидравлические параметры на расчетном участке реки, соответствующие равномерному режиму течения, следующие: Q = 520
м3/с; и0 = 0,912 м/с; Н0 = 3,799 м; т = 0; Ь = 150,0 м; п = 0,02; /0 =
0,00006; у = 1/6.
По приведенным выше формулам определим параметры С0, в, у, П, Ф:
Со =
С0 =^1 9,81x3,799=6,687 м/с,
Во
ß = ß =
2x9,81x0,00006
U
0,912
=0,001
Ф=(1+21)[
2
150
]=-0,334,
6 (150+2x3,799) 150x3,799 у=(2 gi0 -10СЦ Ф), у=(2x9,81x0,00006+0,00006х6,6872 х0,334) = 0,00192,
П=0,000017 .
По формуле (26) вычислим время распространения начального возмущения до расчетного створа:
Г = 10000 =1425 с.
0,912 + 6,687
Изменения расходов воды и глубин потока во времени в фиксированных створах рассчитываются по формулам (28) и (29). Результаты расчетов приведены на рисунках 1 и 2.
Q. м3/с
700 675 650 625 600
550
500 4 "5
425
400
Q=f(T)
Тс
0 400 800 1200 1600 2000 2400 2800 3200 3600 4000
» Расчет [Ю аналитическим формулам, х 10000м —■—Расчет [Ю аналитическим формулами 5000м —*—Расчет [Ю методу характеристик, х 10000м —•— Расчет [Ю методу характеристик, х 5000 м Рис. 1. График изменения расходов воды в зависимости от времени в фиксированных створах
Н, М Н=1(Т)
О 400 800 1200 1600 2000 2400 2800 3200 3600 4000
« Расчет по аналитическим формулам, х 10000м
■ Расчет по аналитическим формулам, х 5000м
* Расчет по методу характеристик, х 10000м
» Расчет по методу характеристик, х 5000 м
Рис. 2. График изменения глубин воды в зависимости от времени в фиксированных створах
Расчеты выполнены по выведенным в работе аналитическим формулам (28), (29) и методу характеристик, который рассматривается как аналоговый. Сравнение результатов расчетов по двум методам позволило определить максимальную относительную погрешность, которая для расчетных расходов воды и глубин не превышает 3,5 %.
Заключение
1. В результате исследований разработан метод гидравлического расчета экстремальных расходов и глубин воды в критических створах водотоков, отводящих воду от водосбросных сооружений, для принятого закона регулирования линейно распределенных попусков сбросных расходов воды из водохранилищ с использованием аналитических решений дифференциальных уравнений Сен-Венана.
2. По полученным аналитическим решениям проведены расчеты изменений расхода воды и глубины потока в двух фиксированных створах нижнего бьефа гидроузла: на расстоянии 5000 и 10000 м от начального створа. При этом аналоговым методом был принят метод характеристик. Сравнение результатов расчетов по двум методам позволило определить максимальную относительную погрешность, которая для расчетных
расходов воды и глубин не превысила 3,5 %.
3. Внедрение разработанного метода гидравлического расчета линейно распределенных попусков сбросных расходов на водотоках, отводящих воду от водосбросных сооружений, позволит оптимизировать холостые и нетехнологические сбросы воды из водохранилищ.
Литература
1. Иваненко, Ю. Г. Гидравлические аспекты устойчивых водных потоков в неразмываемых и размываемых руслах / Ю. Г. Иваненко, А. А. Ткачев, А. Ю. Иваненко. - Новочеркасск: Лик, 2013. - 352 с.
2. Юрченко, И. Ф. Автоматизированное управление водораспределением на межхозяйственных оросительных системах / И. Ф. Юрченко, В. В. Трунин // Известия Нижневолжского агроуниверситетского комплекса: наука и высшее профессиональное образование. - 2012. - № 2. - С. 178-184.
3. Клишин, И. В. Современные системы управления водораспределением на оросительных системах / И. В. Клишин, В. И. Селюков // Мелиорация и водное хозяйство. - 2006. - № 6. - С. 23.
4. Управление водораспределением на открытых оросительных системах на основе гидрологической информации и агрометеопараметров / В. Н. Щедрин, С. М. Васильев, А. В. Акопян, В. В. Слабунов // Известия Нижневолжского агроуниверситетского комплекса: наука и высшее профессиональное образование. -2014. - № 2(34). - С. 152-158.
5. Неймарк, Ю. И. О допустимости линеаризации при исследовании устойчивости / Ю. И. Неймарк // ДАН СССР. - 1959. - Т. 127. - № 5. - С. 961-964.
6. Филиппов, Е. Г. Гидравлика гидрометрических сооружений для открытых потоков / Е. Г. Филлиппов. - Л.: Гидрометеоиздат, 1990. - 288 с.
7. Маковский, Э. Э. Автоматизированные автономные системы трансформации неравномерного стока / Э. Э. Маковский, В. В. Волкова. - Фрунзе: Илим, 1981. - 380 с.
8. Ольгаренко, И. В. Методы теории подобия для решения уравнений Сен-Венана при управлении водораспределением в оросительных системах / И. В. Ольгаренко, О. П. Кисаров, В. И. Ольгаренко // Вестник Российской академии сельскохозяйственных наук. -2012. - № 2. - С. 12.
9. Коваленко, П. И. Автоматизация мелиоративных систем / П. И. Коваленко. -М.: Колос, 1983. - 304 с.
10. Кюнж, Ж. А. Численные методы в задачах речной гидравлики: практическое применение / Ж. А. Кюнж, Ф. М. Холли, А. Вервей: [пер. с англ.]. - М.: Энергоатомиздат, 1985. - 256 с.
11. Dedrick, A. R. Automatic control of irrigation water delivery to and on-farm in open channels / A. R. Dedrick, D. D. Zimbelman // Principles of designing control systems for water resources and irrigation using modern techniques: Symposium, New Delhi. - New Delhi: ICID, 1981. - P. 113-128.
References
1. Ivanenko, Ju. G. Gidravlicheskie aspekty ustojchivyh vodnyh potokov v nerazmyvaemyh i razmyvaemyh ruslah / Ju. G. Ivanenko, A. A. Tkachev, A. Ju. Ivanenko. - Novocherkassk: Lik, 2013. - 352 s.
2. Jurchenko, I. F. Avtomatizirovannoe upravlenie vodoraspredeleniem na mezhhozjajstvennyh orositel'nyh sistemah / I. F. Jurchenko, V. V. Trunin // Izvestija Nizhnevolzhskogo agrouniversitetskogo kompleksa: nauka i vysshee professional'noe obrazovanie. - 2012. - № 2. - S. 178-184.
3. Klishin, I. V. Sovremennye sistemy upravlenija vodoraspredeleniem na orositel'nyh sistemah / I. V. Klishin, V. I. Seljukov // Melioracija i vodnoe hozjajstvo. - 2006. - № 6. - S. 23.
4. Upravlenie vodoraspredeleniem na otkrytyh orositel'nyh sistemah na osnove gidrologicheskoj informacii i agrometeoparametrov / V. N. Shhedrin, S. M. Vasil'ev, A. V. Akopjan, V. V. Slabunov // Izvestija Nizhnevolzhskogo agrouniversitetskogo kompleksa: nauka i vysshee professional'noe obrazovanie. - 2014. - № 2(34). - S. 152-158.
5. Nejmark, Ju. I. O dopustimosti linearizacii pri issledovanii ustojchivosti / Ju. I. Nejmark // DAN SSSR. - 1959. - T. 127. - № 5. - S. 961-964.
6. Filippov, E. G. Gidravlika gidrometricheskih sooruzhenij dlja otkrytyh potokov / E. G. Fillippov. - L.: Gidrometeoizdat, 1990. - 288 s.
7. Makovskij, Je. Je. Avtomatizirovannye avtonomnye sistemy transformacii neravnomernogo stoka / Je. Je. Makovskij, V. V. Volkova. - Frunze: Ilim, 1981. - 380 s.
8. Ol'garenko, I. V. Metody teorii podobija dlja reshenija uravnenij Sen-Venana pri upravlenii vodoraspredeleniem v orositel'nyh sistemah / I. V. Ol'garenko, O. P. Kisarov, V. I. Ol'garenko // Vestnik Rossijskoj akademii sel'skohozjajstvennyh nauk. -2012. - № 2. - S. 12.
9. Kovalenko, P. I. Avtomatizacija meliorativnyh sistem / P. I. Kovalenko. - M.: Kolos, 1983. - 304 s.
10. Kjunzh, Zh. A. Chislennye metody v zadachah rechnoj gidravliki: prakticheskoe primenenie / Zh. A. Kjunzh, F. M. Holli, A. Vervej: [per. s angl.]. - M.: Jenergoatomizdat, 1985. - 256 s.
11. Dedrick, A. R. Automatic control of irrigation water delivery to and on-farm in open channels / A. R. Dedrick, D. D. Zimbelman // Principles of designing control systems for water resources and irrigation using modern techniques: Symposium, New Delhi. - New Delhi: ICID, 1981. - P. 113-128.