Method and apparatus for determining position of cohesive zone in blast furnace, device, and storage medium

Abstract

The present application relates to the technical field of blast furnace smelting, and discloses a method and apparatus for determining the position of a cohesive zone in a blast furnace, a device, and a storage medium. The method comprises: on the basis of a burden-distributing start signal and a burden-distributing end signal for the current burden distribution, performing burden surface scanning to obtain burden surface coordinates and a burden charging rate; determining a first boundary condition on the basis of the burden surface coordinates and blast furnace features, and constructing an initial grid structure; updating the initial grid structure on the basis of the burden surface coordinates and the burden charging rate to obtain a target grid structure; enabling burden in the blast furnace to undergo reaction, when a plurality of burden performance parameters of 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; and determining the position of a cohesive zone in the blast furnace on the basis of the plurality of burden performance parameters corresponding to each grid node in the target grid structure and the temperature field distribution map. The present application improves the accuracy of determining the position of the cohesive zone, and reflects dynamic changes inside the blast furnace.

Description

A method, apparatus, equipment, and storage medium for determining the location of the softening zone in a blast furnace.

[0001] Cross-reference to related applications

[0002] This application claims priority to Chinese Patent Application No. 202411877341.7, filed on December 19, 2024, entitled "A method, apparatus, device and storage medium for determining the position of the blast furnace softening zone", the entire contents of which are incorporated herein by reference. Technical Field

[0003] This application relates to the field of blast furnace smelting technology, specifically to a method, apparatus, equipment, and storage medium for determining the location of the softening zone in a blast furnace. Background Technology

[0004] The softening zone of a blast furnace refers to the area within the furnace where, as the burden descends to a certain depth, the ore begins to soften and partially melt to form primary slag due to the gradual increase in temperature. Within this zone, the gangue in the ore undergoes physicochemical changes, such as moisture evaporation, carbonate decomposition, ore softening, and melting. Above the softening zone is primarily a solid layer for drying, preheating, and initiation of chemical reactions. Below the softening zone forms a flowing molten layer, also known as the drip zone, where iron oxides in the ore continue to be reduced and react with coke, ultimately producing molten iron and slag. The location of the blast furnace softening zone is crucial for blast furnace operation, directly affecting factors such as heat energy utilization, chemical reaction efficiency, and furnace lining life. Therefore, determining the location of the blast furnace softening zone is a problem that needs to be solved.

[0005] Currently, infrared thermal imagers and other equipment can be used to scan the blast furnace shell from the outside and infer the location of the softening zone 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 cannot accurately reflect the actual temperature distribution deep inside the furnace, resulting in poor accuracy in determining the location of the softening zone and failing to reflect the dynamic changes of the blast furnace. Summary of the Invention

[0006] In view of this, this application provides a method, apparatus, equipment and storage medium for determining the location of the blast furnace softening zone, so as to solve the problem of accurately determining the location of the blast furnace softening zone.

[0007] In a first aspect, this application provides a method for determining the location of the softening zone in a blast furnace, the method comprising:

[0008] Based on the start and end signals of this feeding operation, the material surface is scanned to obtain the furnace charge surface coordinates and the charge feeding speed. The furnace charge surface coordinates include the coordinates of each point on the top surface of the blast furnace after the feeding operation ends.

[0009] Based on the coordinates of the charge surface and the characteristics of the blast furnace, the first boundary condition is determined, and an initial mesh structure is constructed based on the first boundary condition. The initial mesh structure is obtained by meshing the physical model of the blast furnace under the constraint of the first boundary condition.

[0010] Based on the furnace charge surface coordinates and furnace charge feeding speed, the initial mesh structure is updated to obtain the target mesh structure for this material distribution.

[0011] The furnace charge in the blast furnace undergoes a three-phase flow of solid, liquid, and gas, as well as chemical reactions. Under the condition that multiple furnace charge 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. The distribution includes the solid temperature, gas temperature, flow rate, and pressure in the region where the grid cell is located.

[0012] Based on the multiple furnace charge performance parameters and temperature field distribution maps corresponding to each grid node in the target grid structure, the location of the blast furnace softening zone after this charging is determined.

[0013] The method for determining the location of the blast furnace softening zone provided in this application involves scanning the charge surface in real time based on the start and end signals of the charge feeding to obtain the coordinates of the charge surface and the charge feeding speed for this charge feeding. The initial mesh structure constructed based on the first boundary condition is then updated using these coordinates and the charge feeding speed to obtain the target mesh structure. By simulating the solid-liquid-gas three-phase flow and chemical reactions within the mesh structure corresponding to this charge feeding, and considering the changes in the charge over time and space, a deeper understanding of the complex physicochemical changes within the furnace can be achieved. This generates a temperature field distribution map. By combining the charge performance parameters and the temperature field distribution map, the location of the softening zone can be determined more accurately, providing strong data support for production scheduling and process optimization, improving the accuracy of determining the location of the softening zone, and reflecting the dynamic changes inside the blast furnace.

[0014] In one optional implementation, the material surface is scanned based on the material feeding start signal and the material feeding end signal of this feeding operation to obtain the furnace charge surface coordinates and the furnace charge feeding speed, including:

[0015] When the blast furnace material flow valve is opened, the material feeding start signal for this feeding is obtained;

[0016] Based on the fabric start signal, the fabric surface is scanned to obtain the starting coordinates of the fabric surface for this fabric operation;

[0017] When the probe of the blast furnace is lifted, the signal indicating the end of the current material placement is obtained;

[0018] Based on the end signal of material feeding, the material surface is scanned to obtain the coordinates of the furnace material surface during this feeding.

[0019] Based on the furnace charge surface coordinates, the starting charge surface coordinates, and the time difference between the charge start signal and the charge end signal, the charge feeding speed for this charge feeding is determined.

[0020] The method for determining the location of the blast furnace softening zone provided in this application embodiment obtains the material surface coordinates by real-time monitoring of the start and end signals of the charging process, and determines the charging speed of the furnace charge by using the material surface coordinates obtained from the two scans, thus providing support for the subsequent determination of the softening zone location.

[0021] In one optional implementation, the initial mesh structure is updated based on the furnace charge surface coordinates and the furnace charge feeding speed to obtain the target mesh structure for this charge distribution, including:

[0022] Determine the mesh parameters for each mesh node in the initial mesh structure. The mesh parameters include radial distance, vertical distance, fluid radial velocity, and fluid vertical velocity.

[0023] For any grid node in the initial grid structure corresponding to the coordinates of the furnace charge surface, the grid node is lowered, and the grid parameters of the lowered grid node are determined.

[0024] The process of descent of the grid node and obtaining the grid parameters of the descent grid node is repeated, with the descent of the grid node as the new grid node, until the grid node or the descent grid node exceeds the physical domain, and the grid parameters of at least one descent grid node corresponding to the grid node are obtained.

[0025] Based on the grid parameters of at least one descending grid node obtained during the descent process for each grid node corresponding to the furnace charge surface coordinates in the initial grid structure, the grid parameters of each grid node in the initial grid structure are updated to obtain the target grid structure.

[0026] The method for determining the location of the blast furnace softening zone provided in this application determines the grid parameters of each descending grid node during the descent process by lowering the grid nodes of the top layer of the initial grid structure. Based on the grid parameters of the descending grid nodes obtained by each grid node in the top layer of the initial grid structure during the descent process, the grid parameters of the grid nodes in the initial grid structure are updated. This simulates the descent process of the furnace charge over time and space, reflects the dynamic changes inside the blast furnace, and obtains a target grid structure that conforms to the actual situation of the current charge distribution, providing support for the subsequent determination of the softening zone location.

[0027] In one alternative implementation, determining the mesh parameters for each mesh node in the initial mesh structure includes:

[0028] The coordinates of the furnace charge surface are transformed to obtain the grid parameters of each point in the furnace charge surface coordinates;

[0029] Based on the grid parameters of each point in the furnace charge surface coordinates, the grid parameters of each grid node in the initial grid structure are determined.

[0030] The method for determining the location of the blast furnace softening zone provided in this application converts the material surface coordinates into grid parameters, ensuring that the physical properties of each point match the actual furnace charge state. This helps to grasp the furnace condition in greater detail. Furthermore, based on the grid parameters of each point in the furnace charge 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.

[0031] In one optional implementation, for any grid node corresponding to the furnace charge surface coordinates in the initial grid structure, the grid node is lowered, and the grid parameters of the lowered grid node are determined, including:

[0032] Determine the descent time of the grid nodes;

[0033] The radial distance of the descending grid nodes is determined based on the radial distance of the grid nodes, the radial velocity of the fluid, and the descent time.

[0034] The vertical distance of the descending grid nodes is determined based on the vertical distance between grid nodes, the vertical velocity of the fluid, and the descent time.

[0035] Based on the radial and vertical distances of the descending grid nodes, determine the grid cell to which the descending grid node belongs in the initial grid structure;

[0036] Based on the mesh parameters of the mesh nodes, the radial and vertical distances of the descending mesh nodes, and the mesh parameters of the four vertex mesh nodes of the mesh cell, determine the four distances between the descending mesh nodes and the four vertex mesh nodes of the mesh cell;

[0037] Based on the mesh parameters of the four vertex mesh nodes and the four distances of the mesh cell, the radial velocity and vertical velocity of the fluid at the descending mesh node are determined.

[0038] The method for determining the location of the blast furnace softening zone provided in this application can simulate the sinking behavior of the furnace charge by calculating the grid parameters of the descending grid nodes obtained after the grid nodes descend, thus providing intuitive visual and quantitative data support for understanding the dynamic changes inside the furnace.

[0039] In one optional implementation, the furnace charge undergoes a three-phase flow and chemical reaction within the blast furnace. When multiple charge 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, generating a temperature field distribution map, including:

[0040] Determine the second boundary conditions;

[0041] The furnace charge in the blast furnace undergoes a solid-liquid-gas three-phase flow and chemical reaction for the first time, resulting in multiple furnace charge performance parameters corresponding to each grid node in the target grid structure after the first solid-liquid-gas three-phase flow and chemical reaction, as well as the distribution of each grid unit.

[0042] Repeat the above process of solid-liquid-gas three-phase flow and chemical reaction in the blast furnace charge until the multiple charge performance parameters corresponding to each grid node converge under the constraint of the second boundary condition, thereby obtaining the distribution of each grid cell in the target grid structure and generating a temperature field distribution map.

[0043] The method for determining the location of the blast furnace softening zone provided in this application involves setting a second boundary condition and repeatedly performing three-phase mass transfer and chemical reaction simulations until the furnace charge performance parameters converge, thereby obtaining the distribution of each grid cell and generating a temperature field distribution map. This map visually displays the spatial distribution of temperature inside the furnace, providing support for subsequent determination of the softening zone location.

[0044] In one optional implementation, the location of the blast furnace softening zone after this charging is determined based on multiple charge performance parameters and temperature field distribution maps corresponding to each grid node in the target grid structure, including:

[0045] Based on multiple furnace charge performance parameters corresponding to each grid node in the target grid structure, high-temperature droplet performance is detected to obtain the softening start temperature and droplet start temperature of the furnace charge.

[0046] Based on the temperature of all grid cells in the target grid structure in the temperature field distribution map, all grid cells whose temperature is between the softening start temperature and the dripping start temperature are taken as the location of the blast furnace softening zone after this charging.

[0047] The method for determining the location of the blast furnace softening zone provided in this application embodiment detects the high-temperature droplet performance by using multiple furnace charge performance parameters corresponding to each grid node in the target grid structure as the detection environment. This allows for the prediction of the softening start temperature and the droplet initiation temperature of the furnace charge, thereby identifying the grid cells in the target grid structure whose temperature is between the two temperatures as the location of the softening zone, thus achieving accurate determination of the softening zone location.

[0048] Secondly, this application provides a device for determining the location of the softening zone in a blast furnace, the device comprising:

[0049] The scanning module is used to scan the material surface based on the material feeding start signal and the material feeding end signal of this feeding to obtain the furnace charge surface coordinates and the furnace charge feeding speed. The furnace charge surface coordinates include the coordinates of each point on the top surface of the blast furnace after this feeding is completed.

[0050] The construction module is used to determine the first boundary conditions based on the furnace charge surface coordinates and blast furnace characteristics, and to construct the initial mesh structure based on the first boundary conditions. The initial mesh structure is obtained by meshing the physical model of the blast furnace under the constraints of the first boundary conditions.

[0051] The update module is used to update the initial mesh structure based on the furnace charge surface coordinates and the furnace charge feeding speed, so as to obtain the target mesh structure for this charge distribution.

[0052] The generation module is used to enable the furnace charge to undergo solid-liquid-gas three-phase flow and chemical reaction. Under the condition that the multiple furnace charge 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. The distribution includes the solid temperature, gas temperature, flow rate and pressure in the region where the grid cell is located.

[0053] The determination module is used to determine the location of the blast furnace softening zone after this charging operation, based on multiple furnace charge performance parameters and temperature field distribution maps corresponding to each grid node in the target grid structure.

[0054] Thirdly, this application 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 perform the method for determining the position of the blast furnace softening zone described in the first aspect or any corresponding embodiment.

[0055] Fourthly, this application provides a computer-readable storage medium storing computer instructions for causing a computer to execute the method for determining the location of the blast furnace softening zone described in the first aspect or any corresponding embodiment.

[0056] Fifthly, this application provides a computer program product, including computer instructions for causing a computer to execute the method for determining the location of the blast furnace softening zone as described in the first aspect or any corresponding embodiment. Attached Figure Description

[0057] To more clearly illustrate the technical solutions in the specific embodiments of this application or the prior art, the drawings used in the description of the specific embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of this application. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0058] Figure 1 is a flowchart of a method for determining the location of the blast furnace softening zone according to an embodiment of this application;

[0059] Figure 2 is a schematic diagram of the initial mesh structure according to an embodiment of this application;

[0060] Figure 3 is a schematic diagram of the distribution of the top burden surface of the blast furnace according to an embodiment of this application;

[0061] Figure 4 is a schematic diagram of a mesh cell according to an embodiment of this application;

[0062] Figure 5 is a schematic diagram of the target mesh structure according to an embodiment of this application;

[0063] Figure 6 is a schematic diagram of the blast furnace softening zone according to an embodiment of this application;

[0064] Figure 7 is a flowchart of another method for determining the location of the blast furnace softening zone according to an embodiment of this application;

[0065] Figure 8 is a structural block diagram of a device for determining the location of the blast furnace softening zone according to an embodiment of this application;

[0066] Figure 9 is a schematic diagram of the hardware structure of a computer device according to an embodiment of this application. Detailed Implementation

[0067] To make the objectives, technical solutions, and advantages of the embodiments of this application clearer, the technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0068] The location of the blast furnace's softening zone is crucial for blast furnace operation, directly impacting factors such as thermal energy utilization, chemical reaction efficiency, and lining life. While infrared thermal imagers and similar devices scan the blast furnace shell externally and infer the softening zone's location based on temperature distribution, this method struggles to accurately reflect the actual temperature distribution deep within the furnace, resulting in poor accuracy in determining the softening zone's location and failing to capture the blast furnace's dynamic changes. This application's embodiment, however, uses real-time scanning of the charge surface based on the start and end signals of the charging start and end to obtain the charge surface coordinates and charging speed for this charge operation. This updates the initial mesh structure to obtain the target mesh structure. By simulating the solid-liquid-gas three-phase flow and chemical reactions, considering the changes in the charge over time and space, a deeper understanding of the complex physicochemical changes within the furnace is achieved. This generates a temperature field distribution map. By combining the charge performance parameters with the temperature field distribution map, the location of the softening zone can be determined more accurately, improving accuracy and reflecting the dynamic changes within the blast furnace.

[0069] According to an embodiment of this application, a method for determining the location of the softening zone in a blast furnace is provided. It should be noted that the steps shown in the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions. Furthermore, although a logical order is shown in the flowchart, in some cases, the steps shown or described may be executed in a different order than that shown here.

[0070] This embodiment provides a method for determining the location of the blast furnace softening zone, which can be used in electronic devices. Figure 1 is a flowchart of the method for determining the location of the blast furnace softening zone according to an embodiment of this application. As shown in Figure 1, the process includes the following steps:

[0071] Step S101: Based on the start and end signals of this charging operation, a surface scan is performed to obtain the coordinates of the charge surface and the charging speed. The charge surface coordinates include the coordinates of each point on the top surface of the blast furnace after the current charging operation. Specifically, the blast furnace charge includes coke and ore. As the coke in the tuyeres at the bottom of the blast furnace burns and molten iron and slag are discharged, the blast furnace charge continuously descends. To ensure continuous operation of the blast furnace, charging must be continuously performed at the top. Each time charging occurs, coke or ore becomes part of the top surface of the blast furnace charge. Optionally, an online laser charge surface scanning device can be used to scan the charge surface, requiring only the top surface of the blast furnace to be scanned each time. By scanning the charge surface of the currently charging operation, the coordinates of each point on the coke or ore surface can be obtained, and the charging speed of coke or ore can be obtained from these coordinates.

[0072] Step S102: Based on the burden surface coordinates and blast furnace characteristics, determine the first boundary conditions and construct an initial mesh structure based on these conditions. The initial mesh structure is obtained by meshing the physical model of the blast furnace under the constraints of the first boundary conditions. Specifically, the mesh structure is obtained by dividing the physical model of the blast furnace into multiple mesh units. The shape and size of each mesh unit are not fixed, but each unit has four mesh nodes as vertices. Each mesh node can be represented by radial and vertical distances, reflecting the radial and vertical velocities of the burden at that mesh node after it becomes fluid through the reaction within the blast furnace. Blast furnace characteristics include the boundary furnace type, the blast furnace centerline, and the shape of the dead coke pile. Figure 2 is a schematic diagram of the initial mesh structure according to an embodiment of this application. As shown in Figure 2, the top layer of material surface is obtained from the coordinates of the furnace charge surface after this charge distribution. The boundary furnace type is composed of the furnace body, furnace waist, furnace belly, furnace hearth and the depth of the swirling zone. The first boundary condition is the boundary formed by the boundary furnace type, the shape of the dead coke pile, the blast furnace centerline and the coordinates of the furnace charge surface. The physical model of the blast furnace located within the boundary range formed by the first boundary condition is meshed to obtain an initial mesh structure that conforms to the shape of the blast furnace and the charge conditions of this charge distribution.

[0073] Step S103: Based on the furnace charge surface coordinates and the furnace charge feeding speed, update the initial mesh structure to obtain the target mesh structure for this charge distribution. Specifically, by updating the initial mesh structure based on the furnace charge conditions for this charge distribution, a target mesh structure that better reflects the current furnace conditions can be obtained, taking into account the dynamic changes in the furnace charge.

[0074] Step S104 involves subjecting the blast furnace charge to a three-phase flow (solid-liquid-gas) and chemical reaction. With the convergence of multiple charge performance parameters corresponding to each grid node in the target grid structure, the distribution of each grid cell in the target grid structure is obtained, generating a temperature field distribution map. This distribution includes the solid temperature, gas temperature, flow rate, and pressure of the region where the grid cell is located. Specifically, when the furnace temperature reaches a certain threshold, the iron ore in the charge begins to soften, forming an initial melt. As the temperature further increases, some of the charge melts and drips, forming the blast furnace softening zone. Therefore, temperature factors need to be carefully considered when determining the location of the blast furnace softening zone. In the blast furnace ironmaking process, the main processes involve solid-liquid-gas three-phase flow and chemical reactions, causing changes in the mass, momentum, and heat of the furnace charge. These mass, momentum, heat, and chemical reactions are coupled and calculated at each grid node. When the charge performance parameters of each grid node converge, the distribution of each grid cell in the current charge distribution grid structure is obtained. Based on the distribution of each grid cell, a temperature field distribution map is generated to reflect the temperature at various points inside the blast furnace. The charge performance parameters include edge solid temperature, gas pressure, silicon composition, and top gas temperature and composition. Optionally, other charge performance parameters can be determined according to actual conditions; this application does not limit this. By generating the temperature field distribution map, the softening zone of the blast furnace can be accurately determined based on the temperature at various points inside the blast furnace.

[0075] Step S105: Based on the multiple charge performance parameters and temperature field distribution maps corresponding to each grid node in the target grid structure, determine the location of the blast furnace softening zone after this charge distribution. Specifically, by using the charge performance parameters and temperature field distribution maps corresponding to each grid node, the furnace internal conditions can be fully understood, thereby enabling a more accurate determination of the blast furnace softening zone location.

[0076] The method for determining the location of the blast furnace softening zone provided in this application involves scanning the material surface in real time based on the start and end signals of the charging start and end to obtain the coordinates of the charge surface and the charging speed of the current charging operation. This allows for updating of the initial mesh structure constructed based on the first boundary condition using the charge surface coordinates and charging speed, resulting in a target mesh structure. By simulating the solid-liquid-gas three-phase flow and chemical reactions within the mesh structure corresponding to the current charging operation, and considering the changes in the charge over time and space, a deeper understanding of the complex physicochemical changes within the furnace can be achieved. This generates a temperature field distribution map. By combining the charge performance parameters and the temperature field distribution map, the location of the softening zone can be determined more accurately, providing strong data support for production scheduling and process optimization. This improves the accuracy of determining the location of the softening zone and reflects the dynamic changes inside the blast furnace.

[0077] This embodiment provides a method for determining the location of the softening zone in a blast furnace, which can be used in the aforementioned electronic equipment. The method specifically includes the following steps:

[0078] Step S201: Based on the start signal and end signal of this feeding, perform a material surface scan to obtain the furnace charge surface coordinates and the furnace charge feeding speed. The furnace charge surface coordinates include the coordinates of each point on the top surface of the blast furnace after this feeding is completed.

[0079] Specifically, step S201 includes:

[0080] Step S2011: When the blast furnace material flow valve is opened, acquire the material feeding start signal for this feeding operation. Specifically, when the blast furnace material flow valve is opened, feeding can be carried out, that is, coke or ore is added to the blast furnace, and at this time, the material feeding start signal is acquired.

[0081] Step S2012: Based on the material placement start signal, a material surface scan is performed to obtain the starting material surface coordinates for this material placement. Specifically, when the material placement start signal is obtained, to ensure that the scanned material surface coordinates are more accurate and comprehensive, a preset time is delayed before performing the material surface scan. At this time, the material placement has stabilized, and the material surface scan is started using the blast furnace online laser material surface scanning equipment to obtain the starting material surface coordinates. Optionally, the preset delay time and the material surface scan time can be set by the user, and this embodiment does not limit this.

[0082] Step S2013: When the probe of the blast furnace is lifted, obtain the end signal of the current material feeding operation. Specifically, when the probe of the blast furnace is lifted, the material feeding operation ends, that is, the addition of coke or ore to the blast furnace stops, and at this time, the end signal of the material feeding operation is obtained.

[0083] Step S2014: Based on the end signal of material distribution, perform a surface scan to obtain the coordinates of the furnace charge surface during this distribution. Specifically, when the end signal of material distribution is obtained, considering the large amount of dust after the distribution ends, which may adversely affect the laser scanning effect of the blast furnace online laser charge surface scanning equipment, a preset time delay is applied after obtaining the end signal of material distribution before performing a surface scan to obtain the coordinates of the furnace charge surface. Optionally, different preset delay times can be set for different furnace charges, for example, a preset delay time of 10 seconds after coke distribution ends and a preset delay time of 20 seconds after ore distribution ends. Optionally, the preset delay times for coke and ore are merely examples, and this application embodiment does not limit this.

[0084] Step S2015: Based on the current charge surface coordinates, the initial charge surface coordinates, and the time difference between the start and end of the charge feeding signals, the charge feeding speed is determined. Specifically, the charge feeding speed is obtained by dividing the height difference between the charge surface coordinates and the initial charge surface coordinates by the time difference between the two signals. By real-time monitoring of the start and end of the charge feeding signals, the charge surface coordinates are scanned and obtained. The charge feeding speed is determined using the charge surface coordinates obtained from the two scans, providing support for subsequently determining the location of the softening zone.

[0085] In some alternative implementations, the central part of the blast furnace is the flame zone, and the laser of the online laser material surface scanning equipment cannot scan the material surface. Therefore, the coordinates of the central part can be calculated using the Lagrange interpolation algorithm. The material surface coordinates obtained directly from the scan and the material surface coordinates of the central part obtained by the interpolation method are then combined to form the complete material surface coordinates.

[0086] In some optional embodiments, Figure 3 is a schematic diagram of the distribution of the top charge surface of the blast furnace according to an embodiment of this application. As shown in Figure 3, the coordinates of the charge surface are obtained each time the blast furnace is charged. For example, for the ore being charged this time, the coordinates of the charge surface at the time of ore charging are obtained by scanning the charge surface, which is the upper boundary of the ore shown in Figure 3.

[0087] Step S202: Based on the coordinates of the blast furnace charge surface and the characteristics of the blast furnace, the first boundary condition is determined, and an initial mesh structure is constructed based on the first boundary condition. The initial mesh structure is obtained by meshing the physical model of the blast furnace under the constraint of the first boundary condition. For details, please refer to step S102 of the embodiment shown in Figure 1, which will not be repeated here.

[0088] Step S203: Based on the furnace charge surface coordinates and the furnace charge feeding speed, update the initial mesh structure to obtain the target mesh structure for this charge distribution.

[0089] Specifically, step S203 includes:

[0090] 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.

[0091] In some optional implementations, step S2031 above includes:

[0092] Step a1: Transform the coordinates of the furnace charge surface to obtain the mesh parameters of each point in the furnace charge surface coordinates. Specifically, transform the coordinates of the furnace charge surface for this charge placement using the Laplace equation, and calculate the mesh parameters of each mesh node corresponding to the furnace charge surface coordinates in the initial mesh structure using the finite difference method in the orthogonal coordinate system according to the following formula (1). Optionally, the above process is all prior art and will not be described in detail here.

[0093] Where ψ represents the stream function; r represents the radial distance; and z represents the vertical distance.

[0094] Step a2 involves determining the mesh parameters of each mesh node in the initial mesh structure based on the mesh parameters of each point in the blast furnace charge surface coordinates. Specifically, since the blast furnace charge surface coordinates are obtained by scanning the top layer of the blast furnace charge surface, they correspond to the top layer mesh nodes in the initial mesh structure. After determining the mesh parameters of each mesh node located at the top layer, existing techniques, such as using the convection trajectory equation, can be employed to calculate the mesh parameters of other mesh nodes in the initial mesh structure, excluding the top layer mesh nodes, based on the physical properties and motion laws of the charge. The specific process will not be elaborated here.

[0095] Step S2032: For any grid node in the initial grid structure corresponding to the coordinates of the furnace charge surface, the grid node is lowered, and the grid parameters of the lowered grid node are determined.

[0096] In some optional implementations, step S2032 above includes:

[0097] Step b1: Determine the descent time for each grid node. Specifically, assuming any grid node is to be descended, the descent time is preset. Optionally, in subsequent repeated descent processes, this embodiment of the application will illustrate by using the same descent time for each descent.

[0098] Step b2: Based on the radial distance of the grid node, the radial velocity of the fluid, and the descent time, determine the radial distance of the descending grid node. Specifically, the radial distance of the descending grid node can be determined using the following equation (2): r(i a j a )=r(i b ,j b)+vr(i b ,j b )*t (2)

[0099] Among them, (i a ,j a ) indicates that it is located at the i-th position. a Line j a The index of the descending grid node of the column; (i b ,j b ) indicates that it is located at the i-th position. b Line j b The index of the grid node of the column; a(i a ,j a ) indicates that it is located at the i-th position. a Line j a The radial distance between the descending grid nodes of the column; r(i b ,j b ) indicates that it is located in row ix and j. b Radial distance between grid nodes of the column; vr(i b ,j b ) indicates that it is located at the i-th position. b Line j b The radial velocity of the fluid at each grid node in the column; t represents the descent time.

[0100] Step b3: Based on the vertical distance of the grid nodes, the vertical velocity of the fluid, and the descent time, determine the vertical distance of the descending grid nodes. Specifically, the vertical distance of the descending grid nodes can be determined using the following formula (3): z(i a ,j a )=z(i b ,j b )+vz(i b ,j b )*t (3)

[0101] Among them, (i a ,j a ) indicates that it is located at the i-th position. a Line j a The index of the descending grid node of the column; (i b i b ) indicates that it is located at the i-th position. b Line j b The index of the grid node of the column; z(i a ,j a ) indicates that it is located at the i-th position. a Line j a Vertical distance of the descending grid nodes of the column; z(i b ,j b ) indicates that it is located at the i-th position. b Line jb Vertical distance between grid nodes in a column; vz(i b ,j b ) indicates that it is located at the i-th position. b Line j b The vertical velocity of the fluid at each grid node in the column; t represents the descent time.

[0102] Step b4: Based on the radial and vertical distances of the descending grid nodes, determine the grid cell to which the descending grid node belongs in the initial grid structure. Specifically, Figure 4 is a schematic diagram of the grid cell according to an embodiment of this application. As shown in Figure 4, assuming point P is a grid node and point M is the descending grid node of point P, determine the grid cell to which point M belongs from the initial grid structure, that is, the grid cell composed of Q1-Q2-Q3-Q4.

[0103] Step b5: Based on the mesh parameters of the mesh nodes, the radial and vertical distances of the descending mesh nodes, and the mesh parameters of the four vertex mesh nodes of the mesh cell, determine the four distances between the descending mesh nodes and the four vertex mesh nodes of the mesh cell. Specifically, the mesh cell where point M is located is Q1-Q2-Q3-Q4. The distances between point M and the four vertices of this mesh cell are determined by the following formulas (4) to (7), that is, the distances between point M and points Q1, Q2, Q3, and Q4 are determined respectively.

[0104] Where s1 represents the distance between point M and point Q1; s2 represents the distance between point M and point Q2; s3 represents the distance between point M and point Q3; s4 represents the distance between point M and point Q4; r(i1,j1) represents the radial distance of point Q1 located in row i1 and column j1; nr(i M ,j M z(i1,j1) represents the x-coordinate of point M located in the Q1-Q2-Q3-Q4 grid cell; z(i1,j1) represents the vertical distance to point Q1 located in the i1th row and j1st column; nz(i M ,j M r(i2,j2) represents the ordinate of point M located in the Q1-Q2-Q3-Q4 grid cell; r(i2,j2) represents the radial distance of point Q2 located in the i2-th row and j2-th column; z(i2,j2) represents the vertical distance of point Q2 located in the i2-th row and j2-th column; r(i3,j3) represents the radial distance of point Q3 located in the i3-th row and j3-th column; z(i3,j3) represents the vertical distance of point Q3 located in the i3-th row and j3-th column; r(i4,j4) represents the radial distance of point Q4 located in the i4-th row and j4-th column; z(i4,j4) represents the vertical distance of point Q4 located in the i4-th row and j4-th column.

[0105] Step b6: Based on the mesh parameters of the four vertex mesh nodes and the four distances of the mesh element, determine the radial and vertical fluid velocities of the descending mesh nodes. Specifically, the radial and vertical fluid velocities of the descending mesh nodes can be determined by constructing an interpolation function using equations (8) and (9). By calculating the mesh parameters of the descending mesh nodes obtained after their descent, the sinking behavior of the furnace charge can be simulated, providing intuitive visual and quantitative data support for understanding the dynamic changes within the furnace.

[0106] Where, vr(i M ,j M ) indicates that it is located at the i-th position. M Line j M The radial velocity of the fluid at point M in the column; vz(i M ,j M ) indicates that it is located at the i-th position. M Line j M The vertical velocity of the fluid at point M in row i1, column j1; vr(i1,j1) represents the radial velocity of the fluid at point Q1 in row i1, column j1; vz(i1,j1) represents the vertical velocity of the fluid at point Q1 in row i1, column j1; s1 represents the distance between point M and point Q1; vr(i2,j2) represents the radial velocity of the fluid at point Q2 in row i2, column j2; vz(i2,j2) represents the vertical velocity of the fluid at point Q2 in row i2, column j2; s2 represents the distance between point M and point Q2. The distance between points; vr(i3,j3) represents the radial velocity of the fluid at point Q3 located in row i3 and column j3; vx(i3,j3) represents the vertical velocity of the fluid at point Q3 located in row i3 and column j3; s3 represents the distance between point M and point Q3; w(i4,j4) represents the radial velocity of the fluid at point Q2 located in row i4 and column j4; vz(i4,j4) represents the vertical velocity of the fluid at point Q2 located in row i4 and column j4; s4 represents the distance between point M and point Q4.

[0107] Step S2033: The descending mesh node is used as a new mesh node. The process of descending the mesh node and obtaining its mesh parameters is repeated until the mesh node or the descending mesh node exceeds the physical domain, thus obtaining the mesh parameters of at least one descending mesh node corresponding to the original mesh node. Specifically, for each mesh node located at the top layer of the initial mesh structure, it is descended within the initial mesh structure. After a descent time, a descending mesh node and its mesh parameters are determined. This descending mesh node is then used as a new mesh node to continue the descent process, and the next descending mesh node and its mesh parameters are determined. This process is repeated until the mesh node or the descending mesh node exceeds the physical domain, thus obtaining the mesh parameters of all descending mesh nodes corresponding to each mesh node at the top layer of the initial mesh structure. The physical domain can be seen in the mesh region shown in Figure 2.

[0108] Step S2034: Based on the mesh parameters of at least one descending mesh node obtained during the descent process of each mesh node corresponding to the furnace charge surface coordinates in the initial mesh structure, the mesh parameters of each mesh node in the initial mesh structure are updated to obtain the target mesh structure. Specifically, for each mesh node corresponding to the furnace charge surface coordinates in the initial mesh structure, that is, each mesh node located at the top layer of the initial mesh structure, based on the mesh parameters of all descending mesh nodes obtained by simulating the descent of the furnace charge, the mesh parameters of the mesh nodes in the initial mesh structure located at the same position as the descending mesh nodes are updated accordingly. This simulates the descent process of the furnace charge over time and space, which can reflect the dynamic changes inside the blast furnace and obtain a target mesh structure that conforms to the actual situation of this charge distribution, providing support for the subsequent determination of the location of the softening zone. Figure 5 is a schematic diagram of the target mesh structure according to an embodiment of this application. The target mesh structure shown in Figure 5 is more consistent with the actual situation of this charge distribution than the initial mesh structure shown in Figure 2.

[0109] Step S204 involves the three-phase flow and chemical reaction of the charge in the blast furnace, which is solid-liquid-gas. When the multiple charge performance parameters corresponding to each grid node in the grid structure corresponding to this charge distribution converge, the distribution of each grid cell in the grid structure corresponding to this charge distribution is obtained, and a temperature field distribution map is generated. The distribution includes the solid temperature, gas temperature, flow rate and pressure in the region where the grid cell is located.

[0110] Specifically, step S204 includes:

[0111] Step S2041: Determine the second boundary conditions. Specifically, the furnace body thermocouple, furnace body static pressure, molten iron composition, furnace top gas temperature, and furnace top gas composition are taken as the second boundary conditions. Among them, the furnace body thermocouple can be obtained by the Fourier heat conduction equation, that is, the following equation (10), which represents the temperature of the solid furnace charge in the furnace wall part; the molten iron composition refers to the actual silicon content in the molten iron.

[0112] Where q represents heat; u represents time; k represents thermal conductivity; A represents furnace area; U represents furnace temperature; and x represents furnace length.

[0113] Step S2042 involves inducing the first solid-liquid-gas three-phase flow and chemical reaction in the blast furnace charge, obtaining multiple charge performance parameters corresponding to each grid node in the target grid structure after the first solid-liquid-gas three-phase flow and chemical reaction, as well as the distribution of each grid cell. Specifically, the charge in the blast furnace undergoes a solid-liquid-gas three-phase flow and chemical reaction during blast furnace ironmaking. The solid-liquid-gas three-phase flow involves the charge's mass, momentum, and heat transfer, while the chemical reaction involves the reduction of iron, carbon, and silicon. The charge performance parameters correspond to the aforementioned second boundary conditions, including edge solid temperature, gas pressure, silicon composition, and the temperature and composition of the furnace top gas. This process induces the first solid-liquid-gas three-phase flow and chemical reaction inside the blast furnace, obtaining the charge performance parameters for each grid node and the distribution of each grid cell after the first reaction. More specifically, the reduction reaction of iron is shown in the chemical reaction equation (11) below, the reduction reaction of carbon is shown in the chemical reaction equation (12) below, and the reduction reaction of silicon is shown in the chemical reaction equation (13) below.

[0114] Wherein, Fe₂O₃ is ferric oxide; CO is carbon monoxide; CO(g) is gaseous carbon monoxide; FeO is iron oxide; CO₂ is carbon dioxide; H₂ is hydrogen; H₂O is water; Fe is iron; C is carbon; Si is silicon; SiO₂ is silicon dioxide; SiO(g) is gaseous silicon monoxide; each R * This indicates the rate of the corresponding chemical reaction.

[0115] Step S2043: Repeat the above process of solid-liquid-gas three-phase flow and chemical reaction in the blast furnace until the multiple charge performance parameters corresponding to each grid node converge under the constraint of the second boundary condition, obtaining the distribution of each grid cell in the target grid structure and generating a temperature field distribution map. Specifically, after the first solid-liquid-gas three-phase flow and chemical reaction occurs inside the blast furnace, determine whether the charge performance parameters of each grid node after the first reaction have converged under the constraint of the second boundary condition. That is, determine whether the errors between the edge solid temperature and the set value of the furnace body thermocouple, the errors between the set values ​​of gas pressure and furnace body static pressure, the errors between the set values ​​of silicon composition and molten iron composition, and the errors between the furnace top gas temperature and composition in the charge performance parameters and the set values ​​of furnace top gas temperature and composition in the second boundary condition are all within the preset error range. If the error between each item of the charge performance parameters of each grid node and the corresponding item in the second boundary condition is within the preset error range, then convergence is considered to have been achieved. If any grid node fails to converge after the first reaction, a second solid-liquid-gas three-phase flow mass transfer and chemical reaction is performed. The convergence of the furnace charge performance parameters after the second reaction is then assessed again until every grid node converges after any solid-liquid-gas three-phase flow mass transfer and chemical reaction. Based on the distribution of each grid cell obtained after this reaction, a temperature field distribution map is generated. This temperature field distribution map visually displays the temperature, flow rate, and pressure at every point inside the blast furnace. By setting a second boundary condition, the three-phase flow mass transfer and chemical reaction simulation is repeated until the furnace charge performance parameters converge, obtaining the distribution of each grid cell and generating the temperature field distribution map. This visually displays the spatial distribution of temperature within the furnace, providing support for subsequently determining the location of the softening zone.

[0116] Step S205: Based on the multiple furnace charge performance parameters and temperature field distribution map corresponding to each grid node in the grid structure corresponding to this charge placement, determine the location of the blast furnace softening zone after this charge placement.

[0117] Specifically, step S205 includes:

[0118] Step S2051 involves performing high-temperature droplet performance testing based on multiple charge performance parameters corresponding to each grid node in the target mesh structure to obtain the charge's softening and dripping initiation temperatures. Specifically, high-temperature droplet performance testing is a materials testing method that focuses on the droplet's morphology, melting characteristics, flow behavior, and mechanical properties at high temperatures; therefore, it can be used to determine the charge's softening and dripping initiation temperatures. The multiple charge performance parameters corresponding to each grid node at the convergence point in step S2043 are used as environmental parameters in the high-temperature droplet performance testing; that is, the testing is performed in the environment corresponding to the charge performance parameters to obtain the charge's softening and dripping initiation temperatures.

[0119] Step S2052: Based on the temperature of all grid cells in the target grid structure in the temperature field distribution map, all grid cells whose temperature is between the softening start temperature and the dripping start temperature are taken as the positions of the blast furnace softening zone after this charging. Specifically, since the temperature field distribution map shows the temperature at every point inside the blast furnace, the grid cells in the target grid structure whose temperature is between two temperatures are determined as the softening zone positions, achieving accurate determination of the softening zone position. Figure 6 is a schematic diagram of the blast furnace softening zone according to an embodiment of this application. As shown in Figure 6, the left side is a center-suppressed blast furnace softening zone, and the right side is an ideal blast furnace softening zone. The shape of the blast furnace softening zone is not fixed; Figure 6 is only an example.

[0120] In some optional embodiments, Figure 7 is a flowchart of another method for determining the location of the blast furnace softening zone according to an embodiment of this application. As shown in Figure 7, firstly, the charge surface is scanned after the current charging to obtain the charge surface coordinates and the charge feeding speed. Then, based on the first boundary condition, the charge surface coordinates, and blast furnace characteristics, an initial grid structure is constructed. Then, based on the charge surface coordinates and the charge feeding speed, the initial grid structure is updated to obtain the target grid structure. Then, the charge in the blast furnace undergoes solid-liquid-gas three-phase flow and chemical reaction to obtain multiple charge performance parameters corresponding to each grid node and the distribution of each grid cell. Then, based on the second boundary condition, it is determined whether the multiple charge performance parameters corresponding to each grid node converge. If converged, a temperature field distribution map is generated based on the distribution of each grid cell. If not converged, the process returns to the step of undergoing solid-liquid-gas three-phase flow and chemical reaction in the blast furnace. Finally, based on the temperature field distribution map, the location of the blast furnace softening zone is determined.

[0121] The method for determining the location of the blast furnace softening zone provided in this application involves scanning the material surface in real time based on the start and end signals of the charging start and end to obtain the coordinates of the charge surface and the charging speed of the current charging operation. This allows for updating of the initial mesh structure constructed based on the first boundary condition using the charge surface coordinates and charging speed, resulting in a target mesh structure. By simulating the solid-liquid-gas three-phase flow and chemical reactions within the mesh structure corresponding to the current charging operation, and considering the changes in the charge over time and space, a deeper understanding of the complex physicochemical changes within the furnace can be achieved. This generates a temperature field distribution map. By combining the charge performance parameters and the temperature field distribution map, the location of the softening zone can be determined more accurately, providing strong data support for production scheduling and process optimization. This improves the accuracy of determining the location of the softening zone and reflects the dynamic changes inside the blast furnace.

[0122] This embodiment also provides a device for determining the location of the blast furnace softening zone. This device is used to implement the above embodiments and optional implementations, and details already described will not be repeated. As used below, the term "module" can refer to a combination of software and / or hardware that performs a predetermined function. Although the device described in the following embodiments is preferably implemented in software, hardware implementation, or a combination of software and hardware, is also possible and contemplated.

[0123] This embodiment provides a device for determining the location of the softening zone in a blast furnace, as shown in Figure 8, including:

[0124] The scanning module 801 is used to scan the material surface based on the material feeding start signal and the material feeding end signal of this feeding, so as to obtain the furnace material surface coordinates and the furnace material feeding speed of this feeding. The furnace material surface coordinates include the coordinates of each point located on the top surface of the blast furnace after this feeding is completed.

[0125] The construction module 802 is used to determine the first boundary conditions based on the furnace charge surface coordinates and blast furnace characteristics, and to construct an initial mesh structure based on the first boundary conditions. The initial mesh structure is obtained by meshing the physical model of the blast furnace under the constraints of the first boundary conditions.

[0126] The update module 803 is used to update the initial mesh structure based on the furnace charge surface coordinates and the furnace charge feeding speed, so as to obtain the target mesh structure for this material distribution.

[0127] The generation module 804 is used to enable the furnace charge in the blast furnace to undergo solid-liquid-gas three-phase flow and chemical reaction. When the multiple furnace charge 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. The distribution includes the solid temperature, gas temperature, flow rate and pressure in the region where the grid cell is located.

[0128] The determination module 805 is used to determine the location of the blast furnace softening zone after this charging based on multiple furnace charge performance parameters and temperature field distribution maps corresponding to each grid node in the target grid structure.

[0129] In some alternative implementations, the scanning module 801 includes:

[0130] The first acquisition unit is used to acquire the material feeding start signal when the material flow valve of the blast furnace is opened.

[0131] The first scanning unit is used to scan the material surface based on the fabric start signal to obtain the starting material surface coordinates of the fabric in this operation.

[0132] The second acquisition unit is used to acquire the end signal of the current material placement when the probe of the blast furnace is lifted.

[0133] The second scanning unit is used to scan the material surface based on the material feeding end signal to obtain the coordinates of the furnace material surface during this feeding.

[0134] The first determining unit is used to determine the charge feeding speed of the current charge feeding based on the coordinates of the charge surface, the starting charge surface coordinates, and the time difference between the charge feeding start signal and the charge feeding end signal.

[0135] In some alternative implementations, the update module 803 includes:

[0136] The second determining unit is used 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.

[0137] The descent element is used to descend any grid node corresponding to the furnace charge surface coordinates in the initial grid structure, and to determine the grid parameters of the descended grid node.

[0138] The third determining unit is used to take the descending grid node as a new grid node, repeat the process of descending the grid node and obtaining the grid parameters of the descending grid node, until the grid node or the descending grid node exceeds the physical domain, and obtain the grid parameters of at least one descending grid node corresponding to the grid node.

[0139] The update unit is used to update the mesh parameters of each mesh node in the initial mesh structure based on the mesh parameters of at least one descending mesh node obtained during the descent process of each mesh node corresponding to the furnace charge surface coordinates in the initial mesh structure, so as to obtain the target mesh structure.

[0140] In some optional implementations, the second determining unit includes:

[0141] The transformation sub-unit is used to transform the coordinates of the furnace charge surface to obtain the mesh parameters of each point in the furnace charge surface coordinates.

[0142] The first determining sub-unit is used to determine the mesh parameters of each mesh node in the initial mesh structure based on the mesh parameters of each point in the furnace charge surface coordinates.

[0143] In some alternative implementations, the descent unit includes:

[0144] The second determining sub-unit is used to determine the descent time of the grid nodes.

[0145] The third determining sub-unit is used to determine the radial distance of the descending grid nodes based on the radial distance of the grid nodes, the radial velocity of the fluid, and the descent time.

[0146] The fourth sub-unit is determined based on the vertical distance of the grid nodes, the vertical velocity of the fluid, and the descent time, to determine the vertical distance of the descending grid nodes.

[0147] The fifth step is to determine the sub-cell, which is used to determine the grid cell to which the descending grid node belongs in the initial grid structure based on the radial and vertical distances of the descending grid node.

[0148] The sixth sub-element is determined based on the mesh parameters of the mesh nodes, the radial and vertical distances of the descending mesh nodes, and the mesh parameters of the four vertex mesh nodes of the mesh element, to determine the four distances between the descending mesh nodes and the four vertex mesh nodes of the mesh element.

[0149] The seventh sub-element is determined based on the mesh parameters of the four vertex mesh nodes and the four distances of the mesh element, to determine the fluid radial velocity and fluid vertical velocity of the descending mesh node.

[0150] In some alternative implementations, the generation module 804 includes:

[0151] The fourth determining unit is used to determine the second boundary conditions.

[0152] The reaction unit is used to enable the furnace charge to undergo a solid-liquid-gas three-phase flow and chemical reaction for the first time, so as to obtain multiple furnace charge performance parameters corresponding to each grid node in the target grid structure after the first solid-liquid-gas three-phase flow and chemical reaction, as well as the distribution of each grid unit.

[0153] The generation unit is used to repeat the above process of solid-liquid-gas three-phase flow and chemical reaction in the furnace charge until the multiple charge performance parameters corresponding to each grid node converge under the constraint of the second boundary condition, so as to obtain the distribution of each grid cell in the target grid structure and generate the temperature field distribution map.

[0154] In some alternative implementations, the determining module 805 includes:

[0155] The detection unit is used to detect the high-temperature droplet performance based on multiple furnace charge performance parameters corresponding to each grid node in the target grid structure, and to obtain the softening start temperature and droplet start temperature of the furnace charge.

[0156] The fifth determining element is used to determine the location of the blast furnace softening zone after this charge placement, based on the temperature of all grid elements in the target grid structure in the temperature field distribution map, and to determine all grid elements whose temperature is between the softening start temperature and the dripping start temperature.

[0157] The optional functional descriptions of the above modules and units are the same as those in the corresponding embodiments described above, and will not be repeated here.

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

[0159] This application also provides a computer device having a device for determining the position of the blast furnace softening zone as shown in FIG8.

[0160] Please refer to Figure 9, which is a schematic diagram of the structure of a computer device provided in an optional embodiment of this application. As shown in Figure 9, the computer device includes: one or more processors 10, a memory 20, and interfaces for connecting the various components, including high-speed interfaces and low-speed interfaces. The various components communicate with each other using different buses and can be installed on a common motherboard or otherwise as needed. The processor can process instructions executed within the computer device, including instructions stored in or on memory to display graphical information of a GUI on an external input / output device (such as a display device coupled to the interface). In some optional embodiments, multiple processors and / or multiple buses can be used with multiple memories and multiple memory modules, if desired. Similarly, multiple computer devices can be connected, each providing some of the necessary operations (e.g., as a server array, a group of blade servers, or a multiprocessor system). Figure 9 uses one processor 10 as an example.

[0161] Processor 10 may be a central processing unit, a network processor, or a combination thereof. Optionally, processor 10 may also include a hardware chip. The hardware chip may be an application-specific integrated circuit (ASIC), a programmable logic device (PLD), or a combination thereof. The programmable logic device may be a complex programmable logic device (CAMP), a field-programmable gate array (FPGA), a general-purpose array logic (GPRS), or any combination thereof.

[0162] The memory 20 stores instructions executable by at least one processor 10 to cause at least one processor 10 to perform the method shown in the above embodiments.

[0163] The memory 20 may include a program storage area and a data storage area. The program storage area may store the operating system and applications required for at least one function; the data storage area may store data created based on the use of the computer device. Furthermore, the memory 20 may include high-speed random access memory and may also include non-transitory memory, such as at least one disk storage device, flash memory device, or other non-transitory solid-state storage device. In some alternative embodiments, the memory 20 may optionally include memory remotely located relative to the processor 10, and these remote memories may be connected to the computer device via a network. Examples of such networks include, but are not limited to, the Internet, intranets, local area networks, mobile communication networks, and combinations thereof.

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

[0165] The computer device also includes a communication interface 30 for communicating with other devices or communication networks.

[0166] This application also provides a computer-readable storage medium. The methods described in this application can be implemented in hardware or firmware, or implemented as recordable on a storage medium, or implemented as computer code downloaded over a network and originally stored on a remote storage medium or a non-transitory machine-readable storage medium and subsequently stored on a local storage medium. Thus, the methods described herein can be processed by software stored on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. The storage medium can be a magnetic disk, optical disk, read-only memory, random access memory, flash memory, hard disk, or solid-state drive, etc.; optionally, the storage medium may also include combinations of the above types of memory. It is understood that computers, processors, microprocessor controllers, or programmable hardware include storage components capable of storing or receiving software or computer code, which, when accessed and executed by the computer, processor, or hardware, implements the methods shown in the above embodiments.

[0167] A portion of this application 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 this application through the operation of the computer. Those skilled in the art will 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 executing the instructions, or the computer compiling the instructions and then executing the corresponding compiled program, or the computer reading and executing the instructions, or the computer reading and installing the instructions and then executing the corresponding installed program. Here, the computer-readable medium can be any available computer-readable storage medium or communication medium accessible to a computer.

[0168] Although embodiments of this application 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 this application, and such modifications and variations all fall within the scope defined by this application.

Claims

1. A method of determining the position of a softening zone of a blast furnace, characterized in that The method comprises: The method comprises: Based on the burden surface coordinates and the blast furnace characteristics, a first boundary condition is determined, and an initial grid structure is constructed based on the first boundary condition, the initial grid structure being obtained by grid division of a physical model of the blast furnace under the constraint of the first boundary condition, the blast furnace characteristics including a boundary furnace type, a blast furnace center line, and a dead coke pile shape; Based on the burden surface coordinates and the burden discharging speed, the initial grid structure is updated to obtain a target grid structure of the current burden distribution; The burden in the blast furnace is subjected to solid-liquid-gas three-phase flow three-way transmission and chemical reaction, and in the case that a plurality of burden performance parameters corresponding to each grid node in the target grid structure converge, the distribution of each grid element in the target grid structure is obtained, and a temperature field distribution map is generated, the distribution including the solid temperature, gas temperature, flow rate, and pressure of the region where the grid element is located; Based on the plurality of burden performance parameters corresponding to each grid node in the target grid structure and the temperature field distribution map, the position of the soft melting zone of the blast furnace after the current burden distribution is determined; The method comprises: When the blast furnace flow valve is opened, a burden start signal of the current burden distribution is obtained; Based on the burden start signal, a delay of a preset time is performed to obtain the start burden surface coordinates of the current burden distribution; When the blast furnace sounding rod is lifted, a burden end signal of the current burden distribution is obtained; Based on the burden end signal, a delay of a preset time is performed to obtain the burden surface coordinates of the current burden distribution; Based on the burden surface coordinates of the current burden distribution, the start burden surface coordinates, and the time difference between the burden start signal and the burden end signal, the burden discharging speed of the current burden distribution is determined, the burden discharging speed being the quotient obtained by dividing the height difference between the burden surface coordinates and the start burden surface coordinates by the time difference; The method comprises: The method comprises: For any grid node corresponding to the burden level coordinate in the initial grid structure, the grid node is lowered to determine the grid parameters of the lowered grid node obtained after lowering; The lowered grid node is taken as a new grid node, and the lowering of the grid node and the obtaining of the grid parameters of the lowered grid node are repeated until the grid node or the lowered grid node exceeds the physical domain, so as to 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 in the lowering process of each grid node corresponding to the burden level coordinate in the initial grid structure, the grid parameters of each grid node in the initial grid structure are updated to obtain the target grid structure; The lowering of the grid node and the obtaining of the grid parameters of the lowered grid node are repeated until the grid node or the lowered grid node exceeds the physical domain, so as to obtain the grid parameters of at least one lowered grid node corresponding to the grid node; The lowering time of the grid node is determined; Based on the radial distance of the grid node, the fluid radial velocity and the lowering time, the radial distance of the lowered grid node is determined; Based on the vertical distance of the grid node, the fluid vertical velocity and the lowering time, the vertical distance of the lowered grid node is determined; Based on the radial distance and the vertical distance of the lowered grid node, the grid cell to which the lowered grid node belongs in the initial grid structure is determined; 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, four distances between the lowered grid node and the four vertex grid nodes of the grid cell are determined; Based on the grid parameters of the four vertex grid nodes of the grid cell and the four distances, the fluid radial velocity and the fluid vertical velocity of the lowered grid node are determined.

2. The method of claim 1, wherein, The determination of the grid parameters of each grid node in the initial grid structure includes: The burden level coordinates are subjected to coordinate conversion to obtain the grid parameters of each point in the burden level coordinates; Based on the grid parameters of each point in the burden level coordinates, the grid parameters of each grid node in the initial grid structure are determined.

3. The method of claim 1, wherein, In the case where the multiple burden performance parameters corresponding to each grid node in the target grid structure converge after the solid-liquid-gas three-phase flow three-transfers and chemical reactions of the burden in the blast furnace, the distribution of each grid cell in the target grid structure is obtained, and a temperature field distribution map is generated, including: A second boundary condition is determined; The burden in the blast furnace is subjected to solid-liquid-gas three-phase flow three-transfers and chemical reactions for the first time, and the multiple burden performance parameters corresponding to each grid node in the target grid structure and the distribution of each grid cell after the first solid-liquid-gas three-phase flow three-transfers and chemical reactions are obtained. The process of solid-liquid-gas three-phase flow and chemical reaction of the burden in the blast furnace is repeated until the multiple burden performance parameters corresponding to each grid node reach convergence under the constraint of the second boundary condition, and the distribution of each grid element in the target grid structure is obtained, and a temperature field distribution map is generated.

4. The method of claim 1, wherein, The soft melting zone position of the blast furnace after the burden distribution is determined based on the multiple burden performance parameters corresponding to each grid node in the target grid structure and the temperature field distribution map, and includes: Based on the multiple burden performance parameters corresponding to each grid node in the target grid structure, high-temperature molten droplet performance detection is performed to obtain the soft melting start temperature and the droplet start temperature of the burden; Based on the temperature of all grid elements in the target grid structure in the temperature field distribution map, all grid elements with a temperature between the soft melting start temperature and the droplet start temperature are taken as the soft melting zone position of the blast furnace after the burden distribution.

5. A device for determining the position of the softening zone of a blast furnace, characterized in that The device includes: The scanning module is configured to perform material surface scanning based on the burden start signal and the burden end signal of the current burden distribution by using an online laser material surface scanning device of the blast furnace to obtain the material surface coordinates and the burden discharging speed of the current burden distribution, wherein the material surface coordinates include the coordinates of each point located on the top surface of the blast furnace after the current burden distribution, the material surface coordinates of the central part of the blast furnace are calculated by Lagrange interpolation algorithm based on the material surface coordinates scanned by the online laser material surface scanning device, and the material surface coordinates include the material surface coordinates scanned by the online laser material surface scanning device and the material surface coordinates of the central part of the blast furnace; The construction module is configured to determine the first boundary condition based on the material surface coordinates and the characteristics of the blast furnace, and construct an initial grid structure based on the first boundary condition, wherein the initial grid structure is obtained by grid division of the physical model of the blast furnace under the constraint of the first boundary condition, and the characteristics of the blast furnace include the boundary furnace type, the center line of the blast furnace, and the shape of the dead coke pile; The updating module is configured to update the initial grid structure based on the material surface coordinates and the burden discharging speed to obtain the target grid structure of the current burden distribution; The generating module is configured to make the burden in the blast furnace perform solid-liquid-gas three-phase flow and chemical reaction, and obtain the distribution of each grid element in the target grid structure when the multiple burden performance parameters corresponding to each grid node in the target grid structure converge, generate a temperature field distribution map, and the distribution includes the solid temperature, gas temperature, flow rate, and pressure of the region where the grid element is located. The determining module is configured to determine the soft melting zone position of the blast furnace after the burden distribution based on the multiple burden performance parameters corresponding to each grid node in the target grid structure and the target grid structure. The scanning module is specifically configured to: acquire the burden start signal of the current burden distribution when the blast furnace material flow valve is opened; perform material surface scanning based on the delay of the burden start signal by a preset time to obtain the starting material surface coordinates of the current burden distribution; acquire the burden end signal of the current burden distribution when the blast furnace probe is lifted; Delaying for a preset time based on the material end signal to perform material surface scanning to obtain a material surface coordinate of the current material; Determine a material discharge speed of the current material based on the material surface coordinate of the current material, the start material surface coordinate, and a time difference between the material start signal and the material end signal, the material discharge speed being a quotient value obtained by dividing a height difference between the material surface coordinate and the start material surface coordinate by the time difference; The updating module is specifically configured to: Determine a grid parameter of each grid node in the initial grid structure, the grid parameter including a radial distance, a vertical distance, a fluid radial velocity, and a fluid vertical velocity; For any grid node corresponding to the material surface coordinate in the initial grid structure, make the grid node descend to determine a grid parameter of a descended grid node obtained after the descent; Take the descended grid node as a new grid node, repeat the process of the descent of the grid node and the acquisition of the grid parameter of the descended grid node until the grid node or the descended grid node exceeds a physical domain, and obtain the grid parameter of at least one descended grid node corresponding to the grid node; The process of making the grid node descend to determine the grid parameter of the descended grid node obtained after the descent for any grid node corresponding to the material surface coordinate in the initial grid structure includes: Determine a descent time of the grid node; Determine a radial distance of the descended grid node based on the radial distance of the grid node, the fluid radial velocity, and the descent time; Determine a vertical distance of the descended grid node based on the vertical distance of the grid node, the fluid vertical velocity, and the descent time; Determine a grid cell to which the descended grid node belongs in the initial grid structure based on the radial distance and the vertical distance of the descended grid node; Determine four distances between the descended grid node and four vertex grid nodes of the grid cell based on the grid parameter of the grid node, the radial distance and the vertical distance of the descended grid node, and the grid parameter of the four vertex grid nodes of the grid cell; Determine the fluid radial velocity and the fluid vertical velocity of the descended grid node based on the grid parameter of the four vertex grid nodes of the grid cell and the four distances.

6. A computer device, comprising: Comprise: A memory and a processor, which are communicatively connected with each other, and the memory stores computer instructions, and the processor executes the computer instructions to perform the method for determining the position of the soft melting zone of the blast furnace according to any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that, The computer readable storage medium stores computer instructions, and the computer instructions are used to make a computer execute the method for determining the position of the soft melting zone of the blast furnace according to any one claim 1 to 4.

Images