РОССИЙСКОЕ АКЦИОНЕРНОЕ ОБЩЕСТВО ЭНЕРГЕТИКИ И ЭЛЕКТРИФИКАЦИИ «ЕЭС РОССИИ»

 

МЕТОДИЧЕСКИЕ УКАЗАНИЯ.
ИНФОРМАЦИОННО-ИЗМЕРИТЕЛЬНЫЕ СИСТЕМЫ.
МЕТОДИКА ОБРАБОТКИ
ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ
МЕТРОЛОГИЧЕСКОЙ АТТЕСТАЦИИ

 

РД 34.11.206-94

 

 

СЛУЖБА ПЕРЕДОВОГО ОПЫТА ОРГРЭС

Москва 1995

 

РАЗРАБОТАНО Акционерным обществом «Фирма по наладке, совершенствованию технологии и эксплуатации электростанций и сетей ОРГРЭС»

ИСПОЛНИТЕЛИ А.Г. АЖИКИН, С.А. СПОРЫХИН

УТВЕРЖДЕНО Департаментом науки и техники РАО «ЕЭС России» 14.01.94

Первый заместитель начальника А.П. БЕРСЕНЕВ

 

Срок действия установлен

с 01.01.95 г.

до 01.01.99 г.

Настоящие Методические указания (МУ) определяют основные положения и устанавливают методы обработки результатов многократных наблюдений, полученных при проведении метрологической аттестации (МА) в соответствии с «РД 34.11.202-87. Методические указания. Информационно-измерительные системы. Метрологическая аттестация. Организация и порядок проведения» (М.: СПО Союзтехэнерго, 1988) измерительных каналов (ИК) информационно-измерительных систем (ИИС), эксплуатируемых на предприятиях электроэнергетической отрасли.

Методические указания предназначены для персонала организаций, проводящих МА ИИС.

В МУ изложены алгоритмы обработки результатов многократных наблюдений при экспериментальных исследованиях метрологических характеристик (MX) ИК ИИС (в условиях эксплуатации) комплектным и поэлементным методами.

Указанные в МУ алгоритмы обработки экспериментальных данных основаны на следующих допущениях:

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

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

закон распределения систематической составляющей погрешности на множестве однотипных ИК является нормальным;

взаимное влияние ИК не имеет места.

С выходом настоящих МУ утрачивают силу «Методические указания. Информационно-измерительные системы. Методика обработки экспериментальных данных метрологической аттестации: РД 34.11.206-88» (М.: СПО Союзтехэнерго, 1988).

1. ОБЩИЕ ПОЛОЖЕНИЯ

1.1. Результаты многократных наблюдений, полученные при экспериментальных исследованиях MX ИК, должны быть представлены в протоколе по форме, приведенной в РД 34.11.202-87.

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

1.2. Результатами обработки экспериментальных данных являются следующие MX:

систематическая составляющая погрешности ИК Δ_s или (и) математическое ожидание систематической составляющей погрешности совокупности однотипных ИК М(Δ_s);

среднее квадратическое отклонение (СКО) систематической составляющей погрешности ИК совокупности однотипных ИК σ(Δ_s);

среднее квадратическое отклонение случайной составляющей погрешности ИК σ() или (и) совокупности однотипных ИК σ();

случайная составляющая погрешности ИК или (и) совокупности однотипных ИК, обусловленная гистерезисом - вариация Н или (и) ;

верхняя и нижняя границы доверительного интервала погрешности измерений ИК или совокупности однотипных ИК Δ_h, Δ_l.

1.3. Результаты МА должны быть представлены в соответствии с МИ 1317-86 в виде Δ_h, Δ_l, Р (Р - доверительная вероятность).

Метрологические характеристики нормируются для групп однотипных ИК или для каждого ИК.

2. АЛГОРИТМ ОБРАБОТКИ РЕЗУЛЬТАТОВ НАБЛЮДЕНИЙ ПРИ КОМПЛЕКТНОМ МЕТОДЕ ИССЛЕДОВАНИЙ MX ИК

2.1. Обработка результатов проводится в три этапа.

На первом этапе полученная совокупность результатов наблюдений для каждого l-го ИК, в каждой j-й точке диапазона измерений проверяется на наличие грубых ошибок (погрешностей). Отдельно проверяются результаты, полученные при прямом и обратном ходах в процессе проведения эксперимента.

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

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

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

Если причина резкого отклонения результата наблюдений не установлена, то следует воспользоваться методом статистической обработки (п. 2.1.2).

2.1.2. Оценивается анормальность результатов наблюдений.

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

                                                     (1)

где x_ljiм(б) - i-e значение результата наблюдений в j-й точке диапазона измерений l-го ИК при прямом (обратном) ходе (i = 1 ... n);

n_м(б) - количество наблюдений при прямом (обратном) ходе.

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

                                         (2)

2.1.2.3. Определяется отношение U по формуле

                                                     (3)

где x*_ljiм(б) - результат наблюдения, по предположению, содержащий грубую погрешность.

2.1.2.4. Полученное значение U сравнивается со значением h по таблице приложения 1.

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

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

Определяется среднее арифметическое значение по формуле (1). Определяется модуль отклонения каждого результата наблюдений от среднего значения

Δx_ljiм(б) = ||.                                                 (4)

Полученные значения отклонений записываются в возрастающий ряд

Δx_lj1м(б) £ ... £ Δx_ljiм(б) £ ... £ Δx_ljnм(б).                                         (5)

Определяется отношение К по формуле

                                                  (6)

Полученное значение К сравнивается со значениями К_д из таблицы.

Если К > К_д, то результат наблюдений, у которого модуль отклонения от среднего значения является наибольшим, исключается из рассматриваемой совокупности и производится его замена, как указано в п. 2.1.2.5.

2.2. На втором этапе определяются следующие характеристики погрешности ИК.

2.2.1. Оценка систематической составляющей погрешности в j-й точке диапазона измерений l-го ИК по формуле

                                                   (7)

где  - среднее арифметическое значение результатов наблюдений в j-й точке диапазона измерений l-го ИК при прямом (обратном) ходе;

x_дj - действительное значение измеряем ой величины в j-й точке диапазона измерений l-го ИК (действительное значение образцового сигнала на входе ИК в единицах измеряемой величины).

Значения  определяются по формулам:

                                                           (8)

                                                          (9)

2.2.2. Оценка дисперсии случайной составляющей погрешности в j-й точке диапазона измерений l-го ИК по формуле

                     (10)

2.2.3. Оценка СКО случайной составляющей погрешности в i-й точке диапазона измерений l-го ИК по формуле

                                                     (11)

2.2.4. Оценка вариаций в j-й точке диапазона измерений l-го ИК по формуле

                                                        (12)

2.2.5. Если обрабатываются результаты наблюдений при исследовании совокупности однотипных ИК, то обработка продолжается и определяются следующие обобщенные характеристики погрешности совокупности ИК.

2.2.5.1. Оценка математического ожидания систематической составляющей погрешности в j-й точке диапазона измерений совокупности однотипных ИК по формуле

                                                  (13)

где m - количество однотипных ИК в обрабатываемой (рассматриваемой) совокупности.

2.2.5.2. Оценка СКО систематической составляющей погрешности в j-й точке диапазона измерений совокупности однотипных ИК по формуле

                                   (14)

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

                                               (15)

Если же дисперсии не однородны, то в качестве  принимается максимальное из т значений , т.е.

                                                 (16)

Алгоритм проверки однородности дисперсий представлен в приложении 2.

2.2.5.4. Оценка вариации в j-й точке диапазона измерений совокупности однотипных ИК в случае допустимости ее расхождения по формуле

                                                       (17)

Если же расхождения значений вариации не допустимо, то в качестве  принимается максимальное из их значений , т.е.

                                                      (18)

Алгоритм проверки допустимости расхождения вариации представлен в приложении 3.

2.3. На третьем этапе обработки в каждой j-й точке диапазона измерений ИК (совокупности однотипных ИК) определяются верхняя (Δ_hl) и нижняя (Δ_lj) границы доверительного интервала погрешности измерений.

2.3.1. Для одного ИК Δ_hl, Δ_lj определяются по формуле

                                             (19)

где t_р - коэффициент Стьюдента, который определяется по таблице приложения 4 в зависимости от доверительной вероятности и количества наблюдений.

2.3.2. Для совокупности однотипных ИК Δ_hl, Δ_lj определяются по формуле

                             (20)

2.4. Полученные значения Δ_hl и Δ_lj записываются как функции Δ_hl = f(x_дj) и Δ_lj = f(x_дj), аппроксимируются по диапазону измерений прямыми (например, методом наименьших квадратов) и представляются в виде:

Δ_h = a₁ + b₁x_изм;                                                       (21)

Δ_l = a₂ + b₂x_изм,                                                       (22)

где a₁(a₂), b₁(b₂) - параметры аппроксимирующих прямых;

х_изм - измеренное значение физической величины.

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

                                              (23)

где q - количество точек диапазона измерений, в которых проводились экспериментальные исследования (j = 1 ... q).

2.5. В качестве нормированных MX ИК (совокупности ИК) во всем диапазоне измерений принимаются приведенные γ_h, γ_l и относительные δ_h, δ_l значения верхней и нижней границ доверительного интервала погрешности измерений, которые определяются по формулам:

;                                                      (24)

;                                                       (25)

δ_h = b₁ · 100 %;                                                           (26)

δ_l = b₂ · 100 %;                                                           (27)

где х_в, х_н - верхний и нижний пределы диапазона измерений.

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

                                                    (28)

                                                    (29)

3. АЛГОРИТМ ОБРАБОТКИ РЕЗУЛЬТАТОВ НАБЛЮДЕНИЙ ПРИ ПОЭЛЕМЕНТНОМ МЕТОДЕ ИССЛЕДОВАНИЯ МХ ИК

3.1. При поэлементном методе исследуются МХ элементов ИК. Обработка результатов наблюдений при исследовании MX элемента ИК проводится согласно пп. 2.1 - 2.3.

3.2. Границы доверительного интервала погрешности измерений определяются следующим образом.

Для одного ИК по формуле

                                 (30)

где  - оценка систематической составляющей погрешности в j-й точке диапазона измерений k-го элемента l-го ИК;

 - оценка СКО случайной составляющей погрешности в j-й точке диапазона измерений k-го элемента l-го ИК;

 - оценка вариации в j-й точке диапазона измерений k-го элемента l-го ИК;

z - количество элементов ИК.

Для совокупности однотипных ИК по формуле

                 (31)

где  - оценка математического ожидания систематической составляющей погрешности в j-й точке диапазона измерений k-го элемента совокупности однотипных ИК;

 - оценка СКО систематической составляющей погрешности в j-й точке диапазона измерений k-го элемента совокупности однотипных ИК;

 - оценка СКО случайной составляющей погрешности в j-й точке диапазона измерений k-го элемента совокупности однотипных ИК;

 - оценка вариации в j-й точке диапазона измерений k-го элемента совокупности однотипных ИК.

3.3. Если поэлементный метод исследования MX ИК является расчетно-экспериментальным, то расчет MX (составляющих погрешности в каждой j-й точке диапазона измерений) элементов ИК, не исследуемых экспериментально, следует проводить по MX агрегатных средств измерений (АСИ), входящих в состав элемента в соответствии с «Методикой определения обобщенных метрологических характеристик измерительных каналов ИИС и АСУ ТП по метрологическим характеристикам агрегатных средств измерений: МТ 34-70-038-87» (M.: СПО Союзтехэнерго, 1987).

3.4. Нормирование MX ИК (совокупности ИК) проводится согласно пп. 2.4 и 2.5.

Приложение 1

ПРЕДЕЛЬНЫЕ ЗНАЧЕНИЯ h ПРИ ОЦЕНКЕ РЕЗУЛЬТАТОВ НАБЛЮДЕНИЙ НА АНОРМАЛЬНОСТЬ

Приложение 2

ПРОВЕРКА ОДНОРОДНОСТИ m ОЦЕНОК ДИСПЕРСИЙ

1. При m = 2 оценки дисперсий считаются однородными, если выполняется неравенство

Значение F находится по табл. П2.1 в зависимости от доверительной вероятности и параметров f₁ и f₂, которые определяются по формуле

f₁ = f₂ = n_м + n_б - 1.

Таблица П2.1

Значения F

2. При т > 2 оценки дисперсий считаются однородными, если выполняется неравенство

G_max £ G.

Значение G_max находится по формуле

Значения G находятся по табл. П2.2 в зависимости от числа наблюдений при доверительной вероятности 0,95.

Таблица П2.2

Значения G

Приложение 3

ПРОВЕРКА ДОПУСТИМОСТИ РАСХОЖДЕНИЯ т ОЦЕНОК ВАРИАЦИИ

1. Определяется оценка межгрупповой дисперсии по формуле

2. Определяется среднее значение внутригрупповых оценок дисперсии по формуле

3. Расхождение вариации считается допустимым, если выполняется неравенство

Значение F находится по табл. П2.1 в зависимости от доверительной вероятности и параметров f₁ и f₂, которые определяются по формулам:

f₁ = m - 1;

f₂ = m(n_м + n_б - 1).

Приложение 4

ЗНАЧЕНИЕ КОЭФФИЦИЕНТА СТЬЮДЕНТА

 

СОДЕРЖАНИЕ