WO2021220901A1 - 制御装置 - Google Patents

制御装置 Download PDF

Info

Publication number
WO2021220901A1
WO2021220901A1 PCT/JP2021/016135 JP2021016135W WO2021220901A1 WO 2021220901 A1 WO2021220901 A1 WO 2021220901A1 JP 2021016135 W JP2021016135 W JP 2021016135W WO 2021220901 A1 WO2021220901 A1 WO 2021220901A1
Authority
WO
WIPO (PCT)
Prior art keywords
amount
value
calculation unit
control
specified
Prior art date
Application number
PCT/JP2021/016135
Other languages
English (en)
French (fr)
Inventor
寿宏 馬場
Original Assignee
株式会社Mhiパワーコントロールシステムズ
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 株式会社Mhiパワーコントロールシステムズ filed Critical 株式会社Mhiパワーコントロールシステムズ
Publication of WO2021220901A1 publication Critical patent/WO2021220901A1/ja

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B11/00Automatic controllers
    • G05B11/01Automatic controllers electric
    • G05B11/32Automatic controllers electric with inputs from more than one sensing element; with outputs to more than one correcting element

Definitions

  • the present disclosure relates to the control of the state quantity (controlled quantity) of the controlled object in the plant, and more particularly to the determination of the manipulated variable at the operating end capable of controlling the controlled quantity.
  • state quantities such as main steam temperature, main steam pressure, reheated steam temperature, and NOx concentration at the outlet of the denitration device are measured so that the measured value of the state quantity becomes the set value.
  • Feedback control is performed on. Further, in such feedback control, the response to a sudden change in the state quantity due to the occurrence of disturbance may be delayed, and the response may deviate from the permissible fluctuation range of the controlled quantity. Therefore, when a disturbance occurs, the preceding amount for canceling the influence before appearing in the control amount is calculated, and the operation amount added to the operation amount at the operation end calculated by the feedback control is given to the operation end, feedforward. Control may also be performed at the same time.
  • the above-mentioned disturbances to the temperature and pressure of the main steam, reheated steam temperature, NOx concentration, etc. are caused by various heat exchanges such as the furnace wall evaporation pipe forming the furnace wall and the superheater and reheater installed in the flue. It is caused by the operation of the steam blower to remove the ash adhering to the surface of the vessel, and events such as ignition and extinguishing of the burner. Specifically, when the soot blower for the furnace wall evaporation pipe is operated, the temperature such as the main steam temperature and the reheat steam temperature is lowered, and when the soot blower for the heat exchanger installed in the flue is operated, These temperatures rise. Therefore, for example, control is performed to reduce or increase the spray amount in the temperature reducer for adjusting the main steam temperature within a certain range in advance (see Patent Document 1).
  • At least one embodiment of the present invention aims to provide a control device capable of more appropriately determining a preceding amount for feedforward control of a controlled amount of a plant.
  • the control device is It is a control device that determines the operation amount of the target operation end that can control the control amount of the plant at each specified timing. Using a learning model that learns the relationship between the process parameter for specifying the operating state of the plant and the fluctuation value of the control amount, the predicted fluctuation value of the control amount according to the input of the process parameter is output. Fluctuation prediction unit configured as With a preceding amount calculation unit configured to calculate a preceding amount for correcting the first operation amount calculated according to the deviation between the measured value of the controlled amount and the control set value based on the predicted fluctuation value.
  • the preceding amount calculation unit When an event that changes the control amount occurs, a first calculation unit configured to calculate a specified preceding amount according to the duration of the event, and a first calculation unit. A second calculation unit configured to calculate a leading index for correcting the specified leading amount based on the predicted fluctuation value and the deviation. It has a third calculation unit configured to calculate the preceding amount based on the specified preceding amount and the preceding index.
  • a control device capable of more appropriately determining a preceding amount for feedforward control of a controlled amount of a plant.
  • expressions such as “same”, “equal”, and “homogeneous” that indicate that things are in the same state not only represent exactly the same state, but also have tolerances or differences to the extent that the same function can be obtained. It shall also represent the existing state.
  • an expression representing a shape such as a quadrangular shape or a cylindrical shape not only represents a shape such as a quadrangular shape or a cylindrical shape in a geometrically strict sense, but also an uneven portion or chamfering within a range in which the same effect can be obtained.
  • the shape including the part and the like shall also be represented.
  • the expressions “equipped”, “equipped”, “equipped”, “included”, or “have” one component are not exclusive expressions that exclude the existence of other components.
  • FIG. 1 is a diagram schematically showing a configuration of a boiler 71 included in a plant 7 according to an embodiment of the present invention.
  • the control device 1 periodically sets the operation amount MV (command value) of the operation end (hereinafter, target operation end 9) capable of controlling the state amount (hereinafter, control amount) of the plant 7 to be controlled by the own device. It is a device that can be determined at each specified timing such as.
  • the plant 7 has a plurality of operating ends capable of changing movable parts such as a valve body (valve opening) according to, for example, an operating amount MV given from the outside, and various types such as temperature, pressure, and flow rate. Multiple sensors are installed to measure the state quantity.
  • the operation amount MV of the target operation end 9, which is one or more of the plurality of operation ends installed in the plant 7, is controlled by the measured value PV of the control amount obtained by using the above sensor. It is determined by feedback control and feedforward control so as to be SV (target value). Further, the control device 1 can control the target operation end 9 by outputting (transmitting) the determined operation amount MV to the outside.
  • the plant 7 shown in FIG. 1 is a power plant equipped with a boiler 71 such as a boiler / turbine power generation facility (BTG), for example.
  • the power generation plant supplies the superheated steam (main steam) generated by the boiler 71 to, for example, the steam turbine T and rotationally drives the rotating shaft thereof, thereby rotating and driving the generator G connected to the rotating shaft of the steam turbine T.
  • the control device 1 may use at least one of state quantities such as the main steam temperature, the main steam pressure, the reheated steam temperature, and the NOx concentration at the outlet of the denitration device as the control amount. ..
  • the target operating ends 9 for controlling these control amounts are the spray valve 7v, the fuel supply device 72v, and the gas recirculation fan inlet damper 87 or burner installed at the inlets of the gas recirculation fan 86, respectively, as described later.
  • An angle adjuster (not shown), an ammonia flow rate control valve 81v, or the like.
  • the boiler 71 exchanges heat generated by burning pulverized coal fuel obtained by crushing coal, for example, inside the furnace 71f with water supply or steam, and supplies the heat to the steam turbine T or the like. It is a pulverized coal-fired boiler configured to generate superheated steam (main steam).
  • the boiler 71 may be a waste heat recovery boiler such as a soda recovery boiler having a different fuel or combustion method, a fluidized bed, a fluidized bed boiler, or the like.
  • the boiler 71 will be described in detail by taking the pulverized coal-cooked boiler shown in FIG. 1 as an example.
  • the boiler 71 (boiler main body) has a furnace 71f, a combustion device 72 installed in the furnace 71f, and a flue 73 connected to the boiler 71.
  • the furnace 71f has a combustion chamber inside, and is installed along the vertical direction, for example, in the shape of a hollow square cylinder (see FIG. 1).
  • the inner wall surface of the wall of the furnace 71f (fire furnace wall) forming the combustion chamber is composed of a furnace wall evaporation tube 71p (heat transfer tube) for generating steam from water and fins connecting the furnace wall evaporation tube 71p.
  • a plurality of furnace wall evaporation pipes 71p arranged in the vertical direction are arranged side by side in the horizontal direction, and the adjacent furnace wall evaporation pipes 71p and the furnace wall in these plurality of furnace wall evaporation pipes 71p are arranged.
  • the fins are arranged so as to close the space between the evaporation pipe 71p and the evaporation pipe 71p.
  • the furnace 71f is provided with an inclined surface at the bottom of the furnace, and a furnace wall evaporation pipe 71p (furnace bottom evaporation pipe) is also provided on the inclined surface to form a bottom surface.
  • an ash removing device such as a soot blower (not shown) configured to remove deposits such as ash and unburned components adhering to the surface of the furnace wall evaporation pipe 71p is provided inside the combustion chamber.
  • the combustion device 72 is a device for supplying fuel to the inside of the above-mentioned furnace 71f for combustion, and is used to supply one or more combustion burners 74 installed on the furnace wall, a wind box 77e, and a boiler 71.
  • a fuel supply device 72v for adjusting the amount is provided. Specifically, it is configured to burn fuel by supplying fuel to the inside of the furnace 71f via a combustion burner 74 and also supplying combustion air (secondary air) through a wind box 77e. NS.
  • the amount of fuel supplied to the combustion burner 74 is adjusted by controlling the fuel supply device 72v (operating end) so that the measured value PV (control amount) of the main steam pressure becomes the control set value SV. .. In the embodiment shown in FIG.
  • the steam generated by the furnace wall evaporation pipe 71p is heated by the boiler 79a when flowing through the steam pipe 7p and then sent to the steam turbine T, but the main steam pressure is the steam pipe 7p. It is detected by the pressure sensor 7s installed at the outlet of the boiler in the above.
  • the combustion device 72 may further include an AA port 77p for supplying additional air (AA) installed on the furnace wall.
  • the combustion air is provided by the transport air (primary air) described later for supplying fuel to the combustion burner 74 and the combustion air (secondary air) supplied to the furnace 71f via the air box 77e.
  • the supply amount is less than the theoretical air amount with respect to the supply amount of the pulverized coal fuel, the inside is maintained in a reducing atmosphere, and after the pulverized coal fuel is reduced and burned, new combustion air is used.
  • (AA) is supplied to complete the oxidative combustion of the pulverized coal fuel, and the amount of NOx generated on the combustion burner 74 side can be reduced (two-stage combustion method).
  • the combustion device 72 includes a plurality of combustion burners 74. Then, one or more (plurality in FIG. 1) combustion burners 74 are installed on the furnace wall in the circumferential direction or opposite positions of the furnace 71f at equal intervals, for example, and are installed in such a circumferential direction or opposite positions. The one or more combustion burners 74 are installed in one or more stages (five stages of 74a to 74e in FIG. 1) along the vertical direction. Each of the plurality of combustion burners 74 is connected to one of a plurality of crushers 76 (76a to 76e in FIG.
  • the pulverized coal produced in (1) is supplied to the combustion burner 74 by the transport air (primary air) via the pulverized coal supply pipe 75 in which the fuel supply device 72v is installed.
  • the crusher 76 is designed to produce pulverized coal by pulverizing coal supplied by a raw coal transport system (not shown) to a predetermined size of pulverized coal.
  • the above-mentioned air box 77e is provided at the installation position of each combustion burner 74 in the furnace 71f, and one end of an air duct 77a for guiding air (combustion air) is connected.
  • the other end of the air duct 77a is connected to the blower 78, and the air pushed by the blower 78 is supplied to the air box 77e through the air duct 77a.
  • the branch duct 77c branches from the branch portion 77b provided in the air duct 77a and is connected to the above AA port 77p, so that the AA corresponding to the opening degree of the AA damper 77d installed in the branch duct 77c can be obtained.
  • the air supplied from the blower 78 to the air duct 77a is heated by the air heater 82 (described later), and then is guided to the air box 77e at the branch portion 77b and the AA port via the branch duct 77c. It is designed to be distributed to AA, which leads to 77p.
  • the flue 73 is a duct that guides the combustion gas generated by the combustion of fuel in the furnace 71f, and is connected to the upper part of the furnace 71f in the vertical direction as shown in FIG.
  • the flue 73 has a plurality of heat exchangers 79 for generating steam and heating (overheating) the steam (in FIG. 1, a heater 79a (SH), a reheater 79b (RH), and a coal saver). 79c (ECO)) and an ash removal device (not shown) provided for each heat exchanger 79 to remove deposits (ash, etc.) adhering to the heat transfer surface of the heat exchanger 79 are installed.
  • 79c (ECO)) and an ash removal device not shown
  • the combustion gas flowing from the furnace 71f into the flue 73 heats the water supply and steam flowing through the furnace wall and the heat exchanger 79 to generate main steam
  • the gas-liquid mixed fluid in which water (saturated water) and steam (saturated steam) generated by boiling the water inside the furnace wall evaporation pipe 71p is moist. It flows into a minute (steam water) separator (not shown), and the gas-liquid mixed fluid is separated into steam (gas phase) and water (liquid phase). After that, the steam (main steam) separated from water by the moisture (steam water) separator is sent to the steam turbine T or the like through the steam pipe 7p (superheated steam pipe) forming the steam flow path.
  • the main steam temperature (control amount) is kept constant at the control set value SV by receiving heating by the superheater installed in the steam pipe 7p and cooling as needed by the temperature reducer 7c (superheater). Adjustments are made to.
  • a cooling medium supply pipe 7L for supplying the cooling medium 7w is connected to the temperature reducing device 7c, and the cooling medium supply pipe 7L is supplied to the temperature reducing device 7c by the spray valve 7v (operating end).
  • the flow rate of the cooling medium 7w can be adjusted.
  • the valve opening degree of the spray valve 7v (operating end) is adjusted so that the measured value PV (controlled amount) of the thermometer (boiler outlet temperature measuring sensor 7t) installed in the steam pipe 7p becomes the control set value SV. It is supposed to be done.
  • the gas recirculation fan inlet damper 87 so that the measured value PV of the reheated steam temperature (controlled amount) becomes the control set value SV. (Operating end) and a burner angle adjuster (operating end) (not shown) that can adjust the burner nozzle angle are controlled.
  • the exhaust gas treatment device is composed of an exhaust gas passage 8 connected to the flue 73 and a plurality of devices installed in the exhaust gas passage 8.
  • An air heater 82 for exchanging heat with the air heater 82, a soot dust treatment device 83, an attractive blower 84, and the like are provided, and a chimney 85 is provided at the downstream end.
  • the denitration device 81 purifies NOx contained in the exhaust gas passing through the catalyst by supplying 8w or the like of ammonia to the catalyst installed inside.
  • a pipe for supplying ammonia 8w (ammonia supply pipe 81p) is connected to the denitration device 81, and by adjusting the valve opening degree of the ammonia flow rate control valve 81v (operating end), The amount of ammonia 8w supplied to the denitration device 81 can be controlled.
  • a sensor (NOx sensor 81s) capable of detecting the NOx concentration is installed at the outlet of the denitration device 81, and the ammonia flow rate control valve 81v so that the measured value PV by the NOx sensor 81s becomes the control set value SV. (Operation end) is controlled.
  • the plant 7 is equipped with a plurality of devices such as a boiler 71, an exhaust gas treatment device, a steam turbine T, an operation end, and a sensor, each of which is composed of a plurality of devices. Then, the control amount is controlled by the control device 1 or the like so that the operation of the plant 7 is optimized while satisfying the output request based on the power demand of the plant 7 (boiler 71), for example.
  • the control device 1 or the like described the control device 1 or the like performs control by feedforward control in addition to feedback control.
  • the feedforward control is performed in order to suppress the influence of the occurrence of an event that causes disturbance in the feedback control of the controlled amount.
  • the operation amount obtained by feedback control (hereinafter, the first operation amount Fb) is determined at a predetermined timing with respect to the spray valve 7v so that the measured value PV of the main steam temperature becomes the control set value SV.
  • the deashing device of the furnace wall evaporation pipe 71p operates (event occurs)
  • the deposits are removed, so that the heat generated by the combustion in the furnace 71f and the fluid (water) in the evaporation pipe 71p are removed.
  • the amount of heat exchanged with and from increases sharply.
  • the first manipulated variable Fb for at least one of the gas recirculation fan inlet damper 87 or the burner angle adjuster is at the specified timing so that the measured value PV of the reheated steam temperature becomes the control set value SV. It is determined.
  • the operation of the deashing device (occurrence of an event) of the furnace wall evaporation pipe 71p becomes the above-mentioned disturbance, and the reheated steam temperature drops sharply.
  • the operation of the deashing device of the reheater 79b becomes the above-mentioned disturbance, and the reheated steam temperature rises sharply.
  • the first operation amount Fb for the fuel supply device 72v is determined at the specified timing so that the measured value PV of the main steam pressure becomes the control set value SV, but the operation of the furnace wall evaporation pipe 71p ( The occurrence of the event) becomes the above-mentioned disturbance, and the main vapor pressure rises sharply.
  • the first manipulated variable Fb with respect to the ammonia flow control valve 81v is determined at a specified timing so that the measured value PV of the NOx concentration becomes the control set value SV.
  • the ignition or extinguishing (occurrence of an event) of the combustion burner 74 becomes the above-mentioned disturbance, and the NOx concentration of the steam rises sharply.
  • the control device 1 calculates (calculates) the preceding amount Ff as described below in order to calculate an appropriate preceding amount Ff according to the operating state of the plant 7 (boiler 71) and the like.
  • FIG. 2 is a block diagram showing a configuration of a control device 1 according to an embodiment of the present invention.
  • FIG. 3 is a diagram showing the configurations of the first calculation unit 31 and the third calculation unit 5 according to the embodiment of the present invention, and is the details of the portion shown by A in FIG.
  • FIG. 4 is a diagram showing an example of the relationship between the specified preceding quantity Qa calculated by the first calculation unit 31 according to the first embodiment of the present invention and the duration t of the event.
  • FIG. 5 is a diagram illustrating the relationship between the specified leading index value Q1 calculated by the specified leading index calculation unit 42 according to the embodiment of the present invention and the predicted fluctuation value V.
  • FIG. 6 is a diagram showing a configuration of an evaluation unit 43 according to an embodiment of the present invention.
  • the control device 1 includes a fluctuation prediction unit 2 and a preceding amount calculation unit 3. Each of these functional parts (circuits) will be described.
  • the control device 1 may be composed of a computer. It includes a processor 11 such as a CPU and a storage unit 12 such as a memory such as a ROM or RAM. Then, the processor 11 operates (data calculation, etc.) according to the instruction of the program (control program 10) loaded in the memory (main storage unit) to realize each of the above functional units.
  • the control program may be realized by an MCU (Micro Controller Unit) or the like, or may be mounted on one integrated circuit (chip) such as a SoC (System on Chip).
  • the above program is software for realizing each functional unit described later in a computer, and may be stored in a storage medium that can be read by the computer.
  • the fluctuation prediction unit 2 uses a learning model M that has learned the relationship between the process parameter P for specifying the operating state of the plant 7 and the fluctuation value of the control amount, and uses the learning model M to learn the control amount according to the input of the process parameter P.
  • This is a functional unit configured to output the predicted fluctuation value V.
  • the process parameter P includes at least one of the measured value of the state quantity and the manipulated variable at the operation end.
  • the measured value of the state quantity can be measured by a sensor (not shown or shown), and the manipulated variable at the operating end can be obtained from the operating end or the control unit at the operating end.
  • the process parameters P are a) main steam pressure, b) exhaust gas temperature at the outlet of the coal saver 79c, c) valve opening of the spray valve 7v, d) fuel coal flow rate (ton / hour), e) water supply flow rate.
  • the learning model M learns the relationship between the process parameter P including the parameter values of one or more parameters described above and the fluctuation value of the controlled variable measured under the process parameter P (machine learning). It is created by doing.
  • This learning may be performed by a well-known machine learning method for learning data (teacher data).
  • the training data is, for example, the measurement result of each parameter value of the process parameter P obtained when the plant 7 (boiler 71) was operated (actual operation, test run) in the past, and the fluctuation value of the control amount when this was obtained. It may be composed of a plurality of data corresponding to (associating) with.
  • the training data may be composed of a plurality of data in which each parameter value of the process parameter P obtained from the literature or numerical analysis (material simulation) is associated with the fluctuation value of the control amount. These data may be combined.
  • the learning model M created in this way outputs the fluctuation value of the control amount obtained from the above relationship as the predicted fluctuation value V with respect to the input process parameter P.
  • the leading amount calculation unit 3 corrects the leading manipulated variable Fb calculated according to the deviation D between the measured value PV of the controlled variable and the control set value SV based on the predicted variable value V described above. It is a functional unit configured to calculate Ff.
  • the first manipulated variable Fb is a feedback manipulated variable calculated by feedback control based on the above deviation D.
  • the value corrected by the preceding quantity Ff is output as the final manipulated variable MV to the target manipulated end 9 operated by this manipulated variable MV (see FIG. 2).
  • the preceding amount calculation unit 3 includes the first calculation unit 31 (see FIG. 3), the second calculation unit 4, and the second calculation unit 4 in order to calculate the preceding amount Ff. It has three functional units (circuits) with a third calculation unit 5 (see FIG. 3).
  • the event is the operation of the deashing device
  • the control amount is the main steam temperature
  • the target operating end 9 is the spray valve 7v will be described as an example of these functional units.
  • the first calculation unit 31 calculates a specified preceding amount Qa according to the duration t of the event when an event such as operation of the deashing device described above, ignition of the combustion burner 74, or fire extinguishing occurs to change the control amount. It is a functional unit (function generator) configured as such. More specifically, as shown in FIG. 3, the operating time of the deashing device (soothing blower in this embodiment) of the superheater 79a is input to the first calculation unit 31 as the duration t of the event. Then, the specified timing is determined by using the relationship such as a function (see FIG. 4) or a table created based on, for example, control theory or experience, which can uniquely calculate the specified preceding quantity Qa from the duration t of the input event. The specified leading amount Qa is calculated according to the duration t of the event that advances each time.
  • the first calculation unit 31 uses a plurality of different functions (three in FIG. 3) with respect to the duration t of the input event.
  • the specified leading amount Qa is calculated respectively.
  • the second calculation unit 4 is based on the predicted fluctuation value V output from the fluctuation prediction unit 2 described above and the deviation D (for example, PV-SV) used for the feedback control described above, and the specified preceding amount Qa described above. It is a functional part configured to calculate the leading index Qb for correcting the above. That is, the second calculation unit 4 outputs the calculation result of the leading index Qb with the predicted fluctuation value V and the deviation D as input values.
  • the deviation D for example, PV-SV
  • the second calculation unit 4 has a value corresponding to the predicted fluctuation value V, and is based on the value of the preceding index Qb.
  • a specified leading index calculation unit 42 (function generator) that calculates the specified leading index value Q1 and an evaluation unit that calculates an adjustment value Q2 for adjusting the above specified leading index value Q1 based on the above deviation D. 43 (see FIG. 6) and a leading index calculation unit 44 that calculates a leading index Qb based on the above-mentioned specified leading index value Q1 and adjustment value Q2 may be provided.
  • the above-mentioned defined leading index calculation unit 42 may perform each calculation using a relationship such as a function capable of calculating a unique value from an input value as shown in FIG.
  • the horizontal axis is the fluctuation value (° C.) of the main steam temperature as the predicted fluctuation value V
  • the other axis is the above-mentioned specified leading index value Q1.
  • the specified leading index value Q1 corresponding to the predicted fluctuation value V is calculated.
  • the adjustment value Q2 is the preceding amount Ff applied immediately before based on the deviation D obtained as a result of actually applying the preceding amount Ff calculated based on the predicted fluctuation value V described above to the target operation end 9. It has a value smaller than the specified leading index value Q1 (Q1> Q2). Then, by adjusting the specified leading index value Q1 with the adjusting value Q2, it has a role of adjusting the leading index Qb to a more appropriate value. In short, instead of directly calculating the leading index Qb based on the predicted fluctuation value V obtained from the learning model M, the specified leading index value Q1 calculated based on this predicted fluctuation value V was obtained immediately before that. It is adjusted in a feedback control manner according to the result of applying the preceding amount Ff to the target operation end 9.
  • the evaluation unit 43 calculates a specified adjustment value Q3 having a value corresponding to the input deviation D, which is the basis of the adjusted value Q2. It is connected to the calculation unit 43a (function generator) and the subsequent stage of this specified adjustment value calculation unit 43a, and outputs the specified adjustment value Q3 to the subsequent stage only when a specified event occurs (when the SH soot blower is operating in FIG. 6).
  • the integrated value calculation unit 43c which is connected to the output switching unit 43b and the subsequent stage of the output switching unit 43b and calculates the integrated value I by integrating the specified adjustment value Q3 input at each specified timing, and the integrated value I. It may have an adjustment value calculation unit 43d (function generator) for calculating the adjustment value Q2.
  • the above-mentioned specified adjustment value calculation unit 43a and adjustment value calculation unit 43d may perform their respective calculations using a relationship such as a function capable of calculating a unique value from the input value, respectively.
  • the output switching unit 43b integrates the output value of the specified adjustment value calculation unit 43a in the subsequent stage. It is output to the value calculation unit 43c, and when it indicates OFF, 0 (%) is output to the subsequent stage. In this way, by calculating the integrated value I, instead of modifying the function of the specified leading index calculation unit 42 (FIG. 2) to obtain the leading index Qb, an appropriate leading through the calculation of the adjusted leading index Q2.
  • the index Qb can be calculated.
  • the second calculation unit 4 further includes a holding unit 41 for holding the input predicted fluctuation value V, and is specified in advance.
  • the index calculation unit 42 calculates the specified leading index value Q1 by inputting the predicted fluctuation value V held by the holding unit 41.
  • the holding unit 41 is provided to hold the value of the predicted fluctuation value V immediately before the occurrence of the event. Since the predicted fluctuation value V also changes due to the transient change in the state quantity during the event occurrence, the holding unit 41 makes it possible to determine the preceding quantity Ff with the predicted fluctuation value V immediately before the occurrence of the event.
  • the third calculation unit 5 is a functional unit for calculating the preceding amount Ff based on the above-mentioned specified leading amount Qa and leading index Qb (see FIG. 3).
  • the leading index Qb obtained based on the predicted fluctuation value V by the second calculation unit 4 is used instead of directly calculating the leading amount Ff based on the predicted fluctuation value V obtained from the learning model M.
  • the leading amount Ff is calculated by modifying the specified leading amount Qa calculated by the first calculation unit 31 to a more appropriate one. That is, the third calculation unit 5 outputs the leading amount Ff with the specified leading amount Qa and the leading index Qb as input values.
  • the third calculation unit 5 has a specified preceding amount Qa calculated by a plurality of functions of the first calculation unit 31, and a second calculation.
  • the leading index Qb calculated by the part 4 is input.
  • one of the plurality of input specified leading amounts Qa is selected according to the value of the input leading index Qb, and the selected specified leading amount Qa is selected. It is output from the third calculation unit 5.
  • the preceding amount Ff is calculated in consideration of the load value of the boiler 71, which will be described later.
  • the preceding quantity Ff output from the preceding quantity calculation unit 3 in this way is input to the feedback control unit 6 and is added to the first manipulated variable Fb described above in the adding unit 63.
  • This addition result is transmitted (output) to the target operation end 9 as the operation amount MV.
  • the target operation end 9 is operated by the operation amount MV
  • the measured value PV of the control amount changes, and the measured value PV of the control amount and the control set value SV after this change are input to the deviation calculation unit 61.
  • the above deviation D is input to the controller 62.
  • the controller 62 calculates the first manipulated variable Fb of the target operating end 9 so that the deviation D becomes 0, and outputs the calculated result to the addition unit 63 described above.
  • the controller 62 performs PI control in FIG.
  • the addition unit 63 calculates the operation amount MV based on the first operation amount Fb and the preceding amount Ff. Feedback control is performed by repeating such control. Further, the deviation D calculated by the deviation calculation unit 61 is input to the third calculation unit 5 when the event generation signal E is on.
  • the preceding amount Ff that corrects the first manipulated variable Fb obtained by the feedback control of the controlled variable of the plant 7 is the operation data obtained in the past operation (actual operation, test run) of the plant 7. It is calculated based on the predicted fluctuation value V obtained by the learning model M that learned the relationship between the process parameter P and the fluctuation value of the controlled variable.
  • the leading amount Ff is calculated by correcting using the leading index Qb calculated based on the deviation D between the measured value PV of the controlled variable and the control set value SV.
  • a more appropriate preceding amount Ff can be obtained according to the operating state of the plant 7 and its conversion over time, and the control amount can be calculated more appropriately through the optimization of the operation amount of the target operation end 9. ..
  • the process parameter P is a measured value of one or more state quantities, a set manipulated quantity of one or more operating ends (c), l), m)) and a leading index Qb. (G)) and.
  • the fluctuation prediction unit 2 described above calculates the predicted fluctuation value V based on such a process parameter P.
  • the process parameter P may include the above b), c), d), f), g), h), and j).
  • the preceding amount calculation unit 3 calculates the preceding amount Ff with respect to the spray valve 7v, which is the target operating end 9, based on the predicted fluctuation value V of the main steam temperature. Further, this calculation result is input to the feedback control unit 6 of the valve opening degree of the spray valve 7v. Thereby, the feedforward control can be performed more appropriately so that the sudden rise or fall of the main steam temperature due to the above-mentioned disturbance is suppressed.
  • the process parameters P are b), c), d), f), g), h), i), l), m), and r) above. May include.
  • the preceding amount calculation unit 3 uses the opening degree of the gas recirculation fan inlet damper 87, which is the target operation end 9, and the burner nozzle by the burner angle adjuster (not shown) based on the predicted fluctuation value V of the reheated steam temperature.
  • the preceding amount Ff with respect to the angle is calculated respectively. Further, this calculation result is input to each feedback control unit 6 of the opening degree of the gas recirculation fan inlet damper 87 and the burner nozzle angle. Thereby, the feedforward control can be performed more appropriately so that the sudden rise or fall of the reheated steam temperature due to the above-mentioned disturbance is suppressed.
  • the process parameters P are the above a), b), c), d), f), g), h), j), k), l), m. ), And p) may be included.
  • the preceding amount calculation unit 3 calculates the preceding amount Ff for the boiler master (fuel supply device 72v), which is the target operating end 9, based on the predicted fluctuation value V of the main steam pressure. Further, this calculation result is input to the feedback control unit 6 of the valve opening degree of the fuel supply device 72v. Thereby, the feedforward control can be performed more appropriately so that the sudden increase in the main steam pressure due to the above-mentioned disturbance is suppressed.
  • the process parameters P are the above d), h), g), l), m), n), o), p), and so on. And q) may be included.
  • the preceding amount calculation unit 3 calculates the preceding amount Ff with respect to the ammonia flow rate control valve 81v, which is the target operation end 9, based on the predicted fluctuation value V of the NOx concentration. Further, this calculation result is input to the feedback control unit 6 of the valve opening degree of the ammonia flow rate control valve 81v. Thereby, the feedforward control can be performed more appropriately so that the rapid increase in the NOx concentration due to the above-mentioned disturbance is suppressed.
  • the fluctuation prediction unit 2 and the preceding amount calculation unit 3 calculate the preceding amount Ff at each specified timing. As shown in FIG. 3, when the duration t of the event is 0, the specified preceding amount Qa is calculated. Since it becomes 0, the preceding amount Ff also becomes 0.
  • the calculation of the preceding amount Ff may be performed in consideration of the load of the boiler 71.
  • the control device 1 further includes a load acquisition unit 32 that acquires the load value of the boiler 71 of the plant 7. Then, the fluctuation prediction unit 2 described above is acquired by the load acquisition unit 32 at the time of prediction based on the predicted fluctuation value V of each of the plurality of learning models M provided according to the plurality of load bands of the boiler 71. It may be configured to output the predicted fluctuation value V corresponding to the load value (hereinafter, the current load value L) of the boiler 71. Further, the above-mentioned advance amount calculation unit 3 has the above-mentioned first calculation unit 31, second calculation unit 4, and third calculation unit 5, which are provided according to the plurality of load zones, respectively, and is described above. It may be configured to calculate the preceding amount Ff according to the load value L.
  • each of the above-mentioned plurality of learning models M has the above-mentioned relationship (learning data) obtained when the boiler 71 is operated with a load value representing each load band (hereinafter, representative load value). ) Is created by learning.
  • the load of the boiler 71 is divided into four load bands, and 30%, 50%, 75%, and 100% belonging to each load band are set as representative load values, so that a total of 4 is obtained.
  • Two learning models M are created.
  • the above-mentioned process parameter P is input to each of these plurality of learning models M in parallel.
  • the fluctuation prediction unit 2 is based on the output values obtained from each of the plurality of learning models M with respect to the input value of the process parameter P and the current load value L acquired by the load acquisition unit 32.
  • the predicted fluctuation value V to be output to the second calculation unit 4 is calculated.
  • the output of the learning model M created with the representative load value having the same value as the current load value L is output as the predicted fluctuation value V.
  • the current load value L and the representative load value are not the same, the current load value L is interpolated between the output values of the two learning models M created at the representative load positions sandwiching the current load value L.
  • the calculated value according to the above is obtained and output as the predicted fluctuation value V.
  • the predicted fluctuation value V In the embodiment shown in FIG. 2, for example, with respect to the current load value L (for example, 90%), two learnings in which the representative load values of the values sandwiching the current load value are L1 (for example, 75%) and L2 (for example, 100%), respectively.
  • the predicted fluctuation value V corresponding to the current load value L is located on a straight line connecting the two coordinate points (V1, L1) and (V2, L2). become.
  • the second calculation unit 4 and the third calculation unit 5 also have the same number as the number of the learning model M according to the load band in which the load value is divided in the same manner as the learning model M.
  • the predicted fluctuation value V output from the fluctuation prediction unit 2 is input to each of the plurality of second calculation units 4 provided for each load band, and the specified leading index value Q1 is obtained from each of them. Is output.
  • These plurality of second calculation units 4 use different functions to calculate the specified leading index value Q1 according to the load band.
  • the evaluation unit 43 is provided for each of the plurality of second calculation units 4, and the specified leading index value Q1 calculated for each load band is calculated for each load band. It is adjusted according to the adjustment value Q2. At this time, each evaluation unit 43 calculates the adjustment value Q2 according to the load band, and the adjusted specified leading index value Q1 and the current load value L obtained for each load band are input to the switching unit 45. .. Then, the leading index Qb is output by performing linear interpolation or the like by the switching unit 45 in the same manner as described above.
  • This leading index Qb is input to each of a plurality of third calculation units 5 provided for each load band. Further, the first calculation unit 31 provided for each third calculation unit 5 calculates the specified preceding amount Qa according to the load band, and is connected to the subsequent stage of the third calculation unit 5. The switching unit 52 outputs the preceding amount Ff according to the load band.
  • the preceding amount Ff is calculated in consideration of the load of the plant 7 (boiler 71) at the time of its calculation. Thereby, an appropriate preceding amount Ff according to the load can be obtained.
  • the present invention is not limited to the above-described embodiment, and includes a modified form of the above-described embodiment and a combination of these embodiments as appropriate.
  • the control device (1) is It is a control device (1) that determines the operation amount (MV) of the target operation end (9) capable of controlling the control amount of the plant (7) at each specified timing.
  • the learning model (M) that learned the relationship between the process parameter (P) for specifying the operating state of the plant (7) and the fluctuation value of the controlled variable.
  • the process parameter (P) is input.
  • a fluctuation prediction unit (2) configured to output the predicted fluctuation value (V) of the controlled variable according to the control amount, and To correct the first manipulated variable (Fb) calculated according to the deviation (D) between the measured value (PV) of the controlled variable and the control set value (SV) based on the predicted variable value (V).
  • the preceding amount calculation unit (3) configured to calculate the preceding amount (Ff) of the above is provided.
  • the preceding amount calculation unit (3) When an event that changes the control amount occurs, a first calculation unit (31) configured to calculate a specified preceding amount (Qa) according to the duration (t) of the event, and a first calculation unit (31).
  • a second calculation unit (4) configured to calculate a leading index (Qb) for correcting the specified leading amount (Qa) based on the predicted fluctuation value (V) and the deviation (D). It has a third calculation unit (5) configured to calculate the preceding amount (Ff) based on the specified preceding amount (Qa) and the preceding index (Qb).
  • the preceding amount (Ff) feed forward operation amount
  • the first operation amount (Fb) obtained by the feedback control of the control amount (state amount of the controlled object) of the plant (7).
  • M learned the relationship between the process parameter (P) and the fluctuation value of the controlled variable from the operation data obtained in the past operation (actual operation, trial operation) of the plant (7). It is calculated based on the prediction result (V) of the fluctuation value of the controlled variable.
  • the operation of the decalcifying device and the occurrence of events such as ignition and extinguishing of the burner are states such as main steam temperature, main steam pressure, reheated steam temperature, and NOx concentration.
  • the predicted fluctuation value (V) obtained based on the above learning model (M) is the preceding quantity (Ff) for suppressing the influence of the disturbance on the control quantity. ) Is calculated.
  • a specified preceding amount (Qa) is calculated according to the operating time of the deashing device (event duration (t)). Then, the leading index (Qb) calculated by calculating this defined leading quantity (Qa) based on the predicted fluctuation value (V) and the deviation (D) between the measured value (PV) of the controlled quantity and the control set value (SV).
  • the preceding amount (Ff) is calculated by modifying it to a more appropriate one using. That is, instead of directly calculating the preceding amount (Ff) based on the predicted fluctuation value (V) obtained from the learning model (M), the calculated value calculated based on this predicted fluctuation value (V) is calculated.
  • the preceding amount (Ff) obtained immediately before that is adjusted in a feedback control manner according to the result of applying the preceding amount (Ff) to the target operation end (9), and the preceding index (Qb) calculated in this way is used as the above-mentioned specified preceding amount (Qa). ) Is modified to calculate the preceding amount (Ff).
  • the second calculation unit (4) The specified leading index calculation unit (42) that calculates the specified leading index value (Q1) according to the predicted fluctuation value (V), and An evaluation unit (43) that calculates an adjustment value (Q2) for adjusting the specified leading index value (Q1) based on the deviation (D). It has a leading index calculation unit (44) that sends out the leading index (Qb) based on the specified leading index value (Q1) and the adjusted value (Q2).
  • the leading index is adjusted by adjusting the uncorrected index calculated according to the predicted fluctuation value (V) with the adjustment value (Q2) calculated based on the deviation (D).
  • Qb is calculated.
  • the process parameter (P) includes a measured value of one or more state quantities, an manipulated variable of one or more operating ends installed in the plant (7), and the leading index (Qb).
  • the measured value of the state quantity at the time of prediction, the manipulated variable at the operation end, and the leading index (Qb) (previous value) are included.
  • the operating state of the plant (7) can be appropriately specified by the process parameter (P).
  • the learning model (M) can output the predicted fluctuation value (V) in consideration of the leading index (Qb) at the previous specified timing, and the learning data collection time can be increased. It is possible to shorten the learning time including the learning time and the test run time of, for example, the boiler (71) for obtaining the learning data.
  • a load acquisition unit (32) for acquiring the load value (L) of the boiler (71) of the plant (7) is further provided.
  • the fluctuation prediction unit (2) is based on the predicted fluctuation value (V) of each of the plurality of learning models (M) provided according to the plurality of load bands of the boiler (71). It is configured to output the predicted fluctuation value (V) according to (L).
  • the preceding amount calculation unit (3) includes the first calculation unit (31), the second calculation unit (4), and the third calculation unit (5), which are provided according to the plurality of load zones, respectively. Then, the preceding amount (Ff) is calculated according to the load value (L).
  • the preceding amount (Ff) is calculated in consideration of the load of the plant (7) (boiler (71)) at the time of its calculation. Thereby, an appropriate leading amount (Ff) according to the load can be obtained.
  • Control device 10 Control program 11 Processor 12 Storage unit 2 Fluctuation prediction unit 3 Preceding amount calculation unit 31 1st calculation unit 32 Load acquisition unit 4 2nd calculation unit 41 Holding unit 42 Specified leading index calculation unit 43 Evaluation unit 43a Specified adjustment value Calculation unit 43b Output switching unit 43c Integrated value calculation unit 43d Adjustment value calculation unit 44 Leading index calculation unit 45 Switching unit 5 Third calculation unit 52 Switching unit 6 Feedback control unit 61 Deviation calculation unit 62 Controller 63 Addition unit 7 Plant 7L Cooling medium Supply pipe 7c Temperature reducer 7p Steam pipe 7s Pressure sensor 7t Boiler outlet temperature measurement sensor 7v Spray valve 7w Cooling medium 71 Boiler 71f Fire furnace 71p Furnace wall evaporation pipe 72 Burning device 72v Fuel supply device 73 Smoke channel 74 Combustion burner 75 Pulverized charcoal supply Pipe 76 Crusher 77a Air duct 77b Branch 77c Branch duct 77d Damper 77e Air box 77p Port 78 Boiler 79

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)
  • Control Of Steam Boilers And Waste-Gas Boilers (AREA)

Abstract

プラントの制御量を制御可能な対象操作端の操作量を決定する制御装置であって、プラントの運転状態を特定するためのプロセスパラメータと制御量の変動値との関係性を学習した学習モデルを用いて、制御量の予測変動値を出力する変動予測部と、予測変動値に基づいて、制御量の計測値と制御設定値との偏差に応じて算出される第1操作量を修正するための先行量を算出する先行量算出部と、を備え、先行量算出部は、制御量を変化させるイベントの発生時に、イベントの継続時間に応じた規定先行量を算出する第1算出部と、予測変動値および偏差に基づいて、規定先行量を修正するための先行指標を算出する第2算出部と、規定先行量および先行指標に基づいて、先行量を算出する第3算出部と、を有する。

Description

制御装置
 本開示は、プラントにおける制御対象の状態量(制御量)の制御に関し、特に、制御量を制御可能な操作端の操作量の決定に関する。
 本願は、2020年5月1日に日本国特許庁に出願された特願2020-081362号に基づき優先権を主張し、その内容をここに援用する。
 発電プラントなどのプラントで用いられるボイラでは、主蒸気温度や主蒸気圧力、再熱蒸気温度、脱硝装置の出口におけるNOx濃度などの状態量を計測し、状態量の計測値が設定値になるようにフィードバック制御が行われる。また、このようなフィードバック制御では外乱の発生による状態量の急変に対して応答が遅れ、制御量の許容変動範囲を逸脱してしまう場合がある。このため、外乱が発生すると、その影響が制御量に表れる前に打ち消すための先行量を算出し、フィードバック制御によって演算される操作端の操作量に加えた操作量を操作端に与える、フィードフォワード制御も一緒に行う場合もある。
 上述した主蒸気などの温度や圧力、再熱蒸気温度、NOx濃度などに対する外乱は、火炉壁を形成する炉壁蒸発管や、煙道に設置された過熱器や再熱器といった各種の熱交換器の表面に付着した灰を除去するためのスートブロアの稼動や、バーナの点火、消火などのイベントによって生じる。具体的には、炉壁蒸発管のためのスートブロアの稼動によって、主蒸気温度や再熱蒸気温度などの温度が低下し、煙道に設置された熱交換器のためのスートブロアを稼働させると、これらの温度は上昇する。このため、例えば主蒸気温度を一定範囲に調整するための温度低減器におけるスプレイ量を先行的に減少あるいは増加する制御が行われる(特許文献1参照)。
特許第6087226号公報
 上述したように、操作端の操作量をフィードバック制御およびフィードフォワード制御により決定する場合であっても、最適な先行量は運転状態やその経時変化により変わる。よって、運転状態等に応じて先行量を算出する必要がある。
 上述の事情に鑑みて、本発明の少なくとも一実施形態は、プラントの制御量のフィードフォワード制御のための先行量をより適切に決定することが可能な制御装置を提供することを目的とする。
 本発明の少なくとも一実施形態に係る制御装置は、
 プラントの制御量を制御可能な対象操作端の操作量を規定タイミング毎に決定する制御装置であって、
 前記プラントの運転状態を特定するためのプロセスパラメータと前記制御量の変動値との関係性を学習した学習モデルを用いて、前記プロセスパラメータの入力に応じた前記制御量の予測変動値を出力するよう構成された変動予測部と、
 前記予測変動値に基づいて、前記制御量の計測値と制御設定値との偏差に応じて算出される第1操作量を修正するための先行量を算出するよう構成された先行量算出部と、を備え、
 前記先行量算出部は、
 前記制御量を変化させるイベントの発生時に、前記イベントの継続時間に応じた規定先行量を算出するよう構成された第1算出部と、
 前記予測変動値および前記偏差に基づいて、前記規定先行量を修正するための先行指標を算出するよう構成された第2算出部と、
 前記規定先行量および前記先行指標に基づいて、前記先行量を算出するよう構成された第3算出部と、を有する。
 本発明の少なくとも一実施形態によれば、プラントの制御量のフィードフォワード制御のための先行量をより適切に決定することが可能な制御装置が提供される。
本発明の一実施形態に係るプラント7が備えるボイラの構成を概略的に示す図である。 本発明の一実施形態に係る制御装置の構成を示すブロック図である。 本発明の一実施形態に係る第1算出部および第3算出部の構成を示す図であり、図2のAで示す部分の詳細である。 本発明の一実施形態に係る第1算出部によって算出される規定先行量とイベントの継続時間との関係を例示する図である。 本発明の一実施形態に係る規定先行指標算出部によって算出される規定先行指標値と予測変動値との関係を例示する図である。 本発明の一実施形態に係る評価部の構成を示す図である。
 以下、添付図面を参照して本発明の幾つかの実施形態について説明する。ただし、実施形態として記載されている又は図面に示されている構成部品の寸法、材質、形状、その相対的配置等は、本発明の範囲をこれに限定する趣旨ではなく、単なる説明例にすぎない。
 例えば、「ある方向に」、「ある方向に沿って」、「平行」、「直交」、「中心」、「同心」或いは「同軸」等の相対的或いは絶対的な配置を表す表現は、厳密にそのような配置を表すのみならず、公差、若しくは、同じ機能が得られる程度の角度や距離をもって相対的に変位している状態も表すものとする。
 例えば、「同一」、「等しい」及び「均質」等の物事が等しい状態であることを表す表現は、厳密に等しい状態を表すのみならず、公差、若しくは、同じ機能が得られる程度の差が存在している状態も表すものとする。
 例えば、四角形状や円筒形状等の形状を表す表現は、幾何学的に厳密な意味での四角形状や円筒形状等の形状を表すのみならず、同じ効果が得られる範囲で、凹凸部や面取り部等を含む形状も表すものとする。
 一方、一の構成要素を「備える」、「具える」、「具備する」、「含む」、又は、「有する」という表現は、他の構成要素の存在を除外する排他的な表現ではない。
 図1は、本発明の一実施形態に係るプラント7が備えるボイラ71の構成を概略的に示す図である。
 制御装置1は、自装置が制御対象とするプラント7の状態量(以下、制御量)を制御可能な操作端(以下、対象操作端9)の操作量MV(指令値)を、例えば周期的などの規定タイミング毎に決定可能な装置である。一般にプラント7には、例えば外部から与えられた操作量MVに従って例えば弁体(弁開度)などの可動部を変化させることが可能な複数の操作端や、温度や圧力、流量などの各種の状態量を計測するための複数のセンサが設置されている。制御装置1は、プラント7に設置された複数の操作端のうちの1以上となる対象操作端9の操作量MVを、上記のセンサを用いて得られる制御量の計測値PVが制御設定値SV(目標値)になるように、フィードバック制御およびフィードフォワード制御により決定する。また、制御装置1は、決定した操作量MVを外部に出力(送信)することで、対象操作端9の制御を行うことが可能となる。
 図1に示すプラント7は、例えばボイラ・タービン発電設備(BTG)などのボイラ71を備えた発電プラントである。発電プラントは、ボイラ71で生成した過熱蒸気(主蒸気)を例えば蒸気タービンTに供給してその回転軸を回転駆動させることで、蒸気タービンTの回転軸に連結した発電機Gを回転駆動して発電を行う。このようなボイラ71を備えるプラント7において、制御装置1は、主蒸気温度や主蒸気圧力、再熱蒸気温度、脱硝装置の出口におけるNOx濃度などの状態量の少なくとも1つを制御量としても良い。これらの制御量を制御するための対象操作端9は、それぞれ、後述するようなスプレイ弁7vや燃料供給装置72v、ガス再循環ファン86の入口に設置されたガス再循環ファン入口ダンパ87あるいはバーナ角度調節器(不図示)、アンモニア流量制御弁81vなどとなる。
 図1に示す実施形態では、ボイラ71は、例えば石炭を粉砕した微粉炭燃料を火炉71fの内部で燃焼させることにより発生させた熱を給水や蒸気と熱交換して、蒸気タービンTなどに供給する過熱蒸気(主蒸気)を生成するよう構成された微粉炭焚きボイラである。ただし、ボイラ71は、使用燃料や燃焼方式の異なるソーダ回収ボイラなどの廃熱回収ボイラ、流動床、流動層ボイラなどであっても良い。
 まず、ボイラ71について、図1に示す微粉炭炊きボイラを例に詳細に説明する。図1に示すように、ボイラ71(ボイラ本体)は、火炉71fと、火炉71fに設置された燃焼装置72と、ボイラ71に接続された煙道73と、を有する。
 火炉71fは、内部に燃焼室を有し、例えば四角筒の中空形状をなして鉛直方向に沿って設置される(図1参照)。燃焼室を形成する火炉71fの壁(火炉壁)の内壁面は、水から蒸気を生成するための炉壁蒸発管71p(伝熱管)と、炉壁蒸発管71pを接続するフィンとで構成される。具体的には、例えば鉛直方向に沿って配置される炉壁蒸発管71pが水平方向に並んで複数配置されると共に、これら複数の炉壁蒸発管71pにおける隣接する炉壁蒸発管71pと炉壁蒸発管71pとの間をフィンが閉塞するように配置される。火炉71fは、炉底に傾斜面が設けられており、傾斜面にも炉壁蒸発管71p(炉底蒸発管)が設けられて底面となる。また、燃焼室の内部には、炉壁蒸発管71pの表面に付着する灰や未燃分など付着物を除去するよう構成された不図示のスートブロアなどの除灰装置が設けられている。
 燃焼装置72は、上記の火炉71fの内部に燃料を供給して燃焼させるための装置であり、火炉壁に設置される1以上の燃焼バーナ74と、風箱77eと、ボイラ71に供給する燃料量を調整するための燃料供給装置72vと、を備える。具体的には、火炉71fの内部に燃焼バーナ74を介して燃料を供給すると共に、風箱77eを介して燃焼用空気(2次空気)を供給することにより、燃料を燃焼させるように構成される。この燃焼バーナ74に供給される燃料量は、主蒸気圧力の計測値PV(制御量)が制御設定値SVになるように、燃料供給装置72v(操作端)が制御されることで調整される。図1に示す実施形態では、炉壁蒸発管71pで生成された蒸気は蒸気管7pを流れる際に過熱器79aで過熱された後に蒸気タービンTに送られるが、主蒸気圧力は、蒸気管7pにおけるボイラ出口に設置された圧力センサ7sによって検出されるようになっている。
 この燃焼装置72は、図1に示すように、火炉壁に設置されるアディショナルエア(AA)を供給するためのAAポート77pを、さらに備えていても良い。この場合、燃焼バーナ74に燃料を供給するための後述する搬送用空気(1次空気)、及び風箱77eを介して火炉71fに供給される燃焼用空気(2次空気)によって、燃焼用空気の供給量が微粉炭燃料の供給量に対して理論空気量未満となるように設定されることで、内部が還元雰囲気に保持されて、微粉炭燃料を還元燃焼した後、新たに燃焼用空気(AA)を供給して微粉炭燃料の酸化燃焼が完結され、燃焼バーナ74側でのNOx発生量の低減が可能となっている(2段燃焼方式)。
 本実施形態では、燃焼装置72は複数の燃焼バーナ74を備えている。そして、火炉壁には、火炉71fの周方向または、対向位置に例えば均等間隔などで1以上(図1では複数)の燃焼バーナ74が設置されると共に、このような周方向または対向位置に設置された1以上の燃焼バーナ74が鉛直方向に沿って1以上の段数(図1では74a~74eの5段)で設置されている。これら複数の燃焼バーナ74は、それぞれ、微粉炭供給管75を介して複数の粉砕機76(図1では76a~76e。微粉炭機/ミル)のいずれかに連結されており、各粉砕機76で生成された微粉炭が搬送用空気(1次空気)によって、上記の燃料供給装置72vが設置された微粉炭供給管75を介して燃焼バーナ74に供給されるようになっている。なお、上記の粉砕機76は、図示しない原炭搬送系統により供給される石炭を所定の微粉の大きさに粉砕することで、微粉炭を生成するようになっている。
 また、上記の風箱77eは、火炉71fにおける各燃焼バーナ74の設置位置に設けられており、空気(燃焼用空気)を導く空気ダクト77aの一端部が連結されている。この空気ダクト77aの他端部は送風機78に接続されており、送風機78によって押し込まれた空気が空気ダクト77aを通って風箱77eに供給されるようになっている。また、空気ダクト77aに設けられた分岐部77bから分岐ダクト77cが分岐して上記のAAポート77pに接続されることで、分岐ダクト77cに設置されたAAダンパ77dの開度に応じたAAをAAポート77p(火炉71fの内部)に供給することが可能に構成されている。つまり、送風機78から空気ダクト77aに供給された空気は、エアヒータ82(後述)で温められた後、分岐部77bにおいて風箱77eへ導かれる2次空気と、分岐ダクト77cを経由してAAポート77pへと導かれるAAとに分配されるようになっている。
 煙道73は、火炉71fにおける燃料の燃焼により生じた燃焼ガスを導くダクトであり、図1に示すように火炉71fの鉛直方向上方に連結される。この煙道73には、蒸気の生成や、蒸気の加熱(過熱)するための複数の熱交換器79(図1では、過熱器79a(SH)、再熱器79b(RH)、節炭器79c(ECO))や、熱交換器79の伝熱面に付着する付着物(灰など)を除去するために各熱交換器79に対してそれぞれ設けられる除灰装置(不図示)が設置される。そして、火炉71fから煙道73に流入した燃焼ガスにより火炉壁及び熱交換器79を流れる給水や蒸気を加熱して、主蒸気などを生成するようになっている。
 より詳細には、例えば貫流型ボイラの場合には、炉壁蒸発管71pの内部の水が沸騰して生成された水(飽和水)と蒸気(飽和蒸気)が混合した気液混合流体が湿分(汽水)分離器(不図示)に流入し、気液混合流体が蒸気(気相)と水(液相)に分離される。その後、湿分(汽水)分離器で水と分離された蒸気(主蒸気)は、蒸気の流路を形成する蒸気管7p(過熱蒸気管)を通って蒸気タービンTなどに送られる際に、蒸気管7pに設置された過熱器による加熱と、温度低減器7c(過熱低減器)による必要に応じた冷却とを受けて、主蒸気温度(制御量)は制御設定値SVで一定になるように調整がなされる。
 図1に示す実施形態では、上記の温度低減器7cには、冷却媒体7wを供給する冷却媒体供給管7Lが接続されており、スプレイ弁7v(操作端)によって、温度低減器7cに供給する冷却媒体7wの流量が調整可能になっている。そして、蒸気管7pに設置された温度計(ボイラ出口温度計測センサ7t)の計測値PV(制御量)が制御設定値SVになるように、スプレイ弁7v(操作端)の弁開度が調整されるようになっている。
 一方、再熱器79bにより加熱された蒸気温度(再熱蒸気温度)については、再熱蒸気温度(制御量)の計測値PVが制御設定値SVになるように、ガス再循環ファン入口ダンパ87(操作端)や、バーナノズル角度を調整可能な不図示のバーナ角度調節器(操作端)が制御される。
 また、煙道73の下流側には排ガス処理装置が接続されることにより、上記の熱交換器79を通過した後の燃焼ガス(排ガス)は、排ガス処理装置により無害化された後、大気に放出される。より詳細には、排ガス処理装置は、煙道73に連結される排ガス通路8と、排ガス通路8に設置される複数の機器で構成される。図1に示す実施形態では、排ガス通路8には、排ガスの上流側から下流側に向けて、順番に、脱硝装置81、送風機78から空気ダクト77aへ供給される空気と排ガス通路8を流れる排ガスとの間で熱交換を行うエアヒータ82、煤塵処理装置83、誘引送風機84などが設けられ、下流端部に煙突85が設けられている。
 上記の脱硝装置81は、内部に設置された触媒にアンモニア8w等を供給することで、触媒を通過する排ガスに含まれるNOxを浄化する。具体的には、脱硝装置81には、アンモニア8wを供給するための配管(アンモニア供給配管81p)が接続されており、アンモニア流量制御弁81v(操作端)の弁開度を調整することで、脱硝装置81へのアンモニア8wの供給量を制御可能になっている。そして、脱硝装置81の出口には、NOx濃度を検出可能なセンサ(NOxセンサ81s)が設置されており、NOxセンサ81sによる計測値PVが制御設定値SVになるように、アンモニア流量制御弁81v(操作端)が制御される。
 以上の通り、プラント7は、各々が複数の機器で構成されるボイラ71や排ガス処理装置、蒸気タービンT、操作端、センサなどの複数の機器を備えている。そして、プラント7(ボイラ71)の例えば電力需要などに基づく出力要求を満たしつつ、プラント7の運転が最適となるように、上記の制御装置1などによる制御量の制御がなされる。具体的には、上記の制御装置1などによって、フィードバック制御に加えて、フィードフォワード制御による制御が行われる。
 ここで、フィードフォワード制御は、制御量のフィードバック制御における外乱を生じされるイベントの発生による影響を抑制するために行われる。例えば主蒸気温度の計測値PVが制御設定値SVとなるように、スプレイ弁7vに対して、フィードバック制御により求まる操作量(以下、第1操作量Fb)が規定タイミングで決定される。この際、炉壁蒸発管71pの除灰装置が稼働すると(イベントの発生)、その付着物が除去されることで、火炉71fでの燃焼により生じた熱と蒸発管71p内の流体(水)との間の熱交換量が急激に増大する。これがフィードバック制御に対する外乱となって、煙道73に向けて流れる燃焼ガスの温度がその分だけ大きく低下し、主蒸気温度の急激な低下及び蒸気圧力の急激な上昇を招来する。また、過熱器79aの除灰装置が稼働すると(イベントの発生)、その熱交換量が急激に増大するため、これがフィードバック制御に対する外乱となって、主蒸気に対する加熱が過剰になり、主蒸気温度の急激な上昇を招来する。
 同様に、再熱蒸気温度の計測値PVが制御設定値SVとなるように、ガス再循環ファン入口ダンパ87あるいはバーナ角度調節器(不図示)の少なくとも一方に対する第1操作量Fbが規定タイミングで決定される。この際、炉壁蒸発管71pの除灰装置の稼働(イベントの発生)が上記の外乱となって、再熱蒸気温度は急激に低下する。また、再熱器79bの除灰装置の稼働が上記の外乱となって、再熱蒸気温度は急激に上昇する。
 主蒸気圧力については、主蒸気圧力の計測値PVが制御設定値SVになるように、燃料供給装置72vに対する第1操作量Fbが規定タイミングで決定されるが、炉壁蒸発管71pの稼動(イベントの発生)が上記の外乱となって、主蒸気圧は急激に上昇する。
 また、脱硝装置81の出口におけるNOx濃度については、そのNOx濃度の計測値PVが制御設定値SVになるように、アンモニア流量制御弁81vに対する第1操作量Fbが規定タイミングで決定されるが、燃焼バーナ74の点火または消火(イベントの発生)が上記の外乱となって、蒸気のNOx濃度が急激に上昇する。
 そして、上述したような制御量のフィードバック制御における外乱に対して、フィードフォワード制御によって外乱の発生に応じて先行的に与える操作量(以下、先行量Ff)を演算するが、その最適値は、プラント7(ボイラ71)の運転状態やその経時変化により変わる。よって、制御装置1は、プラント7(ボイラ71)の運転状態等に応じた適切な先行量Ffを算出するために、次に説明するように先行量Ffを算出(演算)する。
 以下、上記の制御装置1について、図2~図6を用いて詳細に説明する。
 図2は、本発明の一実施形態に係る制御装置1の構成を示すブロック図である。図3は、本発明の一実施形態に係る第1算出部31および第3算出部5の構成を示す図であり、図2のAで示す部分の詳細である。図4は、本発明の一実施形態に係る第1算出部31によって算出される規定先行量Qaとイベントの継続時間tとの関係を例示す図である。図5は、本発明の一実施形態に係る規定先行指標算出部42によって算出される規定先行指標値Q1と予測変動値Vとの関係を例示する図である。また、図6は、本発明の一実施形態に係る評価部43の構成を示す図である。
 図2に示すように、制御装置1は、変動予測部2と、先行量算出部3と、を備える。これらの機能部(回路)についてそれぞれ説明する。
 なお、制御装置1は、コンピュータで構成されても良い。CPUなどのプロセッサ11や、ROMやRAMといったメモリなどの記憶部12を備えている。そして、メモリ(主記憶部)にロードされたプログラム(制御プログラム10)の命令に従ってプロセッサ11が動作(データの演算など)することで、上記の各機能部を実現する。例えば、MCU(Micro Controller Unit)などで実現しても良いし、SoC(System on Chip)など1つの集積回路(チップ)に実装しても良い。また、上記のプログラムは、コンピュータに後述する各機能部を実現させるためのソフトウェアであり、コンピュータによる読み込みが可能な記憶媒体に記憶されても良い。
 変動予測部2は、プラント7の運転状態を特定するためのプロセスパラメータPと制御量の変動値との関係性を学習した学習モデルMを用いて、プロセスパラメータPの入力に応じた制御量の予測変動値Vを出力するよう構成された機能部である。プロセスパラメータPは、状態量の計測値と、操作端の操作量との少なくとも一方を含む。状態量の計測値は図示あるいは不図示のセンサによって計測可能であり、操作端の操作量は、操作端あるいは操作端の制御部から取得可能である。
 例えば、プロセスパラメータPは、a)主蒸気圧力、b)節炭器79cの出口における排ガス温度、c)スプレイ弁7vの弁開度、d)燃料石炭流量(トン/時間)、e)給水流量(トン/時間)、f)主蒸気温度、g)前回の先行指標Qb(後述)、h)発電機出力(負荷指標)、i)再熱蒸気温度、j)過熱器79a(1次過熱器など)の入口蒸気温度、K)過熱器79a(1次過熱器など)の出口蒸気温度、l)ガス再循環ファン入口ダンパ87の開度、m)燃焼バーナ74のバーナノズル角度(バーナ角度調節器による設定角度)、n)燃料ガス流量、o)燃料油流量、p)主蒸気流量、q)全空気流量(上記の分岐部77bに流入する総流量)、r)再熱器79bの出口温度のうちの少なくとも1つのパラメータ値を含む。なお、上記のc)、g)、l)、m)が操作端の操作量であり、それ以外が状態量である。
 また、上記の学習モデルMは、上述した1以上のパラメータのパラメータ値を含むプロセスパラメータPと、このプロセスパラメータPの下で計測された制御量の変動値との関係性を学習(機械学習)することにより作成される。この学習は、学習データ(教師データ)に対して周知な機械学習の手法により行われれば良い。学習データは、例えば過去にプラント7(ボイラ71)を運転(実運転、試運転)した際に得られたプロセスパラメータPの各パラメータ値と、これが得られた際の制御量の変動値の計測結果とを対応(関連)付けた複数のデータで構成されても良い。あるいは、学習データは、文献あるいは数値解析(材料シミュレーション)などから得られるプロセスパラメータPの各パラメータ値と制御量の変動値とを対応付けた複数のデータで構成されても良い。これらのデータを組み合せても良い。こうして作成された学習モデルMは、入力されたプロセスパラメータPに対して、上記の関係性から得られる制御量の変動値を予測変動値Vとして出力する。
 先行量算出部3は、上述した予測変動値Vに基づいて、制御量の計測値PVと制御設定値SVとの偏差Dに応じて算出される第1操作量Fbを修正するための先行量Ffを算出するよう構成された機能部である。第1操作量Fbは、上記の偏差Dに基づいて、フィードバック制御により算出されるフィードバック操作量である。先行量Ffによって修正された値が最終的な操作量MVとして、この操作量MVによって操作される対象操作端9に出力される(図2参照)。
 詳述すると、上記の先行量算出部3は、図2に示すように、上記の先行量Ffを算出するために、第1算出部31(図3参照)と、第2算出部4と、第3算出部5(図3参照)との3つの機能部(回路)を有する。これらの機能部について、イベントが除灰装置の稼働であり、制御量が主蒸気温度、対象操作端9がスプレイ弁7vである場合を例に説明する。
 第1算出部31は、制御量を変化させる例えば上述した除灰装置の稼働、燃焼バーナ74の点火、消火などのイベントの発生時に、イベントの継続時間tに応じた規定先行量Qaを算出するよう構成された機能部(関数発生器)である。より詳細には、図3に示すように、第1算出部31には、イベントの継続時間tとして過熱器79aの除灰装置(本実施形態ではスートブロア)の作動時間が入力される。そして、入力されたイベントの継続時間tから一意に規定先行量Qaを算出可能な例えば制御理論や経験などに基づいて作成された関数(図4参照)やテーブルなどの関係を用いて、規定タイミング毎に進むイベントの継続時間tに応じた規定先行量Qaを算出する。
 図2~図6に示す実施形態では、図3に示すように、第1算出部31は、互いに異なる複数(図3では3つ)の関数を用いて、入力されたイベントの継続時間tに対する規定先行量Qaをそれぞれ算出するようになっている。
 第2算出部4は、上述した変動予測部2から出力される予測変動値Vと、上述したフィードバック制御に用いられる上記の偏差D(例えばPV-SV)に基づいて、上記の規定先行量Qaを修正するための先行指標Qbを算出するよう構成された機能部である。つまり、第2算出部4は、予測変動値Vおよび偏差Dを入力値として、先行指標Qbの演算結果を出力する。
 より詳細には、幾つかの実施形態では、図2に示すように、第2算出部4は、上記の予測変動値Vに応じた値を有し、上記の先行指標Qbの値の基になる規定先行指標値Q1を算出する規定先行指標算出部42(関数発生器)と、上記の偏差Dに基づいて、上記の規定先行指標値Q1を調整するための調整値Q2を算出する評価部43(図6参照)と、上記の規定先行指標値Q1および調整値Q2に基づいて先行指標Qbを算出する先行指標算出部44と、を備えても良い。
 上記の規定先行指標算出部42は、図5に示すような、入力値から一意な値を算出可能な関数などの関係を用いて各々の算出を行っても良い。図5では、横軸が予測変動値Vとしての主蒸気温度の変動値(℃)であり、他軸が上記の規定先行指標値Q1であり、このような関係を示す関数を用いて、入力された予測変動値Vに応じた規定先行指標値Q1を算出する。
 また、上記の調整値Q2は、上述した予測変動値Vに基づいて算出した先行量Ffを実際に対象操作端9に適用した結果として得られる偏差Dに基づいて、直前に適用した先行量Ffを評価したものであり、規定先行指標値Q1よりも小さい値を有する(Q1>Q2)。そして、規定先行指標値Q1を調整値Q2によって調整することで、先行指標Qbをより適切な値に調整する役割を有する。要するに、学習モデルMから求めた予測変動値Vに基づいて直接的に先行指標Qbを算出するのではなく、この予測変動値Vに基づいて算出した規定先行指標値Q1を、その直前に求めた先行量Ffを対象操作端9に適用した結果によりフィードバック制御的に調整する。
 具体的には、上記の評価部43は、図6に示すように、上記の調整値Q2の基になる、入力された偏差Dに応じた値を有する規定調整値Q3を算出する規定調整値算出部43a(関数発生器)と、この規定調整値算出部43aの後段に接続され、規定されたイベントの発生時(図6ではSHスートブロアの稼動時)のみ規定調整値Q3を後段に出力する出力切替部43bと、この出力切替部43bの後段に接続され、規定タイミング毎に入力される規定調整値Q3を積算した積算値Iを算出する積算値算出部43cと、この積算値Iに応じた調整値Q2を算出する調整値算出部43d(関数発生器)と、を有しても良い。上記の規定調整値算出部43aおよび調整値算出部43dは、それぞれ、入力値から一意な値を算出可能な関数などの関係を用いて各々の算出を行っても良い。
 図2~図6に示す実施形態では、図6に示すように、出力切替部43bは、イベントの発生信号EがONを示す場合には、規定調整値算出部43aの出力値を後段の積算値算出部43cに出力し、OFFを示す場合には0(%)を後段に出力するようになっている。このように、積算値Iを算出することで、規定先行指標算出部42(図2)の関数などを修正して先行指標Qbを求めるのではなく、調整値Q2の算出を通した適切な先行指標Qbの算出が可能となる。
 また、図2~図6に示す実施形態では、図2に示すように、第2算出部4は、入力される予測変動値Vを保持するための保持部41をさらに備えており、規定先行指標算出部42は、保持部41に保持された予測変動値Vを入力として、規定先行指標値Q1を算出するようになっている。この保持部41は、イベントの発生直前時の予測変動値Vの値を保持するために設けられている。イベント発生中の過渡的な状態量の変化で予測変動値Vも変化するため、保持部41によって、イベントの発生直前の予測変動値Vで先行量Ffを決定することが可能となる。
 第3算出部5は、上述した規定先行量Qaおよび先行指標Qbに基づいて、先行量Ffを算出するための機能部である(図3参照)。換言すれば、学習モデルMから求めた予測変動値Vに基づいて直接的に先行量Ffを算出するのではなく、第2算出部4により予測変動値Vに基づいて求められた先行指標Qbによって、第1算出部31により演算された規定先行量Qaをより適切なものに修正することで、先行量Ffを算出する。つまり、第3算出部5は、規定先行量Qaおよび先行指標Qbを入力値として、先行量Ffを出力する。
 図2~図6に示す実施形態では、図3に示すように、第3算出部5には、第1算出部31が有する複数の関数によってそれぞれ算出された規定先行量Qaと、第2算出部4によって算出された先行指標Qbとが入力されるようになっている。そして、第3算出部5においては、入力された先行指標Qbの値によって、入力された複数の規定先行量Qaのいずれかが選択されるようになっており、選択された規定先行量Qaが第3算出部5から出力される。なお、本実施形態ではボイラ71の負荷値を考慮して先行量Ffを算出するが、これについては後述する。
 このようにして先行量算出部3から出力された先行量Ffは、フィードバック制御部6に入力されることで、加算部63において上述した第1操作量Fbに加算されるようになっており、この加算結果が操作量MVとして対象操作端9に送信(出力)される。対象操作端9が操作量MVで操作されることで制御量の計測値PVが変化するが、この変化後の制御量の計測値PVおよび制御設定値SVが偏差算出部61に入力されて、上記の偏差Dがコントローラ62に入力される。コントローラ62では、偏差Dが0となるように、対象操作端9の第1操作量Fbを演算し、この演算結果を上述した加算部63に出力する。なお、図2では、コントローラ62はPI制御を行うが、PID制御であっても良い。そして、加算部63において第1操作量Fbと先行量Ffとに基づいて操作量MVが演算される。このような制御が繰り返されることで、フィードバック制御が行われる。また、上記の偏差算出部61によって算出された偏差Dは、イベントの発生信号Eがオンの時に第3算出部5に入力されるようになっている。
 上記の構成によれば、プラント7の制御量のフィードバック制御により求まる第1操作量Fbを修正する先行量Ffを、過去のプラント7の運転(実運転、試運転)等で得られた運転データなどからプロセスパラメータPと制御量の変動値との関係性を学習した学習モデルMにより得られる予測変動値Vに基づいて算出する。具体的には、制御理論や経験などに基づいて定めた、除灰装置の稼働、燃焼バーナ74の点火あるは消火といったイベントの継続時間に応じた規定先行量Qaを、予測変動値V、および制御量の計測値PVと制御設定値SVとの偏差Dに基づいて算出した先行指標Qbを用いて修正することで、先行量Ffを算出する。これによって、プラント7の運転状態やその経時変換に応じたより適切な先行量Ffを求めることができ、対象操作端9の操作量の適正化を通して、制御量の演算をより適切に行うことができる。
 次に、上述したプロセスパラメータPと制御量との関係について詳細に説明する。
 幾つかの実施形態では、プロセスパラメータPは、1以上の状態量の計測値と、1以上の操作端の設定済みの操作量(上記のc)、l)、m))と、先行指標Qb(上記g))と、を含んでいる。そして、上述した変動予測部2はこのようなプロセスパラメータPに基づいて予測変動値Vを算出する。
 具体的には、制御量が主蒸気温度である場合には、プロセスパラメータPは、上記のb)、c)、d)、f)、g)、h)、およびj)を含んでも良い。そして、先行量算出部3は、主蒸気温度の予測変動値Vに基づいて、その対象操作端9となるスプレイ弁7vに対する先行量Ffを算出する。また、この算出結果は、スプレイ弁7vの弁開度のフィードバック制御部6に入力される。これによって、上記の外乱による主蒸気温度の急激な上昇または低下が抑制されるようにフィードフォワード制御をより適切に行うことができる。
 制御量が再熱蒸気温度である場合には、プロセスパラメータPは、上記のb)、c)、d)、f)、g)、h)、i)、l)、m)、およびr)を含んでも良い。そして、先行量算出部3は、再熱蒸気温度の予測変動値Vに基づいて、その対象操作端9となるガス再循環ファン入口ダンパ87の開度およびバーナ角度調節器(不図示)によるバーナノズル角度に対する先行量Ffをそれぞれ算出する。また、この算出結果は、上記のガス再循環ファン入口ダンパ87の開度およびバーナノズル角度の各々のフィードバック制御部6に入力される。これによって、上記の外乱による再熱蒸気温度の急激な上昇または低下が抑制されるようにフィードフォワード制御をより適切に行うことができる。
 制御量が主蒸気圧力である場合には、プロセスパラメータPは、上記のa)、b)、c)、d)、f)、g)、h)、j)、k)、l)、m)、およびp)を含んでも良い。そして、先行量算出部3は、主蒸気圧力の予測変動値Vに基づいて、その対象操作端9となるボイラマスタ(燃料供給装置72v)に対する先行量Ffをそれぞれ算出する。また、この算出結果は、燃料供給装置72vの弁開度のフィードバック制御部6に入力される。これによって、上記の外乱による主蒸気圧力の急激な上昇が抑制されるようにフィードフォワード制御をより適切に行うことができる。
 また、制御量が脱硝装置81の出口におけるNOx濃度である場合には、プロセスパラメータPは、上記のd)、h)、g)、l)、m)、n)、o)、p)、およびq)を含んでも良い。そして、先行量算出部3は、上記のNOx濃度の予測変動値Vに基づいて、その対象操作端9となるアンモニア流量制御弁81vに対する先行量Ffを算出する。また、この算出結果は、アンモニア流量制御弁81vの弁開度のフィードバック制御部6に入力される。これによって、上記の外乱によるNOx濃度の急激な上昇が抑制されるようにフィードフォワード制御をより適切に行うことができる。
 なお、変動予測部2および先行量算出部3は、規定タイミング毎に先行量Ffの算出を行うが、図3に示すように、イベントの継続時間tが0の場合には規定先行量Qaは0となるので、先行量Ffも0となる。
 次に、上記の先行量Ffの算出をボイラ71の負荷を考慮して行う場合の実施形態について説明する。
 幾つかの実施形態では、図2に示すように、上記の先行量Ffの算出はボイラ71の負荷を考慮して行っても良い。
 具体的には、幾つかの実施形態では、図2に示すように、制御装置1は、プラント7が有するボイラ71の負荷値を取得する負荷取得部32をさらに備える。そして、上述した変動予測部2は、ボイラ71の複数の負荷帯に応じて設けられた複数の学習モデルMの各々の予測変動値Vに基づいて、予測時において上記の負荷取得部32によって取得されたボイラ71の負荷値(以下、現負荷値L)に応じた予測変動値Vを出力するように構成されても良い。さらに、上述した先行量算出部3は、上記の複数の負荷帯に応じてそれぞれ設けられた上述した第1算出部31、第2算出部4および第3算出部5を有し、上記の現負荷値Lに応じた先行量Ffを算出するよう構成されても良い。
 より詳細には、上記の複数の学習モデルMは、それぞれ、各負荷帯を代表する負荷値(以下、代表負荷値)でボイラ71が運転されている際に得られる上述の関係性(学習データ)を学習することで作成される。図2に示す実施形態では、ボイラ71の負荷を4つの負荷帯に分けると共に、各負荷帯に属する30%、50%、75%、100%をそれぞれ代表負荷値とすることで、合計で4つの学習モデルMを作成している。
 また、上述したプロセスパラメータPは、これらの複数の学習モデルMの各々にそれぞれ並列に入力される。そして、変動予測部2は、このプロセスパラメータPの入力値に対して複数の学習モデルMの各々からそれぞれ得られる出力値と、負荷取得部32によって取得された現負荷値Lとに基づいて、第2算出部4に対して出力する予測変動値Vを算出する。
 具体的には、現負荷値Lと代表負荷値とが同じ場合には、現負荷値Lと同じ値を有する代表負荷値で作成された学習モデルMの出力が、予測変動値Vとして出力される。現負荷値Lと代表負荷値とが同じではない場合には、現負荷値Lを挟む代表負荷位置で作成された2つの学習モデルMの出力値の間を補間することで、現負荷値Lに応じた算出値を求めて、予測変動値Vとして出力する。図2に示す実施形態において、例えば現負荷値L(例えば90%)に対して、これを挟む値の代表負荷値をそれぞれL1(例えば75%)、L2(例えば100%)とする2つの学習モデルMの出力値がそれぞれV1、V2とすると、現負荷値Lに対応する予測変動値Vは、(V1、L1)、(V2、L2)の2つの座標点を結ぶ直線上に位置することになる。
 同様に、第2算出部4および第3算出部5も、図2に示すように、学習モデルMと同じように負荷値を分割した負荷帯に応じて、学習モデルMの数と同じ数だけそれぞれ設けられる。具体的には、変動予測部2から出力された予測変動値Vは、各負荷帯に対して設けられた複数の第2算出部4にそれぞれ入力されると共に、その各々から規定先行指標値Q1が出力される。これらの複数の第2算出部4は、互いに異なる関数を用いて、負荷帯に応じた規定先行指標値Q1を算出する。
 また、評価部43は、図2に示すように、複数の第2算出部4に対してそれぞれ設けられており、負荷帯毎に算出された規定先行指標値Q1が、負荷帯毎に算出された調整値Q2によってそれぞれ調整される。この際、各評価部43は、負荷帯に応じた調整値Q2を算出すると共に、負荷帯毎にそれぞれ求めた調整後の規定先行指標値Q1および現負荷値Lが切替部45に入力される。そして、切替部45によって上記と同様に線形補間などされることで、先行指標Qbが出力される。
 この先行指標Qbは、負荷帯毎に設けられた複数の第3算出部5にそれぞれ入力される。また、各第3算出部5に対して設けられた第1算出部31は、負荷帯の応じた規定先行量Qaを算出するようになっており、第3算出部5の後段に接続された切替部52によって、負荷帯に応じた先行量Ffが出力される。
 上記の構成によれば、先行量Ffは、その算出時におけるプラント7(ボイラ71)の負荷を考慮して算出される。これによって、負荷に応じた適切な先行量Ffを求めることができる。
 本発明は上述した実施形態に限定されることはなく、上述した実施形態に変形を加えた形態や、これらの形態を適宜組み合わせた形態も含む。
<付記>
(1)本発明の少なくとも一実施形態に係る制御装置(1)は、
 プラント(7)の制御量を制御可能な対象操作端(9)の操作量(MV)を規定タイミング毎に決定する制御装置(1)であって、
 前記プラント(7)の運転状態を特定するためのプロセスパラメータ(P)と前記制御量の変動値との関係性を学習した学習モデル(M)を用いて、前記プロセスパラメータ(P)の入力に応じた前記制御量の予測変動値(V)を出力するよう構成された変動予測部(2)と、
 前記予測変動値(V)に基づいて、前記制御量の計測値(PV)と制御設定値(SV)との偏差(D)に応じて算出される第1操作量(Fb)を修正するための先行量(Ff)を算出するよう構成された先行量算出部(3)と、を備え、
 前記先行量算出部(3)は、
 前記制御量を変化させるイベントの発生時に、前記イベントの継続時間(t)に応じた規定先行量(Qa)を算出するよう構成された第1算出部(31)と、
 前記予測変動値(V)および前記偏差(D)に基づいて、前記規定先行量(Qa)を修正するための先行指標(Qb)を算出するよう構成された第2算出部(4)と、
 前記規定先行量(Qa)および前記先行指標(Qb)に基づいて、前記先行量(Ff)を算出するよう構成された第3算出部(5)と、を有する。
 上記(1)の構成によれば、プラント(7)の制御量(制御対象の状態量)のフィードバック制御により求まる第1操作量(Fb)を修正する先行量(Ff)(フィードフォワード操作量)を、過去のプラント(7)の運転(実運転、試運転)等で得られた運転データなどからプロセスパラメータ(P)と制御量の変動値との関係性を学習した学習モデル(M)により得られる制御量の変動値の予測結果(V)に基づいて算出する。例えばボイラ(71)を有するプラント(7)においては、除灰装置の稼働や、バーナの点火や消火といったイベントの発生は、主蒸気温度や主蒸気圧力、再熱蒸気温度、NOx濃度などの状態量(制御量)のフィードバック制御にとって外乱となるが、この外乱の制御量への影響を抑制するための先行量(Ff)を上記の学習モデル(M)に基づいて得られる予測変動値(V)に基づいて算出する。
 具体的には、制御理論や経験などに基づいて作成した関数などを用いて、例えば除灰装置の作動時間(イベントの継続時間(t))に応じた規定先行量(Qa)を算出する。そして、この規定先行量(Qa)を、予測変動値(V)、および制御量の計測値(PV)と制御設定値(SV)との偏差(D)に基づいて算出した先行指標(Qb)を用いてより適切なものに修正することで、先行量(Ff)を算出する。つまり、学習モデル(M)から求めた予測変動値(V)に基づいて直接的に先行量(Ff)を算出するのではなく、この予測変動値(V)に基づいて算出した算出値を、その直前に求めた先行量(Ff)を対象操作端(9)に適用した結果によりフィードバック制御的に調整すると共に、このようにして算出した先行指標(Qb)で、上記の規定先行量(Qa)を修正するようにして、先行量(Ff)を算出する。
 上述したような学習モデル(M)を用いたとしても予測変動値(V)と実際の偏差(D)との間には誤差が生じる可能性があり、予測変動値(V)からのみ先行指標(Qb)を求める場合には誤差がある状態で予測変動値(V)を求めることになってしまうが、予測変動値(V)および実際の偏差(D)に基づいて先行指標(Qb)を求めることで、予測変動値(V)と実際の偏差(D)との誤差が小さくされたことに相当するような先行指標(Qb)を求めることができる。よって、プラント(7)の運転状態やその経時変換に応じたより適切な先行量(Ff)を求めることができ、対象操作端(9)の操作量(MV)の適正化を通して、制御量の演算をより適切に行うことができる。
(2)幾つかの実施形態では、上記(1)の構成において、
 前記第2算出部(4)は、
 前記予測変動値(V)に応じた規定先行指標値(Q1)を算出する規定先行指標算出部(42)と、
 前記偏差(D)に基づいて、前記規定先行指標値(Q1)を調整するための調整値(Q2)を算出する評価部(43)と、
 前記規定先行指標値(Q1)および前記調整値(Q2)に基づいて、前記先行指標(Qb)を差出する先行指標算出部(44)と、を有する。
 上記(2)の構成によれば、予測変動値(V)に応じて算出される修正前指標を、偏差(D)に基づいて算出される調整値(Q2)で調整することで、先行指標(Qb)を算出する。これによって、より適切な先行量(Ff)を求めることができる。
(3)幾つかの実施形態では、上記(1)~(2)の構成において、
 前記プロセスパラメータ(P)は、1以上の状態量の計測値と、前記プラント(7)に設置された1以上の操作端の操作量と、前記先行指標(Qb)と、を含む。
 上記(3)の構成によれば、予測時における状態量の計測値と、操作端の操作量と、先行指標(Qb)(前回値)とを含む。状態量の計測値と操作端の操作量とを含むことで、プロセスパラメータ(P)によって、プラント(7)の運転状態を適切に特定することができる。また、先行指標(Qb)を含むことで、学習モデル(M)が前回の規定タイミングにおける先行指標(Qb)を考慮した予測変動値(V)を出力することができ、学習データの収集時間を含めた学習時間や、この学習データを得るための例えばボイラ(71)の試運転時間を短縮することができる。
(4)幾つかの実施形態では、上記(1)~(3)の構成において、
 前記プラント(7)が有するボイラ(71)の負荷値(L)を取得する負荷取得部(32)をさらに備え、
 前記変動予測部(2)は、前記ボイラ(71)の複数の負荷帯に応じて設けられた複数の前記学習モデル(M)の各々の前記予測変動値(V)に基づいて、前記負荷値(L)に応じた前記予測変動値(V)を出力するように構成されており、
 前記先行量算出部(3)は、前記複数の負荷帯に応じてそれぞれ設けられた前記第1算出部(31)、前記第2算出部(4)および前記第3算出部(5)を有し、前記負荷値(L)に応じた前記先行量(Ff)を算出するよう構成されている。
 上記(4)の構成によれば、先行量(Ff)は、その算出時におけるプラント(7)(ボイラ(71))の負荷を考慮して算出される。これによって、負荷に応じた適切な先行量(Ff)を求めることができる。
1    制御装置
10   制御プログラム
11   プロセッサ
12   記憶部
2    変動予測部
3    先行量算出部
31   第1算出部
32   負荷取得部
4    第2算出部
41   保持部
42   規定先行指標算出部
43   評価部
43a  規定調整値算出部
43b  出力切替部
43c  積算値算出部
43d  調整値算出部
44   先行指標算出部
45   切替部
5    第3算出部
52   切替部
6    フィードバック制御部
61   偏差算出部
62   コントローラ
63   加算部
7    プラント
7L   冷却媒体供給管
7c   温度低減器
7p   蒸気管
7s   圧力センサ
7t   ボイラ出口温度計測センサ
7v   スプレイ弁
7w   冷却媒体
71   ボイラ
71f  火炉
71p  炉壁蒸発管
72   燃焼装置
72v  燃料供給装置
73   煙道
74   燃焼バーナ
75   微粉炭供給管
76   粉砕機
77a  空気ダクト
77b  分岐部
77c  分岐ダクト
77d  ダンパ
77e  風箱
77p  ポート
78   送風機
79   熱交換器
79a  過熱器
79b  再熱器
79c  節炭器
8    排ガス通路
81   脱硝装置
81p  アンモニア供給配管
81s  NOxセンサ
81v  アンモニア流量制御弁
8w   アンモニア
82   エアヒータ
83   煤塵処理装置
84   誘引送風機
85   煙突
86   再循環ファン
87   ガス再循環ファン入口ダンパ
9    対象操作端
T    蒸気タービン
G    発電機
P    プロセスパラメータ
M    学習モデル
V    予測変動値
Ff   先行量
Fb   第1操作量
PV   制御量の計測値
SV   制御設定値
MV   操作量
D    偏差
E    発生信号
I    積算値
L    現負荷値
Qa   規定先行量
Qb   先行指標
Q1   規定先行指標値
Q2   調整値
Q3   規定調整値
t    イベントの継続時間

Claims (4)

  1.  プラントの制御量を制御可能な対象操作端の操作量を規定タイミング毎に決定する制御装置であって、
     前記プラントの運転状態を特定するためのプロセスパラメータと前記制御量の変動値との関係性を学習した学習モデルを用いて、前記プロセスパラメータの入力に応じた前記制御量の予測変動値を出力するよう構成された変動予測部と、
     前記予測変動値に基づいて、前記制御量の計測値と制御設定値との偏差に応じて算出される第1操作量を修正するための先行量を算出するよう構成された先行量算出部と、を備え、
     前記先行量算出部は、
     前記制御量を変化させるイベントの発生時に、前記イベントの継続時間に応じた規定先行量を算出するよう構成された第1算出部と、
     前記予測変動値および前記偏差に基づいて、前記規定先行量を修正するための先行指標を算出するよう構成された第2算出部と、
     前記規定先行量および前記先行指標に基づいて、前記先行量を算出するよう構成された第3算出部と、を有する制御装置。
  2.  前記第2算出部は、
     前記予測変動値に応じた規定先行指標値を算出する規定先行指標算出部と、
     前記偏差に基づいて、前記規定先行指標値を調整するための調整値を算出する評価部と、
     前記規定先行指標値および前記調整値に基づいて、前記先行指標を差出する先行指標算出部と、を有する請求項1に記載の制御装置。
  3.  前記プロセスパラメータは、1以上の状態量の計測値と、前記プラントに設置された1以上の操作端の操作量と、前記先行指標と、を含む請求項1または2に記載の制御装置。
  4.  前記プラントが有するボイラの負荷値を取得する負荷取得部をさらに備え、
     前記変動予測部は、前記ボイラの複数の負荷帯に応じて設けられた複数の前記学習モデルの各々の前記予測変動値に基づいて、前記負荷値に応じた前記予測変動値を出力するように構成されており、
     前記先行量算出部は、前記複数の負荷帯に応じてそれぞれ設けられた前記第1算出部、前記第2算出部および前記第3算出部を有し、前記負荷値に応じた前記先行量を算出するよう構成されている請求項1~3のいずれか1項に記載の制御装置。
PCT/JP2021/016135 2020-05-01 2021-04-21 制御装置 WO2021220901A1 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2020081362A JP7438014B2 (ja) 2020-05-01 2020-05-01 制御装置
JP2020-081362 2020-05-01

Publications (1)

Publication Number Publication Date
WO2021220901A1 true WO2021220901A1 (ja) 2021-11-04

Family

ID=78300568

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2021/016135 WO2021220901A1 (ja) 2020-05-01 2021-04-21 制御装置

Country Status (2)

Country Link
JP (1) JP7438014B2 (ja)
WO (1) WO2021220901A1 (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7348359B1 (ja) 2022-07-20 2023-09-20 三菱重工業株式会社 ボイラ用のコントローラ、ボイラシステム、および、ボイラ制御プログラム
JP2024124174A (ja) * 2023-03-02 2024-09-12 オムロン株式会社 制御対象に対するフィードフォワード制御に関する機械学習のための装置、教師データの生成プログラム、および教師データの生成方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01503498A (ja) * 1987-03-30 1989-11-22 ザ フォックスボロ カンパニー 多変数適応フィードフォワード制御のためのプロセス制御器及びプロセスの多変数適応制御方法
JPH05204407A (ja) * 1992-01-28 1993-08-13 Matsushita Electric Works Ltd プロセスの特性予測方法ならびにその予測方法を用いたプロセスの監視方法およびプロセスの制御方法
JPH11242503A (ja) * 1998-02-25 1999-09-07 Hitachi Ltd プラント運転制御支援システム
JP2002157003A (ja) * 2000-11-17 2002-05-31 Fuji Electric Co Ltd 多変数モデル予測制御装置、方法、及びその記憶媒体
WO2014208228A1 (ja) * 2013-06-26 2014-12-31 株式会社日立製作所 ボイラ制御装置
JP2020036387A (ja) * 2018-08-27 2020-03-05 株式会社ノーリツ 制御装置およびこれを備えた燃焼装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6032847B2 (ja) 2013-03-18 2016-11-30 三菱日立パワーシステムズ株式会社 蒸気タービン起動制御システム

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01503498A (ja) * 1987-03-30 1989-11-22 ザ フォックスボロ カンパニー 多変数適応フィードフォワード制御のためのプロセス制御器及びプロセスの多変数適応制御方法
JPH05204407A (ja) * 1992-01-28 1993-08-13 Matsushita Electric Works Ltd プロセスの特性予測方法ならびにその予測方法を用いたプロセスの監視方法およびプロセスの制御方法
JPH11242503A (ja) * 1998-02-25 1999-09-07 Hitachi Ltd プラント運転制御支援システム
JP2002157003A (ja) * 2000-11-17 2002-05-31 Fuji Electric Co Ltd 多変数モデル予測制御装置、方法、及びその記憶媒体
WO2014208228A1 (ja) * 2013-06-26 2014-12-31 株式会社日立製作所 ボイラ制御装置
JP2020036387A (ja) * 2018-08-27 2020-03-05 株式会社ノーリツ 制御装置およびこれを備えた燃焼装置

Also Published As

Publication number Publication date
JP7438014B2 (ja) 2024-02-26
JP2021176047A (ja) 2021-11-04

Similar Documents

Publication Publication Date Title
CA2733107C (en) Process temperature control in oxy/fuel combustion system
US8515582B2 (en) Control system for operation of a fossil fuel power generating unit
WO2021220901A1 (ja) 制御装置
EP2107220B1 (en) Steam temperature control in a boiler system using reheater variables
US8973535B2 (en) Steam-generator temperature control and optimization
TWI705316B (zh) 鍋爐之運轉支援裝置、鍋爐之運轉支援方法、及鍋爐之學習模型之作成方法
US9841185B2 (en) Steam temperature control using model-based temperature balancing
US9217565B2 (en) Dynamic matrix control of steam temperature with prevention of saturated steam entry into superheater
CA2747950A1 (en) Dynamic tuning of dynamic matrix control of steam temperature
Gilman et al. Boiler control systems engineering
EP2990617B1 (en) Apparatus and method for controlling at least one operational parameter of a plant
JP7261113B2 (ja) ボイラの制御装置、ボイラシステム、発電プラント、及びボイラの制御方法
CN110094747A (zh) 吹灰器运转控制装置、吹灰器运转控制方法及燃烧系统
TWI735042B (zh) 廠房的運轉支援裝置及廠房的運轉支援方法
Black et al. Challenges when converting coal-fired boilers to natural gas
JPS6375406A (ja) 加圧流動燃焼法
WO2024101074A1 (ja) 脱硝制御装置及び脱硝装置
Braun et al. Dynamic Simulation of a Sub-Critical Coal Fired Power Plant
Huang et al. Control strategy for main steam pressure of combustion system of pulverized coal boiler
GB802032A (en) A steam generator and method of operating the same
Ganapathy Understand boiler performance characteristics
Kellerman et al. Modelling of boiler fireside control in flownex software environment
Leithner Steam generators and steam distribution networks
Smuts Improving Steam Temperature Control with Neural Networks
Bittanti et al. Pressurized oxy-coal mild combustion for clean-coal technology

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 21795966

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 21795966

Country of ref document: EP

Kind code of ref document: A1