CC BY
24
УДК 621.18.021
НЕКОТОРЫЕ РЕЗУЛЬТАТЫ ПРОГНОЗИРОВАНИЯ РЕСУРСА И НАДЕЖНОСТИ ТЕПЛООБМЕННЫХ АППАРАТОВ ТЭС
ДЕРИЙ В.П., инж.
Приведены результаты прогноза ресурса трубчатки сетевых подогревателей ТЭС. Предложена математическая модель оценки надежности теплообменных аппаратов, позволяющая учесть зависимость интенсивности восстановления оборудования от времени ремонта.
Ключевые слова: теплообменные аппараты, надежность, математическая модель, работоспособность системы.
SOME FORECASTING RESULTS OF TPS HEAT-EXCHANGE APPARATUS
RELIABILITY AND RESOURCE
V.P. DERIY, engineer
The work represents the forecasting results of TPS (thermoelectric power station) networked heater pipe heater resource. The author has suggested mathematical model of heat-exchange apparatus reliability evaluation, which allows time dependence of equipment reconstruction intensity to be taken into account.
Key words: heat-exchange apparatus, reliability, mathematical model, system efficiency.
На основе разработанной нами ранее математической модели [1] по данным обследования выполнен прогноз ресурса трубчатки сетевых подогревателей ТЭС. Среднее число заглушенных трубок и дисперсия распределения определялись по формулам < N >= 1- (1-< N0 >)ехр(-д(!)), (1)
g(t) = jx(t)dt = X0t +
at2 ~2
Np
(1- < N >)2
(1-< N0 >)2
(1-< N0 >)(N - N0)
(2)
(3)
(1- < N >)
Здесь Ыр - полное количество трубок аппарата.
В качестве примера (рис. 1) приведены функции прогноза числа заглушенных трубок для двух сетевых подогревателей ИВТЭЦ-2.
Nl(t)
50
1 1 1 \> }
/ /
/ /
7 /
/ /
/ /
/ /
/ /
- / /
/ + /
у
+
+ 1 1 1
Рис. 1. Функции прогноза заглушенных трубок сетевых подогревателей:------ БО-350М;--БО-200
Анализ результатов расчета коэффициентов идентификации функций прогноза по ИВТЭЦ-2 (табл. 1) показал, что параметры
идентификации функций прогноза для отдельных подогревателей значительно отличаются друг от друга, что определяется водно-химическим режимом, конструкцией подогревателей и материалом трубчатки (в третьем столбце указано время достижения аппаратом критического состояния, а в последнем - относительная флуктуация ).
Таблица 1. Параметры функции прогноза
Название МО3 a-103 T N(T) 5
1/год 1/год2 годы - -
БО-2Б 4,9 2,01 5,7 101 0,1
БСВ-1 10 1,22 5,5 130 0.09
БО-5А -7,2 4,3 6 130 0.09
БП-1 -2,34 5,5 4 101 0,1
БП-3 4,3 2,5 5 101 0,1
Условия эксплуатации теплообменного оборудования старых станций таковы, что даже в течение одного года приходится неоднократно выводить в ремонт подогреватели для глушения трубчатки в больших объемах. Это обстоятельство с учетом низкого срока службы исследуемой аппаратуры делает малоэффективным прогноз состояния трубчатки на основе математического моделирования. Математическое моделирование может рассматриваться как дополнение к основному методу, который основывается на измерении фактического недогрева сетевой воды против расчетной величины.
Прогноз состояния трубчатки на основе математического моделирования целесообразен для новых станций с хорошими условиями эксплуатации, когда накопление дефектов происходит медленно. В качестве примера (рис. 2) приведены данные по распределению
0
затушенных трубок за десять лет работы подогревателя сетевой воды ПСВ-500-3-23. Прогнозируемое время достижения трубчаткой критического состояния равно 35 годам.
30
20
N(1)
10
О 5 10
I
Рис. 2. Функция прогноза состояния трубчатки ПСВ-500-3
Для старых и стареющих тепловых станций на первое место встают вопросы оценки надежности аппаратуры. К основным видам работ, проводимых на станциях для восстановления работоспособности теплообменных аппаратов, относятся чистка, заглушка и перепаковка трубчатки. В настоящее время информация по учету различных дефектов в теплоэнергетическом оборудовании паротурбинных установок отсутствует. В таких условиях нельзя выделить роль отдельных составляющих, определяющих надежность аппарата, и можно говорить только об их интегральном действии. Работа теплотехнического оборудования зависит от многих факторов: характера тепловых и динамических режимов отдельных аппаратов, конструктивных особенностей, качества материалов, их усталости, длительности срока службы, действия автоматики и пр. Многие характеристики системы можно считать определенными, детерминированными, но отказы и восстановления - это случайные состояния. Срок эксплуатации системы, определяющийся этими состояниями, тоже будет случайной величиной. Взаимные связи между отказами, восстановлениями, действием противоаварий-ной автоматики реализуют некий стохастический процесс, в результате которого система переходит из рабочего состояния в состояние ремонта с тем или иным сочетанием отказавших элементов. В таких условиях можно говорить о вероятности того, что аппарат в определенный момент времени может находиться в том или ином состоянии. Восстановление работоспособности системы осуществляется средствами противоаварийной автоматики, оперативными переключениями, выполняемыми обслуживающим персоналом и проведением аварийно-восстановительных работ. В первом приближении по своему характеру отказы можно разделить на внезапные и постепенные. Последние возникают при постепенном изме-
нении параметров (старение и износ). Основными задачами теории надежности являются:
• исследование возможности прогноза постепенных отказов оборудования;
• разработка методов проведения испытаний на надежность и обработка полученных результатов;
• выбор оптимальных режимов работы оборудования в ходе холодной и горячей обкатки и в режиме длительной эксплуатации.
При этом прогноз и стратегия эксплуатации должны базироваться как на результатах постоянного контроля состояния оборудования, так и на результатах математического моделирования физико-химических процессов, протекающих в системе и определяющих ее надежность. Таким образом, теория надежности является существенным составным элементом общей концепции по повышению экономичности паротурбинной установки и продления ее ресурса.
Обычно анализ надежности различных аппаратов ведется при одном существенно упрощающем допущении: функция распределения времени ремонта отказавшего оборудования подчиняется экспоненциальному закону [2]. Это равносильно тому, что интенсивность ввода отказавших аппаратов не зависит от длительности их ремонта. Получается, что вероятность ввода отказавшего аппарата непосредственно после выхода его из строя равна вероятности его ввода в любой момент времени в течение всего ремонта. Такое допущение практически никогда не выполняется и вводится оно по причине больших математических трудностей, возникающих при решении задачи. Однако если ввести плотность вероятности состояния аппарата по времени ремонта, то от указанного выше допущения можно отказаться. Далее считаем, что интенсивность восстановления аппарата является произвольной функцией времени восстановления ц = ц(т). Кроме интенсивности ц(т) процесс характеризуется вероятностью отказа аппарата в единицу времени X, вероятностью Р нахождения аппарата в момент времени 1 в рабочем состоянии в ожидании исходного случая и плотностью вероятности по времени ремонта р(т, 1).
Выведем уравнение для вероятности Р(1). Вероятность того, что аппарат к моменту времени находится в рабочем состоянии, складывается из вероятности того, что он находился в рабочем состоянии в момент времени 1 и за время ^ не вышел в ремонт, и вероятности того, что он в течение времени 1 был в ремонте, а за время ^ вышел из ремонта:
Р^ + сК) = Р(1)(1 - ХсК) + сК | ц( т')р( т', 1)с1т'.
0
Из этого следует уравнение ССР = -ХР(1) + } ц(т')р(т',1)с1т'. (4)
Аналогичным образом можно получить уравнение для плотности вероятности нахождения аппарата в ремонте:
ФСМ) ФСМ) , . , ..
„ =--Г---|(ФСМ) (5)
01 от
с граничным условием р(т = 0,1) = <Р(1).
В стационарном режиме при 1 ^ да сро) = 0 др(т,1) = 0 сК , д .
Решение уравнения (5) в стационарном случае имеет вид
р(т) = С exp(-j^(x')dx').
(6)
Подставляя это решение в уравнение (4), с учетом стационарности получим С=<Р. Постоянную Р = Р(1 ^ да) найдем из условия нормировки вероятностей
да
Р + }р(т)с1т = 1. (7)
о
Последний интеграл берется по частям:
=< т > .
(8)
} ехр(-}|(т')с1т') = }тс1Р 0 0 о
т
Здесь Р = 1 - ехр(-}|(т')ст' - функция распре-
о
деления времени ремонта.
Интеграл (8) имеет простой физический смысл - он равен среднему времени ремонта. Подставляя (8) в (7), находим
Р =-1-. (9)
1+ < <т>
Эта величина определяет стационарную вероятность того, что аппарат находится в работе. В теории надежности эту величину называют вероятностью безотказной работы.
Вероятность Р-| того, что аппарат находится в ремонте, определяется соотношением
да <
Р1 = }р(т)с1т = 1 -Р = <<> . (10)
•I 1 + < <т >
Из формул (9) и (10) видно, что искомые вероятности определяются только средним временем ремонта агрегата и не зависят от вида функции распределения времени ремонта. Из этого следует, что обычное допущение о постоянстве ц приводит к правильным результатам относительно надежности аппарата, если среднее время восстановления аппарата определяется экспериментально.
Для расчета указанных вероятностей необходимо знание величин < и <т>. Определение этих величин осуществлено усреднением большого количества данных как по новым, так и по стареющим станциям методом экспертных оценок [3]. Результаты расчета вероятностей Р и Р1 по значениям указанных величин (табл. 2) показывают, что все теплообмен-ное оборудование обладает практически одинаковой высокой вероятностью безотказной работы, тогда как частота повреждаемости и время нахождения оборудования в ремонте оказываются разными.
Таблица 2. Параметры надежности оборудования
Оборудование % <т> P P1
1/год час - -
Конденсаторы 3,0±0,44 6,6±1,51 0,998 0,002
ПВД 2,3±0,43 45,7±9,77 0,988 0,012
ПНД 2,1 ±0,34 16,9±3,24 0,996 0,004
ПСГ 1,1 ±0,20 21,3±3,23 0,997 0,003
ПСП 1,1 ±0,26 18,0±4,2 0,998 0,002
Список литературы
1. К вопросу прогнозирования ресурса теплоэнергетического оборудования тепловых и атомных электрических станций / В.К. Семенов, В.С. Щебнев, В.П. Дерий, В.Ф. Степанов // Вестник ИГЭУ. - Вып. 2. - 2007. - С. 30-33.
2. Гнеденко Б.В., Беляев Ю.К., Соловьев А.Д. Математические методы в теории надежности. - М.: Наука, 1965.
3. Бродов Ю.М., Плотников П.Н. Надежность кожу-хотрубных теплообменных аппаратов паротурбинных установок. - Екатеринбург: УПИ, 2001.
О
да
Дерий Владимир Петрович, ОАО «Атомтехэнерго» (г. Мытищи), заместитель главного инженера, телефон (4932) 26-99-17, e-mail: [email protected]
