CN106759178B - A kind of method of irrigation channel distribution - Google Patents
A kind of method of irrigation channel distribution Download PDFInfo
- Publication number
- CN106759178B CN106759178B CN201611162844.1A CN201611162844A CN106759178B CN 106759178 B CN106759178 B CN 106759178B CN 201611162844 A CN201611162844 A CN 201611162844A CN 106759178 B CN106759178 B CN 106759178B
- Authority
- CN
- China
- Prior art keywords
- channel
- water distribution
- junior
- strip
- period
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02B—HYDRAULIC ENGINEERING
- E02B13/00—Irrigation ditches, i.e. gravity flow, open channel water distribution systems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Abstract
The invention discloses a kind of methods of irrigation channel distribution, belong to irrigation field.The described method includes: establishing encoder matrix P；According to the water distribution initial time t for the j-th strip junior channel that the xth column in the encoder matrix P includejx, distributed water flow ratio ccjx, preset design distributed water flow qjWith the practical water distribution quantity W of preset j-th strip junior channelj, calculate the water distribution end time t of j-th strip junior channeljxy；According to the water distribution initial time t of j-th strip junior channeljxWith water distribution end time tjxy, determine switching variable x of the j-th strip junior channel within the i-th periodijValue；According to the flow ratio of every junior's channel, water distribution initial time, water distribution end time and switching variable, the fitness fit of xth column is calculated；From the maximum column c of encoder matrix P selection fitness, and in water distribution initial time tjcTo water distribution end time tjcyIt is interior according to flow ratio αjcControl j-th strip junior channel corresponds to opening and closing and the lock flow of gate.The present invention can optimize water distribution time and the distributed water flow of channel.
Description
Technical field
The present invention relates to the field of irrigation, in particular to a kind of methods of irrigation channel distribution.
Background technique
Irrigated area water quantity model is mainly the canal engineering and control Project Realization for passing through irrigated area, wherein to meet farmland irrigating
The water for irrigating water demand distributes the main contents distributed rationally for irrigated area water rationally.Water distribution plan is worked out, is exactly being filled entirely
The irrigated area in area, the water intaking time, water intaking water and flow it is fixed under the conditions of draft pouring water every time and distribute to water distribution point
Water, water distribution mode, distributed water flow or water distribution sequence and time.
Irrigation canal water distribution plan is by virtue of experience worked out in current most of irrigated areas, causes biggish canal system water delivery leakage damage
Frequent gate of becoming estranged is adjusted, and with becoming increasingly conspicuous for agricultural water resources imbalance between supply and demand, is proposed to irrigated area water distribution plan higher
Requirement, therefore at present be badly in need of to refine irrigated area water distribution plan, optimize water distribution time and the distributed water flow of channel.
Summary of the invention
In order to solve prior art problem, the present invention provides a kind of methods of irrigation channel distribution.The technical solution is as follows:
The present invention provides a kind of method of irrigation channel distribution, the method is applied to water distribution system, the water distribution system packet
Include K sections of higher level's channels and N junior's channel, K and N are the integer greater than 1, and every section of higher level's canal in K sections of higher level's channel
Road is connected with the part junior channel in N junior channel by gate, which comprises
Establish the encoder matrix P such as (1) formula shown in, the encoder matrix P be include matrix that 2N row and M are arranged, M for greater than
1 integer, jth row include M different flow ratios of j-th strip junior channel, respectively αj1、αj2……αjM, the jth
The flow ratio of junior's channel is used to indicate the actual flow of j-th strip junior channel and the ratio of design discharge, αj1、
αj2……αjMIt is all larger than or is equal to 0.6 and is less than or equal to 1, j=1,2 ... N；N+j row includes the M of j-th strip junior channel
A different water distribution initial time, respectively tj1、tj2……tjM, tj1、tj2……tjMIt is all larger than or equal to 0 and is less than or equal to
T, T are the integer greater than 0；
According to the water distribution initial time t for the j-th strip junior channel that the xth column in the encoder matrix P includejx, water distribution ratio
Value αjxWith the practical water distribution quantity W of preset j-th strip junior channelj, (2) calculate matching for j-th strip junior channel as follows
Water end time tjxy, x=1,2 ... M, qjFor the design distributed water flow of preset j-th strip junior channel；
Wj=αjxqj·(tjxy-tjx)……(2)
According to the water distribution initial time t of j-th strip junior channeljxWith water distribution end time tjxy, determine j-th strip junior channel
Switching variable x within the i-th periodijValue, i=1,2 ... T；
It is determined according to the water distribution initial time for each junior's channel being connected with kth section higher level's channel and water distribution end time
Practical water distribution time t of the kth section higher level channel in the i periodauik, k=1,2 ... K, specific determining method sees below described；
According to every section of higher level's channel the i period practical water distribution time tauik, every junior's channel flow ratio, every
Practical water distribution time t of junior's channel in the i-th periodaijAnd switching variable, (3) to (6) calculate the suitable of xth column as follows
Response fit；
In formula (4), AukFor the default canal bed coefficient of permeability of kth section higher level's channel；mukFor kth section higher level's channel
Default index；QauikPractical distributed water flow for kth section higher level channel in the i-th period,J is in kth section
First junior's channel serial number that grade channel is controlled；LukFor the length of kth section higher level's channel；AjFor j-th strip junior channel
Default canal bed coefficient of permeability；mjFor the default index of j-th strip junior channel；qaijFor j-th strip junior channel the i-th period reality
Border distributed water flow, qaij=αjx·qj；LjFor the length of j-th strip junior channel, taijIt is j-th strip junior channel in the i-th period
The practical water distribution time is the water distribution initial time t according to j-th strip junior channeljxWith water distribution end time tjxyIn the i-th period institute
The real time accounted for determines that specific determining method sees below described；
In formula (5), α is predetermined coefficient, QsFor the preset design discharge of paragraph 1 higher level's channel, j-th strip junior channel
In the practical distributed water flow q of the i-th periodaij=αjxqj, qjFor the preset design discharge of j-th strip junior channel；
In formula (6), β is predetermined coefficient, QskFor the preset design discharge of kth section higher level's channel；
From the maximum column c of encoder matrix P selection fitness, and in water distribution initial time tjcTo the water distribution end time
tjcyIt is interior according to flow ratio αjcThe corresponding gate of j-th strip junior channel is controlled to open.
Optionally, foundation encoder matrix P as shown in (1) formula, comprising:
For j-th strip junior channel, j=1,2 ... N；
It is randomly provided more than or equal to 0.6 and is less than or equal to 1 M different flow ratio, respectively αj1、αj2……
αjM, and by αj1、αj2……αjMJth row as encoder matrix P；
It is randomly provided the M different water distribution initial time more than or equal to 0 and less than or equal to T, respectively tj1、
tj2……tjM, and by tj1、tj2……tjMN+j row as encoder matrix P.
Optionally, the water distribution initial time t according to j-th strip junior channeljxWith water distribution end time tjxy, determine jth
Switching variable x of junior's channel within the i-th periodijValue, comprising:
For the i-th period, i=1,2 ... T, if meeting following first any one condition into third condition,
Switching variable x of the j-th strip junior channel within the i-th period is setijValue be 1；If being unsatisfactory for following first to third item
Switching variable x of the j-th strip junior channel within the i-th period is then arranged in any one condition in partijValue be 0；
First: tjxLess than i and tjxyGreater than i；Second: tjxGreater than i-1 and tjxyLess than i；Third: tjxLess than i-1 and
tjxyGreater than i-1 and it is less than i.
Optionally, the water distribution initial time and water distribution for each junior's channel that the basis is connected with kth section higher level's channel
End time determines kth section higher level channel in the practical water distribution time t of i periodauik, comprising:
When earliest water distribution being selected to originate in the water distribution initial time for each junior's channel being connected with kth section higher level's channel
Between；
At the end of selecting water distribution the latest in the water distribution end time for each junior's channel being connected with kth section higher level's channel
Between；
The real time that the selected water distribution initial time is occupied to the selected water distribution end time in the i period
Section is determined as the practical water distribution time t of kth section higher level's channelauik。
The beneficial effect of the technical scheme provided by the present invention is that:
By the way that encoder matrix is arranged, encoder matrix includes M column, and each column includes the flow ratio and water distribution of every junior's channel
Then initial time calculates the corresponding fitness of every column data that encoder matrix includes, fitness is higher, indicates each channel damage
The water of consumption is just smaller, selects the maximum column of fitness, according to a column of selection, controls the corresponding gate of each junior's channel and exist
When open and opening degree；Further according to the corresponding relationship of each section of higher level's channel and junior's channel, calculate on each section
Practical distributed water flow of the grade channel in the i-th period, so that the preceding paragraph of higher level's channel is controlled to next section of letdown flow, thus
Fining irrigated area water distribution plan is provided, water distribution time and the distributed water flow of channel are optimized.Calculate every level-one junior channel
Water distribution initial time, water distribution end time and after the practical distributed water flow of the i-th period, further according to each section of higher level's channel under
The corresponding relationship of grade channel, calculates each section of higher level's channel in the practical distributed water flow Q of the i-th periodauik, to control higher level
The preceding paragraph of channel is to next section of letdown flow.
Detailed description of the invention
Fig. 1 is a kind of method flow diagram of irrigation channel distribution provided in an embodiment of the present invention；
Fig. 2 is a kind of period structural schematic diagram provided in an embodiment of the present invention；
Fig. 3 is another period structural schematic diagram provided in an embodiment of the present invention；
Fig. 4 is another period structural schematic diagram provided in an embodiment of the present invention；
Fig. 5 is another period structural schematic diagram provided in an embodiment of the present invention；
Fig. 6 is the lateral canal water distribution process provided in an embodiment of the present invention worked out using empirical method；
Fig. 7 is provided in an embodiment of the present invention using fining lateral canal water distribution process；
Fig. 8 is the branch canal water distribution process provided in an embodiment of the present invention worked out using empirical method；
Fig. 9 is provided in an embodiment of the present invention using fining branch canal water distribution process；
Figure 10 is revised branch canal water distribution process provided in an embodiment of the present invention.
Specific embodiment
To make the object, technical solutions and advantages of the present invention clearer, below in conjunction with attached drawing to embodiment party of the present invention
Formula is described in further detail.
It at present include higher level's channel and N junior's channel in irrigated area, N is the integer greater than 1.This bar higher level's channel packet
Include K sections, respectively first segment higher level channel, second segment higher level channel to K sections of higher level's channels；First segment higher level channel, second
The area of section higher level channel to the cross section of K sections of higher level's channels is gradually reduced, and K is the integer greater than 1.The K sections of higher level's channel
In every section of higher level's channel be connected with the part junior channel in N junior's channel by gate, two sections of higher levels of arbitrary neighborhood
Channel is also to be connected by gate.
It is directed to above-mentioned irrigated area, referring to Fig. 1, the embodiment of the invention provides a kind of methods of irrigation channel distribution, comprising:
Step 101: establishing encoder matrix P shown in following (1) formula.
Encoder matrix P is the matrix for including 2N row and M column, and M is the integer greater than 1, and jth row includes j-th strip junior channel
The different flow ratios of M, respectively αj1、αj2……αjM, the flow ratio of j-th strip junior channel is for indicating under j-th strip
The grade actual flow of channel and the ratio of design discharge, αj1、αj2……αjMIt is all larger than or is equal to 0.6 and is less than or equal to 1, j=
1,2……N；N+j row includes M different water distribution initial times of j-th strip junior channel, respectively tj1、tj2……tjM,
tj1、tj2……tjMIt is all larger than or is equal to 0 and is less than or equal to T, T is the integer greater than 0, indicates the default practical water distribution in irrigated area
Time；
Wherein, the corresponding two rows data of every junior's channel in encoder matrix P can be arranged as follows, such as right
In j-th strip junior channel, j=1,2 ... N, j-th strip junior channel corresponds to M data of jth row and M number of N+j row
According to；Detailed setting up procedure is as follows:
It is randomly provided more than or equal to 0.6 and is less than or equal to 1 M different flow ratio, respectively αj1、αj2……
αjM, and by αj1、αj2……αjMJth row as encoder matrix P；
It is randomly provided the M different water distribution initial time more than or equal to 0 and less than or equal to T, respectively tj1、
tj2……tjM, and by tj1、tj2……tjMN+j row as encoder matrix P.
Step 102: according to the water distribution initial time t for the j-th strip junior channel that the xth column in encoder matrix P includejx, match
Water ratio ccjxWith the practical water distribution quantity W of preset j-th strip junior channelj, (2) calculate j-th strip junior channel as follows
Water distribution end time tjxy。
Wj=αjxqj·(tjxy-tjx)……(2)
In above-mentioned formula (2), x=1,2 ... M, qjFor the design distributed water flow of preset j-th strip junior channel.Its
In, j=1,2 ... N, therefore terminated by the water distribution that this step can calculate every junior's channel in N junior's channel
Time.
Step 103: according to the water distribution initial time t of j-th strip junior channeljxWith water distribution end time tjxy, determine j-th strip
Switching variable x of junior's channel within the i-th periodijValue, i=1,2 ... T, j=1,2 ... N.
The detailed process of this step can be with are as follows: for the i-th period, i=1,2 ... T.
(1): working as tjxy> i and tjxWhen < i, then switching variable x of the j-th strip junior channel within the i-th period is setijValue be
1。
tjxy> i and tjx< i refers to the water distribution end time t of j-th strip junior channeljxyNumber i and water distribution greater than the i-th period
Initial time tjxLess than the number i of the i-th period.The situation includes the following two kinds form:
The first, referring to fig. 2, if i-tjx>=1, the practical water distribution time t of the i-th period of j-th strip junior channelaijIt is equal to
One period, i.e. taij=1, x is setij=1.
Second, referring to Fig. 3, if i-tjx< 1, the practical water distribution time t of the i-th period of j-th strip junior channelaij=(i-
t1ij) a period, i.e. taij=(i-t1ij), x is setij=1.
(2): referring to fig. 4, as the water distribution initial time t of j-th strip junior channeljx> i-1 and water distribution end time tjxyWhen < i,
The then practical water distribution time t of the i-th period of j-th strip junior channelaijFor (taij=tjxy-tjx) a period, i.e. taij=tjxy-tjx,
Switching variable x of the j-th strip junior channel within the i-th period is setijValue be 1.
(3): referring to Fig. 5, as water distribution end time i-1 < t of j-th strip junior channeljxy< i and water distribution initial time tjx<i-
When 1: the practical water distribution time t of the i-th period of j-th strip junior channelaijFor (tjxy(i-1)) a period, i.e. taij=tjxy-(i-
1), switching variable x of the setting j-th strip junior channel within the i-th periodijValue be 1.
When (if 1), (2), (3) are all unsatisfactory for: the practical water distribution time t of the i-th period of j-th strip junior channelaij=0,
Switching variable x of the j-th strip junior channel within the i-th period is setijValue be 0.
In this step, due to j=1,2 ... N, every junior's channel is calculated in this step in each period
Interior switching variable xij, all xs of several junior's channels within the i-th period corresponding to kth section higher level's channelijIt is
0, then tauik=0；Otherwise 0 < tauik≤1。
Step 104: according to the water distribution initial time t of j-th strip junior channeljxWith water distribution end time tjxy, determine j-th strip
Junior's channel practical water distribution time t shared by the i-th periodaij, i=1,2 ... T.
This step can be with are as follows: by the water distribution initial time t of j-th strip junior channeljxTo the water distribution end time t of selectionjxy?
The section of real time shared by i-th period is determined as j-th strip junior channel practical water distribution time t shared by the i-th periodaij。
For example, it is assumed that T is 6, that is, there is 6 periods, 4 hours each periods, if the water distribution of j-th strip junior channel rises
Begin time tjxTo the water distribution end time t of selectionjxyIt is the 3rd hour to the 10th hour, then in the 1st period j-th strip junior channel
Shared Duan Weiyi hour of real time is a quarter of entire 1st period, so j-th strip junior channel is in the 1st period
Shared practical water distribution time ta1jIt is 0.25；The Duan Weisi hour of real time shared by the 2nd period j-th strip junior channel,
For entire 2nd period, so j-th strip junior channel practical water distribution time t shared by the 2nd perioda2jIt is 1；In the 3rd period jth
The section of real time shared by junior's channel is two hours, is the 3rd period of half, so j-th strip junior channel exists
Practical water distribution time t shared by 3rd perioda3jIt is 0.5.And the practical water distribution time within the 4th, 5,6 periods is 0.
Step 105: being terminated according to the water distribution initial time for each junior's channel being connected with kth section higher level's channel and water distribution
Time determines kth section higher level channel practical water distribution time t shared by the i-th periodauik, k=1,2 ... K.
This step can be with are as follows: selects most in the water distribution initial time for each junior's channel being connected with kth section higher level's channel
Early water distribution initial time txmin；It is selected most in the water distribution end time for each junior's channel being connected with kth section higher level's channel
The water distribution end time t in eveningxymax；By the water distribution initial time t of selectionxminTo the water distribution end time t of selectionxymaxAt i-th
The real time section of Duan Suozhan is determined as kth section higher level channel practical water distribution time t shared by the i-th periodauik, determination side
The practical water distribution time t of method and the i-th period of j-th strip junior channelaijDetermination method it is identical.
For example, it is assumed that T is still 6, that is, there is 6 periods, 4 hours each periods, if earliest water distribution initial time
txminTo water distribution end time t the latestxymaxIt is the 3rd hour to the 15th hour.So shared by the 1st period kth section grade channel
Duan Weiyi hour of real time, be a quarter of entire 1st period, so kth section higher level channel is shared by the 1st period
Practical water distribution time tau1kIt is 0.25；The Duan Weisi hour of real time shared by the 2nd period kth section higher level's channel, it is
Entire 2nd period, so kth section higher level channel practical water distribution time t shared by the 2nd periodau2kIt is 1；In the 3rd period kth
Duan Weisi hour of real time shared by section grade channel, is entire 3rd period, so kth section higher level channel is in the 3rd period institute
The practical water distribution time t accounted forau3kIt is 1；The Duan Weisan hour of real time shared by the 4th period kth section grade channel is four points
Three the 4th periods, so kth section higher level channel practical water distribution time t shared by the 4th periodau4kIt is 0.75.And the 5th,
The practical water distribution time in 6 periods is 0.
Step 106: according to every section of higher level's channel the i-th period practical water distribution time tauik, every junior's channel reality
Border flow and the flow ratio of design discharge, every junior's channel the i-th period practical water distribution time taijAnd switching variable,
(3) to (6) calculate the fitness fit of xth column as follows；
In formula (4), AukFor the default canal bed coefficient of permeability of kth section higher level's channel；mukFor kth section higher level's channel
Default index；QauikPractical distributed water flow for kth section higher level channel in the i-th period,J be and kth section
First connected junior's channel serial number of higher level's channel；LukFor the length of kth section higher level's channel, k=1,2 ... K；AjFor jth
The default canal bed coefficient of permeability of junior's channel；mjFor the default index of j-th strip junior channel；qaijExist for j-th strip junior channel
The practical distributed water flow of i-th period, qaij=αjx·qj, qjFor the preset design discharge of j-th strip junior channel；LjFor under j-th strip
The length of grade channel；taijPractical water distribution time calculation method for j-th strip junior channel in the i-th period is seen above；tauikIt is
Practical water distribution time of the k sections of higher level's channels within the i-th period refers to that kth section higher level's channel is several met corresponding to it
Water distribution demand of junior's channel in the i-th period, the necessary continuous water distribution time, determination method is seen above；
In formula (5), α is predetermined coefficient, QsFor the preset design discharge of paragraph 1 higher level's channel, qaij=αjxqj, qjFor
The preset design discharge of j-th strip junior channel；
In formula (6), β is predetermined coefficient, QskFor the preset design discharge of kth section higher level's channel；
X=1,2 ... M in this step can calculate the fitness fit of each column in encoder matrix P, fitness
Fit shows that more greatly the canal system water delivery filter loss in irrigated area is smaller, otherwise bigger.
Step 107: using fit as fitness, above-mentioned steps being optimized using genetic algorithm, until obtaining changing for setting
Generation number exports optimal matrix P.
Step 108: from the maximum column c of encoder matrix P selection fitness, and in water distribution initial time tjcTerminate to water distribution
Time tjcyIt is interior according to flow ratio αjcThe corresponding gate of j-th strip junior channel is controlled to open.
Flow ratio αjcBigger, the degree for indicating that the corresponding threshold gate of control j-th strip junior channel is opened is bigger, under j-th strip
Grade channel obtain the actual flow of water just on it is bigger；Conversely, flow ratio αjcIt is smaller, indicate that control j-th strip junior channel is corresponding
The degree opened of threshold gate it is smaller, it is just upper smaller that j-th strip junior channel obtains the actual flow of water.
Next it cites an actual example to select totally 33 winter irrigation of northwest irrigated area junior channel herein to being illustrated herein
Water distribution process data, be shown in Table 1, Fig. 6 and Fig. 8, entire to irrigate the wheel phase be 21 days.
Go down south a grade channel basic parameter in 1 northwest irrigated area of table
The basic document for going down south grade channel is inputted to the water allocation model established herein, wherein canal bed coefficient of permeability A takes 3.4,
Permeable exponent m takes 0.5, and the wheel phase of pouring water is 21 days；Population size takes 1000, and excellent individual takes 100, and the number of iterations 100 makes a variation
Probability takes 0.001, and the water distribution period takes 10h, then water distribution period T is 50；The higher level's channel solved using genetic algorithm is under
Grade irrigation channel distribution process is shown in Fig. 7 and Fig. 9.It is concentrated using the junior's channel and higher level's irrigation channel distribution time of empirical method establishment, is total
Body flow is small, and individual practical distributed water flows of junior's channel are greater than its design discharge or are less than minimum allowable flow；Refine water distribution
The practical distributed water flow of junior's channel of scheme is in the range of flow that it allows, and, water distribution almost the same with design discharge
Rationally, the water distribution time decreases 1d for journey collocation；As shown in Figure 9, it can be seen that the fining practical distributed water flow of higher level's channel is equal
It is even, the adjusting number of higher level's channel gate is effectively reduced, practical water distribution is convenient for.It is computed, the water delivery of optimal water allocation scheme is seeped
The leakage loss mistake water delivery filter loss that relatively water allocation scheme is irrigated in grouping reduces 15%, and show fining water allocation scheme matches water quality
It is high.
Each junior's irrigation channel distribution in fining water allocation scheme starts the period and processing completion time used for them is not quite similar, when practical application
It is required that irrigated area water distribution the degree of automation with higher, therefore, in China, large-scale promotion application this method does not sound feasible at this stage
Border.In view of above situation, herein to refine water allocation scheme as foundation, comprehensive each junior's channel starts water distribution period and water distribution
It lasting, junior's channel is divided into 3 groups of carry out rotation flows, i.e., 1,8,9,11,13,18,22,25, No. 27 junior's channel is first group, 2,
6,7,12,16,17,23,24,26,33,29,32,21, No. 31 junior's channels are second group；28,19,30,5,3,14,20,
15,10, No. 4 junior's channels are third group, wherein No. 4 and No. 25 pouring water for junior's channel last, reality smaller with control area
When water distribution can flexible operating, revised higher level's irrigation channel distribution process is shown in Figure 10.It is computed revised fining water allocation scheme
The more former water allocation scheme of water delivery filter loss reduces 9%, shows that revised fining water allocation scheme is better than former water allocation scheme.
In embodiments of the present invention, by the way that encoder matrix is arranged, encoder matrix includes M column, and each column includes every junior's canal
Then the flow ratio and water distribution initial time in road calculate the corresponding fitness of every column data that encoder matrix includes, adapt to
Degree is higher, indicates that the water of each channel loss is just smaller, selects the maximum column of fitness, according to a column of selection, control is respectively
When the corresponding gate of junior's channel is being opened and opening degree, to provide fining irrigated area water distribution plan, optimizes canal
The water distribution time in road and distributed water flow.
Those of ordinary skill in the art will appreciate that realizing that all or part of the steps of above-described embodiment can pass through hardware
It completes, relevant hardware can also be instructed to complete by program, the program can store in a kind of computer-readable
In storage medium, storage medium mentioned above can be read-only memory, disk or CD etc..
The foregoing is merely presently preferred embodiments of the present invention, is not intended to limit the invention, it is all in spirit of the invention and
Within principle, any modification, equivalent replacement, improvement and so on be should all be included in the protection scope of the present invention.
Claims (3)
1. a kind of method of irrigation channel distribution, which is characterized in that the method is applied to water distribution system, and the water distribution system includes K
Section higher level's channel and N junior's channel, K and N are the integer greater than 1, and every section of higher level's channel in K section higher level's channel and
Part junior channel in N junior channel is connected by gate, which comprises
The encoder matrix P as shown in (1) formula is established, the encoder matrix P is the matrix for including 2N row and M column, and M is greater than 1
Integer, jth row include M different flow ratios of j-th strip junior channel, respectively αj1、αj2……αjM, under the j-th strip
The flow ratio of grade channel is used to indicate the actual flow of j-th strip junior channel and the ratio of design discharge, αj1、
αj2……αjMIt is all larger than or is equal to 0.6 and is less than or equal to 1, j=1,2 ... N；N+j row includes the M of j-th strip junior channel
A different water distribution initial time, respectively tj1、tj2……tjM, tj1、tj2……tjMIt is all larger than or equal to 0 and is less than or equal to
T, T are the integer greater than 0；
According to the water distribution initial time t for the j-th strip junior channel that the xth column in the encoder matrix P includejx, water distribution ratio ccjx
With the practical water distribution quantity W of preset j-th strip junior channelj, (2) calculate the water distribution knot of j-th strip junior channel as follows
Beam time tjxy, x=1,2 ... M, qjFor the design distributed water flow of preset j-th strip junior channel；
Wj=αjxqj·(tjxy-tjx)……(2)
According to the water distribution initial time t of j-th strip junior channeljxWith water distribution end time tjxy, determine j-th strip junior channel i-th
Switching variable x in periodijValue, i=1,2 ... T；
Kth is determined according to the water distribution initial time for each junior's channel being connected with kth section higher level's channel and water distribution end time
Practical water distribution time t of the section higher level's channel in the i-th periodauik, k=1,2 ... K；
According to every section of higher level's channel the i-th period practical water distribution time tauik, every junior's channel flow ratio, under every
Practical water distribution time t of the grade channel in the i-th periodaijAnd switching variable, (3) to (6) calculate the adaptation of xth column as follows
Spend fit；
In formula (4), AukFor the default canal bed coefficient of permeability of kth section higher level's channel；mukFor the default finger of kth section higher level's channel
Number；QauikPractical distributed water flow for kth section higher level channel in the i-th period,J be and kth section higher level's canal
First connected junior's channel serial number of road；LukFor the length of kth section higher level's channel；AjFor the default canal of j-th strip junior channel
Bed coefficient of permeability；mjFor the default index of j-th strip junior channel；qaijFor j-th strip junior channel the i-th period practical water distribution
Flow, qaij=αjx·qj, qjFor the preset design discharge of j-th strip junior channel；LjFor the length of j-th strip junior channel, taijFor
J-th strip junior channel is the water distribution initial time t according to j-th strip junior channel in the practical water distribution time of the i-th periodjxWith match
Water end time tjxyThe real time shared by the i-th period determines；
In formula (5), α is predetermined coefficient, QsFor the preset design discharge of paragraph 1 higher level's channel, j-th strip junior channel is i-th
The practical distributed water flow q of periodaij=αjxqj；
In formula (6), β is predetermined coefficient, QskFor the preset design discharge of kth section higher level's channel；
From the maximum column c of encoder matrix P selection fitness, and in water distribution initial time tjcTo water distribution end time tjcy
It is interior according to flow ratio αjcControl the corresponding gate opening/closing of j-th strip junior channel and lock flow；
Wherein, foundation encoder matrix P as shown in (1) formula, comprising:
For j-th strip junior channel, j=1,2 ... N；
It is randomly provided more than or equal to 0.6 and is less than or equal to 1 M different flow ratio, respectively αj1、αj2……αjM, and
By αj1、αj2……αjMJth row as encoder matrix P；
It is randomly provided the M different water distribution initial time more than or equal to 0 and less than or equal to T, respectively tj1、tj2……tjM,
And by tj1、tj2……tjMN+j row as encoder matrix P.
2. the method as described in claim 1, which is characterized in that the water distribution initial time t according to j-th strip junior channeljx
With water distribution end time tjxy, determine switching variable x of the j-th strip junior channel within the i-th periodijValue, comprising:
For the i-th period, i=1,2 ... T, if meeting following first any one condition into third condition, it is arranged
Switching variable x of the j-th strip junior channel within the i-th periodijValue be 1；If being unsatisfactory for following first into third condition
Any one condition, then switching variable x of the j-th strip junior channel within the i-th period is setijValue be 0；
First: tjxLess than i and tjxyGreater than i；Second: tjxGreater than i-1 and tjxyLess than i；Third: tjxLess than i-1 and tjxyGreatly
In i-1 and it is less than i.
3. the method as described in claim 1, which is characterized in that each junior that the basis is connected with kth section higher level's channel
The water distribution initial time of channel and water distribution end time determine kth section higher level channel in the practical water distribution time t of i periodauik, packet
It includes:
Earliest water distribution initial time is selected in the water distribution initial time for each junior's channel being connected with kth section higher level's channel；
The water distribution end time the latest is selected in the water distribution end time for each junior's channel being connected with kth section higher level's channel；
By the selected earliest water distribution initial time to the selected water distribution end time the latest shared by the i period
Real time is determined as kth section higher level channel in the practical water distribution time t of i periodauik。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611162844.1A CN106759178B (en) | 2016-12-15 | 2016-12-15 | A kind of method of irrigation channel distribution |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611162844.1A CN106759178B (en) | 2016-12-15 | 2016-12-15 | A kind of method of irrigation channel distribution |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106759178A CN106759178A (en) | 2017-05-31 |
| CN106759178B true CN106759178B (en) | 2019-08-13 |
Family
ID=58891430
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611162844.1A Active CN106759178B (en) | 2016-12-15 | 2016-12-15 | A kind of method of irrigation channel distribution |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106759178B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114819735A (en) * | 2022-05-27 | 2022-07-29 | 北京乐水新源智能水务科技有限责任公司 | Irrigation area canal system optimized water distribution method based on incoming water conditions |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107313410B (en) * | 2017-07-13 | 2019-11-29 | 黑龙江省水利水电勘测设计研究院 | Canal water depth design method and device |
| CN114743355B (en) * | 2022-04-12 | 2023-07-21 | 四川水利职业技术学院 | Water diversion engineering leakage alarm method and system |
| CN115981221B (en) * | 2023-03-21 | 2023-09-19 | 北京市农林科学院智能装备技术研究中心 | Channel irrigation gate control method and system optimized step by step |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101622952A (en) * | 2009-08-13 | 2010-01-13 | 中国灌溉排水发展中心 | Water management informatization structure system of irrigation district |
| RU2467560C2 (en) * | 2011-01-14 | 2012-11-27 | Общество с ограниченной ответственностью Научно-производственная фирма "ФИТО" | Irrigation system of nutrient solution |
| CN102094409A (en) * | 2011-03-14 | 2011-06-15 | 中国灌溉排水发展中心 | Method and device for safety identification of canal structures in irrigation areas |
| CN102550370A (en) * | 2012-03-20 | 2012-07-11 | 中国灌溉排水发展中心 | Mobile irrigation method |
| CN105528734A (en) * | 2015-11-27 | 2016-04-27 | 中国科学院寒区旱区环境与工程研究所 | Water volume distribution method based on multilevel water resource management |
-
2016
- 2016-12-15 CN CN201611162844.1A patent/CN106759178B/en active Active
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114819735A (en) * | 2022-05-27 | 2022-07-29 | 北京乐水新源智能水务科技有限责任公司 | Irrigation area canal system optimized water distribution method based on incoming water conditions |
| CN114819735B (en) * | 2022-05-27 | 2023-02-10 | 北京乐水新源智能水务科技有限责任公司 | Irrigation area canal system optimized water distribution method based on incoming water conditions |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106759178A (en) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106759178B (en) | A kind of method of irrigation channel distribution | |
| CN102823402B (en) | Efficient and environment friendly irrigation method for rices in cold region | |
| CN104838868A (en) | Hydrangea cuttage reproduction method | |
| CN105502666A (en) | Construction method of amphibious plant community with water body purification function | |
| CN105660105B (en) | A kind of Rhizoma Gastrodiae implantation methods reducing rotten numb phenomenon | |
| CN107173094A (en) | A kind of warmhouse booth yellow sand substrate grooves and cultural method | |
| CN101248750B (en) | Paddy rice dual-purpose genic male sterile line single-season with regeneration cold water string irrigating seed breeding method | |
| Zhang et al. | Canal delivery and irrigation scheduling optimization based on crop water demand | |
| CN109757341A (en) | A kind of dual tube flow asymmetry drip irrigation method of Extremely arid area grape | |
| CN103741653B (en) | The airtight micro-plot experiment system of a kind of paddy rice | |
| CN105519496A (en) | Natural diffusion myzus persicae vaccination method | |
| CN102084747B (en) | Transplanting method for carex taro | |
| CN205623605U (en) | Controlled drip irrigation system in turn based on drip irrigation technology is store in underground cave | |
| CN206791229U (en) | The irrigation and drainage system of pepper planting | |
| CN104067908A (en) | Water-saving laying method for irrigating and drop irrigation of kiwi fruit | |
| CN105028102B (en) | The cultural method of drought resisting hillside fields banana | |
| CN108450296A (en) | A kind of irrigation method of grassland mining area refuse dump recovering vegetation on slope | |
| CN205320773U (en) | A loop system is irrigated to ground formula of burying for cotton field | |
| CN204811227U (en) | Container seedling loader | |
| CN208047435U (en) | A kind of strawberry cultivating frame | |
| CN206314314U (en) | Cucumber booth plantation irrigation system | |
| CN110036734A (en) | A kind of flue-cured tobacco well cellar formula method for transplanting | |
| Burton | Water management issues at the farm level | |
| CN115222152B (en) | Rotation irrigation system optimization method for improving field drip irrigation uniformity | |
| CN209218848U (en) | A kind of cucumber cultivation system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20221104 Address after: No. 60, Guang'anmen South Street, Xuanwu District, Xicheng District, Beijing 100032 Patentee after: CHINA IRRIGATION AND DRAINAGE DEVELOPMENT CENTER Patentee after: Beijing runhuanong Water Technology Development Co.,Ltd. Address before: 100054 No. 60 Guang'anmen South Street, Xicheng District, Beijing Patentee before: CHINA IRRIGATION AND DRAINAGE DEVELOPMENT CENTER Patentee before: CHINA GREEN WATER INTERNATIONAL CONSULTING CO.,LTD. |
|
| TR01 | Transfer of patent right |