CN109301876B - Constraint condition elasticized electric power day-ahead market clearing method - Google Patents

Constraint condition elasticized electric power day-ahead market clearing method Download PDF

Info

Publication number
CN109301876B
CN109301876B CN201811025518.5A CN201811025518A CN109301876B CN 109301876 B CN109301876 B CN 109301876B CN 201811025518 A CN201811025518 A CN 201811025518A CN 109301876 B CN109301876 B CN 109301876B
Authority
CN
China
Prior art keywords
time period
renewable energy
thermal power
hydropower
reduction
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
CN201811025518.5A
Other languages
Chinese (zh)
Other versions
CN109301876A (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.)
Tsinghua University
State Grid Corp of China SGCC
State Grid Zhejiang Electric Power Co Ltd
Original Assignee
Tsinghua University
State Grid Corp of China SGCC
State Grid Zhejiang 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 Tsinghua University, State Grid Corp of China SGCC, State Grid Zhejiang Electric Power Co Ltd filed Critical Tsinghua University
Priority to CN201811025518.5A priority Critical patent/CN109301876B/en
Publication of CN109301876A publication Critical patent/CN109301876A/en
Application granted granted Critical
Publication of CN109301876B publication Critical patent/CN109301876B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • 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
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Public Health (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Water Supply & Treatment (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Power Engineering (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention relates to a constraint condition-elasticized electric power day-ahead market clearing method, which belongs to the technical field of electric power system scheduling and electric power market trading. The method provides a thermal power generating unit and renewable energy modeling method with elasticated constraint conditions, a modeling method for balancing thermal power deep peak shaving and renewable energy source reduction, and a model for uniformly reducing renewable energy sources, provides technical support for the daily market clearance or daily planning of the elasticated constraint conditions considering thermal power deep peak shaving and renewable energy source orderly reduction, solves the problem of coordination and mutual arrangement of thermal power and renewable energy sources in the daily market, and promotes the consumption of the renewable energy sources.

Description

Constraint condition elasticized electric power day-ahead market clearing method
Technical Field
The invention discloses a constraint condition-elasticized electric power day-ahead market clearing method, and belongs to the technical field of electric power system scheduling and electric power market trading.
Background
In recent years, the renewable development of China is rapid, the installed capacity is continuously improved, but the problem of 'three abandons' of renewable energy sources is more and more prominent, so that the energy conservation and emission reduction of China are seriously restricted, and the serious economic loss is caused. In the traditional power system scheduling, a thermal power generating unit can only perform scheduling between the minimum economic output and installed capacity (Xuqingshan, Dingyifan, Jiyili. a power grid safety optimization scheduling method considering demand response [ P ]. Jiangsu: CN106712005A,2017-05-24.) (Huangqiang, Mengnan wave. a power system scheduling method [ P ]. Guangdong: CN107294126A,2017-10-24.), so that the adjustable space of the system is severely limited.
In order to solve the problem of serious shortage of peak-load-reduction resources of a system, the peak-load-reduction potential of the system is excavated in a market mode, the consumption difficulty of renewable energy sources is relieved, the constraint condition of power system scheduling needs to be elasticized, thermal power deep peak-load-reduction and renewable energy source reduction are combined in the market in the day ahead, thermal power and renewable energy source resources are jointly scheduled, and the problem of coordination and mutual arrangement of thermal power and renewable energy sources in the market in the day ahead is solved. Currently, reports related to peak shaving resource calling for solving renewable energy grid connection (Yuanbo, Zhang jin Fang, Wang Guanghua, Wang Xiao Cheng, Junan, Zhang Fu, Feng Junshu, Zheng Guang, Jiabing Qi, Yan Xiao Qing, Wu Yu, Wang gang, Zhang Jiang Bo, Liujun) exist, but according to the cost and the adjustable capacity of various peak shaving resources, the decision method for peak shaving resource calling is suitable for large-scale renewable energy grid connection [ P ]. Beijing: CN107276122A,2017-10-20 ], but only for various peak shaving resources, the peak-shaving resource calls are sequenced, but the day-ahead planning arrangement and scheduling of the power system are not realized under the condition of considering deep peak shaving, various resources are not modeled in a fine mode, the resources are not subjected to combined optimization in all time periods all day according to the physical characteristics of the resources, and the optimal system and the maximum consumption of renewable energy cannot be realized.
In view of the above, there is no report on the electric power day-ahead market clearing method with the constraint condition elasticization.
Disclosure of Invention
The invention aims to fill the technical blank in the field of current power system scheduling and power market trading, and provides a constraint condition-elasticated power day-ahead market clearing method. The invention provides a thermal power generating unit with an elasticized constraint condition and a renewable energy modeling method, provides a method for linearizing the output of hydropower, provides a modeling method for balancing thermal power deep peak shaving and renewable energy reduction, provides a model for uniformly reducing the renewable energy, can provide technical support for the daily market clearing or daily planning of the elasticized constraint condition considering thermal power deep peak shaving and renewable energy ordered reduction, solves the problem of coordination and mutual arrangement of thermal power and renewable energy in the daily market, promotes the consumption of the renewable energy, and improves the safety of the trading result of the electric power market.
The invention provides a constraint condition elasticized electric power day-ahead market clearing method, which comprises the following steps of:
1) constructing a constraint condition elasticized mathematical model of the clear electric power market at the day before, wherein the involved energy sources are in two forms, namely thermal power and renewable energy sources, the renewable energy sources comprise an uncontrollable renewable energy source and an adjustable hydropower form, the uncontrollable renewable energy sources are respectively radial flow hydropower, wind power and photovoltaic, and the specific process is as follows;
1-1) constructing an objective function of a constraint condition-elasticated electric power day-ahead market clearing mathematical model as follows:
Figure GDA0002952840500000021
the objective function represents the minimum total operating cost of the power system at the maximum renewable energy consumption, where the upper scale th represents thermal power, the upper scale r represents uncontrollable renewable energy, the upper scale h represents regulated hydropower,
Figure GDA0002952840500000022
the price of the k nth energy quotation section of the thermal power generating unit is the price of the electricity supplied to the thermal power generating unit participating in the market at the day beforeThe price of the unit kilowatt-hour electric quantity corresponding to the nth energy quotation section of the second day reported by the trading institution,
Figure GDA0002952840500000023
the nth energy quotation section of the thermal power generating unit is subjected to the medium-winning amount in the time period t,
Figure GDA0002952840500000024
the starting cost of the thermal power generating unit k in the time period t is shown, M is a weight coefficient, and the value range of M is
Figure GDA0002952840500000025
Figure GDA0002952840500000026
Avoiding reducing the price corresponding to the reduction-avoiding claim quantity for the uncontrollable renewable energy k in the nth section of the time period t and reporting the uncontrollable renewable energy participating in the market to the electric power exchange mechanism in the second day section,
Figure GDA0002952840500000027
to avoid the reduction of the reduction declaration quantity in the nth period of the time t for the uncontrollable renewable energy source k,
Figure GDA0002952840500000028
for avoiding reduction of quotations of the regulated hydropower k in the nth section of the time period t, for avoiding reduction of prices corresponding to reduction of reporting quantity in the second day section by section for reporting the regulated hydropower participating in the market in the day ahead to the electric power transaction institution,
Figure GDA0002952840500000029
to avoid the reduction of the declaration quantity in the nth period of the time t,
Figure GDA00029528405000000210
giving out electricity for the thermal power generating unit k in the nth section of the time interval t, and letting out electricity for the thermal power generating unit k participating in the market in the day ahead in the kilowatt-hour unit corresponding to the second day subsection deep peak regulation capacity reported to the electric power trading mechanism by the thermal power generating unitThe price of the amount is such that,
Figure GDA00029528405000000211
the method comprises the steps of obtaining the quantity of the thermal power generating unit k for the nth section depth peak regulation in a time period t;
1-2) constructing constraint conditions of a constraint condition-elasticized electric power day-ahead market clearing mathematical model, wherein expressions are respectively as follows:
1-2-1) power system constraints:
Figure GDA0002952840500000031
Figure GDA0002952840500000032
wherein the content of the first and second substances,
Figure GDA0002952840500000033
indicating that for all of the (-) s that are true,
Figure GDA0002952840500000034
the output of the uncontrollable renewable energy k during the time period t,
Figure GDA0002952840500000035
for regulating the output of the hydropower k over a period of time t, DtFor the total load of the power system in the time period t, the predicted value P of the total load demand of the power system in each time period of the next day is obtained by the power dispatching mechanism by adopting a prediction methodLmaxThe maximum active transmission power of the line L and the maximum active transmission power of each line on the second day given by the dispatching mechanism, G is a generator output power transfer distribution factor matrix of the power system, generator output power transfer distribution factor data of each power transmission line by each node given by the dispatching mechanism, and Gr-k,LThe distribution factor of the output power transfer of the generator from the node where the uncontrollable renewable energy k is located to the line L, r in the subscript represents the uncontrollable renewable energy, Gh-k,LGenerator for regulating the node of hydroelectric k to line LOutput power transfer distribution factor, h in subscript denotes regulated hydropower, Gth-k,LFor a generator output power transfer distribution factor from a node where a thermal power generating unit k is located to a line L, th in subscript represents thermal power, Gb,LFor generator output power transfer distribution factor, D, from node b to line Lb,tFor the load of the node B in the time period t, obtaining a predicted value of the power load demand of each node in each time period of the second day by a power dispatching mechanism by adopting a prediction method, wherein B is the total number of nodes and the total number of power network nodes given by the dispatching mechanism;
the formula (2) is a load balance constraint of the power system, and represents that the total power generation amount of the power system in each period is equal to the total load of the power system;
formula (3) is network constraint, and indicates that the active transmission power of each line must not exceed the active transmission capacity;
1-2-2) uncontrollable renewable energy constraint conditions:
Figure GDA0002952840500000036
Figure GDA0002952840500000037
wherein the content of the first and second substances,
Figure GDA0002952840500000038
predicted output for the uncontrollable renewable energy k in the time period t, predicted output value for each time period of the second day predicted for each uncontrollable renewable energy,
Figure GDA0002952840500000039
the reduction of the reporting amount is avoided for the uncontrollable renewable energy k in the nth period of the time period t, and the reduction of the reporting amount is avoided for the uncontrollable renewable energy participating in the market in the second day and reported to the electric power exchange mechanism in sections; the constraint condition group indicates that the sum of the winning output and the reduction of the uncontrollable renewable energy k in the time period t is equal to the predicted output of the uncontrollable renewable energy k in the time period t, and the reduction of the declaration quantity is avoided in each time periodThe decrement should not be higher than the reduction declaration amount of the renewable energy source in the period;
1-2-3) regulating type hydropower linearization constraint condition
Figure GDA0002952840500000041
Figure GDA0002952840500000042
Figure GDA0002952840500000043
Figure GDA0002952840500000044
Figure GDA0002952840500000045
Figure GDA0002952840500000046
Figure GDA0002952840500000047
Figure GDA0002952840500000048
Wherein the content of the first and second substances,
Figure GDA0002952840500000049
for the lower limit of the applicable water storage capacity range of the regulated hydropower k c output curve,
Figure GDA00029528405000000410
water storage capacity suitable for k (th) power curve of adjustable hydropowerThe upper end of the range is,
Figure GDA00029528405000000411
to regulate the amount of water stored by the hydro-electric k over a period of time t,
Figure GDA00029528405000000412
to represent integer variables applicable to the linear section, if the regulated hydropower k is in the w-th section water discharge amount range of the c-th output curve in the time period t, then
Figure GDA00029528405000000413
If the adjusting type hydropower k is not in the w-th section water discharge amount range of the c-th output curve in the time period t, the water discharge amount of the adjusting type hydropower k is adjusted
Figure GDA00029528405000000414
For regulating the lower limit of the water discharge quantity range of the w section of the c output curve of the hydropower k,
Figure GDA00029528405000000415
the water discharge upper limit of the w-th section of the water discharge range of the c-th output curve of the adjustable hydropower k,
Figure GDA00029528405000000416
to adjust the discharge of the hydropower k at the time period t,
Figure GDA00029528405000000417
in order to adjust the position proportion of the water discharge quantity of the hydropower k in the w-th section water discharge quantity range of the c-th output curve,
Figure GDA00029528405000000418
in order to adjust the output of the water discharge quantity of the hydropower k at the w-th section water discharge quantity lower limit of the c-th output curve,
Figure GDA00029528405000000419
in order to adjust the output of the water discharge quantity of the hydropower k at the w-th section water discharge quantity upper limit of the c-th output curve,
Figure GDA00029528405000000420
Figure GDA00029528405000000421
and
Figure GDA00029528405000000422
the output of the hydroelectric generating set corresponding to different water discharge amounts in different water storage amount sections reported to the dispatching mechanism is a piecewise linear function;
the restriction condition set is obtained by carrying out linear treatment on the nonlinear and non-convex output curve of the hydropower so as to approximately obtain the hydropower output
Figure GDA00029528405000000423
Water discharge amount with water and electricity
Figure GDA00029528405000000424
And the water storage capacity
Figure GDA00029528405000000425
The relationship between;
1-2-4) regulating hydropower constraint conditions:
Figure GDA00029528405000000426
Figure GDA0002952840500000051
Figure GDA0002952840500000052
Figure GDA0002952840500000053
Figure GDA0002952840500000054
Figure GDA0002952840500000055
wherein the content of the first and second substances,
Figure GDA0002952840500000056
and
Figure GDA0002952840500000057
respectively the minimum and maximum water storage capacity of the adjustable hydroelectric power k, the maximum/minimum water storage capacity constraint reported by the hydroelectric generating set,
Figure GDA0002952840500000058
for the predicted inflow of the regulated hydropower k at time period t, for the second water-day prediction data reported by the hydroelectric generating set to the scheduling authority,
Figure GDA0002952840500000059
to adjust the amount of waste water of the hydro-power k during the time period t,
Figure GDA00029528405000000510
the ith hydroelectric station of basin p,
Figure GDA00029528405000000511
for hydropower stations
Figure GDA00029528405000000512
In a period of time
Figure GDA00029528405000000513
The amount of the supplied water is controlled by the amount of the supplied water,
Figure GDA00029528405000000514
hydropower station
Figure GDA00029528405000000515
In a period of time
Figure GDA00029528405000000516
The amount of natural water coming in is predicted,
Figure GDA00029528405000000517
for hydropower stations
Figure GDA00029528405000000518
At the amount of water discharge for the period t,
Figure GDA00029528405000000519
for hydropower stations
Figure GDA00029528405000000520
The amount of reject water at the time period t,
Figure GDA00029528405000000521
for water flow from hydropower stations
Figure GDA00029528405000000522
Flows to downstream hydropower stations
Figure GDA00029528405000000523
The time of (a) is,
Figure GDA00029528405000000524
and
Figure GDA00029528405000000525
for the step hydropower station data acquired by the scheduling authority,
Figure GDA00029528405000000526
the reporting amount is prevented from being reduced for the regulating hydropower k in the nth section of the time period t, and the reporting amount is prevented from being reduced for the regulating hydropower k participating in the day-ahead market in the second day section;
the constraint condition group reflects the operation constraint of the adjustable hydropower station, wherein a formula (13) is a hydropower station water storage quantity constraint, a formula (14) is a reservoir dynamic balance equation constraint, a formula (15) is an upstream and downstream hydropower station coupling constraint, a formula (16) is an adjustable hydropower station output constraint, a formula (17) is an adjustable hydropower station output constraint if water is not abandoned, and a formula (18) is a reduction constraint for avoiding the reduction of the declared quantity in each section;
1-2-5) constraint conditions of the thermal power generating unit:
Figure GDA00029528405000000527
Figure GDA00029528405000000528
Figure GDA00029528405000000529
Figure GDA00029528405000000530
Figure GDA00029528405000000531
Figure GDA0002952840500000061
Figure GDA0002952840500000062
Figure GDA0002952840500000063
Figure GDA0002952840500000064
Figure GDA0002952840500000065
Figure GDA0002952840500000066
wherein the content of the first and second substances,
Figure GDA0002952840500000067
the minimum output of the thermal power generating unit k is the minimum economic output reported by the thermal power generating unit to the power dispatching mechanism if the thermal power generating unit is not the necessary starting group, the minimum economic output is the output of the dispatching mechanism for decomposing the physical contract electric quantity of the thermal power generating unit to each time period of the second day if the thermal power generating unit is the necessary starting group,
Figure GDA0002952840500000068
the unit capacity of the thermal power generating unit k is the unit capacity reported to the power dispatching mechanism by the thermal power generating unit,
Figure GDA0002952840500000069
is an integer variable of the starting state of the thermal power generating unit, if the thermal power generating unit k is in the starting state,
Figure GDA00029528405000000610
if the thermal power generating unit k is not in the on state,
Figure GDA00029528405000000611
the maximum climbing rate of the thermal power generating unit k and the maximum climbing rate constraint of the thermal power generating unit reported to the dispatching mechanism,
Figure GDA00029528405000000612
starting an integer variable for the thermal power generating unit, if the thermal power generating unit k is started in a time period t,
Figure GDA00029528405000000613
if the thermal power generating unit k is not started in the time period t,
Figure GDA00029528405000000614
is an integer variable of the shutdown of the thermal power unit, if the thermal power unit k is in a time intervalt is turned off and stopped, and t is turned off,
Figure GDA00029528405000000615
if the thermal power generating unit k is not shut down in the time period t,
Figure GDA00029528405000000616
Figure GDA00029528405000000617
for the minimum on-state duration of the thermal power generating unit k,
Figure GDA00029528405000000618
is the minimum shutdown state duration of the thermal power generating unit k,
Figure GDA00029528405000000619
the maximum starting times of the thermal power generating unit k in one day,
Figure GDA00029528405000000620
in order to increase the starting cost of the thermal power generating unit k,
Figure GDA00029528405000000621
and
Figure GDA00029528405000000622
the minimum startup state duration, the minimum shutdown state duration, the maximum startup times in a day and the startup cost data of the thermal power unit are reported to a scheduling mechanism by the thermal power unit,
Figure GDA00029528405000000623
determining information whether the thermal power unit is a necessary unit for a scheduling mechanism for necessary constraint parameters of the thermal power unit k, if the thermal power unit k is the necessary unit,
Figure GDA00029528405000000624
if the thermal power generating unit k does not need to be started,
Figure GDA00029528405000000625
is a fireThe output declaration value of the nth energy quotation section of the generator set k is the second day subsection output reported to the electric power trading mechanism by the thermal power generating units participating in the market in the day ahead;
the constraint condition group reflects the operation constraint of the thermal power unit, formulas (19) and (20) are the output range constraint of the thermal power unit, the deep peak shaving of the thermal power unit is embodied in the formula (19), the climbing rate constraint of the thermal power unit is embodied in the formula (21), the formula (22) is the minimum starting state duration constraint of the thermal power unit, the formula (23) is the minimum stopping state duration constraint of the thermal power unit, the formula (24) is the maximum starting time constraint of the thermal power unit, the formulas (25) and (26) are the relation constraint between integer variables, the formula (27) is the starting cost constraint of the thermal power unit, the formula (28) is the necessary starting group constraint, and the formula (29) is the medium output constraint of each energy quotation section of the thermal power unit;
1-3) constructing a mathematical model of the daily market clearing of the uniformly reduced electric power without considering constraint condition elasticization, wherein the expressions are respectively as follows:
1-3-1) objective function of the electric power day-ahead market clearing mathematical model with uniform curtailment elasticated without consideration of constraints:
Figure GDA0002952840500000071
wherein the content of the first and second substances,
Figure GDA0002952840500000072
and
Figure GDA0002952840500000073
are respectively
Figure GDA0002952840500000074
And
Figure GDA0002952840500000075
the weight coefficient of (a) is,
Figure GDA0002952840500000076
and as n is increased,
Figure GDA0002952840500000077
and
Figure GDA0002952840500000078
the size of the composite material is also increased continuously,
Figure GDA0002952840500000079
for the uniform reduction of the uncontrollable renewable energy k in the nth reduction segment of the time period t,
Figure GDA00029528405000000710
the uniform reduction amount of the adjustable hydropower k in the nth section reduction section of the time period t is reduced;
the objective function shows that when the power system can not consume the renewable energy, the renewable energy is evenly reduced;
1-3-2) variable constraints of the electric power day-ahead market clearing mathematical model without consideration of constraint condition elasticization:
Figure GDA00029528405000000711
Figure GDA00029528405000000712
n is large enough to make uncontrollable renewable energy k output predicted in time period t
Figure GDA00029528405000000713
And adjusting the power of the hydropower k in the time period t
Figure GDA00029528405000000714
Can be divided small enough so that renewable energy can be evenly curtailed;
1-3-3) the constraint conditions of the power system of the electric power daily market clearing mathematical model which do not consider the constraint conditions for the elasticization and the uniform reduction are the same as all the constraint conditions of the step 1-2-1);
1-3-4) the constraint conditions of uncontrollable renewable energy sources such as radial flow type hydropower, wind power, photovoltaic and the like of a mathematical model of the market clearing at the day before the electric power with the constraint conditions of the elasticization and the uniform reduction are not considered, and the constraint conditions are the same as all the constraint conditions of the step 1-2-2);
1-3-5) adjusting type hydropower linearization constraint conditions of a power day-ahead market clearing mathematical model, which do not consider constraint conditions for elasticization and uniform reduction, are the same as all the constraint conditions of the step 1-2-3);
1-3-6) considering the adjustable hydropower constraint condition of the power day-ahead market clearing mathematical model with the constraint condition elasticization and uniform reduction, and adding the following constraint conditions on the basis of all the constraint conditions of the step 1-2-4):
Figure GDA00029528405000000715
wherein the content of the first and second substances,
Figure GDA0002952840500000081
calculating the water discharge quantity of the adjustable hydropower k in the time period t in the constraint condition elasticized day-ahead market in the step 1-2),
Figure GDA0002952840500000082
calculating the water abandon amount of the adjustable hydropower k in the time period t in the constraint condition elasticized day-ahead market in the step 1-2);
1-3-7) eliminating constraint conditions and elasticizing the constraint conditions, and uniformly reducing the constraint conditions to obtain the constraint conditions of the thermal power generating unit of the mathematical model in the market before the day:
Figure GDA0002952840500000083
wherein the content of the first and second substances,
Figure GDA0002952840500000084
and
Figure GDA0002952840500000085
respectively calculating results of the output and the depth peak regulation amount of the fire generator set k in the time period t in the constraint condition elasticized day-ahead market in the step 1-2);
the constraint condition ensures that the thermal power generating unit does not carry out deep peak regulation intentionally any more, and the output is fixed on the clear result or the minimum output of the market in the day before the constraint condition is elasticized;
2) solving the constraint condition elasticized electric power day-ahead market clearing mathematical model in the step 1-2) by adopting a mixed integer programming solving method to obtain constraint condition elasticized electric power day-ahead market clearing results, wherein the clearing results comprise the medium-output power of each thermal power generating unit in each quotation section in each time period
Figure GDA0002952840500000086
Output of each uncontrollable renewable energy source in each time period
Figure GDA0002952840500000087
The output of each adjustable water and electricity in each time interval
Figure GDA0002952840500000088
Reduction amount of each uncontrollable renewable energy in each quotation section of each time period
Figure GDA0002952840500000089
Reduction of uncontrollable renewable energy sources in each time period
Figure GDA00029528405000000810
Reduction of each adjusting hydropower in each quotation section of each time period
Figure GDA00029528405000000811
Cutting amount of each adjustable water and electricity in each time period
Figure GDA00029528405000000812
Mid-bid depth peak regulation amount of each thermal power generating unit in each quotation section of each time period
Figure GDA00029528405000000813
3) Solving the uniform reduction daily market clearing mathematical model which does not consider the elasticization of the constraint conditions and is obtained by adopting a mixed integer programming solving method, wherein the uniform reduction daily market clearing mathematical model which does not consider the elasticization of the constraint conditions in the step 1-3) is obtained, and the clearing result comprises the uniform reduction amount of each uncontrollable renewable energy source in each reduction section in each time period
Figure GDA00029528405000000814
Even reduction of uncontrollable renewable energy sources in each time period
Figure GDA00029528405000000815
The even reduction amount of each reduction section of each adjusting type water and electricity in each time period
Figure GDA00029528405000000816
The even reduction of each adjustable water and electricity in each time interval
Figure GDA00029528405000000817
4) Quoting each thermal power generating unit in each time interval according to the step 1) in each period
Figure GDA00029528405000000818
And the peak regulation amount of the thermal power generating units in the step 2) is the mid-bid depth of each quotation section in each period
Figure GDA00029528405000000819
Taking the highest price of the deep peak shaving in the deep peak shaving price quotation section of the winning bid of all the thermal power units in the time period as the deep peak shaving clear price;
5) calculating renewable energy avoids cutting clearing prices, including:
5-1) the reduction amount of each uncontrollable renewable energy source in each time period according to the step 2) above
Figure GDA0002952840500000091
Cutting amount of each adjustable water and electricity in each time period
Figure GDA0002952840500000092
Even reduction of uncontrollable renewable energy sources in each time period
Figure GDA0002952840500000093
And the even reduction amount of each adjustable hydropower in each time period in the step 3) is reduced
Figure GDA0002952840500000094
Subtracting the reduction amount of each renewable energy source in each time period in the step 2) from the uniform reduction amount of each renewable energy source in each time period in the step 3) to obtain the amount of each renewable energy source which is not discarded in each time period;
5-2) avoiding cutting off quotations at each stage of each time interval according to each uncontrollable renewable energy in the step 1) above
Figure GDA0002952840500000095
Avoiding reduction of quotations at each section of each time period by each adjusting type hydropower in step 1)
Figure GDA0002952840500000096
The reduction amount of each uncontrollable renewable energy source obtained in the step 2) in each time period
Figure GDA0002952840500000097
The reduction amount of each adjusting type hydropower obtained in the step 2) in each time period
Figure GDA0002952840500000098
Judging the amount of each renewable energy source to be abandoned in the period of time according to the amount of each renewable energy source to be abandoned obtained in the step 5-1), if the amount of each renewable energy source to be abandoned in the period of time is positive, taking the avoidance reduction quotation corresponding to the reduction amount in the period of time obtained in the step 2) of the renewable energy source as the avoidance reduction quotation of the renewable energy source acting in the period of time, and if the amount of each renewable energy source to be abandoned in the period of time is not positive, avoiding the reduction of the renewable energy source not acting in the period of timeQuoting;
5-3) according to the avoidance reduction quotation of each renewable energy source acting in each time interval obtained in the step 5-2), taking the lowest quotation of the avoidance reduction quotation of each renewable energy source acting in each time interval as the avoidance reduction clearing price of the power system in the time interval, and realizing the clearing of the power market before the day with the elastic constraint condition.
The method for clearing the electric power market in the day before with the constraint condition elasticization, provided by the invention, has the advantages that:
the thermal power generating unit with the elasticized constraint condition and the renewable energy modeling method provide a method for linearizing the output of the hydropower, a modeling method for balancing thermal power deep peak shaving and renewable energy reduction, and a model for uniformly reducing the renewable energy, can provide technical support for the market clearing or planning before the day of the constraint condition elasticized thermal power deep peak shaving and the renewable energy ordered reduction, solve the problem of coordination and mutual arrangement of the thermal power and the renewable energy in the market before the day, promote the consumption of the renewable energy, improve the safety of the trading result of the electric power market, and have great significance for improving the dispatching level of an electric power system and the development of the electric power market.
Drawings
Fig. 1 is a block flow diagram of a method for market clearing of electric power in the day ahead, in which constraint conditions are elasticated, according to the present invention.
FIG. 2 is a linear approximate schematic diagram of hydroelectric power output involved in the method of the present invention.
Detailed Description
The invention provides a constraint condition-elasticated electric power day-ahead market clearing method, a flow chart of which is shown in figure 1, and the method comprises the following steps:
1) constructing a constraint condition elasticized mathematical model of the clear electric power market at the day before, wherein the involved energy sources are in two forms, namely thermal power and renewable energy sources, the renewable energy sources comprise an uncontrollable renewable energy source and an adjustable hydropower form, the uncontrollable renewable energy sources are respectively radial flow hydropower, wind power and photovoltaic, and the specific process is as follows;
1-1) constructing an objective function of a constraint condition-elasticated electric power day-ahead market clearing mathematical model as follows:
Figure GDA0002952840500000101
the objective function represents the minimum total operating cost of the power system at the maximum renewable energy consumption, where the upper scale th represents thermal power, the upper scale r represents uncontrollable renewable energy, the upper scale h represents regulated hydropower,
Figure GDA0002952840500000102
the price of the nth energy quotation section of the thermal power generating unit k is the price of the unit kilowatt-hour electric quantity corresponding to the nth energy quotation section of the second day reported to the electric power trading mechanism by the thermal power generating unit participating in the market in the day ahead,
Figure GDA0002952840500000103
the nth energy quotation section of the thermal power generating unit is subjected to the medium-winning amount in the time period t,
Figure GDA0002952840500000104
for the starting cost of the thermal power generating unit k in a time period t, M is a weight coefficient to ensure that renewable energy is preferentially consumed and deep peak regulation is preferentially not used, and the value range of M is
Figure GDA0002952840500000105
An embodiment of the invention takes the value 1000000000,
Figure GDA0002952840500000106
avoiding reducing the price corresponding to the reduction-avoiding claim quantity for the uncontrollable renewable energy k in the nth section of the time period t and reporting the uncontrollable renewable energy participating in the market to the electric power exchange mechanism in the second day section,
Figure GDA0002952840500000107
for uncontrollable renewable energy k in time periodsthe nth stage of t avoids the reduction of the declaration amount,
Figure GDA0002952840500000108
for avoiding reduction of quotations of the regulated hydropower k in the nth section of the time period t, for avoiding reduction of prices corresponding to reduction of reporting quantity in the second day section by section for reporting the regulated hydropower participating in the market in the day ahead to the electric power transaction institution,
Figure GDA0002952840500000109
to avoid the reduction of the declaration quantity in the nth period of the time t,
Figure GDA00029528405000001010
giving a thermal power generating unit k a deep peak regulation price in the nth section of the time interval t, giving out the price of electric quantity in kilowatt-hours corresponding to the second-day segmented deep peak regulation capacity reported to the electric power trading mechanism by the thermal power generating unit participating in the market in the day ahead,
Figure GDA00029528405000001011
the method comprises the steps of obtaining the quantity of the thermal power generating unit k for the nth section depth peak regulation in a time period t;
1-2) constructing constraint conditions of a constraint condition-elasticized electric power day-ahead market clearing mathematical model, wherein expressions are respectively as follows:
1-2-1) power system constraints:
Figure GDA0002952840500000111
Figure GDA0002952840500000112
wherein the content of the first and second substances,
Figure GDA0002952840500000113
indicating that for all of the (-) s that are true,
Figure GDA0002952840500000114
the output of the uncontrollable renewable energy k during the time period t,
Figure GDA0002952840500000115
for regulating the output of the hydropower k over a period of time t, DtFor the total load of the power system in the time period t, the predicted value P of the total load demand of the power system in each time period of the next day is obtained by the power dispatching mechanism by adopting a prediction methodLmaxThe maximum active transmission power of the line L and the maximum active transmission power of each line on the second day given by the dispatching mechanism, G is a generator output power transfer distribution factor matrix of the power system, generator output power transfer distribution factor data of each power transmission line by each node given by the dispatching mechanism, and Gr-k,LThe distribution factor of the output power transfer of the generator from the node where the uncontrollable renewable energy k is located to the line L, r in the subscript represents the uncontrollable renewable energy, Gh-k,LFor the distribution factor of the generator output power transfer from the node where the regulated hydropower k is located to the line L, h in the subscript denotes the regulated hydropower, Gth-k,LFor a generator output power transfer distribution factor from a node where a thermal power generating unit k is located to a line L, th in subscript represents thermal power, Gb,LFor generator output power transfer distribution factor, D, from node b to line Lb,tFor the load of the node B in the time period t, obtaining a predicted value of the power load demand of each node in each time period of the second day by a power dispatching mechanism by adopting a prediction method, wherein B is the total number of nodes and the total number of power network nodes given by the dispatching mechanism;
the formula (2) is a load balance constraint of the power system, and represents that the total power generation amount of the power system in each period is equal to the total load of the power system;
formula (3) is network constraint, and indicates that the active transmission power of each line must not exceed the active transmission capacity;
1-2-2) uncontrollable renewable energy constraint conditions:
Figure GDA0002952840500000116
Figure GDA0002952840500000117
wherein the content of the first and second substances,
Figure GDA0002952840500000118
predicted output for the uncontrollable renewable energy k in the time period t, predicted output value for each time period of the second day predicted for each uncontrollable renewable energy,
Figure GDA0002952840500000119
the reduction of the reporting amount is avoided for the uncontrollable renewable energy k in the nth period of the time period t, and the reduction of the reporting amount is avoided for the uncontrollable renewable energy participating in the market in the second day and reported to the electric power exchange mechanism in sections; the constraint condition group indicates that the sum of the winning output and the reduction of the uncontrollable renewable energy k in the time period t is equal to the predicted output of the uncontrollable renewable energy k in the time period t, and the reduction of the avoidance reduction declaration quantity of each time period is not higher than the avoidance reduction declaration quantity of the renewable energy in the time period;
1-2-3) regulating type hydropower linearization constraint condition
Figure GDA0002952840500000121
Figure GDA0002952840500000122
Figure GDA0002952840500000123
Figure GDA0002952840500000124
Figure GDA0002952840500000125
Figure GDA0002952840500000126
Figure GDA0002952840500000127
Figure GDA0002952840500000128
Wherein the content of the first and second substances,
Figure GDA0002952840500000129
for the lower limit of the applicable water storage capacity range of the regulated hydropower k c output curve,
Figure GDA00029528405000001210
for the upper limit of the applicable water storage capacity range of the c output curve of the regulating hydropower k,
Figure GDA00029528405000001211
to regulate the amount of water stored by the hydro-electric k over a period of time t,
Figure GDA00029528405000001212
to represent integer variables applicable to the linear section, if the regulated hydropower k is in the w-th section water discharge amount range of the c-th output curve in the time period t, then
Figure GDA00029528405000001213
If the adjusting type hydropower k is not in the w-th section water discharge amount range of the c-th output curve in the time period t, the water discharge amount of the adjusting type hydropower k is adjusted
Figure GDA00029528405000001214
For regulating the lower limit of the water discharge quantity range of the w section of the c output curve of the hydropower k,
Figure GDA00029528405000001215
the water discharge upper limit of the w-th section of the water discharge range of the c-th output curve of the adjustable hydropower k,
Figure GDA00029528405000001216
to adjust the discharge of the hydropower k at the time period t,
Figure GDA00029528405000001217
in order to adjust the position proportion of the water discharge quantity of the hydropower k in the w-th section water discharge quantity range of the c-th output curve,
Figure GDA00029528405000001218
in order to adjust the output of the water discharge quantity of the hydropower k at the w-th section water discharge quantity lower limit of the c-th output curve,
Figure GDA00029528405000001219
in order to adjust the output of the water discharge quantity of the hydropower k at the w-th section water discharge quantity upper limit of the c-th output curve,
Figure GDA00029528405000001220
Figure GDA00029528405000001221
and
Figure GDA00029528405000001222
the output of the hydroelectric generating set corresponding to different water discharge amounts in different water storage amount sections reported to the dispatching mechanism is a piecewise linear function;
the restriction condition set is obtained by carrying out linear treatment on the nonlinear and non-convex output curve of the hydropower so as to approximately obtain the hydropower output
Figure GDA00029528405000001223
Water discharge amount with water and electricity
Figure GDA00029528405000001224
And the water storage capacity
Figure GDA00029528405000001225
The relationship between;
1-2-4) regulating hydropower constraint conditions:
Figure GDA0002952840500000131
Figure GDA0002952840500000132
Figure GDA0002952840500000133
Figure GDA0002952840500000134
Figure GDA0002952840500000135
Figure GDA0002952840500000136
wherein the content of the first and second substances,
Figure GDA0002952840500000137
and
Figure GDA0002952840500000138
respectively the minimum and maximum water storage capacity of the adjustable hydroelectric power k, the maximum/minimum water storage capacity constraint reported by the hydroelectric generating set,
Figure GDA0002952840500000139
for the predicted inflow of the regulated hydropower k at time period t, for the second water-day prediction data reported by the hydroelectric generating set to the scheduling authority,
Figure GDA00029528405000001310
to adjust the amount of waste water of the hydro-power k during the time period t,
Figure GDA00029528405000001311
the ith hydroelectric station of basin p,
Figure GDA00029528405000001312
for hydropower stations
Figure GDA00029528405000001313
In a period of time
Figure GDA00029528405000001314
The amount of the supplied water is controlled by the amount of the supplied water,
Figure GDA00029528405000001315
hydropower station
Figure GDA00029528405000001316
In a period of time
Figure GDA00029528405000001317
The amount of natural water coming in is predicted,
Figure GDA00029528405000001318
for hydropower stations
Figure GDA00029528405000001319
At the amount of water discharge for the period t,
Figure GDA00029528405000001320
for hydropower stations
Figure GDA00029528405000001321
The amount of reject water at the time period t,
Figure GDA00029528405000001322
for water flow from hydropower stations
Figure GDA00029528405000001323
Flows to downstream hydropower stations
Figure GDA00029528405000001324
The time of (a) is,
Figure GDA00029528405000001325
and
Figure GDA00029528405000001326
for the step hydropower station data acquired by the scheduling authority,
Figure GDA00029528405000001327
the reporting amount is prevented from being reduced for the regulating hydropower k in the nth section of the time period t, and the reporting amount is prevented from being reduced for the regulating hydropower k participating in the day-ahead market in the second day section;
the constraint condition group reflects the operation constraint of the adjustable hydropower station, wherein a formula (13) is a hydropower station water storage quantity constraint, a formula (14) is a reservoir dynamic balance equation constraint, a formula (15) is an upstream and downstream hydropower station coupling constraint, a formula (16) is an adjustable hydropower station output constraint, a formula (17) is an adjustable hydropower station output constraint if water is not abandoned, and a formula (18) is a reduction constraint for avoiding the reduction of the declared quantity in each section;
1-2-5) constraint conditions of the thermal power generating unit:
Figure GDA00029528405000001328
Figure GDA00029528405000001329
Figure GDA00029528405000001330
Figure GDA00029528405000001331
Figure GDA0002952840500000141
Figure GDA0002952840500000142
Figure GDA0002952840500000143
Figure GDA0002952840500000144
Figure GDA0002952840500000145
Figure GDA0002952840500000146
Figure GDA0002952840500000147
wherein the content of the first and second substances,
Figure GDA0002952840500000148
the minimum output of the thermal power generating unit k is the minimum economic output reported by the thermal power generating unit to the power dispatching mechanism if the thermal power generating unit is not the necessary starting group, the minimum economic output is the output of the dispatching mechanism for decomposing the physical contract electric quantity of the thermal power generating unit to each time period of the second day if the thermal power generating unit is the necessary starting group,
Figure GDA0002952840500000149
the unit capacity of the thermal power generating unit k is the unit capacity reported to the power dispatching mechanism by the thermal power generating unit,
Figure GDA00029528405000001410
is an integer variable of the starting state of the thermal power generating unit, if the thermal power generating unit k is in the starting state,
Figure GDA00029528405000001411
if the thermal power generating unit k is not in the on state,
Figure GDA00029528405000001412
the maximum climbing rate of the thermal power generating unit k and the maximum climbing rate constraint of the thermal power generating unit reported to the dispatching mechanism,
Figure GDA00029528405000001413
starting an integer variable for the thermal power generating unit, if the thermal power generating unit k is started in a time period t,
Figure GDA00029528405000001414
if the thermal power generating unit k is not started in the time period t,
Figure GDA00029528405000001415
is a thermal power generating unit shutdown integer variable, if the thermal power generating unit k is shutdown in a time period t,
Figure GDA00029528405000001416
if the thermal power generating unit k is not shut down in the time period t,
Figure GDA00029528405000001417
Figure GDA00029528405000001418
for the minimum on-state duration of the thermal power generating unit k,
Figure GDA00029528405000001419
is the minimum shutdown state duration of the thermal power generating unit k,
Figure GDA00029528405000001420
the maximum starting times of the thermal power generating unit k in one day,
Figure GDA00029528405000001421
in order to increase the starting cost of the thermal power generating unit k,
Figure GDA00029528405000001422
and
Figure GDA00029528405000001423
the minimum startup state duration, the minimum shutdown state duration, the maximum startup times in a day and the startup cost data of the thermal power unit are reported to a scheduling mechanism by the thermal power unit,
Figure GDA00029528405000001424
determining information whether the thermal power unit is a necessary unit for a scheduling mechanism for necessary constraint parameters of the thermal power unit k, if the thermal power unit k is the necessary unit,
Figure GDA00029528405000001425
if the thermal power generating unit k does not need to be started,
Figure GDA00029528405000001426
the method comprises the steps that an output declaration value of the nth energy quotation section of the thermal power generating unit k is obtained, and a second day subsection output is reported to a power trading mechanism by the thermal power generating unit participating in the market in the day ahead;
the constraint set reflects the operation constraint of the thermal power unit, formulas (19) and (20) are the output range constraint of the thermal power unit, the deep peak shaving of the thermal power unit is embodied in the formula (19), the climbing rate constraint of the thermal power unit is embodied in the formula (21), the formula (22) is the minimum starting state duration constraint of the thermal power unit, the formula (23) is the minimum stopping state duration constraint of the thermal power unit, the formula (24) is the maximum starting time constraint of the thermal power unit, the formulas (25) and (26) are relationship constraints among integer variables, the formula (27) is the starting cost constraint of the thermal power unit, the formula (28) is the necessary starting group constraint, and the formula (29) is the winning output constraint of each energy quotation section of the thermal power unit;
1-3) constructing a mathematical model of the daily market clearing of the uniformly reduced electric power without considering constraint condition elasticization, wherein the expressions are respectively as follows:
1-3-1) objective function of the electric power day-ahead market clearing mathematical model with uniform curtailment elasticated without consideration of constraints:
Figure GDA0002952840500000151
wherein the content of the first and second substances,
Figure GDA0002952840500000152
and
Figure GDA0002952840500000153
are respectively
Figure GDA0002952840500000154
And
Figure GDA0002952840500000155
the weight coefficient of (a) is,
Figure GDA0002952840500000156
and as n is increased,
Figure GDA0002952840500000157
and
Figure GDA0002952840500000158
the size of the composite material is also increased continuously,
Figure GDA0002952840500000159
for the uniform reduction of the uncontrollable renewable energy k in the nth reduction segment of the time period t,
Figure GDA00029528405000001510
the uniform reduction amount of the adjustable hydropower k in the nth section reduction section of the time period t is reduced;
the objective function shows that when the power system can not consume the renewable energy, the renewable energy is evenly reduced;
1-3-2) variable constraints of the electric power day-ahead market clearing mathematical model without consideration of constraint condition elasticization:
Figure GDA00029528405000001511
Figure GDA00029528405000001512
n is large enough to make uncontrollable renewable energy k output predicted in time period t
Figure GDA00029528405000001513
And adjusting the power of the hydropower k in the time period t
Figure GDA00029528405000001514
Can be divided small enough so that renewable energy can be evenly curtailed;
1-3-3) the constraint conditions of the power system of the electric power daily market clearing mathematical model which do not consider the constraint conditions for the elasticization and the uniform reduction are the same as all the constraint conditions of the step 1-2-1);
1-3-4) the constraint conditions of uncontrollable renewable energy sources such as radial flow type hydropower, wind power, photovoltaic and the like of a mathematical model of the market clearing at the day before the electric power with the constraint conditions of the elasticization and the uniform reduction are not considered, and the constraint conditions are the same as all the constraint conditions of the step 1-2-2);
1-3-5) adjusting type hydropower linearization constraint conditions of a power day-ahead market clearing mathematical model, which do not consider constraint conditions for elasticization and uniform reduction, are the same as all the constraint conditions of the step 1-2-3);
1-3-6) considering the adjustable hydropower constraint condition of the power day-ahead market clearing mathematical model with the constraint condition elasticization and uniform reduction, and adding the following constraint conditions on the basis of all the constraint conditions of the step 1-2-4):
Figure GDA0002952840500000161
wherein the content of the first and second substances,
Figure GDA0002952840500000162
market-by-day reconciliation of step 1-2) constraint elasticizationThe water discharge amount of the water-saving type hydropower k in the time period t is calculated,
Figure GDA0002952840500000163
calculating the water abandon amount of the adjustable hydropower k in the time period t in the constraint condition elasticized day-ahead market in the step 1-2);
1-3-7) eliminating constraint conditions and elasticizing the constraint conditions, and uniformly reducing the constraint conditions to obtain the constraint conditions of the thermal power generating unit of the mathematical model in the market before the day:
Figure GDA0002952840500000164
wherein the content of the first and second substances,
Figure GDA0002952840500000165
and
Figure GDA0002952840500000166
respectively calculating results of the output and the depth peak regulation amount of the fire generator set k in the time period t in the constraint condition elasticized day-ahead market in the step 1-2);
the constraint condition ensures that the thermal power generating unit does not carry out deep peak regulation intentionally any more, and the output is fixed on the clear result or the minimum output of the market in the day before the constraint condition is elasticized;
2) solving the constraint condition elasticized electric power day-ahead market clearing mathematical model in the step 1-2) by adopting a mixed integer programming solving method to obtain constraint condition elasticized electric power day-ahead market clearing results, wherein the clearing results comprise the medium-output power of each thermal power generating unit in each quotation section in each time period
Figure GDA0002952840500000167
Output of each uncontrollable renewable energy source in each time period
Figure GDA0002952840500000168
The output of each adjustable water and electricity in each time interval
Figure GDA0002952840500000169
Reduction amount of each uncontrollable renewable energy in each quotation section of each time period
Figure GDA00029528405000001610
Reduction of uncontrollable renewable energy sources in each time period
Figure GDA00029528405000001611
Reduction of each adjusting hydropower in each quotation section of each time period
Figure GDA00029528405000001612
Cutting amount of each adjustable water and electricity in each time period
Figure GDA00029528405000001613
Mid-bid depth peak regulation amount of each thermal power generating unit in each quotation section of each time period
Figure GDA00029528405000001614
3) Solving the uniform reduction daily market clearing mathematical model which does not consider the elasticization of the constraint conditions and is obtained by adopting a mixed integer programming solving method, wherein the uniform reduction daily market clearing mathematical model which does not consider the elasticization of the constraint conditions in the step 1-3) is obtained, and the clearing result comprises the uniform reduction amount of each uncontrollable renewable energy source in each reduction section in each time period
Figure GDA00029528405000001615
Even reduction of uncontrollable renewable energy sources in each time period
Figure GDA00029528405000001616
The even reduction amount of each reduction section of each adjusting type water and electricity in each time period
Figure GDA00029528405000001617
The even reduction of each adjustable water and electricity in each time interval
Figure GDA00029528405000001618
4) According to the aboveQuoting each thermal power generating unit in step 1) in each period
Figure GDA00029528405000001619
And the peak regulation amount of the thermal power generating units in the step 2) is the mid-bid depth of each quotation section in each period
Figure GDA0002952840500000171
Taking the highest price of the deep peak shaving in the deep peak shaving price quotation section of the winning bid of all the thermal power units in the time period as the deep peak shaving clear price;
5) calculating renewable energy avoids cutting clearing prices, including:
5-1) the reduction amount of each uncontrollable renewable energy source in each time period according to the step 2) above
Figure GDA0002952840500000172
Cutting amount of each adjustable water and electricity in each time period
Figure GDA0002952840500000173
Even reduction of uncontrollable renewable energy sources in each time period
Figure GDA0002952840500000174
And the even reduction amount of each adjustable hydropower in each time period in the step 3) is reduced
Figure GDA0002952840500000175
Subtracting the reduction amount of each renewable energy source in each time period in the step 2) from the uniform reduction amount of each renewable energy source in each time period in the step 3) to obtain the amount of each renewable energy source which is not discarded in each time period;
5-2) avoiding cutting off quotations at each stage of each time interval according to each uncontrollable renewable energy in the step 1) above
Figure GDA0002952840500000176
Avoiding reduction of quotations at each section of each time period by each adjusting type hydropower in step 1)
Figure GDA0002952840500000177
The reduction amount of each uncontrollable renewable energy source obtained in the step 2) in each time period
Figure GDA0002952840500000178
The reduction amount of each adjusting type hydropower obtained in the step 2) in each time period
Figure GDA0002952840500000179
Judging the amount of each renewable energy source to be abandoned in the period of time according to the amount of each renewable energy source to be abandoned obtained in the step 5-1), if the amount of each renewable energy source to be abandoned in the period of time is positive, taking the reduction-avoiding quoted price corresponding to the reduction amount in the period of time obtained in the step 2) of the renewable energy source as the reduction-avoiding quoted price of the renewable energy source acting in the period of time, and if the amount of each renewable energy source to be abandoned in the period of time is not positive, avoiding reduction-avoiding quoted price of the renewable energy source acting in the period of time does not exist;
5-3) according to the avoidance reduction quotation of each renewable energy source acting in each time interval obtained in the step 5-2), taking the lowest quotation of the avoidance reduction quotation of each renewable energy source acting in each time interval as the avoidance reduction clearing price of the power system in the time interval, and realizing the clearing of the power market before the day with the elastic constraint condition.
It is worth mentioning that the objective function in the implementation steps provided by the invention can be flexibly selected and customized according to needs, constraint conditions can be added and deleted according to actual needs, and the expandability is strong; therefore, the above implementation steps are only used for illustrating and not limiting the technical solution of the present invention; any modification or partial replacement without departing from the spirit and scope of the present invention should be covered in the claims of the present invention.

Claims (1)

1. A constraint condition-elasticated electric power day-ahead market clearing method is characterized by comprising the following steps:
1) constructing a constraint condition elasticized mathematical model of the clear electric power market at the day before, wherein the involved energy sources are in two forms, namely thermal power and renewable energy sources, the renewable energy sources comprise an uncontrollable renewable energy source and an adjustable hydropower form, the uncontrollable renewable energy sources are respectively radial flow hydropower, wind power and photovoltaic, and the specific process is as follows;
1-1) constructing an objective function of a constraint condition-elasticated electric power day-ahead market clearing mathematical model as follows:
Figure FDA0002952840490000011
the objective function represents the minimum total operating cost of the power system at the maximum renewable energy consumption, where the upper scale th represents thermal power, the upper scale r represents uncontrollable renewable energy, the upper scale h represents regulated hydropower,
Figure FDA0002952840490000012
the price of the nth energy quotation section of the thermal power generating unit k is the price of the unit kilowatt-hour electric quantity corresponding to the nth energy quotation section of the second day reported to the electric power trading mechanism by the thermal power generating unit participating in the market in the day ahead,
Figure FDA0002952840490000013
the nth energy quotation section of the thermal power generating unit is subjected to the medium-winning amount in the time period t,
Figure FDA0002952840490000014
the starting cost of the thermal power generating unit k in the time period t is shown, M is a weight coefficient, and the value range of M is
Figure FDA0002952840490000015
Figure FDA0002952840490000016
Avoiding reduction of quotes for uncontrollable renewable energy k in the nth period of time t, reporting to the electric power exchange for uncontrollable renewable energy participating in the market in the dayThe daily subsection avoids reducing price corresponding to the reduction of the declaration quantity,
Figure FDA0002952840490000017
to avoid the reduction of the reduction declaration quantity in the nth period of the time t for the uncontrollable renewable energy source k,
Figure FDA0002952840490000018
for avoiding reduction of quotations of the regulated hydropower k in the nth section of the time period t, for avoiding reduction of prices corresponding to reduction of reporting quantity in the second day section by section for reporting the regulated hydropower participating in the market in the day ahead to the electric power transaction institution,
Figure FDA0002952840490000019
to avoid the reduction of the declaration quantity in the nth period of the time t,
Figure FDA00029528404900000110
giving a thermal power generating unit k a deep peak regulation price in the nth section of the time interval t, giving out the price of electric quantity in kilowatt-hours corresponding to the second-day segmented deep peak regulation capacity reported to the electric power trading mechanism by the thermal power generating unit participating in the market in the day ahead,
Figure FDA00029528404900000111
the method comprises the steps of obtaining the quantity of the thermal power generating unit k for the nth section depth peak regulation in a time period t;
1-2) constructing constraint conditions of a constraint condition-elasticized electric power day-ahead market clearing mathematical model, wherein expressions are respectively as follows:
1-2-1) power system constraints:
Figure FDA0002952840490000021
Figure FDA0002952840490000022
wherein the content of the first and second substances,
Figure FDA0002952840490000023
indicating that for all of the (-) s that are true,
Figure FDA0002952840490000024
the output of the uncontrollable renewable energy k during the time period t,
Figure FDA0002952840490000025
for regulating the output of the hydropower k over a period of time t, DtFor the total load of the power system in the time period t, the predicted value P of the total load demand of the power system in each time period of the next day is obtained by the power dispatching mechanism by adopting a prediction methodLmaxThe maximum active transmission power of the line L and the maximum active transmission power of each line on the second day given by the dispatching mechanism, G is a generator output power transfer distribution factor matrix of the power system, generator output power transfer distribution factor data of each power transmission line by each node given by the dispatching mechanism, and Gr-k,LThe distribution factor of the output power transfer of the generator from the node where the uncontrollable renewable energy k is located to the line L, r in the subscript represents the uncontrollable renewable energy, Gh-k,LFor the distribution factor of the generator output power transfer from the node where the regulated hydropower k is located to the line L, h in the subscript denotes the regulated hydropower, Gth-k,LFor a generator output power transfer distribution factor from a node where a thermal power generating unit k is located to a line L, th in subscript represents thermal power, Gb,LFor generator output power transfer distribution factor, D, from node b to line Lb,tFor the load of the node B in the time period t, obtaining a predicted value of the power load demand of each node in each time period of the second day by a power dispatching mechanism by adopting a prediction method, wherein B is the total number of nodes and the total number of power network nodes given by the dispatching mechanism;
the formula (2) is a load balance constraint of the power system, and represents that the total power generation amount of the power system in each period is equal to the total load of the power system;
formula (3) is network constraint, and indicates that the active transmission power of each line must not exceed the active transmission capacity;
1-2-2) uncontrollable renewable energy constraint conditions:
Figure FDA0002952840490000026
Figure FDA0002952840490000027
wherein the content of the first and second substances,
Figure FDA0002952840490000028
predicted output for the uncontrollable renewable energy k in the time period t, predicted output value for each time period of the second day predicted for each uncontrollable renewable energy,
Figure FDA0002952840490000029
the reduction of the reporting amount is avoided for the uncontrollable renewable energy k in the nth period of the time period t, and the reduction of the reporting amount is avoided for the uncontrollable renewable energy participating in the market in the second day and reported to the electric power exchange mechanism in sections; the constraint condition group indicates that the sum of the winning output and the reduction of the uncontrollable renewable energy k in the time period t is equal to the predicted output of the uncontrollable renewable energy k in the time period t, and the reduction of the avoidance reduction declaration quantity of each time period is not higher than the avoidance reduction declaration quantity of the renewable energy in the time period;
1-2-3) regulating type hydropower linearization constraint condition
Figure FDA0002952840490000031
Figure FDA0002952840490000032
Figure FDA0002952840490000033
Figure FDA0002952840490000034
Figure FDA0002952840490000035
Figure FDA0002952840490000036
Figure FDA0002952840490000037
Figure FDA0002952840490000038
Wherein the content of the first and second substances,
Figure FDA0002952840490000039
for the lower limit of the applicable water storage capacity range of the regulated hydropower k c output curve,
Figure FDA00029528404900000310
for the upper limit of the applicable water storage capacity range of the c output curve of the regulating hydropower k,
Figure FDA00029528404900000311
to regulate the amount of water stored by the hydro-electric k over a period of time t,
Figure FDA00029528404900000312
to represent integer variables applicable to the linear section, if the regulated hydropower k is in the w-th section water discharge amount range of the c-th output curve in the time period t, then
Figure FDA00029528404900000313
If the adjusting type hydropower k is not in the w-th section water discharge amount range of the c-th output curve in the time period t, the water discharge amount of the adjusting type hydropower k is adjusted
Figure FDA00029528404900000314
Figure FDA00029528404900000315
For regulating the lower limit of the water discharge quantity range of the w section of the c output curve of the hydropower k,
Figure FDA00029528404900000316
the water discharge upper limit of the w-th section of the water discharge range of the c-th output curve of the adjustable hydropower k,
Figure FDA00029528404900000317
to adjust the discharge of the hydropower k at the time period t,
Figure FDA00029528404900000318
in order to adjust the position proportion of the water discharge quantity of the hydropower k in the w-th section water discharge quantity range of the c-th output curve,
Figure FDA00029528404900000319
in order to adjust the output of the water discharge quantity of the hydropower k at the w-th section water discharge quantity lower limit of the c-th output curve,
Figure FDA00029528404900000320
in order to adjust the output of the water discharge quantity of the hydropower k at the w-th section water discharge quantity upper limit of the c-th output curve,
Figure FDA00029528404900000321
Figure FDA00029528404900000322
and
Figure FDA00029528404900000323
the output of the hydroelectric generating set corresponding to different water discharge amounts in different water storage amount sections reported to the dispatching mechanism is a piecewise linear function;
the restriction condition set is obtained by carrying out linear treatment on the nonlinear and non-convex output curve of the hydropower so as to approximately obtain the hydropower output
Figure FDA00029528404900000324
Water discharge amount with water and electricity
Figure FDA00029528404900000325
And the water storage capacity
Figure FDA00029528404900000326
The relationship between;
1-2-4) regulating hydropower constraint conditions:
Figure FDA0002952840490000041
Figure FDA0002952840490000042
Figure FDA0002952840490000043
Figure FDA0002952840490000044
Figure FDA0002952840490000045
Figure FDA0002952840490000046
wherein the content of the first and second substances,
Figure FDA0002952840490000047
and
Figure FDA0002952840490000048
respectively the minimum and maximum water storage capacity of the adjustable hydroelectric power k, the maximum/minimum water storage capacity constraint reported by the hydroelectric generating set,
Figure FDA0002952840490000049
for the predicted inflow of the regulated hydropower k at time period t, for the second water-day prediction data reported by the hydroelectric generating set to the scheduling authority,
Figure FDA00029528404900000410
to adjust the amount of waste water of the hydro-power k during the time period t,
Figure FDA00029528404900000411
the ith hydroelectric station of basin p,
Figure FDA00029528404900000412
for hydropower stations
Figure FDA00029528404900000413
In a period of time
Figure FDA00029528404900000414
The amount of the supplied water is controlled by the amount of the supplied water,
Figure FDA00029528404900000415
hydropower station
Figure FDA00029528404900000416
In a period of time
Figure FDA00029528404900000417
The amount of natural water coming in is predicted,
Figure FDA00029528404900000418
for hydropower stations
Figure FDA00029528404900000419
At the amount of water discharge for the period t,
Figure FDA00029528404900000420
for hydropower stations
Figure FDA00029528404900000421
The amount of reject water at the time period t,
Figure FDA00029528404900000422
for water flow from hydropower stations
Figure FDA00029528404900000423
Flows to downstream hydropower stations
Figure FDA00029528404900000424
The time of (a) is,
Figure FDA00029528404900000425
and
Figure FDA00029528404900000426
for the step hydropower station data acquired by the scheduling authority,
Figure FDA00029528404900000427
the reporting amount is prevented from being reduced for the regulating hydropower k in the nth section of the time period t, and the reporting amount is prevented from being reduced for the regulating hydropower k participating in the day-ahead market in the second day section;
the constraint condition group reflects the operation constraint of the adjustable hydropower station, wherein a formula (13) is a hydropower station water storage quantity constraint, a formula (14) is a reservoir dynamic balance equation constraint, a formula (15) is an upstream and downstream hydropower station coupling constraint, a formula (16) is an adjustable hydropower station output constraint, a formula (17) is an adjustable hydropower station output constraint if water is not abandoned, and a formula (18) is a reduction constraint for avoiding the reduction of the declared quantity in each section;
1-2-5) constraint conditions of the thermal power generating unit:
Figure FDA00029528404900000428
Figure FDA00029528404900000429
Figure FDA00029528404900000430
Figure FDA00029528404900000431
Figure FDA0002952840490000051
Figure FDA0002952840490000052
Figure FDA0002952840490000053
Figure FDA0002952840490000054
Figure FDA0002952840490000055
Figure FDA0002952840490000056
Figure FDA0002952840490000057
wherein the content of the first and second substances,
Figure FDA0002952840490000058
the minimum output of the thermal power generating unit k is the minimum economic output reported by the thermal power generating unit to the power dispatching mechanism if the thermal power generating unit is not the necessary starting group, the minimum economic output is the output of the dispatching mechanism for decomposing the physical contract electric quantity of the thermal power generating unit to each time period of the second day if the thermal power generating unit is the necessary starting group,
Figure FDA0002952840490000059
the unit capacity of the thermal power generating unit k is the unit capacity reported to the power dispatching mechanism by the thermal power generating unit,
Figure FDA00029528404900000510
is an integer variable of the starting state of the thermal power generating unit, if the thermal power generating unit k is in the starting state,
Figure FDA00029528404900000511
if the thermal power generating unit k is not in the on state,
Figure FDA00029528404900000512
Figure FDA00029528404900000513
the maximum climbing rate of the thermal power generating unit k and the maximum climbing rate constraint of the thermal power generating unit reported to the dispatching mechanism,
Figure FDA00029528404900000514
starting an integer variable for the thermal power generating unit, if the thermal power generating unit k is started in a time period t,
Figure FDA00029528404900000515
if the thermal power generating unit k is not started in the time period t,
Figure FDA00029528404900000516
Figure FDA00029528404900000517
is a thermal power generating unit shutdown integer variable, if the thermal power generating unit k is shutdown in a time period t,
Figure FDA00029528404900000518
if the thermal power generating unit k is not shut down in the time period t,
Figure FDA00029528404900000519
Figure FDA00029528404900000520
for the minimum on-state duration of the thermal power generating unit k,
Figure FDA00029528404900000521
is the minimum shutdown state duration of the thermal power generating unit k,
Figure FDA00029528404900000522
the maximum starting times of the thermal power generating unit k in one day,
Figure FDA00029528404900000523
in order to increase the starting cost of the thermal power generating unit k,
Figure FDA00029528404900000524
and
Figure FDA00029528404900000525
the minimum startup state duration, the minimum shutdown state duration, the maximum startup times in a day and the startup cost data of the thermal power unit are reported to a scheduling mechanism by the thermal power unit,
Figure FDA00029528404900000526
determining information whether the thermal power unit is a necessary unit for a scheduling mechanism for necessary constraint parameters of the thermal power unit k, if the thermal power unit k is the necessary unit,
Figure FDA00029528404900000527
if the thermal power generating unit k does not need to be started,
Figure FDA00029528404900000528
Figure FDA00029528404900000529
the method comprises the steps that an output declaration value of the nth energy quotation section of the thermal power generating unit k is obtained, and a second day subsection output is reported to a power trading mechanism by the thermal power generating unit participating in the market in the day ahead;
the constraint condition group reflects the operation constraint of the thermal power unit, formulas (19) and (20) are the output range constraint of the thermal power unit, the deep peak shaving of the thermal power unit is embodied in the formula (19), the climbing rate constraint of the thermal power unit is embodied in the formula (21), the formula (22) is the minimum starting state duration constraint of the thermal power unit, the formula (23) is the minimum stopping state duration constraint of the thermal power unit, the formula (24) is the maximum starting time constraint of the thermal power unit, the formulas (25) and (26) are the relation constraint between integer variables, the formula (27) is the starting cost constraint of the thermal power unit, the formula (28) is the necessary starting group constraint, and the formula (29) is the medium output constraint of each energy quotation section of the thermal power unit;
1-3) constructing a mathematical model of the daily market clearing of the uniformly reduced electric power without considering constraint condition elasticization, wherein the expressions are respectively as follows:
1-3-1) objective function of the electric power day-ahead market clearing mathematical model with uniform curtailment elasticated without consideration of constraints:
Figure FDA0002952840490000061
wherein the content of the first and second substances,
Figure FDA0002952840490000062
and
Figure FDA0002952840490000063
are respectively
Figure FDA0002952840490000064
And
Figure FDA0002952840490000065
the weight coefficient of (a) is,
Figure FDA0002952840490000066
and as n is increased,
Figure FDA0002952840490000067
and
Figure FDA0002952840490000068
the size of the composite material is also increased continuously,
Figure FDA0002952840490000069
for the uniform reduction of the uncontrollable renewable energy k in the nth reduction segment of the time period t,
Figure FDA00029528404900000610
the uniform reduction amount of the adjustable hydropower k in the nth section reduction section of the time period t is reduced;
the objective function shows that when the power system can not consume the renewable energy, the renewable energy is evenly reduced;
1-3-2) variable constraints of the electric power day-ahead market clearing mathematical model without consideration of constraint condition elasticization:
Figure FDA00029528404900000611
Figure FDA00029528404900000612
n is large enough to make uncontrollable renewable energy k output predicted in time period t
Figure FDA00029528404900000613
And adjusting the power of the hydropower k in the time period t
Figure FDA00029528404900000614
Can be divided small enough so that renewable energy can be evenly curtailed;
1-3-3) the constraint conditions of the power system of the electric power daily market clearing mathematical model which do not consider the constraint conditions for the elasticization and the uniform reduction are the same as all the constraint conditions of the step 1-2-1);
1-3-4) the constraint conditions of uncontrollable renewable energy sources such as radial flow type hydropower, wind power, photovoltaic and the like of a mathematical model of the market clearing at the day before the electric power with the constraint conditions of the elasticization and the uniform reduction are not considered, and the constraint conditions are the same as all the constraint conditions of the step 1-2-2);
1-3-5) adjusting type hydropower linearization constraint conditions of a power day-ahead market clearing mathematical model, which do not consider constraint conditions for elasticization and uniform reduction, are the same as all the constraint conditions of the step 1-2-3);
1-3-6) considering the adjustable hydropower constraint condition of the power day-ahead market clearing mathematical model with the constraint condition elasticization and uniform reduction, and adding the following constraint conditions on the basis of all the constraint conditions of the step 1-2-4):
Figure FDA0002952840490000071
wherein the content of the first and second substances,
Figure FDA0002952840490000072
calculating the water discharge quantity of the adjustable hydropower k in the time period t in the constraint condition elasticized day-ahead market in the step 1-2),
Figure FDA0002952840490000073
calculating the water abandon amount of the adjustable hydropower k in the time period t in the constraint condition elasticized day-ahead market in the step 1-2);
1-3-7) eliminating constraint conditions and elasticizing the constraint conditions, and uniformly reducing the constraint conditions to obtain the constraint conditions of the thermal power generating unit of the mathematical model in the market before the day:
Figure FDA0002952840490000074
wherein the content of the first and second substances,
Figure FDA0002952840490000075
and
Figure FDA0002952840490000076
respectively calculating results of the output and the depth peak regulation amount of the fire generator set k in the time period t in the constraint condition elasticized day-ahead market in the step 1-2);
the constraint condition ensures that the thermal power generating unit does not carry out deep peak regulation intentionally any more, and the output is fixed on the clear result or the minimum output of the market in the day before the constraint condition is elasticized;
2) solving the constraint condition elasticized electric power day-ahead market clearing mathematical model in the step 1-2) by adopting a mixed integer programming solving method to obtain constraint condition elasticized electric power day-ahead market clearing results, wherein the clearing results comprise the medium-output power of each thermal power generating unit in each quotation section in each time period
Figure FDA0002952840490000077
Output of each uncontrollable renewable energy source in each time period
Figure FDA0002952840490000078
The output of each adjustable water and electricity in each time interval
Figure FDA0002952840490000079
Reduction amount of each uncontrollable renewable energy in each quotation section of each time period
Figure FDA00029528404900000710
Reduction of uncontrollable renewable energy sources in each time period
Figure FDA00029528404900000711
Reduction of each adjusting hydropower in each quotation section of each time period
Figure FDA00029528404900000712
Cutting amount of each adjustable water and electricity in each time period
Figure FDA00029528404900000713
Mid-bid depth peak regulation amount of each thermal power generating unit in each quotation section of each time period
Figure FDA00029528404900000714
3) Solving the uniform reduction daily market clearing mathematical model which does not consider the elasticization of the constraint conditions and is obtained by adopting a mixed integer programming solving method, wherein the uniform reduction daily market clearing mathematical model which does not consider the elasticization of the constraint conditions in the step 1-3) is obtained, and the clearing result comprises the uniform reduction amount of each uncontrollable renewable energy source in each reduction section in each time period
Figure FDA00029528404900000715
Even reduction of uncontrollable renewable energy sources in each time period
Figure FDA00029528404900000716
The even reduction amount of each reduction section of each adjusting type water and electricity in each time period
Figure FDA00029528404900000717
The even reduction of each adjustable water and electricity in each time interval
Figure FDA00029528404900000718
4) Quoting each thermal power generating unit in each time interval according to the step 1) in each period
Figure FDA00029528404900000719
And the peak regulation amount of the thermal power generating units in the step 2) is the mid-bid depth of each quotation section in each period
Figure FDA0002952840490000081
Taking the highest price of the deep peak shaving in the deep peak shaving price quotation section of the winning bid of all the thermal power units in the time period as the deep peak shaving clear price;
5) calculating renewable energy avoids cutting clearing prices, including:
5-1) the reduction amount of each uncontrollable renewable energy source in each time period according to the step 2) above
Figure FDA0002952840490000082
Cutting amount of each adjustable water and electricity in each time period
Figure FDA0002952840490000083
Even reduction of uncontrollable renewable energy sources in each time period
Figure FDA0002952840490000084
And the even reduction amount of each adjustable hydropower in each time period in the step 3) is reduced
Figure FDA0002952840490000085
Each renewable energy sourceSubtracting the reduction amount of each renewable energy source in each time period in the step 2) from the uniform reduction amount of each time period in the step 3) to obtain the amount of each renewable energy source which is not abandoned in each time period;
5-2) avoiding cutting off quotations at each stage of each time interval according to each uncontrollable renewable energy in the step 1) above
Figure FDA0002952840490000086
Avoiding reduction of quotations at each section of each time period by each adjusting type hydropower in step 1)
Figure FDA0002952840490000087
The reduction amount of each uncontrollable renewable energy source obtained in the step 2) in each time period
Figure FDA0002952840490000088
The reduction amount of each adjusting type hydropower obtained in the step 2) in each time period
Figure FDA0002952840490000089
Judging the amount of each renewable energy source to be abandoned in the period of time according to the amount of each renewable energy source to be abandoned obtained in the step 5-1), if the amount of each renewable energy source to be abandoned in the period of time is positive, taking the reduction-avoiding quoted price corresponding to the reduction amount in the period of time obtained in the step 2) of the renewable energy source as the reduction-avoiding quoted price of the renewable energy source acting in the period of time, and if the amount of each renewable energy source to be abandoned in the period of time is not positive, avoiding reduction-avoiding quoted price of the renewable energy source acting in the period of time does not exist;
5-3) according to the avoidance reduction quotation of each renewable energy source acting in each time interval obtained in the step 5-2), taking the lowest quotation of the avoidance reduction quotation of each renewable energy source acting in each time interval as the avoidance reduction clearing price of the power system in the time interval, and realizing the clearing of the power market before the day with the elastic constraint condition.
CN201811025518.5A 2018-09-04 2018-09-04 Constraint condition elasticized electric power day-ahead market clearing method Active CN109301876B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811025518.5A CN109301876B (en) 2018-09-04 2018-09-04 Constraint condition elasticized electric power day-ahead market clearing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811025518.5A CN109301876B (en) 2018-09-04 2018-09-04 Constraint condition elasticized electric power day-ahead market clearing method

Publications (2)

Publication Number Publication Date
CN109301876A CN109301876A (en) 2019-02-01
CN109301876B true CN109301876B (en) 2021-08-13

Family

ID=65166209

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811025518.5A Active CN109301876B (en) 2018-09-04 2018-09-04 Constraint condition elasticized electric power day-ahead market clearing method

Country Status (1)

Country Link
CN (1) CN109301876B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110165711B (en) * 2019-02-28 2020-11-24 重庆大学 Wind, water and fire combined scheduling optimization method considering elastic space of power generation operation
CN111476410B (en) * 2020-03-31 2022-10-14 云南电网有限责任公司 Electric quantity settlement method and device for guaranteeing clean energy consumption and unit fair bidding
CN111861162B (en) * 2020-07-06 2023-12-19 大连理工大学 Method for clearing daily current goods of high-proportion hydroelectric system under step electric quantity linkage control
CN112803479B (en) * 2021-01-06 2023-07-21 国网辽宁省电力有限公司 Thermal power generating unit modeling method and model considering depth peak shaving
CN113822707B (en) * 2021-09-10 2024-04-30 国网冀北电力有限公司电力科学研究院 Output decision method and device for electric power market, computer equipment and storage medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106022644A (en) * 2016-06-06 2016-10-12 中国南方电网有限责任公司电网技术研究中心 Construction method of cross-region power generation right transaction-carbon transaction coupling model
US9704206B2 (en) * 2012-10-24 2017-07-11 Causam Energy, Inc. System, method, and apparatus for settlement for participation in an electric power grid
CN107370164A (en) * 2017-07-12 2017-11-21 国电南瑞科技股份有限公司 Consider that the frequency modulation assisted hatching of performance criteria of the response goes out clear decision-making technique
CN107862586A (en) * 2017-12-04 2018-03-30 清华大学 A kind of energy hinge competitive tender method and system towards the transaction of multipotency stream

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7349883B2 (en) * 2003-03-10 2008-03-25 Siemens Power Transmission & Distribution, Inc. Optimized transmission and load security constrained unit commitment dispatch using linear programming for electricity markets

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9704206B2 (en) * 2012-10-24 2017-07-11 Causam Energy, Inc. System, method, and apparatus for settlement for participation in an electric power grid
CN106022644A (en) * 2016-06-06 2016-10-12 中国南方电网有限责任公司电网技术研究中心 Construction method of cross-region power generation right transaction-carbon transaction coupling model
CN107370164A (en) * 2017-07-12 2017-11-21 国电南瑞科技股份有限公司 Consider that the frequency modulation assisted hatching of performance criteria of the response goes out clear decision-making technique
CN107862586A (en) * 2017-12-04 2018-03-30 清华大学 A kind of energy hinge competitive tender method and system towards the transaction of multipotency stream

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"A Convex Model of Risk-Based Unit Commitment for Day-Ahead Market Clearing Considering Wind Power Uncertainty";Ning Zhang,;《IEEE Transactions on Power Systems 》;20141008;1582-1592 *
"约束条件弹性化的月度电力市场机制设计";舒畅;《中国电机工程学报》;20160205;587-594 *

Also Published As

Publication number Publication date
CN109301876A (en) 2019-02-01

Similar Documents

Publication Publication Date Title
CN109301876B (en) Constraint condition elasticized electric power day-ahead market clearing method
CN106300336B (en) It is a kind of meter and load side and source side virtual plant Multiobjective Optimal Operation method
CN110782281B (en) Day-ahead market clearing method for multi-owner cascade power station basin electric quantity transfer
JP2019525705A (en) Long-term federation peaking scheduling method for inter-provincial continuous hydropower plants
CN110336334B (en) Day-ahead peak regulation resource optimal scheduling method based on peak regulation capacity quotation of thermal power generating unit
CN109636674B (en) Large-scale hydropower station group monthly transaction electric quantity decomposition and checking method
CN112039056A (en) Two-stage optimal scheduling method for new energy
CN107451670B (en) Day-ahead hierarchical coordination scheduling method for active power distribution network
Wang et al. Short-term optimal scheduling of cascade hydropower plants shaving peak load for multiple power grids
CN109687534B (en) Power system generator set active power control method based on step water quantity matching
CN113538066B (en) Control method, system, equipment and medium for virtual power plant to participate in power market
CN110649598B (en) Method and system for regulating node electricity price by virtual power plant in area
CN104299072A (en) Security constraint generation schedule planning method based on water and fire coordination
Zhang et al. An ultra-short-term scheduling model for cascade hydropower regulated by multilevel dispatch centers suppressing wind power volatility
Su et al. Optimization model for the short-term operation of hydropower plants transmitting power to multiple power grids via HVDC transmission lines
Peng et al. Optimized economic operation strategy for distributed energy storage with multi-profit mode
CN108879657A (en) Balance of electric power and ener optimization method based on wind-powered electricity generation volume metering
CN117277444B (en) New energy base power capacity optimal configuration method and device
CN110929959B (en) Method and device for allocating local deep peak regulation expense of power grid
CN109726894A (en) Ensure the new energy active command calculation method of spot exchange and medium-term and long-term electricity
CN117374974A (en) Distribution network scheduling method, system, medium and equipment
CN113435673A (en) Deep peak regulation cost allocation calculation method and system based on load rate
CN112348274A (en) Peak shaving-considered power spot market clearing method, system, equipment and medium
Ma Short term hydropower dispatching optimization of cascaded hydropower stations based on two-stage optimization
Chakrabarti et al. Balancing wind intermittency using hydro reserves and demand response

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