Gas turbine control device, gas turbine control method, and program

Abstract

A flow rate per unit time of fuel fed to a gas turbine is calculated. A flow rate per unit time of air fed to the gas turbine is calculated. A turbine inlet temperature is calculated through use of a physical model formula expressing a relationship of input and output of thermal energy relating to a combustor of the gas turbine. A fuel distribution ratio for each of a plurality of fuel supply systems connected to the combustor is calculated based on the turbine inlet temperature.

Claims

1. A gas turbine control device for controlling a gas turbine including a compressor, a combustor, and a turbine, the gas turbine control device comprising a non-transitory memory and a processor configured to execute instructions stored in the memory, the memory and the processor being further configured to: calculate a fuel flow rate per unit time of a fuel fed to the combustor through a plurality of fuel supply systems according to a first formula including (i) a manifold pressure of the combustor, (ii) a pressure of a casing of the combustor, and (iii) a fuel temperature; calculate an air flow rate per unit time of an air flowing into the compressor according to a second formula including an index pressure difference and an inlet temperature of the compressor; calculate a turbine inlet temperature according to a third formula including the fuel flow rate per unit time, the air flow rate per unit time, the fuel temperature, and an air temperature at the casing of the combustor, wherein the turbine inlet temperature reflects a specific enthalpy changed in response to an air humidity; correct the turbine inlet temperature responsive to an exhaust gas temperature at the turbine to obtain a corrected turbine inlet temperature; obtain a plurality of fuel distribution ratios for delivering fuel to the combustor through a plurality of fuel supply systems in response to obtaining the corrected turbine inlet temperature, and individually control a plurality of regulating valves attached to the plurality of fuel supply systems in response to the obtaining the plurality of fuel distribution ratios.

2. The gas turbine control device according to claim 1, wherein the memory and the processor are further configured to correct the turbine inlet temperature is by multiplying a correction factor according to a ratio of the exhaust gas temperature at the turbine and an exhaust gas temperature estimation value estimated from the turbine inlet temperature.

3. The gas turbine control device according to claim 1, wherein the memory and the processor are further configured to calculate the turbine inlet temperature according to the third formula using the air flow rate per unit time to be reduced by an extracted amount of air flow rate per unit time, the extracted amount of air flow being extracted before air flows into the combustor following the compressor.

4. The gas turbine control device according to claim 1, wherein the memory and the processor are further configured to calculate the air temperature at the casing of the combustor according to a fourth formula including an inlet pressure of the compressor and an outlet pressure of the compressor as well as the inlet temperature of the compressor and a measured air temperature of the casing of the combustor.

5. A gas turbine control method to be implemented by a gas turbine control device configured to control a gas turbine including a compressor, a combustor, and a turbine, the gas turbine control device including a non-transitory memory and a processor configured to execute instructions stored in the non-transitory memory, the non-transitory memory and the processor being further configured to perform the gas turbine control method comprising: calculating a fuel flow rate per unit time of a fuel fed to the combustor through a plurality of fuel supply systems according to a first formula including (i) a manifold pressure of the combustor, (ii) a pressure of a casing of the combustor, and (iii) a fuel temperature; calculating an air flow rate per unit time of an air flowing into the compressor according to a second formula including an index pressure difference, and an inlet temperature of the compressor; calculating a turbine inlet temperature according to a third formula including the fuel flow rate per unit time, the air flow rate per unit time, the fuel temperature, and an air temperature at the casing of the combustor, wherein the turbine inlet temperature reflects a specific enthalpy changed in response to an air humidity; correcting the turbine inlet temperature responsive to an exhaust gas temperature at the turbine to obtain a corrected turbine inlet temperature; obtaining a plurality of fuel distribution ratios for delivering fuel to the combustor through a plurality of fuel supply systems in response to obtaining the corrected turbine inlet temperature; and individually controlling a plurality of regulating valves attached to the plurality of fuel supply systems in response to the obtaining the plurality of the fuel distribution ratios.

6. A non-transitory computer-readable storage medium including a stored program causing a computer to implement the gas turbine control method of claim 5.

Description

BRIEF DESCRIPTION OF DRAWINGS

(1) FIG. 1 is a system diagram of a gas turbine power plant.

(2) FIG. 2 is a function block diagram of a gas turbine control device according to a first embodiment.

(3) FIG. 3 is a process flowchart of the gas turbine control device.

(4) FIG. 4 is u function block diagram of a gas turbine control device according to a second embodiment.

DESCRIPTION OF EMBODIMENTS

First Embodiment

(5) In the following, description is made of a gas turbine control device according to a first embodiment of the present invention and a gas turbine power plant including the device with reference to the drawings.

(6) FIG. 1 is a system diagram of the gas turbine power plant according to the present embodiment.

(7) As illustrated in FIG. 1, a gas turbine power plant 100 in the present embodiment includes a gas turbine 10, a generator 16 that generates electricity by being driven by the gas turbine 10, a gas turbine control device 20 that controls the gas turbine 10, and a supply device 40 that supplies fuel. The gas turbine 10 and the generator 16 are connected by a rotor 15. Note that the gas turbine power plant 100 may be a gas turbine combined cycle power plant (GTCC) obtained by further including a steam turbine and the rotor 15 connected to the steam turbine.

(8) The gas turbine 10 includes a compressor 11, a fuel supply system 50, a combustor 12, and a turbine 13.

(9) The compressor 11 generates high-pressure compressed air by compressing external air. More specifically, the compressor 11 includes a compressor rotor 111 that rotates about a main axis line and a compressor casing 112 that covers the compressor rotor 111 from an outer circumferential side. The turbine 13 includes a turbine rotor 131 that rotates about the main axis line and a turbine casing 132 that covers the turbine rotor 131 from an outer circumferential side.

(10) The compressor rotor 111 and the turbine rotor 131 are integrally connected to each other on the main axis line.

(11) The combustor 12 generates high-temperature and high-pressure combustion gas by combusting fuel supplied from the fuel supply system 50 in the compressed air generated by the compressor 11. The combustor 12 includes a combustor internal cylinder 121 and a combustor casing 17 that covers the combustor internal cylinder 121. Further, the compressor casing 112, the combustor casing 17, and the turbine casing 132 are connected to each other.

(12) The generator 16 is connected to one end of the compressor rotor 111. The generator 16 is driven by rotation of the compressor rotor 111, and thus generates power.

(13) The fuel supply system 50 regulates a pressure and a flow rate of fuel supplied from the supply device 40. The fuel supply system 50 includes a regulating valve 18, a nozzle 123, and a manifold pipe 124. The regulating valve 18 is a pressure regulating valve, a fuel flow rate regulating valve, or the like.

(14) Note that FIG. 1 illustrates a mode in which the gas turbine 10 includes fuel supply systems 50a, 50b, 50c, . . . . Further, FIG. 1 illustrates a mode in which, as the regulating valve 18, a first regulating valve 18a, a second regulating valve 18b, and a third regulating valve 18e, which regulate a fuel flow rate and a pressure for the fuel supply systems, respectively, are provided. FIG. 1 illustrates only three systems, including the fuel supply systems 50a, 50b, and 50c. However, the number of the fuel supply systems is not limited to three. For example, the gas turbine 10 may include a top hat fuel supply system that supplies fuel to a top hat portion of the combustor 12, a pilot fuel supply system that supplies fuel to a center portion of the internal cylinder of the combustor 12, and a main fuel supply system that supplies fuel to a portion surrounding the pilot fuel supply system in the internal cylinder of the combustor 12. Further. FIG. 1 illustrates only one combustor 12. However, a plurality of, for example, sixteen combustors 12 are provided, and fuel is supplied to the combustors 12 from the corresponding fuel supply systems 50a, 50b, 50c, . . . . As one example, in the gas turbine power plant 100, sixteen fuel supply systems 50 are connected to the corresponding combustors 12, respectively.

(15) The gas turbine 10 is provided with a plurality of measuring instruments that measure state quantities of the gas turbine 10. As one example, the gas turbine 10 is provided with a fuel flow rate sensor 20a that measures a fuel flow rate at the fuel supply system 50, a nozzle flow rate sensor 20b that measures a passage flow rate (nozzle passage flow rate) at the nozzle 123, an exhaust gas flow rate sensor 20c that measures an exhaust gas flow rate, a manifold pressure sensor 20d that measures the pressure P.sub.in in the manifold pipe 124, a temperature sensor 20c that measures a compressor inlet temperature, an index pressure difference gauge 20f that measures a compressor index pressure difference P.sub.index, and the like. Note that the compressor index pressure difference P.sub.index is a pressure difference between a pressure at a compressor suction port easing portion of the compressor 11 and a pressure in the vicinity of a blade inside the compressor, and is a value being an index of a flow rate of air sucked by the compressor. All of the fuel flow rate, the nozzle passage flow rate or the fuel valve passage flow rate, the exhaust gas flow rate, and the manifold pressure are state quantities of the gas turbine 10. Further, the gas turbine 10 is also provided with other measuring instruments. For example, the gas turbine 10 is provided with measuring instruments that measure a pressure P.sub.out at the combustor casing 17, a fuel temperature T.sub.f, a casing air temperature T.sub.cs, and an exhaust gas temperature T.sub.2T. These are state quantities of the gas turbine 10.

(16) Note that a turbine inlet temperature described below indicates a temperature at an inlet of the turbine 13 of combustion gas which has a high-temperature and is jetted from the combustor 12. The compressor 11 is provided with an inlet guide vane (IGV) 14. The IGV 14 regulates an amount of air flowing in the compressor 11.

(17) FIG. 2 is a function block diagram of the gas turbine control device according to the first embodiment.

(18) The gas turbine control device 20 is a computer, and may have a hardware configuration including a storage unit such as a read-only memory (ROM), a random access memory (RAM), and a hard disk drive (HDD), a central processing unit (CPU), a communication interface, and the like.

(19) Based on a user operation, the CPU of the gas turbine control device 20 executes a stored control program. With this, the gas turbine control device 20 includes a fuel flow rate calculation unit 201, an air flow rate calculation unit 202, a turbine inlet temperature calculation unit 203, a correction unit 204, a fuel distribution ratio calculation unit 205, and a valve opening calculation unit 206. Note that the fuel distribution ratio calculation unit 205 includes a first fuel distribution ratio calculation unit 205a, a second fuel distribution ratio calculation unit 205b, a third fuel distribution ratio calculation unit 205c, and so on that calculate a fuel distribution ratio for the corresponding fuel supply systems. Further, the valve opening calculation unit 206 includes a first valve opening calculation unit 206a, a second valve opening calculation unit 206b, a third valve opening calculation unit 206c, and so on that calculate an opening amount for the corresponding fuel supply systems.

(20) The fuel flow rate calculation unit 201 calculates a flow rate per unit time of fuel fed to the gas turbine 10, based on a measurement value of the gas turbine 10.

(21) The air flow rate calculation unit 202 calculates a flow rate per unit time of air fed to the gas turbine 10, based on a measurement value of the gas turbine 10.

(22) The turbine inlet temperature calculation unit 203 calculates a turbine inlet temperature through use of a physical model formula expressing a relationship of input and output of thermal energy relating to the combustor 12 of the gas turbine 10. In the physical model formula, a flow rate per unit time of the fuel, a flow rate per unit time of the air, a fuel temperature, and an air temperature at the easing of the gas turbine 10 are used.

(23) The correction unit 204 calculates a correction factor for correcting the turbine inlet temperature, based on a ratio of the turbine inlet temperature and the acquired exhaust gas temperature at the gas turbine 10, and corrects the turbine inlet temperature with the correction factor.

(24) Based on the corrected turbine inlet temperature, the fuel distribution ratio calculation unit 205 (205a, 205b, 205c, . . . ) calculates a fuel distribution ratio for each of the fuel supply systems that are correspondingly connected to the plurality of combustors 12.

(25) Based on the fuel distribution ratio for each of the fuel supply systems, the valve opening calculation unit 206 (206a, 206b, 206c, . . . ) calculates an opening amount for each of the fuel supply systems.

(26) In actuality, by executing the control program, the gas turbine control device 20 has functions in addition to the functional units illustrated in FIG. 2. However, for convenience of explanation, description is made only on the functional units illustrated in FIG. 2.

(27) FIG. 3 is a process flowchart of the gas turbine control device according to the present embodiment.

(28) First, the fuel flow rate calculation unit 201 acquires the pressure P.sub.in at the manifold pipe 124, the pressure P.sub.out at the combustor casing 17, and the fuel temperature T.sub.f. The pressure P.sub.in at the manifold pipe 124 is measured by the manifold pressure sensor 20d. Similarly, an outlet pressure P.sub.out and a fuel temperature T.sub.f at the fuel nozzle are measured by sensors. The fuel flow rate calculation unit 201 substitutes the pressure P.sub.in at the manifold pipe 124, the pressure P.sub.out at the combustor casing 17, and the fuel temperature T.sub.f in a fuel flow rate calculation formula (1), and calculates a total flow rate G.sub.f per unit time of the fuel fed to each of the combustors 12 of the gas turbine 10 as a first function of the above-noted parameters, as indicated below in Formula 1 (Step S101).
[Mathematical Formula 1]
G.sub.f=f(P.sub.in,P.sub.out,T.sub.f)  (1)

(29) Further, the air flow rate calculation unit 202 acquires the compressor index pressure difference P.sub.index from the index pressure difference gauge 20f, and acquires a compressor inlet temperature T.sub.1c from the temperature sensor 20e. By measuring the index pressure difference P.sub.index, a flow rate of air flowing through the compressor 11 can be calculated. The air flow rate calculation unit 202 substitutes the index pressure difference P.sub.index and the compressor inlet temperature T.sub.1c in an air flow rate calculation formula (2), and calculates a flow rate Ga per unit time of the air flowing in the compressor 11 as a second function of the above-noted parameters, as indicated below in Formula 2 (Step S102).
[Mathematical Formula 2]
G.sub.a=f(P.sub.index,T.sub.1C)  (2)

(30) The turbine inlet temperature calculation unit 203 acquires the fuel flow rate G.sub.f calculated by the fuel flow rate calculation unit 201. The turbine inlet temperature calculation unit 203 acquires the air flow rate G.sub.a calculated by the air flow rate calculation unit 202. The turbine inlet temperature calculation unit 203 further acquires the fuel temperature T.sub.f. The turbine inlet temperature calculation unit 203 further acquires the casing air temperature T.sub.cs at the combustor casing 17. The air temperature T.sub.cs at the combustor casing is measured by a sensor. The turbine inlet temperature calculation unit 203 substitutes the fuel flow rate G.sub.f, the air flow rate G.sub.a, the fuel temperature T.sub.f, and the casing air temperature T.sub.cs at the combustor casing 17 that are acquired in an abnormality physical model in the periphery of the combustor, which is indicated with a formula (3), and calculates a turbine inlet temperature T.sub.1T as a third function of the above-noted parameters, as indicated below in Formula 3 (Step S103). The turbine inlet temperature calculation unit 203 outputs the calculated turbine inlet temperature T.sub.1T to the correction unit 204. Note that the abnormality physical model expressed with formula (3) is obtained by rearranging a model formula, which indicates that thermal energy flowing in the combustor 12 is equivalent to thermal energy flowing out from the combustor 12, in such manner that the turbine inlet temperature T.sub.1T included in the model formula is put on the left side and other elements are put on the right side. The thermal energy flowing in the combustor 12 in the abnormality physical model is indicated with a total of thermal energy of fuel, thermal energy of air, and heat generation energy of combustion gas. Further, the thermal energy flowing out from the combustor 12 is indicated with the thermal energy at the inlet of the turbine 13. Here, the thermal energy of air is varied in accordance with a specific enthalpy in consideration with influence of humidity of steam contained in the air flowing in the combustor casing 17. Therefore, the turbine inlet temperature calculation unit 203 may acquire humidity of the air flowing in the combustor casing 17 with a sensor, calculate a specific enthalpy based on the humidity, and calculate the thermal energy of the air through use of the specific enthalpy.
[Mathematical Formula 3]
T.sub.1T=f(G.sub.f,G.sub.a,T.sub.f,T.sub.CS)  (3)

(31) Here, in actuality, the air flowing in the compressor 11 is extracted at a stage before flowing in the combustor 12. Thus, the air flow rate calculation unit 202 may calculate the air flow rate G.sub.a per unit time, obtained by subtracting a flow rate per unit time of the extracted air, and the turbine inlet temperature calculation unit 203 may calculate the turbine inlet temperature T.sub.1T with the above-described formula (3) using the air flow rate G.sub.a obtained by subtracting the extracted air flow rate. The extracted air flow rate is monotonically increased in accordance with increase in output value of the gas turbine 10 and value of the turbine inlet temperature T.sub.1T. The air flow rate calculation unit 202 may calculate the extracted air flow rate through use of interpolation calculation with a data table indicating extracted air flow rates corresponding to output values of the gas turbine 10 and values of the turbine inlet temperature T.sub.1T and a predetermined calculation formula. Note that, when the extracted air flow rate is calculated based on the value of the turbine inlet temperature T.sub.1T, the air flow rate calculation unit 202 may use the value of the turbine inlet temperature T.sub.1T that is calculated previously.

(32) Here, immediately after the gas turbine 10 is actuated, the extracted air flow rate described above is also varied based on change of a shape of the compressor 11, the pipe, and the like due to thermal expansion and the like. Therefore, the air flow rate calculation unit 202 may calculate a highly accurate extracted air flow rate after activation through use of a correction formula for an extracted air flow rate per unit time in accordance with an elapsed time after activation of the gas turbine 10. Further, the turbine inlet temperature calculation unit 203 may calculate the turbine inlet temperature T.sub.1T through use of the air flow rate G.sub.a obtained by subtracting the highly accurate extracted air flow rate thus corrected.

(33) The correction unit 204 may further correct the turbine inlet temperature T.sub.1T, based on the exhaust gas temperature T.sub.2T. Specifically, in advance, the gas turbine control device 20 stores information indicating a relationship of a rated exhaust gas temperature T4 at the time of rated operation of the gas turbine 10 and a rated turbine inlet temperature T5 at that time. The correction unit 204 acquires the information on the relationship of the rated exhaust gas temperature T4 and the rated turbine inlet temperature T5. The correction unit 204 calculates a correction value for the turbine inlet temperature T.sub.1T, based on the relationship between the rated exhaust gas temperature T4 and the rated turbine inlet temperature T5. As one example, the correction unit 204 calculates a correction value for the turbine inlet temperature T.sub.1T. The correction value satisfies that a value Tx obtained by subtracting the exhaust gas temperature T.sub.2T from the rated exhaust gas temperature T4 is equal to a value Ty obtained by subtracting the turbine inlet temperature T.sub.1T from the rated turbine inlet temperature T5 or is equal to a predetermined constant multiplication of a. Note that the relationship of the rated exhaust gas temperature T4 and the rated turbine inlet temperature T5 at that time may be varied due to a pressure at the combustor casing 17. In this case, the correction unit 204 acquires a pressure at the combustor casing 17 with a sensor, and uses the relationship of the rated exhaust gas temperature T4 and the rated turbine inlet temperature T5 at that time that are corrected based on the value of the pressure. In this manner, the correction value for the turbine inlet temperature T.sub.1T may be calculated as described above.

(34) Here, the exhaust gas temperature T.sub.2T is a measurement value with slow reactivity with respect to change in output of the gas turbine 10. However, although having slow reactivity, the exhaust gas temperature T.sub.2T is highly accurate value. The exhaust gas temperature and the turbine inlet temperature have a strong correlation. Therefore, the correction unit 204 corrects the turbine inlet temperature T.sub.1T having fast reactivity that is obtained in Step S103 through use of the highly accurate exhaust gas temperature T.sub.2T. Specifically, the correction unit 204 calculates an exhaust gas temperature estimation value T.sub.4T, based on the turbine inlet temperature T.sub.1T. The exhaust gas temperature estimation value T.sub.4T is an estimation value of an exhaust gas temperature in a case of assuming a temperature at the inlet of the turbine 13 as the turbine inlet temperature T.sub.1T.

(35) The correction unit 204 acquires the turbine inlet temperature T.sub.1T, and inputs the turbine inlet temperature T.sub.1T to a predetermined calculation formula set in advance, and calculates an assumed exhaust gas temperature estimation value T.sub.3T. This calculation formula is a formula for converting the turbine inlet temperature T.sub.1T into an exhaust gas temperature.

(36) The correction unit 204 subjects the assumed exhaust gas temperature estimation value T.sub.3T to a process of delaying change of a value for each time, and calculates the exhaust gas temperature estimation value T.sub.4T.

(37) The correction unit 204 calculates a correction factor X4 that corrects the turbine inlet temperature T.sub.1T, based on a ratio of the exhaust gas temperature T.sub.2T and the exhaust gas temperature estimation value T.sub.4T calculated from the turbine inlet temperature T.sub.1T. In the present embodiment, as one example, the correction factor X4 is a value obtained by summing an adjustment ratio X2 and a past adjustment ratio X3. The adjustment ratio X2 is obtained by multiplying an adjustment coefficient α with a ratio X2 obtained by dividing the exhaust gas temperature T.sub.2T by the exhaust gas temperature estimation value T.sub.4T. The past adjustment ratio X3 is obtained by multiplying 1−α with a past correction factor X4′ being the correction factor X4 that is previously calculated. Further, the correction unit 204 calculates a post correction turbine inlet temperature T.sub.1T′ by multiplying the turbine inlet temperature T.sub.1T with the correction factor X4 (Step S104). The correction factor X4 is a ratio of the exhaust gas temperature T.sub.2T to the exhaust gas temperature estimation value T.sub.4T, and hence the post correction turbine inlet temperature T.sub.1T′ is a value having high calculation accuracy in addition to fast reactivity by multiplying the turbine inlet temperature T.sub.1T with the correction factor X4. The gas turbine control device 20 determines a temperature of the combustion gas, based on the post correction turbine inlet temperature T.sub.1T′ having fast reactivity and high calculation accuracy. Therefore, the gas turbine control device 20 is capable of calculating the post correction turbine inlet temperature T.sub.1T′ having fast reactivity and high calculation accuracy. The correction unit 204 outputs the post correction turbine inlet temperature T.sub.1T′ to the fuel distribution ratio calculation unit 205.

(38) The fuel distribution ratio calculation unit 205 calculates a fuel distribution ratio for each of the fuel supply systems 50a, 50b, and 50c corresponding to the first fuel distribution ratio calculation unit 205a, the second fuel distribution ratio calculation unit 205b, and the third fuel distribution ratio calculation unit 205c, respectively (Step S105). The total of each fuel distribution ratio is 100%. A fuel distribution ratio for each of the fuel supply systems 50a, 50b, and 50c corresponding to the first fuel distribution ratio calculation unit 205a, the second fuel distribution ratio calculation unit 205b, and the third fuel distribution ratio calculation unit 205c, respectively, is calculated with a relationship formula of the post correction turbine inlet temperature T.sub.1T′ and a fuel distribution ratio. Note that the fuel supply systems 50a, 50b, and 50c each may have a different corresponding relationship of the post correction turbine inlet temperature T.sub.1T′ and a fuel distribution ratio. The first fuel distribution ratio calculation unit 205a, the second fuel distribution ratio calculation unit 205b, and the third fuel distribution ratio calculation unit 205c output the calculated opening amounts to the valve opening calculation unit 206. That is, the first fuel distribution ratio calculation unit 205a outputs a calculated fuel distribution ratio Da to the first valve opening calculation unit 206a. The second fuel distribution ratio calculation unit 205b outputs a calculated fuel distribution ratio db to the second valve opening calculation unit 206b. The third fuel distribution ratio calculation unit 205c outputs a calculated fuel distribution ratio Dc to the third valve opening calculation unit 206c.

(39) The first valve opening calculation unit 206a, the second valve opening calculation unit 206b, and the third valve opening calculation unit 206c corresponding to the fuel supply systems 50a, 50b, and 50c, respectively, substitute fuel distribution ratios of the corresponding fuel supply systems and a fuel control command value CSO in an opening amount calculation formula, and calculates opening degrees of the first regulating valve 18a, the second regulating valve 18b, and the third regulating valve 18c of the corresponding fuel supply systems, respectively (Step S106). That is, the first valve opening calculation unit 206a calculates an opening amount Oa of the first regulating valve 18a, the second valve opening calculation unit 206b calculates an opening amount Ob of the second regulating valve 18b, and the third valve opening calculation unit 206c calculates an opening amount Oc of the third regulating valve 18c. The first valve opening calculation unit 206a outputs the calculated opening amount Oa to the first regulating valve 18a, the second valve opening calculation unit 206b outputs the calculated opening amount Ob to the second regulating valve 18b, and the third valve opening calculation unit 206c outputs the calculated opening amount Oc to the third regulating valve 18c. The first regulating valve 18a, the second regulating valve 18b, and the third regulating valve 18c are controlled to have the input opening amounts, respectively.

(40) With the above-described process, an output value of the gas turbine 10 is not used for calculating the turbine inlet temperature. Instead, by substituting a measurement value having faster transient response in the abnormality physical model formula in the periphery of the combustor, the turbine inlet temperature is calculated. Therefore, even at the time of transient response in which output of the gas turbine 10 suddenly varies, an appropriate relationship of the turbine inlet temperature and the fuel distribution ratio can be maintained. Thus, stable combustion without causing combustion oscillation can be achieved.

Second Embodiment

(41) FIG. 4 is a function block diagram of a gas turbine control device according to a second embodiment.

(42) The gas turbine control device 20 may calculate a casing air temperature Tcs′, based on a measurement value. In the second embodiment, the casing air temperature Tcs' is calculated by substituting values of a compressor inlet pressure P.sub.1C and a compressor outlet pressure P.sub.2C that are used for calculating the compressor index pressure difference P.sub.index, the compressor inlet temperature T.sub.1ce, and a casing air temperature measurement value T.sub.CS_MEA in a casing air temperature calculation formula. The casing air temperature measurement value T.sub.CS_MEA is generally measured along with time delay. Thus, the measurement value T.sub.CS_MEA is not used as it is, and calculation with the above-described calculation formula for improving reactivity is performed. A casing air temperature calculation unit 207 substitutes the compressor inlet pressure P.sub.1C, the compressor outlet pressure P.sub.2C, the compressor inlet temperature T.sub.1c, and the casing air temperature measurement value T.sub.CS_MEA in a formula (4) being a physical model. With this, the casing air temperature Tcs′ having fast reactivity can be calculated according to a fourth function as noted below (Math.4).
[Math. 4]
T.sub.CS′=f(P.sub.1C,P.sub.2C,T.sub.1C,T.sub.CS_MEA)  (4)

(43) The above-described gas turbine control device 20 internally includes a computer system. Further, the steps of each process described above are stored in a computer readable recording medium in the form of a program, and the processes described above are implemented by the computer reading out and executing this program. Here, the computer readable storing medium refers to a magnetic disk, a magneto-optical disk, a CD-ROM, a DVD-ROM, a semiconductor memory, or the like. This computer program may be distributed to the computer on a communication line, and the computer that receives this distribution may execute the program.

(44) Also, the program as described above may realize some of the functions described above. In addition, the functions as described above may be realized in combination with a program already stored on the computer system, namely, a so-called differential file (differential program).

REFERENCE SIGNS LIST

(45) 10 Gas turbine 20 Gas turbine control device 40 Supply device 201 Fuel flow rate calculation unit 202 Air flow rate calculation unit 203 Turbine inlet temperature calculation unit 204 Correction unit 205 Fuel distribution ratio calculation unit 206 Valve opening calculation unit