Method and device for setting up a micro-irrigation pipe network
By dividing irregular plots into irrigation triangles and optimizing the location of water outlet piles based on hydraulic and economic parameters, and by setting up branch pipes and capillary pipes, the problems of irrigation uniformity and economy in irregular plots are solved, achieving efficient and low-cost irrigation results.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2022-09-30
- Publication Date
- 2026-03-20
AI Technical Summary
In existing micro-irrigation network systems for irregular plots, it is difficult to simultaneously satisfy both irrigation uniformity and economy, and the terrain deviation has a significant impact, resulting in uneven irrigation effects.
Irregular plots are divided into multiple irrigation triangles. The location of water outlet piles is optimized by using hydraulic performance and economic evaluation parameters. Branch pipes and capillary pipes are installed to ensure irrigation uniformity and economy.
It improves irrigation uniformity, reduces costs, and ensures consistent irrigation results in irregular plots by limiting capillary length differences through triangular boundaries.
Smart Images

Figure CN115688334B_ABST
Abstract
Description
TECHNICAL FIELD
[0001] The present application relates to the field of agricultural irrigation technology, and in particular to an irrigation pipe network setting method and device. BACKGROUND
[0002] Micro-irrigation is an irrigation method that delivers water and nutrients to the soil around the plant roots in small flow rates according to the crop requirements, through a pipe network system and emitters on the last stage of the pipe network. Under the influence of the hydraulic deviation of the pipe network system, the manufacturing deviation of the emitters and the topographic deviation, etc., the irrigation effect of the micro-irrigation pipe network system is not uniform to a certain extent.
[0003] In the related art, the flow rate is allocated on the branch pipe and the capillary tube based on the allowable pressure deviation in the micro-irrigation unit of the micro-irrigation pipe network through the multi-hole coefficient method, the energy width line method, the backtracking method or the finite element method, etc., and then the branch pipe position, the branch / capillary tube length and diameter, etc. are optimized and calculated.
[0004] However, the uniformity of the irrigation of the micro-irrigation pipe network system for irregular plots still needs to be improved. SUMMARY
[0005] In view of this, the present application provides a micro-irrigation pipe network setting method and device, which can simplify the pipe network setting process of irregular plots and ensure that the irrigation uniformity is high and the cost is low.
[0006] Specifically, the technical scheme includes the following:
[0007] In a first aspect, a micro-irrigation pipe network setting method is provided, and the method includes:
[0008] dividing an irregular plot into a plurality of irrigation triangles, determining a plurality of optional positions of a water stake in each irrigation triangle, calculating a hydraulic performance evaluation parameter corresponding to each optional position, determining a plurality of candidate positions from the plurality of optional positions based on the hydraulic performance evaluation parameter, and the hydraulic performance evaluation parameter is used to indicate the irrigation uniformity when the plot is irrigated using a micro-irrigation pipe network;
[0009] calculating an economic evaluation parameter corresponding to each candidate position of the water stake in each irrigation triangle, and the economic evaluation parameter is used to indicate the cost required for irrigating the plot using the micro-irrigation pipe network;
[0010] determining a target setting position of the water stake in each irrigation triangle from the plurality of candidate positions based on the economic evaluation parameter;
[0011] Based on the target setting position, the water outlet stake, branch pipes and capillary pipes are set in each irrigation triangle, the branch pipes are connected with the water outlet stake, the capillary pipes are connected with the branch pipes, and a plurality of water emitters are uniformly arranged on the capillary pipes.
[0012] Optionally, the step of dividing the irregular land block into a plurality of irrigation triangles comprises:
[0013] n vertices of the irregular land block are obtained, n being a positive integer greater than 3;
[0014] each vertex of the n vertices is sequentially connected with n-2 vertices adjacent to the vertex, and n division schemes each comprising n-2 triangles sharing one vertex are obtained;
[0015] the area of each triangle in each division scheme is calculated by a first preset relationship, the first preset relationship being:
[0016] wherein, S i represents the area of the ith triangle, (x1, y1), (x2, y2) and (x3, y3) represent the coordinates of the three vertices of the ith triangle respectively;
[0017] the absolute value of the area difference between any two triangles in each division scheme is calculated by a second preset relationship, the second preset relationship being:
[0018] D ij = |S i -S j |, wherein, i = 1:(n-2), j = 1:(n-2), i < j, S j represents the area of the jth triangle, S i represents the area of the ith triangle, D ij represents the area difference between the ith triangle and the jth triangle;
[0019] the area evaluation parameter of each division scheme is calculated by a third preset relationship, the third preset relationship being:
[0020] wherein, i = 1:(n-2), j = 1:(n-2), i < j, D a represents the area evaluation parameter;
[0021] the n-2 triangles corresponding to the division scheme with the minimum area evaluation parameter are determined as first-level triangles;
[0022] the first-level triangles or second-level triangles obtained by dividing the first-level triangles are used as the irrigation triangles.
[0023] Optionally, the hydraulic performance evaluation parameter is a deviation rate of emitter flow in the irrigation triangle.
[0024] The method for determining the candidate position of the water stake in the irregular land, comprising:
[0025] Step a1: setting an initial position of the water stake in the primary triangle, and setting a diameter of the branch pipe as a first diameter, the initial position being any position in the primary triangle, the first diameter being an i-th diameter in optional diameters of the branch pipe, the optional diameters including n optional diameters sequentially increasing with the serial number, i being a positive integer greater than or equal to 1 and less than or equal to n, and the initial value being 1;
[0026] Step a2: determining a plurality of optional positions of the water stake in the primary triangle based on the initial position and a position adjustment step;
[0027] Step a3: calculating a deviation rate of emitter flow in the primary triangle when the water stake is located at each optional position based on each optional position and the first diameter;
[0028] Step a4: determining whether there is at least one candidate position in the plurality of optional positions based on the deviation rate of emitter flow in the primary triangle and a deviation rate condition, if yes, executing step a5, otherwise executing a6;
[0029] Step a5: determining the primary triangle as the irrigation triangle, determining the first diameter as a target diameter of the branch pipe in the primary triangle, and determining all candidate positions from the plurality of optional positions;
[0030] Step a6: adjusting the first diameter of the branch pipe to an i+1-th diameter in the optional diameters, and jumping to step a3 until i+1 is greater than n.
[0031] Optionally, in the step a6, when i+1 is greater than n, the method for determining the candidate position of the water stake in the irregular land, further comprising:
[0032] Step a7: dividing the primary triangle into M secondary triangles with equal areas, M being a positive integer greater than or equal to 2, and the initial value being 2;
[0033] Step a8: determining whether there is at least one candidate position in each of the secondary triangles in sequence, if there is at least one candidate position in each of the secondary triangles, executing step a9, otherwise executing step a10;
[0034] Step a9: determining each of the secondary triangles as the irrigation triangle, and determining all candidate positions in the multiple optional positions of the water outlet stake in each of the secondary triangles in sequence;
[0035] Step a10: assigning M+1 to M, and jumping to step a7.
[0036] Optionally, the economic evaluation parameter is irrigation annual cost, and the calculation of the economic evaluation parameter corresponding to each candidate position of the water outlet stake in each of the irrigation triangles comprises:
[0037] calculating the pipe network system investment of the irrigation triangle through a fourth preset relationship, the fourth preset relationship being:
[0038] wherein i represents annual interest rate, m represents service life of branch pipe and main pipe, L l and L m are length of main pipe and branch pipe in the irrigation triangle respectively, P l and P m are unit price of main pipe and branch pipe in the irrigation triangle respectively, S is area of the irrigation triangle, and C a represents the pipe network system investment;
[0039] multiplying the pipe network system investment C a by a first proportion to obtain pipe network maintenance cost C m of the irrigation triangle;
[0040] calculating the energy consumption cost of the irrigation triangle through a fifth preset relationship, the fifth preset relationship being:
[0041] wherein C e represents the energy consumption cost, N p is pump power, O t is annual irrigation time for meeting irrigation demand, En c is electricity fee, Q s is branch pipe inlet flow, H0 is branch pipe inlet pressure, E p is pump efficiency, R n is net irrigation quota, and EU is irrigation uniformity coefficient, which is calculated through a sixth preset relationship, the sixth preset relationship being:
[0042] wherein CV mN is the number of emitters allocated to each crop, q is the average flow rate of emitters in the irrigation triangle, min qmin is the minimum flow rate of emitters in the irrigation triangle;
[0043] The annual water fee of the irrigation triangle is calculated by a seventh preset relationship, which is:
[0044] C w = R g P w , wherein C w represents the annual water fee, P w is the water price, and R g is the gross irrigation quota during the growth period of crops, which is calculated by an eighth preset relationship, which is:
[0045] The sum of the pipe network system investment C a , the pipe network maintenance cost C m , the energy consumption cost C e , and the annual water fee C w is taken as the annual irrigation cost of the irrigation triangle;
[0046] The target setting position of the water outlet stake in each irrigation triangle is determined from the plurality of candidate positions based on the economic evaluation parameters, which includes:
[0047] The candidate position with the minimum annual irrigation cost in each irrigation triangle is determined as the target setting position of the water outlet stake.
[0048] Optionally, the method further includes:
[0049] A plurality of candidate branch pipe connection points are determined in a first triangle, which is a primary triangle divided into a plurality of secondary triangles;
[0050] A first distance from each candidate branch pipe connection point to any water outlet stake in the first triangle is calculated by a ninth preset relationship, which is:
[0051] wherein i = 1:M, D Noi represents the first distance from the candidate branch pipe connection point (x N , y N ) to the i-th water outlet stake , and M represents the total number of water outlet stakes in the first triangle;
[0052] The absolute value of the difference between any two first distances is calculated by a tenth preset relationship, which is:
[0053] S ij =|D NOi -D NOj | where i = 1:M, j = 1:M, i < j, S ij This represents the absolute value of the difference between the i-th first distance and the j-th first distance;
[0054] The first pressure balance evaluation parameter corresponding to each candidate branch pipe connection point is calculated using the eleventh preset relation, which is:
[0055] Where i = 1:M, j = 1:M, i < j, S a This represents the first pressure equilibrium evaluation parameter;
[0056] The candidate branch pipe connection point with the minimum first pressure equalization evaluation parameter is determined as the target branch pipe connection point;
[0057] A branch pipe is installed between each water outlet pile and the target branch pipe connection point in the first triangle.
[0058] Optionally, the method further includes:
[0059] Multiple candidate trunk connection points were identified in the irregular plot of land;
[0060] The twelfth preset relation is used to calculate multiple second distances corresponding to each candidate main pipe connection point. These second distances include the distance between the candidate main pipe connection point and any branch pipe connection point, or the distance between the candidate main pipe connection point and any outlet pile in a second triangle. The second triangle is a first-level triangle that has not been divided into second-level triangles. The twelfth preset relation is:
[0061] Where i = 1:n-2, D' Noi Representative of candidate trunk connection point (x') N y' N ) to time The second distance, point (x') oi y' oi ) represents the location of any branch pipe connection point or any outlet pile in the second triangle, and n represents the number of vertices of the irregular plot;
[0062] The absolute value of the difference between any two second distances is calculated using the thirteenth preset relation, which is:
[0063] S' ij =|D' NOi -D' NOj|, wherein i=1:n-2, j=1:n-2, i ij an absolute value representing a difference between the ith second distance and the jth second distance;
[0064] calculating the second pressure balance evaluation parameter corresponding to each candidate main pipe connection point by a fourteenth preset relationship, the fourteenth preset relationship being:
[0065] wherein i=1:n-2, j=1:n-2, i a representing the second pressure balance evaluation parameter;
[0066] determining the candidate main pipe connection point with the smallest second pressure balance evaluation parameter as a target main pipe connection point;
[0067] arranging a main pipe in the irregular land plot based on the target main pipe connection point.
[0068] Optionally, the arranging the water outlet stake, the branch pipe and the capillary pipe in each irrigation triangle comprises:
[0069] extending the branch pipe from a position where the water outlet stake is located in the irrigation triangle to the boundary of the irrigation triangle along a first direction, the first direction being perpendicular to an extension direction of the capillary pipe, the extension direction of the capillary pipe being determined according to a planting direction of crops;
[0070] when the branch pipe extends to intersect a first side of the irrigation triangle, extending the branch pipe along a second direction until reaching a vertex of the irrigation triangle, the first side being any side of the irrigation triangle, the second direction being parallel to the first side, and an included angle between the second direction and the first direction being an acute angle;
[0071] when the branch pipe extends to intersect a first vertex of the irrigation triangle, the branch pipe being cut off, the first vertex being any vertex of the irrigation triangle.
[0072] Optionally, the method further comprises:
[0073] obtaining elevation information of a plurality of elevation measurement points in the irregular land plot;
[0074] obtaining elevation information of all water irrigators in the irregular land plot by interpolation based on the elevation information of the plurality of elevation measurement points;
[0075] the calculating the deviation rate of the flow of the water irrigator in the primary triangle when the water outlet stake is located at each optional position comprises:
[0076] calculate a flow rate of each sprinkler in the first triangle based on the elevation information of all sprinklers in the irregular field, a pipe segment energy equation, a pipe network node flow continuity equation, and a sprinkler flow-pressure relationship;
[0077] calculate a flow rate of each sprinkler in the first triangle based on the elevation information of all sprinklers in the irregular field, a pipe segment energy equation, a pipe network node flow continuity equation, and a sprinkler flow-pressure relationship;
[0078] In a second aspect, a micro-irrigation pipe network setting device is provided, and the device comprises:
[0079] a determination module configured to divide an irregular field into a plurality of irrigation triangles, determine a plurality of optional positions of a water outlet post in each of the irrigation triangles, calculate a hydraulic performance evaluation parameter corresponding to each of the optional positions, and determine a plurality of candidate positions from the optional positions based on the hydraulic performance evaluation parameter, the hydraulic performance evaluation parameter being used to indicate irrigation uniformity when the field is irrigated by using a micro-irrigation pipe network;
[0080] a calculation module configured to calculate an economic evaluation parameter corresponding to each of the candidate positions of the water outlet post in each of the irrigation triangles, the economic evaluation parameter being used to indicate a cost required for irrigating the field by using the micro-irrigation pipe network;
[0081] the determination module is further configured to determine a target setting position of the water outlet post in each of the irrigation triangles from the candidate positions based on the economic evaluation parameter;
[0082] a setting module configured to set the water outlet post, a branch pipe, and a micro tube in each of the irrigation triangles based on the target setting position, the branch pipe being connected to the water outlet post, the micro tube being connected to the branch pipe, and the micro tube being uniformly provided with a plurality of sprinklers.
[0083] Embodiments of the present application provide a micro-irrigation pipe network setting method and device. In the method, an irregular field is divided into a plurality of irrigation triangles, and the irrigation triangle is used as a pipe network arrangement unit to effectively simplify the pipe network setting process of an arbitrary irregular field. For each irrigation triangle, the position of the water outlet post is optimized based on a hydraulic performance evaluation parameter and an economic evaluation parameter, so that the pipe network in each irrigation triangle can meet the hydraulic performance requirement and the economic requirement, thereby effectively improving irrigation uniformity and reducing cost. In addition, the micro tube is arranged based on a smaller triangle, which can limit the difference degree of the length of different micro tubes in the field by means of the boundary of the triangle, so as to ensure that the length difference of the micro tube in the irregular field is small, and irrigation uniformity is improved. BRIEF DESCRIPTION OF DRAWINGS
[0084] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings needed to be used in the embodiments description will be briefly introduced.
[0085] Figure 1 A flow chart of a micro-irrigation pipe network setting method provided by the embodiments of the present application;
[0086] Figure 2 A flow chart of a method for determining an irrigation triangle provided by the embodiments of the present application;
[0087] Figure 3 A schematic diagram of an irregular plot provided by the embodiments of the present application;
[0088] Figure 4 Schematic diagrams of different division schemes for dividing an irregular plot provided by the embodiments of the present application;
[0089] Figure 5 A flow chart of a method for arranging branch pipes provided by the embodiments of the present application;
[0090] Figure 6 A schematic diagram of arranging a capillary tube in an irregular plot provided by the embodiments of the present application;
[0091] Figure 7 A schematic diagram of arranging a branch pipe in an irrigation triangle provided by the embodiments of the present application;
[0092] Figure 8 A flow chart of a method for setting a water outlet stake provided by the embodiments of the present application;
[0093] Figure 9 A flow chart of a method for determining a candidate position of a water outlet stake in an irrigation pipe network setting method provided by the embodiments of the present application;
[0094] Figure 10 A flow chart of a method for determining elevation information of an emitter provided by the embodiments of the present application;
[0095] Figure 11 A flow chart of a method for calculating emitter flow deviation rate in a triangle provided by the embodiments of the present application;
[0096] Figure 12 A schematic diagram of emitter flow distribution in a triangle provided by the embodiments of the present application;
[0097] Figure 13 A schematic diagram of M-equal division of a triangle provided by the embodiments of the present application;
[0098] Figure 14 A flow chart of a method for setting a sub-main pipe provided by the embodiments of the present application;
[0099] Figure 15 A schematic diagram of setting a branch pipe provided by an embodiment of the present application;
[0100] Figure 16 A schematic diagram of the arrangement of water outlet piles and branch pipes in an irregular plot provided by an embodiment of the present application;
[0101] Figure 17 A method flowchart of setting a branch pipe provided by an embodiment of the present application;
[0102] Figure 18 A schematic diagram of the distribution of flow deviation rates of emitters corresponding to different candidate positions in a primary triangle provided by an embodiment of the present application;
[0103] Figure 19 A schematic diagram of the elevation distribution of emitters in an irregular plot provided by an embodiment of the present application;
[0104] Figure 20 A schematic diagram of the arrangement of branch pipes in an irregular plot provided by an embodiment of the present application;
[0105] Figure 21 A structural diagram of a micro-irrigation pipe network setting device provided by an embodiment of the present application.
[0106] The specific embodiments of the present application have been shown in the above-described drawings, and will be described in more detail hereinafter. These drawings and textual descriptions are not intended to limit the scope of the concept of the present application in any way, but to illustrate the concept of the present application to those skilled in the art by referring to specific embodiments. DETAILED DESCRIPTION
[0107] The technical solutions in the embodiments of the present application will be described clearly and completely below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only some of the embodiments of the present application, rather than all the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without any creative work fall within the scope of protection of the present application.
[0108] It should be noted that the present application is constructed on the basis that the micro-irrigation pipe network design method in the related art has the following defects:
[0109] (1) When the related art performs hydraulic calculation with the constraint condition of flow deviation rate ≤ 20%, the branch pipe diameter is mainly selected according to experience, although the uniform performance of irrigation is satisfied, but the economic indicators such as micro-irrigation unit pipe network investment, maintenance cost, energy consumption cost and water fee are not comprehensively considered, and it is unknown whether the system is economically optimal;
[0110] (2) Related technologies set the distribution ratio of the micro-irrigation unit allowable pressure deviation on the branch pipe and the capillary pipe as a fixed value, and then calculate the maximum branch / capillary pipe length of the micro-irrigation unit. Affected by the size of the field, the actual branch / capillary pipe length is less than the maximum branch / capillary pipe length and presents more random combination states. After the micro-irrigation unit is arranged, it is unknown how the actual irrigation uniformity is, and the system cannot achieve economic optimization while meeting uniformity;
[0111] (3) Related technologies calculate the maximum branch / capillary pipe length for regular units, but the capillary pipe length in irregular units presents inconsistent states due to the influence of field boundaries. The inconsistency of the capillary pipe length affects the hydraulic performance and uniformity in the micro-irrigation unit. Related technologies do not have clear instructions for the optimal design of the pipe network of irregular micro-irrigation units.
[0112] (4) Related technologies consider the influence of terrain deviation in the form of uniform slope (downhill, flat, and uphill) in the design of micro-irrigation systems. However, the terrain fluctuation characteristics in actual micro-irrigation units are more random, and the slope changes more frequently. Related technologies do not adequately consider terrain deviation under non-uniform slope conditions, making it difficult to consider the actual influence of terrain deviation on the hydraulic performance of the micro-irrigation unit in the hydraulic design of the micro-irrigation unit.
[0113] Therefore, in the field of irrigation, especially in the design of micro-irrigation pipe networks for irregular fields, there is an urgent need for a pipe network optimization arrangement method that considers the irregular terrain characteristics of the field, the hydraulic performance of the unit, and the economy comprehensively.
[0114] In the embodiments of the present application, the micro-irrigation pipe network system can include multiple levels, such as a main pipe, a branch pipe, a branch pipe, and a capillary pipe connected in sequence. The capillary pipe at the end of the pipe network system can be arranged with an irrigation device, which can spray or drip water or fertilizer into the soil around the crops.
[0115] In order to make the technical solutions and advantages of the present application clearer, the embodiments of the present application will be described in further detail below with reference to the drawings.
[0116] In a first aspect, with reference to Figure 1 The embodiments of the present application provide a micro-irrigation pipe network setting method. The method includes steps 101-104.
[0117] Step 101: divide the irregular field into a plurality of irrigation triangles, determine a plurality of optional positions of the water stake in each irrigation triangle, calculate a hydraulic performance evaluation parameter corresponding to each optional position, and determine a plurality of candidate positions from the plurality of optional positions based on the hydraulic performance evaluation parameter. The hydraulic performance evaluation parameter is used to indicate the irrigation uniformity when the field is irrigated using the micro-irrigation pipe network.
[0118] In this step, the irregular plot can be divided to obtain a plurality of irrigation triangles, and the sum of the plurality of irrigation triangles is equal to the irregular plot. In this way, the micro-irrigation pipe network setting problem of an arbitrary irregular plot is converted into a micro-irrigation pipe network setting problem of a triangular plot. On this basis, the hydraulic performance evaluation parameter corresponding to each of a plurality of optional positions in each irrigation triangle is calculated, and a plurality of candidate positions are determined from the plurality of optional positions based on the hydraulic performance evaluation parameter, so as to realize preliminary optimization of the water outlet pile position based on the hydraulic performance evaluation parameter. This can ensure that the final setting position of the water outlet pile in the irrigation triangle is necessarily a position that can meet the hydraulic performance requirement.
[0119] Step 102: Calculate the economic evaluation parameter corresponding to each candidate position of the water outlet pile in each irrigation triangle, which is used to indicate the cost required for irrigating the plot using the micro-irrigation pipe network.
[0120] In this step, the economic evaluation parameter corresponding to each candidate position is calculated to evaluate the economy of different candidate positions, so as to facilitate further optimization of the water outlet pile position.
[0121] Step 103: Determine the target setting position of the water outlet pile in each irrigation triangle from the plurality of candidate positions based on the economic evaluation parameter.
[0122] In this step, the target setting position of the water outlet pile in each irrigation triangle is determined from the candidate positions by means of the economic evaluation parameter, so that the finally determined water outlet pile position can not only meet the hydraulic performance requirement, but also meet the economic requirement.
[0123] Step 104: Based on the target setting position, set the water outlet pile, branch pipe and micro tube in each irrigation triangle, the branch pipe is connected with the water outlet pile, the micro tube is connected with the branch pipe, and a plurality of water emitters are uniformly arranged on the micro tube.
[0124] After obtaining the target setting position of the water outlet pile in each micro-irrigation triangle in step 103, each micro-irrigation triangle can be used as a micro-irrigation pipe network design unit to arrange the water outlet pile, branch pipe and micro tube. The length of the micro tube is limited by the boundary of the micro-irrigation triangle and is kept within a certain range, so that the length difference of the micro tube in the entire irregular plot is also stabilized within a certain range, which can alleviate the problem of poor irrigation uniformity caused by excessive length difference of the micro tube.
[0125] In summary, this application provides a method for setting up a micro-irrigation network, effectively simplifying the network setup process for any irregular plot. Furthermore, for each irrigation triangle obtained from the division of the irregular plot, the position of the outlet piles within each triangle is optimized secondaryly based on hydraulic performance and economic evaluation parameters. This ensures that the network setup within each irrigation triangle meets both hydraulic performance and economic requirements, resulting in high irrigation uniformity and low cost. In addition, this method utilizes the boundaries of the triangles to limit the variation in capillary tube lengths within the plot, ensuring minimal variation in capillary tube lengths in irregular plots and further improving irrigation uniformity.
[0126] In some embodiments, reference Figure 2 Dividing irregular plots into multiple irrigation triangles may include steps 201-207.
[0127] Step 201: Obtain n vertices of the irregular plot, where n is a positive integer greater than 3.
[0128] In this step, obtaining the n vertices of the irregular plot can include obtaining the coordinate values of the n vertices of the irregular plot in a pre-established coordinate system. For example, refer to... Figure 3 Given an irregular plot of land, consisting of five vertices, each corresponding to a set of coordinates, the boundary information of the irregular plot can be represented as an n x 3 matrix B:
[0129] Where n is the number of vertices of the irregular plots, in units of . Figure 3 In the case where n is 5, x i and y i These are the coordinates in the x and y directions, respectively.
[0130] Step 202: Connect each of the n vertices to the n-2 non-adjacent vertices of each vertex in turn to obtain n partition schemes. Each partition scheme includes n-2 triangles with a common vertex.
[0131] refer to Figure 4 It gave Figure 3 The diagram shows three partitioning schemes corresponding to the irregular pentagonal plot. For example, for vertex A, its two non-adjacent vertices are vertices C and D. Connecting vertices A, C, and D respectively yields triangles ABC, ACD, and AED.
[0132] Step 203: Calculate the area of each triangle under each partitioning scheme using the first preset relational formula, which is:
[0133] Among them, S iS i represents the area of the i-th triangle, (x1, y1), (x2, y2) and (x3, y3) represent the coordinates of the three vertices of the i-th triangle respectively.
[0134] In this step, the area of each triangle in the n-2 triangles corresponding to any one division scheme can be calculated by means of the first preset relationship.
[0135] Step 204: Calculate the absolute value of the area difference of any two triangles under each division scheme by a second preset relationship, and the second preset relationship is:
[0136] D ij = |S i -S j |, wherein i = 1: (n-2), j = 1: (n-2), i < j, S j represents the area of the j-th triangle, S i represents the area of the i-th triangle, D ij represents the area difference of the i-th triangle and the j-th triangle.
[0137] In this step, the absolute value of the area difference of any two triangles in the n-2 triangles corresponding to any one division scheme can be calculated by means of the second preset relationship.
[0138] Step 205: Calculate the area evaluation parameter of each division scheme by a third preset relationship, and the third preset relationship is:
[0139] wherein i = 1: (n-2), j = 1: (n-2), i < j, D a represents the area evaluation parameter.
[0140] In this step, the area evaluation parameter of each division scheme can be obtained by means of the third preset relationship, and the area evaluation parameter reflects the difference degree of the n-2 triangles corresponding to each division scheme.
[0141] Step 206: Determine the n-2 triangles corresponding to the division scheme with the smallest area evaluation parameter as the primary triangles.
[0142] In this step, the area evaluation parameters corresponding to a plurality of division schemes are compared, and the n-2 triangles corresponding to the division scheme with the smallest area evaluation parameter are determined as the primary triangles, that is, the scheme with the smallest area difference of the divided triangles is taken as the target scheme, which facilitates to ensure the uniformity of subsequent pipe network irrigation.
[0143] Step 207: Take the primary triangle or the secondary triangle obtained by dividing the primary triangle as the irrigation triangle.
[0144] In this step, the first-order triangle can be taken as an irrigation triangle for subsequent pipe network setting, so as to convert the pipe network setting problem of an arbitrary irregular plot into a pipe network setting problem of a triangular plot, and since the area evaluation parameter of the corresponding division scheme of the first-order triangle is minimum, the area division of the irregular plot is more balanced, which is beneficial to improve the uniformity of irrigation. Alternatively, the first-order triangle can be further divided to obtain smaller second-order triangles as irrigation triangles, which can make the length difference of the emitters in the irregular plot smaller, further improving the uniformity of irrigation.
[0145] The embodiment gives a method for dividing an irregular plot to obtain irrigation triangles, but it should be noted that other ways can also be used for division in practice, and the above embodiments do not constitute a limitation on the present application.
[0146] In some embodiments, the aforementioned hydraulic performance evaluation parameter is the deviation rate of the emitter flow in the irrigation triangle. Referring to Figure 5 , step 101 in the foregoing embodiment can include steps 501-503.
[0147] Step 501: The branch pipe is extended from the position of the water outlet stake in the irrigation triangle along a first direction to the boundary of the irrigation triangle, the first direction is perpendicular to the extension direction of the emitter, and the extension direction of the emitter is determined according to the planting direction of the crops.
[0148] Referring to Figures 6-7 , Figure 6 The schematic diagram of the arrangement of the emitters and the emitters in the irregular plot shown in Figure 3 is given, Figure 7 The schematic diagram of the arrangement of the emitters and the emitters in the triangular plot ABC in Figure 6 is given. Among them, Figures 6-7 The arrangement direction of the emitters in is only an example, and the arrangement direction of the emitters is generally related to the planting direction of the crops. The emitters are arranged at a certain interval, which can be determined according to the row spacing and plant spacing of the crops. The emitters are arranged at a certain interval on the emitters, and the interval of the emitters can also be determined according to the row spacing and plant spacing of the crops (the corresponding relationship between the row spacing and plant spacing of the crops and the arrangement interval of the emitters or the arrangement interval of the emitters can be pre-stored, and the arrangement interval of the emitters or the arrangement interval of the emitters can be determined by looking up the table later).
[0149] In this step, the O point of Figure 7 is taken as an example of the position of the water outlet stake, and when the branch pipe is set, the branch pipe is extended from the O point along a first direction to the boundary of the irrigation triangle, i.e., the OO' pipe section and the OO" pipe section are obtained, and the first direction includes and The first direction is perpendicular to the extension direction of the emitter.
[0150] Step 502: when the branch pipe extends to intersect with a first side of the irrigation triangle, the branch pipe is arranged to extend in a second direction until reaching a vertex of the irrigation triangle, the first side being any side of the irrigation triangle, the second direction being parallel to the first side, and the second direction having an acute angle with the first direction.
[0151] With reference to the foregoing Figure 7 , when the branch pipe extends to intersect with the first side AB or BC of the irrigation triangle, the branch pipe is arranged to extend in a second direction until reaching a vertex of the irrigation triangle, i.e. Figure 7 the middle OO' pipe segment changes the extending direction at point O' to and obtains pipe segment O'A, the middle OO" pipe segment changes the extending direction at point O" to and obtains pipe segment O"C. The second direction has an acute angle with the first direction , the second direction has an acute angle with the first direction , in other words, ∠OO'A and ∠OO"C are obtuse angles. In this way, all the capillary pipes can be connected to the branch pipe, and all the areas in the irrigation triangle can be effectively irrigated.
[0152] Step 503: when the branch pipe extends to intersect with a first vertex of the irrigation triangle, the branch pipe is cut off, the first vertex being any vertex of the irrigation triangle.
[0153] This step can be a parallel step of step 502. In some cases, when the branch pipe extends in the first direction, it can directly intersect with a vertex of the irrigation triangle, and the branch pipe can be directly cut off at the vertex.
[0154] This embodiment provides a way of arranging a branch pipe in an irrigation triangle, which can connect all the capillary pipes to the branch pipe.
[0155] In some embodiments, the economic evaluation parameter is the annual irrigation cost, and with reference to Figure 8 , the foregoing step 102 can include:
[0156] Step 801: calculating the pipe network system investment of the irrigation triangle by a fourth preset relationship formula, the fourth preset relationship formula being:
[0157] wherein i represents the annual interest rate, m represents the service life of the branch pipe and the capillary pipe, L l and L m are the lengths of the capillary pipe and the branch pipe in the irrigation triangle respectively, P l and P m are the unit prices of the capillary pipe and the branch pipe in the irrigation triangle respectively, S is the area of the irrigation triangle, and C a represents the pipe network system investment.
[0158] Step 802: Invest C in the pipeline system a Multiply by the first ratio to obtain the pipeline maintenance cost C of the irrigation triangle. m The first proportion could be, for example, 5%.
[0159] Step 803: Calculate the energy consumption cost of the irrigation triangle using the fifth preset relation. The fifth preset relation is:
[0160] Among them, C e N represents energy consumption costs. p For water pump power, O t To meet the annual irrigation demand, En c For electricity costs, Q s H0 is the inlet flow rate of the branch pipe, and E is the inlet pressure of the branch pipe. p For pump efficiency, R n The net irrigation quota, EU is the irrigation uniformity coefficient, calculated from the sixth preset formula, which is:
[0161] Among them, CV m The manufacturing deviation of the irrigation device is defined by N, where N is the number of irrigation devices allocated to each crop. Let q be the average flow rate of the emitters in the irrigation triangle. min The minimum flow rate of the emitters in the irrigation triangle can be obtained by solving the nodal flow continuity equation. Energy cost C e It can be calculated based on the power consumption required to raise the irrigation water to a specific working pressure.
[0162] Step 804: Calculate the annual water fee for the irrigation triangle using the seventh preset relation, which is: C w =R g P w , where C w P represents the annual water fee. w For water price, R g The gross irrigation quota for the crop growth period is calculated using the eighth preset formula, which is:
[0163] Step 805: Invest C in the pipeline system a Pipeline maintenance costs C m Energy consumption cost C e And annual water fee C w The sum of these values is taken as the annual irrigation cost CT. That is: CT = C a +C m +C e +C w.
[0164] On this basis, the foregoing step 103 can include:
[0165] Step 806: determining the target setting position of the water outlet stake in each irrigation triangle as the candidate position with the minimum irrigation annual cost.
[0166] In steps 801-806 of the embodiment, a specific method for calculating the economic evaluation parameter of the irrigation triangle is provided. After the economic evaluation parameter corresponding to each candidate position is calculated, the economic differences of different candidate positions can be compared. The candidate position with the minimum economic evaluation parameter (i.e., the position with the best economy) can be taken as the target setting position of the water outlet stake in the irrigation triangle, so as to ensure that the cost of irrigation using the micro-irrigation pipe network is lower and the economy of the micro-irrigation pipe network system is improved.
[0167] In some embodiments, the hydraulic performance evaluation parameter is the deviation rate of the flow rate of the irrigator in the irrigation triangle. Correspondingly, the foregoing step 101 in the embodiment can specifically include steps a1-a6. Figure 9
[0168] Step a1: in the primary triangle, setting an initial position of the water outlet stake, setting a diameter of the branch pipe as a first diameter, the initial position being any position in the primary triangle, and the first diameter being an i-th diameter in the selectable diameters of the branch pipe, the selectable diameters including n selectable diameters sequentially increasing with the serial number, i being a positive integer greater than or equal to 1 and less than or equal to n, and the initial value being 1.
[0169] In this step, the initial position of the water outlet stake and the initial diameter of the branch pipe are set in the primary triangle, which facilitates subsequent calculation. The initial diameter of the branch pipe is the first diameter, and the selectable diameters of the branch pipe can include the diameters of conventional branch pipes on the market, such as 40, 50, 63, 75, and 90 mm. The initial diameter of the branch pipe is the smallest diameter in the selectable diameters, and the branch pipe with a smaller diameter also has a lower cost, thereby ensuring that the trial calculation starts from the branch pipe diameter with the best economy. The initial position of the water outlet stake can also be set according to a preset rule, such as being set as the center position of the primary triangle or a position close to the vertex.
[0170] Step a2: determining a plurality of selectable positions of the water outlet stake in the primary triangle based on the initial position and the position adjustment step.
[0171] After the initial position of the outlet stake is determined, a plurality of optional positions of the outlet stake in the primary triangle can be determined based on the initial position and a position adjustment step. For example, from the initial position, a first direction is determined, which is a direction from the initial position to the boundary of the triangle. Then, a first step is moved from the initial position in the first direction, and a new optional position is obtained. The first step is a multiple of the distance between emitters. The process is repeated until the new optional position reaches the boundary of the triangle or exceeds the boundary of the triangle. Then, the first direction is adjusted by a preset angle adjustment step, and the process is repeated until the adjusted first direction coincides with the initial first direction or the angle between the adjusted first direction and the initial first direction is less than an angle threshold. The angle threshold can be equal to the preset angle adjustment step. When the angle between the adjusted first direction and the initial first direction is less than the angle threshold, it can be determined that all optional first directions have been traversed according to the preset angle, and the process is terminated. In other embodiments, other methods can be used to determine the optional positions of the outlet stake according to actual conditions.
[0172] Step a3: Based on each optional position and the first diameter, the deviation rate of the emitter flow in the primary triangle when the outlet stake is located at each optional position is calculated.
[0173] In this step, for any optional position, the deviation rate of the emitter flow in the primary triangle when the outlet stake is located at the optional position can be calculated. Specifically, the branch pipe, the main pipe and the emitter can be arranged based on the optional position. The arrangement method can refer to the method provided in the previous embodiments. Then, the flow-pressure relationship of the emitter is obtained, which is calibrated in advance or fitted by measuring the flow of the emitter under multiple pressure heads. The flow-pressure relationship of the emitter can be expressed as: q=KH x ,
[0174] wherein q is the flow of the emitter, K is the flow coefficient, H is the pressure head, and x is the flow regime index of the emitter. For the measured fitting, K and x are to-be-fitted parameters. The deviation rate of the emitter flow in the primary triangle when the outlet stake is located at each optional position can be calculated based on the flow-pressure relationship of the emitter.
[0175] For irregular plots, the irregularity of the plot boundary is often accompanied by the unevenness of the plot elevation. In related technologies, it is generally assumed that the plot has a uniform slope, but for irregular plots with obvious unevenness of elevation, this assumption can easily lead to poor pipe network design. In this embodiment, the measured elevation information of the partial elevation measurement points of the irregular plot can also be introduced to optimize the pipe network design.
[0176] Specifically, referring to Figure 10The method of the present application further comprises steps 1001-1002.
[0177] Step 1001: Obtain the elevation information of a plurality of elevation measurement points in the irregular plot.
[0178] Continuing with the irregular plot shown in Figure 3 , the elevation information of a plurality of ground elevation measurement points (referred to as ground elevation measurement points) in the irregular plot can be obtained, which can be specifically represented as an m-row and 4-column matrix:
[0179] Where x i and y i are the coordinate values in the x and y directions, z i is the measured ground relative elevation, m is the number of measured ground elevation measurement points, Figure 3 which is 30 in the example.
[0180] Step 1002: Interpolate the elevation information of all sprinklers in the irregular plot based on the elevation information of the plurality of elevation measurement points.
[0181] After obtaining the elevation information of the plurality of elevation measurement points, the elevation information of each sprinkler can be obtained by interpolation based on the coordinates of each sprinkler in the irregular plot. The interpolation method can include inverse distance weighting, spline function, or Kriging method. Referring to Figure 19 , a diagram showing the distribution of the elevation information of the nodes of the sprinklers in the irregular plot is shown in Figure 3 . In this diagram, the distribution of the elevation of the nodes of the sprinklers is shown in the form of contour lines, and the elevations of the nodes of the sprinklers on each contour line are equal.
[0182] On this basis, referring to Figure 11 , in step a3, the deviation rate of the flow of the sprinklers in the primary triangle is calculated, including steps 1101-1102.
[0183] Step 1101: Calculate the flow of each sprinkler in the primary triangle based on the elevation information of all sprinklers in the irregular plot, the pipe segment energy equation, the pipe network node flow continuity equation, and the sprinkler pressure-flow relationship.
[0184] Specifically, for any pipe segment ij in the primary triangle, the pipe segment energy equation is subject to Bernoulli's law, which is specifically embodied as:
[0185] z i +H i =z j +H j +(1+k l )h f-ij ,
[0186] where z i and z j are the interpolated relative ground elevations at the emitter nodes i and j, respectively, H i and H j are the pressure heads at the emitter nodes i and j, respectively, h f-ij is the head loss along the pipe segment ij, k l is the local head loss coefficient.
[0187] The head loss along the pipe segment ij is calculated using the Hazen-Williams formula:
[0188]
[0189] where L ij is the length of the pipe segment ij, d ij is the diameter of the pipe segment ij, Q ij is the flow rate of the pipe segment ij, and C ij is the roughness coefficient of the pipe segment ij, which is related to the material of the pipe.
[0190] The process of calculating the flow rate of each emitter in the primary triangle includes:
[0191] (1) Set the supply pressure head at the outlet stake to H0;
[0192] (2) Assume that the pressure heads of all emitter nodes in the primary triangle are equal to the pressure head at the outlet stake, and calculate the flow rates of the emitters at various positions in the primary triangle using the flow rate-pressure relationship of the emitters;
[0193] (3) Calculate the flow rates of each pipe segment using the backtracking method;
[0194] Specifically, the flow rates of each pipe segment are obtained by continuously backtracking based on the pipe network node flow continuity equation, using the flow rate of the emitter node at the end of the capillary tube, where the pipe network node flow continuity equation is expressed as: Q i = α∑q l where Q i represents the flow rate of the i-th node, α is the pipe segment along-line flow conversion coefficient, and q l represents the flow rate of the pipe segment connected to the i-th node.
[0195] (4) Calculate the head loss along each pipe segment using the aforementioned head loss calculation formula (Hazen-Williams formula) based on the flow rates of each pipe segment;
[0196] (5) Solve the equation set combining the pipe segment energy equation and the pipe network node flow continuity equation to obtain the new pressure heads of each node;
[0197] (6) Based on the calculated new pressure head of each node, the flow rate of all emitters in the primary triangle is calculated by using the flow rate-pressure relationship of the emitter, i.e., one iteration is completed;
[0198] (7) After each iteration, the accuracy of the iteration result is checked, and when the preset calculation accuracy is met, the iteration is ended, and the flow rate of each emitter node in the primary triangular element under the condition of branch inlet pressure H0 is obtained, wherein, in some embodiments, the preset calculation accuracy is defined by the following formula:
[0199]
[0200] wherein, Q ij-iter(n-1) is the flow rate in pipe segment ij in the (n-1)th iteration result, Q ij-iter(n) is the flow rate in pipe segment ij in the nth iteration result;
[0201] (8) When the accuracy is not met, based on the calculated pressure head of each node, the flow rate of the emitters at each position in the primary triangle is calculated by using the flow rate-pressure relationship of the emitter, and steps (3)-(8) are repeated until the iteration result meets the preset calculation accuracy. Reference Figure 12 , an exemplary shows the flow rate distribution of multiple emitter nodes in the primary triangle ABC when the calculation accuracy meets the preset calculation accuracy.
[0202] Step 1102: Based on the flow rate of each emitter in the primary triangle, the deviation rate of the emitter flow rate in the primary triangle is calculated.
[0203] In this step, based on the calculated flow rate of each emitter, the maximum emitter flow rate and the minimum emitter flow rate in the primary triangle are determined, and the deviation rate of the emitter flow rate is calculated based on the maximum emitter flow rate and the minimum emitter flow rate. Specifically, the deviation rate of the emitter flow rate can be calculated by using the following formula:
[0204]
[0205] wherein, q max is the maximum emitter flow rate in the primary triangle, q min is the minimum emitter flow rate in the primary triangle, q d is the design flow rate of the emitter, and q d can be determined according to the soil properties of the plot (the design flow rate value or range corresponding to different soil properties can be preset, and the design flow rate corresponding to the target soil property can be directly determined by referring to the table later).
[0206] Step a4: determining whether there is at least one candidate position in the plurality of optional positions based on the deviation rate of the emitter flow in the primary triangle and the deviation rate condition, if yes, performing step a5, otherwise performing a6.
[0207] In this step, the deviation rate condition can be that the deviation rate of the emitter flow is less than or equal to a deviation rate threshold, which can be set according to requirements, such as 20%. The deviation rates of the emitter flows corresponding to the plurality of optional positions can be compared with the deviation rate threshold respectively, and the optional positions with deviation rates less than the deviation rate threshold are determined as candidate positions. When it is determined that a certain optional position causes there to be at least one candidate position in the primary triangle, step a5 can be jumped to. When it is determined that all optional positions in the primary triangle cannot cause there to be at least one candidate position in the primary triangle to satisfy the deviation rate condition, step a6 can be jumped to. Referring to Figure 18 , a distribution of the emitter deviation rates corresponding to the plurality of candidate positions in the primary triangle under the first diameter is given (the candidate positions contained in Figure 18 are more, so they are displayed in the form of contour lines, and the numbers on each contour line represent the deviation rate).
[0208] Step a5: determining the primary triangle as an irrigation triangle, determining the first diameter as the target diameter of the branch pipe in the primary triangle, and determining all candidate positions from the plurality of optional positions.
[0209] In this step, when there is at least one candidate position in the primary triangle under the condition that the diameter of the branch pipe is the first diameter, the primary triangle can be directly determined as an irrigation triangle, and the first diameter can be determined as the target diameter of the primary triangle. In order to enable the target setting position of the water stake to be determined from different candidate positions subsequently, all candidate positions in the plurality of optional positions are determined.
[0210] Step a6: adjusting the first diameter of the branch pipe to the i+1th diameter in the optional diameters, and jumping to step a3 until i+1 is greater than n.
[0211] Specifically, step a6 can include: adjusting the first diameter of the branch pipe to the i+1th diameter in the optional diameters (step a61); determining whether i+1 is greater than n (step a62), if not, jumping to step a3, if yes, it can be determined that all optional diameters of the branch pipe cannot ensure that a water stake position satisfying the deviation rate condition is found in the plurality of optional positions, at this time, step a7 can be performed.
[0212] The foregoing steps can actually achieve finding a water stake position satisfying the deviation rate condition by adjusting the water stake position and the diameter of the branch pipe, and optimizing the water stake position based on the hydraulic performance parameter.
[0213] Continuing to refer to Figure 9In step a6, when i+1 is greater than n, the aforementioned step 101 may also include steps a7-a9.
[0214] Step a7: Divide the primary triangle into M secondary triangles of equal area, where M is a positive integer greater than or equal to 2 and has an initial value of 2.
[0215] In implementation, refer to Figure 13 The primary triangle can be divided in the following way: (1) When dividing the primary triangle into two secondary triangles, connect the vertex containing the largest interior angle of the primary triangle with the midpoint of its opposite side to form two secondary triangles with equal areas, such as Figure 13 As shown in a; (2) When dividing the primary triangle into 3 secondary triangles, connect the 3 vertices of the primary triangle to the centroid of the primary triangle respectively to obtain 3 secondary triangles with equal areas, such as Figure 13 As shown in b; (3) When dividing the primary triangle into 4 secondary triangles, first divide the primary triangle into 2 triangles according to (1), and then continue to divide the two generated triangles according to (1) to obtain 4 secondary triangles with equal areas, such as Figure 13 As shown in c; (4) When dividing the primary triangle into 5 triangles, connect the vertex of the primary triangle containing the largest interior angle to the point where its opposite side is 2 / 5 of the side length, to obtain two triangles with areas of 2 / 5 and 3 / 5 of the original primary triangle, respectively. Then, use methods (1) and (2) to divide them into 2 and 3 secondary triangles, respectively, as shown in c. Figure 13 As shown in d. A similar rule can be used to divide the area of a triangle into any equal parts.
[0216] Step a8: Determine whether there is at least one candidate position in each secondary triangle. If there is at least one candidate position in each secondary triangle, proceed to step a9; otherwise, proceed to step a10.
[0217] In this step, a method similar to steps a1-a4 can be used to sequentially determine whether there is at least one candidate position in each secondary triangle. When there is at least one candidate position in each secondary triangle, it is determined that the current number of divisions M can make the divided secondary triangles meet the deviation rate requirement, and step a9 is executed. If there is a secondary triangle that, even after adjusting the branch pipe diameter, cannot find at least one candidate position among multiple selectable diameters, it is determined that the current number of divisions M cannot make the current secondary triangle meet the deviation rate requirement, and step a10 is executed.
[0218] Step a9: Determine each secondary triangle as an irrigation triangle, and sequentially determine all candidate positions among the multiple optional positions of the water outlet pile in each secondary triangle.
[0219] In this step, the plurality of secondary triangles obtained by dividing the primary triangle can be determined as irrigation triangles, and all candidate positions of the plurality of optional positions of the water outlet stake in each secondary triangle are determined, so as to determine the target setting position of the water outlet stake in each secondary triangle subsequently. It should be noted that there can be 0 primary triangles that need to be further divided in the n-2 primary triangles, or all the n-2 primary triangles need to be further divided, which needs to be determined according to the calculation result.
[0220] Step a10: M+1 is assigned to M, and jump to step a7.
[0221] In this step, since the original M value cannot make the secondary triangles obtained by division meet the deviation rate requirement, the M value is increased, and jump to step a7 for reevaluation until at least one candidate position is found in the secondary triangles obtained by division. Referring to Figure 16 , the setting position of the water outlet stake in each secondary triangle, the arrangement mode of the branch pipe and the water flow direction when three first triangles are each divided into two secondary triangles are exemplarily given.
[0222] The above steps can ensure that there is a candidate position of the water outlet stake in any irrigation triangle determined finally, which can meet the deviation rate requirement, and further ensure the uniformity of the pipe network system for crop irrigation.
[0223] In some embodiments, referring to Figure 14 , the method provided by the application further includes steps 1401-1406.
[0224] Step 1401: a plurality of candidate branch pipe connection points are determined in the first triangle, and the first triangle is a primary triangle divided into a plurality of secondary triangles.
[0225] Since one water outlet stake needs to be set in each secondary triangle, the first triangle includes a plurality of water outlet stakes, and therefore a branch pipe connection point needs to be set in the first triangle to supply water for each water outlet stake. In order to determine the most suitable position, a plurality of candidate branch pipe connection points are determined in the first triangle in this step. Exemplarily, a position in the first triangle can be taken as an initial position, and the initial position is adjusted by a certain position adjustment step to obtain a plurality of candidate branch pipe connection points, or the first triangle is directly gridded, and a candidate branch pipe connection point is determined in each grid. In implementation, other ways can also be used to determine the candidate branch pipe connection point according to actual needs.
[0226] Step 1402: a first distance from each candidate branch pipe connection point to any water outlet stake in the first triangle is calculated by a ninth preset relationship, and the ninth preset relationship is:
[0227] wherein i = 1 :M, D Noi represents the first distance from the candidate sub-trunk junction point (x N , y N ) to the i-th water outlet pile, and M represents the total number of water outlet piles in the first triangle.
[0228] Step 1403: calculating the absolute value of the difference between any two first distances according to a tenth preset relationship, which is:
[0229] S ij = |D NOi -D NOj |, wherein i = 1 :M, j = 1 :M, i < j, S ij represents the absolute value of the difference between the i-th first distance and the j-th first distance.
[0230] Step 1404: calculating the first pressure balance evaluation parameter corresponding to each candidate sub-trunk junction point according to an eleventh preset relationship, which is:
[0231] wherein i = 1 :M, j = 1 :M, i < j, S a represents the first pressure balance evaluation parameter.
[0232] Steps 1402-1403 give a specific method for calculating the first pressure balance evaluation parameter, and the first pressure balance parameter calculated by this method can reflect the difference degree of the distance between a candidate sub-trunk junction point and different water outlet piles in the first triangle where the candidate sub-trunk junction point is located.
[0233] Step 1405: determining the candidate sub-trunk junction point with the minimum first pressure balance evaluation parameter as the target sub-trunk junction point.
[0234] In this step, the candidate sub-trunk junction point with the minimum S a value is determined as the target sub-trunk junction point. Referring to Figure 15 , the target sub-trunk junction points N determined in the first triangles respectively containing two, three, four and five secondary triangles are shown.
[0235] Step 1406: setting sub-trunks between each water outlet pile and the target sub-trunk junction point in the first triangle. In this step, after the target sub-trunk junction point is determined, sub-trunks can be set between each water outlet pile and the target sub-trunk junction point in the first triangle, i.e., the sub-trunks adopt a star-shaped arrangement mode.
[0236] In this embodiment, the distances from the target branch pipe connection point to each water outlet stake in the first triangle are relatively equal, which can ensure that the water pressure at each water outlet stake is relatively balanced, thereby improving the uniformity of irrigation.
[0237] In some embodiments, with reference to Figure 17 , the method provided in the application further includes steps 1701-1706.
[0238] Step 1701: determining a plurality of candidate trunk pipe connection points in the irregular land plot.
[0239] In this step, the manner of determining a plurality of candidate trunk pipe connection points in the irregular land plot can refer to the manner of determining candidate branch pipe connection points in the first triangle in step 1401, which is not described herein.
[0240] Step 1702: calculating a plurality of second distances corresponding to each candidate trunk pipe connection point by a twelfth preset relationship, the second distances including distances between the candidate trunk pipe connection point and any branch pipe connection point or distances between the candidate trunk pipe connection point and any water outlet stake in a second triangle, the second triangle being a first triangle that is not divided into secondary triangles, and the twelfth preset relationship being:
[0241] wherein i=1:n-2, D' Noi represents the second distance from the candidate trunk pipe connection point (x' N , y' N ) to the point , the point (x' oi , y' oi ) representing the position of any branch pipe connection point or any water outlet stake in a second triangle, and n representing the number of vertices of the irregular land plot.
[0242] Step 1703: calculating the absolute value of the difference between any two second distances by a thirteenth preset relationship, the thirteenth preset relationship being:
[0243] S' ij =|D' NOi -D' NOj |, wherein i=1:n-2, j=1:n-2, i ij represents the absolute value of the difference between the ith second distance and the jth second distance.
[0244] Step 1704: calculating a second pressure balance evaluation parameter corresponding to each candidate trunk pipe connection point by a fourteenth preset relationship, the fourteenth preset relationship being:
[0245] wherein i=1:n-2, j=1:n-2, i aa second pressure balance evaluation parameter.
[0246] In this step, the second pressure balance evaluation parameter can be equal to the mean value of the absolute value of the difference between any two second distances. It should be noted that the specific type of second distance varies depending on the division case. For the case where each primary triangle is divided into multiple secondary triangles, the category of second distance only includes the distance between the candidate main pipe connection point and any branch pipe connection point; for the case where each primary triangle is not divided into secondary triangles, the category of second distance only includes the distance between the candidate main pipe connection point and any water outlet stake in the second triangle; for the case where part of the primary triangles are divided into secondary triangles, the second distance includes the above two types, and both types of second distance are used to calculate the second pressure balance evaluation parameter.
[0247] Step 1705: determining the candidate main pipe connection point with the smallest second pressure balance evaluation parameter as the target main pipe connection point.
[0248] In this step, the candidate main pipe connection point with the smallest second pressure balance evaluation parameter is determined as the target main pipe connection point. Referring to Figure 20 , point P is the determined target main pipe connection point.
[0249] Step 1706: arranging the main pipe in the irregular plot based on the target main pipe connection point.
[0250] Continuing to refer to Figure 20 , a connection diagram of the target main pipe connection point P and each target branch pipe connection point N i is given in the case where each primary triangle is divided into two secondary triangles. For the case where each primary triangle is not divided into secondary triangles, the target main pipe connection point P is directly connected to each water outlet stake O i ; for the case where part of the primary triangles are divided into secondary triangles, the target main pipe connection point P is connected to the target branch pipe connection point in the first triangle and connected to the water outlet stake in the second triangle. That is, the main pipe also adopts a star arrangement.
[0251] The above steps can ensure that the distance between the main pipe connection point and each target branch pipe connection point or water outlet stake in the second triangle is relatively equal, so that the water pressure distributed in each irrigation triangle is relatively balanced, thereby improving the irrigation uniformity of the pipe network.
[0252] In some embodiments, the method provided by the present application further comprises calculating the diameter of the branch pipe or the main pipe based on the pipe flow. Specifically, the following formula can be used for calculation:
[0253]
[0254] Wherein, D is the minimum pipe diameter of the sub-main pipe or the main pipe, Q is the pipe section flow, V is the economic flow rate, which can be 1.5 m / s. After D is calculated, the pipe diameter that is the least more than the minimum pipe diameter can be selected as the optimized value of the sub-main pipe diameter or the main pipe diameter according to the common specifications in the market.
[0255] In another aspect, referring to Figure 21 The application also provides a micro-irrigation pipe network setting device, which comprises:
[0256] The determining module 210 is configured to divide the irregular land into a plurality of irrigation triangles, determine a plurality of optional positions of the water outlet post in each irrigation triangle, calculate a hydraulic performance evaluation parameter corresponding to each optional position, and determine a plurality of candidate positions from the plurality of optional positions based on the hydraulic performance evaluation parameter, wherein the hydraulic performance evaluation parameter is used to indicate irrigation uniformity when the land is irrigated by using the micro-irrigation pipe network.
[0257] The calculating module 220 is configured to calculate an economic evaluation parameter corresponding to each candidate position of the water outlet post in each irrigation triangle, wherein the economic evaluation parameter is used to indicate the cost required for irrigating the land by using the micro-irrigation pipe network.
[0258] The determining module 210 is further configured to determine a target setting position of the water outlet post in each irrigation triangle from the plurality of candidate positions based on the economic evaluation parameter.
[0259] The setting module 230 is configured to set the water outlet post, the branch pipe and the micro-pipe in each irrigation triangle based on the target setting position, wherein the branch pipe is connected with the water outlet post, the micro-pipe is connected with the branch pipe, and a plurality of water emitters are uniformly arranged on the micro-pipe.
[0260] Optionally, the determining module 210 is further configured to perform steps 201-207, 806, a1-a10, 1001-1002, 1401, 1405, 1701 and 1705 in the foregoing method embodiments; the calculating module 220 is further configured to perform steps 801-805, 1101-1102, 1402-1404 and 1702-1704; and the setting module 230 is further configured to perform steps 501-503, 803, 1406 and 1706.
[0261] The device provided in the embodiments of the application is used to implement the foregoing method. In order to implement each step of the foregoing method, the device of the embodiments of the application can have more modules. The specific operation process of the device can be referred to the method embodiments, which will not be described in detail here.
[0262] In this application, the terms "first" and "second" are used only for descriptive purposes and not to connote or imply relative importance. Other embodiments of the present application will be readily apparent to one skilled in the art upon considering the description hereof in light of the disclosure herein. The description and examples are to be considered exemplary only.
[0263] It is to be understood that the application is not limited to the precise construction herein described and as shown in the drawings, and that various modifications and changes can be made by those skilled in the art without departing from the scope thereof. The scope of the application is limited only by the appended claims.
Claims
1. A method for setting up a micro-irrigation network, characterized in that, The method includes: The irregular plot is divided into multiple irrigation triangles, and multiple optional locations of water piles are determined in each irrigation triangle. The hydraulic performance evaluation parameters corresponding to each optional location are calculated. Based on the hydraulic performance evaluation parameters, multiple candidate locations are determined from the multiple optional locations. The hydraulic performance evaluation parameters are used to indicate the irrigation uniformity when the plot is irrigated using a micro-irrigation network. Calculate the economic evaluation parameter corresponding to each candidate location of the outlet pile in each irrigation triangle, the economic evaluation parameter being used to indicate the cost required to irrigate the plot using the micro-irrigation network; Based on the economic evaluation parameters, the target setting position of the water outlet pile in each irrigation triangle is determined from the plurality of candidate positions; Based on the target location, the water outlet pile, branch pipe and capillary pipe are set in each irrigation triangle. The branch pipe is connected to the water outlet pile, the capillary pipe is connected to the branch pipe, and multiple water emitters are evenly arranged on the capillary pipe. The economic evaluation parameter is the annual irrigation cost. The calculation of the economic evaluation parameter corresponding to each candidate position of the outlet pile in each irrigation triangle includes: The investment in the pipeline system of the irrigation triangle is calculated using a fourth preset relation, which is: ,in i Represents the annual interest rate. m This represents the service life of the branch pipe and the capillary pipe. L l and L m These represent the lengths of the capillary tubes and branch tubes within the irrigation triangle. P l and P m These are the unit prices of the capillary tubes and branch tubes in the irrigation triangle, respectively. S To irrigate the area of the triangle, C a This represents the investment in the pipeline network system; Invest in the pipeline system C a Multiply by the first ratio to obtain the pipeline maintenance cost of the irrigation triangle. C m ; The energy consumption cost of the irrigation triangle is calculated using a fifth preset relation, which is: ,in, C e This represents the energy consumption cost. N p For water pump power, O t The annual irrigation time is designed to meet irrigation needs. E n c For electricity bills, Q s For the branch pipe inlet flow rate, H 0 represents the inlet pressure of the branch pipe. E p For pump efficiency, R n For net irrigation quota, EU The irrigation uniformity coefficient is calculated using the sixth preset formula, which is: ,in, CV m To cause deviations in the water dispenser, N The number of irrigation devices allocated to each crop. The average flow rate of the emitters in the irrigation triangle. q min The minimum flow rate of the emitter in the irrigation triangle; The annual water cost of the irrigation triangle is calculated using the seventh preset relation, which is: ,in, C w This represents the annual water fee. P w For water price, R g The gross irrigation quota for the crop growth period is calculated using the eighth preset formula, which is: ; Invest in the pipeline system C a Pipeline maintenance costs C m The energy consumption cost C e and the annual water fee C w The sum of these values is used as the annual cost of irrigation. The step of determining the target location for the outlet pile in each irrigation triangle from the plurality of candidate locations based on the economic evaluation parameters includes: The target location for the outlet pile is determined by selecting the candidate location with the minimum annual irrigation cost within each irrigation triangle.
2. The method according to claim 1, characterized in that, The method of dividing irregular plots into multiple irrigation triangles includes: Obtain the n vertices of the irregular plot, where n is a positive integer greater than 3; Connect each of the n vertices to the n-2 non-adjacent vertices of each vertex in turn to obtain n partition schemes, each partition scheme including n-2 triangles with a common vertex; The area of each triangle under each partitioning scheme is calculated using a first preset relational formula, which is: ,in, S i Representing the i The area of each triangle, ( x 1 ,y 1), ( x 2 ,y 2) and ( x 3 ,y 3) Representing the first i The coordinates of the three vertices of a triangle; The absolute value of the area difference between any two triangles under each partitioning scheme is calculated using a second preset relational formula, which is: ,in, i =1:(n-2) j =1:(n-2) i < j , S j Indicates the first j The area of each triangle, S i Indicates the first i The area of each triangle, D ij Indicates the first i The first triangle and the second j The difference in the areas of the triangles; The area evaluation parameters for each partitioning scheme are calculated using a third preset relation, which is: ,in, i =1:(n-2) j =1:(n-2) i < j , D a This represents the area evaluation parameter; The n-2 triangles corresponding to the partitioning scheme with the minimum area evaluation parameter are determined as first-order triangles; The primary triangle or the secondary triangle obtained by dividing the primary triangle is used as the irrigation triangle.
3. The method according to claim 2, characterized in that, The hydraulic performance evaluation parameter is the deviation rate of the flow rate of the emitters in the irrigation triangle; The process involves dividing the irregular plot of land into multiple irrigation triangles, determining multiple possible locations for water piles within each triangle, calculating hydraulic performance evaluation parameters for each possible location, and determining multiple candidate locations based on these parameters. This includes: Step a1: In the first-order triangle, set the initial position of the water outlet pile, set the diameter of the branch pipe as the first diameter, the initial position is any position in the first-order triangle, the first diameter is the i-th diameter among the optional diameters of the branch pipe, the optional diameters include n optional diameters that increase sequentially with the sequence number, i is a positive integer greater than or equal to 1 and less than or equal to n, and the initial value is 1; Step a2: Based on the initial position and the position adjustment step size, determine multiple possible positions of the water outlet pile in the first-order triangle; Step a3: Based on each selectable position and the first diameter, calculate the deviation rate of the flow rate of the water emitter in the first-order triangle when the water outlet pile is located at each selectable position; Step a4: Based on the deviation rate of the water flow rate of the irrigation device in the first-order triangle and the deviation rate condition, determine whether there is at least one candidate position among the multiple selectable positions. If so, proceed to step a5; otherwise, proceed to step a6. Step a5: Determine the primary triangle as the irrigation triangle, determine the first diameter as the target diameter of the branch pipe in the primary triangle, and determine all candidate positions from the plurality of selectable positions; Step a6: Adjust the first diameter of the branch pipe to the (i+1)th diameter among the selectable diameters, and jump to step a3 until i+1 is greater than n.
4. The method according to claim 3, characterized in that, In step a6, when i+1 is greater than n, the process of dividing the irregular plot into multiple irrigation triangles, determining multiple possible locations for water piles in each irrigation triangle, calculating hydraulic performance evaluation parameters corresponding to each possible location, and determining multiple candidate locations from the multiple possible locations based on the hydraulic performance evaluation parameters further includes: Step a7: Divide the first-order triangle into M second-order triangles of equal area, where M is a positive integer greater than or equal to 2, and the initial value is 2; Step a8: Sequentially determine whether there is at least one candidate position in each of the secondary triangles. If there is at least one candidate position in each of the secondary triangles, then proceed to step a9; otherwise, proceed to step a10. Step a9: Determine each of the secondary triangles as the irrigation triangle, and sequentially determine all candidate positions among the multiple selectable positions of the water outlet pile in each secondary triangle; Step a10: Assign M+1 to M and jump to step a7.
5. The method according to claim 4, characterized in that, The method further includes: Multiple candidate branch pipe connection points are determined in the first triangle, which is a primary triangle divided into multiple secondary triangles; The first distance from each candidate branch pipe connection point to any outlet pile in the first triangle is calculated using the ninth preset relation. The ninth preset relation is: ,in i =1: M , D Noi Representative of candidate branch pipe connection points To the i One water outlet pile The first distance, M This represents the total number of water outlet piles in the first triangle; The absolute value of the difference between any two first distances is calculated using the tenth preset relation, which is: ,in, i = 1:M , j = 1:M , i < j , S ij Representing the i The first distance and the second j The absolute value of the difference between the first distances; The first pressure balance evaluation parameter corresponding to each candidate branch pipe connection point is calculated using the eleventh preset relation, which is: ,in, i = 1:M , j = 1:M , i < j , S a This represents the first pressure equilibrium evaluation parameter; The candidate branch pipe connection point with the minimum first pressure equalization evaluation parameter is determined as the target branch pipe connection point; A branch pipe is installed between each water outlet pile and the target branch pipe connection point in the first triangle.
6. The method according to claim 5, characterized in that, The method further includes: Multiple candidate trunk connection points were identified in the irregular plot of land; The twelfth preset relation is used to calculate multiple second distances corresponding to each candidate main pipe connection point. These second distances include the distance between the candidate main pipe connection point and any branch pipe connection point, or the distance between the candidate main pipe connection point and any outlet pile in a second triangle. The second triangle is a first-level triangle that has not been divided into second-level triangles. The twelfth preset relation is: ,in i =1: n -2, D' Noi Representative of candidate trunk connection points Time The second distance, point This represents the location of any branch pipe connection point or any outlet pile in the second triangle, and n represents the number of vertices of the irregular plot. The absolute value of the difference between any two second distances is calculated using the thirteenth preset relation, which is: ,in, i =1: n -2, j =1: n -2, i < j , S' ij Representing the i The second distance and the first j The absolute value of the difference between the second distances; The second pressure balance evaluation parameter corresponding to each candidate trunk connection point is calculated using the fourteenth preset relation, which is: ,in, i =1: n -2, j =1: n -2, i < j , S' a This represents the second pressure equilibrium evaluation parameter; The candidate trunk connection point with the smallest second pressure equalization evaluation parameter is determined as the target trunk connection point; Based on the target trunk connection point, trunk pipes are laid out in the irregular plot.
7. The method according to claim 1, characterized in that, The method of setting out the outlet pile, branch pipe and capillary pipe in each of the irrigation triangles includes: The branch pipe is set to extend from the position of the water outlet pile in the irrigation triangle along a first direction toward the boundary of the irrigation triangle. The first direction is perpendicular to the extension direction of the capillary pipe, and the extension direction of the capillary pipe is determined according to the crop planting direction. When the branch pipe extends to intersect with the first side of the irrigation triangle, the branch pipe is set to extend along the second direction until it reaches a vertex of the irrigation triangle. The first side is any side of the irrigation triangle, the second direction is parallel to the first side, and the angle between the second direction and the first direction is an acute angle. The branch pipe ends when it intersects with the first vertex of the irrigation triangle, where the first vertex is any vertex of the irrigation triangle.
8. The method according to claim 3, characterized in that, The method further includes: Obtain elevation information from multiple elevation measurement points in the irregular land parcel; Based on the elevation information of the multiple elevation measurement points, the elevation information of all irrigation devices in the irregular plot is obtained by interpolation. The calculation of the deviation rate of the flow rate of the water emitter in the first-order triangle when the water pile is located at each selectable position includes: Based on the elevation information of all irrigation devices in the irregular plot, the pipe segment energy equation, the pipe network node flow continuity equation, and the irrigation device flow-pressure relationship, the flow rate of each irrigation device in the first-order triangle is calculated. Based on the flow rate of each emitter in the first-order triangle, calculate the deviation rate of the emitter flow rate in the first-order triangle.
9. A micro-irrigation network installation device, characterized in that, The apparatus is configured to perform the method as described in any one of claims 1-8, the apparatus comprising: The determination module is used to divide the irregular plot into multiple irrigation triangles, determine multiple optional locations of water piles in each irrigation triangle, calculate the hydraulic performance evaluation parameters corresponding to each optional location, and determine multiple candidate locations from the multiple optional locations based on the hydraulic performance evaluation parameters. The hydraulic performance evaluation parameters are used to indicate the irrigation uniformity when irrigating the plot using a micro-irrigation network. The calculation module is used to calculate the economic evaluation parameters corresponding to each candidate position of the outlet pile in each irrigation triangle, the economic evaluation parameters being used to indicate the cost required to irrigate the plot using the micro-irrigation network; The determining module is further configured to determine the target setting position of the water outlet pile in each irrigation triangle from the plurality of candidate positions based on the economic evaluation parameters; The setting module is used to set the water outlet pile, branch pipe and capillary pipe in each irrigation triangle based on the target setting position. The branch pipe is connected to the water outlet pile, the capillary pipe is connected to the branch pipe, and multiple water emitters are evenly arranged on the capillary pipe.