阅读说明：本技术 一种基于黑匣子试验横向分区数据提高加热炉钢坯温度跟踪模型计算精度的方法 (Method for improving heating furnace billet temperature tracking model calculation accuracy based on black box test transverse partition data ) 是由 韦思丞 巫献华 陈功彬 王�琦 黄育坚 胡俊平 易承钧 赵星星 于 2021-08-30 设计创作，主要内容包括：一种基于黑匣子试验横向分区数据提高加热炉钢坯温度跟踪模型计算精度的方法,属于加热炉控制技术领域。该方法将步进式加热炉的横向均分为N个区,N为大于1的整数,采用N个相同的测试钢坯跟踪测量钢坯在所述步进式加热炉内加热过程中的温度变化趋势,分别获取所述步进式加热炉中N个区的测试过程中钢坯温度、炉气温度随时间变化的测试数据；提取由炉气温度、下炉气温度、钢坯上表面温度、钢坯下表面温度、时间点组成的参照数据；由参照数据求出钢坯各层的C-(p)、λ、q～(u)(i)、q～(d)(i),根据钢坯上表面的热流系数q～(u)(i)、钢坯下表面的热流系数q～(d)(i),带入综合热辐射系数计算公式(式1),求出N个区对应的综合热辐射系数ε-(g)。保证所有炉况下的计算精度。(A method for improving the calculation accuracy of a heating furnace billet temperature tracking model based on black box test transverse partition data belongs to the technical field of heating furnace control. The method comprises the steps of transversely dividing a stepping heating furnace into N regions, wherein N is an integer larger than 1, adopting N identical test steel billets to track and measure the temperature change trend of the steel billets in the heating process of the stepping heating furnace, and respectively obtaining test data of the temperature of the steel billets and the temperature of furnace gas changing along with time in the test process of the N regions in the stepping heating furnace; extracting reference data consisting of furnace gas temperature, lower furnace gas temperature, billet upper surface temperature, billet lower surface temperature and time points; by reference toData determination of C for each layer of the billet p 、λ、q u (i)、q d (i) According to the heat flow coefficient q of the upper surface of the steel billet u (i) Heat flux q of the lower surface of the billet d (i) Substituting the formula (formula 1) into the calculation formula of comprehensive thermal emissivity to obtain the comprehensive thermal emissivity epsilon corresponding to the N regions g . And the calculation accuracy under all furnace conditions is ensured.)

1. A method for improving the calculation accuracy of a heating furnace billet temperature tracking model based on black box test transverse partition data is characterized by comprising the following steps: the method comprises the following steps:

step one, transversely dividing a stepping heating furnace into N regions, wherein N is an integer larger than 1, adopting N identical test steel billets to track and measure the temperature change trend of the steel billets in the heating process of the stepping heating furnace, and respectively obtaining test data of the temperature of the steel billets and the temperature of furnace gas changing along with time in the test process of the N regions in the stepping heating furnace;

step two, processing the test data to obtain trend graphs of the temperature of furnace gas above the test billet, the temperature of furnace gas below the test billet, the temperature of the upper surface of the test billet, the temperature of the lower surface of the test billet, the temperature of the head of the test billet and the temperature of the tail of the test billet along with time and the change of the position in the test billet, and extracting reference data consisting of the temperature of the furnace gas below the test billet, the temperature of the upper surface of the test billet, the temperature of the lower surface of the test billet and time points; determining C of each layer of the billet from the reference data_p、λ、q^u(i)、q^d(i)，C_pSpecific heat of each layer of the steel billet, lambda is the heat conduction coefficient of the steel billet, and the heat flow coefficient q of the upper surface of the steel billet^u(i) Heat flux q of the lower surface of the billet^d(i)；

Thirdly, according to the heat flow coefficient q of the upper surface of the steel billet^u(i) Heat flux q of the lower surface of the billet^d(i) Substituting into a calculation formula (formula 1) of comprehensive thermal emissivity to obtainComprehensive thermal emissivity epsilon corresponding to N regions_g；

In the formula, epsilon_gQ is the surface heat flux of the steel billet, F_sgIs the surface area of the billet, T_gIs the temperature of furnace gas, T_sThe surface temperature of the steel billet is shown, a is thermal diffusivity, delta is the coordinate position of a test point, and y is the coordinate axis of the steel billet along the thickness direction;

step four, when the actual steel billet enters the furnace, the model automatically indexes out the comprehensive thermal emissivity epsilon corresponding to the N areas of the slab according to the chemical components contained in the steel billet data and the temperature calculated in the previous period_gCarrying the temperature data into a slab temperature tracking model to calculate furnace gas temperature, lower furnace gas temperature, upper surface temperature of the steel billet and lower surface temperature of the steel billet at the time of discharging the steel billet corresponding to the N areas;

and fifthly, selecting the lowest values of the upper furnace gas temperature, the lower furnace gas temperature, the upper surface temperature of the steel billet and the lower surface temperature of the steel billet which correspond to the N areas and are not at the same time, and determining control targets of the upper furnace gas temperature, the lower furnace gas temperature, the upper surface temperature of the steel billet and the lower surface temperature of the steel billet.

2. The method of claim 1, wherein C is determined for each layer of the slab from the reference data_pThe specific method of lambda is that the reference data of N areas are respectively substituted into the heat conduction differential equation (formula 2) in the billet,

wherein rho is billet density and C_pThe specific heat of the steel billet, tau, T and lambda are respectively the heating time, the temperature and the heat conduction coefficient of the steel billet.

3. The method of claim 1Characterized in that the heat flow coefficient q of the upper surface of the billet is^u(i) Heat flux q of the lower surface of the billet^d(i) The calculation formula of (2) is as follows:

wherein the content of the first and second substances,

in the formula (I), the compound is shown in the specification,the Fourier coefficient of the j-th layer of the steel billet is obtained after the steel billet is layered in the thickness direction; lambda [ alpha ]_j,j-1Represents the equivalent heat conduction coefficient between the j-1 layer and the j-th layer of the billet, W/(m.K); delta tau is the time interval, s, calculated by the billet temperature model; rho is the density of the steel billet, Kg/m 3;C_pthe specific heat of the steel billet is J/(kg.K); delta y is the billet layering thickness, m;the temperature of the upper surface of the steel billet;the temperature of the lower surface of the steel billet is shown.

4. The method of claim 1, wherein the lateral direction of the walking-beam furnace is divided into two zones.

5. The method of claim 1, wherein the walking beam furnace has an effective furnace length of 28600mm and a width of 16500 mm.

6. The method according to claim 1, wherein the length direction of the walking beam furnace comprises a preheating section, a heating section and a soaking section in sequence.

7. The method of claim 1, wherein the test billet dimensions are: thickness, width, length 150mm 2100mm 7500 mm.

8. The method of claim 1 wherein said test steel blank is selected from the group consisting of Q235 steel, cold blank.

9. The method of claim 1, wherein the temperature of the test billet during the test is obtained by embedding thermocouples in the upper surface and the lower surface of the test billet.

10. The method of claim 9, wherein the position of the embedded thermocouple is 15mm from the surface of the billet.

Technical Field

The application relates to the technical field of heating furnace control, in particular to a method for improving the calculation accuracy of a heating furnace billet temperature tracking model based on black box test transverse partition data.

Background

The steel billet temperature tracking model is a thermal mathematical model and is used for simulating and calculating the change of a temperature field when a steel billet is heated in a heating furnace. The mathematical model technology of the heating furnace mainly comprises functional models of billet heating curve optimization, billet temperature tracking calculation, furnace temperature set value calculation, combustion optimization control, self-adaptive correction and the like. In actual production, the temperature tracking model calculates the temperature field data of the steel billet according to the current incoming material information (including the size, the type, the incoming temperature and the like of the steel billet) and the current furnace condition information (the current air gas flow, the air gas valve position value, the furnace temperature set value and the like), the furnace temperature set value model calculates the furnace temperature set value according to the temperature field data of the steel billet, then the furnace temperature set value is transmitted to the combustion optimization control model to calculate the optimization parameters and the air gas valve position regulating quantity, and the air gas flow is directly regulated and controlled. Therefore, the calculation accuracy of the billet temperature tracking model ultimately affects the heating quality of the billet and the control level of energy consumption.

The development of the temperature tracking model of the steel billet in the stepping heating furnace is mainly a model developed based on data mining and statistical theory, and mainly comprises two modules of data acquisition and data mining analysis. The data processing mainly finishes real-time data acquisition to obtain the current heating system of the heating furnace, then the acquired production data are processed into analysis samples through a data mining analysis module and are stored in a database according to certain rules, and the optimal furnace temperature set value of the current state of the heating furnace is obtained from the database during production and the temperature of steel billets in the furnace is predicted. The method is biased to a thermal engineering professional theory, and the furnace temperature decision and the furnace temperature forecast cannot fully consider the actual furnace condition information, so that the purposes of optimal control of the heating furnace, energy conservation and consumption reduction obviously cannot be reasonably realized in the practical production.

The first prior art discloses a heating furnace billet temperature tracking model correction method, which comprises the following steps: tracking, measuring and recording data of the change of the temperature of the steel billet and the temperature of furnace gas along with time in the test process; processing the test data, and extracting reference data for calculating a radiation coefficient reference value; combining the reference data, and iteratively calculating the radiation coefficient reference value of each furnace area; correcting the radiation coefficient reference value by combining the production historical data of the steel billet; and storing the correction results in a correction coefficient table in a classified manner. The method can improve the calculation precision of the billet temperature tracking model, so that the corrected model can accurately simulate the heating process of the billet in the furnace, and has important significance for improving the control precision of the furnace temperature of the heating furnace. The method mainly corrects the radiation coefficient continuously through a model so as to improve the calculation accuracy continuously, but the method has a complex calculation process, and when the furnace temperature data fluctuation is large along with different steel types and different furnace conditions, the coefficient correction generates a large error, so that the calculation accuracy under all the furnace conditions cannot be ensured.

The invention aims to solve the problem that aiming at the defects of the prior art, the method for improving the calculation accuracy of the heating furnace billet temperature tracking model is obtained by processing data based on the furnace information acquired by the black box test.

Disclosure of Invention

In order to solve the problems in the prior art, the application aims to provide a method for improving the calculation accuracy of a heating furnace billet temperature tracking model based on the horizontal partition data of a black box test, so that the heating furnace is always controlled above the optimal working temperature of the billet, a simple method for improving the calculation accuracy of the model is provided, the problems that different steel types have different furnace conditions, and when the furnace temperature data fluctuates greatly, the coefficient correction generates large errors and the like are solved, and the calculation accuracy under all the furnace conditions is ensured.

The application example provides a method for improving the calculation accuracy of a heating furnace steel billet temperature tracking model based on transverse partition data of black box test, which comprises the following steps:

step one, transversely dividing a stepping heating furnace into N regions, wherein N is an integer larger than 1, adopting N identical test steel billets to track and measure the temperature change trend of the steel billets in the heating process of the stepping heating furnace, and respectively obtaining test data of the temperature of the steel billets and the temperature of furnace gas changing along with time in the test process of the N regions in the stepping heating furnace;

step two, processing the test data to obtain trend graphs of the temperature of furnace gas above the test billet, the temperature of furnace gas below the test billet, the temperature of the upper surface of the test billet, the temperature of the lower surface of the test billet, the temperature of the head of the test billet and the temperature of the tail of the test billet along with time and the change of the position in the test billet, and extracting reference data consisting of the temperature of the furnace gas below the test billet, the temperature of the upper surface of the test billet, the temperature of the lower surface of the test billet and time points; determining C of each layer of the billet from the reference data_p、λ、q^u(i)、q^d(i)，C_pSpecific heat of each layer of the steel billet, lambda is the heat conduction coefficient of the steel billet, and the heat flow coefficient q of the upper surface of the steel billet^u(i) Heat flux q of the lower surface of the billet^d(i)；

Thirdly, according to the heat flow coefficient q of the upper surface of the steel billet^u(i) Heat flux q of the lower surface of the billet^d(i) Substituting the formula (formula 1) into the calculation formula of comprehensive thermal emissivity to obtain the comprehensive thermal emissivity epsilon corresponding to the N regions_g；

In the formula, epsilon_gQ is steel for comprehensive heat radiation coefficientSurface heat flux density of the blank, F_sgIs the surface area of the billet, T_gIs the temperature of furnace gas, T_sThe surface temperature of the steel billet is shown as a, the thermal diffusivity is shown as a, delta is the coordinate position of the test point, and y is the coordinate axis of the steel billet along the thickness direction.

Step four, when the actual steel billet enters the furnace, the model automatically indexes out the comprehensive thermal emissivity epsilon corresponding to the N areas of the slab according to the chemical components contained in the steel billet data and the temperature calculated in the previous period_gCarrying the temperature data into a slab temperature tracking model to calculate furnace gas temperature, lower furnace gas temperature, upper surface temperature of the steel billet and lower surface temperature of the steel billet at the time of discharging the steel billet corresponding to the N areas;

and fifthly, selecting the lowest values of the upper furnace gas temperature, the lower furnace gas temperature, the upper surface temperature of the steel billet and the lower surface temperature of the steel billet which correspond to the N areas and are not at the same time, and determining control targets of the upper furnace gas temperature, the lower furnace gas temperature, the upper surface temperature of the steel billet and the lower surface temperature of the steel billet.

In some examples, C of each layer of the billet is determined from the reference data_pThe specific method of lambda is that the reference data of N areas are respectively substituted into the heat conduction differential equation (formula 2) in the billet,

wherein rho is billet density and C_pThe specific heat of the steel billet, tau, T and lambda are respectively the heating time, the temperature and the heat conduction coefficient of the steel billet.

In some examples, the heat flow coefficient q of the top surface of the steel billet^u(i) Heat flux q of the lower surface of the billet^d(i) The calculation formula of (2) is as follows:

wherein the content of the first and second substances,

in the formula (I), the compound is shown in the specification,the Fourier coefficient of the j-th layer of the steel billet is obtained after the steel billet is layered in the thickness direction; lambda [ alpha ]_j,j-1 represents the equivalent heat conduction coefficient between the j-1 layer and the j-th layer of the billet, W/(m.K); delta tau is the time interval, s, calculated by the billet temperature model; rho is the density of the steel billet, Kg/m 3; c_pThe specific heat of the steel billet is J/(kg.K); delta y is the billet layering thickness, m;the temperature of the upper surface of the steel billet;the temperature of the lower surface of the steel billet is shown.

In some examples, the lateral direction of the walking beam furnace is divided into two zones.

In some examples, the walking beam furnace has an effective furnace length of 28600mm and a width of 16500 mm.

In some examples, the length direction of the walking beam furnace sequentially comprises a preheating section, a first heating section, a second heating section and a soaking section.

In some examples, the test billet dimensions are: thickness, width, length 150mm 2100mm 7500 mm.

In some examples, the test steel blank is selected to be Q235 steel, cold.

In some examples, the temperature of the test billet during the test is obtained by embedding thermocouples in the upper surface and the lower surface of the test billet.

In some examples, the buried thermocouple is located 15mm from the surface of the billet.

The method for improving the calculation accuracy of the heating furnace steel billet temperature tracking model based on the black box test transverse partition data has the advantages that:

the comprehensive thermal radiation coefficient epsilon of each subarea is obtained by taking the furnace information acquired by the black box test of the transverse subarea of the heating furnace as the data base_gThe method obtains the temperatures of different subareas, determines the lowest temperature of the heating furnace by the temperature, ensures that the heating furnace is always controlled above the optimal working temperature of the steel billet, provides a simple method for improving the calculation accuracy of the model, avoids the problems that the coefficient correction generates larger errors when different steel types and different furnace conditions have larger furnace temperature data fluctuation, and the heating uniformity of the steel billet is influenced due to the different distribution positions of burners in the heating furnace, and the like, and ensures the calculation accuracy under all the furnace conditions.

Drawings

In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are required to be used in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered as limiting the scope, and for those skilled in the art, other related drawings can be obtained from the drawings without inventive effort.

FIG. 1 is a coordinate system of differential equation of internal heat conduction of a steel billet;

FIG. 2 is a schematic view of the position of a test billet within a furnace of a heating furnace.

Icon: 1-heating a furnace; 2-testing the steel billet; 3-paragraph number 1; 4-paragraph number 2; 5-paragraph number 3; 6-paragraph number 4; 7-sub-segment number 1; 8-subsection number 2; 9-furnace width; 10-test billet length.

Detailed Description

Embodiments of the present application will be described in detail below with reference to examples, but those skilled in the art will appreciate that the following examples are only illustrative of the present application and should not be construed as limiting the scope of the present application. The examples, in which specific conditions are not specified, were conducted under conventional conditions or conditions recommended by the manufacturer. The reagents or instruments used are not indicated by the manufacturer, and are all conventional products available commercially. The detection method is not particularly described, and the detection is carried out according to national standards or conventional detection methods.

The following specifically explains the method for improving the calculation accuracy of the heating furnace billet temperature tracking model based on the black box test transverse partition data in the embodiment of the present application:

considering that the distribution positions of the burners in the heating furnace are different, the distribution of furnace gas in the same furnace section is different, and finally the heating uniformity of the billet is influenced, so the acquired data can be closer to the real condition of the furnace condition by adopting the method for testing the billet by sections. The following combination is that the heating furnace 1 has four discharging situations in actual production (namely, the heating furnace 1 is transversely divided into two subareas in figure 2, namely, a furnace gas A area 7 (subsection number 1) and a furnace gas B area 8 (subsection number 2) are respectively provided with two test steel billets 10, and the two subareas are uniformly divided from the position of the central line of the furnace), the heating furnace 1 is divided into two areas for experiment, the data comparison is more obvious, the distribution of double discharging and four discharging of the steel billets in the furnace in actual production is also met, and the experiment cost is reduced.

A method for improving the calculation accuracy of a heating furnace billet temperature tracking model based on transverse partition data of black box test tests comprises the following steps:

step one, transversely dividing the walking beam furnace into two areas, tracking and measuring the temperature change trend of the steel billet 2 in the heating process of the walking beam furnace 1 by adopting two identical test steel billets 2, and respectively obtaining the test data of the change of the steel billet temperature and the furnace gas temperature along with time in the test process of the two areas in the walking beam furnace.

Test data: step-by-step heating furnace 1: the effective furnace length is 28600mm, and the furnace width 9 is 16500 mm.

The length direction of the stepping heating furnace sequentially comprises a preheating section, a heating section and a soaking section, and the corresponding sections are (section number 1)3, (section number 2)4, (section number 3)5 and (section number 4) 6.

The material of the test steel billet is Q235 steel and cold billet, and the size of the test steel billet is as follows: thickness, width, test billet length 10-150 mm 2100mm 7500 mm.

The test steel billet obtains the temperature of the steel billet in the test process through black box experiments by embedding thermocouples on the upper surface and the lower surface of the steel billet, wherein the positions of the embedded thermocouples are 15mm away from the surface of the steel billet, and the inventor finds that: the position can reduce the influence of furnace gas temperature, more accurately reflect the temperature of the upper surface and the lower surface of the billet simultaneously, and improve the calculation precision of the model.

And step two, processing the test data to obtain trend graphs of the temperature of furnace gas above the test billet, the temperature of furnace gas below the test billet, the temperature of the upper surface of the test billet, the temperature of the lower surface of the test billet, the temperature of the head of the test billet and the temperature of the tail of the test billet along with time and the change of the position in the test billet, and extracting reference data consisting of the temperature of the furnace gas below the test billet, the temperature of the upper surface of the test billet, the temperature of the lower surface of the test billet and time points from the trend graphs, wherein the reference data refer to table 1.

TABLE 1 Black box experiment table of the change of each temperature of steel billet with time

Carrying out calculation by taking the data into the formula 2 and the formula 3;

determining C of each layer of the billet from the reference data_pThe specific method of lambda is that the reference data of the two areas are respectively substituted into the heat conduction differential equation (formula 2) in the billet,

wherein rho is billet density and C_pThe specific heat of the steel billet, tau, T and lambda are respectively the heating time, the temperature and the heat conduction coefficient of the steel billet.

Heat flow coefficient q of the upper surface of the billet^u(i) Heat flux q of the lower surface of the billet^d(i) The calculation formula of (2) is as follows:

wherein the content of the first and second substances,

in the formula (I), the compound is shown in the specification,the Fourier coefficient of the j-th layer of the steel billet is obtained after the steel billet is layered in the thickness direction; lambda [ alpha ]_j,j-1Represents the equivalent heat conduction coefficient between the j-1 layer and the j-th layer of the billet, W/(m.K); delta tau is the time interval, s, calculated by the billet temperature model; rho is the density of the steel billet, Kg/m 3; c_pThe specific heat of the steel billet is J/(kg.K); delta y is the billet layering thickness, m;the temperature of the upper surface of the steel billet;the temperature of the lower surface of the steel billet is shown.

Calculating C of each layer of the billet by the calculation_p、λ、q^u(i)、q^d(i)，C_pSpecific heat of each layer of the steel billet, lambda is the heat conduction coefficient of the steel billet, and the heat flow coefficient q of the upper surface of the steel billet^u(i) Heat flux q of the lower surface of the billet^d(i)。

Thirdly, according to the heat flow coefficient q of the upper surface of the steel billet^u(i) Heat flux q of the lower surface of the billet^d(i) And substituting the comprehensive emissivity calculation formula (formula 1).

In the formula, epsilon_gQ is the surface heat flux of the steel billet, F_sgIs the surface area of the billet, T_gIs the temperature of furnace gas, T_sIs the surface temperature of the steel billet, a is heatAnd the diffusivity is delta is the coordinate position of the test point, and y is the coordinate axis of the billet along the thickness direction.

Calculating the comprehensive thermal emissivity epsilon corresponding to the two regions_gSee table 2 for emissivity of each furnace section.

The above-mentioned comprehensive thermal emissivity coefficient epsilon_gThe method comprehensively considers the factors such as the coordinates of the test points, the surface area of the steel billet, the heat diffusion rate, the heat flow density and the like, can more comprehensively represent the heat flow characteristics of the steel billet, and improves the calculation precision.

TABLE 2 radiation coefficient of each furnace section

Step four, when the actual steel billet enters the furnace, the model automatically indexes out the comprehensive thermal emissivity epsilon corresponding to the two areas of the slab blank according to the chemical components contained in the steel billet data and the temperature calculated in the previous period_gAnd carrying into a slab temperature tracking model to calculate furnace gas temperature, lower furnace gas temperature, upper surface temperature of the steel billet and lower surface temperature of the steel billet corresponding to the two zones at the moment of discharging the steel billet, taking the following example as model data calculated by the zone A, and referring to a table 3 of time-dependent change table of each temperature of the slab temperature tracking model steel billet. The method for calculating the time variation of each temperature of the slab temperature tracking model billet is the prior art which can be directly adopted by related technicians and is not described in detail herein.

TABLE 3 slab temperature tracking model billet temperature variation with time table

And fifthly, selecting the lowest values of the upper furnace gas temperature, the lower furnace gas temperature, the upper surface temperature of the steel billet and the lower surface temperature of the steel billet at different moments corresponding to the area A and the area B, and determining the control targets of the upper furnace gas temperature, the lower furnace gas temperature, the upper surface temperature of the steel billet and the lower surface temperature of the steel billet.

And (3) checking and calculating the precision:

because 11:41 is divided into the moment of billet tapping, the error between the model calculation value and the actual value at the moment can be compared to obtain the accuracy of the model calculation, and the accuracy is as follows:

upper furnace gas	Bottom furnace gas	Upper surface of	Lower surface	Time of day
						1228.9	1226.1	1206.1	1186.8	11:41	Actual value
1223.5	1216.3	1206.4	1175.6	11:41	Calculated value
						5.4	9.8	0.3	11.2		Error of the measurement

As can be seen from the data in the table, the error range between the calculated value and the actual value of the model in the area A at the tapping time is within 11.2 degrees, and the temperature change trend of the steel billet in the furnace is accurately predicted. The error range of the calculated value and the actual value of the same B zone model is controlled within 15 degrees, so that when the transverse working conditions of the heating furnace are different, a prepared calculation result can be obtained by the model, and the temperature change trend of the steel billet in the heating furnace can be accurately predicted.

The comprehensive thermal radiation coefficient epsilon of each subarea is obtained by taking the furnace information acquired by the black box test of the transverse subarea of the heating furnace as the data base_gThe method obtains the temperatures of different subareas, determines the lowest temperature of the heating furnace by the temperature, ensures that the heating furnace is always controlled above the optimal working temperature of the steel billet, provides a simple method for improving the calculation precision of the model, avoids the problems that the coefficient correction generates larger errors when different steel types and different furnace conditions have different furnace temperature data fluctuation, and the heating uniformity of the steel billet is influenced due to different distribution positions of burners in the heating furnace, and the like, and ensures the calculation precision under all furnace conditions. Obtaining the comprehensive thermal emissivity epsilon of each subarea according to a black box test_gA reference data set is formed, the accuracy fluctuates up and down due to repeated adjustment of the emissivity is not needed, the applicability of the model is increased, the model is suitable for the control of heating furnaces of all steel types, and the problem of inconsistent control standards caused by the steel types is avoided.

The foregoing is merely exemplary of the present application and is not intended to limit the present application, which may be modified or varied by those skilled in the art. Any modification, equivalent replacement, improvement and the like made within the spirit and principle of the present application shall be included in the protection scope of the present application.