Файл: Клемин А.И. Инженерные вероятностные расчеты при проектировании ядерных реакторов.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 11.04.2024

Просмотров: 267

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

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

отказа, обычно

мала. Так, если

вероятность, что каждый из них

(в отдельности)

произойдет в

течение

времени

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

Т'эиспл = 3 0

Л е т >

Р а в н а

Р 1 (Тш<спл) =

Р 2

(^экспл) =

Ю" 3 ,

ТО ВЄ-

роятность,

что

они оба

произойдут за

время Твкйпл,

Р Х Р 2

= Ю - 6 .

Однако вероятность, что эти два отказа произойдут одновременно, совпадут, из выражений (1.2) и (1.7) равна

•Pi (^экспл) Р 2 (Т'эксил) ' зо . 8740 '

где Ат — время, характеризующее длительность протекания и раз­ вития первого отказа, я. Если, например, Дт « 1 ч, то

Pi (Пкспл) Р-2 (Пкспл) ( / 3 0 • 8740) = 0,4 • 10-",

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

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

Проблема исходных данных. Перед инженером, оценивающим радиационную безопасность АЭС, прежде всего возникает вопрос, как вычислить (задать) вероятности отказов элементов реакторной установки в условиях недостаточности эксплуатационных статисти­ ческих данных по таким отказам. Эта задача не такая неразреши­ мая, как может показаться. Во-первых, из перечисленных выше групп оборудования АЭС наиболее специфичными являются группы 1а, 2 и 3. Остальные устройства (как, впрочем, и упомянутые) вклю­ чают в себя оборудование, довольно широко используемое в про­ мышленных установках других отраслей техники, по которому имеется достаточное количество статистических данных (см. на­ пример, табл. П. 10.в приложении). И хотя на некоторое стандартноеоборудование в условиях применения на АЭС влияет фактор облу­ чения, это еще не означает, что опытом, накопленным при эксплуа-


тацни стандартного оборудования в обычных условиях, нельзя вос­ пользоваться. Во-вторых, к настоящему времени у нас и за рубежом уже накоплен немалый опыт эксплуатации ядерных реакторных уста­ новок п собраны определенные статистические данные, в том числе по оборудованию групп 1а, 2и 3 (см. табл. П.10, П.11). Кроме того, часто вероятности отказов соответствующих элементов оборудова­ ния АЭС можно найти, используя формулы и методы для расчета

.надежности, изложенные в гл. 6.

Заметим, что во всех случаях, когда отсутствуют достаточно надежные статистические данные по отказам, оценки радиационной безопасности АЭС могут и должны производиться в запас.

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

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

же табл. П.10

и П.11). Пусть эта вероятность за время эксплуа­

тации Тъкаия

равна рх. Разрыв трубопровода приведет к тому, что

за некоторое время первый контур потеряет практически весь теп­ лоноситель. Поведение определяющего параметра канала (допустим, температуры оболочки твэла) будет зависеть от готовности устройств защиты и, в частности, от того, будут они исправны или откажут (не сработают). Для анализа возможных в данном случае ситуаций удоб­ но построить цепочки предполагаемых событий (первичного, вто­ ричного, третичного и т. д. отказов оборудования) и представить их в форме табл. 8.1. Вероятности различных состояний систем ава­ рийной остановки и охлаждения можно оценить по соответствующим статистическим данным (см., например, табл. П. 11) или из вероят­ ностного анализа схем и элементов этих систем. Заметим, что реше­ ние последней задачи облегчается, так как по отказам элементов радиоэлектроники и автоматики накоплен достаточно большой ста­ тистический материал (см., например, [50]). Полученные вероят­ ности вносим в табл. 8.1.

Для оценки количества радиоактивного вещества, выделяюще­ гося при каждой из рассмотренных восьми возможных аварий (см. табл. 8.1), предварительно необходимо найти среднее число ка­ налов активной зоны, в которых произойдет пережог оболочек твэлов. По средним значениям остаточного тепловыделения в канале, средним значениям проходных сечений и т. д. находим зависимость определяющего параметра канала от времени начиная с момента разрыва трубопровода Г о б . й = Т о б . ) г {і). Возможные случайные от-


Первичный

отказ

Р а з р ы в т р у ­ бопровода первого кон ­ тура (веро­ ятность

этого отка­ за P l )

последую­Номерсобытиящего

1

2

3

4

5

6

7

8

 

 

Х а р а к т е р и с т и к и в о з м о ж н ы х а в а р и й н ы х с и т у а ц и й *

 

Система аварийной

Система аварий­

Система аварийно­

Система локали­

 

зации и подавле­

 

остановки I

ной остановки II

го охлаждения

ния активности

Полная вероятность

Признак состояния

 

Признак состояния

 

Признак состояния

 

Признак состоянпя

 

Вероят­

Вероят­

Вероят­

Вероят­

аварии

 

ность

 

ность

 

ность

 

ность

 

 

состояния

 

состояния

 

состояния

 

состояния

 

+1 — P 2 I

+1 — P 2 I

+

+!—/>21

Р21

Р%\

Р21

Р21

+

+

+

+

1

+

1

+

1

 

1

 

1—/=*211

+

1 — P 2 I I

+

1—Ра п

1~Р2П

1—Рз

1 - Р з

Рз

Рз

! — Р з

1—Рз

Рз

Рз

+

+

+

+

1— Pi

Pi

1— Pi

Pi

I—Pi

Pi

1—Pi

Pi

P i l l — p 2 1 ) (1—Рз) (1—Pi)

P i ( l — P21)

0 — Р з )

Pi

P i O — P21)

Р з ( 1 —

Pi)

Pi ( I — P21) Рз Pi

P1P21 (I—Р211) 0 — Рз) 0 — Pi)

P1.P21 U — P211) (1—Рз) Pi

P1P21 ( 1 — Р 2

ц ) Р з ( 1 — Pi)

P1P21

P211) Р з Р «

Величина выбро­ са Я, кюри

Еі

Ei

Ея

Et

Еь

Еа

Еп

Е%

* Признаки

состоянпя означают следующее: -|

система функционирует нормально;

отказ системы;

состояние системы

в данной ситуации

безразлично для безопасности.

 

 

 


клоненая от средних значений приводят к разбросу значений опреде­ ляющего параметра. Предположим, что закон распределения слу­ чайной величины Г о С . ; і (/)—температуры оболочки в фиксиро­ ванный момент времени t — является нормальным законом, это достаточно близко к истине [61. Вводя определяющую функцию

 

ПО

Z

 

 

 

для /г-й группы

каналов

(работающих

в

одинаковых

режимах)

со среднеквадратичным отклонением а ([,,) =

а Л , по аналогии с вы­

ражением (6.115)

получаем

вероятность

для определяющего пара­

метра выйти в отрицательную область значений (пережог

оболочки)

Р Г = 0 , 5 - Ф (/Г/од).

Тогда среднее количество отказавших каналов /г-й группы (мате­ матическое ожидание) по формуле (3.61)

mh = nh Р* ,

где nk — количество каналов в /г-й группе. Отсюда среднее коли­ чество каналов, в которых произойдет пережог оболочек твэлов, во всей зоне составит

 

 

 

k*

 

т=

Т

m.h=

Ъ

%Pfc

(8-1)

 

А =

1

к=

I

 

(k* — полное число групп

каналов).

 

 

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

ные программы,

позволяющие оценивать

долю 1 3 1 1 (или любого

другого продукта

деления), находящегося

в свободном состоянии

(т. е. в пространстве между оболочкой и топливом и в специальных

пустотах в объеме твэла). Так, по оценке английских специалистов [72[, сделанной для усовершенствованного газоохлаждаемого реак­ тора AGR (топливо U0 2 ), при температуре топлива не больше 1600° С

доля свободного 1 3 1 1 от полного его количества в наиболее

горячих

каналах составляет 1/100, а

в среднем по активной зоне

1/1000.

Основная часть йода удерживается кристаллической решеткой U 0 2 .

В других активных зонах, использующих то же топливо,

но при

более высокой

температуре,

доля свободного

йода

несколько

больше.

 

 

 

 

 

Следует также

учитывать,

что в отдельном отказавшем

канале,

.содержащем несколько твэлов, может произойти

пережог

оболочек


не всех твэлов. Иногда в запас считают, что в отказавшем канале все твэлы имеют нарушение герметичности оболочек.

Пусть после проведенного анализа принято, что величина вы­ хода 1 3 1 1 из отказавшего канала составляете. Тогда, используя фор­ мулу (8.1), можно записать:

 

 

(8.2)

где Е'

— полная

активность ш 1 , выходящего из поврежденных ка­

налов

и первого

контура в реакторное помещение. Символом Е

(в отличие от Е')

будем в дальнейшем обозначать активность радио­

активного вещества, выброшенного в окружающую среду. Эта ве­ личина находится после проведения аналогичного описанному ана­ лиза надежности и эффективности работы устройств локализации радиоактивного вещества на АЭС. При этом для случаев первичных •отказов, подобных разбираемому (нарушение целостности контура высокого давления), необходимо учитывать повышение давления в реакторном помещении и возможность его разгерметизации или других нежелательных явлений. Иными словами, следует учитывать, что отказ системы локализации может быть событием, зависящим от первичного отказа. Полученные в результате анализа величины выбросов Е заносим в табл. 8.1.

Аналогичное рассмотрение должно быть проведено для всех воз­ можных первичных отказов. Полученные в результате такого ана­ лиза таблицы типа 8.1, содержащие вероятности аварийных ситуа­ ций и величины выбросов, будут использоваться на следующих ста­ диях оценки безопасности: при расчете распространения радиоак­ тивного вещества (§ 8.4) и вычисления величин рисков (§ 8.6).

Итак, на первом этапе оценки радиационной безопасности тре­ буется комплексное изучение очень многих вопросов: физических, теплогидравлических, металловедческих и т. д. в их тесной взаимо­ связи и с учетом случайной природы отказов. Поэтому 'оценка ве­ роятностей тех/или иных аварийных ситуаций должна проводиться комплексно с участием различных специалистов, и в том числе спе­ циалиста по количественной оценке надежности. На пути решения этой задачи не существует непреодолимых трудностей. Последова­ тельное рассмотрение, проводимое для конкретной АЭС, обнаружит, что на настоящей стадии развития ядерной энергетики возможно оп­ ределить упомянутые вероятности, основываясь на уже имеющихся статистических данных (см. § 8.3), разработанных математических моделях отказов, а также на принимаемых в запас предположениях.

§ 8.3. Исходные статистические данные по отказам оборудования АЭС

Ниже будут приведены статистические данные по отказам обору­ дования и вероятностям возникновения различных аварийных ситуа­ ций на АЭС, которые рекомендуются в зарубежных работах [42,