UPDATE OF HYDRAULIC CALCULATION DURING SYSTEM CONSTRUCTION DRIP IRRIGATION
Abstract and keywords
Abstract (English):
In this article, the hydraulic calculation was verified by analyzing the calculations of pressure loss with the data of field studies. As a result of the calculations performed, it is possible to more accurately determine the optimal parameters that ensure the uniform distribution of water in the irrigated area with the minimization of losses in the system. Thanks to the conducted studies of the dependence for the hydraulic calculation of drip irrigation systems, it is possible to adjust and set the general aspects that ensure the effective distribution of water in the irrigated area and reduce losses in the system. To ensure constant pressure during operation, it is recommended to determine the pressure loss in pipelines, connections, and drip emitters. To grow crops in favorable climatic conditions, it is necessary to determine the required water flow rate. Proper selection of equipment for different types of drip emitters and pumps will ensure the system's performance. It is recommended to develop a layout plan for the location of pipelines and drip emitters in advance, in order to minimize the cost of materials and installation, while maintaining proper water distribution. When developing the layout plan for the location of pipelines and drip emitters, it is necessary to analyze the impact of external environmental factors, such as temperature, humidity, soil type, and the effectiveness of the chosen drip irrigation system. When using software for modeling hydraulic processes, it is important to conduct a trial test of different system scenarios in advance. It is necessary to assess the initial costs of installing and operating the drip irrigation system compared to other irrigation methods.

Keywords:
drip irrigation, hydraulic calculation, water flow, pipeline throughput, pressure and pressure losses in the system, hydraulic losses, pipeline pressure, transit and travel flow
Text
Text (RU) (PDF): Read Download
Letter of recommendation (RU) (PDF): Read Download

Введение. Капельное орошение позволило повысить продуктивность выращиваемых культур, тем самым обосновав свою эффективность применения в агропроизводстве, что дало ей занять доминирующее место, среди технологий полива точного земледелия. Данная технология позволяет сэкономить около 60-80% воды по сравнению с традиционными системами орошения. Однако, эффективность эксплуатации систем капельного орошения зависит от совокупности различных факторов, в том числе от гидравлических расчётов. Выполненные в советское время натурные измерения потерь напора в участковых и поливных трубопроводах Приднестровской Молдавской Республики на Дубоссарской, Рыбницкой и Вулканештской систем капельного орошения, показало хорошее совпадение при сравнении экспериментальных результатов с расчётными. Разница суммарных потерь напора не превышала 10%, а разница абсолютной величины напора 5%. Примеры данного расчёта подтвердили актуальность применения для участкового и поливного трубопроводов.
Цель исследования. Верификация полученных данных в ходе проведения натурных испытаний полученных в Приднестровье Молдавской Республике в отношении зависимости для гидравлического расчёта при строительстве систем капельного орошения при определении оптимальных параметров, таких как: диаметр трубопроводов, расход и давление воды, потери напора, а также режимы водоподачи, обеспечивающих эффективное распределение воды и минимизации потерь в системе, путём разработки и уточнении математических моделей и алгоритмов, позволяющих обеспечить грамотное распределение влаги в почве. Это необходимо для водосбережения и рационального распределения водных ресурсов, чем поспособствует повышению надёжности подобранного оборудования и эффективности применения данного способа полива, при выращивании выбранной культуры в конкретных климатических условиях.
Методы и объект исследования. Успех применения данной технологии объясняется проведением грамотных гидравлических расчётов систем капельного орошения включает в себя определение различных зависимостей и параметров, таких как давление, расход, потери напора и другие [1]. Объектом исследования в данной теме рассматриваются основные зависимости и формулы, которые должны быть учтены при гидравлическом расчёте:
1. Закон сохранения объемов. Учитывается для системы капельного орошения при расчёте, что объем жидкости сохраняет постоянную величину только в том случае, если в ней нет пустот, и она движется как однородная:
Q вход = Q выход                                               (1)
где (Q) — расход воды.
    2. Потери напора. Могут быть вызваны трением, изменением направления потока и другими факторами. Основные формулы для расчета потерь напора:
Сопротивление трению. Для расчёта сопротивления на трение в трубопроводах можно использовать формулу Дарси-Вейсбаха:
Δh_т = f L/D⋅V^2/2g=f L/D ∫_O^V▒(V^2 d)/2g                                      (2)
где:
∆hт — потери напора (м),
 f — коэффициент трения,
 L — длина трубопровода (м),
 D — диаметр трубопровода (м),
 V — скорость потока (м/с),
 g — ускорение свободного падения (9.81 м/с²).
    2.2. Потери на местные сопротивления, можно рассчитать с помощью, но при использовании эмиттерных капельных лент, этими потерями можно пренебречь:
                                                       (3)
где:
∆h– потери напора на гидравлическом сопротивлении (м);
ξ– коэффициент местного сопротивления (коэффициент потерь);
V – средняя скорость течения жидкости (м/с);
g– ускорение свободного падения (9.81 м/с²). 
Из этого получим полное сопротивление, которое имеет вид:
     h=∆h_т+∆h_м+∆h_0                                                     (4)
где:
∆hт - потери напора на сопротивление трению (м);
∆hм – местные потри (м);
∆hо – общие потери напора (м).

    3. Формула Вейсбаха, определяющая потери давления на гидравлических сопротивлениях, имеет вид:
∆P=∆hρ≥d∆P=ρdh                            (5)
где:
ΔP – потери давления на гидравлическом сопротивлении (Па);
Δh – потери напора на сопротивление (м);
ρ – плотность жидкости (кг/м3).
    4. Давление в системе определяется с помощью расчета гидростатического давления, которое представляет собой сумму атмосферного давления, оказываемого атмосферой на поверхность жидкости, и давления, оказываемого самой жидкостью на глубине (высоте).

dP_n=P_атм+∫▒ρgdh                                (6)
где:
Pатм — атмосферное давление (Па);
ρ — плотность жидкости (кг/м3);
g — ускорение свободного падения (м/с2);
h — глубина, на которой находится жидкость (м).
    5. Расход через капельницы
Расход через капельницы можно рассчитать с использованием уравнения:
Q=k/((P -P_min )^n)                                                    (7)
где:
Q — расход через капельницу (м3/с),
k — коэффициент, зависящий от конструкции капельницы,
P — давление на входе капельницы (кПа),
Pmin - минимальное давление для работы капельницы (кПа),
n — экспонента, зависящая от типа капельницы.
    6. Для расчета оптимизации системы капельного полива необходимо учитывать при:
- распределение капельниц;
- длину и диаметр трубопроводов;
- условия эксплуатации (температура, тип почвы и т.д.).
Данные зависимости и формулы помогут при выполнении гидравлического расчета систем капельного орошения, для конкретных условий применения. 
Поливная сеть капельного орошения (участковые и поливные трубопроводы) характеризуется установившимся движением жидкости с переменным расходом по длине потока. Для расчетов поливной сети капельного орошения, где движение жидкости происходит с переменным расходом по длине потока, важно учитывать уравнения, описывающие динамику жидкости в трубопроводах. Основное уравнение, которое может быть использовано для анализа, основано на уравнении неразрывности и уравнении Бернулли:
6.1 Уравнение неразрывности (постоянства расхода) позволяют определить скорость потока в сечении трубы переменного сечения (если известны скорости в двух сечениях и площадь одного из них, можно найти скорость в другом сечении). Объясняют обратную зависимость между скоростью и площадью. Например, при движении жидкости по напорной трубе скорость в выходном сечении увеличивается обратно пропорционально площади живого сечения потока, будет иметь вид:
dp/dt+divpV=dp/dt+divpV+Vgradp=0                          (8)
6.2 Уравнение Бернулли
(pV^2)/2+pgh+P=const                                            (9)
 где: P — давление (Па); p — плотность жидкости (кг/м3); v — скорость потока (с); g — ускорение свободного падения (м/с2); h — высота (м).
При интегрировании уравнения, получим уравнение движения для учета переменного расхода и динамики жидкости в трубопроводах, которое выглядит так [2]:
h=(Q_тр^2+Q_п Q_тр+(Q_П^2)/3)  L_x/K^2 +(a_0 (a_2-2))/(gω^2 )(Q_тр Q_п+(Q_П^2)/2),                    (10)
где: Qтр — транзитный расход данного участка трубопровода (м3/с); Qп - путевой расход, разбираемый на данном участке трубопровода (м3/с); Lх — длина трубопровода от начального до рассматриваемого сечения (м); К - расходная характеристика (модуль расхода) участкового трубопровода/ коэффициент сопротивления (м3/с); а0 — коэффициент Буссинеска; а2 — коэффициент отсоединения массы жидкости; h - потери напора в рассматриваемом трубопроводе (м); g — ускорение свободного падения (м/с2); ω – угловая скорость (с).
Следует отметить, что результатом интегрирования дифференциального уравнения движения жидкости с переменной массой большинство исследователей делают допущения, что модуль расхода К и коэффициент отсоединения а2 являются постоянными величинами для всего трубопровода. Из этих условий и получено уравнение (1). Однако, как показали раннее проводимые исследования, при движении жидкости с переменным расходом изменяется по длине потока 4 коэффициент гидравлического трения, а следовательно, изменится коэффициент Шези и модуль расхода. 
Для определения численного значения величины коэффициента гидравлического трения λ в перфорированных полиэтиленовых трубопроводах были поставлены специальные опыты. Испытания проводились на трубопроводах разных диаметров с различной степенью перфорации. Изучались трубопроводы полиэтилена низкой плотности диаметром 16,2 и 25,0 мм, (ГОСТ 18599–2001) [3], длиной 2,4 и 20,0 м с калиброванными отверстиями 1,2 и 3,0 мм. Расстояние между отверстиями составляло 0,0625; 0,125; 0,25; 0,5; 1,0; 2,0; 4,0 м, т. е. опыты, проводились с различной интенсивностью отбора жидкости.
Давление жидкости в начальном сечении изменялось в диапазоне работы систем капельного орошения от 0,06 до 0,30 МПа.
Число Рейнольдса при этом варьировало от 1000 до 100 000.
Для сравнения со значением λ при равномерном движении жидкости на этих трубах проведены опыты при Q = const по длине трубопровода  рис. 4-3 стр.34 [2]. При обработке опытных данных получена зависимость
λ=(0,16Re^0.8+115)/R_e ,                                                   (11)
где: R_e=V_d/v — переменная по длине скорость в сечениях трубопровода; d — внутренний диаметр трубопровода; υ — коэффициент кинематической вязкости жидкости.
    Полученная зависимость достаточно хорошо описывает результаты эксперимента.
    Обозначив выражение 
А = 0,16 〖Re〗^0.8+ 115                                      (12)
а значение λ запишем в виде:
 λ=A/R_e                                                      (13)
Коэффициент Шези 
С = √((8g∙R_e)/A)                                                 (14)
И модуля расхода
 K=〖πD〗^2/4∙C√(R_гидр )                                             (15)
где: К – модуль расхода; R_гидр=ω/χ – гидравлический радиус для круглого поперечного сечения трубы рис. 3-3 стр. 21[2].
    Величина коэффициента массы жидкости а2 представляет собой отношение проекции на ось трубы скорости вытекающей струи к скорости основного потока. Для проверки численного значения коэффициента а2 при различном отклонении вытекающей струи нами проведены специальные исследования. В результате математической обработки данных эксперимента получены зависимости для определения коэффициента а2 при ламинарном и переходном режиме 300≤ Re ≤ 50 000                                                                                                          
a2=0,341 -10-4 Re—0,041- 10-8 Re2,                                (16)                                                                                     
    Число Рейнольдса в сечениях трубопровода, где установлены эмиттеры можно записать в виде:
Re=1.27 ∙  〖gn〗_i/Vd                                                    (18)
    Тогда теоретическую зависимость (10) можно преобразовать, если выразить транзитный и путевой расходы формулами:
Q_тр=n_i ∫_x ∫_y ∫_z d_x d_y d_z                                            (19)
Q_п=n_x ∫_x ∫_y ∫_z d_x d_y d_z                                             (20)                                                                       
Расход одного поливного трубопровода (ПТ) или одной капельницы участкового трубопровода (УТ) рассчитывается в зависимости от: n_x— количество отверстий на участке длиной ∆L; nі — количество отверстий от рассматриваемого сечения до конца трубопровода.
    После подстановки величин, преобразований и введений в уравнение (10) коэффициента дискретности, расчетная зависимость будет иметь вид:
Q_тр≥∑▒Q_п ⟹
h_dл=K_g^` q^2 (n_i^2+n_i n_x+(n_x^2)/3)  ∆L/K^2 +(1.67(a_2-2))/〖ga〗^2 ×q^2 (n_i n_x+(n_x^2)/2)                  (21)
где: ∆hдл – потери напора по длине участка перфорированного трубопровода (м); ∆L – длина участка рассматриваемого трубопровода (м); К – модуль расхода (м3/с). 
K_g^`=1+1,7 (L_отв/∆L)^1,04                                              (22)
где: Lотв – расстояние между отверстиями (м); К’g – коэффициент дискретности расхода, учитывающий влияние степени перфорации трубопровода.
    Правильность расчетов можно подтвердить с помощью подобных расчетов при условии: d=16,2 мм (dвн=15,9 мм); L=100 м; Lотв=0,5 м; n= 199 шт; Pатм=0,1 Мпа; g = 9,81 тогда
Δh_т = f L/D⋅V^2/2g=f L/D ∫_O^V▒(V^2 d)/2g
Средняя скорость течения воды в трубопроводе 
V=(4Q_п)/(πd_тр^2 ) =(4∙8,75∙10^(-5))/(79,4∙10^(-5) )=0,44 м/с
Q=g∙n_x=4,4∙10^(-7)∙199=8,75∙10^(-5) м^3/с
〖Re〗_эм=1.27 ∙  〖gn〗_i/(Vd_вн )=1,27∙(4,4∙10^(-7)∙199)/(0,44∙0,0159)=0,015875≈0,0436
Число Рейнольдса для потока воды в трубе
 〖Re〗_тр=(V∙d_вн)/ν=(0,44∙0,0159)/(1∙10^(-6) )=6996
    Так как Число Reтр = 6996, это подтверждает переходный режим течения, который находится в диапазоне от 2000 до 10000.    
    Теперь можно вычислить потери ∆hтр и потери hэ
∆h_тр=λ L/d_вн ∙V^2/2g=64/Re∙100/0,0159∙0,1936/(2∙9,81)=0,56 м
h_э=λ_э  L/d_вн ∙V^2/2g=0,0436∙100/0,0159∙〖0,44〗^2/(2∙9,81)=13,08 м
λ_э=(0,16∙〖Re〗^0,8+115)/Re=(0,16∙6996^0,8+115)/6996=0,0436
∆P=∆hρ=0,56∙1000=560∙9,81=5493,6Па≈5,5кПа
∑▒P=P_атм+∆P=101,3кПа+5,5кПа=106,8 кПа
    При сравнение полученных экспериментальных результатов с расчётными доказано совпадение, а также разница суммарных потерь напора, не превышающая 10% и разница абсолютной величины напора в 5%. Чем подтверждается зависимость (21), которая является расчетной теоретической и может быть применена для вычисления участковых трубопроводов систем капельного орошения.
    При этом учитывая, что величины модуля расхода К и коэффициента отсоединения а2 переменны по длине трубопровода, необходимо весь трубопровод разделить на отдельные участки, определить для каждого участка величину модуля расхода по формуле (15) и коэффициент отсоединения по формулам (16) или (17), а затем вычислить ∆hдл [6]. Примеры расчёта для участкового и поливного трубопроводов приведены в таблице 1.  

Таблица 1 – Сравнение расчета потерь напора в зависимости (21) с данными натурных исследований

Участковый трубопровод ∆l = 12м q2=2,596∙10-8 м3/с    Поливной трубопровод ∆l = 20м q2=0,0121∙10-10 м3/с
dh, мм    hрасч, м    hнат, м    % откл    dh, мм    hрасч, м    hнат, м    % откл
63    0,69    0,70    1,4    20    1,19    1,15    3,4
63    1,22    1,25    2,4    20    2,08    2,00    4,0
63    1,60    1,60    0,0    20    2,80    2,70    3,7
63    1,86    1,90    2,1    20    3,32    3,20    3,8
50    2,39    2,30    3,9    20    3,62    3,50    3,4
50    2,70    2,60    3,8    20    2,81    3,80    0,3
50    2,83    2,70    4,8    20    3,89    3,90    0,2
50    2,87    2,80    2,5    20    3,91    3,95    1,0
                    3,94    4,00    1,5
    
    Местные потери напора для участковых трубопроводов можно учитывать путём прибавления 5% к потерям по длине. Зависимость (21) также применяется для расчёта поливных трубопроводов систем капельного орошения. При этом следует иметь в виду, что q – расход одной капельницы; nх – количество капельницы на длине ∆L поливного трубопровода; nі – количество поливных трубопроводов от рассматриваемого сечения до конца трубопровода. В этом случае принимают а2=0, так как в результате направляющего действия штуцера капельницы струя отсоединяется под углом 90°. Для оценки надёжности предлагаемых зависимостей выполнены натурные измерения потерь напора в участковых и поливных трубопроводах систем капельного орошения. Сравнение экспериментальных результатов с расчётными показало хорошее совпадение. Разница суммарных потерь напора не превышала 10%, а разница абсолютной величины напора 5%. 
    Выводы и рекомендации. При исследовании зависимости для гидравлического расчета систем капельного орошения, можно сделать заключение. Для определения оптимальных параметров, обеспечивающих эффективное распределение воды и минимизацию потерь в системе, необходимо учитывать следующие аспекты:
    - Определение потерь давления в трубопроводах, соединениях и капельницах, для обеспечения необходимого (постоянного) давления при эксплуатации.
    - Определение необходимого расхода воды для выращивания культур в благоприятных климатических условиях.
    - Грамотный подбор оборудования для различных типов капельниц и насосов, при учете их производительности и комплектации с системой.
    - Оптимизация проектирования при разработке схемы расположения трубопроводов и капельниц, с целью минимизировать затраты на материалы и монтаж, с соблюдением водораспределения.
    - Необходимо системное анализирование влияния внешних факторов окружающей среды, таких как: температура, влажность, тип почвы, эффективность влияния подобранной системы орошения.
    - Использование программного обеспечения для моделирования гидравлических процессов и тестирования различных сценариев работы системы.
    - Оценка затрат при монтаже установки и эксплуатации системы капельного орошения в сравнении с другими методами орошения.
    Таким образом, исследование зависимости в гидравлическом расчете систем капельного орошения позволяет разработать эффективные и экономически целесообразные решения для сельского хозяйства и ландшафтного дизайна.
 

References

1. Grushin A.V., Grzhibovsky S.A., Bulgakov D.V., Kolomeets A.V. Aspects of development and features of drip irrigation // Bulletin of Meliorative Science, 2021, No. 3. pp. 57-66.

2. Handbook of hydraulic calculations. Edited by P.G. Kiselyov. 4th Ed., revised. and additional M., "Energy", 1972.

3. GOST 18599-2001 Polyethylene pipes. Technical specifications. : approved and put into effect by the Interstate Council for Standardization, Metrology and Certification (Protocol No. 20 of November 1, 2001) : By Resolution No. 112-st of the State Committee for Standardization and Metrology of the Russian Federation dated March 23, 2002, the interstate standard GOST 18599-2001 was put into effect directly as the state standard of the Russian Federation from January 1 2003 : instead of GOST 18599-83 : Minsk : IPK Publishing House of Standards, 2003. – 38 p. – Text: direct.

4. Safronova T.I., Fiyatkin D.A. Hydraulic calculation of drip irrigation systems // Scientific and educational journal for students and teachers "StudNent", 2020, No. 2, pp. 143-147.

5. Semenenko S.Ya., Lytov M.N., Chushkin A.N., Chushkina E.I. Algorithm for calculating drip irrigation systems with an electrochemical activation module for irrigation water //Scientific Journal of RosNIIPM, 2017, No.4(28), pp. 20-36.

6. Mikityuk A.V., Kazharov V.M., Shugai P.Yu. Hydraulic calculation of an irrigation polyethylene pipeline of a drip irrigation system // Scientific Electronic Journal of KubGAU, 2005, No. 05(13). pp. 105-111.

Login or Create
* Forgot password?