CC BY
84
Власов А.Б., Буев С.А. Оценка показателей надёжности...
УДК 629.123 : 681.518.5
Оценка показателей надёжности корабельного электрооборудования, их влияние на безопасность мореплавания и вероятность страховых рисков
А.Б. Власов, С.А. Буев
Морская академия МГТУ, кафедра электрооборудования судов
Аннотация. Представлен расчёт надёжности электрооборудования с использованием данных тепловизионной диагностики различных морских судов в разные периоды времени. От надёжности оборудования во многом зависят показатели эффективности и безопасности использования судна. Применение систем технической диагностики позволяет предотвращать или существенно сокращать число отказов судового оборудования, вынужденных остановок в море, а следовательно, и число страховых случаев.
Abstract. The authors present the calculation of electrical equipment reliability using the data of IR-thermography diagnostic of various vessels in different dates. The indicators of ships operation efficiency and safety depend essentially on reliability. Application of technical diagnostic's systems allows preventing or substantially reducing the number of ship equipment failures, forced stop at sea, and therefore the number of insurance cases.
Ключевые слова: надёжность, тепловизионная диагностика, морские суда, электрооборудование, страховые случаи, эффективность, безопасность, отказы оборудования
Key words: reliability, IR-diagnostics, vessels, electrical equipment, insurance cases, efficiency, safety, failures
1. Введение
Безопасная эксплуатация водного транспорта обеспечивается безотказной работой судового оборудования. Безотказность характеризует способность объекта (изделия, системы) сохранять работоспособность в течение определённого времени (наработки) при нормальных условиях эксплуатации. Количественными показателями безотказности являются вероятность безотказной работы, средняя наработка до отказа, интенсивность отказов и др. (Рябинин, Парфёнов, 1997).
Тепловизионная диагностика была разработана для системного сопровождения объектов энергетики и доказала свою эффективность (Власов, 2002). Указанная диагностика может применяться на морских судах при условии, что результаты диагностики будут удовлетворять требованиям Российского морского регистра судоходства для освидетельствования судов в эксплуатации и страховых компаний для принятия морских судов на страхование, оценки риска морского страхования.
Несмотря на широкое распространение тепловизоров в электроэнергетике, проблемами остаются оценка достаточности тепловизионного метода как средства определения эксплуатационных параметров надёжности электрооборудования, применимость методов статистической обработки данных, полученных с помощью тепловизионного контроля изделий в процессе их эксплуатации (Власов, Джура, 2002).
2. Статистические методы определения эксплуатационных показателей надёжности электрооборудования морских судов по данным тепловизионной диагностики
Рассмотрим методику статистического анализа данных, полученных в ходе инфракрасной диагностики оборудования морских судов, с целью определения эксплуатационной надёжности контактных соединений различных главных распределительных щитов (ГРЩ).
Статистический параметр потока повреждений т (т) ("точечная оценка" параметра) рассчитывается по соотношению (Власов, 2006):
ю*(т) = п0(т) / АД(т), год-1, (1)
где п0(т) - количество повреждений, накопленных с момента предыдущих испытаний; Д(т) - время, прошедшее между испытаниями, год; N - общее количество элементов (контактных соединений), находящихся под наблюдением на судне.
Статистическая вероятность появления повреждения Q определяется выражением:
Q (т, т + Дт) = п0(т + Дт) / N, (2)
672
Вестник МГТУ, том 16, № 4, 2013 г.
стр.672-680
где п0(т + Ат) - число повреждений, наблюдаемых в интервале (0, т + Ат). При фиксированном потоке повреждений т(т)
т(т) = то. (3)
Функция вероятности безотказной работы Р(т) восстанавливаемых элементов для экспоненциального закона распределения определяется как
Р(т) = exp (-тот). (4)
Статистическая вероятность безотказной работы R (т) рассчитывается как
R (т, т + Ат) = 1 - Q (т, т + Ат). (5)
В табл. 1 приведены статистические данные, полученные в процессе проведения тепловизионной диагностики контактных соединений (КС) на различных морских судах.
Таблица 1. Сведения по контактным соединениям на обследованных судах
Параметр Объект
I II III
Уровень напряжения, В 380 380 380
Общее количество КС на ГРЩ 402 400 430
Верхняя граница вероятности повреждений QB 0,007462 0,0175 0,004651
Нижняя граница вероятности безотказной работы RH 0,992537 0,9825 0,995349
Средний параметр потока повреждений тср, год-1 0,022 0,023 0,007
Примечание: I - транспортно-буксирное судно "Десна", ООО "Газфлот"; II - самоходный плавучий кран "СПК-395", ЗАО "Мурманская судоремонтная компания"; III - рыболовецкий траулер "Альферас", ЗАО "Мурмансельдь-2".
На рисунке приведена зависимость вероятности безотказной работы R (т) контактных соединений на объекте I.
Рис. График R*(6) контактных соединений "СПК-395"
При анализе характеристик надёжности рассчитывается вероятность повреждения Q(т, Ат) за промежуток времени [т, т + Ат]. Для оценки вероятности появления повреждений используется метод доверительных интервалов.
Задачу интервального оценивания в самом общем виде можно сформулировать так: по данным выборки построить числовой интервал, относительно которого с заранее выбранной вероятностью можно сказать, что внутри этого интервала находится оцениваемый параметр. Интервальное оценивание необходимо при малом числе отказов, что характерно для корабельного электрооборудования, когда точечная оценка не может служить устойчивой характеристикой надёжности.
673
Власов А.Б., Буев С.А. Оценка показателей надёжности...
Доверительным интервалом [QH, QJ для исследуемого параметра Q называют такой интервал, относительно которого можно с заранее выбранной вероятностью 5=1 - у, близкой к единице, утверждать, что он содержит неизвестное значение параметра Q, т.е.
P{Qh < Q < Qe} > 1 - у,
где Qe и QH - верхняя и нижняя границы параметра Q; вероятность 5 - доверительная вероятность или коэффициент доверия, а у - уровень значимости.
Абсолютно достоверными границами для вероятности отказа Q(t, Дт) являются QH = 0 и Qe = 1. Объявляя любые другие границы, мы рискуем совершить ошибку. Вероятность такой ошибки характеризует уровень значимости у, а доверительная вероятность 5 гарантирует нам, что в 1005% случаев искомая вероятность отказа будет находиться в указанных границах.
Выбор величины доверительной вероятности зависит от вида изделия, его назначения и целей использования. Желание иметь большую доверительную вероятность (перестраховка) приводит к увеличению интервала, т.е. расширению границ искомого параметра. Для корабельного электрооборудования вполне достаточно величины 5 = 0,95 или 5 = 0,90 (Рябинин, Парфёнов, 1997).
Согласно общей теории интервальных оценок нижний доверительный предел п для неизвестной вероятностиp определяется как решение уравнения (Большее, Смирнов, 1983)
In = (р, n - р + 1) = 1 - P, (6)
где P - заданный коэффициент доверия (0,5 < P < 1); n - общее количество независимых испытаний; р -число испытаний, в которых наблюдался отказ; Ix(a, b) - функция В-распределения.
Функция В-распределения Ix(a, b) зависит от двух положительных параметров а и b и определяется на отрезке 0 < x < 1 формулой
Ix(a, b) = (НЕ (a, b)) |0х f'1 (1 - t)b-1 dt, (7)
где E(a,b) - так называемая В-функция Эйлера:
E(a,b) = Jo1 f'1 (1 - t)b-1 dt = (r(a)r(b)) / r(a + b). (8)
В принятых нами обозначениях формула (6) принимает вид
IQH(n, N - n + 1) = 1 - 5. (9)
Аналогично верхний доверительный предел представляет собой решение уравнения
IQe(n + 1, N - n) = 5, (10)
где n - число отказов на исследуемом интервале времени; N - общее число контактных соединений на обследуемом морском судне.
Пара чисел Qu и Qe, соответствующая одним и тем же значениям n, N и 5, определяет для неизвестной вероятности отказа Q(t, Д) доверительный интервал (QK, Qe) c коэффициентом доверия 25 - 1.
Правомочность применения формул (9) и (10) базируется на равенстве функций В-распределения и функций биноминального распределения, т.е.
Ie(n, N - n + 1) = X*=nN CNkQk(1 - Q)N-k, (11)
т.к. при ограниченном числе контактных соединений N корабельного электрооборудования и малом числе отказов п(т, Дт) оценка частоты отказов имеет биноминальное распределение.
Ещё одним способом определения интервальных границ Qe и Qu является приближённая оценка по формулам (Власов, 2006):
Qe ~ x[100(1 - P)%,2(n + 1)] / (2N- n + 1 + 0,5x[100(1 - P)%, 2(n + 1)]); (12)
QH ~ x[100P%,2n] / (2N- n + 1 + 0,5x[100P%, 2n]), (13)
где x = f(Q, 2a) - Q-процентная точка X-распределения с v = 2a степенями свободы, этот параметр определяется по табл. 2.2а из работы (Большев, Смирнов, 1983).
Рассмотрим третий способ определения доверительных границ, который основан на оценке биноминального распределения с помощью приближённой формулы Пуассона из уравнений:
Qe ~ x[100(1 - P)%,2(n + 1)] / 2N; (14)
QH ~ x[100P%,2n] / 2N, (15)
в которых числитель определяется, как и в формулах (12) и (13).
674
Вестник МГТУ, том 16, № 4, 2013 г.
стр.672-680
Для расчёта доверительного интервала без применения статистических таблиц используют формулу Муавра - Лапласа:
Qe ~ (2n + 1 + Л + V( z4sl + (4n - 4n2/N + 2 - 4n/N- 1/ N) z23l) / (2(N + z281)); (16)
Qh ~ (2n + 1 + z2a - V( z4n + (4n - 4n2/N + 2 - 4n/N - 1/ N) z2^i) / (2(N + z2^i)), (17)
где zS1 = Ф"1(^1) - аргумент функции нормального распределения, при котором функция нормального распределения Ф(д1) = д1.
В табл. 2 приведены аргументы функции для часто используемых значений коэффициентов доверия д1.
Таблица 2
81 0,900 0,950 0,975
z81 1,282 1,645 1,960
На практике нередко приходится сталкиваться с задачей определения доверительного интервала для вероятности повреждения, когда полученная из опыта частота повреждений равна нулю, либо тепловизионная диагностика не проводилась.
При количестве отказов п(т,Дт) = 0 и, соответственно, частоте отказов, равной нулю, нижняя граница будет нулевая. Формула для верхнего значения границы вероятности отказа может быть записана в виде
Qe = 1 - Ч1 = 1 - N(1 - 81). (18)
Например, в течение одного года не обнаружено ни одного дефекта на буксире-спасателе "Мурманрыба", имеющем 410 контактов. Задавая коэффициент доверия 8 = 0,95, определим верхнюю доверительную границу: Qe = 0,00728. Следовательно, с гарантией 95 % можно утверждать, что в течение данного года вероятность безотказной работы на траулере составляет не менее Rh = 0,99272.
От интервальной оценки вероятности отказа легко перейти к построению доверительного интервала I8 = [Rh(t, Дт); Rb(t, Дт)] для вероятности безотказной работы R (т, Дт) путём замены на
функцию надёжности и изменения знаков величин и неравенств:
Qh(t, Дт) < Q(t, Дт) < Qe(T, Дт); (19)
Qh(t, Дт) < 1 - RQ, Дт) < Qfr, Дт); (20)
1 - QeQ, Дт) < RQ, Дт) < 1 - Qh^, Дт). (21)
Таким образом, границы доверительного интервала для искомой вероятности безотказной работы будут следующие:
RhQ, Дт) = 1 - QeQ, Дт); (22)
ReQ, Дт) = 1 - Qh^, Дт). (23)
Например, на рыболовецком траулере "Альферас" в течение 2012 г. обнаружено n0 = 2 дефекта контактных соединений (N = 430). Требуется построить доверительные пределы для неизвестной вероятности повреждений Q (2012;1) контактных соединений в интервале времени, равном 1 году, с коэффициентом доверия 8 = 0,9. Точечную оценку неизвестной вероятности повреждения вычислим по формуле (2):
0‘(т,Дт) = Q(2012; 1) = п0(т,т + Дт) / N = 0,004651.
Коэффициент доверия одностороннего доверительного интервала 81 и уровень у1 значимости вычисляются по формулам:
81 = (1 + 8) / 2 = (1 + 0,9) / 2 = 0,95;
Y = 1 - 81 = 0,05.
По табл. 5.2 (Большее, Смирное, 1983) для 81 = 0,95 находим строку для n = 2 и два соседних столбца с величиной Со = 200, для которой Qe0 = 0,031 и Qh0 = 0,002, а для С = 500 имеем Qe1 = 0,012 и Qh1 = 0,001. Поскольку число С = N - n = 428, в нашем случае применим интерполяцию с коэффициентом
к = (С - С0) / (С - С1) = (428 - 200) / (500 - 200) = 0,76
и найдём Qe и Qh:
Qe = Qe0 Qe1 / (Qe0k + Qe1(1 - к)) = 0,31-0,012 / (0,031-0,76 + 0,012(1 - 0,76) = 0,014;
Qh = Qh0Qh1 / (QH0k + Qh1(1 - к)) = 0,002-0,001 / (0,002-0,76 + 0,001(1 - 0,76) = 0,001.
675
Власов А.Б., Буев С.А. Оценка показателей надёжности...
Таким образом, в силу правила с коэффициентом доверия 28 - 1 = 0,95 можно утверждать, что истинная вероятность повреждения контактного соединения в течение 2012 г. накрывается интервалом Jo,95= [0,001; 0,014] или Р{0,001 <0(2012,1) <0,014} >0,95.
Для оценки интервала воспользуемся формулами (12) и (13) по данным из табл. 2.2а (Большее, Смирнов, 1983):
Qe = x[5%, 2(2 + 1)] / (2-430 - 2 + 1 + 0,5x[5%, 2(2+1)]) = 12,592 / 865,296 = 0,014;
QH = x[95%, 2-2] / (2-430 - 2 + 1 + 0,5x[95%, 2-2]) ~ 0,711 / 859,356 ~ 0,001.
По приближённым формулам Пуассона (14) и (15) имеем:
Qe = x[5%, 2(2 + 1)] / 2-430 = 12,592 / 860 = 0,015;
QH = x[95%, 2-2] / 2-430 = 0,711/860 ~ 0,001.
Доверительный интервал по приближённым формулам Муавра - Лапласа (16) и (17) будет иметь следующие границы:
Qe = (2-2 + 1 + (1,645)2 + 1,645^((1,645)4 + 4 • 2 - 4 • 22 / 430 + 2 - 4 • 2 / 430 - 1 / 430) /
/ (2(430 + (1,645)2)) = (7,71 + 1,645^12,6479) / 865,42 = 13,56 / 865,42 = 0,016;
QH = 7,71 - 1,645^12,6479 / 865,42 = 1,8597 / 865,42 = 0,002.
В табл. 3 представлены сравнительные результаты.
Таблица 3
Результат Наименование способа определения доверительного интервала
По формулам (9), (10) По формулам (12), (13) По прибл. формулам Пуассона По прибл. формулам Муавра - Лапласа
Qe 0,014 0,014 0,015 0,016
Qn 0,001 0,001 0,001 0,002
Результаты показывают, что нижняя граница двустороннего интервала первых трёх методов в пределах данной точности одинакова.
Рассчитаем вероятности безотказной работы R (т) = 0,995349 при точечной вероятности повреждения Q (т) = 0,004651. Для крайних пределов вероятности повреждения получаем: Qe = 0,014; R,* = 0,986; Q* = 0,001; Rв* = 0,999. В итоге имеем: P{0,986 < R*(2012;1) < 0,999} > 0,95.
В результате выполненных расчётов и анализа данных можно сделать некоторые выводы. Во-первых, экспериментальные данные имеют некоторый разброс значений, обусловленный многочисленными факторами. Во-вторых, данные тепловизионной диагностики требуют тщательной обработки для сравнения данных по надёжности, полученных в разное время или с разных объектов. При исследовании надёжности элементов электрических систем кораблей вызывает затруднение совместная обработка всех накопленных статистических данных об их эксплуатации, т.к. определить идентичность условий и режимов работы распределительных устройств путём поверхностного анализа не всегда возможно. Ввиду этого, прежде чем приступить к совместной обработке информации, необходимо определить, случайным или неслучайным является расхождение между статистическими данными. Неслучайность расхождения статистики свидетельствует о наличии существенных различий в условиях или режимах эксплуатации морского оборудования, и, следовательно, последние обладают различными характеристиками надёжности.
3. Проверка достоверности полученных данных на основе статистических гипотез о равенстве двух вероятностей отказа
При анализе данных количественной инфракрасной термографии должна производиться их статистическая обработка, направленная на проверку тех или иных статистических гипотез о виде или параметрах неизвестного закона распределения (Власов, 2006).
Рассматриваемый метод используется для сравнения надёжности однотипного оборудования на разных кораблях или одного и того же оборудования в разные промежутки времени. В первом случае нас интересует, насколько существенны внешние факторы, а во втором - например, можно ли считать, что наступило старение.
Имеем две выборки с разными числами элементов N и N2, в которых отказали соответственно п\ и п2 элементов. Статистические частоты отказов для данных выборок будут:
676
Вестник МГТУ, том 16, № 4, 2013 г.
стр.672-680
Q1 = n / М; (24)
Q2* = «2 / N2. (25)
Они соответствуют неизвестным нам истинным значениям вероятностей отказа Q1 и Q2. Для однозначности предположим, что Q1 > 02.
В общем виде задача проверки статистической гипотезы формулируется следующим образом (Рябинин, Парфёнов, 1997). Имеем статистические данные в виде одной или нескольких выборок. Высказывается предположение о параметрах или о законе распределения выборок. Это основное предположение называют нулевой гипотезой Н0. Альтернативой гипотезе Н0 может быть одна или две гипотезы, называемые альтернативными или конкурирующими гипотезами Н1 и H2. Таким образом, задача проверки гипотезы Н0 для выборки, состоящей из n наблюдений x1 ,x2,...,xn случайной величиныX, заключается в выработке определённых правил, по которым мы решаем, принять гипотезу Н0 или гипотезу Н1. Всё возможное множество выборок объёма n можно разделить на два непересекающихся подмножества S п1 и S n2, т.е. S п1 П S п2 = 0. Если наша выборка (или её характеристики, параметры) попадает в область Sn1, то принимается исходная проверяемая гипотеза Н0, если же в область Sn2, то нулевая гипотеза Н0 отвергается, и принимается альтернативная гипотеза Н1.
Подмножество S n1 называют областью допустимых значений, а Sn2 - критической областью. Выбор одной из областей однозначно определяет и другую область, т.е. задача сводится к определению одной области, а вторая получается автоматически единственным образом. Возникает необходимость формирования критериев, принципов или правил для построения критической области S n2 или области допустимых значений S n1. При выборе критической области и принятии или отклонении гипотезы Н0 по случайной выборке принятое решение соответствует истине с некоторой доверительной вероятностью. Схема возможных исходов представлена в табл. 4.
Таблица 4
Г ипотеза Н0 Принятое решение Критерий Пояснение
Неверна Принято в Ошибка 2-го рода
Не принято 1 - в Мощность критерия
Верна Принято а Доверительная вероятность
Не принято 1 - а Уровень значимости
Рассмотрим нулевую гипотезу Н0 = {Q1 = Q2}, что различие между статистическими частотами Q1 и Q2 случайно и несущественно. Если гипотеза Н0 будет отвергнута, то принимается одна из двух конкурирующих гипотез: Н1 = {Q1 < Q2} или Н2 = {Q1 > Q2}.
Для наглядности и удобства расчётов исходные данные представим в виде сопряжённых признаков 2 х 2, показанных в табл. 5.
Таблица 5
№ выборки Всего Число неисправных Число исправных Q*
элементов элементов элементов
1 N1 n1 r1 = N1 - n1 n1/N1
2 n2 n2 r2 = N2 - n2 nr2iN'2
Всего N = N1 + N2 n = n1 + n2 R = r1 + r2 n/N
Имеем независимые случайные величины G1 и G2, которые могут принимать значения G1 = 0, 1, 2, ., n1, ., N1 и G2 = 0, 1, 2, ., n2, ., N2. Вероятность появления ровно n1 событий при N1 испытаниях в теории вероятности определяется по биноминальному закону
P{n1} = CN1n1 • pn1 • qN1-n1, (26)
где p - вероятность появления одного события; q = 1 - p; Cm"1 - число сочетаний из N1 по n1.
При допущении справедливости Н0 = {Q1 = Q2 = Q} и в силу независимости G1 и G2 вероятность совместного появления n1 отказов в первой выборке и n2 отказов во второй можем записать в виде произведения
P{n1n2} = P^P^} =
__ s~l n1
= CN
n1
CN1
• СN2
• Qn1 • (1 - Q)N1-n1
n2 • Qnl+n2 • (1 - Q)
• CN2n2 • Q"2 • (1 - Q)
N1+N2-n1-n2
N2-n2
(27)
где Q - неизвестная вероятность появления отказа.
Для исключения этой неизвестной вероятности Q найдём вероятность появления n отказов из суммарной выборки объёмом N1 + N2 = N
677
Власов А.Б., Буев С.А. Оценка показателей надёжности...
Т)<-.Л _|_ _ ,„■» _ г'' nl+n2 г\п\+Ш /л /'j\Nl+N2-nl-n2 /лп\
P{ni + П2 — n} — CM+N2 • Q • (1 - Q) (28)
и, разделив (27) на (28), получим величину
P{n1n2} / P{n1 + n2 — n} — СМП1 • Cni2 / CN1+N^^n2, (29)
которая является условной вероятностью появления n1 отказов в первой выборке, при условии, что в обеих выборках было ровно n отказов. Действительно, полную вероятность (27) можно представить в виде произведения вероятностей зависимых событий
P{nin2} — P{ni} P{ni|n,Ho}, (30)
где P{n1 | n,H0} - условная вероятность появления n1 отказов при условии, что в обеих выборках справедлива гипотеза H0. Но эта условная вероятность уже не зависит от неизвестного параметра Q и может быть найдена через факториалы соответствующих величин:
P{n1 | n,Ho} — P{n1n2} / P{n} — CMn1 • CN2n2 / CNn — n1 !N2!n!(N - n)! / n1!(N1 - n1)!n2!(N2 - n2)!N! —
N1!N2!n!r! / n1!r1!n2!r2!N!. (31)
Она имеет гипергеометрическое распределение с интегральной функцией распределения F{n1 | n,H0}, по которой можем определить вероятность попадания случайной величины G1 в критическую область Sn2 в виде суммы условных вероятностей
P1{G1 < щ | n} = Xk—on1P{k | n} — P{0| n} + P{1| n} + ... + P{n1 | n} (32)
или вероятность попадания G2 в критическую область Sn2
P2{G1 > n1 | n} = Tk—n1nP{k | n} — P{n1 | n} + P{n1+1| n} + . + P{n |n}. (33)
Для проверки гипотезы Н0 — {Q1 — Q2} при двусторонней альтернативе Н1 — {Q1 < Q2} или Н2 — {Q1 > Q2} вероятности (32) и (33) сравниваются с величиной а/2. Если обе вероятности P1{G1 < n1 | n} и P2{G1 > n1 | n} окажутся больше уровня значимости а/2, то гипотеза Н0 принимается. При односторонней альтернативе сравнение производится c величиной а. Если одна из этих вероятностей будет меньше а/2, то гипотеза Н0 отклоняется в пользу Н1 и Н2 соответственно.
Вычисление условной вероятности по формуле (31) путём логарифмирования представлено
ниже:
lgP{ nI\n} — (lgN1! + lgN2! + lgn! + lgr! - lgN!) -
- (lgn1! + lgr! + lgn2!) = X1 - X2, (34)
где X1 - постоянная величина при вычислении вероятностей по формулам (32) и (33), а Х2 - переменная величина, вычисляемая для удобства в табличной форме.
Алгоритм действий при проверке гипотезы о равенстве двух вероятностей отказа описан в (Рябинин, Парфёнов, 1997) и заключается в следующем:
1) составляется таблица сопряжённых признаков;
2) вычисляются десятичные логарифмы по формуле (34) для значений G1 = (0, 1, 2, ..., n1, ..., n);
3) определяются вероятности попадания в критическую область по формулам (32) и (33);
4) полученные вероятности сравниваются с уровнем значимости а/2 при двусторонней альтернативе или с величиной а при односторонней альтернативе: при P1{G1 < n1 | n} > а/2 и P2{G1 > n1 | n} > а/2 (а) принимается гипотеза Н0, если P1{G1 < n1 | n} < а/2 (а), то принимается гипотеза Н1 — {Q1 < Q2}; если же P2{G1 > n1 | n} < а/2 (а), то принимается вторая конкурирующая гипотеза Н2 — {Q1 > Q2}.
Перейдём к выполнению практической задачи. В результате тепловизионной диагностики в 2012 г. было зарегистрировано 3 отказа из 402 контактных соединений на объекте I (транспортнобуксирное судно "Десна") и 7 отказов из 400 контактных соединений на объекте II (самоходно-плавучий кран "СПК-395"). Нам необходимо сравнить надёжность электрооборудования на разных типах судов. Можем ли мы считать различие несущественным, а следовательно, статистические данные, полученные с разных морских судов, - однородными?
Представим исходные данные в табл. 6.
Высказываем гипотезу Н0 — {Q1 — Q2} о равенстве вероятностей отказа на объектах I и II при односторонней конкурирующей гипотезе Н1 — {Q1 < Q2}. По формуле (34) вычисляем постоянную сумму (для вычисления воспользуемся табл. 7.6 (Большее, Смирнов, 1983) десятичных логарифмов
факториалов):
678
Вестник МГТУ, том 16, № 4, 2013 г.
стр.672-680
X = lgM! + lgN2! + lgn! + lgr! - lgN = lg402! + lg400! + lgl0! + lg792! - lg802! =
= 1720,3626831.
Таблица 6
№ выборки Всего Число неисправных Число исправных Q*
элементов элементов элементов
1 402 3 399 0,007463
2 400 7 393 0,017500
Всего 802 10 792 0,012468
Переменную сумму X2 вычисляем в табл. 7.
Таблица 7
П1 Г1 0 402 2 401 2 400 3 399
П2 Г2 10 390 9 391 8 392 7 393
1 2 3 4 5
lgn1! 0 0,3010300 0,3010300 0,7781513
lgn! 874,0137846 871,4095586 868,8064142 866,2043542
lgn2! 6,5597630 5,5597630 4,6055205 3,7024305
lgr2! 842,8355064 845,4272406 848,0205267 850,6149192
X2 1723,409054 1722,6975922 1721,7334914 1721,2998552
Вычисляем логарифмы условных вероятностей по (34) и соответствующие вероятности:
lgP{0| 10} = 1720,3626831 - 1723,409054 = - 3,0463708;
P{0| 10} = 0,0008987;
lgP{1| 10} = 1720,3626831 - 1722,6975922 = - 2,3349091;
P{1| 10} = 0,0046247;
lgP{2| 10} = 1720,3626831 - 1721,7334914 = - 1,3708083;
P{2| 10} = 0,0425786;
lgP{3| 10} = 1720,3626831 - 1722,6975922 = - 0,9371720;
P{3| 10} = 0,1155654.
По формуле (32) находим:
< 3 | 10} = Хк=03Р{£ | 10} = P{0| 10} + P{1| 10} + P{2| 10} + P{3| 10} =
= 0,1636675.
При уровне значимости а = 10 % полученная вероятность превосходит а, поэтому принимаем гипотезу Н0 = {Q1 = Q2}. Наблюдаемая разница в числе отказов носит случайный характер.
Проверим равенство статистических вероятностей отказа Q1 и Q2 методом доверительных интервалов по приближённым формулам (12) и (13). Верхняя доверительная граница Q1e при п1 = 3, N1 = 402 для уровня значимости а = 0,10 будет следующей:
Q1e* = х[10%, 2(3 + 1)] / (2 • 402 - 3 + 1 + 0,5х[10%, 2(3 + 1)]) = 13,362 / 867,81 = 0,015397.
Нижняя доверительная граница Q2h при п2 = 7, N1 = 400 для уровня значимости а = 0,10 будет следующей:
Q1e* = х[90%, 2 • 7] / (2 • 400 - 7 + 1 + 0,5х[90%, 2 • 7]) = 7,790 / 797,895 = 0,009763.
Таким образом, имеем область пересечения доверительных интервалов на участке (0,009763; 0,015397), что позволяет принять гипотезу Н0 = {Q1 = Q2}.
4. Вероятность возникновения страховых рисков
При рассматриваемых отказах электрооборудования морского судна появляется возможность возникновения страхового случая по всем рискам страхования (пожар, столкновение с другим судном, штормовое повреждение и т.д.).
По статистике в 18 % случаев причиной возникновения пожара на корабле являлся выход из строя электрооборудования (Подобед, Панкратов, 2011).
679
Власов А.Б., Буев С.А. Оценка показателей надёжности...
По данным Центрального союза морских страховщиков (CEFOR), средняя величина страховой выплаты по риску пожара за первую половину 2012 г. составила 4,5 млн долл. США (Отчёт об убытках..., 2012).
При повышении вероятности безотказной работы электрооборудования риск возникновения уменьшается.
По данным статистики на объекте I (транспортно-буксирное судное "Десна") при применении метода тепловизионной диагностики вероятность безотказной работы R (т) повысилась с 0,982587 в 2009 г. до 0,992537 в 2012 г. При этом риск возникновения пожара Р(пож) в среднем уменьшился на 9 %.
Р2012(пож) = (1 - R*2012) • Лр (пож) / (1 - R*2009) = 9 %.
5. Заключение
Экономический эффект от внедрения метода тепловизионной диагностики для компаний, эксплуатирующих морские суда, обеспечивается за счёт уменьшения эксплуатационных затрат, убытков от простоя судна в связи с ремонтом, уменьшения страховой премии. Обзорные сведения, связанные с экономическим эффектом от применения тепловизионной техники в различных отраслях народного хозяйства, представлены в работе (Власов, 2006).
Применение систем технической диагностики позволяет устанавливать оптимальные сроки выполнения планово-предупредительного ТО и ремонта, предотвращать или существенно сокращать число отказов судового оборудования и вынужденных остановок в море, сокращать время поиска неисправностей. Расчёты показывают: при применении разработанного метода R (т) = 0,995347, при точечной вероятности повреждения Q (т) = 0,004651.
Экономический эффект при избегании ремонтных работ складывается из стоимости простоя судна и стоимости выполнения ремонта и в среднем для рыболовецкого судна доходит до 3 млн руб. в год (Базарное, 2008).
По отношению к морскому страхованию предотвращение возможных страховых случаев, вызванных выходом из строя электрооборудования и возникновением пожаров, является важной задачей, которую способна решать тепловизионная диагностика.
Выявление дефектов электрооборудования на ранней стадии развития является одним из элементов противопожарной профилактики.
Наличие у судовладельца отчёта о проведённой тепловизионной диагностике электрооборудования является документом, который принимается Регистром Судоходства и страховыми организациями как подтверждение технической исправности судна.
Литература
Базарнов В. Безопасность на море и на суше. Морской бизнес, № 11, с. 12-15, 2008.
Большев Л.Н., Смирнов Н.В. Таблицы математической статистики. М., Наука, Главная редакция физико-математической литературы, 416 с., 1983.
Власов А.Б. Методы статистической обработки тепловизионных данных по параметрам надёжности электрооборудования. Методы и средства оценки состояния энергетического оборудования.
Петербургский энергетический ин-т повышения квалификации рук. работников и специалистов, вып. 17, c. 16-31, 2002.
Власов А.Б. Модели и методы термографической диагностики объектов энергетики. М., Колос, 280 с., 2006.
Власов А.Б., Джура А.В. Система сопровождения оборудования по техническому состоянию (тепловизионная диагностика). Методы и средства оценки состояния энергетического оборудования. Петербургский энергетический ин-т повышения квалификации рук. работников и специалистов, вып. 20, c. 294-301, 2002.
Отчёт об убытках Центрального союза морских страховщиков за первую половину 2012 года. 2012. URL: http://www.cefor.no.
Подобед В.А., Панкратов А.А. Анализ пожаров на судах рыбопромыслового флота Северного бассейна. Вестник МГТУ, т. 14, № 4, c. 728-732, 2011.
Рябинин И.А., Парфёнов Ю.М. Надёжность, живучесть и безопасность корабельных электроэнергетических систем. СПб., Изд-во Военно-морской академии им. адмирала флота Советского Союза Н.Г. Кузнецова, 431 с., 1997.
680
