Примеры расчета неопределенности измерений - Школа метрологии
Главная » Полезные статьи » Примеры расчета неопределенности измерений

Примеры расчета неопределенности измерений

Оценивание неопределенности измерений (количественных величин) – одна из важных задач, стоящих перед каждой лабораторией. Требование к оцениванию неопределенности измерений заложено в межгосударственном стандарте ГОСТ ISO/IEC 17025-2019, а также политике ILAC-G17:2002.

Чтобы получить предварительную информацию по теме Неопределенность измерений, посмотрите, пожалуйста, наше обучающее видео:

Международное метрологическое сообщество давно уже разработало и приняло основные принципы концепции неопределенности, закрепив их в серии международных документов JCGM (Joint Committee for Guides in Metrology), а также документах ISO/IEC Guide 98. Разработано много дополнительных руководств по различных подходам оценивания неопределенности измерений в конкретных областях испытаний/измерений (EA, EURACHEM, Nordtest, EUROLAB и т.д.).

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

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

Пример 1. Оценивание неопределенности измерений массовой доли влаги

Исходные данные:

  1. Объект измерений – углекислый барий
  2. Измеряемая величина – массовая доля влаги
  3. Единицы измерений – процент (%)
  4. Методика выполнения измерения – ГОСТ 2149-75 «Барий углекислый технический. Технические условия»
  5. Метод измерений – метод высушивания пробы до постоянной массы

Этап 1. Составление функции измерений

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

Массовую долю влаги бария углекислого Х в процентах вычисляют на основании ГОСТ 2149 (п.3.5) в соответствии с функцией измерений:

1
2

где Хi – массовая доля влаги i-ой пробы углекислого бария, %;

m1 – масса стаканчика для взвешивания с навеской до высушивания, г;

m2 – масса стаканчика для взвешивания с навеской после высушивания, г;

m – масса навески углекислого бария, г;

m – масса стаканчика для взвешивания, г;

i – номер параллельной пробы, i = 1, 2;

F – поправочный множитель, учитывающий допускаемое расхождение между параллельными определениями.

Этап 2. Анализ входных величин

1) Масса стаканчика для взвешивания с навеской до высушивания, m1

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

  • для первой пробы m1(1) = 41,0055 г,
  • для второй пробы m1(2) = 41,3842 г.

Неопределенность, связанную с величиной m1, оцениваем, используя данные производителя на весы. В паспорте на весы лабораторные AV264C для диапазона измерений до 50 г указаны пределы погрешности взвешивания ± 0,001 г. Поскольку значение дано без доверительной вероятности, принимаем прямоугольное распределение значений погрешности взвешивания в этих границах. Стандартная неопределенность массы стаканчика для взвешивания с навеской до высушивания m1 оценивается по типу В и составляет:

2) Масса стаканчика для взвешивания, m

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

  • для первой пробы mст(1) = 21,0034 г,
  • для второй пробы mст(2) = 21,3822 г.

Неопределенность массы стаканчика для взвешивания после высушивания обусловлена двумя факторами:

погрешностью взвешивания стаканчика на весах;
возможными отклонениями массы стаканчика после высушивания вследствие нечеткого определения в методике выполнения измерений момента, в который масса стаканчика после высушивания будет являться постоянной величиной.
Стандартная неопределенность u1(m), связанная с погрешностью взвешивания, оценивается на основании данных производителя на весы лабораторные AV264C, определяется аналогично неопределенности величины m1 и составляет u1(m) = 0,00058 г.

Стандартную неопределенность u2(m), обусловленную отклонениями массы стаканчика для взвешивания после высушивания, можно определить на основании информации о том, что разность между двумя последующими взвешиваниями стаканчика для взвешивания после сушки не должна превышать 0,002 г (принято лабораторией, поскольку ГОСТ 2149-75 четко не устанавливает момент, в который масса стаканчика после высушивания будет являться постоянной величиной). Это значение можно рассматривать как 95-% доверительный интервал для разности двух оценок величины, распределенной по нормальному закону распределения (предел повторяемости по СТБ ИСО 5725-6). Стандартная неопределенность u2(m) будет равна стандартному отклонению, рассчитанному на основании указанного интервала по типу В по формуле:

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

3) Масса стаканчика для взвешивания с навеской после высушивания, m2

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

  • для первой пробы m2(1) = 40,9850 г,
  • для второй пробы m2(2) = 41,3638 г.

Неопределенность массы стаканчика для взвешивания с навеской после высушивания обусловлена двумя факторами:

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

Стандартная неопределенность u1(m2), связанная с погрешностью взвешивания, оценивается на основании данных производителя на весы лабораторные AV264C, определяется аналогично неопределенности величины m2 и составляет u1(m2) = 0,00058 г.

Стандартную неопределенность u2(m2), обусловленную отклонениями массы стаканчика для взвешивания с навеской после высушивания, можно определить на основании информации о том, что разность между двумя последующими взвешиваниями стаканчика для взвешивания с навеской пробы бария углекислого после сушки не должна превышать 0,002 г (принято лабораторией, поскольку ГОСТ 2149-75 четко не устанавливает момент, в который масса пробы после высушивания будет являться постоянной величиной). Это значение можно рассматривать как 95-% доверительный интервал для разности двух оценок величины, распределенной по нормальному закону распределения (предел повторяемости по СТБ ИСО 5725-6). Стандартная неопределенность u2(m2) будет равна стандартному отклонению, рассчитанному на основании указанного интервала по типу В по формуле:

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

4) Поправочный множитель, учитывающий допускаемое расхождение между параллельными определениями, F

Значение оценки величины принимается равным единице: F = 1.

Стандартная неопределенность поправочного множителя рассчитывается на основании информации о допускаемом расхождении между параллельными определениями влажности, приведенными в ГОСТ 2149-75 (п.3.5.2). Приведенное допускаемое относительное расхождение составляет r = 20 % и рассматривается как 95-% доверительный интервал для разности двух оценок величины, распределенной по нормальному закону распределения. Стандартная неопределенность поправочного множителя будет равна стандартному отклонению, рассчитанному на основании указанного интервала, с учетом того, что за результат измерения принимают среднее арифметические определений двух параллельных проб, по типу В по формуле:

Этап 3. Анализ корреляций

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

Этап 4. Измеренное значение величины

Массовая доля влаги i-ой пробы углекислого бария рассчитывается по формуле (2):

1-2

Расхождение между результатами измерений двух параллельных проб не превышает допускаемого значения, установленного в ГОСТ 2149-75 (п.3.5.2):

1-3

Массовая доля влаги бария углекислого Х в процентах вычисляется в соответствии с функцией измерений (1):

1-4

Измеренное значение (оценку измеряемой величины) округляют до четырех знаков после запятой (принято лабораторией, поскольку ГОСТ 2149-75 четко не устанавливает требования к округлению результата измерения).

Этап 5. Суммарная стандартная неопределенность

Стандартную неопределенность измеряемой величины Х получаем по закону распространения неопределенностей путем суммирования квадратов произведений стандартных неопределенностей всех влияющих величин, входящих в функции измерений (1) и (2), на соответствующие коэффициенты чувствительности:

9

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

10 11 12 13

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

Суммарная стандартная неопределенность составит

1-1

Этап 6. Бюджет неопределенности

В таблице представлен бюджет неопределенности для измеряемой величины.

таблица

Вклад в неопределенность от j-ой входной величины рассчитывается как произведение стандартной неопределенности этой величины на соответствующий коэффициент чувствительности.

Этап 7. Расширенная неопределенность

Расширенную неопределенность U получаем умножением суммарной стандартной неопределенности на коэффициент охвата k, который выбирается равным 2 при уровне доверия приблизительно 95 % в предположении нормального распределения вероятностей измеряемой величины:

14

Этап 8. Представление результата измерения

Результат измерения массовой доли влаги углекислого бария представляют в виде:

«Массовая доля влаги бария углекислого составила (0,1022 ± 0,0150) %, где число, следующее за знаком ±, является численным значением расширенной неопределенности, которая получена умножением суммарной стандартной неопределенности на коэффициент охвата k = 2, основанный на предполагаемом нормальном распределении, и определяет интервал, соответствующий вероятности охвата приблизительно 95 %».

Ближайшее обучение

Успейте записаться на ближайшее мероприятие

Популярные услуги