Добавил:
ivanov666
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Оптимизация источников теплоты и систем теплоснабжения. Учебное пособие для обучающихся направления подготовки 08.06.01 «Техника и технология строител
.pdf
61
locate 11,73:?" И "
locate 12,73:?" Т "
locate 13,73:?" Е "
color 3,1
locate 18,5
color 18,8
locate 18,16
locate 18,32
ppmin=7899999*90000000
proc=1
1119
if n=1 then 100 else 200
100
y=(3-1.5)/400*(500-G): qu(1)=3-y: qu(2)=qu(1)
y=(1500-150)/400*(500-g): Q(1)=(1500-y)*1000: Q(2)=Q(1)
y=(1.0-.3)/400*(500-g): Mu(1)=1.0-y: mu(2)=mu(1)
y=(1.0-0.35)/400*(500-g): k(1)=1.0-y
nnn=1: fff=2: k(2)=k(1)
Y=(60-6)/400*(500-g): M1(1)=60-y
y=(700-70)/400*(500-g): M1(2)=700-y
y=(3000-2000)/400*(500-g): MuTB(1)=3000-y: MuTB(2)=MuTB(1)
y=(10*10^6-21*10^5)/400*(500-g): V(1)=(10*10^6)-y
y=(63*10^5-10*10^5)/400*(500-g): V(2)=(63*10^5)-y
y=(5-1)/400*(500-g): zXBO(1)=5-y: zXBO(2)=zxbo(1)
y=(0.8-0.6)/400*(500-g): Pk(1)=0.8-y: Pk(2)=Pk(1)
QTEPH(1)=7.5: QTEPK(1)=75
QTEPH(2)=3.8: QTEPK(2)=45: hag=2
y=(10^6-6*10^5)/400*(500-g): dp(1)=10^6-y: dp(2)=dp(1)
goto 1000
200
if n=2 then 300 else 400
300
y=(3-1.5)/80*(100-G): qu(3)=3-y: qu(4)=qu(3)
y=(300-30)/80*(100-g): Q(3)=(300-y)*1000: Q(4)=q(3)
y=(1.0-.3)/80*(100-g): Mu(3)=1.0-y: mu(4)=mu(3)
62
K(3)=2-(2-1.5)/80*(100-g): k(4)=k(3): nnn=3
Y=(40-4)/80*(100-g): M1(3)=40-y: fff=4
y=(150-15)/80*(100-g): M1(4)=150-y
y=(3000-2000)/80*(100-g): MuTB(3)=3000-y: MuTB(4)=MuTB(3)
y=(21*10^5-4*10^5)/80*(100-g): V(3)=(21*10^5)-y
y=(10*10^5-2*10^5)/80*(100-g): V(4)=(10*10^5)-y
y=(4-1)/80*(100-g): zXBO(3)=4-y: zXBO(4)=zxbo(3)
y=(0.6-0.4)/80*(100-g): Pk(3)=0.6-y
y=(0.5-0.3)/80*(100-g): Pk(4)=0.5-y
QTEPH(3)=4.5: QTEPK(3)=45
QTEPH(4)=2.5: QTEPK(4)=25: hag=1.5
y=(10^6-6*10^5)/80*(100-g): dP(3)=(10^6)-y: dP(4)=dp(3)
goto 1000
400
y=(3-1.5)/48*(50-G): qu(5)=3-y: qu(6)=qu(5)
y=(150-3)/48*(50-g): Q(5)=(150-y)*1000: Q(6)=q(5)
y=(1.0-.3)/48*(50-g): Mu(5)=1.0-y: mu(6)=mu(5)
K(5)=2-(2-1.5)/48*(50-g): k(6)=k(5)
nnn=5: Y=(70-2)/48*(50-g): M1(5)=70-y: fff=6
y=(70-2)/48*(50-g): M1(6)=70-y
y=(3000-2000)/48*(50-g): MuTB(5)=3000-y: MuTB(6)=MuTB(5)
y=(10*10^5-4.2*10^4)/48*(50-g): V(5)=(10*10^5)-y
y=(10*10^5-4.2*10^4)/48*(50-g): V(6)=(10*10^5)-y
y=(3-1)/48*(50-g): zXBO(5)=3-y: zXBO(6)=zxbo(5)
y=(0.5-0.3)/48*(50-g): Pk(5)=0.5-y
y=(0.4-0.3)/48*(50-g): Pk(6)=0.4-y
QTEPH(5)=0.15: QTEPK(5)=7.5
QTEPH(6)=0.15: QTEPK(6)=7.5: hag=0.25
y=(10^6-6*10^5)/48*(50-g): dP(5)=10^6-y: dp(6)=dp(5)
1000 for i=nnn to fff step 1
if g>300 then QTEPH(i)=QTEPH(I)*2
for t2=40 to 80 step 2
if proc<100 then proc=3+proc
if n=1 then 1011 else 1012
1011 for t1=105 to 200 step 2
63
goto 1010
1012 if n=2 then 1013 else 1014
1013 for t1=105 to 150 step 2
goto 1010
1014 for t1=95 to 150 step 2
1010
for QTEP=QTEPH(i) to QTEPK(i) step hag
GCH(i)=Q(i)/(4.2*(t1-t2)): ro=975
NOPT(i)=10000000
for KPD=0.6 to 0.75 step 0.01
NCH(i)=GCH(i)*dP(i)/(ro*1000*KPD)
if NCH(i)<NOPT(i) then NOPT(i)=NCH(i) else NOPT(i)=NOPT(i)
next KPD
NCH(i)=NOPT(i)
NPH(i)=0.1*NCH(i)
NCM(i)=0.1*NCH(i)
NCPK(i)=0.05*NCH(i)
if t2>=70 then koee=0.93 else koee=1
u1(i)=nCH(i)*p(12,i)*NOM*koee
U2(i)=NPH(i)*p(12,i)*NOM*koee
U3(i)=NCM(i)*p(12,i)*NOM*M1(i)*koee
U4(i)=NCPK(i)*p(12,i)*No*M1(i)*koee
U5(i)=NCPK(i)*p(12,i)*M1(i)*MuTB(i)*koee
if i=nnn then u4(i)=0 else u2(i)=0
if i=nnn then u5(i)=0 else u3(i)=0
UPER(i)=U1(i)+u2(I)+U3(i)+U4(i)+U5(i)
koff=0.05: tuh=1000
if t1>150 then 8123 else 8223
8123 koef=(t1-150)/200: koef=1-koef: goto 8323
8223 koef=1
8323
Utop(i)=NOM*k(i)*p(10,i)*(Q(i)+(QTEP*tuh)*p(11,i)+(Q(i)-(QTEP*tuh))*koff)/tuh
Uxbo(i)=v(i)*zxbo(i)
Uobc(i)=1.27*Q(i)*p(16,i)*(pk(i)+p(15,i)+p(14,i))/tuh
PP(i,i)=p(1,i)*(p(5,i)+p(9,i))+p(2,i)*koef*(p(6,i)+p(9,i))+p(3,i)*(p(7,i)+p(9,i))+p(4,i)*(p(8,i)+p(9,i))+
64
Utop(i)+Uper(i)+Uxbo(i)+UPER(i)
if PP(i,i)<PPmin then 1018 else 1019
1018
nomer=i: PPmin=PP(i,i): ppmin(i,i)=pp(i,i): t1min=t1: t2min=t2: qmin(i)=q(i): qumin(i)=qu(i)
for j=1 to 17
pmin(j,i)=p(j,i)
next j
QminTEP=QTEP: zminXBO(i)=zXBO(i): M1min(i)=m1(i): kmin(i)=k(i): vmin(i)=v(i): uper=uper(i)
1019
next qtep
if n=1 then next t1 else 1214
goto 1213
1214 if n=2 then next t1 else next t1
1213
locate 20,68:?" ### %"
locate 20,69:? proc
next t2
next i
1314 CLS
color 7,0
color 12,1
locate 1,2: ? " Паспорт оптимального варианта "
if n>1.1 then 331
if nomer=1 then 1 else 2
1
color 12,1: ii=1
locate 3,2: ? using " Приведенные затраты на ТЭЦ c ЦТП равны ############ руб/год
";ppmin(1,1)
goto 3
2
color 12,1: ii=2
locate 3,2: ? using " Приведенные затраты на ТЭЦ с ИТП равны ############ руб/год
";PPmin(2,2)
goto 3
331 if n>2.1 then 41
65
if nomer=3 then 111 else 21
111
color 12,1: ii=3
locate 3,2: ? using " Приведенные затраты на РК c ЦТП равны ################ руб/год
";ppmin(3,3)
goto 3
21
color 12,1: ii=4
locate 3,2: ? using " Приведенные затраты на РК с ИТП равны ################ руб/год
";PPmin(4,4)
goto 3
41
if nomer=5 then 1111 else 222
1111
color 12,1: ii=5
locate 3,2: ? using "Приведенные затраты на РК с ИТП и БА равны ###############
руб/год";ppmin(5,5)
goto 3
222
color 12,1: ii=6
locate 3,2: ? using "Приведенные затраты на РК с ИТП и ПАУ равны ##############
руб/год";PPmin(6,6)
3
i=ii
locate 5,2: ? using " Число жителей ### тыс человек ";g
locate 7,2: ? using " Удельный расход теплоты #.## кВт/чел ";qumin(i)
locate 9,2:? using " Суммарный расход теплоты #### МВт ";qmin(i)/1000
locate 11,2:? using " Начальная ### C и конечная ## C температура теплоносителя в
сети";t1min,t2min
locate 13,2:? using " Капиталовложения в источник теплоты #### руб/кВт ";pmin(1,i)
locate 15,2:? using " Капиталовложения в тепловые сети ### руб/кВт ";pmin(2,i)
locate 17,2:? using " Капиталовложения в тепловые пункты ### руб/кВт ";pmin(3,i)
locate 19,2:? using " Капиталовложения в местные системы теплопотребления ### руб/кВт
";pmin(4,i)
color 0,0
66
input ,y
cls
color 7,0
color 12,1
locate 1,2:? using " Доля ежегодных отчислений в источник теплоты # 1/год ";pmin(5,i)
locate 3,2:? using " Доля ежегодных отчислений в тепловые сети # 1/год ";pmin(6,i)
locate 5,2:? using " Доля ежегодных отчислений в тепловые пункты #.# 1/год ";pmin(7,i)
if n=>2 then 5
if ii=2 then 4
locate 7,2:? using " Число ЦТП ### шт ";M1min(i)
goto 56
4
locate 7,2:? using " Число ИТП ### шт ";M1min(i)
goto 56
5
if n=>3 then 55
if ii=4 then 44
locate 7,2:? using " Число ЦТП ### шт ";M1min(i)
goto 56
44
locate 7,2:? using " Число ИТП ### шт ";M1min(i)
goto 56
55
locate 7,2:? using " Число ИТП ### шт ";M1min(i)
56 locate 9,2:? using " Доля ежегодных отчислений в местные системы # 1/год ";pmin(8,i)
locate 11,2:? using " Нормативный коэффициент # 1/год ";pmin(9,i)
locate 13,2:? using " Удельн расход ###.## кг/(кВт ч) и затраты ### руб/тут на
топливо";kmin(i),pmin(10,i)
locate 15,2:? using " Длительность работы системы #### ч/год ";NOM
locate 17,2:? using " Коэффициент разветвленности сети #.## ";pmin(11,i)
locate 19,2:? using " Затраты на электрическую энергию #.## руб/(кВт ч) ";pmin(12,i)
color 0,0
input ,y
cls
color 7,0
67
color 12,1
locate 1,2:? using " Длительность работы ГВ #### ч/год и отопления #### ч/год
";pmin(13,i),No
locate 3,2:? using " Количество подпитки на ХВО ######## м3/год ";vmin(i)
locate 5,2:? using " Удельный расход на ХВО # руб/м3 ";zminXBO(i)
locate 7,2:? using " Среднегодовая заработная плата #### руб/год ";pmin(16,i)
locate 9,2:? using " Теплопотери ##.# МВт в тепловой сети с числом трубопроводов #
шт";QminTEP,pmin(17,i)
locate 12,2: ? " оптимальный паспорт записан под именем 2001ИТП.DAT "
color 0,0
input ,y
open "2001итп.dat" for output AS#2
?#2, " Паспорт оптимального варианта "
?#2,
if n>1.1 then 9331
if nomer=1 then 91 else 92
91 ii=1
?#2, using " Приведенные затраты на ТЭЦ c ЦТП равны ############ руб/год; ";ppmin(1,1)
goto 93
92 ii=2
?#2, using " Приведенные затраты на ТЭЦ с ИТП равны ############ руб/год; ";PPmin(2,2)
goto 93
9331 if n>2.1 then 941
if nomer=3 then 9111 else 921
9111 ii=3
?#2, using " Приведенные затраты на РК c ЦТП равны ################ руб/год; ";ppmin(3,3)
goto 93
921 ii=4
?#2, using " Приведенные затраты на РК с ИТП равны ################ руб/год; ";PPmin(4,4)
goto 93
941 if nomer=5 then 991 else 9222
991 ii=5
?#2, using " Приведенные затраты на РК с ИТП и БА равны ###############
руб/год;";ppmin(5,5)
68
goto 93
9222 ii=6
?#2, using " Приведенные затраты на РК с ИТП и ПАУ равны ##############
руб/год;";PPmin(6,6)
93 i=ii
?#2, using " Число жителей ### тыс человек; ";g
?#2, using " Удельный расход теплоты #.## кВт/чел; ";qumin(i)
?#2, using " Суммарный расход теплоты #### МВт; ";qmin(i)/1000
?#2, using " Начальная ### C и конечная ## C температура теплоносителя в сети;";t1min,t2min
?#2, using " Капиталовложения в источник теплоты #### руб/кВт; ";pmin(1,i)
?#2, using " Капиталовложения в тепловые сети ### руб/кВт; ";pmin(2,i)
?#2, using " Капиталовложения в тепловые пункты ### руб/кВт; ";pmin(3,i)
?#2, using " Капиталовложения в местные системы теплопотребления ### руб/кВт; ";pmin(4,i)
?#2, using " Доля ежегодных отчислений в источник теплоты # 1/год; ";pmin(5,i)
?#2, using " Доля ежегодных отчислений в тепловые сети # 1/год; ";pmin(6,i)
?#2, using " Доля ежегодных отчислений в тепловые пункты #.# 1/год; ";pmin(7,i)
if n=>2 then 95
if ii=2 then 94
?#2, using " Число ЦТП ### шт; ";M1min(i)
goto 956
94
?#2, using " Число ИТП ### шт; ";M1min(i)
goto 956
95
if n=>3 then 955
if ii=4 then 944
?#2, using " Число ЦТП ### шт; ";M1min(i)
goto 956
944
?#2, using " Число ИТП ### шт; ";M1min(i)
goto 956
955
?#2, using " Число ИТП ### шт; ";M1min(i)
956
?#2, using " Доля ежегодных отчислений в местные системы # 1/год; ";pmin(8,i)
69
?#2, using " Нормативный коэффициент # 1/год; ";pmin(9,i)
?#2, using " Удельн расход ###.## кг/(кВт ч) и затраты ### руб/тут на
топливо;";kmin(i),pmin(10,i)
?#2, using " Длительность работы системы #### ч/год; ";NOM
?#2, using " Коэффициент разветвленности сети #.##; ";pmin(11,i)
?#2, using " Затраты на электрическую энергию #.## руб/(кВт ч); ";pmin(12,i)
?#2, using " Длительность работы ГВ #### ч/год; ";pmin(13,i)
?#2, using " Количество подпитки на ХВО ######## м3/год; ";vmin(i)
?#2, using " Удельный расход на ХВО # руб/м3; ";zminXBO(i)
?#2, using " Среднегодовая заработная плата #### руб/год; ";pmin(16,i)
?#2, using " Теплопотери ##.# МВт в тепловой сети с числом трубопроводов # шт
";QminTEP,pmin(17,i)
end
70
Приложение Б
Результаты расчетов оптимальных вариантов
Паспорт оптимального варианта
Приведенные затраты на ТЭЦ с ИТП равны 145371968 руб./год;
Число жителей 100 тыс. человек;
Удельный расход теплоты 1.50 кВт/чел;
Суммарный расход теплоты 150 МВт;
Начальная 150 оC и конечная 40 оC температура теплоносителя в сети;
Капиталовложения в источник теплоты 1000 руб./кВт;
Капиталовложения в тепловые сети 360 руб./кВт;
Капиталовложения в тепловые пункты 550 руб./кВт;
Капиталовложения в местные системы теплопотребления 250 руб./кВт;
Доля ежегодных отчислений в источник теплоты 7 1/год;
Доля ежегодных отчислений в тепловые сети 5 1/год;
Доля ежегодных отчислений в тепловые пункты 6.5 1/год;
Число ИТП 70 шт.;
Доля ежегодных отчислений в местные системы 6 1/год;
Нормативный коэффициент 8 1/год;
Удельный расход 0.35 кг/(кВт ч) и затраты 300 руб./т.у.т. на топливо;
Длительность работы системы 8400 ч/год;
Коэффициент разветвленности сети 1.50;
Затраты на электрическую энергию 0.05 руб./(кВт ч);
Длительность работы ГВ 8400 ч/год;
Количество подпитки на ХВО 1000000 м3/год;
Удельный расход на ХВО 1 руб./м3;
Среднегодовая заработная плата 5000 руб./год;
Теплопотери 3.8 МВт в тепловой сети с числом трубопроводов 2 шт.
Паспорт оптимального варианта
Приведенные затраты на ТЭЦ с ИТП равны 374646240 руб./год;
Число жителей 150 тыс. человек;
Удельный расход теплоты 1.69 кВт/чел;
Суммарный расход теплоты 319 МВт;
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]