Приказ Федеральной службы по экологическому, технологическому и атомному надзору от 14 сентября 2011 г. N 535 "Об утверждении Положения о применении методов математической статистики для учета и контроля ядерных материалов"

В целях реализации полномочий, установленных Положением о Федеральной службе по экологическому, технологическому и атомному надзору, утвержденным постановлением Правительства Российской Федерации от 30 июля 2004 г. N 401, приказываю:

Утвердить прилагаемое Положение о применении методов математической статистики для учета и контроля ядерных материалов.

Руководитель

Н.Г. Кутьин

Положение
о применении методов математической статистики для учета и контроля ядерных материалов
(утв. приказом Федеральной службы по экологическому, технологическому и атомному надзору от 14 сентября 2011 г. N 535)

Список сокращений

В настоящем документе используются следующие сокращения:

	МБП	- межбалансовый период
	МВИ	- методика выполнения измерений
	ЗБМ	- зона баланса материалов
	ИР	- инвентаризационная разница
	КТИ	- ключевая точка измерений
	СКД	- средство контроля доступа
	СКО	- среднее квадратическое отклонение
	СФНК	- список фактически наличного количества
	ТВС	- тепловыделяющая сборка
	ТВЭЛ	- тепловыделяющий элемент
	УЕ	- учетная единица
	ФИ	- физическая инвентаризация
	ЯМ	- ядерный материал

I. Общие положения

1. Положение о применении методов математической статистики для учета и контроля ядерных материалов (далее - Положение) входит в число руководств по безопасности, носит рекомендательный характер и не является нормативным правовым актом.

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

3. Настоящее Положение рекомендуется использовать при разработке документов организации по статистической обработке результатов измерения количества и состава ЯМ в каждой ЗБМ, количества и состава ЯМ при передачах между ЗБМ, а также при подведении баланса ЯМ в ЗБМ.

4. Содержание настоящего Положения не охватывает всех существующих методов математической статистики. В нем представлены методы математической статистики, которые наиболее применимы в практике учета и контроля ЯМ, а также рекомендации по их применению для решения таких задач учета и контроля ЯМ как:

установление соответствия фактических параметров УЕ учетным данным;

определение объема подтверждающих измерений и анализ разницы между результатами учетных и подтверждающих измерений количественных параметров ЯМ, УЕ, продуктов;

определение объема случайной выборки при проверке пломб;

анализ расхождения данных по ЯМ организации-отправителя и организации-получателя;

проведение учетных и подтверждающих измерений фактически наличного количества ЯМ с учетом погрешностей измерений;

оценка величины неизмеренных потерь ЯМ и ее погрешности;

определение погрешности ИР, статистический анализ значимости ИР.

5. Рекомендации настоящего Положения касаются только ЯМ, учет и контроль которых осуществляется в форме УЕ.

6. Настоящее Положение содержит практические примеры возможного применения методов математической статистики.

II. Рекомендации по установлению соответствия фактических параметров УЕ существующим учетным данным

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

8. Статистическую обработку результатов однократных измерений для подтверждения количества ЯМ в выборке УЕ рекомендуется осуществлять на основе определенных характеристик.

Так, например, в случае определения массы ЯМ среднее значение массы ЯМ в выборке однородных (изготовленных по одним техническим условиям) УЕ рекомендуется вычислять по формуле:

, (1)

где - значение массы ЯМ в отдельных УЕ выборки;

n - число УЕ в выборке.

Для определения значений с помощью прямых или косвенных методов рекомендуется использовать аттестованные МВИ, определяющие порядок проведения измерений и гарантированную погрешность результата.

Выборочную дисперсию массы ЯМ в УЕ рекомендуется определять в соответствии с выражениями:

, (2)

или:

(2.1)

СКО и коэффициент вариации по выборке рекомендуется определять по формулам:

, (3)

и

. (4)

Расчет перечисленных характеристик рекомендуется выполнять с использованием интегрированных статистических пакетов. В случаях, когда при проведении измерений количества ЯМ в продуктах с использованием имеющихся МВИ обнаружено значимое расхождение с учетными данными рекомендуется убедиться в целесообразности дальнейшего использования перечисленных характеристик.

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

9. Для проверки гипотез рекомендуется подход, основанный на формировании статистического вывода о количестве ЯМ в УЕ с использованием двух типов гипотез: нулевой гипотезы, заключающейся в том, что действительное количество ЯМ равно заявленному и альтернативной гипотезы, заключающейся в том, что количество ЯМ в УЕ отличается от заявленного.

При проверке гипотез возможны ошибка, состоящая в отклонении нулевой гипотезы, когда она истинна (ошибка 1-го рода), и ошибка, состоящая в принятии нулевой гипотезы, когда истинной является альтернативная гипотеза (ошибка 2-го рода).

Вероятности возникновения ошибок 1-го рода и ошибок 2-го рода соответственно обозначаются, как и и .

10. Проверку согласия опытного распределения с нормальным законом рекомендуется проводить в соответствии с Правилами проверки согласия опытного распределения с теоретическим (рекомендации по стандартизации Р 50.1.037-2002). При этом в зависимости от числа проводимых измерений рекомендуется использовать различные критерии.

Например, при числе измерений, превышающем пятьдесят, рекомендуется использовать критерий .

При расчете рекомендуется табличный метод отображения, при котором значения, например, содержания каждого изотопа располагаются в порядке их возрастания и результаты упорядочения заносятся в первый столбец таблиц. Во второй столбец таблиц записываются значения функции распределения проверяемого теоретического распределения . В третий столбец таблиц записываются значения . В четвертый - значения (2j-1)/(2n), j=1, _, n. В пятый - произведения значений в столбцах (3) и (4) для j=1, _, n. В шестой - разность единицы и значений в столбце (4) для j=1, _, n. В седьмой - разность единицы и значений в столбце (2) для j=1, _, n. В восьмой столбец записываются значения . В девятый - произведения значений в столбцах (6) и (8) для j=1,_, n. В десятый - сумма значений в столбцах (5) и (9) для j=1, _, n. По таблицам рассчитываются значения . Здесь суммирование по j=1, _, n. Рассчитанные значения сравниваются с критическим для уровня значимости 0,05 значением, приведенным в Правилах проверки согласия опытного распределения с теоретическим (рекомендации по стандартизации Р 50.1.037-2002).

11. Проверку стохастической независимости результатов параллельных определений рекомендуется проводить на основе критерия отношений квадратов последовательных разностей. Для этого для каждых n одновременных измерений рассчитывается значение

, (5)

где , , .

В случае , где , а - табулированное значение, гипотеза о стохастической независимости результатов параллельных определений не отвергается.

12. Проверку незначимости расхождения результатов параллельных определений

() для каждой серии из n измерений при доверительной вероятности 0,95 рекомендуется оценивать с помощью неравенства:

, (6)

где верхняя граница относительного СКО результатов измерений , доверительные границы (интервал) случайной погрешности результата измерений , среднее арифметическое , СКО . - значение коэффициента Стьюдента с (n-1) степенями свободы при доверительной вероятности 0,95. При определении значения коэффициента используются существующие в этой области стандарты и рекомендации. В случае выполнения неравенства оснований для исключения из анализа i-го значения как аномального нет.

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

. (7)

Здесь , а - наибольшее из . Расчет производится по формуле: , N - число параллельных определений содержания каждого изотопа (элемента) в каждом образце, а m - количество образцов каждого содержания.

Для проверки незначимости различия дисперсий, рассчитанные значения сравниваются с критическим для уровня значимости 0,05 значением , которое определяется в соответствии с существующими в этой области стандартами и рекомендациями.

14. При исследовании однородности образцов содержания элемента (изотопа) для каждого содержания по каждому элементу (изотопу) рассчитывается сумма квадратов отклонений результатов определений внутри проб (под пробами здесь понимаются образцы одного содержания) и между средними арифметическими по пробам . Здесь - среднее арифметическое значение J параллельных определений () каждого элемента (изотопа) в n - ом образце исследуемого содержания, а - среднее арифметическое значение всех NJ определений . Затем вычисляются выборочные средние квадраты отклонений результатов внутри проб и между пробами . Для неделимых образцов СКО, характеризующие неоднородность материала по содержанию каждого элемента (изотопа) в образцах каждого содержания элемента (изотопа), рассчитываются по формуле: .

Пример сравнения двух алгоритмов обработки спектров приведен в Приложении N 1 к настоящему Положению.

III. Определение объема подтверждающих измерений и УЕ для измерений. Анализ разницы между результатами учетных и подтверждающих измерений

15. При наличии возможностей рекомендуется проводить подтверждающие измерения всего ЯМ, находящегося в ЗБМ. В противном случае рекомендуется проводить измерения УЕ, случайным образом выбираемых из общего их количества в ЗБМ. При этом рекомендуется применять дифференцированный подход, учитывающий влияние трех факторов: категории ЯМ, объема применения к нему СКД и массы ЯМ в УЕ.

16. При осуществлении дифференцированного подхода к выбору УЕ для проведения подтверждающих измерений рекомендуется распределить их по стратам (совокупностям отдельных УЕ с одинаковыми или близкими физическими характеристиками и химическим составом ЯМ), и составить перечень УЕ, входящих в каждую страту. При этом объем выборки задается двумя величинами: G - пороговым количеством ЯМ в единицах массы и Р - вероятностью обнаружения недостачи/излишка порогового количества ЯМ.

17. Пороговое количество G рекомендуется определять в зависимости от категории ЯМ.

18. Вероятность обнаружения недостачи/излишка P рекомендуется определять в соответствии с таблицей 1.

Таблица 1

Вероятность обнаружения недостачи/излишка порогового количества ЯМ для расчета объема выборки подтверждающих измерений

Объем применения СКД к ЯМ	Вероятность обнаружения P, %, не менее
	Для категорий 1, 2, 3	Для категории 4
Только пломбы	50	30
Только система наблюдения	50	30
Одновременно два вида разнообразных СКД	25	9
Одновременно m видов разнообразных СКД

19. Если ЯМ представлены в форме УЕ, целостность которых при нахождении в ЗБМ в МБП обеспечена их конструкцией (например, неразборные тепловыделяющие сборки или диски, блочки, изделия с ЯМ, покрытые оболочками, в результате чего ЯМ не могут быть извлечены из них без нарушения целостности оболочек) и может быть подтверждена (например, визуальным осмотром, сравнением образов гамма- или ультразвукового сканирования швов сварки), то для этих ЯМ при расчете объема выборки подтверждающих измерений рекомендуется использовать значения вероятностей обнаружения, аналогичные применяемым к ЯМ с пломбами. При этом рекомендуется принимать во внимание наличие документального подтверждения целостности неразборных УЕ и отсутствия несанкционированного доступа, а также фактической проверки СКД при проведении ФИ.

20. Расчет объема выборки (n) рекомендуется осуществлять в соответствии с Приложением N 2 к настоящему Положению, по формуле:

, (8)

N - количество УЕ в страте;

x - средняя масса ЯМ в одной УЕ;

- в формуле означает округление до ближайшего большего целого числа.

21. Номера подлежащих измерению УЕ рекомендуется выбирать случайным образом из перечня проверяемых УЕ каждой страты.

22. Случайный выбор рекомендуется проводить с помощью генератора случайных чисел с предварительным присвоением каждой УЕ проверяемой страты порядкового номера от 1 до N.

23. После измерения количественных параметров ЯМ, УЕ, продуктов рекомендуется проводить анализ разницы между результатами учетных и подтверждающих измерений.

Для этого рекомендуется рассчитать значение контрольного допуска по формуле:

. (9)

Здесь - СКО подтверждающих измерений, а - СКО учетных измерений.

24. При проведении измерений УЕ нескольких страт в случае необходимости объединения результатов измерений рекомендуется для нормально распределенных случайных величин выполнять проверку гипотез о равенстве средних значений и дисперсий. Это может потребоваться в задачах объединения нескольких страт, либо в задачах разбиения одной страты на несколько для получения представительных результатов. Пример практического решения такой задачи приведен в Приложении N 3 настоящего Положения.

Две выборочные дисперсии рекомендуется сравнивать с помощью критерия Фишера F. Для этого вычисляют отношение большей выборочной дисперсии к меньшей.

Если отношение выборочных дисперсий при больше или равно с числом степеней свободы и , где - объемы выборок, а - уровень значимости, то рекомендуется принять гипотезу о неравенстве двух генеральных дисперсий и, следовательно, нельзя объединять выборки. В случае соблюдения условия рекомендуется принять гипотезу о равенстве генеральных дисперсий, и переходят к проверке равенства выборочных средних значений.

Средние значения нормально распределенных случайных величин рекомендуется сравнивать с помощью t - критерия Стьюдента, значения которого приведены в анализе характеристик выборочного распределения (Приложение N 3 настоящего Положения). Для этого, в случае доказанного равенства выборочных дисперсий, вычисляют сводную дисперсию в соответствии с выражением:

(10)

и значение t согласно выражению:

(11)

Если , то рекомендуется принять гипотезу о равенстве средних. Здесь - уровень значимости, а k - число степеней свободы, определяемое выражением . В случае .

Рекомендуемый порядок применения контрольного допуска, а также примеры сравнения двух выборочных дисперсий и средних приведены в Приложении N 4 к настоящему Положению.

IV. Определение объема случайной выборки при проверке пломб

25. Для определения объема выборки проверяемых пломб рекомендуется рассчитать значение ; где N - количество установленных пломб в ЗБМ. Здесь также означает округление до ближайшего большего целого числа. Значения объема выборки для различных значений N приведены в Приложении N 2 к настоящему Положению.

26. Номера подлежащих проверке пломб рекомендуется выбирать случайным образом.

V. Анализ расхождения данных организации-отправителя и организации-получателя

27. Для анализа расхождения данных по ЯМ организации-отправителя и организации-получателя рекомендуется определить допустимые границы такого расхождения:

, (12)

где - СКО результатов измерений массы с помощью весов отправителя и получателя. В общем случае:

, (13)

где и - значения случайных составляющих погрешностей весов отправителя и получателя соответственно, а

и - значения систематических составляющих погрешностей весов отправителя и получателя соответственно.

На практике в технической документации на весы, как правило, указана полная погрешность, приведенная к доверительному интервалу 0,99, уже включающая в себя как систематическую так и случайную составляющие погрешности (например: в абсолютных единицах г или в относительных %). В этом случае расчет ведется по формуле:

, (14)

или

. (15)

В случае отсутствия данных по весам отправителя соответствующие погрешности принимаются равными нулю и соответственно допустимые границы расхождения данных по ЯМ организации-отправителя и организации-получателя будут более "жесткими". Примеры расчета , а также алгоритм оценки значимости систематических расхождений приведены в Приложении N 5 к настоящему Положению.

28. Если расхождение данных по ЯМ организации-отправителя и организации-получателя , (здесь - масса ЯМ по данным организации-отправителя, а - масса ЯМ по данным организации-получателя), то ЯМ рекомендуется ставить на учет в организации-получателе по данным организации-отправителя.

VI. Проведение учетных и подтверждающих измерений фактически наличного количества ЯМ с учетом погрешностей измерений

29. В учете и контроле ЯМ наиболее распространена задача сопоставления результатов подтверждающих измерений ЯМ с учетными данными. Эта задача описана в разделе III настоящего Положения. Однако иногда возникает необходимость проанализировать данные по партии ЯМ (раздел V настоящего Положения), стратам ЯМ. Фактически наличное количество ЯМ в страте из N УЕ, имеющих массу , рекомендуется определять следующим образом:

, (16)

где - среднее учетное значение массы УЕ.

30. При оценке фактически наличного количества ЯМ по выборке из ограниченной совокупности, где , , - истинное значение массы ЯМ для i-ой УЕ, а - неопределенность этого значения, рекомендуется рассматривать, по крайней мере, три компоненты :

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

- вторая: - функция оценок калибровочной кривой, причем коэффициенты последней обычно зависят от измеряемой массы. При измерении двух разных УЕ каждой из них является причиной ковариации измерений, а ковариация в свою очередь зависит от масс этих УЕ;

- третья: - случайная погрешность, связанная со статистикой счета. Она принимает независимые значения при повторных измерениях одной УЕ, а также для разных УЕ. В общем случае является функцией измеряемой массы.

Выборки из страты обычно являются бесповторными. Случайные величины одинаково распределены, но в отличие от случая выборки из бесконечной генеральной совокупности не являются независимыми. Эта зависимость определяется ковариацией, которая в рассматриваемом случае равна:

. Здесь . (17)

Связь характеристик выборки с характеристиками страты без учета погрешности измерения рекомендуется описывать следующим образом:

, ,

.

Здесь E(x) и Var(x) - математическое ожидание и дисперсия величины x соответственно.

31. Оценка фактически наличного количества ЯМ стремится к истинному значению при уменьшении суммы систематических отклонений, которые вносило бы измерение каждой УЕ страты:

- дисперсия оценки. Здесь , где <m> - среднее значение , а - оценки элементов матрицы дисперсий погрешности измерений. Такие оценки условных значений дисперсии и ковариации зависят от массы и служат оценками математических ожиданий для выборки этих значений дисперсии и ковариации. Последнее позволяет работать с выборками из конечной совокупности с учетом погрешностей измерений, не ограничивая себя конкретной моделью погрешности. Условные значения дисперсии и ковариации обычно содержатся в описании и программном обеспечении большинства современных измерительных систем.

Оценка фактически наличного количества ЯМ в УЕ страты стремится к истинному значению при уменьшении суммы систематических отклонений, которые вносило бы измерение каждой УЕ страты. Для минимизации суммы систематических отклонений рекомендуется проведение более тщательной калибровки, учитывающей зависимость систематических отклонений от измеряемой массы, а также других влияющих на измерение факторов.

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

VII. Определение погрешности ИР. Статистический анализ значимости ИР

32. ИР рекомендуется вычислять по элементу и по изотопу для каждого ЯМ:

ИР = ФК - ДК = ФК - (НК + УВ - УМ), (18)

где ФК - фактически наличное количество ЯМ в ЗБМ, определенное в результате данной ФИ;

ДК - документально зарегистрированное количество ЯМ в ЗБМ на начало ФИ (конец данного МБП);

УВ - документально зарегистрированное увеличение количества ЯМ в ЗБМ за данный МБП в результате всех поступлений, наработок и т.д.;

УМ - документально зарегистрированное уменьшение количества ЯМ в ЗБМ за данный МБП в результате всех отправок из ЗБМ, ядерных превращений, потерь и т.д.;

НК - документально зарегистрированное количество ЯМ в ЗБМ на начало данного МБП.

Уравнение ИР можно представить как сумму ЯМ (элемента или изотопа).

, (19)

где: N1 - количество УЕ в НК

N2 - количество УЕ в УВ

N3 - количество УЕ в УМ

N4 - количество УЕ в ФК

- масса элемента при определении ИР по элементу (массовая доля С):

, (20)

где: - масса нетто;

- масса изотопа при определении ИР по изотопу (массовая доля ):

(21)

33. Для оценки погрешности ИР рекомендуется использовать метод переноса погрешности измерений.

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

(22)

С учетом (22) выражение для определения погрешности ИР может быть представлено в следующем виде:

, (23)

где i - индекс строки матрицы дисперсий и ковариаций;

j - индекс столбца матрицы дисперсий и ковариаций;

- погрешность результатов измерения количества ЯМ в компоненте уравнения баланса с индексом i;

- погрешность результатов измерения количества ЯМ в компоненте уравнения баланса с индексом j;

- коэффициент корреляции между компонентами уравнения баланса с индексами i и j.

Структура погрешности результатов измерения ЯМ рассмотрена в разделе VI настоящего Положения. В этом же разделе приведены аналитические выражения для расчета дисперсий и ковариаций. В таких расчетах рекомендуется использовать данные, содержащиеся в документации на средства измерения, в МВИ.

Причиной корреляции между компонентами уравнения баланса может служить использование при проведении учетных измерений ЯМ одних и тех же средств измерения.

35. Для вычисления СКО величин каждого из компонентов уравнения баланса все УЕ, которые подвергаются измерениям, рекомендуется разбить на несколько (K) независимых друг от друга (некоррелированных) страт. Дисперсия результатов измерений ЯМ в УЕ k-ой страты в общем случае с учетом возможных корреляций определения параметров ЯМ в УЕ, составляющих данную страту, может быть вычислена по формуле:

, (24)

где - дисперсия определения массы ЯМ l-ой УЕ, входящей в k-ую страту, содержащую УЕ, с учетом случайной и систематической составляющих погрешности определения величины массы ЯМ в l-ой УЕ - ;

индексы l и m соответствуют порядковым номерам УЕ, входящих в k-ую страту;

- коэффициент корреляции между значениями масс ЯМ в l-ой УЕ - и в m-ой УЕ - , входящих в k-ую страту.

Величина коэффициентов корреляции может изменяться в пределах от -1 до +1.

Для независимых величин коэффициент корреляции между ними равен нулю.

Определение значения величины коэффициента корреляции при анализе результатов измерений является достаточно трудоемкой задачей. Поэтому в большинстве случаев для упрощения обработки результатов измерений и вычисления величины возможными корреляциями между компонентами уравнения баланса, а также между результатами измерений параметров ЯМ в стратах УЕ рекомендуется пренебречь. При этом значения , определяются в соответствии с выражениями:

, (25)

(26)

Это в конечном итоге приведет к некоторому уменьшению найденного значения и, следовательно, повышению риска ошибки первого рода (фиксации аномалии в учете и контроле при фактическом ее отсутствии). Однако, следует полагать, что при правильной организации учета и контроля ЯМ в ЗБМ данное явление будет происходить достаточно редко, что обусловливает приемлемость допущения об отсутствии корреляций между соответствующими величинами при определении , в практических задачах.

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

36. Для определения дисперсий масс ЯМ в УЕ в организации рекомендуется разработать методики определения массы ЯМ и ее погрешности с учетом оснащения ЗБМ техническими средствами для проведения учетных измерений, МВИ, особенностей технологических процессов.

37. При расчете погрешности определения массы рекомендуется учитывать такие источники погрешности как: статистическая погрешность отбора УЕ, погрешность балк - измерения, погрешность пробоотбора материала, погрешность химического анализа материала и другие в зависимости от особенностей технологических процессов.

38. Для тех страт УЕ с ЯМ в ЗБМ, которые не подвергались каким-либо превращениям в течение МБП, при вычислении ИР и рекомендуется использовать учетные данные, полученные ранее для входящих в их состав УЕ. При этом такие данные какого-либо влияния на величину не оказывают.

39. Для вычисления и анализа результатов на наличие возможных аномалий в учете и контроле ЯМ рекомендуется порядок, изложенный в Стандарте корпорации ОАО "ТВЭЛ". "Система учета и контроля ядерных материалов", в соответствии с которым производится:

установление СКО систематических () и случайных () составляющих погрешностей всех методов измерений ЯМ в ЗБМ;

подготовка данных регистрации движения ЯМ для вычисления ИР и . выделение групп (страт) ЯМ;

вычисление ИР по элементу и по изотопу для каждого ЯМ в ЗБМ;

исключение из расчета "парных учетных записей", которые не должны влиять на ;

вычисление (дисперсий ИР);

проверка по критериям обнаружения аномалий в учете и контроле ЯМ.

40. При определении характеристик погрешностей измерений (СКО систематических и случайных составляющих) рекомендуется:

установить на основании МВИ ЯМ случайные и систематические составляющие характеристики погрешностей измерения ЯМ в ЗБМ;

для вычисления вычислить СКО для каждой случайной и систематической составляющей погрешности всех методов измерения ЯМ в ЗБМ. Если погрешность МВИ регламентируется в абсолютной форме - , то рекомендуется следовать аддитивной модели погрешности. Если погрешность МВИ регламентируется в относительной форме - , то рекомендуется следовать мультипликативной модели. Описание аддитивной и мультипликативной моделей погрешности измерения, а также связь с абсолютной и относительной погрешностями измерения приведены в Приложении N 6 к настоящему Положению.

Переход от интервальных характеристик погрешности к СКО () при P=0,95 рекомендуется осуществлять следующим образом:

- СКО случайной составляющей для абсолютной погрешности (аддитивная модель);

- СКО случайной составляющей для относительной погрешности (мультипликативная модель);

где: , - интервальные характеристики случайной составляющей погрешности измерения при P=0,95;

- СКО систематической составляющей для абсолютной погрешности (аддитивная модель);

- СКО систематической составляющей для относительной погрешности (мультипликативная модель);

где: , - интервальные характеристики систематической составляющей погрешности измерения при P=0,95.

Для процедур взвешивания в расчетах рекомендуется использовать паспортизованные значения погрешностей при Р=1.

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

Пример рекомендуемой формы представления итоговых результатов определения СКО составляющих погрешностей измерения представлен в Приложении N 7 к настоящему Положению.

41. При подготовке данных регистрации движения ЯМ к вычислению ИР и , при выделении групп (страт) ЯМ рекомендуется:

все поступления и другие увеличения ЯМ в ЗБМ, отправления и другие уменьшения ЯМ (например, образование отходов) из ЗБМ, результаты предыдущей и текущей ФИ зафиксировать в отчетных документах в соответствии с документами предприятия по учету и контролю ЯМ. Пример регистрации движения ЯМ для ЗБМ производства ТВЭЛ представлен в таблице П8.1 Приложения N 8 к настоящему Положению;

выделить группы (страты) ЯМ. Каждая группа (страта) характеризуется набором одних и тех же методов измерений (например: взвешивание на одних и тех же весах, определение массовой доли элемента по одной и той же МВИ, определение массовой доли изотопа по одной и той же МВИ) и соответствующими характеристиками погрешности (пример в таблице П8.1 Приложения N 8 к настоящему Положению).

В каждой группе (страте) содержится информация о движении ЯМ в ЗБМ за МБП. Пример выделения групп (страт) для ЗБМ производства ТВЭЛ представлен в таблицах П8.2, П8.3, П8.4 Приложения N 8 к настоящему Положению.

42. При вычислении ИР рекомендуется пользоваться формулой (20). Вычисления проводятся для каждого ЯМ (элемента и изотопа).

43. Для исключения из расчета "парных учетных записей" рекомендуется:

при анализе таблицы регистрации движения ЯМ (пример - таблица П8.1 Приложения N 8 к настоящему Положению выделить такие учетные записи, которым соответствует одна и та же УЕ (с одним и тем же идентификационным номером), которая осталась без изменения и работы с ней за МБП не выполнялись, а достоверность учетных данных для этой УЕ с момента их определения до момента использования подтверждена надлежащим состоянием примененных СКД;

если для каждой такой УЕ существуют две учетные записи, соответствующие разным членам уравнения ИР с разными знаками, то такие "парные учетные записи" необходимо исключить из расчета , чтобы не допустить необоснованного увеличения . Например, возможны следующие типы "парных учетных записей":

Члены уравнения инвентаризационной разницы УВ и ФК.

Поступившие в ЗБМ УЕ (УВ), которые зарегистрированы в соответствии с данными отправителя, затем при проведении ФИ эти же УЕ были занесены в инвентаризационную ведомость (ФК).

Члены уравнения инвентаризационной разницы УВ и УМ.

Поступившие в ЗБМ УЕ (УВ), которые затем были отправлены из ЗБМ (УМ).

Члены уравнения инвентаризационной разницы НК и УМ.

Часть УЕ из предыдущей ФИ (НК) были отправлены из ЗБМ (УМ).

Члены уравнения инвентаризационной разницы НК и ФК.

Часть УЕ из предыдущей ФИ (НК) вошли в список УЕ текущей инвентаризации (ФК).

Пример исключения парных записей приведен в Приложении N 9 к настоящему Положению.

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

Подход с использованием усреднения по потоку заключается в предположении, что в среднем все учетные массы в партии ЯМ являются примерно одинаковыми и масса одной УЕ равна средней массе УЕ по партии. При использовании этого подхода получаемое значение , как правило, определяет целесообразность применения более "жестких" критериев обнаружения аномалий, чем при вычислении без усреднения по потоку.

Подход к вычислению без усреднения по потоку состоит в том, что используются индивидуальные учетные массы ЯМ по каждой УЕ. Для реализации этого подхода рекомендуется разработать специальное программное обеспечение.

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

При принятии решения об использовании того или иного метода рекомендуется сделать оценку величины, на которую уменьшается вычисляемая с усреднением по потоку (Приложение N 10 к настоящему Положению).

45. При вычислении (c усреднением по потоку) рекомендуются аналитические выражения для относительных погрешностей измерения как наиболее часто применяемых в практике при вычислении дисперсии ИР (мультипликативная модель).

При этом вклад систематической составляющей погрешности по одному из методов измерения (взвешивание, измерение объема, других МВИ) в дисперсию ИР для j-той группы (страты) рекомендуется определять с помощью выражения:

, (27)

где - СКО систематической составляющей погрешности одного из методов измерения,

(28)

для из НК, УВ - членов уравнения ИР,

для из ФК, УМ - членов уравнения ИР. Этот подход справедлив, например, при взвешивании всего ЯМ, составляющего члены уравнения ИР, на одних и тех же весах без перекалибровки.

При проведении взвешиваний и определении массовой доли элемента (изотопа) рекомендуется следующий подход.

При вычислении ИР по элементу:

, (29)

где - масса брутто;

- масса тары;

- массовая доля элемента.

При вычислении ИР по изотопу:

, (30)

где: - массовая доля изотопа.

При вычислении вклада систематической составляющей погрешности в определяется общая сумма масс ЯМ в группе (страте) с учетом их знака в уравнении ИР. В связи с тем, что результаты суммы масс как для метода с усреднением по потоку, так и для метода без усреднения, являются одинаковыми, то эти два подхода дают эквивалентные результаты. Порядок вычисления суммарной массы j-той группы (страты) для определения вклада систематических и случайных составляющих погрешностей в приведен ниже.

Вклад случайной составляющей погрешности по одному из методов измерения (взвешивание, измерение объема, МВИ разрушающих и неразрушающих методов анализа) в для j-той группы (страты) рекомендуется вычислять с использованием выражения:

, (31)

где - суммарная масса элемента (изотопа) j-той группы (страты);

- число измерений в j-той группе (страте);

- случайная составляющая погрешности одного из методов измерения j-той группы (страты).

Суммарная масса элемента (изотопа) j-той группы (страты):

, (32)

где R - число партий в j-той группе (страте);

- суммарная масса элемента (изотопа) в i-той партии.

Общее число взвешиваний (измерений объема и т.п.) в j-той группе (страте):

(33)

где R - число партий в j-той группе (страте);

- число контейнеров (измерений объема) в i-той партии.

Общее число анализов по МВИ одним из разрушающих (неразрушающих) методов анализа в j-той группе (страте):

(34)

где R - число партий ЯМ, в которых определяется массовая доля элемента (изотопа) по МВИ для j-той группы (страты);

- число проб на партию.

Аналитические выражения для вычисления без усреднения по потоку рассмотрены в Приложении N 11 к настоящему Положению.

При вычислении суммарной массы в j-той группе (страте) для определения вклада систематических и случайных составляющих погрешностей в (с усреднением по потоку) к его началу рекомендуется провести процедуру исключения "парных учетных записей".

Для вычисления вклада в систематических составляющих погрешностей рекомендуется суммировать массы элемента (изотопа) с учетом знака уравнения ИР для случая, когда при взвешивании всего ЯМ, составляющего члены уравнения ИР используются одни и те же средства измерения без перекалибровки. В примере исключения парных записей, представленном в Приложении N 9 к настоящему Положению суммируются суммы масс элемента (изотопа) и записываются в строку (сист) по элементу и изотопу:

(35)

Для вычисления вклада в случайных составляющих погрешностей рекомендуется суммировать массы элемента (изотопа) без учета знака уравнения ИР для случая, когда при взвешивании всего ЯМ, составляющего члены уравнения ИР, используются одни и те же средства измерения без перекалибровки. В примере исключения парных записей, представленном в Приложении N 9 к настоящему Положению суммируются суммы масс элемента (изотопа) и записываются в строку по элементу и изотопу:

(36)

Пример вычисления суммарной массы материала приведен в Приложении N 12 к настоящему Положению.

При вычислении значения рекомендуется использовать следующие данные:

- характеристики погрешностей ( и );

- количество измерений (количество контейнеров при однократном взвешивании);

- количество проб, которые были взяты из партии ЯМ для определения массовой доли элемента и изотопа.

В Приложении N 13 к настоящему Положению приведен пример вычисления вклада систематических и случайных составляющих погрешностей j-ой группы (страты) в .

Для вычисления вклада в систематических составляющих погрешностей для тех методов измерения, которые выполнялись в j-ой группе (страте) рекомендуется использовать выражение (27). Определяется суммарный вклад систематических составляющих погрешностей j-ой группы (страты) в .

Например:

- - вклад систематической составляющей погрешности взвешивания j-ой группы (страты) в ;

- - вклад систематической составляющей погрешности разрушающего анализа (по элементу и по изотопу) j-ой группы (страты) в .

Определяется общий вклад в систематических составляющих погрешностей как сумма отдельных вкладов (взвешивание и анализ):

(37)

Для вычисления вклада в случайных составляющих погрешностей для тех методов измерения, которые выполнялись в j-ой группе (страте) рекомендуется использовать выражение (28). Определяется суммарный вклад случайных составляющих погрешностей j-ой группы (страты) в .

Например:

- - вклад случайной составляющей погрешности взвешивания j-ой группы (страты) в ;

- - вклад случайной составляющей погрешности разрушающего анализа (по элементу и по изотопу) j-ой группы (страты) в .

Определяется общий вклад в случайных составляющих погрешностей как сумма отдельных вкладов (взвешивание и анализ):

(38)

Суммируются вклады в систематических составляющих погрешностей всех групп (страт):

, (39)

где: - вклад систематических составляющих погрешностей в j-ой группы (страты);

K - число групп (страт).

Суммируются вклады в случайных составляющих погрешностей всех групп (страт):

, (40)

где: - вклад случайных составляющих погрешностей в j-ой группы (страты).

Рекомендуется провести сравнительный анализ вносимой доли каждой из составляющих погрешности в .

Суммарная складывается из систематических и случайных составляющих погрешностей измерений всех групп (страт):

(41)

СКО инвентаризационной разницы:

(42)

46. В качестве критерия обнаружения аномалий в учете и контроле ЯМ за МБП рекомендуется использовать выполнение хотя бы одного из неравенств следующей системы:

(43)

Числовые значения в третьем и четвертом неравенствах системы (43) определяют критические доли от общего количества данного ЯМ, преобразованного и подвергшегося учетным измерениям в данный МБП или в процессе ФИ.

Проверку выполнения второго неравенства системы рекомендуется проводить только в отношении ЯМ, для которых пороговые значения G установлены.

Проверку выполнения третьего и четвертого неравенств рекомендуется проводить отдельно для каждого ЯМ в ЗБМ по элементу и по изотопу.

Приложение N 1

к Положению по применению методов

математической статистики для учета и контроля

ядерных материалов,

утв. приказом Федеральной службы по экологическому,

технологическому и атомному надзору

от 14 сентября 2011 г. N 535

Пример сравнения двух алгоритмов обработки спектров

При разработке стандартных образцов содержания плутония для инспекционных измерений эти образцы измерялись с помощью инспекционной гамма - спектрометрической станции U-Pu InSpector с низкоэнергетическим германиевым детектором (LEGe) активной площадью 5 . Для получения данных об изотопном составе использовалось программное обеспечение MGA. Результаты измерения были обработаны с целью их сопоставления с данными, полученными с помощью разрушающего анализа проб материала, а также для сопоставления двух предусмотренных программным обеспечением станции алгоритмов: с декларацией содержания Pu-242 до измерения образца и без таковой.

Измерениям подверглись 11 комплектов стандартных образцов, в каждый из которых входило по четыре образца разного изотопного состава. Данные о составе были получены с помощью разрушающих методов анализа проб материала до изготовления образцов и представлены до проведения инспекционных измерений. Так, например, процентное содержание Pu-239 в образцах состава N 1 было равно , для состава N 2 - , для состава N 3 - , а для состава N 4 - . Для каждого из 44 образцов было выполнено 8 параллельных определений изотопного состава с обработкой спектра по алгоритму, предусматривающему декларацию содержания Pu-242 до измерения образца. Для обработки каждого девятого спектра использовался алгоритм, не требующий декларировать содержание Pu-242. Разрешение на 122 кэВ при измерениях составляло (590-600) эВ (ПШПВ). Время набора каждого спектра составляло 1000 с. Путем изменения расстояния между детектором и образцом, а также подбором толщины фильтра для снижения интенсивности низкоэнергетических фотонов "мертвое" время для образцов разного содержания регулировалось в диапазоне (18-22)%. При этом общее за время набора каждого спектра число отсчетов составляло , а в энергетических областях при (59,54; 129,29; 208,00) кэВ - не менее при сопоставимых значениях. Смещение центроид, соответствующих перечисленным энергиям, находилось при измерениях в пределах пяти каналов. Необходимости коррекции смещения усиления в период измерений не возникало. Коэффициент, соответствующий усилению в уравнении калибровки, составлял ~ 0,075 кэВ.

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

согласия опытного распределения с теоретическим;

стохастической независимости результатов параллельных определений;

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

незначимости различия дисперсий по критерию Кохрена;

однородности.

Для проверки согласия опытного распределения с нормальным законом значение критерия для малых объемов выборки сравнивалось с критическим для уровня значимости 0,05. В ходе такой проверки несоответствия данных параллельных определений нормальному закону для уровня значимости 0,05 установлено не было.

Стохастическая независимость результатов параллельных определений проверялась с помощью критерия отношений квадратов последовательных разностей. В ходе такой проверки случаев невыполнения неравенства для уровня значимости 0,05 установлено не было.

Проверка незначимости расхождения результатов () для каждой серии из n=8 параллельных определений для доверительной вероятности 0,95 аномальных значений в экспериментальных данных не выявила.

Проверка дисперсий для каждого изотопа образцов каждого содержания значимых различий дисперсий не выявила.

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

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

При проверке согласия опытных распределений объединенных данных с нормальным законом (число измерений превышает пятьдесят) по критерию доказательств несоответствия для уровня значимости 0,05 опытного распределения нормальному закону установлено не было.

На следующем этапе по каждому изотопу формировались таблицы, в каждый из четырех столбцов которых записывались результаты измерения образцов одного содержания плутония. Для анализа таблиц использовался обычный интегрированный статистический пакет, позволяющий, не углубляясь в теорию, выполнить стандартные статистические расчеты. Для таблицы по анализируемому изотопу рассчитывались арифметические средние и дисперсии измерений четырех стандартных образцов содержания плутония, дисперсии сравнивались по критерию Фишера, определялся вид модели погрешностей измерения. Затем с использованием однофакторного дисперсионного анализа определялись случайная и систематическая составляющие погрешностей измерения по следующему алгоритму. Если в ввести следующие обозначения: , , и считать, что - распределено по нормальному закону (индекс "R" относится к случайной составляющей погрешностей измерения, а "S" - к систематической), то получим:

,

, .

Здесь k = 4 - число стандартных образцов в комплекте; n = 88 - количество измерений образцов одного содержания плутония; - для i-го образца j-го содержания разность между измеренным с помощью U-Pu InSpector значением и значением, полученным в аналитической лаборатории; ; .

Аналогично обрабатывались результаты измерений без декларации содержания Pu-242. Только в этом случае анализировались таблицы размером .

Результаты обработки приведены в таблицах ниже (в каждой ячейке таблицы последовательно записаны данные для образцов состава N 1, N 2, N 3 и N 4).

Таблица П1.1

Результаты анализа измерений с декларацией содержания Pu-242

Изотоп Характеристика	Pu-238	Pu-239	Pu-240	Pu-241	Am-241
Среднее квадратическое отклонение, характеризующее неоднородность материала стандартного образца - по каждому изотопу для образцов одного содержания плутония.
Случайная составляющая погрешностей измерения по каждому изотопу для образцов одного содержания плутония.
Систематическая составляющая погрешностей измерения по каждому изотопу для образцов одного содержания плутония.
Доверительные границы (интервал) случайной погрешности результата измерений по каждому изотопу для образцов одного содержания плутония.
Верхняя граница относительного среднего квадратического отклонения результатов определений по каждому изотопу для образцов одного содержания плутония.
Доверительные границы погрешности результата анализа по каждому изотопу для образцов одного содержания плутония.
Доверительные границы погрешности с учетом неоднородности образцов по каждому изотопу для образцов одного содержания плутония.

Таблица П1.2

Результаты анализа измерений без декларации содержания Pu-242

Изотоп Характеристика	Pu-238	Pu-239	Pu-240	Pu-241	Pu-242	Am-241
Среднее квадратическое отклонение, характеризующее неоднородность материала стандартного образца - по каждому изотопу для образцов одного содержания плутония.
Случайная составляющая погрешностей измерения по каждому изотопу для образцов одного содержания плутония.
Систематическая составляющая погрешностей измерения по каждому изотопу для образцов одного содержания плутония.
Доверительные границы (интервал) случайной погрешности результата измерений по каждому изотопу для образцов одного содержания плутония.
Верхняя граница относительного среднего квадратического отклонения результатов определений по каждому изотопу для образцов одного содержания плутония.
Доверительные границы погрешности результата анализа по каждому изотопу для образцов одного содержания плутония.
Доверительные границы погрешности с учетом неоднородности образцов по каждому изотопу для образцов одного содержания плутония.

Рассчитанные значения параметров калибровочных характеристик y=ax+b, где y - измеряемое значение содержания изотопа Pu, а x - полученное в аналитической лаборатории, приведены в таблицах ниже.

Таблица П1.3

Характеристики статистических моделей для измерений с декларацией содержания Pu-242

Изотоп Характеристика	Pu-238	Pu-239	Pu-240	Pu-241
Параметр a (стандартная ошибка параметра).	1,023 ()	1,007 ()	1,007 ()	0,992 ()
Параметр b (стандартная ошибка параметра).	()	()	()	()
Обусловленная регрессией сумма квадратов .	1,72			2,56
Остаточная сумма квадратов .		0,746	0,730
Коэффициент корреляции r.	0,998	0,999	0,999	0,999

Таблица П1.4

Характеристики статистических моделей для измерений без декларации содержания Pu-242

Изотоп Характеристика	Pu-238	Pu-239	Pu-240	Pu-241	Pu-242
Параметр a (стандартная ошибка параметра).	1,021 ()	1,034 ()	1,002 ()	0,987 ()	2,259 (0,338)
Параметр b (стандартная ошибка параметра).	()	-3,032 (0,244)	0,017 (0,017)	()	-0,152 (0,066)
Обусловленная регрессией сумма квадратов .	0,214	1465	1278	0,316	1,473
Остаточная сумма квадратов		0,449	0,085		1,386
Коэффициент корреляции r.	0,999	0,999	0,999	0,999	0,718

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

При сопоставлении характеристик погрешности результатов определения содержания изотопов в стандартных образцах отмечается ухудшение большинства этих характеристик при использовании алгоритма без декларации содержания Pu-242. Это наиболее заметно для Am-241.

Для всех изотопов Pu регрессионные модели для уровня значимости 0,05 оказались адекватны экспериментальным данным как в случае декларации содержания Pu-242, так и без таковой. Для Pu-242 качество регрессионной модели оказалось ниже, чем для других изотопов Pu. Угол наклона калибровочной характеристики для Pu-242 существенно отличается от 45°, а для Pu-239 эта характеристика существенно смещается вниз для случая программного определения содержания Pu-242 в образцах.

При анализе составляющих погрешностей измерения отмечается тенденция к уменьшению систематической составляющей для Pu-238 с увеличением содержания этого изотопа при использовании обоих алгоритмов обработки спектров. Случайная составляющая для Pu-239 уменьшается, а для остальных изотопов Pu увеличивается с увеличением содержания соответствующего изотопа в случае использования алгоритма обработки спектров с декларацией содержания Pu-242. Других тенденций к увеличению (уменьшению) составляющих погрешностей для уровня значимости 0,05 не выявлено.

Характеристики погрешности результатов определения содержания изотопов Pu, приведенные в настоящем приложении, а также выявленные особенности предлагаемых разработчиками программного обеспечения алгоритмов обработки спектров целесообразно учитывать при решении задач учета и контроля ЯМ.

Приложение N 2

к Положению по применению методов

математической статистики для учета и контроля

ядерных материалов,

утв. приказом Федеральной службы по экологическому,

технологическому и атомному надзору

от 14 сентября 2011 г. N 535

Порядок расчета объема выборки

Объем выборки задается двумя величинами: пороговым количеством ЯМ в единицах массы (G) и вероятностью обнаружения недостачи/излишка порогового количества ЯМ в безотносительных единицах (Р). Эти величины определяются в соответствии с Основными правилами учета и контроля ядерных материалов (НП-030-11). Результаты проверки целостности неразборных УЕ, а также состояния СКД документируются на предмет документального подтверждения отсутствия несанкционированного доступа.

Перед расчетом объема выборки определяется отношение G/x.

Расчет объема выборки (n) осуществляется следующим образом:

при величине , объем выборки (n) рассчитывается по формуле

, (П2.1)

N - количество УЕ в страте;

x - средняя масса ЯМ в одной УЕ;

- в формуле означает округление до ближайшего большего целого числа.

при величине G/x > N проверке подлежат все учетные единицы.

Определение объема выборки проверяемых пломб согласно пункту 25 настоящего Положения осуществляется с помощью уравнения ; где N - количество установленных пломб в ЗБМ. Здесь также означает округление до ближайшего большего целого числа. Ниже приведены результаты такого расчета для нескольких значений N.

Всего пломб, штук	Число пломб, работоспособность и состояние которых проверяется
10	10
20	19
50	35
100	46
1000	59
10000	60

Пример расчета объема выборки

Пусть в ЗБМ находятся три страты ЯМ, в каждой из которых содержится 1000 УЕ. Под стратой понимается совокупность отдельных УЕ с одинаковыми или близкими физическими характеристиками и химическим составом ЯМ для осуществления статистической выборки. На практике, в страты обычно выделяют: ТВС одной партии, либо однотипные ТВС; одинаковые в пределах технических условий УЕ для формирования критических (подкритических) сборок; партии одинаковых продуктов, на которые оформляется общий паспорт и т. д. В научно-исследовательских организациях, а также в организациях, где номенклатура используемых ЯМ велика, количество таких страт может оказаться значительным.

Пусть в нашем примере первая и вторая страты представлены одинаковыми УЕ, содержащими металлический высокообогащенный уран в алюминиевой оболочке с номерами, выполненными электроискровым карандашом. Подобным образом выполнены УЕ третьей страты, содержащие металлический низкообогащенный уран. В каждой из УЕ, содержащих высокообогащенный и низкообогащенный уран, находится по 1000 г . Пусть УЕ первой страты в данный МБП после использования в эксперименте были извлечены из критсборки, помещены в трубы как в контейнеры, трубы опечатаны после регистрации соответствующих учетных данных, а УЕ второй и третьей страт в данный МБП не использовались, и находились в опечатанных трубах с момента предыдущей инвентаризации. Таким образом, для первой страты вероятность обнаружения недостачи/излишка порогового количества ЯМ принимается равной 0,5, а для второй и третьей страт - равной 0,25, так как материалы второй и третьей страт находились еще и под системой наблюдения. Такой выбор вероятности обнаружения для каждой из страт рекомендуется делать при условии, что результаты проверки целостности неразборных УЕ, а также состояния СКД документируются на предмет документального подтверждения отсутствия несанкционированного доступа.

Поскольку отношение G/x для первой и второй страт 8000/1000=8, а для третьей 70000/1000=70 и во всех случаях (значение N для каждой из страт равно 1000), то для расчета объема выборки во всех случаях используем формулу (П2.1).

В соответствии с формулой (П2.1) и параметрами G и P для первой страты объем выборки составит:

.

Для второй страты:

, но, округляя до ближайшего большего целого, получаем .

Для третьей страты:

, но, округляя до ближайшего большего целого, получаем .

Таким образом, всего должно быть подвергнуто измерениям:

83 + 36 + 1 = 120 учетных единиц.

Приложение N 3

к Положению по применению методов

математической статистики для учета и контроля

ядерных материалов,

утв. приказом Федеральной службы по экологическому,

технологическому и атомному надзору

от 14 сентября 2011 г. N 535

Анализ характеристик выборочного распределения

В табл. П3.1 приведены значения коэффициента для уровня значимости и различных степеней свободы k.

Если случайная величина x распределена нормально с математическим ожиданием a, то по выборке объема n () можно найти доверительные границы для a следующим образом: (), где k=n-1.

Таблица П3.1

Значения коэффициента Стьюдента

	0,80	0,90	0,95	0,98	0,99	0,995	0,999
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 22 24 26 28 30 40 50 60 80 100 150 200 300 500	1,886 1,638 1,533 1,476 1,440 1,415 1,397 1,383 1,372 1,363 1,356 1,350 1,345 1,341 1,337 1,333 1,330 1,328 1,325 1,321 1,318 1,315 1,313 1,310 1,303 1,299 1,296 1,292 1,290 1,287 1,286 1,284 1,283 1,282	2,920 2,353 2,132 2,015 1,943 1,895 1,860 1,833 1,813 1,796 1,782 1,771 1,761 1,753 1,746 1,740 1,734 1,729 1,725 1,717 1,711 1,706 1,701 1,697 1,684 1,676 1,671 1,664 1,660 1,655 1,653 1,650 1,648 1,645	4,303 3,182 2,776 2,571 2,447 2,365 2,306 2,262 2,228 2,201 2,179 2,160 2,145 2,131 2,120 2,110 2,101 2,093 2,086 2,074 2,064 2,056 2,048 2,042 2,021 2,009 2,000 1,990 1,984 1,976 1,972 1,968 1,965 1,960	6,965 4,541 3,747 3,365 3,143 2,998 2,897 2,821 2,764 2,718 2,681 2,650 2,625 2,603 2,584 2,567 2,552 2,540 2,528 2,508 2,492 2,479 2,467 2,457 2,423 2,403 2,390 2,374 2,364 2,352 2,345 2,339 2,334 2,326	9,925 5,841 4,604 4,032 3,707 3,500 3,355 3,325 3,169 3,106 3,055 3,012 2,977 2,947 2,921 2,898 2,878 2,861 2,845 2,819 2,797 2,770 2,763 2,750 2,705 2,678 2,660 2,639 2,626 2,609 2,601 2,592 2,586 2,576	14,09 7,453 5,598 4,773 4,317 4,029 3,833 3,690 3,581 3,497 3,423 3,373 3,326 3,286 3,252 3,222 3,197 3,174 3,153 3,119 3,091 3,067 3,047 3,030 2,971 2,937 2,915 2,887 2,871 2,849 2,839 2,828 2,820 2,807	31,60 12,92 8,610 6,869 5,959 5,408 5,041 4,781 4,587 4,437 4,318 4,221 4,141 4,073 4,015 3,965 3,922 3,883 3,850 3,792 3,745 3,707 3,674 3,646 3,551 3,496 3,460 3,416 3,391 3,357 3,340 3,323 3,310 3,291

Пример.

Масса изделий подчиняется нормальному закону. Имеется статистический материал для n=25 образцов, по которым найдено г., а s = 600 г. Найти для доверительные границы для математического ожидания массы генеральной совокупности образцов.

Решение. По табл. П3.1 для и k = n - 1 = 24 находим , и определяем доверительные границы для математического ожидания следующим образом:

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

Приложение N 4

к Положению по применению методов

математической статистики для учета и контроля

ядерных материалов,

утв. приказом Федеральной службы по экологическому,

технологическому и атомному надзору

от 14 сентября 2011 г. N 535

Рекомендуемый порядок применения контрольного допуска

С помощью гамма - спектрометрической станции U-Pu InSpector с низкоэнергетическим германиевым детектором измерен образец с процентным содержанием Pu-239, равным по данным разрушающего анализа (учетные данные). Для получения учетных данных согласно требованию Основных правил учета и контроля ядерных материалов (НП-030-11), методики выполнения измерений должны быть метрологически аттестованы в соответствии с требованиями документов в области стандартизации (границы погрешности измерений устанавливаются для доверительной вероятности 95%).

Для получения данных об изотопном составе с помощью гамма - спектрометрической станции U-Pu InSpector с низкоэнергетическим германиевым детектором спектр обрабатывается с помощью программного обеспечения станции MGA. Для определения границ погрешности измерения использовалась методика: Плутоний и его соединения. Методика выполнения измерений массовых долей изотопов плутония и изотопа америций-241 в плутонии гамма-спектрометрическим методом с использованием гамма-спектрометра U-Pu InSpector. МВИ 223.13.17.104/2006 (границы погрешности измерений по этой методике устанавливаются для доверительной вероятности 95%). Результат, полученный с использованием названной методики, оказался равным: .

Поскольку все погрешности документируются и рассчитываются при 95% доверительной вероятности, то для получения среднеквадратических погрешностей значения 0,06 и 0,05 должны быть уменьшены в 1,96 раза. Таким образом, .

Расхождение учетных данных и результата подтверждающих измерений составляет 0,01 и находится в пределах контрольного допуска 0,10.

Примеры сравнения двух выборочных дисперсий и средних

При необходимости объединения данных по двум или нескольким партиям ЯМ рекомендуется проводить сравнение выборочных дисперсий и средних. Случай двух партий рассмотрен ниже. В результате проведения учетных измерений двух партий, в состав которых входят и УЕ, найдены выборочные значения математического ожидания и дисперсии массы учитываемого нуклида в этих УЕ. Значения этих величин составили: кг; кг; ; . Требуется оценить значимость расхождения выборочных дисперсий.

В рассматриваемом примере при числе степеней свободы 29 и 19. Это говорит об отсутствии значимого различия в дисперсиях массы учитываемого нуклида в УЕ первой и второй выборок. То есть можно принять, что данные выборки могут быть объединены по критерию равенства дисперсий после проверки равенства средних.

Для этого определим объединенную оценку выборочной дисперсии:

и значение t статистики:

Т.к. при числе степеней свободы 48, то можно принять, что , а значение объединенной оценки выборочного среднего равно:

кг.

Приложение N 5

к Положению по применению методов

математической статистики для учета и контроля

ядерных материалов,

утв. приказом Федеральной службы по экологическому,

технологическому и атомному надзору

от 14 сентября 2011 г. N 535

Примеры расчета , а также алгоритм оценки значимости систематических расхождений

Примеры расчета :

Допустим, полные погрешности весов, приведенные к доверительной вероятности 0,99 составляют:

а) отправителя г, получателя г; тогда

г;

б) отправителя %, получателя %, а измеряемая масса 25 кг; тогда

кг;

в) отправителя кг, получателя %, а измеряемая масса 25 кг; тогда

кг.

Алгоритм оценки значимости систематических расхождений

Для проверки значимости систематических расхождений между данными отправителя и получателя рекомендуется пользоваться парным критерием Стьюдента. Для этого вычислить:

- наблюдаемые парные разности

- среднее значение парных разностей , где n - количество пар

- среднее квадратическое отклонение среднего значения парных разностей

- наблюдаемый критерий .

При выполнении условия , где табличное значения коэффициента Стьюдента при доверительной вероятности P и числе степеней свободы f=n-1, нулевая гипотеза отвергается.

Пример. В таблице ниже приведены данные масс брутто отправителя и получателя пятнадцати контейнеров с ЯМ. Допустимые расхождения, вычисленные по паспортным данным, не должны превышать 15 условных единиц.

Таблица

N п/п						Z
1	2	3	4	5	6	7
1	282	290	-8	-0,3	0,09	286,0
2	347	350	-3	4,7	22,09	348,5
3	286	296	-10	-2,3	5,29	291,0
4	319	326	-7	0,7	0,49	322,5
5	337	346	-9	-1,3	1,69	341,5
6	290	300	-10	-2,3	5,29	295,0
7	301	300	1	8,7	75,69	300,5
8	314	323	-9	-1,3	1,69	318,5
9	373	380	-7	0,7	0,49	376,5
10	236	246	-10	-2,3	5,29	241,0
11	361	370	--9	-1,3	1,69	365,5
12	333	340	-7	0,7	0,49	336,5
13	292	302	-10	-2,3	5,29	297,0
14	340	347	-7	0,7	0,49	343,5
15	320	330	-10	-2,3	5,29	325,0
	4731	4846	-115	0,5	131,35	4788,5

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

Проверим значимость систематических расхождений.

По данным столбцов 2 и 3 вычисляем и заполняем столбец 5.

По данным столбца 4 вычисляем .

По данным столбца 4 вычисляем и , заполняем столбцы 5 и 6 и вычисляем .

Вычисляем наблюдаемый критерий .

Критическое значение коэффициента Стьюдента при доверительной вероятности 0,99 и числе степеней свободы 14 равно 2,9768.

Так как условие соблюдается, то нулевая гипотеза отвергается. Следовательно, между данными отправителя и получателя наблюдаются значимые систематические расхождения, и до выяснения причин расхождений взаиморасчеты рекомендуется проводить по средним значениям Z, значения которых приведены в столбце 7 таблицы.

Приложение N 6

к Положению по применению методов

математической статистики для учета и контроля

ядерных материалов,

утв. приказом Федеральной службы по экологическому,

технологическому и атомному надзору

от 14 сентября 2011 г. N 535

Аддитивная и мультипликативная модель погрешности измерения. Связь с абсолютной и относительной погрешностью измерения

П6.1 Абсолютная погрешность измерения

По определению абсолютная погрешность есть

(П6.1)

где: - измеренная величина;

- истинное значение.

Абсолютная погрешность выражается в единицах измеряемой величины. Значение установленной погрешности для МВИ или средств измерений является неизменной во всем измеряемом (установленном) диапазоне. Например, грамм для весов означает, что во всем измеряемом диапазоне эта декларируемая погрешность остается постоянной.

Взвешивание 1 кг осуществляется с погрешностью 0,25 г, взвешивание 10 кг также осуществляется с погрешностью 0,25 г.

График зависимости погрешности измерения от измеряемой величины показан на рис. П6.1.

В общем случае значение абсолютной погрешности не зависит от значения измеряемой величины (свойство аддитивности). На практике производитель средства измерения может определить различные значения абсолютной погрешности для разных диапазонов измеряемых значений (например: от 1 до 100 грамм - 0,5 грамм, а от 100 до 500 грамм - 1 грамм).

П6.2 Относительная погрешность измерения

По определению относительная погрешность есть

(П6.2)

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

Значение установленной погрешности для МВИ или средства измерений является неизменной во всем измеряемом (установленном) диапазоне. Например, % для весов означает, что во всем измеряемом диапазоне эта декларируемая погрешность остается постоянной. Взвешивание 1кг осуществляется с погрешностью 0,05%, что составляет 0,5 г, взвешивание 10 кг также осуществляется с погрешностью 0,05%, что составляет 5 г. В данном случае значение абсолютной погрешности прямо пропорционально измеряемой величине (свойство мультипликативности).

График зависимости погрешности измерения от измеряемой величины показан на рис. П6.2.

П6.3 Аддитивная модель погрешности измерения, cвязь с абсолютной погрешностью измерения

Абсолютная погрешность измерения структурируется и представляется как сумма систематической и случайной составляющих погрешности

(П6.3)

где: S - систематическая составляющая погрешности, имеет размерность измеряемой величины;

R - случайная составляющая погрешности, имеет размерность измеряемой величины.

Подставив правую часть выражения (П6.3) в (П6.1) и сделав преобразования, получим:

(П6.4)

Выражение (П6.4) есть аддитивная модель погрешности измерения. Это есть другая запись абсолютной погрешности измерения с учетом выделения систематической и случайной составляющих.

Свойство модели (П6.4) со статистической точки зрения:

- детерминированная величина;

S - случайная величина, как правило подчиняющаяся нормальному закону распределения ;

R - случайная величина, как правило подчиняющаяся нормальному закону распределения ;

Дисперсия погрешности из (П6.4) будет:

(П6.5)

где: - СКО систематической составляющей погрешности (имеет размерность измеряемой величины);

- СКО случайной составляющей погрешности (имеет размерность измеряемой величины).

Значения и необходимо вычислить по известным интервальным оценкам погрешности. Значения и будут использоваться при вычислении .

П6.4 Мультипликативная модель погрешности измерения, cвязь с относительной погрешностью измерения.

Относительная погрешность измерения структурируется и представляется как сумма систематической и случайной составляющих погрешности:

(П6.6)

где: S - систематическая составляющая погрешности, безразмерная величина;

R - случайная составляющая погрешности, безразмерная величина.

Подставив правую часть выражения (П6.6) в (П6.2) и сделав преобразования, получим:

(П6.7)

Выражение (П6.7) есть мультипликативная модель погрешности измерения. Это есть другая запись относительной погрешности измерения с учетом выделения систематической и случайной составляющих.

Свойство модели (П6.7) со статистической точки зрения:

- детерминированная величина;

S - случайная величина, как правило подчиняющаяся нормальному закону распределения ;

R - случайная величина, как правило подчиняющаяся нормальному закону распределения .

Дисперсия погрешности из (П6.7) будет:

(П6.8)

где: - СКО систематической составляющей погрешности (безразмерная величина);

- СКО случайной составляющей погрешности (безразмерная величина).

Значения и необходимо вычислить по известным интервальным оценкам погрешности. Значения и будут использоваться при вычислении .

П6.5 Смешанная модель и приведенная погрешность.

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

Принято считать, что для средства измерения на участке диапазона протяженностью до 1% от его предела модель погрешности можно считать аддитивной. Это допущение позволяет существенно упростить оценки погрешностей с достаточной степенью корректности.

Часто для описания погрешности средства измерения применяют приведенную погрешность, которая равна максимальному значению погрешности по диапазону, отнесенному к пределу измерения (верхнему пределу шкалы). По приведенной погрешности можно оценить максимальную погрешность в диапазоне независимо от ее модели.

Приложение N 7

к Положению по применению методов

математической статистики для учета и контроля

ядерных материалов,

утв. приказом Федеральной службы по экологическому,

технологическому и атомному надзору

от 14 сентября 2011 г. N 535

Характеристики погрешности измерений*

Таблица N П7.1 - Характеристики погрешности измерений ЯМ

Название ЗБМ: Производство таблеток

Группа (страта) ЯМ	Метод измерения	Случайная составляющая погрешности			Систематическая составляющая погрешности
1 Порошок двуокиси урана 4,4%	1 Взвешивание Mettler KE 1500 2 Масс. доля U ОСТ 95 175-90 3 Масс доля ОСТ 95 446-84	- - -	- 0,0014 0,0075	0,0005 отн. ед 0,00071 отн. ед 0,0038 отн. ед		- 0,0016 0,0042	0,0005 отн. ед 0,0008 отн. ед 0,0021 отн. ед
2 Шлиф. отходы 4,4%	1 Взвешивание Mettler KB-60s 2 Масс. доля U ТУ 002.46-82 3 Масс доля ОСТ 95 446-84	- - -	- 0,0071 0,0075	0,0005 отн. ед 0,0036 отн. ед 0,0038 отн. ед	- - -	- 0,0081 0,0042	0,0005 отн. ед 0,0041 отн. ед 0,0021 отн. ед
3 Таблетки 4,4%	1 Взвешивание Mettler KB-32s 2 Масс. доля U ОСТ 95 175-90 3 Масс доля ОСТ 95 446-84	- - -	- 0,0014 0,0075	0,0005 отн. ед 0,00071 отн. ед 0,0038 отн. ед	- - -	- 0,0016 0,0042	0,0005 отн. ед 0,0008 отн. ед 0,0021 отн. ед
4 Отходы спеченные 4,4%	1 Взвешивание Mettler KB-60s 2 Масс. доля U ОСТ 95 175-90 3 Масс доля ОСТ 95 446-84	- - -	- 0,0014 0,0075	0,0005 отн. ед 0,00071 отн. ед 0,0038 отн. ед	- - -	- 0,0016 0,0042	0,0005 отн. ед 0,0008 отн. ед 0,0021 отн. ед

_____________________________

* Статистические методы для контроля качества измерений и анализа инвентаризационной разницы. Методические материалы курса. Учебно-методический центр по учету и контролю ядерных материалов. Обнинск, ГНЦ РФ-ФЭИ, 2010.

Приложение N 8

к Положению по применению методов

математической статистики для учета и контроля

ядерных материалов,

утв. приказом Федеральной службы по экологическому,

технологическому и атомному надзору

от 14 сентября 2011 г. N 535

Таблицы движения материалов в ЗБМ

Таблица П8.1 - Таблица движения материалов в ЗБМ производства ТВЭЛ

Наименование материала	Номер партии	Кол-во УЕ	Масса брутто	Масса тары	Масса нетто	Масс. доля эл-та	Масс. эл-та	Масс. доля изотопа	Масса изотопа
Предыдущая ФИ (СФНК)
Таблетки	Таб05 Таб07 Таб10 Таб12 Таб15	75 39 17 45 53
ТВЭЛ	ТВ07 ТВ09 ТВ11	135 85 37
Брак (таблетки)	Таб07d Таб10d	2 3
Сметки	См01	2
Поступления и др. увеличения ЯМ в ЗБМ (УВ)
Таблетки	Таб16 Таб17 Таб18 Таб19 Таб20 Таб21 Таб22	205 203 200 203 229 211 201
Отправления и др. уменьшения ЯМ из ЗБМ (УМ)
Брак (таблетки)	Таб07d Таб10d Таб35d Таб36d Таб37d	2 3 5 7 3
ТВЭЛ	ТВ07 ТВ09 ТВ11 ТВ15 ТВ17	120 83 37 58 115
Сметки	См01 См02	2 4
Текущая ФИ (СФНК)
Таблетки	Таб12 Таб17 Таб20 Таб22	45 53 29 51
ТВЭЛ	ТВ07 ТВ09 ТВ15	15 2 35
Брак -таб	-	0
Сметки	См03	2

Таблица П8.2 - Группа (страта) ТВЭЛ

Название ЗБМ: Производство ТВЭЛ

Наименование материала	Номер партии	Кол-во УЕ	Масса брутто	Масса тары	Масса нетто	Масс. доля эл-та	Масса эл-та	Масс. доля изотопа	Масса изотопа
Предыдущая ФИ (СФНК)
ТВЭЛ	ТВ07 ТВ09 ТВ11	135 85 37
Поступления и др. увеличения ЯМ в ЗБМ(УВ)
ТВЭЛ	-	0
Отправления и др. уменьшения ЯМ из ЗБМ(УМ)
ТВЭЛ	ТВ07 ТВ09 ТВ11 ТВ15 ТВ17	120 83 37 58 115
Текущая ФИ (СФНК)
ТВЭЛ	ТВ07 ТВ09 ТВ15	15 2 35

Таблица П8.3 - Группа (страта) сметок

Название ЗБМ: Производство ТВЭЛ

Наименование материала	Номер партии	Кол-во УЕ	Масса брутто	Масса тары	Масса нетто	Масс. доля эл-та	Масса эл-та	М. Д. изотопа	Масса изотопа
Предыдущая ФИ (СФНК)
Сметки	См01	2
Поступления и др. увеличения ЯМ в ЗБМ (УВ)
Сметки		0
Отправления и др. уменьшения ЯМ из ЗБМ (УМ)
Сметки	См01 См02	2 4
Текущая ФИ (СФНК)
Сметки	См03	2

Таблица П8.4 - Группа (страта) таблеток (бракованных)

Название ЗБМ: Производство ТВЭЛ

Наименование материала	Номер партии	Кол-во УЕ	Масса брутто	Масса тары	Масса нетто	М.Д. эл-та	Масса эл-та	М. Д. изотопа	Масса изотопа
Предыдущая ФИ (СФНК)
Таблетки	Таб05 Таб07 Таб10 Таб12 Таб15	75 39 17 45 53
Брак (таблетки)	Таб07d Таб10d	2 3
Поступления и др. увеличения ЯМ в ЗБМ (УВ)
Таблетки	Таб16 Таб17 Таб18 Таб19 Таб20 Таб21 Таб22	205 203 200 203 229 211 201
Отправления и др. уменьшения ЯМ из ЗБМ (УМ)
Брак (таблетки)	Таб07d Таб10d Таб35d Таб36d Таб37d	2 3 5 7 3
Текущая ФИ (СФНК)
Таблетки	Таб12 Таб17 Таб20 Таб22	45 53 29 51
Брак -таб	-	0

Приложение N 9

к Положению по применению методов

математической статистики для учета и контроля

ядерных материалов,

утв. приказом Федеральной службы по экологическому,

технологическому и атомному надзору

от 14 сентября 2011 г. N 535

Пример исключения парных записей

Группа (страта) таблеток производство ТВЭЛ представлена в таблице П8.4 приложения N 8 к настоящему Положению.

В данном примере партия Таб12 (45 контейнеров) присутствует в предыдущей (НК) и текущей инвентаризации (ФК), партии Таб17, Таб20, Таб22 поступили в ЗБМ (УВ) и частично были переработаны - собраны в ТВЭЛы, оставшиеся контейнеры из этих партий присутствуют в текущей инвентаризации (ФК). Данные по контейнерам в этих партиях представляют собой "парные учетные записи". Остальные партии таблеток были собраны в ТВЭЛы.

Партии Таб07d, Таб10 (бракованные таблетки), из предыдущей (НК) инвентаризации были отправлены из ЗБМ (УМ). УЕ данных партий представляют собой "парные учетные записи".

После исключения "парных учетных записей" данная группа (страта) представлена в таблице N П9.1.

Таблица П9.1 - Группа(страта) таблеток (исключены "парные учетные записи")

Название ЗБМ: Производство ТВЭЛ

Наименование материала	Номер партии/	Кол-во УЕ	Масса брутто	Масса тары	Масса нетто	Масс. доля эл-та	Масса эл-та	Масс. доля изотопа	Масса изотопа
Предыдущая ФИ (СФНК)
Таблетки	Таб05 Таб07 Таб10 Таб15	75 39 17 53
Поступления и др. увеличения ЯМ в ЗБМ (УВ)
Таблетки	Таб16 Таб17 Таб18 Таб19 Таб20 Таб21 Таб22	205 150 200 203 200 211 150
Отправления и др. уменьшения ЯМ из ЗБМ (УМ)
Брак (таблетки)	Таб35d Таб36d Таб37d	5 7 3
Текущая ФИ (СФНК)
Таблетки

Приложение N 10

к Положению по применению методов

математической статистики для учета и контроля

ядерных материалов,

утв. приказом Федеральной службы по экологическому,

технологическому и атомному надзору

от 14 сентября 2011 г. N 535

Оценка составляющей дисперсии, на которую уменьшается вычисляемая по способу с усреднением по потоку

Если обозначить сумму массы элемента (изотопа) в n УЕ партии как , то средняя масса ЯМ одной УЕ будет .

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

Если найти разность между аналитическими выражениями вычисления по способу без усреднения по потоку и по способу с усреднением по потоку (для мультипликативных моделей погрешности измерения), то составляющая, на которую уменьшается оценка по методу усреднения по потоку, есть:

(П10.1)

где: - масса ЯМ каждой УЕ партии;

- масса ЯМ всей партии;

n - число УЕ в партии;

- СКО случайной составляющей погрешности взвешивания;

- СКО случайной составляющей погрешности МВИ разрушающего метода анализа.

Пример:

Пусть имеется партия ЯМ состоящая из 50 УЕ, измеренные массы которых представлены ниже. Погрешности измерений: %, %.

41.39 40.96 40.66 40.74 41.6  40.87 41.26 41.46 40.77 41.29

41.85 41.22 40.83 40.44 41.24 40.96 41.3  41.87 41.32 41.09

42.33 40.99 41.5  41.02 41.82 40.36 40.9  41.04 40.85 41.34

41.72 41.62 41.92 39.74 41.36 41.08 41.14 40.88 41.62 41.12

40.93 42.2  41.14 41.33 40.7  41.35 41.2  41.1  41.31 41.58

Средняя масса ЯМ одной УЕ = 41.206

Составляющая, на которую уменьшается оценка по методу усреднения по потоку, есть:

Приложение N 11

к Положению по применению методов

математической статистики для учета и контроля

ядерных материалов,

утв. приказом Федеральной службы по экологическому,

технологическому и атомному надзору

от 14 сентября 2011 г. N 535

Аналитические выражения для вычисления без усреднения по потоку

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

П11.1 Вклад систематических составляющих погрешностей в .

Вклад систематической составляющей погрешности измерения по одному из методов измерения (взвешивания, измерения объема, МВИ разрушающих и неразрушающих методов анализа) в в одной из групп (страт) вычисляется (без усреднения по потоку):

(П11.1)

где: для из НК, УВ - членов уравнения ИР;

для из ФК, УМ - членов уравнения ИР;

- СКО систематической составляющей погрешности одного из методов измерения.

Этот подход справедлив, например, при взвешивании всего ЯМ, составляющего члены уравнения баланса, на одних и тех же весах без перекалибровки.

Вклад в дисперсию систематической составляющей погрешности взвешивания будем обозначать с индексом "M": и . При измерении объема в обозначениях вместо "M" следует подставить "V". Вклад в дисперсию систематической составляющей погрешности МВИ разрушающих и неразрушающих методов анализа будем обозначать с индексом "A": и .

При проведении взвешиваний и определении массовой доли элемента (изотопа) рекомендуется следующий подход.

При вычислении ИР по элементу:

, (П11.2)

где: - масса брутто;

- масса тары;

- массовая доля элемента.

При вычислении ИР по изотопу:

, (П11.3)

где: - массовая доля изотопа.

П11.2 Вклад случайных составляющих погрешностей в

П11.2.1 Вклад случайной составляющей погрешности взвешивания в в одной из групп (страт) (без усреднения по потоку):

(П11.4)

где: - количество взвешиваний i-го контейнера в j-той группе (страте) (требование записи среднего результата измерений по взвешиваниям одного и того же контейнера определяется инструкциями предприятия), в большинстве случаев ;

- СКО случайной составляющей погрешности взвешивания в j-той группе (страте);

- массовая доля элемента.

При вычислении по изотопу вместо необходимо в (П11.4) подставить .

Выражение (П11.4) используется при способе расчета без усреднения по потоку, в котором , индивидуальные массы брутто и тары i-го контейнера в j-той группе (страте).

П11.2.2 Вклад случайной составляющей погрешности измерения объема в в одной из групп (страт) (без усреднения по потоку):

(П11.5)

где: - количество измерений i-ого объема (требование записи среднего результата измерений по результатам измерения одного и того же объема определяется инструкциями предприятия), в большинстве случаев ;

- СКО случайной составляющей погрешности измерения объема в j-той группе (страте);

- массовая доля элемента.

При вычислении по изотопу вместо необходимо в (П11.5) подставить .

Выражение (П11.5) используется при способе расчета без усреднения по потоку, в котором индивидуальные измерения объема в j-той группе (страте).

П11.2.3 Вклад в случайной составляющей погрешности МВИ одного из методов разрушающего и неразрушающего анализа в j-той группе (страте) (без усреднения по потоку):

(П11.6)

где: R - количество партий ЯМ, в которых определяется массовая доля элемента (изотопа) по МВИ (для j-той группы (страты);

- число контейнеров в партии;

- массы брутто, тары k-го контейнера из i- той партии;

- число проб на партию;

- массовая доля элемента;

- СКО случайной составляющей погрешности МВИ для j-той группы (страты).

При вычислении по изотопу вместо необходимо в (П11.6) подставить .

Приложение N 12

к Положению по применению методов

математической статистики для учета и контроля

ядерных материалов,

утв. приказом Федеральной службы по экологическому,

технологическому и атомному надзору

от 14 сентября 2011 г. N 535

Пример вычисления суммарной массы материала

В таблице П12.1 представлены учетные данные предыдущего примера по партиям таблеток производства ТВЭЛ и результаты вычислений суммарной массы :

кг (для элемента)

кг (для изотопа)

кг (для элемента)

кг (для изотопа)

Таблица П12.1 - Суммарные массы и , исключены "парные учетные записи"

Название ЗБМ: Производство ТВЭЛ

Группа(страта) таблеток (исключены "парные учетные записи")
Член ур-я ИР	Партия	Кол-во УЕ	Масса нетто кг	м.д. U %	Масса кг	м.д. %	Масса кг
НК	Таб05 Таб07 Таб10 Таб15	75 39 17 53	3098,25 1615,38 706,01 1859,40	87,3 87,5 87,4 87,2	2704,7723 1413,4575 617,05274 1621,3968	2,91 2,87 2,95 2,97	78,708874 40,56623 18,203056 48,155485
	Сумма				6356,6793		185,63364
УВ	Таб16 Таб17 Таб18 Таб19 Таб20 Таб21 Таб22	205 150 200 203 200 211 150	8513,65 6202,50 8286,00 8426,53 8284,45 8754,39 6234,62	87,6 87,7 87,5 87,4 87,3 87,8 87,6	7457,9574 5439,5925 7250,25 7364,7872 7232,3249 7686,3544 5461,5271	2,94 2,95 2,98 2,93 2,95 2,97 2,94	219,26395 160,46798 216,05745 215,78826 213,35358 228,28473 160,5689
	Сумма				47892,794		1413,7849
УМ	Таб35d Таб36d Таб37d	5 7 3	172,85 242,9 105,45	87,5 87,6 87,4	151,24375 212,7804 92,1633	2,98 2,94 2,93	4,5070637 6,2557438 2,7003847
	Сумма	0			456,18745		13,463192
ФК	-	0			0		0
	Сумма	0			0		0
					53793,286		1585,9553
					54705,661		1612,8817

Приложение N 13

к Положению по применению методов

математической статистики для учета и контроля

ядерных материалов,

утв. приказом Федеральной службы по экологическому,

технологическому и атомному надзору

от 14 сентября 2011 г. N 535

Пример вычисления вклада в систематических и случайных составляющих погрешностей группы (страты)

В таблице П13.1 представлены результаты вычисления по данным Приложения N 12 настоящего Положения вклада в систематических и случайных составляющих погрешностей группы (страты).

Для рассматриваемого примера выполняется однократное взвешивание каждого контейнера, поэтому соответствует суммарному числу контейнеров (1518). Массовая доля элемента и изотопа определялась в другой ЗБМ. В среднем отбирались 2 пробы на 200 кг ЯМ для определения массовой доли элемента и 2 пробы на 300 кг ЯМ для определения массовой доли изотопа. Каждая проба анализировалась один раз, поэтому число анализов для определения массовой доли элемента есть , число анализов для определения массовой доли изотопа есть .

Для взвешивания СКО систематической составляющей погрешности % (соответствует международным целевым значениям ITV - 2000). Для МВИ разрушающего анализа СКО систематической составляющей погрешности определения массовой доли элемента %, определения массовой доли изотопа %.

Вклад систематических составляющих погрешностей в :

- - взвешивание (723,4294 - по элементу, 0,62881 - по изотопу);

- - разрушающий анализ (2893,7176 - по элементу, 1,60976 - по изотопу);

- общий вклад систематических составляющих погрешностей в j-ой группы (страты) (3617,147 - по элементу, 2,23857 - по изотопу).

Для взвешивания случайная составляющая погрешности % (соответствует международным целевым значениям ITV - 2000).

Для МВИ разрушающего анализа случайная составляющая погрешности определения массовой доли элемента %, определения массовой доли изотопа %.

Вклад случайных составляющих погрешностей в :

- - взвешивание (0,966 - по элементу; 0,00084 - по изотопу);

- - разрушающий анализ (12,31005 - по элементу; 1,60976 - по изотопу);

- общий вклад систематических составляющих погрешностей j-ой группы (страты) в (3617,147 - по элементу; 0,23649 - по изотопу).

Таблица П13.1 - Вклад систематических и случайных составляющих группы (страты) таблеток в .

Название ЗБМ: Производство ТВЭЛ

Группа (страта) таблеток (исключены "парные учетные записи")
Член ур-я ИР	Партия	Кол-во УЕ	Масса нетто кг	м.д. U %	Масса кг	м.д. %	Масса кг
НК	Таб05 Таб07 Таб10 Таб15	75 39 17 53	3098,25 1615,38 706,01 1859,40	87,3 87,5 87,4 87,2	2704,7723 1413,4575 617,05274 1621,3968	2,91 2,87 2,95 2,97	78,708874 40,56623 18,203056 48,155485
	Сумма				6356,6793		185,63364
УВ	Таб16 Таб17 Таб18 Таб19 Таб20 Таб21 Таб22	205 150 200 203 200 211 150	8513,65 6202,50 8286,00 8426,53 8284,45 8754,39 6234,62	87,6 87,7 87,5 87,4 87,3 87,8 87,6	7457,9574 5439,5925 7250,25 7364,7872 7232,3249 7686,3544 5461,5271	2,94 2,95 2,98 2,93 2,95 2,97 2,94	219,26395 160,46798 216,05745 215,78826 213,35358 228,28473 160,5689
	Сумма				47892,794		1413,7849
ФК	0	0			0
	Сумма	0			0
УМ	Таб35d Таб36d Таб37d	5 7 3	172,85 242,9 105,45	87,5 87,6 87,4	151,24375 212,7804 92,1633	2,98 2,94 2,93	4,5070637 6,2557438 2,7003847
	Сумма	0			456,18745		13,463192
Определение вклада систематических составляющих погрешностей в
					Элемент		Изотоп
					53793,286		1585,9553
Взвешивание					0,05 %		0,05 %
					723,4294		0,62881
Анализ					0,1 %		0,08 %
					2893,7176		1,60976
					3617,147		2,23857
Определение вклада случайных составляющих погрешностей в
					54705,661		1612,8817
Взвешивание					0,07 %		0,07 %
Число контейнеров					1518		1518
					0,966		0,00084
Анализ					0,15 %		0,1 %
Число анализов					547		11
					12,31005		0,23649
					13,27605		0,23733