"ГОСТ Р ИСО 13752-2005. Национальный стандарт Российской Федерации. Качество воздуха. Оценка неопределенности метода измерений в условиях применения с использованием референтного метода"
(утв. и введен в действие Приказом Ростехрегулирования от 22.12.2005 N 321-ст)

Цели и принципы стандартизации в Российской Федерации установлены Федеральным законом от 27 декабря 2002 г. N 184-ФЗ "О техническом регулировании", а правила применения национальных стандартов Российской Федерации - ГОСТ Р 1.0-2004 "Стандартизация в Российской Федерации. Основные положения".

1. Подготовлен Открытым акционерным обществом "Научно-исследовательский центр контроля и диагностики технических систем" (ОАО "НИЦ КД") на основе собственного аутентичного перевода стандарта, указанного в пункте 4.

2. Внесен Техническим комитетом по стандартизации ТК 457 "Качество воздуха".

3. Утвержден и введен в действие Приказом Федерального агентства по техническому регулированию и метрологии от 22 декабря 2005 г. N 321-ст.

4. Настоящий стандарт идентичен международному стандарту ИСО 13752:1998 "Качество воздуха. Оценка неопределенности метода измерений в условиях применения с использованием второго метода в качестве референтного" (ISO 13752:1998 "Air quality - Assessment of uncertainty of a measurement method under field conditions using a second method as reference").

Наименование настоящего стандарта изменено относительно наименования указанного международного стандарта для приведения в соответствие с ГОСТ Р 1.5 - 2004 (подраздел 3.5).

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

5. Введен впервые.

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

Характеристики методов измерений качества воздуха установлены ИСО 6879. Соответствующие процедуры их определения приведены в ИСО 9169, за исключением точности, в настоящем стандарте рассматриваемой как неопределенность измерений в соответствии с требованиями "Руководства по выражению неопределенности измерений".

Неопределенность измерений в условиях применения рассматривается также в стандартах ИСО 7935:1992 "Выбросы стационарных источников - Определение массовой концентрации диоксида серы - Параметры автоматических методов измерений" и ИСО 10849:1996 "Выбросы стационарных источников - Определение массовой концентрации оксидов азота - Рабочие характеристики автоматических измерительных систем". Однако процедуры, приведенные в этих стандартах, ограничиваются определением либо независимой от концентрации систематической погрешности при допущении независимой от концентрации дисперсии, либо пропорциональной концентрации систематической погрешности при допущении относительно пропорциональной концентрации дисперсии.

Настоящий стандарт устанавливает процедуру оценки неопределенности измерений, выполняемых "полевым" <1> методом, подлежащим проверке (далее - проверяемый метод), по другому (референтному) методу. Используемый референтный метод не обязательно должен быть стандартизован.

--------------------------------

<1> Метод измерений в реальных условиях применения (далее - условия применения).

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

Настоящая процедура разработана специально для оценки пригодности (валидации) проверяемого метода.

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

ИСО 6879:1995 Качество воздуха - Характеристики и соответствующие им понятия, относящиеся к методам измерений качества воздуха

ИСО 9169:1994 Качество воздуха - Определение характеристик методов измерений.

В настоящем стандарте применены следующие обозначения и сокращения <1>:

--------------------------------

<1> Определения применяемых терминов можно найти в [1] и [2].

--------------------------------

<1> Данные обозначения в тексте стандарта не употребляются, приведены для сохранения идентичности ИСО 13752:1998.

На основе параллельных измерений в условиях применения получают N пар результатов измерений . Результаты измерений, полученные референтным методом (x-методом), рассматривают как истинные. Разность результатов измерений в паре приписывают отклонению результата измерений, полученному проверяемым методом (y-методом).

При выборе модели зависимости между переменными величинами X и Y принимают допущение об их линейной зависимости:

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

- стандартное отклонение проверяемого метода не зависит от x (т.е. стандартное отклонение является постоянным) и задается формулой

- стандартное отклонение проверяемого метода пропорционально x (т.е. коэффициент вариации постоянен) и задается формулой

Примечание 1. Первое допущение рассматривает флуктуации фона или значения коэффициента (свободного члена - отсекаемого отрезка <1>) без учета флуктуаций коэффициента (коэффициента наклона); второе - флуктуации коэффициента наклона без учета флуктуаций фона или свободного члена.

--------------------------------

<1> См. ИСО 6879.

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

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

Примечание 3. Коэффициенты возведены в квадрат, т.к. коэффициент в большей мере, чем его квадрат, отражает физический смысл.

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

Значения коэффициентов принятой модели дисперсии общего вида [ и в формуле (1) и , и в формуле (4)] не могут быть рассчитаны. Их подбирают на основе критерия максимального правдоподобия итеративным методом. Блок-схема вычисления коэффициентов функции линейной регрессии методом максимального правдоподобия приведена на рисунке 1. После выбора начальных значений коэффициентов, при использовании допущения о нормальности закона распределения, вычисляют вероятность в каждой точке, отвечающей результатам измерений и принадлежащей линии регрессии, по формуле

Функцию правдоподобия L, являющуюся произведением отдельных вероятностей результатов измерений, полученных y-методом, вычисляют по формуле

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

Неопределенность измерений, соответствующую любому полученному значению ХКВ, оценивают по функции регрессии и функции дисперсии.

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

Процедуру, установленную настоящим стандартом, применяют при выполнении следующих условий:

- существует линейная зависимость между сравниваемыми переменными; если зависимость нелинейна, но известно ее математическое выражение, процедура может быть адаптирована;

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

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

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

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

5.2. Проверяемый метод (y-метод)

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

5.3. Референтный метод (x-метод)

Проверяют допущение о незначительности неопределенности x-метода по сравнению с y-методом с учетом условий проведения измерений (в том числе условий окружающей среды), ожидаемых на месте проведения проверки, например наличия мешающих веществ, температуры и т.д. Проверку допущения о неопределенности x-метода проводят на основе анализа принципа измерений, литературных данных или результатов испытаний в лабораторных условиях или в условиях применения. Ориентируясь на детальное описание x-метода, выполняют соответствующие измерения.

5.4. Условия проведения измерений

Условия проведения измерений должны соответствовать условиям, в которых планируют использовать проверяемый метод (продолжительность измерений, диапазон изменений ХКВ, диапазон изменений физических и химических влияющих величин и условия работы). Фиксируют условия, в том числе окружающей среды, на месте проведения измерений.

Средства измерений, используемые в обоих методах, должны быть установлены так, чтобы:

- различия в составе параллельных проб были незначительны;

- средства измерений одного метода не влияли на средства измерений другого метода.

5.5. Обработка данных

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

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

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

- если значения разностей не зависят от , см. 8.2;

- если значения разностей являются пропорциональными , см. 8.3;

- если значения разностей не являются ни независимыми, ни пропорциональными , см. 8.4.

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

--------------------------------

<1> Диапазон измерений может быть разбит на поддиапазоны, в каждом из которых должно быть выполнено конкретное условие из перечисленных.

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

Подразделы 8.2 и 8.3 содержат процедуры проверки, позволяющие выяснить, подтверждается или нет допущение относительно принятой модели дисперсии.

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

Принятая линейная зависимость между переменными X и Y описывается уравнением линейной регрессии

Если дисперсия результатов измерений относительно линии регрессии не зависит от значения ХКВ (см. 8.2), ее определяют по формуле

Если дисперсия относительно линии регрессии пропорциональна значению ХКВ (см 8.3), ее определяют по формуле

В общем случае, если дисперсия относительно линии регрессии является монотонной функцией значения ХКВ (см. 8.4), ее определяют по формуле

Выбирают одну из указанных выше моделей дисперсии.

8.2. Стандартное отклонение постоянно

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

Проверяют, подтверждается или нет допущение относительно постоянства дисперсии:

- выбирают пар результатов измерений вблизи верхнего предела диапазона измерений и пар вблизи нижнего предела при условии ;

- не используют среднюю часть диапазона измерений;

- рассчитывают статистику F по формуле

- если F не превышает табулированное значение для F-распределения в случае одностороннего критерия для уровня значимости , принятое за критическое значение, дисперсию считают постоянной;

- если F превышает табулированное значение F, см. 8.3 или 8.4.

Дисперсию вычисляют по формуле

Стандартные отклонения и коэффициентов и вычисляют по формулам:

8.3. Коэффициент вариации постоянен

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

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

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

- выбирают пар результатов измерений в области верхней границы диапазона измерений и пар в области нижней границы с условием ;

- не используют среднюю часть диапазона;

- статистику F вычисляют по формуле

- если F не превышает табулированное значение для F-распределения в случае одностороннего критерия для уровня значимости , принятое за критическое значение, дисперсию считают постоянной;

- если F превышает табулированное значение F, см. 8.2 или 8.4.

Функцию дисперсии вычисляют по формуле

Проводят обратное преобразование полученных значений по формулам:

Стандартные отклонения и коэффициентов и вычисляют по формулам:

8.4. Модель дисперсии общего вида

Выбирают начальные значения коэффициентов модели и [формула (7)] и , , [формула (10)], где .

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

Логарифм функции правдоподобия l вычисляют по формуле

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

Используют значения коэффициентов , , , и , соответствующие максимальному L.

Стандартные отклонения и коэффициентов и вычисляют по формулам:

где

Образец крупноформатной таблицы для расчета коэффициентов функции регрессии и функции дисперсии приведен в Приложении A, а соответствующий пример крупноформатной таблицы - в Приложении B.

Примечания

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

2. Обычно функция дисперсии имеет вид . Эта функция отражает стремление дисперсии к константе в области нижней границы диапазона измерений и пропорциональность дисперсии значению ХКВ в области верхней границы диапазона измерений. Три члена в выражение для функции дисперсии вводят только при большом числе пар измерений.

3. Если зависимость y = f(x) является нелинейной, но ее математический вид известен, формула (28) может быть заменена соответствующей математической функцией, а коэффициенты регрессии определены аналогично методом максимального правдоподобия.

Значения коэффициентов и существенно отличаются от идеальных значений 0 и 1 соответственно, если:

Если формулы (36) и (37) показывают значимость поправок, систематическая погрешность при X = x может быть рассчитана в пределах диапазона измерений по формуле

В соответствии с процедурами, указанными в [3], систематическая погрешность, возникающая в результате выявленного эффекта, может быть компенсирована путем введения поправки. Поправки могут быть введены, только если значения влияющих переменных, ответственных за эту систематическую погрешность, представительные. Однако неопределенность поправок, приравниваемая к неопределенности систематической погрешности, остается:

где :

a) при модели дисперсии 8.2:

b) при модели дисперсии 8.3:

c) при модели дисперсии 8.4: см. формулы (34) и (35).

При X = x вычисляют стандартное отклонение результата измерений в условиях применения на основе функции дисперсии по формуле

Предполагая, что может быть введена поправка на систематическую погрешность , вычисляют стандартное отклонение исправленного значения по формуле

Неопределенность единичного результата измерений при X = x с коэффициентом охвата k = 2 определяют по формуле

Коэффициент охвата k = 2 соответствует коэффициенту Стьюдента для доверительного интервала 95% и нормального распределения.

В [3] рекомендуется всегда вводить поправку на систематическую погрешность. Обычно в инструкциях по оценке неопределенности измерений не приводят способы определения и введения поправок на систематические погрешности. Если введение поправки на систематические погрешности не является частью y-метода, эти погрешности включают в неопределенность измерений. Применяя согласно [3] принцип выражения погрешностей в показателях дисперсии, неопределенность измерений вычисляют по формуле

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

В настоящем приложении приведен пример расчета коэффициентов функции регрессии и функции дисперсии с помощью программы табличных вычислений <1>. Вводят формулы из таблицы A.1 в соответствующие ячейки пустой формы таблицы A.2.

--------------------------------

<1> - торговое наименование программы, предоставляемой Microsoft. Возможно использование других программ, приводящее к получению аналогичных результатов.

--------------------------------

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

--------------------------------

<1> Руководство по выражению неопределенности измерения: Аутентичный перевод с англ./Под ред. проф. Слаева В.А. - СПб.: Изд-во ВНИИМ им. Д.И. Менделеева, 1999.