Таблица 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;
Анализ полученной зависимости
Одним из применяемых на практике показателей адекватности математической функции для статистического изучения тренда является стандартизированная ошибка аппроксимации - 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,188 |
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. ОПРЕДЕЛЕНИЕ ИЗДЕРЖЕК ПРОИЗВОДСТВА ПРИ ПЕРЕДАЧЕ ТЕПЛОВОЙ ЭНЕРГИИ
Годовые издержки производства (эксплуатационные расходы), связанные с передачей тепловой энергии, включают в себя следующие статьи:
- затраты на подпитку тепловых сетей;
- затраты на перекачку теплоносителя;
- фонд оплаты труда обслуживающего персонала;
- амортизационные отчисления на полное восстановление (реновацию) основных производственных фондов теплосети;
- затраты на текущие и капитальные ремонты основного и вспомогательного оборудования тепловых сетей;
- экономический ущерб, связанный с отказами тепловых сетей;
- общецеховые расходы;
- общесетевые расходы.
Расчет фонда амортизации
По этой статье в себестоимости передачи тепловой энергии учитываются амортизационные отчисления на полное восстановление ( реновацию ) основного и вспомогательного оборудования, а также отчисления на реновацию и капитальный ремонт зданий и сооружений.
Сумма амортизационных отчислений по основным фондам теплосети,
Иам = (Кбал * аср)/100
Где Кбал - общая балансовая стоимость оборудования, зданий и т.п. теплосети;
аср - средняя по основным фондам норма амортизации на полное восстановление (реновацию ), %.
Кбал= 83,3 млн руб/год; аср = 4 %;
Иам = (83,3* 4)/100 = 3,332 млн руб/год;
Общецеховые расходы
К этой статье относятся затраты по обслуживанию и управлению тепловыми сетями, то есть зарплата аппарата управления, затраты по текущему ремонту зданий, расходы по охране труда и т.д.
Цеховые расходы могут быть рассчитаны укрупненно, млн руб/год:
Ицех= g*Иэкс.
Где g - доля общецеховых расходов от суммарных эксплуатационных.
g = 0,1;
Ицех= 0,1 * 33,556378 = 3,3556378 млн руб/год
Общесетевые расходы
В этой статье учитываются расходы по управлению предприятием тепловых сетей:
- административно-управленческие расходы.
- общепроизводственные расходы (НИОКР, текущий ремонт общесетевых средств и т.д.).
Исет = ЗП адм * 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 |
СПИСОК ЛИТЕРАТУРЫ
1. Инструкция по планированию, учету и калькулированию себестоимости электрической и тепловой энергии в энергосистемах и на электростанциях, затрат на передачу и распределение энергии в электрических и тепловых сетях, -М.: Энергоатомиздат, 1970.-68с.
2. Четыркин Е.М. Статистические методы прогнозирования. Изд. 2-е перераб. и доп. -М.: Статистика, 1977. -200 с.
3. Юзбашев М.М., Манелля А.И. Статистический анализ тенденций и колеблемости. -М.: Финансы и статистика, 1983. -216 с.
4. Положение о составе затрат по производству и реализации продукции (работ, услуг), включаемых в себестоимость продукции (работ, услуг), и о порядке формирования финансовых результатов, учитываемых при налогообложении прибыли: Утв. постановлением Правительства РФ от 5 августа 1992 г. №55 // Экономика и жизнь, 1992. -№ЗЗ.
5. Инструкция по расследованию и учету технологических нарушений в работе электростанций, сетей и энергосистем: РД 34.20.801-93. -М.: ОРГРЭС, 1993. -58с.
Таблица 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;
Анализ полученной зависимости
Одним из применяемых на практике показателей адекватности математической функции для статистического изучения тренда является стандартизированная ошибка аппроксимации - 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,188 |
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. ОПРЕДЕЛЕНИЕ ИЗДЕРЖЕК ПРОИЗВОДСТВА ПРИ ПЕРЕДАЧЕ ТЕПЛОВОЙ ЭНЕРГИИ
Годовые издержки производства (эксплуатационные расходы), связанные с передачей тепловой энергии, включают в себя следующие статьи:
- затраты на подпитку тепловых сетей;
- затраты на перекачку теплоносителя;
- фонд оплаты труда обслуживающего персонала;
- амортизационные отчисления на полное восстановление (реновацию) основных производственных фондов теплосети;
- затраты на текущие и капитальные ремонты основного и вспомогательного оборудования тепловых сетей;
- экономический ущерб, связанный с отказами тепловых сетей;
- общецеховые расходы;
- общесетевые расходы.
Дата: 2019-12-22, просмотров: 237.