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

Кривые распределения и кривые обеспеченности стока. Обеспеченностью заданного значения гидрологических характеристик называется вероятность ежегодного его превышения.

ЛЕКЦИЯ N8

Содержание: Изменчивость характеристики стока. Кривые распределения и кривые обеспеченности стока. Параметры кривых. Применение эмпирических и аналитических кривых (трехпараметрического гамма распределения и биноминальной кривой) к расчетам стока.

1.Изменчивость характеристики стока. Величины годового стока непрерывно колеблются под влиянием климатических (осадков и испарения) и физико-географических факторов. Длительными наблюдениями установлено, что колебания годового стока носят циклический характер, выражающийся в последовательной смене многоводных и маловодных лет. Установлено, что появление в том или ином календарном году того или иного годового стока из числа возможных является случайным. Также являются случайными календарные последовательности значений годового стока. Кроме того, как показали исследования, между величинами стока за смежные годы имеется лишь незначительная корреляционная связь, что дает возможность в большинстве случаев пренебрегать ею. Следовательно величины годового стока являются случайными, а ряд образованный ими, вариационным рядом.

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

В инженерной гидрологии для изучения колебаний годового стока и оценки изменений его характеристики в будущем широко используются кривые распределения вероятностей и обеспеченности стока.

Кривая обеспеченности (или кривая превышения) - это интегральная кривая, показывающая обеспеченность или вероятность превышения (в процентах или в долях от единицы) данной величины среди общей совокупности ряда. Следовательно, если обеспеченность годового стока Q_i (или W_i) равна Р%, это означает, что число членов в данном ряду, больших или равных Q_i, составляет Р% от общего числа членов ряда.

Рассмотрим как можно построить кривую повторяемости (частоты) и эмпирическую кривую обеспеченности среднегодовых величин стока.

Пусть имеются данные наблюдений за годовым стоком за период n лет: Q₁,Q₂,...,Q_i,..., Q _n (где Q_{i -} средний годовой сток в i-м году). Выразим эти данные в относительных величинах (модульные характеристики стока): где Q_{0 -} среднее арифметическое значение ряда) и расположим их в убывающем порядке от Кmах до Кmin. Если такой ряд разбить на одинаковые интервалы DК и подсчитать абсолютную n_j частоты повторений каждого интервала модульных коэффициентов Кj (), то можно построить ступенчатый график распределения частоты, или вероятностей, среднегодовых расходов воды. Такой график в математической статистике называется гистограммой распределения (рис. 1).

Рис. 1 Гистограмма (1) и эмпирическая кривая обеспеченности (2) среднегодовых расходов воды.

При неограниченном росте членов ряда (n ® ¥) и стремлении DК ® 0 гистограмма распределения превращается в плавную кривую, которую называют кривой распределения вероятностей (рис. 2). Кривая распределения дает наглядное представление о законе распределения случайной величины. Она характеризует вероятность появления того или иного значения рассматриваемого ряда случайных величин. Из ее рассмотрения, в частности, следует, что чем ближе значение члена ряда к среднему, тем больше вероятность его появления (частота) и, наоборот, чем больше отклонение члена ряда от среднего, тем меньше его вероятность появления.

Рис. 2 Кривая распределения вероятностей.

Кривая распределения имеет три характерные точки на оси абсцисс (рис. 2).

Точка 1 - центр распределения, соответствует среднеарифметическому значению ряда (ордината проходящая через эту точку, называется центральной); точка 2 - медиана, делит ряд на две равные части; точка 3 - мода, соответствует члену ряда, имеющему наибольшую частоту (модальная ордината).

Если центральная, медианная и модальная ординаты совпадают и образуют ось симметрии, кривую распределения называют симметричной. Наоборот, в асимметричных кривых эти ординаты не совпадают. Расстояние d между центральной и модальной ординатами, называемое радиусом асимметрии. характеризует степень асимметричности кривой.

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

Параметры кривых. Параметрами аналитических кривых распределения и кривых обеспеченности являются: средний годовой сток (Q₀), среднеквадратическое отклонение (s _Q), изменчивости (вариации) (C_v)_, коэффициент асимметрии (C_s) и коэффициент автокорреляции.

Методы вычисления C_v и C_s з ависят от коэффициента изменчивости: а) при C _v< 0.5 применяется метод моментов; б) при С _v > 0.5 используется метод наибольшего правдоподобия, в) графо-аналитический метод применяется в случае использования биноминальной кривой обеспеченности при любом значении С_v.

Метод моментов заключается в том, что искомые параметры распределения выражаются через статистические моменты ряда наблюдений. Приближенные параметры аналитических кривых распределения находят по формулам (1), (2),(3), (9) приведенные в лекции № 7.

Метод наибольшого правдоподобия состоит в том, что в качестве приближенных значений параметров распределения принимают такие значения, при которых результаты наблюдений имеют наибольшую вероятность совместного появления. Этот метод доведен до стадии практического использования только применительно к трехпараметрическому гамма-распределению. Составлены номограммы позволяющие находить параметр распределения. Таким образом, для определения параметров распределения C_v и C_s методом наибольшого правдоподобия вычисляют статистики l₁, l₂ и l₃ по формулам:

(1)

(2)

(3)

Как видно из уравнения (1) статистика l₁ полностью совпадает со средним арифметическим значением x₀ в методе моментов. Коэффициенты изменчивости C_v и ассиметрии C_s устанавливаются, как уже отметили, по специально составленным для этой цели номограммам.

Применение эмпирических и аналитических кривых к расчетам стока. Для построения эмпирической кривой обеспеченности необходимо величины годового стока расположить в убывающем порядке от наибольшего значения к наименьшему. Такой ряд дает общее представление о колебании значений годового стока за период наблюдений.

Обеспеченность каждого члена ряда, расположенного в убывающем порядке, выраженная в процентах (или в долях единицы), может быть вычислена по формуле

(4)

где Р - обеспеченность данного расхода в процентах, m -порядковый номер данного расхода в ряду, n- общее число членов ряда.

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

(5)

Эмпирическая кривая обеспеченности может быть построена в прямоугольных координатах или специальной сетке так называемой клетчатке вероятностей. По оси ординат откладываются значения годовых расходов (или модульных коэффициентов), а по оси абсцисс - значения обеспеченности (в процентах) этих расходов, вычисленные по формулам (4) или (5). Нанесенные точки соединяются плавной кривой.

Как видно из формул (4) и (5), обеспеченность крайних точек является наименее достоверной и наименее устойчивой из всех точек. И зависит от длительности ряда наблюденй. Предпочтительнее строить кривую обеспеченности на клетчатке вероятностей. При этом кривые распрямляются.

Эмпирическая кривая обеспеченности позволяет определять расчетные значения гидрологической характеристики Q_р, заданной вероятности превышения Р.

(6)

где К_р - ордината кривой обеспеченности, соответствующая заданному значению Р%; Q₀ - среднее арифметическое значение (норма) гидрологической характеристики.

Зная эмпирическую обеспеченность гидрологической характеристики, можно подсчитать вероятность повторяемости в годах, т.е. число лет N, в течение которых данное значение характеристики встречается в среднем один раз. Обеспеченность Р и повторяемость N связаны между собой следующим образом.

при Р< 50% (7)

при Р> 50% (8)

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

При определении расчетных гидрологических характеристик, широко применяются различные теоретические функции распределения. Среди многочисленных кривых распределения наибольший практический интерес представляют функции гамма распределения, распределение Пирсона III типа и трехпараметрическое гамма-распределение С.Н. Крицкого и М.Ф. Менкеля (1946 г.).

Гамма распределение занимает особое положение при описании существенно положительных случайных величин, варьирующих в интервале

¥ >Х >0. Величины речного стока не могут принимать значений меньше нуля, поэтому гамма распределение широко применяется в гидрологии.

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

Уравнение плотности вероятностей гамма-распределения имеет вид:

(9)

где Х₀ - центр (среднее) распределения; g - параметр, связанный с коэффициентом изменчивости простым соотношением:

(10)

Г(g) - cимвол гамма-функции (интеграла Эйлера II рода);

Г(g)

(11)

Гамма распределение содержит два параметра, характеризующих центр распределения и его рассеяние. Это распределение имеет положительную асимметрию, коэффициент асимметрии С_s функциональной связи с коэффициентом изменчивости С_v соотношением С_{s =} 2 С_v.

Распределение Пирсона III типа. Распределение вероятностей многих гидрологических величин (ливневые осадки, максимумы дождевых паводков в малых реках и т.п.) характеризуются значительной асимметрией, иногда в несколько раз превышающий асимметрию гамма-распределения. В связи с этим возникает необходимость в более гибком - трехпараметрическом распределении вероятностей. До сравнительно недавнего времени в качестве такого распределения применялись почти исключительно кривые Пирсона III типа.

Уравнение плотности вероятностей этих кривых имеет вид

при х>

(12)

где x₀ -среднее распределение; и - параметры распределения.

Связь между этими и стандартными параметрами и выражается соотношениями:

(13)

Cемейство кривых распределения Пирсона III типа допускает свободный выбор трех параметров: среднего x₀, коэффициента изменчивости и коэффициента асимметрии .

Распределение Пирсона является линейным преобразованием гамма распределения: она получается "сдвигом" гамма распределения по оси абсцисс на величину . Для различных значений , мы будем получать распределения с различной относительной асимметрией /. Величины , т.е. начало кривой, связаны с и равенством:

(14)

Таким образом, как видно из уравнений (14) распределение Пирсона III типа ограничено снизу величиной , которая определяет асимметрию распределения. Последнее формально следует из уравнений кривой Пирсона III типа и никак не определяется природой явления речного стока. Это является принципиальным недостатком кривых Пирсона III типа. Кроме того, как видно из рис.2, следует, что

(15)

C другой стороны, при интегрировании биномиальной кривой распределения получено:

(16)

Cоответственно, из (15), (16) получаем:

(17)

Анализ (17) показывает:

если = 2 , то К_min= 0 и, следовательно, кривая выходит из начала координат;

при > 2 , то К_min> 0 и начало кривой отстоит вправо от начала координат на величину К_min;

если же < 2 , то К_min< 0, т.е. имеет отрицательное значение, что противоречит физической сущности гидрологической характеристики.

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

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

Уравнение плотности вероятностей трехпараметрического распределения С.Н.Крицкого и М.Ф.Менкеля имеет вид:

(18)

В этом уравнении x₀ -среднее распределение; и b - параметры распределения, каждому сочетанию которых соответствуют определенные значения коэффициентов изменчивости и коэффициентом асимметрии . Связь между параметрами распределения и b и стандартными параметрами и выражается следующими трансцендентными уравнениями:

(19)

(20)

Аналитические кривые обеспеченности при известных параметрах Q₀, C_v и C_s строят с помощью таблиц, приведенные в работе "Практикум по гидрологии, гидрометрии и регулированию стока. М.; Агропромиздат,1998"(Приложение 2, с.215). В данной работе представлены результаты интегрирования соответствующих кривых распределения. Ординаты кривой обеспеченности, представленные в этих таблицах даны в виде модульных коэффициентов для значений обеспеченности P% от 0,001 до 99,9 % при различных коэффициентах вариации C_v - от 0,1 до 1,0.

Ординаты аналитических кривых обеспеченности определяют следующим образом:

1. Выбирают таблицу соответствующую определенному соотношению C_s/C_v.

2. По выбранной таблице для заданного коэффициента вариации C_v выписывают значения K_P при различных значениях P%. Если заданный коэффициент вариации C_v не совпадает с табличным, то проводят интерполирование.

3. Определяют абсолютные ординаты кривой обеспеченности

Q_P=K_PQ₀

По вычисленным ординатам строят график аналитической кривой обеспеченности трехпараметрического гамма распределения.

Дата добавления: 2014-01-14