CN117153282A - Method, device, medium and equipment for calculating reference temperature of combustion chamber of gas turbine - Google Patents

Method, device, medium and equipment for calculating reference temperature of combustion chamber of gas turbine Download PDF

Info

Publication number
CN117153282A
CN117153282A CN202311164585.6A CN202311164585A CN117153282A CN 117153282 A CN117153282 A CN 117153282A CN 202311164585 A CN202311164585 A CN 202311164585A CN 117153282 A CN117153282 A CN 117153282A
Authority
CN
China
Prior art keywords
candidate
combustion chamber
indexes
preset
mapping relation
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
Application number
CN202311164585.6A
Other languages
Chinese (zh)
Other versions
CN117153282B (en
Inventor
李明
牟影
杨猛
蒋佳华
尹峰
高家春
韩永军
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.)
China Hangfa Gas Turbine Co ltd
Original Assignee
China Hangfa Gas Turbine Co ltd
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 China Hangfa Gas Turbine Co ltd filed Critical China Hangfa Gas Turbine Co ltd
Priority to CN202311164585.6A priority Critical patent/CN117153282B/en
Publication of CN117153282A publication Critical patent/CN117153282A/en
Application granted granted Critical
Publication of CN117153282B publication Critical patent/CN117153282B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/10—Analysis or design of chemical reactions, syntheses or processes
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/70—Machine learning, data mining or chemometrics

Landscapes

  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Health & Medical Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Analytical Chemistry (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Control Of Turbines (AREA)

Abstract

The application discloses a method, a device, a medium and equipment for calculating the reference temperature of a combustion chamber of a gas turbine, comprising the following steps: respectively acquiring actually measured reference temperatures of the combustion chamber at a plurality of first preset moments, and working parameters of a plurality of strong correlation indexes and a plurality of weak correlation indexes; acquiring candidate index sets from the strong association indexes and the weak association indexes; constructing a candidate mapping relation between the reference temperature of the combustion chamber and the candidate index group; and determining a deviation range of the candidate mapping relation, if the deviation range is out of the preset range, updating the candidate index group, and constructing an updated candidate mapping relation until the deviation range of the updated candidate mapping relation is in the preset range, so as to obtain the target mapping relation. In different environments, the index group adopted for constructing the target mapping relation is not fixed, the obtained target mapping relation better accords with the current environment, and the reference temperature error obtained by calculating the target mapping relation is smaller.

Description

Method, device, medium and equipment for calculating reference temperature of combustion chamber of gas turbine
Technical Field
The application relates to the field of gas turbines, in particular to a method, a device, a medium and equipment for calculating the reference temperature of a combustion chamber of a gas turbine.
Background
In the field of gas turbines, pollutants such as carbon monoxide (CO) and nitrogen oxides (NOx) are emitted when fuel in a combustion chamber is combusted, and with the increasing strictness of environmental regulations, it is necessary to strictly control the pollutant equivalent amount of the emissions.
The emission equivalent of pollutants such as carbon monoxide (CO) and nitrogen oxides (NOx) is directly related to the flame temperature during combustion, so the emission equivalent of pollutants can be controlled by controlling the flame temperature (reference temperature) at the time of combustion in the combustion chamber. However, the flame temperature of the combustion chamber of the gas turbine is substantially 1300K or more, and in this high temperature environment, it is difficult for the temperature sensor to directly measure the flame temperature of the combustion chamber for a long period of time.
Currently, there are some methods for calculating the flame temperature of the combustion chamber by fitting, which mostly adopt actually measured and fixed index parameters capable of influencing the flame temperature of the combustion chamber, so as to fit and calculate the flame temperature of the combustion chamber. However, there are many factors that affect the flame temperature of the combustor, and the operating environment of the gas turbine is complex and variable, such as: for the offshore operation environment, the land operation environment and the air operation environment, the influences of various factors on the flame temperature of the combustion chamber of the gas turbine are different, and the fitting indexes in the prior art are relatively fixed, so that the fitting cannot be performed under different environments.
Disclosure of Invention
The application provides a method, a device, a medium and equipment for calculating the reference temperature of a combustion chamber of a gas turbine, and aims to solve the problems that the temperature is higher when the combustion chamber of the gas turbine burns, the flame temperature of the combustion chamber is not easy to directly measure by using a temperature measuring sensor, and the error of the conventional fitting calculation method is larger.
In an embodiment of the present application, a method for calculating a reference temperature of a combustion chamber of a gas turbine is provided, including:
respectively obtaining actually measured reference temperatures of the combustion chamber at a plurality of first preset moments;
respectively acquiring working parameters of a plurality of strong association indexes and a plurality of weak association indexes of the combustion chamber at a plurality of first preset moments, wherein the plurality of strong association indexes comprise: turbine exhaust temperature, fuel flow, compressor outlet pressure, and atmospheric temperature; the plurality of weakly associated candidate indexes includes: compressor outlet temperature, barometric pressure, air intake flow, and turbine exhaust pressure;
obtaining a candidate index group from the strong correlation indexes and the weak correlation indexes, wherein the candidate index group comprises all the strong correlation indexes and at least one weak correlation index;
constructing a candidate mapping relation between the reference temperature of the combustion chamber and the candidate index group based on working parameters of each index in the candidate index group at each first preset time and the actually measured reference temperature of the combustion chamber at each first preset time;
Determining a deviation range of the candidate mapping relation, if the deviation range is out of a preset range, updating the candidate index group, and constructing an updated candidate mapping relation based on working parameters of each index in the updated candidate index group at each first preset time and actual measurement reference temperatures of the combustion chamber at each first preset time until the deviation range of the updated candidate mapping relation is in the preset range;
and taking the candidate mapping relation with the deviation range within the preset range as a target mapping relation of the reference temperature of the combustion chamber.
In an embodiment of the present application, the updating the candidate index set includes:
changing the number of weak association indexes in the candidate index group; or,
and changing weak association indexes in the candidate index group.
In the embodiment of the application, the number of the plurality of first preset moments is not less than four, and the plurality of first preset moments comprise the initial combustion moment of the combustion chamber, the stable moment after stable combustion and the moment between the initial combustion moment and the stable moment.
In the embodiment of the present application, the constructing the candidate mapping relationship of the reference temperature of the combustion chamber based on the working parameters of each index in the candidate index set at each first preset time and the actually measured reference temperature of the combustion chamber at each first preset time includes:
Constructing a parameter matrix based on working parameters of each index in the candidate index group at each first preset time;
constructing a temperature matrix based on the actually measured reference temperatures of the combustion chamber at each first preset time;
determining an influence factor of each index in the candidate index group based on the parameter matrix and the temperature matrix;
and constructing a candidate mapping relation of the reference temperature of the combustion chamber based on the influence factors of the indexes in the candidate index group.
In the embodiment of the present application, the candidate mapping relationship is as follows:
CRT＝(K 1 * 1 + 2 * 2 + 3 * 3 + 4 * 4 )+(K 5 * 1 +…+ n-1 *
Yn-5+Kn*；
wherein the CRT represents the reference temperature of the combustion chamber, X 1 To X 4 Work parameters, K, representing the respective strongly associated index 1 To K 4 An influence factor representing each strong correlation index, Y 1 To Y n-5 Work parameters, K, representing the respective weakly associated index 5 To K n-1 Influence factors representing weak correlation indexes, C being a constant, K n Is the influencing factor of C.
In an embodiment of the present application, the determining the deviation range of the candidate mapping relationship includes:
determining verification indexes, wherein the verification indexes are strong association indexes and weak association indexes in candidate index groups corresponding to the candidate mapping relation;
Respectively acquiring working parameters of verification indexes of the combustion chamber at a plurality of preset second preset moments;
respectively obtaining actual measurement reference temperatures of the combustion chamber at a plurality of preset second preset moments;
determining a theoretical reference temperature of the combustion chamber at each second preset time based on the candidate mapping relation and working parameters of each verification index at each second preset time;
and determining the deviation range of the candidate mapping relation based on the theoretical reference temperature and the actual measurement reference temperature of the combustion chamber at each second preset time.
In the embodiment of the application, at least one second preset time is arranged between two adjacent first preset times.
In the embodiment of the application, the preset range is-1.08K to +0.49K.
The application also provides a reference temperature calculating device of the combustion chamber of the gas turbine, which comprises the following components:
the acquisition module is used for respectively acquiring the actually measured reference temperatures of the combustion chamber at a plurality of first preset moments; and
respectively acquiring working parameters of a plurality of strong association indexes and a plurality of weak association indexes of the combustion chamber at a plurality of first preset moments, wherein the plurality of strong association indexes comprise: turbine exhaust temperature, fuel flow, compressor outlet pressure, and atmospheric temperature; the plurality of weakly associated candidate indexes includes: compressor outlet temperature, barometric pressure, air intake flow, and turbine exhaust pressure;
The processing module is used for acquiring candidate index groups from the strong correlation indexes and the weak correlation indexes, wherein the candidate index groups comprise all the strong correlation indexes and at least one weak correlation index;
constructing a candidate mapping relation between the reference temperature of the combustion chamber and the candidate index group based on working parameters of each index in the candidate index group at each first preset time and the actually measured reference temperature of the combustion chamber at each first preset time;
determining a deviation range of the candidate mapping relation, if the deviation range is out of a preset range, updating the candidate index group, and constructing an updated candidate mapping relation based on working parameters of each index in the updated candidate index group at each first preset time and actual measurement reference temperatures of the combustion chamber at each first preset time until the deviation range of the updated candidate mapping relation is in the preset range;
and taking the candidate mapping relation with the deviation range within the preset range as a target mapping relation of the reference temperature of the combustion chamber.
In an embodiment of the application, the processing module is configured to update the candidate index set by:
Changing the number of weak association indexes in the candidate index group; or,
and changing weak association indexes in the candidate index group.
In the embodiment of the application, the number of the plurality of first preset moments is not less than four, and the plurality of first preset moments comprise the initial combustion moment of the combustion chamber, the stable moment after stable combustion and the moment between the initial combustion moment and the stable moment.
In an embodiment of the application, the processing module is further configured to:
constructing a parameter matrix based on working parameters of each index in the candidate index group at each first preset time;
constructing a temperature matrix based on the actually measured reference temperatures of the combustion chamber at each first preset time;
determining an influence factor of each index in the candidate index group based on the parameter matrix and the temperature matrix;
and constructing a candidate mapping relation of the reference temperature of the combustion chamber based on the influence factors of the indexes in the candidate index group.
In the embodiment of the application, the candidate mapping relation is constructed by the following steps:
CRT＝(K 1 * 1 + 2 * 2 + 3 * 3 + 4 * 4 )+(K 5 * 1 +…+ n-1 *
Yn-5+Kn*；
wherein the CRT represents the reference temperature of the combustion chamber, X 1 To X 4 Work parameters, K, representing the respective strongly associated index 1 To K 4 An influence factor representing each strong correlation index, Y 1 To Y n-5 Work parameters, K, representing the respective weakly associated index 5 To K n-1 Influence factors representing weak correlation indexes, C being a constant, K n Is the influencing factor of C.
In an embodiment of the application, the processing module is further configured to:
determining verification indexes, wherein the verification indexes are strong association indexes and weak association indexes in candidate index groups corresponding to the candidate mapping relation;
respectively acquiring working parameters of verification indexes of the combustion chamber at a plurality of preset second preset moments;
respectively obtaining actual measurement reference temperatures of the combustion chamber at a plurality of preset second preset moments;
determining a theoretical reference temperature of the combustion chamber at each second preset time based on the candidate mapping relation and working parameters of each verification index at each second preset time;
and determining the deviation range of the candidate mapping relation based on the theoretical reference temperature and the actual measurement reference temperature of the combustion chamber at each second preset time.
In the embodiment of the application, at least one second preset time is arranged between two adjacent first preset times.
In the embodiment of the application, the preset range is-1.08K to +0.49K.
The application also proposes a computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, implements a method as claimed in any of the preceding claims.
The application also proposes a computing device comprising a processor for implementing the method of any of the preceding claims when executing a computer program stored in a memory.
According to the embodiment of the application, the indexes capable of influencing the reference temperature of the combustion chamber of the gas turbine are divided into the strong association indexes and the weak association indexes, the candidate index groups are selected based on the strong association indexes and the weak association indexes when the mapping relation is constructed, the mapping relation is constructed according to the parameters of all indexes in the candidate index groups, whether all indexes in the candidate index groups are correctly selected or not is determined according to the deviation range of the mapping relation, and the reference temperature of the combustion chamber can be accurately calculated based on the mapping relation constructed by the finally selected candidate index groups. In contrast, in different environments, the index sets adopted in the construction of the target mapping relation are not fixed, but in the current environment, different candidate index sets are selected, different candidate mapping relations are constructed according to the measured temperatures of the indexes and the measured temperatures of the combustion chamber, then the target mapping relation with the accuracy meeting the conditions is selected according to the preset deviation range, in the selection process, the indexes are screened, so that the index set corresponding to the target mapping relation is the index which is most suitable for the current environment, namely the target mapping relation is more suitable for the current environment, and the reference temperature error calculated by using the target mapping relation is smaller.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, and it is obvious that the drawings in the following description are only some embodiments of the present application, and other drawings may be obtained according to the structures shown in these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a step diagram of a method for calculating a reference temperature for a combustion chamber of a gas turbine in accordance with an embodiment of the present application;
FIG. 2 is a block diagram of a gas turbine combustor reference temperature calculation apparatus in accordance with an embodiment of the present application;
FIG. 3 is a block diagram of a medium in an embodiment of the application;
FIG. 4 is a block diagram of a computing device in an embodiment of the application.
The achievement of the objects, functional features and advantages of the present application will be further described with reference to the accompanying drawings, in conjunction with the embodiments.
Detailed Description
The following description of the embodiments of the present application will be made clearly and fully with reference to the accompanying drawings, in which it is evident that the embodiments described are only some, but not all embodiments of the application. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.
It should be noted that, if directional indications (such as up, down, left, right, front, and rear … …) are included in the embodiments of the present application, the directional indications are merely used to explain the relative positional relationship, movement conditions, etc. between the components in a specific posture (as shown in the drawings), and if the specific posture is changed, the directional indications are correspondingly changed.
In addition, if there is a description of "first", "second", etc. in the embodiments of the present application, the description of "first", "second", etc. is for descriptive purposes only and is not to be construed as indicating or implying a relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include at least one such feature. In addition, the technical solutions of the embodiments may be combined with each other, but it is necessary to base that the technical solutions can be realized by those skilled in the art, and when the technical solutions are contradictory or cannot be realized, the combination of the technical solutions should be considered to be absent and not within the scope of protection claimed in the present application.
As shown in fig. 1, in the embodiment of the present application, the method for calculating the reference temperature of the combustion chamber of the gas turbine includes the following steps S100 to S600:
Step S100: and respectively acquiring actual measurement reference temperatures of the combustion chamber at a plurality of first preset moments.
In the embodiment of the application, a temporary high temperature measuring device can be installed in the combustion chamber of the gas turbine, and the flame temperature of the combustion chamber (namely, the actual reference temperature) is actually measured at a plurality of first preset moments in the debugging stage of the combustion chamber based on the temporary high temperature measuring device.
The first preset times may be times in a period from when the combustion chamber of the gas turbine starts to operate to when the operation is stable, or times in a period from the lowest load to the highest load of the gas turbine.
In addition, each of the first preset times may have an interval therebetween, such as a time at which the gas turbine combustor starts to operate, and a 2 nd, 6 th, 10 th minute or the like after the start of operation. The application does not limit the number of the first preset moments and does not limit the time interval between two adjacent first preset moments.
In addition, the temporary high temperature measuring device can at least ensure that the gas turbine is not damaged in the time period from the lowest load to the highest load or at least ensure that the gas turbine is not damaged in the time period from the start of the operation of the combustion chamber to the stable operation. The actual flame temperature of the combustion chamber at each first preset moment, i.e. the actual measured reference temperature, can be measured and obtained by means of the temporary high temperature measuring device.
Step S200: and respectively acquiring working parameters of a plurality of strong correlation indexes and a plurality of weak correlation indexes of the combustion chamber at a plurality of first preset moments.
In an embodiment of the present application, the plurality of strong association indexes includes: turbine Exhaust Temperature (TET), fuel Flow (FF), compressor outlet pressure (CPD), and Atmospheric Temperature (ATEM);
the plurality of weakly associated candidate metrics includes: compressor outlet temperature (CTD), barometric pressure (BAPO), air intake flow (MAF), and Turbine Exhaust Pressure (TEP).
Among them, the inventors found that the degree of correlation of the temperature of each of the above strong correlation indexes with the temperature of the combustion chamber flame is large, while the degree of correlation of the temperature of the weak correlation index with the temperature of the combustion chamber flame is small. Furthermore, the inventors have found that the degree of influence of each index on the temperature of the combustion chamber flame is not the same in different operating environments of the gas turbine. For example, in working environments such as sea, high altitude areas, extremely cold areas, high temperature desert areas and the like, the atmospheric pressure of the weak correlation index has a larger influence on the flame temperature of the combustion chamber, while in working environments such as low altitude land, the atmospheric pressure basically has no influence on the flame temperature of the combustion chamber, and the air inlet flow and the turbine exhaust pressure have a larger influence on the flame temperature of the combustion chamber. And for the strong correlation index, the flame temperature of the combustion chamber is greatly influenced under any working environment. Therefore, it is obviously inaccurate to select different indexes to construct the mapping relationship between the reference temperature of the combustion chamber and each index in different working environments if a fixed index is selected uniformly for constructing the mapping relationship in any environment.
Therefore, when the reference temperature of the combustion chamber is calculated by fitting based on the above-mentioned indexes, all the indexes may need to be considered in some working environments, and only some indexes may need to be considered in some environments, if the indexes are not distinguished, only a certain index is fixedly used for fitting calculation, and the calculated reference temperature of the combustion chamber has larger deviation compared with the actual reference temperature.
Step S300: and acquiring a candidate index group from the strong correlation indexes and the weak correlation indexes, wherein the candidate index group comprises all the strong correlation indexes and at least one weak correlation index.
In step S200, the strong correlation index and the weak correlation index have been explained, and have different degrees of influence on the standard temperature of combustion of the gas turbine under different working environments, so in the embodiment of the present application, the candidate index group may be selected from the strong correlation index and the weak correlation index, and by verifying the candidate index, it is determined whether the selected index is correct.
The candidate index set includes all strong association indexes and at least one weak association index, that is, the weak association index may include only one weak association index, or may include two, three or all weak association indexes. As shown in table 1 below, the candidate index set may include the following combinations:
TABLE 1
Step S400: and constructing a candidate mapping relation of the reference temperature of the combustion chamber based on the working parameters of each index in the candidate index group at each first preset time and the actually measured reference temperature of the combustion chamber at each first preset time.
In the embodiment of the present application, the candidate mapping relationship may be constructed by the following steps S410 to S440:
step S410: and constructing a parameter matrix based on the working parameters of each index in the candidate index group at each first preset time.
Step S420: and constructing a temperature matrix based on the actually measured reference temperatures of the combustion chamber at each first preset time.
In an embodiment of the present application, as shown in Table 2 below, the parameters of each index from the start of operation to the time when the load reaches a maximum, and the measured reference temperatures of the gas turbine of Table 2 are as follows:
TABLE 2
Wherein N1 is the initial time, and N11 is the time when the load reaches the maximum;
TET (1) -TET (11), representing the respective exhaust temperatures of the turbine at 11 times;
FF (1) -FF (11), representing the fuel flow at 11 hours;
CPD (1) -CPD (11), representing the respective compressor outlet pressures at 11 times;
CTD (1) -CTD (11), representing the respective compressor outlet temperatures at 11 times;
ATEM (1) -ATEM (11), representing the respective atmospheric temperatures at 11 times;
BAPO (1) -BAPO (11), representing the respective atmospheric pressures at 11 times;
MAF (1) -MAF (11), representing the respective air intake flows at 11 times;
TEP (1) -TEP (11), representing the respective exhaust pressures at 11 times;
CPT (1) -CPT (1), represent each measured reference temperature at 11 time instants.
In the embodiment of the application, the number of the first preset moments is not less than four, and the plurality of first preset moments comprise the initial combustion moment of the combustion chamber, the stable moment after stable combustion and the moment between the initial combustion moment and the stable moment. Because the first preset time includes the initial combustion time, the stable time and the working parameters of each index in the candidate index group at other times between the initial combustion time and the stable time, the mapping relationship obtained by constructing the candidate index group based on each parameter at each first preset time can be used for calculating the reference temperature of the combustion chamber at the initial combustion time and the stable time and the reference temperature of the combustion chamber at other times between the initial combustion time and the stable time.
It should be noted that, if the mapping relationship only needs to be able to calculate the reference temperature of a certain section of the combustion section, for example, only needs to calculate the reference temperature between N3 and N9 by using the mapping relationship, then the first preset time may be selected only between N3 and N9, and does not need to include the start time and the stabilization time.
In the embodiment of the present application, N1, N4, N7, N9, N11 may be selected as each first preset time, assuming that the reference temperature of the entire combustion process of the combustion chamber of the gas turbine needs to be calculated using the map. Assume that the candidate index set is the candidate index set 6 shown in the above table 1: TET, FF, CPD, ATEM, CTD, MAF, then in step S410, the following parameter matrix may be constructed:
wherein "C" in the parameter matrix is a constant.
In step S420, the following temperature matrix may be constructed based on the measured reference temperatures at the first preset moments of N1, N4, N7, N9, and N11:
step S430: and determining the influence factors of the indexes in the candidate index group based on the parameter matrix and the temperature matrix.
The influence factors of the indexes can reflect the influence degree of the indexes on the reference temperature of the combustion chamber in the working environment, and in the embodiment of the application, the magnitude of each influence factor can be calculated based on the parameter matrix and the temperature matrix.
For example, the parameter matrix may be inverted based on a least squares method as follows:
after inverting the parameter matrix, the temperature matrix is multiplied by the inverse of the parameter matrix as follows:
the influence factors of the indexes and the constant C can be obtained.
Step S440: and constructing a candidate mapping relation of the reference temperature of the combustion chamber based on the influence factors of the indexes in the candidate index group.
Assume that the impact factors of the index TET, FF, CPD, CPT, CTD, MAF, C calculated in step S430 are K 1 、K 2 、K 3 、K 4 、K 5 、K 6 The influence factors of the constant C are K respectively 7 Then the following candidate mapping relationship can be constructed:
CRT＝K 1 *TET+K 2 *FF+K 3 *CPD+K 4 *CPT+K 5 CTD+K 6 MAF+K 7 *C
the formula is a candidate mapping relationship constructed according to the strong association index TET, FF, CPD, CPT and the weak association indexes CTD and MAF.
Step S500: and determining a deviation range of the candidate mapping relation, if the deviation range is out of a preset range, updating the candidate index group, and constructing an updated candidate mapping relation based on the working parameters of each index in the updated candidate index group at each first preset time and the actually measured reference temperature of the combustion chamber at each first preset time until the deviation range of the updated candidate mapping relation is in the preset range.
In the embodiment of the present application, the deviation range of the mapping relationship may be calculated based on the following steps S510 to S550, as follows:
step S510: and determining verification indexes, wherein the verification indexes are strong association indexes and weak association indexes in candidate index groups corresponding to the candidate mapping relations.
In the embodiment shown in the above step S400, the candidate index set includes TET, FF, CPD, CPT and the weakly associated indexes CTD and MAF, and then the verification index at the time of verification is still TET, FF, CPD, CPT and CTD and MAF.
Step S520: and respectively acquiring working parameters of the verification index of the combustion chamber at a plurality of preset second preset moments.
It should be noted that, the mapping relationship is constructed by using index parameters of each of N1, N4, N7, N9, and N11 in the N1-N11 moments, so that the second preset moment can be acted by the moment in N1-N11 during verification, and it should be noted that parameters of each of N1, N4, N7, N9, and N11 are already involved in constructing the mapping relationship, and during verification, verification moments can be selected from other moments except N1, N4, N7, N9, and N11, for example: time N3, N6, N8, etc.
In addition, the above-mentioned mapping relationship is constructed by fitting based on a plurality of times between N1 to N11, and the combustion states of the combustion chambers represented by different times are also different, and the mapping relationship needs to satisfy the requirement that the reference temperature of the combustion chamber in each state can be calculated, so that when the deviation is verified, each state of the combustion chamber needs to be verified. For example, in the embodiment of the present application, at least one second preset time is set between two first preset times.
For example, the second preset time includes: n3, N5, N8 and N10, wherein N3 is positioned between N1 and N4, verification is performed based on the moment N3, whether the reference temperature of the combustion chamber at the moment N1 and N4 is accurate or not can be calculated based on the mapping relation in a certain degree of reaction, and similarly, the verification moment N5 is used for calculating whether the reference temperature of the combustion chamber at the moment N4 and N7 is accurate or not based on the mapping relation in a certain degree of reaction; the verification moment of N8 can reflect whether the reference temperature of the combustion chamber is accurate or not based on the mapping relation under the moments of N7 and N9 to a certain extent; the verification time of N10 can reflect to a certain extent whether the reference temperature of the combustion chamber at the times N9 and N11 is accurate or not calculated based on the map.
In addition, when the verification time is selected, a plurality of second preset times may be selected between two adjacent first preset times, for example, N2 and N3 between N1 and N4 are both selected as verification times, and N5 and N6 between N4 and N7 are both selected as verification times.
Step S530: and respectively acquiring the actually measured reference temperatures of the combustion chamber at the plurality of preset second preset moments.
In the embodiment of the present application, it is assumed that the second preset time includes: as is clear from table 2, the measured reference temperatures at the times of N2, N3, N5, N6, N8, and N10 are classified into CRT (2), CRT (3), CRT (5), CRT (6), CRT (8), and CRT (10).
Step S540: and determining the theoretical reference temperature of the combustion chamber at each second preset time based on the candidate mapping relation and the working parameters of each verification index at each second preset time.
In the embodiment of the application, the theoretical reference temperature at each second preset time can be calculated and obtained based on the candidate mapping relation brought into by the parameters of each verification index at each time of N2, N3, N5, N6, N8 and N10.
Taking the time of N2 as an example, the theoretical reference temperature CRT (2) at the time of N2'
CRT(2)′＝K 1 *TET(2)+K 2 *FF(2)+K 3 *CPD(2)+K 4 *CPT(2)+K 5 CTD(2)+K 6 MAF(2)+K 7 *C
Similarly, CRT (2) ', CRT (3)', CRT (5) ', CRT (6)', CRT (8) ', CRT (10)', were calculated in this order.
Step S550: and determining the deviation range of the candidate mapping relation based on the theoretical reference temperature and the actual measurement reference temperature of the combustion chamber at each second preset time.
In the embodiment of the application, the deviation between the theoretical reference temperature and the actually measured reference temperature calculated by using the mapping relation at each second preset moment can be obtained by subtracting the actual reference temperature from the theoretical reference temperature at each second preset moment.
For example, regarding the combustion state of the combustion chamber at the time of N2, the deviation of the mapping relationship in each combustion state can be calculated by using the difference between the CRTs (2)' and (2) as the deviation of the mapping relationship in the combustion state of the combustion chamber, and the same can be said to obtain the deviation of the mapping relationship in each combustion state at each second preset time. And the range of the maximum value and the lowest value of the deviation at each second preset time is the deviation range of the mapping relation.
After the deviation range of the mapping relationship is obtained, the deviation range can be compared with a preset range, for example, in the embodiment of the application, the preset range is between-1.08K and +0.49K, and if the deviation range of the candidate mapping relationship is between-1.08K and +0.49K, the mapping relationship obtained by fitting based on each index selected in step S300 is described, so that the reference temperature of the combustion chamber in the combustion state reacted at each moment can be accurately reflected. If the deviation range of the candidate mapping relationship exceeds-1.08K to +0.49K, it indicates that the reference temperature at each time calculated based on the mapping relationship obtained by fitting the indexes selected in step S300 cannot accurately reflect the actual reference temperature of the combustion chamber in the combustion state corresponding to the time at least at a certain time, and the candidate index set needs to be updated at this time.
In the embodiment of the application, when the candidate index group is updated, the number of weak association indexes in the candidate index group can be changed, for example, one weak association index is changed into two weak association indexes; or, changing the weak association index in the candidate index group, such as updating from CTD, TEP to BAPO, MAF.
The above embodiment selects the candidate index set 6 in table 1 at step S300, and at this time, any one of the candidate index sets 1 to 5 or the candidate index sets 7 to 15 may be replaced, and based on the updated candidate index set, the mapping relationship is reconstructed according to step S400, and based on steps S510 to S550, the deviation range of the constructed mapping relationship is re-verified until the deviation range of the constructed mapping relationship is within the preset range based on the final selected candidate index set.
Step S600: and taking the candidate mapping relation with the deviation range within the preset range as a target mapping relation of the reference temperature of the combustion chamber.
In step S500, a mapping relationship can be finally determined, and the deviation range thereof is within the preset deviation, which indicates that the reference temperature of the combustion chamber at each time in the period from the initial combustion to the highest load can be accurately calculated based on the mapping relationship, and then the mapping relationship can be used as the target mapping relationship of the reference temperature of the combustion chamber of the gas turbine, and the reference temperature of the combustion chamber at each time can be accurately calculated by using the target mapping relationship, so that the pollutant emission equivalent can be controlled by using the reference temperature.
In addition, in another embodiment, for the 15 candidate index sets in table 1, if the deviation ranges of the candidate mapping relationships constructed based on the plurality of index sets are all within the preset range, the final target mapping relationship may be determined according to the interval size of the deviation range of each candidate mapping relationship. For example, one of the candidate mappings having the smallest interval of the deviation range may be selected as the target mapping.
According to the embodiment of the application, the indexes capable of influencing the reference temperature of the combustion chamber of the gas turbine are divided into the strong association indexes and the weak association indexes, the candidate index groups are selected based on the strong association indexes and the weak association indexes when the mapping relation is constructed, the mapping relation is constructed according to the parameters of all indexes in the candidate index groups, whether all indexes in the candidate index groups are correctly selected or not is determined according to the deviation range of the mapping relation, and the reference temperature of the combustion chamber can be accurately calculated based on the mapping relation constructed by the finally selected candidate index groups. In contrast, in different environments, the index sets adopted in the construction of the target mapping relation are not fixed, but in the current environment, different candidate index sets are selected, different candidate mapping relations are constructed according to the measured temperatures of the indexes and the measured temperatures of the combustion chamber, then the target mapping relation with the accuracy meeting the conditions is selected according to the preset deviation range, in the selection process, the indexes are screened, so that the index set corresponding to the target mapping relation is the index which is most suitable for the current environment, namely the target mapping relation is more suitable for the current environment, and the reference temperature error calculated by using the target mapping relation is smaller.
Exemplary apparatus
Having described the method of an exemplary embodiment of the present application, an exemplary gas turbine combustor reference temperature calculation apparatus 100 of the present application, as shown in FIG. 2, is described next, which, in an example of the present application, comprises:
the acquisition module is used for respectively acquiring the actually measured reference temperatures of the combustion chamber at a plurality of first preset moments; and
respectively acquiring working parameters of a plurality of strong association indexes and a plurality of weak association indexes of the combustion chamber at a plurality of first preset moments, wherein the plurality of strong association indexes comprise: turbine exhaust temperature, fuel flow, compressor outlet pressure, and atmospheric temperature; the plurality of weakly associated candidate indexes includes: compressor outlet temperature, barometric pressure, air intake flow, and turbine exhaust pressure;
the processing module is used for acquiring candidate index groups from the strong correlation indexes and the weak correlation indexes, wherein the candidate index groups comprise all the strong correlation indexes and at least one weak correlation index;
constructing a candidate mapping relation between the reference temperature of the combustion chamber and the candidate index group based on working parameters of each index in the candidate index group at each first preset time and the actually measured reference temperature of the combustion chamber at each first preset time;
Determining a deviation range of the candidate mapping relation, if the deviation range is out of a preset range, updating the candidate index group, and constructing an updated candidate mapping relation based on working parameters of each index in the updated candidate index group at each first preset time and actual measurement reference temperatures of the combustion chamber at each first preset time until the deviation range of the updated candidate mapping relation is in the preset range;
and taking the candidate mapping relation with the deviation range within the preset range as a target mapping relation of the reference temperature of the combustion chamber.
In an embodiment of the application, the processing module is configured to update the candidate index set by:
changing the number of weak association indexes in the candidate index group; or,
and changing weak association indexes in the candidate index group.
In the embodiment of the application, the number of the plurality of first preset moments is not less than four, and the plurality of first preset moments comprise the initial combustion moment of the combustion chamber, the stable moment after stable combustion and the moment between the initial combustion moment and the stable moment.
In an embodiment of the application, the processing module is further configured to:
constructing a parameter matrix based on working parameters of each index in the candidate index group at each first preset time;
constructing a temperature matrix based on the actually measured reference temperatures of the combustion chamber at each first preset time;
determining an influence factor of each index in the candidate index group based on the parameter matrix and the temperature matrix;
and constructing a candidate mapping relation of the reference temperature of the combustion chamber based on the influence factors of the indexes in the candidate index group.
In the embodiment of the application, the candidate mapping relation is constructed by the following steps:
CRT＝(K 1 * 1 + 2 * 2 + 3 * 3 + 4 * 4 )+(K 5 * 1 +…+ n-1 *
Yn-5+Kn*；
wherein the CRT represents the reference temperature of the combustion chamber, X 1 To X 4 Work parameters, K, representing the respective strongly associated index 1 To K 4 An influence factor representing each strong correlation index, Y 1 To Y n-5 Work parameters, K, representing the respective weakly associated index 5 To K n-1 Influence factors representing weak correlation indexes, C being a constant, K n Is the influencing factor of C.
In an embodiment of the application, the processing module is further configured to:
determining verification indexes, wherein the verification indexes are strong association indexes and weak association indexes in candidate index groups corresponding to the candidate mapping relation;
Respectively acquiring working parameters of verification indexes of the combustion chamber at a plurality of preset second preset moments;
respectively obtaining actual measurement reference temperatures of the combustion chamber at a plurality of preset second preset moments;
determining a theoretical reference temperature of the combustion chamber at each second preset time based on the candidate mapping relation and working parameters of each verification index at each second preset time;
and determining the deviation range of the candidate mapping relation based on the theoretical reference temperature and the actual measurement reference temperature of the combustion chamber at each second preset time.
In the embodiment of the application, at least one second preset time is arranged between two adjacent first preset times.
In the embodiment of the application, the preset range is-1.08K to +0.49K.
The application also proposes a computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, implements a method as claimed in any of the preceding claims.
The application also proposes a computing device comprising a processor for implementing the method of any of the preceding claims when executing a computer program stored in a memory.
According to the embodiment of the application, the indexes capable of influencing the reference temperature of the combustion chamber of the gas turbine are divided into the strong association indexes and the weak association indexes, the processing module is utilized to select the candidate index groups based on the strong association indexes and the weak association indexes when the mapping relation is constructed, the mapping relation is constructed according to the parameters of each index in the candidate index groups, whether each index in the candidate index groups is correctly selected or not is determined according to the deviation range of the mapping relation, and the reference temperature of the combustion chamber can be accurately calculated based on the mapping relation constructed by the finally selected candidate index groups. In contrast, in different environments, the index sets adopted in the construction of the target mapping relation are not fixed, but in the current environment, different candidate index sets are selected, different candidate mapping relations are constructed according to the measured temperatures of the indexes and the measured temperatures of the combustion chamber, then the target mapping relation with the accuracy meeting the conditions is selected according to the preset deviation range, in the selection process, the indexes are screened, so that the index set corresponding to the target mapping relation is the index which is most suitable for the current environment, namely the target mapping relation is more suitable for the current environment, and the reference temperature error calculated by using the target mapping relation is smaller.
Exemplary Medium
Having described the methods and apparatus of exemplary embodiments of the present invention, a computer-readable storage medium of exemplary embodiments of the present invention is described next with reference to FIG. 3.
Referring to fig. 3, a computer readable storage medium is shown as an optical disc 70, on which a computer program (i.e., a program product) is stored, which when executed by a processor, implements the steps described in the above method embodiments, for example: respectively obtaining actually measured reference temperatures of the combustion chamber at a plurality of first preset moments; respectively acquiring working parameters of a plurality of strong association indexes and a plurality of weak association indexes of the combustion chamber at a plurality of first preset moments, wherein the plurality of strong association indexes comprise: turbine exhaust temperature, fuel flow, compressor outlet pressure, and atmospheric temperature; the plurality of weakly associated candidate indexes includes: compressor outlet temperature, barometric pressure, air intake flow, and turbine exhaust pressure; obtaining a candidate index group from the strong correlation indexes and the weak correlation indexes, wherein the candidate index group comprises all the strong correlation indexes and at least one weak correlation index; constructing a candidate mapping relation between the reference temperature of the combustion chamber and the candidate index group based on working parameters of each index in the candidate index group at each first preset time and the actually measured reference temperature of the combustion chamber at each first preset time; determining a deviation range of the candidate mapping relation, if the deviation range is out of a preset range, updating the candidate index group, and constructing an updated candidate mapping relation based on working parameters of each index in the updated candidate index group at each first preset time and actual measurement reference temperatures of the combustion chamber at each first preset time until the deviation range of the updated candidate mapping relation is in the preset range; taking the candidate mapping relation with the deviation range within the preset range as a target mapping relation of the reference temperature of the combustion chamber; the specific implementation of each step is not repeated here.
It should be noted that examples of the computer readable storage medium may also include, but are not limited to, a phase change memory (PRAM), a Static Random Access Memory (SRAM), a Dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), a Read Only Memory (ROM), an Electrically Erasable Programmable Read Only Memory (EEPROM), a flash memory, or other optical or magnetic storage medium, which will not be described in detail herein.
Exemplary computing device
Having described the method, apparatus, and medium of the exemplary embodiments of the present invention, a computing device 80 of the exemplary embodiments of the present invention is next described with reference to FIG. 4.
FIG. 4 illustrates a block diagram of an exemplary computing device 80 suitable for use in implementing embodiments of the invention, the computing device 80 may be a computer system or a server. The computing device 80 shown in fig. 4 is merely an example and should not be taken as limiting the functionality and scope of use of embodiments of the present invention.
As shown in fig. 4, components of computing device 80 may include, but are not limited to: one or more processors or processing units 801, a system memory 802, and a bus 803 that connects the various system components (including the system memory 802 and processing units 801).
Computing device 80 typically includes a variety of computer system readable media. Such media can be any available media that is accessible by computing device 80 and includes both volatile and nonvolatile media, removable and non-removable media.
The system memory 802 may include computer system readable media in the form of volatile memory, such as Random Access Memory (RAM) 8021 and/or cache memory 8022. Computing device 70 may further include other removable/non-removable, volatile/nonvolatile computer system storage media. By way of example only, ROM8023 may be used to read from and write to non-removable, nonvolatile magnetic media (not shown in FIG. 4, commonly referred to as a "hard disk drive"). Although not shown in fig. 4, a magnetic disk drive for reading from and writing to a removable non-volatile magnetic disk (e.g., a "floppy disk"), and an optical disk drive for reading from or writing to a removable non-volatile optical disk (e.g., a CD-ROM, DVD-ROM, or other optical media), may be provided. In such cases, each drive may be coupled to bus 803 via one or more data medium interfaces. The system memory 802 may include at least one program product having a set (e.g., at least one) of program modules configured to carry out the functions of the embodiments of the invention.
A program/utility 8025 having a set (at least one) of program modules 8024 may be stored, for example, in system memory 802, and such program modules 8024 include, but are not limited to: an operating system, one or more application programs, other program modules, and program data, each or some combination of which may include an implementation of a network environment. Program modules 8024 generally perform the functions and/or methods in the embodiments described herein.
The computing device 80 may also communicate with one or more external devices 804 (e.g., keyboard, pointing device, display, etc.). Such communication may be through an input/output (I/O) interface. Moreover, computing device 80 may also communicate with one or more networks such as a Local Area Network (LAN), a Wide Area Network (WAN) and/or a public network, such as the Internet, through network adapter 806. As shown in fig. 4, network adapter 806 communicates with other modules of computing device 80 (e.g., processing unit 801, etc.) over bus 803. It should be appreciated that although not shown in fig. 4, other hardware and/or software modules may be used in connection with computing device 80.
The processing unit 801 executes various functional applications and data processing by running a program stored in the system memory 802, for example, acquiring measured reference temperatures of the combustion chamber at a plurality of first preset times, respectively; respectively acquiring working parameters of a plurality of strong association indexes and a plurality of weak association indexes of the combustion chamber at a plurality of first preset moments, wherein the plurality of strong association indexes comprise: turbine exhaust temperature, fuel flow, compressor outlet pressure, and atmospheric temperature; the plurality of weakly associated candidate indexes includes: compressor outlet temperature, barometric pressure, air intake flow, and turbine exhaust pressure; obtaining a candidate index group from the strong correlation indexes and the weak correlation indexes, wherein the candidate index group comprises all the strong correlation indexes and at least one weak correlation index; constructing a candidate mapping relation between the reference temperature of the combustion chamber and the candidate index group based on working parameters of each index in the candidate index group at each first preset time and the actually measured reference temperature of the combustion chamber at each first preset time; determining a deviation range of the candidate mapping relation, if the deviation range is out of a preset range, updating the candidate index group, and constructing an updated candidate mapping relation based on working parameters of each index in the updated candidate index group at each first preset time and actual measurement reference temperatures of the combustion chamber at each first preset time until the deviation range of the updated candidate mapping relation is in the preset range; taking the candidate mapping relation with the deviation range within the preset range as a target mapping relation of the reference temperature of the combustion chamber; the specific implementation of each step is not repeated here.
Furthermore, although the operations of the methods of the present application are depicted in the drawings in a particular order, this is not required to either imply that the operations must be performed in that particular order or that all of the illustrated operations be performed to achieve desirable results. Additionally or alternatively, certain steps may be omitted, multiple steps combined into one step to perform, and/or one step decomposed into multiple steps to perform.
While the spirit and principles of the present application have been described with reference to several particular embodiments, it is to be understood that the application is not limited to the disclosed embodiments nor does it imply that features of the various aspects are not useful in combination, nor are they useful in any combination, such as for convenience of description. The application is intended to cover various modifications and equivalent arrangements included within the spirit and scope of the appended claims.
The foregoing description is only of the optional embodiments of the present application, and is not intended to limit the scope of the application, and all the equivalent structural changes made by the description of the present application and the accompanying drawings or the direct/indirect application in other related technical fields are included in the scope of the application.
Through the above description, the embodiments of the present application at least provide the following technical solutions, but are not limited thereto:
1. a method of calculating a reference temperature for a combustion chamber of a gas turbine, comprising:
respectively obtaining actually measured reference temperatures of the combustion chamber at a plurality of first preset moments;
respectively acquiring working parameters of a plurality of strong association indexes and a plurality of weak association indexes of the combustion chamber at a plurality of first preset moments, wherein the plurality of strong association indexes comprise: turbine exhaust temperature, fuel flow, compressor outlet pressure, and atmospheric temperature; the plurality of weakly associated candidate indexes includes: compressor outlet temperature, barometric pressure, air intake flow, and turbine exhaust pressure;
obtaining a candidate index group from the strong correlation indexes and the weak correlation indexes, wherein the candidate index group comprises all the strong correlation indexes and at least one weak correlation index;
constructing a candidate mapping relation between the reference temperature of the combustion chamber and the candidate index group based on working parameters of each index in the candidate index group at each first preset time and the actually measured reference temperature of the combustion chamber at each first preset time;
Determining a deviation range of the candidate mapping relation, if the deviation range is out of a preset range, updating the candidate index group, and constructing an updated candidate mapping relation based on working parameters of each index in the updated candidate index group at each first preset time and actual measurement reference temperatures of the combustion chamber at each first preset time until the deviation range of the updated candidate mapping relation is in the preset range;
and taking the candidate mapping relation with the deviation range within the preset range as a target mapping relation of the reference temperature of the combustion chamber.
2. The gas turbine combustor reference temperature calculation method according to claim 1, wherein the updating the candidate index set includes:
changing the number of weak association indexes in the candidate index group; or,
and changing weak association indexes in the candidate index group.
3. The gas turbine combustor reference temperature calculation method according to claim 1 or 2, wherein the number of the plurality of first preset times is not less than four, and the plurality of first preset times includes a start combustion time of the combustor, a stabilization time after combustion stabilization, and a time between the start combustion time and the stabilization time.
4. The method for calculating a reference temperature of a combustion chamber of a gas turbine according to any one of claims 1 to 3, wherein the constructing a candidate mapping relationship of the reference temperature of the combustion chamber based on the operating parameters of each index in the candidate index set at each first preset time and the measured reference temperature of the combustion chamber at each first preset time includes:
constructing a parameter matrix based on working parameters of each index in the candidate index group at each first preset time;
constructing a temperature matrix based on the actually measured reference temperatures of the combustion chamber at each first preset time;
determining an influence factor of each index in the candidate index group based on the parameter matrix and the temperature matrix;
and constructing a candidate mapping relation of the reference temperature of the combustion chamber based on the influence factors of the indexes in the candidate index group.
5. The gas turbine combustor reference temperature calculation method according to any one of claims 1 to 4, wherein the candidate map is as follows:
CRT＝(K 1 * 1 + 2 * 2 + 3 * 3 + 4 * 4 )+(K 5 * 1 +…+ n-1 *
Yn-5+Kn*；
wherein the CRT represents the reference temperature of the combustion chamber, X 1 To X 4 Work parameters, K, representing the respective strongly associated index 1 To K 4 An influence factor representing each strong correlation index, Y 1 To Y n-5 Work parameters, K, representing the respective weakly associated index 5 To K n-1 Influence factors representing weak correlation indexes, C being a constant, K n Is the influencing factor of C.
6. The gas turbine combustor reference temperature calculation method according to any one of claims 1 to 5, wherein the determining the deviation range of the candidate map includes:
determining verification indexes, wherein the verification indexes are strong association indexes and weak association indexes in candidate index groups corresponding to the candidate mapping relation;
respectively acquiring working parameters of verification indexes of the combustion chamber at a plurality of preset second preset moments;
respectively obtaining actual measurement reference temperatures of the combustion chamber at a plurality of preset second preset moments;
determining a theoretical reference temperature of the combustion chamber at each second preset time based on the candidate mapping relation and working parameters of each verification index at each second preset time;
and determining the deviation range of the candidate mapping relation based on the theoretical reference temperature and the actual measurement reference temperature of the combustion chamber at each second preset time.
7. The gas turbine combustor reference temperature calculation method according to any one of claims 1 to 6, wherein at least one second preset time is provided between two adjacent first preset times.
8. The gas turbine combustor reference temperature calculation method as set forth in claim 1, wherein the preset range is-1.08K to +0.49K.
9. A gas turbine combustor reference temperature calculation apparatus comprising:
the acquisition module is used for respectively acquiring the actually measured reference temperatures of the combustion chamber at a plurality of first preset moments; and
respectively acquiring working parameters of a plurality of strong association indexes and a plurality of weak association indexes of the combustion chamber at a plurality of first preset moments, wherein the plurality of strong association indexes comprise: turbine exhaust temperature, fuel flow, compressor outlet pressure, and atmospheric temperature; the plurality of weakly associated candidate indexes includes: compressor outlet temperature, barometric pressure, air intake flow, and turbine exhaust pressure;
the processing module is used for acquiring candidate index groups from the strong correlation indexes and the weak correlation indexes, wherein the candidate index groups comprise all the strong correlation indexes and at least one weak correlation index;
constructing a candidate mapping relation between the reference temperature of the combustion chamber and the candidate index group based on working parameters of each index in the candidate index group at each first preset time and the actually measured reference temperature of the combustion chamber at each first preset time;
Determining a deviation range of the candidate mapping relation, if the deviation range is out of a preset range, updating the candidate index group, and constructing an updated candidate mapping relation based on working parameters of each index in the updated candidate index group at each first preset time and actual measurement reference temperatures of the combustion chamber at each first preset time until the deviation range of the updated candidate mapping relation is in the preset range;
and taking the candidate mapping relation with the deviation range within the preset range as a target mapping relation of the reference temperature of the combustion chamber.
10. The gas turbine combustor reference temperature calculation apparatus of claim 9, the processing module configured to update the candidate index set by:
changing the number of weak association indexes in the candidate index group; or,
and changing weak association indexes in the candidate index group.
11. The gas turbine combustor reference temperature calculation apparatus according to claim 9 or 10, wherein the number of the plurality of first preset times is not less than four, and the plurality of first preset times includes a start combustion time of the combustor, a stabilization time after combustion stabilization, and a time between the start combustion time and the stabilization time.
12. The gas turbine combustor reference temperature calculation apparatus of any one of claims 9-11, the processing module further configured to:
constructing a parameter matrix based on working parameters of each index in the candidate index group at each first preset time;
constructing a temperature matrix based on the actually measured reference temperatures of the combustion chamber at each first preset time;
determining an influence factor of each index in the candidate index group based on the parameter matrix and the temperature matrix;
and constructing a candidate mapping relation of the reference temperature of the combustion chamber based on the influence factors of the indexes in the candidate index group.
13. The gas turbine combustor reference temperature calculation apparatus as set forth in any one of claims 9 to 12, the candidate map is constructed by:
CRT＝(K 1 * 1 + 2 * 2 + 3 * 3 + 4 * 4 )+(K 5 * 1 +…+ n-1 *
Yn-5+Kn*；
wherein the CRT represents the reference temperature of the combustion chamber, X 1 To X 4 Work parameters, K, representing the respective strongly associated index 1 To K 4 An influence factor representing each strong correlation index, Y 1 To Y n-5 Work parameters, K, representing the respective weakly associated index 5 To K n-1 Influence factors representing weak correlation indexes, C being a constant, K n Is the influencing factor of C.
14. The gas turbine combustor reference temperature calculation apparatus of any one of claims 9-13, the processing module further configured to:
Determining verification indexes, wherein the verification indexes are strong association indexes and weak association indexes in candidate index groups corresponding to the candidate mapping relation;
respectively acquiring working parameters of verification indexes of the combustion chamber at a plurality of preset second preset moments;
respectively obtaining actual measurement reference temperatures of the combustion chamber at a plurality of preset second preset moments;
determining a theoretical reference temperature of the combustion chamber at each second preset time based on the candidate mapping relation and working parameters of each verification index at each second preset time;
and determining the deviation range of the candidate mapping relation based on the theoretical reference temperature and the actual measurement reference temperature of the combustion chamber at each second preset time.
15. The gas turbine combustor reference temperature calculation apparatus as set forth in any one of claims 9 to 14, wherein at least one of the second preset times is provided between two adjacent first preset times.
16. The gas turbine combustor reference temperature calculation apparatus as set forth in any one of claims 9 to 15, wherein the preset range is-1.08K to +0.49K.
17. A computer readable storage medium, on which a computer program is stored which, when executed by a processor, implements the method of any of claims 1-8.
18. A computing device comprising a processor for implementing the method of any of claims 1-8 when executing a computer program stored in memory.

Claims (10)

1. A method of calculating a reference temperature for a combustion chamber of a gas turbine, comprising:
respectively obtaining actually measured reference temperatures of the combustion chamber at a plurality of first preset moments;
respectively acquiring working parameters of a plurality of strong association indexes and a plurality of weak association indexes of the combustion chamber at a plurality of first preset moments, wherein the plurality of strong association indexes comprise: turbine exhaust temperature, fuel flow, compressor outlet pressure, and atmospheric temperature; the plurality of weakly associated candidate indexes includes: compressor outlet temperature, barometric pressure, air intake flow, and turbine exhaust pressure;
obtaining a candidate index group from the strong correlation indexes and the weak correlation indexes, wherein the candidate index group comprises all the strong correlation indexes and at least one weak correlation index;
constructing a candidate mapping relation between the reference temperature of the combustion chamber and the candidate index group based on working parameters of each index in the candidate index group at each first preset time and the actually measured reference temperature of the combustion chamber at each first preset time;
Determining a deviation range of the candidate mapping relation, if the deviation range is out of a preset range, updating the candidate index group, and constructing an updated candidate mapping relation based on working parameters of each index in the updated candidate index group at each first preset time and actual measurement reference temperatures of the combustion chamber at each first preset time until the deviation range of the updated candidate mapping relation is in the preset range;
and taking the candidate mapping relation with the deviation range within the preset range as a target mapping relation of the reference temperature of the combustion chamber.
2. The gas turbine combustor reference temperature calculation method of claim 1, the updating the candidate set of metrics comprising:
changing the number of weak association indexes in the candidate index group; or,
and changing weak association indexes in the candidate index group.
3. The gas turbine combustor reference temperature calculation method as set forth in claim 1, wherein the number of the plurality of first preset times is not less than four, and the plurality of first preset times includes a start combustion time of the combustor, a stabilization time after combustion stabilization, and a time between the start combustion time and the stabilization time.
4. The gas turbine combustor reference temperature calculation method as set forth in claim 1, said constructing a candidate map of the combustor reference temperature based on the operating parameters of each index in the candidate index set at each first preset time and the measured reference temperature of the combustor at each first preset time, comprising:
constructing a parameter matrix based on working parameters of each index in the candidate index group at each first preset time;
constructing a temperature matrix based on the actually measured reference temperatures of the combustion chamber at each first preset time;
determining an influence factor of each index in the candidate index group based on the parameter matrix and the temperature matrix;
and constructing a candidate mapping relation of the reference temperature of the combustion chamber based on the influence factors of the indexes in the candidate index group.
5. A gas turbine combustor reference temperature calculation apparatus comprising:
the acquisition module is used for respectively acquiring the actually measured reference temperatures of the combustion chamber at a plurality of first preset moments; and
respectively acquiring working parameters of a plurality of strong association indexes and a plurality of weak association indexes of the combustion chamber at a plurality of first preset moments, wherein the plurality of strong association indexes comprise: turbine exhaust temperature, fuel flow, compressor outlet pressure, and atmospheric temperature; the plurality of weakly associated candidate indexes includes: compressor outlet temperature, barometric pressure, air intake flow, and turbine exhaust pressure;
The processing module is used for acquiring candidate index groups from the strong correlation indexes and the weak correlation indexes, wherein the candidate index groups comprise all the strong correlation indexes and at least one weak correlation index;
constructing a candidate mapping relation between the reference temperature of the combustion chamber and the candidate index group based on working parameters of each index in the candidate index group at each first preset time and the actually measured reference temperature of the combustion chamber at each first preset time;
determining a deviation range of the candidate mapping relation, if the deviation range is out of a preset range, updating the candidate index group, and constructing an updated candidate mapping relation based on working parameters of each index in the updated candidate index group at each first preset time and actual measurement reference temperatures of the combustion chamber at each first preset time until the deviation range of the updated candidate mapping relation is in the preset range;
and taking the candidate mapping relation with the deviation range within the preset range as a target mapping relation of the reference temperature of the combustion chamber.
6. The gas turbine combustor reference temperature calculation apparatus of claim 5, the processing module configured to update the candidate set of metrics by:
Changing the number of weak association indexes in the candidate index group; or,
and changing weak association indexes in the candidate index group.
7. The gas turbine combustor reference temperature calculation apparatus as set forth in claim 5, wherein the number of the plurality of first preset times is not less than four, and the plurality of first preset times includes a start combustion time of the combustor, a stabilization time after combustion stabilization, and a time between the start combustion time and the stabilization time.
8. The gas turbine combustor reference temperature calculation apparatus of claim 5, the processing module further configured to:
constructing a parameter matrix based on working parameters of each index in the candidate index group at each first preset time;
constructing a temperature matrix based on the actually measured reference temperatures of the combustion chamber at each first preset time;
determining an influence factor of each index in the candidate index group based on the parameter matrix and the temperature matrix;
and constructing a candidate mapping relation of the reference temperature of the combustion chamber based on the influence factors of the indexes in the candidate index group.
9. A computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements the method of any of claims 1-4.
10. A computing device comprising a processor for implementing the method of any of claims 1-4 when executing a computer program stored in memory.
CN202311164585.6A 2023-09-11 2023-09-11 Method, device, medium and equipment for calculating reference temperature of combustion chamber of gas turbine Active CN117153282B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311164585.6A CN117153282B (en) 2023-09-11 2023-09-11 Method, device, medium and equipment for calculating reference temperature of combustion chamber of gas turbine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311164585.6A CN117153282B (en) 2023-09-11 2023-09-11 Method, device, medium and equipment for calculating reference temperature of combustion chamber of gas turbine

Publications (2)

Publication Number Publication Date
CN117153282A true CN117153282A (en) 2023-12-01
CN117153282B CN117153282B (en) 2024-04-16

Family

ID=88900556

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311164585.6A Active CN117153282B (en) 2023-09-11 2023-09-11 Method, device, medium and equipment for calculating reference temperature of combustion chamber of gas turbine

Country Status (1)

Country Link
CN (1) CN117153282B (en)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ITCO20090054A1 (en) * 2009-11-27 2011-05-28 Nuovo Pignone Spa METHOD OF MODE CONTROL BASED ON EXHAUST TEMPERATURE FOR GAS TURBINE AND GAS TURBINE
US20160069271A1 (en) * 2014-09-08 2016-03-10 General Electric Company Bulk Flame Temperature Regulator for Dry Low Emission Engines
CN109581870A (en) * 2018-11-27 2019-04-05 中国工程物理研究院化工材料研究所 The temperature in the kettle dynamic matrix control method of energetic material reaction kettle
CN109917829A (en) * 2019-03-20 2019-06-21 华中科技大学 A kind of wind-tunnel temperature gradient control device and its control method
CN111144721A (en) * 2019-12-12 2020-05-12 国网山东省电力公司经济技术研究院 Power grid project demand evaluation model construction method and device and computing equipment
CN113091928A (en) * 2021-04-01 2021-07-09 广东电网有限责任公司佛山供电局 High-voltage chamber equipment temperature fault monitoring method and related device
CN113468762A (en) * 2021-07-22 2021-10-01 广东电网有限责任公司广州供电局 Hot spot temperature calculation method and device, computer equipment and storage medium
CN113851757A (en) * 2021-09-24 2021-12-28 经纬恒润(天津)研究开发有限公司 Power battery thermal management method and device
US20220065127A1 (en) * 2020-09-01 2022-03-03 Purdue Research Foundation Method for reconstructing non-uniform circumferential flow in gas turbine engines
CN114993689A (en) * 2022-04-26 2022-09-02 华电电力科学研究院有限公司 Method, device, equipment and medium for evaluating combustion state of gas turbine
CN115906467A (en) * 2022-11-16 2023-04-04 协鑫电港云科技(海南)有限公司 Data processing method and device based on battery swapping station, electronic equipment and storage medium

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ITCO20090054A1 (en) * 2009-11-27 2011-05-28 Nuovo Pignone Spa METHOD OF MODE CONTROL BASED ON EXHAUST TEMPERATURE FOR GAS TURBINE AND GAS TURBINE
US20160069271A1 (en) * 2014-09-08 2016-03-10 General Electric Company Bulk Flame Temperature Regulator for Dry Low Emission Engines
CN109581870A (en) * 2018-11-27 2019-04-05 中国工程物理研究院化工材料研究所 The temperature in the kettle dynamic matrix control method of energetic material reaction kettle
CN109917829A (en) * 2019-03-20 2019-06-21 华中科技大学 A kind of wind-tunnel temperature gradient control device and its control method
CN111144721A (en) * 2019-12-12 2020-05-12 国网山东省电力公司经济技术研究院 Power grid project demand evaluation model construction method and device and computing equipment
US20220065127A1 (en) * 2020-09-01 2022-03-03 Purdue Research Foundation Method for reconstructing non-uniform circumferential flow in gas turbine engines
CN113091928A (en) * 2021-04-01 2021-07-09 广东电网有限责任公司佛山供电局 High-voltage chamber equipment temperature fault monitoring method and related device
CN113468762A (en) * 2021-07-22 2021-10-01 广东电网有限责任公司广州供电局 Hot spot temperature calculation method and device, computer equipment and storage medium
CN113851757A (en) * 2021-09-24 2021-12-28 经纬恒润(天津)研究开发有限公司 Power battery thermal management method and device
CN114993689A (en) * 2022-04-26 2022-09-02 华电电力科学研究院有限公司 Method, device, equipment and medium for evaluating combustion state of gas turbine
CN115906467A (en) * 2022-11-16 2023-04-04 协鑫电港云科技(海南)有限公司 Data processing method and device based on battery swapping station, electronic equipment and storage medium

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
YANGYI LIU等: "Noncontact temperature estimation method on the actively cooled primary mirror surface of large ground-based solar telescope", JOURNAL OF ASTRONOMICAL TELESCOPES, INSTRUMENTS, AND SYSTEMS, 2 November 2017 (2017-11-02), pages 1 - 8 *
何皑;蒋洪德;谢法;: "燃气轮机燃烧基准温度估算方法", 燃气轮机技术, no. 04, 16 December 2014 (2014-12-16) *
吴彬, 陈海耿: "均匀介质参与体系辐射换热的直接矩阵解法", 东北大学学报(自然科学版), no. 08, 30 August 2002 (2002-08-30) *
田丰等: "基于RBF神经网络的温度场重建算法研究", 仪器仪表学报, 30 November 2006 (2006-11-30), pages 1460 - 1464 *

Also Published As

Publication number Publication date
CN117153282B (en) 2024-04-16

Similar Documents

Publication Publication Date Title
Kwiatkowski et al. Automated generation and assessment of affine LPV models
US7848910B2 (en) Method for analyzing the behavior of complex systems, especially internal combustion engines
Gutjahr et al. Advanced modeling and optimization for virtual calibration of internal combustion engines
Vogt et al. On-line adaptation of grid-based look-up tables using a fast linear regression technique
CN110671219B (en) Method and system for controlling air-fuel ratio of gas engine
JP2006017698A (en) Setting method of experimental planning for engine control parameter, program for making computer execute the setting method of experimental planning, computer-readable recording medium storing the program, and setting apparatus of experimental planning for engine control parameter
CN117153282B (en) Method, device, medium and equipment for calculating reference temperature of combustion chamber of gas turbine
JP2011021518A (en) Virtual adaptation system of engine
Pochet et al. Uncertainty quantification from raw measurements to post-processed data: A general methodology and its application to a homogeneous-charge compression–ignition engine
CN112489734B (en) Simplified method of combustion reaction mechanism model of alternative fuel dimethyl ether of internal combustion engine
JP4272234B2 (en) Method and apparatus for estimating gas concentration in coal fired boiler
CN116258082A (en) Working medium thermodynamic property calculation method and system based on RBF neural network
CN114692525B (en) Combustion simulation dimension reduction and speed acceleration method and device and steady state calculation method
CN114444401A (en) Chemical reaction mechanism optimizing method and device
CN114861126A (en) Method and system for predicting stability in combustion adjustment process of gas turbine
CN113901733A (en) Method for determining thermal environment numerical value of gas flow field
Egorov et al. The methodology of stochastic optimization of parameters and control laws for the aircraft gas-turbine engines flow passage components
Jung et al. In-cylinder pressure-based low-pressure-cooled exhaust gas recirculation estimation methods for turbocharged gasoline direct injection engines
Mai et al. Impact of measurement uncertainty in the characteristic maps of a turbocharger on engine performance
CN111399210A (en) Coarse alignment method, device, equipment and storage medium for large-caliber large-field telescope
Bensouda Extending and formalizing the energy signature method for calibrating simulations and illustrating with application for three California climates
CN110069856A (en) CFD system and implementation method based on adaptive reaction mechanism
CN116522683A (en) Building energy consumption carbon emission determination method
CN113722998B (en) Temperature calibration method, computer and storage medium for passive temperature zone of multispectral pyrometer
CN115680780B (en) Turbine blade plane blade grid inlet axial speed control method

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