| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
межгосударственный совет по стандартизации, метрологии и сертификации (МГС) INTERSTATE COUNCIL FOR STANDARDIZATION, METROLOGY AND CERTIFICATION (ISC)
Государственная система обеспечения ОЦЕНИВАНИЕ
МЕТРОЛОГИЧЕСКИХ
Цели, основные принципы и порядок проведения работ по межгосударственной стандартизации установлены ГОСТ 1.0-92 «Межгосударственная система стандартизации. Основные положения» и ГОСТ 1.2-2009 «Межгосударственная система стандартизации. Стандарты межгосударственные, правила и рекомендации по межгосударственной стандартизации. Правила разработки, принятия, применения, обновления и отмены» Сведения о рекомендациях 1 РАЗРАБОТАНЫ Федеральным государственным унитарным предприятием «Уральский научно-исследовательский институт метрологии» (ФГУП «УНИИМ») 2 ВНЕСЕНЫ Управлением метрологии Федерального агентства по техническому регулированию и метрологии 3 ПРИНЯТЫ Межгосударственным советом по стандартизации, метрологии и сертификации (протокол от 11 ноября 2009 г. № 36-2009) За принятие проголосовали:
4 Приказом Федерального агентства по техническому регулированию и метрологии от 23 декабря 2010 г. № 1001-ст рекомендации по межгосударственной стандартизации РМГ 93-2009 введены в действие в качестве рекомендаций по метрологии Российской Федерации с 1 января 2012 г. 5 ВВЕДЕНЫ ВПЕРВЫЕ Информация об изменениях к настоящим рекомендациям публикуется в ежегодно издаваемом указателе «Руководящие документы, рекомендации и правила», а текст изменений и поправок - в ежемесячно издаваемых информационных указателях «Национальные стандарты». В случае пересмотра (замены) или отмены настоящих рекомендаций соответствующее уведомление будет опубликовано в ежемесячно издаваемом информационном указателе «Национальные стандарты». Соответствующая информация, уведомление и тексты размещаются также в информационной системе общего пользования - на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет СОДЕРЖАНИЕ РЕКОМЕНДАЦИИ ПО МЕЖГОСУДАРСТВЕННОЙ СТАНДАРТИЗАЦИИ Государственная система обеспечения единства измерений ОЦЕНИВАНИЕ МЕТРОЛОГИЧЕСКИХ ХАРАКТЕРИСТИК СТАНДАРТНЫХ ОБРАЗЦОВ State system for ensuring the uniformity of measurements. Estimation of metrological characteristics of reference materials Дата введения - 2012-01-01 1 Область примененияНастоящие рекомендации распространяются на стандартные образцы состава и свойств веществ и материалов (далее - СО) и устанавливают алгоритмы оценивания их метрологических характеристик. Номенклатура метрологических характеристик СО, форма их представления и способы нормирования приведены в приложении А. 2 Нормативные ссылкиВ настоящих рекомендациях использованы нормативные ссылки на следующие нормативные документы: ГОСТ 8.010-90* Государственная система обеспечения единства измерений. Методики выполнения измерений ______________ * Утратил силу на территории РФ. Действует ГОСТ Р 8.563-2009. ГОСТ 8.315-97 Государственная система обеспечения единства измерений. Стандартные образцы состава и свойств веществ и материалов. Основные положения ГОСТ 8.531-2002 Государственная система обеспечения единства измерений. Стандартные образцы состава монолитных и дисперсных материалов. Способы оценивания однородности ГОСТ 8.532-2002 Государственная система обеспечения единства измерений. Стандартные образцы состава веществ и материалов. Межлабораторная метрологическая аттестация. Содержание и порядок проведения работ РМГ 53-2002 Государственная система обеспечения единства измерений. Стандартные образцы. Оценивание метрологических характеристик с использованием эталонов и образцовых средств измерений Примечание - При пользовании настоящими рекомендациями целесообразно проверить действие ссылочных межгосударственных документов на территории государства по соответствующему указателю стандартов, составленному по состоянию на 1 января текущего года, и по соответствующим информационным указателям, опубликованным в текущем году. Если ссылочный документ заменен (изменен), то при пользовании настоящими рекомендациями, следует руководствоваться заменяющим (измененным) документом. Если ссылочный документ отменен без замены, то положение, в котором дана ссылка на него, применяется в части, не затрагивающей эту ссылку. 3 Термины, определения и обозначения3.1 В настоящих рекомендациях применены термины по ГОСТ 8.315, а также следующие термины с соответствующими определениями: 3.1.1
3.1.2
3.1.3
3.1.4
3.1.5
Примечания 1 Расширенную неопределенность U вычисляют исходя из суммарной стандартной неопределенности ис и коэффициента охвата k по формуле U = kuc. 2 Выбор коэффициента k зависит от требуемого уровня достоверности. При доверительной вероятности Р = 0,95 k равно 2. 3 Коэффициент охвата k необходимо указывать для того, чтобы можно было восстановить значение суммарной стандартной неопределенности измеряемой величины для использования в вычислениях суммарной стандартной неопределенности других величин, которые могут зависеть от этой величины. 3.1.6 оценивание типа А: Способ оценивания стандартной неопределенности, основанный на статистическом анализе серии наблюдений. Примечание - Оценивание типа А стандартной неопределенности проводят по Руководству [1], подраздел 4.2. 3.1.7 оценивание типа В: Способ оценивания стандартной неопределенности, не основанный на статистическом анализе серии наблюдений. Примечания 1 Оценивание типа В стандартной неопределенности проводят по Руководству [1], подраздел 4.3. 2 В данных рекомендациях при оценивании стандартной неопределенности величины типа В использованы только предположения относительно функции распределения этой величины. 3.1.8 способ определения аттестованного значения стандартного образца: Метрологически обоснованная процедура установления аттестованного значения СО. 3.1.9 неопределенность от способа определения аттестованного значения стандартного образца: Составляющая неопределенности аттестованного значения СО, обусловленная способом определения аттестованного значения стандартного образца. Примечание - В Руководстве [2] эта составляющая неопределенности трактуется как «неопределенность от характеризации» («uncertainty due to characterisation»). 3.1.10 стандартная неопределенность от способа определения аттестованного значения стандартного образца: Неопределенность от способа определения аттестованного значения стандартного образца, выраженная в виде стандартного отклонения. 3.1.11 неопределенность от неоднородности: Составляющая неопределенности аттестованного значения СО, обусловленная различием значений аттестуемой характеристики СО в разных частях (пробах) материала СО. 3.1.12 стандартная неопределенность от неоднородности: Неопределенность от неоднородности, выраженная в виде стандартного отклонения. 3.1.13 неопределенность от нестабильности: Составляющая неопределенности аттестованного значения СО, обусловленная изменением значения аттестуемой характеристики СО в течение срока годности экземпляра СО. 3.1.14 стандартная неопределенность от нестабильности: Неопределенность от нестабильности, выраженная в виде стандартного отклонения. 3.1.15 влияющие величины: Величины, определяющие значение аттестуемой характеристики СО в использованной модели измерения, в число которых входят поправки, а также все иные величины, которые отражают другие источники изменчивости, вносящие существенную неопределенность в результат измерения. 3.1.16 неопределенность аттестованного значения стандартного образца: Параметр, характеризующий рассеяние значений, которые могут быть обоснованно приписаны аттестуемой характеристике СО. 3.1.17 стандартная неопределенность аттестованного значения стандартного образца: Неопределенность аттестованного значения СО, выраженная в виде стандартного отклонения. 3.1.18 суммарная стандартная неопределенность аттестованного значения стандартного образца: Стандартная неопределенность аттестованного значения СО, оцененная с учетом вклада стандартных неопределенностей от способа определения аттестованного значения стандартного образца, от неоднородности, от нестабильности. 3.1.19 расширенная неопределенность аттестованного значения стандартного образца: Неопределенность аттестованного значения СО, представляющая собой величину, определяющую интервал вокруг аттестованного значения СО, в пределах которого, как можно ожидать, находится большая часть распределения значений, которые с достаточным основанием могут быть приписаны аттестуемой характеристике СО. 3.1.20 неэмпирические (рациональные) методики выполнения измерений: Методики выполнения измерений, результаты которых не зависят от используемого метода. Примечание - В Руководстве [3], подраздел 5.3 использован термин «неэмпирические (рациональные) методы». 3.1.21 эмпирические методики выполнения измерений: Методики выполнения измерений, результаты которых зависят от используемого метода. Примечание - В Руководстве [3], подраздел 5.2 использован термин «эмпирические методы». 3.1.22
3.1.23
3.1.24
3.1.26
3.1.27
3.1.28
3.1.29 стабильность стандартного образца: Свойство материала СО, выражающееся в неизменности значений аттестованной характеристики СО во времени при соблюдении условий хранения и применения. 3.1.30 исследование стабильности стандартного образца: Изучение материала СО с целью установления срока годности экземпляра СО, условий его хранения и применения, при которых оцениваются возможные изменения значений аттестованной характеристики СО под влиянием факторов нестабильности. 3.1.31 факторы нестабильности: Совокупность внешних условий и физических и химических процессов, протекающих в материале СО, вызывающих изменение аттестованной характеристики СО. 3.1.32 ускоренное старение: Хранение материала СО при исследовании стабильности СО в условиях, существенно усиливающих воздействие факторов нестабильности. 3.2 В настоящих рекомендациях использованы следующие обозначения: А - аттестованное значение СО; хi - результат измерения значения аттестуемой характеристики стандартного образца; уi - оценка i-й влияющей величины; и - стандартная неопределенность; uА - стандартная неопределенность при оценивании типа А; иB - стандартная неопределенность при оценивании типа В; ис - суммарная стандартная неопределенность; k - коэффициент охвата; Р - доверительная вероятность; tP(v) - квантиль распределения Стьюдента для доверительной вероятности Р и числа степеней свободы v; р - число лабораторий, участвующих в межлабораторном эксперименте; v - число степеней свободы; Uдоп - допускаемое значение расширенной неопределенности аттестованного значения стандартного образца; Uр - расширенная неопределенность для доверительной вероятности Р; ис(А) - суммарная стандартная неопределенность аттестованного значения СО; UP(A) - расширенная неопределенность аттестованного значения СО для доверительной вероятности Р; uchar - стандартная неопределенность от способа определения аттестованного значения СО; uh - стандартная неопределенность от неоднородности СО; ustab - стандартная неопределенность от нестабильности СО; sА - стандартное отклонение аттестованного значения СО; DА - граница погрешности аттестованного значения СО для доверительной вероятности Р; schar - стандартное отклонение от способа определения аттестованного значения СО; sh - стандартное отклонение от неоднородности; sstab - стандартное отклонение от нестабильности. 4 Метрологические характеристики стандартных образцов4.3 Исходя из положений 4.1 и 4.2 в рекомендации изложены алгоритмы оценивания стандартных неопределенностей аттестованного значения СО. 4.4 При необходимости перехода от стандартных неопределенностей аттестованного значения СО к характеристикам погрешности аттестованного значения СО используют соотношения, приведенные в таблице 4.1. Таблица 4.1 - Соотношения между неопределенностями и характеристиками погрешности аттестованного значения СО
5 Источники неопределенности аттестованного значения стандартного образца5.1 Основными источниками неопределенности аттестованного значения СО являются: - неоднородность материала СО (далее - неоднородность); - нестабильность значений аттестуемой характеристики СО (далее - нестабильность); - способ определения аттестованного значения СО. 5.2 На начальном этапе проводят анализ возможных влияний каждого источника неопределенности, указанного в 5.1, на значения аттестуемой характеристики СО. В общем случае суммарную стандартную неопределенность аттестованного значения СО определяют по уравнению (5.1) Примечание - В зависимости от типа материала СО и условий хранения экземпляров СО некоторые составляющие суммарной стандартной неопределенности аттестованного значения СО либо равны нулю, либо их вклад в неопределенность аттестованного значения СО является незначимым. 6 Оценивание неопределенности от нестабильности6.1 Способы оценивания неопределенности от нестабильностиОценивание неопределенности от нестабильности осуществляют в зависимости от условий проведения измерений двумя способами: - классическое исследование стабильности согласно Руководству [2]; - изохронное исследование стабильности. 6.2 Классическое исследование стабильности6.2.1 Измерения при исследовании стабильности проводят в условиях промежуточной прецизионности. Факторами, влияющими на расхождения результатов, в данном случае являются следующие: а) время; б) градуировка оборудования; в) условия окружающей среды; г) нестабильность измерительной системы; д) изменение значения аттестуемой характеристики; е) повторяемость. Факторы а) - д) варьируют в процессе исследования стабильности. Рассеяние результатов измерений при воздействии факторов характеризуют стандартным отклонением прецизионности (при различиях по фактору «время») sI(T). 6.2.2 Для оценивания неопределенности от нестабильности СО используют методики выполнения измерений (МВИ), аттестованные в соответствии с требованиями ГОСТ 8.010. Рекомендуется использовать МВИ с оцененным значением стандартного отклонения промежуточной прецизионности (при различиях по фактору «время») sI(T). Стандартное отклонение промежуточной прецизионности sI(T) должно удовлетворять условию (6.1) где Uдоп - допускаемое значение расширенной неопределенности аттестованного значения СО в соответствии с техническим заданием на разработку СО. 6.2.3 Продолжительность исследования стабильности t должна быть более половины предполагаемого срока годности экземпляра СО. Примечание - Продолжительность исследования при ускоренном старении материала СО сокращается. Ее определяют по предполагаемому сроку годности экземпляра СО и известной или оцененной зависимости изменений аттестуемой характеристики от факторов нестабильности. 6.2.4 Для оценивания неопределенности от нестабильности за период исследования стабильности СО t получают п результатов измерений значений аттестуемой характеристики хi (i = 1, 2, ..., п). 6.2.5 Число измерений п определяют по таблице 6.1 с учетом предварительно вычисленного значения отношения Таблица 6.1 - Минимальное число измерений при исследовании стабильности СО
Рекомендуется в период исследования стабильности СО t измерения проводить через равные промежутки времени. 6.2.6 Отклонение результата измерения хi в i-й момент времени от первоначального значения х1 оценивают по формуле di = xi - x1, (6.2) где хi - результат измерения значения аттестуемой характеристики СО в i-й момент времени; х1 - первый результат измерения, полученный в период исследования стабильности СО. 6.2.7 Для уменьшения влияния факторов промежуточной прецизионности проводят экспоненциальное сглаживание полученных значений di по формуле Di = adi + (1 - a)Di-1, (6.3) где Di - сглаженное значение разности результатов измерений в i-й момент времени (i = 2, ..., п). Начальным значением D1 считают d1 = 0. Значение a выбирают из интервала от 0,1 до 0,3 в зависимости от значения отношения по таблице 6.2. Таблица 6.2 - Значение коэффициента a для экспоненциального сглаживания 6.2.8 Результаты контроля стабильности записывают в форме, приведенной в таблице 6.3. Таблица 6.3 - Результаты исследования стабильности СО
6.2.9 По вычисленным значениям Di для i = 2, 3, ..., n определяют скользящие размахи Ri по формуле Ri = |Di - Di-1| (6.4) и вписывают полученные значения в графу «Ri» таблицы 6.3. Вычисляют средний скользящий размах по формуле (6.5) 6.2.10 Предполагают линейную модель изменения аттестуемой характеристики D в момент времени t D = at, (6.6) где а - скорость изменения значений результатов измерений аттестуемой характеристики, полученных в условиях промежуточной прецизионности. Оцененное по результатам измерений изменение аттестуемой характеристики Di в момент времени ti представляют в виде уравнения Di = ati + ei, (6.7) где ei - отклонение i-го результата измерения в условиях промежуточной прецизионности. 6.2.11 Определяют коэффициент а в уравнении (6.7) по полученным значениям Di методом наименьших квадратов по формуле (6.8) 6.2.12 Определяют стандартное отклонение sa коэффициента а по формуле (6.9) где sD - стандартное отклонение разностей результатов измерений, вычисленное по среднему скользящему размаху по формуле sD = 0,89. (6.10) 6.2.13 Стандартную неопределенность от нестабильности ustab в момент времени t оценивают по формуле ustab = sat. (6.11) Число степеней свободы для ustab оценивают по формуле vstab = n - 1. (6.12) 6.2.14 Проверяют гипотезу об отсутствии тренда Н0а = 0. (6.13) Для проверки гипотезы Н0 вычисляют отношение по формуле (6.14) Сравнивают полученное значение с квантилем распределения Стьюдента t0,95(n - 1), указанным в приложении Б, таблице Б.2. При выполнении неравенства (6.15) гипотезу об отсутствии тренда принимают. Статистически значимого изменения за период исследования стабильности не обнаружено. Но вклад неопределенности от нестабильности в суммарную неопределенность в этом случае следует оценивать с учетом назначенного срока годности экземпляра СО по формуле (6.11). 6.3 Изохронное исследование стабильности6.3.1 При изохронном исследовании стабильности применяют «метод ускоренного старения». Суть этого метода заключается в том, что усиливают воздействие факторов, вызывающих тренд аттестованной характеристики за период времени намного меньше срока годности СО. Для обоснования метода ускоренного старения обычно применяют правило Вант-Гоффа для медленных реакций: скорость реакции при нагреве на 10 °С увеличивается в 2 - 4 раза. 6.3.2 Материал СО, отобранный для исследования стабильности, разделяют на две части. Одну из этих частей хранят в предполагаемых условиях хранения СО, другую - при повышенной температуре. При фиксированных значениях температуры хранения продолжительность исследования стабильности t оценивают по формуле (6.16) где Т - предполагаемый срок годности экземпляра СО; t0, t1 - температура хранения материала СО и температура хранения СО при ускоренном старении. 6.3.3 В течение периода времени t получают п пар результатов измерений {х0i, х1i} в условиях повторяемости. Рекомендуется проводить измерения в период исследования стабильности СО t через равные промежутки времени. Вычисляют отклонение результатов измерений di в i-й момент времени di = x1i - x0i. (6.17) 6.3.4 После проведения всех измерений оценивают стандартное отклонение повторяемости sr по формуле (6.18) 6.3.5 Предполагают (аналогично классическому методу исследования стабильности по 6.2.10), что разность результатов di в момент времени ti можно представить в виде линейной зависимости di = ati + ei. (6.19) 6.3.6 Определяют коэффициент а в уравнении (6.19) по полученным значениям di методом наименьших квадратов по формуле (6.20) 6.3.7 Определяют стандартное отклонение sa коэффициента а по формуле (6.21) 6.3.8 Стандартную неопределенность от нестабильности ustab в момент времени t и число степеней свободы vstab оценивают по формулам (6.11) и (6.12). 6.3.9 Гипотезу об отсутствии тренда Н0 проверяют в соответствии с требованиями 6.2.14. 7 Оценивание неопределенности от неоднородности материала стандартного образца7.1 Общие требования к оцениванию неоднородности материала стандартного образцаОбщие требования к оцениванию неоднородности материала СО приведены в ГОСТ 8.531. 7.2 Оценивание неопределенности от неоднородности для дисперсных материалов7.2.1 Для оценивания неоднородности материала СО от всей массы материала СО случайным образом отбирают N проб массой М0 каждая. Отбор проб проводят после приготовления материала СО. Масса каждой пробы М0 должна быть достаточной для проведения измерений в соответствии с применяемой МВИ. Для определения числа отбираемых проб рассчитывают отношение Q по формуле (7.1) где Uдоп - допускаемое значение расширенной неопределенности аттестованного значения стандартного образца; SMBИ - стандартное отклонение повторяемости или промежуточной прецизионности методики выполнения измерений по [5]. 7.2.2 Число отбираемых проб N при фиксированном числе многократных измерений J находят по таблице 7.1 для значения Q, определенного в соответствии с требованиями 7.1 Таблица 7.1 - Число отбираемых проб для оценивания однородности дисперсных материалов
7.2.3 В каждой из N проб J раз измеряют значение аттестуемой характеристики. Измерения выполняют либо в одной пробе массой М0 неразрушающим методом, либо в растворе пробы, обеспечивающим ее однородность. 7.2.4 Результаты измерений при оценивании характеристики однородности обрабатывают в следующем порядке. Вычисляют средние арифметические значения всех NJ результатов измерений (7.2) где Xnj - j-й результат измерений в n-й пробе и N средних результатов измерений для каждой пробы (7.3) Вычисляют сумму квадратов отклонений SSe результатов измерений от средних значений каждой пробы (7.4) и сумму квадратов отклонений SSн средних результатов для каждой пробы от среднего арифметического значения всех результатов измерений (7.5) Вычисляют средние квадраты отклонений результатов измерений от средних значений для каждой пробы (7.6) и между пробами (7.7) Стандартную неопределенность от неоднородности оценивают по формуле Если , то (7.9) где М - масса наименьшей представительной пробы СО. Число степеней свободы стандартной неопределенности uh, обозначаемое vh, равно N - 1. 7.3 Оценивание неопределенности от неоднородности для монолитных материалов7.3.1 Оценивание неопределенности от неоднородности проводят после отработки технологии получения материала СО, исключающей регулярные изменения содержаний аттестуемого элемента, порядка приготовления материала СО и разделения его на экземпляры. 7.3.2 Из общего числа экземпляров СО случайным образом отбирают K экземпляров СО (K ³ 25). 7.3.3 Готовят на каждом отобранном экземпляре СО анализируемые поверхности в соответствии с методикой измерений, используемой для оценивания однородности. 7.3.5 После проведения измерений разрезают каждый экземпляр СО по плоскости, параллельной проанализированной поверхности. Положение плоскости разреза на каждом экземпляре СО определяют случайным образом на всей длине (высоте). Подготавливают на срезах анализируемые поверхности и проводят измерения в соответствии с требованиями 7.3.4. 7.3.6 По результатам измерений для каждого аттестуемого элемента вычисляют следующие суммы: - сумму всех результатов S1 - сумму квадратов средних значений результатов измерений, полученных при анализе аналитических поверхностей S2 (7.11) - сумму квадратов средних значений результатов измерений, полученных при анализе экземпляров СО S3 (7.12) - сумму квадратов всех результатов S4 где i - номер экземпляра СО (i = 1, 2, ..., K); j - номер анализируемой поверхности (j = 1, 2); п - номер измерения (п = 1, 2); хijn - результат n-го измерения на j-й поверхности i-го СО. По суммам, вычисленным по формулам (7.10) - (7.13), вычисляют следующие суммы квадратов отклонений: - средних значений результатов измерений, полученных при анализе экземпляров СО, от среднего значения всех результатов измерений SSBL (7.14) - средних значений результатов измерений, полученных при анализе аналитических поверхностей, от средних значений результатов измерений, полученных при анализе экземпляров СО SSBB SSBB = S2 - S3, (7.15) - результатов измерений от средних значений результатов измерений, полученных при анализе аналитических поверхностей SSW SSW = S4 - S2, (7.16) - результатов измерений от среднего значения всех результатов SST (7.17) Для контроля правильности вычислений проверяют соотношение между суммами квадратов. Если вычисления проведены правильно, то должно быть выполнено равенство SSBL + SSBB + SSW = SST. (7.18) В том случае, если суммы квадратов удовлетворяют уравнению (7.18), вычисляют средние квадраты по соответствующим суммам квадратов отклонений, вычисленных по формулам (7.19) - (7.21): (7.20) и стандартное отклонение повторяемости результатов измерений (7.22) Значение SM характеризует неопределенность интегрального (рентгенофлуоресцентного) метода анализа в условиях повторяемости. При оценивании однородности локальным (оптико-эмиссионным) методом SM характеризует суммарную неопределенность повторяемости метода и различие содержания аттестуемого элемента в аналитических объемах. 7.3.7 По средним квадратам, вычисленным по формулам (7.19) - (7.21), оценивают дисперсию внутриэкземплярной неоднородности и дисперсию межэкземплярной неоднородности по формулам (7.23) (7.24) 7.3.8 Оценки стандартных отклонений Sмак, характеризующего межэкземплярную неоднородность, и Sмик, характеризующего внутриэкземплярную неоднородность, проводят в зависимости от соотношений между средними квадратами MSW, MSBB и MSBL по формулам, приведенным в таблице 7.2. В таблице 7.2 буквой п обозначено число измерений для воспроизведения аттестованного значения СО локальным (оптико-эмиссионным) методом. Таблица 7.2 - Оценки стандартных отклонений Sмак и Sмик при различных соотношениях между средними квадратами MSW, MSBB, MSBL
7.3.9 Вычисление стандартной неопределенности от неоднородности проводят по формуле (7.25) Число степеней свободы vh стандартной неопределенности uh в зависимости от соотношений между средними квадратами MSW, MSBB и MSBL определяют по таблице 7.3. Таблица 7.3 - Число степеней свободы при различных соотношениях между средними квадратами MSW, MSBB, MSBL
8 Оценивание неопределенности от способа определения аттестованного значения стандартного образца8.1 Способы определения аттестованного значения стандартного образцаАлгоритм оценивания неопределенности от способа определения аттестованного значения СО зависит от способа определения аттестованного значения СО. Существуют следующие способы определения аттестованного значения СО по ГОСТ 8.315: по результатам, полученным в нескольких лабораториях одним или несколькими методами; по результатам, полученным в одной лаборатории несколькими методами: - с использованием эталонов; - с использованием МВИ, аттестованной в соответствии с требованиями ГОСТ 8.010; - по процедуре приготовления. 8.2 Оценивание неопределенности от способа определения аттестованного значения стандартного образца при межлабораторном эксперименте8.2.1 Виды методик выполнения измерений Методики выполнения измерений, применяемые для измерения значения аттестуемой характеристики стандартного образца, разделяют на рациональные и эмпирические в соответствии с Руководством [2]. Измерения значения аттестуемой характеристики стандартного образца при проведении межлабораторного эксперимента обычно проводят либо с использованием нескольких рациональных МВИ, либо с использованием одной и той же эмпирической МВИ во всех лабораториях. Для оценивания неопределенности от способа определения аттестованного значения стандартного образца, предусматривающего применение эмпирических или рациональных МВИ, используют разные алгоритмы. 8.2.2 Оценивание неопределенности от способа определения аттестованного значения при применении эмпирических методик 8.2.2.1 Применяемая МВИ должна иметь характеристики прецизионности (стандартное отклонение повторяемости sr и стандартное отклонение воспроизводимости sR), установленные по [6]. Примечание - При отсутствии характеристик прецизионности их оценки для использования в дальнейшем могут быть получены по результатам межлабораторного эксперимента по алгоритму, изложенному в [6]. 8.2.2.2 Межлабораторный эксперимент проводят с привлечением лабораторий, имеющих опыт использования применяемой МВИ. Оптимальное число лабораторий р для получения достоверной оценки неопределенности от способа определения аттестованного значения должно быть от 10 до 15. 8.2.2.3 Каждая i-я лаборатория выполняет по МВИ в условиях повторяемости п измерений (число измерений п £ 4). Протокол с результатами всех измерений (, ..., хin) в условиях повторяемости лаборатория направляет разработчику СО. 8.2.2.4 Разработчик СО оценивает приемлемость результатов по [7]. Вычисляют для каждой i-й лаборатории диапазон значений результатов измерений Qi, по формуле Qi = ximax - ximin, (8.1) где ximax и ximin - наибольший и наименьший результаты измерений в i-й лаборатории. Сравнивают диапазон Qi с критическим диапазоном CR0,95(n) CR0,95(n) = f(n)sr, (8.2) где f(2) = 2,8; f(3) = 3,3; f(4) = 3,6. При выполнении неравенства Qi £ CR0,95(n) (8.3) результаты i-й лаборатории принимают для дальнейшей обработки. В противном случае результаты i-й лаборатории исключают и предлагают данной лаборатории провести все измерения повторно. 8.2.2.5 Все результаты измерений, оставленные после отбраковки по 8.2.2.4, обрабатывают следующим образом. Вычисляют: - по результатам измерений каждой i-й лаборатории среднее значение (8.4) - общее среднее значение по всем результатам - стандартное отклонение повторяемости в i-й лаборатории si (8.6) - стандартное отклонение повторяемости sr - межлабораторную дисперсию В случае, если по уравнению (8.8) получено отрицательное значение , то принимают = 0. 8.2.2.6 По стандартным отклонениям воспроизводимости sR и повторяемости sr, приписанным данной МВИ, определяют межлабораторную дисперсию (8.9) 8.2.2.7 Сравнивают оценку межлабораторной дисперсии с дисперсией по критерию хи-квадрат, описываемому неравенством где - квантиль распределения с v = р - 1 степенями свободы, значения которого приведены в таблице Б.1 приложения Б. Если неравенство (8.10) выполнено, то межлабораторная дисперсия является приемлемой. Результаты измерений всех лабораторий в этом случае можно использовать для оценивания аттестованного значения СО и стандартной неопределенности от способа определения аттестованного значения СО в соответствии с требованиями 8.2.2.8. Если неравенство (8.10) не выполнено, то межлабораторная дисперсия является неприемлемой вследствие плохой воспроизводимости МВИ в одной или нескольких лабораториях. Оценивание аттестованного значения СО и стандартной неопределенности uchar в этом случае проводят в соответствии с требованиями 8.2.2.9. 8.2.2.8 Аттестованное значение СО оценивают как среднее арифметическое значение всех результатов, вычисленных по формуле (8.5) (8.11) Стандартную неопределенность от способа определения аттестованного значения СО uchar вычисляют по формуле (8.12) 8.2.2.9 При невыполнении условия (8.10) аттестованное значение СО вычисляют как среднее взвешенное средних результатов измерений, полученных в лабораториях, по формуле (8.13) где W - сумма весовых коэффициентов wi. Весовые коэффициенты вычисляют следующим образом: - располагают средние результаты измерений , полученные во всех лабораториях, в ряд по возрастанию £ ... £ £ ... < , (8.14) где (i) - номер результата в ряду по возрастанию, (р) - общее число результатов измерений; - вычисляют медиану результатов по формуле - вычисляют абсолютные отклонения результатов измерений d0i от медианы по формуле - вычисляют медиану абсолютных ненулевых отклонений по формуле MAD0 = med{d0i}, (8.17) где med{d0i} вычисляют в зависимости от числа ненулевых значений d0i по формуле (8.15); - для каждого результата вычисляют нормированное отклонение Ui от медианы по формуле (8.18) где величины d0i и MAD0 вычислены по формулам (8.16) и (8.17); - весовые коэффициенты wi в зависимости от значений Ui вычисляют по формуле (8.19) Для оценивания стандартной неопределенности uchar вычисляют абсолютные отклонения результатов измерений от среднего взвешенного значения d2i (8.20) и медиану абсолютных ненулевых отклонений MAD2 по формуле MAD2 = med{d2i}, (8.17) где med{d2i} вычисляет в зависимости от числа ненулевых значений d2i по формуле (8.15). Стандартную неопределенность от способа определения аттестованного значения СО uchar вычисляют по формуле uchar = 1,48MAD2. (8.22) Число степеней свободы vchar равно целому значению суммы весовых коэффициентов W. 8.2.3 Оценивание неопределенности от способа определения аттестованного значения СО при применении рациональных методик 8.2.3.1 В межлабораторной аттестации принимают участие р лабораторий, которые применяют q МВИ (q ³ 2), основанных на разных физических и химических принципах. Если лаборатория применяет при межлабораторном эксперименте несколько МВИ, то результаты измерений этой лаборатории по каждой МВИ представляют как результаты измерений отдельной лаборатории. Предполагают, что каждая j-я МВИ имеет оцененные в соответствии с [6] характеристики прецизионности (стандартное отклонение повторяемости srj, стандартное отклонение воспроизводимости sRj и стандартное отклонение систематического смещения метода sс). 8.2.3.2 Разработчик СО оценивает приемлемость результатов в условиях повторяемости по процедуре, приведенной в 8.2.2.3. 8.2.3.3 Для каждой j-й МВИ вычисляют: - среднее значение , полученное в i-й лаборатории по j-й МВИ по формуле (8.23) где xjik - k-й результат измерений, полученный в i-й лаборатории по j-й МВИ; - общее среднее значение всех результатов, полученных по j-й МВИ по формуле (8.24) где рj - число лабораторий, представивших результаты по j-й МВИ; - стандартное отклонение повторяемости sji результатов измерений, полученных в i-й лаборатории по j-й МВИ по формуле (8.25) - оценку стандартного отклонения повторяемости srj результатов измерений, полученных по j-й МВИ по формуле (8.26) - оценку межлабораторной дисперсии результатов измерений, полученных по j-й МВИ по формуле В случае, если по уравнению (8.27) получено отрицательное значение , то принимают равным 0. 8.2.3.4 По характеристикам прецизионности, приписанным j-й МВИ, вычисляют межлабораторную дисперсию (8.28) 8.2.3.5 Сравнивают оценку межлабораторной дисперсии с дисперсией по критерию хи-квадрат, описываемому неравенством в котором - квантиль c2 распределения c v = pj - 1 степенями свободы, значения которого приведены в таблице Б.1 приложения Б. 8.2.3.6 При выполнении неравенства (8.29) межлабораторная дисперсия является приемлемой. Результаты всех лабораторий в этом случае используют для оценивания аттестованного значения СО и стандартной неопределенности способа. Для j-й МВИ стандартную неопределенность среднего значения вычисляют по формуле (8.30) Число степеней свободы стандартной неопределенности равно vj = рj - 1. Если неравенство (8.29) не выполнено, то межлабораторная дисперсия является неприемлемой вследствие недостаточной воспроизводимости результатов измерений, полученных по j-й МВИ в одной или нескольких лабораториях. В этом случае возможны следующие варианты: а) использовать критерий для выявления статистических выбросов (например, критерий Граббса в соответствии с [6]) и повторить процедуру проверки приемлемости результатов по 8.2.2.4; б) забраковать все результаты при их небольшом числе. - располагают все средние значения результатов измерений, полученных по МВИ, по возрастанию £ ... £ ... < (8.31) - вычисляют отношение (8.32) - сравнивают полученное значение с квантилем распределения Стьюдента со степенями свободы v где v оценивают по формуле (8.34) Значения квантилей распределения Стьюдента приведены в таблице Б.2 приложения Б. 8.2.3.8 При выполнении неравенства (8.33) для оценивания аттестованного значения СО принимают результаты, полученные по всем МВИ. При невыполнении неравенства (8.33) проводят проверку совместимости средних значений по 8.2.3.7 после исключения из ряда (8.32) или . 8.2.3.9 В случае принятия положительного решения в соответствии с 8.2.3.7 аттестованное значение оценивают следующим образом: - вычисляют для среднего значения результата j-й МВИ вес по формуле (8.35) и сумму весов по формуле (8.36) Аттестованное значение СО А вычисляют как средневзвешенное значение по формуле (8.37) Стандартную неопределенность от способа определения аттестованного значения СО uchar вычисляют по формуле (8.38) Определяют число степеней свободы vchar (8.39) 8.3 Оценивание неопределенности от способа определения аттестованного значения стандартного образца по методике выполнения измерений в одной лаборатории8.3.1 Способы определения аттестованного значения стандартного образца в одной лаборатории В соответствии с перечислением, приведенным в 8.1, применяют три способа определения аттестованного значения СО в одной лаборатории: - с использованием эталона; - с использованием МВИ, аттестованной по ГОСТ 8.010; - в соответствии с процедурой приготовления. 8.3.2 Оценивание неопределенности от способа определения аттестованного значения стандартного образца по результатам измерений, полученных с помощью эталона 8.3.2.1 Эталоны, используемые для определения аттестованного значения СО, должны иметь действующие сертификаты о калибровке, содержащие значения расширенной неопределенности U(q) с указанием коэффициента охвата k для исключенных систематических эффектов и стандартного отклонения результатов измерений в условиях повторяемости sr. 8.3.2.2 При выборе конкретных эталонов для определения аттестованных значений СО в общем случае следует руководствоваться отношениями U(q) £ Uдоп; sr £ 1,2Uдоп, (8.40) где Uдоп - допускаемое значение расширенной неопределенности аттестованного значения СО в соответствии с техническим заданием на разработку СО. 8.3.2.3 Измерения при определении аттестованного значения СО проводят согласно правилам выполнения измерений на эталоне или МВИ, рекомендуемой разработчиком СО и учитывающей правила применения соответствующего эталона. 8.3.2.4 Для СО монолитных или дисперсных материалов предварительно оценивают стандартную неопределенность от неоднородности uh в соответствии с требованиями раздела 7. При поэкземплярной аттестации СО стандартную неопределенность от неоднородности не оценивают и в дальнейших расчетах принимают uh равным 0. 8.3.2.5 Общее число измерений J, необходимых для оценки аттестованного значения СО методом многократных наблюдений, определяют следующим образом. Вычисляют: - допускаемое значение стандартного отклонения среднего значения результатов в условиях повторяемости sдоп по формуле (8.41) - значения величин x и h (8.42) - определяют число измерений J по таблице 8.1 в зависимости от значений величин x и h. Таблица 8.1 - Число наблюдений J при определении аттестованного значения СО методом многократных наблюдений 8.3.2.6 Аттестованное значение СО и стандартную неопределенность uchar CO оценивают следующим образом: - получают в условиях повторяемости результаты п измерений x1, ..., хп, где п ³ J; - вычисляют среднее значение результатов измерений и их стандартное отклонение s по формулам (8.43) - сравнивают полученное значение стандартного отклонения со стандартным отклонением повторяемости по критерию хи-квадрат по неравенству где v = n - 1; - квантиль c2 распределения с v степенями свободы, значения которого приведены в таблице Б.1 приложения Б; - если неравенство (8.44) выполняется, то аттестованное значение СО вычисляют как среднее значение результатов измерений (8.45) и стандартную неопределенность от способа определения аттестованного значения СО вычисляют по формуле (8.46) Для стандартной неопределенности uchar число степеней свободы определяют по формуле vchar = n - 1. (8.47) 8.3.2.7 При невыполнении неравенства (8.44) полученные результаты измерений бракуют и повторяют процедуру, указанную в 8.3.2.6. 8.3.3 Оценивание неопределенности от способа определения аттестованного значения стандартного образца по результатам измерений с использованием аттестованной методики выполнения измерений 8.3.3.1 Для установления аттестованного значения СО может быть использована МВИ с известными неопределенностями. Предварительно в процессе аттестации МВИ определяют стандартную неопределенность повторяемости sr, стандартную неопределенность систематического смещения результатов измерений sJ, полученных по МВИ для данной лаборатории, и стандартную неопределенность промежуточной прецизионности sI0 (в скобках указывают идентификацию промежуточной ситуации в соответствии с [5]). Составляющими систематического смещения могут быть систематическое смещение результатов МВИ и систематическое смещение результатов измерений, полученных в лаборатории. Стандартная неопределенность промежуточной прецизионности sI(), оцениваемой в соответствии с [5], обусловлена различием между результатами измерений, полученными: - разными операторами; - с использованием разного оборудования; - в разных условиях окружающей среды; - в разное время проведения измерений. В зависимости от имеющихся возможностей лаборатории могут проводить измерения по МВИ с изменением всех вышеперечисленных факторов. 8.3.3.2 Измерения для определения аттестованного значения СО выполняют несколько операторов при различных уровнях промежуточной прецизионности. Каждый оператор на одном уровне промежуточной прецизионности выполняет п измерений в условиях повторяемости. Оптимальное число результатов п измерений находится в интервале от двух до четырех. Результаты измерений, полученные в условиях промежуточной прецизионности, обрабатывают следующим образом. Вычисляют: - для каждого j-го уровня среднее значение в условиях повторяемости по формуле (8.48) - оценку стандартного отклонения повторяемости srj для j-го уровня промежуточной прецизионности по формуле (8.49) - оценку стандартного отклонения повторяемости sr или всех уровней промежуточной прецизионности по формуле (8.50) - среднее значение результатов измерений по всем уровням промежуточной прецизионности по формуле где z - число уровней промежуточной прецизионности; - стандартное отклонение промежуточной прецизионности sI(0), обусловленное участием различных операторов, по формуле (8.52) 8.3.3.3 Сравнивают оценку стандартного отклонения повторяемости sr МВИ со стандартным отклонением повторяемости sr по критерию хи-квадрат, проверяя выполнение неравенства где v - число степеней свободы, v = р(п - 1); - квантиль c2 распределения с v степенями свободы, значения которых приведены в таблице Б.1 приложения Б. При выполнении неравенства (8.53) повторяемость результатов для всех уровней промежуточной прецизионности считают приемлемой. Для дальнейшей обработки используют результаты измерений, полученные на всех уровнях. В противном случае выявляют уровни с набольшими значениями стандартной неопределенности повторяемости и исключают результаты измерений, полученные на этих уровнях. При необходимости проводят дополнительные измерения. 8.3.3.4 Сравнивают оценку стандартного отклонения промежуточной прецизионности sI(0) и стандартное отклонение промежуточной прецизионности sI(0) по критерию хи-квадрат, описываемому неравенством где - квантиль c2 распределения c v = р - 1 степенями свободы, значения которых приведены в таблице Б.1 приложения Б. 8.3.3.5 При выполнении неравенства (8.54) стандартное отклонение sI(0) является приемлемым и результаты измерений на всех уровнях используют для оценивания аттестованного значения СО и стандартной неопределенности от способа определения аттестованного значения СО. При невыполнении неравенства (8.54) стандартное отклонение sI(0) является неприемлемым. В этом случае возможны следующие варианты: а) использовать критерий для выявления статистических выбросов (например, критерий Граббса в соответствии с пунктом 7.3.4 национального стандарта [6]) и повторить процедуру проверки приемлемости результатов в соответствии 8.2.3.3 - 8.2.4; б) забраковать все результаты при их малом числе. 8.3.3.6 При выполнении неравенства (8.54) аттестованное значение СО А вычисляют как среднее арифметическое по всем результатам измерений по формуле (8.51) (8.55) Стандартную неопределенность от способа определения аттестованного значения СО определяют по формуле (8.56) Для стандартной неопределенности uchar число степеней свободы вычисляют по формуле (8.57) где vJ - число степеней свободы стандартного отклонения sj. 8.3.4 Оценивание неопределенности от способа определения аттестованного значения стандартного образца по результатам косвенных измерений 8.3.4.1 Результаты косвенных измерений используют при определении аттестованного значения стандартного образца по процедуре приготовления. В соответствии с этим способом используют при оценивании аттестованного значения СО оценки влияющих величин уi, а при оценивании стандартной неопределенности от способа определения аттестованного значения стандартного образца - оценки стандартных неопределенностей влияющих величин u(уi) 8.3.4.2 Аттестуемую характеристику СО выражают в виде зависимости от влияющих величин уi, где i = 1, 2, ..., n, а аттестованное значение стандартного образца А оценивают по оценкам значений влияющих величин уi, используя зависимость (8.58) (8.59) 8.3.4.3 Стандартную неопределенность способа uchar оценивают по стандартным неопределенностям влияющих величин и их коэффициентам корреляции по формуле (8.60) где r(уi, уj) - коэффициент корреляции влияющих величин уi и уj. В случае статистической независимости влияющих величин значение неопределенности от способа вычисляют по формуле (8.61) Число степеней свободы стандартной неопределенности uchar рассчитывают по формуле 8.3.4.4 Степени свободы влияющих величин вычисляют следующими способами: - если неопределенность влияющей величины уi вычислена с помощью статистического анализа серии из ni измерений (оценивание типа А), то число степеней свободы равно vi = ni - k, (8.63) где k - число параметров, оцененных по данной серии измерений. В частности, если по серии результатов оценено только среднее значение (математическое ожидание), то k = 1 и vi = ni - k; - если неопределенность влияющей величины уi оценена иным способом без проведения статистического анализа серии измерений (оценивание типа В), то число степеней свободы устанавливают по имеющейся информации о неопределенности оценки стандартной неопределенности. В общем случае согласно [1] vi можно вычислить по формуле (8.64) - если известно, что влияющая величина с вероятностью, равной единице, принимает значения в известном интервале и оценивание типа В стандартной неопределенности получено с учетом известного распределения результатов измерений (например, равномерного, треугольного), то число степеней свободы равно бесконечности (¥). В этом случае вычисление числа степеней свободы стандартной неопределенности uchar(А) по формуле (8.62) имеет вид
9 Оценивание суммарной стандартной и расширенной неопределенности аттестованного значения стандартного образца9.1 Суммарную стандартную неопределенность аттестованного значения СО ис(А) вычисляют по формуле (9.1) где составляющие uchar, uh и ustab вычислены по формулам, приведенным в разделах 6 - 8. Число эффективных степеней свободы суммарной стандартной неопределенности ис(А) вычисляют по формуле 9.2 Расширенную неопределенность аттестованного значения СО UP(A) для доверительной вероятности Р вычисляют по формуле UP(A) = tp(vc)uc(A), (9.3) где tp(vc) - квантиль распределения Стьюдента с числом степеней свободы vc при доверительной вероятности Р. Значения квантилей распределения Стьюдента для доверительной вероятности Р = 0,95 приведены в таблице Б.2 приложения Б. Приложение А
|
|
v |
|
v |
|
|
1 |
3,841 |
11 |
19,675 |
22 |
33,924 |
2 |
5,991 |
12 |
21,026 |
24 |
36,415 |
3 |
7,815 |
13 |
22,362 |
26 |
38,885 |
4 |
9,488 |
14 |
23,685 |
28 |
41,337 |
5 |
11,070 |
15 |
24,996 |
30 |
43,773 |
6 |
12,592 |
16 |
26,296 |
32 |
46,194 |
7 |
14,067 |
17 |
27,587 |
34 |
48,602 |
8 |
15,507 |
18 |
28,869 |
36 |
50,998 |
9 |
16,919 |
19 |
30,144 |
38 |
53,384 |
10 |
18,307 |
20 |
31,410 |
40 |
55,758 |
Для v ³ 41 значение вычисляют по формуле
= 1,36v + 10,8.
Таблица Б.2 - Квантили распределения Стьюдента t0,95(v)
t0,95(v) |
v |
t0,95(v) |
v |
t0,95(v) |
|
1 |
12,706 |
11 |
2,201 |
22 |
2,074 |
2 |
4,303 |
12 |
2,179 |
24 |
2,064 |
3 |
3,182 |
13 |
2,160 |
26 |
2,056 |
4 |
2,776 |
14 |
2,145 |
28 |
2,048 |
5 |
2,571 |
15 |
2,131 |
30 |
2,042 |
6 |
2,447 |
16 |
2,120 |
32 |
2,037 |
7 |
2,365 |
17 |
2,110 |
34 |
2,032 |
8 |
2,306 |
18 |
2,101 |
36 |
2,028 |
9 |
2,262 |
19 |
2,093 |
38 |
2,024 |
10 |
2,228 |
20 |
2,086 |
40 |
2,021 |
Для v ³ 41 значение t0,95(v) вычисляют по формуле
Guid to expression of uncertainty in measurement. ISO, Geneva, 1993 (Руководство по выражению неопределенности измерения / Пер. с англ. под ред. В.А. Слаева. С-Пб.: ВНИИМ, 1999) |
||
Руководство ИСО/РЕМКО 35:2006 (ISO/REMCO Guide 35:2006) |
Стандартные образцы. Общие и статистические принципы аттестации (Reference materials. General and statistical principles for certification) |
|
Guid: Quantifying uncertainty in analytical measurement, 2nd edition, 2000, Eurachem/CITAC (Количественное описание неопределенности в аналитических измерениях. Руководство ЕВРАХИМ/СИТАК. 2-е изд. / Пер. с англ. под ред. Л.А. Конопелько С-Пб.: ВНИИМ, 2002 |
||
Точность (правильность и прецизионность) методов и результатов измерений. Часть 1. Основные положения и определения |
||
Точность (правильность и прецизионность) методов и результатов измерений. Часть 3. Промежуточные показатели прецизионности стандартного метода измерений |
||
Точность (правильность и прецизионность) методов и результатов измерений. Часть 2. Основной метод определения повторяемости и воспроизводимости стандартного метода измерений |
||
Точность (правильность и прецизионность) методов и результатов измерений. Часть 6. Использование значений точности на практике |
Ключевые слова: стандартный образец, аттестованное значение, неопределенность, стандартная неопределенность аттестованного значения стандартного образца, однородность стандартного образца, стабильность стандартного образца