Совершенствование способов определения расхода воды в открытых мелиоративных каналах
В.Н. Щедрин, М.В. Вайнберг, А.А. Чураев Российский научно-исследовательский институт проблем мелиорации
Аннотация
Целью исследований являлся анализ точности определения коэффициента шероховатости. Для этого были использованы два способа определения расхода воды по методу «уклон - площадь»: с учетом перепада свободной поверхности воды, без учета перепада. За эталонный расход воды при этом был принят расход, рассчитанный объемным методом. В процессе исследований были проведены опыты по определению влияния глубины потока и расхода воды на коэффициент шероховатости при изменении уклона дна лотка от 0,0002 до 0,0006. Для подтверждения достоверности полученных данных произведена статистическая обработка значений коэффициентов шероховатости, вычисленных без учета и с учетом перепада свободной поверхности воды. В первом опыте коэффициент вариации составил 0,82 % против 1,27 %, во втором - 0,81 % против 1,23 % и в третьем - 0,59 % против 1,64 %. Полученные данные были проанализированы с целью определения влияния перепада свободной поверхности на точность определения расхода воды. В результате проведенных расчетов относительная погрешность расхода воды без учета перепада свободной поверхности составила в первом опыте дQmax=2,51%, во втором дQmax=2,83% и в третьем дQmax=3,24%. С учетом перепада свободной поверхности относительная погрешность по расходу составила в первом опыте дQmax=1,67%, во втором опыте дQmax=1,84%, в третьем дQmax=1,25 %. На основании полученных результатов можно сделать вывод, что определение коэффициента шероховатости и расхода воды по предложенному способу необходимо осуществлять с учетом перепада свободной поверхности воды, при этом предложенный способ позволяет устанавливать и отслеживать величину расхода воды в непрерывном режиме (онлайн).
Ключевые слова: коэффициент шероховатости, метод «уклон - площадь», расход воды, перепад свободной поверхности воды, коэффициент Шези.
Abstract
The objective of the research was to analyze the accuracy of defining roughness coefficient. For this purpose two methods were used to determine water discharge by “slope-area” method: considering drop of free surface of water and excluding the drop. As a reference water discharge the discharge calculating by volume method was assumed. During the investigation experiments on defining the impact of flow depth and water discharge on roughness coefficient at the range of tray bottom slope 0.0002 to 0.0006 were carried out. To confirm reliability of data obtained, statistical treatment of values of roughness coefficient, calculated considering drop of free surface of water and excluding the drop, was done. In the first experiment variation coefficient was 0.82 % against 1.27 %, in the second one - 0.81% against 1.23 %, and in the third - 0.59 % against 1.64 %. Data obtained were analyzed in order to define the impact of free surface drop on the accuracy of determining water discharge. As a result of conducted calculations relative error in the first experiment was дQmax = 2.51 %, in the second - дQmax = 2.83 %, and in the third - дQmax = 3.24 %. Considering the drop of free surface relative error on discharge in the first experiment was дQmax = 1.67 %, in the second - дQmax = 1.84 %, and in the third - дQmax = 1.25 %. Based on the results obtained, the conclusion can be formulated that defining of roughness coefficient and water discharge according to proposed method should be done considering drop of free surface of water. While the method proposed allows determining and monitoring the value of water discharge in a continuous mode (on-line).
Keywords: roughness coefficient, slope-area method, water discharge, drop of free surface of water, Chezy coefficient.
Выбор значения коэффициента шероховатости равносилен определению сопротивления течению, природа которого на современном этапе знаний, в сущности, непостижима [1-6]. Опытные инженеры принимают этот коэффициент, основываясь на своей многолетней практике, а начинающие - наугад, причем результаты значительно различаются.
В каналах оросительных систем значение коэффициента шероховатости весьма изменчиво и зависит от большого количества факторов [7], таких как наличие растительности, неоднородность грунта ложа канала, заиление и размыв, размеры и форма канала и т. д. Оно не может приниматься постоянным и единственным для всех случаев.
От точности определения коэффициента шероховатости зависит точность вычисления расхода воды, так как является элементом, входящим в состав расчетных выражений для определения расхода.
Известен способ определения перепада уровней и расхода воды на открытом водотоке с призматическим руслом по методу «уклон - площадь», разработанный А.Е. Ивахненко (RU 2327114 C 2, G 01 F 1/00, от 20.06.2008) [8]. Сущность способа сводится к использованию одного уровнемера, расположенного в уровнемерном колодце нижнего гидрометрического створа, для измерения уровней в верхнем и нижнем гидрометрических створах, определению перепада уровней между верхним и нижним створами и вычислению расхода. Достоинствами способа являются:
- повышение точности и достоверности определения расхода воды за счет более точного измерения перепада уровней воды в верхнем и нижнем измерительных створах, так как измерения производятся одним измерительным прибором (уровнемером);
- возможность автоматизации процесса измерений;
- возможность уменьшения величины обратного уклона или полный отказ от его устройства;
- снижение затрат на эксплуатацию измерительного участка ввиду близости гидрометрических створов и отсутствия необходимости в строительстве второго уровнемерного колодца.
Недостатками приведенного способа являются:
- отсутствие возможности вычисления расхода воды в режиме непрерывного времени, так как уровни воды в створах и перепад определяются с интервалом 20-30 мин [8];
- низкая надежность измерительного комплекса ввиду наличия движущихся механических элементов (поплавка, задвижки);
- необходимость устройства подводящей трубы с запорным клапаном к успокоительному колодцу с измерительным оборудованием.
С целью устранения недостатков способа измерения, предложенного А.Е. Ивахненко для реализации метода измерения «уклон - площадь», нами предлагается новый способ определения расхода воды на открытых каналах оросительных систем [9]. Сущность его сводится к использованию двух датчиков уровня воды, оснащенных средствами дистанционной передачи показаний уровня и расположенных в уровнемерных колодцах верхнего и нижнего гидрометрических створов, определению уровней воды в створах, перепада уровней между верхним и нижним створами и вычислению расхода воды.
Использование предложенного нами способа измерения расхода воды по методу «уклон - площадь» позволит эксплуатационным организациям, обслуживающим современные оросительные системы, организовать автоматизированный водоучет, который будет характеризоваться следующими достоинствами:
- отслеживание в непрерывном режиме значений уровней воды в створах, перепада уровней между створами;
- полная автоматизация процесса измерений;
- повышение надежности за счет отсутствия движущихся механических частей измерительного оборудования;
- отсутствие необходимости в устройстве подводящей трубы с запорным клапаном.
При использовании предлагаемого способа, измеряя необходимые параметры водотока с помощью соответствующего приборного и программного обеспечения, на выходе получаем значения расхода в режиме онлайн. Однако погрешность использования предложенного нами способа и оборудования для реализации метода «уклон - площадь» необходимо было подтвердить результатами лабораторных исследований.
Материалы и методы. Метод «уклон - площадь» основан на измерении уклона свободной поверхности воды и площади живого сечения с учетом гидравлического сопротивления русла. Предложенный нами способ реализуется следующим образом. Из измерительного канала по соединительным трубам вода поступает в успокоительные колодцы в верхнем и нижнем гидрометрических створах. Когда течение воды установится, в успокоительных колодцах датчиками уровня воды будут непрерывно регистрироваться измеряемые параметры с заданным интервалом, и с помощью средств дистанционной передачи информация будет передаваться на пункт диспетчера, оснащенный средствами ее обработки и вычисления расхода. По полученным данным и при известных параметрах измерительного участка канала (расстояние между сечениями верхнего и нижнего гидрометрического створа , уклон дна канала , коэффициент шероховатости канала , ширина канала по дну , коэффициент заложения откосов и т. д.) вычисляется искомый расход.
Исследования по определению расхода воды методом «уклон - площадь» проводились в лабораторных условиях. В качестве средства измерения расхода воды применялся прямолинейный участок трубы с полукруглым поперечным сечением, являющийся моделью лотка, оснащенный средством измерения уровней воды в двух створах и удовлетворяющий требованиям точности измерения уклона поверхности воды. С помощью насосного агрегата вода поступала в емкость для подачи воды в лоток лабораторной установки. Проводились измерения глубины потока в верхнем и нижнем гидрометрических створах одновременно (рисунок 1). По измеренным данным вычислялись , и , .
Рисунок 1 - Лабораторная установка: 1 - емкость для используемой воды; 2 - задвижка; 3 - насос; 4 - емкость для подачи воды; 5 - затвор; 6 - датчики уровня воды для определения и в верхнем и нижнем створах; 7 - труба; 8 - жалюзи; 9 - емкость для сбросной воды; 10 - контроллер для приема и обработки дистанционного сигнала с датчиков уровня воды; 11 - компьютер с программой вычисления расхода воды
Расход воды определялся по формуле, являющейся решением дифференциального уравнения неравномерного установившегося плавно изменяющегося движения жидкости в открытом канале оросительной сети [10, 11]. Для прямого ( > 0) уклона формула имеет вид:
,(1)
где - расход, м3/с;
- площадь живого сечения в створах 1 и 2, м2;
, - глубина воды в створах 1 и 2, м;
- уклон дна канала;
- расстояние между створами, м;
- коэффициент скорости, равный 1,1;
- ширина канала по свободной поверхности воды в створах 1 и 2, м;
- ускорение свободного падения, равное 9,81 м/с2;
- коэффициент Шези в створах 1 и 2, м0,5/с;
- гидравлические радиусы в створах 1 и 2, м.
Для определения фактического коэффициента шероховатости использовалось основное уравнение расхода воды при равномерном движении жидкости [12]:
.(2)
Используя выражение (2), определяли коэффициент Шези :
.
В практических расчетах открытых русел преимущественно призматической формы наиболее подходящей для определения коэффициента шероховатости является формула Агроскина, которая при проведении лабораторных исследований позволила получить значения п, наиболее близкие к определенным экспериментально:
, отсюда .
Расчеты производились с использованием программы Microsoft Excel.
Результаты и обсуждение. Был проведен анализ точности определения коэффициента шероховатости путем решения обратной задачи двумя выражениями: предлагаемым выражением (1), в состав которого входит перепад свободной поверхности воды, и формулой Шези (2), в состав которой не входит перепад свободной поверхности воды.
Эталонный расход воды определялся объемным методом. Изменение расхода воды осуществлялось с помощью затвора 5 (рисунок 1).
Таблица 1 - Значения измерения расхода воды, подаваемого насосами
|
Объем воды в мерной емкости W, л |
Продолжительность заполнения мерной емкости t, с |
Расход воды Q, л/с |
|
|
При заданном расходе Q = 0,3 л/с |
|||
|
18,6 |
60,5 |
0,307438 |
|
|
18,5 |
60,4 |
0,306291 |
|
|
18,3 |
60,1 |
0,304493 |
|
|
18,1 |
59,9 |
0,302170 |
|
|
18,2 |
60,0 |
0,303333 |
|
|
При заданном расходе Q = 0,5 л/с |
|||
|
30,4 |
60,2 |
0,504983 |
|
|
30,5 |
60,4 |
0,504967 |
|
|
30,3 |
60,2 |
0,503322 |
|
|
29,9 |
60,0 |
0,498333 |
|
|
30,1 |
60,1 |
0,500832 |
|
|
При заданном расходе Q = 0,7 л/с |
|||
|
42,0 |
60,2 |
0,697674 |
|
|
42,7 |
60,6 |
0,704620 |
|
|
42,5 |
60,5 |
0,702479 |
|
|
41,8 |
60,1 |
0,695507 |
|
|
42,3 |
60,3 |
0,701493 |
|
|
При заданном расходе Q = 0,9 л/с |
|||
|
54,0 |
60,4 |
0,894040 |
|
|
53,8 |
60,3 |
0,892206 |
|
|
54,1 |
60,4 |
0,895695 |
|
|
54,5 |
60,5 |
0,900826 |
|
|
54,6 |
60,5 |
0,902479 |
|
|
При заданном расходе Q = 1,1 л/с |
|||
|
66,5 |
60,3 |
1,102819 |
|
|
66,2 |
60,1 |
1,101498 |
|
|
66,4 |
60,2 |
1,102990 |
|
|
66,3 |
60,1 |
1,103161 |
|
|
66,6 |
60,4 |
1,102649 |
Расход, подаваемый насосом, вычисляли по формуле:
Q=W/t,
где W - объем воды в мерной емкости, л;
t - продолжительность заполнения мерной емкости, с.
Значения рассчитанного объемным методом расхода воды приведены в таблице 1. Измерения расхода проводились не менее пяти раз.
При определении расхода воды в лабораторных условиях могут возникать случайные ошибки, которые необходимо учитывать с помощью методов теории вероятности [13], затем следует вносить соответствующие поправки.
В классической теории ошибок доказано, что случайные ошибки подчиняются нормальному закону распределения. Для оценки колеблемости значений признака относительно средней ошибки используются характеристики рассеяния. Они различаются выбранной формой и способами оценки отклонений от нее отдельных вариантов [14]. К таким показателям относятся среднее квадратическое отклонение, дисперсия и коэффициент вариации.
Оценка отклонений у определяется по формуле [13, 15]:
,
где - средняя арифметическая простая, которая определяется по формуле:
,
где - число измерений;
- значение признака в дискретном ряду;
Дисперсией измерений является квадрат величины , определяемый по формуле [16]:
.
Для оценки меры вариации и ее значимости используют коэффициент вариации V, который дает относительную оценку вариации [17]. Коэффициент определяется по следующей формуле и выражается в процентах:
.
Коэффициенты вариации дают относительную характеристику однородности явлений и процессов. Они позволяют сравнивать степень вариации разных признаков. Статистическая обработка значений задаваемого эталонного расхода приведена в таблице 2.
Таблица 2 - Статистическая обработка измерения задаваемого эталонного расхода
|
Значение расхода Qi, л/с |
Дисперсия D |
Коэффициент вариации V, % |
|
|
0,307438 |
4,58256 · 10-6 |
0,702453 |
|
|
0,306291 |
|||
|
0,304493 |
|||
|
0,302170 |
|||
|
0,303333 |
|||
|
Сумма 1,523725 |
|||
|
Среднее значение 0,304745 |
|||
|
0,504983 |
8,2681 · 10-6 |
0,572239 |
|
|
0,504967 |
|||
|
0,503322 |
|||
|
0,498333 |
|||
|
0,500832 |
|||
|
Сумма 2,512437 |
|||
|
Среднее значение 0,502487 |
|||
|
0,697674 |
1,36719 · 10-5 |
0,527954 |
|
|
0,704620 |
|||
|
0,702479 |
|||
|
0,695507 |
|||
|
0,701493 |
|||
|
Сумма 3,501773 |
|||
|
Среднее значение 0,700355 |
|||
|
0,894040 |
1,95232 · 10-5 |
0,49256 |
|
|
0,892206 |
|||
|
0,895695 |
|||
|
0,900826 |
|||
|
0,902479 |
|||
|
Сумма 4,485246 |
|||
|
Среднее значение 0,897049 |
|||
|
1,102819 |
4,32212 · 10-7 |
0,059624 |
|
|
1,101498 |
|||
|
1,102990 |
|||
|
1,103161 |
|||
|
1,102649 |
|||
|
Сумма 5,513117 |
|||
|
Среднее значение 1,102623 |