В качестве экономической задачи в данной работе нужно рассчитать плановую себестоимость передачи тепловой энергии тепловыми сетями с учетом прогнозирования возможных причин отказов, оказывающих большое влияние на величину планируемых расходов (затрат) на ремонт трубопроводов теплосетей и потерь тепловой энергии с утечками. Характерными причинами отказов являются: разрыв трубы, разрыв сварного шва, свищ, отказ арматуры, повреждение сальникового компенсатора.
В качестве метода прогнозирования, возможных причин отказов тепловых сетей в работе используется метод прогнозирования на основе выявления временного тренда, путем аналитического выравнивания уровней динамического ряда.
Глава 1. ПРОГНОЗИРОВАНИЕ ОТКАЗОВ ТЕПЛОВЫХ СЕТЕЙ
Методика статистического прогноза по тренду основана на предположении, что параметры тренда и колебаний сохраняются до прогнозируемого периода.
Данная методика расчета включает в себя несколько видов:
· Выявление типа тренда;
· Вывод уравнения тренда;
· Анализ полученной зависимости;
· Прогноз на требуемый период.
1.1 Выявление типа тренда
В качестве динамического ряда берется изменение количества отказов тепловых сетей за предлагаемый период времени в зависимости от характера отказа. Эмпирические уровни ряда сглаживаются методом скользящей средней. Сглаженные уровни наносят на график.
Характер отказа – Разрыв трубы
Характер отказа – Повреждение сальникового компенсатора
Характер отказа – отказ арматуры
1.2 ВЫВОД УРАВНЕНИЯ ТРЕНДА
На графическом изображении сглаженных уровней отчетливо проявляется основная тенденция изменения количества отказов тепловых сетей по различным причинам.
Следующим этапом работы является получение обобщенной статистической оценки тренда методом аналитического выравнивания.
Основная тенденция развития yt как функция времени:
yti= f (ti).
Определение теоретических (расчетных) уровней yti производится на основе адекватной математической модели, которая наиболее точно отображает основную тенденцию ряда динамики – количество отказов теплосети по данной причине.
Подбор математической функции основан на определении типа развития экономического явления во времени.
Когда тип тренда установлен, необходимо вычислить оптимальные значения параметров тренда, исходя из фактических уровней, т.е. получить значения коэффициентов уравнения. На основе требований метода наименьших квадратов (МНК) составляется система нормальных уравнений, решая которую, вычисляем параметры тренда.
Для линейного тренда нормальные уравнения МКТ имеют вид:
n*a + b*St i= Sy ia*St i + b*St2i = S yi *ti
где yi – уровни исходного ряда динамики;
ti – номера периодов или моментов времени;
Нормальные уравнения МКТ для тренда в виде параболы 2-го порядка:
n*a + b*Sti + c*Sti2 = Syia*Sti + b*Sti2 + c*Sti3 = Syi*ti
a*Sti2 + b*Sti3 + c*Sti4 = Syi*ti2
где yi – уровни ряда динамики;
ti – номера периодов или моментов времени;
a, b, c – константы уравнений.
свищ | 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;
С/к | 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;
Отказ арматуры | 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 |
Для линейного тренда: