Изчисляване на котата на водната повърхност с цел определяне на потенциално опасни участъци
Определянето на потенциално застрашените от наводнения райони е извършено съгласно препоръчаните методи от Методиката за ПОРН.
В тези указания са предложени два основни критерия за дефиниране на потенциално опасните участъци, а именно:
-
хоризонтално отстояние между разглеждания обект и преминаващата покрай него река и
-
вертикално отстояние (дълбочина на водата), т.е. определяне на дълбочина на заливане и нейното съпоставяне с цифровия модел на терена.
Във втория случай може да бъде определена една представителна дълбочина (кота водно ниво) за границите на участъка или да бъдат изчислени различни дълбочини (коти на водните нива по протежение на участъка).
Очевидно вторият критерий, както това е посочено и в самата методика, притежава значителни предимства и е избран като водещ при определяне на потенциално застрашените от заливане участъци в ЗБР за БУ.
В този случай са възможни различни подходи в зависимост от наличната топографска и хидроложка информация:
-
Изчисляване на нормалната дълбочина за условията на стационарно, равномерно движение.
-
Хидравлично моделиране с използване на 1D модели, при стационарен и нестационарен режим на движение на водата.
-
Хидравлично моделиране с използване на 1D+2D модели.
Всеки от цитираните подходи има своите предимства и недостатъци, които са подробно анализирани в рамките на самата методика. От чисто формална гледна точка най-точни са 2D моделите, които освен котите на водната повърхност дават и векторите на скоростта в изчислителните възли на мрежата, т.е. при тях освен котата на водната повърхност може да бъде добавен и критерият за скорост на водното течение. Изработването на такива модели изисква изключително голям обем прецизна и скъпоструваща информация и към тях трябва да се прибягва след съответна обосновка. Очевиден е факта, че използването на 2D моделите само по себе си не може да компенсира липсата на прецизен модел на терена за речното корито и прилежащите заливаеми територии. Следва извода, че използването на 1D и евентуално 2D моделите в условията на стационарно неравномерно движение или нестационарно движение на водата ще бъде приоритетно на следващия етап от изготвянето на Плана за управление на риска от наводнения, когато ще бъде изготвено подробно заснемане на речното корито и прилежащите заливаеми низини, заедно с цялата прилежаща инфраструктура, мостове, речни прагове, водовземни и заустващи съоръжения, защитни диги и др.
Имайки предвид времевата рамка за изготвяне на ПОРН, за предварителните изследвания, като единствено възможен е избран метода с определяне на нормалната дълбочина за условията на стационарно, равномерно движение.
От хидравлична гледна точка този метод е приложим, когато геометричните и хидравличните параметри на течението не се променят в границите на един достатъчно дълъг участък. В този случай се използва формулата на Вези:
средна скорост на водното течение
водно количество
- скоростен множител (формула на Майнинг)
- хидравличен радиус
- площ на намокреното сечение
- наклон на водната повърхност в участъка (съвпада с наклона на дъното и хидравличния наклон)
В много от случаите при наличие на подприщване или наличие на съоръжения точността на метода се влошава, но същото важи и за останалите методи, ако не се разполага с достатъчно точна информация за критичните участъци на подприщване (мостове, прагове, водовземни и др.).
По-надолу е описана последователността на извършване на хидравличните изчисления:
-
Подготовка на входни данни
Геометрични данни - за изчислителните профили. За всеки изчислителен участък са подбрани след анализ на картен материал в мащаб 1:5000, ортофотоснимки и наличните данни за минали наводнения между три и пет изчислителни профила (фигура 1). Въз основа на това е определен средния надлъжен наклон на водната повърхност в участъка. Данните за профилите включват двойки числа: разстояние от постоянно начало и кота дъно (абсолютна или относителна);
Фигура 1. Разположение на изчислителните профили.
Избор на коефициент на грапавина, след анализ на морфологичните характеристики и зърнометрията на дъното в участъка. При хидравличните изчисления са използвани стойности на коефициента на грапавината във формулата на Манинг между n=0.035 и 0,040 за зоната на основното речно легло и n=0.045 и 0,050 за зоната на речната тераса извън основното речно корито, като тези стойности са оценени въз основа на оглед на съществуващото положение “in situ” и сравнителен анализ с литературни данни и резултати аналогични изследвания. (Guide for Selecting Manning's Roughness Coefficients for Natural Channels and Flood Plains United States Geological Survey Water-supply Paper 2339, G.J. Arcement, Jr. and V.R. Schneider, USGS).
Коефициентът на Шези (скоростния множител) е определян по формулата на Манинг. Това изисква определяне на площите на напречното сечение, намокрения периметър и хидравличния радиус за напречния профил, което при неправилна форма на профила изисква значителна изчислителна работа.
-
Провеждане на хидравлични изчисления
По принцип изчисляването на дълбочините (котите на водната повърхност) в условията на равномерно движение е проста задача, която при неправилна форма на напречното сечение се усложнява единствено от необходимостта за определяне на площите на напречното сечение, намокрения периметър и хидравличния радиус при различни коти на водното ниво. Изчисленията могат да се извършат ръчно или с помощта на макроси в EXCEL.
В настоящото изследване изчисленията са извършени с компютърния модел HEC-RAS (Hydrologic Engineering Center - River Analysis System), разработен от корпуса на военните инженери на САЩ (U.S. Army Corps of Engineers). Независимо, че същият се използва главно за провеждане на изчисления в условията на неравномерно движение, той може да бъде пригоден и за случая на равномерно движение, като граничните условия се задават чрез нормалните дълбочини съответстващи на средния надлъжен наклон за изчислителния участък.
Получените резултати включват изчислената кота на водното ниво в профила за оразмерителното водно количество и серия от водни нива съответстващи на различни водни количества в диапазона на възможните такива (ключова крива на профила).
Въз основа на определената кота на водното ниво в изчислителните профили са определени площите на заливане, като през тях е прекарана равнина с надлъжен наклон равен на измерения за представителния участък и е определена пресечницата на тази равнина с модела на терена.
Таблица с оразмерителни водни количества за определяне на петната на заливане при районите с потенциален риск от наводнения
D
|
река
|
Населено място
|
X север
|
Y изток
|
Fi
|
La
|
AF
|
Q1%
|
1
|
Места
|
с. Хаджи Димово
|
4602392,475
|
237541,0821
|
41°31'48.6
|
23°51'14.8
|
2264,36
|
693,60
|
2
|
Кaнинa
|
с. Огняново
|
4612352,907
|
233038,7271
|
41°37'05.7
|
23°47'44.9
|
227,69
|
143,38
|
3
|
Златарица
|
гр. Елешница
|
4640424,324
|
219592,6906
|
41°51'57.8
|
23°37'17.2
|
103,82
|
41,94
|
4
|
Глазне
|
с. Баня
|
4643183,872
|
211376,729
|
41°53'16.4
|
23°31'16.6
|
85,49
|
72,51
|
5
|
Бела река
|
гр. Разлог
|
4642760,479
|
205092,8844
|
41°52'54.4
|
23°26'45.3
|
63,89
|
26,96
|
6
|
Глазне
|
гр. Банско
|
4637200,725
|
207730,2041
|
41°49'58.0
|
23°28'49.3
|
75,40
|
65,26
|
7
|
Места
|
гр. Якоруда
|
4658257,194
|
225603,1692
|
42°01'42.6
|
23°41'07.6
|
264,85
|
98,37
|
8
|
Струма
|
граница
|
4586536,632
|
194074,9772
|
41°22'20.0
|
23°20'31.5
|
10370,92
|
1150,18
|
9
|
Места
|
граница
|
4593064,434
|
253166,3964
|
41°27'04.4
|
24°02'41.9
|
3263,63
|
967,41
|
10
|
Доспат
|
граница
|
4598921,311
|
254963,883
|
41°30'16.0
|
24°03'50.7
|
665,42
|
311,27
|
11
|
Доспат
|
с. Медени поляни
|
4639355,561
|
234787,0935
|
41°51'42.0
|
23°48'17.0
|
11,45
|
5,64
|
12
|
Струма
|
с. Дамяница
|
4602053,145
|
187316,271
|
41°30'32.7
|
23°15'12.3
|
8087,06
|
912,01
|
13
|
Струма
|
с. Струма
|
4604288,559
|
187216,5375
|
41°31'44.9
|
23°15'03.8
|
8042,63
|
907,33
|
14
|
Струма
|
с. Струма
|
4605777,967
|
186222,4209
|
41°32'31.7
|
23°14'18.2
|
7575,73
|
858,10
|
15
|
Сaндaнскa Бистрицa
|
гр. Сандански
|
4604758,344
|
187233,4253
|
41°32'00.1
|
23°15'03.7
|
139,49
|
47,94
|
16
|
Лебница
|
с. Лебница
|
4604173,863
|
186379,5926
|
41°31'40.0
|
23°14'28.0
|
322,81
|
71,50
|
17
|
Струма
|
с. Микрево
|
4615654,906
|
183823,2319
|
41°37'47.9
|
23°12'16.2
|
7407,36
|
840,30
|
18
|
Белишка
|
с. Струмяни
|
4616227,399
|
183333,338
|
41°38'05.7
|
23°11'54.0
|
28,74
|
22,59
|
19
|
Цaпaревскa
|
с. Микрево
|
4615454,266
|
183493,3984
|
41°37'40.9
|
23°12'02.3
|
78,58
|
36,47
|
20
|
Струма
|
с. Струмяни
|
4616530,701
|
182978,7828
|
41°38'15.0
|
23°11'38.1
|
7297,18
|
828,63
|
21
|
Струма
|
гр. Кресна
|
4626664,305
|
180405,9647
|
41°43'39.1
|
23°09'27.6
|
107,70
|
42,38
|
23
|
Брезнишка
|
гр. Кресна
|
4627579,488
|
179356,335
|
41°44'07.2
|
23°08'40.5
|
49,18
|
29,18
|
26
|
Сушичка
|
с. Полена
|
4639599,222
|
176272,1995
|
41°50'31.6
|
23°06'03.7
|
53,28
|
30,31
|
27
|
Струма
|
с. Крупник
|
4640243,101
|
179231,0476
|
41°50'56.7
|
23°08'10.5
|
6785,75
|
774,33
|
28
|
Брежaнскa
|
с. Полето
|
4640540,49
|
179351,0163
|
41°51'06.5
|
23°08'15.1
|
33,64
|
24,35
|
29
|
Струма
|
с. Полето
|
46410,5
|
179311,695
|
41°51'23.3
|
23°0812.4
|
666,57
|
761,63
|
30
|
Струма
|
гр. Симитли
|
4646200,028
|
177791,0834
|
41°54'07.3
|
23°06'56.5
|
6378,66
|
730,91
|
31
|
Грaдевскa
|
гр. Симитли
|
4644992,174
|
178043,7359
|
41°53'28.6
|
23°07'09.8
|
228,83
|
60,69
|
32
|
Струма
|
гр. Симитли
|
4644833,276
|
177498,0586
|
41°53'22.7
|
23°06'46.5
|
6619,47
|
756,61
|
33
|
Струма
|
гр. Благоевград
|
4654652,178
|
175246,8366
|
41°58'37.0
|
23°04'49.6
|
6198,19
|
711,60
|
34
|
Четиркa
|
гр. Благоевград
|
4655065,534
|
174375,6881
|
41°58'49.1
|
23°04'11.0
|
167,87
|
83,34
|
35
|
Струма
|
гр. Благоевград
|
4655694,136
|
174358,3139
|
41°59'09.4
|
23°04'09.0
|
6027,98
|
693,35
|
36
|
Блaгоевгрaдскa Бистрицa
|
гр. Благоевград
|
4656195,094
|
174706,3878
|
41°59'26.1
|
23°04'23.1
|
215,20
|
95,71
|
37
|
Струма
|
гр. Благоевград
|
4656739,663
|
173907,077
|
41°59'42.5
|
23°03'47.3
|
5797,39
|
668,58
|
38
|
Рилска
|
гр. Кочериново
|
4665245,902
|
175059,2029
|
42°04'19.4
|
23°04'20.3
|
356,28
|
126,73
|
39
|
Рилска
|
с. Стоб
|
4668139,432
|
178163,6719
|
42°05'57.6
|
23°06'29.4
|
337,68
|
123,00
|
40
|
Струма
|
гр. Бобошево
|
4674680,241
|
169354,0678
|
42°09'16.0
|
22°59'53.4
|
4334,04
|
509,69
|
41
|
Елешница
|
с. Друмохар
|
4680756,117
|
159072,7549
|
42°12'16.6
|
22°52'13.8
|
326,33
|
120,68
|
42
|
Джерман
|
гр. Дупница
|
4684799,911
|
177638,4122
|
42°14'55.8
|
23°05'33.4
|
471,60
|
148,15
|
45
|
Отовица
|
гр. Дупница
|
4687293,544
|
181234,8374
|
42°16'21.7
|
23°08'05.0
|
32,52
|
33,41
|
46
|
Джерман
|
гр. Дупница
|
4687296,311
|
180674,3901
|
42°16'21.0
|
23°07'40.6
|
399,84
|
135,14
|
47
|
Тополница
|
гр. Дупница
|
4687770,937
|
180802,3787
|
42°16'36.5
|
23°07'45.3
|
255,56
|
105,32
|
48
|
Джерман
|
гр. Дупница
|
4687765,597
|
181019,909
|
42°16'36.7
|
23°07'54.8
|
111,18
|
66,25
|
50
|
Арката
|
с. Бобарци
|
4708786,049
|
168657,3025
|
42°27'38.0
|
22°58'12.9
|
293,14
|
113,68
|
51
|
Струма
|
с. Прибой
|
4712716,329
|
165085,6477
|
42°29'39.6
|
22°55'28.6
|
783,63
|
196,57
|
52
|
Арката
|
с. Прибой
|
4712650,213
|
163715,9639
|
42°29'35.3
|
22°54'28.9
|
349,26
|
125,33
|
53
|
Струма
|
с. Прибой
|
4713556,602
|
163406,7527
|
42°30'04.1
|
22°54'13.5
|
1140,96
|
242,31
|
54
|
Струма
|
с. Жабляно
|
4713089,743
|
154533,8203
|
42°29'35.0
|
22°47'46.7
|
1518,70
|
284,14
|
55
|
Треклянска
|
гр. Земен
|
4712706,362
|
150938,9618
|
42°29'16.8
|
22°45'10.5
|
526,57
|
157,53
|
56
|
Струма
|
гр. Земен
|
4712287,776
|
151063,8447
|
42°29'03.4
|
22°45'16.9
|
1570,81
|
289,53
|
57
|
Струма
|
гр. Земен
|
4710653,918
|
148770,1564
|
42°28'06.9
|
22°43'40.3
|
2110,72
|
341,31
|
58
|
Струма
|
гр. Радомир
|
4719471,494
|
167383,7812
|
42°33'21.6
|
22°56'54.9
|
725,34
|
188,28
|
59
|
Струма
|
гр. Батановци
|
4723873,758
|
169260,4773
|
42°35'46.9
|
22°58'07.8
|
310,83
|
117,45
|
60
|
Струма
|
с. Студена
|
4716193,231
|
183071,0957
|
42°31'59.3
|
23°08'27.8
|
93,00
|
59,98
|
61
|
Струма
|
гр. Перник
|
472448,443
|
176431,5012
|
42°36'17.7
|
23°03'20.5
|
241,75
|
102,11
|
62
|
Мещичка
|
с. Мещица
|
4732419,899
|
172649,8743
|
42°40'28.5
|
23°00'18.5
|
14,23
|
21,09
|
63
|
Конска
|
с. Ярджиловци
|
4728848,805
|
167002,2542
|
42°38'24.3
|
22°56'18.5
|
248,29
|
103,64
|
64
|
Конска
|
гр. Батановци
|
4724494,752
|
167482,9689
|
42°36'04.2
|
22°56'48.7
|
370,75
|
129,57
|
65
|
Струма
|
с. Раждавица
|
4701785,704
|
146485,6861
|
42°23'16.4
|
22°42'20.2
|
2194,03
|
348,75
|
66
|
Глогошка
|
гр. Кюстендил
|
4689959,161
|
143297,7775
|
42°16'48.9
|
22°40'27.4
|
67,89
|
50,34
|
67
|
Рударщица
|
с. Рударци
|
4722136,652
|
184692,8122
|
42°35'13.9
|
23°09'26.9
|
19,51
|
25,14
|
68
|
Рударщица
|
с. Драгичево
|
4723920,06
|
182783,953
|
42°36'08.8
|
23°07'59.8
|
24,87
|
28,78
|
69
|
Арката
|
с. Горна Диканя
|
4704599,142
|
186496,2221
|
42°25'49.3
|
23°11'20.4
|
8,01
|
15,31
|
70
|
Тролско дере
|
с. Сатовча
|
4612233,28
|
247753,6549
|
41°37'19.1
|
23°58'20.0
|
14,80
|
19,88
|
71
|
Дяволски дол
|
с. Плетена
|
4614288,935
|
246293,3716
|
41°38'23.9
|
23°57'13.8
|
22,83
|
27,20
|
72
|
Доспат
|
с. Кочан
|
4608519,278
|
252320,4209
|
41°35'23.9
|
24°01'42.6
|
9,64
|
14,58
|
74
|
Голобрешка
|
с. Драгичево
|
4724138,405
|
182945,0523
|
42°36'16.1
|
23°08'06.4
|
14,73
|
21,50
|
75
|
Рударщица
|
гр. Перник
|
4723764,541
|
181529,018
|
42°36'01.9
|
23°07'05.1
|
45,50
|
40,28
|
77
|
Струма
|
Копиловци
|
4695361,237
|
149177,965
|
42°19'53.1
|
22°44'31.7
|
3098
|
422,63
|
78
|
Струма
|
Гирчевци_Жабокрът
|
4691995,037
|
151152,491
|
42°18'07.5
|
22°46'05.0
|
3609,7
|
460,19
|
79
|
Струма
|
Нови чифлик
|
4688380,136
|
154432,923
|
42°16'15.8
|
22°48'35.7
|
3703,92
|
466,84
|
80
|
Струма
|
Невестино
|
4686916,706
|
157450,932
|
42°15'33.3
|
22°50'50.2
|
3778,48
|
472,05
|
81
|
Елешница
|
с. Четирци
|
4684987,548
|
159672,86
|
42°14'34.4
|
22°52'31.0
|
356,13
|
126,70
|
83
|
Соволянска Бистрица
|
с. Соволяно
|
4696080,204
|
142528,026
|
42°20'05.5
|
22°39'40.3
|
264,05
|
107,25
|
84
|
Струма
|
Шишковци
|
4697570,408
|
147952,104
|
42°21'02.5
|
22°43'33.4
|
3082,71
|
421,47
|
85
|
Драговищица
|
Драговищица
|
4700087,315
|
142411,159
|
42°22'14.8
|
22°39'26.3
|
843,84
|
20484
|
86
|
Драговищица
|
Горановци
|
4702135,727
|
140068,551
|
42°23'17.2
|
22°37'39.5
|
828,79
|
202,80
|
87
|
Новоселска
|
Слокощица
|
4688404,207
|
146243,017
|
42°16'03.4
|
22°42'39.1
|
72,91
|
52,38
|
88
|
Рашка Гращица
|
Неделкова Гращица
|
4684583,097
|
153535,556
|
42°14'11.6
|
22°48'04.7
|
57,02
|
45,68
|
|