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

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

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

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

Добавлен: 15.10.2024

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

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

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

при т = О

 

 

Р(х )= е -«(*-хУ‘

( х о ^ х < о о ) ,

(2.67)

при п = т

 

 

_ т ( х ~ -VnУ”

(*0< х < А Г ).

(2.68)

/>(*)=<?

При использовании уравнений (2.67) и (2.68) непосредственно для аппроксимации различных кривых, т. е. вне связи с оценкой статистических параметров С„ и Cs, обеспеченность Р удобно вы­

разить в процентах, а за основание показательной функции взять

число 10.

 

 

и (2.68)

будут иметь вид:

В этом случае уравнения (2.67)

р — 1Q2

-

(Х~ Х°)П

(2.69)

 

 

 

Х — Хд

 

/>=10

2 -

1 ' (Х - х г

(2.70)

где а' = 0,43а; у' = 0,43у.

Для определения параметров уравнений (2.69) и (2.70) Гудрич предложил специальную клетчатку, названную им клетчаткой асим­ метричной частоты.

Ось абсцисс этой клетчатки представляет шкалу значений (2 —

— 1gP), а ось ординат — логарифмическую шкалу. Смысл построе­

ния данной клетчатки заключается в том, что уравнения

(2.69) и

(2.70) после двухкратного логарифмирования

 

lg (2 — lg До)— lg я'-\-п lg (х — х0),

(2.71)

lg (2 -lg /> 0) = l g T ' + m l g ( ^ ^ )

(2.72)

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

том п (или т) и свободным членом lga'

(или lg y 7),

если на оси

ординат клетчатки откладывать значения

/

X — Xq \

х х0 ^ или —— — у.

Значения х0 и X подбираются при этом

из условия

наилучшего

спрямления эмпирической кривой обеспеченности, определяемой

точками: Р(хi), Р(х2), ..., Р(хп), где Р(хп) = - - - т ^ • 100%; п — ТЬ“г

число всех членов ряда; т — порядковый номер ряда, расположен-

ного по убывающим значениям признака.

Уравнение (2.71) использовал Б. Д. Зайков [55] для спрямления эмпирических кривых обеспеченности максимальных расходов ве­ сенних половодий, а уравнение (2.72) применял В. А. Урываев [139] для аппроксимации кривых обеспеченностей суточных расходов воды.

101


Уравнения (2.67) и (2.68) в дифференциальной форме имеют вид:

dP

по ( х —х 0)п е

( х - х 0) п

(2.73)

dx

 

 

 

( х — х а) \ т

 

=

( х - х 0)

е

(2.74)

Графическое изображение функции

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

(2.73) при

различных значениях параметра п представлено на рис. 2.11.

р(х)

Рис. 2.11. Кривые распределения Гудрича при раз­ личных n= f(Cs, Cv).

l - n - 1,0.

С е = 2,0,

С „ = 1,0;

2 — л = 1,4,

С, —1.19,

С „ = 0,72;

3 — п -2 ,0 ,

Cs= 0,63,

С „ = 0,52;

4 — л = 3,6,

Са=0,0,

С „ = 0,31;

 

5 — /1=6,0, Са«— 0,37, 0^ = 0,19.

 

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

В силу указанного Алексеев рассмотрел детально уравнение (2.73) и составил для него стандартные таблицы нормированных

отклонений ——— ординат кривой обеспеченности Гудрича, анало-

С V

гично тому, как Фостер составил таблицу нормированных отклоне­ ний ординат кривой Пирсона III типа.

Используя обычный прием выражения параметров кривой рас­ пределения через моменты, можно установить связь параметров уравнения (2.73) а, я и хо со средним арифметическим значением

и коэффициентами вариации и асимметрии.

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

1 0 2


фициент асимметрии, связанный с параметром п, следующим соот­

ношением:

г ( 1 + 4

) - з г ( ,

+ Х ) г ( 1+ А ) + гГз(, + Т-)

 

3

г /

2 \

/ 1 \1г/3

‘ *

И '+ тгЬ 'Ф +тг)]

Расчет нормированных значений отклонений х от среднего зна­

чения х осуществляется по формуле

|1п/>|,/я- Г (l + 4-)

(2.76)

^ r (1 + A ) - r2( , + _ L ) ’

которая непосредственно вытекает из уравнения (2.73) и из сле­ дующих соотношений:

3 - * „ = - 1

;гг ( 1 + У ) ,

(2.77)

°2= т Ы Г ( ‘ +

4 - ) - Г !( 1 + Ф ) ] -

(2.78)

Заменяя натуральный логарифм на десятичный и выражая обес­ печенность Р в процентах, формулу (2.76) запишем в виде

J ^ L = =A ( 2 - \ g P ) ' ln- B ,

(2.79)

где

________ (2,3026)11"_________

/ r(1+T)-r2(4 -v) ’

] / r ( 1 + A ) _ r !(, + J_)

причем А и В являются функциями Cs, так как Cs=(/(n). Задаваясь величиной параметра п от 0,5 до 20, на основании выражений (2.75),

(2.78) и (2.79)

Алексеев составил таблицу нормированных отклоне-

.

k — \

зависимости от обеспеченности Р и коэффициента

нии

— —— в

С v

асимметрии Cs = f(n). Эта таблица приведена в работе [10]. Полагая в формуле (2.79) х = х0,Р = 100%, получаем

Х р — X

*0-1

В,

k p = \ - B C v.

о

Су

 

 

103


Из приведенного равенства следует, что нормированное отклоне­

ние нижней границы распределения х0 от середины х равно пара­ метру В, взятому со знаком минус. Из приведенного соотношения

также следует, что при заданной асимметрии распределения (за­ данному Cs отвечает определенное значение В) нижняя граница

кривой распределения Гудрича: положительна при

ka= \ —BCv > 0 , если С» < -g- ,

равна нулю при

k0= \ —BCv=0, если Cv = ,

и, наконец, отрицательна при

&о=1 — BCV< 0, если

В соответствии с этим величину \/B = CVq можно принять в ка­

честве критерия для оценки области положительных значений кри­

вой Гудрича в зоне больших обеспеченностей при заданной

асим­

метрии.

 

[10] построена зависи­

На основании данных таблицы Алексеева

мость Cs= f(Cv0)

Для кривой распределения

Гудрича (рис.

2.12).

Эта зависимость,

приближенно описываемая уравнением

Cs=

= 2,9 Cv — 0,9, определяет границу совпадения наименьшего значе­ ния признака с нулем. Иначе говоря, при соотношениях между Cv

и Cs, соответствующих этой линии или располагающихся выше ее, кривая Гудрича не уходит в отрицательную область. Для сопостав­ ления на рассматриваемом графике нанесена линия, соответствую­ щая равенству CS=2CV, при котором кривая Пирсона III типа не

уходит в отрицательную область. При значении С„ = 1 эти линии пересекаются. Из графика следует, что нижняя граница кривой рас­ пределения Гудрича, в отличие от кривой Пирсона III типа, ос­ тается положительной и при CS<2CV, если только C„<CS=1/,B, или, что то же самое, если (Д>2,9 Cv — 0,9.

Однако при Cv > 1 соотношение между рассматриваемыми кри­

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

Сопоставим распределение Гудрича с распределением Пирсона

III типа при значениях С^ = 0,5; 1,0; 1,5 и при CS = 2C„ и CS= C V.

Как и следовало ожидать, при Сг- = 1,0 и при CS = 2C„ оба распре­

104


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

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

При Си = 0,5 и Cs = 0,5 распределение Пирсона III типа уходит в отрицательную область при Р = 99%, а распределение Гудрича

при Р = 99,9% (рис. 2.9 6

и рис.

2.13). При С^ = 1,5 и CS = 2CV рас­

пределение Пирсона III типа положительно во всем диапазоне обес­

печенностей,

а кривая

Гудрича

уходит в область отрицательных

значений признака при Р = 75%.

 

Таким образом,

при

Cv<

 

<1,0

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

Гудрича

 

в случае, если CS<2C„, может

 

оказаться

предпочтительней

 

по сравнению с распределени­

 

ем Пирсона III типа.

 

 

 

 

При Со>1,0 и при CS<2CV

 

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

Гудрича

ухо­

 

дит в отрицательную

 

область

 

при меньших значениях

обес­

 

печенности

по

сравнению с

 

кривой Пирсона III типа.

 

 

Техника

расчета

 

ординат

 

кривой

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

Гудри­

 

ча по существу

не отличается

 

от расчета ординат

распреде­

 

ления

Пирсона

III типа, так

 

как расчетные таблицы в обо­

 

их случаях

содержат

 

значе­

 

ния

(Р, с.).

С у

 

Пример

расчета величин

годового стока различной обес­ печенности р. Днепра у пгт Лоцманской Каменки приве­ ден в табл. 2.6.

Рис. 2.12. Зависимость соответствен­ ных значений Cs= f (С„), при которых функции распределения проходят че­ рез нулевое значение признака рас­ пределения.

I —биномиальное распределение, 2—рас­

пределение Гудрича.

Расчет величин годового стока р. Днепра у пгт Лоцманской Каменки с использованием распределения Гудрича

Р°10

...................

. 1

3

5

10

50

60

80

90

99

kp =

уpCv -f- 1

. .

1,73

1,58

1,50

1,38

0,98

0,90

0,75

0,65

0,49

х р ~

kpx ■ •

.

6,54

5,97

5,67

5,22

3,70

3,40

2,84

2,46

1,85

1 0 5