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