CN111882181B - Quantitative calculation method for water level constraint parameters in automatic control algorithm of water delivery channel - Google Patents
Quantitative calculation method for water level constraint parameters in automatic control algorithm of water delivery channel Download PDFInfo
- Publication number
- CN111882181B CN111882181B CN202010664671.3A CN202010664671A CN111882181B CN 111882181 B CN111882181 B CN 111882181B CN 202010664671 A CN202010664671 A CN 202010664671A CN 111882181 B CN111882181 B CN 111882181B
- Authority
- CN
- China
- Prior art keywords
- water
- channel
- water level
- constraint parameters
- level
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 190
- 238000004364 calculation method Methods 0.000 title claims abstract description 24
- 230000035945 sensitivity Effects 0.000 claims abstract description 42
- 238000004458 analytical method Methods 0.000 claims abstract description 12
- 238000012512 characterization method Methods 0.000 claims abstract description 6
- 230000004044 response Effects 0.000 claims abstract description 5
- 108091006146 Channels Proteins 0.000 claims description 116
- 230000008859 change Effects 0.000 claims description 15
- 238000000034 method Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 230000002262 irrigation Effects 0.000 description 2
- 238000003973 irrigation Methods 0.000 description 2
- 238000005299 abrasion Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Agronomy & Crop Science (AREA)
- Primary Health Care (AREA)
- Mining & Mineral Resources (AREA)
- Animal Husbandry (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Flow Control (AREA)
Abstract
The invention relates to a quantitative calculation method for water level constraint parameters in an automatic control algorithm of a water delivery channel, which comprises the following steps: when the channel water level constraint parameters are analyzed, water diversion sensitivity indexes are introduced, a fitting formula is constructed through hydraulic response analysis, and a quantitative characterization method of the channel water level constraint parameters is provided. Compared with the currently commonly used fixed water level constraint parameters, the invention has the advantages that: (1) Corresponding water level constraint parameters are determined according to different water supply service requirements, and the calculated water level constraint parameters of channels with high water supply service requirements (the accuracy of controlling the flow of the lower-level water distribution channels is high) are small, so that the flow of the lower-level channels can be ensured to be in the range of user requirements. (2) The influence of the water flow working condition on the constraint parameters of the channel water level is considered, for example, under the high-water-level operation of the channel, the influence of the fluctuation of the water level on the flow of the lower-level channel is smaller, namely, under the same flow control precision of the lower-level channel, the constraint parameters of the channel water level under the high-water-level working condition are larger.
Description
Technical Field
The invention relates to a quantitative calculation method for water level constraint parameters in an automatic control algorithm of a water delivery channel, and belongs to the field of irrigation.
Background
With the continuous advancement of the state to modern agriculture and modern irrigation areas, channel automation control becomes a hot spot for water management of agriculture today. The channel automatic control algorithm forms an instruction for coordinating the operation of the automatic equipment through the processing, analysis and decision-making of the information, and is a core of channel automatic control. In the channel control algorithm, water level constraint parameters are important parameters, and if the water level constraint parameters are set too small, frequent adjustment of channel system buildings such as gates and the like is caused, so that energy waste and gate abrasion are caused; if the setting is too large, the formulated control strategy cannot meet the water supply flow requirement of the channel and the water level safety. And according to the hydraulic response research of the existing channel water delivery process, the water level constraint parameters of the upper channel are influenced by the water supply service requirements (the flow control precision of the lower-level water channel) of the channel and the water flow state of the channel. The water level constraint parameters in the current control algorithm are mainly set to fixed values through experience or macroscopic feature analysis, the space-time difference of the water level constraint parameters is not considered, certain theoretical basis and internal cause explanation are lacked, and the channel control effect cannot completely meet the actual requirements.
Disclosure of Invention
Based on the problems, the research discloses a quantitative calculation method of the water level constraint parameters in the channel control algorithm, and by using the method, the water level constraint parameters can be determined according to the water demand of a user in control time and the water flow working condition, and theoretical basis can be provided for improving the automatic control effect of the channel.
In order to solve the technical problems, the invention adopts the following technical scheme:
a quantitative calculation method for water level constraint parameters in an automatic control algorithm of a water delivery channel comprises the following steps:
when the channel water level constraint parameters are analyzed, water diversion sensitivity indexes are introduced, a fitting formula is constructed through hydraulic response analysis, and then a quantitative characterization method of the channel water level constraint parameters is provided.
Preferably: the quantitative calculation method specifically comprises the following steps:
(1) Constructing a water division sensitivity index:
the water diversion sensitivity of the channel is expressed by the ratio of the relative change of the water diversion flow to the absolute change of the water level of the superior channel, and is shown as follows:
wherein S is out(q-h) The sensitivity of the water diversion channel is expressed as the flow change of the water diversion channel when the water level of the upper channel is changed by 1 m;
delta q/q is the relative change of the lower-fraction water flow, and represents the flow control precision of the water diversion channel in channel control;
Δh is the absolute change of the water level of the superior channel, and represents the channel water level constraint parameter in channel control;
(2) Fitting of the water cut sensitivity:
according to the gate flow calculation formula, the gate flow is mainly equal to the upper channel water level (h u ) Width of gate (g) b ) Gate opening (g) e ) And the water depth of the water diversion channel (water level after gate, h) d ) And on the basis of the analysis rules, analyzing the influence rules of the factors on the water diversion sensitivity, and further providing a water diversion sensitivity fitting formula of a research channel:
S out(q-h) ＝f(h u ，g b ，g e ，h d )
h u g is the water level before the gate b For gate width g e Is the opening degree of the gate, h d Is the water level behind the gate;
(3) Quantitative calculation of channel water level constraint parameters:
according to the definition of the water division sensitivity index and a fitting formula, a quantitative calculation method of channel water level constraint parameters is provided:
the invention has the beneficial effects that:
compared with the currently commonly used fixed water level constraint parameters, the channel water level constraint advantages calculated by the method provided by the technical scheme of the invention comprise the following points:
(1) The calculation method considers the water supply service requirements of the lower-level channel, can determine corresponding channel water level constraint parameters according to different water supply service requirements during channel control, and has small calculated water level constraint parameters for channels with high water supply service requirements (high accuracy of flow control of the lower-level water channel), so that the flow of the lower-level channel can be ensured to be in the range of user requirements.
(2) The calculation method considers the influence of the water flow working condition on the constraint parameters of the channel water level, for example, the influence of the fluctuation of the water level on the flow of the lower-level channel is relatively small under the high-water-level operation of the channel, namely, the constraint parameters of the channel water level under the high-water-level working condition are larger under the same flow control precision of the lower-level channel.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions of the prior art, the drawings which are used in the description of the embodiments or the prior art will be briefly described, it being obvious that the drawings in the description below are only some embodiments of the invention, and that other drawings can be obtained according to these drawings without inventive faculty for a person skilled in the art.
FIG. 1 is a graph of influence factors versus water cut sensitivity;
FIG. 2 is a graph showing the comparison of the fit value of the sensitivity index value and the actual calculated value.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Examples
A quantitative calculation method for water level constraint parameters in an automatic control algorithm of a water delivery channel comprises the following steps:
when the channel water level constraint parameters are analyzed, water diversion sensitivity indexes are introduced, a fitting formula is constructed through hydraulic response analysis, and then a quantitative characterization method of the channel water level constraint parameters is provided.
The method specifically comprises the following steps:
(1) Constructing a water division sensitivity index:
the water diversion sensitivity of the channel is expressed by the ratio of the relative change of the water diversion flow to the absolute change of the water level of the superior channel, and is shown as follows:
wherein S is out(q-h) The sensitivity of the water diversion channel is expressed as the flow change of the water diversion channel when the water level of the upper channel is changed by 1 m;
delta q/q is the relative change of the lower-fraction water flow, and represents the flow control precision of the water diversion channel in channel control;
Δh is the absolute change of the water level of the superior channel, and represents the channel water level constraint parameter in channel control;
(2) Fitting of the water cut sensitivity:
according to a gate flow calculation formula, the gate flow is mainly related to the water level of an upper channel, the width of the gate, the opening of the gate and the water depth of a water diversion channel, the factors are used as analysis variables to set different scenes, water diversion sensitivity index values under the different scenes are calculated, the influence rule of each factor on the water diversion sensitivity is analyzed on the basis, and then a water diversion sensitivity fitting formula of a research channel is provided:
S out(q-h) ＝f(h u ，g b ，g e ，h d )；
h u g is the water level before the gate (the water level of the upper channel) b For gate width g e Is the opening degree of the gate, h d Is the water level (water depth of the diversion channel) behind the sluice;
(3) Quantitative calculation of channel water level constraint parameters:
according to the definition of the water division sensitivity index and a fitting formula, a quantitative calculation method of channel water level constraint parameters is provided:
for the water diversion sensitivity index, a certain channel is taken as an example for analysis, and according to a gate flow calculation formula, the gate flow is mainly related to the water level of the superior channel, the gate width, the gate opening and the water depth of the water diversion channel. The following scenes (table 1, 5016 scene combinations in total) are set by referring to the actual values of various parameters of the example channel, the water diversion sensitivity index value of each scene is calculated on the basis of the scenes, and the influence rule of various influence factors on the water diversion sensitivity is analyzed, wherein the height difference of the bottom of the canal at the joint of the upper-level channel and the lower-level water diversion channel is set to be 0.
TABLE 1 Water division sensitivity index analysis scenario Table
Note that: the following scenarios were culled: (1) A situation that the water depth of the superior channel is smaller than that of the water diversion channel; (2) And the opening of the water diversion gate is larger than the depth of the water diversion channel.
TABLE 2 index value of water division sensitivity under different influence factor changes
And analyzing the influence rule of each influence factor on the water diversion sensitivity by setting a scene scheme so as to further provide a fitting formula of the water diversion sensitivity index. The relationship between the depth of the superior channel, the opening of the water diversion gate, the width of the water diversion gate, and the sensitivity of the water diversion channel is shown in FIG. 1.
According to the analysis between the water diversion sensitivity and each influence factor, the water diversion sensitivity and the water diversion gate opening degree of the water diversion gate are in negative correlation, and in positive correlation with the water depth of the water diversion channel. By Matlab software, based on the obtained series of arguments (h u 、h d 、g e ) And dependent variable (S) out(q-h) ) And (3) fitting the relation between each factor and the water division sensitivity index value by taking the nlinfit function in Matlab as a main, so as to obtain a formula:
S out(q-h) ＝-0.84+0.19*h u -0.037*g e +0.9*(h u -h d ) -0.87
according to the comparison (figure 2) between the sensitivity index value fitting value and the actual calculated value, the overall fitting effect is good, and the method can be used for calculating channel water level constraint parameters.
According to a fitting formula of the channel water diversion sensitivity index value, the quantitative characterization formula of the channel real-time water level constraint parameter is obtained as follows:
wherein, the water supply service requirement for the lower channel is h for the constraint parameters of the channel water level u Is the water level g of the superior channel b Is the width g of the gate e Is the opening degree of the gate, h d The water depth of the water diversion channel.
The characterization formula can be used for obtaining the water supply service requirements of different channels and the channel water level constraint parameters under the water flow working condition, as shown in the following table 3.
TABLE 3 channel Water level constraint parameters under different channel Water supply service requirements and Water flow conditions
| Sequence number | h u | h d | g e | Δq/q | Δh |
| 1 | 3 | 2 | 1 | 0.1 | 0.169 |
| 2 | 3 | 2 | 1 | 0.2 | 0.337 |
| 3 | 3 | 2 | 0.5 | 0.1 | 0.164 |
| 4 | 3 | 1.5 | 1 | 0.1 | 0.307 |
| 5 | 4 | 2 | 1 | 0.1 | 0.266 |
Comparing serial numbers 1 and 2 in table 3, the calculation method of the invention considers the water supply service requirements of the lower-level channel, can determine corresponding channel water level constraint parameters according to different water supply service requirements during channel control, and has small calculated water level constraint parameters for channels with high water supply service requirements (high accuracy of flow control of the lower-level water distribution channel), thereby ensuring that the flow of the lower-level channel is in the range of user requirements.
Comparing serial numbers 1 and 5 in table 3, the calculation method of the invention considers the influence of the water flow working condition on the channel water level constraint parameter, for example, under the high water level operation of the channel, the influence of the fluctuation of the water level on the flow of the lower channel is relatively smaller, namely, under the same lower channel flow control precision, the constraint parameter of the channel water level under the high water level working condition is larger.
The foregoing description of the preferred embodiments of the invention is not intended to be limiting, but rather is intended to cover all modifications, equivalents, alternatives, and improvements that fall within the spirit and scope of the invention.
Claims (1)
1. The quantitative calculation method for the water level constraint parameters in the automatic control algorithm of the water delivery channel is characterized by comprising the following steps of:
introducing a water diversion sensitivity index when analyzing the channel water level constraint parameters, constructing a fitting formula of the water diversion sensitivity index through hydraulic response analysis, and further providing a quantitative characterization method of the channel water level constraint parameters; the quantitative calculation method specifically comprises the following steps:
(1) Constructing a water division sensitivity index:
the water diversion sensitivity of the channel is expressed by the ratio of the relative change of the water diversion flow to the absolute change of the water level of the superior channel, and is shown as follows:
wherein S is out(q-h) The sensitivity of the water diversion channel is expressed as the flow change of the water diversion channel when the water level of the upper channel is changed by 1 m;
delta q/q is the relative change of the lower-fraction water flow, and represents the flow control precision of the water diversion channel in channel control;
Δh is the absolute change of the water level of the superior channel, and represents the channel water level constraint parameter in channel control;
(2) Fitting of the water cut sensitivity:
according to the gate flow calculation formula, the gate flow is mainly equal to the upper channel water level (h u ) Width of gate (g) b ) Gate opening (g) e ) And the water depth of the water diversion channel (water level after gate, h) d ) And on the basis of the analysis rules, analyzing the influence rules of the factors on the water diversion sensitivity, and further providing a water diversion sensitivity fitting formula of a research channel:
S out(q-h) ＝f(h u ，g b ，g e ，h d )
h u g is the water level before the gate b For gate width g e Is the opening degree of the gate, h d Is the water level behind the gate;
(3) Quantitative calculation of channel water level constraint parameters:
according to the definition of the water division sensitivity index and a fitting formula, a quantitative calculation method of channel water level constraint parameters is provided:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010664671.3A CN111882181B (en) | 2020-07-10 | 2020-07-10 | Quantitative calculation method for water level constraint parameters in automatic control algorithm of water delivery channel |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010664671.3A CN111882181B (en) | 2020-07-10 | 2020-07-10 | Quantitative calculation method for water level constraint parameters in automatic control algorithm of water delivery channel |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111882181A CN111882181A (en) | 2020-11-03 |
| CN111882181B true CN111882181B (en) | 2023-09-29 |
Family
ID=73151725
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010664671.3A Active CN111882181B (en) | 2020-07-10 | 2020-07-10 | Quantitative calculation method for water level constraint parameters in automatic control algorithm of water delivery channel |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111882181B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110647039A (en) * | 2019-10-08 | 2020-01-03 | 黄河勘测规划设计研究院有限公司 | Synchronous control self-adaptive balance scheduling method for long-distance open channel water delivery project |
| CN110689173A (en) * | 2019-09-12 | 2020-01-14 | 黄河水利委员会黄河水利科学研究院 | Irrigation area agricultural irrigation water demand decision method and system |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040054564A1 (en) * | 2002-09-17 | 2004-03-18 | Fonseca Adolfo M. | Systems and methods for the optimization of resources in energy markets |
-
2020
- 2020-07-10 CN CN202010664671.3A patent/CN111882181B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110689173A (en) * | 2019-09-12 | 2020-01-14 | 黄河水利委员会黄河水利科学研究院 | Irrigation area agricultural irrigation water demand decision method and system |
| CN110647039A (en) * | 2019-10-08 | 2020-01-03 | 黄河勘测规划设计研究院有限公司 | Synchronous control self-adaptive balance scheduling method for long-distance open channel water delivery project |
Non-Patent Citations (1)
| Title |
|---|
| 南水北调中线工程渠道敏感性研究;穆祥鹏;陈文学;崔巍;郭晓晨;王琦;;水力发电学报(第06期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111882181A (en) | 2020-11-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108108838B (en) | Season regulation reservoir optimal scheduling method with high water utilization rate | |
| Shahdany et al. | Operational performance improvements in irrigation canals to overcome groundwater overexploitation | |
| CN112149984B (en) | Reservoir flood regulation multidimensional uncertainty risk analysis method based on Bayesian network | |
| CN110555239B (en) | River surge sluice optimization scheduling method based on drainage pipe network model | |
| CN111533290B (en) | Method for generating optimal operation plan library of sewage treatment process and applying complex scene | |
| CN114417631B (en) | Modeling method of infusion area water delivery and distribution system based on observation data | |
| CN111882181B (en) | Quantitative calculation method for water level constraint parameters in automatic control algorithm of water delivery channel | |
| CN113361214A (en) | Open channel control model parameter identification method based on water level flow data | |
| CN109325610B (en) | Piecewise linear scheduling method for multi-water-source cascade pump station water transfer system | |
| CN114936742A (en) | Water supply system scheduling agent decision method | |
| CN109447336B (en) | Optimized control method for water level between upstream reservoir and reverse regulation reservoir dam thereof | |
| CN116703134B (en) | Multi-target scheduling method and system for large cross-river basin water diversion reservoir | |
| CN117473876A (en) | SHAP method and gene expression programming fused gate group combined control method | |
| CN106489687A (en) | The control method of irrigation rig and control device | |
| CN107808209B (en) | Wind power plant abnormal data identification method based on weighted kNN distance | |
| CN111738493B (en) | Power grid safety risk pre-control decision method and device for new energy power generation uncertainty | |
| CN113359880A (en) | Control method for constant water level filtration of V-shaped filter tank, computer device and computer readable storage medium | |
| CN104977847A (en) | Stable-state condition discrimination method facing atmospheric and vacuum optimization | |
| CN116362552A (en) | Method for evaluating safety risk level of small reservoir | |
| CN114021902B (en) | Reservoir dispatching method for dynamic planning dimension reduction based on dynamic rope collection and discrete mechanism | |
| CN116777027A (en) | Load prediction method and system for abnormal days | |
| CN110729721A (en) | Method for calculating global reserve capacity of power system | |
| JP3573252B2 (en) | Prediction method of water volume in dam | |
| CN110671340B (en) | Flow self-adaptive adjustment method of single-stage full-adjustment pump station | |
| CN113642256A (en) | Practical gate opening and closing optimization method based on theme database |
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 |