A refined scheduling modeling method based on the topological relationship of graded sections of a water-wind-solar-storage integrated base

Abstract

The present invention provides a refined scheduling modeling method based on the hierarchical cross-section topological relationship of a water-wind-solar-storage integrated base, comprising the following steps: step S1, the division of the scheduling units of the water-wind-solar-storage integrated base; step S2, establishing the scheduling unit coupling relationship of the water-wind-solar-storage integrated base; step S3, based on the scheduling unit coupling relationship of the water-wind-solar-storage integrated base, realizing the refined modeling and solution of the water-wind-solar-storage integrated scheduling, and obtaining the output of each scheduling unit in each time period t during the scheduling period. The present invention provides a refined scheduling modeling method based on the hierarchical cross-section topological relationship of a water-wind-solar-storage integrated base, which accurately analyzes the coupling relationship between large-scale, multi-type power sources and their transmission sections in the water-wind-solar-storage integrated base, and effectively supports the refined scheduling of water-wind-solar-storage integrated. The use of the present invention is not limited to the power source composition, category, quantity and power collection method of the water-wind-solar-storage integrated base, and has strong adaptability and can be widely promoted and applied.

Description

Technical Field

[0001] The present invention belongs to the technical field of clean energy multi-energy complementary scheduling, and specifically relates to a refined scheduling modeling method based on the hierarchical section topological relationship of a water-wind-solar-storage integrated base. Background Art

[0002] The water-wind-solar-storage integrated base is formed by connecting different types of power sources in series through transmission lines of different voltage levels and gathering them into a unified whole. The changes in the operating mode of each power source affect other related power sources through the transmission lines, thereby affecting the overall operating efficiency of the water-wind-solar-storage integrated base. Large-scale new energy clusters in water-wind-solar-storage integrated bases are generally built in remote areas. They need to be stepped up and connected to the 220kV and 500kV main grids through 35kV and 110kV transmission lines or sent out together with other power sources such as large and medium-sized hydropower and energy storage facilities. Therefore, the power access situation in the water-wind-solar-storage integrated base is very complicated. The purpose of the operation and dispatch of the water-wind-solar-storage integrated base is to consider the operating characteristics of all power sources and the coupling relationship between power sources, and to find the operating mode of the power source group that maximizes the base efficiency under the given constraints of transmission lines at all levels. However, the existing dispatch modeling methods for integrated water, wind, solar and storage bases ignore the complex coupling relationship between power sources, and often generalize the integrated base power access mode, network topology and related transmission constraints. It is difficult to effectively reflect the transmission relationship between the output influence of various types of power sources and the impact of the power transmission capacity of various levels of transmission lines on the operation mode of power sources. The water, wind, solar and storage dispatch schemes formulated based on this are of poor practicality and difficult to meet the requirements of refined dispatching. From the perspective of existing technical methods, there is still a lack of effective methods that can accurately describe the coupling relationship and output constraint relationship between multi-category, large-scale power sources and transmission sections of various voltage levels in integrated bases, which has become a key problem to be solved in the multi-energy complementary design and operation dispatch simulation of integrated bases. Summary of the invention

[0003] In view of the defects of the prior art, the present invention provides a refined scheduling modeling method based on the hierarchical section topological relationship of an integrated water, wind, solar and storage base, which can effectively solve the above problems.

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

[0005] The present invention provides a refined scheduling modeling method based on the topological relationship of the hierarchical sections of a water-wind-solar-storage integrated base, comprising the following steps:

[0006] Step S1, division of dispatching units of the integrated water, wind, solar and storage base:

[0007] According to the power source type and power output characteristics of the water-wind-solar-storage integrated base, the hydropower station cluster of the water-wind-solar-storage integrated base is divided into multiple hydropower dispatching units, the wind power station cluster is divided into multiple wind power dispatching units, the photovoltaic power station cluster is divided into multiple photovoltaic dispatching units, and the energy storage power station cluster is divided into multiple energy storage dispatching units;

[0008] Step S2, establishing the scheduling unit coupling relationship of the integrated water, wind, solar and storage base:

[0009] Step S2.1, partitioning the dispatching unit of the integrated water, wind, solar and storage base to form power supply partitions of different levels, each power supply partition having a corresponding level of transmission sections for accessing a higher adjacent level, thereby forming transmission sections of different levels;

[0010] Step S2.2, according to the level of each transmission section and the hierarchical relationship of the associated transmission sections of each level, encode each transmission section to clarify the hierarchical relationship of each transmission section;

[0011] Step S2.3, according to the superior-subordinate relationship of each transmission section, a directed graph is used to construct a transmission section topology relationship;

[0012] Step S2.4, according to the topological relationship of the transmission sections, in order from low level to high level, a depth-first search method is used to determine the calculation order of each transmission section;

[0013] Step S2.5, according to the access mode of each dispatching unit, directly associate it with the corresponding transmission section as the dispatching unit directly associated with the transmission section, so as to construct the dispatching unit coupling relationship of the integrated water, wind, solar and storage base;

[0014] Step S3, based on the coupling relationship of the dispatching units of the integrated water, wind, solar and storage base, implement the refined modeling and solution of the integrated water, wind, solar and storage dispatching, and obtain the output of each dispatching unit in each time period t during the dispatching period, including:

[0015] Step S3.1, determining the output of each hydropower dispatching unit in each time period t during the dispatching period and the output of each energy storage dispatching unit in each time period t during the dispatching period;

[0016] Step S3.2, determining the output of each wind power dispatching unit in each time period t during the dispatching period and the output of each photovoltaic dispatching unit in each time period t during the dispatching period.

[0017] Preferably, in step S1, the method for dividing the hydropower dispatching units is:

[0018] The division is made based on the regulation and storage capacity of the hydropower station and the hydraulic and power connections between cascade hydropower stations, with a single hydropower station or a group of hydropower stations serving as a hydropower dispatching unit. Specifically, a hydropower station with a regulation capacity of daily or weekly or more forms a hydropower dispatching unit on its own; a run-of-river hydropower station and its directly adjacent upstream hydropower stations with the same power consumption direction form a hydropower dispatching unit together.

[0019] Preferably, in step S1, the division method of the wind power dispatching unit and the photovoltaic dispatching unit is the same, and the wind power dispatching unit and the photovoltaic dispatching unit are both called new energy dispatching units, and the division method is:

[0020] Taking into account the output characteristics, output correlation, grid access nodes, wind speed and time-varying characteristics of new energy power stations, the actual output process or simulated output process of each new energy power station is used as the feature input, and the Spearman rank correlation coefficient is used as the output correlation analysis and evaluation index of the new energy power station. The Canopy-K-means comprehensive clustering algorithm combining Canopy clustering and K-means clustering is used to aggregate a single new energy power station to a new energy power station cluster, and a new energy power station cluster is obtained. Each new energy power station cluster forms a new energy dispatching unit. The specific aggregation process is as follows:

[0021] Step S1.1: Rough clustering of new energy power station clusters using the Canopy algorithm

[0022] Step S1.1.1, establish a data set X of the new energy power station output cluster = {X 1 ,X 2 ,...,X n ,...,X N}; where N represents the number of new energy power stations participating in the cluster analysis; X n represents the output sequence of the new energy power station n, n=1,2,…,N; X n ={x 1 ,x 2 ,...,x t ,…,x T}, where x 1 ,x 2 ,...,x t ,…,x T represents the output of the new energy power station n in each time period, and T represents the number of time periods;

[0023] Step S1.1.2, using the cross-validation method to determine the first distance threshold T of the output correlation of the new energy power station 1 and the second distance threshold T 2 , and T 1 >T 2 ;

[0024] Step S1.1.3, randomly select the output sequence X of the new energy power station a from the data set X of the new energy power station output cluster a , a=1,2,…,N, with the output sequence X of new energy power station a a As the center of the canopy class, calculate the output sequence X of the new energy power station a respectively a The output sequence X of each other new energy power station b b Output correlation distance The output correlation distance with the new energy power station a Less than the second distance threshold T 2 All new energy power stations are added to the canopy class where the center is located, marked with a strong mark and deleted from the data set X; the output correlation distance with new energy power station a is Located at the first distance threshold T 1 and the second distance threshold T 2 All new energy power stations between are added to the canopy class where the center is located and marked with a weak mark. The new energy power stations marked with a weak mark will participate in the next round of clustering process and become the center or member of the new canopy class; the output correlation distance with the new energy power station a will be Greater than the first distance threshold T 1 All new energy power stations are individually used as the center of a new canopy class to carry out cluster analysis; this process is repeated until all new energy power stations are assigned to one or more canopy classes, thereby achieving coarse clustering of new energy power station clusters and dividing them into K coarse clusters;

[0025] in: The following formula is used to obtain:

[0026]

[0027] in: and They are the output sequence X of new energy power station a a and the output sequence X of the new energy power station b b The rank after transformation; and They are the output sequence X of new energy power station a a and the output sequence X of the new energy power station b b The standard deviation of is the output sequence X of the new energy power station a a and the output sequence X of the new energy power station b b The Spearman rank correlation coefficient ranges from [-1,1]. The larger the absolute value, the better the correlation. represent and The covariance of

[0028] Step S1.2, based on the K coarse clusters obtained by the Canopy algorithm, considering the output correlation distance between different renewable energy power stations, each coarse cluster is further clustered and analyzed by the K-means algorithm to obtain multiple subsets of renewable energy power stations, each of which forms a renewable energy dispatching unit.

[0029] Preferably, in step S1, the energy storage scheduling unit is divided as follows:

[0030] The energy storage power stations in the integrated base are divided according to their operating characteristics, storage capacity and cycle efficiency, with a single energy storage power station or a group of energy storage power stations serving as the dispatching unit. Specifically, for new energy storage power stations, if the energy storage equipment form, operating characteristics, cycle efficiency and power access point are the same, an energy storage dispatching unit is formed. If the combination conditions are not met, a single new energy storage power station serves as the energy storage dispatching unit. For pumped-storage power stations, due to differences in the development form, upper and lower reservoir capacity conditions, head conditions, unit output characteristics and power station operation mode of pumped-storage power stations, a single pumped-storage power station serves as the energy storage dispatching unit.

[0031] Preferably, in step S2.1, the specific method of power supply partitioning and transmission section classification is:

[0032] According to the spatial distribution of each dispatching unit, the power collection method and the voltage level of the power transmission channel, power partitioning is carried out to form recursive power partitioning and hierarchical transmission sections;

[0033] Specifically, the dispatching units with similar geographical locations, the same power input points, and directly associated with the same power transmission channel, as well as the power supply zones of the next level directly associated with the power transmission channel, are merged into the power supply zone of this level, and the power transmission channel of the power supply zone of this level is defined as the transmission section of the corresponding level;

[0034] Specifically, the transmission section with the highest voltage level of the integrated water, wind, solar and storage base is defined as the first-level transmission section R1; the number of the first-level transmission sections R1 is one or more; for each first-level transmission section R1, a first-level power partition C1 is formed; each of the first-level power partitions C1 includes the first-level transmission section R1 and the first-level dispatching unit U1 directly associated with the power collection end of the first-level transmission section R1, and the second-level power partition C2;

[0035] Each secondary power supply partition C2 includes a secondary transmission section R2 directly connected to the power collection end of the primary transmission section R1, a secondary dispatching unit U2 directly associated with the power collection end of the secondary transmission section R2, and a tertiary power supply partition C3;

[0036] By analogy, a multi-level cascaded power supply partition and a multi-level transmission section are formed, forming a coupling structure of power supply partition and transmission section classification; wherein each transmission section can be directly associated with one or more transmission sections of the next adjacent level, but has only one directly associated transmission section of the previous adjacent level;

[0037] In step S2.2, each transmission section is coded as follows:

[0038] The first-level transmission section is coded as: 1A; where 1 represents the level of the transmission section; A represents the sequence number of the first-level transmission section among all the first-level transmission sections;

[0039] The secondary transmission section is coded as: 1A-2B, which means: the Bth secondary transmission section connected to the power collection end of the primary transmission section 1A;

[0040] The tertiary transmission section is coded as: 1A-2B-3C, which means: the Cth tertiary transmission section connected to the power collection end of the secondary transmission section 1A-2B.

[0041] Preferably, step S2.3 is specifically:

[0042] The vertices of the directed graph represent the transmission sections of each level, the edges of the directed graph represent the direct connection between the transmission sections of adjacent levels, and the direction of the edge of the directed graph represents the superior-subordinate relationship between two transmission sections with direct connection, that is, from the lower transmission section to the upper transmission section, so as to construct a directed graph by combining the superior-subordinate relationship of each transmission section and form the topological relationship of the transmission section;

[0043] Step S2.4 is specifically as follows:

[0044] From the topological relationship of the transmission sections, locate the highest level transmission section, use it as the root s of the directed graph, and put it in the target sorting queue line;

[0045] Then, in the transmission section topology relationship, determine the transmission section set U of the direct subordinate of the root s, traverse the transmission sections u in the transmission section set U in turn, and add each traversed transmission section u to the target sorting queue line in order;

[0046] For each transmission section u in the transmission section set U, locate the transmission section set directly below it, and add each traversed transmission section to the target sorting queue line in order; this cycle continues until all transmission sections are searched;

[0047] The reverse order of the transmission sections in the target sorting queue line is the calculation order of each transmission section determined.

[0048] Preferably, in step S3.1, determining the output of each hydropower dispatching unit in each time period t during the dispatching period and the output of each energy storage dispatching unit in each time period t during the dispatching period specifically includes:

[0049] Step S3.1.1, estimate the average output of each hydropower dispatching unit during the dispatching period:

[0050] According to the coupling relationship of the dispatching units of the integrated water, wind, solar and storage base, the topological relationship of the hydropower dispatching units is determined, so that the hydropower dispatching units are sorted according to the upstream and downstream hydraulic connections; starting from the upstream, the average downstream flow of each hydropower dispatching unit during the dispatching period is calculated in sequence, and the average output of each hydropower dispatching unit during the dispatching period is estimated through the refined allocation method of water-determined current flow;

[0051] Step S3.1.2, generating the initial output process of each hydropower dispatching unit in each time period t during the dispatching period:

[0052] Taking the average output, expected output and forced output of each hydropower dispatching unit in the dispatching period as constraints, the successive load shedding method is adopted. On the basis of the base transmission capacity curve after deducting the designed output of the wind and solar power stations, the power balance calculation is performed for each hydropower dispatching unit separately to obtain the initial output process of each hydropower dispatching unit in each time period t during the dispatching period;

[0053] Step S3.1.3, correct the initial output process of each hydropower dispatching unit in each time period t during the dispatching period, and obtain the corrected initial output process of each hydropower dispatching unit in each time period t during the dispatching period:

[0054] The unit combination vibration zone analysis method is used to determine whether the initial output process of the hydropower dispatching unit in each time period t during the dispatching period meets the stable operation requirements. If there is an output in the time period t that does not meet the requirements, the output in the time period t that does not meet the requirements is corrected to the stable operation boundary, thereby obtaining the corrected initial output process of the hydropower dispatching unit in each time period t during the dispatching period;

[0055] Step S3.1.4, taking the initial output process of the hydropower dispatching unit in each time period t during the dispatching period obtained in step S3.1.3 as the initial solution, taking the output of the hydropower dispatching unit in each time period t during the dispatching period as the decision variable, taking the minimum fluctuation of the output process of the hydro-wind-solar-storage multi-energy complementary power generation system as the dispatching target F, and establishing a dispatching model;

[0056] Among them: The scheduling target F is:

[0057]

[0058] Where: T is the number of time periods in the scheduling period; P tis the output of the hydropower, wind power, photovoltaic and energy storage multi-energy complementary power generation system in time period t, is the sum of the outputs of the hydropower dispatching unit, wind power dispatching unit, photovoltaic dispatching unit and energy storage dispatching unit in time period t; The average output of the hydropower, wind, solar and storage multi-energy complementary power generation system during the dispatch period;

[0059] Step S3.1.5, using the bee colony algorithm to solve and optimize the dispatching model established in step S3.1.4. In the optimization process, the wind and solar power output adopts the wind and solar power design output as a constant value; the extraction state and output of the energy storage dispatching unit are determined according to the hydropower output and wind and solar power output of the hydro-wind-solar-storage multi-energy complementary power generation system, and controlled by the energy storage balance constraint;

[0060] The heuristic constraint processing mechanism is used to modify the output of the hydropower dispatching unit in each time period t during the dispatching period and the extraction state and output of the energy storage dispatching unit in each time period t during the dispatching period in turn, so as to ensure the feasibility of the solution obtained by the algorithm optimization.

[0061] Thus, the output of the hydropower dispatching unit in each time period t during the dispatching period and the extraction state and output of the energy storage dispatching unit in each time period t during the dispatching period are finally solved.

[0062] Preferably, in step S3.1.5, a heuristic constraint processing mechanism is used to sequentially correct the output of the hydropower dispatching unit and the output of the energy storage dispatching unit in each time period t of the dispatching period, including:

[0063] Step S3.1.5.1, for each hydropower dispatching unit i obtained by the bee colony algorithm in this optimization, the output of each time period t during the dispatching period is expressed as The power balance constraint strategy is used for correction to obtain the corrected output of hydropower dispatching unit i in each time period t during the dispatching period. The specific method is:

[0064] ①Through the formula Get the full-time power generation of hydropower dispatching unit i during the dispatching period Where Δt is the length of each time period t;

[0065] ②If Among them, E i is the theoretical power generation of hydropower dispatching unit i during the dispatching period; ε is the deviation margin, then the power balance constraint is corrected, otherwise, no correction is made;

[0066] The correction method is:

[0067] Using formula The power difference Divide into T equal parts, T is the number of scheduling periods, and deltE is the power value of each equal part;

[0068] like That is, the total power generation of hydropower dispatching unit i Greater than the theoretical power generation E i , then set P i h,min is the forced output of hydropower dispatching unit i; if not, the full-time power generation of hydropower dispatching unit i Less than the theoretical power generation E i , then set P i h,max is the expected output limit of hydropower dispatching unit i; the correction is repeated in this cycle until the full-time power generation of hydropower dispatching unit i is And the theoretical power generation E i The difference is controlled within the deviation margin ε;

[0069] Step S3.1.5.2: for the corrected output of hydropower dispatching unit i in each time period t during the dispatching period Adopt the hydropower final water level control constraint repair strategy to Further correction is made to obtain the corrected output of hydropower dispatching unit i in each time period t during the dispatching period:

[0070] Specifically, the final water level of the hydropower dispatching unit i in each time period t during the dispatching period and the final water level of the dispatching period are obtained by using the refined allocation method of power-determined water output; if the final water level of the hydropower dispatching unit i in time period t does not meet the requirements, the output of the time period is re-determined by reverse calculation through flow control or water level control method;

[0071] The specific method is:

[0072] ① If the final water level of hydropower dispatching unit i in each time period t during the dispatching period is greater than the given final water level, it indicates that the power generation is insufficient and a certain power step length ΔE is required. 0 Increase output Get the adjusted output of hydropower dispatching unit i in each time period t during the dispatching period The method is:

[0073] Set the charge step length to ΔE 0 The time period is divided into T parts, and the length of each part is Adjust each period from period 1. If restricted, a single period cannot be If the power is completely consumed, the remaining power will be allocated to its adjacent time periods;

[0074] ② According to the adjusted output of hydropower dispatching unit i in each time period t during the dispatching period Recalculate the final water level of hydropower dispatching unit i in each time period t during the dispatching period by using electricity to determine water, return to step ①, and repeat this cycle until the final water level requirement is met. Finally, the corrected output of hydropower dispatching unit i in each time period t during the dispatching period is obtained, which is expressed as

[0075] ③ If the final water level of hydropower dispatching unit i in each time period t during the dispatching period is lower than the given final water level, the output will be reduced in each time period for correction;

[0076] Step S3.1.5.3, adopt the heuristic constraint processing mechanism to correct the output of the energy storage dispatching unit in each time period t during the dispatching period:

[0077] ① Determine the working status of the energy storage dispatching unit:

[0078] For the kth energy storage dispatching unit put into operation, its working state is jointly determined by the output process of the water, wind, and light dispatching units in the water-wind-solar-storage multi-energy complementary power generation system and the first k-1 energy storage dispatching units put into operation: if the current water-wind-solar-storage multi-energy complementary power generation system has power abandonment in time period t, and the energy storage dispatching unit k still has energy storage space, then the energy storage dispatching unit k is in the energy storage state in time period t; if there is idle capacity in the sending section of the water-wind-solar-storage multi-energy complementary power generation system, and the energy storage dispatching unit k still has the discharge capability, then the energy storage dispatching unit k is in the power generation state in time period t;

[0079] ② Generate energy storage constraint corridor:

[0080] The storage energy E of energy storage scheduling unit k at the beginning of the scheduling period k,0 Based on the forward sequence calculation method, the energy storage constraint corridor is generated, and the formula is:

[0081]

[0082] in: and are the upper limit of energy storage output and power generation output of energy storage dispatching unit k in time period t respectively; E k,t,max1 and E k,t,min1 , are respectively the upper and lower limits of the energy storage capacity of the energy storage dispatching unit k at the end of time period t during the positive sequence calculation; E k,t-1,max1 and E k,t-1,min1 , are the upper and lower limits of the energy storage capacity of the energy storage dispatching unit k at the end of time period t-1 during the positive sequence calculation; E k,max and E k,min , are the upper and lower limits of the designed storage capacity of the energy storage dispatching unit k respectively; is the cycle efficiency of energy storage dispatch unit k; E k,0,max1 and E k,0,min1, are respectively the upper and lower limits of the energy storage capacity of the energy storage dispatching unit k at the beginning of the dispatching period when the forward sequence calculation method is adopted;

[0083] The storage energy E of energy storage dispatch unit k at the end of the dispatch period k,T Based on the reverse calculation method, the energy storage constraint corridor is generated, and the formula is:

[0084]

[0085] Where: E k,t,max2 and E k,t,min2 , are the upper and lower limits of the energy storage capacity of the energy storage dispatching unit k at the end of time period t when calculating in reverse order; E k,t+1,max2 and E k,t+1,min2 , are the upper and lower limits of the energy storage capacity of the energy storage dispatching unit k at the end of time period t+1 when calculating in reverse order; E k,T,max2 and E k,T,min2 , are respectively the upper and lower limits of the energy storage capacity of the energy storage dispatch unit k at the end of the dispatch period when the reverse calculation method is adopted;

[0086] Further take the envelope to generate the energy storage constraint corridor (E k,t,min ,E k,t,max ):

[0087]

[0088] E k,t,max and E k,t,min , are the upper and lower limits of the energy storage capacity of the energy storage scheduling unit k in time period t, respectively;

[0089] ③ According to the energy storage capacity constraint corridor (E k,t,min ,E k,t,max ), the extraction state and output of the energy storage dispatching unit k in time period t obtained by the bee colony algorithm in step S3.1.5 are corrected.

[0090] Preferably, step S3.2, determining the output of each wind power dispatching unit in each time period t during the dispatching period and the output of each photovoltaic dispatching unit in each time period t during the dispatching period, specifically includes:

[0091] Step S3.2.1, in the hydropower, wind, photovoltaic and storage multi-energy complementary power generation system, the output of each hydropower dispatching unit in each time period t during the dispatching period and the output of each energy storage dispatching unit in each time period t during the dispatching period have been solved, and the initial output of the wind power dispatching unit and the photovoltaic dispatching unit is given as the design output;

[0092] Step S3.2.2, traversing each transmission section in sequence according to the calculation order of the transmission sections determined in step S2.4;

[0093] Step S3.2.3: for the current traversed transmission section, it is represented as transmission section ts j , and its corresponding power partition is represented as power partition j; according to the output of the hydropower dispatching unit in each time period t during the dispatching period and the output of the energy storage dispatching unit in each time period t during the dispatching period, as well as the latest output of the wind power dispatching unit and the latest output of the photovoltaic dispatching unit, the power partition output recursive function is used to calculate the output P of power partition j in time period t Sj,t :

[0094]

[0095] in: is the recursive function of power partition j in time period t;

[0096] They are respectively the transmission section ts with power supply partition j in time period t j The sum of the outputs of all directly related hydropower dispatching units, the sum of the outputs of all energy storage dispatching units, the sum of the outputs of all wind power dispatching units, and the sum of the outputs of all photovoltaic dispatching units; D Hj , D Pj , D Qj , D Rj , respectively, the transmission section ts with power supply partition j j The number set of directly related hydropower dispatching units, the number set of energy storage dispatching units, the number set of wind power dispatching units and the number set of photovoltaic dispatching units; m, k, q, r are the numbers of the hydropower dispatching unit, the energy storage dispatching unit, the wind power dispatching unit and the photovoltaic dispatching unit respectively; They are the output of the hydropower dispatching unit m in time period t, the output of the energy storage dispatching unit k in time period t, the latest output of the wind power dispatching unit q in time period t, and the latest output of the photovoltaic dispatching unit r in time period t;

[0097] D is the sum of the outputs of all the next-level power partitions in power partition j after traversal and constraint processing in time period t; Sj Represents the set of sub-level power partition numbers, where l is the sub-level power partition number; Represents the recursive function of power partition l in time period t; represents the output upper limit of the next-level power partition l in time period t;

[0098] Step S3.2.4, determine the output P of power partition j in time period t Sj,t , whether it meets the transmission section ts j The power transmission restriction constraint function is: Sj,t ≤N S,j , where NS,j is the transmission section ts of power supply zone j j If the conditions are met, return to step S3.2.3 and continue to traverse the next transmission section; if the conditions are not met, execute step S3.2.5;

[0099] Step S3.2.5, combining the network topology of transmission sections at all levels and the power transmission restrictions of the transmission sections of each power partition, correct the output of all wind power dispatching units and all photovoltaic dispatching units in power partition j to obtain the latest output of each wind power dispatching unit in time period t and the latest output of each photovoltaic dispatching unit in time period t; then return to step S3.2.3 and continue to traverse the next transmission section.

[0100] Preferably, step S3.2.5 is specifically:

[0101] Step S3.2.5.1, based on the output P of power partition j in time period t Sj,t and transmission section ts j Power supply limit N S,j The difference between the transmission section ts j The curtailment of power in period t

[0102] Step S3.2.5.2: If power partition j does not contain the next-level power partition, the abandoned power Allocate according to the proportion of online electricity of the wind and solar dispatching units directly associated with power partition j, thereby correcting the output of the corresponding wind and solar dispatching units in time period t, and updating the online electricity of the wind and solar dispatching units in time period t; wherein the wind and solar dispatching units include wind power dispatching units and photovoltaic dispatching units;

[0103] If power partition j contains the next-level power partition, the power will be discarded. The online power consumption of all wind and solar dispatching units directly associated with power partition j and indirectly associated with the next-level power partition is allocated to correct the output of the corresponding wind and solar dispatching unit in time period t and update the online power consumption of the wind and solar dispatching unit in time period t.

[0104] The refined scheduling modeling method based on the topological relationship of the hierarchical sections of the integrated water, wind, solar and storage base provided by the present invention has the following advantages:

[0105] The present invention provides a refined scheduling modeling method based on the topological relationship of the hierarchical sections of the water-wind-solar-storage integrated base, which accurately analyzes the coupling relationship between large-scale, multi-type power sources and their transmission sections in the water-wind-solar-storage integrated base, and effectively supports the refined scheduling of water-wind-solar-storage integrated base. The use of the present invention is not limited to the power source composition, category, quantity and power collection method of the water-wind-solar-storage integrated base, and has strong adaptability and can be widely promoted and applied. BRIEF DESCRIPTION OF THE DRAWINGS

[0106] Figure 1 A flow chart of a refined scheduling modeling method based on the topological relationship of graded sections of a water-wind-solar-storage integrated base provided by the present invention;

[0107] Figure 2 This is a schematic diagram of the hierarchical sections and power supply zones of the water-wind-solar integrated base of the present invention;

[0108] Figure 3 It is a directed graph of topological relations of graded sections of the integrated base of the present invention;

[0109] Figure 4 This is the power supply access method of the integrated water, wind, solar and storage base in the basin of the present invention;

[0110] Figure 5 This is a topological diagram of the graded sections and power zoning of the integrated water, wind, solar and storage base in the case basin of the present invention. DETAILED DESCRIPTION

[0111] In order to make the purpose, technical solutions and advantages of the present invention more clearly understood, the present invention is 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 intended to limit the present invention. In addition, the technical features involved in the various embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

[0112] The present invention aims at the fact that the existing integrated water, wind, solar and storage scheduling modeling methods cannot accurately describe the coupling relationship and output constraint relationship between the multi-category, large-scale power sources and transmission sections of various voltage levels in the integrated base, and it is difficult to effectively solve the integrated water, wind, solar and storage scheduling modeling problem affected by the complex power collection method of the base power source. A refined scheduling modeling method based on the hierarchical section topological relationship of the integrated water, wind, solar and storage base is proposed to accurately analyze the coupling relationship between the large-scale, multi-type power sources and their transmission sections in the integrated water, wind, solar and storage base, and effectively support the refined scheduling of integrated water, wind, solar and storage. The use of the present invention is not limited to the power source composition, category, quantity and power collection method of the integrated water, wind, solar and storage base, and it has strong adaptability and can be widely promoted and applied.

[0113] To achieve the above object, the present invention provides a refined scheduling modeling method based on the topological relationship of the hierarchical sections of the integrated water, wind, solar and storage base, which mainly includes the following three steps:

[0114] Step S1, division of dispatching units of the integrated water, wind, solar and storage base;

[0115] Step S2, establishing a scheduling unit coupling relationship of the integrated water, wind, solar and storage base;

[0116] Step S3, based on the coupling relationship of the dispatching units of the integrated water, wind, solar and storage base, realize the refined modeling and solution of the integrated water, wind, solar and storage dispatching.

[0117] The following are the three steps in detail:

[0118] Step S1, division of dispatching units of the integrated water, wind, solar and storage base:

[0119] According to the power source type and power output characteristics of the integrated water, wind, solar and storage base, the hydropower station cluster of the integrated water, wind, solar and storage base is divided into multiple hydropower dispatching units, the wind power station cluster is divided into multiple wind power dispatching units, the photovoltaic power station cluster is divided into multiple photovoltaic dispatching units, and the energy storage power station cluster is divided into multiple energy storage dispatching units.

[0120] The specific division methods of hydropower dispatching unit, wind power and photovoltaic dispatching unit, and energy storage dispatching unit include:

[0121] (1) Division method of hydropower dispatching units

[0122] The division is made based on the regulation and storage capacity of the hydropower station and the complex hydraulic and power connections between cascade hydropower stations, with a single hydropower station or a group of hydropower stations serving as a hydropower dispatching unit. Specifically, a hydropower station with a daily or weekly regulation capacity can form a hydropower dispatching unit independently, and its access method in the base and the power transmission limit of the transmission channel will be set according to the original hydropower settings. Run-of-river hydropower stations can be bundled with their directly adjacent upstream hydropower stations with the same power consumption direction to form a hydropower dispatching unit together. Its access method in the base will be set according to the original regulated hydropower, and the power transmission limit of the transmission channel will be set according to the total installed capacity of the hydropower dispatching unit.

[0123] (2) Division of wind power dispatching units and photovoltaic dispatching units

[0124] The division method of wind power dispatching unit and photovoltaic dispatching unit is the same. Both wind power dispatching unit and photovoltaic dispatching unit are called new energy dispatching unit. The division method is as follows:

[0125] The new energy dispatching unit is mainly divided in the form of new energy power station clusters. In specific operations, the output characteristics, output correlation, grid access nodes, time-varying characteristics of wind speed and light intensity of new energy power stations are comprehensively considered. The actual output process or simulated output process of each new energy power station is used as the feature input, and the Spearman rank correlation coefficient is used as the output correlation analysis and evaluation index of the new energy power station. The Canopy-K-means comprehensive clustering algorithm combining Canopy clustering and K-means clustering is used to aggregate a single new energy power station to a new energy power station cluster, and a new energy power station cluster is obtained. Each new energy power station cluster forms a new energy dispatching unit. The specific aggregation process includes two steps:

[0126] Step S1.1, performing rough clustering of new energy power station clusters using the Canopy algorithm;

[0127] Step S1.2, based on the K coarse clusters obtained by the Canopy algorithm, considering the output correlation distance between different renewable energy power stations, each coarse cluster is further clustered and analyzed by the K-means algorithm to obtain multiple subsets of renewable energy power stations, each of which forms a renewable energy dispatching unit.

[0128] The specific implementation process of step S1.1 and step S1.2 is as follows:

[0129] Step S1.1, perform rough clustering of new energy power station clusters using the Canopy algorithm:

[0130] Step S1.1.1, establish a data set X of the new energy power station output cluster = {X 1 ,X 2 ,...,X n ,...,X N}; where N represents the number of new energy power stations participating in the cluster analysis; X n represents the output sequence of the new energy power station n, n=1,2,…,N; X n ={x 1 ,x 2 ,...,x t ,…,x T}, where x 1 ,x 2 ,...,x t ,…,x T represents the output of the new energy power station n in each time period, and T represents the number of time periods;

[0131] Step S1.1.2, using the cross-validation method to determine the first distance threshold T of the output correlation of the new energy power station 1 and the second distance threshold T 2 , and T 1 >T2 ;

[0132] Step S1.1.3, randomly select the output sequence X of the new energy power station a from the data set X of the new energy power station output cluster a , a=1,2,…,N, with the output sequence X of new energy power station a a As the center of the canopy class, calculate the output sequence X of the new energy power station a respectively a The output sequence X of each other new energy power station b b Output correlation distance b=1,2,…,N, the output correlation distance with the new energy power station a Less than the second distance threshold T 2 All new energy power stations are added to the canopy class where the center is located, marked with a strong mark and deleted from the data set X; the output correlation distance with new energy power station a is Located at the first distance threshold T 1 and the second distance threshold T 2 All new energy power stations between are added to the canopy class where the center is located and marked with a weak mark. The new energy power stations marked with a weak mark will participate in the next round of clustering process and become the center or member of the new canopy class; the output correlation distance with the new energy power station a will be Greater than the first distance threshold T 1 All new energy power stations are individually used as the center of a new canopy class to carry out cluster analysis; this process is repeated until all new energy power stations are assigned to one or more canopy classes, thereby achieving coarse clustering of new energy power station clusters and dividing them into K coarse clusters;

[0133] Among them: Output correlation distance The Spearman rank correlation coefficient is calculated as follows:

[0134]

[0135] in: and They are the output sequence X of new energy power station a a and the output sequence X of the new energy power station b b The rank after transformation; and They are the output sequence X of new energy power station a a and the output sequence X of the new energy power station b b The standard deviation of is the output sequence X of the new energy power station a a and the output sequence X of the new energy power station bb The Spearman rank correlation coefficient ranges from [-1,1]. The larger the absolute value, the better the correlation. represent and The covariance of Indicates the output correlation distance, and its value range is [-1,0].

[0136] Step S1.2, based on the K coarse clusters obtained by the Canopy algorithm, considering the output correlation distance between different renewable energy power stations, each coarse cluster is further clustered and analyzed by the K-means algorithm to obtain multiple subsets of renewable energy power stations, each of which forms a renewable energy dispatching unit.

[0137] After determining the wind / solar new energy dispatching unit, the single wind power or photovoltaic installed capacity in the wind / solar new energy dispatching unit is used as the weight, and the long series output simulation time series process of each wind power and photovoltaic is used as the basis. Combined with the spatial distribution relationship between power stations in the wind / solar new energy dispatching unit, and considering a certain time difference offset, the equivalent output time series process of the wind / solar new energy dispatching unit is obtained by superposing the output process time series, which represents the overall output characteristics of the wind / solar new energy dispatching unit.

[0138] (3) Division method of energy storage dispatching units

[0139] The division method of energy storage dispatching units is as follows: the division is based on the operating characteristics, storage capacity, cycle efficiency and other factors of each energy storage power station in the integrated base, and a single energy storage power station or a group of energy storage power stations is used as the dispatching unit; specifically, for new energy storage power stations, if the energy storage equipment form, operating characteristics, cycle efficiency and power access point are the same, the energy storage dispatching unit can be formed in the form of a group of energy storage power stations; the power transmission limit of the transmission channel is set according to the total installed capacity of the energy storage dispatching unit; if the combination conditions are not met, a single new energy storage power station is used as the energy storage dispatching unit; for pumped storage power stations, due to differences in the development form, upper and lower reservoir capacity conditions, head conditions, unit output characteristics, power station operation mode, etc. of pumped storage power stations, it is advisable to use a single pumped storage power station as the energy storage dispatching unit and no longer combine them.

[0140] Step S2, establishing the scheduling unit coupling relationship of the integrated water, wind, solar and storage base:

[0141] The coupling relationship between the integrated base dispatching units is described by adopting the method of "power source zoning and transmission section classification + establishment of the topological relationship of the classified transmission section", which mainly includes five steps:

[0142] Step S2.1, partitioning the dispatching unit of the integrated water, wind, solar and storage base to form power supply partitions of different levels, each power supply partition having a corresponding level of transmission sections for accessing a higher adjacent level, thereby forming transmission sections of different levels;

[0143] Step S2.2, according to the level of each transmission section and the hierarchical relationship of the associated transmission sections of each level, encode each transmission section to clarify the hierarchical relationship of each transmission section;

[0144] Step S2.3, according to the superior-subordinate relationship of each transmission section, a directed graph is used to construct a transmission section topology relationship;

[0145] Step S2.4, according to the topological relationship of the transmission sections, in order from low level to high level, a depth-first search method is used to determine the calculation order of each transmission section;

[0146] Step S2.5, according to the access mode of each dispatching unit, directly associate it with the corresponding transmission section as the dispatching unit directly associated with the transmission section, thereby constructing the dispatching unit coupling relationship of the integrated water, wind, solar and storage base.

[0147] The specific implementation of steps S2.1 to S2.5 is as follows:

[0148] Step S2.1, power supply zoning and transmission section classification:

[0149] The dispatching units of the integrated water, wind, solar and storage base are divided into different levels of power supply zones. Each level of power supply zone has a corresponding level of transmission sections for accessing higher adjacent levels, thus forming different levels of transmission sections;

[0150] In this step, the specific methods of power supply zoning and transmission section classification are as follows:

[0151] Power zoning is carried out according to the spatial distribution of each dispatching unit of the integrated base, the power collection method and the voltage level of the power transmission channel, forming a recursive power zoning and hierarchical transmission section;

[0152] Specifically, the dispatching units with similar geographical locations, the same power input points, and directly associated with the same power transmission channel, as well as the power supply zones of the next level directly associated with the power transmission channel, are merged into the power supply zone of this level, and the power transmission channel of the power supply zone of this level is defined as the transmission section of the corresponding level;

[0153] Specifically, the transmission section with the highest voltage level of the integrated water, wind, solar and storage base is defined as the first-level transmission section R1; the number of the first-level transmission sections R1 is one or more; for each first-level transmission section R1, a first-level power partition C1 is formed; each of the first-level power partitions C1 includes the first-level transmission section R1 and the first-level dispatching unit U1 directly associated with the power collection end of the first-level transmission section R1, and the second-level power partition C2;

[0154] Each secondary power supply partition C2 includes a secondary transmission section R2 directly connected to the power collection end of the primary transmission section R1, a secondary dispatching unit U2 directly associated with the power collection end of the secondary transmission section R2, and a tertiary power supply partition C3;

[0155] By analogy, a multi-level cascaded power partition and a multi-level transmission section are formed, forming a coupled structure of power partition and transmission section classification, which can be called a recursive partition and classification structure; wherein each transmission section can be directly associated with one or more transmission sections of the next adjacent level, but has only one directly associated transmission section of the previous adjacent level;

[0156] Step S2.2, encode each transmission section according to the level of each transmission section and the hierarchical relationship of its associated transmission sections of each level, clarify the hierarchical relationship of each transmission section, and according to the transmission section coding, the level of the transmission section, the hierarchical relationship of the transmission section and the number of transmission sections of each level can be identified.

[0157] The specific method for encoding each transmission section is as follows:

[0158] The first-level transmission section is coded as: 1A; where 1 represents the level of the transmission section; A represents the sequence number of the first-level transmission section among all the first-level transmission sections;

[0159] The secondary transmission section is coded as: 1A-2B, which means: the Bth secondary transmission section connected to the power collection end of the primary transmission section 1A;

[0160] The tertiary transmission section is coded as: 1A-2B-3C, which means: the Cth tertiary transmission section connected to the power collection end of the secondary transmission section 1A-2B.

[0161] Therefore, the coding principle of transmission section is:

[0162] 1) The transmission section numbers of each level are marked by the number of the directly related upper level transmission section + the level and serial number of the current transmission section, with "-" in the middle. That is, the first level transmission section is 1X, the second level transmission section is 1X-2Y, the third level transmission section is 1X-2Y-3Z, and the codes of other levels of transmission sections are deduced in the same way; among them, X, Y and Z represent the serial numbers of the first level transmission section, the second level transmission section and the third level transmission section in the same level transmission section.

[0163] 2) Each transmission section can only have one directly associated upper transmission section, but can have multiple directly associated lower transmission sections.

[0164] by Figure 2 For example, the base includes 1 primary, 3 secondary and 1 tertiary transmission sections.

[0165] Step S2.3, according to the superior-subordinate relationship of each transmission section, a directed graph is used to construct a transmission section topology relationship;

[0166] The present invention constructs the topological relationship of the hierarchical transmission sections of the integrated water, wind, solar and storage base based on a directed graph. The specific method is: when constructing the topological relationship of the hierarchical transmission sections of the integrated water, wind, solar and storage base based on the directed graph, the vertices of the directed graph represent the transmission sections at each level, the edges of the directed graph represent the direct connection between the adjacent levels of transmission sections, and the direction of the edge of the directed graph represents the superior-subordinate relationship of two transmission sections with direct connection, that is, the direction of the edge of the directed graph is set according to the current collection direction of the transmission sections at each level, that is, from the lower transmission section to the upper transmission section, so as to combine the superior-subordinate relationship of each transmission section to construct a directed graph and form a topological relationship of the transmission section;

[0167] For example, a directed graph D(O, T) is constructed in combination with the hierarchical transmission section relationship, where O represents the starting point number and T represents the end point number. Figure 2 The description of the integrated base hierarchical transmission section in the directed graph is as follows Figure 3 shown.

[0168] Step S2.4, according to the topological relationship of the transmission sections, in order from low level to high level, a depth-first search method is used to determine the calculation order of each transmission section;

[0169] Step S2.4 is specifically as follows:

[0170] When the depth-first search method is used to determine the calculation order of the integrated base transmission section, the highest level transmission section of the entire integrated base hierarchical section topology structure is first located from the topological relationship of the transmission section, and this is used as the root s of the directed graph and placed in the target sorting queue line;

[0171] Then, in the transmission section topology relationship, determine the transmission section set U of the direct subordinate of the root s, traverse the transmission sections u in the transmission section set U in turn, and add each traversed transmission section u to the target sorting queue line in order;

[0172] For each transmission section u in the transmission section set U, locate the transmission section set directly below it, and add each traversed transmission section to the target sorting queue line in order; this cycle continues until all transmission sections are searched;

[0173] The reverse order of the transmission sections in the target sorting queue line is the calculation order of each transmission section determined.

[0174] by Figure 3 Taking the directed graph shown in FIG. 1 as an example, the sections added to the queue line are 0, 1, 2, 3, and 4 in sequence, that is, l = {0, 1, 2, 3, 4}, and the calculation order of the transmission section is 4→3→2→1→0.

[0175] Step S2.5, according to the access mode of each dispatching unit, directly associate it with the corresponding transmission section as the dispatching unit directly associated with the transmission section, so as to construct the dispatching unit coupling relationship of the integrated water, wind, solar and storage base;

[0176] Combined with the access mode of each dispatching unit in the integrated water, wind, solar and storage base, it is directly associated with the aforementioned divided hierarchical transmission section as the power partition directly associated with the hierarchical transmission section. When the hierarchical transmission section topology relationship is constructed, the dispatching unit coupling relationship can be established. The operation mode and output changes of each dispatching unit are transmitted layer by layer to other associated dispatching units through the hierarchical transmission section.

[0177] Figure 2 That is, the coupling relationship of the integrated base dispatching units after considering the topological relationship of the transmission section.

[0178] Step S3, based on the coupling relationship of the dispatching units of the integrated water, wind, solar and storage base, the refined modeling and solution of the integrated water, wind, solar and storage dispatching are realized to obtain the output of each dispatching unit in each time period t during the dispatching period:

[0179] Step S3.1, determining the output of each hydropower dispatching unit in each time period t during the dispatching period and the output of each energy storage dispatching unit in each time period t during the dispatching period;

[0180] Step S3.2, determining the output of each wind power dispatching unit in each time period t during the dispatching period and the output of each photovoltaic dispatching unit in each time period t during the dispatching period.

[0181] The specific process of step S3.1 to step S3.2 is as follows:

[0182] In step S3.1, the output of each hydropower dispatching unit in each time period t during the dispatching period and the output of each energy storage dispatching unit in each time period t during the dispatching period are determined, specifically including:

[0183] Step S3.1.1, estimate the average output of each hydropower dispatching unit during the dispatching period:

[0184] According to the coupling relationship of the dispatching units of the integrated water, wind, solar and storage base, the topological relationship of the hydropower dispatching units is determined, so that the hydropower dispatching units are sorted according to the upstream and downstream hydraulic connections; starting from the upstream, the average downstream flow of each hydropower dispatching unit during the dispatching period is calculated in sequence, and the average output of each hydropower dispatching unit during the dispatching period is estimated through the refined allocation method of water-determined current flow;

[0185] Step S3.1.2, generating the initial output process of each hydropower dispatching unit in each time period t during the dispatching period:

[0186] Taking the average output, expected output and forced output of each hydropower dispatching unit in the dispatching period as constraints, the successive load shedding method is adopted. On the basis of the base transmission capacity curve after deducting the designed output of the wind and solar power stations, the power balance calculation is performed for each hydropower dispatching unit separately to obtain the initial output process of each hydropower dispatching unit in each time period t during the dispatching period;

[0187] In order to improve the calculation efficiency of the model, when the successive load shedding method is used to determine the initial output process of the hydropower dispatching unit, the hydropower water level / flow constraint is used as a relaxation constraint and is not processed temporarily during load shedding.

[0188] Step S3.1.3, correct the initial output process of each hydropower dispatching unit in each time period t during the dispatching period, and obtain the corrected initial output process of each hydropower dispatching unit in each time period t during the dispatching period:

[0189] The unit combination vibration zone analysis method is used to determine whether the initial output process of the hydropower dispatching unit in each time period t during the dispatching period meets the stable operation requirements. If there is an output in the time period t that does not meet the requirements, the output in the time period t that does not meet the requirements is corrected to the stable operation boundary, thereby obtaining the corrected initial output process of the hydropower dispatching unit in each time period t during the dispatching period;

[0190] On the basis of the determination of the initial output process of the hydropower dispatching unit, the operation mode of the energy storage dispatching unit is analyzed by superimposing the design output of wind and solar energy. The dispatching unit stores energy when there is power abandonment in the integrated system and outputs power when the integrated system needs to top load. At the same time, based on the above work, an algorithm is formed to optimize the initial solution and evaluate the dispatching benefit corresponding to the initial solution. The specific process is shown in step S3.1.4.

[0191] Step S3.1.4, taking the initial output process of the hydropower dispatching unit in each time period t during the dispatching period obtained in step S3.1.3 as the initial solution, taking the output of the hydropower dispatching unit in each time period t during the dispatching period as the decision variable, taking the minimum fluctuation of the output process of the hydro-wind-solar-storage multi-energy complementary power generation system as the dispatching target F, and establishing a dispatching model;

[0192] Among them: The scheduling target F is:

[0193]

[0194] Where: T is the number of time periods in the scheduling period; P t is the output of the hydropower, wind power, photovoltaic and energy storage multi-energy complementary power generation system in time period t, is the sum of the outputs of the hydropower dispatching unit, wind power dispatching unit, photovoltaic dispatching unit and energy storage dispatching unit in time period t; The average output of the hydropower, wind, solar and storage multi-energy complementary power generation system during the dispatch period;

[0195] Step S3.1.5, using the bee colony algorithm to solve and optimize the dispatching model established in step S3.1.4. In the optimization process, the wind and solar power output adopts the wind and solar power design output as a constant value; the extraction state and output of the energy storage dispatching unit are determined according to the hydropower output and wind and solar power output of the hydro-wind-solar-storage multi-energy complementary power generation system, and controlled by the energy storage balance constraint;

[0196] During the optimization process, a heuristic constraint processing mechanism is used to correct the output of the hydropower dispatching unit in each time period t during the scheduling period and the extraction state and output of the energy storage dispatching unit in each time period t during the scheduling period in turn, so as to ensure the feasibility of the solution obtained by the algorithm optimization and improve the overall optimization efficiency of the algorithm.

[0197] Thus, the output of the hydropower dispatching unit in each time period t during the dispatching period and the extraction state and output of the energy storage dispatching unit in each time period t during the dispatching period are finally solved.

[0198] In step S3.1.5, the heuristic constraint processing mechanism is used to sequentially correct the output of the hydropower dispatching unit and the output of the energy storage dispatching unit in each time period t of the dispatching period, including:

[0199] Step S3.1.5.1, hydropower dispatching unit power balance constraint processing strategy:

[0200] The output of each hydropower dispatching unit i in each time period t during the dispatching period obtained by the bee colony algorithm is expressed as The power balance constraint strategy is used for correction to obtain the corrected output of hydropower dispatching unit i in each time period t during the dispatching period. The specific method is:

[0201] ①Through the formula Get the full-time power generation of hydropower dispatching unit i during the dispatching period Where Δt is the length of each time period t;

[0202] ②If Among them, E i is the theoretical power generation of hydropower dispatching unit i during the dispatching period; ε is the deviation margin, then the power balance constraint is corrected, otherwise, no correction is made;

[0203] The correction method is:

[0204] Using formula The power difference Divide into T equal parts, T is the number of scheduling periods, and deltE is the power value of each equal part;

[0205] like That is, the total power generation of hydropower dispatching unit i Greater than the theoretical power generation E i , then set P i h,min is the forced output of hydropower dispatching unit i; if not, the full-time power generation of hydropower dispatching unit i Less than the theoretical power generation E i , then set P i h,max is the expected output limit of hydropower dispatching unit i; the correction is repeated in this cycle until the full-time power generation of hydropower dispatching unit i is And the theoretical power generation E i The difference is controlled within the deviation margin ε;

[0206] Step S3.1.5.2, hydropower dispatching unit final water level control constraint repair strategy:

[0207] For the corrected output of hydropower dispatching unit i in each time period t during the dispatching period Adopt the hydropower final water level control constraint repair strategy to Further correction is made to obtain the corrected output of hydropower dispatching unit i in each time period t during the dispatching period:

[0208] Specifically, the final water level of the hydropower dispatching unit i in each time period t during the dispatching period and the final water level of the dispatching period are obtained by using the refined allocation method of power-determined water output; if the final water level of the hydropower dispatching unit i in time period t does not meet the requirements, the output of the time period is re-determined by reverse calculation through flow control or water level control method;

[0209] The specific method is:

[0210] ① If the final water level of hydropower dispatching unit i in each time period t during the dispatching period is greater than the given final water level, it indicates that the power generation is insufficient and a certain power step length ΔE is required. 0 Increase output Get the adjusted output of hydropower dispatching unit i in each time period t during the dispatching period The method is:

[0211] Set the charge step length to ΔE 0 The time period is divided into T parts, and the length of each part is Adjust each period from period 1. If restricted, a single period cannot be If the power is completely consumed, the remaining power will be allocated to its adjacent time periods;

[0212] ② According to the adjusted output of hydropower dispatching unit i in each time period t during the dispatching period Recalculate the final water level of hydropower dispatching unit i in each time period t during the dispatching period by using electricity to determine water, return to step ①, and repeat this cycle until the final water level requirement is met. Finally, the corrected output of hydropower dispatching unit i in each time period t during the dispatching period is obtained, which is expressed as

[0213] ③ If the final water level of hydropower dispatching unit i in each time period t during the dispatching period is lower than the given final water level, the output will be reduced in each time period for correction;

[0214] Step S3.1.5.3, adopt the heuristic constraint processing mechanism to correct the output of the energy storage dispatching unit in each time period t during the dispatching period:

[0215] ① Determine the working status of the energy storage dispatching unit:

[0216] For the kth energy storage dispatching unit put into operation, its working state is jointly determined by the output process of the water, wind, and light dispatching units in the water-wind-solar-storage multi-energy complementary power generation system and the first k-1 energy storage dispatching units put into operation: if the current water-wind-solar-storage multi-energy complementary power generation system (only the first k-1 energy storage dispatching units are put into operation) has power abandonment in time period t, and the energy storage dispatching unit k still has energy storage space, then the energy storage dispatching unit k is in the energy storage state in time period t; if there is idle capacity in the sending section of the water-wind-solar-storage multi-energy complementary power generation system, and the energy storage dispatching unit k still has the discharge capability, then the energy storage dispatching unit k is in the power generation state in time period t;

[0217] ② Generate energy storage constraint corridor:

[0218] The storage energy E of energy storage scheduling unit k at the beginning of the scheduling period k,0 Based on the forward sequence calculation method, the energy storage constraint corridor is generated, and the formula is:

[0219]

[0220] in: and are the upper limit of energy storage output and power generation output of energy storage dispatching unit k in time period t respectively; E k,t,max1 and E k,t,min1 , are respectively the upper and lower limits of the energy storage capacity of the energy storage dispatching unit k at the end of time period t during the positive sequence calculation; E k,t-1,max1 and E k,t-1,min1 , are the upper and lower limits of the energy storage capacity of the energy storage dispatching unit k at the end of time period t-1 during the positive sequence calculation; E k,max and E k,min , are the upper and lower limits of the designed storage capacity of the energy storage dispatching unit k respectively; is the cycle efficiency of energy storage dispatch unit k; E k,0,max1 and E k,0,min1 , are respectively the upper and lower limits of the energy storage capacity of the energy storage dispatching unit k at the beginning of the dispatching period when the forward sequence calculation method is adopted;

[0221] The storage energy E of energy storage dispatch unit k at the end of the dispatch period k,T Based on the reverse calculation method, the energy storage constraint corridor is generated, and the formula is:

[0222]

[0223] Where: E k,t,max2 and E k,t,min2 , are the upper and lower limits of the energy storage capacity of the energy storage dispatching unit k at the end of time period t when calculating in reverse order; E k,t+1,max2 and E k,t+1,min2 , are the upper and lower limits of the energy storage capacity of the energy storage dispatching unit k at the end of time period t+1 when calculating in reverse order; E k,T,max2 and E k,T,min2 , are respectively the upper and lower limits of the energy storage capacity of the energy storage dispatch unit k at the end of the dispatch period when the reverse calculation method is adopted;

[0224] Further take the envelope to generate the energy storage constraint corridor (E k,t,min ,E k,t,max ):

[0225]

[0226] E k,t,max and E k,t,min , are the upper and lower limits of the energy storage capacity of the energy storage scheduling unit k in time period t, respectively;

[0227] ③ According to the energy storage capacity constraint corridor (E k,t,min ,E k,t,max), the extraction state and output of the energy storage dispatching unit k in time period t obtained by the bee colony algorithm in step S3.1.5 are corrected.

[0228] Step S3.2, determining the output of each wind power dispatching unit in each time period t during the dispatching period and the output of each photovoltaic dispatching unit in each time period t during the dispatching period:

[0229] The idea of ​​determining the output of the wind power dispatching unit and the photovoltaic dispatching unit is: after the output of the hydropower dispatching unit and the energy storage dispatching unit is confirmed, the initial output of the wind power dispatching unit and the photovoltaic dispatching unit is used as the design output, and the initial output of the wind power dispatching unit and the photovoltaic dispatching unit is continuously revised in the following way:

[0230] In the process of dispatching the integrated water, wind, solar and storage base, the output change of each dispatching unit in the base takes into account the influence of power transmission restrictions of each level of transmission sections, and transmits the influence to other dispatching units through the hierarchical transmission section topology. The present invention adopts the power partition output recursive function to realize the transmission of the dispatching unit output influence, so as to continuously correct the output of the wind power dispatching unit and the photovoltaic dispatching unit, so that the transmission sections at all levels meet the power transmission restrictions.

[0231] The specific process is as follows:

[0232] Step S3.2.1, in the hydropower, wind, photovoltaic and storage multi-energy complementary power generation system, the output of each hydropower dispatching unit in each time period t during the dispatching period and the output of each energy storage dispatching unit in each time period t during the dispatching period have been solved, and the initial output of the wind power dispatching unit and the photovoltaic dispatching unit is given as the design output;

[0233] Step S3.2.2, traversing each transmission section in sequence according to the calculation order of the transmission sections determined in step S2.4;

[0234] Step S3.2.3: for the current traversed transmission section, it is represented as transmission section ts j , and its corresponding power partition is represented as power partition j; according to the output of the hydropower dispatching unit in each time period t during the dispatching period and the output of the energy storage dispatching unit in each time period t during the dispatching period, as well as the latest output of the wind power dispatching unit and the latest output of the photovoltaic dispatching unit, the power partition output recursive function is used to calculate the output P of power partition j in time period t Sj,t :

[0235]

[0236] in: is the recursive function of power partition j in time period t;

[0237] They are respectively the transmission section ts with power supply partition j in time period tj The sum of the outputs of all directly related hydropower dispatching units, the sum of the outputs of all energy storage dispatching units, the sum of the outputs of all wind power dispatching units, and the sum of the outputs of all photovoltaic dispatching units; D Hj , D Pj , D Qj , D Rj , respectively, the transmission section ts with power supply partition j j The number set of directly related hydropower dispatching units, the number set of energy storage dispatching units, the number set of wind power dispatching units and the number set of photovoltaic dispatching units; m, k, q, r are the numbers of the hydropower dispatching unit, the energy storage dispatching unit, the wind power dispatching unit and the photovoltaic dispatching unit respectively; They are the output of the hydropower dispatching unit m in time period t, the output of the energy storage dispatching unit k in time period t, the latest output of the wind power dispatching unit q in time period t, and the latest output of the photovoltaic dispatching unit r in time period t;

[0238] D is the sum of the outputs of all the next-level power partitions in power partition j after traversal and constraint processing in time period t; Sj Represents the set of sub-level power partition numbers, where l is the sub-level power partition number; Represents the recursive function of power partition l in time period t; represents the output upper limit of the next-level power partition l in time period t;

[0239] Step S3.2.4, determine the output P of power partition j in time period t Sj,t , whether it meets the transmission section ts j The power transmission restriction constraint function is: Sj,t ≤N S,j , where N S,j is the transmission section ts of power supply zone j j If the conditions are met, return to step S3.2.3 and continue to traverse the next transmission section; if the conditions are not met, execute step S3.2.5;

[0240] Step S3.2.5, combining the network topology of transmission sections at all levels and the power transmission restrictions of the transmission sections of each power partition, correct the output of all wind power dispatching units and all photovoltaic dispatching units in power partition j to obtain the latest output of each wind power dispatching unit in time period t and the latest output of each photovoltaic dispatching unit in time period t; then return to step S3.2.3 and continue to traverse the next transmission section.

[0241] Step S3.2.5 is specifically:

[0242] Step S3.2.5.1, based on the output P of power partition j in time period t Sj,t and transmission section ts j Power supply limit N S,j The difference between the transmission section ts j The curtailment of power in period t

[0243] Step S3.2.5.2: If power partition j does not contain the next-level power partition, the abandoned power Allocate according to the proportion of online electricity of the wind and solar dispatching units directly associated with power partition j, thereby correcting the output of the corresponding wind and solar dispatching units in time period t, and updating the online electricity of the wind and solar dispatching units in time period t; wherein the wind and solar dispatching units include wind power dispatching units and photovoltaic dispatching units;

[0244] If power partition j contains the next-level power partition, the power will be discarded. The online power consumption of all wind and solar dispatching units directly associated with power partition j and indirectly associated with the next-level power partition is allocated to correct the output of the corresponding wind and solar dispatching unit in time period t and update the online power consumption of the wind and solar dispatching unit in time period t.

[0245] A specific embodiment is described below:

[0246] (1) Taking a water-wind-solar-storage integrated base in a river basin as an example, the base dispatching units are divided into:

[0247] In order to further explain the specific implementation methods of the present invention in detail, data of a water-wind-solar-storage integrated base in a river basin in the southwest was collected. The total installed capacity of the water-wind-solar-storage integrated base in the river basin is 25.6 million kW, including 3 hydropower projects with a total installed capacity of 4.82 million kW; 16 photovoltaic power generation projects with a total installed capacity of 17.78 million kW; 1 pumped storage power station with a total installed capacity of 3 million kW. The clean electricity of the base is transmitted and consumed through a ±800kV ultra-high voltage direct current transmission project. See the base power supply access method for details. Figure 4 The power supply of the base is shown in Table 1:

[0248] Table 1. Power supply status of the integrated water, wind, solar and storage base in the case basin

[0249]

[0250] The power station cluster division method based on hierarchical clustering and the combined hierarchical iterative method were used to aggregate the photovoltaic clusters of the water-wind-solar integrated base in the basin, and four photovoltaic clusters were obtained. The photovoltaic clusters were used as photovoltaic dispatching units. The three hydropower stations involved in the integrated base all have daily regulation capabilities or above, and the access channels are not shared, so each hydropower station is used as a separate hydropower dispatching unit. The energy storage facilities involved in the base are pumped storage, which is considered as a separate energy storage dispatching unit.

[0251] (2) The coupling relationship of the dispatching unit of the case basin integration base is established by adopting the method of "power source zoning and transmission section classification + establishment of topological relationship of classified sections":

[0252] Considering the spatial distribution of each dispatching unit, the power collection method and the voltage level of the transmission line section, the case basin integration base can be abstracted into one primary and one secondary transmission section. The topology of the hierarchical sections and power partitions of the case base is shown in Figure 5 The coupling relationship of the base dispatching unit is shown in Table 2; the determined transmission section traversal order is 11-21→11.

[0253] Table 2. Coupling relationship of basin dispatching units in the case

[0254]

[0255] (3) Set the power transmission capacity of each level of transmission section, carry out integrated dispatching modeling of water, wind, solar and storage, and use the hierarchical section topology to realize the transmission of dispatching unit output influence and power source output correction in the optimization process to ensure that the power transmission of each level of transmission section meets the corresponding constraints:

[0256] The case basin contains two hierarchical transmission sections numbered "11" and "11-21". According to the above rules, the transmission capacity of the first-level transmission section numbered "11" is 8 million kW, and the transmission capacity of the first-level transmission section numbered "11-21" is 3 million kW. When using the hierarchical transmission section topology to transmit the output influence of the dispatching unit, the sum of the outputs of all dispatching units cannot exceed the transmission capacity of the "11" transmission section, and the sum of the outputs of the energy storage dispatching unit, photovoltaic dispatching unit 2, photovoltaic dispatching unit 3 and photovoltaic dispatching unit 4 cannot exceed the transmission capacity of the "11-21" transmission section. If the output adjustment of a certain dispatching unit causes the "11" or "11-21" transmission section to abandon power, the output influence will be transmitted to each dispatching unit through the transmission restriction constraints of each level of transmission sections and the recursive function of the power supply partition output.

[0257] The power topology relationship construction method of the integrated water, wind, solar and storage base proposed in the present invention is used to carry out short-term scheduling modeling of the integrated water, wind, solar and storage multi-energy complementarity. The hourly scheduling output process of the integrated water, wind, solar and storage base in the case basin on a typical day is shown in Table 3.

[0258] Table 3. Hourly dispatching output process table for a typical day in the case basin

[0259]

[0260] It can be seen from Table 3 that after being processed by the method provided by the present invention, the output of each level of transmission sections and their associated power sources all meet the relevant constraint requirements and meet the needs of refined scheduling of the integrated water, wind, solar and storage base.

[0261] The present invention provides a refined scheduling modeling method based on the topological relationship of the hierarchical sections of the water-wind-solar-storage integrated base, which accurately analyzes the coupling relationship between large-scale, multi-type power sources and their transmission sections in the water-wind-solar-storage integrated base, and effectively supports the refined scheduling of water-wind-solar-storage integrated base. The use of the present invention is not limited to the power source composition, category, quantity and power collection method of the water-wind-solar-storage integrated base, and has strong adaptability and can be widely promoted and applied.

[0262] 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 refined scheduling modeling method based on the topological relationship of the hierarchical sections of a water, wind, solar and storage integrated base, characterized in that: The following steps are involved: Step S1, division of dispatching units of the integrated water, wind, solar and storage base: According to the power source type and power output characteristics of the water-wind-solar-storage integrated base, the hydropower station cluster of the water-wind-solar-storage integrated base is divided into multiple hydropower dispatching units, the wind power station cluster is divided into multiple wind power dispatching units, the photovoltaic power station cluster is divided into multiple photovoltaic dispatching units, and the energy storage power station cluster is divided into multiple energy storage dispatching units; Step S2, establishing the scheduling unit coupling relationship of the integrated water, wind, solar and storage base: Step S2.1, partitioning the dispatching unit of the integrated water, wind, solar and storage base to form power supply partitions of different levels, each power supply partition having a corresponding level of transmission sections for accessing a higher adjacent level, thereby forming transmission sections of different levels; Step S2.2, according to the level of each transmission section and the hierarchical relationship of the associated transmission sections of each level, encode each transmission section to clarify the hierarchical relationship of each transmission section; Step S2.3, according to the superior-subordinate relationship of each transmission section, a directed graph is used to construct a transmission section topology relationship; Step S2.4, according to the topological relationship of the transmission sections, in order from low level to high level, a depth-first search method is used to determine the calculation order of each transmission section; Step S2.5, according to the access mode of each dispatching unit, directly associate it with the corresponding hierarchical transmission section as the dispatching unit directly associated with the hierarchical transmission section, thereby constructing the dispatching unit coupling relationship of the integrated water, wind, solar and storage base; Step S3, based on the coupling relationship of the dispatching units of the integrated water, wind, solar and storage base, implement the refined modeling and solution of the integrated water, wind, solar and storage dispatching, and obtain the output of each dispatching unit in each time period t during the dispatching period, including: Step S3.1, determining the output of each hydropower dispatching unit in each time period t during the dispatching period and the output of each energy storage dispatching unit in each time period t during the dispatching period; In step S3.1, the output of each hydropower dispatching unit in each time period t during the dispatching period and the output of each energy storage dispatching unit in each time period t during the dispatching period are determined, specifically including: Step S3.1.1, estimate the average output of each hydropower dispatching unit during the dispatching period: According to the coupling relationship of the dispatching units of the integrated water, wind, solar and storage base, the topological relationship of the hydropower dispatching units is determined, so that the hydropower dispatching units are sorted according to the upstream and downstream hydraulic connections; starting from the upstream, the average downstream flow of each hydropower dispatching unit during the dispatching period is calculated in sequence, and the average output of each hydropower dispatching unit during the dispatching period is estimated through the refined allocation method of water-determined current flow; Step S3.1.2, generating the initial output process of each hydropower dispatching unit in each time period t during the dispatching period: Taking the average output, expected output and forced output of each hydropower dispatching unit in the dispatching period as constraints, the successive load shedding method is adopted. On the basis of the base transmission capacity curve after deducting the designed output of the wind and solar power stations, the power balance calculation is performed for each hydropower dispatching unit separately to obtain the initial output process of each hydropower dispatching unit in each time period t during the dispatching period; Step S3.1.3, correct the initial output process of each hydropower dispatching unit in each time period t during the dispatching period, and obtain the corrected initial output process of each hydropower dispatching unit in each time period t during the dispatching period: The unit combination vibration zone analysis method is used to determine whether the initial output process of the hydropower dispatching unit in each time period t during the dispatching period meets the stable operation requirements. If there is an output in the time period t that does not meet the requirements, the output in the time period t that does not meet the requirements is corrected to the stable operation boundary, thereby obtaining the corrected initial output process of the hydropower dispatching unit in each time period t during the dispatching period; Step S3.1.4, taking the initial output process of the hydropower dispatching unit in each time period t during the dispatching period obtained in step S3.1.3 as the initial solution, taking the output of the hydropower dispatching unit in each time period t during the dispatching period as the decision variable, taking the minimum fluctuation of the output process of the hydro-wind-solar-storage multi-energy complementary power generation system as the dispatching target F, and establishing a dispatching model; Among them: The scheduling target F is: Where: T is the number of time periods in the scheduling period; P t is the output of the hydropower, wind power, photovoltaic and energy storage multi-energy complementary power generation system in time period t, is the sum of the outputs of the hydropower dispatching unit, wind power dispatching unit, photovoltaic dispatching unit and energy storage dispatching unit in time period t; The average output of the hydropower, wind, solar and storage multi-energy complementary power generation system during the dispatch period; Step S3.1.5, using the bee colony algorithm to solve and optimize the dispatching model established in step S3.1.

4. In the optimization process, the wind and solar power output adopts the wind and solar power design output as a constant value; the extraction state and output of the energy storage dispatching unit are determined according to the hydropower output and wind and solar power output of the hydro-wind-solar-storage multi-energy complementary power generation system, and controlled by the energy storage balance constraint; The heuristic constraint processing mechanism is used to modify the output of the hydropower dispatching unit in each time period t during the dispatching period and the extraction state and output of the energy storage dispatching unit in each time period t during the dispatching period in turn, so as to ensure the feasibility of the solution obtained by the algorithm optimization. Thus, the output of the hydropower dispatching unit in each time period t during the dispatching period and the extraction state and output of the energy storage dispatching unit in each time period t during the dispatching period are finally solved; Step S3.2, determining the output of each wind power dispatching unit in each time period t during the dispatching period and the output of each photovoltaic dispatching unit in each time period t during the dispatching period; Step S3.2, determining the output of each wind power dispatching unit in each time period t during the dispatching period and the output of each photovoltaic dispatching unit in each time period t during the dispatching period, specifically includes: Step S3.2.1, in the hydropower, wind, photovoltaic and storage multi-energy complementary power generation system, the output of each hydropower dispatching unit in each time period t during the dispatching period and the output of each energy storage dispatching unit in each time period t during the dispatching period have been solved, and the initial output of the wind power dispatching unit and the photovoltaic dispatching unit is given as the design output; Step S3.2.2, traversing each transmission section in sequence according to the calculation order of the transmission sections determined in step S2.4; Step S3.2.3: for the current traversed transmission section, it is represented as transmission section ts j , and its corresponding power partition is represented as power partition j; according to the output of the hydropower dispatching unit in each time period t during the dispatching period and the output of the energy storage dispatching unit in each time period t during the dispatching period, as well as the latest output of the wind power dispatching unit and the latest output of the photovoltaic dispatching unit, the power partition output recursive function is used to calculate the output P of power partition j in time period t Sj,t : in: is the recursive function of power partition j in time period t; They are respectively the transmission section ts with power supply partition j in time period t j The sum of the outputs of all directly related hydropower dispatching units, the sum of the outputs of all energy storage dispatching units, the sum of the outputs of all wind power dispatching units, and the sum of the outputs of all photovoltaic dispatching units; D Hj , D Pj , D Qj , D Rj , respectively, the transmission section ts with power supply partition j j The number set of directly related hydropower dispatching units, the number set of energy storage dispatching units, the number set of wind power dispatching units and the number set of photovoltaic dispatching units; m, k, q, r are the numbers of the hydropower dispatching unit, the energy storage dispatching unit, the wind power dispatching unit and the photovoltaic dispatching unit respectively; They are the output of the hydropower dispatching unit m in time period t, the output of the energy storage dispatching unit k in time period t, the latest output of the wind power dispatching unit q in time period t, and the latest output of the photovoltaic dispatching unit r in time period t; D is the sum of the outputs of all the next-level power partitions in power partition j after traversal and constraint processing in time period t; Sj Represents the set of sub-level power partition numbers, where l is the sub-level power partition number; Represents the recursive function of power partition l in time period t; represents the output upper limit of the next-level power partition l in time period t; Step S3.2.4, determine the output P of power partition j in time period t Sj,t , whether it meets the transmission section ts j The power transmission restriction constraint function is: Sj,t ≤N S,j , where N S,j is the transmission section ts of power supply zone j j If the power transmission restriction is satisfied, return to step S3.2.3 and continue to traverse the next power transmission section; if not satisfied, execute step S3.2.5; Step S3.2.5, combining the network topology of transmission sections at all levels and the power transmission restrictions of the transmission sections of each power partition, correct the output of all wind power dispatching units and all photovoltaic dispatching units in power partition j to obtain the latest output of each wind power dispatching unit in time period t and the latest output of each photovoltaic dispatching unit in time period t; then return to step S3.2.3 and continue to traverse the next transmission section.

2. According to claim 1, a refined scheduling modeling method based on the hierarchical cross-section topological relationship of the integrated water, wind, solar and storage base is characterized in that: In step S1, the division method of the hydropower dispatching unit is: The division is made based on the regulation and storage capacity of the hydropower station and the hydraulic and power connections between cascade hydropower stations, with a single hydropower station or a group of hydropower stations serving as a hydropower dispatching unit. Specifically, a hydropower station with a regulation capacity of daily or weekly or more forms a hydropower dispatching unit on its own; a run-of-river hydropower station and its directly adjacent upstream hydropower stations with the same power consumption direction form a hydropower dispatching unit together.

3. According to claim 1, a refined scheduling modeling method based on the hierarchical cross-section topological relationship of the integrated water, wind, solar and storage base is characterized in that: In step S1, the division method of the wind power dispatching unit and the photovoltaic dispatching unit is the same, and the wind power dispatching unit and the photovoltaic dispatching unit are both called new energy dispatching units, and the division method is: Taking into account the output characteristics, output correlation, grid access nodes, wind speed and time-varying characteristics of new energy power stations, the actual output process or simulated output process of each new energy power station is used as the feature input, and the Spearman rank correlation coefficient is used as the output correlation analysis and evaluation index of the new energy power station. The Canopy-K-means comprehensive clustering algorithm combining Canopy clustering and K-means clustering is used to aggregate a single new energy power station to a new energy power station cluster, and a new energy power station cluster is obtained. Each new energy power station cluster forms a new energy dispatching unit. The specific aggregation process is as follows: Step S1.1: Rough clustering of new energy power station clusters using the Canopy algorithm Step S1.1.1, establish a data set X = {X1, X2, ..., X n ,...,X N }; where N represents the number of new energy power stations participating in the cluster analysis; X n represents the output sequence of the new energy power station n, n=1,2,…,N; X n ={x1,x2,...,x t ,…,x T }, where x1, x2, ..., x t ,…,x T represents the output of the new energy power station n in each time period, and T represents the number of time periods; Step S1.1.2, using a cross-validation method to determine the first distance threshold T1 and the second distance threshold T2 of the output correlation of the new energy power station, and T1>T2; Step S1.1.3, randomly select the output sequence X of the new energy power station a from the data set X of the new energy power station output cluster a , a=1,2,…,N, with the output sequence X of new energy power station a a As the center of the canopy class, calculate the output sequence X of the new energy power station a respectively a The output sequence X of each other new energy power station b b Output correlation distance b=1,2,…,N, the output correlation distance with the new energy power station a All new energy power stations with a distance less than the second distance threshold T2 are added to the canopy class where the center is located, marked with a strong mark and deleted from the data set X; the output correlation distance with the new energy power station a is All new energy power stations between the first distance threshold T1 and the second distance threshold T2 are added to the canopy class where the center is located and marked with a weak mark. The new energy power stations marked with a weak mark will participate in the next round of clustering process and become the center or member of the new canopy class; the output correlation distance with the new energy power station a will be All new energy power stations with a distance greater than the first distance threshold T1 are individually used as the center of a new canopy class for cluster analysis; this process is repeated until all new energy power stations are assigned to one or more canopy classes, thereby achieving coarse clustering of new energy power station clusters and dividing them into K coarse clusters; in: The following formula is used to obtain: in: and They are the output sequence X of new energy power station a a and the output sequence X of the new energy power station b b The rank after transformation; and They are the output sequence X of new energy power station a a and the output sequence X of the new energy power station b b The standard deviation of is the output sequence X of the new energy power station a a and the output sequence X of the new energy power station b b The Spearman rank correlation coefficient ranges from [-1,1]. The larger the absolute value, the better the correlation. represent and The covariance of Step S1.2, based on the K coarse clusters obtained by the Canopy algorithm, considering the output correlation distance between different renewable energy power stations, each coarse cluster is further clustered and analyzed by the K-means algorithm to obtain multiple subsets of renewable energy power stations, each of which forms a renewable energy dispatching unit.

4. The refined scheduling modeling method based on the topological relationship of the hierarchical sections of the integrated water, wind, solar and storage base according to claim 1 is characterized in that: In step S1, the energy storage dispatching unit is divided as follows: Combined with the operating characteristics, storage capacity, and cycle efficiency of each energy storage power station in the integrated base, a single energy storage power station or a group of energy storage power stations is used as a dispatching unit; specifically, for new energy storage power stations, if the energy storage equipment form, operating characteristics, cycle efficiency, and power access point are the same, an energy storage dispatching unit is formed; if the combination conditions are not met, a single new energy storage power station is used as an energy storage dispatching unit; For pumped-storage power stations, due to differences in development forms, upper and lower reservoir capacity conditions, head conditions, unit output characteristics, and power station operation modes, a single pumped-storage power station is used as the energy storage dispatching unit.

5. The refined scheduling modeling method based on the topological relationship of the hierarchical sections of the integrated water, wind, solar and storage base according to claim 1 is characterized in that: In step S2.1, the specific method of power source partitioning and transmission section classification is as follows: According to the spatial distribution of each dispatching unit, the power collection method and the voltage level of the power transmission channel, power partitioning is carried out to form recursive power partitioning and hierarchical transmission sections; Specifically, the dispatching units with similar geographical locations, the same power input points, and directly associated with the same power transmission channel, as well as the power supply zones of the next level directly associated with the power transmission channel, are merged into the power supply zone of this level, and the power transmission channel of the power supply zone of this level is defined as the transmission section of the corresponding level; Specifically, the transmission section with the highest voltage level of the integrated water, wind, solar and storage base is defined as the first-level transmission section R1; the number of the first-level transmission sections R1 is one or more; for each first-level transmission section R1, a first-level power partition C1 is formed; each of the first-level power partitions C1 includes the first-level transmission section R1 and the first-level dispatching unit U1 directly associated with the power collection end of the first-level transmission section R1, and the second-level power partition C2; Each secondary power supply partition C2 includes a secondary transmission section R2 directly connected to the power collection end of the primary transmission section R1, a secondary dispatching unit U2 directly associated with the power collection end of the secondary transmission section R2, and a tertiary power supply partition C3; By analogy, a multi-level cascaded power supply partition and a multi-level transmission section are formed, forming a coupling structure of power supply partition and transmission section classification; wherein each transmission section can be directly associated with one or more transmission sections of the next adjacent level, but has only one directly associated transmission section of the previous adjacent level; In step S2.2, each transmission section is coded as follows: The first-level transmission section is coded as: 1A; where 1 represents the level of the transmission section; A represents the sequence number of the first-level transmission section among all the first-level transmission sections; The secondary transmission section is coded as: 1A-2B, which means: the Bth secondary transmission section connected to the power collection end of the primary transmission section 1A; The tertiary transmission section is coded as: 1A-2B-3C, which means: the Cth tertiary transmission section connected to the power collection end of the secondary transmission section 1A-2B.

6. The refined scheduling modeling method based on the topological relationship of the hierarchical sections of the integrated water, wind, solar and storage base according to claim 1 is characterized in that: Step S2.3 is specifically as follows: The vertices of the directed graph represent the transmission sections of each level, the edges of the directed graph represent the direct connection between the transmission sections of adjacent levels, and the direction of the edge of the directed graph represents the superior-subordinate relationship between two transmission sections with direct connection, that is, from the lower transmission section to the upper transmission section, so as to construct a directed graph by combining the superior-subordinate relationship of each transmission section and form the topological relationship of the transmission section; Step S2.4 is specifically as follows: From the topological relationship of the transmission sections, locate the highest level transmission section, use it as the root s of the directed graph, and put it in the target sorting queue line; Then, in the transmission section topology relationship, determine the transmission section set U of the direct subordinate of the root s, traverse the transmission sections u in the transmission section set U in turn, and add each traversed transmission section u to the target sorting queue line in order; For each transmission section u in the transmission section set U, locate the transmission section set directly below it, and add each traversed transmission section to the target sorting queue line in order; this cycle continues until all transmission sections are searched; The reverse order of the transmission sections in the target sorting queue line is the calculation order of each transmission section determined.

7. The refined scheduling modeling method based on the topological relationship of the hierarchical sections of the integrated water, wind, solar and storage base according to claim 1 is characterized in that: In step S3.1.5, the heuristic constraint processing mechanism is used to sequentially correct the output of the hydropower dispatching unit and the output of the energy storage dispatching unit in each time period t of the dispatching period, including: Step S3.1.5.1, for each hydropower dispatching unit i obtained by the bee colony algorithm in this optimization, the output of each time period t during the dispatching period is expressed as The power balance constraint strategy is used for correction to obtain the corrected output of hydropower dispatching unit i in each time period t during the dispatching period. The specific method is: ①Through the formula Get the full-time power generation of hydropower dispatching unit i during the dispatching period Where Δt is the length of each time period t; ②If Among them, E i is the theoretical power generation of hydropower dispatching unit i during the dispatching period; ε is the deviation margin, then the power balance constraint is corrected, otherwise, no correction is made; The correction method is: Using formula The power difference Divide into T equal parts, T is the number of scheduling periods, and deltE is the power value of each equal part; like That is, the total power generation of hydropower dispatching unit i Greater than the theoretical power generation E i , then set t=1,2,...,T,P i h,min is the forced output of hydropower dispatching unit i; if not, the full-time power generation of hydropower dispatching unit i Less than the theoretical power generation E i , then set P i h,max is the expected output limit of hydropower dispatching unit i; the correction is repeated in this cycle until the full-time power generation of hydropower dispatching unit i is And the theoretical power generation E i The difference is controlled within the deviation margin ε; Step S3.1.5.2: for the corrected output of hydropower dispatching unit i in each time period t during the dispatching period Adopt the hydropower final water level control constraint repair strategy to Further correction is made to obtain the corrected output of hydropower dispatching unit i in each time period t during the dispatching period: Specifically, the final water level of the hydropower dispatching unit i in each time period t during the dispatching period and the final water level of the dispatching period are obtained by using the refined allocation method of power-determined water output; if the final water level of the hydropower dispatching unit i in time period t does not meet the requirements, the output of the time period is re-determined by reverse calculation through flow control or water level control method; The specific method is: ① If the final water level of hydropower dispatching unit i in each time period t during the dispatching period is greater than the given final water level, it indicates that the power generation is insufficient and the output needs to be increased according to a certain power step ΔE0 Get the adjusted output of hydropower dispatching unit i in each time period t during the dispatching period The method is: The power step length ΔE0 is divided into T parts according to the number of time periods, and the power length of each part is Adjust each period from period 1. If restricted, a single period cannot be If the power is completely consumed, the remaining power will be allocated to its adjacent time periods; ② According to the adjusted output of hydropower dispatching unit i in each time period t during the dispatching period Recalculate the final water level of hydropower dispatching unit i in each time period t during the dispatching period by using electricity to determine water, return to step ①, and repeat this cycle until the final water level requirement is met. Finally, the corrected output of hydropower dispatching unit i in each time period t during the dispatching period is obtained, which is expressed as ③ If the final water level of hydropower dispatching unit i in each time period t during the dispatching period is lower than the given final water level, the output will be reduced in each time period for correction; Step S3.1.5.3, adopt the heuristic constraint processing mechanism to correct the output of the energy storage dispatching unit in each time period t during the dispatching period: ① Determine the working status of the energy storage dispatching unit: For the kth energy storage dispatching unit put into operation, its working state is jointly determined by the output process of the water, wind, and light dispatching units in the water-wind-solar-storage multi-energy complementary power generation system and the first k-1 energy storage dispatching units put into operation: if the current water-wind-solar-storage multi-energy complementary power generation system has power abandonment in time period t, and the energy storage dispatching unit k still has energy storage space, then the energy storage dispatching unit k is in the energy storage state in time period t; if there is idle capacity in the sending section of the water-wind-solar-storage multi-energy complementary power generation system, and the energy storage dispatching unit k still has the discharge capability, then the energy storage dispatching unit k is in the power generation state in time period t; ② Generate energy storage constraint corridor: The storage energy E of energy storage scheduling unit k at the beginning of the scheduling period k,0 Based on the forward sequence calculation method, the energy storage constraint corridor is generated, and the formula is: in: and are the upper limit of energy storage output and power generation output of energy storage dispatching unit k in time period t respectively; E k,t,max1 and E k,t,min1 , are respectively the upper and lower limits of the energy storage capacity of the energy storage dispatching unit k at the end of time period t during the positive sequence calculation; E k,t-1,max1 and E k,t-1,min1 , are the upper and lower limits of the energy storage capacity of the energy storage dispatching unit k at the end of time period t-1 during the positive sequence calculation; E k,max and E k,min , are the upper and lower limits of the designed storage capacity of the energy storage dispatching unit k respectively; is the cycle efficiency of energy storage dispatch unit k; E k,0,max1 and E k,0,min1 , are respectively the upper and lower limits of the energy storage capacity of the energy storage dispatching unit k at the beginning of the dispatching period when the forward sequence calculation method is adopted; The storage energy E of energy storage dispatch unit k at the end of the dispatch period k,T Based on the reverse calculation method, the energy storage constraint corridor is generated, and the formula is: Where: E k,t,max2 and E k,t,min2 , are the upper and lower limits of the energy storage capacity of the energy storage dispatching unit k at the end of time period t when calculating in reverse order; E k,t+1,max2 and E k,t+1,min2 , are the upper and lower limits of the energy storage capacity of the energy storage dispatching unit k at the end of time period t+1 when calculating in reverse order; E k,T,max2 and E k,T,min2 , are respectively the upper and lower limits of the energy storage capacity of the energy storage dispatch unit k at the end of the dispatch period when the reverse calculation method is adopted; Further take the envelope to generate the energy storage constraint corridor (E k,t,min ,E k,t,max ): E k,t,max and E k,t,min , are the upper and lower limits of the energy storage capacity of the energy storage scheduling unit k in time period t, respectively; ③ According to the energy storage capacity constraint corridor (E k,t,min ,E k,t,max ), the extraction state and output of the energy storage dispatching unit k in time period t obtained by the bee colony algorithm in step S3.1.5 are corrected.

8. The refined scheduling modeling method based on the topological relationship of the hierarchical sections of the integrated water, wind, solar and storage base according to claim 1 is characterized in that: Step S3.2.5 is specifically: Step S3.2.5.1, based on the output P of power partition j in time period t Sj,t and transmission section ts j Power supply limit N S,j The difference between the transmission section ts j The curtailment of power in period t Step S3.2.5.2: If power partition j does not contain the next-level power partition, the abandoned power Allocate according to the proportion of online electricity of the wind and solar dispatching units directly associated with power partition j, thereby correcting the output of the corresponding wind and solar dispatching units in time period t, and updating the online electricity of the wind and solar dispatching units in time period t; wherein the wind and solar dispatching units include wind power dispatching units and photovoltaic dispatching units; If power partition j contains the next-level power partition, the power will be discarded. The online power consumption of all wind and solar dispatching units directly associated with power partition j and indirectly associated with the next-level power partition is allocated to correct the output of the corresponding wind and solar dispatching unit in time period t and update the online power consumption of the wind and solar dispatching unit in time period t.

Images