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

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

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

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

Добавлен: 15.10.2024

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

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

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

Аналогичным образом, вероятность того, что конкретное значе­

ние Xi принадлежит к совокупности

(х ), равна

П2

(•*;)•

Щ+ П2

Наконец, определим вероятность появления конкретного значе­ ния Xi в совокупности Pi(x) или Рг(л:); она в соответствии с теоре­

мой сложения вероятностей равна

Р(хд

щ + л2 Л СО

«2

^2 СО,

(4.14)

Щ+ П2

так как любое конкретное значение лц принадлежит либо к совокуп­ ности Р1 (х), либо к совокупности Рг (х).

Обобщая приведенное доказательство на случай k неоднород­

ных совокупностей, получаем приведенное выше равенство (4.13), которое Бровковичем и Великановым [30] использовалось для оценки соответствия уравнений кривых распределения эмпириче­ ским данным речного стока.

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

В § 11 главы II при изложении типовых эмпирических кривых распределения рассматривались особенности формирования стока весеннего половодья, приводящие к неоднородности рядов стока ве­ сеннего половодья в некоторых районах. Наличие такой неоднород­ ности приводило к своеобразной форме эмпирических кривых рас­ пределения. Технику построения составных аналитических кривых распределения для описания неоднородных эмпирических совокуп­ ностей рассмотрим на примерах расчета колебаний речного стока.

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

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

На рис. 4.2 нанесены эмпирические точки, соответствующие рас­ сматриваемым совокупностям, и аналитические кривые обеспечен­ ности. Построение кривых 3 и 4 выполнено обычными приемами,

исходя из указанных на чертеже параметров.

Рассмотрим порядок построения кривой 5. Исходя из общего объема совокупности (35 членов ряда) и однородных совокупно­ стей (20 и 15 расходов), получаем доли (веса) каждой совокуп­ ности:

А = - § - = 0 ,5 7 , ^2= 4 J -= 0 ,4 3 .

210


Используя кривые обеспеченности (3 и 4), построенные для од­

нородных совокупностей, производим расчет аналитической неод­ нородной кривой обеспеченности по схеме, приведенной в табл. 4.4.

На рис. 4.3 представлена эмпирическая кривая обеспеченности годового стока р. Сакмары у с. Сакмары (эмпирические точки всего ряда) и ее аналитическая аппроксимация в форме интегральной кривой Крицкого—Менкеля при CS = CV (кривая IV). Параметры

этой кривой получены с использованием всей выборки.

Qr 3/C 1000 200100 20 10 5

2

5 10 20 100 200 1000

4

3

2

1

О

0,01 0,1 1 5 10 20 40 00 80 90 95 99 99,9 Р %

Рис. 4:2. Кривая обеспеченности максимальных расходов воды р. Абава у х. Сисени.

/ — эм п и р и ч е ск и е то ч к и , со о т ве тст ву ю щ и е

су м м а р н о й

к р и во й об есп еч ен н о сти

м а к ­

с и м а л ьн ы х

р ас х о д о в

в о д ы ; 2 — э м п и р и ч е ск и е то ч ки ,

со о т ве тст ву ю щ и е кри вой

о б е с ­

п еч ен н ости

сн его вы х

м ак си м у м о в (в е р х н я я

к р и в а я )

и

кр и во й об есп еч ен н о сти

д о ж ­

д ев ы х м ак си м у м о в (н и ж н я я к р и в а я ); 3 — а н а л и т и ч е с к а я к р и в а я К р и ц к о го —М ен к е л я ,

о т в е ч а ю щ а я

со во к у п н о сти м а к с и м а л ь н ы х

р а с х о д о в

в о д ы

в есен н его п о л о в о д ь я при

Q=*184 м 3/с ,

C v =*0,36,

Ca = 2 C v ;

4 — а н а л и т и ч е с к а я

к р и в а я ,

о т в е ч а ю щ а я с о в о к у п н о ­

сти д о ж д е в ы х м а к с и м а л ь н ы х

р а с х о д о в

воды при

Q*=128

м 3/с, (7^ = 0,52, CS= 2 C V \

б — а н а л и т и ч е с к а я к р и в а я о б есп еч е н н о ст и , р а с с ч и т а н н а я

н а

о сн о ван и и ко м п о зи ц и и

с и сп о л ьзо ва н и ем

к р и вы х 3

и

4, о т н о с я щ и х с я к

о д н о р о д н ы м со в о к у п н о стя м .

 

П о

оси о р д и н а т

р а с х о д ы у м ен ьш ен ы в

100 р а з .

Отметим, что принятая аналитическая кривая обеспеченности наилучшим образом осредняет рассматриваемое эмпирическое поле точек по сравнению с другими аналитическими кривыми и различ­ ными соотношениями Cs/Cv. Больше того, интегральная кривая Крицкого—Менкеля при CS=CV испытана на большом материале

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

1 4 *

211


Рис. 4.3. Кривая обеспеченности го­ дового стока р. Сакмары у с. Сакмары,

/ — эм п и р и ч е ск и е т о ч ­

ки о тн о си тел ьн о о д ­

н ородны х с о в о к у п н о ­

стей,

2 — э м п и р и ч е ­

ские точ ки в сего р я д а . К ривы е о б ес п еч е н н о ­

сти

К р и ц к о го —М ен-

к е ля

при

CS= C V: I

п ер в ая

к о м п о н ен та :

лГ=3,64,

С и = 0,46,

п =

= 68, CS= C V\ I I — в т о ­

р а я к о м п о н ен та :

М =

=8,98,

С в - 0 ,1 1 ,

п = 12,

CS = CV;

I I I

с у м ­

м а р н а я и н т е гр а л ь н а я

к р

и вая ,

о с н о в а н н а я

на

су м м е

взвеш е н н ы х

вер о ятн о стей

I к I I

ком пон ент;

I V

— с у м ­

м а р н а я и н т е гр а л ь н а я к р и в а я по в сем у р я д у

н аб лю д ен и й : A f=4.45, C v = 0,56, я = 80, Cs =

Т а б л и ц а 4.4

Схема расчета аналитической неоднородной кривой обеспеченности максимальных расходов воды р. Абавы у х. Сисени

 

С овокупность дож девы х

С овокупность снеговых

А бсцисса неодно­

 

максим ум ов

 

 

максим ум ов

^ м а к с

 

 

родной

кривой

 

 

 

 

 

 

м 3/с

 

 

 

 

 

 

%)

р, %

0,43 Р ,

%

 

 

0,57р 2 %

Р„ = 0 ,4 3 Р ,-|-

 

Рг

%

 

+ 0,57 Р 2

490

0,016

0,007

0,04

0,023

0,03

400

0 ,2 1

0,09

 

0,48

0,27

0,36

350

0,65

0,28

 

1,7

0,97

1,25

300

1,98

0,85

 

5,2

2,97

3,82

250

5,3

2,27

 

14,3

8,17

10,4

2 0 0

13,8

5,92

 

35,8

20,4

26,3

150

31,9

13,7

 

6 6 , 0

37,7

51,4

1 0 0

61,0

26,2

 

91,5

52,2

78,4

50

91,5

39,3

 

99,75

56,96

96,3

40

95,3

40,9

 

99,96

57,08

97,7

10

99,96

42,88

 

1 0 0

 

57,0

99,98

как легко

обнаружить

по

рис.

4.3,

интегральная

кривая

распределения, полученная

по всей

совокупности,

на некоторых

участках существенно отличается от расположения

эмпирических

точек, отклоняясь от них вниз при больших обеспеченностях и вверх в экстраполируемой верхней части этой кривой. Указанное несоот­ ветствие не может быть объяснено только «случайными» отклоне­ ниями эмпирических точек от аналитической кривой. Больше того, аналогичное несоответствие эмпирических данных по стоку анали­ тическим интегральным кривым распределения наблюдалось и для других рек рассматриваемой территории.

Итак, исходя из физического анализа и учитывая систематиче­ ское уклонение эмпирических точек от теоретических кривых, можно сделать предположение, что этот ряд нельзя рассматривать как вполне однородный. Видимо, он состоит из нескольких относи­ тельно однородных совокупностей.

Первые выводы можно сформулировать следующим образом:

1)данный статистический ряд в первом приближении можно разбить на две относительно однородные совокупности;

2)начиная с модулей стока приблизительно 7—8 л/с-км2 и больше действует какая-то причина, нарушающая однородность-

ряда наблюдений; 3) действие этой причины проявляется в увеличении модулей

стока многоводных лет.

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

Признавая неоднородность ряда годового стока р. Сакмары,. произведем разделение разнородного распределения на относи-

21S


тельно однородные совокупности с учетом физического анализа формирования стока и с учетом расположения точек на эмпириче­ ской кривой обеспеченности.

На рис. 4.3 представлены эмпирические точки и интегральные кривые распределения Крицкого и Менкеля при CS=CV относи­

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

Суммарную кривую обеспеченности можно получить по двум относительно однородным совокупностям, включающим соответст­ венно 68 и 12 членов. В таком случае первая совокупность будет участвовать в расчете с весом 0,85, вторая — с весом 0,15.

При выполнении расчета берем несколько произвольных точек на оси ординат, например, соответствующих следующим значениям годовых модулей стока: 12, 10, 8, 6, 4, 2 л/с • км2.

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

Т а б л и ц а 4.5

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

М одуль

1-я совокупность

2-я совокупность

С ум м арное

р ас п р е ­

стока,

 

 

 

 

делен ие

( в % )

р, %

 

 

 

P = 0,85P ,-f-

л / с ‘ К М 2

0.85Р , %

 

0,15 Р 2 %

Рг %

—j—0,15/^2

12

0.01

0,008

0,18

0,027

0,033

10

0.04

0,034

14,5

2,18

2,21

8

1.00

0,85

83,5

12,53

13,38

6

9,3

7,9

99,91

14,87

22,77

4

36,0

30,6

99,99

15,0

45,6

 

2

81.5

69,3

99,99

15,0

84,3

 

На рис. 4.3 видно, что суммарная интегральная кривая обеспе­ ченности III значительно лучше осредняет эмпирическое поле то­ чек, чем аналитическая кривая обеспеченности IV, полученная по

разнородной совокупности.

Разброс точек всего эмпирического ряда относительно суммар­ ной кривой III можно отнести к случайным отклонениям. Однако

и в этом случае в диапазоне модулей 5,0—6,5 л/с • км2 наблюдается некоторое систематическое отклонение точек от суммарной кри­ вой III. В дальнейшем, разбив эту совокупность на две относи­

тельно однородные совокупности, можно было бы исключить и это незначительное несоответствие кривой обеспеченности III наблю­

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

2 1 4


ности, добиваясь более полного соответствия расположения эмпи­ рических точек вычисленной аналитической кривой.

Рассмотрим еще один пример для р. Большой Узень у г. Ново-

узенска. Отметим, что в пределах данного водосбора бессточные депрессии более развиты, чем в условиях р. Сакмары. Кроме того, на водосборе р. Большой Узень имеется большое количество мел­ ких прудов, водохранилищ и больших плёсов, которые дополни­ тельно снижают сток маловодных лет, увеличивая их количество и сводя сток исключительно маловодных лет до нуля (1933 г.).

М л/с км 2

Рис. 4.4. Кривая обеспеченности годового стока р. Большой Узень у г. Нозоузенска.

/ — э м п и р и ч е ск и е

то ч к и о тн о си тел ьн о о д н о р о д н ы х со во к у п н о сте й , 2 — эм п и р и ч е ск и е

точ ки

в сего р я д а ; / ,

I I ,

I I I

— а н а л и т и ч е с к и е

к р и вы е

о б ес п еч ен н о сти

о тн о си т ел ьн о

о д н о р о д н ы х

со во ку п н о стей ;

I V

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

к р и в а я о б есп еч ен н о сти ,

о с н о в а н н а я н а

су м м е

в з в е ­

ш ен н ы х вер о ят н о сте й

к р и вы х / ,

I I ,

/ / / ;

V — а н а л и т и ч е с к а я

к р и в а я о б есп еч ен н о сти ,

р а с ­

 

 

с ч и т а н н а я

по

р а зн о р о д н ы м

д а н н ы м в сего

р я д а .

 

 

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

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

215