Файл: Уриг, Р. Статистические методы в физике ядерных реакторов.pdf

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

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

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

Добавлен: 16.10.2024

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

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

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

ГЛАВА 12. СПЕЦИАЛЬНЫЕ СТАТИСТИЧЕСКИЕ МЕТОДЫ

И ИХ ПРИМЕНЕНИЕ

ВЯДЕРНЫХ РЕАКТОРАХ

§12.1. Введение

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

лезности. Применение многих из этих методов не ограничивает­ ся ядерными реакторами.

§12.2. Оптическое наблюдение корреляции

Вработе [1] доказана возможность применения оптического метода для детектирования корреляции между двумя сигналами. Присутствие некоррелированных и коррелированных событий в цеп­ ной реакции может выражаться в поведении выходных сигналов при визуальном наблюдении. В двухдетекторном корреляционном эксперименте; где детектируются независимые события, не ожи­ дается подобия двух выходных сигналов. Такой случай изображен на рис. 12.1, где показаны выходные сигналы для двух каналов в за­ висимости от времени. Центральная частота fc равна 1 гц, и ширина полосы А/ равна 0,95 гц. Сигналы не проявляют общих характерис­ тик. Иная картина имеет место в присутствии коррелированных со­ бытий, обусловленных ветвящимися процессами цепной реакции.

Если периоды наблюдения достаточно велики, коррелированные события цепной реакции будут проявляться в сигналах. Поэтому корреляция между сигналами должна выражаться в подобии ам­ плитуд и в пересечениях нулевой линии в эквивалентных точках на временной оси. Подобие выходных сигналов двух фильтров должно быть более отчетливым, когда увеличивается вклад коррелирован­ ных шумов. Это явление ясно видно на рис. 12.2. Здесь показаны выходные сигналы фильтров (\ = 1 гц и Af = 0,95 гц) для камер, расположенных рядом с активной зоной реактора.

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

381


Рис. 12.1 Некоррелированные выходные сигналы фильтров в двух­ детекторном взаимном корреляционном эксперимен­ те [1].

Рис. 12.2. Частично коррелированные выходные сигналы фильтров в двухдетекторном взаимном корреляционном экспери­ менте [1].

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

Рис. 12.3. Оптическое наблюдение корреляции [1]:

а — некоррелированные сиг­

налы;

о —

частично кор­

релированные

сигналы

(/=40

гц)\

в — частично

коррелированные

сигналы

(/=15

гц);

г — частично

коррелированные

сигналы

(/=1

гц);

д•—автокорре­

ляционные сигналы.

тремальный случай — автокорреляционный эксперимент

с двумя

идентичными сигналами, который дает прямую линию под углом 45° к осям х и у (см. рис. 12.3). Отношение доли коррелированных к некоррелированным шумам Q увеличивается на рис. 12,3, б, в, г из-за изменения центральной частоты от высоких к более низким

383


величинам. Таким образом, величина корреляции сигналов опреде­ ляется непосредственно из оптического наблюдения.

Впредположении, что распределения амплитуд коррелированных

инекоррелированных частей в сигналах х и у близки к гауссовому, кривые постоянной интенсивности света являются эллипсами. Для случая, в котором основные оси эллипса наклонены под углом 45° к осям х и у, получено следующее соотношение:

a2/62 = (L + l)/(L -l),

(12.1)

где а и Ь — большая и малая полуоси эллипса и

 

 

 

L =V (1+Qr2/3)(1+Qr2/3)-

(12.2)

 

Уравнение (12.1)

проверено в

работе [1]

при

измерениях на

 

реакторе

STARK,

и получено

совпадение

в пределах ± 2 0 % .

 

Если

коэффициент передачи в t-м канале А ;

выбран таким, что

 

А г/А 2

=

некоррелированные сигналы дают круг, и мы имеем

 

простое соотношение:

 

 

 

 

 

 

(a/6)2 = 1 + Q2/3 (cd) + Q!/3 (©).

(12.3)

§ 12.3. Анализ реакторных шумов с применением знаковой корреляции

Использование знака переменной величины в дискретные интер­ валы времени вместо ее численных значений в то же самое время эквивалентно оцифровке перемен­ ной с помощью одного значащего бита. Этот метод был предложен в работе [2] и усовершенствован

Драгтом [3].

 

Предположим, что х (0 = 0. За­

 

меним автокорреляционную

функ­

 

цию

знаковой

автокорреляцион­

 

ной функцией

фр

(т):

 

 

Фрхх СО = sign [х (t)] sign (t + т)].

 

 

 

 

 

(,12.4)

 

В работе [2] показано, что если

 

шумы

являются

гауссовыми, то

 

справедливо соотношение

 

Рис. 12.4. Распределение вероятно­

sin

 

 

Фжх СО

(12.5)

 

 

Ф х х ( 0 )

 

 

 

 

сти [3].

Этот метод проверялся для ана­

 

 

лиза

реакторных

шумов

путем

применения цифрового оборудования в измерениях шумов на реак­ торе LFR [3]. Эксперимент проводился следующим образом. Ток ионизационной камеры в реакторе с небольшой подкритичностью

384


квантовался по времени со скоростью 1000 выборок в секунду, оцифровывался (9 уровней) и вводился в память ЦВМ. Здесь под­ считывалась функция фка. (т) путем усреднения по времени. Затем все числа заменялись их знаками ( + 1 или — 1) относительно средней величины x (t) и опять рассчитывалась автокорреляцион­ ная функция в соответствии с уравнением (12.4). От второй функ­ ции находились синусоидальные значения в соответствии с уравне­ нием (12.5), затем обе функции сравнивались. Эксперимент был

2,Sr

з: Зс

 

 

 

сР00

 

си

2, 0 -

 

оо о~

 

 

<=>5 1<5'

0000°00°0оо0

°

°°Ооо°ООСоСОО

,оо

I *

1,0

TOono6 oQO °

 

°9оо

РоООгР00'

 

 

 

 

 

Рис. 12.5. Сравнение непосредственной и знаковой автокорреляционных функ­ ций шумов реактора LFR [3].

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

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

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

385


отношение cpxa; (х)/флзс (О). Разница между этой функцией и функ­

цией sin

фрхх С’1') показана в увеличенном в десять раз масштабе

(черные кружки) вместе со стандартными отклонениями. Между этими функциями не было обнаружено существенной разницы. В верхней части рисунка показано отношение стандартных откло­

нений функций sin и ер** (т)/фХЛ. (0). Удивительно, что

знаковая корреляционная функция имеет погрешность только приб­ лизительно в 1,5 раза больше, чем обычная корреляционная функ­ ция. Определение корреляционных функций для [sign х (/) ] вместо х (t) может рассматриваться как измерение х (I) с точностью до одного бита. Поскольку задержка, умножение и суммирование чи­ сел, представленных в виде одного бита, очень легко выполняются электронными схемами, метод знаковой корреляции представляет интересные возможности для прямого анализа реакторных шумов

врежиме on-line.

§12.4. Внутриреакторные измерения скорости теплоносителя и вибрации с помощью

внутризонных электродов и взаимной корреляции

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

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

386