CN116362085A - Hearth lining erosion morphology identification method based on cooling wall heat flow intensity - Google Patents
Hearth lining erosion morphology identification method based on cooling wall heat flow intensity Download PDFInfo
- Publication number
- CN116362085A CN116362085A CN202310332273.5A CN202310332273A CN116362085A CN 116362085 A CN116362085 A CN 116362085A CN 202310332273 A CN202310332273 A CN 202310332273A CN 116362085 A CN116362085 A CN 116362085A
- Authority
- CN
- China
- Prior art keywords
- erosion
- cooling wall
- point
- heat flow
- boundary
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000003628 erosive effect Effects 0.000 title claims abstract description 150
- 238000001816 cooling Methods 0.000 title claims abstract description 123
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000004364 calculation method Methods 0.000 claims abstract description 60
- 230000007797 corrosion Effects 0.000 claims abstract description 8
- 238000005260 corrosion Methods 0.000 claims abstract description 8
- 238000012795 verification Methods 0.000 claims abstract description 5
- 230000004907 flux Effects 0.000 claims description 8
- 238000004088 simulation Methods 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 7
- 238000001514 detection method Methods 0.000 abstract description 2
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 14
- 238000009529 body temperature measurement Methods 0.000 description 9
- 229910052742 iron Inorganic materials 0.000 description 7
- 238000012546 transfer Methods 0.000 description 7
- 238000003745 diagnosis Methods 0.000 description 4
- 239000000498 cooling water Substances 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000003211 malignant effect Effects 0.000 description 3
- 230000006378 damage Effects 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000007711 solidification Methods 0.000 description 2
- 230000008023 solidification Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011449 brick Substances 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 238000012625 in-situ measurement Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000009991 scouring Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000008646 thermal stress Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- C—CHEMISTRY; METALLURGY
- C21—METALLURGY OF IRON
- C21B—MANUFACTURE OF IRON OR STEEL
- C21B7/00—Blast furnaces
- C21B7/04—Blast furnaces with special refractories
- C21B7/06—Linings for furnaces
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Manufacturing & Machinery (AREA)
- Computer Hardware Design (AREA)
- Materials Engineering (AREA)
- Metallurgy (AREA)
- Organic Chemistry (AREA)
- Waste-Gas Treatment And Other Accessory Devices For Furnaces (AREA)
Abstract
The invention provides a hearth lining erosion morphology identification method based on cooling wall heat flow intensity, and relates to the technical field of blast furnace hearth erosion detection. According to the method, a reference calculation model and an erosion calculation model are established by using three-dimensional software, and different erosion positions are compared with the reference calculation model to obtain the influence rule of the erosion position change on the heat flow intensity of the cooling wall; calculating the percentage difference delta of the heat flow intensity of the corresponding cooling wall in each erosion calculation model and the reference calculation model d Calibrating delta d The relation between the erosion positions to obtain the approximate position range of the actual concave erosion, further obtain the accurate position, establish the corresponding initial erosion boundary and the erosion boundary movement searchThe method takes the heat flow intensity of the cooling wall as the verification parameter of the corrosion boundary, and reversely calculates the actual corrosion boundary and the minimum residual thickness of the hearth lining.
Description
Technical Field
The invention relates to the technical field of blast furnace hearth erosion detection, in particular to a hearth lining erosion morphology identification method based on cooling wall heat flow intensity.
Background
Blast furnaces are used as large-scale iron-making equipment, wherein the hearth at the lower part is a frequent part of safety accidents. The refractory lining in the hearth is directly contacted with high-temperature molten iron, and irreversible erosion occurs under the comprehensive effects of hot iron scouring, chemical erosion, thermal stress and the like. Because of complex erosion reasons, the erosion morphology is often irregular, local dent often occurs, and when serious, hearth burning-through accidents can be caused, so that high-temperature molten iron in the hearth leaks to cause equipment damage, even severe explosion and other malignant consequences. In addition, the hearth is in a closed state for a long time and is immersed by high-temperature molten iron, so that the erosion state of the hearth cannot be directly detected. Therefore, judging whether the hearth lining is corroded or not and accurately identifying the corrosion position and the corrosion degree are important technologies for improving the safety and the durability of the hearth.
In the prior art, the calculation of the erosion boundary of the blast furnace hearth is generally to perform the thermal conductivity calculation on thermocouple temperature measurement data pre-embedded in the lining of the hearth to obtain the lining erosion boundary and the hot iron solidification boundary.
The invention patent is a two-dimensional rapid calculation method for erosion boundary of longitudinal section of blast furnace hearth, and the application publication number is: CN 110765623a, solving the erosion boundary position of the hearth by using the detected temperature value of the thermocouple disposed in the hearth, passing through the outer layer thermocouple P i,out The temperature of the outer boundary layer is determined by the temperature and the coordinates of the outer boundary layer, the initial inner wall erosion line is calculated preliminarily by utilizing a linear interpolation method, and then the initial inner wall erosion line is connected with an actual inner layer thermocouple P i,in The temperature values are compared and corrected for a plurality of times until the error meets the range requirement.
The invention patent 'a finite element method for calculating erosion boundary of longitudinal section of hearth', application publication number: CN 114996997A, by combining the temperature values detected by thermocouple groups arranged in the hearth carbon bricks with the particle swarm algorithm, implements calculation of the position of the hearth erosion boundary, and also considers the curvature constraint of the erosion boundary to obtain a smoother erosion boundary. The invention patent is a method for detecting the erosion condition of a blast furnace hearth, which has the application publication number: CN 109929955A, a thermocouple temperature measurement historical data database and a furnace wall internal temperature field calculation model are established, and a series of thermocouple historical temperature measurement data are calculated and analyzed through an objective function, so as to obtain the erosion conditions of the hearth in different periods.
Feng proposes a hearth erosion monitoring model which can carry out inversion solution on unknown boundary problems of a hearth and monitor the residual thickness of a hearth lining based on thermocouple feedback temperature in the hearth to obtain a hearth erosion boundary.
The brannbackup creates a mathematical model containing erosion boundaries based on embedded thermocouple data in the hearth. And the erosion boundary of the model is adjusted, so that the temperature distribution of the mathematical model approaches to reality, and the liner erosion morphology diagnosis is realized. Based on two-dimensional heat transfer theory such as Zagaria, solving of erosion boundary of hearth lining is realized based on thermocouple temperature measurement data, and a solution for interpolating new thermocouples is provided for solving the problem of thermocouple damage.
In the above technology, the calculation of the erosion boundary of the blast furnace hearth is generally that the lining erosion boundary and the hot iron solidification boundary are obtained after the thermal conductivity calculation is carried out on thermocouple temperature measurement data pre-embedded in the lining of the hearth. However, the embedded thermocouple in the hearth lining often fails in a considerable proportion due to long-term working at high temperature, and the inserted thermocouple has danger, and particularly, the iron leakage accident is often caused by the thinner residual thickness of the lining.
Secondly, the number of thermocouples pre-buried in the prior middle and small blast furnace lining is generally small, so that a temperature measurement blind area exists in the hearth lining, the erosion morphology in the blind area cannot be accurately diagnosed, and the erosion morphology is important to the safety of the hearth structure, especially in the situation of concave erosion.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a hearth lining erosion morphology identification method based on the heat flow intensity of the cooling wall, and provides a moving boundary search method based on the heat flow intensity of the cooling wall on the basis, which realizes the identification of the lining erosion morphology of the thermocouple temperature measurement blind area and has important significance for the hearth structure safety assessment work. Solves the difficult problem of erosion diagnosis caused by the failure phenomenon of the thermocouple or the temperature measurement blind area of the thermocouple in the hearth, avoids the occurrence of malignant accidents of hearth burning, and ensures the safety and durability of the blast furnace.
A hearth lining erosion morphology identification method based on cooling wall heat flow intensity comprises the following steps:
step 1: based on the blast furnace hearth structure and the cooling wall distribution, a reference calculation model and an erosion calculation model are established by using three-dimensional simulation software, and the influence rule of erosion position change on the cooling wall heat flow intensity is obtained by comparing the calculation results of different erosion positions and the reference calculation model;
the reference calculation model is a calculation model when no erosion occurs;
the erosion calculation model is a calculation model with local concave erosion at five typical positions A1-A5 in the longitudinal direction respectively; five positions A1-A5 are positioned on the same reference line and are arranged from high to low according to the sequence, wherein A1 and A2 are positioned at the left upper part of the cooling wall, A3 is positioned at the juncture of the cooling wall, and A4 and A5 are positioned at the right lower part of the cooling wall.
Step 2: calculating the percentage difference delta of the heat flow intensity of the corresponding cooling wall in each erosion calculation model and the reference calculation model d Calibrating the relation between the percentage difference and each erosion position;
wherein the percentage difference delta d The calculation formula of (2) is as follows:
wherein q is iE To erode the cooling wall i heat flow intensity value, q in the calculation model iS The heat flow intensity value of the cooling wall i in the model is calculated as a reference.
Step 3: equivalent average heat flow intensity values of all cooling walls in the same section are used as a reference calculation model, and the erosion position range is judged according to the influence rule in the step 1 and the actually measured heat flow intensity change of the cooling walls;
step 4: and (2) after determining the erosion position range, calculating the heat flow intensity percentage difference between the cooling wall where the erosion position range is positioned and the reference model, and judging the accurate position of the actual pit erosion according to the calibration result in the step (2).
Step 5: and (3) establishing a corresponding initial erosion boundary and an erosion boundary moving search method by combining the accurate position of the actual concave erosion in the step (4), taking the heat flow intensity of the cooling wall as the verification parameter of the required erosion boundary, and obtaining the actual erosion boundary of the hearth lining by using the erosion boundary moving search method.
According to the erosion boundary moving search method, 7 boundary control points are sequentially connected through a broken line or spline curve to construct an initial erosion boundary, and 3 control moving points are set to be strongly related to the heat flow intensity of the three-section cooling wall respectively. If the simulation value of the heat flow intensity of the cooling wall is smaller than the actual measurement value, the corresponding boundary control point moves to the outer side of the hearth, otherwise, moves to the inner side, and stops searching when the convergence judgment type boundary is met, and the obtained boundary is the actual erosion boundary of the hearth.
The 7 boundary control points are boundary control points 1-7 respectively, and point 3, point 4 and point 5 are control moving points which are strongly related to the heat flow intensity of the three-section cooling wall; determining a point 1 by calculating the temperature of a thermocouple at the bottom of the furnace; the passing point 1 is a horizontal straight line L 1 An extension line passing through the inner furnace corner point as a side wall boundary is connected with a straight line L 1 Crossing at point 2'; taking the midpoint of the point 1 and the point 2' as the point 2; connecting the passing inner furnace corner point with the midpoint of the first section cooling wall and a straight line L 1 Crossing at point 3; the temperature line position of the liner at 1150 ℃ is calculated according to the cooling conditions and the heat flow intensity of the second and third sections of cooling walls and is marked as a straight line L 2 And L 3 Taking a straight line L 3 Is point 6. If erosion occurs on the second section cooling wall, making a horizontal line determination point 4 at the center position of the second section cooling wall; making horizontal lines at the junction positions of the second and third sections of cooling walls and respectively connecting with the straight lines L 2 And L 3 Intersecting point 5 'and point 5 ", taking the midpoint between point 5' and point 5" as point 5; making a straight line with the passing point 5 and the point 6, and intersecting with the top at a point 7; point 6 and point 7 move synchronously with point 5And (5) moving. If corrosion occurs at the junction of the second cooling wall and the third cooling wall, taking the midpoint of the straight line L2 as a point 4; making a horizontal line L at the junction position of the second and third sections of cooling walls 4 The method comprises the steps of carrying out a first treatment on the surface of the The passing point 6 is a straight line forming an angle theta with the horizontal direction and intersects with the top of the hearth at a point 7 and L 4 Crossing at point 5; the points 6 and 7 move synchronously with the point 5, and the three points are always collinear and keep an angle theta with the horizontal direction.
The beneficial effects of adopting above-mentioned technical scheme to produce lie in:
the invention provides a hearth lining erosion morphology identification method based on cooling wall heat flow intensity, which can judge local erosion occurrence positions according to cooling wall heat flow intensity changes, and reversely calculate erosion boundaries based on an erosion boundary moving search method to obtain the minimum residual thickness of the hearth lining, so that the recognition of the erosion morphology of the lining in a thermocouple temperature measurement blind area is realized, further, the hearth lining erosion diagnosis method is enriched, the occurrence of malignant accidents caused by hearth burning can be effectively avoided, and the method has important significance in prolonging the service life of a blast furnace and improving the safety and durability of the hearth.
Drawings
FIG. 1 is a flowchart of a method for identifying erosion morphology of a hearth lining according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a stave position distribution and a pit erosion position distribution according to an embodiment of the present invention;
FIG. 3 is a graph showing the results of the differential calibration of the heat flow intensity percentages of the cooling wall for different calculation models provided by the embodiment of the invention;
FIG. 4 is a schematic illustration of an initial boundary layout for erosion occurring on a second stage stave according to an embodiment of the present invention;
FIG. 5 is a schematic diagram illustrating an initial boundary planning of the boundary between two or three cooling walls according to the embodiment of the present invention;
FIG. 6 is a schematic diagram of calculated erosion boundaries provided by an embodiment of the present invention.
Detailed Description
The following describes in further detail the embodiments of the present invention with reference to the drawings and examples. The following examples are illustrative of the invention and are not intended to limit the scope of the invention.
A hearth lining erosion morphology identification method based on cooling wall heat flow intensity, as shown in figure 1, comprises the following steps:
step 1: based on the blast furnace hearth structure and the cooling wall distribution, a reference calculation model and an erosion calculation model are established by using three-dimensional simulation software, and the influence rule of erosion position change on the cooling wall heat flow intensity is obtained by comparing the calculation results of different erosion positions and the reference calculation model;
the reference calculation model is a calculation model when no erosion occurs;
the erosion calculation model is a calculation model with local concave erosion at five typical positions (A1-A5) in the longitudinal direction respectively;
as shown in FIG. 2, the hearth part has three cooling walls, the cooling walls are staggered, each cooling wall consists of 32 cooling walls, and the cooling walls are uniformly distributed along the circumferential direction. A reference calculation model in which no erosion occurs and an erosion calculation model in which recess erosion occurs at five positions (A1 to A5) in the longitudinal direction in fig. 2, respectively, are established. Five positions A1-A5 are positioned on the same reference line and are arranged from high to low according to the sequence, wherein A1 and A2 are positioned at the left upper part of the 2-2 cooling wall, A3 is positioned at the junction of the 2-2 cooling wall and the 3-1 cooling wall, and A4 and A5 are positioned at the right lower part of the 3-1 cooling wall.
The heat transfer process of two calculation models is simulated based on a finite element method, the average heat flow intensity value of each cooling wall in each calculation model is extracted, the number of the cooling wall is shown in a table 1, the number of the cooling wall is shown in fig. 2, and the erosion position and the heat flow intensity change rule are compared as follows:
(1) If the heat flux intensity of a stave is significantly increased while the heat flux intensity of its neighboring staves is increased less, localized dishing erosion occurs in the region corresponding to the stave.
(2) If the heat flow intensity of two or three adjacent staves is increased at the same time and there is no abnormality in other staves adjacent to them, localized pitting corrosion is present at the boundary positions of these staves.
(3) If the increase in the heat flux intensity of the staves is different, the local pit erosion position is biased toward the staves with the larger increase.
TABLE 1 average heat flow intensity simulation values (W/m for staves 2 )
Step 2: calculating the percentage difference delta of the heat flow intensity of the corresponding cooling wall in each erosion calculation model and the reference calculation model d The relationship between the calibrated percentage difference and each erosion location is shown in FIG. 3;
wherein the percentage difference delta d The calculation formula of (2) is as follows:
wherein q is iE To erode the cooling wall i heat flow intensity value, q in the calculation model iS The heat flow intensity value of the cooling wall i in the model is calculated as a reference.
Calculating the percentage difference between the cooling wall of the erosion position range and the reference model, wherein the average heat flow intensity value of all cooling walls in the same section in actual measurement is equivalent to the reference model to calculate the percentage difference;
after determining the approximate range of erosion, the heat flow intensity percentage difference of all the cooling walls around the approximate range is calculated, the calibration result of the heat flow intensity percentage difference of each cooling wall in the models A1-A5 is shown in figure 3, and the rule of the calibration result is as follows:
(1) If the erosion site occurs in a single stave, the heat flux intensity of the stave is maximized.
(2) The increase in heat flux intensity increases gradually as the localized pit erodes toward the center of the stave.
(3) If erosion occurs at the juncture of two staves, heat is commonly carried by the two staves, resulting in a percent increase in heat flow intensity per stave of less than 10%.
Step 3: equivalent average heat flow intensity values of all cooling walls in the same section are used as a reference calculation model, and the erosion position range is judged according to the influence rule in the step 1 and the actually measured heat flow intensity change of the cooling walls;
in the actual test, the change of the cooling wall accords with the rule (1), and the partial concave erosion can be primarily judged to be approximately generated on the 2-2 cooling wall.
Step 4: and (2) after determining the erosion position range, calculating the heat flow intensity percentage difference between the cooling wall where the erosion position range is positioned and the reference model, and judging the accurate position of the actual pit erosion according to the calibration result in the step (2).
The stave change conforms to rule (1) in step 1, and a preliminary determination is made that localized pitting corrosion generally occurs on the 2-2 stave. And (3) equivalent average heat flow intensity values of all cooling walls in the same section in actual measurement are used as a reference model to calculate a percentage difference, the calculated result is compared with a calibration result, the percentage difference of the heat flow intensity of 2-2 is far greater than that of the other cooling walls, the heat flow intensity of the cooling walls of 1-1 and 1-2 is greater than that of 3-1 and 3-2, the erosion position is shown to occur on the cooling walls of the second section and is closer to the cooling walls of the first section, and the A1 position of the erosion in the cooling walls of the second section is determined by comparison.
Step 5: and (3) establishing a corresponding initial erosion boundary and an erosion boundary moving search method by combining the accurate position of the actual concave erosion in the step (4), taking the heat flow intensity of the cooling wall as the verification parameter of the required erosion boundary, and obtaining the actual erosion boundary of the hearth lining by using the erosion boundary moving search method.
The initial erosion boundary is obtained by sequentially connecting boundary control points 1-7 through broken lines or spline curves, wherein points 3, 4 and 5 are corresponding control moving points which are strongly related to the heat flow intensity of the first, second and third sections of cooling walls and can move along the directions shown in fig. 4 and 5; the furnace bottom center boundary control point 1 is determined by calculating the heat flow intensity of a furnace bottom water pipe according to one-dimensional heat transfer; the passing point 1 is a horizontal straight line L 1 The corner point of the passing inner furnace is taken as an extension line of the inner boundary of the original side wall and is connected with a straight line L 1 Intersection point 2'; taking the midpoint of the point 1 and the point 2' as a boundary control point 2; the point 1 and the point 2 are far away from the cooling wall, and fixed-point treatment is carried out; connecting the passing inner furnace corner point with the midpoint of the first section cooling wall and a straight line L 1 Crossing the boundary control point 3; respectively calculating the 1150 ℃ isotherm position of the lining according to the cooling conditions and the heat flow intensity of the second and third sections of cooling walls, and marking as straightLine L 2 And L 3 Taking a straight line L 3 Is the boundary control point 6.
If erosion occurs in the second-stage cooling wall, as shown in the initial boundary planning in fig. 4, determining a boundary control point 4 by making a horizontal line at the center of the second-stage cooling wall; making horizontal lines at the junction positions of the second and third sections of cooling walls and respectively connecting with the straight lines L 2 And L 3 Intersecting with the point 5 'and the point 5', taking the midpoint between the point 5 'and the point 5' as a boundary control point 5; the passing point 5 and the point 6 are made into a straight line and are intersected with the top part at a boundary control point 7; points 6 and 7 move synchronously with point 5.
If erosion occurs at the junction of the second and third cooling walls, as shown in the initial boundary plan in FIG. 5, a straight line L is taken 2 Is the boundary control point 4; making a horizontal line L at the junction position of the second and third sections of cooling walls 4 The method comprises the steps of carrying out a first treatment on the surface of the The passing point 6 is a straight line forming an angle theta with the horizontal direction, and is intersected with the top of the hearth at the boundary control point 7 and a horizontal line L 4 Intersecting with the boundary control point 5; the points 6 and 7 move synchronously with the point 5, the three points are always collinear and keep an angle theta with the horizontal direction, the size of the angle theta is determined according to specific conditions, and theoretically, the smaller the theta is, the smaller the residual thickness of the lining at the position of the point 5 is.
After determining the erosion occurrence position and the initial search boundary, the heat flow intensity q is measured by the cooling wall of the second section and the third section 2T 、q 3T Calculating an isothermal line L at 1150 ℃ according to a one-dimensional flat wall heat transfer theory 2 And L 3 The boundary control point moves according to a preset movement strategy, if the simulation value of the heat flow intensity of the cooling wall is smaller than the actual measurement value, the corresponding boundary control point moves to the outside of the hearth, otherwise, the boundary control point moves to the inside. When the relative error delta i of the heat flow intensity simulation value and the actual measurement value of the three-section cooling wall belongs to a reasonable interval e, the erosion boundary can be obtained, namely, the conditions are satisfied:
|δ i |≤e(i＝1,2,3)
the calculation accuracy required for different blast furnaces varies, and e=1 to 2% is generally desirable depending on the circumstances.
Measured heat flow intensity q of cooling wall iT The calculation formula of (2) is as follows:
wherein, c p J/(kg.K) is the specific heat capacity of cooling water; m is M s For cooling water flow, kg/s; delta T is the temperature difference of water and DEG C; s is cooling area, m 2 。
The calculation formula of the relative error delta i between the cooling wall heat flow intensity simulation value and the measured value is as follows:
wherein q is iC (i=1, 2, 3) and q iT (i=1, 2, 3) are respectively the simulation value and the actual measurement value of the heat flow intensity of the cooling wall, W/m 2 。
After determining that the erosion occurs at the A1 position, selecting an initial erosion boundary in fig. 4, and performing heat transfer process simulation calculation on the initial erosion boundary by using a finite element method. In the embodiment, the convection heat exchange condition of the inner wall surface of the cooling water pipe is equivalent to the inner lining cold surface, and the equivalent convection heat exchange coefficient h is calculated e 67.7W/(m) 2 K) equivalent convective heat transfer coefficient h of bottom surface of furnace bottom lining d 40W/(m) 2 K) the coefficient of thermal conductivity of the lining is 12.0W/(mK), the average temperature T of the air f 30℃was taken. In situ measurement q for actual erosion diagnosis 1T 、q 2T 、q 3T As a verification parameter for solving the erosion boundary, the field test value is shown in FIG. 6, q 1T ＝6130W/m 2 ，q 1T ＝26457W/m 2 ，q 3T ＝18039W/m 2 。
Calculating an isothermal line L at 1150 ℃ according to a one-dimensional flat wall heat transfer theory by the cooling conditions and the heat flow intensity of the second cooling wall and the third cooling wall 2 And L is equal to 3 331mm and 528mm from the outer edge of the hearth lining, respectively. The boundary movement strategy is that the point 6 and the point 7 move synchronously with the point 5, and the initial search step size is set to be 40mm. The iterative search process is shown in table 2, and after the 6 th step search is completed, |δi| (i=1, 2, 3) is less than 1%, the convergence judgment formula is satisfied, and the graph shown in fig. 6 is obtainedThe erosion boundary is shown. The radius of the concave position boundary control point 4 is 3930mm, the hearth radius is 4200mm, and the remaining thickness of the liner is 270mm.
Table 2 iterative search procedure
The foregoing description is only of the preferred embodiments of the present disclosure and description of the principles of the technology being employed. It will be appreciated by those skilled in the art that the scope of the invention in the embodiments of the present disclosure is not limited to the specific combination of the above technical features, but encompasses other technical features formed by any combination of the above technical features or their equivalents without departing from the spirit of the invention. Such as the above-described features, are mutually substituted with (but not limited to) the features having similar functions disclosed in the embodiments of the present disclosure.
Claims (7)
1. A hearth lining erosion morphology identification method based on cooling wall heat flow intensity is characterized by comprising the following steps:
step 1: based on the blast furnace hearth structure and the cooling wall distribution, a reference calculation model and an erosion calculation model are established by using three-dimensional simulation software, and the influence rule of erosion position change on the cooling wall heat flow intensity is obtained by comparing the calculation results of different erosion positions and the reference calculation model;
step 2: calculating the percentage difference delta of the heat flow intensity of the corresponding cooling wall in each erosion calculation model and the reference calculation model d Calibrating the relation between the percentage difference and each erosion position;
step 3: equivalent average heat flow intensity values of all cooling walls in the same section are used as a reference calculation model, and the erosion position range is judged according to the influence rule in the step 1 and the actually measured heat flow intensity change of the cooling walls;
step 4: after determining the erosion position range, calculating the heat flow intensity percentage difference between the cooling wall where the erosion position range is located and the reference model, and judging the accurate position of the actual concave erosion according to the calibration result in the step 2;
step 5: and (3) establishing a corresponding initial erosion boundary and an erosion boundary moving search method by combining the accurate position of the actual concave erosion in the step (4), taking the heat flow intensity of the cooling wall as the verification parameter of the required erosion boundary, and obtaining the actual erosion boundary of the hearth lining by using the erosion boundary moving search method.
2. The method for recognizing the erosion morphology of the hearth lining based on the heat flow intensity of the cooling wall according to claim 1, wherein the reference calculation model in the step 1 is a calculation model when no erosion occurs.
3. The method for recognizing the erosion morphology of the hearth lining based on the heat flow intensity of the cooling wall according to claim 1, wherein the erosion calculation model in the step 1 is a calculation model in which partial dent erosion occurs in five typical positions A1 to A5 in the longitudinal direction respectively; five positions A1-A5 are positioned on the same reference line and are arranged from high to low according to the sequence, wherein A1 and A2 are positioned at the left upper part of the cooling wall, A3 is positioned at the juncture of the cooling wall, and A4 and A5 are positioned at the right lower part of the cooling wall.
4. The method for recognizing erosion morphology of hearth lining based on heat flow intensity of cooling wall according to claim 1, wherein the percentage difference δ in step 2 d The calculation formula of (2) is as follows:
wherein q is iE To erode the cooling wall i heat flow intensity value, q in the calculation model iS The heat flow intensity value of the cooling wall i in the model is calculated as a reference.
5. The method for recognizing the erosion morphology of the hearth lining based on the heat flux intensity of the cooling wall according to claim 1, wherein in the erosion boundary moving search method in the step 5, 7 boundary control points are sequentially connected by a broken line or spline curve to construct an initial erosion boundary, and 3 control moving points are set to be strongly related to the heat flux intensity of the three-section cooling wall respectively.
6. The hearth lining erosion morphology identification method based on the cooling wall heat flux intensity according to claim 5, wherein the searching method specifically comprises the following steps: if the simulation value of the heat flow intensity of the cooling wall is smaller than the actual measurement value, the corresponding boundary control point moves to the outer side of the hearth, otherwise, moves to the inner side, and stops searching when the convergence judgment type boundary is met, and the obtained boundary is the actual erosion boundary of the hearth.
7. The method for recognizing the erosion morphology of the hearth lining based on the heat flow intensity of the cooling wall according to claim 4, wherein the 7 boundary control points are boundary control points 1-7 respectively, and the points 3, 4 and 5 are control moving points which are strongly related to the heat flow intensity of the three sections of cooling wall; determining a point 1 by calculating the temperature of a thermocouple at the bottom of the furnace; the passing point 1 is a horizontal straight line L 1 An extension line passing through the inner furnace corner point as a side wall boundary is connected with a straight line L 1 Crossing at point 2'; taking the midpoint of the point 1 and the point 2' as the point 2; connecting the passing inner furnace corner point with the midpoint of the first section cooling wall and a straight line L 1 Crossing at point 3; the temperature line position of the liner at 1150 ℃ is calculated according to the cooling conditions and the heat flow intensity of the second and third sections of cooling walls and is marked as a straight line L 2 And L 3 Taking a straight line L 3 Is point 6; if erosion occurs on the second section cooling wall, making a horizontal line determination point 4 at the center position of the second section cooling wall; making horizontal lines at the junction positions of the second and third sections of cooling walls and respectively connecting with the straight lines L 2 And L 3 Intersecting point 5 'and point 5 ", taking the midpoint between point 5' and point 5" as point 5; making a straight line with the passing point 5 and the point 6, and intersecting with the top at a point 7; point 6 and point 7 move synchronously with point 5; if corrosion occurs at the junction of the second cooling wall and the third cooling wall, taking the midpoint of the straight line L2 as a point 4; making a horizontal line L at the junction position of the second and third sections of cooling walls 4 The method comprises the steps of carrying out a first treatment on the surface of the The passing point 6 is a straight line forming an angle theta with the horizontal direction and intersects with the top of the hearth at a point 7 and L 4 Intersection point5, a step of; the points 6 and 7 move synchronously with the point 5, and the three points are always collinear and keep an angle theta with the horizontal direction.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310332273.5A CN116362085B (en) | 2023-03-31 | 2023-03-31 | Hearth lining erosion morphology identification method based on cooling wall heat flow intensity |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310332273.5A CN116362085B (en) | 2023-03-31 | 2023-03-31 | Hearth lining erosion morphology identification method based on cooling wall heat flow intensity |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116362085A true CN116362085A (en) | 2023-06-30 |
| CN116362085B CN116362085B (en) | 2024-01-30 |
Family
ID=86935905
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310332273.5A Active CN116362085B (en) | 2023-03-31 | 2023-03-31 | Hearth lining erosion morphology identification method based on cooling wall heat flow intensity |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116362085B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3585605A (en) * | 1968-07-04 | 1971-06-15 | Ibm | Associative memory data processor |
| CN101812559A (en) * | 2009-12-28 | 2010-08-25 | 中冶南方工程技术有限公司 | Method for analyzing and monitoring erosion of furnace lining of blast furnace |
| CN103866061A (en) * | 2012-12-18 | 2014-06-18 | 宝山钢铁股份有限公司 | Method for monitoring erosion deviation of blast furnace hearth |
| CN110765623A (en) * | 2019-10-29 | 2020-02-07 | 南京航空航天大学 | Two-dimensional rapid calculation method for erosion boundary of longitudinal section of blast furnace hearth |
| CN113139353A (en) * | 2021-05-11 | 2021-07-20 | 东北大学 | Dynamic calculation and online monitoring early warning analysis method for steam pipe network |
| CN113971359A (en) * | 2021-10-28 | 2022-01-25 | 中冶南方工程技术有限公司 | Method for identifying hearth erosion, electronic equipment and storage medium |
| US20220072593A1 (en) * | 2019-01-17 | 2022-03-10 | Jfe Steel Corporation | Production specification determination method, production method, and production specification determination apparatus for metal material |
-
2023
- 2023-03-31 CN CN202310332273.5A patent/CN116362085B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3585605A (en) * | 1968-07-04 | 1971-06-15 | Ibm | Associative memory data processor |
| CN101812559A (en) * | 2009-12-28 | 2010-08-25 | 中冶南方工程技术有限公司 | Method for analyzing and monitoring erosion of furnace lining of blast furnace |
| CN103866061A (en) * | 2012-12-18 | 2014-06-18 | 宝山钢铁股份有限公司 | Method for monitoring erosion deviation of blast furnace hearth |
| US20220072593A1 (en) * | 2019-01-17 | 2022-03-10 | Jfe Steel Corporation | Production specification determination method, production method, and production specification determination apparatus for metal material |
| CN110765623A (en) * | 2019-10-29 | 2020-02-07 | 南京航空航天大学 | Two-dimensional rapid calculation method for erosion boundary of longitudinal section of blast furnace hearth |
| CN113139353A (en) * | 2021-05-11 | 2021-07-20 | 东北大学 | Dynamic calculation and online monitoring early warning analysis method for steam pipe network |
| CN113971359A (en) * | 2021-10-28 | 2022-01-25 | 中冶南方工程技术有限公司 | Method for identifying hearth erosion, electronic equipment and storage medium |
Non-Patent Citations (5)
| Title |
|---|
| SHILIN等: "Numerical simulation of Erosion boundary identification of blast furnace lining", JOURNAL OF IRON AND STEEL RESEARCH, no. 04, pages 4 - 8 * |
| YANG LI等: "Numerical Study on the Relationship Between the Localized Depression Erosion of a Commercial Blast Furnace Hearth Lining and the Heat Flux of Cooling Staves", 《IEEE ACCESS》, vol. 7, pages 60984 - 60994, XP011726491, DOI: 10.1109/ACCESS.2019.2915915 * |
| 李强;冯明霞;储文;邹宗树;: "基于边界移动法的高炉炉缸侵蚀监测模型", 东北大学学报(自然科学版), no. 01, pages 60 - 65 * |
| 石琳;程素森;左海滨;: "高炉炉衬侵蚀边界识别的数值模拟", 钢铁研究学报, no. 04, pages 4 - 8 * |
| 陈良玉;李玉;王子金;张明;: "传热边界逆解在高炉炉缸侵蚀诊断中的应用", 东北大学学报(自然科学版), no. 08, pages 73 - 76 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116362085B (en) | 2024-01-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103866061B (en) | A kind of blast furnace crucible corrosion deviation monitoring method | |
| CN101812559B (en) | Method for analyzing and monitoring erosion of furnace lining of blast furnace | |
| CN101319256A (en) | Intelligent monitoring method for cooling wall of blast furnace | |
| CN107622154A (en) | The determination method of blast furnace heat load region cooling wall heat flow rate per unit area and thickness of the slag crust | |
| CN108517384A (en) | A kind of monitoring method of blast furnace crucibe side wall | |
| CN109929955A (en) | A kind of detection method of blast furnace crucible corrosion situation | |
| Zagaria et al. | Monitoring erosion and skull profile in blast furnace hearth | |
| CN116362085B (en) | Hearth lining erosion morphology identification method based on cooling wall heat flow intensity | |
| CN110765623A (en) | Two-dimensional rapid calculation method for erosion boundary of longitudinal section of blast furnace hearth | |
| KR102531803B1 (en) | Method for monitoring wear of refractory linings of blast furnaces | |
| CN113111549A (en) | Method and system for modeling erosion model after pouring and repairing blast furnace hearth | |
| CN105463142A (en) | Method for measuring temperature of molten iron inside blast-furnace hearth | |
| Ding et al. | Multi‐Objective Optimization of Slab Heating Process in Walking Beam Reheating Furnace Based on Particle Swarm Optimization Algorithm | |
| CN114896546A (en) | High-precision calculation method for residual thickness of carbon brick in blast furnace hearth | |
| KR101185300B1 (en) | Method for estimating position bordered to furnace wall of softening zone | |
| TWI728743B (en) | Monitoring method for agglomeration on inner wall of blast furnace | |
| CN114756820A (en) | Method for reducing damage of cooling wall water pipe in high-temperature zone of blast furnace | |
| CN117540583B (en) | Blast furnace bottom lining erosion morphology identification method based on water temperature difference of water cooling pipe | |
| CN116817603B (en) | High-temperature smelting furnace molten pool temperature monitoring and inverting method based on heat conduction inverse problem | |
| Ghorbani et al. | Thermal assessment and identification of wear zones in a blast furnace hearth and tap-holes | |
| CN114395651A (en) | Intelligent management method and system for furnace hearth | |
| US20220298591A1 (en) | Tuyere stock, method of fabricating the same, and air blast system for melting furnace | |
| JPH0978113A (en) | Method for designing structure of furnace bottom in blast furnace | |
| CN104313224A (en) | Method for alarming leaching and thickening of hearth of vanadium titano-magnetite blast furnace without ceramic cup | |
| JP3740108B2 (en) | Blast furnace hearth structure |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |