A hydraulic calculation method and device for a self-pressure irrigation system

Abstract

The present application relates to the technical field of self-pressure irrigation pipe network hydraulic, disclose a kind of hydraulic calculation method and device of self-pressure irrigation system, the actual water flow and actual working pressure of the multiple water outlets of first branch pipe and second branch pipe are calculated by different formula, and the actual working pressure of the second intersection of second branch pipe and main pipe, and the optimal actual water flow and actual working pressure are determined, finally the present application can ensure that pipe network pressure is uniform, prevent pipeline explosion accident, and in self-pressure irrigation system, it can provide system and comprehensive pipe network hydraulic calculation mode.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of self-pressure irrigation pipe network hydraulics, and particularly relates to a self-pressure irrigation system hydraulic calculation method and device. BACKGROUND

[0002] Drip irrigation is one of the efficient water-saving irrigation technologies in the world today, but the problems of resource waste and energy consumption in the application process of drip irrigation technology have attracted much attention, and environmental protection and energy saving have become the theme of the development of drip irrigation technology.

[0003] There are many large irrigation areas in mountain basin systems in China. Due to the large evaporation and seepage loss of the existing surface water canal system and plain reservoir water, combined with the natural topographic conditions of the mountain basin system, which provides unique water sources and energy conditions for the development of pipeline water distribution and the implementation of self-pressure irrigation, therefore, constructing a large self-pressure pipeline water distribution system from the mountain front to the basin can effectively solve the problem of water resource utilization.

[0004] In the related art, the pressure regulation technology of the current self-pressure irrigation system is not mature, such as the pressure regulation equipment and the pressure regulation technology cannot well guarantee the stable operation of the large-diameter pipeline, in addition, the overall pressure of the pipe network is uneven, combined with the weakening of the aging performance of the pipeline, thereby causing the local point pressure to be too large, which eventually easily causes the main pipe to burst, and in addition, the current research focuses on urban pipe networks, and there is less research on agricultural self-pressure irrigation pipe networks, therefore, in the self-pressure irrigation system, there is a lack of systematic and comprehensive pipe network hydraulic calculation method. SUMMARY

[0005] Therefore, the present application provides a self-pressure irrigation system hydraulic calculation method and device to solve the problem of lack of systematic and comprehensive pipe network hydraulic calculation method in the self-pressure irrigation system.

[0006] In the first aspect, the present application provides a self-pressure irrigation system hydraulic calculation method applied to a self-pressure irrigation system, the self-pressure irrigation system comprising: a water source, a first branch pipe, a second branch pipe and a main pipe with the same pipe diameter, wherein the first branch pipe and the second branch pipe are parallel, the main pipe is perpendicular to the first branch pipe and the second branch pipe, a plurality of water outlets with the same caliber and height are arranged on the first branch pipe and the second branch pipe according to a preset interval, a first electromagnetic valve is arranged at a first intersection of the main pipe and the first branch pipe, and a second electromagnetic valve is arranged at a second intersection of the main pipe and the second branch pipe, water flows from the water source into the main pipe, the first branch pipe and the second branch pipe, and the method comprises:

[0007] Based on the preset range value, the initial working pressure of the end water outlet of the first branch pipe and the second branch pipe is set respectively;

[0008] Based on the branch pipe hydraulic calculation parameters and the initial working pressure, the actual water outlet flow and the actual working pressure of the end water outlet of the first branch pipe and the second branch pipe are calculated respectively through the first hydraulic calculation formula;

[0009] Based on the branch pipe hydraulic calculation parameters, the actual water outlet flow and the actual working pressure of the first water outlet adjacent to the end water outlet of the first branch pipe and the second branch pipe are calculated respectively through the second hydraulic calculation formula;

[0010] The actual water outlet flow and the actual working pressure of multiple water outlets between the first water outlet and the first electromagnetic valve of the first branch pipe and between the first water outlet and the second electromagnetic valve of the second branch pipe are calculated gradually in a recursive manner through the second hydraulic calculation formula;

[0011] Based on the branch pipe hydraulic calculation parameters, the actual working pressure at the second intersection between the second branch pipe and the main pipe is calculated through the third hydraulic calculation formula;

[0012] Based on the branch pipe hydraulic calculation parameters, the target working pressure at the second intersection between the first branch pipe and the main pipe is calculated through the fourth hydraulic calculation formula;

[0013] If the absolute value of the difference between the actual working pressure and the target working pressure is greater than the first preset threshold value, return to the step of setting the initial working pressure of the end water outlet of the first branch pipe and the second branch pipe based on the preset range value;

[0014] If the absolute value of the difference between the actual working pressure and the target working pressure is less than or equal to the first preset threshold value, the total gravitational potential energy of the main pipe is calculated based on the main pipe hydraulic calculation parameters and the branch pipe hydraulic calculation parameters through the fifth calculation formula and the sixth calculation formula;

[0015] If the absolute value of the difference between the actual gravitational potential energy and the total gravitational potential energy between the water source and the first end inlet of the main pipe is less than or equal to the second preset threshold value, the actual water outlet flow and the actual working pressure of the multiple water outlets of the first branch pipe and the second branch pipe are output one by one;

[0016] If the absolute value of the difference between the actual gravitational potential energy and the total gravitational potential energy between the water source and the first end inlet of the main pipe is greater than the second preset threshold value, return again to the step of setting the initial working pressure of the end water outlet of the first branch pipe and the second branch pipe based on the preset range value.

[0017] The embodiments of the present disclosure execute the above-mentioned mode, calculate the actual water flow and the actual working pressure of the multiple water outlets of the first branch pipe and the second branch pipe and the actual working pressure at the second intersection point of the second branch pipe and the main pipe through different formulas, and determine the optimal actual water flow and actual working pressure, so that the present application can finally ensure uniform pipe network pressure, prevent pipe explosion accidents, and provide a comprehensive pipe network hydraulic calculation mode in a self-pressure irrigation system.

[0018] In some optional embodiments, the actual water flow and the actual working pressure of the terminal water outlet of the first branch pipe are calculated through the following first hydraulic calculation formula;

[0019]

[0020] Q u (a,n a )=Q s (a,n a )

[0021]

[0022]

[0023] wherein Q s (a,n a ) is the actual water flow of the terminal water outlet n a of the branch pipe a, μ is the water outlet flow coefficient of the branch pipe a, A s is the water section area of the terminal water outlet n a of the branch pipe a, H s (a,n a ) is the initial working pressure of the terminal water outlet n a of the branch pipe a, Q u (a,n a ) is the upstream branch pipe inter-pipe section flow of the terminal water outlet n a of the branch pipe a, H u (a,n a ) is the upstream section working pressure of the terminal water outlet n a of the branch pipe a, V s (a,n a ) is the outlet average flow velocity of the terminal water outlet n a of the branch pipe a, g is the gravity acceleration, f is the pipe friction coefficient, L s is the vertical pipe length of the branch pipe a, m is the flow index, D S is the vertical pipe inner diameter of the branch pipe a, ζ1 is the variable diameter right-angle bend local resistance coefficient, b is the pipe diameter index, V u (a,n a ) is the terminal water outlet na the average flow rate of the inlet.

[0024] The embodiments of the present disclosure can accurately calculate the actual water flow rate and the actual working pressure of the end water outlet of the first branch pipe and the second branch pipe by the first hydraulic calculation formula through the implementation of the above-mentioned embodiments.

[0025] In some optional embodiments, the actual water flow rate and the actual working pressure of the first water outlet adjacent to the end water outlet of the first branch pipe and the second branch pipe are respectively calculated by the following second hydraulic calculation formula;

[0026]

[0027]

[0028] V u (a,n a )=Q u (a,n a ) / A z

[0029] V s (a,n a-1 )=Q s (a,n a-1 ) / A s

[0030]

[0031]

[0032] Q u (a,n a-1 )=Q s (a,n a-1 )+Q u (a,n a )

[0033] wherein, H d (a,n a-1 ) is the working pressure of the first water outlet n a-1 downstream section of the branch pipe a, H u (a,n a ) is the working pressure of the upstream section of the end water outlet n a of the branch pipe a, f is the pipe friction coefficient, Q u (a,n a ) is the upstream branch pipe inter-pipe segment flow rate of the end water outlet n a of the branch pipe a, L z is the branch pipe pipe segment length between the adjacent two water outlets of the branch pipe a, b is the pipe diameter index, D z is the inner diameter of the branch pipe a, and Zs (a, n a ) is the height difference between the two ends of the branch pipe section between adjacent water outlets, H u (a, n a-1 ) is the first water outlet n a-1 of the branch pipe a, ζ2 is the local resistance coefficient of the straight tee of the branch pipe a, V u (a, n a ) is the inlet average flow velocity of the terminal water outlet n a of the branch pipe a, g is the acceleration of gravity, A z is the water passing section area of the branch pipe a, V s (a, n a-1 ) is the outlet average flow velocity of the first water outlet n a-1 of the branch pipe a, A s is the water passing section area of the terminal water outlet n a of the branch pipe a, Q s (a, n a-1 ) is the actual water discharge flow of the first water outlet n a-1 of the branch pipe a, μ is the water outlet flow coefficient of the branch pipe, H s (a, n a-1 ) is the actual working pressure of the first water outlet n a-1 of the branch pipe a, D S is the inner diameter of the vertical pipe of the branch pipe a, L s is the length of the vertical pipe of the branch pipe a, ζ3 is the local resistance coefficient of the deflection tee of the branch pipe a, Q u (a, n a-1 ) is the upstream branch pipe section flow of the first water outlet n a-1 of the branch pipe a.

[0034] The embodiment of the present disclosure can accurately calculate the actual water discharge flow and the actual working pressure of the first water outlet adjacent to the terminal water outlet of the first branch pipe and the second branch pipe by the second hydraulic calculation formula by executing the above-mentioned implementation manner.

[0035] In some optional implementation manners, the actual water discharge flow and the actual working pressure at the second intersection point of the second branch pipe and the main pipe are calculated by a third hydraulic calculation formula;

[0036]

[0037] Q z (c) = Q u (c, 1)

[0038]

[0039] Wherein, H u(c,1) is the working pressure of the upstream section of the first branch pipe outlet, f is the pipe friction coefficient, Q u (c,1) is the flow rate of the upstream section of the first branch pipe outlet, m is the flow index, L c is the horizontal distance from the first branch pipe outlet to the main pipe, D z is the internal diameter of the branch pipe, Z c0 is the height difference from the first branch pipe outlet to the main pipe, ζ4 is the local resistance coefficient of the first branch pipe valve, V u (c,1) is the average flow velocity of the first branch pipe outlet, g is the acceleration of gravity, Q z (c) is the total flow rate at the second electromagnetic valve, H u (g) is the actual working pressure at the second intersection of the second branch pipe and the main pipe, H z (c) is the total working pressure at the second electromagnetic valve, ζ6 is the local resistance coefficient of the second branch pipe at the three-way valve, V z (c) is the average flow velocity at the second electromagnetic valve.

[0040] The embodiments of the present disclosure can accurately calculate the actual flow rate and the actual working pressure at the second intersection of the second branch pipe and the main pipe by the third hydraulic calculation formula by executing the above-mentioned implementation.

[0041] In some optional implementations, the target flow rate and the target working pressure at the second intersection of the first branch pipe and the main pipe are calculated by the following fourth hydraulic calculation formula;

[0042]

[0043] Q z (a) = Q u (a,1)

[0044]

[0045] Q g (a) = Q z (a)

[0046] Again Substituting, we get:

[0047]

[0048] where H u (c,1) is the working pressure of the upstream section of the first branch pipe outlet, f is the pipe friction coefficient, Q u (c,1) is the flow rate of the upstream section of the first branch pipe outlet, m is the flow index, L a is the horizontal distance from the first branch pipe outlet to the main pipe, Dz is the inner diameter of the branch pipe, Z a0 is the height difference from the first branch pipe outlet to the main pipe, ζ4 is the local resistance coefficient of the first branch pipe valve, V u is the average inlet flow rate of the first branch pipe outlet, g is the acceleration of gravity, Q z is the total outlet flow rate at the first electromagnetic valve, H d is an intermediate parameter, H z is the total working pressure at the first electromagnetic valve, Q g is the total outlet flow rate at the first branch pipe, L g is the distance between the first intersection point and the second intersection point, D g is the inner diameter of the main pipe, b is the pipe diameter index, Z0 is the height difference between the two ends of the main pipe between the first intersection point and the second intersection point, ζ5 is the sum of the local resistance coefficients of the straight-through tees on the main pipe between the first intersection point and the second intersection point, V z is the average flow rate at the first electromagnetic valve, V z is the average flow rate at the second electromagnetic valve, ζ7 is the local resistance coefficient of the straight-through tee at the second branch pipe, H u is the target working pressure at the second intersection point of the first branch pipe and the main pipe.

[0049] The embodiments of the present disclosure can accurately calculate the target outlet flow rate and the target working pressure at the second intersection point of the first branch pipe and the main pipe by the fourth hydraulic calculation formula by executing the above-mentioned embodiments.

[0050] In some optional embodiments, the total gravitational potential energy of the main pipe is calculated by the following fifth calculation formula and sixth calculation formula;

[0051] When there is an intersection point of the second branch pipe and the main pipe inlet, the total gravitational potential energy of the main pipe is calculated by the following fifth calculation formula;

[0052] Q g = Q z (a) + Q z (c)

[0053]

[0054] When there is no intersection point of the second branch pipe and the main pipe inlet, the total gravitational potential energy of the main pipe is calculated by the following sixth calculation formula;

[0055] Q g = Q z (a) + Q z (c)

[0056]

[0057]

[0058] wherein, Q g is the total water flow of the main pipe, Q z (a) is the total water flow at the first electromagnetic valve, Q z (c) is the total water flow at the second electromagnetic valve, Z g is the total gravitational potential energy of the main pipe, H u (g) is the actual working pressure at the second intersection point of the second branch pipe and the main pipe, f is the pipe friction coefficient, m is the flow index, L gs is the length of the pipe section between the water source and the main pipe inlet, D gs is the inner diameter of the pipe between the water source and the main pipe inlet, b is the pipe diameter index, ζ8 is the local resistance coefficient of the main pipe head valve, V g is the average flow velocity of the main pipe, g is the gravitational acceleration, L0 is the length of the pipe section between the second branch pipe and the main pipe inlet, Z1 is the height difference between the second branch pipe and the main pipe inlet, ζ9 is the sum of the local resistance coefficients of the straight-through three-way valves at the main pipe inlet.

[0059] The embodiment of the present disclosure can accurately calculate the total gravitational potential energy of the main pipe by the fifth calculation formula and the sixth calculation formula by executing the above-mentioned implementation.

[0060] In a second aspect, the present application provides a hydraulic calculation device of a self-pressure irrigation system, which is applied to a self-pressure irrigation system, and the self-pressure irrigation system comprises a water source, a first branch pipe, a second branch pipe and a main pipe with the same pipe diameter, wherein the first branch pipe and the second branch pipe are parallel, the main pipe is perpendicular to the first branch pipe and the second branch pipe, a plurality of water outlets with the same caliber and height are arranged on the first branch pipe and the second branch pipe at a preset interval, a first electromagnetic valve is arranged at a first intersection point of the main pipe and the first branch pipe, and a second electromagnetic valve is arranged at a second intersection point of the main pipe and the second branch pipe, water flows from the water source into the main pipe, the first branch pipe and the second branch pipe, and the device comprises:

[0061] A parameter setting module is configured to set initial working pressures of end water outlets of the first branch pipe and the second branch pipe based on preset range values.

[0062] A first calculation module is configured to calculate actual water flow and actual working pressures of the end water outlets of the first branch pipe and the second branch pipe based on branch pipe hydraulic calculation parameters and the initial working pressures by a first hydraulic calculation formula.

[0063] The second calculation module is configured to calculate actual water outlet flow and actual working pressure of the first water outlet adjacent to the end water outlet of the first branch pipe and the second branch pipe based on the branch pipe hydraulic calculation parameters and through a second hydraulic calculation formula.

[0064] The third calculation module is configured to calculate actual water outlet flow and actual working pressure of multiple water outlets between the first water outlet and the first electromagnetic valve of the first branch pipe and between the first water outlet and the second electromagnetic valve of the second branch pipe in a recursive manner through the second hydraulic calculation formula.

[0065] The fourth calculation module is configured to calculate actual working pressure at the second intersection between the second branch pipe and the main pipe based on the branch pipe hydraulic calculation parameters and through a third hydraulic calculation formula.

[0066] The fifth calculation module is configured to calculate target working pressure at the second intersection between the first branch pipe and the main pipe based on the branch pipe hydraulic calculation parameters and through a fourth hydraulic calculation formula.

[0067] The first optimization module is configured to return to a module for setting initial working pressure of the end water outlet of the first branch pipe and the second branch pipe based on preset range values if an absolute value of a difference between the actual working pressure and the target working pressure is greater than a first preset threshold value.

[0068] The sixth calculation module is configured to calculate total gravitational potential energy of the main pipe based on the main pipe hydraulic calculation parameters and the branch pipe hydraulic calculation parameters through a fifth calculation formula and a sixth calculation formula if the absolute value of the difference between the actual working pressure and the target working pressure is less than or equal to the first preset threshold value.

[0069] The parameter output module is configured to output actual water outlet flow and actual working pressure of the multiple water outlets of the first branch pipe and the second branch pipe one by one if an absolute value of a difference between actual gravitational potential energy between the water source and the first end inlet of the main pipe and total gravitational potential energy is less than or equal to a second preset threshold value.

[0070] The second optimization module is configured to return to the module for setting initial working pressure of the end water outlet of the first branch pipe and the second branch pipe based on preset range values again if the absolute value of the difference between the actual gravitational potential energy between the water source and the first end inlet of the main pipe and the total gravitational potential energy is greater than the second preset threshold value.

[0071] In a third aspect, the present application provides a computer device, comprising a memory and a processor, the memory and the processor are communicatively connected with each other, the memory stores computer instructions, and the processor executes the computer instructions to perform the hydraulic calculation method of the self-pressure irrigation system according to the first aspect or any one of the corresponding embodiments.

[0072] In a fourth aspect, the present application provides a computer readable storage medium, having stored thereon computer instructions for causing a computer to execute the hydraulic calculation method of the self-pressure irrigation system according to the first aspect or any one of the corresponding embodiments thereof.

[0073] In a fifth aspect, the present application provides a computer program product, characterized by comprising computer instructions for causing a computer to execute the hydraulic calculation method of the self-pressure irrigation system according to the first aspect or any one of the corresponding embodiments thereof. BRIEF DESCRIPTION OF DRAWINGS

[0074] In order to more clearly illustrate the specific embodiments or prior art technical solutions of the present application, the drawings required to be used in the specific embodiments or prior art description will be briefly introduced as follows. Obviously, the drawings described below are some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor on the basis of these drawings.

[0075] Figure 1 is a schematic diagram of a pipe network of a self-pressure irrigation system according to an embodiment of the present application;

[0076] Figure 2 is a flowchart of a hydraulic calculation method of a self-pressure irrigation system according to an embodiment of the present application;

[0077] Figure 3 is a schematic diagram of a pipe network of a self-pressure irrigation system according to an embodiment of the present application;

[0078] Figure 4 is a schematic diagram of a pipe network of a self-pressure irrigation system according to an embodiment of the present application;

[0079] Figure 5 is a structural block diagram of a hydraulic calculation device of a self-pressure irrigation system according to an embodiment of the present application;

[0080] Figure 6 is a schematic diagram of a hardware structure of a computer device according to an embodiment of the present application. DETAILED DESCRIPTION

[0081] In order to make the purpose, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions of the embodiments of the present application will be described clearly and completely below with reference to the drawings of the embodiments of the present application. Obviously, the described embodiments are some of the embodiments of the present application, but not all the embodiments. Based on the embodiments of the present application, all other embodiments obtained by those skilled in the art without creative labor fall within the scope of protection of the present application.

[0082] According to an embodiment of the present application, a self-pressure irrigation system is provided, as shown inFigure 1 As shown, the self-pressurized irrigation system includes: a water source 10, a first branch pipe 11, a second branch pipe 12, and a main pipe 13, all with the same diameter. The first branch pipe 11 and the second branch pipe 12 are parallel, and the main pipe 13 is perpendicular to both the first branch pipe 11 and the second branch pipe 12. Multiple outlets (n) of the same diameter and height are arranged on the first branch pipe 11 and the second branch pipe 12 at predetermined intervals. a n a-1 ...n a-n A first solenoid valve 14 is installed at the first intersection point D of the main pipe 13 and the first branch pipe 11, and a second solenoid valve 15 is installed at the second intersection point E of the main pipe 13 and the second branch pipe 12. Water flows from the water source 10 into the main pipe 13, the first branch pipe 11 and the second branch pipe 12.

[0083] Specifically, the water source can be a reservoir or a river. Figure 1 In this embodiment, a self-pressurized irrigation system is formed by selecting a first branch pipe 11, a second branch pipe 12, a main pipe 13, a first solenoid valve 14, and a second solenoid valve 15. Multiple outlets of the same diameter and height are arranged on the first branch pipe 11 and the second branch pipe 12 at preset intervals. Using the self-pressurized irrigation system in this embodiment, the hydraulic parameters of the self-pressurized irrigation system are further calculated to prevent uneven pressure distribution in the overall pipe network. Combined with the weakening performance of aging pipes, this can lead to excessive pressure at local points, ultimately causing a main pipe burst accident.

[0084] According to an embodiment of the present invention, this disclosure provides an embodiment of a hydraulic calculation method for a self-pressurized irrigation system. It should be noted that the steps shown in the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions. Furthermore, although a logical order is shown in the flowchart, in some cases, the steps shown or described may be executed in a different order than that shown here.

[0085] This embodiment provides a hydraulic calculation method for a self-pressurized irrigation system, which can be used on mobile terminals such as mobile phones and tablets. Figure 2 This is a flowchart of a hydraulic calculation method for a self-pressurized irrigation system according to an embodiment of the present invention, such as... Figure 2 As shown, the process includes the following steps:

[0086] Step S201: Based on a preset range of values, set the initial working pressure of the end outlets of the first branch pipe and the second branch pipe respectively.

[0087] Specifically, the preset range value is n at the end outlet of the first branch pipe and the second branch pipe. a Let H be the reasonable range of initial working pressure (h1, h2). s (a,n a) = (h1 + h2) / 2, H s is the initial working pressure of the end water outlet n a of the first branch pipe and the second branch pipe, and the preset range value can be flexibly set in combination with the actual application scene.

[0088] Step S202, based on the branch pipe hydraulic calculation parameters and the initial working pressure, the actual water outlet flow and the actual working pressure of the end water outlet of the first branch pipe and the second branch pipe are calculated respectively through the first hydraulic calculation formula.

[0089] Specifically, the branch pipe hydraulic calculation parameters are hydraulic calculation parameters related to the first branch pipe or the second branch pipe, such as Figure 3 As shown, the actual water outlet flow Q a (a, n s ) and the actual working pressure H a (a, n s ) of the end water outlet n a of the branch pipe a, where the branch pipe a includes the first branch pipe and the second branch pipe.

[0090] In a specific example, the first hydraulic calculation formula is expressed by the following formula:

[0091]

[0092] Q u (a, n a ) = Q s (a, n a )

[0093]

[0094]

[0095] Wherein, Q s (a, n a ) is the actual water outlet flow of the end water outlet n a of the branch pipe a, μ is the water outlet flow coefficient of the branch pipe a, A s is the water section area of the end water outlet n a of the branch pipe a, H s (a, n a ) is the initial working pressure of the end water outlet n a of the branch pipe a, Q u (a, n a ) is the upstream branch pipe inter-pipe segment flow of the end water outlet n a of the branch pipe a, H u (a, n a ) is the upstream section working pressure of the end water outlet n a of the branch pipe a, and V s(a, n a ) is the average flow velocity of the outlet of the terminal water outlet n a of the branch pipe a, g is the acceleration of gravity, f is the pipe friction coefficient, L s is the vertical pipe length of the branch pipe a, m is the flow index, D S is the vertical pipe inner diameter of the branch pipe a, ζ1 is the local resistance coefficient of the variable diameter right-angle bend, b is the pipe diameter index, V u (a, n a ) is the average flow velocity of the inlet of the terminal water outlet n a of the branch pipe a.

[0096] In the above formula, except for the actual water discharge Q s (a, n a ) and the actual working pressure H s (a, n a ), the rest of the parameters belong to the branch pipe hydraulic calculation parameters.

[0097] Step S203, based on the branch pipe hydraulic calculation parameters, the actual water discharge and the actual working pressure of the first water outlet adjacent to the terminal water outlet of the first branch pipe and the second branch pipe are respectively calculated through the second hydraulic calculation formula.

[0098] Specifically, in Figure 3 , the first water outlet adjacent to the terminal water outlet a is n a-1 , the actual water discharge Q a-1 (a, n s ) and the actual working pressure H s (a, n a-1 ) of the first water outlet n a-1 .

[0099] In a specific example, the second hydraulic calculation formula is expressed by the following formula:

[0100]

[0101]

[0102] V u (a, n a ) = Q u (a, n a ) / A z

[0103] V s (a, n a-1 ) = Q s (a, n a-1 ) / A s

[0104]

[0105]

[0106] Q u (a,n a-1 )=Q s (a,n a-1 )+Q u (a,n a )

[0107] where H d (a,n a-1 ) is the working pressure of the downstream section of the first water outlet n a-1 of the branch pipe a, H u (a,n a ) is the working pressure of the upstream section of the last water outlet n a of the branch pipe a, f is the pipe friction coefficient, Q u (a,n a ) is the upstream branch pipe inter-pipe segment flow of the last water outlet n a of the branch pipe a, L z is the length of the branch pipe segment between the adjacent two water outlets of the branch pipe a, b is the pipe diameter index, D z is the inner diameter of the branch pipe a, Z s (a,n a ) is the height difference between the two ends of the branch pipe segment between the adjacent water outlets, H u (a,n a-1 ) is the working pressure of the upstream section of the first water outlet n a-1 of the branch pipe a, ζ2 is the local resistance coefficient of the straight tee of the branch pipe a, V u (a,n a ) is the inlet average flow velocity of the last water outlet n a of the branch pipe a, g is the gravitational acceleration, A z is the water passing section area of the branch pipe a, V s (a,n a-1 ) is the outlet average flow velocity of the first water outlet n a-1 of the branch pipe a, A s is the water passing section area of the last water outlet n a of the branch pipe a, Q s (a,n a-1 ) is the actual water flow of the first water outlet n a-1 of the branch pipe a, μ is the water outlet flow coefficient of the branch pipe, H s (a,n a-1 ) is the actual working pressure of the first water outlet n a-1 of the branch pipe a, D S is the inner diameter of the vertical pipe of the branch pipe a, L s is the length of the vertical pipe of the branch pipe a, ζ3 is the local resistance coefficient of the baffle tee of the branch pipe a, Qu (a, n a-1 ) is the first outlet n a-1 of the branch pipe a

[0108] In the above formula, in addition to the actual water flow Q a-1 of the first outlet n s (a, n a-1 ) and the actual working pressure H s (a, n a-1 ), the remaining parameters all belong to the branch pipe hydraulic calculation parameters.

[0109] Step S204, the actual water flow and the actual working pressure of the multiple outlets between the first outlet of the first branch pipe and the first electromagnetic valve and between the first outlet of the second branch pipe and the second electromagnetic valve are gradually calculated in a recursive manner through the second hydraulic calculation formula.

[0110] Specifically, in Figure 1 , since there are multiple outlets between the first outlet of the first branch pipe and the first electromagnetic valve and between the first outlet of the second branch pipe and the second electromagnetic valve, the actual water flow and the actual working pressure of the previous outlet are taken as the input parameters of the subsequent outlet, in Figure 1 , the terminal outlet is the subsequent outlet n a , and the previous outlet is n a-1 , the actual water flow and the actual working pressure of the multiple outlets between the first outlet of the first branch pipe and the first electromagnetic valve and between the first outlet of the second branch pipe and the second electromagnetic valve are gradually calculated through the above-mentioned second formula.

[0111] Exemplarily, according to the calculation results of the first hydraulic calculation parameters of the first branch pipe and the terminal outlet, the calculation is carried out in a recursive manner, the second hydraulic calculation formula is taken as a calling formula, the second hydraulic calculation formula is repeatedly called by using the for loop in the python environment, and the variable is initialized each time the second hydraulic calculation formula is called, that is, the input value of the subsequent outlet is the calculation value of the previous outlet, the loop is called until the first end outlet of the first branch pipe is calculated, and finally H and Q of all outlets of the first branch pipe except the terminal outlet can be obtained, and the calculation results are saved.

[0112] Step S205, based on the branch pipe hydraulic calculation parameters, the actual working pressure at the second intersection between the second branch pipe and the main pipe is calculated through the third hydraulic calculation formula.

[0113] Specifically, in Figure 1 , the actual working pressure at the second intersection E between the second branch pipe 12 and the main pipe 13 is H u (g). As Figure 4 ​​​​​​​​As shown, the first branch pipe first end water outlet and the second branch pipe first end water outlet are shown, and Figure 4 wherein the a branch pipe is the first branch pipe, the c branch pipe is the second branch pipe, L gs is the length of the pipe section of the second branch pipe to the water source, L a is the horizontal distance from the first branch pipe first end water outlet to the main pipe, L g is the distance between the first intersection point D and the second intersection point E, L c is the distance between the first intersection point E and the second branch pipe first end water outlet.

[0114] In a specific example, the actual working pressure at the second intersection point of the second branch pipe and the main pipe is calculated by the following third hydraulic calculation formula:

[0115]

[0116] Q z (c) = Q u (c,1)

[0117]

[0118] wherein H u (c,1) is the working pressure of the upstream section of the second branch pipe first end water outlet, f is the pipe friction coefficient, Q u (c,1) is the flow rate of the upstream section of the second branch pipe first end water outlet, m is the flow index, L c is the horizontal distance from the second branch pipe first end water outlet to the main pipe, D z is the internal diameter of the branch pipe, Z c0 is the height difference from the second branch pipe first end water outlet to the main pipe, ζ4 is the local resistance coefficient of the branch pipe first end valve, V u (c,1) is the average flow velocity of the inlet of the second branch pipe first end water outlet, g is the gravitational acceleration, Q z (c) is the total water outlet flow rate at the second electromagnetic valve, H u (g) is the actual working pressure at the second intersection point of the second branch pipe and the main pipe, H z (c) is the total working pressure at the second electromagnetic valve, ζ6 is the local resistance coefficient of the branch pipe at the second electromagnetic valve, V z (c) is the average flow velocity at the second electromagnetic valve.

[0119] In the above formula, in addition to the actual working pressure H u (g) at the second intersection point of the second branch pipe and the main pipe, the remaining parameters are all branch pipe hydraulic calculation parameters.

[0120] Step S206, based on the branch pipe hydraulic calculation parameters, the target working pressure at the second intersection point of the first branch pipe and the main pipe is calculated through the fourth hydraulic calculation formula.

[0121] Specifically, in Figure 1 , the actual working pressure at the second intersection point E of the second branch pipe 12 and the main pipe 13 is H u ′(g). In Figure 4 , the first end water outlet of the first branch pipe and the first end water outlet of the second branch pipe are shown, L gs is the pipe segment length of the second branch pipe to the water source, L a is the horizontal distance from the first branch pipe first end water outlet to the main pipe, L g is the distance between the first intersection point D and the second intersection point E, L c is the distance between the first intersection point E and the first end water outlet of the second branch pipe c.

[0122] In a specific example, the target water outlet flow rate and the target working pressure at the second intersection point of the first branch pipe and the main pipe are calculated through the following fourth hydraulic calculation formula;

[0123]

[0124] Q z (a) = Q u (a,1)

[0125]

[0126] Q g (a) = Q z (a)

[0127] Again Substituting, we get:

[0128]

[0129] Where H u (a,1) is the working pressure of the first end water outlet of the first branch pipe upstream section, f is the pipe friction coefficient, Q u (a,1) is the water flow rate of the first end water outlet of the first branch pipe upstream section, m is the flow index, L a is the horizontal distance from the first end water outlet of the first branch pipe to the main pipe, D z is the internal diameter of the branch pipe, Z a0 is the elevation difference from the first end water outlet of the first branch pipe to the main pipe, ζ4 is the local resistance coefficient of the first end valve of the branch pipe, V u (a,1) is the inlet average flow rate of the first end water outlet of the first branch pipe, g is the gravitational acceleration, Q z (a) is the total water outlet flow rate at the first electromagnetic valve, H d(g) is an intermediate parameter, H z (a) is the total working pressure at the first electromagnetic valve, Q g (a) is the total water outlet flow at the first branch pipe, L g is the distance between the first intersection point and the second intersection point, D g is the internal diameter of the main pipe, b is the pipe diameter index, Z0 is the height difference between the two ends of the main pipe between the first intersection point and the second intersection point, ζ5 is the sum of the local resistance coefficients of each straight tee on the main pipe between the first intersection point and the second intersection point, V z (a) is the average flow rate at the first electromagnetic valve, V z (c) is the average flow rate at the second electromagnetic valve, ζ7 is the local resistance coefficient of the straight tee at the second branch pipe, H u ′(g) is the target working pressure at the second intersection point of the first branch pipe and the main pipe.

[0130] In the above formula, in addition to the target working pressure H u ′(g), the remaining parameters are branch pipe hydraulic calculation parameters.

[0131] Step S207, if the absolute value of the difference between the actual working pressure and the target working pressure is greater than the first preset threshold value, return to the step of setting the initial working pressure of the end water outlet of the first branch pipe based on the preset range value.

[0132] Specifically, the actual working pressure is H u (g), the target working pressure H u ′(g), the first preset threshold value is E1, and further judgment is made on |H u (g)-H u ′(g)|≦E1=0.01, and the first preset threshold value E1 can be flexibly set in combination with the actual scene. If |H u (g)-H u ′(g)|>E1=0.01, return to the step S201 of setting the initial working pressure of the end water outlet of the first branch pipe based on the preset range value, if |H u (g)-H u ′(g)|≦E1=0.01, then proceed to the next calculation, otherwise, return to step S201 using the optimization principle and cycle calculation.

[0133] Step S208, if the absolute value of the difference between the actual working pressure and the target working pressure is less than or equal to the first preset threshold value, the total gravitational potential energy of the main pipe is calculated based on the main pipe hydraulic calculation parameters and the branch pipe hydraulic calculation parameters through the fifth calculation formula and the sixth calculation formula.

[0134] Specifically, if |H u (g)-Hu (g) | ≦ E1 = 0.01, total gravitational potential energy Z of the main pipe is calculated g .

[0135] In a specific example, the total gravitational potential energy of the main pipe is calculated by the following fifth calculation formula and sixth calculation formula.

[0136] When there is an intersection point of the second branch pipe and the main pipe water inlet, the total gravitational potential energy of the main pipe is calculated by the following fifth calculation formula;

[0137] Q g = Q z (a) + Q z (c)

[0138]

[0139] When there is no intersection point of the second branch pipe and the main pipe water inlet, the total gravitational potential energy of the main pipe is calculated by the following sixth calculation formula;

[0140] Q g = Q z (a) + Q z (c)

[0141]

[0142]

[0143] wherein Q g is the total water outlet flow of the main pipe, Q z (a) is the total water outlet flow at the first electromagnetic valve, Q z (c) is the total water outlet flow at the second electromagnetic valve, Z g ' is the total gravitational potential energy of the main pipe, H u (g) is the actual working pressure at the second intersection point of the second branch pipe and the main pipe, f is the pipe friction coefficient, m is the flow index, L gs is the pipe segment length between the water source and the main pipe water inlet, D gs is the pipe inner diameter between the water source and the main pipe water inlet, b is the pipe diameter index, ζ8 is the local resistance coefficient of the main pipe head valve, V g is the average flow rate of the main pipe, g is the gravitational acceleration, L0 is the pipe segment length between the second branch pipe and the main pipe water inlet, Z1 is the elevation difference between the second branch pipe and the main pipe water inlet, ζ9 is the sum of the local resistance coefficients of the straight-through three-way valves of the main pipe water inlet.

[0144] In the above parameters, the parameter associated with the first branch pipe is a first hydraulic calculation parameter of the first branch pipe, the parameter associated with the second branch pipe is a second hydraulic calculation parameter of the second branch pipe, and the parameter associated with the main pipe is a third hydraulic calculation parameter of the main pipe.

[0145] In step S209, if the absolute value of the difference between the actual gravity potential energy and the total gravity potential energy between the water source and the first end inlet of the main pipe is less than or equal to a second preset threshold value, the actual water outlet flow rate and the actual working pressure of each of the plurality of water outlets of the first branch pipe and the second branch pipe are output.

[0146] In step S210, if the absolute value of the difference between the actual gravity potential energy and the total gravity potential energy between the water source and the first end inlet of the main pipe is greater than the second preset threshold value, the step of setting the initial working pressure of the last water outlet of the first branch pipe based on the preset range value is returned again.

[0147] Specifically, the actual gravity potential energy is Z g , the total gravity potential energy is Z g , and the second preset threshold value is E2, which can be flexibly applied in combination with actual application scenarios.

[0148] Specifically, the main pipe section connected between the reservoir and the irrigation area is taken as a carrier, and the fifth calculation formula and the sixth calculation formula are used to calculate the total output potential energy Z g ′ of the main pipe section. The self-pressure irrigation technology is to convert the gravity potential energy accumulated by the pipe into the working pressure water head required for irrigation by using the natural terrain elevation difference, and by using this feature, the calculated total gravity potential energy Z g ′ of the main pipe section is compared with the accumulated actual gravity potential energy Z g to determine whether the difference between them meets the accuracy requirement of |Z g -Z g ′|≦E2=0.01. If it meets the requirement, the actual working pressure and the actual water outlet flow rate values of all the water outlets of the first branch pipe and the second branch pipe are output one by one. If it does not meet the requirement, i.e., |Z g -Z g ′|﹥E2=0.01, the pressure value range is narrowed by using the optimization principle, and the appropriate initial working pressure is selected for cyclic calculation until the accuracy requirement is met and the result is output.

[0149] In the present embodiment, a hydraulic calculation device for a self-pressure irrigation system is also provided, which is used to implement the above embodiments and preferred embodiments, and will not be described again. As used below, the term "module" can be a combination of software and / or hardware that implements a predetermined function. Although the device described in the following embodiments is preferably implemented in software, hardware, or a combination of software and hardware is also possible and contemplated.

[0150] The embodiment provides a hydraulic calculation device of a self-pressure irrigation system, which comprises the following modules. Figure 5

[0151] The parameter setting module 51 is configured to set initial working pressures of end water outlets of the first branch pipe and the second branch pipe respectively based on preset range values.

[0152] The first calculation module 52 is configured to calculate actual water flow rates and actual working pressures of the end water outlets of the first branch pipe and the second branch pipe respectively by a first hydraulic calculation formula based on branch pipe hydraulic calculation parameters and the initial working pressures.

[0153] The second calculation module 53 is configured to calculate actual water flow rates and actual working pressures of first water outlets adjacent to the end water outlets of the first branch pipe and the second branch pipe respectively by a second hydraulic calculation formula based on the branch pipe hydraulic calculation parameters.

[0154] The third calculation module 54 is configured to calculate actual water flow rates and actual working pressures of a plurality of water outlets between the first water outlets and the first electromagnetic valve of the first branch pipe and between the first water outlets and the second electromagnetic valve of the second branch pipe gradually by the second hydraulic calculation formula in a recursive manner.

[0155] The fourth calculation module 55 is configured to calculate an actual working pressure at a second intersection between the second branch pipe and the main pipe by a third hydraulic calculation formula based on the branch pipe hydraulic calculation parameters.

[0156] The fifth calculation module 56 is configured to calculate a target working pressure at the second intersection between the first branch pipe and the main pipe by a fourth hydraulic calculation formula based on the branch pipe hydraulic calculation parameters.

[0157] The first optimization module 57 is configured to return to the module of setting the initial working pressures of the end water outlets of the first branch pipe and the second branch pipe respectively based on the preset range values if an absolute value of a difference between the actual working pressure and the target working pressure is greater than a first preset threshold value.

[0158] The sixth calculation module 58 is configured to calculate total gravitational potential energy of the main pipe by a fifth calculation formula and a sixth calculation formula based on main pipe hydraulic calculation parameters and the branch pipe hydraulic calculation parameters if the absolute value of the difference between the actual working pressure and the target working pressure is less than or equal to the first preset threshold value.

[0159] The parameter output module 59 is configured to output the actual water flow rates and the actual working pressures of the plurality of water outlets of the first branch pipe and the second branch pipe one by one if an absolute value of a difference between actual gravitational potential energy between the water source and the first end inlet of the main pipe and the total gravitational potential energy is less than or equal to a second preset threshold value.

[0160] ​The second optimization module 60 is configured to return to the preset range-based value setting module if the absolute value of the difference between the actual gravity potential energy and the total gravity potential energy between the water source and the first end inlet of the main pipe is greater than the second preset threshold value.

[0161] Further function descriptions of the above-mentioned modules and units are the same as those of the corresponding embodiments, and will not be described here.

[0162] The hydraulic calculation device of the self-pressure irrigation system in the embodiment is in the form of a functional unit, where the unit refers to an ASIC (Application Specific Integrated Circuit) circuit, a processor and a memory executing one or more software or fixed programs, and / or other devices capable of providing the above functions.

[0163] The embodiment of the present application also provides a computer device having the hydraulic calculation device of the self-pressure irrigation system.

[0164] Please refer to Figure 6 , Figure 6 is a structural schematic diagram of a computer device provided by an optional embodiment of the present application, as shown in Figure 6 , the computer device comprises one or more processors 10, a memory 20, and an interface for connecting various components, including a high-speed interface and a low-speed interface. Various components are communicatively connected with each other by different buses, and can be installed on a common mainboard or in other ways as needed. The processor can process instructions executed in the computer device, including instructions stored in the memory or graphics information of the memory to display a GUI on an external input / output device (such as a display device coupled to the interface). In some optional embodiments, multiple processors and / or multiple buses can be used with multiple memories and multiple memories as needed. Similarly, multiple computer devices can be connected, and each device provides part of the necessary operations (for example, as a server array, a group of blade servers, or a multi-processor system). Figure 6 In the embodiment, the processor 10 is taken as an example.

[0165] The processor 10 can be a central processor, a network processor or a combination thereof. The processor 10 can further include a hardware chip. The hardware chip can be an application specific integrated circuit, a programmable logic device or a combination thereof. The programmable logic device can be a complex programmable logic device, a field programmable logic gate array, a general array logic or any combination thereof.

[0166] The memory 20 stores instructions executable by the at least one processor 10 to cause the at least one processor 10 to perform the methods illustrated by the above embodiments.

[0167] The memory 20 can include a program storage area and a data storage area. The program storage area can store an operating system, application programs required by at least one function, and the like. The data storage area can store data created according to the use of the computer device, and the like. In addition, the memory 20 can include a high-speed random access memory, and can further include a non-transitory memory such as at least one magnetic disk storage device, a flash memory device, or other non-transitory solid-state memory device. In some alternative embodiments, the memory 20 can optionally include a memory disposed remotely from the processor 10, which can be connected to the computer device through a network. Examples of the network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and a combination thereof.

[0168] The memory 20 can include a volatile memory such as a random access memory, and can further include a non-volatile memory such as a flash memory, a hard disk, or a solid state disk, and a combination of the above-mentioned kinds of memories.

[0169] The computer device further includes a communication interface 30 for communication of the computer device with other devices or communication networks.

[0170] The embodiments of the present application also provide a computer readable storage medium, and the above-mentioned method according to the embodiments of the present application can be implemented in hardware, firmware, or recorded in a storage medium, or stored in a remote storage medium or a non-transitory machine readable storage medium and downloaded to a local storage medium through network downloading of computer code, so that the method described herein can be processed by such software on a storage medium using a general purpose computer, a special purpose processor, or programmable or special purpose hardware. The storage medium can be a magnetic disk, an optical disk, a read-only memory, a random access memory, a flash memory, a hard disk, or a solid state disk, and the like. Further, the storage medium can also include a combination of the above-mentioned kinds of memories. It can be understood that the computer, the processor, the microprocessor controller, or the programmable hardware includes a storage component that can store or receive software or computer code, which, when accessed and executed by the computer, the processor, or the hardware, implements the methods illustrated by the above embodiments.

[0171] Part of the present application can be applied as a computer program product, for example, computer program instructions, when executed by a computer, through the operation of the computer, can invoke or provide the method and / or technical solutions according to the present application. Those skilled in the art should understand that the form of computer program instructions in computer readable medium includes but is not limited to source file, executable file, installation package file and the like, and accordingly, the way of computer program instructions executed by computer includes but is not limited to: the computer directly executes the instructions, or the computer compiles the instructions and then executes the corresponding compiled program, or the computer reads and executes the instructions, or the computer reads and installs the instructions and then executes the corresponding installed program. Here, the computer readable medium can be any available computer readable storage medium or communication medium accessible to the computer.

[0172] Although the embodiments of the present application are described in conjunction with the drawings, various modifications and changes can be made by those skilled in the art without departing from the spirit and scope of the present application, and such modifications and changes fall within the scope defined by the appended claims.

Claims

1. A hydraulic calculation method for a self-pressurized irrigation system, characterized in that, The method is applied to the self-pressurized irrigation system, which includes: a water source, a first branch pipe, a second branch pipe, and a main pipe of the same diameter. The first and second branch pipes are parallel, and the main pipe is perpendicular to both the first and second branch pipes. Multiple outlets of the same diameter and height are arranged on the first and second branch pipes at predetermined intervals. A first solenoid valve is installed at a first intersection of the main pipe and the first branch pipe, and a second solenoid valve is installed at a second intersection of the main pipe and the second branch pipe. Water flows from the water source into the main pipe, the first branch pipe, and the second branch pipe. Based on a preset range of values, the initial working pressure of the end outlets of the first branch pipe and the second branch pipe is set respectively. Based on the hydraulic calculation parameters of the branch pipe and the initial working pressure, the actual water flow rate and actual working pressure of the end outlet of the first branch pipe and the second branch pipe are calculated respectively using the first hydraulic calculation formula. Based on the hydraulic calculation parameters of the branch pipe, the actual water flow rate and actual working pressure of the first outlet adjacent to the end outlet of the first branch pipe and the second branch pipe are calculated respectively using the second hydraulic calculation formula. Using the second hydraulic calculation formula, the actual water flow and actual working pressure of multiple outlets between the first outlet of the first branch pipe and the first solenoid valve, and between the first outlet of the second branch pipe and the second solenoid valve, are calculated step by step in a recursive manner. Based on the hydraulic calculation parameters of the branch pipe, the actual working pressure at the second intersection of the second branch pipe and the main pipe is calculated using the third hydraulic calculation formula. Based on the hydraulic calculation parameters of the branch pipe, the target working pressure at the second intersection of the first branch pipe and the main pipe is calculated using the fourth hydraulic calculation formula. If the absolute value of the difference between the actual working pressure and the target working pressure is greater than the first preset threshold, return to the step of setting the initial working pressure of the end outlet of the first branch pipe and the second branch pipe respectively based on the preset range value; If the absolute value of the difference between the actual working pressure and the target working pressure is less than or equal to the first preset threshold, the total gravitational potential energy of the main pipe is calculated based on the hydraulic calculation parameters of the main pipe and the hydraulic calculation parameters of the branch pipe, using the fifth calculation formula and the sixth calculation formula. If the absolute value of the difference between the actual gravitational potential energy between the water source and the first inlet of the main pipe and the total gravitational potential energy is less than or equal to the second preset threshold, the actual water flow rate and actual working pressure of multiple outlets of the first branch pipe and the second branch pipe are output one by one. If the absolute value of the difference between the actual gravitational potential energy between the water source and the first inlet of the main pipe and the total gravitational potential energy is greater than the second preset threshold, the process returns to the step of setting the initial working pressure of the end outlets of the first branch pipe and the second branch pipe respectively based on the preset range value.

2. The method according to claim 1, characterized in that, The actual flow rate and actual working pressure of the end outlets of the first branch pipe and the second branch pipe are calculated using the following first hydraulic calculation formula. Q u (a,n a )=Q s (a,n a ) Among them, Q s (a,n a (n) is the outlet at the end of branch pipe a. a The actual outflow rate, μ is the outflow coefficient of branch pipe a, A s The outlet n at the end of branch pipe a a The cross-sectional area of ​​the water passage, H s (a,n a (n) is the outlet at the end of branch pipe a. a Initial working pressure, Q u (a,n a (n) is the outlet at the end of branch pipe a. a The flow rate of the upstream branch pipe section, H u (a,n a (n) is the outlet at the end of branch pipe a. a The upstream section working pressure, V s (a,n a The outlet n at the end of branch pipe a a The average outlet velocity, g is the acceleration due to gravity, f is the pipe friction coefficient, and L s Let m be the length of the vertical pipe of branch a, and m be the flow rate index. S Let V be the inner diameter of the vertical pipe of branch a, ζ1 be the local resistance coefficient of the reducing right-angle elbow, b be the pipe diameter index, and V be the vertical pipe diameter index. u (a,n a (n) is the outlet at the end of branch pipe a. a The average inflow velocity.

3. The method according to claim 1, characterized in that, The actual flow rate and actual working pressure of the first outlet adjacent to the end outlet of the first branch pipe and the second branch pipe are calculated respectively using the following second hydraulic calculation formula. V u (a,n a )=Q u (a,n a ) / A z V s (a,n a-1 )=Q s (a,n a-1 ) / A s Q u (a,n a-1 )=Q s (a,n a-1 )+Q u (a,n a ) Among them, H d (a,n a-1 (n) is the first outlet of branch pipe a. a-1 Downstream section working pressure, H u (a,n a (n) is the outlet at the end of branch pipe a. a The upstream section working pressure, f is the pipe friction coefficient, Q u (a,n a (n) is the outlet at the end of branch pipe a. a The flow rate of the upstream branch pipe section, L z Let D be the length of the branch pipe segment between two adjacent outlets of branch pipe a, b be the pipe diameter index, and D be the diameter index. z Z is the inner diameter of branch pipe a. s (a,n a H represents the elevation difference between the two ends of the branch pipe section between adjacent outlets. u (a,n a-1 (n) is the first outlet of branch pipe a. a-1 The upstream section working pressure, ζ2 is the local resistance coefficient of the DC tee of branch pipe a, V u (a,n a (n) is the outlet at the end of branch pipe a. a The average inlet velocity, g is the acceleration due to gravity, A z V is the cross-sectional area of ​​branch pipe a. s (a,n a-1 (n) is the first outlet of branch pipe a. a-1 The average outlet velocity, A s The outlet n at the end of branch pipe a a The cross-sectional area of ​​the water passage, Q s (a,n a-1 (n) is the first outlet of branch pipe a. a-1 The actual outflow rate, μ is the outflow coefficient of the branch pipe, H s (a,n a-1 (n) is the first outlet of branch pipe a. a-1 The actual work pressure, D S L is the inner diameter of the vertical pipe of branch a. s Let ζ3 be the vertical pipe length of branch a, and ζ3 be the local resistance coefficient of the deflector tee of branch a. u (a,n a-1 (n) is the first outlet of branch pipe a. a-1 The flow rate of the upstream branch pipe section.

4. The method according to claim 1, characterized in that, The actual outflow rate and actual working pressure at the second intersection of the second branch pipe and the main pipe are calculated using the following third hydraulic calculation formula. Q z (c)=Q u (c,1) Among them, H u (c,1) represents the working pressure at the upstream section of the outlet of the second branch pipe, f is the pipe friction coefficient, and Q is the pressure at the upstream section of the outlet of the second branch pipe. u (c,1) represents the flow rate at the upstream section of the outlet at the first end of the second branch pipe, where m is the flow rate index and L is the flow rate. c D is the horizontal distance from the outlet of the second branch pipe to the main pipe. z Z is the inner diameter of the branch pipe. c0 ζ4 is the elevation difference between the outlet of the second branch pipe and the main pipe, ζ4 is the local resistance coefficient of the valve at the beginning of the branch pipe, and V u (c,1) represents the average inlet velocity of the outlet at the first end of the second branch pipe, g is the acceleration due to gravity, and Q is the average velocity of the outlet. z (c) represents the total outflow rate at the second solenoid valve, H. u (g) represents the actual working pressure at the second intersection of the second branch pipe and the main pipe, H. z (c) represents the total working pressure at the second solenoid valve, and ζ6 represents the local resistance coefficient of the deflector tee valve at the second branch pipe. z (c) represents the average flow rate at the second solenoid valve.

5. The method according to claim 1, characterized in that, The target outflow rate and target working pressure at the second intersection of the first branch pipe and the main pipe are calculated using the following fourth hydraulic calculation formula. Q z (a)=Q u (a,1) Q g (a)=Q z (a) again Substituting, we get: Among them, H u (a,1) represents the working pressure at the upstream section of the outlet of the first branch pipe, f is the pipe friction coefficient, and Q u (a,1) represents the flow rate at the upstream section of the outlet of the first branch pipe, m is the flow rate index, and L... a D is the horizontal distance from the outlet of the first branch pipe to the main pipe. z Z is the inner diameter of the branch pipe. a0 Let V be the elevation difference between the outlet of the first branch pipe and the main pipe, ζ4 be the local resistance coefficient of the valve at the beginning of the branch pipe, and V be the local resistance coefficient of the valve at the beginning of the branch pipe. u (a,1) represents the average inlet velocity of the outlet at the first end of the first branch pipe, g is the acceleration due to gravity, and Q z (a) represents the total outflow rate at the first solenoid valve, H. d (g) is an intermediate parameter, H z (a) represents the total working pressure at the first solenoid valve, Q. g (a) represents the total outflow rate at the first branch pipe, L g D is the distance between the first intersection point and the second intersection point. g Z0 is the inner diameter of the main pipe, Z0 is the height difference between the two ends of the main pipe between the first and second intersections, ζ5 is the sum of the local resistance coefficients of each DC tee on the main pipe between the first and second intersections, and V is the inner diameter of the main pipe. z (a) represents the average flow velocity at the first solenoid valve, V z (c) represents the average flow velocity at the second solenoid valve, ζ7 represents the local resistance coefficient of the DC tee at the second branch pipe, and H u ′(g) is the target working pressure at the second intersection of the first branch pipe and the main pipe.

6. The method according to claim 1, characterized in that, The total gravitational potential energy of the main pipe is calculated using the fifth and sixth calculation formulas below. When there is an intersection between the second branch pipe and the main pipe inlet, the total gravitational potential energy of the main pipe is calculated using the following fifth calculation formula. Q g =Q z (a)+Q z (c) When there is no intersection between the second branch pipe and the main pipe inlet, the total gravitational potential energy of the main pipe is calculated using the following sixth calculation formula. Q g =Q z (a)+Q z (c) Among them, Q g The total outflow rate of the main pipe, Q z (a) represents the total outflow rate at the first solenoid valve, Q. z (c) represents the total outflow rate at the second solenoid valve, Z g ′ is the total gravitational potential energy of the main pipe, H u (g) represents the actual working pressure at the second intersection point of the second branch pipe and the main pipe, f is the pipe friction coefficient, m is the flow index, and L gs D is the length of the pipe section between the water source and the main pipe inlet. gs The inner diameter of the pipe between the water source and the main pipe inlet is given by , b is the pipe diameter index, ζ8 is the local resistance coefficient of the valve at the head of the main pipe, and V is the... g ζ is the average flow velocity of the main pipe, g is the acceleration due to gravity, L0 is the length of the pipe section between the second branch pipe and the inlet of the main pipe, Z1 is the height difference between the second branch pipe and the inlet of the main pipe, and ζ9 is the sum of the local resistance coefficients of each DC three-way valve at the inlet of the main pipe.

7. A hydraulic calculation device for a self-pressurized irrigation system, characterized in that, The device is applied to the self-pressurized irrigation system, which includes: a water source, a first branch pipe, a second branch pipe, and a main pipe of the same diameter. The first and second branch pipes are parallel, and the main pipe is perpendicular to both the first and second branch pipes. Multiple outlets of the same diameter and height are arranged on the first and second branch pipes at preset intervals. A first solenoid valve is installed at a first intersection of the main pipe and the first branch pipe, and a second solenoid valve is installed at a second intersection of the main pipe and the second branch pipe. Water flows from the water source into the main pipe, the first branch pipe, and the second branch pipe. The device includes: The parameter setting module is used to set the initial working pressure of the end outlet of the first branch pipe and the second branch pipe respectively based on a preset range of values. The first calculation module is used to calculate the actual outflow rate and actual working pressure of the end outlet of the first branch pipe and the second branch pipe respectively, based on the hydraulic calculation parameters of the branch pipe and the initial working pressure, and by using the first hydraulic calculation formula. The second calculation module is used to calculate the actual flow rate and actual working pressure of the first outlet adjacent to the end outlet of the first branch pipe and the second branch pipe respectively, based on the hydraulic calculation parameters of the branch pipe and the second hydraulic calculation formula. The third calculation module is used to calculate the actual water flow and actual working pressure of multiple outlets between the first outlet of the first branch pipe and the first solenoid valve, and between the first outlet of the second branch pipe and the second solenoid valve in a recursive manner. The fourth calculation module is used to calculate the actual working pressure at the second intersection of the second branch pipe and the main pipe based on the hydraulic calculation parameters of the branch pipe and using the third hydraulic calculation formula. The fifth calculation module is used to calculate the target working pressure at the second intersection of the first branch pipe and the main pipe based on the hydraulic calculation parameters of the branch pipe and using the fourth hydraulic calculation formula. The first optimization module is used to return the initial working pressure of the end outlet of the first branch pipe and the second branch pipe respectively based on the preset range value if the absolute value of the difference between the actual working pressure and the target working pressure is greater than the first preset threshold. The sixth calculation module is used to calculate the total gravitational potential energy of the main pipe based on the hydraulic calculation parameters of the main pipe and the hydraulic calculation parameters of the branch pipe, using the fifth calculation formula and the sixth calculation formula, if the absolute value of the difference between the actual working pressure and the target working pressure is less than or equal to the first preset threshold. The parameter output module is used to output the actual water flow rate and actual working pressure of multiple outlets of the first branch pipe and the second branch pipe one by one if the absolute value of the difference between the actual gravitational potential energy between the water source and the first inlet of the main pipe and the total gravitational potential energy is less than or equal to a second preset threshold. The second optimization module is used to return to the module that sets the initial working pressure of the end outlets of the first branch pipe and the second branch pipe respectively based on the preset range value if the absolute value of the difference between the actual gravitational potential energy between the water source and the first inlet of the main pipe and the total gravitational potential energy is greater than the second preset threshold.

8. A computer device, characterized in that, include: A memory and a processor are interconnected, the memory stores computer instructions, and the processor executes the computer instructions to perform the hydraulic calculation method for the self-pressurized irrigation system according to any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer instructions for causing the computer to execute the hydraulic calculation method for the self-pressurized irrigation system according to any one of claims 1 to 6.

10. A computer program product, characterized in that, Includes computer instructions for causing a computer to execute the hydraulic calculation method for the self-pressurized irrigation system according to any one of claims 1 to 6.

Images