CC BY
65
УДК 519.6 + 519.8 ББК 22.19+20.1
ПРОЕКТИРОВАНИЕ СИСТЕМЫ ЭКОЛОГО-ЭКОНОМИЧЕСКОГО УПРАВЛЕНИЯ ТЕРРИТОРИЕЙ ВОЛГО-АХТУБИНСКОЙ ПОЙМЫ НА ОСНОВЕ
ГИДРОДИНАМИЧЕСКОГО И ГЕОИНФОРМАЦИОННОГО МОДЕЛИРОВАНИЯ
Воронин А. А.1, Васильченко А. А.2, Писарева М. В.3, Писарев А. В.4, Хоперсков А. В.5, Храпов С. С.6, Подщипкова Ю. Е.7
(Волгоградский государственный университет, Волгоград)
Описана система эколого-экономического управления территорией Волго-Ахтубинской поймы, построенная на основе гидродинамического и геоинформационного моделирования. Приведены результаты оптимизации паводкового гидрографа Волжской ГЭС, пространственного распределения работ по проекту восстановления русел Волго-Ахтубинской поймы, механизмов управления хозяйствующими субъектами территории.
Ключевые слова: компьютерное моделирование, гидрологический режим, гидрограф, цифровая модель рельефа, Вол-го-Ахтубинская пойма, эколого-экономическое управление.
1 Александр Александрович Воронин, доктор физико-математических наук, профессор (voronin.prof@gmail.com).
2 Анна Анатольевна Васильченко, старший преподаватель (aa-vasilchenko@mail.ru).
3 Мария Владимировна Писарева, аспирант (^. eliseeva. vlg@gmail.com).
4 Андрей Владимирович Писарев, кандидат физико-математических наук, доцент (andrew.pisarev@gmail.com).
5 Александр Валентинович Хоперсков, доктор физико-математических наук, профессор (ahoperskov@gmail.com).
6 Сергей Сергеевич Храпов, кандидат физико-математических наук, доцент (xss-ip@mail.ru).
7 Юлия Евгеньевна Подщипкова, студент (podschipkova1993@mail.ru).
Управление большими системами. Специальный выпуск 55 1. Введение
Волго-Ахтубинская пойма (ВАП) - уникальное природное образование (протяженность около 450 км, территория - свыше 20 тыс. кв. км), жизнь которого полностью определяется весенним паводком р. Волги. Волжская ГЭС (ВГЭС) регулирует подчиненный интересам гидроэнергетики гидрологический режим р. Волги, характеризуемый увеличенным объемом ее меженного стока и значительно сокращенным объемом паводкового, что служит причиной прогрессирующей деградации русел и обезвоживания ВАП [5]. Факторами обезвоживания являются также дополнительное ограничение переменных паводков требованиями гидрологической безопасности выходящих из зоны затопления осваиваемых территорий, а также расширение и углубление русла р. Волги. Результатом предпринимаемых в последние годы усилий по улучшению гидрологического режима северной части ВАП (далее - ВАП), находящейся на грани экологической катастрофы, является лишь замедление деградационных процессов [6]. Для коренного изменения ситуации необходима разработка и реализация научно обоснованной комплексной системы эколого-экономического управления, создание которой на территории ВАП - пространственно распределенной высокодинамичной многоагентной слабоустойчивой социоприродохозяй-ственной системы с высоким уровнем техногенной опасности, социальной и экологической значимости - возможно только на основе междисциплинарного подхода, интегрирующего методы и результаты геоинформационного, гидродинамического, оптимизационного, теоретико-игрового и сценарно-имитационного моделирования.
Главной целью проектируемой комплексной системы эко-лого-экономического управления ВАП является создание математической модели ее устойчивого развития на основе достижения и поддержки устойчивости ее экосистемы в условиях оптимального гидрологического режима. Она включает в себя несколько постоянно актуализируемых модулей - математических моделей и отвечающих им программных комплексов, реализующих основной функционал: когнитивный анализ проблемной ситуации; многокритериальную оптимизацию паводкового гид-
рографа ВГЭС и природовосстановительных проектов; проектирование механизмов эколого-экономического управления; сце-нарно-имитационное моделирование социоприрохозяйственной динамики и развития.
Функционирование главной части комплекса обеспечивается вспомогательным функционалом, реализуемым программным комплексом «ЭКОГИС», включающим в себя построение и актуализацию цифровой модели рельефа (ЦМР) на основе ГИС-технологий, численное гидродинамическое моделирование динамики поверхностных вод. Первые версии основных моделей и результатов моделирования описаны в [4, 5, 7-9, 11].
2. Геоинформационное моделирование рельефа
Разработка цифровой модели рельефа (ЦМР) ВАП b(x, y) была начата в 2007 г. Ее основой служат спутниковые данные ASTER [14] и SRTM [12] с разрешением до 20 м в плоскости Земли и до 0,5 м по вертикали. Программа Landsat [13] позволяет получить актуальные данные высокого разрешения (15-30 метров/пиксель). Для добавления объектам атрибутивной информации использовалась топографическая карта 1:50000 открытого пользования. Цифровая карта пойменных русел и речного дна строилась совместной векторизацией топографических и лоцманских карт русел ВАП, Волги и Ахтубы [1] с использованием результатов экспедиционных исследований [2]. Актуализация ЦМР производится ежегодно с использованием регулярно обновляемых космических снимков, размещенных на сайте геологической службы США [15] в открытом режиме (спутник ДЗЗ Landsat 8), а также с помощью GPS-измерений границ паводковых затоплений. Векторная карта ВАП 2014 года включает в себя слой гидросистемы, состоящий из 1542 русловых объектов общей протяженностью 887 км, слой инфраструктуры, включающий в себя 118 населенных пунктов, и слой рельефа местности, включающий более 15 тыс. рельефных объектов.
Рис. 1. Карта русловых зон северной части ВАП
Естественной основой структуризации пойменной территории, необходимой для проектирования системы управления, служит ее гидрологическая сеть. Магистральную гидросистему ВАП образуют три водотока: р. Волга, р. Ахтуба и ерик «Гнилой» («Пахотный»), из которых берут начало соответственно 14, 15 и 57 русел - входов на ее территорию. Эти русла, практически равномерно (с интервалом в 1,5-2 км) распределенные вдоль магистральных водотоков, разветвляются на территории поймы, образуя (в версии 2014 г.) 86 гидросистем второго уровня и соответствующих им территориальных зон (рис. 1) с функциями рельефа ь(%^, у}) = 1, ..., 86).
3. Гидродинамическое моделирование паводков
Численная модель динамики поверхностных вод в версии 2014 г. учитывает все основные факторы затопления территории: поверхностные и подземные источники воды - плотины, осадки, ключи, выход грунтовых вод на поверхность суши; рельеф местности с учетом антропогенной застройки территорий и рельефа дна водоемов; свойства подстилающей поверхности -придонное трение, инфильтрация (новая многослойная нелинейная модель); внутреннее вязкое трение; ветровое воздействие - нагонные волны; вращение Земли - сила Кориолиса; испарение. Адекватность модели подтверждается результатами сравнения с данными наблюдений уровня воды на четырех гидропостах в северной части ВАП и площади водной поверхности в различные годы с использованием данных спутника Landsat-7.
Реализованы параллельные OpenMP, CUDA-, OpenMP-CUDA-версии расчетного модуля метода CSPH-TVD, снизившие в 700 раз вычислительную сложность задачи и позволившие довести время расчета 100-суточного паводка до 2 часов.
Вид функций расхода (гидрографа) Q(t) паводкового попуска ВГЭС разных лет (рис. 2) демонстрирует его относительную хаотичность в 1962-1996 гг., обусловленную следованием требованиям безопасности ВГЭС и переменному спросу на электроэнергию. В 2003-2015 гг. в нем выделяются две ступени с почти постоянными значениями Q1 и Q2, а меженный гидрограф характеризуется меньшей суточной изменчивостью. Анализ космических снимков и результатов гидродинамических расчетов показал, что перед весенним попуском воды относительная площадь занятой водой территории не превышает 1%, площадь затопления территории двухступенчатым гидрографом определяется величиной Q1 и длительностью ft его первой ступени, при Q1 < 20 000 м3/с затопления территории не происходит, а при Q1 > 30 000 м3/с и t1 > 5 суток соответствуют катастрофическому затоплению.
На рис. 3 показана расчетная динамика площади территории паводкового затопления iSb(t) с плановыми гидрографами 2006-2011 гг.
На рис. 4, 5 показана расчетная динамика затопления всей территории ВАП и указанных на рис. 1 ее зон модельными (Q1 = const) паводковыми гидрографами. Различия в динамике зон связаны с характеристиками их рельефа и гидросистемы.
Используемая в настоящее время гидродинамическая модель не учитывает прогрессирующей деградации русел ВАП и рек Волги и Ахтубы, снижающей эффективность паводков. По данным гидропоста в г. Волгограде снижение высоты паводка вследствие изменений в русле р. Волги составляет 1 м (рис. 6). Деградация русла р. Ахтубы и русел территории ВАП добавляет к этой величине значительную долю, для учета которой необходимы широкомасштабные полевые исследования.
01.04 14.04 28.04 05.05 26.05 09.06 23.06 I Рис. 2. Паводковые гидрографы ВГЭСразных лет
6x1 о8
30 40 50 60 70 80 /
Рис. 3. Динамика площади затопления Бо (м2) для гидрографов за 2006-2011 гг. ([¿] = сутки)
Рис. 4. Динамика относительной площади затопления Б (м2) для постоянных модельных гидрографов ([¿] = сутки)
4. Моделирование системы управления
4.1. ОБЩАЯ МОДЕЛЬ УПРАВЛЕНИЯ
Постановка задачи управления, с одной стороны, является результатом когнитивного анализа ситуации и экспертного поиска баланса между экологическим, социальным и экономическим краткосрочным и долгосрочным приоритетами различных групп акторов в условиях требований безопасности, деградаци-онной природной динамики, изменчивой гидрологии волжского бассейна, гидротехнических ограничений, с другой - определяется возможностью построения соответствующей математической модели.
10 8 6 4
2 0 о
Рис. 5. Динамика затопления зон ВАП для (\ = 23000 м3/с
H, m BS 1 l 1961 . t 2010
- l l l * / í""""^ 2011 2012
- r .¡Mr
i ♦
- 1 1
- / « 1 l 1 00
1 l 1
- 1 1 l
- i 1 1 i Q, m3/s
1 5000 10 1 . . 1 , 1 1 . I 00 15000 20 , ... I..... 00 25000 30
Рис. 6. Зависимость уровня воды от паводкового гидрографа ВГЭС в разные годы (по гидропосту г. Волгограда)
В проектируемой модели управления состояние системы описывается кортежем S(x, y, l, t) = [A(l, t), B(x, y, l, t), C(A, B, Ext), U(A, B, Ext), F(A, B, U), (Ext(l, t))] 6 Sдоп, где l = 1, 2, ... - «медленное» время (шаг по l равен одному году); A(l, t) = (Ao(l, t), Aij(l)) (i = 1, .., n; j = 1, .., k), Ао = {fo, wo} - модель агента ВГЭС: fo - целевая функция, wo(l) = (Оф, ro(l)) - действие, Оф = 2ф(1, t) - фактический паводковый гидрограф, t -«быстрое» время паводка, ro - объем софинансирования проектов Центра; A,j = fj, Wj(l)} - модели агентов (хозяйствующих субъектов) зон ВАП: f - целевые функции, w,j - действия, Wij(l) = (dij(l), rij(l)), dij - изменение рельефа (число созданных агентами дамб, удерживающих после окончания паводка необходимую для орошения воду в руслах зон), rj - объем софинансирования проектов Центра, n - число зон (1 < n < 86), ki - число агентов зоны i; B(x,y, l) = (bi, ..., be) - пространственная сеточная вектор-функция с координатами: bi(x, y, l) - ЦМР; b2 = ф(Оф, bi(x, y, l)) - максимальная (по t) высота паводкового затопления, определяемая гидродинамическими расчетами, b3(x, y, l) = sgn(b2(x, y, l)) - индикатор затопления территории;
Ьл(х, y, l), bs(x, y, l), Ьб(х, y, l) - векторные коэффициенты соответственно экологической (гидрологической, ихтиологической, орнитологической и др.), социальной (урбанистической, рекреационной и др.) и экономической (сельское, лесное хозяйство, рекреационный бизнес и пр.) значимости, рассчитываемые на основе данных мониторинга и экспертных оценок (число координат каждого вектора определяется наличием данных и целями моделирования); C(A, B, Ext) = {U(A, B, Ext), F(A, B, U)} - модель Центра (органа федеральной власти), F - целевая функция Центра, U = (Gm, W, P, M) - действие Центра (управление): Опл(Т) = [Qim, t\m, Qim, t2rn), - плановый двухступенчатый паводковый гидрограф; W = (Wi,..., Wn) - проекты русловосстановле-ния в зонах ВАП, P - штрафы за сверхнормативные действия, M - стимулирование участия в проектах; Бдоп - множество допустимых состояний системы S, определяемое природными, техническими, социально-экономическими и финансовыми ограничениями.
В рамках построенной модели S различные задачи эколого-экономического управления описываются следующими соотношениями:
(1) Bt(l + 1) = f(B„ Wi, Wi, Ext),
(2) F ^ maxu,F = Z(i,k,0e aC%(l)Blk(l) " (r, U(l))),
(3) fo = ао(Уш - Уф ) - P ( Gm - ) ^ max{ ,
(4) f = aijSi3 - P(wv ) - rt] (M) ^ maxw¡j,
(5) S(xy,T) g Ss (g2,g3,g4),
(6) S( X, y, l ) G Ss ( g2, g3) при всех I
(7) V(l) < Vo(l), V(l) = ZaM (l)By при всех I,
(i = 1, ..., n; j = 1, ..., k) Q = {l = 1, ..., T; i = 1, ..., n; k = 1, ..., 6}, 1 < n < 86, где b, = (Bii, ..., Bi6), f = (fi, ...,fu), Уф = Ql t)dt, Si3(l) = maxtBi3(l), Bik(l) = \bk(x, y, l)dx,dy, - интегральные функции зоны i; ao, aj - коэффициенты экономической эффективности агентов, Wi = {Wij} - вектор действий агентов зоны i, Ext - вектор-функция сценариев состояния внешней среды (уровень наполненности Волгоградского водохранилища, социохозяй-ственные воздействия и ограничения), Упл = Qimti пл + Q2Knt2Kn;
V - экологический ущерб (Уо - его нормативное значение); & -состояние стабильности системы (в котором действия Центра компенсируют природную и антропогенную деградацию, и средние значения показателей паводкового затопления и пойменных популяций определяются средним значением паводкового гидрографа, т.е. существуют однозначные функции gk (к = 2, 3, 4), определяемые соотношениями
Ьк(х,у, I) = gk(Gпл(l), х,у)); г = (г\, ..., г4), г\ - стоимость упущенной гидроэнергетической полезности от реализации Gпл, Г2 - цена работ по проектам; г2, г2 - затраты на мониторинг, необходимый для синтеза механизмов, Т - горизонт планирования, Цк (к = \, ..., 6) - весовые коэффициенты.
В зависимости от степени адекватности каждого из соотношений (\)-(7) и целей управления Центр может решать следующие задачи:
- (\)-(4) - задача оптимального управления и ее частные случаи (см. далее);
- (\)-(5) - оптимальная стабилизация экосистемы;
- (\), (3), (4), (5) - поддержание стабильности экосистемы;
- (\), (3), (4), (7) - управление гидрологическим риском;
- (\)-(6) ((\)-(7)) - устойчивое развитие (с экологическим нормативом).
На каждом шаге I динамической оптимизации поиск решения производится на множестве решений соответствующих иерархических игр и условной параметрической оптимизации, что требует численного построения функций Б,к(1, А, V) в широком диапазоне параметров. Функции gk определяются в ходе решения соответствующих задач. Ниже приведены результаты исследования некоторых частных случаев задачи (1)-(4) - параметрической оптимизации паводкового гидрографа и проекта расчистки русел с пассивными (уравнения (2)) и активными (система (2),(3)) агентами.
4.2. МОДЕЛЬ ОПТИМИЗАЦИИ ПАВОДКОВОГО ГИДРОГРАФА ВГЭС
Оптимизация паводкового гидрографа ВГЭС может исследоваться и реализовываться как в составе комплексной задачи управления, так и в виде самостоятельной задачи. Она состоит в
поиске функции Gm как решения задачи оптимизации с несколькими критериями: площадь нерестилищ (Spx > Sonm), длительность их затопления (Трх > Топт), максимальная площадь затопления
(S3 = ТЦ=1 Sl3 ) ^ Sonm , объем гидрографа (Vm ^ min) в условиях гидрологической безопасности (Qi < Qm).
Когнитивный анализ ситуации показывает, что фактическая иерархия критериев S3^V<Spx<Tpx
сводит многокритериальную задачу к параметрической одно-критериальной с параметром Si и заданными величинами Qm,
ТРХ'
(8) Упшн ^ minGm , Spx = S1,
решение которой существует при Si < Smax, Smax = lim(t) S3(Qm, t). Ее решением на множестве одно- и двухступенчатых гидрографов является двухступенчатый гидрограф G* = {Q*i(=Qm), t*i(Si), Q*2(Si), t*2(Si)} [5]. График функции Smin(Si) = Sз(Sl, Gorm) проведен жирной линией на рис. 1.
Эколого-экономическая иерархия критериев V<S3<Spx<Tpx
переводит (8) в задачу:
(9) VmaH ^ mnGm , SPx = SU S = ^
решение которой существует при выполнении неравенств Si < S2 < Smax. На рис. 1 штриховкой и цифрами выделены области решения задачи (9) в виде одноступенчатого (1) и двухступенчатого (2) гидрографов. (Вопрос о свойствах гидрографа с большим числом ступеней остается открытым.) Сравнение рис. 2 и рис. 7 показывает, что фактический переход от одноступенчатого к двухступенчатому гидрографу можно объяснить целевым уменьшением отношения S2 / Si.
4.3. МОДЕЛЬ ОПТИМИЗАЦИИ ПРОЕКТА ВОССТАНОВЛЕНИЯ РУСЕЛ
Задача параметрической оптимизации проекта восстановления русел имеет вид
^(й А,яц М,Ь(Х,У )) ^ тахб(лу,),
= щ(вх,ь,яц,ь2,ь(хг,уг)),
где ASi - проектное приращение площади паводкового затопления зоны г; Яц - затраты Центра на реализацию проекта,
£ = Иги=1 , где Ь - общая длина восстановленных русел в зоне г;
Лг - проектное приращение глубины г русел.
Рис. 7. График Smin(Spx) и области решений задачи (9), отвечающие одноступенчатому (1) и двухступенчатому (2) гидрографам
В созданной ЦМР ширина всех русел ВАП равна 20 м, глубина - z = 2 м. Расчеты [3] показали, что для i = 1, ..., 20 при различных Qi, Яц, t\ решением задачи Д5г ^ max^ , (Az = 0; 0,5; 1; ...), Az L = const, является значение Az* = 1 м, которое использовалось в дальнейшем. (Далее считаем, что Яц = rrL). Предположение о независимом затоплении гидрологических зон ВАП позволяет свести оптимизацию по простран-
ственно распределенной переменной Ъ(х, у) к двухуровневой векторной оптимизации с параметрами Ql, ¿1, Ь: (10) ЛСД, К Ь) ^ тахь, Ь = (Ц, .., Ь„),
АБ = АБ1(вх,1х,г2,Ь1).
Функции ЛХ в (10) являются решением задач оптимизации работ внутри каждой зоны с параметром Ьг (общей длиной восстановленных русел в зоне г) и числом альтернатив ](Ы) в соответствии с числом русловых разветвлений. В связи с большой трудоемкостью построения виртуальных ЦМР с различными глубинами русел (Ьвап = 887 км) ниже приведены результаты решения частной задачи (10): п = 20, Ьг = 1 км и Ьг = 2 км (г =1,...,20). Зависимости ЛХ^ь 11, Ь) (г = 1, ..., 20) строились численно для серии параметров Ql, 11 и различных цифровых рельефов ВАП, отличающихся высотами дна русел выбранных зон. На рис. 9-12 представлены результаты расчета динамики ЛБг(Ьг) и ЛХ(0 (г = 1, ..., 20) для нескольких значений параметров Ql и ¿1. Представленные рисунки демонстрируют общую тенденцию начального роста с последующей стабилизацией или падением значений функций ЛХ(0 с различиями в амплитудах и временной локализации. При Ql = 17 тыс. м3/с в отдельных зонах наблюдаются отрицательные значения ЛХ, растущие по абсолютной величине с ростом Ь, что показывает доминирование эффекта накопления вод в углубленных руслах над их разливом по территории. С ростом Ql это доминирование исчезает. Нерегулярный характер зависимости ЛХ(^) для некоторых зон связан с нерегулярностью рельефа. Путем упорядочения зон ВАП по ЛХ были построены приведенные на рис. 13-16 решения частной задачи (10) АБ(Ь) для Ьг = 1 км; 2км (г = 1, ..., 20) и различных Ql, ¿1.
Рис. 8. Зависимость А(г = 1, ..., 20) при Q\ = 20 000 мъ/с,
¿1 = 6 дней
Рис. 9. Зависимость А%(0 (г = 1, ..., 20) при Ql = 17 000 мъ/с,
Ьг = 1 км
Л5г,кв. км
8 - —9 -^—10—^—11 —•—12—13-14
15 16—•—17-18 19 20
Рис. 10. Зависимость Д£г(0 при Ь, = 2 км (/ = 1, ..., 20) Ql = 17 000 м3/с
8 -9 —10—11-*-12-в-13 14
15—16——17-18 19 20
Рис. 11. Зависимость Д$() при Ь, = 1 км, Ql = 23 000 м3/с (/ = 1, ..., 20)
-*-15-*<-16-^17—18 19 20
Рис. 12. Зависимость АБ() при Ь = 2 км, Ql = 23 000 м3/с
(г = 1, .. 20)
Рис. 13. Зависимость А^(Ь) при Ь = 1 км;2 км, Ql = 17 000 м3/с, ¿1 = 3 дня (г = 1, ..., 20)
ЛЯ КЕ_ КМ
1 0 2 Э 3 0
1
—♦—1км -»-2 км
Рис. 14. Зависимость Д£(Ь) при Ь, = 1 км;2 км, Ql = 17 000 м3/с, ¿1 = 6 дней (/ = 1, 20)
КБ. КМ
О 4 8 12 16 20 24 28 32 3« 40 —♦—1м *2н
Рис. 15. Зависимость Д£(Ь) при Ь, = 1 км; 2 км, Ql = 23 000 м3/с, Н = 3 дня (/ = 1, ..., 20)
Л 5; ке. км
О 4 3 12 16 20 24 28 32 36 40 —♦—1км -в-т^,
Рис. 16. Зависимость Д£(Ь) при Ь, = 1 км;2 км, Ql = 23 000 м3/с, ¿1 = 6 дней (/ = 1, 20)
4.4. МОДЕЛИ СИНТЕЗА ОПТИМАЛЬНЫХ МЕХАНИЗМОВ УПРАВЛЕНИЯ
Целями механизмов управления агентами ВАП в проектируемой системе управления являются поддержание планового состояния системы (&, В™, S3(Sрх, Gпл), Gпл), плановой траектории ее управляемой динамики, снижение затрат г на действия Центра. В [5] представлен механизм управления агентом ВГЭС Po, поддерживающий равенство Sз(SрX, Gфакт, Po) = Sз(SрX, Gm).
Опишем процедуру синтеза механизма софинансирования Агентами ВАП проекта русловосстановления, включение которого в задачу (10) увеличивает общее финансирование R = Rц + RА = г^Ъ, эквивалентное уменьшению для Центра стоимости работ в каждой из зон, потенциально изменяет ее решение. (Горизонт планирования предполагается единым для всех участников проекта.)
Стратегией Центра является выбор величин ду - долей агентов Лу в финансировании работ по восстановлению русла каждой из зон длиной Ьг. Агент Лу, информированный о величине Цу и функции АSг(Gm, Ьг), находит оптимальную для него величину
т *
Ь у восстановленного русла как решение задачи
= а щ (Опя^,) - д13 г ^ ^ тахЬу. Оптимальная стратегия агента Лгу:
Ц (Ц )=(№,Г%г2/а„),
где Ц (дг ^(Л^ /А) - оптимальная стратегия коллектив-
ного агента зоны г - решение задачи
Рг = Е 11г] = А ( Ц ) " 1, Г2 Ц ^ таХЦ ,
А =Еа,
Цг > 0 - определяемая Центром доля участия коллективного Агента зоны г в финансировании работ по расчистке русла длиной Ьг. Согласование интересов агентов Лу (Ь у = Ь г, у = 1, ..., к) достигается при выполнении равенств
1* = ау1*/Аг (г = 1, п,у = 1, ..., кг),
где Цг*(Ьг) - оптимальная стратегия Центра, определяемая из уравнения Ь (д,) = Ьг.
Обозначим
q* = argmax^ (qtr2L] (qt)),
тогда L* =(ASi)- (q*r2 /Д.) - максимальное значение желательной для Агентов зоны i величины Li. Таким образом, оптимальные для Центра цены работ в зонах имеют вид
* fr2(1 - q*X L, ^ L*,
, ^ i > L*;
(i = 1,..., n) .
Далее Центр решает задачу (10) с условием Кц (L) = L.
5. Анализ результатов
Как показали расчеты, действия Центра образуют единую систему и эффективность каждого из них обусловлена всей системой действий, что доказывает актуальность комплексной задачи эколого-экономического управления. Проектирование адекватной системы управления ограничивается алгоритмической сложностью и неопределенностью задачи. Вычислительная сложность гидродинамических расчетов и оптимизационных задач в значительной мере преодолевается распараллеливанием алгоритмов. Объективную неопределенность комплексной модели управления придает ее зависимость от внешних сценариев Ext. Субъективная неопределенность обусловлена недостаточной разработанностью теоретических основ динамики и рисков в слабоустойчивых распределенных социоприродохозяйствен-ных системах, отсутствием данных пространственного мониторинга и социально-экономической статистики в ВАП [6, 10].
В настоящее время обоснована адекватность нескольких (описанных здесь) частных моделей - составных частей проектируемой комплексной системы: моделей рельефа и паводкового затопления, оптимизации гидрографа и восстановления русел. Однако обоснование практической значимости их результатов требует обширных исследований по определению и повышению их точности. Сравнение с данными спутниковых измерений показывает, что относительная погрешность вычисления
r2i =
величины S3 для всей территории ВАП равна 4%-6%, тогда как для небольших по размеру равнинных зон и больших гидрографов она может достигать 30%.
6. Заключение
В силу ограниченного объема настоящей статьи в ней не описаны модели когнитивного анализа (SWOT- , PEST-анализ, когнитивные карты) и риска, сценарно-имитационная модель социоприродохозяйственной динамики ВАП. Недостаточная практическая адекватность уравнений (1) является сдерживающим фактором проектирования комплексной модели управления социально-экономической системой, концентрирующей в себе большое число видов методологической и алгоритмической сложности и неопределенности, что актуализирует новые теоретические и практические задачи моделирования динамики и механизмов управления, параллельных вычислений, мониторинга и социально-экономической статистики территорий.
Работа выполнена при частичной финансовой поддержке гранта РГНФ 14-12-34019 (комплексная модель эколого-экономического управления, оптимизация гидрографа, проекта русловосстановления, механизмов управления), гранта РФФИ 13-01-97062 (модель динамики поверхностных вод), гранта РФФИ 14-17-97030 (использование ГИС для актуализации цифровой модели рельефа), гранта РФФИ 13-07-97056 (программный комплекс моделирования динамики поверхностных вод). Авторы благодарят суперкомпьютерный комплекс МГУ «Ломоносов» за предоставленные вычислительные ресурсы.
Литература
1. Атлас единой глубоководной системы Европейской части РФ. Т. 7: Река Волга. От Волгоградского гидроузла до Астрахани. - СПб.: Б.и., 2009. - 34 л.
2. БАРМИН АН., ГОЛУБ В.Б., ИОЛИН М.М. и др. Индикация изменений условий среды в северной части ВАП при использовании шкал Л.Г. Раменского и DCA-
ординации // Геодезия и аэросъемка. - 2010. - №5. -С.21-24.
3. ВОРОНИН А.А., ВАСИЛЬЧЕНКО А.А., ХРАПОВ С.С.
и др. Анализ эффективности природовосстановитель-ных проектов в эколого-экономической системе «Волжская ГЭС - Волго-Ахтубинская пойма» // Управление большими системами. -2014. - №52. - С. 133-147.
4. ВОРОНИН А.А., ЕЛИСЕЕВА М.В., ПИСАРЕВ А.В. и др. Имитационные модели динамики поверхностных вод с использованием данных дистанционного зондирования: влияние рельефа местности // Прикаспийский журнал: управление и высокие технологии. - 2012. - №3(19). -С.54-62.
5. ВОРОНИН А.А, ЕЛИСЕЕВА М.В., ХРАПОВ С.С. и др.
Задача управления гидрологическим режимом в эколого-экономической системе «Волжская ГЭС - Волго-Ахтубинская пойма». Ч2. Синтез системы управления // Проблемы управления. - 2012. - №6. - С. 19-25.
6. ЗЕМЛЯНОВ И.В., ГОРЕЛИЦ О.В., ПАВЛОВСКИЙ А.Е. и др. Анализ экологических последствий эксплуатации Волгоградского водохранилища для сохранения биоразнообразия основных водно-болотных территорий Нижней Волги // Отчет о НИР ФГУ «ГОИН», 2010. - 675 с.
7. ХОПЕРСКОВ А.В., ХРАПОВ С.С., ПИСАРЕВ А.В. и др. Задача управления гидрологическим режимом в эколого-экономической системе «Волжская ГЭС - Волго-Ахтубинская пойма». Ч. 1. Моделирование динамики поверхностных вод в период весеннего паводка // Проблемы управления. - 2012. - №5. - С. 18-25.
8. ХОПЕРСКОВ А.В., ХРАПОВ С.С., ПИСАРЕВ А.В. Прямое моделирование динамики поверхностных вод на территории Волго-Ахтубинской поймы // Суперкомпьютерные технологии в науке, образовании и промышленности. - 2012. - С. 177-181.
9. ХРАПОВ С.С., ПИСАРЕВ А.В., ВОРОНИН А.А. Особенности динамики затопления Волго-Ахтубинской поймы в зависимости от режимов испарения и инфиль-
трации // Вестник ВолГУ. Серия 1: Математика. Физика. - 2012. - Т. 15, №1.5. - С. 36-42.
10. ШЕВАНДРИН А.В., ПЕТРОВА Е.А., ВОРОНИН А.А.
Методика мониторинга и оценки уровня социально-экономического развития социохозяйственных слабоустойчивых природных систем (на примере территории Волго-Ахтубинской поймы) // Фундаментальные исследования. - 2014.- № 11-10.- С. 2233-2237.
11. KHRAPOV S., PISAREV A., KOBELEV I. ETC. The Numerical Simulation of Shallow Water: Estimation of the Roughness Coefficient on the Flood Stage // Advances in Mechanical Engineering. - 2013. - Vol. 2013. - Article ID 787016. - 11 p.
12. RABUS B., EINEDER M., ROTH A., ET AL. The shuttle radar topography mission - a new class of digital elevation models acquired by spaceborne radar // ISPRS Journal of Photogrammetry and Remote Sensing. - 2003. -Vol. 57(4). - P. 241-262.
13. ROY DP., WULDER M.A., LOVELAND T.R. ETC. Land-sat-8: Science and product vision for terrestrial global change research // Remote Sens. Environ. - 2014. -No. 145. - P. 154-172.
14. SUWANDANA E., KAWAMURA K., SAKUNO Y. ETC. Thematic information content assessment of the ASTER GDEM: A case study of watershed delineation in West Java, Indonesia // Remote Sensing Letters. - 2012. - Vol. 3(5). -P. 423-432.
15. [Электронный ресурс] - URL: http://earthexplorer.usgs.gov (дата обращения: 20.02.2012).
DESIGNING A SYSTEM FOR ECOLOGICAL-ECONOMICAL MANAGEMENT OF THE VOLGA-AKHTUBA FLOODPLAIN ON BASIS OF HYDRODYNAMIC AND GEOINFORMATIONAL SIMULATION
Alexander Voronin, Volgograd State University, Volgograd, Doctor of Science, professor (voronin.prof@gmail.com). Anna Vasilchenko, Volgograd State University, Volgograd, Senior Lecture (aa-vasilchenko@mail.ru).
Maria Pisareva, Volgograd State University, Volgograd, postgraduate (m.eliseeva.vlg@gmail.com).
Andrei Pisarev, Volgograd State University, Volgograd, cand. of Science, associate professor (andrew.pisarev@gmail.com). Alexander Khoperskov, Volgograd State University, Volgograd, Doctor of Science, professor (ahoperskov@gmail.com). Sergey Khrapov, Volgograd State University, Volgograd, cand. of Science, associate professor (xss-ip@mail.ru). Julia Podschipkova, Volgograd State University, Volgograd,student (podschipkova1993@mail.ru).
Abstract: A system is introduced for ecological-economical management of the territory of the Volga-Akhtuba floodplain. The system is built on the basis of hydrodynamic computer simulations and geo-informational technologies. We optimize the flood hydrograph for the Volga hydroelectric station, suggest the optimal plan of the system of channels for Volga-Akhtuba floodplain recovery, and develop control mechanisms for local economic entities.
Keywords: computer simulation, hydrological mode, hydrograph, digital terrain model, the Volga-Akhtuba floodplain, ecological-economical management.
Статья представлена к публикации членом редакционной коллегии Д.А. Новиковым
Поступила в редакцию 17.11.2014.
Опубликована 31.05.2015.
