Файл: Рождественский, А. В. Статистические методы в гидрологии.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 15.10.2024
Просмотров: 199
Скачиваний: 0
Из теоретического |
анализа, выполненного, в частности, |
Г. А. Алексеевым [11], |
следует, что предельное (при л-»-оо) значе |
ние третьего центрального момента (цз) |
равно 2,404, в таком случае |
||||
С |
= с , |
2,404 |
|
1,14. |
|
1,2823 |
|||||
|
|
|
Рассмотрим пример вычисления максимальных расходов воды дождевых паводков р. Москвы у г. Звенигорода за период наблю-
Рис. 2.14. Эмпирическая и аналитические кривые обеспеченности годового стока р. Москвы — г. Звенигород.
1 — р а с п р е д е л е н и е Г у м б е л я , 2 — б и н о м и а л ь н о е р а с п р е д е л е н и е .
Т а б л и ц а 2.8
Средние значения параметров у ( п) |
и а у(п) |
при разном числе членов ряда п |
|||||||
|
|
|
|
(по |
Гумбелю) |
|
|
|
|
п |
У (л) |
° у («) |
п |
|
У (л) |
° у <"> |
.п |
У (л) |
° у (") |
20 |
0,524 |
1,063 |
40 |
0,544 |
1,141 |
60 |
0,552 |
1,175 |
|
22 |
0,527 |
1,076 |
42 |
0,545 |
1,146 |
65 |
0,554 |
1,180 |
|
24 |
0,530 |
1,086 |
44 |
0,546 |
1,150 |
70 |
0,555 |
1,185 |
|
26 |
0,532 |
1,096 |
46 |
0,547 |
1,154 |
75 |
0,556 |
1,190 |
|
28 |
0,534 |
1,105 |
48 |
0,548 |
1,157 |
80 |
0,557 |
1,194 |
|
30 |
0,536 |
1,112 |
50 |
0,548 |
1,161 |
85 |
0,558 |
1,197 |
|
32 |
0,538 |
1,119 |
52 |
0,549 |
1,164 |
90 |
0,559 |
1,201 |
|
34 |
0,540 |
1,126 |
54 |
0,550 |
1,167 |
95 |
0,559 |
1,204 |
|
36 |
0,541 |
1,131 |
56 |
0,551 |
1,170 |
100 |
0,560 |
1,206 |
|
38 |
0,542 |
1,136 |
58 |
0,552 |
1,172 |
С О |
0,577 |
1,282 |
дений 1924— 1970 гг. с применением кривой распределения Гумбеля, Исходные данные представлены в виде эмпирической кривой обес печенности на рис. 2.14. Стандартные параметры распределения оп
1 1 0
ределены графоаналитическим методом, рассматриваемым в гла
ве III. Они равны: х = Н 9 м3/с; CUjc = I,14; Cs = 2,5.
По табл. 2.8 находим среднее значение и среднее квадратическое
отклонение вспомогательной переменной у при п = 47: г/(47)=0,55;
о У(47)=1,16. Далее, используя уравнения (2.95) и (2.96), находим:
д = х — ^ - ^ - = 1 4 9 - 8 0 ,5 = 6 8 ,5 . |
||
1 |
а |
' |
Окончательные расчеты с использованием распределения Гумбеля приведены в табл. 2.9, значения у Р взяты из табл. 2.7.
Т а б л и ц а 2.9
Расчет максимальных расходов воды дождевого стока различной обеспеченности с использованием кривой Гумбеля
|
|
р. М осква — г. |
Звенигород |
р . А р х а р а - - с . |
А ркадьевка |
р |
УР |
1 |
V |
1 |
у/> |
|
|||||
|
|
|
х р = — + ч |
уР— |
х р = — |
0,1 |
6,90 |
1016 |
1084 |
2037 |
2682 |
0,5 |
5,30 |
779 |
848 |
1562 |
2206 |
1,0 |
4,60 |
677 |
745 |
1356 |
2000 |
3,0 |
3,49 |
514 |
582 |
1029 |
1674 |
5,0 |
2,97 |
437 |
505 |
876 |
1520 |
10,0 |
2,25 |
331 |
399 |
670 |
1314 |
20,0 |
1,50 |
221 |
289 |
442 |
1087 |
30,0 |
1,03 |
152 |
220 |
304 |
948 |
40,0 |
0,672 |
98,8 |
167 |
198 |
843 |
50,0 |
0,367 |
51,0 |
122 |
108 |
753 |
60,0 |
0,087 |
12,8 |
81,3 |
25,6 |
670 |
70,0 |
-0 ,186 |
-2 7 ,4 |
41,1 |
-5 4 ,8 |
590 |
80,0 |
-0,476 |
-7 0 ,0 |
-1 ,5 3 |
-1 4 0 |
504 |
90,0 |
-0 ,8 3 4 |
-123 |
-5 4 ,2 |
-246 |
400 |
95,0 |
- 1 ,1 0 |
-1 6 2 |
-9 2 ,9 |
-3 2 3 |
320 |
99,0 |
-1 ,5 3 |
-225 |
-156 |
-4 5 0 |
194 |
99,9 |
-1 ,9 3 |
-293 |
-2 9 3 |
-5 7 0 |
75,0 |
Сопоставление величии, рассчитанных с использованием кривой распределения Гумбеля, с эмпирическими данными (рис. 2.14) по зволяет обнаружить определенное несоответствие в их расположе нии. Больше того, начиная с обеспеченности 78% аналитическая кривая Гумбеля уходит в отрицательную область. Это объясняется тем, что у кривой распределения Гумбеля коэффициент асиммет рии, как указано выше, является величиной постоянной (Cs=l,14), в то время как в данном случае Cs~ 2,5.
I l l
Биномиальная |
кривая распределения, представленная на |
рис. 2.14, вполне |
удовлетворительно соответствует эмпирическим |
данным и во всем диапазоне обеспеченностей положительна. Интересно сопоставить распределение Гумбеля с данными фак
тических наблюдений, для которых Cs < 1,14. Для этой цели исполь зовались максимальные расходы воды р. Архары у с. Аркадьевки с 1941 по 1968 г.
Стандартные статистические параметры этого ряда равны: х =
= 801; Cv =0,40; Cs . = 0,50 при п= 24. Как и ранее, у (24) —0,530;
Рис. 2.15. Эмпирическая и аналитические кривые обеспеченности годового стока р. Архары — с. Аркадьевна.
1 — р а с п р е д е л е н и е Г у м б е л я , 2 — б и н о м и а л ь н о е р а с п р е д е л е н и е .
0 i/(24)~ 1,09; 1/а = 295 и ^= 645. Последующие расчеты представлены
в табл. 2.9, а сопоставление рассчитанных и наблюденных величин выполнено на рис. 2.15.
Оценивая возможности рассматриваемого распределения для описания статистических совокупностей гидрометеорологических величин, следует иметь в виду указанные выше основные допуще ния, принятые при его построении. Эти допущения применительно
ксовокупности расходов воды следующие:
1)годовая совокупность суточных расходов воды имеет распре
деление типа экспоненциального;
2)число элементов выборки, из которой извлекается наиболь ший максимальный расход воды, равное 365 дням, достаточно для применения асимптотической теории (n -э-оо);
3)элементы выборки (суточные расходы воды) взаимно неза висимы.
Указанное и некоторые дополнительные допущения приводят к получению в итоге распределения, характеризующегося постоян ным значением коэффициента асимметрии Cs= 1,14.
Рассмотрим, насколько указанные ограничения приемлемы для рассматриваемого нами класса задач.
112
Прежде всего |
отметим, |
что положение, сформулированное |
в п. 3, находится |
в очевидном |
противоречии с особенностями ре |
жима большинства рек. Наличие достаточно высокой корреляции между суточными расходами воды приводит к резкому уменьшению объема независимой информации, которая может оказаться совер шенно недостаточной для применения асимптотической теории при
п —у- О О .
Возможность использования первого допущения не может быть заранее отвергнута, но одновременно следует учитывать, что оно принято па основе общих соображений, а не выявлено на основе анализа статистических совокупностей гидрометеорологических ха рактеристик.
Указанные соображения, как отмечает Ю. Б. Виноградов [36], еще в большей мере относятся к случаю применения кривой Гумбеля для вычисления годовых максимумов суточных осадков раз личной обеспеченности.
Важной особенностью рассматриваемого распределения, суще ственно уменьшающей пределы его использования в гидрологии, является постоянство коэффициента асимметрии (С*=1,14). В силу этого кривая Гумбеля применительно к статистическим совокупно стям, обладающим меньшей асимметрией, в области малых обеспе ченностей будет давать завышенные значения, а для рядов с коэф фициентом асимметрии, существенно большим, чем Cs=l,14, — за ниженные, что иллюстрируется на приведенных выше примерах.
В зоне больших обеспеченностей, очевидно, будет наблюдаться обратное соотношение между эмпирической кривой и распределе нием Гумбеля. Только в частном случае, когда коэффициент асим метрии статистической совокупности не существенно уклоняется от единицы, распределение Гумбеля покажет достаточное согласова ние с эмпирическими данными.
Вследствие указанных особенностей кривая распределения Гум беля не получила широкого применения в советской гидрологии. Отмеченные свойства кривой Гумбеля, вытекающие из ее теорети ческой схемы, на ряде примеров подтверждены в статье Д. Л. Соко ловского и В. А. Шелутко [135].
§ 8
нормальный закон распределения
Нормальный закон распределения, часто называемый распреде лением Гаусса или Гаусса—Лапласа, находит широкое применение при решении многих вопросов, связанных с исследованием законо мерностей случайных величин. Прежде всего это наиболее часто встречающийся закон распределения случайных величин. Кроме того, к нормальному закону могут быть сведены, как к пределу, многие другие законы распределения при некоторых условиях, ти пичных для формирования различных совокупностей случайных величин. Следовательно, по отношению к таким законам
8 З ак . № 88 |
113 |
распределения нормальный закон выступает в форме предельного закона.
Нормальный закон возникает в том случае, когда переменная величина формируется под действием суммы большого числа неза висимых (или слабо зависимых) факторов при условии, что каждый из этих факторов не оказывает на изучаемое явление превалирую щего влияния. Если это условие не выполняется, то закон распре деления превалирующего фактора окажет заметное влияние на сумму и определит уклонение закона распределения суммы от нор мального распределения. Возникнув одним из первых, рассматри ваемый закон к настоящему времени наиболее разработан в теоре тическом отношении и потому широко применяется при решении многих задач статистики и теории вероятностей.
Нормальное распределение было получено в связи с анализом ошибок измерений, но впоследствии нашло применение для описа ния статистических совокупностей многих природных явлений. В ча стности, нормальный закон распределения был использован А. Хазеном [152, 153] для описания многолетних колебаний речного стока.
Применительно к исследованию статистических закономерно стей речного стока в последующем стали использоваться иные за коны распределения. Однако при описании ряда других статисти ческих совокупностей гидрологических величин нормальный закон находит непосредственное применение.
Например, это имеет место при изучении турбулентных пульса ций скоростей течения и гидродинамических усилий, возникающих
впотоке, распределения высот снежного покрова на маршруте или
впределах некоторой территории, ошибок измерений параметров
гидрологического режима и т. д. М. А. Великанов [33] использовал нормальный закон распределения для обоснования теории форми рования донных гряд, И. Ф. Карасев [60] — для оценки вероятности срыва частиц грунта со дна турбулентного потока, К. И. Россий ский [113—115] — при исследовании закономерностей движения донных наносов.
Подробный вывод нормального закона распределения содер жится почти во всех курсах теории вероятностей и математической статистики и поэтому здесь не излагается. Укажем лишь, что он, например, может быть получен из приведенного ранее дискретного биномиального закона распределения при оо.
Уравнение нормальной кривой распределения в дифференциаль ном виде или в форме кривой распределения плотности вероятно стей имеет вид
,U-3H2
Я(*) = - 7 7 й Г е |
" |
• |
<2-97> |
где х — математическое ожидание (среднее |
значение) |
переменной |
|
х; а — среднее^квадратическое отклонение. |
|
|
|
Величины х и а являются параметрами нормальной кривой. Пре дел простирания нормальной кривой распределения — от минус до плюс бесконечности (—о о < х < о о ) .
114
Кривая нормального закона распределения располагается сим
метрично относительно максимальной ординаты, равной у = ---- —,
а^2п
находящейся в точке х = х. Следовательно, среднее арифметическое
значение (х) |
и мода |
(Л10) в нормальном распределении совпадают. |
В силу этого |
с ними |
совпадает и медиана (Me). По мере удаления |
от точки х ординаты кривой распределения уменьшаются и при
х ztоо стремятся к нулю. При х = х ± а функция Р(х) имеет точки
перегиба.
Непосредственно из форму лы (3.97) следует, что величи
на х в нормальном законе рас |
|
|||||
пределения |
является |
центром |
|
|||
симметрии. Действительно, при |
|
|||||
изменении |
знака |
разности |
|
|||
(.V — х) на обратный |
выраже |
|
||||
ние (2.97) |
не меняется. Функ |
|
||||
ция Р (х) является четной, т. е. |
|
|||||
Р(х) = Р ( —х). |
Из |
симметрич |
|
|||
ности нормального закона сле |
|
|||||
дует, что все. нечетные |
цент |
|
||||
ральные моменты равны нулю. |
|
|||||
Соответственно |
равен |
нулю и |
|
|||
коэффициент асимметрии. |
Для |
Рис. 2.16. Нормальные кривые распре |
||||
четных моментов |
имеем |
соот |
деления при различных средних квад |
|||
ношения Ц2= о2; |
Ц4= Зо4. |
Из |
ратических отклонениях. . |
|||
выражения |
четвертого |
мо |
|
|||
мента следует, |
что эксцесс нормального распределения равен нулю |
£ = - § - - 3 = 0 .
При изменении величины математического ожидания (х) кривая
распределения будет смещаться вдоль оси абсцисс, не изменяя
своей формы. При х = 0 имеем семейство центрированных нормаль ных кривых (т. е. с центром в начале координат). Параметр о яв
ляется характеристикой размаха отклонений переменной величины
х от центра рассеяния (х) и, следовательно, характеристикой формы
кривой распределения. Действительно, максимальная ордината, как указано выше, обратно пропорциональна а и, следовательно, умень шается при возрастании о. Это при условии неизменности площади ограниченной кривой приводит к тому, что с возрастанием а кривая
распределения становится все более плоской, растягиваясь вдоль оси абсцисс. При уменьшении а, наоборот, кривая распределения сжимается с боков и вытягивается вверх вдоль оси симметрии
(рис. 2.16).
8* |
115 |