Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

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

..pdf
Скачиваний:
6
Добавлен:
19.11.2023
Размер:
30.92 Mб
Скачать

расстояние между модой и центром распределения; расстояние между медианой и центром распределения равно приблизительно

Тг• Если радиус асимметрии г = 0 , то кривая распределения стано­

вится симметричной. В этом случае точки а, b, с совпадают.

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

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

уровней — весьма часто отрицательную.

вполне

Кривая обеспеченности выражается уравнением, которое

определяется тремя параметрами:

 

1) средним арифметическим значением ряда Q0,

 

2) коэффициентом вариации (изменчивости) ряда Cv,

 

3) коэффициентом асимметрии ряда С*.

 

Рассмотрим эти параметры.

 

Среднее арифметическое значение ряда Q0 находится

по фор­

муле

 

Qo —

(IV-1)

где EQi— сумма всех максимальных расходов данного ряда; п — число членов ряда.

Одно среднее арифметическое значение Q0 не может достаточно

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

на какое-то число, а вторую половину членов уменьшим на это же число, то величина Q0 не изменится, тогда как характер распреде­

ления максимальных расходов, безусловно, изменится.

Для оценки изменчивости ряда определяют среднее квадратиче­ ское отклонение

о =

£ (Qi Q o)a

(IV-2)

п 1

 

 

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

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

Итак, коэффициент вариации

 

г

— _ L- —

1

л /

£ (Qi— QaV

 

 

 

 

v ~

Qa ~~

Qa

У

п -

1

 

 

Но Qi KlQo-

 

 

 

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

.

=

л / U K j Q o - W

_ 1 /

Е (Kj-

1)а

(IV-3)

*

Q0

У

п — 1

 

У

п — 1

 

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

Коэффициент вариации Cv не полностью характеризует измен­

чивость ряда. В самом деле, два ряда с одинаковыми значениями Q0 и Cv могут иметь различное число отклонений Q,— Q0 в поло­

жительную и отрицательную сторону, т. е. разную асимметричность. Для оценки асимметричности ряда необходимо определить

коэффициент асимметрии, который находится по формуле

Е (Ki - 1)3

(IV-4)

("■ - 1 ) c l

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

Для точного вычисления коэффициента асимметрии по формуле (IV-4) необходимо иметь весьма продолжительный ряд наблюдений, содержащий 100 и более членов. Такой ряд на практике в боль­ шинстве случаев отсутствует. Поэтому коэффициент асимметрии приходится определять по приближенным формулам, приведенным в § 16.

Построение кривой обеспеченности модульных коэффициентов по трем параметрам Q0, С,0 и Cs производится с помощью составлен­

ной С. И. Рыбкиным специальной таблицы, в которой приведены отклонения ординат кривой обеспеченности от середины, т. е. от /С/=1 при Cv— 1,0, С /= 0 -КЗ и при разной обеспеченности. Эти

отклонения обозначаются через Ф*. Таблица С. И. Рыбкина при­ ведена в приложении За.

Для получения ординат кривой обеспеченности модульных коэф­ фициентов нужно взятую из таблицы величину Ф* умножить на коэффициент вариации Cv (так как в таблице отклонения Ф1 да-

ются для С „=1,0) и к полученной величине прибавить единицу

(так как в таблице отклонения Фг даются от /Сг=1) — см. рис. 44. Таким образом, ординаты кривой обеспеченности подсчитыва­

ются по формуле

 

*1 = ^ + 1 .

(IV-5)

В гидрологических расчетах наряду с методом математической статистики применяется теория корреляции, которая является частью теории вероятностей. Теория корреляции используется для

а)

Ф

Ф

Рис. 45. Возможные случаи связи между величинами

Уи X:

а— функциональная зависимость У от X; б—отсутствие связи между У и X; а — коррелятивная связь между У и X

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

Поясним сущность метода корреляции.

Пусть мы имеем два ряда изменяющихся гидрологических вели­ чин X и Y, причем каждому значению X соответствует определен­ ное значение У:

*1, х * х 3,

к*.

Если в системе координат X Y нанести точки, определяемые координатами Х г и Y v Х г и У2, Х 3 и У3 и т. д., то получится неко­

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

1. Точки располагаются таким образом, что по ним можно построить некоторую кривую (рис. 45, а). Если подобрать уравне­ ние этой кривой, то получится зависимость между У и X в матема­ тической форме, т. е. Y —f{X) . Следовательно, в данном случае между величинами Y и X имеет место функциональная зависимость.

2.Точки рассеяны по всему полю (рис. 45, б). В этом случае между величинами Y и X никакой связи не существует.

3.Точки хотя и рассеяны по полю и не Дают функциональной зависимости между величинами У и X, но разбросаны не так силь-

но, как во втором случае (рис. 45, в). Они располагаются таким

образом, что видна некоторая зависимость между величинами Y и X . Эта зависимость тем теснее, чем ближе находятся точки от некоторой средней линии ab. В данном случае между величинами Y и X имеет место коррелятивная связь.

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

 

г =

 

Е

-

Х0) (К, -

Y0)

,

 

 

 

 

(IV-6)

 

- ------

 

 

 

 

Уо)2

 

 

 

 

 

 

V Е ( * i - *о)2 Е (Yi -

 

 

 

 

 

 

 

где Х 0 и Y 0— средние

арифметические

значения

рядов величин

X

и

Y

 

 

 

 

 

 

 

 

 

± 1. Если г О,

Коэффициент корреляции г изменяется от 0 до

то между величинами

Y

и X

никакой связи

нет.Если

r = ± 1, то

 

 

 

между

величинами

Y

и X существует

 

 

 

функциональная

связь.

В

 

остальных

 

 

 

случаях

связь

 

коррелятивная,

причем

 

 

 

при

—0 ,6 <

г < + 0 ,6

эта

связь

очень

 

 

 

слабая.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поэтому

практически

считается, что

 

 

 

коррелятивная

связь

между величина­

 

 

 

ми

У и X

существует

при

значениях

 

 

 

коэффициента

корреляции

г > + 0 ,6 или

 

 

 

г < —0,6, причем чем

ближе

коэффици­

 

 

 

ент корреляции

к ± 1, тем связь теснее.

 

 

 

Если

с увеличением X

увеличивает­

 

 

 

ся

К,

то коэффициент

корреляции

г

Рис. 46. Прямые регрессии:

имеет

положительные

значения.

Если

I — прямая регрессии

X по

Y;

же с увеличением X

величина Y

умень­

Л — прямая регрессии К по X

шается,

то

коэффициент

корреляции

г

 

 

 

имеет отрицательные

значения.

 

 

Широкое применение в теории корреляции имеют так называе­ мые прямые регрессии. Различают две прямые регрессии: прямую регрессии Y по X и прямую регрессии X по Y (рис. 46). Прямой регрессии Y по X называется прямая, проведенная по нанесенным

на графике точкам так, чтобы сумма квадратов отклонений от нее ординат Y отдельных точек была бы наименьшей. Прямой регрес­ сии X no Y называется прямая, проведенная по нанесенным на

графике точкам так, чтобы сумма квадратов отклонений от нее

абсцисс X отдельных точек была бы наименьшей.

 

Уравнение прямой

регрессии Y по

X

имеет

следующий вид:

Y

- Y a= r ^ - ( X

- X

о).

(IV-7)

QX

Уравнение прямой регрессии X по У

имеет следующий вид:

X - X 0= r ^ . ( Y - Y

0).

(IV-8)

В этих уравнениях оЛ. и оу — средние квадратичные

отклоне­

ния X и У от средних величин Х 0 и Y 0.

 

 

Величины ох и оу определяются по формулам

 

-

л /

Ъ(Хг-Хо)а .

Л' ~

У

п — 1

'

„ _

I f

ЫГ1-Уо)1

аУ-

V

„ _ 1

■»

(IV-9)

(IV-10)

где п — число членов ряда.

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

_ £ (* < -* ,) (Г ,-Г ,)

(/Z 1) Cj. Zy

(IV-11)

Прямые регрессии Y по X и X по Y пересекаются в точке с координатами Х = Х 0 и К=К„ (рис. 46). Чем ближе коэффициент

корреляции г к ± 1,

тем меньше угол а, образованный этими пря­

мыми. При г 0,

т.

е. в случае функциональной связи» угол а=0

и прямые регрессии

сливаются.

 

 

Вероятная ошибка в определении коэффициента корреляции

находится по формуле

 

 

 

 

£ = ± 0,674 1 —г2

 

 

 

Vn

'

 

Предельная ошибка в определении коэффициента корреляции

принимается равной

£ тах= 4 £ .

 

 

Вероятные значения коэффициента

корреляции заключаются

в пределах г ± Е ,

а

предельные значения — в

пределах r ± £ mox.

Коррелятивная связь между величинами X

и Y считается тес­

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

1) коэффициент корреляции г и все его вероятные и предельные

значения больше

+ 0,6 или меньше —0,6;

 

2)величина г ± Е тйХ сохраняет в обоих случаях знак коэффи­

циента корреляции;

3)предельная ошибка £ гаах значительно меньше коэффициента корреляции.

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

Вэтом отношении все реки можно подразделить на две группы.

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

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

§ 16. Расчеты максимальных расходов и уровней воды хорошо изученных рек

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

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

После составления ряда наивысших за каждый год уровней воды по кривой расхода Q = f (Н) находятся соответствующие этим

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

Затем определяются среднее арифметическое значение ряда максимальных расходов Q0 по формуле (IV-1), коэффициент вариа­ ции Cv по формуле (IV-3) и коэффициент асимметрии Cs. Общая формула для коэффициента асимметрии Cs (IV-4) при сравнительно

непродолжительном ряде наблюдений, обычно имеющем место на практике, дает очень большие ошибки. Например, при числе лет

наблюдений л = 20

ошибка составляет 55%;

даже

при я=100 лет

эта ошибка равна

24,5%. Поэтому

коэффициент

асимметрии Cs

обычно определяют по другим формулам.

С.

 

М. Ф. Менкеля

с , = -Г-2^ — ,

(IV-12)

Широкое применение получила

формула

Н. Крицкого и

*—Amin

где К т\п— наименьший модульный коэффициент ряда;

К min

__ Qmln

Qo

 

ГОСТ СН 2—57 рекомендует применять следующие соотношения между коэффициентами Cs и Cv:

для снеговых паводков

С ,= 2С„;

(IV-13)

для дождевых паводков

 

CS = Z + 4C9.

(IV-14)

После установления величин Q0, Cv и Cs определяются ординаты

кривой обеспеченности по формуле (IV-5), и по полученным данным строится кривая обеспеченности.

Для проверки построенной теоретической кривой обеспечен­ ности модульных коэффициентов нужно нанести на полученный график имеющиеся данные фактических наблюдений. Для этого необходимо определить эмпирическую вероятность превышения в процентах фактических модульных коэффициентов по формуле Н. Н. Чегодаева:

 

? =

(IV-16)

где р — эмпирическая

вероятность превышения каждого

члена

ряда модульных коэффициентов в процентах;

 

п — число членов

ряда модульных коэффициентов;

 

т — порядковый номер члена данного ряда при расположе­

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

Данные, полученные по формуле Н. Н. Чегодаева, наносятся на график (на рис. 44 показаны кружочками). Если эти данные располагаются близко от теоретической кривой обеспеченности, то это означает, что параметр кривой обеспеченности С5 выбран

правильно. Если же данные фактических наблюдений располагают­

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

то в этом слу­

чае необходимо изменить величину коэффициента

асимметрии Cs.

При известных параметрах Q0, Cv и Cs максимальный

расход

заданной вероятности превышения находится по формуле

 

С =

«о<<рС„+1).

 

(IV-16)

Величина Ф, входящая

в эту формулу, берется из

таблицы

С. И. Рыбкина (см. приложение За) в зависимости от принятого значения коэффициента асимметрии Cs и заданной вероятности

превышения паводка.

Для построения кривой обеспеченности необходимо иметь ряд, содержащий не менее 15 максимальных годовых расходов, из кото­ рых 1—2 расхода должны относиться к многоводным годам. Необ­ ходимо также иметь надежную кривую расхода Q = f (Н), построен­

ную по данным гидрометрических измерений. Кривая расхода имеет

характерное петлеобразное очертание (рис. 47). Это объясняется тем, что при движении паводочной волны уклон водной поверх­ ности на лобовой ее части больше, чем на тыльной (см. рис. 58). Следовательно, при одной и той же отметке Я (рис. 58) уклон вод­ ной поверхности в период подъема паводка будет больше, чем

Рис. 47. Кривая расхода Q=f (Я):

Qn — расход на подъеме паводка; Qc — расход на спаде паводка

в период его спада. В результате этого при постоянном значении отметки Я скорости течения, а значит, и расход будут иметь боль­ шие значения на подъеме паводка и меньшие на спаде его. Это и приводит к петлеобразному очертанию кривой расхода Q = f(H ).

Рис. 48. Графическая экстраполяция кривой расхода до ГВВ с контролем по зависимости

Фгвв ~ шгвв °гвв

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

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

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

(рис.

49). При

наличии

н

 

функциональной

или кор-

 

релятивной

связи между

 

 

максимальными уровнями

 

 

на

этих

водомерных по­

 

 

стах нужно по кривой свя­

 

 

зи

максимальных уровней

 

 

найти

все

недостающие

 

 

уровни, а затем по кривой

 

 

расхода —

соответствую­

 

 

щие этим уровням макси­

 

 

мальные

расходы.

 

 

При определении макси­

 

 

мальных

расходов задан­

Рк

нб

ной вероятности превыше-

49. Кривая связи годовых макси­

ния

иногда

встречаются

 

 

мальных уровней

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

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

Предположим, что имеется непрерывный ряд годовых макси­ мальных расходов за п лет Qlf Q2> Q& > Qn и один расход редкой повторяемости QH, который является наибольшим за N лет, при­ чем N^>n. В этом случае параметры Q0 и Cv кривой обеспеченности

расходов определяются по формулам С. Н. Крицкого и М. Ф. Менкеля:

(IV-17)

(IV-18)

где K t и K N— модульные коэффициенты, которые находятся по

формулам

Коэффициент асимметрии

Cs

определяется в этом случае по

обычным формулам, приведенным

выше.

После определения параметров

Q0, Cv и Cs максимальный рас­

ход заданной

вероятности

превышения находится по формуле

(IV-16).

 

 

 

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

максимального

расхода заданной

вероятности превышения, когда

От

Ш

I

.5 W

20 X 40 50 SO 70

80

90 95

99

99.9

99.39

 

 

 

 

05еспечг,тост.

%

 

 

 

 

Рис. 50.

Клетчатка вероятностей нормального

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

равно­

 

 

 

мерной

вертикальной

шкалой

 

 

 

этот расход находится непосредственно по формуле (IV-16), в практике проектирования мостовых переходов большое распро­ странение получил графо-аналитический способ определения ука­ занного расхода. По этому способу кривая обеспеченности, построен­ ная по данным фактических наблюдений, графически экстраполи­ руется до расчетного ‘значения вероятности превышения паводка. Если кривую обеспеченности построить в обычных координатах (см. рис. 44), то верхняя и нижняя ветви кривой будут иметь очень крутые подъемы, в результате чего надежная графическая экстра­ поляция кривой обеспеченности до расчетного значения вероят­ ности превышения паводка невозможна. Поэтому кривая обеспе­ ченности строится не в обычных координатах, а в так называемой клетчатке вероятностей, которая позволяет значительно умень­ шить кривизну кривой обеспеченности. Наибольшее распростране­ ние получила клетчатка нормального (симметричного) распределе­ ния. Неравномерная горизонтальная шкала этой клетчатки стро-

Соседние файлы в папке книги