CN101881563A - Multi-area intelligent online optimizing control method for thermal efficiency of heating furnace - Google Patents
Multi-area intelligent online optimizing control method for thermal efficiency of heating furnace Download PDFInfo
- Publication number
- CN101881563A CN101881563A CN 201010217097 CN201010217097A CN101881563A CN 101881563 A CN101881563 A CN 101881563A CN 201010217097 CN201010217097 CN 201010217097 CN 201010217097 A CN201010217097 A CN 201010217097A CN 101881563 A CN101881563 A CN 101881563A
- Authority
- CN
- China
- Prior art keywords
- oxygen content
- thermal efficiency
- heating furnace
- combustion chamber
- 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.)
- Granted
Links
- 238000010438 heat treatment Methods 0.000 title claims abstract description 80
- 238000000034 method Methods 0.000 title claims abstract description 63
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims abstract description 94
- 229910052760 oxygen Inorganic materials 0.000 claims abstract description 94
- 239000001301 oxygen Substances 0.000 claims abstract description 94
- 238000005457 optimization Methods 0.000 claims abstract description 49
- 230000008569 process Effects 0.000 claims abstract description 38
- 238000002485 combustion reaction Methods 0.000 claims description 68
- 239000000446 fuel Substances 0.000 claims description 50
- OIGNJSKKLXVSLS-VWUMJDOOSA-N prednisolone Chemical compound O=C1C=C[C@]2(C)[C@H]3[C@@H](O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 OIGNJSKKLXVSLS-VWUMJDOOSA-N 0.000 claims description 29
- 239000007789 gas Substances 0.000 claims description 21
- 230000004044 response Effects 0.000 claims description 21
- 230000008859 change Effects 0.000 claims description 17
- 238000012360 testing method Methods 0.000 claims description 15
- 239000002737 fuel gas Substances 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 9
- 230000009123 feedback regulation Effects 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 235000019504 cigarettes Nutrition 0.000 claims description 3
- 239000000470 constituent Substances 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 230000007115 recruitment Effects 0.000 claims description 3
- 229910052717 sulfur Inorganic materials 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 claims description 2
- 238000009412 basement excavation Methods 0.000 abstract 1
- 239000000779 smoke Substances 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 6
- 230000007774 longterm Effects 0.000 description 4
- 230000001276 controlling effect Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000001105 regulatory effect Effects 0.000 description 3
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 2
- 206010013082 Discomfort Diseases 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000003546 flue gas Substances 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000004939 coking Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005292 vacuum distillation Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Landscapes
- Regulation And Control Of Combustion (AREA)
Abstract
The invention relates to a multi-area intelligent online optimizing control method for thermal efficiency of a heating furnace, which belongs to the technical field of the control of the thermal efficiency of heating furnaces. The method is characterized in that a thermal efficiency optimized off-line part carries out working area division on the working condition of a heating furnace according to a thermal load and obtains optimizing values of oxygen content and negative pressure in each working area by historical data excavation. An online part confirms the working area of the heating furnace according to process real-time data, can be maintained in a better working state in time when the load is changed by starting with the optimizing values of the oxygen content and the negative pressure obtained by the off-line part and carries out thermal efficiency self-optimizing control when operation is steady, thereby achieving the aims of rapid optimization and long-time working near an optimizing state. A smoke air system carries out area control on the oxygen content and the negative pressure of a hearth by adopting a control method based on dynamic feedforward and steady feedback, and can ensure that the heating furnace works near an optimized working point given by a thermal efficiency optimizing method so as to realize the optimal thermal efficiency.
Description
Technical field
The present invention relates to the on-line optimizing and controlling method of thermal efficiency of heating furnace, belong to petrochemical industry tubular heater and production process automation field.
Background technology
Tubular heater (to call heating furnace in the following text) is the important technological equipment of production process extensive uses such as oil refining and petrochemical industry, almost in every suit commercial plant heating furnace is arranged all.Heating furnace provides thermal source for device, is the capital equipment of device energy-wasting.For example, the atmospheric and vacuum distillation fuel consumption accounts for the device total energy consumption at 82-92%, and delayed coking is about 90%.Heating furnace fuel can produce CO in combustion process
2, gases such as CO.Heating furnace to operation carries out thermal efficiency optimization, makes its long-term work near optimum combustion regime, and under the situation that does not change technology, realize saving energy and reduce the cost, and reduce the environmental pollution that causes because of imperfect combustion, be a urgent task.
The thermal efficiency of heating furnace is the percentage that the effective heat that utilizes of heated medium accounts for the aggregate supply heat.For heating furnace, if air quantity is too small, can cause fuel combustion insufficient, waste fuel also produces gaseous contamination environment such as CO; If air quantity is excessive, it is too much then to discharge the heat that flue gas takes away, and the thermal efficiency reduces.Therefore, there is the optimum thermal efficiency in corresponding certain operating mode.The adjustment means of thermal efficiency of heating furnace are air quantity, by adjusting air quantity, keep rational excess air coefficient, reach optimum thermal efficiency.
Control and optimization about thermal efficiency of heating furnace have had some researchs and application result.The at present industrial method that more generally adopts is that oxygen content is controlled efficiency of combustion indirectly in the control flue gas, but load, when fuel changes, respective change also can take place in optimal oxygen content.The bibliographical information internal model control and the thermal efficiency are used to the real-time control and the optimization of thermal efficiency of heating furnace from methods such as optimizing.But the control to the thermal efficiency can not guarantee optimum thermal efficiency, also exists the controllability problem.Adopt the method for the thermal efficiency from optimizing, when the heating furnace thermic load changed greatly, optimizing speed was slower.
PID control is generally adopted in the control of industrial air and gas system.Owing to the heating furnace air quantity is regulated reasons such as executing agency's sensitivity is relatively poor, the dead band is big, general discomfort is fit to do continuous adjusting, and above-mentioned control loop major part is in manual state for a long time.Owing to the heating furnace air quantity is regulated reasons such as executing agency's sensitivity is relatively poor, the dead band is big, general discomfort is fit to do continuous adjusting, and above-mentioned control loop is in manual state for a long time, let alone the optimization of the thermal efficiency.Therefore, need the automatic control of the suitable control strategy realization of exploitation, for the optimization that realizes the thermal efficiency lays the foundation to the heating furnace air and gas system.
Summary of the invention
Purpose of the present invention: the multi-area intelligent online optimizing control method for thermal efficiency of heating furnace that provides a kind of practicality.Method comprises thermal efficiency optimal control and air and gas system control two parts.Thermal efficiency optimal control partly is divided into off-line and online two parts.The off-line part is carried out the working region division according to thermic load to the heating furnace operating mode, excavates by historical data, obtains the optimization operating point in each working region.Online part is determined the working region of heating furnace according to the process real time data, carries out the thermal efficiency then from seeking optimum control in each working region, reaches near quick optimizing and the long-term work purpose the optimization state.The control method based on " dynamic Feedforward, stable state feedback " is adopted in the control of air and gas system, and oxygen content and combustion chamber draft are controlled, and guarantees that heating furnace has stable fired state, for optimization lays the foundation.
The invention is characterized in: described method realizes in host computer successively according to the following steps:
Steps A: host computer initialization:
In described host computer, set up with lower module: on-line optimization module, air and gas system control module and real-time data base/OPC bitcom module, wherein:
The on-line optimization module, pass through the real time data that the OPC bitcom is gathered heating furnace from heating furnace controlled device and Distributed Control System by described operation control, and described real-time data base is sent in the current optimization operating point of the history in the thermic load zone of heating furnace optimization operating point and thermic load uses for described air and gas system control module;
The air and gas system control module, under the effect of described OPC bitcom, gather numerical value such as oxygen content and combustion chamber draft in real time, and read the determined optimization of on-line optimization module operating point in the described real-time data base, the combustion chamber draft in the air-fuel ratio when fuel is increased, thermic load zone oxygen content and thermic load zone is controlled, oxygen content and combustion chamber draft are maintained optimize near the operating point, the control action of being calculated is sent to described heating furnace controlled device and Distributed Control System by described OPC bitcom;
Step B: described on-line optimization module, carry out thermal efficiency offline optimization and thermal efficiency on-line optimization successively according to the following steps:
Step B1: thermal efficiency offline optimization, its step is as follows:
Step B1.1: step test obtains the steady-state response time of the thermal efficiency:
Under the steady situation of process, the oxygen content setting value is applied the step test signal, record thermal efficiency change curve, the steady-state response time T of the acquisition thermal efficiency
r
Step B1.2: set sampling period T, T ∈ [0.25T
r, T
r], and gather following heating furnace field data: heating furnace oxygen content, heating furnace burner hearth negative pressure, furnace outlet temperature, heating furnace inlet temperature, furnace charge flow, fuel flow rate and intake, and carry out the off-line modeling analysis according to the following steps;
Step B1.3: heated medium is not had the heating furnace of phase transformation, calculate the effective heat duty Q (k) of heating furnace:
Q(k)=F(k)C
p[T
out(k)-T
in(k)]
Wherein: k is sampling instant,
Q (k) is the k effective heat duty of heating furnace constantly,
F (k) is the k flow of heated medium constantly,
T
Out(k) be the k outlet temperature of heated medium constantly,
T
In(k) be the k inlet temperature of heated medium constantly,
C
pSpecific heat for heated medium;
Step B1.4: utilize thermic load territorial classification device to be divided into N zone according to the residing working region of big young pathbreaker's heating furnace thermic load of thermic load, as N=5, zone Ω
iExpression, i=1,2 ..., 5;
The working range of loading when the heating furnace normal heat is designing effective heat duty Q
00.75 times when between 1.25 times, changing, the zone in 5 zones limit is respectively:
[0.75Q
0 0.85Q
0),[0.85Q
0 0.95Q
0),[0.95Q
0 1.05Q
0),
[1.05Q
0 1.15Q
0),[1.15Q
0 1.25Q
0]
Use Ω
IL, Ω
IHRepresent regional Ω respectively
iThe lower limit and the upper limit;
Step B1.5: calculate the thermal efficiency according to the positive balance method:
Wherein: η (k) calculates thermal efficiency value constantly for k,
F
f(k) be that k is fuel flow rate constantly,
H
fBe the fuel combustion calorific value;
Step B1.6:, seek thermic load working region Ω according to the historical data that obtains among the step B1.2
iThe history of internally heated oven is optimized the operating point:
Wherein: η
i(k) be thermic load working region Ω
iInterior k calculates thermal efficiency value constantly,
O
2i(k) be the corresponding constantly measurement of oxygen content value of k,
P
i(k) be the corresponding constantly combustion chamber draft measured value of k;
The oxygen content that makes thermal efficiency maximum and combustion chamber draft value as described thermic load working region Ω
iThe history of internally heated oven is optimized the operating point, uses Ω
I, optExpression:
Ω
i,opt={O
2i,opt P
i,opt} i=1,2,…,5
Wherein: Ω
I, optBe thermic load working region Ω
iThe historical operating point of optimizing,
O
2i, optBe thermic load working region Ω
iThe historical optimal value of interior oxygen content,
P
I, optBe thermic load working region Ω
iThe historical optimal value of internal furnace negative pressure;
Step B2: thermal efficiency on-line optimization, its step is as follows:
Step B2.1: the real time data of gatherer process comprises: heating furnace oxygen content, heating furnace burner hearth negative pressure, furnace outlet temperature, heating furnace inlet temperature, furnace charge flow and fuel flow rate;
Step B2.2: utilize the thermic load working region grader deterministic process present located thermic load working region Ω among the step B1.4
i, and the history of obtaining corresponding calculated off-line gained is optimized point value Ω
I, opt
Step B2.3: the described method of B1.5 is calculated current time thermal efficiency η (k) set by step;
Step B2.4: by following decision criteria deterministic process whether for stable state:
Wherein: y
Jl(j=1,2,3) are respectively l value of characteristic variable (furnace outlet temperature, fire box temperature and inlet amount),
L is whether deterministic process is in stable historical data length, and L*T=30min, T are the sampling periods,
ε is preassigned stable state decision threshold, and span is (0,0.1);
Step B2.5: determine to optimize under the current thermic load operating point according to the following steps:
Step B2.5.1: if process is in stable state, then the online optimizing of the thermal efficiency is adopted from seeking method for optimally controlling, is the tuning variable with the oxygen content setting value, optimizes operating point O with history
2i, optBe initial value, online searching makes the highest oxygen content setting value O of the thermal efficiency
2s(k), step is as follows:
Step B2.5.1.1: the changing value Δ η (k) that calculates the thermal efficiency of going up a moment k-1 relatively:
Δη(k)=η(k)-η(k-1);
Step B2.5.1.2: if | Δ η (k) |<Δ η
Min, then stop optimizing and record oxygen content setting value at this moment, wherein Δ η
MinFor the default thermal efficiency is adjusted the dead band;
Step B2.5.1.3: if | Δ η (k) | 〉=Δ η
Min, then with Δ O
2s(k) be described oxygen content setting value O
2s(k) optimizing step-length automatic optimal:
O
2s(k)=O
2s(k-1)+ΔO
2s(k)
Wherein: λ is for adjusting coefficient, span be (0,1];
Step B2.5.2:, optimize operating point O with history if process is not in stable state
2i, optBe setting value O
2s(k), guarantee that heating furnace is in the duty of suboptimum;
Step B2.6: the oxygen content of calculating setting value O
2s(k) be sent to described real-time data base, come piece to implement for the Control for Oxygen Content device in the described air and gas system control module;
Step C: the control object of described air and gas system control module comprises the air-fuel ratio when oxygen content, hearth load and fuel increase, the control target is that described oxygen content and combustion chamber draft are controlled in the optimization set point or given range that is drawn jointly by described offline optimization and on-line optimization, and step is as follows:
Step C1: step test obtains the steady-state response time of oxygen content and combustion chamber draft:
Under the steady situation of process, intake is applied the step test signal, record oxygen content change curve, the steady-state response time T of acquisition oxygen content
O2ss
Under the steady situation of process, the blower inlet baffle plate is applied the step test signal, record combustion chamber draft change curve, the steady-state response time T of acquisition combustion chamber draft
Pss
Step C2: the control cycle of setting cigarette wind control system module is T
c, T
c=min (T
O2ss, T
Pss)/40;
Step C3: when fuel increases air-fuel ratio is controlled according to the following steps:
Step C3.1: if fuel recruitment Δ F
f(k
c) last relatively intake adjustment moment k
c-1 fuel quantity F
F0Ratio surpasses predetermined threshold value β, and β ∈ (0,0.2], promptly
The feedforward variation delta F of air intake when then fuel increases
AF(k
c) be
ΔF
aF(k
c)=α·AFR·ΔF
f(k
c)
Wherein: Δ F
f(k
c)=F
f(k
c)-F
F0Be fuel change amount, F
f(k
c) be fuel flow rate,
α is an excess air coefficient, is [1.05,1.15] to the fuel gas span;
AFR is a stoichiometric air/, and to fuel gas, it is calculated as:
Wherein: CO, H
2, C
mH
n, H
2S, O
2Being each constituent content in the fuel gas, is unit with %;
Step C4: according to the following steps oxygen content is carried out Region control:
Step C4.1: if the setting value of oxygen content is O
2s, the regional extent of its permission is:
[O
2s-δ?O
2s+δ]
Wherein: δ is that oxygen content departs from setting value O
2sZone limit, δ ∈ (00.5);
Step C4.2: at Control for Oxygen Content moment k
c,
If O
2(k
c)>O
2HPerhaps O
2(k
c)<O
2L, O wherein
2H=O
2s+ δ, O
L=O
2s-δ, and described Control for Oxygen Content device is not in the stand-by period, then is calculated as follows required intake and changes:
Wherein: Δ f
a(k
c) be the variable quantity of intake,
O
2sBe the setting value of oxygen content,
O
2(k
c) be the measured value of oxygen content;
If O
2L≤ O
2(k
c)≤O
2H, Δ F then
a(k
c)=0;
According to Control for Oxygen Content moment k
cCombustion chamber draft P (k
c) by following principle correction intake variation delta F
a(k
c):
If P (k
c)>P
HS, P
HSBe the upper safety limit of combustion chamber draft, and the variation delta F of intake
a(k
c)>0 then keeps intake constant, makes Δ F
a(k
c)=0,
If P (k
c)<P
Ls, P
LsBe the lower safety limit of combustion chamber draft, and the variation delta F of intake
a(k
c)<0 then keeps intake constant, makes Δ F
a(k
c)=0;
Step C4.3: be calculated as follows required intake:
F
a(k
c)=F
a(k
c-1)+ΔF
a(k
c)
Wherein: F
a(k
c) be intake;
After process made feedback regulation, wait for the steady-state response time T of an oxygen content
O2ss
Step C5: according to the following steps combustion chamber draft is carried out Region control:
Step C5.1: if the setting value of combustion chamber draft is P
s, and be in thermic load working region Ω
i, P then
s=P
I, opt, i=1,2 ..., 5, it allows the regional extent of change to be:
[P
s-σ?P
s+σ]
Wherein: σ is the zone limit that combustion chamber draft departs from setting value, δ ∈ (010);
Step C5.2: at combustion chamber draft control moment k
c,
If P (k
c)>P
HPerhaps P (k
c)<P
L, P wherein
H=P
s+ σ, P
L=P
s-σ, and described combustion chamber draft controller is not in the stand-by period, and then the deviation of combustion chamber draft is e
P(k
c)=P
s-P (k
c), air-introduced machine inlet baffle variation delta MV then
2(k
c) be
Wherein: K
1Proportionality coefficient for air-introduced machine inlet baffle and negative pressure variation;
After process made feedback regulation, wait for the steady-state response time T of a combustion chamber draft
Pss
Under other situations, Δ MV
2(k
c)=0;
Step C5.3: FEEDFORWARD CONTROL is carried out in the variation of blower variable frequency by following formula:
ΔMV
2F(k
c)=K
2ΔMV
1(k
c)
Wherein: Δ MV
2F(k
c) be the feedforward variation of air-introduced machine inlet baffle,
Δ MV
1(k
c)=MV
1(k
c)-MV
1(k
c-1) variation of exporting for blower variable frequency,
K
2Be the feed-forward coefficients between baffle plate variation and the frequency conversion output.
Use proof: the heating furnace offline optimization carries out the working region according to thermic load to the heating furnace operating mode to be divided, and excavates by historical data, can access the optimization operating point in each working region.On-line optimization is determined the working region of heating furnace according to the process real time data, and optimizing the operating point with history in each working region is that starting point is carried out the thermal efficiency from seeking optimum control, reaches near quick optimizing and the long-term work purpose the optimization state.The control of air and gas system can effectively be removed the dangerous influence that the intake continuous closed-loop is regulated based on the control method of " stable state feedback, dynamic Feedforward " thought, removes the stability influence of Dynamic Coupling to control.
Description of drawings
Fig. 1 is heating furnace flow process and thermal efficiency multi-area intelligent optimal control structural representation.
Fig. 2 is that the thermal efficiency is optimized off-line part steps schematic diagram.
Fig. 3 is that the thermal efficiency is optimized online optimizing and air and gas system control general steps schematic diagram.
Fig. 4 is that the thermal efficiency is optimized thermic load territorial classification device.
Fig. 5 is a thermal efficiency on-line optimization program flow diagram.
The specific embodiment
Below in conjunction with accompanying drawing and case study on implementation, the specific embodiment of the present invention is described in further detail.Wherein, Fig. 1 is heating furnace flow process and thermal efficiency multi-area intelligent optimal control structural representation; Fig. 2 is that the thermal efficiency is optimized off-line part steps schematic diagram; Fig. 3 is that the thermal efficiency is optimized online optimizing and air and gas system control general steps schematic diagram; Fig. 4 is that the thermal efficiency is optimized thermic load territorial classification device; Fig. 5 is a thermal efficiency on-line optimization program flow diagram.
Thermal efficiency offline optimization
Step B1.1: step test obtains the steady-state response time of the thermal efficiency:
Under the steady situation of process, the oxygen content setting value is applied the step test signal, record thermal efficiency change curve, the steady-state response time T of the acquisition thermal efficiency
r
Step B1.2: set sampling period T, T ∈ [0.25T
r, T
r], and gather following heating furnace field data: heating furnace oxygen content, heating furnace burner hearth negative pressure, furnace outlet temperature, heating furnace inlet temperature, furnace charge flow, fuel flow rate and intake;
Gathering above data a period of time (as one month) is used for the off-line modeling analysis;
Step B1.3: heated medium is not had the heating furnace of phase transformation, calculate the effective heat duty Q (k) of heating furnace:
Q(k)=F(k)C
p[T
out(k)-T
in(k)]
Wherein: k is sampling instant,
Q (k) is the k effective heat duty of heating furnace constantly,
F (k) is the k flow of heated medium constantly,
T
Out(k) be the k outlet temperature of heated medium constantly,
T
In(k) be the k inlet temperature of heated medium constantly,
C
pSpecific heat for heated medium;
Step B1.4: utilize thermic load territorial classification device to be divided into N zone according to the residing working region of big young pathbreaker's heating furnace thermic load of thermic load, as N=5, zone Ω
iExpression, i=1,2 ..., 5;
The working range of loading when the heating furnace normal heat is designing effective heat duty Q
00.75 times when between 1.25 times, changing, the zone in 5 zones limit is respectively:
[0.75Q
0?0.85Q
0),[0.85Q
0?0.95Q
0),[0.95Q
0?1.05Q
0),
[1.05Q
0?1.15Q
0),[1.15Q
0?1.25Q
0]
Use Ω
IL, Ω
IHRepresent regional Ω respectively
iThe lower limit and the upper limit;
Step B1.5: calculate the thermal efficiency according to the positive balance method:
Wherein: η (k) calculates thermal efficiency value constantly for k,
F
f(k) be that k is fuel flow rate constantly,
H
fBe the fuel combustion calorific value;
Step B1.6:, seek thermic load working region Ω according to the historical data that obtains among the step B1.2
iThe history of internally heated oven is optimized the operating point:
Wherein: η
i(k) be thermic load working region Ω
iInterior k calculates thermal efficiency value constantly,
O
2i(k) be the corresponding constantly measurement of oxygen content value of k,
P
i(k) be the corresponding constantly combustion chamber draft measured value of k;
The oxygen content that makes thermal efficiency maximum and combustion chamber draft value as described thermic load working region Ω
iThe history of internally heated oven is optimized the operating point, uses Ω
I, optExpression:
Ω
i,opt={O
2i,opt?P
i,opt} i=1,2,…,5
Wherein: Ω
I, optBe thermic load working region Ω
iThe historical operating point of optimizing,
O
2i, optBe thermic load working region Ω
iThe historical optimal value of interior oxygen content,
P
I, optBe thermic load working region Ω
iThe historical optimal value of internal furnace negative pressure.
Thermal efficiency on-line optimization
In host computer, the optimization timer is set, the cycle is T.The thermal efficiency is optimized online part and is carried out the following step at each optimization cycle:
Step B2.1: the real time data of gatherer process comprises: heating furnace oxygen content, heating furnace burner hearth negative pressure, furnace outlet temperature, heating furnace inlet temperature, furnace charge flow and fuel flow rate;
Step B2.2: utilize the thermic load working region grader deterministic process present located thermic load working region Ω among the step B1.4
i, and the history of obtaining corresponding calculated off-line gained is optimized point value Ω
I, opt
Step B2.3: the described method of B1.5 is calculated current time thermal efficiency η (k) set by step;
Step B2.4: by following decision criteria deterministic process whether for stable state:
Wherein: y
Jl(j=1,2,3) are respectively l value of characteristic variable (furnace outlet temperature, fire box temperature and inlet amount),
L is whether deterministic process is in stable historical data length, and L*T=30min, T are the sampling periods,
Be the mean value of j characteristic variable of selection,
ε is preassigned stable state decision threshold, and span is (0,0.1);
Step B2.5: determine to optimize under the current thermic load operating point according to the following steps:
Step B2.5.1: if process is in stable state, then the online optimizing of the thermal efficiency is adopted from seeking method for optimally controlling, is the tuning variable with the oxygen content setting value, optimizes operating point O with history
2i, optBe initial value, online searching makes the highest oxygen content setting value O of the thermal efficiency
2s(k), step is as follows:
Step B2.5.1.1: the changing value Δ η (k) that calculates the thermal efficiency of going up a moment k-1 relatively:
Δη(k)=η(k)-η(k-1);
Step B2.5.1.2: if | Δ η (k) |<Δ η
Min, then stop optimizing and record oxygen content setting value at this moment, wherein Δ η
MinFor the default thermal efficiency is adjusted the dead band;
Step B2.5.1.3: if | Δ η (k) | 〉=Δ η
Min, then with Δ O
2s(k) be described oxygen content setting value O
2s(k) optimizing step-length automatic optimal:
O
2s(k)=O
2s(k-1)+ΔO
2s(k)
Wherein: λ is for adjusting coefficient, span be (0,1];
Step B2.5.2:, optimize operating point O with history if process is not in stable state
2i, optBe setting value O
2s(k), guarantee that heating furnace is in the duty of suboptimum;
Step B2.6: the oxygen content of calculating setting value O
2s(k) be sent to described real-time data base, come piece to implement for the Control for Oxygen Content device in the described air and gas system control module;
By the classification of thermic load multizone, reach near the purpose the optimization state at quick optimizing of specific operation and long-term work, especially be fit to the heating furnace that thermic load often changes.
Air and gas system control
The control of air and gas system comprises air-fuel ratio control, oxygen content Region control and the combustion chamber draft Region control when fuel increases.Its target is that oxygen content and combustion chamber draft are controlled in the optimization set point or given range that above-mentioned offline optimization and on-line optimization determine.
Step C1: step test obtains the steady-state response time of oxygen content and combustion chamber draft:
Under the steady situation of process, intake is applied the step test signal, record oxygen content change curve, the steady-state response time T of acquisition oxygen content
O2ss
Under the steady situation of process, the blower inlet baffle plate is applied the step test signal, record combustion chamber draft change curve, the steady-state response time T of acquisition combustion chamber draft
Pss
Step C2: the control cycle of setting cigarette wind control system module is T
c, T
c=min (T
O2ss, T
Pss)/40;
Step C3: the air-fuel ratio control when fuel increases.The purpose of air-fuel ratio control is to carry out dynamic Feedforward when fuel increases, and prevents that burning is not enough.
Step C3.1: if fuel recruitment Δ F
f(k
c) last relatively intake adjustment moment k
c-1 fuel quantity F
F0Ratio surpasses predetermined threshold value β, and β ∈ (0,0.2], promptly
Wherein: Δ F
f(k
c)=F
f(k
c)-F
F0Be fuel change amount, F
f(k
c) be fuel flow rate;
The feedforward variation delta F of air intake when then fuel increases
AF(k
c) be
ΔF
aF(k
c)=α·AFR·ΔF
f(k
c)
Wherein: α is an excess air coefficient, is [1.05,1.15] to the fuel gas span;
AFR is a stoichiometric air/, and it is calculated as to fuel gas:
Wherein: CO, H
2, C
mH
n, H
2S, O
2Being each constituent content in the fuel gas, is unit with %;
Step C4: oxygen content Region control:
The controlled variable of oxygen content Region control is an oxygen content, and performance variable is intake (or calculating intake).The target of oxygen content Region control is with near the zone of Control for Oxygen Content setting value.
Step C4.1: if the setting value of oxygen content is O
2s, the regional extent of its permission is:
[O
2s-δ?O
2s+δ]
Wherein: δ is that oxygen content departs from setting value O
2sZone limit, δ ∈ (00.5);
Step C4.2: at Control for Oxygen Content moment k
c,
If O
2(k
c)>O
2HPerhaps O
2(k
c)<O
2L, O wherein
2H=O
2s+ δ, O
L=O
2s-δ, and described Control for Oxygen Content device is not in the stand-by period, then is calculated as follows required intake and changes:
Wherein: Δ F
a(k
c) be the variable quantity of intake,
O
2sBe the setting value of oxygen content,
O
2(k
c) be the measured value of oxygen content;
If O
2L≤ O
2(k
c)≤O
2H, Δ F then
a(k
c)=0;
According to Control for Oxygen Content moment k
cCombustion chamber draft P (k
c) by following principle correction intake variation delta F
a(k
c):
If P (k
c)>P
HS, P
HSBe the upper safety limit of combustion chamber draft, and the variation delta F of intake
a(k
c)>0 then keeps intake constant, makes Δ F
a(k
c)=0,
If P (k
c)<P
LS, P
LsBe the lower safety limit of combustion chamber draft, and Δ F
a(k
c)<0 then keeps intake constant, makes Δ F
a(k
c)=0;
Step C4.3: be calculated as follows required intake:
F
a(k
c)=F
a(k
c-1)+ΔF
a(k
c)
Wherein: F
a(k
c) be intake;
After process made feedback regulation, wait for the steady-state response time T of an oxygen content
O2ss
Step C5: combustion chamber draft Region control:
The controlled variable of combustion chamber draft Region control is a combustion chamber draft, and performance variable is the blower inlet baffle plate.The target of combustion chamber draft Region control is that combustion chamber draft is controlled near the zone of setting value.
Step C5.1: if the setting value of combustion chamber draft is P
s, and be in thermic load working region Ω
i, P then
s=P
I, opt, i=1,2 ..., 5, it allows the regional extent of change to be:
[P
s-σ?P
s+σ]
Wherein: σ is the zone limit that combustion chamber draft departs from setting value, δ ∈ (010);
Step C5.2: at combustion chamber draft control moment k
c,
If P (k
c)>P
HPerhaps P (k
c)<P
L, P wherein
H=P
s+ σ, P
L=P
s-σ, and described combustion chamber draft controller is not in the stand-by period, and then the deviation of combustion chamber draft is e
P(k
c)=P
s-P (k
c), air-introduced machine inlet baffle variation delta MV then
2(k
c) be
Wherein: K
1Proportionality coefficient for air-introduced machine inlet baffle and negative pressure variation;
After process made feedback regulation, wait for the steady-state response time T of a combustion chamber draft
Pss
Under other situations, Δ MV
2(k
c)=0;
Step C5.3: FEEDFORWARD CONTROL is carried out in the variation of blower variable frequency by following formula:
ΔMV
2F(k
c)=K
2ΔMV
1(k
c)
Wherein: Δ MV
2F(k
c) be the feedforward variation of air-introduced machine inlet baffle,
Δ MV
1(k
c)=MV
1(k
c)-MV
1(k
c-1) variation of exporting for blower variable frequency,
K
2Be the feed-forward coefficients between baffle plate variation and the frequency conversion output.
Claims (1)
1. multi-area intelligent online optimizing control method for thermal efficiency of heating furnace is characterized in that, described method realizes in host computer successively according to the following steps:
Steps A: host computer initialization:
In described host computer, set up with lower module: on-line optimization module, air and gas system control module and real-time data base/OPC bitcom module, wherein:
The on-line optimization module, pass through the real time data that the OPC bitcom is gathered heating furnace from heating furnace controlled device and Distributed Control System by described operation control, and described real-time data base is sent in the current optimization operating point of the history in the thermic load zone of heating furnace optimization operating point and thermic load uses for described air and gas system control module;
The air and gas system control module, under the effect of described OPC bitcom, gather numerical value such as oxygen content and combustion chamber draft in real time, and read the determined optimization of on-line optimization module operating point in the described real-time data base, the combustion chamber draft in the air-fuel ratio when fuel is increased, thermic load zone oxygen content and thermic load zone is controlled, oxygen content and combustion chamber draft are maintained optimize near the operating point, the control action of being calculated is sent to described heating furnace controlled device and Distributed Control System by described OPC bitcom;
Step B: described on-line optimization module, carry out thermal efficiency offline optimization and thermal efficiency on-line optimization successively according to the following steps:
Step B1: thermal efficiency offline optimization, its step is as follows:
Step B1.1: step test obtains the steady-state response time of the thermal efficiency:
Under the steady situation of process, the oxygen content setting value is applied the step test signal, record thermal efficiency change curve, the steady-state response time T of the acquisition thermal efficiency
r
Step B1.2: set sampling period T, T ∈ [0.25T
r, T
r], and gather following heating furnace field data: heating furnace oxygen content, heating furnace burner hearth negative pressure, furnace outlet temperature, heating furnace inlet temperature, furnace charge flow, fuel flow rate and intake, and carry out the off-line modeling analysis according to the following steps;
Step B1.3: heated medium is not had the heating furnace of phase transformation, calculate the effective heat duty Q (k) of heating furnace:
Q(k)=F(k)C
p[T
out(k)-T
in(k)]
Wherein: k is sampling instant,
Q (k) is the k effective heat duty of heating furnace constantly,
F (k) is the k flow of heated medium constantly,
T
Out(k) be the k outlet temperature of heated medium constantly,
T
In(k) be the k inlet temperature of heated medium constantly,
C
pSpecific heat for heated medium;
Step B1.4: utilize thermic load territorial classification device to be divided into N zone according to the residing working region of big young pathbreaker's heating furnace thermic load of thermic load, as N=5, zone Ω
iExpression, i=1,2 ..., 5;
The working range of loading when the heating furnace normal heat is designing effective heat duty Q
00.75 times when between 1.25 times, changing, the zone in 5 zones limit is respectively:
[0.75Q
0?0.85Q
0),[0.85Q
0?0.95Q
0),[0.95Q
0?1.05Q
0),
[1.05Q
0?1.15Q
0),[1.15Q
0?1.25Q
0]
Use Ω
IL, Ω
IHRepresent regional Ω respectively
iThe lower limit and the upper limit;
Step B1.5: calculate the thermal efficiency according to the positive balance method:
Wherein: η (k) calculates thermal efficiency value constantly for k,
F
f(k) be that k is fuel flow rate constantly,
H
fBe the fuel combustion calorific value;
Step B1.6:, seek thermic load working region Ω according to the historical data that obtains among the step B1.2
iThe history of internally heated oven is optimized the operating point:
Wherein: η
i(k) be thermic load working region Ω
iInterior k calculates thermal efficiency value constantly,
O
2i(k) be the corresponding constantly measurement of oxygen content value of k,
P
i(k) be the corresponding constantly combustion chamber draft measured value of k;
The oxygen content that makes thermal efficiency maximum and combustion chamber draft value as described thermic load working region Ω
iThe history of internally heated oven is optimized the operating point, uses Ω
I, optExpression:
Ω
i,opt={O
2i,opt?P
i,opt} i=1,2,…,5
Wherein: Ω
I, optBe thermic load working region Ω
iThe historical operating point of optimizing,
O
2i, optBe thermic load working region Ω
iThe historical optimal value of interior oxygen content,
P
I, optBe thermic load working region Ω
iThe historical optimal value of internal furnace negative pressure;
Step B2: thermal efficiency on-line optimization, its step is as follows:
Step B2.1: the real time data of gatherer process comprises: heating furnace oxygen content, heating furnace burner hearth negative pressure, furnace outlet temperature, heating furnace inlet temperature, furnace charge flow and fuel flow rate;
Step B2.2: utilize the thermic load working region grader deterministic process present located thermic load working region Ω among the step B1.4
i, and the history of obtaining corresponding calculated off-line gained is optimized point value Ω
I, opt
Step B2.3: the described method of B1.5 is calculated current time thermal efficiency η (k) set by step;
Step B2.4: by following decision criteria deterministic process whether for stable state:
Wherein: y
Jl(j=1,2,3) are respectively l value of characteristic variable (furnace outlet temperature, fire box temperature and inlet amount),
L is whether deterministic process is in stable historical data length, and L*T=30min, T are the sampling periods,
ε is preassigned stable state decision threshold, and span is (0,0.1);
Step B2.5: determine to optimize under the current thermic load operating point according to the following steps:
Step B2.5.1: if process is in stable state, then the online optimizing of the thermal efficiency is adopted from seeking method for optimally controlling, is the tuning variable with the oxygen content setting value, optimizes operating point O with history
2i, optBe initial value, online searching makes the highest oxygen content setting value O of the thermal efficiency
2s(k), step is as follows:
Step B2.5.1.1: the changing value Δ η (k) that calculates the thermal efficiency of going up a moment k-1 relatively:
Δη(k)=η(k)-η(k-1);
Step B2.5.1.2: if | Δ η (k) |<Δ η
Min, then stop optimizing and record oxygen content setting value at this moment, wherein Δ η
MinFor the default thermal efficiency is adjusted the dead band;
Step B2.5.1.3: if | Δ η (k) | 〉=Δ η
Min, then with Δ O
2s(k) be described oxygen content setting value O
2s(k) optimizing step-length automatic optimal:
O
2s(k)=O
2s(k-1)+ΔO
2s(k)
Wherein: λ is for adjusting coefficient, span be (0,1];
Step B2.5.2:, optimize operating point O with history if process is not in stable state
2i, optBe setting value O
2s(k), guarantee that heating furnace is in the duty of suboptimum;
Step B2.6: the oxygen content of calculating setting value O
2s(k) be sent to described real-time data base, come piece to implement for the Control for Oxygen Content device in the described air and gas system control module;
Step C: the control object of described air and gas system control module comprises the air-fuel ratio when oxygen content, hearth load and fuel increase, the control target is that described oxygen content and combustion chamber draft are controlled in the optimization set point or given range that is drawn jointly by described offline optimization and on-line optimization, and step is as follows:
Step C1: step test obtains the steady-state response time of oxygen content and combustion chamber draft:
Under the steady situation of process, intake is applied the step test signal, record oxygen content change curve, the steady-state response time T of acquisition oxygen content
O2ss
Under the steady situation of process, the blower inlet baffle plate is applied the step test signal, record combustion chamber draft change curve, the steady-state response time T of acquisition combustion chamber draft
Pss
Step C2: the control cycle of setting cigarette wind control system module is T
c, T
c=min (T
O2ss, T
Pss)/40;
Step C3: when fuel increases air-fuel ratio is controlled according to the following steps:
Step C3.1: if fuel recruitment Δ F
f(k
c) last relatively intake adjustment moment k
c-1 fuel quantity F
F0Ratio surpasses predetermined threshold value β, and β ∈ (0,0.2], promptly
The feedforward variation delta F of air intake when then fuel increases
AF(k
c) be
ΔF
aF(k
c)=α·AFR·ΔF
f(k
c)
Wherein: Δ F
f(k
c)=F
f(k
c)-F
F0Be fuel change amount, F
f(k
c) be fuel flow rate,
α is an excess air coefficient, is [1.05,1.15] to the fuel gas span;
AFR is a stoichiometric air/, and to fuel gas, it is calculated as:
Wherein: CO, H
2, C
mH
n, H
2S, O
2Being each constituent content in the fuel gas, is unit with %;
Step C4: according to the following steps oxygen content is carried out Region control:
Step C4.1: if the setting value of oxygen content is O
2s, the regional extent of its permission is:
[O
2s-δ?O
2s+δ]
Wherein: δ is that oxygen content departs from setting value O
2sZone limit, δ ∈ (00.5);
Step C4.2: at Control for Oxygen Content moment k
c,
If O
2(k
c)>O
2HPerhaps O
2(k
c)<O
2L, O wherein
2H=O
2s+ δ, O
L=O
2s-δ, and described Control for Oxygen Content device is not in the stand-by period, then is calculated as follows required intake and changes:
Wherein: Δ F
a(k
c) be the variable quantity of intake,
O
2sBe the setting value of oxygen content,
O
2(k
c) be the measured value of oxygen content;
If O
2L≤ O
2(k
c)≤O
2H, Δ F then
a(k
c)=0;
According to Control for Oxygen Content moment k
cCombustion chamber draft P (k
c) by following principle correction intake variation delta F
a(k
c):
If P (k
c)>P
HS, P
HSBe the upper safety limit of combustion chamber draft, and the variation delta F of intake
a(k
c)>0 then keeps intake constant, makes Δ F
a(k
c)=0,
If P (k
c)<P
LS, P
LSBe the lower safety limit of combustion chamber draft, and the variation delta F of intake
a(k
c)<0 then keeps intake constant, makes Δ F
a(k
c)=0;
Step C4.3: be calculated as follows required intake:
F
a(k
c)=F
a(k
c-1)+ΔF
a(k
c)
Wherein: F
a(k
c) be intake;
After process made feedback regulation, wait for the steady-state response time T of an oxygen content
O2ss
Step C5: according to the following steps combustion chamber draft is carried out Region control:
Step C5.1: if the setting value of combustion chamber draft is P
s, and be in thermic load working region Ω
i, P then
s=P
I, opt, i=1,2 ..., 5, it allows the regional extent of change to be:
[
Ps-σ?Px+σ]
Wherein: σ is the zone limit that combustion chamber draft departs from setting value, δ ∈ (0 10);
Step C5.2: at combustion chamber draft control moment k
c,
If P (k
c)>P
HPerhaps P (k
c)<P
L, P wherein
H=P
s+ σ, P
L=P
s-σ, and described combustion chamber draft controller is not in the stand-by period, and then the deviation of combustion chamber draft is e
P(k
c)=P
s-P (k
c), air-introduced machine inlet baffle variation delta MV then
2(k
c) be
Wherein: K
1Proportionality coefficient for air-introduced machine inlet baffle and negative pressure variation;
After process made feedback regulation, wait for the steady-state response time T of a combustion chamber draft
Pss
Under other situations, Δ MV
2(k
c)=0;
Step C5.3: FEEDFORWARD CONTROL is carried out in the variation of blower variable frequency by following formula:
ΔMV
2F(k
c)=K
2ΔMV
1(k
c)
Wherein: Δ MV
2F(k
c) be the feedforward variation of air-introduced machine inlet baffle,
Δ MV
1(k
c)=MV
1(k
c)-MV
1(k
c-1) variation of exporting for blower variable frequency,
K
2Be the feed-forward coefficients between baffle plate variation and the frequency conversion output.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102170973A CN101881563B (en) | 2010-07-02 | 2010-07-02 | Multi-area intelligent online optimizing control method for thermal efficiency of heating furnace |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102170973A CN101881563B (en) | 2010-07-02 | 2010-07-02 | Multi-area intelligent online optimizing control method for thermal efficiency of heating furnace |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101881563A true CN101881563A (en) | 2010-11-10 |
CN101881563B CN101881563B (en) | 2012-05-23 |
Family
ID=43053635
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010102170973A Active CN101881563B (en) | 2010-07-02 | 2010-07-02 | Multi-area intelligent online optimizing control method for thermal efficiency of heating furnace |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101881563B (en) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102032801A (en) * | 2010-12-22 | 2011-04-27 | 莱芜钢铁股份有限公司 | Air feed control equipment and method for heating furnace |
CN102540883A (en) * | 2010-12-20 | 2012-07-04 | 北京西玛通科技有限公司 | Intelligent optimized energy saving and emission reducing controlling system |
CN102841983A (en) * | 2012-08-10 | 2012-12-26 | 南京南瑞继保电气有限公司 | On-line monitoring method of thermal efficiency of industrial kiln |
CN102865752A (en) * | 2012-08-15 | 2013-01-09 | 北京世纪隆博科技有限责任公司 | Branch temperature balancing and load controlling method for heating furnace |
CN103279652A (en) * | 2013-05-13 | 2013-09-04 | 浙江大学 | Complex chemical process modeling method of DNA genetic algorithm based on swarm behavior |
CN103939940A (en) * | 2014-04-04 | 2014-07-23 | 青岛科技大学 | Boiler combustion optimization and online energy efficiency evaluation system |
CN107977736A (en) * | 2017-11-15 | 2018-05-01 | 浙江中控软件技术有限公司 | The optimization method and system of fuel gas system based on calorific value balance |
CN108036299A (en) * | 2018-01-08 | 2018-05-15 | 中国恩菲工程技术有限公司 | Intelligent waste heat boiler Optimal Control System |
CN108224393A (en) * | 2018-01-08 | 2018-06-29 | 中国恩菲工程技术有限公司 | Intelligent waste heat boiler optimal control method |
CN109269117A (en) * | 2018-08-10 | 2019-01-25 | 中国石油天然气股份有限公司 | Heating furnace thermal efficiency dynamic operation calculation method |
CN114001561A (en) * | 2021-09-22 | 2022-02-01 | 拓正化工工程(上海)有限公司 | Intelligent digital operation and maintenance platform based on Internet of things |
CN114756068A (en) * | 2022-04-11 | 2022-07-15 | 深圳市佳运通电子有限公司 | Method for controlling trusteeship temperature under condition of furnace condition mutation of oil field heating furnace |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4435149A (en) * | 1981-12-07 | 1984-03-06 | Barnes Engineering Company | Method and apparatus for monitoring the burning efficiency of a furnace |
CN101256418A (en) * | 2008-03-28 | 2008-09-03 | 清华大学 | Combination control method for exit temperature of heating furnace |
-
2010
- 2010-07-02 CN CN2010102170973A patent/CN101881563B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4435149A (en) * | 1981-12-07 | 1984-03-06 | Barnes Engineering Company | Method and apparatus for monitoring the burning efficiency of a furnace |
CN101256418A (en) * | 2008-03-28 | 2008-09-03 | 清华大学 | Combination control method for exit temperature of heating furnace |
Non-Patent Citations (3)
Title |
---|
《仪器仪表学报》 20090630 吕文祥等 加热炉热效率的智能优化控制研究 第30卷, 第6期 2 * |
《化工自动化及仪表》 19931231 赵育才 加热炉多变量热效率优化控制 第20卷, 第6期 2 * |
《燃烧科学与技术》 19951231 唐林等 工业加热炉最佳燃烧状况自动寻优控制系统的研究 第1卷, 第2期 2 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102540883A (en) * | 2010-12-20 | 2012-07-04 | 北京西玛通科技有限公司 | Intelligent optimized energy saving and emission reducing controlling system |
CN102032801A (en) * | 2010-12-22 | 2011-04-27 | 莱芜钢铁股份有限公司 | Air feed control equipment and method for heating furnace |
CN102841983B (en) * | 2012-08-10 | 2015-11-18 | 南京南瑞继保电气有限公司 | A kind of thermal efficiency of industrial kiln on-line monitoring method |
CN102841983A (en) * | 2012-08-10 | 2012-12-26 | 南京南瑞继保电气有限公司 | On-line monitoring method of thermal efficiency of industrial kiln |
CN102865752A (en) * | 2012-08-15 | 2013-01-09 | 北京世纪隆博科技有限责任公司 | Branch temperature balancing and load controlling method for heating furnace |
CN102865752B (en) * | 2012-08-15 | 2014-07-16 | 北京世纪隆博科技有限责任公司 | Branch temperature balancing and load controlling method for heating furnace |
CN103279652B (en) * | 2013-05-13 | 2016-03-30 | 浙江大学 | A kind of complicated chemical procedure modelling approach of the DNA genetic algorithm based on bee colony behavior |
CN103279652A (en) * | 2013-05-13 | 2013-09-04 | 浙江大学 | Complex chemical process modeling method of DNA genetic algorithm based on swarm behavior |
CN103939940A (en) * | 2014-04-04 | 2014-07-23 | 青岛科技大学 | Boiler combustion optimization and online energy efficiency evaluation system |
CN107977736A (en) * | 2017-11-15 | 2018-05-01 | 浙江中控软件技术有限公司 | The optimization method and system of fuel gas system based on calorific value balance |
CN107977736B (en) * | 2017-11-15 | 2021-04-09 | 浙江中控软件技术有限公司 | Fuel gas system optimization method and system based on heat value balance |
CN108036299A (en) * | 2018-01-08 | 2018-05-15 | 中国恩菲工程技术有限公司 | Intelligent waste heat boiler Optimal Control System |
CN108224393A (en) * | 2018-01-08 | 2018-06-29 | 中国恩菲工程技术有限公司 | Intelligent waste heat boiler optimal control method |
CN108036299B (en) * | 2018-01-08 | 2019-11-08 | 中国恩菲工程技术有限公司 | Intelligent waste heat boiler Optimal Control System |
CN109269117A (en) * | 2018-08-10 | 2019-01-25 | 中国石油天然气股份有限公司 | Heating furnace thermal efficiency dynamic operation calculation method |
CN109269117B (en) * | 2018-08-10 | 2020-10-02 | 中国石油天然气股份有限公司 | Method for determining operating state of heating furnace |
CN114001561A (en) * | 2021-09-22 | 2022-02-01 | 拓正化工工程(上海)有限公司 | Intelligent digital operation and maintenance platform based on Internet of things |
CN114756068A (en) * | 2022-04-11 | 2022-07-15 | 深圳市佳运通电子有限公司 | Method for controlling trusteeship temperature under condition of furnace condition mutation of oil field heating furnace |
Also Published As
Publication number | Publication date |
---|---|
CN101881563B (en) | 2012-05-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101881563B (en) | Multi-area intelligent online optimizing control method for thermal efficiency of heating furnace | |
CN201737965U (en) | Automatic optimizing combustion intelligent control system of hot blast stove | |
CN101892338B (en) | Constant wind temperature control system for hot blast stove | |
CN101556038B (en) | Optimization control system for stable operation and economical combustion of circulating fluidized-bed boiler | |
CN101749730B (en) | Coal burning heat medium furnace high-precision temperature regulation and combustion control method | |
WO2020015694A1 (en) | Combustion air volume control system and method for biomass boiler | |
CN102453792B (en) | Method and equipment for controlling furnace pressure of continuous annealing furnace | |
CN103423750B (en) | A kind of method for controlling combustion of incinerator | |
CN112795716B (en) | Efficient and practical hot blast stove burning control method | |
CN106774521A (en) | The temprature control method of heat accumulation type radiant tube fast pyrogenation stove | |
CN104850151A (en) | Temperature control method for airflow type cut tobacco dryer combustion chamber | |
WO2022252414A1 (en) | Coal-air synchronous dynamic coordinated control method for coal-fired unit | |
CN105605609A (en) | Thermal power plant boiler combustion oxygen amount optimizing method | |
CN101749731A (en) | Automatic optimal combustion intelligent control system for hot blast stove | |
CN103593578A (en) | Flue suction force feedback setting method in coke oven heating combustion process | |
CN114721253A (en) | Heating furnace temperature fractional order PID control system and method based on artificial bee colony algorithm | |
CN109539359A (en) | The phase transformation electric heat-storage heating system and method for divided working status PID+ adaptive feedforward compensation | |
CN203052713U (en) | Efficient double thermal storage type heating furnace pulse combustion device | |
CN103939880B (en) | A kind of self-adapting intelligent control method of gasification of biomass thermal power transfer system | |
CN202512382U (en) | Chain boiler burning rolling self-optimization - proportion integration differentiation (PID) compound control system | |
CN115371055A (en) | Automatic control method and system for waste incineration | |
CN113741200A (en) | Intelligent optimization calcination control system for lime sleeve kiln | |
CN109579044B (en) | Air-fuel ratio dynamic feedforward combustion control method for walking beam furnace | |
CN109055640B (en) | Method and device for accurately controlling temperature of furnace burning flue of blast furnace hot blast stove | |
CN111780155A (en) | Control system and control method for high-efficiency combustion hot blast stove |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20180913 Address after: 519085 block A, Tsinghua Science and Technology Park (Zhuhai), No. 101, University Road, Tangjia Bay, Zhuhai high tech Zone, Guangdong Patentee after: Zhuhai Czech Pioneer Technology Co., Ltd. Address before: 100084 mailbox 100084-82, Beijing City Patentee before: Tsinghua University |
|
TR01 | Transfer of patent right |