Combustion self-optimization control method, system and equipment for header system circulating fluidized bed boiler and medium

Abstract

The invention discloses a header system circulating fluidized bed boiler combustion self-optimization control method, system, equipment and medium, and relates to the technical field of automatic control, and the method comprises the following steps: S01, collecting the steam pressure of a high-pressure header and time domain data of process values of coal quantity, steam flow, bed temperature, primary air quantity, oxygen quantity and secondary air quantity of N boilers, and constructing a process value time sequence, calculating process value stability based on the time series; step S02, calculating optimized basic values and constraint values of process values of coal feeding quantity, steam flow, bed temperature, primary air quantity, oxygen quantity and secondary air quantity of the N boilers at the current moment; and S03, calculating the oxygen amount of the N boilers at the current moment, the upper and lower optimization weights of the ratio of the primary air volume to the secondary air volume and the like. According to the method, the stability of the combustion process can be guaranteed, the ton steam coal consumption can be reduced, and the benefits of the production process are maximized.

Description

Technical Field

[0001] This invention relates to the field of automatic control technology, and in particular to a method, system, equipment and medium for self-optimization control of combustion in a main-pipe circulating fluidized bed boiler. Background Technology

[0002] Steam supply from multiple boilers is delivered to a single steam header, from which all steam-consuming units obtain steam. This header operation method facilitates a stable steam supply from the boilers to the steam-consuming units. Boilers can be classified according to their combustion method, such as circulating fluidized bed boilers and direct-fired pulverized coal boilers. Circulating fluidized bed boilers are widely used in combined heat and power (CHP) and chemical production due to their wide fuel adaptability and low emissions. In the combustion process of a circulating fluidized bed boiler, bed temperature, primary air volume, oxygen content, and secondary air volume are crucial control parameters. These parameters affect both the stability and economic efficiency of combustion. System stability can be ensured by constructing an automatic control loop that controls bed temperature with primary air volume and oxygen with secondary air volume. Economic efficiency requires optimizing the oxygen setpoint and the ratio of primary to secondary air volume to maximize combustion economic benefits, i.e., reducing coal consumption per ton of steam. Since equipment performance and fuel characteristics change over time, failure to optimize and adjust these factors will reduce the economic efficiency of the production process. Therefore, it is necessary to design a self-optimizing control method to automatically find the oxygen and primary air volume that reduce the consumption per ton of steam and coal when the combustion process is stable. Summary of the Invention

[0003] To address the aforementioned issues, this invention provides a method, system, equipment, and medium for self-optimizing combustion control of a mother-pipe circulating fluidized bed boiler. The aim is to maximize production efficiency by automatically identifying and reducing the oxygen and primary air volume per ton of steam and coal consumption.

[0004] A method and system for self-optimizing combustion control of a circulating fluidized bed boiler with a main control pipe includes the following main steps:

[0005] Step S01: Collect time-domain data of the process values ​​of high-pressure main pipe steam pressure and N boilers coal quantity, steam flow rate, bed temperature, primary air volume, oxygen quantity, and secondary air volume, and construct a time series of process values. Calculate the stability of process values ​​based on the time series.

[0006] Step S02: Calculate the optimized baseline and constraint values ​​for the process values ​​of coal feed rate, steam flow rate, bed temperature, primary air volume, oxygen volume, and secondary air volume of N boilers at the current time.

[0007] Step S03: Calculate the upper and lower optimization weights of the ratio of oxygen content, primary air volume and secondary air volume of N boilers at the current time.

[0008] Step S04: Construct oxygen and primary air volume optimization controllers for N boilers, determine the optimization direction using the ton-coal consumption rotation factor, and calculate the optimized setpoints for oxygen and primary air volume that reduce ton-coal consumption at the current time by performing optimization within the constraints.

[0009] Step S05: Compare the setpoints of oxygen and primary air volume of N boilers in the previous period with the time series of process values ​​at the current time to determine whether to optimize the setpoint output. If the condition is met, the optimized setpoint can be added to the current setpoint; otherwise, the optimized setpoint will not be added to the current setpoint.

[0010] Further, step S01 specifically includes the following steps:

[0011] Step S101: Obtain the process values ​​of the K high-pressure main pipe steam pressure, N boiler coal quantity, steam flow rate, bed temperature, primary air volume, oxygen quantity, and secondary air volume before the current time, and construct them into time series data.

[0012] Process value time series: (Formula 1) in, PV represents a process-value time series. This represents the sampled value of the first process value. This represents the sampled value of the Kth process value. Indicates superscript, when "main" represents the high-voltage main pipe. Let BN represent the Nth boiler, where BN is the boiler-number. Indicates a subscript, when STPR represents high-pressure steam pressure, where STPR stands for steam-pressure. Let coal represent the amount of coal, when STFL represents steam flow rate. BDT stands for bed-temperature. air1 represents the primary air volume, when O2 represents the amount of oxygen. air2 indicates the secondary air volume. The following... and Both have the same meaning.

[0013] Step S102: Calculate the stability of the high-pressure main pipe steam pressure, coal feed rate of N boilers, steam flow rate, bed temperature, primary air volume, oxygen content, and secondary air volume based on the process value mean, oscillation number, and deviation.

[0014] For the process value time series established in "step S101", the mean of the process value, the rate of change and the number of oscillations are calculated first, and then the stability is calculated.

[0015] (1) Calculation of the mean of process values: (Formula 2) in, AVE represents the average value of the process. This represents a time series of process values, where K represents the number of samples and m represents the current period. (Subscript) and The meaning is the same as described in "Step S101".

[0016] (2) Calculation of process value change rate: (Formula 3) in, RT represents the rate of change of the process value. This represents a time series of process values, where K represents the number of samples. This represents the average value of the process, and m represents the current period. and The meaning is the same as described in "Step S101".

[0017] (3) Calculation of the number of oscillations in the process value like ,but ,otherwise Finally, it was concluded that = (Formula 4) in, Represents a time series of process values. This represents the average value of the process values. This represents an intermediate variable used in calculating the number of oscillations during the process. , This indicates the number of oscillations in the process value. OSCI stands for oscillating, K represents the number of samples, and m represents the current period. (Subscript) and The meaning is the same as described in "Step S101".

[0018] (4) Calculation of process value deviation: (Formula 5) in, DEV stands for deviation from process value. Represents a time series of process values. This represents the average value of the process, K represents the number of samples, and m represents the current period. (Subscript) and The meaning is the same as described in "Step S101".

[0019] (5) Calculation of process value stability: (Formula 6) in, STA stands for stability, representing process value stability. Indicates the rate of change of process value. Indicates the number of oscillations in the process value. This indicates the preset value for the number of oscillations in the process value. For set-oscillating, , , , , , , , This indicates the deviation of the process value, `max` indicates finding the maximum value, `K` indicates the number of samples, and `m` indicates the current period. (Subscript) and The meaning is the same as described in "Step S101".

[0020] Furthermore, step S02 specifically includes the following steps:

[0021] Step S201: Calculate the optimized baseline values ​​for the process values ​​of coal feed rate, steam flow rate, bed temperature, primary air volume, oxygen volume, and secondary air volume of N boilers at the current moment.

[0022] To smooth the mean of process values ​​and avoid the impact of noise and abnormal disturbances on the optimization, the specific calculation method for the optimization base value is as follows: = (Formula 7) in, This represents the base value for process value optimization. For base-process-value, This represents the average value of the process, where m represents the current period, m-1 represents the previous period, m-2 represents the previous two periods, and m-3 represents the previous three periods. (Subscript) and The meaning is the same as described in "Step S101".

[0023] Step S202: Calculate the optimized constraint values ​​for the process values ​​of coal feed rate, steam flow rate, bed temperature, primary air volume, oxygen volume, and secondary air volume of N boilers at the current time.

[0024] By utilizing the set process value constraints and the stability of the process values, the upper and lower constraints of each process value optimization constraint are calculated comprehensively, serving as the constraints for optimization.

[0025] (1) Calculation of upper constraint limits for process value optimization: (Formula 8) in, Upper constraint limit for process value optimization For optimal-high-limit, SETHL indicates the upper limit of the process value constraint. For set-high-limit, Indicates process value stability, where m represents the current period. Subscript and The meaning is the same as described in "Step S101".

[0026] The method for calculating the upper limit of the primary air volume to secondary air volume ratio is as follows: (Formula 9) in, This represents the upper constraint limit for optimizing the ratio of primary air volume to secondary air volume. This indicates the upper limit of the set primary air volume to secondary air volume ratio constraint. This indicates the stability of the primary air volume process value. This indicates the stability of the secondary air volume process value, where m represents the current period. (Subscript) The meaning is the same as described in "Step S101".

[0027] (2) Calculation of lower constraint limits for process value optimization: (Formula 10) in, This represents the lower constraint limit for process value optimization. SETLL, meaning optimal-low-limit, indicates the lower limit of the process value constraint. For set-low-limit, Indicates process value stability, where m represents the current period. Subscript and The meaning is the same as described in "Step S101".

[0028] The method for calculating the lower constraint limit of the ratio of primary air volume to secondary air volume is as follows: (Formula 11) in, This represents the lower constraint limit for optimizing the ratio of primary air volume to secondary air volume. This indicates the lower limit of the set constraint on the ratio of primary air volume to secondary air volume. This indicates the stability of the primary air volume process value. This indicates the stability of the secondary air volume process value, where m represents the current period. (Subscript) The meaning is the same as described in "Step S101".

[0029] Furthermore, step S03 specifically includes the following steps:

[0030] Step S301: Calculate the optimized weight of oxygen based on the process value time series constructed in step S01.

[0031] The time series of oxygen process values ​​is weighted and summed with the optimization constraints to calculate the optimization weight of oxygen. The larger the weight value, the larger the optimization output, and vice versa.

[0032] (1) Calculation of the optimization weight for oxygen content: (Formula 12) in, To optimize the weighting of oxygen content in the Nth boiler, For high-weight, This represents the Nth boiler. Indicates oxygen content. This represents the upper constraint limit for oxygen optimization of the Nth boiler, where m represents the current period. This represents the oxygen content process value of the Nth boiler, and K represents the number of samples.

[0033] (2) Calculation of optimization weights under oxygen content: (Formula 13) in, To optimize the weights for the oxygen content of the Nth boiler, For low-weight, This represents the Nth boiler. Indicates oxygen content. This represents the lower constraint limit for oxygen optimization in the Nth boiler. This represents the oxygen content process value of the Nth boiler, and K represents the number of samples.

[0034] Step S302: Calculate the optimized weight of the ratio of primary air volume to secondary air volume based on the process value time series constructed in step S01.

[0035] The optimization weight of the primary air volume is calculated by weighting and summing the time series ratios of the primary and secondary air volume process values ​​with the optimization constraints. The larger the weight value, the larger the optimization output, and vice versa.

[0036] (1) Optimize the weighting of the ratio of primary air volume to secondary air volume: (Formula 14) in, Optimize the weighting of the ratio of primary air volume to secondary air volume for the Nth boiler. For high-weight, This indicates the Nth boiler, and AIRR represents the ratio of primary air volume to secondary air volume. AIRR is... - , This represents the upper constraint limit for optimizing the ratio of primary air volume to secondary air volume of the Nth boiler. This represents the process value of the primary air volume of the Nth boiler. This represents the process value of the secondary air volume of the Nth boiler, where m represents the current period and K represents the number of samples.

[0037] (2) Optimization weights based on the ratio of primary air volume to secondary air volume: (Formula 15) in, The optimal weights are determined for the ratio of primary air volume to secondary air volume of the Nth boiler. For low-weight, This represents the Nth boiler. This represents the ratio of primary air volume to secondary air volume. This represents the lower constraint limit for optimizing the ratio of primary air volume to secondary air volume of the Nth boiler. This represents the process value of the primary air volume of the Nth boiler. This represents the process value of the secondary air volume of the Nth boiler, where m represents the current period and K represents the number of samples.

[0038] Furthermore, step S04 specifically includes the following steps:

[0039] Step S401: Calculate the incremental factor, incremental rotation factor, and number of consecutive flips for the ton of steam coal consumption.

[0040] (1) Calculation of incremental factor for coal consumption per ton of steam: (Formula 16) in, This represents the incremental factor of coal consumption per ton of steam for the Nth boiler. For factor, This indicates the Nth boiler, where COST represents the coal-to-steam ratio. This represents the baseline value for optimizing the coal quantity process of the Nth boiler. This represents the optimized baseline value of the steam flow process for the Nth boiler, where m represents the current cycle and m-1 represents the previous cycle.

[0041] (2) Calculation of the rotation factor for the incremental coal consumption per ton of steam: (Formula 17) in, This represents the rotational factor for the incremental coal consumption per ton of steam in the Nth boiler. For optimal-factor, This represents the incremental factor of coal consumption per ton of steam for the Nth boiler. Let N represent the Nth boiler, COST represent the coal-to-steam ratio, m represent the current cycle, and m-1 represent the previous cycle.

[0042] (3) Continuous flip counting calculation: If , +1 (Formula 18), otherwise , ,in, This represents the rotational factor for the incremental coal consumption per ton of steam in the Nth boiler. This indicates the Nth boiler, and COST represents the coal-to-steam ratio. This represents the number of times the ton of steam kerosene is turned over in the Nth boiler. For count-factor, m represents the current period, and m-1 represents the previous period.

[0043] Step S402: Construct an oxygen optimization controller and perform oxygen optimization calculations.

[0044] when At that time, the oxygen content optimization calculation is as follows: If , At that time, the oxygen optimization setpoint is calculated as follows: (Formula 19), otherwise

[0045] in, This represents the rotational factor for the incremental coal consumption per ton of steam in the Nth boiler. This indicates the Nth boiler, and COST represents the coal-to-steam ratio. This represents the number of times the ton of steam kerosene is turned over in the Nth boiler. This represents the preset value for determining the number of times the ton of steam kerosene is rotated in the Nth boiler. Greater than or equal to Stop optimizing. This represents the oxygen setting value for the Nth boiler. This represents the upper limit of oxygen optimization for the Nth boiler. This represents the optimized oxygen setting value for the Nth boiler. This represents the baseline value for optimizing the coal quantity process of the Nth boiler. This represents the optimized baseline value for the steam flow rate process of the Nth boiler. This represents the optimization weight of oxygen quantity for the Nth boiler, where m represents the current period and m-1 represents the previous period.

[0046] when At that time, the oxygen content optimization calculation is as follows: If , At that time, the oxygen optimization setpoint is calculated as follows: ,otherwise ;

[0047] in, This represents the rotational factor for the incremental coal consumption per ton of steam in the Nth boiler. This indicates the Nth boiler, and COST represents the coal-to-steam ratio. This represents the number of times the ton of steam kerosene is turned over in the Nth boiler. This represents the preset value for determining the number of times the ton of steam kerosene is rotated in the Nth boiler. Greater than or equal to Stop optimizing. This represents the oxygen setting value for the Nth boiler. This represents the lower limit of oxygen optimization for the Nth boiler. This represents the optimized oxygen setting value for the Nth boiler. This represents the baseline value for optimizing the coal quantity process of the Nth boiler. This represents the optimized baseline value for the steam flow rate process of the Nth boiler. This represents the optimization weight for the oxygen content of the Nth boiler, where m represents the current period and m-1 represents the previous period.

[0048] Step S403: Construct a primary airflow optimization controller and perform a primary airflow optimization calculation.

[0049] when When the primary air volume is optimized, the calculation method is as follows: If , , , , Primary air volume optimization setpoint calculation: ,otherwise ; in, This represents the rotational factor for the incremental coal consumption per ton of steam in the Nth boiler. This indicates the Nth boiler, and COST represents the coal-to-steam ratio. This represents the number of times the ton of steam kerosene is turned over in the Nth boiler. This represents the preset value for determining the number of times the ton of steam kerosene is rotated in the Nth boiler. Greater than or equal to Stop optimizing. This represents the optimized baseline value for the primary air volume process of the Nth boiler. This represents the optimized baseline value for the secondary air volume process of the Nth boiler. This represents the upper limit of the optimization ratio of primary air volume to secondary air volume for the Nth boiler. This represents the optimized baseline value for the bed temperature process of the Nth boiler. This represents the lower limit of the optimal bed temperature for the Nth boiler. This represents the upper limit of the primary air volume optimization for the Nth boiler. This represents the lower limit for optimizing the secondary air volume of the Nth boiler. This represents the optimized setpoint for the primary air volume of the Nth boiler. This represents the baseline value for optimizing the coal quantity process of the Nth boiler. This represents the optimized baseline value for the steam flow rate process of the Nth boiler. Optimize the weighting of the ratio of primary air volume to secondary air volume for the Nth boiler. This represents the setpoint for the secondary air volume of the Nth boiler, where m represents the current cycle and m-1 represents the previous cycle.

[0050] when When the primary air volume is optimized, the calculation method is as follows: If , , , , Primary air volume optimization setpoint calculation: ,otherwise ; in, This represents the rotational factor for the incremental coal consumption per ton of steam in the Nth boiler. This indicates the Nth boiler, and COST represents the coal-to-steam ratio. This represents the number of times the ton of steam kerosene is turned over in the Nth boiler. This represents the preset value for determining the number of times the ton of steam kerosene is rotated in the Nth boiler. Greater than or equal to Stop optimizing. This represents the optimized baseline value for the primary air volume process of the Nth boiler. This represents the optimized baseline value for the secondary air volume process of the Nth boiler. This represents the lower limit of the optimization ratio of primary air volume to secondary air volume for the Nth boiler. This represents the optimized baseline value for the bed temperature process of the Nth boiler. This represents the upper limit of the optimized bed temperature for the Nth boiler. This represents the lower limit of the primary air volume optimization for the Nth boiler. This represents the upper limit of the optimized secondary air volume for the Nth boiler. This represents the optimized setpoint for the primary air volume of the Nth boiler. This represents the baseline value for optimizing the coal quantity process of the Nth boiler. This represents the optimized baseline value for the steam flow rate process of the Nth boiler. The optimal weights are determined for the ratio of primary air volume to secondary air volume of the Nth boiler. This represents the setpoint for the secondary air volume of the Nth boiler, where m represents the current cycle and m-1 represents the previous cycle.

[0051] Furthermore, step S05 specifically includes the following steps:

[0052] Step S501: Compare the current oxygen process value time series with the previous cycle set value. If the condition is met, the optimized set value can be added to the previous cycle oxygen set value and used as the current cycle oxygen set value; otherwise, the set value remains unchanged.

[0053] The time series of oxygen quantity process values ​​is analyzed, and when each time series value is... When (i=1, ..., K), the oxygen optimization setting value is accumulated to the oxygen setting value of the previous cycle and used as the oxygen setting value of the current cycle; otherwise, the oxygen setting value of the current cycle remains unchanged.

[0054] (Formula 23) in, This represents the time series of oxygen quantity process values ​​for the Nth boiler. This represents the oxygen setting value for the Nth boiler. This represents the oxygen optimization setting value for the Nth boiler, m represents the current cycle, m-1 represents the previous cycle, and K represents the number of samples.

[0055] Step S502: Compare the time series of the primary air volume process value at the current time with the set value of the previous period. If the condition is met, the optimized set value can be added to the primary air volume set value of the previous period and used as the primary air volume set value of the current period. Otherwise, the set value remains unchanged.

[0056] The time series of primary air volume values ​​is analyzed, and when each value in the time series... When i=1, ..., K, the primary air volume optimization setting value is accumulated to the primary air volume setting value of the previous cycle and used as the primary air volume setting value of the current cycle; otherwise, the primary air volume setting value of the current cycle remains unchanged.

[0057] (Formula 24) in, This represents the time series of primary air volume values ​​for the Nth boiler. This represents the primary air volume setpoint for the Nth boiler. This represents the optimized setpoint for the primary air volume of the Nth boiler, where m represents the current cycle, m-1 represents the previous cycle, and K represents the number of samples.

[0058] A self-optimizing control system for combustion in a circulating fluidized bed boiler with a main pipe includes a data acquisition module for acquiring process data of main pipe pressure and coal quantity, steam flow rate, bed temperature, primary air volume, oxygen quantity, and secondary air volume of N boilers.

[0059] The control module is used to implement the self-optimization control method for combustion in the main pipe circulating fluidized bed boiler, and to write back the calculation results to the field control system.

[0060] The data storage module is used to store the data generated by the data acquisition module and the control module.

[0061] The data is displayed by the data storage module.

[0062] A computing device includes a memory, a processor, and a program stored in the memory and running on the processor, wherein the processor executes the program to implement the steps of the self-optimizing combustion control method for a mother-pipe circulating fluidized bed boiler.

[0063] A computer-readable storage medium storing instructions that, when executed on a terminal device, cause the terminal device to perform the steps of the self-optimizing combustion control method for a main-pipe circulating fluidized bed boiler.

[0064] Compared with existing technologies, this technical solution has the following characteristics:

[0065] 1. Perform stability calculations for high-pressure header steam pressure, boiler coal quantity, steam flow rate, bed temperature, primary air volume, oxygen quantity, and secondary air volume every cycle, and incorporate stability into the optimization calculations to ensure that optimization is carried out under stable operation during the production process.

[0066] 2. Use the stability of the main pipe pressure as the upper and lower limits of the process value optimization constraint for each boiler to ensure the self-optimization stability under the operation of the main pipe.

[0067] 3. Based on the upper and lower limits of the optimization constraints, dynamically calculate the weights of the ratio of oxygen, primary air volume, and secondary air volume, and output the optimization setpoints according to the weights to reduce the disturbance of the optimization setpoints to the production process.

[0068] 4. By using a self-optimization method to reduce coal consumption per ton of steam, the oxygen content and primary air volume of the boiler are optimized to maximize the efficiency of the production process. Attached Figure Description

[0069] Figure 1 This is a schematic diagram of the combustion process of the circulating fluidized bed boiler with a main control system according to the present invention;

[0070] Figure 2 This is a schematic diagram of the self-optimizing combustion control system for a circulating fluidized bed boiler according to the present invention. Detailed Implementation

[0071] The specific implementation process of the present invention will be described in detail below with reference to the accompanying drawings and specific embodiments.

[0072] like Figure 1-2As shown, this invention provides a method, system, equipment, and medium for self-optimizing combustion control of a main-pipe circulating fluidized bed boiler. The method includes the following steps: Step S01: Collect time-domain data of the high-pressure main pipe steam pressure and the process values ​​of coal feed, steam flow, bed temperature, primary air volume, oxygen content, and secondary air volume of N boilers, and construct a time series of process values. Calculate the stability of the process values ​​based on the time series. Step S02: Calculate the optimized baseline and constraint values ​​of the process values ​​of coal feed, steam flow, bed temperature, primary air volume, oxygen content, and secondary air volume of N boilers at the current time. Step S03: Calculate the upper and lower optimization weights of the oxygen content and the ratio of primary air volume to secondary air volume of N boilers at the current time. Step S04: Construct an optimization controller for the oxygen content and primary air volume of N boilers, determine the optimization direction using a ton-coal consumption rotation factor, and calculate the optimized setpoints for oxygen content and primary air volume that reduce ton-coal consumption at the current time by performing optimization within the constraints. Step S05: Compare the setpoints of oxygen and primary air volume of N boilers in the previous period with the time series of process values ​​at the current time to determine whether to optimize the setpoint output. If the condition is met, the optimized setpoint can be added to the current setpoint; otherwise, the optimized setpoint will not be added to the current setpoint.

[0073] This invention's method includes performing stability calculations for high-pressure header steam pressure, boiler coal quantity, steam flow rate, bed temperature, primary air volume, oxygen quantity, and secondary air volume every cycle. Stability is incorporated into the optimization calculations to ensure optimization occurs under stable production conditions. The stability of the header pressure is used as the upper and lower limits of the optimization constraints for each boiler process value, ensuring self-optimization stability under header operation. Based on the upper and lower limits of the optimization constraints, the weights of the oxygen quantity, primary air volume, and secondary air volume ratio are dynamically calculated, and the optimization setpoints are output according to these weights, reducing the disturbance of the optimization setpoints to the production process. This invention enables a self-optimization method to reduce coal consumption per ton of steam by optimizing the boiler's oxygen quantity and primary air volume, maximizing production efficiency.

[0074] This invention also provides a self-optimizing combustion control system for a main-pipe circulating fluidized bed boiler, comprising: a data acquisition module for acquiring process data on main pipe pressure and coal quantity, steam flow rate, bed temperature, primary air volume, oxygen quantity, and secondary air volume of N boilers; a control module for implementing the self-optimizing combustion control method of the main-pipe circulating fluidized bed boiler and writing back the calculation results to the field control system, which is an existing DCS (Distributed Control System); a data storage module for storing the data generated by the data acquisition module and the control module; and a data display module for displaying the data from the data storage module.

[0075] The system of this invention communicates with the field control system via OPC (OLE for Process Control) to acquire and write back data. The acquired data can be stored in databases such as SQLite and MySQL.

[0076] The present invention also provides a computing device, including a memory, a processor, and a program stored in the memory and running on the processor, wherein the processor executes the program to implement the steps of the self-optimizing combustion control method for the main-pipe circulating fluidized bed boiler.

[0077] The present invention also provides a computer-readable storage medium storing instructions that, when executed on a terminal device, cause the terminal device to perform the steps of the self-optimizing combustion control method for a mother-pipe circulating fluidized bed boiler. Example

[0078] Suppose a power plant is equipped with 2×280t / h circulating fluidized bed boilers.

[0079] According to the method proposed in this invention, the following five main processes are calculated respectively:

[0080] Step S01: Collect time-domain data of the process values ​​of high-pressure main pipe steam pressure and coal quantity, steam flow rate, bed temperature, primary air volume, oxygen quantity, and secondary air volume of N=2 boilers, construct the process value time series, and calculate the process value stability based on the time series.

[0081] Step S02: Calculate the optimized baseline and constraint values ​​for the process values ​​of coal feed rate, steam flow rate, bed temperature, primary air volume, oxygen volume, and secondary air volume of N=2 boilers at the current time.

[0082] Step S03: Calculate the upper and lower optimization weights of the oxygen content, primary air volume and secondary air volume ratio of N=2 boilers at the current time.

[0083] Step S04: Construct an optimization controller for oxygen and primary air volume of N=2 boilers, determine the optimization direction using the ton-coal consumption rotation factor, and calculate the optimized setpoints for oxygen and primary air volume that reduce ton-coal consumption at the current time by performing optimization within the constraints.

[0084] Step S05: Compare the setpoints for oxygen and primary air volume of N=2 boilers in the previous period with the time series of process values ​​at the current time to determine whether to optimize the setpoint output. If the condition is met, the optimized setpoint can be added to the current setpoint; otherwise, the optimized setpoint will not be added to the current setpoint.

[0085] Specifically, step S01 includes the following functions:

[0086] Step S101: Obtain the process values ​​of steam pressure of K=3600 high-pressure main pipes, coal quantity of N=2 boilers, steam flow rate, bed temperature, primary air volume, oxygen quantity, and secondary air volume before the current time, and construct them into time series data.

[0087] Given that historical data from the control system is acquired via OPC communication, a time series can be constructed using "Formula 1":

[0088] The time series of steam pressure process values ​​in the high-pressure header is as follows (unit: MPa):

[0089]

[0090] The time series of process values ​​for boiler 1 are as follows:

[0091] Coal quantity (unit: t / h): .

[0092] Steam flow rate (unit: t / h): .

[0093] Bed temperature (unit: °C): .

[0094] Primary air volume (unit: Nm3 / h): .

[0095] Oxygen content (unit: %): .

[0096] Secondary air volume (unit: Nm3 / h): .

[0097] The time series of process values ​​for the second boiler are as follows:

[0098] Coal quantity (unit: t / h): .

[0099] Steam flow rate (unit: t / h): .

[0100] Bed temperature (unit: °C): .

[0101] Primary air volume (unit: Nm3 / h): .

[0102] Oxygen content (unit: %): .

[0103] Secondary air volume (unit: Nm3 / h): .

[0104] Step S102: Calculate the stability of the high-pressure main pipe steam pressure, coal feed rate of N=2 boilers, steam flow rate, bed temperature, primary air volume, oxygen content, and secondary air volume based on the process value mean, oscillation number, and deviation.

[0105] For the process value time series established in "step S101", the mean of the process value, the rate of change and the number of oscillations are calculated first, and then the stability is calculated.

[0106] (1) Calculate the mean of the process values ​​using "Formula 2"

[0107] Average value of steam pressure process in high-pressure header:

[0108] The average process values ​​for boiler 1 are as follows: Coal quantity: .

[0109] Steam flow rate: .

[0110] Bed temperature: .

[0111] Primary air volume: .

[0112] Oxygen content: .

[0113] Secondary air volume: .

[0114] The average process values ​​for the second boiler are as follows:

[0115] Coal quantity: .

[0116] Steam flow rate: .

[0117] Bed temperature: .

[0118] Primary air volume: .

[0119] Oxygen content: .

[0120] Secondary air volume: .

[0121] (2) Calculate the rate of change of process value using "Formula 3"

[0122] Rate of change of steam pressure in high-pressure header: .

[0123] The rate of change of process values ​​for boiler 1 are as follows:

[0124] Coal quantity: .

[0125] Steam flow rate: .

[0126] Bed temperature: .

[0127] Primary air volume: .

[0128] Oxygen content: .

[0129] Secondary air volume: .

[0130] The rate of change of process values ​​for the second boiler are as follows:

[0131] Coal quantity: .

[0132] Steam flow rate: .

[0133] Bed temperature: .

[0134] Primary air volume: .

[0135] Oxygen content: .

[0136] Secondary air volume: .

[0137] (3) Use "Formula 4" to calculate the number of oscillations in the process value.

[0138] Number of oscillations in the high-pressure header steam pressure process value:

[0139] When i=1, , And so on, i=2, ..., 3600, finally obtaining = .

[0140] The number of oscillations in the process values ​​of the first boiler are as follows:

[0141] Coal quantity: When i=1, , And so on, i=2, ..., 3600, finally obtaining = .

[0142] Steam flow rate: When i=1, , And so on, i=2, ..., 3600, finally obtaining = .

[0143] Bed temperature: When i=1, , And so on, i=2, ..., 3600, finally obtaining = .

[0144] Primary air volume: When i=1, , And so on, i=2, ..., 3600, finally obtaining = .

[0145] Oxygen content: When i=1, , And so on, i=2, ..., 3600, finally obtaining = .

[0146] Secondary air volume: When i=1, , And so on, i=2, ..., 3600, finally obtaining = .

[0147] The number of oscillations in the process values ​​of the second boiler are as follows:

[0148] Coal quantity: When i=1, , And so on, i=2, ..., 3600, finally obtaining = .

[0149] Steam flow rate: When i=1, , And so on, i=2, ..., 3600, finally obtaining = .

[0150] Bed temperature: When i=1, , And so on, i=2, ..., 3600, finally obtaining = .

[0151] Primary air volume: When i=1, , And so on, i=2, ..., 3600, finally obtaining = .

[0152] Oxygen content: When i=1, , And so on, i=2, ..., 3600, finally obtaining = .

[0153] Secondary air volume: When i=1, <0.05, And so on, i=2, ..., 3600, finally obtaining = .

[0154] (4) Calculate the deviation of process values ​​using "Formula 5".

[0155] Deviation of steam pressure process value in high-pressure header: .

[0156] The deviations of the process values ​​for boiler 1 are as follows: Coal quantity: .

[0157] Steam flow rate: .

[0158] Bed temperature: .

[0159] Primary air volume: .

[0160] Oxygen content: .

[0161] Secondary air volume: .

[0162] The deviations of the process values ​​for the second boiler are as follows:

[0163] Coal quantity: .

[0164] Steam flow rate: .

[0165] Bed temperature: .

[0166] Primary air volume: .

[0167] Oxygen content: .

[0168] Secondary air volume: .

[0169] (5) Calculate the stability of process values ​​using "Formula 6"

[0170] Stability of steam pressure process value in high-pressure header:

[0171]

[0172] The process stability values ​​for the first boiler are as follows:

[0173] Coal quantity: .

[0174] Steam flow rate: .

[0175] Bed temperature: .

[0176] Primary air volume: .

[0177] Oxygen content: .

[0178] Secondary air volume: .

[0179] The process stability values ​​for the second boiler are as follows:

[0180] Coal quantity: .

[0181] Steam flow rate: .

[0182] Bed temperature: .

[0183] Primary air volume: .

[0184] Oxygen content: .

[0185] Secondary air volume: .

[0186] Specifically, step S02 includes the following functions:

[0187] Step S201: Calculate the optimized baseline values ​​for the process values ​​of coal feed rate, steam flow rate, bed temperature, primary air volume, oxygen volume, and secondary air volume of N=2 boilers at the current time.

[0188] A smoothing transformation is performed on the mean of the process values ​​to avoid the impact of noise and abnormal disturbances on the optimization. The optimization base value is calculated using "Formula 7":

[0189] The optimized baseline values ​​for the process values ​​of the first boiler are as follows:

[0190] Coal quantity: known from the previous cycle The first two cycles The first three cycles , .

[0191] Steam flow rate: known for the previous cycle The first two cycles The first three cycles , .

[0192] Bed temperature: known from the previous cycle The first two cycles The first three cycles , .

[0193] Primary air volume: known from the previous cycle The first two cycles The first three cycles , .

[0194] Oxygen content: known from the previous cycle The first two cycles The first three cycles , .

[0195] Secondary air volume: known from the previous cycle The first two cycles The first three cycles , .

[0196] The optimized baseline values ​​for the process values ​​of the second boiler are as follows:

[0197] Coal quantity: known from the previous cycle The first two cycles The first three cycles , .

[0198] Steam flow rate: known for the previous cycle The first two cycles The first three cycles , .

[0199] Bed temperature: known from the previous cycle The first two cycles The first three cycles , .

[0200] Primary air volume: known from the previous cycle The first two cycles The first three cycles , .

[0201] Oxygen content: known from the previous cycle The first two cycles The first three cycles , .

[0202] Secondary air volume: known from the previous cycle The first two cycles The first three cycles , .

[0203] Step S202: Calculate the optimized constraint values ​​for the process values ​​of coal feed rate, steam flow rate, bed temperature, primary air volume, oxygen volume, and secondary air volume of N=2 boilers at the current time.

[0204] By utilizing the set process value constraints and the stability of the process values, the upper and lower constraints of each process value optimization constraint are calculated comprehensively, serving as the constraints for optimization.

[0205] (1) Calculate the process value using "Formula 8" and "Formula 9" to optimize the upper constraint limit.

[0206] The optimization upper constraints for the first boiler are as follows:

[0207] Coal quantity: Known upper limit of preset coal quantity process value constraint , .

[0208] Steam flow rate: Known preset upper limit of steam flow rate process value constraint , .

[0209] Bed temperature: Known preset bed temperature process value constraint upper limit , .

[0210] Primary air volume: Known preset primary air volume process value constraint upper limit , .

[0211] Oxygen content: Known preset oxygen content process value constraint upper limit , .

[0212] Secondary air volume is known; preset secondary air volume process value constraint upper limit. , .

[0213] Ratio of primary air volume process value to secondary air volume process value: The upper limit of the preset constraint on the ratio of primary air volume process value to secondary air volume process value is known. , .

[0214] The optimized upper constraints for the second boiler are as follows:

[0215] Coal quantity: Known upper limit of preset coal quantity process value constraint , .

[0216] Steam flow rate: Known preset upper limit of steam flow rate process value constraint , .

[0217] Bed temperature: Known preset bed temperature process value constraint upper limit , .

[0218] Primary air volume: Known preset primary air volume process value constraint upper limit , .

[0219] Oxygen content: Known preset oxygen content process value constraint upper limit , .

[0220] Secondary air volume: Known preset secondary air volume process value constraint upper limit , .

[0221] Ratio of primary air volume process value to secondary air volume process value: The upper limit of the preset constraint on the ratio of primary air volume process value to secondary air volume process value is known. , .

[0222] (2) Calculate the process value optimization lower constraint limit using "Formula 10" and "Formula 11".

[0223] The optimization lower constraints for the first boiler are as follows:

[0224] Coal quantity: Lower limit of known preset coal quantity process value constraint , .

[0225] Steam flow rate: Known preset steam flow rate process value constraint lower limit , .

[0226] Bed temperature: Lower limit of known preset bed temperature process value constraint , .

[0227] Primary air volume: Lower limit of the known preset primary air volume process value constraint , .

[0228] Oxygen content: Known preset oxygen content process value constraint lower limit , .

[0229] Secondary air volume: Lower limit of known preset secondary air volume process value constraint , .

[0230] Ratio of primary air volume process value to secondary air volume process value: Given a preset lower limit constraint on the ratio of primary air volume process value to secondary air volume process value. , .

[0231] The optimized lower constraints for the second boiler are as follows:

[0232] Coal quantity: Lower limit of known preset coal quantity process value constraint , .

[0233] Steam flow rate: Known preset steam flow rate process value constraint lower limit , .

[0234] Bed temperature: Lower limit of known preset bed temperature process value constraint , .

[0235] Primary air volume: Lower limit of the known preset primary air volume process value constraint , .

[0236] Oxygen content: Known preset oxygen content process value constraint lower limit , .

[0237] Secondary air volume: Lower limit of known preset secondary air volume process value constraint , .

[0238] Ratio of primary air volume process value to secondary air volume process value: Given a preset lower limit constraint on the ratio of primary air volume process value to secondary air volume process value. , .

[0239] Specifically, step S03 includes the following functions:

[0240] Step S301: Calculate the optimized weight of oxygen based on the process value time series constructed in step S01.

[0241] The time series of oxygen process values ​​is weighted and summed with the optimization constraints to calculate the optimization weight of oxygen. The larger the weight value, the larger the optimization output, and vice versa.

[0242] (1) Calculate the oxygen content optimization weight using "Formula 12" Optimize the oxygen content weight of the first boiler: .

[0243] Optimize the oxygen content weight for the second boiler: .

[0244] (2) Calculate the optimal weight under oxygen content using "Formula 13" Optimization weights for the oxygen content of boiler 1: .

[0245] Optimization weights for oxygen levels in the second boiler: .

[0246] Step S302: Calculate the optimized weight of the ratio of primary air volume to secondary air volume based on the process value time series constructed in step S01.

[0247] The optimization weight of the primary air volume is calculated by weighting and summing the time series ratios of the primary and secondary air volume process values ​​with the optimization constraints. The larger the weight value, the larger the optimization output, and vice versa.

[0248] (1) Calculate the ratio of primary air volume to secondary air volume using "Formula 14" and optimize the weight.

[0249] The weighting of the primary air volume to secondary air volume ratio for the first boiler is optimized as follows: .

[0250] The weighting of the primary air volume to secondary air volume ratio for the second boiler was optimized. .

[0251] (2) Calculate the optimization weights based on the ratio of primary air volume to secondary air volume using Formula 15.

[0252] Optimization weights for the ratio of primary air volume to secondary air volume of boiler 1: .

[0253] Optimization weights for the ratio of primary air volume to secondary air volume of the second boiler: .

[0254] Specifically, step S04 includes the following functions:

[0255] Step S401: Calculate the incremental factor, incremental rotation factor, and number of consecutive flips for the ton of steam coal consumption.

[0256] (1) Calculate the incremental factor of coal consumption per ton of steam using "Formula 16".

[0257] Incremental factor for coal consumption per ton of steam for boiler #1:

[0258] The optimized baseline value of the steam flow rate process in the previous cycle is known. Coal quantity process value optimization base value .

[0259]

[0260] Incremental factor for coal consumption per ton of steam for the second boiler:

[0261] The optimized baseline value of the steam flow rate process in the previous cycle is known. Coal quantity process value optimization base value .

[0262]

[0263] (2) Calculate the rotation factor of ton-to-ton coal consumption increment using "Formula 17".

[0264] Rotation factor for coal consumption per ton of steam for boiler #1:

[0265] Given the rotation factor of ton gasoline consumption in the previous period .

[0266]

[0267] Rotation factor for coal consumption per ton of steam for boiler #2:

[0268] Given the rotation factor of ton gasoline consumption in the previous period .

[0269]

[0270] (3) Calculate the number of consecutive flips using Formula 18

[0271] Continuous rotation counting of coal consumption per ton of steam for boiler #1:

[0272] Given the continuous reversal count of ton gasoline consumption in the previous cycle .

[0273] because ,but

[0274] Continuous rotation counting of coal consumption per ton of steam for boiler #2:

[0275] Given the continuous reversal count of ton gasoline consumption in the previous cycle .

[0276] because ,but .

[0277] Step S402: Construct an oxygen optimization controller and perform oxygen optimization calculations.

[0278] Using "Formula 19" and "Formula 20", an oxygen optimization controller is constructed to perform oxygen optimization calculations.

[0279] Oxygen Optimization Calculation for Boiler No. 1

[0280] Given the number of rotations per ton of gasoline, determine the preset value. oxygen setpoint for the previous cycle Steam flow rate process value optimization base value Coal quantity process value optimization base value .because , , Calculation of oxygen setpoint increment: .

[0281] Oxygen Optimization Calculation for the Second Boiler

[0282] Given the number of rotations per ton of gasoline, determine the preset value. oxygen setpoint for the previous cycle Steam flow rate process value optimization base value Coal quantity process value optimization base value .because , , Calculation of oxygen setpoint increment: .

[0283] Step S403: Construct a primary airflow optimization controller and perform a primary airflow optimization calculation.

[0284] A primary airflow optimization controller is constructed using "Formula 21" and "Formula 22" to perform a primary airflow optimization calculation.

[0285] Primary air volume optimization calculation for the first boiler

[0286] Given the number of rotations per ton of gasoline, determine the preset value. The set value of the primary air volume in the previous cycle Steam flow rate process value optimization base value Coal quantity process value optimization base value .because , , Not valid. , , Calculation of primary air volume setpoint increment: Since the above conditions are not met, therefore, for this cycle... .

[0287] Primary air volume optimization calculation for the second boiler

[0288] Given the number of rotations per ton of gasoline, determine the preset value. The set value of the primary air volume in the previous cycle Steam flow rate process value optimization base value Coal quantity process value optimization base value .because , , , Not valid. , Calculation of primary air volume setpoint increment: Since the above conditions are not met, therefore, for this cycle... .

[0289] Specifically, step S05 includes the following functions:

[0290] Step S501: Compare the current oxygen process value time series with the previous cycle set value. If the condition is met, the optimized set value can be added to the previous cycle oxygen set value and used as the current cycle oxygen set value; otherwise, the set value remains unchanged.

[0291] Use "Formula 23" to calculate the oxygen optimization setpoint update.

[0292] Oxygen setpoint update calculation for boiler #1

[0293] Known oxygen setpoint for the previous cycle When i=1, And so on, the oxygen setpoint is updated when i=2, ..., 3600.

[0294] Oxygen setpoint update calculation for boiler #2

[0295] Known oxygen setpoint for the previous cycle When i=1, And so on, the oxygen setpoint is updated when i=2, ..., 3600.

[0296] Step S503: Compare the time series of the primary air volume process value at the current time with the set value of the previous period. If the condition is met, the optimized set value can be added to the primary air volume set value of the previous period and used as the primary air volume set value of the current period. Otherwise, the set value remains unchanged.

[0297] Use "Formula 24" to perform an update calculation of the air volume optimization setpoint.

[0298] Calculation of primary air volume setpoint update for boiler #1

[0299] Knowing the primary air volume setpoint of the previous cycle When i=1, And so on, the oxygen setpoint is updated when i=2, ..., 3600. .

[0300] Calculation of primary air volume setpoint update for boiler #2

[0301] Knowing the primary air volume setpoint of the previous cycle When i=1, And so on, the oxygen setpoint is updated when i=2, ..., 3600. .

Claims

1. A method for self-optimizing combustion control of a circulating fluidized bed boiler with a main control pipe, characterized in that, Includes the following steps: Step S01: Collect time-domain data of the process values ​​of high-pressure main pipe steam pressure and N boilers coal quantity, steam flow rate, bed temperature, primary air volume, oxygen quantity, and secondary air volume, construct the process value time series, and calculate the process value stability based on the time series. Step S02: Calculate the optimized baseline and constraint values ​​of the process values ​​of coal feed rate, steam flow rate, bed temperature, primary air volume, oxygen volume, and secondary air volume for N boilers at the current time. Step S03: Calculate the upper and lower optimization weights of the oxygen content, primary air volume and secondary air volume ratio of N boilers at the current time; Step S04: Construct oxygen and primary air volume optimization controllers for N boilers, determine the optimization direction using the ton-coal consumption rotation factor, and calculate the optimized setpoints for oxygen and primary air volume that reduce ton-coal consumption at the current time by performing optimization within the constraints. Step S05: Compare the setpoints of oxygen and primary air volume of N boilers in the previous period with the time series of process values ​​at the current time to determine whether to optimize the setpoint output. If the condition is met, the optimized setpoint is added to the current setpoint; otherwise, the optimized setpoint is not added to the current setpoint.

2. The self-optimizing combustion control method for a circulating fluidized bed boiler with a main control pipe according to claim 1, characterized in that, Step S01 includes the following steps: Step S101: Obtain the sampled process values ​​of the K previous high-pressure header steam pressures, N boiler coal quantities, steam flow rates, bed temperatures, primary air volumes, oxygen quantities, and secondary air volumes at the current time, and construct a time series of process values. (Formula 1) in, Represents a time series of process values. This represents the sampled value of the first process value. This represents the sampled value of the Kth process value, where K represents the number of samples. Superscripts include the high-pressure main pipe (main) and the Nth boiler (BN). Subscripts indicate the following parameters: high-pressure steam pressure (STPR), coal quantity (coal), steam flow rate (STFL), bed temperature (BDT), primary air volume (air1), oxygen quantity (O2), and secondary air volume (air2). Step S102: Calculate the stability of high-pressure main pipe steam pressure, coal feed rate of N boilers, steam flow rate, bed temperature, primary air volume, oxygen content, and secondary air volume based on the process value mean, number of oscillations, and deviation. (1) Calculate the mean of the process values ​​according to Formula 1. (Formula 2) in, This represents the average value of the process, m represents the current period, and i is the i-th data point. (2) Calculate the rate of change of process value (Formula 3) in, Indicates the rate of change of process value; (3) Calculate the number of oscillations in the process value like ,but ,otherwise , = (Formula 4) in, This represents an intermediate variable used in calculating the number of oscillations during the process. , Indicates the number of oscillations in the process value; (4) Calculate the deviation of process values (Formula 5) in, Indicates the deviation of process values; (5) Calculation of process value stability (Formula 6) in, Indicates the stability of process values. The preset number of oscillations in the process value is indicated by 'max', which indicates the maximum value to be calculated.

3. The combustion self-optimization control method for a circulating fluidized bed boiler with a main control pipe according to claim 1, characterized in that, Step S02 includes the following steps: Step S201: Calculate the optimized baseline values ​​for the process values ​​of coal feed rate, steam flow rate, bed temperature, primary air volume, oxygen content, and secondary air volume of N boilers at the current time. (Formula 7) in, This represents the base value for process value optimization. This represents the average value of the process, and m represents the current period. Superscripts include the high-pressure main pipe (main) and the Nth boiler (BN). Subscripts indicate the following parameters: high-pressure steam pressure (STPR), coal quantity (coal), steam flow rate (STFL), bed temperature (BDT), primary air volume (air1), oxygen quantity (O2), and secondary air volume (air2). Step S202: Calculate the optimal constraint values ​​for the process values ​​of coal feed rate, steam flow rate, bed temperature, primary air volume, oxygen volume, and secondary air volume of N boilers at the current time. (1) Calculation of upper constraint limits for process value optimization (Formula 8) in, The upper constraint limit for process value optimization, SETHL, represents the upper limit of the set process value constraint. Indicates the stability of process values; This indicates the stability of the high-pressure steam pressure process value of the high-pressure main pipe in the current cycle; The formula for calculating the upper limit of the primary air volume to secondary air volume ratio is as follows: (Formula 9) in, This represents the upper constraint limit for optimizing the ratio of primary air volume to secondary air volume. This indicates the upper limit of the set constraint on the ratio of primary air volume to secondary air volume. This indicates the stability of the primary air volume process value. This indicates the stability of the secondary air volume process value; (2) Calculation of lower constraint limits for process value optimization (Formula 10) in, SETLL represents the lower limit of process value optimization constraints; The lower limit of the primary air volume to secondary air volume ratio is calculated as follows: (Formula 11) in, This represents the lower constraint limit for optimizing the ratio of primary air volume to secondary air volume. This indicates the lower limit of the set constraint on the ratio of primary air volume to secondary air volume.

4. The combustion self-optimization control method for a circulating fluidized bed boiler with a main pipe according to claim 1, characterized in that, Step S03 includes the following steps: Step S301: Calculate the optimized weight of oxygen based on the process value time series constructed in step S01; (1) Calculation of the weight for oxygen content optimization (Formula 12) in, To optimize the weighting of oxygen content in the Nth boiler, This represents the upper constraint limit for oxygen optimization of the Nth boiler, where m represents the current period. Let K represent the oxygen content process value of the Nth boiler, K represent the number of samples, and i represent the i-th data point. (2) Calculation of optimization weights under oxygen content (Formula 13) in, To optimize the weights for the oxygen content of the Nth boiler, This represents the lower constraint limit for oxygen optimization of the Nth boiler; Step S302: Calculate the optimized weight of the ratio of primary air volume to secondary air volume based on the process value time series constructed in step S01; (1) Optimize the weight of the ratio of primary air volume to secondary air volume (Formula 14) in, Optimize the weighting of the ratio of primary air volume to secondary air volume for the Nth boiler. This represents the upper constraint limit for optimizing the ratio of primary air volume to secondary air volume of the Nth boiler. This represents the process value of the primary air volume of the Nth boiler. This represents the process value of the secondary air volume of the Nth boiler; (2) Optimization weight under the ratio of primary air volume to secondary air volume (Formula 15) in, The optimal weights are determined for the ratio of primary air volume to secondary air volume of the Nth boiler. This represents the lower constraint limit for optimizing the ratio of primary air volume to secondary air volume of the Nth boiler. This represents the process value of the primary air volume of the Nth boiler. This represents the process value of the secondary air volume of the Nth boiler.

5. The combustion self-optimization control method for a circulating fluidized bed boiler with a main pipe according to claim 1, characterized in that, Step S04 includes the following steps: Step S401: Calculate the incremental factor, incremental rotation factor, and number of consecutive flips for the ton of gasoline coal consumption; (1) Calculation of incremental factor for coal consumption per ton of steam (Formula 16) in, This represents the incremental factor of coal consumption per ton of steam for the Nth boiler. This represents the baseline value for optimizing the coal quantity process of the Nth boiler. represents the optimized baseline value of the steam flow process value of the Nth boiler, and m represents the current period; (2) Calculation of the rotation factor for the incremental coal consumption per ton of steam (Formula 17) in, This represents the rotation factor for the incremental coal consumption per ton of steam for the Nth boiler; (3) Continuous flip count calculation like , +1 (Formula 18), otherwise , in, This indicates the number of times the ton of steam kerosene is turned over in the Nth boiler; Step S402: Construct an oxygen optimization controller and perform oxygen optimization calculations. when At that time, the oxygen content optimization calculation is as follows: If , At that time, the oxygen optimization setpoint is calculated as follows: (Formula 19), otherwise ; in, This represents the preset value for determining the number of times the ton of steam kerosene is rotated in the Nth boiler. Greater than or equal to Stop optimizing. This represents the oxygen setting value for the Nth boiler. This represents the upper limit of oxygen optimization for the Nth boiler. This represents the optimized oxygen setting value for the Nth boiler. This represents the baseline value for optimizing the coal quantity process of the Nth boiler. This represents the optimized baseline value for the steam flow rate process of the Nth boiler. This indicates the optimization weight for the oxygen content of the Nth boiler; when At that time, the oxygen content optimization calculation is as follows: If , At that time, the oxygen optimization setpoint is calculated as follows: (Formula 20), otherwise Among them, when Greater than or equal to Stop optimizing. This represents the lower limit of oxygen optimization for the Nth boiler. This represents the optimization weight under the oxygen content of the Nth boiler; Step S403: Construct a primary airflow optimization controller and perform primary airflow optimization calculations. when When the primary air volume is optimized, the calculation method is as follows: If , , , , Primary air volume optimization setpoint calculation: (Formula 21), otherwise ; when Greater than or equal to Stop optimizing. This represents the optimized baseline value for the primary air volume process of the Nth boiler. This represents the optimized baseline value for the secondary air volume process of the Nth boiler. This represents the upper limit of the optimization ratio of primary air volume to secondary air volume for the Nth boiler. This represents the optimized baseline value for the bed temperature process of the Nth boiler. This represents the lower limit of the optimal bed temperature for the Nth boiler. This represents the upper limit of the primary air volume optimization for the Nth boiler. This represents the lower limit for optimizing the secondary air volume of the Nth boiler. This represents the optimized setpoint for the primary air volume of the Nth boiler. This represents the setpoint for the secondary air volume of the Nth boiler; when When the primary air volume is optimized, the calculation method is as follows: If , , , , Primary air volume optimization setpoint calculation: (Formula 22), otherwise ; in, This represents the preset value for determining the number of times the ton of steam kerosene is rotated in the Nth boiler. Greater than or equal to Stop optimizing. This represents the lower limit of the optimization ratio of primary air volume to secondary air volume for the Nth boiler. This represents the upper limit of the optimized bed temperature for the Nth boiler. This represents the lower limit of the primary air volume optimization for the Nth boiler. This represents the upper limit of the optimized secondary air volume for the Nth boiler. The optimal weights are determined by the ratio of primary air volume to secondary air volume for the Nth boiler.

6. The combustion self-optimization control method for a circulating fluidized bed boiler with a main pipe according to claim 1, characterized in that, Step S05 includes the following steps: Step S501: The time series of oxygen process value at the current time is compared with the set value of the previous period. If the condition is met, the optimized set value is added to the oxygen set value of the previous period and used as the oxygen set value of the current period. Otherwise, the set value remains unchanged. The time series of oxygen quantity process values ​​is analyzed, and when each time series value is... When (i=1, ..., K), the oxygen optimization setting value is accumulated to the oxygen setting value of the previous cycle and used as the oxygen setting value of the current cycle; otherwise, the oxygen setting value of the current cycle remains unchanged. (Formula 23) in, This represents the time series of oxygen quantity values ​​for the Nth boiler. This represents the oxygen setting value for the Nth boiler. This represents the optimized oxygen setting value for the Nth boiler, where m represents the current cycle and K represents the number of samples. Step S502: The time series of the primary air volume process value at the current time is compared with the set value of the previous period. If the condition is met, the optimized set value can be added to the primary air volume set value of the previous period and used as the primary air volume set value of the current period. Otherwise, the set value remains unchanged. The time series of primary air volume values ​​is analyzed, and when each value in the time series... When i=1, ..., K, the primary air volume optimization setting value is accumulated to the primary air volume setting value of the previous cycle and used as the primary air volume setting value of the current cycle; otherwise, the primary air volume setting value of the current cycle remains unchanged. (Official 24) in, This represents the time series of primary air volume values ​​for the Nth boiler. This represents the primary air volume setpoint for the Nth boiler. This represents the optimized primary air volume setting for the Nth boiler.

7. A self-optimizing combustion control system for a circulating fluidized bed boiler with a main control pipe, characterized in that: Includes a data acquisition module, used to collect process data on main pipe pressure and coal quantity, steam flow, bed temperature, primary air volume, oxygen quantity, and secondary air volume of N boilers; The control module is used to implement the self-optimization control method for combustion of a circulating fluidized bed boiler with a main pipe as described in any one of claims 1-6, and to write back the calculation results to the field control system. The data storage module is used to store the data generated by the data acquisition module and the control module; The data is displayed by the data storage module.

8. A computing device comprising a memory, a processor, and a program stored in the memory and running on the processor, characterized in that, When the processor executes the program, it implements the steps of the self-optimizing combustion control method for a circulating fluidized bed boiler as described in any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores instructions that, when executed on a terminal device, cause the terminal device to perform the steps of the self-optimizing combustion control method for a circulating fluidized bed boiler as described in any one of claims 1 to 6.

Images