A material distribution method for blast furnace flat material layer

Abstract

The present invention relates to the technical field of blast furnace flat material layers, and specifically discloses a charging method for blast furnace flat material layers, comprising: dividing furnace charge into large, medium and small particles according to particle diameter, filtering through sieves with different apertures, and then distributing the charge to the blast furnace through a charging chute in order of particle size, dividing the blast furnace charging layer into multiple areas, obtaining sub-area height and wrinkle area information by laser radar scanning, calculating a height difference ratio mean and a charging wrinkle area ratio, thereby obtaining a charging smoothing coefficient, comparing the height difference ratio with a threshold value to generate a smoothing analysis signal, monitoring the charging smoothing coefficient based on the coefficient, plotting a time-charging smoothing coefficient variation curve and calculating a material layer fluctuation ratio, generating a flattening analysis signal if the material layer fluctuation ratio exceeds the threshold, and finally calculating a charging chute rotation angle adjustment value based on the signal and related parameters, so as to realize monitoring and adjustment of the blast furnace charging situation, effectively improve the blast furnace charging uniformity and stability, and ensure stable production of the blast furnace.

Description

Technical Field

[0001] The present invention relates to the technical field of blast furnace flat material layers, and in particular to a material distribution method for blast furnace flat material layers. Background Art

[0002] During the blast furnace charge distribution process, the uniformity and stability of the charge distribution are crucial to the efficient and stable operation of the blast furnace. Traditional blast furnace charge distribution methods have many problems in actual application. On the one hand, it is difficult to control the distribution of the charge. On the other hand, there is a lack of effective monitoring and adjustment methods, and it is impossible to optimize and adjust the charge layer in the blast furnace.

[0003] A Chinese patent application with publication number CN116497166A discloses an intelligent control method for a blast furnace top charging process flow, comprising: establishing a blast furnace top charging geometric model, a blast furnace body geometric model, and a blast furnace top charging layer structure model; processing the blast furnace top charging layer structure model to obtain the result as an input condition for the blast furnace body geometric model; establishing a blast furnace body model through the blast furnace body geometric model to obtain key economic and technical indicators of blast furnace ironmaking; obtaining data through the blast furnace body model to solve the functional expressions of various key economic and technical indicators of blast furnace ironmaking, performing multi-objective optimization on the functional expressions of the key economic and technical indicators of blast furnace ironmaking, and obtaining the optimal key economic and technical indicators of blast furnace ironmaking; and using standard mean square error to evaluate the blast furnace top charging layer structure model to obtain the optimal blast furnace top charging layer structure model.

[0004] Existing technologies offer relatively simple charge distribution control during the charge distribution process, lacking refined distribution strategies. Some traditional blast furnace charge distribution methods fail to fully consider the impact of varying particle size distribution within the furnace on the flatness of the charge bed, resulting in a lack of real-time monitoring of the charge bed's status. This results in poor charge distribution uniformity and an inability to promptly detect changes in the flatness of the charge bed during the distribution process. The present invention classifies charge particles by diameter and distributes them to the blast furnace via a distribution chute. Using a laser radar (LiDAR) to scan the charge bed, the system acquires height and wrinkle area information, calculates the charge flatness coefficient, and monitors parameter changes in real time to optimize the flatness of the charge bed.

[0005] To this end, the present invention provides a method for distributing material for a blast furnace flat material layer. Summary of the Invention

[0006] The object of the present invention is to provide a method for distributing material for a blast furnace flat layer, so as to solve the above-mentioned problems.

[0007] The purpose of the present invention can be achieved through the following technical solutions:

[0008] A method for distributing material for a blast furnace flat layer, comprising:

[0009] Step 1: The charge is divided into large particles, medium particles and small particles according to the size of the particle diameter, and the charge is filtered through different hole sieves. When the same charge is distributed in the blast furnace, it is distributed in the order of particle size through the distribution chute;

[0010] Step 2: Divide the charge layer in the blast furnace into multiple areas, use LiDAR scanning to obtain the height information of the high charge layer in each sub-area and the wrinkle area information of each sub-area, calculate the average height difference ratio Gd and the charge wrinkle area ratio BI, and calculate the charge smoothing coefficient Pz based on the average height difference ratio Gd and the charge wrinkle area ratio BI. If the charge smoothing coefficient Pz is less than or equal to the charge smoothing coefficient threshold Pzy, generate a smoothing analysis signal.

[0011] Step 3: Based on the flattening analysis signal, monitor the change of the fabric flattening coefficient Pz. Use the monitoring time as the X-axis and the fabric flattening coefficient Pz as the Y-axis. Draw a time-fabric flattening coefficient Pz change curve in the coordinate system. Use the flattening coefficient threshold Pzy as the flattening coefficient reference value and draw a flattening coefficient reference line. Calculate the vertical distance between the point on the curve and the flattening coefficient reference line to obtain the vertical distance mean Jl and distance standard deviation Bz. Calculate the material layer fluctuation ratio Bd of the time-fabric flattening coefficient Pz change curve based on the vertical distance mean Jl and distance standard deviation Bz. If the material layer fluctuation ratio Bd ≥ the material layer fluctuation ratio threshold Bdz, generate a flattening analysis signal.

[0012] Step 4: Based on the tiling analysis signal, obtain the speed and angle of the distribution chute rotation, calculate the chute rotation angular speed ratio Js, and draw a change curve with the chute rotation angular speed ratio Js as the X-axis and the material layer fluctuation ratio Bd as the Y-axis. Draw a material layer fluctuation reference line with the material layer fluctuation ratio threshold Bdz. Divide the angular speed ratio Js into multiple target segments, calculate the midpoint coordinates of the target segment change curve, connect adjacent midpoints to obtain target segment straight segments, calculate the slope difference ratio Xl and angle difference ratio Jc of all target segment straight segments, calculate the curve linearity Qx of the change curve based on the slope difference ratio Xl and the angle difference ratio Jc, compare the curve linearity Qx with the curve linearity threshold Qxy to determine whether the curve changes linearly or nonlinearly;

[0013] Step 5. Obtain the linear change state of the change curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd. If it is a linear change, use the least squares method in the rectangular coordinate system to construct the optimal angular speed ratio fitting line, and obtain the vertical coordinate value of the intersection of the fitting line and the material layer fluctuation reference line as the chute angular speed adjustment ratio Tz. If it is a nonlinear change, calculate the average of the fluctuation difference of all time periods in the monitoring period to obtain the fluctuation difference mean Bdc, calculate the average of the leveling coefficient difference to obtain the leveling coefficient difference mean Pzc, and then calculate the chute angular speed adjustment ratio Tz based on the height difference ratio mean Gd, the fluctuation difference mean Bdc, and the leveling coefficient difference mean Pzc.

[0014] As a further technical solution of the present invention: the cloth flatness coefficient Pz is obtained as follows:

[0015] Calculate the fabric flatness coefficient Pz based on the height difference ratio mean Gd and the fabric wrinkle area ratio BI;

[0016] By formula: The fabric flatness coefficient Pz is calculated, where a1 and a2 are preset proportional coefficients;

[0017] The smoothing analysis signal is generated as follows:

[0018] Compare the fabric flatness coefficient Pz with the fabric flatness coefficient threshold Pzy;

[0019] If the fabric flatness coefficient Pz is greater than the fabric flatness coefficient threshold Pzy, continue to monitor the subsequent changes in the fabric flatness coefficient Pz;

[0020] If the fabric flatness coefficient Pz ≤ the fabric flatness coefficient threshold Pzy, a flatness analysis signal is generated.

[0021] As a further technical solution of the present invention: the height difference ratio mean value Gd is obtained as follows:

[0022] Divide the blast furnace distribution area into multiple areas and mark the divided areas as sub-areas;

[0023] Use LiDAR to scan the height information of blast furnace distribution in the sub-area, calculate the height difference between each sub-area and the plane where the sub-area is located, and obtain the regional height difference;

[0024] Based on the height information of each sub-area, the heights of all sub-areas are summed and averaged to obtain the average height of the fabric layer;

[0025] Calculate the ratio of the regional height difference to the average height of the fabric layer to obtain the height difference ratio;

[0026] The height difference ratios of all sub-regions are summed and averaged to obtain the mean height difference ratio, which is marked as Gd;

[0027] The method for obtaining the fabric wrinkle area ratio BI is as follows:

[0028] Get the area of uneven fabric in each sub-region and mark the uneven area as the wrinkle area;

[0029] The sub-region wrinkle area is processed as a ratio to the sub-region area to obtain the sub-region wrinkle area ratio;

[0030] The wrinkle area ratios of all sub-regions are summed and averaged to obtain the fabric wrinkle area ratio, which is marked as B1.

[0031] As a further technical solution of the present invention: the material layer fluctuation ratio Bd is obtained as follows:

[0032] Calculate the material layer fluctuation ratio Bd of the time-material flatness coefficient Pz change curve based on the vertical distance mean J l and standard deviation Bz;

[0033] By formula: The material layer fluctuation ratio Bd is calculated, where b1 and b2 are preset proportional coefficients;

[0034] The tiling analysis signal is generated as follows:

[0035] Compare the material layer fluctuation ratio Bd with the material layer fluctuation ratio threshold Bdz;

[0036] If the material layer fluctuation ratio Bd ≥ the material layer fluctuation ratio threshold Bdz, a tiling analysis signal is generated;

[0037] If the material layer fluctuation ratio Bd is less than the material layer fluctuation ratio threshold Bdz, the change of the material layer fluctuation ratio will continue to be monitored.

[0038] As a further technical solution of the present invention: the vertical distance mean J l is obtained as follows:

[0039] Monitor the change of the fabric flatness coefficient Pz in real time, use the monitoring time as the X-axis and the fabric flatness coefficient Pz as the Y-axis, and draw a time-fabric flatness coefficient Pz change curve in a two-dimensional rectangular coordinate system;

[0040] The fabric flattening coefficient threshold value Pzy is selected as the flattening coefficient reference value. Based on the flattening coefficient reference value, a flattening coefficient reference line parallel to the X-axis is drawn in the rectangular coordinate system.

[0041] Select several points from the time-fabric smoothness coefficient Pz variation curve as calculation points, and calculate the vertical distance between the smoothness coefficient reference line and the calculation points;

[0042] Sum and average the vertical distances of all calculation points on the time-fabric smoothness coefficient Pz change curve to obtain the vertical distance mean J l;

[0043] The standard deviation Bz is obtained as follows:

[0044] The standard deviation of the vertical distances of all calculation points on the time-fabric flatness coefficient Pz change curve is calculated as a further technical solution of the present invention: the method for obtaining the curve linearity Qx is:

[0045] The linearity of the curve of the change curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd is calculated based on the slope difference ratio Xl and the angle difference ratio Jc, and the curve linearity is marked as Qx;

[0046] The method to obtain the curve linearity Qx is: The curve linearity Qx is calculated, where c1 and c2 are preset proportional coefficients;

[0047] The linear change judgment method of the change curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd is:

[0048] Compare the curve linearity Qx with the curve linearity threshold Qxy;

[0049] If the curve linearity Qx ≥ the curve linearity threshold Qxy, it means that the change curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd changes linearly;

[0050] If the curve linearity Qx is less than the curve linearity threshold Qxy, it means that the change curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd is nonlinear.

[0051] As a further technical solution of the present invention: the slope difference ratio Xl is obtained as follows:

[0052] During the monitoring period, obtain the speed and angle of rotation of the material distribution chute;

[0053] The angle and speed of the distribution chute are ratioed to obtain the chute rotation angular speed ratio, which is marked as Js.

[0054] With the chute rotation angular speed ratio Js as the X-axis and the material layer fluctuation ratio Bd as the Y-axis, a curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd is drawn in the rectangular coordinate system;

[0055] Taking the material layer fluctuation ratio threshold value Bdz as the material layer fluctuation reference value, a material layer fluctuation reference line based on the material layer fluctuation reference value is drawn in a rectangular coordinate system;

[0056] Divide the chute rotation angular speed ratio Js into multiple sections, and mark each section as a target section;

[0057] Obtain the midpoint coordinates of the curve of the chute rotation angular speed ratio Js and the material layer fluctuation ratio Bd in the target section;

[0058] Connect the two adjacent midpoints with a straight line to obtain the target segment straight line segment;

[0059] Perform difference processing on the abscissa and ordinate of the two midpoints to obtain the difference between the abscissa and ordinate;

[0060] The difference between the ordinate and the abscissa is processed as a ratio to obtain the slope of the straight line at both ends of the midpoint, that is, the slope of the straight line segment of the target section;

[0061] In the coordinate system, the slopes of the two adjacent target segments are subtracted to obtain the straight line segment slope difference. In the coordinate system, the slope differences of all straight line segments are summed and averaged to obtain the mean of the straight line segment slope differences.

[0062] The slope difference ratio is obtained by ratioing the mean value of the straight line segment slope difference with the straight line segment slope difference threshold value, and the slope difference ratio is marked as Xl;

[0063] The angle difference ratio Jc is obtained as follows:

[0064] Extend the straight line segment of the target segment to intersect with the X-axis, obtain the angle between the straight line segment of the target segment and the X-axis, and obtain the angle between the straight line segment of the target segment;

[0065] The angle difference between the straight line segments of two adjacent target segments is processed to obtain the straight line segment angle difference;

[0066] Sum and average the angle differences of all straight line segments to obtain the mean of the angle differences of the straight line segments;

[0067] The angle difference ratio is obtained by ratioing the mean of the straight line segment angle difference with the straight line segment angle difference threshold, and the angle difference ratio is marked as Jc.

[0068] As a further technical solution of the present invention: the variation curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd is a linear variation, and the chute angular speed adjustment ratio Tz is obtained as follows:

[0069] The variation curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd is linear, and the optimal angular speed ratio fitting line is constructed using the least squares method in the rectangular coordinate system;

[0070] Obtain the coordinates of the intersection of the fitting straight line and the material layer fluctuation reference line, use the ordinate value of the intersection as the chute angular velocity adjustment ratio, and mark the chute angular velocity adjustment ratio as Tz.

[0071] As a further technical solution of the present invention: the variation curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd is nonlinear, and the chute angular speed adjustment ratio Tz is obtained as follows:

[0072] Calculate the chute angular velocity adjustment ratio Tz based on the height difference ratio mean Gd, the fluctuation difference mean Bdc, and the leveling coefficient difference mean Pzc;

[0073] The rotation angle adjustment value Tz is constructed as follows;

[0074] By formula: The rotation angle adjustment value is calculated, where c1, c2, and c3 are preset proportional coefficients.

[0075] As a further technical solution of the present invention: the height difference ratio mean Gd, the fluctuation difference mean Bdc, and the leveling coefficient difference mean Pzc are obtained as follows:

[0076] Calculate the difference between the material layer fluctuation ratio Bd and the material layer fluctuation ratio threshold Bdz in each period during the monitoring period to obtain the fluctuation difference;

[0077] During the monitoring period, the fluctuation differences of all periods are summed and averaged to obtain the average fluctuation difference, which is marked as Bdc;

[0078] The difference between the flattening coefficient Pz and the fabric flattening coefficient threshold Pzy in each time period is calculated to obtain the flattening coefficient difference; within the monitoring period, the flattening coefficient differences of all time periods are summed and averaged to obtain the flattening coefficient difference mean, which is marked as Pzc.

[0079] (1) The charge is divided into three types of particles: large, medium and small according to the particle diameter, and is distributed to the blast furnace through the distribution chute in the order of particle size. This method makes the distribution of the charge in the blast furnace more reasonable. The large-particle charge first enters the bottom of the blast furnace, which can form a relatively stable bottom material layer, providing good support for the subsequent accumulation of charge, while ensuring a certain degree of air permeability, which is conducive to the rise and uniform distribution of air flow in the furnace. The lidar is used to scan the distribution material layer in the blast furnace, obtain the sub-area height information and fold area information, calculate the height difference ratio mean and the distribution material fold area ratio, and then obtain the distribution material smoothness coefficient. Through real-time monitoring and analysis of the distribution material smoothness coefficient, the unevenness of the material layer during the distribution process can be discovered in time. If the distribution material smoothness coefficient exceeds the set threshold, the system will generate a corresponding signal so that timely measures can be taken to adjust it;

[0080] (2) By real-time monitoring of the distribution smoothness coefficient and drawing a time-distribution smoothness coefficient change curve, combined with the calculation of the material layer fluctuation ratio, the unevenness and instability factors in the distribution process can be discovered in a timely manner. When the material layer fluctuation ratio is too large, by appropriately adjusting the rotation angle of the distribution chute, the charge can be distributed more evenly in the blast furnace, reducing the segregation of the charge. This not only reduces the consumption of raw materials, but also reduces the fluctuation of furnace conditions caused by charge segregation, reduces production costs, and reduces production interruptions and resource waste caused by uneven distribution. BRIEF DESCRIPTION OF THE DRAWINGS

[0081] The present invention will be further described below with reference to the accompanying drawings.

[0082] Figure 1 This is a flow chart of a material distribution method for a blast furnace flat material layer according to the present invention;

[0083] Figure 2 The present invention is a flow chart of a method for generating a paving analysis signal in a blast furnace paving material layer distribution method. DETAILED DESCRIPTION

[0084] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only 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 ordinary technicians in this field without making any creative efforts shall fall within the scope of protection of the present invention.

[0085] Example 1

[0086] See also Figure 1 As shown, the present invention is a method for distributing material for a blast furnace flat layer, comprising:

[0087] Step 1: The charge is divided into large particles, medium particles and small particles according to the size of the particle diameter, and the charge is filtered through different hole sieves. When the same charge is distributed in the blast furnace, it is distributed in the order of particle size through the distribution chute;

[0088] Specifically, the charge is divided into three grades according to different particle size ranges, namely large particles, medium particles and small particles according to the size of the particles. Before being distributed to the blast furnace, the charge is filtered and screened through sieves with different apertures. The screened charge is distributed into the blast furnace through the rotation of the distribution chute;

[0089] For example, the charge particle diameter is divided into 20mm and 50mm as the reference value of the particle diameter. The reference value of the particle diameter is set by professionals in this field based on experience. The charge particles with a diameter greater than 50mm are marked as large particles, the charge particles with a diameter between 50mm and 20mm are marked as medium particles, and the charge particles with a diameter less than 20mm are marked as small particles.

[0090] Step 2: Divide the charge layer in the blast furnace into multiple areas, use LiDAR scanning to obtain the height information of the high charge layer in each sub-area and the wrinkle area information of each sub-area, calculate the average height difference ratio Gd and the charge wrinkle area ratio BI, and calculate the charge smoothing coefficient Pz based on the average height difference ratio Gd and the charge wrinkle area ratio BI. If the charge smoothing coefficient Pz is less than or equal to the charge smoothing coefficient threshold Pzy, generate a smoothing analysis signal.

[0091] Specifically, the blast furnace burden distribution system is divided into multiple areas, and the divided areas are marked as sub-areas;

[0092] Use LiDAR to scan the height information of blast furnace distribution in the sub-area, calculate the height difference between each sub-area and the plane where the sub-area is located, and obtain the regional height difference;

[0093] Based on the height information of each sub-area, the heights of all sub-areas are summed and averaged to obtain the average height of the fabric layer;

[0094] Calculate the ratio of the regional height difference to the average height of the fabric layer to obtain the height difference ratio;

[0095] The height difference ratios of all sub-regions are summed and averaged to obtain the mean height difference ratio, which is marked as Gd;

[0096] Get the area of uneven fabric in each sub-region and mark the uneven area as the wrinkle area;

[0097] The sub-region wrinkle area is processed as a ratio to the sub-region area to obtain the sub-region wrinkle area ratio;

[0098] The wrinkle area ratios of all sub-regions are summed and averaged to obtain the fabric wrinkle area ratio, which is marked as B1.

[0099] Calculate the fabric flatness coefficient Pz based on the height difference ratio mean Gd and the fabric wrinkle area ratio BI;

[0100] By formula: The fabric flatness coefficient Pz is calculated, where a1 and a2 are preset proportional coefficients;

[0101] Compare the fabric flatness coefficient Pz with the fabric flatness coefficient threshold Pzy;

[0102] If the distribution flatness coefficient Pz is greater than the distribution flatness coefficient threshold Pzy, it indicates that the distribution of the blast furnace charge layer is relatively flat, and the subsequent changes in the distribution flatness coefficient Pz will continue to be monitored;

[0103] If the distribution flatness coefficient Pz ≤ the distribution flatness coefficient threshold Pzy, it indicates that the distribution flatness of the blast furnace material layer is poor, and further analysis of the flatness of the distribution material layer is required to generate a flatness analysis signal;

[0104] Step 3: Based on the flattening analysis signal, monitor the change of the fabric flattening coefficient Pz. Use the monitoring time as the X-axis and the fabric flattening coefficient Pz as the Y-axis. Draw a time-fabric flattening coefficient Pz change curve in the coordinate system. Use the flattening coefficient threshold Pzy as the flattening coefficient reference value and draw a flattening coefficient reference line. Calculate the vertical distance between the point on the curve and the flattening coefficient reference line to obtain the vertical distance mean Jl and distance standard deviation Bz. Calculate the material layer fluctuation ratio Bd of the time-fabric flattening coefficient Pz change curve based on the vertical distance mean Jl and distance standard deviation Bz. If the material layer fluctuation ratio Bd ≥ the material layer fluctuation ratio threshold Bdz, generate a flattening analysis signal.

[0105] like Figure 2 As shown in the figure, the tiled analysis signal is generated as follows:

[0106] Specifically, the change of the fabric flatness coefficient Pz is monitored in real time, with the monitoring time as the X-axis and the fabric flatness coefficient Pz as the Y-axis, and a time-fabric flatness coefficient Pz change curve is drawn in a two-dimensional rectangular coordinate system;

[0107] The fabric flattening coefficient threshold value Pzy is selected as the flattening coefficient reference value. Based on the flattening coefficient reference value, a flattening coefficient reference line parallel to the X-axis is drawn in the rectangular coordinate system.

[0108] Select several points from the time-fabric smoothness coefficient Pz variation curve as calculation points, and calculate the vertical distance between the smoothness coefficient reference line and the calculation points;

[0109] For example, the vertical distance is marked as d, the coordinates of the calculation point are (3,2), the equation of the leveling coefficient reference line is: Y=6, and the vertical distance can be obtained by the distance formula between the two points: d= | Y-6 | , that is, the distance of the calculated point (3, 2) is 4;

[0110] Sum and average the vertical distances of all calculation points on the time-fabric smoothness coefficient Pz change curve to obtain the vertical distance mean J l;

[0111] Calculate the standard deviation of the vertical distances of all calculation points on the time-fabric flatness coefficient Pz change curve to obtain the distance standard deviation, which is marked as Bz;

[0112] It should be noted that when the mean vertical distance Jl is small, it indicates that the average deviation between the smoothing coefficient reference line and the calculation point is small, that is, the overall change in the flatness coefficient Pz is closer to the flatness coefficient reference line; when the standard deviation Bz is small, the material layer fluctuation ratio Bd will be relatively small, indicating that the fluctuation of the flatness coefficient in the time dimension is relatively gentle, the material distribution uniformity is relatively good, and the distribution process is relatively stable without large fluctuations;

[0113] Calculate the material layer fluctuation ratio Bd of the time-material flatness coefficient Pz change curve based on the vertical distance mean Jl and standard deviation Bz;

[0114] By formula: The material layer fluctuation ratio Bd is calculated, where b1 and b2 are preset proportional coefficients;

[0115] Compare the material layer fluctuation ratio Bd with the material layer fluctuation ratio threshold Bdz;

[0116] If the material layer fluctuation ratio Bd ≥ the material layer fluctuation ratio threshold Bdz, it indicates that the distribution stability of the flat material layer is poor, and a flat distribution analysis signal is generated. The stability of the flat material layer is optimized by adjusting the rotation angle of the distribution chute in the blast furnace.

[0117] If the material layer fluctuation ratio Bd is less than the material layer fluctuation ratio threshold Bdz, it indicates that the material layer is laid smoothly, but the changes in the material layer fluctuation ratio are still monitored;

[0118] It should be noted that the material layer fluctuation ratio Bd reflects the degree of fluctuation of the fabric smoothness coefficient over time. The higher the material layer fluctuation ratio Bd, the more drastic the fluctuation of the fabric smoothness coefficient Pz over time, the poorer the fabric uniformity, and the problems of uneven fabric distribution and unstable fabric process. The material layer fluctuation ratio Bd is an indicator to measure the stability and uniformity of the fabric process.

[0119] The technical solution of this embodiment is as follows: the charge is divided into large particles, medium particles and small particles according to the size of the particle diameter, the charge is filtered through sieves with different holes, and the charge is distributed into the blast furnace through the distribution chute in the order of the particle size, the distribution layer in the blast furnace is divided into multiple areas, and the height information of the high distribution layer of the sub-area and the wrinkle area information of the sub-area are obtained by using laser radar scanning, and the height difference ratio mean Gd and the distribution wrinkle area ratio BI are calculated, and the distribution smoothing coefficient Pz is calculated according to the height difference ratio mean Gd and the distribution wrinkle area ratio BI, and the distribution smoothing coefficient Pz is compared with the distribution smoothing coefficient threshold Pzy to monitor the distribution smoothing coefficient Pz, and a time-distribution distribution smoothing coefficient Pz variation curve is drawn in the coordinate system, and the flatness coefficient threshold Pzy is used as the flatness coefficient reference value and a flatness coefficient reference line is drawn, and the vertical distance between the point on the curve and the flatness coefficient reference line is calculated to obtain the vertical distance mean J l, distance standard deviation Bz, based on the vertical distance mean Jl and the distance standard deviation Bz, calculate the material layer fluctuation ratio Bd of the time-material flatness coefficient Pz change curve. If the material layer fluctuation ratio Bd ≥ the material layer fluctuation ratio threshold Bdz, it indicates that the smoothness of the flat material layer is poor, and a flattening analysis signal is generated.

[0120] Example 2

[0121] Step 4: Based on the tiling analysis signal, obtain the speed and angle of the distribution chute rotation, calculate the chute rotation angular speed ratio Js, and draw a change curve with the chute rotation angular speed ratio Js as the X-axis and the material layer fluctuation ratio Bd as the Y-axis. Draw a material layer fluctuation reference line with the material layer fluctuation ratio threshold Bdz. Divide the angular speed ratio Js into multiple target segments, calculate the midpoint coordinates of the target segment change curve, connect adjacent midpoints to obtain the target segment straight line segment, calculate the slope difference ratio Xl and angle difference ratio Jc of all target segment straight line segments, calculate the curve linearity Qx of the change curve based on the slope difference ratio Xl and the angle difference ratio Jc, compare the curve linearity Qx with the curve linearity threshold Qxy to determine whether the curve changes linearly or nonlinearly.

[0122] Specifically, during the monitoring period, the speed and angle of rotation of the material distribution chute are obtained;

[0123] The angle and speed of the distribution chute are ratioed to obtain the chute rotation angular speed ratio, which is marked as Js.

[0124] With the chute rotation angular speed ratio Js as the X-axis and the material layer fluctuation ratio Bd as the Y-axis, a curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd is drawn in the rectangular coordinate system;

[0125] Taking the material layer fluctuation ratio threshold value Bdz as the material layer fluctuation reference value, a material layer fluctuation reference line based on the material layer fluctuation reference value is drawn in a rectangular coordinate system;

[0126] Divide the chute rotation angular speed ratio Js into multiple sections, and mark each section as a target section;

[0127] Obtain the midpoint coordinates of the curve of the chute rotation angular speed ratio Js and the material layer fluctuation ratio Bd in the target section;

[0128] Connect the two adjacent midpoints with a straight line to obtain the target segment straight line segment;

[0129] Perform difference processing on the abscissa and ordinate of the two midpoints to obtain the difference between the abscissa and ordinate;

[0130] The difference between the ordinate and the abscissa is processed as a ratio to obtain the slope of the straight line at both ends of the midpoint, that is, the slope of the straight line segment of the target section;

[0131] Extend the straight line segment of the target segment to intersect with the X-axis, obtain the angle between the straight line segment of the target segment and the X-axis, and obtain the angle between the straight line segment of the target segment;

[0132] For example, the coordinates of the two endpoints of the target segment are (X1, Y1) and (X2, Y2), and the coordinates of the midpoint are

[0133] In the coordinate system, the slopes of the two adjacent target segments are processed by difference to obtain the slope difference of the straight line segments;

[0134] The angle difference between the straight line segments of two adjacent target segments is processed to obtain the straight line segment angle difference;

[0135] In the coordinate system, the slope differences of all straight line segments are summed and averaged to obtain the mean value of the slope differences of the straight line segments;

[0136] Sum and average the angle differences of all straight line segments to obtain the mean of the angle differences of the straight line segments;

[0137] The slope difference ratio is obtained by ratioing the mean value of the straight line segment slope difference with the straight line segment slope difference threshold value, and the slope difference ratio is marked as Xl;

[0138] The angle difference ratio is obtained by ratioing the mean of the straight line segment angle difference with the straight line segment angle difference threshold, and the angle difference ratio is marked as Jc;

[0139] The linearity of the curve of the change curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd is calculated based on the slope difference ratio Xl and the angle difference ratio Jc, and the curve linearity is marked as Qx;

[0140] The method to obtain the curve linearity Qx is: The curve linearity Qx is calculated, where c1 and c2 are preset proportional coefficients;

[0141] Compare the curve linearity Qx with the curve linearity threshold Qxy;

[0142] If the curve linearity Qx ≥ the curve linearity threshold Qxy, it means that the change curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd changes linearly;

[0143] If the curve linearity Qx is less than the curve linearity threshold Qxy, it means that the change curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd is nonlinear;

[0144] Step 5. Obtain the linear change state of the variation curve of the chute rotation angular speed ratio Js - material layer fluctuation ratio Bd. If it is a linear change, use the least squares method in the rectangular coordinate system to construct the optimal angular speed ratio fitting line. Obtain the vertical coordinate value of the intersection of the fitting line and the material layer fluctuation reference line as the chute angular speed adjustment ratio Tz. If it is a nonlinear change, calculate the average of the fluctuation difference of all time periods in the monitoring period to obtain the fluctuation difference mean value Bdc, and calculate the average of the leveling coefficient difference to obtain the leveling coefficient difference mean value Pzc. Then, based on the height difference ratio mean value Gd, the fluctuation difference mean value Bdc, and the leveling coefficient difference mean value Pzc, calculate the chute angular speed adjustment ratio Tz;

[0145] If the variation curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd changes linearly, the least square method is used in the rectangular coordinate system to construct the optimal angular speed ratio fitting line;

[0146] Obtain the coordinates of the intersection of the fitting straight line and the material layer fluctuation reference line, use the ordinate value of the intersection as the chute angular velocity adjustment ratio, and mark the chute angular velocity adjustment ratio as Tz;

[0147] If the variation curve of chute rotation angular speed ratio Js-material layer fluctuation ratio Bd is nonlinear;

[0148] Calculate the difference between the material layer fluctuation ratio Bd and the material layer fluctuation ratio threshold Bdz in each period during the monitoring period to obtain the fluctuation difference;

[0149] During the monitoring period, the fluctuation differences of all periods are summed and averaged to obtain the average fluctuation difference, which is marked as Bdc;

[0150] Calculate the difference between the flatness coefficient Pz and the fabric flatness coefficient threshold Pzy in each time period to obtain the flatness coefficient difference;

[0151] During the monitoring period, the leveling coefficient differences of all time periods are summed and averaged to obtain the average leveling coefficient difference, which is marked as Pzc;

[0152] Calculate the chute angular velocity adjustment ratio Tz based on the height difference ratio mean Gd, the fluctuation difference mean Bdc, and the leveling coefficient difference mean Pzc;

[0153] The rotation angle adjustment value Tz is constructed as follows;

[0154] By formula: The rotation angle adjustment value is calculated, where c1, c2, and c3 are preset proportional coefficients;

[0155] The technical solution of this embodiment is as follows: based on the tiling analysis signal, the speed and angle of the distribution chute rotation are obtained, the chute rotation angular speed ratio Js is calculated, and a change curve is drawn with the chute rotation angular speed ratio Js as the X-axis and the material layer fluctuation ratio Bd as the Y-axis. A material layer fluctuation reference line is drawn with the material layer fluctuation ratio threshold Bdz. The angular speed ratio Js is divided into multiple target segments, the midpoint coordinates of the target segment change curve are calculated, adjacent midpoints are connected to obtain the target segment straight line segment, the slope difference ratio Xl and the angle difference ratio Jc of all the target segment straight line segments are calculated, and the slope difference ratio Xl and the angle difference ratio Jc are calculated based on the slope difference ratio Xl. The linearity Qx of the change curve is calculated based on the angle difference ratio Jc. This linearity Qx is compared with the linearity threshold Qxy to determine whether the curve changes linearly or nonlinearly. The linear change state of the change curve of the chute rotation angular speed ratio Js and the material layer fluctuation ratio Bd is determined. If the change is linear, the least squares method is used in a rectangular coordinate system to construct the optimal angular speed ratio fitting line. The ordinate value of the intersection of the fitting line and the material layer fluctuation reference line is obtained as the chute angular speed adjustment ratio Tz. If the change is nonlinear, the average of the fluctuation differences for all time periods within the monitoring period is calculated to obtain the fluctuation difference mean Bdc, and the average of the leveling coefficient differences is calculated to obtain the leveling coefficient difference mean Pzc. The chute angular speed adjustment ratio Tz is then calculated based on the height difference ratio mean Gd, the fluctuation difference mean Bdc, and the leveling coefficient difference mean Pzc.

[0156] The above is a detailed description of an embodiment of the present invention. However, the content described is only a preferred embodiment of the present invention and should not be considered to limit the scope of the present invention. All equivalent changes and improvements made within the scope of the present invention should still fall within the scope of the patent coverage of the present invention.

Claims

1. A method for distributing material in a blast furnace flat layer, characterized in that: include: Step 1: filter the charge through sieves with different holes. When distributing the same charge in the blast furnace, prioritize the order of particle size and distribute it through the distribution chute. Step 2: Divide the charge layer in the blast furnace into multiple areas, use LiDAR scanning to obtain the height information of the high charge layer in each sub-area and the wrinkle area information of each sub-area, calculate the average height difference ratio Gd and the charge wrinkle area ratio BI, and calculate the charge smoothing coefficient Pz based on the average height difference ratio Gd and the charge wrinkle area ratio BI. If the charge smoothing coefficient Pz is less than or equal to the charge smoothing coefficient threshold Pzy, generate a smoothing analysis signal. The cloth flatness coefficient Pz is obtained as follows: Calculate the fabric flatness coefficient Pz based on the height difference ratio mean Gd and the fabric wrinkle area ratio BI; By formula: The fabric flatness coefficient Pz is calculated, where a1 and a2 are preset proportional coefficients; The smoothing analysis signal is generated as follows: Compare the fabric flatness coefficient Pz with the fabric flatness coefficient threshold Pzy; If the fabric flatness coefficient Pz is greater than the fabric flatness coefficient threshold Pzy, continue to monitor the subsequent changes in the fabric flatness coefficient Pz; If the fabric flatness coefficient Pz ≤ the fabric flatness coefficient threshold Pzy, a flatness analysis signal is generated; Step 3: Based on the flattening analysis signal, monitor the change of the fabric flattening coefficient Pz, draw a time-fabric flattening coefficient Pz change curve, use the flattening coefficient threshold Pzy as the flattening coefficient reference value and draw a flattening coefficient reference line, calculate the vertical distance between the point on the curve and the flattening coefficient reference line, and obtain the vertical distance mean Jl and distance standard deviation Bz. Based on the vertical distance mean Jl and distance standard deviation Bz, calculate the material layer fluctuation ratio Bd of the time-fabric flattening coefficient Pz change curve. If the material layer fluctuation ratio Bd ≥ the material layer fluctuation ratio threshold Bdz, generate a flattening analysis signal; The material layer fluctuation ratio Bd is obtained as follows: Calculate the material layer fluctuation ratio Bd of the time-fabric smoothness coefficient Pz change curve based on the vertical distance mean Jl and standard deviation Bz; By formula: , calculate the material layer fluctuation ratio Bd, where b1 and b2 are preset proportional coefficients; The tiling analysis signal is generated as follows: Compare the material layer fluctuation ratio Bd with the material layer fluctuation ratio threshold Bdz; If the material layer fluctuation ratio Bd ≥ the material layer fluctuation ratio threshold Bdz, a tiling analysis signal is generated; If the material layer fluctuation ratio Bd is less than the material layer fluctuation ratio threshold Bdz, the change of the material layer fluctuation ratio will continue to be monitored; Step 4: Based on the tiling analysis signal, obtain the speed and angle of the distribution chute rotation, calculate the chute rotation angular speed ratio Js, and draw a change curve with the chute rotation angular speed ratio Js as the X-axis and the material layer fluctuation ratio Bd as the Y-axis. Draw a material layer fluctuation reference line with the material layer fluctuation ratio threshold Bdz. Divide the angular speed ratio Js into multiple target segments, calculate the midpoint coordinates of the target segment change curve, connect adjacent midpoints to obtain the target segment straight line segment, calculate the slope difference ratio Xl and angle difference ratio Jc of all target segment straight line segments, calculate the curve linearity Qx of the change curve based on the slope difference ratio Xl and the angle difference ratio Jc, compare the curve linearity Qx with the curve linearity threshold Qxy to determine whether the curve changes linearly or nonlinearly; Step 5: Obtain the linear change state of the change curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd, and calculate the chute angular speed adjustment ratio Tz under different states.

2. A method for distributing material for a blast furnace flat layer according to claim 1, characterized in that: The height difference ratio mean value Gd is obtained as follows: Divide the blast furnace charge layer into multiple areas, and mark the divided areas as sub-areas; Use LiDAR to scan the height information of the blast furnace charge layer in the sub-area, calculate the height difference between each sub-area and the plane where the sub-area is located, and obtain the regional height difference; Based on the height information of each sub-area, the heights of all sub-areas are summed and averaged to obtain the average height of the fabric layer; Calculate the ratio of the regional height difference to the average height of the fabric layer to obtain the height difference ratio; The height difference ratios of all sub-areas are summed and averaged to obtain the mean height difference ratio, which is marked as Gd; The method for obtaining the fabric wrinkle area ratio BI is as follows: Get the area of uneven fabric in each sub-region and mark the uneven area as the wrinkle area; The sub-region wrinkle area is processed as a ratio to the sub-region area to obtain the sub-region wrinkle area ratio; The wrinkle area ratios of all sub-regions are summed and averaged to obtain the fabric wrinkle area ratio, which is marked as Bl.

3. The method for distributing material for a blast furnace flat layer according to claim 1, characterized in that: The vertical distance mean value J1 is obtained as follows: Monitor the change of the fabric flatness coefficient Pz in real time, use the monitoring time as the X-axis and the fabric flatness coefficient Pz as the Y-axis, and draw a time-fabric flatness coefficient Pz change curve in a two-dimensional rectangular coordinate system; The fabric flattening coefficient threshold value Pzy is selected as the flattening coefficient reference value. Based on the flattening coefficient reference value, a flattening coefficient reference line parallel to the X-axis is drawn in the rectangular coordinate system. Select several points from the time-fabric smoothness coefficient Pz variation curve as calculation points, and calculate the vertical distance between the smoothness coefficient reference line and the calculation points; Sum and average the vertical distances of all calculation points on the time-fabric smoothness coefficient Pz change curve to obtain the vertical distance mean Jl; The standard deviation Bz is obtained as follows: Calculate the standard deviation of the vertical distances of all calculation points on the time-fabric flatness coefficient Pz change curve to obtain the distance standard deviation, and mark the distance standard deviation as Bz.

4. The method for distributing material for a blast furnace flat layer according to claim 1, characterized in that: The method for calculating the curve linearity Qx is: The linearity of the curve of the change curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd is calculated based on the slope difference ratio Xl and the angle difference ratio Jc, and the curve linearity is marked as Qx; The method for calculating the curve linearity Qx is: , calculate the curve linearity Qx, where c1 and c2 are preset proportional coefficients; The linear change judgment method of the change curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd is: Compare the curve linearity Qx with the curve linearity threshold Qxy; If the curve linearity Qx ≥ the curve linearity threshold Qxy, it means that the change curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd changes linearly; If the curve linearity Qx is less than the curve linearity threshold Qxy, it means that the change curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd is nonlinear.

5. The method for distributing material for a blast furnace flat layer according to claim 1, characterized in that: The slope difference ratio X1 is obtained as follows: During the monitoring period, obtain the speed and angle of rotation of the material distribution chute; The angle and speed of the distribution chute are ratioed to obtain the chute rotation angular speed ratio, which is marked as Js. With the chute rotation angular speed ratio Js as the X-axis and the material layer fluctuation ratio Bd as the Y-axis, a curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd is drawn in the rectangular coordinate system; Taking the material layer fluctuation ratio threshold value Bdz as the material layer fluctuation reference value, a material layer fluctuation reference line based on the material layer fluctuation reference value is drawn in a rectangular coordinate system; Divide the chute rotation angular speed ratio Js into multiple sections, and mark each section as a target section; Obtain the midpoint coordinates of the curve of the chute rotation angular speed ratio Js and the material layer fluctuation ratio Bd in the target section; Connect the two adjacent midpoints with a straight line to obtain the target segment straight line segment; Perform difference processing on the abscissa and ordinate of the two midpoints to obtain the difference between the abscissa and ordinate; The difference between the ordinate and the abscissa is processed as a ratio to obtain the slope of the straight line at both ends of the midpoint, that is, the slope of the straight line segment of the target section; In the coordinate system, the slopes of the two adjacent target segments are subtracted to obtain the straight line segment slope difference. In the coordinate system, the slope differences of all straight line segments are summed and averaged to obtain the mean of the straight line segment slope differences. The slope difference ratio is obtained by ratioing the mean of the straight line segment slope differences with the straight line segment slope difference threshold, and the slope difference ratio is marked as X1; The angle difference ratio Jc is obtained as follows: Extend the straight line segment of the target segment to intersect with the X-axis, obtain the angle between the straight line segment of the target segment and the X-axis, and obtain the angle between the straight line segment of the target segment; The angle difference between the straight line segments of two adjacent target segments is processed to obtain the straight line segment angle difference; Sum and average the angle differences of all straight line segments to obtain the mean of the angle differences of the straight line segments; The angle difference ratio is obtained by ratioing the mean of the straight line segment angle difference with the straight line segment angle difference threshold, and the angle difference ratio is marked as Jc.

6. The method for distributing material for a blast furnace flat layer according to claim 1, characterized in that: The variation curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd is a linear variation, and the chute angular speed adjustment ratio Tz is obtained as follows: The variation curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd is linear, and the optimal angular speed ratio fitting line is constructed using the least squares method in the rectangular coordinate system; Obtain the coordinates of the intersection of the fitting straight line and the material layer fluctuation reference line, use the ordinate value of the intersection as the chute angular velocity adjustment ratio, and mark the chute angular velocity adjustment ratio as Tz.

7. The method for distributing material for a blast furnace flat layer according to claim 1, characterized in that: The variation curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd is nonlinear, and the chute angular speed adjustment ratio Tz is obtained as follows: Calculate the chute angular velocity adjustment ratio Tz based on the height difference ratio mean Gd, the fluctuation difference mean Bdc, and the leveling coefficient difference mean Pzc; The rotation angle adjustment value Tz is constructed as follows; By formula: , calculate the rotation angle adjustment value, where c1, c2, and c3 are preset proportional coefficients.

8. A method for distributing material for a blast furnace flat layer according to claim 7, characterized in that: The method for obtaining the height difference ratio mean Gd, the fluctuation difference mean Bdc, and the leveling coefficient difference mean Pzc is as follows: Calculate the difference between the material layer fluctuation ratio Bd and the material layer fluctuation ratio threshold Bdz in each period during the monitoring period to obtain the fluctuation difference; During the monitoring period, the fluctuation differences of all periods are summed and averaged to obtain the average fluctuation difference, which is marked as Bdc; Calculate the difference between the flatness coefficient Pz in each time period and the fabric flatness coefficient threshold Pzy to obtain the flatness coefficient difference; within the monitoring period, sum and average the flatness coefficient differences of all time periods to obtain the average flatness coefficient difference, and mark the average flatness coefficient difference as Pzc.

Images