Исследование характеристик речного бассейна в районе чрезвычайной ситуации

Размещено на http://www.allbest.ru/

Новосибирская государственная академия

Водного транспорта

Кафедра водных изысканий и гидрогеологии

Курсовая работа

По дисциплине «Русловые изыскания»

Тема: «Исследование характеристик речного бассейна в районе чрезвычайной ситуации»

Студент гр. Чекалина Д.С.

Работа зачтена:

Преподаватель: Шамова В.В.

Новосибирск-2011

СОДЕРЖАНИЕ

ВВЕДЕНИЕ

1. ФИЗИКО-ГЕОГРАФИЧЕСКАЯ ХАРАКТЕРИСТИК ИССЛЕДУЕМОГО РАЙОНА

1.1 Описание рельефа местности

1.2 Характеристика ландшафта и геологические характеристики

2. ХАРАКТЕРИСТИКА СТРАТЕГИЧЕСКИХ ОБЪЕКТОВ

3. ХАРАКТЕРИСТИКА ВОДНЫХ ОБЪЕКТОВ

3.1 Физико-географические характеристики бассейна

3.2 Определение морфологических характеристик бассейна реки

3.3 Антропогенное влияние на речной сток

3.4 Построение продольного профиля реки и определение зоны затопления

3.5Вычисление средней скорости и оценка размываемости русла

3.6 Изменение скорости по вертикали

3.7 Определение скорости поперечной циркуляции

3.7.1 На повороте русла, от действия сил инерции

3.7.2 От действия сил Кориолиса

3.7.3 Определение суммарной скорости, построение эпюры

4. ВЫЧИСЛЕНИЕ РАСХОДА ВОДЫ

5. ВЫЧИСЛЕНИЕ РАСХОДА ВЗВЕШЕННЫХ НАНОСОВ

6. ХАРАКТЕРИСТИКА ДОРОЖНОЙ СЕТИ

6.1 Общее описание дорог

6.2 Выбор оптимального подъезда к стратегическому объекту при возникновении чрезвычайной ситуации

7. ИССЛЕДОВАНИЕ ВОЗДЕЙСТВИЯ ВЕТРА

7.1 Построение розы ветров

7.2 Расчет сгонно-нагонных течений

ЗАКЛЮЧЕНИЕ

СПИСОК ИСПОЛЬЗУЕМОЙ ЛИТЕРАТУРЫ

ВВЕДЕНИЕ

бассейн река чрезвычайная ситуация

В данной работе будет выполнен подробный анализ физико-географических условий, строительных и водных объектов, населенных пунктов и автодорог, проверка вероятности их затопления. Будет произведен расчет характеристик подъездных дорог к стратегическим объектам, выбор наиболее оптимального подъезда в зону ЧС, расчеты скоростного режима русла реки, вычисления расходов воды и наносов, вычисление сгонно-нагонных течений.

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

1. ФИЗИКО-ГЕОГРАФИЧЕСКАЯ ХАРАКТЕРИСТИКА ИССЛЕДУЕМОГО РАЙОНА

1.1 Описание рельефа местности

Рельеф участка представлен холмисто - горной местностью, левый берег р. Соть, приток р. Тихая заболочен. Р. Соть течет с северо-запада на юго-восток. На юго-западной части рельеф более холмистый, на северо-востоке более крутой с горами и крутыми обрывами. В северо-западной и восточной части располагается цепь холмов и гор, направленных от водного объекта. В средней части карты преобладает горный рельеф.

Таблица 1.1 - Характерные точки

Название точки	Высота z,м	Географические координаты	Прямоугольные координаты
		Долгота л (восточная)	Широта ц (северная)	X	У
Наивысшая	259,4	18?8ґ37Ѕ	54?49ґ6Ѕ	6079км 900м	4316км 400м
Наинизшая	108	18?5ґ45Ѕ	54?46ґ52Ѕ	6071км 900м	4313км 459м

z_max-z_min;

151,4м.

1.2 Характеристика ландшафта и геологические характеристики

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

В юго-западной части рельеф горно-холмистый с большим количеством оврагов, обрывов и седловин. Грунт черноземный. Протекают р. Соть с песчаными берегами и впадающая в неё реки Тихая и Калитинка. Берега р. Тихой окружены заболоченными пастбищами с отдельно растущими кустарниками. Имеется оз. Щучье, заболоченное с зарослями камыша и тростника. Вдоль реки Калитинка растет сосново-дубовый лес. Располагаются большие населенные пункты Ивановка, Борисово. Находится аэропорт. Проходят железная дорога и автомобильные дороги без покрытия и с шоссейным и усовершенствованным покрытиями.

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

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

2. ХАРАКТЕРИСТИКА СТРАТЕГИЧЕСКИХ ОБЪЕКТОВ

Таблица 2.1- Характеристика стратегического объекта

Название объекта	Высота z, м	Долгота л (восточная)	Широта ц (северная)	Координаты	Площадь км2	Расстояние L, км
				X	Y
Ивановка	120	18?10ґ00Ѕ	54?73ґ00Ѕ	6072-6073	4309-4312	4	Lmin=0,1 Lmax=0,85

3. ХАРАКТЕРИСТИКА ВОДНЫХ ОБЪЕКТОВ

Правильная организация работ при ЧС на гидротехнических объектах требует знания морфометрических, физико-географических, скоростных и других характеристик речного бассейна.

Бассейн - часть земной поверхности, включая толщу почво-грунтов, откуда происходит сток воды в отдельную реку или речную систему. Площадь каждого речного бассейна ограничена водоразделом. Бассейн, определяемый по топографической карте для главной реки до её устья, состоит из бассейнов отдельных притоков, каждый из которых в свою очередь имеет бассейны притоков следующего порядка.

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

3.1 Физико-географические характеристики речного бассейна

1) Географическое положение бассейна.

Таблица 3.1- Характеристики бассейна реки.

Название	Порядок	Длинна, км	Координаты истока	Координаты устья
			X,Y	л , ц	X,Y	л , ц
Соть	главная	7,25 L=6	6072км 850м 4306км 850м	18?00?00? 54?44?55?	6071км 900м 4313км 450м	18?5?45? 54?46?52?
Тихая	1	12,5 L=10	4309км 450м 6082км 10м		6072км 550м 4312км 750м

К_изв=l_р/L,

где l_р, км - длина по руслу реки; L, км - длина по прямой.

K_изв(Соть)=,

К_изв(Тихая)=,

К_разв. = (Уli + 1с) / 1с,

где Уli - сумма длин несудоходных рукавов,

1с - длина судоходного рукава.

К_разв=1.

2) Климатические условия бассейна.

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

Главнейшими климатическими факторами являются:

Среднее годовое количество осадков: x,

Средний годовой слой испарения: z,

Средний годовой сток «y», находится как разность y=z-x,

Средний годовой сток: y.

3) Геологическое строение бассейна.

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

4) Рельеф бассейна.

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

Z_cp. = w / F,

где w - площадь в масштабах графика, ограниченная гипсографической кривой и осями координат, выраженная в км³;

F - общая площадь бассейна в квадратных километрах.

Для построения гипсографической кривой произведём расчет в таблице:

Таблица 3.2 - Данные для построения гипсографической кривой бассейна р. Соть.

Высотные зоны, м	Площадь высотной зоны бассейна, км2	Высота бассейна, м zi	Площадь бассейна, соответствующая высоте, км2	Нарастание площади бассейна, км2
233-200	12,55	233	0.00	0,00
200-150	26,5	200	12,55	12,55
150-108	50	150 108	26,5 50	39,05 89,05

Нарастание площади бассейна с уменьшением высоты получаем путем последовательного суммирования площадей между горизонталями. Также строим график распределения площадей бассейна по высотным зонам. На этом графике значение площадей между горизонталями относим к середине высотной зоны (интервала). Точки гипсографической кривой относим к нижней границе высотной зоны. Под горизонтальной шкалой площадей дополнительно наносим шкалу процентов из расчета F=89,05 км² и составляет 100%.

По гипсографической кривой, устанавливаем среднюю высоту бассейна. Для этого определяем площадь, ограниченную гипсографической кривой и осями координат, которая в масштабе графика составляет

км³.

Подставляя вычисленные значения в формулу, имеем:

Z_cp. ==91,97м.

Второй способ определения средней высоты бассейна выполняем по карте с горизонталями по формуле:

Z_cp. = (Z₁f₁ + Z₂f₂ + Z₃f₃+ ...+Z_nf_n) / F,

где f₁, f₂, f₃ ,... f_n - площади между соседними горизонтами (или в пределах высотных интервалов), км²;

Z₁, Z₂, Z₃ , ... Z_n - средние высоты этих интервалов, м;

F - общая площадь бассейна, км².

Для выполнения этой работы удобно пользоваться палеткой. Расчет произведём в таблице 3.3.

Таблица 3.3 - Данные к вычислению средней высоты бассейна р. Соть

№	Отметки горизонталей, м	Площадь бассейна, fi, тыс. км2	Средние высоты между горизонталями Zi, м	Произведение Zi fi м тыс. км2
1	233-200
2	200-150
3	150-108
Сумма

Третий способ определения средней высоты заключается в следующем: делением площади бассейна на квадраты равного размера, и для каждого из которых определяем среднюю высоту z _ср.к.

Средняя высота бассейна равна сумме средних высот квадратов, деленной на количество этих квадратов n :

Z_cp.= УZ_ср.к/n

Z_ср.к= (170+150+179+150+180+150+188+175)/6=167,7м.

Для определения среднего уклона поверхности бассейна применяем следующую формулу:

I_ср. = h(1 ₀/2 + 1₁+1₂+1₃+…+1_n2)/F,

где h- высота сечения рельефа местности в метрах;

I₀, I₁, I₂,....длина горизонталей в метрах;

F - площадь бассейна реки в км². Вычисление уклона производится в таблице 3.4.

Таблица 3.4 - Вычисление уклона бассейна реки Соть

№	Отметки горизонталей, м	Длина горизонталей, км	Полусумма длин горизонталей, км	Превышение между горизонталями, h, км	Произведение
1
2
3
4
Icp=0,23725/5,425=0,0437

5) Растительный покров. Для его характеристики необходимо определить площади, занятые лесом и другой растительностью. Степень залесенности речного бассейна определяется коэффициентом лесистости, представляющим собой частное от деления площади лесов F_лес, находящихся на территории бассейна, на общую площадь бассейна F, т.е.

6) Озерность бассейна. Озерность бассейна определяется коэффициентом озерности, т. е. отношением площади зеркала всех водоемов ?F_оз (озер, водохранилищ) к общей площади бассейна

7) Заболоченность речного бассейна. Определяется отношением площади, занятой в бассейне болотами F_бол, ко всей его площади:

8) Наличие ледников и вечномерзлых пород. Ледники и вечномерзлые породы на местности отсутствуют.

3.2 Определение морфометрических характеристик бассейна реки

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

1) Площадь бассейна F, км². Площадь бассейна определяется несколькими способами: при помощи палетки, планиметра, геодезических таблиц, графическим методом путем разбиения всей площади на правильные геометрические фигуры с вычислением их площадей и последующим суммированием.

F=89,05 км²

2) Длина бассейна L, км - расстояние по прямой от устья реки до наиболее удаленной точки бассейна.

L=12,5км

3) Наибольшая ширина бассейна В, км - длина линии, перпендикулярной длине бассейна в наиболее широком месте.

В=15 км

4) Средняя ширина бассейна В_ср., км - определяется путем деления площади бассейна на его длину L, т.е.

В_ср. = F / L,

В_ср. = 89,05 / 12,5=7,12 км.

Иногда определяют ширину отдельно для левой и правой частей бассейна.

где F_л. - площадь левобережной части бассейна в км²;
F_пp. - площадь правобережной части бассейна в км².

5) Коэффициент ассиметрии бассейна. Главная река может занимать симметричное положение (посреди бассейна) или боковое, т.е. подходить к одному из водоразделов. Мерой ассиметричности является коэффициент, определяемый по формуле

а = (F_л -F_пр )/0,5(F_л + F_пр ),

а=(82,8-6,25)/0,5(82,8+6,25)=1,7,

где F_л - площадь левобережной части бассейна в км²;

F_пр - площадь правобережной части бассейна в км².

Иногда коэффициент ассиметрии речного бассейна определяют отношением левой части площади бассейна к правой:

6) График нарастания площади бассейна реки характеризует постепенное увеличение площади бассейна реки по длине от истока к устью. Выделяются главная река и притоки, проводятся водораздельные линии их бассейнов, определяются площади бассейнов притоков, главной реки и межприточных участков. Результаты вычислений заносятся в таблицу 3.5. На основании данных таблицы строится график, на котором по горизонтальной оси откладывается длина главной реки в принятом масштабе, а по вертикальной - площадь межприточных участков и площади бассейнов притоков в створах впадения их в главную реку. Наклонные линии графика показывают постепенное нарастание площадей межприточных участков главной реки. В местах впадения притоков в главную реку откладываются в масштабе отрезки вертикальных линий, показывающие увеличение площади бассейна за счет площади бассейна притока. Суммарный график строится последовательным суммированием ординат площадей графиков нарастания левой и правой частей бассейна.

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

Таблица 3.5 - Данные к построению графика нарастания площади бассейна р.Соть

Название участка	Расст. от истока, км	Площадь, F	Название участка	Расст. От истока, км	Площадь, F
Правый берег	Левый берег
Межприточный участок 1	-	0,75	0,75	Межприточный участок 3	-	8,5	8,5
Приточный 1 р.Андора	1,5	1,5	2,25	Приточный 2 р. Тихая	6,25	51,5	60
Межприточный участок 2	-	3,5	5,75	Межприточный участок 4	-	23	83

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

Таблица 3.6 - Данные к построению кругового графика р. Соть

Название притока	ПЛОЩАДЬ	Название притока	ПЛОЩАДЬ
	F, кв.км	F,%	градусы		F, кв.км	F,%	градусы
ПРАВЫЙ БЕРЕГ	ЛЕВЫЙ БЕРЕГ
М1	0,75	0,8	3	М3	8,5	9	32
П1	1,5	1,7	6	П2	51,5	54,3	195
М2	3,5	4	14	М4	23	30,2	110

График нарастания площади бассейна по длине и круговой график позволяют правильно распределить ресурсы (питание, оборудование, технику, обмундирование и т.д.), рационально организовать спасательные работы и экстренную ликвидацию чрезвычайной ситуации на объекте.

3.3 Антропогенное влияние на речной сток

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

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

3.4 Построение продольного и поперечного профилей реки и определение зоны затопления

Построение продольного профиля с вычислением отметок уровней воды весеннего половодья или паводка необходимо при определении зон затопления и организации работ по ликвидации чрезвычайных ситуаций.

Построение ведется в следующей последовательности:

1) Отложим длину реки в километрах в масштабе по горизонтали. Измеряем и откладываем в масштабе расстояния между горизонталями;

2) Вычисляем отметки максимального уровня воды в начале и конце исследуемого участка реки по формуле

Z_Hmax = Z₀ + H_max,

где Z₀ - отметка нуля гидропоста, м;

H_max - максимальный уровень, берём из таблиц гидрологического ежегодника в один день для верхнего и нижнего гидропостов (в метрах);

Отметки нулей верхнего и нижнего гидропостов определяются расположением плоскости отчета уровней воды ниже самых низких уровней за многолетний период. Отметку нуля Z₀ гидропоста ориентировочно (в первом приближении) принимаем на три метра ниже отметки берега.

Z_{0 истока} = 109,8 - 3= 106,8 м,

Z_{0 устья} = 108 - 3 = 105 м,

H_max= 6,80 м - для нижнего гидропоста,

H_max = 6,35 м - для верхнего гидропоста,

Z_Hmax = 106,8 + 6,35 = 113,15 м - для истока,

Z_Hmax = 105 + 6,80 = 111,8 м - для устья.

3) По вычисленным отметкам максимального уровня воды проводим линию водной поверхности и выполняем анализ зон затопления; на плане эти зоны закрашиваем голубым цветом;

4) Вычисляем отметки проектных уровней Н_пр. для верхнего и нижнего гидропостов:

Z_{Нпр в} = Z_Ов+ Н_{пр в ,}

Z_{Нпр н} = Z_Он+ Н_{пр н,}

Z_{Нпр в} = 106,8 + 0,35 = 107,15 м - для истока,

Z_{Нпр н} = 105 + 0,70 = 105,7 м - для устья.

Линия, проведенная по отметкам проектного уровня, оформляется красным цветом.

5) Вычисляем уклон водной поверхности при проектном уровне воды, по формуле:

J =(Z_{Нпр истока} -Z_{Нпр устья})/L,

J =(107,15 -105,7)/7,25=0,2,

Таблица 3.7 - Вычисление уклона водной поверхности

Г/пост	Отметка нуля г/поста, м	Проектный уровень, м	Отметка проектного уровня, м	Общее падение, м	Расстояние м/у г/постами, км	Осредненное падение м/км	Уклон
Верхний	106,8	0,35	107,15	1,45	7,25	0,2	0,2
Нижний	105	0,7	105,7

Вывод: Так как отметки максимальных уровней воды превышают отметки берега и проектные уровни воды, то произойдет затопление прибрежных территорий. Стратегическому объекту опасность затопления не угрожает, так как в зону наводнения он не попадает.

3.5. Вычисление средней скорости и оценка размываемости русла

Для вычисления средней скорости расхода воды и средних наносов, строится поперечный профиль реки по заданным значениям. Для построения берём масштабированную бумагу и по горизонтали откладываем B - ширину русла, а по вертикали h - высоту. Ширина разделяется на четыре участка B₁=40м, B₂=30м, B₃=30, B₄=30м. По этим сечениям откладывается глубина h₁=1.2м, h₂=3.3м, h₃=2.8м. Получаем четыре участка, площадь которых будет равна w-площади поперечного сечения русла. Вычисляем эти площади как площади геометрических фигур.

₁=; ₁==24 м²,

₂=; ₂==69м²,

₃=; ₃==153 м²,

₄=; ₄==45 м²,

=,

=24+69+153+45=291 м².

Средняя скорость течения в поперечном сечении русла определяется по формуле Шези:

где С - коэффициент Шези;

R - гидравлический радиус, м;

I - продольный уклон водной поверхности.

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

С = R 1/6 / n,

где n - коэффициент шероховатости, равный n=0,027,

R = w / - гидравлический радиус, м

w - площадь поперечного сечения русла, м²;

- длина смоченного периметра русла, м;

Для равнинных рек гидравлический радиус равен средней глубине:

R = h_ср

где h_ср = w / В =291/130,2=2,2м;

В - ширина русла, м;

Коэффициент шероховатости n можно определить по формуле:

,

С = 2,2 1/6 / 0,03=37,9 ,

V _ср. = 37,9=0,83 м/с.

Средняя скорость течения V_ср. на вертикали, отвечающая состоянию предельного равновесия донных частиц, когда отдельные частицы срываются с места и перемещаются, но общего движения наносов еще нет, называется неразмывающей V_нр. Она является предельной скоростью, отвечающей началу сдвига отдельных частиц, т.е. скоростью начала влечения и определяется с помощью следующих зависимостей.

Формула В.Н.Гончарова, которая по данным исследований дает на песчаных перекатах наилучшие результаты.

где h - средняя глубина потока, м;

d - средний диаметр донных частиц, м;

- диаметр частиц с обеспеченностью 95% по кривой гранулометрического состава, т.е. такой диаметр, который оказывается превзойденным лишь у 5% частиц, м.

=3,74

Скорость течения , при которой движение донных наносов становится массовым, называется размывающей. Соотношение между размывающей и неразмывающей скоростями равно 1,30.

;

м/с.

Вывод:

- русло размываться не будет,

- русло будет размываться.

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

0,83<4,42 м/с.

3.6 Изменение скорости по вертикали

Для аналитического выражения закономерности изменения осредненных скоростей течения по глубине (на вертикали) воспользуемся формулой А.В. Краушева для вычисления поверхностной скорости течения:

V=V_п , м/с,

h - глубина на вертикали, принимается равной h₂, м.

у/ h - относительная глубина точки наблюдения,

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

Vп - поверхностная скорость течения, м/с;

V_п = 1,11C V _ср. / (С - 1),

V_п = = 0,95 м/с.

Р - безразмерный параметр, определяемый по выражению.

Р= 0,57 +3,3/С при 10<C<60,

Р= 0,57 +3,3/37,9 = 0,1.

Порядок расчета следующий:

Используя исходные данные, вычисляем скорости течения в шести точках на вертикали у/ h =0.02; у/ h =0,2; у/ h =0,4; у/ h =0,6; у/ h =0,8; у/ h =1.

V_0.0=0,95=0,95 м/с,

V_0.2=0,95=0,94 м/с,

V_0.4=0,95=0,93 м/с,

V_0.6=0,95=0,92 м/с,

V_0.8=0,95=0,91 м/с,

V₁=0,95=0,9 м/с.

По вычисленным значениям строится эпюра скорости, принимая масштаб: для глубины в 1см-0,5 м, для скоростей в 1см-1 м/c.



3.7 Определение скорости поперечной циркуляции

3.7.1 На повороте русла

При изменении направления потока возникают центробежные силы инерции, направленные по нормали к криволинейным линиям тока в плане. Они приводят к образованию поперечного уклона I_поп. и поперечных составляющих скорости, перпендикулярных к основному направлению потока.

Поперечный уклон определяется по формуле

I_поп.= V_ср² / rg,

где: скорость Vср определяется по формуле Шези;

g = 9,8 м/с²,

I_поп= 0,83²/ 5009,8=0,00014.

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

H= I_поп.В ;

H= 0,00014130=0,018 м.

Для определения величины поперечной составляющей скорости V_поп. на вертикали, находящейся на повороте русла, воспользуемся формулой К.И. Россинского и И.А. Кузьмина

V_поп= (1,53C² V_срh) / g r (у1 / h) ^0,15 (у1 /h) ^0,3 - 0,80, м/ c,

где у1- высота точки над дном, м;

h - глубина потока по судовому ходу на повороте русла, м.

Задавая относительную глубину у1/h, получаем значение поперечной скорости V_поп в этой точке на вертикали.. Величины у1/h принимаем следующую:

у1/h = 1 (точка находится на поверхности у1 = h);

у1/h =0,8 (вторая точка);

у1/h = 0,6 (третья точка);

у1/h = 0,4 (четвертая точка);

у1/h = 0,2 (пятая точка).

у1/h =0,02 - на дне (шестая точка).

V_поп.= (1,53 0,833,4) / 9,8 500 (у1 / h) ^0,15 (у1 /h) ^0,3 - 0,80=

= м/ c

V₁= =0,254 м/с,

V_0.8= =0,172 м/с,

V_0.6= =0,07 м/с,

V_0.4= = - 0,044 м/с,

V_0.2= = - 0,18 м/с,

V_0.01= = - 0,35 м/с.

По вычисленным значениям строится эпюра скорости, принимая масштаб: для глубины в 1см-0,5 м, для скоростей в 1см-1 м/c.

3.7.2 От действия сил Кориолиса

1) Определяется поперечный уклон:

I_k=(0,0001456Vsin )/g;

I_k =0,00014560,83sin54/9.8=0,0000099

Соответственно перепад уровней воды:

м.

Для определения составляющей скорости Кориолиса V_к, воспользуемся формулой К.И. Россинского и И.А. Кузьмина:

V_к.= (1,53C²wh sin) (у1 / h) ^0,15 (у1 /h) ^0,65 - 0,89, м/ c,

где w =0.0000728- угловая скорость вращения Земли, рад/с;

h=3,3-глубина потока, м;

sin=0,8,

С= 40,6 - коэффициент Шези;

V_к.= (2,65² sin54) (у1 / h) ^0,15 (у1 /h) ^0,65 - 0,89=

=0,22(у1 / h) ^0,15 (у1 /h) ^0,65 - 0,89,

V_к1= 0,221 ^0,15 1 ^0,65 - 0,89= 0,024 м/ c ,

V_к0.8= 0,22(0,8) ^0,15 (0,8) ^0,65 - 0,89= - 0,006 м/ c,

V_к0.6= 0,220,6 ^0,15 0,6 ^0,65 - 0,89= - 0,035 м/ c,

V_к0.4= 0,220,4 ^0,15 0,4 ^0,65 - 0,89= - 0,065м/ c,

V_к0.2= 0,220,2 ^0,15 0,2 ^0,65 - 0,89= - 0,092 м/ c,

V_к0.01= 0,220,01 ^0,15 [0,01 ^0,65 - 0,89= - 0,093 м/ c.

По вычисленным значениям строится эпюра скорости, принимая масштаб: для глубины в 1см-0,5 м, для скоростей в 1см-1 м/c.

3.7.3 Определение суммарной скорости, построение эпюры

Суммарная скорость вычисляется путем суммирования V_поп - скорости поперечной циркуляции и V_к скорости Кориолиса на одной глубине.

V_с= V_поп+ V_к, м/с.

V_с1=0,254+0,024= 0,478 м/с.

V_с0.8=0,172-0,006= 0,166 м/с.

V_с0.6=0,07 -0,035= 0,035 м/с.

V_с0.4=-0,044 -0,065= - 0,109 м/с.

V_с0.2=-0,18-0,092= - 0,272 м/с.

V_с0.02=-0,35-0,093= - 0,443 м/с.

По вычисленным значениям строится эпюра скорости, принимая масштаб: для глубины в 1см-0,5 м, для скоростей в 1см-1 м/c.

Вывод: Превышение уровня воды у вогнутого берега над выпуклым больше от действия центробежной силы, чем от действия силы Кориолиса.

4. ВЫЧИСЛЕНИЕ РАСХОДА ВОДЫ

Порядок выполнения работы:

1) Вычисляем скорости течения воды в каждой из пяти точек на трех скоростных вертикалях и определяем тных вертикалях и определяютем средние скорости на каждой вертикали V ср.

Вычисления будем проводить в расчетной табл. 4.1 .

Уравнение скорости течения воды имеет вид:

V = a + b n, м / c,

где n - число оборотов в 1 секунду,

a=0,036 и b=0,755 - коэффициенты уравнения вертушки по варианту.

Для первой вертикали:

V_п = 0,036 + 0,755 0,74=0,59 м / c,

V_0.2 = 0,036 + 0,7550,72=0,58 м / c,

V_0.6 = 0,036 + 0,755 0,69=0,56 м / c,

V_0.8 = 0,036 + 0,755 0,54=0,44 м / c,

V_дна = 0,036 + 0,755 0,42 =0,35 м / c.

Для второй вертикали:

V_п = 0,036 + 0,7551,03=0,81 м / c,

V_0.2 = 0,036 + 0,7551,02 =0,80 м / c,

V_0.6 = 0,036 + 0,7550,92=0,73 м / c,

V_0.8 = 0,036 + 0,7550,79=0,63 м / c,

V_дна = 0,036 + 0,7350,60 =0,49 м / c.

Для третей вертикали:

V_п = 0,036 + 0,7550,90=0,72 м / c,

V_0.2 = 0,036 + 0,7550,95=0,75 м / c,

V_0.6 = 0,036 + 0,7550,92 =,73м / c,

V_0.8 = 0,036 + 0,55,74=0,59 м / c,

V_дна = 0,036 + 0,7550,50 =0,41 м / c.

Средняя скорость на каждой вертикали вычисляется по формуле:

V ср. = 0,1(Vпов + 3V + 3V + 2V + Vдно), м/с,

V _ср. = 0,1(0,59+ 30,58+ 30,56+ 20,44+0,35)= 0,52 м/с;

V _ср. = 0,1(0,81+ 30,80+ 3,73+ 20,63+0,49)=0,72 м/с;

V _ср. = 0,1(0,72+ 30,75+ 30,73 + 20,59 +0,41)=0,68 м/с.

Таблица 4.1 - Вычисление скоростей в живом сечении русла р. Соть

№ вертикали	Расстояние от постоянного начала l, м	Глубина на вертикали h, м	Точки Измерения на вертикали	Суммарное Число Оборотов в точке N	Время измерения в точке, с t	Число оборотов в 1 сек. n=N / t	Скорость течения в точке, м/с v =a+ b* n	Средняя скорость течения на вертикали м/с v
1	40	1,2	Пов-ть 0,2 0,6 0,8 Дно	120 120 80 80 60	163 167 116 148 144	0,74 0,72 0,69 0,54 0,42	0,59 0,58 0,56 0,44 0,35	0,52
2	70	3,4	Пов-ть 0,2 0,6 0,8 Дно	140 140 120 100 80	136 137 130 126 134	1,03 1,02 0,92 0,79 0,60	0,81 0,80 0,73 0,63 0,49	0,72
3	100	3,0	Пов-ть 0,2 0,6 0,8 Дно	160 140 120 100 80	178 147 130 135 160	0,90 0,95 0,92 0,74 0,50	0,72 0,75 0,73 0,59 0,41	0,68

2) Производим аналитическое вычисление расхода воды и площади живого сечения в таблице 4.2.

Средние скорости течения между вертикалями vⁿ_ср определяются полусуммами средних скоростей vср.i на смежных вертикалях, вычисленных в таблице 4.1

v_ср.1 = v_ср.10,7,м/с;

v_ср.2 = (v_ср.1 + v_ср.2) / 2 м/с;

v_ср.3 = (v_ср.2 + v_ср.3) / 2 м/с;

v_ср.4 = v_ср.30,7 м/с.

При этом для первой и последней части живого сечения среднюю скорость определяем умножением средней скорости на ближайшей вертикали на коэффициент, учитывающий распределение скоростей течения у берегов, коэффициент принимается равным 0,7 - для холмистого.

v_ср.1 = 0,520,7 =0,36 м/с,

v_ср.2 = (0,52+0,72) / 2=0,62 м/с,

v_ср.3 = (0,72+0,68) / 2 =0,7 м/с,

v_ср.4 = 0,680,7=0,48 м/с.

Расходы воды через части живого сечения qi находим как произведение их площадей wi на средние скорости vⁿ_ср между вертикалями:

q_i=w_i v_срi;

q₁=240,36=8,64 м/ с,

q₂ =690,62=42,78 м/ с,

q₃ =1530,7=107,1 м/ с,

q₄=450,48=21,6 м/ с.

Общий расход воды Q определяем суммированием расходов через элементы сечения:

Q =q_i м/с;

Q=8,64+42,78+107,1+21,6=180,12 м/с.

3) По полученным данным дополнительно вычисляем:

а) среднюю скорость в живом сечении:

V_ср. =Q / W , м/с;

V ср. = 180,12/ 291 =0,62 м/с.

б) среднюю глубину живого сечения:

h_ср = W / B м;

где B-ширина русла

h_ср = 291 / 130=2,24 м,

Таблица 4.2 - Вычисление расхода воды аналитическим способом

№ вертикалей	Расстояние от постоянного начала l,м	Площадь сечения между вертикалями w, м	Средняя скорость	Расход между вертикалями расчет вручную q , м/ с
			На вертикали V, м/с	Между вертикалями V , м/с
Урез левого берега	0
1	40	24	0,52	0,36	8,64
2	70	69	0,72	0,62	42,78
3	100	153	0,68	0,7	107,1
Урез правого берега	130	45		0,48	21,6
	В=130	W=291			Q =180,12

5. ВЫЧИСЛЕНИЕ РАСХОДА ВЗВЕШАНЫХ НАНОСОВ

1) Определяем мутность , единичные расходы , взвешенные наносы в каждой точке отбора проб, затем вычисляем аналитическим способом средние единичные расходы наносов на каждой вертикали. При этом принимаем, что пробы на мутность отбирались в тех же точках, в которых проводилось измерение скоростей течения. Вычисления произведём в табл. 5.1.

Величины мутностей в точках отбора проб находим делением весового количества наносов на объем пробы:

= (P 10⁶) / A , г/м3

где P - вес наносов, г;

A - объем пробы, см.

На первой вертикали:

_п = (0,054 106) / 3000=18,00г/м³

_0.2 = (0,052 106) / 3000=17,33 г/м³

_0.6 = (0,057 106) / 3000=19,00г/м³

_0.8 = (0,066106) / 3000=22,00г/м³

_дно = (0,096106) / 3000=32,00 г/м³

На второй вертикали:

_п = (0,071106) / 3000=23,70г/м³

_0.2 = (0,066106) / 3000=22,00г/м³

_0.6 = (0,072106) / 3000=24,00г/м³

_0.8 = (0,085106) / 3000=28,33г/м³

_дно = (0,095106) / 3000=31,66 г/м³

На третьей вертикали:

_п = (0,057106) / 3000=19,00г/м³

_0.2 = (0,057106) / 3000=19,00г/м³

_0.6 = (0,064106) / 3000=21,33г/м³

_0.8 = (0,068106) / 3000=22,66г/м³

_дно = (,073106) / 3000=24,33г/м³

Единичные расходы наносов в точках вычисляем по формуле:

= V , (г/ м²с),

где V- скорость течения в точке, м /c;

- мутность в той же точке, г / м³.

На первой вертикали:

= 0,5918,00=10,62 (г/ м²с),

= 0,5817,33=10,03 (г/ м²с),

= 0,5619,00= 10,64(г/ м²с),

= 0,44322,00= 9,68 (г/ м²с),

= 0,3532,00=11,2 (г/ м²с).

На второй вертикали:

= 0,8123,7=19,2 (г/ м²с),

= 0,8022,00=17,6 (г/ м²с),

= 0,73 24,00 =17,52 (г/ м²с),

= 0,6328,33 =17,83 (г/ м²с),

= 0,4931,66 =15,53 (г/ м²с).

На третьей вертикали:

= 0,7219,00=13,68 (г/ м²с),

= 0,7519,00=14,25 (г/ м²с),

= 0,7321,33=15,55 (г/ м²с),

= 0,5922,66 = 13,39 (г/ м²с),

= 0,41 24,33= 9,96 (г/ м²с).

Средний единичный расход наносов на вертикали определяем по формуле:

_ср.i =0,1(_пов.+ 3+3+ 2+_дно), г/ м²с;

_ср.1 =0,1(10,62+ 3 10,03+ 310,64+ 2 9,68 + 11,2)= 10,32 г/ м²с,

_ср.2 = 0,1(19,2+ 317,6+ 317,52+ 217,83 + 14,53)= 17,58 г/ м²с,

_ср.3 = 0,1(13,68+ 314,25 + 315,55+ 213,39 + 9,96)=13,98 г/ м²с.

2) Средние единичные расходы наносов между вертикалями определяем как полусуммы средних единичных расходов наносов на смежных вертикалях, которые вычислены в табл. 5.1. При этом для первой и последней частей живого сечения средний единичный расход определяем умножением среднего единичного расхода на ближайшей вертикали на коэффициент, равный 0.7 (аналогично определению скоростей течения между вертикалями в табл. 5.2).

Таблица 5.1 - Вычисление мутностей и единичных расходов взвешенных наносов в точках живого сечения

№	Расстояние от постоянного начала l, м	Глубина на вертикали h, м	Точки отбора проб на вертикали	Скорость течения в точке, м/с v	Объем пробы cм А	Вес наносов в пробе, г Р	Мутность точке, г / м	Единичный расход наносов в точке, , г / мс	Сред. един. Расход наносов на вертикали ,г/мc
1	40	1,2	Пов-ть 0.2 0.6 0.8 Дно	0,59 0,58 0,56 0,44 0,35	3000	0,054 0,052 0,057 0,066 0,096	18,00 17,33 19,00 22,00 32,00	10,62 10,03 10,64 9,68 11,2	10,32
2	70	3,4	Пов-ть 0.2 0.6 0.8 Дно	0,81 0,80 0,73 0,63 0,49	3000	0,071 0,066 0,072 0,085 0,095	23,70 22,00 24,00 28,33 31,66	19,2 17,6 17,52 17,83 15,53	17,58
3	100	3,0	Пов-ть 0.2 0.6 0.8 Дно	0,72 0,75 0,73 0,59 0,41	3000	0,057 0,057 0,064 0,068 0,073	19,00 19,00 21,33 22,66 24,33	13,68 14,25 15,55 13,39 9.,96	13,98

_ср1 =_ср1 0,7 =10,320,7=7,22 г/м²с,

_ср2 = (_ср1 +_ср2) / 2=(10,32+17,58) / 2=13,95 г/м²с,

_ср3 = (_ср2 +_ср3)/2 =(17,58+13,98) / 2 =15,78 г/м²с,

_ср4 =_ср3 0,7= 13,98 0,7=9,79 г/м²с.

Расходы наносов через части живого сечения r_i вычисляем умножением значений их площадей на средние единичные расходы между вертикалями с последующим делением на 10² для приведения размерности расхода наносов в кг/с.

r_i=_ср1 w/1000, кг/с;

r₁=7,2224 /1000= 0,17 кг/с,

r ₂=13,9569 /1000= 0,96кг/с,

r ₃=15,78 153/1000=2,41 кг/с,

r ₄=9,7945/1000=0,44 кг/с.

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

R= r_i, кг/ с;

R=0,17+0,96+2,41+ 0,44=3,98 кг/ с.

3) По полученным данным вычисляем среднюю мутность в сечении:

_ср. = 1000R / Q , г/ м,³

где R - расход взвешенных наносов, кг / c,

Q - расход воды, вычисленный аналитическим методом, м³/ c.

_ср. = 10003,98/180,12=22г/ м³

Таблица 5.2 - Вычисление расхода взвешенных наносов в живом сечении русла

№ вертикалей	Расстояние от постоянного начала l, м	Глубина На Вертикали h,м	Средний единичный расход наносов на вертикали, ,г/мс	Средний единичный расход наносов между вертикалями , г/мс	Площадь сечения между вертикалям w, м	Расход наносов между вертикалями (вручную) r,кг/с
Урез левого берега	0	0
1	40	1,2	10,32	7,22	24	0,17
2	70	3,4	17,58	13,95	69	0,96
3	100	3,0	13,98	15,78	153	2,41
Урез правого берега	130	0		9,79	45 W=291	0,44 R=3,98

6. ХАРАКТЕРИСТИКА ДОРОЖНОЙ СЕТИ

6.1 Общее описание дорог

К стратегическому объекту Ивановка подходят три шоссейные дороги в северном, западном, юго-восточном направлениях. Дорога в северном направлении пересекает р. Соть.

Дорога в западном направлении пересекает р. Андога по мосту, далее переходит в грунтовую дорогу в населенном пункте Коровино.

С северо-востока на юго-запад к проходит железная дорога.

6.2 Выбор оптимального подъезда к стратегическому пункту при возникновении чрезвычайной ситуации

Для подъезда к стратегическому объекту во время чрезвычайной ситуации лучше всего использовать две автомобильные дороги №1 и №2 (смотри карту) н.п. Дворики и н.п. Калитино.

Описание подъездных дорог.

Дорога №1.Проходит в северном направлении.

Тип покрытия - шоссейное улучшенное.

Уклон: ,

где - перепад высот в начале и конце дороги, L- длина дороги.

Дорога №2. Проходит в юго-восточном направлении.

Тип покрытия - шоссейное улучшенное покрытия

Уклон:

Вывод: На основе анализа наиболее оптимальным подъездным путем к н.п. Шустики является дорога №1, так как она короче, чем дорога №2, хотя имеет и незначительно больший уклон, но наиболее удобный для проезда рельеф. По ней в случае возникновения чрезвычайной ситуации на р. Соть в районе стратегического пункта можно наиболее оперативно доставить группы спасателей, инструменты, оборудование, питание, технику и так далее. Также по этой дороге можно в более короткие сроки провести эвакуацию населения в н.п. Дворики и Борисово. В любом случае решение выбора подъезда к стратегическому объекту следует принимать исходя из сложившейся обстановки и прогнозов развития ЧС.

7. ИССЛЕДОВАНИЕ ДЕЙСТВИЯ ВЕТРА

7.1 Построение розы ветров

Для построения розы ветров на исследуемом участке используются данные гидрологического ежегодника или водного кадастра по гидропосту, где проводились измерения силы ветра на высоте 10 метров анеморумбометром. Приводятся сведения о распределении ветра по направлению и скорости в целом за период, свободный ото льда. Повторяемость ветра по градациям направления и скорости, выраженная в процентах от и числа наблюдений, определяется по формуле:

,

где n - число случаев повторяемости ветра по данным интервалов скорости и румбов направлений ветра;

- число всех случаев измерений ветра за период.

Данные сводятся в таблицу 7.1.

Таблица 7.1 - Повторяемость направлений ветра по румбам в процентах

Скорость ветра	Повторяемость направлений ветра по румбам, %
	С	СВ	В	ЮВ	Ю	ЮЗ	З	СЗ
1 - 3	6,3	2,9	3,1	3,3	2,6	2,5	2,1	1,8
4 - 5	3,8	1,2	1,9	2,1	1,8	1,8	1,4	0,7
6 - 7	1,0	0,1	0,1	0,7	1,3	1,5	0,8	0,3
8 - 9	0,3	-	-	0,1	1,2	1,1	0,8	0,1
10-11	0,2	-	-	-	0,4	0,8	0,3	0,1
12-13	-	-	-	-	0,1	0,3	-	-
14-15	-	-	-	-	-	-	0,1	-
16-17	-	-	-	-	-	0,1	-	-
Сумма	11,6	4,2	5,1		7,4	8,1	6,5	3,0

7.2 Расчет сгонно-нагонных течений

При наличии на карте водохранилища, большого пруда, приустьевого участка произведём расчет волн для этих водных объектов, в нашем случае такие объекты отсутствуют, поэтому расчет выполняем в соответствии с данными в приложении. Исходные данные для выполнения задания берём по вариантам в приложении № 6, из методического пособия.

Порядок выполнения задания.

Вычисляются безразмерные величины:

1)

2)

3)

где g= 9,81 м/с² - ускорение силы тяжести;

t= 6ч =21600с-продолжительность действия ветра;

D=22 км=22000 м - длина разгона волны;

W=10 м/с - скорость ветра на высоте 10 м над поверхностью воды;

Т=17 м - глубина.

1)

2)

3)

По вычисленным значениям (1), (2), (3) с помощью номограммы (Приложение №6) из методического пособия находим величины g/W² и g/W. В нашем случае левее значение gD/W², это значит, что волнение установившееся и его параметры определяем в зависимости от длины разгона волны, т.е. по величине gD/W². Находим величины g/W²=0,06 и g/W=3,2 , из которых определяем среднюю высоту волны в метрах и средний период волны в секундах.

=0,61 м,

=3,26 с.

Среднюю длину волны вычисляем по формуле:

м;

Определяем высоту волны обеспеченностью р % в системе волн по формуле:

, м,

где К_р - коэффициент, для определения которого используем (приложение) и предварительно вычисленные безразмерные величины gD/W²и gT/W², а также значение обеспеченности Р, которое в наше случае равно Р=0,1% . По номограмме находим соответствующие им значения коэффициента, которые соответственно равны 2,5 и 2,7, из которых за расчетное принимаем меньшее по величине, то есть К_р=2,5.

м.

Анализируя полученные данные можно сделать вывод:

Наш водный объект имеет установившееся волнение со средней длиной волны 16,6 м и средней высотой 0,61 м при 3,26 с. При обеспеченности равной 2,5% высота волны достигает 1,5 м. На нашем участке местности преобладают ветра северного направления. Наката волны на стратегический пункт от действия ветра не будет.

ЗАКЛЮЧЕНИЕ

В данной курсовой работе произведен анализ бассейна реки Соть по топографической карте. В курсовой работе расчеты сопровождаются рисунками.

В ходе исследования определены следующие характеристики речного бассейна:

1) Общая площадь бассейна- 89,05 км²

2)Средняя ширина бассейна - 7,12 км;

3) Самая высокая точка участка имеет высоту Z_max=259,4 метра;

Самая низкая точка участка имеет высоту Z_min=108 метров;

Перепад высот :151,4м;

4) Географическое положение бассейна:

ц: от 54є00'с.ш. до 54є19'с.ш;

л: от 18є50'в.д до 18є5'в.д ;

5) Длина р. Соть 7,25 км;

6) Длина бассейна- 12,5 км;

7)Уклон водной поверхности при проектном уровне воды 0,23;

8) Средняя глубина 2,2 метра;

9)Средняя скорость течения в поперечном сечении русла 0,83 м/с;

10) Общий расход воды 180,12 м/с;³

11) Общий расход взвешенных наносов 3,98 кг/ с;

12) Средняя мутность в сечении 19,99 г/м;³

Стратегическим объектом на данном участке местности является населенный пункт Ивановка.

На основании проведенных расчетов и анализа ситуации сделаны следующие выводы:

1) Во время весеннего половодья произойдет затопление прибрежных территорий. Населенный пункт Ивановка в зону затопления не попадает.

2) Русло реки не подвергается размыванию.

3) Наиболее оптимальным подъездным путем к стратегическому объекту в случае возникновения там чрезвычайной ситуации является автомобильная дорога №1 (Дворики-Ивановка).

4) На участке преобладают северные ветра.

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

СПИСОК ЛИТЕРАТУРЫ

1. Шамова В.В. Методические указания по курсовой работе по дисциплине “Русловые изыскания” для студентов ГТФ специальности «Защита в чрезвычайных ситуациях»: - Новосибирск: Новосиб. гос. акад. водн. трансп., 2007.

Размещено на Allbest.ru