Method for evaluating intraday hydrological regime alterations under integrated operation of hydropower, wind power and photovoltaic power in river basin

Abstract

Provided is a method for evaluating intraday hydrological regime alterations under integrated operation of hydropower, wind power and photovoltaic power in a river basin, the method comprising the following steps: on the basis of river basin hydrological characteristics and a reservoir operation mode, determining scheduling periods after integration of hydropower, wind power and photovoltaic power into river basin operation, and, for each scheduling period, determining a corresponding evaluation baseline period without integration of hydropower, wind power and photovoltaic power into river basin operation; and evaluating hydrological regime alteration degrees and hydrological pattern alteration degrees of intraday hydrological regime evaluation indicators during each scheduling period, so as to obtain comprehensive intraday hydrological regime alteration degrees, thereby achieving comprehensive and accurate evaluation of the impact of hydropower station scheduling on the intraday hydrological regime of rivers.

Description

A method for assessing intraday hydrological changes in integrated watershed operation involving water, wind, and solar power. Technical Field

[0001] This invention relates to the field of hydrological situation assessment technology, specifically to a method for assessing intraday hydrological situation changes during the integrated operation of watershed water, wind and solar power. Background Technology

[0002] Natural hydrological conditions are a key factor in maintaining the health and integrity of river ecosystems, and are crucial for preserving the ecological balance and biodiversity of rivers. However, with the increasing demand for energy from humans, the construction and operation of hydropower stations have had a significant impact on the natural hydrological conditions of rivers. In particular, with the emergence of integrated hydropower-wind-solar power scheduling, water resources in the basin are being redistributed in time and space, resulting in significant changes in the natural hydrological conditions of rivers. These changes not only lead to the degradation of river ecosystems and the fragmentation of fish spawning habitats, but also significantly restrict the health and integrity of river ecosystems. Changes in intraday hydrological conditions, such as a sharp drop in water level, can easily cause adhesive fish eggs to dry out and die when exposed to air, and juvenile fish to become stranded; conversely, a sharp rise in water level can easily wash away the gravel and pebbles at the bottom of the spawning grounds, as well as benthic invertebrates that serve as food for fish, while also causing fish eggs to be dragged into deeper waters with lower oxygen levels and die. Therefore, how to identify the impact of integrated hydropower-wind-solar power scheduling on the natural hydrological conditions of rivers, especially the impact on intraday hydrological conditions, has become an important issue for the sustainable development of hydropower in the future.

[0003] Although various methods exist for assessing changes in river hydrological conditions, such as the Index of Hydrologic Alternation (IHA), these methods use a minimum daily timescale for data. Therefore, they often define the impact of run-of-river hydropower stations or peak-shaving hydropower stations on natural hydrological conditions as having no impact, neglecting intraday changes in hydrological conditions. In recent years, with the development of multi-energy complementary systems such as water-wind-solar power, the impact of intraday hydrological changes on river ecosystems has intensified. In particular, after hydropower stations operate during intraday peak-shaving, they cause "sharp rises and falls" in runoff, significantly altering the relatively gentle intraday hydrological processes, leading to the erosion or stranding of aquatic organisms and damaging the health and integrity of their habitats.

[0004] Therefore, how to effectively quantify and assess the impact of different integrated water, wind, and solar power scheduling operation conditions on intraday hydrological changes, and thus adjust the integrated water, wind, and solar power scheduling 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. Summary of the Invention

[0005] To address the shortcomings of existing technologies, this invention provides a method for assessing intraday hydrological changes in integrated watershed operation involving water, wind, and solar power, which can effectively solve the aforementioned problems.

[0006] The technical solution adopted in this invention is as follows:

[0007] This invention provides a method for assessing intraday hydrological changes in a watershed's integrated water, wind, and solar power operation, comprising the following steps:

[0008] Step S1: Determine the intraday hydrological situation assessment index vector X = {x1, x2, ..., x...} n};x1,x2,...,x n These represent intraday hydrological situation assessment index 1, intraday hydrological situation assessment index 2, ..., intraday hydrological situation assessment index n, respectively.

[0009] Step S2: Combining the watershed hydrological characteristics and reservoir scheduling methods, determine the scheduling period after N types of watershed water, wind and solar power participate in operation. For each scheduling period k, k = 1, 2, ..., N, determine the corresponding evaluation benchmark period s for watershed water, wind and solar power not participating in operation. Both the scheduling period k and the evaluation benchmark period s have m days.

[0010] Step S3: Collect the intraday hydrological situation assessment index x i The daily intraday hydrological situation assessment index values ​​for the scheduling period k are expressed as follows: Where i = 1, 2, ..., n, j = 1, 2, ..., m, the specific meaning is: intraday hydrological situation assessment index x i The intraday hydrological situation assessment index value on the j-th day of the scheduling period k;

[0011] Obtain intraday hydrological situation assessment indicators x i The daily intraday hydrological situation assessment index values ​​for the baseline period s are expressed as:

[0012] Step S4: Using formula (1), the intraday hydrological situation assessment index x is obtained. i Changes in hydrological conditions during the scheduling period k

[0013] in: The index representing the intraday hydrological situation assessment is x i Daily intraday hydrological situation assessment index values ​​during scheduling period k The number of items within the RVA threshold range;

[0014] The index representing the intraday hydrological situation assessment is x iDaily intraday hydrological situation assessment index values ​​during scheduling period k The expected number within the RVA threshold range; Determined by formula (2):

[0015] Where: p k For intraday hydrological situation assessment indicators x i Daily intraday hydrological situation assessment index values ​​during scheduling period k The expected proportion within the RVA threshold range;

[0016] Step S5: Evaluate and obtain the intraday hydrological situation assessment index x i Changes in hydrological morphology during the scheduling period k

[0017] Step S5.1, according to the level conversion rules, adjust the intraday hydrological situation assessment index x i Daily intraday hydrological situation assessment index values ​​during scheduling period k Converted to intraday hydrological situation assessment index level The intraday hydrological situation assessment index x i Daily intraday hydrological situation assessment index values ​​for the baseline period s Converted to intraday hydrological situation assessment index level

[0018] Step S5.2, since j = 1, 2, ..., m, therefore, m intraday hydrological situation assessment index levels are obtained. Daily hydrological situation assessment index levels By performing pairwise comparisons and based on the matrix creation rules, an intraday hydrological situation assessment index x corresponding to the scheduling period k is established. i The Hasse matrix is ​​represented as: matrix It is an m*m matrix;

[0019] Similarly, the daily intraday hydrological situation assessment index levels for the evaluation baseline period s. By performing pairwise comparisons and based on the matrix creation rules, an intraday hydrological situation assessment index x corresponding to the evaluation baseline period s is established. i The Hasse matrix is ​​represented as: matrix It is an m*m matrix;

[0020] matrix sum matrix The m elements on the main diagonal are all 0;

[0021] Step S5.3, for the matrix After correction, the measured values ​​of the daily reservoir scheduling characteristic data for scheduling period k are incorporated to obtain the corrected matrix, which is represented as follows:

[0022] For matrix After correction, the measured values ​​of daily reservoir operation characteristic data from the evaluation baseline period s are incorporated to obtain the corrected matrix, which is represented as follows:

[0023] Step S5.4: Calculate the corrected matrix. and the corrected matrix The distance between them is the intraday hydrological situation assessment index x. i Changes in hydrological morphology during the scheduling period k

[0024] Step S6: Using formula (3), the intraday hydrological situation assessment index x is obtained. i The overall change in the intraday hydrological situation during the scheduling period k

[0025] Daily hydrological situation assessment indicators x i The overall change in the intraday hydrological situation during the scheduling period k The impact of water, wind, and solar power on the hydrological situation after their operation is reflected.

[0026] The preferred intraday hydrological situation assessment index refers to an index that quantifies the degree of intraday hydrological situation change based on hourly data of inflow and outflow from the hydropower station.

[0027] Preferably, the intraday hydrological situation assessment indicators include the hourly variation difference of hydrological situation, the intraday variation 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.

[0028] Preferably, in step S2, for each scheduling period k, k = 1, 2, ..., N, the corresponding evaluation baseline period s for which watershed water, wind, and solar power do not participate in operation is determined, specifically as follows:

[0029] The average annual water inflow during the time period k is statistically analyzed to determine the characteristics of the time period k, including the flood season, normal water season, and dry season.

[0030] The average value of multiple historical benchmark periods that are concurrent with scheduling period k, have the same characteristics, and whose watershed water, wind and solar power are not involved in the operation are obtained as the evaluation benchmark period s.

[0031] Preferably, in step S4, the method for determining the RVA threshold range is as follows:

[0032] Obtain intraday hydrological situation assessment indicators x i The 25%-75% quantiles of multiple historical data points from the same period in scheduling period k are used as the RVA threshold range RF. i .

[0033] Preferably, in step S5.1, the intraday hydrological situation assessment index x is adjusted according to the level conversion rules. i Daily intraday hydrological situation assessment index values ​​during scheduling period k Converted to intraday hydrological situation assessment index level Specifically:

[0034] If the intraday hydrological situation assessment index value In RF i The range of [0, 25%) indicates the level of the intraday hydrological situation assessment index. =1;

[0035] If the intraday hydrological situation assessment index value In RF i The range of [25%, 75%) indicates the level of the intraday hydrological situation assessment index. It is 2;

[0036] If the intraday hydrological situation assessment index value In RF i The range of [75%, 100%] indicates the level of the intraday hydrological situation assessment index. The value is 3.

[0037] Preferably, in step S5.2, according to the matrix creation rules, the intraday hydrological situation assessment index x corresponding to the scheduling period k is established. i The Hasse matrix is ​​represented as: matrix Specifically:

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

[0039] in: The index representing the intraday hydrological situation assessment is x i The intraday hydrological situation assessment index level for the l-th day of the scheduling period k, where l = 1, 2, ..., m;

[0040] Representative matrix The value of the element in the j-th row and l-th column;

[0041] Therefore, the matrix The m elements on the main diagonal all have a value of 0.

[0042] Preferably, in step S5.3, the matrix is... After correction, the measured values ​​of the daily reservoir scheduling characteristic data for scheduling period k are incorporated to obtain the corrected matrix, which is represented as follows: Specifically:

[0043] For matrix The m elements of the main diagonal are modified, where the j-th element of the main diagonal is... The correction method is as follows:

[0044] make Where Feature(j) represents the reservoir scheduling characteristics data on day j of scheduling period k, specifically three-dimensional feature data, including the intraday difference index of the reservoir's intraday outflow on day j of scheduling period k. Instantaneous indicators of reservoirs The maximum daily water level fluctuation of the reservoir

[0045] Preferably, step S5.4 specifically includes:

[0046] Step S5.4.1, the corrected matrix For an m*m matrix, its elements are re-represented as: Represents the corrected matrix The element in the j-th row and l-th column of the array, where j = 1, 2, ..., m; l = 1, 2, ..., m;

[0047] Corrected matrix For an m*m matrix, its elements are re-represented as: Represents the corrected matrix The element in the j-th row and l-th column of the array, where j = 1, 2, ..., m; l = 1, 2, ..., m;

[0048] Step S5.4.2: Using formula (5), the corrected matrix is ​​calculated. and the corrected matrix Distance between off-diagonal elements:

[0049] Step S5.4.3: Using formula (6), the corrected matrix is ​​calculated. and the corrected matrix Distance between diagonal elements:

[0050] Step S5.4.4: Using formula (7), the corrected matrix is ​​calculated. and the corrected matrix The combined distance between them is

[0051] in: The intraday hydrological situation assessment index x is located in [0,1]. i The weighted parameter of k during the scheduling period.

[0052] Preferred weighted parameters The method for determining it is as follows:

[0053] In step S5.4.4.1, formulas (8) and (9) are used to calculate the intraday hydrological situation assessment index x for scheduling period k. i Diagonal element proportion and overall proportion

[0054] In step S5.4.4.2, formula (10) is used to obtain the intraday hydrological situation assessment index x for scheduling period k. i Information entropy

[0055] In step S5.4.4.3, the weighted parameters are obtained using formula (11).

[0056] This step is now complete.

[0057] The method for assessing intraday hydrological changes in an integrated watershed operation based on water, wind, and solar power, provided by this invention, has the following advantages:

[0058] This invention provides a method for assessing intraday hydrological changes during integrated water, wind, and solar power operations. It can effectively quantify the impact of different integrated water, wind, and solar power scheduling and operation conditions on intraday hydrological changes, thereby providing a scientific basis for ecological scheduling of hydropower stations participating in integrated water, wind, and solar power operations and realizing coordinated development of hydropower operation and ecological protection. Attached Figure Description

[0059] Figure 1 is a flowchart of an assessment method for intraday hydrological situation changes in a watershed integrated operation of water, wind and solar power provided by the present invention. Detailed Implementation

[0060] To make the technical problems solved, the technical solutions, and the beneficial effects of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and are not intended to limit the invention.

[0061] This invention provides a method for assessing intraday hydrological changes during integrated operation of water, wind, and solar power in a river basin. This method can comprehensively and accurately assess the impact of hydropower station scheduling on intraday river hydrological conditions, providing a scientific basis for ecological scheduling of hydropower stations and contributing to the harmonious development of hydropower stations and the ecological environment.

[0062] Referring to Figure 1, the present invention provides a method for assessing intraday hydrological situation changes in an integrated watershed operation involving water, wind, and solar power, comprising the following steps:

[0063] Step S1: Determine the intraday hydrological situation assessment index vector X = {x1, x2, ..., x...} n};x1,x2,...,x n These represent intraday hydrological situation assessment index 1, intraday hydrological situation assessment index 2, ..., intraday hydrological situation assessment index n, respectively.

[0064] In this invention, the intraday hydrological situation assessment index refers to an index that quantifies the degree of intraday hydrological situation change based on hourly data of inflow and outflow from a hydropower station. For example, intraday hydrological situation assessment indicators include hourly variation difference of hydrological situation, intraday variation difference, daily average runoff, time of flow increase, time of flow decrease, duration of flow increase, and duration of flow decrease. This invention does not limit the specific intraday hydrological situation assessment index used.

[0065] Specifically, intraday hydrological situation assessment indicators can include indicators such as intraday hydrological situation magnitude, rate of change, duration, frequency, and timing. Among them, magnitude indicators include: hourly variation difference, intraday variation difference, and daily average runoff; rate of change indicators include: average flow increase rate, average flow decrease rate, maximum flow increase rate, maximum flow decrease rate, and instantaneous index; duration indicators include: times of flow increase and times of flow decrease; timing indicators include: duration of flow increase and duration of flow decrease; and frequency indicators include: frequency of flow reversal, frequency of flow increase, and frequency of flow decrease.

[0066] Furthermore, referring to Table 1, the indicators and their meanings are as follows:

[0067] (1) Scale indicator. This indicator mainly assesses the difference between high and low flow rates during the day, including the daily average runoff, the difference between the maximum and minimum flow rates during the day, and the maximum difference in runoff between consecutive hours.

[0068] (2) Rate of change index. This index mainly assesses the dynamic characteristics of runoff during the day, including the rate of increase in flow, the rate of decrease in flow, and the instantaneous index.

[0069] (3) Duration index. This index is used to assess the duration of the increase or decrease in runoff during the day caused by the water storage and release scheduling of a hydropower station.

[0070] (4) Frequency index. This index assesses the frequency of runoff increase or decrease during the day. It is calculated by dividing the time of runoff increase or decrease by the total time of day (24h). In addition, the flow reversal frequency is defined as the number of times runoff changes between high and low flow.

[0071] (5) Timing indicator. Because the water storage and release scheduling of hydropower stations alters the intraday distribution of runoff, the timing of runoff increases and decreases also changes. Therefore, this indicator is used to assess the moment when runoff increases to its maximum or decreases to its minimum.

[0072] Table 1. Daily Hydrological Situation Assessment Indicators and Their Meanings

[0073] Step S2: Combining the watershed hydrological characteristics and reservoir scheduling methods, determine the scheduling period after N types of watershed water, wind and solar power participate in operation. For each scheduling period k, k = 1, 2, ..., N, determine the corresponding evaluation benchmark period s for watershed water, wind and solar power not participating in operation. Both the scheduling period k and the evaluation benchmark period s have m days.

[0074] In this step, different scheduling periods are divided according to the changing trends of the watershed hydrological situation and the reservoir scheduling rules.

[0075] In this step, for each scheduling period k, k = 1, 2, ..., N, the corresponding evaluation baseline period s for which watershed water, wind, and solar power do not participate in operation is determined, specifically as follows:

[0076] The average annual water inflow during the time period k is statistically analyzed to determine the characteristics of the time period k, including the flood season, normal water season, and dry season.

[0077] The average value of multiple historical benchmark periods that are concurrent with scheduling period k, have the same characteristics, and whose watershed water, wind and solar power are not involved in the operation are obtained as the evaluation benchmark period s.

[0078] Step S3: Collect the intraday hydrological situation assessment index x i The daily intraday hydrological situation assessment index values ​​for the scheduling period k are expressed as follows: Where i = 1, 2, ..., n, j = 1, 2, ..., m, the specific meaning is: intraday hydrological situation assessment index x i The intraday hydrological situation assessment index value on the j-th day of the scheduling period k;

[0079] Obtain intraday hydrological situation assessment indicators x i The daily intraday hydrological situation assessment index values ​​for the baseline period s are expressed as:

[0080] Step S4: Using formula (1), the intraday hydrological situation assessment index x is obtained. i Changes in hydrological conditions during the scheduling period k

[0081] in: The index representing the intraday hydrological situation assessment is x i Daily intraday hydrological situation assessment index values ​​during scheduling period k The number of items within the RVA threshold range;

[0082] The index representing the intraday hydrological situation assessment is x i Daily intraday hydrological situation assessment index values ​​during scheduling period k The expected number within the RVA threshold range; Determined by formula (2):

[0083] Where: p k For intraday hydrological situation assessment indicators x i Daily intraday hydrological situation assessment index values ​​during scheduling period k The expected proportion within the RVA threshold range;

[0084] In this step, the method for determining the RVA threshold range can be: obtaining the intraday hydrological situation assessment index x i The 25%-75% quantiles of multiple historical data points from the same period in scheduling period k are used as the RVA threshold range RF. i .

[0085] Step S5: Evaluate and obtain the intraday hydrological situation assessment index x i Changes in hydrological morphology during the scheduling period k

[0086] Step S6: Using formula (3), the intraday hydrological situation assessment index x is obtained. i The overall change in the intraday hydrological situation during the scheduling period k

[0087] Daily hydrological situation assessment indicators x i The overall change in the intraday hydrological situation during the scheduling period k The impact of water, wind, and solar power on the hydrological situation after their operation is reflected.

[0088] In this invention, step S5 involves evaluating and obtaining the intraday hydrological situation assessment index x. iChanges in hydrological morphology during the scheduling period k This is an important innovative step in the present invention, and its specific implementation steps are as follows:

[0089] Step S5.1, according to the level conversion rules, adjust the intraday hydrological situation assessment index x i Daily intraday hydrological situation assessment index values ​​during scheduling period k Converted to intraday hydrological situation assessment index level The intraday hydrological situation assessment index x i Daily intraday hydrological situation assessment index values ​​for the baseline period s Converted to intraday hydrological situation assessment index level

[0090] Daily hydrological situation assessment indicators x i Daily intraday hydrological situation assessment index values ​​during scheduling period k Converted to intraday hydrological situation assessment index level The method, and the intraday hydrological situation assessment index x i Daily intraday hydrological situation assessment index values ​​for the baseline period s Converted to intraday hydrological situation assessment index level The method and principle are the same; therefore, only the intraday hydrological situation assessment index x is used. i Daily intraday hydrological situation assessment index values ​​during scheduling period k Converted to intraday hydrological situation assessment index level For example, the conversion is performed in the following way:

[0091] If the intraday hydrological situation assessment index value In RF i The range of [0, 25%) indicates the level of the intraday hydrological situation assessment index. =1;

[0092] If the intraday hydrological situation assessment index value In RF i The range of [25%, 75%) indicates the level of the intraday hydrological situation assessment index. It is 2;

[0093] If the intraday hydrological situation assessment index value In RF i The range of [75%, 100%] indicates the level of the intraday hydrological situation assessment index. The value is 3.

[0094] Therefore, by using the above methods, the intraday hydrological situation assessment index values ​​are... Convert to level 1, 2, or 3, corresponding to low, medium, and high levels.

[0095] Step S5.2, since j = 1, 2, ..., m, therefore, m intraday hydrological situation assessment index levels are obtained. Daily hydrological situation assessment index levels By performing pairwise comparisons and based on the matrix creation rules, an intraday hydrological situation assessment index x corresponding to the scheduling period k is established. i The Hasse matrix is ​​represented as: matrix It is an m*m matrix;

[0096] Similarly, the daily intraday hydrological situation assessment index levels for the evaluation baseline period s. By performing pairwise comparisons and based on the matrix creation rules, an intraday hydrological situation assessment index x corresponding to the evaluation baseline period s is established. i The Hasse matrix is ​​represented as: matrix It is an m*m matrix; matrix sum matrix The m elements on the main diagonal are all 0.

[0097] In this step, the matrix sum matrix The creation principle is the same, only using a matrix. For example, its creation method is as follows:

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

[0099] in: The index representing the intraday hydrological situation assessment is x i The intraday hydrological situation assessment index level for the l-th day of the scheduling period k, where l = 1, 2, ..., m; Representative matrix The value of the element in the j-th row and l-th column;

[0100] Therefore, the matrix The m elements on the main diagonal have a value of 0, while the elements off the diagonal have a value of +1 or -1.

[0101] Step S5.3, for the matrix After correction, the measured values ​​of the daily reservoir scheduling characteristic data for scheduling period k are incorporated to obtain the corrected matrix, which is represented as follows:

[0102] For matrix After correction, the measured values ​​of daily reservoir operation characteristic data from the evaluation baseline period s are incorporated to obtain the corrected matrix, which is represented as follows:

[0103] In this step, the matrix Perform corrections and adjust the matrix The principle for making corrections is the same, only the matrix is ​​modified. Taking the correction method as an example, the correction is made in the following way:

[0104] For matrix The m elements of the main diagonal are modified, where the j-th element of the main diagonal is... The correction method is as follows:

[0105] make Where Feature(j) represents the reservoir scheduling characteristics data on day j of scheduling period k, specifically three-dimensional feature data, including the intraday difference index of the reservoir's intraday outflow on day j of scheduling period k. Instantaneous indicators of reservoirs The maximum daily water level fluctuation of the reservoir

[0106] Therefore, in this invention, the uncorrected matrix sum matrix All are antisymmetric matrices, consisting only of 0 and ±1, with all elements on the main diagonal being 0. By adding the reservoir's operational characteristics data to the main diagonal of the matrix, the correlation between intraday hydrological situation assessment indicators and reservoir operational conditions is enhanced.

[0107] Intra-day difference index of reservoir intra-day outflow on day j of scheduling period k Instantaneous indicators of reservoirs The maximum daily water level fluctuation of the reservoir The specific meaning is:

[0108] Intraday Difference Indicator The difference between the maximum and minimum flow rates on day j of scheduling period k;

[0109] Maximum daily water level fluctuation in the reservoir It is the difference between the maximum and minimum reservoir water levels on the j-th day of the scheduling period k;

[0110] Instantaneous indicators of reservoirs It is the instantaneous index after Z-score normalization; specifically, in order to smooth short-term abnormal fluctuations in traffic and capture the trend of traffic changes, the weighted moving average method (EWMA) is used to calculate the instantaneous index, and Z-score normalization is used to normalize the instantaneous index. The calculation formula is as follows: I(t)=α·ΔQ(t)+(1-α)·I(t-1)

[0111] Where: ΔQ(t) is the flow difference over a continuous period; I(t) is the instantaneous exponent at time t; I(t-1) is the instantaneous exponent at time t-1; α is the smoothing coefficient; I norm (t) is the instantaneous exponent after Z-score normalization; μ I It is the average of the instantaneous index at all times during the day; σ I It is the standard deviation of the instantaneous exponent.

[0112] Furthermore, in order to calculate the Hasse matrix, the reservoir scheduling condition feature data is normalized to be located in the interval [0,1]. The normalization method is as follows: determine the maximum value of the m reservoir scheduling condition feature data on the main diagonal, and divide the reservoir scheduling condition feature data of each j-th day on the main diagonal by the maximum value to achieve normalization.

[0113] Step S5.4: Calculate the corrected matrix. and the corrected matrix The distance between them is the intraday hydrological situation assessment index x. i Changes in hydrological morphology during the scheduling period k

[0114] Step S5.4 specifically includes:

[0115] Step S5.4.1, the corrected matrix For an m*m matrix, its elements are re-represented as: Represents the corrected matrix The element in the j-th row and l-th column of the array, where j = 1, 2, ..., m; l = 1, 2, ..., m;

[0116] Corrected matrix For an m*m matrix, its elements are re-represented as: Represents the corrected matrix The element in the j-th row and l-th column of the array, where j = 1, 2, ..., m; l = 1, 2, ..., m;

[0117] Step S5.4.2: Using formula (5), the corrected matrix is ​​calculated. and the corrected matrix Distance between off-diagonal elements:

[0118] Step S5.4.3: Using formula (6), the corrected matrix is ​​calculated. and the corrected matrix Distance between diagonal elements:

[0119] Step S5.4.4: Using formula (7), the corrected matrix is ​​calculated. and the corrected matrix The combined distance between them is

[0120] in: The intraday hydrological situation assessment index x is located in [0,1]. i The weighted parameter of k during the scheduling period.

[0121] To accurately quantify changes in hydrological conditions during reservoir operation and comprehensively reflect changes in hydrological condition assessment indicators across different days, the entropy weight method is used to dynamically calculate the weights of diagonal and off-diagonal elements within each operation period k. This comprehensively evaluates the Hasse distances of the diagonal and off-diagonal elements. Specifically, the weighting parameters... The method for determining it is as follows:

[0122] Step S5.4.4.1, Entropy weight calculation: Since the diagonal and off-diagonal elements of the Hasse matrix are all in the interval [0,1], no standardization is required. Formulas (8) and (9) are directly used to calculate the intraday hydrological situation assessment index x for scheduling period k. i Diagonal element proportion and overall proportion

[0123] Step S5.4.4.2, calculate information entropy: using formula (10), obtain the intraday hydrological situation assessment index x for scheduling period k. i Information entropy

[0124] In step S5.4.4.3, the weighted parameters are obtained using formula (11).

[0125] This step is now complete.

[0126] This invention provides a method for assessing intraday hydrological changes in an integrated watershed operation involving water, wind, and solar power. The underlying concept can be summarized as follows:

[0127] Step 1: Collect hourly data on inflow and outflow of hydropower stations, water levels, etc., during each scheduling period before and after the integrated operation of water, wind and solar power, and statistically analyze the intraday changes in reservoir operation status and watershed hydrological conditions.

[0128] Step 2: Define intraday hydrological situation assessment indicators based on the watershed's eco-hydrological needs;

[0129] Step 3: Based on the hydrological characteristics of the watershed and the reservoir operation mode, divide the reservoir operation scheduling period and determine the corresponding evaluation benchmark period;

[0130] Based on the defined intraday hydrological situation assessment index, the RVA method was used to assess the degree of change in hydrological situation during the scheduling period.

[0131] Step 4: Based on the range of changes in the daily hydrological situation assessment indicators during the scheduling period determined by the RVA method, classify the hydrological situation change levels, thereby converting the daily daily hydrological situation assessment indicator values ​​of the scheduling period and the evaluation baseline period into daily hydrological situation assessment indicator levels.

[0132] Step 5: Compare the daily intraday hydrological situation assessment index levels for each day during the scheduling period pairwise to construct the Hasse matrix corresponding to the scheduling period; similarly, compare the daily intraday hydrological situation assessment index levels for each day during the evaluation baseline period pairwise to construct the Hasse matrix corresponding to the evaluation baseline period.

[0133] The Hasse matrices corresponding to the scheduling period and the evaluation benchmark period are corrected by adding the reservoir scheduling condition feature data of the corresponding day to the diagonal elements of the matrices, so as to obtain the corrected Hasse matrices corresponding to the scheduling period and the corrected Hasse matrices corresponding to the evaluation benchmark period.

[0134] Step 6: Combine the entropy weight method to calculate the distance between the Hasse matrix corresponding to the modified scheduling period and the Hasse matrix corresponding to the modified evaluation baseline period, and obtain the degree of hydrological morphological change for each intraday hydrological situation assessment index.

[0135] Step 7: Combine the degree of change in hydrological situation and the degree of change in hydrological morphology of the intraday hydrological situation assessment indicators to obtain the comprehensive degree of change in intraday hydrological situation for each intraday hydrological situation assessment indicator in each scheduling period.

[0136] This invention provides a method for assessing intraday hydrological changes during integrated water, wind, and solar power operations. It can effectively quantify the impact of different integrated water, wind, and solar power scheduling and operation conditions on intraday hydrological changes, thereby providing a scientific basis for ecological scheduling of hydropower stations participating in integrated water, wind, and solar power operations and realizing coordinated development of hydropower operation and ecological protection.

[0137] The above description is only a preferred embodiment of the present invention. It should be noted that for those skilled in the art, several improvements and modifications can be made without departing from the principle of the present invention, and these improvements and modifications should also be considered within the scope of protection of the present invention.

Claims

1. A method for assessing intraday hydrological situation changes in a watershed integrated water, wind, and solar power operation, characterized in that, Includes the following steps: Step S1: Determine the intraday hydrological situation assessment index vector X = {x1, x2, ..., x...} n };x1,x2,...,x n These represent intraday hydrological situation assessment index 1, intraday hydrological situation assessment index 2, ..., intraday hydrological situation assessment index n, respectively. Step S2: Combining the watershed hydrological characteristics and reservoir scheduling methods, determine the scheduling period after N types of watershed water, wind and solar power participate in operation. For each scheduling period k, k = 1, 2, ..., N, determine the corresponding evaluation benchmark period s for watershed water, wind and solar power not participating in operation. Both the scheduling period k and the evaluation benchmark period s have m days. Step S3: Collect the intraday hydrological situation assessment index x i The daily intraday hydrological situation assessment index values ​​for the scheduling period k are expressed as follows: Where i = 1, 2, ..., n, j = 1, 2, ..., m, the specific meaning is: intraday hydrological situation assessment index x i The intraday hydrological situation assessment index value on the j-th day of the scheduling period k; Obtain intraday hydrological situation assessment indicators x i The daily intraday hydrological situation assessment index values ​​for the baseline period s are expressed as: Step S4: Using formula (1), the intraday hydrological situation assessment index x is obtained. i Changes in hydrological conditions during the scheduling period k in: The index representing the intraday hydrological situation assessment is x i Daily intraday hydrological situation assessment index values ​​during scheduling period k The number of items within the RVA threshold range; The index representing the intraday hydrological situation assessment is x i Daily intraday hydrological situation assessment index values ​​during scheduling period k The expected number within the RVA threshold range; Determined by formula (2): Where: p k For intraday hydrological situation assessment indicators x i Daily intraday hydrological situation assessment index values ​​during scheduling period k The expected proportion within the RVA threshold range; Step S5: Evaluate and obtain the intraday hydrological situation assessment index x i Changes in hydrological morphology during the scheduling period k Step S5.1, according to the level conversion rules, adjust the intraday hydrological situation assessment index x i Daily intraday hydrological situation assessment index values ​​during scheduling period k Converted to intraday hydrological situation assessment index level The intraday hydrological situation assessment index x i Daily intraday hydrological situation assessment index values ​​for the baseline period s Converted to intraday hydrological situation assessment index level Step S5.2, since j = 1, 2, ..., m, therefore, m intraday hydrological situation assessment index levels are obtained. Daily hydrological situation assessment index levels By performing pairwise comparisons and based on the matrix creation rules, an intraday hydrological situation assessment index x corresponding to the scheduling period k is established. i The Hasse matrix is ​​represented as: matrix It is an m*m matrix; Similarly, the daily intraday hydrological situation assessment index levels for the evaluation baseline period s. By performing pairwise comparisons and based on the matrix creation rules, an intraday hydrological situation assessment index x corresponding to the evaluation baseline period s is established. i The Hasse matrix is ​​represented as: matrix It is an m*m matrix; matrix sum matrix The m elements on the main diagonal are all 0; Step S5.3, for the matrix After correction, the measured values ​​of the daily reservoir scheduling characteristic data for scheduling period k are incorporated to obtain the corrected matrix, which is represented as follows: For matrix After correction, the measured values ​​of daily reservoir operation characteristic data from the evaluation baseline period s are incorporated to obtain the corrected matrix, which is represented as follows: Step S5.4: Calculate the corrected matrix. and the corrected matrix The distance between them is the intraday hydrological situation assessment index x. i Changes in hydrological morphology during the scheduling period k Step S6: Using formula (3), the intraday hydrological situation assessment index x is obtained. i The overall change in the intraday hydrological situation during the scheduling period k Daily hydrological situation assessment indicators x i The overall change in the intraday hydrological situation during the scheduling period k The impact of water, wind, and solar power on the hydrological situation after their operation is reflected.

2. The method for assessing intraday hydrological situation changes in an integrated watershed operation based on water, wind, and solar power, as described in claim 1, is characterized in that... Intraday hydrological situation assessment indicators refer to indicators that quantify the degree of intraday hydrological situation change based on hourly data of inflow and outflow from hydropower stations.

3. The method for assessing intraday hydrological situation changes in a watershed integrated water, wind, and solar power system according to claim 1, characterized in that, Intraday hydrological condition assessment indicators include hourly variation of hydrological conditions, intraday variation, daily average runoff, times of flow increase, times of flow decrease, duration of flow increase, and duration of flow decrease.

4. The method for assessing intraday hydrological situation changes in a watershed integrated water, wind, and solar power system according to claim 1, characterized in that, In step S2, for each scheduling period k, k = 1, 2, ..., N, determine the corresponding evaluation baseline period s for which watershed water, wind, and solar power do not participate in operation, specifically as follows: The average annual water inflow during the time period k is statistically analyzed to determine the characteristics of the time period k, including the flood season, normal water season, and dry season. The average value of multiple historical benchmark periods that are concurrent with scheduling period k, have the same characteristics, and whose watershed water, wind and solar power are not involved in the operation are obtained as the evaluation benchmark period s.

5. The method for assessing intraday hydrological situation changes in an integrated watershed operation based on water, wind, and solar power, as described in claim 1, is characterized in that... In step S4, the method for determining the RVA threshold range is as follows: Obtain intraday hydrological situation assessment indicators x i The 25%-75% quantiles of multiple historical data points from the same period in scheduling period k are used as the RVA threshold range RF. i .

6. The method for assessing intraday hydrological situation changes in an integrated watershed operation based on water, wind, and solar power, as described in claim 5, is characterized in that... In step S5.1, according to the level conversion rules, the intraday hydrological situation assessment index x is... i Daily intraday hydrological situation assessment index values ​​during scheduling period k Converted to intraday hydrological situation assessment index level Specifically: If the intraday hydrological situation assessment index value In RF i The range of [0, 25%) indicates the level of the intraday hydrological situation assessment index. =1; If the intraday hydrological situation assessment index value In RF i The range of [25%, 75%) indicates the level of the intraday hydrological situation assessment index. It is 2; If the intraday hydrological situation assessment index value In RF i The range of [75%, 100%] indicates the level of the intraday hydrological situation assessment index. The value is 3.

7. The method for assessing intraday hydrological situation changes in an integrated watershed operation based on water, wind, and solar power, as described in claim 1, is characterized in that... In step S5.2, according to the matrix creation rules, the intraday hydrological situation assessment index x corresponding to the scheduling period k is established. i The Hasse matrix is ​​represented as: matrix Specifically: According to formula (4), determine the matrix. The value of each element: in: The index representing the intraday hydrological situation assessment is x i The intraday hydrological situation assessment index level for the l-th day of the scheduling period k, where l = 1, 2, ..., m; Representative matrix The value of the element in the j-th row and l-th column; Therefore, the matrix The m elements on the main diagonal all have a value of 0.

8. The method for assessing intraday hydrological situation changes in a watershed integrated water, wind, and solar power system according to claim 1, characterized in that, In step S5.3, the matrix is... After correction, the measured values ​​of the daily reservoir scheduling characteristic data for scheduling period k are incorporated to obtain the corrected matrix, which is represented as follows: Specifically: For matrix The m elements of the main diagonal are modified, where the j-th element of the main diagonal is... The correction method is as follows: make Where Feature(j) represents the reservoir scheduling characteristics data on day j of scheduling period k, specifically three-dimensional feature data, including the intraday difference index of the reservoir's intraday outflow on day j of scheduling period k. Instantaneous indicators of reservoirs The maximum daily water level fluctuation of the reservoir 9. The method for assessing intraday hydrological situation changes in an integrated watershed operation based on water, wind, and solar power, as described in claim 1, is characterized in that... Step S5.4 specifically includes: Step S5.4.1, the corrected matrix For an m*m matrix, its elements are re-represented as: Represents the corrected matrix The element in the j-th row and l-th column of the array, where j = 1, 2, ..., m; l = 1, 2, ..., m; Corrected matrix For an m*m matrix, its elements are re-represented as: Represents the corrected matrix The element in the j-th row and l-th column of the array, where j = 1, 2, ..., m; l = 1, 2, ..., m; Step S5.4.2: Using formula (5), the corrected matrix is ​​calculated. and the corrected matrix Distance between off-diagonal elements: Step S5.4.3: Using formula (6), the corrected matrix is ​​calculated. and the corrected matrix Distance between diagonal elements: Step S5.4.4: Using formula (7), the corrected matrix is ​​calculated. and the corrected matrix The combined distance between them is in: The intraday hydrological situation assessment index x is located in [0,1]. i The weighted parameter of k during the scheduling period.

10. The method for assessing intraday hydrological situation changes in an integrated watershed operation based on water, wind, and solar power, as described in claim 9, is characterized in that... Weighted parameters The method for determining it is as follows: In step S5.4.4.1, formulas (8) and (9) are used to calculate the intraday hydrological situation assessment index x for scheduling period k. i Diagonal element proportion and overall proportion In step S5.4.4.2, formula (10) is used to obtain the intraday hydrological situation assessment index x for scheduling period k. i Information entropy In step S5.4.4.3, the weighted parameters are obtained using formula (11). This step is now complete.

Images