Файл: Рождественский, А. В. Статистические методы в гидрологии.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 15.10.2024
Просмотров: 185
Скачиваний: 0
В уравнениях интегральной биномиальной кривой распределе ния в виде (2.45) и (2.47) величина обеспеченности признака рас
пределения |
связывается с величинами коэффициентов |
|
х |
вариации и асимметрии.
Параметры более высоких моментов, чем третий, функционально связаны с коэффициентами асимметрии и вариации. Воспользу емся этим для определения эксцесса через величины параметров Cv и Cs, для чего в рекуррентное соотношение для моментов (2.29)
подставим (2.30)
|
|
|
|
|
|
2mi |
-т4-{-■2тп |
|
|
|
|
|
|
|
Зт2-\-- т2 |
|
|||
Из полученного уравнения определен четвертый момент |
|||||||||
|
|
|
|
6пё, + 4 |
— тз |
3 (2т\ -)- m^j |
|
||
|
|
nli |
|
|
2т2 |
|
2т2 |
|
|
Следовательно, |
выражение эксцесса можно представить в виде |
||||||||
р _ |
т 4 |
з - |
3 |
|
2т\ + т\ \ |
о |
3 (2т\ + т\ — 2яг®) |
2>т\ |
|
С |
- |
1 - |
|
|
|
|
|||
„ 2 |
й |
2 |
1 |
А |
|
|
2т\ |
2т\ |
|
|
т2 |
|
|
|
) |
|
|||
Наконец, выражая моменты через |
параметры Cv и |
Cs, полу |
|||||||
чаем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
^ = - |- С*> |
|
(2.48) |
|
а при |
CS = 2CV |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
E — 6Cl. |
|
(2.49) |
Аналогичным образом через Cs и С„ могут быть выражены лю бые моменты (т 3 ).
Численное решение уравнений (2.45) и (2.47) в каждом конкрет ном случае, т. е. при изменении параметров эмпирического ряда, является достаточно сложной задачей. Поэтому возникла необхо димость составления стандартной расчетной таблицы рассматри ваемого семейства кривых распределения, пользуясь которой можно определить ординаты кривой обеспеченности в зависимости
от величин параметров х, Сv и Cs. Впервые такую таблицу соста
вил А. Фостер [151]. В СССР она была уточнена С. И. Рыбкиным [117], а затем и другими авторами.
Эта таблица в современном ее виде опубликована в «Руковод стве по определению расчетных гидрологических характеристик». Указанная таблица1 представляет собой численное решение урав-
1 В дальнейшем при упоминании этой таблицы указание на необходимость использования Руководства не повторяется.
90
нения (2.45) при C\,= l и различных р (выраженных в процентах) и Cs. В таблице приведены величины отклонения модульных коэф
фициентов { ^ = ~^J от еДиницы (т- |
е- среднего значения |
при |
||
знака), выраженные в долях коэффициента вариации1 (СД, |
т. е. |
|||
|
|
Xп |
|
|
yP= f ( c s, р)-~ |
1 |
- 1 |
(2.50) |
|
|
||||
СV |
Су |
|||
|
|
Следовательно, для определения ординат биномиальной кривой
обеспеченности по заданным параметрам х, Cv и Cs табличные зна чения величины ур, полученные при заданном Cs для различных
обеспеченностей р%, необходимо преобразовать в соответствии с равенством (2.50) к виду
kp—УрСу-\-^>
(2.51)
х р= { урС„-{-\)х.
Если среднее значение рассматриваемого ряда мало отличается от нуля (или равно ему), коэффициент вариации становится не определенной характеристикой. В таком случае для расчета вели чин различной обеспеченности необходимо использовать величину среднего квадратического отклонения. Сказанное вытекает из сле дующих очевидных соотношений:
kP- 1
Ур Со
или
хр= у рз-\-х.
Приведем пример. Для статистической совокупности величин годового стока р. Днепра у пгт Лоцманской Каменки получены сле
дующие значения параметров: <7= 3,78 л/с-км2; = 0,28; Cs = 0,46.
Рассчитаем модули стока различной обеспеченности (т. е. значения ординат биномиальной кривой обеспеченности) при указанных зна чениях параметров. Исходные величины у р при различных р для
Cs= 0,46 взяты из таблицы биномиального закона. Порядок рас чета дан в табл. 2.4.
Полученные значения, определяющие очертание теоретической кривой обеспеченности, целесообразно сопоставить с эмпирической кривой. Такое построение обычно выполняется на клетчатках веро ятностей, рассматриваемых в главе III. При этом эмпирическая
! Табличное значение величин |
kp— 1 |
в гидрологической литературе часто |
||
Cv |
||||
|
kP- |
|
||
обозначают символом <!>(*) = |
1 |
|
||
|
С, |
|
|
Т а б л и ц а 2.4
Порядок расчета величин годовых модулей стока различной обеспеченности (qp) р. Днепра у пгт Лоцманской Каменки
О беспеченность Р %
1 |
3 |
5 |
10 |
50 |
60 |
80 |
90 |
99 |
Ур УpCv
k — УpCv + О ex II '£
|
2,64 |
2,06 |
1,76 |
1,32 |
-0,075 |
-0 ,3 2 |
—0,85 |
-1 ,2 2 |
- 2 ,0 |
|
0,74 |
0,58 |
0,49 |
0,37 |
-0,021 |
-0 ,0 9 |
-0 ,2 4 |
-0 ,3 4 |
-0 ,5 6 |
I |
1,74 |
1,58 |
1,49 |
1,37 |
0,98 |
0,91 |
0,76 |
0,66 |
0,44 |
|
6,58 |
5,20 |
5,63 |
5,18 |
3,70 |
3,44 |
2,87 |
2,49 |
1,66 |
•обеспеченность членов статистического ряда обычно определяется по формуле
■100°/о, |
(2.52) |
где т = 1 , 2 , . . « — порядковые номера членов статистического ряда, расположенных в убывающем порядке; п — общее число чле
нов статистической совокупности.
При сопоставлении аналитической и эмпирических кривых обес печенностей (рис. 2 .8 ) обнаруживается достаточно хорошая их
сходимость. Рядам годового стока и некоторым другим характери
стикам гидрологического |
режима свойственно соотношение Cs — |
= 2Cv. Для этого частного, |
но часто используемого соотношения |
с целью облегчения расчетов составлена специальная таблица ор динат кривых обеспеченностей, выраженных в модульных коэффи циентах в зависимости от Cv и р \
yp=*kp= ^ £ - = f ( C B, р ). |
(2.53) |
Рассмотренная биномиальная кривая распределения представ ляет частный случай (III тип) обширной системы кривых К. Пир сона, состоящей из 13 типов распределения (включая нормальную кривую).
Совокупность кривых Пирсона получается в результате реше ния дифференциального уравнения вида
1 . |
аУ |
— ____ г + х____ |
/ |
/о щи |
у |
dx |
Ь0 b\X + b-2 X2 ’ |
|
V • / |
которое представляет собой более общий вид рассмотренного ра нее дифференциального уравнения биномиальной кривой распре деления (2.23). Уравнение (2.54) соответствует статистической1
1 См. «Руководство по определению расчетных гидрологических характери стик» (приложение 3, случай Cs= 2Cr ).
9 2
3,0
0,2
0,1----------- |
— ---- |
— L J — _J-------------------------— — ______ _________ ____ _________ |
|||||||
0,01 |
0,1 0,2 |
0,5 |
1 2 |
5 |
10 |
20 30 40 50 60 70 80 |
90 95 |
98 |
99 99,5 99,8P % |
Рис. 2.8. |
Аналитическая |
и эмпирическая кривые обеспеченности модульных |
коэффициентов |
годового стока |
|||||
|
|
|
|
р. |
Днепра у пгт Лоцманской Каменки. |
|
|
|
/ —эмпирические точки, 2—биномиальная кривая.
схеме, описываемой так называемым гипергеометрическим распре делением. Это распределение получается, если оценивается вероят ность появления белого шара т раз в S извлечениях, когда в рас сматриваемой совокупности объемом п имеется пр белых шаров и nq черных ( р > 0, <7 > 0, p + q=\), причем извлеченный шар снова
не возвращается в исходную совокупность (схема невозвращенного шара). В этой схеме, в отличие от схемы Бернулли, испытания становятся зависимыми и потому выводы более общие.
В геометрическом смысле уравнение (2.54), так же как и урав нение (2.23), отражает характерные особенности очертания кривых распределения, которые обычно начинаются от нуля (в отношении значений вероятности), достигают максимума и затем снижаются с той же или большей частью с иной интенсивностью.
Использование в уравнении (2.54) третьего параметра (62) при
водит к необходимости в его оценке опираться на значение четвер того момента, который не может быть установлен с необходимой точностью вследствие ограниченности статистических совокупно стей гидрологических характеристик. Область простирания боль шинства типов кривых Пирсона выходит за пределы, характерные для гидрологических характеристик (О ^ хсоо ). Поэтому в гидро логии преимущественное применение получила кривая Пирсона III типа, которая здесь и рассмотрена достаточно подробно.
Сведения о всей совокупности кривых Пирсона можно получить из работ по статистике [111, 144].
Очертание биномиальной кривой обеспеченности при различ ных Cs и Cv показано на рис. 2.9.
§ 5
кривая распределения вероятностей С. Н. Крицкого и М. Ф. Менкеля
Кривая распределения вероятностей Пирсона III типа, широко используемая в практике гидрологических расчетов, обладает, как было показано, одним существенным недостатком, заключающимся в том, что она при CS<2CV уходит в область отрицательных зна
чений гидрологических характеристик, образующих статистическую совокупность. Применительно к тем достаточно многочисленным характеристикам гидрологического режима, которые по их физиче ской сущности не могут принимать отрицательных значений, такая схематизация оказывается принципиально противоречивой.
Лишь в одном частном случае, при CS = 2C„, кривая Пирсона III типа охватывает область изменения переменной от нуля до бес конечности, что не входит в логическое противоречие с представ лениями о законах колебания речного стока и многих других гидро логических величин.
94