Файл: Рождественский, А. В. Статистические методы в гидрологии.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