ГЛАВА 12. СПЕЦИАЛЬНЫЕ СТАТИСТИЧЕСКИЕ МЕТОДЫ
И ИХ ПРИМЕНЕНИЕ
ВЯДЕРНЫХ РЕАКТОРАХ
§12.1. Введение
Вэтой главе представлены некоторые из специальных шумовых методов и их применения, которые логически не связаны с вопро сами, обсуждаемыми в других главах книги. Однако их расположе ние в конце книги отнюдь не имеет отношения к их важности и по
лезности. Применение многих из этих методов не ограничивает ся ядерными реакторами.
§12.2. Оптическое наблюдение корреляции
Вработе [1] доказана возможность применения оптического метода для детектирования корреляции между двумя сигналами. Присутствие некоррелированных и коррелированных событий в цеп ной реакции может выражаться в поведении выходных сигналов при визуальном наблюдении. В двухдетекторном корреляционном эксперименте; где детектируются независимые события, не ожи дается подобия двух выходных сигналов. Такой случай изображен на рис. 12.1, где показаны выходные сигналы для двух каналов в за висимости от времени. Центральная частота fc равна 1 гц, и ширина полосы А/ равна 0,95 гц. Сигналы не проявляют общих характерис тик. Иная картина имеет место в присутствии коррелированных со бытий, обусловленных ветвящимися процессами цепной реакции.
Если периоды наблюдения достаточно велики, коррелированные события цепной реакции будут проявляться в сигналах. Поэтому корреляция между сигналами должна выражаться в подобии ам плитуд и в пересечениях нулевой линии в эквивалентных точках на временной оси. Подобие выходных сигналов двух фильтров должно быть более отчетливым, когда увеличивается вклад коррелирован ных шумов. Это явление ясно видно на рис. 12.2. Здесь показаны выходные сигналы фильтров (\ = 1 гц и Af = 0,95 гц) для камер, расположенных рядом с активной зоной реактора.
Корреляция сигналов в двухдетекторном эксперименте наблю дается оптически более четко, если выходной сигнал одного фильтра подается на пластины осциллоскопа, отклоняющие луч в направле нии оси х, а выходной сигнал другого фильтра на пластины, откло-
Рис. 12.1 Некоррелированные выходные сигналы фильтров в двух детекторном взаимном корреляционном эксперимен те [1].
Рис. 12.2. Частично коррелированные выходные сигналы фильтров в двухдетекторном взаимном корреляционном экспери менте [1].
няющне луч в направлении оси у, и результирующий сигнал на эк ране фотографируется в течение времени, достаточного для получе ния геометрических фигур, показанных на рис. 12.3. Другой экс-
Рис. 12.3. Оптическое наблюдение корреляции [1]:
а — некоррелированные сиг
|
|
|
|
налы; |
о — |
частично кор |
релированные |
сигналы |
(/=40 |
гц)\ |
в — частично |
коррелированные |
сигналы |
(/=15 |
гц); |
г — частично |
коррелированные |
сигналы |
(/=1 |
гц); |
д•—автокорре |
ляционные сигналы. |
тремальный случай — автокорреляционный эксперимент |
с двумя |
идентичными сигналами, который дает прямую линию под углом 45° к осям х и у (см. рис. 12.3). Отношение доли коррелированных к некоррелированным шумам Q увеличивается на рис. 12,3, б, в, г из-за изменения центральной частоты от высоких к более низким
величинам. Таким образом, величина корреляции сигналов опреде ляется непосредственно из оптического наблюдения.
Впредположении, что распределения амплитуд коррелированных
инекоррелированных частей в сигналах х и у близки к гауссовому, кривые постоянной интенсивности света являются эллипсами. Для случая, в котором основные оси эллипса наклонены под углом 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]. Эксперимент проводился следующим образом. Ток ионизационной камеры в реакторе с небольшой подкритичностью
квантовался по времени со скоростью 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, где кружками показано
отношение cpxa; (х)/флзс (О). Разница между этой функцией и функ
цией sin |
фрхх С’1') показана в увеличенном в десять раз масштабе |
(черные кружки) вместе со стандартными отклонениями. Между этими функциями не было обнаружено существенной разницы. В верхней части рисунка показано отношение стандартных откло
нений функций sin и ер** (т)/фХЛ. (0). Удивительно, что
знаковая корреляционная функция имеет погрешность только приб лизительно в 1,5 раза больше, чем обычная корреляционная функ ция. Определение корреляционных функций для [sign х (/) ] вместо х (t) может рассматриваться как измерение х (I) с точностью до одного бита. Поскольку задержка, умножение и суммирование чи сел, представленных в виде одного бита, очень легко выполняются электронными схемами, метод знаковой корреляции представляет интересные возможности для прямого анализа реакторных шумов
врежиме on-line.
§12.4. Внутриреакторные измерения скорости теплоносителя и вибрации с помощью
внутризонных электродов и взаимной корреляции
В работе [4] сообщалось о методе измерения локальных скоро стей теплоносителя и распределения расходов в гидравлических мо делях реакторных активных зон с целью оптимизации теплопередачи и относительных характеристик потока теплоносителя. Проводящие электроды, смонтированные на стенках внутризонных каналов теп лоносителя, использовались для измерения изменений проводимости без нарушения потока теплоносителя. Изменения проводимости индуцировались путем инжекции малых количеств раствора соли в поток, входящий в модель активной зоны. Взаимные корреляцион ные измерения дают время запаздывания между сигналами двух детекторов.
Основой для измерений профиля потока является дальнейшее распространение рассмотренного выше метода, в котором раствор соли инжектируется против течения в секции, где проводится испы тание, и измеряются изменения сопротивлений в этой секции с по мощью проводящих элементов. Корреляционный метод включает в себя инжекцию раствора соли при медленной постоянной скоро сти. Гидравлическая турбулентность превращает концентрирован ный раствор соли в малые вихревые образования, которые распро страняются через систему и проявляются вблизи электродов как случайные изменения проводимости. Основная проблема заключает ся в измерении запаздывания переноса изменений проводимости между двумя электродами вдоль одного и того же канала теплоно сителя. Поскольку взаимная корреляционная функция между дву мя случайными когерентными сигналами имеет максимальную положительную величину при временном сдвиге, соответствующем