Аналитические методы контроля величины инвентаризационной разницы при подведении баланса ядерных материалов тема диссертации и автореферата по ВАК РФ 05.11.13, кандидат технических наук Ишков, Юрий Геннадьевич

• Ишков, Юрий Геннадьевич
• кандидат технических наук
• 2007

Оглавление диссертации кандидат технических наук Ишков, Юрий Геннадьевич

Введение.

1. Анализ существующих отечественных и зарубежных подходов к оцениванию величины инвентаризационной разницы.

1.1 Требования нормативных документов.

1.2 Российские методы статистической обработки результатов измерений.

1.3 Нормативные документы, применяемые в США для измерений в области учета и контроля ядерных материалов.

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

1.5 Подход МАГАТЭ к методам оценивания инвентаризационной разницы

1.6 Статистическое оценивание инвентаризационной разницы, представленное в нормативных документах США.

1.7 Сравнение национальных и международных подходов к оцениванию результатов измерений и инвентаризационных разностей.

1.8 Выводы.

2 Аналитические методы контроля инвентаризационных разниц ядерных материалов на основании погрешностей средств и методов измерений методом «переноса погрешностей».

2.1 Измерение массы, аддитивная модель погрешности.

2.2 Перекалибровка.

2.3 Мультипликативная модель.

2.4 Аддитивная и мультипликативные модели.

2.5 Учет корреляций.

2.6 Смешанная модель.

2.7 Измерения N учетных единиц.

2.8 Обобщение формул.

2.9 Анализ линеаризации дисперсии инвентаризационной разницы методом Монте - Карло.

2.10 Выводы.

3 Методы контроля инвентаризационных разниц на основании прогнозирования временных рядов.

3.1 Постановка задачи.

3.2 Выявление характеристик временных рядов.

3.2.1 Анализ на стационарность, сезонность.

3.2.2 Анализ рядов с использованием автокорреляционного анализа.

3.3 Методы прогнозирования временных рядов.

3.3.1 Методы, использующие скользящее среднее.

3.3.2 Экспоненциальное сглаживание.

3.3.3 Регрессионный анализ.

3.3.4 Метод Бокса-Дженкинса.

3.3.5 Расчет предикторов в методе Бокса-Дженкинса.

3.3.6 Сравнение методов прогнозирования.

3.5 Моделирование временных рядов методом экспоненциального сглаживания.

3.6 Методы адаптивного прогнозирования.

3.7 Выводы.

4 Определение универсальной модели прогнозирования исследуемых временных рядов.

4.1 Выбор модели.

4.2 Расчет предикторов моделей и выбор наилучшей модели методологии Бокса - Дженкинса.

4.3 Проверка выбранной модели на адекватность.

4.4 Расчет предикторов методами экспоненциального сглаживания и адаптивного прогнозирования.

4.5 Выводы.

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

Введение диссертации (часть автореферата) на тему «Аналитические методы контроля величины инвентаризационной разницы при подведении баланса ядерных материалов»

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

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

Актуальность осуществления контроля за величиной инвентаризационной разницы, направленного в конечном счете на обеспечение сохранности ядерных материалов, подтверждается наличием в федеральных нормах и правилах учета и контроля ядерных материалов НП-030-05, являющимися обязательными для выполнения всеми эксплуатирующими организациями и предприятиями требований к определению наличия или отсутствия аномалии в системе УиК ЯМ по результатам сравнения ИР и ее погрешности с заданными пороговыми значениями.

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

Значение ИР для каждого ЯМ (или каждого вида ЯМ) рассчитывается по формуле:

ИР = КК-ДК = КК-(НК + УВ-УМ), (1) где КК - фактически наличное количество ЯМ в зоне баланса ядерного материала (ЗБМ), определенное в результате данной физической инвентаризации на конец межбалансового периода (МБП);

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

НК-наличное количество ЯМ в ЗБМ, определенное и документально зарегистрированное на начало данного МБП в результате предыдущей физической инвентаризации;

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

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

Появление ИР объясняется тем, что в процессе технологической переработки ядерных материалов происходит их преобразование из одной физической формы в другую. Например, на сублиматных заводах получение гек-сафторида урана [1, стр.71] производится прямым фторированием в газопламенном реакторе закиси-окиси урана с последующей десублимацией урана и очисткой газов, на радиохимических заводах для переработки облученного топлива [1, стр.73] используется экстракционный пюрекс-процесс. Обогащение урана производится на заводах разделения изотопов [1,стр.97] центробежным методом.

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

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

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

Таким образом, объектом исследования являются: а) индивидуальные инвентаризационные разницы, определенные уравнением (1) и представляющие собой алгебраическую сумму компонент. Отдельные компоненты ИР является в общем случае произведением измерений массы/объема на концентрацию химического элемента и изотопного состава ЯМ. Компоненты ИР входят в уравнение баланса с различными знаками. Концентрация и изотопный состав, как правило, определяется от партии ЯМ.

Погрешности измерений массы, объема, концентрации имеют в общем случае случайные и неисключенные систематические составляющие, изменяющиеся по аддитивной или мультипликативной модели. Измерения могут проводиться одними и теми же или различными средствами измерений. Между измерениями существуют корреляции, которые необходимо учитывать; б) Последовательность значений инвентаризационных разниц, полученных с заданной периодичностью. Физические инвентаризации ядерных материалов, по результатам которых подводятся балансы ЯМ и определяются инвентаризационные разницы, осуществляются с установленной периодичностью от одного месяца до одного года. Наиболее важные с точки зрения контроля зоны баланса ядерных материалов, инвентаризируются ежемесячно. Это означает, что для анализа инвентаризационных разниц имеется дополнительный статистический материал, представляющий собой временной ряд значений инвентаризационных разниц. Объемы таких временных рядов достигают сотен значений.

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

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

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

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

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

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

Для реализации поставленной цели решались следующие задачи:

• вывод аналитических формул расчета дисперсии компонент и результирующей дисперсии значения инвентаризационной разницы ИР( в произвольном инвентаризационном периоде и Расчет интервальной оценки погрешности инвентаризационной разницы. Сравнение результатов, полученных по аналитическим формулам, с результатами моделирования методом Монте-Карло;

• анализ поведения последовательностей инвентаризационных разниц ИР( различных видов технологических процессов переработки ядерных материалов, как членов временных рядов. Выбор модели прогнозирования, которая адекватно описывает все исследуемые временные ряды;

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

• разработка решающих правил для анализа значения ИР1 в инвентаризационном периоде и используя найденные доверительные интервалы в задачах 2 и 3, рассматривая величину ИР в качестве функции случайных величин и в качестве члена временного ряда. При построении решающих правил необходимо использовать требования и пороговые значения, установленные федеральными нормами и правилами государственной системы учета и контроля ядерных материалов;

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

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

• исследованы модели поведения временных рядов инвентаризационных разниц с учетом их нестационарности и сезонности всех технологических процессов переработки ЯМ на Сублиматном, Разделения изотопов, Радиохимическом и Химико-металлургическом заводах СХК;

• впервые найдена автором модель прогнозирования (авторегрессия первого порядка метода авторегрессии - интегрированного скользящего среднего Бокса-Дженкинса), являющейся универсальной и адекватно описывающей имеющиеся на комбинате временные ряды инвентаризационных разниц;

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

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

Практическая ценность. На основании результатов, полученных в результате диссертационной работы, автором разработан и введен в действие на ФГУП «СХК» стандарт организации СТО 76-2007, в котором для выдачи заключения о том, выявлена или нет аномалия в учете ядерных материалов по результатам физической инвентаризации, приведена методика анализа ИР, как случайного значения, являющегося следствием коррелированных погрешностей средств и методов измерений и найденный универсальный метод прогнозирования временного ряда инвентаризационных разниц.

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

Полученные аналитические методы расчетов реализованы автором в виде шаблонов компьютерных программ, используемых на заводах ФГУП «СХК» при регламентных работах по подведению балансов ядерных материалов.

Стандарт организации введен в действие приказом по комбинату, акт внедрения приведен в приложении к диссертации.

Апробация результатов работы и публикации. Апробация результатов работы проводилась при подготовке учебных материалов и чтении лекций в Учебно-методическом центре подготовки специалистов по учету и контролю ЯМ (УМЦУК) Федерального агентства по атомной энергии в г. Обнинске для специалистов отрасли, занимающихся учетом и контролем ядерных материалов. Лекции изданы в сборниках методических материалов УМЦУК в 2001,2005 и в 2007 годах.

Методы расчета погрешности инвентаризационной разницы, использующие погрешности средств и методов измерений, были оформлены в виде нормативного документа комбината «Рекомендации. Оценивание статистической значимости инвентаризационной разницы на основании погрешностей средств и методов измерений. Методы расчета», введены приказом по Сибирскому химическому комбинату и используются на заводах комбината с 2002 года.

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

Материалы работы отдельными частями докладывались и изданы в материалах 7 научно-технической конференции СХК в 2003 году и конференции ТААЭП-2007, организованной СГТА.

Результаты, полученные автором, были опубликованы в журнале «Метрология» в 2003г. №11, с.33-39 и в 2005г. №12, с.39-45.

Автор входит в состав рабочей группы Росатома по разработке федеральных норм и правил в области государственного учета и контроля ядерных материалов. Результаты диссертационной работы и опубликованные совместные статьи с сотрудниками ГНЦ РФ-ФЭИ по оценке результатов физических инвентаризаций ядерных материалов были использованы при разработке федеральных норм и правил государственного учета и контроля ядерных материалов.

Структура и объем диссертации. Диссертационная работа изложена на 148 страницах машинописного текста, иллюстрируется 33 рисунками, 18 таблицами и состоит из введения, пяти глав, заключения, списка литературы из 83 наименований и 2 приложений.

Заключение диссертации по теме «Приборы и методы контроля природной среды, веществ, материалов и изделий», Ишков, Юрий Геннадьевич

4.5 Выводы

Исследованы десять временных рядов, являющихся инвентаризационными разницами основных технологических процессов переработки ЯМ на комбинате: сублиматного, разделения изотопов, радиохимического и химико- металлургического. Определена оптимальная модель - авторегрессия первого порядка метода Бокса-Дженкинса, адекватно прогнозирующая все анализируемые временные ряды.

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

• экспоненциального сглаживания для полиномиальной модели нулевого порядка (EWMA);

• экспоненциального сглаживания с учетом тренда -метод Хольта

Holt);

• экспоненциального сглаживания с учетом сезонности -метод Винтерса (Winters);

• двойного экспоненциального сглаживания (2xEWMA);

• адаптивного экспоненциального сглаживания -метод Тригга-Лича

Шона (TLS).

При сравнениях моделей и проверок на адекватность использовались автокорреляционные функции. Временные ряды разбивались на «обучающие» и «экзаменующие» выборки. Для проверки, изменяются ли модели'рядов со временем и сравнения различных рядов, применялся метод оценки точности называемый «методом складного ножа» (jackknife method) или «перекрестным анализом дееспособности метода» (cross-validation method).

Выбор универсальной модели основывался на анализе поведения временных рядов в последние 48-60 периодов, являющиеся наиболее значимыми в характере их поведения для целей учета. Модель проверялась на близость прогнозов в последние 6-12 периодов.

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

5 Разработка процедур и программного обеспечения для контроля инвентаризационных разниц методом авторегрессии первого порядка

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

Пусть имеется 1 значений временного ряда. Тогда алгоритм прогнозирования выглядит следующим образом:

• Находятся параметры модели. Для рассчитанного на предыдущем шаге значения коэффициента автокорреляции первого порядка строится ряд остатков модели;

• Модель проверяется на адекватность. Для первых 20 задержек (лагов) ряда рассчитываются коэффициенты автокорреляций и С* статистика Бокса-Льюнга. С* статистика сравнивается с коэффициентом % . Если статистика меньше табличного значения с доверительной вероятностью 95%, то принимается гипотеза о том, что остатки модели являются случайными некоррелированными и модель адекватна.

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

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

• Полученное в следующем 1+1 периоде значение временного ряда сравнивается с рассчитанными в предыдущем периоде доверительным интервалом. Если это значение находится в доверительном интервал, то делается заключение о том, что характер поведения временного ряда не изменился, следовательно аномалия не выявлена. В случае, если появляется необходимость в применении данного программного обеспечения для новых коротких временных рядов (п<20), то для них в алгоритм прогнозирования включается метод экспоненциального сглаживания с определением оптимального коэффициента сглаживания. Далее проверка на адекватность производится аналогично. Если модель становится не адекватной, то пересчитывается оптимальный коэффициент сглаживания, если повторно модель не адекватна, то работа программного обеспечения останавливается и выводится сообщение о необходимости пересматривать применяемые модели.

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

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

Заключение

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

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

При разложении в ряд Тейлора использовались величины первого порядка. Методом Монте-Карло проверялась не значимость округления, вызванного отбрасыванием членов порядка выше первого при разложении ИР в ряд Тейлора.

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

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

• сглаживания скользящими средними;

• регрессионные;

• экспоненциального сглаживания для полиномиальной модели нулевого порядка (EWMA);

• экспоненциального сглаживания с учетом тренда -метод Хольта (Holt);

• экспоненциального сглаживания с учетом сезонности -метод Винтерса (Winters);

• двойного экспоненциального сглаживания (2xEWMA);

• адаптивного прогнозирования -метод Тригга-Лича-Шона (TLS);

• авторегрессии - интегрированного скользящего среднего (Бокса-Дженкинса).

Проведена оценка методов и для дальнейшего исследования выведены аналитические формулы расчетов предикторов следующих методов прогнозирования: авторегрессии - интегрированного скользящего среднего (Бокса-Дженкинса) на основании решения уравнений автокорреляций Юла-Уокера, метода экспоненциального сглаживания Брауна и адаптивного прогнозирования Тригга-Лича-Шона.

Для определения оптимальных предикторов метода Бокса-Дженкинса использован метод анализа графиков функций автокорреляций и функции частной автокорреляции. Выбран метод проверки адекватности модели с помощью проверки теста «хи-квадрат» Q -статистики Льюнга-Бокса.

Для проверки качества построения моделей предложены показатели оценок ошибок модели: среднее абсолютное отклонение (MAD), среднеквадрати-ческая ошибка (MSE), средняя абсолютная ошибка в процентах (МАРЕ) и средняя процентная ошибка (МРЕ). Для дальнейшего поиска оптимального метода, указанные методы реализованы в виде компьютерных программ.

В результате проведенных в диссертации исследований определена оптимальная модель - авторегрессия первого порядка, адекватно прогнозирующая все анализируемые временные ряды. При сравнениях моделей и проверок на адекватность использовались автокорреляционные функции.

Временные ряды разбивались на «обучающие» и «экзаменующие» выборки. Для проверки, изменяются ли модели рядов со временем и сравнении различных рядов, применялся метод оценки точности называемый «методом складного ножа» (jackknife method) или «перекрестным анализом дееспособности метода» (cross-validation method).

Выбор универсальной модели основывался на анализе поведения временных рядов в последние 48-60 периодов, являющиеся наиболее значимыми в характере их поведения для целей учета. Модель проверялась на близость прогнозов в последние 6-12 периодов. Также был рассмотрен класс адаптивных моделей, в которых массы, с которыми в модели учитываются последние наблюдения, убывают по экспоненциальному закону. Определено, что модель - авторегрессия второго порядка является наиболее близкой к оптимальной.

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

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

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

Разработано программное обеспечение, использующее анализ временных рядов инвентаризационных разниц методом авторегрессии первого порядка Бокса-Дженкинса.

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

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

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

Список литературы диссертационного исследования кандидат технических наук Ишков, Юрий Геннадьевич, 2007 год

1. В.М.Лебедев Ядерный топливный цикл. Технологии, безопасность, экономика. - М.: Энергоатомиздат, 2005. - 316с.

2. Федеральные нормы и правила в области использования атомной энергии. Основные правила учета и контроля ядерных материалов. НП-030-05. Москва, 2005. - 27с.

3. ГОСТ 8.009-84 нормируемые метрологические характеристики средств измерений.

4. РД 50-453-84 Руководящий нормативный документ. Методические указания. Характеристики погрешности средств измерений в реальных условиях эксплуатации. Методы расчета.

5. ГОСТ 8-207-76. Прямые измерения с многократными наблюдениями. Методы обработки результатов наблюдений. Основные положения.

6. МИ 2083-90. Рекомендации. Государственная система обеспечения единства измерений. Измерения косвенные. Определение результатов измерений и оценивание их погрешности. Москва, 1991

7. РНЦ «Курчатовский институт». Модельная автоматизированная система учета, контроля и гарантий ядерныхз материалов для сложных ядерных установок (система «НУМАКС»). Технический проект.

8. IAEA Safeguards Statistical Concepts and Techniques. Fourth Revised Edition. IAEA/SG/SCT/4, International Atomic Energy Agency, 1989.

9. ANSI N15-36-1994. American National Standard for Nuclear Materials -Nondestructive Assay Measurement Control and Assurance.

10. ANSI N15.19-1989. Volume Calibration Methods for Nuclear Materials Control. 1990.

11. Taylor B.N. and C.E.Kuyatt. Guidelines for evaluating and Expressing the Uncertainty of Measurement Results. NIST Technical Note 1297.

12. Measurement Control Guide. U.S. Department of Energy. April 1993.

13. ANSI N 15/51-90 "Measurement Control Program Analytical Chemistry Laboratory". American Nation Standards Institute. INC., 1430 Broudway, New York, NY, 10018 (1984).

14. R. Picard, J.Hater, "MAWST" Matirials accounting with sequential testing, Version 2.0 June 1991. N-4/91-633.

15. Bowen, W.Michael and Carl A.Bennet. Statistical Methods for Nuclear Matirial Management. NUREG/CR-4604, PNL-5849, 1988.- 1000p.

16. Kern,E.A. and K.K.S.Pillay. September 1992. Variance Propagation Code User's Manual, N-4/91-978.

17. Albert M. Liebetrau, Dennis R. Weier. Statistical Methods for Calculation of Inventory Differences and Uncertainty Estimation. Joint meeting of Russian and American Statistical Working Groups. Moscow. VNIIA. 1998.

18. МИ 2552-99 Рекомендация. Применение «Руководства по выражению неопределенности измерений».- Санкт-Петербург, 1999

19. Ю.Г.Ишков «Оценивание статистической значимости инвентаризационной разницы на основании погрешностей средств и методов измерений. Методы расчета»// Рекомендации предприятия, ФГУП «СХК».26с.

20. Ю.Г.Ишков, Б.М.Кербель «Влияние неисключенных систематических погрешностей на дисперсию инвентаризационной разницы в учете ядерных материалов»// Журн. «Метрология».- 2003.- №11.- С.33-39.

21. Ю.Г.Ишков, Б.М.Кербель «Оценка погрешности инвентаризационной разницы в учете ядерных материалов методом Монте-Карло» // Журн. «Метрология»,-2005.-№12.-С.39-45.

22. Система учета и контроля ядерных материалов. Инвентаризационная разница. Определение и оценивание статистической значимости. СТО 76-2007// Северск, ФГУП «СХК», 2007.- 50с.

23. А.Я.Карпенко «Статистический анализ инвентаризационной разности в учете и контроле ядерных материалов».- Обнинск: ГЦИПК, 2001.- 151с.

24. Д.Э.Ханк, Д.У Уичерн, А.Дж. Райте «Бизнес прогнозирование», 7-е издание. М.: Изд. дом «Вильяме», 2003. 656с.

25. Лукашин Ю.П. Адаптивные методы краткосрочного прогнозирования.- М.: «Статистика», 1979. 245с.

26. Статистические методы для ЭВМ/ Под ред. К.Энслейна.- М.: «Наука», 1986. 151с.

27. Дж.Мердок Контрольные карты,- М.: «Финансы и статистика», 1986. 120с.

28. В.А.Колемаев, В.Н.Калинина Теория вероятностей и математическая статистика.- М.: ЮНИТИ-ДАНА, 2003. 252с.

29. Ю.В.Чуев Ю.Б.Михайлов, В.И.Кузьмин Прогнозирование количественных характеристик процессов.- М.: «Наука», 1975.- 362с.

30. П.В.Новицкий, И.А.Зограф Динамика погрешностей средств измерений.- Ленинград: Энергоатомиздат, 1990.-192с.

31. Статистические методы анализа экономической динамики.- М.: «Наука», 1983.-391с.

32. Статистика и управление случайными процессами/ Сб.статей под ред. А.Н.Ширяева.- М.: «Наука», 1989. -232с.

33. К.Д.Льюис Методы прогнозорования экономических показателей.-М.: Финансы и статистика, 1986

34. Дж.Теннант-Смит Бейсик для статистиков,- М.: «Мир», 1988. -276с.

35. Э.Леман Проверка статистических гипотез.- М.: «Наука», 1979. -349с.

36. С.А. Айвазян Прикладная статистика. Основы эконометрики, том 1, Том 2: М.: ЮНИТИ-ДАНА, 2001. -432с.

37. Статистические методы повышения качества: Пер. с англ. под ред. Х.Куме.- М.: «Финансы и статистика», 1990. 304с.

38. М.Дж.Кендалл, А.Стьюарт Многомерный статистический анализ и временные ряды.- М.: «Наука», 1976. -736с.

39. Е.М.Четыркин, И.Л.Калихман Вероятность и статистика.- М.: «Финансы и статистика», 1982. 319с.

40. С.А.Айвазян Статистическое исследование зависимостей.- М.: «Наука»,1968. 380с.

41. Дж.Бокс, Г.Дженкинс Анализ временных рядов. Прогноз и управление, Вып.1 и 2.- М.: Мир, 1974.-416с.

42. И.В.Никифоров Последовательное обнаружение изменения свойств временных рядов.- М.: «Наука», 1983. 320с.

43. Statistical Methods in Nuclear Material Control, by John L. Jaech. -US AEC/TID-26298,1973/ 227p.

44. Ермаков C.M. Метод Монте-Карло и смежные вопросы.- М.: Наука, 1971.-232с.

45. М.Кендэл Временные ряды: Пер. с англ. и предисловие Ю.П.Лукашина.- М.: «Финансы и статистика», 1981,- 199с.

46. Reid DJ. Forecasting in action: A comparison of forecasting techniques in economic time-series, 1971

47. Дж.Мэйндоналд Вычислительные алгоритмы в прикладной статистике.- М.: Финансы и статистика, 1988. 302с.

48. С.А. Айвазян, И.С. Енюков, Л.Д. Мешалкин, Прикладная статистика, Классификация и снижение размерности. -М.: Финансы и статистика, 1989.-312с.

49. А. Афифи, С. Эйзен, Статистический анализ. Подход с использованием ЭВМ. М.: Мир, 1982. - 215с.

50. Таблицы математической статистики: Л.Н. Большев, Н.В. Смирнов. -М. .-Наука, 1983.-416с.

51. Е.С. Вентцель, Л.А. Овчаров, Теория вероятностей и её инженерные приложения.- М.: Наука, 1988. 317с.

52. В.А. Грановский, Т.Н. Сирая, Методы обработки экспериментальных данных при измерениях. -JL: Энергоатомиздат, 1990. 129с.

53. Н. Джонсон, Ф. Лион, Статистика и планирование эксперимента в технике и науке, Методы планирования эксперимента.- М.: Мир, 1981.- 119с.

54. Н. Дрейпер, Г. Смит, Прикладной регрессионный анализ, -М.: Финансы и статистика, Книга 1, 1986. Книга 2, 1988. 392с.

55. А.Н. Зайдель, Ошибки измерений физических величин.- М.: Наука, 1974.- 108с.

56. С.Р. Рао, Линейные статистические методы и их применение. М.: Наука, 1968.-218с.

57. Д. Химмельблау, Анализ процессов статистическими методами. -М.: Мир, 1973.-957 с.

58. Д. Химмельблау, Обнаружение и диагностика неполадок в химических и нефтехимических процессах.- Л.: Химия, 1983, -351 с.

59. Douglas С. Montgomery, Introduction to Statistical Quality Control, Second Edition, New York, John Wiley & Sons.

60. James M. Lucas, Ronald B. Crosier, Fast Initial Response for CUSUM Quality Control Shames: Give Your CUSUM A Head Start, Technomet-rics, August 1982, vol. 24, No. 3, p. 199-205.

61. James M. Lucas, Counted Data CUSUM's, Technometrics, May 1985, vol. 27, No. 2, p. 129- 144.

62. James M. Lucas and Michael S. Saccuci, Response, Technometrics, February 1990 vol. 32, No., p. 27-28.

63. U.S. DOE Order № 5633.3B. "Control and Accountability of Nuclear Materials". -1994.

64. The Guide on implemtntation of U.S. DOE Order № 5633.3B. "Control and Accountability of Nuclear Materials". -1995.

65. S. Deron, et al., @1993 International Taret Values for Uncertainty Components in Fissile Isotope and Element Accountancy for the Effective Safeguarding of Nuclear Materials", Jornal of Nuclear Materials Management, Vol. XXII, No. II.

66. T. Burr, A Coulter, J. Prommel, "Variance Propagation by Simulation (VPSim)", LA-UR-97-2549,presented at 38th annual meeting of INMM, Phoenix Arizona, July 20-24 1997.

67. M.T. Franklin. Statictical Methods for cmaf Computations. 1996.

68. Andrews D.W.K. Test for Parameter Instability and Structural Change With Unknown Change Point. Econometrica, Vol 6 № 4. Jul, 1993.

69. Brown R.L., Durbin J., Evans J.M. Techniques for Testing the Constancy of Regression Relationship over Tame. Jornal of Royal Statistical Society Series В (Metodological). Vol.37, №2. 1995.

70. A.H. Ширяев Статистический последовательный анализ. М. :Наука, 1969. -275с.

71. А.Н.Горбань, Д.А.Россиев Нейронные сети на персональном компьютере.- Новосибирск.: Наука. Сибирская издательская фирма РАН, 1996.-218с.

72. Бриллинджер Временные ряды. -М.: Наука, 1980. 326с.

73. Хеннан Многомерные временные ряды. -М. :МИР, 1984. 575с.

74. С.А.Айвазян, И.С.Енюков, Л.Д.Мешалкин Прикладная статистика-Исследование зависимостей. -М.: Финансы и статистика, 1985. 312с.

75. Кремер Н.Ш., Путко Б.А. Эконометрика: Учебное пособие для вузов/ Под ред.Н.Ш.Кремера. М.: ЮНИТИ-ДАНА, 2003.- 311с.ной работы