CN111623369B - Control method for adjusting boiler fuel feeding quantity by using smoke oxygen content signal - Google Patents

Control method for adjusting boiler fuel feeding quantity by using smoke oxygen content signal Download PDF

Info

Publication number
CN111623369B
CN111623369B CN202010594628.4A CN202010594628A CN111623369B CN 111623369 B CN111623369 B CN 111623369B CN 202010594628 A CN202010594628 A CN 202010594628A CN 111623369 B CN111623369 B CN 111623369B
Authority
CN
China
Prior art keywords
interval
load
oxygen content
formula
control
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN202010594628.4A
Other languages
Chinese (zh)
Other versions
CN111623369A (en
Inventor
沈赫男
王承亮
牛树森
周怀春
孙铭泽
周远科
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huadian International Power Co ltd Technical Service Branch
Wuhan Jiuzhou 3d Combustion Technology Co ltd
Northeast Electric Power University
Huadian Weifang Power Generation Co Ltd
Original Assignee
Huadian International Power Co ltd Technical Service Branch
Wuhan Jiuzhou 3d Combustion Technology Co ltd
Northeast Dianli University
Huadian Weifang Power Generation Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huadian International Power Co ltd Technical Service Branch, Wuhan Jiuzhou 3d Combustion Technology Co ltd, Northeast Dianli University, Huadian Weifang Power Generation Co Ltd filed Critical Huadian International Power Co ltd Technical Service Branch
Priority to CN202010594628.4A priority Critical patent/CN111623369B/en
Publication of CN111623369A publication Critical patent/CN111623369A/en
Application granted granted Critical
Publication of CN111623369B publication Critical patent/CN111623369B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F23COMBUSTION APPARATUS; COMBUSTION PROCESSES
    • F23NREGULATING OR CONTROLLING COMBUSTION
    • F23N1/00Regulating fuel supply
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F23COMBUSTION APPARATUS; COMBUSTION PROCESSES
    • F23NREGULATING OR CONTROLLING COMBUSTION
    • F23N2239/00Fuels
    • F23N2239/02Solid fuels
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F23COMBUSTION APPARATUS; COMBUSTION PROCESSES
    • F23NREGULATING OR CONTROLLING COMBUSTION
    • F23N2241/00Applications
    • F23N2241/10Generating vapour

Landscapes

  • Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Control Of Steam Boilers And Waste-Gas Boilers (AREA)
  • Regulation And Control Of Combustion (AREA)

Abstract

The invention relates to a control method for fuel supply, in particular to a control method for adjusting the amount of fuel fed into a boiler by using a smoke oxygen content signal. The method comprises the following steps: establishing a control rule which is superior to an optimization target according to the load and air temperature distribution, acquiring the optimal smoke oxygen content from the optimization control rule, and performing deviation operation according to the current smoke oxygen content signal and the optimal smoke oxygen content to obtain a total fuel correction signal. And accurately controlling the amount of fuel entering the furnace according to the total fuel correction signal. And the control operation adopts an incremental PID algorithm. According to the invention, the fuel quantity is subjected to advanced and accurate intervention by the optimal change rule of the oxygen content of the flue gas along with the load and the steam temperature, so that the response speed of the load of the coal-fired unit is increased; the safety and the stability of the boiler in the operation process are improved, the heat efficiency of a power plant is effectively improved, and the energy consumption is reduced.

Description

Control method for adjusting boiler fuel feeding quantity by using smoke oxygen content signal
Technical Field
The invention relates to a control method for fuel supply, in particular to a control method for adjusting the amount of fuel fed into a boiler by using a smoke oxygen content signal.
Background
The effective improvement of the operating thermal efficiency of a power plant and the reduction of the energy consumption are inevitable requirements for the development of the economy and the power industry of China, wherein whether the control and adjustment methods of the air quantity and the fuel quantity are effective or not is an important factor influencing the combustion thermal efficiency of the boiler. For a coal-fired power plant, the main tasks of a combustion control system of the coal-fired power plant are to control the fuel quantity in the operation process to meet the requirements of steam pressure or load output, adjust the air supply quantity through the oxygen content deviation of flue gas and simultaneously ensure the safety, stability and economy of the operation process of a boiler.
At present, the control of air supply quantity is mostly based on the measurement of oxygen content of flue gas, the numerical value of the oxygen content of operating flue gas is firstly detected, and then the control of the air supply quantity and the control of the oxygen content of operating flue gas are completed by an air supply quantity control device. The control mode is based on a control strategy of 'wind-coal ratio' and 'water-coal ratio', and realizes the regulation of fuel entering the furnace by adopting a regulation means of relevant parameter deviation-feedback-output, and has the disadvantages of hysteresis, simple model, small amount of detection signals, single control loop, no analysis and utilization of 'mass data' in each operation condition of the unit, incapability of performing overall optimization control on fuel quantity, air door opening degree, water supply quantity and the like according to the change conditions of coal quality and load in time, incapability of adapting to the quick load change condition, incapability of adapting to the regulation performance change caused by the characteristic change of external equipment such as a fan and the like, and poor disturbance resistance of the system.
Disclosure of Invention
Aiming at the defects, the invention discloses a control method for adjusting the amount of fuel fed into a boiler by using a smoke oxygen content signal, which can accurately control the amount of fuel according to the change conditions of load and steam temperature.
The invention relates to a control method for adjusting the amount of fuel fed into a boiler by using a smoke oxygen content signal, which comprises the following steps of:
1) establishing a control rule which is superior to an optimization target according to the load and air temperature distribution:
acquiring historical operating data of a coal-fired unit:
1a, determining a parameter x extracted from a DCS historical databasekThe parameter types include: actual load, total fuel quantity, fuel quantity of each layer of coal feeder, primary air quantity of each layer of coal mill, total secondary air quantity, secondary air door opening degree of each layer, rotating speed of each layer of coal feeder, total air quantity, water supply flow, main steam pressure, regulation level pressure, flue gas oxygen content, air temperature of air feeder inlet, NOXConcentration and superheater wall temperature, wherein the parameter classes are divided by K, K being 1,2, …, K;
1b, setting a data time interval value;
and 1c, extracting the DCS historical data of the coal-fired unit in the past year by taking the data unit of the parameter at the same moment as a data packet according to the parameter listed in the step 1a and the data time interval value set in the step 1 b.
Dividing two-dimensional intervals according to the load and the temperature, and classifying the historical data according to the two-dimensional intervals:
2a, determining a load optimization interval according to a main annual load operation section: the highest operation load of the coal-fired unit in the main operation section is SmaxThe lowest running load is SminThen the load optimization interval is [ S ]min,Smax];
2b, determining the number of load intervals: setting the load division interval to LSThe divided load interval is represented by i, i is 1,2, …, m, and the number of load intervals m is obtained by the formula (1):
m=(Smax-Smin)/LS (1);
2c, determining an air temperature optimization interval according to the annual air temperature change of an air feeder inlet: the maximum air temperature at the inlet of the annual air feeder is set as TmaxThe lowest wind temperature is TminIf the temperature is within the optimum range [ T ]min,Tmax];
2d, determining the number of temperature intervals: air temperature divisionInterval LTAnd j represents a divided air temperature section, and j is 1,2, …, n, the number of air temperature sections n is obtained by the formula (2):
n=(Tmax-Tmin)/LT (2);
2e, the two-dimensional interval of the i-th load interval and the j-th air temperature interval obtained in the steps 2a to 2d is formula (3)
{[Smin+(i-1)×LS,Smin+i×LS],[Tmin+(j-1)×LT,Tmin+j×LT]} (3)。
Calculating the average value of various parameters, the average value of boiler efficiency and NO of each two-dimensional intervalXEmission mean value:
3a, classifying the data packets of the historical data of the coal-fired unit obtained in the first step according to the two-dimensional interval divided in the second step, and discarding the data packets exceeding the two-dimensional interval of the load and the air temperature;
3b, counting the number L of the data packets in each two-dimensional interval in the step 3a, wherein L is 1,2, …, and Li,j(ii) a The history data in the two-dimensional interval expressed by equation (3) is expressed as: x is the number ofi,j,k,lThe mean value of various parameters of the two-dimensional interval is shown in formula (4):
Figure BDA0002557003560000021
in the formula:
Figure BDA0002557003560000022
means for expressing the mean value of the kth parameter in the ith load interval and the jth air temperature interval;
and 3c, defining the average value of the boiler efficiency as the average value of the actual load/the total fuel quantity, wherein the average value of the boiler efficiency is shown in formula (5):
Figure BDA0002557003560000023
in the formula:
Figure BDA0002557003560000024
and
Figure BDA0002557003560000025
is the average value of the total fuel quantity of the two-dimensional interval,
Figure BDA0002557003560000026
and
Figure BDA0002557003560000027
is the actual load mean value of the two-dimensional interval, wherein
Figure BDA0002557003560000028
3d, calculating NOXMean value of emissions
Figure BDA0002557003560000029
And is
Figure BDA00025570035600000210
Fourthly, screening the data packets in the two-dimensional interval according to the optimization target:
in the two-dimensional interval formed by the ith load interval and the jth gas temperature interval, the average value of the efficiency is superior to that of the boiler
Figure BDA00025570035600000211
And NOXMean value of emissions
Figure BDA00025570035600000212
To optimize the data packet, the screening condition is shown in formula (6):
Figure BDA00025570035600000213
in the formula: x is the number ofi,j,E,lAnd
Figure BDA00025570035600000216
cooker for indicating ith load interval, jth air temperature interval and ith data packetFurnace efficiency and NOXIs discharged, wherein
Figure BDA00025570035600000214
Each two-dimensional interval is superior to the average efficiency and NO of the boilerXProcessing of data packets of emission mean:
screening the data packets meeting the optimization target, storing the data packets in the original two-dimensional interval, and counting the number L 'of the data packets meeting the optimization target, wherein L' is 1,2, … and Li,jThe' packets that do not meet the optimization objective are discarded.
Calculating the average value of various parameters in each two-dimensional interval from the data set superior to the optimization target:
and (3) calculating the average value of various parameters of the two-dimensional interval again for the data packets meeting the optimization target as shown in formula (7):
Figure BDA00025570035600000215
in the formula:
Figure BDA0002557003560000031
representing the mean value of the screened ith load interval, the screened jth air temperature interval and the screened kth parameter; x'i,j,k,l′And (4) showing the ith load interval, the jth air temperature interval and the ith data of the kth parameter after screening, thereby obtaining an optimized control rule superior to an optimized target according to the load and air temperature distribution.
2) Obtaining the optimal oxygen content of the smoke from an optimization control law:
the optimal control rule comprises the optimal values of all analysis parameters, including the optimal oxygen content of the flue gas under different loads and temperatures
Figure BDA0002557003560000032
And is
Figure BDA0002557003560000033
3) And (3) performing deviation calculation according to the current oxygen content signal of the flue gas and the optimal oxygen content of the flue gas to obtain a total fuel correction signal:
i, measuring the current oxygen content signal O of the smoke by a smoke oxygen content measuring device2(t);
II, using the current oxygen content signal O of the smoke2(t) and optimum flue gas oxygen content value
Figure BDA0002557003560000034
Obtaining an oxygen deviation signal e (t) by deviation calculation:
Figure BDA0002557003560000035
III, performing control operation on the oxygen amount deviation signal e (t) to obtain a fuel amount control correction signal delta w (t);
IV, when a total fuel correction signal delta w (t) with a positive value is output after the oxygen amount deviation signal e (t) is input, the input of the fuel amount is increased; when the total fuel correction signal delta w (t) with a negative value is output after the oxygen amount deviation signal e (t) is input, the input of the fuel amount is reduced.
Preferably, the control operation in step iii is performed by using an incremental PID algorithm, and the PID algorithm used in step iii needs to be discretized, where the PID algorithm formula is as follows:
Figure BDA0002557003560000036
wherein u is a control signal, wherein KPIs a proportionality coefficient, T1To integrate the time constant, TDIs a differential time constant, u0Are control constants. The specific control operation steps are as follows:
taking T as a sampling period and s as a sampling number (s is 1,2, …, q), corresponding discrete sampling time qT and T, replacing integration by summation and differentiation by increment, and performing approximate transformation:
Figure BDA0002557003560000037
a discrete PID expression can be obtained:
Figure BDA0002557003560000038
in the formula: u. ofsControl number of s-th sampling time, esDeviation value input for s-th sampling time, es-1Is the deviation of the input at the s-1 th sampling instant.
The control signal output by the controller at the s-1 th sampling moment is as follows:
Figure BDA0002557003560000041
and (3) performing difference operation on the formula (11) and the formula (12) to obtain an incremental PID control formula:
Figure BDA0002557003560000042
in the formula
Figure BDA0002557003560000043
Determining the values of A, B and C by using a constant sampling period T, and calculating the control increment delta u by using the deviation of three measurement valuessThe control increment Δ usFurther calculation is performed to generate a total fuel correction signal Δ w (t).
The invention has the beneficial effects that: according to the invention, the fuel quantity is subjected to advanced and accurate intervention by the optimal change rule of the oxygen content of the flue gas along with the load and the steam temperature, so that the response speed of the load of the coal-fired unit is increased; the safety and the stability of the boiler in the operation process are improved, the heat efficiency of a power plant is effectively improved, and the energy consumption is reduced.
Drawings
The invention will be further described with reference to the accompanying drawings, which are simplified schematic drawings and illustrate only the basic structure of the invention in a schematic manner and do not constitute a limitation on the invention:
FIG. 1 is a general flow chart of a control method for adjusting the amount of fuel fed into a boiler by using a flue gas oxygen content signal according to the present invention;
FIG. 2 is a diagram of an incremental PID control system.
Detailed Description
It should be noted that the following detailed description is intended to provide further explanation of the disclosure, and the terminology used is for the purpose of describing the particular embodiments only, and is not intended to be limiting of exemplary embodiments based on the disclosure.
Referring to fig. 1 and 2, the control method for adjusting the amount of fuel fed into a boiler by using a flue gas oxygen content signal according to the present invention includes the following steps:
1) establishing a control rule which is superior to an optimization target according to the load and air temperature distribution:
acquiring historical operating data of a coal-fired unit:
1a, determining a parameter x extracted from a DCS historical databasekThe parameter types include: actual load, total fuel quantity, fuel quantity of each layer of coal feeder, primary air quantity of each layer of coal mill, total secondary air quantity, secondary air door opening degree of each layer, rotating speed of each layer of coal feeder, total air quantity, water supply flow, main steam pressure, regulation level pressure, flue gas oxygen content, air temperature of air feeder inlet, NOXConcentration and superheater wall temperature, wherein the parameter classes are divided by K, K being 1,2, …, K;
1b, setting a data time interval value;
and 1c, extracting the DCS historical data of the coal-fired unit in the past year by taking the data unit of the parameter at the same moment as a data packet according to the parameter listed in the step 1a and the data time interval value set in the step 1 b.
Dividing two-dimensional intervals according to the load and the temperature, and classifying the historical data according to the two-dimensional intervals:
2a, determining a load optimization interval according to a main annual load operation section: setting the maximum operating load of the coal-fired unit in the main operating sectionIs SmaxThe lowest running load is SminThen the load optimization interval is [ S ]min,Smax];
2b, determining the number of load intervals: setting the load division interval to LSI represents the divided load section, i is 1,2, …, m, and the number of load sections m is obtained by the formula (1),
m=(Smax-Smin)/LS (1);
2c, determining an air temperature optimization interval according to the annual air temperature change of an air feeder inlet: the maximum air temperature at the inlet of the annual air feeder is set as TmaxThe lowest wind temperature is TminIf the temperature is within the optimum range [ T ]min,Tmax];
2d, determining the number of temperature intervals: setting the temperature division interval to LTAnd j represents a divided air temperature section, and j is 1,2, …, n, the number of air temperature sections n is obtained by the formula (2):
n=(Tmax-Tmin)/LT (2);
2e, the two-dimensional interval of the i-th load interval and the j-th air temperature interval obtained in the steps 2a to 2d is formula (3)
{[Smin+(i-1)×LS,Smin+i×LS],[Tmin+(j-1)×LT,Tmin+j×LT]} (3)。
Calculating the average value of various parameters, the average value of boiler efficiency and NO of each two-dimensional intervalXEmission mean value:
3a, classifying the data packets of the historical data of the coal-fired unit obtained in the first step according to the two-dimensional interval divided in the second step, and discarding the data packets exceeding the two-dimensional interval of the load and the air temperature;
3b, counting the number L of the data packets in each two-dimensional interval in the step 3a, wherein L is 1,2, …, and Li,j(ii) a The history data in the two-dimensional interval expressed by equation (3) is expressed as: x is the number ofi,j,k,lThe mean value of various parameters of the two-dimensional interval is shown in formula (4):
Figure BDA0002557003560000051
in the formula:
Figure BDA0002557003560000052
means for expressing the mean value of the kth parameter in the ith load interval and the jth air temperature interval;
and 3c, defining the average value of the boiler efficiency as the average value of the actual load/the total fuel quantity, wherein the average value of the boiler efficiency is shown in formula (5):
Figure BDA0002557003560000053
in the formula:
Figure BDA0002557003560000054
and
Figure BDA0002557003560000055
is the average value of the total fuel quantity of the two-dimensional interval,
Figure BDA0002557003560000056
and
Figure BDA0002557003560000057
is the actual load mean value of the two-dimensional interval, wherein
Figure BDA0002557003560000058
3d, calculating NOXMean value of emissions
Figure BDA0002557003560000059
And is
Figure BDA00025570035600000510
Fourthly, screening the data packets in the two-dimensional interval according to the optimization target:
in the two-dimensional interval formed by the ith load interval and the jth gas temperature interval, the average value of the efficiency is superior to that of the boiler
Figure BDA00025570035600000511
And NOXMean value of emissions
Figure BDA00025570035600000512
To optimize the data packet, the screening condition is shown in formula (6):
Figure BDA00025570035600000513
in the formula: x is the number ofi,j,E,lAnd
Figure BDA00025570035600000515
indicating the boiler efficiency and NO of the ith load interval, the jth air temperature interval and the ith data packetXIs discharged, wherein
Figure BDA00025570035600000514
Each two-dimensional interval is superior to the average efficiency and NO of the boilerXProcessing of data packets of emission mean:
screening the data packets meeting the optimization target, storing the data packets in the original two-dimensional interval, and counting the number L 'of the data packets meeting the optimization target, wherein L' is 1,2, … and Li,jThe' packets that do not meet the optimization objective are discarded.
Calculating the average value of various parameters in each two-dimensional interval from the data set superior to the optimization target:
and (3) calculating the average value of various parameters of the two-dimensional interval again for the data packets meeting the optimization target as shown in formula (7):
Figure BDA0002557003560000061
in the formula:
Figure BDA0002557003560000062
representing the mean value of the screened ith load interval, the screened jth air temperature interval and the screened kth parameter; x'i,j,k,l′And (4) showing the ith load interval, the jth air temperature interval and the ith data of the kth parameter after screening, thereby obtaining an optimized control rule superior to an optimized target according to the load and air temperature distribution.
2) Obtaining the optimal oxygen content of the smoke from an optimization control law:
the optimal control rule comprises the optimal values of all analysis parameters, including the optimal oxygen content of the flue gas under different loads and temperatures
Figure BDA0002557003560000063
And is
Figure BDA0002557003560000064
3) And (3) performing deviation calculation according to the current oxygen content signal of the flue gas and the optimal oxygen content of the flue gas to obtain a total fuel correction signal:
i, measuring the current oxygen content signal O of the smoke by a smoke oxygen content measuring device2(t);
II, using the current oxygen content signal O of the smoke2(t) and optimum flue gas oxygen content value
Figure BDA0002557003560000071
Obtaining an oxygen deviation signal e (t) by deviation calculation:
Figure BDA0002557003560000072
III, performing control calculation on the oxygen amount deviation signal e (t) to obtain a fuel amount control correction signal delta w (t):
the control operation adopts an incremental PID algorithm, the used PID algorithm needs discretization treatment, and the formula of the PID algorithm is as follows:
Figure BDA0002557003560000073
wherein u is a control signal, wherein KPIs a proportionality coefficient, T1Is an integration time constant,TDIs a differential time constant, u0Are control constants.
The specific control operation steps are as follows:
taking T as a sampling period and s as a sampling number (s is 1,2, …, q), corresponding discrete sampling time qT and T, replacing integration by summation and differentiation by increment, and performing approximate transformation:
Figure BDA0002557003560000074
a discrete PID expression can be obtained:
Figure BDA0002557003560000075
in the formula: u. ofsControl number of s-th sampling time, esDeviation value input for s-th sampling time, es-1Is the deviation of the input at the s-1 th sampling instant.
The control signal output by the controller at the s-1 th sampling moment is as follows:
Figure BDA0002557003560000076
and (3) performing difference operation on the formula (11) and the formula (12) to obtain an incremental PID control formula:
Figure BDA0002557003560000077
Figure BDA0002557003560000081
in the formula
Figure BDA0002557003560000082
By constant miningDetermining values of A, B and C, and calculating control increment delta u by using the deviation of three measurement values before and after sampling period TsThe control increment Δ usFurther operation is performed (here, operation is Δ u)sAnd a transfer function matched with the boiler combustion) to generate a total fuel correction signal delta w (t).
IV, when a total fuel correction signal delta w (t) with a positive value is output after the oxygen amount deviation signal e (t) is input, the input of the fuel amount is increased; when the total fuel correction signal delta w (t) with a negative value is output after the oxygen amount deviation signal e (t) is input, the input of the fuel amount is reduced.
According to the invention, the fuel quantity is subjected to advanced and accurate intervention by the optimal change rule of the oxygen content of the flue gas along with the load and the steam temperature, so that the response speed of the load of the coal-fired unit is increased; the safety and the stability of the boiler in the operation process are improved, the heat efficiency of a power plant is effectively improved, and the energy consumption is reduced.
It is worth mentioning that the fuel quantity control of the present invention is performed under the condition that the air supply system is under the optimal control.
In light of the foregoing description of the preferred embodiment of the present invention, it is to be understood that various changes and modifications may be made by one skilled in the art without departing from the spirit and scope of the invention. The technical scope of the present invention is not limited to the content of the specification, and must be determined according to the scope of the claims.

Claims (2)

1. A control method for adjusting the amount of fuel fed into a boiler by using a smoke oxygen content signal is characterized by comprising the following steps of: the method comprises the following steps:
1) establishing a control rule which is superior to an optimization target according to the load and air temperature distribution:
acquiring historical operating data of a coal-fired unit:
1a, determining a parameter x extracted from a DCS historical databasekThe parameter types include: actual load, total fuel quantity, fuel quantity of each layer of coal feeders, primary air quantity of each layer of coal mills, total secondary air quantity, opening degree of each layer of secondary air doors and rotating speed of each layer of coal feedersTotal air quantity, feed water flow, main steam pressure, regulation stage pressure, flue gas oxygen content, air temperature of air feeder inlet, NOXConcentration and superheater wall temperature, where the parameter classes are divided by K, K is 1,2, …, K,
1b, setting a data time interval value,
1c, extracting the DCS historical data of the coal-fired unit in the past year by taking the data unit of the parameter at the same moment as a data packet according to the parameter listed in the step 1a and the data time interval value set in the step 1b,
dividing two-dimensional intervals according to the load and the temperature, and classifying the historical data according to the two-dimensional intervals:
2a, determining a load optimization interval according to a main annual load operation section: the highest operation load of the coal-fired unit in the main operation section is SmaxThe lowest running load is SminThen the load optimization interval is [ S ]min,Smax],
2b, determining the number of load intervals: setting the load division interval to LSI represents the divided load section, i is 1,2, …, m, and the number of load sections m is obtained by the formula (1),
m=(Smax-Smin)/LS (1)
2c, determining an air temperature optimization interval according to the annual air temperature change of an air feeder inlet: the maximum air temperature at the inlet of the annual air feeder is set as TmaxThe lowest wind temperature is TminIf the temperature is within the optimum range [ T ]min,Tmax],
2d, determining the number of temperature intervals: setting the temperature division interval to LTAnd j represents a divided air temperature section, and j is 1,2, …, n, the number of air temperature sections n is obtained by the formula (2):
n=(Tmax-Tmin)/LT (2)
2e, the two-dimensional interval of the i-th load interval and the j-th air temperature interval obtained in the steps 2a to 2d is formula (3)
{[Smin+(i-1)×LS,Smin+i×LS],[Tmin+(j-1)×LT,Tmin+j×LT]} (3)
Calculating the average value of various parameters, the average value of boiler efficiency and NO of each two-dimensional intervalXEmission mean value:
3a, classifying the data packets of the historical data of the coal-fired unit obtained in the step one according to the two-dimensional interval divided in the step two, discarding the data packets exceeding the two-dimensional interval of the load and the air temperature,
3b, counting the number L of the data packets in each two-dimensional interval in the step 3a, wherein L is 1,2, …, and Li,j(ii) a The history data in the two-dimensional interval expressed by equation (3) is expressed as: x is the number ofi,j,k,lThe mean value of various parameters of the two-dimensional interval is shown in formula (4):
Figure FDA0002557003550000011
in the formula:
Figure FDA0002557003550000012
represents the mean value of the kth parameter in the ith load interval and the jth air temperature interval,
and 3c, defining the average value of the boiler efficiency as the average value of the actual load/the total fuel quantity, wherein the average value of the boiler efficiency is shown in formula (5):
Figure FDA0002557003550000021
in the formula:
Figure FDA0002557003550000022
and
Figure FDA0002557003550000023
is the average value of the total fuel quantity of the two-dimensional interval,
Figure FDA0002557003550000024
and
Figure FDA0002557003550000025
is the actual load mean value of the two-dimensional interval, wherein
Figure FDA0002557003550000026
3d, calculating NOXMean value of emissions
Figure FDA0002557003550000027
And is
Figure FDA0002557003550000028
Fourthly, screening the data packets in the two-dimensional interval according to the optimization target:
in the two-dimensional interval formed by the ith load interval and the jth gas temperature interval, the average value of the efficiency is superior to that of the boiler
Figure FDA0002557003550000029
And NOXMean value of emissions
Figure FDA00025570035500000210
To optimize the data packet, the screening condition is shown in formula (6):
Figure FDA00025570035500000211
in the formula: x is the number ofi,j,E,lAnd
Figure FDA00025570035500000219
indicating the boiler efficiency and NO of the ith load interval, the jth air temperature interval and the ith data packetXIs discharged, wherein
Figure FDA00025570035500000212
All two dimensionsInterval is superior to average boiler efficiency and NOXProcessing of data packets of emission mean:
the data packets conforming to the optimization target are screened and stored in the original two-dimensional interval, and the number l' of the data packets conforming to the optimization target is counted, wherein
l′=1,2,…,Li,j' packets that do not meet the optimization objective are eliminated,
calculating the average value of various parameters in each two-dimensional interval from the data set superior to the optimization target:
and (3) calculating the average value of various parameters of the two-dimensional interval again for the data packets meeting the optimization target as shown in formula (7):
Figure FDA00025570035500000213
in the formula:
Figure FDA00025570035500000214
representing the mean value of the screened ith load interval, the screened jth air temperature interval and the screened kth parameter; x'i,j,k,l′The ith data of the screened ith load interval, the jth air temperature interval and the kth parameter are expressed, so that the optimized control rule superior to the optimized target according to the load and air temperature distribution is obtained,
2) obtaining the optimal oxygen content of the smoke from an optimization control law:
the optimal control rule comprises the optimal values of all analysis parameters, including the optimal oxygen content of the flue gas under different loads and temperatures
Figure FDA00025570035500000215
And is
Figure FDA00025570035500000216
3) And (3) performing deviation calculation according to the current oxygen content signal of the flue gas and the optimal oxygen content of the flue gas to obtain a total fuel correction signal:
i, measuring the current oxygen content signal O of the smoke by a smoke oxygen content measuring device2(t),
II, using the current oxygen content signal O of the smoke2(t) and optimum flue gas oxygen content value
Figure FDA00025570035500000217
Obtaining an oxygen deviation signal e (t) by deviation calculation:
Figure FDA00025570035500000218
III, performing control operation on the oxygen amount deviation signal e (t) to obtain a fuel amount control correction signal delta w (t),
and IV, when a total fuel correction signal delta w (t) with a positive value is output after the oxygen amount deviation signal e (t) is input, the input of the fuel amount is increased, and when the total fuel correction signal delta w (t) with a negative value is output after the oxygen amount deviation signal e (t) is input, the input of the fuel amount is reduced.
2. The control method for regulating the amount of fuel fed into a boiler by using a flue gas oxygen content signal according to claim 1, wherein the control method comprises the following steps: and step III, selecting an incremental PID algorithm for the control operation, wherein the PID algorithm needs to be discretized, and the PID algorithm formula is as follows:
Figure FDA0002557003550000031
wherein u is a control signal, wherein KPIs a proportionality coefficient, T1To integrate the time constant, TDIs a differential time constant, u0In order to control the constant amount of the liquid,
the specific control operation steps are as follows:
taking T as a sampling period and s as a sampling number (s is 1,2, …, q), corresponding discrete sampling time qT and T, replacing integration by summation and differentiation by increment, and performing approximate transformation:
Figure FDA0002557003550000032
a discrete PID expression can be obtained:
Figure FDA0002557003550000033
in the formula: u. ofsControl number of s-th sampling time, esDeviation value input for s-th sampling time, es-1For the deviations input at the s-1 th sampling instant,
the control signal output by the controller at the s-1 th sampling moment is as follows:
Figure FDA0002557003550000034
and (3) performing difference operation on the formula (11) and the formula (12) to obtain an incremental PID control formula:
Figure FDA0002557003550000035
in the formula
Figure FDA0002557003550000036
Determining the values of A, B and C by using a constant sampling period T, and calculating the control increment delta u by using the deviation of three measurement valuessThe control increment Δ usFurther calculation is performed to generate a total fuel correction signal Δ w (t).
CN202010594628.4A 2020-06-28 2020-06-28 Control method for adjusting boiler fuel feeding quantity by using smoke oxygen content signal Expired - Fee Related CN111623369B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010594628.4A CN111623369B (en) 2020-06-28 2020-06-28 Control method for adjusting boiler fuel feeding quantity by using smoke oxygen content signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010594628.4A CN111623369B (en) 2020-06-28 2020-06-28 Control method for adjusting boiler fuel feeding quantity by using smoke oxygen content signal

Publications (2)

Publication Number Publication Date
CN111623369A CN111623369A (en) 2020-09-04
CN111623369B true CN111623369B (en) 2022-01-25

Family

ID=72257614

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010594628.4A Expired - Fee Related CN111623369B (en) 2020-06-28 2020-06-28 Control method for adjusting boiler fuel feeding quantity by using smoke oxygen content signal

Country Status (1)

Country Link
CN (1) CN111623369B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112856478A (en) * 2021-01-14 2021-05-28 新智数字科技有限公司 Method, device, equipment and medium for adjusting air-fuel ratio of gas boiler
CN113236387B (en) * 2021-04-23 2022-11-29 国网浙江省电力有限公司电力科学研究院 Gas turbine coal machine combined power generation system and method
CN113834092B (en) * 2021-10-12 2023-02-14 中国矿业大学 Combustion optimization control method for gas-fired boiler
CN113834091B (en) * 2021-10-12 2022-11-11 中国矿业大学 Control method for combustion optimization air supply system of gas boiler
CN116481033B (en) * 2023-05-19 2024-03-05 北京中科润宇环保科技股份有限公司 Automatic combustion control method for universal household garbage incinerator

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102183015A (en) * 2011-04-01 2011-09-14 北京三博中自科技有限公司 Combustion optimizing control system of circulating fluidized bed boiler in wide-range change of load
CN102425807A (en) * 2011-11-23 2012-04-25 华北电力大学(保定) Combustion feedforward and feedback composite optimization controlling method for pulverized coal fired boiler
CN103576655A (en) * 2013-11-06 2014-02-12 华北电力大学(保定) Method and system for utility boiler combustion subspace modeling and multi-objective optimization
CN107023825A (en) * 2016-08-31 2017-08-08 西安艾贝尔科技发展有限公司 Fluidized-bed combustion boiler is controlled and combustion optimizing system
CN109507910A (en) * 2018-10-31 2019-03-22 国网河南省电力公司电力科学研究院 A kind of fired power generating unit modeling and control method based on data-driven

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102183015A (en) * 2011-04-01 2011-09-14 北京三博中自科技有限公司 Combustion optimizing control system of circulating fluidized bed boiler in wide-range change of load
CN102425807A (en) * 2011-11-23 2012-04-25 华北电力大学(保定) Combustion feedforward and feedback composite optimization controlling method for pulverized coal fired boiler
CN103576655A (en) * 2013-11-06 2014-02-12 华北电力大学(保定) Method and system for utility boiler combustion subspace modeling and multi-objective optimization
CN107023825A (en) * 2016-08-31 2017-08-08 西安艾贝尔科技发展有限公司 Fluidized-bed combustion boiler is controlled and combustion optimizing system
CN109507910A (en) * 2018-10-31 2019-03-22 国网河南省电力公司电力科学研究院 A kind of fired power generating unit modeling and control method based on data-driven

Also Published As

Publication number Publication date
CN111623369A (en) 2020-09-04

Similar Documents

Publication Publication Date Title
CN111623369B (en) Control method for adjusting boiler fuel feeding quantity by using smoke oxygen content signal
CN107561941B (en) Full-working-condition standard-reaching emission control method for thermal power generating unit denitration system
US8554706B2 (en) Power plant control device which uses a model, a learning signal, a correction signal, and a manipulation signal
CN103019097B (en) Optimal control system for steel rolling heating furnace
CN110263452B (en) Flue gas time distribution characteristic analysis method and system in flue and denitration system
CN111881554B (en) Optimization control method for boiler changing along with air temperature
WO2009080282A2 (en) Optimizing the operation of a power plant
CN112628712A (en) Secondary air closed-loop optimization control system based on air door resistance coefficient
CN113033873B (en) Method for predicting content of sulfur oxides at inlet of desulfurization system based on measurement of coal quality entering furnace
CN110762540A (en) Combustion decision method and automatic combustion control method and system for garbage incinerator
CN115111594A (en) Intelligent regulation and control system and method for heat accumulating type thermal oxidation furnace
CN111061149A (en) Circulating fluidized bed coal saving and consumption reduction method based on deep learning prediction control optimization
CN111457392B (en) Intelligent setting method for air quantity in urban domestic garbage incineration process
CN114721263B (en) Intelligent regulation and control method for cement decomposing furnace based on machine learning and intelligent optimization algorithm
CN104729888B (en) Flue fly ash constant-speed sampling system based on BP neural network and control method
CN115685743A (en) Intelligent control coal-fired boiler and intelligent prediction regulation and control flue gas emission method thereof
CN108488831B (en) Boiler combustion control system and method
CN104561514A (en) Double heat accumulating type steel rolling heating furnace oxidizing atmosphere adjustment method and automatic control method thereof
CN109882883B (en) Method and device for optimizing boiler coal burning efficiency based on artificial intelligence
CN115198047B (en) Hot blast stove combustion monitoring system and method based on big data analysis
CN111855953A (en) Method and system for quickly feeding back calorific value of coal as fired in coal-fired power plant
CN114187976A (en) Oxygen amount prediction method based on real-time measurement data
CN115145152A (en) Boiler combustion and denitration process collaborative optimization control method
CN211526443U (en) Automatic combustion control system of garbage incinerator
CN111765478A (en) Automatic control method for primary air volume of garbage incinerator

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220125