Файл: Уриг, Р. Статистические методы в физике ядерных реакторов.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 16.10.2024
Просмотров: 111
Скачиваний: 0
при Д = t. Величина S является мощностью нейтронного источника в подкритической системе. Два вида уравнений (3.69) записаны для подкритической и критической систем соответственно.
Уравнение (3.66) имеет некоторые особенности, представляющие интерес. Если положить е < р£, т. е. у -v 1, и увеличить мощность источника S так, что С3 (/) -> 0, то в результате получится выражеНие
рсс (t)dt = Ее ехр (— Fet) dt, |
(3.70) |
которое идентично уравнению (3.64) для пуассоновского распре деления отсчетов, т. е. некоррелируемого процесса.
Вероятность рсс (0 зависит как от уровня мощности (мощности источника), так и от эффективности детектора, но ее использование имеет преимущества перед другими статистическими методами. При высоких уровнях мощности, когда метод росси-альфа неприме ним для измерения параметров, С2->- 0 , получаем выражение
Рсс (t)dt = Сг {t) dt, |
|
(3.71) |
которое может быть использовано для измерений параметров. |
||
Если эффективность детектора очень низка (т. е. е |
р£ и у |
1) |
все методы счета слабых интенсивностей непригодны. Однако, если 2eDv/Pp С 1, уравнение (3.65) записывается в виде
Рсс (t) — ехр |
(—Fet) [Ее + |
(eDJ2Ap) ехр |
(—ocf)] = |
= ехр |
(—At) [А + |
В ехр (— at)], |
(3.72) |
которое, исключая член ехр (—Fet) сопоставимо с выражением рос си-альфа. Для быстрых реакторов, когда эффективность детекторов низка, низка скорость счета и временные интервалы коротки, эк споненциальный член ехр (—Fet) примерно равен единице и урав нение (3.72) переходит в уравнение
Рсс (0 dt = (А + В ехр (—at))dt, |
(3.73) |
идентичное уравнению (3.20) для эксперимента росси-альфа. Это объясняет успешное применение метода росси-альфа Брансоном и др. [16], когда они фактически измеряли счетно-счетные времена.
Экспериментальная методика была описана в § 3.3. Первый им пульс запускает анализатор, каналы которого подключаются пре цизионным таймером (отметчиком времени). Второй отсчет останав ливает анализатор, отсчет помещается в ячейку памяти соответст вующего канала, на котором анализатор был остановлен, и систе ма восстанавливается для следующего отсчета. Такое устройство регистрирует только половину данных, т. е. интервал времени меж ду двумя разными импульсами. Если анализатор автоматически за пускается схемой остановки и восстановления, то могут быть за регистрированы все данные. Однако этот прием укорачивает интер вал времени на величину, равную мертвому времени (На время, не обходимое для остановки анализатора, запоминания отсчета, вос
51
становления и пуска анализатора). В современных анализаторах это мертвое время может быть очень коротким, однако, как обычно, необходимо сделать соответствующую поправку.
Метод распределения случайно начатых интервалов. Метод рас пределения случайно начатых интервалов близок методу распре деления счетно-счетных интервалов. Главное отличие его состоит в том, что начало интервала выбирается случайно процессом, ко торый не коррелирует с изучаемым ядерным явлением. Бабала [5] получил выражение для вероятности распределения случайно нача тых интервалов:
(V+ 1) + (Y 1) exp ( ayl) 1 Pro (0 * = 2FePo (О (V+1)2—(У—1)2ехр(—ayt) J
где у и р0 (0 определяются уравнениями (3.69) и (3.54) соответст венно. Как и в случае распределения счетно-счетных интервалов, процесс становится пуассоновским и
Pro (t) dt = Fe exp (—Fst) dt, |
(3.75) |
когда у -v 1 из-за уменьшения эффективности e или значительной
подкритичности системы.
(
Пачилио [32] отметил, |
что J p«o(0 dt представляет вероятность |
того, что после времени |
о |
t = 0 , выбранного случайно, первый им |
пульс придет в интервале от 0 до t. Поскольку р0 (0 есть вероят ность того, что то же самое событие произойдет в интервале между t и бесконечностью, то получается выражение
Po(t) + {pRo(t)dl=h |
(3-76) |
о |
|
из которого следует, что |
|
Pro (0 = — др0 (t) /dt. |
(3.77) |
Мы также можем считать pRo(t) вероятностью того, что после импульса, пришедшего в интервале между произвольным началом и dt последует пустой интервал. Эта вероятность может быть выраже на как произведение вероятности одного отсчета между 0 и dt и веро ятности следующего отсчета, происходящего за пределами времени, большего t. Отсюда
i
Pro W dt = (F*dt) |
j) Pcc W > |
(3.78) |
|
0 |
|
где интеграл дает вероятность того, что после отсчета при t = О следующий отсчет произойдет в интервале между 0 и /. Следова-
тельно,
1 &Pro (О |
дг Ро (О |
(3.79) |
Рсс М — Ге Jt |
Fe дР |
|
52
Экспериментальный подход по существу тот же, что и в методе рас пределения счетно-счетного интервала, за исключением того, что анализатор после восстановления запускается случайным импуль сом. Этот метод эффективен при изучении систем с тепловыми ней тронами, но имеет существенные ограничения и не может применять ся при изучении систем на быстрых нейтронах. Остин и др. [36] использовали этот метод, который они назвали «выжидающий вре менной альфа-метод», и получили хорошее согласие с измерениями росси-альфа и измерениями с импульсным нейтронным источником.
§ 3.8. Метод мертвого времени (метод Сринивасана)
Другой метод измерения а был предложен в работе Сринива сана 137]. Он основывается на том факте, что при введении искусст венно изменяемого мертвого времени в измерительную аппаратуру оно оказывает влияние на корреляцию между отсчетами. Обсудим это влияние для случая нечувствительности прибора, определенной Сринивасаном следующим образом. Предположим, что последова тельность входных импульсов (реальных отсчетов) от нейтронного детектора поступает в прибор, который дает последовательность выходных импульсов (выходных отсчетов). Если прибор передал реальный отсчет на выход, то появляется запрет на второй выходной сигнал, если прошел временной интервал, меньший d (мертвого вре мени) между двумя последующими реальными отсчетами. Таким образом, этот прибор регистрирует число интервалов длиннее, чем d между реальными отсчетами.
Для некоррелированных отсчетов отношение между скоростью счета Cdна выходе нечувствительного прибора и реальной скоростью
С дается формулой |
|
Cd = С ехр (—Cd), |
(3.80) |
где экспоненциальная функция является просто вероятностью того» что интервал между двумя реальными отсчетами будет больше» чем d. Дисперсия выходных отсчетов за время А для такой системы при пуассоновском распределении входных импульсов по Сринивасану определяется следующим образом:
сг = с + |
с2 ( A _d ] \ |
(A |
d )> 0. |
(3.81) |
||
откуда |
|
|
|
|
|
|
с5—? |
. |
- г , |
/ Д—d у |
(3.82) |
||
-с |
1 |
Т |
( А |
) } |
||
|
||||||
где |
|
|
|
|
|
с — C dА. |
(3.83) |
Для коррелированных отсчетов, т. е. таких, которые проис ходят в ядерном реакторе нулевой мощности, дисперсия отсчетов для выключающегося прибора равна
|
|
|
с2 = с + с2 |
Д—d \ ~n+ |
|
|
|
|||
|
+ - а д |
А—d \ |
1—ехр [— а (Д—d)] |
exp (—ad), |
(3.84) |
|||||
|
|
|
Д У1 |
|
|
а (Д■—d) |
|
|
|
|
где для |
мертвых |
времен короткой продолжительности |
|
|||||||
|
|
|
В = 1—Cd |
|
Dv |
у 3' |
|
|
(3.85) |
|
|
|
|
1 + vPy |
/ У— |
|
|
||||
|
|
|
|
|
|
2SA |
VY+1 ' |
|
|
|
Преобразование |
уравнения |
(3.84) |
дает |
|
|
|
|
|||
|
|
ГИ— “С2 |
1 |
|
Д—d\- |
|
|
|
||
|
|
|
1— С |
|
|
|
|
|
|
|
х — |
а |
1—ехр [—а ( Д —d)] |е х р ( — a d ) = l — с 1 — |
-ь |
|||||||
|
|
|
а (Д—d) |
|
|
|
|
|
|
|
f |
Рсо (d) Вех р (—ad) ■ |
{1 |
1—ехр [ —а(Д—d)] |
(3.86) |
||||||
а (Д—d) |
|
|||||||||
|
|
|
|
|
|
|
|
|
||
где pco(d)— вероятность того, |
что интервал |
между двумя |
отсчет |
|||||||
тами будет |
больше, чем |
|
|
|
|
|
|
|
||
|
Pco(d) = 2Po (d) |
(V + 1) + (V— 1)ехр(— ayd) |
1 |
(3.87) |
||||||
|
(7+ l)2—(Y—1)2ехр (— ayd) |
J’ |
||||||||
|
|
причем p0 (d) вычисляется по формуле (3.54) при A = d, а остальные обозначения прежние.
Уравнение (3.86) является довольно сложным для практического определения а при вариации мертвого времени d. Однако оно может быть использовано для оценки влияния мертвого времени прибора на дисперсию. Легко увидеть, например, что влияние мертвого вре мени может быть незначительным, если Cd 1, ad < 1 и d « А.
Подобный анализ вюлучае невыключающегося прибора является более сложным и здесь не приводится.
§ 3.9. Методы корреляционного анализа
Взаимная корреляционная функция стационарного процесса cpKI/ (т), частным случаем которой является автокорреляционная функция (т) при х = у, определяется уравнениями (2.74) и (2.77):
00 |
оо |
фж„ = Е [х (t) y(t + т)] = ^ |
^ X (ti) У (t2) Р [* (к), у (/2)] dxdy, (3.88) |
54
где р [х (У , у (/2)] — функция совместной вероятности наступления события х в момент времени tx и события у в момент t2при определен ном
т = / 2 — tx. |
(3.89) |
Если событие х — регистрация нейтрона детектором 1 и у — де тектором 2 (или детектором 1 в частном случае, когда х — у), то легко показать, что корреляционная функция будет вероятностью пары отсчетов, происходящих в интервале Аг при tx и Л2 при t2, т. е. происходящих с интервалом т. Эта же величина исследовалась в § 3.3 при обсуждении метода росси-альфа. Отсюда
Ф*» 00 = РИъ t2) = Рс *s) + Pr (*i>*2). (3-90)
где р (tlt ■— вероятность отсчета в момент tlt за которым следует отсчет при /2; индексы С и R относятся к коррелированным и случай ным событиям. Если используется один детектор, то уравнение (3.90) становится подобным уравнению (3.17) при условии, если исключить
дельта-функцию Дирака при т |
= 0: |
|
|
Фта (т) = Р е2 |
Fe2 — |
— |
ехр (— ат) +F&б (т) = |
' |
2 (1 |
—/гр) 1 |
^ |
= А2 + АВ ехр( —ат) + Л8 (т) = А [Л + Д ехр(— ат)] + Л8 (т), (3.91)
где Л и В определяются уравнениями (3.25) и (3.26). Этот подход известен как метод автокорреляционного анализа. С теоретической точки зрения, метод по существу является методом росси-альфа, однако методика измерения совершенно другая. Следует отметить, что случайные или фоновые величины зависят от квадрата скорости делений F2(или мощности), в то время как амплитуда перед экспо нентой зависит только от F. Следовательно, этот метод ограничивает ся очень низкими скоростями делений. Член с дельта-функцией Ди рака отсутствует в экспериментах росси-альфа вследствие запазды вания по фронту первого канала совпадений (см. рис. 3.2).
При использовании двух детекторов одинаковой эффективности е накапливаемые случайные отсчеты независимы и, следовательно, некоррелированны, так как детектируемые нейтроны поглощаются, и уравнение (3.91) принимает вид:
Fe2 D ft2 |
|
|
Фжу 0) = F2е2 + |
ехр (—ат) = А2+ АВ ехр (— ат) = |
|
2. (1— ftp) I |
|
|
= |
Л (Л + Д ехр (—ат)), |
(3.92) |
Этот метод известен как метод взаимно-корреляционного анализа.
Исключение члена с дельта-функцией Дирака в уравнении (3.92) является принципиальным отличием при использовании метода вза имной корреляции во времени с двумя детекторами. Как будет по казано ниже, это соответствует исключению постоянного фона в ча стотной области и позволяет проводить измерения детекторами с от
55