Анализ почасовых архивов

AGL
VIP персона
Posts in topic: 32
Сообщения: 3771
Зарегистрирован: 17 окт 2011, 14:43
Благодарил (а): 14 раз
Поблагодарили: 19 раз

Анализ почасовых архивов

#41

Непрочитанное сообщение AGL » 20 апр 2014, 03:17

Adrenalin Splash » Вчера, 22:47 писал(а):как Вы это увидели: "...тут не работает (или вовсе отсутствует) регулятор температуры воды, подаваемой в систему ГВС. Т.е. подмес обратной (более холодной) воды в систему ГВС не производится...".
Это из-за вот этих скачков между d12 и d34, которые отмечены стрелками на рисунке, или не в этом дело?
Если бы в данном теплоцентре регулятор температуры ГВС был в наличии и исправно функционировал, то температура в подающем трубопроводе системы ГВС t3 (Т1гвс) была бы близка к константе на уровне, близком к заданном в регуляторе. Например, t3 = Т1гвс = 60 - 65 °С, независимо от значений Т1. Здесь же t3 (Т1гвс) практически равна Т1. Т.е. вся вода, подаваемая на вход системы ГВС, подаётся только из подающего трубопровода (см. схему измерений № 3 download/file.php?id=5016" onclick="window.open(this.href);return false;), а из обратного трубопровода вода в систему ГВС не поступает.
Можно ожидать, что в сильные морозы Т1 достигнет 110 - 115 °С, и в систему ГВС будет подана перегретая вода, что смертельно опасно и потому категорически недопустимо. Следовательно, надзорные органы (например, жилинспекция) обязаны выдать виновным грозное предписание на немедленную установку регулятора Т1гвс. В противном случае (не дай Бог, конечно!) рано или поздно в какой-то квартире случится разрыв трубы, перегретая вода хлынет на кухню и мгновенно превратится в пар, и жильцы тут же сварятся на высокотемпературном пару.

Скачки перепадов температур dt12 и dt34 никак не связаны с отсутствием регулирования t3 (Т1гвс), а обусловлены изменением напора на входе систем отопления и ГВС. Видимо, кто-то в теплоцентре периодически крутит задвижку. При этом изменяется напор (разность давлений DP12 = P1 - P2) и, как следствие, изменяются циркрасходы в системах отопления и вентиляции. Уменьшение расхода приводит к уменьшению Т2 (Т2гвс) и росту dT12 и dT34. И наоборот.

Аватара пользователя
Токарев Виктор
Ветеран форума
Posts in topic: 19
Сообщения: 1245
Зарегистрирован: 23 июн 2012, 05:21
Благодарил (а): 3 раза
Поблагодарили: 15 раз

Анализ почасовых архивов

#42

Непрочитанное сообщение Токарев Виктор » 20 апр 2014, 11:33

AGL » Вчера, 11:11 писал(а): но не хватает важных графиков двух разностей масс G1 - G2 и G3 - G4.
Согласен М1-М2 и М3-М4, это важные элементы анализа, по этому я их предпочел для наглядности отсортировать по возрастанию М1-М2, и посмотреть что вышло.
А вышло следующее:
Вложения
2 график.png

Adrenalin Splash
Начинающий
Posts in topic: 18
Сообщения: 18
Зарегистрирован: 11 апр 2014, 23:25

Анализ почасовых архивов

#43

Непрочитанное сообщение Adrenalin Splash » 26 апр 2014, 17:02

Извините, что вsпал из обсуждений. Работа. Я подредактировал графики и сделал покрасивше для глаз! Ну по крайней мере постарался так сделать, думаю - получилось. Какие шаги нужно проделать дальше? Виктор, какой тип диаграммы нужно выбрать, чтобы получилась диаграмма как у вас? Все перепробовал, не нашел такую. Пологий график не выходит...

AGL
VIP персона
Posts in topic: 32
Сообщения: 3771
Зарегистрирован: 17 окт 2011, 14:43
Благодарил (а): 14 раз
Поблагодарили: 19 раз

Анализ почасовых архивов

#44

Непрочитанное сообщение AGL » 26 апр 2014, 17:55

Если построить графики изменения во времени разностей масс dM12 = М1 - М2 и Мгвс = М3 - М4, то получим следующее:
749_4-канальный ТС (6).png
Такие графики хороши для оценки динамики изменения разностей масс с течением времени.
Но эти графики не очень наглядны с точки зрения оценки соотношения часовых величин dM12 и Мгвс. Поэтому иногда для таких оценок необходимо отсортировать таблицу архива по возрастанию (или по убыванию) архивных значений одного из параметров.
В данном случае Виктор отсортировал архив по возрастанию значений dM12. В результате график dM12 стал монотонно возрастающим (каждое последующее значение dM12 равно или больше предыдущего значения), и на фоне возрастающего графика dM12 хорошо видно, как изменяются значения разности масс Мгвс.[/color]

AGL
VIP персона
Posts in topic: 32
Сообщения: 3771
Зарегистрирован: 17 окт 2011, 14:43
Благодарил (а): 14 раз
Поблагодарили: 19 раз

Анализ почасовых архивов

#45

Непрочитанное сообщение AGL » 26 апр 2014, 18:26

Здесь исходная таблица архива отсортирована по убыванию часовых разностей масс Мгвс = М3 - М4.
750_4-канальный ТС (7).png
Видим, что пурпурный график Мгвс монотонно убывает, а зелёный график разности масс dM12 = М1 - М2 расположился чуть выше графика Мгвс. При этом наблюдается заметная "вибрация" dM12 по отношению к Мгвс (случайные расхождения dM12 и Мгвс на часовых интервалах времени довольно существенны).
В дальнейшем мы получим численные характеристики систематического и случайного расхождения измеренных разностей масс dM12 и Мгвс и сравним эти характеристики с метрологическими допусками. В конечном итоге мы вынесем вердикт "годен - не годен" для этого теплосчётчика.[/color]

Аватара пользователя
Токарев Виктор
Ветеран форума
Posts in topic: 19
Сообщения: 1245
Зарегистрирован: 23 июн 2012, 05:21
Благодарил (а): 3 раза
Поблагодарили: 15 раз

Анализ почасовых архивов

#46

Непрочитанное сообщение Токарев Виктор » 26 апр 2014, 18:27

Adrenalin Splash » Сегодня, 10:02 писал(а):Извините, что вsпал из обсуждений. Работа. Я подредактировал графики и сделал покрасивше для глаз! Ну по крайней мере постарался так сделать, думаю - получилось. Какие шаги нужно проделать дальше? Виктор, какой тип диаграммы нужно выбрать, чтобы получилась диаграмма как у вас? Все перепробовал, не нашел такую. Пологий график не выходит...
Простой график. Выше Александр Григорьевич пояснил, как я его построил.
Работа у всех, тут дело не в работе, тут больше в желании все заложено, как я понимаю у Вас нет рвения чего-то изучать.
Кстати, кем Вы работаете? Представитель ТСО или обслуживающая организация? Может кем другим, для чего вообще заикнулись, что хотите прикоснуться к анализу почасовок, может не стоит время тратить то?

Аватара пользователя
Токарев Виктор
Ветеран форума
Posts in topic: 19
Сообщения: 1245
Зарегистрирован: 23 июн 2012, 05:21
Благодарил (а): 3 раза
Поблагодарили: 15 раз

Анализ почасовых архивов

#47

Непрочитанное сообщение Токарев Виктор » 26 апр 2014, 18:41

AGL » 13 минут назад писал(а):Здесь исходная таблица архива отсортирована по убыванию часовых разностей масс Мгвс = М3 - М4.
Вспомнилось обсуждение схождения масс на одном из источников теплоснабжения моего города, надо сказать дела в обсуждаемом архиве немного лучше чем тут: viewtopic.php?p=3693#p3693" onclick="window.open(this.href);return false;

AGL
VIP персона
Posts in topic: 32
Сообщения: 3771
Зарегистрирован: 17 окт 2011, 14:43
Благодарил (а): 14 раз
Поблагодарили: 19 раз

Анализ почасовых архивов

#48

Непрочитанное сообщение AGL » 26 апр 2014, 18:49

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

Аватара пользователя
Токарев Виктор
Ветеран форума
Posts in topic: 19
Сообщения: 1245
Зарегистрирован: 23 июн 2012, 05:21
Благодарил (а): 3 раза
Поблагодарили: 15 раз

Анализ почасовых архивов

#49

Непрочитанное сообщение Токарев Виктор » 26 апр 2014, 18:56

AGL » 5 минут назад писал(а):Всем, кто интересуется таким анализом, рекомендую посетить эту тему и ознакомиться с её содержимым.
Это точно, там путь анализа проложен во всех красках, хотя все зависит от цели анализа, что нас смущает и т.п.

AGL
VIP персона
Posts in topic: 32
Сообщения: 3771
Зарегистрирован: 17 окт 2011, 14:43
Благодарил (а): 14 раз
Поблагодарили: 19 раз

Анализ почасовых архивов

#50

Непрочитанное сообщение AGL » 26 апр 2014, 19:13

На картинке - фрагмент исходного архива после сортировки строк по убыванию значений столбца "Н" (в этом столбце размещены значения Мгвс).
751_4-канальный ТС (8).png
Чтобы выполнить такую сортировку, необходимо поставить курсор в заголовок столбца и нажать инструмент "Я - А" (т.е. сортировать от бОльшего к меньшему).
После такой сортировки строки архива перемешались так, что значения Мгвс стали убывать от бОльшего к меньшему.
Именно после такой сортировки строк построены графики Мгвс и dM12 на последней картинке.
Чтобы вернуть таблицу архива в исходное состояние, надо поставить курсор в заголовок столбца "А" (№ строки) и нажать кнопку "А - Я" (сортировать по возрастанию). В результате все строки вернутся на свои места. [/color]

AGL
VIP персона
Posts in topic: 32
Сообщения: 3771
Зарегистрирован: 17 окт 2011, 14:43
Благодарил (а): 14 раз
Поблагодарили: 19 раз

Анализ почасовых архивов

#51

Непрочитанное сообщение AGL » 26 апр 2014, 21:33

По данным исследуемого архива построим точечную диаграмму, связывающую измеренные разности масс Мгвс и dM12.
752_4-канальный ТС (9).png
Наносим на функцию Мгвс = f(dM12) линейный тренд и уравнение этого тренда:
Мгвс = 1,0013*dM12 - 0,0147, т за час.
Из уравнения линейного тренда видно, что результаты измерений двух разностей масс имеют как мультипликативное (относительное), так и аддитивное (абсолютное) рассогласования.
Среднее относительное рассогласование составляет 0,13 % (наклон функции Мгвс = f(dM12) на 0,13 % отличается от единицы).
Среднее абсолютное рассогласование (иначе - смещение функции Мгвс = f(dM12)) двух разностей масс составило 14,7 кг/ч.
Статистическая функция Мгвс = f(dM12) имеет заметное рассеивание, из-за чего параметр R2 имеет значение, довольно далеко удалённое от 1 (R2 = 0,984; при стабильных измерениях R2 не менее 0,9998).[/color]

AGL
VIP персона
Posts in topic: 32
Сообщения: 3771
Зарегистрирован: 17 окт 2011, 14:43
Благодарил (а): 14 раз
Поблагодарили: 19 раз

Анализ почасовых архивов

#52

Непрочитанное сообщение AGL » 26 апр 2014, 23:29

Построим временнЫе графики изменения фактической (измеренной) и допускаемой утечек и сравним факт с допуском.
753_4-канальный ТС (10).png
Фактическая утечка Мут (иначе - рассогласование двух разностей масс) рассчитана для каждого часа по формуле
Мут = (М1 - М2) - (М3 - М4) = dM12 - Мгвс.
Допускаемая утечка зависит от допускаемой погрешности измерений каждого их четырёх расходомеров и определяется по формуле
Мутдоп = +/-(DM12 + DM22 + DM32 + DM42)0,5,
где DMi - допускаемая абсолютная погрешность измерения расхода i-м расходомером (i = 1...4).
Допускаемая абсолютная погрешность i-го расходомера определятся по формуле
DMi = 0,01*dMi*Mi,
где dMi - допускаемая относительная погрешность i-го расходомера (в процентах);
Mi - часовая масса, измеренная i-м расходомером.
При расчётах принято, что в данном узле учета применяются 1-процентные расходомеры.
Как следует из графиков, фактическое рассогласование показаний двух пар расходомеров в 7 раз меньше допускаемого рассогласования (в среднем Мут = 14 кг/ч, Мутдоп = +/-101 кг/ч). Следовательно, мы можем говорить о высокой степени согласования всех четырёх каналов измерения расхода, благодаря которому рассогласование (иначе - расхождение) разностей масс dM12 и Мгвс не должно вызывать метрологического беспокойства.[/color]

AGL
VIP персона
Posts in topic: 32
Сообщения: 3771
Зарегистрирован: 17 окт 2011, 14:43
Благодарил (а): 14 раз
Поблагодарили: 19 раз

Анализ почасовых архивов

#53

Непрочитанное сообщение AGL » 27 апр 2014, 00:27

На этой диаграмме значения фактической и допускаемой утечек построены как функции разности масс dM12.
754_4-канальный ТС (11).png
Как видим, и Мут, и Мутдоп почти не зависят от измеренной разности масс dM12 - средние функции Мут = f(dM12) и Мутдоп = f(dM12) почти горизонтальны. Независимость Мут от dM12 стала возможной благодаря незначительной (0,13 %) непропорциональности Мгвс и dM12 (см. download/file.php?id=5107" onclick="window.open(this.href);return false;).[/color]

AGL
VIP персона
Posts in topic: 32
Сообщения: 3771
Зарегистрирован: 17 окт 2011, 14:43
Благодарил (а): 14 раз
Поблагодарили: 19 раз

Анализ почасовых архивов

#54

Непрочитанное сообщение AGL » 28 апр 2014, 02:26

Часто на практике возникает необходимость оценить допускаемую погрешность измерения разности масс на вводе в теплоцентр (в наших обозначениях это dM12 = М1 - М2) и/или разности масс, измеренной в системе ГВС (у нас это Мгвс = М3 - М4).
Рассчитаем допускаемые погрешности измерения двух разностей масс по данным рабочего архива, который можно скачать здесь: download/file.php?id=5090" onclick="window.open(this.href);return false;.
При расчетах будем предполагать, что допускаемая отн. погрешность всех четырёх расходомеров dMi = +/-1 %.
Допускаемую отн. погрешность измерения разности масс dM12 для каждого архивного часа рассчитаем по формуле
ddM12 = +/-[(DdM12/dM12)*100 %],
где DdM12 - допускаемая абсолютная погрешность измерения разности масс dM12, т.
Значения DdM12 для каждого часа рассчитаем по формуле
DdM12 = +/-(DM12 + DM22)0,5.
В этой формуле DM1 и DM2 - допускаемые абсолютные погрешности измерения часовых масс М1 и М2 соответственно.
DM1 = 0,01*dM1*M1, DM2 = 0,01*dM2*M2,
где dM1 = dM2 = +/-1 %.
По результатам расчётов построим временнОй график ddM12.

AGL
VIP персона
Posts in topic: 32
Сообщения: 3771
Зарегистрирован: 17 окт 2011, 14:43
Благодарил (а): 14 раз
Поблагодарили: 19 раз

Анализ почасовых архивов

#55

Непрочитанное сообщение AGL » 28 апр 2014, 02:58

Так выглядит график изменения во времени часовых значений допускаемой относительной погрешности измерения разности масс dM12:
754_4-канальный ТС (12).png
Выглядит этот график, прямо скажем, неважно: в часы активного водопотребления доп. погрешность разности масс составляет десятки процентов, а в ночные и околоночные часы, когда потребление воды снижается до минимума, доп. погрешность dM12 составляет сотни процентов.
За 960 часов работы измеренная разность масс dM12 = 300 т. Эти 300 т измерены 1-процентными расходомерами с допускаемой погрешностью +/-95 т.
Средневзвешенная (за 960 часов) допускаемая относительная погрешность измерения разности масс dM12 = М1 - М2 составила +/-32 %.
На мой взгляд это очень плохо - считать товар (горячую воду) и деньги за этот товар с погрешностью в 32 %. [/color]

Adrenalin Splash
Начинающий
Posts in topic: 18
Сообщения: 18
Зарегистрирован: 11 апр 2014, 23:25

Анализ почасовых архивов

#56

Непрочитанное сообщение Adrenalin Splash » 29 апр 2014, 02:04

Запнулся на построении точечной диаграммы. Диаграмму построить получилось, совсем как у Вас, но не совсем понял куда вводить формулу: Мгвс = 1,0013*dM12 - 0,0147? Это нужно сделать отдельный столбец с расчетами и на основании этого построить диаграмму? И вот момент не понятен - как Вы по графику определили:Среднее абсолютное рассогласование (иначе - смещение функции Мгвс = f(dM12)) двух разностей масс составило 14,7 кг/ч. :smu:sche_nie:
Виктор, моя работа связана с ЖКХ. Помимо анализа работы приборов учета в мою компетенцию входит много различных вопросов. Часто приходится заниматься работой в личное время. Если остальные вопросы я решаю самостоятельно, то анализ для меня - это дебри, в которых без специалистов мне не разобраться и поэтому я обратился сюда. Так как анализ работы узлов учета не является основой моей работы, то я не всегда могу уделять этому столько времени, сколько мне самому хотелось бы. А порой и физически не хватает сил смотреть в монитор от напряжения, да и голова порой дуреет за день и иногда лучше обучение отложить до лучших времен. Надеюсь, что тема не загнется, ведь все что здесь обсуждается не происходит зря! Спасибо, что помогаете разобраться.

AGL
VIP персона
Posts in topic: 32
Сообщения: 3771
Зарегистрирован: 17 окт 2011, 14:43
Благодарил (а): 14 раз
Поблагодарили: 19 раз

Анализ почасовых архивов

#57

Непрочитанное сообщение AGL » 29 апр 2014, 03:25

Adrenalin Splash писал(а):куда вводить формулу: Мгвс = 1,0013*dM12 - 0,0147?
Эту формулу нам выдал Excel по нашей же просьбе. В дальнейшем эта формула нам пригодится для расчёта линейной (линеаризованной) зависимости Мгвс от dM12.
Формула ценна тем, что она показывает нам размер мультипликативного (0,13 %) и аддитивного (0,0147 т за час) рассогласований двух разностей масс - dM12 = M1 - M2 и Мгвс = М3 - М4.
Добывать эту формулу будем в следующей последовательности.
После построения точечной статистической функции Мгвс = f(dM12) ставим курсор на полученную функцию, кликаем правую кнопку мыши и в появившемся меню выбираем строку "добавить линию тренда":
1 Добавить линию тренда.png
[/color]

AGL
VIP персона
Posts in topic: 32
Сообщения: 3771
Зарегистрирован: 17 окт 2011, 14:43
Благодарил (а): 14 раз
Поблагодарили: 19 раз

Анализ почасовых архивов

#58

Непрочитанное сообщение AGL » 29 апр 2014, 03:45

Из шести типов аппроксимации статистической функции Мгвс = f(dM12) мы выбираем линейную аппроксимацию, т.к. по сути своей функция Мгвс = f(dM12) - это есть линейная функция (в идеальном случае Мгвс = dM12). Выбрав линейную аппроксимацию, идем в меню "Параметры".
2 Выбираем линейный тренд.png
[/color]

AGL
VIP персона
Posts in topic: 32
Сообщения: 3771
Зарегистрирован: 17 окт 2011, 14:43
Благодарил (а): 14 раз
Поблагодарили: 19 раз

Анализ почасовых архивов

#59

Непрочитанное сообщение AGL » 29 апр 2014, 04:02

3 Параметры линии тренда.png
В меню "Параметры" вписываем название линейного тренда (аппроксимирущей функции) - "Средняя функция Мгвс = f(dM12)".
Ставим галочку в окошке "показывать уравнение на диаграмме" (имеется ввиду линейное уравнение, описывающее заказанный нами линейный тренд).
Ставим галочку в окошке "поместить на диаграмме величину достоверности аппроксимации R2".
Кликаем ОК и смотрим, что получилось.[/color]

AGL
VIP персона
Posts in topic: 32
Сообщения: 3771
Зарегистрирован: 17 окт 2011, 14:43
Благодарил (а): 14 раз
Поблагодарили: 19 раз

Анализ почасовых архивов

#60

Непрочитанное сообщение AGL » 29 апр 2014, 04:21

А получилось то, что нам и надо:
4 Тренд и его уравнение.png
На диаграмме появился линейный аналог размытой статистической функции и уравнение этой линии: у = 1,0013х - 0,0147, R2 = 0,9842.
Теперь для лучшего восприятия меняем "у" на "Мгвс", "х" на "dM12", форматируем линию тренда и уравнение и получаем диаграмму, которая приведена здесь: download/file.php?id=5107" onclick="window.open(this.href);return false;.[/color]

Ответить