CN111428922B - Electric power market clearing method and system for rapid posterior condition section - Google Patents

Electric power market clearing method and system for rapid posterior condition section Download PDF

Info

Publication number
CN111428922B
CN111428922B CN202010196782.6A CN202010196782A CN111428922B CN 111428922 B CN111428922 B CN 111428922B CN 202010196782 A CN202010196782 A CN 202010196782A CN 111428922 B CN111428922 B CN 111428922B
Authority
CN
China
Prior art keywords
section
constraint
unit
scuc
model
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
CN202010196782.6A
Other languages
Chinese (zh)
Other versions
CN111428922A (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.)
Yunnan Power Grid Co Ltd
Original Assignee
Beijing Tsintergy Technology Co ltd
Yunnan Power Grid 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 Beijing Tsintergy Technology Co ltd, Yunnan Power Grid Co Ltd filed Critical Beijing Tsintergy Technology Co ltd
Priority to CN202010196782.6A priority Critical patent/CN111428922B/en
Publication of CN111428922A publication Critical patent/CN111428922A/en
Application granted granted Critical
Publication of CN111428922B publication Critical patent/CN111428922B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • 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
    • G06Q30/00—Commerce
    • G06Q30/02—Marketing; Price estimation or determination; Fundraising
    • G06Q30/0201—Market modelling; Market analysis; Collecting market data
    • G06Q30/0206—Price or cost determination based on market factors
    • 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas or water supply

Abstract

The invention discloses a quick electric power market clearing method for posterior condition sections, which comprises the following steps: establishing an SCUC model according to preset constraint conditions and an SCUC optimization target, and solving the SCUC model to obtain the output and start-stop states of each unit at different time periods; calculating the section flow according to the output and start-stop states of each unit at different time intervals to obtain branch flow and section flow results; judging whether the conditional section constraints or the control section constraints falling into the conditional sections are respectively locked into the effective conditional sections according to the branch flow and section flow results, and respectively adding the effective conditional sections into the SCUC model for iterative solution; when any one of the new SCUC models is successfully solved, the other SCUC model is subjected to termination operation, and the obtained result is used as an optimal mode; and determining constraint conditions and an optimization target according to the optimal mode, establishing an SCED model, solving to obtain Lagrange multiplier values corresponding to constraints of each unit at different time intervals, and calculating to obtain the electricity prices of each node at different time intervals.

Description

Electric power market clearing method and system for rapid posterior condition section
Technical Field
The invention relates to the technical field of power markets, in particular to a method and a system for clearing a power market of a section under a rapid posterior condition.
Background
With the gradual advance of the construction of the electric power spot market in China, the mode of the traditional electric power dispatching and power generation planning is changed, the constraint conditions of electric power supply and demand balance, unit physical characteristics, power grid safety and the like need to be comprehensively considered, and the day-ahead and real-time power generation planning of the marketized unit is completed with the aim of minimum cost on the power generation side or maximum social benefits. Modeling and solving of a Unit combination (SCUC) and an economic Dispatch (SCED) of safety constraint are core problems of spot market clearing, how to process the safety constraint of a power grid in the clearing process is directly related to the coordination level of system economy and reliability, and direct influence is generated on a market clearing result.
For a long time, domestic and foreign research scholars have developed a great deal of research on power generation scheduling plan optimization under market conditions, and provide reliable theoretical basis and engineering practice method for power generation scheduling plan optimization under market conditions, but power grid safety constraint is mainly expressed by normalized branch or section tidal current limit values, and the problem of condition control section is not fully considered. Based on the existing research results, the SCUC solving process in the existing market clearing calculation framework in the day before can only control section constraints without considering complex conditions, and the starting and effective conditions of the conditional section are judged according to the power generation plan obtained by solving; and then, adjusting the power generation plan through a closed-loop correction optimization decision to ensure that the condition control section constraint is met.
However, in the closed-loop correction optimization decision mode according to the latest out-of-limit result, although the calculation scale of each iteration is reduced, endless iteration times are easy to generate, and the overall efficiency of the solution is greatly reduced. Therefore, a scientific, effective and rapid power market clearing method and system considering the condition section of the power system are urgently needed in the actual industry.
Disclosure of Invention
The invention provides a method and a system for clearing an electric power market of a rapid posterior condition section, which can reduce the loss of an optimization effect, reduce the working difficulty of decision correction and improve the efficiency of iterative computation.
In order to solve the above technical problem, an embodiment of the present invention provides a method for clearing a power market of a section under a rapid posterior condition, including:
establishing an SCUC model according to preset constraint conditions and an SCUC optimization target, and solving the SCUC model by a branch-and-bound method to obtain the output and start-stop states of each unit at different time intervals;
calculating the section tidal current by a direct current tidal current method according to the output and start-stop states of the units at different time intervals to obtain branch tidal current and section tidal current results;
judging the conditional section constraint or the control section constraint falling into the conditional section according to the branch flow and section flow results, respectively locking the effective conditional sections, respectively adding the effective conditional sections into the SCUC model to obtain two new SCUC models, and respectively carrying out iterative solution on the two new SCUC models;
when any one new SCUC model is successfully solved, carrying out termination operation on the other new SCUC model, and taking a result obtained after the solution is successfully solved as an optimal mode;
determining SCED constraint conditions and an SCED optimization target according to the optimal mode, establishing an SCED model, solving the SCED model through a dual simplex method to obtain Lagrangian multiplier values corresponding to constraints of each unit in different time periods, and calculating the electrovalence of each node in different time periods according to the Lagrangian multiplier values corresponding to the constraints in different time periods.
As a preferred scheme, the SCUC optimization target is:
Figure BDA0002417919190000021
wherein: n represents the total number of the units; t represents the total number of time segments considered;
Pi,trepresenting the output of the unit i in the time period t;
Figure BDA0002417919190000022
the running cost, the starting cost and the shutdown cost of the unit i in the time period t are respectively, wherein the unit running cost Ci,t(Pi,t) Is a multi-segment linear function of the unit output.
Preferably, the preset constraint condition includes: the method comprises the following steps of load balance constraint, unit output upper and lower limit constraint, unit climbing constraint, unit minimum continuous start-stop time constraint and unit maximum start-stop times constraint.
Preferably, the formula for calculating the cross-sectional power flow by the direct current power flow method is as follows:
Figure BDA0002417919190000031
wherein G isl-iOutputting a power transfer distribution factor for a generator of a line or a section l by a node where a unit i is located; k is the number of nodes of the system; gl-kThe distribution factor of the output power transfer of the generator of the node k to the line or the section l; dk,tIs the bus load value of node k in the period t.
Preferably, the SCED constraint includes: load balance constraint, unit output upper and lower limit constraint, unit climbing constraint and conditional section constraint.
As a preferred scheme, the SCED optimization target is:
Figure BDA0002417919190000032
wherein: n represents the total number of the units; t represents the total number of time segments considered;
Pi,trepresenting the output of the unit i in the time period t; ci,t(Pi,t) The operating cost of the unit i in the time period t is a multi-segment linear function of the unit output.
As a preferred scheme, the calculation formula for calculating the electricity prices of the nodes in different time periods according to the lagrangian multiplier values corresponding to the constraints in different time periods is as follows:
Figure BDA0002417919190000033
wherein λ ist: time period t system loadA balanced constrained lagrange multiplier;
Figure BDA0002417919190000034
a lagrange multiplier for the maximum forward power flow constraint of the line l;
Figure BDA0002417919190000035
a lagrangian multiplier constrained by the maximum reverse power flow of the line l;
Figure BDA0002417919190000036
a lagrange multiplier with section s maximum forward power flow constraint;
Figure BDA0002417919190000037
a lagrange multiplier with section s maximum reverse power flow constraint; gl-k: the node k transfers the distribution factor to the generator output power of the line l; gs-k: the node k transfers the distribution factor to the generator output power of the section s.
The embodiment of the invention also provides a rapid posterior condition section electric power market clearing system, which comprises:
the model establishing module is used for establishing an SCUC model according to preset constraint conditions and an SCUC optimization target, and solving the SCUC model through a branch-and-bound method to obtain the output and start-stop states of each unit at different time intervals;
the section power flow module is used for calculating the section power flow through a direct current power flow method according to the output and start-stop states of the units at different time intervals to obtain branch power flow and section power flow results;
the iteration solving module is used for judging the conditional section constraint or the control section constraint falling into the conditional section according to the branch flow and section flow results, respectively locking the effective conditional section, respectively adding the effective conditional section into the SCUC model to obtain two new SCUC models, and respectively carrying out iteration solving on the two new SCUC models;
the termination operation module is used for carrying out termination operation on another new SCUC model after any new SCUC model is successfully solved, and taking a result obtained after the solution is successfully solved as an optimal mode;
and the electricity price calculation module is used for determining the SCED constraint condition and the SCED optimization target according to the optimal mode, establishing an SCED model, solving the SCED model through a dual simplex method to obtain Lagrange multiplier values corresponding to the constraints of each unit in different time periods, and calculating the electricity price of each node in different time periods according to the Lagrange multiplier values corresponding to the constraints in different time periods.
As a preferred scheme, the SCUC optimization target is:
Figure BDA0002417919190000041
wherein: n represents the total number of the units; t represents the total number of time segments considered;
Pi,trepresenting the output of the unit i in the time period t; ci,t(Pi,t)、
Figure BDA0002417919190000042
The running cost, the starting cost and the shutdown cost of the unit i in the time period t are respectively, wherein the unit running cost Ci,t(Pi,t) Is a multi-segment linear function of the unit output.
Preferably, the preset constraint condition includes: the method comprises the following steps of load balance constraint, unit output upper and lower limit constraint, unit climbing constraint, unit minimum continuous start-stop time constraint and unit maximum start-stop times constraint.
Preferably, the formula for calculating the cross-sectional power flow by the direct current power flow method is as follows:
Figure BDA0002417919190000043
wherein G isl-iOutputting a power transfer distribution factor for a generator of a line or a section l by a node where a unit i is located; k is the number of nodes of the system; gl-kFor node k to line or sectionl generator output power transfer distribution factor; dk,tIs the bus load value of node k in the period t.
Preferably, the SCED constraint includes: load balance constraint, unit output upper and lower limit constraint, unit climbing constraint and conditional section constraint.
As a preferred scheme, the SCED optimization target is:
Figure BDA0002417919190000051
wherein: n represents the total number of the units; t represents the total number of time segments considered;
Pi,trepresenting the output of the unit i in the time period t; ci,t(Pi,t) The operating cost of the unit i in the time period t is a multi-segment linear function of the unit output.
As a preferred scheme, the calculation formula for calculating the electricity prices of the nodes in different time periods according to the lagrangian multiplier values corresponding to the constraints in different time periods is as follows:
Figure BDA0002417919190000052
wherein λ ist: lagrangian multipliers of system load balance constraint in the period t;
Figure BDA0002417919190000053
a lagrange multiplier for the maximum forward power flow constraint of the line l;
Figure BDA0002417919190000054
a lagrangian multiplier constrained by the maximum reverse power flow of the line l;
Figure BDA0002417919190000055
a lagrange multiplier with section s maximum forward power flow constraint;
Figure BDA0002417919190000056
a lagrange multiplier with section s maximum reverse power flow constraint; gl-k: the node k transfers the distribution factor to the generator output power of the line l; gs-k: the node k transfers the distribution factor to the generator output power of the section s.
An embodiment of the present invention further provides a computer-readable storage medium, where the computer-readable storage medium includes a stored computer program; wherein the computer program, when executed, controls an apparatus in which the computer-readable storage medium is located to perform the method of rapid posterior conditional profiling for electricity market clearing as defined in any of the above.
An embodiment of the present invention further provides a terminal device, which includes a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, where the processor, when executing the computer program, implements the power market clearing method for a rapid posterior condition section as described in any one of the above.
Compared with the prior art, the embodiment of the invention has the following beneficial effects:
according to the technical scheme, the SCUC model is subjected to iterative solution to obtain an optimal mode, the SCED model is established for solution, the electricity price of each node in different time periods is calculated, the loss of the optimization effect is reduced, the working difficulty of decision correction is reduced, and the efficiency of iterative calculation is improved.
Drawings
FIG. 1: the invention provides a flow chart of a method for clearing a power market of a rapid posterior condition section in an embodiment.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Referring to fig. 1, a preferred embodiment of the present invention provides a method for clearing a power market with a rapid posterior conditional section, including:
and S1, establishing an SCUC model according to preset constraint conditions and an SCUC optimization target, and solving the SCUC model by a branch-and-bound method to obtain the output and start-stop states of each unit at different time intervals.
Specifically, the objective function of the future market clearing SCUC is to minimize the electricity purchase cost:
Figure BDA0002417919190000061
wherein:
n represents the total number of the units; t represents the total number of considered time periods, and assuming 24 time periods are considered in a day, T is 24; pi,tRepresenting the output of the unit i in the time period t; ci,t(Pi,t)、
Figure BDA0002417919190000062
The operating costs, the start-up costs and the shutdown costs (optionally considered) for the unit i in the time period t, respectively, wherein the unit operating costs Ci,t(Pi,t) Is a multi-segment linear function of the unit output.
Specifically, the constraints include:
1. a load balancing constraint, which for each time period t may be described as:
Figure BDA0002417919190000071
wherein, Pi,tRepresenting the output of the unit i in the time period t, DtFor the time t period, the system load has deducted the tie-line net injected power.
2. The upper and lower limits of the unit output force are restricted, the output force of the unit should be within the maximum/minimum technical output range, and the restriction conditions can be described as follows:
Figure BDA0002417919190000072
if the unit is shut down, αi,tIf the output power of the unit is 0, the output power of the unit can be limited to 0 by the constraint condition; when the unit is started, alphai,tThe constraint is a conventional upper and lower force limit constraint.
3. The unit climbing restriction is that the unit climbing speed requirement should be met when climbing up or down. The hill climbing constraint can be described as:
Pi,t-Pi,t-1≤ΔPi Uαi,t-1
Pi,t-1-Pi,t≤ΔPi Dαi,t
wherein, Δ Pi UThe maximum upward slope climbing rate of the unit i,
Figure BDA0002417919190000073
the maximum downward climbing rate of the unit i.
4. The minimum continuous start-stop time of the thermal power generating unit is constrained, and the thermal power generating unit is required to meet the minimum continuous start-stop time due to the physical properties and actual operation requirements of the thermal power generating unit. The minimum continuous on-off time constraint can be described as:
Figure BDA0002417919190000074
Figure BDA0002417919190000075
wherein alpha isi,tStarting and stopping a unit i at a time t; t isU、TDThe minimum continuous start-up time and the minimum continuous stop time of the unit.
5. The maximum start-stop times of the unit are constrained, firstly, the switching variable of start and stop is defined, and eta is definedi,tWhether the unit i is switched to a starting state in a time period t or not is judged; definition of gammai,tAnd indicating whether the unit i is switched to a shutdown state in the period t. Starting and stopping of corresponding unit iThe frequency limit can be expressed as follows:
Figure BDA0002417919190000076
Figure BDA0002417919190000077
wherein the content of the first and second substances,
Figure BDA0002417919190000081
the maximum starting and stopping times of the unit i are respectively.
And then, solving the models by using a branch-and-bound method to obtain the output and the start-stop state of each unit in different time periods.
And S2, calculating the section flow through a direct current flow method according to the output and start-stop states of the units at different time intervals to obtain branch flow and section flow results.
Specifically, the output and the start-stop state of each unit obtained in the step 1 at different time periods are input. Calculating branch flow and section flow by using the output power transfer distribution factor of the generator:
Figure BDA0002417919190000082
wherein G isl-iOutputting a power transfer distribution factor for a generator of a line or a section l by a node where a unit i is located; k is the number of nodes of the system; gl-kThe distribution factor of the output power transfer of the generator of the node k to the line or the section l; dk,tIs the bus load value of node k in the period t.
And S3, judging the conditional section constraint or the control section constraint falling in the conditional section according to the branch flow and section flow results, respectively locking the effective conditional section, respectively adding the effective conditional section into the SCUC model to obtain two new SCUC models, and respectively carrying out iterative solution on the two new SCUC models.
Specifically, the form is a general expression form of the conditional section as in formula (1). When the "conditional part" section ST is in the range [ STj,min,STj,max]The limit of the "control section" section S is [ S ]j,min,Sj,max]：
Figure BDA0002417919190000083
The following two modes are simultaneously calculated in parallel:
mode (1): inputting the branch flow and section flow results obtained in the step 2, locking the section with the effective condition according to the section condition of the 'condition part' falling in the condition section, adding the section with the effective condition to the model of 1.1, and re-iterating and solving the SCUC;
mode (2): and (3) inputting the branch flow and section flow results obtained in the step (2), locking the section with the effective condition according to the section condition of the 'control part' falling into the condition section, adding the section with the effective condition into the model of 1.1, and re-iterating and solving the SCUC.
And S4, when any one new SCUC model is successfully solved, carrying out termination operation on the other new SCUC model, and taking the result obtained after the solution is successfully solved as an optimal mode.
Specifically, in step 3, once the solution of any mode is successful, the calculation of the other mode is terminated, and the mode with the successful solution is selected as the optimal mode.
S5, determining SCED constraint conditions and SCED optimization targets according to the optimal mode, establishing an SCED model, solving the SCED model through a dual simplex method to obtain Lagrangian multiplier values corresponding to constraints of each unit in different time periods, and calculating the electrovalence of each node in different time periods according to the Lagrangian multiplier values corresponding to the constraints in different time periods.
Specifically, the SCED model (containing optimization objectives and constraints) is constructed.
Wherein, the objective function with the optimization target of coming out of SCED in the market at the day before is the minimization of the electricity purchasing cost:
Figure BDA0002417919190000091
in the above formula: n represents the total number of the units; t represents the total number of considered time periods, and assuming 24 time periods are considered in a day, T is 24; pi,tRepresenting the output of the unit i in the time period t; ci,t(Pi,t) The operating cost of the unit i in the time period t is a multi-segment linear function of the unit output.
Wherein the constraint condition comprises:
1. a load balancing constraint, which for each time period t may be described as:
Figure BDA0002417919190000092
wherein, Pi,tRepresenting the output of the unit i in the time period t, DtFor the time t period, the system load has deducted the tie-line net injected power.
2. The upper and lower limits of the unit output force are restricted, the output force of the unit should be within the maximum/minimum technical output range, and the restriction conditions can be described as follows:
Figure BDA0002417919190000093
if the unit is shut down, αi,tIf the output power of the unit is 0, the output power of the unit can be limited to 0 by the constraint condition; when the unit is started, alphai,tThe constraint is a conventional upper and lower force limit constraint.
3. The unit climbing restriction is that the unit climbing speed requirement should be met when climbing up or down. The hill climbing constraint can be described as:
Pi,t-Pi,t-1≤ΔPi Uαi,t-1
Pi,t-1-Pi,t≤ΔPi Dαi,t
wherein, Δ Pi UFor the unit i maximum climbing rate, Δ Pi DThe maximum downward climbing rate of the unit i.
4. And (4) adding corresponding conditional section constraints according to the optimal mode selected in the step 4.
Pl,min≤Pl≤Pl,max
Wherein, PlIs the tide of a line or a section l; pl,minThe lower current limit of the line or the section l; pl,maxIs the upper current limit of the line or section l.
And then, solving the model by adopting a dual simplex method to obtain the output of each unit in different periods and Lagrangian multiplier values corresponding to each constraint.
And calculating the electricity price of each node in different periods based on the following formula by combining the Lagrange multiplier value corresponding to each constraint in each period obtained in the steps. The node electricity price of the node i in the t period is as follows:
Figure BDA0002417919190000101
in the formula, λt: lagrangian multipliers of system load balance constraint in the period t;
Figure BDA0002417919190000102
a lagrange multiplier for the maximum forward power flow constraint of the line l;
Figure BDA0002417919190000103
a lagrangian multiplier constrained by the maximum reverse power flow of the line l;
Figure BDA0002417919190000104
a lagrange multiplier with section s maximum forward power flow constraint;
Figure BDA0002417919190000105
lagrange multiplication of maximum reverse power flow constraint of section sA seed; gl-k: the node k transfers the distribution factor to the generator output power of the line l; gs-k: the node k transfers the distribution factor to the generator output power of the section s.
It should be noted that all lagrange multipliers are equal to or greater than 0.
The establishment of the alternating current-direct current interconnected power system enables the optimal configuration of power resources to be realized in a large range, and simultaneously, the requirement on the operation control of a complex power grid is also improved. For example, a strong coupling relationship is formed between a part of ultrahigh-voltage direct current lines and a matching unit due to a short electrical distance, and when the limit of a control section related to the tie line is set, the operation condition (the number of power plant starting units, the output power of the power plant unit, and the like) of the corresponding matching unit needs to be considered. Therefore, whether the section constraint is started or not and the effective limit range after the section constraint is started need to be judged according to certain preconditions; whether the precondition is satisfied or not is usually related to decision variables of the SCUC model to be solved, such as unit start and stop, unit output and the like, and a constraint relation exists between the solution of the SCUC model and the setting of the condition control section constraint. The invention provides a method and a system for clearing an electric power market of a rapid posterior conditional section, which can lock the conditional part constraint and the control part constraint of the conditional section through the first SCUC calculation and carry out iterative calculation again in order to reduce the loss of an optimization effect and reduce the working difficulty of decision correction, and can complete the link of clearing calculation with the longest time only by twice SCUC.
Correspondingly, the embodiment of the invention also provides a rapid posterior conditional section electric power market clearing system, which comprises:
the model establishing module is used for establishing an SCUC model according to preset constraint conditions and an SCUC optimization target, and solving the SCUC model through a branch-and-bound method to obtain the output and start-stop states of each unit at different time intervals;
the section power flow module is used for calculating the section power flow through a direct current power flow method according to the output and start-stop states of the units at different time intervals to obtain branch power flow and section power flow results;
the iteration solving module is used for judging the conditional section constraint or the control section constraint falling into the conditional section according to the branch flow and section flow results, respectively locking the effective conditional section, respectively adding the effective conditional section into the SCUC model to obtain two new SCUC models, and respectively carrying out iteration solving on the two new SCUC models;
the termination operation module is used for carrying out termination operation on another new SCUC model after any new SCUC model is successfully solved, and taking a result obtained after the solution is successfully solved as an optimal mode;
and the electricity price calculation module is used for determining the SCED constraint condition and the SCED optimization target according to the optimal mode, establishing an SCED model, solving the SCED model through a dual simplex method to obtain Lagrange multiplier values corresponding to the constraints of each unit in different time periods, and calculating the electricity price of each node in different time periods according to the Lagrange multiplier values corresponding to the constraints in different time periods.
An embodiment of the present invention further provides a computer-readable storage medium, where the computer-readable storage medium includes a stored computer program; wherein the computer program, when running, controls the device on which the computer readable storage medium is located to execute the method for clearing a power market of a rapid posterior conditional section according to any of the above embodiments.
The embodiment of the present invention further provides a terminal device, where the terminal device includes a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, and the processor, when executing the computer program, implements the method for clearing an electric power market of a rapid posterior condition section according to any of the above embodiments.
Preferably, the computer program may be divided into one or more modules/units (e.g., computer program) that are stored in the memory and executed by the processor to implement the invention. The one or more modules/units may be a series of computer program instruction segments capable of performing specific functions, which are used for describing the execution process of the computer program in the terminal device.
The Processor may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), an off-the-shelf Programmable Gate Array (FPGA) or other Programmable logic device, a discrete Gate or transistor logic device, a discrete hardware component, etc., the general purpose Processor may be a microprocessor, or the Processor may be any conventional Processor, the Processor is a control center of the terminal device, and various interfaces and lines are used to connect various parts of the terminal device.
The memory mainly includes a program storage area and a data storage area, wherein the program storage area may store an operating system, an application program required for at least one function, and the like, and the data storage area may store related data and the like. In addition, the memory may be a high speed random access memory, may also be a non-volatile memory, such as a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash Card (Flash Card), and the like, or may also be other volatile solid state memory devices.
It should be noted that the terminal device may include, but is not limited to, a processor and a memory, and those skilled in the art will understand that the terminal device is only an example and does not constitute a limitation of the terminal device, and may include more or less components, or combine some components, or different components.
The above-mentioned embodiments are provided to further explain the objects, technical solutions and advantages of the present invention in detail, and it should be understood that the above-mentioned embodiments are only examples of the present invention and are not intended to limit the scope of the present invention. It should be understood that any modifications, equivalents, improvements and the like, which come within the spirit and principle of the invention, may occur to those skilled in the art and are intended to be included within the scope of the invention.

Claims (10)

1. A method for clearing a power market of a rapid posterior condition section is characterized by comprising the following steps:
establishing an SCUC model according to preset constraint conditions and an SCUC optimization target, and solving the SCUC model by a branch-and-bound method to obtain the output and start-stop states of each unit at different time intervals;
calculating the section tidal current by a direct current tidal current method according to the output and start-stop states of the units at different time intervals to obtain branch tidal current and section tidal current results;
judging the conditional section constraint or the control section constraint falling into the conditional section according to the branch flow and section flow results, respectively locking the effective conditional sections, respectively adding the effective conditional sections into the SCUC model to obtain two new SCUC models, and respectively carrying out iterative solution on the two new SCUC models;
when any one new SCUC model is successfully solved, carrying out termination operation on the other new SCUC model, and taking a result obtained after the solution is successfully solved as an optimal mode;
determining SCED constraint conditions and an SCED optimization target according to the optimal mode, establishing an SCED model, solving the SCED model through a dual simplex method to obtain Lagrangian multiplier values corresponding to constraints of each unit in different time periods, and calculating the electrovalence of each node in different time periods according to the Lagrangian multiplier values corresponding to the constraints in different time periods.
2. The rapid a posteriori conditional section electricity market clearing method of claim 1 wherein said SCUC optimization objective is:
Figure FDA0002417919180000011
wherein: n represents the total number of the units; t represents the total number of time segments considered;
Pi,tto representThe output of the unit i in the time period t;
Figure FDA0002417919180000012
the running cost, the starting cost and the shutdown cost of the unit i in the time period t are respectively, wherein the unit running cost Ci,t(Pi,t) Is a multi-segment linear function of the unit output.
3. The method for rapid a posteriori conditional cut-out of the electricity market of claim 1 wherein said predetermined constraints comprise: the method comprises the following steps of load balance constraint, unit output upper and lower limit constraint, unit climbing constraint, unit minimum continuous start-stop time constraint and unit maximum start-stop times constraint.
4. The method for rapidly clearing an electric power market of a posterior condition section according to claim 1, wherein the formula for calculating the section power flow by the direct current power flow method is as follows:
Figure FDA0002417919180000021
wherein G isl-iOutputting a power transfer distribution factor for a generator of a line or a section l by a node where a unit i is located; k is the number of nodes of the system; gl-kThe distribution factor of the output power transfer of the generator of the node k to the line or the section l; dk,tIs the bus load value of node k in the period t.
5. The method for rapid a posteriori profiling electricity market clearing according to claim 1 wherein said SCED constraints comprise: load balance constraint, unit output upper and lower limit constraint, unit climbing constraint and conditional section constraint.
6. The rapid a posteriori profiling electric power market clearing method according to claim 1 wherein said SCED optimization objectives are:
Figure FDA0002417919180000022
wherein: n represents the total number of the units; t represents the total number of time segments considered;
Pi,trepresenting the output of the unit i in the time period t; ci,t(Pi,t) The operating cost of the unit i in the time period t is a multi-segment linear function of the unit output.
7. The method for clearing the electric power market of the rapid posterior conditional section as recited in claim 1, wherein the calculation formula for calculating the electricity prices of the nodes in different periods according to the lagrangian multiplier values corresponding to the constraints in different periods is as follows:
Figure FDA0002417919180000023
wherein λ ist: lagrangian multipliers of system load balance constraint in the period t;
Figure FDA0002417919180000031
a lagrange multiplier for the maximum forward power flow constraint of the line l;
Figure FDA0002417919180000032
a lagrangian multiplier constrained by the maximum reverse power flow of the line l;
Figure FDA0002417919180000033
a lagrange multiplier with section s maximum forward power flow constraint;
Figure FDA0002417919180000034
a lagrange multiplier with section s maximum reverse power flow constraint; gl-k: the node k transfers the distribution factor to the generator output power of the line l; gs-k: node k to section s generator outputPower transfer profile factor.
8. The utility model provides a quick posterior conditional section electric power market goes out clear system which characterized in that includes:
the model establishing module is used for establishing an SCUC model according to preset constraint conditions and an SCUC optimization target, and solving the SCUC model through a branch-and-bound method to obtain the output and start-stop states of each unit at different time intervals;
the section power flow module is used for calculating the section power flow through a direct current power flow method according to the output and start-stop states of the units at different time intervals to obtain branch power flow and section power flow results;
the iteration solving module is used for judging the conditional section constraint or the control section constraint falling into the conditional section according to the branch flow and section flow results, respectively locking the effective conditional section, respectively adding the effective conditional section into the SCUC model to obtain two new SCUC models, and respectively carrying out iteration solving on the two new SCUC models;
the termination operation module is used for carrying out termination operation on another new SCUC model after any new SCUC model is successfully solved, and taking a result obtained after the solution is successfully solved as an optimal mode;
and the electricity price calculation module is used for determining the SCED constraint condition and the SCED optimization target according to the optimal mode, establishing an SCED model, solving the SCED model through a dual simplex method to obtain Lagrange multiplier values corresponding to the constraints of each unit in different time periods, and calculating the electricity price of each node in different time periods according to the Lagrange multiplier values corresponding to the constraints in different time periods.
9. A computer-readable storage medium, characterized in that the computer-readable storage medium comprises a stored computer program; wherein the computer program, when executed, controls an apparatus in which the computer-readable storage medium is located to perform the fast a posteriori conditional cut-off for electricity market clearing method according to any of claims 1 to 7.
10. A terminal device comprising a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, the processor implementing the fast a posteriori conditional cut-out for electricity market clearing method according to any of claims 1 to 7 when executing the computer program.
CN202010196782.6A 2020-03-19 2020-03-19 Electric power market clearing method and system for rapid posterior condition section Active CN111428922B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010196782.6A CN111428922B (en) 2020-03-19 2020-03-19 Electric power market clearing method and system for rapid posterior condition section

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010196782.6A CN111428922B (en) 2020-03-19 2020-03-19 Electric power market clearing method and system for rapid posterior condition section

Publications (2)

Publication Number Publication Date
CN111428922A CN111428922A (en) 2020-07-17
CN111428922B true CN111428922B (en) 2022-04-05

Family

ID=71546486

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010196782.6A Active CN111428922B (en) 2020-03-19 2020-03-19 Electric power market clearing method and system for rapid posterior condition section

Country Status (1)

Country Link
CN (1) CN111428922B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111950800A (en) * 2020-08-21 2020-11-17 中国电力科学研究院有限公司 Method, system, equipment and readable storage medium for market clearing before electric power spot goods day
CN112993974B (en) * 2020-12-26 2023-03-14 广东电网有限责任公司电力调度控制中心 Method, device, terminal and medium for calculating market discharge price before current electric power spot date
CN112906918B (en) * 2021-03-10 2023-04-07 广东电网有限责任公司电力调度控制中心 Method and device for checking start-stop rationality of clear result unit combination
CN113609439B (en) * 2021-06-17 2024-04-12 国电南瑞南京控制系统有限公司 Spot market clearing method and device considering clearing of power generation side with quotation proportion
CN115907140A (en) * 2022-11-18 2023-04-04 中国南方电网有限责任公司 Method and device for optimizing power spot shipment scheme, computer equipment and medium
CN116757388A (en) * 2023-04-20 2023-09-15 中国南方电网有限责任公司 Electric power market clearing method and device based on redundancy constraint screening

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101917024A (en) * 2010-03-29 2010-12-15 福建省电力有限公司 Generating method of universality cost space in security-constrained dispatch
CN105303256A (en) * 2015-10-29 2016-02-03 西安交通大学 Power inter-provincial and inter-district trade path analysis method
WO2018049737A1 (en) * 2016-09-18 2018-03-22 国电南瑞科技股份有限公司 Safe correction calculation method based on partition load control
CN109980631A (en) * 2017-12-28 2019-07-05 广东电网有限责任公司电力调度控制中心 A kind of power spot market a few days ago go out clearly with node electricity price calculation method
CN110061510A (en) * 2019-03-11 2019-07-26 中国电力科学研究院有限公司 A kind of period decoupling security constrained economic dispatch fast solution method and system
CN110555590A (en) * 2019-07-31 2019-12-10 云南电网有限责任公司 secondary clearing method for watershed cascade upstream and downstream power stations participating in electric power spot market
CN110766480A (en) * 2019-10-31 2020-02-07 中国南方电网有限责任公司 Real-time market clearing optimization method and device considering condition control section
CN110782281A (en) * 2019-10-23 2020-02-11 四川大学 Day-ahead market clearing method for multi-owner cascade power station basin electric quantity transfer
CN110874692A (en) * 2019-11-05 2020-03-10 中国南方电网有限责任公司 Alternating current-direct current hybrid regional power grid spot market clearing calculation optimization method and device

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100107173A1 (en) * 2008-09-29 2010-04-29 Battelle Memorial Institute Distributing resources in a market-based resource allocation system
US10991041B2 (en) * 2014-10-03 2021-04-27 Open Access Technology International, Inc. Next-generation energy market design and implementation

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101917024A (en) * 2010-03-29 2010-12-15 福建省电力有限公司 Generating method of universality cost space in security-constrained dispatch
CN105303256A (en) * 2015-10-29 2016-02-03 西安交通大学 Power inter-provincial and inter-district trade path analysis method
WO2018049737A1 (en) * 2016-09-18 2018-03-22 国电南瑞科技股份有限公司 Safe correction calculation method based on partition load control
CN109980631A (en) * 2017-12-28 2019-07-05 广东电网有限责任公司电力调度控制中心 A kind of power spot market a few days ago go out clearly with node electricity price calculation method
CN110061510A (en) * 2019-03-11 2019-07-26 中国电力科学研究院有限公司 A kind of period decoupling security constrained economic dispatch fast solution method and system
CN110555590A (en) * 2019-07-31 2019-12-10 云南电网有限责任公司 secondary clearing method for watershed cascade upstream and downstream power stations participating in electric power spot market
CN110782281A (en) * 2019-10-23 2020-02-11 四川大学 Day-ahead market clearing method for multi-owner cascade power station basin electric quantity transfer
CN110766480A (en) * 2019-10-31 2020-02-07 中国南方电网有限责任公司 Real-time market clearing optimization method and device considering condition control section
CN110874692A (en) * 2019-11-05 2020-03-10 中国南方电网有限责任公司 Alternating current-direct current hybrid regional power grid spot market clearing calculation optimization method and device

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A Hierarchical Method for Robust SCUC of Multi-Area Power Systems With Novel Uncertainty Sets;Xiaodong Zheng 等;《IEEE Transactions on Power Systems》;20190904;第35卷(第2期);第1364-1375页 *
An optimal pricing scheme in electricity markets by parallelizing security constrained optimal power flow based market-clearing model;M.K.Kim 等;《International Journal of Electrical Power & Energy Systems》;20130630;第48卷;第161-171页 *
区域互联电网电能量与备用辅助服务联合优化模型研究;周华锋 等;《电网技术（CNKI网络首发）》;20200120;第991-1001页 *
考虑条件控制断面的跨区高压直流发输电系统发电计划模型;徐帆 等;《电力系统自动化》;20190910;第43卷(第17期);第176-183页 *

Also Published As

Publication number Publication date
CN111428922A (en) 2020-07-17

Similar Documents

Publication Publication Date Title
CN111428922B (en) Electric power market clearing method and system for rapid posterior condition section
Wang et al. Robust co-planning of energy storage and transmission line with mixed integer recourse
Rezvani et al. Environmental/economic scheduling of a micro-grid with renewable energy resources
Ju et al. A bi-level stochastic scheduling optimization model for a virtual power plant connected to a wind–photovoltaic–energy storage system considering the uncertainty and demand response
Hesamzadeh et al. Computation of extremal-Nash equilibria in a wholesale power market using a single-stage MILP
CN107392395A (en) A kind of power distribution network and micro electric network coordination optimization method based on price competition mechanism
CN113541166B (en) Distributed energy storage optimal configuration method, system, terminal and storage medium
CN106844916A (en) A kind of generating and electricity consumption method for organizing and device based on peak regulation assisted hatching
CN111291985B (en) Day-ahead power market clearing method and device based on conditional section modeling
CN110232583B (en) Electric power market marginal price planning method considering carbon emission right
Maffei et al. Optimal Power Flow model with energy storage, an extension towards large integration of renewable energy sources.
CN115271397A (en) Calculation method and device for electric power day-ahead market clearing and storage medium
CN111667136A (en) Clearing method and device for regional power market and storage medium
CN115841187A (en) Method, device, equipment and storage medium for optimizing operation strategy of flexible power distribution network
Malysz et al. MILP-based rolling horizon control for microgrids with battery storage
CN113659576B (en) Cross-region interaction method, device, equipment and medium under alternating current-direct current hybrid power grid
Song et al. Bidding strategy for virtual power plants with the day-ahead and balancing markets using distributionally robust optimization approach
CN113538066B (en) Control method, system, equipment and medium for virtual power plant to participate in power market
Jia et al. Dynamic pricing strategy and regional energy consumption optimization based on different stakeholders
Yuan et al. Transmission expansion planning with optimal transmission switching considering uncertain nk contingency and renewables
KR102319625B1 (en) Microgrid operation method and device to provide operational flexibility in main grid operation
CN116402223A (en) Cooperative scheduling method, system and equipment for power distribution network
CN115760194A (en) Method, device and equipment for optimizing peak-valley electricity price and time interval and storage medium
CN113139271B (en) Aggregation resource modeling method and device, electronic device and storage medium
CN111404154A (en) Power distribution network power supply capacity optimization method, equipment and storage medium

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20231120

Address after: No.73, Tuodong Road, Kunming, Yunnan 650051

Patentee after: YUNNAN POWER GRID Co.,Ltd.

Address before: No.73, Tuodong Road, Kunming, Yunnan 650051

Patentee before: YUNNAN POWER GRID Co.,Ltd.

Patentee before: BEIJING TSINTERGY TECHNOLOGY CO.,LTD.