ФЕДЕРАЛЬНОЕ
АГЕНТСТВО Статистические методы ВЫБОРОЧНЫЙ
КОНТРОЛЬ Часть 1 Общие принципы ISO 11648-1:2003 Statistical aspects of sampling from bulk (IDT)
Предисловие Цели и принципы стандартизации в Российской Федерации установлены Федеральным законом от 27 декабря 2002 г. № 184-ФЗ «О техническом регулировании», а правила применения национальных стандартов Российской Федерации - ГОСТ Р 1.0-2004 «Стандартизация в Российской Федерации. Основные положения» Сведения о стандарте 1 ПОДГОТОВЛЕН Автономной некоммерческой организацией «Научно-исследовательский центр контроля и диагностики технических систем» (АНО «НИЦ КД») на основе собственного аутентичного перевода стандарта, указанного в пункте 4 2 ВНЕСЕН Техническим комитетом по стандартизации ТК 125 «Статистические методы в управлении качеством продукции» 3 УТВЕРЖДЕН И ВВЕДЕН В ДЕЙСТВИЕ Приказом Федерального агентства по техническому регулированию и метрологии от 7 декабря 2009 г. № 583-ст 4 Настоящий стандарт идентичен международному стандарту ИСО 11648-1:2003 «Статистические аспекты выборочного контроля нештучной продукции. Часть 1. Общие принципы» (ISO 11648-1:2003 «Statistical aspects of sampling from bulk materials - Part 1: General principles»). Наименование настоящего стандарта изменено относительно наименования указанного международного стандарта для приведения в соответствие с ГОСТ Р 1.5-2004 (подраздел 3.5). При применении настоящего стандарта рекомендуется использовать вместо ссылочных международных стандартов соответствующие им национальные стандарты, сведения о которых приведены в приложении ДА 5 ВВЕДЕН ВПЕРВЫЕ Информация об изменениях к настоящему стандарту публикуется в ежегодно издаваемом информационном указателе «Национальные стандарты», а текст изменений и поправок - в ежемесячно издаваемых информационных указателях «Национальные стандарты». В случае пересмотра (замены) или отмены настоящего стандарта соответствующее уведомление будет опубликовано в ежемесячно издаваемом информационном указателе «Национальные стандарты». Соответствующая информация, уведомление и тексты размещаются также в информационной системе общего пользования - на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет СОДЕРЖАНИЕ Введение Настоящий стандарт устанавливает статистические методы выборочного контроля нештучной продукции. Существуют стандарты методов выборочного контроля конкретных видов нештучной продукции такой, как твердое топливо, железная руда и т.п. Настоящий стандарт устанавливает термины и методы выборочного контроля нештучной продукции, на которые не распространяются упомянутые стандарты. Применяемый в настоящем стандарте международный стандарт разработан техническим комитетом ИСО/ТК 69 «Применение статистических методов». ГОСТ Р ИСО 11648-1-2009 НАЦИОНАЛЬНЫЙ СТАНДАРТ РОССИЙСКОЙ ФЕДЕРАЦИИ Статистические методы ВЫБОРОЧНЫЙ КОНТРОЛЬ НЕШТУЧНОЙ ПРОДУКЦИИ Часть 1 Общие принципы Statistical methods. Sampling from bulk materials. Part 1. General principles Дата введения - 2010-12-01 1 Область определенияНастоящий стандарт устанавливает общие принципы применения статистической обработки данных при выборочном контроле нештучной продукции. Стандарт также дает общие рекомендации и примеры оценки дисперсии, проверки прецизионности и наличия смещения при исследовании среднего характеристики качества продукции. Кроме того, настоящий стандарт применим при проведении статистического анализа серий данных с помощью вариограмм и коррелограмм. Настоящий стандарт устанавливает основные термины и определения для выборочного контроля нештучной продукции. Эти термины обеспечивают лучшее понимание методов выборочного контроля, облегчают выполнение необходимых требований. Примечание - ИСО 11648-2 применим к сыпучим материалам. 2 Нормативные ссылкиВ настоящем стандарте использованы ссылки на следующие стандарты: ИСО 565 Серии испытаний. Металлическая проволочная сетка, перфорированные металлические пластины и листы, изготовленные методом гальванопластики. Номинальные размеры отверстий (ISO 565, Test sieves - Metal wire cloth, perforated metal plate and electroformed sheet - Nominal sizes of openings) ИСО 3534-1 (все части) Статистика. Словарь и условные обозначения (ISO 3534 (all parts), Statistics - Vocabulary and symbols) ИСО 5725 (все части) Точность (правильность и прецизионность) методов и результатов измерений (ISO 5725 (all parts), Accuracy (trueness and precision) of measurement methods and results) 3 Термины, определения, обозначения и сокращения3.1 Термины и определения В настоящем стандарте применены термины и определения по серии стандартов ИСО 3534, а также следующие термины с соответствующими определениями. Примечания 1 Слова «в случае нештучной продукции» в определении термина означают, что данное определение относится только к выборочному контролю нештучной продукции. 2 Более подробная информация приведена в приложении А. 3.1.1 нештучная продукция (bulk material): Вещество или материал, в котором выборочные единицы изначально не различимы на макроскопическом уровне. 3.1.2 выборка (sample): Часть общей совокупности нештучной продукции, состоящая из одной или большего количества отобранных выборочных единиц. 3.1.3 отбор выборки (sampling): Процесс извлечения или составления выборки. 3.1.4 отбор простой случайной выборки (simple random sampling): Отбор, при котором выборка из п выборочных единиц отобрана из совокупности таким образом, что все возможные комбинации п выборочных единиц совокупности имеют равную вероятность быть отобранными в выборку. Примечание - При отборе выборки нештучной продукции, если выборочной единицей является разовая проба, место отбора, выделение и извлечение разовой пробы должны гарантировать, что все выборочные единицы имеют равную вероятность быть отобранными. 3.1.5 слой (stratum): Часть совокупности, рассматриваемая как более однородная по отношению к исследуемым характеристикам качества (слои являются взаимно непересекающимися и полными). Примеры - В нештучной продукции слои на основе времени, массы и пространства связаны с периодами времени (например в 15 мин), изготовленной массой (например 100 т), грузоподъемностью транспортного средства или объемом контейнера. 3.1.6 отбор расслоенной выборки (stratified sampling): Отбор выборки, выполненный таким образом, что выборочные единицы извлекают из различных слоев, и при этом из каждого слоя отбирают хотя бы одну выборочную единицу. Примечание - В некоторых случаях пропорции, соответствующие слоям, устанавливают заранее. Однако в пострасслоенной выборке пропорции, соответствующие слоям, заранее неизвестны. 3.1.7 отбор простой расслоенной случайной выборки (stratified simple random sampling): Отбор, при котором из каждого слоя отбирают простую случайную выборку. 3.1.8 систематический отбор выборки (systematic sampling): Отбор выборки в соответствии со схемой. Примечание 1 - В случае нештучной продукции систематический отбор может быть основан на отборе выборочных единиц через установленные расстояния или установленные интервалы времени. Периодичность отбора может быть задана на основе массы или времени. Если периодичность отбора основана на массе, выборочные единицы или разовые пробы отбирают из объемов равной массы. Если периодичность отбора основана на времени, выборочные единицы или разовые пробы должны быть отобраны из потока или с конвейера в одинаковые интервалы времени. В этом случае масса каждой выборочной единицы или разовой пробы должна быть пропорциональна скорости потока материала в момент отбора образца или разовой пробы. Примечание 2 - Если партию делят на слои, систематический отбор расслоенной выборки может быть выполнен путем отбора разовой пробы из каждого слоя. 3.1.9 выборочная единица (sampling unit): В случае нештучной продукции один из элементов генеральной совокупности, объединяющей всю рассматриваемую нештучную продукцию, каждый из которых с равной вероятностью может быть отобран в выборку. Примечание 1 - В случае нештучной продукции выборочные единицы имеют равную вероятность отбора в выборку. После отбора выборочная единица становится частью выборки. Примечание 2 - При отборе выборки нештучной продукции выборочной единицей является первоначальная разовая проба. 3.1.10 прецизионность (precision): Степень близости независимых результатов наблюдений, полученных в конкретных условиях. Примечание 1 - Прецизионность зависит только от распределения случайных ошибок и погрешностей и не имеет отношения к истинному или установленному значению контролируемой величины. Примечание 2 - В качестве прецизионности обычно используют стандартное отклонение результатов наблюдений. Чем больше стандартное отклонение, тем меньше прецизионность. Примечание 3 - Количественные значения прецизионности зависят от установленных условий. Условия повторяемости и воспроизводимости представляют собой два крайних случая установленных условий. 3.1.11 смещение, систематическая погрешность (bias): Разность между математическим ожиданием результатов наблюдений и истинным, или (в его отсутствие) принятым опорным значением. Примечание 1 - Смещение характеризует систематическую ошибку или погрешность в противоположность случайным ошибке или погрешности. Смещение может иметь одну или несколько составляющих. Большее систематическое отклонение от истинного или принятого опорного значения соответствует большему значению смещения. Примечание 2 - Смещение средства измерений обычно оценивают на основе среднего арифметического погрешностей в показаниях средства измерения, найденных по соответствующему количеству повторных измерений. Погрешность показаний средства измерений - отклонение измеренного значения входной величины от ее истинного значения. 3.1.12 партия (lot): Для нештучной продукции это часть совокупности, для которой должны быть определены установленные характеристики качества. Примечание - В коммерческих операциях предметом сделки часто является единственная партия. В этом случае партия совпадает с генеральной совокупностью. 3.1.13 часть партии (sub-lot): Установленная часть нештучной продукции, относящаяся к партии. 3.1.14 разовая проба (increment): Количество нештучной продукции, отбираемое за одно действие устройством для отбора выборки. Примечание 1 - Место отбора, выделение и извлечение разовой пробы должны гарантировать, что все части нештучной продукции в партии имеют равную вероятность быть отобранными. Примечание 2 - Если отбор выборки или пробы выполняют в несколько этапов, необходимо различать первоначальную разовую пробу, которую отбирают из партии на первом этапе, и вторичную разовую пробу, которую отбирают из первичной разовой пробы на втором этапе и так далее. 3.1.15 сложная проба (composite sample): В случае нештучной продукции проба, объединяющая две или большее количество разовых проб, отобранных из партии. 3.1.16 объединенная проба (gross sample): В случае нештучной продукции проба, объединяющая все разовые пробы, отобранные из части партии или из всей партии с применением процедур стандартного отбора выборки. 3.1.17 исследуемая проба (test sample): Для нештучной продукции проба, подготовленная к испытаниям или анализу, которую целиком или частично используют для исследований. Примечание - Термин «исследуемая проба» может быть заменен на такие термины как «проба для химического анализа», «проба для определения влажности», «проба для определения размеров частиц» и «проба для физических испытаний». 3.1.18 исследуемая порция (test portion): Для нештучной продукции часть исследуемой пробы, единовременно используемая для анализа или испытаний. 3.1.19 многоэтапный отбор выборки (multi-stage sampling): Для нештучной продукции отбор выборки, при котором выборку отбирают поэтапно, а выборочные единицы на каждом этапе выбирают из выборочных единиц, отобранных на предыдущем этапе. 3.1.20 стандартный отбор выборки (routine sampling): Для нештучной продукции отбор выборки для коммерческих целей, выполняемый в соответствии с процедурой, установленной в одном из стандартов и предназначенной для оценки характеристики качества партии. Примечание - Термин «регламентированный отбор» иногда используют как синоним к термину «стандартный отбор выборки». 3.1.21 экспериментальный отбор выборки (experimental sampling): Для нештучной продукции нестандартный отбор выборки, выполняемый в соответствии с установленным планом эксперимента, применяемым для исследования источников отклонений и/или смещения. 3.1.22 отбор сквозной выборки (interpenetrating sampling): Для нештучной продукции повторный отбор проб из партии или части партии, при котором первоначально отобранные разовые пробы для i-й партии или i-й части партии собирают поочередно в разные контейнеры для составления совокупности сложных проб (Аi, Вi, ..., Сi), используемых при исследовании дисперсии, характеризующей разброс между разовыми пробами в партии или части партии. Примечание 1 - Термин «отбор межуровневой выборки» иногда используют как синоним термину «отбор сквозной выборки». Примечание 2 - Наиболее часто в схемах отбора сквозных выборок используют метод двойного отбора выборок с формированием пар сложных проб (Ai, Вi) для каждой i-й партии или части партии. 3.1.23 повторный отбор проб (replicate sampling): Для нештучной продукции отбор выборки, при котором разовые пробы отбирают попарно (одновременно или последовательно) для формирования совокупности из нескольких сложных проб. 3.1.24 двойной отбор проб (duplicate sampling): Для нештучной продукции это повторный отбор выборки, при котором разовые пробы отбирают попарно (одновременно или последовательно) для формирования двух сложных проб. Примечание - Двойной отбор проб - частный случай повторного отбора проб. 3.1.25 мануальный отбор пробы (manual sampling): Для нештучной продукции отбор разовых проб вручную. 3.1.26 автоматизированный отбор пробы (mechanical sampling): Для нештучной продукции отбор разовых проб с применением технических средств. 3.1.27 резка (cut): Для нештучной продукции один ход технического средства при автоматизированном отборе пробы из потока нештучной продукции. 3.1.28 подготовка пробы (sample preparation): Для нештучной продукции совокупность действий для формирования из отобранных проб исследуемой пробы. Пример - Измельчение, смешивание и деление пробы. Примечание - Для сыпучих материалов завершение каждой операции деления пробы является началом следующей стадии подготовки пробы. Таким образом, количество стадий при подготовке пробы равно количеству выполненных действий деления. 3.1.29 измельчение образца (sample reduction): Для нештучной продукции стадия процесса подготовки пробы, при которой размер частиц материала пробы уменьшают путем дробления, размалывания или растирания. 3.1.30 деление пробы (sample division): Для нештучной продукции стадия процесса подготовки пробы, при которой пробу материала нештучной продукции делят на части, одну или более из которых сохраняют. Пример - Нарезание, механическое деление или деление на четыре равные части (квартование). 3.1.31 деление на части с фиксированным процентом массы (fixed ratio division): Для нештучной продукции деление пробы, при котором сохраняемые части отдельных проб представляют собой фиксированную долю массы исходной пробы. 3.1.32 деление на части с фиксированной массой (fixed mass division): Для нештучной продукции деление пробы, при котором сохраняемые части имеют почти одинаковую массу, независимо от различий в массе делимых проб. 3.1.33 сушка пробы (sample drying): Для нештучной продукции стадия процесса подготовки пробы, включающая ее сушку, предназначенную для доведения влажности пробы до уровня, который не будет давать смещение результатов дальнейших исследований или влиять на дальнейшую подготовку пробы. 3.1.34 стандартная подготовка пробы (routine sample preparation): Для нештучной продукции подготовка пробы, выполняемая в соответствии с процедурами, установленными в одном из стандартов и предназначенными для определения выборочного среднего характеристики качества партии. 3.1.35 нестандартная подготовка пробы (поп routine sample preparation): Для нештучной продукции подготовка пробы, выполняемая для эксперимента, предназначенного для определения параметров схемы выборочного контроля. 3.1.36 номинальный надрешетный размер (nominal top size): Для нештучной продукции размер частиц, выраженный в виде апертурного размера испытательного сита (с квадратными отверстиями в соответствии с ИСО 565), на котором остается не более 5 % пробы. 3.1.37 номинальный подрешетный размер (nominal bottom size): Для нештучной продукции размер частиц, выраженный в виде размера отверстия испытательного сита (с квадратными отверстиями в соответствии с ИСО 565), через которое проходит не более 5 % пробы. 3.1.38 изменчивость характеристики качества (quality variation): Для нештучной продукции стандартное отклонение характеристики качества, определяемое на основе оценки дисперсии, полученной либо по сквозной выборке, отобранной из партии или части партии, либо полученной с помощью вариографического анализа отклонений между разовыми пробами, отстоящими друг от друга на разные интервалы времени. 3.1.39 процедура отбора выборки (пробы) (sampling procedure): Для нештучной продукции требования и/или инструкции, устанавливающие порядок отбора разовых проб и формирования пробы или выборки. 3.1.40 процедура подготовки пробы (sampling preparation procedure): Для нештучной продукции требования и/или инструкции, устанавливающие методы и критерии деления пробы. 3.1.41 схема выборочного контроля (sampling scheme): Для нештучной продукции требования к типу выборочного контроля, а также требования к контролируемому объекту, разовым пробам, формированию выборки (пробы) и выполняемым измерениям. Пример - Схема выборочного контроля может устанавливать, например, проведение последовательного или двухступенчатого контроля. В комбинации с требованиями к типу контроля схема может устанавливать количество разовых проб, отбираемых из партии, количество объединенных проб или сквозных выборок, формируемых для партии, количество исследуемых проб, подготавливаемых из сложной пробы и количество измерений, выполняемых на исследуемой пробе. 3.1.42 система для отбора проб (sampling system): Для нештучной продукции механизм и/или механическая установка для отбора разовых проб и подготовки пробы. 3.2 Условные обозначения и сокращения В настоящем стандарте применены следующие обозначения и сокращения. Обозначения приведены в таблице 1, сокращения - в таблице 2. Таблица 1 - Условные обозначения
Таблица 2 - Сокращения, используемые в качестве нижних индексов
4 Цель и применение статистических методов в выборочном контроле нештучной продукцииДля оценки характеристик или свойств нештучной продукции отбирают выборки нескольких типов. Выборки могут быть отобраны из непрерывного потока вещества отдельной партии или последовательности партий. Необходимость в настоящем стандарте вызвана наличием большого количества причин отклонений, связанных с особенностями партии, методами отбора выборки и подготовки сложных проб, погрешностями измерений. Стандарты на выборочный контроль нештучной продукции (уголь, железная руда и сырая нефть) изданы и пересматриваются в соответствующих технических комитетах. Однако в этих стандартах существуют различия в использовании терминов и применении статистических методов. Поэтому одной из целей настоящего стандарта является установление единых терминов и определений, необходимых для применения выборочного контроля нештучной продукции. Второй целью настоящего стандарта является регламентация требований по применению статистических методов. Например, в вышеупомянутых стандартах приведены различные методы исследований для оценки смещения, и пользователи часто не могут принять решение, какой из них лучше. Настоящий стандарт устанавливает альтернативный метод испытаний для оценки смещения. Математическая модель вышеупомянутых методов не может быть применена для большинства автоматизированных систем отбора проб, существующих на сегодняшний день. Нормальные рабочие условия метода испытаний не могут быть точно смоделированы, если система отбора проб не предназначена для работы в нормальных условиях. Предложенный метод испытаний является модификацией обычного метода испытаний для контроля смещения на основе парных данных. Метод испытаний предусматривает прямую оценку дисперсии ошибки посредством двойных измерений каждого элемента пар данных. Это обеспечивает получение большего количества информации о дисперсии ошибок по сравнению с другими методами оценки смещения. Кроме того, рекомендуется включать последовательный анализ данных на основе вариограммы в планы выборочного контроля для нештучной продукции. В настоящем стандарте описано несколько вариантов применения анализа данных различного вида. Основным направлением применения процедур выборочного контроля нештучной продукции является торговая и коммерческая деятельность. Выборочный контроль нештучной продукции выполняют с помощью процедур двух типов. В одном случае выборочный контроль нештучной продукции применяют для оценки среднего значения характеристики качества в партии с приемлемой точностью, в другом случае - для принятия решения о приемке партии. Стандарты для процедур первого типа применимы к выборочному контролю угля, железной руды и иной нештучной продукции. Настоящий стандарт является основополагающим серии стандартов ИСО 11648. Стандарт, устанавливающий процедуры второго типа, - это ИСО 10725 [30]. Выборочный контроль нештучной продукции различается по категории его применения. К одной категории относится выборочный контроль из совокупности нештучной продукции, как описано выше, а к другой - выборочный контроль на предприятии. Цель применения выборочного контроля на предприятии - управление производственным процессом и подтверждение качества продукции. Например, при производстве стали кислородно-конверторным способом выборочный контроль жидкой стали выполняют для управления производственным процессом, а результаты контроля используют для подтверждения того, что химический состав стали соответствует установленным требованиям. Таким образом, применяемые методы выборочного контроля предприятие выбирает самостоятельно, но они должны соответствовать методам выборочного контроля, установленным в стандартах серии ИСО 11648. 5 Особые задачи для выборочного контроля нештучной продукцииЕсли партия состоит из сотен ламп или болтов, случайный выбор ламп или болтов дает представительную выборку из партии. В случае выборочного контроля нештучной продукции вместо отдельных ламп или болтов из партии отбирают разовые пробы. При выборочном контроле нештучной продукции очень важно определить минимальную массу разовой пробы. Например, при последовательном плане выборочного контроля нештучной продукции, упакованной в мешки по 50 килограмм (например, муки или цемента), необходимо: - выбрать мешки для отбора выборки; - определить массу разовой пробы; - отобрать разовые пробы из мешков устройством для отбора пробы и составить представительную выборку (т.е. избежать смещения, соответствующего слоям продукции с различными свойствами); - выполнить необходимую подготовку и исследование пробы. При выборе устройства по отбору проб следует учесть, что слишком маленькое устройство может внести смещение из-за пропуска больших частиц в партии, а слишком большое устройство может привести к чрезмерной нагрузке при подготовке проб. Соответственно, устройство для отбора проб должно быть определено на основе компромисса между этими крайними вариантами. Однако при выборочном контроле сыпучих материалов необходимо рассмотреть влияние окружающей среды и удобство работы с разовыми пробами, в том числе массу разовой пробы, которая при вычислении по приведенной ниже формуле может быть слишком мала для ручного метода отбора проб. На практике применяются ручные (мануальные) и автоматизированные методы отбора проб. При контроле сыпучих материалов минимальная масса разовой пробы при мануальном отборе проб соответствует объему с размерами (3´3´3)d, где d - номинальный максимальный размер частиц в партии, выраженный в миллиметрах. Мануальный отбор разовой пробы основан на предположении о случайном отборе проб из партии. При отборе выборки с остановленной ленты конвейера профилированную рамку для отбора проб с минимальными внутренними размерами, равными утроенному номинально максимальному размеру частиц в партии или 30 мм (какой из них больше), размещают на неподвижном конвейере, вставляя ее в материал в поперечном сечении до упора с лентой конвейера по всей ее ширине. Удаляют материал в пределах рамки для отбора проб, перемещая его с ленты, и ссыпают каждую разовую пробу так, чтобы все ее частицы были включены в разовую пробу. Отбор выборки с остановленной ленты конвейера не всегда удобен, но часто является предпочтительным методом. Минимальная масса разовой пробы, отобранной типовым пробоотборником из отгружаемой с конвейера нештучной продукции, получается при минимальном отверстии пробоотборника и максимальной скорости отбора пробы пробоотборником. Максимальную скорость пробоотборника ограничивают так, чтобы избежать появления смещения, связанного с пропуском больших частиц материала. Масса разовой пробы, отобранной типовым пробоотборником, обычно в 10 - 50 раз больше массы разовой пробы, отобранной вручную. Поперечные ленточные пробоотборники отбирают разовую пробу из материала, движущегося по конвейеру. Пробоотборник должен вырезать поток нештучной продукции в перпендикулярной к поверхности конвейера плоскости. При отборе выборки из зернистого материала математическое ожидание оценки дисперсии выборочного среднего характеристики качества в партии имеет следующий вид где Nite - количество элементов в совокупности; пite - количество элементов в выборке; - оценка дисперсии показателя качества, характеризующей разброс между элементами выборки. В выражении (1) сомножитель (1 - nite/Nite) называют «поправкой на конечность совокупности». Если значение nite/Nite меньше 0,1, то поправкой можно пренебречь. В выборочном контроле нештучной продукции значение, соответствующее nite/Nite, т.е. nsub/Nsub, обычно меньше 0,1, и поправка на конечность совокупности может быть опущена. Здесь nsub - количество разовых проб, отобранных из части партии, а Nsub - общее количество возможных разовых проб в части партии. Этот вывод применим не только к стадии отбора выборки (отбора разовых проб), но также и к стадии подготовки пробы (отбор исследуемой пробы из объединенной пробы), стадии анализа (отбор исследуемой порции из исследуемой пробы). Он также применим к жидкостям и газам. Поправку на конечность совокупности необходимо применять при выборочном контроле нештучной продукции из вагона поезда, металлических бочек из кузова грузовика и т.п. Характеристики качества, которые подлежат контролю, устанавливают в договоре. Влажность вещества в партии определяют путем вычисления сухой массы на основе взвешивания. При этом анализируют различные виды химического состава веществ. При выборе весового оборудования для определения сухой массы компонента необходимо учитывать точность весов, влажность вещества в партии и состав композиции. Иногда определяют распределение частиц по размерам, физические и химические свойства материала. Установленные методы выборочного контроля должны удовлетворять всем требованиям по каждой характеристике качества. 6 Различия для сыпучих веществ, жидкостей и газовПроцесс выборочного контроля сыпучих материалов обычно делится на три стадии: a) процесс отбора разовых проб; b) процесс подготовки пробы; c) процесс измерений. Каждой стадии процесса соответствует своя дисперсия: - дисперсия отбора выборки, характеризующая изменчивость, возникающую при отборе разовых проб; - дисперсия подготовки пробы, характеризующая изменчивость, возникающую при подготовке исследуемой пробы; - дисперсия измерений, характеризующая прецизионность используемого метода измерений (аналитического метода); Если при систематическом выборочном контроле отобрано п разовых проб из партии сыпучего материала, объединенная проба состоит из п разовых проб. Из объединенной пробы подготавливают исследуемую пробу, на которой выполняют пM измерений. Для дисперсии оценки среднего1) характеристики качества в партии можно использовать аппроксимацию (2) где - дисперсия, характеризующая изменчивость разовых проб в пределах слоя, включая каждую разовую пробу партии; ____________ 1) Среднее - математическое ожидание (ред.). - дисперсия, характеризующая изменчивость между исследуемыми пробами, подготовленными для анализа из объединенной пробы; - дисперсия, характеризующая изменчивость результатов измерений, полученных на исследуемой пробе; п - количество разовых проб, отобранных из партии; пМ - количество измерений, выполненных на исследуемой пробе. Примечание - Теория систематического выборочного контроля изложена в [1] и [2]. Если необходимо, чтобы дисперсия была меньше предельного значения, тогда путем подбора соответствующей комбинации количества разовых проб п и количества измерений пМ уменьшают значения первого и третьего членов в выражении (2), поскольку второй член этого выражения сложно изменить. Если дисперсия, характеризующая изменчивость в исследуемой пробе , составляет большую часть выражения (2), а должно быть меньше заданного предельного значения, дальнейшее уменьшение только за счет увеличения п и пм может быть невозможно. В частности, уменьшение дисперсии, соответствующей исследуемой пробе и характеризующей изменчивость, связанную с приготовлением пробы, труднодостижимо при подготовке пробы с зернистостью, близкой к естественной. Единственное решение в этом случае - разбиение партии на соответствующее количество частей партии. Если партия разделена на иlot частей с отбором nsub разовых проб из каждой части партии, объединенную пробу составляют для каждой части партии и выполняют пМ повторных измерений на каждой объединенной пробе, то дисперсия оценки среднего характеристики качества партии будет иметь вид (3) Таким образом, дисперсия оценки среднего характеристики качества партии может быть уменьшена путем выбора соответствующего количества частей партии utot. Часть партии должна представлять собой известное количество нештучной продукции для вычисления оценки характеристики качества партии в виде взвешенного среднего. В процессе отбора выборки из жидкости изменчивость в пределах объединенной пробы сравнительно мала и процесс подготовки пробы обычно отсутствует. При необходимости, объединенная проба может быть перемешана, что позволяет сделать эту изменчивость незначительной. В процессе отбора выборки из газообразного вещества разовую пробу, отобранную из партии, сразу подвергают анализу, и процесс подготовки пробы обычно отсутствует. При отборе выборки из сыпучего материала по возможности весь материал должен быть гомогенизирован, например, путем выделения нескольких партий до отбора разовых проб. Системы отбора проб для сыпучих материалов должны быть подготовлены до начала погрузки материала на суда. Это позволяет уменьшить изменчивость характеристики качества в пределах партии. Отбор разовых проб из слоев, на которые разделена партия, снижает изменчивость характеристики качества на стадии подготовки пробы. Уменьшение размеров частиц является еще одним способом гомогенизации. На стадии исследуемой пробы в лаборатории выполняют механическое перемешивание. Однако специальные операции по гомогенизации на данном этапе могут иногда приводить к разделению материала по свойствам. 7 Экспериментальные методы оценки дисперсии на различных стадиях отбора выборки7.1 Дисперсии, соответствующие различным стадиям отбора выборки Используемые планы контроля нештучной продукции должны быть установлены так, чтобы была достигнута установленная полная прецизионность для партии с учетом прошлого опыта и результатов специальной серии экспериментов. Составляющие дисперсии при стандартном отборе выборки обычно делятся на дисперсию отбора выборки (отбора разовых проб), дисперсию подготовки пробы и дисперсию измерений. Для оценки этих дисперсий, по отдельности или совместно, используют следующие три способа отбора выборки: - иерархический эксперимент; - отбор сквозной выборки; - систематический отбор выборок на основе массы с измерениями каждой разовой пробы поочередно. 7.2 Иерархические эксперименты В начале выборочного контроля, когда нет предыдущего опыта, должен быть выполнен эксперимент для оценки составляющих полной дисперсии на различных стадиях отбора выборки, т.е. дисперсии, характеризующей разброс между партиями; дисперсии, характеризующей разброс между разовыми пробами; дисперсии, характеризующей разброс между пробами, и дисперсии, характеризующей погрешности (ошибки) измерений. Самый простой план эксперимента - это полностью сгруппированный эксперимент с двумя пробами или измерениями на каждой стадии, как показано на рисунке 1. 1 - партия; 2 - стадия отбора выборки; 3 - стадия подготовки пробы; 4 - стадия выполнения измерений Рисунок 1 - Полностью сгруппированный эксперимент Чтобы получить достаточную информацию о составляющих дисперсии на стадиях контроля, необходимо проверить приблизительно 20 партий (хотя в большинстве случаев может потребоваться несколько пар на каждой стадии из одной партии). Недостатком такого плана эксперимента является то, что на каждой стадии необходимо выполнить четыре измерения, а это достаточно трудоемко. Число степеней свободы и математическое ожидание дисперсии для этого случая приведены в таблице 3. Таблица 3 - Таблица ANOVA для полностью сгруппированного эксперимента
Число степеней свободы 4р для дисперсии, соответствующей измерениям, является слишком большим и более предпочтителен план эксперимента, которому соответствует более равномерное распределение числа степеней свободы. Это может быть сделано за счет использования нерегулярного иерархического эксперимента (см. рисунок 2). 1 - партия; 2 - стадия отбора выборки; 3 - стадия подготовки пробы; 4 - стадия выполнения измерений Рисунок 2 - Нерегулярный иерархический эксперимент В случае нерегулярного иерархического эксперимента количество измерений сокращается с 8р до 4р. Число степеней свободы и математическое ожидание дисперсии приведены в таблице 4. В таком эксперименте оценки дисперсий и находят, приравнивая каждую формулу в последней графе таблицы к соответствующей наблюдаемой дисперсии. Детали полностью сгруппированного эксперимента и нерегулярного иерархического эксперимента в соответствии с ИСО 5725-3. Примеры иерархического эксперимента приведены в приложении В. Таблица 4 - Таблица ANOVA для нерегулярного иерархического эксперимента
Отбор сквозной выборки применяют в тех случаях, когда дисперсия отбора выборки существенно больше дисперсии подготовки пробы и дисперсии измерений. Кроме того, отбор сквозной выборки применим в случае, когда допустимо объединение или накопление разовых проб, т.е. к сыпучим материалам или жидкостям. При систематическом выборочном контроле железной руды периодически проводят анализ изменчивости характеристики качества, включая отбор двух разовых проб в пределах слоя. Партию делят более чем на 10 частей, и количество разовых проб, отбираемых из каждой части партии, определяют путем деления общего количества разовых проб, определенных в соответствии с массой партии, на количество частей. Разовые пробы отбирают через установленные объемы массы вещества. Разовые пробы с нечетными и четными номерами, отобранные из каждой части, объединяют в две сложные пробы соответственно. (В следующем примере эти сложные пробы обозначены Ai и Вi соответственно, где i - номер части). Для каждой сложной пробы определяют характеристики качества, а также их изменчивость в пределах слоя на основе двух разовых проб. Применяемые методы иллюстрированы на следующих примерах: Пример 1 (количество разовых проб в сложной пробе) × (количество сложных проб в части) × (количество частей) = 3 × 2 × 10.
Рисунок 3 - Отбор сквозной выборки Пример определения общего содержания железа показан в таблице 5. В этом примере из партии отобрано 60 разовых проб. Разовые пробы № 1, № 3 и № 5 объединены в сложную пробу А1, а разовые пробы № 2, № 4 и № 6 объединены в сложную пробу В1. Таким образом, получены сложные пробы А1, ..., А10 и В1, ..., В10, а затем для каждой сложной пробы определено содержание железа после подготовки каждой пробы отдельно. Размах между аi и bi обозначают Ri. Изменчивость характеристики качества в пределах слоя с двумя отобранными разовыми пробами, включающая дисперсии подготовки пробы и выполнение измерений, оценивают по следующей формуле (среднее арифметическое размаха равно 0,23) (4) где d2 - коэффициент для оценки дисперсии по размаху пар данных, подчиняющихся нормальному распределению, d2 = 1,128; 3 - количество разовых проб в каждой сложной пробе. Другие примеры отбора сквозной выборки приведены в С.7. Таблица 5 - Пример отбора сквозной выборки
7.4 Систематический отбор выборки продукции с последовательным измерением каждой разовой пробы Систематический отбор выборки часто применяют для отбора разовых проб нештучной продукции вместо отбора простой случайной выборки на стационарной установке, поскольку его проще выполнить и механизировать. После отбора разовых проб при систематическом отборе выборки выполняют подготовку исследуемых проб для каждой разовой пробы отдельно и затем измеряют характеристику качества по каждой исследуемой пробе. Данные, полученные таким способом, анализируют с помощью коррелограммы или вариограммы. Вариограмма - график, изображающий дисперсию как функцию интервала между исходными данными. Расстояние между соседними данными называют единичным интервалом, расстояние между каждым вторым значением данных называют двойным интервалом и т.д. Значение дисперсии Vexp(t), соответствующее интервалу в k разовых проб, может быть вычислено по следующей формуле где хi - значение характеристики качества для i-ой разовой пробы (i = 1, 2, ..., n); (п - k) - количество пар разовых проб в интервале, кратном k единичных интервалов; t - значение длины интервала для вычисления вариограммы в единицах времени или массы. Коррелограмма - график, изображающий коэффициент корреляции как функцию интервала между исходными данными. Значение коэффициента корреляции r(t), соответствующее интервалу с k разовыми пробами, может быть вычислено по следующей формуле Вариограмма и коррелограмма для одних и тех же данных зависят одна от другой (см. С.3). В зависимости от ситуации применяют одну из них или обе. Отбор выборки железной руды при разгрузке (например, в порту) обычно выполняют с помощью механического оборудования с количеством разовых проб, соответствующим систематическому отбору выборки, основанному на массе, с учетом массы партии, попадающей на основную ленту конвейера, двигающуюся от судна до места складирования. Выборки из разовых проб для анализа размера частиц обычно просеивают поочередно с помощью специального механического устройства. В последнее время при выборочном контроле железной руды разовые пробы из выборки для определения влажности часто измеряют друг за другом по мере поступления после подготовки, чтобы избежать потери влажности при хранении. Массы разовых проб, отобранных при систематическом отборе выборки на основе времени, пропорциональны интенсивности потока продукции, а соответствующая масса партии неизвестна. Соответственно, характеристика качества разовой пробы, отобранной на основе времени, не может быть мерой интенсивности потока продукции. Таким образом, для последовательно поступающих данных погрешности, связанные с подготовкой пробы и выполнением измерений, определяют стандартным образом. Статистический анализ последовательных данных приведен в приложении С. Метод вариограммы разработан, прежде всего, для оценки составляющих выборочной дисперсии для каждого плана выборочного контроля и нескольких различных периодов отбора выборки, например для периода отбора выборки, увеличенного вдвое. Значение вариограммы для единицы интервала соответствует изменчивости характеристики качества в пределах слоя с двумя соседними разовыми пробами при сквозном отборе выборки. Однако метод вариограммы используют для выявления особенностей последовательных данных, а не для непосредственной оценки составляющих выборочной дисперсии. 8 Выбор плана контроля для получения необходимой прецизионностиЕсли при выборочном контроле сыпучих материалов партию делят на ulot частей партии, nsub разовых проб отбирают из каждой части партии, для каждой части партии составляют объединенную пробу и пМ повторных измерений выполняют на каждой объединенной пробе, то дисперсия оценки среднего характеристики качества партии имеет вид (7) Здесь - необходимая прецизионность, и соответствуют материалу выборки, - соответствует применяемому методу измерений характеристики качества. Однако в большинстве случаев достаточно малая величина по сравнению с и . Следовательно, ulot и nsub должны быть основными параметрами, с помощью которых можно обеспечить необходимую прецизионность. В большинстве случаев при выборочном контроле жидкостей дисперсия, соответствующая исследуемой пробе, подготовленной из объединенной пробы , является малой величиной из-за перемешивания объединенной пробы. Соответственно, дисперсия оценки среднего характеристики качества партии имеет вид (8) где - дисперсия характеристики качества разовых проб в партии; п - количество разовых проб, отобранных из партии. В этом случае для регулировки используют п и nМ. Однако если партия разделена на ulot контейнеров равной массы (части партии), дисперсия оценки среднего характеристики качества партии имеет вид (9) где - дисперсия характеристики качества разовых проб в пределах контейнера. При выборочном контроле газов накопление разовых проб является достаточно трудным делом. Соответственно, дисперсия оценки среднего характеристики качества имеет вид (10) где - дисперсия характеристики качества разовых проб в партии. В этом случае для обеспечения необходимой прецизионности используют n и пМ. 9 Оценка прецизионностиПрецизионность, соответствующую стандартному отбору выборки, подготовке пробы и методу выполнения измерений, необходимо периодически проверять на основе двойного отбора проб. В экспериментах при систематическом отборе выборки должно быть отобрано вдвое больше разовых проб, чем при стандартном отборе выборки за половину интервала времени, или объема массы и составлены две сложные пробы, каждая из которых состоит из п разовых проб соответственно, которые должны быть объединены. Две сложные пробы из каждой партии должны быть подготовлены и на них должны быть выполнены измерения в соответствии со стандартными процедурами. Рекомендуется выполнять эксперименты не меньше чем для двадцати партий одного и того же материала. Независимо от количества двойных данных для обнаружения неконтролируемых точек и оценки полученной прецизионности можно применять контрольные карты размаха в соответствии с ИСО 8258 [29]. Практические рекомендации приведены в приложении D. 10 Проверка смещенияДанные, полученные при стандартном отборе выборки, обычно используют для вычисления стоимости товара. Присутствие в полученных данных смещения искажает результаты расчетов. Наличие смещения очень важно для обеих заинтересованных сторон - и покупателя, и поставщика. Смещение является результатом всех воздействий в системе выборочного контроля от отбора разовых проб до измерения характеристики качества. Смещение возникает при появлении отклонений от критериев, планов и установленного выполнения элементов системы выборочного контроля. Для предупреждения появления смещения все элементы выборочного контроля должны быть проверены путем сопоставления их с критериями плана эксперимента. Например, устройство по отбору разовых проб в системе выборочного контроля сыпучих материалов должно захватывать все поперечное сечение потока материала. В противном случае, например при использовании ложечного пробоотборника, смещение обязательно будет присутствовать, даже если смещение может быть не обнаружено. Детали критериев, плана эксперимента, систем отбора проб для сыпучих материалов приведены в ИСО 11648-2. Смещение - это разность между математическим ожиданием результатов испытаний и истинным значением. Однако на практике истинное значение, как правило, неизвестно. Если применение «несмещенного метода» доступно в качестве «принятого эталонного метода», смещение обычно исследуют на основе сравнения результатов испытаний со значениями, полученными несмещенным методом. При автоматизированном отборе пробы из потока сыпучих материалов примером несмещенного метода может быть метод отбора выборки с остановкой конвейера. Пусть xi1, xi2 - результаты измерений двойных разовых проб, полученных с помощью механического устройства, а уi1, уi2 - результаты измерений двойных разовых проб, полученных с остановкой конвейера соответственно. Разовые пробы одного набора должны быть отобраны как можно ближе друг к другу, k - количество разовых проб (рекомендуется больше двадцати). Вычисляют gi = xi1 - xi2, i =1, 2, ..., k (11) (12) hi = yi1 - yi2, (13) (14) Большее из , помещают в числитель F0 и вычисляют отношение F0 дисперсий ошибок или . (15) Если F0 > Fa/2(v1, v2), то нулевую гипотезу отклоняют, следовательно, нет оснований считать, что две группы данных принадлежат одной совокупности с общей дисперсией. Уровень значимости a обычно устанавливают равным 0,05, число степеней свободы и равными k (v1 = v2 = k). Если F0 < Fa/2(v1, v2), можно считать, что две группы данных имеют общую дисперсию. Границы доверительных интервалов T1(х), T2(х) и T1(у), T2(y) вычисляют следующим образом: (16) (17) где - общее среднее арифметическое хi1 и хi2; - общее среднее арифметическое уi1 и уi2. Вычисляют (18) (19) (20) Если то гипотезу о равенстве двух средних отклоняют. Оценкой смещения является . Если абсолютная величина больше максимально допустимого смещения d, для сокращения или устранения смещения необходимо проанализировать фактическое влияние смещения на оценку качества партии. Существует достаточно много стандартизованных статистических методов для исследования смещения. Однако любой применяемый метод должен быть оценен на пригодность и эффективность. Метод, предложенный в настоящем стандарте, является основным для анализа смещения и случайных ошибок, связанных с особенностями исследуемого материала и примененной методикой выполнения измерений. Более детальные рассмотрения и практические рекомендации приведены в приложении Е. 11 Прецизионность и смещение на стадии выполнения измеренийПрецизионность и смещение на стадии выполнения измерений необходимо анализировать в соответствии с требованиями серии стандартов ИСО 5725 и других серий стандартов ИСО 11648. Приложение А
|
Зольность, % |
||||||||
x111 |
x112 |
x121 |
x122 |
x211 |
x212 |
x221 |
x222 |
|
1 |
9,38 |
9,24 |
9,02 |
8,98 |
9,22 |
9,32 |
8,40 |
8,38 |
2 |
9,76 |
9,80 |
9,92 |
9,92 |
9,36 |
9,36 |
9,72 |
9,54 |
3 |
7,40 |
7,26 |
7,32 |
7,40 |
7,55 |
7,61 |
7,57 |
7,49 |
4 |
8,62 |
8,76 |
8,82 |
8,84 |
9,20 |
9,34 |
10,00 |
10,00 |
5 |
9,16 |
9,18 |
8,72 |
8,68 |
8,89 |
8,75 |
9,51 |
9,47 |
6 |
9,08 |
9,08 |
9,06 |
8,86 |
8,80 |
8,84 |
8,76 |
8,60 |
7 |
8,77 |
8,69 |
8,77 |
8,75 |
9,16 |
8,92 |
9,06 |
8,94 |
8 |
8,62 |
8,68 |
8,80 |
8,42 |
8,78 |
9,02 |
8,62 |
8,94 |
9 |
8,60 |
8,74 |
7,10 |
7,22 |
8,88 |
9,10 |
9,08 |
9,00 |
10 |
6,96 |
7,20 |
7,32 |
7,40 |
8,59 |
8,89 |
7,55 |
7,43 |
11 |
8,44 |
8,26 |
7,92 |
7,70 |
8,65 |
8,45 |
8,37 |
8,15 |
12 |
8,24 |
8,00 |
8,38 |
8,12 |
8,42 |
8,26 |
8,78 |
8,72 |
13 |
7,21 |
7,25 |
6,85 |
7,03 |
7,21 |
7,31 |
7,31 |
7,39 |
14 |
8,84 |
9,00 |
8,96 |
8,90 |
9,24 |
9,16 |
9,20 |
9,38 |
15 |
8,45 |
8,51 |
8,91 |
8,79 |
9,00 |
9,06 |
8,86 |
8,96 |
16 |
9,02 |
9,08 |
9,16 |
9,08 |
8,75 |
8,83 |
8,65 |
8,75 |
17 |
8,71 |
8,77 |
8,75 |
8,75 |
8,98 |
8,96 |
9,00 |
9,18 |
18 |
8,77 |
8,92 |
9,24 |
9,32 |
8,82 |
8,64 |
8,32 |
8,42 |
19 |
7,37 |
7,39 |
7,13 |
7,25 |
7,10 |
6,92 |
6,64 |
6,74 |
20 |
10,12 |
10,02 |
9,96 |
9,94 |
10,72 |
10,78 |
10,30 |
10,30 |
Примечание - Значения х111, х112, ..., х222 соответствуют результатам измерений, показанным на рисунке В.1. |
В.4 Статистический анализ с применением контрольной карты
В.4.1 Контрольная карта
На рисунке В.1 в качестве примера приведена контрольная карта части данных на стадии выполнения измерений. Аналогично могут быть построены контрольные карты на стадии исследуемой пробы и сложной выборки.
Рисунок В.1 - Контрольная карта на стадии выполнения измерений
B.4.2 Интерпретация контрольной карты
Точка на карте средних - среднее арифметическое двух измерений на исследуемой пробе. Точка на карте размахов - размах двух измерений на исследуемой пробе.
На карте размахов точки выходят за допустимую область на стадии выполнения измерений. На карте среднего на стадии выполнения измерений в этом примере 14 точек из 20 лежат вне границы допустимой области. Данные карт размаха и среднего показывают, что прецизионность, представленная на карте размахов, является удовлетворительной и устойчивой и позволяет обнаружить изменчивость по точкам, лежащим за границами допустимой области на карте среднего.
B.4.3 Вычисление дисперсии на каждой стадии
С помощью карты размаха могут быть получены следующие значения (см. таблицу В.2):
- на стадии выполнения измерений;
- на стадии исследуемой пробы;
- на стадии сложной пробы.
Таблица В.2 - Вычисление дисперсий, соответствующих различным стадиям
Результаты |
||||||||||||||||||||||
x111 |
x112 |
x11 |
R1 |
x121 |
x122 |
x12 |
R1 |
x1 |
R2 |
x211 |
x212 |
x21 |
R1 |
x221 |
x222 |
x22 |
R1 |
x2 |
R2 |
x |
R3 |
|
1 |
9,38 |
9,24 |
9,31 |
0,14 |
9,02 |
8,98 |
9 |
0,04 |
9,155 |
0,31 |
9,22 |
9,32 |
9,27 |
0,1 |
8,4 |
8,38 |
8,39 |
0,02 |
8,83 |
0,88 |
8,9925 |
0,325 |
2 |
9,76 |
9,8 |
9,78 |
0,04 |
9,92 |
9,92 |
9,92 |
0 |
9,85 |
0,14 |
9,36 |
9,36 |
9,36 |
0 |
9,72 |
9,54 |
9,63 |
0,18 |
9,495 |
0,27 |
9,6725 |
0,355 |
3 |
7,4 |
7,26 |
7,33 |
0,14 |
7,32 |
7,4 |
7,36 |
0,08 |
7,345 |
0,03 |
7,55 |
7,61 |
7,58 |
0,06 |
7,57 |
7,49 |
7,53 |
0,08 |
7,555 |
0,05 |
7,45 |
0,21 |
4 |
8,62 |
8,76 |
8,69 |
0,14 |
8,82 |
8,84 |
8,83 |
0,02 |
8,76 |
0,14 |
9,2 |
9,34 |
9,27 |
0,14 |
10 |
10 |
10 |
0 |
9,635 |
0,73 |
9,1975 |
0,875 |
5 |
9,16 |
9,18 |
9,17 |
0,02 |
8,72 |
8,68 |
8,7 |
0,04 |
8,935 |
0,47 |
8,89 |
8,75 |
8,82 |
0,14 |
9,51 |
9,47 |
9,49 |
0,04 |
9,155 |
0,67 |
9,045 |
0,22 |
6 |
9,08 |
9,08 |
9,08 |
0 |
9,06 |
8,86 |
8,96 |
0,2 |
9,02 |
0,12 |
8,8 |
8,84 |
8,82 |
0,04 |
8,76 |
8,6 |
8,68 |
0,16 |
8,75 |
0,14 |
8,885 |
0,27 |
7 |
8,77 |
8,69 |
8,73 |
0,08 |
8,77 |
8,75 |
8,76 |
0,02 |
8,745 |
0,03 |
9,16 |
8,92 |
9,04 |
0,24 |
9,06 |
8,94 |
9 |
0,12 |
9,02 |
0,04 |
8,8825 |
0,275 |
8 |
8,62 |
8,68 |
8,65 |
0,06 |
8,8 |
8,42 |
8,61 |
0,38 |
8,63 |
0,04 |
8,78 |
9,02 |
8,9 |
0,24 |
8,62 |
8,94 |
8,78 |
0,32 |
8,84 |
0,12 |
8,735 |
0,21 |
9 |
8,6 |
8,74 |
8,67 |
0,14 |
7,1 |
7,22 |
7,16 |
0,12 |
7,915 |
1,51 |
8,88 |
9,1 |
8,99 |
0,22 |
9,08 |
9 |
9,04 |
0,08 |
9,015 |
0,05 |
8,465 |
1,1 |
10 |
6,96 |
7,2 |
7,08 |
0,24 |
7,32 |
7,4 |
7,36 |
0,08 |
7,22 |
0,28 |
8,59 |
8,89 |
8,74 |
0,3 |
7,55 |
7,43 |
7,49 |
0,12 |
8,115 |
1,25 |
7,6675 |
0,895 |
11 |
8,44 |
8,26 |
8,35 |
0,18 |
7,92 |
7,7 |
7,81 |
0,22 |
8,08 |
0,54 |
8,65 |
8,45 |
8,55 |
0,2 |
8,37 |
8,15 |
8,26 |
0,22 |
8,405 |
0,29 |
8,2425 |
0,325 |
12 |
8,24 |
8 |
8,12 |
0,24 |
8,38 |
8,12 |
8,25 |
0,26 |
8,185 |
0,13 |
8,42 |
8,26 |
8,34 |
0,16 |
8,78 |
8,72 |
8,75 |
0,06 |
8,545 |
0,41 |
8,365 |
0,36 |
13 |
7,21 |
7,25 |
7,23 |
0,04 |
6,85 |
7,03 |
6,94 |
0,18 |
7,085 |
0,29 |
7,21 |
7,31 |
7,26 |
0,1 |
7,31 |
7,39 |
7,35 |
0,08 |
7,305 |
0,09 |
7,195 |
0,22 |
14 |
8,84 |
9 |
8,92 |
0,16 |
8,96 |
8,9 |
8,93 |
0,06 |
8,925 |
0,01 |
9,24 |
9,16 |
9,2 |
0,08 |
9,2 |
9,38 |
9,29 |
0,18 |
9,245 |
0,09 |
9,085 |
0,32 |
15 |
8,45 |
8,51 |
8,48 |
0,06 |
8,91 |
8,79 |
8,85 |
0,12 |
8,665 |
0,37 |
9 |
9,06 |
9,03 |
0,06 |
8,86 |
8,96 |
8,91 |
0,1 |
8,97 |
0,12 |
8,8175 |
0,305 |
16 |
9,02 |
9,08 |
9,05 |
0,06 |
9,16 |
9,08 |
9,12 |
0,08 |
9,085 |
0,07 |
8,75 |
8,83 |
8,79 |
0,08 |
8,65 |
8,75 |
8,7 |
0,1 |
8,745 |
0,09 |
8,915 |
0,34 |
17 |
8,71 |
8,77 |
8,74 |
0,06 |
8,75 |
8,75 |
8,75 |
0 |
8,745 |
0,01 |
8,98 |
8,96 |
8,97 |
0,02 |
9 |
9,18 |
9,09 |
0,18 |
9,03 |
0,12 |
8,8875 |
0,285 |
18 |
8,77 |
8,92 |
8,845 |
0,15 |
9,24 |
9,32 |
9,28 |
0,08 |
9,0625 |
0,435 |
8,82 |
8,64 |
8,73 |
0,18 |
8,32 |
8,42 |
8,37 |
0,1 |
8,55 |
0,36 |
8,80625 |
0,5125 |
19 |
7,37 |
7,39 |
7,38 |
0,02 |
7,13 |
7,25 |
7,19 |
0,12 |
7,285 |
0,19 |
7,1 |
6,92 |
7,01 |
0,18 |
6,64 |
6,74 |
6,69 |
0,1 |
6,85 |
0,32 |
7,0675 |
0,435 |
20 |
10,12 |
10,02 |
10,07 |
0,1 |
9,96 |
9,94 |
9,95 |
0,02 |
10,01 |
0,12 |
10,72 |
10,78 |
10,75 |
0,06 |
10,3 |
10,3 |
10,3 |
0 |
10,525 |
0,45 |
10,2675 |
0,515 |
Сумма |
2,07 |
2,12 |
5,235 |
2,6 |
2,24 |
6,54 |
8,3525 |
|||||||||||||||
Примечание - = (2,07 + 2,12 + 2,6 + 2,24)/80 = 0,112875; = (5,235 + 6,54)/40 = 0,294375; = 8,3525/20 = 0,417625. |
На стадии выполнения измерений
На стадии исследуемой пробы
На стадии сложной выборки
(В.1)
(B.2)
(B.3)
где - оценка дисперсии, соответствующей стадии выполнения измерений;
- оценка дисперсии, характеризующей изменчивость характеристики качества между исследуемыми пробами (дисперсия, соответствующая подготовке пробы);
- оценка дисперсии, характеризующей изменчивость характеристики качества между сложными пробами (дисперсия, соответствующая отбору выборки);
d2 - коэффициент, используемый для оценки стандартного отклонения по размаху пар данных, подчиняющихся нормальному распределению для п = 2, d2 = 1,128.
В.5 Использование дисперсионного анализа (ANOVA)
Данные периодического эксперимента могут быть исследованы с применением ANOVA (дисперсионного анализа). Таблица ANOVA приведена в таблице В.3.
Таблица В.3 - Таблица ANOVA для полностью сгруппированного эксперимента по определению зольности
Сумма квадратов отклонений |
Число степеней свободы |
Средний квадрат отклонений |
Математическое ожидание среднего квадрата отклонений (дисперсия) |
|
Партии |
96,172 |
19 |
5,062 |
|
Сложная проба |
9,373 |
20 |
0,469 |
|
Исследуемая проба |
7,679 |
40 |
0,192 |
|
Выполнение измерений |
0,777 |
80 |
0,010 |
|
Сумма |
114,00 |
159 |
||
Примечание - - дисперсия, характеризующая изменчивость характеристики качества между партиями. |
В таблице В.3 средние квадраты отклонений являются несмещенными оценками соответствующих параметров. Следовательно
Таким образом
Итоговая оценка дисперсии полученная с помощью контрольной карты, совпадает со значением, полученным с применением дисперсионного анализа ANOVA. Однако оценки дисперсии на промежуточных стадиях и полученные этими методами, различны. Одной из возможных причин этого могут быть точки, выходящие за допустимые границы в картах размахов (отклонение от нормального распределения). Другой возможной причиной могут быть ошибки при оценке дисперсии.
Доверительные интервалы для дисперсии могут быть найдены методами Саттервейта [5] с использованием c2-распределения, Андерсона-Банкрофта [6] или Моридичи [7] с использованием F-распределения.
a) Метод Саттервейта
В случае однофакторного эксперимента
(B.4)
где
vA и vE - число степеней свободы для дисперсий VА и VE соответственно.
b) Метод Андерсона-Банкрофта
В случае однофакторного эксперимента
(В.5)
где F0 = VA/VE,
F2 = F(vA, vE; a/2),
F1 = 1/F(vE, vA; a/2),
c) Метод Моридичи
В случае однофакторного эксперимента с односторонним планом, если доверительный интервал для с доверительной вероятностью 100(1 - a) % имеет вид
(B.6)
где vA - число степеней свободы для VA;
Значения bU, bL, F(¥, vA) и 1/F(¥, vA) определяют по таблицам (см. [7]).
В эксперименте с двойными измерениями химического состава х % на 24 партиях сыпучего материала дисперсия, соответствующая отбору выборки и разбросу между партиями равна 10,17. Ниже приведены границы доверительного интервала для этой дисперсии с уровнем доверия 90 %
Нижняя граница |
Верхняя граница |
|
Метод Саттервейта |
6,47 |
18,71 |
Метод Андерсона-Банкрофта |
4,81 |
20,35 |
Метод Моридичи |
6,40 |
17,98 |
Пример определения доверительных интервалов для дисперсий по экспериментальным данным на стадии деления пробы см. в [5].
Статистический анализ последовательных данных
С.1 Общие положения
Настоящее приложение содержит рекомендации по практическому применению статистического анализа последовательных данных. Некоторые из примеров статистического анализа последовательных данных приведены в [8], [9], [10], [11], [12]. Метод вариограмм был разработан, прежде всего, для оценки дисперсии при различных планах выборочного контроля и различных интервалах отбора выборки, например для удвоенного интервала отбора выборки. Статистический анализ последовательных данных для нештучной продукции рассмотрен в [13] - [15].
Основная цель статистического анализа последовательно поступающих данных состоит в том, чтобы определить тип изменения данных. Существует четыре типа данных:
a) случайные;
b) циклические;
c) изменяющиеся в соответствии с некоторой тенденцией;
d) коррелированные.
На практике обычно имеет место комбинация двух или трех перечисленных типов данных. Если вариограмма или коррелограмма горизонтальные, данные являются случайными. Если вариограмма или коррелограмма циклические, данные изменяются также циклически с той же периодичностью. Если вариограмма или коррелограмма представляют собой квадратичную кривую, изменения данных подчиняются некоторой тенденции, или имеет место корреляция. В этом случае удобно использовать коррелограмму, чтобы выявить наличие тренда или автокорреляции, которой соответствуют существенные значения коэффициента корреляции.
Если изменения исходных данных не являются случайными для выбора интервала отбора выборки, используемого при систематическом отборе выборок и для принятия решения о применении систематического отбора выборок или систематического случайного отбора выборок, необходимо провести специальные исследования.
В настоящем приложении даны практические рекомендации, относящиеся к:
a) данным доменного производства;
b) толщине бумаги;
c) содержанию влаги в железной руде;
d) распределению размеров частиц и содержанию влаги в железной руде;
e) тенденциям характеристик в железной руде.
На примере последовательного анализа данных доменного производства и толщины бумаги показаны модели вариограмм и коррелограмм.
С.2 Основные принципы
Вариограмма - график дисперсии, как функции интервала между первоначальными данными. Расстояние между соседними данными называют единичным интервалом, расстояние между каждыми вторыми данными называют двойным интервалом и т.д. (см. рисунок С.1).
Последовательные данные |
xi |
xi+1 |
xi+2 |
xi+3 |
xi+4 |
Рисунок С.1 - Графическое представление данных, используемых для построения вариограммы и коррелограммы
Значения дисперсии Vexp(t) и коэффициента корреляции rexp(t) для интервала t могут быть вычислены в соответствии с формулами (5) и (6) в 7.4.
С.3 Данные доменного производства
В таблицах С.1 и С.2 показаны последовательные данные некоторых из важных показателей доменного процесса. Данные собраны в мае (серия а в таблице С.1) и в октябре месяце 1987 г. (серия b в таблице С.2). Производительностью является масса жидкой стали в тоннах, поступившей из летки при одном выпуске плавки. Летка доменной печи открывается в точно указанное время с интервалом в 3 ч, например, в 700, 1000, 1300, 1600 и т.д. Точное соблюдение этого времени очень важно. Летка закрывается перед началом плавки жидким припоем под давлением воздуха в печи. Производительность (Pmi), содержание кремния (wSi) и серы (wSu) в каждой плавке представлены в таблицах С.1 и С.2.
Таблица С.1 - Данные доменной печи, серия а
Pmi(t) |
wSi(%) |
wSu(%) |
№ |
Pmi(t) |
wSi(%) |
wSu(%) |
|
1 |
1183,7 |
0,13 |
0,037 |
37 |
1065,5 |
0,15 |
0,029 |
2 |
1139,8 |
0,10 |
0,052 |
38 |
1187,2 |
0,18 |
0,030 |
3 |
952,4 |
0,19 |
0,033 |
39 |
957,1 |
0,15 |
0,028 |
4 |
1234,4 |
0,22 |
0,055 |
40 |
713,2 |
0,11 |
0,037 |
5 |
491,8 |
0,45 |
0,013 |
41 |
1107,5 |
0,20 |
0,023 |
6 |
949,1 |
0,25 |
0,030 |
42 |
880,7 |
0,13 |
0,045 |
7 |
881,0 |
0,13 |
0,049 |
43 |
1534,2 |
0,15 |
0,033 |
8 |
1308,5 |
0,14 |
0,053 |
44 |
942,9 |
0,12 |
0,035 |
9 |
795,3 |
0,25 |
0,020 |
45 |
992,7 |
0,21 |
0,031 |
10 |
1133,4 |
0,23 |
0,028 |
46 |
1091,3 |
0,22 |
0,028 |
11 |
467,6 |
0,15 |
0,049 |
47 |
1191,6 |
0,23 |
0,026 |
12 |
1158,8 |
0,22 |
0,026 |
48 |
679,0 |
0,34 |
0,018 |
13 |
1246,2 |
0,21 |
0,032 |
49 |
1611,1 |
0,32 |
0,016 |
14 |
1177,7 |
0,24 |
0,025 |
50 |
937,0 |
0,18 |
0,030 |
15 |
943,5 |
0,19 |
0,035 |
51 |
1233,8 |
0,20 |
0,025 |
16 |
1233,7 |
0,22 |
0,022 |
42 |
834,5 |
0,16 |
0,036 |
17 |
648,9 |
0,15 |
0,045 |
53 |
1164,7 |
0,22 |
0,021 |
18 |
1483,9 |
0,17 |
0,022 |
54 |
856,5 |
0,31 |
0,021 |
19 |
810,5 |
0,22 |
0,033 |
55 |
985,3 |
0,18 |
0,025 |
20 |
1011,2 |
0,25 |
0,017 |
56 |
840,0 |
0,16 |
0,029 |
21 |
757,9 |
0,19 |
0,046 |
57 |
938,7 |
0,14 |
0,028 |
22 |
1101,5 |
0,29 |
0,016 |
58 |
1010,0 |
0,16 |
0,026 |
23 |
672,9 |
0,28 |
0,020 |
59 |
1308,1 |
0,19 |
0,023 |
24 |
1197,5 |
0,21 |
0,026 |
60 |
652,2 |
0,34 |
0,019 |
25 |
1169,3 |
0,29 |
0,018 |
61 |
1109,8 |
0,18 |
0,027 |
26 |
756,5 |
0,15 |
0,040 |
62 |
939,3 |
0,28 |
0,023 |
27 |
806,9 |
0,32 |
0,019 |
63 |
1137,2 |
0,23 |
0,026 |
28 |
906,0 |
0,20 |
0,032 |
64 |
776,6 |
0,24 |
0,018 |
29 |
1179,6 |
0,30 |
0,018 |
65 |
654,1 |
0,18 |
0,021 |
30 |
1048,1 |
0,19 |
0,042 |
66 |
1379,0 |
0,22 |
0,035 |
31 |
1133,4 |
0,22 |
0,022 |
67 |
1275,5 |
0,31 |
0,019 |
32 |
1361,2 |
0,18 |
0,038 |
68 |
1376,2 |
0,25 |
0,028 |
33 |
879,4 |
0,25 |
0,023 |
69 |
1023,3 |
0,31 |
0,018 |
34 |
961,5 |
0,19 |
0,032 |
70 |
827,0 |
0,45 |
0,019 |
35 |
1196,0 |
0,22 |
0,023 |
71 |
1322,0 |
0,26 |
0,025 |
36 |
999,5 |
0,17 |
0,027 |
72 |
913,3 |
0,34 |
0,015 |
73 |
1405,4 |
0,32 |
0,020 |
82 |
1130,6 |
0,20 |
0,027 |
74 |
1522,0 |
0,30 |
0,022 |
83 |
1058,0 |
0,20 |
0,034 |
75 |
926,7 |
0,31 |
0,017 |
84 |
1099,8 |
0,42 |
0,018 |
76 |
952,7 |
0,31 |
0,018 |
85 |
889,0 |
0,31 |
0,025 |
77 |
884,4 |
0,27 |
0,022 |
86 |
946,1 |
0,36 |
0,021 |
78 |
996,3 |
0,12 |
0,044 |
87 |
1118,0 |
0,27 |
0,026 |
79 |
1061,8 |
0,24 |
0,023 |
88 |
928,1 |
0,36 |
0,018 |
80 |
826,0 |
0,40 |
0,018 |
89 |
1136,7 |
0,27 |
0,023 |
81 |
893,3 |
0,40 |
0,018 |
90 |
1605,5 |
0,21 |
0,030 |
Таблица С.2 - Данные доменной печи, серия b
Pmi(t) |
wSi(%) |
wSu(%) |
№ |
Pmi(t) |
wSi(%) |
wSu(%) |
|
1 |
1088,1 |
0,32 |
0,023 |
30 |
1197,0 |
0,19 |
0,021 |
2 |
769,3 |
0,32 |
0,014 |
31 |
1093,1 |
0,20 |
0,025 |
3 |
1665,8 |
0,30 |
0,025 |
32 |
1178,4 |
0,16 |
0,027 |
4 |
1191,8 |
0,34 |
0,016 |
33 |
1098,7 |
0,16 |
0,026 |
5 |
1194,9 |
0,20 |
0,020 |
34 |
769,4 |
0,12 |
0,028 |
6 |
1040,6 |
0,27 |
0,024 |
35 |
1071,7 |
0,19 |
0,031 |
7 |
1542,2 |
0,16 |
0,027 |
36 |
787,6 |
0,22 |
0,023 |
8 |
1012,1 |
0,20 |
0,027 |
37 |
1035,2 |
0,24 |
0,020 |
9 |
1770,4 |
0,19 |
0,027 |
38 |
1017,2 |
0,38 |
0,025 |
10 |
913,4 |
0,32 |
0,014 |
39 |
1087,0 |
0,35 |
0,014 |
11 |
976,6 |
0,16 |
0,036 |
40 |
848,3 |
0,28 |
0,018 |
12 |
1033,9 |
0,30 |
0,017 |
41 |
1324,5 |
0,38 |
0,015 |
13 |
1299,4 |
0,25 |
0,027 |
42 |
991,1 |
0,35 |
0,014 |
14 |
1156,7 |
0,37 |
0,012 |
43 |
1017,5 |
0,28 |
0,018 |
15 |
1260,5 |
0,22 |
0,026 |
44 |
1086,9 |
0,27 |
0,015 |
16 |
940,7 |
0,35 |
0,013 |
45 |
1329,2 |
0,17 |
0,031 |
17 |
1488,2 |
0,32 |
0,017 |
46 |
1041,4 |
0,30 |
0,014 |
18 |
1630,9 |
0,22 |
0,019 |
47 |
1188,3 |
0,18 |
0,027 |
19 |
1127,1 |
0,29 |
0,020 |
48 |
1159,7 |
0,17 |
0,022 |
20 |
1123,2 |
0,14 |
0,035 |
49 |
1423,8 |
0,24 |
0,026 |
21 |
1087,2 |
0,26 |
0,022 |
50 |
1262,6 |
0,27 |
0,017 |
22 |
963,2 |
0,18 |
0,022 |
51 |
937,5 |
0,29 |
0,016 |
23 |
1446,9 |
0,25 |
0,019 |
52 |
846,2 |
0,14 |
0,026 |
24 |
1038,7 |
0,16 |
0,027 |
53 |
1066,8 |
0,21 |
0,020 |
25 |
956,0 |
0,23 |
0,025 |
54 |
991,2 |
0,21 |
0,024 |
26 |
846,7 |
0,21 |
0,022 |
55 |
985,6 |
0,15 |
0,025 |
27 |
732,4 |
0,27 |
0,030 |
56 |
877,0 |
0,13 |
0,028 |
28 |
1007,1 |
0,18 |
0,026 |
57 |
1201,4 |
0,27 |
0,017 |
29 |
1039,6 |
0,23 |
0,019 |
58 |
974,9 |
0,23 |
0,023 |
59 |
1071,0 |
0,28 |
0,020 |
75 |
1053,2 |
0,16 |
0,029 |
60 |
726,2 |
0,21 |
0,027 |
76 |
1218,9 |
0,19 |
0,022 |
61 |
642,6 |
0,20 |
0,021 |
77 |
1449,1 |
0,22 |
0,026 |
62 |
1098,2 |
0,18 |
0,022 |
78 |
888,2 |
0,18 |
0,026 |
63 |
1188,4 |
0,15 |
0,028 |
79 |
1172,1 |
0,12 |
0,035 |
64 |
1016,9 |
0,20 |
0,023 |
80 |
1273,2 |
0,32 |
0,017 |
65 |
911,7 |
0,22 |
0,025 |
81 |
1554,6 |
0,18 |
0,029 |
66 |
1038,1 |
0,32 |
0,015 |
82 |
1287,5 |
0,25 |
0,019 |
67 |
746,2 |
0,21 |
0,026 |
83 |
1526,3 |
0,26 |
0,019 |
68 |
1231,7 |
0,18 |
0,025 |
84 |
1224,6 |
0,30 |
0,013 |
69 |
1261,2 |
0,21 |
0,021 |
85 |
1370,8 |
0,22 |
0,023 |
70 |
1102,5 |
0,28 |
0,022 |
86 |
1128,1 |
0,16 |
0,021 |
71 |
1222,3 |
0,19 |
0,030 |
87 |
1043,3 |
0,21 |
0,024 |
72 |
1118,5 |
0,29 |
0,023 |
88 |
952,4 |
0,34 |
0,020 |
73 |
1051,8 |
0,27 |
0,030 |
89 |
1550,1 |
0,26 |
0,022 |
74 |
824,0 |
0,24 |
0,020 |
90 |
1192,6 |
0,40 |
0,014 |
В таблице С.3 приведены значения вариограммы для общих серий Pmi(a) и Pmi(b), wSi(a), wSi(b), wSu(a) и wSu(b) для интервала t (t пробегает значения от 0 до 18).
Таблица С.3 - Значения вариограммы для данных доменной печи
Pmi(a)10-4 |
Pmi(b)10-4 |
wSi(a)102 |
wSi(b)102 |
wSu(a)104 |
wSu(b)104 |
|
1 |
6,618 |
4,448 |
0,398 |
0,336 |
0,922 |
0,326 |
2 |
4,642 |
3,978 |
0,490 |
0,278 |
0,639 |
0,183 |
3 |
6,590 |
4,922 |
0,537 |
0,410 |
0,875 |
0,327 |
4 |
4,907 |
3,621 |
0,447 |
0,381 |
0,510 |
0,235 |
5 |
5,870 |
4,601 |
0,436 |
0,500 |
0,710 |
0,340 |
6 |
3,816 |
3,430 |
0,466 |
0,477 |
0,717 |
0,315 |
7 |
6,494 |
5,589 |
0,451 |
0,466 |
0,695 |
0,286 |
8 |
4,571 |
4,009 |
0,490 |
0,422 |
0,770 |
0,280 |
9 |
7,150 |
4,037 |
0,520 |
0,512 |
0,611 |
0,332 |
10 |
5,070 |
4,249 |
0,482 |
0,470 |
0,702 |
0,310 |
11 |
6,275 |
4,915 |
0,469 |
0,435 |
0,648 |
0,299 |
12 |
4,469 |
4,195 |
0,512 |
0,456 |
0,810 |
0,276 |
13 |
6,696 |
4,526 |
0,464 |
0,406 |
0,533 |
0,285 |
14 |
5,126 |
3,518 |
0,446 |
0,501 |
0,800 |
0,265 |
15 |
5,897 |
4,664 |
0,497 |
0,404 |
0,646 |
0,287 |
16 |
4,565 |
5,354 |
0,495 |
0,466 |
0,836 |
0,281 |
17 |
4,880 |
5,170 |
0,517 |
0,476 |
0,677 |
0,308 |
18 |
5,180 |
4,577 |
0,523 |
0,511 |
0,879 |
0,282 |
В таблице С.4 приведены значения коррелограммы Pmi(a), Pmi(b) и wSi(a), wSi(b), wSu(a) и wSu(b) (×104) для интервала t (t пробегает значения от 0 до 18) и уровней значимости 5 % и 1 %.
Таблица С.4 - Значения коррелограммы для данных доменной печи
Pmi(a) |
Pmi(b) |
wSi(a) |
wSi(b) |
wSu(a) |
wSu(b) |
|
0 |
1,000 |
1,000 |
1,000 |
1,000 |
1,000 |
1,000 |
1 |
-0,256а |
0,103 |
0,344b |
0,214а |
-0,030 |
-0,142 |
2 |
0,127 |
0,179 |
0,188 |
0,351b |
0,266а |
0,358b |
3 |
-0,228а |
-0,046 |
0,106 |
0,033 |
-0,002 |
0,137 |
4 |
0,091 |
0,238а |
0,264а |
0,098 |
0,397b |
0,185 |
5 |
-0,110 |
0,043 |
0,244а |
-0,183 |
0,152 |
-0,167 |
6 |
0,285b |
0,288b |
0,196 |
-0,117 |
0,153 |
-0,071 |
7 |
-0,206 |
-0,173 |
0,197 |
-0,093 |
0,159 |
0,019 |
8 |
0,154 |
0,149 |
0,129 |
0,020 |
0,029 |
0,043 |
9 |
-0,317b |
0,097 |
0,083 |
-0,178 |
0,242а |
-0,127 |
10 |
0,075 |
0,029 |
0,137 |
-0,083 |
0,132 |
-0,067 |
11 |
-0,181 |
-0,116 |
0,142 |
-0,010 |
0,182 |
-0,061 |
12 |
0,169 |
0,059 |
0,074 |
-0,073 |
-0,025 |
-0,010 |
13 |
-0,237а |
-0,015 |
0,159 |
0,053 |
0,334b |
-0,039 |
14 |
0,060 |
0,207 |
0,203 |
-0,186 |
-0,015 |
0,017 |
15 |
-0,069 |
-0,045 |
0,116 |
0,053 |
0,196 |
-0,052 |
16 |
0,179 |
-0,188 |
0,126 |
-0,109 |
-0,054 |
-0,052 |
17 |
0,092 |
-0,173 |
0,091 |
-0,130 |
0,152 |
-0,145 |
18 |
0,006 |
-0,073 |
0,084 |
-0,202 |
-0,118 |
-0,053 |
а Уровень значимости 5 %. b Уровень значимости 1 %. |
На рисунках С.2, С.4 и С.6 показаны вариограммы для Pmi(a), Pmi(b), wSi(a), wSi(b), wSu(a) и wSu(b) соответственно. На рисунках С.3, С.5 и С.7 показаны коррелограммы для Pmi(a), Pmi(b), wSi(a), wSi(b), wSu(a) и wSu(b) соответственно.
Рисунок С.2 - Вариограммы Ртi(а) и Ртi(b)
Рисунок С.3 - Коррелограммы Ртi(а) и Ртi(b)
Рисунок С.4 - Вариограммы для wSi(a) и wSi(b)
Рисунок С.5 - Коррелограммы для wSi(a) и wSi(b)
Рисунок С.6 - Вариограммы для wSu(a) и wSu(b)
Рисунок С.7 - Коррелограммы для wSu(a) и wSu(b)
На рисунке С.2 показана особая ситуация, когда график похож на лезвие пилы. На рисунке С.6 показана аналогичная ситуация только для малых интервалов. Вариограмму на рисунке С.2 можно разделить следующим образом:
Vexp(t) = Vг = Vc, (С.1)
или Vexp(t) = (-1)t+1Va, (C.2)
где Vг - дисперсия случайных изменений;
Vc - дисперсия циклических изменений;
Va - дисперсия, соответствующая амплитуде циклических изменений.
В этом случае Vexp(t) для серии Pmi может быть аппроксимирована с использованием значений, приведенных на рисунке С.2 следующим уравнением
Vexp(t) = 5,2 + (-1)t+1Va×1,0. (C.3)
Однако в общем случае разделить эту дисперсию на составляющие сложно.
Обычно предполагают, что циклическое перемещение, показанное на рисунке С.2, вызвано смещением времени открытия летки (в часах), независимо от времени предыдущего закрытия летки. Однако дальнейшее исследование этого перемещения и действий, которые необходимо предпринять, должны проводить специалисты по процессам плавки. Этот пример показывает, что статистический анализ последовательных данных дает полезную информацию, которая не может быть получена на основе традиционных методов статистического анализа.
На рисунке С.8 показана зависимость между вариограммой и коррелограммой Ртi(а). На этом рисунке коррелограмма отражена на оси r = +1. Рисунок показывает поведение графика [1 - rexp(t)], аналогичное вариограмме Vexp(t). Между вариограммой и коррелограммой существует следующая зависимость
(C.4)
где - общая дисперсия первоначальных серий (константа).
Что касается серии wSu(a), вариограмма Vexp(t) и график [1 - rexp(t)] демонстрируют аналогичное поведение (см. рисунок С.9).
Рисунок С.8 - Соотношение между Vexp(t) и [1 - rexp(t)] для Pmi(a)
Рисунок С.9 - Соотношение между Vexp(t) и [1 - rexp(t)] для wSu(a)
С.4 Толщина бумаги
Настоящее приложение содержит пример коррелированных данных. Толщину бумаги (мм) измеряют в интервалах постоянной длины (5 м) вдоль листа бумаги, 208 результатов измерений представлены в таблице С.5.
Таблица С.5 - Результаты измерений толщины бумаги
Значения в микронах
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|
0 |
572 |
575 |
568 |
558 |
555 |
564 |
576 |
570 |
570 |
575 |
10 |
587 |
562 |
580 |
573 |
576 |
573 |
585 |
576 |
580 |
576 |
20 |
586 |
571 |
575 |
582 |
570 |
568 |
588 |
560 |
580 |
590 |
30 |
596 |
564 |
605 |
576 |
579 |
579 |
581 |
573 |
570 |
576 |
40 |
573 |
576 |
577 |
575 |
573 |
572 |
562 |
580 |
565 |
572 |
50 |
558 |
569 |
559 |
558 |
556 |
572 |
552 |
560 |
565 |
561 |
60 |
553 |
562 |
555 |
560 |
552 |
561 |
559 |
562 |
550 |
553 |
70 |
556 |
548 |
556 |
536 |
560 |
540 |
548 |
539 |
552 |
543 |
80 |
546 |
541 |
542 |
551 |
538 |
547 |
542 |
541 |
546 |
548 |
90 |
543 |
545 |
547 |
559 |
541 |
543 |
553 |
546 |
550 |
559 |
100 |
548 |
557 |
555 |
559 |
565 |
552 |
555 |
552 |
560 |
562 |
110 |
563 |
563 |
562 |
559 |
581 |
574 |
560 |
583 |
563 |
573 |
120 |
577 |
582 |
573 |
571 |
589 |
588 |
587 |
592 |
582 |
589 |
130 |
589 |
584 |
584 |
599 |
580 |
589 |
575 |
573 |
581 |
576 |
140 |
568 |
570 |
586 |
560 |
570 |
574 |
570 |
556 |
569 |
578 |
150 |
569 |
562 |
576 |
570 |
560 |
572 |
570 |
554 |
568 |
558 |
160 |
562 |
562 |
560 |
547 |
550 |
569 |
547 |
561 |
548 |
556 |
170 |
558 |
549 |
553 |
559 |
564 |
536 |
555 |
558 |
539 |
542 |
180 |
564 |
546 |
557 |
550 |
558 |
550 |
549 |
551 |
541 |
552 |
190 |
542 |
553 |
551 |
558 |
574 |
562 |
561 |
560 |
565 |
569 |
200 |
575 |
568 |
569 |
572 |
565 |
562 |
582 |
570 |
- |
- |
Значения вариограммы и коррелограммы для интервала t = 25 приведены в таблице С.6, а графики вариограммы и коррелограммы - на рисунках С.10 и С.11 соответственно.
Таблица С.6 - Вариограмма и коррелограмма для данных по толщине бумаги
Значения вариограммы |
Значения коррелограммы |
t |
Значения вариограммы |
Значения коррелограммы |
|
0 |
0 |
1000а |
13 |
99,928 |
0,496а |
1 |
62,435 |
0,674а |
14 |
113,487 |
0,429а |
2 |
49,638 |
0,741а |
15 |
117,096 |
0,413а |
3 |
48,324 |
0,749а |
16 |
122,034 |
0,389а |
4 |
48,931 |
0,747а |
17 |
132,296 |
0,337а |
5 |
70,569 |
0,636а |
18 |
137,608 |
0,309а |
6 |
58,223 |
0,701а |
19 |
142,077 |
0,286а |
7 |
64,995 |
0,667а |
20 |
152,202 |
0,233а |
8 |
77,583 |
0,604а |
21 |
154,944 |
0,217а |
9 |
78,638 |
0,600а |
22 |
165,304 |
0,166b |
10 |
81,748 |
0,586а |
23 |
170,103 |
0,144 |
11 |
90,165 |
0,542а |
24 |
183,158 |
0,079 |
12 |
94,781 |
0,521а |
25 |
184,011 |
0,077 |
а Очень высокая корреляция. b Высокая корреляция. |
Рисунок С.10 - Вариограмма для данных по толщине бумаги
Рисунок С.11 - Коррелограмма для данных по толщине бумаги
Коррелограмма показывает существенную корреляцию для t = 21 и существенную корреляцию для t = 22. Так как более 20 последовательных данных коррелированы друг с другом, данные разделены на девять групп по 20 последовательных данных, начиная с № 11 для удобства. Анализ ANOVA показывает существенные различия между группами. Доверительные интервалы с 95 %-ным уровнем доверия для среднего в каждой группе показаны на рисунке С.12. Рекомендуется регулировка и более частый контроль отверстия смесителя.
Рисунок С.12 - Доверительные интервалы для каждой границы
Количество разовых проб, отобранных для определения влажности в партии, обычно составляет в среднем приблизительно 50 на партию. Некоторое время назад процедура предусматривала объединение нескольких разовых проб (например, 5 разовых проб) в пробы и измерение влажности по этим пробам. Эту процедуру применяет большинство сталелитейных заводов.
В последнее время влажность железной руды определяют путем последовательного отбора разовых проб, что позволяет избежать потери влажности при хранении проб. Данные, полученные таким способом, рассмотрены ниже.
Данные для железной руды, поставленной на сталелитейные заводы А и В из одной и той же шахты, группируют в зависимости от размера кусков руды на крупную и мелкую фракции. Значения влажности разовых проб с 12 судов составили для крупной фракции А (от S1 до S3), для крупной фракции В (от S4 до S6), для мелкой фракции А (от S7 до S9) и для мелкой фракции В (от S10 до S12) (см. таблицы С.7 и С.8). Проведен анализ этих данных.
Таблица С.7 - Содержание влаги в руде (%) для S1 - S6
S1 |
S2 |
S3 |
S4 |
S5 |
S6 |
||
1 |
2,230 |
2,573 |
3,123 |
2,743 |
2,744 |
2,895 |
|
2 |
2,109 |
2,532 |
3,040 |
2,719 |
2,763 |
2,966 |
|
3 |
2,263 |
2,394 |
2,870 |
2,912 |
2,881 |
2,964 |
|
4 |
1,993 |
2,679 |
2,914 |
2,672 |
2,988 |
2,769 |
|
5 |
2,091 |
2,567 |
3,118 |
2,656 |
2,756 |
2,833 |
|
6 |
2,012 |
2,174 |
2,871 |
2,582 |
2,813 |
2,945 |
|
7 |
2,213 |
2,421 |
2,934 |
2,780 |
2,795 |
2,852 |
|
8 |
2,049 |
2,350 |
2,515 |
3,090 |
2,948 |
2,826 |
|
9 |
2,080 |
2,381 |
2,705 |
2,861 |
2,943 |
2,728 |
|
10 |
2,111 |
2,557 |
2,755 |
2,801 |
3,022 |
2,665 |
|
11 |
2,004 |
2,465 |
2,975 |
2,559 |
3,012 |
2,569 |
|
12 |
1,993 |
2,373 |
2,921 |
2,567 |
2,74 |
2,808 |
|
13 |
2,110 |
2,417 |
2,918 |
2,680 |
2,723 |
2,775 |
|
14 |
2,211 |
2,449 |
2,730 |
2,788 |
2,329 |
3,020 |
|
15 |
2,116 |
2,489 |
3,631 |
2,694 |
2,803 |
2,850 |
|
16 |
2,018 |
2,528 |
2,798 |
2,481 |
2,602 |
2,679 |
|
17 |
1,970 |
2,394 |
2,841 |
2,434 |
2,635 |
2,776 |
|
18 |
1,732 |
2,297 |
2,812 |
2,435 |
2,603 |
2,824 |
|
19 |
2,515 |
2,306 |
2,639 |
2,445 |
2,551 |
2,644 |
|
20 |
1,941 |
2,447 |
2,495 |
2,732 |
2,572 |
2,555 |
|
21 |
1,769 |
2,379 |
2,571 |
2,674 |
2,507 |
2,600 |
|
22 |
1,963 |
2,316 |
2,646 |
2,738 |
2,575 |
2,622 |
|
23 |
1,607 |
2,242 |
2,447 |
2,739 |
2,531 |
2,442 |
|
24 |
1,829 |
2,310 |
2,484 |
2,573 |
2,488 |
2,353 |
|
25 |
1,891 |
2,378 |
2,528 |
2,435 |
2,612 |
2,274 |
|
26 |
1,691 |
2,331 |
2,874 |
2,760 |
2,602 |
2,273 |
|
27 |
1,969 |
2,342 |
2,524 |
2,549 |
2,625 |
2,396 |
|
28 |
2,240 |
2,368 |
2,504 |
2,547 |
2,756 |
2,282 |
|
29 |
2,000 |
2,394 |
2,603 |
2,522 |
2,723 |
2,333 |
|
30 |
2,000 |
2,438 |
2,395 |
2,815 |
2,572 |
2,715 |
|
31 |
1,857 |
2,395 |
2,109 |
2,363 |
2,504 |
2,564 |
|
32 |
1,874 |
2,330 |
2,261 |
2,347 |
2,315 |
2,516 |
|
33 |
1,836 |
2,392 |
2,412 |
2,393 |
2,460 |
2,597 |
|
34 |
1,776 |
2,569 |
2,426 |
2,804 |
2,425 |
5,729 |
|
35 |
1,858 |
2,174 |
2,398 |
2,584 |
2,354 |
2,554 |
|
36 |
1,906 |
2,386 |
2,365 |
2,501 |
2,463 |
2,497 |
|
37 |
1,919 |
2,180 |
2,669 |
2,367 |
2,590 |
2,723 |
|
38 |
2,018 |
2,276 |
2,697 |
2,217 |
2,212 |
2,762 |
|
39 |
1,749 |
2,179 |
2,711 |
2,668 |
2,304 |
2,278 |
|
40 |
1,982 |
2,398 |
2,591 |
2,485 |
2,271 |
2,432 |
|
41 |
1,986 |
2,339 |
2,734 |
2,306 |
2,454 |
2,527 |
|
42 |
1,693 |
2,507 |
2,351 |
2,326 |
2,368 |
2,580 |
|
43 |
1,643 |
2,281 |
2,551 |
2,076 |
2,117 |
- |
|
44 |
1,889 |
2,144 |
2,297 |
2,136 |
2,165 |
- |
|
45 |
1,898 |
1,999 |
2,269 |
2,701 |
2,157 |
- |
|
46 |
2,003 |
1,999 |
2,510 |
- |
1,979 |
- |
|
47 |
1,891 |
2,190 |
2,409 |
- |
- |
- |
|
48 |
1,849 |
2,201 |
- |
- |
- |
- |
|
49 |
1,928 |
2,265 |
- |
- |
- |
- |
|
50 |
1,961 |
2,203 |
- |
- |
- |
- |
|
51 |
1,760 |
1,984 |
- |
- |
- |
- |
|
52 |
1,734 |
2,239 |
- |
- |
- |
- |
|
53 |
1,511 |
2,429 |
- |
- |
- |
- |
|
54 |
1,822 |
2,280 |
- |
- |
- |
- |
|
55 |
1,701 |
1,977 |
- |
- |
- |
- |
|
56 |
1,592 |
2,141 |
- |
- |
- |
- |
|
57 |
1,482 |
2,012 |
- |
- |
- |
- |
|
58 |
1,505 |
- |
- |
- |
- |
- |
|
59 |
1,493 |
- |
- |
- |
- |
- |
|
60 |
1,540 |
- |
- |
- |
- |
- |
|
Таблица С.8 - Содержание влаги в руде (%) для S7 - S12
S7 |
S8 |
S9 |
S10 |
S11 |
S12 |
|
1 |
4,625 |
5,554 |
4,535 |
4,595 |
4,710 |
4,992 |
2 |
4,617 |
5,303 |
4,505 |
4,579 |
4,660 |
5,066 |
3 |
4,558 |
5,267 |
4,605 |
4,559 |
4,696 |
4,911 |
4 |
4,643 |
5,261 |
4,462 |
4,544 |
4,681 |
4,998 |
5 |
4,623 |
5,144 |
4,530 |
4,585 |
4,656 |
5,024 |
6 |
4,656 |
5,197 |
4,482 |
4,574 |
4,684 |
5,061 |
7 |
4,586 |
5,193 |
4,603 |
4,572 |
4,679 |
4,995 |
8 |
4,638 |
5,171 |
4,626 |
4,633 |
4,689 |
4,976 |
9 |
4,474 |
5,159 |
4,503 |
4,616 |
4,633 |
4,982 |
10 |
4,687 |
5,161 |
4,426 |
4,693 |
4,685 |
4,977 |
11 |
4,600 |
5,000 |
4,432 |
4,483 |
4,758 |
4,976 |
12 |
4,659 |
5,292 |
4,437 |
4,559 |
4,691 |
5,042 |
13 |
4,659 |
5,000 |
4,507 |
4,528 |
4,709 |
1,912 |
14 |
4,553 |
5,127 |
4,542 |
4,543 |
4,690 |
4,893 |
15 |
4,591 |
5,108 |
4,569 |
4,571 |
4,649 |
4,982 |
16 |
4,630 |
4,828 |
4,564 |
4,500 |
4,668 |
4,883 |
19 |
4,590 |
4,951 |
4,535 |
4,530 |
4,722 |
4,936 |
20 |
4,617 |
5,019 |
4,505 |
4,488 |
4,594 |
4,956 |
21 |
4,603 |
5,148 |
4,564 |
4,404 |
4,605 |
4,880 |
22 |
4,576 |
5,123 |
4,573 |
4,610 |
4,669 |
4,899 |
23 |
4,585 |
5,079 |
4,564 |
4,474 |
4,615 |
4,992 |
24 |
4,603 |
4,991 |
4,456 |
4,456 |
4,659 |
5,115 |
25 |
4,635 |
5,095 |
4,480 |
4,418 |
4,601 |
4,959 |
26 |
5,285 |
4,972 |
4,503 |
4,427 |
4,572 |
4,997 |
27 |
5,026 |
5,017 |
4,385 |
4,486 |
4,537 |
4,971 |
28 |
4,929 |
5,005 |
4,454 |
4,467 |
4,588 |
4,950 |
29 |
5,000 |
5,050 |
4,478 |
4,464 |
4,564 |
4,960 |
30 |
4,000 |
4,979 |
4,411 |
4,450 |
4,566 |
4,997 |
31 |
4,868 |
5,029 |
4,434 |
4,634 |
4,577 |
4,988 |
32 |
4,733 |
4,977 |
4,521 |
4,535 |
4,637 |
4,963 |
33 |
4,597 |
5,011 |
4,478 |
4,528 |
4,458 |
4,996 |
34 |
4,733 |
5,030 |
- |
4,450 |
4,714 |
4,945 |
35 |
4,496 |
5,019 |
- |
4,458 |
4,648 |
5,041 |
36 |
4,584 |
5,025 |
- |
4,424 |
4,644 |
4,981 |
37 |
4,671 |
5,113 |
- |
4,571 |
4,698 |
5,029 |
38 |
4,723 |
4,849 |
- |
4,456 |
4,689 |
4,885 |
39 |
4,688 |
5,069 |
- |
4,501 |
4,638 |
4,975 |
40 |
4,570 |
- |
- |
4,544 |
4,627 |
4,866 |
41 |
4,522 |
- |
- |
4,390 |
4,659 |
4,659 |
42 |
4,332 |
- |
- |
4,466 |
4,599 |
4,788 |
43 |
4,584 |
- |
- |
4,467 |
4,637 |
4,891 |
44 |
4,557 |
- |
- |
4,391 |
4,677 |
4,778 |
45 |
4,607 |
- |
- |
4,405 |
4,647 |
4,969 |
46 |
4,537 |
- |
- |
- |
4,649 |
4,843 |
47 |
4,534 |
- |
- |
- |
4,515 |
4,907 |
48 |
4,459 |
- |
- |
- |
4,579 |
4,942 |
49 |
4,571 |
- |
- |
- |
4,553 |
4,945 |
50 |
4,684 |
- |
- |
- |
4,561 |
4,915 |
51 |
4,548 |
- |
- |
- |
- |
4,913 |
52 |
4,568 |
- |
- |
- |
- |
4,937 |
53 |
4,368 |
- |
- |
- |
- |
4,947 |
54 |
4,524 |
- |
- |
- |
- |
4,839 |
55 |
4,320 |
- |
- |
- |
- |
4,950 |
56 |
4,479 |
- |
- |
- |
- |
4,818 |
57 |
4,643 |
- |
- |
- |
- |
В качестве примера ниже показан результат статистического анализа серии № 1. Данные серии № 1 показаны на рисунке С.13. На рисунке показана тенденция уменьшения значений влажности с увеличением номера разовой пробы. Регрессионный анализ этих данных показывает существенную зависимость влажности от номера разовой пробы. Уравнение линии регрессии имеет вид
wm = 2,164 - 0,0085i, (С.5)
где i - номер разовой пробы;
wm - доля массы, представляющая собой процент содержания влаги в разовых пробах.
Рисунок С.13 - Результаты измерений влажности в серии № 1
Поскольку количество разовых проб в партии меняется от судна к судну, переход к относительным номерам разовых проб дает следующее уравнение регрессии
wm = 2,164 - 0,51u, (С.6)
где и - относительный номер разовой пробы.
Если общее количество разовых проб, отобранных из партии n, относительный номер для i-й разовой пробы равен i/n. Уравнение (С.6) можно представить в следующей общей форме
wm = b0 + b1u. (C.7)
Значения b1 для каждой серии приведены в таблице С.9.
Таблица С.9 - Результаты статистического анализа влажности железной руды
Сталелитейные заводы |
№ серии данных |
Номера разовых проб |
Угловой коэффициент b1 линии регрессии |
Значение вариограммы Vexp (t = 1) |
|
Крупные фракции |
А |
1 |
60 |
-0,51а |
0,012 |
А |
2 |
57 |
-0,36а |
0,012 |
|
А |
3 |
47 |
-0,37а |
0,016 |
|
В |
4 |
45 |
-0,27а |
0,022 |
|
В |
5 |
46 |
-0,46а |
0,012 |
|
В |
6 |
42 |
-0,47а |
0,011 |
|
С |
13 |
48 |
-0,09а |
0,134 |
|
С |
14 |
54 |
0 |
0,110 |
|
С |
15 |
50 |
0 |
0,166 |
|
Мелкие фракции |
А |
7 |
57 |
0 |
0,023 |
А |
8 |
40 |
0 |
0,007 |
|
А |
9 |
33 |
-0,07b |
0,002 |
|
В |
10 |
45 |
-0,09а |
0,003 |
|
В |
11 |
50 |
-0,10а |
0,002 |
|
В |
12 |
56 |
-0,11а |
0,004 |
|
С |
16 |
33 |
0 |
0,075 |
|
С |
17 |
48 |
0 |
0,017 |
|
С |
18 |
47 |
0 |
0,206 |
|
а Очень высокая корреляция. b Высокая корреляция. |
Вариограмма S1 показана на рисунке С.14. На этом рисунке «первоначальной серии» соответствует вариограмма, обозначенная буквой «b», буквой «а» обозначена вариограммы для серии данных № 19 без выбросов. Следует заметить, что вариограмма с выбросами выше вариограммы без выбросов.
а - без выброса; b - исходная серия
Рисунок С.14 - Вариограмма для серии № 1
Значения вариограмм для S1 - S6 и S7 - S12 указаны в таблицах С.10 и С.11. Графики вариограмм показаны на рисунках С.15, С.16 соответственно. Так как серии фракции S1 - S6 и S7 - S12 представляют собой одну и ту же железную руду соответственно, важно, что вариограммы на этих рисунках в большой степени согласованы.
Таблица С.10 - Значения вариограммы для S1 - S6
S1 |
S1Ra |
S2 |
S3 |
S3Ra |
S4 |
S5 |
S6 |
S6Ra |
|
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
1 |
0,020 |
0,012 |
0,012 |
0,032 |
0,016 |
0,022 |
0,012 |
0,254 |
0,011 |
2 |
0,022 |
0,014 |
0,015 |
0,033 |
0,022 |
0,032 |
0,015 |
0,278 |
0,019 |
3 |
0,021 |
0,017 |
0,017 |
0,039 |
0,026 |
0,034 |
0,020 |
0,262 |
0,018 |
4 |
0,026 |
0,018 |
0,016 |
0,043 |
0,027 |
0,031 |
0,022 |
0,254 |
0,020 |
5 |
0,024 |
0,021 |
0,015 |
0,056 |
0,033 |
0,024 |
0,023 |
0,339 |
0,025 |
6 |
0,025 |
0,022 |
0,011 |
0,059 |
0,038 |
0,020 |
0,024 |
0,342 |
0,026 |
7 |
0,032 |
0,024 |
0,015 |
0,064 |
0,040 |
0,031 |
0,025 |
0,334 |
0,031 |
8 |
0,026 |
0,023 |
0,016 |
0,057 |
0,035 |
0,035 |
0,028 |
0,352 |
0,037 |
9 |
0,025 |
0,026 |
0,015 |
0,058 |
0,034 |
0,037 |
0,035 |
0,219 |
0,037 |
10 |
0,030 |
0,025 |
0,014 |
0,050 |
0,032 |
0,036 |
0,038 |
0,218 |
0,040 |
11 |
0,032 |
0,027 |
0,014 |
0,041 |
0,028 |
0,031 |
0,036 |
0,213 |
0,042 |
12 |
0,032 |
0,026 |
0,014 |
0,056 |
0,031 |
0,027 |
0,039 |
0,200 |
0,040 |
13 |
0,033 |
0,026 |
0,017 |
0,058 |
0,039 |
0,030 |
0,043 |
0,210 |
0,041 |
14 |
0,033 |
0,027 |
0,019 |
0,060 |
0,043 |
0,036 |
0,044 |
0,224 |
0,043 |
15 |
0,039 |
0,034 |
0,021 |
0,066 |
0,049 |
0,037 |
0,051 |
0,217 |
0,042 |
16 |
0,033 |
0,028 |
0,018 |
0,084 |
0,059 |
0,039 |
0,052 |
0,202 |
0,044 |
17 |
0,036 |
0,032 |
0,021 |
0,083 |
0,067 |
0,043 |
0,054 |
0,224 |
0,052 |
18 |
0,035 |
0,031 |
0,019 |
0,095 |
0,072 |
0,042 |
0,051 |
0,248 |
0,054 |
а R после обозначения серии указывает на обработку с исключением выбросов. |
Таблица C.11 - Значения вариограммы для S7 - S12
S7 |
S7Ra |
S8 |
S8Ra |
S9 |
S10 |
S11 |
S12 |
|
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
1 |
0,025 |
0,023 |
0,008 |
0,007 |
0,002 |
0,003 |
0,002 |
0,004 |
2 |
0,025 |
0,021 |
0,006 |
0,005 |
0,003 |
0,003 |
0,002 |
0,004 |
3 |
0,028 |
0,025 |
0,006 |
0,005 |
0,004 |
0,003 |
0,002 |
0,004 |
4 |
0,039 |
0,023 |
0,010 |
0,008 |
0,004 |
0,004 |
0,003 |
0,005 |
5 |
0,025 |
0,019 |
0,008 |
0,007 |
0,004 |
0,004 |
0,003 |
0,005 |
6 |
0,028 |
0,021 |
0,009 |
0,008 |
0,003 |
0,004 |
0,003 |
0,005 |
7 |
0,031 |
0,022 |
0,011 |
0,009 |
0,003 |
0,004 |
0,003 |
0,006 |
8 |
0,033 |
0,027 |
0,010 |
0,008 |
0,003 |
0,004 |
0,003 |
0,006 |
9 |
0,036 |
0,026 |
0,012 |
0,010 |
0,004 |
0,003 |
0,004 |
0,007 |
10 |
0,034 |
0,024 |
0,014 |
0,009 |
0,004 |
0,004 |
0,004 |
0,006 |
11 |
0,032 |
0,024 |
0,011 |
0,010 |
0,004 |
0,005 |
0,004 |
0,007 |
12 |
0,031 |
0,021 |
0,017 |
0,012 |
0,004 |
0,005 |
0,004 |
0,006 |
13 |
0,036 |
0,029 |
0,014 |
0,011 |
0,003 |
0,005 |
0,004 |
0,005 |
14 |
0,040 |
0,029 |
0,018 |
0,015 |
0,003 |
0,005 |
0,004 |
0,007 |
15 |
0,044 |
0,032 |
0,020 |
0,009 |
0,003 |
0,005 |
0,004 |
0,006 |
16 |
0,041 |
0,030 |
0,017 |
0,012 |
0,003 |
0,006 |
0,004 |
0,006 |
17 |
0,045 |
0,032 |
0,017 |
0,012 |
0,003 |
0,005 |
0,004 |
0,007 |
18 |
0,041 |
0,028 |
0,019 |
0,012 |
0,003 |
0,005 |
0,003 |
0,006 |
a R после обозначения серии указывает на обработку с исключением выбросов. |
Рисунок С.15 - Вариограмма для S1 - S6
Рисунок С.16 - Вариограмма для S7 - S12
Результаты статистического анализа каждой серии приведены в таблице С.9. Из таблицы С.9 видно, что для большинства серий характерны существенные тренды, и линии регрессии всегда имеют отрицательный наклон. Угловой коэффициент линии регрессии для крупных фракций по абсолютной величине больше, чем для мелких фракций, тогда как уровни содержания влаги крупных фракций ниже, чем мелких. Полезно рассмотреть часть партии, соответствующую двум смежным разовым пробам как слой. Изменчивость характеристики качества в пределах слоя можно оценить по значениям вариограммы для t = 1, как показано в таблице С.9.
Дополнительные данные представлены сталелитейным заводом С, который расположен в области усиленного контроля требований запыленности. Импортированные железные руды распыляются обычно при разгрузке. Значения влажности разовых проб, взятых с шести судов, из которых три судна с крупной фракцией (от S13 до S15) и три судна с мелкой фракцией (от S16 до S18), приведены в таблице С.12. Регрессионный анализ данных с исключением выбросов, при необходимости, не дает результатов, характерных для данных заводов А и В. Значения вариограммы приведены в таблице С.13. Значения вариограммы для интервала t = 1 более чем в 10 раз больше, чем таковые же для руды, распыленной без воды. Вариограммы для S13 - S18 приведены на рисунке С.17. Результаты статистического анализа на данных, полученных с завода С, также приведены в таблице С.9.
Таблица С.12 - Содержание влаги в руде (%) для S13 - S18
S13 |
S14 |
S15 |
S16 |
S17 |
S18 |
|
1 |
3,299 |
3,432 |
3,532 |
5,293 |
5,080 |
6,442 |
2 |
3,451 |
2,585 |
2,885 |
5,363 |
4,897 |
6,184 |
3 |
2,833 |
3,226 |
3,270 |
5,406 |
5,040 |
5,157 |
4 |
2,708 |
3,049 |
3,051 |
5,368 |
4,949 |
5,036 |
5 |
3,646 |
2,880 |
2,404 |
5,417 |
4,912 |
5,435 |
6 |
3,231 |
2,964 |
2,592 |
5,837 |
5,125 |
4,995 |
7 |
2,205 |
2,424 |
2,603 |
5,083 |
4,986 |
5,008 |
8 |
2,763 |
2,403 |
2,320 |
5,002 |
5,074 |
5,230 |
9 |
3,294 |
2,331 |
2,413 |
5,415 |
5,032 |
6,290 |
10 |
2,423 |
2,402 |
2,570 |
5,467 |
5,016 |
4,942 |
11 |
2,671 |
2,602 |
2,421 |
5,126 |
5,054 |
4,716 |
12 |
3,506 |
2,517 |
2,792 |
5,262 |
5,024 |
5,014 |
13 |
3,417 |
2,757 |
2,658 |
5,038 |
4,894 |
5,188 |
14 |
4,022 |
2,453 |
2,383 |
5,326 |
4,938 |
5,064 |
15 |
2,239 |
2,669 |
2,375 |
5,091 |
5,008 |
4,968 |
16 |
2,070 |
2,598 |
2,566 |
5,005 |
5,023 |
5,006 |
17 |
2,571 |
2,728 |
2,886 |
4,938 |
4,949 |
4,980 |
18 |
2,520 |
2,440 |
2,919 |
4,942 |
4,923 |
5,057 |
19 |
2,432 |
2,862 |
3,565 |
4,949 |
4,993 |
4,646 |
20 |
2,707 |
3,018 |
3,484 |
4,894 |
5,235 |
5,205 |
21 |
2,221 |
2,686 |
3,765 |
4,949 |
5,132 |
5,038 |
22 |
2,166 |
2,777 |
2,500 |
5,008 |
5,099 |
5,031 |
23 |
2,710 |
2,339 |
4,177 |
4,963 |
4,995 |
5,007 |
24 |
2,777 |
2,294 |
3,139 |
5,333 |
5,197 |
5,407 |
25 |
2,306 |
2,401 |
3,800 |
4,928 |
5,170 |
5,868 |
26 |
2,344 |
2,882 |
3,840 |
5,253 |
4,736 |
5,313 |
27 |
2,112 |
2,762 |
2,411 |
5,747 |
4,978 |
5,093 |
28 |
1,735 |
2,654 |
2,452 |
5,064 |
5,078 |
5,723 |
29 |
1,772 |
2,474 |
3,273 |
5,205 |
5,063 |
6,378 |
30 |
2,004 |
2,640 |
3,588 |
5,121 |
5,018 |
5,163 |
31 |
2,821 |
2,562 |
2,642 |
5,181 |
5,055 |
6,000 |
32 |
2,780 |
2,926 |
3,481 |
6,296 |
4,889 |
4,688 |
33 |
2,889 |
3,177 |
3,578 |
5,239 |
4,997 |
4,900 |
34 |
2,746 |
3,028 |
2,898 |
- |
5,152 |
5,403 |
35 |
2,725 |
2,519 |
3,444 |
- |
4,942 |
5,416 |
36 |
3,493 |
3,405 |
2,906 |
- |
5,270 |
5,430 |
37 |
2,813 |
2,462 |
2,550 |
- |
3,949 |
6,753 |
38 |
2,315 |
2,483 |
2,749 |
- |
4,983 |
5,932 |
39 |
2,614 |
4,468 |
2,435 |
- |
5,377 |
5,426 |
40 |
2,163 |
3,370 |
2,389 |
- |
5,528 |
5,134 |
41 |
2,151 |
2,726 |
2,703 |
- |
5,515 |
6,154 |
42 |
2,326 |
2,814 |
2,963 |
- |
5,027 |
5,351 |
43 |
2,626 |
2,847 |
2,542 |
- |
4,907 |
6,138 |
44 |
2,438 |
2,261 |
3,455 |
- |
4,946 |
5,093 |
45 |
2,452 |
2,516 |
3,195 |
- |
5,015 |
6,058 |
46 |
2,955 |
2,842 |
3,331 |
- |
5,257 |
5,094 |
47 |
2,381 |
2,788 |
3,959 |
- |
4,820 |
5,086 |
48 |
2,380 |
2,721 |
3,038 |
- |
5,014 |
- |
49 |
- |
3,017 |
3,020 |
- |
- |
- |
50 |
- |
2,506 |
3,251 |
- |
- |
- |
51 |
- |
2,989 |
- |
- |
- |
- |
52 |
- |
2,868 |
- |
- |
- |
- |
53 |
- |
3,069 |
- |
- |
- |
- |
54 |
- |
2,753 |
- |
- |
- |
- |
Таблица С.13 - Значение вариограммы для S13 - S18
Интервал t |
S13 |
S14 |
S15 |
S16 |
S17 |
S18 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
1 |
0,134 |
0,110 |
0,166 |
0,075 |
0,017 |
0,206 |
2 |
0,218 |
0,131 |
0,169 |
0,067 |
0,024 |
0,223 |
3 |
0,192 |
0,107 |
0,169 |
0,055 |
0,026 |
0,246 |
4 |
0,222 |
0,143 |
0,232 |
0,069 |
0,025 |
0,202 |
5 |
0,248 |
0,158 |
0,250 |
0,049 |
0,019 |
0,274 |
6 |
0,252 |
0,123 |
0,240 |
0,071 |
0,026 |
0,222 |
7 |
0,247 |
0,122 |
0,283 |
0,088 |
0,029 |
0,219 |
8 |
0,239 |
0,164 |
0,328 |
0,067 |
0,028 |
0,184 |
9 |
0,230 |
0,162 |
0,256 |
0,090 |
0,021 |
0,272 |
10 |
0,235 |
0,155 |
0,307 |
0,094 |
0,019 |
0,302 |
11 |
0,222 |
0,178 |
0,306 |
0,107 |
0,020 |
0,250 |
12 |
0,192 |
0,148 |
0,263 |
0,107 |
0,017 |
0,199 |
13 |
0,232 |
0,128 |
0,315 |
0,108 |
0,023 |
0,255 |
14 |
0,278 |
0,143 |
0,318 |
0,119 |
0,023 |
0,253 |
15 |
0,298 |
0,166 |
0,331 |
0,118 |
0,021 |
0,288 |
16 |
0,330 |
0,126 |
0,345 |
0,127 |
0,016 |
0,238 |
17 |
0,302 |
0,119 |
0,341 |
0,106 |
0,019 |
0,290 |
18 |
0,252 |
0,141 |
0,347 |
0,084 |
0,022 |
0,362 |
Рисунок С.17 - Вариограмма для S13 - S18
С.6 Распределение по размерам и содержание влаги в железной руде
В С.5 показано, что уменьшение влажности с увеличением номера разовой пробы железной руды не соблюдается в большинстве случаев, особенно большим является уменьшение для руды в крупной фракции. Чтобы узнать причину этих трендов, приведены парные данные влажности и распределения по размерам руды с шести судов. Разовые пробы отобраны в половине стандартного интервала систематическим отбором выборок на основе массы. Разовые пробы с нечетными номерами используют для определения влажности, а четные разовые пробы используют для определения распределения по размеру частиц. Полученные данные приведены в таблице С.14, где wm - влажность в виде процента массы, a wsf - распределение по размерам частиц (фракция железной руды с размерами частиц менее 6 мм).
Таблица С.14 - Содержание влаги (%) и распределение частиц по размерам (% фракции < 6 мм) для S19 - S24
S19 |
S20 |
S21 |
S22 |
S23 |
S24 |
|||||||
wm |
wsf |
wm |
wsf |
wm |
wsf |
wm |
wsf |
wm |
wsf |
wm |
wsf |
|
1 |
2,573 |
12,606 |
3,123 |
2,863 |
11,965 |
2,863 |
2,586 |
18,235 |
2,430 |
21,843 |
2,74 |
9,85 |
2 |
2,532 |
14,046 |
3,040 |
2,704 |
12,244 |
2,704 |
2,521 |
13,347 |
1,936 |
25,294 |
2,76 |
11,30 |
3 |
2,394 |
12,500 |
2,870 |
2,610 |
11,173 |
2,610 |
2,723 |
15,169 |
2,384 |
11,764 |
2,88 |
10,72 |
4 |
2,679 |
9,594 |
2,914 |
2,681 |
12,796 |
2,681 |
2,517 |
13,773 |
2,152 |
9,973 |
2,99 |
15,90 |
5 |
2,567 |
6,315 |
3,188 |
2,604 |
10,887 |
2,604 |
2,377 |
12,361 |
2,466 |
7,455 |
2,76 |
12,06 |
6 |
2,174 |
7,331 |
2,871 |
2,497 |
6,934 |
2,497 |
2,588 |
18,433 |
2,444 |
8,806 |
2,81 |
11,59 |
7 |
2,421 |
6,681 |
2,934 |
2,581 |
8,928 |
2,581 |
2,491 |
13,191 |
2,443 |
8,707 |
2,80 |
16,03 |
8 |
2,350 |
10,728 |
2,515 |
2,570 |
9,090 |
2,570 |
2,365 |
14,143 |
2,307 |
6,117 |
2,95 |
13,87 |
9 |
2,381 |
8,823 |
2,705 |
2,776 |
10,114 |
2,776 |
2,398 |
12,828 |
2,294 |
18,449 |
2,94 |
15,18 |
10 |
2,557 |
9,146 |
2,755 |
2,739 |
11,567 |
2,739 |
2,651 |
6,692 |
2,558 |
6,020 |
3,02 |
13,32 |
11 |
2,373 |
10,697 |
2,975 |
2,548 |
6,933 |
2,548 |
2,403 |
18,534 |
1,965 |
13,192 |
3,01 |
14,45 |
12 |
2,417 |
12,078 |
2,921 |
2,638 |
8,146 |
2,638 |
2,525 |
15,957 |
2,247 |
22,743 |
2,77 |
15,34 |
13 |
2,449 |
11,397 |
2,918 |
2,622 |
14,598 |
2,622 |
2,502 |
14,470 |
1,813 |
7,853 |
2,72 |
12,33 |
14 |
2,528 |
8,872 |
2,730 |
2,613 |
6,470 |
2,613 |
2,230 |
6,231 |
2,002 |
6,084 |
3,33 |
13,52 |
15 |
2,394 |
10,272 |
3,631 |
2,437 |
4,910 |
2,437 |
2,262 |
11,548 |
2,107 |
10,810 |
2,80 |
12,03 |
16 |
2,297 |
8,085 |
2,798 |
2,471 |
3,448 |
2,471 |
2,588 |
10,516 |
2,116 |
9,023 |
2,60 |
13,27 |
17 |
2,306 |
10,119 |
2,841 |
2,851 |
9,819 |
2,851 |
2,469 |
7,536 |
2,071 |
17,785 |
2,64 |
7,67 |
18 |
2,447 |
4,680 |
2,812 |
2,544 |
7,784 |
2,544 |
2,348 |
14,114 |
1,989 |
10,424 |
2,60 |
11,97 |
19 |
2,379 |
5,128 |
2,639 |
2,564 |
10,204 |
2,564 |
2,332 |
5,090 |
1,913 |
5,741 |
2,55 |
9,00 |
20 |
2,316 |
10,714 |
2,495 |
2,656 |
9,202 |
2,656 |
2,165 |
10,762 |
2,062 |
6,862 |
2,57 |
10,43 |
21 |
2,242 |
11,974 |
2,646 |
2,589 |
6,115 |
2,589 |
2,232 |
4,815 |
2,155 |
6,854 |
2,51 |
10,67 |
22 |
2,310 |
6,896 |
2,447 |
2,670 |
7,391 |
2,670 |
2,237 |
9,333 |
2,141 |
6,034 |
2,58 |
8,88 |
23 |
2,378 |
7,162 |
2,484 |
2,566 |
5,042 |
2,566 |
2,046 |
16,949 |
2,378 |
9,795 |
2,53 |
8,39 |
24 |
2,331 |
5,023 |
2,528 |
2,506 |
8,043 |
2,506 |
2,548 |
8,385 |
2,432 |
4,954 |
2,49 |
7,62 |
25 |
2,342 |
5,544 |
2,874 |
2,610 |
8,695 |
2,610 |
2,340 |
7,575 |
1,901 |
6,637 |
2,61 |
6,95 |
26 |
2,641 |
4,516 |
2,542 |
2,765 |
16,560 |
2,765 |
2,722 |
4,469 |
1,939 |
7,848 |
2,60 |
17,97 |
27 |
2,394 |
4,871 |
2,504 |
2,441 |
11,326 |
2,441 |
2,326 |
7,079 |
1,826 |
8,787 |
2,63 |
9,09 |
28 |
2,438 |
3,614 |
2,603 |
2,397 |
5,389 |
2,397 |
2,290 |
3,562 |
1,809 |
7,692 |
2,76 |
6,56 |
29 |
2,395 |
7,908 |
2,395 |
2,931 |
2,764 |
2,931 |
2,313 |
12,121 |
1,729 |
6,591 |
2,72 |
7,78 |
30 |
2,330 |
7,808 |
2,109 |
2,896 |
4,304 |
2,896 |
2,588 |
15,300 |
1,680 |
8,158 |
2,57 |
6,74 |
31 |
2,392 |
3,719 |
2,412 |
3,871 |
6,578 |
3,871 |
2,459 |
13,533 |
2,156 |
5,384 |
2,50 |
5,49 |
32 |
2,174 |
5,291 |
2,426 |
2,823 |
3,484 |
2,823 |
2,674 |
7,077 |
1,947 |
12,703 |
2,32 |
5,10 |
33 |
2,386 |
3,977 |
2,398 |
3,432 |
2,380 |
3,432 |
2,547 |
7,155 |
1,648 |
6,230 |
2,46 |
9,63 |
34 |
2,180 |
6,788 |
2,365 |
2,611 |
6,578 |
2,611 |
2,517 |
14,742 |
1,386 |
5,276 |
2,43 |
6,68 |
35 |
2,276 |
2,912 |
2,559 |
3,498 |
0,800 |
3,498 |
2,556 |
7,730 |
1,895 |
19,458 |
2,35 |
7,37 |
36 |
2,179 |
2,836 |
2,697 |
3,245 |
2,222 |
3,245 |
2,403 |
6,796 |
1,788 |
7,950 |
2,46 |
7,16 |
37 |
2,398 |
2,811 |
2,711 |
3,588 |
2,727 |
3,588 |
- |
- |
1,699 |
4,359 |
2,59 |
6,86 |
38 |
2,339 |
4,699 |
2,591 |
4,574 |
11,818 |
4,574 |
- |
- |
1,695 |
4,249 |
2,21 |
9,42 |
39 |
2,507 |
3,676 |
2,734 |
2,917 |
16,772 |
2,917 |
- |
- |
1,913 |
5,357 |
2,30 |
8,45 |
40 |
2,281 |
6,527 |
2,351 |
2,540 |
1,785 |
2,540 |
- |
- |
2,249 |
7,762 |
2,27 |
5,79 |
41 |
2,144 |
4,407 |
2,551 |
2,598 |
2,659 |
2,598 |
- |
- |
2,101 |
7,829 |
2,45 |
4,99 |
42 |
2,190 |
2,739 |
2,397 |
2,174 |
3,278 |
2,174 |
- |
- |
1,813 |
5,555 |
2,37 |
6,47 |
43 |
1,999 |
4,000 |
2,269 |
2,502 |
9,696 |
2,502 |
- |
- |
- |
- |
2,12 |
6,26 |
44 |
2,201 |
2,447 |
2,510 |
2,539 |
7,446 |
2,539 |
- |
- |
- |
- |
2,17 |
9,42 |
45 |
2,265 |
6,229 |
2,409 |
2,334 |
6,341 |
2,334 |
- |
- |
- |
- |
2,16 |
7,49 |
46 |
2,203 |
10,666 |
- |
2,378 |
6,748 |
2,378 |
- |
- |
- |
- |
1,98 |
8,62 |
47 |
1,984 |
3,902 |
- |
2,640 |
5,056 |
2,640 |
- |
- |
- |
- |
- |
- |
48 |
2,239 |
4,332 |
- |
2,849 |
3,401 |
2,849 |
- |
- |
- |
- |
- |
- |
49 |
2,429 |
11,337 |
- |
2,533 |
4,511 |
2,533 |
- |
- |
- |
- |
- |
- |
50 |
2,280 |
13,151 |
- |
2,475 |
4,761 |
2,475 |
- |
- |
- |
- |
- |
- |
51 |
1,977 |
4,310 |
- |
2,441 |
1,785 |
2,441 |
- |
- |
- |
- |
- |
- |
52 |
2,141 |
3,954 |
- |
3,120 |
0,626 |
3,120 |
- |
- |
- |
- |
- |
- |
53 |
2,012 |
7,051 |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
54 |
2,569 |
2,564 |
- |
- |
- |
- |
- |
- |
- |
- |
- |
- |
Приведен регрессионный анализ содержания влаги и распределения по размерам частиц в зависимости от номера разовой пробы, а также регрессионный анализ содержания влаги от распределения по размерам частиц. Распределения по размерам выражают как фракцию железной руды с величиной частиц менее 6 мм. Результаты статистического анализа для этой серии приведены в таблице С.15. Значения вариограммы вычислены для содержания влаги и распределения частиц по размерам (см. таблицы С.16 и С.17). Вариограммы содержания влаги и распределения частиц по размерам показаны на рисунках С.18 и С.19 соответственно.
Таблица С.15 - Результаты статистического анализа для S19 - S24
Номер разовой пробы |
Выбросы |
Регрессия |
Vexp (t = 1) |
|||||
wm |
wsf |
wm - i |
wsf - i |
wm - wsf |
wm |
wsf |
||
S19 |
54 |
0 |
0 |
-а |
-а |
+b |
0,016 |
4,622 |
S20 |
45 |
1 |
0 |
-а |
-b |
0 |
0,034 |
7,142 |
S21 |
52 |
1 |
0 |
0 |
-а |
0 |
0,096 |
8,273 |
S22 |
36 |
0 |
0 |
0 |
-а |
0 |
0,020 |
14,289 |
S23 |
42 |
0 |
1 |
-а |
-а |
0 |
0,033 |
19,242 |
S24 |
46 |
0 |
0 |
-а |
-а |
+а |
0,015 |
4,762 |
wm - массовая доля влаги, выраженная в процентах (% влажности); wsf - фракция железной руды с размерами частиц мене 6 мм; i - номер разовой пробы. |
||||||||
а Очень высокая корреляция. b Высокая корреляция. |
Таблица С.16 - Значения вариограммы содержания влаги для S19 - S24
S19 |
S20 |
S21 |
S22 |
S23 |
S24 |
|
0 |
0 |
0 |
0 |
0 |
0 |
0 |
1 |
0,016 |
0,034 |
0,096 |
0,020 |
0,033 |
0,015 |
2 |
0,017 |
0,034 |
0,105 |
0,020 |
0,038 |
0,021 |
3 |
0,019 |
0,040 |
0,142 |
0,024 |
0,048 |
0,020 |
4 |
0,016 |
0,044 |
0,162 |
0,025 |
0,044 |
0,022 |
5 |
0,015 |
0,060 |
0,155 |
0,026 |
0,049 |
0,026 |
6 |
0,016 |
0,056 |
0,175 |
0,020 |
0,061 |
0,028 |
7 |
0,019 |
0,067 |
0,175 |
0,025 |
0,055 |
0,034 |
8 |
0,017 |
0,055 |
0,210 |
0,025 |
0,056 |
0,035 |
9 |
0,014 |
0,059 |
0,221 |
0,026 |
0,048 |
0,045 |
10 |
0,015 |
0,037 |
0,243 |
0,033 |
0,071 |
0,043 |
11 |
0,016 |
0,054 |
0,254 |
0,038 |
0,065 |
0,041 |
12 |
0,018 |
0,059 |
0,245 |
0,041 |
0,066 |
0,048 |
13 |
0,018 |
0,059 |
0,254 |
0,034 |
0,066 |
0,051 |
14 |
0,021 |
0,065 |
0,238 |
0,025 |
0,058 |
0,046 |
15 |
0,019 |
0,085 |
0,200 |
0,030 |
0,059 |
0,054 |
16 |
0,020 |
0,076 |
0,185 |
0,028 |
0,054 |
0,060 |
17 |
0,022 |
0,089 |
0,175 |
0,036 |
0,057 |
0,065 |
18 |
0,024 |
0,091 |
0,182 |
0,034 |
0,075 |
0,069 |
Таблица С.17 - Значения вариограммы распределения частиц по размерам для S19 - S24
S19 |
S20 |
S21 |
S22 |
S23 |
S24 |
|
0 |
0 |
0 |
0 |
0 |
0 |
0 |
1 |
4,622 |
7,142 |
8,273 |
14,289 |
19,242 |
4,762 |
2 |
7,297 |
7,257 |
13,739 |
15,689 |
23,312 |
4,980 |
3 |
6,693 |
8,041 |
15,001 |
16,430 |
19,360 |
5,167 |
4 |
6,414 |
10,029 |
12,404 |
16,357 |
26,790 |
5,376 |
5 |
7,178 |
8,421 |
11,341 |
12,479 |
20,857 |
5,538 |
6 |
7,091 |
9,407 |
13,438 |
17,028 |
22,840 |
6,375 |
7 |
6,940 |
7,722 |
12,358 |
10,500 |
25,784 |
5,863 |
8 |
6,726 |
6,227 |
13,612 |
18,821 |
19,231 |
6,733 |
9 |
5,264 |
7,545 |
16,007 |
18,471 |
21,019 |
7,806 |
10 |
6,394 |
9,316 |
17,775 |
17,994 |
17,724 |
7,881 |
11 |
7,342 |
8,996 |
16,070 |
12,758 |
19,784 |
7,996 |
12 |
7,345 |
9,156 |
12,561 |
14,314 |
24,406 |
7,706 |
13 |
8,150 |
9,736 |
12,954 |
26,553 |
24,739 |
9,732 |
14 |
9,781 |
10,081 |
14,041 |
18,035 |
20,668 |
9,973 |
15 |
10,185 |
8,542 |
13,968 |
24,037 |
17,418 |
11,249 |
16 |
10,100 |
6,340 |
12,958 |
21,422 |
21,631 |
12,257 |
17 |
9,494 |
4,824 |
10,605 |
19,795 |
25,140 |
11,569 |
18 |
9,945 |
7,389 |
12,154 |
20,286 |
23,557 |
12,213 |
Рисунок С.18 - Вариограмма содержания влаги для S19 - S24
Рисунок С.19 - Вариограмма распределения частиц для S19 - S24
С.7.1 Сбор данных
Для железной руды изменения характеристики качества в пределах слоя, из которого отбирают две разовые пробы, проверяют периодически с помощью систематического отбора выборок на основе массы в соответствии с ИСО 3084 [16]. Методом наблюдения является отбор сквозной выборки (см. 7.3). На основе С.5 и С.6 отобраны и обработаны данные последней декады. Первоначальные данные приведены в таблице 5 (см. 7.3) и повторены для удобства в таблице С.18. Данные для крупных и мелких фракций объединены в таблице С.19.
Таблица С.18 - Пример отбора сквозной выборки
Общее содержание железа |
||||
А |
В |
Среднее арифметическое |
Размах |
|
1 |
65,37 |
64,36 |
64,865 |
1,01 |
2 |
64,82 |
64,82 |
64,82 |
0 |
3 |
64,81 |
65,10 |
64,955 |
0,29 |
4 |
64,96 |
65,06 |
65,01 |
0,10 |
5 |
65,23 |
65,20 |
65,215 |
0,03 |
6 |
65,34 |
65,22 |
65,28 |
0,12 |
7 |
65,54 |
65,80 |
65,67 |
0,26 |
8 |
65,41 |
65,34 |
65,375 |
0,07 |
9 |
65,16 |
65,22 |
65,19 |
0,06 |
10 |
65,34 |
65,69 |
65,515 |
0,35 |
Среднее арифметическое |
65,19 |
0,23 |
||
swst |
0,35 |
|||
Общее содержание железа: метод - 3´2´10; сталелитейный завод - G; дата - 1985-05-19; масса - 97101 т. |
С.7.2 Распределение по размерам частиц
Исключим L04, для которого определяют долю железной руды с размером частиц менее 8 мм вместо доли руды с размером частиц менее 6 мм, L08, для которого распределение по размерам частиц не рассматривается, и L10, для которого включено 12 частей. Другие серии данных для крупных фракций, приведенные в таблице С.19, включают одну и ту же модель эксперимента для фракции с размером частиц менее 6 мм, которая содержит 10 частей, Для каждой из этих частей составлено две сложных пробы, независимо от того, что каждая сложная проба содержит различное количество разовых проб. Соответственно, девять партий проанализированы вместе. Каждое значение данных для сложных проб хВС можно представить следующим образом
хВС = m + fBV + fBP + fBC, (С.8)
где m - математическое ожидание совокупности;
fBV - описывает изменчивость между судами;
fBP - описывает изменчивость между частями судна;
fBC - описывает изменчивость между сложными пробами в пределах части.
Таблица С.19 - Обзор изменчивости характеристик качества
Масса партии |
Метод |
t |
Рассмотренные характеристики качества |
|
|
L01 |
69597 |
2×2×10 = 40 |
1740 |
wsf, wFe |
|
L02 |
74141 |
5×2×10 = 100 |
741 |
wsf, wFe |
|
L03 |
77749 |
5×2×10 = 100 |
777 |
wsf, wFe |
|
L04 |
60820 |
3×2×10 = 60 |
1014 |
wsf, wFe |
|
L05 |
144583 |
3×2×10 = 60 |
2410 |
wsf, wm, wFe |
|
L06 |
135626 |
4×2×10 = 80 |
1695 |
wsf, wm, wFe |
|
L07 |
93884 |
3×2×10 = 60 |
1565 |
wsf, wm, wFe |
|
L08 |
78334 |
3×2×10 = 60 |
1306 |
wFe |
|
L09 |
97101 |
3×2×10 = 60 |
1618 |
wsf, wFe |
|
L10 |
50532 |
4×2×2 = 96 |
526 |
wsf, wm, wFe |
|
L11 |
111132 |
3×2×10 = 60 |
1852 |
wsf, wm, wFe |
|
L12 |
90892 |
5×2×10 = 100 |
909 |
wsf, wm, wFe |
|
f01 |
149081 |
3×2×10 = 60 |
2485 |
wFe |
|
f02 |
109759 |
4×2×10 = 80 |
1372 |
wsf, wm, wFe, wSi, wAl |
|
f03 |
44467 |
2×2×10 = 40 |
1112 |
wsf, wFe |
|
f04 |
55725 |
3×2×10 = 60 |
929 |
wFe |
|
f05 |
120073 |
4×2×10 = 80 |
1501 |
wsf, wm, wFe, wSi, wAl |
|
f06 |
103051 |
3×2×10 = 60 |
1718 |
wFe |
|
f07 |
66691 |
2×2×10 = 40 |
1667 |
wm, wFe |
|
f08 |
103956 |
3×2×10 = 60 |
1733 |
wFe |
|
f09 |
102398 |
3×2×10 = 60 |
1707 |
wsf, wm, wFe |
|
f10 |
60418 |
4×2×10 = 80 |
755 |
wsf, wm, wFe |
|
f11 |
30130 |
4×2×11 = 88 |
342 |
wsf, wm, wFe |
|
f12 |
31446 |
2×2×12 = 48 |
655 |
wsf, wm, wFe |
|
L - соответствует крупной фракции и сопровождается номером партии; f - соответствует мелкой фракции и сопровождается номером партии; wsf - распределение по размерам (фракциям (%) с размером частиц > 6 мм для крупной фракции и < 6 мм для мелкой фракции); wm - доля массы влаги, выраженная в процентах (%), (влажность); wFe - доля массы железа, выраженная в процентах (%), (содержание железа); wsi - доля массы SiO2, выраженная в процентах (%), (содержание SiO2); wAl - доля массы Al2O3, выраженная в процентах (%), (содержание Al2O3). |
|||||
В этом случае количество судов - девять, количество частей - 10 и количество сложных проб - две. Таблица ANOVA приведена в таблице С.20.
Таблица С.20 - Таблица ANOVA для данных изменчивости характеристики качества
Сумма квадратов отклонений |
Число степеней свободы |
Средний квадрат отклонений |
Математическое ожидание среднего квадрата отклонений |
|
Суда |
SBV |
v - 1 |
SBV/fBV |
|
Части |
SBP |
v(p - 1) |
SBP/fBP |
|
Сложные пробы |
SBC |
vp(c - 1) |
SBC/fBC |
|
Сумма |
St |
vpc - 1 |
Таблица С.21 получена из таблицы С.20 путем подстановки числовых значений величин.
Таблица С.21 - Таблица ANOVA для распределения по размерам крупных фракций
Сумма квадратов отклонений |
Число степеней свободы |
Средний квадрат отклонений |
Математическое ожидание среднего квадрата отклонений |
|
Суда |
572,14 |
8 |
71,52 |
|
Части |
1289,59 |
81 |
15,92 |
|
Сложные пробы |
153,82 |
90 |
1,71 |
|
Сумма |
2015,55 |
179 |
Приравнивание выражения в последней колонке таблицы С.21 к значению среднего квадрата отклонений, позволяет определить соответствующие дисперсии
Следовательно,
Таким образом, дисперсия, соответствующая частям судна, дает наибольший вклад в исследуемую изменчивость.
Результаты ANOVA для данных мелких фракций пяти судов, включая f02, f03, f05, f09 и f10, приведены в таблице С.22.
Таблица С.22 - Таблица ANOVA для распределения по размерам мелких фракций
Сумма квадратов отклонений |
Число степеней свободы |
Средний квадрат отклонений |
Математическое ожидание среднего квадрата отклонений |
|
Суда |
462,65 |
4 |
115,66 |
|
Части |
142,83 |
45 |
3,17 |
|
Сложные пробы |
52,59 |
50 |
1,05 |
|
Сумма |
658,07 |
99 |
Таким образом:
Наибольшей является дисперсия, характеризующая изменчивость.
Линия регрессии, описывающая распределение по размерам для относительного номера части, построена для каждого судна. Относительный номер равен отношению номера части к общему количеству частей. Если общее количество частей равно единице, угловой коэффициент линии регрессии указывает разность между процентом фракции в начале и в конце разгрузки. Угловой коэффициент линии регрессии распределения по размерам приведен в таблице С.23. В таблице С.23 приведены дисперсии и полученные с помощью ANOVA для каждого судна.
Таблица С.23 - Угловой коэффициент линии регрессии для распределения по размерам
Угловой коэффициент линии регрессии |
|
|
|
L01 |
-11,895а |
3,122 |
13,19 |
L02 |
-7,058а |
1,936 |
4,57 |
L03 |
-4,150а |
0,717 |
2,30 |
L05 |
-11,468а |
1,996 |
12,29 |
L06 |
-6,684а |
0,742 |
5,00 |
L07 |
-9,373а |
2,205 |
8,88 |
L09 |
-6,371а |
2,878 |
8,54 |
L10 |
-9,924а |
0,879 |
9,83 |
L11 |
-7,098а |
0,547 |
4,63 |
L12 |
-5,318а |
1,239 |
4,56 |
f02 |
+1,71а |
0,437 |
0,159 |
f03 |
+2,69 |
0,599 |
3,314 |
f02 |
+1,14 |
0,647 |
0,428 |
f02 |
-0,06 |
2,494 |
0,992 |
f10 |
-1,86b |
1,082 |
0,413 |
f11 |
+1,55 |
1,911 |
0,250 |
f12 |
+1,85а |
0,271 |
0,551 |
a Очень высокая корреляция. b Высокая корреляция. |
С.7.3 Влажность
Результаты ANOVA для влажности крупных фракций, включая L05, L06, L07, L11 и L12, приведены в таблице С.24.
Таблица С.24 - Таблица ANOVA для влажности крупных фракций
Сумма квадратов отклонений |
Число степеней свободы |
Средний квадрат отклонений |
Математическое ожидание среднего квадрата отклонений |
|
Суда |
0,44 |
4 |
0,11 |
|
Части |
5,91 |
45 |
0,13 |
|
Сложные пробы |
1,2 |
50 |
0,02 |
|
Сумма |
7,55 |
99 |
Приравнивая выражения из последней колонки таблицы С.24 к полученным значениям среднего квадрата отклонений, определяют дисперсии
Наибольшая дисперсия соответствует изменчивости между частями судна.
Результаты ANOVA для влажности мелких фракций, включая f02, f05, f07, f09 и f10, приведены в таблице С.25.
Таблица С.25 - Таблица ANOVA для влажности мелких фракций
Сумма квадратов отклонений |
Число степеней свободы |
Средний квадрат отклонений |
Математическое ожидание среднего квадрата отклонений |
|
Суда |
3,43 |
4 |
0,86 |
|
Части |
0,66 |
45 |
0,01 |
|
Сложные пробы |
0,34 |
50 |
0,01 |
|
Сумма |
4,43 |
99 |
Таким образом:
Наибольшая дисперсия соответствует изменчивости между судами. Угловой коэффициент линии регрессии и приведен в таблице С.26.
Таблица С.26 - Угловой коэффициент линии регрессии для влажности
Угловой коэффициент линии регрессии |
|
|
|
L05 |
+0,061 |
0,074 |
0,019 |
L06 |
-0,181 |
0,004 |
0,012 |
L07 |
-0,455a |
0,007 |
0,050 |
L10 |
+0,032 |
0,006 |
0,011 |
L11 |
+0,731a |
0,032 |
0,104 |
L12 |
-0,915a |
0,003 |
0,084 |
f02 |
+0,022 |
0,004 |
0,006 |
f05 |
+0,076 |
0,006 |
0 |
f07 |
-0,336a |
0,011 |
0,017 |
f09 |
+0,030 |
0,004 |
0 |
f10 |
+0,011 |
0,009 |
0 |
f11 |
+0,110 |
0,025 |
0,278 |
f12 |
-0,287a |
0,003 |
0,010 |
a Очень высокая корреляция. |
C.7.4 Общее содержание железа
Результаты ANOVA для общего содержания железа по данным крупной фракции 10 судов, исключая L04 (другой тип руды) и L10 (другое количество частей), приведены в таблице С.27.
Таблица С.27 - Таблица ANOVA для общего содержания железа в крупных фракциях
Сумма квадратов отклонений |
Число степеней свободы |
Средний квадрат отклонений |
Математическое ожидание среднего квадрата отклонений |
|
Суда |
21,98 |
9 |
2,44 |
|
Части |
9,89 |
90 |
0,11 |
|
Сложные пробы |
8,43 |
100 |
0,08 |
|
Сумма |
40,30 |
199 |
Таким образом:
Наибольшая дисперсия соответствует изменчивости между судами.
Результаты ANOVA для общего содержания железа в мелкой фракции для 10 судов, исключая f11 и f12 (оба имеют различное количество частей), приведены в таблице С.28.
Таблица С.28 - Таблица ANOVA для общего содержания железа в мелких фракциях
Сумма квадратов отклонений |
Число степеней свободы |
Средний квадрат отклонений |
Математическое ожидание среднего квадрата отклонений |
|
Суда |
25,20 |
9 |
2,80 |
|
Части |
5,71 |
90 |
0,06 |
|
Сложные пробы |
2,70 |
100 |
0,03 |
|
Сумма |
33,61 |
199 |
Таким образом:
Наибольшая дисперсия соответствует изменчивости между судами.
Угловой коэффициент линии регрессии для общего содержания железа приведен в таблице С.29. Данные для L07 и L11 относятся к одному заводу D и дают линии регрессии с отрицательным угловым коэффициентом. Бульдозеры обычно используют на самой последней стадии для быстрой разгрузки железной руды с судна. Завод D использует бульдозеры на более ранней стадии для разгрузки судна, что приводит к измельчению частиц руды, перемещению их в верхнюю часть массива при разгрузке судна и последующему перемешиванию руды.
Таблица С.29 - Угловой коэффициент линии регрессии для общего содержания железа
Угловой коэффициент линии регрессии |
|
|
|
L01 |
+0,455 |
0,116 |
0,004 |
L02 |
+0,658а |
0,040 |
0,042 |
L03 |
+0,344 |
0,073 |
0,044 |
L04 |
+0,868а |
0,005 |
0,083 |
L05 |
+0,573 |
0,289 |
0 |
L06 |
+0,266а |
0,013 |
0,002 |
L07 |
-0,140 |
0,083 |
0 |
L08 |
+0,849а |
0,098 |
0,033 |
L09 |
+0,763а |
0,066 |
0,047 |
L10 |
+0,066 |
0,048 |
0 |
L11 |
-0,231 |
0,044 |
0,002 |
L12 |
-0,102 |
0,020 |
0,019 |
f01 |
+0,395b |
0,023 |
0,039 |
f02 |
-0,226b |
0,021 |
0,003 |
f03 |
+0,432 |
0,070 |
0,062 |
f04 |
-0,128 |
0,007 |
0,002 |
f05 |
-0,122 |
0,014 |
0 |
f06 |
-0,441а |
0,019 |
0,029 |
f07 |
-0,045 |
0,022 |
0,003 |
f08 |
-0,249 |
0,041 |
0,026 |
f09 |
+0,050 |
0,023 |
0,029 |
f10 |
+0,050 |
0,030 |
0 |
f11 |
-0,360а |
0,018 |
0,009 |
f12 |
+ 0,113 |
0,022 |
0 |
a Очень высокая корреляция. b Высокая корреляция. |
С.7.5 Смещение при отборе проб железной руды
Тенденция распределения по размерам (фракция с размером частиц менее 6 мм) составляет минус 12 % от максимума для крупной фракции и минус 2 % для мелкой фракции с размером частиц менее 6 мм (в таблице С.23, около плюс 2 % для фракции крупнее 6 мм).
Тенденция изменения влажности из таблицы С.26 не ясна. Однако в соответствии с таблицей С.9 она составляет не менее 1,0 % для крупной фракции и менее 0,3 % для мелкой фракции.
Тенденция общего содержания железа составляет менее 1,0 % для крупной фракции (см. таблицу С.29), в то время как тенденция изменений для мелкой фракции не наблюдается.
В случае максимального тренда (- 12 %) распределения по размерам для крупной фракции возможно смещение, которое вызвано систематическим отбором выборок на основе массы из 50 разовых проб, что составляет 12/(2×50) = 0,12 % фракции с размером частиц менее 6 мм. Также возможно смещение для влажности и общего содержания железа 1,0/(2×50) = 0,01 %. Эти возможные смещения незначительны по сравнению с общим стандартным отклонением оценки среднего характеристики качества партии.
D.1 Общие положения
В настоящем приложении установлена процедура проведения эксперимента и описано ее практическое применение для оценки прецизионности при выполнении стандартного отбора выборки, подготовки пробы и измерений. Прецизионность, соответствующая стандартному отбору выборки, подготовке пробы и выполнению измерений, необходимо периодически проверять на основе двойного отбора проб. Если метод стандартного отбора выборки (включая подготовку пробы и выполнение измерений) вводят впервые или он изменен, или прошло значительное время с его последней проверки, должен быть проведен специальный эксперимент. Этот эксперимент должен соответствовать одному из стандартов на выборочный контроль нештучной продукции.
Стандартный отбор выборки нештучной продукции часто выполняют путем систематического отбора выборки на основе массы. Если систематический отбор выборки на основе массы применяют к сыпучим материалам, стандартный отбор выборки выполняют в соответствии с одним из следующих планов эксперимента:
План 1. Партию делят на ulot частей партии, все разовые пробы отбирают из каждой части партии и соединяют в ulot объединенных проб, подготавливают ulot исследуемых проб, которые анализируют отдельно. Взвешенное среднее результатов исследований в соответствии с массой частей партии является оценкой среднего характеристики качества партии. Когда масса партии меньше установленной массы части партии, то партия не может быть разделена на части. В этом случае все разовые пробы отбирают из неразделенной партии и соединяют в объединенную пробу, подготавливают исследуемую пробу и выполняют ее анализ. Результат анализа - оценка среднего качества партии.
План 2. Несколько разовых проб, отобранных из партии, соединяют в пробы, подготавливают исследуемые пробы и анализируют их отдельно. Взвешенное среднее результатов исследований в соответствии с количеством разовых проб, объединенных в выборки, является оценкой среднего характеристики качества партии. Обычно в выборки объединяют равное количество разовых проб.
План 3. Для всех разовых проб, отобранных из партии, последовательно измеряют характеристики качества (после подготовки, при необходимости). Выборочное среднее результатов исследований является оценкой среднего характеристики качества партии.
Каждая из описанных процедур позволяет получить оценку среднего характеристики качества партии с соответствующей плану прецизионностью. Достигнутую прецизионность можно проверить с помощью двойного отбора проб.
В случае систематического отбора выборки на основе времени применимы процедуры плана 1, но процедуры планов 2 и 3 не могут быть применены при пропорциональном изменении разовой пробы по мере поступления материала.
D.2 Общие условия
D.2.1 Количество экспериментов
Следует выполнять эксперименты не менее чем на 20 партиях материала. Если это невозможно, следует выполнять эксперимент не менее чем на 10 партиях.
D.2.2 Количество разовых сложных проб
Минимальное количество разовых проб, необходимое для эксперимента, должно быть вдвое больше количества разовых проб, установленного для стандартного отбора выборки. Например, если партия разделена на ulot частей партии, количество разовых проб, отобранных из части партии, равно nsub. Из этих разовых проб составляют одну объединенную пробу. При стандартном отборе выборки минимальное количество разовых проб, необходимое для эксперимента, должно быть 2ulot, из которых составляют две сложные пробы (каждая сложная проба соответствует объединенной пробе при стандартном отборе выборки) по nsub разовых проб в каждой.
D.2.3 Подготовка пробы и выполнение измерений
Подготовка пробы и выполнение измерений должны соответствовать методу, установленному в одном из стандартов.
D.3 Порядок проведения эксперимента
В экспериментах с использованием систематического отбора выборок на основе массы должно быть отобрано вдвое большее количество разовых проб, чем при стандартном отборе выборки, т.е. 2nsub разовых проб из части партии через половину объема массы для стандартного отбора, из которых составляют две сложные пробы по nsub разовых проб в каждой. При этом разовые пробы включают в сложные пробы попеременно по мере отбора.
В экспериментах с использованием отбора простой случайной выборки должно быть отобрано вдвое большее количество разовых проб, чем при стандартном отборе выборки, т.е. 2n разовых проб. Разовые пробы отбирают из партии случайным образом и составляют из них две сложные пробы по п разовых проб в каждой. Разовые пробы для включения в сложную пробу также отбирают случайным образом.
В эксперименте с использованием расслоенной выборки должно быть отобрано вдвое большее количество разовых проб, чем при стандартном отборе выборки. Разовые пробы отбирают случайным образом из каждого слоя с равным количеством разовых проб в слое, равным количеству разовых проб при стандартном отборе выборки. Затем из разовых проб формируют две сложные пробы для частей партии или всей партии.
В экспериментах с использованием двухэтапного отбора выборки из партии случайным образом должно быть отобрано сначала usub выборочных единиц как при стандартном отборе выборки. Затем таким же способом отбирают еще usub выборочных единиц. Выборочные единицы, отобранные на втором этапе, могут отбираться одновременно с выборочными единицами первого этапа. Объединяют все разовые пробы каждого этапа отдельно и составляют две сложные пробы для партии.
В экспериментах каждого типа должны быть подготовлены две сложные пробы, на которых выполняют измерения отдельно в соответствии с процедурами, установленными в соответствующем стандарте.
D.4 Статистический анализ полученных данных
Независимо от количества двойных данных могут быть применены верхняя контрольная граница UCL и центральная линия карты размахов, описанные в ИСО 8258 [29], для обнаружения точек, выходящих за границу и оценки достигнутой прецизионности. Для двойного отбора проб коэффициент D3 карты размахов равен нулю, а нижняя контрольная граница отсутствует. Коэффициент D4 карты размахов равен 3,267, а коэффициент d2 равен 1,128 для двойных данных. Если ни одно из значений размаха не превышает границу UCL, то можно считать, что процесс находится в состоянии статистической управляемости.
На основе двойных данных вычисляют средний размах . Верхняя контрольная граница для карты размахов равна По расположению размахов на карте размахов можно судить о стабильности процесса. Если процесс находится в состоянии статистической управляемости, дает оценку стандартного отклонения по двойным данным, т.е. оценку достигнутой прецизионности.
D.5 Практическое применение
Влажность (%) и распределение по размерам (фракция (%) железной руды с размером частиц менее 6 мм) оценивали по 2n разовым пробам, отобранных из половины интервала массы для стандартного отбора выборки, как показано в таблице D.1, где измерения п разовых проб выполнялись последовательно при стандартном систематическом отборе выборок на основе массы.
Таблица D.1 - Влажность и распределение по размерам по двойным данным
Влажность, % |
Распределение по размерама), % |
|
1 |
2,573 2,532 |
12,606 14,046 |
2 |
2,394 2,679 |
12,500 9,594 |
3 |
2,567 2,174 |
6,315 7,331 |
4 |
2,421 2,350 |
6,681 10,728 |
5 |
2,381 2,557 |
8,823 9,146 |
6 |
2,373 2,417 |
10,697 12,078 |
7 |
2,449 2,528 |
11,397 8,872 |
8 |
2,394 2,297 |
10,272 8,085 |
9 |
2,306 2,447 |
10,119 4,680 |
10 |
2,379 2,316 |
5,128 10,714 |
11 |
2,242 2,310 |
11,974 6,896 |
12 |
2,378 2,331 |
7,162 5,023 |
13 |
2,342 2,641 |
5,544 4,516 |
14 |
2,394 2,438 |
4,871 3,614 |
15 |
2,395 2,330 |
7,908 7,808 |
16 |
2,392 2,174 |
3,719 5,291 |
17 |
2,386 2,180 |
3,977 6,788 |
18 |
2,276 2,179 |
2,912 2,836 |
19 |
2,398 2,339 |
2,811 4,699 |
20 |
2,507 2,281 |
3,676 6,527 |
21 |
2,144 2,190 |
4,407 2,739 |
22 |
1,999 2,201 |
4,000 2,447 |
23 |
2,265 2,203 |
6,229 10,666 |
24 |
1,984 2,239 |
3,902 4,332 |
25 |
2,429 2,280 |
11,337 13,151 |
26 |
1,977 2,141 |
4,310 3,954 |
а) Фракция железной руды с размером частиц менее 6 мм. |
Обе контрольные карты размаха (на влажность и распределение частиц по размерам) демонстрируют состояние статистической управляемости. Средние размахи равны 0,138 и 2,150 соответственно. Математическое ожидание стандартного отклонения оценок среднего характеристики качества партии при стандартном отборе выборки имеет вид.
Для влажности
где - стандартное отклонение в пределах слоя двойных данных;
- стандартное отклонение оценки среднего влажности партии.
Для распределения по размерам
В приложении для железной руды партия разделена на 10 частей и пять разовых проб отобраны из каждой части партии на основе массы при систематическом отборе выборок в соответствии со стандартной процедурой. Для оценки прецизионности 10 разовых проб отобраны из каждой части партии через половину интервала массы для стандартного отбора выборки при систематическом отборе выборок на основе массы и объединены в две сложные пробы поочередно по мере отбора разовых проб. Каждую сложную пробу готовят отдельно и проводят анализ для определения общего содержания железа. Результаты приведены в таблице D.2.
Таблица D.2 - Общее содержание железа для оценки прецизионности
Общее содержание железа, % |
Ri |
|
1 |
65,17 65,54 |
0,37 |
2 |
65,19 65,42 |
0,23 |
3 |
65,23 65,38 |
0,15 |
4 |
65,58 65,63 |
0,05 |
5 |
65,33 65,51 |
0,18 |
6 |
65,38 65,59 |
0,21 |
7 |
65,48 65,23 |
0,25 |
8 |
65,80 65,57 |
0,23 |
9 |
65,18 65,13 |
0,05 |
10 |
65,14 65,16 |
0,02 |
Сумма |
1,74 |
Верхняя контрольная граница карты размахов: Карта размахов демонстрирует состояние статистической управляемости процесса. Математическое ожидание стандартного отклонения оценки среднего характеристики качества партии при стандартном отборе выборки имеет вид
Примечание - Экспериментальные методы и числовые примеры для проверки прецизионности отбора выборки из железной руды приведены в ИСО 3085 [17].
Е.1 Общие положения
В настоящем приложении приведены методы и рекомендации для проверки смещения в процессе отбора выборки нештучной продукции. Методы оценки случайных ошибок, возникающих в исследуемой системе*), позволяют проверить смещение исследуемой системы по отношению к эталонной системе с помощью единственного эксперимента.
_____________
* Система выборочного контроля.
Приведенные методы применимы к отдельным компонентам исследуемой системы отбора проб включая магистральные квантизаторы, дробилки, делители и средства измерений, а также применимы к исследуемой системе в целом.
В случае двойных экспериментов непосредственно по результатам эксперимента могут быть получены среднее арифметическое характеристики качества партии, математическое ожидание дисперсии оценки среднего характеристики качества партии и другие дисперсии. Однако в случае двойного отбора проб по результатам эксперимента могут быть получены только дисперсии, соответствующие отдельным компонентам системы отбора проб. Соответственно, математическое ожидание дисперсии оценки среднего характеристики качества партии на основе стандартного отбора выборки, подготовки пробы, выполнения измерений должно быть получено методами, приведенными в разделе 8 приложения D.
Различные методы для проверки смещения приведены в ИСО 3086 [18], ИСО 3301 [19], а также в [20] - [26]. В настоящем стандарте предложен метод на основе двойных наблюдений исследуемой и эталонной системы.
Е.2 Принцип проверки наличия смещения
Обычно проверка смещения системы отбора выборки основана на отборе серии пар проб по существу одного и того же материала. Один элемент каждой пары отбирается исследуемой системой или ее компонентом. Другой элемент пары отбирают несмещенным методом посредством эталонной системы.
В методах, предложенных в настоящем приложении, два элемента каждого набора отбирают с помощью системы или компонента, а другие два элемента отбирают с помощью несмещенного эталонного метода. Дисперсии ошибки исследуемой системы оценивают по двум элементам, отобранным исследуемой системой. Проверка смещения, соответствующего системе отбора выборки, основана на сравнении среднего арифметического данных исследуемой и эталонной систем. Количество данных должно быть определено экспериментально.
Е.3 Процедуры проверки смещения
Отбирают не менее 20 наборов проб для исследований. При этом в каждом наборе два элемента отбирают с применением исследуемой системы и два элемента отбирают несмещенным эталонным методом. Рекомендуемые действия при отборе проб следующие:
a) отбирают первую пробу с применением исследуемой системы;
b) отбирают первую пробу эталонным методом;
c) отбирают вторую пробу с применением исследуемой системы;
d) отбирают вторую пробу эталонным методом и т.д.
Однако при исследовании механического устройства для отбора проб, допускается порядок, при котором первую и вторую пробы отбирают два раза, сначала с применением исследуемой системы, а затем - эталонным методом. Пробы в наборе должны быть отобраны как можно ближе друг к другу насколько позволяет конвейер.
Для четырех элементов каждого набора результаты измерений определяют отдельно и затем проводят статистический анализ.
Номер набора данных |
Исследуемая система |
Эталонный метод |
1 |
x1,1 x1,2 |
y1,1 y1,2 |
2 |
x2,1 x2,2 |
y2,1 y2,2 |
... |
... ... |
... ... |
20 |
x20,1 x20,2 |
y20,1 y20,2 |
Рисунок Е.1 - Схема эксперимента для проверки смещения
Е.4 Анализ смещения
Если i - номер набора данных, хi1, хi2 - результаты измерений проб i-го набора проб, отобранных с применением исследуемой системы, уi1, уi2 - результаты измерений проб i-го набора проб, отобранных эталонным методом, где i = 1, 2, ..., k (k = 20).
Вычисляют следующие значения
gi = xi1 - xi2, (E.1)
hi = yi1 - yi2, (Е.3)
Вычисляют отношение F0 дисперсий ошибки, помещая большее из и в числитель
или (E.5)
Если F0 > Fa/2(v1, v2), то нулевую гипотезу отклоняют, следовательно две группы данных нельзя считать принадлежащими одной совокупности с общей дисперсией. Уровень значимости а обычно устанавливают равным 0,05, a v1 и v2 - число степеней свободы и и соответственно, и оба в рассматриваемом случае равны k.
Если F0 ≤ Fa/2(v1, v2), можно считать, что эти две группы данных имеют общую дисперсию. Доверительные границы Т1(х), T2(х) и Т1(у), Т2(у) для уровня доверия 95 % вычисляют следующим образом
где - общее среднее арифметическое хi1 и хi2;
- общее среднее арифметическое уi1 и уi2.
Вычисляют
(E.8)
(E.9)
(E.10)
Если то гипотезу равенства двух способов отбора проб отклоняют. Оценка смещения равна .
Е.5 Содержание остаточного углерода
Характеристикой качества является остаточное содержание углерода после перегонки 10 % тяжелой нефти сорта А. Эталонный метод - метод микроперегонки требует около трех часов для получения результатов исследований.
Метод исследуемой системы - составление мультиуравнений регрессии по результатам спектрофотометрии. Результат таких исследований может быть получен в течение 10 мин. Результаты эксперимента приведены в таблице Е.1. Результаты расчета и приведены в таблице Е.2.
Таблица Е.1 - Результаты эксперимента по определению остаточного углерода в тяжелой нефти
Данные представляют собой процент массы, умноженный на 103
Исследуемая система |
Эталонный метод |
|||
xi1 |
xi2 |
yi1 |
yi2 |
|
1 |
335 |
335 |
318 |
343 |
2 |
343 |
342 |
310 |
307 |
3 |
340 |
342 |
356 |
300 |
4 |
355 |
357 |
372 |
329 |
5 |
346 |
347 |
318 |
332 |
6 |
348 |
346 |
302 |
320 |
7 |
356 |
359 |
325 |
301 |
8 |
351 |
353 |
344 |
351 |
9 |
343 |
343 |
326 |
317 |
10 |
351 |
351 |
320 |
315 |
11 |
315 |
314 |
296 |
276 |
12 |
329 |
329 |
282 |
283 |
13 |
334 |
334 |
280 |
292 |
14 |
330 |
330 |
353 |
317 |
15 |
331 |
331 |
261 |
312 |
16 |
345 |
345 |
276 |
300 |
17 |
345 |
344 |
280 |
304 |
18 |
333 |
334 |
283 |
317 |
19 |
359 |
359 |
313 |
350 |
20 |
353 |
354 |
299 |
328 |
Сумма |
6842 |
6849 |
6214 |
6294 |
Среднее арифметическое |
342,275 |
312,700 |
Таблица Е.2 - Вычисление и
gi = xi1 - xi2 |
|
hi = yi1 - yi2 |
|
|
1 |
0 |
0 |
-25 |
625 |
2 |
1 |
1 |
3 |
9 |
3 |
-2 |
4 |
56 |
3136 |
4 |
-2 |
4 |
43 |
1849 |
5 |
-1 |
1 |
-14 |
196 |
6 |
2 |
4 |
-18 |
324 |
7 |
-3 |
9 |
24 |
576 |
8 |
-2 |
4 |
-7 |
49 |
9 |
0 |
0 |
9 |
81 |
10 |
0 |
0 |
5 |
25 |
11 |
1 |
1 |
20 |
400 |
12 |
0 |
0 |
-1 |
1 |
13 |
0 |
0 |
-12 |
144 |
14 |
0 |
0 |
36 |
1296 |
15 |
0 |
0 |
-51 |
2601 |
16 |
0 |
-24 |
576 |
|
17 |
0 |
1 |
-24 |
576 |
18 |
1 |
1 |
-34 |
1156 |
19 |
-1 |
0 |
-37 |
1369 |
20 |
0 |
1 |
-29 |
841 |
Сумма |
31 |
15830 |
В соответствии с уравнениями (Е.2) и (Е.4):
se(x) = 0,880,
se(y) = 19,893,
Fa/2(20, 20) = 2,46.
Поскольку Fa/2 < F0, две группы данных нельзя считать принадлежащими к одной совокупности с общей дисперсией.
Доверительные границы Т1(х), Т2(х) и Т1(у), Т2(у) с уровнем доверия 95 % вычисляют в соответствии с уравнениями (Е.6) и (Е.7):
Т1(х) = 340,439,
Т2(х) = 344,111,
Т1(y) = 271,203,
Т2(y) = 354,197,
Вычисление смещения:
Таким образом, гипотезу о равенстве двух способов контроля отклоняют. Оценка смещения равна
Тем не менее, дисперсия ошибки исследуемой системы существенно меньше дисперсии ошибки эталонного метода. Рекомендуется провести еще один эксперимент для проверки того, что между двойными измерениями исследуемой системы нет статистической зависимости, например, выполнить измерения в разные дни.
Е.6 Набивка табака в сигареты
Для обеспечения плотности при производстве сигарет проверяют заполнение сигареты табаком. После разработки нового прибора по проверке набивки сигареты было проведено исследование нового и старого приборов. Для этого из партии в 200 кг случайным образом отобрано 200 г исследуемых проб. Исследуемые порции для старого испытательного прибора составляют 15 г, для нового испытательного прибора - 10 г. Результаты двойных измерений старым и новыми приборами приведены в таблице Е.3.
Таблица Е.3 - Результаты измерений набивки сигарет старым и новым приборами
Значения приведены в кубических сантиметрах на грамм
Старый прибор |
Новый прибор |
|||
xi1 |
xi2 |
yi1 |
yi2 |
|
1 |
4,590 |
4,509 |
4,940 |
4,969 |
2 |
4,399 |
4,346 |
4,792 |
4,760 |
3 |
4,700 |
4,700 |
5,259 |
5,312 |
4 |
4,562 |
4,619 |
5,131 |
5,114 |
5 |
4,619 |
4,700 |
4,846 |
4,731 |
6 |
4,399 |
4,481 |
4,771 |
4,763 |
7 |
4,183 |
4,318 |
4,520 |
4,516 |
8 |
4,644 |
4,562 |
4,910 |
4,874 |
9 |
5,214 |
5,189 |
5,988 |
6,168 |
10 |
5,405 |
5,270 |
6,215 |
6,066 |
11 |
5,026 |
4,998 |
5,751 |
5,755 |
12 |
4,916 |
4,916 |
5,581 |
5,600 |
13 |
5,107 |
5,026 |
5,704 |
5,943 |
14 |
4,427 |
4,456 |
4,756 |
4,758 |
15 |
4,807 |
4,782 |
5,738 |
5,591 |
16 |
4,863 |
4,888 |
5,581 |
5,568 |
17 |
5,107 |
5,132 |
5,864 |
5,804 |
18 |
4,835 |
4,888 |
5,630 |
5,572 |
19 |
5,026 |
4,970 |
5,628 |
5,570 |
20 |
5,894 |
5,866 |
6,946 |
6,856 |
Сумма |
96,723 |
96,616 |
108,551 |
108,290 |
Среднее арифметическое |
4,833 |
5,421 |
В соответствии с уравнениями (Е.2) и (Е.4)
se(x) = 0,047,
se(y) = 0,066,
Fa/2(20, 20) = 2,46.
Поскольку Fa/2(20, 20) > F0, то можно считать, что рассматриваемые две группы данных имеют общую дисперсию.
Доверительные границы Т1(х), Т2(х) и Т1(у), Т2(у) для уровня доверия 95 % вычисляют в соответствии с уравнениями (Е.6) и (Е.7)
Т1(х) = 4,735,
Т2(х) = 4,931,
Т1(y) = 5,283,
Т2(y) = 5,559.
Вычисляют
Поскольку гипотезу о равенстве двух способов контроля отклоняют. Оценка смещения отрицательна и равна (-0,588).
Е.7 COD сточных вод
Степень насыщения кислородом COD*) воды, спускаемой предприятием, проверяют двумя методами в зависимости от расположения предприятия. Один метод контролирует количество вещества, связанного с COD. Другой метод контролирует массовую концентрацию COD в сточной воде rCOD. Исследуемый химический завод контролируют с применением второго метода.
_____________
*) COD - Chemical oxygen demand (параметр качества воды).
Пробы отбирают из резервуара для хранения сточных вод химического завода перед разбавлением их чистой водой. Затем выполняют измерения COD. Степень разведения сточных вод чистой водой определяют на основе полученных значений COD. Для выполнения стандартного метода необходимо 2,5 ч, тогда как для быстрого метода требуется не более 20 мин.
Для сравнения прецизионности и смещения быстрого и стандартного методов, определяют COD в двойных выборках быстрым и стандартным методом. Вторые определения быстрым методом двойных данных получены в другие дни для обеспечения независимости данных. Полученные данные приведены в таблице Е.4.
Таблица Е.4 - COD сточных вод до разбавления
Значения в миллиграммах кислорода на литр
Исследуемая система |
Эталонный метод |
|||
xi1 |
xi2 |
yi1 |
yi2 |
|
2 |
4050 |
4300 |
3902 |
4513 |
3 |
3525 |
3925 |
3576 |
4028 |
4 |
2575 |
3400 |
3004 |
3491 |
5 |
4300 |
4450 |
4794 |
5004 |
6 |
3250 |
1825 |
3357 |
3174 |
7 |
2875 |
2215 |
2699 |
2656 |
8 |
3725 |
2975 |
3634 |
3482 |
9 |
3750 |
4325 |
4380 |
4500 |
10 |
201 |
229 |
240 |
229 |
11 |
37 |
33 |
45 |
42 |
12 |
655 |
655 |
1928 |
1838 |
13 |
1788 |
1650 |
2001 |
1987 |
14 |
1675 |
1215 |
1555 |
1476 |
15 |
2050 |
2050 |
2308 |
2384 |
16 |
1180 |
1200 |
1764 |
1740 |
17 |
2325 |
1750 |
2543 |
2385 |
18 |
2100 |
2588 |
3022 |
3040 |
19 |
2650 |
2913 |
2815 |
2823 |
20 |
3050 |
3163 |
3526 |
3230 |
21 |
2488 |
2738 |
2781 |
2554 |
22 |
1340 |
1345 |
1517 |
1554 |
23 |
2575 |
2113 |
2296 |
1874 |
24 |
1800 |
1520 |
1802 |
1838 |
25 |
2400 |
2400 |
2895 |
2745 |
26 |
3150 |
3488 |
3708 |
3945 |
27 |
2875 |
2463 |
3021 |
2963 |
28 |
2413 |
2388 |
2570 |
2555 |
29 |
1925 |
1463 |
2065 |
1948 |
30 |
1975 |
2488 |
2795 |
2861 |
Среднее арифметическое |
2395,9 |
2702,7 |
В соответствии с уравнениями (Е.2) и (Е.4):
se(x) = 327,6;
se(y) = 156,4;
и Fa/2(30, 30) = 2,07.
Поскольку Fa/2(30, 30) < F0, данные нельзя считать принадлежащими к одной совокупности с общей дисперсией.
Доверительные границы Т1(х), Т2(х) и Т1(у), Т2(у) для уровня доверия 95 % вычисляют в соответствии с уравнениями (Е.6) и (Е.7):
Т1(х) = 1726,9;
Т2(х) = 3064,9;
Т1(y) = 2383,3;
Т2(y) = 3022,1.
Вычисляют смещение
Поскольку гипотезу о равенстве двух способов контроля отклоняют. Оценка смещения отрицательна и равна (-306,8).
Задача корректировки недостаточной прецизионности быстрого метода должна решаться за счет увеличения количества выполняемых измерений на партию, но схему установки более точного и дорогостоящего оборудования, обеспечивающего выполнение измерений альтернативным методом, следует рассмотреть.
Е.8 Сравнение портативного и стационарного эмиссионного спектрофотометров
Портативный эмиссионный спектрофотометр предназначен для быстрого анализа химического состава алюминиевых слитков на основе прямой эмиссии с поверхности слитка. Принцип работы портативного спектрофотометра тот же, что и у стационарного спектрофотометра.
Анализируемый портативный спектрофотометр исследуют, используя метод, установленный в ИСО 3301 [19]. При этом анализируют Si и Fe, но для исследований важен только Si. Результаты приведены в таблице Е.5, где xi1 и уi1 - результаты анализа i-й пробы портативным и стационарным спектрофотометрами соответственно.
Таблица Е.5 - Результаты определения содержания Si (%) с использованием портативного и стационарного спектрофотометров
xi1 |
yi1 |
di1 = xi1 - yi2 |
|
|
1 |
1,117 |
0,895 |
0,222 |
0,049284 |
2 |
3,270 |
2,418 |
0,852 |
0,725904 |
3 |
0,720 |
0,620 |
0,100 |
0,010000 |
4 |
2,563 |
2,026 |
0,537 |
0,288369 |
5 |
1,087 |
0,922 |
0,165 |
0,027225 |
6 |
2,627 |
1,898 |
0,729 |
0,531441 |
7 |
0,933 |
0,783 |
0,150 |
0,022500 |
8 |
1,977 |
1,518 |
0,459 |
0,210681 |
9 |
0,677 |
0,627 |
0,050 |
0,002500 |
10 |
2,060 |
1,677 |
0,383 |
0,146689 |
11 |
0,453 |
0,487 |
-0,034 |
0,001156 |
12 |
1,613 |
1,410 |
0,203 |
0,041209 |
13 |
1,097 |
0,965 |
0,132 |
0,017424 |
14 |
0,493 |
0,470 |
0,023 |
0,000529 |
15 |
2,187 |
1,688 |
0,499 |
0,249001 |
16 |
0,567 |
0,482 |
0,085 |
0,007225 |
17 |
2,023 |
1,656 |
0,367 |
0,134689 |
18 |
0,457 |
0,386 |
0,071 |
0,005041 |
19 |
1,643 |
1,491 |
0,152 |
0,023104 |
20 |
1,073 |
0,946 |
0,127 |
0,016129 |
21 |
1,120 |
0,925 |
0,195 |
0,038025 |
Сумма |
29,757 |
24,290 |
5,467 |
2,548125 |
Среднее арифметическое |
1,4170 |
1,1567 |
0,2603 |
- |
Поскольку t1-a/2(v) = t0,975(20) = 2,086, нулевую гипотезу: отклоняют. Между портативным и стационарным спектрофотометрами существует большое смещение. Однако ИСО 3301 [19] лишь описывает сравнение двух методов и содержит рекомендации по проверке их смещения.
Через месяц был проведен эксперимент по исследованию промежуточной прецизионности, в котором были выполнены измерения на одних и тех же пробах с интервалом в один месяц, но пробы, использованные в вышеупомянутом эксперименте, не были сохранены. Тогда был разработан новый эксперимент, предусматривающий выполнение двух последовательных измерений с разницей в один день на одной и той же пробе портативным и стационарным спектрофотометрами. Результаты измерений приведены в таблице Е.6, где хi2 и хi3 - результаты измерений портативным спектрофотометром, а уi2 и уi3 - результаты измерений стационарным спектрофотометром.
Таблица Е.6 - Результаты двойных измерений Si (%) с применением портативного и стационарного спектрофотометров
xi2 |
xi3 |
yi2 |
yi3 |
|
1 |
1,160 |
1,085 |
1,041 |
1,059 |
2 |
1,525 |
1,450 |
1,487 |
1,455 |
3 |
0,340 |
0,315 |
0,342 |
0,334 |
4 |
1,200 |
1,100 |
1,132 |
1,067 |
5 |
0,455 |
0,405 |
0,435 |
0,428 |
6 |
1,820 |
1,730 |
1,657 |
1,611 |
7 |
0,405 |
0,355 |
0,382 |
0,381 |
8 |
2,030 |
1,820 |
1,750 |
1,717 |
9 |
0,465 |
0,410 |
0,435 |
0,432 |
10 |
2,130 |
1,925 |
1,807 |
1,812 |
11 |
1,205 |
1,290 |
1,126 |
1,112 |
12 |
0,915 |
0,850 |
0,909 |
0,897 |
13 |
1,915 |
1,630 |
1,633 |
1,640 |
14 |
2,825 |
2,610 |
2,410 |
2,317 |
15 |
1,955 |
1,930 |
1,680 |
1,726 |
16 |
0,755 |
0,740 |
0,724 |
0,726 |
17 |
2,580 |
2,200 |
1,682 |
1,653 |
18 |
1,945 |
1,895 |
1,743 |
1,774 |
19 |
1,450 |
1,320 |
1,275 |
1,298 |
20 |
1,895 |
1,815 |
2,081 |
2,030 |
21 |
1,990 |
1,710 |
1,630 |
1,686 |
Сумма |
30,960 |
28,585 |
27,361 |
27,155 |
Среднее арифметическое |
1,4177 |
1,2980 |
В соответствии с (Е.2) и (Е.4):
se(x) = 0,1105;
se(y) = 0,0259;
Fa/2(40, 40) = 1,88.
Поскольку Fa/2 < F0, две группы данных нельзя считать принадлежащими к одной совокупности с общей дисперсией.
Доверительные границы Т1(х), Т2(х) и Т1(у), Т2(у) с уровнем доверия 95 % вычислены в соответствии с (Е.6) и (Е.7)
Т1(х) = 1,1877;
Т2(х) = 1,6477;
Т1(y) = 1,2441;
Т2(y) = 1,3519.
Вычисляют
Поскольку гипотезу о равенстве двух способов исследований отклоняют. Оценка смещения равна
Т.к. доверительные интервалы Т1(х), Т2(х) и Т1(у), Т2(у) перекрываются, остается неясным, существует ли смещение между портативным и стационарным спектрофотометрами.
Для улучшения прецизионности портативного спектрофотометра должны быть предприняты специальные действия, например, направленные на снижение шероховатости поверхности исследуемого слитка.
Таблица F.1
Обозначение ссылочного международного стандарта |
Степень соответствия |
Обозначение и наименование соответствующего национального стандарта |
ИСО 565:1990 |
- |
* |
ИСО 3534-1:2006 |
IDT |
* |
ИСО 3534-2:2006 |
IDT |
* |
ISO 3534-3:1999 |
IDT |
Р 50.1.040-2002 Статистические методы. Планирование экспериментов. Термины и определения |
ИСО 5725-1:1994 |
IDT |
ГОСТ Р ИСО 5725-1-2002 Точность (правильность и прецизионность) методов и результатов измерений. Часть 1. Основные положения и определения |
ИСО 5725-2:1994 |
IDT |
ГОСТ Р ИСО 5725-2-2002 Точность (правильность и прецизионность) методов и результатов измерений. Часть 2. Основной метод определения повторяемости и воспроизводимости стандартного метода измерений |
ИСО 5725-3:1994 |
IDT |
ГОСТ Р ИСО 5725-3-2002 Точность (правильность и прецизионность) методов и результатов измерений. Часть 3. Промежуточные показатели прецизионности стандартного метода измерений |
ИСО 5725-4:1994 |
IDT |
ГОСТ Р ИСО 5725-4-2002 Точность (правильность и прецизионность) методов и результатов измерений. Часть 4. Основные методы определения правильности стандартного метода измерений |
ИСО 5725-5:1998 |
IDT |
ГОСТ Р ИСО 5725-5-2002 Точность (правильность и прецизионность) методов и результатов измерений. Часть 5. Альтернативные методы определения прецизионности стандартного метода измерений |
ИСО 5725-6:1994 |
IDT |
ГОСТ Р ИСО 5725-6-2002 Точность (правильность и прецизионность) методов и результатов измерений. Часть 6. Использование значений точности на практике |
* Соответствующий национальный стандарт отсутствует. До его утверждения рекомендуется использовать перевод на русский язык данного международного стандарта. Перевод данного международного стандарта находится в Федеральном информационном фонде технических регламентов и стандартов. Примечание - В настоящей таблице использованы следующие условные обозначения степени соответствия стандартов: IDT - идентичные стандарты |
Cochran, W.G. The relative accuracy of systematic and stratified random samples from a certain class of populations. Annals of Mathematical Statistics, 17, 1946, p. 164
Cochran, W.G. Sampling Techniques. John Wiley, New York, 1953
Aoki, Shigeo. Sample division for size determination. Reports of statistical application research, JUSE, 29 (4), 1982
Aoki, Shigeo and Yoneda, T. Sample division of bulk materials. Reports of statistical application research, JUSE, 34 (4), 1987
Satterthwaite, RE. Biometrics. Bull. 2, 1946
Anderson, R.L. and Bancroft, ТА. Statistical theory in research, McGraw-Hill Book Co., Inc., New York, 1952
Confidence limit of variance component, in JUSE Statistical Table B, 1982 ed.
Jowett, G.H. The accuracy of systematic sampling from conveyor belts. Applied Statistics, 1952
Hebden, J. and Jowett, G.H. The accuracy of sampling coal, Applied Statistics, 1952
Jowett, G.H. and Scott, J.F. Simple graphical techniques for calculating serial and spatial correlations, and mean semi-squared differences. J.R. Statist. Soc, B, 15, 1953, p. 81
Jowett, G.H. The comparison of means of industrial time series. Applied Statistics, 1953 Ward, D.H. Weekly, monthly, and quarterly tolerances for coke quality. Applied Statistics, 1959
GY, P.M. Sampling of Particulate Materials - Theory and Practice. Elsevier Scientific Publishing Co., Amsterdam, 2nd ed., 1982
GY, P.M. Sampling of heterogeneous and dynamic material systems - Theories of heterogeneity, sampling and homogenizing. Elsevier Scientific Publishing Co., Amsterdam, 1992
GY, P.M. L¢echantillonnage des lots de matiere en vue de leur analyse. Masson publ., 1996 ISO 3084, Iron ores - Experimental methods for evaluation of quality variation
ISO 3085, Iron ores - Experimental methods for checking the precision of sampling, sample preparation and measurement (ИСО 3085:2002 Руды железные. Экспериментальные методы контроля точности отбора проб, приготовления образцов и измерения)*
ISO 3086, Iron ores - Experimental methods for checking the bias of sampling
ISO 3301, Statistical interpretation of data - Comparison of two means in the case of paired observations (ИСО 3301:1975 Статистическое представление данных. Сравнение двух средних в парных наблюдениях)*
Snedecor, G.W and Cochran, W.G. Statistical Methods. IOWA State University Press, 7th ed., 1980
Aoki, Shigeo. Case studies for bias testing. Bulletin for sampling research JUSE, 237, February 1991
Grubbs, F.E. On estimating precision of measuring instruments and product variability. J.A.S.A., 1948
Thomson, G.J. JR. Precision of simultaneous measurement procedures, J.A.S.A., 1963
Hahn, G.J. and Nelson, W.A. Problem in the statistical comparison of measuring devices. Technometrics, 1970
Maloney, C.J. and Rastogi, S.C. Significance test for Grubbs estimators. Biometrics, 1970
Jaech, John L. Further tests of significance for Grubbs estimators. Biometrics, December 1971
Grubbs, F.E. Errors of measurement, precision, accuracy and the statistical comparison of measuring instruments. Technometrics, 1973
ISO 2854 Statistical interpretation of data - Techniques of estimation and tests relating to means and variances (ИСО 2854 Статистическое представление данных. Методы оценки и проверки гипотез о средних значениях и дисперсиях)*
____________
* Официальный перевод этого стандарта находится в Федеральном информационном фонде.
ISO 8258 Shewhart control charts
ISO 10725 Acceptance sampling plans and procedures for the inspection of bulk materials
Guide to the expression of uncertainty in measurement (GUM), BIPM, IEC, IFCC, ISO, IUPAC, IUPAP, OIML, 1st ed., 1993
Ключевые слова: статистический приемочный контроль, план выборочного контроля, контроль по альтернативному признаку, выборка, партия, приемлемый уровень качества, единица продукции, несоответствие, несоответствующая единица продукции