До початку сторінка 12 з 14

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

Завантажити документ

8 W9=SYS(7, 0)

9 OPEN «LP:» FOR OUTPUT AS FILE #1

10 PRINT #1, «ICE3»\ PRINT #1,

15 DIM F(30), G1(30), F2(30)

22 DIM O0(30), O1(30), O2(30), O3(30), O(30), G2(30)

23 DIM X(30), Y(30), U(30), V(30), V1(30)

24 DIM Y1(30), Y2(30)

25 READ N, G0

26 DATA 25, 1961

27 REM ВВОД ПРОДОЛЖИТЕЛЬНОСТИ МОРОЗНОГО ПЕРИОДА КАЖДОГО ГОДА

60 FOR I=1 ТО N

70 READ F(I)

364 REM ФОРМИРОВАНИЕ РЯДА ЛЕТ НАБЛЮДЕНИЙ ПРОДОЛЖИТЕЛЬНОСТИ МОРОЗНОГО ПЕРИОДА

366 G1(I)=G0+1

390 NEXT I

399 REM РАНЖИРОВАНИЕ РЯДА ПРОДОЛЖИТЕЛЬНОСТИ МОРОЗНОГО ПЕРИОДА

400 К=0

402 FOR I=1 ТО N-1

404 IF F(I)>F(I+1) THEN 420

406 R1=F(I+1)

408 R2=G1(I+1)

410 F(I+1)F(I)

412 G1(I+1)=G1(I)

414 F(I)=R1

416 G1(I)=R2

418 K=1

420 NEXT I

422 IF K=1 THEN 400

449 REM

450 DATA 226, 229, 226, 230, 216

451 DATA 219, 225, 228, 214, 222

452 DATA 222, 227, 225, 221, 223

453 DATA 232, 232, 220, 232, 230

454 DATA 226, 229, 226, 229, 226

1000 REM РАСЧЕТ СУММ ЖИДКИХ ОСАДКОВ

1043 REM ВВОД КОЛ-ВА ЖИДКИХ ОСАДКОВ ЗА ИЮЛЬ КАЖДОГО ГОДА

1044 FOR I=1 ТО N

1046 READ O0(I)

1048 NEXT I

1049 REM ВВОД КОЛ-ВА ЖИДКИХ ОСАДКОВ ЗА АВГУСТ КАЖДОГО ГОДА

1050 FOR I=1 ТО N

1052 READ O1(I)

1054 NEXT I

1055 REM ВВОД КОЛ-ВА ЖИДКИХ ОСАДКОВ ЗА СЕНТЯБРЬ КАЖДОГО ГОДА

1056 FOR I=1 ТО N

1058 READ O2(I)

1060 NEXT I

1061 REM ВВОД КОЛ-ВА ЖИДКИХ ОСАДКОВ ЗА ОКТЯБРЬ КАЖДОГО ГОДА

1062 FOR I=1 ТО N

1064 READ O3(I)

1066 NEXT I

1499 REM РАСЧЕТ СУММ ЖИДКИХ ОСАДКОВ ЗА 3 МЕСЯЦА

1500 FOR I=1 TO N

1510 O(I)=0

1520 NEXT I

1630 FOR I=1 TO N

1540 O(I) =O0(I)+O1(I)+O2(I)+O3(I)

1560 NEXT I

1609 REM ФОРМИРОВАНИЕ РЯДА ЛЕТ НАБЛЮДЕНИЙ ЖИДКИХ ОСАДКОВ

1610 FOR I=1 TO N

1620 G2(I)=G0-1+I

1630 NEXT I

1639 REM РАНЖИРОВАНИЕ РЯДА СУММ ЖИДКИХ ОСАДКОВ

1640 K=0

1650 FOR I=1 TO N-1

1660 IF O(I)>O(I+1) THEN 1740

1670 R1=O(I+1)

1680 R2=G2(I+1)

1690 O(I+1)=O(I)

1700 G2(I+1)=G2(I)

1710 O(I)=R1

1720 G2(I)=R2

1730 K=1

1740 NEXT I

1750 IF K=1 THEN 1640

1760 FOR I=1 TO N

1764 REM ВЫЧИСЛЕНИЕ ЭМПИРИЧЕСКОЙ ВЕРОЯТНОСТИ

1765 F2(I)=I/(N+1)

1780 NEXT I

1949 REM ЖИДКИЕ ОСАДКИ ЗА ИЮЛЬ ПО КАЖДОМУ ГОДУ НАБЛЮДЕНИЙ

1950 DATA 83.1, 53.8, 50.0, 45.1, 5.7

1951 DATA 63.0, 68.8, 48.3, 90.7, 58.3

1952 DATA 23.9, 68.4, 17.4, 48.8, 60.5

1953 DATA 39.1, 54.3, 99.4, 86.5, 67.7

1954 DATA 62.7, 34.9, 94.3, 93.8, 163.2

1959 REM ЖИДКИЕ ОСАДКИ ЗА АВГУСТ ПО КАЖДОМУ ГОДУ НАБЛЮДЕНИЙ

1960 DATA 37.6, 61.9, 20.6, 75.2, 71.9

1961 DATA 90.3, 28.6, 12.5, 16.9, 94.9

1962 DATA 14.3, 75.6, 17.6, 33.9, 96.7

1963 DATA 34.3, 66.2, 37.8, 91.6, 33.7

1964 DATA 24.3, 111.4, 103.0, 84.6, 64.7

1979 REM ЖИДКИЕ ОСАДКИ ЗА СЕНТЯБРЬ ПО КАЖДОМУ ГОДУ НАБЛЮДЕНИЙ

1980 DATA 80.4, 0.4, 16.6, 12.1, 8.8

1981 DATA 25.9, 61.6, 42.2, 30.3, 32.5

1982 DATA 29.8, 19.0, 42.3, 76.3, 64.5

1983 DATA 1.3, 49.1, 2.5, 28.9, 42.4

1984 DATA 38.9, 47.7, 51.4, 12.9, 28.1

1998 REM ЖИДКИЕ ОСАДКИ ЗА ОКТЯБРЬ ПО КАЖДОМУ ГОДУ НАБЛЮДЕНИЙ

1999 DATA 0, 0, 10.7, 0, 0

2000 DATA 0, 0, 0, 0, 0

2001 DATA 0, 0, 30.1, 2.6, 0

2002 DATA 0, 0, 0, 0, 0

2003 DATA 16.9, 76.4, 0, 19.6, .5

2100 REM АППРОКСИМАЦИЯ КРИВЫХ

2110 FOR I=1 TO N

2120 X(I)=F2(I)

2130 Y(I)=O(I)

2140 NEXT I

2145 PRINT 1, «A, U - КОЭФФИЦИЕНТЫ УРАВНЕНИЯ КРИВОЙ Y=1/(A*X+B) СУМ. Ж. ОСАД.»

2146 PRINT #1. «Е – ТОЧНОСТЬ - УСРЕДНЕННОЕ СРЕДНЕКВАДРАТИЧНОЕ ОТКЛОНЕНИЕ

2150 GOSUB 3000

2160 А1=А9

2170 В1=В9

2180 Е1=Е

2190 FOR I=1 TO N

2200 Y(I)=F(I)

2210 NEXT I

2215 PRINT #1, «А, В - КОЭЭФИЦИЕНТЫ УРАВНЕНИЯ КРИВОЙ Y=1/(A*X+B) ПРОД. М. П.»

2216 PRINT #1, «Е - ТОЧНОСТЬ - УСРЕДНЕННОЕ СРЕДНЕКВАДРАТИЧНОЕ ОТКЛОНЕНИЕ

2220 GOSUB 3000

2230 А2=А9

2240 В2=В9

2250 Е2=Е

2255 REM ПЕЧАТЬ ЗАГОЛОВКА ТАБЛИЦЫ РЕЗУЛЬТАТОВ

2260 PRINT #1, «Р ЭМП, «СУМ. Ж. ОСАД.», «СУМ. Ж. ОСАД.», «ГОД», «ПРОД. М.П.», «ПРОД. М.П.», «ГОД»

2265 PRINT #1, « », «НАБЛ.». «АППР.», « », «РАСЧ.», «АППР.»

2270 FOR I=1 ТО N

2280 Y1(I)=1/(A1*F2(I)+B1)

2290 Y2(I)=1/(A2*F2(I)+B2)

2295 REM ПЕЧАТЬ ТАБЛИЦЫ РЕЗУЛЬТАТОВ

2300 PRINT #1, F2(I), O(I), YI(I), G2(I), F(I), Y2(I), G1(I)

2310 NEXT I

2315 CLOSE #1

2350 STOP

3000 REM ПОДПРОГРАММА АППРОКСИМАЦИИ КРИВОЙ ПО МЕТОДУ НАИМЕНЬШИХ КВАДРАТОВ

3020 S=0

3021 S1=0

3022 S2=0

3023 S3=0

3024 S4=0

3070 FOR I=1 TO N

3075 U(I)=X(I)

3080 V(I)=1/Y(I)

3085 NEXT I

3195 FOR I=1 TO N

3170 S1=S1+U(I)*U(I)

3175 S2=S2+U(I)

3180 S3=S3+U(1)*V(I)

3185 S4=S4+V(I)

3190 NEXT I

3195 D=S1*N-S2*S2

3200 D1=S3*N-S4*S2

3205 D2=S1*S4-S2*S3

3210 A9=D1/D

3211 B9=D2/D

3215 FOR I=1 TO N

3220 V1(I)=A9*U(I)+B9

3225 S=S+(V1(I)-V(I))*(V1(I)-V(I))

3230 NEXT I

3235 E=SQR(S/N)

3255 PRINT #1, «A=»; A9, «B=»; B9,«E=»; E

3260 RETURN

3300 END

ПРИЛОЖЕНИЕ 7 Пример прогнозирования расчетных параметров наледи ключевых вод на реке Д Горная река находится в I дорожно-климатической зоне. Ширина долины 100-140 м, с обеих сторон ограничена крутыми склонами. Русло сложено гравийно-галечниковыми отложениями, которые на глубине 2,4 м подстилаются глинистыми сланцами. На высоте 2,5 м от поверхности русла из разлома коренных пород правого склона в течение всего года изливается источник, расход которого год от года изменяется в зависимости от количества осадков, выпадающих в летне-осеннее время. При наступлении морозного периода источник образует наледь.

Обследование наледи производилось в конце апреля 1981 г. Ее объем составлял 7140 м³, площадь 14600 м, средняя толщина 0,49 м, максимальная толщина 0,93 м.

На водотоке проектируется четырехпролетный железобетонный мост отверстием 80 м на столбчатых опорах. Дорога IV категории.

Определение расчетных параметров ключевой наледи, согласно п. 1.25 СНиП 2.05.03-84, производим с вероятностью превышения 2 %.

Данные ближайшей метеорологической станции о жидких осадках и продолжительности морозного периода за 25 лет (табл. 1) вводим в программу Лед-III. ЭВМ выдает расположенные в ранжированные ряды ординаты аппроксимирующих кривых обеспеченности сумм жидких осадков и продолжительности морозного периода, которые приводятся в табл. 2. На клетчатке вероятностей строим кривые обеспеченности сумм жидких осадков и продолжительности морозного периода (рис. 1).

Таблица 1. Суммы жидких осадков и продолжительности морозного периода в районе мостового перехода через реку Д

Годы

Суммы жидких осадков

Морозный период

Июль

Август

Сентябрь

Октябрь

Переходы температуры

Количество морозных дней

Осень

Весна

1961

83,1

37,6

80,4

28,09

1962

53,8

61,9

0,4

29,09

11,05

226

1963

50,0

20,6

16,6

10,7

28,09

14,05

226

1964

45,1

75,2

12,1

28,09

11,05

230

1965

5,7

71,9

8,8

5,10

14,05

216

1966

63,0

90,3

25,9

26,09

8,05

219

1967

68,8

28,6

61,6

1,10

2,05

225

1968

48,3

12,5

42,2

l,10

13,05

228

1969

90,7

16,9

30,3

26,09

15,05

214

1970

58,3

94,9

32,5

2,10

27,04

222

1971

23,9

14,3

29,8

21,09

11,05

229

1972

68,4

75,6

19,0

16,09

30,04

232

1973

17,4

17,6

42,3

30,1

29,09

29,04

227

1974

48,8

33,9

76,3

2,6

29,09

11,05

225

1975

60,5

96,7

64,5

3,10

7,05

221

1976

39,1

34,3

1,3

21,09

13,05

223

1977

54,3

66,2

49,1

22,09

9,05

232

1978

99,4

37,8

2,5

26,09

11,05

232

1979

86,5

91,6

28,9

22,09

3,05

220

1980

67,7

33,7

42,4

28,09

11,05

232

1981

62,7

24,3

38,9

16,9

28,09

14,05

230

1982

34,9

111,4

47,7

76,4

28,09

11,05

226

1983

94,3

103,0

51,4

28,09

14,05

229

1984

93,8

81,6

12,9

19,6

28,09

11,05

226

1985

163,2

61,7

28,1

5,0

28,09

14,05

229

11,05

226

Рис. 1. Кривые обеспеченности сумм жидких осадков и продолжительности морозного периода на реке Д (питание ключевыми водами)

Таблица 2

Ординаты аппроксимирующих кривых обеспеченности суммы жидких осадков и продолжительности морозного периода на водотоке Д

ICE3

А, Б - КОЭФФИЦИЕНТЫ УРАВНЕНИЯ КРИВОЙ Y=1/(A*X+B) СУМ. Ж. ОСАД.

Е - ТОЧНОСТЬ - УСРЕДНЕННОЕ СРЕДНЕКВАДРАТИЧНОЕ ОТКЛОНЕНИЕ

А=9.60166Е-03 В=2.38243Е-3 Е=1.02441Е-03

А, В - КОЭФФИЦИЕНТЫ УРАВНЕНИЯ КРИВОЙ Y=1/(A*X+B) ПРОД. М. П.

Е - ТОЧНОСТЬ - УСРЕДНЕННОЕ СРЕДНЕКВАДРАТИЧНОЕ ОТКЛОНЕНИЕ

А=3.26907Е-04

Р ЭМП

СУМ. Ж. ОСАД. НАБЛ.

В=4.27509Е-03

СУМ. Ж. ОСАД. АППР.

ГОД

Е=2.68221Е-05

ПРОД. М. П. НАБЛ.

ПРОД. М. П. АППР.

ГОД

.0384615

270.4

363.409

1982

932

233.227

1977

.0769231

256.5

320.408

1985

232

232.545

1978

.115385

248.7

286.507

1983

232

231.867

1980

.153846

221.7

259.094

1975

230

231.193

1965

.192308

210.9

236.468

1984

230

230.523

1981

.230769

207

217.477

1979

229

229.857

1963

.269231

201.1

201.309

1961

229

229.195

1983

.307692

185.7

187.379

1970

229

228.536

1985

.346154

179.2

175.252

1966

228

227.881

1969

.384615

169.6

164.599

1977

227

227.23

1973

.123077

163

155.167

1972

226

226.583

1962

.461538

161.6

146.757

1974

226

225.938

1964

159

139.213

1967

226

225.299

1982

.538462

143.8

132.406

1980

226

224.663

1984

.576923

142.8

126.233

1981

226

224.03

1986

.615385

139.7

120.611

1978

225

223.401

1968

.653846

137.9

115.468

1969

225

222.775

1974

.692308

132.4

110.745

1964

223

222.153

1976

.730769

116.1

106.394

1962

222

221.534

1971

.769231

107.4

102.372

1973

222

220.919

1972

.807692

103

98.6425

1968

221

220.307

1975

.846154

97.9

95.1755

1963

220

219.698

1979

.884615

86.4

91.9439

1965

219

219.093

1967

.923077

74.7

88.9245

1976

216

218.491

1966

.961538

86.0971

1971

214

217.892

1970

Из табл. 2 в год обследования наледи (1981) продолжительность морозного периода была 231 день с эмпирической вероятностью 0,192. В предыдущий год (1980) сумма жидких осадков за четыре летне-осенних месяца составляла 132 мм с эмпирической вероятностью 0,538.