CC BY
28
grafiya /; Federanoe gos. Unitarnoe pred-priyatie "Kaspijskij nauch.-issled. in-t rybnogo khoz-va" (FGUP "KaspNIRKH"). - Astrakhan: FGUP "KaspNIRKH", 2014 g, S. 56-99.
7. Chugunov N.L. Biologiya molodi pro-myslovykh ryb Volgo-Kaspijskogo rajona. / LGr. Astrakhanskoj nauchnoj Rybokhozyajstvennoj Stantsii. T. VI, vyp. 4, 1928, 282 s.
8. Koblitskaya A.F. Izuchenie neresta presnovodnykh ryb. - M.: Pishchevaya pro-myshlennost. 1966. - 109 s.
9. Alekhina R.P., Kurapov A.A., Fina-eva V.G. Rezhim vesennikh popuskov i ego vliyanie na rannij ontogenez poluprokhodnykh ryb Volgo-Kaspijskogo rajona // 4-ya Vsesoyuz. konf. Po rannemu ontogenezu ryb. Tez. dokl. -Murmansk, 1988. - S. 16-18.
10. Koblitskaya A.F. Vliyanie dlitelno-go zaregulirovaniya stoka reki i kolebanij urovnya Kaspijskogo morya na estestvennoe razmnozhenie promyslovykh ryb v ustjevoj oblasti Volgi // Ekologiya molodi i problem
vosproizvodstva kaspijskikh ryb: sb. nauch. Trudov FGUP «KaspNIRKH». M.: VNIRO, 2001. S. 126-138.
The material was received at the editorial office
04.06.2018 g.
Information about the authors Zhezmer Valentin Borisovich, candidate of agricultural sciences, the leading researcher, head of the department of safety of hydraulic structures of irrigation and drainage complex (FGBNU «VNIIGiM named after A.N. Kostyakov»). Tel / Fax: +7(499)1532133, e-mail: v1532133@yandex.ru
Buber Alexander Leonidovich, Leading Researcher, Head of the Department of Melioration and Water Management Complex (FGBNU VNIIGIM named after A.N. Kostyakov), 127550, Moscow, ul. Bolshaya Akademi-cheskaya, 44, building. 2; tel.: +7(499)1531678; e-mail: buber49@yandex.ru, buber@vniigim.ru
УДК 502/504:532.592.7 DOI 10.26897/1997-6011/2018-5-30-35
О.Я. МАСЛИКОВА, В.К. ДЕБОЛЬСКИй
Институт водных проблем РАН, г. Москва, Российская Федерация
ВОЗДЕЙСТВИЕ СУДОВЫХ ВОЛН НА ДИНАМИКУ БЕРЕГОВОГО СКЛОНА
Воздействие судовых волн на берег, сложенный рыхлыми песчаными отложениями, обусловливает деформации подводного склона. Перенос размытого грунта волновым потоком в продольном и поперечном направлениях приводит к значительным изменениям морфологических характеристик русла и, в результате, к уменьшению глубин на судовом ходу. С другой стороны, для полноценного проектирования берегозащитных работ требуется знание поведения береговой линии в естественных условиях, установление взаимосвязи между волнами, подходящими к берегу и характеристиками произведенного профиля берега. Предлагается аналитический метод расчета, позволяющий оценить ряд характеристик, достаточно полно отражающий условия динамики прибрежной зоны: параметры поля высот волн, средний волновой уровень, скорость вдольбереговых течений, продольные расходы песчаных наносов с целью выявления оптимальных параметров судового движения. Рассматривается возникновение волн под действием движущихся судов, рассчитывается амплитуда и угол подхода волн к берегу. В результате, зная состав пород, слагающих берег, можно предложить максимальные скорости движения судов как по течению, так и против течения, при которых профиль равновесия берегового склона останется ненарушенным. Можно решать и обратную задачу: зная необходимые скорости движения водного транспорта, создавать искусственные вдольбереговые насыпи, наименее чувствительные к данным воздействиям.
Геофизика, прибрежная зона, русловые процессы, физические процессы, деформации
профиля дна.
Введение. При проектировании стро- при реконструкции существующих в числе ящихся гидротехнических объектов, а также прочего следует учитывать влияние судовых
волн, которые, подходя к берегам, разрушают их. Навигационный транспорт генерирует определенный вид береговых волн, которые имеют особое влияние на топографические границы водного объекта и локальные гидравлические условия у берега посредством более выраженного по сравнению с ветровыми волнами воздействия. Характеристики этого воздействия, такие как: направления распространения, высоты и периоды волн определяются характеристиками судна (разные типы, тяга, скорость)и направлением их движения. Разрушительная сила судовых волн наиболее полно проявляется на берегах, сложенных рыхлыми песчаными отложениями. При приближении судна, идущего полным ходом, уровень воды возле берега заметно повышается. Как только судно поравняется с берегом, начинаются отлив воды и резкое понижение ее уровня. После того как судно пройдет мимо, на берег вкатывается большая судовая волна, иногда с пенистым гребнем, что может вызвать большие разрушения. При подходе одиночного судна по течению или судна с составом влияние судовых волн на берега и стоящие суда оказывается несколько меньшим. Чтобы оценить разрушительную силу судовых волн и их воздействие на береговой размыв, надо выявить основные факторы, влияющие на возникновение и распространение судового волнения. В настоящей работе объединены теоретические исследования, касающиеся условий возникновения и распространения судовых волн, в том числе с учетом угла подхода к берегу, исследования существования и сохранения профиля динамического равновесия берега, а также теория транспорта наносов в потоке.
Результаты и их обсуждение. Параметры судовых волн. Так как судовые волны приблизительно можно считать плоскими, то высота волны, исходящей от судна в момент возникновения, равна высоте волны, подходящей к прибрежной зоне. Можно рассчитать все параметры рассматриваемых волн, используя заданные параметры судна. Значения высоты КеН, м, и длины Ан, м, волн, возникающих при движении водоизмещающих судов, следует определять по формулам [1]:
Кн =
2V2 Ш
ё
и '
V2
Ан = 5,65 V- (1)
учитывать, что скорость судна не должна превышать допускаемое значение определенное нормативными актами для рассматриваемого водного пути. При отсутствии нормативных данных Уш1т обычно считают по формуле:
уаат = 0,9.
п + агссов (1 - ка) . , . беев-^-'а! - 2 (1-ка)
4 (2)
Ь
где k — отношение подводной площади поперечного сечения судна к площади живого сечения канала А, м; Ь — ширина канала, м, по урезу воды. При двухстороннем движении однотипных судов численное значение k в формуле (2) следует принимать увеличенным вдвое.
Высоту наката Н н, м, судовой волны на откос (рисунок 1) следует определять по формуле
0,5кн + 0,05со^
(при У < Ут): =в
V2 ё
1 - 0,05со^
(3)
где Рв1 — коэффициент (для откосов, облицованных сплошными плитами, равным 1,4, каменным мощением — 1,0 и каменнойнаброской — 0,8).
При отсутствии облицовки волну наката можно считать [2]
ь =
ё1апЬТ2 2п
(4)
Т — период волн, Ь — угол заложения откоса.
Рис. 1. Эпюра давления судовых волн на крепления берегов каналов при накате волны на откос
Транспорт наносов. Количество материала, которое будет эродировано, или транспортирующая способность потока, зависит от средней по глубине скорости течения воды, глубины потока и гидравлической крупности частиц [2]:
где ( и 1и - осадка и длина судна, м; V - скорость судна, м/с; 8 - коэффициент полноты водоизмещения судна. При определении параметров волны, возникающей при движении во-доизмещающих судов в водотоках (реках, каналах), следует
и3
(5)
где 8 — транспортирующая способность потока; и - средняя скорость потока; Н - глубина потока; w - гидравлическая крупность материала, слагающего берег.
Скорости противотечения, действующие на откосе при глубине к, определяются соотношением:
U =
5п2 H2 4TLsh2kh
(6)
Н — высота волн, Ь — длина наката волны
Умножая выражение (5) на время прохождения судна, которое равно отношению длины судна к его скорости, получаем объем размыва от воздействия одного судна.
Угол распространения волн за судном. От носа и кормы образуются самостоятельные волновые системы, которые интерферируют друг с другом. По [3]:
r =1 V t cos 6> 2 c
(7)
где V — скорость судна; t — текущее время; r — путь, пройденный волной от судна до точки интерференции на берегу; y = rsind — расстояние от судна до берега
(или половина ширины канала), или —y— = 1 ct c0Sí?,
sin в 2
откуда угол, под которым распространяются судовые волны в момент возникновения:
„ 1 ■ (4 У в = — arcsin I —
2 1 ct
(8)
Рис. 2. Волны судовые:
1 - расходящиеся; 2 - поперечные
Фото с сайта http://udivitelno.com/mesta/item/396-gejranger-ford-norvegija-foto
Область допустимых значений арксинуса:
1 4у 1 4у тх
-1 < — < 1, или t > —. Именно начиная с этого
^ с
времени (не раньше) волны от судна начинают достигать берега.
В [3] геометрически показано, что волны чувствительной (интерферирующей друг с другом амплитуды) будут образовываться только в угле раствором 38°56', биссектрисой которого является ось движения судна. Значит, в обе стороны и к берегу волны будут распространяться под углом 19°28'. На рис. 2 показано, что через любую точку внутри этого раствора проходит два вида кривых одинаковой фазы.
Влияние течения реки на распространение судовых волн. В [3] было доказано, что скорость распространения образующихся при движении корабля поперечных волн равна скорости самого корабля. Посчитаем скорость распространения V диагональных (расходящихся) волн (рис. 3). У-компонента этих волн будет равна Vс — скорости судна. Х-компонента диагональной волны в момент возникновения равна соответственно Vcx0 = Vtg19028' = 0,35Ус. При подходе к берегу эта компонента будет изменяться с изменением скорости течения реки.
Рис. 3. Распространение диагональных волн
Определим поверхностную скорость течения реки в виде:
U y = Ur 311 -
( V
У
V У max J
(9)
Здесь ик — максимальная скорость течения реки, утах — расстояние от стрежня до берега. Линия стрежня совпадает с траекторией судна.
Координату у = к , начала обруше-tan г
ния волн определим из [Лонге-Хиг-гинс] Гв = = 0,8.
h
У =
H,
sh
0,8 tan i
(10)
где Нак — высота подходящих волн, в нашем случае — судовых, I — береговой уклон. Значит, х-компо-нента пришедшей к берегу диагональной волны будет
V = V . + и(у), у-компонента останется неизменной,
ох охи / > ^ ?
или равной скорости судна. Таким образом, в зависимости от взаимного направления течения и движения судна, х-компонента рассчитывается как:
Vrr = 0,3SVp ± U(y'),
(11)
Угол подхода волн к прибойной зоне:
V
tanp = -^-
(12)
впадают, скорость максимальна, а там, где они противоположны — минимальна. Если
в и в — характерные уклоны дна по нормах у
ли и вдоль берега, а ф - угол подхода волн, то, как показано И.О. Леонтьевым [6], энергетическое течение доминирует над градиентным при условии
srn^ >
1 в 2 в
Таким образом, система уравнений (9-12) описывает расчет угла подхода судовых волн к берегу, который зависит от уклона дна в прибойной зоне, ширины реки, скорости течения реки и скорости судна. При приближении к берегу х-компонента скорости диагональной волны уменьшается вместе с поверхностным течением при неизменной у-компоненте, что изменяет угол подхода волны. На широких реках этот угол меньше, чем на узких. Данный параметр является важным для расчета динамики берегового профиля равновесия.
Вдольбереговые течения. В динамике береговой (или прибрежной) зоны, определяемой как зона взаимодействия волн с дном, главную роль играют волновые течения. Они возникают отчасти благодаря стоксову переносу масс, но главным образом как результат процессов диссипации энергии и потери импульса в волнах, распространяющихся к берегу, и образуют в совокупности систему прибрежных течений (nearshorecurrentsystem- в англоязычной литературе).
Ее элементами являются как поперечные по отношению к изобатам, так и продольные (вдольбереговые) течения. К поперечным могут быть причислены потоки, связанные со стоксовым переносом и направленные к берегу, а также компенсационное противотечение (undertow), действующее в противоположном направлении. Сюда же относятся и разрывные течения (ripcurrents) в виде узких струй, отходящих от берега. Вдольбереговые течения условно подразделяются на энергетические, генерируемые при косом подходе волн, и градиентные, обусловленные изменениями среднего уровня [4]. При наблюдениях в природных условиях бывает довольно трудно идентифицировать тот или иной вид течений, поскольку все они, как правило, действуют одновременно [5]. На участках, где направления энергетического и градиентного течений со-
Вдольбереговое энергетическое течение, генерируемое косоподходящим волнением, служит одним из важнейших механизмов адвекции растворенного и взвешенного вещества в прибрежной зоне и играет главную роль в формировании потока наносов вдоль берега. Протяженность течения зависит от характера контура береговой линии, ширина определяется масштабом ширины прибойной зоны, в которой течение зарождается и откуда оно диффундирует в соседние области.
М. Лонге-Хиггинс [7] нашел теоретическое распределение скоростей течения для регулярных волн, распространяющихся под небольшим углом к берегу на мелкой воде над плоским наклонным дном. Коэффициент горизонтального обмена (вихревой вязкости) в прибойной зоне с монотонно повышающимся дном определялся с помощью соотношения и1 ~ ёНТ, где Ни Т - высота и период волн.
Масштаб скорости на линии обрушения при подходе волн под углом ф [5]:
Vo = ТЛ si
16' C
sin (p
(13)
H
С. - коэффициент трения, ув = —, Н — полная глубина,
f п
включающая и отклонение среднего уровня от штилевой отметки, в - коэффициент порядка единицы.
Соотношение показывает, что скорость течения прямо пропорциональна уклону дна, определяющему скорость диссипации, синусу угла подхода волн и квадратному корню из глубины обрушения, зависящей в основном от исходной высоты волн.
На рисунке показан поперечный профиль вдольберегового течения над плоским наклонным дном по Лонге-Хиггинсу [7]. Параметр Р характеризует отношение вкладов бокового перемешивания и донного трения. При максимальной относительной высоте разрушающихся волн у = 0.8 (критерий разрушения одиночных волн) имеем приблизи-
тельно Р = 4в, у - расстояние от берега, ув ширина прибойной зоны, V - переносная скорость.
торое равно отношению длины судна к его скорости, можно получить объем размыва от воздействия одного судна.
Рис. 4. Поперечный профиль вдольберегового течения над плоским наклонным дном по Лонге-Хиггинсу при различных значениях параметра перемешивания Р. Точки относятся к значениям Р = 0.14-0.16
Профиль динамического равновесия. Уравнение профиля динамического-равновесияпри нормальном подходе волн к берегу [6]:
dh + 4Wf h» 12 - A = О,
dx у C0 ^ h
(14)
где h — локальная глубина, h0 — глубина начала разрушения волн, С0 = ^ gh0 — фазовая скорость на глубокой воде, близкая к фазовой скорости в мелководной прибрежной зоне, у = H/h— отношение высоты волн
dh
к глубине в данной точке,_= tan i — средний уклон,
dx
W — гидравлическая крупность твердых частиц, слага-
С»
ющих дно, — — характеристика подвижности наносов, W
поскольку возможности перемещения твердых частиц увеличиваются как с ростом силы волнения (с ростом hg), так и с уменьшением их крупности.
Величина А включает в себя ряд констант и имеет смысл осредненной по профилю характеристики асимметрии волновых скоростей. Анализ данных [8] показал четкую линейную зависимость, определяемую уравнением:
W
A = 22,8— С»
Подставляя в выражение (14) масштаб скорости на линии обрушения (13) и интегрируя по времени прохождения судна, ко-
Выводы
Исходя из теории волнообразования под воздействием движущихся судов, рассчитывается амплитуда и угол подхода волн к берегу. На основе модели деформации профиля дна в прибрежной зоне, предложенной И.О. Леонтьевым, была выведена зависимость размыва берегового склона от судового волнения. Зная состав пород, слагающих берег, можно предложить максимальные скорости движения судов как по течению, так и против течения, при которых профиль равновесия берегового склона останется ненарушенным. Можно решать и обратную задачу - зная необходимые скорости движения водного транспорта создавать искусственные вдольбереговые насыпи, наименее чувствительные к данным воздействиям.
Интегрирование воздействий этих скоростей во времени позволит определить интенсивность переформирования берегов водных объектов и разработать мероприятия по их укреплению.
Таким образом, предложена теория воздействия судовых волн на береговые насыпи и сооружения, объединяющая описание возникновения и распространения судовых волн, исследования существования и сохранения профиля динамического равновесия берега, а также теорию транспорта наносов в потоке.
Библиографический список
1. СП 38.13330.2012. Нагрузки и воздействия на гидротехнические сооружения (волновые, ледовые и от судов) СНИП строительные нормы и правила. https://docplan. ги/Оа!а2/1/4293793/4293793649.Ы;т
2. Динамика русловых потоков и лито-динамика прибрежной зоны моря. / под ред. Дебольского В.К. - М.: Наука, 1994. - 303 с.
3. Кочин Н.Е., Кибель И.А., Розе Н.В. Теоретическая гидромеханика. Ч. I. Под ред. И.А. Кибеля., 6-е изд., испр. и доп. — М.: Физматлит, 1963. - 583 с.
4. Шадрин И.Ф. Течения береговой зоны бесприливного моря. - М.: Наука. 1972. -128 с.
5. Динамические процессы береговой зоны моря. / Под ред. Р.Д. Косьяна, И.С. Подымова, Н.В. Пыхова). - М.: Научный мир, 2003. - 320 с.
6. Леонтьев И.О. Динамика прибойной зоны. — М.: Ин-т океанологии им. П.П. Ширшова АН СССР, 1989. - 184 с.
7. Longuet-Higgins M.S. Longshore currents generated by obliquely incident waves. Parts 1 and 2. J. of Geophys. Res. 1970. V. 75. № 33. P. 6778-6801.
8. Леонтьев И.О. Профиль динамического равновесия: проверка теории. // Океанология. - 1992. - т. 32, вып. 2. - С. 355-361.
Материал поступил в редакцию 24.06.2018 г.
O.YA. MASLIKOVA, V.K. DEBOLJSKY
Ras institute of water problems, Moscow, Russian Federation
Сведения об авторах
Масликова Оксана Яковлевна, кандидат технических наук, старший научный сотрудник ИВП РАН;119333, г. Москва, ул. Губкина, д. 3; тел.: +7(499)1357201, e-mail: oksana68@mail.ru
Дебольский Владимир Кириллович, доктор технических наук, профессор, заведующий лабораторией ИВП РАН; 119333, г. Москва, ул. Губкина, д. 3; тел.: +7(499)1357201, e-mail: vdebolsky@mail.ru
IMPACT OF SHIP WAVES ON THE DYNAMICS OF COASTAL SLOPE
The impact of ship waves on folded by loose or sand deposits shore, causes deformations of the underwater slope. The transfer of the eroded soil by the wave flow in the longitudinal and transverse direction leads to a significant change in the morphological characteristics of the channel and, as a result, to a decrease in depth on the ship's run. On the other hand, for the rational design of bank protection works, knowledge of the behavior of the shoreline in natural conditions is required, establishing the relationship between the waves approaching the shore and the characteristics of the produced shore profile. The analytical method of calculation is proposed, which makes it possible to evaluate a number of characteristics, parameters of wave completeness, average wave level, speed of alongshore currents, longitudinal expense of sand with tasks of optimal parameters of ship velocities. The appearance of waves under the action of moving vessels is considered, the amplitude and angle of approach of the waves to the shore are calculated. As a result, knowing the composite soil of the shore, it is possible to offer the maximum speed of the vessels' movement along the flow and against the flow, at which the equilibrium profile of the coastal slope will remain permanent. It is also possible to solve the inverse problem - knowing the necessary speed of water transport, create artificial coastal embankments, the least sensitive to these impacts.
Geophysics, coastal zone, channel processes, physical processes, deformation of the bottom profile.
References
1. SP 38.13330.2012. Nagruzki i vozde-jstviya na gidrotehnicheskie sooruzheniya (volnovye, ledovye i ot sudov) SNIP stroi-telnye normy i pravila. https://docplan.ru/ Data2/1/4293793/4293793649.htm
2. Dinamika ruslovyh potokov i litodi-namika pribrezhnoj zony morya. / pod red. Debolskogo V.K. - M.: Nauka, 1994. - 303 s.
3. Kochin N.E., Kibel I.A., Roze N.V. Teoreticheskaya gidromehanika. Ch. I. Pod red. I.A. Kibelya, 6-e izd., ispr. i dop. - M.: Fizmatlit, 1963. - 583 s.
4. Shadrin I.F. Techeniya beregov-oj zony ^eHaa besprilivnogo morya. - M.: Nauka. 1972. - 128 s.
5. Dinamicheskie protsessy beregovoj zony morya. / Pod red. R.D. Kosjana, I.S. Podymova, N.V. Pyhova). - M.: Nauchny mir, 2003. - 320 s.
6. Leontjev I.O. Dinamika pribojnoj zony. - M.: in-t okeanologii im. P.P. Shirshova AN SSSR, 1989. - 184 s.
7. Longuet-Higgins M.S. Longshore currents generated by obliquely incident waves. Parts 1 and 2. J. of Geophys. Res. 1970. V. 75. № 33. P. 6778-6801.
8. Leontjev I.O. Profil dinamicheskogo ravnovesiya: proverka teorii. // Okeanologiya. -1992. - t. 32, vyp. 2. - S. 355-361.
The material was received at the editorial office
24.06.2018 g.
Information about the authors Maslikova Oksana Yakovlevna, candidate of technical sciences, senior researcher IVP RAN; 119333, Moscow, ul. Gubkina, d. 3; tel.: +7(499)1357201, e-mail: oksana68@mail.ru
DebolskyVladimirKirillovich, doctor of technical sciences, professor, head of the chair IVP RAN; 119333, Moscow, ul. Gubkina, d. 3; tel.: +7(499)1357201, e-mail: vdebolsky@mail.ru
