Исходный ряд гидрометрических наблюдений, относящийся к максимальным уровням воды, теперь можно представить в виде ряда наблюдений за максимальными расходами воды на водотоке. Эти ежегодные максимальные расходы принимаются в качестве исходного статистического ряда наблюдений для дальнейшего расчета.
3.2. Определение расходов и уровней воды требуемой вероятности превышения
Нормы проектирования мостовых переходов регламентируют максимальные значения основных гидрологических характеристик пересекаемого водотока – расхода Q и соответствующего ему уровня воды H – через вероятность их превышения р.
На железных дорогах России мостовые переходы, как и малые водопропускные сооружения, проектируются на пропуск двух максимальных расходов: расчетного и наибольшего, значения вероятностей превышения которых устанавливаются нормами проектирования сооружений [3] (табл. 3).
Таблица 3
Вероятность превышения расходов и соответствующих им уровней воды
|
Вероятность превышения расходов |
||
|
и соответствующих им уровней воды |
||
Категория железной дороги |
|
на пике паводка (%) |
|
|
|
|
|
|
Расчетного |
|
Наибольшего |
|
|
|
|
I–III |
1:100 (1) |
|
1:300 (0,33) |
|
|
|
|
IV (общей сети) |
1:50 (2) |
|
1:100 (1) |
|
|
|
|
В рассматриваемом примере категория проектируемой железной дороги – III, следовательно, расчетный расход водотока определяется с вероятностью превышения 1 %, а наибольший – с вероятностью превышения 0,33%.
Расход и соответствующий ему уровень воды заданной вероятности превышения определяется по теоретической кривой распределения вероятностей, для построения которой нужен статистический ряд многолетних наблюдений за расходами. Таким статистическим рядом являются расходы, соответствующие годовым максимумам уровня воды, зафиксированным на водомерном посту в течение 20 лет.
15
Вначале на клетчатке вероятностей строится эмпирическая кривая распределения ежегодных вероятностей превышения максимальных расходов Q, наблюдавшихся на водотоке (рис. 4). Клетчатки используют для выпрямления кривых вероятностей превышения. Клетчатку можно построить и самостоятельно, пользуясь рекомендациями учебника [2].
Для построения эмпирической кривой распределения вероятностей имеющийся ряд наблюдений Нi следует расположить в убывающем порядке. В табл. 4 в убывающем порядке записаны наблюдавшиеся значения отметок уровня воды Нi (графа 2) и соответствующих значений расхода Qi (графа 3).
Для каждого члена ряда наблюдений за гидрологической величиной определяется эмпирическая вероятность его превышения рэ (в процентах) по формуле
р = |
m |
100 % , |
(7) |
э |
n +1 |
|
где m – порядковый номер члена ряда гидрологической характеристики, расположенного в убывающем порядке (графа 1 табл. 4); n – общее число членов ряда наблюдений, в рассматриваемом примере n = 20.
Например, для первого значения (Нi = 317,6 м, Qi = 7674 м3/с):
р |
= |
m |
|
|
100 % = |
1 |
100 % = 5 %. |
|
|
|
|
||||
э |
|
n + |
1 |
|
20 +1 |
||
|
|
|
|||||
Модульные коэффициенты расходов воды определяют по формуле
ki = |
Qi |
|
, |
(8) |
|
Q |
|||||
|
|
|
|||
где Qi – значения максимального годового значения расхода воды для i-го члена ряда; Q – среднее арифметическое значение максимальных расходов воды, определяемое по формуле
n
∑Qi
|
= |
i=1 |
. |
|
|
Q |
(9) |
||||
|
|||||
|
|
n |
|
||
n
Значение ∑Qi соответствует итоговому значению графы 3 табл. 4. Для
i=1
рассматриваемого примера
16
n
Q = ∑i=1Qi = 107 084 = 5354,2 м3/с. n 20
На клетчатке вероятности по оси абцисс указывают значения вероятности р в процентах, а на оси ординат – модульные коэффициенты k. При этом следует учитывать диапазон изменений величин модульных коэффициентов ki (в рассматриваемом примере ki см. в графе 5 табл. 4), который изменяется от 0,6729 до 1,4333. Максимальное значение на шкале модульных
Таблица 4
Определение параметров эмпирической кривой распределения вероятностей и значений статистик λ2 и λ3
Поряд- |
Отметка |
|
|
Эмпириче- |
Модульный |
|
|
ковый |
Расход воды |
lg ki |
ki lg ki |
||||
номер |
уровня |
3 |
/с |
ская вероят- |
коэффици- |
||
воды Нi, м |
Qi, м |
ность рэ,% |
ент ki |
|
|
||
m |
|
|
|
|
|||
1 |
2 |
3 |
|
4 |
5 |
6 |
7 |
|
|
|
|
|
|
|
|
1 |
317,6 |
7674 |
5 |
1,4333 |
0,1563 |
0,2241 |
|
|
|
|
|
|
|
|
|
2 |
317,4 |
7331 |
10 |
1,3692 |
0,1365 |
0,1869 |
|
|
|
|
|
|
|
|
|
3 |
317,2 |
7026 |
14 |
1,3122 |
0,1180 |
0,1549 |
|
|
|
|
|
|
|
|
|
4 |
317,0 |
6707 |
19 |
1,2527 |
0,0978 |
0,1226 |
|
|
|
|
|
|
|
|
|
5 |
316,9 |
6557 |
24 |
1,2246 |
0,0880 |
0,1078 |
|
|
|
|
|
|
|
|
|
6 |
316,8 |
6405 |
29 |
1,1963 |
0,0778 |
0,0931 |
|
|
|
|
|
|
|
|
|
7 |
316,5 |
5985 |
33 |
1,1178 |
0,0484 |
0,0541 |
|
|
|
|
|
|
|
|
|
8 |
316,2 |
5565 |
38 |
1,0394 |
0,0168 |
0,0174 |
|
|
|
|
|
|
|
|
|
9 |
316,1 |
5429 |
43 |
1,0140 |
0,0060 |
0,0061 |
|
|
|
|
|
|
|
|
|
10 |
316,0 |
5297 |
48 |
0,9893 |
–0,0047 |
–0,0046 |
|
|
|
|
|
|
|
|
|
11 |
315,8 |
5061 |
52 |
0,9452 |
–0,0245 |
–0,0231 |
|
|
|
|
|
|
|
|
|
12 |
315,6 |
4812 |
57 |
0,8987 |
–0,0464 |
–0,0417 |
|
|
|
|
|
|
|
|
|
13 |
315,4 |
4569 |
62 |
0,8533 |
–0,0689 |
–0,0588 |
|
|
|
|
|
|
|
|
|
14 |
315,3 |
4453 |
67 |
0,8317 |
–0,0800 |
–0,0666 |
|
|
|
|
|
|
|
|
|
15 |
315,2 |
4342 |
71 |
0,8110 |
–0,0910 |
–0,0738 |
|
|
|
|
|
|
|
|
|
16 |
315,1 |
4230 |
76 |
0,7900 |
–0,1024 |
–0,0809 |
|
|
|
|
|
|
|
|
|
17 |
315,0 |
4118 |
81 |
0,7691 |
–0,1140 |
–0,0877 |
|
|
|
|
|
|
|
|
|
18 |
314,9 |
4015 |
86 |
0,7499 |
–0,1250 |
–0,0937 |
|
|
|
|
|
|
|
|
|
19 |
314,8 |
3905 |
90 |
0,7293 |
–0,1371 |
–0,1000 |
|
|
|
|
|
|
|
|
|
20 |
314,5 |
3603 |
95 |
0,6729 |
–0,1720 |
–0,1158 |
|
|
|
|
|
|
|
|
|
|
ИТОГО |
107084 |
|
|
–0,2202 |
0,2202 |
|
|
|
|
|
|
|
|
|
17
коэффициентов k должно быть приблизительно в 1,5 раза больше, чем максимальное из значений, определенных в табл. 4. Для рассматриваемого примера 1,4333 × 1,5 = 2,1499, на клетчатке вероятностей значение k изменяется в диапазоне от 0 до 2,0.
Точки с координатами (рэ,ki) наносят на клетчатку вероятностей (рис. 4, кривая 1), строят эмпирическую кривую распределения вероятностей максимальных расходов воды.
2 |
1 |
Рис. 4. Кривые распределения вероятностей превышения максимальных расходов воды:
1 – эмпирическая; 2 – аналитическая
Для сглаживания и экстраполяции эмпирических кривых распределения ежегодных вероятностей превышения максимальных расходов, как правило, применяется трехпараметрическое гамма-распределение.
Параметрами аналитической трехпараметрической кривой гаммараспределения являются:
• среднее многолетнее значение максимального расхода Q ;
18
•коэффициент вариации Cv;
•отношение коэффициента асимметрии к коэффициенту вариации
Cs/Cv. Параметры аналитической кривой гамма-распределения устанавливаются по гидрологическому ряду наблюдений ежегодных максимальных расходов методом приближенно наибольшего правдоподобия в зависимости от статистик λ2 и λ3, вычисляемых по формулам
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
∑lg ki |
|
|
|
|
|
λ |
2 |
= |
i=1 |
|
; |
|
(10) |
|
|
n |
−1 |
|
||||||
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
||
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
∑ki |
lg ki |
|
|
||
λ |
3 |
= |
i=1 |
|
|
. |
(11) |
||
|
−1 |
|
|||||||
|
|
|
|
n |
|
|
|
||
|
|
|
|
|
|
|
|
||
Значения ki, lgki и ki lgki рассчитываются в табличной форме (см., соот-
n
ветственно, графы 5, 6 и 7 в табл. 4). Для рассматриваемого примера ∑lg ki =
i=1
n
= –0,2202, а ∑ki lg ki = 0,2202 (см. итоговые значения граф 6 и 7 табл. 4).
i=1
Следовательно,
|
|
|
n |
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
∑lg ki |
|
−0,2202 |
|
|
|
∑ki |
lg ki |
|
0,2202 |
|
|
λ |
|
= |
i=1 |
= |
= −0,0116; λ |
|
= |
i=1 |
|
= |
= 0,0116. |
|||
2 |
|
3 |
|
|
||||||||||
n −1 |
20 −1 |
n |
−1 |
20 −1 |
||||||||||
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
Значение коэффициента вариации Cv и отношение Cs/Cv определяются по значениям статистик λ2 и λ3 по номограммам, приведенным в Приложении 1.
В рассматриваемом примере, учитывая значения λ2 и λ3, следует использовать номограмму, изображенную в приложении на рис.П.1.2, на которую наносится точка с координатами (λ3, λ2) (точка должна попасть на поле номограммы), по положению этой точки на номограмме определяются значения Cv и соотношение Cs/Cv.
Для рассматриваемого примера (рис. 5) Cv = 0,23, Cs = 3,0 Cv.
По значениям Cv и Cs/Cv по таблицам Приложения 2 определяются ординаты (модульные коэффициенты) аналитической кривой распределения ежегодных вероятностей превышения максимальных расходов воды. При необходимости значения, приведенные в Приложении 2, интерполируются. В рассматриваемом примере значения ординат получены интерполяцией значений, указанных в графах для Cv = 0,1 и Cv = 0,3 для случая Cs = 3,0 Cv.
19