В качестве экономической задачи в данной работе нужно рассчитать плановую себестоимость передачи тепловой энергии тепловыми сетями с учетом прогнозирования возможных причин отказов, оказывающих большое влияние на величину планируемых расходов (затрат) на ремонт трубопроводов теплосетей и потерь тепловой энергии с утечками. Характерными причинами отказов являются: разрыв трубы, разрыв сварного шва, свищ, отказ арматуры, повреждение сальникового компенсатора.
В качестве метода прогнозирования, возможных причин отказов тепловых сетей в работе используется метод прогнозирования на основе выявления временного тренда, путем аналитического выравнивания уровней динамического ряда.
Глава 1. ПРОГНОЗИРОВАНИЕ ОТКАЗОВ ТЕПЛОВЫХ СЕТЕЙ
Методика статистического прогноза по тренду основана на предположении, что параметры тренда и колебаний сохраняются до прогнозируемого периода.
Данная методика расчета включает в себя несколько видов:
· Выявление типа тренда;
· Вывод уравнения тренда;
· Анализ полученной зависимости;
· Прогноз на требуемый период.
1.1 Выявление типа тренда
Для определения основной тенденции развития значения эмпирических уровней их изображают на графике y = f (t).
В качестве динамического ряда берется изменение количества отказов тепловых сетей за предлагаемый период времени в зависимости от характера отказа. Эмпирические уровни ряда сглаживаются методом скользящей средней. Сглаженные уровни наносят на график.
Характер отказа – Свищ
Характер отказа – Разрыв трубы
Характер отказа – Разрыв сварного шва
Характер отказа – Повреждение сальникового компенсатора
Характер отказа – отказ арматуры
1.2 ВЫВОД УРАВНЕНИЯ ТРЕНДА
На графическом изображении сглаженных уровней отчетливо проявляется основная тенденция изменения количества отказов тепловых сетей по различным причинам.
Следующим этапом работы является получение обобщенной статистической оценки тренда методом аналитического выравнивания.
Основная тенденция развития yt
как функция времени:
yti
= f (ti
).
Определение теоретических (расчетных) уровней yti
производится на основе адекватной математической модели, которая наиболее точно отображает основную тенденцию ряда динамики – количество отказов теплосети по данной причине.
Подбор математической функции основан на определении типа развития экономического явления во времени.
Когда тип тренда установлен, необходимо вычислить оптимальные значения параметров тренда, исходя из фактических уровней, т.е. получить значения коэффициентов уравнения. На основе требований метода наименьших квадратов (МНК) составляется система нормальных уравнений, решая которую, вычисляем параметры тренда.
Для линейного тренда нормальные уравнения МКТ имеют вид:
n*a + b*St i
= Sy i
a*St i
+ b*St2
i
= S yi
*ti
где yi
– уровни исходного ряда динамики;
ti
– номера периодов или моментов времени;
a, b – константы уравнений.
Нормальные уравнения МКТ для тренда в виде параболы 2-го порядка:
n*a + b*Sti
+ c*Sti
2
= Syi
a*Sti
+ b*Sti
2
+ c*Sti
3
= Syi
*ti
a*Sti
2
+ b*Sti
3
+ c*Sti
4
= Syi
*ti
2
где yi
– уровни ряда динамики;
ti
– номера периодов или моментов времени;
a, b, c – константы уравнений.
Таблица 1.2.1 Формирование данных для вывода уравнения тренда – характер отказа свищ
| свищ | Yi | скользящая средняя(3 г.) | Ti | Yi*Ti | Ti*Ti | Yi*Ti 2 | Ti 3 | Ti 4 | Ti 2 |
| 1991 | 25 | - | 1 | 25 | 1 | 25 | 1 | 1 | 1 |
| 1992 | 31 | 30,33333333 | 2 | 62 | 4 | 124 | 8 | 16 | 4 |
| 1993 | 35 | 31,66666667 | 3 | 105 | 9 | 315 | 27 | 81 | 9 |
| 1994 | 29 | 31,66666667 | 4 | 116 | 16 | 464 | 64 | 256 | 16 |
| 1995 | 31 | 30 | 5 | 155 | 25 | 775 | 125 | 625 | 25 |
| 1996 | 30 | 32,66666667 | 6 | 180 | 36 | 1080 | 216 | 1296 | 36 |
| 1997 | 37 | 33,33333333 | 7 | 259 | 49 | 1813 | 343 | 2401 | 49 |
| 1998 | 33 | 34,33333333 | 8 | 264 | 64 | 2112 | 512 | 4096 | 64 |
| 1999 | 33 | 33,66666667 | 9 | 297 | 81 | 2673 | 729 | 6561 | 81 |
| 2000 | 35 | - | 10 | 350 | 100 | 3500 | 1000 | 10000 | 100 |
| итого | 319 | 257,6666667 | 55 | 1813 | 385 | 12881 | 3025 | 25333 | 385 |
Для линейного тренда:
10*a +55*b = 319; a = 28;
55*a + 385*b = 1813; b = 0,71;
y = 28 + 0,71*t;
Для тренда в виде параболы второго порядка:
10*a + 55*b + 385*c = 319; a = 26,02;
55*a + 385*b + 3025*c = 1813; b = 1,7;
385*a + 3025*b + 25333*c = 12881; c = -0,09;
y = 26,02 + 1,7*t –0,09*t2
;
Таблица 1.2.2 Формирование данных для вывода уравнения тренда – характер отказа разрыв трубы
| разрыв трубы | Yi | скольз. ср. за 3 г. | Ti | Yi*Ti | Ti*Ti | Yi*Ti 2 | Ti 3 | Ti 4 | Ti 2 |
| 1991 | 12 | 1 | 12 | 1 | 12 | 1 | 1 | 1 | |
| 1992 | 8 | 10 | 2 | 16 | 4 | 32 | 8 | 16 | 4 |
| 1993 | 10 | 9,666667 | 3 | 30 | 9 | 90 | 27 | 81 | 9 |
| 1994 | 11 | 11,33333 | 4 | 44 | 16 | 176 | 64 | 256 | 16 |
| 1995 | 13 | 11,33333 | 5 | 65 | 25 | 325 | 125 | 625 | 25 |
| 1996 | 10 | 11,66667 | 6 | 60 | 36 | 360 | 216 | 1296 | 36 |
| 1997 | 12 | 11,33333 | 7 | 84 | 49 | 588 | 343 | 2401 | 49 |
| 1998 | 12 | 13 | 8 | 96 | 64 | 768 | 512 | 4096 | 64 |
| 1999 | 15 | 12,66667 | 9 | 135 | 81 | 1215 | 729 | 6561 | 81 |
| 2000 | 11 | 10 | 110 | 100 | 1100 | 1000 | 10000 | 100 | |
| итого | 114 | 91 | 55 | 652 | 385 | 4666 | 3025 | 25333 | 385 |
Для линейного тренда:
y = 9,75 + 0,3*t;
Для тренда в виде параболы второго порядка:
y = 9,82 + 0,26*t + 0,004*t2
;
Таблица 1.2.3 Формирование данных для вывода уравнения тренда – характер отказа разрыв сварного шва.
| разрыв св.шва | Yi | ср. скользящая за 3 г. | Ti | Yi*Ti | Ti*Ti | Yi*Ti 2 | Ti 3 | Ti 4 | Ti 2 |
| 1991 | 10 | 1 | 10 | 1 | 10 | 1 | 1 | 1 | |
| 1992 | 6 | 7,666667 | 2 | 12 | 4 | 24 | 8 | 16 | 4 |
| 1993 | 7 | 6,666667 | 3 | 21 | 9 | 63 | 27 | 81 | 9 |
| 1994 | 7 | 8,333333 | 4 | 28 | 16 | 112 | 64 | 256 | 16 |
| 1995 | 11 | 9 | 5 | 55 | 25 | 275 | 125 | 625 | 25 |
| 1996 | 9 | 10,33333 | 6 | 54 | 36 | 324 | 216 | 1296 | 36 |
| 1997 | 11 | 11,66667 | 7 | 77 | 49 | 539 | 343 | 2401 | 49 |
| 1998 | 15 | 12 | 8 | 120 | 64 | 960 | 512 | 4096 | 64 |
| 1999 | 10 | 12,33333 | 9 | 90 | 81 | 810 | 729 | 6561 | 81 |
| 2000 | 12 | 10 | 120 | 100 | 1200 | 1000 | 10000 | 100 | |
| итого | 98 | 78 | 55 | 587 | 385 | 4317 | 3025 | 25333 | 385 |
Для линейного тренда:
y = 6,61+ 0,58*t;
Для тренда в виде параболы второго порядка:
y = 7,27+ 0,25*t + 0,03*t2
;
Таблица 1.2.4 Формирование данных для вывода уравнения тренда – характер отказа повреждение сальникового компенсатора
| С/к | Yi | ср скользящая за 3 г. | Ti | Yi*Ti | Ti*Ti | Yi*Ti 2 | Ti 3 | Ti 4 | Ti 2 |
| 1991 | 7 | 1 | 7 | 1 | 7 | 1 | 1 | 1 | |
| 1992 | 12 | 9,666667 | 2 | 24 | 4 | 48 | 8 | 16 | 4 |
| 1993 | 10 | 11 | 3 | 30 | 9 | 90 | 27 | 81 | 9 |
| 1994 | 11 | 12 | 4 | 44 | 16 | 176 | 64 | 256 | 16 |
| 1995 | 15 | 12,66667 | 5 | 75 | 25 | 375 | 125 | 625 | 25 |
| 1996 | 12 | 13 | 6 | 72 | 36 | 432 | 216 | 1296 | 36 |
| 1997 | 12 | 10,33333 | 7 | 84 | 49 | 588 | 343 | 2401 | 49 |
| 1998 | 7 | 9,666667 | 8 | 56 | 64 | 448 | 512 | 4096 | 64 |
| 1999 | 10 | 10 | 9 | 90 | 81 | 810 | 729 | 6561 | 81 |
| 2000 | 13 | 10 | 130 | 100 | 1300 | 1000 | 10000 | 100 | |
| итого | 109 | 88,33333 | 55 | 612 | 385 | 4274 | 3025 | 25333 | 385 |
y = 9,69 + 0,22*t;
Для тренда в виде параболы второго порядка:
y = 7,66 + 1,36*t - 0,11*t2
;
Таблица 1.2.5 Формирование данных для вывода уравнения тренда – характер отказа отказ арматуры
| Отказ арматуры | Yi | ср. скользящая за 3 г. | Ti | Yi*Ti | Ti*Ti | Yi*Ti 2 | Ti 3 | Ti 4 | Ti 2 |
| 1991 | 6 | 1 | 0 | 0 | 0 | 0 | 0 | 1 | |
| 1992 | 8 | 7,333333 | 2 | 8 | 1 | 8 | 1 | 1 | 4 |
| 1993 | 8 | 8,666667 | 3 | 16 | 4 | 32 | 8 | 16 | 9 |
| 1994 | 10 | 9 | 4 | 30 | 9 | 90 | 27 | 81 | 16 |
| 1995 | 9 | 10,33333 | 5 | 36 | 16 | 144 | 64 | 256 | 25 |
| 1996 | 12 | 10,33333 | 6 | 60 | 25 | 300 | 125 | 625 | 36 |
| 1997 | 10 | 13 | 7 | 60 | 36 | 360 | 216 | 1296 | 49 |
| 1998 | 17 | 13,66667 | 8 | 119 | 49 | 833 | 343 | 2401 | 64 |
| 1999 | 14 | 14,66667 | 9 | 112 | 64 | 896 | 512 | 4096 | 81 |
| 2000 | 13 | 10 | 117 | 81 | 1053 | 729 | 6561 | 100 | |
| 107 | 87 | 55 | 558 | 285 | 3716 | 2025 | 15333 | 385 |
Для линейного тренда:
y = 12,74 - 0,37*t;
Для тренда в виде параболы второго порядка:
y = 9,9 + 1,05*t - 0,129*t2
;
1.3 Анализ полученной зависимости
Одним из применяемых на практике показателей адекватности математической функции для статистического изучения тренда является стандартизированная ошибка аппроксимации - syt
:
где уti
- теоретические уровни;
yi
- экспериментальные уровни;
n - число уровней ряда.
За наиболее адекватную принимается функция, у которой syt
минимальная.
Правильность окончательных расчетов проверяется по равенству:
åyi
= å уti
Таблица 1.3.1 Расчет параметров тренда (характер отказа - свищ)
| Yi | Ti | прям | параб | отклонения | прям | параб | |
| Yti-Yi | Yti-Yi 2 | Yti-Yi | Yti-Yi 2 | ||||
| 25 | 1 | 28,71 | 27,63 | 3,71 | 13,7641 | 2,63 | 6,9169 |
| 31 | 2 | 29,42 | 29,06 | -1,58 | 2,4964 | -1,94 | 3,7636 |
| 35 | 3 | 30,13 | 30,31 | -4,87 | 23,7169 | -4,69 | 21,9961 |
| 29 | 4 | 30,84 | 31,38 | 1,84 | 3,3856 | 2,38 | 5,6644 |
| 31 | 5 | 31,55 | 32,27 | 0,55 | 0,3025 | 1,27 | 1,6129 |
| 30 | 6 | 32,26 | 32,98 | 2,26 | 5,1076 | 2,98 | 8,8804 |
| 37 | 7 | 32,97 | 33,51 | -4,03 | 16,2409 | -3,49 | 12,1801 |
| 33 | 8 | 33,68 | 33,86 | 0,68 | 0,4624 | 0,86 | 0,7396 |
| 33 | 9 | 34,39 | 34,03 | 1,39 | 1,9321 | 1,03 | 1,0609 |
| 35 | 10 | 35,1 | 34,02 | 0,1 | 0,01 | -0,98 | 0,9604 |
| 319 | 55 | 319,05 | 319,05 | 0,05 | 67,4185 | 0,05 | 63,7753 |
sпрямая
= 2,5965; sпарабола
= 2,5254;
Уравнение тренда имеет вид:
y = 26,02 +1,7*t-0,09 t2
;
Таблица 1.3.2 Расчет параметров тренда (характер отказа – разрыв трубы)
| Yi | Ti | прям | параб | отклонения | прям | параб | |
| Yti-Yi | Yti-Yi 2 | Yti-Yi | Yti-Yi 2 | ||||
| 12 | 1 | 10,05 | 10,08 | -1,950 | 3,803 | -1,916 | 3,671 |
| 8 | 2 | 10,35 | 10,36 | 2,350 | 5,523 | 2,356 | 5,551 |
| 10 | 3 | 10,65 | 10,64 | 0,650 | 0,423 | 0,636 | 0,404 |
| 11 | 4 | 10,95 | 10,92 | -0,050 | 0,003 | -0,076 | 0,006 |
| 13 | 5 | 11,25 | 11,22 | -1,750 | 3,063 | -1,780 | 3,168 |
| 10 | 6 | 11,55 | 11,52 | 1,550 | 2,403 | 1,524 | 2,323 |
| 12 | 7 | 11,85 | 11,84 | -0,150 | 0,023 | -0,164 | 0,027 |
| 12 | 8 | 12,15 | 12,16 | 0,150 | 0,023 | 0,156 | 0,024 |
| 15 | 9 | 12,45 | 12,48 | -2,550 | 6,503 | -2,516 | 6,330 |
| 11 | 10 | 12,75 | 12,82 | 1,750 | 3,063 | 1,820 | 3,312 |
| 114 | 55 | 114,00 | 114,04 | 0,000 | 24,825 | 0,04 | 24,817 |
sпрямая
= 1,5756;
sпарабола
= 1,5753;
Уравнение тренда имеет вид:
y = 9,82 + 0,26*t + 0,004*t2
;
Таблица 1.3.3 Расчет параметров тренда (характер отказа – разрыв сварного шва)
| Yi | Ti | прям | параб | отклонения | прям | параб | |
| Yti-Yi | Yti-Yi 2 | Yti-Yi | Yti-Yi 2 | ||||
| 10 | 1 | 7,19 | 7,550 | -2,8100 | 7,8961 | -2,4500 | 6,0025 |
| 6 | 2 | 7,77 | 7,890 | 1,7700 | 3,1329 | 1,8900 | 3,5721 |
| 7 | 3 | 8,35 | 8,290 | 1,3500 | 1,8225 | 1,2900 | 1,6641 |
| 7 | 4 | 8,93 | 8,750 | 1,9300 | 3,7249 | 1,7500 | 3,0625 |
| 11 | 5 | 9,51 | 9,270 | -1,4900 | 2,2201 | -1,7300 | 2,9929 |
| 9 | 6 | 10,09 | 9,850 | 1,0900 | 1,1881 | 0,8500 | 0,7225 |
| 11 | 7 | 10,67 | 10,490 | -0,3300 | 0,1089 | -0,5100 | 0,2601 |
| 15 | 8 | 11,25 | 11,190 | -3,7500 | 14,0625 | -3,8100 | 14,5161 |
| 10 | 9 | 11,83 | 11,950 | 1,8300 | 3,3489 | 1,9500 | 3,8025 |
| 12 | 10 | 12,41 | 12,770 | 0,4100 | 0,1681 | 0,7700 | 0,5929 |
| 98 | 55 | 98,00 | 98,00 | 0,0000 | 37,67 | 0,0000 | 37,18
8 |
sпрямая
= 1,94087;
sпарабола
= 1,92842;
Уравнение тренда имеет вид:
y = 7,27+ 0,25*t + 0,03*t2
;
Таблица 1.3.4 Расчет параметров тренда (характер отказа – повреждение сальникового компенсатора)
| Yi | Ti | прям | параб | отклонения | прям | параб | |
| Yti-Yi | Yti-Yi 2 | Yti-Yi | Yti-Yi 2 | ||||
| 7 | 1 | 9,910 | 8,910 | 2,9100 | 8,4681 | 1,9100 | 3,6481 |
| 12 | 2 | 10,130 | 9,940 | -1,8700 | 3,4969 | -2,0600 | 4,2436 |
| 10 | 3 | 10,350 | 10,750 | 0,3500 | 0,1225 | 0,7500 | 0,5625 |
| 11 | 4 | 10,570 | 11,340 | -0,4300 | 0,1849 | 0,3400 | 0,1156 |
| 15 | 5 | 10,790 | 11,710 | -4,2100 | 17,7241 | -3,2900 | 10,8241 |
| 12 | 6 | 11,010 | 11,860 | -0,9900 | 0,9801 | -0,1400 | 0,0196 |
| 12 | 7 | 11,230 | 11,790 | -0,7700 | 0,5929 | -0,2100 | 0,0441 |
| 7 | 8 | 11,450 | 11,500 | 4,4500 | 19,8025 | 4,5000 | 20,2500 |
| 10 | 9 | 11,670 | 10,990 | 1,6700 | 2,7889 | 0,9900 | 0,9801 |
| 13 | 10 | 11,890 | 10,260 | -1,1100 | 1,2321 | -2,7400 | 7,5076 |
| 109 | 55 | 109,00 | 109,05 | 0,0000 | 55,3930 | 0,0500 | 48,1953 |
sпрямая
= 2,3536;
sпарабола
= 2,1953;
Уравнение тренда имеет вид:
y = 7,66 + 1,36*t - 0,11*t2
;
Таблица 1.3.5 Расчет параметров тренда (характер отказа – отказ арматуры)
| Yi | Ti | прям | параб | отклонения | прям | параб | |
| Yti-Yi | Yti-Yi 2 | Yti-Yi | Yti-Yi 2 | ||||
| 6 | 1 | 12,370 | 10,821 | 6,3700 | 40,5769 | 4,8210 | 23,2420 |
| 8 | 2 | 12,000 | 11,484 | 4,0000 | 16,0000 | 3,4840 | 12,1383 |
| 8 | 3 | 11,630 | 11,889 | 3,6300 | 13,1769 | 3,8890 | 15,1243 |
| 10 | 4 | 11,260 | 12,036 | 1,2600 | 1,5876 | 2,0360 | 4,1453 |
| 9 | 5 | 10,890 | 11,925 | 1,8900 | 3,5721 | 2,9250 | 8,5556 |
| 12 | 6 | 10,520 | 11,556 | -1,4800 | 2,1904 | -0,4440 | 0,1971 |
| 10 | 7 | 10,150 | 10,929 | 0,1500 | 0,0225 | 0,9290 | 0,8630 |
| 17 | 8 | 9,780 | 10,044 | -7,2200 | 52,1284 | -6,9560 | 48,3859 |
| 14 | 9 | 9,410 | 8,901 | -4,5900 | 21,0681 | -5,0990 | 25,9998 |
| 13 | 10 | 9,040 | 7,500 | -3,9600 | 15,6816 | -5,5000 | 30,2500 |
| 107 | 107,05 | 107,085 | 0,0500 | 166,0045 | 0,0850 | 168,9015 |
sпрямая
= 4,0744;
sпарабола
= 4,1097;
Уравнение тренда имеет вид:
y = 12,74 - 0,37*t;
1.4 ПРОГНОЗ НА ТРЕБУЕМЫЙ ПЕРИОД
В результате анализа трендов получим математические функции, характеризующие изменение количества отказов тепловых сетей по различным причинам от времени.
Количество полученных уравнений должно соответствовать количеству причин отказов тепловых сетей.
Для прогнозирования на требуемый период времени необходимо рассчитать, какому периоду времени соответствует выбранный год.
Номер периода времени определяется следующим образом:
t = ( Tпрог
– Tно
) +1
где Tпрог
- время (дата), на которое делается прогноз, год:
Tно
- время (дата) начального отсчета, год.
В уравнение тренда надо подставить полученное значение t =... В результате получаем количество отказов теплосети по выбранной причине отказа.
t= (2001 – 1991)+1=11;
Таблица 1.4.1 Прогнозируемое количество отказов теплосети на 2000 год
| Причины отказов теплосети | Количество отказов |
| Свищ | 34 |
| Разрыв трубы | 14 |
| Разрыв сварного шва | 14 |
| Повреждение сальникового компенсатора | 10 |
| Отказ арматуры | 9 |
Глава 2. ОПРЕДЕЛЕНИЕ ИЗДЕРЖЕК ПРОИЗВОДСТВА ПРИ ПЕРЕДАЧЕ ТЕПЛОВОЙ ЭНЕРГИИ
Годовые издержки производства (эксплуатационные расходы), связанные с передачей тепловой энергии, включают в себя следующие статьи:
- затраты на подпитку тепловых сетей;
- затраты на перекачку теплоносителя;
- фонд оплаты труда обслуживающего персонала;
- амортизационные отчисления на полное восстановление (реновацию) основных производственных фондов теплосети;
- затраты на текущие и капитальные ремонты основного и вспомогательного оборудования тепловых сетей;
- экономический ущерб, связанный с отказами тепловых сетей;
- общецеховые расходы;
- общесетевые расходы.
2.1 Расчет затрат на подпитку тепловых сетей
Иподт
= Qпот
/(Cтн
*tтн
)*Цхов,
Где Qпот
- количество потерянной тепловой энергии в связи с отказами теплосети, Гкал/год;
Cтн
- теплоемкость теплоносителя, Ккал/кг*С; для воды принимается равной Ккал/кг*С;
tтн
- температура теплоносителя. С; принимается в пределах 60+70 С;
Цхов
- стоимость химически очищенной воды предприятия, руб/кг; принимается по действующим тарифам.
Годовое нормированное количество потерянной тепловой энергии (утечки), Гкал/год,
Qпот
= 0,025*Qпер
,
где 0,025 - нормированная доля утечек теплоносителя, связанная с отказами в работе теплосети для Ивгортеплосети;
Qпер
- годовое количество теплоты, передаваемой потребителям, Гкал/год.
Qпот
= 0,025* Qпер
=0,025*3,38*106
= 0,0845*106
Гкал/год
Cтн
= 1 Ккал/кг*С; tтн
= 60 С; Цхов
= 0,0045 руб/кг;
Иподт
= 0,0845*106
/ (1*65)*0,0045*103
= 5,85 млн руб/год;
2.2 Затраты на перекачку теплоносителя
Ипер
=Wн
*Цээ
*10-6
где Wн
- количество электрической энергии, затрачиваемой на перекачку теплоносителя насосом, КВт-ч/год;
ЦЭЭ
- стоимость КВт-ч электрической энергии по двухставочному тарифу, руб/КВт*ч; принимается по действующим тарифам.
а =106,97 руб/кВт*ч – в месяц
б = 0,42 руб/кВт*ч
б = 0,316 руб/кВт*ч
h = 6000 часов
Цээ =0,71
2.2.1 Количество электрической энергии, потребляемой насосом
Wн
=( Qпот
+Qпер
)Dpн
*106
/ Cтн
*tтн
*3600 *hн
*rтн
где Dpн
- полное гидравлическое сопротивление насоса, Па,
rтн
- плотность теплоносителя ( воды ), кг/м3
,
hн
- коэффициент полезного действия насоса.
Dpн
=500 Па;
rтн
= 975 кг/м3
;
hн
= 0,6;
Wн
= ((0,0845+3,38)*500*1012
)/(1*65*975*3600*0,6)=15 432 098,77 кВт*ч/год
Ипер
= 15 432 098,77 * 0,71 = 10,95679013 млн руб/год
2.3 Расчет фонда оплаты труда обслуживающего персонала
В этой статье учитывается основная, дополнительная заработная плата и начисления на нее для производственных рабочих и других категорий работников, непосредственно участвующих в процессе передачи и распределения энергии ( определяется укрупненно ).
Фонд оплаты труда обслуживающего персонала тепловых сетей, млн руб/год,
Изп
= (Иозп
+ Идзп
+ Ин
)
Где Иозп
- основная заработная плата обслуживающего персонала теплосети, руб/год;
Идзп
дополнительная заработная плата работников теплосети, руб/год;
Ин
- начисления на зарплату, руб/год.
2.3.1 Основная зарплата обслуживающего персонала тепловых сетей
Иозп
= ЗП ср
* 12 * n
где ЗП р
- среднемесячная зарплата одного рабочего, обслуживающего тепловые сети, руб;
n - списочное количество работников, обслуживающих тепловые сети.
ЗПср
=1930 руб/мес ; n= 138 человек;
И озп
= 1930 * 12 * 163 = 3636120 руб/год;
2.3.2 Дополнительная зарплата
И дзп
= (И озп
/(12 * 25,25)) * Nср
Где Nср
- среднее количество дней отдыха, болезни и выполнения государственных обязанностей работников теплосети..
Nср
= 24 дня;
И дзп
= 3636120 / (12*25,25)*24 = 288009,5 руб/год;
2.3.3 Начисления на заработную плату
Ин
= j * (Иозп
+ Идзп
)/100,
где j - процент начислений на основную и дополнительную зарплату.
j = 35,6 %
Ин
= 35,6 * (3636120+ 288009,5)/100 = 1396990,102руб/год;
И зп
= 5,32112 млн. руб/год
2.4 Расчет фонда амортизации
По этой статье в себестоимости передачи тепловой энергии учитываются амортизационные отчисления на полное восстановление ( реновацию ) основного и вспомогательного оборудования, а также отчисления на реновацию и капитальный ремонт зданий и сооружений.
Сумма амортизационных отчислений по основным фондам теплосети,
Иам
= (Кбал
* аср
)/100
Где Кбал
- общая балансовая стоимость оборудования, зданий и т.п. теплосети;
аср
- средняя по основным фондам норма амортизации на полное восстановление (реновацию ), %.
Кбал
= 83,3 млн руб/год; аср
= 4 %;
Иам
= (83,3* 4)/100 = 3,332 млн руб/год;
2.5 Расчет затрат на текущие и капитальные ремонты оборудования
Текущие и капитальные ремонты финансируются из единого ремонтного фонда предприятия. Для определения затрат на ремонты составляется годовой календарный график ремонтов оборудования тепловых сетей с учетом норм периодичности имеющегося оборудования. При расчете затрат на ремонты используются нормативы трудоемкости ремонтов, нормы расхода вспомогательных материалов и энергетических ресурсов.
Ирем
= 1,5*Иам
= 1,5 * 5,321122 = 7,981683 млн руб/год;
2.6 Расчет экономического ущерба, связанного с отказами тепловых сетей
Общий экономический ущерб,
Иущ
= Иу1
+Иу2
+Иуз
,
где Иу1
- потери тепловой энергии в трубопроводах, связанные с утечками теплоносителя (упущенная выгода), тыс руб/год;
Иу2
- возмещение убытков потребителям тепловой энергии, тыс руб/год;
Иуз
- затраты на внеплановые ремонты оборудования тепловых сетей, тыс руб/год.
2.6 1 Потери тепловой энергии в трубопроводах, тыс руб/год
Иу1
= Qпот.отк
* Цтэ
где Qпот.отк
- количество потерянной тепловой энергии, связанной с отказами теплосети, Гкал в год;
Цтэ
- цена 1 Гкал теплоты, отпущенной потребителю, руб/Гкал;
Qпот
.
отк
= åqпот
i
* l i
* ni
* Ткал
i
где qпот
i
- удельное количество утечек тепловой энергии при данном характере отказа теплосети, Гкал /пог.м-день-единицу характера отказа, определяется по данным табл.4;
ni
- количество отказов теплосети данного характера, определяется по результатам обработки статистических данных;
Ткал i
- время, прошедшее с момента обнаружения отказа до его устранения, дни;
li
- длина замененного участка теплосети при устранении отказа, пог.м.
Таблица 2.6.1 Расчет Qпот.отк
для каждого характера отказа
| Характер отказа | qпот
|
n | Ткал
|
l | Qпот.отк
|
| Свищ | 0,05 | 34 | 6 | 0,4 | 4,08 |
| Разрыв трубы | 0,279 | 14 | 5 | 5 | 97,65 |
| Разрыв сварного шва | 0,015 | 14 | 5 | 3 | 3,15 |
| Повреждение сальникового компенсатора | 0,005 | 10 | 3 | 0,1 | 0,015 |
| Отказ арматуры | 0,001 | 9 | 3 | 0,2 | 0,0054 |
Qпот.отк
= 104,9 Гкал/год
Цтэ
= 227,04 руб/Гкал
Иу1
= 227,04 * 104,9 = 23816,5 руб/год
Иу1
=0,0238165 млн руб/год
2.6.2 Возмещение убытков потребителям тепловой энергии
Возмещение убытков потребителям за недоотпуск или снижение параметров теплоносителя определяется в соответствии с условиями коммерческих договоров между энергоснабжающей организацией и потребителем тепловой энергии, млн руб/год,
Иу2
= Цдог
* Qпот
. отк
.
где Цдог
- договорная цена промышленного ущерба, приходящаяся на 1 Гкал недоотпущенной тепловой энергии, тыс.руб/Гкал.
Цдог
= 227,64 руб/год
Иу2
= Иу1
=0,0238165 млн руб/год
2.6.3 Затраты на внеплановые ремонты оборудования
В эту статью включаются затраты на обнаружение места утечки, на земляные работы, непосредственное проведение ремонта (установка заплат, сварочные работы и т.п.) и не входит стоимость строительно-монтажных работ, связанных с заменой оборудования тепловых сетей.
Иуз
= åЦрем
i
* n i
* l i
Где Ц рем i
- удельная стоимость ремонтных работ по устранению утечек воды в теплосетях при данном характере отказа, тыс.руб/пог.м;
ni
- число отказов теплосети данного характера;
li
- длина замененного участка теплосети, пог.м.
Таблица 2.6.3 Расчет Иуз
для каждого характера отказа
| Характер отказа | n | l | Ц рем
|
Иуз
|
| Свищ | 34 | 0,5 | 0,9 | 15,3 |
| Разрыв трубы | 14 | 4 | 0,5 | 28 |
| Разрыв сварного шва | 14 | 3 | 0,5 | 21 |
| Повреждение сальникового компенсатора | 10 | 0,3 | 0,5 | 1,5 |
| Отказ арматуры | 9 | 0,3 | 0,5 | 1,35 |
Иуз
= 67,15 тыс руб/год
Иущ
=0,0238165 + 0,0238165 + 0,06715= 0,114783 млнруб/год
Общие эксплуатационные расходы, тыс руб/год,
Иэкс
= Иподп
+Ипер
+Изп
+Иам
+Ирем
+Иущ
Иэкс
= 5,85 +10,95679013 + 5,321122 + 3,332 + 7,981683 + 0,114783 = 33,556378 млн руб/ год
2.7 Общецеховые расходы
К этой статье относятся затраты по обслуживанию и управлению тепловыми сетями, то есть зарплата аппарата управления, затраты по текущему ремонту зданий, расходы по охране труда и т.д.
Цеховые расходы могут быть рассчитаны укрупненно, млн руб/год:
Ицех
= g*Иэкс
.
Где g - доля общецеховых расходов от суммарных эксплуатационных.
g = 0,1;
Ицех
= 0,1 * 33,556378 = 3,3556378 млн руб/год
2.8 Общесетевые расходы
В этой статье учитываются расходы по управлению предприятием тепловых сетей:
- административно-управленческие расходы.
- общепроизводственные расходы (НИОКР, текущий ремонт общесетевых средств и т.д.).
Исет
= ЗП адм
* nупр
+ k * (И экс
. + И цех
)
где ЗП адм
- среднегодовая заработная плата административно-управленческого персонала теплосетей, тыс руб/год;
nупр
- численность административно-управленческого персонала теплосетей, чел;
k - коэффициент, учитывающий общепроизводственные расходы.
ЗП адм
= 0,036 млн руб/год nупр
= 19 чел , k = 0,05;
Исет
= 0,036* 19 + 0,05* (3,3556378 +33,556378 ) = 2,529601 млн руб/год
Полные издержки производства при передаче тепловой энергии, тыс руб/год,
Иполн
= Иэкс
+ Ццэх
+ Цсет
Иполн
= 3,3556378 + 33,556378 +2,529601 = 39,4416168 млн руб/год
Глава 3. РАСЧЕТ СЕБЕСТОИМОСТИ ПЕРЕДАВАЕМОЙ ТЕПЛОВОЙ ЭНЕРГИИ И ПЛАНИРОВАНИЕ КОЛИЧЕСТВА ПОЛУЧАЕМОЙ ПРИБЫЛИ
Плановая себестоимость передачи тепловой энергии (без учета отказов тепловых сетей ), руб/Гкал,
Спл
= (Иполн
– Иущ
)*106
/Qпер
Спл
= (39,4416168 –0,114783 )* 106
/ 3,38*106
= 11,63 руб/Гкал
Плановая себестоимость передачи тепловой энергии (с учетом отказов тепловых сетей), руб/Гкал,
Cпл
от
= Иполн
* 106
/ Qпер
Cпл
от
= 39,4416168 * 106
/ 3,38*106
= 11,67 руб/Гкал
Планируемое снижение прибыли тепловых сетей в связи с отказами, млн руб/год.
DП = П2
– П1
где П2
, П1
- прибыль тепловых сетей при различных значениях себестоимости передаваемой тепловой энергии.
П2
= (Ттэ
- Скот
- Спл
)*Qпер
П1
= (Ттэ
- Скот
- Cпл
от
)*Qпер
где Ттэ
- тариф на тепловую энергию, установленный региональной энергетической комиссией на рассматриваемый период времени, руб/Гкал;
Скот
- себестоимость тепловой энергии, отпускаемой котельной иди ТЭЦ теплосети, руб/Гкал.
Скот
= 60 руб/Гкал
П2
= (189,7 – 162,7 - 11,64 )*3,38 * 106
= 51,933166 млн. руб. /год
П1
= (189,7 – 162,7 - 11,67 )*3,38 * 106
= 51,818385 млн. руб. /год
DП = 51,933166 – 51,818385 = 0,114781 млн.руб./год
Таблица 3.1 Сводная таблица отказов
| Состав затрат | Сумма затрат, млн. руб. | Удельный вес, % | Себестоимость |
| 1. На подпитку тепловых сетей | 5,85 | 14,83 | 1,73 |
| 2. На перекачку теплоносителя | 10,95679 | 27,76 | 3,24 |
| 3. На зарплату | 5,32112 | 13,49 | 1,57 |
| 4. На амортизацию | 3,332 | 8,45 | 0,99 |
| 5. На ремонт оборудования | 7,981683 | 20,23 | 2,36 |
| 6. Экономический ущерб | 0,114783 | 0,29 | 0,03 |
| 7. На общецеховые расходы | 3,355638 | 8,52 | 0,99 |
| 8. На общесетевые расходы | 2,529601 | 6,41 | 0,75 |
| Итого | 39,441615 | 100 | 11,67 |
4. АНАЛИЗ ПОЛУЧЕННЫХ РЕЗУЛЬТАТОВ
В ходе проведения расчетов было получено, что в 2001 году намечается снижение количества отказов по повреждениям сальникового компенсатора и отказа арматуры. В других случаях отказов (свищей, разрывов трубы, разрывов сварного шва) прогнозируется увеличение числа отказов.
Таблица 4.1 Прогнозируемое количество отказов теплосети на 2000 год
| Причины отказов теплосети | Количество отказов |
| Свищ | 34 |
| Разрыв трубы | 14 |
| Разрыв сварного шва | 14 |
| Повреждение сальникового компенсатора | 10 |
| Отказ арматуры | 9 |
Таким образом, зная ожидаемое количество отказов по различным причинам, я рассчитал затраты на ремонт трубопроводов и другие издержки теплосети, а также себестоимость передачи тепловой энергии. В результате получили плановые себестоимости тепловой энергии с учетом отказов теплосети, и без учета отказов.
С пл
= 11,63 руб/Гкал
Cпл
от
= 11,67 руб/Гкал
Плановая себестоимость сети с учетом отказов сети превышает плановую себестоимость сети без учета отказов. Это означает снижение прибыли теплосети за счет отказов. Фактическая прибыль теплосети составляет П1
= 51,818385 млн. руб. /год, плановая - П2
= 51,818385 млн. руб. /год .Таким образом предприятие недополучает DП = 0,114781 млн.руб./год.
На мой взгляд, необходимо поднять тарифы на тепловою энергию, передаваемую теплосетью. Конечно, это ударит по карману потребителей тепловой энергии, но это позволит увеличить прибыль и, следовательно, качество трубопроводов, что в свою очередь приведет к улучшению качества обслуживания потребителей. Так же считается целесообразным увеличить ремонтный фонд, что поможет проводить ремонтные работы не только по мере надобности, но и в опережающем порядке.
СПИСОК ЛИТЕРАТУРЫ
1. Инструкция по планированию, учету и калькулированию себестоимости электрической и тепловой энергии в энергосистемах и на электростанциях, затрат на передачу и распределение энергии в электрических и тепловых сетях, -М.: Энергоатомиздат, 1970.-68с.
2. Четыркин Е.М. Статистические методы прогнозирования. Изд. 2-е перераб. и доп. -М.: Статистика, 1977. -200 с.
3. Юзбашев М.М., Манелля А.И. Статистический анализ тенденций и колеблемости. -М.: Финансы и статистика, 1983. -216 с.
4. Положение о составе затрат по производству и реализации продукции (работ, услуг), включаемых в себестоимость продукции (работ, услуг), и о порядке формирования финансовых результатов, учитываемых при налогообложении прибыли: Утв. постановлением Правительства РФ от 5 августа 1992 г. №55 // Экономика и жизнь, 1992. -№ЗЗ.
5. Инструкция по расследованию и учету технологических нарушений в работе электростанций, сетей и энергосистем: РД 34.20.801-93. -М.: ОРГРЭС, 1993. -58с.