Файл: Рождественский, А. В. Статистические методы в гидрологии.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 15.10.2024

Просмотров: 185

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

В уравнениях интегральной биномиальной кривой распределе­ ния в виде (2.45) и (2.47) величина обеспеченности признака рас­

пределения

связывается с величинами коэффициентов

 

х

вариации и асимметрии.

Параметры более высоких моментов, чем третий, функционально связаны с коэффициентами асимметрии и вариации. Воспользу­ емся этим для определения эксцесса через величины параметров Cv и Cs, для чего в рекуррентное соотношение для моментов (2.29)

подставим (2.30)

 

 

 

 

 

 

2mi

-т4-{-■2тп

 

 

 

 

 

 

Зт2-\-- т2

 

Из полученного уравнения определен четвертый момент

 

 

 

 

6пё, + 4

тз

3 (2т\ -)- m^j

 

 

 

nli

 

 

2т2

 

2т2

 

Следовательно,

выражение эксцесса можно представить в виде

р _

т 4

з -

3

 

2т\ + т\ \

о

3 (2т\ + т\ 2яг®)

2>т\

С

-

1 -

 

 

 

 

„ 2

й

2

1

А

 

 

2т\

2т\

 

т2

 

 

 

)

 

Наконец, выражая моменты через

параметры Cv и

Cs, полу

чаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^ = - |- С*>

 

(2.48)

а при

CS = 2CV

 

 

 

 

 

 

 

 

 

 

 

 

 

E — 6Cl.

 

(2.49)

Аналогичным образом через Cs и С„ могут быть выражены лю­ бые моменты (т 3 ).

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

от величин параметров х, Сv и Cs. Впервые такую таблицу соста­

вил А. Фостер [151]. В СССР она была уточнена С. И. Рыбкиным [117], а затем и другими авторами.

Эта таблица в современном ее виде опубликована в «Руковод­ стве по определению расчетных гидрологических характеристик». Указанная таблица1 представляет собой численное решение урав-

1 В дальнейшем при упоминании этой таблицы указание на необходимость использования Руководства не повторяется.

90


нения (2.45) при C\,= l и различных р (выраженных в процентах) и Cs. В таблице приведены величины отклонения модульных коэф­

фициентов { ^ = ~^J от еДиницы (т-

е- среднего значения

при­

знака), выраженные в долях коэффициента вариации1 (СД,

т. е.

 

 

Xп

 

yP= f ( c s, р)-~

1

- 1

(2.50)

 

СV

Су

 

 

Следовательно, для определения ординат биномиальной кривой

обеспеченности по заданным параметрам х, Cv и Cs табличные зна­ чения величины ур, полученные при заданном Cs для различных

обеспеченностей р%, необходимо преобразовать в соответствии с равенством (2.50) к виду

kpУрСу-\-^>

(2.51)

х р= { урС„-{-\)х.

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

kP- 1

Ур Со

или

хр= у рз-\-х.

Приведем пример. Для статистической совокупности величин годового стока р. Днепра у пгт Лоцманской Каменки получены сле­

дующие значения параметров: <7= 3,78 л/с-км2; = 0,28; Cs = 0,46.

Рассчитаем модули стока различной обеспеченности (т. е. значения ординат биномиальной кривой обеспеченности) при указанных зна­ чениях параметров. Исходные величины у р при различных р для

Cs= 0,46 взяты из таблицы биномиального закона. Порядок рас­ чета дан в табл. 2.4.

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

! Табличное значение величин

kp— 1

в гидрологической литературе часто

Cv

 

kP-

 

обозначают символом <!>(*) =

1

 

 

С,

 

 


Т а б л и ц а 2.4

Порядок расчета величин годовых модулей стока различной обеспеченности (qp) р. Днепра у пгт Лоцманской Каменки

О беспеченность Р %

1

3

5

10

50

60

80

90

99

Ур УpCv

k — УpCv + О ex II '£

 

2,64

2,06

1,76

1,32

-0,075

-0 ,3 2

—0,85

-1 ,2 2

- 2 ,0

 

0,74

0,58

0,49

0,37

-0,021

-0 ,0 9

-0 ,2 4

-0 ,3 4

-0 ,5 6

I

1,74

1,58

1,49

1,37

0,98

0,91

0,76

0,66

0,44

 

6,58

5,20

5,63

5,18

3,70

3,44

2,87

2,49

1,66

•обеспеченность членов статистического ряда обычно определяется по формуле

100°/о,

(2.52)

где т = 1 , 2 , . . « — порядковые номера членов статистического ряда, расположенных в убывающем порядке; п — общее число чле­

нов статистической совокупности.

При сопоставлении аналитической и эмпирических кривых обес­ печенностей (рис. 2 .8 ) обнаруживается достаточно хорошая их

сходимость. Рядам годового стока и некоторым другим характери­

стикам гидрологического

режима свойственно соотношение Cs —

= 2Cv. Для этого частного,

но часто используемого соотношения

с целью облегчения расчетов составлена специальная таблица ор­ динат кривых обеспеченностей, выраженных в модульных коэффи­ циентах в зависимости от Cv и р \

yp=*kp= ^ £ - = f ( C B, р ).

(2.53)

Рассмотренная биномиальная кривая распределения представ­ ляет частный случай (III тип) обширной системы кривых К. Пир­ сона, состоящей из 13 типов распределения (включая нормальную кривую).

Совокупность кривых Пирсона получается в результате реше­ ния дифференциального уравнения вида

1 .

аУ

— ____ г + х____

/

/о щи

у

dx

Ь0 b\X + b-2 X2 ’

 

V • /

которое представляет собой более общий вид рассмотренного ра­ нее дифференциального уравнения биномиальной кривой распре­ деления (2.23). Уравнение (2.54) соответствует статистической1

1 См. «Руководство по определению расчетных гидрологических характери­ стик» (приложение 3, случай Cs= 2Cr ).

9 2


3,0

0,2

0,1-----------

— ----

— L J — _J-------------------------— — ______ _________ ____ _________

0,01

0,1 0,2

0,5

1 2

5

10

20 30 40 50 60 70 80

90 95

98

99 99,5 99,8P %

Рис. 2.8.

Аналитическая

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

коэффициентов

годового стока

 

 

 

 

р.

Днепра у пгт Лоцманской Каменки.

 

 

 

/ —эмпирические точки, 2—биномиальная кривая.

схеме, описываемой так называемым гипергеометрическим распре­ делением. Это распределение получается, если оценивается вероят­ ность появления белого шара т раз в S извлечениях, когда в рас­ сматриваемой совокупности объемом п имеется пр белых шаров и nq черных ( р > 0, <7 > 0, p + q=\), причем извлеченный шар снова

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

В геометрическом смысле уравнение (2.54), так же как и урав­ нение (2.23), отражает характерные особенности очертания кривых распределения, которые обычно начинаются от нуля (в отношении значений вероятности), достигают максимума и затем снижаются с той же или большей частью с иной интенсивностью.

Использование в уравнении (2.54) третьего параметра (62) при­

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

Сведения о всей совокупности кривых Пирсона можно получить из работ по статистике [111, 144].

Очертание биномиальной кривой обеспеченности при различ­ ных Cs и Cv показано на рис. 2.9.

§ 5

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

Кривая распределения вероятностей Пирсона III типа, широко используемая в практике гидрологических расчетов, обладает, как было показано, одним существенным недостатком, заключающимся в том, что она при CS<2CV уходит в область отрицательных зна­

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

Лишь в одном частном случае, при CS = 2C„, кривая Пирсона III типа охватывает область изменения переменной от нуля до бес­ конечности, что не входит в логическое противоречие с представ­ лениями о законах колебания речного стока и многих других гидро­ логических величин.

94