CN109936164A - Multiple-energy-source electric power system optimization operation method based on the analysis of power supply complementary characteristic - Google Patents

Abstract

A kind of multiple-energy-source electric power system optimization operation method based on the analysis of power supply complementary characteristic of the invention, its main feature is that, it include: to establish to provide multiple forms of energy to complement each other in multiple-energy-source electric system to coordinate the complementary mechanisms of power generation, mathematical model based on complementary mechanisms building complementary index system and complementary demand, define renewable energy complementary power supply, it calculates with matching water power capacity ratio required in the renewable energy complementary power supply of the minimum target of complementary demand, the hierarchy optimization operation reserve of polyisocyanate mass-energy source current is formulated according to complementary index, the corresponding optimization object function of each optimization layer is solved using particle swarm algorithm, calculating can make complementary index be optimal corresponding polyisocyanate mass-energy source current in the output power value of day part, this method has science, rationally, simply, it is practical, it is able to ascend renewable energy Dissolve horizontal advantage.

Description

Optimal operation method of multi-energy power system based on analysis of complementary characteristics of power sources

technical field

The invention relates to the field of optimal operation of electric power systems, and is an optimal operation method of a multi-energy electric power system based on the analysis of complementary characteristics of power sources.

Background technique

With the outbreak of the energy crisis and the prominence of environmental pollution, the development and utilization of renewable energy has received unprecedented attention. However, due to the intermittent and fluctuating output power of wind power and photovoltaics, it brings certain challenges to the optimal operation of the power grid. , Through the complementary characteristics between heterogeneous power sources, the impact of wind power and photovoltaic grid connection can be weakened, which plays an important role in reducing environmental pollution, improving energy utilization, and ensuring the stable operation of the power system. Therefore, utilizing the space-time complementary characteristics between power sources is an important means to improve the capacity of renewable energy consumption.

Existing researches on the optimal operation of power systems with new energy include using energy storage to cope with the uncertainty of renewable energy and optimizing the power system, but do not consider the complementary characteristics of other power sources; some only consider two Or the complementary characteristics between the three power sources to optimize the power system for energy saving, and to optimize the power system on multiple time scales from the perspective of source-grid-load coordination, but there are fewer types of power sources to consider. The power system is optimized from the perspective of peak regulation of the power system. In the power system optimization involving multiple energy sources, the existing research fails to make good use of the complementary characteristics of multiple power sources to improve the consumption of renewable energy.

SUMMARY OF THE INVENTION

The purpose of the present invention is to provide a scientific, reasonable, simple and practical method for optimizing the operation of a multi-energy power system based on the analysis of the complementary characteristics of the power source.

The technical solution adopted to realize the purpose of the present invention is a method for optimizing operation of a multi-energy power system based on the analysis of the complementary characteristics of the power supply, which is characterized in that it comprises the following steps:

1) Establish a complementary mechanism for multi-energy complementary and coordinated power generation in a multi-energy power system

① In a multi-energy power system, a variety of heterogeneous energy sources with different output power characteristics are included, which is a prerequisite for multi-energy complementary and coordinated power generation. The power supply complementary characteristic is based on the characteristics of mutual aid and mutual power generation between different heterogeneous energy sources. Formula (1) is used to express the power supply complementary characteristic as the characteristic that the output power of various power supplies meets the system load.

In the formula, is the load value of the multi-energy power system in the t-th period; is the output power value of the i-th thermal power unit in the multi-energy power system in the t-th period, i=1,2,3...N _th , N _th is the in-service thermal power unit in the multi-energy power system in the t-th period quantity; is the output power value of the j-th hydroelectric unit in the multi-energy power system in the t-th period, j=1,2,3...N _hy , N _hy is the multi-energy power system in the t-th period of the hydroelectric unit number of service; is the output power value of the k-th photovoltaic unit in the multi-energy power system in the t-th period, k=1,2,3...N _pv , where N _pv is the output power of the photovoltaic unit in the t-th period of the multi-energy power system number of service; is the output power value of the g-th wind turbine in the multi-energy power system in the t-th period, g=1, 2, 3...N _w , N _w is the multi-energy power system in the t-th period of the wind turbine number of times; t=1, 2, 3...T, T is the number of time periods, Δt is the time step,

② Establish a complementary mechanism for multi-energy complementary and coordinated power generation in a multi-energy power system: based on the natural complementary characteristics of various power sources, relying on the good regulation capabilities of thermal power, gas power, and adjustable hydropower to stabilize wind power and photovoltaics. The fluctuation of the output power of resource-constrained power sources, making full use of clean renewable energy sources, reducing the proportion of thermal power in the power system, realizing the rational allocation of power system resources, and finally achieving a real-time balance between the total power generation and total load of the system. To optimize the operation purpose,

2) Constructing complementary index system and mathematical model of complementary demand respectively

①Build a complementary index system

The complementarity index is defined as the quantitative index of the complementary effect pursued by the multi-energy power system, that is, the direction of its optimization. Combined with the complementary mechanism of multi-energy complementary and coordinated power generation in the multi-energy power system, it can improve the system's renewable energy consumption capacity and energy saving. From the perspective of efficiency enhancement, build a complementary index system,

A. Calculate the penetration rate index of renewable energy in the multi-energy power system

Define r _pe (t) as the ratio of the output power of the renewable energy source to the load in the t-th period of the multi-energy power system, which mainly reflects the status of the renewable energy source in the power system and the penetration rate of the output power of the renewable energy source. The larger the value is, the cleaner the multi-energy power system is. Formula (2) is used to calculate the penetration rate of renewable energy in the multi-energy power system.

In the formula, is the penetration rate index of renewable energy in the multi-energy power system, represents the output power value of the renewable energy power source in the multi-energy power system in the t-th period, represents the load value of the multi-energy power system in the t-th period, T is the number of periods, Δt is the time step, t=1, 2, 3...T,

B. Calculate the coal consumption index of thermal power

In a multi-energy power system, the smaller the coal consumption of thermal power, the better the economy and environmental protection of thermal power. Formula (3) is used to calculate the coal consumption index of thermal power units,

In the formula, a _i , bi , c _i are the three fuel consumption characteristic coefficients of the _ith thermal power unit, i=1, 2, 3...N _th , N _th is the thermal power in the t-th time period in the multi-energy power system The number of units in service, f _it is the coal consumption index of the i-th thermal power unit in the t-th period, P _th.it is the output power value of the i-th thermal power unit in the t-th period, t=1,2, 3...T, Δt is the time step, T is the number of time periods,

C. Calculate the fluctuation range index of the load borne by thermal power

In order to reduce the frequent changes of output power of thermal power units, reduce coal consumption and improve utilization efficiency, the fluctuation range of load borne by thermal power should be minimized. The standard deviation of load value is used to represent the fluctuation range of load, and formula (4) is used to calculate thermal power The fluctuation range index of the unit's load,

In the formula, δ _th is the fluctuation range index of the load borne by thermal power units in the multi-energy power system, t = 1, 2, 3...T, Δt is the time step, T is the number of time periods, is the load value borne by all thermal power in the multi-energy power system in the t-th period, is the average load of all thermal power in the multi-energy power system during T periods,

D. Use formula (5) to calculate the power generation of hydropower,

In the formula, W _H is the total power generation of all hydropower in the multi-energy power system in T periods; is the water consumption of the jth hydroelectric unit for power generation in the tth period, is the head height of the j-th hydroelectric unit in the t-th period, j=1,2,3...N _hy , N _hy is the number of in-service hydro-electric units in the multi-energy power system in the t-th period, η is the hydroelectric conversion efficiency ;t=1,2,3...T, Δt is the time step, T is the number of time periods,

E. Use formula (6) to calculate the water abandonment index of hydropower,

In the formula, ΔQ is the index of water abandonment of all hydropower in T periods, Q _jmax is the maximum allocated water consumption of the jth hydropower unit in T periods, is the water consumption of the j-th hydroelectric unit in the t-th period, t=1, 2, 3...T, Δt is the time step, T is the number of periods, j=1, 2, 3...N _hy , N _hy is The number of in-service hydropower units in the multi-energy power system in the t-th period,

②Construct a mathematical model describing the complementary requirements of multiple heterogeneous power sources

The complementary demand of multi-heterogeneous power supplies is defined as the degree of matching between the output powers of the various heterogeneous power supplies after complementing each other with the load within a certain period of time. The two elements in the complementary demand are the power supply and the load. quantified between the complementary demand index and the complementary demand index between the power source and the load,

A. Calculate the complementary demand indicators between various heterogeneous power sources

A1. Use formula (7) to calculate the rate of change of the output power of the power supply,

r _i ^t =(P _i ^t -P _i ^t-1 )/Δt, (7)

In the formula, r _i ^t is the rate of change of the output power of the i-th power supply in the t-1th period to the t-th period, P _i ^t is the output power value of the i-th power source in the t-th period, P _i ^{t -1} is the output power value of the i-th power supply in the t-1th period, i=1,2,3...n, n is the number of power supply types under investigation, t=1,2,3...T; Δt is the time step size, T is the total number of epochs,

A2. Use formula (8) to calculate the set of absolute values of the sum of the output power change rates of thermal power, hydropower, photovoltaic, and wind power in each period,

In the formula, S _s is the set of absolute values of the sum of the output power change rates of thermal power, hydropower, photovoltaic, and wind power in each time period of T, and β _t is the thermal power, hydropower, photovoltaic, and wind power in the t-1th time period. The absolute value of the sum of the output power change rates in the t-th period, is the output power change rate of thermal power from the t-1th period to the tth period, is the output power change rate of hydropower from the t-1th period to the tth period, is the rate of change of photovoltaic output power from the t-1th period to the tth period, is the output power change rate of wind power from the t-1th period to the tth period, t=1, 2, 3...T, Δt is the time step, T is the number of periods,

A3. Use formula (9) to calculate the complementary demand index between various heterogeneous power sources,

In the formula, D _ss is the complementary demand index between various heterogeneous power sources in T periods, and β _t is the output power change of thermal power, hydropower, photovoltaic, and wind power from the t-1th period to the tth period. The absolute value of the sum of the rates, t=1,2,3...T, Δt is the time step, T is the number of time periods,

The smaller the value of the complementary demand index between various heterogeneous power sources, the stronger the mutual support between the wind-solar, hydro-thermal power sources in the time scale under investigation, that is, the better the complementary effect. The weaker the mutual support, the

B. Calculate the complementary demand index between the power supply and the load

B1. Use formula (10) to calculate the relative rate of change of the output power of the power supply,

In the formula, is the relative rate of change of the total output power of all power supplies from the t-1th period to the tth period, is the output power change rate of the total output power of all power sources from the t-1th period to the tth period, P _sc is the installed capacity of the generator sets in service in all power sources, T is the number of periods, t=1,2,3 ...T, Δt is the time step,

B2. Use formula (11) to calculate the relative rate of change of the load,

In the formula, is the relative change rate of the system load from the t-1th period to the tth period, P _lmax.T is the maximum load value in the T period, is the change rate of the system load from the t-1th period to the tth period, t=1, 2, 3...T, Δt is the time step, T is the number of periods,

B3. Use formula (12) to calculate the set of absolute values of the sum of the relative change rates in each period between the output power of all power supplies and the system load,

In the formula, S _l is the set of absolute values of the relative rate of change in each period between the output power of all power supplies and the system load in T time periods, and α _t is the difference between the output power of the power supply and the system load at the t-th The absolute value of the sum of the relative rate of change from 1 period to the t period, is the relative rate of change of the total output power of all power supplies from the t-1th period to the tth period, is the relative change rate of the system load from the t-1th period to the tth period, t=1, 2, 3...T, Δt is the time step, T is the number of periods,

B4. Use formula (13) to calculate the complementary demand index between the power supply and the load,

In the formula, D _sl is the complementary demand index between the power supply and the load of the multi-energy power system in T time periods, α _t is the difference between the output power of the power supply and the system load from the t-1th period to the tth period. The absolute value of the relative change rate of the time period and the sum, T is the number of time periods, t=1, 2, 3...T, Δt is the time step,

The smaller the value of the complementary demand index between the power source and the load, the more similar the change trend of the power source and the load in the investigated time scale; otherwise, the more different the change trend of the power source and the load,

The better the complementarity between power sources and between power sources and loads, the smaller the complementary demand, that is, the closer the complementary demand index value is to zero;

3) Define renewable energy complementary power sources

The renewable energy power sources that can meet the complementary needs after complementation are aggregated into one power source, which is defined as Renewable Energy Complementary Power Supply (RECPS). The main purpose of aggregation is to reduce the power fluctuation brought by wind and solar power to the multi-energy power system. The basic rule of aggregation is to meet the complementary needs of the system under the complementary mechanism of multi-energy complementary and coordinated power generation in the multi-energy power system. Based on energy conservation and environmental protection, the renewable energy power sources are aggregated, and the hydropower capacity ratio is calculated to minimize the complementary demand. After the renewable energy complementary power source is formed, its output power can follow the load fluctuation. In a multi-energy power system , the renewable energy complementary power supply is regarded as a kind of power supply, and it operates optimally together with other conventional power supplies, and when the system load value remains unchanged, the output power of the renewable energy complementary power supply also remains unchanged;

4) Formulate a hierarchical optimization operation strategy for multi-heterogeneous energy power sources

The hierarchical optimization operation strategy of multi-heterogeneous energy power sources is based on the complementary mechanism of multi-energy complementary and coordinated power generation in the multi-energy power system, giving priority to full acceptance of wind power and photovoltaic power generation; making full use of adjustable hydropower to deal with these uncertain power sources such as wind and solar. The randomness, intermittency and anti-peak regulation characteristics brought about by the present invention aggregate hydropower, wind power, and photovoltaics into a renewable energy complementary power supply. The output power of the renewable energy complementary power supply is relatively stable and can follow the load fluctuation, which can improve the multi-energy power supply. The ability of the system to absorb the wind and the wind and reduce the adverse impact of the uncertainty of wind and wind resources on the stable operation of the system,

The selected optimization objective of the multi-energy power system is to optimize the complementary indexes of the multi-energy power system and realize the rational allocation of power system resources. Including complementary power optimization layer, residual hydropower optimization layer, thermal power optimization layer,

①Complementary power optimization layer

In the optimization operation of the multi-energy power system, the renewable energy complementary power supply is firstly optimized, and the renewable energy complementary power supply is obtained by the aggregation of wind power, photovoltaic and hydropower. Obtain the aggregate capacity ratio of matching hydropower and wind power in the renewable energy complementary power supply, and then determine the output power of wind power, photovoltaic, and matching hydropower in each time period. The main goal of the renewable energy complementary power optimization layer is to determine the wind power, wind power The aggregate proportion of

In the formula, D _sl is the complementary demand between the power source and the load of the multi-energy power system in T time periods, is the relative rate of change of the load from the t-1th period to the tth period, is the relative change rate of the output power of the renewable energy complementary power supply from the t-1th period to the tth period, t=1, 2, 3...T, Δt is the time step, T is the number of periods,

②Remaining hydropower optimization layer

The net load curve is obtained by subtracting the output power of the renewable energy complementary power supply from the load value of the multi-energy power system. Under the condition that the net load is stable, the residual hydropower is used to generate electricity with the goal of minimizing the amount of waste water, and is calculated by formula (15). The objective function of the minimum amount of water waste from hydropower,

In the formula, ΔQ is the discarded water volume of the remaining hydropower in T periods, Q _jmax is the maximum allocated water consumption of the jth hydropower unit in T periods, is the water consumption of the jth hydroelectric unit in the tth period, j=1,2,3...N _rhy , N _rhy is the total number of remaining hydroelectric units, t=1,2,3...T, Δt is the time step , T is the number of time periods,

③ Thermal power optimization layer

The residual load is obtained by subtracting the residual hydropower output power from the net load curve. Under the dual action of the renewable energy complementary power supply and residual hydropower, the fluctuation range of the residual load is small, and the output power of the thermal power unit is arranged with the goal of minimum coal consumption. When the output power of the thermal power unit reaches the minimum, and the total output power of various heterogeneous power sources is still greater than the load, part of the output power of the renewable energy power source needs to be discarded, and formula (16) is used to calculate the objective function of the minimum coal consumption for thermal power generation,

In the formula, F is the total coal consumption of the thermal power unit, u _it is the on-off coefficient of thermal power, which is 1 when starting up and 0 when shutting down, f _it is the coal consumption of the i-th thermal power unit in the t-th period, i=1 ,2,3…N _th , N _th is the total number of thermal power units, t=1,2,3…T, Δt is the time step, T is the number of time periods,

5) Determine the constraints

The optimal operation of the multi-energy power system needs to satisfy the constraint equations of equations (17) to (23),

①Determine power balance constraints

The power balance constraint is expressed by Equation (17),

In the formula, is the total thermal power output power in the t-th period, is the remaining hydropower output power in the t-th period, is the total hydropower output power in the t-th period, is the total output power of the hydroelectric unit in the renewable energy complementary power supply in the t-th period, is the output power of the renewable energy complementary power supply in the t-th period, P _l ^t is the load value of the multi-energy power system in the t-th period, t=1, 2, 3...T, Δt is the time step, T is the number of periods,

② Determine the active output power constraints of the unit

The active output power constraint of the unit is expressed by formula (18),

In the formula, is the output power value of the i-th thermal power unit in the multi-energy power system in the t-th period, i=1,2,3...N _th , N _th is the in-service thermal power unit in the multi-energy power system in the t-th period quantity; is the output power value of the j-th hydroelectric unit in the multi-energy power system in the t-th period, j=1,2,3...N _hy , N _hy is the multi-energy power system in the t-th period of the hydroelectric unit number of service; is the output power value of the k-th photovoltaic unit in the multi-energy power system in the t-th period, k=1,2,3...N _pv , where N _pv is the output power of the photovoltaic unit in the t-th period of the multi-energy power system number of service; is the output power value of the g-th wind turbine in the multi-energy power system in the t-th period, g=1, 2, 3...N _w , N _w is the multi-energy power system in the t-th period of the wind turbine P _max.i is the upper limit of the output power of the ith thermal power unit in the multi-energy power system, P _max.j is the upper limit of the output power of the j-th hydropower unit in the multi-energy power system, and P _max.k is the multi-energy power system. The upper limit of the output power of the kth photovoltaic unit in the power system, P _max.g is the upper limit of the output power of the gth wind turbine in the multi-energy power system; P _min.i is the output of the ith thermal power unit in the multi-energy power system. Power lower limit, P _min.j is the output power lower limit of the jth hydroelectric unit in the multi-energy power system, t=1, 2, 3...T, Δt is the time step, T is the number of time periods,

③ Determine the system spinning reserve capacity constraints

The system spinning reserve capacity constraint is expressed by formula (19),

where: is the spinning reserve capacity of the system in the t-th period, is the rotating reserve capacity of the i-th thermal power unit in the t-th period, i=1,2,3...N _th , where N _th is the number of thermal power units in service in the t-th period in the multi-energy power system; is the rotating reserve capacity of the j-th hydroelectric unit in the t-th period, j=1, 2, 3...N _hy , N _hy is the number of in-service hydro-electric units in the multi-energy power system in the t-th period; α is The demand coefficient of the system load forecast error to the spinning reserve, β is the demand coefficient of the wind power output power prediction error to the spinning reserve; γ is the demand coefficient of the photovoltaic output power prediction error to the spinning reserve, is the load value of the multi-energy power system in the t-th period; is the output power value of the k-th photovoltaic unit in the multi-energy power system in the t-th period, k=1,2,3...N _pv , where N _pv is the output power of the photovoltaic unit in the t-th period of the multi-energy power system number of service; is the output power value of the g-th wind turbine in the multi-energy power system in the t-th period, g=1, 2, 3...N _w , N _w is the multi-energy power system in the t-th period of the wind turbine number of times, t=1, 2, 3...T, Δt is the time step, T is the number of time periods,

④ Determine the up-slope constraint and down-slope constraint of the unit

The up-climbing constraints and down-climbing constraints of the unit are expressed by formulas (20)~(21),

P _i ^{t+1 -} P _i ^t ≤ΔP _i ^up (20)

P _i ^t -P _i ^t+1 ≤ΔP _i ^down (21)

In the formula, P _i ^t+1 is the output power of the i-th unit in the multi-energy power system in the t+1-th period, P _i ^t is the output power of the i-th unit in the multi-energy power system in the t-th period, ΔP _i ^up is the maximum value of the ith unit uphill in the multi-energy power system, ΔP _i ^down is the maximum value of the ith unit in the multi-energy power system downhill, i=1,2,3...N,N is the number of active units in the multi-energy power system, t=1, 2, 3...T, Δt is the time step, T is the number of time periods,

⑤ Determine the constraints of hydropower generation

The constraint of hydropower generation is expressed by formula (22),

where: is the output power value of the j-th hydroelectric unit in the multi-energy power system in the t-th period, j=1,2,3...N _hy , N _hy is the number of in-service hydro-electric units in the multi-energy power system in the t-th period , W _H is the total power generation of all hydropower units in the multi-energy power system in T time periods, t=1, 2, 3...T, Δt is the time step, T is the number of time periods,

⑥ Determine hydroelectric flow constraints

The hydroelectric flow constraint is expressed by Equation (23),

In the formula: Q _jmin is the minimum allocated water consumption of the jth hydropower unit in the T period, Q _jmax is the maximum allocated water consumption of the jth hydropower unit in the T period, is the water consumption of the j-th hydropower unit in the t-th period, j=1,2,3...N _hy , N _hy is the number of in-service hydro-power units of the multi-energy power system in the t-th period, t=1,2 ,3…T, Δt is the time step, T is the number of time periods,

6) Through steps 1) to 5), using the constructed complementary index system and the mathematical model of complementary demand, combined with the actual operation constraints of the multi-energy power system, the output power of the renewable energy complementary power source is most matched with the load, and the remaining The objective function of the minimum water waste of hydropower and the minimum coal consumption of thermal power is solved by particle swarm algorithm, and finally the output power value of the multi-heterogeneous energy power source in each time period can be calculated to optimize the complementary index.

The features of the method for optimizing operation of a multi-energy power system based on the analysis of the power supply complementary characteristics of the present invention include the following steps: first, establishing a complementary mechanism for multi-energy complementary and coordinated power generation in the multi-energy power system, and constructing a complementary mechanism based on the complementary mechanism Then, define the renewable energy complementary power source, which is composed of all wind and solar power sources and the required matching hydropower, and calculate the hydropower capacity ratio aiming at the minimum complementary demand; Secondly, the hierarchical optimization operation strategy of multi-heterogeneous energy power sources is formulated according to the complementary indicators, including the complementary power source optimization layer, the remaining hydropower optimization layer, and the thermal power optimization layer. Solve the problem and calculate the output power value of the multi-heterogeneous energy power source in each period that can make the complementary index reach the optimum. This method is scientific, reasonable, simple and practical, and can improve the level of renewable energy consumption. .

Description of drawings

1 is a flowchart of a method for optimizing operation of a multi-energy power system based on the analysis of complementary characteristics of power sources according to the present invention;

Figure 2 is a schematic diagram of the change curve of the complementary demand of the multi-energy power system with the aggregation ratio of hydropower and wind and solar power capacity;

Figure 3 is a comparison diagram of the output power and load of the renewable energy complementary power supply;

Figure 4 is a graph showing the results of the optimized operation of the multi-energy power system;

Figure 5 is a comparison chart of the power generation of various heterogeneous power sources before and after optimization.

Detailed ways

The present invention will be further described below with reference to the accompanying drawings and embodiments.

Referring to Figures 1-5, Figure 1 shows the basic data processing, the construction of the complementary index system and the mathematical model of complementary requirements to the formulation of the hierarchical optimization operation strategy for multi-heterogeneous energy power sources. Finally, the particle swarm algorithm is used for each optimization layer. The corresponding optimization objective function is solved, and the technical roadmap of the output power value of the multi-heterogeneous energy power source in each time period corresponding to the optimal complementarity index is calculated. Figure 3 shows the comparison between the output power and load of the optimized renewable energy complementary power supply; Figure 4 shows the results of the optimal operation of the multi-energy power system; Figure 5 shows the comparison of the power generation of various heterogeneous power sources before and after optimization.

In a method for optimizing operation of a multi-energy power system based on the analysis of the complementary characteristics of the power supply of the present invention, the parameter values of the embodiment are set as follows:

The parameters of the thermal power unit are shown in Table 1

Table 1 embodiment thermal power unit parameters

parameter TH1 TH2 TH3 TH4 TH5 Capacity (MW) 6020 270 270 1320 700 Maximum output power (MW) 6020 270 270 1320 700 Minimum output power (MW) 300 135 135 560 320 Ramp rate (MW/min) 10 2.7 2.7 6 7 a 0.0088 0.0028 0.0085 0.0022 0.0058 b 0.17 0.52 0.32 0.35 0.45 c 10.8 14.7 12.3 13.8 11.9

The total installed capacity of wind power is 1,920MW; the total installed capacity of photovoltaics is 7,954MW; Requirement coefficient for spinning reserve γ=9%; T=24; ΔT=1h.

A method for optimizing operation of a multi-energy power system based on the analysis of the complementary characteristics of the power supply of the present invention includes the following steps:

1) Establish a complementary mechanism for multi-energy complementary and coordinated power generation in a multi-energy power system

① In a multi-energy power system, a variety of heterogeneous energy sources with different output power characteristics are included, which is a prerequisite for multi-energy complementary and coordinated power generation. The power supply complementary characteristic is based on the characteristics of mutual aid and mutual aid power generation between different heterogeneous energy sources. Formula (1) is used to express the power supply complementary characteristic as the characteristic that the output power of various power sources meets the system load.

In the formula, is the load value of the multi-energy power system in the t-th period; is the output power value of the i-th thermal power unit in the multi-energy power system in the t-th period, i=1,2,3...N _th , N _th is the in-service thermal power unit in the multi-energy power system in the t-th period quantity; is the output power value of the j-th hydroelectric unit in the multi-energy power system in the t-th period, j=1,2,3...N _hy , N _hy is the multi-energy power system in the t-th period of the hydroelectric unit number of service; is the output power value of the k-th photovoltaic unit in the multi-energy power system in the t-th period, k=1,2,3...N _pv , where N _pv is the output power of the photovoltaic unit in the t-th period of the multi-energy power system number of service; is the output power value of the g-th wind turbine in the multi-energy power system in the t-th period, g=1, 2, 3...N _w , N _w is the multi-energy power system in the t-th period of the wind turbine number of times; t=1, 2, 3...T, T is the number of time periods, Δt is the time step,

② Establish a complementary mechanism for multi-energy complementary and coordinated power generation in a multi-energy power system: based on the natural complementary characteristics of various power sources, relying on the good regulation capabilities of thermal power, gas power, and adjustable hydropower to stabilize wind power and photovoltaics. The fluctuation of the output power of resource-constrained power sources, making full use of clean renewable energy sources, reducing the proportion of thermal power in the power system, realizing the rational allocation of power system resources, and finally achieving a real-time balance between the total power generation and total load of the system. To optimize the operation purpose,

2) Constructing complementary index system and mathematical model of complementary demand respectively

①Build a complementary index system

The complementarity index is defined as the quantitative index of the complementary effect pursued by the multi-energy power system, that is, the direction of its optimization. Combined with the complementary mechanism of multi-energy complementary and coordinated power generation in the multi-energy power system, it can improve the system's renewable energy consumption capacity and energy saving. From the perspective of efficiency enhancement, build a complementary index system,

A. Calculate the penetration rate index of renewable energy in the multi-energy power system

Define r _pe (t) as the ratio of the output power of the renewable energy source to the load in the t-th period of the multi-energy power system, which mainly reflects the status of the renewable energy source in the power system and the penetration rate of the output power of the renewable energy source. The larger the value is, the cleaner the multi-energy power system is. Formula (2) is used to calculate the penetration rate of renewable energy in the multi-energy power system.

In the formula, is the penetration rate index of renewable energy in the multi-energy power system, represents the output power value of the renewable energy power supply in the multi-energy power system in the t-th period, P _l ^t represents the load value of the multi-energy power system in the t-th period, T is the number of periods, Δt is the time step, t = 1,2,3…T,

B. Calculate the coal consumption index of thermal power

In a multi-energy power system, the smaller the coal consumption of thermal power, the better the economy and environmental protection of thermal power. Formula (3) is used to calculate the coal consumption index of thermal power units,

In the formula, a _i , bi , c _i are the three fuel consumption characteristic coefficients of the _ith thermal power unit, i=1, 2, 3...N _th , N _th is the thermal power in the t-th time period in the multi-energy power system The number of units in service, f _it is the coal consumption index of the i-th thermal power unit in the t-th period, P _th.it is the output power value of the i-th thermal power unit in the t-th period, t=1,2, 3...T, Δt is the time step, T is the number of time periods,

C. Calculate the fluctuation range index of the load borne by thermal power

In order to reduce the frequent changes of output power of thermal power units, reduce coal consumption and improve utilization efficiency, the fluctuation range of load borne by thermal power should be minimized. The standard deviation of load value is used to represent the fluctuation range of load, and formula (4) is used to calculate thermal power The fluctuation range index of the unit's load,

In the formula, δ _th is the fluctuation range index of the load borne by thermal power units in the multi-energy power system, t = 1, 2, 3...T, Δt is the time step, T is the number of time periods, is the load value borne by all thermal power in the multi-energy power system in the t-th period, is the average load of all thermal power in the multi-energy power system during T periods,

D. Use formula (5) to calculate the power generation of hydropower,

In the formula, W _H is the total power generation of all hydropower in the multi-energy power system in T periods; is the water consumption of the jth hydroelectric unit for power generation in the tth period, is the head height of the j-th hydroelectric unit in the t-th period, j=1,2,3...N _hy , N _hy is the number of in-service hydro-electric units in the multi-energy power system in the t-th period, η is the hydroelectric conversion efficiency ;t=1,2,3...T, Δt is the time step, T is the number of time periods,

E. Use formula (6) to calculate the water abandonment index of hydropower,

In the formula, ΔQ is the index of water abandonment of all hydropower in T periods, Q _jmax is the maximum allocated water consumption of the jth hydropower unit in T periods, is the water consumption of the j-th hydroelectric unit in the t-th period, t=1, 2, 3...T, Δt is the time step, T is the number of periods, j=1, 2, 3...N _hy , N _hy is The number of in-service hydropower units in the multi-energy power system in the t-th period,

②Construct a mathematical model describing the complementary requirements of multiple heterogeneous power sources

The complementary demand of multi-heterogeneous power supplies is defined as the degree of matching between the output powers of the various heterogeneous power supplies after complementing each other with the load within a certain period of time. The two elements in the complementary demand are the power supply and the load. quantified between the complementary demand index and the complementary demand index between the power source and the load,

A. Calculate the complementary demand indicators between various heterogeneous power sources

A1. Use formula (7) to calculate the rate of change of the output power of the power supply,

r _i ^t =(P _i ^t -P _i ^t-1 )/Δt, (7)

In the formula, r _i ^t is the rate of change of the output power of the i-th power supply in the t-1th period to the t-th period, P _i ^t is the output power value of the i-th power source in the t-th period, P _i ^{t -1} is the output power value of the i-th power supply in the t-1th period, i=1,2,3...n, n is the number of power supply types under investigation, t=1,2,3...T; Δt is the time step size, T is the total number of epochs,

A2. Use formula (8) to calculate the set of absolute values of the sum of the output power change rates of thermal power, hydropower, photovoltaic, and wind power in each period,

In the formula, S _s is the set of absolute values of the sum of the output power change rates of thermal power, hydropower, photovoltaic, and wind power in each time period of T, and β _t is the thermal power, hydropower, photovoltaic, and wind power in the t-1th time period. The absolute value of the sum of the output power change rates in the t-th period, is the output power change rate of thermal power from the t-1th period to the tth period, is the output power change rate of hydropower from the t-1th period to the tth period, is the rate of change of photovoltaic output power from the t-1th period to the tth period, is the output power change rate of wind power from the t-1th period to the tth period, t=1, 2, 3...T, Δt is the time step, T is the number of periods,

A3. Use formula (9) to calculate the complementary demand index between various heterogeneous power sources,

In the formula, D _ss is the complementary demand index between various heterogeneous power sources in T periods, and β _t is the output power change of thermal power, hydropower, photovoltaic, and wind power from the t-1th period to the tth period. The absolute value of the sum of the rates, t=1,2,3...T, Δt is the time step, T is the number of time periods,

The smaller the value of the complementary demand index between various heterogeneous power sources, the stronger the mutual support between the wind-solar, hydro-thermal power sources in the time scale under investigation, that is, the better the complementary effect. The weaker the mutual support, the

B. Calculate the complementary demand index between the power supply and the load

B1. Use formula (10) to calculate the relative rate of change of the output power of the power supply,

In the formula, is the relative rate of change of the total output power of all power supplies from the t-1th period to the tth period, is the output power change rate of the total output power of all power sources from the t-1th period to the tth period, P _sc is the installed capacity of the generator sets in service in all power sources, T is the number of periods, t=1,2,3 ...T, Δt is the time step,

B2. Use formula (11) to calculate the relative rate of change of the load,

In the formula, is the relative change rate of the system load from the t-1th period to the tth period, P _lmax.T is the maximum load value in the T period, is the change rate of the system load from the t-1th period to the tth period, t=1, 2, 3...T, Δt is the time step, T is the number of periods,

B3. Use formula (12) to calculate the set of absolute values of the sum of the relative change rates in each period between the output power of all power supplies and the system load,

In the formula, S _l is the set of absolute values of the relative rate of change in each period between the output power of all power supplies and the system load in T time periods, and α _t is the difference between the output power of the power supply and the system load at the t-th The absolute value of the sum of the relative rate of change from 1 period to the t period, is the relative rate of change of the total output power of all power supplies from the t-1th period to the tth period, is the relative change rate of the system load from the t-1th period to the tth period, t=1, 2, 3...T, Δt is the time step, T is the number of periods,

B4. Use formula (13) to calculate the complementary demand index between the power supply and the load,

In the formula, D _sl is the complementary demand index between the power supply and the load of the multi-energy power system in T time periods, α _t is the difference between the output power of the power supply and the system load from the t-1th period to the tth period. The absolute value of the relative change rate of the time period and the sum, T is the number of time periods, t=1, 2, 3...T, Δt is the time step,

The smaller the value of the complementary demand index between the power source and the load, the more similar the change trend of the power source and the load in the investigated time scale; otherwise, the more different the change trend of the power source and the load,

The better the complementarity between power sources and between power sources and loads, the smaller the complementary demand, that is, the closer the complementary demand index value is to zero;

3) Define renewable energy complementary power sources

The renewable energy power sources that can meet the complementary needs after complementation are aggregated into one power source, which is defined as Renewable Energy Complementary Power Supply (RECPS). The main purpose of aggregation is to reduce the power fluctuation brought by wind and solar power to the multi-energy power system. The basic rule of aggregation is to meet the complementary needs of the system under the complementary mechanism of multi-energy complementary and coordinated power generation in the multi-energy power system. Based on energy conservation and environmental protection, the renewable energy power sources are aggregated, and the hydropower capacity ratio is calculated to minimize the complementary demand. After the renewable energy complementary power source is formed, its output power can follow the load fluctuation. In a multi-energy power system , the renewable energy complementary power supply is regarded as a kind of power supply, and it operates optimally together with other conventional power supplies, and when the system load value remains unchanged, the output power of the renewable energy complementary power supply also remains unchanged;

4) Formulate a hierarchical optimization operation strategy for multi-heterogeneous energy power sources

The hierarchical optimization operation strategy of multi-heterogeneous energy power sources is based on the complementary mechanism of multi-energy complementary and coordinated power generation in the multi-energy power system, giving priority to full acceptance of wind power and photovoltaic power generation; making full use of adjustable hydropower to deal with these uncertain power sources such as wind and solar. The randomness, intermittency and anti-peak regulation characteristics brought about by the present invention aggregate hydropower, wind power, and photovoltaics into a renewable energy complementary power supply. The output power of the renewable energy complementary power supply is relatively stable and can follow the load fluctuation, which can improve the multi-energy power supply. The ability of the system to absorb the wind and the wind and reduce the adverse impact of the uncertainty of wind and wind resources on the stable operation of the system,

The selected optimization objective of the multi-energy power system is to optimize the complementary indexes of the multi-energy power system and realize the rational allocation of power system resources. Including complementary power optimization layer, residual hydropower optimization layer, thermal power optimization layer,

①Complementary power optimization layer

In the optimization operation of the multi-energy power system, the renewable energy complementary power supply is firstly optimized, and the renewable energy complementary power supply is obtained by the aggregation of wind power, photovoltaic and hydropower. Obtain the aggregate capacity ratio of matching hydropower and wind power in the renewable energy complementary power supply, and then determine the output power of wind power, photovoltaic, and matching hydropower in each time period. The main goal of the renewable energy complementary power optimization layer is to determine the wind power, wind power The aggregate proportion of

In the formula, D _sl is the complementary demand between the power source and the load of the multi-energy power system in T time periods, is the relative rate of change of the load from the t-1th period to the tth period, is the relative change rate of the output power of the renewable energy complementary power supply from the t-1th period to the tth period, t=1, 2, 3...T, Δt is the time step, T is the number of periods,

②Remaining hydropower optimization layer

The net load curve is obtained by subtracting the output power of the renewable energy complementary power supply from the load value of the multi-energy power system. Under the condition that the net load is stable, the residual hydropower is used to generate electricity with the goal of minimizing the amount of waste water, and is calculated by formula (15). The objective function of the minimum amount of water waste from hydropower,

In the formula, ΔQ is the discarded water volume of the remaining hydropower in T periods, Q _jmax is the maximum allocated water consumption of the jth hydropower unit in T periods, is the water consumption of the jth hydroelectric unit in the tth period, j=1,2,3...N _rhy , N _rhy is the total number of remaining hydroelectric units, t=1,2,3...T, Δt is the time step , T is the number of time periods,

③ Thermal power optimization layer

The residual load is obtained by subtracting the residual hydropower output power from the net load curve. Under the dual action of the renewable energy complementary power supply and residual hydropower, the fluctuation range of the residual load is small, and the output power of the thermal power unit is arranged with the goal of minimum coal consumption. When the output power of the thermal power unit reaches the minimum, and the total output power of various heterogeneous power sources is still greater than the load, part of the output power of the renewable energy power source needs to be discarded, and formula (16) is used to calculate the objective function of the minimum coal consumption for thermal power generation,

In the formula, F is the total coal consumption of the thermal power unit, u _it is the on-off coefficient of thermal power, which is 1 when starting up and 0 when shutting down, f _it is the coal consumption of the i-th thermal power unit in the t-th period, i=1 ,2,3…N _th , N _th is the total number of thermal power units, t=1,2,3…T, Δt is the time step, T is the number of time periods,

5) Determine the constraints

The optimal operation of the multi-energy power system needs to satisfy the constraint equations of equations (17) to (23),

①Determine power balance constraints

The power balance constraint is expressed by Equation (17),

In the formula, is the total thermal power output power in the t-th period, is the remaining hydropower output power in the t-th period, is the total hydropower output power in the t-th period, is the total output power of the hydroelectric unit in the renewable energy complementary power supply in the t-th period, is the output power of the renewable energy complementary power supply in the t-th period, P _l ^t is the load value of the multi-energy power system in the t-th period, t=1, 2, 3...T, Δt is the time step, T is the number of periods,

② Determine the active output power constraints of the unit

The active output power constraint of the unit is expressed by formula (18),

In the formula, is the output power value of the i-th thermal power unit in the multi-energy power system in the t-th period, i=1,2,3...N _th , N _th is the in-service thermal power unit in the multi-energy power system in the t-th period quantity; is the output power value of the j-th hydroelectric unit in the multi-energy power system in the t-th period, j=1,2,3...N _hy , N _hy is the multi-energy power system in the t-th period of the hydroelectric unit number of service; is the output power value of the k-th photovoltaic unit in the multi-energy power system in the t-th period, k=1,2,3...N _pv , where N _pv is the output power of the photovoltaic unit in the t-th period of the multi-energy power system number of service; is the output power value of the g-th wind turbine in the multi-energy power system in the t-th period, g=1, 2, 3...N _w , N _w is the multi-energy power system in the t-th period of the wind turbine P _max.i is the upper limit of the output power of the ith thermal power unit in the multi-energy power system, P _max.j is the upper limit of the output power of the j-th hydropower unit in the multi-energy power system, and P _max.k is the multi-energy power system. The upper limit of the output power of the kth photovoltaic unit in the power system, P _max.g is the upper limit of the output power of the gth wind turbine in the multi-energy power system; P _min.i is the output of the ith thermal power unit in the multi-energy power system. Power lower limit, P _min.j is the output power lower limit of the jth hydroelectric unit in the multi-energy power system, t=1, 2, 3...T, Δt is the time step, T is the number of time periods,

③ Determine the system spinning reserve capacity constraints

The system spinning reserve capacity constraint is expressed by formula (19),

where: is the spinning reserve capacity of the system in the t-th period, is the rotating reserve capacity of the i-th thermal power unit in the t-th period, i=1,2,3...N _th , where N _th is the number of thermal power units in service in the t-th period in the multi-energy power system; is the rotating reserve capacity of the j-th hydroelectric unit in the t-th period, j=1, 2, 3...N _hy , N _hy is the number of in-service hydro-electric units in the multi-energy power system in the t-th period; α is The demand coefficient of the system load forecast error to the spinning reserve, β is the demand coefficient of the wind power output power prediction error to the spinning reserve; γ is the demand coefficient of the photovoltaic output power prediction error to the spinning reserve, is the load value of the multi-energy power system in the t-th period; is the output power value of the k-th photovoltaic unit in the multi-energy power system in the t-th period, k=1,2,3...N _pv , where N _pv is the output power of the photovoltaic unit in the t-th period of the multi-energy power system number of service; is the output power value of the g-th wind turbine in the multi-energy power system in the t-th period, g=1, 2, 3...N _w , N _w is the multi-energy power system in the t-th period of the wind turbine number of times, t=1, 2, 3...T, Δt is the time step, T is the number of time periods,

④ Determine the up-slope constraint and down-slope constraint of the unit

The up-climbing constraints and down-climbing constraints of the unit are expressed by formulas (20)~(21),

P _i ^{t+1 -} P _i ^t ≤ΔP _i ^up (20)

P _i ^t -P _i ^t+1 ≤ΔP _i ^down (21)

In the formula, P _i ^t+1 is the output power of the i-th unit in the multi-energy power system in the t+1-th period, P _i ^t is the output power of the i-th unit in the multi-energy power system in the t-th period, ΔP _i ^up is the maximum value of the ith unit uphill in the multi-energy power system, ΔP _i ^down is the maximum value of the ith unit in the multi-energy power system downhill, i=1,2,3...N,N is the number of active units in the multi-energy power system, t=1, 2, 3...T, Δt is the time step, T is the number of time periods,

⑤ Determine the constraints of hydropower generation

The constraint of hydropower generation is expressed by formula (22),

where: is the output power value of the j-th hydroelectric unit in the multi-energy power system in the t-th period, j=1,2,3...N _hy , N _hy is the number of in-service hydro-electric units in the multi-energy power system in the t-th period , W _H is the total power generation of all hydropower units in the multi-energy power system in T time periods, t=1, 2, 3...T, Δt is the time step, T is the number of time periods,

⑥ Determine hydroelectric flow constraints

The hydroelectric flow constraint is expressed by Equation (23),

In the formula: Q _jmin is the minimum allocated water consumption of the jth hydropower unit in the T period, Q _jmax is the maximum allocated water consumption of the jth hydropower unit in the T period, is the water consumption of the j-th hydropower unit in the t-th period, j=1,2,3...N _hy , N _hy is the number of in-service hydro-power units of the multi-energy power system in the t-th period, t=1,2 ,3…T, Δt is the time step, T is the number of time periods,

6) Through steps 1) to 5), using the constructed complementary index system and the mathematical model of complementary demand, combined with the actual operation constraints of the multi-energy power system, the output power of the renewable energy complementary power source is most matched with the load, and the remaining The objective function of the minimum water waste of hydropower and the minimum coal consumption of thermal power is solved by particle swarm algorithm, and finally the output power value of the multi-heterogeneous energy power source in each time period can be calculated to optimize the complementary index.

Combined with the multi-energy power system optimization operation model, the particle swarm algorithm program is written. The parameters in the algorithm program are set as: the particle swarm size is 20, the number of iterations is 500 times, the particle motion speed range is [-10, 10], and the learning factor is set. is 2, the linear inertia weight is used, the maximum value is set to 0.9, and the minimum value is set to 0.4.

The results of the optimization operation verification using the multi-energy power system optimization operation method based on the analysis of the complementary characteristics of the power supply of the present invention show that the acceptance capacity of the power grid for the wind and solar grid connection is enhanced, the grid connection ratio of the renewable energy is increased, and the natural resource constraints are reduced. Based on the influence of uncontrollable power generation on the power grid, the proportion of renewable energy generation increased from 72.76% before optimization to 76.57%. It can also be seen from the optimized operation results that the mean square error of thermal power output power is reduced from 196.775MW to 22.56MW, which reduces the fluctuation of thermal power output power and reduces the start, stop and adjustment of thermal power.

The specific embodiment of the present invention has made a detailed description of the content of the present invention, but it is not limited to this embodiment, and any obvious changes made by those skilled in the art according to the inspiration of the present invention belong to the scope of the right protection of the present invention .

Claims (1)

1. a kind of multiple-energy-source electric power system optimization operation method based on the analysis of power supply complementary characteristic, which is characterized in that it includes Following steps:

1) it establishes to provide multiple forms of energy to complement each other in multiple-energy-source electric system and coordinates the complementary mechanisms of power generation

1. being to provide multiple forms of energy to complement each other comprising the heterogeneous energy power supply that a variety of characteristics of output power are different in multiple-energy-source electric system Coordinate the precondition of power generation, wherein the complementary characteristic between each power supply has the characteristics that multiple-energy-source, multi-space, various dimensions, electricity Source complementary characteristic, will with formula (1) for being able to carry out mutually Ji mutual assistance power generation feature between different mass-energy source current Power supply complementary characteristic is expressed as the characteristic that various output power of power supply meet system loading,

In formula,For multiple-energy-source electric system t-th of period load value；For i-th of thermoelectricity in multiple-energy-source electric system Output power value of the unit t-th of period, i=1,2,3 ... N_th, N_thFor in multiple-energy-source electric system in t-th of period thermoelectricity The in-service quantity of unit；For output power value of j-th of Hydropower Unit t-th of period, j=in multiple-energy-source electric system 1,2,3...N_hy, N_hyFor multiple-energy-source electric system t-th of period Hydropower Unit in-service quantity；For multiple-energy-source power train Output power value of k-th of photovoltaic unit t-th of period in system, k=1,2,3...N_pv, N_pvExist for multiple-energy-source electric system The in-service quantity of t-th of period photovoltaic unit；It is g-th of Wind turbines in multiple-energy-source electric system defeated in t-th period Performance number out, g=1,2,3...N_w, N_wFor multiple-energy-source electric system t-th of period Wind turbines in-service quantity；T=1,2, 3 ... T, number of segment when T is, Δ t are time step,

Coordinate the complementary mechanisms of power generation 2. establishing and providing multiple forms of energy to complement each other in multiple-energy-source electric system are as follows: special with the complementation that various power supplys are natural Based on property, these natures of wind-powered electricity generation, photovoltaic are stabilized by the well-tuned ability that thermoelectricity, pneumoelectric, controllability water power have The fluctuation of resource limitation area output power of power supply makes full use of clean renewable energy source current, reduces thermoelectricity in power train Proportion in system realizes the reasonable disposition of power system resource, is finally reached system and always generates electricity and puts down in real time between total load The optimization of weighing apparatus runs purpose,

2) mathematical model of complementary index system, complementary demand is constructed respectively

1. constructing complementary index system

The quantizating index for the complementary effect that complementary index is pursued by multiple-energy-source electric system, the i.e. direction of its optimization are defined, In conjunction with the complementary mechanisms for coordinating power generation of providing multiple forms of energy to complement each other in multiple-energy-source electric system, from lifting system renewable energy digestion capability and From the point of view of energy efficiency, complementary index system is constructed,

A, renewable energy penetrance index in multiple-energy-source electric system is calculated

Define r_pe(t) ratio for multiple-energy-source electric system in t-th period renewable energy output power of power supply and load, master Reflect the status of renewable energy source current in the power system, the penetrance of renewable energy output power of power supply is bigger, card Bright multiple-energy-source electric system is more cleaned, and calculates renewable energy penetrance index in multiple-energy-source electric system with formula (2),

In formula,For renewable energy penetrance index in multiple-energy-source electric system,Indicating can be again in multiple-energy-source electric system Output power value of the raw energy power supply t-th of period, P_l ^tIndicate that load value of the multiple-energy-source electric system t-th of period, T are When number of segment, Δ t be time step, t=1,2,3 ... T,

B, the coal consumption figureofmerit of thermoelectricity is calculated

In multiple-energy-source electric system, the smaller economy for illustrating thermoelectricity of the coal consumption amount of thermoelectricity and the feature of environmental protection are better, with formula (3) The coal consumption figureofmerit of fired power generating unit is calculated,

In formula, a_i、b_i、c_iFor three fuel consumption characteristic coefficients of i-th of fired power generating unit, i=1,2,3 ... N_th, N_thFor multiple-energy-source In the in-service quantity of t-th of period fired power generating unit, f in electric system_i.tFor i-th of fired power generating unit t-th of period coal consumption amount Index, P_th.i.tOutput power value for i-th of fired power generating unit t-th of period, t=1,2,3 ... T, Δ t are time step, T For when number of segment,

C, the fluctuating range index of the undertaken load of thermoelectricity is calculated

In order to reduce fired power generating unit output power frequent change, reduce coal consumption amount, improve utilization efficiency, fire should be reduced to the greatest extent The fluctuating range of the undertaken load of electricity, the fluctuating range of load is indicated using the standard deviation of load value, calculates fire with formula (4) Motor group undertakes the fluctuating range index of load,

In formula, δ_thThe fluctuating range index of load, t=1,2,3 ... T, Δ t are undertaken by fired power generating unit in multiple-energy-source electric system For time step, number of segment when T is,For the load value that thermoelectricitys whole in multiple-energy-source electric system undertake t-th of period,The average value of load is undertaken within T period for thermoelectricitys whole in multiple-energy-source electric system,

D, the generated energy of water power is calculated with formula (5),

In formula, W_HFor within T period in multiple-energy-source electric system whole water power total power generation；For j-th of Hydropower Unit T-th of period power generation water consumption,Head height for j-th of Hydropower Unit t-th of period, j=1,2,3 ... N_hy, N_hyIt is multiple-energy-source electric system in the in-service quantity of t-th of period Hydropower Unit, η is water power transfer efficiency；T=1,2,3 ... T, Δ t be time step, number of segment when T is,

E, the abandoning water index of water power is calculated with formula (6),

In formula, Δ Q is abandoning water index of whole water power within T period, Q_jmaxIt is j-th of Hydropower Unit within T period Maximum allocated water consumption,Water consumption for j-th of Hydropower Unit t-th of period, t=1,2,3 ... T, Δ t are time step It is long, number of segment when T is, j=1,2,3 ... N_hy, N_hyFor multiple-energy-source electric system t-th of period Hydropower Unit in-service quantity,

2. the mathematical model of building description polyisocyanate matter power supply complementarity demand

Define polyisocyanate matter power supply complementarity demand be within a certain period of time, the output power between each heterogeneous power supply be complementary to one another after with The matching degree of load, two elements in complementary demand are power supply and load, pass through the complementation between various heterogeneous power supplys Complementary demand parameter between property demand parameter and power supply and load quantifies,

A, the complementary demand parameter between various heterogeneous power supplys is calculated

A1, the change rate that output power of power supply is calculated with formula (7),

In formula,It is i-th kind of power supply in the t-1 period to the output power change rate of t-th of period, P_i ^tExist for i-th kind of power supply The output power value of t-th of period, P_i ^t-1For the i-th kind of output power value of power supply in the t-1 period, i=1,2,3 ... n, n are The power type number investigated, t=1,2,3 ... T；Δ t be time step, number of segment when T is total,

A2, calculated with formula (8) thermoelectricity, water power, photovoltaic, the sum of wind-powered electricity generation day part output power change rate absolute value collection It closes,

In formula, S_sFor the absolute value of the sum of thermoelectricity, water power, photovoltaic, wind-powered electricity generation day part output power change rate within T period Set, β_tIt is thermoelectricity, water power, photovoltaic, wind-powered electricity generation in the t-1 period to the absolute of the sum of t-th of period output power change rate Value,It is thermoelectricity in the t-1 period to the output power change rate of t-th of period,It is water power in the t-1 period to The output power change rate of t period,It is photovoltaic in the t-1 period to the output power change rate of t-th of period,For Wind-powered electricity generation is in the t-1 period to the output power change rate of t-th of period, and t=1,2,3 ... T, Δ t is time step, when T is Number of segment,

A3, the complementary demand parameter between various heterogeneous power supplys is calculated with formula (9),

In formula, D_ssFor the complementary demand parameter between heterogeneous power supplys various in T period, β_tFor thermoelectricity, water power, photovoltaic, Wind-powered electricity generation is in the t-1 period to the absolute value of the sum of the output power change rate of t-th of period, t=1,2,3 ... T, when Δ t is Between step-length, number of segment when T is,

The value of complementary demand parameter between various heterogeneous power supplys is smaller, shows the honourable extreme misery in the time scale investigated Mutual supporting role between power supply is stronger, i.e., complementary effect is better, otherwise the mutual supporting role between honourable extreme misery power supply It is weaker,

B, the complementary demand parameter between power supply and load is calculated

B1, the relative change rate that output power of power supply is calculated with formula (10),

In formula,Relative change rate for all power supply gross outputs in the t-1 period to t-th of period,For all electricity Source gross output is in the t-1 period to the output power change rate of t-th of period, P_scFor in-service generator in all power supplys The installed capacity of group, number of segment when T is, t=1,2,3 ... T, Δ t are time step,

B2, with the relative change rate of formula (11) calculated load,

In formula,For in the relative change rate of the t-1 period to t-th of period system loading, P_lmax.TFor in T period most Big load value,For in the change rate of the t-1 period to t-th of period system loading, t=1,2,3 ... T, Δ t is the time Step-length, number of segment when T is,

B3, calculated with formula (12) relative change rate of day part between the output power and system loading of all power supplys and it is exhausted To the set of value,

In formula, S_lFor the relative change rate of day part between the output power and system loading of power supplys all in T period and The set of absolute value, α_tBecome in the t-1 period to t-th of the opposite of period between the output power and system loading of power supply The absolute value of rate sum,The relative change rate of gross output for all power supplys in the t-1 period to t-th of period,Relative change rate for system loading in the t-1 period to t-th of period, t=1,2,3 ... T, Δ t are time step, T For when number of segment,

B4, the complementary demand parameter between power supply and load is calculated with formula (13),

In formula, D_slFor the complementary demand parameter between the power supply and load of multiple-energy-source electric system in T period, α_tFor electricity Between the output power and system loading in source the relative change rate of the t-1 period to t-th period and absolute value, T is When number of segment, t=1,2,3 ... T, Δ t be time step,

The value of complementary demand parameter between power supply and load is smaller, shows power supply and load in the time scale investigated Variation tendency is more close；Otherwise, the variation tendency of power supply and load is more different,

Complementarity between power supply, between power supply and load is better, and complementary demand is smaller, i.e., complementary demand parameter value more connects It is bordering on zero；

3) renewable energy complementary power supply is defined

The renewable energy source current that can satisfy complementary demand after complementation is polymerized to a kind of power supply, is defined as renewable energy Complementary power supply (Renewable Energy Complementary Power Supply, RECPS), renewable energy complementary electrical Source is polymerized by whole wind-solar power supplies and required matching water power, and main purpose is to reduce wind-solar power supply and give multiple-energy-source electric power System bring power swing, the primitive rule of polymerization are the mutual assisting vehicle provided multiple forms of energy to complement each other in multiple-energy-source electric system and coordinate power generation Under system, meet system complementarity demand, using energy conservation and environmental protection as criterion, renewable energy source current is polymerize, calculates with complementation Property the minimum target of demand water power capacity ratio, formed renewable energy complementary power supply after, output power can follow load Fluctuation, in multiple-energy-source electric system, renewable energy complementary power supply sees a kind of power supply as, together with other normal power supplies Operation is optimized, and when system load value is constant, the output power of renewable energy complementary power supply is also remained unchanged；

4) the hierarchy optimization operation reserve of polyisocyanate mass-energy source current is formulated

The hierarchy optimization operation reserve of polyisocyanate mass-energy source current coordinates the complementation of power generation to provide multiple forms of energy to complement each other in multiple-energy-source electric system Based on mechanism, preferential receiving wind-powered electricity generation in full and photovoltaic power generation；Make full use of these not true property electricity of adjustable water power reply scene Water power, wind-powered electricity generation, photovoltaic are aggregated into renewable energy by randomness brought by source and intermittent and anti-tune peak character, the present invention Complementary power supply, renewable energy complementary power supply output power is relatively stable and can follow load fluctuation, can be improved multiple-energy-source Electric system gives system stable operation bring to adversely affect the digestion capability of scene, the honourable resource uncertainty of reduction,

The optimization aim of selected multiple-energy-source electric system is that the complementary index of multiple-energy-source electric system is allowed to be optimal simultaneously The reasonable disposition for realizing power system resource runs plan according to the hierarchy optimization that complementary index formulates polyisocyanate mass-energy source current Slightly, wherein comprising complementary power supply optimization layer, remaining optimization of hydroelectric generation layer, thermoelectricity optimization layer,

1. complementary power supply optimization layer

Multiple-energy-source electric power system optimization operation in renewable energy complementary power supply is optimized first, using wind-powered electricity generation, photovoltaic, Water power polymerize to obtain renewable energy complementary power supply, with the minimum target of complementary demand, the output work based on wind-powered electricity generation and photovoltaic Rate predicted value obtains the aggregate capacity of required matching water power and scene in renewable energy complementary power supply and matches, and then determines each The output power of wind-powered electricity generation, photovoltaic, required matching water power in period, the main target of renewable energy complementary power supply optimization layer is true The polymerization ratio of fixed scene water calculates the smallest objective function of complementary demand with formula (14),

In formula, D_slFor complementary demand of the multiple-energy-source electric system between the power supply and load in T period,For load In the relative change rate of the t-1 period to t-th of period,It is renewable energy complementary power supply output power in t-1 Relative change rate of a period to t-th of period, t=1,2,3 ... T, Δ t be time step, number of segment when T is,

2. remaining optimization of hydroelectric generation layer

Net load curve is obtained after the load value of multiple-energy-source electric system is subtracted renewable energy complementary power supply output power, It under the conditions of guaranteeing that net load is stable, is generated electricity with abandoning the minimum target exploitation residue water power of water, is calculated with formula (15) Water power abandons the smallest objective function of water,

In formula, Δ Q is abandoning water of the remaining water power within T period, Q_jmaxFor maximum of j-th of Hydropower Unit within T period Water consumption is distributed,Water consumption for j-th of Hydropower Unit t-th of period, j=1,2,3 ... N_rhy, N_rhyFor remaining hydroelectric machine Group total quantity, t=1,2,3 ... T, Δ t be time step, number of segment when T is,

3. thermoelectricity optimization layer

Remaining water power output power is subtracted on net load curve and obtains remaining load, in renewable energy complementary power supply and residue Under water power double action, the fluctuating range of remaining load is smaller, and the output work of fired power generating unit is arranged with the minimum target of coal consumption amount Rate when the total output power of various heterogeneous power supplys is still greater than load, then needs to give up when the output power of fired power generating unit reaches minimum Part renewable energy output power of power supply calculates the smallest objective function of thermoelectricity coal consumption amount with formula (16),

In formula, F is fired power generating unit total consumption of coal amount, u_i.tFor thermoelectricity startup-shutdown coefficient, the duration that is switched on is 1, and when shutdown is 0, f_i.tFor Coal consumption amount of i-th of fired power generating unit t-th of period, i=1,2,3 ... N_th, N_thFor fired power generating unit total quantity, t=1,2,3 ... T, Δ t be time step, number of segment when T is,

5) constraint condition is determined

The optimization operation of multiple-energy-source electric system needs to meet formula (17)~formula (23) constraint equation,

1. determining power-balance constraint

Power-balance constraint formula (17) expression,

In formula,For whole thermoelectricity output powers of t-th of period,For t-th of period residue water power output power, For t-th of period whole water power output power,For in t-th of period renewable energy complementary power supply Hydropower Unit gross output,For t-th of period renewable energy complementary power supply output power, P_l ^tFor multiple-energy-source electric power The load value of t-th of period of system, t=1,2,3 ... T, Δ t be time step, number of segment when T is,

2. determining the active output power constraint of unit

The active output power constraint of unit is indicated with formula (18),

In formula,For output power value of i-th of fired power generating unit t-th of period in multiple-energy-source electric system, i=1,2,3 ... N_th, N_thFor in multiple-energy-source electric system in the in-service quantity of t-th of period fired power generating unit；It is in multiple-energy-source electric system Output power value of the j Hydropower Unit t-th of period, j=1,2,3...N_hy, N_hyIt is multiple-energy-source electric system at t-th The in-service quantity of section Hydropower Unit；It is k-th of photovoltaic unit in multiple-energy-source electric system in the output power of t-th of period Value, k=1,2,3...N_pv, N_pvFor multiple-energy-source electric system t-th of period photovoltaic unit in-service quantity；For multiple-energy-source Output power value of g-th of Wind turbines t-th of period in electric system, g=1,2,3...N_w, N_wFor multiple-energy-source power train It unites in the in-service quantity of t-th of period Wind turbines；P_max.iFor the output power of i-th of fired power generating unit in multiple-energy-source electric system The upper limit, P_max.jFor the output power upper limit of j-th of Hydropower Unit in multiple-energy-source electric system, P_max.kFor multiple-energy-source electric system In k-th of photovoltaic unit the output power upper limit, P_max.gOn output power for g-th of Wind turbines in multiple-energy-source electric system Limit；P_min.iFor the output power lower limit of i-th of fired power generating unit in multiple-energy-source electric system, P_min.jFor in multiple-energy-source electric system The output power lower limit of j-th of Hydropower Unit, t=1,2,3 ... T, Δ t be time step, number of segment when T is,

3. determining the constraint of system spinning reserve capacity

System spinning reserve capacity constraint formula (19) expression,

In formula:Spinning reserve capacity for system t-th of period,For i-th of fired power generating unit t-th of period rotation Turn spare capacity, i=1,2,3 ... N_th, N_thFor in multiple-energy-source electric system in the in-service quantity of t-th of period fired power generating unit；Spinning reserve capacity for j-th of Hydropower Unit t-th of period, j=1,2,3...N_hy, N_hyFor multiple-energy-source electric system In the in-service quantity of t-th of period Hydropower Unit；α is that system loading predicts error to the service demand factor of spinning reserve, and β is wind-powered electricity generation Output power predicts error to the service demand factor of spinning reserve；γ is that photovoltaic output power predicts demand of the error to spinning reserve Coefficient,For multiple-energy-source electric system t-th of period load value；For k-th of photovoltaic machine in multiple-energy-source electric system Output power value of the group t-th of period, k=1,2,3...N_pv, N_pvIt is multiple-energy-source electric system in t-th of period photovoltaic machine The in-service quantity of group；For output power value of g-th of Wind turbines t-th of period in multiple-energy-source electric system, g=1, 2,3...N_w, N_wIn-service quantity for multiple-energy-source electric system in t-th of period Wind turbines, t=1,2,3 ... T, Δ t are the time Step-length, number of segment when T is,

4. determining Climing constant, lower Climing constant on unit

Climing constant, lower Climing constant are indicated with formula (20)~(21) on unit,

P_i ^t+1-P_i ^t≤ΔP_i ^up (20)

P_i ^t-P_i ^t+1≤ΔP_i ^down (21)

In formula, P_i ^t+1For output power of i-th of unit in the t+1 period, P in multiple-energy-source electric system_i ^tFor multiple-energy-source electric power Output power of i-th of unit t-th of period in system, Δ P_i ^upTo climb most on i-th of unit in multiple-energy-source electric system Big value, Δ P_i ^downFor for maximum value of climbing under i-th of unit in multiple-energy-source electric system, i=1,2,3 ... N, N are multiple-energy-source electricity In-service unit number in Force system, t=1,2,3 ... T, Δ t be time step, number of segment when T is,

5. determining the constraint of water power generated energy

Water power generated energy constraint formula (22) expression,

In formula:For output power value of j-th of Hydropower Unit t-th of period in multiple-energy-source electric system, j=1,2,3 ... N_hy, N_hyIn-service quantity for multiple-energy-source electric system in t-th of period Hydropower Unit, W_HFor all water in multiple-energy-source electric system Total power generation of the motor group within T period, t=1,2,3 ... T, Δ t be time step, number of segment when T is,

6. determining water power traffic constraints

Formula (23) expression of water power traffic constraints,

In formula: Q_jminThe smallest allocation water consumption for being j-th of Hydropower Unit within T period, Q_jmaxIt is j-th of Hydropower Unit in T Maximum allocated water consumption in a period,Water consumption for j-th of Hydropower Unit t-th of period, j=1,2,3 ... N_hy, N_hyIn-service quantity for multiple-energy-source electric system in t-th of period Hydropower Unit, t=1,2,3 ... T, Δ t are time step, and T is When number of segment,

6) by 1)~5) step, using the mathematical model of constructed complementary index system and complementary demand, in conjunction with more The constraint of electricity power running most matches the output power of renewable energy complementary power supply with load, remaining water power Abandoning water is minimum, the smallest objective function of coal consumption amount of thermoelectricity, solved using particle swarm algorithm, finally calculating can Complementary index is set to be optimal corresponding polyisocyanate mass-energy source current in the output power value of day part.