A method, device, equipment and storage medium for determining the position of the cohesive zone in a blast furnace

Abstract

The present invention relates to the technical field of blast furnace smelting, and discloses a method, device, equipment and storage medium for determining the position of the cohesive zone in a blast furnace, including: performing burden surface scanning based on the burden charging start signal and the burden charging end signal of the current burden charging to obtain the burden surface coordinates and the burden charging speed; determining the first boundary conditions based on the burden surface coordinates and the blast furnace characteristics, and constructing an initial grid structure; updating the initial grid structure based on the burden surface coordinates and the burden charging speed to obtain a target grid structure; enabling the blast furnace burden to react, and when the multiple burden performance parameters at each grid node in the target grid structure converge, obtaining the distribution of each grid cell in the target grid structure, and generating a temperature field distribution map; determining the position of the cohesive zone in the blast furnace based on the multiple burden performance parameters corresponding to each grid node in the target grid structure and the temperature field distribution map. The present invention improves the accuracy of determining the position of the cohesive zone and reflects the dynamic changes inside the blast furnace at the same time.

Description

Technical Field

[0001] The present invention relates to the technical field of blast furnace smelting, and particularly to a method, device, equipment and storage medium for determining the position of the cohesive zone in a blast furnace. Background Art

[0002] The cohesive zone in a blast furnace refers to the area in the blast furnace where, when the burden descends to a certain depth, due to the gradually increasing temperature, the ore begins to soften and partially melt to form the initial slag. In this area, the gangue in the ore begins to undergo physical and chemical changes, such as water evaporation, carbonate decomposition, ore softening, melting, etc. The area above the cohesive zone is mainly a solid layer for drying, preheating and starting chemical reactions, and below the cohesive zone, a flowing molten layer is formed, also known as the dropping zone, where the iron oxides in the ore continue to be reduced and react with the coke, ultimately producing hot metal and slag. The position of the cohesive zone in a blast furnace is crucial for blast furnace operation because it directly affects multiple aspects such as heat energy utilization, chemical reaction efficiency and lining life in the blast furnace. Therefore, how to determine the position of the cohesive zone in a blast furnace is a problem to be solved.

[0003] Currently, devices such as infrared thermal imagers can be used to scan the blast furnace shell from the outside, and the position of the cohesive zone is inferred based on the temperature distribution. However, infrared detection mainly relies on the heat radiated from the furnace shell, and its ability to penetrate the furnace shell is limited. It can only provide the apparent temperature and is difficult to accurately reflect the actual temperature distribution deep inside the furnace, resulting in poor accuracy in determining the position of the cohesive zone and inability to reflect the dynamic changes of the blast furnace. Summary of the Invention

[0004] In view of this, the present invention provides a method, device, equipment and storage medium for determining the position of the cohesive zone in a blast furnace to solve the problem of accurately determining the position of the cohesive zone in a blast furnace.

[0005] In a first aspect, the present invention provides a method for determining the position of the cohesive zone in a blast furnace, the method comprising:

[0006] Performing a burden surface scan based on the burden charging start signal and the burden charging end signal of the current charging to obtain the burden surface coordinates and the burden charging speed of the current charging, where the burden surface coordinates include the coordinates of each point on the top surface of the blast furnace after the current charging is completed;

[0007] Determining a first boundary condition based on the burden surface coordinates and the blast furnace characteristics, and constructing an initial grid structure based on the first boundary condition, where the initial grid structure is obtained by dividing the physical model of the blast furnace into grids under the constraint of the first boundary condition;

[0008] Updating the initial grid structure based on the burden surface coordinates and the burden charging speed to obtain the target grid structure of the current charging;

[0009] Enable the burden materials in the blast furnace to undergo three-phase flow mass, heat, and momentum transfer and chemical reactions of solid-liquid-gas. When the multiple burden material performance parameters corresponding to each grid node in the target grid structure converge, obtain the distribution of each grid cell in the target grid structure, and generate a temperature field distribution map. The distribution includes the solid temperature, gas temperature, flow rate, and pressure in the area where the grid cell is located.

[0010] Based on the multiple burden material performance parameters corresponding to each grid node in the target grid structure and the temperature field distribution map, determine the position of the cohesive zone in the blast furnace after this batch of burden distribution.

[0011] The method for determining the position of the cohesive zone in the blast furnace provided by the embodiments of the present invention scans the burden surface in real time according to the burden distribution start and end signals, obtains the burden surface coordinates and the burden discharging speed of this batch of burden distribution. Thus, based on the burden surface coordinates and the burden discharging speed, update the initial grid structure constructed based on the first boundary conditions to obtain the target grid structure. By simulating the three-phase flow mass, heat, and momentum transfer and chemical reactions in the grid structure corresponding to this batch of burden distribution, considering the changes of the burden materials over time and space, it is possible to understand the complex physical and chemical changes in the furnace more deeply. Thereby, generate a temperature field distribution map. Combining the burden material performance parameters and the temperature field distribution map can more accurately determine the position of the cohesive zone, provide strong data support for production scheduling and process optimization, improve the accuracy of determining the position of the cohesive zone, and at the same time reflect the dynamic changes inside the blast furnace.

[0012] In an optional embodiment, based on the burden distribution start signal and the burden distribution end signal of this batch of burden distribution, perform burden surface scanning to obtain the burden surface coordinates and the burden discharging speed of this batch of burden distribution, including:

[0013] When the material flow valve in the blast furnace is opened, obtain the burden distribution start signal of this batch of burden distribution;

[0014] Based on the burden distribution start signal, perform burden surface scanning to obtain the starting burden surface coordinates of this batch of burden distribution;

[0015] When the sounding rod in the blast furnace is lifted, obtain the burden distribution end signal of this batch of burden distribution;

[0016] Based on the burden distribution end signal, perform burden surface scanning to obtain the burden surface coordinates of this batch of burden distribution;

[0017] Based on the burden surface coordinates of this batch of burden distribution, the starting burden surface coordinates, and the time difference between the burden distribution start signal and the burden distribution end signal, determine the burden discharging speed of this batch of burden distribution.

[0018] The method for determining the position of the cohesive zone in the blast furnace provided by the embodiments of the present invention supports the subsequent determination of the position of the cohesive zone by real-time monitoring of the burden distribution start and end signals, respectively scanning to obtain the burden surface coordinates, and determining the burden discharging speed based on the burden surface coordinates obtained from the two scans.

[0019] In an alternative embodiment, based on the burden surface coordinates and the burden feeding speed, the initial grid structure is updated to obtain the target grid structure for this burden distribution, including:

[0020] Determine the grid parameters of each grid node in the initial grid structure, where the grid parameters include radial distance, vertical distance, fluid radial velocity, and fluid vertical velocity;

[0021] For any grid node corresponding to the burden surface coordinates in the initial grid structure, lower the grid node and determine the grid parameters of the lowered grid node obtained after lowering;

[0022] Take the lowered grid node as the new grid node, and repeat the process of lowering the grid node and obtaining the grid parameters of the lowered grid node until the grid node or the lowered grid node exceeds the physical domain, to obtain the grid parameters of at least one lowered grid node corresponding to the grid node;

[0023] Based on the grid parameters of at least one lowered grid node obtained by each grid node corresponding to the burden surface coordinates in the initial grid structure during the lowering process, update the grid parameters of each grid node in the initial grid structure to obtain the target grid structure.

[0024] The method for determining the position of the cohesive zone in the blast furnace provided by the embodiments of the present invention determines the grid parameters of each lowered grid node during the lowering process by lowering the grid nodes on the top burden surface in the initial grid structure. Thus, according to the grid parameters of the lowered grid nodes obtained by each grid node on the top burden surface during the lowering process, the grid parameters of the grid nodes in the initial grid structure are updated, simulating the lowering process of the burden over time and space, being able to reflect the dynamic changes inside the blast furnace, obtaining the target grid structure that conforms to the actual situation of this burden distribution, and providing support for subsequent determination of the position of the cohesive zone.

[0025] In an alternative embodiment, determining the grid parameters of each grid node in the initial grid structure includes:

[0026] Perform coordinate transformation on the burden surface coordinates to obtain the grid parameters of each point in the burden surface coordinates;

[0027] Based on the grid parameters of each point in the burden surface coordinates, determine the grid parameters of each grid node in the initial grid structure.

[0028] The method for determining the position of the cohesive zone in a blast furnace provided by the embodiments of the present invention converts the burden surface coordinates into grid parameters, ensuring that the physical properties of each point match the actual burden state, which helps to more precisely understand the furnace conditions. Based on the grid parameters of each point in the burden surface coordinates, the grid parameters of each grid node in the initial grid structure are determined, thereby obtaining the initial grid structure with complete grid node information.

[0029] In an alternative embodiment, for any grid node corresponding to the burden surface coordinates in the initial grid structure, the grid node is lowered, and the grid parameters of the lowered grid node obtained after lowering are determined, including:

[0030] Determine the lowering time of the grid node;

[0031] Based on the radial distance, fluid radial velocity, and lowering time of the grid node, determine the radial distance of the lowered grid node;

[0032] Based on the vertical distance, fluid vertical velocity, and lowering time of the grid node, determine the vertical distance of the lowered grid node;

[0033] Based on the radial distance and vertical distance of the lowered grid node, determine the grid cell to which the lowered grid node belongs in the initial grid structure;

[0034] Based on the grid parameters of the grid node, the radial distance and vertical distance of the lowered grid node, and the grid parameters of the four vertex grid nodes of the grid cell, determine the four distances between the lowered grid node and the four vertex grid nodes of the grid cell;

[0035] Based on the grid parameters of the four vertex grid nodes of the grid cell and the four distances, determine the fluid radial velocity and fluid vertical velocity of the lowered grid node.

[0036] The method for determining the position of the cohesive zone in a blast furnace provided by the embodiments of the present invention can simulate the sinking behavior of the burden in the furnace by calculating the grid parameters of the lowered grid node obtained after the grid node is lowered, providing intuitive visual and quantitative data support for understanding the dynamic changes in the furnace.

[0037] In an alternative embodiment, the solid-liquid-gas three-phase flow, heat and mass transfer, and chemical reactions of the burden in the blast furnace are carried out. When the multiple burden performance parameters corresponding to each grid node in the target grid structure converge, the distribution of each grid cell in the target grid structure is obtained, and a temperature field distribution map is generated, including:

[0038] Determine the second boundary condition;

[0039] Enable the burden in the blast furnace to undergo the three - transfer of solid - liquid - gas and chemical reactions for the first time, and obtain multiple burden property parameters corresponding to each grid node in the target grid structure and the distribution of each grid cell after the first three - transfer of solid - liquid - gas and chemical reactions;

[0040] Repeat the above process of enabling the burden in the blast furnace to undergo the three - transfer of solid - liquid - gas and chemical reactions until the multiple burden property parameters corresponding to each grid node converge under the constraint of the second boundary condition, obtain the distribution of each grid cell in the target grid structure, and generate a temperature field distribution map.

[0041] The method for determining the position of the cohesive zone in the blast furnace provided by the embodiments of the present invention, by setting the second boundary condition, repeating the simulation of three - phase mass transfer and chemical reactions until the burden property parameters converge, obtains the distribution of each grid cell, thereby generating a temperature field distribution map, visually showing the spatial distribution of the temperature in the furnace, and providing support for subsequent determination of the position of the cohesive zone.

[0042] In an alternative embodiment, based on the multiple burden property parameters corresponding to each grid node in the target grid structure and the temperature field distribution map, determine the position of the cohesive zone in the blast furnace after this burden charging, including:

[0043] Based on the multiple burden property parameters corresponding to each grid node in the target grid structure, conduct high - temperature droplet performance detection to obtain the softening start temperature and dripping start temperature of the burden;

[0044] Based on the temperatures of all grid cells in the temperature field distribution map in the target grid structure, regard all grid cells with temperatures between the softening start temperature and the dripping start temperature as the position of the cohesive zone in the blast furnace after this burden charging.

[0045] The method for determining the position of the cohesive zone in the blast furnace provided by the embodiments of the present invention, by conducting high - temperature droplet performance detection when using the multiple burden property parameters corresponding to each grid node in the target grid structure as the detection environment, can predict the softening start temperature and dripping start temperature of the burden, and thus determine the grid cells with temperatures between the two temperatures in the target grid structure as the position of the cohesive zone, realizing the accurate determination of the position of the cohesive zone.

[0046] In a second aspect, the present invention provides a device for determining the position of the cohesive zone in a blast furnace, and the device includes:

[0047] A scanning module, configured to perform burden surface scanning based on the burden charging start signal and the burden charging end signal of this burden charging, obtain the burden surface coordinates and the burden discharging speed of this burden charging, and the burden surface coordinates include the coordinates of each point on the top surface of the blast furnace after this burden charging ends;

[0048] A construction module, configured to determine a first boundary condition based on the burden surface coordinates and blast furnace characteristics, and construct an initial grid structure based on the first boundary condition. The initial grid structure is obtained by meshing the physical model of the blast furnace under the constraint of the first boundary condition;

[0049] An update module, configured to update the initial grid structure based on the burden surface coordinates and the burden charging speed to obtain the target grid structure for this burden distribution;

[0050] A generation module, configured to enable the burden in the blast furnace to undergo three-phase flow heat, mass, and momentum transfer and chemical reactions of solid-liquid-gas phases. When the multiple burden performance parameters corresponding to each grid node in the target grid structure converge, obtain the distribution of each grid cell in the target grid structure, and generate a temperature field distribution diagram. The distribution includes the solid temperature, gas temperature, flow rate, and pressure in the area where the grid cell is located;

[0051] A determination module, configured to determine the position of the cohesive zone in the blast furnace after this burden distribution based on the multiple burden performance parameters corresponding to each grid node in the target grid structure and the temperature field distribution diagram.

[0052] In a third aspect, the present invention provides a computer device, including: a memory and a processor, which are communicatively connected to each other. The memory stores computer instructions, and the processor executes the computer instructions to execute the method for determining the position of the cohesive zone in the blast furnace according to the first aspect or any corresponding embodiment thereof.

[0053] In a fourth aspect, the present invention provides a computer-readable storage medium, on which computer instructions are stored. The computer instructions are used to cause a computer to execute the method for determining the position of the cohesive zone in the blast furnace according to the first aspect or any corresponding embodiment thereof.

[0054] In a fifth aspect, the present invention provides a computer program product, including computer instructions, which are used to cause a computer to execute the method for determining the position of the cohesive zone in the blast furnace according to the first aspect or any corresponding embodiment thereof. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0056] Figure 1 is a flowchart of the method for determining the position of the cohesive zone in the blast furnace according to an embodiment of the present invention;

[0057] Figure 2 It is a schematic diagram of the initial grid structure according to an embodiment of the present invention;

[0058] Figure 3 It is a schematic diagram of the distribution of the burden surface at the top of the blast furnace according to an embodiment of the present invention;

[0059] Figure 4 It is a schematic diagram of the grid unit according to an embodiment of the present invention;

[0060] Figure 5 It is a schematic diagram of the target grid structure according to an embodiment of the present invention;

[0061] Figure 6 It is a schematic diagram of the cohesive zone of the blast furnace according to an embodiment of the present invention;

[0062] Figure 7 It is a flowchart of a method for determining the position of the cohesive zone of the blast furnace according to another embodiment of the present invention;

[0063] Figure 8 It is a structural block diagram of a device for determining the position of the cohesive zone of the blast furnace according to an embodiment of the present invention;

[0064] Figure 9 It is a schematic diagram of the hardware structure of a computer device according to an embodiment of the present invention. Detailed implementation manners

[0065] To make the objectives, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Apparently, the described embodiments are some but not all of the embodiments of the present invention. All other embodiments obtained by those skilled in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0066] The position of the cohesive zone in a blast furnace is crucial for blast furnace operation as it directly affects multiple aspects such as the heat energy utilization, chemical reaction efficiency, and lining life inside the blast furnace. By using equipment such as an infrared thermal imager to scan the blast furnace shell from the outside and inferring the position of the cohesive zone based on the temperature distribution, this method is difficult to accurately reflect the actual temperature distribution deep inside the furnace, resulting in poor accuracy in determining the position of the cohesive zone and being unable to reflect the dynamic changes of the blast furnace. In the embodiments of the present invention, by scanning the burden surface in real time according to the charging start and end signals, the burden surface coordinates and the burden discharging speed of this charging are obtained, the initial grid structure is updated to obtain the target grid structure, and by simulating the three transports of solid-liquid-gas and chemical reactions, considering the changes of the burden over time and space, a deeper understanding of the complex physical and chemical changes inside the furnace can be achieved, thereby generating a temperature field distribution map. By integrating the burden performance parameters and the temperature field distribution map, the position of the cohesive zone can be determined more accurately, improving the accuracy of determining the position of the cohesive zone and reflecting the dynamic changes inside the blast furnace at the same time.

[0067] According to an embodiment of the present invention, there is provided an embodiment of a method for determining the position of the cohesive zone in a blast furnace. It should be noted that the steps shown in the flowchart of the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although the logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in a different order than here.

[0068] In this embodiment, a method for determining the position of the cohesive zone in a blast furnace is provided, which can be used in an electronic device. Figure 1 It is a flowchart of the method for determining the position of the cohesive zone in a blast furnace according to an embodiment of the present invention, as Figure 1 shown, and this process includes the following steps:

[0069] Step S101, perform burden surface scanning based on the charging start signal and the charging end signal of this charging to obtain the burden surface coordinates and the burden discharging speed of this charging. The burden surface coordinates include the coordinates of each point on the top surface of the blast furnace after the end of this charging. Specifically, the burden of blast furnace charging includes coke and ore. As the coke in the tuyere raceway at the bottom of the blast furnace burns and the molten iron and slag are discharged out of the furnace, the blast furnace burden continuously descends. To keep the blast furnace working continuously, charging needs to be continuously carried out at the top of the blast furnace. Each time during charging, coke or ore will become the top burden surface of the blast furnace. Optionally, an on-line laser burden surface scanning device of the blast furnace can be used for burden surface scanning, and only the top burden surface of the blast furnace needs to be scanned each time. By performing burden surface scanning on the burden currently being charged, the coordinates of each point on the coke burden surface or the ore burden surface can be obtained, and the discharging speed of coke or ore can be obtained through the coordinates.

[0070] Step S102: Based on the burden surface coordinates and the blast furnace characteristics, determine the first boundary conditions, and construct an initial grid structure based on the first boundary conditions. The initial grid structure is obtained by meshing the physical model of the blast furnace under the constraints of the first boundary conditions. Specifically, the grid structure is obtained by dividing the physical model of the blast furnace into multiple grid cells. The shape and size of each grid cell are not fixed, but each has four grid nodes as vertices. Each grid node can be represented by a radial distance and a vertical distance, and can reflect the fluid radial velocity and fluid vertical velocity of the burden at this grid node after the burden in the furnace becomes a fluid through the reactions in the blast furnace. The blast furnace characteristics include the boundary furnace profile, the blast furnace centerline, and the shape of the dead coke bed. Figure 2 is a schematic diagram of the initial grid structure according to an embodiment of the present invention, as Figure 2 shown. The top burden surface is obtained from the burden surface coordinates after this batch of charging. The boundary furnace profile is composed of the furnace shaft, furnace waist, furnace belly, hearth, and the depth of the raceway. The first boundary conditions are the boundaries formed by the boundary furnace profile, the shape of the dead coke bed, the blast furnace centerline, and the burden surface coordinates. Mesh the physical model of the blast furnace within the boundaries formed by the first boundary conditions to obtain an initial grid structure that conforms to the shape of the blast furnace and the burden situation of this batch of charging.

[0071] Step S103: Based on the burden surface coordinates and the burden charging speed, update the initial grid structure to obtain the target grid structure for this batch of charging. Specifically, by the burden situation of this batch of charging, updating the initial grid structure can obtain a target grid structure that better fits the actual situation in the furnace currently under the consideration of the dynamic changes of the burden in the furnace.

[0072] Step S104: Cause the burden in the blast furnace to undergo three-phase flow mass, momentum, and heat transfer and chemical reactions of solid-liquid-gas. When the multiple burden property parameters corresponding to each grid node in the target grid structure converge, obtain the distribution of each grid cell in the target grid structure, and generate a temperature field distribution map. The distribution includes the solid temperature, gas temperature, flow rate, and pressure in the area where the grid cell is located. Specifically, when the temperature in the furnace reaches a certain threshold, the iron ore in the burden begins to soften and form initial melting. As the temperature further increases, part of the burden melts and drips, forming the softening-melting zone in the blast furnace. Therefore, the temperature factor needs to be considered emphatically when determining the position of the softening-melting zone in the blast furnace. During the ironmaking process in the blast furnace, three-phase flow mass, momentum, and heat transfer and chemical reactions mainly occur, causing changes in the mass, momentum, and heat of the burden in the blast furnace. Couple and calculate the above mass, momentum, heat, and chemical reactions at each grid node. When the burden property parameters of each grid node converge, obtain the distribution of each grid cell in the grid structure corresponding to this burden distribution, and thus generate a temperature field distribution map based on the distribution of each grid cell to reflect the temperature at various parts inside the blast furnace. Among them, the burden property parameters include the edge solid temperature, gas pressure, silicon content, top gas temperature and composition. Optionally, other burden property parameters can also be determined according to the actual situation, and the embodiments of the present invention do not limit this. By generating the temperature field distribution map, it is possible to accurately determine the softening-melting zone of the blast furnace based on the temperature at various parts inside the blast furnace.

[0073] Step S105: Based on the multiple burden property parameters corresponding to each grid node in the target grid structure and the temperature field distribution map, determine the position of the softening-melting zone of the blast furnace after this burden distribution. Specifically, through the burden property parameters corresponding to each grid node and the temperature field distribution map, the situation inside the furnace can be fully understood, and thus the position of the softening-melting zone of the blast furnace can be determined more accurately.

[0074] The method for determining the position of the softening-melting zone of the blast furnace provided by the embodiments of the present invention scans the burden surface in real time according to the start and end signals of the burden distribution, obtains the burden surface coordinates and the burden feeding speed of this burden distribution, and thus updates the initial grid structure constructed based on the first boundary conditions through the burden surface coordinates and the burden feeding speed to obtain the target grid structure. By simulating three-phase flow mass, momentum, and heat transfer and chemical reactions in the grid structure corresponding to this burden distribution, considering the changes of the burden over time and space, the complex physical and chemical changes inside the furnace can be understood more deeply, and thus a temperature field distribution map is generated. Combining the burden property parameters and the temperature field distribution map, the position of the softening-melting zone can be determined more precisely, providing strong data support for production scheduling and process optimization, achieving the improvement of the accuracy of determining the position of the softening-melting zone, and at the same time reflecting the dynamic changes inside the blast furnace.

[0075] In this embodiment, a method for determining the position of the softening-melting zone of the blast furnace is provided, which can be used in the above-mentioned electronic device. The method specifically includes the following steps:

[0076] Step S201: Based on the cloth start signal and the cloth end signal of this cloth, perform a burden surface scan to obtain the burden surface coordinates and the burden charging speed of this cloth. The burden surface coordinates include the coordinates of each point on the top surface of the blast furnace after the end of this cloth.

[0077] Specifically, the above step S201 includes:

[0078] Step S2011: When the material flow valve of the blast furnace is opened, obtain the cloth start signal of this cloth. Specifically, when the material flow valve of the blast furnace is opened, cloth can be carried out, that is, coke or ore is added to the blast furnace. At this time, the cloth start signal is obtained.

[0079] Step S2012: Based on the cloth start signal, perform a burden surface scan to obtain the start burden surface coordinates of this cloth. Specifically, when the cloth start signal is obtained, in order to make the scanned burden surface coordinates more accurate and comprehensive, delay for a preset time before performing the burden surface scan. At this time, the cloth has stabilized, and the burden surface scan starts through the on-line laser burden surface scanning device of the blast furnace to obtain the start burden surface coordinates. Optionally, the preset delay time and the burden surface scan time can be set by oneself, and the embodiments of the present invention do not limit this.

[0080] Step S2013: When the sounding rod of the blast furnace is lifted, obtain the cloth end signal of this cloth. Specifically, when the sounding rod of the blast furnace is lifted, the cloth ends, that is, the addition of coke or ore to the blast furnace stops. At this time, the cloth end signal is obtained.

[0081] Step S2014: Based on the cloth end signal, perform a burden surface scan to obtain the burden surface coordinates of the burden of this cloth. Specifically, when the cloth end signal is obtained, considering that there is a large amount of dust after the cloth ends, which has an adverse effect on the laser scanning effect of the on-line laser burden surface scanning device of the blast furnace, so delay for a preset time after obtaining the cloth end signal before performing the burden surface scan to obtain the burden surface coordinates. Optionally, different preset delay times can be set for different burdens. For example, the preset delay time after the end of coke cloth is set to 10 seconds, and the preset delay time after the end of ore cloth is set to 20 seconds. Optionally, the above-mentioned preset delay times corresponding to coke and ore are only examples, and the embodiments of the present invention do not limit this.

[0082] Step S2015: Determine the charging speed of the burden materials for this charging based on the burden surface coordinates of this batch of materials, the starting burden surface coordinates, and the time difference between the charging start signal and the charging end signal. Specifically, divide the height difference between the burden surface coordinates and the starting burden surface coordinates by the time difference between the two signals to obtain the charging speed of the burden materials. By real-time monitoring the charging start and end signals, scan to obtain the burden surface coordinates respectively, and determine the charging speed of the burden materials based on the burden surface coordinates obtained from the two scans, providing support for subsequent determination of the position of the cohesive zone.

[0083] In some alternative embodiments, the central part of the blast furnace is a flame zone, and the laser of the on-line laser burden surface scanning device of the blast furnace cannot scan out the burden surface. Therefore, the coordinates of the central part can be calculated using the Lagrange interpolation algorithm. Concatenate the burden surface coordinates directly obtained by scanning and the burden surface coordinates of the central part obtained by the difference method to form the complete burden surface coordinates of the burden materials.

[0084] In some alternative embodiments, Figure 3 is a schematic diagram of the distribution of the burden surface at the top of the blast furnace according to an embodiment of the present invention, as Figure 3 shown. The burden surface coordinates are obtained each time the blast furnace is charged. For example, when charging ore this time, the burden surface coordinates of the ore during this charging are obtained through burden surface scanning, that is, Figure 3 the upper boundary of the ore shown.

[0085] Step S202: Determine the first boundary conditions based on the burden surface coordinates and the characteristics of the blast furnace, and construct an initial grid structure based on the first boundary conditions. The initial grid structure is obtained by meshing the physical model of the blast furnace under the constraint of the first boundary conditions. For details, please refer to Figure 1 Step S102 of the embodiment shown, which will not be elaborated here.

[0086] Step S203: Update the initial grid structure based on the burden surface coordinates and the charging speed of the burden materials to obtain the target grid structure for this charging.

[0087] Specifically, the above-mentioned Step S203 includes:

[0088] Step S2031: Determine the grid parameters of each grid node in the initial grid structure. The grid parameters include radial distance, vertical distance, fluid radial velocity, and fluid vertical velocity.

[0089] In some alternative embodiments, the above-mentioned Step S2031 includes:

[0090] Step a1: Convert the coordinates of the burden surface to obtain the grid parameters of each point in the burden surface coordinates. Specifically, convert the coordinates of the burden surface of this batch of burden distribution through the Laplace equation. In the orthogonal coordinate system, according to the finite difference method, calculate the grid parameters of each grid node corresponding to the burden surface coordinates in the initial grid structure through the following formula (1). Optionally, the above processes are all prior arts and will not be elaborated here.

[0091] (1)

[0092] Where, represents the stream function; represents the radial distance; represents the vertical distance.

[0093] Step a2: Based on the grid parameters of each point in the burden surface coordinates, determine the grid parameters of each grid node in the initial grid structure. Specifically, since the burden surface coordinates are obtained by scanning the top surface of the blast furnace, they correspond to the top grid nodes of the initial grid structure. After determining the network parameters of each grid node located at the top layer, prior arts can be used. For example, by using the convective trajectory equation, starting from the known top grid nodes, according to the physical properties and movement laws of the burden, calculate the grid parameters of other grid nodes in the initial grid structure except the top grid nodes. The specific process will not be elaborated here.

[0094] Step S2032: For any grid node corresponding to the burden surface coordinates in the initial grid structure, lower the grid node and determine the grid parameters of the lowered grid node obtained after lowering.

[0095] In some alternative embodiments, the above Step S2032 includes:

[0096] Step b1: Determine the lowering time of the grid node. Specifically, assume that any grid node is lowered and the lowering time is preset. Optionally, in the subsequent repeated lowering process, the embodiments of the present invention take the same lowering time for each lowering as an example for illustration.

[0097] Step b2: Based on the radial distance, fluid radial velocity, and lowering time of the grid node, determine the radial distance of the lowered grid node. Specifically, the radial distance of the lowered grid node of this grid node can be determined through the following formula (2).

[0098] (2)

[0099] Where, represents the index of the lowered grid node located in the th row and the th column; Indicates that it is located at Rank The index of the grid node of the column; Indicates that it is located at Rank The radial distance of the descending grid node of the column; Indicates that it is located at Rank The radial distance of the grid nodes of the column; Indicates that it is located at Rank The radial velocity of the fluid at the grid nodes of the column; Indicates the fall time.

[0100] Step b3: Determine the vertical distance of the descending grid node based on the vertical distance of the grid node, the vertical velocity of the fluid, and the descending time. Specifically, the vertical distance of the descending grid node of the grid node can be determined by the following formula (3).

[0101] (3)

[0102] in, Indicates that it is located at Rank The index of the descending grid node of the column; Indicates that it is located at Rank The index of the grid node of the column; Indicates that it is located at Rank The vertical distance of the descending grid node of the column; Indicates that it is located at Rank The vertical distance of the grid nodes of the column; Indicates that it is located at Rank The vertical velocity of the fluid at the grid nodes of the column; Indicates the fall time.

[0103] Step b4, based on the radial distance and vertical distance of the descending grid node, determine the grid unit to which the descending grid node belongs in the initial grid structure. Specifically, Figure 4 is a schematic diagram of a grid unit according to an embodiment of the present invention, such as Figure 4 As shown, assuming that point P is a grid node and point M is a descending grid node of point P, the grid unit where point M is located is determined from the initial grid structure, that is, the grid unit composed of Q1-Q2-Q3-Q4.

[0104] In step b5, based on the grid parameters of the grid node, the radial distance and vertical distance of the descending grid node, and the grid parameters of the four vertex grid nodes of the grid cell, the four distances between the descending grid node and the four vertex grid nodes of the grid cell are determined. Specifically, the grid cell where point M is located is Q1-Q2-Q3-Q4. The distances between point M and the four vertices of the grid cell are determined using the following equations (4) to (7), that is, the distances between point M and points Q1, Q2, Q3, and Q4 are determined respectively.

[0105] (4)

[0106] (5)

[0107] (6)

[0108] (7)

[0109] in, Indicates the distance between point M and point Q1; Represents the distance between point M and point Q2; Indicates the distance between point M and point Q3; Indicates the distance between point M and point Q4; Indicates that it is located at Rank The radial distance of the Q1 point of the column; represents the horizontal coordinate of point M located in the Q1-Q2-Q3-Q4 grid cell; Indicates that it is located at Rank The vertical distance of point Q1 of the column; represents the ordinate of point M located in the Q1-Q2-Q3-Q4 grid cell; Indicates that it is located at Rank The radial distance of the Q2 point of the column; Indicates that it is located at Rank The vertical distance of point Q2 of the column; Indicates that it is located at Rank The radial distance of the Q3 point of the column; Indicates that it is located at Rank The vertical distance of point Q3 of the column; Indicates that it is located at Rank The radial distance of the Q4 point of the column; Indicates that it is located at Rank The vertical distance of the Q4 point in the column;

[0110] Step b6: Determine the fluid radial velocity and fluid vertical velocity of the descending grid node based on the grid parameters of the four vertex grid nodes of the grid cell and the four distances. Specifically, the interpolation function can be constructed by the following equations (8) and (9) to determine the fluid radial velocity and fluid vertical velocity of the descending grid node. By calculating the grid parameters of the descending grid node obtained after the grid node descends, the sinking behavior of the burden in the furnace can be simulated, providing intuitive visual and quantitative data support for understanding the dynamic changes in the furnace.

[0111] (8)

[0112] (9)

[0113] Wherein, represents the fluid radial velocity of point M located in the th row and the th column; represents the fluid vertical velocity of point M located in the th row and the th column; represents the fluid radial velocity of point Q1 located in the th row and the th column; represents the fluid vertical velocity of point Q1 located in the th row and the th column; represents the distance between point M and point Q1; represents the fluid radial velocity of point Q2 located in the th row and the th column; represents the fluid vertical velocity of point Q2 located in the th row and the th column; represents the distance between point M and point Q2; represents the fluid radial velocity of point Q3 located in the th row and the th column; represents the fluid vertical velocity of point Q3 located in the th row and the th column; represents the distance between point M and point Q3; represents the fluid radial velocity of point Q2 located in the th row and the th column; represents the fluid vertical velocity of point Q2 located in the th row and the th column; Indicates the distance between point M and point Q4.

[0114] Step S2033: Take the descending grid nodes as new grid nodes, and repeat the process of descending the grid nodes and obtaining the grid parameters of the descending grid nodes until the grid nodes or the descending grid nodes exceed the physical domain, so as to obtain the grid parameters of at least one descending grid node corresponding to the grid nodes. Specifically, for each grid node located at the top layer in the initial grid structure, make it descend in the initial grid structure. After the descending time, determine a descending grid node and its grid parameters, then take this descending grid node as a new grid node and continue to descend, and determine the next descending grid node and its grid parameters. Repeat the above process until the grid nodes or the descending grid nodes exceed the physical domain, and obtain the grid parameters of all descending grid nodes corresponding to each grid node at the top layer of the initial grid structure. Among them, the physical domain can be referred to as Figure 2 the grid area shown.

[0115] Step S2034: Update the grid parameters of each grid node in the initial grid structure based on the grid parameters of at least one descending grid node obtained during the descent of each grid node corresponding to the burden surface coordinates in the initial grid structure, so as to obtain the target grid structure. Specifically, for each grid node corresponding to the burden surface coordinates in the initial grid structure, that is, each grid node located at the top layer of the initial grid structure, based on the grid parameters of all descending grid nodes obtained by simulating the descent of the burden for each grid node, correspondingly update the grid parameters of the grid node in the initial grid structure that is in the same position as the descending grid node. The descent process of the burden over time and space is simulated, which can reflect the dynamic changes inside the blast furnace, and obtain the target grid structure that conforms to the actual situation of this burden distribution, providing support for subsequent determination of the softening-melting zone position. Figure 5 is a schematic diagram of the target grid structure according to an embodiment of the present invention, Figure 5 The target grid structure shown relative to Figure 2 the initial grid structure shown is more in line with the actual situation of this burden distribution.

[0116] Step S204: Make the burden in the blast furnace undergo three-phase flow heat, mass and momentum transfer and chemical reactions. When the multiple burden performance parameters corresponding to each grid node in the grid structure corresponding to this burden distribution converge, obtain the distribution of each grid cell in the grid structure corresponding to this burden distribution, and generate a temperature field distribution diagram. The distribution includes the solid temperature, gas temperature, flow rate and pressure in the area where the grid cell is located.

[0117] Specifically, the above step S204 includes:

[0118] Step S2041, determine the second boundary conditions. Specifically, the furnace body thermocouple, furnace body static pressure, hot metal composition, top gas temperature, and top gas composition are used as the second boundary conditions. Among them, the furnace body thermocouple can be obtained through the Fourier heat conduction equation, that is, the following formula (10), which represents the temperature of the solid burden in the furnace wall part; the hot metal composition refers to the actual content of silicon in the hot metal.

[0119] (10)

[0120] Among them, represents heat; represents time; represents the thermal conductivity; represents the furnace body area; represents the furnace body temperature; represents the furnace body length.

[0121] Step S2042, make the burden in the blast furnace undergo the first solid-liquid-gas three-phase flow mass, momentum, and heat transfer and chemical reactions for the first time, and obtain the multiple burden performance parameters corresponding to each grid node and the distribution of each grid cell in the target grid structure after the first solid-liquid-gas three-phase flow mass, momentum, and heat transfer and chemical reactions. Specifically, the burden in the blast furnace undergoes solid-liquid-gas three-phase flow mass, momentum, and heat transfer and chemical reactions during the blast furnace ironmaking process. Among them, the solid-liquid-gas three-phase flow mass, momentum, and heat transfer involve the mass, momentum, and heat transfer of the burden, and the chemical reactions involve the reduction reactions of iron, carbon, and silicon. The burden performance parameters correspond to the above second boundary conditions and include the edge solid temperature, gas pressure, silicon composition, and top gas temperature and composition. Make the first solid-liquid-gas three-phase flow mass, momentum, and heat transfer and chemical reactions occur inside the blast furnace, and obtain the burden performance parameters of each grid node and the distribution of each grid cell after the first reaction. More specifically, the reduction reaction of iron refers to the chemical reaction equation shown in the following formula (11), the reduction reaction of carbon refers to the chemical reaction equation shown in the following formula (12), and the reduction reaction of silicon refers to the chemical reaction equation shown in the following formula (13).

[0122] (11)

[0123] (12)

[0124] (13)

[0125] Among them, is ferric oxide; is carbon monoxide; is gaseous carbon monoxide; is iron oxide; is carbon dioxide; is hydrogen; is water; is iron; is carbon; is silicon; is silicon dioxide; is gaseous silicon monoxide; each represents the rate of the corresponding chemical reaction.

[0126] Step S2043: Repeat the above process of enabling the burden in the blast furnace to undergo three-phase mass, heat, and momentum transfer and chemical reactions in the solid-liquid-gas phases until the multiple burden property parameters corresponding to each grid node converge under the constraints of the second boundary conditions, obtain the distribution of each grid cell in the target grid structure, and generate a temperature field distribution map. Specifically, after the first occurrence of three-phase mass, heat, and momentum transfer and chemical reactions inside the blast furnace, determine whether the burden property parameters of each grid node after the first reaction converge under the constraints of the second boundary conditions, that is, determine whether the errors between the edge solid temperature and the set value of the furnace body thermocouple, the errors between the gas pressure and the set value of the furnace shaft static pressure, the errors between the silicon composition and the set value of the hot metal composition, and the errors between the top gas temperature and composition in the burden property parameters and the set values of the top gas temperature and composition in the second boundary conditions are all within the preset error range. If the error between each item of the burden property parameters of each grid node and the corresponding item in the second boundary conditions is within the preset error range, it is considered to have converged. If any grid node does not converge after the first reaction, perform the second three-phase mass, heat, and momentum transfer and chemical reactions, and continue to judge whether the burden property parameters after the second reaction converge until, after any occurrence of three-phase mass, heat, and momentum transfer and chemical reactions, each grid node converges. Then, based on the distribution of each grid cell obtained after this reaction, generate a temperature field distribution map. Among them, this temperature field distribution map can visually display the temperature, flow rate, and pressure at each location inside the blast furnace. By setting the second boundary conditions, repeating the three-phase mass transfer and chemical reaction simulation until the burden property parameters converge, obtaining the distribution of each grid cell, and thus generating a temperature field distribution map, it visually shows the spatial distribution of the temperature inside the furnace, providing support for determining the position of the cohesive zone in the subsequent process.

[0127] Step S205: Based on the multiple burden property parameters corresponding to each grid node in the grid structure corresponding to this burden distribution and the temperature field distribution map, determine the position of the cohesive zone in the blast furnace after this burden distribution.

[0128] Specifically, the above Step S205 includes:

[0129] Step S2051: Based on multiple burden performance parameters corresponding to each grid node in the target grid structure, perform high-temperature dripping performance detection to obtain the softening start temperature and dripping start temperature of the burden. Specifically, high-temperature dripping performance detection is a material testing method that focuses on the morphology, melting characteristics, flow behavior of the droplets, and their mechanical properties at high temperatures. Therefore, it can be used to determine the softening start temperature and dripping start temperature of the burden. Use the multiple burden performance parameters corresponding to each grid node when convergence is reached in the above step S2043 as the environmental parameters in the high-temperature dripping performance detection, that is, perform the detection in the environment corresponding to the burden performance parameters to obtain the softening start temperature and dripping start temperature of the burden.

[0130] Step S2052: Based on the temperatures of all grid cells in the temperature field distribution map of the target grid structure, take all grid cells with temperatures between the softening start temperature and the dripping start temperature as the position of the cohesive zone of the blast furnace after this burden distribution. Specifically, since the temperature field distribution map shows the temperature at each location inside the blast furnace, the grid cells in the target grid structure with temperatures between the two temperatures are determined as the position of the cohesive zone, realizing the accurate determination of the position of the cohesive zone. Figure 6 is a schematic diagram of the cohesive zone of the blast furnace according to an embodiment of the present invention, as Figure 6 shown. The left side is a center-suppressed cohesive zone of the blast furnace, and the right side is an ideal cohesive zone of the blast furnace. The shape of the cohesive zone of the blast furnace is not fixed, Figure 6 and the shown is only an example.

[0131] In some alternative embodiments, Figure 7 is a flowchart of another method for determining the position of the cohesive zone of the blast furnace according to an embodiment of the present invention, as Figure 7 shown. First, perform a burden surface scan on the burden after this burden distribution to obtain the burden surface coordinates and the burden charging speed. Then, based on the first boundary condition, the burden surface coordinates, and the blast furnace characteristics, construct an initial grid structure. Then, based on the burden surface coordinates and the burden charging speed, update the initial grid structure to obtain the target grid structure. Then, let the burden in the blast furnace undergo three-phase flow mass, heat, and momentum transfer and chemical reactions to obtain multiple burden performance parameters corresponding to each grid node and the distribution of each grid cell. Then, based on the second boundary condition, determine whether the multiple burden performance parameters corresponding to each grid node converge. If they converge, generate a temperature field distribution map based on the distribution of each grid cell. If they do not converge, return to the step of letting the burden in the blast furnace undergo three-phase flow mass, heat, and momentum transfer and chemical reactions. Finally, based on the temperature field distribution map, determine the position of the cohesive zone of the blast furnace.

[0132] The method for determining the position of the cohesive zone in a blast furnace provided by the embodiments of the present invention scans the burden surface in real time according to the charging start and end signals, obtains the burden surface coordinates and the burden discharging speed of the current charging, and then updates the initial grid structure constructed based on the first boundary conditions through the burden surface coordinates and the burden discharging speed to obtain the target grid structure. By simulating the three-phase mass, heat, and momentum transfer and chemical reactions of solid-liquid-gas in the grid structure corresponding to the current charging, considering the changes of the burden over time and space, it is possible to understand the complex physical and chemical changes in the furnace more deeply, thereby generating a temperature field distribution map. Combining the burden performance parameters and the temperature field distribution map can more accurately determine the position of the cohesive zone, providing strong data support for production scheduling and process optimization, achieving the improvement of the accuracy of determining the position of the cohesive zone, and reflecting the dynamic changes inside the blast furnace at the same time.

[0133] In this embodiment, a device for determining the position of the cohesive zone in a blast furnace is also provided. This device is used to implement the above embodiments and preferred implementation manners, and those that have been described will not be repeated. As used hereinafter, the term "module" can be a combination of software and / or hardware that can achieve a predetermined function. Although the devices described in the following embodiments are preferably implemented in software, implementation in hardware, or a combination of software and hardware is also possible and contemplated.

[0134] This embodiment provides a device for determining the position of the cohesive zone in a blast furnace, as Figure 8 shown, including:

[0135] A scanning module 801, configured to perform burden surface scanning based on the charging start signal and the charging end signal of the current charging, and obtain the burden surface coordinates and the burden discharging speed of the current charging. The burden surface coordinates include the coordinates of each point on the top surface of the blast furnace after the current charging ends.

[0136] A construction module 802, configured to determine the first boundary conditions based on the burden surface coordinates and the blast furnace characteristics, and construct an initial grid structure based on the first boundary conditions. The initial grid structure is obtained by dividing the physical model of the blast furnace into grids under the constraints of the first boundary conditions.

[0137] An updating module 803, configured to update the initial grid structure based on the burden surface coordinates and the burden discharging speed to obtain the target grid structure of the current charging.

[0138] A generating module 804, configured to cause the burden in the blast furnace to perform three-phase mass, heat, and momentum transfer and chemical reactions. When the multiple burden performance parameters corresponding to each grid node in the target grid structure converge, obtain the distribution of each grid cell in the target grid structure, and generate a temperature field distribution map. The distribution includes the solid temperature, gas temperature, flow rate, and pressure in the area where the grid cell is located.

[0139] A determination module 805, configured to determine the position of the cohesive zone in the blast furnace after this batch of burden distribution based on multiple burden performance parameters and the temperature field distribution map corresponding to each grid node in the target grid structure.

[0140] In some alternative embodiments, the scanning module 801 includes:

[0141] A first acquisition unit, configured to acquire the burden distribution start signal for this batch of burden distribution when the material flow valve of the blast furnace is opened.

[0142] A first scanning unit, configured to perform burden surface scanning based on the burden distribution start signal to obtain the start burden surface coordinates for this batch of burden distribution.

[0143] A second acquisition unit, configured to acquire the burden distribution end signal for this batch of burden distribution when the sounding rod of the blast furnace is lifted.

[0144] A second scanning unit, configured to perform burden surface scanning based on the burden distribution end signal to obtain the burden surface coordinates of the burden for this batch of burden distribution.

[0145] A first determination unit, configured to determine the burden discharging speed for this batch of burden distribution based on the burden surface coordinates of the burden for this batch of burden distribution, the start burden surface coordinates, and the time difference between the burden distribution start signal and the burden distribution end signal.

[0146] In some alternative embodiments, the updating module 803 includes:

[0147] A second determination unit, configured to determine the grid parameters of each grid node in the initial grid structure, where the grid parameters include the radial distance, the vertical distance, the fluid radial velocity, and the fluid vertical velocity.

[0148] A descending unit, configured to, for any grid node corresponding to the burden surface coordinates in the initial grid structure, lower the grid node to determine the grid parameters of the lowered grid node obtained after lowering.

[0149] A third determination unit, configured to use the lowered grid node as a new grid node, repeat the process of lowering the grid node and obtaining the grid parameters of the lowered grid node until the grid node or the lowered grid node exceeds the physical domain, to obtain the grid parameters of at least one lowered grid node corresponding to the grid node.

[0150] An updating unit, configured to update the grid parameters of each grid node in the initial grid structure based on the grid parameters of at least one lowered grid node obtained by each grid node corresponding to the burden surface coordinates in the initial grid structure during the lowering process, to obtain the target grid structure.

[0151] In some alternative embodiments, the second determination unit includes:

[0152] A conversion subunit for converting the coordinates of the burden surface to obtain the grid parameters of each point in the burden surface coordinates.

[0153] A first determination subunit for determining the grid parameters of each grid node in the initial grid structure based on the grid parameters of each point in the burden surface coordinates.

[0154] In some alternative embodiments, the descent unit includes:

[0155] A second determination subunit for determining the descent time of the grid node.

[0156] A third determination subunit for determining the radial distance of the descending grid node based on the radial distance of the grid node, the radial fluid velocity, and the descent time.

[0157] A fourth determination subunit for determining the vertical distance of the descending grid node based on the vertical distance of the grid node, the vertical fluid velocity, and the descent time.

[0158] A fifth determination subunit for determining the grid cell in the initial grid structure to which the descending grid node belongs based on the radial distance and vertical distance of the descending grid node.

[0159] A sixth determination subunit for determining the four distances between the descending grid node and the four vertex grid nodes of the grid cell based on the grid parameters of the grid node, the radial distance and vertical distance of the descending grid node, and the grid parameters of the four vertex grid nodes of the grid cell.

[0160] A seventh determination subunit for determining the radial fluid velocity and vertical fluid velocity of the descending grid node based on the grid parameters of the four vertex grid nodes of the grid cell and the four distances.

[0161] In some alternative embodiments, the generation module 804 includes:

[0162] A fourth determination unit for determining the second boundary condition.

[0163] A reaction unit for enabling the burden in the blast furnace to undergo the three transfers and chemical reactions of solid-liquid-gas three-phase flow for the first time, to obtain multiple burden property parameters corresponding to each grid node in the target grid structure after the first solid-liquid-gas three-phase flow three transfers and chemical reactions and the distribution of each grid cell.

[0164] A generation unit for repeating the above process of enabling the burden in the blast furnace to undergo the three transfers and chemical reactions of solid-liquid-gas three-phase flow until the multiple burden property parameters corresponding to each grid node converge under the constraint of the second boundary condition, to obtain the distribution of each grid cell in the target grid structure and generate a temperature field distribution map.

[0165] In some alternative embodiments, the determining module 805 includes:

[0166] A detection unit configured to perform high-temperature droplet performance detection based on multiple burden performance parameters corresponding to each grid node in the target grid structure, and obtain the softening start temperature and dripping start temperature of the burden.

[0167] A fifth determining unit configured to, based on the temperatures of all grid cells in the temperature field distribution map of the target grid structure, use all grid cells with temperatures between the softening start temperature and the dripping start temperature as the position of the BF cohesive zone after this burden distribution.

[0168] The further function descriptions of the above various modules and units are the same as those in the corresponding above embodiments, and will not be elaborated here.

[0169] The device for determining the position of the BF cohesive zone in this embodiment is presented in the form of functional units. Here, the unit refers to an ASIC (Application Specific Integrated Circuit) circuit, a processor and a memory that execute one or more software or fixed programs, and / or other devices that can provide the above functions.

[0170] An embodiment of the present invention further provides a computer device having the Figure 8 device for determining the position of the BF cohesive zone as shown above.

[0171] Please refer to Figure 9 , Figure 9 which is a schematic structural diagram of a computer device provided by an alternative embodiment of the present invention. As Figure 9 shown, the computer device includes: one or more processors 10, a memory 20, and interfaces for connecting various components, including a high-speed interface and a low-speed interface. Each component communicates with each other using different buses and can be installed on a common main board or installed in other ways as needed. The processor can process instructions executed within the computer device, including instructions stored in the memory or on the memory to display graphical information of the GUI on an external input / output device (such as a display device coupled to the interface). In some alternative embodiments, if necessary, multiple processors and / or multiple buses can be used together with multiple memories and multiple memories. Similarly, multiple computer devices can be connected, and each device provides some necessary operations (such as a server array, a set of blade servers, or a multi-processor system). Figure 9 In

[0172] The processor 10 may be a central processing unit, a network processor, or a combination thereof. Among them, the processor 10 may further include a hardware chip. The above-mentioned hardware chip may be an application specific integrated circuit, a programmable logic device, or a combination thereof. The above-mentioned programmable logic device may be a complex programmable logic device, a field programmable gate array, a generic array logic, or any combination thereof.

[0173] Among them, the memory 20 stores instructions executable by at least one processor 10, so that at least one processor 10 executes the method shown in the above embodiments.

[0174] The memory 20 may include a program storage area and a data storage area. Among them, the program storage area may store an operating system and application programs required for at least one function; the data storage area may store data created according to the use of the computer device, etc. In addition, the memory 20 may include a high-speed random access memory, and may further include a non-transitory memory, such as at least one disk storage device, a flash memory device, or other non-transitory solid-state storage devices. In some alternative embodiments, the memory 20 may optionally include a memory remotely provided with respect to the processor 10, and these remote memories may be connected to the computer device through a network. Examples of the above-mentioned network include but are not limited to the Internet, an enterprise intranet, a local area network, a mobile communication network, and combinations thereof.

[0175] The memory 20 may include a volatile memory, such as a random access memory; the memory may also include a non-volatile memory, such as a flash memory, a hard disk, or a solid-state drive; the memory 20 may further include a combination of the above types of memories.

[0176] The computer device further includes a communication interface 30 for the computer device to communicate with other devices or communication networks.

[0177] The embodiments of the present invention further provide a computer-readable storage medium. The method according to the embodiments of the present invention can be implemented in hardware, firmware, or be implemented as computer code that can be recorded on a storage medium, or be implemented as computer code originally stored in a remote storage medium or a non-temporary machine-readable storage medium and downloaded through a network and to be stored in a local storage medium, so that the method described herein can be stored in such software processing on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. Among them, the storage medium may be a magnetic disk, an optical disk, a read-only memory, a random access memory, a flash memory, a hard disk, or a solid-state drive, etc.; further, the storage medium may further include a combination of the above types of memories. It can be understood that a computer, a processor, a microprocessor controller, or programmable hardware includes a storage component that can store or receive software or computer code, and when the software or computer code is accessed and executed by the computer, the processor, or the hardware, the method shown in the above embodiments is implemented.

[0178] A part of the present invention can be applied as a computer program product, such as computer program instructions, which, when executed by a computer, can invoke or provide the methods and / or technical solutions according to the present invention through the operations of the computer. Those skilled in the art should understand that the forms in which computer program instructions exist in a computer-readable medium include, but are not limited to, source files, executable files, installation package files, etc. Correspondingly, the ways in which computer program instructions are executed by a computer include, but are not limited to: the computer directly executes the instructions, or the computer compiles the instructions and then executes the corresponding compiled program, or the computer reads and executes the instructions, or the computer reads and installs the instructions and then executes the corresponding installed program. Herein, the computer-readable medium can be any available computer-readable storage medium or communication medium accessible by the computer.

[0179] Although the embodiments of the present invention have been described in conjunction with the accompanying drawings, those skilled in the art can make various modifications and variations without departing from the spirit and scope of the present invention, and such modifications and variations all fall within the scope defined by the present invention.

Claims

1. A method for determining the position of the cohesive zone in a blast furnace, characterized in that, The method includes: Using an on-line laser burden surface scanning device for a blast furnace to perform burden surface scanning based on the burden charging start signal and the burden charging end signal of this burden charging, obtaining the burden surface coordinates of the burden materials and the burden discharging speed of the burden materials in this burden charging. The burden surface coordinates include the coordinates of each point located on the top surface of the blast furnace after the completion of this burden charging. The burden surface coordinates of the central part of the blast furnace burden surface are calculated by the Lagrange interpolation algorithm from the burden surface coordinates scanned by the on-line laser burden surface scanning device for the blast furnace. The burden surface coordinates of the burden materials include the burden surface coordinates scanned by the on-line laser burden surface scanning device for the blast furnace and the burden surface coordinates of the central part of the blast furnace burden surface; Based on the burden surface coordinates of the burden materials and the blast furnace characteristics, determining the first boundary conditions, and constructing an initial grid structure based on the first boundary conditions. The initial grid structure is obtained by performing grid division on the physical model of the blast furnace under the constraints of the first boundary conditions. The blast furnace characteristics include the boundary furnace type, the blast furnace center line, and the shape of the dead coke bed; Based on the burden surface coordinates of the burden materials and the burden discharging speed, updating the initial grid structure to obtain the target grid structure of this burden charging; Enabling the burden materials in the blast furnace to undergo three-phase flow heat, mass, and momentum transfer and chemical reactions of solid-liquid-gas. When the multiple burden material performance parameters corresponding to each grid node in the target grid structure converge, obtaining the distribution of each grid unit in the target grid structure, and generating a temperature field distribution map. The distribution includes the solid temperature, gas temperature, flow rate, and pressure in the area where the grid unit is located; Based on the multiple burden material performance parameters corresponding to each grid node in the target grid structure and the temperature field distribution map, determining the position of the cohesive zone in the blast furnace after this burden charging; Among them, the performing burden surface scanning based on the burden charging start signal and the burden charging end signal of this burden charging, and obtaining the burden surface coordinates of the burden materials and the burden discharging speed of the burden materials includes: When the material flow valve of the blast furnace is opened, obtaining the burden charging start signal of this burden charging; Performing burden surface scanning after delaying a preset time based on the burden charging start signal to obtain the starting burden surface coordinates of this burden charging; When the sounding rod of the blast furnace is lifted, obtaining the burden charging end signal of this burden charging; Performing burden surface scanning after delaying a preset time based on the burden charging end signal to obtain the burden surface coordinates of the burden materials in this burden charging; Based on the burden surface coordinates of the burden materials in this burden charging, the starting burden surface coordinates, and the time difference between the burden charging start signal and the burden charging end signal, determining the burden discharging speed of the burden materials in this burden charging. The burden discharging speed is the quotient obtained by dividing the height difference between the burden surface coordinates of the burden materials and the starting burden surface coordinates by the time difference; The updating the initial grid structure based on the burden surface coordinates of the burden materials and the burden discharging speed to obtain the target grid structure of this burden charging includes: Determining the grid parameters of each grid node in the initial grid structure. The grid parameters include the radial distance, the vertical distance, the fluid radial velocity, and the fluid vertical velocity; For any grid node corresponding to the burden surface coordinates in the initial grid structure, lower the grid node and determine the grid parameters of the lowered grid node obtained after lowering; Take the lowered grid node as a new grid node, repeat the process of lowering the grid node and obtaining the grid parameters of the lowered grid node until the grid node or the lowered grid node exceeds the physical domain, and obtain the grid parameters of at least one lowered grid node corresponding to the grid node; Based on the grid parameters of at least one lowered grid node obtained during the lowering process of each grid node corresponding to the burden surface coordinates in the initial grid structure, update the grid parameters of each grid node in the initial grid structure to obtain the target grid structure; The step of, for any grid node corresponding to the burden surface coordinates in the initial grid structure, lowering the grid node and determining the grid parameters of the lowered grid node obtained after lowering, includes: Determine the lowering time of the grid node; Based on the radial distance of the grid node, the radial velocity of the fluid, and the lowering time, determine the radial distance of the lowered grid node; Based on the vertical distance of the grid node, the vertical velocity of the fluid, and the lowering time, determine the vertical distance of the lowered grid node; Based on the radial distance and the vertical distance of the lowered grid node, determine the grid cell in which the lowered grid node is located in the initial grid structure; Based on the grid parameters of the grid node, the radial distance and the vertical distance of the lowered grid node, and the grid parameters of the four vertex grid nodes of the grid cell, determine the four distances between the lowered grid node and the four vertex grid nodes of the grid cell; Based on the grid parameters of the four vertex grid nodes of the grid cell and the four distances, determine the radial velocity and the vertical velocity of the fluid of the lowered grid node.

2. The method according to claim 1, wherein The step of determining the grid parameters of each grid node in the initial grid structure includes: Perform coordinate transformation on the burden surface coordinates to obtain the grid parameters of each point in the burden surface coordinates; Based on the grid parameters of each point in the burden surface coordinates, determine the grid parameters of each grid node in the initial grid structure.

3. The method according to claim 1, characterized in that, The step of making the burden in the blast furnace undergo three-phase flow mass, heat, and momentum transfer and chemical reactions, and obtaining the distribution of each grid cell in the target grid structure and generating a temperature field distribution map when the multiple burden performance parameters corresponding to each grid node in the target grid structure converge, includes: Determine the second boundary condition; Make the burden in the blast furnace first undergo three-phase flow mass, heat, and momentum transfer and chemical reactions to obtain the multiple burden performance parameters corresponding to each grid node in the target grid structure after the first three-phase flow mass, heat, and momentum transfer and chemical reactions and the distribution of each grid cell; Repeat the above process of enabling the burden materials in the blast furnace to undergo three-phase flow mass, heat, and momentum transfer and chemical reactions in the solid-liquid-gas phases until the multiple burden material property parameters corresponding to each grid node converge under the constraints of the second boundary condition, obtain the distribution of each grid cell in the target grid structure, and generate the temperature field distribution map.

4. The method according to claim 1, wherein Determining the position of the cohesive zone in the blast furnace after this batch of burden distribution based on the multiple burden material property parameters corresponding to each grid node in the target grid structure and the temperature field distribution map includes: Performing high-temperature droplet property detection based on the multiple burden material property parameters corresponding to each grid node in the target grid structure to obtain the softening start temperature and dripping start temperature of the burden materials. Based on the temperatures of all grid cells in the target grid structure in the temperature field distribution map, take all grid cells with temperatures between the softening start temperature and the dripping start temperature as the position of the cohesive zone in the blast furnace after this batch of burden distribution.

5. A device for determining the position of the cohesive zone in a blast furnace, characterized in that, The device includes: A scanning module, configured to perform burden surface scanning using an on-line laser burden surface scanning device for the blast furnace based on the burden distribution start signal and the burden distribution end signal of this batch of burden distribution to obtain the burden surface coordinates and the burden discharging speed of this batch of burden distribution. The burden surface coordinates include the coordinates of each point on the top surface of the blast furnace after the end of this batch of burden distribution. The burden surface coordinates at the central part of the blast furnace burden surface are calculated by the Lagrange interpolation algorithm using the burden surface coordinates scanned by the on-line laser burden surface scanning device for the blast furnace. The burden surface coordinates include the burden surface coordinates scanned by the on-line laser burden surface scanning device for the blast furnace and the burden surface coordinates at the central part of the blast furnace burden surface. A construction module, configured to determine the first boundary condition based on the burden surface coordinates and the blast furnace characteristics, and construct an initial grid structure based on the first boundary condition. The initial grid structure is obtained by performing grid division on the physical model of the blast furnace under the constraints of the first boundary condition. The blast furnace characteristics include the boundary furnace shape, the blast furnace center line, and the shape of the dead coke bed. An update module, configured to update the initial grid structure based on the burden surface coordinates and the burden discharging speed to obtain the target grid structure of this batch of burden distribution. A generation module, configured to enable the burden materials in the blast furnace to undergo three-phase flow mass, heat, and momentum transfer and chemical reactions. When the multiple burden material property parameters corresponding to each grid node in the target grid structure converge, obtain the distribution of each grid cell in the target grid structure and generate a temperature field distribution map. The distribution includes the solid temperature, gas temperature, flow rate, and pressure in the area where the grid cell is located. A determination module, configured to determine the position of the cohesive zone in the blast furnace after this batch of burden distribution based on the multiple burden material property parameters corresponding to each grid node in the target grid structure and the temperature field distribution map. Among them, the scanning module is specifically configured to: Obtain the burden distribution start signal of this batch of burden distribution when the material flow valve of the blast furnace is opened. Perform burden surface scanning with a time delay of a preset time based on the burden distribution start signal to obtain the start burden surface coordinates of this batch of burden distribution. Obtain the burden distribution end signal of this batch of burden distribution when the sounding rod of the blast furnace is lifted. Perform a burden surface scan based on the cloth end signal delayed by a preset time to obtain the burden surface coordinates of the current cloth. Based on the burden surface coordinates of the current cloth, the start burden surface coordinates, and the time difference between the cloth start signal and the cloth end signal, determine the burden discharging speed of the current cloth. The burden discharging speed is the quotient obtained by dividing the height difference between the burden surface coordinates and the start burden surface coordinates by the time difference. The update module is specifically configured to: Determine the grid parameters of each grid node in the initial grid structure. The grid parameters include radial distance, vertical distance, fluid radial velocity, and fluid vertical velocity. For any grid node corresponding to the burden surface coordinates in the initial grid structure, lower the grid node and determine the grid parameters of the lowered grid node obtained after lowering. Take the lowered grid node as a new grid node, and repeat the process of lowering the grid node and obtaining the grid parameters of the lowered grid node until the grid node or the lowered grid node exceeds the physical domain, to obtain the grid parameters of at least one lowered grid node corresponding to the grid node. The step of, for any grid node corresponding to the burden surface coordinates in the initial grid structure, lowering the grid node and determining the grid parameters of the lowered grid node obtained after lowering includes: Determine the lowering time of the grid node. Based on the radial distance of the grid node, the fluid radial velocity, and the lowering time, determine the radial distance of the lowered grid node. Based on the vertical distance of the grid node, the fluid vertical velocity, and the lowering time, determine the vertical distance of the lowered grid node. Based on the radial distance and vertical distance of the lowered grid node, determine the grid cell in the initial grid structure to which the lowered grid node belongs. Based on the grid parameters of the grid node, the radial distance and vertical distance of the lowered grid node, and the grid parameters of the four vertex grid nodes of the grid cell, determine the four distances between the lowered grid node and the four vertex grid nodes of the grid cell. Based on the grid parameters of the four vertex grid nodes of the grid cell and the four distances, determine the fluid radial velocity and fluid vertical velocity of the lowered grid node.

6. A computer device, characterized in that, It includes: A memory and a processor, which are communicatively connected to each other. The memory stores computer instructions, and the processor executes the computer instructions to execute the method for determining the position of the cohesive zone in a blast furnace according to any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that, Computer instructions are stored on the computer-readable storage medium, and the computer instructions are used to cause a computer to execute the method for determining the position of the cohesive zone in a blast furnace according to any one of claims 1 to 4.

Images