CN107784375B - Coordination optimization method for bilateral power contract participating in day-ahead power and electric quantity balance - Google Patents

Coordination optimization method for bilateral power contract participating in day-ahead power and electric quantity balance Download PDF

Info

Publication number
CN107784375B
CN107784375B CN201610739177.2A CN201610739177A CN107784375B CN 107784375 B CN107784375 B CN 107784375B CN 201610739177 A CN201610739177 A CN 201610739177A CN 107784375 B CN107784375 B CN 107784375B
Authority
CN
China
Prior art keywords
contract
bilateral
power
power plant
unit
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
CN201610739177.2A
Other languages
Chinese (zh)
Other versions
CN107784375A (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.)
State Grid Corp of China SGCC
State Grid Zhejiang Electric Power Co Ltd
China Electric Power Research Institute Co Ltd CEPRI
Original Assignee
State Grid Corp of China SGCC
State Grid Zhejiang Electric Power Co Ltd
China Electric Power Research Institute Co Ltd CEPRI
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, State Grid Zhejiang Electric Power Co Ltd, China Electric Power Research Institute Co Ltd CEPRI filed Critical State Grid Corp of China SGCC
Priority to CN201610739177.2A priority Critical patent/CN107784375B/en
Publication of CN107784375A publication Critical patent/CN107784375A/en
Application granted granted Critical
Publication of CN107784375B publication Critical patent/CN107784375B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06Electricity, gas or water supply
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS 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/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems 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)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

A coordination optimization method for participating in day-ahead power electric quantity balance in bilateral power contract comprises the following steps: acquiring the same electric quantity of three unions of each power plant in the future day, acquiring bilateral power contract data of the power plants or bilateral power contract data of the units, and determining contract reduction factors of each power plant; acquiring system load prediction data in the future 24 hours, basic parameters of each thermal power generating unit and coal consumption characteristics; establishing an optimization model according to the type of the bilateral power contract, determining a target function, setting constraint conditions and solving a mixed integer programming problem; and finally, determining a starting and stopping plan and a unit output plan of the unit in the future 24 hours and the reduction of the power plant. According to the technical scheme, the unit output is divided into a bilateral output part, a three-output part and other output parts, reasonable constraint is established, so that the bilateral actual output is not greater than the corresponding contractual output, the contractual adjustment amount in a target function is guaranteed to be non-negative, and the solving difficulty is reduced.

Description

Coordination optimization method for bilateral power contract participating in day-ahead power and electric quantity balance
Technical Field
The invention relates to the technical field of power markets, in particular to a coordination optimization method for participating in day-ahead power and electric quantity balance in a bilateral power contract.
Background
The electric quantity of the international typical electric power market is market electricity, and market electricity and planned electricity coexist in many regions and countries; unlike the PJM (Pennsylvania-New Jersey-Maryland) power market, the initial bilateral agreement of many regional and national power markets is a physical power contract with execution constraints, and day-ahead planning requires consideration of execution of bilateral transactions; different from the balance mechanism of developed countries, the market maturity of main bodies in the initial market of electric power markets in many regions and countries is not high, the price of increasing and decreasing output is not provided, and the adjustment of bilateral trading according to the price of increasing and decreasing output cannot be carried out.
The difference between the bilateral power contract and the bilateral power contract is that the power contract specifies a bilateral output curve, which will be the basis for bilateral settlement, which will affect the trend of the unit output curve, and the shape of the unit output curve will affect the completion rate of the bilateral contract. At present, research is mainly focused on wholesale competitive markets, bilateral transaction research, particularly quantitative analysis research, has few results, and is more focused on game models, quotation strategies, annual monthly decomposition and other problems of bilateral transaction. Therefore, the coordination optimization method for the day-ahead power and electric quantity balance participated by the bilateral power contract and the unit bilateral power contract of the power plant is provided, is suitable for day-ahead power optimization scheduling under various power contract types, is suitable for a market mode with coexistence of market power and planned power, and provides reference for realizing smooth transition of a power market.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention provides a coordination optimization method for participating in day-ahead power and electric quantity balance in a bilateral power contract.
The adopted solution for realizing the purpose is as follows:
a coordinated optimization method for enabling bilateral power contracts to participate in day-ahead power electricity quantity balance includes the following steps:
(1) acquiring the three-fair-contract same electric quantity of each power plant in the future one day, acquiring bilateral power contract data of the power plants or bilateral power contract data of the units according to the bilateral contract types, and determining contract reduction factors of each power plant;
(2) acquiring system load prediction data in the future 24 hours, basic parameters of each thermal power generating unit and coal consumption characteristics;
(3) if the bilateral power contract is the bilateral power contract of the power plant, turning to the step (4); if the bilateral power contract is the unit bilateral power contract, turning to the step (5);
(4) the method for establishing the day-ahead power and electric quantity balance coordination optimization model considering the bilateral power contract of the power plant comprises the following steps: determining an objective function and setting constraint conditions, and turning to the step (6);
(5) the method for establishing the day-ahead power and electric quantity balance coordination optimization model considering the unit bilateral power contract comprises the following steps: determining an objective function and setting constraint conditions, and turning to the step (6);
(6) solving a mixed integer programming problem;
(7) and determining a starting and stopping plan and a unit output plan of the unit in the future 24 hours and a power plant reduction amount.
Preferably, the step (1) includes:
a. the bilateral contract is the bilateral power contract of the power plant, and the bilateral contract quantity of the power plant in one day is counted according to a curve specified by the bilateral power contract of the power plant;
b. and the bilateral contract is the bilateral electric power contract of the units, the bilateral electric quantity of each unit in one day is counted according to a curve specified by the bilateral electric power contract of the units, and the bilateral contract quantities of all the units in one power plant are added to obtain the bilateral contract quantity of one power plant.
Preferably, in step (1), the reduction factor includes: the bilateral reduction factor is the proportion of the bilateral contract reduction amount of the power plant to the total bilateral reduction amount of the system; the three-fair reduction factor is the proportion of the three-fair reduction amount of the power plant to the total three-fair reduction amount of the system.
Preferably, in step (4), the optimization objective includes: minimizing the amount of contract adjustments and minimizing the cost of power generation; the minimized contract adjustment amount has higher priority;
the objective function is:
F=M1*(W1*ΔS+W2*ΔJ)+M2*C(P)
Figure BDA0001092896220000031
Figure BDA0001092896220000032
Figure BDA0001092896220000033
wherein F is the value of the objective function, M1Reducing penalty factors for contracts, M2For power generation cost weighting, W1Weight reduction for bilateral contracts, W2The weight is reduced for three fair contracts, Δ S is reduced for two-sided contracts, Δ J is reduced for three fair contracts, C (P) is the system power generation cost, Pg,0,tThe method is characterized in that the power is provided for bilateral contract regulation of the power plant in g time period T, H is the total number of the power plant, N is the total number of thermal power generating units, T is the total number of time periods, and P is the total number of the thermal power generating unitsi,S,tRepresenting the double-sided part of the unit output, Jg,0Is the same electric quantity of the power plant g in three commons, Pi,J,tRepresenting the three-revolution-force part of the unit force, fi,tFor the operating costs of thermal power units, SUi,t、SDi,tThe starting cost and the shutdown cost of the thermal power generating unit are respectively.
Preferably, in the step (4), the constraint condition includes:
and power balance constraint:
Figure BDA0001092896220000034
and (3) constraint of starting variables and stopping dynamic variables of the thermal power generating unit:
Ii,t-Ii,t-1=ui,t-vi,t,ui,t+vi,t≤1
minimum on-off time constraint:
Figure BDA0001092896220000041
vertical rotation standby restraint:
Figure BDA0001092896220000042
and (3) climbing restraint: pi,t-Pi,t-1≤Ri(1+Ii,t-1-Ii,t)+Pi,min(2-Ii,t-1-Ii,t)
Landslide restraint: pi,t-1-Pi,t≤Di(1-Ii,t-1+Ii,t)+Pi,min(2-Ii,t-1-Ii,t)
And (3) bilateral contract reduction constraint of the power plant g:
Figure BDA0001092896220000043
and (3) the relation constraint of the bilateral contract reduction amount of the power plant g and the total bilateral contract reduction amount of the system is as follows:
ΔSg=SXg·ΔS
the t-period bilateral contract output reduction of the power plant g is non-negative constraint:
Figure BDA0001092896220000044
and (3) the three commons of the power plant g share the reduction constraint:
Figure BDA0001092896220000045
and (3) the relationship constraint of the three-fair simultaneous reduction amount of the power plant g and the total three-fair simultaneous reduction amount of the system is as follows:
ΔJg=JXg·ΔJ
the three commons of the power plant g are reduced by the same amount and are not restricted to negative:
ΔJg≥0
and (3) constraining the relationship between the unit output and the unit transaction output:
Pi,t=Pi,S,t+Pi,J,t+Pi,z,t
limiting and restricting the transaction output part of the unit: pi,S,t≥0
And (3) limiting and restricting a third output part of the unit: pi,J,t≥0
And (3) limiting and restricting other components of unit output: pi,z,t≥0
Thermal powerAnd (3) unit output limit constraint: pi,minIi,t≤Pi,t≤Pi,maxIi,t
Wherein, Ii,tIs the running state of the thermal power generating unit at t time period Ii,t-1The operating state of the thermal power generating unit in the t-1 period is shown,
Figure BDA0001092896220000046
for the powered on time and the powered off time to the end of the t-1 period,
Figure BDA0001092896220000047
minimum boot time and minimum downtime, P, respectivelyi,tPlanned output, P, of thermal power generating unit for time period ti,max、Pi,minRespectively an upper limit and a lower limit of output of the thermal power generating unit, Pi,S,tIs the double-side output part of the unit output, Pi,J,tIs the part of three public forces, P, in the unit forcei,z,tFor other components of the unit, excluding the trade force and the three-force portions, LtFor t period system load, RUt、RDtUp and down rotation standby requirements, R, respectively, for a period of ti、DiThe climbing speed and the landslide speed u of the thermal power generating unit i respectivelyi,tStarting variable v for thermal power generating unit ii,tFor the i outage variable of the thermal power generating unit, delta S is the total bilateral contract reduction amount of the system, delta J is the total three-common contract reduction amount of the system, and delta SgFor a bilateral contract reduction of the power plant g, Δ JgFor the three-fair simultaneous reduction of the g of the power plant, Pg,0,tFor a bilateral contract for the g period t of the power plant, Jg,0Is the same electric quantity of the power plant g in three commons, NgNumbering all units of the power plant g; SXgFor the power plant g bilateral contract quantity reduction to account for the proportion factor of the total contract quantity reduction, JXgAnd reducing the scale factor accounting for the total contract reduction amount for the three-common-contract equivalent amount of the power plant g.
Preferably, in the step (5), the optimization objective includes: minimizing the amount of contract adjustments and minimizing the cost of power generation; the minimized contract adjustment amount has higher priority;
the objective function is:
F=M1*(W1*ΔS+W2*ΔJ)+M2*C(P)
Figure BDA0001092896220000051
Figure BDA0001092896220000052
Figure BDA0001092896220000053
wherein, Pi,S_0,tAnd (5) providing power for the unit i in the bilateral contract of the time t.
Preferably, in the step (5), the constraint condition includes:
and power balance constraint:
Figure BDA0001092896220000054
and (3) constraint of starting variables and stopping dynamic variables of the thermal power generating unit:
Ii,t-Ii,t-1=ui,t-vi,t,ui,t+vi,t≤1
minimum on-off time constraint:
Figure BDA0001092896220000061
vertical rotation standby restraint:
Figure BDA0001092896220000062
and (3) climbing restraint: pi,t-Pi,t-1≤Ri(1+Ii,t-1-Ii,t)+Pi,min(2-Ii,t-1-Ii,t)
Landslide restraint: pi,t-1-Pi,t≤Di(1-Ii,t-1+Ii,t)+Pi,min(2-Ii,t-1-Ii,t)
And (3) bilateral contract reduction constraint of the power plant g:
Figure BDA0001092896220000063
and (3) the relation constraint of the bilateral contract reduction amount of the power plant g and the total bilateral contract reduction amount of the system is as follows:
ΔSg=SXg·ΔS
the t-period bilateral contract output reduction non-negative constraint of the unit i:
Pi,S_0,t≥Pi,S,t
and (3) the three commons of the power plant g share the reduction constraint:
Figure BDA0001092896220000064
and (3) the relationship constraint of the three-fair simultaneous reduction amount of the power plant g and the total three-fair simultaneous reduction amount of the system is as follows:
ΔJg=JXg·ΔJ
the three commons of the power plant g are reduced by the same amount and are not restricted to negative:
ΔJg≥0
and (3) constraining the relationship between the unit output and the unit transaction output:
Pi,t=Pi,S,t+Pi,J,t+Pi,z,t
limiting and restricting the transaction output part of the unit: pi,S,t≥0
And (3) limiting and restricting a third output part of the unit: pi,J,t≥0
And (3) limiting and restricting other components of unit output: pi,z,t≥0
And (3) output limit constraint of the thermal power generating unit: pi,minIi,t≤Pi,t≤Pi,maxIi,t
Wherein, Ii,tIs the running state of the thermal power generating unit at t time period Ii,i-1The operating state of the thermal power generating unit in the t-1 period is shown,
Figure BDA0001092896220000071
is started up at the end of t-1 timeThe time and the time of the down time,
Figure BDA0001092896220000072
minimum boot time and minimum downtime, P, respectivelyi,tPlanned output, P, of thermal power generating unit for time period ti,max、Pi,minRespectively an upper limit and a lower limit of output of the thermal power generating unit, Pi,S,tIs the double-side output part of the unit output, Pi,J,tIs the part of three public forces, P, in the unit forcei,z,tFor other components of the unit, excluding the trade force and the three-force portions, LtFor t period system load, RUt、RDtUp and down rotation standby requirements, R, respectively, for a period of ti、DiThe climbing speed and the landslide speed u of the thermal power generating unit i respectivelyi,tStarting variable v for thermal power generating unit ii,tFor the i outage variable of the thermal power generating unit, delta S is the total bilateral contract reduction amount of the system, delta J is the total three-common contract reduction amount of the system, and delta SgFor a bilateral contract reduction of the power plant g, Δ JgFor the three-fair simultaneous reduction of the g of the power plant, Pi,S_0,tThe output is specified for the bilateral contract of the unit g time period t, Jg,0Is the same electric quantity of the power plant g in three commons, NgNumbering all units of the power plant g; SXgFor the power plant g bilateral contract quantity reduction to account for the proportion factor of the total contract quantity reduction, JXgAnd reducing the scale factor accounting for the total contract reduction amount for the three-common-contract equivalent amount of the power plant g.
Compared with the closest prior art, the technical scheme of the invention has the following beneficial effects:
the invention considers a coordination optimization method of double-side power contract participating in day-ahead power electric quantity balance, the double-side power contract and the three-public power contract belong to two types of contracts with different properties, and the two types of contracts are respectively modeled in an optimization model.
Two types of bilateral power contracts are considered, namely a power plant bilateral power contract and a unit bilateral power contract, and different bilateral power contract models are adopted when the contract types are different.
The unit output is decomposed into a bilateral output part, a three-output part and other output parts, and reasonable constraint is established to enable the bilateral actual output to be not greater than the corresponding contractual specified output, so that the non-negative contract adjustment quantity in the objective function is ensured, the objective function is prevented from being changed into non-linearity by adopting an absolute value or a square value, and the solving difficulty is reduced.
The two types of optimization targets are effectively combined, and the two targets with different priorities are combined together through penalty factors, so that the transaction contract can be executed maximally, the power generation cost can be minimized on the basis, the transaction can be completed maximally at the lowest cost, and the economy of transaction execution and power generation is well considered.
Drawings
Fig. 1 is a comparison graph of the completion of the bilateral transaction under the output curves of the two units provided by the invention.
Fig. 2 is a flow chart of a day-ahead power electric quantity balancing technical scheme considering a bilateral power contract provided by the invention.
FIG. 3 is a result of a day-ahead crew portfolio provided by the present invention that considers a bilateral power contract for a power plant.
Fig. 4 is a result of the day-ahead unit combination considering the bilateral electricity contract provided by the present invention.
FIG. 5 is a comparison of the actual bilateral output of the power plant provided by the present invention with a bilateral contract output curve.
Fig. 6 shows the combination result of the day ahead unit considering the unit bilateral power contract provided by the present invention.
Detailed Description
The following describes embodiments of the present invention in further detail with reference to the accompanying drawings.
The first step is as follows: and obtaining the three-fair electricity quantity of each power plant in the next day according to the electricity quantity decomposed to the day by the three-fair electricity quantity of the year and the month. According to a bilateral power contract submitted by a power generation and utilization party, if the power contract specifies a bilateral power curve of a power plant in the next day, obtaining bilateral output data of the power plant; and if the power contract specifies a bilateral power curve of the unit in the future day, obtaining bilateral output data of the unit so as to determine constraint conditions related to transactions in the day-ahead power electric quantity balance.
And obtaining three electric quantities required to be completed by each power plant. For example, according to the medium-and-long-term contract decomposition result of the schedule, the power plant a will have three common electric quantities of 2000MWH in 24 hours in the future, and then three common electric quantities need to be completed as far as possible when the electric power electric quantities are balanced in the future.
And obtaining a bilateral power contract curve. If the power plant output curve of 24 hours in the future is stipulated in the contract, the power plant bilateral curve needs to be completed as far as possible in the day-ahead power and electricity balance; if a 24-hour future output curve of a unit in a power plant is specified in a contract, the day-ahead power and electricity balance needs to complete the unit bilateral curve as far as possible.
(3) And determining a contract reduction factor of each power plant. That is, a scaling factor is determined for the contract curtailment amount for each power plant to the total contract adjustment amount for the system. The bilateral reduction factor is the proportion of the bilateral contract reduction of the power plant to the total bilateral reduction of the system, and the three-party reduction factor is the proportion of the three-party contract reduction of the power plant to the total three-party reduction of the system. Bilateral contract reduction factor SX for each power plant on the same dayiThree-common same reduction factor JXiComprises the following steps:
Figure BDA0001092896220000091
wherein S isiBilateral contract power quantity, J, determined for the i bilateral power contract of the power plantiThe three-fair electric quantity of the power plant i is decomposed to the same day.
And according to the contract quantity of each power plant, each power plant performs contract reduction in proportion to the contract quantity of each power plant. The power plant three-public contract is an electric quantity contract. The statistical method of the bilateral contract electric quantity of the power plant comprises the following steps:
1) and the bilateral contract is the power contract of the power plant, the electric quantity of one day is counted according to a curve specified by the power contract of the power plant, and if the given curve is a 24-point curve, the 24-point data is always added to the bilateral contract quantity of the power plant for one day.
2) The double-side contract is a power contract of the unit power plant, and meanwhile, the contract electric quantity of the power plant is calculated through the following steps. Step 1: and (3) counting the daily electric quantity of each unit according to a curve specified by the unit electric power contract, and if the given curve is a 24-point curve, adding 24 points of data to obtain the double-side contract quantity of the unit for one day. Step 2: and (4) summing the bilateral contract quantities (obtained in the step 1) of all the units in one power plant, and counting the bilateral contract quantities of one power plant.
Assuming A, B, C three power plants, two types of contract reduction factors of each power plant are counted according to two types of contract electric quantity.
Power plant Bilateral contract (MWH) Three public contract (MWH) Bilateral reduction factor Three common reduction factor
A 6000 2000 2/5 1/2
B 5000 1000 1/3 1/4
C 4000 1000 4/15 1/4
According to the information in the table, a scaling factor 6 of the total bilateral reduction amount occupied by bilateral contract reduction of each power plant in the constraint related to bilateral contract reduction in the optimization problem can be determined: 5: 4, the scaling factors for the bilateral contract reductions of the plant A, B, C are SXA、SXB、SXCThen, then
Figure BDA0001092896220000101
And (3) a scaling factor 2 of the total bilateral reduction amount of the three-fair equivalent reduction of each power plant in the constraint related to the three-fair equivalent: 1: 1, the scaling factors of the bilateral contract reductions of the plant A, B, C are JX respectivelyA、JXB、JXCThen, then
Figure BDA0001092896220000102
The second step is that: and obtaining the predicted data of the system load in the future 24 hours. And acquiring basic parameters and coal consumption characteristics of each thermal power generating unit so as to determine power balance constraint and related constraint conditions of the thermal power generating units. The related constraints of the thermal power generating unit comprise unit output limit and minimum on-off time constraint.
The third step: and according to the bilateral contract type determined in the first step, if the bilateral power contract is a power plant power contract, jumping to the fourth step. And if the bilateral power contract is the unit power contract, jumping to the fifth step.
The fourth step: and establishing a power and electricity balance coordination optimization model considering the power plant power contract day ahead according to the data of the first step and the data of the second step, and determining an objective function and constraint conditions according to the operation requirements and purposes, wherein the constraint conditions comprise thermal power unit related constraints, contract related constraints and system operation constraints.
The optimization objective includes two: the method comprises the steps of minimizing a contract adjustment amount and minimizing the power generation cost, wherein the minimized contract adjustment amount is higher in priority, and the minimized power generation cost is the second time, namely, the power generation cost is minimized on the premise of executing the transaction to the maximum extent possible. To achieve this, a large penalty factor will be introduced to coordinate the contribution of the two objectives in the overall objective. When the bilateral contract is the electric contract, the bilateral contract and the three contracts are used as two independent electric quantity components for modeling respectively. In order to establish a linear model, the contract adjustment quantity in the objective function needs to be positive, the implementation method is that the unit output is decomposed into an output part corresponding to bilateral transaction, an output part corresponding to three-centimeter electricity and other parts, and the bilateral output of the power plant is restrained to be not more than the bilateral contract regulation output of the power plant in the restraint process.
F=M1*(W1*ΔS+W2*ΔJ)+M2*C(P)
Figure BDA0001092896220000111
Figure BDA0001092896220000112
Figure BDA0001092896220000113
H is the total number of the power plant, N is the total number of the thermal power generating unit, T is the total number of the time period, Delta S is the bilateral contract reduction amount, Delta J is the three-fair contract reduction amount, C (P) is the system power generation cost, M1Reducing penalty factors for contracts, M2For power generation cost weighting, W1Weight reduction for bilateral contracts, W2For three co-ordinates to reduce the weight, Pg,0,tFor a bilateral contract-defined output, P, of the power plant at time period gi,S,tRepresents the output of the unit corresponding to the trade contract part, Jg,0Is the same electric quantity of the power plant g in three commons, Pi,J,tRepresenting the output of the unit corresponding to the trade contract part, fi,tFor the operating costs of thermal power units, SUi,t、SDi,tThe starting cost and the shutdown cost of the thermal power generating unit are respectively.
M1、M2Is selected fromThe selection is based on the priority of transaction execution, and when the transaction needs to be executed maximally, the contract reduction weight is selected to be larger, such as M1Get 104,M2And 1 is taken, namely when the contract reduction is reduced by 1, the objective function can be reduced by 10000, the influence is far greater than the influence of the power generation cost change on the objective function, the minimized contract reduction can be guaranteed preferentially, and then the power generation cost is reduced on the premise of the minimum contract reduction.
W1、W2The selection is based on the priority of bilateral and three-public execution, when bilateral needs to be executed preferentially, W is1Take the larger value, otherwise, W2Take the larger value. E.g. M1Get 104,M2When taking 1, W1Take 4, W 21 is selected, when the bilateral contract reduction is reduced by 1, the target function can be reduced by 40000, when the three-fair contract reduction is reduced by 1, the target function can be reduced by 10000, so that the influence of the bilateral contract reduction is far greater than the influence of the three-fair contract reduction, the influence of the contract reduction is far greater than the influence of the power generation cost change on the target function, the minimized contract reduction can be guaranteed preferentially, and the power generation cost is reduced on the premise.
The constraint conditions include power balance, minimum on-off time, standby, ramp-up, unit output limit and contract-related constraints as follows:
and power balance constraint:
Figure BDA0001092896220000121
and (3) constraint of starting variables and stopping dynamic variables of the thermal power generating unit:
Ii,t-Ii,t-1=ui,t-vi,t,ui,t+vi,t≤1
minimum on-off time constraint:
Figure BDA0001092896220000122
vertical rotation standby restraint:
Figure BDA0001092896220000123
and (3) climbing restraint: pi,t-Pi,t-1≤Ri(1+Ii,t-1-Ii,t)+Pi,min(2-Ii,t-1-Ii,t)
Landslide restraint: pi,t-1-Pi,t≤Di(1-Ii,t-1+Ii,t)+Pi,min(2-Ii,t-1-Ii,t)
And (3) bilateral contract reduction constraint of the power plant g:
Figure BDA0001092896220000124
and (3) the relation constraint of the bilateral contract reduction amount of the power plant g and the total bilateral contract reduction amount of the system is as follows:
ΔSg=SXg·ΔS
the t-period bilateral contract output reduction of the power plant g is non-negative constraint:
Figure BDA0001092896220000125
and (3) the three commons of the power plant g share the reduction constraint:
Figure BDA0001092896220000126
and (3) the relationship constraint of the three-fair simultaneous reduction amount of the power plant g and the total three-fair simultaneous reduction amount of the system is as follows:
ΔJg=JXg·ΔJ
the three commons of the power plant g are reduced by the same amount and are not restricted to negative:
ΔJg≥0
and (3) constraining the relationship between the unit output and the unit transaction output:
Pi,t=Pi,S,t+Pi,J,t+Pi,z,t
limiting and restricting the transaction output part of the unit: pi,S,t≥0
And (3) limiting and restricting a third output part of the unit: pi,J,t≥0
And (3) limiting and restricting other components of unit output: pi,z,t≥0
And (3) output limit constraint of the thermal power generating unit: pi,minIi,t≤Pi,t≤Pi,maxIi,t
Wherein, Ii,tIs the running state of the thermal power generating unit at t time period Ii,i-1The operating state of the thermal power generating unit in the t-1 period is shown,
Figure BDA0001092896220000131
for the powered on time and the powered off time to the end of the t-1 period,
Figure BDA0001092896220000132
minimum boot time and minimum downtime, P, respectivelyi,tPlanned output, P, of thermal power generating unit for time period ti,max、Pi,minRespectively an upper limit and a lower limit of output of the thermal power generating unit, Pi,S,tIs the double-side output part of the unit output, Pi,J,tIs the part of three public forces, P, in the unit forcei,z,tFor other components of the unit, excluding the trade force and the three-force portions, LtFor t period system load, RUt、RDtUp and down rotation standby requirements, R, respectively, for a period of ti、DiThe climbing speed and the landslide speed u of the thermal power generating unit i respectivelyi,tStarting variable v for thermal power generating unit ii,tFor the i outage variable of the thermal power generating unit, delta S is the total bilateral contract reduction amount of the system, delta J is the total three-common contract reduction amount of the system, and delta SgFor a bilateral contract reduction of the power plant g, Δ JgFor the three-fair simultaneous reduction of the g of the power plant, Pg,0,tFor a bilateral contract for the g period t of the power plant, Jg,0Is the same electric quantity of the power plant g in three commons, NgNumbering all units of the power plant g; SXgFor the power plant g bilateral contract quantity reduction to account for the proportion factor of the total contract quantity reduction, JXgAnd reducing the scale factor accounting for the total contract reduction amount for the three-common-contract equivalent amount of the power plant g.
Contract reduction between power plants needs to be coordinated, so the relation constraint of the power plant reduction and the system contract reduction is added in the constraint. The non-negative constraint of the output reduction of the bilateral contract of the power plant g and the non-negative constraint of the reduction of the three-fair contract can prevent the reduction of part of the power plant from being negative and the reduction of part of the power plant from being positive, which leads to the cancellation of the reduction of the contract between the power plants when the total reduction of the contract of the power plant is increased. For a power plant, the contract completion rate of the whole power plant is concerned, the contract reduction is distributed among the units in the power plant under the influence of the economy of the generator sets, the aim of minimizing the power generation cost is fulfilled, the power plant is guaranteed to complete the power generation contract to the maximum extent at the lowest cost, the benefit of the power plant is met, and the power generation cost is reduced.
The fifth step: and establishing a day-ahead power and electric quantity balance coordination optimization model considering the power contract of the unit according to the data of the first step and the data of the second step, and determining an objective function and constraint conditions according to the operation requirement and the purpose, wherein the constraint conditions comprise thermal power unit related constraints, contract related constraints and system operation constraints.
The optimization objective includes two: the method comprises the steps of minimizing a contract adjustment amount and minimizing the power generation cost, wherein the minimized contract adjustment amount is higher in priority, and the minimized power generation cost is the second time, namely, the power generation cost is minimized on the premise of executing the transaction to the maximum extent possible. To achieve this, a large penalty factor will be introduced to coordinate the contribution of the two objectives in the overall objective. When the bilateral contract is the electric contract, the bilateral contract and the three contracts are used as two independent electric quantity components for modeling respectively. In order to establish a linear model, the contract adjustment quantity in the objective function needs to be positive, the method is implemented by decomposing the unit output into an output part corresponding to bilateral transaction, an output part corresponding to three-centimeter electricity and other parts, and constraining the unit bilateral output part not to be larger than the unit bilateral contract specified output in constraint.
F=M1*(W1*ΔS+W2*ΔJ)+M2*C(P)
Figure BDA0001092896220000141
Figure BDA0001092896220000142
Figure BDA0001092896220000143
H is the total number of the power plant, N is the total number of the thermal power generating unit, T is the total number of the time period, Delta S is the bilateral contract reduction amount, Delta J is the three-fair contract reduction amount, C (P) is the system power generation cost, M1Reducing penalty factors for contracts, M2For power generation cost weighting, W1Weight reduction for bilateral contracts, W2For three co-ordinates to reduce the weight, Pi,S_0,tFor the bilateral contract regulation of the unit in the period i t, the output Pi,S,tRepresents the output of the unit corresponding to the trade contract part, Jg,0Is the same electric quantity of the power plant g in three commons, Pi,J,tRepresenting the output of the unit corresponding to the trade contract part, fi,tFor the operating costs of thermal power units, SUi,t、SDi,tThe starting cost and the shutdown cost of the thermal power generating unit are respectively.
M1、M2The selection is based on the priority of transaction execution, and when the transaction needs to be executed maximally, the contract-reduction weight is greater, such as M1Get 104,M2And 1 is taken, namely when the contract reduction is reduced by 1, the objective function can be reduced by 10000, the influence is far greater than the influence of the power generation cost change on the objective function, the minimized contract reduction can be guaranteed preferentially, and then the power generation cost is reduced on the premise of the minimum contract reduction.
W1、W2The selection is based on the priority of bilateral and three-public execution, when bilateral needs to be executed preferentially, W is1Take the larger value, otherwise, W2Take the larger value. E.g. M1Get 104,M2When taking 1, W1Take 4, W2If 1 is selected, the target function can be reduced by 40000 when the reduction of the bilateral contract is reduced by 1, and the target function can be reduced by 10000 when the reduction of the three-common contract is reduced by 1, so that the influence of the reduction of the bilateral contract is far greater than that of the bilateral contractThe influence of the three-common contract is reduced, the influence of the contract reduction is far greater than the influence of the power generation cost change on the objective function, the minimized contract reduction can be guaranteed preferentially, and then the power generation cost is reduced on the premise.
The constraint conditions include power balance, minimum on-off time, standby, ramp-up, unit output limit and contract-related constraints as follows:
and power balance constraint:
Figure BDA0001092896220000151
and (3) constraint of starting variables and stopping dynamic variables of the thermal power generating unit:
Ii,t-Ii,t-1=ui,t-vi,t,ui,t+vi,t≤1
minimum on-off time constraint:
Figure BDA0001092896220000152
vertical rotation standby restraint:
Figure BDA0001092896220000153
and (3) climbing restraint: pi,t-Pi,t-1≤Ri(1+Ii,t-1-Ii,t)+Pi,min(2-Ii,t-1-Ii,t)
Landslide restraint: pi,t-1-Pi,t≤Di(1-Ii,t-1+Ii,t)+Pi,min(2-Ii,t-1-Ii,t)
And (3) bilateral contract reduction constraint of the power plant g:
Figure BDA0001092896220000154
and (3) the relation constraint of the bilateral contract reduction amount of the power plant g and the total bilateral contract reduction amount of the system is as follows:
ΔSg=SXg·ΔS
the t-period bilateral contract output reduction non-negative constraint of the unit i:
Pi,S_0,t≥Pi,S,t
and (3) the three commons of the power plant g share the reduction constraint:
Figure BDA0001092896220000161
and (3) the relationship constraint of the three-fair simultaneous reduction amount of the power plant g and the total three-fair simultaneous reduction amount of the system is as follows:
ΔJg=JXg·ΔJ
the three commons of the power plant g are reduced by the same amount and are not restricted to negative:
ΔJg≥0
and (3) constraining the relationship between the unit output and the unit transaction output:
Pi,t=Pi,S,t+Pi,J,t+Pi,z,t
limiting and restricting the transaction output part of the unit: pi,S,t≥0
And (3) limiting and restricting a third output part of the unit: pi,J,t≥0
And (3) limiting and restricting other components of unit output: pi,z,t≥0
And (3) output limit constraint of the thermal power generating unit: pi,minIi,t≤Pi,t≤Pi,maxIi,t
Wherein, Ii,tIs the running state of the thermal power generating unit at t time period Ii,t-1The operating state of the thermal power generating unit in the t-1 period is shown,
Figure BDA0001092896220000162
for the powered on time and the powered off time to the end of the t-1 period,
Figure BDA0001092896220000163
minimum boot time and minimum downtime, P, respectivelyi,tPlanned output, P, of thermal power generating unit for time period ti,max、Pi,minRespectively an upper limit and a lower limit of output of the thermal power generating unit, Pi,S,tIs the double-side output part of the unit output, Pi,J,tIs the part of three public forces, P, in the unit forcei,z,tFor other components of the unit, excluding the trade force and the three-force portions, LtFor t period system load, RUt、RDtUp and down rotation standby requirements, R, respectively, for a period of ti、DiThe climbing speed and the landslide speed u of the thermal power generating unit i respectivelyi,tStarting variable v for thermal power generating unit ii,tFor the i outage variable of the thermal power generating unit, delta S is the total bilateral contract reduction amount of the system, delta J is the total three-common contract reduction amount of the system, and delta SgFor a bilateral contract reduction of the power plant g, Δ JgFor the three-fair simultaneous reduction of the g of the power plant, Pi,S_0,tThe output is specified for the bilateral contract of the unit g time period t, Jg,0Is the same electric quantity of the power plant g in three commons, NgNumbering all units of the power plant g; SXgFor the power plant g bilateral contract quantity reduction to account for the proportion factor of the total contract quantity reduction, JXgAnd reducing the scale factor accounting for the total contract reduction amount for the three-common-contract equivalent amount of the power plant g.
Contract reduction between power plants needs to be coordinated, so the relation constraint of the power plant reduction and the system contract reduction is added in the constraint. The output reduction non-negative constraint of the unit bilateral contract and the three-common contract reduction non-negative constraint of the power plant can prevent the reduction of part of the power plant from being negative and the reduction of part of the power plant from being positive, so that the reduction of the contract between the power plants is cancelled when the contract reduction of the power plant is totally increased. For a power plant, the contract completion rate of the whole power plant is concerned, the contract reduction is distributed among the units in the power plant under the influence of the economy of the generator sets, the aim of minimizing the power generation cost is fulfilled, and the power plant is guaranteed to complete the power generation contract at the lowest cost, load the benefits of the power plant and reduce the power generation cost.
And a sixth step: and according to the type of the bilateral power contract, simultaneously solving the mixed integer programming problem determined in the fourth step for the power contract of the power plant, and simultaneously solving the mixed integer programming problem determined in the fourth step for the unit power contract.
The seventh step: and taking the result obtained in the sixth step as a power grid unit scheduling scheme, determining a starting and stopping plan and a unit output plan of the unit and a power plant reduction amount in 24 hours in the future, ensuring bilateral transaction and execution of three public electric quantities as far as possible, and improving the safety and the economy of power grid operation.
As shown in fig. 1: in the figure, the completion electric quantity of the output curve 1 is the same as that of the output curve 2, but the bilateral transaction completion conditions are different due to different curve shapes, the output curve 1 completes the bilateral power contract, but the output of the output curve 2 is lower than that of the bilateral power contract in some time periods, so that the bilateral contract is not completely completed.
3 thermal power plants, 7 thermal power generating units, the parameters of the power generating units are shown in a table 1, the load prediction data is shown in a table 2, the bilateral power contract information of the power plant is shown in a table 3, and the three public contract information of the power plant is shown in a table 4.
TABLE 1 Power plant Unit parameters
Figure BDA0001092896220000181
TABLE 2 load forecast data
Time period load/MW Time period load/MW
1 575.19 13 642.18
2 565.15 14 643.6
3 558.67 15 648.86
4 554.73 16 655.79
5 555.06 17 656
6 560.48 18 638.74
7 573.39 19 645.97
8 590.4 20 632.35
9 605.56 21 637.31
10 617.2 22 627.14
11 628.61 23 601.05
12 636.1 24 596.75
TABLE 3 bilateral Power contract information for Power plants
Figure BDA0001092896220000182
Figure BDA0001092896220000191
Statistics of plant contract information from Table 3 As shown in Table 4, the contract curtailment factor for each plant is calculated. Selecting M1M 210000 and 1, W respectively1、W2And 4 and 1, preferably ensuring bilateral execution, solving the day-ahead power and electric quantity balance optimization problem considering the bilateral power contract of the power plant, and obtaining a unit combination result as shown in the figure 3, comparing an actual bilateral output curve of the power plant with a power plant contract curve as shown in the figure 5, a unit bilateral output result as shown in the table 5, a unit output as shown in the table 6 and a contract reduction result as shown in the table 7.
TABLE 4 contract reduction factor for Power plants
Power plant Bilateral contract/MW Three public contract/MW Bilateral reduction factor Three common reduction factor
A 6447 1611.75 0.5 0.5
B 1074.5 268.63 0.083 0.083
C 5372.5 1343.13 0.417 0.417
Accumulation 12894 3223.51 1 1
Bilateral contract reduction factor:
Figure BDA0001092896220000192
Figure BDA0001092896220000201
Figure BDA0001092896220000202
similarly, three commonalities share the reduction factor:
Figure BDA0001092896220000203
Figure BDA0001092896220000204
Figure BDA0001092896220000205
TABLE 5 double-side output results of the unit
Figure BDA0001092896220000206
TABLE 6 Unit output results
Figure BDA0001092896220000211
TABLE 7 contract reduction of Power plants and Power Generation cost results
Figure BDA0001092896220000212
Comparing a unit combination result (shown in figure 3) considering a bilateral power contract of a power plant with a unit combination result (shown in figure 4) considering a bilateral power contract, wherein the former is increased by 3 time intervals compared with the latter by the unit G2, the unit G4 is increased by 6 time intervals, the unit G5 is increased by 6 time intervals, and the unit G6 is decreased by 1 time interval.
According to the attached figure 5, it can be seen that when the bilateral contract stipulates the bilateral output curve of the power plant, the actual bilateral output curve of the power plant is as close to the bilateral contract curve as possible, and the optimization result is that the bilateral contract is executed as far as possible.
The invention effectively combines two types of optimization targets, hopes to realize the transaction completion with the lowest cost maximization, and considers the economy of transaction execution and power generation in the optimization, so compared with the power generation cost under two conditions, one condition is that the minimized contract adjustment amount and the minimized power generation cost are optimized in a combined manner, and the other condition is that the objective function only comprises the minimized contract adjustment amount, and the optimization results are shown in a table 8.
TABLE 8 comparison of Power Generation costs
Figure BDA0001092896220000221
Table 8 shows that, compared with the single-target optimization using the contract adjustment amount, the joint optimization using the contract adjustment amount and the power generation cost can ensure that the trade power is executed to the maximum extent, and the reduction of the trade is consistent, but the two different optimizations have significant influence on the power generation cost, the joint optimization considers the power generation economy, and the power generation cost is reduced by 27.68% compared with the single-target optimization, which indicates that the joint optimization is effective.
It should be noted that the above-mentioned embodiments are only for illustrating the technical solutions of the present application and not for limiting the scope of protection thereof, and although the present application is described in detail with reference to the above-mentioned embodiments, those skilled in the art should understand that after reading the present application, they can make various changes, modifications or equivalents to the specific embodiments of the application, but these changes, modifications or equivalents are all within the scope of protection of the claims to be filed.

Claims (5)

1. A coordination optimization method for enabling a bilateral power contract to participate in day-ahead power electric quantity balance is characterized by comprising the following steps:
(1) acquiring the three-fair-contract same electric quantity of each power plant in the future one day, acquiring bilateral power contract data of the power plants or bilateral power contract data of the units according to the bilateral contract types, and determining contract reduction factors of each power plant;
(2) acquiring system load prediction data in the future 24 hours, basic parameters of each thermal power generating unit and coal consumption characteristics;
(3) if the bilateral power contract is the bilateral power contract of the power plant, turning to the step (4); if the bilateral power contract is the unit bilateral power contract, turning to the step (5);
(4) the method for establishing the day-ahead power and electric quantity balance coordination optimization model considering the bilateral power contract of the power plant comprises the following steps: determining an objective function and setting constraint conditions, and turning to the step (6);
(5) the method for establishing the day-ahead power and electric quantity balance coordination optimization model considering the unit bilateral power contract comprises the following steps: determining an objective function and setting constraint conditions, and turning to the step (6);
(6) solving a mixed integer programming problem;
(7) determining a starting and stopping plan and a unit output plan of a unit and a power plant reduction amount in 24 hours in the future;
in the step (4), the optimization objective includes: minimizing the amount of contract adjustments and minimizing the cost of power generation; the minimized contract adjustment amount has higher priority;
the objective function is:
F=M1*(W1*ΔS+W2*ΔJ)+M2*C(P)
Figure FDA0003214299270000011
Figure FDA0003214299270000021
Figure FDA0003214299270000022
wherein F is the value of the objective function, M1Reducing penalty factors for contracts, M2For power generation cost weighting, W1Weight reduction for bilateral contracts, W2The weight is reduced for three fair contracts, Δ S is reduced for two-sided contracts, Δ J is reduced for three fair contracts, C (P) is the system power generation cost, Pg,0,tThe bilateral contract specified output of the power plant in the period of T is g, H is the total number of the power plant, N is the total number of the thermal power generating units, T is the total number of the period of T, P isi,S,tRepresenting the double-sided part of the unit output, Jg,0Is g three common and same electric quantity of the power plant, Pi,J,tRepresenting the three-revolution-force part of the unit force, fi,tFor the operating costs of thermal power units, SUi,t、SDi,tRespectively representing the starting cost and the shutdown cost of the thermal power generating unit;
in the step (5), the optimization objective includes: minimizing the amount of contract adjustments and minimizing the cost of power generation; the minimized contract adjustment amount has higher priority;
the objective function is:
F=M1*(W1*ΔS+W2*ΔJ)+M2*C(P)
Figure FDA0003214299270000023
Figure FDA0003214299270000024
Figure FDA0003214299270000025
wherein,Pi,S_0,tAnd (5) providing power for the unit i in the bilateral contract of the time t.
2. The coordinated optimization method according to claim 1, wherein the step (1) comprises:
a. the bilateral contract is a bilateral power contract of the power plant, and the bilateral contract amount of one day is counted according to a power plant output curve specified by the bilateral power contract of the power plant;
b. and the bilateral contract is the bilateral electric contract of the units, the electric quantity of each unit in one day is counted according to the unit output curve specified by the bilateral electric contract of the units, and the bilateral contract quantities of all the units in one power plant are added to obtain the bilateral contract quantity of one power plant.
3. The coordinated optimization method according to claim 1, wherein in the step (1), the reduction factor includes: the bilateral reduction factor is the proportion of the bilateral contract reduction amount of the power plant to the total bilateral reduction amount of the system; the three-fair reduction factor is the proportion of the three-fair reduction amount of the power plant to the total three-fair reduction amount of the system.
4. The coordinated optimization method according to claim 1, wherein in the step (4), the constraint conditions include:
and power balance constraint:
Figure FDA0003214299270000031
and (3) constraint of starting variables and stopping dynamic variables of the thermal power generating unit:
Ii,t-Ii,t-1=ui,t-vi,t,ui,t+vi,t≤1
minimum on-off time constraint:
Figure FDA0003214299270000032
vertical rotation standby restraint:
Figure FDA0003214299270000033
and (3) climbing restraint: pi,t-Pi,t-1≤Ri(1+Ii,t-1-Ii,t)+Pi,min(2-Ii,t-1-Ii,t)
Landslide restraint: pi,t-1-Pi,t≤Di(1-Ii,t-1+Ii,t)+Pi,min(2-Ii,t-1-Ii,t)
And (3) bilateral contract reduction constraint of the power plant g:
Figure FDA0003214299270000034
and (3) the relation constraint of the bilateral contract reduction amount of the power plant g and the total bilateral contract reduction amount of the system is as follows: delta Sg=SXg·ΔS
The t-period bilateral contract output reduction of the power plant g is non-negative constraint:
Figure FDA0003214299270000041
and (3) the three commons of the power plant g share the reduction constraint:
Figure FDA0003214299270000042
and (3) the relationship constraint of the three-fair simultaneous reduction amount of the power plant g and the total three-fair simultaneous reduction amount of the system is as follows:
ΔJg=JXg·ΔJ
the three commons of the power plant g are reduced by the same amount and are not restricted to negative:
ΔJg≥0
and (3) constraining the relationship between the unit output and the unit transaction output:
Pi,t=Pi,S,t+Pi,J,t+Pi,z,t
limiting and restricting the transaction output part of the unit: pi,S,t≥0
And (3) limiting and restricting a third output part of the unit: pi,J,t≥0
And (3) limiting and restricting other components of unit output: pi,z,t≥0
And (3) output limit constraint of the thermal power generating unit: pi,minIi,t≤Pi,t≤Pi,maxIi,t
Wherein, Ii,tIs the running state of the thermal power generating unit at t time period Ii,t-1The operating state of the thermal power generating unit in the t-1 period is shown,
Figure FDA0003214299270000043
for the powered on time and the powered off time to the end of the t-1 period,
Figure FDA0003214299270000044
Figure FDA0003214299270000045
minimum boot time and minimum downtime, P, respectivelyi,tPlanned output, P, of thermal power generating unit for time period ti,max、Pi,minRespectively an upper limit and a lower limit of output of the thermal power generating unit, Pi,S,tIs the double-side output part of the unit output, Pi,J,tIs the part of three public forces, P, in the unit forcei,z,tFor other components of the unit, excluding the trade force and the three-force portions, LtFor t period system load, RUt、RDtUp and down rotation standby requirements, R, respectively, for a period of ti、DiThe climbing speed and the landslide speed u of the thermal power generating unit i respectivelyi,tStarting variable v for thermal power generating unit ii,tFor the i outage variable of the thermal power generating unit, delta S is the total bilateral contract reduction amount of the system, delta J is the total three-common contract reduction amount of the system, and delta SgFor a bilateral contract reduction of the power plant g, Δ JgFor the three-fair simultaneous reduction of the g of the power plant, Pg,0,tFor a bilateral contract for the g period t of the power plant, Jg,0Is the same electric quantity of the power plant g in three commons, NgNumbering all units of the power plant g; SXgReducing total contract for bilateral contract amount of power plant gReduced scale factor, JXgAnd reducing the scale factor accounting for the total contract reduction amount for the three-common-contract equivalent amount of the power plant g.
5. The coordinated optimization method according to claim 1, wherein in the step (5), the constraint conditions include:
and power balance constraint:
Figure FDA0003214299270000051
and (3) constraint of starting variables and stopping dynamic variables of the thermal power generating unit:
Ii,t-Ii,t-1=ui,t-vi,t,ui,t+vi,t≤1
minimum on-off time constraint:
Figure FDA0003214299270000052
vertical rotation standby restraint:
Figure FDA0003214299270000053
and (3) climbing restraint: pi,t-Pi,t-1≤Ri(1+Ii,t-1-Ii,t)+Pi,min(2-Ii,t-1-Ii,t)
Landslide restraint: pi,t-1-Pi,t≤Di(1-Ii,t-1+Ii,t)+Pi,min(2-Ii,t-1-Ii,t)
And (3) bilateral contract reduction constraint of the power plant g:
Figure FDA0003214299270000054
and (3) the relation constraint of the bilateral contract reduction amount of the power plant g and the total bilateral contract reduction amount of the system is as follows: delta Sg=SXg·ΔS
The t-period bilateral contract output reduction non-negative constraint of the unit i:
Pi,S_0,t≥Pi,S,t
and (3) the three commons of the power plant g share the reduction constraint:
Figure FDA0003214299270000061
and (3) the relationship constraint of the three-fair simultaneous reduction amount of the power plant g and the total three-fair simultaneous reduction amount of the system is as follows:
ΔJg=JXg·ΔJ
the three commons of the power plant g are reduced by the same amount and are not restricted to negative:
ΔJg≥0
and (3) constraining the relationship between the unit output and the unit transaction output:
Pi,t=Pi,S,t+Pi,J,t+Pi,z,t
limiting and restricting the transaction output part of the unit: pi,S,t≥0
And (3) limiting and restricting a third output part of the unit: pi,J,t≥0
And (3) limiting and restricting other components of unit output: pi,z,t≥0
And (3) output limit constraint of the thermal power generating unit: pi,minIi,t≤Pi,t≤Pi,maxIi,t
Wherein, Ii,tIs the running state of the thermal power generating unit at t time period Ii,t-1The operating state of the thermal power generating unit in the t-1 period is shown,
Figure FDA0003214299270000062
for the powered on time and the powered off time to the end of the t-1 period,
Figure FDA0003214299270000063
Figure FDA0003214299270000064
minimum boot time and minimum downtime, P, respectivelyi,tPlanned output, P, of thermal power generating unit for time period ti,max、Pi,minRespectively an upper limit and a lower limit of output of the thermal power generating unit, Pi,S,tIs the double-side output part of the unit output, Pi,J,tIs the part of three public forces, P, in the unit forcei,z,tFor other components of the unit, excluding the trade force and the three-force portions, LtFor t period system load, RUt、RDtUp and down rotation standby requirements, R, respectively, for a period of ti、DiThe climbing speed and the landslide speed u of the thermal power generating unit i respectivelyi,tStarting variable v for thermal power generating unit ii,tFor the i outage variable of the thermal power generating unit, delta S is the total bilateral contract reduction amount of the system, delta J is the total three-common contract reduction amount of the system, and delta SgFor a bilateral contract reduction of the power plant g, Δ JgFor the three-fair simultaneous reduction of the g of the power plant, Pi,S_0,tThe output is specified for the bilateral contract of the unit g time period t, Jg,0Is the same electric quantity of the power plant g in three commons, NgNumbering all units of the power plant g; SXgFor the power plant g bilateral contract quantity reduction to account for the proportion factor of the total contract quantity reduction, JXgAnd reducing the scale factor accounting for the total contract reduction amount for the three-common-contract equivalent amount of the power plant g.
CN201610739177.2A 2016-08-26 2016-08-26 Coordination optimization method for bilateral power contract participating in day-ahead power and electric quantity balance Active CN107784375B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610739177.2A CN107784375B (en) 2016-08-26 2016-08-26 Coordination optimization method for bilateral power contract participating in day-ahead power and electric quantity balance

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610739177.2A CN107784375B (en) 2016-08-26 2016-08-26 Coordination optimization method for bilateral power contract participating in day-ahead power and electric quantity balance

Publications (2)

Publication Number Publication Date
CN107784375A CN107784375A (en) 2018-03-09
CN107784375B true CN107784375B (en) 2021-12-03

Family

ID=61439516

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610739177.2A Active CN107784375B (en) 2016-08-26 2016-08-26 Coordination optimization method for bilateral power contract participating in day-ahead power and electric quantity balance

Country Status (1)

Country Link
CN (1) CN107784375B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111898272A (en) * 2020-07-31 2020-11-06 国网青海省电力公司电力科学研究院 Electric quantity decomposition method and device

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103066626A (en) * 2012-10-31 2013-04-24 中国电力科学研究院 Multi-source coordinating day-ahead generation scheduling method
CN103855718A (en) * 2014-03-10 2014-06-11 东南大学 Scheduling method for pumped storage power station to participate in electric power system with wind power
CN105048491A (en) * 2015-06-29 2015-11-11 国电南瑞科技股份有限公司 Multi-stage wind power accepted range calculating method based on unit combination and economic dispatching

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103066626A (en) * 2012-10-31 2013-04-24 中国电力科学研究院 Multi-source coordinating day-ahead generation scheduling method
CN103855718A (en) * 2014-03-10 2014-06-11 东南大学 Scheduling method for pumped storage power station to participate in electric power system with wind power
CN105048491A (en) * 2015-06-29 2015-11-11 国电南瑞科技股份有限公司 Multi-stage wind power accepted range calculating method based on unit combination and economic dispatching

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
主动配电网的两阶段优化调度模型;庄慧敏等;《西南交通大学学报》;20151015;第50卷(第05期);正文第929-934页 *

Also Published As

Publication number Publication date
CN107784375A (en) 2018-03-09

Similar Documents

Publication Publication Date Title
CN107785885B (en) Coordination optimization method for participating in day-ahead power and electric quantity balance by bilateral electric quantity contract
Hirth et al. Balancing power and variable renewables: Three links
Riveros et al. Bidding strategies for virtual power plants considering CHPs and intermittent renewables
Valenzuela et al. Modeling and simulation of consumer response to dynamic pricing with enabled technologies
Wolak Measuring the competitiveness benefits of a transmission investment policy: The case of the Alberta electricity market
Rövekamp et al. Renewable electricity business models in a post feed-in tariff era
Kozlova et al. Modeling the effects of the new Russian capacity mechanism on renewable energy investments
Craine et al. Clean energy services for all: financing universal electrification
Zakeri et al. Pricing wind: a revenue adequate, cost recovering uniform price auction for electricity markets with intermittent generation
CN101794424A (en) Self-adaptive unified power market transaction method
CN106355285B (en) Power load prediction method based on parameter correction
Ding et al. Promoting energy-saving and environmentally friendly generation dispatching model in China: Phase development and case studies
Simoglou et al. Electricity market models and RES integration: The Greek case
Liu et al. Extending the bidding format to promote demand response
Fraunholz et al. Agent-based generation and storage expansion planning in interconnected electricity markets
Lee et al. Development of energy storage system scheduling algorithm for simultaneous self-consumption and demand response program participation in South Korea
CN110021932B (en) Bilateral participation peak shaving auxiliary service capacity selection and economic model construction method
Shen et al. Facilitating power grid decarbonization with distributed energy resources: Lessons from the united states
Bokkisam et al. Effective community energy management through transactive energy marketplace
Ali Development and Improvement of Renewable Energy Integrated with Energy Trading Schemes based on Advanced Optimization Approaches
Khajeh et al. Blockchain-based demand response using prosumer scheduling
CN107784375B (en) Coordination optimization method for bilateral power contract participating in day-ahead power and electric quantity balance
CN108695851A (en) A kind of generation optimization dispatching method based on overloaded partition electricity price
CN108075494B (en) Day-ahead unit combination method considering new energy consumption and transaction execution
Wong et al. Coordination of investor-owned DG capacity growth in distribution systems

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