Статьи
English version
О компании
История
Руководство
Сотрудники
Разработки
Объекты
Наши публикации
Наши патенты
Контакты
Обратная связь


Метод прогнозирования срока службы алюминиевых электролизеров

В.А Козлов, С.Н.Ахмедов, Б.С.Громов, Р.В Пак

Прогнозирование срока службы электролизеров имеет большое значение в практике работы заводов по производству алюминия, так как позволяет оценить эффективность предпринимаемых усовершенствований в конструкции электролизера и технологии его эксплуатации, а также делать более обоснованный выбор величины объема капитального ремонта на планируемый период. Для прогнозирования срока службы электролизеров можно применять нормальное, логарифмическое, биноминальное распределения. Однако на практике наибольшее распространение получило применение трехпараметрического распределения Вейбулла [1,2,3]. Преимущество применения функции распределения трёх параметров Вэйбулла состоит в получении качественной возможности прогнозирования выхода из строя, средних сроков службы, медианы отказов на базе ограниченного числа ранних выходов из строя электролизеров данной группы. Другое существенное преимущество распределения Вейбулла состоит в его гибкости, т.е. возможности моделировать практически любые функции риска.

Нарастающая функция распределения Вейбулла выражается как:

(1)

где F(t) - часть электролизеров, выходящая из строя за время t

 

α - масштабный множитель

β - коэффициент формы

Θ - параметр положения

Плотность функции распределения выражается формулой:

(2)

Влияние масштабного множителя α аналогично влиянию стандартного отклонения нормального распределения, т.е. с увеличением α происходит выравнивание распределения. Коэффициент β влияет на асимметрию распределения. Когда величина этого коэффициента находится в пределах 3-4, что чаще всего наблюдается при аппроксимации данных по сроку службы электролизеров, распределение плотности функции близко к симметричному. При уменьшении коэффициента наблюдается перекос распределения в сторону ранних поломок, при увеличении коэффициента имеется перекос к электролизерам с длительным сроком службы. Параметр положения Θ позволяет сдвинуть распределение от нулевого значения на некоторую величину, т.е. при этом подразумевается, что спустя некоторое время после ввода в эксплуатацию вероятность выхода из строя практически равна нулю.

Рис. 1. Вид функции распределения P(t) и плотности распределения f(t) Вейбулла

Пример функции распределения, плотности распределения, для распределения Вейбулла с масштабным множителем 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



Рис. 2. Параметры аппроксимационной кривой зависимости ln(ln1/R) от ln(t- Θ) при Θ = 150

На основании полученных данных, исходя из условия максимума значения R2, выбираем значение Θ = 150. Уравнение аппроксимационной кривой (рис.2):

y = 2,8148x - 20,691

Отсюда:

β = 2,8148

α = ехр (20,691/2,8148) = 1557

Рис. 3. Фактические данные по отключению электролизеров и расчетная кривая функции Вейбулла

График полученной функции распределения с фактическими точками и рассчитанной кривой представлен на рис. 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.


Реферат

Рассмотрено использование трехпараметрического распределения Вейбулла для использования с целью прогнозирования среднего срока службы электролизеров для получения алюминия. Приведены расчетные формулы для расчета прогнозируемого срока службы, моды и медианы отключений. Показано хорошее согласование результатов расчета и фактических данных по выходу электролизеров из строя. Использование метода расчета на компьютерной программе позволит выполнить прогнозы по сроку службы электролизеров различных типов и модификаций, определяя эффективные решения в конструкции электролизеров и технологии их эксплуатации, и рассчитать обоснованные объемы капитального ремонта электролизеров.




designed by

Sagitta


ALCORUS ENGINEERING Ltd,
Copyright ©  2004-2019