CN103146905A - Heating furnace temperature decision-making method based on billet optimizing heating curve - Google Patents

Abstract

The invention provides a heating furnace temperature decision-making method based on a billet optimizing heating curve. The method comprises a first step of initialization of a calculation parameter, a second step of calculation segment division of a heating furnace, a third step of position tracking and temperature tracking of all billets in the heating furnace, a fourth step of judging whether the position tracking and the temperature tracking of all the billets in the heating furnace are finished, returning to the third step if the positioning tracking and the temperature tracking are not finished and continuing if the positioning tracking and the temperature tracking are finished, and a fifth step of judging whether a furnace temperature decision-making cycle arrives, returning to the third step if the furnace temperature decision-making cycle does not arrive and conducting the furnace temperature decision-making if the furnace temperature decision-making cycle arrives. The method is simple in correcting algorithm, fast in calculating speed, high in efficiency and beneficial for achieving online control of the heating furnace, saves energy and reduces consumption.

Description

A kind of furnace temperature of heating furnace decision-making technique of optimizing heating curve based on steel billet

Technical field

The present invention relates to process furnace mathematical Model Development and process furnace Optimized-control Technique field, relate in particular to a kind of furnace temperature of heating furnace decision-making technique of optimizing heating curve based on steel billet.

Background technology

Process furnace is the visual plant during hot rolling is produced, also that one of larger equipment of line power consumption is produced in hot rolling, along with power supply shortage and energy demand increase, the contradiction of excessive resources and energy consumption amount and limited resources environmental carrying capacity outstanding, the quick upgrading of hot-rolled product and tooling in addition, high added value steel plate make process furnace face severe challenge energy-saving and cost-reducing and raising self control level to improving constantly that the process furnace heating process requires.

Various equivalent modifications, carried out a large amount of research work for the problems referred to above, especially aspect process furnace optimal control mathematical Model Development, most typically mainly contain two kinds: the one,, the furnace temperature of heating furnace set(ting)value Optimized model that based on data excavates and statistical theory is developed, it mainly comprises data gathering and two modules of data mining analysis.Data processing is mainly completed real-time data acquisition and is obtained the current warm system of process furnace, then by the data mining analysis module, the production data that gathers is processed into analyzing samples and according to certain rale store in database, obtain the optimal furnace temperature set value of process furnace current state during production from database; The 2nd,, based on furnace temperature of heating furnace mathematical model of decision and the Steel In Reheating Furnace base temperature prediction mathematical model of thermal technology's Specialized Theory exploitation.

The former has bypassed thermal technology's Specialized Theory, and the good and bad of sample directly affects furnace temperature decision-making and the forecast of steel temperature, thereby affects the Heating Furnace Control effect, and it is used for instructing actual production obviously can't rationally realize process furnace optimum control and energy-saving and cost-reducing aim; The latter more pays attention to furnace temperature of heating furnace decision-making and the forecast of steel temperature, and it mainly completes the calculating of process furnace optimal furnace temperature set value by mathematical model based on thermal technology's Specialized Theory.But its furnace temperature set(ting)value generates mainly to reach the steel billet tapping temperature as target, and the steel billet tapping temperature of calculating by mathematical model and the deviation of target tapping temperature constantly revise the furnace temperature set(ting)value, improves and consider intensification situation in the middle of process furnace energy consumption and steel billet.

Summary of the invention

The present invention is directed to above-mentioned technical problem, heat-transfer mechanism in the process furnace burner hearth, a kind of furnace temperature of heating furnace decision-making technique of optimizing heating curve based on steel billet is proposed, optimizing heating curve by virtual heating take steel billet revises a corresponding furnace temperature of steel billet current position by foundation, then the corresponding furnace temperature of all steel billets in each control section of process furnace is weighted on average, thereby calculates the decision-making furnace temperature that obtains each control section of process furnace.The present invention has considered the impact of all steel billets in process furnace, is beneficial to and realizes the process furnace optimal control.A corresponding furnace temperature of steel billet is revised, be based on the deviation of the optimization heating curve target temperature of the accounting temperature of this steel billet current position model and current position, and correction algorithm is simple, and computing velocity is fast, efficient is high, is beneficial to realize the process furnace on-line Control.Simultaneously, will be heated to take steel billet steel billet that the required minimum fuel consumption amount of target temperature set up as objective function and optimize heating curve and be introduced in the furnace temperature of heating furnace decision making algorithm in process furnace, and be beneficial to and realize Energy Saving of Heating Furnaces.

The present invention solve the technical problem and mainly takes following technical scheme: a kind of furnace temperature of heating furnace decision-making technique of optimizing heating curve based on steel billet comprises the following steps:

S1) calculating parameter initialize;

S2) process furnace is carried out the division of compute segment;

S3) all steel billets in process furnace are carried out location tracking and temperature tracking;

S4) judge whether to complete location tracking and the temperature tracking of all steel billets in process furnace, do not return to step S3 if complete; Enter next step if complete;

S5) judge whether to arrive furnace temperature decision-making period, if do not arrive furnace temperature decision-making period, return to step S3; As arrive furnace temperature decision-making period, carry out the furnace temperature decision-making.

Described method, the initialized parameter of step S1 comprises: the heating furnace structure parameter, add Steel In Reheating Furnace base position and temperature field, each control section furnace temperature of process furnace, the initialize of steel billet counter, Steel In Reheating Furnace base sum, timer parameter.

Described method, the method that step S2 divides comprises: according to the number of thermopair in process furnace and the layout in process furnace thereof, each control section of process furnace is divided into several compute segment.

Described method, the method that step S3 temperature is followed the tracks of comprises:

Find the solution steel billet heat conduction governing equation, calculate and obtain the inside steel billet temperature distribution, its heat conduction governing equation is as follows:

$\mathrm{\ρ}\left(T\right)\·c\left(T\right)\·\frac{\∂T(y,\mathrm{\τ})}{\∂\mathrm{\τ}}=\frac{\∂}{\∂y}[\mathrm{\λ}\left(T\right)\·\frac{\∂T(y,\mathrm{\τ})}{\∂y}]$

In formula, ρ (T) is steel billet density, and c (T) is steel billet specific heat, and λ (T) is the steel billet thermal conductivity, and τ is the time, and T=T (y, τ) is steel billet temperature, and y is the steel billet thickness coordinate;

Its final condition is:

$q_u=\mathrm{\λ}\left(T\right)\·\frac{\∂T(y,\mathrm{\τ})}{\∂y}{|}_{y=d}$

$q_b=\mathrm{\λ}\left(T\right)\·\frac{\∂T(y,\mathrm{\τ})}{\∂y}{|}_{y=0}$

Its starting condition is:

T(y,τ)| _τ＝0＝f(y)0≤y≤d

In formula, q _uBe steel billet upper surface heat flow density, q _bBe steel billet lower surface heat flow density, f (y) is steel billet initial temperature field, and d is steel billet thickness.

Described method, the method for step S5 furnace temperature decision-making comprises:

S51) steel billet counter j is initially 1;

S52) the timing register τ of j piece steel billet in process furnace _jBe initialized as zero, generate the optimization heating curve of j piece steel billet, the inferior counter i of correction of the corresponding furnace temperature of steel billet institute is initialized as 1;

S53) calculate the deviation △ T of j piece steel billet Current Temperatures and its optimization heating curve target temperature _j, according to this deviation, the corresponding furnace temperature of steel billet j is revised, obtain revised furnace temperature Tf _{I, j, k}

S54) with furnace temperature Tf _{I, j, k}Steel billet j is carried out virtual heating, the temperature after the i time virtual heating of calculating steel billet j and the deviation △ T of its optimization heating curve target temperature _{I, j}

S55) judgement △ T _{I, j}Whether less than or equal to allowable value ξ, if greater than allowable value ξ, with △ T _{I, j}Assignment is to △ T _jGetting back to step S53 continues to calculate; If △ is T _{I, j}Less than or equal to allowable value ξ, with furnace temperature Tf _{I, j, k}As the corresponding decision-making furnace temperature of steel billet j;

S56) judge whether steel billet counter j equals Steel In Reheating Furnace base overall number N, if j is not equal to N, steel billet counter j adds 1, then get back to step S52 and continue to calculate; If steel billet counter j equals Steel In Reheating Furnace base overall number N, continue following steps;

S57) according to the furnace temperature decision value Tf of every block of steel billet in process furnace _{I, j, k}Calculate the furnace temperature decision value Tf that obtains each control section of process furnace _k

Described method, step S53 obtains revised furnace temperature Tf _{I, j, k}Method obtain by following formula:

${\mathrm{Tf}}_{i,j,k}={\mathrm{Tf}}_{i,j,k}\·(1-c_i\·\frac{\mathrm{\Δ}{T}_{i,j}}{T_{\mathrm{Opt},j}})$

In formula, the Tf on equation the right _{I, j, k}The corresponding furnace temperature of steel billet j when being the i time virtual heating, the Tf on the equation left side _{I, j, k}Be the corresponding furnace temperature of steel billet j of revising after the i time virtual heating, it is the corresponding furnace temperature of steel billet j during as the i+1 time virtual heating, c _iBe correction factor, T _{Opt, j}Be the optimization heating curve target temperature of steel billet j, Δ T _i,jBe temperature and its optimization heating curve target temperature T of steel billet j after the i time virtual heating _{Opt, j}Deviation.

Described method, the method for the virtual heating of step S54 comprises:

S541) with the initial temperature field T of steel billet j _{0, j}And the corresponding furnace temperature Tf of steel billet j _{I, j, k}Steel billet is carried out location tracking and temperature tracking;

S542) temperature of steel billet j and the deviation △ T of its optimization heating curve target temperature after accounting temperature is followed the tracks of _{I, j}

S543) judgement △ T _{I, j}Whether (value that ξ generally gets between-5 ℃ to 5 ℃ gets final product less than or equal to allowable value ξ; In fact ξ is equivalent to a convergence, and what ξ can get when computational accuracy is had relatively high expectations is smaller, but program runtime is longer like this; Otherwise ξ can suitably relax, to improve the program computing velocity):

S5431) if △ is T _{I, j}Less than or equal to allowable value ξ, finish to the virtual heating of steel billet j and with Tf _{I, j, k}As the corresponding decision-making furnace temperature of steel billet j;

S5432) if greater than allowable value ξ, according to △ T _{I, j}To Tf _{I, j, k}Revise, obtain revising the corresponding furnace temperature Tf of rear steel billet j _{I+1j, k}, get back to step S541 still with T _{0, j}As the initial temperature field of steel billet j, but with the corresponding furnace temperature Tf of revised steel billet j _{I+1, j, k}Continue to calculate.

Described method, step S57 obtains furnace temperature decision value Tf _kMethod comprise:

Be calculated as follows the furnace temperature decision value Tf that obtains process furnace k control section _k:

${\mathrm{Tf}}_k=\underset{j=1}{\overset{N_k}{\mathrm{\Σ}}}({\mathrm{\ω}}_j\·{\mathrm{Tf}}_{i,j,k})$

In formula, N _kBe the number of steel billet in k control section of process furnace, ω _jBe the furnace temperature decision value Tf of j piece steel billet _{I, j, k}Corresponding weights (weights ω _jCan get the Arbitrary Digit between 0 to 1, comprise 0 or 1, but should satisfy

Its weights of steel billet in general control section mid-way are got high, and both sides are got low, and be parabolic shape and be advisable, such as: there are 5 blocks of steel billets to export from the control section entrance to control section in control section and are followed successively by B1, B2, B3, B4, B5, desirable ω ₁=0.15, ω ₂=0.2, ω ₃=0.3, ω ₄=0.2, ω ₅=0.15), Tf _{I, j, k}Be i the furnace temperature decision value of j piece steel billet in k control section of process furnace.

The present invention compared with prior art has following main beneficial effect:

1. its furnace temperature decision-making is optimized heating curve based on steel billet, it is to be heated to the required minimum fuel consumption amount of target temperature to be objective function take steel billet in process furnace that steel billet is optimized heating curve, the motion of simulation calculation steel billet in process furnace and the heating steel billet curve of heat temperature raising process acquisition, be introduced into to the furnace temperature of heating furnace decision making algorithm, be beneficial to and realize Energy Saving of Heating Furnaces.

2. optimizing heating curve by virtual heating take steel billet revises a corresponding furnace temperature of steel billet by foundation, then to all steel billets in process furnace corresponding furnace temperature be weighted the furnace temperature that average computation obtains each control section of process furnace, consider the impact of all steel billets in process furnace, be beneficial to and realize the process furnace optimal control.

3. a kind of furnace temperature of heating furnace decision-making technique of optimizing heating curve based on steel billet provided by the present invention, its inside there is no complicated algorithm and more iteration, and computing velocity is fast, efficient is high, is beneficial to realize the process furnace on-line Control.

Description of drawings

Fig. 1 is general flow chart of the present invention.

Fig. 2 is the furnace temperature decision flow diagram.

Embodiment

A kind of furnace temperature of heating furnace decision-making technique of optimizing heating curve based on steel billet provided by the present invention, it comprises: the calculating parameter initialization step, process furnace is carried out compute segment partiting step, steel billet location tracking and temperature tracking step, furnace temperature decision-making judgement and performing step.

Described calculating parameter initialization step mainly comprises timing register τ and steel billet counter j initialize, Steel In Reheating Furnace base position and temperature field initialize, each control section furnace temperature initialize of process furnace.

Described process furnace is carried out the compute segment partiting step, main according to the number of thermopair in process furnace and the layout in process furnace thereof, each control section of process furnace is divided into several compute segment.

Described steel billet location tracking and temperature tracking step, its steel billet location tracking are mainly to obtain Steel In Reheating Furnace base current location information; It is mainly to obtain the inside steel billet temperature field by finding the solution steel billet heat conduction governing equation, calculating that its steel billet temperature is followed the tracks of.

Described furnace temperature decision-making judgement and performing step, its furnace temperature decision-making determining step are mainly to judge whether to arrive furnace temperature decision-making period according to timing register timing situation; Its furnace temperature decision-making performing step mainly refers to according to when front furnace internal temperature field, each control section furnace temperature set(ting)value of process furnace being calculated.

Its furnace temperature decision-making performing step mainly comprises: the optimization heating curve of steel billet generates, the virtual heating of correction, steel billet of the corresponding furnace temperature of steel billet institute, each control section furnace temperature set(ting)value generation of process furnace.

The optimization heating curve of described steel billet generates, and generates the optimization heating curve of Steel In Reheating Furnace base.Described steel billet is optimized heating curve, it comprises the position array POS[M-1 of steel billet in process furnace] and the steel billet node temperature array TEMP[M-1 corresponding with each position, N-1], wherein M is the interstitial content of steel billet thickness direction for position number, the N that divides by the heating-furnace length direction.

A described steel billet the correction of corresponding furnace temperature, calculate the deviation △ T of j piece steel billet Current Temperatures and its optimization heating curve target temperature in process furnace _j, according to this deviation, the corresponding furnace temperature of steel billet j is revised, obtain revised furnace temperature Tf _{I, j, k}

The virtual heating of described steel billet is with furnace temperature Tf _{I, j, k}Steel billet j is carried out virtual heating, the steel billet j after virtual heating is carried out location tracking and temperature tracking, the temperature after the virtual heating of calculating steel billet j and the deviation △ T of its optimization heating curve target temperature _i,j

It is characterized in that, the virtual heating of steel billet mainly comprises the steps:

(1) with the initial temperature field T of steel billet j _{0, j}And the corresponding furnace temperature Tf of steel billet j _{I, j, k}Steel billet is carried out location tracking and temperature tracking;

(2) temperature of steel billet j and the deviation △ T of its optimization heating curve target temperature after accounting temperature is followed the tracks of _i,j

(3) judgement △ T _i,jWhether less than or equal to allowable value ξ:

If △ Ti A., j finishes to the virtual heating of steel billet j and with Tf less than or equal to allowable value ξ _{I, j, k}As the corresponding decision-making furnace temperature of steel billet j;

If B. greater than allowable value ξ, according to △ Ti, j is to Tf _{I, j, k}Revise, obtain revising the corresponding furnace temperature Tf of rear steel billet j _{I+1j, k}, getting back to step (1) still with T0, j is as the initial temperature field of steel billet j, but with the corresponding furnace temperature Tf of revised steel billet j _{I+1, j, k}Continue to calculate.

It is characterized in that, when steel billet j is carried out virtual heating for the first time with T0, j carries out temperature as the initial temperature field of steel billet j to be followed the tracks of and location tracking, if steel billet j is carried out repeatedly virtual heating, each virtual heating is all with T0, j carries out temperature as the initial temperature field of steel billet j to be followed the tracks of and location tracking, but not with the temperature of steel billet j after the virtual heating last time initial temperature as this virtual heating.

It is characterized in that, by following formula correction Tf _i,j:

${\mathrm{Tf}}_{i,j,k}={\mathrm{Tf}}_{i,j,k}\·(1-c_i\·\frac{\mathrm{\Δ}{T}_{i,j}}{T_{\mathrm{Opt},j}})$

In formula, the Tf on equation the right _{I, j, k}The corresponding furnace temperature of steel billet j when being the i time virtual heating; The Tf on the equation left side _{I, j, k}Be the corresponding furnace temperature of steel billet j of revising after the i time virtual heating, it is the corresponding furnace temperature of steel billet j during as the i+1 time virtual heating; c _iBe correction factor, generally get the inverse of the virtual heating times i of steel billet j, namely

T _{Opt, j}Optimization heating curve target temperature for steel billet j; Δ T _i,jBe temperature and its optimization heating curve target temperature T of steel billet j after the i time virtual heating _{Opt, j}Deviation.

Each control section furnace temperature set(ting)value of described process furnace generates, according to the furnace temperature decision value Tf of every block of steel billet in process furnace _{I, j, k}Calculate the furnace temperature decision value Tf that obtains each control section of process furnace _k

It is characterized in that Tf _kBe calculated as follows:

${\mathrm{Tf}}_k=\underset{j=1}{\overset{N_k}{\mathrm{\Σ}}}({\mathrm{\ω}}_j\·{\mathrm{Tf}}_{i,j,k})$

In formula, Tf _kDecision-making furnace temperature for k control section of process furnace; N _kNumber for steel billet in k control section of process furnace; ω _jBe its furnace temperature decision value of j piece steel billet Tf _{I, j, k}Corresponding weights; Tf _{I, j, k}I the furnace temperature decision value of j piece steel billet in k control section of process furnace.

A kind of furnace temperature of heating furnace decision-making technique of optimizing heating curve based on steel billet, its step:

1. calculating parameter initialize, the heating furnace structure parameter initialization, add Steel In Reheating Furnace base position and temperature field initialize, each control section furnace temperature initialize of process furnace, steel billet counter initialize (j=1), Steel In Reheating Furnace base sum is initialized as N, timer parameter initialize (τ=0).

2. process furnace is carried out the division of compute segment, according to the number of thermopair in process furnace and the layout in process furnace thereof, each control section of process furnace is divided into several compute segment.

3. all steel billets in process furnace are carried out location tracking and temperature tracking.

4. judge whether to complete location tracking and the temperature tracking of all steel billets in process furnace according to steel billet counter j.Do not follow the tracks of if complete location tracking and the temperature of all steel billets in process furnace, steel billet counter j adds 1, gets back to step 3; Follow the tracks of step below continuing if complete location tracking and the temperature of all steel billets in process furnace.

5. judge whether to arrive furnace temperature decision-making period, if do not arrive furnace temperature decision-making period, timing register continues timing, i.e. τ=τ+Δ τ, and wherein τ is that timing register time, Δ τ are time step, then gets back to step 3; As arrive furnace temperature decision-making period, carry out according to the following steps the furnace temperature decision-making:

A. steel billet counter j is initially 1;

B. the timing register τ of j piece steel billet in process furnace _jBe initialized as zero, generate the optimization heating curve of j piece steel billet, the inferior counter i of correction of the corresponding furnace temperature of steel billet institute is initialized as 1;

C. calculate the deviation △ Tj of j piece steel billet Current Temperatures and its optimization heating curve target temperature, according to this deviation, the corresponding furnace temperature of steel billet j is revised, obtain revised furnace temperature Tf _{I, j, k}

D. with furnace temperature Tf _{I, j, k}Steel billet j is carried out virtual heating, the temperature after the i time virtual heating of calculating steel billet j and the deviation △ Ti of its optimization heating curve target temperature, j;

E. judge △ Ti, whether j is less than or equal to allowable value ξ, if greater than allowable value ξ, with △ Ti, the j assignment is got back to step B to △ Tj and continued to calculate; If △ is Ti, j is less than or equal to allowable value ξ, with furnace temperature Tf _{I, j, k}As the corresponding decision-making furnace temperature of steel billet j;

F. judge whether steel billet counter j equals Steel In Reheating Furnace base overall number N, if j is not equal to N, steel billet counter j adds 1, then get back to step a and continue to calculate; If steel billet counter j equals Steel In Reheating Furnace base overall number N, continue following steps;

G. according to the furnace temperature decision value Tf of every block of steel billet in process furnace _{I, j, k}Calculate the furnace temperature decision value Tf that obtains each control section of process furnace _k

Describe embodiments of the present invention in detail below in conjunction with accompanying drawing 1, accompanying drawing 2.

1. calculating parameter initialize (step 1) mainly comprises: the heating furnace structure parameter initialization, add Steel In Reheating Furnace base position and temperature field initialize, each control section furnace temperature initialize of process furnace, steel billet counter initialize (j=1), Steel In Reheating Furnace base sum is initialized as N, timer parameter initialize (τ=0).

2. process furnace is carried out the division (step 2) of compute segment, according to the number of thermopair in process furnace and the layout in process furnace thereof, each control section of process furnace is divided into several compute segment, total compute segment number that process furnace is divided is designated as M.

3. j piece steel billet in process furnace is carried out location tracking and temperature tracking (step 3).At first j piece steel billet in process furnace is carried out location tracking (step 4), namely read and obtain the positional information of steel billet in process furnace; J piece steel billet in process furnace is carried out temperature follow the tracks of (step 5), namely find the solution steel billet heat conduction governing equation, calculate and obtain the inside steel billet temperature distribution, its heat conduction governing equation is as follows:

$\mathrm{\ρ}\left(T\right)\·c\left(T\right)\·\frac{\∂T(y,\mathrm{\τ})}{\∂\mathrm{\τ}}=\frac{\∂}{\∂y}[\mathrm{\λ}\left(T\right)\·\frac{\∂T(y,\mathrm{\τ})}{\∂y}]$

In formula, ρ is steel billet density; C is steel billet specific heat; λ is the steel billet thermal conductivity.

Its final condition is:

$q_u=\mathrm{\λ}\left(T\right)\·\frac{\∂T(y,\mathrm{\τ})}{\∂y}{|}_{y=d}$

$q_b=\mathrm{\λ}\left(T\right)\·\frac{\∂T(y,\mathrm{\τ})}{\∂y}{|}_{y=0}$

Its starting condition is:

T(y,τ)| _τ＝0＝f(y)0≤y≤d

In formula, q _uBe steel billet upper surface heat flow density; q _bBe steel billet lower surface heat flow density; F (y) is steel billet initial temperature field; D is steel billet thickness.

4. judge whether to complete location tracking and the temperature tracking (step 6) of all steel billets in process furnace according to steel billet counter j.Do not follow the tracks of if complete location tracking and the temperature of all steel billets in process furnace, steel billet counter j adds 1 (step 7), gets back to step 3; Follow the tracks of step below continuing if complete location tracking and the temperature of all steel billets in process furnace.

5. according to timing register timing situation, judge whether to arrive furnace temperature decision-making period (step 8), if do not arrive furnace temperature decision-making period, timing register continues timing, be τ=τ+Δ τ (step 9), wherein τ is that timing register time, Δ τ are time step, then gets back to step 3; As arrive furnace temperature decision-making period, carry out following steps and carry out the furnace temperature decision-making:

A. counter initialize (step 11), steel billet counter j is initially 1;

B. the timing register τ of j piece steel billet in process furnace _jBe initialized as zero, generate the optimization heating curve of j piece steel billet, the inferior counter i of correction of the corresponding furnace temperature of steel billet institute is initialized as 1 (step 12);

C. calculate the deviation △ T of j piece steel billet Current Temperatures and its optimization heating curve target temperature _j, the corresponding furnace temperature of steel billet j is revised by following formula according to this deviation, obtain revised furnace temperature Tf _{I, j, k}(step 13):

${\mathrm{Tf}}_{i,j,k}={\mathrm{Tf}}_{i,j,k}\·(1-c_i\·\frac{\mathrm{\Δ}{T}_{i,j}}{T_{\mathrm{Opt},j}})$

In formula, the Tf on equation the right _{I, j, k}The corresponding furnace temperature of steel billet j when being the i time virtual heating; The Tf on the equation left side _{I, j, k}Be the corresponding furnace temperature of steel billet j of revising after the i time virtual heating, it is the corresponding furnace temperature of steel billet j during as the i+1 time virtual heating; c _iBe correction factor, generally get the inverse of the virtual heating times i of steel billet j, namely T _{Opt, j}Optimization heating curve target temperature for steel billet j; Δ T _i,jBe temperature and its optimization heating curve target temperature T of steel billet j after the i time virtual heating _{Opt, j}Deviation.

D. with furnace temperature Tf _{I, j, k}Steel billet j is carried out virtual heating (step 14), the temperature after the i time virtual heating of calculating steel billet j and the deviation △ Ti of its optimization heating curve target temperature, j (step 15);

E. judge △ Ti, whether j is less than or equal to allowable value ξ (step 16), if greater than allowable value ξ, with △ Ti, the j assignment is to △ Tj, and virtual heating times counter i adds 1 (step 17), gets back to step c and continues to calculate; If △ is Ti, j is less than or equal to allowable value ξ, with furnace temperature Tf _{I, j, k}, as the corresponding decision-making furnace temperature of steel billet j (step 18) (Tf _{I, j, k}Be corresponding k the furnace temperature of steel billet j, it obtains through after i virtual heating);

F. judge whether steel billet counter j equals Steel In Reheating Furnace base overall number N (step 19), if j is not equal to N, steel billet counter j adds 1, then get back to step a and continue to calculate; If steel billet counter j equals Steel In Reheating Furnace base overall number N, continue following steps;

G. according to (the furnace temperature decision value Tf of j=1 ~ N) of every block of steel billet in process furnace _{I, j, k}Be calculated as follows the furnace temperature decision value Tf that obtains each control section of process furnace _k:

${\mathrm{Tf}}_k=\underset{j=1}{\overset{N_k}{\mathrm{\Σ}}}({\mathrm{\ω}}_j\·{\mathrm{Tf}}_{i,j,k})$

In formula, Tf _kDecision-making furnace temperature for k control section of process furnace; N _kNumber for steel billet in k control section of process furnace; ω _jBe its furnace temperature decision value of j piece steel billet Tf _{I, j, k}Corresponding weights; Tf _{I, j, k}I the furnace temperature decision value of j piece steel billet in k control section of process furnace.

Embodiment

Certain process furnace, it comprises: heat one section, heat two sections, heating three sections, soaking zone, the total furnace superintendent of process furnace is 50.5m.By method provided by the present invention, heat one section be divided into two sections of 2 compute segment, heating be divided into 3 compute segment, heating divided in three sections minutes 3 compute segment, soaking zones be divided into 2 compute segment with each stove section itself as 3 compute segment, i.e. compute segment number M=10.

The process furnace initial furnace temperature of each stove section such as table 1.

The initial furnace temperature of each stove section of table 1

The stove section	Heat one section	Heat two sections	Heat three sections	Soaking zone
					Furnace temperature (℃)	1017	1043	1136	1193

Steel billet initial parameter information such as table 2 (in table 2, the temperature of steel billet and positional information are that 3 pairs of Steel In Reheating Furnace bases of performing step carry out location tracking and temperature is followed the tracks of acquisition).

Table 2 steel billet initial information

The every corresponding furnace temperature decision value of steel billet Tf _{I, j, k}Weights ω _j=1 (j=1,2,3).By method provided by the present invention, namely step 1 ~ 5, carry out the furnace temperature decision-making, the furnace temperature result of decision such as table 3.

Each stove section furnace temperature decision value of table 3

The stove section	Heat one section	Heat two sections	Heat three sections	Soaking zone
					The furnace temperature decision value (℃)	873	1111	1125	1202

[0123]By this decision-making furnace temperature, the Steel In Reheating Furnace base is heated and to obtain good heats.

Claims (8)

1. a furnace temperature of heating furnace decision-making technique of optimizing heating curve based on steel billet, is characterized in that, comprises the following steps:

S1) calculating parameter initialize;

S2) process furnace is carried out the division of compute segment;

S3) all steel billets in process furnace are carried out location tracking and temperature tracking;

S4) judge whether to complete location tracking and the temperature tracking of all steel billets in process furnace, do not return to step S3 if complete; Enter next step if complete;

S5) judge whether to arrive furnace temperature decision-making period, if do not arrive furnace temperature decision-making period, return to step S3; As arrive furnace temperature decision-making period, carry out the furnace temperature decision-making.

2. method according to claim 1, it is characterized in that, the initialized parameter of step S1 comprises: the heating furnace structure parameter, add Steel In Reheating Furnace base position and temperature field, each control section furnace temperature of process furnace, the initialize of steel billet counter, Steel In Reheating Furnace base sum, timer parameter.

3. method according to claim 1, is characterized in that, the method that step S2 divides comprises: according to the number of thermopair in process furnace and the layout in process furnace thereof, each control section of process furnace is divided into several compute segment.

4. method according to claim 1, is characterized in that, the method that step S3 temperature is followed the tracks of comprises:

Find the solution steel billet heat conduction governing equation, calculate and obtain the inside steel billet temperature distribution, its heat conduction governing equation is as follows:

$\mathrm{\ρ}\left(T\right)\·c\left(T\right)\·\frac{\∂T(y,\mathrm{\τ})}{\∂\mathrm{\τ}}=\frac{\∂}{\∂y}[\mathrm{\λ}\left(T\right)\·\frac{\∂T(y,\mathrm{\τ})}{\∂y}]$

In formula, ρ (T) is steel billet density, and c (T) is steel billet specific heat, and λ (T) is the steel billet thermal conductivity, and τ is the time, and T=T (y, τ) is steel billet temperature, and y is the steel billet thickness coordinate;

Its final condition is:

$q_u=\mathrm{\λ}\left(T\right)\·\frac{\∂T(y,\mathrm{\τ})}{\∂y}{|}_{y=d}$

$q_b=\mathrm{\λ}\left(T\right)\·\frac{\∂T(y,\mathrm{\τ})}{\∂y}{|}_{y=0}$

Its starting condition is:

T(y,τ)| _τ＝0＝f(y)0≤y≤d

In formula, q _uBe steel billet upper surface heat flow density, q _bBe steel billet lower surface heat flow density, f (y) is steel billet initial temperature field, and d is steel billet thickness.

5. method according to claim 1, is characterized in that, the method for step S5 furnace temperature decision-making comprises:

S51) steel billet counter j is initially 1;

S52) the timing register τ of j piece steel billet in process furnace _jBe initialized as zero, generate the optimization heating curve of j piece steel billet, the inferior counter i of correction of the corresponding furnace temperature of steel billet institute is initialized as 1;

S53) calculate the deviation △ T of j piece steel billet Current Temperatures and its optimization heating curve target temperature _j, according to this deviation, the corresponding furnace temperature of steel billet j is revised, obtain revised furnace temperature Tf _{I, j, k}

S54) with furnace temperature Tf _{I, j, k}Steel billet j is carried out virtual heating, the temperature after the i time virtual heating of calculating steel billet j and the deviation △ T of its optimization heating curve target temperature _{I, j}

S55) judgement △ T _{I, j}Whether less than or equal to allowable value ξ, if greater than allowable value ξ, with △ T _{I, j}Assignment is to △ T _jGetting back to step S53 continues to calculate; If △ is T _{I, j}Less than or equal to allowable value ξ, with furnace temperature Tf _{I, j, k}As the corresponding decision-making furnace temperature of steel billet j;

S56) judge whether steel billet counter j equals Steel In Reheating Furnace base overall number N, if j is not equal to N, steel billet counter j adds 1, then get back to step S52 and continue to calculate; If steel billet counter j equals Steel In Reheating Furnace base overall number N, continue following steps;

S57) according to the furnace temperature decision value Tf of every block of steel billet in process furnace _{I, j, k}Calculate the furnace temperature decision value Tf that obtains each control section of process furnace _k

6. method according to claim 5 is characterized in that:

Step S53 obtains revised furnace temperature Tf _{I, j, k}Method obtain by following formula:

${\mathrm{Tf}}_{i,j,k}={\mathrm{Tf}}_{i,j,k}\·(1-c_i\·\frac{\mathrm{\Δ}{T}_{i,j}}{T_{\mathrm{Opt},j}})$

In formula, the Tf on equation the right _{I, jk}The corresponding furnace temperature of steel billet j when being the i time virtual heating, the Tf on the equation left side _{I, j, k}Be the corresponding furnace temperature of steel billet j of revising after the i time virtual heating, it is the corresponding furnace temperature of steel billet j during as the i+1 time virtual heating, c _iBe correction factor, T _{Opt, j}Be the optimization heating curve target temperature of steel billet j, Δ T _i,jBe temperature and its optimization heating curve target temperature T of steel billet j after the i time virtual heating _{Opt, j}Deviation.

7. method according to claim 5, is characterized in that, the method for the virtual heating of step S54 comprises:

S541) with the initial temperature field T of steel billet j _{0, j}And the corresponding furnace temperature Tf of steel billet j _{I, j, k}Steel billet is carried out location tracking and temperature tracking;

S542) temperature of steel billet j and the deviation △ T of its optimization heating curve target temperature after accounting temperature is followed the tracks of _{I, j}

S543) judgement △ T _{I, j}Whether less than or equal to allowable value ξ:

S5431) if △ is T _{I, j}Less than or equal to allowable value ξ, finish to the virtual heating of steel billet j and with Tf _{I, j, k}As the corresponding decision-making furnace temperature of steel billet j;

S5432) if greater than allowable value ξ, according to △ T _{I, j}To Tf _{I, j, k}Revise, obtain revising the corresponding furnace temperature Tf of rear steel billet j _{I+1, j, k}, get back to step S541 still with T _{0, j}As the initial temperature field of steel billet j, but with the corresponding furnace temperature Tf of revised steel billet j _{I+1, j, k}Continue to calculate.

8. method according to claim 5, is characterized in that, step S57 obtains furnace temperature decision value Tf _kMethod comprise:

Be calculated as follows the furnace temperature decision value Tf that obtains process furnace k control section _k:

${\mathrm{Tf}}_k=\underset{j=1}{\overset{N_k}{\mathrm{\Σ}}}({\mathrm{\ω}}_j\·{\mathrm{Tf}}_{i,j,k})$

In formula, N _kBe the number of steel billet in k control section of process furnace, ω _jBe the furnace temperature decision value Tf of j piece steel billet _{I, j, k}Corresponding weights, Tf _{I, j, k}Be i the furnace temperature decision value of j piece steel billet in k control section of process furnace.

Images