Страница 2 из 35

Где подвох?

Добавлено: 24 дек 2017, 18:45
Токарев Виктор
AGL » 21 минуту назад писал(а):
Я в двух картинках разницы не вижу... И где тут подвох?
Разницы в уравнениях нет, а в линейных трендах есть, моя линия Мут разместилась на уровне 6 т по оси М2, а Ваша на ровне 1 т по оси М2, почему это произошло не знаю, ведь уравнения при этом одинаковые.

Находим показатели функционирования узла учета.

Добавлено: 24 дек 2017, 18:54
AGL
Теперь понятно, где подвох. Ось М2 не имеет отношения к функции Мут = f(М1). Пересечение линии тренда этой функции надо смотреть на оси "Мут", а не на оси "М2".
А расхождение в точках пересечения функции Мут = f(М1) с осями "М2" вызвано различиями в масштабах осей "М2" и "Мут". Если масштабы всех осей (левой и правой) сделать такими же, как и на моей картинке, то подвох тут же исчезнет.

Находим показатели функционирования узла учета.

Добавлено: 24 дек 2017, 18:58
Токарев Виктор
Все верно при построении Мут задействована ось Х = М1 и вспомогательная ось У = Мут, почему я привязал ее к М2 не знаю, ведь мозг то говорил что нужно смотреть на другую ось.

Находим показатели функционирования узла учета.

Добавлено: 24 дек 2017, 19:07
AGL
Ради интереса можно масштаб оси "М2" сделать 0...45, а "Мут" - 0...-0,9. И две картинки будут, как близнецы-братья :-):

Находим показатели функционирования узла учета.

Добавлено: 24 дек 2017, 19:43
Токарев Виктор
AGL » 35 минут назад писал(а):Ради интереса можно масштаб оси "М2" сделать 0...45, а "Мут" - 0...-0,9. И две картинки будут, как близнецы-братья :-):
2.png
Так оно и есть, что-то я заплутал маленько, ну ничего, главное во время поправиться.

Находим показатели функционирования узла учета.

Добавлено: 24 дек 2017, 19:56
AGL
Теперь всё встало на свои места - две картинки почти совпали. "Почти" - потому что ось "У" получилась короткая (0...45), а ось "Х" - длинная (0...50). В результате появилось пустое место на диаграмме справа. Желательно (если нет противопоказаний), чтобы информативная часть изображаемых функций занимала максимум полезной площади области построения диаграммы.

Открыто-закрытая система ГВС

Добавлено: 27 дек 2017, 01:36
AGL
Теперь мы можем оценить качество измерений в системе ГВС.
На прищепке - часовой архив за ноябрь-декабрь из открыто-закрытой системы ГВС детского садика. Система названа открыто-закрытой потому, что по рабочим дням в рабочее время система открытая, а по ночам и в выходные дни система ГВС становится закрытой.

Находим показатели функционирования узла учета.

Добавлено: 27 дек 2017, 16:47
AGL
По данным архива-4 строим традиционные временные графики М1гвс, М2гвс и Мгвс:
1093_М1гвс, М2гвс, Мгвс.png
Здесь видно, что в часы максимального водоразбора двухтрубная система ГВС превращается в однотрубную - периодически в линии циркуляции М2гвс = 0, т.к. вся вода, поступившая в систему ГВС, расходуется на водоразбор и возвращать в теплоцентр уже нечего.
Также видно, что по ночам и в выходные дни имеется утечка Мут = М1 - М2. Что это за утечка - физическая или метрологическая - нам предстоит выяснить. [/color]

Изменение Мгвс с течением времени

Добавлено: 27 дек 2017, 17:10
AGL
Так выглядит измеренная утечка по ночам и в выходные дни:
1094_Мгвс по ночам.png
Эти "горные ландшафты" изменяются в диапазоне от 20 до 100 кг/ч. На фоне этих "ландшафтов" заметен небольшой водоразбор в выходные дни. Видимо, охрана садика изредка пользуется горячей водой. [/color]

Зависимость "утечки" от М1гвс

Добавлено: 27 дек 2017, 18:19
AGL
На этой картинке показана зависимость Мгвс = М1гвс - М2гвс от М1гвс в области наименьших значений Мгвс.
1095_Мгвс = f(М1гвс).png
Видим, что линия нулевого водоразбора в ночные часы и выходные дни расположена под большим углом к оси "Х": чем больше часовые массы М1гвс, тем больше измеренная "утечка" Мут = М1 - М2. Такая взаимосвязь Мгвс и М1гвс указывает на наличие положительного мультипликативного рассогласования погрешностей двух расходомеров (наклон ГХ расходомера М2гвс меньше наклона ГХ расходомера М1гвс). Именно по этой причине "утечки" ночью и в выходные дни строго следят за циркрасходом.
Насколько велико мультипликативное рассогласование погрешностей dМ1гвс и dМ2гвс? Попробуем это выяснить...[/color]

Распределение Мгвс по часам суток

Добавлено: 27 дек 2017, 18:58
AGL
1096_Мгвс по часам суток.png
Из этой диаграммы следует, что горячеводная жизнь в садике начинается на 8-м часу суток и заканчивается на 21-м. Всё остальное время система ГВС является закрытой. Правда, по ночам и в выходные охранники изредка открывают кран на несколько секунд.
Для получения уравнения зависимости утечки Мут = f(М1гвс) из общего архива отберём данные часов с 22-го по 7-й и построим функцию Мут = f(М1гвс).[/color]

Зависимость ночной "утечки" от М1

Добавлено: 27 дек 2017, 20:31
AGL
Получаем вот такую зависимость метрологической утечки от часовых масс М1:
1097_Мут = f(М1гвс).png
Несколько точек на этой диаграмме выпадают из общего стройного ряда. Скорее всего, это небольшие водоразборы увеличивают "утечку" (охраняющий садик человек помыл руки). Чтобы не вносить дополнительную погрешность в статистическое уравнение Мут = f(М1), эти нехарактерные точки нужно исключить из массива данных.[/color]

Окончательные функции Мут = f(M1) и М2 = f(M1)

Добавлено: 28 дек 2017, 00:30
AGL
После удаления подозрительных значений Мут получаем окончательные зависимости Мут = f(M1) и М2 = f(М1) для тех часов, в течение которых система ГВС являлась закрытой.
1098_Мут и М2 от М1.png
Теперь мы для каждого часа из исходного архива-4 можем рассчитать метрологическую утечку по формуле
Мут = 0,02848*М1 + 0,0220, т за час. (1)
А формула
М2 = 0,9715*М1 - 0,0220, т за час (2)
нам понадобится для расчёта показателей функционирования данной пары расходомеров.[/color]

Графики Мгвс и Мут

Добавлено: 28 дек 2017, 01:12
AGL
Рассчитав для каждого архивного часа "утечку" по формуле (1), можно добавить значения Мут на график изменения Мгвс.
1099_Мгвс и Мут.png
Видим, что в часы, когда потребления воды нет, график Мут совпадает с графиком Мгвс. А благодаря формуле (1) мы имеем точные значения Мут и для тех часовых интервалов, когда система ГВС является открытой.
Всего за 1248 архивных часов метрологическая утечка составила +64 тонны (в среднем +51 кг/ч), что привело к завышению результата учёта Мгвс на +30,0 %.[/color]

Находим показатели функционирования узла учета.

Добавлено: 28 дек 2017, 16:25
AGL
Имея результаты измерений "утечки" в ночные часы в закрытой системе ГВС, рассчитаем относительную (в процентах по отношению к М1гвс) утечку dМут, которая в закрытой системе представляет из себя разность относительных погрешностей расходомеров М1 и М2. Получаем такую картину:
1100_dМут = f(M1).png
Здесь аддитивная и мультипликативная составляющие расхождений ГХ двух расходомеров имеют положительный знак, из-за чего хотя бы частичная взаимокомпенсация этих двух составляющих расхождения показаний М1 и М2 не происходит. Поэтому нет никакой надежды на то, что расхождения погрешностей dM1 и dM2 хотя бы на наибольших расходах (близких к Qmax = 12 т/ч) не превысят допуска +/-1,4 %.
Эта статистическая функция dМут = f(М1гвс) указывает на то, что фактический диапазон измерений расхода у этой пары равен нулю при нормируемом диапазоне 1:150. [/color]

Закрытые часы из открытой системы ГВС

Добавлено: 30 дек 2017, 21:39
AGL
В прилагаемом архиве-4 исключены все часы, в течение которых система ГВС была открытой. Оставлены только данные, накопленные в ночные часы, когда водоразбор в дет. садике отсутствует. Для получившейся закрытой системы можно провести оценку показателей работы пары расходомеров так, как это было подробно проделано с архивом-2.

Показатели для системы ГВС (архив-4)

Добавлено: 30 дек 2017, 22:34
AGL
По данным "закрытого" архива-4 (система ГВС в дет. садике) получаем следующие "показатели функционирования" пары расходомеров:
1101_Архив-4. Показатели.png
В общем, всё (кроме критерия R2) получилось не очень хорошо при нулевом диапазоне измерений. Надо бы, конечно, сдать эту пару изделий в реанимацию. Но внеочередно поверять у нас не принято, а очередная (плановая) поверка ещё не скоро. Так что придётся садику потерпеть. Тем более, что 30-процентное завышение учёта (на радость продавцу) оплачивает какой-то там бюджет/комитет, а не директор садика из собственной зарплаты. [/color]

Находим показатели функционирования узла учета.

Добавлено: 03 янв 2018, 21:22
Токарев Виктор
AGL » 22 дек 2017, 08:25 писал(а): А здесь мы видим график изменения случайных относительных расхождений часовых масс М1 и М2, выраженных в процентах от М2. Взаимная относительная нестабильность М1 и М2 (ВзНо), как и случайная погрешность, характеризуется среднеквадратическим отклонением (СКО).
Из графика следует, что по результатам учёта за 744 часа СКО ВзНо = +/-0,446 % при допуске +/-0,18 % (для 1-процентных расходомеров). Превышение допуска - в 2,5 раза, и по этому показателю результат тестирования получился "НЕ ГОДЕН".
Коль скоро взаимная "болтанка" часовых масс М1 и М2 существенна, у расходомеров данного типа нужно нормировать как систематическую, так и случайную погрешность. Но о необходимости нормирования случайной составляющей погрешности давно уже никто не вспоминает... Видимо, считается, что случайная составляющая несущественна и потому не стоит о ней вспоминать.
Александр Григорьевич, а как найти Взаимную относительную нестабильность ВзНо и среднеквадратическое отклонение СКО?

Находим взаимную нестабильность ВзНо

Добавлено: 03 янв 2018, 22:15
AGL
Эту задачу будем решать в такой последовательности.
1. Ранее (см. сообщ. № 19) мы для архива-2 получили среднюю линейную функцию М2лин = f(М1) следующего вида:
М2лин = 0,985*М1 + 0,883, т за час.
2. Создаём в исходном архиве-2 столбец с заголовком "М2лин" и в каждой ячейке этого столбца рассчитываем значения М2лин по вышеозначенной формуле, подставляя в формулу архивные часовые массы М1.
3. В соседнем столбце с заголовком "М2сл" находим абсолютные отклонения измеренных часовых масс М2 от М2лин, т.е абсолютную взаимную нестабильность пары расходомеров ВзНа: ВзНа = М2 - М2лин.
4. В новом столбце с заголовком "ВзНо" (Взаимная Нестабильность относительная) в каждой ячейке рассчитываем значения ВзНо по формуле
ВзНо = (ВзНа/М2лин)*100 %.
ВременнОй график ВзНо показан в сообщении № 25.
На этом графике также показаны границы допускаемого СКО ВзНод = +/-0,18 % и границы фактического СКО ВзНо = +/-0,446 %.

Расчет СКО ВзНо

Добавлено: 03 янв 2018, 22:50
AGL
Для расчёта среднеквадратического отклонения (СКО) взаимной относительной нестабильности (ВзНо) работы пары расходомеров следует воспользоваться стандартной функцией "СТАНДОТКЛОНП", встроенной в Excel.
1104_Расчет СКО ВзНо.png
В нашем случае СКО ВзНо = +/-0,446 % при допуске +/-0,18 %. Т.е. случайная болтанка измеренных часовых масс М2 относительно М1 в 2,5 раза превысила допуск для стабильно работающих расходомеров.
Одной из причин повышенной взаимной болтанки М1 и М2 является известная болезнь, именуемая как температуробоязнь...[/color]