Файл: Уриг, Р. Статистические методы в физике ядерных реакторов.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. (1ftp) I

 

=

Л (Л + Д ехр (—ат)),

(3.92)

Этот метод известен как метод взаимно-корреляционного анализа.

Исключение члена с дельта-функцией Дирака в уравнении (3.92) является принципиальным отличием при использовании метода вза­ имной корреляции во времени с двумя детекторами. Как будет по­ казано ниже, это соответствует исключению постоянного фона в ча­ стотной области и позволяет проводить измерения детекторами с от­

55