Method, device and equipment for determining position of soft melting zone of blast furnace and storage medium

Abstract

The invention relates to the technical field of blast furnace smelting, and discloses a method, device and equipment for determining the position of a blast furnace soft melting zone and a storage medium, and the method comprises the following steps: carrying out charge level scanning based on a material distribution starting signal and a material distribution ending signal of current material distribution to obtain a charge level coordinate and a charge discharging speed; determining a first boundary condition based on the furnace charge level coordinates and the blast furnace characteristics, and constructing an initial grid structure; the initial grid structure is updated based on the furnace charge level coordinates and the furnace charge discharging speed, and a target grid structure is obtained; enabling the blast furnace burden to react, obtaining the distribution condition of each grid unit in the target grid structure under the condition that a plurality of burden performance parameters of each grid node in the target grid structure converge, and generating a temperature field distribution diagram; and determining the position of the soft melting zone of the blast furnace based on the plurality of furnace charge performance parameters and the temperature field distribution diagram corresponding to each grid node in the target grid structure. According to the method, the accuracy of determining the position of the soft melting zone is improved, and dynamic changes in the blast furnace are reflected.

Description

Technical Field

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

[0002] The soft melting zone of a blast furnace refers to the area in the blast furnace where, when the charge descends to a certain depth, the ore begins to soften and partially melt to form primary slag due to the gradual increase in temperature. 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 soft melting zone is mainly a solid layer for drying, preheating and the beginning of chemical reactions, and a flowing molten layer is formed below the soft melting zone, also known as the dripping zone, in which the iron oxides in the ore continue to be reduced and react with coke to eventually generate molten iron and slag. The location of the blast furnace soft melting zone is crucial to the operation of the blast furnace because it directly affects the utilization of heat energy in the blast furnace, the efficiency of chemical reactions, and the life of the lining. Therefore, how to determine the location of the blast furnace soft melting zone is a problem that needs to be solved.

[0003] At present, infrared thermal imagers and other equipment can be used to scan the blast furnace shell from the outside and infer the location of the soft melting zone based on the temperature distribution. However, infrared detection mainly relies on the heat radiated by the furnace shell, and its ability to penetrate the furnace shell is limited. It can only provide the apparent temperature and it is difficult to accurately reflect the actual temperature distribution deep inside the furnace. As a result, the accuracy of determining the location of the soft melting zone is poor and it cannot 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 soft melting zone of a blast furnace, so as to solve the problem of accurately determining the position of the soft melting zone of a blast furnace.

[0005] In a first aspect, the present invention provides a method for determining a position of a blast furnace soft melting zone, the method comprising: Based on the feeding start signal and feeding end signal of this feeding, the material surface is scanned to obtain the material surface coordinates and the material unloading speed of this feeding, and the material surface coordinates include the coordinates of each point on the top surface of the blast furnace after the feeding is completed; Based on the charge 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, wherein the initial grid structure is obtained by gridding a physical model of the blast furnace under the constraint of the first boundary condition; Based on the charge surface coordinates and charge unloading speed, the initial grid structure is updated to obtain the target grid structure for this charge placement; The charge in the blast furnace is subjected to solid-liquid-gas three-phase flow three-transfer and chemical reaction. When multiple charge performance parameters corresponding to each grid node in the target grid structure converge, the distribution of each grid unit in the target grid structure is obtained, and a temperature field distribution diagram is generated. The distribution includes the solid temperature, gas temperature, flow rate and pressure in the area where the grid unit is located. Based on multiple charge performance parameters and temperature field distribution diagrams corresponding to each grid node in the target grid structure, the position of the blast furnace soft melting zone after this charge distribution is determined.

[0006] The method for determining the position of the soft melting zone of a blast furnace provided in an embodiment of the present invention obtains the charge surface coordinates and charge unloading speed of this charge by scanning the charge surface in real time according to charge start and end signals, thereby updating the initial grid structure constructed based on the first boundary condition through the charge surface coordinates and charge unloading speed to obtain the target grid structure, and simulating the solid-liquid-gas three-phase flow transmission and chemical reaction in the grid structure corresponding to this charge, taking into account the changes of charge with time and space, and more deeply understanding the complex physical and chemical changes in the furnace, thereby generating a temperature field distribution map, and combining the charge performance parameters and the temperature field distribution map to more accurately determine the position of the soft melting zone, providing strong data support for production scheduling and process optimization, improving the accuracy of determining the position of the soft melting zone, and reflecting the dynamic changes inside the blast furnace.

[0007] In an optional implementation, a material surface scan is performed based on a material distribution start signal and a material distribution end signal of this material distribution to obtain the material surface coordinates and the material discharge speed of the furnace material distribution of this material distribution, including: When the material flow valve of the blast furnace is opened, the material distribution start signal of this material distribution is obtained; Scan the material surface based on the cloth start signal to obtain the starting material surface coordinates of this cloth; When the probe of the blast furnace is lifted, the signal for the end of the current material distribution is obtained; Scan the material surface based on the material placement end signal to obtain the material surface coordinates of the current material placement; The material unloading speed of the current burden is determined based on the material surface coordinates of the current burden, the starting material surface coordinates, and the time difference between the burden start signal and the burden end signal.

[0008] The method for determining the position of the soft melting zone of a blast furnace provided in an embodiment of the present invention monitors the start and end signals of the material distribution in real time, scans the material surface coordinates respectively, and determines the material discharge speed through the material surface coordinates obtained by the two scans, thereby providing support for the subsequent determination of the position of the soft melting zone.

[0009] In an optional implementation, based on the charge surface coordinates and the charge unloading speed, the initial grid structure is updated to obtain the target grid structure for this charge distribution, including: Determine the grid parameters of each grid node in the initial grid structure, the grid parameters including radial distance, vertical distance, fluid radial velocity and fluid vertical velocity; For any grid node corresponding to the charge surface coordinates in the initial grid structure, the grid node is lowered, and the grid parameters of the lowered grid node obtained after the lowering are determined; The descending grid node is taken as a new grid node, and the process of descending the grid node and obtaining the grid parameters of the descending grid node is repeated until the grid node or the descending grid node exceeds the physical domain, and the grid parameters of at least one descending grid node corresponding to the grid node are obtained; Based on the grid parameters of at least one descending grid node obtained during the descending process of each grid node corresponding to the charge surface coordinates in the initial grid structure, the grid parameters of each grid node in the initial grid structure are updated to obtain a target grid structure.

[0010] The method for determining the position of the soft melting zone of a blast furnace provided in an embodiment of the present invention determines the grid parameters of each descending grid node during the descending process by descending the grid nodes of the top material surface in the initial grid structure, thereby updating the grid parameters of the grid nodes in the initial grid structure according to the grid parameters of the descending grid nodes obtained by each grid node in the top material surface during the descending process, simulating the descending process of the charge with time and space, reflecting the dynamic changes inside the blast furnace, and obtaining the target grid structure that meets the actual situation of this charge distribution, providing support for the subsequent determination of the position of the soft melting zone.

[0011] In an optional implementation, determining the grid parameters of each grid node in the initial grid structure includes: The coordinates of the charge surface are transformed to obtain the grid parameters of each point in the charge surface coordinates; Based on the mesh parameters of each point in the charge surface coordinates, the mesh parameters of each mesh node in the initial mesh structure are determined.

[0012] The method for determining the position of the soft melting zone of a blast furnace provided in an embodiment of the present invention ensures that the physical properties of each point match the actual charge state by converting the charge surface coordinates into grid parameters, which helps to grasp the furnace conditions more carefully, and based on the grid parameters of each point in the charge surface coordinates, determines the grid parameters of each grid node in the initial grid structure, thereby obtaining an initial grid structure with complete grid node information.

[0013] In an optional embodiment, for any grid node corresponding to the charge surface coordinates in the initial grid structure, the grid node is lowered, and the grid parameters of the lowered grid node obtained after the lowering are determined, including: Determine the fall time of the grid nodes; Determine the radial distance of the descending grid node based on the radial distance of the grid node, the radial velocity of the fluid and the descending time; 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; Based on the radial distance and the vertical distance of the descending grid node, determining the grid unit to which the descending grid node belongs in the initial grid structure; Determine four distances between the descending mesh node and the four vertex mesh nodes of the mesh unit based on the mesh parameters of the mesh node, the radial distance and the vertical distance of the descending mesh node, and the mesh parameters of the four vertex mesh nodes of the mesh unit; Based on the mesh parameters of the four vertex mesh nodes of the mesh unit and the four distances, the fluid radial velocity and the fluid vertical velocity of the descending mesh node are determined.

[0014] The method for determining the position of the soft melting zone of a blast furnace provided in an embodiment of the present invention can simulate the sinking behavior of the charge in the furnace by calculating the grid parameters of the descending grid nodes obtained after the grid nodes are descended, thereby providing intuitive visual and quantitative data support for understanding the dynamic changes in the furnace.

[0015] In an optional embodiment, the charge in the blast furnace is subjected to solid-liquid-gas three-phase flow three-transfer and chemical reaction, and when a plurality of charge performance parameters corresponding to each grid node in the target grid structure converge, the distribution of each grid unit in the target grid structure is obtained, and a temperature field distribution diagram is generated, including: Determine the second boundary condition; The charge in the blast furnace is subjected to solid-liquid-gas three-phase flow three-transfer and chemical reaction for the first time, and a plurality of charge performance parameters corresponding to each grid node in the target grid structure and the distribution of each grid unit are obtained after the first solid-liquid-gas three-phase flow three-transfer and chemical reaction; The above process of subjecting the charge in the blast furnace to solid-liquid-gas three-phase flow, three transfers and chemical reactions is repeated until the multiple charge performance parameters corresponding to each grid node converge under the constraint of the second boundary condition, and the distribution of each grid unit in the target grid structure is obtained to generate a temperature field distribution diagram.

[0016] The method for determining the position of the soft melting zone of a blast furnace provided in an embodiment of the present invention sets a second boundary condition, repeats the three-phase mass transfer and chemical reaction simulation until the charge performance parameters converge, and obtains the distribution of each grid unit, thereby generating a temperature field distribution diagram, which intuitively displays the spatial distribution of the temperature in the furnace, and provides support for the subsequent determination of the position of the soft melting zone.

[0017] In an optional implementation, based on a plurality of charge performance parameters and a temperature field distribution diagram corresponding to each grid node in the target grid structure, determining the position of the blast furnace soft melting zone after this charge distribution includes: Based on multiple charge performance parameters corresponding to each grid node in the target grid structure, high-temperature droplet performance detection is performed to obtain the charge's soft melting start temperature and dripping start temperature; Based on the temperatures of all grid cells in the target grid structure in the temperature field distribution diagram, all grid cells with temperatures between the soft melting start temperature and the dripping start temperature are taken as the positions of the blast furnace soft melting zone after this charging.

[0018] The method for determining the position of the soft melting zone of a blast furnace provided in an embodiment of the present invention can predict the soft melting start temperature and dripping start temperature of the charge by performing high-temperature molten droplet performance detection with multiple charge performance parameters corresponding to each grid node in the target grid structure as the detection environment, thereby determining the grid unit in the target grid structure whose temperature is between the two temperatures as the position of the soft melting zone, thereby achieving accurate determination of the position of the soft melting zone.

[0019] In a second aspect, the present invention provides a device for determining the position of a blast furnace soft melting zone, the device comprising: A scanning module is used to scan the material surface based on the material distribution start signal and the material distribution end signal of this material distribution, and obtain the material surface coordinates and the material discharge speed of the furnace material distribution of this material distribution. The material surface coordinates of the furnace material include the coordinates of each point on the top surface of the blast furnace after the end of this material distribution; A construction module is used to determine a first boundary condition based on the charge surface coordinates and the blast furnace characteristics, and to construct an initial grid structure based on the first boundary condition, wherein the initial grid structure is obtained by gridding the physical model of the blast furnace under the constraint of the first boundary condition; An updating module is used to update the initial grid structure based on the coordinates of the charge surface and the charge unloading speed to obtain the target grid structure for this charge placement; A generation module is used to make the charge in the blast furnace undergo solid-liquid-gas three-phase flow three-transfer and chemical reaction. When multiple charge performance parameters corresponding to each grid node in the target grid structure converge, the distribution of each grid unit in the target grid structure is obtained, and a temperature field distribution diagram is generated. The distribution includes the solid temperature, gas temperature, flow rate and pressure in the area where the grid unit is located. The determination module is used to determine the position of the blast furnace soft melting zone after this charging based on multiple charge performance parameters and temperature field distribution diagrams corresponding to each grid node in the target grid structure.

[0020] In a third aspect, the present invention provides a computer device, comprising: a memory and a processor, the memory and the processor being communicatively connected to each other, the memory storing computer instructions, and the processor executing the method for determining the position of the soft melting zone of a blast furnace according to the first aspect or any corresponding embodiment thereof by executing the computer instructions.

[0021] In a fourth aspect, the present invention provides a computer-readable storage medium having computer instructions stored thereon, the computer instructions being used to enable a computer to execute the method for determining the position of the soft melting zone of a blast furnace according to the first aspect or any corresponding embodiment thereof.

[0022] In a fifth aspect, the present invention provides a computer program product, comprising computer instructions for causing a computer to execute the method for determining the position of a blast furnace soft melting zone according to the first aspect or any corresponding embodiment thereof. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] In order to more clearly illustrate the specific implementation methods of the present invention or the technical solutions in the prior art, the drawings required for use in the specific implementation methods or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some implementation methods of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0024] Figure 1 is a flow chart of a method for determining a position of a blast furnace soft melting zone according to an embodiment of the present invention; Figure 2 is a schematic diagram of an initial grid structure according to an embodiment of the present invention; Figure 3 is a schematic diagram of the distribution of the top material surface of a blast furnace according to an embodiment of the present invention; Figure 4 is a schematic diagram of a grid unit according to an embodiment of the present invention; Figure 5 is a schematic diagram of a target grid structure according to an embodiment of the present invention; Figure 6 is a schematic diagram of a blast furnace soft melting zone according to an embodiment of the present invention; Figure 7 is a flow chart of another method for determining the position of a blast furnace soft melting zone according to an embodiment of the present invention; Figure 8 is a structural block diagram of a device for determining a position of a blast furnace soft melting zone according to an embodiment of the present invention; Fig. 9 It is a schematic diagram of the hardware structure of a computer device according to an embodiment of the present invention. DETAILED DESCRIPTION

[0025] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of the present invention.

[0026] The position of the soft melting zone of a blast furnace is crucial to the operation of the blast furnace because it directly affects the utilization of heat energy, chemical reaction efficiency and life of the lining in the blast furnace. Using infrared thermal imagers and other equipment, the blast furnace shell is scanned from the outside, and the position of the soft melting zone is inferred 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 soft melting zone, and cannot reflect the dynamic changes of the blast furnace. The embodiment of the present invention obtains the charge surface coordinates and charge unloading speed of this charge by scanning the charge surface in real time according to the start and end signals of the charge distribution, updates the initial grid structure, obtains the target grid structure, and simulates the solid-liquid-gas three-phase flow, three transmissions and chemical reactions, taking into account the changes of the charge over time and space, and can more deeply understand the complex physical and chemical changes in the furnace, thereby generating a temperature field distribution map, and comprehensively combining the charge performance parameters and the temperature field distribution map. The position of the soft melting zone can be determined more accurately, which improves the accuracy of determining the position of the soft melting zone and reflects the dynamic changes inside the blast furnace.

[0027] According to an embodiment of the present invention, an embodiment of a method for determining the position of a soft melting zone in a blast furnace is provided. 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 a logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in an order different from that shown here.

[0028] In this embodiment, a method for determining the position of a blast furnace soft melting zone is provided, which can be used in electronic equipment. Figure 1 is a flow chart of a method for determining the position of a blast furnace soft melting zone according to an embodiment of the present invention, such as Figure 1 As shown, the process includes the following steps: Step S101, based on the material distribution start signal and the material distribution end signal of this material distribution, the material surface coordinates and the material discharge speed of the furnace material distribution of this material distribution are obtained. The material surface coordinates of the furnace material distribution include the coordinates of each point on the surface of the top layer of the blast furnace after the end of this material distribution. Specifically, the charge of the blast furnace distribution includes coke and ore. As the coke in the tuyere whirlpool area at the bottom of the blast furnace burns and the molten iron and slag are discharged from the furnace, the blast furnace charge continues to fall. In order to keep the blast furnace working, it is necessary to continue to distribute the material on the top of the blast furnace. Each time the material is distributed, the coke or ore will become the top material surface of the blast furnace. Optionally, the blast furnace online laser material surface scanning equipment can be used to scan the material surface, and only the top material surface of the blast furnace needs to be scanned each time. By scanning the material surface of the current distribution charge, the coordinates of each point on the coke material surface or the ore material surface can be obtained, and the discharge speed of the coke or ore can be obtained through the coordinates.

[0029] Step S102, based on the charge surface coordinates and blast furnace characteristics, determine the first boundary condition, 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. Specifically, the grid structure is obtained by dividing the physical model of the blast furnace into a plurality of grid units, each of which has a non-fixed shape and size, but has four grid nodes as vertices, each of which can be represented by radial distance and vertical distance, and can reflect the radial velocity and vertical velocity of the fluid at the grid node after the charge becomes a fluid in the furnace through the reaction in the blast furnace. Blast furnace characteristics include boundary furnace type, blast furnace centerline, and dead coke pile shape. Figure 2 is a schematic diagram of an initial grid structure according to an embodiment of the present invention. Figure 2 As shown, the top charge surface is obtained by the charge surface coordinates after this charging, the boundary furnace shape is composed of the furnace body, furnace waist, furnace belly, hearth and the depth of the whirlpool area, and the first boundary condition is the boundary formed by the boundary furnace shape, the shape of the dead coke pile, the center line of the blast furnace and the charge surface coordinates. 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 blast furnace shape and the charge conditions of this charging.

[0030] Step S103, based on the charge surface coordinates and charge feeding speed, the initial grid structure is updated to obtain the target grid structure for this charge feeding. Specifically, the initial grid structure is updated according to the charge situation for this charge feeding, so that a target grid structure that is more in line with the actual situation in the current furnace can be obtained while considering the dynamic changes of the charge in the furnace.

[0031] Step S104, the charge in the blast furnace undergoes solid-liquid-gas three-phase flow three-transmission and chemical reaction, and when the multiple charge performance parameters corresponding to each grid node in the target grid structure converge, the distribution of each grid unit in the target grid structure is obtained, and a temperature field distribution diagram is generated, wherein the distribution includes the solid temperature, gas temperature, flow rate and pressure in the area where the grid unit is located. Specifically, when the temperature in the furnace reaches a certain threshold, the iron ore in the charge begins to soften and forms an initial melt. As the temperature further increases, part of the charge melts and drips, forming a blast furnace soft melting zone. Therefore, the temperature factor needs to be considered when determining the location of the blast furnace soft melting zone. During the blast furnace ironmaking process, solid-liquid-gas three-phase flow, three transfers and chemical reactions mainly occur, which cause changes in mass, momentum and heat of the charge in the blast furnace. The above 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 unit in the grid structure corresponding to this distribution is obtained, so that according to the distribution of each grid unit, a temperature field distribution diagram is generated to reflect the temperature of various places inside the blast furnace. Among them, the charge performance parameters include edge solid temperature, gas pressure, silicon composition, top gas temperature and composition. Optionally, other charge performance parameters can also be determined according to actual conditions, and the embodiment of the present invention does not limit this. By generating a temperature field distribution diagram, the blast furnace soft melting zone can be accurately determined based on the temperature of various places inside the blast furnace.

[0032] Step S105, based on the multiple charge performance parameters and temperature field distribution diagrams corresponding to each grid node in the target grid structure, determine the position of the blast furnace soft melting zone after this charge distribution. Specifically, through the charge performance parameters and temperature field distribution diagrams corresponding to each grid node, the situation in the furnace can be fully understood, so that the position of the blast furnace soft melting zone can be determined more accurately.

[0033] The method for determining the position of the soft melting zone of a blast furnace provided in an embodiment of the present invention obtains the charge surface coordinates and charge unloading speed of this charge by scanning the charge surface in real time according to charge start and end signals, thereby updating the initial grid structure constructed based on the first boundary condition through the charge surface coordinates and charge unloading speed to obtain the target grid structure, and simulating the solid-liquid-gas three-phase flow transmission and chemical reaction in the grid structure corresponding to this charge, taking into account the changes of charge with time and space, and more deeply understanding the complex physical and chemical changes in the furnace, thereby generating a temperature field distribution map, and combining the charge performance parameters and the temperature field distribution map to more accurately determine the position of the soft melting zone, providing strong data support for production scheduling and process optimization, and realizing the improvement of the accuracy of determining the position of the soft melting zone, while reflecting the dynamic changes inside the blast furnace.

[0034] In this embodiment, a method for determining the position of a blast furnace soft melting zone is provided, which can be used for the above-mentioned electronic device. The method specifically includes the following steps: Step S201, based on the burden start signal and burden end signal of this burden distribution, the burden surface coordinates and burden unloading speed of this burden distribution are obtained, and the burden surface coordinates include the coordinates of each point on the top surface of the blast furnace after the burden distribution is completed.

[0035] Specifically, the above step S201 includes: Step S2011, when the material flow valve of the blast furnace is opened, a material distribution start signal of this material distribution is obtained. Specifically, when the material flow valve of the blast furnace is opened, material distribution can be carried out, that is, coke or ore can be added to the blast furnace, and at this time, a material distribution start signal is obtained.

[0036] Step S2012, based on the material distribution start signal, the material surface is scanned to obtain the starting material surface coordinates of this material distribution. Specifically, when the material distribution start signal is obtained, in order to make the scanned material surface coordinates more accurate and comprehensive, the material surface scan is performed again after a preset delay time. At this time, the material distribution is stable, and the material surface scan is started by the blast furnace online laser material surface scanning device to obtain the starting material surface coordinates. Optionally, the preset delay time and the material surface scanning time can be set voluntarily, and the embodiment of the present invention does not limit this.

[0037] Step S2013, when the probe of the blast furnace is lifted, a signal for the end of the current feeding is obtained. Specifically, the feeding is ended when the probe of the blast furnace is lifted, that is, the addition of coke or ore to the blast furnace is stopped, and the signal for the end of the feeding is obtained at this time.

[0038] Step S2014, based on the distribution end signal, the material surface is scanned to obtain the material surface coordinates of the furnace charge for this distribution. Specifically, when the distribution end signal is obtained, considering that the dust is large after the distribution is completed, it has an adverse effect on the laser scanning effect of the blast furnace online laser material surface scanning equipment. Therefore, after the distribution end signal is obtained, the material surface is scanned again after a preset time delay to obtain the material surface coordinates of the furnace charge. Optionally, different preset delay times can be set for different charges. For example, the preset delay time after the coke distribution is completed is set to 10 seconds, and the preset delay time after the ore distribution is completed 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 are not limited to this.

[0039] Step S2015, based on the charge surface coordinates of this charge distribution, the starting charge surface coordinates, and the time difference between the charge start signal and the charge end signal, determine the charge unloading speed of this charge distribution. Specifically, the height difference between the charge surface coordinates and the starting charge surface coordinates is divided by the time difference between the two signals to obtain the charge unloading speed. By real-time monitoring of the charge distribution start and end signals, the charge surface coordinates are scanned and obtained respectively, and the charge unloading speed is determined by the charge surface coordinates obtained by the two scans, which provides support for the subsequent determination of the position of the soft melting zone.

[0040] In some optional implementations, the center of the blast furnace is the flame zone, and the laser of the blast furnace online laser material surface scanning device cannot scan the discharge surface, so the coordinates of the center can be calculated using the Lagrange interpolation algorithm. The material surface coordinates directly obtained by scanning and the material surface coordinates of the center obtained by the difference method are spliced ​​into the complete furnace charge material surface coordinates.

[0041] In some optional embodiments, Figure 3 Schematic diagram of the distribution of the top material surface of a blast furnace according to an embodiment of the present invention, such as Figure 3 As shown in the figure, the coordinates of the charge surface are obtained every time the blast furnace is charged. For example, the charge surface coordinates of the ore charged this time are obtained by scanning the charge surface, that is, Figure 3 Upper boundary of the ore shown.

[0042] Step S202: Determine the first boundary condition based on the charge 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 gridding the physical model of the blast furnace under the constraint of the first boundary condition. Figure 1 Step S102 of the illustrated embodiment will not be described in detail here.

[0043] Step S203, based on the charge surface coordinates and charge feeding speed, the initial grid structure is updated to obtain the target grid structure for this charge feeding.

[0044] Specifically, the above step S203 includes: Step S2031, determining the grid parameters of each grid node in the initial grid structure, the grid parameters including radial distance, vertical distance, fluid radial velocity and fluid vertical velocity.

[0045] In some optional implementations, the above step S2031 includes: Step a1, coordinate transformation of the charge surface coordinates to obtain the grid parameters of each point in the charge surface coordinates. Specifically, the charge surface coordinates of this batch are transformed by the Laplace equation, and the grid parameters of each grid node corresponding to the charge surface coordinates in the initial grid structure are calculated by the following formula (1) in the orthogonal coordinate system according to the finite difference method. Optionally, the above process is all prior art and will not be repeated here.

[0046] (1) in, represents the stream function; represents radial distance; Indicates vertical distance.

[0047] Step a2, based on the grid parameters of each point in the charge surface coordinates, determine the grid parameters of each grid node in the initial grid structure. Specifically, since the charge surface coordinates are obtained by scanning the top surface of the blast furnace, the initial grid structure corresponds to the top grid node of the initial grid structure. After determining the network parameters of each grid node at the top layer, the existing technology can be used, such as using the convection trajectory equation, starting from the known top grid nodes, according to the physical properties and motion laws of the charge, to calculate the grid parameters of other grid nodes in the initial grid structure except the top grid nodes, and the specific process will not be repeated here.

[0048] Step S2032, for any grid node corresponding to the charge surface coordinates in the initial grid structure, the grid node is lowered, and the grid parameters of the lowered grid node obtained after the lowering are determined.

[0049] In some optional implementations, the above step S2032 includes: Step b1, determining the descending time of the grid node. Specifically, assuming that any grid node is descended, the descending time is preset. Optionally, in the subsequent repeated descending process, the embodiment of the present invention is described by taking the same descending time for each descent as an example.

[0050] Step b2, based on the radial distance of the grid node, the radial velocity of the fluid and the falling time, determine the radial distance of the descending grid node. Specifically, the radial distance of the descending grid node of the grid node can be determined by the following formula (2).

[0051] (2) in, Indicates that it is located at Line The index of the descending grid node of the column; Indicates that it is located at Line The index of the grid node of the column; Indicates that it is located at Line The radial distance of the descending grid node of the column; Indicates that it is located at Line The radial distance of the grid nodes of the column; Indicates that it is located at Line The radial velocity of the fluid at the grid nodes of the column; Indicates the fall time.

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

[0053] (3) in, Indicates that it is located at Line The index of the descending grid node of the column; Indicates that it is located at Line The index of the grid node of the column; Indicates that it is located at Line The vertical distance of the descending grid node of the column; Indicates that it is located at Line The vertical distance of the grid nodes of the column; Indicates that it is located at Line The vertical velocity of the fluid at the grid nodes of the column; Indicates the fall time.

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

[0055] Step b5, based on the grid parameters of the grid node, the radial distance and the vertical distance of the descending grid node, and the grid parameters of the four vertex grid nodes of the grid unit, determine the four distances between the descending grid node and the four vertex grid nodes of the grid unit. Specifically, the grid unit where point M is located is Q1-Q2-Q3-Q4, and the distances between point M and the four vertices of the grid unit are determined respectively by the following equations (4) to (7), that is, the distances between point M and point Q1, point Q2, point Q3, and point Q4 are determined respectively.

[0056] (4) (5) (6) (7) in, Represents the distance between point M and point Q1; Represents the distance between point M and point Q2; Represents the distance between point M and point Q3; Indicates the distance between point M and point Q4; Indicates that it is located at Line 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 Line The vertical distance of the Q1 point of the column; represents the ordinate of point M located in the Q1-Q2-Q3-Q4 grid cell; Indicates that it is located at Line The radial distance of the Q2 point of the column; Indicates that it is located at Line The vertical distance of point Q2 of the column; Indicates that it is located at Line The radial distance of the Q3 point of the column; Indicates that it is located at Line The vertical distance of the Q3 point of the column; Indicates that it is located at Line The radial distance of the Q4 point of the column; Indicates that it is located at Line The vertical distance of the Q4 point of the column; Step b6, based on the grid parameters of the four vertex grid nodes of the grid unit and the four distances, determine the radial velocity and vertical velocity of the fluid at the descending grid node. Specifically, the interpolation function can be constructed by the following equations (8) and (9) to determine the radial velocity and vertical velocity of the fluid at 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 charge in the furnace can be simulated, providing intuitive visual and quantitative data support for understanding the dynamic changes in the furnace.

[0057] (8) (9) in, Indicates that it is located at Line The radial velocity of the fluid at point M in the column; Indicates that it is located at Line The vertical velocity of the fluid at point M in the column; Indicates that it is located at Line The radial velocity of the fluid at point Q1 of the column; Indicates that it is located at Line The vertical velocity of the fluid at point Q1 of the column; Represents the distance between point M and point Q1; Indicates that it is located at Line The radial velocity of the fluid at point Q2 of the column; Indicates that it is located at Line The vertical velocity of the fluid at point Q2 of the column; Represents the distance between point M and point Q2; Indicates that it is located at Line The radial velocity of the fluid at point Q3 of the column; Indicates that it is located at Line The vertical velocity of the fluid at point Q3 of the column; Represents the distance between point M and point Q3; Indicates that it is located at Line The radial velocity of the fluid at point Q2 of the column; Indicates that it is located at Line The vertical velocity of the fluid at point Q2 of the column; Indicates the distance between point M and point Q4.

[0058] Step S2033, taking the descending grid node as a new grid node, repeating the above 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 obtaining the grid parameters of at least one descending grid node corresponding to the grid node. Specifically, for each grid node at the top layer in the initial grid structure, make it descend in the initial grid structure, determine a descending grid node and its grid parameters after the descending time, then continue to descend the descending grid node as a new grid node, and determine the next descending grid node and its grid parameters, repeat the above process until the grid node or the descending grid node exceeds 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 refer to Figure 2 The grid area shown.

[0059] Step S2034, based on the grid parameters of at least one descending grid node obtained during the descent process of each grid node corresponding to the charge surface coordinates in the initial grid structure, the grid parameters of each grid node in the initial grid structure are updated to obtain a target grid structure. Specifically, for each grid node corresponding to the charge 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 charge at each grid node, the grid parameters of the grid nodes located at the same position as the descending grid nodes in the initial grid structure are correspondingly updated, simulating the descent process of the charge over time and space, being able to reflect the dynamic changes inside the blast furnace, obtaining a target grid structure that meets the actual situation of this charge distribution, and providing support for the subsequent determination of the soft melting zone position. Figure 5 is a schematic diagram of a target grid structure according to an embodiment of the present invention, Figure 5 The target grid structure shown is relative to Figure 2 The initial grid structure shown is more in line with the actual situation of this cloth.

[0060] Step S204, causes the charge in the blast furnace to undergo solid-liquid-gas three-phase flow three-transfer and chemical reaction, and when multiple charge performance parameters corresponding to each grid node in the grid structure corresponding to this charge distribution converge, obtain the distribution of each grid unit in the grid structure corresponding to this charge distribution, and generate a temperature field distribution diagram, the distribution including the solid temperature, gas temperature, flow rate and pressure in the area where the grid unit is located.

[0061] Specifically, the above step S204 includes: Step S2041, determine the second boundary condition. Specifically, the furnace body couple, furnace body static pressure, molten iron composition, furnace top gas temperature and furnace top gas composition are used as the second boundary condition. Among them, the furnace body couple can be obtained by the Fourier heat conduction equation, that is, the following formula (10), which represents the temperature of the solid charge at the furnace wall; the molten iron composition refers to the actual content of silicon in the molten iron.

[0062] (10) in, Indicates heat; Indicates time; represents thermal conductivity; Indicates the furnace area; Indicates furnace temperature; Indicates the length of the furnace body.

[0063] Step S2042, the charge in the blast furnace is subjected to solid-liquid-gas three-phase flow three-transfer and chemical reaction for the first time, and a plurality of charge performance parameters corresponding to each grid node in the target grid structure after the first solid-liquid-gas three-phase flow three-transfer and chemical reaction and the distribution of each grid unit are obtained. Specifically, the charge in the blast furnace is subjected to solid-liquid-gas three-phase flow three-transfer and chemical reaction during the blast furnace ironmaking process. Among them, solid-liquid-gas three-phase flow three-transfer involves the mass, momentum and heat transfer of the charge, and the chemical reaction involves the reduction reaction of iron, carbon and silicon. The charge performance parameters correspond to the above-mentioned second boundary condition, including edge solid temperature, gas pressure, silicon composition, and furnace top gas temperature and composition. The first solid-liquid-gas three-phase flow three-transfer and chemical reaction are caused inside the blast furnace, and the charge performance parameters of each grid node and the distribution of each grid unit after the first reaction are obtained. More specifically, the reduction reaction of iron is shown in the chemical reaction equation shown in the following formula (11), the reduction reaction of carbon is shown in the chemical reaction equation shown in the following formula (12), and the reduction reaction of silicon is shown in the chemical reaction equation shown in the following formula (13).

[0064] (11) (12) (13) in, It is ferric oxide; is carbon monoxide; It is gaseous carbon monoxide; For iron oxide; for carbon dioxide; for hydrogen; For water; for iron; for carbon; For silicon; is silicon dioxide; is gaseous silicon monoxide; each Represents the rate of the corresponding chemical reaction.

[0065] Step S2043, repeat the above process of making the charge in the blast furnace perform solid-liquid-gas three-phase flow three-transfer and chemical reaction, until the multiple charge performance parameters corresponding to each grid node reach convergence under the constraint of the second boundary condition, obtain the distribution of each grid unit in the target grid structure, and generate a temperature field distribution diagram. Specifically, after the first solid-liquid-gas three-phase flow three-transfer and chemical reaction occurs inside the blast furnace, it is judged whether the charge performance parameters of each grid node after the first reaction reach convergence under the constraint of the second boundary condition, that is, it is judged whether the error between the edge solid temperature and the set value of the furnace body electric couple, the error between the set value of the gas pressure and the furnace body static pressure, the error between the set value of the silicon component and the molten iron component, the error between the top gas temperature and the composition in the charge performance parameters and the set value of the 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 parameter of each grid node and the corresponding item in the second boundary condition is within the preset error range, it is considered that convergence is reached. If any grid node does not reach convergence after the first reaction, the second solid-liquid-gas three-phase flow three-transfer and chemical reaction will be carried out, and the performance parameters of the furnace charge after the second reaction will continue to be judged whether they have converged, until each grid node reaches convergence after any solid-liquid-gas three-phase flow three-transfer and chemical reaction, then a temperature field distribution diagram is generated based on the distribution of each grid unit obtained after this reaction. Among them, the temperature field distribution diagram can intuitively display the temperature, flow rate and pressure at each location inside the blast furnace. By setting the second boundary condition, repeating the three-phase mass transfer and chemical reaction simulation until the furnace charge performance parameters converge, the distribution of each grid unit is obtained, and thus a temperature field distribution diagram is generated, which intuitively shows the spatial distribution of the temperature in the furnace, providing support for the subsequent determination of the soft melting zone position.

[0066] Step S205, based on a plurality of furnace charge performance parameters and a temperature field distribution diagram corresponding to each grid node in the grid structure corresponding to the current charge distribution, determine the position of the blast furnace soft melting zone after the current charge distribution.

[0067] Specifically, the above step S205 includes: Step S2051, based on the multiple charge performance parameters corresponding to each grid node in the target grid structure, high-temperature molten droplet performance detection is performed to obtain the softening start temperature and dripping start temperature of the charge. Specifically, high-temperature molten droplet performance detection is a material testing method that focuses on the morphology, melting characteristics, flow behavior and mechanical properties of the molten droplet under high temperature conditions, so it can be used to determine the softening start temperature and dripping start temperature of the charge. The multiple charge performance parameters corresponding to each grid node corresponding to the convergence reached in the above step S2043 are used as environmental parameters in the high-temperature molten droplet performance detection, that is, the detection is performed in the environment corresponding to the charge performance parameters to obtain the softening start temperature and dripping start temperature of the charge.

[0068] Step S2052, based on the temperature of all grid units in the target grid structure in the temperature field distribution map, all grid units with temperatures between the soft melting start temperature and the dripping start temperature are used as the soft melting zone position of the blast furnace after this batching. Specifically, since the temperature field distribution map shows the temperature of each location inside the blast furnace, the grid units with temperatures between the two temperatures in the target grid structure are determined as the soft melting zone position, thereby achieving accurate determination of the soft melting zone position. Figure 6 Schematic diagram of a blast furnace soft melting zone according to an embodiment of the present invention, Figure 6 As shown in the figure, the left side is the center-suppressed blast furnace soft melting zone, and the right side is the ideal blast furnace soft melting zone. The shape of the blast furnace soft melting zone is not fixed. Figure 6 Shown are examples only.

[0069] In some optional embodiments, Figure 7 is a flow chart of another method for determining the position of a blast furnace soft melting zone according to an embodiment of the present invention, such as Figure 7 As shown, firstly, the charge surface after this distribution is scanned to obtain the charge surface coordinates and charge unloading speed. Then, based on the first boundary condition, the charge surface coordinates and the blast furnace characteristics, an initial grid structure is constructed. Then, based on the charge surface coordinates and the charge unloading speed, the initial grid structure is updated to obtain the target grid structure. Then, the charge in the blast furnace is subjected to solid-liquid-gas three-phase flow three-transmission and chemical reaction to obtain multiple charge performance parameters corresponding to each grid node and the distribution of each grid unit. Then, based on the second boundary condition, it is judged whether the multiple charge performance parameters corresponding to each grid node converge. If converged, a temperature field distribution diagram is generated based on the distribution of each grid unit. If not converged, return to the step of subjecting the charge in the blast furnace to solid-liquid-gas three-phase flow three-transmission and chemical reaction. Finally, based on the temperature field distribution diagram, the position of the blast furnace soft melting zone is determined.

[0070] The method for determining the position of the soft melting zone of a blast furnace provided in an embodiment of the present invention obtains the charge surface coordinates and charge unloading speed of this charge by scanning the charge surface in real time according to charge start and end signals, thereby updating the initial grid structure constructed based on the first boundary condition through the charge surface coordinates and charge unloading speed to obtain the target grid structure, and simulating the solid-liquid-gas three-phase flow transmission and chemical reaction in the grid structure corresponding to this charge, taking into account the changes of charge with time and space, and more deeply understanding the complex physical and chemical changes in the furnace, thereby generating a temperature field distribution map, and combining the charge performance parameters and the temperature field distribution map to more accurately determine the position of the soft melting zone, providing strong data support for production scheduling and process optimization, and realizing the improvement of the accuracy of determining the position of the soft melting zone, while reflecting the dynamic changes inside the blast furnace.

[0071] In this embodiment, a device for determining the position of the soft melting zone of a blast furnace is also provided, and the device is used to implement the above-mentioned embodiments and preferred embodiments, and the descriptions that have been made will not be repeated. As used below, the term "module" can implement a combination of software and / or hardware of a predetermined function. Although the device described in the following embodiments is preferably implemented in software, the implementation of hardware, or a combination of software and hardware, is also possible and conceivable.

[0072] This embodiment provides a device for determining the position of a blast furnace soft melting zone. Figure 8 As shown, including: Scanning module 801 is used to scan the material surface based on the material distribution start signal and the material distribution end signal of this material distribution, and obtain the material surface coordinates and material unloading speed of this material distribution. The material surface coordinates include the coordinates of each point on the top surface of the blast furnace after the end of this material distribution.

[0073] The construction module 802 is used to determine the first boundary condition based on the charge surface coordinates and the blast furnace characteristics, and to construct an initial grid structure based on the first boundary condition. The initial grid structure is obtained by gridding the physical model of the blast furnace under the constraint of the first boundary condition.

[0074] The updating module 803 is used to update the initial grid structure based on the coordinates of the charge surface and the charge feeding speed to obtain the target grid structure for this charge feeding.

[0075] Generation module 804 is used to make the charge in the blast furnace undergo solid-liquid-gas three-phase flow three-transfer and chemical reaction. When multiple charge performance parameters corresponding to each grid node in the target grid structure converge, the distribution of each grid unit in the target grid structure is obtained, and a temperature field distribution diagram is generated. The distribution includes the solid temperature, gas temperature, flow rate and pressure of the area where the grid unit is located.

[0076] The determination module 805 is used to determine the position of the soft melting zone of the blast furnace after the current charge distribution based on a plurality of charge performance parameters and a temperature field distribution diagram corresponding to each grid node in the target grid structure.

[0077] In some optional implementations, the scanning module 801 includes: The first acquisition unit is used to acquire a material distribution start signal for this material distribution when the material flow valve of the blast furnace is opened.

[0078] The first scanning unit is used to scan the material surface based on the cloth starting signal to obtain the starting material surface coordinates of this cloth.

[0079] The second acquisition unit is used to acquire a material distribution end signal of this material distribution when the probe of the blast furnace is lifted.

[0080] The second scanning unit is used to scan the material surface based on the material distribution end signal to obtain the material surface coordinates of the furnace material for this distribution.

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

[0082] In some optional implementations, the updating module 803 includes: The second determining unit is used to determine the grid parameters of each grid node in the initial grid structure, and the grid parameters include radial distance, vertical distance, fluid radial velocity and fluid vertical velocity.

[0083] The descending unit is used to descend any grid node corresponding to the charge surface coordinates in the initial grid structure, and determine the grid parameters of the descending grid node obtained after the descending.

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

[0085] The updating unit is used to 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 descending process of each grid node corresponding to the charge surface coordinates in the initial grid structure to obtain a target grid structure.

[0086] In some optional implementations, the second determining unit includes: The conversion subunit is used to convert the coordinates of the charge surface to obtain the grid parameters of each point in the charge surface coordinates.

[0087] The first determination subunit is used to determine the grid parameters of each grid node in the initial grid structure based on the grid parameters of each point in the coordinates of the charge surface.

[0088] In some optional embodiments, the descending unit comprises: The second determining subunit is used to determine the falling time of the grid node.

[0089] The third determining subunit is used to determine the radial distance of the descending grid node based on the radial distance of the grid node, the radial velocity of the fluid and the descending time.

[0090] The fourth determining subunit is used to 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.

[0091] The fifth determining subunit is used to determine the grid unit to which the descending grid node belongs in the initial grid structure based on the radial distance and the vertical distance of the descending grid node.

[0092] The sixth determining subunit is used to determine four distances between the descending mesh node and the four vertex mesh nodes of the mesh unit based on the mesh parameters of the mesh node, the radial distance and the vertical distance of the descending mesh node, and the mesh parameters of the four vertex mesh nodes of the mesh unit.

[0093] The seventh determination subunit is used to determine the fluid radial velocity and the fluid vertical velocity of the descending mesh node based on the mesh parameters of the four vertex mesh nodes of the mesh unit and the four distances.

[0094] In some optional implementations, the generating module 804 includes: The fourth determining unit is used to determine the second boundary condition.

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

[0096] A generation unit is used to repeat the above process of subjecting the charge in the blast furnace to solid-liquid-gas three-phase flow three-transfer and chemical reaction until a plurality of charge performance parameters corresponding to each grid node converge under the constraint of the second boundary condition, obtain the distribution of each grid unit in the target grid structure, and generate a temperature field distribution diagram.

[0097] In some optional implementations, the determination module 805 includes: The detection unit is used to perform high-temperature molten droplet performance detection based on multiple charge performance parameters corresponding to each grid node in the target grid structure to obtain the soft melting start temperature and dripping start temperature of the charge.

[0098] The fifth determination unit is used to determine all grid cells with temperatures between the soft melting start temperature and the dripping start temperature as the soft melting zone position of the blast furnace after this charging based on the temperatures of all grid cells in the target grid structure in the temperature field distribution diagram.

[0099] The further functional description of each of the above modules and units is the same as that of the above corresponding embodiments and will not be repeated here.

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

[0101] The embodiment of the present invention also provides a computer device having the above Figure 8 The device for determining the position of the soft melting zone of a blast furnace is shown.

[0102] See also Fig. 9 , Fig. 9 is a schematic diagram of the structure of a computer device provided by an optional embodiment of the present invention, such as Fig. 9 As shown, the computer device includes: one or more processors 10, a memory 20, and interfaces for connecting various components, including high-speed interfaces and low-speed interfaces. Various components are connected to each other using different buses for communication, and can be installed on a common mainboard or installed in other ways as needed. The processor can process the instructions executed in the computer device, including instructions stored in or on the memory to display the graphical information of the GUI on an external input / output device (such as, a display device coupled to the interface). In some optional 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 (for example, as a server array, a group of blade servers, or a multi-processor system). Fig. 9 A processor 10 is taken as an example.

[0103] The processor 10 may be a central processing unit, a network processor or a combination thereof. The processor 10 may further include a hardware chip. The hardware chip may be a dedicated integrated circuit, a programmable logic device or a combination thereof. The programmable logic device may be a complex programmable logic device, a field programmable gate array, a general purpose array logic or any combination thereof.

[0104] 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 embodiment.

[0105] The memory 20 may include a program storage area and a data storage area, wherein the program storage area may store an operating system, an application 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 also include a non-transient memory, such as at least one disk storage device, a flash memory device, or other non-transient solid-state storage device. In some optional embodiments, the memory 20 may optionally include a memory remotely arranged relative to the processor 10, and these remote memories may be connected to the computer device via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0106] 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 also include a combination of the above types of memory.

[0107] The computer device further comprises a communication interface 30 for the computer device to communicate with other devices or a communication network.

[0108] The embodiment of the present invention also provides a computer-readable storage medium. The method according to the embodiment of the present invention can be implemented in hardware, firmware, or can be implemented as a computer code that can be recorded in a storage medium, or can be implemented as a computer code that is originally stored in a remote storage medium or a non-temporary machine-readable storage medium and will be stored in a local storage medium through a network download, 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 can be a magnetic disk, an optical disk, a read-only storage memory, a random access memory, a flash memory, a hard disk or a solid-state hard disk, etc.; further, the storage medium can also 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. When the software or computer code is accessed and executed by a computer, a processor, or hardware, the method shown in the above embodiment is implemented.

[0109] A part of the present invention may be applied as a computer program product, such as a computer program instruction, which, when executed by a computer, can call or provide the method and / or technical solution according to the present invention through the operation of the computer. Those skilled in the art should understand that the existence of the computer program instruction in a computer-readable medium includes, but is not limited to, a source file, an executable file, an installation package file, etc., and accordingly, the way in which the computer program instruction is executed by the computer includes, but is not limited to: the computer directly executes the instruction, or the computer compiles the instruction and then executes the corresponding compiled program, or the computer reads and executes the instruction, or the computer reads and installs the instruction and then executes the corresponding installed program. Here, the computer-readable medium may be any available computer-readable storage medium or communication medium accessible to the computer.

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

Claims

1. A method for determining the position of a blast furnace soft melting zone, characterized in that: The method comprises: Based on the feeding start signal and the feeding end signal of this feeding, the material surface is scanned to obtain the material surface coordinates and the material unloading speed of the furnace charge of this feeding, wherein the material surface coordinates of the furnace charge include the coordinates of each point on the top surface of the blast furnace after the feeding of this time is completed; Based on the charge surface coordinates and blast furnace characteristics, a first boundary condition is determined, and an initial grid structure is constructed based on the first boundary condition, wherein the initial grid structure is obtained by gridding the physical model of the blast furnace under the constraint of the first boundary condition; Based on the charge surface coordinates and the charge unloading speed, the initial grid structure is updated to obtain a target grid structure for this charge distribution; The charge in the blast furnace is subjected to solid-liquid-gas three-phase flow three-transfer and chemical reaction, and when a plurality of charge performance parameters corresponding to each grid node in the target grid structure converge, the distribution of each grid unit in the target grid structure is obtained, and a temperature field distribution diagram is generated, wherein the distribution includes solid temperature, gas temperature, flow rate and pressure in the area where the grid unit is located; Based on a plurality of furnace charge performance parameters corresponding to each grid node in the target grid structure and the temperature field distribution diagram, the position of the blast furnace soft melting zone after this charging is determined.

2. The method according to claim 1, characterized in that The material surface scanning is performed based on the material distribution start signal and the material distribution end signal of this material distribution to obtain the material surface coordinates and the material discharge speed of the furnace material distribution of this material distribution, including: When the material flow valve of the blast furnace is opened, obtaining a material distribution start signal for this material distribution; Scan the material surface based on the cloth start signal to obtain the starting material surface coordinates of this cloth; When the probe of the blast furnace is lifted, a signal for the end of material distribution is obtained; Scan the material surface based on the material distribution end signal to obtain the material surface coordinates of the furnace material for this distribution; The material unloading speed of the current material distribution is determined based on the material surface coordinates of the current material distribution, the starting material surface coordinates, and the time difference between the material distribution start signal and the material distribution end signal.

3. The method according to claim 1, characterized in that The updating of the initial grid structure based on the charge surface coordinates and the charge unloading speed to obtain the target grid structure for this charge distribution includes: Determining grid parameters of each grid node in the initial grid structure, wherein the grid parameters include radial distance, vertical distance, fluid radial velocity, and fluid vertical velocity; For any grid node corresponding to the charge surface coordinates in the initial grid structure, the grid node is lowered, and a grid parameter of the lowered grid node obtained after the lowering is determined; The descending grid node is used as a new grid node, and the process of descending the grid node and obtaining the grid parameters of the descending grid node is repeated until the grid node or the descending grid node exceeds the physical domain, and the grid parameters of at least one descending grid node corresponding to the grid node are obtained; Based on the grid parameters of at least one descending grid node obtained during the descending process of each grid node corresponding to the 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.

4. The method according to claim 3, characterized in that The step of determining a grid parameter of each grid node in the initial grid structure comprises: Performing coordinate transformation on the charge surface coordinates to obtain grid parameters of each point in the charge surface coordinates; Based on the grid parameters of each point in the charge surface coordinates, the grid parameters of each grid node in the initial grid structure are determined.

5. The method according to claim 3, characterized in that: For any grid node corresponding to the charge surface coordinates in the initial grid structure, the grid node is lowered, and the grid parameters of the lowered grid node obtained after the lowering are determined, including: determining a fall time of the grid node; Determining the radial distance of the descending grid node based on the radial distance of the grid node, the radial velocity of the fluid and the descending time; Determining 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; Determine, based on the radial distance and the vertical distance of the descending grid node, the grid unit to which the descending grid node belongs in the initial grid structure; Determining four distances between the descending mesh node and the four vertex mesh nodes of the mesh unit based on the mesh parameters of the mesh node, the radial distance and the vertical distance of the descending mesh node, and the mesh parameters of the four vertex mesh nodes of the mesh unit; Based on the mesh parameters of the four vertex mesh nodes of the mesh unit and the four distances, the fluid radial velocity and the fluid vertical velocity of the descending mesh node are determined.

6. The method according to claim 1, characterized in that The method of causing the charge in the blast furnace to undergo solid-liquid-gas three-phase flow three-transfer and chemical reaction, and obtaining the distribution of each grid unit in the target grid structure when multiple charge performance parameters corresponding to each grid node in the target grid structure converge, and generating a temperature field distribution diagram, comprises: Determine the second boundary condition; Allowing the charge in the blast furnace to undergo solid-liquid-gas three-phase flow three-transfer and chemical reaction for the first time, and obtaining a plurality of charge performance parameters corresponding to each grid node in the target grid structure and the distribution of each grid unit after the first solid-liquid-gas three-phase flow three-transfer and chemical reaction; The above process of subjecting the charge in the blast furnace to solid-liquid-gas three-phase flow three-transfer and chemical reaction is repeated until a plurality of charge performance parameters corresponding to each grid node converge under the constraint of the second boundary condition, the distribution of each grid unit in the target grid structure is obtained, and the temperature field distribution diagram is generated.

7. The method according to claim 1, characterized in that The determining of the position of the blast furnace soft melting zone after the current charge distribution based on the multiple charge performance parameters corresponding to each grid node in the target grid structure and the temperature field distribution diagram comprises: Based on a plurality of furnace charge performance parameters corresponding to each grid node in the target grid structure, a high-temperature droplet performance test is performed to obtain a remelting start temperature and a dripping start temperature of the furnace charge; Based on the temperatures of all grid units in the target grid structure in the temperature field distribution diagram, all grid units with temperatures between the soft melting start temperature and the dripping start temperature are used as the soft melting zone positions of the blast furnace after this charging.

8. A device for determining the position of a blast furnace soft melting zone, characterized in that: The device comprises: A scanning module is used to scan the material surface based on the material distribution start signal and the material distribution end signal of this material distribution, and obtain the material surface coordinates and the material discharge speed of the furnace charge of this material distribution, wherein the material surface coordinates of the furnace charge include the coordinates of each point on the top surface of the blast furnace after the end of this material distribution; A construction module, for determining a first boundary condition based on the charge surface coordinates and the blast furnace characteristics, and constructing an initial grid structure based on the first boundary condition, wherein the initial grid structure is obtained by gridding the physical model of the blast furnace under the constraint of the first boundary condition; An updating module, used for updating the initial grid structure based on the charge surface coordinates and the charge feeding speed, to obtain a target grid structure for this charge feeding; A generation module is used to make the charge in the blast furnace perform solid-liquid-gas three-phase flow three-transfer and chemical reaction, and obtain the distribution of each grid unit in the target grid structure when multiple charge performance parameters corresponding to each grid node in the target grid structure converge, and generate a temperature field distribution diagram, wherein the distribution includes the solid temperature, gas temperature, flow rate and pressure of the area where the grid unit is located; The determination module is used to determine the position of the blast furnace soft melting zone after this charging based on multiple charge performance parameters corresponding to each grid node in the target grid structure and the temperature field distribution diagram.

9. A computer device, characterized in that: include: A memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the method for determining the position of the soft melting zone of a blast furnace as described in any one of claims 1 to 7 by executing the computer instructions.

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

Images