Method for determining operation level of hydraulic device on basis of kriging interpolation algorithm, and device

Abstract

Provided in the present disclosure are a method and apparatus for determining the operation level of a hydraulic device on the basis of a Kriging interpolation algorithm, and a device and a medium. The method comprises: acquiring three-dimensional data corresponding to operation parameter information of a hydraulic device, and a plurality of target operation parameter values; mapping the three-dimensional data to a two-dimensional grid corresponding to a unit rotational speed of the hydraulic device and a unit flow of the hydraulic device, so as to obtain an initial operation parameter grid; performing interpolation processing on the initial operation parameter grid, so as to obtain interpolation points, calculating interpolation-point operation parameter values on the basis of a Kriging interpolation algorithm, and repeating the interpolation process until the interpolation-point operation parameter values meet a preset parameter condition, so as to obtain an operation parameter grid; on the basis of the interpolation points and each target operation parameter value, determining a plurality of equal-value points, and performing curve fitting to obtain a contour line corresponding to each target operation parameter value; and on the basis of the contour line corresponding to each target operation parameter value, determining the operation level of the hydraulic device. The present disclosure realizes the accurate statistics and analysis of operation parameters of a hydraulic device, and thus realizes the accurate determination of the operation situation of the hydraulic device.

Description

Method and device for determining operating level of hydraulic equipment based on Kriging interpolation algorithm Technical Field

[0001] The present disclosure relates to the field of hydraulic equipment operation evaluation, and in particular to a method, device, equipment and medium for determining the operation level of hydraulic equipment. Background Art

[0002] The operating parameters of hydraulic equipment are constantly changing during operation. By conducting statistical analysis on the operating parameters, the operating status of the hydraulic equipment can be monitored.

[0003] Currently, statistics and analysis of the operating parameters of hydraulic equipment usually require plotting the test data of each operating point in a plane coordinate diagram. However, the two-dimensional plotting method has low accuracy, which leads to inaccurate judgment of the operating status of the hydraulic equipment.

[0004] In view of this, how to accurately count and analyze the operating parameters of hydraulic equipment and then accurately judge the operating conditions of hydraulic equipment has become an important technical issue. Summary of the Invention

[0005] In view of this, the purpose of the present disclosure is to propose a method, device, equipment and medium for determining the operating level of hydraulic equipment, so as to solve or partially solve the above problems.

[0006] Based on the above objectives, a first aspect of the present disclosure provides a method for determining an operating level of a hydraulic device, the method comprising:

[0007] Acquire three-dimensional data corresponding to operating parameter information of the hydraulic equipment and multiple target operating parameter values, wherein the three dimensions of the three-dimensional data are a unit speed of the hydraulic equipment, a unit flow rate of the hydraulic equipment, and an operating parameter value of the hydraulic equipment;

[0008] Mapping the three-dimensional data to a two-dimensional grid corresponding to a unit speed of the hydraulic equipment and a unit flow rate of the hydraulic equipment to obtain at least one initial operating parameter grid, wherein each grid vertex in the initial operating parameter grid represents the three-dimensional data;

[0009] For each initial operating parameter grid, interpolation processing is performed on the initial operating parameter grid to obtain an interpolation point, and the operating parameter value of the interpolation point corresponding to the interpolation point is calculated based on the Kriging interpolation algorithm. The interpolation process is repeated until the operating parameter value of the interpolation point meets the preset parameter condition, thereby obtaining an operating parameter grid;

[0010] For each target operating parameter value, multiple isovalue points are determined based on the interpolation points on all operating parameter grids and the target operating parameter value, and curve fitting is performed on the multiple isovalue points to obtain isovalue lines corresponding to the target operating parameter value;

[0011] The operation level of the hydraulic equipment is determined according to the contour lines corresponding to each target operation parameter value.

[0012] Based on the same inventive concept, the second aspect of the present disclosure provides a device for determining an operating level of a hydraulic device, comprising:

[0013] a parameter information acquisition module configured to acquire three-dimensional data corresponding to operating parameter information of the hydraulic equipment and a plurality of target operating parameter values, wherein the three dimensions of the three-dimensional data are a unit speed of the hydraulic equipment, a unit flow rate of the hydraulic equipment, and an operating parameter value of the hydraulic equipment;

[0014] a parameter grid determination module configured to map the three-dimensional data to a two-dimensional grid corresponding to a unit speed of the hydraulic equipment and a unit flow rate of the hydraulic equipment to obtain at least one initial operating parameter grid, wherein each grid vertex in the initial operating parameter grid represents the three-dimensional data;

[0015] an interpolation processing module configured to perform interpolation processing on each initial operating parameter grid to obtain an interpolation point, calculate an interpolation point operating parameter value corresponding to the interpolation point based on a Kriging interpolation algorithm, and repeat the interpolation process until the interpolation point operating parameter value meets a preset parameter condition, thereby obtaining an operating parameter grid;

[0016] an isovalue point determination module configured to determine, for each target operating parameter value, a plurality of isovalue points based on interpolation points on all operating parameter grids and the target operating parameter value, and perform curve fitting on the plurality of isovalue points to obtain an isovalue line corresponding to the target operating parameter value;

[0017] The operation level determination module is configured to determine the operation level of the hydraulic equipment according to the contour line corresponding to each target operation parameter value.

[0018] Based on the same inventive concept, the third aspect of the present disclosure proposes an electronic device, comprising a memory, a processor, and a computer program stored on the memory and executable by the processor, wherein the processor implements the above-mentioned method for determining the operating level of the hydraulic equipment when executing the computer program.

[0019] Based on the same inventive concept, a fourth aspect of the present disclosure proposes a non-transitory computer-readable storage medium, which stores computer instructions for causing a computer to execute the method for determining the operating level of hydraulic equipment as described above.

[0020] As can be seen from the above, the present disclosure proposes a method, apparatus, equipment, and medium for determining the operating level of hydraulic equipment, obtaining three-dimensional data corresponding to the operating parameter information of the hydraulic equipment, wherein the three dimensions of the three-dimensional data are the unit speed of the hydraulic equipment, the unit flow rate of the hydraulic equipment, and the operating parameter value of the hydraulic equipment, mapping the obtained three-dimensional data to a two-dimensional grid corresponding to the unit speed of the hydraulic equipment and the unit flow rate of the hydraulic equipment, and obtaining at least one initial operating parameter grid for subsequent interpolation processing of each initial operating parameter grid. For each initial operating parameter grid, the initial operating parameter grid is interpolated to obtain an interpolation point, and the interpolation point operating parameter value corresponding to the interpolation point is calculated based on the Kriging interpolation algorithm. The interpolation process is repeated until the interpolation point operating parameter value meets the preset parameter conditions, thereby obtaining an operating parameter grid. By comparing the interpolation point operating parameter value with the preset parameter conditions, better control of the interpolation accuracy and effectiveness is achieved, reducing invalid or inefficient investment. Multiple target operating parameter values ​​are obtained. For each target operating parameter value, multiple isovalue points are determined based on the interpolation points on the operating parameter grid and the target operating parameter value. Curve fitting is performed on these multiple isovalue points to obtain isovalue lines corresponding to the target operating parameter value. This provides more accurate identification of isovalue points, and thus more accurate fitting of isovalue lines. The operating level of the hydraulic equipment is determined based on the isovalue lines corresponding to each target operating parameter value, enabling accurate judgment of the hydraulic equipment's operating status, facilitating subsequent determination of whether the hydraulic equipment needs optimization. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] In order to more clearly illustrate the technical solutions in the present disclosure or related technologies, the following briefly introduces the drawings required for use in the embodiments or related technical descriptions. Obviously, the drawings described below are only embodiments of the present disclosure. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0022] FIG1 is a flow chart of a method for determining an operating level of hydraulic equipment based on a Kriging interpolation algorithm according to an embodiment of the present disclosure;

[0023] FIG2 is a schematic diagram of three-dimensional data corresponding to operating parameter information of hydraulic equipment based on a Kriging interpolation algorithm according to an embodiment of the present disclosure;

[0024] FIG3 is a schematic diagram of label information according to an embodiment of the present disclosure;

[0025] FIG4 is a schematic diagram of interpolation points according to an embodiment of the present disclosure;

[0026] FIG5 is a structural block diagram of a device for determining an operating level of hydraulic equipment based on a Kriging interpolation algorithm according to an embodiment of the present disclosure;

[0027] FIG6 is a schematic structural diagram of an electronic device according to an embodiment of the present disclosure. DETAILED DESCRIPTION

[0028] In order to make the objectives, technical solutions and advantages of the present disclosure more clearly understood, the present disclosure is further described in detail below in conjunction with specific embodiments and with reference to the accompanying drawings.

[0029] It should be noted that, unless otherwise defined, the technical terms or scientific terms used in the embodiments of the present disclosure should have the usual meanings understood by people with ordinary skills in the field to which the present disclosure belongs. The "first", "second" and similar words used in the embodiments of the present disclosure do not indicate any order, quantity or importance, but are only used to distinguish different components. "Include" or "comprise" and similar words mean that the elements or objects appearing before the word include the elements or objects listed after the word and their equivalents, without excluding other elements or objects. "Connect" or "connected" and similar words are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect. "Up", "down", "left", "right" and the like are only used to indicate relative position relationships. When the absolute position of the described object changes, the relative position relationship may also change accordingly.

[0030] Based on the above description, this embodiment proposes a method for determining the operating level of hydraulic equipment based on the Kriging interpolation algorithm, as shown in FIG1 , the method comprising:

[0031] Step 101: Acquire three-dimensional data and multiple target operating parameter values ​​corresponding to the operating parameter information of the hydraulic equipment, wherein the three dimensions of the three-dimensional data are the unit speed of the hydraulic equipment, the unit flow of the hydraulic equipment, and the operating parameter value of the hydraulic equipment.

[0032] In a specific implementation, operating parameter information of hydraulic equipment is obtained. The hydraulic equipment includes at least one of the following: a water pump, a water turbine, or a pump-turbine. The water pump includes at least one of the following types: a piston pump, a plunger pump, or a rotor pump. The water turbine includes at least one of the following types: a tangential-impact turbine, a skew-impact turbine, a double-impact turbine, a Francis turbine, a diagonal-flow turbine, or an axial-flow turbine. The operating parameter information includes at least one of the following: energy conversion efficiency of the hydraulic equipment, pressure pulsation of the hydraulic equipment, or incipient cavitation value of the hydraulic equipment.

[0033] The unit speed and unit flow of the hydraulic equipment are obtained, and three-dimensional data corresponding to the operating parameter information is determined based on the unit speed, unit flow and operating parameter information of the hydraulic equipment, as shown in Figure 2. Multiple preset target operating parameter values ​​are obtained.

[0034] Exemplarily, the operating parameter information is the energy conversion efficiency of the hydraulic equipment, and the three dimensions of the three-dimensional data are the unit speed of the hydraulic equipment, the unit flow rate of the hydraulic equipment, and the energy conversion efficiency of the hydraulic equipment.

[0035] Step 102: Map the three-dimensional data to a two-dimensional grid corresponding to a unit speed and a unit flow of the hydraulic equipment to obtain at least one initial operating parameter grid, wherein each grid vertex in the initial operating parameter grid is the three-dimensional data.

[0036] In a specific implementation, the obtained three-dimensional data is mapped to a two-dimensional grid corresponding to the unit speed and unit flow rate of the hydraulic equipment. Each grid vertex in the initial operating parameter grid represents the three-dimensional data. In other words, the three-dimensional data corresponding to the operating parameter information of every four hydraulic equipment constitutes an initial operating parameter grid.

[0037] Exemplarily, the three-dimensional data corresponding to the operating parameter information of four hydraulic equipment constitutes an initial operating parameter grid, the three-dimensional data corresponding to the operating parameter information of six hydraulic equipment constitutes two initial operating parameter grids, and the three-dimensional data corresponding to the operating parameter information of nine hydraulic equipment constitutes three output operating parameter grids.

[0038] After the three-dimensional data is reflected in the two-dimensional grid, the value corresponding to the three-dimensional data is the operating parameter value of the operating parameter information.

[0039] Step 103: For each initial operating parameter grid, interpolation processing is performed on the initial operating parameter grid to obtain interpolation points, and the interpolation point operating parameter values ​​corresponding to the interpolation points are calculated based on the Kriging interpolation algorithm. The interpolation process is repeated until the interpolation point operating parameter values ​​meet the preset parameter conditions, thereby obtaining an operating parameter grid.

[0040] In a specific implementation, for each initial operating parameter grid, interpolation processing is performed on the initial operating parameter grid to obtain interpolation points. The interpolation point operating parameter value corresponding to each interpolation point is calculated using the Kriging interpolation algorithm.

[0041] The initial operating parameter grid is repeatedly interpolated until the obtained operating parameter value at the interpolation point meets the preset parameter condition. The interpolation process ends and the initial operating parameter grid after interpolation processing, that is, the operating parameter grid, is obtained.

[0042] Step 104 , for each target operating parameter value, determine multiple isovalue points based on the interpolation points on all operating parameter grids and the target operating parameter value, perform curve fitting on the multiple isovalue points, and obtain isovalue lines corresponding to the target operating parameter value.

[0043] In a specific implementation, after interpolation processing is performed on each initial operating parameter grid, at least one operating parameter grid is obtained, and all interpolation points in all operating parameter grids are obtained. For each preset target operating parameter value, multiple equivalent points are determined based on all interpolation points, i.e., the target operating parameter value. Each equivalent point is equal to the target operating parameter value. Curve fitting is performed on these multiple equivalent points to obtain the corresponding contour lines of the target operating parameter value.

[0044] Step 105 : determining the operating level of the hydraulic equipment according to the contour lines corresponding to each target operating parameter value.

[0045] In specific implementation, the above steps are used to respectively determine the contour lines corresponding to each target operating parameter value, and the operating level of the hydraulic equipment is determined based on the multiple contour lines.

[0046] Contour lines highlight areas with similar data values ​​and allow users to easily view the distribution and trends of hydraulic equipment operating parameters. Multiple contour lines can form a contour map. By analyzing the distribution of operating parameter values ​​along each contour line, the operating level of the hydraulic equipment can be determined.

[0047] Obtain the minimum value of the preset operating parameter value, obtain the number of data points on the contour line whose operating parameter value is less than the minimum value of the operating parameter value, record it as the first number, obtain the number of data points on all contour lines, record it as the second number. Calculate the ratio of the first number to the second number, determine that if the ratio is less than the first preset threshold, the operating level of the hydraulic equipment is the first level; determine that if the ratio is greater than or equal to the first preset threshold and less than the second preset threshold, the operating level of the hydraulic equipment is the second level; determine that if the ratio is greater than or equal to the second preset threshold and less than the third preset threshold, the operating level of the hydraulic equipment is the third level; determine that if the ratio is greater than or equal to the third preset threshold, the operating level of the hydraulic equipment is the fourth level. The above method of dividing the operating levels of hydraulic equipment is only an example, and the specific number of divided levels and the size of the preset threshold are not limited.

[0048] Exemplarily, the operating parameter information is the energy conversion efficiency of the hydraulic equipment, and the minimum value of the preset operating parameter value is 90.5%. There are four contour lines, corresponding to energy conversion efficiencies of 90.5%, 90%, 89%, and 88%, respectively. It is determined that the number of data points on the contour lines where the operating parameter value is less than the minimum value of the operating parameter value is 8, that is, the number of data points on the contour lines 89% and 88% is 8, and the number of data points on all contour points is 20, so the ratio is 0.4. If the first preset threshold is 0.2, the second preset threshold is 0.35, and the third preset threshold is 0.5, then the current operating level of the hydraulic equipment is the third level.

[0049] Through the above scheme, three-dimensional data corresponding to the operating parameter information of the hydraulic equipment is obtained, where the three dimensions of the three-dimensional data are the unit speed of the hydraulic equipment, the unit flow rate of the hydraulic equipment, and the operating parameter value of the hydraulic equipment. The obtained three-dimensional data is mapped to a two-dimensional grid corresponding to the unit speed of the hydraulic equipment and the unit flow rate of the hydraulic equipment, thereby obtaining at least one initial operating parameter grid, so as to facilitate subsequent interpolation processing of each initial operating parameter grid. For each initial operating parameter grid, the initial operating parameter grid is interpolated to obtain an interpolation point. The operating parameter value of the interpolation point corresponding to the interpolation point is calculated based on the Kriging interpolation algorithm. The interpolation process is repeated until the operating parameter value of the interpolation point meets the preset parameter conditions, thereby obtaining the operating parameter grid. By comparing the operating parameter value of the interpolation point with the preset parameter conditions, better control of the accuracy and effectiveness of the interpolation is achieved, thereby reducing invalid or inefficient investment. Multiple target operating parameter values ​​are obtained. For each target operating parameter value, multiple isovalue points are determined based on the interpolation points on the operating parameter grid and the target operating parameter value. Curve fitting is performed on these multiple isovalue points to obtain isovalue lines corresponding to the target operating parameter value. This provides more accurate identification of isovalue points, and thus more accurate fitting of isovalue lines. The operating level of the hydraulic equipment is determined based on the isovalue lines corresponding to each target operating parameter value, enabling accurate judgment of the hydraulic equipment's operating status, facilitating subsequent determination of whether the hydraulic equipment needs optimization.

[0050] In some embodiments, for each initial operating parameter grid, step 103 specifically includes:

[0051] Step 1031: Obtain the interpolation point operating parameter value corresponding to each interpolation point in the initial operating parameter grid.

[0052] Step 1032 : Calculate the difference between the interpolation point operating parameter values ​​of any two interpolation points. In response to the difference being less than a preset difference threshold, determine whether the interpolation point operating parameter value meets a preset parameter condition.

[0053] In a specific implementation, for each initial operating parameter grid, all interpolation points in the initial operating parameter grid and the interpolation point operating parameter value corresponding to each interpolation point are obtained. The difference in the operating parameter values ​​between any two interpolation points is calculated. When it is determined that the difference in the operating parameter between any two interpolation points is less than a preset difference threshold, it is determined that the interpolation point operating parameter value meets the preset parameter condition.

[0054] For example, an initial operating parameter grid A has four interpolation points, B, C, D, and E, with a preset difference threshold of 0.2. The differences in the operating parameter values ​​at the interpolation points are calculated for each of the four interpolation points, yielding a total of six differences. If all six differences are less than 0.2, the operating parameter values ​​at the interpolation points meet the preset parameter conditions, meaning no further interpolation processing is required, and the operating parameter grid is obtained.

[0055] or,

[0056] Step 1033 : In response to the interpolation point operating parameter value corresponding to each interpolation point being smaller than the target operating parameter value, it is determined that the interpolation point operating parameter value meets a preset parameter condition.

[0057] In a specific implementation, for each initial operating parameter grid, all interpolation points in the initial operating parameter grid and the interpolation point operating parameter value corresponding to each interpolation point are obtained. If the interpolation point operating parameter value corresponding to each interpolation point is determined to be smaller than the target operating parameter value, then the interpolation point operating parameter value is determined to meet the preset parameter condition.

[0058] For example, the initial operating parameter grid F has four interpolation points G, H, I, and J, and the target operating parameter value is 75. The interpolation point operating parameter value corresponding to interpolation point G is 60, the interpolation point operating parameter value corresponding to interpolation point H is 65, the interpolation point operating parameter value corresponding to interpolation point I is 50, and the interpolation point operating parameter value corresponding to interpolation point J is 68. These interpolation point operating parameter values ​​meet the preset parameter conditions, meaning that no interpolation processing is required to obtain the operating parameter grid.

[0059] or,

[0060] Step 1034 : In response to the interpolation point operating parameter value corresponding to each interpolation point being greater than the target operating parameter value, determining whether the interpolation point operating parameter value meets a preset parameter condition.

[0061] In a specific implementation, for each initial operating parameter grid, all interpolation points in the initial operating parameter grid and the interpolation point operating parameter value corresponding to each interpolation point are obtained. If it is determined that the interpolation point operating parameter value corresponding to each interpolation point is greater than the target operating parameter value, then it is determined that the interpolation point operating parameter value meets the preset parameter condition.

[0062] For example, the initial operating parameter grid J has four interpolation points K, L, M, and N, and the target operating parameter value is 80. The interpolation point operating parameter value corresponding to interpolation point K is 88, the interpolation point operating parameter value corresponding to interpolation point L is 95, the interpolation point operating parameter value corresponding to interpolation point M is 92, and the interpolation point operating parameter value corresponding to interpolation point N is 90. These interpolation point operating parameter values ​​meet the preset parameter conditions, meaning that no interpolation processing is required to obtain the operating parameter grid.

[0063] Through the above scheme, by determining the preset parameter conditions corresponding to the operating parameter values ​​of the interpolation points, better control of the interpolation accuracy is achieved. When the preset parameter conditions are met, it indicates that the interpolation accuracy has met the requirements, and the interpolation process can be stopped, thereby ensuring the effectiveness of the interpolation.

[0064] In some embodiments, step 104 specifically includes:

[0065] Step 1041 : for each operating parameter grid: obtain all interpolation points included in the operating parameter grid, and select any two interpolation points from the all interpolation points as a first target interpolation point and a second target interpolation point.

[0066] Step 1042: Determine the magnitude relationship between the first interpolation point operating parameter value corresponding to the first target interpolation point and the second interpolation point operating parameter value corresponding to the second target interpolation point, and the target operating parameter value.

[0067] In a specific implementation, for each operating parameter grid, all interpolation points in the operating parameter grid are obtained, and any two interpolation points are selected as the first target interpolation point and the second target interpolation point. The operating parameter value of the first interpolation point corresponding to the first target interpolation point, that is, the operating parameter value of the second interpolation point corresponding to the second target interpolation point, is obtained. The magnitude relationship between the operating parameter value of the first interpolation point and the target operating parameter value, and the magnitude relationship between the operating parameter value of the second interpolation point and the target operating parameter value are compared respectively.

[0068] Step 1043: In response to the fact that, among the first interpolation point operating parameter value and the second interpolation point operating parameter value, there exists an interpolation point operating parameter value greater than the target operating parameter value and an interpolation point operating parameter value less than the target operating parameter value, it is determined that there is an equal value point between the first target interpolation point and the second target interpolation point.

[0069] In a specific implementation, when there is an interpolation point operating parameter value greater than the target operating parameter value and an interpolation point operating parameter value less than the target operating parameter value among the first interpolation point operating parameter value and the second interpolation point operating parameter value, it is determined that there is an equal value point between the first target interpolation point and the second target interpolation point.

[0070] That is, if the operating parameter value of the first interpolation point is greater than the target operating parameter value, and the operating parameter value of the second interpolation point is less than the target operating parameter value, it is determined that there is an equal value point between the first target interpolation point and the second target interpolation point.

[0071] or,

[0072] If the operating parameter value of the first interpolation point is less than the target operating parameter value, and the operating parameter value of the second interpolation point is greater than the target operating parameter value, it is determined that an equal value point exists between the first target interpolation point and the second target interpolation point.

[0073] For example, the target operating parameter value is 65, the acquired operating parameter value of the first interpolation point is 40, and the acquired operating parameter value of the second interpolation point is 89, then there is an equal value point between the first target interpolation point and the second target interpolation point.

[0074] Step 1044 : Calculate a first proportional relationship according to the first interpolation point operating parameter value, the second interpolation point operating parameter value, and the target operating parameter value.

[0075] During specific implementation, the proportional relationship between the acquired operating parameter value at the first interpolation point and the acquired operating parameter value at the second interpolation point and the target operating parameter value is calculated to obtain the first proportional relationship.

[0076] Step 1045: Obtain the first unit speed and the first unit flow corresponding to the first target interpolation point, and the second unit speed and the second unit flow corresponding to the second target interpolation point. According to the first unit speed, the first unit flow, the second unit speed and the second unit flow, the target unit speed and the target unit flow corresponding to the equivalent point are calculated based on the first proportional relationship.

[0077] Step 1046: Obtain an equivalent point according to the target unit speed, the target unit flow rate, and the target operating parameter value.

[0078] In specific implementation, since the equivalent points are all within the initial operating parameter grid, the first unit speed and first unit flow corresponding to the first target interpolation point and the second unit speed and second unit flow corresponding to the second target interpolation point can be calculated.

[0079] For example, if the unit rotational speeds corresponding to the three-dimensional data on the left and right sides of the first target interpolation point are 80 and 90, respectively, then the first unit rotational speed corresponding to the first target interpolation point is 85. If the unit flow rates corresponding to the three-dimensional data on the left and right sides of the first target interpolation point are 40 and 50, respectively, then the first unit flow rate corresponding to the first target interpolation point is 45.

[0080] Based on the calculated first proportional relationship and the acquired first unit speed, first unit flow, second unit speed and second unit flow, the target unit speed and target unit flow corresponding to the equivalent point are calculated.

[0081] The equivalent points and interpolation points are all three-dimensional data, so the corresponding equivalent points can be determined based on the target unit speed and target unit flow rate obtained by calculation, as well as the target operating parameter values.

[0082] In some embodiments, after step 103, the method further includes:

[0083] Step 10A: obtaining the origin of the two-dimensional grid corresponding to the unit speed and unit flow of the hydraulic equipment, and labeling the operating parameter grid according to the horizontal axis and vertical axis of the two-dimensional grid to obtain the label information of the operating parameter grid.

[0084] During specific implementation, the origin of the two-dimensional grid corresponding to the unit speed and unit flow of the hydraulic equipment is obtained, and the operating parameter grid is labeled according to the horizontal axis and vertical axis of the two-dimensional grid to obtain the label information of each operating parameter grid, which is two-dimensional information.

[0085] Exemplarily, as shown in Figure 3, the origin of the two-dimensional grid is O. According to the horizontal and vertical axes of the two-dimensional grid, the label information corresponding to the operating parameter grid A is (x1, y1), the label information corresponding to the operating parameter grid B is (x2, y1), the label information corresponding to the operating parameter grid C is (x2, y1), and the label information corresponding to the operating parameter grid D is (x2, y2).

[0086] In some embodiments, step 104 specifically includes:

[0087] Step 1041: For each operating parameter grid, obtain the equivalent points contained in the operating parameter grid.

[0088] Step 1042 : Connect all the isovalue points contained in each operating parameter grid in sequence along a preset direction according to the label information of the operating parameter grid to obtain an isovalue line corresponding to the target operating parameter value.

[0089] In a specific implementation, when performing curve fitting on the isovalue points, all isovalue points contained in the operating parameter grid are first obtained. According to the preset direction, the isovalue points contained in each parameter operating grid are sequentially connected according to the label information of the operating parameter grid to obtain the isovalue line corresponding to the target operating parameter value.

[0090] For example, taking the label information in step 10A as an example, when the preset direction is horizontal, the isovalue points included in operating parameter grid A are sequentially connected horizontally to the isovalue points included in operating parameter grid B. After all the isovalue points in operating parameter grids A and B are connected, the isovalue points included in operating parameter grid C are obtained and sequentially connected horizontally to the isovalue points included in operating parameter grid D, ultimately obtaining the isovalue line corresponding to the target operating parameter value.

[0091] In another example, using the label information in step 10A as an example, when the preset direction is vertical, the isovalue points included in operating parameter grid A are sequentially connected vertically to the isovalue points included in operating parameter grid C. After all the isovalue points in operating parameter grids A and C are connected, the isovalue points included in operating parameter grid B are obtained and sequentially connected vertically to the isovalue points included in operating parameter grid D, ultimately obtaining the isovalue line corresponding to the target operating parameter value.

[0092] The above scheme, by labeling the operating parameter grid, allows for accurate determination of the corresponding fitting order for equivalent point fitting, avoiding omissions of equivalent point fitting, and also avoiding duplicate fitting, overfitting, and invalid fitting. Furthermore, the accurate determination of the corresponding contour lines for the operating parameter information of the hydraulic equipment has important guiding significance for the selection of power station operating areas and the operating efficiency of power stations. It can also provide a reference for evaluating the level of hydraulic development of turbines and a method for runner hydraulic optimization, which is of great significance to the technological development of the hydropower industry.

[0093] In some embodiments, for each initial operating parameter grid, step 103 specifically includes:

[0094] Step 103A: Obtain the corresponding four grid vertices in the initial operating parameter grid.

[0095] Step 103B: Select the middle point of the line connecting any two adjacent mesh vertices as the first initial interpolation point.

[0096] Step 103C: Select the intersection point of the lines connecting two diagonal mesh vertices and use the intersection point as the second initial interpolation point.

[0097] Step 103D: Count the first initial interpolation points and the second initial interpolation points to obtain interpolation points corresponding to the initial operating parameter grid.

[0098] During specific implementation, the four corresponding grid vertices in the initial operating parameter grid are obtained, the middle point of the line connecting any two adjacent grid vertices is selected as the first initial interpolation point, the intersection point of the line connecting two diagonal grid vertices is selected, and the intersection point is used as the second initial interpolation point. The first initial interpolation point and the second initial interpolation point are counted to obtain the interpolation point corresponding to the initial operating parameter grid.

[0099] Exemplarily, as shown in Figure 4, the corresponding four grid vertices in the initial operating parameter grid are A, B, C, and D, and the midpoint between A and B is E, the midpoint between B and C is F, the midpoint between C and D is G, the midpoint between D and A is H, and the intersection point O of the line between A and C and the line between B and D are obtained respectively. The midpoints E, F, G, and H are the first initial interpolation points, and the intersection point O is the second initial interpolation point. The midpoints E, F, G, H and the intersection point O are the interpolation points corresponding to the initial operating parameter grid.

[0100] In some embodiments, for each initial operating parameter grid, step 103 specifically includes:

[0101] Step 103a: Calculate the first plane distance and the first semi-variance value between any two mesh vertices in each initial operating parameter mesh obtained according to the three-dimensional data to obtain a first plane distance set and a first semi-variance value set.

[0102] In a specific implementation, the first plane distance and the first semi-variance value between any two grid vertices in each initial operating parameter grid obtained according to the three-dimensional data are calculated to obtain a first plane distance set and a first semi-variance value set.

[0103] The first plane distance is expressed by the formula:

[0104] Among them, the coordinates of the grid vertex A are (N i , Q i ), the coordinates of mesh vertex B are (N i , Q i ), the grid vertex A and the grid vertex B are both grid vertices in the initial operation parameter grid, the d ij is the first plane distance between mesh vertex A and mesh vertex B.

[0105] The first semivariance value is expressed by the formula:

[0106] Among them, γ ij is the first semivariance between mesh vertex A and mesh vertex B, η iis the operating parameter value in the three-dimensional data corresponding to the grid vertex A, η j is the operating parameter value in the three-dimensional data corresponding to the mesh vertex B, and E is the expectation.

[0107] Step 103b: Fit the first plane distance set and the first semi-variance value set using the least square method to obtain a first functional relationship.

[0108] Step 103c: Calculate the second plane distance between the interpolation point and each grid vertex in each initial operating parameter grid to obtain a second plane distance set.

[0109] Step 103d: Based on the first functional relationship and the second plane distance set, obtain a second semi-variance value set between the interpolation point and each grid vertex in each initial operating parameter grid.

[0110] Specifically, a first plane distance set and a first semi-variance value set are fitted using the least squares method to obtain a first functional relationship. A second plane distance set is calculated between the interpolation point and each mesh vertex in each initial operating parameter grid to obtain a second plane distance set. Based on the first functional relationship and the second plane distance set, a second semi-variance value set is obtained between the interpolation point and each mesh vertex in each initial operating parameter grid.

[0111] Step 103e: Calculate the interpolation point operating parameter value corresponding to the interpolation point based on the Kriging interpolation algorithm and the second semi-variance value set.

[0112] In specific implementation, according to the definition of the Kriging algorithm, the value of the unknown point is estimated by weighted summation of the data of all known points in space. Therefore, the interpolation point operating parameter value corresponding to the interpolation point is calculated based on the Kriging interpolation algorithm and the second semi-variance value set. The interpolation point operating parameter value is expressed by the formula:

[0113] in, is the interpolation point (N m , Q m ) corresponds to the estimated value of the operating parameter, λ i is the weight coefficient of the i-th mesh vertex, η i is the running parameter value of the i-th mesh vertex, k is the number of mesh vertices. The weight coefficient is the one that can satisfy the interpolation point (N m , Q m ) The set of coefficients that minimizes the variance from the true value.

[0114] In some embodiments, step 103e specifically includes:

[0115] Step 103e1 , obtaining an operating parameter value corresponding to each mesh vertex in each initial operating parameter mesh obtained based on the three-dimensional data, to obtain an operating parameter set.

[0116] Step 103e2, calculating the expected value and variance value of the operating parameter set.

[0117] Step 103e3, obtaining a random deviation value of a preset interpolation point, and calculating the estimation error of the interpolation point based on the Kriging interpolation algorithm according to the random deviation value, the expected value and the variance value.

[0118] In a specific implementation, the operating parameter value corresponding to each grid vertex in each initial operating parameter grid obtained based on the three-dimensional data is obtained to obtain an operating parameter set. The expected value and variance value of the operating parameter set are calculated, and the random deviation value of the preset interpolation point is obtained. The estimated error of the interpolation point is calculated based on the random deviation value, the expected value and the variance value, and the Kriging interpolation algorithm. E[η(Nm, Qm)] = E[η] = c Var[η(Nm, Qm)] = σ 2 η(Nm, Qm)=c+R(Nm, Qm)

[0119] Among them, c is the expected value of the operating parameter set, σ 2 is the variance value of the operating parameter set, R(N, Q) is the random deviation value, δ is the estimation error, and η(Nm, Qm) is the operating parameter value of the interpolation point.

[0120] Step 103e4: Calculate the weight coefficient of each mesh vertex using the estimation error and the second semivariance value set.

[0121] Step 103e5: According to the weight coefficient of each mesh vertex and the operating parameter value corresponding to each mesh vertex, the interpolation point operating parameter value corresponding to the interpolation point is calculated based on the Kriging interpolation algorithm.

[0122] In specific implementation, the estimated error and the second semivariance value set are substituted into the pre-constructed weight coefficient equation group to obtain the weight coefficient of each mesh vertex, wherein the weight coefficient equation group is expressed by the formula:

[0123] According to the weight coefficient of each mesh vertex and the operating parameter value corresponding to each mesh vertex, the interpolation point operating parameter value corresponding to the interpolation point is calculated based on the Kriging interpolation algorithm, that is, substituted into the formula shown in step 103e to obtain the interpolation point operating parameter value corresponding to the interpolation point.

[0124] It should be noted that the method of the embodiments of the present disclosure can be performed by a single device, such as a computer or server. The method of the embodiments of the present disclosure can also be applied in a distributed scenario, where multiple devices cooperate to perform the method. In such a distributed scenario, one of the multiple devices may only perform one or more steps of the method of the embodiments of the present disclosure, and the multiple devices will interact with each other to complete the method.

[0125] It should be noted that the above description is limited to some embodiments of the present disclosure. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recited in the claims may be performed in an order different from that described in the above embodiments and still achieve the desired results. Furthermore, the processes depicted in the accompanying drawings do not necessarily require the specific order or sequential order shown to achieve the desired results. In certain embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0126] Based on the same inventive concept, corresponding to any of the above-mentioned embodiments and methods, the present disclosure also provides an operating level determination device for hydraulic equipment.

[0127] Referring to FIG5 , FIG5 is a diagram illustrating an operating level determination device for hydraulic equipment according to an embodiment, comprising:

[0128] The parameter information acquisition module 501 is configured to acquire three-dimensional data corresponding to the operating parameter information of the hydraulic equipment and a plurality of target operating parameter values, wherein the three dimensions of the three-dimensional data are the unit speed of the hydraulic equipment, the unit flow rate of the hydraulic equipment, and the operating parameter value of the hydraulic equipment;

[0129] a parameter grid determination module 502 configured to map the three-dimensional data to a two-dimensional grid corresponding to a unit speed of the hydraulic equipment and a unit flow rate of the hydraulic equipment to obtain at least one initial operating parameter grid, wherein each grid vertex in the initial operating parameter grid represents the three-dimensional data;

[0130] The interpolation processing module 503 is configured to perform interpolation processing on each initial operating parameter grid to obtain an interpolation point, calculate the operating parameter value of the interpolation point corresponding to the interpolation point based on the Kriging interpolation algorithm, and repeat the interpolation process until the operating parameter value of the interpolation point meets the preset parameter conditions, thereby obtaining an operating parameter grid;

[0131] The isovalue point determination module 504 is configured to determine, for each target operating parameter value, a plurality of isovalue points based on the interpolation points on all operating parameter grids and the target operating parameter value, and perform curve fitting on the plurality of isovalue points to obtain an isovalue line corresponding to the target operating parameter value;

[0132] The operation level determination module 505 is configured to determine the operation level of the hydraulic equipment according to the contour line corresponding to each target operation parameter value.

[0133] In some embodiments, for each initial operating parameter grid, the interpolation processing module 503 is specifically configured to:

[0134] Obtaining an interpolation point operating parameter value corresponding to each interpolation point in the initial operating parameter grid;

[0135] Calculating a difference between interpolation point operating parameter values ​​of any two interpolation points, and in response to the difference being less than a preset difference threshold, determining that the interpolation point operating parameter value meets a preset parameter condition;

[0136] or,

[0137] In response to the interpolation point operating parameter value corresponding to each interpolation point being less than the target operating parameter value, determining that the interpolation point operating parameter value meets a preset parameter condition;

[0138] or,

[0139] In response to the interpolation point operating parameter value corresponding to each interpolation point being greater than the target operating parameter value, it is determined that the interpolation point operating parameter value meets the preset parameter condition.

[0140] In some embodiments, the equivalent point determination module 504 specifically includes:

[0141] The interpolation point acquisition unit is configured to: for each operating parameter grid: acquire all interpolation points included in the operating parameter grid, and select any two interpolation points from all the interpolation points as a first target interpolation point and a second target interpolation point;

[0142] an operating parameter value determining unit configured to determine a magnitude relationship between a first interpolation point operating parameter value corresponding to the first target interpolation point and a second interpolation point operating parameter value corresponding to the second target interpolation point, and the target operating parameter value;

[0143] an equivalue point determination unit configured to determine that an equivalue point exists between the first target interpolation point and the second target interpolation point in response to the presence of an interpolation point operating parameter value greater than the target operating parameter value and an interpolation point operating parameter value less than the target operating parameter value among the first interpolation point operating parameter value and the second interpolation point operating parameter value;

[0144] a proportional relationship determining unit, configured to calculate a first proportional relationship according to the first interpolation point operating parameter value, the second interpolation point operating parameter value, and the target operating parameter value;

[0145] a unit speed determining unit configured to obtain a first unit speed and a first unit flow rate corresponding to the first target interpolation point, and a second unit speed and a second unit flow rate corresponding to the second target interpolation point, and calculate, based on the first proportional relationship, a target unit speed and a target unit flow rate corresponding to the equivalent point according to the first unit speed, the first unit flow rate, the second unit speed, and the second unit flow rate;

[0146] The equivalent point determination unit is configured to obtain an equivalent point according to the target unit speed, the target unit flow rate and the target operating parameter value.

[0147] In some embodiments, the device also includes a labeling module, which is specifically configured to obtain the origin of the two-dimensional grid corresponding to the unit speed of the hydraulic equipment and the unit flow of the hydraulic equipment, and label the operating parameter grid according to the horizontal axis and vertical axis of the two-dimensional grid to obtain the labeling information of the operating parameter grid.

[0148] In some embodiments, the isovalue point determination module 504 is further configured to obtain, for each operating parameter grid, isovalue points contained in the operating parameter grid;

[0149] Along a preset direction, all the isovalue points contained in each operating parameter grid are connected in sequence according to the label information of the operating parameter grid to obtain the isovalue line corresponding to the target operating parameter value.

[0150] In some embodiments, for each initial operating parameter grid, the interpolation processing module 503 is specifically configured to:

[0151] Obtaining four corresponding grid vertices in the initial operating parameter grid;

[0152] Select the middle point of the line connecting any two adjacent mesh vertices as the first initial interpolation point;

[0153] Selecting an intersection point of a line connecting two diagonal mesh vertices, and using the intersection point as a second initial interpolation point;

[0154] The first initial interpolation points and the second initial interpolation points are counted to obtain interpolation points corresponding to the initial operating parameter grid.

[0155] In some embodiments, for each initial operating parameter grid, the interpolation processing module 503 specifically further includes:

[0156] A first plane distance calculation unit is configured to calculate a first plane distance and a first semi-variance value between any two mesh vertices in each initial operating parameter mesh obtained according to the three-dimensional data, to obtain a first plane distance set and a first semi-variance value set;

[0157] a functional relationship determining unit, configured to fit the first plane distance set and the first semivariance value set using a least squares method to obtain a first functional relationship;

[0158] A second plane distance determination unit is configured to calculate a second plane distance between the interpolation point and each grid vertex in each initial operating parameter grid to obtain a second plane distance set;

[0159] a second semivariance value set determining unit configured to obtain, based on the first functional relationship and the second plane distance set, a second semivariance value set between the interpolation point and each mesh vertex in each initial operating parameter mesh;

[0160] The interpolation point operation parameter value determining unit is configured to calculate the interpolation point operation parameter value corresponding to the interpolation point based on the Kriging interpolation algorithm and the second semi-variance value set.

[0161] In some embodiments, the interpolation point operation parameter value determination unit is specifically configured to:

[0162] Obtaining an operating parameter value corresponding to each grid vertex in each initial operating parameter grid obtained according to the three-dimensional data to obtain an operating parameter set;

[0163] Calculating the expected value and variance value of the operating parameter set;

[0164] Obtaining a random deviation value of a preset interpolation point, and calculating an estimation error of the interpolation point based on the random deviation value, the expected value, and the variance value using a Kriging interpolation algorithm;

[0165] Calculating a weight coefficient of each mesh vertex using the estimation error and the second semivariance value set;

[0166] According to the weight coefficient of each mesh vertex and the operating parameter value corresponding to each mesh vertex, the interpolation point operating parameter value corresponding to the interpolation point is calculated based on the Kriging interpolation algorithm.

[0167] For the convenience of description, the above devices are described as being functionally divided into various modules. Of course, when implementing the present disclosure, the functions of each module can be implemented in the same or multiple software and / or hardware.

[0168] The device of the above embodiment is used to implement the method for determining the operating level of the corresponding hydraulic equipment in any of the above embodiments, and has the beneficial effects of the corresponding method embodiment, which will not be repeated here.

[0169] Based on the same inventive concept, corresponding to any of the above-mentioned embodiments and methods, the present disclosure also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the method for determining the operating level of the hydraulic equipment described in any of the above embodiments is implemented.

[0170] FIG6 shows a more specific schematic diagram of the hardware structure of an electronic device provided in this embodiment. The device may include: a processor 1010, a memory 1020, an input / output interface 1030, a communication interface 1040, and a bus 1050. The processor 1010, the memory 1020, the input / output interface 1030, and the communication interface 1040 are communicatively connected to each other within the device via the bus 1050.

[0171] The processor 1010 can be implemented using a general-purpose CPU (Central Processing Unit), a microprocessor, an application-specific integrated circuit (ASIC), or one or more integrated circuits, and is used to execute relevant programs to implement the technical solutions provided in the embodiments of this specification.

[0172] The memory 1020 can be implemented in the form of ROM (Read Only Memory), RAM (Random Access Memory), static storage devices, dynamic storage devices, etc. The memory 1020 can store an operating system and other application programs. When the technical solutions provided in the embodiments of this specification are implemented through software or firmware, the relevant program code is stored in the memory 1020 and is called and executed by the processor 1010.

[0173] The input / output interface 1030 is used to connect input / output modules to implement information input and output. The input / output modules can be configured as components within the device (not shown in the figure) or can be externally connected to the device to provide corresponding functions. Input devices may include a keyboard, mouse, touch screen, microphone, various sensors, etc., and output devices may include a display, speaker, vibrator, indicator light, etc.

[0174] The communication interface 1040 is used to connect to a communication module (not shown) to enable communication between the device and other devices. The communication module can communicate via a wired method (such as USB, network cable, etc.) or a wireless method (such as mobile network, WiFi, Bluetooth, etc.).

[0175] The bus 1050 comprises a path for transmitting information between the various components of the device (eg, the processor 1010 , the memory 1020 , the input / output interface 1030 , and the communication interface 1040 ).

[0176] It should be noted that although the above device only shows the processor 1010, the memory 1020, the input / output interface 1030, the communication interface 1040, and the bus 1050, in a specific implementation, the device may also include other components necessary for normal operation. In addition, it will be understood by those skilled in the art that the above device may only include the components necessary to implement the embodiments of this specification, and does not necessarily include all the components shown in the figure.

[0177] The electronic device of the above embodiment is used to implement the method for determining the operating level of the corresponding hydraulic equipment in any of the above embodiments, and has the beneficial effects of the corresponding method embodiment, which will not be repeated here.

[0178] Based on the same inventive concept, corresponding to any of the above-mentioned embodiment methods, the present disclosure also provides a non-transitory computer-readable storage medium, wherein the non-transitory computer-readable storage medium stores computer instructions, and the computer instructions are used to enable the computer to execute the method for determining the operating level of the hydraulic equipment as described in any of the above embodiments.

[0179] The computer-readable media of this embodiment include permanent and non-permanent, removable and non-removable media that can be used to store information by any method or technology. The information can be computer-readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, read-only compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device.

[0180] The computer instructions stored in the storage medium of the above embodiment are used to enable the computer to execute the method for determining the operating level of hydraulic equipment as described in any of the above embodiments, and have the beneficial effects of the corresponding method embodiments, which will not be repeated here.

[0181] It is understandable that before using the technical solutions of each embodiment of the present disclosure, the type, scope of use, usage scenarios, etc. of the personal information involved will be informed to the user in an appropriate manner, and the user's authorization will be obtained.

[0182] For example, in response to a user's active request, a prompt message is sent to the user to clearly inform the user that the requested operation will require the acquisition and use of the user's personal information. This allows the user to independently choose whether to provide personal information to the electronic device, application, server, storage medium, or other software or hardware that performs the operation of the disclosed technical solution based on the prompt message.

[0183] As an optional but non-limiting implementation, in response to a user's active request, the prompt information may be sent to the user in the form of a pop-up window, in which the prompt information may be presented in text form. Furthermore, the pop-up window may also contain a selection control for the user to select "agree" or "disagree" to provide personal information to the electronic device.

[0184] It is understandable that the above notification and user authorization process are merely illustrative and do not constitute a limitation on the implementation of the present disclosure. Other methods that comply with relevant laws and regulations may also be applied to the implementation of the present disclosure.

[0185] Those skilled in the art should understand that the discussion of any of the above embodiments is merely illustrative and is not intended to imply that the scope of the present disclosure (including the claims) is limited to these examples. Within the scope of the present disclosure, the technical features in the above embodiments or different embodiments may be combined, the steps may be implemented in any order, and there are many other variations of the different aspects of the embodiments of the present disclosure as described above, which are not provided in detail for the sake of simplicity.

[0186] In addition, to simplify the description and discussion, and so as not to obscure the embodiments of the present disclosure, known power / ground connections to integrated circuit (IC) chips and other components may or may not be shown in the provided figures. In addition, devices may be shown in the form of block diagrams to avoid obscuring the embodiments of the present disclosure, and this also takes into account the fact that the details of the implementation of these block diagram devices are highly dependent on the platform on which the embodiments of the present disclosure are to be implemented (i.e., these details should be fully within the purview of those skilled in the art). Where specific details (e.g., circuits) are set forth to describe exemplary embodiments of the present disclosure, it will be apparent to those skilled in the art that the embodiments of the present disclosure may be implemented without these specific details or with variations in these specific details. Therefore, these descriptions should be considered illustrative rather than restrictive.

[0187] Although the present disclosure has been described in conjunction with specific embodiments thereof, many alternatives, modifications, and variations of these embodiments will be apparent to those skilled in the art based on the foregoing description. For example, other memory architectures (e.g., dynamic RAM (DRAM)) may use the embodiments discussed.

[0188] The embodiments of the present disclosure are intended to cover all such substitutions, modifications, and variations that fall within the broad scope of the appended claims. Therefore, any omissions, modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the embodiments of the present disclosure should be included in the scope of protection of the present disclosure.

Claims

1. A method for determining the operating level of hydraulic equipment based on Kriging interpolation algorithm, characterized in that: include: Acquire three-dimensional data and multiple target operating parameter values ​​corresponding to the operating parameter information of the hydraulic equipment, wherein the three dimensions of the three-dimensional data are the unit speed of the hydraulic equipment, the unit flow of the hydraulic equipment and the operating parameter value of the hydraulic equipment; Mapping the three-dimensional data to a two-dimensional grid corresponding to a unit rotation speed of the hydraulic equipment and a unit flow rate of the hydraulic equipment to obtain at least one initial operation parameter grid, wherein each grid vertex in the initial operation parameter grid is the three-dimensional data; For each initial operating parameter grid, interpolation processing is performed on the initial operating parameter grid to obtain an interpolation point, and the interpolation point operating parameter value corresponding to the interpolation point is calculated based on the Kriging interpolation algorithm, and the interpolation process is repeated until the interpolation point operating parameter value meets the preset parameter condition to obtain an operating parameter grid; For each target operating parameter value, multiple equal value points are determined according to the interpolation points on all operating parameter grids and the target operating parameter value, and curve fitting is performed on the multiple equal value points to obtain the equal value line corresponding to the target operating parameter value; The operation level of the hydraulic equipment is determined according to the contour lines corresponding to each target operation parameter value.

2. The method according to claim 1, characterized in that For each initial run parameter grid, The interpolation point operation parameter value meets the preset parameter conditions, including: Obtaining an interpolation point operating parameter value corresponding to each interpolation point in the initial operating parameter grid; Calculating the difference between the interpolation point operating parameter values ​​of any two interpolation points, and in response to the difference being less than a preset difference threshold, determining that the interpolation point operating parameter value meets a preset parameter condition; or, In response to the interpolation point operating parameter value corresponding to each interpolation point being less than the target operating parameter value, determining that the interpolation point operating parameter value meets a preset parameter condition; or, In response to the interpolation point operating parameter value corresponding to each interpolation point being greater than the target operating parameter value, it is determined that the interpolation point operating parameter value meets a preset parameter condition.

3. The method according to claim 1, characterized in that The step of determining a plurality of equal value points according to the interpolation points on all the operating parameter grids and the target operating parameter values ​​comprises: For each run parameter grid: Acquire all interpolation points included in the operating parameter grid, and select any two interpolation points from all the interpolation points as the first target interpolation point and the second target interpolation point; Determine the magnitude relationship between a first interpolation point operating parameter value corresponding to the first target interpolation point and a second interpolation point operating parameter value corresponding to the second target interpolation point and the target operating parameter value; In response to the existence of an interpolation point operating parameter value greater than the target operating parameter value and an interpolation point operating parameter value less than the target operating parameter value among the first interpolation point operating parameter value and the second interpolation point operating parameter value, determining that an equal value point exists between the first target interpolation point and the second target interpolation point; Calculating a first proportional relationship according to the first interpolation point operating parameter value, the second interpolation point operating parameter value and the target operating parameter value; Obtaining a first unit speed and a first unit flow rate corresponding to the first target interpolation point, and a second unit speed and a second unit flow rate corresponding to the second target interpolation point, and calculating the target unit speed and the target unit flow rate corresponding to the equivalent point based on the first proportional relationship according to the first unit speed, the first unit flow rate, the second unit speed and the second unit flow rate; An equivalue point is obtained according to the target unit rotation speed, the target unit flow rate and the target operating parameter value.

4. The method according to claim 1, characterized in that: After getting the running parameter grid, it also includes: Obtaining the origin of the two-dimensional grid corresponding to the unit speed of the hydraulic equipment and the unit flow of the hydraulic equipment, and labeling the operation parameter grid according to the horizontal axis and the vertical axis of the two-dimensional grid to obtain the label information of the operation parameter grid; The performing curve fitting on the multiple isovalue points to obtain isovalue lines corresponding to the target operating parameter values ​​includes: For each operating parameter grid, obtaining equal value points contained in the operating parameter grid; Along a preset direction, all the equal-value points contained in each operating parameter grid are connected in sequence according to the label information of the operating parameter grid to obtain the equal-value line corresponding to the target operating parameter value.

5. The method according to claim 1, characterized in that For each initial run parameter grid, The interpolation processing is performed on the initial operation parameter grid to obtain interpolation points, including: Obtaining four corresponding grid vertices in the initial operating parameter grid; Select the middle point of the line connecting any two adjacent mesh vertices as the first initial interpolation point; Selecting an intersection point of lines connecting two diagonal mesh vertices, and using the intersection point as a second initial interpolation point; The first initial interpolation points and the second initial interpolation points are counted to obtain interpolation points corresponding to the initial operating parameter grid.

6. The method according to claim 1, characterized in that For each initial run parameter grid, The calculating the interpolation point operation parameter value corresponding to the interpolation point based on the Kriging interpolation algorithm includes: Calculating the first plane distance and the first semi-variance value between any two mesh vertices in each initial operating parameter mesh obtained according to the three-dimensional data to obtain a first plane distance set and a first semi-variance value set; Fitting the first plane distance set and the first semi-variance value set by using the least square method to obtain a first functional relationship; Calculating the second plane distance between the interpolation point and each grid vertex in each initial operating parameter grid to obtain a second plane distance set; Based on the first functional relationship and the second plane distance set, obtaining a second semivariance value set between the interpolation point and each mesh vertex in each initial operating parameter mesh; An interpolation point operating parameter value corresponding to the interpolation point is calculated based on the Kriging interpolation algorithm and the second semi-variance value set.

7. The method according to claim 6, characterized in that The operation parameter value of the interpolation point is calculated based on the Kriging interpolation algorithm and the second semi-variance value set, including: Obtaining an operating parameter value corresponding to each mesh vertex in each initial operating parameter mesh obtained according to the three-dimensional data to obtain an operating parameter set; Calculating the expected value and variance value of the operating parameter set; Obtaining a random deviation value of a preset interpolation point, and calculating an estimation error of the interpolation point based on the Kriging interpolation algorithm according to the random deviation value, the expected value and the variance value; The weight coefficient of each mesh vertex is calculated using the estimation error and the second semi-variance value set; According to the weight coefficient of each mesh vertex and the operating parameter value corresponding to each mesh vertex, the interpolation point operating parameter value corresponding to the interpolation point is calculated based on the Kriging interpolation algorithm.

8. A device for determining the operating level of hydraulic equipment based on Kriging interpolation algorithm, characterized in that: include: A parameter information acquisition module is configured to acquire three-dimensional data and a plurality of target operating parameter values ​​corresponding to the operating parameter information of the hydraulic equipment, wherein the three dimensions of the three-dimensional data are a unit speed of the hydraulic equipment, a unit flow rate of the hydraulic equipment, and an operating parameter value of the hydraulic equipment; A parameter grid determination module is configured to map the three-dimensional data to a two-dimensional grid corresponding to a unit rotation speed of the hydraulic equipment and a unit flow rate of the hydraulic equipment to obtain at least one initial operation parameter grid, wherein each grid vertex in the initial operation parameter grid is the three-dimensional data; The interpolation processing module is configured to perform interpolation processing on each initial operating parameter grid to obtain an interpolation point, calculate the interpolation point operating parameter value corresponding to the interpolation point based on the Kriging interpolation algorithm, and repeat the interpolation process until the interpolation point operating parameter value meets the preset parameter condition to obtain an operating parameter grid; An isovalue point determination module is configured to determine, for each target operating parameter value, a plurality of isovalue points according to interpolation points on all operating parameter grids and the target operating parameter value, and perform curve fitting on the plurality of isovalue points to obtain an isovalue line corresponding to the target operating parameter value; The operation level determination module is configured to determine the operation level of the hydraulic equipment according to the contour lines corresponding to each target operation parameter value.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the method for determining the operating level of a hydraulic device as claimed in any one of claims 1 to 7 is implemented.

10. A non-transitory computer-readable storage medium, characterized in that: The non-transitory computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a computer to execute the method for determining the operating level of a hydraulic equipment according to any one of claims 1 to 7.

Images