CN108827643A - A kind of high-temperature component of gas turbine fault early warning method for considering to arrange warm temperature field rotation - Google Patents
A kind of high-temperature component of gas turbine fault early warning method for considering to arrange warm temperature field rotation Download PDFInfo
- Publication number
- CN108827643A CN108827643A CN201810646418.8A CN201810646418A CN108827643A CN 108827643 A CN108827643 A CN 108827643A CN 201810646418 A CN201810646418 A CN 201810646418A CN 108827643 A CN108827643 A CN 108827643A
- Authority
- CN
- China
- Prior art keywords
- exhaust temperature
- measuring point
- temperature
- 4pre
- rotation
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000012544 monitoring process Methods 0.000 claims abstract description 14
- 238000012549 training Methods 0.000 claims description 32
- 238000012937 correction Methods 0.000 claims description 26
- 238000005259 measurement Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 2
- 230000002159 abnormal effect Effects 0.000 abstract description 8
- 238000009826 distribution Methods 0.000 abstract description 8
- 238000003745 diagnosis Methods 0.000 abstract description 4
- 230000007423 decrease Effects 0.000 abstract description 2
- 238000002485 combustion reaction Methods 0.000 description 22
- 238000010586 diagram Methods 0.000 description 5
- 238000001514 detection method Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000005856 abnormality Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000009987 spinning Methods 0.000 description 1
- 230000007704 transition Effects 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
- G01M15/14—Testing gas-turbine engines or jet-propulsion engines
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Combustion & Propulsion (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Temperature Or Quantity Of Heat (AREA)
- Control Of Turbines (AREA)
Abstract
一种考虑排温温度场旋转的燃气轮机高温部件故障预警方法,本发明涉及燃气轮机高温部件故障预警方法。本发明为了解决现有方法存在事后诊断的现象,当系统报警时,部件已经发生较严重的损坏以及故障预警的准确性低的问题。与目前已知的同类实时监测方法相比,本发明方法充分利用排温各个测点的数据信息。考虑不同工况对排温的影响,消除了燃气轮机高温气体旋转作用对出口温度场分布的干扰,能够更准确检测出异常演变过程,在故障发生的早期就能及时发现,提高故障预警的准确性。当排温以0.008的斜率下降恶化时,本发明方法相比于传统方法提前了50分钟检测到故障。本发明用于燃气轮机故障预测领域。
A fault early warning method for high temperature components of a gas turbine considering the rotation of exhaust temperature field, the invention relates to a fault early warning method for high temperature components of a gas turbine. The present invention solves the problems of after-the-fact diagnosis in the existing method, when the system alarms, the components have already been seriously damaged and the accuracy of fault early warning is low. Compared with the known real-time monitoring method of the same kind at present, the method of the invention makes full use of the data information of each measuring point of exhaust temperature. Considering the influence of different working conditions on the exhaust temperature, the interference of the high-temperature gas rotation of the gas turbine on the outlet temperature field distribution can be eliminated, and the abnormal evolution process can be detected more accurately, and the fault can be detected in the early stage of the fault, improving the accuracy of fault warning . When the exhaust temperature decreases and deteriorates with a slope of 0.008, the method of the invention detects the fault 50 minutes earlier than the traditional method. The invention is used in the field of gas turbine failure prediction.
Description
技术领域technical field
本发明涉及燃气轮机故障预测领域,考虑排温温度场旋转的燃气轮机高温部件故障预警方法。The invention relates to the field of failure prediction of gas turbines, and relates to a method for early warning of failures of high temperature components of gas turbines in consideration of the rotation of exhaust temperature field.
背景技术Background technique
燃气轮机作为新型的动力设备,具有结构紧凑、运行平稳、安全可靠、可以快速启动并带动负载,具有较高的热效率等优点,在航空、地面和舰船等方面得到了广泛的应用,因此燃气轮机的异常检测对生产实际有着重要意义。在燃气轮机机组运行过程中,对于燃气轮机燃烧系统异常检测方面,50%以上的故障都与燃烧室有关。由于燃烧室燃烧筒等部件长期工作在1600℃的高温区域,工作环境恶劣,设备一旦出现缺陷将可能会对下游的喷嘴和动叶部件安全构成威胁。因此,有必要对燃烧室的工作状况进行监控。As a new type of power equipment, gas turbine has the advantages of compact structure, stable operation, safety and reliability, quick start and load driving, and high thermal efficiency. It has been widely used in aviation, ground and ships. Therefore, the gas turbine Anomaly detection is of great significance to production practice. During the operation of the gas turbine unit, more than 50% of the faults in the abnormal detection of the gas turbine combustion system are related to the combustion chamber. Since the combustor, combustor and other components work in the high temperature area of 1600 ° C for a long time, the working environment is harsh, and once the equipment is defective, it may pose a threat to the safety of downstream nozzles and moving blade components. Therefore, it is necessary to monitor the working condition of the combustion chamber.
燃烧系统一旦出现故障,会使燃烧室出口温度产生差异。因此我们可以通过检测燃烧室出口温度来监测燃烧系统的运行状况。但是,常规的温度测量元件无法在如此高温的区域长期工作,因此,在机组透平排气通道中周向均匀布置了若干个排气测温热电偶,热电偶所测的温度就是燃气轮机的排温。用户通过监视燃气轮机的排气温度来间接监视燃烧室内的工作状况。在实际运行时,当燃烧筒出现异常的时候,排温的结果也会出现异常,所以就可以通过排温的情况来判断燃烧筒的工作情况是否出现异常。热电偶分布情况如图1所示。Once the combustion system fails, it will cause a difference in the outlet temperature of the combustion chamber. Therefore, we can monitor the operation status of the combustion system by detecting the outlet temperature of the combustion chamber. However, conventional temperature measuring elements cannot work in such a high temperature area for a long time. Therefore, several exhaust temperature measuring thermocouples are evenly arranged circumferentially in the turbine exhaust passage of the unit. The temperature measured by the thermocouples is the exhaust temperature of the gas turbine. temperature. The user indirectly monitors the working conditions in the combustor by monitoring the exhaust temperature of the gas turbine. In actual operation, when the combustion tube is abnormal, the result of exhaust temperature will also be abnormal, so it can be judged whether the working condition of the combustion tube is abnormal according to the exhaust temperature. The distribution of thermocouples is shown in Figure 1.
GE公司开发的MARK VI燃烧监测系统定义S为排气温度的允许排温分散度,认为S是燃气轮机出口的平均排气温度T4 *、压气机出口温度T4 *的函数,具体函数是个经验公式:The MARK VI combustion monitoring system developed by GE defines S as the allowable dispersion of exhaust gas temperature, and considers that S is a function of the average exhaust temperature T 4 * at the outlet of the gas turbine and the temperature T 4 * at the outlet of the compressor. The specific function is an experience formula:
在该公式里,温度均是以℉为计量单位的。公式右端的100带有括号,表示变工况条件下才加入该项。In this formula, temperature is measured in °F. The 100 at the right end of the formula has brackets, which means that this item is only added under variable working conditions.
此外,MARK VI燃烧监测系统还定义:S1为排气温度热电偶的最高读数与最低读数之间的差;S2为排气温度热电偶的最高读数与第2个低读数之间的差;S3为排气温度热电偶的最高读数与第3个低读数之间的差。In addition, the MARK VI combustion monitoring system also defines: S1 is the difference between the highest reading and the lowest reading of the exhaust temperature thermocouple; S2 is the difference between the highest reading and the second lowest reading of the exhaust temperature thermocouple; S3 is the difference between the highest reading and the 3rd lowest reading of the discharge temperature thermocouple.
基于上述的公式和定义,MARK VI燃烧监测保护系统的判别原理见图2。图2中,K1,K2,K3是三个依据经验定义的参数。典型情况下:K1=1.0;K2=5.0;K3=0.8;Based on the above formulas and definitions, the discrimination principle of MARK VI combustion monitoring and protection system is shown in Figure 2. In Fig. 2, K 1 , K 2 , and K 3 are three parameters defined based on experience. Typically: K 1 =1.0; K 2 =5.0; K 3 =0.8;
在实际应用中发现,该种方法检测不出异常演变的过程,无法对燃烧状态变化趋势做出判断,存在严重的“事后”诊断现象,即当检测系统发出报警时燃烧系统已经损坏较严重。In practical application, it is found that this method cannot detect the abnormal evolution process, and cannot make judgments on the changing trend of the combustion state. There is a serious "post-event" diagnosis phenomenon, that is, the combustion system has been seriously damaged when the detection system sends out an alarm.
同时,现有的技术都是假设机组的透平流道是均匀的,但是在实际运行过程中,不同工况会使高温燃气产生旋转的现象,如图3所示,从而使燃气轮机出口温度场分布的产生差异。旋转产生的差异对于燃烧系统发生故障时燃烧室出口温度产生的差异是一种干扰因素,在早期预警中,旋转造成的影响可能比燃烧筒异常产生的影响还要大,如图4所示,导致无法检测燃烧室出现异常情况。除此之外,不同工况还会导致不同测点温度的缩放。因此,需要考虑不同工况对排温的影响,消除旋转对燃气轮机出口温度分布的影响。At the same time, the existing technologies all assume that the turbine flow path of the unit is uniform, but in the actual operation process, different working conditions will cause the high-temperature gas to rotate, as shown in Figure 3, so that the gas turbine outlet temperature field distribution of difference. The difference caused by the rotation is an interference factor for the difference in the outlet temperature of the combustion chamber when the combustion system fails. In the early warning, the influence caused by the rotation may be greater than the abnormality of the combustion cylinder, as shown in Figure 4. As a result, abnormalities in the combustion chamber cannot be detected. In addition, different working conditions will also lead to the scaling of the temperature of different measuring points. Therefore, it is necessary to consider the influence of different working conditions on the exhaust temperature and eliminate the influence of rotation on the gas turbine outlet temperature distribution.
发明内容Contents of the invention
本发明的目的是为了解决现有方法存在事后诊断的现象,当系统报警时,部件已经发生较严重的损坏以及故障预警的准确性低的缺点,而提出一种考虑排温温度场旋转的燃气轮机高温部件故障预警方法。The purpose of the present invention is to solve the problem of after-the-fact diagnosis in the existing method. When the system alarms, the components have been seriously damaged and the accuracy of fault early warning is low, and a gas turbine considering the rotation of the exhaust temperature field is proposed. High temperature component failure early warning method.
一种考虑排温温度场旋转的燃气轮机高温部件故障预警方法包括以下步骤:A fault early warning method for high temperature components of a gas turbine considering the rotation of the exhaust temperature field includes the following steps:
步骤一:燃气轮机在透平出口周向均匀分布n个热电偶测量燃气轮机排气温度,燃气轮机在正常运行时间[t1,tm],得到了m个时刻n个测点的排温数据序列,分别构成排温数据序列和T4,1,T4,2,…,T4,n,将t1时刻的排温数据序列作为训练集的参考序列;tm时刻的排温数据序列作为训练集的待处理对比序列;T4为排气温度;Step 1: The exhaust temperature of the gas turbine is measured by n thermocouples evenly distributed in the circumferential direction of the turbine outlet. During the normal operation time of the gas turbine [t 1 , t m ], the exhaust temperature data sequence of n measuring points at m times is obtained. Respectively constitute the exhaust temperature data series and T 4,1 , T 4,2 ,..., T 4,n , the exhaust temperature data sequence at time t 1 As the reference sequence of the training set; the exhaust temperature data sequence at time t m As the comparison sequence to be processed in the training set; T 4 is the exhaust temperature;
所述为t2时刻的排温数据序列,T4,1为第一个测点的排温数据序列,T4,2为第一个测点的排温数据序列,T4,n为第n个测点的排温数据序列;said is the exhaust temperature data sequence at time t 2 , T 4,1 is the exhaust temperature data sequence of the first measuring point, T 4,2 is the exhaust temperature data sequence of the first measuring point, T 4,n is the nth Exhaust temperature data series of measuring points;
步骤二:定义训练集的待处理对比序列与训练集的参考序列的偏差为e,采用欧氏距离判别法确定e的最小值,得到e取最小值时的m个时刻的旋转修正后的排温数据序列和n个测点的旋转修正后的排温数据序列T′4,1,T′4,2,…,T′4,n;Step 2: Define the comparison sequence to be processed in the training set Reference sequence with training set The deviation of e is e, and the minimum value of e is determined by the Euclidean distance discriminant method, and the exhaust temperature data sequence after rotation correction at m moments when e takes the minimum value is obtained and the rotation-corrected exhaust temperature data sequence T′ 4,1 , T′ 4,2 ,…, T′ 4,n of n measuring points;
步骤三:确定n个测点的旋转修正后的排温数据序列T′4,1,T′4,2,…,T′4,n与n个测点的平均温度之间的关系;Step 3: Determine the relationship between the rotation-corrected exhaust temperature data series T′ 4,1 , T′ 4,2 , ..., T′ 4,n of n measuring points and the average temperature of n measuring points;
步骤四:根据步骤三确定的m个时刻的旋转修正后的排温数据序列与n个测点的平均温度之间的关系,得到n个测点的温度预测值T4pre,1,T4pre,2,…,T4pre,i,…,T4pre,n;Step 4: Exhaust temperature data sequence corrected by rotation at m moments determined in step 3 and the relationship between the average temperature of n measuring points to obtain the predicted temperature values T 4pre,1 ,T 4pre,2 ,...,T 4pre,i ,...,T 4pre,n of n measuring points;
步骤五:根据步骤四得到的n个测点的温度预测值T4pre,1,T4pre,2,…,T4pre,i,…,T4pre,n,得到每一个测点的误差带ΔT1,ΔT2,…,ΔTi,…,ΔTn;其中ΔTi为第i个测点的误差带;Step 5: According to the temperature prediction values T 4pre,1 ,T 4pre,2 ,…,T 4pre,i ,…,T 4pre,n of n measuring points obtained in Step 4, get the error band ΔT 1 of each measuring point ,ΔT 2 ,…,ΔT i ,…,ΔT n ; where ΔT i is the error band of the i-th measuring point;
步骤六:监测时,重新执行步骤一至步骤五,若得到的第i个测点的误差带在ΔTi的范围内,则说明燃气轮机机组无故障,否则燃气轮机机组故障。Step 6: When monitoring, re-execute steps 1 to 5. If the error band of the i-th measuring point is within the range of ΔT i , it means that the gas turbine unit is not faulty, otherwise the gas turbine unit is faulty.
在监测时,将每个时刻测得的排温数据代入,得到的测点1的预测值T4pre,1。监测每个时刻机组无故障时测点1的预测值T4pre,1与测点1的实际值(修正值)T4',1t之差为ΔT1,则说明机组无故障,若超出该范围,则说明机组发生故障。During monitoring, the exhaust temperature data measured at each moment is substituted into the predicted value T 4pre,1 of measuring point 1 . The difference between the predicted value T 4pre,1 of measuring point 1 and the actual value (corrected value) T 4 ' ,1t of measuring point 1 when the unit is not faulty at each moment of monitoring is ΔT 1 , which means that the unit is not faulty. If it exceeds this range , it means that the unit is malfunctioning.
分别得到若机组无故障时的测点2至n的预测值与实际值之差ΔT2,ΔT3,…,ΔTn,若偏差在ΔTi范围内则说明机组无故障,若超出该范围,则说明机组发生故障。The difference between the predicted value and the actual value ΔT 2 , ΔT 3 ,...,ΔT n of the measuring points 2 to n when the unit has no faults is respectively obtained. If the deviation is within the range of ΔT i , it means that the unit has no faults. If it exceeds this range, It means that the unit is malfunctioning.
所述高温部件是指燃气轮机燃烧室中的火焰筒、过渡段及首级透平喷嘴、动叶片等部件。The high-temperature components refer to the flame cylinder, the transition section, the first stage turbine nozzle, the moving blade and other components in the combustion chamber of the gas turbine.
本发明的有益效果为:The beneficial effects of the present invention are:
目前已有的方法存在事后诊断的现象,当系统报警时,部件已经发生较严重的损坏。为实现燃气轮机高温部件故障的早期预警,及早发现故障,本发明提供了一种考虑排温温度场旋转的燃气轮机高温部件故障早期预警方法。The current existing methods have the phenomenon of post-event diagnosis. When the system alarms, the components have been seriously damaged. In order to realize the early warning of the failure of the high temperature component of the gas turbine and find the failure early, the present invention provides an early warning method for the failure of the high temperature component of the gas turbine considering the rotation of the exhaust temperature field.
与目前已知的同类实时监测方法相比,本发明方法可以充分利用排温各个测点的数据信息。同时,考虑不同工况对排温的影响,消除了燃气轮机高温气体旋转作用对出口温度场分布的干扰,能够更准确检测出异常演变过程,在故障发生的早期就能及时发现,提高故障预警的准确性。因此,可以降低因为燃气轮机产生故障不能及时发现所造成的损失和可能。当排温以0.008的斜率下降恶化时,本发明方法相比于传统方法提前了50分钟检测到故障。Compared with the known real-time monitoring methods of the same kind at present, the method of the invention can make full use of the data information of each measuring point of exhaust temperature. At the same time, considering the influence of different working conditions on the exhaust temperature, the interference of the high-temperature gas rotation of the gas turbine on the outlet temperature field distribution can be eliminated, and the abnormal evolution process can be detected more accurately, which can be detected in the early stage of the fault and improve the effectiveness of fault early warning. accuracy. Therefore, it can reduce the loss and possibility caused by the gas turbine failing to be discovered in time. When the exhaust temperature decreases and deteriorates with a slope of 0.008, the method of the invention detects the fault 50 minutes earlier than the traditional method.
附图说明Description of drawings
图1为热电偶分布情况图;Figure 1 is a diagram of the distribution of thermocouples;
图2为燃烧监测的判别原理图;Fig. 2 is the discrimination schematic diagram of combustion monitoring;
图3为不同工况旋转的现象图;图中C0为气流在静叶进口的绝对速度(米/秒),C1为气流在动叶进口的绝对速度(米/秒),C2为气流在动叶出口的绝对速度(米/秒),u为动叶牵连速度(米/秒),α1为C1绝对速度方向角,β1为相对速度方向角,α2为C2绝对速度方向角,β2为相对速度方向角;Figure 3 is the phenomenon diagram of rotation under different working conditions; C 0 in the figure is the absolute velocity of the airflow at the inlet of the stator blade (m/s), C 1 is the absolute velocity of the airflow at the inlet of the moving blade (m/s), and C 2 is The absolute velocity of the airflow at the outlet of the rotor blade (m/s), u is the implicated velocity of the rotor blade (m/s), α 1 is the absolute velocity direction angle of C 1 , β 1 is the relative velocity direction angle, α 2 is the absolute velocity direction angle of C 2 Velocity direction angle, β2 is the relative velocity direction angle;
图4为旋转前后排温分布图;Figure 4 is a diagram of exhaust temperature distribution before and after rotation;
图5为训练集数据图,图中横坐标为时间,纵坐标为功率;Fig. 5 is a training set data diagram, in which the abscissa is time, and the ordinate is power;
图6为测试集正常结果图;Figure 6 is a normal result graph of the test set;
图7为测试集异常结果图。Figure 7 is a graph of the abnormal results of the test set.
具体实施方式Detailed ways
具体实施方式一:一种考虑排温温度场旋转的燃气轮机高温部件故障预警方法包括以下步骤:Embodiment 1: A method for early warning of failure of high temperature components of a gas turbine considering the rotation of the exhaust temperature field includes the following steps:
步骤一:燃气轮机在透平出口周向均匀分布n个热电偶测量燃气轮机排气温度,燃气轮机在正常运行时间[t1,tm],得到了m个时刻n个测点的排温数据序列,分别构成排温数据序列和T4,1,T4,2,…,T4,n,将t1时刻的排温数据序列作为训练集的参考序列;tm时刻的排温数据序列作为训练集的待处理对比序列;Step 1: The exhaust temperature of the gas turbine is measured by n thermocouples evenly distributed in the circumferential direction of the turbine outlet. During the normal operation time of the gas turbine [t 1 , t m ], the exhaust temperature data sequence of n measuring points at m times is obtained. Respectively constitute the exhaust temperature data series and T 4,1 , T 4,2 ,..., T 4,n , the exhaust temperature data sequence at time t 1 As the reference sequence of the training set; the exhaust temperature data sequence at time t m As the comparison sequence to be processed in the training set;
为m个时刻的排温数据序列,T4,1,T4,2,…,T4,n为n个测点的排温数据序列; is the exhaust temperature data sequence of m moments, T 4,1 , T 4,2 ,..., T 4,n is the exhaust temperature data sequence of n measuring points;
所述为t2时刻的排温数据序列,T4,1为第一个测点的排温数据序列,T4,2为第一个测点的排温数据序列,T4,n为第n个测点的排温数据序列;said is the exhaust temperature data sequence at time t 2 , T 4,1 is the exhaust temperature data sequence of the first measuring point, T 4,2 is the exhaust temperature data sequence of the first measuring point, T 4,n is the nth Exhaust temperature data series of measuring points;
步骤二:定义训练集的待处理对比序列与训练集的参考序列的偏差为e,采用欧氏距离判别法确定e的最小值,得到e取最小值时的m个时刻的旋转修正后的排温数据序列和n个测点的旋转修正后的排温数据序列T′4,1,T′4,2,…,T′4,n;Step 2: Define the comparison sequence to be processed in the training set Reference sequence with training set The deviation of e is e, and the minimum value of e is determined by the Euclidean distance discriminant method, and the exhaust temperature data sequence after rotation correction at m moments when e takes the minimum value is obtained and the rotation-corrected exhaust temperature data sequence T′ 4,1 , T′ 4,2 ,…, T′ 4,n of n measuring points;
步骤三:确定n个测点的旋转修正后的排温数据序列T′4,1,T′4,2,…,T′4,n与n个测点的平均温度之间的关系;Step 3: Determine the relationship between the rotation-corrected exhaust temperature data series T′ 4,1 , T′ 4,2 , ..., T′ 4,n of n measuring points and the average temperature of n measuring points;
步骤四:根据步骤三确定的m个时刻的旋转修正后的排温数据序列与n个测点的平均温度之间的关系,得到n个测点的温度预测值T4pre,1,T4pre,2,…,T4pre,i,…,T4pre,n;Step 4: Exhaust temperature data sequence corrected by rotation at m moments determined in step 3 and the relationship between the average temperature of n measuring points to obtain the predicted temperature values T 4pre,1 ,T 4pre,2 ,...,T 4pre,i ,...,T 4pre,n of n measuring points;
步骤五:根据步骤四得到的n个测点的温度预测值T4pre,1,T4pre,2,…,T4pre,i,…,T4pre,n,得到每一个测点的误差带ΔT1,ΔT2,…,ΔTi,…,ΔTn;其中ΔTi为第i个测点的误差带;Step 5: According to the temperature prediction values T 4pre,1 ,T 4pre,2 ,…,T 4pre,i ,…,T 4pre,n of n measuring points obtained in Step 4, get the error band ΔT 1 of each measuring point ,ΔT 2 ,…,ΔT i ,…,ΔT n ; where ΔT i is the error band of the i-th measuring point;
步骤六:监测时,重新执行步骤一至步骤五,若得到的第i个测点的误差带在ΔTi的范围内,则说明燃气轮机机组无故障,否则燃气轮机机组故障。Step 6: When monitoring, re-execute steps 1 to 5. If the error band of the i-th measuring point is within the range of ΔT i , it means that the gas turbine unit is not faulty, otherwise the gas turbine unit is faulty.
在监测时,将每个时刻测得的排温数据代入,得到的测点1的预测值T4pre,1。监测每个时刻机组无故障时测点1的预测值T4pre,1与测点1的实际值(修正值)T′4,1t之差为ΔT1,则说明机组无故障,若超出该范围,则说明机组发生故障。During monitoring, the exhaust temperature data measured at each moment is substituted into the predicted value T 4pre,1 of measuring point 1 . The difference between the predicted value T 4pre,1 of measuring point 1 and the actual value (corrected value) T′ 4,1t of measuring point 1 when there is no fault in the monitoring unit at each moment is ΔT 1 , which means that the unit has no fault. If it exceeds this range , it means that the unit is malfunctioning.
分别得到若机组无故障时的测点2至n的预测值与实际值之差ΔT2,ΔT3,…,ΔTn,若偏差在ΔTi范围内则说明机组无故障,若超出该范围,则说明机组发生故障。The difference between the predicted value and the actual value ΔT 2 , ΔT 3 ,...,ΔT n of the measuring points 2 to n when the unit has no faults is respectively obtained. If the deviation is within the range of ΔT i , it means that the unit has no faults. If it exceeds this range, It means that the unit is malfunctioning.
训练集数据如图5所示,测试集结果如图6和图7所示。The training set data is shown in Figure 5, and the test set results are shown in Figures 6 and 7.
具体实施方式二:本实施方式与具体实施方式一不同的是:所述步骤一中得到m个时刻的排温数据序列以及n个测点的排温数据序列T4,1,T4,2,…,T4,n具体为:Specific embodiment 2: The difference between this embodiment and specific embodiment 1 is that in the step 1, the temperature exhaust data series at m times are obtained And the exhaust temperature data series T 4,1 , T 4,2 ,..., T 4,n of n measuring points are specifically:
m个时刻的排温数据序列为:Exhaust temperature data series at m time points for:
……...
其中分别为t1时刻第1个测点到第n个测点的温度值,in are the temperature values from the first measuring point to the nth measuring point at time t1 respectively,
分别为tm时刻第1个测点到第n个测点的温度值; are the temperature values from the first measuring point to the nth measuring point at time t m respectively;
n个测点的排温数据序列T4,1,T4,2,…,T4,n具体为:The exhaust temperature data series T 4,1 , T 4,2 ,..., T 4,n of n measuring points are specifically:
……...
其中分别为第1个测点t1到tm时刻的温度值,分别为第2个测点t1到tm时刻的温度值,分别为第n个测点t1到tm时刻的温度值。in are the temperature values of the first measuring point from t1 to tm , respectively, are the temperature values at the second measuring point from t 1 to t m respectively, are the temperature values at the nth measuring point from t1 to tm , respectively.
其它步骤及参数与具体实施方式一相同。Other steps and parameters are the same as those in Embodiment 1.
具体实施方式三:本实施方式与具体实施方式一或二不同的是:所述步骤二中得到e取最小值时的m个时刻的旋转修正后的排温数据序列和n个测点的旋转修正后的排温数据序列T′4,1,T′4,2,…,T′4,n的具体过程为:Specific embodiment 3: The difference between this embodiment and specific embodiment 1 or 2 is that in the step 2, the rotation-corrected exhaust temperature data sequence at m moments when e takes the minimum value is obtained The specific process of the exhaust temperature data sequence T′ 4,1 , T′ 4,2 , ..., T′ 4,n after rotation correction of n measuring points is:
当不考虑旋转影响时,训练集的待处理对比序列与训练集的参考序列的偏差e0为:When the effect of rotation is not considered, the sequence of comparisons to be processed in the training set Reference sequence with training set The deviation e 0 of is:
当旋转一个测点时,训练集的待处理对比序列与训练集的参考序列的偏差e1为:When rotating a survey point, the sequence of comparisons to be processed in the training set Reference sequence with training set The deviation e1 of is:
当旋转 when spinning
当旋转两个测点时,训练集的待处理对比序列与训练集的参考序列的偏差e2为:The sequence of comparisons to be processed for the training set when rotating two survey points Reference sequence with training set The deviation e2 of is:
当旋转N个测点时,训练集的待处理对比序列与训练集的参考序列的偏差eN为:When rotating N measurement points, the comparison sequence to be processed in the training set Reference sequence with training set The deviation e N of is:
其中e=[e0,e1,…,eN];where e=[e 0 ,e 1 ,...,e N ];
取e中的最小值为:Take the smallest value in e:
emin=argmin(e1,e2,…,eN)e min =argmin(e 1 ,e 2 ,…,e N )
argmin代表取最小值;argmin represents the minimum value;
当旋转了N个测点时,e取最小值,对不同时刻排温序列都进行如上操作,可以m个时刻的旋转修正后的排温数据序列:When N measuring points are rotated, e takes the minimum value, and the above operation is carried out for the temperature exhaust sequence at different times, and the corrected exhaust temperature data sequence can be rotated at m times:
……...
其中分别为t1时刻第1个测点到第n个测点旋转修正后的温度值,分别为tm时刻第1个测点到第n个测点旋转修正后的温度值;in Respectively, the temperature values after rotation correction from the first measuring point to the nth measuring point at time t1 , Respectively, the temperature values after rotation correction from the first measuring point to the nth measuring point at time t m ;
如t1时刻,训练集的待处理(对比)序列旋转了N个测点时,偏差e取最小值,此时可知tm时刻, For example, at time t1 , when the sequence to be processed (comparison) of the training set rotates N measurement points, the deviation e takes the minimum value, at this time It can be seen that at time t m ,
n个测点的旋转修正后的排温数据序列T4',1,T4',2,…,T4',n具体为:The exhaust temperature data series T 4 ' ,1 , T 4 ' ,2 ,..., T 4 ' ,n of n measuring points after rotation correction are specifically:
……...
其中分别为第1个测点t1时刻到tm时刻旋转修正后的温度值,分别为第2个测点t1时刻到tm时刻旋转修正后的温度值,分别为第n个测点t1时刻到tm时刻旋转修正后的温度值。in are the temperature values after rotation correction of the first measuring point from time t 1 to time t m respectively, are the temperature values after rotation correction of the second measuring point from time t 1 to time t m respectively, Respectively, the temperature values after rotation correction of the nth measuring point from time t 1 to time t m .
其它步骤及参数与具体实施方式一或二相同。Other steps and parameters are the same as those in Embodiment 1 or Embodiment 2.
具体实施方式四:本实施方式与具体实施方式一至三之一不同的是:所述步骤三中确定n个测点的旋转修正后的排温数据序列T4',1,T4',2,…,T4',n与n个测点的平均温度之间的关系具体为:Embodiment 4: The difference between this embodiment and one of Embodiments 1 to 3 is that in the step 3, the exhaust temperature data sequence T 4 ' ,1 and T 4 ' ,2 after rotation correction of n measuring points are determined ,..., the relationship between T 4 ' ,n and the average temperature of n measuring points is specifically:
在燃烧室中,其出口排温的分布情况反映了燃烧室的运转情况。同时,为了消除不同工况对不同测点温度的缩放,利用平均温度进行计算。出口处各个测点温度满足线性关系,即对于每个测点而言,其旋转修正后的排温与平均温度之间存在线性关系:In the combustion chamber, the distribution of the outlet exhaust temperature reflects the operation of the combustion chamber. At the same time, in order to eliminate the scaling of different measuring point temperatures under different working conditions, the average temperature is used for calculation. The temperature of each measuring point at the outlet satisfies a linear relationship, that is, for each measuring point, there is a linear relationship between the exhaust temperature after rotation correction and the average temperature:
其中T′4,i为第i个测点的旋转修正后的排温数据序列,i=[1,2,…,n];为n个测点的平均温度,αi为平均温度拟合参数,βi为常数项拟合参数;Where T′ 4,i is the exhaust temperature data sequence after rotation correction of the i-th measuring point, i=[1,2,…,n]; is the average temperature of n measuring points, α i is the average temperature fitting parameter, and β i is the constant term fitting parameter;
i为排温测点号,总共为n个测点。αi,βi为拟合参数,αi,βi确定方法如下:设燃气轮机在排气端均匀分布着n个热电偶,对于任一测点i,i=1,2,…n选取机组正常运行时一段时间内的数据作为训练集;i is the exhaust temperature measuring point number, and there are n measuring points in total. α i , β i are the fitting parameters, and the determination method of α i , β i is as follows: Assume that the gas turbine has n thermocouples uniformly distributed at the exhaust end, and for any measuring point i, i=1,2,...n select the unit The data during normal operation for a period of time is used as the training set;
用该时间段内各个时间点的由最小二乘法确定该测点处对应的参数αi,βi的值。at each time point in the time period The values of the parameters α i and β i corresponding to the measuring point are determined by the least square method.
其它步骤及参数与具体实施方式一至三之一相同。Other steps and parameters are the same as those in Embodiments 1 to 3.
具体实施方式五:本实施方式与具体实施方式一至四之一不同的是:所述步骤四中得到n个测点的温度预测值T4pre,1,T4pre,2,…,T4pre,i,…,T4pre,n的计算公式为:Embodiment 5: The difference between this embodiment and Embodiment 1 to Embodiment 4 is that the temperature prediction values T 4pre,1 , T 4pre,2 ,..., T 4pre,i of n measuring points are obtained in the step 4 ,...,T 4pre,n is calculated as:
根据步骤三的函数关系将t时刻正常运行条件下的平均排温代入,可以确定各个测点的温度预测值:According to the functional relationship in step 3, the average exhaust temperature under normal operating conditions at time t is calculated as Substituting in, the temperature prediction value of each measuring point can be determined:
其它步骤及参数与具体实施方式一至四之一相同。Other steps and parameters are the same as in one of the specific embodiments 1 to 4.
具体实施方式六:本实施方式与具体实施方式一至五之一不同的是:所述步骤五中根据步骤四得到的n个测点的温度预测值T4pre,1,T4pre,2,…,T4pre,i,…,T4pre,n,得到每一个测点的误差带ΔT1,ΔT2,…,ΔTi,…,ΔTn的计算公式为:Embodiment 6: The difference between this embodiment and one of Embodiments 1 to 5 is that the temperature prediction values T 4pre,1 , T 4pre,2 ,..., T 4pre,i ,…,T 4pre,n , the calculation formula of the error band ΔT 1 ,ΔT 2 ,…,ΔT i ,…,ΔT n of each measuring point is:
ΔTi=T4pre,i-T′4,1t ΔT i =T 4pre,i -T′ 4,1t
其中T′4,1t为第1个测点t时刻旋转修正后的温度值,t=[t1,t2,…,tm]。Where T′ 4,1t is the temperature value after rotation correction at the first measuring point at time t, t=[t 1 ,t 2 ,…,t m ].
可以求出t时刻机组无故障时测点i的理论值T4pre,i与测点i的实测值T′4,1t之差为ΔTi,ΔTi=T4pre,i-T′4,1t。将这一误差称为i测点的误差带。在计算时,把每个时刻的排温数据均代入到步骤四中,得到机组无故障时各个测点的排温预测值的时间序列,再利用ΔTi=T4pre,i-T′4,1t计算测点的误差带。The difference between the theoretical value T 4pre,i of measuring point i and the measured value T′ 4,1t of measuring point i at time t can be calculated as ΔT i , ΔT i =T 4pre,i -T′ 4,1t . This error is called the error band of the i measuring point. When calculating, the exhaust temperature data at each moment are substituted into step 4 to obtain the time series of exhaust temperature prediction values of each measuring point when the unit has no faults, and then use ΔT i =T 4pre,i -T′ 4, 1t calculates the error band of the measuring point.
其它步骤及参数与具体实施方式一至五之一相同。Other steps and parameters are the same as one of the specific embodiments 1 to 5.
本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。The present invention can also have other various embodiments, without departing from the spirit and essence of the present invention, those skilled in the art can make various corresponding changes and deformations according to the present invention, but these corresponding changes and deformations are all Should belong to the scope of protection of the appended claims of the present invention.
Claims (6)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810646418.8A CN108827643B (en) | 2018-06-21 | 2018-06-21 | Gas turbine high-temperature component fault early warning method considering exhaust temperature field rotation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810646418.8A CN108827643B (en) | 2018-06-21 | 2018-06-21 | Gas turbine high-temperature component fault early warning method considering exhaust temperature field rotation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108827643A true CN108827643A (en) | 2018-11-16 |
CN108827643B CN108827643B (en) | 2020-04-07 |
Family
ID=64143034
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810646418.8A Active CN108827643B (en) | 2018-06-21 | 2018-06-21 | Gas turbine high-temperature component fault early warning method considering exhaust temperature field rotation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108827643B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109556876A (en) * | 2018-11-07 | 2019-04-02 | 国网浙江省电力有限公司电力科学研究院 | A kind of diagnostic method for distinguishing gas turbine combustion failure and passage of heat equipment fault |
CN110672332A (en) * | 2019-09-10 | 2020-01-10 | 上海电力大学 | A Gas Turbine Fault Early Warning System Based on SARIMA Model |
CN111639401A (en) * | 2020-05-13 | 2020-09-08 | 中国航发贵阳发动机设计研究所 | Method for calculating turbine front temperature field by using turbine rear temperature field |
CN112116157A (en) * | 2020-09-18 | 2020-12-22 | 广东能源集团科学技术研究院有限公司 | Method and system for constructing combustion state comprehensive early warning system of gas turbine |
CN112460634A (en) * | 2020-11-23 | 2021-03-09 | 西安热工研究院有限公司 | Method for determining fault combustion chamber of gas turbine |
CN116380474A (en) * | 2023-02-20 | 2023-07-04 | 中国船舶重工集团公司第七0三研究所 | A Method for Predicting Average Temperature of Marine Gas Turbine Power Turbine Inlet Based on Energy Balance Method |
CN117171517A (en) * | 2023-11-02 | 2023-12-05 | 无锡尚航数据有限公司 | Dynamic early warning method for operation fault risk of data center |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0632268A1 (en) * | 1993-07-01 | 1995-01-04 | Johnson Service Company | Apparatus and method for determining the indoor air quality within an enclosed space |
CN105067275A (en) * | 2015-07-24 | 2015-11-18 | 哈尔滨工业大学 | Gas turbine combustion system online monitoring and diagnosis method based on exhaust temperature deviation index |
CN105134386A (en) * | 2015-09-02 | 2015-12-09 | 哈尔滨工业大学 | On-line monitoring method for gas turbine combustion system based on measuring-point weighted value |
-
2018
- 2018-06-21 CN CN201810646418.8A patent/CN108827643B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0632268A1 (en) * | 1993-07-01 | 1995-01-04 | Johnson Service Company | Apparatus and method for determining the indoor air quality within an enclosed space |
CN105067275A (en) * | 2015-07-24 | 2015-11-18 | 哈尔滨工业大学 | Gas turbine combustion system online monitoring and diagnosis method based on exhaust temperature deviation index |
CN105134386A (en) * | 2015-09-02 | 2015-12-09 | 哈尔滨工业大学 | On-line monitoring method for gas turbine combustion system based on measuring-point weighted value |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109556876A (en) * | 2018-11-07 | 2019-04-02 | 国网浙江省电力有限公司电力科学研究院 | A kind of diagnostic method for distinguishing gas turbine combustion failure and passage of heat equipment fault |
CN109556876B (en) * | 2018-11-07 | 2020-09-04 | 国网浙江省电力有限公司电力科学研究院 | A diagnostic method for distinguishing combustion faults of gas turbines and hot aisle equipment faults |
CN110672332A (en) * | 2019-09-10 | 2020-01-10 | 上海电力大学 | A Gas Turbine Fault Early Warning System Based on SARIMA Model |
CN111639401A (en) * | 2020-05-13 | 2020-09-08 | 中国航发贵阳发动机设计研究所 | Method for calculating turbine front temperature field by using turbine rear temperature field |
CN111639401B (en) * | 2020-05-13 | 2023-10-27 | 中国航发贵阳发动机设计研究所 | Method for calculating temperature field before turbine by using temperature field after turbine |
CN112116157A (en) * | 2020-09-18 | 2020-12-22 | 广东能源集团科学技术研究院有限公司 | Method and system for constructing combustion state comprehensive early warning system of gas turbine |
CN112116157B (en) * | 2020-09-18 | 2023-05-30 | 广东能源集团科学技术研究院有限公司 | Method and system for constructing comprehensive early warning system of combustion state of gas turbine |
CN112460634A (en) * | 2020-11-23 | 2021-03-09 | 西安热工研究院有限公司 | Method for determining fault combustion chamber of gas turbine |
CN116380474A (en) * | 2023-02-20 | 2023-07-04 | 中国船舶重工集团公司第七0三研究所 | A Method for Predicting Average Temperature of Marine Gas Turbine Power Turbine Inlet Based on Energy Balance Method |
CN116380474B (en) * | 2023-02-20 | 2025-07-15 | 中国船舶重工集团公司第七0三研究所 | A method for predicting the average temperature at the inlet of a marine gas turbine based on the energy balance method |
CN117171517A (en) * | 2023-11-02 | 2023-12-05 | 无锡尚航数据有限公司 | Dynamic early warning method for operation fault risk of data center |
CN117171517B (en) * | 2023-11-02 | 2024-01-26 | 无锡尚航数据有限公司 | Dynamic early warning method for operation fault risk of data center |
Also Published As
Publication number | Publication date |
---|---|
CN108827643B (en) | 2020-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108827643A (en) | A kind of high-temperature component of gas turbine fault early warning method for considering to arrange warm temperature field rotation | |
US9650909B2 (en) | Multi-stage compressor fault detection and protection | |
US8770913B1 (en) | Apparatus and process for rotor creep monitoring | |
CN105241669B (en) | On-line Monitoring Method of Gas Turbine Combustion System Based on Comparison Coding | |
US6526358B1 (en) | Model-based detection of leaks and blockages in fluid handling systems | |
JP6628734B2 (en) | Method of monitoring abnormal combustion in gas turbomachine and gas turbomachine including abnormal combustion detection system | |
US20090228230A1 (en) | System and method for real-time detection of gas turbine or aircraft engine blade problems | |
WO2012084453A1 (en) | Method of detecting a predetermined condition in a gas turbine and failure detection system for a gas turbine | |
US20210131354A1 (en) | Methods and systems for detection of control sensor override | |
US20130104516A1 (en) | Method of monitoring an operation of a compressor bleed valve | |
CN106525442B (en) | Method and device for monitoring gas path performance of gas turbine | |
CN109556876B (en) | A diagnostic method for distinguishing combustion faults of gas turbines and hot aisle equipment faults | |
CN105067275A (en) | Gas turbine combustion system online monitoring and diagnosis method based on exhaust temperature deviation index | |
RU2636602C2 (en) | Method for monitoring engine start cycle of gas-turbine plant | |
RU2733150C2 (en) | Clogging control in blower circuit of starter nozzle for turbomachine | |
CN105510045B (en) | Gas turbine combustion system on-line monitoring method based on coefficient matrix | |
CN105134386B (en) | On-line monitoring method for gas turbine combustion system based on measuring-point weighted value | |
US20120141251A1 (en) | Method and device for predicting the instability of an axial compressor | |
CN105114977B (en) | Gas turbine combustion system online monitoring method based on exhaust temperature measuring point correlation | |
CN108518285B (en) | Gas turbine temperature monitoring method and gas turbine temperature monitoring system | |
JP5164928B2 (en) | Gas turbine abnormality diagnosis device | |
JP7682723B2 (en) | Control of a power generation system during online maintenance using multiple maintenance modes | |
JP2024070523A (en) | Temperature measurement system and temperature measurement method | |
CN115144186A (en) | A continuous high-precision diagnosis method for gas path faults of gas turbine engine | |
CN115045807A (en) | A method for detecting abnormal rotation speed of wind turbine generators |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210408 Address after: Room 206-10, building 16, 1616 Chuangxin Road, Songbei District, Harbin City, Heilongjiang Province Patentee after: Harbin jizuo technology partnership (L.P.) Patentee after: Harbin Institute of Technology Asset Management Co.,Ltd. Address before: 150001 No. 92 West straight street, Nangang District, Heilongjiang, Harbin Patentee before: HARBIN INSTITUTE OF TECHNOLOGY Patentee before: NANJING POWER HORIZON INFORMATION TECHNOLOGY Co.,Ltd. |
|
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210611 Address after: Room 206-12, building 16, 1616 Chuangxin Road, Songbei District, Harbin City, Heilongjiang Province Patentee after: Harbin Institute of Technology Institute of artificial intelligence Co.,Ltd. Address before: Room 206-10, building 16, 1616 Chuangxin Road, Songbei District, Harbin City, Heilongjiang Province Patentee before: Harbin jizuo technology partnership (L.P.) Patentee before: Harbin Institute of Technology Asset Management Co.,Ltd. |