Файл: Пояг, М. А. Комплексное использование и охрана водных ресурсов Молдавии.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 16.10.2024

Просмотров: 58

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

полнительный так называемый фиктивный водопотребитель, которому в качестве объема спроса приписывает­ ся разница между общим объемом воды и фактической потребностью в ней, то есть

вп+ 1 = 2

а 1 — S Bi — 5 161 — 4 967= 194 млн. м3.

i=i

j=i J

После введения фиктивного водопотребителя открытая водотранспортирующая задача преобразуется в закры­ тую модель.

Выполним первое распределение объемов воды мето­ дом северо-западного угла.

У первого водоисточника имеется 2 235 млн. м3, а пер­ вому водопотребителю требуется 3 798 млн. м3 воды, поэтому все 2 235 млн. м3 записываем первому водопо­ требителю. При этом для удобства водоподачу записы­ ваем в первом нижнем углу клетки.

Так как все запасы первого водоисточника распреде­ лены, приступаем к распределению запасов воды второ­ го водоисточника. Во втором водоисточнике имеется 993

млн. м3 воды,

а первому водопотребителю требуется еще

1 563 млн. м3

(3798—2235). Записываем первому водо­

потребителю 993 млн. м3 воды от второго водоисточника. При этом ему еще недостает 570 млн. м3 (1563—993). Снова записываем первому водопотребителю недостаю­ щие 570 млн. м3 воды от третьего водоисточника, а у по­ следнего остается еще 19 м3 (589—570). Эти 19 млн. м3 от третьего водоисточника передаем второму водопотре­ бителю, которому не хватает еще 837 млн. м3 (856—19). Они доставляются от четвертого водоисточника. У чет­ вертого водоисточника осталось еще 313 млн. м3 (1150—— 837), которые доставляются третьему водопотребителю, и тем самым полностью распределяются запасы воды четвертого водоисточника. Этим заканчивается распреде­ ление водных ресурсов;" все водопотребители удовлет­ ворены и вся вода распределена. Следовательно, полу­ чен первоначальный опорный план распределения потребностей, Удельные затраты воды в клетках с потреб­ ностями, то есть в занятых клетках, обведены квадрати­ ками.

Прежде чем проверить на оптимальность полученную схему распределения поставок воды, необходимо прове­ рить количество занятых клеток. Их должно быть

151


га+ П'—1, так как допустимое распределение поставок во­ ды можно выполнить только в том случае, если число за­ нятых клеток не меньше, чем m + n— 1. Однако допусти­ мое распределение по произвольному набору такого чис­ ла занятых клеток выполнить не всегда можно. Необхо­ димо выполнение и других условий:

а) в каждом столбце или строке должно быть не ме­ нее одной занятой клетки, а любая занятая клетка мо­ жет оказаться единственной в своей строке или столбце, но не может быть единственной одновременно и в столб­ це, и в строке;

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

строке (столбце) имеется

единственная занятая клетка

и если ее зачеркнуть, то

в столбце (строке), в .котором

находится эта клетка, должна остаться только одна за­ нятая клетка и т. д. Если последовательно проходя по столбцам и строкам матрицы удается вычеркнуть все без исключения занятые клетки, значит их комбинация вычеркиваемая.

Проверим, составляют ли занятые клетки в табл. 20 вычеркиваемую комбинацию. Начнем просмотр матри­ цы со строк. Вычеркнем занятые клетки 1,1 и 2,1, явля­ ющиеся единственными в своих строках. Переходим за­ тем на просмотр столбцов и вычеркиваем клетки 3,1 и 4,3 — единственные в своих столбцах. Снова просматри­ ваем матрицу по строкам и вычеркиваем клетки 3,2 и 4,2. Следовательно, занятые клетки расположены в порядке вычеркиваемой комбинации.

Если в матрице число занятых клеток будет меньше гп+ п—1, то имеют место случаи вырождения задачи и ее решение затрудняется. Для устранения случая вырож­ дения необходимо поменять местами строки или столб­ цы или же свести в некоторые клетки так называемые нулевые водоподачи. Для этого вписывается в опреде­ ленную клетку нулевая поставка и клетка считается за­ нятой. В нашем случае заполненных клеток 6 , то есть равно гп+ п— 1 (3+4—1). Следовательно, план не явля­ ется вырожденным и можно продолжить решение за­

дачи.

152

После первого распределения поставок подсчитыва­ ется общая стоимость транспортирования воды

ШШ

F(x) = 2 ^С,,-X,j = 0,1199x2235+0,1184x993 +

wJ=i

+0,1231 х 570+0,0211X 19+ 0,147X 837+0,0278 х 313 =

=268+118 + 72 + 0,4+12 + 8,6 — 479 млн. руб.

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

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

вает номер строки,

а второй индекс — номер столбца

данной клетки.

цепи

всех свободных клеток, по­

Ниже приводятся

строенные по данным табл.

2 0 , и их характеристики.

Для клетки 1,2

 

Для клетки 2,2

1-0,11991------+ 0,0181

|-0,11841------- + 0,0175 .

1+0,1231 |------ 1-0,0211 [

| + 0,12311------ 1-0,0211 |

д1>2=

+0,0181 — 0,0211 +

Л2о=: + °,0175 - 0,0211 +

 

+ 0,01231-0,1199=

+ 0,1231 -0 ,1 1 8 4 =

=

+0,0002.

 

= + 0 ,0 0 1 1 .

153


Д л я клетки

3,3

 

Для

клетки 4,1

1 -0 ,0 2 1 1 1------

+0,0256

-0,1231

| - — | + 0,0211 |

I

 

1

 

I

____ ]____

1+ 0,01471------1-

0,02781

+ 0,1201

------ 1 - 0,01471

Д2,3=

+0,0256 — 0,0278 +

Дм =

+0,0211 -0,0147 +

 

+0,0147-0,0211 =

 

+ 0,1201 -0,1231 =

=

-0,0086.

 

=

+0,0034.

 

Для клетки 2,3

 

 

 

1-0,1184 I---------------------

+ 0,0323

 

I

_________

 

 

 

1+ 0,123Ц----- 1-0,02111

 

 

 

 

 

I

| - 0,0278 1

 

1+ 0,01471

д, 13= 0,1231 — 0,1184 + 0,0147 — 0,0211 +0,0323 -

-0.0278 = +0,0028.

Теперь рассмотрим, что собой представляют характе­ ристики, вычисленные для каждой цепи.

Очевидно, для улучшения исходного распределения водотранспортирования необходимо проверить целесооб­ разность перемещения части запасов воды из занятых клеток в свободные. Если в клетку 1,2 запишем транс­ портирование хотя бы одного кубометра воды, то для сохранения общей балансировки водных ресурсов по строкам и столбцам придется в клетке 1,1 по той же строке уменьшить поставку на кубометр воды, и баланс воды по первой строке сохранится. Для сохранения же баланса воды по первому и второму столбцам придет­ ся в клетке 2,1 увеличить, а в клетке 2 ,2 уменьшить на кубометр воды. В результате такого перераспределения кубометра воды по первой цепи мы уменьшаем удельные затраты на 0,0009 руб. (0,0181—0,1199+0,1184—0,0175).

Следовательно, характеристика данной клетки показыва­ ет, что при включении в нее каждого кубометра воды удельные затраты снижаются на 0,0009 руб.

154


Характеристика клетки 1,3 показывает, что перенесе­ ние поставок в эту клетку вызывает уменьшение удель­ ных затрат в размере 0,0115 руб. ( + 0,0124—0,0181— —0,0323 + 0,0175) в расчете на кубометр воды.

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

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

Внашем случае отрицательная характеристика при­ надлежит клетке 2,3 (Л 2,3 = —0,0086). Поэтому по дан­ ной цепи и выполняется перераспределение поставок так, чтобы в свободной клетке появились поставки. Для этого вдоль всей цепи перемещается такой объем воды, который равен наименьшему объему водоподачи одной из ее занятых клеток, имеющей в цепи отрицательные удельные затраты. Этот объем вычитается из объемов водоподачи всех клеток данной цепи с отрицательными затратами и прибавляется к объему водоподачи всех клеток цепи с положительными затратами. В результате такого перемещения сохранился баланс воды по всем строкам и столбцам матрицы.

Вцепи, построенной для свободной клетки 2,3 пере­

мещается поставка, равная 313 млн. м3 воды. После этой поставки вдоль всей цепи получается следующее но­ вое распределение воды в клетках этой цепи: Х41=313— —313 = 0; Х23=0+313=313; X,, = 993—313=680; Х42= = 1150—294=856.

После этого составляется новая таблица распределе­ ния поставок воды с учетом выполненного перераспре­ деления по цепи 2 , 3(табл. 2 1 ).

В результате полученной новой схемы .распределения водных ресурсов приведенные затраты уменьшились на 2,7 млн. руб. (31ЗХ 0,0086). Все затраты по транспорти­ рованию воды водопотребителям составят 476,3 млн. руб. (479,0—2,7). Однако следует проверить, является ли полученная схема оптимальной. Для этого, снова со-

155

 

 

 

 

 

 

 

Т а б л и ц а

21

 

 

Второе базисное распределение водных ресурсов

 

 

 

 

 

Водопотребителн н нх потребности

 

.

Водоисточники

 

 

 

 

 

 

 

 

 

и их объемы

п,

 

 

 

п„

 

п ,

 

 

 

 

 

 

 

 

 

 

 

 

3 798

 

 

 

856

 

313

 

 

 

 

1 0,1199

1

 

0,0181

0,0214

 

м ,

2 2 3 5

2 235

 

 

 

 

 

 

 

 

 

 

1 0,1184

|

 

0,0175

1 0,0323

| .

Mj

993

680

 

 

 

 

 

313

 

 

 

 

| 0,1231

|

 

0,0211

0,0256

 

М»

589

589

 

 

 

 

 

 

 

 

 

 

1 0,1201

|

|

0,0147 |

0,0278

 

 

 

1 150

294

 

 

 

\

 

 

 

 

 

 

 

 

856

 

 

 

ставляем для каждой свободной клетки табл.

21

цепи и

определяем их характеристики.

 

 

 

 

 

 

Для клетки

1,2.

 

 

 

Для клетки 2,2.

 

| -0,1199]------

+0,0181

 

| — 0,11841------

+0,0175

I

 

1

 

...— Г.......

 

1

 

1+0,12011----- 1-0,01471

 

1+0,1201 1---- - 1-0,01471

Д,г =

+0,1201-0,1199 +

 

Д2,2 == 0,0175 —0,0147 +

 

 

+0,0181 -0,0147=

 

 

 

+6,1201 -0,1184 =

 

=

+0,0036.

 

 

 

=

+0,0045.

 

 

 

 

Для клетки 3,2

 

 

 

Для клетки 3,3

 

|-0,1231|------- +0,0211

 

1 +

0,1184|------ 1 — 0,0323

1+

I

 

I

 

1-

О,*!2311-------

+0.1)256

0,1201|------1— 0,01471

 

д3>2=

+ 0,0211 -0,0147+

 

д3>3=

+0,1184 -

0,1231 +

 

 

+ 0,1201-0,1231 =

 

 

 

+0,0256 -

0,0323 =

 

=

+ 0,0034.

 

 

 

=

-0,0114.

 

 

 

155


Д ля клетки

1,3

Д л я

клетки

4,3

1— 0,11991—----

+0,0214

| + 0,11841------ 1-

0,0323 |

I

 

____ !____

'

I

| +0,11841------1- 0,03231

[ -0,1201

|------ +0,0278

д13= + 0 ,1 1 8 4 - 0,1199 +

д4,3 =

+0,1184 — 0,0323 +

+ 0,0214 -0,0323=

 

+0,0278 - 0,1201 =

= - 0,0124. '

=

- 0,0062. .

Наличие отрицательных характеристик указывает на неоптимальность схемы и необходимость дальнейшего

ееулучшения.

Вцепи, построенной для клетки 1,3, перемещаем

объем наименьшей поставки из вершины с отрицательным тарифом, то есть 313. После перемещения этой поставки

вдоль цепи получим: Х2,3 = 313—313 = 0; Xi,з=

0 + 313=

= 313; Х2.1 = 6 8 0 + 313 —993; Хы = 2 235 - 3 1 3 = I 922.

Составляем новую таблицу распределения

поставок

с учетом перераспределения водных ресурсов в цепи клетки 1,3 (табл. 22).

 

 

 

 

 

Т а б л и ц а

22

 

Третье базисное распределение водных ресурсов

 

 

 

Водопотребнтели н их потребности

 

Водоисточники

 

 

 

 

 

и их объемы

п ,

 

п ,

п ,

 

 

 

3 793

 

856

313

 

 

|

0,1199

1

0.0181

| 0,0214

|

М,

2 235

I 922

 

 

313

 

 

1

0,1184

|

0,0175

0,0323

 

м 2

993

993

 

 

 

 

 

|

0,1231

|

0,0211

0,0256

 

М3

589

589

 

 

 

 

 

|

0,1201

|

1 0,0147 |

’ 0,0278

 

м<

1 150

294.

 

856

 

 

157