CN113343439B - Accident identification method for open channel water delivery system - Google Patents

Accident identification method for open channel water delivery system Download PDF

Info

Publication number
CN113343439B
CN113343439B CN202110555791.4A CN202110555791A CN113343439B CN 113343439 B CN113343439 B CN 113343439B CN 202110555791 A CN202110555791 A CN 202110555791A CN 113343439 B CN113343439 B CN 113343439B
Authority
CN
China
Prior art keywords
channel
accident
water level
gate
pool
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
Application number
CN202110555791.4A
Other languages
Chinese (zh)
Other versions
CN113343439A (en
Inventor
管光华
樊甲伟
冯晓波
刘王嘉仪
贾梦洁
朱哲立
郭兴涛
黄一飞
曹佳婷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN202110555791.4A priority Critical patent/CN113343439B/en
Publication of CN113343439A publication Critical patent/CN113343439A/en
Application granted granted Critical
Publication of CN113343439B publication Critical patent/CN113343439B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas or water supply
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling

Abstract

The invention discloses an accident identification method for an open channel water delivery system, which comprises the steps of firstly establishing a simulation model of an actual channel; then, using the water plan as an input quantity, and obtaining the gate target opening, the downstream target water level and the target gate passing flow of each channel pool at different moments through simulation calculation of a simulation model; and taking the water level, the gate opening and the passing flow of the front and back of the gate of each channel pond obtained through actual measurement as input quantities, obtaining four judgment indexes of absolute maximum gate opening deviation, maximum water level deviation speed, end state flow increment and end state water level increment by combining with a target value obtained through simulation calculation, and judging whether an accident, an accident type, an accident channel pond and other information occur according to the four judgment indexes. The obtained accident source information provides basis for emergency scheduling, is favorable for quickly starting an emergency plan, reduces adverse impact of accidents, prevents the accidents from further spreading or causing secondary disasters, and is favorable for guaranteeing the structural safety of high-efficiency water supply and canal system or diversion projects.

Description

Accident identification method for open channel water delivery system
Technical Field
The invention belongs to the field of water conservancy, relates to an automatic channel control technology, and particularly relates to an accident identification method for an open channel water delivery system.
Background
Water resources are very important natural resources and strategic economic resources, and are an organic component of national comprehensive strength. The water resource problem of China is very outstanding, and the per capita water resource quantity of China is 2200m3At present, the per-capita water resource of 16 provinces (cities and districts) is lower than the serious water shortage standard, and the per-capita water resource of 6 provinces and districts is lower than 500m3. Meanwhile, water resources in China also show that annual allocation is uneven in the year, drought and waterlogging disasters occur frequently, and supply and demand contradictions are prominent; secondly, the areas are distributed unevenly, and water and soil resources are not matched; thirdly, the utilization rate of water resources is low, the pollution is serious and the like.
Aiming at the problem of uneven spatial and temporal distribution of water resources in China, China invests a large amount of funds, and numerous water conservancy workers contribute and pay a lot of wisdom and sweat, so that a lot of irrigation channel systems for irrigation and long-distance water diversion projects for relieving uneven regional water resource space distribution are designed and built in China. General survey data of irrigation areas in China show that the designed irrigation area of 30 ten thousand mu and more than 456 parts in the irrigation areas in China reaches 2.80 hundred million mu; the designed irrigation area is 1-30 ten thousand mu irrigation area 7316, and the designed irrigation area is 2.33 hundred million mu; 205.82 ten thousands of irrigated areas of 50-1 ten thousand mu, the designed irrigation area is 3.42 hundred million mu, and the overall scale is very large. At present, China has already built a large number of large diversion projects for solving the contradiction between supply and demand of urban water supply, such as the project of south-to-north water diversion, the project of ' dao ' 28390, the project of ' jin ', the project of ' qing ji qin ', the project of ' huang ' -jin ', the project of Dong-deep water supply and the like, the length of the project is mostly 100-200 km, and the central line of the south-to-north water diversion breaks through 1200 km.
Such irrigation canals and diversion works are mainly open channels as water delivery structures, and are defined as open channel water delivery systems. The open channel water delivery system often causes sudden accidents due to poor operation management, or special geological conditions or adverse meteorological conditions, brings huge challenges to the smooth operation of the system, and even causes great secondary disasters and great economic losses. Taking the project of "guiding river Jihan" as an example, it is recommended that 70% of channels of the line-Longgao I line pass through the expansive soil distribution area, during the excavation and operation of the channels, the side slope may be damaged by expansion and contraction deformation, and further cause instability accidents, and then like the project of south-to-north water transfer, the total main canal coke is made to the cang river, the section about 100km long and the sentry head tunnel of the main line channel section of the south-to-north water transfer have the debris flow, the total main canal high filling channel accumulation 48 sections of the central line channel in the project of south-to-north water transfer and the total length 53.3km, and the high filling channel of the type is easy to have breach accidents under extreme working conditions or rainstorm conditions. Because the irrigation areas in China are widely distributed, the seepage damage accidents such as piping and the like also happen occasionally due to poor geological conditions of part of the irrigation areas. Channel overflow accidents may also occur due to improper operation of gates or influx of rainfall into the channel.
As mentioned above, an emergency accident is inevitable for the open channel water delivery system, and it often causes great loss of lives and properties of people after the emergency accident occurs, so it is very necessary to research the scheduling of the open channel water delivery system under the accident condition, and the current research on the open channel accident mainly focuses on emergency scheduling and hydraulic response characteristics of the accident, and on the problem of tracing to the source of pollutants. Accident emergency dispatch plans are primarily co-dispatch plans for buildings including gates, sluices, weirs, etc., with the expectation of minimizing losses. The concept of incident identification is defined and explained herein for open channel incident identification. The accident identification is to judge the accident type, the accident channel pool and the accident occurrence approximate time according to the measuring point data of the open channel system and the hydraulic response characteristic judgment index under the accident working condition. The method is characterized in that the acquisition of accident source information is a necessary condition for starting an emergency plan or emergency scheduling pertinently, because a proper emergency plan can be provided for emergency administration only by knowing an accident type, an accident channel pool, an accident specific point position and an accident severity, and accident identification is a basic and precondition for emergency scheduling of accidents. Thus, a large amount of time is wasted, the best opportunity for reducing the accident loss through accident scheduling is missed, and the risk that the accident continues to develop into other more serious accidents or secondary disasters are caused is increased. Therefore, the establishment of an on-line accident inversion system of the open channel water delivery system has very important significance for quickly starting an emergency plan, reducing adverse effects of accidents and ensuring the structural safety of the channel system or the diversion project. And secondly, an identification algorithm is established and incorporated into a channel control system, so that the working efficiency is improved, the accuracy and the real-time performance are improved, and the method is also the necessary requirement of intelligent water network construction and water conservancy digital construction in China.
Disclosure of Invention
The invention aims to solve the technical problem of providing an accident identification method of an open channel system aiming at the technical defect of water delivery scheduling of the open channel water delivery system. According to the accident identification method, water level and flow data of a measuring point in an open channel system are used as input, and relevant information of an open channel accident source is identified and obtained, wherein the relevant information mainly comprises information such as whether an accident occurs, the type of the accident, the channel pool where the accident occurs and the like. The obtained accident source information provides basis for emergency scheduling, is favorable for quickly starting an emergency plan, reduces adverse impact of accidents, prevents the accidents from further spreading or causing secondary disasters, and is favorable for guaranteeing the structural safety of high-efficiency water supply and canal system or diversion projects.
The method is realized in such a way that in the daily scheduling process of the open channel water delivery system, the values of target flow, target gate opening and target water level are calculated through a water plan, and then the values are analyzed and calculated by combining with actual measurement data of the open channel to obtain specific numerical values (delta ABE max, delta Zspeed max, delta Qend and delta Zend) of 4 judgment indexes. Due to different types of accidents, hydraulic response characteristics are greatly different and reflect large numerical difference on judgment indexes. And then, by further combining an accident identification method and the numerical value of a specific judgment index, whether an accident occurs or not, the type of the accident and the accident channel pool can be identified and judged. Necessary accident source information is provided for emergency accident scheduling and starting of emergency scheduling plans.
In order to solve the technical problems, the invention adopts the following technical scheme:
an accident identification method for an open channel water delivery system is characterized by comprising the following steps:
step 1, inputting channel modeling parameters and establishing a simulation model of an actual channel;
step 2, taking the water plan as an input quantity, and calculating a gate target opening process, a downstream target water level process and a target gate flow process of each channel pool under a normal operation condition through simulation of a simulation model to obtain the gate target opening, the downstream target water level and the target gate flow of each channel pool at different moments;
step 3, reading the water level, the gate opening and the passing gate flow of each channel pool, which are obtained by actual measurement, and calculating four judgment indexes of each channel pool by combining the target opening, the downstream target water level and the target passing gate flow of each channel pool at different moments, which are obtained in the step 2, wherein the four judgment indexes are respectively an absolute maximum gate opening deviation, a maximum water level deviation speed, a final state flow increment and a final state water level increment, and the method specifically comprises the following steps:
(1) maximum gate opening deviation delta ABE max
Δ ABE max refers to the absolute value of the maximum value of the difference between the actually measured gate opening of all the canal ponds at all the sampling moments and the target opening under the water use plan, and the maximum gate opening deviation calculation formula is as follows:
ΔABE max＝max(max(|Erij-Etij|))
in the above formula, ErijThe opening degree of a gate is measured actually at the moment i, the channel pond j; etijPlanning the opening degree of a target gate at the moment i, of the channel pool j under the water use plan;
(2) maximum water level deviation speed delta Zspeed max
The delta Zspeed max refers to the maximum offset of the water level of the downstream measuring point compared with the water level of the downstream water plan at the unit sampling time step of all the ditch pools at all the sampling moments, and the maximum water level offset speed comprises the positive maximum water level offset speed delta Zspeed max1And a negative maximum water level bias speed Δ Zspeed max2The calculation formula is as follows:
Figure BDA0003077149600000031
Figure BDA0003077149600000032
in the above formula, DT is the calculated time interval, yrijThe measured downstream water level of the i channel pool at the moment j is obtained; y istijPlanning a downstream target water level of the i-channel pool at the time j for water utilization;
(3) end state flow increment delta Qend
ΔQendiThe difference value between the inlet channel flow at the last moment of the channel pool i and the target flow under the water use plan is obtained, and the calculation formula is as follows:
ΔQendi＝Qrend,i-Qtend,i
in the above formula, Qrend,iThe measured flow is the flow into the canal measured at the end of the canal pool i or indirectly calculated by using other measured data; qtend,iPlanning the flow rate of the ditch in the ditch pool i at the last moment for water use;
(4) end state water level increment delta Zend
ΔZendiThe difference value between the downstream water level actually measured at the last moment of the channel pond i and the target water level under the water use plan is obtained, and the calculation formula is as follows:
ΔZendi＝Zrend,i-Ztend,i
in the above formula, Zrend,iIs the downstream water level Z actually measured at the end of the ditch pool itend,iThe water level of the downstream of the ditch pool i at the last moment of water use plan;
step 4, judging the accident channel pool and the accident type according to the four judgment indexes, wherein the specific method comprises the following steps:
step 4.1, judging the size of the maximum gate opening deviation delta ABE max, if the maximum gate opening deviation is larger than a threshold value, determining that the accident type is a gate fault, finding out an accident gate through gate opening and opening change, determining that a corresponding channel pool is an accident channel pool, and if the maximum gate opening deviation is smaller than the threshold value, executing the following steps:
step 4.2, setting a water level fluctuation speed threshold value
Figure BDA0003077149600000041
Judging the maximum water level deviation speed, if the maximum water level deviation speed meets the requirement
Figure BDA0003077149600000042
And is
Figure BDA0003077149600000043
Judging the position of the accident channel pool according to the final-state flow increment index, specifically, calculating the final-state flow increment of each channel pool, and when the upstream final-state flow increment of the n channel pools is positive and the downstream final-state flow increment of the n channel pools is negative, the n channel pools are accident channel pools; otherwise, executing the following steps;
step 4.3, if the maximum water level deviation speed meets the requirement
Figure BDA0003077149600000044
And is
Figure BDA0003077149600000045
Judging the channel overflow accident caused by rainfall inflow, judging the position of the accident channel pool through the final state flow increment index, specifically, calculating the final state flow increment of each channel pool, and when the final state flow increment of the upstream of the n channel pools isThe increment is negative, and the end-state flow increment of the downstream channel pool is positive, then the n channel pool is an accident channel pool; otherwise, calculating the final water level increment Z of each channel pool for the channel silting accidenttendWhen the end state water level increment of the upstream channel pool is positive and the end state water level increment of the downstream channel pool is negative, the n channel pool is an accident channel pool; and completing the accident identification of the open channel water delivery system.
Further, in step 4.1, the method for judging the gate fault includes: traversing all channel ponds gate's of channel last time aperture Erj, judging whether the gate aperture is zero, if be zero, then drop for the gate, corresponding channel pond is accident channel pond, if the gate aperture is nonzero, judge whether the gate aperture changes, if unchangeable, the gate card dies, corresponding channel pond is accident channel pond, if the difference of gate aperture and gate aperture target value is greater than the threshold value, then the gate can not be opened to appointed aperture, corresponding channel pond is accident channel pond, all the other channel ponds are normal channel ponds.
Further, in step 4.2, a downstream water level fluctuation speed critical value epsilon is set to further judge the flow outflow type accident, and when the maximum negative water level deviation speed meets | Delta Zspeedmax2| ≧ epsilon, the accident is a breach structure damage accident; when | zspedmax 2| < epsilon, it is a piping accident.
Further, the water level fluctuation speed threshold value
Figure BDA0003077149600000051
In the range of 2 to 6 x 10-4m/min。
Further, the range of the downstream water level fluctuation speed critical value epsilon is 2-6 multiplied by 10-2m/min。
Further, in step 4.1, the threshold value range of the maximum gate opening deviation is judged to be 1.5-3 cm.
The invention has the following beneficial effects:
(1) the identification algorithm applied to various accident types of the open channel water delivery system does not exist, the established identification algorithm can be used for identifying common open channel accidents such as gate faults, piping accidents, breach or structure damage accidents, channel overflow accidents and channel silting accidents, the technical blank in the field of industry is filled, and reference is provided for the development maturity of the subsequent technology;
(2) the established accident identification method has the functions of automatically identifying the accident type and judging the channel pool in which the accident is positioned, has high identification precision and can provide reference for accident emergency scheduling;
(3) the established recognition algorithm has strong mechanization and small calculation force (the calculation time is less than 1 minute), and is favorable for realizing online recognition;
(4) the established online accident identification method has very important significance for quickly starting an emergency plan, reducing adverse effects of accidents, ensuring the structural safety of a canal system or a diversion project, and being beneficial to improving the working efficiency and the accuracy and the real-time performance.
Drawings
FIG. 1 is a schematic flow diagram of the present invention.
FIG. 2 is an algorithmic logic diagram of a simulation model employed in an embodiment of the present invention.
Note: nac is the number of the accident channel pool, QP1 is the flow boundary at the T moment at the upstream of the channel pool, and QPn is the flow boundary at the T moment at the downstream of the channel pool
Fig. 3 is a logic diagram of the accident recognition method in step 4 according to the embodiment of the present invention.
Detailed Description
The embodiments of the present invention will be described in further detail with reference to the drawings and examples. The following examples are intended to illustrate the invention but are not intended to limit the scope of the invention.
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in further detail below with reference to the accompanying drawings, but the present invention is not limited thereto.
As shown in fig. 1, an accident identification method for an open channel water delivery system includes the following steps:
step 1, inputting channel modeling parameters, and establishing a simulation model of an actual channel, wherein the channel modeling parameters comprise the number of channel pools, the size of a gate, the gradient of the channel, buildings and the roughness, and the model is established on the basis of water delivery channel system operation simulation and control software developed in 2011 of Wangshengde and Guanguanghua.
Step 2, taking the water plan as an input quantity, and calculating a gate target opening process, a downstream target water level process and a target gate flow process of each channel pool under a normal operation condition through simulation of a simulation model to obtain the gate target opening, the downstream target water level and the target gate flow of each channel pool at different moments;
step 3, reading the water level, the gate opening and the passing flow of each channel pond before and after the gate is obtained through actual measurement, and calculating four judgment indexes of each channel pond by combining the gate target opening, the downstream target water level and the target passing flow of each channel pond at different moments obtained in the step 1, wherein the four judgment indexes are respectively the absolute maximum gate opening deviation, the maximum water level deviation speed, the final state flow increment and the final state water level increment, and the method specifically comprises the following steps:
(1) maximum gate opening deviation delta ABE max
Δ ABE max is the absolute value of the maximum of the difference between the actual measured gate opening of all the canals and the target opening of the water plan at all sampling times, and reflects the deviation degree of the gate relative to the target opening under the accident condition, and the value is maintained near zero under the condition of no gate fault without considering feedback. If the value is larger than 2cm, the gate accident is considered to occur, and the maximum gate opening deviation calculation formula is as follows:
ΔABE max＝max(max(|Erij-Etij|))
in the above formula, ErijThe opening (m) of a gate is measured in the I channel pond j at the moment; etijPlanning the target gate opening (m) of the i-channel pool j moment under the water use plan;
(2) maximum water level deviation speed delta Zspeed max
The delta Zspeed max refers to the maximum offset of the water level of the downstream measuring point relative to the downstream water level of the water plan at the unit sampling time step length of all the channel pools at all the sampling moments, and reflects the deviation degree of the downstream water level relative to the target water level under the accident condition; the maximum water level deviation speed comprises a positive maximum water level deviation speed delta Zspeed max1And a negative maximum water level bias speed Δ Zspeed max2The calculation formula is as follows:
Figure BDA0003077149600000061
Figure BDA0003077149600000062
in the above formula, DT is the calculated time interval, yrijThe water level (m) is actually measured at the time j for the i channel pool; y istijPlanning a downstream target water level (m) of the i-channel pool at the time j for water utilization;
(3) end state flow increment delta Qend
ΔQendiThe difference value between the actual measured flow into the channel at the end of the channel pool i or the indirectly calculated flow into the channel through the actual measured water level and the gate opening and the target flow under the water plan reflects the deviation degree of the flow into the channel under the accident condition compared with the target flow, and the calculation formula is as follows:
ΔQendi＝Qrend,i-Qtend,i
in the above formula, Qrend,iThe flow (m) of the inlet channel measured at the end of the channel pool i or indirectly calculated by other measured data3/s)；Qtend,iPlanning the flow (m) of the canal at the end of the water supply3/s)；
(3) End state water level increment delta Zend
ΔZendiThe difference between the downstream water level actually measured at the last moment of the channel pond and the target water level under the water use plan reflects the deviation degree of the channel pond water level relative to the target water level under the accident condition, and the calculation formula is as follows:
ΔZendi＝Zrend,i-Ztend,i
in the above formula, ZrendIs the downstream water level (m), Z measured at the end of the ditch pool itend,iThe water level (m) at the downstream of the ditch pool i at the last time of the water use plan;
step 4, judging the accident channel pool and the accident type according to the four judgment indexes, wherein the specific method comprises the following steps:
step 4.1, judging the size of the maximum gate opening deviation delta ABE max, if the maximum gate opening deviation is larger than a threshold value, determining that the accident type is a gate fault, finding out an accident gate through gate opening and opening change, determining that a corresponding channel pool is an accident channel pool, and if the maximum gate opening deviation is smaller than the threshold value, executing the following steps:
step 4.2, setting a water level fluctuation speed threshold value
Figure BDA0003077149600000071
Judging the maximum water level deviation speed, if the maximum water level deviation speed meets the requirement
Figure BDA0003077149600000072
And is
Figure BDA0003077149600000073
If the water level is reduced, the accident is a flow outflow type accident; and meanwhile, judging the position of the accident channel pool through a final state flow increment index, specifically, calculating the final state flow increment of each channel pool, and for the final state flow increment of each channel pool, when the final state flow increment at the upstream of n channel pools is positive (delta Qend)>0) And the end state flow increment of the downstream ditch pool is negative (delta Qend)<0) The n channel pond is an accident channel pond; otherwise, executing the following steps;
step 4.3, if the maximum water level deviation speed meets the requirement
Figure BDA0003077149600000074
And is
Figure BDA0003077149600000075
If the water level rises, the channel overflow accident caused by rainfall inflow is caused, the position of the accident channel pool is judged according to the final state flow increment index, and specifically, the final state flow of each channel pool is calculated to be positive (delta Qend)>0) When the end state flow increment at the upstream of the n-channel pool is negative (delta Qend)<0) If the end state flow increment of the downstream channel pool is positive, the n channel pools are accident channel pools; otherwise, the water level of the whole open channel water delivery system is increased or decreased (
Figure BDA0003077149600000081
And is
Figure BDA0003077149600000082
) Or (a)
Figure BDA0003077149600000083
And is
Figure BDA0003077149600000084
) Calculating the final water level increment Z of each channel pool for channel silting accidentstendWhen the increment of the final water level at the upstream of the n-channel pool is positive, the upstream is choked (delta Zend)>0) If the end state water level increment of the downstream channel pool is negative, the n channel pools are accident channel pools; and completing the accident identification of the open channel water delivery system.
As a preferred embodiment, in step 4.1, the gate failure determination method is: traversing opening Ern at the end of all canal pond gates of the channel, judging whether the gate opening is zero or not, if zero, falling off the gate, wherein the corresponding canal pond is an accident canal pond, if the gate opening is not zero, judging whether the gate opening is changed or not, if not, blocking the gate, wherein the corresponding canal pond is an accident canal pond, if the difference value between the gate opening and the target value of the gate opening is greater than the threshold value, the gate cannot be opened to the specified opening, the corresponding canal pond is an accident canal pond, and the rest canal ponds are normal canal ponds.
As a preferred embodiment, in step 4.2, a downstream water level fluctuation speed critical value epsilon is set to further judge the flow outflow type accident, and when the maximum negative water level deviation speed meets | Δ Zspeedmax2| ≧ epsilon, the accident is a breach structure damage accident; when | zspedmax 2| < epsilon, it is a piping accident.
The present invention provides two decision thresholds that are,
Figure BDA0003077149600000085
a threshold value of the speed of fluctuation of the water level is indicated,
Figure BDA0003077149600000086
the method indicates that the water level fluctuation is caused by the accident of the open channel water delivery system, and the wrong judgment of the algorithm caused by small fluctuation such as wind waves is avoided.
Figure BDA0003077149600000087
Take 5X 10-4m/min, it is considered that the small fluctuation caused by swell etc. does not exceed 7.2 cm/day. Epsilon represents the critical value between the velocity of the downstream water level fluctuation caused by the small flow change of piping and the large flow fluctuation of breach and structural damage, and the value is an empirical value and is obtained by accumulating more data. After a plurality of attempts and summaries, the invention takes the following values: 5X 10-2The m/min, epsilon values can also be further improved by learning the accident data.
In summary, the recognition algorithm can recognize that: 4 types of accidents of channel overflow and channel silting caused by gate faults, piping, breach or structural damage and rainfall convergence flowing into the channel can be further determined.
And judging the accident channel pool and the accident type by the aid of the established accident identification algorithm. Accident identification is based on significant differences between hydraulic response characteristics of different accident types (see table 1 for details) and can be used to identify channel fouling, piping accidents, breach or structural failure accidents, gate failure accidents, and overflow accidents due to the influx of rainfall into the channel. The difference of the hydraulic response characteristics can also be measured by using the calculated discrimination indexes.
TABLE 1 summary table of simulation rules of different types of accidents
Figure BDA0003077149600000088
Figure BDA0003077149600000091
Note: ↓ represents a small descent or a decrease; ×. either indicates a small rise or increase; ↓ represents a large decrease or a reduction; ↓ [ c ]. means a large rise or increase.
The simulation model is shown in fig. 3, and is based on water delivery channel system operation simulation and control software developed in 2011 by Wang Changde and guan GuanghuaEstablished (Wang Changde, guan Guanghua. water delivery channel system operation simulation and control software [ P ]]China, 2011SR034392, 2011) setting different boundary conditions for different types of accidents can simulate and analyze the hydraulic response characteristics of the channel pond under the accident working condition, and can calculate the gate target passing gate flow, the gate target opening and the downstream target water level under the normal working condition. Under the condition of lacking accident data, the accident condition can be simulated and analyzed by using the model, and then an accident recognition algorithm judgment threshold (a) which can be used for domesticating a more adaptive target open channel water delivery system is obtained
Figure BDA0003077149600000092
Epsilon), the accuracy of the accident identification method is improved.
The accident identification method can identify and obtain the accident source information (accident type, accident channel pool and accident occurrence time) of the open channel according to the water plan and the channel pool measuring point data under the working condition that the open channel has an accident, provides reference for accident emergency scheduling, can effectively avoid time waste caused by on-site accident confirmation, can improve the accuracy and real-time performance of accident emergency scheduling, and has important significance for guaranteeing the structure safety of efficient water supply and engineering.
Appendix description
1. Basic assumptions established by accident recognition algorithm
(1) Considering only the effect of changes in the water plan on the system, i.e., feedforward, when the water plan changes, the target flow rate of the trench gate will also change accordingly, but not the real-time feedback of the check gate.
(2) One type of incident is not considered to develop into other incidents, such as a channel blockage is not considered to further develop into a channel flood incident.
(3) The whole open channel system does not consider additional disturbance, namely the influence of a downstream boundary, and the change of the water taking flow of the water diversion port.
The above embodiments are merely illustrative of the present invention and are not to be construed as limiting the invention. Although the present invention has been described in detail with reference to the embodiments, it should be understood by those skilled in the art that various combinations, modifications or equivalents may be made to the technical solution of the present invention without departing from the spirit and scope of the technical solution of the present invention, and the technical solution of the present invention is covered by the claims of the present invention.

Claims (6)

1. An accident identification method for an open channel water delivery system is characterized by comprising the following steps:
step 1, inputting channel modeling parameters and establishing a simulation model of an actual channel;
step 2, taking the water plan as an input quantity, and calculating a gate target opening process, a downstream target water level process and a target gate flow process of each channel pool under a normal operation condition through simulation of a simulation model to obtain the gate target opening, the downstream target water level and the target gate flow of each channel pool at different moments;
step 3, reading the water level, the gate opening and the passing gate flow of each channel pool, which are obtained by actual measurement, and calculating four judgment indexes of each channel pool by combining the target opening, the downstream target water level and the target passing gate flow of each channel pool at different moments, which are obtained in the step 2, wherein the four judgment indexes are respectively an absolute maximum gate opening deviation, a maximum water level deviation speed, a final state flow increment and a final state water level increment, and the method specifically comprises the following steps:
(1) maximum gate opening deviation delta ABE max
Δ ABE max refers to the absolute value of the maximum value of the difference between the actually measured gate opening of all the canal ponds at all the sampling moments and the target opening under the water use plan, and the maximum gate opening deviation calculation formula is as follows:
ΔABE max＝max(max(|Erij-Etij|))
in the above formula, ErijThe opening degree of a gate is measured actually at the moment i, the channel pond j; etijPlanning the opening degree of a target gate at the moment i, of the channel pool j under the water use plan;
(2) maximum water level deviation speed delta Zspeed max
The delta Zspeed max refers to the maximum offset of the water level of the downstream measuring point compared with the water level of the downstream water plan under the unit sampling time step of all the ditch pools at all the sampling moments, and the maximum water level offset speed comprisesPositive maximum water level deviation delta Zspeed max1And a negative maximum water level bias speed Δ Zspeed max2The calculation formula is as follows:
Figure FDA0003463040120000011
Figure FDA0003463040120000012
in the above formula, DT is the calculated time interval, yrijThe measured downstream water level of the i channel pool at the moment j is obtained; y istijPlanning a downstream target water level of the i-channel pool at the time j for water utilization;
(3) end state flow increment delta Qend
ΔQendiThe difference value between the inlet channel flow at the last moment of the channel pool i and the target flow under the water use plan is obtained, and the calculation formula is as follows:
ΔQendi＝Qrend,i-Qtend,i
in the above formula, Qrend,iThe measured flow is the flow into the canal measured at the end of the canal pool i or indirectly calculated by using other measured data; qtend,iPlanning the flow rate of the ditch in the ditch pool i at the last moment for water use;
(4) end state water level increment delta Zend
ΔZendiThe difference value between the downstream water level actually measured at the last moment of the channel pond i and the target water level under the water use plan is obtained, and the calculation formula is as follows:
ΔZendi＝Zrend,i-Ztend,i
in the above formula, Zrend,iIs the downstream water level Z actually measured at the end of the ditch pool itend,iThe water level of the downstream of the ditch pool i at the last moment of water use plan;
step 4, judging the accident channel pool and the accident type according to the four judgment indexes, wherein the specific method comprises the following steps:
step 4.1, judging the size of the maximum gate opening deviation delta ABE max, if the maximum gate opening deviation is larger than a threshold value, determining that the accident type is a gate fault, finding out an accident gate through gate opening and opening change, determining that a corresponding channel pool is an accident channel pool, and if the maximum gate opening deviation is smaller than the threshold value, executing the following steps:
step 4.2, setting a water level fluctuation speed threshold value
Figure FDA0003463040120000021
Judging the maximum water level deviation speed, if the maximum water level deviation speed meets the requirement
Figure FDA0003463040120000022
And is
Figure FDA0003463040120000023
Judging the position of the accident channel pool according to the final-state flow increment index, specifically, calculating the final-state flow increment of each channel pool, and when the upstream final-state flow increment of the n channel pools is positive and the downstream final-state flow increment of the n channel pools is negative, the n channel pools are accident channel pools; otherwise, executing the following steps;
step 4.3, if the maximum water level deviation speed meets the requirement
Figure FDA0003463040120000024
And is
Figure FDA0003463040120000025
The channel overflow accident caused by rainfall inflow is determined, the position of the accident channel pool is judged according to the final state flow increment index, specifically, the final state flow increment of each channel pool is calculated, and when the final state flow increment of the upstream channel pool of the n channel pools is negative and the final state flow increment of the downstream channel pool is positive, the n channel pools are the accident channel pools; otherwise, calculating the final water level increment Z of each channel pool for the channel silting accidenttendWhen the end state water level increment of the upstream channel pool is positive and the end state water level increment of the downstream channel pool is negative, the n channel pool is an accident channel pool; and completing the accident identification of the open channel water delivery system.
2. The open channel water delivery system accident identification method according to claim 1, wherein in step 4.1, the gate fault judgment method comprises the following steps: traversing opening Ern at the end of all channel pond gates of the channel, judging whether the gate opening is zero or not, if the gate opening is zero, the gate falls off, the corresponding channel pond is an accident channel pond, if the gate opening is not zero, judging whether the gate opening changes or not, if the gate opening does not change, the gate is blocked, the corresponding channel pond is the accident channel pond, if the difference value between the gate opening and the target value of the gate opening is greater than a threshold value, the gate cannot be opened to a specified opening, the corresponding channel pond is the accident channel pond, and the rest channel ponds are normal channel ponds.
3. The method for identifying the accident of the open channel water delivery system according to claim 1, wherein in step 4.2, a downstream water level fluctuation speed critical value epsilon is set to further judge the flow outflow type accident, and when the maximum negative water level deviation speed meets | - Δ Zspeedmax2| ≧ epsilon, the accident is a breach structure damage accident; when | zspedmax 2| < epsilon, it is a piping accident.
4. The open channel water delivery system incident identification method of claim 1, wherein the water level fluctuation speed threshold value
Figure FDA0003463040120000031
In the range of 2 to 6 x 10-4m/min。
5. The method for identifying an accident in an open channel water delivery system according to claim 3, wherein the downstream water level fluctuation speed threshold value ε is in the range of 2-6 x 10-2m/min。
6. The method for identifying the accident of the open channel water delivery system according to claim 1, wherein in step 4.1, the threshold value range of the maximum gate opening deviation is judged to be 1.5-3 cm.
CN202110555791.4A 2021-05-21 2021-05-21 Accident identification method for open channel water delivery system Active CN113343439B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110555791.4A CN113343439B (en) 2021-05-21 2021-05-21 Accident identification method for open channel water delivery system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110555791.4A CN113343439B (en) 2021-05-21 2021-05-21 Accident identification method for open channel water delivery system

Publications (2)

Publication Number Publication Date
CN113343439A CN113343439A (en) 2021-09-03
CN113343439B true CN113343439B (en) 2022-03-04

Family

ID=77470375

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110555791.4A Active CN113343439B (en) 2021-05-21 2021-05-21 Accident identification method for open channel water delivery system

Country Status (1)

Country Link
CN (1) CN113343439B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117824788A (en) * 2024-03-05 2024-04-05 河海大学 Water level monitoring and analyzing system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106223257A (en) * 2016-07-25 2016-12-14 中国水利水电科学研究院 A kind of long distance open channel accident section upstream emergency response gate group control method
CN106647289A (en) * 2017-03-06 2017-05-10 武汉大学 Feedforward control method for opening and closing time of open channel water dispatching gate
CN111784536A (en) * 2020-07-01 2020-10-16 中国水利水电科学研究院 Method for estimating water level overrun time in open channel according to actually measured water level change condition

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106223257A (en) * 2016-07-25 2016-12-14 中国水利水电科学研究院 A kind of long distance open channel accident section upstream emergency response gate group control method
CN106647289A (en) * 2017-03-06 2017-05-10 武汉大学 Feedforward control method for opening and closing time of open channel water dispatching gate
CN111784536A (en) * 2020-07-01 2020-10-16 中国水利水电科学研究院 Method for estimating water level overrun time in open channel according to actually measured water level change condition

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
Research on emergency control mode of sluice gates in water delivery canal;NIE Yan-hua et al.;《MATEC Web of Conferences 246》;20181231;第1-13页 *
基于实测水位和流量数据的渠道控制模型参数辨识与验证;管光华，等;《农业工程学报》;20201208;第36卷(第23期);第92-98页 *
大型调水工程突发事件及应急调度预案初探;段文刚，等;《中国水利学会2008学术会议论文集（下册）》;20081001;第252-256页 *
大型输水明渠运行控制模式研究;崔巍等;《南水北调与水利科技》;20091015;第7卷(第05期);第6-10页 *
跨流域调水工程突发事件及应急调度措施研究;段文刚等;《长江科学院院报》;20100415;第27卷(第04期);第24-27页 *
长距离明渠突发事件应急调度策略设计及应用;韩黎明;《水资源开发与管理》;20181125(第11期);第71-75页 *
长距离渠道闸门故障扰动及小影响应急调度研究;万蕙等;《人民长江》;20180714;第49卷(第13期);第74-78页 *

Also Published As

Publication number Publication date
CN113343439A (en) 2021-09-03

Similar Documents

Publication Publication Date Title
WO2021073192A1 (en) Forecasting and dispatching method by lowering reservoir flood initial dispatch water level in consideration of forecast error
CN108108838B (en) Season regulation reservoir optimal scheduling method with high water utilization rate
CN102032935B (en) Soft measurement method for sewage pumping station flow of urban drainage converged network
CN108764573B (en) Inland arid region-oriented multidimensional equilibrium configuration system for water resources
CN104460577A (en) Water quality safety early warning decision making system
CN112101818B (en) Sponge city flood optimal scheduling method suitable for complex hydraulic connection
CN106355274A (en) Auxiliary system for flood defense dispatching decision
CN113343595B (en) Inversion model of open channel water delivery system accident and method for determining accident flow and position
CN113701060B (en) Intelligent diagnosis method for urban drainage pipe network based on drainage model
CN111046574B (en) Flood control and drainage calculation method for lake and river gate pump system in plain lake region
CN111460742B (en) Method for constructing urban river and lake water network flood control and rainfall flood resource utilization model
CN104615845A (en) Supercritical flood risk analysis method
CN113343439B (en) Accident identification method for open channel water delivery system
CN115271304A (en) Cross-basin water transfer engineering optimal scheduling method based on climate change
CN109002946B (en) Water resource optimization scheduling method for &#39;two-reservoir-two-station&#39; system for replenishing water in rivers and lakes
CN117332908B (en) Multi-objective optimization scheduling method and system for cascade reservoir of coupling set forecast
CN105404760A (en) Rainwater design discharge calculation method in constant and non-uniform flow condition
CN111160763B (en) Safety risk assessment method for hydraulic absorption basin
CN112926786B (en) Shallow lake target water level reverse prediction method and system based on association rule model and numerical simulation
CN102776872B (en) Optimization method of efficient flood preventing power generation in hydropower station
Vereide et al. Case study: Damaging effects of increasing the installed capacity in an existing hydropower plant
CN108959816B (en) Detection method for high arch dam construction early-middle diversion risk rate
CN113239642A (en) Method for calculating reservoir warehousing flow
Wu et al. Multi-objective model of irrigation water distribution based on particle swarm optimization
Zheng et al. Study on Countermeasures for Risks of Flood Resources Utilization in the Three Gorges Project

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