Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Клемин А.И. Инженерные вероятностные расчеты при проектировании ядерных реакторов

.pdf
Скачиваний:
34
Добавлен:
25.10.2023
Размер:
14.61 Mб
Скачать

цим вынести из-под знака интеграла, взяв их в точке максимума функции ехр [—(а/0)2 ], т. е. при а = 0. В итоге находим

| ¥ -

[/• (а)] а ехр

 

 

 

da :

 

 

\V2)

dr

r = R/)R2

—д

 

 

 

 

 

 

 

 

 

так как

г (0)

R/V2.

 

 

 

 

 

 

 

 

 

 

 

Аналогичным образом вычисляем интеграл в знаменателе фор­

мулы (5.74):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

' Ч'2

[г (а)] ехр

 

 

 

а_

\ 2

 

 

 

 

 

 

 

 

 

а

da:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y2 [/-(0)] j

ехр

 

 

 

da « У2noW

[

^

 

 

 

—д

 

 

 

 

 

 

 

 

 

 

 

Воспользовавшись полученными выражениями, формулу (5.74)

можно преобразовать к виду

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d4

 

 

 

 

Ар = 2 / г

\ І

і

dr

r = R/Y2

 

R

dr

r=R/Y2

Vі 4 rr

(^i 77\

 

]г = к/ r

 

 

 

'\

"* \г=к/у 'і

В частном

случае,

когда W (г) =

J0

 

2,405

г),

где

#э —-эффек-

тивныи

радиус

активной

зоны,

получаем

 

 

 

 

 

 

 

 

 

j

(

2 ' 4

0 5

 

 

 

 

 

 

 

Ар

R

 

2,405

 

у 2

-

 

 

 

 

 

(5.78)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

УТл

 

Rs

(

2,405

 

R

 

 

 

 

 

 

 

 

 

 

 

У2

'

Rs

 

 

 

 

Здесь J0

и /х —функции

Бесселя

 

1-го рода соответственно нулевого

и первого порядков.

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом,

формулы (5.77)

и

(5.78)

устанавливают

связь

между характеристиками

рассеяния

конструкционных

параметров

(в пределах допусков) и максимальным отклонением запаса реактив­ ности, возможным при специальном размещении технологических каналов в активной зоне. Средние квадратические отклонения а г можно вычислить по результатам замеров конструкционных пара­ метров xt на определенной партии каналов, см. гл. П . Если партия еще не изготовлена или недостаточно исследована, то можно оценить величины ot приближенно, положив at = Aj/З или же (это будет в запас) приняв ot = Аг / |/2я , как для треугольного распределения Симпсона [6].

Г л а в а

6.

О Ц Е Н К А Н А Д Е Ж Н О С Т И Р Е А К Т О Р Н О Й У С Т А Н О В К И *

§ 6.1.

Основные количественные показатели

общей

надежности реакторной установки

Говорить о надежности, повышать ее, поддерживать в эксплуата­ ции можно, и не интересуясь ее количественной стороной. Такой ка­ чественный подход к надежности существует с момента появления инженерной деятельности, и в наши дни он вовсе не исключается, а иногда остается единственно возможным. Однако более рациональ­ ным, эффективным, а потому и более прогрессивным, является под­ ход, базирующийся на количественной основе. Он позволяет ква­ лифицированно сравнивать изделия по надежности друг с другом, оп­ тимизировать их по критерию надежности, выбирать действительно лучший вариант, оценивать эффективность того или иного способа повышения надежности, выяснять в количественном выражении вли­ яние на надежность самых различных факторов, схемных решений, резервирования, материалов, технологии, отдельных деталей, узлов ит. д.

Основными понятиями современного количественного подхода к проблеме надежности ядерных реакторных установок являются следующие.

Н а д е ж н о с т ь р е а к т о р н о й у с т а н о в к и (в широ­ ком смысле) есть свойство, обусловленное ее безотказностью, ремон­ топригодностью и долговечностью и обеспечивающее нормальное выполнение установкой требуемой задачи в заданном объеме и в за­ данных условиях эксплуатации.

О т к а з р е а к т о р н о й у с т а н о в к и — событие, в резуль­ тате которого установка полностью или частично утрачивает рабо­ тоспособность. П о л н ы м называется отказ, приводящий к вынуж­ денной срочной остановке реактора, а следовательно, и всей установ­ ки в целом. Под ч а с т и ч н ы м понимают отказ, приводящий к не­ обходимости частичного снижения рабочих параметров реакторной установки. Без отказов проблемы надежности вообще и надежности ядерных реакторных установок в частности не существовало бы. Про­ блемы надежности не существовало бы даже в том случае, когда отка­ зы и возникали бы, но заранее был бы известен момент наступления каждого из них, другими словами, если бы отказы носили не случай­ ный, а строго детерминированный характер.

В отличие от отказа под н е и с п р а в н о с т ь ю

будем пони­

мать такой дефект оборудования, который не требует

немедленного

прекращения его эксплуатации, т. е. с которым оборудование в тече-

*

П о д

р е а к т о р н о й у с т а н о в к о й

д л я о п р е д е л е н н о с т и т а м , где это не

огово­

р е н о о с о б о ,

б у д е м п о н и м а т ь л и б о

р е а к т о р н у ю п а р о п р о и з в о д и т е л ь н у ю

уста­

н о в к у , л и б о о т д е л ь н ы й б л о к А Э С , л и б о м о н о б л о ч н у ю А Э С в ц е л о м .

 

4 Зак.

1282

 

81

ниє определенного времени еще может выполнять свои основные функции.'

Б е з о т к а з н о с т ь р е а к т о р н о й у с т а н о в к и (или надежность в узком смысле) количественно характеризуется вероят­ ностью Р (і) того, что в заданных режиме и условиях эксплуатации установка проработает без отказов в течение времени і:

 

Р ( 0 = Р { т > 0 .

 

(6-і)

гдет — время безотказной работы

реакторной

установки, или нара­

ботка на отказ, — случайная величина.

 

 

Р е м о н т о п р и г о д н о с т ь р е а к т о р н о й у с т а н о в -

к и

количественно характеризуется суммарным

временем аварий­

ных

ремонтов 9р и суммарным временем планово-предупредитель­

ных

ремонтов Ґ1ПР (профилактических работ) 6П

за весь период ее

эксплуатации

 

 

 

 

8р = 2 6р ;,

6П = 2 0D ;.

 

/

/

 

 

Д о л г о в е ч н о с т ь р е а к т о р н о й

у с т а н о в к и ко­

личественно характеризуется ее сроком службы т0 . с

и полным тех­

ническим ресурсом

 

т т . р = X Xj,

(6.2)

т. е. полным временем безотказной работы реакторной установки до момента вывода ее из эксплуатации (т; — интервалы работоспособ­ ности). Срок службы

т с . о = Т т . р + Є = 2 (т, + BPJ- + 0 n J + Q'j)

(6.3)

і

 

есть полная календарная продолжительность эксплуатации реактор­

ной установки (6 — полное

время простоя установки,

включающее

в себя периоды аварийного

ремонта — 0р ;, ППР —

0п ; - и прочие

остановки — Q'j, например, для перегрузки ядерного топлива и т. д.). Срок службы стационарных энергетических реакторных установок обычно составляет 20—30 лет [41—43].

Количественные показатели общей надежности реакторной уста­ новки за период времени t следующие.

1. Коэффициент технического использования (КТИ), называемый иногда реакторостроителями коэффициентом использования кален­ дарного времени (КИКВ), равен отношению «чистого» времени рабо­

ты т (І) установки за некоторый

календарный период эксплуатации

к величине этого периода t:

 

 

/Ст. и О = т {t)lt

= х (0ЛЄ (0 + т (01,

(6.4)

где G (t) — полное время простоя установки за период t. Для установки, снятой с эксплуатации, Кт.,, = тт > р с , с . Для доста­ точно большого периода t КТИ численно равен вероятности застать

реакторную установку работающей в произвольно выбранный мо­ мент времени (внутри этого периода).

КТИ характеризует в основном надежность установки в отноше­ нии полных отказов, приводящих к ее остановке, и плановых ремон­ тов. Чем больше таких отказов, чем больше времени тратится на их устранение и на проведение плановых ремонтов, тем ниже КТИ. Обычно для АЭС Кг. „ ^ 75% [44].

2. Коэффициент использования установленной мощности(КИУМ) равен отношению фактически выработанной мощности реакторной

установкой за период эксплуатации t

к мощности,

которую бы она

выработала за время t, работая без

остановок

на

номинальном

уровне

мощности:

 

 

 

 

 

 

 

 

 

Кп.у.

м if) —

I Q (0 dt

//Q4

=

т (О Q> - Q H

=

 

 

 

 

= Я т . и (0 QVQh>

 

 

(6-5)

где Q(t) — мощность

установки

в произвольный

момент времени;

Qt — средний

уровень

мощности

за период времени

t

(исключая

время

остановок); QH

— номинальная

(установленная)

мощность.

Если реакторная установка эксплуатируется только на

номиналь­

ном уровне мощности QH, то Ки. у . M.(t)

= Кт. и (О-

 

 

 

Таким образом,

КИУМ характеризует

надежность

реакторной

установки не только в отношении полных,

но и частичных отказов,

которые не приводят к ее остановке, а требуют снижения мощности.

Чем

ниже

мощность работающей установки по сравнению с кана­

лом,

тем

ниже КИУМ при постоянном КТИ. Обычно для АЭС

Л и у . м ^

70% [41, 44]. Например, за 9 месяцев 1969 г. КИУМ для

БАЭС составил 71% [44].

3.Коэффициент готовности равен отношению чистого времени работы x(t) реакторной установки за календарный период эксплуа­ тации t к сумме этого времени и продолжительности аварийных ре­ монтов Qp(t) за период t:

К? (0 =

т (t)/[x (t) + 9Р

(01.

(6.6)

Коэффициент готовности

является характеристикой

надежности

установки за период, когда

не проводятся

плановые ремонты и пе­

регрузки горючего. Поэтому он удобен для оценки надежности транс­ портных реакторных установок, ибо численно равен вероятности застать реакторную установку работающей в произвольный момент

времени 0<t'<.t

в период ее нахождения, например, в плавании

(когда нет плановых

остановок для ремонтов

и перегрузки

или

они очень редки). Для АЭС, как правило, Кг ^

90% [45, 46].

4. Вероятность исправной работы реакторной установки в ин­

тервале времени от 4 до t2 (в общем случае, когда состояние

уста­

новки в момент tx

неизвестно):

 

 

Яі ih

- *i) = К,. „ (0 Р (t2 -

h),

(6.7)

4*

 

 

 

83

где Р(/2 — /і) — вероятность безотказной работы установки в те­ чение времени (t2 — tx), причем 0 < tx < t2 < t.

Иногда, интересуясь лишь одной стороной надежности — безот­ казностью, т. е. рассматривая надежность в узком смысле, в каче­ стве ее количественного показателя принимают вероятность безот­ казной работы (6.1) в течение заданного периода времени t:

R (0 = Р (0.

(6.8)

Показатели надежности реакторной установки (6.4)—(6.8) мо­ гут быть рассчитаны, если известны числовые характеристики, во-первых, безотказности: параметр потока отказов или интенсив­ ность отказов X, средняя наработка на отказ Т0 ; во-вторых, ре­ монтопригодности: интенсивность восстановления работоспособ­ ности установки- после отказа А.п, средняя продолжительность аварийного ремонта Г в . Остановимся подробнее на этих характе­ ристиках надежности.

§ 6.2. Характеристики безотказности реакторной

установки и ее оборудования

Параметр потока отказов. В период эксплуатации реакторной установки отказы возникают всегда в случайные моменты времени tx, t2, ... Эти моменты образуют случайный поток отказов. Можно рассматривать поток полных отказов реакторной установки (ава­ рийных остановок), поток частичных отказов, поток неисправно­ стей, поток отказов каналов активной зоны и т. д. Без ущерба для общности рассмотрения в этом параграфе будем предполагать, что длительность остановок реакторной установки много меньше чи­ стого времени ее работы. Разобьем срок службы или рассматривае­

мый интервал

эксплуатации установки на /V равных промежутков

А? = rc_c/N.

Будем регистрировать количество каких-то опреде­

ленных отказов установки, попавших в каждый из этих интервалов

тх, т2,

ты. Отложим на каждом из интервалов ординаты, рав­

ные mh/At,

и проведем через них горизонтальные прямые в пределах

интервалов. Получим ступенчатую кривую (рис. 19). При

достаточ­

но большом N

эта. эмпирическая

кривая будет близка

к функции

X(t), называемой п а р а м е т р о м

п о т о к а о т к а з о в р е ­

а к т о р н о й

у с т а н о в к и ,

 

 

 

X (t) = m (t,

t + At)/At,

(6.9)

где m \t, t + At) — число отказов в достаточно малом интервале времени At, построенном около момента t. Параметр потока отказов представляет собой среднее число отказов за достаточно малую единицу времени около момента t. Зная X(t), можно найти среднее количество отказов в любой отрезок времени в течение эксплуа­ тации:

. rn (tx, t2) = \X(t) dt = X(t2 — tx).

(6.10)

' l

 

Функцию X(t) можно рассматривать не только для реакторной установки в целом, но и для любого элемента ее оборудования. Од­ нако надо иметь в виду одно принципиальное обстоятельство. О по­ токе отказов и о его параметре можно говорить только применитель­ но к тем изделиям, которые ремонтируются (восстанавлііЕаются) в процессе эксплуатации. Если же рассматривать перемонтируемые (иевосстанавливаемые) изделия, например большинство типов

каналов

или твэлов реактора, то для них понятия параметра потока

отказов

Ці) не существует

(поскольку самого потока

нет). Вместо

него используется понятие

и н т е н с и в н о с т и

о т к а з о в

Л

 

 

 

 

 

—1

1

I

1

л

I

ш!

0

 

 

 

tec. *

Р и с . 19. П а р а м е т р п о т о к а

о т к а з о в

д л я р е м о н т и р у е м ы х и з д е ­

л и й и л и и н т е н с и в н о с т ь о т к а з о в д л я н е р е м о н т и р у е м ы х .

Щ). Подробнее

об этом см.

ниже.

Часто смешивают два упомя­

нутых (принципиально различных) понятия, употребляя во всех случаях последнее. Путаницу и ошибки в этом вопросе исключает специально введенный ГОСТ 13377—67 [47].

Зависимость % от времени. Обработка статистических данных по отказам реакторных установок и их элементов, а также теорети­ ческие соображения свидетельствуют о том, что кривая %(t) для них имеет вид, изображенный на рис. 19 [15, 48,- 49]. Подобная же зависимость параметра потока отказов и интенсивности отказов

от

времени наблюдается для многих изделий смежных отраслей

машиностроения

[7, 8,

18,

30, 50—55]. Характер этой зависи­

мости таков, что

позволяет

выделить три периода работы изделия

на

которых

X(t)

ведет

себя

принципиально

различным образом

 

I . Период

приработки

начальный период

работы устройства;

когда отмечается повышенная частота отказов (или вероятностьт отдельного отказа), связанная с выявлением скрытых дефектов проектирования, изготовления (включая дефекты материалов), монтажа и наладки. Для больших энергетических реакторных установок он может составлять несколько лет.

По окончании периода приработки отказы установки опреде­ ляются, в основном, только эксплуатационными факторами, так как дефектные элементы конструкций и систем к этому времени уже вышли из строя (условно говоря, отбраковалнсь) и остались только работоспособные, качественно выполненные. Поэтому период при­ работки изделия может служить объективной мерой для установ­ ления гарантийного срока, в течение которого завод-изготовитель должен нести финансовые издержки за отказы изделия. Разумеется, этот гарантийный срок не имеет ничего общего со сроком службы изделия (как, например, гарантия на автомобиль 1 год не определяет его долговечности).

I I .

Период нормальной

эксплуатации, характеризующийся

прак­

тически постоянной Х((), связанной с окончанием периода

прира­

ботки

и со стабилизацией

эксплуатационных факторов,

влияющих

на надежность

изделия.

 

 

 

I I I .

Период

старения

(износа)—заключительный

период ра­

боты, когда начинают сказываться необратимые физико-химические изменения в материалах и конструкциях, связанные с износом (ста­ рением), приводящие к возрастанию частоты отказов (или вероят­ ности отдельного отказа). Это период постепенных (т. е. происхо­ дящих в результате постепенного изменения определяющих пара­ метров) отказов. Обычно ядерные реакторные установки выводят из эксплуатации, когда чрезмерно возрастает количество износных отказов основной части оборудования. Предотвращают износные отказы путем выбора оптимальных сроков для капитального ремон­

та установки, ППР, путем регулярных

осмотров и контроля состоя­

ния оборудования,

замены элементов

установки до выхода

их из

строя и т. д.

 

 

 

 

Итак, весь срок

службы реакторной установки можно разбить

на интервалы, где

параметр потока

отказов (число отказов

в еди­

ницу времени)

 

 

 

 

kt

(t) « const,

<t<zti.

(6.11)

В качестве таких интервалов можно рассматривать перечисленные три с параметрами:

 

 

^•прираб

^-нор < С ^-стар-

 

 

(6.. 12)

Пуассоновский

поток

отказов.

В

пределах

интервала,

где

X(t) = X =

const,

поток

отказов реакторной

установки

с

боль­

шой достоверностью можно считать простейшим

стационарным

потоком,

поскольку он

удовлетворяет

всем

свойствам

послед­

него: 1) стационарности

(X =

const);

2) отсутствию последействия

(отказы установки

обычно устраняются

полностью, наиболее

часто

выходящие из строя элементы заменяются новыми, т. е. можно счи­ тать, что надежность реакторной установки возвращается практи­ чески к исходному уровню после каждого восстановления и после­ дующие отказы практически не зависят от предыдущих); 3) орди-

нарности (отказы установки — события редкие, и возникновение одновременно двух или более независимых отказов маловероятно, зависимые же отказы, происходящие один за другим, следует счи­ тать одним отказом реакторной установки).

Кроме отмеченных есть еще одно обстоятельство, которое под­ тверждает, что поток отказов реакторной установки близок, к про­ стейшему. Оно заключается в том, что этот поток складывается из большого количества обычно независимых потоков отказов от­ дельных составных частей, устройств и элементов установки. Каж­ дый из этих потоков может быть и не простейшим, однако их сумма (при условии, что параметр А,,, для отдельного потока отказов мал по сравнению с параметром общего потока отказов реактора) в со­ ответствии с предельной теоремой Пальма будет простейшим (пуассоновским) потоком [18]. Следовательно, вероятность возникнове­ ния т отказов реакторной установки за время /, на протяжении которого % = const, будет описываться законом Пуассона (3.71).

В общем случае, когда % Ф const на рассматриваемом интервале времени t, но остальные перечисленные выше условия для потока выполняются; можно показать, что поток отказов будет простей­ шим нестационарным потоком или пуассоновским потоком с пере­

менным параметром [7, 18]. Для такого

потока

 

 

t

 

t

 

P (m) =

jо X (t) dt

e x P

nI' %{t)dt

(6.13)

В этом случае для расчета P{tn) не обязательно иметь точную за­ висимость X(t). Достаточно знать среднее значение X на рассматри­ ваемом интервале і:

X=j

. ^%(t).dt

(6.14)

 

о

 

 

или среднее число отказов \t

за

время t. Из выражений

(6.13),

(6.14) получаем

 

 

 

Р(т)=Ш1

ml

ехр( — Xt).

(6.15)

Поскольку вид формул (3.71) и (6.15) совершенно идентичен, даль­ нейшее рассмотрение будем вести для случая К = const.

Функция (закон) надежности. Вероятность безотказной работы реакторной установки в течение любого периода времени t (6.1) будет равна вероятности Р(т) при т = 0:

Р (0 = ехр. (—U).

(6.16)

Эту зависимость называют функцией или законом надежности. Фор­ мула (6.16) предполагает, что в момент отсчета времени t — 0 Р(0) = 1, т. е. установка находится в работоспособном состоянии.

Если в качестве момента / = 0 рассматривать момент пуска уста­ новки, для которого в общем случае может быть

Р (0) = / ? „ < ! ,

(6.17)

то вероятность безотказной работы реакторной установки в тече­ ние времени і (включающего момент пуска t= 0) следует записать в виде

R (/) - R0P (0,

(6.18)

где R0 — вероятность, что отказ в момент пуска не наступит. Из формул (6.13) и (6.15) вытекают следствия, полезные для практики вычисления функции надежности Р(/). Рассмотрим два изделия:

Р и с . 20. В е р о я т н о с т ь

б е з о т к а з н о м

р а б о т ы

Р (/) (а)

д л я

и з д е л и й

с в о з р а с т а ю щ е й

(1) и у б ы в а ю щ е й (2)

ф у н к ц и е й

А (/)

(б).

одно с возрастающей функцией

в интервале времени 0 <I t ^

^ Э, а другое с убывающей X2{t).

Конкретный вид функций Xx(t)

и X.2(t) может быть любым, например, как на рис. 20. Сравнивая

функции надежности P(t)

для этих двух изделий,

вычисленные по

точной формуле, полученной из выражения (6.13),

 

Р (t)

= ехр {— lX{t)di)

(6.19)

 

о

 

ипо формуле (6.15), можно сделать следующие выводы:

1.Как и следовало ожидать (см. рис. 20, а), в момент t = й точные выражения Р,(^) и P2(i) соответственно для первого и вто­

рого изделий равны значению P(t) — ехр(—Xgt), где XQ— средняя на интервале времени 6 величина параметра потока отказов.

2. Для всех моментов времени t, не выходящих за интервал О

(на котором, найдено среднее значение XQ), истинная

функция на­

дежности для изделий с возрастающей X(t)

 

Р (0 > ехр (—Хе(),

(6.20)

для изделий с убывающей X{t)

Р

( 0 < ехр

(6.21)

Таким образом, для t <

0 вычисление Р {() —

ехр (— Xt) будет

с запасом при возрастающей X{t),n не будет в запас при убывающей

Щ.

Средняя наработка на отказ. Вероятность безотказной работы (6.16), или функция надежности P(t), однозначно связана с законом распределения случайной величины т — времени работы изделия до отказа (между двумя соседними отказами). В соответствии с вы­ ражением (6.1) имеем

 

Р (/) =

Р {т >

t) = 1 — Р {т < і) = 1 — F (0,

(6.22) *

где F(i) — интегральный

закон

(2.1) распределения

случайной ве­

личины т>, называемой по ГОСТ

13377—67 [47] наработкой

на от­

каз. Следовательно,

из

выражений (6.22)

и (6.16) F(t) — 1 —

— ехр

{—yd).

Отсюда

дифференциальный

закон для т в

соответ­

ствии

с формулой

(2.3)

 

 

 

 

 

 

 

f (t)

=

F'

(0

=

X ехр {-Xt);

0 < t <

оо.

(6.23)

Получили экспоненциальную плотность распределения (3.4). По­

этому функцию (6.16) называют

э к с п о н е н ц и а л ь н ы м

з а-

к о н о м н а д е ж н о е т и.

 

 

Зная f(t) для реакторной установки, легко найти среднее время

ее безотказной работы по формуле (2.10)

 

оо

оо

 

Т0 = М (т) = оf tX ехр {—Xt)dt = —t ехр {—Xt)оI - f

 

со

со

 

+ [ ехр (—Xt) dt

= J Р (0 dt = 1/1.

(6.24)

6

о

 

В работе [7] показано, что если для изделия, имеющего возра­ стающую Х{(), известна истинная средняя наработка на отказ Т0, то, принимая в качестве функции надежности изделия экспонен­ циальный закон с X = 1/7"0 — const

 

 

Р {t) =

ехр (—t/T0 ),

(6.25)

занизим фактическую величину. Р(^) для всех

моментов времени

t <С Т0

(т. е. получим результаты с запасом).

.

г

Интенсивность отказов

неремонтируемых

изделий. Такие

из­

делия

работают

только до первого отказа и далее заменяются

на

новые

(потока

отказов нет, отказ всего один). Их надежность ха­

рактеризуют интенсивностью отказов. Ее часто обозначают той же буквой X, что и параметр потока отказов для ремонтируемых из­ делий, так как она аналогично ведет себя во времени (см. рис. 19)

Соседние файлы в папке книги из ГПНТБ