CN113270901B

CN113270901B - Multi-energy power supply control method giving consideration to battery state and active power dynamic stability

Info

Publication number: CN113270901B
Application number: CN202110667723.7A
Authority: CN
Inventors: 胡林; 李红刚; 李铁山; 师碧; 钟东明; 陶卿; 王军; 王昱倩; 魏超; 张攀全
Original assignee: Huaneng Lancang River Hydropower Co Ltd
Current assignee: Huaneng Lancang River Hydropower Co Ltd
Priority date: 2021-06-16
Filing date: 2021-06-16
Publication date: 2022-04-08
Anticipated expiration: 2041-06-16
Also published as: CN113270901A

Abstract

The invention discloses a multi-energy power supply control method giving consideration to both battery state and active power dynamic stability, which carries out coordination control on a thermal power supply, an energy storage power supply and a photovoltaic energy source through a multi-energy complementary integrated power supply centralized control center: the complementary integration unit sends out instructions including an instruction for distributing unit active power target values of the thermal power supply unit and the energy storage power supply unit, an instruction for setting a primary frequency modulation adjustment coefficient of the thermal power supply unit and an instruction for generating an on-off operation suggestion of the photovoltaic power supply unit to various types of power supply units. According to the invention, flexible charging and discharging of the energy storage power supply battery are realized by introducing charging and discharging correction power into the active power target value of the thermal power supply unit. On one hand, the battery state of each unit of the energy storage power supply is introduced into the calculation of the adjustment coefficient of the energy storage unit, on the other hand, a control strategy for preventing the adjustment coefficient of each energy storage unit from changing violently is designed, and meanwhile, the requirements of the battery state of each unit and the dynamic stability of active power in the adjustment process are considered.

Description

Multi-energy power supply control method giving consideration to battery state and active power dynamic stability

Technical Field

The invention belongs to the technical field of automatic control of power systems, and relates to a multi-energy power supply control method giving consideration to both battery state and active power dynamic stability.

Background

With the implementation of new energy strategies, the proportion of photovoltaic power generation is continuously increased, but the photovoltaic power generation mainly depends on weather resources which cannot be regulated and stored, the power generation capacity strongly depends on the weather resources which cannot be regulated and stored, the characteristics of strong randomness and volatility are provided, and the safety of a power grid is seriously threatened.

The thermal power supply takes combustion heat energy of coal and natural gas as a motive power source of the generator, so that the thermal power supply has good adjustability and storability (depending on coal storage amount and gas storage amount) compared with photovoltaic power generation, and is a core support power supply of the power system up to now together with the hydraulic power supply.

The unbalance between the power generation power and the power consumption of the power grid is represented by the frequency of the power grid and the rated frequency (50H)_z) When the deviation between the power grid frequency and the rated frequency exceeds a threshold value, the output active power of each grid-connected power station in the control range is regulated by scheduling, so that the power generation power and the consumed power of the power grid are restored to a balanced state, the difference between the power grid frequency and the rated frequency is ensured to be in an allowable range, and the whole process is called secondary frequency modulation. The secondary frequency modulation comprises the following steps: 1) the dispatching mechanism calculates the generating power variation required for enabling the power grid frequency to return to the rated frequency according to the power grid frequency deviation and the power grid frequency-power sensitivity coefficient; 2) the dispatching corrects the active power set value of each grid-connected power station in the control area according to the calculation result, and sends a power regulation instruction; 3) after each power station receives the new active power set value, the AGC distributes the total active power set value of the power station to each unit controlled by the AGC; 4) and the active power control system of each unit performs closed-loop feedback regulation on the active power of the unit according to the new single-unit active power set value.

When the deviation of the power grid frequency and the rated frequency exceeds a primary frequency modulation threshold value (most of domestic power grids are thermal power 0.03H)_z) Then, the speed regulator system of each unit regulates the active power of the unit according to a preset 'frequency-power' regulation coefficient so as to regulate the active power of the unitAnd the unbalance between the generated power and the consumed power of the power grid is compensated to a certain extent. Compared with secondary frequency modulation, because a unified control center is not provided for performing coordinated control on each unit participating in primary frequency modulation and is related to a calculation mechanism of an adjustment amount, the primary frequency modulation cannot enable the frequency of a power grid to be completely recovered to a rated frequency, so that the primary frequency modulation is also called as differential adjustment, but the primary frequency modulation has the advantages that: 1) because a uniform control center is not arranged, the risk of complete failure like secondary frequency modulation (for example, abnormal exit of a secondary frequency modulation function module is scheduled) is avoided, and thus extremely high overall reliability is obtained; 2) the regulating instruction is directly calculated by the unit, and processes of scheduling calculation, instruction transmission, AGC distribution of a power station and the like of secondary frequency modulation are omitted, so that the response speed to the power grid frequency abnormality is far higher than that of the secondary frequency modulation.

The photovoltaic power supply and the thermal power supply are taken as an organic whole, and the power regulation task is executed aiming at the dynamic balance of consumption and supply of the power system. Compared with a single photovoltaic power supply or a single thermal power supply, the solar photovoltaic power supply has the advantages that the solar photovoltaic power supply has the adjusting capacity equivalent to the scale of the thermal power supply, and meanwhile, under the condition of abundant sunlight, the active power output of the thermal power supply can be correspondingly reduced, so that the coal-saving and energy-saving aim is realized; however, the limitation is that the photovoltaic and thermal power supply has the performance disadvantage of secondary frequency modulation as the thermal power supply, and due to the existence of the intrinsic active power regulation delay of the thermal power supply, the photovoltaic and thermal power supply can only be restrained to a certain extent and cannot solve the problem of random fluctuation of the output power of the photovoltaic power supply, and in an extreme case, when the output power of the photovoltaic power supply oscillates in a similar simple harmonic wave, the thermal power supply may even have resonant regulation of the active power due to the regulation delay, so that the overall output power oscillation of the photovoltaic and thermal power supply multifunctional complementary integrated power supply is intensified.

Disclosure of Invention

The invention solves the technical problem of providing a multi-energy power supply control method giving consideration to both the battery state and the dynamic stability of active power, and the power supply with better regulation performance is utilized to carry out compensation regulation on the power supply with poorer regulation performance or the power supply without regulation capacity so as to improve the total regulation performance of the active power of the multi-energy complementary integrated power supply.

The invention is realized by the following technical scheme:

a multi-energy power supply control method giving consideration to both battery state and active power dynamic stability carries out coordination control on a thermal power supply, an energy storage power supply and a photovoltaic energy source through a multi-energy complementary integrated power supply centralized control center:

the multi-energy complementary integrated power supply centralized control center is provided with a complementary integrated unit, a firepower power supply unit, an energy storage power supply unit and a photovoltaic power supply unit; the complementary integration unit sends an instruction for distributing a unit active power target value of the firepower power supply unit and an instruction for setting a primary frequency modulation adjustment coefficient of the firepower supply unit to the firepower supply unit; the complementary integration unit sends an instruction for distributing the unit active power target value of the energy storage power supply unit to the energy storage power supply unit; the complementary integration unit sends an instruction of a start-up and shut-down operation suggestion of the photovoltaic power supply unit to the photovoltaic power supply unit;

the complementary integrated unit distributes the unit active power target value of the firepower power supply unit as follows: the unit active power target value of the thermal power supply is equal to the sum of the total active power set value of the multi-energy complementary integrated power supply minus the unit active power real-emitting value of the photovoltaic power supply unit and the calculated value filtering value, and the charging and discharging correction power of the energy storage power supply unit;

the real unit active power value of the photovoltaic power supply unit is involved in the calculated quantity filtering value and is updated according to a fixed period based on the real unit active power value of the photovoltaic power supply unit, a filtering threshold and a dead output area of the photovoltaic power supply unit;

the charging and discharging correction power of the energy storage power supply unit is periodically updated by the energy storage power supply unit according to the battery electric quantity state, the battery capacity and the charging and discharging coefficient of each energy storage unit;

complementary integrated unit all transfers the primary frequency modulation task of photovoltaic power supply unit to firepower power supply unit, and its setting to firepower power supply unit's primary frequency modulation adjustment coefficient is: multiplying a primary frequency modulation adjustment coefficient of a thermal power supply unit issued by a power grid by a primary frequency modulation scaling coefficient; the primary frequency modulation scaling coefficient is equal to (rated capacity of active power of the photovoltaic power supply unit + rated capacity of active power of the fire power supply unit) ÷ rated capacity of active power of the fire power supply unit;

the complementary integrated unit distributes the unit active power target value of the energy storage power supply unit as follows: adding a total active power set value of the multi-energy complementary integrated power supply to a primary frequency modulation correction quantity of a unit of a thermal power supply unit, then subtracting a real unit active power value of the photovoltaic power supply unit from a calculated quantity, and then subtracting a real unit active power value of the thermal power supply unit to obtain a total active power adjustment deviation; updating the unit active power target value of the energy storage power supply unit according to the total active power regulation deviation and a fixed period;

the real unit active power value of the photovoltaic power supply unit is involved in the calculated quantity, and is updated according to a fixed period based on the real unit active power value of the photovoltaic power supply unit and the output dead zone of the photovoltaic power supply unit;

the complementary integrated unit sends out the operation suggestion of starting and stopping the machine of the photovoltaic power supply unit to the photovoltaic power supply unit: the complementary integration unit generates an operation suggestion on the startup and shutdown of the photovoltaic power supply unit according to the mismatch quantitative value of the total active power set value of the multi-energy complementary integrated power supply and by combining the current startup and shutdown sequence of the photovoltaic power supply unit and the possible active power fluctuation range sequence respectively corresponding to the startup and shutdown sequence;

the complementary integration unit acquires a mismatching degree quantization value of the total active power set value of the multi-energy complementary integrated power supply according to the possible active power fluctuation range of the photovoltaic power supply unit, the total active power set value of the multi-energy complementary integrated power supply and the unit combined operation area of the thermal power supply unit;

the thermal power supply unit obtains thermal power supply control intermediate parameters according to basic parameters of a thermal power supply and sends the intermediate parameters to the complementary integrated unit, and performs unit-level AGC distribution and unit active power closed-loop regulation on the thermal power supply according to the received unit active power target value and primary frequency modulation regulation coefficient, and sends an operation suggestion of the thermal power supply unit;

the energy storage power supply unit obtains energy storage power supply control intermediate parameters including a battery state and an energy storage power supply regulation coefficient according to basic parameters of the energy storage power supply and sends the energy storage power supply control intermediate parameters to the complementary integration unit, and unit-level AGC distribution and unit active power closed-loop regulation of the energy storage power supply are carried out according to the received unit active power target value;

the photovoltaic power supply unit sends the photovoltaic power supply control intermediate parameters to the complementary integration unit; and sending a suggested command of the start-up and shutdown operation of the photovoltaic generator set.

Compared with the prior art, the invention has the following beneficial technical effects:

the invention takes the emerging energy storage power supply, particularly the shallow charging and shallow discharging problem of the energy storage power supply battery as the focus of attention, and when the charging and discharging strategy of the battery is designed, the firepower supply is used as the charging and discharging source of the energy storage power supply battery by introducing charging and discharging correction power into the active power target value of the firepower supply unit; meanwhile, the battery states of each unit of the energy storage power supply are introduced into the calculation of the adjustment coefficients of the energy storage units, and a control strategy for preventing the adjustment coefficients of each energy storage unit from changing violently is designed so as to simultaneously take the requirements of the battery states of each unit and the dynamic stability of the active power in the adjustment process into consideration;

further using a fire power supply to perform compensation adjustment on the photovoltaic power supply, using an energy storage power supply to perform compensation adjustment on the fire power supply and the photovoltaic power supply unit, and using a power supply with better adjustment performance to perform compensation adjustment on a power supply with poorer adjustment performance or a power supply without adjustment capacity so as to improve the active power overall adjustment performance of the multi-energy complementary integrated power supply;

the invention also aims at the nonideal of the adjusting process and the adjusting result caused by the problems of time delay, precision and the like of the adjustment of various power supplies, and a large number of parameters such as operation dead zones and the like are introduced into the active power control strategy so as to inhibit the overall sensitivity of the control strategy and prevent the problems of overhigh calculation frequency, frequent change of an adjusting target, excessive compensation and the like; and further, by introducing the method that the real active power value of the photovoltaic power supply unit participates in the calculated amount filtering value into the multi-energy complementary integrated power supply, the sensitivity of the fire power supply to the random fluctuation of the active power of the photovoltaic power supply unit is reduced, and the fire power supply corrects the large-amplitude deviation of the real active power value of the photovoltaic power supply unit.

Drawings

FIG. 1 is a simulation modeling diagram of a multi-energy complementary integrated power supply of the invention, namely a thermal power supply, a photovoltaic power supply and an energy storage power supply;

FIG. 2 is a block diagram of the logic framework for computing and controlling the energy storage power supply unit according to the present invention;

FIG. 3 is a logic flow diagram of calculating the adjustment coefficients of the energy storage units of the energy storage power supply unit according to the present invention;

FIG. 4 is a schematic diagram showing the relationship between the proportional variation of the battery SOC capacity for the upward and downward adjustment validation threshold parameters of each energy storage unit according to the present invention;

FIG. 5 is a logic diagram illustrating the dead zone processing of the active power target value of the energy storage power supply unit in the multi-energy complementary integrated power supply according to the present invention;

FIG. 6 is a diagram of the adjusting effect of the multi-energy complementary integrated power supply of the invention, namely a thermal power supply, a photovoltaic power supply and an energy storage power supply;

FIG. 7 is a simulation modeling diagram II of the multi-energy complementary integrated power supply of "thermal power supply + photovoltaic + energy storage power supply" according to the present invention;

FIG. 8 is an adjustment effect diagram of a simulation modeling diagram II of a multi-energy complementary integrated power supply of a thermal power supply, a photovoltaic power supply and an energy storage power supply.

Detailed Description

The present invention will now be described in further detail with reference to the following examples, which are intended to be illustrative, but not limiting, of the invention.

The following describes each unit in detail.

S1000), the parameters acquired by the complementary integration unit comprise:

s1100), parameters input by a complementary integration unit:

s1111) directly inputting a total active power set value of the multi-energy complementary integrated power supply;

s1112) unit active power rated capacity, wherein the unit active power rated capacity of the thermal power supply and the photovoltaic power supply is equal to the sum of single machine active power rated capacities of the units which are generating power by the power supply units, and the unit active power rated capacity of the energy storage power supply depends on the rated capacity of each energy storage unit and the charge state of a battery;

s1113) the real active power value of the unit is respectively equal to the sum of the real active power values of the units of the thermal power supply unit, the energy storage power supply unit and the photovoltaic power supply unit;

s1114) unit active power regulation dead zones are respectively equal to the sum of single-machine active power regulation dead zones of the units in which the thermal power supply unit and the energy storage power supply unit are running;

s1120) input parameters transmitted by the thermal power supply unit:

s1121) the unit primary frequency modulation target regulating quantity of the thermal power supply unit is equal to the sum of the single-machine primary frequency modulation target regulating quantities of the generating set;

s1122) a unit joint operation area of the thermal power supply unit;

s1123) carrying out primary frequency modulation actual adjustment quantity on the unit of the thermal power supply unit;

s1124) correcting the primary frequency modulation of the unit of the thermal power supply unit, when the primary frequency modulation actual regulating quantity of each unit of the thermal power supply unit can be measured, the correction is equal to the primary frequency modulation actual regulating quantity of the unit of the thermal power supply unit, otherwise, the correction is equal to the primary frequency modulation target regulating quantity of the unit of the thermal power supply unit in S1121;

s1130) parameters sent by the energy storage power supply unit: the charging and discharging correction power of the energy storage power supply unit is calculated by the energy storage power supply unit according to parameters such as battery states of the energy storage units;

s1140) input parameters sent by the photovoltaic power supply unit:

s1141) the real unit active power value of the photovoltaic power supply unit participates in the calculated quantity, and the photovoltaic power supply unit updates according to the real unit active power value and the output dead zone of each photovoltaic unit according to a fixed period;

s1142) the real unit active power value of the photovoltaic power supply unit is involved in the calculated value of the filtered value, and the photovoltaic power supply unit updates according to the real unit active power value, the scaling coefficient and the dead output area of each photovoltaic unit according to a fixed period; (ii) a

S1143) the possible fluctuation range of the active power of the photovoltaic power supply unit is a prediction result of the fluctuation range of the active power of the photovoltaic power supply unit within a certain time in the future;

s1144) a starting sequence and a stopping sequence of the photovoltaic power supply unit, and active power possible fluctuation range sequences respectively corresponding to the starting sequence and the stopping sequence are used for generating a starting and stopping operation suggestion for the photovoltaic unit;

s1145) the unit primary frequency modulation target regulating quantity of the photovoltaic power supply unit is equal to the sum of the single-machine primary frequency modulation target regulating quantity of the photovoltaic unit generating electricity;

the instruction of the operation suggestion of the start-up and shutdown of the photovoltaic power supply unit is obtained according to a total active power set value of the multi-energy complementary integrated power supply, a unit combined operation area of the thermal power supply unit, a possible active power fluctuation range of the photovoltaic power supply unit, a start-up and shutdown sequence of the photovoltaic power supply unit and a possible active power fluctuation range sequence respectively corresponding to the start-up and shutdown sequence, and the operation suggestion of the start-up and shutdown of the photovoltaic power supply unit, which is referred by an operator, is generated.

S2000) operation of the fire power supply is specifically given below.

S2100) determining the unit type of the thermal power supply unit, including:

s2110) dividing the generator set and the non-generator set according to different generator set states, wherein the non-generator set comprises a generator set in a shutdown state, an idle state and an unsteady state;

s2120) according to the difference of the active power regulation controlled state of the generator set, further dividing the generator set:

s2131) a single-machine open-loop unit, namely a unit of which the single-machine active power actual output value is not regulated by any source;

s2132) a single-machine closed-loop unit, namely, a single-machine active power real output value of the unit is subjected to closed-loop regulation according to a single-machine active power set value or an execution value, so that the single-machine active power real output value of the unit continuously tends to the single-machine active power set value or the execution value and is finally stabilized in the unit within a single-machine active power set value or execution value regulation dead zone range;

s2133) putting into an AGC unit, namely a unit closed loop, wherein the unit active power set value of the unit is distributed and set by unit-level AGC;

s2134) the unit which is not put into AGC, namely the generator units except the unit which is put into AGC, comprise a single-machine open-loop unit and a single-machine closed-loop unit which does not accept unit-level AGC distribution and setting of the active power set value of the single machine;

s2200) establishing a combined output model for each unit of the AGC, and calculating a combined operation area and a combined recommended operation area, wherein the combined output model comprises the following steps:

s2210) determining a single machine recommended operation area and a single machine forbidden operation area which are put into each unit of the AGC, comprising the following steps:

s2211) a stand-alone operation forbidden area refers to a load area in which the set value of the stand-alone active power of the unit is forbidden to be set (between the upper limit and the lower limit of the stand-alone operation forbidden area); the real value of the single-machine active power of the unit is allowed to pass through or pass through the single-machine forbidden operation area, but is not allowed to reside or stay in the single-machine forbidden operation area for a long time;

s2212) the single machine suggested operation area is a load area with high unit operation efficiency and stable operation when the single machine active power actual value of the unit is between the upper limit and the lower limit of the single machine suggested operation area; under the condition that the conditions allow, the single machine active power set value of the unit is preferably set in a single machine suggested operation area;

s2215) the low-load area of the conventional thermal power unit is a single-machine forbidden operation area, the single-machine forbidden operation area of the thermal power unit is about 0-50% of rated capacity, and the rest part of the rated capacity minus the single-machine forbidden operation area is a single-machine suggested operation area;

s2220) establishing a suggested combined output model of the unit which is put into the AGC, and calculating a combined suggested operation area which is put into the AGC unit, wherein the method comprises the following steps: adding the single machine suggested operation areas put into the AGC unit to obtain a combined suggested operation area put into the AGC unit;

s2240) determining the current single-machine AGC active power distribution value of each unit, including:

s2241) for the unit which is put into the AGC, the unit AGC active power distribution value is distributed by the unit-level AGC;

s2242) for a single-machine closed-loop unit which is not put into AGC, tracking a single-machine active power set value by a single-machine AGC active power distribution value;

s2243) for the single-machine open-loop unit which is not put into the AGC, the single-machine AGC active power distribution value tracks the single-machine active power set value, and the single-machine active power set value is assigned by the single-machine active power real sending value, namely when the single-machine active power set value is not equal to the single-machine active power real sending value and the absolute value of the difference between the single-machine active power set value and the single-machine active power set value is larger than the single-machine active power regulation dead zone, the single-machine active power real sending value is written into the single-machine active power set value.

S2250) adding the joint recommended operation area obtained in S2220 and input into the AGC unit and the active power distribution values of all stand-alone AGC units not input into the AGC unit to obtain a unit joint recommended operation area of the thermal power supply and a unit joint operation area of the thermal power supply, and providing reference for automatic active power control of the thermal power supply unit and comprehensive control of the multi-energy complementary integrated power supply;

s2300) comparing a unit active power target value of the thermal power supply with the unit combined operation area in the S2250, and skipping the rest step of the S2300 if the unit active power target value is feasible when the unit active power target value is included in the unit combined operation area; when the unit active power target value is not included in the unit joint operation area and the unit active power target value is not feasible, searching an operation proposal for enabling the unit active power target value:

s2320) finding an operation proposal for making a unit active power target value of the thermal power supply feasible by putting a unit not put into AGC control, including:

s2321) setting a loop variable i₁，i₁Is set to 1;

s2322) for i₁Making a judgment if i₁If the number of the units not put into the AGC is larger than the number of the units not put into the AGC, the S2320 is terminated, otherwise, the following steps are continuously executed to find the number of the units i₁The unit which is not put into AGC is put into AGC control so as to make the unit active power target value of the thermal power supply feasible;

s2323) listing and selecting i from all the units which are not put into AGC₁All combinations of stages, C (j)₁,i₁) Wherein C () is a combination number function, j₁The number of the units which are not put into AGC;

s2324) C (j) listed respectively as S2323₁,i₁) In the combination mode, a unit which is selected in various modes and is not put into AGC is assumed to be put into AGC, a unit joint operation area and a unit joint suggested operation area are calculated by adopting the S2200 method, and then the feasibility of the unit active power target value is judged again by adopting the S2300 method according to the newly calculated unit joint operation area;

s2325) according to the calculation result of S2324, if there is and only 1 unit joint operation area regenerated in the mode can make the unit active power target value feasible, generating operation suggestions, namely 'putting the unit selected in the mode and not put into AGC', if there are unit joint operation areas regenerated in the modes can make the unit active power target value feasible, respectively generating operation suggestions according to the modes, namely 'putting the unit selected in the corresponding mode and not put into AGC', and jumping to step S2326 to continue execution, and if there is no unit joint operation area regenerated in any mode can make the unit active power target value feasible, i₁＝i₁+1, then go to step S2322 for i₁And judging whether the number of the units not put into the AGC is larger than that of the units not put into the AGC, and determining whether to execute the subsequent steps according to the judgment result.

S2326) carrying out priority sequencing on the multiple operation suggestions generated in the step S2325, wherein the sequencing is based on the absolute value of the difference value (the larger the difference value, the better the difference value) between the unit active power target value and the unit combined operation zone boundary.

S2330) finding a running operation recommendation that makes a unit active power target value of the thermal power source feasible by turning the non-power generating unit to the power generating state and putting into AGC, including:

s2331) setting a circulation variable i₂，i₂Is set to 1;

s2332) pairs of i₂Making a judgment if i₂If the number of the units which are available and do not generate electricity is larger than the number of the units which are available and do not generate electricity, the step S2330 is terminated, otherwise, the following steps are continuously executed to search for the unit i₂The unit which can be used and does not generate power is converted into a power generation state and is put into AGC to make the unit active power target value of the thermal power supply feasible;

s2333) enumerating the selection of i from all available and unenergized units₂All combinations of stages, C (j)₂,i₂) Wherein j is₂The number of units which are available and not generating electricity;

s2334) C (j) listed according to S2333, respectively₂,i₂) A combination mode is adopted, available and non-power generation units selected in various modes are assumed to be in a power generation state and are put into AGC, a unit joint operation area and a unit joint recommended operation area are calculated by adopting the S2200 method again, and then the feasibility of the unit active power target value is judged again by adopting the S2300 method according to the newly calculated unit joint operation area;

s2335) according to the calculation result of S2334, if there are and only 1 unit joint operation area regenerated by 1 mode to enable the unit active power target value, generating operation suggestions, namely converting the available and non-power generation unit selected by the mode into the power generation state and putting the unit into AGC, if there are unit joint operation areas regenerated by multiple modes to enable the unit active power target value, respectively generating operation suggestions according to the modes, namely converting the available and non-power generation unit selected by the corresponding mode into the power generation state and putting the unit into AGC, and jumping to step S2336 to continue execution, if there is no unit joint operation area regenerated by any mode to enable the unit active power target value, i₂＝i₂+1, and then go to step S2332 for i₂Whether greater than available and not generating electricityAnd judging the number of the units, and determining whether to execute the subsequent steps according to the judgment result.

S2336) carrying out priority ranking on the plurality of operation suggestions generated in the S2335, wherein the ranking is carried out according to the absolute value of the difference value of the unit active power target value and the unit combined operation area boundary (the larger the absolute value is, the better the absolute value is).

S2340) finding a running operation recommendation that makes a unit active power target value of the thermal power source feasible by turning the generating unit to the non-generating state, includes:

s2341) setting a Loop variable i₃，i₃Is set to 1;

s2342) pairs of i₃Making a judgment if i₃If the number of the generating units is larger than the number of the generating units, S2340 is ended, otherwise, the following steps are continuously executed to find the number i of the generating units₃The unit of the platform power generation is changed into a non-power generation state, so that the unit active power target value of the thermal power supply becomes feasible;

s2343) listing and selecting i from all power generation units₃All combinations of stages, C (j)₃,i₃) Wherein j is₃The number of generating units;

s2344) C (j) listed according to S2343, respectively₃,i₃) In the combination mode, the unit for generating power selected by various modes is assumed to be in a non-power generation state, the unit combined operation area and the unit combined suggested operation area are calculated by adopting the S2200 method, and then the feasibility of the unit active power target value is judged again by adopting the S2300 method according to the newly calculated unit combined operation area;

s2345) according to the calculation result of S2344, if the unit active power target value can be enabled by the unit joint operation area regenerated in only 1 mode, an operation suggestion that 'the generating set selected in the mode is converted into the non-generating state' is generated, if the unit active power target value can be enabled by the unit joint operation area regenerated in multiple modes, operation suggestions that 'the generating set selected in the corresponding mode is converted into the non-generating state' are respectively generated according to the modes, and the operation is continued by jumping to the step S2346, if no one mode existsThe regenerated unit joint operation area can enable the unit active power target value to be feasible, i₃＝i₃+1, and then go to step S2342 for i₃And judging whether the number of the units is larger than the number of the generating sets or not, and determining whether to execute the subsequent steps or not according to a judgment result.

S2346) carrying out priority ranking on the multiple operation suggestions generated in the S2345 according to the fact that the operation suggestions are selected from the generating set i correspondingly₃The combination mode of the station set and the range of the unit joint operation area and the unit joint recommended operation area after the operation suggestions obtained in step S2344 are changed correspondingly, and the ranking is respectively from high to low according to the importance degree: and selecting the number of the units which are not subjected to AGC (the more units are the better) and the number of the units which are subjected to AGC (the less units are the better), and the absolute value of the difference value of the unit active power target value and the unit combined operation area boundary (the larger unit is the better).

S2350) classifying the operation suggestions generated by the S2320, the S2330 and the S2340, and orderly displaying the operation suggestions according to the priorities (when more than 1 operation suggestion in a certain class is obtained) obtained by the S2326, the S2336 and the S2346 so as to assist the decision of an operator.

S2400) calculating a single AGC active power distribution value which is put into an AGC unit, wherein the calculation comprises the following steps:

s2410) calculating a unit AGC active power distribution value of the thermal power supply, including:

s2411) calculating the active power distribution values of all single AGC units which are not put into the AGC unit, wherein the obtaining mode of the active power distribution values of the single AGC units is as the mode of S2240;

s2412) subtracting all single AGC active power distribution values which are not put into the AGC unit from the unit active power target value to obtain a unit AGC active power distribution value.

S2450) according to the unit AGC active power distribution value, carrying out AGC active power average distribution on each unit which is put into the AGC;

s2500) active power regulation of each single closed-loop unit of the thermal power supply unit, comprising:

s2510) determining the single-machine active power setting value of each single-machine closed-loop unit, including:

s2511) for a stand-alone closed-loop unit which is not put into AGC, setting a stand-alone active power set value manually by an operator;

s2512) for the thermal power generating unit which is put into the AGC, the single-machine active power set value is equal to the single-machine AGC active power distribution value.

S2530) the active power control system of each single-machine closed-loop unit of the thermal power supply unit takes a single-machine active power set value as a target, calculates the deviation between the single-machine active power actual output value and the single-machine active power set value, and outputs continuous signals according to the calculation result to adjust the single-machine active power actual output value of the unit so that the single-machine active power actual output value of the unit tends to the single-machine active power set value and is finally stabilized in the adjustment dead zone range of the single-machine active power set value.

S3000) the operation of the energy storage power supply unit is given below, and the calculation and control logic thereof is as shown in fig. 2, and includes:

s3100) calculating the charged capacity proportion of each energy storage unit battery of the energy storage power supply unit and the charged total capacity proportion of the energy storage power supply unit batteries, and the method comprises the following steps:

s3110) calculating the charge capacity ratio of each energy storage unit battery,

in the formula r_iFor the battery state of charge capacity ratio, SOC, of the energy storage unit i_iIs the battery charge state of the energy storage unit i,

and

respectively maximum and minimum of battery charge of energy storage unit i, e.g. when certain energy storage unit SOC_i＝50，

And

are respectively 100 and10, then

S3120) calculating the overall capacity ratio of the energy storage power source unit cell charge,

where r is the total charge capacity ratio of the energy storage power supply units, for example, if a certain energy storage power supply unit comprises 3 energy storage units, the charge states of the batteries are 40, 50 and 60 respectively, the maximum charge states of the batteries are 100, 110 and 120 respectively, and the minimum charge states of the batteries are 0,5 and 10 respectively, then

S3200) setting judgment threshold value R of total capacity proportion of energy storage power source unit battery state of charge₁’～R₆' the setting principle comprises:

S3210)0＜R₁’＜R₂’＜R₃’＜R₄’＜R₅’＜R₆’＜1；

S3220)R₁’+R₆’＝1；

S3230)R₂’+R₅’＝1；

S3230)R₃’+R₄’＝1。

in this example, R is₁’～R₆' is set to 20%, 30%, 45%, 55%, 70%, 80%, respectively.

S3300) judging the overall electric quantity state of the battery of the energy storage power supply unit, including:

s3310) when S3120 is reached, the overall capacity ratio of the energy storage power source unit cell charge is not less than 0_r＜R₁When the battery of the energy storage power supply unit is in an extremely low power state;

s3320) when R₁’≤_r＜R₂When the battery of the energy storage power supply unit is in a lower state of charge;

s3330) when R₂’≤_r＜R₃' or R₄’＜_r≤R₅When the battery is in a more ideal electric quantity state, the whole battery of the energy storage power supply unit is in a more ideal electric quantity state;

s3340) when R₃’≤_r≤R₄When the battery of the energy storage power supply unit is in an extremely ideal electric quantity state;

s3350) when R₅’＜_r≤R₆When the battery of the energy storage power supply unit is in a higher state of charge;

s3360) when R₆’＜_rWhen the battery capacity is less than or equal to 1, the battery of the energy storage power supply unit is in an extremely high electric quantity state.

S3400) setting judgment threshold R of battery state-of-charge capacity ratio of energy storage unit₁～R₄The setting principle comprises the following steps:

S3210)0＜R₁＜R₂＜R₃＜R₄＜1；

S3220)R₁+R₄＝1；

S3230)R₂+R₃＝1。

this example will show that₁～R₄Set to 20%, 40%, 60%, 80%, respectively.

S3500) auxiliary calculation parameters of the adjustment coefficients of the energy storage units of the energy storage power supply unit are set, and the method comprises the following steps:

s3510) setting 4 threshold parameters K₁、K₂、K₃、K₄Wherein 0 < K₁＜K₂＜K₃＜K₄This example compares K with₁～K₄Respectively set to 0.5, 1, 1.5 and 2;

s3520) setting a variable gradient parameter delta K of an energy storage unit adjusting coefficient, wherein delta K is more than 0 and less than min K₁,K₂－K₁,K₃－K₂,K₄－K₃]Wherein min 2]In order to take a minimum function, setting Δ K is to prevent the dynamic stability of the real power value of the unit from being reduced due to too severe change of the adjustment coefficient of the energy storage unit in the adjustment process, and in this embodiment, Δ K is set to 0.1;

s3600) calculating an adjustment coefficient of each energy storage unit of the energy storage power supply unit, as shown in fig. 3, including:

s3610) calculating upward adjustment coefficients of each energy storage unit of the energy storage power supply unit, and the method comprises the following steps:

s3611) upward adjustment coefficients of energy storage units of the energy storage power supply unit are initially set

In the formula

The upward adjustment coefficient of the energy storage unit i is obtained;

s3612) correcting the upward adjustment coefficients of the energy storage units according to a fixed period, namely continuously and circularly operating the subsequent steps according to the fixed period;

s3613) calculating effective threshold parameters of upward adjustment of each energy storage unit

When the content is less than or equal to 0_ri＜R₁Time of flight

When R is₁≤_ri＜R₂Time of flight

When R is₂≤_ri≤R₃Time of flight

When R is₃＜_ri≤R₄Time of flight

When R is₄＜_riWhen the temperature is less than or equal to 1

S3614) comparison

And

when the absolute value of the difference between the two is less than or equal to delta K

When the absolute value of the difference between the two is greater than delta K and

time of flight

time of flight

E.g. certain energy storage units

And

all are equal to 1 originally, and the proportion of the charge capacity of the battery is reduced to R in the discharging process₁And R₂In between, thus

With a consequent decrease of 0.5, so that in the next few cycles,

corrected to 0.9, 0.8, 0.7, 0.6, 0.5, respectively.

S3620) calculating downward adjustment coefficients of each energy storage unit of the energy storage power supply unit, and the method comprises the following steps:

s3621) initializing and setting downward adjustment coefficients of energy storage units of the energy storage power supply unit

In the formula

The downward adjustment coefficient of the energy storage unit i is obtained;

s3622) correcting the downward adjustment coefficients of the energy storage units according to a fixed period, namely continuously and circularly operating the subsequent steps according to the fixed period;

s3623) calculating effective threshold parameters of downward adjustment of each energy storage unit _ikWhen 0 is less than or equal to r_i＜R₁Time of flight _ik＝K₄When R is₁≤r_i＜R₂Time of flight _ik＝K₃When R is₂≤r_i≤R₃Time of flight _ik＝K₂When R is₃＜r_i≤R₄Time of flight _ik＝K₁When R is₄＜r_iWhen the temperature is less than or equal to 1 _ik＝0；

S3624) comparison

And _ikwhen the absolute value of the difference between the two is less than or equal to Δ K

time of flight

time of flight

In the above embodiment, the capacity ratio r is determined according to the state of charge of the battery_iEffective threshold value parameter for upward and downward adjustment of energy storage units

_ikAs shown in fig. 4, as the soc capacity ratio of the battery increases, the upward adjustment effective threshold parameter of the energy storage unit increases, and the downward adjustment effective threshold parameter decreases, and since the upward adjustment coefficient and the downward adjustment coefficient of the energy storage unit respectively tend to change according to the upward adjustment effective threshold parameter and the downward adjustment effective threshold parameter, the upward adjustment coefficient and the downward adjustment coefficient of the energy storage unit also increase and decrease according to the increase of the soc capacity ratio of the battery.

S3700) performing unit-level AGC distribution on the unit active power target value of the energy storage power supply unit, including:

s3710) when the unit active power target value of the energy storage power supply unit is equal to 0, the single-machine active power set value of each energy storage unit is equal to 0;

s3720) when the unit active power target value of the energy storage power supply unit is larger than 0, the single active power set value of each energy storage unit is distributed according to the mutual proportion of the product of the upward adjustment coefficient and the battery capacity of each energy storage unit, namely the single active power set value of each energy storage unit is equal to

In the formula

The active power target value of the unit of the energy storage power supply unit is obtained, if the calculation result is larger than the positive single-machine active power rated capacity of the energy storage unit, the positive single-machine active power rated capacity of the energy storage unit is used as a single-machine active power set value, and if the active power target value of the energy storage power supply unit is 300MW, 3 energy storage units are provided

0.5, 1, 1.5, respectively, cell capacity

200MW, 150 MW, 220MW respectively, the single active power set values of 3 energy storage units are respectively

S3730) when the unit active power target value of the energy storage power supply unit is less than 0, the single active power set value of each energy storage unit is distributed according to the mutual proportion of the product of the downward adjustment coefficient and the battery capacity of each energy storage unit, namely the single active power set value of each energy storage unit is equal to

If the calculation result is smaller than the negative single-machine active power rated capacity of the energy storage unit, taking the negative single-machine active power rated capacity of the energy storage unit as a single-machine active power set value, and assuming that the active power target value of the energy storage power supply unit is-300 MW, 3 energy storage units exist

0.5, 1, 1.5, respectively, cell capacity

200, 150 and 220MW respectively, and the single active power set values of the 3 energy storage units are-51.7, -77.6 and-170.7 MW respectively.

As described in S3600, the upward adjustment coefficient and the downward adjustment coefficient of the energy storage unit respectively increase and decrease with the increase of the battery soc capacity ratio, so according to the calculation methods of S3720 and S3730, when the active power target value of the energy storage unit is greater than 0, that is, the energy storage unit is generally in a discharge state, the energy storage unit with a higher battery soc capacity ratio tends to discharge, and when the active power target value of the energy storage unit is less than 0, that is, the energy storage unit is generally in a charge state, the energy storage unit with a lower battery soc capacity ratio tends to charge, so that the soc capacity ratios of the energy storage units can be kept consistent, and the battery overcharge or overdischarge of one or several energy storage unit batteries compared with other energy storage unit batteries can be avoided.

S3800) an active power control system of each energy storage unit of the energy storage power supply unit calculates a deviation between a single-machine active power actual output value and the single-machine active power set value by taking a single-machine active power set value as a target, and outputs continuous signals according to a calculation result to adjust the single-machine active power actual output value of the energy storage unit so that the single-machine active power actual output value of the energy storage unit tends to the single-machine active power set value and is finally stabilized within an adjustment dead zone range of the single-machine active power set value.

S3900) calculating unit active power rated capacity of the energy storage power supply unit, and the method comprises the following steps:

s3910) calculating upward adjusting capacity of each energy storage unit of the energy storage power supply unit, wherein the upward adjusting capacity comprises the following steps:

s3911) when the energy storage unit is adjusted upwards as calculated in S3613, the effective threshold parameter

Then, the upward regulating capacity of the unit is the positive single-machine active power rated capacity of the unit;

s3912) when the energy storage unit is adjusted upwards as calculated in S3613, the effective threshold parameter is obtained

The upward regulating capacity of the unit is the product of the positive single-machine active power rated capacity of the unit and the rated capacity

Then divided by K₂；

For example, when the forward single-machine active power rated capacity of the unit is 50MW, if K is₂When 1, then

Are respectively 1.5, 1 and 0And 5, the upward adjusting capacity of the unit is respectively 50MW, 50MW and 25 MW.

S3920) accumulating the upward regulating energy of each energy storage unit obtained in the S3910 to obtain the active power rated capacity of the forward unit of the energy storage power supply unit;

s3930) calculating the downward adjusting capacity of each energy storage unit of the energy storage power supply unit, wherein the downward adjusting capacity comprises the following steps:

s3931) when the energy storage unit is adjusted downward as calculated in S3623, the effective threshold parameter _ik≥K₂Then, the downward regulating capacity of the unit is the negative single-machine active power rated capacity of the unit;

s3932) when the energy storage unit is adjusted downward as calculated in S3623, the effective threshold value parameter _ik＜K₂The downward regulating capacity of the unit is the product of the negative single-machine active power rated capacity of the unit and the rated capacity _ikThen divided by K₂。

S3940) accumulating the downward regulating energy of each energy storage unit obtained in the S3930 to obtain the negative direction unit active power rated capacity of the energy storage power supply unit.

S4000) the operation of the photovoltaic power supply unit is given below; the method specifically comprises the following steps:

s4100) aiming at the characteristics of non-adjustable active power, output power fluctuation and intermittence of the photovoltaic power supply, generating future T for each unit₁Possible fluctuation range of active power in time is calculated, and possible fluctuation range of unit active power of the photovoltaic power supply is calculated, wherein T₁For the manual setting parameter, the purpose is in order to reserve sufficient time for the perhaps start-up and shut-down operation of photovoltaic unit, include:

s4110) if a power prediction system is deployed, adopting future T of each photovoltaic unit output by a power prediction function₁The power prediction system is a system which adopts a physical method, a regression method, a time series method, a neural network method, a deep learning method and the like to establish a prediction model according to the past power, the contemporaneous historical data, the seasonal variation, the weather forecast and the like and predicts the future active power variation trend of the photovoltaic power supply so as to improve the accuracy of the prediction result and improve the accuracy of the prediction resultThe availability, the prediction system usually adopts the interval prediction method, namely the maximum value and the minimum value which the active power change may reach are predicted;

s4120) if the power prediction system is not deployed, employing a method comprising:

s4121) for the photovoltaic set of power generation, using the current power multiplied by the upper prediction parameter as the future T₁Using the current power multiplied by a lower limit prediction parameter as the lower limit value of the possible fluctuation range of the active power, wherein the upper limit prediction parameter is more than 1 and the lower limit prediction parameter is more than 0;

s4122) for photovoltaic plants not generating electricity, future T using plants with performance consistent or similar to that of the photovoltaic plants (in particular, consistent with the capacity of a single machine)₁The possible fluctuation range of the active power in time is used as the future T of the unit₁The possible fluctuation range of active power in time;

s4123) for the upper limit prediction parameter and the lower limit prediction parameter described in S4121, fixed values may be used, or different parameters may be used at different time points, and the latter is more suitable for photovoltaic power stations with obvious time regularity in the year and day, for example, a higher prediction parameter is used in a period after sunrise, and a lower prediction parameter is used in a period before sunset.

S4130) calculating future T₁The unit active power possible fluctuation range of the photovoltaic power supply unit in time comprises the following steps:

s4131) future T₁Accumulating and summing the upper limits of possible fluctuation ranges of the active power of all the generator sets of the photovoltaic power supply unit within the time, namely the sum is the future T₁The upper limit of the possible fluctuation range of the unit active power of the photovoltaic power supply unit within the time;

s4132) will be T in the future₁Accumulating and summing the lower limits of possible fluctuation ranges of the active power of all the generator sets of the photovoltaic power supply unit in time, namely obtaining the T in the future₁And the lower limit of the possible fluctuation range of the unit active power of the photovoltaic power supply unit in time.

S4200) respectively generating a startup and shutdown sequence for the photovoltaic units, including:

s4210) generating a shutdown sequence of the photovoltaic unit for power generation, wherein the priority is calculated according to the duration of the unit in the power generation state, and the longer the duration of the unit in the power generation state is, the higher the priority is;

s4220) generating a startup sequence of available and unenergy photovoltaic units, the priority being calculated according to the duration of the units in the non-electricity generating state, the longer the duration of the units in the non-electricity generating state, the higher the priority, the so-called available and unenergy units being relative to the unavailable units which cannot be converted into the electricity generating state due to equipment failure or maintenance work.

S4300) respectively generating possible fluctuation range sequences of active power corresponding to the startup and shutdown sequences aiming at the photovoltaic unit, wherein the possible fluctuation range sequences comprise:

s4310) generating an active power possible fluctuation range sequence corresponding to the starting sequence aiming at the photovoltaic unit:

s4311) setting variable u₁，u₁Is 1;

s4312) adding the possible fluctuation range of the active power of the photovoltaic power supply unit obtained in the step S4130 to the sequence u in the starting sequence of the photovoltaic unit₁The possible fluctuation range of the active power of the photovoltaic unit is obtained, and the sequence u in the possible fluctuation range sequence of the active power corresponding to the starting sequence of the photovoltaic unit is obtained₁In which u is ordered₁The upper limit of the range of (a) is equal to the upper limit of the possible fluctuation range of the active power of the photovoltaic power supply unit obtained in the step S4130 and the sequencing u in the starting sequence of the photovoltaic unit₁The upper limit of the possible fluctuation range of the active power of the unit is sorted u₁The lower limit of the range is equal to the sum of the lower limit of the possible fluctuation range of the active power of the photovoltaic power supply unit obtained in the step S4130 and the sequencing u in the starting sequence of the photovoltaic unit₁The lower limit of the possible fluctuation range of the active power of the unit;

s4313) determination of u₁Whether it is equal to the starting sequence length of the photovoltaic unit, if u₁If the length of the starting sequence of the photovoltaic unit is equal to the length of the starting sequence of the photovoltaic unit, the step S4310 is terminated, otherwise u is executed₁＝u₁+1, and then continuing to perform the subsequent steps;

s4314) converting the possible fluctuation range sequence of the active power corresponding to the startup sequence of the photovoltaic unitRank u in column₁Range of-1, plus sequence u in the startup sequence of the photovoltaic set₁The possible fluctuation range of the active power of the photovoltaic unit is obtained, and the sequence u in the possible fluctuation range sequence of the active power corresponding to the starting sequence of the photovoltaic unit is obtained₁In which u is ordered₁Is equal to the rank u₁Upper limit of range of-1 plus sequence u in the photovoltaic set startup sequence₁The upper limit of the possible fluctuation range of the active power of the unit is sorted u₁Is equal to the rank u₁-1 lower limit of range plus the sequence u in the photovoltaic set startup sequence₁The lower limit of the possible fluctuation range of the active power of the unit;

s4315) jumping to step S4313 until u₁Equals to the length of the startup sequence of the photovoltaic set, and ends S4310.

E.g. future T₁The possible fluctuation range of the active power of the photovoltaic power supply unit in time is 310-360 MW, and the starting sequence of the photovoltaic unit is No. 1 machine, No. 3 machine and No. 2 machine]Wherein, the possible fluctuation range of the active power of the

photovoltaic units

1, 2 and 3 is 40-60, 50-70 and 40-80, and the possible fluctuation range sequence of the active power corresponding to the startup sequence of the photovoltaic units is [ (350,420), (390,500) and (440,570)]。

S4320) generating a possible fluctuation range sequence of active power corresponding to the shutdown sequence aiming at the photovoltaic unit, wherein the possible fluctuation range sequence of the active power comprises the following steps:

s4321) setting variable u₂，u₂Is 1;

s4322) subtracting the possible fluctuation range of the active power of the photovoltaic power supply unit obtained in the step S4130 from the sequence u in the photovoltaic shutdown sequence₂The possible fluctuation range of the active power of the photovoltaic unit is obtained, and the order u in the possible fluctuation range sequence of the active power corresponding to the photovoltaic shutdown sequence is obtained₂In which u is ordered₂Is equal to the upper limit of the possible fluctuation range of the active power of the photovoltaic power supply unit obtained in S4130 minus the sequence u in the photovoltaic shutdown sequence₂The upper limit of the possible fluctuation range of the active power of the photovoltaic unit is sorted u₂The lower limit of the range of (1) is equal to the possible active power wave of the photovoltaic power supply unit obtained in S4130Sequence u in dynamic range lower limit minus photovoltaic shutdown sequence₂The lower limit of the possible fluctuation range of the active power of the photovoltaic unit;

s4323) judgment of u₂Whether it is equal to the photovoltaic shutdown sequence length, if u₂Equal to the photovoltaic shutdown sequence length, step S4320 is terminated, otherwise u is executed₂＝u₂+1, and then continuing to perform the subsequent steps;

s4324) sorting u in the possible fluctuation range sequence of the active power corresponding to the photovoltaic shutdown sequence₂Range of-1, minus the order u in the photovoltaic shutdown sequence₂The possible fluctuation range of the active power of the photovoltaic unit is obtained, and the order u in the possible fluctuation range sequence of the active power corresponding to the photovoltaic shutdown sequence is obtained₂In which u is ordered₂Is equal to the rank u₂Upper limit of range of-1 minus the order u in the photovoltaic shutdown sequence₂The upper limit of the possible fluctuation range of the active power of the photovoltaic unit is sorted u₂Is equal to the rank u₂Lower bound of range of-1 minus the order u in the photovoltaic shutdown sequence₂The lower limit of the possible fluctuation range of the active power of the photovoltaic unit;

s4325) jumping to step S4323 until u₂Equals the photovoltaic shutdown sequence length and ends step S4320.

S4400) calculating the real unit active power value sending parameter of the photovoltaic power supply unit, including:

s4410) initially setting the active power real-emission value parameter calculation quantity of the photovoltaic power supply unit to be equal to the unit active power real-emission value;

s4420) accumulating the output dead zones of each set of the photovoltaic power supply unit which is given by scheduling or manually set to obtain the unit output dead zone of the photovoltaic power supply unit;

s4430) comparing the real active power value of the photovoltaic power supply unit with the calculated quantity and the real active power value of the photovoltaic power supply unit at the current period according to a fixed period, wherein the comparing comprises the following steps:

s4431) if the absolute value of the difference value of the two is less than or equal to the output dead zone of the photovoltaic power supply unit, the real output parameter of the active power of the photovoltaic power supply unit and the calculated quantity are kept unchanged;

s4432) if the absolute value of the difference value of the two is larger than the output dead zone of the photovoltaic power supply unit, the active power real-time value parameter calculation amount of the photovoltaic power supply unit is equal to the active power real-time value of the photovoltaic power supply unit in the current period.

For example, the dead zone of the output of the photovoltaic power supply unit is 20MW, the real active power value of the photovoltaic power supply unit is 300MW both in the calculated quantity and the real active power value of the unit, the real power value of the unit is changed into 305MW due to power fluctuation, the absolute value of the difference value between the real power value of the photovoltaic power supply unit and the calculated quantity 300MW and the real power value of the unit 305MW is 5MW which is smaller than the dead output zone 20MW, therefore, the real active power parameter of the photovoltaic power supply unit keeps 300MW unchanged, and later, due to further power fluctuation, the real unit active power parameter changes to 321MW, so that the absolute value of the difference value between the real active power parameter of the photovoltaic power supply unit 300MW and the real unit active power parameter 321MW changes to 21MW, which is larger than the output dead zone 20MW, therefore, the active power real-emitting value of the photovoltaic power supply unit is changed into 321MW according to the unit active power real-emitting value.

S4500) calculating a unit active power real-time value of the photovoltaic power supply unit and a calculated value filtering value, including:

s4510) initially setting the active power real-emission value of the photovoltaic power supply unit and the calculated amount filter value to be equal to the unit active power real-emission value;

s4520) calculating a filtering threshold of an active power real-time value of the photovoltaic power supply unit, including:

s4521) setting a scaling coefficient lambda, wherein lambda is larger than 1;

s4522) a filtering threshold of the active power real output value of the photovoltaic power supply unit is equal to the unit output dead zone multiplied by λ in S4420, and in this embodiment, assuming that λ is 3, the filtering threshold is equal to 3 times of the unit output dead zone.

S4530) the real active power value of the photovoltaic power supply unit is compared with the calculated value of the filter and the real active power value of the photovoltaic power supply unit in the current period according to a fixed period, and the method comprises the following steps:

s4531) if the absolute value of the difference value between the active power parameter and the calculated value is less than or equal to the filtering threshold obtained in S4522, keeping the active power actual value parameter and the calculated value filtering value of the photovoltaic power supply unit unchanged;

s4532) if the absolute value of the difference value between the real power value and the real power value is larger than the filtering threshold obtained in S4522, the real power value of the photovoltaic power supply unit and the calculated value are equal to the real power value of the real power of the photovoltaic power supply unit in the current period.

S4600) calculating a unit primary frequency modulation target adjustment amount of the photovoltaic power supply unit, including:

s4610) calculating a power grid frequency deviation, wherein the power grid frequency deviation is equal to the power grid rated frequency (50Hz) minus the real-time frequency of the power grid;

s4620) if the absolute value of the power grid frequency deviation is smaller than or equal to a primary frequency modulation threshold (given by scheduling), the unit primary frequency modulation target regulating quantity of the photovoltaic power supply unit is equal to 0;

s4630) if the absolute value of the power grid frequency deviation is larger than a primary frequency modulation threshold, the primary frequency modulation target regulating quantity of the unit of the photovoltaic power supply unit is equal to the real power output value of the photovoltaic power supply unit multiplied by the power grid frequency deviation and then multiplied by a photovoltaic primary frequency modulation regulating coefficient (power grid given parameter).

S5000) the operation of the complementary integration unit is given below, which specifically includes:

the method comprises the following steps of distributing unit active power target values of a thermal power supply unit and an energy storage power supply unit, setting a primary frequency modulation adjusting coefficient of the thermal power supply unit, calculating a start-stop operation suggestion of a photovoltaic power supply unit to meet the adjusting requirements of a total active power set value and primary frequency modulation of a multi-energy complementary integrated power supply and the charging and discharging requirements of an energy storage power supply battery, wherein a control model is shown in figure 1, and in order to visually display an adjusting effect, the influence of the primary frequency modulation is eliminated in the control model, but technicians in the industry can easily know that the implementation effect of the method cannot be influenced even the primary frequency modulation is introduced, and the method comprises the following steps:

s5100) calculating charging and discharging correction power of the energy storage power supply unit;

s5110) manually setting charge-discharge parameter alpha₁And an emergency charge-discharge parameter alpha₂Wherein 0 < alpha₁＜α₂，α₁And alpha₂The units of (a) and (h) are all/h, and in practical engineering, the energy storage unit battery is generally configured according to the supporting rated power for charging or discharging for 30 minutes, so that the embodiment can set alpha₁And alpha₂The power is respectively 0.6/h and 1.2/h, namely, the battery is charged and discharged according to 30 percent and 60 percent of rated power respectively;

s5120) calculating the charge-discharge coefficient alpha every fixed period according to the battery total electric quantity state of the energy storage power supply unit calculated in the S3300, wherein the calculation comprises the following steps:

s5121) when the total amount of the battery is in an extremely ideal state of charge, the charge-discharge coefficient α is 0;

s5122) when the total amount of the battery is in a low state of charge, the charge-discharge coefficient α is α₁；

S5123) when the total amount of the battery is in the extremely low state of charge, the charge-discharge coefficient α is α₂；

S5124) when the total amount of the battery is in a higher state of charge, the charge-discharge coefficient α is- α₁；

S5125) when the total amount of the battery is in an extremely high state of charge, the charge-discharge coefficient α is- α₂；

S5126) when the total battery amount is in the more ideal state, the charge/discharge coefficient remains unchanged, the more ideal state of the total battery amount is made to be a buffer area for changing the charge/discharge state to prevent the charge/discharge correction power from changing frequently, i.e. the charge/discharge coefficient in the more ideal state is determined by the previous total state of charge of the battery, when the total battery amount is changed from the extremely ideal state of charge to the more ideal state of charge, the charge/discharge coefficient α is 0, and when the total battery amount is changed from the lower state of charge to the more ideal state of charge, the charge/discharge coefficient α is α₁When the total quantity of the battery is changed from a higher state of charge to a more ideal state of charge, the charge-discharge coefficient alpha is-alpha₁；

S5130) calculating the charging and discharging correction power of the energy storage power supply unit according to the charging and discharging coefficient obtained in the S5120, wherein the charging and discharging correction power is equal to

S5200) calculating a unit active power target value of the thermal power supply by the complementary integrated unit, wherein the unit active power target value of the thermal power supply is equal to a calculated quantity filter value obtained by subtracting the unit active power actual value of the photovoltaic power supply unit obtained by the S4500 from the total active power set value of the multifunctional complementary integrated power supply, and adding the charging and discharging correction power obtained by the S5100; in the step, the real unit active power value of the photovoltaic power supply unit is used for participating in the calculated value filtering value, and the energy storage power supply unit is considered to be included in the multi-energy complementary integrated power supply, so that the sensitivity of the firepower power supply unit to the random fluctuation of the real active power value of the photovoltaic power supply unit is reduced appropriately under the condition, and according to the S4500 example, the sensitivity of the firepower supply unit to the random fluctuation of the real active power value of the photovoltaic power supply unit is one third of that of the energy storage power supply unit;

s5210) comparing the unit active power target value of the thermal power supply with the unit combined operation zone described in S2250, there are two possible results: when the unit active power target value is included in the unit combined operation area, the unit active power target value is feasible, and then according to the S2000 method, unit-level AGC distribution is performed on the unit active power target value of the thermal power supply obtained in S5200, and primary frequency modulation and secondary frequency modulation are performed according to the active power of the thermal power supply unit, including:

s5211) calculating a primary frequency modulation adjustment coefficient of the thermal power supply unit, wherein each unit of the thermal power supply unit is executed according to the primary frequency modulation adjustment coefficient when actually executing a primary frequency modulation task;

the complementary integration unit calculates a primary frequency modulation scaling coefficient of the thermal power supply unit, and the primary frequency modulation scaling coefficient is equal to (rated capacity of active power of the photovoltaic power supply unit + rated capacity of active power of the thermal power supply unit) ÷ rated capacity of active power of the thermal power supply unit; assuming that the active power rated capacity of the thermal power supply unit is 200MW and the active power rated capacity of the photovoltaic power supply unit is 100MW, the primary frequency modulation scaling coefficient of the thermal power supply unit is (200+100)/200 ═ 1.5;

the primary frequency modulation adjustment coefficient of the thermal power supply unit is equal to the primary frequency modulation adjustment coefficient of the thermal power supply unit issued by the power grid multiplied by a primary frequency modulation scaling coefficient;

when each unit of the thermal power supply unit actually executes primary frequency modulation regulation, the adjustment is carried out according to the scaled primary frequency modulation regulation coefficient; assuming that the primary frequency modulation adjustment amount of a unit of the original thermal power supply unit is 40MW when a certain deviation occurs in the grid frequency, the primary frequency modulation adjustment amount of the unit is amplified to 40 × 1.5 — 60MW in order to undertake the primary frequency modulation task of the photovoltaic power supply.

S5220) when the unit active power target value is not included in the unit combined operation area, the unit active power target value is not feasible, and then it is necessary to find an operation proposal for making the unit active power target value feasible through the following steps:

s5221) finding an operation advice for the fire power supply unit, including:

1) according to the S2320 method, searching for operation suggestions for enabling unit active power target values of the thermal power supply by putting the units which are not put into AGC control, and sequencing the priority of the operation suggestions;

2) according to the method of S2330, searching operation suggestions for enabling the unit active power target value of the thermal power supply by converting the unit which does not generate power into a power generation state and putting AGC into operation, and sequencing the priority of the operation suggestions;

3) according to the S2340 method, a running operation advice that makes the unit active power target value of the thermal power supply feasible by turning the generating unit to the non-generating state is found, and the priority of the running operation advice is ranked.

S5300) generating a photovoltaic unit start-up and shutdown operation suggestion and assigning:

s5310) calculating future T₁Unit active power accommodation range of photovoltaic power supply unit in time, wherein T₁Setting parameters for the human described in S4100:

s5311) calculating future T₁The photovoltaic power supply unit's of each time point in time unit active power accommodation range lower limit or each continuous interval lower limit of accommodation range includes:

if the scheduling issues the multi-energy complementary integrated power supply in advancePower plan curve, then future T₁Subtracting the upper limit of the combined operation area of the thermal power supply unit obtained by S2250 from the total active power set value of the multi-energy complementary integrated power supply at each time point in time is the upper limit of the combined operation area of the thermal power supply unit in the future T₁The lower limit of the unit active power accommodation range of the photovoltaic power supply unit at each time point in time;

if the active power plan curve of the multi-energy complementary integrated power supply is not issued in advance in the scheduling process, the upper limit of the combined operation area of the thermal power supply unit obtained by subtracting the S2250 from the total active power set value of the current multi-energy complementary integrated power supply is the future T₁The lower limit of the unit active power accommodation range of the photovoltaic power supply unit at each time point in time;

s5312) calculating future T₁The upper limit of the unit active power accommodation range of the photovoltaic power supply unit at each time point in time comprises:

if the active power plan curve of the multi-energy complementary integrated power supply is issued in advance by scheduling, the future T is determined₁Subtracting the lower limit of the combined operation area of the thermal power supply unit from the total active power set value of the multi-energy complementary integrated power supply at each time point in time by S2250 to obtain the lower limit of the combined operation area of the thermal power supply unit, which is the future T₁The upper limit of the unit active power accommodation range of the photovoltaic power supply unit at each time point in time;

if the active power plan curve of the multi-energy complementary integrated power supply is not issued in advance in the scheduling process, the lower limit of the combined operation area of the thermal power supply unit obtained by subtracting the S2250 from the total active power set value of the current multi-energy complementary integrated power supply is the future T₁The upper limit of the unit active power accommodation range of the photovoltaic power supply unit at each time point in time;

s5313) future T₁The unit active power accommodation range of the photovoltaic power supply unit in time is T in the future₁The unit active power accommodation ranges of the photovoltaic power supply units at each time point in time are intersected, namely T in the future₁The upper limit of the unit active power accommodation range of the photovoltaic power supply unit in time is equal to the minimum value of the upper limit of the accommodation range at each time point, T in the future₁The lower limit of the unit active power accommodation range of the photovoltaic power supply unit in time is equal to the maximum value of the lower limit of the accommodation range at each time pointSuppose future T₁The corresponding accommodation ranges of each time point in time are (150,300), (120,270), (100,250), (160,310) and (200,350), and the intersection is obtained to obtain the future T₁The unit active power accommodation range of the photovoltaic power supply unit in time is (200, 250).

S5320) calculating the on-off state and the future T of the current photovoltaic power supply unit₁The mismatch quantization value of the total active power set value of the multi-energy complementary integrated power supply in time is as follows:

s5321) calculating future T of S5310₁Unit active power accommodation range and future T of photovoltaic power supply unit in time₁The upper limit mismatching degree of the possible fluctuation range of the active power of the photovoltaic power supply unit in time is T in the future₁Subtracting future T from the upper limit of the possible fluctuation range of the active power of the photovoltaic power supply unit in time₁Judging the upper limit of the unit active power accommodation range of the photovoltaic power supply unit within time, if the upper limit active power accommodation range is larger than 0, determining that the upper limit mismatching degree is equal to the calculation result, and otherwise, determining that the upper limit mismatching degree is equal to 0;

s5322) calculating future T of S5310₁The unit active power accommodation range of the photovoltaic power supply unit in time includes and is T in future₁The lower limit mismatching degree of the possible fluctuation range of the active power of the photovoltaic power supply unit in time is T in the future₁Subtracting the future T from the lower limit of the unit active power accommodation range of the photovoltaic power supply unit in time₁Judging the calculation result according to the lower limit of the possible fluctuation range of the active power of the photovoltaic power supply unit in time, wherein if the lower limit is larger than 0, the lower limit mismatching degree is equal to the calculation result, and otherwise, the lower limit mismatching degree is equal to 0;

s5323) subtracting the lower limit mismatching degree obtained by S5322 from the upper limit mismatching degree obtained by S5321, and taking the absolute value of the result to obtain the on-off state and the future T of the current photovoltaic power supply unit set₁The mismatch degree quantization value of the total active power set value of the multi-energy complementary integrated power supply in time; e.g. future T₁The unit active power accommodation range of the photovoltaic power supply unit in time is (200,250), and when the possible fluctuation range of the photovoltaic power supply unit active power is (100,130), the upper limit of the range is not matchedThe distribution degree is max [0, 130-]0, the mismatch at the lower end of the range is max [0,200-]100, and then the mismatch quantization value | 0-100 | ═ 100, where max [, [ m ] is]Is a function of the maximum.

S5330) finding an operation proposal for stopping the photovoltaic unit for generating power, which specifically comprises the following steps:

s53310) manually setting a judgment threshold parameter for the suggested shutdown operation;

s53320) setting variable v₃，v₃Is 1;

s53330) if v₃If the length of the photovoltaic shutdown sequence is smaller than or equal to the length of the photovoltaic shutdown sequence, setting an original mismatching degree quantization value variable, wherein the original mismatching degree quantization value variable is equal to the mismatching degree quantization value obtained in the step S5320, otherwise, skipping to the step S53360;

s53340) calculating the sequence v in the possible active power fluctuation range sequence corresponding to the photovoltaic shutdown sequence₃Range and future T of₁The mismatch degree quantization value of the total active power set value of the multi-energy complementary integrated power supply in time comprises the following steps:

s53341) calculating future T of S5310₁Sequencing v in the possible fluctuation range sequence of the active power corresponding to the photovoltaic shutdown sequence and the unit active power accommodation range of the photovoltaic power supply unit in time₃The upper limit mismatching degree of the range, and sorting v in the possible fluctuation range sequence of the active power corresponding to the photovoltaic shutdown sequence₃Upper range limit minus the future T₁Judging the calculation result according to the upper limit of the active power accommodation range of the photovoltaic power supply unit in time, wherein if the upper limit is larger than 0, the upper limit mismatching degree is equal to the calculation result, and otherwise, the upper limit mismatching degree is equal to 0;

s53342) calculating future T₁Sequencing v in the possible fluctuation range sequence of the active power corresponding to the photovoltaic shutdown sequence and the unit active power accommodation range of the photovoltaic power supply unit in time₃The lower limit mismatch of the range of (1), will be T in the future₁Sequencing v in the sequence of subtracting the possible fluctuation range of the active power corresponding to the photovoltaic shutdown sequence from the lower limit of the active power accommodation range of the unit of the photovoltaic power supply unit in time₃And to the calculation resultJudging, if the lower limit mismatching degree is larger than 0, the lower limit mismatching degree is equal to the calculation result, otherwise, the lower limit mismatching degree is equal to 0;

s53343) subtracting the lower limit mismatching degree obtained by S53342 from the upper limit mismatching degree obtained by S53341, and taking an absolute value of the result to obtain a sequence v in the possible fluctuation range sequence of the active power corresponding to the photovoltaic shutdown sequence₃Range and future T of₁And the mismatching degree quantization value of the total active power set value of the multi-energy complementary integrated power supply in time.

S53350) subtracting the quantization value of mismatch degree obtained in S53343 from the original quantization value variable of mismatch degree, and performing the following operations according to the calculation result, including:

s53351) if the calculation result is equal to or greater than the judgment threshold parameter set in S53310, v is₃＝v₃+1 if v is present at this time₃If the length of the photovoltaic shutdown sequence is larger than the length of the photovoltaic shutdown sequence, jumping to a step S53360, otherwise, updating the original mismatch quantization value variable into the mismatch quantization value obtained in the step S53343, and jumping to the step S53340 to continue execution;

s53352) if the calculation result is less than the judgment threshold parameter set in S53310, jumping to S53360 and continuing execution.

S53360) according to the variable v₃Generates an operation recommendation, comprising:

s53361) if v₃If 1, no operation suggestion is generated;

s53362) if v₃If the number of the photovoltaic shutdown sequences is more than 1, a shutdown operation suggestion is generated, and the suggestions are used for sequencing 1 to v in the photovoltaic shutdown sequence₃-1 the corresponding photovoltaic unit executes shutdown operation;

s5340) finding an operation recommendation for starting up a photovoltaic power plant which is available and not generating power, comprising:

s53410) manually setting a judgment threshold parameter of the suggested startup operation;

s53420) setting variable v₄，v₄Is 1;

s53430) if v₄If the length of the starting sequence of the photovoltaic unit is less than or equal to the length of the starting sequence of the photovoltaic unit, setting an original mismatching degree quantization value variable which is equal to that obtained in S5320Otherwise, jumping to step S53460;

s53440) calculating the sequence v in the possible fluctuation range sequence of the active power corresponding to the startup sequence of the photovoltaic unit₄Range and future T of₁The mismatch degree quantization value of the total active power set value of the multi-energy complementary integrated power supply in time comprises the following steps:

s53441) calculating future T of S5310₁Sequencing v in the possible fluctuation range sequence of the active power corresponding to the starting sequence of the photovoltaic unit and the unit active power accommodation range of the photovoltaic power supply unit in time₄The upper limit mismatching degree of the range, and sorting v in the possible fluctuation range sequence of the active power corresponding to the starting sequence of the photovoltaic unit₄Upper range limit minus the future T₁Judging the calculation result according to the upper limit of the active power accommodation range of the photovoltaic power supply unit in time, wherein if the upper limit is larger than 0, the upper limit mismatching degree is equal to the calculation result, and otherwise, the upper limit mismatching degree is equal to 0;

s53442) calculating the future T of S5310₁Sequencing v in the possible fluctuation range sequence of the active power corresponding to the starting sequence of the photovoltaic unit and the unit active power accommodation range of the photovoltaic power supply unit in time₄The lower limit mismatch of the range of (1), will be T in the future₁In time, the lower limit of the unit active power accommodation range of the photovoltaic power supply unit is subtracted by the sequence v in the possible active power fluctuation range sequence corresponding to the starting sequence of the photovoltaic unit₄If the lower limit is larger than 0, the lower limit mismatching degree is equal to the calculation result, otherwise, the lower limit mismatching degree is equal to 0;

s53443) subtracting the lower limit mismatching degree obtained by the S53442 from the upper limit mismatching degree obtained by the S53441, and taking an absolute value of the result to obtain a sequence v in a possible fluctuation range sequence of active power corresponding to the startup sequence of the photovoltaic unit₄Range and future T of₁And the mismatching degree quantization value of the total active power set value of the multi-energy complementary integrated power supply in time.

S53450) subtracting the quantization value of mismatch degree obtained in S53443 from the original quantization value variable of mismatch degree, and performing the following operations according to the calculation result, including:

s53451) if the calculation result is equal to or greater than the judgment threshold parameter set in S53410, v is₄＝v₄+1 if v is present at this time₄If the length of the starting sequence of the photovoltaic unit is larger than the length of the starting sequence of the photovoltaic unit, jumping to a step S53460, otherwise, updating the original mismatching degree quantization value variable into the mismatching degree quantization value obtained in the step S53443, and jumping to the step S53440 to continue execution;

s53452) if the calculation result is less than the judgment threshold parameter set in S53410, jumping to step S53460 and continuing execution.

S53460) according to the variable v₄Generates an operation recommendation, comprising:

s53461) if v₄If 1, no operation suggestion is generated;

s53462) if v₄If the number of the photovoltaic units is more than 1, generating a starting operation suggestion, and sequencing 1 to v in a starting sequence of the photovoltaic units according to the suggestion₄The photovoltaic unit corresponding to the-1 executes the starting operation.

Then, orderly displaying the shutdown operation suggestions of the photovoltaic units generated in the step S5330 respectively, and assigning the shutdown operation suggestions to the photovoltaic power supply units;

and respectively and orderly displaying the power-on operation suggestions of the photovoltaic units generated in the step S5340 and assigning the power-on operation suggestions to the photovoltaic power supply units.

S5400) calculating a unit active power target value of the energy storage power supply unit, including:

s5410) adding a unit primary frequency modulation correction quantity of the fire power supply unit to a total active power set value of the multi-energy complementary integrated power supply, then subtracting a unit active power actual emission value of the photovoltaic power supply unit from a calculated quantity, and then subtracting a unit active power actual emission value of the fire power supply unit to obtain an active power total adjustment deviation of the fire power supply unit and the photovoltaic power supply unit;

s5420) initially setting the compensation adjustment quantity of the energy storage power supply unit to be the active power total adjustment deviation obtained in S5410, and then comparing the compensation adjustment quantity of the energy storage power supply unit with the current active power total adjustment deviation according to a fixed period:

s5421) when the absolute value of the difference value of the two is larger than the unit active power regulation dead zone of the energy storage power supply unit, the single compensation regulation quantity of the energy storage power supply unit is equal to the total active power regulation deviation of the firepower power supply unit and the photovoltaic power supply unit in the current period;

s5422) when the absolute value of the difference value of the two is less than or equal to the unit active power regulation dead zone of the energy storage power supply unit, the compensation regulation quantity of the energy storage power supply unit is kept unchanged.

S5430) performing dead zone processing on the compensation adjustment quantity of the energy storage power supply unit with reference to the method of S5430 to obtain a unit active power target value of the energy storage power supply unit, where the processing logic is as shown in fig. 5, and includes:

s5431) manually setting the timer and the time parameter T₄；

S5432) when the absolute value of the total active power regulation deviation obtained in the S5410 is smaller than or equal to the active power regulation dead zone of the thermal power supply unit, starting timing by a timer set in the S5431;

s5433) resetting and clearing a timer set in S5431 when the absolute value of the total active power regulation deviation obtained in S5410 is larger than the active power regulation dead zone of the thermal power supply unit;

s5434) when the timer time is less than the time parameter T₄When the active power target value of the energy storage power supply unit is equal to the compensation adjustment quantity obtained in S5420;

s5435) when the timer time is more than or equal to the time parameter T₄And when the active power target value of the energy storage power supply unit is equal to 0.

And S5500) the energy storage power supply unit performs unit-level AGC distribution on the unit active power target value obtained in the S5430, and adjusts the active power of each energy storage unit.

The regulation effect of the multi-energy complementary integrated power supply in the control model shown in fig. 1 is shown in fig. 6, and it can be easily seen that: 1) the compensation function of the thermal power supply for the substantial deviation of the real active power value of the photovoltaic power supply unit and the compensation function of the energy storage power supply for the random fluctuation of the real active power value of the photovoltaic power supply unit are benefited, the total real active power value of the multi-energy complementary integrated power supply always keeps extremely high stability, and the energy storage power supply is basically in a low load state except for individual time periods, so that the requirement of 'shallow charging and shallow discharging' of a battery is met;

2) due to the excellent adjusting performance of the energy storage power supply, when the total active power set value of the multi-energy complementary integrated power supply is changed from 300MW to 400MW within 70s, the real active power value responsiveness of the multi-energy complementary integrated power supply is very good, and indexes such as adjusting delay, adjusting speed and adjusting precision are all at a high level;

3) when the charging and discharging correction power of the energy storage power supply battery is changed from 0 to 100M in 140s, the real active power value of the thermal power supply unit is increased, so that the energy storage power supply battery can enter a charging state, and the process does not cause adverse effects on the stability of the real total active power value of the multi-energy complementary integrated power supply.

In order to further show the characteristics of shallow charging and shallow discharging of the energy storage power supply battery in the method of the invention, a simulation model of a thermal power supply, an energy storage power supply and a photovoltaic is further constructed, wherein 3 energy storage units are arranged in an energy storage power supply unit, the battery capacity ratio of the 3 energy storage units is 5:8:10, a control model is shown in fig. 7, wherein the relationship graphs of the actual value of the total active power of the integrated power supply, the actual value of the active power of the thermal power supply unit, the actual value of the active power of the photovoltaic unit, the actual value of the active power of each unit of the energy storage power supply unit, the battery charge state of each unit of the energy storage power supply unit, the battery charge capacity ratio of each unit of the energy storage power supply unit, the total battery charge capacity ratio of the energy storage power supply unit, the charging and discharging correction power of the energy storage power supply unit and the like are respectively shown in fig. 8, and the regulation effect of fig. 8 shows that:

1) the adjusting amplitude of the energy storage unit during active power adjustment is related to the battery capacity and the battery charge state, although the battery capacity of the energy storage unit 3 is twice that of the energy storage unit 1, the discharging amplitude of the energy storage unit 3 is smaller than that of the energy storage unit 1 on the contrary because the initial charge capacity proportion of the battery of the energy storage unit 3 is far lower than that of the energy storage unit 1;

2) although there is a large difference in the battery charge capacity ratios of the 3 energy storage units artificially set in the initial stage of simulation, under the control of the "shallow charging and shallow discharging" strategy of the present invention, the capacity ratios of the battery charges of all the energy storage units gradually tend to be consistent, and meanwhile, as described above, the charging and discharging strategy of the present invention can maintain the total charge capacity of the unit cells of the energy storage units in a better balance, so the batteries of all the energy storage units are naturally in a more balanced state (neither overcharging nor overdischarging).

The embodiments given above are preferable examples for implementing the present invention, and the present invention is not limited to the above-described embodiments. Any non-essential addition and replacement made by the technical characteristics of the technical scheme of the invention by a person skilled in the art belong to the protection scope of the invention.

Claims

1. A multi-energy power supply control method giving consideration to both battery state and active power dynamic stability is characterized in that a thermal power supply, an energy storage power supply and a photovoltaic energy source are coordinately controlled through a multi-energy complementary integrated power supply centralized control center:

2. The method according to claim 1, wherein the parameters obtained by the complementary integrated unit include:

s1100), parameters input by a complementary integration unit:

s1113) the real active power value of the unit is respectively equal to the sum of the real active power values of the units of the thermal power supply, the energy storage power supply and the photovoltaic power supply unit;

s1114) unit active power regulation dead zones are respectively equal to the sum of single-machine active power regulation dead zones of the units in which the thermal power supply and the energy storage power supply unit are running;

s1120) input parameters transmitted by the thermal power supply unit:

s1122) a unit joint operation area of the thermal power supply unit;

s1130) parameters sent by the energy storage power supply unit: the charging and discharging correction power of the energy storage power supply unit is obtained by the energy storage power supply unit according to the battery charge state and the charging and discharging coefficient of each energy storage unit;

s1140) input parameters sent by the photovoltaic power supply unit:

s1142) the real unit active power value of the photovoltaic power supply unit is involved in the calculated value of the filtered value, and the photovoltaic power supply unit updates according to the real unit active power value, the scaling coefficient and the dead output area of each photovoltaic unit according to a fixed period;

s1145) the unit primary frequency modulation target regulating quantity of the photovoltaic power supply unit is equal to the sum of the single-machine primary frequency modulation target regulating quantity of the photovoltaic unit generating power.

3. The method of claim 1, wherein the operation of the thermal power supply unit comprises:

s2100) determining the unit type of the thermal power supply unit: according to the difference of the active power regulation controlled states of the generator sets, dividing a single-machine open-loop unit, a single-machine closed-loop unit, a unit which is put into AGC (automatic gain control) and a unit which is not put into AGC;

s2200) establishing a combined output model for each unit of AGC, calculating a combined operation area and a combined recommended operation area, and determining the current single-machine AGC active power distribution value of each unit;

s2300) comparing a unit active power target value of the thermal power supply with a unit combined operation area, wherein when the unit active power target value is included in the unit combined operation area, the unit active power target value is feasible; when the unit active power target value is not included in the unit combined operation area, the unit active power target value is not feasible, and an operation suggestion enabling the unit active power target value to be feasible is searched; classifying the generated running operation suggestions, and displaying the running operation suggestions in order according to the obtained priority;

s2400) calculating a single AGC active power distribution value which is put into an AGC unit: calculating unit AGC active power distribution values of the thermal power supply, and distributing AGC active power to each unit which is put into the AGC;

s2500) adjusting the active power of each single-machine closed-loop unit of the thermal power supply unit, determining a single-machine active power set value of each single-machine closed-loop unit, calculating the deviation between a single-machine active power actual output value and a single-machine active power set value by taking the single-machine active power set value as a target, and outputting a continuous signal according to a calculation result to adjust the single-machine active power actual output value of the unit so that the single-machine active power actual output value of the unit tends to the single-machine active power set value and is finally stabilized in the adjustment dead zone range of the single-machine active power set value.

4. The multi-energy power supply control method considering both battery status and active power dynamic stability as claimed in claim 3, wherein the steps of the operation of the thermal power supply unit specifically include the following operations:

s2200) establishing a combined output model for each unit of the AGC, and calculating a combined operation area and a combined recommended operation area:

s2210) determining a single machine recommended operation area and a single machine forbidden operation area which are put into each unit of the AGC:

s2211) a single machine operation forbidding area refers to a load area for forbidding the single machine active power set value of the unit to be set in the load area; the real value of the single-machine active power of the unit is allowed to pass through or pass through the single-machine forbidden operation area, but is not allowed to reside or stay in the single-machine forbidden operation area for a long time;

s2212) a single machine suggested operation area is a load area with high unit operation efficiency and stable operation when the real single machine active power value of the unit is in the middle; under the condition that the conditions allow, the single machine active power set value of the unit is preferably set in a single machine suggested operation area;

s2215) the low-load area of the thermal power generating unit is a single-machine operation forbidden area, the single-machine operation forbidden area is about 0-50% of rated capacity, and the rest part of the rated capacity after deducting the single-machine operation forbidden area is a single-machine operation recommended area;

s2220) establishing a suggested combined output model of the unit which is put into the AGC, and calculating a combined suggested operation area which is put into the AGC unit: adding the single machine suggested operation areas put into the AGC unit to obtain a combined suggested operation area put into the AGC unit;

s2243) for the stand-alone open-loop unit which is not put into the AGC, the stand-alone AGC active power distribution value tracks the stand-alone active power set value, and the stand-alone active power set value is assigned by the stand-alone active power real sending value, namely when the stand-alone active power set value is not equal to the stand-alone active power real sending value and the absolute value of the difference between the stand-alone active power set value and the stand-alone active power real sending value is larger than the stand-alone active power regulation dead zone, the stand-alone active power real sending value is written into the stand-alone active power set value;

s2250) adding the joint recommended operation area obtained in S2220 into the single AGC active power distribution value of the AGC unit and obtaining a unit joint recommended operation area of the thermal power supply and a unit joint operation area of the thermal power supply;

s2321) setting a loop variable i₁，i₁Is set to 1;

s2324) C (j) listed respectively as S2323₁,i₁) A combination mode is selected from various modesThe unit which is put into the AGC is assumed to be put into the AGC, a unit joint operation area and a unit joint suggested operation area are calculated by adopting the S2200 method, and then the feasibility of the unit active power target value is judged again by adopting the S2300 method according to the newly calculated unit joint operation area;

s2325) according to the calculation result of S2324, if the unit active power target value is feasible by the unit joint operation zone regenerated in only 1 mode, generating an operation proposal; if the unit active power target value can be enabled by the unit joint operation zone regenerated in multiple ways, respectively generating operation suggestions according to the ways, and jumping to the step S2326 to continue executing; if the unit active power target value is feasible without the unit joint operation zone regenerated in any way, i₁＝i₁+1, then go to step S2322 for i₁Judging whether the number of the units not put into the AGC is larger than that of the units not put into the AGC, and determining whether to execute the subsequent steps or not according to the judgment result;

s2326) carrying out priority ordering on the multiple operation suggestions generated in the S2325, wherein the ordering is based on the following steps: the difference absolute value of the unit active power target value and the unit combined operation area boundary is large or small;

s2331) setting a circulation variable i₂，i₂Is set to 1;

s2335) generating an operation proposal according to the calculation result of S2334 if the unit active power target value is feasible by the unit joint operation area regenerated in 1 way or only by the unit joint operation area regenerated in 1 way; if the unit combined operation area regenerated by multiple modes can enable the unit active power target value to be feasible, respectively generating operation suggestions according to the modes, namely converting the available and non-generating units selected by the corresponding modes into a generating state and putting the generating state into AGC (automatic gain control), and jumping to the step S2336 to continue to execute; if the unit active power target value is feasible without the unit joint operation zone regenerated in any way, i₂＝i₂+1, and then go to step S2332 for i₂Judging whether the number of the units is larger than the number of the available and non-power-generating units, and determining whether to execute the subsequent steps according to the judgment result;

s2336) carrying out priority ordering on the multiple operation suggestions generated in the S2335 according to the absolute value of the difference value between the unit active power target value and the unit combined operation zone boundary;

s2341) setting a Loop variable i₃，i₃Is set to 1;

s2343) listing and selecting i from all power generation units₃Of a tableAll combinations, C (j)₃,i₃) Wherein j is₃The number of generating units;

s2345) generating an operation proposal according to the calculation result of S2344 if the unit active power target value is feasible by the unit joint operation area regenerated in 1 way or only regenerated in 1 way; if the unit combined operation area regenerated in multiple modes can enable the unit active power target value to be feasible, respectively generating operation suggestions for converting the generating set selected in the corresponding mode into a non-generating state according to the modes, and jumping to the step S2346 to continue execution; if the unit active power target value is feasible without the unit joint operation zone regenerated in any way, i₃＝i₃+1, and then go to step S2342 for i₃Judging whether the number of the units is larger than the number of the generating sets or not, and determining whether to execute the subsequent steps or not according to the judgment result;

s2346) carrying out priority ranking on the multiple operation suggestions generated in the S2345 according to the fact that the operation suggestions are selected from the generating set i correspondingly₃The combination mode of the station set and the range of the unit joint operation area and the unit joint recommended operation area after the operation suggestions obtained in step S2344 are changed correspondingly, and the ranking is respectively from high to low according to the importance degree: selecting the number of the units which are not put into AGC and the units which are put into AGC in the units, and the absolute value of the difference value between the unit active power target value and the unit combined operation zone boundary;

s2350) classifying the operation suggestions generated by the S2320, the S2330 and the S2340, and displaying the operation suggestions in order according to the priority;

s2400) calculating a single AGC active power distribution value which is put into an AGC unit:

s2411) calculating the distribution values of the active power of all single AGC units which are not put into the AGC unit;

s2412) subtracting all single AGC active power distribution values which are not put into the AGC unit from the unit active power target value to obtain a unit AGC active power distribution value;

s2510) determining the single-machine active power set value of each single-machine closed-loop unit:

5. The method of claim 1, wherein the operation of the energy storage power supply unit comprises:

s3100) calculating the charged capacity proportion of each energy storage unit of the energy storage power supply unit_riAnd the ratio of the total capacity of the energy storage power unit cell_r；

S3200) setting judgment threshold value R of total capacity proportion of energy storage power source unit battery state of charge₁’～R₆'; wherein, R is more than 0₁’＜R₂’＜R₃’＜R₄’＜R₅’＜R₆’＜1、R₁’+R₆’＝1、R₂’+R₅’＝1、R₃’+R₄’＝1；

S3300) judging the total battery electric quantity state of the energy storage power supply unit according to the judgment threshold;

s3400) setting judgment threshold R of battery state-of-charge capacity ratio of energy storage unit₁～R₄(ii) a Wherein, R is more than 0₁＜R₂＜R₃＜R₄＜1、R₁+R₄＝1、R₂+R₃＝1；

S3500) auxiliary calculation parameters of the adjustment coefficients of the energy storage units of the energy storage power supply unit are set: s3510) setting 4 threshold parameters K₁、K₂、K₃、K₄Wherein 0 < K₁＜K₂＜K₃＜K₄(ii) a S3520) setting a variable gradient parameter delta K of an energy storage unit adjusting coefficient, wherein delta K is more than 0 and less than min K₁,K₂－K₁,K₃－K₂,K₄－K₃]Wherein min 2]Setting delta K for taking a minimum function to prevent the dynamic stability reduction of the real power value of the unit active power caused by the excessively severe change of the adjusting coefficient of the energy storage unit in the adjusting process;

s3600) calculating the adjusting coefficient of each energy storage unit of the energy storage power supply unit: calculating the upward adjustment coefficient of each energy storage unit of the energy storage power supply unit and calculating the downward adjustment coefficient of each energy storage unit of the energy storage power supply unit;

s3700) performing unit-level AGC distribution on the unit active power target value of the energy storage power supply unit: s3710) when the unit active power target value of the energy storage power supply unit is equal to 0, the single-machine active power set value of each energy storage unit is equal to 0;

s3720) when the unit active power target value of the energy storage power supply unit is larger than 0, distributing the single machine active power set value of each energy storage unit according to the mutual proportion of the product of the upward adjustment coefficient of each energy storage unit and the battery capacity, and if the calculation result is larger than the positive single machine active power rated capacity of the energy storage unit, taking the positive single machine active power rated capacity of the energy storage unit as the single machine active power set value;

s3730) when the unit active power target value of the energy storage power supply unit is less than 0, distributing the single-machine active power set value of each energy storage unit according to the mutual proportion of the product of the downward adjustment coefficient and the battery capacity of each energy storage unit; if the calculation result is smaller than the negative single-machine active power rated capacity of the energy storage unit, taking the negative single-machine active power rated capacity of the energy storage unit as a single-machine active power set value;

s3800) an active power control system of each energy storage unit of the energy storage power supply unit, aiming at a single-machine active power set value, outputting continuous signals to adjust the single-machine active power real value of the energy storage unit according to the deviation between the single-machine active power real value and the single-machine active power set value, so that the single-machine active power real value of the energy storage unit tends to the single-machine active power set value and is finally stabilized in the adjustment dead zone range of the single-machine active power set value;

s3900) calculating the unit active power rated capacity of the energy storage power supply unit:

s3910) calculating the upward adjusting capacity of each energy storage unit of the energy storage power supply unit;

s3920) accumulating the upward adjusting energy of each energy storage unit to obtain the active power rated capacity of the forward unit of the energy storage power supply unit;

s3930) calculating the downward adjusting capacity of each energy storage unit of the energy storage power supply unit;

s3940) accumulating the downward regulating energy of each energy storage unit to obtain the negative direction unit active power rated capacity of the energy storage power supply unit.

6. The method for controlling the multi-energy power supply with consideration of both the battery state and the active power dynamic stability according to claim 5, wherein the adjustment coefficients of the energy storage units of the energy storage power supply unit are modified as follows:

s1610) correcting upward adjustment coefficients of energy storage units of the energy storage power supply unit:

s1611) initializing upward adjustment coefficients of energy storage units of the energy storage power supply unit

In the formula

The upward adjustment coefficient of the energy storage unit i is obtained;

s1612) correcting the upward adjustment coefficients of the energy storage units in a fixed period: in the cycle period, the effective threshold parameter of upward adjustment of each energy storage unit is calculated first

When the content is less than or equal to 0_ri＜R₁Time of flight

When R is₁≤_ri＜R₂Time of flight

When R is₂≤_ri≤R₃Time of flight

When R is₃＜_ri≤R₄Time of flight

When R is₄＜_riWhen the temperature is less than or equal to 1

Then compare

And

time of flight

time of flight

S1620) correcting downward adjustment coefficients of energy storage units of the energy storage power supply unit:

s1621) initializing downward adjustment coefficients of energy storage units of the energy storage power supply unit

In the formula

The downward adjustment coefficient of the energy storage unit i is obtained;

s1622) correcting the downward adjustment coefficients of the energy storage units according to a fixed period: in the cycle period, the effective threshold parameter of each energy storage unit which is adjusted downwards is calculated first _ikWhen 0 is less than or equal to r_i＜R₁Time of flight _ik＝K₄When R is₁≤r_i＜R₂Time of flight _ik＝K₃When R is₂≤r_i≤R₃Time of flight _ik＝K₂When R is₃＜r_i≤R₄Time of flight _ik＝K₁When R is₄＜r_iWhen the temperature is less than or equal to 1 _ik＝0；

Then compare

time of flight

time of flight

The energy storage power supply unit also monitors the unit active power rated capacity of the energy storage power supply unit in real time:

s3910) calculating the upward adjusting capacity of each energy storage unit of the energy storage power supply unit, and when the upward adjusting effective threshold value parameter of the energy storage unit

effective threshold parameter when upward adjustment of energy storage unit

Then divided by K₂；

s3930) calculating the downward adjustment capability of each energy storage unit of the energy storage power supply unit, and when the downward adjustment of the energy storage units becomes effective, calculating the downward adjustment effective threshold parameter _ik≥K₂Then, the downward regulating capacity of the unit is the negative single-machine active power rated capacity of the unit;

when the effective threshold value parameter of the energy storage unit is adjusted downwards _ik＜K₂The downward regulating capacity of the unit is the negative single-machine active power of the unitRated capacity multiplied by _ikThen divided by K₂；

7. The method of claim 1 for controlling a multi-energy power supply with consideration of both battery state and active power dynamic stability, wherein the operation of the photovoltaic power supply unit comprises:

s4100) generating future T for each photovoltaic unit₁Possible fluctuation range of active power in time is calculated, and possible fluctuation range of unit active power of the photovoltaic power supply is calculated, wherein T₁The method is a parameter set for reserving sufficient time for possible startup and shutdown operations of the photovoltaic unit:

s4200) respectively generating a startup and shutdown sequence of the photovoltaic unit:

s4220) generating a startup sequence of available photovoltaic units without power generation, wherein the priority is calculated according to the duration of the units in a non-power generation state, and the longer the duration of the units in the non-power generation state is, the higher the priority is;

s4300) respectively generating possible active power fluctuation range sequences corresponding to the startup and shutdown sequences of the photovoltaic unit: s4310) generating an active power possible fluctuation range sequence corresponding to the starting sequence aiming at the photovoltaic unit; s4320) generating an active power possible fluctuation range sequence corresponding to the shutdown sequence of the photovoltaic unit;

s4400) calculating real unit active power value sending parameters of the photovoltaic power supply unit;

s4500) calculating a unit active power real-emitting value of the photovoltaic power supply unit and a calculated value filtering value:

s4520) calculating a filtering threshold of an active power real-sending value of the photovoltaic power supply unit;

s4530) comparing the real active power value of the photovoltaic power supply unit with the calculated value filter and the real active power value of the photovoltaic power supply unit according to a fixed period, and calculating to obtain the real active power value of the photovoltaic power supply unit and the calculated value filter;

s4600) calculating a unit primary frequency modulation target adjustment amount of the photovoltaic power supply unit:

s4610) calculating a power grid frequency deviation;

s4620) if the absolute value of the power grid frequency deviation is smaller than or equal to a given primary frequency modulation threshold, the unit primary frequency modulation target regulating quantity of the photovoltaic power supply unit is equal to 0;

s4630) if the absolute value of the power grid frequency deviation is larger than a primary frequency modulation threshold, the unit primary frequency modulation target regulating quantity of the photovoltaic power supply unit is equal to the real power output value of the photovoltaic power supply unit multiplied by the power grid frequency deviation and then multiplied by a given photovoltaic primary frequency modulation regulating coefficient.

8. The multi-energy power supply control method giving consideration to both battery state and active power dynamic stability according to claim 7, wherein the operation of the photovoltaic power supply unit is specifically:

s4100) if the photovoltaic power supply unit is provided with the power prediction system, adopting the future T of each photovoltaic unit output by the power prediction function₁The possible fluctuation range of the active power over time;

if the power prediction system is not deployed, the following method is adopted:

s4121) for the photovoltaic set of power generation, using the current power multiplied by the upper prediction parameter as the future T₁Using the current power multiplied by a lower limit prediction parameter as the lower limit value of the possible fluctuation range of the active power, wherein the upper limit prediction parameter is more than 1 and the lower limit prediction parameter is more than 0; the upper limit prediction parameter and the lower limit prediction parameter adopt fixed values or set dynamic parameters according to prior experience;

s4122) for photovoltaic units not generating electricity, using future T of generator units with performance consistent with or similar to that of photovoltaic units₁The possible fluctuation range of the active power in time is used as the future T of the unit₁The possible fluctuation range of active power in time;

s4130) calculating future T₁The unit active power of the photovoltaic power supply unit in time may fluctuate within a range: will be T in future₁Accumulating and summing the upper limits of possible fluctuation ranges of the active power of all the generator sets of the photovoltaic power supply unit within time to serve as the upper limit of the possible fluctuation ranges; will be T in future₁Accumulating and summing the lower limits of possible fluctuation ranges of the active power of all the generator sets of the photovoltaic power supply unit within the time to serve as the lower limit of the possible fluctuation ranges;

s4200) generating the startup and shutdown sequence for each photovoltaic module includes:

s4210) generating a shutdown sequence of the power generation photovoltaic unit, wherein the priority is calculated according to the duration of the unit in the power generation state, and the longer the duration of the unit in the power generation state is, the higher the priority is;

s4300) respectively generating possible fluctuation range sequences of active power corresponding to the startup and shutdown sequences for the photovoltaic unit comprises:

s4311) setting variable u₁，u₁Is 1;

s4312) adding the possible fluctuation range of the active power of the photovoltaic power supply unit to the sequence u in the photovoltaic startup sequence₁The possible fluctuation range of the active power of the unit is obtained, and the sequence u in the possible fluctuation range sequence of the active power corresponding to the photovoltaic startup sequence is obtained₁In which u is ordered₁The upper limit of the range of (a) is equal to the upper limit of the possible fluctuation range of the active power of the photovoltaic power supply unit plus the sequence u in the photovoltaic startup sequence₁The upper limit of the possible fluctuation range of the active power of the unit is sorted u₁The lower limit of the range of (1) is equal to the lower limit of the possible fluctuation range of the active power of the photovoltaic power supply unit plus the upper lightSorting u in the boot sequence₁The lower limit of the possible fluctuation range of the active power of the unit;

s4313) determination of u₁Whether it is equal to the length of the photovoltaic power-on sequence, if u₁Equal to the length of the photovoltaic power-on sequence, terminate step S4310, otherwise execute u₁＝u₁+1, and then continuing to perform the subsequent steps;

s4314) sorting u in the possible fluctuation range sequence of the active power corresponding to the photovoltaic startup sequence₁Range of-1, plus sequence u in the photovoltaic power-on sequence₁The possible fluctuation range of the active power of the photovoltaic unit is obtained, and the sequence u in the possible fluctuation range sequence of the active power corresponding to the photovoltaic startup sequence is obtained₁In which u is ordered₁Is equal to the rank u₁Upper limit of range of-1 plus the order u in the photovoltaic power-on sequence₁The upper limit of the possible fluctuation range of the active power of the unit is sorted u₁Is equal to the rank u₁Lower bound of range of-1 plus the order u in the photovoltaic boot sequence₁The lower limit of the possible fluctuation range of the active power of the unit;

s4315) jumping to step S4313 until u₁Ending when the length of the photovoltaic startup sequence is equal to the length of the photovoltaic startup sequence;

s4320) generating a possible fluctuation range sequence of active power corresponding to the shutdown sequence for the photovoltaic unit comprises:

s4321) setting variable u₂，u₂Is 1;

s4322) subtracting the sequencing u in the photovoltaic shutdown sequence from the possible fluctuation range of the active power of the photovoltaic power supply unit₂The possible fluctuation range of the active power of the photovoltaic unit is obtained, and the order u in the possible fluctuation range sequence of the active power corresponding to the photovoltaic shutdown sequence is obtained₂In which u is ordered₂Is equal to the upper limit of the possible fluctuation range of the active power of the photovoltaic power supply unit minus the sequence u in the photovoltaic shutdown sequence₂The upper limit of the possible fluctuation range of the active power of the unit is sorted u₂Is equal to the lower limit of the possible fluctuation range of the active power of the photovoltaic power supply unit minus the sequence u in the photovoltaic shutdown sequence₂The lower limit of the possible fluctuation range of the active power of the unit;

s4324) sorting u in the possible fluctuation range sequence of the active power corresponding to the photovoltaic shutdown sequence₂Range of-1, minus the order u in the photovoltaic shutdown sequence₂The possible fluctuation range of the active power of the unit is obtained, and the sequence u in the possible fluctuation range sequence of the active power corresponding to the photovoltaic shutdown sequence is obtained₂In which the order u is₂Is equal to the rank u₂Upper limit of range of-1 minus the order u in the photovoltaic shutdown sequence₂The upper limit of the possible fluctuation range of the active power of the unit is sorted u₂Is equal to the rank u₂Lower bound of range of-1 minus the order u in the photovoltaic shutdown sequence₂The lower limit of the possible fluctuation range of the active power of the unit;

s4325) jumping to step S4323 until u₂Ends equal to the photovoltaic shutdown sequence length;

s4400) calculating the real unit active power value sending parameter calculation quantity of the photovoltaic power supply unit as follows:

s4420) output dead zones of all the units of the photovoltaic power supply unit are set and accumulated to obtain unit output dead zones of the photovoltaic power supply unit;

s4430) comparing the real active power value of the photovoltaic power supply unit with the calculated quantity and the real active power value of the photovoltaic power supply unit at the current period according to a fixed period:

s4432) if the absolute value of the difference value of the two is larger than the output dead zone of the photovoltaic power supply unit, the active power real-time value parameter of the photovoltaic power supply unit is equal to the active power real-time value of the photovoltaic power supply unit in the current period;

s4500) calculating a unit active power real-emission value parameter and a calculated quantity filtering value of the photovoltaic power supply unit as follows:

s4521) setting a scaling coefficient lambda, wherein lambda is larger than 1;

s4522) the filtering threshold of the active power real output value of the photovoltaic power supply unit is equal to the unit output dead zone multiplied by lambda in S4420;

s4530) comparing the real active power value of the photovoltaic power supply unit with the calculated value of the filter and the real active power value of the photovoltaic power supply unit at the current period according to a fixed period:

s4532) if the absolute value of the difference value between the real power value and the real power value is greater than the filtering threshold obtained in S4522, the real power value of the photovoltaic power supply unit and the calculated value of the filter are equal to the real power value of the real power of the photovoltaic power supply unit in the current period;

s4600) the unit primary frequency modulation target adjustment amount of the photovoltaic power supply unit is:

s4610) the power grid frequency deviation is equal to the power grid rated frequency minus the power grid real-time frequency;

s4620) if the absolute value of the power grid frequency deviation is smaller than or equal to a primary frequency modulation threshold, the unit primary frequency modulation target regulating quantity of the photovoltaic power supply unit is equal to 0;

s4630) if the absolute value of the power grid frequency deviation is larger than a primary frequency modulation threshold, the unit primary frequency modulation target regulating quantity of the photovoltaic power supply unit is equal to the real power output value of the photovoltaic power supply unit multiplied by the power grid frequency deviation and then multiplied by a photovoltaic primary frequency modulation regulating coefficient given by the power grid.

9. The method of claim 1, wherein the operation of the complementary integrated units comprises:

s5110) setting charge-discharge parameter alpha₁And an emergency charge-discharge parameter alpha₂Wherein 0 < alpha₁＜α₂；

S5120) calculating a charge-discharge coefficient alpha every fixed period according to the total electric quantity state of the battery of the energy storage power supply unit: s5121) when the total amount of the battery is in an extremely ideal state of charge, the charge-discharge coefficient α is 0;

S5126) when the total amount of the battery is in a more ideal electric quantity state, keeping the original value of the charge-discharge coefficient unchanged;

s5130) calculating the charging and discharging correction power of the energy storage power supply unit according to the charging and discharging coefficient;

the charging and discharging corrected power is

S5200) calculating a unit active power target value of the thermal power supply, wherein the unit active power target value of the thermal power supply is equal to a total active power set value of the multi-energy complementary integrated power supply minus a unit active power actual value of the photovoltaic power supply unit to participate in a calculated quantity filter value, and then charging and discharging correction power of the energy storage power supply unit;

comparing the unit active power target value of the thermal power supply with the unit joint operation area of the thermal power supply unit, wherein when the unit active power target value is contained in the unit joint operation area, the unit active power target value is feasible; the complementary integration unit sends the power to a firepower power supply unit, and the firepower supply unit carries out unit-level AGC distribution on a unit active power target value;

s5210) calculating a primary frequency modulation adjustment coefficient of the thermal power supply unit, wherein the primary frequency modulation adjustment coefficient is used when each unit of the thermal power supply unit actually executes a primary frequency modulation task; the thermal power supply unit adjusts the active power of each thermal power supply single-machine closed-loop unit;

s5220) when the unit active power target value is not included in the unit joint operation region, the unit active power target value of the fire power source unit is not feasible, finding an operation proposal that makes the unit active power target value feasible:

s5221) sequentially finding a running operation advice that makes the unit active power target value of the thermal power source feasible by putting the unit not put into AGC control, finding a running operation advice that makes the unit active power target value of the thermal power source feasible by turning the unit not put into AGC into a generating state and putting into AGC, finding a running operation advice that makes the unit active power target value of the thermal power source feasible by turning the unit generating power into a non-generating state; classifying the found operation suggestions, and displaying the operation suggestions in order according to the priority to generate operation suggestions aiming at the thermal power supply unit;

the complementary integration unit sends the generated operation suggestions of the thermal power supply unit to the thermal power supply unit, and the thermal power supply unit adjusts the active power of each thermal power supply unit closed-loop unit;

s5300) generating a proposal for startup and shutdown of the photovoltaic power supply unit and assigning the proposal to the photovoltaic power supply unit:

after the unit active power target value of the thermal power supply unit is determined, the on-off state and the future T of the current photovoltaic power supply unit are obtained₁The mismatch degree quantization value of the total active power set value of the multi-energy complementary integrated power supply in time;

then the current photovoltaic power supply unit is started or stopped and the future T₁Multiple functions within timeThe mismatching degree quantization value of the total active power set value of the complementary integrated power supply, the possible active power fluctuation range sequence corresponding to the start-up and shut-down sequence of the photovoltaic unit and the future T₁Comparing the mismatch quantitative values of the total active power set values of the multi-energy complementary integrated power supply within the time, and generating an operation suggestion according to the comparison result; orderly displaying the generated operating suggestions of the photovoltaic units during the startup and shutdown and assigning the operating suggestions to the photovoltaic power supply units;

s5400) the complementary integrated unit calculates the unit active power target value of the energy storage power supply unit and allocates:

s5410) adding a total active power set value of the multi-energy complementary integrated power supply to a unit primary frequency modulation correction quantity of the thermal power supply unit, then subtracting a unit active power actual emission value of the photovoltaic power supply unit from a calculated quantity, and then subtracting the unit active power actual emission value of the thermal power supply unit to obtain an active power total adjustment deviation of the thermal power supply unit and the photovoltaic power supply unit;

s5420) initially setting the compensation adjustment quantity of the energy storage power supply unit as the active power total adjustment deviation, and then comparing the compensation adjustment quantity of the energy storage power supply unit with the current active power total adjustment deviation according to a fixed period, wherein the method comprises the following steps:

s5421) when the absolute value of the difference value of the two is larger than the unit active power regulation dead zone of the energy storage power supply unit, the compensation regulation quantity of the energy storage power supply unit is equal to the total regulation deviation of the current active power;

s5422) when the absolute value of the difference value of the two is less than or equal to the unit active power regulation dead zone of the energy storage power supply unit, the compensation regulation quantity of the energy storage power supply unit is kept unchanged;

s5430) dead zone processing is carried out on the compensation adjustment quantity of the energy storage power supply unit:

s5431) manually setting the timer and the time parameter T₄；

S5432) when the absolute value of the active power total regulation deviation is less than or equal to the active power regulation dead zone of the thermal power supply unit, the timer set in S5431 starts to time;

s5433) resetting and clearing a timer set in S5431 when the absolute value of the total active power regulation deviation obtained in S5410 is larger than the active power regulation dead zone of the fire unit;

s5435) when the timer time is more than or equal to the time parameter T₄When the active power target value of the energy storage power supply unit is equal to 0;

the complementary integration unit distributes the unit active power target value after the dead zone processing to the energy storage power supply unit; and the energy storage power supply unit performs unit-level AGC distribution according to the allocated unit active power target value, and adjusts the active power of each energy storage unit.

10. The method of claim 9 for controlling a multi-energy power supply with consideration of both battery state and active power dynamic stability, wherein generating a photovoltaic unit start-up and shut-down operation recommendation comprises:

s5310) calculating future T₁Unit active power accommodation range of photovoltaic power supply unit in time, wherein T₁Time parameters set for manual:

s5311) calculating future T₁The unit active power of photovoltaic power supply unit of each time point in time holds the scope lower limit, includes:

if the active power plan curve of the multi-energy complementary integrated power supply is issued in advance by scheduling, the future T is determined₁Subtracting the upper limit of the combined operation area of the thermal power supply unit from the total active power set value of the multi-energy complementary integrated power supply at each time point in time to obtain the future T₁The lower limit of the unit active power accommodation range of the photovoltaic power supply unit at each time point in time;

if the active power plan curve of the multi-energy complementary integrated power supply is not issued in advance in the scheduling process, subtracting the upper limit of the combined operation area of the thermal power supply unit from the total active power set value of the current multi-energy complementary integrated power supply to obtain the future T₁The lower limit of the unit active power accommodation range of the photovoltaic power supply unit at each time point in time;

s5312) calculating future T₁Photovoltaic power supply unit at each time point in timeThe upper limit of the unit active power accommodation range of (1), comprising:

if the active power plan curve of the multi-energy complementary integrated power supply is issued in advance by scheduling, the future T is determined₁Subtracting the lower limit of the combined operation area of the thermal power supply unit from the total active power set value of the multi-energy complementary integrated power supply at each time point in time to obtain the future T₁The upper limit of the unit active power accommodation range of the photovoltaic power supply unit at each time point in time;

if the active power plan curve of the multi-energy complementary integrated power supply is not issued in advance in the scheduling process, subtracting the lower limit of the combined operation area of the thermal power supply unit from the total active power set value of the current multi-energy complementary integrated power supply to obtain the future T₁The upper limit of the unit active power accommodation range of the photovoltaic power supply unit at each time point in time;

s5313) future T₁The unit active power accommodation range of the photovoltaic power supply unit in time is T in the future₁Taking intersection of unit active power accommodation ranges of the photovoltaic power supply units at each time point in time;

s5321) calculating future T₁Unit active power accommodation range and future T of photovoltaic power supply unit in time₁The upper limit mismatching degree of the possible fluctuation range of the active power of the photovoltaic power supply unit in time is T in the future₁Subtracting future T from the upper limit of the possible fluctuation range of the active power of the photovoltaic power supply unit in time₁Judging the upper limit of the unit active power accommodation range of the photovoltaic power supply unit within time, if the upper limit active power accommodation range is larger than 0, determining that the upper limit mismatching degree is equal to the calculation result, and otherwise, determining that the upper limit mismatching degree is equal to 0;

s5322) calculating future T₁The unit active power accommodation range of the photovoltaic power supply unit in time includes and is T in future₁The lower limit mismatching degree of the possible fluctuation range of the active power of the photovoltaic power supply unit in time is T in the future₁Lower limit of unit active power accommodation range of photovoltaic power supply unit in timeMinus future T₁Judging the calculation result according to the lower limit of the possible fluctuation range of the active power of the photovoltaic power supply unit in time, wherein if the lower limit is larger than 0, the lower limit mismatching degree is equal to the calculation result, and otherwise, the lower limit mismatching degree is equal to 0;

s5323) subtracting the lower limit mismatching degree obtained by S5322 from the upper limit mismatching degree obtained by S5321, and taking the absolute value of the result to obtain the on-off state and the future T of the current photovoltaic power supply unit set₁The mismatch degree quantization value of the total active power set value of the multi-energy complementary integrated power supply in time;

s5330) finding an operational recommendation to shutdown the photovoltaic power generating set:

s53320) setting variable v₃，v₃Is 1;

s53341) calculating future T₁Sequencing v in the possible fluctuation range sequence of the active power corresponding to the photovoltaic shutdown sequence and the unit active power accommodation range of the photovoltaic power supply unit in time₃The upper limit mismatching degree of the range, and sorting v in the possible fluctuation range sequence of the active power corresponding to the photovoltaic shutdown sequence₃Upper range limit minus the future T₁Judging the calculation result according to the upper limit of the active power accommodation range of the photovoltaic power supply unit in time, wherein if the upper limit is larger than 0, the upper limit mismatching degree is equal to the calculation result, and otherwise, the upper limit mismatching degree is equal to 0;

s53342) calculating future T₁Active power corresponding to photovoltaic shutdown sequence in unit active power accommodation range of photovoltaic power supply unit in timeOrdering v in a sequence of possible fluctuation ranges of the rate₃The lower limit mismatch of the range of (1), will be T in the future₁Sequencing v in the sequence of subtracting the possible fluctuation range of the active power corresponding to the photovoltaic shutdown sequence from the lower limit of the active power accommodation range of the unit of the photovoltaic power supply unit in time₃If the lower limit is larger than 0, the lower limit mismatching degree is equal to the calculation result, otherwise, the lower limit mismatching degree is equal to 0;

s53343) subtracting the lower limit mismatching degree obtained by S53342 from the upper limit mismatching degree obtained by S53341, and taking an absolute value of the result to obtain a sequence v in the possible fluctuation range sequence of the active power corresponding to the photovoltaic shutdown sequence₃Range and future T of₁The mismatch degree quantization value of the total active power set value of the multi-energy complementary integrated power supply in time;

s53352) if the calculation result is less than the judgment threshold parameter set in S53310, jumping to S53360 and continuing to execute;

s53361) if v₃If 1, no operation suggestion is generated;

s5340) finding an operational recommendation to start up an available and non-generating photovoltaic power plant:

s53420) setting variable v₄，v₄Is 1;

s53430) if v₄If the length of the starting sequence of the photovoltaic unit is less than or equal to the length of the starting sequence of the photovoltaic unit, setting an original mismatching degree quantization value variable, wherein the original mismatching degree quantization value variable is equal to the mismatching degree quantization value obtained in the step S5320, otherwise, skipping to the step S53460;

s53443) subtracting the lower limit mismatching degree obtained by the S53442 from the upper limit mismatching degree obtained by the S53441, and taking an absolute value of the result to obtain a sequence v in a possible fluctuation range sequence of active power corresponding to the startup sequence of the photovoltaic unit₄Range and future T of₁The mismatch degree quantization value of the total active power set value of the multi-energy complementary integrated power supply in time;

s53452) if the calculation result is less than the judgment threshold parameter set in S53410, jumping to S53460 and continuing to execute;

s53461) if v₄If 1, no operation suggestion is generated;

s53462) if v₄If the number of the photovoltaic units is more than 1, generating a starting operation suggestion, and sequencing 1 to v in a starting sequence of the photovoltaic units according to the suggestion₄-1, the corresponding photovoltaic unit executes a startup operation;

then orderly displaying the shutdown operation suggestions of the photovoltaic units generated by the S5330 respectively;

and respectively and orderly displaying the photovoltaic unit startup operation suggestions generated by the S5340.