CN115423202A - Method and system for clearing cascade hydropower day-ahead spot goods based on time interval decoupling - Google Patents

Method and system for clearing cascade hydropower day-ahead spot goods based on time interval decoupling Download PDF

Info

Publication number
CN115423202A
CN115423202A CN202211139082.9A CN202211139082A CN115423202A CN 115423202 A CN115423202 A CN 115423202A CN 202211139082 A CN202211139082 A CN 202211139082A CN 115423202 A CN115423202 A CN 115423202A
Authority
CN
China
Prior art keywords
unit
time period
hydropower station
constraint
hydropower
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.)
Pending
Application number
CN202211139082.9A
Other languages
Chinese (zh)
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.)
State Grid Corp of China SGCC
State Grid Hunan Electric Power Co Ltd
Economic and Technological Research Institute of State Grid Hunan Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Hunan Electric Power Co Ltd
Economic and Technological Research Institute of State Grid Hunan Electric Power Co Ltd
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 State Grid Corp of China SGCC, State Grid Hunan Electric Power Co Ltd, Economic and Technological Research Institute of State Grid Hunan Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202211139082.9A priority Critical patent/CN115423202A/en
Publication of CN115423202A publication Critical patent/CN115423202A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION 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/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION 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/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315—Needs-based resource requirements planning or analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Development Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Educational Administration (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a method for clearing cascade hydropower station daily spot goods based on time interval decoupling, which comprises the steps of obtaining working data of a target cascade hydropower station system; constructing a stock clearing objective function of the cascade hydropower day ahead; establishing corresponding constraint conditions; and solving the objective function by adopting a variable-scale optimization method to finish the daily spot shipment clearing of the target cascade hydropower system. According to the method and the system for clearing the cascade hydropower daily spot goods based on the time-interval-decoupling, provided by the invention, multi-stage solution is carried out by utilizing variable-scale optimization, the constraint consideration of the original problem can be kept by a model by increasing the step length to weaken or even clear the time-interval-coupling constraint, the practicability of a clearing result is ensured, the optimization and solving space of an optimization model is promoted, the solving efficiency is promoted, and the method and the system are high in reliability, good in accuracy and good in practicability.

Description

Method and system for clearing cascade hydropower daily spot goods based on time interval decoupling
Technical Field
The invention belongs to the field of electrical automation, and particularly relates to a method and a system for clearing cascade hydropower daily spot goods based on time interval decoupling.
Background
With the development of economic technology and the improvement of living standard of people, electric energy becomes essential secondary energy in production and life of people, and brings endless convenience to production and life of people. Therefore, ensuring a stable and reliable supply of electric energy is one of the most important tasks of an electric power system.
At present, hydropower in China has been developed for a long time. In the electric power market of high-proportion hydropower, hydropower becomes an important power generation main body and has great influence on clearing results, and clearing algorithms are also influenced by the close space-time hydraulic power and electric power coupling relationship among cascade hydropower stations. Therefore, in the day-ahead clearance model with hydropower participating in the spot market, besides traditional time-interval coupling constraints such as unit climbing constraints, the time-interval coupling relationship of the model is further complicated by considering the relevant constraints of hydraulic power and electric power connection between upstream and downstream of the cascade hydropower, and how to weaken or remove the time-interval coupling relationship of the constraints so as to simplify the clearance model and improve the solution efficiency is a difficult problem to be solved urgently.
At present, few documents propose a method of decoupling a model into a plurality of submodels by relaxation constraint, and performing subsequent correction after parallel solution on all the submodels to perform time period decoupling. However, these methods are optimization models for which only a single period coupling constraint, the hill climbing constraint, exists. For the day-ahead emerging model considering the cascade hydropower constraint, because the time interval coupling relation is not limited to the adjacent time intervals any more, but exists in a plurality of time intervals before and after the delay time of the water flow of the upstream and downstream cascade power stations, and the output relation coupling of the upstream and downstream power stations is tight, the solving method of firstly relaxing the constraint and then adjusting the model afterwards is obviously not applicable any more. On the other hand, after the cascade hydroelectric relationship is subjected to linearization processing, a large number of variables 0 and 1 are inevitably added to a current stock-in-stock model, and the solution efficiency of a large-scale mixed integer programming model is relatively poor, so that the requirement of the current power grid is difficult to meet.
Disclosure of Invention
The invention aims to provide a method for clearing cascade hydropower station daily spot goods based on decoupling time period, which has high reliability, good accuracy and good practicability.
The invention also aims to provide a system for realizing the step hydropower station spot shipment method based on the decoupling of the time interval.
The invention provides a method for clearing cascade hydropower station spot goods before day based on time interval decoupling, which comprises the following steps:
s1, acquiring working data of a target step hydroelectric system;
s2, constructing a step hydropower station spot shipment clearing objective function in the day ahead according to the working data obtained in the step S1;
s3, establishing corresponding constraint conditions according to the objective function constructed in the step S2;
s4, solving the objective function constructed in the step S2 by adopting a variable scale optimization method under the constraint condition constructed in the step S3;
and S5, according to the solving result of the step S4, finishing the daily spot shipment of the target step hydropower system.
S2, constructing a step hydropower station spot shipment clearing objective function in the day ahead according to the working data acquired in the step S1, and specifically comprising the following steps:
the following equation is used as the objective function:
Figure BDA0003852692730000031
n is the total number of the units in the system; t is the optimized total time period number; p i,t The output of the unit i in the time period t is obtained; c i,t (P i,t ) The operation cost of the unit i in the time period t is calculated;
Figure BDA0003852692730000032
starting cost of the unit i in a time period t; NL is the total number of lines in the system; m is a group of 1 A penalty value for a network constraint;
Figure BDA0003852692730000033
a slack variable that is the upper limit of the transmission capacity of line l;
Figure BDA0003852692730000034
a slack variable that is the lower limit of the transmission capacity of line l; NS is the total number of sections in the system;
Figure BDA0003852692730000035
a relaxation variable which is the upper limit of the transmission capacity of the section s;
Figure BDA0003852692730000036
a relaxation variable which is a lower limit of a transmission capacity of the section s; omega C Is a set of hydropower stations in the system; m C A hydraulic power station hydraulic energy abandoning penalty factor is given;
Figure BDA0003852692730000037
the power of the hydropower station h is abandoned in the time period t.
Step S3, establishing a corresponding constraint condition according to the objective function established in step S2, specifically including the following steps:
the following equation is adopted as the system load balance constraint:
Figure BDA0003852692730000038
wherein N is the total number of the units in the system; p is i,t The output of the unit i in the time period t is obtained; NT is the total number of links in the system; t is a unit of j,t Planned power for tie j for time period t; d t System load for time period t;
the following equation is used as the system standby constraint:
Figure BDA0003852692730000039
Figure BDA00038526927300000310
in the formula of alpha i,t Is the start-stop state of the unit i in the time period t, and alpha i,t =1 indicates that the unit set i is in a start-up state in time t, α i,t =0 represents that the unit set i is in a shutdown state in the time period t;
Figure BDA00038526927300000311
the maximum output of the unit i in the time period t is obtained;
Figure BDA0003852692730000041
the minimum output of the unit i in the time period t is obtained;
Figure BDA0003852692730000042
the system positive spare capacity requirement for time period t;
Figure BDA0003852692730000043
the system negative spare capacity requirement of the time interval t;
the following formula is adopted as the output constraint of the unit:
Figure BDA0003852692730000044
in the formula P i,t The output of the unit i in the time period t is obtained;
the following formula is adopted as the unit climbing constraint:
Figure BDA0003852692730000045
Figure BDA0003852692730000046
in the formula, delta P i U The maximum upward slope climbing rate of the unit i is set; delta P i D The maximum downward climbing speed of the unit i;
the following formula is adopted as the minimum start-stop time constraint of the unit:
Figure BDA0003852692730000047
Figure BDA0003852692730000048
in the formula
Figure BDA0003852692730000049
For the time that the unit i has been continuously shut down during the time period t, an
Figure BDA00038526927300000410
Figure BDA00038526927300000411
Is the time that the unit i has been continuously started up in the time period t, an
Figure BDA00038526927300000412
T D Minimum continuous down time for the unit; t is U The minimum continuous starting time of the unit is set;
the following formula is adopted as the network security constraint:
Figure BDA00038526927300000413
Figure BDA00038526927300000414
in the formula P l min A lower limit for power flow transmission for line l; p is l max An upper limit for the power flow transmission of line l; g l-i Outputting a power transfer distribution factor for the generator of the line l by the node where the unit i is located; g l-j Outputting a power transfer distribution factor for the generator of the link line l by the node where the link line j is located; g l-k A generator output power transfer distribution factor for node k to line l; k is the total number of the network nodes; d k,t The load for node k at time period t;
Figure BDA0003852692730000051
is the lower limit of the power flow transmission of the section s;
Figure BDA0003852692730000052
is the upper limit of the power flow transmission of the section s; g s-i The generator output power of the section s is transferred to a distribution factor for the node where the unit i is located; g s-j The generator output power transfer distribution factor of the section s is calculated for the node where the tie line j is located; g s-k The distribution factor of the output power transfer of the generator of the node k to the section s is obtained;
the following equation is adopted as the determination constraint of the water disposal electric quantity:
Figure BDA0003852692730000053
in the formula
Figure BDA0003852692730000054
The hydraulic power of the hydropower station h in the time period t is obtained; a is h,t Is a variable from 0 to 1, and
Figure BDA0003852692730000055
when a h,t ＝1，
Figure BDA0003852692730000056
When a h,t ＝0；
Figure BDA0003852692730000057
The water discharge of the hydropower station h in the time period t; h h The water consumption rate of the hydropower station h; p h,max The maximum adjustable output of the hydropower station h; p is h,t The output of the hydropower station h in the time period t is obtained; m is a set maximum positive number and is used for carrying out linear transformation on the water abandonment electric quantity judgment constraint;
the following formula is adopted as the hydropower water level control constraint:
Figure BDA0003852692730000058
in the formula
Figure BDA0003852692730000059
Is the dead water level of the hydropower station h;
Figure BDA00038526927300000510
limiting the water level for the hydropower station h in the flood season; z h,0 The initial water level of the hydropower station h at zero point of the next day; t is the current time period; p is h,τ The output of the hydropower station h in the time period tau is obtained; h h Is the water consumption rate of the hydropower station h; s h The water surface area of the reservoir of the hydropower station h; f. of h,τ The amount of incoming water for hydropower station h in time period τ, and
Figure BDA00038526927300000511
I h,τ the natural inflow water flow of the hydropower station h in a time interval tau is shown as up (h), the upper hydropower station of the hydropower station h is shown as s (h), and the lag time of the upper hydropower station of the hydropower station h is shown as s (h);
the following formula is adopted as the constraint of the hydroelectric vibration region:
Figure BDA0003852692730000061
Figure BDA0003852692730000062
Figure BDA0003852692730000063
in the formula
Figure BDA0003852692730000064
The upper limit of the 1 st vibration area of the ith unit is set;
Figure BDA0003852692730000065
the lower limit of the 1 st vibration area of the ith unit; p i,1,t The output value of the operation interval is a continuous variable and represents the output value of the operation interval when the ith unit is in the 1 st operable interval in time period t; delta i,l,t Is a variable of 0 to 1 and indicates whether the ith unit is in the t period
Figure BDA0003852692730000066
And if the ith unit is in the t period
Figure BDA0003852692730000067
Delta between i,l,t =1, if the ith unit is not in t time interval
Figure BDA0003852692730000068
Is between delta i,l,t ＝0。
Step S4, solving the objective function constructed in step S2 by using the variable scale optimization method under the constraint condition constructed in step S3, specifically including the following steps:
A. determining a sequence of time interval steps;
B. b, establishing a plurality of optimization model sets with the same optimization period and control requirements but different step sizes according to the time period step length sequence selected in the step A; sequencing the optimization models in the sequence from large step to small step, and setting the serial number of the current optimization model to be solved as N = N Δ Corresponding time interval step is delta n (ii) a The optimization model comprises a cascade hydropower day-ahead spot shipment clearing objective function and a corresponding constraint condition;
C. for step size of Δ n Re-describing the time interval coupling type constraint in the optimization model:
equivalent conversion is carried out on the unit climbing constraint by adopting the following formula:
Figure BDA0003852692730000071
Figure BDA0003852692730000072
in the formula P i,t The output of the unit i in the time period t is obtained; delta P i U The maximum upward climbing speed is set i; delta 1 Step length of the original optimization model;
Figure BDA0003852692730000073
the maximum climbing upper limit between adjacent time periods of the current optimization model is set; delta P i D The maximum downward climbing speed is set for the unit i;
Figure BDA0003852692730000074
the maximum landslide upper limit between adjacent time periods of the current optimization model is set; if it is
Figure BDA0003852692730000075
Or
Figure BDA0003852692730000076
If the current optimization model is larger than the rated capacity of the unit, determining that the current optimization model has no corresponding constraint;
the hydropower water level control constraint is equivalently converted by adopting the following formula:
Figure BDA0003852692730000077
in the formula
Figure BDA0003852692730000078
Is the dead water level of the hydropower station h;
Figure BDA0003852692730000079
limiting the water level for the hydropower station h in the flood season; z h,0 The initial water level of the hydropower station h at the zero point of the next day; p h,τ The output of the hydropower station h in the time period tau is obtained; h h Is the water consumption rate of the hydropower station h;
Figure BDA00038526927300000710
the water discharge of the hydropower station h in the time period t; i is h,τ The natural incoming water flow of the hydropower station h in the time period tau is obtained; s h The water surface area of the reservoir of the hydropower station h;
Figure BDA00038526927300000711
converting the water flow delay time of the upstream power station of the hydropower station h in the current optimization model, and if so
Figure BDA00038526927300000712
Determining that the influence of the ex-warehouse flow of the upstream hydropower station on the downstream in the current optimization model is immediate, and not considering the influence of time periods;
D. if N = N Δ Then solving the corresponding optimization model to obtain an initial solution; otherwise, setting the initial output process of each set in the optimization model n as the last oneAs a result of the optimization process of stage n +1,
Figure BDA0003852692730000081
and is provided with
Figure BDA0003852692730000082
Calculating to obtain the corresponding outlet flow and water level of the hydropower station;
after determining the initial solution of part of time period in the optimization model n, solving again;
E. updating the value of n to be n-1; if n is larger than or equal to 1, skipping to the step C to continuously solve the next long optimization model; otherwise, ending the variable-scale optimization solving process to obtain a corresponding solving result.
The time interval step length sequence in the step A is specifically 6h, 3h, 1.5h and 15min, or 8h, 4h, 2h, 1h and 15min.
The invention also discloses a system for realizing the step hydropower station daily spot shipment method based on the decoupling of the time interval, which comprises a data acquisition module, an objective function construction module, a constraint condition construction module, a model solving module and a data output module; the data acquisition module, the objective function construction module, the constraint condition construction module, the model solving module and the data output module are sequentially connected in series; the data acquisition module is used for acquiring the working data of the target step hydropower system and uploading the working data to the target function construction module; the target function construction module is used for constructing a gradient hydropower station spot shipment and clearing target function in the day ahead according to the acquired working data and uploading the gradient hydropower station spot shipment and clearing target function to the constraint condition construction module; the constraint condition construction module is used for establishing corresponding constraint conditions according to the constructed objective function and uploading the constraint conditions to the model solving module; the model solving module is used for solving the objective function by adopting a variable-scale optimization method and uploading the result to the data output module; and the data output module is used for outputting the final current spot shipment and clearing result of the target step hydropower system.
According to the method and the system for clearing the cascade hydropower day-ahead spot goods based on time-interval-decoupling, provided by the invention, multi-stage solution is carried out by using variable-scale optimization, and the time-interval-coupling constraint is weakened by increasing the step length and even eliminated, so that the constraint consideration of the original problem can be kept for the model, the practicability of a clearing result is ensured, the optimization solution space of the optimization model is promoted, the solution efficiency is promoted, and the method and the system are high in reliability, good in accuracy and good in practicability.
Drawings
FIG. 1 is a schematic process flow diagram of the process of the present invention.
FIG. 2 is a functional block diagram of the system of the present invention.
Detailed Description
FIG. 1 is a schematic flow chart of the method of the present invention: the invention provides a method for clearing cascade hydropower station spot goods before day based on time interval decoupling, which comprises the following steps:
s1, acquiring working data of a target step hydropower system;
s2, constructing a step hydropower station spot shipment clearing objective function in the day ahead according to the working data acquired in the step S1; the method specifically comprises the following steps:
the following equation is used as the objective function:
Figure BDA0003852692730000091
n is the total number of the units in the system; t is the optimized total time period number; p is i,t The output of the unit i in the time period t is obtained; c i,t (P i,t ) The operation cost of the unit i in the time period t is calculated;
Figure BDA0003852692730000092
starting cost of the unit i in a time period t; NL is the total number of lines in the system; m is a group of 1 A penalty value for a network constraint;
Figure BDA0003852692730000093
a slack variable that is the upper limit of the transmission capacity of line l;
Figure BDA0003852692730000094
a slack variable which is the lower limit of the transmission capacity of the line l; NS is the total number of sections in the system;
Figure BDA0003852692730000095
a relaxation variable which is the upper limit of the transmission capacity of the section s;
Figure BDA0003852692730000096
a relaxation variable which is a lower limit of a transmission capacity of the section s; omega C Is a set of hydropower stations in the system; m is a group of C A hydraulic power station hydraulic energy abandoning penalty factor is given;
Figure BDA0003852692730000097
the power of the abandoned water of the hydropower station h in the time period t;
s3, establishing corresponding constraint conditions according to the objective function constructed in the step S2; the method specifically comprises the following steps:
the following equation is adopted as the system load balance constraint:
Figure BDA0003852692730000101
wherein N is the total number of the units in the system; p is i,t The output of the unit i in the time period t is obtained; NT is the total number of links in the system; t is j,t Planned power for tie j for time period t; d t System load for time period t;
the following equation is used as the system standby constraint:
Figure BDA0003852692730000102
Figure BDA0003852692730000103
in the formula of alpha i,t Is the start-stop state of the unit i in the time period t, and alpha i,t =1 indicates that the unit set i is in a start-up state in time t, α i,t =0 represents that the unit set i is in a shutdown state in the time period t;
Figure BDA0003852692730000104
the maximum output of the unit i in the time period t is obtained;
Figure BDA0003852692730000105
the minimum output of the unit i in the time period t is obtained;
Figure BDA0003852692730000106
the system positive spare capacity requirement for time period t;
Figure BDA0003852692730000107
the system negative spare capacity requirement for time period t;
the following formula is adopted as the output constraint of the unit:
Figure BDA0003852692730000108
in the formula P i,t The output of the unit i in the time period t is obtained;
the following formula is adopted as the unit climbing constraint:
Figure BDA0003852692730000109
Figure BDA00038526927300001010
in the formula,. DELTA.P i U The maximum upward slope climbing rate of the unit i is set; delta P i D The maximum downward climbing speed of the unit i;
the following formula is adopted as the minimum start-stop time constraint of the unit:
Figure BDA00038526927300001011
Figure BDA0003852692730000111
in the formula
Figure BDA0003852692730000112
For the time that the unit i has been continuously shut down during the time period t, an
Figure BDA0003852692730000113
Figure BDA0003852692730000114
Is the time that the unit i has been continuously started up in the time period t, an
Figure BDA0003852692730000115
T D Minimum continuous down time for the unit; t is U The minimum continuous starting time of the unit is set;
the following equation is used as the network security constraint:
Figure BDA0003852692730000116
Figure BDA0003852692730000117
in the formula P l min A lower limit for power flow transmission for line l; p is l max The upper limit of power flow transmission of the line l; g l-i Outputting a power transfer distribution factor for a generator of a line l by a node where a unit i is located; g l-j Outputting a power transfer distribution factor for the generator of the link line l by the node where the link line j is located; g l-k A generator output power transfer distribution factor for node k to line l; k is the total number of the network nodes; d k,t Load for node k at time period t;
Figure BDA0003852692730000118
is the lower limit of the power flow transmission of the section s;
Figure BDA0003852692730000119
is the upper limit of the power flow transmission of the section s; g s-i The generator output power of the section s is transferred to a distribution factor for the node where the unit i is located; g s-j The generator output power transfer distribution factor of the section s is calculated for the node where the tie line j is located; g s-k The distribution factor of the output power transfer of the generator is the node k to the section s;
the following formula is adopted as the discarded water electric quantity judgment constraint:
Figure BDA00038526927300001110
in the formula
Figure BDA00038526927300001111
The power of the hydropower station h is the power of the water abandoned in the time period t; a is h,t Is a variable from 0 to 1, and
Figure BDA00038526927300001112
when a h,t ＝1，
Figure BDA0003852692730000121
When a h,t ＝0；
Figure BDA0003852692730000122
The water discharge of the hydropower station h in the time period t; h h Is the water consumption rate of the hydropower station h; p is h,max The maximum adjustable output of the hydropower station h; p h,t The output of the hydropower station h in the time period t is obtained; m is a set maximum positive number and is used for carrying out linear transformation on the water abandonment electric quantity judgment constraint;
the following formula is adopted as the hydropower water level control constraint:
Figure BDA0003852692730000123
in the formula
Figure BDA0003852692730000124
For a hydropower station hThe dead water level of (c);
Figure BDA0003852692730000125
limiting the water level for the flood season of the hydropower station h; z h,0 The initial water level of the hydropower station h at the zero point of the next day; t is the current time period; p h,τ The output of the hydropower station h in the time period tau is obtained; h h Is the water consumption rate of the hydropower station h; s h The water surface area of the reservoir of the hydropower station h; f. of h,τ The amount of incoming water for the hydropower station h in the time period tau, and
Figure BDA0003852692730000126
I h,τ the natural inflow flow of the hydropower station h in the time period tau is defined as up (h), the upper hydropower station of the hydropower station h is defined as s (h), and the lag time of the upper hydropower station of the hydropower station h is defined as s (h);
the following formula is adopted as the constraint of the hydroelectric vibration region:
Figure BDA0003852692730000127
Figure BDA0003852692730000128
Figure BDA0003852692730000129
in the formula
Figure BDA00038526927300001210
The upper limit of the 1 st vibration area of the ith unit is set;
Figure BDA00038526927300001211
the lower limit of the 1 st vibration area of the ith unit; p is i,1,t The continuous variable represents the output of the operation interval when the time t of the ith unit is in the 1 st operable interval; delta i,l,t Is a variable of 0 to 1 and indicates whether the ith unit is in the t period
Figure BDA0003852692730000131
And if the ith unit is in the t period
Figure BDA0003852692730000132
Is between delta i,l,t =1, if the ith unit is not in t time interval
Figure BDA0003852692730000133
Is between delta i,l,t ＝0；
S4, solving the objective function constructed in the step S2 by adopting a variable scale optimization method under the constraint condition constructed in the step S3; the method specifically comprises the following steps:
A. determining a time interval step sequence, preferably 6h, 3h, 1.5h and 15min or 8h, 4h, 2h, 1h and 15min;
B. b, establishing a plurality of optimization model sets with the same optimization period and control requirements but different step sizes according to the time interval step length sequence selected in the step A; sequencing the optimization models in the sequence from large step to small step, and setting the serial number of the current optimization model to be solved as N = N Δ Corresponding time interval step is delta n (ii) a The optimization model comprises a cascade hydropower day-ahead spot shipment clearing objective function and a corresponding constraint condition;
C. for step size of Δ n The optimization model needs to describe the time interval coupling type constraint in the basic model again, N is more than or equal to 1 and less than or equal to N Δ ，N Δ To optimize the total number of problem sets, the larger n represents the larger step size of the problem, so Δ 1 Namely the original problem step length; re-describing the time-interval coupled constraints in the optimization model:
equivalent conversion is carried out on the unit climbing constraint by adopting the following formula:
Figure BDA0003852692730000134
Figure BDA0003852692730000135
in the formula P i,t The output of the unit i in the time period t is obtained; delta P i U The maximum upward slope climbing rate of the unit i is set; delta 1 Step length of the original optimization model;
Figure BDA0003852692730000136
the maximum climbing upper limit between adjacent time periods of the current optimization model is set; delta P i D The maximum downward climbing speed is set for the unit i;
Figure BDA0003852692730000137
setting the maximum landslide upper limit between adjacent time periods of the current optimization model; if it is
Figure BDA0003852692730000141
Or
Figure BDA0003852692730000142
If the current optimization model is larger than the rated capacity of the unit, determining that the current optimization model has no corresponding constraint;
the hydropower water level control constraint is equivalently converted by adopting the following formula:
Figure BDA0003852692730000143
in the formula
Figure BDA0003852692730000144
Is the dead water level of the hydropower station h;
Figure BDA0003852692730000145
limiting the water level for the flood season of the hydropower station h; z is a linear or branched member h,0 The initial water level of the hydropower station h at zero point of the next day; p is h,τ The output of the hydropower station h in the time interval tau is obtained; h h The water consumption rate of the hydropower station h;
Figure BDA0003852692730000146
the water discharge of the hydropower station h in the time period t; i is h,τ For hydropower stationsh natural incoming water flow at time period τ; s h The water surface area of the reservoir of the hydropower station h;
Figure BDA0003852692730000147
converting the water flow delay time of the upstream power station of the hydropower station h in the current optimization model, and if so
Figure BDA0003852692730000148
Determining that the influence of the ex-warehouse flow of the upstream hydropower station on the downstream in the current optimization model is immediate, and not considering the influence of time periods;
D. if N = N Δ Then solving the corresponding optimization model to obtain an initial solution; otherwise, setting the initial output process of each unit in the optimization model n as the optimization processing result of the last stage n +1,
Figure BDA0003852692730000149
and is provided with
Figure BDA00038526927300001410
Calculating to obtain the corresponding outlet flow and water level of the hydropower station;
after determining the initial solution of part of time period in the optimization model n, solving again;
E. updating the value of n to be n-1; if n is more than or equal to 1, skipping to the step C to continue to solve the optimization model of the next step; otherwise, ending the variable-scale optimization solving process to obtain a corresponding solving result;
and S5, according to the solving result of the step S4, finishing the daily spot shipment of the target step hydropower system.
FIG. 2 is a schematic diagram of functional modules of the system of the present invention: the invention discloses a system for realizing the stepped hydropower station daily spot shipment clearing method based on time interval decoupling, which comprises a data acquisition module, an objective function construction module, a constraint condition construction module, a model solving module and a data output module; the data acquisition module, the objective function construction module, the constraint condition construction module, the model solving module and the data output module are sequentially connected in series; the data acquisition module is used for acquiring the working data of the target step hydropower system and uploading the working data to the target function construction module; the target function construction module is used for constructing a step hydropower station spot shipment and clearing target function in the day ahead according to the acquired working data and uploading the target function to the constraint condition construction module; the constraint condition construction module is used for establishing corresponding constraint conditions according to the constructed objective function and uploading the constraint conditions to the model solving module; the model solving module is used for solving the objective function by adopting a variable-scale optimization method and uploading the result to the data output module; and the data output module is used for outputting the final daily spot shipment result of the target step hydropower system.

Claims (6)

1. A step hydropower station daily spot shipment clearing method based on time interval decoupling comprises the following steps:
s1, acquiring working data of a target step hydropower system;
s2, constructing a step hydropower station spot shipment clearing objective function in the day ahead according to the working data acquired in the step S1;
s3, establishing corresponding constraint conditions according to the objective function constructed in the step S2;
s4, solving the objective function constructed in the step S2 by adopting a variable scale optimization method under the constraint condition constructed in the step S3;
and S5, according to the solving result of the step S4, finishing the daily spot shipment of the target step hydropower system.
2. The method for removing the coupling of the time periods based on the step hydropower station daily spot shipment according to claim 1, wherein the step S2 is used for constructing a step hydropower station daily spot shipment objective function according to the work data acquired in the step S1, and comprises the following steps:
the following equation is used as the objective function:
Figure FDA0003852692720000011
n is the total number of the units in the system; t is the optimized total time period number; p i,t The output of the unit i in the time period t is obtained; c i,t (P i,t ) The operation cost of the unit i in the time period t is calculated;
Figure FDA0003852692720000012
starting cost of the unit i in a time period t; NL is the total number of lines in the system; m 1 A penalty value for network constraints;
Figure FDA0003852692720000013
a slack variable that is the upper limit of the transmission capacity of line l;
Figure FDA0003852692720000014
a slack variable which is the lower limit of the transmission capacity of the line l; NS is the total number of sections in the system;
Figure FDA0003852692720000015
a relaxation variable which is the upper limit of the transmission capacity of the section s;
Figure FDA0003852692720000016
a relaxation variable that is a lower limit of a transmission capacity of the section s; omega C Is a set of hydropower stations in the system; m is a group of C A hydraulic power station hydraulic energy abandoning penalty factor is given;
Figure FDA0003852692720000017
the power of the hydropower station h is abandoned in the time period t.
3. The method for releasing the time-interval coupling based daily spot shipment of the step hydropower station as claimed in claim 2, wherein the step S3 establishes corresponding constraint conditions according to the objective function established in the step S2, and specifically comprises the following steps:
the following equation is adopted as the system load balance constraint:
Figure FDA0003852692720000021
wherein N is the total number of the units in the system; p i,t The output of the unit i in the time period t is obtained; NT is the total number of links in the system; t is a unit of j,t Planned power for tie j for time period t; d t System load for time period t;
the following equation is used as the system standby constraint:
Figure FDA0003852692720000022
Figure FDA0003852692720000023
in the formula of alpha i,t Is the start-stop state of the unit i in the time period t, and alpha i,t =1 represents that the unit set i is in a starting state in a time period t, and alpha represents that the unit set i is in a starting state in a time period t i,t =0 represents that the unit set i is in a shutdown state in the time period t;
Figure FDA0003852692720000024
the maximum output of the unit i in the time period t is obtained;
Figure FDA0003852692720000025
the minimum output of the unit i in the time period t is obtained;
Figure FDA0003852692720000026
the system positive spare capacity requirement for time period t;
Figure FDA0003852692720000027
the system negative spare capacity requirement of the time interval t;
the following formula is adopted as the output constraint of the unit:
Figure FDA0003852692720000028
in the formula P i,t The output of the unit i in the time period t is obtained;
the following formula is adopted as the unit climbing constraint:
Figure FDA0003852692720000029
Figure FDA00038526927200000210
in the formula
Figure FDA00038526927200000211
The maximum upward climbing rate of the unit i; delta P i D The maximum downward climbing speed of the unit i is obtained;
the following formula is adopted as the minimum start-stop time constraint of the unit:
Figure FDA0003852692720000031
Figure FDA0003852692720000032
in the formula
Figure FDA0003852692720000033
For the time that the unit i has been continuously shut down during the time period t, an
Figure FDA0003852692720000034
Figure FDA0003852692720000035
Is the time that the unit i has been continuously started up in the time period t, an
Figure FDA0003852692720000036
T D Minimum continuous down time for the unit; t is U The minimum continuous starting time of the unit is set;
the following formula is adopted as the network security constraint:
Figure FDA0003852692720000037
Figure FDA0003852692720000038
in the formula P l min A lower limit for power flow transmission of line l; p is l max An upper limit for the power flow transmission of line l; g l-i Outputting a power transfer distribution factor for the generator of the line l by the node where the unit i is located; g l-j Outputting a power transfer distribution factor for the generator of the link line l by the node where the link line j is located; g l-k A generator output power transfer distribution factor for node k to line l; k is the total number of the network nodes; d k,t The load for node k at time period t; p is s min The lower limit of power flow transmission of the section s; p s max Is the upper limit of the power flow transmission of the section s; g s-i The generator output power of the section s is transferred to a distribution factor for the node where the unit i is located; g s-j The generator output power transfer distribution factor of the section s is calculated for the node where the tie line j is located; g s-k The distribution factor of the output power transfer of the generator is the node k to the section s;
the following formula is adopted as the discarded water electric quantity judgment constraint:
Figure FDA0003852692720000039
Figure FDA00038526927200000310
in the formula
Figure FDA0003852692720000041
The power of the hydropower station h is the power of the water abandoned in the time period t; a is a h,t Is a variable from 0 to 1, and
Figure FDA0003852692720000042
when a h,t ＝1，
Figure FDA0003852692720000043
When a h,t ＝0；
Figure FDA0003852692720000044
The water discharge of the hydropower station h in the time period t is determined; h h The water consumption rate of the hydropower station h; p is h,max The maximum adjustable output of the hydropower station h; p is h,t For the hydropower station h during the time period t; m is a set maximum positive number and is used for carrying out linear transformation on the water abandonment electric quantity judgment constraint;
the following formula is adopted as the hydropower water level control constraint:
Figure FDA0003852692720000045
in the formula
Figure FDA0003852692720000046
Is the dead water level of the hydropower station h;
Figure FDA0003852692720000047
limiting the water level for the flood season of the hydropower station h; z h,0 The initial water level of the hydropower station h at zero point of the next day; t is the current time period; p is h,τ The output of the hydropower station h in the time period tau is obtained; h h The water consumption rate of the hydropower station h; s h The water surface area of the reservoir of the hydropower station h; f. of h,τ The amount of incoming water for hydropower station h in time period τ, and
Figure FDA0003852692720000048
I h,τ the natural inflow water flow of the hydropower station h in a time interval tau is shown as up (h), the upper hydropower station of the hydropower station h is shown as s (h), and the lag time of the upper hydropower station of the hydropower station h is shown as s (h);
the following formula is adopted as the constraint of the hydroelectric vibration region:
Figure FDA0003852692720000049
Figure FDA00038526927200000410
Figure FDA00038526927200000411
in the formula
Figure FDA00038526927200000412
The upper limit of the 1 st vibration area of the ith unit is set;
Figure FDA00038526927200000413
the lower limit of the 1 st vibration area of the ith unit; p i,1,t The output value of the operation interval is a continuous variable and represents the output value of the operation interval when the ith unit is in the 1 st operable interval in time period t; delta i,l,t Is a variable of 0 to 1 and indicates whether the ith unit is in the t period
Figure FDA0003852692720000051
And if the ith unit is in the t period
Figure FDA0003852692720000052
Is between delta i,l,t =1, if the ith unit is not in t time interval
Figure FDA0003852692720000053
Delta between i,l,t ＝0。
4. The off-site daily delivery method of step hydropower based on decoupling of time period coupling according to claim 3, characterized in that the step S4 is implemented by adopting a variable scale optimization method, and solving the objective function constructed in the step S2 under the constraint condition constructed in the step S3, and specifically comprises the following steps:
A. determining a sequence of interval step sizes;
B. b, establishing a plurality of optimization model sets with the same optimization period and control requirements but different step sizes according to the time interval step length sequence selected in the step A; sequencing the optimization models in the sequence from large step to small step, and setting the serial number of the current optimization model to be solved as N = N Δ Corresponding time interval step is delta n (ii) a The optimization model comprises a cascade hydropower day-ahead spot shipment clearing objective function and corresponding constraint conditions;
C. for step size of Δ n The time-interval coupling type constraint in the optimization model is described again:
equivalent conversion is carried out on the unit climbing constraint by adopting the following formula:
Figure FDA0003852692720000054
Figure FDA0003852692720000055
in the formula P i,t The output of the unit i in the time period t is obtained; delta P i U The maximum upward climbing rate of the unit i; delta 1 Step length of the original optimization model;
Figure FDA0003852692720000056
the maximum climbing upper limit between adjacent time periods of the current optimization model is set; delta P i D The maximum downward climbing speed of the unit i is obtained;
Figure FDA0003852692720000057
the maximum landslide upper limit between adjacent time periods of the current optimization model is set; if it is
Figure FDA0003852692720000061
Or
Figure FDA0003852692720000062
If the current optimization model is larger than the rated capacity of the unit, determining that the current optimization model has no corresponding constraint;
the hydropower water level control constraint is equivalently converted by adopting the following formula:
Figure FDA0003852692720000063
in the formula
Figure FDA0003852692720000064
Is the dead water level of the hydropower station h;
Figure FDA0003852692720000065
limiting the water level for the flood season of the hydropower station h; z is a linear or branched member h,0 The initial water level of the hydropower station h at zero point of the next day; p is h,τ The output of the hydropower station h in the time interval tau is obtained; h h Is the water consumption rate of the hydropower station h;
Figure FDA0003852692720000066
the water discharge of the hydropower station h in the time period t is determined; i is h,τ The natural incoming water flow of the hydropower station h in the time period tau is obtained; s h The water surface area of the reservoir of the hydropower station h;
Figure FDA0003852692720000067
for the water flow delay time conversion of the upstream power station of the hydropower station h in the current optimization model, and if so
Figure FDA0003852692720000068
Determining that the influence of the ex-warehouse flow of the upstream hydropower station on the downstream in the current optimization model is immediate, and not considering the influence of time periods;
D. if N = N Δ Then solving the corresponding optimization model to obtain an initial solution; otherwise, setting the initial output process of each unit in the optimization model n as the optimization processing result of the last stage n +1,
Figure FDA0003852692720000069
and is provided with
Figure FDA00038526927200000610
Calculating to obtain the corresponding outlet flow and water level of the hydropower station;
when the initial solution of part of time intervals in the optimization model n is determined, solving again;
E. updating the value of n to be n-1; if n is larger than or equal to 1, skipping to the step C to continuously solve the next long optimization model; otherwise, ending the variable-scale optimization solving process to obtain a corresponding solving result.
5. The method for off-the-shelf delivery of step hydropower plants on day basis based on decoupling of time periods according to claim 4, wherein the time period step sequence of step A is 6h, 3h, 1.5h and 15min or 8h, 4h, 2h, 1h and 15min.
6. A system for realizing the stepped hydropower day-ahead spot shipment method based on decoupling of the time period coupling in claims 1-5 is characterized by comprising a data acquisition module, an objective function construction module, a constraint condition construction module, a model solving module and a data output module; the data acquisition module, the target function construction module, the constraint condition construction module, the model solving module and the data output module are sequentially connected in series; the data acquisition module is used for acquiring the working data of the target step hydropower system and uploading the working data to the target function construction module; the target function construction module is used for constructing a step hydropower station spot shipment and clearing target function in the day ahead according to the acquired working data and uploading the target function to the constraint condition construction module; the constraint condition construction module is used for establishing corresponding constraint conditions according to the constructed objective function and uploading the constraint conditions to the model solving module; the model solving module is used for solving the objective function by adopting a variable-scale optimization method and uploading the result to the data output module; and the data output module is used for outputting the final daily spot shipment result of the target step hydropower system.
CN202211139082.9A 2022-09-19 2022-09-19 Method and system for clearing cascade hydropower day-ahead spot goods based on time interval decoupling Pending CN115423202A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211139082.9A CN115423202A (en) 2022-09-19 2022-09-19 Method and system for clearing cascade hydropower day-ahead spot goods based on time interval decoupling

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211139082.9A CN115423202A (en) 2022-09-19 2022-09-19 Method and system for clearing cascade hydropower day-ahead spot goods based on time interval decoupling

Publications (1)

Publication Number Publication Date
CN115423202A true CN115423202A (en) 2022-12-02

Family

ID=84205197

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211139082.9A Pending CN115423202A (en) 2022-09-19 2022-09-19 Method and system for clearing cascade hydropower day-ahead spot goods based on time interval decoupling

Country Status (1)

Country Link
CN (1) CN115423202A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115983733A (en) * 2023-01-03 2023-04-18 中国南方电网有限责任公司 Power market clearing data processing method and device based on water level control

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115983733A (en) * 2023-01-03 2023-04-18 中国南方电网有限责任公司 Power market clearing data processing method and device based on water level control
CN115983733B (en) * 2023-01-03 2023-10-20 中国南方电网有限责任公司 Electric power market clear data processing method and device based on water level control

Similar Documents

Publication Publication Date Title
CN107276122B (en) Peak-shaving resource calling decision method suitable for large-scale renewable energy grid connection
CN112701687B (en) Robust optimization operation method of gas-electricity distribution network system considering price type combined demand response
CN110400056B (en) Method and device for cascade hydropower day-ahead optimized scheduling based on wind-light space-time characteristics
CN110826773B (en) Thermal power generating unit monthly power generation plan optimization method considering new energy access
CN108322139A (en) It is a kind of meter and overall cost wind-powered electricity generation-photovoltaic-photo-thermal power generation combined scheduling method
CN109687534B (en) Power system generator set active power control method based on step water quantity matching
CN114400698B (en) Optimal carbon reduction operation method for high-proportion clean energy power grid power supply
CN109149656B (en) Gas-electricity interconnection comprehensive energy system unit combination method
CN115423202A (en) Method and system for clearing cascade hydropower day-ahead spot goods based on time interval decoupling
CN112072710A (en) Source network load integrated economic dispatching method and system considering demand response
CN114243691A (en) Low-carbon optimized scheduling method for power system
CN116979578A (en) Electric and thermal triple generation optimal scheduling method and system for wind, light, water and fire storage
CN115455709A (en) Low-carbon comprehensive energy system simulation and configuration method considering carbon capture equipment
CN114757509A (en) Day-ahead optimized scheduling method for hydropower enrichment power grid in spot market environment
CN110867907A (en) Power system scheduling method based on multi-type power generation resource homogenization
CN117013522A (en) Comprehensive energy system scheduling optimization method considering distributed power supply and gas-electricity cooperation
CN116540545A (en) Photovoltaic power generation hydrogen production cluster random optimization scheduling method based on ember process
CN114884101A (en) Pumped storage dispatching method based on self-adaptive model control prediction
CN116231767A (en) Multi-energy complementary scheduling method and system for cascade hydropower station
CN111799842B (en) Multi-stage power transmission network planning method and system considering flexibility of thermal power generating unit
CN115936336A (en) Virtual power plant capacity configuration and regulation operation optimization method
CN110854929B (en) Day-ahead scheduling method considering uncertainty in time period
CN114465226A (en) Method for establishing multi-level standby acquisition joint optimization model of power system
CN114362255A (en) Multi-target day-ahead scheduling optimization method and system for source-network charge storage power system
CN113240546A (en) Monthly scheduling method for units in dense hydropower region

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