Real-time automatic allocation method for flood discharge gate opening degree of centralized control center

Abstract

The invention discloses a flood discharge gate opening degree real-time automatic distribution method for a centralized control center, and relates to the technical field of water conservancy gate control, and the method comprises the steps: setting a monitoring point in each sub-region, and collecting the mechanical state, the vibration stress parameter and the surrounding geological parameter of a flood discharge gate, so as to construct a gate wear index Msz; therefore, the state of the gate can be monitored in real time, and abrasion and fault problems can be found and early warned in time. And establishing a sub-region prediction model, extracting water flow characteristics and river volume characteristics, calculating an average water flow index Slzs and a regional river clogging index Ydzs of each sub-region, correlating the obtained predicted river load index Fzz, and performing evaluation. And setting a first difference threshold value C1 and a second difference threshold value C2, calculating and obtaining a load difference value FCz according to the predicted river load index Fzz and the load threshold value R, generating an opening instruction, and setting a priority distribution mechanism. The opening degree of the gate is optimized, and the accuracy and response speed of flood discharge control are ensured.

Description

Technical Field

[0001] The invention relates to the technical field of water conservancy gate control, and in particular to a method for real-time automatic distribution of flood discharge gate openings used in a centralized control center. Background Art

[0002] Flood gates are key components in reservoir management. They are used to regulate the water level of the reservoir, ensure the safety of the dam, and manage the water flow in the downstream river. Proper control of the opening of the flood gates is essential for preventing floods, managing water resources, and protecting the downstream ecological environment. The opening of the flood gates directly affects the flow and water level of the downstream river, which in turn affect the sedimentation, siltation, and riverbed stability of the river.

[0003] In traditional technology, the control of flood gates usually relies on manual operation and empirical judgment. This method has many shortcomings, including slow response speed, low control accuracy, and inability to monitor and adjust in real time. Especially in the face of sudden floods or other emergencies, traditional manual operations are often difficult to respond in a timely and accurate manner, which can easily lead to flooding or drought in downstream rivers. In addition, the traditional flood gate control system lacks real-time monitoring and evaluation of the gate's mechanical state, vibration stress parameters, and surrounding geological parameters, making it difficult to detect and handle gate wear and failure in a timely manner, increasing the maintenance cost and safety risks of water conservancy projects.

[0004] Therefore, it is urgent to propose a real-time automatic allocation method for the flood discharge gate opening for the centralized control center. The centralized control center plays a vital role in the management of modern water conservancy projects. By integrating various monitoring, control and management functions, the centralized control center has greatly improved the operation efficiency and safety of water conservancy projects. The centralized control center can monitor and manage various parameters of the reservoir dam, flood discharge valves, flow channels and downstream rivers in real time, and realize precise control of the flood discharge gates through the automation system. This not only improves the operation efficiency of water conservancy projects, but also can respond quickly to emergencies and reduce disaster losses. Summary of the invention

[0005] In view of the deficiencies in the prior art, the present invention provides a method for real-time automatic allocation of flood discharge gate openings for a centralized control center to solve the problems mentioned in the background technology.

[0006] To achieve the above objectives, the present invention is implemented through the following technical solutions: A method for real-time automatic allocation of flood discharge gate openings for a centralized control center, comprising the following steps:

[0007] S1. After pre-processing the selected reservoir dam body, flood discharge valve structure and downstream river data, a valve centralized control 3D model is obtained by drawing using 3D construction software AutoCAD or 3ds-Max. The valve centralized control 3D model includes at least one reservoir dam body and a plurality of gates connected to the reservoir dam body, a plurality of flow path paths corresponding to the gates and a corresponding terminal river channel;

[0008] S2, dividing the valve centralized control three-dimensional model into several sub-areas, each sub-area includes: a gate and a corresponding flow path and a corresponding terminal river channel, and numbering each sub-area;

[0009] S3. Set monitoring points in each sub-area to collect and obtain the mechanical state of the flood discharge gate, vibration stress parameters and geological parameters around the gate, and establish a gate stress data set; and extract the mechanical state of the flood discharge gate, vibration stress parameters and geological parameters around the gate to construct the gate wear index Msz. If the gate wear index Msz exceeds the wear threshold value X, send a first warning instruction to the outside;

[0010] S4. After obtaining and executing the first warning instruction, the gate flood discharge flow-related data is collected through the monitoring points in the sub-area to establish a water flow data set; and the corresponding flow path and the corresponding terminal river data are collected to establish a third data set, and the metadata collected is pre-processed and sent to the database of the centralized control center. After training, the sub-area prediction model is established and trained through neural network technology or time series analysis technology, and the flood discharge flow characteristics and downstream river volume characteristics are extracted. After in-depth analysis, the following are obtained: the average water flow index Slzs and the regional river channel siltation index Ydzs, and the average water flow index Slzs and the regional river channel siltation index Ydzs are associated to calculate and obtain the predicted river channel load index Fzz, and the predicted river channel load index Fzz is compared and analyzed with the load threshold R to obtain the corresponding second warning instruction;

[0011] S5. Preset the first difference threshold C1 and the second difference threshold C2, and calculate the load difference FCz according to the predicted river load index Fzz and the load threshold R. If the load difference FCz is 0, maintain the current sub-area gate opening value. If the load difference FCz is a positive value, compare and analyze the load difference FCz with the first difference threshold C1 to obtain a first corresponding evaluation result, and generate a corresponding opening instruction based on the first corresponding evaluation result; if the load difference FCz is a negative value, compare and analyze the load difference FCz with the second difference threshold C2 to obtain a second corresponding evaluation result, and generate a corresponding opening instruction based on the second corresponding evaluation result, and set the priority allocation mechanism from large to small according to the size of the load difference FCz.

[0012] Preferably, S1 comprises:

[0013] S11. Obtain topographic data of the reservoir area and downstream river channels through GIS or LiDAR scanning, and collect engineering drawings and design parameter data of the reservoir dam, flood discharge gates and related structures;

[0014] S12. Generate a terrain model based on terrain data and engineering drawings using 3D construction software AutoCAD or 3ds-Max, and use 3ds-Max's "Box", "Cylinder" and "Extrude" tools to create the dam body and gate structure based on engineering drawings; use the "Spline" tool to draw the flood discharge channel path, and use the "Loft" or "Extrude" tool to create a 3D channel model;

[0015] S13, integrating the terrain model and the three-dimensional flow channel model to obtain a three-dimensional valve centralized control model;

[0016] S14. Divide the valve centralized control 3D model into several sub-areas according to the positions of gates and flow channels; each sub-area should contain a gate and the corresponding flow channel path and terminal river channel. The valves include spillway gates, bottom outlet flood discharge gates, surface outlet flood discharge gates and power plant flood discharge gates.

[0017] Preferably, in the valve centralized control three-dimensional model, several sub-areas are labeled as: “Fm1, Fm2, Fm3, ..., Fm n "n" represents the number of sub-regions.

[0018] Preferably, S3 includes:

[0019] S31, setting up monitoring points in each sub-area to collect and obtain the mechanical state of the flood discharge gate, vibration stress parameters and geological parameters around the gate, and establish a gate stress data set;

[0020] The gate stress data set includes: current valve opening angle θ, gate operating torque T, gate hydraulic system pressure value P, vibration amplitude A generated during gate operation, gate surface defect area K, gate area geological density M, gate area rock stress N and gate area groundwater level W;

[0021] S32. Preprocessing the gate stress data set, including data cleaning, data calibration, feature extraction and dimensionless processing.

[0022] Preferably, S3 further includes:

[0023] S33, extracting the mechanical state of the flood discharge gate, the vibration stress parameters and the geological parameter characteristics around the gate from the pre-processed gate stress data set to obtain: a mechanical state index Jxz, a vibration stress index Zdy and a gate area geological index Dzz;

[0024] The mechanical state index Jxz, vibration stress index Zdy and gate area geological index Dzz are calculated by the following formula:

[0025] θ f =θ*(1+β1*sin(θ));

[0026] T f =T*(1+β2*log(T));

[0027] P f =P*(1+β3*e -P );

[0028] Jxz=γ1*θ f +γ2*T f +γ3*P f ;

[0029] In the formula, θ f Expressed as the opening angle index, T f Indicates the operating torque index, P f represents the hydraulic system pressure index; θ represents the current valve opening angle, T represents the gate operating torque, P represents the gate hydraulic system pressure value, β1 represents the angle adjustment coefficient, which is used to correct the difference between the actual situation and the ideal state; β2 represents the operating torque adjustment coefficient, β3 represents the gate hydraulic system pressure adjustment coefficient, γ1, γ2 and γ3 are the opening angle index θ respectively. f , operating torque index T f And the hydraulic system pressure index P f The weight coefficient of

[0030] A f =A*(1+β4*A 2 );

[0031] K f =K*(1+β5*e -K );

[0032] Jxz=δ1*A f +δ2*K f ;

[0033] In the formula, A f Represents the vibration amplitude index, K f represents the gate defect index, A represents the vibration amplitude generated during gate operation, K represents the gate surface defect area, β4 represents the vibration amplitude adjustment coefficient, which is used to correct the difference between the actual situation and the ideal state; β5 represents the gate surface defect area adjustment coefficient, δ1 and δ2 are the vibration amplitude index A respectively. f And the gate defect index Kf The weight coefficient of

[0034] M f =M*(1+β6*M);

[0035] N f =N*(1+β7*N 2 );

[0036] W f =W*(1+β8*log(W+1));

[0037] Dzz=α1*M f +α2*N f +α3*W f ;

[0038] Where M f Represents the geological density index, N f Represents rock stress index; W f represents the groundwater level index, M represents the geological density of the gate area, N represents the rock stress of the gate area, W represents the groundwater level in the gate area, β6 represents the geological density adjustment coefficient, which is used to correct the difference between the actual situation and the ideal state; β7 represents the rock stress adjustment coefficient, β8 represents the groundwater level adjustment coefficient, α1, α2 and α3 are the geological density index M respectively. f , rock stress index N f And the groundwater level index W f The weight coefficient; log represents the logarithmic function, and e represents the Sigmoid function;

[0039] S34, the mechanical state index Jxz, the vibration stress index Zdy and the gate area geological index Dzz are correlated with each other according to the following formula to generate the gate wear index Msz:

[0040] Msz=α4*Jxz+α5*Zdy+α6*Dzz;

[0041] Wherein, α4, α5 and α6 are the weight values ​​of the mechanical state index Jxz, the vibration stress index Zdy and the gate area geological index Dzz, respectively, which are adjusted and set by the user, and α4+α5+α6=1;

[0042] S35, preset a wear threshold X, compare and judge the gate wear index Msz with the wear threshold X, and obtain a first judgment result, including:

[0043] When the gate wear index Msz>wear threshold X, it means that the gate wear in the sub-area is abnormal, triggering the first warning instruction, including welding, grinding, lubrication and repair of the gate wear parts, detecting the leakage of hydraulic pipelines, repairing or replacing damaged pipelines, and strengthening the geology of the gate area, including foundation grouting and building retaining walls;

[0044] When the gate wear index Msz ≤ the wear threshold X, it means that the gate wear in the sub-area is normal and continuous monitoring is required.

[0045] Preferably, S4 includes:

[0046] S41, acquiring gate flood discharge flow-related data through monitoring points in the sub-area, and establishing a water flow data set; the water flow data set includes: sub-area water flow velocity Vs, real-time gate discharge flow Q and reservoir real-time water level Sw;

[0047] S42, dividing the river path corresponding to the sub-area and the corresponding terminal river into several river areas at equal intervals according to the length, and marking them as "Hd1, Hd2, Hd3, ..., Hd m ” m represents the number of river areas;

[0048] S43, and collect the corresponding flow path and the corresponding terminal river data to establish a third data set, which includes: regional river water level Xyhd, river cross-sectional area Amj, river sediment concentration Nsnd, river longitudinal slope Spd, river depth Dsd and real-time rainfall JyR.

[0049] Preferably, the S4 further includes:

[0050] S44, after pre-processing the acquired second data set and third data set and sending them to the centralized control center database, a sub-region prediction model is established, and the flood discharge flow characteristics, downstream river volume characteristics and regional river siltation characteristics are extracted. After in-depth analysis, the following are obtained: average water flow index Slzs and regional river siltation index Ydzs:

[0051] The average water flow index Slzs and the regional river siltation index Ydzs are calculated by the following formula:

[0052]

[0053] In the formula, θ f Expressed as the opening angle index, Vs i is the water flow velocity in the ith river channel area of ​​the sub-area, Amj i is the cross-sectional area of ​​the ith river channel region in the sub-region, Spd i represents the longitudinal slope of the ith river channel area in the sub-area, Q iis the real-time outflow of the ith river area in the sub-area, Sw represents the real-time water level of the reservoir, and m is the number of river areas; Nsnd i is the sediment concentration of the ith river channel area in the sub-area; Dsd i is the depth of the ith channel area in the sub-area; Xyhd i is the water level of the ith river channel area in the sub-area; JyR i is the real-time rainfall in the ith river channel area of ​​the sub-area; β9 is the adjustment coefficient of the longitudinal slope of the river channel on the water flow velocity, β 10 is the coefficient that adjusts the effect of water level on sediment deposition;

[0054] S45, and generating a predicted river load index Fzz by associating the average water flow index Slzs and the regional river siltation index Ydzs through the following formula:

[0055] Fzz=α7*Slzs+α8*Ydzs;

[0056] Wherein, α7 and α8 are the weight values ​​of the average water flow index Slzs and the regional river siltation index Ydzs, respectively, which are adjusted and set by the user, and α7+α8=1.

[0057] Preferably, the S4 further comprises:

[0058] S46, preset a load threshold R, and compare and evaluate the predicted river load index Fzz with the load threshold R to obtain a first evaluation result, including:

[0059] When the predicted river load index Fzz is less than the load threshold R, the comparison value between the current valve opening angle θ and the minimum opening threshold is determined. If the current valve opening angle θ is higher than the minimum opening threshold and the real-time reservoir water level Sw is not lower than the minimum threshold, it indicates that there is a risk of siltation in the river in the sub-region, and a second siltation risk warning is issued; if the current valve opening angle θ is higher than the minimum opening threshold, but the real-time reservoir water level Sw is lower than the minimum threshold, it indicates that there is a risk of drought in the reservoir water flow, and a second drought risk warning is issued;

[0060] When the predicted river load index Fzz = load threshold R, it means that the gate opening of the sub-area is normal and the current opening allocation value is maintained;

[0061] When the predicted river load index Fzz>load threshold R, it means that the gate opening in the sub-area is abnormal, and there is a flood risk in the downstream river flow, and the second flood risk warning is sent.

[0062] Preferably, S5 includes:

[0063] S51, when receiving the second drought risk warning and the second flood risk warning, the load difference FCz is calculated according to the predicted river load index Fzz and the load threshold R by the following formula:

[0064] FCz = Fzz-R;

[0065] S52, presetting a first difference threshold C1, if the load difference FCz is a positive value, comparing and analyzing the load difference FCz with the first difference threshold C1, and generating a first gate opening adjustment strategy, including:

[0066] When the load difference FCz is less than the first difference threshold C1, the current gate opening angle is adjusted to be reduced by 3-5%;

[0067] When the first difference threshold C1≤load difference FCz≤first difference threshold C1*130%, the current gate opening angle is adjusted to be reduced by 6-9%;

[0068] When the load difference FCz> the first difference threshold C1*130%, the current gate opening angle is adjusted to decrease by 10-15%;

[0069] S53, presetting a second difference threshold C2, if the load difference FCz is a negative value, comparing and analyzing the load difference FCz with the second difference threshold C2, and generating a second gate opening adjustment strategy, including:

[0070] When the load difference FCz> the second difference threshold C2, the current gate opening angle is adjusted to increase by 3-5%;

[0071] If the second difference threshold C2 ≤ load difference FCz ≤ second difference threshold C2*130%, the current gate opening angle is adjusted to increase by 6-9%;

[0072] When the load difference FCz is less than the second difference threshold C2, the current gate opening angle is adjusted to increase by 10-15%.

[0073] Preferably, the S5 further includes:

[0074] S54, according to steps S52 and S53, according to the size of the load difference FCz, a priority processing mechanism is set from large to small, including:

[0075] The first priority processing sequence is to adjust the current gate opening angle by reducing 10-15% and to increase the current gate opening angle by 10-15%. The second priority processing sequence is to adjust the current gate opening angle by reducing 6-9% and to increase the current gate opening angle by 6-9%. The third priority processing sequence is to adjust the current gate opening angle by reducing 3-5% and to increase the current gate opening angle by 3-5%.

[0076] and sequentially processing the plurality of sub-regions according to a first priority processing sequence, a second priority processing sequence, and a third priority processing sequence;

[0077] S55. When receiving the second siltation risk warning, locate and mark the river area where the regional river siltation index Ydzs>siltation threshold YD in the valve centralized control three-dimensional model, and dispatch an engineering team to clean up the silt in real time until the regional river siltation index Ydzs≤siltation threshold YD.

[0078] The present invention provides a method for real-time automatic allocation of flood discharge gate openings for a centralized control center. It has the following beneficial effects:

[0079] (1) Traditional technologies lack real-time monitoring of the gate's mechanical state, vibration stress parameters, and surrounding geological parameters, making it difficult to detect and handle wear and failure in a timely manner. Monitoring points are set up in each sub-area to collect the mechanical state, vibration stress parameters, and surrounding geological parameters of the flood discharge gate, establish a gate stress data set, and extract features to construct the gate wear index Msz. This allows real-time monitoring of the gate state, timely detection and warning of wear and failure problems, reducing maintenance costs, and improving system safety.

[0080] (2) Traditional flood discharge control relies on manual operation and lacks real-time data analysis and prediction, resulting in slow response to emergencies. Water flow-related data and downstream river data are collected through monitoring points to establish a data set and perform preprocessing. Using neural network or time series analysis technology, a sub-region prediction model is established to extract water flow characteristics and river volume characteristics, calculate the average water flow index Slzs and regional river siltation index Ydzs of each sub-region, and correlate the obtained predicted river load index Fzz and evaluate it. This can achieve accurate control and prediction of flood discharge flow, quickly respond to emergencies, reduce disaster losses, and improve the operating efficiency and safety of the system.

[0081] (3) Traditional methods lack an accurate opening adjustment mechanism and cannot optimize control based on real-time data. Set the first difference threshold C1 and the second difference threshold C2, calculate the load difference FCz based on the predicted river load index Fzz and the load threshold R, generate the opening instruction and set the priority allocation mechanism. In this way, the gate opening can be optimized based on real-time data, ensuring the accuracy and response speed of flood discharge control, improving the operating efficiency of the system, and reducing the uncertainty and risk caused by human intervention.

[0082] (4) This method for real-time automatic allocation of flood discharge gate openings for the centralized control center sets up a priority processing mechanism, and processes several sub-areas in sequence from large to small according to the size of the load difference FCz. By giving priority to key areas, the system can ensure that important areas are adjusted and processed first, thereby improving resource utilization efficiency and emergency response capabilities. This priority processing mechanism not only improves the flexibility of the system, but also reduces resource waste and improves overall operating efficiency. It also locates and identifies siltation areas in real time, and quickly dispatches engineering teams for cleaning. Through this automated siltation management method, it is possible to ensure that the river channel is unobstructed, reduce the impact of siltation on water flow and water conservancy facilities, and ensure the normal operation of the water conservancy system. This not only reduces dependence on manual intervention, improves the automation level of the system, but also reduces errors and risks caused by human operation. BRIEF DESCRIPTION OF THE DRAWINGS

[0083] Figure 1 The present invention is a schematic diagram of the steps of a method for real-time automatic allocation of flood discharge gate openings in a centralized control center. DETAILED DESCRIPTION

[0084] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0085] Example 1

[0086] See also Figure 1 The present invention provides a method for real-time automatic allocation of flood discharge gate opening for a centralized control center, comprising the following steps:

[0087] S1. After pre-processing the selected reservoir dam body, flood discharge valve structure and downstream river data, a valve centralized control 3D model is obtained by drawing using 3D construction software AutoCAD or 3ds-Max. The valve centralized control 3D model includes at least one reservoir dam body and a plurality of gates connected to the reservoir dam body, a plurality of flow path paths corresponding to the gates and a corresponding terminal river channel;

[0088] Traditional flood discharge gate control lacks precise modeling and real-time data integration, making it difficult to monitor and control accurately. : By preprocessing the selected reservoir dam, flood discharge valve structure and downstream river data, and using 3D construction software to create a 3D model, an accurate system representation is provided. This not only intuitively displays the flood discharge system, but also provides an accurate data basis for subsequent monitoring and analysis, improving the overall control accuracy and visual management capabilities of the system.

[0089] S2. Divide the valve centralized control 3D model into several sub-areas. Each sub-area includes: a gate and a corresponding flow path and a corresponding terminal river channel, and each sub-area is numbered. This allows independent monitoring and control of each sub-area, helps with refined management and precise control, and improves the response speed and control accuracy of the system.

[0090] S3. Set monitoring points in each sub-area to collect and obtain the mechanical state of the flood discharge gate, vibration stress parameters and geological parameters around the gate, and establish a gate stress data set; and extract the mechanical state of the flood discharge gate, vibration stress parameters and geological parameters around the gate to construct the gate wear index Msz. If the gate wear index Msz exceeds the wear threshold value X, send a first warning instruction to the outside;

[0091] Traditional technologies lack real-time monitoring of the gate's mechanical state, vibration stress parameters, and surrounding geological parameters, making it difficult to detect and handle wear and failure in a timely manner. Monitoring points are set up in each sub-area to collect the mechanical state, vibration stress parameters, and surrounding geological parameters of the flood discharge gate, establish a gate stress data set, and extract features to construct the gate wear index Msz. This allows real-time monitoring of the gate state, timely detection and warning of wear and failure problems, reducing maintenance costs, and improving system safety.

[0092] S4. After obtaining and executing the first warning instruction, the gate flood discharge flow-related data is collected through the monitoring points in the sub-area to establish a water flow data set; and the corresponding flow path and the corresponding terminal river data are collected to establish a third data set, and the metadata collected is pre-processed and sent to the database of the centralized control center. After training, the sub-area prediction model is established and trained through neural network technology or time series analysis technology, and the flood discharge flow characteristics and downstream river volume characteristics are extracted. After in-depth analysis, the following are obtained: the average water flow index Slzs and the regional river channel siltation index Ydzs, and the average water flow index Slzs and the regional river channel siltation index Ydzs are associated to calculate and obtain the predicted river channel load index Fzz, and the predicted river channel load index Fzz is compared and analyzed with the load threshold R to obtain the corresponding second warning instruction;

[0093] Traditional flood discharge control relies on manual operation and lacks real-time data analysis and prediction, resulting in slow response to emergencies. Water flow-related data and downstream river data are collected through monitoring points to establish a data set and preprocess it. Using neural network or time series analysis technology, a sub-region prediction model is established to extract water flow characteristics and river volume characteristics, calculate the average water flow index Slzs and regional river siltation index Ydzs of each sub-region, and correlate the obtained predicted river load index Fzz and evaluate it. This can achieve accurate control and prediction of flood discharge flow, quickly respond to emergencies, reduce disaster losses, and improve the operating efficiency and safety of the system.

[0094] S5. Preset the first difference threshold C1 and the second difference threshold C2, and calculate the load difference FCz according to the predicted river load index Fzz and the load threshold R. If the load difference FCz is 0, maintain the current sub-area gate opening value. If the load difference FCz is a positive value, compare and analyze the load difference FCz with the first difference threshold C1 to obtain a first corresponding evaluation result, and generate a corresponding opening instruction based on the first corresponding evaluation result; if the load difference FCz is a negative value, compare and analyze the load difference FCz with the second difference threshold C2 to obtain a second corresponding evaluation result, and generate a corresponding opening instruction based on the second corresponding evaluation result, and set the priority allocation mechanism from large to small according to the size of the load difference FCz.

[0095] Traditional methods lack an accurate opening adjustment mechanism and cannot optimize control based on real-time data. Set the first difference threshold C1 and the second difference threshold C2, calculate the load difference FCz based on the predicted river load index Fzz and the load threshold R, generate the opening instruction and set the priority allocation mechanism. In this way, the gate opening can be optimized based on real-time data, ensuring the accuracy and response speed of flood discharge control, improving the operating efficiency of the system, and reducing the uncertainty and risks caused by human intervention.

[0096] In this embodiment, the method of the present invention can realize the real-time automatic allocation of the flood gate opening, improve the accuracy and response speed of flood gate control, and enhance the safety and operation efficiency of water conservancy projects. At the same time, through intelligent analysis and prediction, early warning instructions can be issued in time to reduce flood disasters and waste of water resources. Real-time monitoring and evaluation of the gate mechanical state and surrounding geological parameters can timely detect gate wear and failure, reduce maintenance costs, and improve the safety and operation efficiency of water conservancy projects.

[0097] Embodiment 2, S1 comprises:

[0098] S11. Obtain topographic data of the reservoir area and downstream river channels through geographic information system (GIS) or laser radar (LiDAR) scanning, and collect engineering drawings and design parameter data of the reservoir dam, flood discharge gates and related structures; this step ensures the accuracy and comprehensiveness of the data, and provides a reliable data basis for subsequent model construction.

[0099] S12. Generate a terrain model based on terrain data and engineering drawings using 3D construction software AutoCAD or 3ds-Max, and use 3ds-Max's "Box", "Cylinder" and "Extrude" tools to create the dam and gate structures based on engineering drawings; use the "Spline" tool to draw the flood discharge channel path, and use the "Loft" or "Extrude" tool to create a 3D channel model; this can accurately simulate the actual terrain and structure, providing an intuitive and accurate 3D representation of the system.

[0100] S13. Integrate the terrain model and the three-dimensional flow channel model to obtain the valve centralized control three-dimensional model; through this step, a comprehensive and accurate system model can be formed, making the relationship between each part clearer, providing a good foundation for subsequent monitoring and control.

[0101] S14. Divide the valve centralized control 3D model into several sub-areas according to the location of the gates and flow channels; each sub-area should contain a gate and the corresponding flow channel path and terminal river channel, and the valves include spillway gates, bottom hole flood discharge gates, surface hole flood discharge gates and power plant flood discharge gates. This can achieve independent monitoring and control of each sub-area, which is conducive to refined management and precise control, and improves the response speed and control accuracy of the system.

[0102] In the valve centralized control 3D model, several sub-areas are labeled as: “Fm1, Fm2, Fm3, ..., Fm n "n" represents the number of sub-regions.

[0103] In this embodiment, the valve centralized control 3D model is divided into multiple sub-areas and numbered, so that each sub-area includes a gate and its corresponding flow path. In this way, each sub-area can be independently monitored and controlled, which is conducive to refined management and precise control, and improves the response speed and control accuracy of the system.

[0104] Embodiment 3, S3 comprises:

[0105] S31, setting up monitoring points in each sub-area to collect and obtain the mechanical state of the flood discharge gate, vibration stress parameters and geological parameters around the gate, and establish a gate stress data set;

[0106] The gate stress data set includes: current valve opening angle θ, gate operating torque T, gate hydraulic system pressure value P, vibration amplitude A generated during gate operation, gate surface defect area K, gate area geological density M, gate area rock stress N and gate area groundwater level W;

[0107] The current valve opening angle θ is measured and obtained through the angle sensor, and the gate operating torque T is measured and obtained through the torque sensor; the gate hydraulic system pressure value P is measured and obtained through the pressure sensor; the vibration amplitude A generated during the gate operation is measured and obtained through the acceleration sensor; the gate surface defect area K is measured and obtained through the visual sensor or laser scanner; the geological density M of the gate area is measured and obtained through the geological density sensor; the rock stress N in the gate area is measured and obtained through the strain sensor; the groundwater level W in the gate area is measured and obtained through the water level sensor;

[0108] S32. Preprocessing the gate stress data set, including data cleaning, data calibration, feature extraction and dimensionless processing.

[0109] In this embodiment, by setting up a variety of sensors, the mechanical state of the gate, vibration stress parameters and surrounding geological parameters can be fully and real-time monitored to provide accurate data support for the operation of the gate. Through data cleaning and calibration, the accuracy and consistency of the data can be improved, providing a high-quality data basis for subsequent analysis. Through feature extraction and dimensionless processing, the complexity of the data can be simplified, and the efficiency and effect of data analysis and model training can be improved.

[0110] In embodiment 4, S3 further comprises:

[0111] S33, extracting the mechanical state of the flood discharge gate, the vibration stress parameters and the geological parameter characteristics around the gate from the pre-processed gate stress data set to obtain: a mechanical state index Jxz, a vibration stress index Zdy and a gate area geological index Dzz;

[0112] The mechanical state index Jxz, vibration stress index Zdy and gate area geological index Dzz are calculated by the following formula:

[0113] θ f =θ*(1+β1*sin(θ));

[0114] T f =T*(1+β2*log(T));

[0115] P f =P*(1+β3*e -P );

[0116] Jxz=γ1*θ f +γ2*Tf +γ3*P f ;

[0117] In the formula, θ f Expressed as the opening angle index, T f Indicates the operating torque index, P f represents the hydraulic system pressure index; θ represents the current valve opening angle, T represents the gate operating torque, P represents the gate hydraulic system pressure value, β1 represents the angle adjustment coefficient, which is used to correct the difference between the actual situation and the ideal state; β2 represents the operating torque adjustment coefficient, β3 represents the gate hydraulic system pressure adjustment coefficient, γ1, γ2 and γ3 are the opening angle index θ respectively. f , operating torque index T f And the hydraulic system pressure index P f The weight coefficient of

[0118] A f =A*(1+β4*A 2 );

[0119] K f =K*(1+β5*e -K );

[0120] Jxz=δ1*A f +δ2*K f ;

[0121] In the formula, A f Represents the vibration amplitude index, K f represents the gate defect index, A represents the vibration amplitude generated during gate operation, K represents the gate surface defect area, β4 represents the vibration amplitude adjustment coefficient, which is used to correct the difference between the actual situation and the ideal state; β5 represents the gate surface defect area adjustment coefficient, δ1 and δ2 are the vibration amplitude index A respectively. f And the gate defect index K f The weight coefficient of

[0122] M f =M*(1+β6*M);

[0123] N f =N*(1+β7*N 2 );

[0124] W f =W*(1+β8*log(W+1));

[0125] Dzz=α1*M f +α2*N f +α3*W f ;

[0126] Where M f Represents the geological density index, N f Represents rock stress index; W f represents the groundwater level index, M represents the geological density of the gate area, N represents the rock stress of the gate area, W represents the groundwater level in the gate area, β6 represents the geological density adjustment coefficient, which is used to correct the difference between the actual situation and the ideal state; β7 represents the rock stress adjustment coefficient, β8 represents the groundwater level adjustment coefficient, α1, α2 and α3 are the geological density index M respectively. f , rock stress index N f And the groundwater level index W f The weight coefficient of the gate is logarithmic function, and e represents the sigmoid function. Areas with higher geological density usually have better stability and can provide sufficient support to reduce the settlement and displacement of the gate. On the contrary, low-density geology may cause unstable foundations and increase the risk of settlement of the gate. The geological density will also affect the wear and corrosion of the gate. High-density geological materials may produce greater friction on the gate and increase wear. Low-density geological materials may be more easily eroded, leading to changes in the surrounding geological environment and affecting the service life of the gate. : Changes in rock stress will directly affect the stress distribution of the gate structure. When the rock stress is large, the gate structure may be subjected to additional pressure, resulting in structural stress concentration and increasing the risk of structural damage to the gate. Excessive rock stress may cause geological disasters such as landslides and collapses, which will pose a serious threat to the gate and its surrounding structures and increase safety risks.

[0127] S34, the mechanical state index Jxz, the vibration stress index Zdy and the gate area geological index Dzz are correlated with each other according to the following formula to generate the gate wear index Msz:

[0128] Msz=α4*Jxz+α5*Zdy+α6*Dzz;

[0129] Wherein, α4, α5 and α6 are the weight values ​​of the mechanical state index Jxz, the vibration stress index Zdy and the gate area geological index Dzz, respectively, which are adjusted and set by the user, and α4+α5+α6=1;

[0130] The calculation of the mechanical state index Jxz, vibration stress index Zdy and gate area geological index Dzz takes into account the mechanical state, vibration stress and geological factors, and conducts multi-dimensional analysis, which can more comprehensively evaluate the operating state and potential risks of the gate. It can accurately evaluate the operating state of the gate and provide a basis for decision-making. By calculating the gate wear index Msz, the wear of the gate can be comprehensively evaluated to provide data support for maintenance and maintenance. According to the changes in the gate wear index, the maintenance plan can be dynamically adjusted to improve maintenance efficiency and reduce failure rate.

[0131] S35, preset a wear threshold X, compare and judge the gate wear index Msz with the wear threshold X, and obtain a first judgment result, including:

[0132] When the gate wear index Msz>wear threshold X, it means that the gate wear in the sub-area is abnormal, triggering the first warning instruction, including welding, grinding, lubrication and repair of the gate wear parts, detecting the leakage of hydraulic pipelines, repairing or replacing damaged pipelines, and strengthening the geology of the gate area, including foundation grouting and building retaining walls;

[0133] When the gate wear index Msz ≤ the wear threshold X, it means that the gate wear in the sub-area is normal and continuous monitoring is required.

[0134] In this embodiment, when the gate wear index Msz exceeds the wear threshold X, the system will automatically trigger the first early warning instruction. In this way, problems can be discovered in time before a major failure of the gate occurs, sudden accidents caused by abnormal gate wear can be prevented, and the safety of gate operation can be improved. By timely warning and taking maintenance measures, the sudden failure of the gate due to severe wear can be avoided, thereby ensuring the safety of the reservoir and downstream river. Even if the gate wear index Msz ≤ the wear threshold X, indicating that the gate wear is normal, the system will continue to monitor to ensure that preventive maintenance is carried out at the early stage of wear, reduce the degree of wear, and extend the service life of the gate. When the gate wear index Msz> the wear threshold X, the system will trigger specific maintenance instructions, such as welding, grinding, lubrication repair, etc., to ensure the pertinence and effectiveness of maintenance work and reduce the waste of resources caused by blind maintenance.

[0135] Embodiment 5, S4 comprises:

[0136] S41, acquiring gate flood discharge flow-related data through monitoring points in the sub-area, and establishing a water flow data set; the water flow data set includes: sub-area water flow velocity Vs, real-time gate discharge flow Q and reservoir real-time water level Sw;

[0137] Water velocity Vs: measured using a Doppler flow meter or a current meter.

[0138] Real-time outflow water flow Q: measured using ultrasonic flowmeter or turbine flowmeter.

[0139] Reservoir real-time water level Sw: measured using a water level sensor or laser rangefinder.

[0140] Step S41 can obtain key parameters such as water flow velocity, outflow water flow and reservoir water level in the sub-area in real time, ensuring that the system has a comprehensive grasp of the current water flow conditions.

[0141] S42. Divide the river path corresponding to the sub-area and the corresponding terminal river into several river areas at equal intervals according to their lengths, and mark them in the valve centralized control 3D model as: "Hd1, Hd2, Hd3, ..., Hdm" m represents the number of several river areas; dividing the river into several river areas makes the management of the river more detailed and accurate. Independent analysis and monitoring can be carried out for different river areas, because each river area is different, and the specific conditions of each area are obtained, which is convenient for formulating regional management and maintenance measures. Through fine division, the monitoring accuracy and control accuracy of the overall system are improved, which is conducive to better water flow management and flood control scheduling.

[0142] S43, and collect the corresponding flow path and the corresponding terminal river data to establish a third data set, which includes: regional river water level Xyhd, river cross-sectional area Amj, river sediment concentration Nsnd, river longitudinal slope Spd, river depth Dsd and real-time rainfall JyR.

[0143] Regional river water level Xyhd: measured using a water level sensor or a floating water level gauge.

[0144] River cross-sectional area Amj: calculated by topographic survey instruments (such as GPS and total station) combined with flow measurement data.

[0145] River sediment concentration Nsnd: measured using an optical turbidity meter or an acoustic suspended matter concentration meter.

[0146] River channel longitudinal slope Spd: calculated after measuring terrain data using a laser rangefinder or total station.

[0147] River channel depth Dsd: measured using an echo sounder or laser depth sounder.

[0148] Real-time rainfall JyR: measured using a rain gauge or radar rainfall sensor.

[0149] By collecting various data of the river, a comprehensive third data set is established, covering important parameters such as water level, cross-sectional area, sediment concentration, longitudinal slope, depth and rainfall. Based on the comprehensive data set, scientific analysis and calculation can be carried out to help predict and evaluate the flow changes and sediment deposition of the river. Through the analysis of the third data set, the river management strategy can be optimized, the sediment deposition can be reduced, the river blockage can be prevented, and the river can be kept unobstructed and stable.

[0150] In embodiment 6, the S4 further comprises:

[0151] S44, after pre-processing the acquired second data set and third data set and sending them to the centralized control center database, a sub-region prediction model is established, and the flood discharge flow characteristics, downstream river volume characteristics and regional river siltation characteristics are extracted. After in-depth analysis, the following are obtained: average water flow index Slzs and regional river siltation index Ydzs:

[0152] The average water flow index Slzs and the regional river siltation index Ydzs are calculated by the following formula:

[0153]

[0154] In the formula, θ f Expressed as the opening angle index, Vs i is the water flow velocity in the ith river channel area of ​​the sub-area, Amj i is the cross-sectional area of ​​the ith river channel region in the sub-region, Spd i represents the longitudinal slope of the ith river channel area in the sub-area, Q i is the real-time outflow of the ith river area in the sub-area, Sw represents the real-time water level of the reservoir, and m is the number of river areas; Nsnd i is the sediment concentration of the ith river channel area in the sub-area; Dsd i is the depth of the ith channel area in the sub-area; Xyhd i is the water level of the ith river channel area in the sub-area; JyR i is the real-time rainfall in the ith river channel area of ​​the sub-area; β9 is the adjustment coefficient of the longitudinal slope of the river channel on the water flow velocity, β 10 is the coefficient that adjusts the effect of water level on sediment deposition;

[0155] S45, and generating a predicted river load index Fzz by associating the average water flow index Slzs and the regional river siltation index Ydzs through the following formula:

[0156] Fzz=α7*Slzs+α8*Ydzs;

[0157] Wherein, α7 and α8 are the weight values ​​of the average water flow index Slzs and the regional river siltation index Ydzs, respectively, which are adjusted and set by the user, and α7+α8=1.

[0158] In this embodiment, by extracting flood discharge flow characteristics, downstream river volume characteristics and regional river siltation characteristics, the system can comprehensively analyze the water flow and sediment deposition conditions to provide a scientific basis for decision-making. Based on the pre-processed data and sub-region prediction model, the system can predict future water flow conditions and river siltation conditions to improve the foresight of emergency treatment and resource management. The average water flow index Slzs and the regional river siltation index Ydzs provide quantitative indicators of water flow and sediment conditions, helping to optimize the opening control of the flood discharge gate and avoid excessive siltation and excessive water flow caused by sub-regional gate opening.

[0159] In embodiment 7, the S4 further comprises:

[0160] S46, preset a load threshold R, and compare and evaluate the predicted river load index Fzz with the load threshold R to obtain a first evaluation result, including:

[0161] When the predicted river load index Fzz is less than the load threshold R, the comparison value between the current valve opening angle θ and the minimum opening threshold is determined. If the current valve opening angle θ is higher than the minimum opening threshold and the real-time reservoir water level Sw is not lower than the minimum threshold, it indicates that there is a risk of siltation in the river in the sub-region, and a second siltation risk warning is issued; if the current valve opening angle θ is higher than the minimum opening threshold, but the real-time reservoir water level Sw is lower than the minimum threshold, it indicates that there is a risk of drought in the reservoir water flow, and a second drought risk warning is issued;

[0162] When the predicted river load index Fzz = load threshold R, it means that the gate opening of the sub-area is normal and the current opening allocation value is maintained;

[0163] When the predicted river load index Fzz>load threshold R, it means that the gate opening in the sub-area is abnormal, and there is a flood risk in the downstream river flow, and the second flood risk warning is sent.

[0164] In this embodiment, by evaluating the predicted river load index Fzz, the system can accurately determine the current river load status and issue targeted warnings according to different situations, including siltation risk, drought risk and flood risk, thereby improving the accuracy and timeliness of the warning. According to the real-time reservoir water level and valve opening angle, the system can dynamically adjust the opening distribution value to ensure that the river load is within a reasonable range, effectively prevent and respond to problems such as river siltation, drought and flood, and improve the flexibility and response speed of regulation.

[0165] In embodiment 8, the S5 comprises:

[0166] S51, when receiving the second drought risk warning and the second flood risk warning, the load difference FCz is calculated according to the predicted river load index Fzz and the load threshold R by the following formula:

[0167] FCz = Fzz-R;

[0168] S52, presetting a first difference threshold C1, if the load difference FCz is a positive value, comparing and analyzing the load difference FCz with the first difference threshold C1, and generating a first gate opening adjustment strategy, including:

[0169] When the load difference FCz is less than the first difference threshold C1, the current gate opening angle is adjusted to be reduced by 3-5%;

[0170] When the first difference threshold C1≤load difference FCz≤first difference threshold C1*130%, the current gate opening angle is adjusted to be reduced by 6-9%;

[0171] When the load difference FCz> the first difference threshold C1*130%, the current gate opening angle is adjusted to decrease by 10-15%;

[0172] S53, presetting a second difference threshold C2, if the load difference FCz is a negative value, comparing and analyzing the load difference FCz with the second difference threshold C2, and generating a second gate opening adjustment strategy, including:

[0173] When the load difference FCz> the second difference threshold C2, the current gate opening angle is adjusted to increase by 3-5%;

[0174] If the second difference threshold C2 ≤ load difference FCz ≤ second difference threshold C2*130%, the current gate opening angle is adjusted to increase by 6-9%;

[0175] When the load difference FCz is less than the second difference threshold C2, the current gate opening angle is adjusted to increase by 10-15%.

[0176] Data example: Current gate opening angle: 50%;

[0177] Load difference FCz: calculated load difference;

[0178] First difference threshold C1: 10;

[0179] Second difference threshold C2: -10;

[0180] Example 1: When the load difference is positive;

[0181] Assume that the calculated load difference FCz is 15.

[0182] When the load difference FCz=15:

[0183] FCz<first difference threshold C1: not satisfied;

[0184] First difference threshold C1≤FCz≤first difference threshold C1*130%(10≤15≤13): not satisfied;

[0185] FCz>first difference threshold C1*130%(15>13): satisfied;

[0186] According to the adjustment strategy:

[0187] Adjust to reduce the current gate opening angle by 7-10%.

[0188] Current opening angle: 50%;

[0189] Opening angle after adjustment: 50%-(7-10)%=45%-43%;

[0190] Example 2: When the load difference is negative;

[0191] Assume that the calculated load difference FCz is -20.

[0192] When the load difference FCz=-20:

[0193] FCz> second difference threshold C2 (-20>-10): not satisfied;

[0194] Second difference threshold C2≤FCz≤second difference threshold C2*130%(-10≤-20≤-13): not satisfied;

[0195] FCz<second difference threshold C2(-20<-10): satisfied;

[0196] According to the adjustment strategy:

[0197] Adjustment increases the current gate opening angle by 10-15%.

[0198] Current opening angle: 50%;

[0199] Opening angle after adjustment: 50% + (10-15)% = 55% - 57.5%;

[0200] S54, according to steps S52 and S53, according to the size of the load difference FCz, a priority processing mechanism is set from large to small, including:

[0201] The first priority processing sequence is to adjust the current gate opening angle by reducing 10-15% and to increase the current gate opening angle by 10-15%. The second priority processing sequence is to adjust the current gate opening angle by reducing 6-9% and to increase the current gate opening angle by 6-9%. The third priority processing sequence is to adjust the current gate opening angle by reducing 3-5% and to increase the current gate opening angle by 3-5%.

[0202] and sequentially processing the plurality of sub-regions according to a first priority processing sequence, a second priority processing sequence, and a third priority processing sequence;

[0203] S55. When receiving the second siltation risk warning, locate and mark the river area where the regional river siltation index Ydzs>siltation threshold YD in the valve centralized control three-dimensional model, and dispatch an engineering team to clean up the silt in real time until the regional river siltation index Ydzs≤siltation threshold YD.

[0204] Priority Processing Example

[0205] Assume that the system needs to process multiple sub-areas, and the load difference of each sub-area is as follows:

[0206] Sub-area A: load difference FCz = 25

[0207] Sub-area B: load difference FCz = -5

[0208] Sub-area C: load difference FCz = 12

[0209] According to the priority processing mechanism:

[0210] First priority processing sequence:

[0211] Sub-area A: Adjust to reduce the current gate opening angle by 10-15%.

[0212] Second priority processing sequence:

[0213] Sub-area C: Adjust to reduce the current gate opening angle by 4-6%.

[0214] Third priority processing sequence:

[0215] Sub-area B: Adjustment increases the current gate opening angle by 3-5%.

[0216] In this embodiment, by calculating the load difference FCz and the preset difference thresholds C1 and C2, the system can accurately adjust the gate opening angle to ensure the reasonable distribution and control of water flow, and avoid the risk of drought or flood caused by excessive or insufficient opening adjustment. The method provides a dynamic response mechanism, which adjusts the gate opening angle in real time according to the actual situation, improves the flexibility and response speed of the system, and can more effectively respond to emergencies and ensure the safety of reservoirs and downstream rivers. By setting a priority processing mechanism, the system can process each sub-area according to its importance and urgency, ensuring that key areas are regulated and processed first, and improving resource utilization efficiency and emergency response capabilities. When the risk of siltation is detected, the system can locate and identify the siltation area in real time, quickly dispatch an engineering team to clean it up, ensure the smooth flow of the river, reduce the impact of siltation on water flow and water conservancy facilities, and ensure the normal operation of the water conservancy system. Through automated and intelligent adjustment strategies, the dependence on manual intervention is reduced, the automation level and operation efficiency of the system are improved, and the errors and risks caused by human operation are reduced. The system provides comprehensive decision-making support by comprehensively analyzing multiple parameters (such as water flow velocity, water level, siltation index, etc.), helping managers make more scientific and reasonable decisions and ensuring the safe and stable operation of water conservancy facilities.

[0217] The threshold is set to facilitate comparison. The size of the threshold depends on the amount of sample data and the number of bases set by technicians in this field for each group of sample data; as long as it does not affect the proportional relationship between the parameter and the quantized value.

[0218] The above formulas are obtained by collecting a large amount of data for software simulation and selecting a formula that is close to the actual value. The coefficients in the formula are set by technical personnel in this field according to actual conditions. The above is only a preferred specific implementation manner of the present invention, but the protection scope of the present invention is not limited to this. Any technical personnel familiar with the technical field within the technical scope disclosed by the present invention, according to the technical solution and the inventive concept of the present invention, make equivalent replacement or change, which should be covered within the protection scope of the present invention.

Claims

1. A method for real-time automatic allocation of flood gate opening for a centralized control center, characterized in that: The following steps are included: S1. After pre-processing the selected reservoir dam body, flood discharge valve structure and downstream river data, a valve centralized control 3D model is obtained by drawing using 3D construction software AutoCAD or 3ds-Max. The valve centralized control 3D model includes at least one reservoir dam body and a plurality of gates connected to the reservoir dam body, a plurality of flow path paths corresponding to the gates and a corresponding terminal river channel; S2, dividing the valve centralized control three-dimensional model into several sub-areas, each sub-area includes: a gate and a corresponding flow path and a corresponding terminal river channel, and numbering each sub-area; S3. Set monitoring points in each sub-area to collect and obtain the mechanical state of the flood discharge gate, vibration stress parameters and geological parameters around the gate, and establish a gate stress data set; and extract the mechanical state of the flood discharge gate, vibration stress parameters and geological parameters around the gate to construct the gate wear index Msz. If the gate wear index Msz exceeds the wear threshold value X, send a first warning instruction to the outside; S4. After obtaining and executing the first warning instruction, the gate flood discharge flow-related data is collected through the monitoring points in the sub-area to establish a water flow data set; and the corresponding flow path and the corresponding terminal river data are collected to establish a third data set, and the metadata collected is pre-processed and sent to the database of the centralized control center. After training, the sub-area prediction model is established and trained through neural network technology or time series analysis technology, and the flood discharge flow characteristics and downstream river volume characteristics are extracted. After in-depth analysis, the following are obtained: the average water flow index Slzs and the regional river channel siltation index Ydzs, and the average water flow index Slzs and the regional river channel siltation index Ydzs are associated to calculate and obtain the predicted river channel load index Fzz, and the predicted river channel load index Fzz is compared and analyzed with the load threshold R to obtain the corresponding second warning instruction; S5. Preset the first difference threshold C1 and the second difference threshold C2, and calculate the load difference FCz according to the predicted river load index Fzz and the load threshold R. If the load difference FCz is 0, maintain the current sub-area gate opening value. If the load difference FCz is a positive value, compare and analyze the load difference FCz with the first difference threshold C1 to obtain a first corresponding evaluation result, and generate a corresponding opening instruction based on the first corresponding evaluation result; if the load difference FCz is a negative value, compare and analyze the load difference FCz with the second difference threshold C2 to obtain a second corresponding evaluation result, and generate a corresponding opening instruction based on the second corresponding evaluation result, and set the priority allocation mechanism from large to small according to the size of the load difference FCz.

2. A method for real-time automatic allocation of flood discharge gate opening for a centralized control center according to claim 1, characterized in that: S1 includes: S11. Obtain topographic data of the reservoir area and downstream river channels through GIS or LiDAR scanning, and collect engineering drawings and design parameter data of the reservoir dam, flood discharge gates and related structures; S12. Generate a terrain model based on terrain data and engineering drawings using 3D construction software AutoCAD or 3ds-Max, and use 3ds-Max's "Box", "Cylinder" and "Extrude" tools to create the dam body and gate structure based on engineering drawings; use the "Spline" tool to draw the flood discharge channel path, and use the "Loft" or "Extrude" tool to create a 3D channel model; S13, integrating the terrain model and the three-dimensional flow channel model to obtain a three-dimensional valve centralized control model; S14. Divide the valve centralized control 3D model into several sub-areas according to the positions of gates and flow channels; each sub-area should contain a gate and the corresponding flow channel path and terminal river channel. The valves include spillway gates, bottom outlet flood discharge gates, surface outlet flood discharge gates and power plant flood discharge gates.

3. A method for real-time automatic allocation of flood discharge gate opening for a centralized control center according to claim 1, characterized in that: In the valve centralized control 3D model, several sub-areas are labeled as: "Fm1, Fm2, Fm3, ..., Fm n "n" represents the number of sub-regions.

4. A method for real-time automatic allocation of flood discharge gate opening for a centralized control center according to claim 1, characterized in that: S3 includes: S31, setting up monitoring points in each sub-area to collect and obtain the mechanical state of the flood discharge gate, vibration stress parameters and geological parameters around the gate, and establish a gate stress data set; The gate stress data set includes: current valve opening angle θ, gate operating torque T, gate hydraulic system pressure value P, vibration amplitude A generated during gate operation, gate surface defect area K, gate area geological density M, gate area rock stress N and gate area groundwater level W; S32. Preprocessing the gate stress data set, including data cleaning, data calibration, feature extraction and dimensionless processing.

5. The method for real-time automatic allocation of flood gate opening for a centralized control center according to claim 1, characterized in that: S3 also includes: S33, extracting the mechanical state of the flood discharge gate, the vibration stress parameters and the geological parameter characteristics around the gate from the pre-processed gate stress data set to obtain: a mechanical state index Jxz, a vibration stress index Zdy and a gate area geological index Dzz; The mechanical state index Jxz, vibration stress index Zdy and gate area geological index Dzz are calculated by the following formula: i f =θ*(1+β1*sin(θ)); T f =T*(1+β2*log(T)); P f =P*(1+β3*e -P ); Jxz=γ1*θ f +γ2*T f +γ3*P f ; In the formula, θ f Expressed as the opening angle index, T f Indicates the operating torque index, P f represents the hydraulic system pressure index; θ represents the current valve opening angle, T represents the gate operating torque, P represents the gate hydraulic system pressure value, β1 represents the angle adjustment coefficient, which is used to correct the difference between the actual situation and the ideal state; β2 represents the operating torque adjustment coefficient, β3 represents the gate hydraulic system pressure adjustment coefficient, γ1, γ2 and γ3 are the opening angle index θ respectively. f , operating torque index T f And the hydraulic system pressure index P f The weight coefficient of A f =A**(1+β4*A 2 ); K f =K*(1+β5*e -K ); Jxz=δ1*A f +δ2*K f ; In the formula, A f Represents the vibration amplitude index, K f represents the gate defect index, A represents the vibration amplitude generated during gate operation, K represents the gate surface defect area, β4 represents the vibration amplitude adjustment coefficient, which is used to correct the difference between the actual situation and the ideal state; β5 represents the gate surface defect area adjustment coefficient, δ1 and δ2 are the vibration amplitude index A respectively. f And the gate defect index K f The weight coefficient of M f =M*(1+β6*M); N f =N*(1+β7*N 2 ); W f =W*(1+β8*log(W+1)) Dzz=α1*M f +α2*N f +α3*W f ; Where M f Represents the geological density index, N f Represents rock stress index; W f represents the groundwater level index, M represents the geological density of the gate area, N represents the rock stress of the gate area, W represents the groundwater level in the gate area, β6 represents the geological density adjustment coefficient, which is used to correct the difference between the actual situation and the ideal state; β7 represents the rock stress adjustment coefficient, β8 represents the groundwater level adjustment coefficient, α1, α2 and α3 are the geological density index M respectively. f , rock stress index N f And the groundwater level index W f The weight coefficient; log represents the logarithmic function, and e represents the Sigmoid function; S34, the mechanical state index Jxz, the vibration stress index Zdy and the gate area geological index Dzz are correlated with each other according to the following formula to generate the gate wear index Msz: Msz=α4*Jxz+α5*Zdy+α6*Dzz; Wherein, α4, α5 and α6 are the weight values ​​of the mechanical state index Jxz, the vibration stress index Zdy and the gate area geological index Dzz, respectively, which are adjusted and set by the user, and α4+α5+α6=1; S35, preset a wear threshold X, compare and judge the gate wear index Msz with the wear threshold X, and obtain a first judgment result, including: When the gate wear index Msz>wear threshold X, it means that the gate wear in the sub-area is abnormal, triggering the first warning instruction, including welding, grinding, lubrication and repair of the gate wear parts, detecting the leakage of hydraulic pipelines, repairing or replacing damaged pipelines, and strengthening the geology of the gate area, including foundation grouting and building retaining walls; When the gate wear index Msz ≤ the wear threshold X, it means that the gate wear in the sub-area is normal and continuous monitoring is required.

6. A method for real-time automatic allocation of flood discharge gate opening for a centralized control center according to claim 1, characterized in that: S4 includes: S41, acquiring gate flood discharge flow-related data through monitoring points in the sub-area, and establishing a water flow data set; the water flow data set includes: sub-area water flow velocity Vs, real-time gate discharge flow Q and reservoir real-time water level Sw; S42, dividing the river path corresponding to the sub-area and the corresponding terminal river into several river areas at equal intervals according to the length, and marking them in the valve centralized control three-dimensional model as: "Hd1, Hd2, Hd3, ..., Hd m ” m represents the number of river areas; S43, and collect the corresponding flow path and the corresponding terminal river data to establish a third data set, which includes: regional river water level Xyhd, river cross-sectional area Amj, river sediment concentration Nsnd, river longitudinal slope Spd, river depth Dsd and real-time rainfall JyR.

7. The method for real-time automatic allocation of flood gate opening for a centralized control center according to claim 1, characterized in that: The S4 further comprises: S44, after pre-processing the acquired second data set and third data set and sending them to the centralized control center database, a sub-region prediction model is established, and the flood discharge flow characteristics, downstream river volume characteristics and regional river siltation characteristics are extracted. After in-depth analysis, the following are obtained: average water flow index Slzs and regional river siltation index Ydzs: The average water flow index Slzs and the regional river siltation index Ydzs are calculated by the following formula: In the formula, θ f Expressed as the opening angle index, Vs i is the water flow velocity in the ith river channel area of ​​the sub-area, Amj i is the cross-sectional area of ​​the ith river channel region in the sub-region, Spd i represents the longitudinal slope of the ith river channel area in the sub-area, Q i is the real-time outflow of the ith river area in the sub-area, Sw represents the real-time water level of the reservoir, and m is the number of river areas; Nsnd i is the sediment concentration of the ith river channel area in the sub-area; Dsd i is the depth of the ith channel area in the sub-area; Xyhd i is the water level of the ith river channel area in the sub-area; JyR i is the real-time rainfall in the ith river channel of the sub-region; β9 is the adjustment coefficient of the longitudinal slope of the river channel on the water flow velocity, β 10 is the coefficient that adjusts the effect of water level on sediment deposition; S45, and generating a predicted river load index Fzz by associating the average water flow index Slzs and the regional river siltation index Ydzs through the following formula: Fzz=α7*Slzs+α8*Ydzs; Wherein, α7 and α8 are the weight values ​​of the average water flow index Slzs and the regional river siltation index Ydzs, respectively, which are adjusted and set by the user, and α7+α8=1.

8. A method for real-time automatic allocation of flood gate opening for a centralized control center according to claim 7, characterized in that: Said S4 also includes: S46, preset a load threshold R, and compare and evaluate the predicted river load index Fzz with the load threshold R to obtain a first evaluation result, including: When the predicted river load index Fzz is less than the load threshold R, the comparison value between the current valve opening angle θ and the minimum opening threshold is determined. If the current valve opening angle θ is higher than the minimum opening threshold and the real-time reservoir water level Sw is not lower than the minimum threshold, it indicates that there is a siltation risk in the river of the sub-region, and a second siltation risk warning is sent; If the current valve opening angle θ is higher than the minimum opening threshold, but the real-time reservoir water level Sw is lower than the minimum threshold, it means that there is a drought risk in the reservoir flow, and the second drought risk warning is sent; When the predicted river load index Fzz = load threshold R, it means that the gate opening of the sub-area is normal and the current opening allocation value is maintained; When the predicted river load index Fzz>load threshold R, it means that the gate opening in the sub-area is abnormal, and there is a flood risk in the downstream river flow, and the second flood risk warning is sent.

9. A method for real-time automatic allocation of flood discharge gate opening for a centralized control center according to claim 8, characterized in that: The S5 includes: S51, when receiving the second drought risk warning and the second flood risk warning, the load difference FCz is calculated according to the predicted river load index Fzz and the load threshold R by the following formula: FCz = Fzz-R; S52, presetting a first difference threshold C1, if the load difference FCz is a positive value, comparing and analyzing the load difference FCz with the first difference threshold C1, and generating a first gate opening adjustment strategy, including: When the load difference FCz is less than the first difference threshold C1, the current gate opening angle is adjusted to be reduced by 3-5%; When the first difference threshold C1≤load difference FCz≤first difference threshold C1*130%, the current gate opening angle is adjusted to be reduced by 6-9%; When the load difference FCz> the first difference threshold C1*130%, the current gate opening angle is adjusted to decrease by 10-15%; S53, presetting a second difference threshold C2, if the load difference FCz is a negative value, comparing and analyzing the load difference FCz with the second difference threshold C2, and generating a second gate opening adjustment strategy, including: When the load difference FCz> the second difference threshold C2, the current gate opening angle is adjusted to increase by 3-5%; If the second difference threshold C2 ≤ load difference FCz ≤ second difference threshold C2*130%, the current gate opening angle is adjusted to increase by 6-9%; When the load difference FCz is less than the second difference threshold C2, the current gate opening angle is adjusted to increase by 10-15%.

10. A method for real-time automatic allocation of flood discharge gate opening for a centralized control center according to claim 9, characterized in that: The S5 further includes: S54, according to steps S52 and S53, according to the size of the load difference FCz, a priority processing mechanism is set from large to small, including: The first priority processing sequence is to adjust the current gate opening angle by reducing 10-15% and to increase the current gate opening angle by 10-15%. The second priority processing sequence is to adjust the current gate opening angle by reducing 6-9% and to increase the current gate opening angle by 6-9%. The third priority processing sequence is to adjust the current gate opening angle by reducing 3-5% and to increase the current gate opening angle by 3-5%. and sequentially processing the plurality of sub-regions according to a first priority processing sequence, a second priority processing sequence, and a third priority processing sequence; S55. When receiving the second siltation risk warning, locate and mark the river area where the regional river siltation index Ydzs>siltation threshold YD in the valve centralized control three-dimensional model, and dispatch an engineering team to clean up the silt in real time until the regional river siltation index Ydzs≤siltation threshold YD.

Images