CC BY
20
www.volsu.ru
DOI: https://doi.org/10.15688/mpcm.jvolsu.2017.4.3
УДК 519.6 ББК 22.19 + 20.1
ОПТИМИЗАЦИЯ И СЦЕНАРНО-ИМИТАЦИОННОЕ МОДЕЛИРОВАНИЕ ДИНАМИКИ СТРУКТУРЫ МАЛЫХ
_ ___ _ __о ________о л
ИСКУССТВЕННЫХ РУСЕЛ ПОИМЕННЫХ ТЕРРИТОРИИ 1
Владислав Сергеевич Бакулин
Студент кафедры фундаментальной информатики и оптимального управления, Волгоградский государственный университет gitara_93@mail.ru, fiou@volsu.ru
просп. Университетский, 100, 400062 г. Волгоград, Российская Федерация
Анна Анатольевна Васильченко
Старший преподаватель кафедры фундаментальной информатики и оптимального управления, ^ Волгоградский государственный университет
aa-vasilchenko@mail.ru, fiou@volsu.ru < просп. Университетский, 100, 400062 г. Волгоград, Российская Федерация
о к о Е-S
а га
X
Александр Александрович Воронин
Доктор физико-математических наук, профессор, заведующий кафедрой фундаментальной информатики и оптимального управления, < Волгоградский государственный университет
^ voronin.prof@gmail.com, fiou@volsu.ru
® просп. Университетский, 100, 400062 г. Волгоград, Российская Федерация
т
^ Михаил Алексеевич Харитонов
Младший научный сотрудник кафедры фундаментальной информатики
о
х и оптимального управления,
2 Волгоградский государственный университет
к kharitonov.mihail@gmail.com, fiou@volsu.ru
§ просп. Университетский, 100, 400062 г. Волгоград, Российская Федерация
^ Аннотация. Целью настоящей работы является исследование свойств
® оптимальности локальных искусственных русловых систем пойменных тер-
^ риторий. Предложен метод поиска топологической структуры и площади по-
рэ перечных сечений искусственных русел, максимизирующих суммарную взве-
@ шенную площадь затапливаемой территории. Для адекватности постановки и
решения задачи управления проведено адаптивное зонирование цифровой модели территории поймы на участки с незначительным эффектом межзонного затопления при малых и средних паводках. Оптимальная пропускная способность малых русел разыскивалась с помощью многошагового рекурсивного метода динамического программирования. Представлены результаты численной реализации алгоритма с использованием компьютерных гидродинамических имитаций паводковой динамики и геоинформационным моделированием пространственно распределенных параметров рельефа нескольких участков территории северной части Волго-Ахтубинской поймы с искусственными локальными русловыми системами.
Построен взвешенный ориентированный граф межфакторных взаимодействий, определяющих паводковую динамику северной части Волго-Ахтубинской поймы, при помощи которого проведена формализация сценариев развития ситуации и динамики параметров, непосредственно определяющих оптимальную локальную русловую структуру. Приведены субъективные оценки вероятности сценариев и нечеткие оценки динамики параметров русловой структуры.
Ключевые слова: структурная оптимизация, динамическое программирование, производственная функция, геоинформационное моделирование, гидродинамическое моделирование, системы с потоками.
Введение
Волго-Ахтубинская пойма (ВАП), расположенная в нижнем течении р. Волги, имеет протяженность около 450 км, площадь ее территории превышает 20 тыс. кв. км. На ее территории расположено несколько тысяч малых русел, через которые происходит весеннее паводковое затопление. В северной части ВАП расположены природный парк с ценными представителями флоры и фауны, значительное число населенных пунктов, ведется сельскохозяйственный, лесной и рекреационный бизнес. Жизнь экосистемы ВАП полностью определяется весенним паводковым режимом р. Волги, который с 1959 г. регулируется Волжской ГЭС (ВГЭС) и поэтому существенно отличается от природного. В последние десятилетия наблюдается деградация пойменной экосистемы ВАП вследствие прогрессирующего снижения уровня паводкового затопления [2; 5].
Одной из мер противодействия природной деградации ВАП является повышение эффективности территориального распределения паводковых вод путем строительства каналов с управляемой структурой и поперечным сечением русел. В данной работе предложен метод оптимизации структуры локальных искусственных русловых систем пойменных территорий, а также представлены результаты его численной реализации для нескольких участков территории ВАП при использовании компьютерных гидродинамических имитаций паводковой динамики с геоинформационным моделированием пространственно распределенных параметров русловой системы. Проанализированы варианты динамики оптимальной структуры русел при различных сценариях изменения системных параметров. Использованная в работе цифровая модель рельефа (ЦМР) ВАП и численная гидродинамическая модель (программно-алгоритмический комплекс <^еЬ-ЭКОГИС») описаны в [3;6-9].
1. Системный анализ и имитационное моделирование динамики паводковых вод
на территории северной части ВАП
Анализ космических снимков и результатов гидродинамических расчетов [6; 7] показал, что площадь затопления территории определяется величиной постоянного расхода 25 и длительностью Л первой фазы паводка, средние значения которых за последние 30 лет равны соответственно 25 000мъ/с и 10 суток. Карта территории паводкового затопления северной части ВАП при Ц = 25 000 м^/с и £ =10 суток представлена на рисунке 1.
Рис. 1. Распределение воды на территории ВАП при Q = 25 000 м^/с
и £ = 10 суток
Топологическая структура русел северной части ВАП представима в виде бинарного дерева с многими десятками уровней ветвления. Однако существенное различие в величине и характере расходов воды в ветвях этого дерева обусловливают представление структуры с учетом способа ее функционирования и существа решаемой задачи управления. С учетом этого будем считать, что первый (высший) уровень иерархической гидрологической структуры территории северной части ВАП образуют три магистральных водотока (р. Волга, р. Ахтуба и ерик «Пахотный»). Второй — связанные с ними 58 русел глубиной 2-2,5 м, образующих локальные русловые системы. Третий уровень образуют русла глубиной 1,5-2 м, связанные с руслами второго уровня. Русла четвертого уровня меньшей глубины имеют небольшое число локальных русловых систем.
Такая гидрологическая структура создает возможность разбиения территории на зоны преимущественно независимого затопления локальными русловыми системами или даже отдельными ее руслами при малых паводках или в их начальных фазах. Мера погрешности такого затопления, растущая вместе с ростом числа зон, характеризуя потенциальную управляемость процесса паводкового затопления, ставит тем самым предел сложности системе пространственно распределенного управления.
Границы зон, минимизирующие погрешности межзонного затопления при заданных параметрах паводка, находились в ходе гидродинамического и геоинформационного моделирования. Для построения нулевого приближения границ зон на основе иерархиче-
ской цифровой модели русловой структуры использовался модифицированный алгоритм Вороного [10]. Последующее адаптивное «квазинепрерывное» изменение границ проводилось в ходе анализа временного ряда матриц цифровых карт высот водного слоя при = 25 000 мМ/с. В случае превышения предельно допустимой объемной доли трансграничных водных потоков производилось объединение соответствующих зон. Результат моделирования представлен в таблице 1.
Таблица 1
Зависимость параметров зон территории ВАП от величины предельно допустимой объемной доли трансграничных потоков при Q = 25000 м3/с, £ = 5 ■ 105 с
Предельно допустимая объемная доля межзонных расходов паводковых вод (%) Число зон (п) Средняя площадь зоны (км2) Средний расход воды через главное русло зоны (м3/с)
20 54 16,1 55,6
15 41 21,1 73,2
10 32 27,1 93,8
5 28 30,1 107,1
2 17 51,0 176,5
Зависимость числа зон от времени затопления при предельно допустимой объемной доле трансграничных водных потоков 5% представлена в таблице 2.
Таблица 2
Зависимость числа зон от времени затопления при предельно допустимой объемной доле
трансграничных водных потоков 5%
t х 105 с 2 4 6 8 10 12
Число зон 58 49 27 24 19 11
Как показало компьютерное гидродинамическое имитационное моделирование динамики паводковых вод на территории северной части ВАП, непосредственное затопление большей части зон при малых паводках осуществляется из русел третьего иерархического уровня при несущественных трансграничных эффектах. С другой стороны, характер затопления территории при больших паводках характеризуются низкой чувствительностью к структурным параметрам русел. Поэтому управление пропускной способностью русел для максимизации общей площади затапливаемой территории актуально и потенциально осуществимо в большей части зон при значениях параметров первой фазы паводка , £ и меньших, чем их средние величине. При этом общая площадь затопленной территории является суммой (взвешенной суммой) площадей независимых зонных затоплений, расчитываемых по результатам численного гидродинамического моделирования.
Расчет расходов паводковых вод на элементах локальных русловых структур показал примерно равное их распределение между боковыми руслами третьего уровня при их числе от 2 до 4.
2. Задача оптимизации структуры искусственных малых русел пойменной территории с управляемым поперечным сечением
Представленные в п. 1 результаты имитационного моделирования динамики паводковых вод в ВАП позволяют сделать вывод о потенциальной возможности управления их пространственным распределением в некоторой области параметров паводка при адаптивном зонировании территории с учетом локальных особенностей рельефа и русловой структуры. Однако структура локальных русловых систем ВАП, сформированная в условиях природных паводков, может быть неоптимальной в изменившихся за последние десятилетия паводковых условиях. Поэтому представляет интерес поиск зависимости оптимальной структуры локальных русловых систем от параметров паводка и рельефа территории. Для исследования этой проблемы представим следующую формальную задачу оптимизации структуры искусственных малых русел пойменной территории с управляемым поперечным сечением.
Модель пойменной территории (квадрат со стороной Ь) описывается цифровой моделью рельефа, представленной сеточной функцией Ь(г,высот. На территории задана серия (т, п)-разбиений: п-кратных рекурсивных дроблений на т равных частей, дающих в результате тп одинаковых участков (зон). Каждому (т, п)-разбиению территории соответствует структура (дерево) искусственных русел, схематично представленная для нескольких случаев на рисунке 2.
а б в
Рис. 2: а — (3,1)-дерево; б — (3,2)-дерево; в — (3,3)-дерево
Объем паводковых вод У0, поступивших на территорию за время ¿, определяется формулой У0 = , где Я0 — расход воды через входное русло территории. Обозначим у объем вод, локализованных в руслах. Для объемов вод в зонах справедливы балансовые соотношения
т
Уо = У0 - У, У- = £ У-,, д = Щ1 ,к = 1л
3 = 1
Здесь и далее комплексный индекс %\% 2... гп (г к = 1,т, к = 1,п) для краткости записывается как п. Площадь затопленной территории каждой из зон описывается вогнутыми функциями = фгп (Угп), полученными путем кусочно-линейных аппроксимаций временных рядов Б^^к), Уггпаг^к), полученных агрегированием соответствующих участков массива карт затоплений К(¿к), ¿к = 103к, к = 1, 500, полученных в ходе численного гидродинамического моделирования (см. п. 1). С учетом результатов п. 1
расходы воды в русловой системе можно задать формулами = ~0 = 1,т,
к = 1,п). Управлениями пропускной способностью русла с площадью сечения
будем считать доли открытых русел и € [0; 1], так что Б8есН = и8$есН . Тогда
Ягпи3 ^гп-к иу
Q
^ п — к )3
U
Т U
и = ^>г.
г=1
В каждой зоне задан коэффициент относительной экологической ценности а
п.. , (аг3 0,атах).
_а^
Задача управления водными потоками в русловой (т, п)-структуре имеет вид
фг„—1 ( Vin—1) = ап—d Vin—li ( Vin—d)
3=1
—у max
{Vin — li >j = 1 j=l
'У ] Vin — 1,j Vin — 1 , Vin — lj £
0, ^
m"
, Vin-i. (1)
Win — k
(V,
п — к ) = а
п — к 3 *гп — к 3 ( Vin — k 3 ) ^
3 = 1
m
^ vaaxYl Vn—k,3
Vin-k , Vin-k3 £
-kl
3 = 1
2 Vi
го
п-к+1
m
, Угп-к, к = 2, п. (2)
Решение задачи (1)-(2), которое далее находится рекурсивным (по глубине (т,п)-дерева) численным многошаговым (дискретизацией величин объема паводковых вод) методом динамического программирования (рекурсией по ширине (т, п)-дерева), обозначим функциями
Щ-кз( )}, ( Угд_к_1)}, ( Угд_к) = тах^ ( Угд_к3). Оптималь-
ное управление и% 3(Угд—к—1) определяется следующими формулами:
j* = arg max Vi ,V
га — кЗ' гп—1,3
V
u
п 1 ,j
in—1 m
, un-1,j
3 = 1
Uin — 1,j*
1, u
n— 1 ,j
V* ■ _
b- ,J = 1,m. (3)
V
-1,3*
Оптимальные значения т и п находятся методом полного перебора.
Численное решение этой задачи проводилось в рамках серии численных гидродинамических имитаций паводкового затопления пяти ВАП с искусственными (т,п)-системами русел при т = 2, 3, 4 ит =1, 2, 3 и следующими выбранными на основании результатов п. 1 характеристиками параметров паводка и территории: размера (квадрат со стороной Ь = 5 км), Я = 100м3/с, Ь = (0 — 100) х 104с при уклоне русел 0, 03°, и следующих значениях диаметров (¿(п)) и глубин (И(п)) русел: (¿(1) = 50м, Л,(1) = 3м, ¿(2) = 30 м, К(2) = 2 м, ¿(3) = 10 м, К(3) = 1, 5 м. Коэффициенты паводковой неоднородности каждой из тп зон выбирались случайно в диапазоне 1,атпх.
0
и
«3
О
X
о
э!
А в с Б Е
у -
// /У '''
ж у //
// > /X ^
// ^ 1/ #
ж г
1 /
(2;1)
(2; 2) (2;3) (3;2)
0 1 2 3 4 5 6 Уо х 107(м3)
Рис. 3. Динамика затопления одного из участков территории ВАП при различных параметрах русловой структуры. Указан вид оптимальных структур и области их оптимальности {А : (т*,п*) = (2;1), В : (т*,п*) = (2; 2), С : (т*,п*) = (2;3), И : (т*,п*) = (3; 2), Е : (т*,п*) = (2 - 3;2 - 3))
При расчете зависимостей площади затопления зон от объема поступивших в них паводковых вод = (Уп) проводился контроль эффекта межзонного затопления. Если при значении Ь максимальный {по зонам) относительный эффект превышал 5%, то при У > процесс затопления считался неуправляемым.
Результаты численного решения задачи (1)-(3) и поиска оптимальных значений структурных параметров для всех пяти участков территории ВАП незначительно отличаются друг от друга. Типичная динамика затопления одной из территорий при различных параметрах русловой структуры представлена на рисунке 3. Участки А-И соответствуют управляемым, участок Е — неуправляемым сечениям русел и весьма слабой зависимости площади затопленной территории от русловой структуры. Реальной структуре русел данного участка территории ВАП отвечает участок В.
Результаты численного решения задачи (1)-(3) для серии значений объема паводковых вод 1,0 х 107 < У0(ж3) < 7,5 х 107 и переменного параметра паводковой неоднородности рельефа представлены в таблице 3. Ее строки помечены значениями оптимальных параметров структуры т*, п*, столбцы — значениями атах = 3,6. Данные ее первого столбца соответствуют кривым на рисунке 3. Последняя строка таблицы 3 отвечает случаю неуправляемого затопления.
8
6
4
2
0
3. Когнитивный анализ факторов и сценарное моделирование долгосрочной динамики оптимальной структуры искусственных малых русел ВАП
Продолжающаяся деградация природных весенних паводков вместе с социохозяй-ственным освоением и экологической деградацией территории ВАП актуализируют исследование зависимости оптимальной структуры ее локальных русловых систем от долгосрочных тенденций социоприродохозяйственной динамики.
Таблица 3
Значения объемов паводка У0 х 107(м3) , отвечающих параметрам оптимальной локальной русловой структуры т*, п* и паводковой неоднородности рельефа атах = 3, 6
т*; П* \flmax 3 4 5 6
2;1 (A) 1,0-3,0 1,0-2,6 1,0-2,4 1,0-2,1
2;2 (B) 3,0-3,9 2,7-3,8 2,5-3,7 2,2-3,2
2;3 (C) 3,9-4,5 3,8-4,4 3,7-4,3 3,3-3,8
3;2 (D) 4,5-5,5 4,4-5,0 4,3-4,5 3,8-4,2
3;3 — 5,0-5,5 4,5-5,5 4,3-5,5
2-3;2-3 (E) 5,5-7,5 5,5-7,5 5,5-7,5 5,5-7,5
Результаты п. 1 показывают, что оптимальная структура русел с управляемым сечением определяется объемом паводковых вод У и параметром паводковой неоднородности рельефа атпх, поэтому прогноз ее изменения требует дополнения данных таблицы 3 данными прогноза изменения этих параметров.
Анализ результатов гидрологических исследований волжского бассейна [1;4] и комплексного обследования северной части ВАП [2; 5] позволяет формально описать взаимодействие главных факторов ее паводковой динамики взвешенным ориентированным графом, вершины которого имеют следующий содержательный смысл: средний объем паводковых вод, поступающих на территорию; максимальный объем паводковых вод, поступающих на территорию (2); средний объем природного весеннего паводка (3); средний объем весеннего попуска Волжской ГЭС (4); максимальный объем весеннего попуска Волжской ГЭС (5); доля паводковых вод, поступающих в р. Ахтуба (6); доля паводковых вод, поступающих в малые русла (7); пропускная способность малых русел (8); доля социохозяйственной территории (9); уровень развития дорожной сети (10); уровень гидрологической безопасности (11).
В матрице связности ||агз-1| 1,] = 1,11 графа 4 обозначим слабые (медленные) взаимодействия элементами агз = ±1, сильные (быстрые) — элементами агз = ±2. Тогда ненулевые элементы этой матрицы имеют вид: а 18 = а19 = а21 = а28 = а9,ю = +1,
а34 = а41 = а52 = а61 = а71 = а72 = а87 = +2, а9,11 = — 1 а2,11 = а98 = а10,8 = —2.
Анализ построенного графа позволяет выявить тенденции изменения главных факторов динамики паводковой ситуации в ВАП и субъективные оценки их вероятности (за десятилетний период), представленные в таблице 4. Введем следующую простейшую лингвистическую шкалу силы тенденции каждого фактора и ее цифровое обозначение: «сильно (2)», «слабо (1)», «рост (+1;+2)», «падение (-1;-2)», «стабильность (0)». Используя данные таблицы 4, получим 8 трехмерных векторов — сценариев развития комплексной тенденции изменения факторов паводковой динамики. Каждому сценарию с помощью построенного графа межфакторных взаимодействий поставим в соответствие вектор (А У0, А атпх) сценария развития ситуации, непосредственно определяющей динамику оптимальной локальной структуры русел, и функции принадлежности нечеткого оценивания его относительных изменений (рис. 5).
Результаты анализа, представленные в таблице 5, показывают, что существует 4 сценария развития ситуации, субъективные оценки вероятности которых соответственно равны: р(—1; —1) = 0, 729; р(—2; — 1) = 0, 089; р(—1; —2) = 0, 08, р(—2; —2) = 0, 019.
Таблица 4
Тенденции изменения главных факторов паводковой динамики ВАП
№ п/п Фактор Тенденция {оценка вероятности)
1 Средний объем весеннего паводка Сильное (0,1) или слабое (0,9) снижение
2 Доля паводковых вод, поступающих в р. Ахтуба Сильное (0,1) или слабое (0,9) снижение
3 Доля социохозяйственной территории ВАП Сильное (0,9) или слабое (0,1) повышение
Таблица 5
Сценарии развития паводковой ситуации в ВАП
Фактор\комплексная 1 2 3 4 5 6 7 8
тенденция
1 -2 -1 -1 -2 -2 -1 -2
2 -1 -2 -1 -2 -1 -2 -2
3 -1 -1 -2 -1 -2 -2 -2
ДУо -2 -2 -1 -2 -2 -2 -2
Дйтах -1 -1 -2 -1 -2 -2 -2
Субъективная оценка 0,729 0,08 0,08 0,08 0,009 0,009 0,009 0,001
вероятности
о ж
ад
ö ж
з £
ж
— «слабое снижение» «сильное снижение»
1
0.5
0 0.2 0.4 0.6
—AVq
Vo
о ж
ад
ö ж
з £
ж
«слабый рост» «сильный рост»
0.5
0.5
Рис. 5. Функции принадлежности нечеткого оценивания относительных изменений параметров
развития ситуации
1
0
0
0
1
а
О Ж
ад
Ж £
Ж
0.4
0.2
10
An*(T)
о ж
ад
ж £
ж
0.5
1
0
Am* (Т)
Рис. 6. Нечеткие оценки динамики параметров русловой структуры т*, п*
С учетом этих данных, функций, представленных на рисунке 5, и данных таблицы 3 были построены нечеткие оценки динамики параметров русловой структуры т*, п*, представленные на рисунке 6.
0
0
1
2
1
Заключение
Проведенные исследования показали потенциальную возможность управления структурой локальных русловых систем пойменных территорий с целью повышения эффективности использования паводковых вод для их орошения. Главными факторами, определяющими оптимальную структуру этих систем, являются объем паводковых вод и мера неоднородности территории. При малых и средних паводках управление русловой системой может дать существенный эффект. При больших паводках затопление территории определяется главным образом ее рельефом и практически не зависит от структуры ее локальных русловых систем.
ПРИМЕЧАНИЕ
1 Работа выполнена при финансовой поддержке РФФИ и Администрации Волгоградской области в рамках научного проекта № 16-48-340147.
СПИСОК ЛИТЕРАТУРЫ
1. Асарин, А. Е. Речные наводнения: причины и последствия. Что можно и нужно сделать? / А. Е. Асарин // Стратегия гражданской защиты: проблемы и исследования. — 2013. — № 1. — C. 416-427.
2. Землянов, И. В. Анализ экологических последствий эксплуатации Волгоградского водохранилища для сохранения биоразнообразия основных водно-болотных территорий Нижней Волги / И. В. Землянов, О. В. Горелиц, А. Е. Павловский. — М. : Отчет о НИР ФГУ «ГОИН», 2010. — 675 с.
3. Имитационные модели динамики поверхностных вод с использованием данных дистанционного зондирования: влияние рельефа местности / А. А. Воронин, М. В. Елисеева, А. В. Писарев, А. В. Хоперсков, С. С. Храпов // Прикаспийский журнал: управление и высокие технологии. — 2012. — № 3. — C. 54-62.
4. Георгиевский, В. Ю. Гидрологический режим и водные ресурсы / В. Ю. Георгиевский, А. Л. Шалыгин // Методы оценки последствий изменения климата для физических и биологических систем. — М. : Росгидромет, 2012. — C. 53-86.
5. Иванов, В. В. Влияние гидроузлов на деформации пойменных берегов и русловых форм в низовьях реки Волги и Кубани / В. В. Иванов, В. Н. Коротаев // Эрозия почв и русловые процессы. — 2008. — № 16. — C. 224-242.
6. Задача управления гидрологическим режимом в эколого-экономической системе «Волжская ГЭС - Волго-Ахтубинская пойма». Ч. 1. Моделирование динамики поверхностных вод в период весеннего паводка / А. В. Хоперсков, С. С. Храпов, А. В. Писарев, А. А. Воронин, М. В. Елисеева, И. А. Кобелев // Проблемы управления. — 2012. — № 5. — C. 18-25.
7. Особенности динамики затопления Волго-Ахтубинской поймы в зависимости от режимов испарения и инфильтрации / С. С. Храпов, А. В. Писарев, А. А. Воронин, Т. А. Дьяконова, Е. А. Циркова // Вестник Волгоградского государственного университета. Серия 1, Математика. Физика. — 2012. — Т. 16, № 1. — C. 43-47.
8. Численная схема для моделирования динамики поверхностных вод на основе комбинированного SPH-TVD-подхода / С. С. Храпов, А. В. Хоперсков, Н. М. Кузьмин, А. В. Писарев, И. А. Кобелев // Вычислительные методы и программирование. — 2011. — Т. 12, № 1. — C. 282-297.
9. The Numerical Simulation of Shallow Water: Estimation of the Roughness Coefficient on the Flood Stage / S. Khrapov, A. Pisarev, I. Kobelev, A. Zhumaliev, E. Agafonnikova, A. Losev, A. Khoperskov // Advances in Mechanical Engineering. — 2013. — Vol. 5. — Article ID: 787016.
10. Karimipour, F. Watershed delineation from the medial axis of river networks / F. Karimipour, M. Ghandehari, H. Ledoux // Computers and Geosciences. — 2013. — Vol. 53. — P. 132-147.
REFERENCES
1. Asarin A.E. Rechnye navodneniya: prichiny i posledstviya. Chto mozhno i nuzhno sdelat? [River Floods: Causes and Consequences. What Can and Should Be Done?]. Strategiya grazhdanskoy zashchity: problemy i issledovaniya, 2013, no. 1, pp. 416-427.
2. Zemlyanov I.V., Gorelits O.V., Pavlovskiy A.E. Analiz ekologicheskikh posledstviy ekspluatatsii Volgogradskogo vodokhranilishcha dlya sokhraneniya bioraznoobraziya
osnovnykh vodno-bolotnykh territoriy Nizhney Volgi [Analysis of the Environmental Impacts of Operating the Volgograd Reservoir for Biodiversity Major Wetlands of the Lower Volga]. Moscow, Otchet o NIR FGU «GOIN» Publ., 2010. 675 p.
3. Voronin A.A., Eliseeva M.V., Pisarev A.V., Khoperskov A.V., Khrapov S.S. Imitatsionnye modeli dinamiki poverkhnostnykh vod s ispolzovaniem dannykh distantsionnogo zondirovaniya: vliyanie relyefa mestnosti [Simulation Models of Surface Water Dynamics Using Remote Sensing Data]. Prikaspiyskiy zhurnal: upravlenie i vysokie tekhnologii, 2012, no. 3, pp. 54-62.
4. Georgievskiy V.Yu., Shalygin A.L. Gidrologicheskiy rezhim i vodnye resursy [Hydrology and Water Resources]. Metody otsenki posledstviy izmeneniya klimata dlya fizicheskikh i biologicheskikh sistem. Moscow, Rosgidromet Publ., 2012, pp. 53-86.
5. Ivanov V.V., Korotaev V.N. Vliyanie gidrouzlov na deformatsii poymennykh beregov i ruslovykh form v nizovyakh reki Volgi i Kubani [Effect of Hydro Strain on the Shores of the Floodplain and Channel Forms in the Lower Reaches of the Volga and Kuban]. Eroziya pochv i ruslovye protsessy, 2008, no. 16, pp. 224-242.
6. Khoperskov A.V., Khrapov S.S., Pisarev A.V., Voronin A.A., Eliseeva M.V., Kobelev I.A. Zadacha upravleniya gidrologicheskim rezhimom v ekologo-ekonomicheskoy sisteme «Volzhskaya GES - Volgo-Akhtubinskaya poyma». Ch. 1. Modelirovanie dinamiki poverkhnostnykh vod v period vesennego pavodka [The Regimen Control Task in the Eco-Economic System «Volzhskaya Hydroelectric Power Station — the Volga-Akhtuba Floodplain». I. Simulation of Dynamics of Surface Water During Spring Floods]. Problemy upravleniya [Control sciences], 2012, no. 5, pp. 18-25.
7. Khrapov S.S., Pisarev A.V., Voronin A.A., Dyakonova T.A., Tsirkova E.A. Osobennosti dinamiki zatopleniya Volgo-Akhtubinskoy poymy v zavisimosti ot rezhimov ispareniya i infiltratsii [The Role of Infiltration and Evaporation in the Flooding Dynamics of the Volga-Akhtuba Floodplain]. Vestnik Volgogradskogo gosudarstvennogo universiteta, Seriya 1, Matematika. Fizika [Science Journal of Volgograd State University. Mathematics. Physics], 2012, vol. 16, no. 1, pp. 43-47.
8. Khrapov S.S., Khoperskov A.V., Kuzmin N.M., Pisarev A.V., Kobelev I.A. Chislennaya skhema dlya modelirovaniya dinamiki poverkhnostnykh vod na osnove kombinirovannogo SPH-TVD-podkhoda [A Numerical Scheme for Modeling the Dynamics of Surface Water Based on the Combined SPH-TVD-Approach]. Vychislitelnye metody i programmirovanie, 2011, vol. 12, no. 1, pp. 282-297.
9. Khrapov S., Pisarev A., Kobelev I., Zhumaliev A., Agafonnikova E., Losev A., Khoperskov A. The Numerical Simulation of Shallow Water: Estimation of the Roughness Coefficient on the Flood Stage. Advances in Mechanical Engineering, 2013, vol. 5, article ID: 787016.
10. Karimipour F., Ghandehari M., Ledoux H. Watershed Delineation From the Medial Axis of River Networks. Computers and Geosciences, 2013, vol. 53, pp. 132-147.
OPTIMIZATION AND SCENARIO-SIMULATION MODELLING OF THE DYNAMICS OF THE STRUCTURE OF SMALL ARTIFICIAL CHANNELS
IN FLOODPLAIN AREAS
Vladislav Sergeevich Bakulin
Student,
Department of Fundamental Informatics and Optimal Control, Volgograd State University gitara_93@mail.ru, fiou@volsu.ru
Prosp. Universitetsky, 100, 400062 Volgograd, Russian Federation
Anna Anatolyevna Vasilchenko
Senior Lecturer, Department of Fundamental Informatics and Optimal Control, Volgograd State University aa-vasilchenko@mail.ru, fiou@volsu.ru
Prosp. Universitetsky, 100, 400062 Volgograd, Russian Federation
Alexander Aleksandrovich Voronin
Doctor of Physical and Mathematical Sciences, Professor, Head of Department of Fundamental Informatics and Optimal Control, Volgograd State University voronin.prof@gmail.com, fiou@volsu.ru
Prosp. Universitetsky, 100, 400062 Volgograd, Russian Federation
Mikhail Alekseevich Kharitonov
Junior Researcher, Department of Fundamental Informatics and Optimal Control,
Volgograd State University
kharitonov.mihail@gmail.com, fiou@volsu.ru
Prosp. Universitetsky, 100, 400062 Volgograd, Russian Federation
Abstract. The purpose of this paper is to study the properties of optimality of local artificial channel systems of floodplain territories. A method for searching the topological structure and cross-sectional area of artificial channels maximizing the total weighted area of the flooded territory is proposed. For the adequacy of the formulation and solution of the control problem, an adaptive zoning of the digital model of the floodplain territory to areas with an insignificant effect of interband flooding at small and medium floods has been carried out. The optimal throughput of small channels was sought with the help of a multi-step recursive dynamic programming method. The results of numerical realization of the algorithm using computer hydrodynamic imitations of flood dynamics and geoinformation modeling of spatially distributed relief parameters of several sections of the territory of the northern part of the Volga-Akhtuba floodplain with artificial local channel systems are presented.
A weighted oriented graph of interfactor interactions determining the flood dynamics of the northern part of the Volga-Akhtuba floodplain was constructed, with the help of which the scenarios for the development of the situation and the dynamics of parameters directly determining the optimal local channel structure were formalized. Subjective estimates of the probability of scenarios and fuzzy estimates of the dynamics of the channel structure parameters are given.
Key words: structural optimization, dynamic programming, production function, geoinformation modeling, hydrodynamic modeling, flow systems.
