Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

Федеральное государственное образовательное учреждение

высшего профессионального образования

кафедра

КУРСОВОЙ ПРОЕКТ

Расчет принципиальной тепловой схемы энергоблока

Красноярск 2011

Исходные данные

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

Дано:

W_э=135 МВт; Р₀=14МПа; t₀=570 ⁰С;P₁=2,5 МПа;P₂=1,7МПа;Р₃=1,1 МПа;Р₄=0,12 МПа;Р_д=0,6;Р_к=0,003МПа;Р_дв=0,12 МПа;D_п=16 кг/с; дD_п=0,6; Р_ок=1,2; Q_m=15МВт; t_хов=32^оС; з_0i=0,89; з_эм=0,99; з_т=0,98; з_н=0,7; з_к=0,91, k_рег=1,2.

Рис 1. принципиальная тепловая схема ТЭС

1. Расчет процесса расширения пара в турбине

Определяем энтальпию острого пара перед турбиной:

1) Теоретический теплоперепад в первом отсеке:

действительный теплоперепад в первом отсеке:

Энтальпия пара в первом отборе:

2) Теоретический теплоперепад во втором отсеке:

Действительный теплоперепад во втором отсеке:

Энтальпия пара во втором отборе:

3) Теоретический теплоперепад в третьем отсеке:

Действительный теплоперепад в третьем отсеке:

Энтальпия пара в третьем отборе:

4) Теоретический теплоперепад в четвертом отсеке:

Действительный теплоперепад в четвертом отсеке:

Энтальпия пара в четвертом отборе:

Теоретический теплоперепад в пятом отсеке:

Действительный теплоперепад в пятом отсеке:

Энтальпия пара в пятом отсеке:

Рис.2 Процесс расширения пара в турбине

2. Расчет расхода пара на турбину

где - коэффициент регенерации,

- электрическая мощность,

- электромеханический КПД,

- коэффициент недовыработки паром теплофикационного отбора,

- расход пара на отопление,

- коэффициент недовыработки паром отопительного отбора,

- расход пара на производство.

где - энтальпия в теплофикационном отборе.

Где - тепловая энергия, - КПД теплообменника, - энтальпия конденсата пара теплофикационного отбора.

3. Определение расхода питательной воды на котельный агрегат

Для прямоточного котла:

где - внутренние потери пара и конденсата,

- расход пара на собственные нужды станции.

/с

Расход химически очищенной воды:

Энтальпия химически очищенной воды:

4. Расчет регенеративной схемы

Расчет ПВД-1.

Расчетная схема подогревателя высокого давления 1.

1) Расчет ПВД-2.

2) Расчет питательного насоса

Расчетная схема питательного насоса

3) Расчет деаэратора

Расчетная схема деаэрационной установки

4) Расчет дренажного насоса

Расчетная схема дренажного насоса

5) Расчет деаэратора добавочной воды

Расчетная схема деаэрационной установки добавочной воды

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

6) Расчетная схема подогревателя низкого давления , насоса и двух точек смешения

Расчетная схема подогревателя низкого давления, насоса и две точки смешения.



5. Определение расчетной нагрузки

По найденной расчетной мощности найдем погрешность расчета.

Расчетная погрешность больше 0,5%поэтому необходимо выполнить перерасчет тепловой схемы заново, начиная с момента определения расхода острого пара на турбину, при этом новое значение коэффициента регенерации находится по следующей зависимости:

6. Расчет технико-экономических показателей тепловой схемы

Полный расход топлива:

где - расход острого пара на турбину ;

- соответственно энтальпии острого пара и питательной воды;

- КПД котельного агрегата;

- теплота сгорания условного топлива.

Расход топлива на выработку тепловой энергии.

где - расход тепла с турбины на теплофикационные нужды.

- полезное тепло, полученное в котельной установке.

КПД брутто по выработке тепловой энергии.

КПД брутто по выработке электроэнергии.

Удельный расход топлива на выработку тепловой энергии:

Удельный расход топлива на выработку электроэнергии:



7. Текст программы

#include <stdio.h>

#include <conio.h>

#include <math.h>

#include <eheat.h>

#include <emath.h>

#include <stdlib.h>

void main ()

{

clrscr();

float P0=140,T0=570,P1=25,P2=17,P3=11,P4=1.2,Pd=6,Pk=0.03,Pdv=1.2,W=135000,

Dn=16,qDn=0.6,Tvk=90,Qt=15000,Txov=32,noi=0.89,nYm=0.97,nka=0.91,nt=0.98;

float kreg=1.2,Tdt,Dm,Dyt,Dcn,Gpv,Dxov,D0,D3,D4;

float Ho1,Ho2,Ho3,Ho4,Ho5,hH1,hH2,hH3,hH4,hH5,H0,H1,H2,H3,H4,H5,Hd,ym,yn;

float Sigma,Z=0.5,Qnr=29330;

float O=4,PN,nm=0.99,tpv,Hnv,Hd1,t1pv,H1nv,td1,Hvk;

float B,nbr,bkec,Qk,Qmm,Btetec,Beetec,nbrtetec,nbreetec,beetec,btetec;

float nd=0.7,td,Vd,H111nv,Ht1,H11nv,Hd2,td2,PN1;

FILE*pfi;

pfi=fopen("C:\\rezult.txt","w");

if (pfi==NULL)

{printf ("err");

exit(1);}

H0=hpt(P0,T0);

Ho1=h0(P0,H0,P1);

hH1=Ho1*noi;

H1=H0-hH1;

Ho2=h0(P1,H1,P2);

hH2=Ho2*noi;

H2=H1-hH2;

Ho3=h0(P2,H2,P3);

hH3=Ho3*noi;

H3=H2-hH3;

Ho4=h0(P3,H3,P4);

hH4=Ho4*noi;

H4=H3-hH4;

Ho5=h0(P4,H4,Pk);

hH5=Ho5*noi;

H5=H4-hH5;

Hd=h1(ts(P4));

Dm=Qt/((H4-Hd)*nYm);

ym=(H4-H5)/(H0-H5);

yn=(H3-H5)/(H0-H5);

do

{

D0=kreg*(W/(hH1+hH2+hH3+hH4+hH5)*nYm)+yn*Dn+yn*Dm;

Dyt=0.015*D0;

Dcn=0.025*D0;

Gpv=D0+Dyt+Dcn;

Dxov=Dyt+Dcn;

PN1=1.4*Pd;

PN=P0*1.4;

tpv=ts(P1)-O;

Hnv=hpt(PN,tpv);

td1=ts(P1)+10;

Hd1=h1(td1);

t1pv=ts(P2)-O;

H1nv=hpt(PN,t1pv);

float D1;

D1=(Gpv*(Hnv-H1nv))/((H1-Hd1)*nt);

td=ts(Pd);

Vd=v1(td);

H111nv=h1(td);

Ht1=((PN1-Pd)*Vd/nd)*100;

H11nv=H111nv+Ht1;

td2=ts(P2);

Hd2=h1(td2);

float D2;

D2=(Gpv*(H1nv-H11nv)-D1*(Hd1-Hd2))*nYm/(H2-Hd2)*nYm;

float Pok=13,K=5,t1ok,Gxov,H1ok,PH1,Hvc,tvc,Vdv,Hdv,Ht2,Dxov,Hxov,Cp=4.19,Hvk;

t1ok=ts(P4)-K;

H1ok=hpt(Pok,t1ok);

PH1=1.4*Pd;

tvc=ts(Pdv);

Hvc=h1(tvc);

Vdv=v1(tvc);

Ht2=((PH1-Pdv)*Vdv*100)/nd;

Hdv=Hvc+Ht2;

Dxov=Dyt+Dcn;

Hxov=Txov*Cp;

Hvk=Tvk*Cp;

Gxov=Dyt+Dcn;

float det,det1,det2,det3,det4;

float a[4][4]={{1,1,1,0},

{H3*nYm,H1ok,Hdv,0},

{0,0,-1,1},

{0,0,-Hvc,H4}};

det=MatDeter(&a[0][0],4);

float b[4][4]={{Gpv-D1-D2,1,1,0},

{Gpv*H111nv-(D1+D2)*Hd2,H1ok,Hdv,0},

{-Gxov-Dn*qDn,0,-1,1},

{-Gxov*Hxov-(Dn*qDn)*Hvc,0,-Hvc,H4}};

det1=MatDeter(&b[0][0],4);

float c[4][4]={{1,Gpv-D1-D2,1,0},

{H3*nYm,Gpv*H111nv-(D1+D2)*Hd2,Hdv,0},

{0,-Gxov-Dn*qDn,-1,1},

{0,-Gxov*Hxov-(Dn*qDn)*Hvc,-Hvc,H4}};

det2=MatDeter(&c[0][0],4);

float d[4][4]={{1,1,Gpv-D1-D2,0},

{H3*nYm,H1ok,-Gxov-Dn*qDn,0},

{0,0,-Gxov-Dn*qDn,1},

{0,0,-Gxov*Hxov-(Dn*qDn)*Hvc,H4}};

det3=MatDeter(&d[0][0],4);

float e[4][4]={{1,1,1,Gpv-D1-D2},

{H3*nYm,H1ok,Hdv,-Gxov-Dn*qDn},

{0,0,-1,-Gxov-Dn*qDn},

{0,0,-Hvc,-Gxov*Hxov-(Dn*qDn)*Hvc}};

det4=MatDeter(&e[0][0],4);

float D113,Gok,Gdv,Dm1;

D113=det1/det;

Gok=det2/det;

Gdv=det3/det;

Dm1=det4/det;

float Vcm=0.0011,td3,Hd3,Ht3,Hk,tk;

td3=ts(P4);

tk=ts(Pk);

Hk=h1(tk);

float dett,det11,det22,det33,det44;

float z[4][4]={{1,1,0,0,},

{0,-Hd3,Gok,-Hk},

{-1,0,0,1},

{0,H4-Hd3,Gok,0}};

dett=MatDeter(&z[0][0],4);

float x[4][4]={{Gok,1,0,0,},

{Dm*Hd,-Hd3,Gok,-Hk},

{-Dm,0,0,1},

{Gok*H1ok+1.949*Gok,H4-Hd3,Gok,0}};

det11=MatDeter(&x[0][0],4);

float v[4][4]={{1,Gok,0,0,},

{0,Dm*Hd,Gok,-Hk},

{-1,-Dm,0,1},

{0,Gok*H1ok+1.949*Gok,Gok,0}};

det22=MatDeter(&v[0][0],4);

float f[4][4]={{1,1,Gok,0,},

{0,-Hd3,Dm*Hd,-Hk},

{-1,0,-Dm,1},

{0,H4-Hd3,Gok*H1ok+1.949*Gok,0}};

det33=MatDeter(&f[0][0],4);

float l[4][4]={{1,1,0,Gok},

{0,-Hd3,Gok,Dm*Hd},

{-1,0,0,-Dm},

{0,H4-Hd3,Gok,Gok*H1ok+1.949*Gok}};

det44=MatDeter(&l[0][0],4);

float Gcm,Dm2,Hcm,Dk;

Gcm=det11/dett;

Dm2=det22/dett;

Hcm=det33/dett;

Dk=det44/dett;

float Wep,D11,D22,D33,D44,D55;

D11=D0;

D22=D0-D1;

D33=D0-D1-D2;

D44=D0-D1-D2-Dn-D113;

D55=D44-Dm-Dm1-Dm2;

Wep=(D11*hH1+D22*hH2+D33*hH3+D44*hH4+D55*hH5)*nYm;

Sigma=fabs((Wep-W)/W)*100;

kreg=kreg*W/Wep;

D3=D113+(Dn*qDn);

D4=Dm+Dm1+Dm2;

int i=0;

i++;

if (Sigma<Z||i==1)

{

fprintf (pfi,"\n\n h0=%.2f \th1=%.2f \th2=%.2f\n h3=%.2f \th4=%.2f \thk=%.2f",H0,H1,H2,H3,H4,H5);

fprintf (pfi,"\n\n D0=%.2f \tD1=%.2f \tD2=%.2f ",D0,D1,D2);

fprintf (pfi,"\n\n D3=%.2f \tD4=%.2f \tDk=%.2f",D3,D4,Dk);

fprintf (pfi,"\n\n Gpv=%.2f\tGok=%.2f\tGdv=%.2f",Gpv,Gok,Gdv);

fprintf (pfi,"\n\n kreg=%.5f\tWep=%.2f\tSigma=%.2f",kreg,Wep,Sigma);

}

}

while (Sigma >Z);

B=D0*(H0-Hnv)/(Qnr*nka);

nbr=W/(B*Qnr);

bkec=0.123/nbr;

Qk=D0*(H0-Hnv);

Qmm=(Dn*(H3-Hvk))+(Dm*(H4-Hd));

Btetec=B*Qmm/Qk;

Beetec=B-Btetec;

nbreetec=W/(Beetec*Qnr);

nbrtetec=Qmm/(Btetec*Qnr);

beetec=0.123/nbreetec;

btetec=0.0342/nbrtetec;

fprintf(pfi,"\n\n\t\t Tehnico-economicheskii raschet");

fprintf(pfi,"\n Rashod yslovnogo topliva\t\t\t\t\t\t= %.2f",B);

fprintf(pfi,"\n Rashod topliva na virabotky teplovoi energii\t\t\t\t=%.2f",Btetec);

fprintf(pfi,"\n Rashod topliva na virabotky electricheskoi energii\t\t\t=%.2f",Beetec);

fprintf(pfi,"\n KPD Brutto po virabotke elektroenergii\t\t\t\t\t=%.2f",nbreetec);

fprintf(pfi,"\n KPD Brutto po virabotke tepla\t\t\t\t\t\t=%.2f",nbrtetec);

fprintf(pfi,"\n Ydelnii rashod yslovnogo topliva (kg/kWt*h) na virabotky 1 kWt*h\t=%.2f",beetec);

fprintf(pfi,"\n Ydelnii rashod yslovnogo topliva (kg/MDj) na virabotky 1 MDj\t\t=%.4f",btetec);

fclose(pfi);

getch();

}

h0=3511.49 h1=3068.79 h2=2982.82

h3=2893.60 h4=2535.15 hk=2101.75

D0=124.60 D1=5.60 D2=10.67

D3=20.02 D4=18.60 Dk=69.72

Gpv=129.58 Gok=87.59 Gdv=15.31

kreg=1.26657 Wep=127904.27 Sigma=5.26

h0=3511.49 h1=3068.79 h2=2982.82

h3=2893.60 h4=2535.15 hk=2101.75

D0=130.78 D1=5.88 D2=11.20

D3=20.54 D4=19.27 Dk=73.90

Gpv=136.01 Gok=92.41 Gdv=15.59

kreg=1.26690 Wep=134965.27 Sigma=0.03

Tehnico-economicheskii raschet

Rashod yslovnogo topliva = 12.56

Rashod topliva na virabotky teplovoi energii =1.71

Rashod topliva na virabotky electricheskoi energii =10.85

KPD Brutto po virabotke elektroenergii =0.42

KPD Brutto po virabotke tepla =0.91

Ydelnii rashod yslovnogo topliva (kg/kWt*h) na virabotky 1 kWt*h =0.29

Ydelnii rashod yslovnogo topliva (kg/MDj) na virabotky 1 MDj =0.037

Список использованных источников

тепловая схема энергоблок

1. Математическое моделирование теплоэнергетических задач на ЭВМ: Методические указания по лабораторным работам и расчетно графическому заданию для студентов специальностей 140101-«Тепловые Электрические Станции»/ Составили Е.А. Бойко, Д.Г. Дидичин, П.В. Шишмарев.Красноярск: ИПЦ КГТУ, 2002.120с.

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