Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Резниковский, А. Ш. Управление режимами водохранилищ гидроэлектростанций

.pdf
Скачиваний:
23
Добавлен:
19.10.2023
Размер:
7.24 Mб
Скачать

Т а б л и ц а 3-4

Сопоставление параметров ряда среднегодовых расходов

 

По ряду наблюдений

По искусственному ряду

Параметр

 

 

 

в

1 000 лет

 

 

 

 

 

 

 

р.

Нарын

р. Енисей

р. Кура

р. Нарын р.

Енисей

р. Кура

Q , м 3/ с е к

431

1457

434

433

1460

435

c v

0,19

0,13

0,17

0,18

0,13

0,16

с .

0,34

0,30

0,48

0,03

0,10

0,18

Г

0,29

0,22

0,00

—0,01

0,03

0,06

метрии, как следует из работы [Л. 21], приводит к завы­ шению необходимой емкости водохранилищ, а (преумень­ шение коэффициента корреляции — к занижению ее ве­ личины, т. е. действие этих несоответствий параметров взаимно компенсирующее. Тем не менее для осторожно­ сти, видимо, целесообразно следовать рекомендации [Л. 21] и использовать данный метод разложения про­ цесса стока и его моделирования для расчетов сезонного регулирования, а для расчетов многолетнего регулиро­ вания— лишь в случае отсутствия коррелятивной связи между смежными членами ряда годового стока.

Для более общего случая рассмотрим еще один воз­ можный вариант разложения процесса стока на состав­

ляющие.

q(t)

принимаем величины,

зави­

В качестве базиса

сящие от предшествующих состояний процесса:

 

9(/) = Q '< = Q i - 5

] ( ^ ^

“ & -i)

(3-13)

Случайную часть ДQi, так же как и в предыдущем варианте, будем нормировать с помощью условных стан­ дартных отклонений для данного месяца

6(0 =

Qi Q'i

(3-14)

»i Vd/du ‘

 

В формулах (3-13) и

 

(3-14) D — определитель квад­

ратной симметричной

автокорреляционной матрицы

[см. уравнение (3-2)];

 

т — количество учитываемых

звеньев, т, е. длина автокорреляционной функции гидро­ логического ряда; Du и Дщ -д— алгебраические^ допол­ нения элементов ги и гш-п в определителе D\ Qi и о,- —

60

Норма и стандарт расходов данного месяца; Q i - j

месяцев.

норма и стандарт

расходов предшествующих

Способ моделирования непосредственно вытекает из

выражения (3-13)

 

 

Qi =

Q'i + 6 (0 з,

(3-15)

При наличии датчика нормальных чисел с параметра­ ми (0,1) его реализация на ЦВМ не представляет, боль­ ших трудностей. Для учета наблюденной асимметрично­ сти распределений остатков стока В' отдельные месяцы может быть использован прием, аналогичный предло­ женному в [Л. 83] и описанному выше (вывод зависимо­ стей, применимых для многозвенных цепей Маркова, произведен А. Ю. Александровским при консультации В. В. Атурина).

Выражение поправки на асимметричность имеет вид:

 

т

 

_

i=i

(3-16)

 

 

(D/Duf l 2

Здесь все обозначения аналогичны приведенным вы­ ше. Одним из основных в рассматриваемой схеме выде­ ления базисной части процесса является вопрос о коли­ честве учитываемых предшествующих состояний процес­ са т. Его выбор может осуществляться, например, с по­ мощью методов, описанных в [Л. 60], исходя из критерия минимума дисперсии прогноза. Вполне очевидным яв­ ляется стремление выбирать число т по возможности небольшим, т. е. начинать расчет для т —2 и лишь при недостаточно хорошем совпадении параметров перехо­ дить к большим т.

Этот вопрос для других целей и иными методами исследовался Н. А. Картвелишвили [Л. 37]. Он указыва­ ет, что «. . . увеличение т сверх четырех, а иногда сверх трех приводит не к улучшению, а к ухудшению совпаде­ ния. Это объясняется тем, что функции распределения вероятностей, по которым производится розыгрыш, не удовлетворяют точно так называемому обобщенному уравнению Маркова, т. е., строго говоря, являются функ­ циями перехода марковского процесса лишь с некоторым приближением».

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

61

8 результате разложения остатки (случайная состав­ ляющая) являлись с вероятностью 99,9%' некоррелиро­ ванной выборкой из одной генеральной совокупности. Они имели распределение, близкое к нормальному, с ну­ левым математическим ожиданием, единичной диспер­ сией и коэффициентом асимметрии, близким единице. Указанные результаты получены для нескольких значе­ ний величины т, т. е. при разном количестве учитывае­ мых в разложении предшествующих состояний процесса.

Рассмотрим результаты сопоставления параметров искусственных рядов и параметров исходных. Этот ана­ лиз проведем на примере трех значений величины т: 2, 3 и 12. Как и в предыдущем варианте, параметры единого ряда среднемесячных расходов совпадают с ис­ ходными. Для сечений одноименных месяцев получен­ ные результаты близки результатам предыдущего вари­ анта выделения составляющих и поэтому для сокраще­ ния места здесь не приводятся. Для ряда среднегодовых величин стока параметры даны в табл. 3-5.

 

 

 

 

 

Т а б л и ц а 3-5

 

Сопоставление параметров ряда среднегодовых

 

 

расходов р. Нарына

 

 

 

(без учета поправки на асимметричность)

Параметр

По наблю­

По одно­

По двух­

По трех-

По двенад­

денному

звенной

звенной

звеннсй

цатизвенной

 

ряду

цепи

цепи

цепи

цепи

Q

Ч з 1

433

431

431

427

С\

0 , 1 9

0 , 1 8

0 , 1 8

0 , 1 8

0 ,2 0

С.

0 , 3 4

0 , 0 3

0,21

0 ,2 2

0 , 1 5

Г

0 , 2 9

— 0,01

0 , 1 7

0 , 3 2

0 , 7 8

Приведенные данные, полученные при моделирова­ нии рядов без учета поправки (3-16) на асимметрич­ ность, свидетельствуют: 1) о нецелесообразности исполь­ зования в разложении большого числа звеньев в мар­ ковской цепи (см. выше); 2) о возможности с помощью выбора соответствующего количества учитываемых при моделировании предшествующих состояний процесса по­ лучить достаточно хорошее совпадение с заданными сле­ дующих параметров искусственных рядов годового сто­ ка: Q, Cv, г. Несколько хуже совпадают характеристики асимметрии. При использовании поправки (3-16), а для схемы простой цапи — поправки (3-11) результаты моде-

02

лирования 1000-летних рядов для рек Нарына и Енисея представлены в табл. 3-6. Приведенные данные свиде­ тельствуют о полном совпадении всех параметров годо­ вого стока при описании месячного стока рек трехзвен­ ной марковской цепью. Совпадение с исходными пара­ метров месячного стока также полное (для сокращения

места эти характеристики здесь

не приводятся). Для

р. Нарына не было

отмечено отрицательных

величин

в искусственном ряду

месячного

стока при его

модели-

Параметр

 

 

 

 

 

 

Т а б л и ц а

3-6

Сопоставление

параметров годового стока

 

 

(с учетом поправки на асимметричность)

 

 

По наблю­ денному ряду

Р- Нарын

 

По наблю­ денному ряду

Р-

Енисей

 

 

Простая цепь Маркова

Двухзвен­ ная цепь Маркова

Трехзвен­ ная цепь Маркова

Простая цепь Маркова

Двухзвен­ ная цепь Маркова

Трехзвен­ ная цепь

Маркова

Q , ж3! с е к

431

432

433

430

1457

1440

1460

 

1452

с„

0,19

0,19

0,18

0,18

0,13

0,13

0,12

0,13

C s

0,34

0,33

0,32

0,33

0,30

0,31

0,30

0,31

Г

0,29

О'.ОО

0,15

0,25

0,22

0,00

0,09

 

0,20

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

и[Л. 38]).

Ометодах моделирования речного стока для каска­

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

63

иые некоторой постоянной положительной величине (на­ пример, 7з Qш ш ) для практических расчетов вполне до­

пустимы.

Исследование закономерностей и моделирование слу­ чайно изменяющегося водопотребления и испарения

вэтом случае полностью аналогично изложенному выше

в§ 3-2, и дополнительных пояснений не требует.

3-4. Сопоставление водноэнергетических расчетов на основе разных методов моделирования

Для оценки различных способов моделирования гидрологических рядов проведем сопоставление результатов водноэнергетических и водохозяйственных расчетов, выполненных с использованием разных искусственных рядов, и сопоставим их с расчетами по календарному ряду наблюдений за речным стоном. Результаты таких сопоставле­ ний для Мингечаурского гидроузла приведены в [Л. 21], где был использован для моделирования искусственных рядов описанный в § 3:2 метод фрагментов.

Для определения влияния способа моделирования гидрологиче­ ских рядов на результаты водноэнергетических расчетов произведем сравнительные расчеты по искусственным рядам, смоделированным методом Монте-Карло с применением метода фрагментов и метода разностей.

Результаты расчетов показывают (табл. 3-7), что различия меж­ ду методами имеются, однако они сравнительно малы. В данном случае метод моделирования непосредственно месячных величин дает более надежный результат. Это объясняется тем, что при моделиро­ вании методом фрагментов в качестве исходных было взято сравни­ тельно небольшое количество фрагментов (10 моделей внутригодово­ го распределения из имевшихся 46). Это количество фрагментов, видимо, не является достаточным для полного совпадения результа­ тов расчета.

Такие же близкие результаты были получены при сопоставлении среднемноголетннх выработрк энергии ГЭС (табл. 3-8).

Таким образом, проведенное сопоставление показало, что исполь­ зование метода фрагментов и метода разностей при моделировании гидрологических рядов дает близкие результаты расчетов регулиро­ вания речного стока водохранилищем Мипгечаурской ГЭС. Распро­ странить указанный вывод па любые другие объекты, видимо, было бы преждевременным.

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

при моделировании

рядов, исходя из выборочных оценок (Q 0 и Cv0),

а также исходя из

этих оценок плюс или

минус их стандартная

ошибка (подробней — см. [Л. 2]).

 

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

раметров распределения число

перебоев нормальной отдачи может

изменяться очень существенно

(1727 и 530),

а обеспеченность нор-

64

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 3-7

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

 

 

 

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

 

 

 

 

 

 

Энергетическая отдача

 

 

 

 

Ирригационная отдача

 

 

 

 

нормальная

 

 

 

сниженная

 

 

нормальная

 

 

сниженная

 

Номера

Число

Обеспечен­

Число

Обеспечен­

Число

Обеспечен­

Число

Обеспечен­

столеток

 

перебоев

ность,

%1

перебоев

ность, %1

перебоев

ность, %1

перебоев

ность,

%1

.

А* 1 -

А

 

В

А

В

А '

В

А

А

А

В

А

В

А

В

1

20

34

98,3

97,2

0

0

99,9

99,9

140

176

88,3

85,3

27

35

97,8

97,1

2

40

13

96,7

98,9

0

0

99,9

99,9

185

108

84,6

91,0

51

18

95,8

98,5

3

20

16

98,3

98,7

0

0

99,9

99,9

164

103

86,3

91,4

30

27

9 7 ,5 ,

97,8

4

16

42

98,7

96,5

0

8

99,9

99,3

116

129

90,3

89,3

24

49

98,(Г 95,9

5

0

40

99,9

96,7

0

0

99,9

99,9

104

121

91,3

89,9

5

39

99,6

96,8

6

30

14

97,5

98,8

2

0

99,8

99,9

102

137

91,5

88,6

33

21-

97,3

98,3

7

29

3

97,6

99,8

0

0

99,9

99,9

129

103

89,3

91,4

34

14

97,2

98,8

8

3

2

99,8

99,8

0

0

99,9

99,9

108

103

91,0

91,4

12

12

99,0

99,0

9

23

27

98,1

 

97,8

0

0

99,9

99,9

121

150

89,9

87,5

26

42

97,8

96,5

10

21

7

98,3

99,4

0

0

99,9

99,9

163

109

86,4

90,9

41

14

96,6

98,8

Всего за

202

198

98,3

98,4

2

8

99,9

99,9

1332

1239

88,9

89,7

283

271

97,6

97,7

1 000 лет

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Наблюден­

 

0

99,8

 

0

99,8

 

46

91,5

 

0

99,8

ный 46-лет­

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ний ряд

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 Обеспеченность в процентах

определялась по числу бесперебойных месяцев.

,фрагментот‘ ;

В—ряды,

смоделированные разностным

* Л—ряды,

смоделированные

методом

Монте-Карло

с применением

метода

методом.

Т а б л и ц а 3-8

Сопоставление среднемноголетних выработок энергии

П о к а з а т е л ь

э , м л р д , к в т - ч

а*

 

 

 

 

 

 

а*

£

 

 

 

S

 

 

5

03

С

(Э)

С, О)

03

 

 

03

а

М

 

 

М

4

V

' '

 

"

2

**

 

 

 

£ К

(Ту о.

 

 

 

 

£

^

о 3

 

 

 

я й.

 

н

а

 

 

 

СТ) S

 

сту

£

Наблюденный ряд —

1,4341

46 лет

лет,

метод

1,4318

0,0210

0,015

—0,73

1,4641

1,3847

1 000

фрагментов

метод

1,4349

0,0200

0,013

—0,09

1,4655

1,4045

I 000

лет,

разностей

 

 

 

 

 

 

 

 

малыюй отдачи соответственно изменяться от 85,6 до 95,5%, т. е. значительно (расхождения выходят за пределы случайного рассея­ ния оценок обеспеченности отдач).

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

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

Г л а в а ч е т в е р т а я

ОСНОВНЫЕ ПОКАЗАТЕЛИ РЕЖИМА РАБОТЫ ГИДРОУЗЛОВ И ИХ ОЦЕНКА

4-1. Гарантированная отдача гидроузлов

Под г а р а н т и р о в а н н о й о т д а ч е й г и д р о у з л а подразумевается такая его отдача, которая с заданной надежностью (обеспеченностью) удовлетворяет требова­ ния водопотребителей и водопользователей в многолет­ нем разрезе [Л. 40].

66

Под г а р а н т и р о в а н н о й о т д а ч е й

в о д о х о з я й ­

с т в е н н о й с и с т е м ы подразумевается

некоторая -со­

вокупная отдача ее гидроузлов, которая с определенной надежностью (иной, -чем у отдельных гидроузлов) мо­ жет удовлетворять требования водо-потребителей и водо­ пользователей системы в многолетнем разрезе. При со­ хранении неизменной надежности (обеспеченности) абсо­ лютная величина отдачи системы возрастает по сравне­ нию с суммой гарантированных отдач отдельных гидро­ узлов при изолированной работе [Л. 21, 55, 67, 76].

В зависимости от назначения гидроузла его гаранти­ рованная отдача может быть выражена: 1) некоторой величиной и обеспеченностью поступающего в нижний бьеф расхода воды для нужд водного транспорта, рыб­ ного, лесного хозяйства и др.; 2) величиной забираемого из верхнего бьефа объема воды для ирригации, водо­

снабжения и пр.; 3) при наличии в

составе гидроузла

гидроэлектростанции — количеством

гидравлической

энергии, используемой в расчетных маловодных условиях для участия в покрытии графика нагрузки энергосисте­ мы во время прохождения ее максимума; 4) совокупно­ стью указанных величин в комплексных гидроузлах; 5) общей величиной гидравлической энергии ГЭС, полу­ чаемой от каскадов гидроэлектростанций, расположенных на разных реках и работающих в одной энергосистеме. При этом роль каждого гидроузла в системе и взаимо-'- связь режимов их работы регламентируется (или опти­ мизируется) с целью достижения максимальной величи­ ны совокупной гарантированной отдачи ГЭС системы. Оптимальность такого решения обосновывается тем, что величина совокупной гарантированной -отдачи ГЭС энер­ госистемы определяет собой необходимую величину за­ меняемых в ней мощностей тепловых (или атомных) электростанций. Сокращение же установленных мощно­ стей последних ведет к общему уменьшению капиталь­ ных вложений в системе. Это происходит, потому что указанное повышение' гарантированных мощностей ГЭС в системе достигается только за счет оптимизации режи­ мов их работы, т. е. без дополнительных вложений

вустановленные мощности гидроэлектростанций.

Вобщем случае гарантированная отдача гидроузлов определяется величиной полезной емкости их водо­

хранилища, размером и распределением притока

воды

в реке во времени. При сезонном регулировании

вели­

5*

67

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

При многолетнем регулировании стока величина га­

рантированной отдачи

определяется

притоком воды

к гидроузлам за расчетный я-летний

период (критиче­

ский период сработки

и наполнения

водохранилищ) и

полезной емкостью водохранилищ, срабатываемой и за­ тем наполняемой в течение этого периода времени. При известных параметрах речного стока и заданной надеж­ ности (обеспеченности) величина гарантированной отда­ чи может быть определена однозначно для каждого из рассматриваемых значений полезной емкости водохрани­ лища. При многолетнем регулировании стока для этой цели могут быть' использованы известные номограммы для определения полезной емкости водохранилищ мно­ голетнего регулирования стока [Л. 21, 40, 51 и др.]. Ког­ да в составе ‘ гидроузла имеется гидроэлектростанция, размер гарантированной отдачи, помимо указанных вы­ ше факторов, зависит также и от напора ГЭС. Расчет гарантированной энергоотдачи, особенно при многолет­ нем регулировании стока, производится по календарным гидрологическим рядам, наблюденным или искусствен­ ным, смоделированным методом Монте-Карло (Л. 21, 64]. Это дает возможность учесть снижение напора при сра­ ботке водохранилища и связанное с ним нарастание ин­ тенсивности расходования воды, аккумулированной в его полезной емкости, а также сложное взаимодействие гид­ роэлектростанций в энергосистеме.

При проектировании гидроузла его гарантированная отдача оценивается и экономически обосновывается для

каждого

варианта нормального

подпорного

уровня

(НПУ),

а экономический эффект от использования этой

отдачи в системе является исходным

материалом

при вы­

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

68

которого общего расчетного маловодного периода

[Л. 57, 76].

Для условий, когда выбраны основные параметры гидроузла (НПУ, полезная емкость водохранилища, установленная мощность ГЭС) и когда задана обеспе­ ченность его гарантированной отдачи, абсолютная вели­ чина ее, рассматриваемая как объем зарегулированного стока или энергии, определяется однозначно. Однако величина гарантированной отдачи, как правило, не рас­ ходуется равномерно в течение года, а изменяется в со­ ответствии с графиком потребления и в развивающейся системе может быть переменной также и из года в год.

Под понятием «развитие водохозяйственной системы»

подразумеваются

рост водопотребления

н связанные

с ним изменения

во времени (чаще всего

увеличение)

требований водо- и энергопотребителей к использованию водных ресурсов рек, а та^же изменения, происходящие в самой водохозяйственной системе: строительство но­ вых гидроузлов, начальное заполнение пх водохранилищ, ввод новых агрегатов на существующих или строящихся ГЭС и т. п. Для удовлетворения изменяющихся во вре­ мени требований потребителей гидроузлы работают в ре­ жиме, обеспечивающем необходимое перераспределение во времени'гарантированной отдачи. При проектирова­ нии гарантированная отдача распределяется во времени на основании данных о перспективе развития энерго­ системы в условиях эксплуатации в соответствии с кон­ кретной обстановкой в системе.

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

В развивающейся водохозяйственной и энергетической системе установленная мощность гидроэлектростанций

69

Соседние файлы в папке книги из ГПНТБ