Математическая модель распределения водных ресурсов трансграничных рек

Установлен характер изменения высоты и длины образующихся гряд в процессе их развития.

Определена скорость перемещения гряд и степень ее снижения по мере развития гряд.

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

1. Знаменская Н. С. Донные наносы и русловые процессы. - Л.: Гидрометеоиз-дат, 1976. - 192 с.

2. Волынов М. А. Пропускная способность саморегулирующихся речных русел // При-родообустройство. - 2011. - № 5. - С. 66-71.

3. Kennedy J. F. The formation of sediment ripples, dunes and antidunes // Annual revive of fluid mech. - 1969. - Vol. 1. - Р. 147-168.

4. Гришанин К. В. Устойчивость русел рек и каналов. - Л.: Гидрометеоиздат, 1974.

5. Волынов М. А. Внутрирусловые грядовые образования на речных излучинах // Научное обозрение. - 2012. - № 3. - С. 60-71.

6. Михайлова Н. А. Перенос твердых частиц турбулентными потоками воды. -Л.: Гидрометеоиздат, 1966. - 232 с.

7. Fudjiko Isaya. An experimental study dune development and its effect on s8. edi-ment suspension // Environm. Res. Center Univ. of Tsukuba. - 1984. - № 5.

8. Волынов М. А. О плоскости раздела радиальных осредненных скоростей в поперечном сечении речного потока на излучине: Мелиорация и окружающая среда: сб. трудов ВНИИГиМ. - М.: ВНИИГиМ, 2004. - Т. 2. - С. 99-101.

Материал поступил в редакцию 26.09.12. Волынов Михаил Анатольевич, кандидат технических наук, доцент Тел. 8 (499) 153-21-33 E-mail: v1532133@yandex.ru

УДК 502/504 : 556.16

Г. Х. ИСМАЙЫЛОВ, И. В. ПРОШЛЯКОВ, А. В. ПЕРМИНОВ

Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Московский государственный университет природообустройства»

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ РАСПРЕДЕЛЕНИЯ ВОДНЫХ РЕСУРСОВ ТРАНСГРАНИЧНЫХ РЕК*

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

Трансграничные реки, водные ресурсы, водноресурсная система, двухуровневая

модель, имитационная модель, каскад водохранилищ.

There is proposed a methodology of building a two-level model of water resources management of transboundary rivers on the basis of optimization and simulation systems. For the problem decision of distribution of water resources between sovereign states in view of interests of natural complexes the principle of decomposition is used.

Transboundary rivers, water resources, water resource system, two-level model,

simulation model, cascade of water reservoirs.

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

а Работа выполнена при финансовой поддержке (ВРТР), текуЩих в естественных руслах РФФИ, грант 12 - 05 - 00193а. и питающихся за счет поверхностного и

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

В мире насчитывается около 250 крупнейших трансграничных рек, воды которых используются нескольким государствами. Это Дунай, Нил, Нигер, Меконг и др. Россия, располагающая богатейшими водными ресурсами, также имеет около 70 полноводных трансграничных рек. Так, по Вуоксе проходит граница с Финляндией, по Неману - с Литвой, по Днепру - с Белоруссией и Украиной, по Западной Двине - с Белоруссией и Латвией, по Самуру - с Азербайджаном, по Волге, Уралу, Иртышу - с Казахстаном, по Селенге - с Монголией, по Амуру - с КНР, по Туманной - с КНР и КНДР.

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

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

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

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

Как видно из рисунка, первый тип взаимосвязей - координирующие воздействия V., поступающие от вышестоящей к нижестоящим управляющим системам. Под влиянием координирующих воздействий V. системы нижнего уровня вырабатывают локальные плановые решения х. и сигналы обратной связи X. . Передача наверх этих сигналов порождает второй тип взаимосвязей между управляющими системами. В двухуровневой системе управления управляющее решение х вырабатывается композицией локальных решений х = {х.} (¿ = 1,л) отдельных суверенных государств или субъектов этих государств, расположенных в бассейне реки. Вышестоящая система непосредственно не участвует в разработке решения х. , она влияет на него с помощью координирующих воздействий V . Таким образом, в рассматриваемых двухуровневых системах нахождение планов использования водных ресурсов х. и сигналов обратной связи X каждой нижестоящей локальной системой производится с помощью решения некоторой оптимизационной задачи Ь. V) , которую будем называть локальной оптимизационной задачей. Условимся также, что для выработки координирующих воздействий V = (у.1) используется координирующая модель

L0 (l). Следовательно, для управлении водными ресурсами трансграничных рек двухуровневая система моделей строится с целью расчленения сложной глобальной оптимизационной задачи на ряд локальных подзадач, увязку которых осуществляет координационная модель [2].

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

прямая задача - двойственная задача -

(С, Xmax; (Y, b) ^ min;

АХ < b;

DX > B; YA > C;

X > 0; Y > 0,

где X =

- n - мерный вектор-строка (эффективность); А

- матрица условий (удельные затраты всех видов производственных и природных ресурсов) размерности m х n; b - m - мерный вектор-столбец (производственные и природные ресурсы, необходимые для ведения сельскохозяйственного производства); Y = (y1,...,ym )- m - мерный вектор-столбец (переменная двойственной задачи - система расчетных цен, т. е. оценка производственных и природных ресурсов); D - матрица интенсивности (урожайности) производства, размерностью mD х n; b - mD -мерный вектор-столбец (объем необходимой сельскохозяйственной продукции, который предприятие получит за год).

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

Пусть Q - множество допустимых

*

ПИ-планов; Q - множество оптимальных ПИ-планов; Т* - множество допу-

*

стимых ПИ-систем расчетных цен; Т -множество оптимальных ПИ-систем расчетных цен. Тогда задача (1) приводится к следующему виду:

Q = {x / АХ < b,x > 0};

* I * * * / * \ / \1

Q =|x /x ей ,

Т = {у / YA > C, y > 0};

Т* = |y* / y* е Т,

Предположим, что ПИ-задача разрешима, т. е. оптимальный ПИ-план существует, Q^{0} ({0} - пустое множество). Известно, что в этом случае также существуют оптимальные ПИ-системы расчетных цен Т ^{0}. Более того, согласно первой теореме двойственности, максимальная величина целевой функции прямой задачи и минимальная величина целевой функции двойственной задачи равны и их общая величина есть оптимум F ПИ-задачи:

F = max (с, x) = min (у, b) =

xeQ

уеТ

=

Разобьем матрицу А на подматрицы Ар..., А.,..., АJ, где при каждом j е J матрица А. имеет размерность т х п. Речь идет о делении бассейна реки на отдельные подсистемы, каждая из которых является относительно самостоятельной; в то же время она взаимосвязана с некоторыми общими ресурсами, прежде всего с водными. Соответственно разбиваются: матрица Б = (Б,..., Б.,..., DJ ), векторы С = ( С. _ С.,..., CJ ) и X = ( _ X.,..., XJ ), причем размерность С.,X. для j е J равняется п. [2]. Тогда ПИ-задача преобразуется к следующему виду:

Е С Xj)

^ max

Е АХ < b

, j j

Е D.X. > B

1 j j

Xj > 0, j е[1, j]

г и

(Yj, bj min

Y1A1 > C1

YA > C

j j j

yj aj > cj

y> 0, j е[1, J].

Вводим m-мерный вектор-столбец u.,

j е [1, J] , удовлетворяющий условию

j

Еu. < b.

Скомпонованный план называется

1

ПИ-планом (центральным планом). Вектор и есть план j-го подсистемного компонента центрального плана и. В дальнейшем задачи (4) и (5) с ограничением (6) сводятся к задаче, для решения которой можно использовать теоретико-игровую модель [3]. Для численного решения рассматриваемых задач используется симплексный алгоритм линейного программирования.

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

Имитационная модель функционирования ВРС. Общая постановка задачи функционирования водно-ресурсной системы с водохранилищами формулируется следующим образом [1].

Рассмотрим водно-ресурсную

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

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

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

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

Математическая постановка рассматриваемой задачи такова: требуется минимизировать функционал -

Ф(V, и, £) = шшМ

V = АЖ + Ви;

V < V < V;

— —

т и - и

I Ь орЬ

и

1 = 0 орЬ

и > 0 при г = 0, V = У0,

где V - вектор наполнения; и - вектор попусков из водохранилищ; иорЬ - оптимальное значение по-пуско=; Ж - вектор водных ресурсов; £ - текущее время; А и В - матрицы системных условий.

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

ко)

№ 1' 2013

функционирования каскада водохранилищ многоцелевого назначения. Ограничения (9) и (10) применительно к ВРС с N взаимосвязанными водохранилищами и накопителями высокоминерализованного стока возвратных вод описывают:

динамику воды в водохранилище в момент времени £

ак.

dt

- Q.(t) + p.(t) + U. (t) + r.(t) =

= 0, V.(t0) = V.0; i = 1, N

V = К (Ч, Н), i = 1, N

где £ - время; К,(£) - объем воды в г -м водохранилище в момент £; <£,(£) = I - объем главных и боковых притоков к водохранилищу, включая сток возвратных вод из накопителей в момент £; Г(£), и(£), г(£) - соответственно водопотребление, попуск из водохранилища, потеря воды на испарение, фильтрацию и ледообразование из водохранилища в момент £; - начальный объем водохранилища в момент £0; н - глубина водохранилища;

кинетику процессов смешения солей в водохранилище в момент времени £:

а

—1— - Е юсч («> - ю+р, ащ а)+

а* 5=1

+ = М1

где М,(£) = С.(£) • К,(£) - масса солей в ¿-м водохранилище в момент £; С. - минерализация притока; С. - минерализация воды, вытекающей из водохранилища; тйп(£) - количество солей, поступающих в момент времени в водохранилище с подземным стоком через донные отложения; а - параметр, учитывающий самоочищение; М0 -масса солей в водохранилище в момент £0;

динамику воды в накопителе в момент £ можно описать следующим дифференциальным балансовым уравнением -

dV.c

dt

0, V.c (to) = V.c°;

- qC(t) + U.c(t) + r.c(t) =

= 1, N, V.c = V.c (qC, H),

Уравнения (11).(13) записываются в конечно-разностном виде.

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

Заключение

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

i = 1, N

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

1. Воропаев Г. В., Исмайылов Г. Х., Федоров В. М. Проблемы управления водными ресурсами Арало-Каспийского региона. - М.: Наука, 2003. - 427 с.

2. Волконский В. А. Оптимальное планирование в условиях большой размерности (итеративные методы и принципы декомпозиции) // Экономика и математические методы. - 1956. - Т. 1. - № 2. - С. 195-219.

3. Робинсон Дж. Итеративный метод решения игр // Матричные игры. - 1961. - С. 110-117.

Материал поступил в редакцию 03.05.12. Исмайылов Габил Худуш оглы, профессор, доктор технических наук, зав. кафедрой «Гидрология, метеорология и регулирование стока». Тел. 8 (499) 976-23-68 E-mail: gabil-1937@mail.ru Прошляков Игорь Валентинович, кандидат технических наук, профессор кафедры «Гидрология, метеорология и регулирование стока» Тел. 8 (499) 976-17-45 E-mail: iprosh@mail.ru Перминов Алексей Васильевич, кандидат технических наук, доцент кафедры «Гидрология, метеорология и регулирование стока» Тел. 8 (499) 976-17-45 E-mail: alexperminov@gmail.com