Material distribution method for flat material layer of blast furnace

Abstract

The invention relates to the technical field of blast furnace tiled material layers, and particularly discloses a material distribution method for a blast furnace tiled material layer, which comprises the following steps: dividing furnace charge into large, medium and small particles according to particle diameters, filtering through different pore sieves, then distributing to a blast furnace through a material distribution chute according to a particle size sequence, and dividing a blast furnace material distribution layer into a plurality of areas, a laser radar is used for scanning to obtain sub-region height and wrinkle area information, a height difference value ratio mean value and a cloth wrinkle area ratio are calculated to obtain a cloth flatness coefficient, the cloth flatness coefficient is compared with a threshold value to generate a flatness analysis signal, the cloth flatness coefficient is monitored accordingly, a time-cloth flatness coefficient change curve is drawn, and a material layer fluctuation ratio is calculated. And if the fluctuation ratio of the material layer exceeds a threshold value, generating a tiling analysis signal, and finally, calculating the adjustment value of the rotation angle of the material distribution chute based on the signal and related parameters, so as to monitor and adjust the material distribution condition of the blast furnace, effectively improve the material distribution uniformity and stability of the blast furnace and guarantee stable production of the blast furnace.

Description

Technical Field

[0001] The 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] In the process of burden distribution in a blast furnace, the uniformity and stability of burden distribution are crucial to the efficient and stable operation of the blast furnace. The traditional burden distribution method for blast furnaces has many problems in actual application. On the one hand, it is difficult to control the distribution of burden, and on the other hand, it lacks effective monitoring and adjustment methods, and cannot optimize the burden distribution 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 geometry model, a blast furnace body geometry 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 geometry model; establishing a blast furnace body model through the blast furnace body geometry model to obtain key economic and technical indicators of blast furnace ironmaking; obtaining data through the blast furnace body model, solving various key economic and technical indicators of blast furnace ironmaking function formulas, performing multi-objective optimization on the key economic and technical indicators of blast furnace ironmaking function formulas, and obtaining the optimal key economic and technical indicators of blast furnace ironmaking; 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] In the prior art, the charge control is relatively simple during the charge distribution process, and there is a lack of refined charge distribution strategies. In some traditional blast furnace charge distribution processes, the influence of the distribution of different particle sizes in the blast furnace on the paving of the material layer is not fully considered, and the material layer state cannot be monitored in real time, resulting in low charge distribution uniformity and failure to timely detect changes in the flatness of the material layer during the charge distribution process. The present invention divides the charge particle diameters, distributes the charge to the blast furnace through a charge distribution chute, obtains height information and wrinkle area information by scanning the blast furnace charge distribution layer using a laser radar, calculates the charge distribution flatness coefficient, monitors parameter changes in real time, and optimizes the flatness of the blast furnace charge distribution.

[0005] To this end, the present invention provides a material distribution method 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 materials for a blast furnace flat material layer to solve the above-mentioned problems.

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

[0008] A material distribution method for a blast furnace flat material 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 placed in the blast furnace, it is preferably placed in the blast furnace through the distribution chute in the order of the particle size;

[0010] Step 2: Divide the charge layer in the blast furnace into multiple areas, use laser radar scanning to obtain the height information of the high charge layer in the sub-area and the wrinkle area information of the sub-area, calculate the height difference ratio mean Gd and the charge wrinkle area ratio BI, calculate the charge smoothing coefficient Pz according to the height difference ratio mean Gd and the charge wrinkle area ratio BI, and generate a smoothing analysis signal if the charge smoothing coefficient Pz ≤ the charge smoothing coefficient threshold Pzy;

[0011] Step 3: Based on the flattening analysis signal, monitor the change of the fabric flattening coefficient Pz, take 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, take 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, obtain the vertical distance mean Jl and the 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 the 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 sections, calculate the midpoint coordinates of the target section change curve, connect adjacent midpoints to obtain the target section straight line segment, calculate the slope difference ratio Xl and the angle difference ratio Jc of all the target section 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, and judge 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 best angular speed ratio fitting line, and obtain the ordinate 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 differences in all time periods within 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 in the following manner:

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

[0019] If the fabric flatness coefficient Pz>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 by:

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

[0023] Use laser radar to scan the height information of blast furnace burden 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-areas are summed and averaged to obtain the mean height difference ratio, which is marked as Gd;

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

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

[0029] The sub-region wrinkle area is processed as a ratio with respect 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 cloth wrinkle area ratio, which is marked as B l.

[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-fabric flatness coefficient Pz change curve according to 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 in the following manner:

[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 value J l is obtained by:

[0039] Monitor the change of fabric flatness coefficient Pz in real time, take 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, and 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 flatness coefficient Pz variation curve as calculation points, and calculate the vertical distance between the flatness coefficient reference line and the calculation point;

[0042] Sum and average the vertical distances of all calculation points on the time-fabric flatness 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 material distribution chute are processed as a ratio 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 variation 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 variation curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd in the target section;

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

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

[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 segment;

[0061] In the coordinate system, the slopes of the straight line segments of two adjacent target segments are processed by difference 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 straight line segment slope difference mean is processed by ratio with the straight line segment slope difference threshold to obtain the slope difference ratio, 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 angles of the straight line segments of two adjacent target segments are processed as a difference to obtain the angle difference of the straight line segments;

[0066] The angle differences of all straight line segments are summed and averaged to obtain the average angle difference of the straight line segments;

[0067] The straight line segment angle difference mean is ratioed to the straight line segment angle difference threshold to obtain the angle difference ratio, which 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 changes linearly, 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 a nonlinear variation, and the chute angular speed adjustment ratio Tz is obtained as follows:

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

[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 time periods are summed and averaged to obtain the average fluctuation difference, which is marked as Bdc;

[0078] The difference between the flattening coefficient Pz in each time period and the fabric flattening coefficient threshold Pzy 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 the charge, while ensuring a certain degree of air permeability, which is conducive to the rise and uniform distribution of airflow 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 make adjustments;

[0080] (2) By real-time monitoring of the material distribution smoothness coefficient and drawing a time-material distribution smoothness coefficient change curve, combined with the calculation of the material layer fluctuation ratio, the unevenness and instability factors in the material distribution process can be discovered in time. When the material layer fluctuation ratio is too large, by properly adjusting the rotation angle of the material 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 furnace condition fluctuations caused by charge segregation, reduces production costs, and reduces production interruptions and resource waste caused by uneven material distribution. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0082] Figure 1 This is a flow chart of a material distribution method for flattening a material layer in a blast furnace 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 be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. 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 creative work are 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 materials for a blast furnace flat material 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 placed in the blast furnace, it is preferably placed in the blast furnace through the distribution chute in the order of the particle size;

[0088] Specifically, the charge is divided into three levels according to different particle size ranges, and is divided into 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, and the screened charge is distributed to the blast furnace through the rotation of the distribution chute;

[0089] For example, the charge particle diameter is divided into 20 mm and 50 mm 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 50 mm are marked as large particles, the charge particles with a diameter between 50 mm and 20 mm are marked as medium particles, and the charge particles with a diameter less than 20 mm are marked as small particles.

[0090] Step 2: Divide the charge layer in the blast furnace into multiple areas, use laser radar scanning to obtain the height information of the high charge layer in the sub-area and the wrinkle area information of the sub-area, calculate the height difference ratio mean Gd and the charge wrinkle area ratio BI, calculate the charge smoothing coefficient Pz according to the height difference ratio mean Gd and the charge wrinkle area ratio BI, and generate a smoothing analysis signal if the charge smoothing coefficient Pz ≤ the charge smoothing coefficient threshold Pzy;

[0091] Specifically, the blast furnace burden distribution is divided into a plurality of regions, and the divided regions are marked as sub-regions;

[0092] Use laser radar to scan the height information of blast furnace burden 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-areas 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 unevenness area as the wrinkle area;

[0097] The sub-region wrinkle area is processed as a ratio with respect 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 B l;

[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 material distribution flatness coefficient Pz>the material distribution flatness coefficient threshold Pzy, it indicates that the material distribution of the blast furnace material layer is relatively flat, and the subsequent changes in the material distribution flatness coefficient Pz will continue to be monitored;

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

[0104] Step 3: Based on the flattening analysis signal, monitor the change of the fabric flattening coefficient Pz, take 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, take 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, obtain the vertical distance mean Jl and the 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 the 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 cloth flatness coefficient Pz is monitored in real time, the monitoring time is used as the X-axis, the cloth flatness coefficient Pz is used as the Y-axis, and a time-cloth 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, and 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 flatness coefficient Pz variation curve as calculation points, and calculate the vertical distance between the flatness coefficient reference line and the calculation point;

[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 flatness 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, and mark the distance standard deviation as Bz;

[0112] It should be noted that when the vertical distance mean J l is small, it means that the average deviation between the flatness coefficient reference line and the calculation point is small, that is, the overall change of the fabric 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 fabric flatness coefficient fluctuates relatively slowly in the time dimension, the fabric uniformity is relatively good, and the fabric process is relatively stable without large fluctuations;

[0113] Calculate the material layer fluctuation ratio Bd of the time-material flatness coefficient Pz change curve according to the vertical distance mean J l 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 level of stability of the flat material layer is poor, and a flat 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 change of the material layer fluctuation ratio is still monitored;

[0118] It should be noted that the material layer fluctuation ratio Bd reflects the degree of fluctuation of the cloth flatness coefficient over time. The higher the material layer fluctuation ratio Bd, the more violent the cloth flatness coefficient Pz fluctuates over time, the poorer the uniformity of the cloth, and there are problems such as uneven cloth spreading and unstable cloth process. The material layer fluctuation ratio Bd is an indicator to measure the stability and uniformity of the cloth 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 different hole sieves, 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 smoothness coefficient Pz is calculated according to the height difference ratio mean Gd and the distribution wrinkle area ratio BI, and the distribution smoothness coefficient Pz is compared with the distribution smoothness coefficient threshold Pzy to monitor the distribution smoothness coefficient Pz, and a time-distribution smoothness coefficient Pz change curve is drawn in the coordinate system, and the smoothness coefficient threshold Pzy is used as the smoothness coefficient reference value and a smoothness coefficient reference line is drawn, and the vertical distance between the point on the curve and the smoothness 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 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 material 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 sections, calculate the midpoint coordinates of the target section change curve, connect adjacent midpoints to obtain the target section straight line segment, calculate the slope difference ratio Xl and the angle difference ratio Jc of all target section 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, and judge 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 material distribution chute are processed as a ratio 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 variation 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 variation curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd in the target section;

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

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

[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 segment;

[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), then the coordinates of the midpoint are

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

[0134] The angles of the straight line segments of two adjacent target segments are processed as a difference to obtain the angle difference of the straight line segments;

[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] The angle differences of all straight line segments are summed and averaged to obtain the average angle difference of the straight line segments;

[0137] The straight line segment slope difference mean is processed by ratio with the straight line segment slope difference threshold to obtain the slope difference ratio, and the slope difference ratio is marked as Xl;

[0138] The angle difference mean of the straight line segments is processed by ratio with the angle difference threshold of the straight line segments to obtain the angle difference ratio, which 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 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 best angular speed ratio fitting line, and obtain the ordinate 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 in all time periods during the monitoring period to obtain the fluctuation difference mean Bdc, and calculate the average of the leveling coefficient difference to obtain the leveling coefficient difference mean Pzc. Then, based on the height difference ratio mean Gd, the fluctuation difference mean Bdc, and the leveling coefficient difference mean 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 the 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 time periods are summed and averaged to obtain the average fluctuation difference, which is marked as Bdc;

[0150] Calculate the difference between the flattening coefficient Pz in each time period and the fabric flattening coefficient threshold Pzy to obtain the flattening 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] Based on the height difference ratio mean Gd, the fluctuation difference mean Bdc, and the leveling coefficient difference mean Pzc, the chute angular velocity adjustment ratio Tz is calculated;

[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 material 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, and 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 sections, the midpoint coordinates of the target section change curve are calculated, adjacent midpoints are connected to obtain the target section straight line segment, the slope difference ratio Xl and the angle difference ratio Jc of all the target section 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. l and the angle difference ratio Jc calculate the curve linearity Qx of the change curve, compare the curve linearity Qx with the curve linearity threshold Qxy, determine whether the curve changes linearly or nonlinearly, and 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 best angular speed ratio fitting line, and obtain the ordinate 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 value of the fluctuation difference in all time periods during the monitoring period to obtain the fluctuation difference mean value Bdc, and calculate the average value 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, the chute angular speed adjustment ratio Tz is calculated.

[0156] The above is a detailed description of an embodiment of the present invention, but the content is only a preferred embodiment of the present invention and cannot be considered to limit the scope of implementation 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 material distribution method for a blast furnace flat material layer, characterized in that: include: Step 1: filter the charge through different hole sieves. When the same charge is placed in the blast furnace, the charge is placed in the blast furnace through the distribution chute in the order of particle size. Step 2: Divide the charge layer in the blast furnace into multiple areas, use laser radar scanning to obtain the height information of the high charge layer in the sub-area and the wrinkle area information of the sub-area, calculate the height difference ratio mean Gd and the charge wrinkle area ratio BI, calculate the charge smoothing coefficient Pz according to the height difference ratio mean Gd and the charge wrinkle area ratio BI, and generate a smoothing analysis signal if the charge smoothing coefficient Pz ≤ the charge smoothing coefficient threshold Pzy; Step 3: Based on the flattening analysis signal, monitor the change of the fabric flattening coefficient Pz, draw the time-fabric flattening coefficient Pz change curve, use the flattening coefficient threshold Pzy as the flattening coefficient reference value and draw the flattening coefficient reference line, calculate the vertical distance between the point on the curve and the flattening coefficient reference line, obtain the vertical distance mean Jl and the 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 the distance standard deviation Bz, if the material layer fluctuation ratio Bd≥the material layer fluctuation ratio threshold Bdz, generate a flattening analysis signal; 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 sections, calculate the midpoint coordinates of the target section change curve, connect adjacent midpoints to obtain the target section straight line segment, calculate the slope difference ratio Xl and the angle difference ratio Jc of all the target section 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, and judge 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 material distribution method for a blast furnace flat material layer according to claim 1, characterized in that: 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 in the following manner: Compare the fabric flatness coefficient Pz with the fabric flatness coefficient threshold Pzy; If the fabric flatness coefficient Pz>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.

3. A material distribution method for a blast furnace flat material layer according to claim 2, characterized in that: The height difference ratio mean value Gd is obtained as follows: Divide the blast furnace distribution into multiple areas, and mark the divided areas as sub-areas; Use laser radar to scan the height information of blast furnace burden 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; 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 cloth wrinkle area ratio BI is as follows: Get the area of ​​uneven fabric in each sub-region and mark the unevenness area as the wrinkle area; The sub-region wrinkle area is processed as a ratio with respect 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 cloth wrinkle area ratio, which is marked as Bl.

4. A material distribution method for a blast furnace flat material layer according to claim 1, characterized in that: The material layer fluctuation ratio Bd is obtained as follows: Calculate the material layer fluctuation ratio Bd of the time-fabric flatness coefficient Pz change curve according to the vertical distance mean Jl and standard deviation Bz; By formula: The material layer fluctuation ratio Bd is calculated, where b1 and b2 are preset proportional coefficients; The tiling analysis signal is generated in the following manner: 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.

5. A material distribution method for a blast furnace flat material layer according to claim 4, characterized in that: The vertical distance mean value J1 is obtained as follows: Monitor the change of fabric flatness coefficient Pz in real time, take 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, and 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 flatness coefficient Pz variation curve as calculation points, and calculate the vertical distance between the flatness coefficient reference line and the calculation point; Sum and average the vertical distances of all calculation points on the time-fabric flatness 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.

6. A material distribution method for a blast furnace flat material layer according to claim 1, characterized in that: The method for obtaining 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 to obtain the curve linearity Qx is: The curve linearity Qx is calculated, 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.

7. A material distribution method for a blast furnace flat material layer according to claim 5, 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 material distribution chute are processed as a ratio 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 variation 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 variation curve of the chute rotation angular speed ratio Js-material layer fluctuation ratio Bd in the target section; Connect the two ends of two adjacent midpoints with a straight line to obtain the straight line segment of the target segment; Perform difference processing on the abscissas and ordinates of the two midpoints to obtain the abscissa difference and ordinate difference; 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 segment; In the coordinate system, the slopes of the straight line segments of two adjacent target segments are processed by difference 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 straight line segment slope difference mean is processed by ratio with the straight line segment slope difference threshold to obtain the slope difference ratio, 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 angles of the straight line segments of two adjacent target segments are processed as a difference to obtain the angle difference of the straight line segments; The angle differences of all straight line segments are summed and averaged to obtain the average angle difference of the straight line segments; The straight line segment angle difference mean is ratioed to the straight line segment angle difference threshold to obtain the angle difference ratio, which is marked as Jc.

8. A material distribution method for a blast furnace flat material 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 changes linearly, 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.

9. A material distribution method for a blast furnace flat material 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: Based on the height difference ratio mean Gd, the fluctuation difference mean Bdc, and the leveling coefficient difference mean Pzc, the chute angular velocity adjustment ratio Tz is calculated; The rotation angle adjustment value Tz is constructed as follows; By formula: The rotation angle adjustment value is calculated, where c1, c2, and c3 are preset proportional coefficients.

10. A material distribution method for a blast furnace flat material layer according to claim 9, 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 time periods are summed and averaged to obtain the average fluctuation difference, which is marked as Bdc; Calculate the difference between the smoothing coefficient Pz in each time period and the fabric smoothing coefficient threshold Pzy to obtain the smoothing coefficient difference; within the monitoring period, sum and average the smoothing coefficient differences of all time periods to obtain the mean of the smoothing coefficient differences, and mark the mean of the smoothing coefficient differences as Pzc.

Images