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

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

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

А.И. КЛЕМИН

ИНЖЕНЕРНЫЕ ВЕРОЯТНОСТНЫЕ РАСЧЕТЫ ПРИ ПРОЕКТИРОВАНИИ ЯДЕРНЫХ РЕАКТОРОВ

0

А Т О М И З Д А Т М О С К В А 1973

У Д К 621 . 039 . 5:519 . 21

1. . . I ,•

Клемин А . И . И н ж е н е р н ы е в е р о я т н о с т н ы е р а с ч е т ы п р и п р о е к т и р о в а н и и я д е р н ы х р е а к т о р о в . М . , А т о м ­

и з д а т , 1973 г., 304

с.

К н и г а п р е д с т а в

л я е т собой р у к о в о д с т в о по р е ш е н и ю

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

н и я и з д е л и й

р е а к т о р о с т р о е и и я ,

в ы б о р у

р е ш е н и й

в у с л о в и я х н е о п р е д е л е н н о с т и .

 

 

Р и с у н к о в

32, т а б л и ц 53, б и б л и о г р а ф и я

131 н а з в а ­

н и е .

 

 

 

.03315 — 113

К, j47—72 Атомиздат, 1973 034(01)—73'

ПРЕДИСЛОВИЕ

Широкое развитие атомной энергетики в нашей стране выдви­ гает новые задачи перед создателями ядерных реакторов, в первую очередь, перед конструкторами, а также перед всеми специалиста­ ми, связанными с проектированием, изготовлением и эксплуата­ цией ядерных энергетических установок. Среди этих задач можно выделить в качестве главных следующие: 1) повышение научнотехнического уровня разработок; 2) сокращение сроков проектных работ; 3) создание экономичных, надежных и безопасных в эксплу­ атации ядерных реакторов. Решение их немыслимо без совершен­ ствования инженерных расчетов ядерных реакторов, без широкого внедрения в практику проектных и исследовательских работ новых математических методов, в частности, вероятностно-статистических методов современной теории надежности, теории планирования эксперимента, теории ошибок, статистических решений, теории игр и др.

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

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

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

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

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

Книга содержит много конкретных примеров, доводимых до численных результатов, которые иллюстрируют применение изла­ гаемых методов в реакторостроении. Некоторые примеры имеют

самостоятельное практическое значение. В основу книги положены

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

исследователей,

монографии последних лет по теории вероятностей,

математической

статистике, планированию эксперимента, теории

надежности, игр

и т . п., а также опыт работы автора по данной тематике.

В книге § 5.4 и 11.2 написаны В. Н. Орловым и А. И. Клеминым, глава 8 и таблицы приложения П. 10, П. 11 — Е. Ф. Поляковым и А. И. Клеминым, гл. 9, 13, 14, табл. П. 1—П. 9, П. 12 и предмет­ ный указатель—Л. Г. Клеминой и А. И. Клеминым. На стр. 72, 77 — 80 книги использованы материалы, любезно предостав­ ленные Г. Е. Романцовым, за что автор выражает ему благодар­ ность.

Автор глубоко признателен академику Н. А. Доллежалю, одоб­ рившему идею написания настоящей книги и принявшему участие в обсуждении ее плана. Автор особо признателен проф. И. Я- Емель­ янову, просмотревшему рукопись и сделавшему ряд ценных заме­ чаний, проф. С. М. Фейнбергу, взявшему на себя труд рецензиро­ вания рукописи, а также докт. техн. наук И. С. Коченову и канд. техн. наук М. М. Стригулину за постоянный интерес к работе. Автор считает своим приятным долгом выразить благодарность канд. техн. наук Л. Г. Клеминой и ст. технику Т. Г. Коробовой, оказав­ шим большую техническую помощь в подготовке рукописи, а так­ же Л. И. Скоромной, Л. М. Бортниковой, Е. А. Шиверскому,

В.С. Бутнику и Н. А. Ковальчук за участие в оформлении работы.

А.И. КЛЕМИН.

Раздел I . ЭЛЕМЕНТЫ ТЕОРИИ ВЕРОЯТНОСТЕЙ И МАТЕМАТИЧЕСКОЙ СТАТИСТИКИ, НЕОБХОДИМЫЕ ДЛЯ ИНЖЕНЕРНЫХ ПРИЛОЖЕНИЙ В РЕАКТОРОСТРОЕНИИ

Г л а в а 1.

М Е Т О Д Ы В Ы Ч И С Л Е Н И Я В Е Р О Я Т Н О С Т Е Й

СО Б Ы Т И Й

§1.1. Прямые методы

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

название

в е р о я т н о с т и . Можно

сказать,

что

вероятность

Р{А]

есть

мера степени объективной возможности события

А.

Ее

условились измерять в долях единицы:

1 > Р { Л )

>

О, Р{А}

=

О,

если

А — невозможное событие, Р{Л} =

1, если

А —достоверное

событие.

 

 

 

 

 

 

Все события в природе можно разбить на два

класса: 1)

собы­

тия, вероятности которых вычисляются точно; 2) события, вероят­ ности которых могут быть вычислены только приближенно.

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

тремя свойствами: несовместностью,

равновозможностью и образо­

ванием полной группы. Н е с о в м е

с т н ы м и называются такие

события, которые не могут появиться вместе в результате одного опыта. Наступление одного из них исключает возможность наступ­ ления другого. Р а в н о в о з м о ж н ы м и (или равновероят­ ными) называются события, обладающие объективно одинаковыми возможностями для осуществления в результате конкретного опыта. П о л н у ю г р у п п у образуют события, связанные с конкрет­ ным опытом, из которых хотя бы одно обязательно происходит в ре­ зультате этого опыта. Классическими примерами событий, обла­ дающих тремя приведенными свойствами, являются: выпадение герба и выпадение цифры при подбрасывании правильной монеты; появление 1, 2, 3, 4, 5, 6 очков при бросании игральной кости (симметричного кубика).

События, обладающие тремя упомянутыми свойствами, называ­ ются в теории вероятностей с л у ч а я м и или шансами, а про

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

задачи, связанные с такими событиями,

говорят,

что они сводятся

к схеме случаев II—3]. Вероятность таких событий

вычисляется

по формуле

 

 

 

 

Р{А}

= т/п,

 

 

(1.1)

где-гс — общее число случаев; т — число

случаев, благоприятных

событию А. Случай называют

б л а ' г о п р и я т н

ы м

некоторому

событию, если появление случая означает осуществление этого со­

бытия. Например, событию А — появлению

нечетной цифры

при

бросании

игральной кости — благоприятны

три случая: выпадения

очков 1, 3 и 5; остальные три случая

(2,

4

и 6)

неблагоприятны в

соответствии с формулой (1.1) Р[А}

=

3 / 0

=

72 .

 

п и т в

кон­

Трудности, возникающие при определении чисел

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

фор­

мулами

комбинаторики — раздела

математики,

в

котором

изу­

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

тов [4,

5].

Формула (1.1) имеет геометрическую модификацию— формулу,

которая

позволяет вычислять так называемые г е о м е т р и ч е ­

с к и е

в е р о я т н о с т и . Если событие А есть попадание точки

в результате случайного бросания ее в область площадью 5 на часть

этой области, имеющую площадь SA, то

 

Р {А} = SA/S.

(1.2)

Аналогичные формулы можно записать для одномерной и трех­ мерной областей.

Для вычисления вероятностей событий, которые не укладываются

в схему случаев (принадлежат ко второму классу

событий), необ­

ходимо организовать серию испытаний. Допустим,

необходимо вы­

числить вероятность события А,

которое либо происходит, либо нет

в результате некоторого опыта.

С этой целью данный опыт необхо­

димо повторить п раз и найти частоту события А

в проведенной се­

рии опытов по формуле

 

 

 

 

Р с т

{ Л }

= т/п,

(1.3)

где т — число опытов, в которых

появилось событие А; п — пол­

ное количество опытов. Эту

частоту

называют

с т а т и с т и ч е ­

с к о й в е р о я т н о с т ь ю

события

А.

 

Дальнейшие рассуждения базируются на законе больших чи­ сел, который в форме Бернулли гласит: если Р{Л) —вероятность

появления

события А

в результате некоторого опыта, то как

бы ни

было мало

число є >

0, при п. -»- оо

 

 

Р { | Р с т { Л } - Р { Л } | > е } - > 0 ,

( ] 4 )

т. е. при неограниченном увеличении числа независимых опытов частота события будет сколь угодно мало отличаться от его вероят­ ности: Р[А] ^ Р С т И ) -

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

§ 1.2. Вычисление вероятностей событий

с помощью теорем сложения и умножения

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

Часто на практике вероятностьпоявления события А зависит от того, произошло или нет некоторое другое событие В. Такие со­

бытия

называются

з а в и с и м ы м и .

Для учета зависимости со­

бытий

вводится понятие у с л о в н о й

в е р о я т н о с т и

собы­

тия А при наличии В Р{А/В]

или вероятности события А при усло­

вии, что событие

В произошло. Величину

Р{А)

иногда

называют

б е з у с л о в н о й

в е р о я т н о с т ь ю .

Для

зависимых

собы­

тий А и В

 

 

 

 

 

 

 

 

 

 

 

 

 

Р{А)фР{А1В),

 

 

 

(1.5)

если же они независимы, то

 

 

 

 

 

 

 

 

 

 

Р {А} = Р {А/В}. .

 

 

 

(1.6)

Формулы

для вычисления

вероятностей

суммы и произведения

двух событий

приведены в табл. 1.1.

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

1.1

 

Теоремы сложения и умножения вероятностей для двух событий

 

События

 

 

Зависимые

 

 

Независимые

 

С о в м е с т н ы е

 

 

 

 

 

 

 

 

 

 

 

 

 

Р{А+В)=Р{А)

+

Р ( Л + В } = Р { Л } +

 

 

 

 

+

Р{В}—Р{А-В]

 

+ Р І В ) —

Р{А-В]

 

 

 

 

Р{А-В}=Р{А}.р{В/А}

 

Р{А-В}=Р{А}-Р{В]

 

Несовместные

 

 

 

 

 

 

 

 

 

•® ©

Р{А +

В)=Р{А)+Р{В)

Т а к и х событий

нет

 

Р{А-В)=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для п независимых событий

РІА.-А,-...

п)

=

П Р { Л , } ,

(1.7)

 

 

 

 

(=1

 

для п несовместных

событий

 

 

 

 

Р{А,+

А2+

... +

Ап)

= S P { A - } .

(1.8)

Формула полной вероятности. Пусть в результате некоторого опыта обязательно наступает только одно (заранее неизвестно какое) из событий ННО, НП (иными словами, имеем несовместных событий, образующих полную группу). В свою очередь, совместно с каждым из этих событий Н( может произойти интересующее нас событие А с вероятностью Р{ А/Я; \. Тогда вероятность, что событие А произойдет, может быть вычислена по формуле полной вероятности:

 

 

Р{А}

=

f Р {Ht}

Р {АІНі).

(1.9)

Например, если Я х — работа реактора

на мощности, Я 2 — простой

реактора, А — повышение

активности

в реакторном зале выше не­

которого-уровня,

причем

известно, что P ( # j } = 0 , 9 ;

Р ( Я 2 } = 0 , Г ,

РІА/Н^

= 0,01;

Р\А/Н2]

=

0,005, то вероятность события А для

такого

реактора

 

 

 

 

 

Р {А} = Р Г} Р {А1НЛ}

+ Р {Я2 } Р {Л/Я,} =

0,0095.

Теорема гипотез (формула Бейеса). Вероятность гипотезы* Ht

при условии, что событие А произошло (или вероятность, что собы­ тие А произошло именно совместно с реализацией гипотезы Я г ) , записывают в виде

P{Ht/A} = * { Н і ] Р { А / Н і } ,

(1.10)

где #!, Я 2 , Я п — как и в формуле (1.9), взаимоисключающие события, одно из которых обязательно реализуется.

*

З д е с ь под г и п о т е з о й НІ п о н и м а е т с я п р е д п о л о ж е н и е о т н о с и т е л ь н о п о я в ­

л е н и я

с о б ы т и я Ні.

Г л а в а 2.

Ф У Н К Ц И О Н А Л Ь Н Ы Е И Ч И С Л О В Ы Е Х А Р А К Т Е Р И С Т И К И С Л У Ч А Й Н Ы Х В Е Л И Ч И Н

§2.1. Законы распределения случайной величины

ивычисление вероятности попадания ее

взаданный интервал

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

н о й, например

число остановок реактора АЭС в течение года. Ко­

гда возможные

значения случайной величины непрерывно запол­

няют один (в том числе бесконечный) или несколько промежутков

на числовой оси, ее называют

н е п р е р ы в н о й ,

например по­

грешность эксперимента.

 

 

 

 

 

 

И н т е г р а л ь н ы м

з а к о н о м

(или функцией)

распреде­

ления случайной величины X называется

функция

F(x),

численно

равная вероятности, что

случайная

величина X примет

значение

меньшее, чем х:

 

 

 

 

 

 

F

(х) =

Р{Х<х);

 

 

О < F (х) < 1, F (—

оо)

= О,

F (оо) =

1.

 

Интегральный закон несложно найти опытным путем. Для этого необходимо провести п независимых наблюдений над случайной ве­ личиной^ и получитьп ее возможных значений хх2, Если обозначить т{х) число тех наблюдений из п, в которых величина X оказалась меньше некоторого фиксированного х, то в силу опре­ деления вероятности (1.3) и функции (2.1) эмпирический интеграль­ ный закон можно записать в виде

 

 

 

F3 (х)

= т (х)/п.

 

 

 

(2.2)

Для

дискретных

случайных

 

величин F(x)

представляет

собой

ступенчатую функцию.

Для

непрерывных

случайных

величин-

F(x) — неубывающая функция,

как правило, монотонно

возраста­

ющая в области возможных значений.

 

 

 

 

 

Фиксированное значение хр

 

случайной величины

X,

для

кото­

рого F(xp) = р, в теории вероятностей называют к в а н т и л ь

ю.

Стало быть, квантиль хр,

отвечающая вероятности р,

есть

значе­

ние X,

вероятность

попадания

левее которого Р{Х

хр]

=

р.

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