Reservoir evaporation control method integrating evaporation rule analysis and dispatching scene simulation

Abstract

The invention discloses a reservoir evaporation control method integrating evaporation rule analysis and dispatching scene simulation, and relates to the field of reservoir control. A main dam area remote sensing point, a water inlet area remote sensing point and a branched branch shallow bay area remote sensing point serve as representatives to achieve microclimate three-dimensional heterogeneous coverage, and high frequency of water evaporation driving factors can be achieved in combination with a total radiometer, an aerosol laser radar, a layered temperature and humidity probe and other multi-type sensor sets. And carrying out time sequence sampling once in 5 minutes and carrying out spatial distinguishing identification on local disturbance behaviors. Compared with an intermittent measurement mode depending on a single-point water surface evaporation dish or a meteorological station in the prior art, the method has the advantages that the environmental basic data integrity of evaporation behavior prediction is remarkably improved, high-resolution input conditions are provided for subsequent heat flux modeling and feedback control, and the method has the dual technical advantages of precision and real-time performance.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of reservoir control, in particular to a reservoir evaporation control method combining evaporation rule analysis and dispatching scenario simulation. BACKGROUND

[0002] In a large reservoir system in a plateau arid and semi-arid area, there are typical meteorological characteristics such as significant diurnal temperature gradient difference, large fluctuation range of radiant energy, and frequent air disturbance, which make the water body evaporation process present nonlinear, high-frequency disturbance, and nested with complex condensation compensation behavior. As a closed hydrological system, the evaporation behavior of the reservoir is not only affected by thermal driving such as solar radiation and wind disturbance, but also coupled with multiple complex factors such as microclimate stability and aerosol nucleation potential, so the traditional linear scheduling mode cannot reflect the thermal, humid, and dynamic multi-field coupling mechanism in such a system. How to build a closed-loop control method with evaporation potential prediction, condensation compensation identification, and dispatching response feedback as the main axis has become a key challenge in the regulation and control scenarios of plateau lakes and water conservation reservoir groups.

[0003] The current mainstream reservoir dispatching system often performs rigid water distribution or planned regulation and control according to water level lines, inflow, and demand load, and its response to evaporation loss only stays at the level of periodically correcting water balance tables or regressive compensation of empirical loss coefficients, without dynamically identifying the immediate change trend of water body evaporation heat flux, especially lacking modeling and utilization of the night condensation energy compensation behavior. The one-way push defect of this evaporation regulation results in the phenomenon of synchronous amplification of water release and evaporation in sunny high-temperature days in the dispatching system, that is, the maximum dispatching action is implemented in the period with the strongest energy input, while the window with condensation capacity at night is completely not activated, causing the water body heat output peak to lag behind the energy input peak in the daily scale, and failing to build a dynamic water vapor balance mechanism of day release and night compensation. At the same time, the existing models mostly use single-parameter prediction models of wind speed, humidity, and temperature to predict potential evaporation, without effective collection and modeling of key microphysical variables such as disturbance frequency and nucleation particle size, so the accuracy of the regulation and control response cannot penetrate to the bottom logic of the condensation mechanism.

[0004] Due to the lack of control strategy based on the dynamic complementarity principle of evaporation and condensation, the existing reservoir scheduling system is prone to heat flux over-consumption type out-of-control phenomenon, that is, the evaporation heat energy release corresponding to the unit water release behavior is greater than the night energy that can be supplemented by natural condensation, which will lead to secondary effects such as abnormal water level drop, microclimate instability and water body edge drying frequency increase. At the same time, the scheduling response window is fixed and rigid without considering the condensation potential, so that the water release threshold cannot be converged in time during the day with high evaporation, and the energy storage opportunity is lost at night with high potential for recovery, resulting in the double imbalance of "excessive evaporation during the day and waste of cooling at night". Once these systematic deviations accumulate in the high-frequency scheduling period, not only will it cause the mismatch of regional water supply plan, but also it may interfere with the surrounding microclimate structure, causing water vapor disturbance, unstable aerosol inversion path and other feedback abnormalities, and further aggravating the system heat disturbance effect. SUMMARY

[0005] In view of the deficiencies of the prior art, the present application provides a reservoir evaporation control method combining evaporation rule analysis and scheduling scenario simulation, which solves the problems mentioned in the background art.

[0006] To achieve the above purpose, the present application is implemented by the following technical scheme: a reservoir evaporation control method combining evaporation rule analysis and scheduling scenario simulation, comprising the following steps:

[0007] S1, arranging remote sensing microclimate collection points in the reservoir area, and collecting evaporation environment data using a sensor group, and transmitting the evaporation environment data to an edge computing terminal;

[0008] S2, performing time series data regularization and preprocessing on the evaporation environment data in the edge computing terminal to obtain a dimensionless data set, and establishing evaporation heat flux qevp and condensation recovery heat flux qcnd based on the dimensionless data set;

[0009] S3, calculating the day-night water vapor energy difference factor ΔEdn based on the evaporation heat flux qevp and the condensation recovery heat flux qcnd, and then performing preliminary comparative evaluation according to the day-night water vapor energy difference factor ΔEdn to judge the condensation compensation potential;

[0010] S4, triggering the adjustment response mechanism based on the preliminary comparative evaluation result, outputting the responsive scheduling window function Amod based on the day-night water vapor energy difference factor ΔEdn, and adjusting the water release capacity based on the responsive scheduling window function Amod;

[0011] S5, associating and calculating the responsive scheduling window function Amod and the evaporation heat flux qevp to output the behavior energy consumption matching degree index, and performing secondary comparative evaluation based on the behavior energy consumption matching degree index Uctrl to execute the corresponding scheduling strategy.

[0012] Preferably, S1 comprises S11 and S12.

[0013] S11, arranging remote sensing microclimate collection points in the reservoir area, the remote sensing microclimate collection points comprising main dam area remote sensing points, intake area remote sensing points, and branch channel shallow bay area remote sensing points; and configuring a sensor group at each remote sensing microclimate collection point, and setting the sampling frequency of the sensor group to 5 minutes each time to collect evaporation environment data of each remote sensing microclimate collection point in real time;

[0014] The main dam area remote sensing points are arranged in the main dam area of the reservoir.

[0015] The intake area remote sensing points are arranged in the intake area of the reservoir.

[0016] The branch channel shallow bay area remote sensing points are arranged in the branch channel shallow bay area of the reservoir.

[0017] The sensor group comprises a total radiation table, a layered temperature and humidity probe, an aerosol laser radar, a water vapor density sensor, and a temperature and humidity sensor.

[0018] The evaporation environment data comprises a thermal radiation intensity parameter Rsol, a night-time convective gradient factor Xconv, an upper air aerosol condensation potential factor Rgsa, and a water vapor disturbance frequency Vwv.

[0019] S12, transmitting the evaporation environment data collected by each remote sensing microclimate collection point to an edge computing terminal in a LoRa communication protocol through the communication module of the sensor group, and synchronizing the evaporation environment data through a GNSS timestamp.

[0020] Preferably, S2 comprises S21.

[0021] S21, receiving the evaporation environment data in real time at the edge computing terminal, and performing time sequence alignment on all parameters through a sliding window mechanism to obtain a structured data matrix Xsync={x ij}n×m, where x ij represents the jth evaporation environment data at the ith time point, n represents the number of time points, and m represents the number of evaporation environment parameter types.

[0022] Preferably, S2 further comprises S22.

[0023] S22, preprocessing the structured data matrix Xsync, wherein the preprocessing converts the evaporation environment data in each column of the structured data matrix Xsync into standard dimensionless parameters with a standard deviation of 1 and a mean of 0 through a Z-score standardization method for each column of evaporation environment data, to obtain a dimensionless data set.

[0024] Preferably, S2 further comprises S23.

[0025] S23, combining the solar radiation intensity parameter Rsol and the water vapor disturbance frequency parameter Vwv in the dimensionless data set through an evaporation driving factor fitting model to calculate an evaporation heat flux parameter qevp for representing the evaporation latent heat release capacity;

[0026] combining the nighttime convection gradient factor Xconv and the high-altitude aerosol condensation potential energy factor Rgsa with the water vapor disturbance frequency parameter Vwv in the dimensionless data set through a condensation nucleus potential modeling process to construct a condensation re-supply heat flux parameter qcnd for representing the nighttime water vapor condensation behavior.

[0027] Preferably, the S3 includes S31.

[0028] S31, based on the time period integration operation of the evaporation heat flux parameter qevp and the condensation re-supply heat flux parameter qcnd, integrating the unit area heat released by the evaporation process and the unit area heat absorbed by the condensation process in the respective time period to obtain a day-night water vapor energy difference factor ΔEdn, which is calculated by subtracting the integral value of the condensation heat flux in the nighttime period tn from the integral value of the evaporation heat flux in the daytime period td.

[0029] Preferably, the S3 further includes S32.

[0030] S32, a condensation compensation potential reference threshold Eth is preset, which is based on the maximum recoverable energy under the optimal condensation capacity condition in the historical monitoring data as the boundary.

[0031] After the condensation compensation potential reference threshold Eth and the day-night water vapor energy difference factor ΔEdn are preliminarily compared and evaluated, the re-supply potential of nighttime condensation and the condensation compensation potential of evaporation loss in the reservoir water system are judged, and the specific evaluation content is as follows.

[0032] When the day-night water vapor energy difference factor ΔEdn is greater than or equal to the condensation compensation potential reference threshold Eth, it indicates that the condensation is insufficient and the evaporation is abnormal, at which time the adjustment response mechanism is triggered immediately.

[0033] When the day-night water vapor energy difference factor ΔEdn is less than the condensation compensation potential reference threshold Eth, it indicates that the day-night energy balance is sufficient, and no intervention is needed at this time.

[0034] Preferably, the S4 includes S41.

[0035] S41. After preliminary comparative evaluation of the trigger regulation response mechanism, based on the numerical results of the diurnal water vapor energy difference factor ΔEdn and combined with the principle of thermal balance regulation, a type of scheduling window function mapping model is constructed. This type of scheduling window function mapping model uses the diurnal water vapor energy difference factor ΔEdn as the input variable and the regulation time window width as the output target to dynamically generate a responsive scheduling window function Amod. The responsive scheduling window function Amod is the regulation capacity function. The responsive scheduling window function Amod is passed as the control input parameter to the variable water release capacity regulation module for preliminary actual water level adjustment.

[0036] Preferably, S5 includes S51;

[0037] S51. After the initial actual water level adjustment, the variable water release capacity control module of the reservoir scheduling system extracts the energy consumption per unit area at time t and the actual water release flow at time t corresponding to the actual scheduling behavior, qactual(t). Within the set observation period T, the ratio of the evaporation heat flux parameter qevp(t) at time t to the energy consumption per unit area at time t corresponding to the actual scheduling behavior is calculated to form an instantaneous matching index. Combined with the weight of the ratio of the responsive scheduling window function Amod(t) at the current time point t to the actual water release flow at time t, a combined kernel function is formed. Then, the combined kernel function is processed by integral averaging within the time interval [0,T] to output the behavior energy consumption matching degree index Uctrl.

[0038] Preferably, S5 further includes S52;

[0039] S52. Based on the behavioral energy consumption matching index Uctrl, a secondary comparative evaluation is performed to determine the initial control status of the actual water level adjustment. Based on the results of the secondary comparative evaluation, the corresponding scheduling strategy is executed. The specific evaluation content is as follows:

[0040] When the behavior energy consumption matching index Uctrl < 0.85, it indicates that the control state is ideal, the behavior is reasonable, and the current scheduling strategy should be maintained.

[0041] When 0.85 ≤ behavior energy consumption matching degree index Uctrl < 1, it indicates that the control state is critically deviating and there is a slight deviation. At this time, the tightening strategy is executed, and the responsive scheduling window function Amod(t) at time t is iteratively regenerated until the control state is ideal and the iteration stops.

[0042] When the behavior energy consumption matching index Uctrl≥1, the expression control state is severely over-consumed and the strategy is out of control. At this time, forced compression is started, compressing the responsive scheduling window function Amod(t) at time t to no more than 35%, freezing 12:00-16:00, and opening window permissions.

[0043] The application provides a reservoir evaporation control method combining evaporation rule analysis and scheduling scenario simulation. The method has the following beneficial effects:

[0044] (1) The method can realize high-frequency, 5-minute time-series sampling and spatial resolution identification of local disturbance behavior of water body evaporation driving factors by deploying a multi-zone remote sensing microclimate collection point system in the reservoir area, taking the main dam area remote sensing point, the water inlet area remote sensing point and the branch bay area remote sensing point as representatives to realize three-dimensional heterogeneous coverage of microclimate, and combining a multi-type sensor group such as a total radiation table, an aerosol laser radar, and a layered temperature and humidity probe. Compared with the intermittent measurement method of the prior art which relies on a single-point water surface evaporation pan or a meteorological station, the scheme has the advantages of high throughput, low lag, wide regional distribution, and various parameter configurations, significantly improving the environmental basic data integrity of evaporation behavior prediction, providing high-resolution input conditions for subsequent heat flux modeling and feedback control, and having dual technical advantages of precision and real-time.

[0045] (2) The method introduces an evaporation heat flux parameter qevp constructed based on the solar radiation intensity parameter Rsol and the water vapor disturbance frequency parameter Vwv, and at the same time, a condensation back-supply heat flux parameter qcnd constructed by combining the night-time convection gradient factor Xconv and the high-altitude aerosol condensation potential energy factor Rgsa, generates a day-night water vapor energy difference factor ΔEdn through time period integral analysis, and for the first time, incorporates evaporation released latent heat and condensation absorbed latent heat into a unified heat energy conservation calculation model, breaking through the problem in the prior art that only considers energy dissipation from the evaporation side. The method establishes a new modeling framework of "heat, humidity, and disturbance coupling driving + day-night asymmetric mechanism analysis + disturbance suppression term deduction", which not only improves the scientific accuracy of water body energy loss estimation, but also provides a decision-making intermediate variable with a physical basis for dynamic scheduling, filling the modeling gap of lacking quantification of condensation contribution in current regulation models.

[0046] (3) The method realizes a two-stage comparative evaluation mechanism of scheduling behavior by constructing a responsive scheduling window function Amod and a behavior energy matching index Uctrl: first, a type of scheduling window function mapping model is used to convert the day-night water vapor energy difference factor AEdn into a dynamic controllable scheduling capability input, forming a closed-loop response path of "prediction → water release capability convergence"; second, the behavior energy matching index Uctrl is used to comprehensively depict the "scheduling behavior energy matching" and "behavior strategy capacity matching", so as to realize the self-iteration and saturation protection mechanism of the fine adjustment strategy. Compared with the rigid threshold water release scheduling mode in the prior art, the method can automatically compress the high evaporation period window according to the dynamic change of water vapor heat flux, or freeze the release authority in the key period, so as to realize the prospective suppression of abnormal evaporation scene and the strategic release of condensation back-up potential, has strong dynamic adaptability and control precision advantage, and significantly improves the regulation stability and water saving ability of the reservoir in the extreme weather environment BRIEF DESCRIPTION OF DRAWINGS

[0047] Figure 1 A step schematic diagram of a reservoir evaporation control method fusing evaporation rule analysis and scheduling scenario simulation according to the present application;

[0048] Figure 2 High evaporation period window freeze control diagram;

[0049] Figure 3 Daytime evaporation heat flux qevp and nighttime condensation back-up heat flux qcnd integral comparison diagram. DETAILED DESCRIPTION

[0050] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor are within the scope of protection of the present application.

[0051] Embodiment 1

[0052] The present application provides a reservoir evaporation control method fusing evaporation rule analysis and scheduling scenario simulation, please refer to Figure 1 , comprising the following steps:

[0053] S1, arranging remote sensing microclimate collection points in the reservoir area, and collecting evaporation environment data by using a sensor group, and transmitting the evaporation environment data to an edge computing terminal;

[0054] S2, time series data regularization and preprocessing of evaporation environment data in the edge computing terminal to obtain a dimensionless data set, and establishing evaporation heat flux qevp and condensation heat flux qcnd based on the dimensionless data set;

[0055] S3, calculating the day-night water vapor energy difference factor ΔEdn based on the evaporation heat flux qevp and the condensation heat flux qcnd, and preliminarily comparing and evaluating according to the day-night water vapor energy difference factor ΔEdn to judge the condensation compensation potential;

[0056] S4, triggering the adjustment response mechanism based on the preliminary comparison and evaluation result, outputting the responsive scheduling window function Amod based on the day-night water vapor energy difference factor ΔEdn, and adjusting the water release capacity based on the responsive scheduling window function Amod;

[0057] S5, associating and calculating the responsive scheduling window function Amod with the evaporation heat flux qevp to output the behavior energy consumption matching degree index, and performing secondary comparison and evaluation based on the behavior energy consumption matching degree index Uctrl to execute the corresponding scheduling strategy.

[0058] In this embodiment, the method can cover key hydrological heterogeneous areas such as the main dam, the water inlet and the branch bay in space and realize minute-level evaporation environment data sampling in time through the multi-point remote sensing microclimate collection layout of step S1 and the high-frequency sensor group. Then, the dimensionless processing and parameter fusion modeling are completed by the edge computing terminal in step S2, and the evaporation heat flux parameter qevp and the condensation heat flux parameter qcnd that can dynamically respond to the perturbed environment are constructed. Subsequently, in step S3, the diurnal water vapor energy difference factor ΔEdn is generated by quantitatively evaluating the diurnal energy difference, the dynamic balance relationship between condensation potential and evaporation loss is judged, and the evaporation anomaly recognition from the perspective of "energy conservation" is realized. In step S4, the responsive scheduling window function Amod is generated based on ΔEdn, and the mapping bridge between the diurnal energy deviation and the actual water level adjustment is constructed by controlling the scheduling ability to compress or relax according to the energy deviation. Finally, in step S5, the dynamic matching evaluation between the discharge behavior and the evaporation demand is realized by constructing the behavior energy consumption matching degree index Uctrl, and the scheduling strategy iteration optimization is completed through the Uctrl secondary discrimination mechanism, and the complete control logic chain from identification, adjustment to closed loop rebalancing is established. The method not only breaks through the traditional evaporation control idea of "one-way passive inhibition", but also introduces the night condensation potential as an active scheduling factor for the first time, and realizes the fine quantitative modeling of water energy flow in the heat flux space, effectively improving the early warning ability and adjustment accuracy of extreme evaporation behavior. Finally, the method realizes three aspects of significant improvement. In terms of prediction accuracy, the regional adaptability of evaporation evaluation is significantly improved through the microclimate distributed sensing network and the perturbation sensitive modeling. In terms of regulation ability, the abnormal energy consumption behavior is pre-emptively constrained by dynamically compressing the scheduling window through the responsive scheduling window function Amod. In terms of strategy closed loop, the precise alignment between scheduling behavior and energy demand is realized through the feedback mechanism driven by the behavior energy consumption matching degree index Uctrl. In summary, the technical scheme of the present application not only solves the core problems of existing reservoir evaporation regulation model response lag, condensation potential non-dimensional modeling, and scheduling behavior and energy consumption disconnection, but also significantly improves the system adaptability, scheduling water saving and model material interpretability, and has significant technical progress value and application potential.

[0059] Embodiment 2

[0060] Please refer to Figure 1 , specifically: S1 includes S11 and S12;

[0061] S11. Remote sensing microclimate collection points are set up in the reservoir area, including remote sensing points in the main dam area, remote sensing points in the intake area, and remote sensing points in the shallow bay area of ​​the tributaries; and a sensor group is configured at each remote sensing microclimate collection point, with the sampling frequency of the sensor group set to 5 minutes each time, to collect the evaporation environment data of each remote sensing microclimate collection point in real time.

[0062] Remote sensing points in the main dam area are set up in the main dam area of ​​the reservoir, where the water body heat capacity response is the strongest, the wind field is stable, and the light reflection on the water surface is the most concentrated, reflecting the main trend of overall evaporation of the reservoir.

[0063] Remote sensing points in the intake area are set up in the reservoir's intake area to detect common water disturbances and disturbances to the atmosphere structure caused by water diversion paths, and to identify the impact of water diversion activities on local evaporation disturbances.

[0064] Remote sensing points in the tributary shallow bay area are set up in the tributary shallow bay area of ​​the reservoir. The water is shallow and has low thermal inertia, making it easy for condensation core zones to form, thus capturing microscale nighttime condensation events.

[0065] The sensor group includes a total radiation meter, a stratified temperature and humidity probe, an aerosol lidar, a water vapor density sensor, and a temperature and humidity sensor;

[0066] Evaporation environment data include thermal radiation intensity parameter Rsol, nighttime convection gradient factor Xconv, upper-air aerosol condensation potential energy factor Rgsa, and water vapor disturbance frequency Vwv;

[0067] The thermal radiation intensity parameter Rsol is calculated as the energy flux per unit area by measuring the temperature rise of the radiant heat on the thermistor element of the total radiation meter installed at the unobstructed main dam area remote sensing point, including direct and diffuse light.

[0068] The nighttime convection gradient factor Xconv is obtained by using stratified temperature and humidity probes at three heights (0.5m, 2m, and 4m above the water surface) to collect temperature and relative humidity in real time and calculate the gradient. The larger the gradient, the stronger the nighttime air convection and the easier it is to produce condensation. The superposition of convective heat transfer and humidity migration describes the possibility of cold air sinking and humidified air rising.

[0069] The upper-level aerosol condensation potential energy factor Rgsa was determined by using aerosol lidar to observe the concentrations of aerosols in the atmosphere (PM0.3-PM2.5) and their average particle size. Based on the equation meteorological physics, the concentrations of atmospheric aerosols (PM0.3-PM2.5) and average particle size are fitted to obtain the expression, which characterizes the ability of fine particulate matter in the atmosphere, such as salt particles and dust, to provide condensation nuclei at night and is a triggering factor for non-explicit condensation.

[0070] The water vapor disturbance frequency Vwv is obtained by continuously recording the water vapor content in the air using a water vapor density sensor. The number of second derivative changes is analyzed within a fixed time window, and a disturbance judgment threshold is set. The number of times the disturbance judgment threshold is exceeded is considered a sudden change in water vapor concentration, and frequency analysis is performed. Frequent disturbances can disrupt the balance of the water vapor layer and induce abnormal evaporation; they can also cause the condensation layer to be broken down.

[0071] S12. The evaporation environment data collected by each remote sensing microclimate collection point is transmitted to the edge computing terminal via the LoRa communication protocol through the communication module deployed in the sensor group, and the evaporation environment data is synchronized through GNSS timestamp to ensure the temporal consistency and regional resolution effectiveness of the evaporation environment data.

[0072] In this embodiment, the method achieves precise coverage of water body evaporation heterogeneity, disturbance, and condensation-sensitive areas by scientifically deploying remote sensing microclimate acquisition points in the main dam area, the intake area, and the shallow bay areas of tributaries in S11. The functional differences in the acquisition points in each area, such as capturing the overall evaporation trend of the reservoir in the main dam area, identifying the impact of water diversion disturbances at the intake, and detecting nighttime condensation formation conditions in the shallow bay areas of tributaries, ensure a comprehensive understanding of the spatial heterogeneity of evaporation behavior. Combined with a multi-source sensor group consisting of a total radiation meter, stratified temperature and humidity probes, aerosol lidar, and water vapor density sensor, the system can acquire key evaporation control factors such as thermal radiation intensity parameter Rsol, nighttime convection gradient factor Xconv, upper-level aerosol condensation potential energy factor Rgsa, and water vapor disturbance frequency Vwv at high frequency, realizing a quantitative characterization of the core driving mechanism of the evaporation process. In S12, by uploading data using the LoRa communication protocol through the communication modules configured at each acquisition point, and combining this with GNSS timestamps for cross-point data timing alignment, the time drift and spatial distortion problems under multi-point concurrent acquisition are effectively solved, ensuring the spatiotemporal consistency and processing reliability of evaporation environment data. This implementation scheme achieves the following: improved environmental perception coverage: by deploying a structured acquisition point network, taking into account the spatial diversity of reservoir evaporation behavior and microclimate interference characteristics; improved physical targeting of environmental parameter perception: by selectively using sensor parameters with physical interpretability, such as Rgsa and Xconv, multi-dimensional capture of condensation conditions and evaporation intensity is achieved; enhanced data transmission stability and timing consistency: through LoRa low-power wide-area network and GNSS synchronization mechanism, continuous, reliable, and low-latency transmission capabilities are ensured over a large area of ​​water. Overall, step S1 constructs a solid data foundation and environmental perception front-end in the entire reservoir evaporation control system, providing highly reliable and spatially resolved environmental information support for subsequent flux modeling, scheduling simulation, and feedback evaluation, greatly improving the responsiveness, accuracy, and predictability of reservoir fine-tuning.

[0073] Example 3

[0074] Please see Figure 1 Specifically: S2 includes S21;

[0075] S21. Receive evaporation environment data in real time at the edge computing terminal, and perform time-series alignment of all parameters using a sliding window mechanism to obtain the structured data matrix Xsync = {x ij}n×m, where x ij This represents the evaporation environment data at time point i, n represents the number of time points, i.e. the length of the time series obtained by continuous sampling. For example, if data is collected every 5 minutes, there will be 288 time points in 24 hours. m represents the number of types of evaporation environment parameters, i.e. the evaporation environment data collected at each time point, such as the four core parameters collected in this method.

[0076] S2 also includes S22;

[0077] S22. Preprocess the structured data matrix Xsync. The preprocessing method is to use the Z-score standardization method to transform the evaporation environment data of each column in the structured data matrix Xsync into a standard dimensionless parameter with a standard deviation of 1 and a mean of 0, thereby obtaining a dimensionless data set.

[0078] S2 also includes S23;

[0079] S23. Combine the solar radiation intensity parameter Rsol and the water vapor disturbance frequency parameter Vwv in the dimensionless dataset through the evaporation driving factor fitting model to construct the evaporation heat flux parameter qevp to characterize the evaporation latent heat release capacity. The evaporation driving factor fitting model is based on the variable simplification of the classic Penman evaporation model. The total solar radiation intensity parameter Rsol is used to replace the net radiation index to optimize the physical interpretation under clear sky conditions. At the same time, the water vapor disturbance frequency parameter Vwv is used to replace the traditional wind speed function to optimize the response sensitivity to evaporation behavior in the disturbed environment.

[0080] The net radiation factor in the traditional model is replaced by the solar radiation intensity parameter Rsol. The reason is that in cloudless or partly cloudy weather, the actual net radiation of the earth's surface is mainly determined by solar radiation. Therefore, the two are highly correlated, and the replacement simplifies the requirements for measurement equipment.

[0081] The wind speed function is replaced by the water vapor disturbance frequency parameter Vwv to reflect the reinforcing effect of frequent disturbances on evaporation behavior in areas such as plateaus, wind gaps, or canyons. The disturbance frequency, as a micro-disturbance excitation index, has a higher physical interpretation capability.

[0082] The nighttime convection gradient factor Xconv, the upper-level aerosol condensation potential energy factor Rgsa, and the water vapor disturbance frequency parameter Vwv from the dimensionless dataset are combined using a condensation nucleation potential modeling process to construct the condensation return heat flux parameter qcnd, which characterizes nighttime water vapor condensation behavior. The condensation nucleation potential modeling process is based on... Based on the theory and principles of environmental disturbance, this paper considers the condensation nucleation capacity and atmospheric vertical stability as positive contributing factors, and introduces disturbance frequency as a nucleation inhibition term to form a comprehensive modeling method for nighttime latent heat release behavior. By considering the temperature and humidity gradients in the vertical direction at night, the stability of the condensation environment is reflected. At the same time, the potential of fine particles and aerosols in the atmosphere as condensation nuclei at night is combined as triggering factors for the possibility of condensation. Furthermore, the adverse effects of water vapor disturbance on nucleation stability are eliminated, thus constructing a condensation heat flux model that can simultaneously reflect nucleation conditions, condensation driving forces, and disturbance disturbances.

[0083] Throughout the modeling process, all parameters are normalized via edge computing terminals before use to ensure their dimensionlessness, thus ensuring comparability and stability in the modeling process. When each physical quantity is substituted into the model, the unit system remains consistent. The evaporation rate and condensation rate are based on the liquid surface loss / replenishment thickness per unit length per unit time. The flux output is uniformly represented as the energy value per unit area per unit time, ensuring the consistency and interpretability of the model output.

[0084] In this embodiment, the method introduces a sliding window mechanism in the edge computing terminal to perform temporal alignment of the original evaporation environment data, constructing a structured data matrix Xsync. This ensures that each parameter has a consistent temporal semantic relationship at different time points, improving the temporal accuracy of data processing and the robustness of model input. S22 further uses the Z-score normalization method to perform dimensionless processing on each column of parameter data, eliminating parameter scale differences and constructing a dimensionless data set under a unified dimensional system, laying a mathematical foundation for subsequent multi-source parameter fusion modeling. In S23, the solar radiation intensity parameter Rsol and the water vapor disturbance frequency parameter Vwv from the dimensionless data set are innovatively introduced into the classic Penman evaporation model framework to construct an evaporation driving factor fitting model for calculating the evaporation heat flux parameter qevp. This effectively simplifies the dependence on factors that are difficult to measure precisely, such as net radiation and wind speed, improving the physical adaptability of evaporation behavior modeling under clear skies and disturbed environments. Simultaneously, it employs... Based on theory and disturbance principles, a condensation nucleus potential energy model was established by combining the nighttime convection gradient factor Xconv, the upper-air aerosol condensation potential energy factor Rgsa, and the water vapor disturbance frequency parameter Vwv. This model was used to calculate the condensation recovery heat flux parameter qcnd, comprehensively characterizing the nighttime condensation heat release behavior from three perspectives: nucleation, driving force, and disturbance. Both flux models were normalized to ensure that the physical quantity output is the energy value per unit area per unit time (J / m²). 2 The S2 step ( / s) enhances the model's versatility and interpretability across time and location. It achieves high-precision temporal alignment and standardized normalization of multi-source evaporation environment data, ensuring data modeling consistency and dimensional uniformity. A heat flux modeling mechanism with physical foundations and engineering simplification is constructed, significantly improving the availability and prediction accuracy of the two key energy factors, qevp and qcnd. It overcomes the limitations of traditional models' dependence on meteorological stations, enhancing sensitivity and discriminative ability to local microclimate disturbances.

[0085] Example 4

[0086] Please see Figure 1 and Figure 3 Specifically: S3 includes S31;

[0087] S31. The time-span integration operation based on the evaporation heat flux parameter qevp and the condensation heat flux parameter qcnd is constructed based on the basic principles of energy conservation and heat flux accumulation over time in physics. The heat released per unit area during the evaporation process and the heat absorbed per unit area during the condensation process are integrated and compared in their respective time periods to obtain the diurnal water vapor energy difference factor ΔEdn. ​​The diurnal water vapor energy difference factor ΔEdn is calculated by subtracting the integral value of the condensation heat flux in the nighttime time period tn from the integral value of the evaporation heat flux in the daytime time period td.

[0088] The specific calculation formula is: △Edn=(∫tdqevp(t)dt)-(∫tnqcnd(t)dt), where td represents the daytime period, tn represents the nighttime period, qebp(t) represents the evaporation heat flux parameter at time t, qcnd(t) represents the condensation compensation heat flux parameter at time t, and dt represents the time integral function;

[0089] (∫tdqevp(t)dt) represents the cumulative value of the evaporative heat flux per unit area qevp during the daytime period. Its physical meaning is that it represents the total latent heat absorbed by water bodies during the daytime evaporation process, which converts liquid water into water vapor.

[0090] (∫tnqcnd(t)dt) represents the cumulative value of the condensation return heat flux parameter qcnd per unit area during the nighttime phase. Its physical meaning is that it represents the total latent heat released during the nighttime environment when water vapor is condensed into droplets.

[0091] Specifically, the units of the evaporation heat flux parameter qevp and the condensation compensation heat flux parameter qcnd are both joules per square meter per second (J / m²). 2 / s), after time integration, is converted into the cumulative energy per unit area during the day and night (J / m²). 2 To ensure dimensional consistency during calculations, the upper and lower limits of integration are determined based on the local "daytime" and "nighttime" durations as defined in the environmental monitoring system, typically using local sunrise and sunset times for separation.

[0092] In this embodiment, the evaporation heat flux parameter qevp and the condensation compensation heat flux parameter qcnd can be numerically integrated using the sliding time window method. By using heat flux data with a resolution of five minutes or ten minutes, a fixed-width rectangular integral is performed on the daytime and nighttime time periods to obtain the cumulative value of evaporation energy and the cumulative value of condensation energy, and then the value of the diurnal water vapor energy difference factor ΔEdn is calculated.

[0093] S3 also includes S32;

[0094] S32. Preset condensation compensation potential reference threshold Eth. The condensation compensation potential reference threshold Eth is based on the maximum recoverable energy under the optimal condensation capacity condition in historical monitoring data as the boundary. It is used to evaluate whether the replenishment potential of nighttime condensation in the current water system can meet the evaporation loss.

[0095] A preliminary comparative assessment was conducted by comparing the condensation compensation potential reference threshold Eth with the diurnal water vapor energy difference factor ΔEdn to determine the nighttime condensation replenishment potential and evaporation loss condensation compensation potential in the reservoir water system. The specific assessment content is as follows:

[0096] When the diurnal water vapor energy difference factor ΔEdn is greater than or equal to the condensation compensation potential reference threshold Eth, it indicates insufficient condensation and abnormal evaporation, and the regulation response mechanism is immediately triggered.

[0097] When the diurnal water vapor energy difference factor ΔEdn is less than the condensation compensation potential reference threshold Eth, it indicates that the diurnal energy is balanced and condensation is sufficient, and no intervention is needed.

[0098] In this embodiment, in S31, the method employs the principles of heat flux conservation and time integration from physics to integrate the daytime evaporation heat flux parameter qevp and the nighttime condensation replenishment heat flux parameter qcnd over time periods, calculating the difference in total latent heat release per unit area between day and night, and generating a day-night water vapor energy difference factor ΔEdn. ​​This factor not only reflects the balance between evaporation and condensation dynamics under current climate and hydrological conditions but also provides a clear basis for quantitative decision-making in subsequent scheduling and regulation. In particular, the use of sliding time window integration and high-frequency data input, such as a 5-minute resolution, gives the day-night water vapor energy difference factor ΔEdn good timeliness and trend capture capabilities, significantly improving its sensitivity to abnormal evaporation or insufficient condensation. In S32, a dynamic heuristic evaluation system is constructed by constructing a condensation compensation potential reference threshold Eth and comparing it in real time with the day-night water vapor energy difference factor ΔEdn. ​​This reference threshold Eth is derived from the energy recovery capacity under optimal condensation conditions in historical monitoring data, possessing practical rationality in an engineering context. Compared to traditional evaporation models that only provide static calculation results, this scheme can make adjustment decisions based on whether the diurnal water vapor energy difference factor ΔEdn exceeds the condensation compensation potential reference threshold Eth, achieving real-time judgment on whether "condensation is sufficient to compensate for evaporation". When the diurnal water vapor energy difference factor ΔEdn exceeds the threshold, the system can immediately trigger the adjustment response mechanism; while under energy balance conditions, excessive intervention is avoided, reflecting the ideas of differentiated control and energy consumption optimization. Therefore, the purpose and beneficial effects achieved by this implementation method are mainly reflected in the following aspects: based on the principle of energy conservation, a two-way quantitative model for the latent heat behavior of evaporation and condensation is constructed, forming a complete energy budget framework; the dynamic accuracy and real-time performance of the diurnal water vapor energy difference factor ΔEdn calculation are significantly improved through the sliding integral method, enabling the evaluation mechanism to have immediate judgment capabilities; and the introduction of the condensation compensation potential reference threshold Eth as a dynamic trigger threshold is used to construct an active judgment mechanism, avoiding misjudgments in traditional passive scheduling.

[0099] Example 5

[0100] Please see Figure 1 Specifically: S4 includes S41;

[0101] S41. After preliminary comparative evaluation of the trigger regulation response mechanism, based on the numerical results of the diurnal water vapor energy difference factor ΔEdn and combined with the principle of heat balance regulation, a type of scheduling window function mapping model is constructed. This type of scheduling window function mapping model uses the diurnal water vapor energy difference factor ΔEdn as the input variable and the regulation time window width as the output target to dynamically generate a responsive scheduling window function Amod. The responsive scheduling window function Amod is a regulation capacity function that reflects the degree to which the water release window should be expanded. The responsive scheduling window function Amod is then passed as a control input parameter to the variable water release capacity regulation module for preliminary actual water level adjustment.

[0102] The responsive scheduling window function Amod calculates its output through the following type of scheduling window function mapping model;

[0103]

[0104] In the formula, Amod(t) represents the responsive scheduling window function at time t, Qref(t) represents the daily reference water supply at time t, Eref represents the upper limit of historical condensation replenishment energy, used to normalize the day-night water vapor energy difference factor ΔEdn, and α represents the response sensitivity coefficient, which controls the compression intensity and is set by the user according to experience to be between 0.5 and 1.0.

[0105] The design principle of the responsive scheduling window function Amod draws on environmental feedback scheduling theory and proportional weakening feedback mechanism in control engineering. It performs dynamic compression based on the default reference flow. The compression factor is dynamically adjusted according to the ratio of the diurnal water vapor energy difference factor ΔEdn to the reference energy value. The larger the diurnal water vapor energy difference factor ΔEdn is, the stronger the compression. A first-order linear compression function is used, which is convenient for algorithm implementation and has a simple physical interpretation. This algorithm belongs to the nonlinear scaling first-order weakening model and is widely used in energy balance control and flow control strategies.

[0106] One type of scheduling window function mapping model originates from the Logistic suppression nonlinear mapping model in mathematical physics and is widely used in control modeling of ecological, hydrological, and energy systems with characteristics of "lag, response, and saturation".

[0107] △Edn

[0108] Since Eref is dimensionless, the entire bracket is also dimensionless. After multiplying it by the daily reference water supply Qref(t) at time t, the output responsive scheduling window function Amod(t) at time t still maintains m. 3 / s, so the units on both sides of the formula are consistent, which is logically and physically reasonable and has a rigorous dimension.

[0109] In this embodiment, the method, through step S4 and sub-step S41, constructs and outputs a responsive scheduling window function Amod based on the real-time analysis results of the diurnal water vapor energy difference factor ΔEdn, achieving precise linkage between energy difference perception and scheduling behavior control. Specifically, this implementation method, based on the principle of thermal balance regulation, uses the diurnal water vapor energy difference factor ΔEdn as the adjustment input variable, introduces a type of scheduling window function mapping model, and combines it with the historical condensation replenishment energy upper limit value Eref to complete the normalization process, dynamically outputting a window adjustment capability that matches the degree of coupling with the diurnal water vapor energy difference factor ΔEdn. ​​In the implementation of this method, the responsive scheduling window function Amod is not only constructed as a water release capacity control function, but also controls its compression intensity by setting the response sensitivity coefficient α to an empirical range of 0.5 to 1.0, thereby achieving a sensitive response to changes in the diurnal water vapor energy difference factor ΔEdn. In terms of functional form, this model borrows from the Logistic suppression nonlinear mapping mechanism, which is widely used in behavioral modeling with "hysteresis, response, and saturation" characteristics in fields such as ecosystems and hydrological energy control, and has good physical interpretability and dynamic regulation performance. Through this function, the reservoir scheduling system can automatically determine the intensity of the current energy imbalance and dynamically compress the reference water supply Qref(t) to generate the responsive scheduling window function Amod(t) at time t, thereby realizing real-time water level regulation behavior. Therefore, this method establishes a functional bridge from energy difference determination to water body scheduling capacity compression, bridging the technical gap between previous evaporation compensation assessment and water release scheduling control. By employing a first-order linear compression algorithm incorporating the proportional weakening feedback principle, the control response process is both simple and feasible, and possesses good engineering adaptability. While maintaining the original physical dimensions of Qref(t), it achieves dimensional consistency and energy consistency in the output of the responsive scheduling window function Amod(t) at time t, ensuring the physical rigor of the model. Through dynamic mapping response to the diurnal water vapor energy difference factor ΔEdn, a synchronization mechanism between evaporation anomalies and control capacity compression is constructed, significantly improving the water resource system's adaptability to extreme climates and short-term fluctuations. Overall, this step effectively improves the dynamic accuracy and energy efficiency of reservoir scheduling, providing a central control module support for the entire evaporation control strategy system.

[0110] Example 6

[0111] Please see Figure 1 and Figure 2 Specifically: S5 includes S51;

[0112] S51. After the initial actual water level adjustment, the variable water release capacity control module of the reservoir scheduling system extracts the energy consumption per unit area at time t and the actual water release flow at time t corresponding to the actual scheduling behavior. Within the set observation period T, the ratio of the evaporation heat flux parameter qevp(t) at time t to the energy consumption per unit area at time t corresponding to the actual scheduling behavior is calculated to form an instantaneous matching index. Combined with the ratio weight of the responsive scheduling window function Amod(t) at the current time point t and the actual water release flow at time t, a combined kernel function is formed. Then, the combined kernel function is integrated and averaged over the time interval [0,T] to output the behavior energy consumption matching degree index Uctrl, which is used to characterize the dynamic fit between scheduling behavior and energy consumption demand.

[0113] The behavior energy consumption matching index Uctrl is calculated and output using the following algorithm formula;

[0114]

[0115] In the formula, dt represents the time integral function;

[0116] This represents an instantaneous matching indicator, used to analyze whether behavior is excessive or lacking in energy release;

[0117] This represents the ratio of the compressed scheduling window capacity to the actual behavior, used to analyze whether scheduling closely follows the policy capacity.

[0118] Since both indicators are fractions, both sides are dimensionless parameters. The first ratio reflects behavior-energy consumption matching; the second ratio reflects behavior-strategy capacity matching; their product forms a composite matching factor.

[0119] For a specific process example, let qactual(t) = 60.27 J / m 2 / s; Qactual(t) = 40m 3 / s, Amod(t)=48m 3 / s, qevp(t)=85J / m 2 ;

[0120] Single-point integral function value The average integral over all points within one hour yields a behavior energy consumption matching index Uctrl = 0.59.

[0121] S5 also includes S52;

[0122] S52. Based on the behavioral energy consumption matching index Uctrl, a secondary comparative evaluation is performed to determine the initial control status of the actual water level adjustment. Based on the results of the secondary comparative evaluation, the corresponding scheduling strategy is executed. The specific evaluation content is as follows:

[0123] When the behavior energy consumption matching index Uctrl < 0.85, it indicates that the control state is ideal, the behavior is reasonable, and the current scheduling strategy should be maintained.

[0124] When 0.85 ≤ behavior energy consumption matching degree index Uctrl < 1, it indicates that the control state is critically deviating and there is a slight deviation. At this time, the tightening strategy is executed, and the responsive scheduling window function Amod(t) at time t is iteratively regenerated until the control state is ideal and the iteration stops.

[0125] The tightening strategy reduces the upper limit of the release function of the responsive scheduling window function Amod(t) at time t by tightening the response sensitivity coefficient α by 10%.

[0126] The response sensitivity coefficient α determines the response strength of the responsive scheduling window function Amod to the change of the diurnal water vapor energy difference factor ΔEdn, which is reflected as the scaling factor between the diurnal water vapor energy difference factor ΔEdn and the responsive scheduling window function Amod.

[0127] When the response sensitivity coefficient α is large, even slight changes in the diurnal water vapor energy difference factor ΔEdn will be amplified into a large change in the responsive scheduling window function Amod, resulting in excessive adjustment of the system's water release range.

[0128] When the behavior energy consumption matching index Uctrl enters the 0.85–1 range, it means that the system is “approaching overconsumption” or “behavior is following too fast” energy consumption signal, and there has been an excessive energy release or a lag in regulation.

[0129] If the response sensitivity coefficient α is not adjusted at this point, the response will still be amplified in the next strategy iteration, which may lead to further deviation and cause "strategy positive feedback to go out of control".

[0130] When the behavior energy consumption matching degree index Uctrl≥1, the expression control state is severely over-consumed and the strategy is out of control. At this time, forced compression is started, compressing the responsive scheduling window function Amod(t) at time t to no more than 35%, and freezing 12:00-16:00, i.e. the high evaporation period, and opening the window.

[0131] In this embodiment, step S5 of the method constructs a closed-loop matching evaluation and dynamic correction mechanism between scheduling behavior and energy consumption demand through sub-steps S51 and S52, realizing full-process linkage control from energy efficiency identification of water release behavior to adaptive optimization of strategy. Specifically, S51 introduces the "behavior energy consumption matching degree index Uctrl" to comprehensively analyze the energy adaptability between the actual energy consumption per unit area (qactual) and the evaporation heat flux (qevp), and simultaneously combines the strategy fit between the responsive scheduling window function (Amod) and the actual water release flow rate (Qactual). A composite matching factor is formed through the two dimensionless ratios, and then the integral average is processed within a set time interval to accurately quantify whether the current scheduling behavior is consistent with the energy demand and strategy objectives of the evaporation environment. Based on this, S52 constructs a targeted secondary comparison evaluation decision logic by performing hierarchical judgment on the value of the behavior energy consumption matching degree index Uctrl. When the value is less than 0.85, it indicates that the scheduling behavior is well controlled and the energy consumption is well matched, and the current strategy is maintained. When it is between 0.85 and 1, it indicates that the control is approaching the critical point, and the system automatically tightens the response sensitivity coefficient α to iteratively compress the scheduling window function Amod to improve the matching degree. When it is greater than or equal to 1, it is considered that the strategy is significantly out of control, and the system will implement a dual intervention of "forced compression + time freeze" to close the window during the high evaporation period, immediately curb excessive behavior, and prevent energy waste and ineffective loss of water resources. This implementation introduces a comprehensive evaluation index, Uctrl, that integrates behavior, energy consumption, and strategy. It establishes a quantitative basis for the behavioral energy efficiency of the scheduling system, addressing the lack of perception of actual energy consumption in traditional water resource scheduling. A differentiated adaptive adjustment mechanism of disturbance-free iteration and strategy freezing is constructed, enabling the system to respond to both minor deviations and severe loss of control, ensuring both flexibility and rigidity in regulation. The behavioral energy consumption matching index Uctrl is highly coupled with actual behavioral parameters, reflecting whether behavior closely follows the strategy capacity boundary, enhancing the system's dynamic identification and rapid correction capabilities for strategy execution deviations. By freezing open window permissions during high evaporation periods, energy waste and evaporation loss during critical periods are effectively suppressed, strengthening the control stability against extreme weather or abnormal water levels.

[0132] Although embodiments of the invention have been shown and described, it will be understood by those skilled in the art that various changes, modifications, substitutions and alterations can be made to these embodiments without departing from the principles and spirit of the invention, the scope of which is defined by the appended claims and their equivalents.

Claims

1. A reservoir evaporation control method that fuses evaporation law analysis and dispatch scenario simulation, characterized in that: The method comprises the following steps: S1, arranging remote sensing microclimate collection points in the reservoir area, and collecting evaporation environment data by using a sensor group, and transmitting the evaporation environment data to an edge computing terminal; S2, performing time series data regularization and preprocessing on the evaporation environment data in the edge computing terminal, obtaining a dimensionless data set, and establishing evaporation heat flux qevp and condensation compensation heat flux qcnd based on the dimensionless data set; S3, calculating the day-night water vapor energy difference factor ΔEdn based on the evaporation heat flux qevp and the condensation compensation heat flux qcnd, and preliminarily comparing and evaluating according to the day-night water vapor energy difference factor ΔEdn to judge the condensation compensation potential; S4, based on the preliminary comparison and evaluation result, triggering the adjustment response mechanism, outputting the responsive scheduling window function Amod based on the day-night water vapor energy difference factor ΔEdn, and adjusting the water release capacity based on the responsive scheduling window function Amod; S5, associating and calculating the responsive scheduling window function Amod and the evaporation heat flux qevp, outputting the behavior energy consumption matching degree index, and performing secondary comparison and evaluation based on the behavior energy consumption matching degree index Uctrl to execute the corresponding scheduling strategy.

2. The reservoir evaporation control method of claim 1, wherein: The S1 comprises S11 and S12; S11, arranging remote sensing microclimate collection points in the reservoir area, wherein the remote sensing microclimate collection points comprise main dam area remote sensing points, intake area remote sensing points and branch channel shallow bay area remote sensing points; and a sensor group is arranged at each remote sensing microclimate collection point, and the sampling frequency of the sensor group is set to 5 minutes each time, and the evaporation environment data of each remote sensing microclimate collection point is collected in real time; The main dam area remote sensing points are arranged in the main dam area of the reservoir; The intake area remote sensing points are arranged in the intake area of the reservoir; The branch channel shallow bay area remote sensing points are arranged in the branch channel shallow bay area of the reservoir; The sensor group comprises a total radiation table, a layered temperature and humidity probe, an aerosol laser radar, a water vapor density sensor and a temperature and humidity sensor; The evaporation environment data comprises a thermal radiation intensity parameter Rsol, a night-time convection gradient factor Xconv, an upper air aerosol condensation potential energy factor Rgsa and a water vapor disturbance frequency Vwv; S12, transmitting the evaporation environment data collected by each remote sensing microclimate collection point to the edge computing terminal by using the LoRa communication protocol through the communication module arranged in the sensor group, and synchronizing the evaporation environment data by using the GNSS time stamp.

3. The reservoir evaporation control method of claim 1, wherein: The S2 comprises S21; S21, receiving the evaporation environment data in real time at the edge computing terminal, and performing time sequence alignment on all parameters through a sliding window mechanism to obtain a structured data matrix Xsync={x ij}n×m, wherein x ij represents the jth evaporation environment data at the ith time point, n represents the number of time points, and m represents the number of types of evaporation environment parameters.

4. The reservoir evaporation control method of claim 3, wherein: The S2 further comprises S22; S22, preprocessing the structured data matrix Xsync, wherein the preprocessing is performed by respectively converting the evaporation environment data of each column in the structured data matrix Xsync into a standard dimensionless parameter with a standard deviation of 1 and a mean value of 0 by using the Z-score standardization method, to obtain a dimensionless data set.

5. The reservoir evaporation control method of claim 4, wherein: The S2 further comprises S23; S23, combining and calculating the solar radiation intensity parameter Rsol and the water vapor disturbance frequency parameter Vwv in the dimensionless data set by using an evaporation driving factor fitting model to construct the evaporation heat flux parameter qevp; The night-time convection gradient factor Xconv, the high-altitude aerosol condensation potential energy factor Rgsa, and the water vapor disturbance frequency parameter Vwv in the dimensionless data set are combined and processed through a condensation nucleus potential modeling process to construct a condensation re-supply heat flux parameter qcnd.

6. The reservoir evaporation control method of claim 4, wherein: The S3 includes S31. In the S31, the heat per unit area released by the evaporation process and the heat per unit area absorbed by the condensation process are integrated and compared in the respective time periods based on the time period integration operation of the evaporation heat flux parameter qevp and the condensation re-supply heat flux parameter qcnd, to obtain a day-night water vapor energy difference factor ΔEdn, which is calculated by subtracting the integral value of the condensation heat flux in the night-time period tn from the integral value of the evaporation heat flux in the day-time period td.

7. The reservoir evaporation control method of claim 6, wherein: The S3 also includes S32. In the S32, a condensation compensation potential reference threshold Eth is preset, which is based on the maximum recoverable energy under the optimal condensation capacity condition in the historical monitoring data as the boundary. In the preliminary comparison and evaluation of the condensation compensation potential reference threshold Eth and the day-night water vapor energy difference factor ΔEdn, the re-supply potential of the night-time condensation and the condensation compensation potential of the evaporation loss in the reservoir water system are judged, and the specific evaluation content is as follows: When the day-night water vapor energy difference factor ΔEdn is greater than or equal to the condensation compensation potential reference threshold Eth, it indicates that the condensation is insufficient and the evaporation is abnormal, and at this time, the adjustment response mechanism is triggered immediately. When the day-night water vapor energy difference factor ΔEdn is less than the condensation compensation potential reference threshold Eth, it indicates that the day-night energy balance is balanced and the condensation is sufficient, and at this time, no intervention is required.

8. The reservoir evaporation control method of claim 7, wherein: The S4 includes S41. In the S41, after the preliminary comparison and evaluation triggers the adjustment response mechanism, a type of scheduling window function mapping model is constructed according to the numerical results of the day-night water vapor energy difference factor ΔEdn and combined with the heat balance regulation principle, the type of scheduling window function mapping model takes the day-night water vapor energy difference factor ΔEdn as the input variable and the adjustment time window width as the output target, and dynamically generates a responsive scheduling window function Amod; the responsive scheduling window function Amod is an adjustment capacity function, and the responsive scheduling window function Amod is taken as a control input parameter and transmitted to the variable water release capacity regulation module for preliminary actual water level adjustment.

9. The reservoir evaporation control method of claim 8, wherein: The S5 includes S51. S51, after the preliminary actual water level adjustment, the variable water release capacity control module of the reservoir dispatching system extracts the unit area energy consumption qactual(t) corresponding to the actual dispatching behavior at time t and the actual water release flow Qactual(t) at time t. Within the set observation period T, the evaporation heat flux parameter qevp(t) at time t is calculated by ratio calculation with the unit area energy consumption qactual(t) corresponding to the actual dispatching behavior at time t, to form an instantaneous matching index. Combined with the ratio weight of the responsive dispatching window function Amod(t) at the current time point t and the actual water release flow Qactual(t) at time t, a combined kernel function is formed. Then the combined kernel function is integrated and averaged in the time interval [0, T] to output the behavior energy consumption matching degree index Uctrl.

10. The reservoir evaporation control method of claim 9, wherein: The S5 further comprises S52; S52, based on the behavior energy consumption matching degree index Uctrl, secondary comparison and evaluation are carried out to judge the control state of the preliminary actual water level adjustment, and based on the secondary comparison and evaluation result, the corresponding dispatching strategy is executed. The specific evaluation contents are as follows: When the behavior energy consumption matching degree index Uctrl is less than 0.85, it indicates that the control state is ideal and the behavior is reasonable, and the current dispatching strategy is maintained. When 0.85≤ behavior energy consumption matching degree index Uctrl<1, it indicates that the control state deviates from the critical state, and there is a slight deviation. At this time, the tightening strategy is executed, and the responsive dispatching window function Amod(t) at time t is iteratively regenerated until the control state is ideal and the iteration is stopped. When the behavior energy consumption matching degree index Uctrl is greater than or equal to 1, it indicates that the control state is seriously overconsumed, and the strategy is out of control. At this time, the responsive dispatching window function Amod(t) at time t is compressed to not more than 35%, and the window opening permission is frozen from 12:00 to 16:00.

Images