Evaluation method for intraday hydrological situation change of watershed water-wind-solar integrated operation

Abstract

The invention provides a method for evaluating hydrological situation change in a watershed water-wind-light integrated operation day, and the method comprises the following steps: determining a scheduling period after watershed water-wind-light operation in combination with watershed hydrological characteristics and a reservoir scheduling mode, and determining an evaluation reference period when watershed water-wind-light does not participate in operation for each scheduling period; the hydrological situation change degree and the hydrological form change degree of the intraday hydrological situation evaluation index in the scheduling period are obtained through evaluation, and therefore the intraday hydrological situation comprehensive change degree is obtained. According to the evaluation method for the intraday hydrological situation change of watershed water-wind-light integrated operation provided by the invention, the influence of hydropower station scheduling on the intraday hydrological situation of a river can be comprehensively and accurately evaluated, a scientific basis is provided for hydropower station ecological scheduling, and harmonious development of the hydropower station and the ecological environment can be realized.

Description

Technical Field

[0001] The present invention relates to the technical field of hydrological situation assessment, and in particular to a method for assessing changes in hydrological situation during a day of integrated water-wind-solar operation in a river basin. Background Art

[0002] The natural hydrological regime is a key factor in maintaining the health and integrity of river ecosystems, and is essential for maintaining the ecological and environmental balance and biodiversity of rivers. However, with the growing demand for energy, the construction and operation of hydropower stations have had a significant impact on the natural hydrological regime of rivers. In particular, with the demand for integrated water, wind and light dispatching and operation, the water resources in the basin have been redistributed in time and space, and the natural hydrological regime of rivers has changed significantly. This change not only leads to the degradation of river ecosystems and the fragmentation of fish spawning habitats, but also significantly restricts the health and integrity of river ecosystems. Changes in the intraday hydrological regime, if the water level drops sharply, it is easy to cause the sticky fish eggs to be exposed to the air and dry and die, and the young fish to be stranded; if the water level rises sharply, it is easy to wash away the gravel at the bottom of the spawning ground and the benthic invertebrates that serve as fish bait, and at the same time cause the fish eggs to be dragged into the deep water area with low oxygen content and die. Therefore, how to identify the impact of integrated water, wind and light dispatching on the natural hydrological regime of rivers, especially the impact of intraday hydrological regime, has become an important issue for the sustainable development of hydropower in the future.

[0003] Although there are many methods for evaluating changes in river hydrological regimes, such as the Index of Hydrological Alternation (IHA), the data time scale used by these methods is at least daily, so the impact of run-of-river hydropower stations or peak-shaving hydropower stations on natural hydrological regimes is often defined as no impact, ignoring the changes in intraday hydrological regimes. In recent years, with the development of multi-energy complementary systems such as water-wind-light, the impact of intraday hydrological regime changes on river ecosystems has intensified, especially after the intraday peak-shaving operation of hydropower stations, which causes the "steep rise and fall" of runoff, significantly changing the relatively gentle intraday hydrological process, causing aquatic organisms to be washed away or stranded, and destroying the health and integrity of aquatic habitats.

[0004] Therefore, how to effectively quantify and evaluate the impact of different integrated water-wind-solar scheduling and operation conditions on daily hydrological changes, and thus adjust the integrated water-wind-solar scheduling and operation conditions according to the degree of impact to reduce the impact on the river ecosystem, is a key issue that needs to be addressed at present. Summary of the invention

[0005] In view of the defects of the prior art, the present invention provides a method for evaluating the changes in hydrological conditions during the day of integrated water-wind-solar operation in a river basin, which can effectively solve the above problems.

[0006] The technical solution adopted by the present invention is as follows:

[0007] The present invention provides a method for evaluating the change of hydrological situation during the day of integrated operation of water, wind and solar power in a river basin, comprising the following steps:

[0008] Step S1: Determine the intraday hydrological situation assessment index vector ; They represent the intraday hydrological situation assessment index 1, intraday hydrological situation assessment index 2, ..., intraday hydrological situation assessment index ;

[0009] Step S2: Determine the The scheduling period after the water, wind and solar power in the basin participate in the operation, for each scheduling period , , determine the evaluation benchmark period for the corresponding basin water, wind and solar power without participation in operation , where the scheduling period and evaluation base period Both have day;

[0010] Step S3: Collect and obtain the daily hydrological situation assessment indicators During the scheduling period The daily intraday hydrological situation assessment index value is expressed as: ,in, , , specifically: daily hydrological situation assessment index During the scheduling period No. The daily hydrological situation assessment index value of the day;

[0011] Obtain daily hydrological situation assessment indicators During the evaluation period The daily intraday hydrological situation assessment index value is expressed as: ;

[0012] Step S4: Use formula (1) to evaluate the daily hydrological situation assessment index During the scheduling period The degree of change in hydrological regime :

[0013] (1)

[0014] in: Represents the daily hydrological situation assessment index During the scheduling period Daily intraday hydrological situation assessment index values , the number of items within the RVA threshold range;

[0015] Represents the daily hydrological situation assessment index During the scheduling period Daily intraday hydrological situation assessment index values , the expected number within the RVA threshold range; Determine by formula (2):

[0016] (2)

[0017] in: Indicators for daily hydrological situation assessment During the scheduling period Daily intraday hydrological situation assessment index values , the expected proportion within the RVA threshold range;

[0018] Step S5: Evaluate and obtain the daily hydrological situation assessment index During the scheduling period The degree of hydrological change :

[0019] Step S5.1: According to the grade conversion rules, the daily hydrological situation assessment index During the scheduling period Daily intraday hydrological situation assessment index values Converted to daily hydrological situation assessment index level ; Daily hydrological situation assessment index During the evaluation period Daily intraday hydrological situation assessment index values Converted to daily hydrological situation assessment index level ;

[0020] Step S5.2, due to , so we get Hydrological situation assessment index level within a day , the daily hydrological situation assessment index level of each day Perform pairwise comparisons and establish scheduling periods based on matrix creation rules Corresponding daily hydrological situation assessment indicators The Hasse matrix is ​​expressed as: ,for matrix;

[0021] Similarly, for the evaluation base period Daily intraday hydrological situation assessment index level Make pairwise comparisons and establish the evaluation benchmark period according to the matrix creation rules. Corresponding daily hydrological situation assessment indicators The Hasse matrix is ​​expressed as: ,for matrix;

[0022] matrix and matrix The main diagonal The value of each element is 0;

[0023] Step S5.3, matrix Make corrections and integrate into the scheduling period The measured values ​​of the daily reservoir operation condition characteristic data are used to obtain the corrected matrix, which is expressed as: ;

[0024] Pair Matrix Make corrections and incorporate them into the evaluation benchmark period The measured values ​​of the daily reservoir operation condition characteristic data are used to obtain the corrected matrix, which is expressed as: ;

[0025] Step S5.4, calculate the corrected matrix And the corrected matrix The distance between the two is the daily hydrological situation assessment index. During the scheduling period The degree of hydrological change ;

[0026] Step S6: Use formula (3) to evaluate the daily hydrological situation assessment index During the scheduling period The comprehensive change of hydrological conditions within a day :

[0027] (3)

[0028] Daily hydrological situation assessment indicators During the scheduling period The comprehensive change of hydrological conditions within a day , reflecting the impact of the basin’s water, wind and solar power systems on the hydrological situation.

[0029] Preferably, the daily hydrological situation assessment index refers to an index that quantifies the degree of change in the daily hydrological situation based on the hourly data of the inflow and outflow runoff of the hydropower station.

[0030] Preferably, the daily hydrological situation assessment indicators include hourly change difference of the hydrological situation, intraday change difference, daily average runoff, flow increase time, flow decrease time, flow increase duration and flow decrease duration.

[0031] Preferably, in step S2, for each scheduling period , , determine the evaluation benchmark period for the corresponding basin water, wind and solar power without participation in operation , specifically:

[0032] Scheduling period The average water inflow over the period of time is counted to determine the scheduling period. characteristics, including flood season, normal water season and dry season;

[0033] Acquisition and dispatch period The average value of multiple historical benchmark periods with the same characteristics and without the participation of water, wind and light in the basin is used as the evaluation benchmark period. .

[0034] Preferably, in step S4, the RVA threshold range is determined by:

[0035] Obtain daily hydrological situation assessment indicators During the scheduling period The 25%-75% quantiles of multiple historical data of the same period are used as the RVA threshold range .

[0036] Preferably, in step S5.1, according to the grade conversion rule, the intraday hydrological situation assessment index During the scheduling period Daily intraday hydrological situation assessment index values Converted to daily hydrological situation assessment index level , specifically:

[0037] If the daily hydrological situation assessment index value In of range, the level of the daily hydrological situation assessment index is 1;

[0038] If the daily hydrological situation assessment index value In of range, the level of the daily hydrological situation assessment index is 2;

[0039] If the daily hydrological situation assessment index value In of range, the level of the daily hydrological situation assessment index is 3.

[0040] Preferably, in step S5.2, according to the matrix creation rule, the scheduling period is established Corresponding daily hydrological situation assessment indicators The Hasse matrix is ​​expressed as: , specifically:

[0041] According to formula (4), determine the matrix The value of each element of :

[0042] (4)

[0043] in: Represents the daily hydrological situation assessment index During the scheduling period No. The level of the daily hydrological situation assessment index of the day, ;

[0044] Representative Matrix No. Line The element value of the column;

[0045] It can be seen that the matrix The main diagonal The element value of each element is 0.

[0046] Preferably, in step S5.3, the matrix Make corrections and integrate into the scheduling period The measured values ​​of the daily reservoir operation condition characteristic data are used to obtain the corrected matrix, which is expressed as: , specifically:

[0047] Pair Matrix The main diagonal The elements of the main diagonal are modified. The element is The correction method is:

[0048] make = ,in, Representative scheduling period No. The reservoir operation condition characteristic data of the day, specifically three-dimensional characteristic data, including the operation period No. The intraday difference index of the reservoir daily discharge flow , reservoir daily instantaneous indicators The maximum daily water level variation of the reservoir .

[0049] Preferably, step S5.4 is specifically:

[0050] Step S5.4.1, the modified matrix for The matrix, whose elements are reformulated as: , represents the modified matrix Middle Line The elements of the column, where ; ;

[0051] The corrected matrix for The matrix, whose elements are reformulated as: , represents the modified matrix Middle Line The elements of the column, where ; ;

[0052] Step S5.4.2, use formula (5) to calculate the corrected matrix And the corrected matrix The distance between off-diagonal elements is:

[0053] (5)

[0054] Step S5.4.3, use formula (6) to calculate the corrected matrix And the corrected matrix The distance between diagonal elements is:

[0055] (6)

[0056] Step S5.4.4, use formula (7) to calculate the corrected matrix And the corrected matrix The comprehensive distance between :

[0057] (7)

[0058] in: It is an intraday hydrological situation assessment index between [0,1]. During the scheduling period The weighting parameter of .

[0059] Preferably, the weighted parameter The method to determine is:

[0060] Step S5.4.4.1, use formula (8) and formula (9) to calculate the scheduling period Indicators for daily hydrological situation assessment The proportion of diagonal elements And comprehensive weight :

[0061] (8)

[0062] (9)

[0063] Step S5.4.4.2, using formula (10), obtain the scheduling period Indicators for daily hydrological situation assessment Information entropy :

[0064] (10)

[0065] Step S5.4.4.3, using formula (11), obtain the weighted parameter :

[0066] (11)

[0067] This step is finished.

[0068] The method for evaluating the changes in hydrological conditions during the day of integrated water, wind and solar operation in a river basin provided by the present invention has the following advantages:

[0069] The present invention provides an evaluation method for changes in hydrological conditions during the day of integrated water-wind-solar operation in a river basin, which can effectively quantify and evaluate the impact of different integrated water-wind-solar scheduling and operation conditions on changes in hydrological conditions during the day, thereby providing a scientific basis for hydropower stations to participate in ecological scheduling under integrated water-wind-solar operation, and achieving coordinated development of hydropower operation and ecological protection. BRIEF DESCRIPTION OF THE DRAWINGS

[0070] Figure 1 A flow chart of a method for evaluating changes in hydrological conditions during a day of integrated water-wind-solar operation in a river basin provided by the present invention. DETAILED DESCRIPTION

[0071] In order to make the technical problems, technical solutions and beneficial effects solved by the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0072] The present invention provides a method for evaluating changes in the hydrological situation during the day of integrated water-wind-solar operation in a river basin, which can comprehensively and accurately evaluate the impact of hydropower station scheduling on the daily hydrological situation of the river, provide a scientific basis for the ecological scheduling of hydropower stations, and help achieve the harmonious development of hydropower stations and the ecological environment.

[0073] See also Figure 1 The present invention provides a method for evaluating the changes in hydrological conditions during the day of integrated water, wind and solar operation in a river basin, comprising the following steps:

[0074] Step S1: Determine the intraday hydrological situation assessment index vector ; They represent the intraday hydrological situation assessment index 1, intraday hydrological situation assessment index 2, ..., intraday hydrological situation assessment index ;

[0075] In the present invention, the intraday hydrological situation evaluation index refers to an index that quantifies the degree of change of the intraday hydrological situation based on the hourly data of the inflow and outflow runoff of the hydropower station. For example, the intraday hydrological situation evaluation index includes the hourly change difference of the hydrological situation, the intraday change difference, the daily average runoff, the time of flow increase, the time of flow decrease, the duration of flow increase and the duration of flow decrease, etc. The present invention does not limit the specific intraday hydrological situation evaluation index adopted.

[0076] Specifically, the intraday hydrological situation assessment indicators may include indicators on the magnitude, change rate, duration, frequency and timing of the intraday hydrological situation, among which magnitude indicators: hourly change difference, intraday change difference and daily average runoff; change rate indicators: average flow increase rate, average flow decrease rate, maximum flow increase rate, maximum flow decrease rate and instantaneous index; duration indicators: the moment of flow increase and the moment of flow decrease; timing indicators: the duration of flow increase and the duration of flow decrease; frequency indicators: flow reversal frequency, flow increase frequency and flow decrease frequency.

[0077] Further, combined with Table 1, the indicators and their meanings are as follows:

[0078] (1) Magnitude index. This index mainly evaluates the difference between high and low flow in a day, including the daily average runoff, the difference between the maximum and minimum flow in a day, and the maximum difference in runoff between consecutive hours.

[0079] (2) Change rate index. This index mainly evaluates the dynamic change characteristics of the daily increase and decrease of runoff, including the flow increase rate, flow decrease rate and instantaneous index.

[0080] (3) Duration index: This index is used to evaluate the duration of the process of runoff increase or decrease caused by the water storage and release scheduling of the hydropower station.

[0081] (4) Frequency index. This index evaluates the frequency of runoff increase or decrease within a day. It is calculated by dividing the time of runoff increase or decrease by the total time of the day (24 hours). In addition, the flow reversal frequency is defined as the number of times the runoff switches between high and low flows.

[0082] (5) Timing index. Since the water storage and release scheduling of the hydropower station changes the daily distribution process of runoff, the timing of runoff increase and decrease also changes. Therefore, this index is used to evaluate the time when the runoff increases to the maximum or decreases to the minimum.

[0083] Table 1 Indicators and meanings of intraday hydrological situation assessment

[0084]

[0085] Step S2: Determine the The scheduling period after the water, wind and solar power in the basin participate in the operation, for each scheduling period , , determine the evaluation benchmark period for the corresponding basin water, wind and solar power without participation in operation , where the scheduling period and evaluation base period Both have day;

[0086] In this step, different dispatching periods are divided according to the changing trend of the basin hydrological situation and the law of reservoir dispatching.

[0087] In this step, for each scheduling period , , determine the evaluation benchmark period for the corresponding basin water, wind and solar power without participation in operation , specifically:

[0088] Scheduling period The average water inflow over the period of time is counted to determine the scheduling period. characteristics, including flood season, normal water season and dry season;

[0089] Acquisition and dispatch period The average value of multiple historical benchmark periods with the same characteristics and without the participation of water, wind and light in the basin is used as the evaluation benchmark period. .

[0090] Step S3: Collect and obtain the daily hydrological situation assessment indicators During the scheduling period The daily intraday hydrological situation assessment index value is expressed as: ,in, , , specifically: daily hydrological situation assessment index During the scheduling period No. The daily hydrological situation assessment index value of the day;

[0091] Obtain daily hydrological situation assessment indicators During the evaluation period The daily intraday hydrological situation assessment index value is expressed as: ;

[0092] Step S4: Use formula (1) to evaluate the daily hydrological situation assessment index During the scheduling period The degree of change in hydrological regime :

[0093] (1)

[0094] in: Represents the daily hydrological situation assessment index During the scheduling period Daily intraday hydrological situation assessment index values , the number of items within the RVA threshold range;

[0095] Represents the daily hydrological situation assessment index During the scheduling period Daily intraday hydrological situation assessment index values , the expected number within the RVA threshold range; Determine by formula (2):

[0096] (2)

[0097] in: Indicators for daily hydrological situation assessment During the scheduling period Daily intraday hydrological situation assessment index values , the expected proportion within the RVA threshold range;

[0098] In this step, the RVA threshold range determination method can be: Obtain the intraday hydrological situation assessment index During the scheduling period The 25%-75% quantiles of multiple historical data of the same period are used as the RVA threshold range .

[0099] Step S5: Evaluate and obtain the daily hydrological situation assessment index During the scheduling period The degree of hydrological change .

[0100] Step S6: Use formula (3) to evaluate the daily hydrological situation assessment index During the scheduling period The comprehensive change of hydrological conditions within a day :

[0101] (3)

[0102] Daily hydrological situation assessment indicators During the scheduling period The comprehensive change of hydrological conditions within a day , reflecting the impact of the basin’s water, wind and solar power systems on the hydrological situation.

[0103] In the present invention, step S5, the daily hydrological situation evaluation index is evaluated. During the scheduling period The degree of hydrological change , which is an important innovative step of the present invention, and its specific implementation steps are:

[0104] Step S5.1: According to the grade conversion rules, the daily hydrological situation assessment index During the scheduling period Daily intraday hydrological situation assessment index values Converted to daily hydrological situation assessment index level ; Daily hydrological situation assessment index During the evaluation period Daily intraday hydrological situation assessment index values Converted to daily hydrological situation assessment index level ;

[0105] Daily hydrological situation assessment indicators During the scheduling period Daily intraday hydrological situation assessment index values Converted to daily hydrological situation assessment index level Methods and indicators for daily hydrological situation assessment During the evaluation period Daily intraday hydrological situation assessment index values Converted to daily hydrological situation assessment index level The method and principle are the same, so only the intraday hydrological situation assessment index is used. During the scheduling period Daily intraday hydrological situation assessment index values Converted to daily hydrological situation assessment index level For example, convert in the following way:

[0106] If the daily hydrological situation assessment index value In of range, the level of the daily hydrological situation assessment index is 1;

[0107] If the daily hydrological situation assessment index value In of range, the level of the daily hydrological situation assessment index is 2;

[0108] If the daily hydrological situation assessment index value In of range, the level of the daily hydrological situation assessment index is 3.

[0109] Therefore, through the above method, the daily hydrological situation assessment index value Convert to level 1, 2, or 3, corresponding to low, medium, and high levels of change.

[0110] Step S5.2, due to , so we get Hydrological situation assessment index level within a day , the daily hydrological situation assessment index level of each day Perform pairwise comparisons and establish scheduling periods based on matrix creation rules Corresponding daily hydrological situation assessment indicators The Hasse matrix is ​​expressed as: ,for matrix;

[0111] Similarly, for the evaluation base period Daily intraday hydrological situation assessment index level Make pairwise comparisons and establish the evaluation benchmark period according to the matrix creation rules. Corresponding daily hydrological situation assessment indicators The Hasse matrix is ​​expressed as: ,for matrix; matrix and matrix The main diagonal The value of each element is 0.

[0112] In this step, the matrix and matrix The creation principle is the same as that of For example, the creation method is:

[0113] According to formula (4), determine the matrix The value of each element of :

[0114] (4)

[0115] in: Represents the daily hydrological situation assessment index During the scheduling period No. The level of the daily hydrological situation assessment index of the day, ; Representative Matrix No. Line The element value of the column;

[0116] It can be seen that the matrix The main diagonal The element values ​​of all elements are 0, and the element values ​​of non-diagonal elements are +1 or -1.

[0117] Step S5.3, matrix Make corrections and integrate into the scheduling period The measured values ​​of the daily reservoir operation condition characteristic data are used to obtain the corrected matrix, which is expressed as: ;

[0118] Pair Matrix Make corrections and incorporate them into the evaluation benchmark period The measured values ​​of the daily reservoir operation condition characteristic data are used to obtain the corrected matrix, which is expressed as: ;

[0119] In this step, the matrix Make corrections and adjust the matrix The principle of correction is the same, only the matrix Take the correction method as an example, and use the following method to correct it:

[0120] Pair Matrix The main diagonal The elements of the main diagonal are modified. The element is The correction method is:

[0121] make = ,in, Representative scheduling period No. The reservoir operation condition characteristic data of the day, specifically three-dimensional characteristic data, including the operation period No. The intraday difference index of the reservoir daily discharge flow , reservoir daily instantaneous indicators The maximum daily water level variation of the reservoir .

[0122] Therefore, in the present invention, the uncorrected matrix and matrix , are antisymmetric matrices, consisting only of 0 and ±1, and the main diagonal elements are all 0. By adding the reservoir's dispatching condition characteristic data to the main diagonal of the matrix, the connection between the daily hydrological situation evaluation indicators and the reservoir dispatching conditions is enhanced.

[0123] Scheduling period No. The intraday difference index of the reservoir daily discharge flow , reservoir daily instantaneous indicators The maximum daily water level variation of the reservoir , the specific meaning is:

[0124] Intraday Difference Index :For the scheduling period No. The flow difference between the maximum and minimum flow of the day;

[0125] Maximum daily water level variation of reservoir :For the scheduling period No. The difference between the maximum and minimum reservoir water levels on a given day;

[0126] Reservoir intraday instantaneous indicators : is the instantaneous index after Z-score normalization; specifically, in order to smooth the short-term abnormal fluctuation of traffic and capture the trend of traffic change, the weighted moving average method (EWMA) is used to calculate the instantaneous index, and the instantaneous index is normalized by Z-score normalization. The calculation formula is as follows:

[0127]

[0128]

[0129] in: is the flow difference between consecutive periods; is the instantaneous index at time t; is the instantaneous index at time t-1; is the smoothing coefficient; is the instantaneous index after normalization of Z-score; μI is the mean of the instantaneous index at all times of the day; σ I is the standard deviation of the instantaneous index.

[0130] Furthermore, in order to calculate the Hasse matrix, the reservoir operation condition characteristic data is normalized to be in the interval [0,1]. The normalization method is: determine the main diagonal The maximum value of the characteristic data of the reservoir operation condition is The daily reservoir operation condition characteristic data are divided by the maximum value for normalization.

[0131] Step S5.4, calculate the corrected matrix And the corrected matrix The distance between the two is the daily hydrological situation assessment index. During the scheduling period The degree of hydrological change ;

[0132] Step S5.4 is specifically as follows:

[0133] Step S5.4.1, the modified matrix for The matrix, whose elements are reformulated as: , represents the modified matrix Middle Line The elements of the column, where ; ;

[0134] The corrected matrix for The matrix, whose elements are reformulated as: , represents the modified matrix Middle Line The elements of the column, where ; ;

[0135] Step S5.4.2, use formula (5) to calculate the corrected matrix And the corrected matrix The distance between off-diagonal elements is:

[0136] (5)

[0137] Step S5.4.3, use formula (6) to calculate the corrected matrix And the corrected matrix The distance between diagonal elements is:

[0138] (6)

[0139] Step S5.4.4, use formula (7) to calculate the corrected matrix And the corrected matrix The comprehensive distance between :

[0140] (7)

[0141] in: It is an intraday hydrological situation assessment index between [0,1]. During the scheduling period The weighting parameter of .

[0142] In order to accurately quantify the changes in the hydrological situation during reservoir operation and comprehensively reflect the changes in the hydrological situation evaluation indicators on different days, the entropy weight method is used to dynamically calculate the weights of the diagonal elements and off-diagonal elements in each operation period k, and comprehensively evaluate the Hasse distance of the diagonal elements and off-diagonal elements. Specifically, the weighted parameters The method to determine is:

[0143] Step S5.4.4.1, entropy weight calculation: Since the diagonal elements and off-diagonal elements of the Hasse matrix are in the interval [0,1], no standardization is required. Formula (8) and formula (9) are directly used to calculate the scheduling period respectively. Indicators for daily hydrological situation assessment The proportion of diagonal elements And comprehensive weight :

[0144] (8)

[0145] (9)

[0146] Step S5.4.4.2, calculate information entropy: Use formula (10) to get the scheduling period Indicators for daily hydrological situation assessment Information entropy :

[0147] (10)

[0148] Step S5.4.4.3, using formula (11), obtain the weighted parameter :

[0149] (11)

[0150] This step is finished.

[0151] The present invention provides a method for evaluating the changes in hydrological conditions during the day of integrated water, wind and solar operation in a river basin. The idea can be summarized as follows:

[0152] Step 1: Collect hourly data on the inflow and outflow of hydropower stations, water levels, etc. in each dispatching period before and after the water-wind-solar integration operation, and count the daily changes in reservoir operation conditions and basin hydrological conditions;

[0153] Step 2: Define the daily hydrological situation assessment indicators based on the watershed eco-hydrological needs;

[0154] Step 3: Based on the basin hydrological characteristics and reservoir operation mode, divide the reservoir operation scheduling period and determine the evaluation benchmark period corresponding to the scheduling period;

[0155] Based on the defined daily hydrological situation assessment indicators, the RVA method is used to evaluate the hydrological situation change degree during the dispatch period.

[0156] Step 4: Based on the range of changes in the daily hydrological situation assessment index during the dispatching period determined by the RVA method, the hydrological situation change level is divided, so as to convert the daily hydrological situation assessment index values ​​of the dispatching period and the evaluation base period into the daily hydrological situation assessment index level;

[0157] Step 5: Compare the levels of the daily intraday hydrological situation assessment indicators during the dispatch period in pairs, and construct the Hasse matrix corresponding to the dispatch period; similarly, compare the levels of the daily intraday hydrological situation assessment indicators during the evaluation base period in pairs, and construct the Hasse matrix corresponding to the evaluation base period;

[0158] The Hasse matrix corresponding to the dispatching period and the Hasse matrix corresponding to the evaluation benchmark period are respectively modified, and the diagonal elements of the matrix are added to the reservoir dispatching condition characteristic data of the corresponding day, so as to obtain the modified Hasse matrix corresponding to the dispatching period and the modified Hasse matrix corresponding to the evaluation benchmark period;

[0159] Step 6: Combine the entropy weight method to calculate the distance between the Hasse matrix corresponding to the revised scheduling period and the Hasse matrix corresponding to the revised evaluation benchmark period, and obtain the hydrological morphological change degree of each intraday hydrological situation assessment indicator;

[0160] Step 7: Integrate the hydrological regime change degree and hydrological morphological change degree of the intraday hydrological regime assessment indicators to obtain the comprehensive intraday hydrological regime change degree of each intraday hydrological regime assessment indicator in each scheduling period.

[0161] The present invention provides an evaluation method for changes in hydrological conditions during the day of integrated water-wind-solar operation in a river basin, which can effectively quantify and evaluate the impact of different integrated water-wind-solar scheduling and operation conditions on changes in hydrological conditions during the day, thereby providing a scientific basis for hydropower stations to participate in ecological scheduling under integrated water-wind-solar operation, and achieving coordinated development of hydropower operation and ecological protection.

[0162] The above is only a preferred embodiment of the present invention. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principle of the present invention. These improvements and modifications should also be considered as the scope of protection of the present invention.

Claims

1. A method for evaluating the changes in hydrological conditions during the day of integrated water, wind and solar power operation in a river basin, characterized in that: The following steps are involved: Step S1: Determine the intraday hydrological situation assessment index vector ; They represent the intraday hydrological situation assessment index 1, intraday hydrological situation assessment index 2, ..., intraday hydrological situation assessment index ; Step S2: Determine the The scheduling period after the water, wind and solar power in the basin participate in the operation, for each scheduling period , , determine the evaluation benchmark period for the corresponding basin water, wind and solar power without participation in operation , where the scheduling period and evaluation base period Both have day; Step S3: Collect and obtain the daily hydrological situation assessment indicators During the scheduling period The daily intraday hydrological situation assessment index value is expressed as: ,in, , , specifically: daily hydrological situation assessment index During the scheduling period No. The daily hydrological situation assessment index value of the day; Obtain daily hydrological situation assessment indicators During the evaluation period The daily intraday hydrological situation assessment index value is expressed as: ; Step S4: Use formula (1) to evaluate the daily hydrological situation assessment index During the scheduling period The degree of change in hydrological regime : （1） in: Represents the daily hydrological situation assessment index During the scheduling period Daily intraday hydrological situation assessment index values , the number of items within the RVA threshold range; Represents the daily hydrological situation assessment index During the scheduling period Daily intraday hydrological situation assessment index values , the expected number within the RVA threshold range; Determine by formula (2): （2） in: Indicators for daily hydrological situation assessment During the scheduling period Daily intraday hydrological situation assessment index values , the expected proportion within the RVA threshold range; Step S5: Evaluate and obtain the daily hydrological situation assessment index During the scheduling period The degree of hydrological change : Step S5.1: According to the grade conversion rules, the daily hydrological situation assessment index During the scheduling period Daily intraday hydrological situation assessment index values Converted to daily hydrological situation assessment index level ; Daily hydrological situation assessment index During the evaluation period Daily intraday hydrological situation assessment index values Converted to daily hydrological situation assessment index level ; Step S5.2, due to , so we get Hydrological situation assessment index level within a day , the daily hydrological situation assessment index level of each day Perform pairwise comparisons and establish scheduling periods based on matrix creation rules Corresponding daily hydrological situation assessment indicators The Hasse matrix is ​​expressed as: ,for matrix; Similarly, for the evaluation base period Daily intraday hydrological situation assessment index level Make pairwise comparisons and establish the evaluation benchmark period according to the matrix creation rules. Corresponding daily hydrological situation assessment indicators The Hasse matrix is ​​expressed as: ,for matrix; matrix and matrix The main diagonal The value of each element is 0; Step S5.3, matrix Make corrections and integrate into the scheduling period The measured values ​​of the daily reservoir operation condition characteristic data are used to obtain the corrected matrix, which is expressed as: ; Pair Matrix Make corrections and incorporate them into the evaluation benchmark period The measured values ​​of the daily reservoir operation condition characteristic data are used to obtain the corrected matrix, which is expressed as: ; Step S5.4, calculate the corrected matrix And the corrected matrix The distance between the two is the daily hydrological situation assessment index. During the scheduling period The degree of hydrological change ; Step S6: Use formula (3) to evaluate the daily hydrological situation assessment index During the scheduling period The comprehensive change of hydrological conditions within a day : （3） Daily hydrological situation assessment indicators During the scheduling period The comprehensive change of hydrological conditions within a day , reflecting the impact of the basin’s water, wind and solar power systems on the hydrological situation.

2. The method for evaluating the changes in hydrological conditions during the day of integrated water, wind and solar operation in a river basin according to claim 1 is characterized in that: The daily hydrological situation assessment index refers to an index that quantifies the degree of change in the daily hydrological situation based on the hourly data of inflow and outflow runoff of hydropower stations.

3. The method for evaluating the changes in hydrological conditions during the day of integrated water, wind and solar operation in a river basin according to claim 1 is characterized in that: The daily hydrological situation assessment indicators include the hourly change difference of the hydrological situation, the intraday change difference, the daily average runoff, the time of flow increase, the time of flow decrease, the duration of flow increase and the duration of flow decrease.

4. The method for evaluating the changes in hydrological conditions during the day of integrated water, wind and solar operation in a river basin according to claim 1 is characterized in that: In step S2, for each scheduling period , , determine the evaluation benchmark period for the corresponding basin water, wind and solar power without participation in operation , specifically: Scheduling period The average water inflow over the period of time is counted to determine the scheduling period. characteristics, including flood season, normal water season and dry season; Acquisition and dispatch period The average value of multiple historical benchmark periods with the same characteristics and without the participation of water, wind and light in the basin is used as the evaluation benchmark period. .

5. The method for evaluating the changes in hydrological conditions during the day of integrated water, wind and solar operation in a river basin according to claim 1 is characterized in that: In step S4, the RVA threshold range is determined by: Obtain daily hydrological situation assessment indicators During the scheduling period The 25%-75% quantiles of multiple historical data of the same period are used as the RVA threshold range .

6. The method for evaluating the changes in hydrological conditions during the day of integrated water, wind and solar operation in a river basin according to claim 5 is characterized in that: In step S5.1, according to the grade conversion rules, the daily hydrological situation assessment index During the scheduling period Daily intraday hydrological situation assessment index values Converted to daily hydrological situation assessment index level , specifically: If the daily hydrological situation assessment index value In of range, the level of the daily hydrological situation assessment index is 1; If the daily hydrological situation assessment index value In of range, the level of the daily hydrological situation assessment index is 2; If the daily hydrological situation assessment index value In of range, the level of the daily hydrological situation assessment index is 3.

7. The method for evaluating the changes in hydrological conditions during the day of integrated water, wind and solar operation in a river basin according to claim 1 is characterized in that: In step S5.2, according to the matrix creation rule, the scheduling period is established Corresponding daily hydrological situation assessment indicators The Hasse matrix is ​​expressed as: , specifically: According to formula (4), determine the matrix The value of each element of : （4） in: Represents the daily hydrological situation assessment index During the scheduling period No. The level of the daily hydrological situation assessment index of the day, ; Representative Matrix No. Line The element value of the column; It can be seen that the matrix The main diagonal The element value of each element is 0.

8. The method for evaluating the changes in hydrological conditions during the day of integrated water, wind and solar operation in a river basin according to claim 1 is characterized in that: In step S5.3, the matrix Make corrections and integrate into the scheduling period The measured values ​​of the daily reservoir operation condition characteristic data are used to obtain the corrected matrix, which is expressed as: , specifically: Pair Matrix The main diagonal The elements of the main diagonal are modified. The element is The correction method is: make = ,in, Representative scheduling period No. The reservoir operation condition characteristic data of the day, specifically three-dimensional characteristic data, including the operation period No. The intraday difference index of the reservoir daily discharge flow , reservoir daily instantaneous indicators The maximum daily water level variation of the reservoir .

9. The method for evaluating the changes in hydrological conditions during the day of integrated water, wind and solar operation in a river basin according to claim 1 is characterized in that: Step S5.4 is specifically as follows: Step S5.4.1, the modified matrix for The matrix, whose elements are reformulated as: , represents the modified matrix Middle Line The elements of the column, where ; ; The corrected matrix for The matrix, whose elements are reformulated as: , represents the modified matrix Middle Line The elements of the column, where ; ; Step S5.4.2, use formula (5) to calculate the corrected matrix And the corrected matrix The distance between off-diagonal elements is: （5） Step S5.4.3, use formula (6) to calculate the corrected matrix And the corrected matrix The distance between diagonal elements is: （6） Step S5.4.4, use formula (7) to calculate the corrected matrix And the corrected matrix The comprehensive distance between : （7） in: It is an intraday hydrological situation assessment index between [0,1]. During the scheduling period The weighting parameter of .

10. The method for evaluating the changes in hydrological conditions during the day of integrated water, wind and solar operation in a river basin according to claim 9 is characterized in that: Weighting parameters The method to determine is: Step S5.4.4.1, use formula (8) and formula (9) to calculate the scheduling period Indicators for daily hydrological situation assessment The proportion of diagonal elements And comprehensive weight : （8） （9） Step S5.4.4.2, using formula (10), obtain the scheduling period Indicators for daily hydrological situation assessment Information entropy : （10） Step S5.4.4.3, using formula (11), obtain the weighted parameter : （11） This step is finished.

Images