阅读说明：本技术 一种高炉料面形状测量方法、终端设备及存储介质 (Blast furnace burden surface shape measuring method, terminal equipment and storage medium ) 是由 严晗 叶理德 欧燕 秦涔 崔伟 方明新 闫朝付 于 2021-07-23 设计创作，主要内容包括：本发明涉及一种高炉料面形状测量方法、终端设备及存储介质,该方法中包括：S1：通过安装于高炉炉顶的测距雷达对高炉内料面的位置进行测量；S2：对测量得到的料面的雷达测距数据进行预处理,剔除其内的异常数据；S3：将雷达测距数据转换为对应的XY直角坐标系中的坐标数据；S4：通过克里金插值法,使用坐标数据对料面形状进行拟合。本发明相对应现有技术,不需要复杂的传感器布置,利于安装维护,且相对于一般单点雷达测量,具有更好的测量效果。(The invention relates to a blast furnace burden surface shape measuring method, terminal equipment and a storage medium, wherein the method comprises the following steps: s1: measuring the position of the charge level in the blast furnace by a ranging radar arranged on the top of the blast furnace; s2: preprocessing radar ranging data of the measured charge level, and removing abnormal data in the radar ranging data; s3: converting radar ranging data into coordinate data in a corresponding XY rectangular coordinate system; s4: and fitting the shape of the charge level by using coordinate data through a Krigin interpolation method. Compared with the prior art, the invention does not need complex sensor arrangement, is beneficial to installation and maintenance, and has better measuring effect compared with the common single-point radar measurement.)

1. A blast furnace burden surface shape measuring method is characterized by comprising the following steps:

s1: measuring the position of the charge level in the blast furnace by a ranging radar arranged on the top of the blast furnace;

s2: preprocessing radar ranging data of the measured charge level, and removing abnormal data in the radar ranging data;

s3: converting radar ranging data into coordinate data in a corresponding XY rectangular coordinate system;

s4: and fitting the shape of the charge level by using coordinate data through a Krigin interpolation method.

2. The blast furnace burden surface shape measuring method according to claim 1, characterized in that: the method for measuring the position of the charge level in the blast furnace by the ranging radar comprises the following steps: and controlling the rotation of the ranging radar so that the ranging radar measures the distance between the charge level in the radius range of the furnace throat and the ranging radar.

3. The blast furnace burden surface shape measuring method according to claim 1, characterized in that: the range radar includes two, and both install in furnace throat top and the both sides of same height.

4. The blast furnace burden surface shape measuring method according to claim 1, characterized in that: the radar ranging data comprises the distance between the charge level and the ranging radar and the included angle between the direction of the ranging radar and the horizontal plane.

5. The blast furnace burden surface shape measuring method according to claim 1, characterized in that: the pretreatment method comprises the following steps: setting the range of radar ranging data according to the installation position of the ranging radar, thereby eliminating the radar ranging data which is not in the range; and according to the set echo intensity threshold value, eliminating radar ranging data with the echo intensity smaller than the echo intensity threshold value.

6. The blast furnace burden surface shape measuring method according to claim 1, characterized in that: the calculation formula for converting the radar ranging data into the coordinate data in the corresponding XY rectangular coordinate system is as follows:

y₀＝d₀ sinα-h_r

wherein x is₀、y₀Are coordinate values of X-axis and Y-axis in the coordinate data, r represents furnace throat radius, h_rIndicates the height h from the installation position of the range radar to the zero stockline_lThe height of the boundary on the furnace throat from a zero material line is shown, gamma is the included angle between the furnace wall on the furnace throat and the horizontal plane, alpha is the included angle between the direction of the range radar and the horizontal plane, and d₀Representing the distance between the level and the range radar.

7. The blast furnace burden surface shape measuring method according to claim 1, characterized in that: step S4 specifically includes the following steps:

s41: carrying out radial segmentation on the charge level, setting an X-axis coordinate range corresponding to each segment, and setting each segmentation point as a point to be inserted;

s42: calculating any two coordinate data x in the coordinate data_iAnd x_jA distance d between_ij＝x_i-x_jGrouping all the calculated distances according to the distance; where i is 0,1, …, n, j is 0,1, …, n, n denotes the total number of coordinate data, i and j denote the serial numbers of the coordinate data, x_iAnd x_jRespectively representing the coordinates of the ith and the j th coordinate data in the X-axis direction;

s43: the mean value d of all the distances in each group is calculated and the estimated value γ of the corresponding half-variance function of each group is calculated according to the following formula^*(d)：

Wherein, y (x)_i) And y (x)_i+ d) respectively representing X as the coordinate in the X direction_iAnd (x)_iThe Y-axis direction coordinate corresponding to + d);

s44: the half-variance function model is constructed as follows: γ (d) ═ c (1-e)^-d/r) Wherein c and r are both parameters in the half variance function model, gamma (d) represents the half variance operation, and the estimated value gamma of each group of corresponding half variance functions is determined according to the half variance function model^*(d) Fitting to obtain a simulationA combined half-variance function model;

s45: the following linear system of equations was constructed:

according to the fitted half variance function model and the formula r_ij＝γ(d_ij) Calculating any two coordinate data x_iAnd x_jHalf-variance of an attribute r_ijAnd further obtaining a coefficient matrix in the linear equation system:

s46: for each point x to be inserted₀Calculating x₀With arbitrary coordinate data x_iHalf-variance of an attribute r_i0Further obtain the point x to be inserted₀Right-hand array in the corresponding linear system of equations:

s47: solving the linear equation set according to the coefficient matrix and the right-side array in the linear equation set to obtain the point x to be inserted₀Corresponding weight array:

s48: for each point x to be inserted₀According to the point x to be inserted₀Corresponding weight array, calculating the point x to be interpolated by₀Is estimated value of

Wherein, y_iRepresenting x in coordinate data_iThe corresponding Y-axis direction coordinates.

8. A blast furnace burden surface shape measuring terminal device is characterized in that: comprising a processor, a memory and a computer program stored in said memory and running on said processor, said processor implementing the steps of the method according to any one of claims 1 to 7 when executing said computer program.

9. A computer-readable storage medium storing a computer program, characterized in that: the computer program when executed by a processor implementing the steps of the method as claimed in any one of claims 1 to 7.

Technical Field

The invention relates to the technical field of blast furnace ironmaking detection, in particular to a blast furnace burden surface shape measuring method, terminal equipment and a storage medium.

Background

The blast furnace is used as a black box-like closed container, the internal environment is complex, and information is difficult to directly obtain. In the traditional technology, a mechanical stock rod is generally used for acquiring data about the charge level of a blast furnace, but the measurement range is small, the data error is large, and the distribution condition of the charge level in the furnace is difficult to analyze. Then, the blast furnace charge level is researched and developed to be monitored by an infrared imaging or other imaging means, but the method is greatly interfered due to large temperature change and much dust on the upper part of the charge level, and the shape of the charge level is difficult to accurately obtain. In some measuring methods, a certain number of distance measuring radars are installed on holes at multiple positions on the furnace top to obtain the approximate shape of the charge level, but the method has the defects that the number of sensors is large, a sensing system is complex, the installation and maintenance of multiple devices bring difficulties to field operators, the real-time performance of the system operation is difficult to guarantee, or the sensing number is too small to reflect the charge level shape information.

Disclosure of Invention

In order to solve the above problems, the present invention provides a blast furnace burden surface shape measuring method, a terminal device, and a storage medium.

The specific scheme is as follows:

a blast furnace burden surface shape measuring method comprises the following steps:

s1: measuring the position of the charge level in the blast furnace by a ranging radar arranged on the top of the blast furnace;

s2: preprocessing radar ranging data of the measured charge level, and removing abnormal data in the radar ranging data;

s3: converting radar ranging data into coordinate data in a corresponding XY rectangular coordinate system;

s4: and fitting the shape of the charge level by using coordinate data through a Krigin interpolation method.

Further, the method for measuring the position of the charge level in the blast furnace by the ranging radar comprises the following steps: and controlling the rotation of the ranging radar so that the ranging radar measures the distance between the charge level in the radius range of the furnace throat and the ranging radar.

Further, the range radar includes two, and both install in furnace throat top and the both sides of same height.

Further, the radar ranging data comprise the distance between the charge level and the ranging radar and the included angle between the direction of the ranging radar and the horizontal plane.

Further, the pretreatment method comprises the following steps: setting the range of radar ranging data according to the installation position of the ranging radar, thereby eliminating the radar ranging data which is not in the range; and according to the set echo intensity threshold value, eliminating radar ranging data with the echo intensity smaller than the echo intensity threshold value.

Further, a calculation formula for converting the radar ranging data into coordinate data in a corresponding XY rectangular coordinate system is as follows:

y₀＝d₀sinα-h_r

wherein x is₀、y₀Are coordinate values of X-axis and Y-axis in the coordinate data, r represents furnace throat radius, h_rIndicates the height h from the installation position of the range radar to the zero stockline_lThe height of the boundary on the furnace throat from a zero material line is shown, gamma is the included angle between the furnace wall on the furnace throat and the horizontal plane, alpha is the included angle between the direction of the range radar and the horizontal plane, and d₀Representing the distance between the level and the range radar.

Further, step S4 specifically includes the following steps:

s41: carrying out radial segmentation on the charge level, setting an X-axis coordinate range corresponding to each segment, and setting each segmentation point as a point to be inserted;

s42: calculating any two coordinate data x in the coordinate data_iAnd x_jA distance d between_ij＝x_i-x_jGrouping all the calculated distances according to the distance; where i is 0,1, …, n, j is 0,1, …, n, n denotes the total number of coordinate data, i and j denote the serial numbers of the coordinate data, x_iAnd x_jRespectively represent the ith and the j th coordinate dataCoordinates in the X-axis direction;

s43: the mean value d of all the distances in each group is calculated and the estimated value γ of the corresponding half-variance function of each group is calculated according to the following formula^*(d)：

Wherein, y (x)_i) And y (x)_i+ d) respectively representing X as the coordinate in the X direction_iAnd (x)_iThe Y-axis direction coordinate corresponding to + d);

s44: the half-variance function model is constructed as follows: γ (d) ═ c (1-e)^-d/r) Wherein c and r are both parameters in the half variance function model, gamma (d) represents the half variance operation, and the estimated value gamma of each group of corresponding half variance functions is determined according to the half variance function model^*(d) Fitting to obtain a fitted half variance function model;

s45: the following linear system of equations was constructed:

according to the fitted half variance function model and the formula r_ij＝γ(d_ij) Calculating any two coordinate data x_iAnd x_jHalf-variance of an attribute r_ijAnd further obtaining a coefficient matrix in the linear equation system:

s46: for each point x to be inserted₀Calculating x₀With arbitrary coordinate data x_iHalf-variance of an attribute r_i0Further obtain the point x to be inserted₀Right-hand array in the corresponding linear system of equations:

s47: solving the linear equation set according to the coefficient matrix and the right-side array in the linear equation setTo obtain the point x to be inserted₀Corresponding weight array:

s48: for each point x to be inserted₀According to the point x to be inserted₀Corresponding weight array, calculating the point x to be interpolated by₀Is estimated value of

Wherein, y_iRepresenting x in coordinate data_iThe corresponding Y-axis direction coordinates.

A blast furnace burden surface shape measuring terminal device comprises a processor, a memory and a computer program stored in the memory and capable of running on the processor, wherein the processor executes the computer program to realize the steps of the method of the embodiment of the invention.

A computer-readable storage medium, in which a computer program is stored, which, when being executed by a processor, carries out the steps of the method as described above for an embodiment of the invention.

The technical scheme is adopted, and the device does not need complex sensor arrangement corresponding to the prior art, is beneficial to installation and maintenance, and has better measuring effect compared with the common single-point radar measurement.

Drawings

Fig. 1 is a flowchart illustrating a first embodiment of the present invention.

Fig. 2 is a schematic diagram illustrating a ranging radar measurement process in this embodiment.

Detailed Description

To further illustrate the various embodiments, the invention provides the accompanying drawings. The accompanying drawings, which are incorporated in and constitute a part of this disclosure, illustrate embodiments of the invention and, together with the description, serve to explain the principles of the embodiments. Those skilled in the art will appreciate still other possible embodiments and advantages of the present invention with reference to these figures.

The invention will now be further described with reference to the accompanying drawings and detailed description.

The first embodiment is as follows:

the embodiment of the invention provides a blast furnace burden surface shape measuring method, as shown in fig. 1, the method comprises the following steps:

s1: and measuring the position of the charge level in the blast furnace by a distance measuring radar arranged on the top of the blast furnace.

The embodiment of the range radar adopts a millimeter wave radar, and in other embodiments, a person skilled in the art may also adopt other range radars, which is not limited herein.

As shown in fig. 2, the method for measuring the position of the burden surface in the blast furnace by the ranging radar in this embodiment includes: and controlling the rotation of the ranging radar so that the ranging radar measures the distance between the charge level in the radius range of the furnace throat and the ranging radar.

In this embodiment, it is preferable to set the range radar to include two, R₁、R₂The two are arranged above the furnace throat and on two sides with the same height. In other embodiments, only one distance measuring radar can be installed, and then the shape of the material surface on the other side (on the same diameter) is obtained through symmetrical processing.

The scanning measurement of the range radar to the charge level is a plurality of times of measurement in the radial direction of the furnace throat within the radius range of the furnace throat, and in the process, scanning is not required to be carried out with the step length of equal distance or equal angle, and only enough data need to be sampled and the distribution is not concentrated (random).

S2: and preprocessing radar ranging data of the measured charge level, and rejecting abnormal data in the radar ranging data.

In this embodiment the radar ranging data comprises the distance d between the level and the ranging radar₀And the angle alpha between the direction of the ranging radar and the horizontal plane.

The method of preprocessing in this embodiment includes: setting the range of radar ranging data according to the installation position of the ranging radar, thereby eliminating the radar ranging data which is not in the range; and according to the set echo intensity threshold value, eliminating radar ranging data with the echo intensity smaller than the echo intensity threshold value.

S3: converting the radar ranging data into coordinate data (x0, y) in a corresponding XY rectangular coordinate system₀)。

In this embodiment, the conversion formula is set as:

y₀＝d₀ sinα-h_r

wherein x is₀、y₀Are coordinate values of X-axis and Y-axis in the coordinate data, r represents furnace throat radius, h_rIndicates the height h from the installation position of the range radar to the zero stockline_lThe height of the upper boundary of the furnace throat from a zero material line is shown, and gamma represents the included angle between the upper furnace wall of the furnace throat and the horizontal plane.

In the XY rectangular coordinate system, the X axis is parallel to the scanning direction of the range radar (radial direction of the throat), and the Y axis is vertical.

S4: and fitting the shape of the charge level by using coordinate data through a Krigin interpolation method.

Step S4 specifically includes the following steps:

s41: carrying out radial segmentation on the charge level, setting an X-axis coordinate range corresponding to each segment, and setting each segmentation point as a point to be inserted;

s42: calculating any two coordinate data x in the coordinate data_iAnd x_jA distance d between_ij＝x_i-x_jGrouping all the calculated distances according to the distance; where i is 0,1, …, n, j is 0,1, …, n, n denotes the total number of coordinate data, i and j denote the serial numbers of the coordinate data, x_iAnd x_jRespectively representing the coordinates of the ith and the j th coordinate data in the X-axis direction;

s43: the average d of all distances in each group is calculated and calculated according toThe estimated value gamma of each group of corresponding half-variance functions^*(d)：

Wherein, y (x)_i) And y (x)_i+ d) respectively representing X as the coordinate in the X direction_iAnd (x)_iThe Y-axis direction coordinate corresponding to + d);

s44: the half-variance function model is constructed as follows: γ (d) ═ c (1-e)^-d/r) Wherein c and r are both parameters in the half variance function model, gamma (d) represents the half variance operation, and the estimated value gamma of each group of corresponding half variance functions is determined according to the half variance function model^*(d) Fitting to obtain a fitted half variance function model;

in the embodiment, the fitting method adopts a Levenberg-Marquardt method;

s45: the following linear system of equations was constructed:

according to the fitted half variance function model and the formula r_ij＝γ(d_ij) Calculating any two coordinate data x_iAnd x_jHalf-variance of an attribute r_ijAnd further obtaining a coefficient matrix in the linear equation system:

s46: for each point x to be inserted₀Calculating x₀With arbitrary coordinate data x_iHalf-variance of an attribute r_i0Further obtain the point x to be inserted₀Right-hand array in the corresponding linear system of equations:

s47: solving the linear equation set according to the coefficient matrix and the right-side array in the linear equation set to obtain the point x to be inserted₀Corresponding weight array:

s48: for each point x to be inserted₀According to the point x to be inserted₀Corresponding weight array, calculating the point x to be interpolated by₀Is estimated value of

Wherein, y_iRepresenting x in coordinate data_iThe corresponding Y-axis direction coordinates.

The method is also suitable for multi-radar sensing three-dimensional charge level shape measurement, and when multiple radars are used for measuring the charge level from multiple angles, after data are collected, a data dimension and Euclidean distance measurement mode in a Krigin method are set, and three-dimensional interpolation fitting can be carried out on the charge level.

Compared with the prior art, the embodiment of the invention does not need complex sensor arrangement, is beneficial to installation and maintenance, and has better measuring effect compared with the common single-point radar measurement.

Example two:

the invention also provides blast furnace burden surface shape measuring terminal equipment which comprises a memory, a processor and a computer program stored in the memory and capable of running on the processor, wherein the processor executes the computer program to realize the steps of the method embodiment of the first embodiment of the invention.

Further, as an executable scheme, the blast furnace burden surface shape measurement terminal device may be a desktop computer, a notebook computer, a palm computer, a cloud server, or other computing devices. The blast furnace burden surface shape measuring terminal equipment can comprise, but is not limited to, a processor and a memory. It will be understood by those skilled in the art that the above-mentioned blast furnace burden surface shape measurement terminal device is only an example of the blast furnace burden surface shape measurement terminal device, and does not constitute a limitation on the blast furnace burden surface shape measurement terminal device, and may include more or less components than the above, or may combine some components, or may be different components, for example, the blast furnace burden surface shape measurement terminal device may further include an input/output device, a network access device, a bus, etc., which is not limited by the embodiment of the present invention.

Further, as an executable solution, the Processor may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other Programmable logic device, a discrete Gate or transistor logic device, a discrete hardware component, and the like. The general purpose processor may be a microprocessor or the processor may be any conventional processor or the like, the processor being the control center of the blast furnace burden surface shape measurement terminal equipment, various interfaces and lines connecting the various parts of the entire blast furnace burden surface shape measurement terminal equipment.

The memory can be used for storing the computer program and/or the module, and the processor can realize various functions of the blast furnace burden surface shape measurement terminal equipment by operating or executing the computer program and/or the module stored in the memory and calling the data stored in the memory. The memory can mainly comprise a program storage area and a data storage area, wherein the program storage area can store an operating system and an application program required by at least one function; the storage data area may store data created according to the use of the mobile phone, and the like. In addition, the memory may include high speed random access memory, and may also include non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), at least one magnetic disk storage device, a Flash memory device, or other volatile solid state storage device.

The invention also provides a computer-readable storage medium, in which a computer program is stored, which, when being executed by a processor, carries out the steps of the above-mentioned method of an embodiment of the invention.

The integrated module/unit of the blast furnace burden surface shape measuring terminal device can be stored in a computer readable storage medium if it is implemented in the form of a software functional unit and sold or used as a separate product. Based on such understanding, all or part of the flow of the method according to the embodiments of the present invention may also be implemented by a computer program, which may be stored in a computer-readable storage medium, and when the computer program is executed by a processor, the steps of the method embodiments may be implemented. Wherein the computer program comprises computer program code, which may be in the form of source code, object code, an executable file or some intermediate form, etc. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, usb disk, removable hard disk, magnetic disk, optical disk, computer Memory, Read-Only Memory (ROM), Random Access Memory (RAM), software distribution medium, and the like.

While the invention has been particularly shown and described with reference to a preferred embodiment, it will be understood by those skilled in the art that various changes in form and detail may be made therein without departing from the spirit and scope of the invention as defined by the appended claims.