A method for selecting hydraulic units based on vibration zone limitation and technical and economic indicators

Abstract

The present invention discloses a method for selecting a hydraulic unit based on vibration zone limitation and technical and economic indicators, which establishes a double-layer dynamic programming unit load distribution model; establishes technical and economic evaluation indicators; and obtains an optimal unit selection scheme based on the double-layer dynamic programming unit load distribution model and the technical and economic evaluation indicators. From the perspective of the operation of the hydraulic unit, during the process of load distribution of each unit selection in the plant, by considering the vibration zone limitation and multiple technical and economic indicators, the operation safety and economy are taken into account, the long-term operation of the unit in the vibration zone is avoided, and at the same time, the load distribution results of the unit are quantified from multiple angles such as additional start-stop times and output fluctuations, providing a significant comparison between indicators and providing favorable support for the optimal unit selection scheme.

Description

Technical Field

[0001] The present invention belongs to the field of power generation operation of hydraulic units, and particularly relates to a method for selecting hydraulic units based on vibration zone limitation and technical and economic indexes. Background Art

[0002] As an important part of automatic generation control, the load distribution of hydropower station units plays an important role in the dispatching automation of power systems, especially in frequency modulation auxiliary services. The goal of unit load distribution is to automatically distribute the total load target value of the whole plant among units under the premise of ensuring the safe and reliable operation of the units, taking into account various factors such as water conditions, unit capacity, unit non-operable area (cavitation vibration area), and operating conditions, so as to achieve the purpose of economic operation.

[0003] When a hydropower station undertakes the power grid frequency modulation task, the load task it receives is at the second level, and the load shows characteristics such as small change amplitude and large change rate in a short period of time. At this time, in addition to pursuing the economy of the operation mode, the safety and real-time performance of the hydropower station need to be emphasized more. Especially for some cascade power stations with annual regulation, in the face of multi-fluctuating loads, there is a high probability that the units will encounter large-scale load transfer problems among units in a short time, and the units frequently cross the vibration zone, which is not recommended in engineering practice.

[0004] Therefore, how to carry out a rehearsal for the unit load distribution problem at the stage of unit selection, based on multiple objectives such as the water consumption of the hydropower station and the source-load matching situation, by comparing the load distribution results of the same operating day under each selection scheme, statistically analyzing the situation of the unit in the stable operating power range, and comparing the unit in the vibration zone, crossing the vibration zone or related economic evaluation indexes, is an urgent problem to be solved in the current hydropower station unit load distribution work. Summary of the Invention

[0005] In view of the above-mentioned defects existing in the prior art, the present invention provides a method for selecting hydraulic units based on vibration zone limitation and technical and economic indexes, taking into account both operation safety and economy, in order to consider vibration zone limitation and related technical and economic evaluation indexes in the load distribution of hydraulic units.

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

[0007] The present invention discloses a method for selecting hydraulic units based on vibration zone limitation and technical and economic indexes, the method comprising:

[0008] Establishing a two-layer dynamic programming unit load distribution model;

[0009] Establishing technical and economic evaluation indexes;

[0010] Obtain the optimal unit selection scheme based on the above-mentioned two-layer dynamic programming unit load distribution model and technical and economic evaluation indicators.

[0011] Furthermore, the two-layer dynamic programming unit load distribution model takes the minimum consumption flow as the objective function.

[0012] Furthermore, the definition formula of the objective function is:

[0013]

[0014] In the formula, T is the number of time periods; L is the number of units put into operation by the power station at time t; l is the unit number of the power station; P l,t is the load borne by the l-th unit at time t; H t is the working head at time t; Q l,t is the working flow required for the l-th unit to bear the power plant load at time t.

[0015] Furthermore, the constraint conditions of the two-layer dynamic programming unit load distribution model include: power balance constraint, unit output constraint, unit flow constraint, unit start-stop constraint, and unit vibration zone constraint.

[0016] Furthermore, the definition formula of the constraint conditions is:

[0017] Power balance constraint:

[0018]

[0019] Unit output constraint:

[0020] N l,min ≤P l,t ≤N l,max

[0021] Unit flow constraint:

[0022] Q l,min ≤Q l,t ≤Q l,max

[0023] Unit start-stop constraint:

[0024]

[0025] Unit vibration zone constraint:

[0026] P l,t ≤P l,min |P l,t ≥P l,max

[0027] In the formula, N l,min , N l,maxare the minimum and maximum possible outputs of the l-th unit respectively, Q l,min , Q l,max are the minimum and maximum possible reference flows of the l-th unit respectively, N t is the total load borne by the power station at time t, P l,min is the lower limit of the vibration zone of the l-th unit, P l,max is the upper limit of the vibration zone of the l-th unit, P l,t is the output of the l-th unit at the t-th time period, is the allowable continuous start-stop time of the l-th unit.

[0028] Furthermore, in the double-layer dynamic programming unit load distribution model, the dynamic programming algorithm is first used to calculate the optimal load distribution values and corresponding flow consumptions corresponding to all unit selection combinations within a single time period, and all time periods are repeatedly calculated to obtain the optimal solutions for unit selection combination load distribution in different time periods; the dynamic programming algorithm is then used to calculate the unit start-stop optimization between different time periods, and all units are repeatedly calculated to obtain the optimal unit selection combination between time periods.

[0029] Furthermore, the specific process of using the dynamic programming algorithm to calculate the optimal load distribution values and corresponding flow consumptions corresponding to all unit selection combinations within a single time period is as follows:

[0030] Set the variables of the dynamic programming algorithm. Among them, the number of started units in the current unit selection combination is used as the stage variable, the total input load is used as the state variable, and the load borne by a single unit in the current unit selection combination is used as the decision variable;

[0031] Obtain the corresponding power generation flow consumption according to the decision variable;

[0032] Set the objective function. Among them, the objective function is set to find the optimal load distribution method in the load discrete space to minimize the total power generation flow consumption in the current time period, and the value range excludes the vibration zone of the unit.

[0033] Furthermore, the specific process of using the dynamic programming algorithm to calculate the unit start-stop optimization between different time periods is as follows:

[0034] Set the variables of the dynamic programming algorithm. Among them, different time periods are used as the stage variable, the total number of unit selection combinations is used as the state space, and the current unit selection combination is used as the decision variable;

[0035] Set the objective function. Among them, the objective function is set to minimize the total power generation flow consumption in the current time period.

[0036] Furthermore, the technical and economic evaluation indicators include: total water consumption, number and proportion of non-recommended operation zone time periods, additional start-stop times, output fluctuation ratio, and average efficiency.

[0037] The hydraulic unit selection method based on vibration zone limitation and technical and economic indicators proposed by the present invention, compared with the existing hydraulic unit selection methods that usually only consider economy, starting from the operation perspective of hydraulic units, during the process of load distribution in the plant for each unit selection, by considering vibration zone limitation and multiple technical and economic indicators, taking into account both operation safety and economy, avoiding long-term operation of the unit in the vibration zone, and at the same time quantifying the unit load distribution results from multiple perspectives such as additional start-stop times and output fluctuations, providing significant comparisons between indicators, and providing favorable support for the optimal unit selection scheme. Brief Description of the Drawings

[0038] The drawings forming a part of this application are used to provide a further understanding of the present invention. The schematic embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention. In the drawings:

[0039] Figure 1 It is a schematic diagram of the working process of Embodiment 1 of the present invention;

[0040] Figure 2 It is a schematic diagram of the load distribution result of Unit Selection Scheme 1 shown in Table 3 on a certain distribution day during the flood season shown in Table 1 of Embodiment 1 of the present invention;

[0041] Figure 3 It is a schematic diagram of the load distribution result of Unit Selection Scheme 2 shown in Table 3 on a certain distribution day during the flood season shown in Table 1 of Embodiment 1 of the present invention;

[0042] Figure 4 It is a schematic diagram of the load distribution result of Unit Selection Scheme 3 shown in Table 3 on a certain distribution day during the flood season shown in Table 1 of Embodiment 1 of the present invention;

[0043] Figure 5 It is a schematic diagram of the load distribution result of Unit Selection Scheme 4 shown in Table 3 on a certain distribution day during the flood season shown in Table 1 of Embodiment 1 of the present invention;

[0044] Figure 6 It is a schematic diagram of the load distribution result of Unit Selection Scheme 1 shown in Table 3 on a certain distribution day during the dry season shown in Table 2 of Embodiment 1 of the present invention;

[0045] Figure 7 It is a schematic diagram of the load distribution result of Unit Selection Scheme 2 shown in Table 3 on a certain distribution day during the dry season shown in Table 2 of Embodiment 1 of the present invention;

[0046] Figure 8 It is a schematic diagram of the load distribution result of Unit Selection Scheme 3 shown in Table 3 on a certain distribution day during the dry season shown in Table 2 of Embodiment 1 of the present invention;

[0047] Figure 9It is a schematic diagram of the load distribution result of the unit selection scheme 4 shown in Table 3 on a certain allocation day during the dry period in Example 1 of the present invention as shown in Table 2. Detailed implementation mode

[0048] To enable those skilled in the art to better understand the technical solution of the present invention, the present invention will be further described in detail below in conjunction with the accompanying drawings and specific implementation modes.

[0049] Example 1

[0050] Refer to Figure 1 , this embodiment provides a hydraulic unit selection method based on vibration zone limitation and technical and economic indicators. The method includes the following steps:

[0051] S1: Establish a two-layer dynamic programming load distribution model considering the vibration zone;

[0052] Specifically, the model takes the minimum consumption flow as the objective function, and the established objective function is as follows:

[0053]

[0054] In the formula, T is the number of time periods; L is the number of units put into operation in the power station at time t; l is the unit number of the power station; P l,t is the load borne by the l-th unit at time t, with the unit of MW; H t is the working head at time t, with the unit of m; Q l,t is the working flow required for the l-th unit to bear the power plant load N t at time t.

[0055] According to the actual situation, the constraint conditions of the model are simplified. Only the upper and lower limits of the power station output are considered; the comprehensive utilization requirements above the reservoir dam (power generation, irrigation, inundation, and flood control) are not considered, and only the highest and lowest water level constraints of the reservoir are considered; the comprehensive utilization requirements below the reservoir dam (navigation below the dam, river channel scouring) are not considered, and only the minimum discharge flow constraint, ecological flow, and maximum discharge flow amplitude constraint are considered; the established model constraint conditions are as follows:

[0056] Power (load) balance constraint:

[0057]

[0058] Unit output constraint:

[0059] N l,min ≤P l,t ≤N l,max

[0060] Unit flow constraint:

[0061] Q l,min≤Q l,t ≤Q l,max

[0062] Unit start-stop constraint:

[0063]

[0064] Unit vibration zone constraint:

[0065] P l,t ≤P l,min |P l,t ≥P l,max

[0066] In the formula, N l,min , N l,max are the minimum and maximum possible outputs of the l-th unit respectively. They are restricted by the single-unit capacity and are also functions of the working head H t and the working flow rate Q l,t ; Q l,min , Q l,max are the minimum and maximum possible diverted flow rates of the l-th unit respectively. They are functions of the working head H t ; N t is the total load borne by the power station in the t-th time period; P l,min is the lower limit of the vibration zone of the l-th unit, P l,max is the upper limit of the vibration zone of the l-th unit, P l,t is the output of the l-th unit in the t-th time period, t is the minimum scheduling time period, is the allowable continuous start-stop time of the l-th unit.

[0067] The algorithm used in the two-layer dynamic programming load distribution model is as follows:

[0068] For the first time, use the dynamic programming algorithm to calculate the optimal load distribution value and the corresponding flow consumption corresponding to all unit selection combinations within a single time period, repeat the calculation for all time periods, and obtain the optimal load distribution solutions for different unit selection combinations within different time periods; for the second time, use the dynamic programming algorithm to calculate the unit start-stop optimization between different time periods, repeat the calculation for all units, and obtain the optimal unit selection combination between time periods.

[0069] (1) Optimal load distribution algorithm for units in a single time period

[0070] Given the power station time period load and the upstream incoming water flow rate, solve the unit power generation load and power generation flow rate based on the dynamic programming method. Use the number of units started in the current unit selection combination as the stage variable, the total input load as the state variable, and the load borne by a single unit in the current unit selection combination as the decision variable. For each unit start-stop combination, use the dynamic programming method to determine the corresponding optimal unit load distribution plan with the goal of minimizing the flow consumption.

[0071] The specific algorithm steps are as follows:

[0072] a. Given the initial water level, the in-plant load Pt (t = 1, …, T) is allocated period by period starting from the first period; there are L units in the plant, and each unit is represented by l (l = 1, …, L); since the unit numbers are fixed, the combination of their start-stop states is 2^L, and each unit combination is represented by a positive integer J (J = 1, …, 2^L). This positive integer can be conveniently expressed in binary form to represent the off and on states of the units (0 - off, 1 - on). For example, the binary expression of the integer 15 is 1111, representing that all 4 units are in the on state. In each period, this 2^L combination of unit states is looped through, and the load decomposition dynamic programming algorithm in step b is executed to decompose the total period load Pt to the units;

[0073] b. Within a single period, for the Jth unit combination, determine the number of operating units m (0 ≤ m ≤ L). Using the dynamic programming algorithm, with the number of available generating units m as the stage variable, the total available load Pt as the state variable, and the load borne by a single unit P l,t as the decision variable, determine the decision variable P l,t and then solve for the corresponding power generation flow Q l,t , perform a water balance calculation to obtain the water level and reservoir capacity at the end of the period, which are used as the water level and reservoir capacity inputs for the next period. The objective function is to find the optimal load distribution plan within the load discrete space to minimize the total power generation consumption flow Qt of the period. During the search within the load discrete space at each stage of the dynamic programming, the vibration zones of the units {P l,min ≤P l,t ≤P l,max} are excluded from the range of values. Record the optimal allocated load P l,t corresponding to each combination and the total power generation consumption flow Qt of the period;

[0074] c. At this point, the load distribution for one period is completed. Return to step a to start allocating the next period, and calculate until the end of the last period to obtain a set of 2^L sets of in-plant load distribution solutions P l,t and Q l,t .

[0075] (2) Unit start-stop optimization algorithm between different periods

[0076] Based on the optimization results of the first stage, the dynamic programming algorithm is used again. Taking different optimization periods t as the stage variable, the total number of unit combinations 2^L as the state space, the decision variable as the unit combination J, and the objective function as minimizing the total power generation consumption flow Q.

[0077] The specific algorithm steps are as follows:

[0078] a. The decision variable J for each optimization period t tThe decision variable J at time period t-1 t-1 are all converted into binary numbers, compared bit by bit, and the start-stop flow AcrossQ of the l-th unit at time period t is calculated according to the following rules l,t :

[0079]

[0080] where Q s and Q e are the flow consumption for starting and shutting down the unit respectively, and are generally fixed values.

[0081] Meanwhile, if the unit is in the on state at both time period t and t-1, check whether the unit load crosses the vibration zone from time period t-1 to t. Calculate the penalty flow AcrossQ for crossing the vibration zone according to the following rules l,t :

[0082]

[0083] where Q c is the flow consumption generated by the unit crossing the vibration zone, and is generally a fixed value.

[0084] Similarly, corresponding penalty flow consumption PunishQ can be designed for constraints such as unit operation time and start-stop times l,t .

[0085] b. After determining the time period decision variables, calculate the total flow consumption after adding costs such as start-stop

[0086]

[0087] c. Return to step a to start calculating the total flow consumption after adding costs such as start-stop for different time periods of the next unit combination. Traverse a total of T*2^L combinations, and use dynamic programming to find the optimal unit combination for each time period with the goal of minimizing the total flow consumption for each time period.

[0088] This double dynamic programming algorithm first forms a set of solutions mainly based on single-time period load distribution to obtain the optimal load distribution solutions of units within different time periods; then, on this basis, perform combinations between time periods, and at the same time use penalty flow to meet the constraint conditions to obtain the optimal unit combination between time periods.

[0089] During the solution process of dynamic programming, pay attention to accelerating the state space search process through conditional judgment. For example, in the calculation of the optimal load distribution of 2^L unit combinations for each time period, first judge the minimum number of units to be started. If the number of units started in the combination is lower than the minimum number of units to be started, then directly jump out without performing subsequent calculations of the unit power generation load and power generation flow, and continue to calculate the next combination.

[0090] S2: Establish technical and economic evaluation indicators

[0091] Specifically, when the unit load adjustment is small, the load distribution strategy should reduce the total adjustment frequency of the units, thereby reducing equipment wear and system energy consumption, and minimizing the number of units participating in the adjustment as much as possible. This strategy can also reduce the cumulative load deviation caused by multiple units participating in small load fluctuations, so as to better track the active power set value. Therefore, in order to compare the impacts of various unit selection schemes on the load distribution of the hydropower station, the present invention proposes five evaluation indicators from technical and economic perspectives, namely total water consumption, non-recommended operation area and ratio, additional start-stop times, output fluctuation, and unit evaluation efficiency, to evaluate the impacts of the scheme on the safe and stable operation of the units and the water energy utilization situation.

[0092] The defined technical and economic evaluation indicators are as follows:

[0093] A. Total water consumption

[0094]

[0095] In the formula, L is the number of units; T is the number of time periods; Q l,t (H t , P l,t ) represents the water consumption generated by the l-th unit in the t-th time period according to the water head and output in the t-th time period; Δt is the minimum scheduling scale; λ l,t represents the start-stop state of the l-th unit in the t-th time period, represented in binary, 1 for starting up and 0 for shutting down; n cross represents the number of times the unit crosses the vibration area; Q on_off , Q cross represent the additional start-stop water flow (excluding the initial start-up and the last shut-down) and the water flow for crossing the vibration area of the unit during the scheduling process. The calculation method of start-stop loss refers to the calculation method in "Application of Dynamic Programming Method in the Economic Operation of Hydropower Station".

[0096] B. Number of time periods and ratio of non-recommended operation area

[0097] This method sets the output range of 50% - 70% of the unit as the non-recommended operation area. Specifically, it is shown in the following formula:

[0098]

[0099] In the formula, T l,t is the running time when the unit operates within the non-recommended operation area range, with 15 minutes as a time period and the minimum scheduling scale. The ratio of the non-recommended operation area is obtained by dividing the non-recommended operation area by the start-up time periods of all units.

[0100] C. Additional start-stop times

[0101] For a single unit, only the initial startup action and the final shutdown action of this unit are considered throughout the scheduling process. If there are additional startup and shutdown operations in between, they are the additional startup and shutdown counts of the unit. When there is one additional startup and shutdown operation, the additional startup and shutdown count of the unit increases by two.

[0102] D. Output fluctuation ratio

[0103]

[0104] In the formula, P l,t represents the output of the l-th unit in the t-th time period, and represents the average output of the l-th unit throughout the scheduling process.

[0105] E. Average efficiency

[0106] For each stage of load distribution, after the working head H, working flow Q, and output N of the unit are determined, the efficiency of the unit's current operating condition can be approximately solved as shown in the following formula:

[0107]

[0108] When there are multiple units, the average operating efficiency can be expressed as:

[0109]

[0110] In the formula, represents the average efficiency of the unit operation in the t-th time period, and η l,t represents the operating efficiency of the l-th unit in the t-th time period.

[0111] S3: Obtain the data of the head and demand load changes of a certain typical year of the power station;

[0112] Specifically, select the data of the demand load changes of the power station at 96 time periods on a certain distribution day during the flood season of a normal water year and display it in Table 1. The average head of that day is 272.11 m;

[0113] Table 1

[0114]

[0115]

[0116] Select the data of the demand load changes of the power station at 96 time periods on a certain distribution day during the dry season of the same typical year and display it in Table 2. The average head of that day is 266.51 m;

[0117] Table 2

[0118]

[0119] S4: Select the optimal unit selection plan according to the double-layer dynamic programming load distribution model and technical and economic evaluation indicators;

[0120] Specifically, this embodiment shows the load distribution calculation and technical and economic index statistics of the selection plan shown in Table 3 under the two distribution days shown in step S3.

[0121] Table 3

[0122]

[0123] The load distribution calculation results of Plans 1 to 4 on the load distribution day shown in Table 1 are as Figures 2 to 5 shown, and the comparison of technical and economic indicators is shown in Table 4:

[0124] Table 4

[0125]

[0126] The load distribution calculation results of Plans 1 to 4 on the load distribution day shown in Table 2 are as Figures 6 to 9 shown, and the comparison of technical and economic indicators is shown in Table 5:

[0127] Table 5

[0128]

[0129] From Figures 2 to 9 and the information given in Tables 3 to 4, it can be seen that under the current flood season distribution day, when the number of the same units is the same, different runner models do not bring obvious differences. When there are four units, the selection of runner model a has lower water consumption compared with b; while when there are six units, the runner of model b has lower energy consumption and higher efficiency. Comparing the plans with different numbers of units under the same runner longitudinally, it can be known that the six-unit plan has a lower single-unit capacity, and can better coordinate the best matching between the unit load and the number of units when encountering short-term fluctuations in the demand load, which can bring a lower unit output fluctuation ratio. However, this also leads to more additional starts and stops of units and more operation hours in the non-recommended operation area during the flood season distribution day, increasing the loss during the power generation process. During the dry season distribution day, the number of non-recommended operation hours in the four-unit plan is more than that in the six-unit plan. After analysis, it is because the demand load of the power station during the dry season is lower, and most units are easy to work in the power generation range of 50% to 70%. And the six units have a lower single-unit power generation capacity, so they can have a better power generation range during the dry season, but this also brings higher power generation water consumption. Through the above analysis, it is concluded that among the four unit selection plans, the two plans with four units have better performance, and there are no significant differences between the different runner models of the two plans.

[0130] In addition, an embodiment of the present invention further provides a computer-readable storage medium, wherein the computer-readable storage medium may store a program, and when the program is executed, it includes some or all of the steps of any of the hydraulic unit selection methods based on vibration zone limitation and technical and economic indicators described in the above method embodiments.

[0131] In addition, in each embodiment of the present invention, each functional unit may be integrated in a processing unit, or each unit may exist physically alone, or two or more units may be integrated in one unit. The above integrated unit may be implemented in the form of hardware or in the form of a software functional unit.

[0132] If the above integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it may be stored in a computer-readable memory. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or all or part of the technical solution, may be embodied in the form of a software product. The computer software product is stored in a memory and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in each embodiment of the present invention. The aforementioned memory includes: various media such as USB flash drives, read-only memories (ROMs), random access memories (RAMs), mobile hard disks, magnetic disks, or optical discs that can store program codes.

[0133] The specific embodiments described above further elaborate on the purpose, technical solution, and beneficial effects of the present application. It should be understood that the above description is only the specific embodiments of the present application and is not used to limit the protection scope of the present application. Any modifications, equivalent replacements, improvements, etc. made on the basis of the technical solution of the present application shall be included in the protection scope of the present application.

Claims

1. A method for selecting a hydraulic unit based on vibration zone limitation and technical and economic indicators, characterized in that, It includes the following steps: Establish a double - layer dynamic programming unit load distribution model; Establish technical and economic evaluation indicators; Obtain the optimal unit selection scheme based on the double - layer dynamic programming unit load distribution model and technical and economic evaluation indicators; The double - layer dynamic programming unit load distribution model takes the minimum consumption flow as the objective function; The definition formula of the objective function is: Wherein, T is the number of time periods; L is the number of units put into operation at the power station in the t-th time period; l is the unit number of the power station; P l,t is the load borne by the l-th unit in the t-th time period; H t is the working head in the t-th time period; Q l,t is the working flow required for the l-th unit to bear the power plant load in the t-th time period; The constraint conditions of the double - layer dynamic programming unit load distribution model include: power balance constraint, unit output constraint, unit flow constraint, unit start - stop constraint, unit vibration zone constraint; The definition formula of the constraint conditions is: Power balance constraint: Unit output constraint: N l,min ≤P l,t ≤N l,max Unit flow constraint: Q l,min ≤Q l,t ≤Q l,max Unit start - stop constraint: Unit vibration zone constraint: P l,t ≤P l,min |P l,t ≥P l,max Wherein, N l,min , N l,max are respectively the minimum and maximum possible outputs of the l-th unit, Q l,min , Q l,max are respectively the minimum and maximum possible reference flows of the l-th unit, N t is the total load borne by the power station in the t period, P l,min is the lower limit of the vibration zone of the l-th unit, P l,max is the upper limit of the vibration zone of the l-th unit, P l,t is the output of the l-th unit in the t period, is the allowable continuous start-stop time of the l-th unit; The double - layer dynamic programming unit load distribution model first uses the dynamic programming algorithm to calculate the optimal load distribution value and the corresponding consumption flow corresponding to all unit selection combinations within a single time period, repeats the calculation for all time periods, and obtains the optimal load distribution solution of the unit selection combination within different time periods; second, it uses the dynamic programming algorithm to calculate the unit start - stop optimization between different time periods, repeats the calculation for all units, and obtains the optimal unit selection combination between time periods; The specific process of using the dynamic programming algorithm to calculate the optimal load distribution value and the corresponding consumption flow corresponding to all unit selection combinations within a single time period is as follows: Set the variables of the dynamic programming algorithm. Among them, the number of operating units in the current unit selection combination is used as the stage variable, the total input load is used as the state variable, and the load borne by a single unit in the current unit selection combination is used as the decision variable; Obtain the corresponding power generation consumption flow according to the decision variable; Set the objective function. Among them, the objective function is set to find the optimal load distribution method in the load discrete space to minimize the total power generation consumption flow within the current time period, and the value range excludes the unit vibration zone; The specific process of using the dynamic programming algorithm to calculate the unit start - stop optimization between different time periods is as follows: Set the variables of the dynamic programming algorithm. Among them, different time periods are used as the stage variable, the total number of unit selection combinations is used as the state space, and the current unit selection combination is used as the decision variable; Set the objective function. Among them, the objective function is set to minimize the total power generation consumption flow within the current time period; The technical and economic evaluation indicators include: total water consumption, the number and proportion of non - recommended operation area time periods, additional start - stop times, output fluctuation ratio, average efficiency.

2. A storage medium, characterized in that, Stores a computer program or instruction, and when the computer program or instruction is run, it implements the hydraulic unit selection method based on vibration zone limitation and technical and economic indicators as described in claim 1.

Images