CC BY
66
УДК 004.94+556.56
Ан.А. Лыгин, К.В. Марусин, А.Ш. Хабидов, А.А. Лыгин,
Е.А. Федорова, А.Л. Хомчановский
Расчетные программные компоненты прототипа информационной системы мониторинга и прогнозирования динамики берегов*
An.A. Lygin, K.V. Marusin, A.Sh. Khabidov, A.A. Lygin,
E.A. Fedorova, A.L. Khomchanovsky
Software Components of the Information System Prototype for Coastal Zone Monitoring
Экспериментальный образец информационной системы мониторинга и прогнозирования динамики берегов содержит семь расчетных программных компонентов для решения прогностических задач на основе полученных данных: «Программа расчета параметров ветровых волн вне береговой зоны», «Программа расчета параметров ветровых волн в береговой зоне», «Программа расчета профиля динамического равновесия берегового склона», «Программа расчета вдоль-берегового потока наносов», «Динамика профиля», «Абразия», «Программа расчета изменения положения береговой линии». В статье представлены первые четыре программных компонента.
Ключевые слова: информационная система, прогностическое программное обеспечение, мониторинг береговой зоны, морфология берегов и дна искусственных водоемов, Новосибирское водохранилище.
Введение. Согласно Постановлению Правительства Российской Федерации от 10 апреля 2007 г. № 219 [1] одной из основных целей государственного мониторинга водных объектов, их побережий и береговых зон является информационное обеспечение управления в области использования и охраны морей, озер и водохранилищ. Инструментарием информационного обеспечения должна стать разрабатываемая в рамках Федерального целевого проекта № 16.515.11.5075 информационная система мониторинга состояния берегов естественных и искусственных водоемов.
При выполнении Государственного контракта были изучены и обобщены результаты современных исследований по математическому моделированию и прогнозированию природных процессов, протекающих в береговой зоне водоемов. Реализованные математические модели базируются на эмпирических уравнениях, положенных в основу нормативных документов, действующих на территории Российской Федерации.
Prototype of the Information System for Coastal Zone Monitoring includes seven software components for solving predicted problems based on the obtained data: program for calculating wind wave parameters outside the coastal zone, program for calculating wave shoaling and refraction, program for calculating equilibrium beach profile and nourishment volume, program for calculating the net annual alongshore transport, program for modeling storm-induced sand beach profile changes, program for modeling erosion of the clay coastal profile, program for modeling shoreline changes. The paper presents the first four software components.
Key words: information system, predictive software, coastal zone monitoring, man-made lake shores and bottom morphology, Novosibirsk Reservoir.
Организация сбора, хранения и предоставления информации, получаемой от пространственно-распределенной сети станций Федеральной службы по гидрометеорологии и мониторингу окружающей среды России и из наблюдений, ведущихся с использованием разработанного в рамках обозначенного выше проекта мобильного аппаратного комплекса мониторинга береговой зоны, является отдельной большой темой и в данной статье не рассматривается.
Общие положения реализации расчетных программных компонентов. Графический пользовательский интерфейс приложения, оформленный в виде исполняемого файла, реализован отдельно от динамической библиотеки расчета и организации доступа к данным.
Входные и выходные данные каждого программного продукта объединены в наборы-проекты. Для хранения проектов и данных используются бинарные файлы, полученные средствами сериализации.МБТ. Разновидности всех наборов сохраняемых данных
* Работа выполнена при поддержке РФФИ (проекты 11 -05-00615-а и 11-05-10046-к) и федеральной целевой программы «Исследования и разработки по приоритетным направлениям развития научно-технологического комплекса России на 2007-2013 гг.» (проект 16.515.11.5075). Заказчик — Министерство образования и науки РФ.
оформлены в виде сущностей для удобной организации контрактов пересылки по сетевым протоколам. В качестве стандарта взаимодействия с внешними приложениями (например, Microsoft Excel) выбран XML 1.0. Для удобства пользования реализованы расшифровки обозначений и допустимых значений входных параметров в виде всплывающих подсказок. Язык программирования: С#.
Платформа: .NETFrameworkвер сии4.0.30319.
Среда разработки: SharpDevelop версии 4.2.0.8783.
CC01 — Программа расчета параметров ветровых волн вне береговой зоны/
Задача генерации волнения. Ветер, дующий над поверхностью водоема, генерирует нерегулярные волны, параметры которых в общем случае зависят от длины разгона (расстояния от подветренного берега до расчетной точки), глубины водоема, скорости ветра и времени его действия. Нами рассматривается только установившийся режим волнообразования, в котором время действия ветра не влияет на параметры волн.
Средняя высота и период волны рассчитываются по эмпирическим зависимостям, полученным Ю. М. Крыловым с соавторами [2].
В случае если глубина дна существенно изменяется, путь волны до расчетной точки разбивается на ряд участков с примерно постоянным уклоном
дна. При этом расчет ведется для конечных точек каждого участка. При уменьшении глубин на некотором участке разгон к концу этого участка корректируется с учетом так называемого условного разгона. В случае превышения отношения уменьшения глубин на концах некоторого участка к длине этого участка значения 0,002 возникают трансформации волн, вызванные в большей степени влиянием изменения глубин, а не ветра (следующая расчетная задача). Кроме того, высота волны Н в каждой расчетной точке должна проверяться по критерию обрушения
H / h < b.
(l)
где Ь — коэффициент обрушения, обычно принимаемый 0,78; h — глубина.
Входными данными задачи является массив характеристик последовательных участков (участка) разгона — длины и глубины в конце участка и набор желаемых скоростей ветра. Пользователь также имеет возможность изменять коэффициент обрушения.
Выходными данными программы являются значения статистических (средних, среднеквадратических и значимых) параметров волн (высоты и периода) в точке обрушения и во всех конечных точках допустимых участков разгона. Пример работы программы СС01 представлен на рисунке 1.
СС01 - Программа расчета параметров ветровых волн вне береговой зоны
Проект Данные О программе - Область данных
ІІ hi
: ) 11.0
S30.0 U.O
EEC.5 13.11
vi
7. С
-.4
5.1
3.0
16.0
11.с
Параметры волн в конечной точке
Расчет параметров
vi 1 н т Hrms Hsig Тзід Тр
и. оаз 1,25 0.1 С С 0.V2 1.51 1.67
5.1 D.ICE 1.35 С.11Э 0.163 1.Є2 1.30
7.0 СЛЕЗ 1.57 0.173 С.2^5 1.35 2,05
3.0 3,173 | г 3.2С1 С.23Е 2.С1 2.2?
11.С С,25^ 1J3 3.237 .. 37 2.?2 2.57
16.0 H.3S2 2,27 3/31 3.611 2.72 3,31
Фактический разгон, м
Параметры волн на участках разгона
vi ІІ hi X н Т Hrms Hsig Тзід Тр -
7,0 Е5~ 13.0 1 75 С С. 15? 1.57 С,173 С.2^5 1.39 2,0S
4,4 Е5:" 13.0 1 75 С С, 388 1.25 С. 100 С. 142 1.51 1,67
5.1 Е5:" 13.0 1 75 С С. IDE 1.35 С. 11В C.16S 1.Е2 1,30
3.0 Е5:" 13.0 1 75 С С. 173 1 ЕЗ D.201 0.2 ЗЕ 2.01 2.2?
16.0 Е5~ 13.0 1 75 С С.332 2.27 С. 431 D.E11 2.72 3,01
11.С Е5~ 13.0 1 75 С С .25^ 1.33 С.237 0.4 С 7 2.32 2.5“
7.0 ЕСС 14.D 1 130 С.125 1.33 с.нг 0.201 1.Е7 1.35
4,4 ЕСС К.С 1 133 С, С 73 1.12 OLOB3 С.113 1.3^ 1,45
5.1 ЕСС К.С 1 133 С.С57 1.2 С D.C'55 C.14D 1.44 1,60
3,0 ЕСС КС 1 133 C.U5 1/7 С. 164 с.;:: 1.77 1,56 -
Рис. 1. Пример работы программы СС01; параметры волн в конечной точке отсортированы по возрастанию скоростей ветра, параметры волн на участках разгона отсортированы по убыванию фактического разгона волн к концу участков
l07
В программе СС01 предусмотрена проверка введенного значения параметра по типу и значению («на лету», по событию) и полная проверка всех входных данных перед расчетом или по команде пользователя.
^02 — Программа расчета параметров ветровых волн в береговой зоне.
Трансформация и рефракция волн. Классическая задача в рамках линейной теории волн [3]. В береговой зоне, двигаясь над уменьшающимися глубинами, волны испытывают процессы трансформации — изменения их высоты и рефракции — изменения угла подхода к берегу, т. е. угла между фронтом волны и линией берега (или между лучом волны и нормалью к берегу). При достижении критической высоты волна начинает обрушаться и далее распространяется в виде буруна или прибоя. В целом с приближением к берегу высота волн возрастает, а угол подхода уменьшается, периоды волны остаются практически постоянными.
В инженерных расчетах параметров волн в береговой зоне полагают, что изобаты на данном участке береговой зоны прямолинейны и параллельны, волнение является регулярным и описывается линейной теорией волн, донное трение весьма мало. Также пренебрегают генерирующим действием ветра ввиду относительной малости глубин и протяжения береговой зоны по сравнению с остальной долей разгона.
С учетом этих предположений и допущений справедливы следующие законы сохранения:
— закон постоянства периода волны;
CC03 — Программа расчета профиля динамического равновесия берегового склона
Программа предназначена для расчета профиля динамического равновесия подводного берегового склона и объемов отсыпки (питания) искусственно-
— закон сохранения потока энергии между двумя волновыми лучами (линиями нормальными к фронту волны);
— закон рефракции Снелла.
В линейной теории волн высота волны Н никак не связана с ее длиной L, т. е. могут существовать волны любой крутизны Н/Ь, что, конечно, физически неверно. Поэтому начальные параметры волнения предварительно оцениваются по критерию максимальной крутизны (нелинейная теория Стокса 1-го порядка).
По известным параметрам волнения в некоторой начальной точке, используя законы сохранения и дисперсионное уравнение, решенное методом аппроксимации Ханта [4], находятся следующие интересующие нас параметры волн в расчетных точках:
— высота волны;
— длина волны или период волны;
— угол подхода волны;
— фазовая скорость волны.
Расчет ведется в точках (значениях глубин), заполненных автоматически с некоторым настраиваемым шагом либо, используя заданные пользователем глубины, до достижения критической высоты волны, определяемой через равенство левой и правой части соотношения (1) и в самой точке обрушения.
Предусмотрено графическое представление зависимостей обозначенных параметров от глубины при подходе к берегу. Пример работы программы СС02 представлен на рисунке 2.
го пляжа, необходимых для компенсации дефицита наносов в целях защиты морских берегов и берегов крупных внутренних водоемов.
Энергия волн при обрушении частично рассеивается и частично расходуется на перемещение слагающе-
И СС02 - Программа расчета параметров ветровых волн в береговой зоне
Проект Данные
О программе
Параметры волн в начальной точке
Глубина.м 10,4 “ і
Высота волны, м 2.2С - і
Период волны, с 3.15 “ :! г
Длина волны, м 1Є.&4 “
Угол подхода волны. град 7&.75 -
Коэффициент обрушения 0.8Є - 1
Расчет параметров
Точки расчета
N
1:7 і
Е,С"
:.сс
Заполнять автоматически с шагом. м
17 | с.1 ±|
Параметры воли на расчетных глубинах
ы Н Т і_ А С ^1
10/с 2,2и '■.21 1Є.С4 іг~і '.СС □
1С.?С 2.13 121 1Є.СЗ 72.51 Е.""
1 С.2С 2.1 В ;.21 1Е.СЗ 78/3 Е.ии
1 С.1 с 2.17 121 1Е.С2 73/7 Е.""
1..., 2.17 ;.21 1Е.С2 7В/4 5.ии
э.Э" 2.17 121 1Е.С2 73/2 Е.СС
5.30 2.17 ;.21 1Е.С2 73.?; Е.ии
э.7" 2.1 Є 121 1Е.С2 ■;.?£ Е.СС
Е.ЕС 2.1Є 3.21 1Е.С2 7;.?3 Е.ии
Построить график
Рис. 2. Пример работы программы СС02; в области результатов выделены столбцы высот волн (Н) и скоростей волн (С)
для построения их зависимостей от глубины
го дно несвязного (рыхлого) материала. При условно постоянном во временном плане внешнем воздействии (скорость диссипации энергии волн на единицу объема водной толщи постоянна) профиль подводного берегового склона стремится к состоянию динамического равновесного [5], по достижении которого никаких дальнейших изменений как по профилю, так и в плане, в среднем, не происходит. Динамически равновесный профиль не остается неизменным. Под действием текущих волнений материал все время перемещается по склону к берегу и от берега, но не покидает профиль. Форма склона на профиле при динамически равновесном состоянии аналитически описывается функцией Дина [5]. Параметр крутизны профиля, зависящий от медианной крупности материала, определяется по соотношениям, представленным в [6].
Для расчета параметров необходимы следующие исходные данные:
— таблица координат исходного профиля дна;
СС04 — Программа расчета вдольберегового потока наноса. Программа предназначена для расчета потоков наносов (по направлениям, за период и результирующего).
Подходящие под некоторым углом к берегу волны, обрушаясь, создают вдольбереговой дрейф и, следовательно, перенос материала (транспорт наносов). Полный (через весь створ береговой зоны) расход наносов зависит от высоты, периода, угла подхода волны в точке ее обрушения и крупности донного материала.
В данной расчетной задаче для песчаных материалов используется формула Леонтьева [7]. Для пес-
— линейные параметры бермы;
— набор значений медианной крупности материала;
— высота расчетной волны;
— набор координат расчетных точек вдоль поперечной к береговой линии оси (линии распространения волны) или шаг приращения координат;
— коэффициент обрушения;
— угол естественного откоса грунта.
Результатом работы программы является таблица координат естественного профиля, динамического равновесного состояния рельефа дна и их разности в расчетных точках. Объемы отсыпки материала на погонный метр береговой линии для каждого значения медианной крупности материала сводятся в отдельную таблицу. Пример работы программы СС03 представлен на рисунке 3. Дополнительно предусмотрена прорисовка профиля с возможностью экспорта изображения в растровом формате (рис. 4).
чаных наносов в интервале медианных крупностей от 0,01 до 1 мм дополнительно предусмотрено использование формулы CERC [8]. В случае переноса гравийных и галечных материалов применяется формула Камфиуса [9]. Для галечных наносов значение коэффициента емкости потока принимают равным 0,4 [10].
Входными данными программы являются:
— экспозиция береговой линии — угол между направлением на север и внутренней нормалью к береговой линии, отсчитываемый по часовой стрелке;
— медианная крупность наносов;
— средний уклон дна в прибрежной зоне;
^ССОЗ: Программа расчета профиля динамического равновесия берегового склона
ЛПІх]
Проект
Данные
О программе
“Берма — Выдвижение, м
|!-С ±1
Высота, м
Высота волны, м
3.26
Расчет
Нарисовать
профиль
- Исходный профиль —
Коэффициент і г ——тп обрушения I - - ^11
Глубина замыкания, м
/І л
<с С. 1С
8.0 С,40
2С.С 1 .ЕС
4С.І г.ес
ЕС.С з.сс
зс.с 4, ЕС
_ Материал -
Угол естественного откоса.град
[Ттс 3
СІ50
1.0С
2. СО
С.ЕС
“Точки расчета
Ы
11. С С
Заполнить автоматически с шагом, м
11.00 3
“Сопоставление профилей ■
СІ50 у. ТАе 2л ТА е-Тп
1. С С 11. С С О.ЗЕЕ С.Е7Е с.зю
2, С С 11. с с С/ЕЕ С.Е7Е 0.21 Э
0. с 11. с с 0,310 ІС.Е7Е С.ЗЕЕ
■Объемы отсыпки искусственного профиля ■
(150 хйс V ЧЬ ч* 46
1.0С Е2.3С Е/33 Н.2ІЕ 3,313 с.ссс
2. С С ^ 7.2 С 13.23Е и .2^6 -Е.С40 с.ссс
С.ЕС 77.1 С -17.Е72 и,2^Е 31.І7Й с.з^с
Рис. 3. Пример работы программы СС03
— начальная расчетная глубина на входе в береговую зону;
— коэффициент обрушения;
— коэффициент пористости материала;
— плотность материала;
— плотность воды;
— таблица из строк параметров волнения в начальной точке для каждого шторма с указанием продолжительности воздействия.
Результаты расчета представлены в таблице суммарных значений расхода наносов за период, распределенных по штормам, с возможностью сортировки вывода по объему перенесенного материала. Результирующее значение потоков (влево, вправо относительно наблюдателя на берегу и итогового) отображается в отдельной области. Пример работы программы СС04 представлен на рисунке 5.
•^ССОЗ: Программа расчета профиля динамического равновесия
Диаграмма профилей 560 х 297
Проект
Данные
Точки расчета п і Сопоставление профилей
Берма
Выдвижение, м
Ця ^
Высота, м
ГйГЭ
Высота і---------л
волны, м |2.00 З]
Коэффициент |—------~7]
обрушени я 10. / 8 з]
Глубина замыкания, м
3.2Є
Расчет
Нарисовать
профиль
сі 50
1,00
2,00
о.єо
(15С хОс V № Ур
— 1.СС е;.зс Е/33 8.313 ::::
;.сс ±7 2" 15.205 \ -2--І -ем)
С'.ЄО 77.1 С -17.572 и .5^6 31.-73 0.3^1
Рис. 4. Диаграмма сопоставления профилей с указанием ключевых точек
СС04 - Программа расчета вдольберегового потока наносов
Проект Данные О программе Входные параметры
Экспозиция береговой линии [шнГ"3
.__________. Мешанная крупность I О.ДО
Глубина в начальной точке 110.00 ’
КоэФФинт пористости 10.40
Плотность 12650
Уклон дна
Коэффициент обрушения Плотность воды
Суммарные значения расхода наносов [куб. м за период ]
158734.В <-х-> |45546.3
Результирующий
-13Ш.З
Расчет
Шторма Шторма: расходы наносов и параметры волн
Напр. НІ) тс Р Напр. р о НС тс ІіЬ нь АЬ
22.; 1Л 11.4Є |Продол -кит 5ЛЬ ■юсть што зма,час| ^4В1В,0 1/1 11.4Є 2.37 1.ВЕ 11.19
-І5.С 1.ЕС 14, [С 4ё,0 -■-її. 43.0 -и 2Е 1.В 1,5С 14, ОС 2.31 1,81 -2иМ
-Е7.5 1.ЁС 15.СС 61Л -Е7 .5 61. С -12 ЭВС'.Э 1,60 15.СС 1.53 1,51 -2І.ЕЙ
І-.С 2.СС 12. С С 23,о І-.С 23.с 13 31Є.Є 2.сс 12.СО 2.91 2.27 23.25
1.41) і;.:; ^7,0 -^Е.Э І 7. С -11 ііЗ.1 1/С 13.СС 2.17 1.7С -19.58
-22.5 1.2 С 12.CC 5Е.С -22.5 ЕЁ.С -Е £Сі/ Л 2" 12.СС 2.0В 1.ЕЗ -1 С/Е
Рис. 5. Пример работы программы СС04
Заключение. Расчетные программные компоненты информационной системы мониторинга и прогнозирования динамики берегов, используемые для анализа накапливаемых данных, позволяют составлять как оперативные, так и средне- и долгосрочные прогнозы развития участков береговой зоны в целях защиты морских берегов и берегов крупных внутренних водоемов.
Система мониторинга и прогнозирования, испытанная на Новосибирском водохранилище, мо-
жет стать прототипом других аналогичных систем, намеченных к развертыванию на береговых зонах морей, озер и искусственных внутренних водоемах.
Программные компоненты СС01, СС02, СС03, СС04 прошли процедуру государственной регистрации и внесены в Реестр программ для ЭВМ под номерами 2012660411, 2012660412, 2012660413, 2012660489 соответственно.
Библиографический список
1. Постановление Правительства Российской Федерации от 10 апреля 2007 г. № 219.
2. Крылов Ю. М., Стрекалов С. С., Цыплухин В. Ф. Ветровые волны и их воздействие на сооружения. — Л., 1976.
3. Dean R. G., Dalrymple R. A. Water wave mechanics for engineers and scientists. — World Scientific, 1992.
4. Hunt J. F. Direct solution of the wave dispersion equation // J. waterway, port, coastal and ocean eng. — 1979. — Vol. 105, №WW4.
5. Dean R. G. Beach nourishment. Theory and practice. — World Scientific, 2002.
6. Hanson H., Kraus N. C. GENESIS: Generalized model for simulating shoreline change. Tech. Report CERC-89-
19. — Coastal Engineering Research Center. US Army Corps of Engineers, 1989.
7. Леонтьев И. О. Прибрежная динамика: волны, течения, потоки наносов. — М., 2001.
8. Coastal Engineering Manual. — US Army Corp. of Engineers, 2000.
9. Kamphuis J. W. Alongshore sediment transport rate // J. of Waterways, Port, Coastal and Ocean Eng. ASCE. — 1991. — V. 117.
10. Van Wellen E., Chadwick A. J., Mason T. A review and assessment of longshore sediment transport equations for coarse-grained beaches // Coastal Eng. — 2000. — V. 40.
