CN108009941B - Nesting optimization method for solving water-light complementary power station unit combination problem - Google Patents

Nesting optimization method for solving water-light complementary power station unit combination problem Download PDF

Info

Publication number
CN108009941B
CN108009941B CN201711221609.1A CN201711221609A CN108009941B CN 108009941 B CN108009941 B CN 108009941B CN 201711221609 A CN201711221609 A CN 201711221609A CN 108009941 B CN108009941 B CN 108009941B
Authority
CN
China
Prior art keywords
water
unit
constraint
power station
units
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201711221609.1A
Other languages
Chinese (zh)
Other versions
CN108009941A (en
Inventor
明波
刘攀
高仕达
谢艾利
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN201711221609.1A priority Critical patent/CN108009941B/en
Publication of CN108009941A publication Critical patent/CN108009941A/en
Application granted granted Critical
Publication of CN108009941B publication Critical patent/CN108009941B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Development Economics (AREA)
  • Operations Research (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Engine Equipment That Uses Special Cycles (AREA)

Abstract

The invention provides a nesting optimization method for solving a water-light complementary power station unit combination problem, which is characterized by comprising the following steps of: the method comprises the following steps: predicting a photoelectric output process, generating a plurality of photoelectric output scenes by considering prediction uncertainty, and calculating the probability corresponding to each scene; step two: establishing a water-light complementary power station unit combination mathematical model, and determining a target function, a constraint condition and a decision variable of the model; the optimization target of the water-light complementary power station unit combination mathematical model is that the average water consumption of the hydroelectric generating set is minimum under multiple scenarios, and the calculation formula is as follows:
Figure DDA0001480110500000011
step three: and constructing a nested optimization method, wherein the outer layer adopts an intelligent algorithm to optimize the number of the starting units of the unit, and the inner layer adopts a dynamic planning method to determine the optimal load distribution strategy under the given number of the starting units of the unit.

Description

Nesting optimization method for solving water-light complementary power station unit combination problem
Technical Field
The invention belongs to the cross field of renewable energy utilization and reservoir scheduling, and particularly relates to a nesting optimization method for solving a water-light complementary power station unit combination problem.
Technical Field
The solar photovoltaic is a non-schedulable renewable energy source, the output of the solar photovoltaic is influenced by factors such as day and night alternation, weather change, cloud layer thickness and the like, and the solar photovoltaic has obvious intermittence, volatility and randomness. The large-scale direct grid connection of photoelectricity can bring great pressure to peak regulation and stable operation of the system. The method is characterized in that water-light complementary power generation is implemented, photoelectricity is transmitted to a hydropower station, the photoelectricity is compensated by utilizing the quick adjustment capacity of a hydroelectric generating set, and the photoelectricity and the hydropower are packed together to be networked, so that the method is an effective way for promoting photoelectricity absorption.
However, due to the strong randomness of the photoelectric, the current prediction method cannot accurately predict the photoelectric. When uncertain photoelectricity is accessed into the hydropower station, the hydropower scheduling decision becomes uncertain, and the difficulty of safe and economic operation of the hydropower station is further increased. The hydropower station unit combination problem is the problem to be solved firstly when a short-term power generation plan of a hydropower station is compiled. The reasonable starting and stopping scheme can further increase the economic benefit of the hydropower station or save the operation cost. Due to the complexity of the hydroelectric generating set combination problem, the optimization solution is very difficult. Especially for large hydropower stations (a large number of units and large installed capacity), when the existing optimization technology is adopted for solving, the unification of calculation precision and efficiency can not be realized: if a dynamic programming method is adopted, a global optimal solution of the problem can be obtained, but the method occupies a large number of calculation inner layers and cannot effectively process time period coupling constraints (such as minimum startup and shutdown constraints); the time interval coupling constraint can be processed by adopting an intelligent algorithm, but the time interval coupling constraint can be trapped into a local optimal solution when the optimization variables are more, and premature convergence is shown. Therefore, the solution for researching the water-light complementary power station unit combination problem not only can provide technical support for economic operation of the water-light complementary power station unit, but also can further enrich and expand the hydropower optimization scheduling theory.
Disclosure of Invention
The invention aims to solve the problem of safe and economic operation of a water-light complementary power station, and provides a nested optimization method for solving the problem of unit combination of the water-light complementary power station.
In order to achieve the purpose, the invention adopts the following scheme:
the invention provides a nesting optimization method for solving a water-light complementary power station unit combination problem, which is characterized by comprising the following steps of: the method comprises the following steps: predicting a photoelectric output process, generating a plurality of photoelectric output scenes by considering prediction uncertainty, and calculating the probability corresponding to each scene; step two: establishing a water-light complementary power station unit combination mathematical model, and determining a target function, a constraint condition and a decision variable of the model; the optimization target of the water-light complementary power station unit combination mathematical model is that the average water consumption of the hydroelectric generating set is minimum under multiple scenarios, and the calculation formula is as follows:
Figure BDA0001480110480000021
in the formula: f is the average water consumption of the water-light complementary power station in the whole scheduling period; n is the number of hydropower station units; t is the number of scheduling period; s is lightNumber of power-out scenarios; n, t and s are respectively the unit number, the scheduling time interval number and the photoelectric scene number; u. ofn,tThe on-off state (0-1 variable) of the unit is obtained;
Figure BDA0001480110480000022
the machine flow of the machine set is measured; delta t is the scheduling period length; step three: the method comprises the steps of constructing a nested optimization method, optimizing the number of the units for starting by adopting an intelligent algorithm (such as a genetic algorithm and a cuckoo algorithm) on the outer layer, and determining the optimal load distribution strategy by adopting a dynamic planning method on the inner layer under the condition of the given number of the units for starting.
The nesting optimization method for solving the water-light complementary power station unit combination problem provided by the invention can also have the following characteristics: in the first step: firstly, a mathematical statistical method or a physical method is adopted to predict the photoelectric output process (P) of the next daytT1, …, T); second, the predicted contribution process is subtracted by a different prediction error (e)1,e2,e3) A variety of scenarios can be generated; finally, it is assumed that the prediction error of the photo follows a normal distribution N (μ, σ)2) The discrete probability distribution is used to replace the continuous probability distribution, and the probability (rho) corresponding to each scene can be calculated123) The calculation formula is as follows:
Figure BDA0001480110480000031
the nesting optimization method for solving the water-light complementary power station unit combination problem provided by the invention can also have the following characteristics: in the second step:
Figure BDA0001480110480000032
Figure BDA0001480110480000033
in the formula: f. ofrphThe relation among the unit flow, output and water head in the power characteristic curve of the hydroelectric generating set; f. ofvzThe relation between water level and reservoir capacity is formed; f. ofqzThe relationship between the downward discharge flow and the tail water level is adopted;
Figure BDA0001480110480000034
the output of the nth unit in the s scene in the t time period;
Figure BDA0001480110480000035
Figure BDA0001480110480000036
the water purification head, the dam front water level, the tail water level and the head loss in the t time period in the s type of scenes are respectively;
Figure BDA0001480110480000037
and
Figure BDA0001480110480000038
the storage capacity of the reservoir is at the beginning of the t-th time period and at the end of the t-th time period in the s-th scene.
The nesting optimization method for solving the water-light complementary power station unit combination problem provided by the invention can also have the following characteristics: in the second step, the constraint conditions considered by the model are as follows: water quantity balance constraint, reservoir capacity constraint, flow constraint of passing machine, unit output constraint, electric power balance constraint, load standby constraint, output lifting constraint, minimum start-stop constraint and vibration area constraint,
Figure BDA0001480110480000039
Figure BDA00014801104800000310
Figure BDA00014801104800000311
Figure BDA00014801104800000312
Figure BDA00014801104800000313
in the formula: i istThe reservoir warehousing flow is the t-th time period;
Figure BDA00014801104800000314
for reservoir in the s th sceneThe total generating flow in the t period; WPt sThe water discharge rate of the reservoir in the t th time period in the s th scene is obtained;
Figure BDA00014801104800000315
and
Figure BDA00014801104800000316
respectively the lower limit and the upper limit of the storage capacity;
Figure BDA00014801104800000317
and
Figure BDA00014801104800000318
respectively is a lower limit value and an upper limit value of the generating flow of the nth unit in the t period;
Figure BDA00014801104800000319
and
Figure BDA00014801104800000320
respectively the lower limit and the upper limit of the unit output; pt sPhotoelectric force generation values for the t-th time period of the s-th scene; dtThe total load assigned to the system; LRtA load reserve value at the t-th time period of the hydropower station; Δ pdAnd Δ puRespectively the upper limit values of the output descending and ascending speeds of the hydropower station; SUnAnd SDnThe minimum duration time of the starting and stopping states of the hydroelectric generating set; sun,tIndicating the state for the starting process of the unit (1 is starting, and 0 is non-starting); sdn,tIndicating the state (1 is closed, and 0 is not closed) for the shutdown process of the unit;
Figure BDA0001480110480000041
and
Figure BDA0001480110480000042
respectively the lower limit and the upper limit of the vibration area of the unit.
The nesting optimization method for solving the water-light complementary power station unit combination problem provided by the invention can also have the following characteristics: in the third step, considering that most of the existing hydropower stations are composed of the same units, the determined on-off state of the units in the unit combination optimization problem is converted into the determined number of the units to be started, so that the dimension reduction is realized.
The nesting optimization method for solving the water-light complementary power station unit combination problem provided by the invention can also have the following characteristics: in the third step, when the outer layer adopts an intelligent algorithm to optimize the number of the units for starting up, in order to effectively process the minimum start-up and shut-down constraints in the unit combination model and simultaneously realize dimension reduction, the encoding mode of the intelligent algorithm solution is as follows:
Figure BDA0001480110480000043
m is more than or equal to 2 and less than T, wherein:
Figure BDA0001480110480000044
is the solution of the intelligent algorithm; m is the number of time segments for dividing the whole scheduling period according to the number of the online units; t is the number of time segments in the whole scheduling period; x is the number of1,x2,…,xm-1Respectively numbering the time interval before the number of the units started in the whole scheduling period is changed;
Figure BDA0001480110480000045
are respectively 1 to x1,x1+1~x2,…,xm-1The number of the units started in the period of +1 to T.
The nesting optimization method for solving the water-light complementary power station unit combination problem provided by the invention can also have the following characteristics: in the third step, when the inner layer adopts dynamic programming to perform optimal load distribution, the two-stage recursion equation is as follows:
Figure BDA0001480110480000046
in the formula:
Figure BDA0001480110480000047
as a total load of
Figure BDA0001480110480000048
Optimal water consumption distributed among the d units; f. ofrph(pd,t,ht) Is negativeThe lotus is pd,tHead of water htThe water consumption of the d-th unit is calculated;
Figure BDA0001480110480000049
as a total load of
Figure BDA00014801104800000410
Optimal water consumption is allocated among d-1 units.
The nesting optimization method for solving the water-light complementary power station unit combination problem provided by the invention can also have the following characteristics: before executing the nested optimization, a dynamic programming method is adopted in advance to calculate the optimal load distribution strategy of the unit, and the optimal load distribution strategy is stored in a database; and directly calling an optimization result when the nested optimization is executed so as to improve the optimization efficiency.
Action and Effect of the invention
The nesting optimization method provided by the invention has theoretical global optimality, the number of optimized variables is small, and the minimum startup and shutdown constraint in the unit combination model can be effectively processed. In addition, uncertainty of photoelectric output prediction is considered in the optimization model, so that the scheduling scheme formulated according to the invention can still realize safe and economic operation of the water-light complementary power station under the condition of inaccurate photoelectric prediction.
Drawings
FIG. 1 is a flowchart of a nested optimization method for solving a water-light complementary power station unit combination problem in an embodiment of the present invention;
fig. 2 is a schematic diagram of an encoding mode of an intelligent algorithm during outer layer optimization.
Detailed Description
The following describes in detail a specific embodiment of the nested optimization method for solving the problem of the water-photovoltaic complementary power plant unit combination according to the present invention with reference to the accompanying drawings.
< example >
As shown in fig. 1, the nested optimization method for solving the water-photovoltaic complementary power station unit combination problem provided by this embodiment includes the following steps:
1. predicting a photoelectric output process; considering the prediction uncertainty, generating a plurality of photoelectric output scenes, and calculating the probability corresponding to each scene;
firstly, a mathematical statistical method or a physical method is adopted to predict the photoelectric output process (P) of the next dayt,t=1,…,T);
Second, the predicted contribution process is subtracted by a different prediction error (e)1,e2,e3) A variety of scenarios may be generated.
Finally, it is assumed that the prediction error of the photo follows a normal distribution N (μ, σ)2) The discrete probability distribution is used to replace the continuous probability distribution, and the probability (rho) corresponding to each scene can be calculated123) The calculation formula is as follows:
Figure BDA0001480110480000061
2. and establishing a water-light complementary power station unit combination mathematical model, and determining a target function, a constraint condition and a decision variable of the model.
The optimization goals of the model are: the average water consumption of the hydroelectric generating set under multiple scenarios is minimum, and the calculation formula is as follows:
Figure BDA0001480110480000062
Figure BDA0001480110480000063
Figure BDA0001480110480000064
Figure BDA0001480110480000065
Figure BDA0001480110480000066
in the formula: f. ofrphFor the machine in the power characteristic curve of hydroelectric generating setThe relationship of the flow, the output and the water head of the cross machine is formed; f. ofvzThe relation between water level and reservoir capacity is formed; f. ofqzThe relationship between the downward discharge flow and the tail water level is adopted;
Figure BDA0001480110480000067
the output of the nth unit in the s scene in the t time period;
Figure BDA0001480110480000068
the water purification head, the dam front water level, the tail water level and the head loss in the t time period in the s type of scenes are respectively;
Figure BDA0001480110480000069
and
Figure BDA00014801104800000610
the storage capacity of the reservoir is at the beginning of the t-th time period and at the end of the t-th time period in the s-th scene.
The constraints considered by the model are: the system comprises a water quantity balance constraint, a storage capacity constraint, an over-machine flow constraint, a unit output constraint, an electric power balance constraint, a load standby constraint, an output lifting constraint, a minimum start-stop constraint and a vibration area constraint.
Figure BDA0001480110480000071
Figure BDA0001480110480000072
Figure BDA0001480110480000073
Figure BDA0001480110480000074
Figure BDA0001480110480000075
Figure BDA0001480110480000076
Figure BDA0001480110480000077
Figure BDA0001480110480000078
Figure BDA0001480110480000079
Figure BDA00014801104800000710
In the formula: i istThe reservoir warehousing flow is the t-th time period;
Figure BDA00014801104800000711
the total generating flow of the reservoir in the t-th scene is obtained; WPt sThe water discharge rate of the reservoir in the t th time period in the s th scene is obtained;
Figure BDA00014801104800000712
and
Figure BDA00014801104800000713
respectively the lower limit and the upper limit of the storage capacity;
Figure BDA00014801104800000714
and
Figure BDA00014801104800000715
respectively is a lower limit value and an upper limit value of the generating flow of the nth unit in the t period;
Figure BDA00014801104800000716
and
Figure BDA00014801104800000717
lower limit and upper limit of unit outputA limit value; pt sPhotoelectric force generation values for the t-th time period of the s-th scene; dtThe total load assigned to the system; LRtA load reserve value at the t-th time period of the hydropower station; Δ pdAnd Δ puRespectively the upper limit values of the output descending and ascending speeds of the hydropower station; SUnAnd SDnThe minimum duration time of the starting and stopping states of the hydroelectric generating set; sun,tIndicating the state for the starting process of the unit (1 is starting, and 0 is non-starting); sdn,tIndicating the state (1 is closed, and 0 is not closed) for the shutdown process of the unit;
Figure BDA00014801104800000718
and
Figure BDA00014801104800000719
respectively the lower limit and the upper limit of the vibration area of the unit.
The model inputs are: total load given by the system, reservoir inflow and different photoelectric prediction scenarios.
The variables to be optimized in the above model are: the on-off state of the unit, the output of the unit, the flow rate of the unit passing through the unit and the water head.
3. And constructing a nested optimization method, wherein the outer layer adopts an intelligent algorithm to optimize the number of the starting units of the unit, and the inner layer adopts a dynamic planning method to determine the optimal load distribution strategy under the given number of the starting units of the unit.
Considering that most of the existing hydropower stations are composed of the same units, the on-off state of the optimized units in the model can be further converted into the number of the started optimized units.
Due to the existence of the minimum startup and shutdown constraint, the number of the startup units of the unit is required to be kept unchanged in a plurality of continuous adjacent time intervals. As shown in fig. 2, optimizing the number of unit boots during the entire scheduling period may further be converted into optimizing the time when the number of unit boots changes and the number of unit boots. Therefore, the encoding method of the intelligent algorithm solution can be represented by the following formula:
Figure BDA0001480110480000081
the number of the sets started in the whole scheduling period can be obtained by the following decoding mode:
Figure BDA0001480110480000082
in the formula:
Figure BDA0001480110480000083
is an individual (solution) of the intelligent algorithm; y is a decision variable set in the whole scheduling period;
Figure BDA0001480110480000084
is the solution of the intelligent algorithm; m is the number of time segments for dividing the whole scheduling period according to the number of the online units; t is the number of time segments in the whole scheduling period; x is the number of1,x2,…,xm-1Respectively numbering the time interval before the number of the units started in the whole scheduling period is changed;
Figure BDA0001480110480000085
are respectively 1 to x1,x1+1~x2,…,xm-1The number of the units started in the period of +1 to T.
When the dynamic planning method is adopted in the inner layer to optimize the load distribution strategy, the calculation of the dynamic planning can be completed in advance due to the fact that the calculation time consumption of the dynamic planning is large. Namely, the optimal distribution strategy (the load born by each unit and the power generation flow) of all possible water head loads under different starting numbers is calculated. When the nested optimization is carried out, related calculation results are directly called. When the optimal load distribution is carried out by adopting dynamic planning, a two-stage recursion equation is as follows:
Figure BDA0001480110480000091
in the formula:
Figure BDA0001480110480000092
as a total load of
Figure BDA0001480110480000093
Optimal water consumption distributed among the d units; f. ofrph(pd,t,ht) As a load of pd,tHead of water htThe water consumption of the d-th unit is calculated;
Figure BDA0001480110480000094
as a total load of
Figure BDA0001480110480000095
Optimal water consumption is allocated among d-1 units.
The above embodiments are merely illustrative of the technical solutions of the present invention. The nested optimization method for solving the problem of the water-light complementary power station unit combination is not limited to the contents described in the above embodiments, but is subject to the scope defined by the claims. Any modification or supplement or equivalent replacement made by a person skilled in the art on the basis of this embodiment is within the scope of the invention as claimed in the claims.

Claims (5)

1. A nesting optimization method for solving a water-light complementary power station unit combination problem is characterized by comprising the following steps:
the method comprises the following steps: the photoelectric output process is predicted, the uncertainty of prediction is considered, a plurality of photoelectric output scenes are generated, and the probability corresponding to each scene is calculated:
firstly, a mathematical statistical method or a physical method is adopted to predict the photoelectric output process (P) of the next daytT1, …, T); second, the predicted contribution process is subtracted by a different prediction error (e)1,e2,e3) A variety of scenarios can be generated; finally, it is assumed that the prediction error of the photo follows a normal distribution N (μ, σ)2) The discrete probability distribution is used to replace the continuous probability distribution, and the probability (rho) corresponding to each scene can be calculated123) The calculation formula is as follows:
Figure FDA0002269074030000011
step two: establishing a water-light complementary power station unit combination mathematical model, and determining a target function, a constraint condition and a decision variable of the model;
the optimization target of the water-light complementary power station unit combination mathematical model is that the average water consumption of the hydroelectric generating set is minimum under multiple scenarios, and the calculation formula is as follows:
Figure FDA0002269074030000012
in the formula: f is the average water consumption of the water-light complementary power station in the whole scheduling period; n is the number of hydropower station units; t is the number of scheduling period; s is the number of photoelectric output scenes; n, t and s are respectively the unit number, the scheduling time interval number and the photoelectric scene number; u. ofn,tThe on-off state (0-1 variable) of the unit is obtained;
Figure FDA0002269074030000013
△ t is the dispatching time period length;
step three: constructing a nested optimization method, optimizing the number of the units for starting up by adopting an intelligent algorithm on the outer layer, determining the optimal load distribution strategy by adopting a dynamic programming method on the inner layer under the condition of the given number of the units for starting up,
when the outer layer adopts an intelligent algorithm to optimize the number of the units for starting up, in order to effectively process the minimum start-up and shut-down constraint in the unit combination model and simultaneously realize dimension reduction, the coding mode of the intelligent algorithm solution is as follows:
Figure FDA0002269074030000021
in the formula:
Figure FDA0002269074030000022
is the solution of the intelligent algorithm; m is the number of time segments for dividing the whole scheduling period according to the number of the online units; t is the number of time segments in the whole scheduling period; x is the number of1,x2,…,xm-1Respectively changing the number of the units in the whole scheduling periodChanging the number of the previous time interval;
Figure FDA0002269074030000023
are respectively 1 to x1,x1+1~x2,…,xm-1The number of the units started in the time interval of +1 to T,
when the inner layer adopts dynamic planning to carry out optimal load distribution, the two-stage recursion equation is as follows:
Figure FDA0002269074030000024
in the formula:
Figure FDA0002269074030000025
as a total load of
Figure FDA0002269074030000026
Optimal water consumption distributed among the d units; f. ofrph(pd,t,ht) As a load of pd,tHead of water htThe water consumption of the d-th unit is calculated;
Figure FDA0002269074030000027
as a total load of
Figure FDA0002269074030000028
Optimal water consumption is allocated among d-1 units.
2. The method for nested optimization for solving the problem of the combination of the water-light complementary power station set according to claim 1, is characterized in that:
wherein, in the step two:
Figure FDA0002269074030000029
Figure FDA00022690740300000210
Figure FDA00022690740300000211
Figure FDA00022690740300000212
in the formula: f. ofrphThe relation among the unit flow, output and water head in the power characteristic curve of the hydroelectric generating set; f. ofvzThe relation between water level and reservoir capacity is formed; f. ofqzThe relationship between the downward discharge flow and the tail water level is adopted;
Figure FDA00022690740300000213
the output of the nth unit in the s scene in the t time period;
Figure FDA00022690740300000214
the water purification head, the dam front water level, the tail water level and the head loss in the t time period in the s type of scenes are respectively;
Figure FDA00022690740300000215
and
Figure FDA00022690740300000216
the storage capacity of the reservoir is at the beginning of the t-th time period and at the end of the t-th time period in the s-th scene.
3. The method for nested optimization for solving the problem of the combination of the water-light complementary power station set according to claim 2, characterized in that:
in the second step, the constraint conditions considered by the model are as follows: water quantity balance constraint, reservoir capacity constraint, flow constraint of passing machine, unit output constraint, electric power balance constraint, load standby constraint, output lifting constraint, minimum start-stop constraint and vibration area constraint,
Figure FDA0002269074030000031
Figure FDA0002269074030000032
Figure FDA0002269074030000033
Figure FDA0002269074030000034
Figure FDA0002269074030000035
Figure FDA0002269074030000036
Figure FDA0002269074030000037
Figure FDA0002269074030000038
Figure FDA0002269074030000039
Figure FDA00022690740300000310
in the formula: i istThe reservoir warehousing flow is the t-th time period;
Figure FDA00022690740300000311
the total generating flow of the reservoir in the t-th scene is obtained; WPt sThe water discharge rate of the reservoir in the t th time period in the s th scene is obtained;
Figure FDA00022690740300000312
and
Figure FDA00022690740300000313
respectively the lower limit and the upper limit of the storage capacity;
Figure FDA00022690740300000314
and
Figure FDA00022690740300000315
respectively is a lower limit value and an upper limit value of the generating flow of the nth unit in the t period;
Figure FDA00022690740300000316
and
Figure FDA00022690740300000317
respectively the lower limit and the upper limit of the unit output; pt sPhotoelectric output value is obtained for the t-th time period in the s-th situation; dtThe total load assigned to the system; LRtLoad reserve value of hydropower station t period △ pdAnd △ puRespectively the upper limit values of the output descending and ascending speeds of the hydropower station; SUnAnd SDnThe minimum duration time of the starting and stopping states of the hydroelectric generating set; sun,tIndicating the state for the starting process of the unit (1 is starting, and 0 is non-starting); sdn,tIndicating the state (1 is closed, and 0 is not closed) for the shutdown process of the unit;
Figure FDA0002269074030000041
and
Figure FDA0002269074030000042
respectively the lower limit and the upper limit of the vibration area of the unit.
4. The method for nested optimization for solving the problem of the combination of the water-light complementary power station set according to claim 1, is characterized in that:
in the third step, the determined on-off state of the unit is converted into the determined number of the units for starting up, so that the dimension reduction is realized.
5. The method for nested optimization for solving the problem of the combination of the water-light complementary power station set according to claim 1, is characterized in that:
before executing nested optimization, a dynamic programming method is adopted in advance to calculate a unit optimal load distribution strategy, and the optimal load distribution strategy is stored in a database; the optimization results are directly invoked when performing nested optimizations.
CN201711221609.1A 2017-11-23 2017-11-23 Nesting optimization method for solving water-light complementary power station unit combination problem Active CN108009941B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711221609.1A CN108009941B (en) 2017-11-23 2017-11-23 Nesting optimization method for solving water-light complementary power station unit combination problem

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711221609.1A CN108009941B (en) 2017-11-23 2017-11-23 Nesting optimization method for solving water-light complementary power station unit combination problem

Publications (2)

Publication Number Publication Date
CN108009941A CN108009941A (en) 2018-05-08
CN108009941B true CN108009941B (en) 2020-05-26

Family

ID=62054625

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711221609.1A Active CN108009941B (en) 2017-11-23 2017-11-23 Nesting optimization method for solving water-light complementary power station unit combination problem

Country Status (1)

Country Link
CN (1) CN108009941B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112762576A (en) * 2020-12-29 2021-05-07 广东美的白色家电技术创新中心有限公司 Air conditioning system control method, temperature reaching time prediction model training method and equipment

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103745023A (en) * 2013-11-22 2014-04-23 华中科技大学 Coupling modeling method for hydropower station power generated output scheme making and optimal load distribution
CN103855718A (en) * 2014-03-10 2014-06-11 东南大学 Scheduling method for pumped storage power station to participate in electric power system with wind power
CN104239967A (en) * 2014-08-29 2014-12-24 华北电力大学 Multi-target economic dispatch method for power system with wind farm
CN104600713A (en) * 2014-12-25 2015-05-06 国家电网公司 Device and method for generating day-ahead reactive power dispatch of power distribution network containing wind/photovoltaic power generation
CN104795846A (en) * 2015-04-24 2015-07-22 清华大学 Optimized operation method of pumped-storage power station and wind power combined system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103745023A (en) * 2013-11-22 2014-04-23 华中科技大学 Coupling modeling method for hydropower station power generated output scheme making and optimal load distribution
CN103855718A (en) * 2014-03-10 2014-06-11 东南大学 Scheduling method for pumped storage power station to participate in electric power system with wind power
CN104239967A (en) * 2014-08-29 2014-12-24 华北电力大学 Multi-target economic dispatch method for power system with wind farm
CN104600713A (en) * 2014-12-25 2015-05-06 国家电网公司 Device and method for generating day-ahead reactive power dispatch of power distribution network containing wind/photovoltaic power generation
CN104795846A (en) * 2015-04-24 2015-07-22 清华大学 Optimized operation method of pumped-storage power station and wind power combined system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
水-光电联合运行短周期调度可行性分析;明波 等;《太阳能学报》;20151130;第36卷(第11期);正文第2731-2737页 *

Also Published As

Publication number Publication date
CN108009941A (en) 2018-05-08

Similar Documents

Publication Publication Date Title
CN105046395B (en) Method for compiling day-by-day rolling plan of power system containing multiple types of new energy
CN108009672B (en) Water-light complementary power station daily power generation planning method based on double-layer optimization model
CN103151803B (en) Method for optimizing wind power system-contained unit and backup configuration
CN111555281B (en) Method and device for simulating flexible resource allocation of power system
CN113193547B (en) Day-ahead-day cooperative scheduling method and system for power system considering uncertainty of new energy and load interval
CN103825280B (en) The autonomous on off control method and apparatus of pump-storage generator
CN104795846A (en) Optimized operation method of pumped-storage power station and wind power combined system
CN104299173B (en) It is a kind of to optimize dispatching method a few days ago suitable for the robust that various energy resources are accessed
CN107800153B (en) Electric heat energy rolling robust scheduling method for electric heat storage and wind power consumption
CN109347152B (en) Random production simulation method considering participation of multi-type power supply in peak shaving and application
CN106026184A (en) Pumped storage power station and wind power plant combination system for power grid peak regulation and optimized scheduling method thereof
CN112381375B (en) Rapid generation method for power grid economic operation domain based on tide distribution matrix
CN105305423A (en) Method for determining optimal error boundary considering intermittent energy uncertainty
CN109149631A (en) It is a kind of to consider that wind-light storage provides the two stages economic load dispatching method of flexible climbing capacity
CN111682531B (en) PL-IMOCS-based wind, light, water and fire primary energy complementary short-term optimization scheduling method and device
CN111934366A (en) Power grid multivariate optimization scheduling method for improving wind power receiving capacity
CN102593855B (en) Method for stabilizing fluctuation of output power of renewable energy power supply in power system
CN113659620A (en) Water-wind complementary power generation system day-ahead scheduling method based on dynamic frequency constraint
CN111130145B (en) Wind-solar unit assembly capacity optimization planning method based on wind and light discarding
CN108009941B (en) Nesting optimization method for solving water-light complementary power station unit combination problem
CN106505637B (en) Method for optimizing active scheduling conservative degree of power system
CN116205345A (en) Optimal scheduling method and system considering wind-solar prediction and maintenance plan
CN110729721A (en) Method for calculating global reserve capacity of power system
CN106655282B (en) Power system economic dispatching method realized by adopting 0-1 integer programming model
CN110854929B (en) Day-ahead scheduling method considering uncertainty in time period

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant