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

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

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

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

Добавлен: 15.10.2024

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

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

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

(1940 г.) по критерию Стьюдента, использование которого предпо­ лагает равенство средних квадратических отклонений рассматри­ ваемых рядов. Напомним, что за период наблюдений с 1877 по 1940 г. среднее квадратическое отклонение оказалось равным ау=

= У 75 600, а с

1941 по 1955 г. — 0Ж= У721ОО. Требуется опреде­

лить, являются

ли полученные различия в выборочных значениях

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

Для оценки однородности средних квадратических отклонений воспользуемся критерием F. По формуле (4.11) определяем

F = 75 600

 

 

 

 

72 100 ~ ' 1*'

 

 

 

 

Этот критерий подчиняется

распределению

Фишера со степе­

нями свободы &i= 64— 1=63 и /е2= 15— 1 = 14.

Задаваясь

уров­

нями значимости 10 и 2% и полученными

значениями ki и

k2, по

таблице распределения F [89]

при <7= 5 и

1% определяем

крити­

ческие области Лер, которые для &i= 63

интерполированы

между

значениями &i= 50 и &i = 75 и равны при

уровне

значимости 10%

Рю%=2,23, при уровне значимости 2% F

 

=3,18. Следовательно,

выборочное значение Z7 = 1,05 лежит в области допустимых значе­

ний при любом из выбранных нами уровне

значимости. Поэтому

предположение о равенстве средних квадратических отклонений не противоречит наблюденным данным о годовом стоке р. Волги у г. Ярославля.

Используем критерий F для оценки однородности средних квад­

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

Среднее квадратическое отклонение максимальных

расходов

воды за 1877—-1940 гг. составило 0 i= У 3 354 000, а

за 1941—

1955 гг. 02 = У 795 200 . По формуле (4.11) получаем

 

3 354 000

4,22.

 

 

795 200

 

 

 

 

 

В качестве нулевой гипотезы примем 01 = 02,

а в качестве аль­

тернативной 0 1^=02. Так как число членов ряда

в каждом из рас­

сматриваемых периодов точно такое, как и для годового стока, по­ лучаем ki = 63 и *2 = 14. Следовательно, критические значения F1(p

при 10 и 2%-ных уровнях значимости будут такими же, как и в пре­ дыдущем примере, и равны F 10%= 2,23 и F 2%=3,18.

Полученное значение /■’-критерия, равное 4,22, даже при 2%-ном уровне значимости лежит в критической области (FKp< F ).

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

201


быть признаны однородными. Иначе говоря, принятая ранее нуле­ вая гипотеза отвергается, а признается альтернативная гипотеза о неоднородности средних квадратических отклонений максималь­ ного стока до и после создания Рыбинской ГЭС.

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

В качестве упрощенного критерия проверки однородности ряда средних квадратических отклонений можно использовать G2-Kpirre- рий, который представлен отношением

О2 =■

+

(4.12)

°1 + °2 +

 

где Омане — наибольшее из эмпирических средних

квадратических

отклонений; оi, 02....... ok — рассчитанные по наблюденным данным

средние квадратические отклонения.

Данный критерий однородности применим, когда рассматривае­ мые наблюденные ряды имеют один и тот же объем. В монографии

[89]приведено распределение критерия G2Kp при числе выборок k

ичисле наблюдений в каждой выборке п при 5 и 1%-ных уровнях

значимости. Критическая область будет G2p < G 2.

Вкачестве примера использования критерия (4.12) оценим од­ нородность дисперсий запасов воды в снежном покрове по данным,, приведенным в табл. 4.1,

О2 = __________ ________ ____ .—о 98

и

590 + 218 +

967 + 676 +

1129

’ '

При 5%-ном уровне значимости критическое значение G2=0,46;

полученное значение критерия

G2 = 0,28

попадает в область допу­

стимых значений, и следовательно, дисперсии рассматриваемых ря­ дов признаются однородными.

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

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

202


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

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

Для примера рассмотрим максимальные годовые расходы воды р. Стрый у г. Турка и р. Абава у х. Сисени, приведенные в табл. 4.2. Из этих данных видно, что максимальные расходы воды в рассмат­ риваемых створах в одни годы сформированы в результате весен­ него снеготаяния (отмеченные скобками), а в другие годы — за счет дождей. По этому признаку каждая выборка может быть разде­ лена на две части. На р. Стрый максимальные годовые расходы воды в течение 43 лет наблюдений 24 раза формировались в весен­ ний период и 19 раз — в летне-осенний; на р. Абава весенние рас­ ходы наблюдались 20 раз, а дождевые — 15 раз.

Для выделенных совокупностей весенних и летне-осенних мак­ симальных расходов воды были вычислены параметры распределе­ ния, приведенные в табл. 4.3.

Выясним, относятся ли полученные значения параметров к од­ нородным или неоднородным совокупностям.

Применяя критерий Фишера, оценим в отношении однородности величины средних квадратических отклонений:

для р. Стрый

с-

авес

9988

, АГ1

 

 

 

F = - r~ = - W f T ^

1,02,

 

 

 

 

дожд

 

 

 

 

 

для р. Абава

 

 

 

 

 

 

 

F

=

-

1,00.

 

 

 

Критическая область применительно к ряду

по

р. Стрый

при

5%-ном уровне значимости

равна F 5%= 2,07;

она

получена

по

203


Т а б л и ц а 4.2

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

р . С тры й — г. Т у р к а , У7 = 897 км 2

р . А бава — х. С исени,

F =

1990 км 2

ГОД

^ м а к с

дата

г о д

^ м а к с

 

 

 

М 3 / С

 

 

 

 

 

 

 

 

1907

(248)

30/IV

1927

(137)

 

8 /111

1908

(215)

4/111

1928

183

20/V; 22/VI

1900

(296)

24/III

1929

105

 

20/XI

1910

(108)

24/11

1930

179

 

8 /VI11

1913

335

3/V III

1931

(316)

 

24/IV

1914

(143)

11/111

1932

(207)

 

9/IV

1916

239

19, 22/VII

1933

(146)

 

17/III

1917

(114)

1/1V

1934

103

 

21/1

1918

281

9/1

1935

291

 

23/11

1922

(174)

3/111

1936

(107)

 

14/1II

1923

88,4

6/VI

1937

(155)

 

25/III

1924

(446)

25/1II

1938

57,8

 

6/II

1925

268

28/VI

1939

89,1

 

11/11

1926

385

24/X

1940

(191)

 

15/1V

1927

515

31 /V III

1941

(188)

 

16/IV

1928

(238)

з о /ш

1942

(303)

 

12/ IV

1929

222

14/VII

1943

61,0

 

18/11

1931

222

26/X

1944

(125)

 

11/V

1932

(330)

6 /IV

1945

(119)

 

4 /IV

1933

239

8/V II; 5/XI

1946

(164)

 

24/III

1934

(198)

12/1II

1947

(146)

 

1 /IV

1935

(150)

12/IV

1948

235

 

4/II

1936

219

23/1

1950

89,4

 

16/XI

1937

239

20/VI11

1951

(306)

 

9/ IV

1938

(222)

10/III

1952

(149)

 

12/ IV

1939

198

13/VI

1953

(158)

 

27/1II

1942

(114)

20' IV

1954

62,6

 

22/Х

1944

(222)

18/IV

1955

172

 

1/V

1946

(93)

6/IV

1956

(258)

 

26/IV

1947

(299)

23/II

1957

94,6

 

21/ 1X

1948

395

17/VII

1958

(174)

 

18/IV

1949

(359)

8/IV

1959

(79,6)

 

9/111

1950

(77,3)

28/II

1960

88,6

 

7/XII

1951

(137)

15/1II

1961

101

 

7/XII

1952

(396)

2/IV

1962

(242)

 

8 /IV .

1953

155

1/1

 

 

 

 

1954

(108)

5/111

 

 

 

 

1955

(206)

25/1II

 

 

 

 

1956

(143)

15/1V

 

 

 

 

1957

190

23/VI

 

 

 

 

1958

213

23/IX

 

 

 

 

1959

121

4/V III

 

 

 

 

1960

108

28/VII

 

 

 

 

П р и м е ч а н и е .

В скобках приведены

максимальные расходы

воды весен­

него половодья.

204


Т а б л и ц а 4.3

Параметры статистических выборок снеговых и дождевых максимальных расходов воды

Фаза режима

р. Стрый—г. Турка

р. Абава--х. Сисени

X

а

X

о

 

 

Весеннее

половодье

210

98,9

184

66,7

Дождевые

паводки

244

99,8

128

66,6

таблице,

приведенной в работе [89],

с учетом того, что число степе­

ней свободы для большей дисперсии равно

ki — tii— 1 = 18, для

меньшей дисперсии кг=П2 — 1 =23.

 

 

 

Применительно к ряду по р. Абава критическая область равна

F5%=2,26 при ki — Ы и

= 19. В обоих случаях в критическую об­

ласть при 5%-ном уровне значимости попадают значения критерия F, большие, чем полученные нами. Это значит, что в рассматривае­ мом случае критерий однородности F попадает в область допусти­

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

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

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

Условие однородности дисперсий было установлено выше. Усло­

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

может быть проверено

с помощью критериев согласия

эмпирического и аналитического

распределений, рассматриваемых

в § 3

настоящей главы. Здесь

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

205