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

Abstract

To provide a gas turbine control device capable of optimizing a control input by accurately predicting a future state quantity of a gas turbine even when an extracted steam flow rate from a compressor changes.SOLUTION: A gas turbine control device comprises: a prediction unit that predicts, using a prediction model, a future state quantity of a gas turbine corresponding to a control input to the gas turbine in a prediction horizon; an optimization unit that optimizes the control input in at least a part of the prediction horizon, using the prediction result; a storage unit that stores sensitivity information indicating sensitivity of the control input to a change speed of a state quantity for each of extracted steam flow rates from a compressor; and an update unit that reads, from the storage unit, the sensitivity information corresponding to the extracted steam flow rate assumed in the prediction horizon, and updates one or more coefficients of a prediction equation of the state quantity used in the prediction model, on the basis of the sensitivity information corresponding to the extracted steam flow rate.SELECTED DRAWING: Figure 7

Description

[Technical field]

[0001] The present disclosure relates to a gas turbine control device, a gas turbine control input determination method, and a gas turbine control input determination program. [Background technology]

[0002] Patent Document 1 discloses a gas turbine control device. This control device is configured to use a prediction model to predict future state quantities of the gas turbine corresponding to control inputs to the gas turbine in a prediction horizon, and to optimize the control inputs based on the prediction results. Specifically, it describes optimizing a fuel flow command and an inlet guide vane opening control command, which are control inputs to the gas turbine, while updating sensitivity information of the prediction model according to state quantities such as the turbine inlet temperature. [Prior art documents] [Patent documents]

[0003] [Patent Document 1] JP 2021-156260 A Summary of the Invention [Problem to be solved by the invention]

[0004] Incidentally, a gas turbine may be provided with an extraction pipe that returns a portion of the compressed air compressed by the compressor to the upstream side of the compressor, a turbine bypass pipe that bypasses the turbine and supplies a portion of the compressed air compressed by the compressor to the exhaust side of the turbine, etc. When a portion of the compressed air compressed by the compressor is extracted from the compressor using the above-mentioned extraction pipe or turbine bypass pipe, etc., even if, for example, the fuel flow rate or the compressor intake flow rate is the same, the amount of combustion air for the gas turbine changes depending on the extraction flow rate, and therefore, in the control device described in Patent Document 1, the prediction accuracy of the prediction model decreases.

[0005] In consideration of the above-described circumstances, at least one embodiment of the present disclosure has an object to provide a gas turbine control device, a gas turbine control input determination method, and a gas turbine control input determination program that are capable of accurately predicting future state variables of a gas turbine and optimizing a control input even if the extraction air flow rate from a compressor changes. [Means for solving the problem]

[0006] In order to achieve the above object, a gas turbine control device according to at least one embodiment of the present disclosure includes: A gas turbine control device for controlling a gas turbine including a compressor, a combustor, and a turbine, comprising: a prediction unit configured to predict a future state quantity of the gas turbine corresponding to a control input to the gas turbine in a prediction horizon using a prediction model; an optimization unit configured to optimize the control input over at least a portion of the prediction horizon using a prediction result from the prediction unit; a storage unit for storing sensitivity information indicating a sensitivity of the control input to a rate of change of the state quantity for each of the extracted air flow rates, where the flow rate of the compressed air extracted from the compressor is defined as an extracted air flow rate; an update unit configured to read out from the storage unit the sensitivity information corresponding to the bleed air flow rate assumed in the prediction horizon, and to update one or more coefficients of a prediction equation for the state quantity used in the prediction model, based on the sensitivity information corresponding to the bleed air flow rate; Equipped with.

[0007] In order to achieve the above object, a gas turbine control input determination method according to at least one embodiment of the present disclosure includes: a prediction step in which a prediction model predicts future state quantities of the gas turbine corresponding to control inputs to the gas turbine in a prediction horizon; an optimization step of optimizing the control input in at least a part of the prediction horizon using a prediction result in the prediction step; a step of storing, in a storage unit, sensitivity information indicating a sensitivity of the control input to a rate of change of the state quantity for each of the extracted air flow rates, where the flow rate of compressed air extracted from a compressor of the gas turbine is defined as an extracted air flow rate; reading out from the storage unit the sensitivity information corresponding to the bleed air flow rate assumed in the prediction horizon, and updating one or more coefficients of a prediction equation for the state quantity used in the prediction model based on the sensitivity information corresponding to the bleed air flow rate; Equipped with.

[0008] In order to achieve the above object, a gas turbine control input determination program according to at least one embodiment of the present disclosure includes: On the computer, predicting future state quantities of the gas turbine corresponding to control inputs to the gas turbine over a prediction horizon using a prediction model; a step of optimizing the control input over at least a portion of the prediction horizon using a prediction result from the step of predicting; a step of storing, in a storage unit, sensitivity information indicating a sensitivity of the control input to a rate of change of the state quantity for each of the extracted air flow rates, where the flow rate of compressed air extracted from a compressor of the gas turbine is defined as an extracted air flow rate; and a step of reading out from the storage unit the sensitivity information corresponding to the bleed air flow rate assumed in the prediction horizon, and updating one or more coefficients of a prediction equation for the state quantity used in the prediction model based on the sensitivity information corresponding to the bleed air flow rate; Execute the command. Effect of the Invention

[0009] According to at least one embodiment of the present disclosure, there are provided a gas turbine control device, a gas turbine control input determination method, and a gas turbine control input determination program that are capable of accurately predicting future state variables of a gas turbine and achieving high control performance even if the extraction air flow rate from a compressor changes. [Brief description of the drawings]

[0010] [Figure 1] FIG. 1 is a diagram illustrating a schematic configuration of a gas turbine system 1 according to an embodiment. [Diagram 2] 2 is a diagram illustrating an example of a hardware configuration of a control device 100 illustrated in FIG. [Diagram 3] 3 is a block diagram showing an example of a functional configuration of a control device 100 shown in FIGS. 1 and 2. FIG. [Figure 4] 1 is a schematic diagram for explaining the principle of optimization performed by a control device 100 according to an embodiment. [Diagram 5] 2 is a diagram showing an example of a performance table stored in a control device 100 according to an embodiment. FIG. [Figure 6] 4 is a graph showing an example of a relationship between a turbine inlet temperature T1T and a combustion efficiency η of a gas turbine 2 according to an embodiment. [Figure 7] FIG. 7 is a block diagram for illustrating an example of the function of the control device 100 described with reference to FIGS. 1 to 6. [Figure 8] FIG. 4 is a block diagram for explaining a modified example of the function of the control device 100. [Figure 9] FIG. 11 is a block diagram for explaining another modified example of the function of the control device 100. [Figure 10] 10 is a diagram for explaining an example of a method for correcting the target value trajectories of the gas turbine output GTMW and the turbine inlet temperature T1T by the target value trajectory correction unit 180. FIG. [Figure 11] 10 is a diagram for explaining an example of a method for correcting the target value trajectories of the gas turbine output GTMW and the turbine inlet temperature T1T by the target value trajectory correction unit 180. FIG. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0011] Hereinafter, some embodiments of the present disclosure will be described with reference to the accompanying drawings. However, the dimensions, materials, shapes, relative arrangements, etc. of components described as the embodiments or shown in the drawings are merely illustrative examples and are not intended to limit the scope of the invention. For example, expressions expressing relative or absolute configuration, such as "in a certain direction," "along a certain direction," "parallel," "orthogonal," "center," "concentric," or "coaxial," not only strictly express such a configuration, but also express a state in which there is a relative displacement with a tolerance or an angle or distance to the extent that the same function is obtained. For example, expressions indicating that things are in an equal state, such as "identical," "equal," and "homogeneous," not only indicate a state of strict equality, but also indicate a state in which there is a tolerance or a difference to the extent that the same function is obtained. For example, expressions describing shapes such as a rectangular shape or a cylindrical shape do not only refer to rectangular shapes, cylindrical shapes, etc. in the strict geometric sense, but also refer to shapes that include uneven portions, chamfered portions, etc., to the extent that the same effect is obtained. On the other hand, the expressions "comprise," "include," "have," "includes," or "have" of one element are not exclusive expressions excluding the presence of other elements.

[0012] (Gas Turbine System) Fig. 1 is a diagram illustrating a schematic configuration of a gas turbine system 1 according to one embodiment. As illustrated in Fig. 1, the gas turbine system 1 includes a gas turbine 2, a generator 7, and a control device 100 (gas turbine control device).

[0013] The gas turbine 2 is a power generating gas turbine and includes a compressor 3, a combustor 4, and a turbine 5 as main components.

[0014] The compressor 3 is connected to the turbine 5 via a rotating shaft 8A. The compressor 3 is rotationally driven by the rotational energy of the turbine 5 to generate compressed air. The compressed air generated by the compressor 3 is supplied to the combustor 4. An intake facility 10 (e.g., an intake duct provided with a filter) and an inlet guide vane 9 provided in the intake facility 10 are provided on the inlet side of the compressor 3, and an air flow rate G1, which is the intake flow rate of the compressor 3, is controlled by adjusting the opening of the inlet guide vane 9 by an actuator 35. The opening of the inlet guide vane 9 is controlled based on an inlet guide vane opening control command IGVCSO from the control device 100, and specifically, the actuator 35 adjusts the opening of the inlet guide vane 9 based on the inlet guide vane opening control command IGVCSO from the control device 100. Upstream of the inlet guide vane 9 in the intake equipment 10, a flowmeter 30 for measuring the air flow rate G1 and a thermometer 32 for measuring the compressor inlet temperature, i.e., the intake temperature T1C of the compressor 3, are provided.

[0015] The combustor 4 is connected to a fuel supply line 12, and burns the fuel supplied from the fuel supply line 12 with compressed air supplied from the compressor 3 to generate combustion gas, which is a working fluid for the turbine 5. The combustion gas generated in the combustor 4 is sent from the combustor 4 to the downstream turbine 5. The flow rate of the fuel supplied to the combustor 4 is adjusted by a fuel flow control valve 14 provided in the fuel supply line 12. The opening of the fuel flow control valve 14 is controlled based on a fuel flow command CSO (fuel flow control valve opening control command) from the control device 100. A flow meter 13 is provided in the fuel supply line 12 to measure a fuel flow rate Gf, which is the flow rate of the fuel supplied to the combustor 4.

[0016] The turbine 5 is driven by the combustion gas generated in the combustor 4. The turbine 5 is connected to the generator 7 by a rotating shaft 8B, and the generator 7 generates electricity using the rotational energy of the turbine 5. The generator 7 is provided with a power meter 34 configured to measure a gas turbine output GTMW, which is the output of the gas turbine 2. The power meter 34 may measure the amount of electricity generated by the generator 7 as the gas turbine output GTMW.

[0017] As shown in FIG. 1, the gas turbine 2 includes an extraction pipe 15 and a turbine bypass pipe 18. The extraction pipe 15 is configured to extract a part of the compressed air compressed by the compressor 3 from the compressor 3 and return it to the upstream side of the compressor 3 (the upstream side of the inlet guide vane 9 in the intake facility 10 in the illustrated example). The extraction pipe 15 is provided with an extraction valve 16 and a flow meter 17. The extraction valve 16 is a flow rate control valve that adjusts an extraction valve flow rate Gai, which is the flow rate of compressed air flowing through the extraction valve 16. The extraction valve flow rate Gai is the flow rate of compressed air extracted from the compressor 3 by the extraction pipe 15 (i.e., the flow rate of compressed air flowing through the extraction pipe 15). The opening of the extraction valve 16 is controlled based on an extraction valve opening control command AICSO (anti-icing command) from the control device 100. The flow meter 17 is configured to measure the extraction valve flow rate Gai.

[0018] Since the intake temperature T1C of the compressor 3 can be increased by increasing the opening of the bleed valve 16, the bleed valve 16 is also called an anti-icing valve. When anti-icing operation is performed by operating the gas turbine 2 with the bleed valve 16 open, high-temperature, high-pressure air at the outlet of the compressor 3 is bled through the bleed pipe 15 and injected into the inlet of the compressor 3, thereby controlling the intake temperature T1C of the compressor 3 to a target temperature. The bleed valve flow rate Gai is controlled by the bleed valve opening control command AICSO, but the combustion air flow rate, which is the flow rate of air supplied to the combustor 4, is reduced by the bleed from the compressor 3, so that the turbine inlet temperature T1T increases with the bleed unless the inlet guide vane opening control command IGVCSO is changed. Therefore, the control device 100 uses model predictive control as described below to control the intake air temperature T1C to a target temperature while simultaneously controlling the inlet guide vane opening control command IGVCSO, thereby achieving control of the desired gas turbine output GTMW, the turbine inlet temperature T1T, and the intake air temperature T1C.

[0019] The turbine bypass piping 18 is configured to supply a portion of the compressed air compressed by the compressor 3 to the downstream side of the turbine 5 (in the illustrated example, a turbine exhaust piping 20 through which the exhaust gas from the turbine 5 flows) while bypassing the combustor 4 and the turbine 5. The turbine bypass piping 18 is provided with a turbine bypass valve 22 and a flow meter 24. The turbine bypass valve 22 is a flow rate control valve that adjusts a turbine bypass valve flow rate Gtb, which is the flow rate of the compressed air flowing through the turbine bypass valve 22. The turbine bypass valve flow rate Gtb is the flow rate of the compressed air extracted from the compressor 3 by the turbine bypass piping 18 (i.e., the flow rate of the compressed air flowing through the turbine bypass piping 18). The opening of the turbine bypass valve 22 is controlled based on a turbine bypass valve opening control command TBCSO from the control device 100. The flow meter 24 is configured to measure the turbine bypass valve flow rate Gtb.

[0020] Since the turbine inlet temperature T1T can be increased by increasing the opening of the turbine bypass valve 22, the turbine bypass valve 22 is used for the purpose of increasing the turbine inlet temperature T1T and improving efficiency and emissions when the gas turbine 2 is at a low load, particularly when the opening of the inlet guide vane 9 is at or near the minimum opening. The control device 100 improves the controllability of the turbine inlet temperature T1T in the low load range of the gas turbine 2 by controlling the turbine bypass valve opening control command TBCSO together with the inlet guide vane opening control command IGVCSO using model predictive control as described later.

[0021] In addition, a thermometer 26 is provided in the final stage of the turbine 5 for measuring the blade pass temperature TBP, which is the temperature of the combustion gas that has passed through the final stage blades, and a thermometer 28 is provided in the turbine exhaust piping 20 downstream of the thermometer 26 for measuring the exhaust gas temperature TEX, which is the temperature of the turbine exhaust gas.

[0022] FIG. 2 is a diagram illustrating an example of a hardware configuration of the control device 100 illustrated in FIG. As shown in FIG. 2, the control device 100 includes, for example, a processor 72, a RAM (Random Access Memory) 74, a ROM (Read Only Memory) 76, a HDD (Hard Disk Drive) 78, an input I / F 80, and an output I / F 82, and is configured using a computer in which these are connected to each other via a bus 84. The hardware configuration of the control device 100 is not limited to the above, and may be configured by a combination of a control circuit and a storage device. The control device 100 is configured by a computer executing a program that realizes each function of the control device 100. The functions of each part in the control device 100 described below are realized by, for example, loading a program (for example, a control input determination program described later) held in the ROM 76 into the RAM 74 and executing it with the processor 72, and reading and writing data from and to the RAM 74 and the ROM 76. The hardware that configures the control device 100 may be concentrated in one place, or may be distributed across multiple places.

[0023] FIG. 3 is a block diagram showing an example of a functional configuration of the control device 100 shown in FIGS. As shown in FIG. 3, the control device 100 includes a prediction unit 110, an optimization unit 120, a storage unit 130, an update unit 140, a state quantity acquisition unit 150, a target value trajectory acquisition unit 160, and a table update unit 170.

[0024] The prediction unit 110 is configured to predict future state quantities of the gas turbine 2 corresponding to control inputs to the gas turbine 2 in a prediction horizon using a prediction model. The control inputs in this prediction model include a fuel flow rate command CSO, an inlet guide vane opening control command IGVCSO, an extraction valve opening control command AICSO, and a turbine bypass valve opening control command TBCSO. The state quantities used in this prediction model include, for example, information indicating each of the fuel flow rate Gf, the air flow rate G1, the extraction valve flow rate Gai, the turbine bypass valve flow rate Gtb, the power output GTMW, the turbine inlet temperature T1T, the blade path temperature TBP, the exhaust gas temperature TEX, and the intake temperature T1C in the gas turbine 2.

[0025] The optimization unit 120 is configured to optimize the control input in at least a part of the prediction horizon using the prediction result in the prediction unit 110. The optimization unit 120 inputs the control input in the prediction horizon to the prediction unit 110 to obtain the prediction result.

[0026] If the flow rate of compressed air bleed from the compressor 3 is defined as the bleed flow rate, the memory unit 130 stores sensitivity information indicating the sensitivity of the control input to the rate of change of the state quantity for each bleed flow rate. The memory unit 130 also stores constraint conditions, a prediction model, and a control input determination program, which will be described later. The sensitivity information includes first sensitivity information and second sensitivity information, as will be described later. The sensitivity information may be a performance table, which will be described later. In this embodiment, an example will be described in which the "bleed flow rate" is the sum of the bleed valve flow rate Gai and the turbine bypass valve flow rate Gtb.

[0027] The update unit 140 reads out sensitivity information according to the bleed air flow rate assumed in the prediction horizon from the storage unit 130, and updates one or more coefficients of a prediction equation for a state quantity used in the prediction model based on the sensitivity information according to the bleed air flow rate. The sensitivity information according to the bleed air flow rate assumed in the prediction horizon may be sensitivity information corresponding to the current bleed air flow rate.

[0028] The state quantity acquisition unit 150 is configured to acquire state quantities of the gas turbine 2 from sensors and the like. The state quantity acquisition unit 150 acquires, for example, a fuel flow rate Gf measured by the flowmeter 13, an air flow rate G1 measured by the flowmeter 30, an extraction valve flow rate Gai measured by the flowmeter 17, a turbine bypass valve flow rate Gtb measured by the flowmeter 24, a gas turbine output GTMW measured by the output meter 34, a turbine inlet temperature T1T measured by the thermometer 36, a blade path temperature TBP measured by the thermometer 26, an exhaust gas temperature TEX measured by the thermometer 28, and an intake air temperature T1C acquired by the thermometer 32.

[0029] The target value trajectory acquisition unit 160 is configured to acquire the target value trajectories of each of the gas turbine output GTMW, the turbine inlet temperature T1T, and the intake air temperature T1C, among the above-mentioned state quantities.

[0030] The table update unit 170 is configured to update the performance table stored in the storage unit 130 based on performance data acquired during actual operation. For example, the initial performance table is created based on design information. After that, the performance table is updated based on performance data acquired during actual operation.

[0031] (Prediction model) We will now explain the prediction model used by the prediction unit 110. In some embodiments, the prediction model may be configured to use a state space model described by a linear model, as described below.

[0032] First, the definitions of the terms in the prediction formula will be explained. x(k) is a state vector at the current time k. x(k) may be considered as a matrix whose elements are multiple state quantities. u(k) is a control input vector at the current time k. y(k) is an observation output vector at the current time k. y(k) may be considered as a matrix whose elements are control quantities for some of the multiple state quantities. z(k) is an evaluation output vector at the current time k. A, B, C y , C z are matrices that make up the state equation. Matrices A and B have elements that are one or more coefficients included in the prediction formula. Matrices A and B are updated when one or more coefficients are updated.

[0033] With these definitions, the following relational expressions can be considered to hold as functions of time. These expressions can be used as prediction expressions. x(k+1)=Ax(k)+Bu(k) (1) y(k)=C y x(k) (2) z(k)=C z x(k) (3)

[0034] In formula (1), x(k+1) is a future state vector. The unit of time may be a control period. In other words, x(k+1) may be a future state vector expected when the current time k is advanced by one control period.

[0035] According to formula (1), it is possible to calculate the future state vector x(k+1) based on the state vector x(k) at the current time k, the control input vector u(k) at the current time k, and the matrices A and B. By repeating such calculations, it is also possible to obtain the future state vector x(k+i) expected when the current time k is advanced by any control period i.

[0036] According to equation (2), the state vector x(k) and matrix C y According to the formula (3), the observation output vector y(k) at the current time k can be calculated based on the state vector x(k) and the matrix C z Based on this, it is possible to calculate the evaluation output vector z(k) at the current time k. In addition, by combining equations (1) to (3), it is also possible to calculate the observation output vector y(k+i) and the evaluation output vector z(k+i) corresponding to the future state quantity vector x(k+i).

[0037] Here, the difference from the equilibrium point (which may be close to the equilibrium point) of the thermal equilibrium in the operation of the gas turbine 2 may be considered as the state quantity x. In this case, the state quantity vector x may be defined as a matrix having elements dGf, dG1, dGai, dGtb, dGTMW, dT1T, dTBP, dTEX, and dT1C, which are the differences from the equilibrium points in each of the fuel flow rate Gf, the air flow rate G1, the bleed valve flow rate Gai, the turbine bypass valve flow rate Gtb, the output GTMW, the turbine inlet temperature T1T, the blade path temperature TBP, the exhaust gas temperature TEX, and the intake temperature T1C. The control input u may be defined as a matrix having elements dCSO, dIGVCSO, dAICSO, and dTBCSO, which correspond to the state quantities (differences from the equilibrium points) of the fuel flow rate command CSO, the inlet guide vane opening control command IGVCSO, the bleed valve opening control command AICSO, and the turbine bypass valve opening control command TBCSO.

[0038] More specifically, dGf, dG1, dGai, dGtb, dGTMW, dT1T, dTBP, dTEX, and dT1C can be expressed by the following differential equations. For example, Gf-Gf0 corresponds to the difference of the fuel flow rate Gf from the equilibrium point, and dGf corresponds to d(Gf-Gf0). dx / dt indicates the rate of change of the state quantity x. For example, d(Gf-Gf0) / dt indicates the rate of change of the fuel flow rate Gf. The same can be considered for other state quantities.

[0039] (Gf-Gf0)+τ Gf ×d(Gf-Gf0) / dt=K_Gf_CSO×(CSO-CSO0)...(4) (G1-G10)+τ G1 ×d(G1-G10) / dt=K_G1_IGVCSO×(IGVCSO-IGVCSOo)···(5) (Gai-Gai0)+τ Gai ×d(Gai-Gai0) / dt=K_Gai_AICSO×(AICSO-AICSO0)···(6) (Gtb-Gtb0)+τ Gtb ×d(Gtb-Ttb0) / dt=K_Gtb_TBCSO×(TBCSO-TBCSO0)...(7) (GTMW-GTMW0)+τ GTMW ×d(GTMW-GTMW0) / dt=K_GTMW_Gf×(Gf-Gf0)+K_GTMW_G1×(G1-G10)+K_GTMW_Gai×(Gai-Gai0)+K_GTMW_Gtb×(Gtb-Gtb0)...(8) (T1T-T1T0)+τ T1T ×d(T1T-T1T0) / dt=K_T1T_Gf×(Gf-Gf0)+K_T1T_G1×(G1-G10)+K_T1T_Gai×(Gai-Gai0)+K_T1T_Gtb×(Gtb-Gtb0)...(9) (TBP-TBP0)+τ TBP ×d(TBP-TBP0) / dt=K_TBP_Gf×(Gf-Gf0)+K_TBP_G1×(G1-G10)+K_TBP_Gai×(Gai-Gai0)+K_TBP_Gtb×(Gtb-Gtb0)...(10) (TEX-TEX0)+τ TEX ×d(TEX-TEX0) / dt=K_TEX_Gf×(Gf-Gf0)+K_TEX_G1×(G1-G10)+K_TEX_Gai×(Gai-Gai0)+K_TEX_Gtb×(Gtb-Gtb0)...(11) (T1C-T1C0)+τ T1C ×d(T1C-T1C0) / dt=K_T1C_Gai×(Gai-Gai0)···(12)

[0040] Equations (4) to (7) are first sensitivity information indicating the relationship between the rate of change of the state quantity x and the control input u. From these equations, the sensitivity of the control input u to the rate of change of the state quantity x can be seen. Equations (8) to (11) are second sensitivity information indicating the relationship between the first state quantities Gf, G1, Gai, and Gtb and the rate of change of the second state quantities GTMW, T1T, TBP, and TEX. Equation (12) is second sensitivity information indicating the relationship between the first state quantity Gai and the rate of change of the second state quantity T1C. τ Gf , τ G1 , τ Gai , τ Gtb , τ GTMW , τ T1T , τTBP and τ TEX and τ T1C are the respective time constants.

[0041] These sensitivity information include sensitivity coefficients K_Gf_CSO, K_G1_IGVCSO, K_Gai_AICSO, K_Gtb_TBCSO, K_GTMW_Gf, K_GTMW_G1, K_GTMW_Gai, K_GTMW_Gtb, K_T1T_Gf, K_T1T_G1, K_T1T_Gai, K_T1T_Gtb, K_TBP_Gf, K_TBP_G1, K_TBP_Gai, K_TBP_Gtb, K_TEX_Gf, K_TEX_G1, K_TEX_Gai, K_TEX_Gtb, and K_T1C_Gai. For example, the sensitivity coefficient K_Gf_CSO is a coefficient that indicates the sensitivity of the fuel flow rate Gf to a change in the fuel flow rate command CSO. Similarly, the other sensitivity coefficients are coefficients that indicate the sensitivity of a previous parameter to a change in a subsequent parameter. As described later, among these sensitivity coefficients, the second sensitivity coefficients K_GTMW_Gf, K_T1T_Gf, K_TBP_Gf, and K_TEX_Gf may be corrected by taking into account the combustion efficiency η.

[0042] When equations (4) to (12) are rearranged into a matrix form, the state equation obtained by time-differentiating the state quantity vector x is expressed by the following equation (13). When the observation output vector y is put into a matrix form, it is expressed by the following equation (14). When the evaluation output vector z is put into a matrix form, it is expressed by the following equation (15). TIFF2025034955000002.tif65170TIFF2025034955000003.tif54170TIFF2025034955000004.tif33170

[0043] (optimization) An example of optimization by the optimization unit 120 will be specifically described below. Fig. 4 is a schematic diagram for explaining the principle of optimization executed by the control device 100 according to one embodiment. The upper graph in Fig. 4 shows the time progression of the controlled variables and the like, and the lower graph shows the time progression of the control input. The controlled variables are, for example, the output GTMW of the gas turbine 2, the turbine inlet temperature T1T, and the intake air temperature T1C. Note that the upper graph in Fig. 4 shows the progress of the output GTMW as an example of the controlled variable.

[0044] In the upper graph of Fig. 4, y(t) shown by a solid line indicates the observed output, and z(t) shown by a dotted line indicates the evaluation output. In other words, y(t) indicates the progress of the actual results up to the current time k, and z(t) indicates the progress of the predicted results in the future after the current time k (i.e., time k+i). r(k+i|k) indicates the operating line, i.e., the target value trajectory. This target value trajectory is acquired by the target value trajectory acquisition unit 160. It is preferable that the target value trajectory and the trajectory of the controlled variable are close to each other.

[0045] The horizontal axis is the time axis, and the time axis is the same for the upper and lower graphs in Figure 4. On the time axis, the section from the evaluation start time k+Hw to time k+HP is the prediction horizon. The section from the current time k to time k+Hu is the control horizon. The control horizon is a shorter time than the prediction horizon.

[0046] In the lower graph of Fig. 4, the actual progress of the control input u(t) is shown by a solid line, and the progress of the control input u(t) after the current time k is shown by a dashed line. The progress of the evaluation output z(t) in the upper graph is determined as the controlled variable predicted according to the control input u(t) after the current time k in the lower graph.

[0047] The optimization unit 120 determines the control input u(t) so that the value of the control input u(t) in the control horizon changes, and determines the control input u(t) so that the value of the control input u(t) in the section after the control horizon in the prediction horizon is constant. For example, as shown in Fig. 4, the value of the control input u(t) changes in the section from the current time k to time k+Hu, whereas the value of the control input u(t) is constant in the section from time k+Hu to time k+Hp.

[0048] This constant value may be the value of the immediately preceding control input u(k+Hu). In other words, calculation of the constant value may be omitted. Note that the constant value may be the value of the control input u(t) calculated so that the evaluation output z(t) matches the reference r(k+i|k) indicating the target value trajectory at the end of the prediction horizon (time k+Hp), as shown in FIG. 4.

[0049] The optimization unit 120 inputs the control input in the determined prediction horizon to the prediction unit 110. As a result, the transition of the evaluation output z(t) is obtained. Note that the control input determined before input to the prediction unit 110 is not determined as the optimal control input, but is provisionally determined as a candidate for the optimal control input. That is, the dashed line portion in the graph at the bottom of FIG. 4 indicates the control input u(t) during the optimization calculation.

[0050] For example, if the control input u(t) actually used is for one control cycle immediately following the current time (i.e., from the current time k to time k+Hw), the prediction horizon is from one control cycle after the current time to 11 control cycles later, and the control horizon is for five control cycles from the current time. In this case, calculated values ​​are set for the control inputs up to five moves ahead, and constant values ​​are set for the control inputs for the next five moves ahead. Based on these control inputs, state quantities and control quantities for up to 11 moves ahead are predicted.

[0051] The optimization unit 120 may be configured to optimize the control input u(t) based on an evaluation function that depends on the deviation between the target value trajectory and the control amount (parameter). For example, the optimization unit 120 may determine the optimal control input u(t) by solving an optimization problem that minimizes the evaluation function while referring to the above-mentioned prediction result.

[0052] Here, the evaluation function may be a function including a first term indicating an error between the target value trajectory and the controlled variable in the prediction horizon, and a second term indicating an accumulated value of the change in the control input u(t) in the control horizon. For example, the evaluation function may be V shown in the following formula (16). TIFF2025034955000005.tif11170

[0053] The evaluation function V shown in equation (16) is a weighting matrix in which each l2 norm (distance) is weighted. Q is a state weighting matrix, and R is a control weighting matrix. Specifically, Q is a matrix in which the weights for the output GTMW of the gas turbine 2, the weights for the fuel-air ratio, i.e., the turbine inlet temperature T1T, the weights for the blade path temperature TBP, the weights for the exhaust gas temperature TEX, and the intake temperature T1C are arranged on the diagonal elements, and the other matrix elements are set to zero. R is a matrix in which the weights for the fuel flow command CSO change, the weights for the inlet guide vane opening control command IGVCSO change, the weights for the bleed valve opening control command AICSO change, and the weights for the turbine bypass valve opening control command TBCSO change are arranged on the diagonal elements, and the other matrix elements are set to zero. These Q and R are parameters that affect the controllability. If Q is made larger relative to R, the responsiveness in the control improves, but the stability decreases. Conversely, if R is made larger relative to R, the responsiveness in the control decreases, but the stability improves. In addition, since the turbine bypass valve 22 is operated when the opening of the inlet guide vane 9 is at the minimum opening or an opening close to the optimum opening, when the inlet guide vane 9 is not at the minimum opening or an opening close to the minimum opening, the control weight of the TBCSO may be set large so that the turbine bypass valve 22 does not operate.

[0054] The optimization unit 120 may be configured to optimize the control input u(t) based on constraint conditions for constraining one or more of the rate of change of the control input u(t), the absolute value of the control input u(t), and the absolute value of the state quantity x. In this case, for example, the control input u(t) can be determined so as to satisfy the constraint conditions on the operation of the gas turbine 2.

[0055] As a constraint condition for the rate of change of the control input u(t), for example, a constraint condition for satisfying −10(% / sec)<rate of change of fuel flow rate command CSO<10(% / sec) may be set. As a constraint condition for the rate of change of the control input u(t), for example, a constraint condition for satisfying −10(deg / sec)<rate of change of inlet guide vane opening control command IGVCSO<10(deg / sec) may be set.

[0056] As a constraint condition for the absolute value of the control input u(t), for example, a constraint condition for satisfying 0(%)<absolute value of fuel flow rate command CSO<100(%) may be set. As a constraint condition for the absolute value of the control input u(t), for example, a constraint condition for satisfying −8(deg)<absolute value of inlet guide vane opening control command IGVCSO<46(deg). As a constraint condition for the absolute value of the state quantity x, for example, a constraint condition for satisfying the condition that the turbine inlet temperature is less than 1650° C. may be set. These constraint conditions may be applied in combination.

[0057] The control device 100 may be configured to execute, for each control period, a process of predicting the state quantity x, a process of updating one or more coefficients included in the matrices A and B, and a process of optimizing the control input u(t). That is, the control device 100 may execute these processes every time one control period elapses.

[0058] (Coefficient update) Hereinafter, a specific example of updating the coefficients of the prediction equation by the update unit 140 will be described. As described above, the update unit 140 reads out sensitivity information corresponding to the bleed air flow rate assumed in the prediction horizon (i.e., sensitivity information determined according to the bleed valve flow rate Gai and the turbine bypass valve flow rate Gtb) from the storage unit 130, and updates one or more coefficients of the prediction equation of the state quantity used in the prediction model based on the sensitivity information corresponding to the bleed air flow rate.

[0059] The sensitivity information according to the bleed air flow rate assumed in the prediction horizon may be sensitivity information according to the current bleed air flow rate. The sensitivity information may be a performance table of the gas turbine 2. Fig. 5 is a diagram showing an example of the performance table stored in the control device 100 according to one embodiment.

[0060] For example, the performance table may be an information table in which information indicating the inlet guide vane position control command IGVCSO, information indicating the output GTMW of the gas turbine 2, information indicating the turbine inlet temperature T1T, information indicating the exhaust gas temperature TEX, information indicating the bleed valve flow rate Gai, information indicating the turbine bypass valve flow rate Gtb, the intake air temperature T1C, and the atmospheric temperature are associated with each other as state quantities. If the sum of the bleed valve flow rate Gai and the turbine bypass valve flow rate Gtb is defined as the bleed air flow rate Gaitb, as shown in Fig. 5, the storage unit 130 may store an information table for each bleed air flow rate Gaitb, and multiple state quantities using the bleed air flow rate Gaitb as a condition may be associated with each other in the information table.

[0061] In the example shown in Fig. 5, a plurality of state quantities are associated when the bleed air flow rate Gaitb is 10% of the air flow rate. In addition to the state quantities shown in Fig. 5, other state quantities (e.g., air flow rate G1, fuel flow rate Gf, etc.) may also be associated in the performance table. Although omitted in Fig. 5, the performance table makes it possible to read the value of another state quantity when a certain state quantity is changed. In other words, the values ​​are arranged vertically in the table.

[0062] The performance table may be discrete data of values ​​indicating each state quantity as shown in Fig. 5. However, as the extracted air flow rate Gaitb, not only the discrete values ​​shown in the performance table but also a linear interpolation value between them may be used. As described above, the performance table stored in the storage unit 130 may be updated by the table update unit 170 based on the actual data during operation.

[0063] The performance table includes, for example, a plurality of sensitivity coefficients K_Gf_CSO, K_G1_IGVCSO, K_Gai_AICSO, K_Gtb_TBCSO, K_GTMW_Gf, K_GTMW_G1, K_GTMW_Gai, K_GTMW_Gtb, K_T1T_Gf, K_T1T_G1, K_T1T_Gai, K_T1T_Gtb, K_TBP_Gf, K_TBP_G1, K_TBP_Gai, K_TBP_Gtb, K_TEX_Gf, K_TEX_G1, K_TEX_Gai, K_TEX_Gtb, and K_T1C_Gai. As described above, the sensitivity coefficient is a coefficient indicating the sensitivity of a previous parameter to a change in a subsequent parameter indicated by a symbol. These sensitivity coefficients can be obtained by dividing the numerical values ​​in the performance table.

[0064] The update unit 140 can obtain one or more coefficients included in the prediction formula from such a performance table, i.e., sensitivity information. The update unit 140 can update the coefficients of the prediction formula (i.e., matrices A and B) using the obtained coefficients.

[0065] Fig. 6 is a graph showing an example of a relationship between the turbine inlet temperature T1T and the combustion efficiency η of the gas turbine 2 according to one embodiment. As shown in Fig. 6, the combustion efficiency η varies depending on the turbine inlet temperature T1T. The variation in the combustion efficiency η significantly affects the sensitivity coefficient (second sensitivity coefficient) related to the fuel flow rate Gf.

[0066] Therefore, it is preferable that the update unit 140 corrects a coefficient (second sensitivity coefficient) indicating the sensitivity of one or more state quantities among the power GTMW, the turbine inlet temperature T1T, the blade path temperature TBP, and the exhaust gas temperature TEX to a change in the fuel flow rate Gf based on the rate of change of the combustion efficiency η with respect to the turbine inlet temperature T1T. For example, the update unit 140 may be configured to correct the second sensitivity coefficients K_GTMW_Gf, K_T1T_Gf, K_TBP_Gf, and K_TEX_Gf among the sensitivity coefficients included in the performance table by taking into account the combustion efficiency η of the gas turbine 2.

[0067] (Functions of the control device 100) The following describes the functions of the control device 100. Fig. 7 is a block diagram for explaining an example of the functions of the control device 100 described with reference to Figs.

[0068] First, the control device 100 acquires, as observed outputs, the fuel flow rate Gf measured by the flowmeter 13, the air flow rate G1 measured by the flowmeter 30, the bleed valve flow rate Gai measured by the flowmeter 17, the turbine bypass valve flow rate Gtb measured by the flowmeter 24, the gas turbine output GTMW measured by the output meter 34, the turbine inlet temperature T1T measured by the thermometer 36, the blade path temperature TBP measured by the thermometer 26, the exhaust gas temperature TEX measured by the thermometer 28, and the intake temperature T1C measured by the thermometer 32. The state quantity acquisition unit 150 acquires state quantities of the fuel flow rate Gf, the air flow rate G1, the bleed valve flow rate Gai, the turbine bypass valve flow rate Gtb, the gas turbine output GTMW, the turbine inlet temperature T1T, the blade path temperature TBP, the exhaust gas temperature TEX, and the intake temperature T1C from these observed outputs. The observed outputs may also be used as input data of a disturbance observer in the optimization unit 120. The output data of the disturbance observer may be used in formulation in optimization (i.e., updating the coefficients of the prediction formula and calculating the evaluation function).

[0069] The control device 100 also acquires the current target value and the final target value for each of the gas turbine output GTMW, the turbine inlet temperature T1T, and the intake air temperature T1C. These may be values ​​based on values ​​determined by a user operation. The target value trajectory acquisition unit 160 acquires the target value trajectories of the gas turbine output GTMW, the turbine inlet temperature T1T, and the intake air temperature T1C by calculating the difference between the final target value and the current target value for each of the gas turbine output GTMW, the turbine inlet temperature T1T, and the intake air temperature T1C. The prediction unit 110 predicts the future state quantity using a prediction model based on the state quantity acquired by the state quantity acquisition unit 150. At this time, the update unit 140 updates the coefficient of the prediction formula of the prediction model based on the sensitivity information stored in the storage unit 130. The update unit 140 also formulates the calculation formula of the evaluation function V using the target value trajectory acquired by the target value trajectory acquisition unit 160 and the prediction model. The optimization unit 120 calculates an evaluation function V while acquiring the prediction result of the prediction model, and optimizes the control inputs so as to determine the control inputs that minimize the evaluation function V (here, a combination of the above four control commands, i.e., a combination of the fuel flow command CSO, the inlet guide vane opening control command IGVCSO, the bleed valve opening control command AICSO, and the turbine bypass valve opening control command TBCSO) as optimal control inputs. When the optimal control inputs are determined, those control inputs (i.e., CSO, IGVCSO, AICSO, and TBCSO) are output.

[0070] In this way, the control device 100 determines the optimal control input. Such a control input determination method may be executed based on a control input determination program. Furthermore, the control device 100 may be configured to execute control of the gas turbine 2 in accordance with the control input optimized by the optimization unit 120.

[0071] The following describes the effects of the above-described control device 100. According to the control device 100, one or more coefficients of the prediction formula can be updated according to the bleed valve flow rate Gai and the turbine bypass valve flow rate Gtb. Therefore, even if the bleed valve flow rate Gai or the turbine bypass valve flow rate Gtb changes, the future state quantities of the gas turbine 2 can be accurately predicted using a prediction model that uses coefficients according to the bleed valve flow rate Gai and the turbine bypass valve flow rate Gtb, making it possible to realize high control performance.

[0072] In addition, by using not only the CSO and the IGVCSO but also the AICSO and the TBCSO as control inputs for bringing the gas turbine output GTMW and the turbine inlet temperature T1T closer to their target value trajectories, respectively, it is possible to calculate optimal control inputs even when the extraction valve 16 and the turbine bypass valve 22 are in operation, and high controllability can be achieved.

[0073] Furthermore, when controlling the bleed valve 16, the CSO, the IGVCSO, the AICSO, and the TBCSO can be optimized in consideration of the target value trajectory of the intake air temperature T1C of the compressor 3 as well as the gas turbine output GTMW and the turbine inlet temperature T1T, and therefore the control accuracy of the intake air temperature T1C can be improved.

[0074] Fig. 8 is a block diagram for explaining a modified example of the functions of the control device 100. In the following, of the functions of the control device 100 shown in Fig. 8, explanations of parts that overlap with the functions explained using Figs. 3 to 7 etc. will be omitted.

[0075] 8, the control device 100 determines the extraction valve position control command AICSO and the turbine bypass valve position control command TBCSO by a means other than the prediction model outside the model predictive control. The control device 100 may calculate the AICSO and TBCSO by, for example, feedforward control or feedback control, separately from the model predictive control, and the calculation method of the AICSO and TBCSO is not particularly limited.

[0076] 8, the update unit 140 reads out sensitivity information corresponding to the bleed valve flow rate Gai and the turbine bypass valve flow rate Gtb from the storage unit 130 based on the bleed valve flow rate Gai and the turbine bypass valve flow rate Gtb acquired as observation outputs by the control device 100, and updates one or more coefficients of a prediction equation for a state quantity used in a prediction model based on the sensitivity information corresponding to the bleed valve flow rate Gai and the turbine bypass valve flow rate Gtb. The prediction unit 110 predicts a future state quantity using a prediction model including the updated coefficients, thereby taking into account the influences of AICSO and TBCSO. In this case, the state quantities input to the prediction model do not include the bleed valve flow rate Gai and the turbine bypass valve flow rate Gtb, and the state quantity acquisition unit 150 acquires, from the fuel flow rate Gf, air flow rate G1, bleed valve flow rate Gai, turbine bypass valve flow rate Gtb, gas turbine output GTMW, turbine inlet temperature T1T, blade path temperature TBP, exhaust gas temperature TEX, and intake temperature T1C acquired by the control device 100 as observation outputs, only the fuel flow rate Gf, the air flow rate G1, the gas turbine output GTMW, the turbine inlet temperature T1T, the blade path temperature TBP, and the exhaust gas temperature TEX as state quantities to be input to the prediction model.

[0077] 8, the target value trajectory acquiring unit 160 acquires the target value trajectories of the gas turbine output GTMW and the turbine inlet temperature T1T, but does not acquire the target value trajectory of the intake air temperature T1C. The updating unit 140 formulates an arithmetic expression for the evaluation function V by using the target value trajectory acquired by the target value trajectory acquiring unit 160 and the updated prediction model.

[0078] The optimization unit 120 calculates an evaluation function V while acquiring the prediction result of the prediction model, and optimizes the control input so as to determine the control input (here, a combination of the fuel flow rate command CSO and the inlet guide vane opening control command IGVCSO) that minimizes the evaluation function V as the optimal control input. When the optimal control inputs are determined, those control inputs (CSO and IGVCSO) are output from the optimizer 120.

[0079] In the embodiment described using Fig. 7 etc., equations (13) to (15) were shown as examples of prediction formulas for the prediction model, but in the embodiment shown in Fig. 8, as described above, the state quantities and control inputs used in the prediction model are different from those in the embodiment described using Fig. 7 etc. For this reason, in the embodiment shown in Fig. 8, the control device 100 may use, for example, the following equations (17) to (19) as prediction formulas for the prediction model.

[0080] TIFF2025034955000006.tif71170

[0081] TIFF2025034955000007.tif47170

[0082] TIFF2025034955000008.tif36170

[0083] In this case, as shown in equation (17), the performance table as sensitivity information read out from memory unit 130 may include K_Gf_CSO, K_G1_IGVCSO, K_GTMW_Gf, K_GTMW_G1, K_T1T_Gf, K_T1T_G1, K_TBP_Gf, K_TBP_G1, K_TEX_Gf and K_TEX_G1 as multiple sensitivity coefficients included in the prediction equation of the prediction model.

[0084] 8, even if the bleed valve flow rate Gai or the turbine bypass valve flow rate Gtb changes, the future state quantities of the gas turbine can be accurately predicted using a prediction model using coefficients according to the bleed valve flow rate Gai and the turbine bypass valve flow rate Gtb, and high control performance can be realized. In particular, when the movement of the AICSO is relatively small compared to the movements of the other control commands (CSO and IGVCSO), the behavior of the other control commands within the prediction horizon can be ignored, and therefore high control performance can be expected.

[0085] Fig. 9 is a block diagram for explaining another modified example of the functions of the control device 100. In the following, among the functions of the control device 100 shown in Fig. 9, the description of parts that overlap with the functions explained using Fig. 8 etc. will be omitted.

[0086] In the embodiment shown in FIG. 9, the basic functions of the control device 100 are the same as those described using FIG. 8 etc., but differ from the embodiment shown in FIG. 8 in that the control device 100 further includes a target value trajectory correction unit 180 in addition to the units shown in FIG. 3.

[0087] 9, the target value trajectory acquisition unit 160 acquires the target value trajectories of the gas turbine output GTMW and the turbine inlet temperature T1T, but does not acquire the target value trajectory of the intake air temperature T1C. The target value trajectory correction unit 180 acquires a bleed valve flow rate future prediction indicating the transition of a predicted value for the bleed valve flow rate Gai in a prediction horizon, and a turbine bypass valve flow rate future prediction indicating the transition of a predicted value for the turbine bypass valve flow rate Gtb in a prediction horizon. The bleed valve flow rate future prediction here may be, for example, a future predicted value of the bleed valve flow rate Gai, a future predicted value of the bleed valve opening control command AICSO, or a future predicted value of the ratio of the bleed valve flow rate Gai to the air flow rate G1. Further, the future turbine bypass valve flow rate prediction may be, for example, a future predicted value of a turbine bypass valve flow rate Gtb, a future predicted value of a turbine bypass valve opening control command TBCSO, or a future predicted value of a turbine bypass ratio (hereinafter, referred to as "TB ratio") which is a ratio of the turbine bypass valve flow rate Gtb to the air flow rate G1. The target value trajectory correction unit 180 corrects the target value trajectories of the gas turbine output GTMW and the turbine inlet temperature T1T over the prediction horizon based on the acquired future extraction valve flow rate prediction and future turbine bypass valve flow rate prediction.

[0088] 10 and 11 are diagrams for explaining an example of a method for correcting the target value trajectories of the gas turbine output GTMW and the turbine inlet temperature T1T by the target value trajectory correction unit 180.

[0089] As shown in FIG. 10, in the prediction horizon, when the gas turbine output GTMW is constant and the turbine bypass valve flow rate Gtb is constant (when the TB ratio is constant), the target value trajectory of the turbine inlet temperature T1T is also constant.

[0090] In contrast to this, as shown in FIG. 11 , in the case where the TB ratio indicated by the future prediction of the turbine bypass valve flow rate (future prediction of the TB ratio in the illustrated example) increases over time in the prediction horizon, the target value trajectory correction unit 180 may perform a correction by adding a positive value to the value of the target value trajectory of the gas turbine output GTMW acquired by the target value trajectory acquisition unit 160, and may perform a correction by adding a negative value to the value of the target value trajectory of the turbine inlet temperature T1T acquired by the target value trajectory acquisition unit 160.

[0091] In this case, the target value trajectory correction unit 180 may correct the target value trajectory of the gas turbine output GTMW based on, for example, the following equation (20). GTMWb(k+i)=GTMWa(k+i)-K_GTMW_Gtb×((Gtb(k+i)-Gtb(k))...(20)

[0092] Here, GTMWb(k+i) is the target value of the gas turbine output GTMW at a time advanced by any control period i from the current time k, and is the target value corrected by the target value trajectory correction unit 180. GTMWe(k+i) is the target value of the gas turbine output GTMW at a time advanced by any control period i from the current time k, and is the target value before correction by the target value trajectory correction unit 180 (i.e., the target value of the gas turbine output GTMW acquired by the target value trajectory acquisition unit 160). K_GTMW_Gtb is the above-mentioned sensitivity coefficient stored in the storage unit 130, and is a coefficient indicating the sensitivity of the gas turbine output to the change in the turbine bypass valve flow rate. Gtb(k+i) is the turbine bypass valve flow rate at a time advanced by any control period i from the current time k in the future turbine bypass valve flow rate prediction. Gtb(k) is the turbine bypass valve flow rate at the current time k.

[0093] Furthermore, the target value trajectory correction unit 180 may correct the target value trajectory of the turbine inlet temperature T1T based on, for example, the following equation (21). T1Tb(k+i)=T1Ta(k+i)-K_T1T_Gtb×((Gtb(k+i)-Gtb(k)) ···(21)

[0094] Here, T1Tb(k+i) is the target value of the turbine inlet temperature T1T at a time advanced by any control period i from the current time k, and is the target value corrected by the target value trajectory correction unit 180. T1Ta(k+i) is the target value of the turbine inlet temperature T1T at a time advanced by any control period i from the current time k, and is the target value corrected by the target value trajectory correction unit 180. K_T1T_Gtb is the above-mentioned sensitivity coefficient stored in the storage unit 130, and is a coefficient indicating the sensitivity of the gas turbine output to the change in the turbine bypass valve flow rate. Gtb(k+i) is the turbine bypass valve flow rate at a time advanced by any control period i from the current time k in the future turbine bypass valve flow rate prediction. Gtb(k) is the turbine bypass valve flow rate at the current time k.

[0095] In the embodiment described with reference to FIGS. 9 to 11 , the optimization unit 120 may optimize the control inputs CSO and IGVCSO based on an evaluation function V (see equation (16)) that depends on the deviation from the corresponding parameter (i.e., the deviation between the target value trajectory after the correction of the gas turbine output and the predicted value of the gas turbine output in the prediction horizon, and the deviation between the target value trajectory after the correction of the turbine inlet temperature and the predicted value of the turbine inlet temperature in the prediction horizon) for each of the target value trajectories after the correction by the target value trajectory correction unit 180 (i.e., each of the target value trajectory after the correction of the gas turbine output and the target value trajectory after the correction of the turbine inlet temperature).

[0096] 9 to 11, even if the bleed valve flow rate Gai or the turbine bypass valve flow rate Gtb changes, the future state quantities of the gas turbine 2 can be accurately predicted using a prediction model using coefficients according to the bleed valve flow rate Gai and the turbine bypass valve flow rate Gtb, making it possible to realize high control performance. Moreover, even if the variation of the AICSO or TBCSO is relatively large, high control performance can be expected by correcting the target value trajectory of the gas turbine output and the target value trajectory of the turbine inlet temperature based on the future bleed valve flow rate prediction and the future turbine bypass valve flow rate prediction.

[0097] The present disclosure is not limited to the above-described embodiments, and includes modifications to the above-described embodiments and appropriate combinations of these modifications. 1 is provided with the extraction piping 15 and the turbine bypass piping 18, the gas turbine 2 may be provided with only one of the extraction piping 15 and the turbine bypass piping 18. In other words, the gas turbine 2 may be provided with only the extraction piping 15 or only the turbine bypass piping 18 out of the extraction piping 15 and the turbine bypass piping 18.

[0098] For example, in the embodiment described with reference to Figs. 9 to 11, the target value trajectory correction unit 180 corrects each of the target value trajectory of the gas turbine output GTMW and the target value trajectory of the turbine inlet temperature T1T based on the future prediction of the bleed valve flow rate and the future prediction of the turbine bypass valve flow rate, but the target value trajectory correction unit 180 may correct the target value trajectory for only one of the gas turbine output GTMW and the turbine inlet temperature T1T based on the future prediction of the bleed valve flow rate and the future prediction of the turbine bypass valve flow rate. That is, the target value trajectory correction unit 180 may correct the target value trajectory for only the gas turbine output GTMW out of the gas turbine output GTMW and the turbine inlet temperature T1T, or may correct the target value trajectory for only the turbine inlet temperature T1T, based on the future prediction of the bleed valve flow rate and the future prediction of the turbine bypass valve flow rate. In this case, the target value trajectory acquisition unit 160 may acquire only the target value trajectory to be corrected out of the target value trajectory of the gas turbine output GTMW and the target value trajectory of the turbine inlet temperature T1T.

[0099] The contents described in each of the above embodiments can be understood, for example, as follows.

[0100] [1] A gas turbine control device according to at least one embodiment of the present disclosure (e.g., the control device 100 described above) includes: A gas turbine control device for controlling a gas turbine (e.g., the above-mentioned gas turbine 2) including a compressor (e.g., the above-mentioned compressor 3), a combustor (e.g., the above-mentioned combustor 4), and a turbine (e.g., the above-mentioned turbine 5), a prediction unit (e.g., the prediction unit 110 described above) configured to predict future state quantities of the gas turbine corresponding to control inputs to the gas turbine over a prediction horizon using a prediction model; an optimization unit (e.g., the optimization unit 120 described above) configured to optimize the control input over at least a portion of the prediction horizon using a prediction result from the prediction unit; A storage unit (e.g., the storage unit 130 described above) for storing sensitivity information indicating the sensitivity of the control input to the rate of change of the state quantity for each of the extracted air flow rates, where the flow rate of the compressed air extracted from the compressor is defined as an extracted air flow rate; an update unit (e.g., the update unit 140 described above) configured to read out the sensitivity information corresponding to the bleed air flow rate assumed in the prediction horizon from the storage unit, and update one or more coefficients of a prediction equation for the state quantity used in the prediction model based on the sensitivity information corresponding to the bleed air flow rate; Equipped with.

[0101] According to the gas turbine control device described in [1] above, one or more coefficients of the prediction formula can be updated according to the bleed air flow rate from the compressor. Therefore, even if the bleed air flow rate from the compressor changes, the future state quantities of the gas turbine can be accurately predicted using a prediction model using coefficients according to the bleed air flow rate, making it possible to realize high control performance.

[0102] [2] In some embodiments, in the gas turbine control device according to [1] above, The gas turbine comprises: A bleed pipe (for example, the above-mentioned bleed pipe 15) that returns a part of the compressed air generated by the compressor to the upstream side of the compressor; A bleed valve (for example, the above-mentioned bleed valve 16) provided in the bleed pipe; Further comprising: the bleed flow rate includes a bleed valve flow rate which is a flow rate of compressed air bled from the compressor by the bleed piping, The update unit is configured to update one or more coefficients of a prediction equation for the state quantity used in the prediction model, based on the sensitivity information corresponding to the bleed valve flow rate.

[0103] According to the gas turbine control device described in [2] above, one or more coefficients of the prediction formula can be updated according to the bleed valve flow rate (i.e., the flow rate of compressed air compressed by the compressor and returned to the upstream side of the compressor). Therefore, even if the bleed valve flow rate changes, the future state quantity of the gas turbine can be predicted with high accuracy, and high control performance can be realized.

[0104] [3] In some embodiments, in the gas turbine control device according to [2] above, The gas turbine further includes a target value trajectory acquisition unit (e.g., the target value trajectory acquisition unit 160) configured to acquire a target value trajectory of an output of the gas turbine, the control input includes an air bleed valve opening control command for controlling an opening of the air bleed valve; The optimization unit is configured to optimize the control input based on an evaluation function that depends on a deviation between an output of the gas turbine and the target value trajectory.

[0105] According to the gas turbine control device described in [3] above, by using the bleed valve position command as one of the control inputs for bringing the gas turbine output closer to the target value trajectory, it is possible to calculate the optimal control input for bringing the gas turbine output closer to the target value trajectory when the bleed valve is operating, and high controllability can be realized.

[0106] [4] In some embodiments, in the gas turbine control device according to any one of [3] above, the target value trajectory acquisition unit is configured to acquire target value trajectories of an output of the gas turbine, a turbine inlet temperature, and an intake temperature of the compressor, the state quantities include information indicating an output of the gas turbine, information indicating the turbine inlet temperature, information indicating an intake temperature of the compressor, and information indicating the extraction valve flow rate, the control input includes a fuel flow rate command for controlling a flow rate of fuel supplied to the combustor, an inlet guide vane opening control command for controlling an opening of an inlet guide vane of the compressor, and an air bleed valve opening control command for controlling an opening of the air bleed valve, The optimization unit is configured to optimize the control input based on an evaluation function that depends on deviations of each of the gas turbine output, the turbine inlet temperature, and the compressor intake temperature from the target value trajectory.

[0107] According to the gas turbine control device described in [4] above, when controlling the opening of the bleed valve, the fuel flow rate command, the air flow rate command, and the bleed valve opening command can be optimized by taking into consideration not only the gas turbine output and the turbine inlet temperature but also the target value trajectory of the compressor intake air temperature, thereby improving the control accuracy of the intake air temperature.

[0108] [5] In some embodiments, in the gas turbine control device according to [2] above, The state quantity does not include information indicating the bleed valve flow rate, The gas turbine control device is configured to determine a bleed valve opening control command for controlling the opening of the bleed valve by a means different from the prediction model.

[0109] According to the gas turbine control device described in [5] above, even if the bleed air flow rate from the compressor changes, the future state quantities of the gas turbine can be accurately predicted using a prediction model using a coefficient according to the bleed air flow rate, and high control performance can be realized. In particular, when the movement of the bleed valve opening control command is relatively small compared to the movements of other control commands (e.g., a fuel flow command and an inlet guide vane opening control command), the behavior of the other control commands within the prediction horizon can be ignored, so high control performance can be expected.

[0110] [6] In some embodiments, in the gas turbine control device according to any one of [2] or [5] above, a target value trajectory acquisition unit (for example, the above-mentioned target value trajectory acquisition unit 160) configured to acquire a target value trajectory of at least one of the parameters of the output of the gas turbine and the turbine inlet temperature among the state quantities; a target value trajectory correction unit (for example, the above-mentioned target value trajectory correction unit 180) configured to obtain a bleed valve flow rate future prediction indicating a future predicted value of the bleed valve flow rate, and to correct a target value trajectory of the at least one parameter based on the bleed valve flow rate future prediction; Further comprising: The optimization unit is configured to optimize the control input based on an evaluation function that depends on a deviation between the target value trajectory corrected by the target value trajectory correction unit and the parameter.

[0111] According to the gas turbine control device described in [6] above, even if the fluctuation of the extraction valve position control command is relatively large, high control performance can be expected by correcting the target value trajectory of at least one of the parameters, the gas turbine output and the turbine inlet temperature, based on the future prediction of the extraction valve flow rate.

[0112] [7] In some embodiments, in the gas turbine control device according to [1] or [2] above, The gas turbine comprises: a turbine bypass pipe (for example, the above-mentioned turbine bypass pipe 18) that bypasses the turbine and supplies a portion of the compressed air generated by the compressor to an exhaust side of the turbine; A turbine bypass valve (for example, the above-mentioned turbine bypass valve 22) provided in the turbine bypass piping; Further comprising: the bleed flow rate includes a turbine bypass valve flow rate which is a flow rate of compressed air bled from the compressor by the turbine bypass piping, The update unit is configured to update one or more coefficients of a prediction equation for the state quantity used in the prediction model, based on the sensitivity information corresponding to the turbine bypass valve flow rate.

[0113] According to the gas turbine control device described in [7] above, one or more coefficients of the prediction formula can be updated according to the turbine bypass valve flow rate (i.e., the flow rate of compressed air that bypasses the combustor and turbine from the compressor and is supplied to the exhaust side of the turbine). Therefore, even if the turbine bypass valve flow rate changes, the future state quantities of the gas turbine can be predicted with high accuracy, making it possible to achieve high control performance.

[0114] [8] In some embodiments, in the gas turbine control device according to [7] above, The gas turbine further includes a target value trajectory acquisition unit (e.g., the target value trajectory acquisition unit 160) configured to acquire a target value trajectory of an output of the gas turbine, the control input includes a turbine bypass valve opening control command for controlling an opening of the turbine bypass valve, The optimization unit is configured to optimize the control input based on an evaluation function that depends on a deviation between an output of the gas turbine and the target value trajectory.

[0115] According to the gas turbine control device described in [8] above, by using a turbine bypass valve opening command as one of the control inputs for bringing the gas turbine output closer to the target value trajectory, it is possible to calculate an optimal control input for bringing the gas turbine output closer to the target value trajectory when the turbine bypass valve is operating, and high controllability can be realized.

[0116] [9] In some embodiments, in the gas turbine control device according to [8] above, the target value trajectory acquisition unit is configured to acquire target value trajectories of each of an output and a turbine inlet temperature of the gas turbine, the state quantities include information indicating an output of the gas turbine, information indicating the turbine inlet temperature, and information indicating the turbine bypass valve flow rate, the control input includes a fuel flow rate command for controlling a flow rate of fuel supplied to the combustor, an inlet guide vane opening control command for controlling an opening of an inlet guide vane of the compressor, and a turbine bypass valve opening control command for controlling an opening of the turbine bypass valve, The optimization unit is configured to optimize the control input based on an evaluation function that depends on a deviation between each of an output power and a turbine inlet temperature of the gas turbine and the target value trajectory.

[0117] According to the gas turbine control device described in [9] above, when controlling the aperture of the turbine bypass valve, the fuel flow rate command, the air flow rate command, and the turbine bypass valve aperture command can be optimized in consideration of the target value trajectories of the gas turbine output and the turbine inlet temperature, respectively, so that the control accuracy of the turbine inlet temperature can be improved.

[0118]

[10] In some embodiments, in the gas turbine control device according to [7] above, the state quantity does not include information indicating the turbine bypass valve flow rate, The gas turbine control device is configured to determine a turbine bypass valve opening control command for controlling the opening of the turbine bypass valve by a means different from the prediction model.

[0119] According to the gas turbine control device described in

[10] above, even if the bleed air flow rate from the compressor changes, the future state quantities of the gas turbine can be accurately predicted using a prediction model using a coefficient according to the bleed air flow rate, and high control performance can be realized. In particular, when the movement of the turbine bypass valve opening control command is relatively small compared to the movements of other control commands (e.g., the fuel flow command and the inlet guide vane opening control command), the behavior of the other control commands within the prediction horizon can be ignored, so high control performance can be expected.

[0120]

[11] In some embodiments, in the gas turbine control device according to [7] or

[10] above, a target value trajectory acquisition unit (for example, the above-mentioned target value trajectory acquisition unit 160) configured to acquire a target value trajectory of at least one parameter of the gas turbine output and the turbine inlet temperature among the state quantities; a target value trajectory correction unit (for example, the above-mentioned target value trajectory correction unit 180) configured to obtain a turbine bypass valve flow rate future prediction indicating a future predicted value related to the turbine bypass valve flow rate, and to correct a target value trajectory of the at least one parameter based on the turbine bypass valve flow rate future prediction; Further comprising: The optimization unit is configured to optimize the control input based on an evaluation function that depends on a deviation between the target value trajectory corrected by the target value trajectory correction unit and the parameter.

[0121] According to the gas turbine control device described in

[11] above, even if the fluctuation of the turbine bypass valve position control command is relatively large, high control performance can be expected by correcting the target value trajectory of at least one of the parameters, the gas turbine output and the turbine inlet temperature, based on the future prediction of the turbine bypass valve flow rate.

[0122]

[12] In some embodiments, in the gas turbine control device according to any one of [1] to

[11] above, The process of predicting the state quantity, the process of updating the one or more coefficients, and the process of optimizing the control input are configured to be executed for each control period.

[0123] According to the gas turbine control device described in

[11] above, the control input is reviewed for each control period, thereby improving prediction accuracy.

[0124]

[13] A method for determining control inputs for a gas turbine according to at least one embodiment of the present disclosure includes: a prediction step in which a prediction model predicts future state quantities of a gas turbine (e.g., the above-mentioned gas turbine 2) corresponding to control inputs to the gas turbine in a prediction horizon; an optimization step of optimizing the control input in at least a part of the prediction horizon using a prediction result in the prediction step; When a flow rate of compressed air extracted from a compressor of the gas turbine is defined as an extracted flow rate, a storage unit (e.g., the above-mentioned storage unit 130) stores sensitivity information indicating a sensitivity of the control input to a change rate of the state quantity for each of the extracted flow rates; reading out from the storage unit the sensitivity information corresponding to the bleed air flow rate assumed in the prediction horizon, and updating one or more coefficients of a prediction equation for the state quantity used in the prediction model based on the sensitivity information corresponding to the bleed air flow rate; Equipped with.

[0125] According to the gas turbine control input determination method described in

[13] above, one or more coefficients of the prediction formula can be updated according to the bleed air flow rate from the compressor. Therefore, even if the bleed air flow rate from the compressor changes, the future state variables of the gas turbine can be accurately predicted using a prediction model using coefficients according to the bleed air flow rate, making it possible to achieve high control performance.

[0126]

[14] A gas turbine control input determination program according to at least one embodiment of the present disclosure, On the computer, predicting future state quantities of the gas turbine corresponding to control inputs to the gas turbine over a prediction horizon using a prediction model; a step of optimizing the control input over at least a portion of the prediction horizon using a prediction result from the step of predicting; a step of storing, in a storage unit (e.g., the storage unit 130 described above), sensitivity information indicating the sensitivity of the control input to the rate of change of the state quantity for each of the extracted air flow rates, where the flow rate of compressed air extracted from a compressor of the gas turbine is defined as an extracted air flow rate; and a step of reading out from the storage unit the sensitivity information corresponding to the bleed air flow rate assumed in the prediction horizon, and updating one or more coefficients of a prediction equation for the state quantity used in the prediction model based on the sensitivity information corresponding to the bleed air flow rate; Execute the command.

[0127] According to the gas turbine control input determination program described in

[14] above, one or more coefficients of the prediction formula can be updated according to the bleed air flow rate from the compressor. Therefore, even if the bleed air flow rate from the compressor changes, the future state variables of the gas turbine can be accurately predicted using a prediction model using coefficients according to the bleed air flow rate, making it possible to achieve high control performance. [Explanation of symbols]

[0128] 1 Gas Turbine System 2. Gas turbine 3. Compressor 4. Combustor 5. Turbine 7. Generator 8A, 8B Rotating shaft 9 Inlet guide vane 10. Intake equipment 12 Fuel supply line 13,17,24,30 Flow meter 14 Fuel flow control valve 15 Extraction piping 16 Bleed valve 18 Turbine bypass piping 20 Turbine exhaust piping 22 Turbine bypass valve 26,28,32,36 thermometer 34 Output meter 35 Actuator 72 processors 74 RAM 76 ROM 78 HDD 80 Input I / F 82 Output I / F 84 Bus 100 Control device 110 Prediction Department 120 Optimization Department 130 Storage section 140 Update Department 150 State quantity acquisition unit 160 Target value trajectory acquisition unit 170 Table Update Section 180 Target value trajectory correction unit

Claims

1. A gas turbine control device for controlling a gas turbine including a compressor, a combustor, and a turbine, comprising: a prediction unit configured to predict a future state quantity of the gas turbine corresponding to a control input to the gas turbine in a prediction horizon using a prediction model; an optimization unit configured to optimize the control input over at least a portion of the prediction horizon using a prediction result from the prediction unit; a storage unit for storing sensitivity information indicating a sensitivity of the control input to a rate of change of the state quantity for each of the extracted air flow rates, where the flow rate of the compressed air extracted from the compressor is defined as an extracted air flow rate; an update unit configured to read out from the storage unit the sensitivity information corresponding to the bleed air flow rate assumed in the prediction horizon, and to update one or more coefficients of a prediction equation for the state quantity used in the prediction model, based on the sensitivity information corresponding to the bleed air flow rate; A gas turbine control device comprising:

2. The gas turbine comprises: a bleed pipe for returning a portion of the compressed air generated by the compressor to the upstream side of the compressor; A bleed valve provided in the bleed pipe; Further comprising: the bleed flow rate includes a bleed valve flow rate which is a flow rate of compressed air bled from the compressor by the bleed piping, 2. The gas turbine control device according to claim 1, wherein the update unit is configured to update one or more coefficients of a prediction equation for the state quantity used in the prediction model, based on the sensitivity information corresponding to the extraction valve flow rate.

3. a target value trajectory acquisition unit configured to acquire a target value trajectory of an output of the gas turbine, the control input includes an air bleed valve opening control command for controlling an opening of the air bleed valve; 3. The gas turbine control device according to claim 2, wherein the optimization unit is configured to optimize the control input based on an evaluation function that depends on a deviation between an output of the gas turbine and the target value trajectory.

4. the target value trajectory acquisition unit is configured to acquire target value trajectories of an output of the gas turbine, a turbine inlet temperature, and an intake temperature of the compressor, the state quantities include information indicating an output of the gas turbine, information indicating the turbine inlet temperature, information indicating an intake temperature of the compressor, and information indicating the extraction valve flow rate, the control input includes a fuel flow rate command for controlling a flow rate of fuel supplied to the combustor, an inlet guide vane opening control command for controlling an opening of an inlet guide vane of the compressor, and an air bleed valve opening control command for controlling an opening of the air bleed valve, 4. The gas turbine control device according to claim 3, wherein the optimization unit is configured to optimize the control input based on an evaluation function that depends on a deviation between an output of the gas turbine, a turbine inlet temperature, and an intake temperature of the compressor and the target value trajectory.

5. The state quantity does not include information indicating the bleed valve flow rate, 3. The gas turbine control device according to claim 2, wherein the gas turbine control device is configured to determine an air bleed valve opening control command for controlling an opening of the air bleed valve by a means other than the prediction model.

6. a target value trajectory acquisition unit configured to acquire a target value trajectory of at least one parameter of an output of the gas turbine and a turbine inlet temperature among the state quantities; a target value trajectory correction unit configured to obtain a bleed valve flow rate future prediction indicating a future predicted value of the bleed valve flow rate, and to correct a target value trajectory of the at least one parameter based on the bleed valve flow rate future prediction; Further comprising:

6. The gas turbine control device according to claim 2, wherein the optimization unit is configured to optimize the control input based on an evaluation function that depends on a deviation between the target value trajectory corrected by the target value trajectory correction unit and the parameter.

7. The gas turbine comprises: a turbine bypass pipe that bypasses the turbine and supplies a portion of the compressed air generated by the compressor to an exhaust side of the turbine; a turbine bypass valve provided in the turbine bypass piping; Further comprising: The extraction flow rate includes a turbine bypass valve flow rate, which is a flow rate of compressed air extracted from the compressor by the turbine bypass piping.

2. The gas turbine control device according to claim 1, wherein the update unit is configured to update one or more coefficients of a prediction equation of the state quantity used in the prediction model, based on the sensitivity information corresponding to the turbine bypass valve flow rate.

8. a target value trajectory acquisition unit configured to acquire a target value trajectory of an output of the gas turbine, the control input includes a turbine bypass valve opening control command for controlling an opening of the turbine bypass valve, The gas turbine control device according to claim 7 , wherein the optimization unit is configured to optimize the control input based on an evaluation function that depends on a deviation between an output of the gas turbine and the target value trajectory.

9. the target value trajectory acquisition unit is configured to acquire target value trajectories of each of an output and a turbine inlet temperature of the gas turbine, the state quantities include information indicating an output of the gas turbine, information indicating the turbine inlet temperature, and information indicating the turbine bypass valve flow rate, the control input includes a fuel flow rate command for controlling a flow rate of fuel supplied to the combustor, an inlet guide vane opening control command for controlling an opening of an inlet guide vane of the compressor, and a turbine bypass valve opening control command for controlling an opening of the turbine bypass valve, 9. The gas turbine control device according to claim 8, wherein the optimization unit is configured to optimize the control input based on an evaluation function that depends on a deviation between an output of the gas turbine and a turbine inlet temperature of the gas turbine and the target value trajectory.

10. the state quantity does not include information indicating the turbine bypass valve flow rate, 8. The gas turbine control device according to claim 7, wherein the gas turbine control device is configured to determine a turbine bypass valve opening control command for controlling an opening of the turbine bypass valve by a means other than the prediction model.

11. a target value trajectory acquisition unit configured to acquire a target value trajectory of at least one parameter of an output of the gas turbine and a turbine inlet temperature among the state quantities; a target value trajectory correction unit configured to acquire a turbine bypass valve flow rate future prediction indicating a future predicted value related to the turbine bypass valve flow rate, and to correct a target value trajectory of the at least one parameter based on the turbine bypass valve flow rate future prediction; Further comprising: the optimization unit is configured to optimize the control input based on an evaluation function that depends on a deviation between the target value trajectory corrected by the target value trajectory correction unit and the parameter. The gas turbine control device according to claim 7 or 10.

12. 2. The gas turbine control device according to claim 1, configured to execute the process of predicting the state quantity, the process of updating the one or more coefficients, and the process of optimizing the control input for each control period.

13. a prediction step in which a prediction model predicts future state quantities of the gas turbine corresponding to control inputs to the gas turbine in a prediction horizon; an optimization step of optimizing the control input in at least a part of the prediction horizon using a prediction result in the prediction step; a step of storing, in a storage unit, sensitivity information indicating a sensitivity of the control input to a rate of change of the state quantity for each of the extracted air flow rates, where the flow rate of compressed air extracted from a compressor of the gas turbine is defined as an extracted air flow rate; reading out from the storage unit the sensitivity information corresponding to the bleed air flow rate assumed in the prediction horizon, and updating one or more coefficients of a prediction equation for the state quantity used in the prediction model based on the sensitivity information corresponding to the bleed air flow rate; A method for determining control inputs for a gas turbine comprising:

14. On the computer, predicting future state quantities of the gas turbine corresponding to control inputs to the gas turbine over a prediction horizon using a prediction model; a step of optimizing the control input over at least a portion of the prediction horizon using a prediction result from the step of predicting; a step of storing, in a storage unit, sensitivity information indicating a sensitivity of the control input to a rate of change of the state quantity for each of the extracted air flow rates, where the flow rate of compressed air extracted from a compressor of the gas turbine is defined as an extracted air flow rate; and a step of reading out from the storage unit the sensitivity information corresponding to the bleed air flow rate assumed in the prediction horizon, and updating one or more coefficients of a prediction equation for the state quantity used in the prediction model based on the sensitivity information corresponding to the bleed air flow rate; A gas turbine control input determination program for executing the above.

Images