Комплексный расчет по строительству парогазового блока мощность 400 МВт

15. Цыганок А.П. Проект ТЭС (Часть 2): нормативные материалы к дипломному и курсовому проектированию для студентов специальностей 0301, 0305 - "Электрические станции", "Тепловые электрические станции"/ А.П. Цыганок, Н.А. Сеулин; КрПИ - Красноярск, 1981. - 36 с.

16. Михайленко С.А. Тепловые электрические станции: учеб. пособие. 2-е изд. испр. / С.А. Михайленко, А.П. Цыганок. - Красноярск: ИПЦ КГТУ, 2005. - 302 с.

17. Цыганок А.П. Тепловые и атомные электрические станции: учеб. пособие: в 2 ч./ А.П. Цыганок. - Красноярск: ИПЦ КГТУ, 2000. - 123 с.

18. Трухний А.Д. Стационарные паровые турбины. - 2-е изд., перераб. и доп./ А.Д. Трухний. - М.: Энергоатомиздат, 1990. - 640 с.: ил.

19. Костюк, А.Г. Турбины тепловых и атомных электрических станций: учебник для вузов. - 2-е изд., перераб. и доп. / А.Г. Костюк, В.В. Фролов, А.Е. Булкин и др.; Ред. А.Г. Костюк. - М.: Издательство МЭИ, 2001. -488 с.: ил.

20. Стерман Л.С. Тепловые и атомные электрические станции: учебник для вузов. - 3-е изд., перераб. / Л.С. Стерман. - М.: Издательство МЭИ, 2004. - 424 с., ил.

21. Цанев, С.В. Газотурбинные и парогазовые установки тепловых электростанций: Учеб. пособие для вузов / С.В. Цанев, В.Д. Буров, А.Н. Ремезов. - М.: Издательство МЭИ, 2002. - 584 с.

22. Трухний, А.Д. Расчет тепловых схем парогазовых установок утилизационного типа / А.Д . Трухний, С.В. Петрунин // МЭИ.- 2001.-21 с.

23. Куликов С.М., Бойко Е.А. Расчет содержания вредных веществ в дымовых газах при проектировании котлов и энергетических установок: Методические указания по дипломному проектированию для студентов специальности 1005 - "Тепловые электрические станции", 1007 - "Промтеплоэнергетика". - Красноярск, КГТУ, 1995.

24. Бойко Е. А. Котельные установки и парогенераторы (тепловой расчет парового котла): Учебное пособие / Е. А. Бойко, И. С. Деринг, Т. И. Охорзина. Красноярск: ИПЦ КГТУ, 2005. 96 с.

25. Орлов, К. А. Исследование схем парогазовых установок на основе разработанных прикладных программ по свойствам рабочих тел [Текст]: автореф. дис. …канд. техн. наук: 05.14.14 / Орлов Константин Александрович. - М., 2004. - 32 с.

26. Подбельский В. В. Язык С++: Учеб. пособие. - 5-е изд. - М.: Финансы и статистика, 2005. - 560 с.: ил.

27. Цыганок, А.П. Проектирование тепловых электрических станций: учебное пособие / А.П. Цыганок, С.А. Михайленко. Красноярск: ИПЦ КГТУ, 2006, 136 с.

28. Энергетические и теплотехнические процессы и оборудование: сб. науч. тр. / ред. А. В. Бойко. - Харьков: НТУ ХПИ, 2008. - 196 с.

29. Кругликов, П. А. Технико-экономические основы проектирования ТЭС и АЭС: письменные лекции. - СПб, СЗТУ, 2003. -118 с.

30. Энергетическое машиностроение - новые решения: сб. конф. / ред. А. А. Бельтюков. - Екатеринбург, 2007. - 96 с.

31. Сокращение потребления природного газа и перспективы электроэнергетики: "атомный" и "парогазовый" сценарии: доклад. / И. В. Бабанин, В. А. Чупров. - М. , 2005 - 18 с.

32. Паровые турбины: номенклатурный каталог. / ПО "Силовые машины" СПб, 2005 - 78 с.



Приложение

/*Подпрограмма теплового расчета газотурбинной установки*/

//---------------------------------------------------------------------------

#include <stdio.h>

#include <conio.h>

#include <math.h>

#include <eheat.h>

#include <emath.h>

void main()

{

float

/*расход сухого воздуха на входе в камеру сгорания

с учетом охлаждения, кг\с*/G_k_ohl=309.431,

/*теоретически необходимое количество воздуха, кг/кг*/ L_0=15,

/*расход топлива на ГТ, кг/с*/ B=6.697,

/*расход газов на выходе из ГТ с учетом охлаждения, кг/с*/G_t_ohl_=359.463,

/* расход газов на выходе из ГТ без учета охлаждения, кг/с */G_t_ohl=316.652,

/*температура воздуха за компрессором, ?C*/ t_b=389.52,

/*температура газов за ГТ,? C*/ t_d=584.834,

/*мощность паровой турбины, кВт*/ N_e_pt=150000,

/*термический КПД ПТУ*/ n_t_ptu=0.418,

/*внутренний относительный КПД паровой турбины*/ n_oi=0.85,

/*механический КПД*/n_m=0.98,

/*КПД генератора*/ n_g=0.99,

/*КПД НПГ*/n_pg=0.89,

/*КПД транспорта*/n_tr=0.98,

/*теплота сгорания газа, kDzh/kg*/Q_n_r_gas=44300,

/*содержание кислорода в окислителе после добавки

воздуха от вентилятора в топку НПГ, %*/O2_t_pg=15,

/*коэффициент избытка воздуха в топке НПГ*/alfa_t_pg=1.2,

/*давление острого пара, бар*/P0=130,

/*температура острого пара, ?C*/t0=540,

/*давление в конденсаторе, бар*/Pk=0.035,

/*коэффициент регенерации*/K_p=1.01,

/*температура питательной воды, ?C*/ t_pv=270,

/*теплоёмкость воды, кДж/(кг*?C)*/ Cp=4.186,

/*температура газов за экономайзером, ?C*/t11_ek=352,

/*теплоёмкость топлива НПГ, кДж/(кг*?C)*/c_tl=2.5,

/*температура топлива НПГ, ?C*/t_tl=70,

/*потеря теплоты от наружного охлаждения*/q5=0.005,

/*КПД камеры сгорания*/ n_ks=0.995,

/*теплоёмкость добавочного воздуха за

вентилятором, кДж/(кг*?C)*/c_dv=1.00464,

/*температура добавочного воздуха за

вентилятором, ?C*/ t_dv=150,

/*энтальпия добавочного воздуха, кДж/кг*/i_dv=2582.762,

/*давление в деаэраторе, бар*/P_d=7,

/*KPD teploobmennika*/n_to=0.98,

/*температура добавочного воздуха до ВЗП, ?C*/t_dv0=45,

/*температура основного конденсата за ГВПНД, C*/t1_ok=160,

/*давление пара в отборе на деаэратор, бар*/P_otb=12.9,

/*температура основного конденсата после точки смешения, C*/t_ok=60,

/*температура ХОВ, ?C*/t_xov=30,

/*КПД питательного насоса*/n_pn=0.75,

/*перепад давления в ПН, бар*/dP_pn=173,

/*электрическая мощность ГТ, кВт*/N_el_gt=110000,

/*электромеханический КПД*/n_em=0.98;

float alfa_ux_gt, G_oxl, t_ux_gt, v0_ks, O2_ux_gt, n_ptu, B_r_pg_ut, B_r_pg_nt, G_pg, G1_pg, G_dv, q_dv, Q_pg, D_ne, i_ne, i_k, D_np, q_pv, H_pt, P_b, i_np, I4_sv, s_para, I4_sg, I_ek_sg, q_tl, I_ek_sv, I_ek_sg_pg, alfa_yx_pg, alfa_yx_pg_dv, Bp1_pg, h_dv, t1_vzp, I1_vzp_sv, I1_vzp_sg, h_d, D_pv, Q_gvpvd_g, Q_gvpvd_v, h_dv0, Q_vzp_g, Q_vzp_vozd, t11_vzp, I11_vzp_sg, I11_vzp_sv, D_d, D_ok, i0_otb, i_otb, t_ux, I_ux_sv, I_ux_sg, h_ok0, Q_gvpnd_g, Q_gvpnd_v, D_rec, D1_ok, D_ezh, D_upl, D_ut, D_ne0, i1_np, i_c1, D_c1, D_np1, D_sn, D_xov, t1_xov, v_pv, h_d0, Wep, d, N_sn_pt, N_sn_gt, N_sn_pgu, N_pgu_br, N_pgu_nt, n_pgu_br, n_pgu_nt;

/*массив зависимости удельного объёма воздуха от температуры*/

float mass_v0[2][21] = { {0, 50, 100, 150, 200, 250, 300, 350, 400, 450, 500, 550, 600, 650, 700, 750, 800, 850, 900, 950, 1000}, {0.774, 0.916, 1.057, 1.199, 1.340, 1.484, 1.626, 1.767, 1.908, 2.049, 2.193, 2.336, 2.475, 2.618, 2.762, 2.899, 3.040, 3.185, 3.322, 3.472, 3.61}};

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

float mass_I_sv[2][21] = { {0, 50, 100, 150, 200, 250, 300, 350, 400, 450, 500, 550, 600, 650, 700, 750, 800, 850, 900, 950, 1000},

{0, 50.19, 100.56, 151.21, 202.22, 253.67, 305.61, 358.13, 411.22, 464.91, 519.2, 574.09, 629.58, 685.58, 742.13, 799.17, 856.68, 914.62, 972.97, 1031.68, 1090.75}};

/*массив зависимости энтальпии продуктов сгорания от температуры */

float mass_I_sg[2][21] = { {0, 50, 100, 150, 200, 250, 300, 350, 400, 450, 500, 550, 600, 650, 700, 750, 800, 850, 900, 950, 1000},

{0, 53.59, 107.78, 162.60, 218.10, 274.33, 331.31, 389.08, 447.65, 507.03, 567.20, 628.19, 689.93, 752.41, 815.60, 879.47, 943.98, 1009.10, 1074.80, 1141.04, 1207.81}};

alfa_ux_gt=(G_k_ohl*0.992)/(L_0*B); /*коэффициент избытка воздуха в сбросных газах ГТ*/

G_oxl=G_t_ohl_-G_t_ohl; /*расход воздуха на охлаждение*/

t_ux_gt=(G_t_ohl*t_d+G_oxl*t_b)/G_t_ohl_; /*температура смеси выхлопных газов и воздуха за ГТ*/

int i;

for (i=0; i<21; i++)

{

if (t_ux_gt>mass_v0[0][i])

{v0_ks=mass_v0[1][i]+((mass_v0[1][i+1]-mass_v0[1][i])/50)*(t_ux_gt-mass_v0[0][i]);} /*удельный объём воздуха*/

}

O2_ux_gt=(0.21*v0_ks*(alfa_ux_gt-1))/(v0_ks*alfa_ux_gt+1); /*объёмное содержание кислорода в окислителе после ГТ*/

n_ptu=n_t_ptu*n_oi*n_m*n_g*n_pg*n_tr; /*КПД ПТУ*/

G1_pg=0.992*G_k_ohl*0.5; /*расход окислителя, идущего в топку НПГ от ГТ*/

G_dv=G1_pg*(O2_t_pg-O2_ux_gt*100)/(21-O2_t_pg); /*расход добавочного воздуха*/

G_pg=G1_pg+G_dv; /*суммарный расход окислителя, идущего на горение в топку НПГ*/

q_dv=G_dv/(G_k_ohl*0.5*0.992); /*относительная добавка воздуха*/

i_ne=hpt(P0, t0); /*энтальпия перегретого пара*/

s_para=s(P0,i_ne); /*энтропия пара при идеальном процессе расширения*/

i_k=hps(Pk,s_para); /*энтальпия в конденсаторе*/

H_pt=i_ne-i_k;/*теплоперепад ПТ*/

/*цикл уточнения коэффициента регенерации*/

do

{

D_ne=(N_e_pt/(H_pt*n_oi*n_m*n_g))*K_p; /*расход острого пара на ПТ*/

q_pv=t_pv*Cp; /*энтальпия питательной воды*/

D_np=0.015*D_ne; /*расход продувочной воды*/

P_b=1.1*P0; /*давление в барабане НПГ*/

i_np=h1(ts(P_b)); /*энтальпия продувочной воды*/

D_ezh=0.005*D_ne; /*расход пара на эжектора*/

D_upl=0.01*D_ne; /*расхлд пара на уплотнения ПТ*/

D_ut=0.015*D_ne; /*расход пара на утечки*/

D_sn=0.024*D_ne;/*расход пара на собственные нужды*/

D_ne0=D_ne+D_ezh+D_upl+D_sn; /*расход острого пара с НПГ*/

Q_pg=D_ne0*(i_ne-q_pv)+D_np*(i_np-q_pv); /*теплота, использованная в НПГ*/

i1_np=h1(ts(P_d));/*энтальпия продувочной воды за РНП*/

i_c1=h11(ts(P_d)); /*энтальпия выпара из РНП*/

D_c1=D_np*(i_np-i1_np)/(i_c1-i1_np); /*расход выпара из РНП*/

D_np1=D_np-D_c1; /*расход продувочной воды, сбрасываемой в технологическую канализацию*/

D_xov=D_np1+D_ut+D_sn; /*расход ХОВ*/

t1_xov=t_xov+D_np1*((i1_np/Cp)-60)/D_xov;/*температура ХОВ на входе в деаэратор*/

for (i=0; i<21; i++)

{

if (t_d>mass_I_sv[0][i])

{I4_sv=mass_I_sv[1][i]+((mass_I_sv[1][i+1]-mass_I_sv[1][i])/50)*(t_d-mass_I_sv[0][i]);}

/*энтальпия сухого воздуха при температуре за ГТ*/

}

for (i=0; i<21; i++)

{

if (t_d>mass_I_sg[0][i])

{I4_sg=mass_I_sg[1][i]+((mass_I_sg[1][i+1]-mass_I_sg[1][i])/50)*(t_d-mass_I_sg[0][i]);}

/*энтальпия продуктов сгорания при температуре за ГТ*/

}

for (i=0; i<21; i++)

{

if (t11_ek>mass_I_sg[0][i])

{I_ek_sg=mass_I_sg[1][i]+((mass_I_sg[1][i+1]-mass_I_sg[1][i])/50)*(t11_ek-mass_I_sg[0][i]);}

/*энтальпия газов ГТ при температуре за экономайзером*/

}

q_tl=c_tl*t_tl; /*теплота, внесённая с топливом в НПГ*/

for (i=0; i<21; i++)

{

if (t11_ek>mass_I_sv[0][i])

{I_ek_sv=mass_I_sv[1][i]+((mass_I_sv[1][i+1]-mass_I_sv[1][i])/50)*(t11_ek-mass_I_sv[0][i]);}

/*энтальпия воздуха при температуре за экономайзером*/

}

I_ek_sg_pg=I_ek_sg; /* энтальпия газов НПГ при температуре за экономайзером */

for (i=0; i<21; i++)

{

if (t_dv>mass_I_sv[0][i])

{h_dv=mass_I_sv[1][i]+((mass_I_sv[1][i+1]-mass_I_sv[1][i])/50)*(t_dv-mass_I_sv[0][i]);}

/* энтальпия добавочного воздуха после ВЗП */

}

for (i=0; i<21; i++)

{

if (t_dv0>mass_I_sv[0][i])

{h_dv0=mass_I_sv[1][i]+((mass_I_sv[1][i+1]-mass_I_sv[1][i])/50)*(t_dv0-mass_I_sv[0][i]);}

/* энтальпия добавочного воздуха перед ВЗП */

}

B_r_pg_nt=((D_ne*(i_ne-q_pv)+G_dv*(I_ek_sv-h_dv)+G_t_ohl_*(I_ek_sg-I4_sg))/((Q_n_r_gas-I_ek_sg)*n_pg))*3.6;/*предварительный расход топлива*/

alfa_yx_pg=(G_k_ohl*0.992*0.5-B*L_0*n_ks*0.5)/(B_r_pg_nt/3.6*L_0);

/*коэффициент избытка воздуха в уходящих газах после экономайзера НПГ без учёта добавочного воздуха*/

alfa_yx_pg_dv=(G_k_ohl*(0.992+q_dv)*0.5-B*L_0*n_ks*0.5)/(B_r_pg_nt/3.6*L_0);

/* коэффициент избытка воздуха в уходящих газах после экономайзера НПГ с учётом добавочного воздуха */

Bp1_pg=((Q_pg-B*n_ks*(I4_sg+(alfa_ux_gt-1)*I4_sv-I_ek_sg)))*3.6/((1-0.01*(q5))*(Q_n_r_gas+q_tl)+(alfa_yx_pg_dv-alfa_yx_pg)*L_0*c_dv*t_dv-(I_ek_sg+I_ek_sv*(alfa_yx_pg_dv-1)));

/*уточнённый расход топлива в топку НПГ*/

/*Тепловой баланс ГВПВД*/

t1_vzp=t11_ek-100; /*температура газов за ГВПВД, ?C*/

for (i=0; i<21; i++)

{

if (t1_vzp>mass_I_sv[0][i])

{I1_vzp_sv=mass_I_sv[1][i]+((mass_I_sv[1][i+1]-mass_I_sv[1][i])/50)*(t1_vzp-mass_I_sv[0][i]);}

/*энтальпия воздуха при температуре за ГВПВД*/

}

for (i=0; i<21; i++)

{

if (t1_vzp>mass_I_sg[0][i])

{I1_vzp_sg=mass_I_sg[1][i]+((mass_I_sg[1][i+1]-mass_I_sg[1][i])/50)*(t1_vzp-mass_I_sg[0][i]);}

/*энтальпия газов ГТ при температуре за ГВПВД*/

}

h_d0=h1(ts(P_d)); /*энтальпия воды за деаэратором*/

v_pv=v1(ts(P_d)); /*удельный объём питательной воды*/

h_d=h_d0+dP_pn*pow(10,5)*v_pv/(n_pn*1000); /*энтальпия питательной воды за ПН*/

D_pv=D_ne0+D_np; /*расходпитательной воды*/

Q_gvpvd_g=((G_t_ohl_+Bp1_pg/3.6)*(I_ek_sg-I1_vzp_sg)+G_dv*(I_ek_sv-I1_vzp_sv))*n_to;

/*тепловой баланс ГВПВД по газу*/

Q_gvpvd_v=D_pv*(q_pv-h_d); /*тепловой баланс ГВПВД по воде*/

if (Q_gvpvd_g>Q_gvpvd_v)

{do

{for (i=0; i<21; i++)

{

if (t1_vzp>mass_I_sv[0][i])

{I1_vzp_sv=mass_I_sv[1][i]+((mass_I_sv[1][i+1]-mass_I_sv[1][i])/50)*(t1_vzp-mass_I_sv[0][i]);}

}

for (i=0; i<21; i++)

{

if (t1_vzp>mass_I_sg[0][i])

{I1_vzp_sg=mass_I_sg[1][i]+((mass_I_sg[1][i+1]-mass_I_sg[1][i])/50)*(t1_vzp-mass_I_sg[0][i]);}

}

h_d0=h1(ts(P_d));

v_pv=v1(ts(P_d));

h_d=h_d0+dP_pn*pow(10,5)*v_pv/(n_pn*1000);

D_pv=D_ne0+D_np;

Q_gvpvd_g=((G_t_ohl_+Bp1_pg/3.6)*(I_ek_sg-I1_vzp_sg)+G_dv*(I_ek_sv-I1_vzp_sv))*n_to;

Q_gvpvd_v=D_pv*(q_pv-h_d);

t1_vzp=t1_vzp+0.1;

}

while ((fabs(Q_gvpvd_g-Q_gvpvd_v)/Q_gvpvd_g)*100>=0.5);}

else

{do

{for (i=0; i<21; i++)

{

if (t1_vzp>mass_I_sv[0][i])

{I1_vzp_sv=mass_I_sv[1][i]+((mass_I_sv[1][i+1]-mass_I_sv[1][i])/50)*(t1_vzp-mass_I_sv[0][i]);}

}

for (i=0; i<21; i++)

{

if (t1_vzp>mass_I_sg[0][i])

{I1_vzp_sg=mass_I_sg[1][i]+((mass_I_sg[1][i+1]-mass_I_sg[1][i])/50)*(t1_vzp-mass_I_sg[0][i]);}

}

h_d0=h1(ts(P_d));

v_pv=v1(ts(P_d));

h_d=h_d0+dP_pn*pow(10,5)*v_pv/(n_pn*1000);

D_pv=D_ne0+D_np;

Q_gvpvd_g=((G_t_ohl_+Bp1_pg/3.6)*(I_ek_sg-I1_vzp_sg)+G_dv*(I_ek_sv-I1_vzp_sv))*n_to;

Q_gvpvd_v=D_pv*(q_pv-h_d);

t1_vzp=t1_vzp-0.1;}

while ((fabs(Q_gvpvd_g-Q_gvpvd_v)/Q_gvpvd_g)*100>=0.5);}

/*Тепловой баланс ВЗП*/

t11_vzp=t1_vzp-50; /*температура газов за ВЗП*/

for (i=0; i<21; i++)

{

if (t11_vzp>mass_I_sv[0][i])

{I11_vzp_sv=mass_I_sv[1][i]+((mass_I_sv[1][i+1]-mass_I_sv[1][i])/50)*(t11_vzp-mass_I_sv[0][i]);}

/*энтальпия воздуха при температуре за ВЗП*/

}

for (i=0; i<21; i++)

{

if (t11_vzp>mass_I_sg[0][i])

{I11_vzp_sg=mass_I_sg[1][i]+((mass_I_sg[1][i+1]-mass_I_sg[1][i])/50)*(t11_vzp-mass_I_sg[0][i]);}

/* энтальпия газов ГТ при температуре за ВЗП */

}

Q_vzp_g=((G_t_ohl_+Bp1_pg/3.6)*(I1_vzp_sg-I11_vzp_sg)+G_dv*(I1_vzp_sv-I11_vzp_sv))*n_to;

/*тепловой баланс ВЗП по газу*/

Q_vzp_vozd=G_dv*(h_dv-h_dv0); /*тепловой баланс ВЗП по воздуху*/

if (Q_vzp_g>Q_vzp_vozd)

{do

{for (i=0; i<21; i++)

{

if (t11_vzp>mass_I_sv[0][i])

{I11_vzp_sv=mass_I_sv[1][i]+((mass_I_sv[1][i+1]-mass_I_sv[1][i])/50)*(t11_vzp-mass_I_sv[0][i]);}

}

for (i=0; i<21; i++)

{

if (t11_vzp>mass_I_sg[0][i])

{I11_vzp_sg=mass_I_sg[1][i]+((mass_I_sg[1][i+1]-mass_I_sg[1][i])/50)*(t11_vzp-mass_I_sg[0][i]);}

}

Q_vzp_g=((G_t_ohl_+Bp1_pg/3.6)*(I1_vzp_sg-I11_vzp_sg)+G_dv*(I1_vzp_sv-I11_vzp_sv))*n_to;

Q_vzp_vozd=G_dv*(h_dv-h_dv0);

t11_vzp=t11_vzp+0.1;}

while ((fabs(Q_vzp_g-Q_vzp_vozd)/Q_vzp_g)*100>=0.5);}

else

{do

{for (i=0; i<21; i++)

{

if (t11_vzp>mass_I_sv[0][i])

{I11_vzp_sv=mass_I_sv[1][i]+((mass_I_sv[1][i+1]-mass_I_sv[1][i])/50)*(t11_vzp-mass_I_sv[0][i]);}

}

for (i=0; i<21; i++)

{

if (t11_vzp>mass_I_sg[0][i])

{I11_vzp_sg=mass_I_sg[1][i]+((mass_I_sg[1][i+1]-mass_I_sg[1][i])/50)*(t11_vzp-mass_I_sg[0][i]);}

}

Q_vzp_g=((G_t_ohl_+Bp1_pg/3.6)*(I1_vzp_sg-I11_vzp_sg)+G_dv*(I1_vzp_sv-I11_vzp_sv))*n_to;

Q_vzp_vozd=G_dv*(h_dv-h_dv0);

t11_vzp=t11_vzp-0.1;}

while ((fabs(Q_vzp_g-Q_vzp_vozd)/Q_vzp_g)*100>=0.5);}

/*Деаэратор*/

i0_otb=hps(P_otb, s_para); /*теоретическая энтальпия отборного пара*/

i_otb=i_ne-(i_ne-i0_otb)*n_oi; /*действительнаяая энтальпия отборного пара */

D_d=((D_pv+D_ut)*(h_d0-t1_ok*Cp)-D_c1*(i_c1-t1_ok*Cp)-D_xov*(t1_xov*Cp-t1_ok*Cp)-(D_ezh+D_upl)*(h_d0-t1_ok*Cp))/(i_otb-t1_ok*Cp); /*расход пара на деаэратор*/

D_ok=D_pv+D_ut-D_d-D_c1-D_xov-D_ezh-D_upl; /*расход основного конденсата*/

/*Точка смешения (расчёт рециркуляции для обеспечениянеобходимой температуры основного конденсата перед ГВПНД)*/

h_ok0=h1(ts(Pk)); /*энтальпия воды за конденсатором*/

D_rec=D_ok*(h_ok0-t_ok*Cp)/(t_ok*Cp-t1_ok*Cp); /*расход рециркуляции*/

D1_ok=D_ok+D_rec; /*расход основного конденсата через ГВПНД*/

/*Тепловой баланс ГВПНД*/

t_ux=t11_vzp-90; /*температура газов за ГВПНД*/

for (i=0; i<21; i++)

{

if (t_ux>mass_I_sv[0][i])

{I_ux_sv=mass_I_sv[1][i]+((mass_I_sv[1][i+1]-mass_I_sv[1][i])/50)*(t_ux-mass_I_sv[0][i]);}

/*энтальпия воздуха при температуре за ГВПНД*/

}

for (i=0; i<21; i++)

{

if (t_ux>mass_I_sg[0][i])

{I_ux_sg=mass_I_sg[1][i]+((mass_I_sg[1][i+1]-mass_I_sg[1][i])/50)*(t_ux-mass_I_sg[0][i]);}

/*энтальпия газов ГТ при температуре за ГВПНД*/

}

Q_gvpnd_g=((G_t_ohl_+Bp1_pg/3.6)*(I11_vzp_sg-I_ux_sg)+G_dv*(I11_vzp_sv-I_ux_sv))*n_to;

/*тепловой баланс ГВПНД по газу*/

Q_gvpnd_v=D1_ok*(t1_ok*Cp-t_ok*Cp); /*тепловой баланс ГВПНД по воде*/

if (Q_gvpnd_g>Q_gvpnd_v)

{do

{for (i=0; i<21; i++)

{

if (t_ux>mass_I_sv[0][i])

{I_ux_sv=mass_I_sv[1][i]+((mass_I_sv[1][i+1]-mass_I_sv[1][i])/50)*(t_ux-mass_I_sv[0][i]);}

}

for (i=0; i<21; i++)

{

if (t_ux>mass_I_sg[0][i])

{I_ux_sg=mass_I_sg[1][i]+((mass_I_sg[1][i+1]-mass_I_sg[1][i])/50)*(t_ux-mass_I_sg[0][i]);}

}

Q_gvpnd_g=((G_t_ohl_+Bp1_pg/3.6)*(I11_vzp_sg-I_ux_sg)+G_dv*(I11_vzp_sv-I_ux_sv))*n_to;

Q_gvpnd_v=D1_ok*(t1_ok*Cp-t_ok*Cp);

t_ux=t_ux+0.1;}

while ((fabs(Q_gvpnd_g-Q_gvpnd_v)/Q_gvpnd_g)*100>=0.5);}

else

{do

{for (i=0; i<21; i++)

{

if (t_ux>mass_I_sv[0][i])

{I_ux_sv=mass_I_sv[1][i]+((mass_I_sv[1][i+1]-mass_I_sv[1][i])/50)*(t_ux-mass_I_sv[0][i]);}

}

for (i=0; i<21; i++)

{

if (t_ux>mass_I_sg[0][i])

{I_ux_sg=mass_I_sg[1][i]+((mass_I_sg[1][i+1]-mass_I_sg[1][i])/50)*(t_ux-mass_I_sg[0][i]);}

}

Q_gvpnd_g=((G_t_ohl_+Bp1_pg/3.6)*(I11_vzp_sg-I_ux_sg)+G_dv*(I11_vzp_sv-I_ux_sv))*n_to;

Q_gvpnd_v=D1_ok*(t1_ok*Cp-t_ok*Cp);

t_ux=t_ux-0.1;}

while ((fabs(Q_gvpnd_g-Q_gvpnd_v)/Q_gvpnd_g)*100>=0.5);}

/*Проверка по балансу мощности паровой турбины*/

Wep=(D_d*(i_ne-i_otb)+D_ok*(H_pt*n_oi))*n_em;/*мощность паровой турбины*/

d=fabs((N_e_pt-Wep)/N_e_pt)*100; /*погрешность расчёта*/

K_p=K_p*(N_e_pt/Wep); /*уточнение коэффициента регенерации*/

}

while(d>=0.5);

/*Расчёт технико-экономических показателей*/

N_sn_pt=2*0.038*Wep; /*мощность механизмов собственных нужд паровой турбины*/

N_sn_gt=0.01*N_el_gt; /* мощность механизмов собственных нужд газовой турбины */

N_sn_pgu=N_sn_pt+N_sn_gt; /* мощность механизмов собственных нужд ПГУ*/

N_pgu_br=2*Wep+N_el_gt; /*мощность ПГУ брутто*/

N_pgu_nt=N_pgu_br-N_sn_pgu; /*мощность ПГУ нетто*/

n_pgu_br=(N_pgu_br*100)/(Q_n_r_gas*(B+2*Bp1_pg/3.6)); /*КПД ПГУ брутто*/

n_pgu_nt=(N_pgu_nt*100)/(Q_n_r_gas*(B+2*Bp1_pg/3.6)); /*КПД ПГУ нетто*/

Размещено на Allbest.ru