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