Решение простейшей задачи ВЭР1

Дано: основные исходные данные п.5.1.; Dt=2,63·10⁶ с; Z_вб^н = 103,0 м;`Q<sub>нб</sub><sup>треб</sup> = 500 м<sup>3</sup>/с;`Q_пр = 200 м³/с.

Требуется найти: Z_вб^к,`Z<sub>нб</sub>,`Н_а,`Q<sub>ГЭС</sub>,`N_ГЭС.

Задача решается безитерационным путем по алгоритму, представленному на рис.5.5. в последовательности, представленной в табл.5.5.

 

Решение простейшей задачи ВЭР2

Дано: основные исходные данные п.5.1.; Dt=2,63·10⁶ с; Z_вб^н = 103,0 м; Z_вб = 102,0 м;`Q_пр = 200 м³/с.

Требуется найти:`Z<sub>вб</sub>,`Z_нб,`Н<sub>а</sub>,`Q_ГЭС,`N_ГЭС.

Задача решается безитерационным путем по алгоритму, представленному на рис.5.6. в последовательности, представленной в табл.5.6.

 

Дано: Dt, Z_вб^н,`Q<sub>нб</sub><sup>треб</sup>,`Q_пр при`Q<sub>нб</sub><sup>треб</sup> >`Q_пр.

¯

`Z<sub>нб</sub> =`Z_нб(`Q_нб^треб)

¯

`Q<sub>в</sub> =`Q_нб^треб –`Q_пр;

¯

DV_срб =`Q_в·Dt

¯

V_в^к = V_в^н(Z_вб^н) – DV_срб

¯

Z_вб^к = Z_вб^к(V_в^к)

¯

`Н<sub>а</sub> = 0,5·(Z<sub>вб</sub><sup>н</sup>+Z<sub>вб</sub><sup>к</sup>)–`Z_нб

¯

Q_ГЭС^пред =

¯

`Q_ГЭС =

¯

`N<sub>ГЭС</sub> = k<sub>N</sub>·`Q_ГЭС·`Н_а

 

Рис.5.5

 

 

 

Дано: Dt, Z_вб^н, Z_вб^к,`Q_пр при Z_вб^н > Z_вб^к.

¯

DV_срб = V_вб^н(Z_вб^н)– V_вб^к(Z_вб^к)

¯

`Q_в =

¯

`Q<sub>нб</sub> =`Q_в +`Q_пр

¯

`Z<sub>нб</sub> =`Z_нб(`Q_нб)

¯

`Н<sub>а</sub> = 0,5·(Z<sub>вб</sub><sup>н</sup>+Z<sub>вб</sub><sup>к</sup>)–`Z_нб

¯

Q_ГЭС^пред =

¯

`Q_ГЭС =

¯

`N<sub>ГЭС</sub> = k<sub>N</sub>·`Q_ГЭС·`Н_а

 

Рис.5.6

 

Решение простейшей задачи ВЭР3

Дано: основные исходные данные п.5.1.; Dt=2,63·10⁶ с; Z_вб^н = 103,0 м;`N<sub>ГЭС</sub> = 100 МВт;`Q_пр = 200 м³/с; e_N=5%.

Требуется найти: Z_вб^к,`Z<sub>нб</sub>,`Q_в,`Q<sub>ГЭС</sub>,`Q_нб,`Н<sub>а</sub>,`N_ГЭС.

Задача решается итерационным путем с использованием алгоритма на рис.5.6 для задаваемых вариантов Z_вб^к, рассчитываемых значений`N<sub>ГЭС</sub><sup>расч</sup> и сопоставлением их с`N_ГЭС с точностью до e_N.

Вначале решается задача 5.2.4.: Z_вб^н = 103,0 м;`Q<sub>пр</sub> =`Q_нб = 200 м³/с;`Z<sub>нб</sub> = 81,00 м; Н<sub>ГЭС</sub> =Z<sub>вб</sub><sup>н</sup> –`Z_нб = 22,0 м;`N_ГЭС^расч = 8,2·22,0´ ´ 0,2 = 36 МВт, следовательно, требуется сработка водохранилища.

Расчет для нескольких вариантов Z_вб^к проводится в форме табл.5.6. с добавлением необходимых расчетных переменных и сопоставлением`N<sub>ГЭС</sub><sup>расч</sup>(Z<sub>вб</sub><sup>к</sup>) с заданным значением`N_ГЭС (см. табл.5.7).

Для поиска решения задачи используется графический способ, реализованный ранее в решении водохозяйственных задач на рис.4.1., 4.2. (см.рис.5.7).

На рис. 5.7 получена т.А с Z^к_вб= 101,85 м для которой выполнены все расчеты в табл.5.7, соответствующее искомому решению задачи ВЭР3.

 

 

Рис. 5.7

 

Обобщение задач ВЭР1 и ВЭР2

5.6.1. ВЭР1.1: работа ГЭС с Q_нб^треб > Q_ГЭС^пред (Zвб) с обеспечением N_ГЭС Þ max.

Дано: основные исходные данные приведены в п.5.1; Dt=2,63×10⁶ c; Z_вб^н=104,0 м; Q_пр=2200 м³/с; Q_нб^треб=2800 м³/с; e_Q=4%.

Требуется найти: Z_вб^к, Q_ГЭС, H_аi, N_ГЭС, Q_х.сбр., Q_в. Решение задачи реализуется итерационным путем и поясняется в табл.5.8 и на рис. 5.8.

Решение задачи: Z_вб^к=103,2 м; Q_ГЭС^пред=2335 м³/с; Q_в=623 м³/с; Q_х.сбр=465 м³/с; N_ГЭС^max= N_ГЭС^уст=345 МВт; Q_нб^факт=2800 м³/с.

 

5.6.2. ВЭР2.1: работа ГЭС предельными мощностями с Q_х.сбр¹0.

Дано: основные исходные данные приведены в п.5.1; Dt=2,63×10⁶ c; Z_вб^н=102,0 м; Z_вб^к=100,0 м;Q_пр=2100 м³/с.

Требуется найти: Q_вi, Q_ГЭС^пред, Q_х.сбр, N_ГЭС^пред, H_а.

Решение задачи реализуется безитерационным путем и поясняется в табл.5.9.

 

Рис. 5.8

 

Обобщение задачи ВЭР3

5.7.1. ВЭР3.1: сработка водохранилища по критерию N_ГЭС Þ max в зонах I и II.

Для зоны II дано: основные исходные данные приведены в п.5.1; Dt=2,63×10⁶ c; Z_вб^н=104,0 м; Q_пр=2000 м³/с.

Требуется найти: Z_вб^к, Q_в, Q_ГЭС, N_ГЭС, Q_х.сбр., Q_в по критерию N_ГЭС Þ max при e_Q=2%.

Поиск решения осуществляется в зоне II методом простых итераций и поясняется в табл.5.10 и на рис.5.9.

Координаты решения: Z_вб^к=103,65 м; Q_в=296 м³/с; Q_ГЭС^пред=2300 м³/с; N_ГЭС^пред= N_ГЭС^уст=345 МВт.

Для зоны I дано: основные исходные данные приведены в п.5.1; Dt=2,63×10⁶ c; Z_вб^н=101,7 м; Q_пр=2200 м³/с.

Требуется найти: Z_вб^к, Q_в, Q_ГЭС, N_ГЭС по критерию N_ГЭС Þ max при e_Q=2%.

Поиск решения осуществляется в зоне I методом простых итераций и поясняется в табл.5.11.

Координаты точки решения: Z_вб^к=100,5 м; Q_в=365 м³/с; Q_ГЭС^пред= Q_нб=2570 м³/с; N_ГЭС^пред= 322 МВт, Н_а=15,32 м.

Рис. 5.9

 

5.7.2. ВЭР3.2: Расчет режима наполнения водохранилища по критерию N_ГЭС Þ max в зоне I.

Дано: основные исходные данные приведены в п.5.1; Dt=2,63×10⁶ c; Z_вб^н=100,0 м; Q_пр=3000 м³/с; Z_вб⁰=101,83 м.

Требуется найти: Z_вб^к, Q_в, Q_ГЭС, Q_нб, H_а, N_ГЭС по критерию N_ГЭС Þ max при e_Q=2%.

Поиск решения осуществляется итерационным путем методом простых итераций и поясняется расчетами в табл.5.12 и на рис.5.10.

 

 

Рис. 5.10

 

Так как Q_пр > Q_ГЭС^max, то реализуется режим наполнения в зоне I (Z_вб^н < Z_вб⁰) с переменными значениями Q_ГЭС^пред(Z_вб) и N_ГЭС^пред(Z_вб).

Предварительно определяются предельные параметры ГЭС при Z_вб^н=100,0 м: Q_ГЭС^пред=2485 м³/с; Z_нб(Q_ГЭС^пред)=85,70 м; Н_а=14,3 м; N_ГЭС^пред = 8,2× Q_ГЭС^пред × Н_а=291 МВт. Соответственно: Q_в^н = Q_пр - Q_ГЭС^пред =515 м³/с. Это значение и принимается в табл.5.12 за Q_в¹, для которого и выполняются все расчеты для первой итерации в следующей последовательности:

 

(DV_нап¹ = Q_в¹×Dt) ® (V_в^к1 = V_в^н(Z_вб^н) + DV_нап¹) ® Z_вб^к1 ®

 

® Z_вб¹® Q_ГЭС^пред1(Z_вб¹) ® (Q_в^1расч = Q_пр - Q_ГЭС^пред1) ®

 

® [ dQ=| (Q_в¹- Q_в^1расч)/ Q_в¹| ] £ e_Q.

 

Если dQ £ e_Q, то решение найдено. В противном случае следует задаваться новым значением Z_вб^к2¹ Z_вб^к1 и продолжать поиск решения в последовательности, указанной в табл.5.12. Из результатов первой итерации видно, что Q_в>Q_в^расч, т.е. Z_вб^к2 должно быть меньше, чем Z_вб^к1=101,79 м, что и сделано в табл.5.12. Поиск Z_вб^к2 можно даже реализовать и для Q_в² = 0,5×(Q_в¹+ Q_в^расч1) = 480 м³/с.

Для двух итераций в табл.5.12 получены условия, говорящие о том, что решение задачи находится в интервале Z_вб^к от 101,50 до 101,79 м. Поиск решения (т.А) показан на рис.5.10 и в табл.5.12. Координаты точки решения: Z_вб^к=101,62 м; Q_в=452 м³/с; Q_ГЭС=2548 м³/с; Z_нб(Q_нб)=85,76 м; Z_вб=100,81 м; Н_а=15,05 м и N_ГЭС =314 МВт. При этом Z_вб^к < Z_вб⁰, т.е. решение задачи реализовано в зоне I корректно для заданного Dt и предположением, что линия ограничения Q_ГЭС^пред(Z_вб) в зоне от Z_вб^н до Z_вб^к - линейна. Если же указанная характеристика - нелинейна, то процесс поиска решения остается тот же, однако при этом следует проводить расчеты с Dt₁< Dt. Например, вместо месячного интервала Dt=2,63×10⁶ c следует брать декадные или Dt₁=1/3×Dt.

 

 

5.7.3. ВЭР3.3: Расчет режима наполнения водохранилища по критерию N_ГЭС Þ max в зоне II.

Дано: основные исходные данные приведены в п.5.1; Dt=2,63×10⁶ c; Z_вб^н=102,0 м; Q_пр=3000 м³/с; Z_вб⁰=101,83 м

Требуется найти: Z_вб^к, Q_в, Q_ГЭС, Q_нб, H_а, N_ГЭС по критерию N_ГЭС Þ max при e_Q=2%.

Поиск решения осуществляется итерационным путем методом простых итераций и поясняется расчетами в табл.5.13 и на рис.5.11.

 

Рис. 5.11

 

Как и в предыдущем случае при Q_пр > Q_ГЭС^max реализуется режим наполнения, но уже в зоне II (Z_вб^н > Z_вб⁰) с переменными значениями Q_ГЭС^пред(Z_вб) при N_ГЭС^пред(Z_вб)= N_ГЭС^уст=сonst при Н_а>H_p^N.

Предварительно определяются предельные параметры ГЭС при Z_вб^н=102,0 м: Q_ГЭС^пред=2600 м³/с; Z_нб(Q_ГЭС^пред)=85,80 м; Н_а=16,19 м; N_ГЭС^пред = N_ГЭС^уст=345 МВт. Соответственно: Q_в^н = = Q_пр - Q_ГЭС^пред = 3000-2600 = 400 м³/с. Это значение и принимается в табл.5.12 за Q_в¹, для которого выполняются все расчеты в последовательности, указанной выше в п.5.7.2.

По результатам первой итерации получено, что Q_в¹< Q_в^1расч, т.е в качестве Z_вб^к2 следует брать значение, больше Z_вб^к1=103,1 м, что и реализовано в табл.5.13.

5.7.4. ВЭР3.4: Расчет режима наполнения водохранилища по критерию N_ГЭС Þ max при переходе из зоны I в зону II.

Дано: основные исходные данные приведены в п.5.1; Dt=2,63×10⁶ c; Z_вб^н=101,5 м; Q_пр=3000 м³/с; Z_вб⁰=101,83 м

Требуется найти: Z_вб^к, Q_в, Q_ГЭС, Q_нб, H_а, N_ГЭС по критерию N_ГЭС Þ max при e_Q=2%.

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

Действительно, для Z_вб^н=101,5 м имеем следующие предельные параметры режима ГЭС: Q_ГЭС^пред=2600 м³/с; Z_нб(Q_ГЭС^пред)=85,80 м; Н_а=16,19 м; N_ГЭС^пред = N_ГЭС^уст=345 МВт. Соответственно: Q_в¹ = Q_пр - Q_ГЭС^пред = 400 м³/с, DV_нап¹= Q_в¹×Dt=1,05×10⁹ м³; V_вб^к1 = V_вб^н(Z_вб^н) + DV_нап¹ =(4,92+1,05) ×10⁹ =5,97×10⁹ м³; Z_вб^к1=102,58 м.

Полученное значение Z_вб^к1=102,58 м больше, чем Z_вб⁰=101,83 м, т.е. произошел переход режима ГЭС из зоны I в зону II, что требует изменения порядка расчета.

В начале рассматривается расчет ВЭР для заданных значений Z_вб^н=101,5 м и Z_вб^к1= Z_вб⁰=101,83 м при Q_пр =3000 м³/с в последовательности рассмотренной выше и представленной в табл.5.14.

В табл.5.14 значение Dt₁ определяется по формуле:

 

Dt₁= [V_вб^к(Z_вб⁰) - V_вб^н(Z_вб^н)]/ Q_в, (5.5)

 

где

Q_в = Q_пр - Q_ГЭС^пред , (5.6)

при

Q_ГЭС^пред =0,5× [Q_ГЭС^пред.к(Z_вб⁰)+ Q_ГЭС^пред.н(Z_вб^н)]. (5.7)

 

Здесь: Q_ГЭС^пред.к(Z_вб⁰)= Q_ГЭС^max=2630 м³/с.

 

Во время режима наполнения водохранилища предельными мощностями в зоне I получаем следующие параметры режима:

 

Q_ГЭС^пред(Z_вб^н)=2600 м³/с; Z_нб(Q_ГЭС^пред)=85,90 м;

Н_а= Z_вб^н - Z_нб(Q_ГЭС^пред)=15,6 м; N_ГЭС^пред(Z_вб^н)=332 МВт; соответственно для Z_вб^к1= Z_вб⁰=101,83 м получаем, что с Q_ГЭС^пред= =Q_ГЭС^max = 2630 м³/с; Н_а= Н_р^N=16 м и N_ГЭС^пред= N_ГЭС^уст=345 МВт.

Далее рассматривается новая задача, где Z_вб^н1= Z_вб⁰=101,83 м; Dt₂=Dt - Dt₁=1,8×10⁶ с; Q_пр=3000 м³/с и требуется найти оконча-тельные параметры режима наполнения.

Расчет производится по алгоритму, аналогичному рассмотренному в п.5.7.3 для ВЭР3.3 с Dt₂¹Dt (см.табл.5.15, рис.5.12).

Координаты окончательного решения: Z_вб^к=102,66 м; Dt₂=1,8×10⁶ с; Q_ГЭС^пред = 2555 м³/с; Z_нб(Q_ГЭС^пред) = 85,77 м; Н_а= Z_вб - Z_нб=102,24 - 85,77 =16,47 м; N_ГЭС^пред= N_ГЭС^уст=345 МВт. Графики изменения Z_вб(t) и N_ГЭС^пред(t) показаны на рис.5.13.

 

Рис. 5.12

 

 

Рис. 5.13

 

5.7.5. ВЭР3.5: Расчет режима наполнения водохранили-ща с Z_вб^н до НПУ по двум способам: равномерный режим наполнения до НПУ за Dt и максимально быстрое наполнение до НПУ с Q_х.сбр.=0 за Dt₁<Dt.

Дано: основные исходные данные приведены в п.5.1; Dt=2,63×10⁶ c; Z_вб^н=103,0 м; Q_пр=4000 м³/с.

Требуется найти: Z_вб(t), Q_в(t), Q_ГЭС(t), N_ГЭС(t) по двум указанным выше критериям.

Режим равномерного заполнения от Z_вб^н=103,0 м до НПУ=104,0 м: последовательность безитерационного расчета представлена в табл.5.16. В связи с тем, что в этих расчетах при Q_х.сбр.¹0 и Q_нб(t)= Q_ГЭС(t)+ Q_х.сбр.(t) нельзя пользоваться данными табл.5.3.

Режим максимально быстрого наполнения с Z_вб^н=103,0 м до НПУ с Q_х.сбр.=0 за Dt₂<Dt: по данным табл.5.16 имеем DV_нап^треб=1,84×10⁹ м³ и Q_ГЭС^пред(Z_вб=103,5 м)=2345 м³/с при Q_х.сбр.=0 по данным табл. 5.3.

Соответственно: Q_в = Q_пр - Q_ГЭС^пред=4000-2345=1655 м³/с и Dt₁=DV_нап^треб/ Q_в= 1,84×10⁹/1655=1,11×10⁶ с =12,84 суток.

Параметры полученного режима: Z_вб^н=103,0 м Z_вб^к= = НПУ=104,0 м; Q_нб=Q_ГЭС^пред= 2345 м³/с; Z_нб(Q_ГЭС^пред) =85,58 м; Н_а= Z_вб - Z_нб=17,92 м; N_ГЭС= N_ГЭС^уст=345 МВт; Dt₁=1,11×10⁶ с =12,84 суток.

В оставшийся период времениDt₂=Dt - Dt₁ = 1,52×10⁶ с ГЭС работает по водотоку с Q_нб= Q_пр и Q_в=0 при Q_х.сбр.¹0. Расчет этого режима работы представлен в табл.5.17.

 

Таблица 5.17

Расчет режима работы ГЭС "по-водотоку" при НПУ=соnst=104,0 м

Zвб, м	Qпр, м3/с	Qнб, м3/с	Zнб, м	Ha, м	QГЭСпред, м3/с	NГЭСпред, МВт	Qх.сбр, м3/с
104,0			86,85	17,15

В графическом виде полученный режим работы ГЭС по двум рассмотренным способам представлен на рис. 5.14 а) и б) соответственно.

 

Рис. 5.14

 

⇐ Предыдущая12345Следующая ⇒

Поделиться:

Воспользуйтесь поиском по сайту: