Hearth lining erosion line online prediction method

Abstract

The invention relates to the field of blast furnace ironmaking, in particular to a hearth and lining erosion line online prediction method. The method comprises the following steps: reading thermocouple data arranged on a hearth, carrying out three-dimensional visualization on the thermocouple data, and extracting internal temperature data of a closed area; establishing a 2D geometric model of the hearth and discretizing the 2D geometric model to obtain a discretized grid; interpolating the internal boundary condition by using the internal temperature data of the closed area to obtain the internal boundary condition; based on the pre-obtained real-time tapping temperature and cooling water temperature difference, external boundary conditions of the three-dimensional visualization model are set, internal boundary conditions are set, and a heat conduction model is established through the external boundary conditions and the internal boundary conditions; and solving the discretized grid by using a heat conduction model, and integrating the obtained results. According to the method, the temperatures at the two ends serve as input, the intermediate temperature isotherm is directly solved, and therefore the calculation efficiency is much higher.

Description

Technical Field

[0001] The present invention relates to the field of blast furnace ironmaking, and particularly to an on-line prediction method for the erosion line of the hearth lining. Background Art

[0002] The remaining thickness of the hearth lining directly affects the safety of blast furnace operation. Especially in the later stage of the furnace campaign, the lining has been severely eroded, and major safety accidents such as the burning through of the hearth base are likely to occur. To avoid such accidents and extend the furnace campaign life as much as possible, blast furnace practitioners need to closely monitor the lining status in real time. Limited by the high temperature, high pressure, and toxic environment inside the hearth, it is quite difficult to directly observe the lining status. A relatively common method is to use thermocouples to collect data or information such as the temperature difference of cooling water to infer the internal thermal state of the hearth, and use the solidification temperature of hot metal as a reference to judge the erosion position of the hearth. Research believes that when the lining temperature is lower than the solidification temperature of hot metal, the hot metal will gradually solidify to form a solidified iron layer, thus avoiding the direct contact between the lining and the slag, and the lining erosion stops at this time. Therefore, the prediction of the lining erosion line can be realized by calculating the isotherm of the hot metal solidification temperature (usually 1150°C).

[0003] For the prediction of the hearth thermal state, the existing technologies can be roughly divided into two categories, namely "forward prediction" and "inverse prediction". The so-called forward prediction is to directly predict the thermal state inside the hearth (usually the whole hearth) based on the energy conservation equation, and compare the obtained results with the thermocouple measurement data, and continuously adjust parameters such as the lining thickness and thermal conductivity according to the comparison situation until the calculation results match the measurement results. Obviously, the forward prediction calculation has problems such as low calculation efficiency and random parameter adjustment. The inverse prediction refers to directly using the thermocouple data as the boundary condition to calculate the lining temperature field. Generally, the data of the first layer of thermocouples are used as the boundary, and the lining thickness is continuously adjusted so that the outermost temperature of the lining is just equal to the solidification temperature of hot metal. The inverse prediction is more intuitive than the forward prediction. And because the measured values are introduced, the credibility of the prediction is guaranteed, and the calculation efficiency is also relatively higher. However, the inverse prediction generally only considers the heat transfer in the area between the first layer of thermocouples and the erosion line, and the thermocouples are buried deep inside the lining, so the temperature distribution of the lining cannot be fully reflected. In addition, during the calculation process, the adjustment of the lining thickness is usually carried out in a fixed direction (such as the vertical direction), but the heat transfer process inside the hearth is not isotropic, so the rationality and accuracy of this adjustment method also need to be verified. The most crucial thing is that whether it is forward prediction or inverse prediction, the existing calculation methods only consider one working condition. Although the input parameters can be continuously adjusted as production progresses, it is essentially a "steady-state" model and cannot reflect the dynamic heat transfer process inside the hearth due to the periodic discharge and accumulation characteristics. Summary of the Invention

[0004] To solve the above technical problems existing in the prior art, the present invention provides an online prediction method for the erosion line of the hearth lining.

[0005] To achieve the above object, the embodiments of the present invention provide the following technical solutions: In a first aspect, in an embodiment provided by the present invention, an online prediction method for the erosion line of the hearth lining is provided. The method includes the following steps: Read the thermocouple data arranged on the hearth, perform three-dimensional visualization on the thermocouple data, and extract the internal temperature data of the closed area; Establish a 2D geometric model of the hearth and discretize it to obtain a discretized grid; Interpolate the internal boundary conditions using the internal temperature data of the closed area to obtain the internal boundary conditions; Based on the pre-acquired real-time tapping temperature and the temperature difference of the cooling water, set the external boundary conditions of the three-dimensional visualization model. Use the set internal boundary conditions, and establish a heat conduction model using the external boundary conditions and the internal boundary conditions; Solve the discretized grid using the heat conduction model and integrate the obtained results.

[0006] As a further solution of the present invention, the thermocouple data includes position and temperature information.

[0007] As a further solution of the present invention, the step of reading the thermocouple data arranged on the hearth, performing three-dimensional visualization on the thermocouple data, and extracting the internal temperature data of the closed area includes: Display the coordinates of the thermocouple group based on the thermocouple data, enclose the area including all the acquisition points, and perform interpolation display on the temperature in the enclosed area.

[0008] Specifically, consider Figure 3 the heat transfer process shown. Molten iron and slag drip from the upper part into the hearth and gradually accumulate and are periodically discharged from the taphole. The hearth molten pool can be divided into a dynamic area and a static area according to the height of the taphole. In the dynamic area, due to the flow of molten iron, the heat transfer is mainly by heat convection; while in the static area, since it is below the taphole and the internal flow is relatively static, the internal heat transfer is mainly by heat conduction. Due to the relative stillness in the static area, the molten iron is likely to solidify at the bottom to form a solidified iron layer. Under the action of chemical erosion and mechanical scouring, the hearth lining is continuously eroded to form an erosion line. In the static area, the chemical erosion is the main, and generally 1150°C is used as the reference line. This temperature is the solidification temperature of molten iron. Below this temperature, due to the formation of the solidified iron layer, the erosion will not continue.

[0009] Based on the above characteristics, the hearth can be regarded as Figure 4The calculation domain only considers the static zone and the furnace lining. Heat conduction is the main mode in this calculation domain, so it can be regarded as a whole, and the influence of the furnace lining properties and the condensed iron layer on heat transfer is mainly reflected in its thermal conductivity. The dynamic zone is simplified as a boundary condition, assuming that there is no temperature gradient inside the dynamic zone, and the temperature at the static and dynamic interface is the molten iron temperature in the furnace. Considering the heat loss in the iron discharge process, it is generally believed that the molten iron temperature inside the furnace is 50°C lower than the iron discharge temperature. The molten iron temperature in the furnace is considered to be between 1480-1520°C. It should be pointed out that since different thermal conductivity coefficients are assigned to the molten iron, furnace lining and condensed iron layer, and this coefficient is a function of position and temperature, there is no need to dynamically adjust the thickness of the furnace lining. The predicted 1150°C isotherm is the erosion line under the current calculation conditions. The biggest difference between the method described in the present invention and the existing "reverse calculation" method is that the existing method uses the furnace lining thickness and the temperature at one end as input to solve the temperature at the other end, and continuously adjusts the furnace lining thickness to make the solved temperature reach a preset value; while the existing method uses the temperatures at both ends as input to directly solve the intermediate temperature isotherm, so the calculation efficiency is much higher.

[0010] As a further solution of the present invention, the reading of the thermocouple data arranged on the furnace, three-dimensional visualization of the thermocouple data, and extraction of the temperature data inside the closed area include: Based on the thermocouple data, the coordinates of the thermocouple group are displayed by 3D visualization software, and the area including all the acquisition points is closed, and the temperature in the closed area is interpolated and displayed using the griddata interpolation algorithm.

[0011] As a further solution of the present invention, the step S20, establishing and discretizing a 2D geometric model of the furnace, comprises: Depending on the needs of on-site monitoring, select a specific 2D section for geometric modeling and discretize it according to the calculation accuracy requirements. The maximum grid size cannot exceed 1 / 2 of the minimum furnace lining material thickness; The geometric models under normal and limit conditions are established respectively.

[0012] Specifically, consider Figure 5 The extreme working condition shown is used to simulate the phenomenon of enhanced convective heat transfer at the bottom of the molten pool caused by factors such as the flow of molten iron "annular gap", natural convection caused by temperature difference of molten iron, and liquid level fluctuation caused by drastic pressure fluctuation in the furnace during iron discharge (such as large discharge). Since the above phenomenon cannot be quantified, it is considered in the most extreme way, that is, the entire molten pool is regarded as a dynamic area, and its temperature is the molten iron temperature. The erosion line predicted under this working condition is the maximum erosion line that can be achieved in theory. The actual erosion line should be between Figure 4 and Figure 5 Between the two working conditions.

[0013] As a further solution of the present invention, the interpolation of the internal boundary conditions by using the internal temperature data of the closed region to obtain the internal boundary conditions includes: Interpolate the internal boundary conditions by using the internal temperature data of the closed region to obtain the temperature data of all nodes in the overlapping region.

[0014] As a further solution of the present invention, the step S50 of solving the discretized grid by using the heat conduction model and integrating the obtained results includes: Solve the thermal state of the hearth based on the discretized grid.

[0015] As a further solution of the present invention, solving the thermal state of the hearth includes calculating the node temperature. The calculation formula of the node temperature is as follows:

[0016] In the formula, is the node temperature, i and j respectively represent the nodes at the positions of the i th row and j th column, k is the thermal conductivity of the material.

[0017] The technical solution provided by the present invention has the following beneficial effects: The method of the present invention can completely display the temperature field distribution of the furnace lining from the inside to the outside; and different from the existing prediction methods, the method of the present invention mixes two calculation ideas of "forward prediction" and "reverse prediction", directly regards the thermocouple-covered area as the internal boundary, and the temperature value inside it is directly determined by the measurement data, while the remaining furnace lining areas are calculated based on the heat conduction model. In addition, since the calculation domain is extended and the heat conduction of the static area at the bottom of the molten pool is considered, the isotherm of the molten iron condensation temperature can be obtained without dynamically adjusting the furnace lining thickness, and the calculation efficiency is improved compared with the existing methods; by superimposing the results of "normal" and "extreme" two working conditions, the concept of "erosion zone" is used to replace the "erosion line" in the prior art to predict the erosion state of the furnace lining, reducing the influence of the dynamic heat transfer process of the hearth on the calculation results and enhancing the confidence level of the prediction results. The present invention takes the temperatures at both ends as the input and directly solves the intermediate temperature isotherm, so the calculation efficiency is much higher.

[0018] These aspects or other aspects of the present invention will be more clearly understood in the following description of the embodiments. It should be understood that the above general description and the following detailed description are only exemplary and explanatory, and cannot limit the present invention. Description of the Drawings

[0019] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required for the description of the embodiments or the prior art. Obviously, the drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other embodiments can be obtained based on these drawings.

[0020] Figure 1 It is a flowchart of an online prediction method for the erosion line of the hearth lining according to an embodiment of the present invention.

[0021] Figure 2 It is a specific flowchart example of an online prediction method for the erosion line of the hearth lining according to an embodiment of the present invention.

[0022] Figure 3 It is a schematic diagram of heat conduction inside the hearth according to an embodiment of the present invention.

[0023] Figure 4 It is a diagram of the division of the calculation domain under normal conditions according to an embodiment of the present invention.

[0024] Figure 5 It is a diagram of the division of the calculation domain under extreme conditions according to an embodiment of the present invention.

[0025] Figure 6 It is a diagram of the interpolation method for the thermocouple coverage area according to an embodiment of the present invention.

[0026] Figure 7 It is a schematic diagram of the grid nodes of the finite difference method according to an embodiment of the present invention. Detailed implementation manners

[0027] Next, the technical solutions in the embodiments of the present invention will be clearly and completely described in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are some, but not all, of the embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts belong to the scope of protection of the present invention.

[0028] The flowcharts shown in the drawings are only illustrative examples, and do not necessarily include all the contents and operations / steps, nor do they necessarily need to be executed in the described order. For example, some operations / steps can be decomposed, combined, or partially merged, so the actual execution order may change according to the actual situation.

[0029] It should be understood that the terms used in the specification of the present invention are only for the purpose of describing specific embodiments and are not intended to limit the present invention. As used in the specification of the present invention and the appended claims, unless the context clearly indicates otherwise, the singular forms "a", "an" and "the" are intended to include the plural forms.

[0030] Specifically, in combination with the accompanying drawings, embodiments of the present invention will be further described below.

[0031] Please refer to Figure 1 , Figure 1 which is a flowchart of an online prediction method for the erosion line of the hearth lining provided by an embodiment of the present invention. As Figure 1 shown, the online prediction method for the erosion line of the hearth lining includes steps S10 to S50.

[0032] S10. Read the thermocouple data arranged on the hearth, perform three-dimensional visualization on the thermocouple data, and extract the temperature data inside the closed area.

[0033] The thermocouple data includes position and temperature information.

[0034] In an embodiment of the present invention, the step S10 of reading the thermocouple data arranged on the hearth, performing three-dimensional visualization on the thermocouple data, and extracting the temperature data inside the closed area includes: Display the coordinates of the thermocouple group based on the thermocouple data, enclose the area including all the acquisition points, and perform interpolation display on the temperature inside the closed area. Thus, the obtained 3D area data can be used to grab 2D data by slicing along any direction for backup.

[0035] In an embodiment of the present invention, the step S10 of reading the thermocouple data arranged on the hearth, performing three-dimensional visualization on the thermocouple data, and extracting the temperature data inside the closed area includes: Display the coordinates of the thermocouple group based on the thermocouple data through 3D visualization software, enclose the area including all the acquisition points, and perform interpolation display on the temperature inside the closed area by using the griddata interpolation algorithm. Thus, the obtained 3D area data can be used to grab 2D data by slicing along any direction for backup.

[0036] Specifically, considering Figure 3The heat transfer process shown. Molten iron and slag drop from the upper part into the hearth, gradually accumulate, and are periodically discharged from the taphole. The molten pool in the hearth can be divided into a dynamic zone and a static zone according to the height of the taphole. In the dynamic zone, due to the flow of molten iron, heat transfer is mainly by thermal convection; while in the static zone, since it is below the taphole and the flow inside is relatively static, heat transfer inside is mainly by thermal conduction. Due to the relative stillness in the static zone, molten iron is likely to solidify at the bottom to form an iron freezing layer. Under the action of chemical erosion and mechanical scouring, the furnace lining is continuously eroded, forming an erosion line. Chemical erosion is mainly in the static zone. Generally, 1150°C is used as the reference line, which is the solidification temperature of molten iron. Below this temperature, due to the formation of the iron freezing layer, the erosion will not continue.

[0037] Based on the above characteristics, the hearth can be modeled as shown Figure 4 in the following way. Its computational domain only considers the static zone and the furnace lining. Heat transfer in this computational domain is mainly by thermal conduction, so it can be regarded as a whole. The influence of the furnace lining properties and the iron freezing layer on heat transfer is mainly reflected in its thermal conductivity. The dynamic zone is simplified as a boundary condition, assuming that there is no temperature gradient inside the dynamic zone, and the temperature at the interface between the dynamic and static zones is the temperature of the molten iron in the hearth. Considering the heat loss during the iron discharging process, it is generally considered that the temperature of the molten iron inside the hearth is 50°C lower than the tapping temperature. Therefore, the temperature of the molten iron in the furnace is considered to be about 1480 - 1520°C. It should be noted that since different thermal conductivities are assigned to the molten iron, the furnace lining, and the iron freezing layer, and this coefficient is a function of position and temperature, there is no need to dynamically adjust the thickness of the furnace lining. The predicted 1150°C isotherm is the erosion line under the current calculation conditions. The biggest difference between the method described in the present invention and the existing "reverse calculation" method is that the existing method takes the furnace lining thickness and one-end temperature as inputs, solves for the other-end temperature, and continuously adjusts the furnace lining thickness to make the solved temperature reach the preset value; while the existing method takes the two-end temperatures as inputs and directly solves for the intermediate temperature isotherm, so the calculation efficiency is much higher.

[0038] S20. Establish a 2D geometric model of the hearth and discretize it to obtain a discretized grid; In the embodiment of the present invention, the step S20 of establishing a 2D geometric model of the hearth and discretizing it includes: According to the on-site monitoring requirements, select a specific 2D cross-section for geometric modeling and discretize it according to the calculation accuracy requirements. The maximum grid scale cannot exceed 1 / 2 of the minimum furnace lining material thickness.

[0039] Establish geometric models under normal and extreme conditions respectively. The three-dimensional asymmetry of the hearth can be considered by calculating different cross-sections.

[0040] Specifically, consider as Figure 5The extreme conditions shown are used to simulate the phenomena of enhanced convective heat transfer at the bottom of the molten pool caused by factors such as the "annular gap" flow of hot metal, natural convection caused by the temperature difference of hot metal, and liquid level fluctuations (such as large discharges) caused by severe fluctuations in the furnace pressure during iron tapping. Since the above phenomena cannot be quantified, they are considered in the most extreme way, that is, the entire interior of the molten pool is regarded as a dynamic region, and the temperature inside is the hot metal temperature. The erosion line predicted under this condition is the maximum erosion line that can be theoretically achieved, and the actual erosion line should be between Figure 4 and Figure 5 the two working conditions.

[0041] S30. Interpolate the internal boundary conditions using the internal temperature data of the closed region to obtain the internal boundary conditions.

[0042] In the embodiment of the present invention, the S30. Interpolate the internal boundary conditions using the internal temperature data of the closed region to obtain the internal boundary conditions, including: Interpolate the internal boundary conditions using the internal temperature data of the closed region to obtain the temperature data of all nodes in the overlapping region, and the temperature data of these nodes will be used as the internal boundary conditions to participate in the subsequent calculation process.

[0043] Specifically, as Figure 6 shown, taking the two-dimensional condition as an example, the closed region formed inside the thermocouple coverage area will partially overlap with the established discretized grid. By using Delaunay triangulation for interpolation, the temperature data of all nodes in the overlapping region can be obtained. And the temperature data of these nodes will be used as the internal boundary conditions to participate in the subsequent calculation process.

[0044] S40. Set the external boundary conditions of the three-dimensional visualization model based on the pre-obtained real-time tapping temperature and the temperature difference between cooling water, and use the set internal boundary conditions. Establish a heat conduction model using the external boundary conditions and the internal boundary conditions.

[0045] S50. Solve the discretized grid using the heat conduction model and integrate the obtained results. Since the calculation conditions in the region between the thermocouple and the external boundary are the same under the two working conditions, the calculation results are also the same, while there will be obvious differences in the temperature of the furnace lining inside the thermocouple. Therefore, when superimposing, directly superimpose the erosion line under the "extreme" condition on the result under the "normal" condition, and thus the temperature distribution of the entire furnace lining and the position of the "erosion zone" can be obtained.

[0046] In the embodiment of the present invention, as Figure 7 shown, the S50. Solve the discretized grid using the heat conduction model and integrate the obtained results, including: Based on the discretized grid, the thermal state of the hearth is solved. The heat conduction process inside the hearth lining satisfies Equation [1]. [1] For simplicity, considering the two-dimensional case, we have: [2] In the formula, k is the thermal conductivity of the material, which is a function of position and temperature. Since different lining materials and frozen iron layers are involved, it cannot be directly taken outside the parentheses. After expanding the above formula, we get: [3] The finite difference method is used to solve the above formula, where the first derivative uses forward difference and the second derivative uses central difference. That is: ; ; ; ; ;

[0047] Among them, i and j respectively represent the nodes at the positions of the i th row and j th column, as shown in Figure 5 . Substituting the above difference form into Equation [3], we have: [4] After sorting, we get: [5] Then the temperature expression at the node of the i-th row and j-th column can be derived as: [6] From Equation [6], we know that solving requires knowing the temperatures of the four adjacent nodes, that is: , , , , as well as the thermal conductivities of the nodes on its own, right, and upper sides, that is , , . Under the given boundary conditions, can be solved through iterative loops.

[0048] The method of the present invention can completely display the temperature field distribution of the furnace lining from the inside to the outside; different from the existing prediction methods, the method of the present invention combines two calculation ideas of "forward prediction" and "reverse prediction", directly regards the thermocouple coverage area as the internal boundary, and the temperature value inside it is directly determined by the measurement data, while the remaining furnace lining areas are calculated based on the heat conduction model. In addition, since the calculation domain is extended and the heat conduction of the static area at the bottom of the molten pool is considered, the isotherm of the molten iron condensation temperature can be obtained without dynamically adjusting the furnace lining thickness, and the calculation efficiency is improved compared with the existing methods; by superimposing the results of two working conditions of "normal" and "limit", the concept of "erosion zone" is used to replace the "erosion line" in the prior art to predict the erosion state of the furnace lining, reducing the influence of the dynamic heat transfer process in the hearth on the calculation results and enhancing the confidence level of the prediction results. The present invention takes the temperatures at both ends as inputs and directly solves the intermediate temperature isotherm, so the calculation efficiency is much higher.

[0049] It should be understood that although the above is described in a certain order, these steps are not necessarily executed in the above order. Unless there is a clear description in this article, there is no strict order limit for the execution of these steps, and these steps can be executed in other orders. Moreover, a part of the steps in this embodiment may include multiple steps or multiple stages. These steps or stages are not necessarily executed at the same time, but can be executed at different times. The execution order of these steps or stages is not necessarily sequential, but can be executed alternately or alternately with at least a part of other steps or steps or stages in other steps.

[0050] It should be understood that, as used in this article, unless the context clearly supports an exception, the singular form "a" is also intended to include the plural form. It should also be understood that "and / or" as used in this article refers to any and all possible combinations of one or more of the related listed items. The serial numbers of the disclosed embodiments of the present invention are only for description and do not represent the advantages and disadvantages of the embodiments.

[0051] Those of ordinary skill in the art should understand that: the discussion of any of the above embodiments is only exemplary and is not intended to imply that the scope of the disclosure of the embodiments of the present invention (including the claims) is limited to these examples; under the idea of the embodiments of the present invention, the technical features in the above embodiments or different embodiments can also be combined, and there are many other variations in different aspects of the embodiments of the present invention as above, and they are not provided in detail for the sake of brevity. Therefore, any omission, modification, equivalent replacement, improvement, etc. made within the spirit and principle of the embodiments of the present invention shall be included in the protection scope of the embodiments of the present invention.

Claims

1. An online prediction method for furnace lining erosion line, characterized in that: The method includes: Reading thermocouple data arranged on the furnace, performing three-dimensional visualization on the thermocouple data, and extracting temperature data inside the closed area; Establish a 2D geometric model of the furnace and discretize it to obtain a discretized grid; The internal boundary conditions are interpolated using the internal temperature data of the closed area to obtain the internal boundary conditions; The external boundary conditions of the three-dimensional visualization model are set based on the pre-acquired real-time iron tapping temperature and cooling water temperature difference, and the internal boundary conditions are set, and the heat conduction model is established using the external boundary conditions and the internal boundary conditions; The discretized grid is solved using a heat conduction model and the results are integrated.

2. The method for online prediction of furnace lining erosion line according to claim 1, characterized in that: Thermocouple data includes location and temperature information.

3. The method for online prediction of furnace lining erosion line according to claim 1, characterized in that: The step of reading the thermocouple data arranged on the furnace, performing three-dimensional visualization on the thermocouple data, and extracting the temperature data inside the closed area includes: Thermocouple group coordinates are displayed based on the thermocouple data.

4. The method for online prediction of furnace lining erosion line according to claim 3, characterized in that: Also includes: The area including all the acquisition points is closed, and the temperature in the closed area is interpolated and displayed.

5. The method for online prediction of furnace lining erosion line according to claim 4, characterized in that: The step of reading the thermocouple data arranged on the furnace, performing three-dimensional visualization on the thermocouple data, and extracting the temperature data inside the closed area includes: The coordinates of the thermocouple group are displayed by 3D visualization software based on the thermocouple data.

6. The method for online prediction of furnace lining erosion line according to claim 5, characterized in that: Also includes: The area including all the collected points is closed, and the temperature in the closed area is interpolated and displayed using the griddata interpolation algorithm.

7. The method for online prediction of furnace lining erosion line according to claim 4, characterized in that: The step of establishing and discretizing the 2D geometric model of the furnace includes: Depending on the needs of on-site monitoring, select a specific 2D section for geometric modeling and discretize it according to the calculation accuracy requirements. The maximum grid size cannot exceed 1 / 2 of the minimum furnace lining material thickness; The geometric models under normal and limit conditions are established respectively.

8. The method for online prediction of furnace lining erosion line according to claim 1, characterized in that: The method of interpolating the internal boundary conditions using the internal temperature data of the closed area to obtain the internal boundary conditions includes: The internal boundary conditions are interpolated using the internal temperature data of the closed area to obtain the temperature data of all nodes in the overlapping area.

9. The method for online prediction of furnace lining erosion line according to claim 1, characterized in that: The heat conduction model is used to solve the discretized grid and the obtained results are integrated, including: Based on the discretized grid, the thermal state of the furnace is solved.

10. The method for online prediction of furnace lining erosion line according to claim 9, characterized in that: Solve the furnace thermal state, including calculating node temperatures, The node temperature calculation formula is as follows: In the formula, is the node temperature, i and j Respectively represent i OK, j The node where the column is located, k is the thermal conductivity of the material.

Images