ES2334684T3 - Diagnostico y pronosticos de fallos de sensores que utilizan un modelo de componente y expansiones ortogonales a escala temporal. - Google Patents

Diagnostico y pronosticos de fallos de sensores que utilizan un modelo de componente y expansiones ortogonales a escala temporal. Download PDF

Info

Publication number
ES2334684T3
ES2334684T3 ES05732341T ES05732341T ES2334684T3 ES 2334684 T3 ES2334684 T3 ES 2334684T3 ES 05732341 T ES05732341 T ES 05732341T ES 05732341 T ES05732341 T ES 05732341T ES 2334684 T3 ES2334684 T3 ES 2334684T3
Authority
ES
Spain
Prior art keywords
sensor
compressor
component model
controller
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
ES05732341T
Other languages
English (en)
Inventor
Pengju Kang
Mohsen Farzad
Slaven Stricevic
Payman Sadegh
Alan M. Finn
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Carrier Corp
Original Assignee
Carrier Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Carrier Corp filed Critical Carrier Corp
Application granted granted Critical
Publication of ES2334684T3 publication Critical patent/ES2334684T3/es
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B9/00Safety arrangements
    • G05B9/02Safety arrangements electric
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/31From computer integrated manufacturing till monitoring
    • G05B2219/31359Object oriented model for fault, quality control
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Air Conditioning Control Device (AREA)
  • Testing And Monitoring For Control Systems (AREA)
  • Testing Of Engines (AREA)
  • Air Bags (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

Un método de diagnosis de fallos en sensores para un sistema de tratamiento de aire, que comprende las etapas de: a) crear un modelo de componente y determinar un valor esperado para una medición de un sensor; b) detectar una medición real por parte del sensor; c) comparar el valor esperado con el valor real; y d) determinar la diferencia, caracterizado por que crear el modelo de componente comprende producir una serie de valores de medición, determinar el valor esperado comprende producir una distribución de probabilidades basándose en la serie de valores de medición y determinar una función de densidad de probabilidad para la serie de mediciones, y comparar el valor esperado con el valor real comprende comparar la función de densidad de probabilidad con una función de distribución normal, y por que el método comprende producir una salida de paso bajo y una salida de paso alto con una transformación de escala temporal y/o de frecuencia temporal, comparar la salida de paso bajo y la de paso alto con un valor de umbral predefinido, y determinar un estado de fallo como respuesta al hecho de que una diferencia entre el valor de umbral y la salida de paso bajo y de paso alto es mayor que un valor deseado.

Description

Diagnósticos y pronósticos de fallos de sensores que utilizan un modelo de componente y expansiones ortogonales a escala temporal.
Antecedentes de la invención
Esta invención se refiere a la detección de fallos en sensores de temperatura y de presión para un sistema de calentamiento, ventilación y acondicionamiento de aire.
Típicamente, un sistema de calentamiento, ventilación y acondicionamiento de aire (HVAC -"heating, ventilation and air conditioning") incluye sensores de temperatura y de presión que son cruciales para garantizar el funcionamiento adecuado del sistema. Un controlador se sirve de los datos procedentes de los sensores para ajustar el funcionamiento de los componentes tales como un evaporador, un compresor y un condensador. Si los sensores fallan o son imprecisos, el sistema de HVAC no funcionará como se pretende. Por otra parte, si los sensores no logran proporcionar datos precisos, los componentes del interior del sistema de HVAC pueden funcionar fuera de parámetros deseables, con la posibilidad de provocar un fallo prematuro. Por estas razones, es deseable supervisar la precisión de los sensores e indicar un estado de fallo cuando la precisión de los sensores se encuentra fuera de un intervalo de funcionamiento deseado.
Las técnicas de detección de fallos de sensores de la técnica anterior incluyen el uso de sensores redundantes para medir el mismo parámetro dentro del sistema. Se detecta la diferencia de los datos medidos entre los dos sensores redundantes y se compara esta diferencia con un cierto valor de umbral. Si el valor de umbral es superado por la diferencia real, entonces se detecta un fallo y se inicia una acción correctiva.
Los sistemas que utilizan sensores redundantes son de un coste prohibitivo. Por añadidura, los sensores redundantes hacen que se dupliquen los requisitos de programación y los requisitos del controlador, lo que aumenta la complejidad global del sistema.
Otro sistema conocido de la técnica anterior detecta los fallos en los sensores de presión restando la presión del evaporador de la presión del condensador. Se indica un sensor defectuoso si la diferencia entre el sensor de presión en el evaporador y el sensor de presión en el condensador es cero o negativa. La presión en el condensador debería ser siempre superior a la presión dentro del evaporador, de manera que una desviación de esto indica un fallo. Desventajosamente, esta solución únicamente puede detectar un fallo completo de un sensor. La desviación o sesgo en un sensor según el cual el sensor proporciona una indicación imprecisa, no puede ser detectado por tal sistema.
Otro sistema conocido de la técnica anterior utiliza datos de sensores para determinar una relación lineal que representa un sistema de componentes relacionados estadísticamente. El fallo en un sensor se detecta por una desviación detectada de esa relación lineal. Esta solución fundamentada en datos requiere una gran cantidad de datos relacionados linealmente, así como de un análisis estadístico que es complicado y computacionalmente arduo. Las grandes cantidades de computación pueden conducir a imprecisiones y al incremento de los costes.
De acuerdo con ello, es deseable desarrollar un método y un sistema de detección de fallos en sensores, que trabaje con respuesta inmediata o en tiempo real, utilice características localizadas así como características de tendencia a largo plazo, obtenidas de datos para detectar defectos en sensores, y se adapte a los sistemas actualmente en funcionamiento.
El documento US-B1-6.804.600 divulga un método como el del preámbulo de la reivindicación 1.
Sumario de la invención
La presente invención proporciona un método según se expone en la reivindicación 1, y un sistema según se expone en la reivindicación 15.
Esta invención consiste en un método de diagnóstico de fallos en sensores para sistemas de manejo o tratamiento del aire, que utilizan modelos de componentes y análisis estadístico para evaluar de un modo fiable las condiciones de funcionamiento y los fallos de los sensores.
Un método de este sistema diagnostica fallos en sensores al crear modelos de componentes para predecir y determinar un valor esperado para la medición de un sensor. Los datos reales del sensor son entonces comparados con el modelo de componente y se indica un fallo si una diferencia entre el valor de umbral y las mediciones reales supera un valor predeterminado.
El método de esta invención comprende también las etapas de producir una serie de valores de medición y producir una distribución probabilística de valores de sensor basándose en la serie de mediciones. La función de distribución de probabilidades que se crea basándose en la serie de valores medidos, se compara entonces con una función de distribución normal. La diferencia entre la función de distribución de probabilidades y una función de distribución normal produce puntos de datos que se encuentran fuera de la función de distribución normal. Se utiliza un procedimiento de transformación de escala temporal que produce una salida de pasos bajo y alto. La salida de pasos bajo y alto es entonces comparada con un umbral predeterminado con el fin de diagnosticar fallos en el sensor.
Este sistema proporciona una técnica diagnóstica y predictiva de fallos en sensores eficaz y económica. Esta técnica resulta particularmente útil para tratar los fallos en los sensores que se producen lentamente, tales como la desviación y el sesgo en las lecturas del sensor. Se utilizan modelos de componente del sistema de tratamiento del aire como base para la detección y la estimación del sesgo en el sensor. La combinación de las dos soluciones de detección de fallos, que son la comparación con modelos de componentes y la manipulación del sensor y de los datos utilizando técnicas estadísticas, proporcionan un método fiable para evaluar el estado de un sensor.
De acuerdo con ello, el método y el sistema de esta invención proporcionan un método fiable y económico para determinar el estado de un sensor, que es adaptable a los sistemas actuales.
Breve descripción de los dibujos
Las diversas características y ventajas de esta invención se pondrán de manifiesto para los expertos de la técnica por la siguiente descripción detallada de la realización actualmente preferida. Los dibujos que acompañan la descripción detallada pueden ser brevemente descritos como sigue:
la Figura 1 es una representación esquemática de un sistema de bomba de calor que incluye varios sensores; y
la Figura 2 es un diagrama de bloques que ilustra las etapas del método de esta invención.
Descripción detallada de la realización preferida
Haciendo referencia a la Figura 1, se ilustra en ella un sistema de bomba de calor 10 que incluye un condensador 14, un compresor 12 y un evaporador 16. Una válvula de expansión 18 controla el flujo de refrigerante entre el evaporador 16 y el condensador 14. El compresor 12 impulsa el refrigerante hasta el condensador 14. Existen varios sensores dispuestos dentro del sistema para medir la temperatura y la presión del refrigerante, así como la temperatura del aire que se desplaza al interior del sistema 10 de bomba de calor y que sale de él. El comprensor 12 incluye un sensor de temperatura 24 y un sensor de presión 26 que miden la temperatura y la presión a la salida del compresor 12. En el lado de aspiración o lado de entrada del compresor 12 existe un sensor de presión 30 y un sensor de temperatura 28.
El evaporador 16 incluye sensores de temperatura 32, 34 que miden el aire de entrada y de salida que se desplaza a través del evaporador 16. Por otra parte, unos sensores de temperatura 36, 38 miden la temperatura del refrigerante que se desplaza al interior del evaporador 16 y que sale de éste.
La presión y la temperatura del refrigerante dentro del sistema 10 son descritas por un sistema de ecuaciones que describe la termodinámica del sistema. El funcionamiento termodinámico de un sistema 10 de bomba de calor es bien conocido y las ecuaciones algebraicas que describen tal sistema y que relacionan cada uno de los parámetros del sistema con los otros son bien conocidas.
El compresor 12 es, preferiblemente, un compresor alternativo o de vaivén en el que el proceso de recorrido del vapor refrigerante a través del compresor 12, se supone que es politrópico. Esto significa que las temperaturas de aspiración y de descarga del compresor, conjuntamente con las presiones de aspiración y de descarga, vienen representadas por la relación:
Ecuación 1
1
donde T_{suc}, T_{dis}, P_{suc} y P_{dis} son, respectivamente, la temperatura de entrada al compresor, la temperatura de salida, la presión de entrada y la presión de salida, y n es el índice politrópico. Un valor típico para n es 1,03. El valor exacto de n puede calcularse a partir de los datos del fabricante. Esta relación proporciona un modelo del funcionamiento del compresor 12. Los cuatro sensores 24, 26, 28 y 30 supervisan el estado del compresor 12. El funcionamiento adecuado de los cuatro sensores 24, 26, 28 y 30 es esencial para la fiabilidad del sistema 10 de bomba de calor. Los errores de desviación en los sensores 24, 26, 30 y 28 se detectan de acuerdo con este método, de tal manera que puede emprenderse una acción correctiva. Si los sensores 24, 26, 28 y 30 comienzan a sufrir desviación o deriva, el modelo de componente de compresor proporcionado por la Ecuación 1 ya no es válido.
De acuerdo con ello, la diferencia entre la temperatura de descarga calculada, que se calcula utilizando la Ecuación 1, y el valor real dictado por el sensor 24 indica si se ha producido un fallo. Si los cuatro sensores están libres de defectos, la discrepancia debería ser próxima a 0.
A medida que las lecturas de un sensor comienzan a desviarse alejándose del valor correspondiente al modelo o predicho, se determina un fallo y se da inicio a una señal de instancia al mantenimiento.
La estimación de las desviaciones en los sensores 24, 26, 28 y 30 se lleva a cabo utilizando el siguiente procedimiento de optimización. El procedimiento de optimización viene representado por la ecuación:
Ecuación 3
2
donde \theta_{1}, \theta_{2}, \theta_{3}, \theta_{4} son el valor de sesgo de los sensores. Pueden utilizarse técnicas de optimización no lineal conocidas en la técnica para encontrar la solución para los valores de sesgo.
Si los dos sensores de presión 26, 30 que miden la presión de descarga y la presión de succión están leyendo valores correctos, la ecuación del modelo debería ser válida para cualquier velocidad del compresor. A continuación, se indican las ecuaciones del modelo para el compresor 12, específicamente, de los dos sensores de presión 26, 30 que miden las presiones de descarga y de succión en torno al compresor 12:
Ecuación 4
3
donde V_{suc} es el caudal de flujo volumétrico de succión del compresor, P_{r} es la relación de compresión, P_{suc} es la presión de succión, \dot{\mathit{W}}_{Comp} es la potencia del compresor, y A, B, C, d, e, f y g son constantes proporcionadas por los datos del fabricante.
Los valores medidos de los sensores de presión de succión y de presión de descarga, 26, 30, se utilizan en la Ecuación 4, y arrojan valores para el caudal de flujo volumétrico V_{suc} y \dot{\mathit{W}}_{Comp}. La \dot{\mathit{W}}_{Comp} calculada se compara con el trabajo actual medido por un sensor de potencia o por un medidor de potencia de mano. Si la diferencia entre los dos valores es mayor que un valor de umbral predeterminado, entonces se hace aparecer una indicación de que uno de los sensores de presión es defectuoso. Estas discrepancias se miden a lo largo del tiempo y los valores computados pueden seguirse a lo largo del tiempo. Una tendencia hacia arriba o hacia abajo en los datos para unas condiciones de funcionamiento habituales o estables indicaría una desviación o deriva en los datos.
Para el sistema 10 de bomba de calor, se dispone siempre de sensores capaces de supervisar el comportamiento o rendimiento de los intercambiadores de calor (condensador 14, evaporador 16). Dependiendo del modo de funcionamiento del sistema 10 de bomba de calor, el condensador 14 y el evaporador 16 estarán, en última instancia, expulsando calor al aire circundante o absorbiéndolo de éste. El funcionamiento del condensador 14 y del evaporador 16 puede ser descrito como modelo como un intercambiador de calor de acuerdo con las ecuaciones siguientes:
Ecuación 5
4
Ecuación 6
5
donde Q = la velocidad o ritmo de transferencia de calor, \dot{\mathit{m}}_{1} = caudal de flujo másico de aire, \dot{\mathit{m}}_{2} = caudal de flujo másico de refrigerante, c_{p1} = calor específico del aire seco, T = la temperatura. SHR = la relación de calor sensible ("sensible heat ratio"), y h_{r1}, h_{r2} = entalpías específicas del vapor de refrigerante en la entrada y en la salida del intercambiador de calor.
Las entalpías del refrigerante que fluyen a través de los intercambiadores de calor pueden ser obtenidas a partir de las propiedades del refrigerante, utilizando las mediciones de temperatura y de presión. En un estado en el que se conocen la relación de calor sensible y el caudal de flujo másico del aire, es posible deducir el caudal de flujo de refrigerante por las Ecuaciones 5 y 6, de las que la relación obtenida como resultado es:
Ecuación 7
6
El caudal de flujo másico de refrigerante, por otra parte, puede estimarse utilizando un modelo de componente de compresor. El modelo de componente de compresor puede obtenerse de los datos del fabricante. Un modelo de componente de compresor proporciona una aproximación del caudal de flujo volumétrico teórico que pasa a través del compresor. Esta ecuación viene dada como:
Ecuación 8
7
donde A, B y C son constantes que se estiman a partir de los datos del fabricante. P_{r} es la relación de compresión del compresor, que es la relación que existe entre la presión de descarga P_{dis} y la presión de succión P_{suc}. El caudal de flujo volumétrico puede obtenerse utilizando la densidad del refrigerante, de acuerdo con la ecuación:
Ecuación 9
8
donde \rho es la densidad del refrigerante. La estimación del caudal de flujo de refrigerante de acuerdo con el modelo de compresor debe ser cercana al valor que se ha calculado utilizando las entalpías del refrigerante, según se evalúa en la Ecuación 7. Si hay un fallo en los sensores, existirá una gran discrepancia entre los dos caudales de flujo volumétrico. De acuerdo con ello, si uno de los sensores de temperatura 36, 38 es defectuoso, esta diferencia será mayor que un cierto valor umbral deseado.
Por otra parte, con el fin de determinar si los sensores se están desviando, se determina un valor de estimación de sesgo de sensor mediante el uso del siguiente procedimiento de optimización:
Ecuación 10
9
donde \theta_{5}, \theta_{6}, \theta_{7}, \theta_{8} son valores de sesgo en los cuatro sensores de temperatura 32, 34, 36 y 38 instalados en el intercambiador de calor. Se supone que los sensores de presión no están sesgados, basándose en la suposición de que las lecturas de presión sesgadas han sido corregidas. En consecuencia, únicamente los sesgos en las lecturas de temperatura afectan a las entalpías h_{1r} y h_{2r}.
La solución de la Ecuación 10 proporciona un valor de estimación del sesgo en los sensores. Los valores de sesgo de los sensores son comparados con umbrales predeterminados que, en caso de superarse, indican que se han producido fallos de sesgo en los sensores de presión 26, 30. El seguimiento de los valores de sesgo estimados a lo largo del tiempo hace posible la predicción de fallos en los sensores. Si un fallo de sesgo no es lo suficientemente serio para poner en peligro la integridad del sistema 10, esto es, si el valor de sesgo se encuentra por debajo del valor de umbral para el que se determina un fallo, el valor de sesgo es introducido como entrada en una base de datos estadística con el fin de mejorar y optimizar el sistema en su conjunto.
Un sensor carente de fallos proporciona datos medidos que aportan valores reales conjuntamente con una cierta cantidad de ruido blanco o aleatorio. El ruido blanco es indicativo de las interferencias que acompañan el funcionamiento y la medición de datos por medio de sensores electrónicos.
El hecho de que una serie de datos incluya datos medidos además de ruido blanco, se utiliza en el análisis estadístico de la serie de datos medidos. Para un sensor carente de defectos, una función de densidad de probabilidad (pdf -"probability density function") de valores de datos medidos deberá ser cercana a una distribución normal. Esto es así tras la eliminación de factores deterministas mediante una transformación de tiempo-escala. A medida que los defectos comienzan a desarrollarse, el sensor producirá características o rasgos únicos que aparecen en las series de datos medidos y que pueden utilizarse para detectar el tipo específico de fallo. Por ejemplo, los fallos de desviación o deriva y de sesgo tienen características de baja frecuencia, en tanto que los fallos intermitentes, tales como un cortocircuito ocasionalmente abierto, exhiben propiedades de fallo a frecuencias elevadas. Se utilizan técnicas de transformación de escala para analizar los datos y aislar los rasgos, tanto de baja frecuencia como de alta, que se producen por una serie de datos originada en un sensor defectuoso. Preferiblemente, se produce una transformación de onda pequeña de la serie de datos con el fin de proporcionar una serie de elementos de datos de paso bajo y de paso alto. Estos elementos de datos de paso bajo y de paso alto son revelados dentro de la transformación de onda pequeña y utilizados para detectar los tipos específicos de fallo.
Si se aplica la transformación de onda pequeña a las mediciones de un sensor carente de defectos, la salida de paso alto de la transformación de onda pequeña será próxima a una distribución normal con un valor medio cero. Debido a que la transformación de onda pequeña es un procedimiento de transformación lineal que conserva la función de densidad de probabilidad de la señal original, la distribución deberá ser normal con un valor medio cero. La salida de paso bajo de la transformación de onda pequeña deberá ser una distribución normal con un valor medio que se corresponde con el valor medio de la magnitud física del sensor. En otras palabras, la salida de paso bajo deberá ser, para un sensor carente de defectos y que funciona adecuadamente, una distribución normal con un valor medio que es el mismo que el valor medio de los datos reales medidos. Cualquier varianza de la salida de paso alto con respecto a la distribución normal y con respecto al valor medio que es común con el valor medio de los datos medidos, indica que el sensor es defectuoso.
La desviación de la salida de paso bajo con respecto a una distribución normal y un valor medio común con el valor medio medido, indica un sesgo, congelación o pérdida de precisión dentro del sensor.
Se identifica dentro de la señal un defecto del sensor con fallos intermitentes que generan puntos de datos de frecuencia elevada periódicos. La transformación de escala temporal es sensible a cambios locales en las señales y a los datos de medición procedentes de los sensores. El resultado de esta alta sensibilidad es que se producirán más coeficientes y de una amplitud mayor. La magnitud incrementada de amplitud superior o salida de paso alto que se produce en virtud de la transformación de escala temporal, hace que la función de distribución de probabilidades se desvíe de una distribución normal. Cuando la función de distribución de probabilidades de la transformación de escala temporal se desvía de una distribución normal, este cambio en la salida de paso alto indica que un sensor está experimentando un defecto intermitente. Este defecto intermitente dentro del sensor puede indicar que el sensor tiene algún tipo de fallo eléctrico, o bien que la señal está sufriendo interferencias con un ruido electrónico externo, o sencillamente que se está congelando físicamente.
Haciendo referencia a la Figura 2, el método de diagnosis de fallos en sensores se indica e ilustra esquemáticamente por el diagrama de flujo 50. El método de detección y diagnosis de fallos en sensores para un sistema de HVAC incluye las etapas de recoger datos de un sensor, representadas por el bloque 52. Estos datos son recogidos e introducidos, a continuación, en un controlador de microprocesador. El controlador de microprocesador puede ser tal y como se conoce por un profesional experto en la técnica. Por otra parte, un profesional experto en la técnica comprenderá, con el beneficio de esta descripción, la manera de programar el microcontrolador de que se dispone en el mercado para implementar esta invención y leer datos de sensores.
Cada uno de los sensores es analizado para proporcionar una indicación del estado de ese sensor. Los datos de sensor, representados por el bloque 52, se suministran como entrada al sistema y son almacenados en una base de datos. Los datos almacenados dentro del sistema se someten a continuación al procedimiento de transformación de escala temporal. Preferiblemente, la transformación de escala temporal que se utiliza es una transformación de onda pequeña que descompone los datos reunidos en una salida de paso bajo y de paso alto. Se lleva a cabo entonces un análisis estadístico en la salida, según se indica en el bloque 56. Este análisis estadístico se lleva a cabo en la salida de paso bajo con el fin de determinar si las mediciones estadísticas tales como la media, la varianza y los momentos de orden superior han cambiado con respecto a los datos que se habían reunido con anterioridad. En otras palabras, a medida que se acumula cada punto de datos o conjunto de puntos de datos, éste se compara con algunos de los conjuntos de puntos de datos previos que se han compilado o con todos ellos. A medida que se van acumulando datos adicionales y se utilizan para confeccionar los datos estadísticos del sistema, los datos nuevos son comparados con los datos previamente reunidos, más antiguos, para determinar cualesquiera diferencias. Se lleva a cabo entonces la detección de los cambios en el contenido de paso bajo, con el fin de determinar si ha habido un cambio.
Se lleva también a cabo, además, un análisis estadístico en una salida de paso alto con el fin de determinar si ésta ha cambiado estadísticamente por lo que respecta a la media, la varianza o los momentos de orden superior que constituyeron los datos compilados previamente reunidos.
Una vez realizado el análisis estadístico de la salida de paso bajo y de paso alto que se ha producido a partir de los datos medidos, un dispositivo reconocedor de patrones lleva a cabo un auto-diagnóstico. El dispositivo reconocedor de patrones es un componente interior al controlador que se ha programado de acuerdo con criterios predeterminados para identificar características o rasgos estadísticos que están presentes dentro de las funciones de distribución de salida para cada uno de los sensores. Se identifican los tipos de fallos concretos como un patrón reconocible dentro de los datos reunidos. El tipo de fallo puede incluir fallos eléctricos tales como el sesgo o la desviación o deriva del sensor, o bien fallos mecánicos tales como la congelación o el funcionamiento intermitente, a menudo causado por un cortocircuito o un circuito abierto. De acuerdo con ello, puede identificarse el fallo concreto del sensor según se indica en el bloque 62, a fin de dar instrucciones con respecto a qué acción correctiva se requiere para corregir el fallo.
El análisis estadístico se lleva a cabo utilizando técnicas de expansión según se indica por la referencia 56. Se utiliza, preferiblemente, una técnica de expansión estadística ortogonal. Técnicas conocidas de expansión estadística ortogonal incluyen la expansión ortogonal de Gram-Charlier (GC), que se utiliza preferentemente en este método. Se encuentran también dentro del ámbito contemplado para la invención otras técnicas de expansión según se conoce en la técnica. La idea básica de la expansión GC es que cualquier función de distribución de probabilidades puede ser aproximada por una función de distribución normal estándar o normalizada y términos residuales según se muestra en la ecuación:
Ecuación 11
10
donde \phi(x) es la función de distribución normal estándar o normalizada y r(x) es el residuo de la expansión. H_{i}(x) es la función de Hermite de grado n-ésimo, y c_{i} son los coeficientes de expansión GC, que pueden calcularse por:
11
donde i! significa i factorial.
Si el sensor se encuentra en buen estado de funcionamiento, los datos medidos y obtenidos de ese sensor deberán estar distribuidos de forma normal. Los defectos en el sensor hacen que los datos se desvíen de una función de distribución de probabilidades normal. En la aplicación práctica de la expansión GC, se utiliza un número limitado de términos residuales. Para una expansión de 8 términos, los coeficientes de la expansión GC pueden calcularse como:
Ecuación 12
12
\vskip1.000000\baselineskip
donde m_{k} es el momento de orden k de x en torno a la media. El momento de orden k se define por la relación que se indica a continuación:
Ecuación 13
13
donde \overline{\mathit{x}} es el valor medio de las muestras de datos de sensor x(i), i = 1, 2,..., N.
En general, los momentos estadísticos impares indican información acerca de la densidad de picos de la distribución. Los momentos pares indican características de la extensión de la distribución. Los coeficientes de la expansión GC contienen la información acerca tanto de los momentos impares como de los pares de la función de distribución de probabilidades de los valores de datos, y, como resultado de ello, la expansión GC es más sensible a los cambios que si se usan expansiones ortogonales convencionales. Una propiedad adicional de utilidad de la expansión GC es que es posible formular una estadística a partir de los coeficientes de la expansión GC proporcionando un único índice distribuido (T^{2}) en Chi al cuadrado (\chi^{2}) que puede utilizarse como una indicación de la degradación del estado del sensor. Éste se representa por la ecuación:
Ecuación 14
14
donde los coeficientes individuales c_{i} son \chi^{2} distribuida con n grados de libertad. Si se especifica un nivel de confianza \alpha, entonces se pueden tener los siguientes criterios para la detección de fallos. La detección de los fallos puede determinarse por comparación de los valores obtenidos dentro de las relaciones que se indican a continuación:
Ecuación 15
15
Esto proporciona una relación entre los diferentes tipos de fallos y los coeficientes de GC correspondientes. Ello hace posible el uso de los parámetros de coeficientes de GC para la clasificación del tipo concreto de fallo del sensor.
El sistema de esta invención proporciona un sistema y un método eficaces y económicos para la detección de fallos en sensores. La técnica es útil para hacerse cargo de fallos en sensores que se producen lentamente, tales como la desviación o deriva y el sesgo. Los modelos de componente de los diversos componentes dentro del sistema se utilizan como base para detectar el sesgo del sensor a la hora de estimar la desviación del sensor. Estos modelos de componente se crean utilizando los datos del fabricante que se obtienen para cada uno de los componentes. Se utilizan entonces los modelos para estimar el sesgo que se produce en las lecturas del sensor. Los datos obtenidos de las lecturas del sensor son entonces manipulados por medio de una transformación de escala temporal, y se utilizan expansiones estadísticas de función de distribución de probabilidades para detectar los fallos en el sensor que no son del tipo de desviación o deriva. La combinación de la utilización de modelos de componente y aproximaciones de detección de dos fallos usando la transformación de escala temporal y las expansiones estadísticas de función de distribución de probabilidades, proporciona una forma fiable y adaptable de evaluar la salud de los sensores en ese momento.
La anterior descripción es a título de ejemplo y no sólo una memoria material. La invención se ha descrito de una manera ilustrativa, y debe comprenderse que es la intención que la terminología empleada sea de naturaleza descriptiva en lugar de limitativa. Son posibles muchas modificaciones y variaciones de la presente invención a la luz de las anteriores enseñanzas. Se han descrito las realizaciones preferidas de esta invención; sin embargo, una persona con conocimientos ordinarios en la técnica constatará que ciertas modificaciones se encuentran dentro del ámbito de esta invención. Se comprende que, dentro del ámbito de las reivindicaciones que se acompañan, la invención puede llevarse a la práctica de una forma diferente a la específicamente descrita. Por esta razón, deben estudiarse las siguientes reivindicaciones para determinar el verdadero ámbito y contenido de esta invención.

Claims (20)

1. Un método de diagnosis de fallos en sensores para un sistema de tratamiento de aire, que comprende las etapas de:
a) crear un modelo de componente y determinar un valor esperado para una medición de un sensor;
b) detectar una medición real por parte del sensor;
c) comparar el valor esperado con el valor real; y
d) determinar la diferencia,
caracterizado por que crear el modelo de componente comprende producir una serie de valores de medición, determinar el valor esperado comprende producir una distribución de probabilidades basándose en la serie de valores de medición y determinar una función de densidad de probabilidad para la serie de mediciones, y comparar el valor esperado con el valor real comprende comparar la función de densidad de probabilidad con una función de distribución normal, y por que el método comprende producir una salida de paso bajo y una salida de paso alto con una transformación de escala temporal y/o de frecuencia temporal, comparar la salida de paso bajo y la de paso alto con un valor de umbral predefinido, y determinar un estado de fallo como respuesta al hecho de que una diferencia entre el valor de umbral y la salida de paso bajo y de paso alto es mayor que un valor deseado.
2. El método de acuerdo con la reivindicación 1, que comprende la etapa de llevar a cabo una transformación estadística ortogonal en la función de densidad de probabilidad.
3. El método de acuerdo con la reivindicación 2, en el cual la transformación estadística ortogonal es una transformación de Gram-Charlier.
4. El método de acuerdo con la reivindicación 1, que comprende las etapas de determinar las funciones de densidad de probabilidad de las salidas a partir de una transformación de escala temporal, comparar la función de densidad de probabilidad determinada con una función de densidad de probabilidad predefinida cuando un sensor está sano, y determinar el estado de salud del sensor basándose en los resultados de la comparación.
5. El método de acuerdo con cualquiera de las reivindicaciones precedentes, en el cual dicho modelo de componente representa el funcionamiento de un compresor (12).
6. El método de acuerdo con la reivindicación 5, en el cual dicho modelo de componente que representa el compresor comprende la relación:
16
donde T_{suc}, T_{dis}, P_{suc} y P_{dis} son, respectivamente, la temperatura de entrada al compresor, la temperatura de salida, la presión de entrada y la presión de salida, y n es el índice politrópico.
7. El método de acuerdo con la reivindicación 5, en el cual dicho modelo de componte que representa el compresor (12) comprende la relación:
17
donde V_{suc} es el caudal de flujo volumétrico de succión del compresor, P_{r} es la relación de compresión, P_{suc} es la presión de succión, es la potencia del compresor, y A, B, C, d, e, f y g son constantes proporcionadas por los datos del fabricante.
8. El método de acuerdo con cualquiera de las reivindicaciones precedentes, en el cual dicho sensor es un sensor de presión (26, 30), que detecta la presión en una entrada o en una salida del compresor (12).
9. El método de acuerdo con cualquiera de las reivindicaciones 1 a 4, en el cual dicho sensor es un sensor de temperatura (32, 34, 36, 38) que detecta la temperatura en una entrada o en una salida de un intercambiador de calor (14, 16).
10. El método de acuerdo con cualquiera de las reivindicaciones 1 a 4, en el cual dicho modelo de componente representa el funcionamiento de un intercambiador de calor (14, 16).
11. El método de acuerdo con la reivindicación 10, en el cual dicho funcionamiento de dicho intercambiador de calor (14, 16) viene representado por las ecuaciones:
18
donde Q = la velocidad o ritmo de transferencia de calor, \dot{\mathit{m}}_{1} = caudal de flujo másico de aire, \dot{\mathit{m}}_{2} = caudal de flujo másico de refrigerante, c_{p1} = el calor específico del aire seco, T = la temperatura, SHR = la relación de calor sensible, y h_{r1}, h_{r2} = entalpías específicas del vapor de refrigerante en la entrada y en la salida del evaporador, en todas las unidades compatibles.
12. El método de acuerdo con la reivindicación 10, en el cual se estima un sesgo como la diferencia entre un caudal de flujo determinado por un modelo de intercambiador de calor y un caudal de flujo determinado por un modelo de compresor.
13. El método de acuerdo con cualquiera de las reivindicaciones precedentes, que comprende la etapa de reunir una pluralidad de las diferencias determinadas y hacer un seguimiento de esas diferencias a lo largo del tiempo.
14. El método de acuerdo con la reivindicación 13, que comprende la etapa de predecir un estado futuro de un sensor basándose en la pluralidad de diferencias determinadas.
15. Un sistema para supervisar sensores para un sistema de tratamiento de aire, que comprende: un controlador, que comprende un modelo de componente para determinar un valor esperado de una medición de un sensor, de tal manera que dicho controlador recibe un valor medido de un sensor, que es comparado con dicho valor esperado con el fin de determinar un estado de dicho sensor, de tal modo que dicho controlador una base de datos para producir una distribución de probabilidades basándose en una pluralidad de valores medidos procedentes de dicho sensor, de manera que dicho controlador determina una función de densidad de probabilidad basándose en dicha pluralidad de valores medidos procedentes de dicho sensor, y compara dicha distribución de probabilidades con dicha distribución de probabilidades, de tal forma que dicho controlador lleva a cabo una transformación estadística de escala temporal en dicha función de densidad de probabilidad, y de manera que dicho controlador produce una salida de paso bajo y una salida de paso alto con dicha transformación estadística de escala temporal y determina dicho estado de dicho sensor basándose en una comparación de dicha salida de paso bajo y dicha salida de paso alto con un valor de umbral deseado.
16. El sistema de acuerdo con la reivindicación 15, en el cual dicho controlador determina una diferencia entre dicho valor esperado y dicho valor medido y determina que existe un estado de fallo como respuesta al hecho de que dicha diferencia satisface un criterio deseado.
17. El sistema de acuerdo con la reivindicación 15 ó la reivindicación 16, en el cual dicho modelo de componente comprende una representación del funcionamiento de un compresor.
18. El sistema de acuerdo con la reivindicación 15 ó la reivindicación 16, en el cual dicho modelo de componente comprende una representación del funcionamiento de un intercambiador de calor (14, 16).
19. El sistema de acuerdo con cualquiera de las reivindicaciones 15 a 18, en el cual dicho controlador es independiente de un controlador del sistema de tratamiento de aire.
20. El sistema de acuerdo con cualquiera de las reivindicaciones 15 a 18, en el cual dicho controlador es una parte integral de un controlador para el sistema de tratamiento de aire.
ES05732341T 2004-05-06 2005-04-07 Diagnostico y pronosticos de fallos de sensores que utilizan un modelo de componente y expansiones ortogonales a escala temporal. Active ES2334684T3 (es)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US840566 2001-04-23
US10/840,566 US7200524B2 (en) 2004-05-06 2004-05-06 Sensor fault diagnostics and prognostics using component model and time scale orthogonal expansions

Publications (1)

Publication Number Publication Date
ES2334684T3 true ES2334684T3 (es) 2010-03-15

Family

ID=35240498

Family Applications (1)

Application Number Title Priority Date Filing Date
ES05732341T Active ES2334684T3 (es) 2004-05-06 2005-04-07 Diagnostico y pronosticos de fallos de sensores que utilizan un modelo de componente y expansiones ortogonales a escala temporal.

Country Status (7)

Country Link
US (1) US7200524B2 (es)
EP (1) EP1763754B1 (es)
JP (1) JP2007536490A (es)
AT (1) ATE450830T1 (es)
DE (1) DE602005018056D1 (es)
ES (1) ES2334684T3 (es)
WO (1) WO2005111806A2 (es)

Families Citing this family (60)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI276954B (en) * 2004-12-30 2007-03-21 Ind Tech Res Inst Key performance index and monitoring method thereof
EP2998894B1 (en) 2005-07-11 2021-09-08 Brooks Automation, Inc. Intelligent condition monitoring and fault diagnostic system
US9104650B2 (en) 2005-07-11 2015-08-11 Brooks Automation, Inc. Intelligent condition monitoring and fault diagnostic system for preventative maintenance
DE102006052321A1 (de) * 2005-11-24 2007-06-06 Danfoss A/S Verfahren zum Analysieren einer Kühlanlage und Verfahren zur Regelung einer Kühlanlage
WO2008133715A2 (en) 2006-11-03 2008-11-06 Air Products And Chemicals, Inc. System and method for process monitoring
US7496472B2 (en) * 2007-01-25 2009-02-24 Johnson Controls Technology Company Method and system for assessing performance of control systems
US7802114B2 (en) * 2007-03-16 2010-09-21 Spansion Llc State change sensing
US8024610B2 (en) * 2007-05-24 2011-09-20 Palo Alto Research Center Incorporated Diagnosing intermittent faults
US7800412B2 (en) * 2008-02-12 2010-09-21 Honeywell International Inc. Fault detection and isolation of redundant signals
KR101470631B1 (ko) * 2008-03-12 2014-12-08 엘지전자 주식회사 공기 조화기의 제어방법
US20090252845A1 (en) * 2008-04-03 2009-10-08 Southwick Kenneth J Collider chamber apparatus and method of use
US7827008B2 (en) * 2008-10-21 2010-11-02 General Electric Company System including phase signal saving during anomaly and related method
US20100187320A1 (en) * 2009-01-29 2010-07-29 Southwick Kenneth J Methods and systems for recovering and redistributing heat
US8935037B2 (en) * 2009-08-24 2015-01-13 Robert Bosch Gmbh Good checking for vehicle steering angle sensor
US8738219B2 (en) * 2009-08-24 2014-05-27 Robert Bosch Gmbh Good checking for vehicle longitudinal acceleration sensor
US8754764B2 (en) * 2009-08-24 2014-06-17 Robert Bosch Gmbh Good checking for vehicle pressure sensor
US8401730B2 (en) * 2009-08-24 2013-03-19 Robert Bosch Llc Good checking for vehicle lateral acceleration sensor
US8467929B2 (en) * 2009-08-24 2013-06-18 Robert Bosch Gmbh Good checking for vehicle wheel speed sensors
US8494708B2 (en) * 2009-08-24 2013-07-23 Robert Bosch Gmbh Good checking for vehicle yaw rate sensor
US20110141912A1 (en) * 2009-12-10 2011-06-16 Xerox Corporation Field corrective action based on part life durations
US8386849B2 (en) * 2010-01-29 2013-02-26 Honeywell International Inc. Noisy monitor detection and intermittent fault isolation
US8854226B2 (en) 2010-06-02 2014-10-07 Trane International Inc. System and process for monitoring and scheduling maintenance on mechanical equipment
WO2011153486A1 (en) 2010-06-03 2011-12-08 Cummins Inc. Fresh air flow estimation
DE102010031298A1 (de) * 2010-07-13 2012-01-19 Krones Aktiengesellschaft Verfahren zum Überprüfen der Ausleitung eines Gegenstandes
US20120031985A1 (en) * 2010-08-09 2012-02-09 Terry Lien Do Fault tolerant appliance
WO2012118550A1 (en) * 2011-03-02 2012-09-07 Carrier Corporation Spm fault detection and diagnostics algorithm
JP5812653B2 (ja) * 2011-03-31 2015-11-17 三菱重工業株式会社 熱媒流量推定装置、熱源機、及び熱媒流量推定方法
US8302625B1 (en) * 2011-06-23 2012-11-06 General Electric Company Validation of working fluid parameter indicator sensitivity in system with centrifugal machines
CN203170019U (zh) 2012-02-10 2013-09-04 欧格运动有限责任公司 一种具有可拆卸式弹性垫的飞盘
US9417638B2 (en) * 2012-12-20 2016-08-16 Automotive Research & Testing Center Intelligent thermostatic control method and device for an air conditioner blowing cold and hot air
WO2014105031A2 (en) * 2012-12-28 2014-07-03 Schneider Electric It Corporation Method for air flow fault and cause identification
US9958190B2 (en) 2013-01-24 2018-05-01 Advantek Consulting Engineering, Inc. Optimizing energy efficiency ratio feedback control for direct expansion air-conditioners and heat pumps
US9261542B1 (en) * 2013-01-24 2016-02-16 Advantek Consulting Engineering, Inc. Energy efficiency ratio meter for direct expansion air-conditioners and heat pumps
DE112014000324B4 (de) 2013-08-20 2021-10-07 Jiangsu Tenesun Electrical Appliance Co., Ltd. Online-Selbstdiagnosesystem für einen Wärmepumpen-Wassererhitzer und Steuerverfahren und -gerät dafür
CN104566762B (zh) * 2013-10-12 2018-03-20 珠海格力电器股份有限公司 除湿机的出风口异常监控方法及系统
KR102317340B1 (ko) * 2014-01-14 2021-10-26 삼성전자주식회사 공기 조화기 및 그의 고장 진단 방법
US10316690B2 (en) 2014-07-16 2019-06-11 Siemens Aktiengesellschaft Method for validation of an investigated sensor and corresponding machine
US10296670B2 (en) 2014-09-17 2019-05-21 Infineon Technologies Ag Age estimator for safety monitoring based on local detectors
US10378936B2 (en) 2015-09-04 2019-08-13 International Business Machines Corporation Identification of failed sensors in a system of interconnected devices
CN105423482B (zh) * 2015-10-30 2019-01-22 四川长虹电器股份有限公司 一种温度检测方法及空调
US9892012B2 (en) * 2015-12-30 2018-02-13 International Business Machines Corporation Detecting anomalous sensors
WO2018095785A1 (en) * 2016-11-22 2018-05-31 Danfoss A/S A method for handling fault mitigation in a vapour compression system
WO2018095786A1 (en) 2016-11-22 2018-05-31 Danfoss A/S A method for controlling a vapour compression system during gas bypass valve malfunction
CN109923356B (zh) 2016-11-22 2020-10-13 丹佛斯有限公司 在气体旁通阀故障期间控制蒸气压缩系统的方法
US10337753B2 (en) * 2016-12-23 2019-07-02 Abb Ag Adaptive modeling method and system for MPC-based building energy control
CN107861915B (zh) * 2017-11-09 2021-06-25 东软集团股份有限公司 获取预警阈值的方法、装置以及存储介质
JP2019100687A (ja) * 2017-12-08 2019-06-24 パナソニックIpマネジメント株式会社 空調制御方法及び空調制御装置
US10989427B2 (en) 2017-12-20 2021-04-27 Trane International Inc. HVAC system including smart diagnostic capabilites
CN108829535A (zh) * 2018-06-08 2018-11-16 上海擎创信息技术有限公司 数据处理方法、终端及计算机可读存储介质
CN112074691B (zh) * 2018-06-29 2021-12-14 日立江森自控空调有限公司 空调管理系统、空调管理方法及程序
CN109407640B (zh) * 2018-12-13 2021-03-09 宁波大学 一种基于动态正交成分分析的动态过程监测方法
JP6595139B1 (ja) * 2019-04-05 2019-10-23 日立ジョンソンコントロールズ空調株式会社 空調管理システム、空調管理方法、及びプログラム
CN109974242B (zh) * 2019-04-10 2021-10-08 山东建筑大学 基于热成像的空调系统智能调温方法与系统
CN110425694B (zh) * 2019-08-09 2021-05-25 长江慧控科技(武汉)有限公司 基于phm的高铁智能车站暖通空调能效控制管理方法
CN110513958B (zh) * 2019-08-29 2021-05-11 青岛聚好联科技有限公司 一种确定设备健康状况方法及装置
CN111123885B (zh) * 2019-11-29 2022-05-10 南京航空航天大学 一种基于自适应技术的高超声速飞行器间歇故障诊断方法
US11562137B2 (en) 2020-04-14 2023-01-24 Bank Of America Corporation System to correct model drift for natural language understanding
CN112254270B (zh) * 2020-09-02 2022-05-31 海信(山东)空调有限公司 一种空调及空调故障分类处理方法
US11796201B2 (en) * 2021-04-20 2023-10-24 Lennox Industries Inc. HVAC sensor validation while HVAC system is off
CN113435593B (zh) * 2021-06-29 2022-11-04 广州博通信息技术有限公司 基于传感器时序数据分析的制冷设备结霜预测方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4653280A (en) 1985-09-18 1987-03-31 Hansen John C Diagnostic system for detecting faulty sensors in a refrigeration system
JPH0627598B2 (ja) * 1986-08-13 1994-04-13 三菱重工業株式会社 冷凍装置における圧力センサの故障診断方法
US5267277A (en) * 1989-11-02 1993-11-30 Combustion Engineering, Inc. Indicator system for advanced nuclear plant control complex
GB2255838A (en) * 1991-05-13 1992-11-18 Gen Electric Filtered signal validation.
US5703777A (en) * 1994-10-20 1997-12-30 Anr Pipeline Company Parametric emissions monitoring system having operating condition deviation feedback
SE9600334D0 (sv) * 1996-01-30 1996-01-30 Radi Medical Systems Combined flow, pressure and temperature sensor
US5987399A (en) * 1998-01-14 1999-11-16 Arch Development Corporation Ultrasensitive surveillance of sensors and processes
EP1112532B1 (en) * 1998-08-17 2003-04-02 Aspen Technology, Inc. Sensor validation apparatus and method
JP3411841B2 (ja) * 1998-12-04 2003-06-03 リオン株式会社 故障診断方法および故障診断器
US6917845B2 (en) * 2000-03-10 2005-07-12 Smiths Detection-Pasadena, Inc. Method for monitoring environmental condition using a mathematical model
US6598195B1 (en) * 2000-08-21 2003-07-22 General Electric Company Sensor fault detection, isolation and accommodation
US6804600B1 (en) * 2003-09-05 2004-10-12 Honeywell International, Inc. Sensor error detection and compensation system and method

Also Published As

Publication number Publication date
US20050251364A1 (en) 2005-11-10
EP1763754A4 (en) 2007-08-08
WO2005111806A2 (en) 2005-11-24
DE602005018056D1 (de) 2010-01-14
US7200524B2 (en) 2007-04-03
EP1763754A2 (en) 2007-03-21
WO2005111806A3 (en) 2006-03-16
ATE450830T1 (de) 2009-12-15
EP1763754B1 (en) 2009-12-02
JP2007536490A (ja) 2007-12-13

Similar Documents

Publication Publication Date Title
ES2334684T3 (es) Diagnostico y pronosticos de fallos de sensores que utilizan un modelo de componente y expansiones ortogonales a escala temporal.
Breuker et al. Evaluating the performance of a fault detection and diagnostic system for vapor compression equipment
Wang et al. A system-level fault detection and diagnosis strategy for HVAC systems involving sensor faults
Wang et al. A robust fault detection and diagnosis strategy for centrifugal chillers
US7444251B2 (en) Detecting and diagnosing faults in HVAC equipment
Li et al. Development, evaluation, and demonstration of a virtual refrigerant charge sensor
JP3746729B2 (ja) 機器の劣化を検出する方法
JP2007536073A (ja) エアフィルタ状態の検出および予測技術
CN113654182A (zh) 检测冷媒泄漏的方法和计算机可读存储介质以及空调器
Katipamula et al. Rooftop unit embedded diagnostics: Automated fault detection and diagnostics (AFDD) development, field testing and validation
JP2010025475A (ja) 冷凍サイクル機器に用いられる故障診断装置
Wang et al. Sensor fault detection and diagnosis of air-handling units using a condition-based adaptive statistical method
EP1946066A1 (en) Diagnostic method for proper refrigerant valve operation
JP4301085B2 (ja) 熱源機器の劣化診断システム
US20240142125A1 (en) Air conditioning system, abnormality estimation method for air conditioning system, air conditioner, and abnormality estimation method for air conditioner
CN110046662B (zh) 基于模型与数据融合的冷水机组复合故障诊断方法及系统
Hu Development and demonstration of a method to determine residential air-conditioning system fault prevalence
CN115638507A (zh) 空调系统
JP2004286279A (ja) 空調機の故障診断装置及びその診断方法
JP4231024B2 (ja) 吸収式冷凍機の異常診断方法及びその装置
JP7361935B2 (ja) 空気調和機診断装置
JPH0493567A (ja) 冷凍機の性能診断装置
JP4253648B2 (ja) 吸収冷温水機の性能評価方法と診断システム
CN117782475A (zh) 利用大气压力变化定量监测密闭空间气密性的方法和装置
JP2021188962A (ja) センサドリフト診断装置、センサドリフト診断方法およびセンサドリフト診断プログラム