|
|
Метод прогнозирования срока службы алюминиевых электролизеров
В.А Козлов, С.Н.Ахмедов, Б.С.Громов, Р.В Пак
Прогнозирование срока службы электролизеров имеет большое значение в практике работы заводов по производству алюминия, так как позволяет оценить эффективность предпринимаемых усовершенствований в конструкции электролизера и технологии его эксплуатации, а также делать более обоснованный выбор величины объема капитального ремонта на планируемый период. Для прогнозирования срока службы электролизеров можно применять нормальное, логарифмическое, биноминальное распределения. Однако на практике наибольшее распространение получило применение трехпараметрического распределения Вейбулла [1,2,3]. Преимущество применения функции распределения трёх параметров Вэйбулла состоит в получении качественной возможности прогнозирования выхода из строя, средних сроков службы, медианы отказов на базе ограниченного числа ранних выходов из строя электролизеров данной группы. Другое существенное преимущество распределения Вейбулла состоит в его гибкости, т.е. возможности моделировать практически любые функции риска.
Нарастающая функция распределения Вейбулла выражается как:
|
(1) |
где F(t) - часть электролизеров, выходящая из строя за время t
|
α - масштабный множитель
β - коэффициент формы
Θ - параметр положения |
Плотность функции распределения выражается формулой:
|
(2) |
Влияние масштабного множителя α аналогично влиянию стандартного отклонения нормального распределения, т.е. с увеличением α происходит выравнивание распределения. Коэффициент β влияет на асимметрию распределения. Когда величина этого коэффициента находится в пределах 3-4, что чаще всего наблюдается при аппроксимации данных по сроку службы электролизеров, распределение плотности функции близко к симметричному. При уменьшении коэффициента наблюдается перекос распределения в сторону ранних поломок, при увеличении коэффициента имеется перекос к электролизерам с длительным сроком службы. Параметр положения Θ позволяет сдвинуть распределение от нулевого значения на некоторую величину, т.е. при этом подразумевается, что спустя некоторое время после ввода в эксплуатацию вероятность выхода из строя практически равна нулю.
Пример функции распределения, плотности распределения, для распределения Вейбулла с масштабным множителем 1600, коэффициентом формы 2,5 и параметром положения 40 показан на риc.1.
Функция надежности, отражающая количество остающихся в работе электролизеров через время t вычисляется по формуле:
|
(3) |
Математическое ожидание (среднее значение срока службы) распределения Вейбулла можно вычислить по формуле:
|
(4) |
где:
|
(5) |
Г(х) - гамма функция, значение которой табулировано. Для упрощения расчетов значения гамма-функции в диапазоне значений от 2 до 6 можно аппроксимировать с точностью R2=0,9998 полиномом второй степени следующего вида:
|
(6) |
Величина моды (при условии β>1) подсчитывается из уравнения:
|
(7) |
Величина медианы определяется из условия равенства функции распределения 0,5 и рассчитывается в соответствии с уравнением:
|
(8) |
Как уже отмечалось выше, в диапазоне значений коэффициента β между 3 и 5 функция распределения приближена к симметричному виду, т.е. величины математического ожидания близки к значению моды и практически совпадают со значением медианы. Разница в значении срока службы и медианы в указанном диапазоне β не превышает 1%. Отсюда следует, что время с момента пуска до момента, при наступлении которого будет отключено 50% электролизеров исследуемой группы, практически равно ожидаемому среднему сроку службы всех электролизеров.
Оценку параметров распределения с удовлетворительной точностью можно производить при выходе из строя не менее 10% электролизеров исследуемой группы /2/. При большом объеме выборок наилучшие результаты показывает применение метода максимального правдоподобия. Однако данный метод требует применения сложного математического аппарата и не всегда обеспечивает точные результаты при малом объеме выборки (n<15). Для практики алюминиевого завода более простым и достаточно точным методом определения параметров распределения может служить графический метод. После преобразования и логарифмирования уравнения функции надежности можно получить следующее выражение:
|
(9) |
Если на графике по оси абсцисс отложить значение ln(t-Θ), по оси ординат - ln(ln(1/R(t))) и аппроксимировать полученную зависимость прямой линией, то угол наклона будет соответствовать значению β, а значение α можно будет подсчитать, разделив значение свободного члена уравнения на β и взяв экспоненту полученного значения.
Наибольшую трудность вызывает определение значения коэффициента положения Θ. Определить его величину можно итерационным методом, последовательно наращивая величину ? от нуля на некоторое значение и рассчитывая аппроксимационные прямые зависимости ln(ln1/R) от ln(t- ?). Для каждой полученной прямой необходимо вычислить значение квадрата коэффициента корреляции Пирсона R2 и выбрать величину ?, при которой значение R2 проходит через максимум. На неправильно выбранную величину ? косвенно указывает слишком большая величина параметра формы ? (>6), что вызывает чрезмерный перекос в сторону ванн с малым сроком службы.
Проведение статистических расчетов срока службы требует постоянства условий эксплуатации анализируемой группы электролизеров во времени. Изменение технологических параметров в процессе эксплуатации исследуемой группы (увеличение силы тока, изменение состава электролита, условий питания ванн глиноземом и т.п.) может существенно изменить динамику выхода электролизеров из строя. Также необходим тщательный анализ выхода из строя первых электролизеров исследуемой группы. Механизм их разрушений может отличаться от такового для остальных ванн группы. Причиной раннего отключения могут быть брак по материалам футеровки, грубые нарушения технологии пуска и эксплуатации, и другие воздействия, которым не были подвержены остальные электролизеры группы. Поэтому данные о досрочном отключении подобных электролизеров следует исключить из анализа.
Ниже приведен пример расчета ожидаемого среднего срока службы группы из 75 электролизеров, после момента выхода из строя 17 штук или 23% от общего количества. Необходимые исходные данные приведены в таблице 1.
Исходные данные для расчета ожидаемого срока службы
Время с момента пуска, t
| Действует электролизеров
| Остановлено электролизеров
|
дни
| шт.
| % от пущенных, P
| шт.
| % от пущенных, R
|
167 |
74 |
98,7% |
1 |
1,3% |
590 |
73 |
97,3% |
2 |
2,7% |
645 |
72 |
96,0% |
3 |
4,0% |
688 |
71 |
94,7% |
4 |
5,3% |
745 |
70 |
93,3% |
5 |
6,7% |
806 |
69 |
92,0% |
6 |
8,0% |
850 |
68 |
90,7% |
7 |
9,3% |
867 |
67 |
89,3% |
8 |
10,7% |
925 |
66 |
88,0% |
9 |
12,0% |
930 |
65 |
86,7% |
10 |
13,3% |
941 |
64 |
85,3% |
11 |
14,7% |
988 |
63 |
84,0% |
12 |
16,0% |
1010 |
62 |
82,7% |
13 |
17,3% |
1040 |
61 |
81,3% |
14 |
18,7% |
1049 |
60 |
80,0% |
15 |
20,0% |
1080 |
59 |
78,7% |
16 |
21,3% |
1120 |
58 |
77,3% |
17 |
22,7% |
|
Из представленных данных для проведения дальнейших расчетов исключаем первую запись, т.к. срок службы первого отключенного электролизера несопоставимо меньше сроков отключения других электролизеров, и можно предположить, что механизм его выхода из строя существенно отличается от других электролизеров исследуемой группы. Последовательно изменяем значение коэффициента Θ и строим таблицу зависимости величины ln(ln1/R) от ln(t- Θ) (табл. 2). Находим для каждой зависимости методом наименьших квадратов коэффициенты аппроксимационной прямой и значение квадрата коэффициента корреляции Пирсона R2.
Расчет параметров распределения
ln(t - Θ)
| ln(ln1/R)
|
Θ = 0
| Θ = 50
| Θ = 100
| Θ = 150
| Θ = 200
| Θ = 250
|
3,611 |
6,3801 |
6,2916 |
6,1944 |
6,0868 |
5,9661 |
5,8289 |
3,199 |
6,4693 |
6,3886 |
6,3008 |
6,2046 |
6,0981 |
5,9789 |
2,904 |
6,5338 |
6,4583 |
6,3767 |
6,2879 |
6,1903 |
6,0822 |
2,674 |
6,6134 |
6,5439 |
6,4693 |
6,3886 |
6,3008 |
6,2046 |
2,484 |
6,6921 |
6,628 |
6,5596 |
6,4862 |
6,4069 |
6,3208 |
2,323 |
6,7452 |
6,6846 |
6,6201 |
6,5511 |
6,477 |
6,3969 |
2,182 |
6,765 |
6,7056 |
6,6425 |
6,5751 |
6,5028 |
6,4249 |
2,057 |
6,8298 |
6,7742 |
6,7154 |
6,6529 |
6,5862 |
6,5147 |
1,944 |
6,8352 |
6,7799 |
6,7214 |
6,6593 |
6,593 |
6,5221 |
1,841 |
6,8469 |
6,7923 |
6,7346 |
6,6733 |
6,608 |
6,5381 |
1,747 |
6,8957 |
6,8437 |
6,789 |
6,731 |
6,6695 |
6,6039 |
1,659 |
6,9177 |
6,8669 |
6,8134 |
6,7569 |
6,697 |
6,6333 |
1,577 |
6,947 |
6,8977 |
6,8459 |
6,7912 |
6,7334 |
6,672 |
1,5 |
6,9556 |
6,9068 |
6,8554 |
6,8013 |
6,7441 |
6,6834 |
1,427 |
6,9847 |
6,9373 |
6,8876 |
6,8352 |
6,7799 |
6,7214 |
1,359 |
7,0211 |
6,9754 |
6,9276 |
6,8773 |
6,8244 |
6,7685 |
коэф. α |
3,4562 |
3,2435 |
3,0298 |
2,8148 |
2,5982 |
2,3795 |
коэф. β |
25,578 |
23,942 |
22,314 |
20,691 |
19,073 |
17,458 |
R2 |
0,9941 |
0,9943 |
0,9944 |
0,9945 |
0,9944 |
0,9941 |
|
На основании полученных данных, исходя из условия максимума значения R2, выбираем значение Θ = 150. Уравнение аппроксимационной кривой (рис.2):
y = 2,8148x - 20,691
Отсюда:
β = 2,8148
α = ехр (20,691/2,8148) = 1557
График полученной функции распределения с фактическими точками и рассчитанной кривой представлен на рис. 3.
Подставляя коэффициенты в уравнения 4-8, получаем, что значение моды для полученного распределения равно 1483 (сут.), значение медианы - 1517 (сут.), средний срок службы - 1537 (сут.).
Как показывает опыт авторов, точность прогноза срока службы электролизеров при выходе из строя около 20% от общего количества пущенных в анализируемой группе составляет не более 15%.
Все вышеприведенные расчеты реализованы ООО "АЛКОРУС ИНЖИНИРИНГ" в виде расчетной программы обработки данных по отключению электролизеров. Использование программы на алюминиевых заводах позволит выполнить прогнозы по сроку службы электролизеров различных типов и модификаций, определяя эффективные решения в конструкции электролизеров и технологии их эксплуатации, и рассчитать обоснованные объемы капитального ремонта электролизеров.
Библиографический список
|
1. C.A. Brown // Light Metals. 1985. P. 705-710.
2. A. Weatherly // Light Metals. 1985. P. 711-717.
3. M.A. Doheim, S.M. El-Raghy, M.A. Mohamed // Light Metals. 2001. P. 519-526.
|
Реферат
Рассмотрено использование трехпараметрического распределения Вейбулла для использования с целью прогнозирования среднего срока службы электролизеров для получения алюминия. Приведены расчетные формулы для расчета прогнозируемого срока службы, моды и медианы отключений. Показано хорошее согласование результатов расчета и фактических данных по выходу электролизеров из строя. Использование метода расчета на компьютерной программе позволит выполнить прогнозы по сроку службы электролизеров различных типов и модификаций, определяя эффективные решения в конструкции электролизеров и технологии их эксплуатации, и рассчитать обоснованные объемы капитального ремонта электролизеров.
|
|