ES2286184T3 - Procedimiento de evaluacion de las condiciones de operacion de una maquina o de una instalacion. - Google Patents
Procedimiento de evaluacion de las condiciones de operacion de una maquina o de una instalacion. Download PDFInfo
- Publication number
- ES2286184T3 ES2286184T3 ES02019455T ES02019455T ES2286184T3 ES 2286184 T3 ES2286184 T3 ES 2286184T3 ES 02019455 T ES02019455 T ES 02019455T ES 02019455 T ES02019455 T ES 02019455T ES 2286184 T3 ES2286184 T3 ES 2286184T3
- Authority
- ES
- Spain
- Prior art keywords
- extrapolation
- values
- database
- curve
- quality
- 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.)
- Expired - Lifetime
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000011156 evaluation Methods 0.000 title claims abstract description 18
- 238000009434 installation Methods 0.000 title claims abstract description 16
- 238000013213 extrapolation Methods 0.000 claims abstract description 61
- 230000002123 temporal effect Effects 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims description 18
- 238000009826 distribution Methods 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 5
- 230000006978 adaptation Effects 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims description 2
- 238000005303 weighing Methods 0.000 claims 1
- 238000010606 normalization Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000018109 developmental process Effects 0.000 description 3
- 239000006185 dispersion Substances 0.000 description 3
- 230000011218 segmentation Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 239000000446 fuel Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000003723 Smelting Methods 0.000 description 1
- 241000251131 Sphyrna Species 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000005242 forging Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 238000003801 milling Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 238000012731 temporal analysis Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M15/00—Testing of engines
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Complex Calculations (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Procedimiento para la evaluación de las condiciones de operación de una máquina (1) o instalación, para la que se mide en múltiples ocasiones, por lo menos, un parámetro, para constituir una base de datos (6), que consiste en valores (x1, y1)....(xn, yn) del parámetro, calculándose una medida de la calidad (K) de una extrapolación sobre la base de la base de datos (6), caracterizado porque en el procedimiento la medida de la calidad (K) es una función de por lo menos dos variables del grupo V, deltaI, S, C; siendo V una relación del rango de valores de la base de datos (6) para el rango de extrapolación xs, que se determina mediante xs > x1, xn; siendo deltaI la incertidumbre x de la curva de ajuste (21) en dirección x, siendo S la continuidad como una medida de la variación de los valores y de la base de datos (6) y C la constancia temporal de la extrapolación.
Description
Procedimiento de evaluación de las condiciones
de operación de una máquina o de una instalación.
La presente invención hace referencia a un
procedimiento para la evaluación de las condiciones de operación de
una máquina o de una instalación.
Hasta ahora, para la evaluación de las
condiciones de operación de una máquina o de una instalación se
efectuaban extrapolaciones visuales de una evolución, por ejemplo,
una evolución crítica, del valor medido y/o evaluaciones subjetivas
de los efectos sobre los parámetros de la máquina o de la
instalación, para estimar las condiciones de operación de una
máquina o de una instalación y reaccionar correspondientemente,
modificándose, o influyéndose de otro forma, el propio parámetro
correspondientemente a la evaluación, por ejemplo, por la ausencia
de un valor nominal paramétrico. Se ajustan, además, las funciones
de ajuste/regresión a partir de una base de datos y se emprende una
optimización mediante selección interactiva de funciones de curvas,
que representan los máximos coeficientes de correlación. La
evolución de la curvatura de una curva ajustada de este modo no
tiene que coincidir obligatoriamente con la de la base de datos. El
coeficiente de correlación r (valor máximo = 1, valor mínimo = - 1)
sólo puede emplearse condicionalmente como criterio de calidad del
ajuste, ya que este valor no sólo depende de la calidad del ajuste
de la función de la curva utilizada, sino también de la pendiente
de la función de la curva utilizada. Si la pendiente de, por
ejemplo, un ajuste lineal tiende a cero, r sigue también esta
tendencia, independiente de la dispersión de los puntos individuales
de la curva. Por consiguiente, r no puede emplearse como medida de
la calidad para efectuar una extrapolación.
Los instrumentos para la generación para la
evaluación de las condiciones de operación de una máquina o
instalación tienen que satisfacer determinados requisitos mínimos.
Por tanto, es posible una evaluación del resultado de la
extrapolación en lo que se refiere a
\bullet La fiabilidad de la predicción
\bullet Los factores que influyen sobre una
seguridad de la predicción
\bullet La reproducibilidad de la
predicción.
Naturalmente, cada predicción está aquejada de
incertidumbres y su medida de la calidad es extraordinariamente
importante para la evaluación de superaciones pronosticadas de la
alarma/valor límite, así como base de decisión para, por ejemplo,
las acciones automatizadas resultantes de esto. Por otra parte, en
un procedimiento generado cíclicamente para la evaluación de las
condiciones de operación de una máquina o instalación de una
determinada relación, por ejemplo, técnica de medición, puede
considerarse de nuevo como tendencia el desarrollo temporal de una
medida de la calidad y proporciona, por tanto, explicaciones
adicionales sobre el momento de aparición del acontecimiento
pronosticado, por ejemplo, de la superación del valor límite, para
garantizar una seguridad operacional de la máquina o de la
instalación. La estimación de la fiabilidad de una extrapolación
posee importancia central, ya que apenas son posibles las
evaluaciones de los síntomas de daños, superaciones del valor límite
y optimizaciones operacionales, por ejemplo,
reemplazo de piezas desgastadas, sin conocimiento de la fiabilidad de la predicción, o sea, de la medida de la calidad.
reemplazo de piezas desgastadas, sin conocimiento de la fiabilidad de la predicción, o sea, de la medida de la calidad.
Una evaluación de las condiciones de operación
de una máquina o de una instalación, en la que se analiza una
tendencia del desarrollo ulterior de los parámetros, es muy
importante en relación a las soluciones del problema, en caso de
establecimiento de objetivos como supervisión operacional o análisis
del sistema. Estos procedimientos se emplean en todos aquellos
lugares, donde las extrapolaciones en rangos desconocidos de valores
a partir de regularidades conocidas sean de importancia, o sea,
esencialmente en el caso de predicciones de incidencias
(reconocimiento de futuras alarmas/superación de límites de carga) y
prevenciones de daños (reconocimiento de futuros daños).
De manera similar al caso de los valores
medidos, estos son prácticamente inútiles sin indicación de la
tolerancia de medida, se tiene que asignar una variable a una
extrapolación, que proporcione información acerca de la calidad de
la predicción (medida de la calidad). Sin esta información, los
resultados de la extrapolación son apenas evaluables y, por tanto,
las decisiones para acciones subsiguientes (por ejemplo,
interferencias en el proceso) son sumamente inciertas y
posiblemente contraproductivas.
Es, por tanto, objetivo de la presente invención
subsanar este problema.
Un procedimiento, tal y como se representa en el
término genérico de la Reivindicación 1, se conoce gracias al
documento US-A-4 590 562.
El objetivo se resuelve con un procedimiento
conforme a la Reivindicación 1.
En las subreivindicaciones se recogen ulteriores
ordenaciones favorables.
El ajuste óptimo de la función de extrapolación
a la base de datos no se lleva a cabo con un coeficiente máximo de
correlación r, sino con ayuda de una medida de la calidad K, que es
independiente de la pendiente de una curva de ajuste a una base de
datos.
La medida de la calidad K, que da información
sobre la calidad de la predicción de la extrapolación, se calcula a
partir de por lo menos dos subcriterios diferentes. La existencia de
K posibilita al menos dos modos de extrapolación:
En este caso, la amplitud de extrapolación, es
decir, la distancia entre el último valor de la base de datos y un
punto final de la extrapolación, es variable.
Este procedimiento es importante durante la
eliminación de alarmas de fallo o de poca importancia. Aquí, puede
informarse a un operario de una instalación industrial, durante qué
intervalo de tiempo no va a aparecer ninguna alarma, es decir, se
da una determinada seguridad operacional.
En este caso, K es variable. La evolución
temporal de K puede emplearse de nuevo como base para un análisis
de tendencias, para estimar el probable momento de aparición del
acontecimiento para K = 1.
Los análisis de tendencias, para evaluar las
condiciones de operación de una máquina o disposición, pueden
extrapolar completamente, en líneas generales, la evolución de dos
variables correlativas en rangos de valores aún desconocidos. Su
empleo puede aplicarse, tanto a evoluciones temporales
Ejemplo: P_{el} = f(t); P_{el} =
rendimiento del generador
y a asignaciones formalmente
independientes
temporalmente
Ejemplo: P_{el} = f(flujo másico de
combustible).
El empleo del análisis de tendencias supone de
manera implícita, que la regularidad, derivable del pasado, de la
evolución de dos variables correlativas prosigue también en el
futuro y, de este modo, pueden generarse asignaciones de valores
aún desconocidas. Esta condición presupone para los análisis
razonables de tendencias, que éstos sólo puedan generarse con
evoluciones que hayan sido originadas en idénticas
("normalizadas") condiciones de contorno- condiciones físicas
del entorno, combustible, rendimiento, modo de operación, quemador,
etc. La provisión de tales asignaciones normalizadas de valores
puede realizarse, o bien mediante selección de datos con ayuda de
condiciones lógicas (normalización discreta) o con ayuda de un
modelo de simulación del proceso (normalización analítica), que
puede retrocalcular, por ejemplo, todos los valores dependientes
del entorno en condiciones ISO. Muestran:
Figura 1 una base de datos y curva de
extrapolación,
Figura 2 curvas para la determinación de
\DeltaI
Figura 3 una segmentación de una base de
datos,
Figura 4 una evolución temporal de K,
Figura 5 una turbina con piezas,
Figura 6 un corte longitudinal de una turbina en
representación simplificada, y
Figura 7 un álabe de turbina como pieza.
La medida de la calidad K es una medida de la
calidad de una extrapolación. Relaciona varias influencias de fallo
y es, por tanto, una función de al menos dos variables, por ejemplo,
V, \DeltaI. En este caso, son, por ejemplo, una relación V, una
incertidumbre x \DeltaI, una continuidad S y una constancia
temporal C:
Estas variables se relacionan, por ejemplo, de
forma que la provisión de valores de K se encuentre dentro de un
rango normalizado, por ejemplo, del 0% al 100%.
V se calcula, por ejemplo, como sigue:
Si V contribuye considerablemente a una
reducción de la medida de la calidad K, la base de datos de la
extrapolación puede extenderse gradualmente.
Cada ajuste de la curva (condición para la
extrapolación) posee naturalmente también una incertidumbre (zona
de confianza) de los parámetros calculados de la curva. Para
determinar la incertidumbre de la pendiente, se utiliza en la
Figura 2, por ejemplo, una función lineal de regresión 21.
Generalmente pueden emplearse como funciones de extrapolación todas
las funciones, que puedan conducir a estructuras lineales, o
sea:
entre
otras.
La incertidumbre \DeltaI se calcula entonces
como sigue.
La Figura 2 muestra una base de datos 6 ejemplar
con una curva de extrapolación 21 y con otras curvas 24, 27 para el
cálculo de \DeltaI.
La incertidumbre \DeltaI puede representarse
mediante el giro de la curva de ajuste 21 (en representación
lineal) alrededor de un fulcro (x, y). Los valores medios
\overline{x} e \overline{y} se determinaron, por ejemplo, por
medios geométricos o aritméticos. El punto (\overline{x},
\overline{y}) no tiene que encontrarse obligatoriamente sobre la
curva lineal de ajuste 21 (curva de extrapolación). La
incertidumbre de una pendiente de la curva lineal 21 corresponde a
un giro angular, en el que la pendiente m varía \pm\Deltam y,
por tanto, a una incertidumbre en x \DeltaI de la curva de
ajuste, entre otras, en dirección x (Fig. 2). Mediante la variación
de las pendientes resultan dos curvas adicionales 24, 27, que
forman, en cada caso, un punto de corte 30, 33 con una paralela 18
(y_{s}=constante), que corresponde a un valor límite. Para cada
punto de corte 30, 33 existe un valor apropiado de x I_{max} e
I_{min}, siendo I_{max} > I_{min} y \DeltaI = I_{max} -
I_{min}. Esta incertidumbre en x se correlaciona con un ángulo de
corte \alpha, que se determina a partir del ángulo entre la
paralela 18 y la curva en el punto de corte de ambas. La
incertidumbre en x \DeltaI crece al decrecer el ángulo de corte
\alpha. Si en una extrapolación no hay ningún valor límite y/o no
se supera ninguno, resulta válido lo comentado anteriormente para
el punto final en x dado (x_{s}) de la extrapolación y una
paralela al eje x pensada, que recorre el punto (x_{s}, y_{s}).
Como la función de extrapolación existe siempre analíticamente,
puede calcularse con ayuda de la derivada primera del ángulo de
corte \alpha entre una horizontal o función de alarma y la
función de extrapolación.
El comienzo del cálculo de la incertidumbre de
la pendiente \pm\Deltam es, por ejemplo, como sigue: se
determina una pendiente b mediante un intervalo de confianza. Un
procedimiento de este tipo se conoce gracias a Kreyszig, Erwin:
"Métodos Estadísticos y sus Aplicaciones", Editorial
Vandenhoeck y Ruprecht, Göttingen, página 270.
Para la determinación de la pendiente b se
procede como sigue:
a) Se establece una base de datos 6. La base de
datos 6 comprende n valores correlativos de x e y (Fig. 2)
b) cálculo de \overline{x} e \overline{y} de
la base de datos 6, así como de la variable
\Sigmax_{i}y_{i}
c) cálculo de 7
d) cálculo de 8
e) cálculo de 9
f) la pendiente se origina a partir de
10
g) cálculo de 11
Por consiguiente, se origina la ecuación de la
recta de regresión 21
A continuación, se fija un número de
confidencialidad de, por ejemplo, el 95%, es decir \gamma= 0,95,
a partir del cual se calcula una variable F (c):
h) 13
i) Con F(c) = 0,975 y n-2
(n = número de valores medidos) grados de libertad, la distribución
de t (distribución de estudiantes) proporciona un valor c (0,975
corresponde a la integral de la distribución de t hasta la posición
x=c).
\Deltam se origina a partir de
j) 14
Por consiguiente, se origina una incertidumbre
de la pendiente m:
Se calcula, por tanto, la incertidumbre en x
\DeltaI de la pendiente m en el fulcro (Fig. 2) en función de las
características de dispersión de la base de datos actual 6. Por
tanto, \DeltaI = f(\Deltam, \alpha) puede calcularse
por los teoremas de Tales o goniometría. Puede realizarse asimismo
una normalización de \DeltaI.
Una segmentación de una base de datos 6
proporciona, ante todo, información sobre la dirección actual de
curvatura de la base de datos 6 utilizada. La base de datos 6 es una
cantidad parcial verdadera o falsa de una provisión de datos 54, que
comprende todos los puntos de medida 15 disponibles. Una
segmentación significa, que los puntos de medida, o sea, las parejas
de datos (x, y) de la base de datos 6 se dividen en al menos tres
segmentos 45, 48, 52. Para cada segmento 44, 48, 52 se determina una
curva lineal de ajuste 36, 39, 42. La dirección de curvatura de la
base de datos 6 se determina por formación de los valores medios de
los segmentos g_{1}, g_{2} y g_{3} y cálculo de la segunda
derivada numérica.
En este procedimiento puede asumirse que el
cálculo de la curva de ajuste (curva de extrapolación) se lleva a
cabo cíclicamente a partir de una base de datos dinámica y flexible
6, es decir, un rango de datos 3 (Fig. 3) de longitud constante o
variable se desplaza un intervalo \Deltax dado en cada caso
(variable independiente) tras cada ciclo aislado de extrapolación en
la dirección de crecimiento de la variable x. En la práctica, el
avance \Deltax del rango de datos 3, es decir, de la base de datos
6, puede ocurrir como máximo en el ciclo temporal, generándose
nuevos valores medidos 15 (por ejemplo, \Deltax = 5 s). El término
de la continuidad S no ha de entenderse en este contexto como
definición matemática, sino como medida de la variación de la
evolución del punto de la base de datos 6 en relación al último paso
de extrapolación. Para poder hacer una afirmación sobre la
continuidad de la evolución de los valores "y" y, por tanto, de
manera implícita, una afirmación sobre la extrapolabilidad de la
evolución, se subdivide el rango de datos 3 (Fig. 3) en por lo menos
tres segmentos 45, 48, 52. Para cada segmento 45, 48, 52 se calcula,
en cada caso, un valor medio y, en cada caso, una función lineal de
adaptación y_{1}, y_{2}, y_{3} con las pendientes c_{1},
c_{2} y c_{3} (Fig. 3), en un paso intermedio. Si los valores
medios de las tres rectas de regresión 36, 39, 42 se designan por
g_{1}, g_{2} (segmento medio) y g_{3}, puede determinarse la
dirección actual de curvatura de la variación de la pendiente con la
medida numérica de la curvatura p
En el ejemplo de la Figura 3 q es negativa. Esto
significa una curvatura derecha. Para p = 0 los tres valores medios
se encuentran sobre una recta. A partir de una familia de curvas,
ese tipo de curva se ajusta interactivamente al rango de valores de
toda la base actual de datos.
Las funciones de ajuste (funciones de regresión)
de la familia de curvas para la curva de extrapolación de las
correlaciones de datos tienen que satisfacer la condición del
transcurso monótono, porque las funciones
no-monótonas pueden conducir a considerables
incertidumbres de extrapolación. Generalmente pueden emplearse todas
las funciones, que puedan transferirse en estructuras lineales de la
forma y = a_{0} + a_{1}*x, o sea:
entre
otras.
El tipo de curva seleccionado de la familia de
curvas tiene que satisfacer las siguientes condiciones:
La dirección de curvatura de la curva, que se
determina a partir de la extrapolación, tiene que coincidir con la
de p y el cociente Q_{k} de numerador = (dado el caso, valor medio
ponderado del cuadrado de la distancia entre los valores medidos y
la curva de extrapolación) y denominador = cuadrado del valor medio
del rango de valores y de la curva de extrapolación en la zona del
rango de datos (además, por ejemplo, y^{2}_{mitt_{k}} =
[(y_{max_{k}} + y_{max_{k}})/2]^{2}, siendo y_{max}
el valor máximo de Y y correspondiendo y_{min} al valor mínimo de
Y de la curva k-ésima, y sirve para la normalización) tiene que ser
mínimo:
k es, además, un numerador de los
tipos de curva de extrapolación disponibles (familia de curvas);
y_{i}(x_{i}) es el valor medido en la posición x_{i};
f_{k}(x_{i}) es el valor funcional del tipo k-ésimo de
curva de extrapolación en la posición x_{i}; w_{i} es un factor
de ponderación para cada valor medido individual o para todos los
valores medidos de un
segmento.
La evaluación de la continuidad S se efectúa
como sigue: comparación de las tres pendientes de las rectas
C_{1}, C_{2}, C_{3} con las pendientes de la curva de
extrapolación en los respectivos puntos centrales (x_{s1},
x_{s2}, x_{s3}) de los tres segmentos del rango de datos. Aquí
es posible y razonable una compensación diferente de las tres
diferencias de pendiente. Así se pueden evaluar más fuertemente los
valores más actuales en el último segmento, para reconocer más
rápidamente una modificación del transcurso de la curva. Las tres
diferencias de pendiente son una medida de la continuidad de
evolución del transcurso de la curva.
S se calcula como sigue: O_{1} a O_{3} son
las pendientes de la k-ésima curva de ajuste 36, 39, 42 seleccionada
en los, en cada caso, semianchos de segmento; C_{1} a C_{3} las
pendientes de los ajustes de los segmentos lineales (Fig. 3).
Por consiguiente, resulta válido:
con \gammai: coeficientes de
compensación
1...n.
El rango de valores de S es como sigue: 0
(continuidad abs.), es decir las pendientes de la función de
extrapolación en los centros de los segmentos son idénticas a la
recta de regresión de los segmentos individuales 45, 48, 52 a +
\infty, es decir, no existe ninguna coincidencia de las
pendientes.
La varianza temporal de una serie de
extrapolaciones es un indicador importante y adicional de la
seguridad de la predicción, ya que se pueden generar altos valores
absolutos de la varianza, por ejemplo, de la medida de la calidad de
condiciones inestables de extrapolación y, por tanto, dudosa fuerza
expresiva. La consulta de la constancia temporal de una
extrapolación puede realizarse, en falta de una medida firme de la
calidad, con la variable amplitud de extrapolación (x_{s}) o, en
caso de fijación de la amplitud de extrapolación - por ejemplo,
x_{s} puede corresponder al sobrepaso de un valor límite firme -a
través del número confidencial entonces variable. La Figura 4
muestra una evolución temporal ejemplar de la medida de la calidad
K(t) con diversos valores 58. Para esta evolución de los
valores se determina una función de adaptación q a partir de una
familia de curvas. El cálculo de las respectivas variables se lleva
a cabo interactivamente, por ejemplo, con el procedimiento de
Newton. El análisis temporal de la varianza se lleva a cabo con
ayuda de una curva de regresión (por ejemplo, mediante un polinomio
de 1^{er} orden o mayor, siendo el orden función de n), en cada
caso, a los últimos n (por ejemplo, n=10) valores temporales de la
amplitud de extrapolación y/o números confidenciales. De este modo
pueden evaluarse también los valores variables linealmente con el
tiempo de las variables a analizar, con precisión en su
variabilidad. La evaluación de la varianza temporal de la amplitud
de extrapolación se calcula como sigue:
i = número de
interacciones.
También aquí, los coeficientes de compensación
\gamma_{i} pueden evaluar los valores más actuales
diferentes/mayores que los más antiguos. Además,
q^{2}_{mitt_{k}} corresponde a y^{2}_{mitt_{k}}. El rango
de valores de C se encuentra entre C= 0 (ninguna dispersión
temporal de las variables, todos los puntos variables se encuentran
sobre la recta de regresión) y + \infty. (sucesión dispersa de
variables).
La medida de la calidad K relaciona las
variables V, \Delta1_{norm}, S y C, por ejemplo, como sigue:
El rango de valores de K se encuentra entre 0
(extrapolación anulada) y + \infty (extrapolación exacta). Puede
realizarse una normalización de K, por ejemplo, a través de una
relación no lineal:
La instalación es, por ejemplo, un rotor de una
turbina de gas. Mediante la determinación de K del parámetro
frecuencia del rotor se constata que, en un determinado periodo de
tiempo (x_{s}-x_{n}), se verifica una
superación inadmisible del valor límite. Por ello se regulará la
frecuencia hacia abajo.
La Figura 5 muestra esquemáticamente una turbina
de gas 1 en un corte longitudinal.
Una turbina de gas 1 se selecciona como pieza
ejemplar para una máquina, en la que se miden parámetros como la
temperatura, vibración, rendimiento eléctrico u otros parámetros. A
lo largo de un eje 4 se disponen un compresor 7, una cámara de
combustión 10 y una pieza de turbina 13 en serie. La pieza de
turbina 13 presenta un canal de gas recalentado 16. En el canal de
gas recalentado 16 se disponen álabes de la turbina de gas 20. Los
límites de los álabes propulsores y de los álabes móviles se prevén
consecutivamente variables. Los álabes de la turbina de gas 20 se
enfrían, por ejemplo, a través de un enfriamiento combinado de aire
y/o vapor. A tal efecto se extrae, por ejemplo, aire del compresor
del compresor 7 y se alimenta a través de una alimentación de aire
23 a los álabes de la turbina de gas 22. A través de una
alimentación de vapor 26 se alimenta, por ejemplo, también vapor a
los álabes de la turbina de gas 20.
En la Figura 6 se muestra un detalle de una
pieza de turbina 13 en representación simplificada, en un corte
longitudinal. La pieza de turbina 13 presenta un eje 4, que se
extiende a lo largo de un eje de rotación 41. La pieza de turbina
13 presenta adicionalmente una zona de entrada 49, una zona de
paletaje 51, así como una zona de salida 53, consecutivamente a lo
largo del eje de rotación 41. En la zona de paletaje 51 se disponen
álabes móviles rotatorios 20 y álabes directrices fijos 20. Los
álabes móviles 20 se fijan, además, al eje 4, mientras que los
álabes directrices 20 se disponen en un portador de álabes
directrices 47 inmediato al eje 4.
Un canal anular de flujo para un medio
circulante A, por ejemplo, vapor recalentado, está formado por el
eje 4, la zona de paletaje 51, así como el portador de álabes
directrices 47. La zona de entrada 49 que sirve para la
alimentación del medio circulante A queda limitada en dirección
radial por una carcasa de entrada 55 dispuesta aguas arriba del
portador de álabes directrices 47. Una carcasa de salida 57 se
dispone es aguas abajo en el portador de álabes directrices 47 y
limita la zona de salida 53 en dirección radial, es decir,
perpendicular al eje de rotación 41. Durante la operación de la
turbina de gas 1, el medio circulante A circula desde la zona de
entrada 49 hasta la zona de paletaje 51, donde el medio circulante
efectúa trabajo bajo expansión, y abandona posteriormente la
turbina de gas 1 a través de la zona de salida 53. El medio
circulante A se acumula, a continuación, en un condensador aguas
abajo de la carcasa de salida 57, no representado en detalle en la
Figura 5, para una turbina de vapor.
Durante el paso por la zona de paletaje 51 se
expande el medio circulante A y efectúa trabajo en los álabes
móviles 20, desplazándose estos en rotación.
La Figura 7 muestra una vista en perspectiva de
un álabe móvil 20, que se extiende a lo largo de un eje radial 60.
El álabe móvil presenta a lo largo del eje radial 60,
consecutivamente, una zona de fijación 63, una plataforma de álabes
66 adyacente, así como una zona de paletas de álabes 69. En la zona
de fijación 63 se diseña un pie de álabe 72, que sirve para la
fijación del álabe 20 al eje 4 de una turbina de gas 1. El pie de
álabe 72 se configura, por ejemplo, como cabeza de martillo. En los
álabes móviles 20 habituales se utilizan materiales metálicos
macizos en todas las zonas 63, 66, 69. El álabe móvil 20 puede estar
fabricado, en este contexto, mediante un procedimiento de
fundición, de forja, de fresado o combinaciones de estos. Además, la
pieza presenta ya múltiples fallos directamente después de la
elaboración.
Claims (14)
1. Procedimiento para la evaluación de las
condiciones de operación de una máquina (1) o instalación, para la
que se mide en múltiples ocasiones, por lo menos, un parámetro, para
constituir una base de datos (6), que consiste en valores
(x_{1},y_{1})....(x_{n}, y_{n}) del parámetro, calculándose
una medida de la calidad (K) de una extrapolación sobre la base de
la base de datos (6), caracterizado porque en el
procedimiento la medida de la calidad (K) es una función de por lo
menos dos variables del grupo V, \DeltaI, S, C; siendo V una
relación del rango de valores de la base de datos (6) para el rango
de extrapolación x_{s}, que se determina mediante x_{s} >
x_{1}, x_{n}; siendo \DeltaI la incertidumbre x de la curva de
ajuste (21) en dirección x, siendo S la continuidad como una medida
de la variación de los valores y de la base de datos (6) y C la
constancia temporal de la extrapolación.
2. Procedimiento acorde a la Reivindicación 1,
caracterizado porque la evaluación de las condiciones de
operación se emplea, para influir apropiadamente sobre el
parámetro, en función de la medida de la calidad (K).
3. Procedimiento acorde a la Reivindicación 1 ó
2, caracterizado porque mediante la evaluación de las
condiciones de operación se eleva la seguridad operacional de la
máquina (1) o de la instalación, influyéndose apropiadamente sobre
el parámetro, en función de la medida de la calidad (K).
4. Procedimiento acorde a la Reivindicación 1 ó
2, caracterizado porque la evaluación de las condiciones de
operación se utiliza para una optimización del funcionamiento de la
máquina (1) o de la instalación.
5. Procedimiento acorde a la Reivindicación 3,
caracterizado porque para el parámetro se prescribe un valor
límite (18), y porque se determina un intervalo de tiempo, en el que
no ocurre ninguna superación del valor límite (18) del
parámetro.
6. Procedimiento acorde a la Reivindicación 1,
caracterizado porque la selección de las variables se efectúa
de forma que la medida de la calidad (K) sea independiente de la
pendiente de una curva de ajuste a la base de datos
(6).
(6).
7. Procedimiento acorde a la Reivindicación 1,
caracterizado porque la medida de la calidad (K) está
normalizada, particularmente por 1-e^{-K}.
8. Procedimiento acorde a la Reivindicación 7,
caracterizado porque la medida de la calidad (K) se normaliza
a un rango de valores del 0 al 100%.
9. Procedimiento acorde a la Reivindicación 1,
caracterizado porque la medida de la calidad (K) es dada
por:
10. Procedimiento acorde a la Reivindicación 1 ó
9, caracterizado porque la relación (V) del rango de valores
de la base de datos (6) es dada por
(x_{n}-x_{1})/(x_{s}-x_{1}).
11. Procedimiento acorde a la Reivindicación 1,
caracterizado porque la base de datos (6) se subdivide en por
lo menos tres segmentos (45, 48, 52); porque para cada segmento (45,
48, 52) de la base de datos (6) se calcula, en cada caso, un valor
medio g_{1}, g_{2}, g_{3}, y, en cada caso, una función lineal
de adaptación y_{1}, y_{2}, y_{3} (36, 39, 42) con las
pendientes c_{1}, c_{2} y c_{3}; porque se calcula una medida
numérica de la curvatura p
que expresa la dirección actual de
curvatura de la variación de la pendiente; que a partir de una
familia de curvas de tipos de curva de al menos el
grupo:
\newpage
se selecciona iterativamente el
tipo de curva de la función de adaptación y se adapta al rango de
valores de toda la base de datos actual (6), teniendo que cumplir el
tipo de curva seleccionado a partir de la familia de curvas las
siguientes
condiciones:
la dirección de curvatura de la curva, que se
determina a partir de la extrapolación, tiene que coincidir con la
de p y el cociente Q_{k} de numerador (= dado el caso, valor medio
ponderado del cuadrado de la distancia entre los valores medidos y
la curva de extrapolación) y denominador (= cuadrado del valor medio
del rango de valores y de la curva de extrapolación en la zona del
rango de datos) tiene que ser mínimo:
siendo k un numerador de los tipos
de curva disponibles (familia de curvas) para la extrapolación,
particularmente, y^{2}_{mitt_{k}} = [(y_{max_{k}} +
y_{min_{k}})/2]^{2}, siendo y_{i}(x_{i}) el
valor medido en la posición x_{i}, siendo f_{k}(x_{i})
el valor funcional del tipo k-ésimo de curva para la extrapolación
en la posición x_{i}, siendo w_{i} un factor de ponderación para
cada valor medido individual o para todos los valores medidos de un
segmento; de forma que la continuidad (S) se calcula como
sigue:
siendo i = 1...3 una enumeración
para las zonas de segmento, siendo \gammai: coeficientes de
compensación 1...n; O_{1} a O_{3}, las pendientes de la curva
k-ésima seleccionada (36, 39, 42) para la extrapolación en los, en
cada caso, semianchos del segmento; y C_{1} a C_{3}, las
pendientes de los ajustes de los segmentos
lineales.
12. Procedimiento acorde a la Reivindicación 1,
caracterizado porque la incertidumbre x se determina como
sigue:
Selección de una función de extrapolación, que
puede convertirse en estructuras lineales, o sea, se selecciona al
menos del grupo
establecimiento de una base de
datos (6), comprendiendo la base de datos (6) n valores de x e y
correlativos
cálculo de \overline{x} e \overline{y} de la
base de datos (6) así como de la variable \Sigmax_{i}y_{i}
cálculo de 28
cálculo de 29
cálculo de 30
cálculo de una pendiente
31
cálculo de 32
determinación de una ecuación de una recta de
regresión
con un número confidencial \gamma
se calcula una variable F
(c)
con F(c) y
n-2 (n = número de valores medidos) grados de
libertad, la distribución de t (distribución de estudiantes)
proporciona un valor
c,
determinación de \Deltam
resultando una incertidumbre de la
pendiente
m:
determinación de las ecuaciones
lineales (27, 21) con las pendientes b - \Deltam, b +
\Deltam,
determinación de los puntos de corte (I_{min},
constante) e (I_{max}, constante) de la recta con una paralela
(18)
(y = constante), que corresponde a un valor límite (18),
(y = constante), que corresponde a un valor límite (18),
determinación de los correspondientes valores de
x I_{max} y I_{min}, siendo I_{max} > I_{min},
cálculo de \DeltaI = I_{max} -
I_{min}.
13. Procedimiento acorde a la Reivindicación 11,
caracterizado porque el rango de valores de la continuidad
(S) se encuentra en el rango de 0 a + \infty.
14. Procedimiento acorde a la Reivindicación 1,
caracterizado porque la medida de la calidad (K) varía con
el tiempo, y porque una evaluación de la varianza temporal de la
amplitud de extrapolación (C) se calcula como sigue:
siendo i un número de
interacciones, siendo, q^{2}_{mitt_{k}} = [(q_{max_{k}} +
q_{min_{k}})/2]^{2} siendo \Deltai un factor de
ponderación.
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EP02019455A EP1394526B1 (de) | 2002-08-30 | 2002-08-30 | Verfahren zur Bewertung der Betriebsbedingungen einer Maschine oder einer Anlage |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| ES2286184T3 true ES2286184T3 (es) | 2007-12-01 |
Family
ID=31197874
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| ES02019455T Expired - Lifetime ES2286184T3 (es) | 2002-08-30 | 2002-08-30 | Procedimiento de evaluacion de las condiciones de operacion de una maquina o de una instalacion. |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US7194381B2 (es) |
| EP (1) | EP1394526B1 (es) |
| JP (1) | JP4295045B2 (es) |
| CN (1) | CN100420930C (es) |
| DE (1) | DE50210278D1 (es) |
| ES (1) | ES2286184T3 (es) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102004054271A1 (de) * | 2004-11-09 | 2006-05-11 | Fev Motorentechnik Gmbh | Vorhersageverfahren für Betriebszustände einer Verbrennungskraftmaschine |
| EP2012209A1 (de) * | 2007-07-02 | 2009-01-07 | Siemens Aktiengesellschaft | Verfahren zum Ermitteln der Lebensdauer einer Kraftwerkskomponente |
| US20100095232A1 (en) * | 2008-10-15 | 2010-04-15 | Daniel Measurement And Control, Inc. | Calculating and plotting statistical data |
| CN102400902B (zh) * | 2010-09-15 | 2014-12-24 | 中国石油天然气股份有限公司 | 往复式压缩机性能状态可靠性评价方法 |
| CN107609230B (zh) * | 2017-08-24 | 2020-10-16 | 沈阳鼓风机集团股份有限公司 | 一种压缩机调节性能预测方法和系统 |
| RU2681320C1 (ru) * | 2017-12-15 | 2019-03-06 | Общество с ограниченной ответственностью "Западно-Сибирский Центр по сертификации объектов Воздушного транспорта" | Способ оценки остаточного ресурса объектов аттракционной техники в условиях эксплуатации |
| CN111413926A (zh) * | 2020-03-31 | 2020-07-14 | 成都飞机工业(集团)有限责任公司 | 一种持续超限的故障预警方法 |
| RU2756827C1 (ru) * | 2020-06-05 | 2021-10-06 | Общество с ограниченной ответственностью «Западно-Сибирский Центр по сертификации объектов Воздушного транспорта» | Способ продления срока эксплуатации аттракционов |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4590562A (en) * | 1982-12-28 | 1986-05-20 | United Technologies Corporation | Statistically correlated electrostatic engine diagnostics |
| DE4321448C2 (de) * | 1993-06-29 | 1995-08-31 | Kapp Werkzeugmasch | Verfahren zur Ermittlung des Summenteilungsfehlers eines zu bearbeitenden Zahnrades |
| US5847259A (en) * | 1995-12-01 | 1998-12-08 | Ford Motor Company | Computer program, system and method to specify sinusoidal vibration tests for product durability validation |
| CN1287075C (zh) * | 2001-06-18 | 2006-11-29 | 株式会社日立制作所 | 燃气轮机的状态诊断方法及诊断系统 |
-
2002
- 2002-08-30 ES ES02019455T patent/ES2286184T3/es not_active Expired - Lifetime
- 2002-08-30 EP EP02019455A patent/EP1394526B1/de not_active Expired - Lifetime
- 2002-08-30 DE DE50210278T patent/DE50210278D1/de not_active Expired - Lifetime
-
2003
- 2003-08-27 US US10/648,930 patent/US7194381B2/en not_active Expired - Fee Related
- 2003-08-29 JP JP2003306817A patent/JP4295045B2/ja not_active Expired - Fee Related
- 2003-09-01 CN CNB031557724A patent/CN100420930C/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN100420930C (zh) | 2008-09-24 |
| JP2004093567A (ja) | 2004-03-25 |
| US7194381B2 (en) | 2007-03-20 |
| EP1394526A1 (de) | 2004-03-03 |
| JP4295045B2 (ja) | 2009-07-15 |
| DE50210278D1 (de) | 2007-07-19 |
| EP1394526B1 (de) | 2007-06-06 |
| US20040073416A1 (en) | 2004-04-15 |
| CN1492366A (zh) | 2004-04-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20210287122A1 (en) | Prediction device, prediction method, and program | |
| US7827803B1 (en) | Method and apparatus for an aerodynamic stability management system | |
| US8177474B2 (en) | System and method for turbine engine clearance control with rub detection | |
| BR112014001652B1 (pt) | Método com base em computador para determinar uma expectativa de vida residual de um rotor de uma turbina a gás | |
| CN1038333A (zh) | 汽轮机叶片围带间隙监测器 | |
| EP3757356A1 (en) | Control of power generation system by visually monitoring gauge during operation | |
| EP3757696A1 (en) | Control of power generation system by visually monitoring valve during operation | |
| US20210131354A1 (en) | Methods and systems for detection of control sensor override | |
| ES2425875T3 (es) | Método de análisis basado en una técnica de optimización de la función de error | |
| JP2023025354A (ja) | ブレード損傷評価装置、方法及びプログラム | |
| US7899647B2 (en) | Method and device for identifying the operating condition of a turbine | |
| US20120141251A1 (en) | Method and device for predicting the instability of an axial compressor | |
| US7194381B2 (en) | Method for evaluating the operating conditions of a machine or an installation | |
| EP3974646A1 (en) | Method and system for wind speed determination using vibration data | |
| EP3757357A1 (en) | Control of power generation system by visually monitoring component during operation | |
| WO2022239612A1 (ja) | プラント監視方法、プラント監視装置及びプラント監視プログラム | |
| KR20190037643A (ko) | 회전날개의 상태 평가장치 및 그 방법 | |
| KR102775770B1 (ko) | 균열들을 확인하기 위한 터빈 블레이드 상태 모니터링 시스템 | |
| US11105318B2 (en) | Method for calibrating nominal frequencies | |
| EP3990830B1 (en) | Gas turbine swirl detection | |
| JP2001289006A (ja) | タービン車室のキャットバック変形量計測方法 | |
| JP6807950B2 (ja) | 軸流ターボ機械のブレードを輪郭付けするための方法 | |
| JP2018200215A (ja) | 流体機械監視システム、流体機械監視装置、および流体機械監視方法 | |
| US20250103671A1 (en) | Systems and methods of training and using a reduced order model to estimate turbomachine clearances | |
| HK40069811A (en) | Gas turbine swirl detection |