5.3.4 Для вычисления среднего суточного уровня за предшествующие сутки необходимо проанализировать полученные данные срочных наблюдений за уровнями воды.
5.5.4.1 Если в полученной таблице А.2 (приложение А) за прошедшие сутки нет пропусков наблюдений и наблюдения выполнены через равные промежутки времени в течение всех суток, то средний суточный уровень воды за прошедшие сутки определяется как среднее арифметическое из всех наблюденных значений уровня воды.
5.5.4.2 Если для рассматриваемых суток специалист-гидролог принял решение об игнорировании неполных данных АГК и использовании занесенных вручную данных, то для расчета среднего суточного уровня воды используются только данные ручных измерений. Средний суточный уровень воды вычисляется как среднее арифметическое из значений, наблюденных в 8 и 20 ч.
При отсутствии реального значения уровня воды хотя бы в один из двух стандартных сроков средний суточный уровень воды не определяется, формируется признак пропуска наблюдений.
5.5.4.3 Если интервал между измерениями уровня воды с помощью АГК или ручным способом изменялся в течение прошедших суток, то определение среднего суточного уровня воды производится по всем имеющимся в файле данных АГК измерениям уровня воды следующим образом.
Средний суточный уровень воды Hср, см определяется по формуле среднего взвешенного значения. При этом, если имеется уровень в 0 ч текущих суток, то он также участвует в расчете и принимается за 24 ч предшествующих суток.
Hср = ((H1 + H2)(T2 - T1)/2 + (H2 + H3)(T3 - T2)/2 +...
+ (Hn-1 + Hn)(Tn - Tn-1)/2)/(Tn - T1), (7)
где Hср - средний суточный уровень воды, см;
T1, T2, ..., Tn - сроки наблюдений, следующие в хронологически возрастающем порядке, в сутках;
H1, H2, ..., Hn - значения уровня воды в сантиметрах в соответствующие сроки T1, T2, ..., Tn.
Вычисление производится по неокругленным данным, затем средний суточный уровень воды округляется до целых сантиметров.
Таким же образом рассчитывается средний суточный уровень воды по данным ручных наблюдений в стандартные и нестандартные сроки. При этом, значение уровня воды на 1 ч предшествующих суток и 0 ч текущих (или 24 ч предшествующих) суток вычисляются методом линейной интерполяции между смежными наблюденными значениями, в интервал между которыми попадает срок 1 или 24 ч предшествующих суток.
5.5.5 Высший Hвысш и низший Hнизш уровни воды за сутки определяются из тех же значений уровня воды, прошедших первичный контроль, что и средний суточный уровень. Выбранные значения округляются до целых сантиметров.
5.5.6 По занесенным кодам СВО в рассматриваемые сутки проверяется наличие затора, зажора, сгона, нагона, а также подпора от плотины, озера (реки), засорения русла, мостовых переправ. Если хоть одно из этих явлений имело место, то выбранным экстремальным за сутки уровням присваивается признак подпора, отметка о котором делается в графе "Примечание" таблицы А1 (приложение А). Этот же признак присваивается среднему суточному уровню воды.
6.1.1 Задача оперативного учета стока заключается в получении срочных и ежедневных значений расходов воды (ЕРВ) по данным актуальных (оперативных) наблюдений за уровнем воды и совокупности ИРВ, выполненных в период, предшествующий моменту расчета. Отличие оперативного учета стока от режимного заключается в следующем. Режимные характеристики стока вычисляют по истечении года, когда специалист-гидролог имеет в распоряжении всю годовую совокупность измеренных расходов и уровней воды, дополненную сведениями о состоянии реки и ее русла в течение всего годового цикла. При этом специалист-гидролог располагает сведениями о состоянии пропускной способности русла (ПСР) реки на моменты измерений расходов воды и должен воспроизвести эти характеристики в интервалах между измерениями, т.е. выполнить их интерполяцию в ретроспективном варианте. Оперативный же учет стока ведется в условиях существенно меньшей определенности, связанной с необходимостью прогнозировать (экстраполировать) характеристики ПСР в промежутке времени от даты состоявшегося измерения до текущего момента времени. Поэтому важно, чтобы оба варианта учета стока опирались на единую объективную основу, учитывающую особенности изменения ПСР в конкретном гидрометрическом створе. Это позволит преодолеть разрыв между оперативными и режимными данными. До настоящего времени нередко оперативные значения расходов воды до 30% - 50%, а иногда и до 100% отличаются от данных, публикуемых в ежегодниках.
6.1.2 Как в режимном, так и в оперативном вариантах учета стока в качестве первого приближения к значению расхода воды при фиксированном уровне принимается значение расхода воды, полученное по опорной кривой расходов (КР) (зависимости расходов от уровней воды Q(Ht) для текущего уровня воды Ht. Для уточнения расхода, вычисленного по опорной КР, используются поправки к ней
6.1.3 В общем случае значение оперативного расхода воды на момент (срок) t Qоп, м3/с, вычисляется по формуле
В качестве характеристик изменения ПСР во времени используются относительные отклонения ИРВ
где Qи - ИРВ, м3/с;
Q(Hи) - расход воды, полученный по КР Q(H) для уровня Hи, зафиксированного во время измерения расхода воды Qи, м3/с.
Изменения ПСР происходят под воздействием ряда сезонных, гидроморфологических и динамических факторов: таких как, ледовые явления зимой и зарастание в летнюю межень, деформации неустойчивого русла, переменный подпор от впадающих ниже притоков или размещенных в русле гидротехнических сооружений и т.п. Влияние этих факторов обусловливает разброс точек ИРВ относительно аппроксимирующей зависимости Q(H) и составляет физическую причину отклонений
Наиболее обоснованным и простым при расчетах оказывается оперативный учет стока в условиях свободного состояния русла при наличии многолетней кривой расходов (МКР) Q(H), освещенной измерениями в широком диапазоне изменения уровней воды, включая их минимальные и максимальные многолетние значения, и подтвержденной измерениями расходов воды расчетного года. В этом случае оперативные значения ЕРВ определяются непосредственно по многолетней кривой расходов (МКР) или соответствующему ей аналитическому уравнению регрессии (УРМ), аппроксимирующему многолетнюю кривую расходов воды по уровням. Поправка
Нарушение однозначности КР происходит также в условиях сезонного изменения ПСР при наличии ледовых явлений в зимний и переходные периоды и зарастания в теплое время года.
6.2.1 В соответствии с основными положениями методики, изложенными в 6.1, прежде всего для исследуемого гидрометрического створа необходимо выбрать опорную КР, которая будет использоваться для расчета оперативных расходов воды. В качестве опорной КР для оперативного учета стока рекомендуется принимать:
- МКР в случаях, когда в течение последних нескольких лет не происходило существенных изменений русла и режима реки;
- КР, полученную по данным ИРВ одного-двух последних лет, если такие изменения прослеживались.
Об изменениях, происходящих в русле реки, можно судить на основе анализа положения годовых КР на графике Q(H) и хронологических графиков изменения средних и минимальных отметок дна в исследуемом гидрометрическом створе в течение ряда лет.
6.2.2 Для удобства выполнения расчетов необходимо представить зависимость в виде аналитического выражения, то есть выполнить аппроксимацию зависимости расходов от уровней воды Q(H). Уравнение аналитической зависимости расходов от уровней воды может быть получено непосредственно по данным совокупности ИРВ за многолетие при использовании в качестве опорной МКР или по совокупности данных одного или двух последних лет. Аналитическое уравнение может быть получено также по данным (опорным точкам) таблицы координат выбранной для использования опорной КР.
Для аналитической аппроксимации зависимости расходов Q от уровней воды H рекомендуется использовать следующие уравнения [5]:
параболу Глушкова
или полиномиальную зависимость в виде
где Q - расход воды, м3/с;
H - уровень воды над "нулем" поста, м;
H0 - уровень, при котором расход воды в исследуемом гидрометрическом створе равен нулю, м;
a и b - параметры уравнения (10), определяемые при аппроксимации по данным совокупности ИРВ;
a0, ..., am - параметры уравнений регрессии (УР) расходов по уровням воды в форме полинома (11).
Аппроксимация зависимости расходов от уровней воды выполняется с использованием метода наименьших квадратов (МНК) [8]. Для удобства представления коэффициентов аналитических уравнений при выполнении расчетов по МНК значения уровней рекомендуется выражать не в сантиметрах над "нулем" поста, как это принято в гидрологии, а в метрах над "нулем" поста.
Подробный алгоритм аппроксимации изложен в приложении Б, примеры такой аппроксимации даны в приложениях Г и Д.
6.2.3 Для принятия решения об использовании того или иного метода учета стока необходимо выполнить оценку однозначности полученной опорной зависимости расходов от уровней воды. Такая оценка выполняется на основе анализа рядов относительных отклонений ИРВ
6.3.1 Как уже отмечалось выше, наиболее обоснованным и простым при расчетах оказывается оперативный учет стока при наличии МКР, освещенной измерениями во всем диапазоне изменения уровней воды, включая минимальные и максимальные исторические значения. МКР может считаться обоснованной, если она получена при использовании данных ИРВ не менее чем за три года наблюдений на ГП. Для ГП с длительными рядами наблюдений рекомендуется использовать сведения об ИРВ за последние 10 лет; данные предшествующего периода привлекаются для аппроксимации МКР только в случаях, когда они освещают исключительно высокие паводки и могут уточнить ее положение в верхней части.
6.3.2 МКР может считаться установленной, если:
- обоснована 50 - 60 измерениями расходов в течение не менее трех лет в периоды свободного состояния русла;
- ее экстраполяция не превышает вверх 15%, вниз 5% от полной амплитуды уровней за прошлые годы соответственно в верхнем и нижнем диапазонах уровней воды;
- в соответствии с приложением В установлена ее однозначность по методике оценки однозначности;
- точки ИРВ в пределах полосы рассеяния не образуют хронологически обособленных групп, обусловленных, например, прохождением паводочных волн или однонаправленных деформаций русла.
6.3.3 При установленном факте наличия однозначной МКР оперативные значения ЕРВ для данного ГП вычисляются непосредственно по МКР или соответствующему ей УРМ при подстановке в него соответствующего срочного оперативного значения уровня воды. Поправка
6.3.4 Для вычисления оперативных расходов воды рекомендуется использовать уравнение, полученное в ходе аппроксимации МКР, но допускается также использование координат этой зависимости, рассчитанных по всему диапазону изменения уровня воды в расчетном створе, представленных в табличной форме. Примеры аппроксимации МКР приведены в приложениях Г и Д.
6.3.5 Каждый вновь измеренный расход воды Qи текущего года используется для контроля и подтверждения положения многолетней зависимости Q(H). Вычисляется относительное значение его отклонения от принятой для расчета МКР -
Здесь Hи - уровень воды, наблюденный при измерении расхода воды.
МКР считается подтвержденной, если относительное отклонение
где
, (14)где
N - количество ИРВ в многолетней совокупности данных, принятых для получения МКР.
6.3.6 При анализе отклонений каждого вновь поступившего измерения от МКР может наблюдаться несколько случаев:
- соотношение (13) выполняется. Многолетняя зависимость Q(H) признается подтвержденной, значения оперативных расходов воды вычисляются непосредственно по МКР или аналитическому уравнению (УРМ);
- соотношение (13) не соблюдается, следовательно, МКР признается неподтвержденной. В этом случае требуется назначить дополнительное незапланированное измерение расхода воды в гидрометрическом створе и после поступления данных вновь выполнить проверку в соответствии с условиями (12), (13). Если по очередному ИРВ критерий (13) вновь не соблюдается, МКР признается неподтвержденной и не может быть принята в качестве опорной для оперативного учета стока. В этом случае при расчете оперативных расходов воды следует использовать метод, описанный в 6.5 для оперативного учета стока при отсутствии однозначной КР.
6.3.7 Каждый вновь измеренный расход воды для дополнительного визуального контроля наносится на график зависимости Q(H), представленный на рисунке 1.
![]() 6.4.1 Если для расчетного ГП устойчивая МКР отсутствует, для оперативного учета стока рекомендуется использовать КР за один или два предшествующих года.
6.4.2 Годовая зависимость расходов от уровней воды, как правило, если это не выдающийся многоводный год, строится в ограниченном диапазоне уровней воды и для ее использования в целях оперативного учета стока в расчетном году требуется ее экстраполяция вверх и/или вниз за пределы диапазона уровней, освещенных измерениями расходов воды.
6.4.3 Экстраполяцию как вверх, так и вниз за пределы диапазона измерения рекомендуется выполнять с использованием данных ИРВ прошлых лет и расчетных значений расходов, полученных по уравнению, предложенному И.Ф. Карасевым [5]
где Qэ - расход воды в любой фазе водного режима за пределами диапазона измерений, который необходимо рассчитать, м3/с;
Qi - расход воды, вычисленный по УР Q(H) в пределах диапазона измерений, м3/с;
Параметр mi для каждого интервала между измерениями расходов воды определяется по формуле
где Qi и Qi-1 - смежные во времени на фазе подъема в верхнем диапазоне уровней или фазе спада в нижнем диапазоне уровней значения расходов воды, но не измеренные, а вычисленные по УР для значений уровней воды, наблюдавшихся при измерении, м3/с;
6.4.4 Для расчета параметра mi, кроме кривой расходов, требуется использовать зависимость площадей поперечного сечения в исследуемом гидростворе от уровня воды, полученную по данным промеров в гидростворе и нивелировки поперечного профиля до отметок наивысших наблюденных исторических уровней воды. Аппроксимация этой зависимости выполняется по алгоритму, описанному в приложении Б.
6.4.5 Расчет по формуле (16) будет корректен при условии, что расходы Qi и Qi-1 не слишком близки между собой - модуль их разности должен существенно превышать погрешность измерения расходов воды, а его относительное значение составлять не менее 20% от среднего значения расходов Qi и Qi-1.
Параметр mi рекомендуется определять для нескольких пар расходов воды в пределах одной фазы водного режима и принимать осредненное в ее пределах значение. В диапазоне экстраполяции по формуле (15) рекомендуется вычислить не менее 7 - 10 расходов воды Qэi.
6.4.6 Процедура экстраполяции зависимости Q(H) выполняется с использованием метода кусочно-полиномиальной регрессии [5], [7] и состоит в следующем. По уравнению Q(H), полученному в пределах диапазона уровней, освещенных измерениями расходов воды, вычисляется значение Q в верхней точке этой КР Q(Hвгр). За пределами этого диапазона уровней вплоть до наивысшего наблюденного уровня воды формируется совокупность, состоящая из отобранных ИРВ прошлых лет в этом диапазоне и рассчитанных по формуле (15) значений. Затем на основе метода полиномиальной регрессии с ограничениями определяются параметры зависимости Q(H) в диапазоне экстраполяции.
Приведенный метод экстраполяции может быть использован как в верхней, так и в нижней частях диапазона изменения уровней и расходов воды.
6.5 Алгоритм оперативного учета стока в условиях свободного состояния русла при отсутствии однозначной кривой расходов
6.5.1 Основой для выбора метода гидрометрического учета стока служит оценка однозначности зависимости расходов от уровней воды, которая выполняется на основе сравнения дисперсии относительных отклонений ИРВ от осредняющей КР и дисперсии случайных погрешностей ИРВ [7].
6.5.2 При получении оперативных значений расходов воды следует исходить из того, что принятая опорная КР Q(H) осредняет поле точек ИРВ за предшествующий период и, следовательно, лишь приближенно отражает характеристики ПСР. Поэтому значение расхода воды, полученное по опорной КР Q(H), принимается в качестве детерминированной характеристики ПСР и является первым приближением к искомому значению расхода воды в оперативном варианте. Его можно существенно уточнить, если учесть изменения ПСР, нарушающие однозначность связи Q(H) в соответствии с основными положениями методики, изложенными в 6.1.
6.5.3 При оперативном учете стока необходимо получить значение поправки
6.5.4 При нарушении однозначности связи Q(H) временной ряд относительных отклонений
где
- относительное значение дисперсии случайной погрешности измерения расходов воды;6.5.5 При использовании метода оптимальной экстраполяции значение
где pi - весовые коэффициенты;
Схема индексных обозначений для временного ряда отклонений
![]() ряда отклонений
6.5.6 Для применения метода оптимальной экстраполяции требуется надежная оценка статистических характеристик процесса
![]() 1 - маловодных; 2 - средних по водности; 3 - многоводных
в створе ГП р. Оскол - с. Ниновка
6.5.7 При использовании метода оптимальной экстраполяции элементов установлено, что для рядов гидрологических данных в расчетах достаточно опираться на два последних измерения расхода воды.
Расчет весовых коэффициентов оптимальной экстраполяции в этом случае производится по формулам:
где
6.5.8 Оценка математического ожидания (среднего значения) элементов относительных отклонений
6.5.9 Иногда в расчетах приходится опираться только на одно последнее значение ИРВ. Это случается, например, при редких измерениях расходов воды, которые обусловливают малую корреляционную связанность отклонений от опорной КР. В этом случае необходимые характеристики вычисляются по формуле:
где
p1 - коэффициент оптимальной экстраполяции, который вычисляется по формуле
где
6.5.10 При невозможности надежной оценки АКФ, а, следовательно, и использования метода оптимальной экстраполяции, поправки
Вычисление поправки для случая, когда в расчет принимаются данные только двух последних измерений расходов воды, производится по формуле
где
6.5.11 Оперативные расходы воды вычисляются ежедневно после поступления сведений о срочных уровнях воды для каждого измеренного значения срочного уровня воды по формуле (8) и заносятся в таблицу А.2 (приложение А).
Таким образом, значение поправки
6.5.12 Такой алгоритм используется в фазово-однородных диапазонах уровня воды. Для каждого гидрометрического створа на КР должны быть выделены фазово-репрезентативные уровни воды, соответствующие скачкообразному изменению характеристик ПСР. Первый из таких уровней соответствует уровню затопления русловых мезоформ, второй - выходу воды на пойму. Необходимо так планировать измерения расходов воды, чтобы они были приурочены к диапазонам фазово-репрезентативных уровней воды. Смена фазы происходит также при достижении максимального значения половодья или дождевого паводка. Это означает, что при смене фазы режима поправку в формуле (8) следует принимать равной значению относительного отклонения лишь последнего ИРВ, вычисленного по формуле (9).
6.5.13 Необходимость корректировки, вычисленной на начало расчетных суток поправки
6.5.14 В ходе расчетов каждый очередной ИРВ Qи должен быть подвергнут объективной проверке его достоверности. Для этого согласно алгоритму, изложенному выше, вычисляется значение оперативного расхода воды на момент измерения, а именно, расход Qоп рассчитывается при использовании формулы (8) по значению среднего за время измерения уровня воды и данных двух последних ИРВ. Вычисленный расход Qоп используется для проверки достоверности вновь измеренного расхода воды Qи посредством сравнения их относительной разности
с относительной погрешностью оперативного учета стока При сравнении значений Qи и Qоп могут наблюдаться следующие случаи:
- при выполнении условия
, (25)последний ИРВ Qи признается достоверным. По его данным вычисляется относительное отклонение
- при выполнении условия
, (26)ИРВ Qи бракуется и не принимается для вычисления стока в последующие дни расчеты выполняются так, как если бы этого измерения не было. При необходимости измерение расхода воды должно быть произведено повторно;
- при выполнении условия
, (27)ИРВ признается вызывающим сомнение, и выполняется дополнительный анализ по оценке его достоверности. В случае, если в ходе анализа найдены подтверждения тому, что произошли значительные изменения ПСР, например, при изменении фазы водности с фазы подъема на фазу спада и наоборот, отклонению
6.5.15 Примеры расчетов оперативных значений расходов воды с использованием вышеизложенного алгоритма представлены в приложениях Е и Ж.
6.5.16 Выбор наилучшего способа экстраполяции рядов
6.6.1.1 Наибольшие сложности возникают при расчете оперативных значений расходов воды в весенний и осенний переходные периоды, когда измерения расходов воды практически не выполняются по соображениям техники безопасности.
6.6.1.2 К переходным периодам относятся фазы предледоставного (осеннего) ледообразования и вскрытия реки ото льда весеннего ледохода. С гидрологической точки зрения для них характерна неопределенность характеристик ПСР. С появлением ледяных образований гидравлическая система "поток - открытое русло" теряет определенность таких элементов как смоченный периметр, гидравлический радиус, шероховатость обтекаемых поверхностей и уклон свободной поверхности потока - градиент потерь энергии. При ледоходе тормозящее действие на поток оказывают забереги и льдины, остановившиеся на отмелых участках русла. Особенно сложная картина складывается при заторах, когда оценка ПСР становится крайне затрудненной. В связи с этим в эти периоды особенно необходимы измерения расходов воды, но они, как правило, не выполняются по техническим причинам и условиям безопасности.
6.6.1.3 Для вычисления стока в зимний период в практике используются переходные коэффициенты Kзим [12]
Kзим = Qзим/QУР, (28)
где Qзим - расход воды, измеренный при ледоставе, м3/с;
QУР - расход воды, полученный на основе КР или УР Q(H) при уровне, к которому отнесены ИРВ, м3/с.
Переходные коэффициенты связаны с отклонениями
, (29)поэтому для зимнего и переходных периодов может быть использован тот же алгоритм оперативного учета стока, что и при свободном русле, т.е. формула (8).
6.6.1.4 В практике для учета стока в зимний и переходные периоды используются коэффициенты: начальный Kпн (для предледоставного периода) и конечный Kвк (для периода вскрытия и очищения реки ото льда), значения которых при открытом русле принимаются равными единице. В граничные даты ледостава значения Kп и Kв устанавливаются на основе некоторых качественных соображений и даже "импровизаций", как это отмечено в работе [13].
По аналогии с коэффициентами Kпер значения относительных отклонений (начальное
При оперативном учете стока в начале осеннего переходного периода
6.6.1.4.1 Промежуточные значения
. (30)6.6.1.4.2 В период стаивания льда этот процесс по аналогии со стаиванием наледей [14] развивается более интенсивно: уменьшение массы ледяных образований
. (31)6.6.1.4.3 Таким образом для осенне-зимнего переходного периода по данным прошлых лет следует построить график связи
, а для весеннего периода стаивания график связи .В этом случае для расчета оперативных расходов воды используются аналитические зависимости
где
aп, aв - коэффициенты, полученные для зависимостей (32) и (33) по данным прошлых лет;
ti - температура воздуха в i-тый день после перехода ее через нулевое значение осенью или весной;
Для расчетов используется данные о средних суточных температурах воздуха в весенний или осенний переходные периоды.
На следующем этапе для расчета оперативных значений расходов воды в переходные периоды используется формула (8).
Примеры расчета оперативных значений расходов воды за переходные периоды с использованием данных о температуре воздуха представлены в приложении И.
6.6.2 Модель учета стока, основанная на закономерности спада зимних расходов и истощения стока
6.6.2.1 В период осеннего ледообразования и установления устойчивого ледостава часто наблюдаются заторные и зажорные явления, обусловливающие подпор уровней воды. При режимном учете стока принято выполнять срезку подпорных уровней воды, посредством соединения прямой линией точек в начале и конце подпора на хронологическом графике уровней воды. В случае оперативного учета стока такая операция невозможна, поэтому рекомендуется опираться на физические закономерности изменения стока в этот период.
6.6.2.2 Одна из особенностей гидрографа в устойчивые зимы без оттепелей - спад расходов воды по мере убывания запасов влаги в бассейне реки. На основе этой физической тенденции построена одна из моделей гидрометрического учета стока.
6.6.2.3 Уменьшение расходов воды с наступлением межени, в том числе и зимней, происходит по экспоненциальному закону, который впервые был установлен Буссинеском как частное решение линеаризованного уравнения неустановившегося движения грунтовых вод при большой мощности водоносного пласта [16]. В действительности же подземное питание рассредоточено по различным водоносным горизонтам, располагающимся на разных высотах (ярусах). Применительно к этим условиям кривая спада гидрографа стока эффективно описывается формулой, выведенной Буссинеском как решение нелинейного дифференциального уравнения фильтрации [16]
где Q0 - расход, при котором река переходит на подземное питание, м3/с;
T - время от начала ледовых явлений до момента расчета, сутки.
На самом деле уравнение (34) относится к отдельному водоносному горизонту. Исследованиями ГГИ [17] на основе анализа 1500 кривых спада для 60 гидрометрических створов на реках Сибири и Дальнего Востока установлено, что величины
, (35)где Q0 - значение расхода в последний день открытого русла перед началом ледообразования, м3/с;
QT - первый ИРВ при ледоставе, м3/с;
T - время от последнего дня при свободном состоянии русла до даты первого измерения расхода воды при ледоставе, сутки.
В случаях, когда в осенний переходный период удается выполнить измерения расходов воды методом поверхностных поплавков с помощью плывущих льдин, данные измерений используются для уточнения константы истощения грунтового питания
При низкой водности рек в предледоставный период и наличии подпорных явлений без большой потери точности можно рекомендовать при оперативном учете стока выполнять экстраполяцию значения расхода воды в последний день открытого русла, как константы, вплоть до первого измерения расхода воды при ледоставе.
Примеры расчетов оперативных расходов в период осеннего ледообразования, осуществленных с учетом закономерностей спада зимних расходов и истощения стока, показаны в приложении К.
6.7.1 Зарастание, как фактор, влияющий на изменение ПСР - одно из наиболее распространенных явлений почти во всех климатических зонах. На территории Северо-Запада, например, около 40% рек относятся к зарастающим.
6.7.2 Традиционным методом учета стока в условиях зарастающего русла остается применение переходного коэффициента Kзар = Qзар/Q0, где Qзар и Q0 - расходы воды, соответствующие зарастающему и свободному руслу при определенном уровне воды [12]. Коэффициенту Kзар приписывается роль основной характеристики ПСР, причем заведомо считается, что она снижается при зарастании, имеется в виду, что значения Kзар при этом становятся меньше единицы (Kзар < 1). Это условие и предопределяет положение КР свободного русла Q0(H) правее точек измеренных при зарастании расходов воды.
На самом же деле положение КР Q0(H) в этой зоне всегда остается неопределенным, так как именно при зарастании русла оно не бывает свободным. КР Q0(H) в зоне низких уровней, при которых, собственно, и наблюдается развитие водной растительности, строят, экстраполируя верхнюю ее ветвь, полученную до начала зарастания, и она может не отражать истинных характеристик пропускной способности свободного русла. Таким образом, Kзар в действительности является искусственной характеристикой, не отражающей реальные условия движения потока в зарастающем русле.
6.7.3 В отделе гидрометрии ФГБУ "ГГИ" разработаны более совершенные модели учета стока в условиях зарастающего русла. Одна из них - интерполяционно-гидравлическая [18], в которой в качестве характеристики изменения пропускной способности русла используется параметр Великанова m3, который определяется непосредственно по данным ИРВ
где
B - ширина русла, м;
I - уклон свободной поверхности;
n3 - коэффициент шероховатости заросшего русла.
При разработке этого метода гидрометрического учета стока было использовано свойство "параметра Великанова" сохранять относительное постоянство при однородной шероховатости русла и малых изменениях уклонов водной поверхности независимо от колебаний расходов воды. В режимном варианте учета стока выполняется интерполяция значений "параметра Великанова" в интервалах между измерениями расходов воды.
При расчете оперативных расходов воды необходимо выполнять экстраполяцию значений "параметра Великанова" в интервале от состоявшегося до следующего измерения расхода воды, при этом метод из интерполяционно-гидравлического соответственно становится экстраполяционно-гидравлическим.
Для расчета оперативного расхода воды Qt в м3/с на момент времени t в этом случае используется формула
где
Bt - ширина русла, м, полученная по зависимости B3(H) заросшего русла для значения уровня воды H на момент времени t.
6.7.4 В период зарастания может быть использован также алгоритм оперативного учета стока при нарушении однозначности зависимости Q(H), описанный в 6.5 с некоторыми дополнениями. Как уже отмечалось в пределах отдельной фазы гидрологического режима хронологическую последовательность относительных отклонений
где
- осредненные по модулю наибольшие значения Tн - продолжительность периода, в течение которого связь Q(H) неоднозначна;
При таком представлении нормированная АКФ случайного процесса, представленного косинусоидой, также может быть выражена в виде тригонометрической функции, имеющей вид
Это дает возможность оценивать АКФ временного ряда
Для периода зарастания русла в летний период хронологический ход
![]() с 18.06.2002 по 25.09.2002 включительно,
р. Олонка - г. Олонец,
Уравнение косинусоиды, представленной на рисунке 4, в виде:
. (40)6.7.5 Алгоритм оперативного учета стока при таком представлении заключается в следующем:
- по данным одного или нескольких предшествующих лет выполняется аппроксимация рядов отклонений
- вычисляются значения АКФ по формуле (39);
- вычисление оперативных расходов воды выполняется по алгоритму, описанному в 6.5 с использованием метода оптимальной экстраполяции и результатов аппроксимации ряда
6.7.6 Если в изменении пропускной способности русла в период зарастания не выявлено периодичности, согласно рисунку 5, для вычисления оперативных расходов воды в период зарастания для оперативного учета стока используется алгоритм, описанный в 6.5.
![]() в период зарастания русла, р. Сосна - г. Елец
Примеры оперативного учета стока в условиях зарастающего русла приведены в приложении Л.
7.1 Перед началом сезона оперативного учета стока по всем ГП, по которым будет выполняться расчет оперативных значений расходов воды необходимо подготовить всю мало изменяющуюся в течение сезона информацию, необходимую для проверки достоверности исходной информации об ИРВ и выполнения расчетов:
- иметь оценки относительной случайной погрешности измерения расходов воды
- выбрать в соответствии с 6.2 опорную КР, получить ее аналитическое выражение в виде УР и выполнить оценку однозначности;
- в случае, если для расчетов выбрана зависимость Q(H) одного или двух последних лет, выполнить ее экстраполяцию в соответствии с алгоритмом, изложенным в 6.4; для этой цели получить аналитическое выражение для зависимости площади поперечного сечения
- в случае, если зависимость Q(H) признана неоднозначной, необходимо установить причины неоднозначности и по данным прошлых лет получить статистические оценки рядов
- перед началом расчетов иметь совокупность состоявшихся в текущем году ИРВ (если измерения состоялись) и последних ИРВ предшествующего года;
- по данным прошлых лет выполнить все необходимые проработки, чтобы иметь возможность использовать алгоритмы, описанные в 6.6 и 6.7, для оперативного учета стока в зимний и переходные периоды, а также при зарастании русел.
7.2 Текущая работа по оперативной обработке гидрологических данных и учету стока заключается в следующем:
- получение из базы данных ЦСД данных срочных наблюдений за уровнями воды как АГП, так и ГП, на которых наблюдения выполняются ручным способом;
- проверка качества данных наблюдений за уровнями воды и, при необходимости, их корректировка в соответствии с 5.3 и 5.4;
- формирование таблицы проверенных и принятых для расчета срочных данных наблюдений за уровнем воды за предшествующие и текущие сутки; занесение исходных данных в таблицы оперативной обработки гидрологических данных и оперативного учета стока А.1 и А.2 (приложение А);
- формирование таблицы ИРВ в виде таблицы А.4 (приложение А);
- непосредственное вычисление оперативных расходов воды.
8 Ежедневная обработка данных об уровнях и расходах воды в оперативном режиме и получение обобщенных характеристик
8.1 Порядок ежедневной оперативной обработки гидрологических данных различается в зависимости от того автоматизированы или нет наблюдения за уровнем воды. При оперативном учете стока необходимо учитывать следующее:
- если учет стока ведется по данным ручных наблюдений, вычисление оперативных расходов воды выполняется только в стандартные сроки наблюдений (8 и 20 ч) в обычном режиме и дополнительно в нестандартные сроки в период паводков и половодий. Таблица оперативного учета стока в форме А.2 (приложение А) заполняется только для этих сроков. Вычисленные оперативные значения расходов воды заносятся также в таблицу оперативной обработки данных формы А.1 (приложение А);
- по данным наблюдений за уровнями воды, полученным с АГП, расчет оперативных значений расходов воды выполняется для всех наблюденных значений уровней воды. Таблица А.2 (приложение А) заполняется для всех сроков. В таблицу А.1 (приложение А) заносятся только данные в стандартные сроки;
- для контроля и визуализации данных каждое вновь полученное (при необходимости откорректированное) значение уровня воды и вычисленное значение оперативного расхода воды наносится на совмещенный хронологический график.
8.2 Порядок ежедневной оперативной обработки уровней воды осуществляется по 5.5. Аналогичный порядок обработки и обобщения данных выполняется и по данным оперативного учета стока. По данным об оперативных расходах воды необходимо получить значения следующих элементов оперативного учета стока:
- срочный расход воды за 8-ми часовой срок местного времени текущих суток;
- срочный расход воды за 20-ти часовой срок местного времени предшествующих суток;
- среднесуточный расход воды за предшествующие сутки;
- наименьший расход воды за предшествующие сутки;
- срок наблюдения наименьшего расхода воды за предшествующие сутки;
- наибольший расход воды за предшествующие сутки;
- срок наблюдения наибольшего расхода воды за предшествующие сутки.
8.3 По окончании декады и месяца вычисляются обобщенные за эти периоды характеристики уровней (средние, высшие и низшие) и расходов воды (средние, наибольшие и наименьшие), которые помещаются в таблицу формы А.6 (приложение А).
В графе "Примечание" таблицы А.6 (приложение А) делаются отметки о наличии дополнительных факторов, влияющих на изменение ПСР: заторах, зажорах, переменном подпоре, попусках вышерасположенной ГЭС и т.п.
Информацию о средних и характерных за многолетний период расходов и уровней воды за соответствующую декаду или месяц в таблице А.6 (приложение А) формируют специалисты-гидрологи.
В последней строке таблицы помещаются отклонения данных текущего года от их многолетних характеристик.
(рекомендуемое)
ДЛЯ ОПЕРАТИВНОГО УЧЕТА СТОКА
обработки гидрологических данных в оперативном режиме
результатов сравнения данных АГК и стандартных наблюдений
обобщенных характеристик оперативного учета стока
(обязательное)
ОТ УРОВНЕЙ ВОДЫ
Б.1 Общие положения построения алгоритмов аппроксимации зависимостей расходов от уровней воды
Б.1.1 Гидравлические и геометрические элементы руслового потока в конечном счете зависят от положения уровня воды H. Поэтому в гидрометрии принято использовать зависимость расходов Q от уровней воды H, в аналитической форме или в виде так называемой кривой расходов Q = f(H). Хотя с точки зрения физики процесса независимой переменной в этой связи является расход воды, обратная зависимость Q(H) удобна тем, что она позволяет вычислить расход воды непосредственно по заданному значению уровня.
Б.1.2 Различают однозначные и неоднозначные зависимости расходов от уровней воды. Понятие однозначности связи Q(H) и алгоритм ее оценки изложены в приложении В. В общем случае зависимость Q(H) не может считаться однозначной. Следует исходить именно из этой предпосылки с тем, чтобы исследовать факторы, нарушающие однозначность связи Q(H), и обосновать, с одной стороны, наиболее эффективные методы учета стока, с другой - требования к точности измерений расходов воды.
Однозначная зависимость Q(H) наблюдается в гидростворах с устойчивыми руслами при равномерном движении потока. В некоторых гидростворах она может сохраняться в течение нескольких лет (трех и более). В этом случае ее принято называть многолетней.
Б.1.3 Для получения зависимостей расхода воды Q(H) и его элементов: площади живого сечения
Кроме того, для анализа характера связи Q(H) привлекаются сведения об ИРВ и наблюденных при этом уровнях, помещаемых в книжках КГ-3М, в которых содержатся данные о состоянии русла и влияющих на изменения его ПСР. В исходную совокупность следует включать все данные, полученные при свободном состоянии русла, при отсутствии зарастания, переменного подпора и прочих влияющих на изменение ПСР факторов, в том числе следующие ряды данных:
- ИРВ Qи, м3/с;
- наблюденные при измерениях расходов воды уровни Hи, м;
- зафиксированные при измерениях площади живого сечения
- средние скорости течения в гидростворе во время измерений расходов воды Vи, м/с.
Б.1.4 В графическом варианте связь расходов и уровней воды представляется КР в поле точек (Q, H), которая строится на миллиметровке (проводится на глаз посередине полосы рассеивания точек) [12]. Такой способ представления зависимости Q(H) при современном развитии вычислительных методов не может быть рекомендован. Вопросы аналитической интерпретации гидрометрических данных к настоящему времени достаточно разработаны. Ниже излагаются основные методы, рекомендуемые для аналитической аппроксимации зависимости расходов от уровней воды, равно как и для аналитического представления зависимостей элементов расхода (площадей водного сечения
Б.1.5 При аналитической аппроксимации зависимостей Q(H),
- стандарта рассеяния аналитической зависимости
;- множественного коэффициента корреляции R, которые определяются по формулам.
Стандарт рассеяния зависимости Y(H) определяется по формуле
где Yи - измеренная характеристика, для которой подбирается аналитическое УР;
Y(Hи) - характеристика, вычисленная по УР для наблюденного при измерении уровня Hи;
n - количество ИРВ;
p - число степеней свободы подбираемого уравнения.
Выборочный множественный коэффициент корреляции R определяется как коэффициент между Yиi и Y(Hиi) по формуле
Для аппроксимации Q(H) рекомендуется применять представленные в 6.2 аналитические выражения (10) и (11), из которых должно быть выбрано наиболее эффективное (оптимальное).
Для удобства представления коэффициентов аналитических уравнений при выполнении расчетов по МНК значения уровней рекомендуется выражать не в сантиметрах над "нулем" поста, как это принято в гидрологии, а в метрах над "нулем" поста.
Б.2.1 Предложенная В.Г. Глушковым аналитическая форма аппроксимации КР имеет вид
где Q - расход воды, м3/с;
H - наблюденный уровень воды, м;
H0 - уровень, при котором расход воды в исследуемом гидростворе равен нулю, м;
a и b - параметры УР, определяемые по данным измерений.
Б.2.2 Для получения параметров a и b выражение (Б.3) преобразуется к логарифмическому виду:
ln (Q) = ln (a) + b ln (H - H0). (Б.4)
Б.2.3 Далее на основе совокупности ИРВ, преобразованной к виду (Б.3) МНК получают параметры ln(a) и b, затем уравнение вновь преобразуют к виду (Б.3). Хотя понятие об уровне нулевого расхода H0 имеет вполне явный гидравлический смысл, его определение в действительности затруднено и требует дополнительных морфометрических данных, которые, как правило, отсутствуют в исходных материалах. В связи с этим величина H0 на практике оказывается таким же статистическим параметром, как и a или b и определяется методом решения оптимизационной задачи: принимается такое значение H0, при котором рассеяние связи Q(H) наименьшее.
При расчетах следует задавать некоторый диапазон изменения H0. В качестве верхней границы этого диапазона рекомендуется принимать значение низшего наблюденного в исследуемом гидрометрическом створе уровня воды - Hнизш. Нижняя граница определяется вычитанием из Hнизш максимальной глубины потока hmax, соответствующей низшему наблюденному уровню воды. Таким образом, H0 располагается в диапазоне:
(Hнизш - hmax) <= H0 <= Hнизш. (Б.5)
Выбирается значение H0, отвечающее лучшим характеристикам аппроксимации связи Q(H): наименьшему стандарту рассеяния УР и наибольшему значению множественного коэффициента корреляции. Значения параметров УР оцениваются при использовании МНК.
Б.2.4 Парабола Глушкова не всегда пригодна для аппроксимации КР во всем диапазоне изменения уровня воды. Лишь в отдельных гидрометрических створах с правильной параболической формой поперечного сечения при отсутствии выраженного изменения крутизны боковых склонов русла по всей амплитуде колебания уровней воды аппроксимация Q(H) в виде уравнения (Б.3) может давать хорошие результаты.
Б.3 Алгоритм II. Аппроксимация методом полиномиальной регрессии
Б.3.1 В соответствии с алгоритмом II аналитически КР представляется в виде уравнения:
Исходная совокупность данных включает ряды:
- ИРВ Qи;
- наблюденные при измерениях уровни воды Hи;
- H2и - значения уровней воды в квадрате;
- Hmи - значения уровней воды в степени m.
Использование при аппроксимации КР степеней полинома вида (Б.6) выше 4-ой приводит, как правило, к плохой обусловленности матрицы (H, Hm), что многократно увеличивает погрешности расчета параметров УР.
Б.3.2 Параметры УР вида (Б.6) определяются с использованием МНК для полиномов 2-ой, 3-ей, 4-ой степени, а также неполных многочленов, например:
Q(H) = a0 + a1H2 + a2H3;
Q(H) = a0 + a1H + a2H3; (Б.7)
и т.д.
Для аппроксимации выбирается уравнение, имеющее наибольший множественный коэффициент корреляции (Б.2) и наименьший стандарт рассеяния (Б.1) связи Q(H).
Б.3.3 Исследования показывают, что уравнение вида (Б.6) применимо для аппроксимации связи Q(H) лишь в ограниченном числе случаев, а именно: оно хорошо аппроксимирует зависимость Q(H), если диапазон изменения уровня воды на гидростворе невелик (например, в маловодные годы) и значения минимальных и максимальных расходов воды различаются друг от друга не более, чем на порядок. В случаях, когда не удается получить приемлемую аппроксимацию КР методом обычной полиномиальной регрессии, следует использовать другие способы, например, полиномиальную регрессию с ограничениями или кусочно-полиномиальную регрессию.
Б.4 Алгоритм III. Аппроксимация зависимости Q(H) методом полиномиальной регрессии с введением ограничений в экстремальных точках амплитуды уровня воды
Б.4.1 Для того чтобы получить удовлетворительную аппроксимацию во всем диапазоне изменения уровней воды, а также для выполнения условия выпуклости кривой dQ/dH > 0 предлагается использовать метод кусочно-полиномиальной регрессии и его частный случай - метод полиномиальной регрессии с заданием ограничений на значение функции и ее производных в точке одного из экстремумов аргумента [8]. В соответствии с этим способом расчеты выполняются по следующему алгоритму.
Б.4.2 Первоначально для аппроксимирующей функции задаются ограничения в точке одного из наблюденных экстремумов уровня воды (для КР расходов - в точке низшего наблюденного уровня воды, для кривых площадей и средних скоростей течения - наивысшего наблюденного уровня воды). Для этого необходимо установить положение КР в нижней, а кривых площадей и средних скоростей течения в верхней частях диапазона изменения уровня воды, что может быть выполнено посредством графической или аналитической аппроксимации (с использованием МНК) данных ИРВ в ограниченном диапазоне изменения уровней воды. Оценку граничных значений рекомендуется производить по данным многолетних измерений. Чем больше ИРВ в нижней части КР будет использовано для определения граничного условия, тем с меньшей погрешностью будет определено положение КР при Hнизш.
Для оценки Qгр желательно получить уравнение аппроксимирующей зависимости в виде алгебраического полинома, что облегчает вычисление граничных условий: значений расхода воды и производных функции Q(H) в точке низшего наблюденного уровня воды. Аппроксимирующая функция, удовлетворяющая заданным ограничениям, записывается в следующем виде
, (Б.8)где
g(H) - полином, подбираемый по МНК.
В порядке первого приближения первым двум полиномам придается вид
, (Б.9) . (Б.10)Б.4.3 Подбор аналитической формы полинома g(H), производится на основе МНК с учетом тех же критериев качества, которые описаны выше. При этом предварительно осуществляется преобразование исходных данных об ИРВ следующим образом
, (Б.11)где Qиi - i-ый ИРВ,
Hиi - соответствующий ему уровень воды.
Таким образом, расчет выполняется по совокупности преобразованных данных в виде Yi и соответствующих им уровней Hи.
Б.4.4 После подбора полинома g(H) окончательное УР Q(H) для всего диапазона уровней воды приобретает вид
Q(H) = Q(Hнизш) + (H - Hнизш)(a0 + a1H +
+ a2H2 + + amHm). (Б.12)
Б.4.5 Качество аппроксимации Q(H), V(H),
![]() а - регрессия с ограничениями; б - полиномиальная регрессия.
в створе ГП р. Дон - х. Беляевский, за 1971 год
Б.4.6 Такая картина может наблюдаться из-за резкого изменения ПСР в некоторых диапазонах уровней воды, например, при выходе воды на пойму. В практике гидрологических исследований в таких случаях выделяются как минимум две фазовые кривые (меженная и паводочная). При этом неизбежно возникает вопрос об их сопряжении. При использовании метода кусочно-полиномиальной регрессии эта задача решается на основе алгоритма III. Зависимость Q(H) в виде параболы Глушкова (Б.3) или полинома (Б.6), полученная на первом шаге аппроксимации, принимается для представления связи Q(H) в нижней части диапазона уровней воды, где ее можно считать удовлетворительной. Далее, вышележащую амплитуду уровня воды следует разбить на отрезки для кусочного представления аппроксимирующей функции. Эти отрезки должны соответствовать морфометрии русла. Однако, в этом случае необходимо было бы иметь наибольшее количество ИРВ именно в окрестности точек, ограничивающих заданные отрезки, а именно, при переходе уровня воды через точки резкого изменения формы поперечного сечения русла.
В верхней точке Hгр1 нижнего отрезка кривой, уравнение, для которого уже получено, задаются ограничения, согласно алгоритму III, т.е. по УР вычисляется значение расхода воды Q(Hгр1). Затем по данным ИРВ, попадающим в диапазон, расположенный выше граничной точки Hгр1, подбирается аппроксимирующий полином. Для получения удовлетворительной аппроксимации связи Q(H) по всей амплитуде уровня воды может быть сделано несколько таких разбиений в соответствии с морфометрическими особенностями русла или поймы при их различных наполнениях (изменения крутизны берегов, рельефа русла, характера растительности).
Б.4.7 При аппроксимации зависимостей
При аппроксимации зависимостей средней скорости течения от уровня воды V(H), как годовых, так и многолетних, целесообразнее ограничения на функцию задавать в точке Hвысш. Это объясняется тем, что малые скорости течения измеряются с более низкой точностью, что снижает надежность оценки положения кривой V(H) в ее нижней части.
Б.4.8 Пример использования алгоритма кусочно-полиномиальной аппроксимации приводится в приложении Д.
Б.5 Выбор оптимального алгоритма аппроксимации
Б.5.1 Известно, что одна и та же кривая на каком-то конкретном отрезке может быть описана самыми различными аналитическими выражениями. Рациональный выбор того или иного уравнения должен опираться не только на статистические оценки и критерии, но еще учитывать и другие характеристики вычислительного алгоритма. Решение об окончательном выборе алгоритма вычисления ЕРВ должен принимать специалист-гидролог на основе аттестации выбранного алгоритма учета стока.
Б.5.2 Применительно к условиям каждого конкретного гидроствора выбор оптимального алгоритма аппроксимации зависимостей Q(H),
Б.5.3 Характеристики точности предназначены для оценивания погрешностей результатов, получаемых при использовании выбранного алгоритма. Они содержат как методические, так и трансформированные погрешности, обусловленные неточностью исходных данных и погрешностями алгоритмов. К основным характеристикам точности относятся:
- статистические показатели для случайных погрешностей (среднее квадратическое значение отклонений (СКО); смещение оценок - математическое ожидание отклонений);
- границы погрешностей, как случайных, так и систематических, которые могут быть доверительными или детерминированными.
Б.5.4 Перечень характеристик, которыми следует руководствоваться при выборе оптимальной аппроксимирующей функции Q(H), должен включать:
- множественный коэффициент корреляции (должен быть наибольшим);
- дисперсия и СКО абсолютных остатков регрессии Dq и
- дисперсия и СКО относительных остатков регрессии
- значения математического ожидания (среднего арифметического) остатков регрессии, как абсолютных, так и относительных, должны быть близки к нулю (
Б.5.5 Адекватность полученной модели исходным данным проверяется также на основании анализа графиков остатков регрессии, представленных на рисунке Б.1. При этом могут быть обнаружены следующие типы отклонений:
- наличие выделяющихся измерений;
- систематические отклонения остатков на отдельных отрезках;
- непостоянство дисперсии по амплитуде изменения расходов и уровней воды.
При оптимальной аппроксимации связи Q(H) относительные остатки регрессии на графике зависимости от H или Qрасч должны располагаться симметричными полосами по обе стороны от ординаты H или Q, (см рисунок Б.1). Рекомендуется анализировать поведение остатков в различные фазы режима, например, в зависимости от времени t или от порядкового номера ИРВ n, как показано на рисунке Б.2б. Подобный анализ позволит выявить периоды или фазы режима, в которые нарушается однозначность связи Q(H). На рисунке Б.1б приведен пример, для которого в верхнем диапазоне уровней воды все отклонения
![]() а - МКР Q(H),
б - график относительных остатков регрессии
в створе ГП р. Енисей - г. Енисейск, 1960 - 1968 годы
Б.5.6 По степени сложности все описанные выше алгоритмы (I - III) достаточно близки и физически обоснованы, а вот что касается физической интерпретируемости, то вопреки ранее существовавшему представлению, парабола Глушкова в большинстве случаев не заслуживает предпочтения из-за неопределенности уровня нулевого расхода, H0.
Б.5.7 При необходимости аналитическое представления зависимости Q(H) в виде (Б.3) или (Б.6) может быть получено по данным таблицы координат многолетней или годовой кривой расходов воды или площадей поперечного сечения. В этом случае вместо данных совокупности ИРВ (Hи, Qи) в расчетах участвуют данные в узлах таблицы координат (H, Q).
Б.5.8 Описанные алгоритмы аппроксимации зависимостей расходов воды Q(H), площадей живого сечения
(обязательное)
ОТ УРОВНЕЙ ВОДЫ
В.1 В общем случае КР или УР расходов от уровней воды отражают осредненные за исследуемый период характеристики ПСР и поэтому УР в виде (10), (11) правомерно принять в качестве детерминированной характеристики ПСР, на основе которой можно получить значение расхода воды в первом приближении. Отклонения ИРВ исходной совокупности от УР оказываются обусловленными двумя причинами: случайными погрешностями измерения расходов воды и изменениями условий протекания потока под воздействием динамических и морфологических факторов. Количественная оценка совокупного влияния этих факторов на изменение ПСР возможна лишь на основе совместного решения уравнения неустановившегося движения потока и уравнения деформации русла. На практике такое решение затруднено недостаточностью необходимых данных. Вместе с тем изменения ПСР, обусловленные неучтенными факторами, можно рассматривать как случайную функцию времени и оценивать ее с помощью статистических методов [5], [6].
В.2 В качестве характеристики изменения ПСР могут служить относительные отклонения
, (В.1)где Qи, и Q(Hи) - ИРВ и расход воды, полученный по связи Q(H) для уровня Hи, зафиксированного во время измерения расхода воды Qи, соответственно, м3/с.
В.3 Для оценки однозначности зависимости Q(H) рекомендуется использовать критериальное отношение дисперсий F Фишера [19]. Связь Q(H) признается однозначной, если
где
Fкр - критическое значение отношения дисперсий Фишера при заданном уровне значимости.
Величина
. (В.3)В.4 Погрешность этой оценки в значительной степени зависит от объема исходной выборки. Относительное значение погрешности измерения расхода воды
В.5 Критические значения отношения Фишера для 5%-ного уровня значимости и различного числа измерений n при условии, если принято, что длина выборки для оценки
Таблица В.1
для 5%-ного уровня значимости
В.6 Согласно предложенному критерию для 5%-ного уровня значимости и числа измерений n = 31 значение Fкр равно 1,46. Это означает, что зависимость Q(H) можно признать однозначной, если
.В.7 Критериальное отношение (В.2) дает возможность судить о близости или различии дисперсий
В.8 Поэтому наряду с критериальным отношением (В.2) для оценки однозначности КР рекомендуется использовать непараметрический критерий оценки случайности временных рядов, в частности, наиболее простой для этого случая - критерий оценки рядов на случайность с учетом знаков отклонений [20]. Выполняется анализ временных рядов отклонений
В.9 Временной ряд
В.10 Зависимость Q(H) признается однозначной в тех случаях, когда отношение дисперсий удовлетворяет критерию Фишера при заданном уровне значимости, а относительные отклонения ИРВ от КР или УР носят случайный характер. Проверка временных рядов
Такая процедура должна осуществляться при оценке однозначности зависимости расходов воды и от других, кроме уровней воды, определяющих факторов. Если связь Q(H) на основе соответствующих УР признается неоднозначной, устанавливаются причины ее обусловившие и в соответствии с ними выбирается адекватный способ учета стока.
В.11 В случае нарушения неравенства (В.2) и неслучайности группировки последовательности
В.12 На основе анализа временных рядов отклонений
(рекомендуемое)
ЗАВИСИМОСТИ РАСХОДОВ ОТ УРОВНЕЙ ВОДЫ
Г.1 Рассмотрим алгоритм выбора аналитической зависимости, аппроксимирующей МКР для ГП м. ст. Кантегир на реке Кантегир (приток Саяно-Шушенского водохранилища).
Г.2 Для получения аналитической зависимости расходов от уровней воды использовались данные наблюдений за многолетний период по 2011 год включительно. Для исследуемого створа был выполнен анализ многолетних данных об ИРВ и получены зависимости расходов от уровней воды и их аналитические представления, рекомендуемые для расчета оперативных значений расходов воды. Анализ многолетних данных показал, что для ГП р. Кантегир - м. ст. Кантегир все годовые зависимости Q(H) остаются однозначными за весь период с 2000 по 2011 год. (рисунок Г.1).
![]() р. Кантегир - м. ст. Кантегир
Г.3 Наилучшая аппроксимация МКР получена методом полиномиальной регрессии с ограничениями в виде полинома 4-й степени. Его параметры приведены в таблице Г.1. Ограничение задано на значение расхода воды при уровне Hнизш = 4,13 м над "нулем" поста. Качество аппроксимации оценивается параметрами, представленными в последних трех графах таблицы Г.1.
Таблица Г.1
ГП р. Кантегир - м. ст. Кантегир за 2000 - 2011 годы
Основной предпосылкой для выбора метода гидрометрического учета стока служит оценка однозначности зависимости расходов от уровней воды, которая выполняется на основе сравнения дисперсии относительных отклонений ИРВ от осредняющей кривой и дисперсии случайных погрешностей измерения расходов воды по алгоритму, описанному в приложении В. В соответствии с этим алгоритмом полученная МКР признана однозначной. При длине совокупности ИРВ n = 178 дисперсия относительных отклонений составила
, дисперсия случайных погрешностей измерения расходов воды при свободном русле в створе ГП р. Кантегир - м. ст. Кантегир (рекомендуемое)
ЗАВИСИМОСТИ РАСХОДОВ ОТ УРОВНЕЙ ВОДЫ
Д.1 Рассмотрим пример использования алгоритма кусочно-полиномиальной аппроксимации для получения аналитической зависимости Q(H). С этой целью возьмем МКР для створа ГП р. Лаба - х. Догужиев за период с 2003 по 2006 год. Удовлетворительную аппроксимацию Q(H) одной кривой во всем диапазоне изменения уровня воды получить не удалось, как видно из рисунка Д.1.
![]() за 2003 - 2006 годы
Д.2 В связи с этим следует выделить два диапазона уровней, для каждого из которых подбирается своя полиномиальная зависимость. Первоначально по многолетним данным в граничной точке низшего наблюденного уровня воды Hнизш = 210 см (2,1 м) определено значение расхода воды, равное 15,0 м3/с, которое и было принято в качестве ограничения на положение кривой Q(H) в этой точке.
Д.3 Методом полиномиальной регрессии с ограничениями было получено уравнение Q(H) в диапазоне от 2,4 до 3,20 м над "нулем" поста, параметры которого приведены в первой строке таблицы Д.1. Затем по полученному УР вычислено значение расхода воды в точке H = 3,20 м над "нулем" поста (Q3,20 = 156 м3/с), которое было задано как ограничение при аппроксимации кривой в диапазоне от 3,21 до 4,95 м над "нулем" поста. Далее методом полиномиальной регрессии с ограничениями было получено уравнение Q(H) в диапазоне от 3,21 до 4,95 м над "нулем" поста, параметры которого приведены во второй строке таблицы Д.1.
Таблица Д.1
Зависимость Q(H), полученная методом кусочно-полиномиальной регрессии с использованием алгоритма, описанного в приложении Б, представлена на рисунке Д.2. Ее качество оценивается параметрами, представленными в последних трех графах таблицы Д.1.
![]() от уровней в створе р. Лаба - х. Догужиев, полученная
методом кусочно-полиномиальной регрессии
Д.4 На рисунке Д.3 приведен график распределения относительных отклонений
![]() 1 - для МКР, аппроксимирующей точки диапазона
от 2,42 до 3,20 м, 2 - для МКР, аппроксимирующей точки
диапазона от 3,21 до 4,95 м
относительных отклонений
ГП р. Лаба - х. Догужиев за период с 2003 по 2006 год
(рекомендуемое)
В УСЛОВИЯХ СВОБОДНОГО СОСТОЯНИЯ РУСЛА
Е.1 Расчет оперативных значений расходов воды с использованием относительных отклонений ИРВ текущего года от принятой опорной КР рассмотрим на примере створа ГП станица (ст.) Ладожская на реке Кубань. На рисунке Е.1 представлена МКР для этого створа за период с 1994 по 1999 год. Как видно, эта зависимость может считаться однозначной и аппроксимируется кривой, описываемой уравнением
Q = -0,5948H3 + 23,426H2 - 53,057H - 75,829 (Е.1)
Диапазон уровней воды от 3,7 до 7,05 м над "нулем" поста.
![]() за период 1994 - 1999 годы
Е.2 Эта зависимость используется для расчетов оперативных значений расходов воды за период открытого русла в 2000 году.
Е.3 В таблице Е.1 представлен фрагмент расчета оперативных значений ЕРВ за период открытого русла 2000 года.
Таблица Е.1
в створе ГП р. Кубань - ст. Ладожская
за период открытого русла 2000 года
Е.4 Относительные погрешности оперативного учета стока для примера, представленного в таблице Е.1, не превышают 10%. Случайная средняя ошибка по итогам 2000 года составила 10%. Полученные в результате расчетов режимные и оперативные гидрографы стока показаны на рисунке Е.2.
![]() 1 - режимный; 2 - оперативный; 3 - измеренные расходы воды
ст. Ладожская за период свободного русла 2000 года
(рекомендуемое)
С ИСПОЛЬЗОВАНИЕМ МЕТОДА ОПТИМАЛЬНОЙ ЭКСТРАПОЛЯЦИИ
Ж.1 Рассмотрим способ оперативного учета стока в условиях деформации русла с использованием оптимальной экстраполяции рядов относительных отклонений
Ж.2 МКР Q(H) в створе ГП за предшествующие годы представлена на рисунке Ж.1 уравнение МКР для ГП р. Оскол - с. Ниновка имеет вид:
Q(H) = 4,4308 + 28,805H - 0,9984H2 + 1,7847H3. (Ж.1)
![]() Параметры МКР приведены в таблице Ж.1.
Таблица Ж.1
Ж.3 Для оценки характеристик АКФ по формуле (9) рассчитываются относительные отклонения от МКР измеренных расходы воды за различные по водности годы. Пример такого расчета приведен в таблице Ж.2.
Таблица Ж.2
в створе ГП р. Оскол - с. Ниновка
Ж.4 По данным о среднегодовых расходах воды предшествующих лет определяется величина среднего расхода воды за многолетие, и в совокупности лет выделяются маловодные, средние по водности и многоводные года. Составляются таблицы относительных отклонений для каждой группы лет по форме таблицы Ж.2, после чего для каждого из трех рядов относительных отклонений производится расчет АКФ. На рисунке Ж.2 представлены графики АКФ для групп лет разной водности.
![]() 1 - маловодных; 2 - средних по водности; 3 - многоводных
для групп лет разной водности
Ж.5 Для каждого из трех рядов относительных отклонений вычисляются значения меры погрешности по формуле (17).
Для групп маловодных, средних по водности и многоводных лет значения дисперсии рядов относительных отклонений
, , .Ж.6 Расчет оперативных значений расходов воды р. Оскол в створе с. Ниновка за период свободного русла 2010 года представлен в таблице Ж.3.
Таблица Ж.3
в створе с. Ниновка за период свободного русла
На каждую дату в графу 2 заносятся значения средних суточных уровней воды. В графу 3 на даты измерений заносятся значения ИРВ. В графу 4 помещаются значения Q(Hи), рассчитанные по МКР для значений средних суточных уровней воды. Параметры МКР представлены в таблице Ж.1. В графу 5 занесены значения относительных отклонений
Ж.7 Для расчета весовых коэффициентов оптимальной экстраполяции необходимо определить интервалы
а) время в сутках от даты t0, на которую производится расчет оперативного значения расхода воды до даты последнего ИРВ t1
; (Ж.2)б) время в сутках от даты последнего ИРВ t1 до даты предпоследнего ИРВ t2
; (Ж.3)в) время в сутках от даты t0, на которую производится расчет оперативного значения расхода воды до даты предпоследнего ИРВ t2
. (Ж.4)Ж.8 Следующим этапом является расчет АКФ для каждого значения
Водность года должна определяться по прогнозу. В начальный период расчета водность определяется по предыдущему году. В нашем случае в начальный период год признан маловодным (первоначально по предыдущему 2009 году, а затем - средним по водности по данным расчетного 2010 года).
Ж.9 Весовые коэффициенты p1 и p2 определяются по формулам (18) и (19) оптимальной экстраполяции с учетом значений АКФ
Ж.10 Выполняется оценка математического ожидания относительных отклонений
Ж.11 Вычисляется экстраполированное значение относительного отклонения по формуле (18), которое используется для расчета оперативных значений расходов воды вплоть до следующего измерения расхода воды.
Иногда в расчетах приходится опираться только на одно последнее значение ИРВ. Это случается, например, при редких измерениях расходов воды, которые обусловливают малую корреляционную связанность отклонений от УР. В этом случае необходимые характеристики вычисляются по формулам (23) и (24).
Ж.12 Расчет оперативных значений расходов воды на каждые сутки (графа 17) производится по формуле (8).
Ж.13 Относительные погрешности ЕРВ (графа 18) оцениваются путем сравнения рассчитанных оперативных значений расходов воды и режимных данных о ЕРВ, полученных по истечении года (графа 16) по формуле
, (Ж.5)где
Qоп t - значение оперативного расхода воды, рассчитанного на дату t;
Qр t - значение ЕРВ на дату t, полученное по истечении года (из режимных данных).
Случайная среднеквадратическая погрешность оперативного учета за весь период рассчитывается по формуле
, (Ж.6)где N - количество членов ряда.
Случайная среднеквадратическая погрешность оперативного учета стока р. Оскол в створе с. Ниновка за период свободного русла 2010 г. составила 6,5%.
Ж.14 Пример оперативного учета стока за период свободного русла 2010 года для ГП р. Оскол в створе с. Ниновка представлен в таблице Ж.3.
На рисунке Ж.4 представлены режимный и оперативный гидрографы р. Оскол в створе с. Ниновка за период свободного русла 2010 г.
![]() 1 - режимный; 2 - оперативный; 3 - измеренные расходы воды
за период свободного русла 2010 года
(рекомендуемое)
В ОСЕННИЙ И ВЕСЕННИЙ ПЕРЕХОДНЫЕ ПЕРИОДЫ
И.1 Оперативный учет стока в осенне-зимний переходный период
И.1.1 Методика расчета оперативных значений расходов воды в переходные периоды с использованием зависимости изменения ПСР от накопившейся суммы температуры воздуха, изложенная в 6.6.1, была опробована для ряда ГП на реках Сибири.
И.1.2 Рассмотрим подробнее алгоритм расчета для ГП с. Локоть на реке Алей (Алтайский ЦГМС) для осенне-зимнего периода 2005 года, результаты расчета по которому представлены в таблице И.1.
Таблица И.1
за осенне-зимний период 2004 года и расчеты, необходимые
для построения зависимости (32)
И.1.3 В качестве опорной для оперативного учета стока в 2005 году принята КР 2004 года, аппроксимированная аналитической зависимостью
И.1.4 Для получения расчетных зависимостей использовались данные наблюдений в осенне-зимний период 2004 года, представленные в графах 1, 2, 3, 4, 7 таблицы И.1.
И.1.5 По УР (И.1) рассчитывались расходы, занесенные в 5-ю графу таблицы И.1. Это же УР в дальнейшем использовалось и для расчетов оперативных расходов воды за осенне-зимний период 2005 года.
И.1.6 Относительные отклонения, занесенные в 6-ю графу таблицы И.1, рассчитывались по формуле (9) для каждого ИРВ.
И.1.7 Нарастающие суммы отрицательных температур рассчитывались, начиная от первой отрицательной температуры, предшествующей появлению ледяных образований и заносились в 8-ю графу таблицы И.1. В графе 9 производились расчеты, необходимые для построения зависимости (32).
И.1.8 Для построения линейной зависимости по формуле (32) необходимо выбрать две точки. Первая точка должна соответствовать относительному отклонению
В качестве второй точки для расчета зависимости (32) берется относительное отклонение, рассчитанное для расхода воды, измеренного уже при ледоставе (16 декабря). В результате построения получаем линейную зависимость в виде прямой линии (далее - прямая), представленную на рисунке И.1, которая описывается уравнением (И.2).
![]() от корня квадратного из модуля нарастающей суммы
отрицательной температуры воздуха для ГП р. Алей - с. Локоть
за осенне-зимний период 2004 года
И.1.9 Зависимость (И.2) используется для расчетов оперативных значений расходов воды за осенне-зимний период 2005 года.
Таблица И.2
Результаты расчета оперативных значений расходов воды для ГП
р. Алей - с. Локоть за осенне-зимний период 2005 года
И.1.10 Поэтапный алгоритм расчета, как если бы он производился в режиме реального времени, представлен в перечислении а) - е):
а) по уравнению КР за предыдущий 2004 год после получения срочных значений уровней воды рассчитываем срочное значение оперативного расхода воды Q(H), которое записывается в графу 5 таблицы И.2;
б) после каждого очередного ИРВ следует вычислить значение относительного отклонения
; (И.3)в) с момента появления отрицательных среднесуточных температур воздуха начинаем рассчитывать их нарастающие суммы (в нашем случае это 8 ноября). Если ледяные образования не появляются, а среднесуточные температуры остаются положительными на протяжении нескольких суток, расчет начинают заново. До появления ледяных образований оперативные значения расходов воды рассчитываются по алгоритмам, изложенным в 6.3, 6.4, 6.5.
С того дня, когда замечены первые ледяные образования (16 ноября), относительные отклонения
г) со второго дня наличия ледяных образований (17 ноября) оперативные значения расходов воды можно получать с использованием
; (И.4)д) расчет по этому алгоритму можно вести до следующего измерения расхода воды (в нашем случае это 30 ноября);
е) с получением очередного ИРВ 30 ноября и рассчитанного для него относительного отклонения можно построить новую зависимость по формуле (15), используя данные уже текущего года. Для построения уточненной зависимости используем значения
, равные соответственно 0 и 10. По этим данным получена линейная зависимость в виде прямой 1 на рисунке И.2, которая описывается уравнением . (И.5)Это уравнение можно использовать для расчета оперативных значений расходов воды с 1 декабря вплоть до следующего измерения расхода.
![]() от корня квадратного из модуля нарастающей суммы
отрицательной температуры воздуха для ГП р. Алей - с. Локоть
за осенне-зимний период 2005 года
И.1.11 После получения следующего ИРВ (15 декабря) зависимость можно уточнить, добавив еще одну точку на график и получить уточненную зависимость в виде прямой 2 (см. рисунок И.2), которая описывается уравнением (И.6). До следующего измерения расхода воды для расчетов можно использовать уточненное уравнение (И.6).
После следующего измерения (18 декабря) добавляем на график связи еще одну точку и получаем новое уточненное уравнение
. (И.7)Зависимость, представленная на рисунке И.2, существенно не изменялась в период с 1 по 18 декабря (прямые 1, 2, 3), что подтверждает правомерность принятого решения об использовании дополнительного фактора - среднесуточной температуры воздуха (см. 6.6.1.4.1) при расчете поправок
И.1.12 Уточненная с использованием очередного ИРВ 27 декабря зависимость (И.8), которой соответствует прямая 4 на рисунке И.2, существенно отклоняется от прямых 1, 2, 3. Вероятнее всего, это связано с установлением более устойчивого ледостава. В этом случае дальнейший расчет оперативных расходов воды следует осуществлять по алгоритму, описанному в 6.5 без привлечения данных о температуре воздуха.
И.1.13 Для оценки эффективности предложенного алгоритма рассчитанные оперативные расходы воды сравнивались с режимными данными, полученными по годовой КР 2005 года с учетом изменения Kзим (графа 11 таблицы И.2). Относительные погрешности оперативного учета стока на каждый день даны в графе 12. Как показывают результаты расчета, большая часть погрешностей в большинстве случаев не превышают 10%. Среднеквадратическая погрешность оперативного учета стока за весь рассмотренный период составила 8,6%.
И.2 Оперативный учет стока в зимне-весенний переходный период
И.2.1 Расчет оперативных значений расходов воды в зимне-весенний период показан на примере ГП г. Хабаровск на реке Амур для 2012 года. Данные наблюдений в зимне-весенний период предыдущего 2011 года представлены в графах 1, 2, 3, 4, 7 таблицы И.3. По итогам года была построена годовая КР, аппроксимированная аналитической зависимостью:
Q = 482,62H2 - 3653H + 10502. (И.9)
По этому уравнению для каждого ежедневного уровня рассчитываются ЕРВ, записанные в графу 5 таблицы И.3. По формуле (9) рассчитываются относительные отклонения
Таблица И.3
Данные наблюдений на ГП р. Амур - г. Хабаровск
за зимне-весенний период 2011 года и расчеты, необходимые
для построения зависимости (33)
И.2.2 Для построения зависимости
используется значение относительного отклонения и . Значение Полученная зависимость представлена на рисунке И.3.
![]() от нарастающей суммы положительных температур воздуха для ГП
р. Амур - г. Хабаровск в зимне-весенний период 2011 года
Зависимость, представленная на рисунке И.3, описывается уравнением
. (И.10)Эта зависимость используется для расчетов оперативных значений расходов воды в зимне-весенний период 2012 года.
И.2.3 Дальнейший алгоритм расчета аналогичен представленному в предыдущем примере для ГП р. Алей - с. Локоть, за исключением расчета нарастающей суммы положительных температур воздуха вместо нарастающей суммы отрицательных температур в осенне-зимний переходный период.
Результаты расчетов заносятся в таблицу И.4.
Таблица И.4
р. Амур - г. Хабаровск в зимне-весенний период 2012 года
И.2.4 Поэтапный алгоритм расчета, как если бы он производился в режиме реального времени представлен в перечислении а) - ж):
а) по уравнению КР за предыдущий (2011) год после получения среднесуточных значений уровней воды рассчитываем ЕРВ, которое помещается в графу 5 таблицы И.4;
б) после каждого ИРВ вычисляется значение относительного отклонения
; (И.11)в) с момента устойчивого перехода температуры воздуха через ноль начинаем рассчитывать нарастающие суммы их среднесуточных значений (в нашем случае это 14 апреля). При возврате устойчивых отрицательных температур расчет начинают заново. До появления устойчивых положительных среднесуточных температур воздуха оперативные значения расходов воды рассчитываются по алгоритму, изложенному в 5.2;
г) расчет поправок
д) со второго дня наличия положительных температур воздуха (15 апреля) оперативные значения оперативных расходов воды рассчитываются по формуле (8) с использованием
. (И.12)Расчет по этому алгоритму выполняется до следующего измерения расхода воды (в нашем случае это 3 мая);
е) с 30 апреля русло полностью освободилось от ледяных образований, и расчет оперативных значений расходов воды с этой даты может опять производиться по алгоритму, изложенному в 6.4, а именно, непосредственно по опорной КР, в качестве которой принята КР 2011 года;
ж) в графе 11 таблицы И.4 приведены относительные погрешности оперативного учета стока, полученные путем сравнения оперативных значений расходов воды с режимными данными, полученными в конце года. В большинстве случаев погрешности для периода весеннего таяния не превышают 15%, а среднеквадратическая погрешность за весь переходный зимне-весенний период с 14 по 29 апреля оказалась равной 13,5%.
(рекомендуемое)
В ОСЕННЕ-ЗИМНИЙ ПЕРЕХОДНЫЙ ПЕРИОД, ОСНОВАННЫЙ
НА ЗАКОНОМЕРНОСТЯХ СПАДА ЗИМНИХ РАСХОДОВ И ИСТОЩЕНИЯ СТОКА
К.1 Расчет оперативных расходов в период осеннего ледообразования с учетом закономерностей спада зимних расходов и истощения стока выполнен по данным ГП г. Алейск на реке Алей (Алтайский ЦГМС) в осенне-зимний переходный период 2005 года. Расчет оперативных расходов воды с учетом закономерностей спада зимних расходов выполнен в период от даты появления первых ледяных образований на участке ГП до даты первого измерения расхода воды в период устойчивого ледостава.
Первые ледяные образования в 2005 году на участке ГП г. Алейск на реке Алей отмечены 14 ноября, устойчивый ледостав установился 20 ноября, а первый расход при устойчивом ледоставе измерен 9 декабря.
К.2 Перед началом расчетов по данным предыдущего (2004) года по формуле (34) необходимо оценить значение параметра
Рассчитанное по (К.1) значение
К.3 Для выполнения расчетов оперативных расходов воды в осенне-зимний переходный период 2005 года формируется таблица К.1.
Таблица К.1
в период осенне-зимнего ледообразования 2005 года для ГП
р. Алей - г. Алейск, осуществленный с учетом закономерностей
спада зимних расходов и истощения стока
К.3.1 В графы 1 - 4 таблицы К.1 заносятся данные наблюдений на ГП р. Алей - г. Алейск за осенне-зимний период 2005 года.
К.3.2 В графу 5 таблицы К.1 заносится расход воды, полученный по опорной КР, в качестве которой принята КР 2004 года, аппроксимированная уравнением (К.2)
К.3.3 В графу 6 таблицы К.1 заносятся значения относительных отклонений, полученные по формуле (9).
К.4 В соответствии 6.6.2.3 для расчета оперативного расхода воды по формуле (34) используется значение расхода воды Q0 в последний день открытого русла перед началом ледообразования - в 2005 году это 13 ноября. Для оценки Q0 не может быть использован последний измеренный при свободном состоянии русла, так как он был измерен 17 октября (почти за месяц до появления ледяных образований). Поэтому Q0 определяется по опорной КР для значения уровня воды, наблюденного в день, предшествующий дате наступления ледовых явлений - 13 ноября. В графу 7 таблицы К.1 заносятся значения количества суток, прошедших от первого появления ледяных образований, а в графу 8 таблицы К.1 заносится оперативный расход, вычисленный по формуле (34). Например, на дату 20 ноября расчет оперативного расхода воды записывается в виде
. (К.3)Такой порядок расчета оперативных расходов воды сохраняется до даты первого измерения расхода воды при устойчивом ледоставе, которое в 2005 году состоялось 9 декабря.
К.5 Вычисление оперативных расходов воды в период устойчивого ледостава следует производить по алгоритму, изложенному в 6.5 для оперативного учета стока при отсутствии однозначной КР.
К.6 По окончании года для оценки эффективности принятого алгоритма оперативного учета стока необходимо выполнить сравнение оперативных и режимных расходов воды. Вычисленные ежедневные относительные погрешности оперативного учета стока заносятся в графу 10 таблицы К.1. Среднеквадратическая погрешность оперативного учета стока с учетом закономерностей спада зимних расходов и истощения стока за расчетный период составила 4,1%.
(рекомендуемое)
Л.1 Оперативный учет стока в условиях зарастания русла с использованием интерполяционно-гидравлической модели
Л.1.1 Алгоритм оперативного учета стока в период зарастания русла на основе интерполяционно-гидравлической модели заключается в следующем:
- с каждым новым измерением расхода воды, площади поперечного сечения и ширины русла вычисляется значение параметра Великанова m3 по формуле (36);
- полученное значение m3 экстраполируется - принимается для расчета оперативного расхода воды Qоп по формуле (37) вплоть до следующей даты измерения расхода воды.
Л.1.2 Для надежной оценки параметра Великанова по данным прошлых лет строятся зависимости площади живого сечения и ширины русла от уровня воды при заросшем русле
и , представленные на рисунках Л.1 и Л.2. Для построения этих зависимостей необходимо выбирать значения ![]() русла от уровня воды за период зарастания русла
2004 - 2011 годы в створе N 2 ГП р. Матыра - с. Крутое
![]() от уровня воды за период зарастания русла 2004 - 2011 годы
в створе N 2 ГП р. Матыра - с. Крутое
Л.1.3 Значения оперативных расходов воды определяются на каждую новую дату по формуле (37). При этом значения площади живого сечения и ширины русла определяются по зависимостям
и (рисунки Л.1 и Л.2) для значения среднего суточного уровня воды.Пример расчета на дату 10 июня представлен в перечислении а) - в):
а) последнее измерение расхода воды, площади живого сечения и ширины русла состоялось 31 мая. Измеренные значения составили: Qи 31 мая = 5,57 м3/с,
, Bи 31 мая = 18,3 м. По этим данным рассчитано значение параметра Великанова: ; (Л.1)б) значение среднего суточного уровня воды 10 июня H10 июн = 1,69 м. По уравнениям зависимостей, представленным на рисунках Л.1 и Л.2, определяем значения площади живого сечения и ширины русла:
;в) по формуле (37) вычисляем значение оперативного расхода воды на дату 10 июня:
. (Л.2)Л.2 Пример оперативного учета стока в условиях зарастания с использованием усовершенствованного метода оптимальной экстраполяции относительных отклонений
Л.2.1 Рассмотрим пример расчета оперативных расходов воды на р. Матыра в створе с. Крутое в период зарастания русла 2012 года.
В качестве опорной КР для расчета оперативных расходов воды принята МКР, полученная по данным расходов воды, измеренных при свободном состоянии русла в 1994, 2000 - 2011 годы, представленная на рисунке Л.3. Параметры полученной МКР приведены в таблице Л.1.
![]() при свободном состоянии русла
Таблица Л.1
р. Матыра - с. Крутое за периоды свободного русла
1994, 2000 - 2011 годы
Л.2.2 Хронологический ход относительных отклонений
![]() от опорной МКР расходов воды, измеренных в периоды
зарастания русла 1994, 2000 - 2011 годов
Значения относительных отклонений начинают уменьшаться после перехода среднесуточной температуры воды через 8 - 10 °C незадолго до явного появления водной растительности. По ходу развития растительности
В рассматриваемом примере аппроксимацию хронологического хода относительных отклонений
Сложнее определить даты окончания нарастания и начала отмирания. С этой целью необходимо построить хронологические графики отклонений
![]() отклонений
(период зарастания с 4 мая по 29 ноября. 2010 г.)
Л.2.3 Значения параметра
, (Л.3)где n - последний по счету день нарастания/спада растительности.
Длительности периодов нарастания и спада вычисляются как разности дат окончания и начала соответствующих периодов.
Л.2.4 Значение параметра a вычисляется для периода нарастания как среднее из максимальных (по модулю) значений
, (Л.4) , (Л.5)где n - число лет в совокупности.
Л.2.5 С использованием полученных параметров
- для периода нарастания водной растительности:
- для периода спада (отмирания) водной растительности:
Л.2.6 Уравнения для расчета АКФ приобретают следующий вид:
- для периода нарастания водной растительности
- для периода спада (отмирания) водной растительности:
На рисунке Л.6 представлены: 1 - АКФ, рассчитанная по формуле (Л.8) для периода нарастания, 2 - АКФ, рассчитанная по формуле (Л.9) для периода спада водной растительности. Как видно из рисунка Л.6, в период нарастания связанность между значениями ряда
![]() 1 - для периода нарастания водной растительности;
2 - для периода спада водной растительности
Л.2.7 Следующие шаги соответствуют расчетам по методу оптимальной экстраполяции, изложенному в приложении Ж.
Л.2.8 Представленные на рисунке Л.4 отрезки косинусоид 2 и 3 отражают осредненные за многолетние характеристики ПСР в створе ГП р. Матыра - с. Крутое в период его зарастания водной растительностью. Как видно из графика (см. рисунок Л.4), разброс точек
Значения поправок
- для периода нарастания водной растительности:
- для периода спада (отмирания) водной растительности:
где
Исследования показали, что в период приблизительно стабильного состояния водной растительности, когда не происходит значительного изменения ПСР, поправки
Л.2.9 Процедура расчета оперативных расходов воды в период зарастания русла 2012 г. для ГП р. Матыра - с. Крутое, результаты которого представлены в таблице Л.2, представлена в перечислении а) - ж):
а) после завершения периода ледовых явлений на каждые сутки по значению среднего суточного уровня воды вычисляется значение расхода по опорной МКР. В то же время необходимо отслеживать рост температур воздуха и воды;
б) с каждым новым измерением расхода воды по формуле (9) получаем значение отклонения
в) если при приближении средней многолетней даты начала периода зарастания и установлении устойчивых температур воды выше 10 °C - 12 °C значение
г) расчет значений АКФ осуществляется по формуле (Л.8). Расчет весовых коэффициентов p1 и p2 - по формулам (19) и (20);
д) с каждым новым измерением расхода воды производится расчет значений остаточных отклонений
е) экстраполированные значения
ж) расходы воды вычисляются по формуле (8).
Таблица Л.2
р. Матыра в створе с. Крутое в период зарастания русла
Л.2.10 Например, на дату 3 июня, 23-й день от начала зарастания, расчет выглядит следующим образом:
- предварительно по данным прошлых лет рассчитывается значение меры погрешности измерений расходов воды
. (Л.13)- по значению уровня воды, наблюденного 3 июня (H = 1,52 м) вычисляется расход воды (Q(H) = 11,0 м3/с);
- по значениям относительных отклонений последнего и предпоследнего ИРВ, равным соответственно:
, и , и значениям отклонений и , по формуле (Л.10) вычисляются остаточные отклонения ; (Л.14) ; (Л.15)- с учетом значений
, по формуле (Л.8) вычисляются значения АКФ: r01 = 0,993, r12 = 0,908, r02 = 0,853;- затем рассчитываются весовые коэффициенты формулы оптимальной экстраполяции:
; (Л.16) . (Л.17)- по формуле (Л.6) рассчитывается значение отклонения
;- далее вычисляется экстраполированное значение поправки по формуле (Л.11)
(Л.18)и значение оперативного расхода воды на дату 3 июня по формуле (8)
Qоп = 11,0·(1 - 0,518) = 5,31 м3/с. (Л.19)
Дата начала периода приблизительно стабильного состояния растительности определяется как средняя по многолетним данным. Для этого периода значения поправок
Дату начала периода спада развития водной растительности также можно определить как среднюю многолетнюю. Для периода спада расчет оперативных расходов воды выполняется аналогично описанному в перечислении а) - ж). При этом расчет значений АКФ осуществляется по формуле (Л.9), расчет отклонений
Л.2.11 На рисунке Л.7 представлены режимный и оперативный гидрографы стока реки Матыра в створе с. Крутое за период зарастания русла 2012 г.
![]() 1 - режимный; 2 - оперативный; 3 - ИРВ
в створе села Крутое за период зарастания русла
с 18 мая по 10 ноября 2012 г.
Л.2.12 Для оценки эффективности предложенных методов выполнен расчет оперативных расходов воды в периоды зарастания русла 2008 - 2013 годов. Результаты сравнительных расчетов оперативных расходов воды с использованием экстраполяционно-гидравлической модели и усовершенствованного метода оптимальной экстраполяции приведены в таблице Л.3.
Таблица Л.3
в период зарастания русла
Относительные погрешности оперативных ЕРВ определялись по формуле
, (Л.20)где Qр t - значение ЕРВ по данным, опубликованным в справочниках водного кадастра "Ежедневные данные о режиме и ресурсах поверхностных вод суши".
Случайная среднеквадратическая погрешность оперативного учета за весь период зарастания русла рассчитывается по формуле
, (Л.21)где N - количество членов ряда.
Систематические погрешности оперативного учета стока определялись как средние значения по всей совокупности относительных погрешностей оперативного учета стока в период зарастания.
Значение максимально допустимой относительной погрешности оперативного учета стока принято равным 15%.
Как видно из таблицы Л.3, расчет оперативных расходов воды в период зарастания русла с использованием усовершенствованного метода оптимальной экстраполяции относительных отклонений позволяет существенно снизить случайные средние квадратические, систематические и наибольшие погрешности оперативного учета стока по сравнению с экстраполяционно-гидравлической моделью.
Следует отметить, что отдельные случаи превышения допустимой погрешности наблюдаются в каждом из исследованных периодов зарастания. Установлено, что в период приблизительно стабильного состояния водной растительности снижение ПСР, которому соответствует увеличение отрицательных значений
Лист регистрации изменений
Вернуться в "Каталог нормативных документов"
Источник информации: https://internet-law.ru/documents/prod/prikaz/110/r_52.html
На правах рекламы:
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||