METHOD FOR ESTIMATING THROAT TEMPERATURE OF BLAST FURNACE BASED ON MULTILAYER ORE-TO-COKE RATIO DISTRIBUTION MODEL

Abstract

Disclosed is a method for estimating a blast furnace throat temperature based on a multilayer ore-to-coke ratio distribution model. According to the method, blast furnace equipment parameters and a burden distribution matrix are utilized, the burden layer profile of each layer is calculated according to the burden distribution movement process, a burden layer distribution model is established in combination with the descending process, and the ore-to-coke ratio of each burden layer is obtained. According to the method, the ore-to-coke ratio distribution of multiple layers and main parameters of a blast furnace are used as input, a generalized regression neural network is used for estimating the temperature at the corresponding measurement position of throat temperature, so as to realize the monitoring of throat temperature in the blast furnace process.

Claims

1. A method for estimating a blast furnace throat temperature based on a multilayer ore-to-coke ratio distribution model, comprising the following steps: step (1) obtaining equipment parameters of the blast furnace, burden distribution process parameters, burden parameters, a burden distribution matrix, main operation parameters, main indication parameters of the blast furnace and throat temperature measurement data; step (2) calculating a burden layer profile of each layer according to a burden distribution law, comprising processes of the burden moving from a storage tank to a chute, moving on the chute, descending from the chute to the burden surface, and forming the burden layer profile, which is as follows: step (2.1) discharging the burden from the storage tank and reaching the chute through a central throat, comprising: calculating an initial speed of the burden along the direction of the chute when reaching the chute, according to a known chute length, a chute inclination angle, a central throat length and other parameters based on a law of free fall; step (2.2) calculating a speed of the burden when leaving the chute based on the initial speed in step (2.1) through stress analysis according to the known chute length, a chute rotation speed and a friction coefficient of the burden on the chute; step (2.3) calculating a coordinate position of a tip of a burden pile, which is formed when the burden reaches the burden layer surface, in a radius direction of the blast furnace according to a known chute inclination angle and a burden line height; step (2.4) determining the burden layer profile according to known internal and external burden pile angles and the coordinate position of the tip of the burden pile, wherein an abscissa of the tip of burden pile is determined in step (2.3), and an ordinate is calculated according to the principle that a single-loop burden volume in the burden distribution matrix is equal to a volume between two successive burden layer profiles; step (2.5) defining a burden layer surface formed by a previous inclination angle as a new initial burden layer profile, calculating a burden layer profile function from a second chute to a last chute inclination angle in turn according to the burden distribution matrix, and completing a burden distribution cycle of the burden distribution matrix, and obtaining a final result as the burden layer profile of a certain layer; step (3) calculating the burden layer profile of each layer in the furnace according to layer descending rules to implement an iterative cycle of burden layers, wherein the layers have different descending rules in different positions in the blast furnace, and the descent starts from a lowest layer and goes up layer by layer, and wherein the volume of each layer is calculated by a method of segmentation until a top layer descending is completed, a current burden layer distribution of the blast furnace is recorded, and then a next top layer burden distribution is carried out to prepare for the next descent, which is as follows: step (3.1) carrying out the burden descending process if a height of the burden line of the top layer is higher than a set value, and taking the burden surface of the current top layer as an initial burden layer surface, and calculating the burden layer profile of the top layer, if a height of the burden line of the top layer is not higher than a set value, according to the burden distribution matrix and the step (2); step (3.2) acquiring a burden layer descending trajectory: the burden has different descending trajectories at different positions in the blast furnace, and vertical descent occurs at the throat and bosh of the furnace, while the radial coordinates remain unchanged; the radial and axial movement laws of the burden at a furnace shaft and a furnace waist are calculated according to the principle of similar triangles and a uniform descending mode; step (3.3) calculating a descending volume of each layer: dividing each layer into several triangles according to the shapes of upper and lower interfaces thereof, calculating an area of each triangle, and taking a result of the accumulation of the volumes enclosed by rotation around a center line of the blast furnace as the volume of each layer; step (3.4) defining a descent volume as the burden volume of the last layer when the layer above last descends after the last layer of the burden descends, applying the falling rules in steps (3.2) and (3.3) to each burden layer one by one upwards till the top layer, and defining, after the top layer descends, an upper interface thereof as a new initial burden layer surface and reading the new initial burden layer surface into the next burden distribution matrix, and then returning to step (3.1); building a blast burden furnace layer burden distribution model through step (2) and step (3); step (4) calculating the distribution of the ore-to-coke ratio in each layer by the following calculation formula: ${\mathrm{OCR}}_k\left(x\right)=\begin{array}{cc}\frac{{{\gamma }_o\left(x\right)}_k-{{\gamma }_c\left(x\right)}_{k-1}}{{{\gamma }_c\left(x\right)}_{k-1}-{{\gamma }_o\left(x\right)}_{k-2}}& 1<k\le K\end{array}$ where x represents a distance between a certain point and the center line of the blast furnace, γ(x).sub.k represents a burden surface distribution function of a k.sup.th layer, and the subscripts o and c are used to distinguish an ore layer from a coke layer, and K represents a number of selected burden layers, OCR.sub.k (x) represents the ore-to-coke ratio of the k.sup.th layer, and wherein the selected multilayer burden layers ranges from the throat position to the bosh position; step (5) establishing a throat temperature estimation model based on a generalized regression neural network by taking the ore-to-coke ratio of each layer and the main parameters of the blast furnace as inputs and the measurement point data of the throat temperature as outputs; inputting, after finishing training the model, the current ore-to-coke ratio of each layer and the main parameters of the blast furnace to obtain a throat temperature estimation, which is as below: step (5.1) time-registering the main parameters of the blast furnace and the measurement data of the throat temperature with the burden distribution process, and selecting the main parameters of the blast furnace and the measurement data of the throat temperature which are consistent with a time of the burden distribution matrix; step (5.2) data-preprocessing, comprising data cleaning and normalization; step (5.3) the generalized regression neural network being composed of an input layer, a pattern layer, a summation layer and an output layer, wherein the relationship between an input and an output thereof is expressed by the following formula: $E\left[T.Math.U\right]=\frac{{\int }_{-\infty }^{\infty }\mathrm{Tg}\left(U,T\right)\mathrm{dT}}{{\int }_{-\infty }^{\infty }g\left(U,T\right)\mathrm{dT}}$ where T represents an output result of GRNN, and the input vector U represents a N×1-dimensional vector composed of the ore-to-coke ratio of each layer and the main parameters of the blast furnace, E[T|U] represents an expected value of an output T of a given input vector U and g(U,T) represents a joint probability density function of U and T; step (5.4) training the model by taking the main parameters of the blast furnace and ore-to-coke ratios of every burden layer in the multilayer in a training set as an output vector and a throat temperature value as an output vector; inputting, after finishing training the model, the current main parameters of the blast furnace and ore-to-coke ratios of every burden layer in the multilayer to obtain the estimated value of the throat temperature.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0035] FIG. 1 is a flow chart of a method for estimating a blast furnace throat temperature based on a multilayer ore-to-coke ratio distribution model according to an embodiment of the present disclosure.

[0036] FIG. 2 is a schematic diagram of measuring points of a throat temperature measuring device (a cross temperature measuring device) according to an embodiment of the present disclosure.

[0037] FIG. 3 is a schematic diagram of the dropping process of burden according to an embodiment of the present disclosure.

[0038] FIG. 4 is a schematic diagram of construction of a burden surface function according to an embodiment of the present disclosure.

[0039] FIG. 5 is a schematic diagram of burden descent according to an embodiment of the present disclosure.

[0040] FIG. 6 is a schematic diagram of the method of calculating the volume of the burden layer according to an embodiment of the present disclosure.

[0041] FIG. 7 is a flow chart of establishing a burden layer model according to an embodiment of the present disclosure.

[0042] FIG. 8 is a schematic diagram of a GRNN network structure according to an embodiment of the present disclosure.

[0043] FIG. 9 shows the estimation results of two temperature measurement points by two methods according to an embodiment of the present disclosure.

DESCRIPTION OF EMBODIMENTS

[0044] The present disclosure will be further explained with reference to the drawings and specific embodiments.

[0045] FIG. 1 shows the overall flow of a method for estimating a blast furnace throat temperature based on a multilayer ore-to-coke ratio distribution model according to an embodiment of the present disclosure. The burden layer distribution model of a blast furnace is established by using equipment parameters of the blast furnace, a burden distribution matrix and state parameters of the blast furnace. Data processing of state parameters of the blast furnace and burden distribution matrix includes sampling synchronization, data cleaning and feature selection. The processed state parameters are combined with the multilayer ore-to-coke ratio burden distribution to obtain the required input data. A GRNN algorithm is used to obtain the estimated throat temperature.

[0046] In order to realize the modeling process of the burden layer structure, the following assumptions are made:

[0047] A. the volume change of the burden in the movement process is ignored, and the collapse and deformation of the burden surface is not considered.

[0048] B. the burden surface is central symmetrically distributed in the center.

[0049] C. the burden surface and burden flow remain continuous.

[0050] D. the burden flow speed of the burden distribution in each rotating circle of the chute is kept stable.

[0051] E. the chute and burden flow shall be implemented strictly according to the burden distribution matrix.

[0052] The embodiment of the present disclosure relates to a method for estimating a blast furnace throat temperature based on a multilayer ore-to-coke ratio distribution model, and each step is specifically described as follows:

[0053] (1) Equipment parameters of the blast furnace, burden distribution process parameters, burden parameters, a burden distribution matrix, main operation parameters, main indication parameters of the blast furnace and throat temperature measurement data are obtained as follows:

[0054] The equipment parameters of the blast furnace include a total blast furnace height, a throat height, a bosh height, a body height, a waist height, a throat radius, a body inclination angle, a bosh radius and a waist inclination angle.

[0055] The burden distribution process parameters include a central throat length h.sub.0, a throttle opening S, a chute length l, a chute tilting distance b, a chute rotating speed ω, a chute friction coefficient μ and a burden line depth H.

[0056] The burden parameters include the average particle size D.sub.o,c, the average density ρ, the inner stacking angle φ.sub.in and the outer stacking angle φ.sub.out of the ore coke.

[0057] The burden distribution matrix includes a chute inclination angle the number of rotation turns corresponding to each inclination angle, and the single-circle ore coke volume.

[0058] The main parameters of the blast furnace include an oxygen enrichment rate (%), a permeability index (%), a cold gas flow rate (m3/h), an oxygen enrichment flow rate (m.sup.3/h), a top pressure (kPa), hot and cold air temperatures (° C.), a top temperature (° C.) and a blast humidity (%).

[0059] The schematic diagram of the throat temperature measuring device is shown in FIG. 2, which is installed at the blast furnace throat and detected by a thermocouple sensor.

[0060] (2) The burden layer profile of each layer is calculated according to the burden distribution law and the principle of equal volume, and the specific steps include:

[0061] FIG. 3 shows the concrete schematic diagram of the falling process of the burden.

[0062] (2.1) The burden (ore or coke) is discharged from the hopper and enters the rotary chute through the central throat. Assuming that the movement of the burden before it reaches the chute falls freely, and considering the collision process, the initial velocity of particles entering the chute is:

[00003]$v_1=K_f\cos \beta \sqrt{v_0^2+2g\left(h_0+\frac{b}{\sin \beta }\right)}$

[0063] where v.sub.0 is the speed of particles leaving the hopper and K.sub.f is an impact attenuation coefficient.

[0064] (2.2) When the burden with a mass of m is moving on the chute, it is subjected to various forces, including a gravity mg, a chute reaction force F.sub.N, a friction force F.sub.f, an inertial centrifugal force F.sub.c, an inertial Coriolis force F.sub.Coli, and the force F.sub.L between the burden and the chute sidewall caused by rotation. According to the force analysis of particles and Newton's second law, the speed of the burden leaving the chute is calculated as follows:

[00004]$v_2={\left[2\mathrm{gl}\left(\sin \beta -\mu \cos \beta \right)+v_1^2+4{\pi }^2{\omega }^2{l}^2\cos \beta \left(\cos \beta +\mu \sin \beta \right)\right]}^{\frac{1}{2}}$

[0065] (2.3) The burden undergoes a projectile motion in the empty area, and is affected by its own gravity and drag force of the rising gas. Because the drag force of the gas is negligible, it can be considered that the burden moves with an initial speed of v.sub.2, and a vertical acceleration of g. The distance d between the tip point of the burden pile and the main shaft of the blast furnace in the horizontal direction is calculated by the following formula:

[00005]$t=\frac{L_x}{v_2\cos \beta }=\frac{L_y}{\omega l\cos \beta }H=L_x\tan \beta +\frac{\mathrm{mg}}{2m{v}_2^2{\cos }^2\beta }{L}_x^2-l\left(1-\sin \beta \right)d=\sqrt{{\left(L_x+l\cos \beta \right)}^2+L_y^2}$

[0066] where t is a falling time of the burden, L.sub.x is a projection distance of the burden trajectory in the radial direction and L.sub.y is a projection distance of the burden trajectory in the tangential direction.

[0067] (2.4) The shape of the burden pile is shown in FIG. 4. According to the internal and external stacking angles φ.sub.in and φ.sub.out of the ore coke, the slopes of two line segments of the burden surface can be obtained. The burden layer profile function form is as follows:

[00006]$\gamma \left(x\right)=\{\begin{array}{c}\tan {\phi }_{\mathrm{in}}\left(x-x_{\mathrm{peak}}\right)+y_{\mathrm{peak}},x_{\mathrm{Left}}\le x<x_{\mathrm{peak}}\\ \tan {\phi }_{\mathrm{out}}\left(x-x_{\mathrm{peak}}\right)+y_{\mathrm{peak}},x_{\mathrm{peak}}\le x\le x_{R\mathrm{ight}}\end{array}$

[0068] where (X.sub.peak, Y.sub.peak) is the coordinate of the tip of the burden pile, X.sub.peak is d calculated in (2.3), and X.sub.Left and X.sub.Right radial coordinates of the left and right end points of the pile. y.sub.peak is obtained through calculation according to the principle that the volume of single-circle burden in the burden distribution matrix is equal to the volume between two successive burden layer profiles, thereby obtaining the burden layer profile function.

[0069] (2.5) After burden distribution is completed for all inclination angles of a burden distribution matrix, the result is the burden layer profile of a certain layer (an ore layer or a coke layer). According to FIG. 4, the burden layer profile function of the multi-circle burden distribution mode can be expressed as:

[00007]${\gamma }_i\left(x\right)=\{\begin{array}{cc}{\gamma }_{i-1}\left(x\right),& 0\le x<x_{A_i}\\ \tan {\phi }_{\mathrm{in}}\left(x-x_{C_i}\right)+y_{C_i}& x_{A_i}\le x<x_{C_i}\\ \tan {\phi }_{\mathrm{out}}\left(x-x_{C_i}\right)+y_{C_i}& x_{C_i}\le x<x_{B_i}\\ {\gamma }_{i-1}\left(x\right),& x_{B_i}\le x\le D_0\end{array}$

[0070] where (x.sub.ci, y.sub.ci) is the coordinate of the pile tip C.sub.i with the i.sup.th inclination angle, x.sub.Ai and x.sub.Bi are the radial coordinates of the intersection points of the new burden line and the original burden surface, and D.sub.0 is the radius of the throat.

[0071] (3) The burden layer profile of each layer in the furnace is calculated according to the descending law, and the iterative cycle of the burden layer is realized. The specific steps are as below:

[0072] (3.1) If the height of the burden line of the top layer is higher than the set value, the descending process is carried out. Otherwise, the current burden surface of the top layer is taken as the initial burden layer surface, and the burden layer profile of the top layer is calculated according to the burden distribution matrix and the step (2).

[0073] (3.2) FIG. 5 shows the lowering mode of the burden at the throat and shaft. Point “O” is the intersection of the central axis of the blast furnace and the extension line of the shaft wall, a is the angle of the shaft of the blast furnace and L.sub.throat is the length of the throat area. In the throat area, it is assumed that the original position of the burden is (r,y), and the position after drop of a unit volume ΔV is

[00008]$\left(r,\frac{4\times \Delta V}{\pi D_0^2}+y\right).$

In the furnace shaft area, the burden trajectory becomes along the ray from point O, and the new position (r′,y′) is:

[00009]$y^{\prime }={\left[{\left(\frac{D_0}{2\tan \alpha }-L_{\mathrm{throat}}+y\right)}^3+\frac{3\Delta V}{{\pi \left(\tan \alpha \right)}^2}\right]}^{\frac{1}{3}}-\frac{D_0}{2\tan \alpha }+L_{\mathrm{throat}}$$r^{\prime }=r\times \frac{\frac{D_0}{2\tan \alpha }-L_{\mathrm{throat}}+y^{\prime }}{\frac{D_0}{2\tan \alpha }-L_{\mathrm{throat}}+y}$

[0074] The law of the waist part is consistent with that of the shaft part, and the law of the bosh part is consistent with that of the throat part, so the formula will not be described in detail.

[0075] (3.3) The cross section of each burden layer is shown in FIG. 6 (a), and the calculation method of the burden volume of each layer is: dividing the burden of each layer into several triangles according to the shapes of its upper and lower interfaces, as shown in FIG. 6 (b), calculating the area of each triangle, and taking the accumulated results of the volumes enclosed by its rotation around the centerline of the blast furnace as the volume of each layer.

[0076] (3.4) The layer above last descends after the last layer of the burden descends, and the descent volume is the burden volume of the last layer. By applying the falling rules in (3.2) and (3.3) to each burden layer one by one upwards till the top layer, after the top layer descends, an upper interface thereof is taken as a new initial burden layer surface and read into the next burden distribution matrix, return to (3.1).

[0077] (4) The ore-to-coke ratio distribution in each layer is calculated.

[0078] The ore-to-coke ratio is a parameter describing the radial coke thickness of the blast furnace. Assuming that the last batch of blast furnace raw materials is ore and the penultimate batch is coke, the calculation formula of the ore-to-coke ratio is obtained by combining the above burden layer profile model:

[00010]${\mathrm{OCR}}_k\left(x\right)=\begin{array}{cc}\frac{{{\gamma }_o\left(x\right)}_k-{{\gamma }_c\left(x\right)}_{k-1}}{{{\gamma }_c\left(x\right)}_{k-1}-{{\gamma }_o\left(x\right)}_{k-2}}& 1<k\le K\end{array}$

[0079] where x represents a distance between a certain point and the center line of the blast furnace, γ(x).sub.k represents a burden surface distribution function of a k.sup.th layer, and the subscripts o and c are used to distinguish an ore layer from a coke layer, K represents a number of selected burden layers, OCR.sub.k (x) represents the ore-to-coke ratio of the k.sup.th layer; the selected multiple burden layers ranges from the throat position to the bosh position.

[0080] FIG. 7 shows the flow chart of building the burden layer distribution model, including steps (2), (3) and (4).

[0081] (5) The ore-to-coke ratio of each layer and the main parameters of the blast furnace are taken as inputs and the temperature measurement point data of the throat as outputs, the throat temperature estimation model is established based on the generalized regression neural network. This step consists of the following substeps:

[0082] (5.1) The main parameters of the blast furnace and the measurement data of the throat temperature are time-registered with the burden distribution process, and the main parameters of the blast furnace and the temperature measurement data of the throat which are consistent with a time of the burden distribution matrix are selected.

[0083] (5.2) Data preprocessing, including data cleaning and normalization. In order to ensure the validity and reliability of the data, input/output combinations with invalid or missing values (data such as sensor failure or malfunction) are eliminated from the data.

[0084] (5.3) The structure of a generalized regression neural network is shown in FIG. 8, which is divided into four layers. The first layer is the input neuron, and the input vector U is a N×1-dimensional vector composed of the ore-to-coke ratio of each layer and main parameters of the blast furnace. These variables are expressed in a vector form (u.sub.1, u.sub.2, . . . , u.sub.N).

[0085] After receiving the information, the neurons in the second pattern layer systematically process and combine the data. The number of neurons in the pattern layer is equal to the number of samples included in the selected blast furnace period, and the transfer function of the i.sup.th neuron to input and output processing is:

θ=e.sup.−(U−Ui)T(U−Ui)/2σ.sup.2

[0086] where θ.sub.i is the output of the neuron in the pattern layer, U.sub.i is the input sample vector corresponding to the i.sup.th neuron, and a is a smoothing factor.

[0087] The neurons in the third layer strengthen the output of the second layer, and perform arithmetic summation and weighted summation among the outputs, with the following formula:

W.sub.1=Σ.sub.iθ.sub.i

W.sub.2=Σ.sub.ip.sub.iθ.sub.i

[0088] where W.sub.1 is the arithmetic sum result, W.sub.2 is the weighted sum result and p.sub.i is the weight value corresponding to θ.sub.i.

[0089] After the results are summed and transmitted to the last output neuron, the estimated value of the output measured temperature of the throat is obtained.

T=W.sub.2/W.sub.1

where T is an the estimated value of the output throat temperature.

[0090] The internal logical relationship between the input and output of the GRNN method is represented by the following formula:

[00011]$E\left[T.Math.U\right]=\frac{{\int }_{-\infty }^{\infty }\mathrm{Tg}\left(U,T\right)\mathrm{dT}}{{\int }_{-\infty }^{\infty }g\left(U,T\right)\mathrm{dT}}$

[0091] where E[T|U] is the expected value of the output T of a given input vector U, and g(U,T) is the joint probability density function of U and T.

[0092] (5.4) The model is trained by taking the main parameters of the blast furnace and ore-to-coke ratios of every burden layer in the multilayer in a training set as an output vector and a throat temperature value as an output vector; after the training of the model is completed, the current main parameters of the blast furnace and the ore-to-coke ratios of every burden layer in the multilayer are input to obtain the throat temperature.

[0093] In this embodiment, the main parameter data, throat temperature measuring point data and blast furnace burden distribution matrix of a blast furnace of 2650 m.sup.3 in two months of a certain year and in China are selected. The radius length of the throat part of the blast furnace is 4.15 m, and the temperature measuring device shown in FIG. 2 is used. The method is evaluated by selecting two measuring points A and B at different positions. By establishing the model, the burden layer distribution and the ore-to-coke ratio distribution of each layer are obtained. The ore-to-coke ratio distribution of multiple layers and main parameters are used as inputs. A total of 360 samples are selected, the first 300 samples are used for training and the last 60 samples are used for testing.

[0094] For a more comprehensive analysis and discussion, two methods are used to estimate the temperature:

[0095] (1) Classical method: only the OCR distribution and main state parameters of the top layer are adopted.

[0096] (2) The method provided by the present disclosure: the OCR distribution and main state parameters of multiple layers are adopted.

[0097] Statistic indicators including MAPE (mean absolute percentage error), MAE (Mean Absolute Error) and RMSE (Root Mean Square Error) are used to evaluate the estimation effect of the model. The three indicators are calculated as follows:

[00012]$\mathrm{RMSE}=\sqrt{\frac{1}{M}{.Math.}_{i=1}^M{\left(y\left(i\right)-\tilde{y}\left(i\right)\right)}^2}$$\mathrm{MAE}=\frac{1}{M}{.Math.}_{i=1}^M.Math.y\left(i\right)-\tilde{y}\left(i\right).Math.$$\mathrm{MAPE}=\frac{1}{M}{.Math.}_{i=1}^M.Math.\frac{y\left(i\right)-\tilde{y}\left(i\right)}{y\left(i\right)}.Math.\times 100\%$

[0098] where y(i) is the measured value of throat temperature, {tilde over (y)}(i) is the estimated value of throat temperature and M is the number of samples.

[0099] The result of that two methods are shown in FIG. 9, in which OCR_TOP is the result of the traditional method, and OCR_Multi is the method propose by the present disclosure. The evaluation results are as follows:

TABLE-US-00001 GRNN Position A B The present Improvement The present Improvement Algorithm Traditional disclosure rate (%) Traditional disclosure rate (%) MAPE 7.84% 4.42% 43.62% 5.71% 5.27% 7.71% MAE 4.36° C. 2.45° C. 43.81% 1.55° C. 1.43° C. 7.74% RMSE 5.05° C. 3.02° C. 40.20% 2.07° C. 1.85° C. 10.63%

[0100] As shown in the figure, the ore-to-coke ratio of multiple layers is obviously closer to the measured value than only considering the top layer distribution. The overall trend of the results of the method provided by the present disclosure is consistent with the measured temperature. By changing the dimension of the input vector from a single layer to multiple layers, a better temperature estimation result and a higher accuracy are obtained without over-fitting. The validity and reliability of the proposed method are proved.

[0101] It can be seen from the table that taking the distribution of the OCR of multiple layers as input features can promote the prediction results. The accuracy of all evaluation criteria is the highest at point B, with MAPE of 5.27%, MAE of 1.43° C. and RMSE of 1.85° C. The column about the improvement rates of these two points shows the same results as those in FIG. 9, which intuitively proves that the temperature estimation method proposed by the present disclosure is superior to the traditional method.

[0102] The above specific embodiments have explained the technical solution and beneficial effects of the present disclosure in detail. It should be understood that the above embodiments are only the most preferred embodiments of the present disclosure, and are not intended to limit the present disclosure. Any modification, supplement, equivalent substitution, etc. made within the scope of the principles of the present disclosure should be included in the scope of protection of the present disclosure.