Прогноз ущерба при аварийных перерывах в подаче теплоэнергии
Аварии или технические отказы, приводящие к полному или частичному отключению теплоэнергии, могут вызвать простои производства, брак продукции, а также нарушение параметров воздушной среды, способное снизить работоспособность, производительность труда, повышение заболеваемости персонала.
Ущерб Y, руб/год, от аварий или отказов в денежном выражении за год составит
Y = аN τв + bN
а = Y зарп + Y зав + ( Ц – С)В; b = Y бр + Y об + Y техн, (16)
где – У зарп - зарплата работников при отсутствии выпуска продукции, руб/ч;
Узав – заводские затраты, связанные с управлением и обслуживанием производства при отсутствии выпуска продукции, руб/ч;
Ц – цена выпускаемой продукции, руб/ед.прод.;
С – себестоимость продукции, руб/ед.прод.;
В – производительность (выработка продукции) при нормальном функционировании производства, ед.прод./ч.;
(Ц – С)В – потеря прибыли от простоя, руб/ч;
Yбp, Уоб, Ymexн – ущерб от брака, поломки оборудования и от восстановления технологического процесса, руб/отказ;
N – число отказов в год, отк/год,;
τв – время восстановления каждого отказа (продолжительность каждого простоя производства при одном отказе), ч.
Число отказов (простоев) и длительность простоев являются случайными величинами со своими законами распределения. Поэтому ущерб также представляет собой случайную величину – функцию двух случайных аргументов: числа отказов и длительностей восстановления системы. Среднюю величину ущерба можно найти как функцию матожиданий отказов в год
и длительностей восстановлений системы Тв
(17)
гдe а, руб/ч; b, руб/отк – константы.
Полагаем, что распределение отказов подчиняется закону Пуассона, а длительность простоя (восстановления) – экспоненциальному. Величины и Тв можно определить по значениям интенсивностей отказов λcucm и восстановления системы μcucm. Однако, из статистики отказов эти параметры известны только для элементов, из которых состоят системы, подсистемы, блоки (т.е. известны λi и μi). Поэтому, требуется вначале составить структурно-логическую схему надежности системы теплоснабжения, а затем, пользуясь аппаратом теории надежности, подобрать соответствующие формулы расчета вычислить
и Тв системы.
Рассмотрим одну из типовых систем теплоснабжения npомпредприятия, состоящую из подсистем: котельной и теплотрассы. Соответствующая структурно-логическая схема показана на рис. 3.
Рис. 3. Структурно-логическая схема надежности системы теплоснабжения промпредприятия: А - элемент, узел (1 - котел, 2 - вентилятор, 3 - дымосос, 4 - экономайзер, 5 - трубопровод, 6 - задвижка); Б - агрегаты, блоки (I - рабочий котлоагрегат, II - резервный котлоагрегат, III, IV - блоки питательных насосов, V - блок химводоподготовки); В - подсистемы (котельная, теплотрасса)
По параметрам надежности элементов λi и μi проводим расчет от элементов к агрегату, от агрегата к подсистеме, от подсистемы к системе по формулам табл. 6, 7, в которых приняты обозначения: l – длина трубопровода, км; – интенсивность отказов трубопровода на 1 км длины; пз, – числ задвижек;
– интенсивность отказов задвижки; Твт – среднее время восстановления теплотрассы, ч; d – диаметр трубопровода, м; Твт – среднее время восстановления системы.
Определить величину ущерба при заданной вероятности ее превышения можно двумя способами: построением закон распределения ущерба как композиции законов двух случайных аргументов (методом Монте-Карло) или вычисление дисперсии распределения ущерба и, в предположении, что композиция распределена нормально, использованием таблицы нормального распределения.
Таблица 6
Формулы для определения интенсивности отказов λ
Котлоагрегат | λка = ∑ λi | (7) |
Блок котлоагрегатов из п рабочих котлов и отсутствий резервных | λк = n λка | (9) |
Блок котлоагрегатов из п рабочих котлов и т резервных (ненагруженный скользящий резерв с ограниченным восстановлением) | ![]() | (12) |
Теплотрасса | ![]() | (18) |
Система теплоснабжения | ![]() | (19) |
Матожидание числа отказов системы | ![]() | (20) |
Таблица 7
Формулы для определения интенсивности восстановления μ
Котлоагрегат | ![]() | (8) |
Блок котлоагрегатов при отсутствий резервных | ![]() ![]() | (10) |
Блок котлоагрегатов из п рабочих котлов и т резервных (ненагруженный скользящий резерв с ограниченным восстановлением) | ![]() ![]() | (13) |
Теплотрасса | ![]() ![]() | (21) ф-ла МИСИ |
Система теплоснабжения | ![]() ![]() | (22) |
Второй способ более простой. Дисперсия ущерба
(23)
а для конкретных законов распределения (Пуассона и экспоненциального), т.е. при и
(24)
Среднеквадратичное отклонение ущерба
(25)
Задаваясь вероятностью непревышения ущерба Ф(Y), найдем по табл. 9 значение t, представляющее число сигм, отстоящих от среднего значения ущерба Y. Тогда величина ущерба с заданной вероятностью его непревышения составит
(26)
Таблица 8
Варианты заданий к лабораторной работе № 2
№ Варианта | Завод | Котлы | Длинна теплотрассы км. | Количество задвижек шт. | |
Марка | Количество n | ||||
«ЖБИ-1» | НР-18 | 2,5 | |||
«Тесмостепс» | ДКВР 2,5-13 | 5,0 | |||
«Красные Баррикады» | ГВ 900 | 2,0 | |||
«им. Ленина» | ДКВР 2,5-13 | 5,0 | |||
«Ахтуба» | Хопер 100 | 7,0 | |||
«Аврора» | ДКВР 2,5-13 | 2,5 | |||
«Медоборудование» | ДЕ10-14 | 11,7 | |||
«ЖБИ-1» | НР-18 | 2,7 | |||
«Тесмостепс» | ДКВР 2,5-13 | ||||
«Красные Баррикады» | ГВ 900 | 4,8 | |||
«им. Ленина» | ДКВР 2,5-13 | ||||
«Ахтуба» | Хопер 100 | ||||
«Аврора» | ДКВР 2,5-13 | 12,1 | |||
«Медоборудование» | ДЕ10-14 | 2,5 | |||
«ЖБИ-1» | НР-18 | 6,6 | |||
«Тесмостепс» | ДКВР 2,5-13 | 2,0 | |||
«Красные Баррикады» | ГВ 900 | 3,0 | |||
«им. Ленина» | ДКВР 2,5-13 | 6,0 | |||
2 | |||||
«Ахтуба» | Хопер 100 | 3,5 | |||
«Аврора» | ДКВР 2,5-13 | 9,7 | |||
«Медоборудование» | ДЕ10-14 | 6,7 | |||
«ЖБИ-1» | НР-18 | 3,3 | |||
«Тесмостепс» | ДКВР 2,5-13 | 4,8 | |||
«Красные Баррикады» | ГВ 900 | ||||
«им. Ленина» | ДКВР 2,5-13 | ||||
«Ахтуба» | Хопер 100 | 8,1 | |||
«Аврора» | ДКВР 2,5-13 | 3,7 | |||
«Медоборудование» | ДЕ10-14 | 2,5 | |||
«ЖБИ-1» | ДКВР 2,5-13 | 4,8 | |||
«Тесмостепс» | Хопер 100 | ||||
«Красные Баррикады» | ДКВР 2,5-13 | ||||
«им. Ленина» | ДЕ10-14 | 8,1 |
Таблица 9