ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ УНИТАРНОЕ ПРЕДПРИЯТИЕ
ВСЕРОССИЙСКИЙ НАУЧНО-ИССЛЕДОВАТЕЛЬСКИЙ
ИНСТИТУТ МЕТРОЛОГИЧЕСКОЙ СЛУЖБЫ

(ФГУП ВНИИМС)

ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО
ТЕХНИЧЕСКОМУ РЕГУЛИРОВАНИЮ И МЕТРОЛОГИИ

УТВЕРЖДАЮ
Зам. директора ФГУП «ВНИИМС»
В.Н. Яншин
24 мая 2006 г.

РЕКОМЕНДАЦИЯ

Государственная система обеспечения единства измерений.

Объем и масса нефти и нефтепродуктов.
Методика контроля величины погрешности
счетчиков ультразвуковых Altosonic V (VM)
по результатам сличений
двух равноточных приборов

 

МИ 2987-2006

 

 

Москва

2006

 

ПРЕДИСЛОВИЕ

1 РАЗРАБОТАНА ФГУП ВНИИМС

ИСПОЛНИТЕЛИ: Вересков А.И. (рук. темы);

Беляев Б.М.;

Дудыкин А.А.

2 УТВЕРЖДЕНА ФГУП ВНИИМС 24 мая 2006 г.

3 ЗАРЕГИСТРИРОВАНА ФГУП ВНИИМС 25 мая 2006 г.

4 ВВЕДЕНА ВПЕРВЫЕ

 

СОДЕРЖАНИЕ

 

РЕКОМЕНДАЦИЯ

1 ОБЛАСТЬ ПРИМЕНЕНИЯ

1.1 Настоящая рекомендация распространяется на объем и массу нефти и нефтепродуктов (далее - продукт) и устанавливает методику контроля метрологических характеристик канала измерений объема (массы) счетчиков ультразвуковых Altosonic V (VM) (далее - прибор), применяемых в составе систем измерения количества и показателей качества нефти (СИКН) или нефтепродуктов (СИКНП) по результатам сличений двух однотипных приборов по объему (массе) продукта.

1.2 Методика предназначена для контроля метрологических характеристик (далее - MX) приборов в межповерочный период.

1.3 Контроль метрологических характеристик по настоящей методике не заменяет поверки в соответствии с ПР 50.2.006.

1.4 Рекомендация разработана с учетом требований ГОСТ 8.207 и МИ 1317.

2 НОРМАТИВНЫЕ ССЫЛКИ

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

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

ПР 50.2.006-94 ГСИ. Порядок проведения поверки средств измерений

МИ 1317-2004 ГСИ. Результаты и характеристики погрешности измерений.

Формы представления. Способы использования при испытаниях образцов продукции и контроле их параметров.

3 ОПИСАНИЕ МЕТОДА КОНТРОЛЯ

Контроль MX проводится в два этапа:

Этап I. Формирование исходных данных для контроля.

Этап II. Обработка результатов сличений.

3.1 Этап I. Формирование исходных данных для контроля.

Формируют исходные данные для контроля. При этом применяют указания изложенные в п.3.1.1 - 3.1.3.

3.1.1 Результаты измерений объема (массы) продукта (далее - количество продукта) при последовательной установке в трубопроводе (измерительной линии) двух приборов (далее - исходные данные):

Q_1i, Q_2i, t_i, i = 1, ..., N,

где Q_1i, Q_2i - результаты измерений количества продукта при i-м сличении соответственно первым и вторым прибором за время (продолжительность) сличения t_i;

N - число сличений, N ³ 1.

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

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

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

При проведении контроля MX при постоянном расходе будут получены результаты контроля MX приборов, соответствующие этому значению расхода. Расход считают постоянным, если его значение изменяется при формировании исходных данных не более чем на 5 %.

При проведении контроля MX в диапазоне расходов рекомендуется ограничивать диапазон расходов таким образом, чтобы отношение максимального значения расхода к минимальному было не более 1,5.

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

3.1.2 Характеристики погрешности годного прибора

В качестве характеристик погрешности годного прибора могут быть применены следующие характеристики.

3.1.2.1 Пределы допускаемой относительной погрешности при измерении объема или массы: ± δ, %.

3.1.2.2 Пределы неисключенной систематической относительной погрешности прибора δ_CL, % и среднего квадратического отклонения (СКО) случайной составляющей погрешности прибора σ_L, %.

3.1.3 Число α - вероятность получить отрицательный результат контроля, когда MX приборов удовлетворяют существующим требованиям (далее - вероятность ошибки первого рода). Используется для формирования критерия принятия решения о результате контроля. Рекомендации по выбору значения α даны в Приложении А.

3.2 Этап II. Обработка результатов сличений.

3.2.1 Контроль MX приборов проводят в соответствии с методом, изложенным в Приложении Б.

3.2.2 В результате обработки исходных данных, сформированных по п. 3.1, получают:

3.2.2.1 Ряд результатов сличений d_i, основанных на результатах измерений количества продукта Q_1i, Q_2i (п. 3.1.1):

d_i, %, I = 1, ..., N

3.2.2.2 При N ³ 10 - заключение об обнаружении дрейфа разности систематических составляющих погрешностей приборов (далее - дрейфа) или о том, что дрейф не обнаружен.

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

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

3.2.2.4 Оценку СКО разности случайных составляющих погрешностей сличаемых приборов.

3.2.2.5 Критическое значение С, %, рассчитанное на основании исходных данных п.п. 3.1.2, 3.1.3 с учетом числа сличений N. Критическое значение С характеризует максимальное допускаемое расхождение результатов измерений количества продукта двумя приборами.

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

В противоположном случае результат контроля MX приборов считают отрицательным.

3.2.2.7 При N > 1 в случае отрицательного результата контроля рекомендуется провести контроль MX повторно только по последнему результату сличения из имеющегося ряда результатов сличений. По полученному в соответствии с процедурой п. 3.2.2 результату делают окончательное заключение о результате контроля MX приборов.

3.2.2.8 В случае обнаружения дрейфа погрешности - ожидаемый момент, когда результат контроля MX приборов станет отрицательным.

3.2.3 Методика обработки результатов сличений, приведена в Приложении Б.

4 ПРОВЕДЕНИЕ И ОФОРМЛЕНИЕ РЕЗУЛЬТАТОВ КОНТРОЛЯ

Контроль MX приборов проводят в автоматическом режиме с помощью программы для ЭВМ «Контроль погрешности счетчиков по результатам сличений», разработанной ФГУП ВНИИМС. При работе программы формируется отчет, включающий в себя исходные данные для проведения и результаты контроля MX приборов. Пример отчета приведен в Приложении Г.

ПРИЛОЖЕНИЕ А

РЕКОМЕНДАЦИИ ПО ВЫБОРУ ЗНАЧЕНИЯ ВЕРОЯТНОСТИ ОШИБКИ ПЕРВОГО РОДА

Вероятность ошибочно признать прибор негодным, когда в действительности он годен, принято называть вероятностью ошибки первого рода и использовать для нее обозначение a.

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

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

Используют для формирования критерия принятия решения о положительном или отрицательном результате контроле MX приборов. При выборе значения a руководствуются следующими соображениями. При увеличении a критическое значение С (см. п. 3.2.2.5) уменьшается, т.е. критерий становится более жестким. При этом вероятность ошибки второго рода b (см. Приложение В) также уменьшается.

Рекомендован выбор a из следующего ряда:

0,01; 0,02; 0,03; 0,04; 0,05; 0,07; 0,10; 0,15; 0,20; 0,25; 0,30

или любое число в диапазоне 0,01 £ a £ 0,30.

ПРИЛОЖЕНИЕ Б

АЛГОРИТМ ПРОВЕДЕНИЯ КОНТРОЛЯ MX ПРИБОРОВ

Б.1 Статистический анализ данных для получения модели разности относительных погрешностей приборов

Б.1.1 Результаты сличений d_i рассчитывают, используя результаты измерений количества продукта (п. 3.1.1), по формуле

                                            (Б.1)

Б.1.2 В качестве модели дрейфа разности систематических составляющих погрешностей приборов принимают зависимость вида d = A + B´Q_Σ, где Q_Σ - величина, пропорциональная количеству продукта нарастающим итогом с момента первого сличения.

Примечание - Для расчета Q_Σ используют значения 0,5´(Q_1i + Q_2i), = 1, ..., N.

Б.1.3 Получают оценки параметров А, В по методу наименьших квадратов, при этом:

при N < 10 принимают В = 0,

при N ³ 10 проверяют значимость полученной оценки параметра В. В случае отрицательного результата проверки значимости принимают В = 0.

Результат статистического анализа ряда сличений получают в зависимости от того, обнаружен дрейф или нет: если обнаружен (коэффициент В отличен от нуля), принимают d = A + B´Q_ΣN, а если не обнаружен (B = 0), принимают d = A.

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

Б.2 Расчет максимального допускаемого расхождения показаний приборов

Б.2.1 Для расчета максимального допускаемого расхождения показаний приборов используют формулу для плотности совместного распределения погрешностей приборов:

f(x₁, y₁, x₂, y₂) = f^C(x₁)´f°(y₁)´f^C(x₂)´f°(y₂)                               (Б.2)

При этом функции f^C(x₁), f°(y₁), f^C(x₂), f°(y₂) определены следующими формулами в соответствии с обозначениями п.3.1.2:

f^C(х) = 0 при x £ - δ_CL

f^C(x) =  при -δ_CL < x £ δ_CL                                                                        (Б.3)

f^C(x) = 0 при δ_CL < х

                                       (Б.4)

Б.2.2 Рассчитывают критическое значение С, %.

Функцию f(x₁, y₁, x₂, y₂), определенную в п. Б.2.1, интегрируют по области G_С, в которой переменные удовлетворяют неравенству

|x₁ + y₁ - x₂ - y₂| > С                                                       (Б.5)

при этом С выбирают таким, чтобы интеграл был равен заданному значению a:

                                       (Б.6)

Б.3 Заключение о положительном либо отрицательном результате контроля MX принимают в соответствии с п. 3.2.2.6 путем сопоставления значения результата статистического анализа п. Б.1.3 и критического значения С п. Б.2.2.

Б.4 В случае обнаружения дрейфа получают ожидаемый момент, когда результат контроля MX прибора станет отрицательным. Используют статистический прогноз модели дрейфа по результатам анализа данных п. Б.1.

ПРИЛОЖЕНИЕ В

(СПРАВОЧНОЕ)

РАСЧЕТ ВЕРОЯТНОСТИ ОШИБКИ ВТОРОГО РОДА ПРИ ПРОВЕДЕНИИ КОНТРОЛЯ MX ПРИБОРОВ

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

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

В.1 Для расчета вероятности ошибки второго рода в случае одного негодного прибора задают характеристики его погрешности: конкретные значения систематической составляющей погрешности δ_cw, %, СКО случайной составляющей погрешности σ_w, %.

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

f(x₁, y₁, x₂, y₂) = f₁^C(x₁) ´ f_l⁰(y₁) ´ f^C(x₂) ´ f⁰(y₂)                                  (В.1)

При этом функции f₁^C(x₁), f₁⁰(y₁), f₂^C(x₂), f₂⁰(y₂) определены следующими формулами:

f^C(х₂), f⁰(у₂) - формулами (Б.3), (Б.4);

f₁^C(x₁) = δ(x₁ - δ_cw)                                                             (B.2)

где δ(x₁ - δ_cw) - дельта-функция, равная нулю при х₁, отличном от δ_cw;

f₁⁰(y₁) - формулой (Б.4), в которой σ_N принимают равной

                                                                   (B.3)

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

для первого - значения систематической составляющей δ_1cw, %, СКО случайной составляющей σ_1w, %;

для второго - значения систематической составляющей δ_2cw, %, СКО случайной составляющей σ_2w, %.

Используют формулу для плотности совместного распределения погрешностей приборов

f(x₁, y₁, x₂, y₂) = f₁^C(x₁) ´ f₁⁰(y₁) ´ f₂^C(x₂) ´ f₂⁰(y₂)                               (В.4)

При этом функции f₁^C(x₁), f_l⁰(y₁), f₂^C(x₂), f₂⁰(y₂) определены следующими формулами

f₁^C(x₁) = δ(x₁ - δ_1cw)                                                           (B.5)

f₂^C(x₂) = δ(x₁ - δ_2cw)                                                           (B.6)

f₁⁰(y₁) - формулой (Б.4), в которой σ_N принимают равной

                                                                   (В.7)

f₂⁰(y₁) - формулой (Б.4), в которой σ_N принимают равной

                                                                  (B.8)

В.3 Рассчитывают вероятность ошибки второго рода

В случае одного негодного прибора используют функцию f(x₁, y₁, x₂, y₂), определенную в п. В.1, а в случае двух негодных приборов - в п. В.2.

Вероятность ошибки второго рода b берут равной интегралу от функции f(x₁, y₁, x₂, y₂) по области G_A, в которой переменные удовлетворяют неравенству

|x₁ + y₁ - x₂ - у₂| £ С                                                             (В.9)

                                            (В.10)

ПРИЛОЖЕНИЕ Г

ПРИМЕР РАСЧЕТА ПО ПРОГРАММЕ «КОНТРОЛЬ ПОГРЕШНОСТИ СЧЕТЧИКОВ ПО РЕЗУЛЬТАТАМ СЛИЧЕНИЙ». МОДЕЛЬНЫЕ ДАННЫЕ

ПРОГРАММА Контроль погрешности счетчиков по результатам сличений

ДАННЫЕ по сличениям. Всего сличений 14

	Показание 1	Показание 2	Расхождение %
1	1607,153	1607,221	-0,004
2	4,000	4,010
3	2576,605	2576,288	0,012
4	110,198	110,193	0,005
5	1381,612	1382,455	-0,061
6	3441,138	3441,491	-0,010
7	1591,410	1591,002	0,026
8	5010,410	5010,300	0,002
9	716,755	716,472	0,039
10	1105,009	1104,316	0,063
11	3000,170	2998,008	0,072
12	3,000	3,010
13	4771,122	4769,016	0,044
14	2825,194	2824,627	0,020

ДАННЫЕ для формирования критерия годности прибора

Метрологические характеристики годного прибора:

граница допускаемых значений погрешности 0,150 %

Вероятность признать годный прибор негодным (ошибка I рода) 0,01

РЕЗУЛЬТАТЫ анализа данных по сличениям

Обнаружен дрейф расхождения показаний приборов.

Модель дрейфа расхождения: А + B ´ Q

Оценки параметров модели (А, В) - 0,01601 0,00251

Стандартные отклонения оценок параметров 0,01618 0,00103

СКО отклонений от модели 0,030

№	Сумма/1000	Расхождение	Модель	Отклонение от модели	Границы полосы		Норм. отклон.
					Нижняя	Верхняя
1	1,607	-0,004	-0,012	0,008	-0,042	0,018	0,3
2	4,188	0,012	-0,006	0,017	-0,031	0,020	0,6
3	4,298	0,005	-0,005	0,010	-0,031	0,020	0,3
4	5,680	-0,061	-0,002	-0,059	-0,025	0,021	-2,0
5	9,121	-0,010	0,007	-0,017	-0,012	0,026	-0,6
6	10,712	0,026	0,011	0,015	-0,007	0,029	0,5
7	15,723	0,002	0,023	-0,021	0,005	0,041	-0,7
8	16,439	0,039	0,025	0,014	0,007	0,044	0,5
9	17,544	0,063	0,028	0,035	0,009	0,047	1,2
10	20,544	0,072	0,035	0,037	0,013	0,058	1,2
11	25,319	0,044	0,047	-0,004	0,017	0,078	-0,1
12	28,144	0,020	0,055	-0,034	0,019	0,090	-1,2

РЕЗУЛЬТАТЫ КОНТРОЛЯ

Критическое значение:

для заданной вероятности ошибки I рода, равной 0,01,

предел допускаемого расхождения показаний двух приборов 0,180 %

Заключение. Результат контроля положительный.

Оперативная характеристика критерия, соответствующая критическому значению, равному 0,180 %:

При одном годном приборе и заданных значениях характеристик

погрешности негодного прибора

(систематическая сост. 0,150, СКО случайной сост. 0,030)

вероятность признать негодный прибор годным (ошибка II рода) равна 0,67

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

(первый - систематическая сост. 0,150, СКО случайной сост. 0,030)

(второй - систематическая сост. - 0,100, СКО случайной сост. 0,040)

вероятность признать их годными (ошибка II рода) равна 0,01

Превышение предела допускаемого расхождения ожидается при поступлении дополнительно 50 тыс. т (м³) в режиме проведения сличений.