CN104573875B - A kind of method of the power generating facilities and power grids optimization planning of low-carbon - Google Patents
A kind of method of the power generating facilities and power grids optimization planning of low-carbon Download PDFInfo
- Publication number
- CN104573875B CN104573875B CN201510041210.XA CN201510041210A CN104573875B CN 104573875 B CN104573875 B CN 104573875B CN 201510041210 A CN201510041210 A CN 201510041210A CN 104573875 B CN104573875 B CN 104573875B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- munder
- msubsup
- mtd
- 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
Links
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 165
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000005457 optimization Methods 0.000 title claims abstract description 20
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims abstract description 115
- 230000005611 electricity Effects 0.000 claims abstract description 93
- 238000011161 development Methods 0.000 claims abstract description 29
- 230000005540 biological transmission Effects 0.000 claims abstract description 28
- 230000008901 benefit Effects 0.000 claims abstract description 5
- 238000012546 transfer Methods 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims abstract description 4
- 238000010276 construction Methods 0.000 claims description 19
- 238000004519 manufacturing process Methods 0.000 claims description 10
- 230000033228 biological regulation Effects 0.000 claims description 9
- 238000005265 energy consumption Methods 0.000 claims description 7
- 230000005619 thermoelectricity Effects 0.000 claims description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 5
- 240000002853 Nelumbo nucifera Species 0.000 claims description 4
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 4
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 4
- 230000018109 developmental process Effects 0.000 description 24
- 208000028659 discharge Diseases 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000002146 bilateral effect Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000012010 growth Effects 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000033772 system development Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000003763 carbonization Methods 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 239000002737 fuel gas Substances 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 230000033001 locomotion Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 238000013433 optimization analysis Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- 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
-
- 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/80—Management or planning
- Y02P90/84—Greenhouse gas [GHG] management systems
-
- 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/80—Management or planning
- Y02P90/84—Greenhouse gas [GHG] management systems
- Y02P90/845—Inventory and reporting systems for greenhouse gases [GHG]
-
- 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)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Secondary Cells (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The present invention relates to a kind of power generating facilities and power grids Method for optimized planning of low-carbon, this method includes:All kinds of power supply capacitys in each region are divided into local capacity and the class of capacity two is sent outside, interregional point-to-point direct sending circuit is set up in equivalent network, with the flowing of clearly interregional electricity and the transfer of carbon emission, is easy to the metering of electricity consumption carbon emission;The decision variable of low-carbon power generating facilities and power grids Optimized model is set up, the low-carbon power generating facilities and power grids Optimized model being made up of object function and constraints is built with this, the constraints in Optimized model is converted into available matrix computations form;By solver processed, the low-carbon power generating facilities and power grids Optimized model is solved, the power generating facilities and power grids optimization planning of the low-carbon of all kinds of power supplys and transmission line of electricity in power system is tried to achieve so that power system has optimal economic benefit under low-carbon development model.The present invention provides technical support to improve the low-carbon level of Power System Planning.
Description
Technical field
The invention belongs to electric power system optimization analysis technical field, more particularly to power generating facilities and power grids Method for optimized planning.
Background technology
As energy problem increasingly protrudes with climate change problem, realize low carbon development, reduce excessively disappearing for fossil energy
Consumption is increasingly becoming the common objective of human society.The core of low carbon development is the change of technological innovation, system innovation and the view of development,
This will be related to production model, life style, values and readjusts, with national rights and interests closely related.
Power industry is as basic energy sector of China, while being also the maximum industry of CO2 emissions.Extremely
2011, the whole society of China carbon emission amount broke through 8,000,000,000 tons, per capita CO2Discharge also has been over the average level in the whole world, together
When power industry carbon emission amount break through 4,000,000,000 tons, account for the ratio of national carbon emission amount from 2006 37% rises to 50%.Electricity
Lixing industry is faced with politics and public opinion pressure huge in the world no matter on total emission volumn or under discharge development trend
Power.Power system has " carbon locking " effect, i.e., due to longer Years Of Services of equipment such as generating, transmissions of electricity so that power system
Carbon emission situation will be " locked " within a very long time.Therefore, low-carbon planning is carried out, for realizing electric power
The low carbon development of system is extremely important.
Power System Planning includes power source planning and two aspects of Electric Power Network Planning, is that one of power system development is important
Previous work, the basic goal of Power System Planning is, according to the result predicted in a certain period internal loading a certain region, to seek
A most economical electric power development scheme is sought, is allowed to enough meet the requirement of operational reliability.Existing Power System Planning method
Fail fully meter and the carbon emission and its corresponding financial cost of power system, it is impossible to fully meet electric power under low-carbon economy situation
The growth requirement of system.
For China's low carbon development target, the planning to power industry is also brought shadow by the difference of carbon emission metering method
Ring.Different carbon emissions calculate bore to the statistical discrepancy of region electric power carbon emission close to 50%, and this is completed to low carbon development target
The accounting of situation will produce tremendous influence.Because electric energy belongs to the secondary energy sources of cleaning, carbon is not produced in use
Discharge, carbon emission of the tradition based on macrostatistical approach, which is calculated, is adjusting the low-carbon target that such as per GDP carbon intensity declines
When will be unfair to electric energy output region, while being difficult to the enthusiasm for transferring energy input province energy-saving and emission-reduction work.Thus, compel
The carbon emission to be realized of being essential metering is from generating link to the transformation of electricity consumption link, and this transformation is also by for towards low-carbon target
Power generating facilities and power grids planning brings new challenge.The core concept of carbon emission metering based on consumption side is by the carbon emission of generating link
Share to electricity consumption link, the metering of electricity consumption correspondence carbon emission is realized from user side.It is existing to consider the carbon based on consumption side
Emission measurement angle and the trans-regional power transmission model of electricity consumption carbon emission constraint are strong with average generating (or electricity consumption) carbon emission
Spend the carbon emission density as electricity sent outside stream.Such a analysis method have ignored the actual conditions that Energy Base sends electric energy outside, difficult
To consider the bilateral electricity contract in multi area interconnection network between different zones.
When power supply architecture and interregional interconnection are not determined completely, the carbon emission metering bore based on electricity consumption link will
New change is brought for Power System Planning problem.It is multiple it is interregional there is electricity exchange, there is respective electricity consumption carbon emission
During constraint, allow whole system to carry out the extension of power generating facilities and power grids in the way of the lowest cost, will be low-carbon generation network optimization rule
Draw problem to be solved.
The content of the invention
The purpose of the present invention is the power system development demand being directed under low-carbon economy, proposes a kind of power supply electricity of low-carbon
Net Method for optimized planning.
A kind of power generating facilities and power grids Method for optimized planning of low-carbon proposed by the present invention, it is characterised in that:
This method comprises the following steps:
1) the low-carbon power generating facilities and power grids Optimized model being made up of object function and constraints is built, is specifically included:
All kinds of power supply capacitys in each region 1-1) are divided into two major classes, a class is local capacity, for balancing this
The capacity of ground load;Another kind of is to send capacity outside, the capacity for supplying other region loads;
Interregional point-to-point direct sending circuit is set up in equivalent network, is arranged with the flowing of clearly interregional electricity and carbon
The transfer put, is easy to the metering of electricity consumption carbon emission;
1-2) set up the decision variable of low-carbon power generating facilities and power grids Optimized model
The model is using year as decision-making time unit, and the decision variable in model mainly includes following three class:The first kind isRepresent that a and b represents area in the extension capacity of u class power supplys powered in a of y regions for region b loads, footnote
Domain;When a and b represent the same area, such asDuring for local capacity, show that such unit increases capacity newly and is used for equilibrium region a sheets
Ground load;When a and b represent different zones, show that such unit increases capacity newly and is used for equilibrium region b load, to send appearance outside
Amount;Equations of The Second Kind isRepresent y regions a to the newly-increased capacity of power transmission passage between the b of region.3rd class isFor
It is exerting oneself for the u class power supplys that region b loads are powered in y, typical day d under period t, in a of region;
1-3) set up low-carbon power generating facilities and power grids Optimized model object function:
In formula (1),WithRepresent respectively system y gross investment construction cost, fixed run
Cost, variable operation cost and interregional interconnection construction cost,Represent carbon transaction cost;Y represents all years to be planned
Part, r represents inflation rate;
1-3-1) shown in total Installed capital cost such as formula (2):
Model pairUsing etc. year value decompose form, the cost of investment of all kinds of power plant is divided into power plant by discount rate
In Years Of Service, and ignore residual value of the power supply after project period.In addition, the unit of all kinds of power plant is held in model acquiescence FX
Amount cost of investment is definite value, therefore in formulaThere is no footnote b;Wherein A represents the set of all subregion in model, and U represents to treat
Plan the set of power type;
1-3-2) fixed operating cost
In formula, Pu,y,aWithRepresent that the active volume of u class power supplys and unit capacity fix operation in a of y regions
Cost.Active volume Pu,y,aIt is expressed from the next:
In formula, Pu,y,a,bThe extension capacity of u class power supplys powered in a of y regions for region b loads is represented, is had:
In formula (5), Pu,0,a,bThe appearance for the u class power supplys powered in region a at the beginning of representing project period for region b loads
Amount,Represent retired capacity of such power supply in i-th of planning year;
1-3-3) interregional interconnection construction cost is represented by formula (6):
In formula,WithRepresent that y is built from region a to the newly-built capacity of region b power transmission passages and unit capacity
This is set as, B represents the set of optional power transmission passage;
1-3-4) variable operation cost
Variable operation cost is system cost of electricity-generating, can be represented by formula (7):
In formula,The unit cost of electricity-generating of u class power supplys in a of y regions is represented, D is all typical cases of power system
Day set under the method for operation, TD is period set of the power system in typical day D, and Δ T is period lasts duration;
1-3-5) carbon emission transaction cost
It is embodied as:
Wherein Ey,aThe electricity consumption carbon emission amount of y in a of region is represented,Y's is arranged with electrical carbon in expression region a
Put quota, πy.aRepresent the price of y unit carbon emission quotas in a of region;Ey,aIt can be expressed from the next:
In formula, eu,bRepresent the carbon potential of u class units in sending end region b;
The constraints of low-carbon power generating facilities and power grids Optimal Planning Model 1-4) is set up, is specifically included:
1-4-1) power supply annual electricity generating capacity is constrained
This is constrained to all kinds of power supplys, and the constraint of its annual electricity generating capacity is equivalent to its year gas-to electricity hourage and must not exceed its year
Highest can utilize hourage Tu,y,a max, can not also be less than certain minimum utilization hour T in yearu,y,a min, therefore have:
In formula, NtAnd NdThe quantity of quantity of the typical period of time t in typical day d and typical day d in y is represented respectively;
1-4-2) subregion power supply and demand is constrained
The total capacity for being constrained to the annual all kinds of power supplys of power system should keep balancing with this year system load demand, together
When, being used to balance the power supply capacity of local load in each region should be not less than with other regions to the feeding capacity sum in the region
The peak load in the region, the latter is the former adequate condition, i.e.,:
In formula, Pu,y,b,aThe existing capacity of u class power supplys powered in the first term area b of planning for region a loads is represented,
Py,a maxFor y regions a maximum predicted load;
1-4-3) region power supply maximum can development capacity constraint
The consumed primary energy consumption amount sum of thermoelectricity generating for being constrained to each region must not exceed primary energy most
The big amount of can be supplied to, to region a, has:
In formula, fu,y,aRepresent the energy consumption of region a moderate heat individual item generated energy, UnThen represent all with primary energy n to make
Gather for the power plant of energy input, Fy,n,aRepresent Maximum Supply Quantitys of such primary energy n in a of region;
For the constraint of wind-powered electricity generation, water power, should limit its maximum year by year can development capacity, the limitation depends on the year of resource
Can development capacity, the or manufacturing maximum production of power-supply device have:
In formula, Pu,y,a maxRepresent that u classes power supply can development capacity, P by y maximum in a of regionu,0,aThen represent
In the in-service installed capacity of such initial power supply of planning;
1-4-4) interregional electricity contract constraint
The electricity that the constraint can behave as the u class power supplys powered in a of region for region b loads is no less than Contract Energy;
I.e.:
In formula,For the Contract Energy of the u class power supplys from region a to region b;
1-4-5) interregional interconnection capacity-constrained
This is constrained to the demand that interconnection channel capacity must assure that interregional exchange of electric power, i.e.,:
P in formula(a,b),yFor the active volume of the point-to-point power transmission passages of y regions a in equivalent network to region b, specifically
For:
P in formula(a,b),0For the existing capacity of the point-to-point power transmission passages of region a to region b in planning equivalent network at initial stage;
In addition, the extension of interconnection capacity also needs to consider that two interregional interconnection capacity have the upper limit, i.e.,:
P(a,b),y≤Pab max a,b∈A (29)
For the interconnection of point-to-point transmission in real network, interconnection capacity-constrained between above-mentioned zone can be directly applied;
For the interconnection in virtual equivalent network, formula (17) is constrained then equivalent to profile constraints, (a, b) represent region a and region b it
Between virtual interconnection set;
1-4-6) interregional interconnection Constraint
Hour can be utilized by being equivalent to its annual utilization hours to the year Constraint of each interconnection and must not exceed its year highest
Number T(a,b),y max, i.e.,:
1-4-7) typical day operation peak regulation and Reserve Constraint
This is constrained under given load curve of typical day, and group need that have been switched in each region meet the maximum of one's respective area and born
Lotus and peak regulation and stand-by requirement.On the premise of in a few days Unit Commitment is not considered, typical day operation constraint is embodied as:
In formula, UdThe set of the start unit in typical day d is represented, D is typical day set, αa,uRepresent type in a of region
The accounting of rated capacity is accounted for for u unit minimum output;Pb,y,d maxAnd Pb,y,d minRepresent that typical day d is most in the b of y regions
Big and minimum load,WithRRepresent positive percentage reserve and negative percentage reserve;
1-4-8) typical day unit output constraint
This was constrained under given typical day, and all kinds of units of power system are exerted oneself no more than the available appearance of the type unit
Amount:
Meanwhile, exerting oneself for all kinds of units of system can not be less than the minimum output of the type unit:
In formula,Represent y, the active volume of the u classes unit in a regions in typical day d accounts for Pu,y,a,bRatio
Value.Represent unit capacity minimum output ratio of the u classes unit in a regions in typical day d;
1-4-9) power generating facilities and power grids investment construction is constrained
This is constrained to the limitation of the power supply gone into operation every year each department and interzone interconnection capacity, should meet such as formula (22) institute
Show:
1-4-10) carbon emission transaction constraint
The constraint is used as carbon emission measurement criteria using electricity consumption rather than generating;What power industry can purchase in one's respective area matches somebody with somebody
There is the upper limit in volume, total electricity consumption carbon emission has the upper limit, i.e.,:
E in formulamax y,aThe electricity consumption carbon emission maximum permissible value of y after consideration carbon transaction in a of region is represented, can be by each
Regional government independently formulates, when there is the carbon emission control targe clearly towards power industry in the government department of system region
When, Emaxy,aDirectly give;
1-4-11) newly-built power supply nonnegativity restriction
I.e. all kinds of power supply operation total amounts in region need to meet nonnegativity restriction, as shown in formula (24):
1-4-12) electricity and interconnection decision variable nonnegativity restriction, as shown in formula (25):
1-4-13) constraint of contact line justification, is stated as shown in formula (26):
In formula,Represent the set in region pair exchanged in the absence of the energy;
2) to the solution of low-carbon power generating facilities and power grids Optimized model
The low-carbon power generating facilities and power grids Optimized model is expressed as linear model:
min cT·x
s.t.A·x≤b
(29)
AeqX=beq
x≥0
Wherein A and b characterizes the inequality constraints in optimization problem, AeqAnd beqCharacterize the equality constraint in optimization problem, x
For decision variable, cTFor system cost term coefficient;
2-1) by step 1) each class cost is represented as in the object functionForm, by formula (2)~(9)
Construct ck.Order:
K represents the number of cost in object function in formula, thus the c in constructive formula (27);
According to step 1-4) in all kinds of constraintss, each class equality constraint and inequality constraints are represented as Ai·x
≤biOr AeqjX=beqjForm;By step 1-4) construction Ai、biAnd Aeqi、beqi:Order:
In formula (29), I and J represent the number of inequality constraints and equality constraint respectively, thus by the pact in Optimized model
Beam condition is converted to available matrix computations form;
Solver is worked out by MATLAB, calls CPLEX12.5 to carry out the low-carbon power generating facilities and power grids Optimized model
Solve, try to achieve the power generating facilities and power grids optimization planning of the low-carbon of all kinds of power supplys and transmission line of electricity in power system so that power system
There is optimal economic benefit under low-carbon development model.
The features of the present invention:
Accurate consideration carbon emission constraint of the invention is transferred to electricity consumption link from generating link, that is, introduces the constraint of electricity consumption carbon emission
When low-carbon power generating facilities and power grids planning.Model will attempt carbon of the application based on consumption side in the investment problem that conventional electric power is planned
Emission measurement angle is carried out to key elements such as carbon emission constraint, the transaction of inter-trade carbon emission, the interregional electricity contracts of electricity consumption link
Description, and consider that carbon price differential between different regions is different, system operation peak regulation Reserve Constraint, form complete planing method.
It is multiple it is interregional there is electricity exchange, when there is the constraint of respective electricity consumption carbon emission, it is contemplated that interregional carbon emission transfer
Source is generally the interregional fossil capacity supported each other, and this method describes interregional carbon row by sending unit and direct sending circuit outside
Put sharing for responsibility.And then the concept by setting up local capacity He sending capacity outside, the electric power networks of multi area interconnection are modeled
Power generating facilities and power grids optimization planning is carried out, allows whole system to carry out the extension of power generating facilities and power grids in the way of the lowest cost.
Beneficial effects of the present invention:
The inventive method accurately considers that carbon emission constraint is transferred to electricity consumption link from generating link, that is, introduces electricity consumption carbon emission
Low-carbon power generating facilities and power grids plan model during constraint, will attempt application based on consumption side in the investment problem that conventional electric power is planned
Carbon emission measuring angle to key elements such as carbon emission constraint, the transaction of inter-trade carbon emission, the interregional electricity contracts of electricity consumption link
It is described, and considers that carbon price differential between different regions is different, system operation peak regulation Reserve Constraint, forms complete planning.
It is multiple it is interregional there is electricity exchange, when there is the constraint of respective electricity consumption carbon emission, allow whole power system with totle drilling cost most
Low mode carries out the extension of power generating facilities and power grids.
The inventive method can provide practical advice for the optimization planning of power system, introduce the basis of electricity consumption side carbon emission
On, the development of power generating facilities and power grids is planned by objective optimization of the lowest cost, is that the low-carbonization operation of system lays a good foundation.
Embodiment
A kind of method of the power generating facilities and power grids optimization planning of low-carbon proposed by the present invention further illustrate in conjunction with the embodiments as
Under:
The method and specific implementation of the power generating facilities and power grids optimization planning of the low-carbon of the present invention, comprise the following steps:
1) the low-carbon power generating facilities and power grids Optimized model being made up of object function and constraints is built
Specifically include:
All kinds of power supply capacitys in each region 1-1) are divided into two major classes, a class is local capacity, for balancing this
The capacity of ground load;Another kind of is to send capacity outside, the capacity for supplying other region loads;
Interregional point-to-point direct sending circuit is set up in equivalent network, is arranged with the flowing of clearly interregional electricity and carbon
The transfer put, is easy to the metering of electricity consumption carbon emission;
1-2) set up the decision variable of low-carbon power generating facilities and power grids Optimized model
The model is using year as decision-making time unit, and the decision variable in model mainly includes following three class:The first kind isRepresent that a and b represents area in the extension capacity of u class power supplys powered in a of y regions for region b loads, footnote
Domain, N represents power extension capacity;When a and b represent the same area, such asDuring for local capacity, show that such unit is increased newly
Capacity is used for the local loads of equilibrium region a;When a and b represent different zones, show that such unit increases capacity newly and is used for equilibrium area
Domain b load, to send capacity outside.Equations of The Second Kind isY regions a is represented to the newly-increased capacity of power transmission passage between the b of region,
I represents power transmission passage and increases capacity newly.3rd class isFor y, in typical day d under period t, born in a of region for region b
U class power supplys that lotus is powered are exerted oneself.
1-3) set up low-carbon power generating facilities and power grids Optimized model object function:
In formula (1),WithRepresent respectively system y gross investment construction cost, fixed run
Section 5 in cost, variable operation cost and interregional interconnection construction cost, object functionRepresent carbon transaction cost;Its
InWithDescribed by investment problem,By running subproblem description;Y represents all times to be planned, and r represents logical
Goods expansion rate;The expression of each cost can be calculated by formula (2), (3), (6), (7) and (8) and obtained in object function;
1-3-1) total Installed capital cost
Model pairUsing etc. year value decompose form, the cost of investment of all kinds of power plant is divided into power plant by discount rate
Years Of Service in, and ignore residual value of the power supply after project period.In addition, model gives tacit consent to the unit of all kinds of power plant in FX
Capacity cost of investment is definite value, therefore in formulaThere is footnote b again;Wherein A represents the set of all subregion in model, and U represents to treat
Plan the set of power type;
1-3-2) fixed operating cost
In formula, Pu,y,aWithRepresent that the active volume of u class power supplys and unit capacity fix operation in a of y regions
Cost.Active volume Pu,y,aIt can be expressed from the next:
In formula, Pu,y,a,bThe extension capacity of u class power supplys powered in a of y regions for region b loads is represented, is had:
In formula, Pu,0,a,bThe existing capacity for the u class power supplys powered in region a at the beginning of representing project period for region b loads,Represent retired capacity of such power supply in i-th of planning year;
1-3-3) interregional interconnection construction cost
Interregional interconnection construction cost can be expressed from the next:
In formula,WithRepresent y from region a to the newly-built capacity of region b power transmission passages and unit capacity
Construction cost, B represents the set of optional power transmission passage, and the direct sending that power transmission passage herein refers in improved equivalent network leads to
Road;
1-3-4) variable operation cost
Total variable operation cost is system cost of electricity-generating, can be represented by formula (7):
In formula,The unit cost of electricity-generating of u class power supplys in a of y regions is represented, D is all typical case's operations of system
Day set under mode, TD is period set of the system in typical day D, and Δ T is period lasts duration;
1-3-5) carbon emission transaction cost
Carbon emission transaction cost is introduced to be meant that:When system, each the electricity consumption total release of subdomain is higher than and always allowed
During quota of discharge, system needs to buy extra quota of discharge to cause cost;And when system total release is matched somebody with somebody less than discharge
During volume, remaining quota can be sold and earned a profit.Specifically it is represented by:
Wherein Ey,aThe electricity consumption carbon emission amount of y in a of region is represented,Y's is arranged with electrical carbon in expression region a
Put quota, πy.aRepresent the price of y unit carbon emission quotas in a of region.This model thinks that carbon emission quota gives to be outside
Value.EY, aIt can be expressed from the next:
In formula, eu,bRepresent the carbon potential of u class units in sending end region b;
1-4) set up the constraints of low-carbon power generating facilities and power grids Optimal Planning Model
1-4-1) power supply annual electricity generating capacity is constrained
To all kinds of power supplys, the constraint of its annual electricity generating capacity be equivalent to its year gas-to electricity hourage must not exceed its year highest can profit
Use hourage Tu,y,a max(highest is mainly derived from primary energy constraint, unit maintenance etc. using hourage), can not also be less than one
Minimum utilization hour T in fixed yearu,y,a min(minimum utilization hourage forces from Policy Conditions, water power and exerts oneself or maintain power plant
Operation needs), therefore have:
In formula, NtAnd NdThe quantity of quantity of the typical period of time t in typical day d and typical day d in y is represented respectively;
1-4-2) subregion power supply and demand is constrained
The total capacity of the annual all kinds of power supplys of system should keep balancing with this year system load demand, meanwhile, used in each region
The maximum in the region should be not less than to the feeding capacity sum in the region in the power supply capacity and other regions for balancing local load
Load, the latter is the former adequate condition, i.e.,:
In formula, Pu,y,b,aThe existing capacity of u class power supplys powered in the first term area b of planning for region a loads is represented,
Py,a maxFor y regions a maximum predicted load;
1-4-3) region power supply maximum can development capacity constraint
For fired power generating unit, power supply, which generates electricity, to be needed to consume primary energy, predominantly coal, natural gas, oil etc.;Obviously,
The primary energy consumption amount sum that the thermoelectricity generating in each region is consumed must not exceed the primary energy maximum amount of can be supplied to, and be presented as
The maximum Constraint of system thermoelectricity year, to region a, has:
In formula, fu,y,aRepresent the energy consumption of region a moderate heat individual item generated energy, UnThen represent all with primary energy n to make
Gather for the power plant of energy input, Fy,n,aRepresent that Maximum Supply Quantitys of such primary energy n in a of region be not (straight in this model
Connect consideration primary energy transport constraint and cost);
The power supply of fuel need not be expended for wind-powered electricity generation, water power, it is not required to the constraint in accordance with above formula, and it should be limited year by year
Maximum can development capacity, the limitation depend on resource year can development capacity, or the manufacturing maximum production of power-supply device, i.e.,
Have:
In formula, Pu,y,a maxRepresent that u classes power supply can development capacity, P by y maximum in a of regionu,0,aThen represent
In the in-service installed capacity of such initial power supply of planning;
1-4-4) interregional electricity contract constraint
It is different it is interregional can freely sign bilateral electricity contract, exchanged as interregional electricity or even carbon emission share
Foundation.By taking the contract that region a and region b is present as an example, the constraint can behave as the u classes powered in a of region for region b loads
The electricity of power supply is no less than Contract Energy;I.e.:
In formula,For the Contract Energy of the u class power supplys from region a to region b;
1-4-5) interregional interconnection capacity-constrained
Interconnection channel capacity must assure that the demand of interregional exchange of electric power, i.e.,:
P in formula(a,b),yFor the active volume of the point-to-point power transmission passages of y regions a in equivalent network to region b, specifically
For:
P in formula(a,b),0For the existing capacity of the point-to-point power transmission passages of region a to region b in planning equivalent network at initial stage;
In addition, the extension of interconnection capacity also needs to consider the limitation of the factors such as geographical and execution conditions, two interregional contacts
There is the upper limit in line capacity, i.e.,:
P(a,b),y≤Pab max a,b∈A (46)
For the interconnection of point-to-point transmission in real network, interconnection capacity-constrained between above-mentioned zone can be directly applied;
For the interconnection in virtual equivalent network, formula (17) is constrained then equivalent to profile constraints, (a, b) represent region a and region b it
Between virtual interconnection set;
1-4-6) interregional interconnection Constraint
To each interconnection, its year Constraint, which is equivalent to its annual utilization hours and must not exceed its year highest, can utilize hour
Number T(a,b),y max, i.e.,:
1-4-7) typical day operation peak regulation and Reserve Constraint
This constraint facing to manufacture is simulated, under the load curve that typical day gives, and the group that has been switched in each region needs satisfaction
The peak load and peak regulation of one's respective area and stand-by requirement.On the premise of in a few days Unit Commitment is not considered, typical day operation is about
Beam is embodied as:
In formula, UdThe set of the start unit in typical day d is represented, D is typical day set, αa,uRepresent type in a of region
The accounting of rated capacity is accounted for for u unit minimum output;Pb,y,d maxAnd Pb,y,d minRepresent that typical day d is most in the b of y regions
Big and minimum load,WithRRepresent positive percentage reserve and negative percentage reserve;
(for the new energy unit (based on wind-powered electricity generation) in each region, its anti-peak-shaving capability need to be considered, it is just standby in analysis system
During with constraint, by the new energy unit start capacity (P of correspondence Wind turbinesu,y,a,b) it is considered as zero, and born in analysis system standby
During constraint, by the corresponding P that exerts oneself of wind-powered electricity generationu,y,a,bIt is taken as EIAJ under certain confidence level.It is processed as may insure system
Peak regulation meets the requirement of actual motion with Reserve Constraint;)
1-4-8) typical day unit output constraint
This constraint facing to manufacture is simulated, under given typical day, no more than the type machine of exerting oneself of all kinds of units of system
The active volume of group:
Meanwhile, exerting oneself for all kinds of units of system can not be less than the minimum output of the type unit:
In formula,Represent y, the active volume of the u classes unit in a regions in typical day d accounts for Pu,y,a,bRatio
Value.Represent unit capacity minimum output ratio of the u classes unit in a regions in typical day d;
1-4-9) power generating facilities and power grids investment construction is constrained
The stationarity of the expansion scheme of construction ability and each region power supply and power network in view of power generating facilities and power grids, table
The power supply and interzone interconnection capacity now gone into operation every year for each department should meet a definite limitation, as shown in formula (22):
1-4-10) carbon emission transaction constraint
In view of the industrial structure of different zones, economic development planning in system and to National Macroscopic low carbon development target
Response, each region to the carbon emission amount of itself by set constrain.Before the formation of trans-regional carbon emission mechanism of exchange, each region
Power network can seek the carbon emission transaction cooperation between other industry in region;In order to promote each region can mode and production to itself using
The optimization of industry structure, rationally recognizes itself electricity consumption carbon emission, and the think of of the carbon emission metering based on consumption side is introduced into this constraint
Think, carbon emission measurement criteria is used as using electricity consumption rather than generating;In view of carbon emission quota in each region and the limitation of transaction, electricity
There is the upper limit in the quota that Lixing industry can purchase in one's respective area, be presented as that total electricity consumption carbon emission constraint has the upper limit, i.e.,:
E in formulamax y,aThe electricity consumption carbon emission maximum permissible value of y after consideration carbon transaction in a of region is represented, can be by each
Regional government independently formulates, and is used as outside known conditions.When there is clear and definite face in the government department (or country) of system region
To power industry carbon emission control targe when, Emaxy,aCan directly it give;
1-4-11) newly-built power supply nonnegativity restriction
In power source planning problem, generally there is the nonnegativity restriction of newly-built power supply capacity, set forth herein multizone
In power source planning problem, it is allowed to decision variableThere is negative value, but all kinds of power supply operation total amounts in all regions need to meet non-negative
Property constraint, i.e.,:
1-4-12) electricity and interconnection decision variable nonnegativity restriction
In addition to newly-built power constraints, remaining decision variable meets nonnegativity restriction, as shown in formula (25):
1-4-13) contact line justification constraint
According to ENERGY PLANNING strategy or system actual conditions, not any two each interregional can enter in multizone system
Row electrical energy transportation, for the region pair exchanged in the absence of the energy, its interconnection capacity need to be set to zero, rational to obtain
Optimum results.With reference to decision variable nonnegativity restriction, contact line justification constraint can be stated as shown in formula (26):
In formula,Represent the set in region pair exchanged in the absence of the energy;
The low-carbon power generating facilities and power grids Optimized model is linear model;
2) to the solution of low-carbon power generating facilities and power grids Optimized model
The general type of linear model is expressed as:
min cT·x
s.t.A·x≤b
(56)
AeqX=beq
x≥0
Wherein A and b characterizes the inequality constraints in optimization problem, AeqAnd beqCharacterize the equality constraint in optimization problem, x
For decision variable, cTFor system cost term coefficient;
2-1) by step 1) each class cost is represented as in the object functionForm, by formula (2)~(9)
Construct ck.Order:
K represents the number of cost in object function in formula, thus the c in constructive formula (27);
According to step 1-4) in all kinds of constraintss, each class equality constraint and inequality constraints are represented as Ai·x
≤biOr AeqjX=beqjForm;By step 1-4) construction Ai、biAnd Aeqi、beqi:Order:
In formula (29), I and J represent the number of inequality constraints and equality constraint respectively, thus can be by Optimized model
Constraints is converted to available matrix computations form;
Solver is worked out by MATLAB, calls CPLEX12.5 to solve above-mentioned Optimized model, can be in the hope of electricity
The power generating facilities and power grids optimization planning of the low-carbon of all kinds of power supplys and transmission line of electricity in Force system so that power system is in low carbon development mould
There is optimal economic benefit under formula.
The planning implementation example that the following power generating facilities and power grids of the province of south China five extend prospect is described as follows with the above method:
The project period of the present embodiment be 2015 to the year two thousand twenty, to show the application effect of this method.South electric network is pressed
Each province is considered as five node systems, does not repartition its internal different zones to each node, only distinguishes its internal inhomogeneity
The unit of type.It is main to consider water power, thermoelectricity, four kinds of power supply types of nuclear power and new energy.According to《South electric network " 12 " develops
Program results collect》" 12 " and medium-term and long-term electric power development project study with each province under south electric network region, required base
Plinth data are as follows:
(1) electricity needs
Each peak loads inside the province in 2015 and prediction " 13 " growth rate are as shown in the table:
1 south electric network of table 2015 and " 13 " load prediction
(2) power parameter
The installed capacity of each province's different type power supply in 2015 refers to following table:
Each region installed capacity in 2015 of 2 south electric network of table and composition (ten thousand kW)
For all kinds of units in different provinces using hourage, maximum can the parameter such as development capacity, there is difference different regions
It is different, obtained according to south electric network entirety and " 12 " of each province and medium-term and long-term electric power industry development planning file.Construction of Unit
Cost, operating cost etc. use canonical parameter.
(3) equivalent network parameter
Set up southern five provinces and regions simplification equivalent network and acquire power supply and interconnection that model is available for programmed decision-making.
Each interregional interconnection situation is as shown in table 3, table 4 within 2015:
3 south electric network of table each provinces and regions interconnection capacity (ten thousand kW) in 2015
4 south electric network of table each provinces and regions interconnection length (km) in 2015
To simplify the analysis, all interconnectors transprovincially planned in the present embodiment are 500kV AC or DC power transmission lines
Road, range estimation result reference table 4.
Each provinces and regions interconnection unit capacity cost of the south electric network of table 5 (member/kW)
(4) analogue data is produced
One-year age is set as 8760 hours, the present embodiment takes two typical days, is respectively peak day then
With minimum load day, every kind of typical day continue 4380 hours.Each typical case's day includes two periods, when being respectively the peak of this day
Section and low-valley interval, each period lasts 12 hours.To 2015 to the planning years of the year two thousand twenty six, 12 typical days totally 24 allusion quotations
The associative simulation of type period.
(5) electricity contract data
According to south electric network actual conditions, set Yunnan Province and export electricity as water power, sent in the total electricity of Guizhou Province's output
Electricity toward Guangdong Province is thermoelectricity.If it is as follows that what each year Guizhou project period and Yunnan need to be met always sends electricity contract outside:
The province of 6 Yunnan-Guizhou of table two sends contract demand (hundred million kWh) outside
Remaining Transmission Corridor is constrained without electricity contract.
(6) carbon emission data
Carbon emission quota is based on compared to 2005 per GDP carbon intensities of China's the year two thousand twenty and declines 40%~45% system
It is fixed.Compared to 2005 2015 per GDP carbon intensities can be limited and decline 30%.Decline 32% respectively within 2016 to 2020
To 42%.On the other hand, to determine carbon consumption constraint, the pass between clear and definite power industry carbon emission and whole society's carbon emission is also needed
System.Guarded in this example and assume that 2015 ratios for accounting for whole society's total carbon emission to power industry carbon emission between the year two thousand twenty are maintained at
50%.Thus, the electricity consumption carbon emission quota for showing that each province is linked up with for 2015 to the year two thousand twenty with GDP phases can be estimated, as shown in table 7,
Unit is ten thousand tons.Concurrently set each province's carbon transaction total amount must not exceed the whole province's carbon emission quota total amount 10%.
Each time electricity consumption carbon emission quota in each provinces and regions of table 7 (ten thousand tCO2)
In view of each province carbon market difference, it is assumed that price is not changed over, setting each province's carbon quota price is as follows:
Each time electricity consumption carbon emission quota price in each provinces and regions of table 8
Above-mentioned model is programmed with MATLAB, calls CPLEX12.5 to be solved, result is obtained as follows:
(1) power source planning
Power source planning is as follows year by year for each province after arrangement:
Each provinces and regions each time of table 9 increases hydropower installed capacity (ten thousand kW) newly
Each provinces and regions each time of table 10 increases thermoelectricity installed capacity (ten thousand kW) newly
Each provinces and regions each time of table 11 increases nuclear power installed capacity (ten thousand kW) newly
Each provinces and regions each time of table 12 increases fuel gas generation installed capacity (ten thousand kW) newly
Each provinces and regions each time of table 13 increases generation of electricity by new energy installed capacity (ten thousand kW) newly
On the whole, south electric network seeks big water power under low-carbon constraint and first developing for big nuclear power is imitated with obtaining low-carbon
Benefit, in the province that hydroelectric resources is not enough or regulating power is not enough, Gas Generator Set turns into preferred.And the new energy based on wind-powered electricity generation
On the premise of small be also towards low-carbon planning preferred pair as.
(2) power network expansion scheme
Inter-provincial each interconnection expansion scheme is as follows after arrangement:
Each provinces and regions each time of table 14 increases interconnection capacity (ten thousand kW) newly
Finally it should be noted that:The above embodiments are merely illustrative of the technical scheme of the present invention and are not intended to be limiting thereof, to the greatest extent
The present invention is described with reference to above-described embodiment for pipe, those of ordinary skills in the art should understand that:Still can be right
The embodiment of the present invention is modified or equivalent substitution, and any modification without departing from spirit and scope of the invention or
Person's equivalent substitution, it all should cover among scope of the presently claimed invention.
Claims (1)
1. the method for the power generating facilities and power grids optimization planning of a kind of low-carbon, it is characterised in that this method comprises the following steps:
1) the low-carbon power generating facilities and power grids Optimized model being made up of object function and constraints is built, is specifically included:
All kinds of power supply capacitys in each region 1-1) are divided into two major classes, a class is local capacity, it is local negative for balancing
The capacity of lotus;Another kind of is to send capacity outside, the capacity for supplying other region loads;
Interregional point-to-point direct sending circuit is set up in equivalent network, with the flowing of clearly interregional electricity and carbon emission
Transfer, is easy to the metering of electricity consumption carbon emission;
1-2) set up the decision variable of low-carbon power generating facilities and power grids Optimized model
The model is using year as decision-making time unit, and the decision variable in model mainly includes following three class:The first kind is
Represent that a and b represents region in the extension capacity of u class power supplys powered in a of y regions for region b loads, footnote;Work as a
The same area is represented with b, such asDuring for local capacity, show that such power supply increases capacity newly and locally born for equilibrium region a
Lotus;When a and b represent different zones, show that such power supply increases capacity newly and is used for equilibrium region b load, to send capacity outside;The
Two classes areRepresent y regions a to the newly-increased capacity of power transmission passage between the b of region;3rd class isFor y,
It is exerting oneself for the u class power supplys that region b loads are powered in typical day d under period t, in a of region;
1-3) set up low-carbon power generating facilities and power grids Optimized model object function:
<mrow>
<mi>min</mi>
<mi> </mi>
<mi>C</mi>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>y</mi>
<mo>&Element;</mo>
<mi>Y</mi>
</mrow>
</munder>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mi>y</mi>
</mrow>
</msup>
<mo>&CenterDot;</mo>
<mo>&lsqb;</mo>
<msubsup>
<mi>C</mi>
<mi>y</mi>
<mi>I</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>C</mi>
<mi>y</mi>
<mi>O</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>C</mi>
<mi>y</mi>
<mi>G</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>C</mi>
<mi>y</mi>
<mi>T</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>C</mi>
<mi>y</mi>
<mi>E</mi>
</msubsup>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula (1),WithRespectively represent system y gross investment construction cost, fixation operating cost,
Variable operation cost and interregional interconnection construction cost,Represent carbon transaction cost;Y represents all times to be planned, r generations
Table inflation rate;
1-3-1) shown in total Installed capital cost such as formula (2):
<mrow>
<msubsup>
<mi>C</mi>
<mi>y</mi>
<mi>I</mi>
</msubsup>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>a</mi>
<mo>&Element;</mo>
<mi>A</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>b</mi>
<mo>&Element;</mo>
<mi>A</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>u</mi>
<mo>&Element;</mo>
<mi>U</mi>
</mrow>
</munder>
<msubsup>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
<mi>N</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>c</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
<mi>I</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
Model pairUsing etc. year value decompose form, the cost of investment of all kinds of power plant is divided to the military service of power plant by discount rate
In the time limit, and ignore residual value of the power supply after project period;In addition, the unit capacity of all kinds of power plant is thrown in model acquiescence FX
Money cost is definite value, therefore in formulaThere is no footnote b;Wherein A represents the set of all subregion in model, and U represents to be planned
The set of power type;
1-3-2) fixed operating cost
<mrow>
<msubsup>
<mi>C</mi>
<mi>y</mi>
<mi>O</mi>
</msubsup>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>a</mi>
<mo>&Element;</mo>
<mi>A</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>u</mi>
<mo>&Element;</mo>
<mi>U</mi>
</mrow>
</munder>
<msub>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>c</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
<mi>O</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, Pu,y,aWithRepresent that the active volume of u class power supplys and unit capacity fix operating cost in a of y regions;
Active volume Pu,y,aIt is expressed from the next:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>b</mi>
<mo>&Element;</mo>
<mi>A</mi>
</mrow>
</munder>
<msub>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, Pu,y,a,bThe extension capacity of u class power supplys powered in a of y regions for region b loads is represented, is had:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mn>0</mn>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>y</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>i</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
<mi>I</mi>
</msubsup>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>y</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>i</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
<mi>D</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula (5), Pu,0,a,bThe existing capacity for the u class power supplys powered in region a at the beginning of representing project period for region b loads,Represent retired capacity of such power supply in i-th of planning year;
1-3-3) interregional interconnection construction cost is represented by formula (6):
<mrow>
<msubsup>
<mi>C</mi>
<mi>y</mi>
<mi>T</mi>
</msubsup>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
<mo>&Element;</mo>
<mi>B</mi>
</mrow>
</munder>
<msubsup>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
<mo>,</mo>
<mi>y</mi>
</mrow>
<mi>I</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>c</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<mi>I</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula,WithRepresent that y is built into from region a to the newly-built capacity of region b power transmission passages and unit capacity
This, B represents the set of optional power transmission passage;
1-3-4) variable operation cost
Variable operation cost is system cost of electricity-generating, can be represented by formula (7):
<mrow>
<msubsup>
<mi>C</mi>
<mi>y</mi>
<mi>G</mi>
</msubsup>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>d</mi>
<mo>&Element;</mo>
<mi>D</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>&Element;</mo>
<msub>
<mi>T</mi>
<mi>D</mi>
</msub>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>a</mi>
<mo>&Element;</mo>
<mi>A</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>b</mi>
<mo>&Element;</mo>
<mi>A</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>u</mi>
<mo>&Element;</mo>
<mi>U</mi>
</mrow>
</munder>
<msubsup>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
<mrow>
<mi>d</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>c</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
<mi>G</mi>
</msubsup>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>T</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula,The unit cost of electricity-generating of u class power supplys in a of y regions is represented, D is all typical operation sides of power system
Day set under formula, TD is period set of the power system in typical day D, and Δ T is period lasts duration;
1-3-5) carbon emission transaction cost
It is embodied as:
<mrow>
<msubsup>
<mi>C</mi>
<mi>y</mi>
<mi>E</mi>
</msubsup>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>a</mi>
<mo>&Element;</mo>
<mi>A</mi>
</mrow>
</munder>
<msub>
<mi>&pi;</mi>
<mrow>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>E</mi>
<mi>a</mi>
<mrow>
<mi>c</mi>
<mi>a</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein Ey,aThe electricity consumption carbon emission amount of y in a of region is represented,Represent that the electricity consumption carbon emission of y in a of region is matched somebody with somebody
Volume, πy.aRepresent the price of y unit carbon emission quotas in a of region;Ey,aIt can be expressed from the next:
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>d</mi>
<mo>&Element;</mo>
<mi>D</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>&Element;</mo>
<msub>
<mi>T</mi>
<mi>D</mi>
</msub>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>b</mi>
<mo>&Element;</mo>
<mi>A</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>u</mi>
<mo>&Element;</mo>
<mi>U</mi>
</mrow>
</munder>
<msubsup>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>b</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
<mrow>
<mi>d</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<msub>
<mi>e</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>T</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, eu,bRepresent the carbon potential of u class units in sending end region b;
The constraints of low-carbon power generating facilities and power grids Optimal Planning Model 1-4) is set up, is specifically included:
1-4-1) power supply annual electricity generating capacity is constrained
This is constrained to all kinds of power supplys, and the constraint of its annual electricity generating capacity is equivalent to its year gas-to electricity hourage and must not exceed its year highest
Using hourage Tu,y,a max, can not also be less than certain minimum utilization hour T in yearu,y,a min, therefore have:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>&le;</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>d</mi>
<mo>&Element;</mo>
<mi>D</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>&Element;</mo>
<msub>
<mi>T</mi>
<mi>D</mi>
</msub>
</mrow>
</munder>
<msub>
<mi>N</mi>
<mi>t</mi>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>N</mi>
<mi>d</mi>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
<mrow>
<mi>d</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>T</mi>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, NtAnd NdThe quantity of quantity of the typical period of time t in typical day d and typical day d in y is represented respectively;
1-4-2) subregion power supply and demand is constrained
The total capacity for being constrained to the annual all kinds of power supplys of power system should keep balancing with this year system load demand, meanwhile, respectively
The area should be not less than with other regions to the feeding capacity sum in the region by being used to balance the power supply capacity of local load in region
The peak load in domain, the latter is the former adequate condition, i.e.,:
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>u</mi>
<mo>&Element;</mo>
<mi>U</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>b</mi>
<mo>&Element;</mo>
<mi>A</mi>
</mrow>
</munder>
<msub>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>b</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>&GreaterEqual;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, Pu,y,b,aRepresent the existing capacity of u class power supplys powered in the first term area b of planning for region a loads, Py,a max
For y regions a maximum predicted load;
1-4-3) region power supply maximum can development capacity constraint
The consumed primary energy consumption amount sum of thermoelectricity generating for being constrained to each region must not exceed primary energy maximum can
Quantity delivered, to region a, has:
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>u</mi>
<mo>&Element;</mo>
<msub>
<mi>U</mi>
<mi>n</mi>
</msub>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>b</mi>
<mo>&Element;</mo>
<mi>A</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>d</mi>
<mo>&Element;</mo>
<mi>D</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>&Element;</mo>
<msub>
<mi>T</mi>
<mi>D</mi>
</msub>
</mrow>
</munder>
<msubsup>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
<mrow>
<mi>d</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>T</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>f</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>F</mi>
<mrow>
<mi>y</mi>
<mo>,</mo>
<mi>n</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, fu,y,aRepresent the energy consumption of region a moderate heat individual item generated energy, UnThen represent all and energy is used as using primary energy n
Power plant's set of amount input, Fy,n,aRepresent Maximum Supply Quantitys of such primary energy n in a of region;
For the constraint of wind-powered electricity generation, water power, should limit its maximum year by year can development capacity, the limitation depends on the Nian Kekai of resource
Capacity, or the manufacturing maximum production of power-supply device are sent out, that is, is had:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mn>0</mn>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>y</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>i</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
<mi>N</mi>
</msubsup>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>,</mo>
<mi>u</mi>
<mo>&Element;</mo>
<mi>U</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, Pu,y,a maxRepresent that u classes power supply can development capacity, P by y maximum in a of regionu,0,aThen represent in rule
Draw the in-service installed capacity of such initial power supply;
1-4-4) interregional electricity contract constraint
The electricity that the constraint can behave as the u class power supplys powered in a of region for region b loads is no less than Contract Energy;I.e.:
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>d</mi>
<mo>&Element;</mo>
<mi>D</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>&Element;</mo>
<msub>
<mi>T</mi>
<mi>D</mi>
</msub>
</mrow>
</munder>
<msubsup>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
<mrow>
<mi>d</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>T</mi>
<mo>&GreaterEqual;</mo>
<msubsup>
<mi>G</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
<mi>C</mi>
</msubsup>
<mo>,</mo>
<mi>u</mi>
<mo>&Element;</mo>
<mi>U</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula,For the Contract Energy of the u class power supplys from region a to region b;
1-4-5) interregional interconnection capacity-constrained
This is constrained to the demand that interconnection channel capacity must assure that interregional exchange of electric power, i.e.,:
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>u</mi>
<mo>&Element;</mo>
<mi>U</mi>
</mrow>
</munder>
<msub>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
<mo>,</mo>
<mi>y</mi>
</mrow>
</msub>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<mo>&Element;</mo>
<mi>B</mi>
<mo>,</mo>
<mi>a</mi>
<mo>&NotEqual;</mo>
<mi>b</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
P in formula(a,b),yFor the active volume of the point-to-point power transmission passages of y regions a in equivalent network to region b, it is specially:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
<mo>,</mo>
<mi>y</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
<mo>,</mo>
<mn>0</mn>
</mrow>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>y</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>I</mi>
</msubsup>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
<mo>&Element;</mo>
<mi>B</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
P in formula(a,b),0For the existing capacity of the point-to-point power transmission passages of region a to region b in planning equivalent network at initial stage;
In addition, the extension of interconnection capacity also needs to consider that two interregional interconnection capacity have the upper limit, i.e.,:
P(a,b),y≤Pab max a,b∈A (17)
For the interconnection of point-to-point transmission in real network, interconnection capacity-constrained between above-mentioned zone can be directly applied;For
Interconnection in virtual equivalent network, formula (17) is constrained then equivalent to profile constraints, and (a, b) is represented between region a and region b
Virtual interconnection set;
1-4-6) interregional interconnection Constraint
Hourage can be utilized by being equivalent to its annual utilization hours to the year Constraint of each interconnection and must not exceed its year highest
T(a,b),y max, i.e.,:
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>u</mi>
<mo>&Element;</mo>
<mi>U</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>d</mi>
<mo>&Element;</mo>
<mi>D</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>&Element;</mo>
<msub>
<mi>T</mi>
<mi>D</mi>
</msub>
</mrow>
</munder>
<msubsup>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
<mrow>
<mi>d</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>T</mi>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
<mo>,</mo>
<mi>y</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
<mo>,</mo>
<mi>y</mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
1-4-7) typical day operation peak regulation and Reserve Constraint
This is constrained under given load curve of typical day, in each region have been switched on group need meet the peak load of one's respective area with
And peak regulation and stand-by requirement;On the premise of in a few days Unit Commitment is not considered, typical day operation constraint is embodied as:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>a</mi>
<mo>&Element;</mo>
<mi>A</mi>
</mrow>
</munder>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>u</mi>
<mo>&Element;</mo>
<msub>
<mi>U</mi>
<mi>d</mi>
</msub>
</mrow>
</munder>
<msub>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
</msub>
<mo>&GreaterEqual;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>b</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>d</mi>
<mi>max</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>a</mi>
<mo>&Element;</mo>
<mi>A</mi>
</mrow>
</munder>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>u</mi>
<mo>&Element;</mo>
<msub>
<mi>U</mi>
<mi>d</mi>
</msub>
</mrow>
</munder>
<msub>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>u</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>b</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>d</mi>
<mi>min</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<munder>
<mi>R</mi>
<mo>&OverBar;</mo>
</munder>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>y</mi>
<mo>&Element;</mo>
<mi>Y</mi>
<mo>,</mo>
<mi>d</mi>
<mo>&Element;</mo>
<mi>D</mi>
<mo>,</mo>
<mi>b</mi>
<mo>&Element;</mo>
<mi>A</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>19</mn>
<mo>)</mo>
</mrow>
</mrow>
3
In formula, UdThe set of the start unit in typical day d is represented, D is typical day set, αa,uRepresent that type is u's in a of region
Unit minimum output accounts for the accounting of rated capacity;Pb,y,d maxAnd Pb,y,d minRepresent in the b of y regions typical day d maximum and
Minimum load,WithRRepresent positive percentage reserve and negative percentage reserve;
1-4-8) typical day unit output constraint
This was constrained under given typical day, and all kinds of units of power system are exerted oneself no more than the active volume of the type unit:
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
<mrow>
<mi>d</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>&lambda;</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
<mi>d</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>20</mn>
<mo>)</mo>
</mrow>
</mrow>
Meanwhile, exerting oneself for all kinds of units of system can not be less than the minimum output of the type unit:
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
<mrow>
<mi>d</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&GreaterEqual;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msubsup>
<mi>&lambda;</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
<mi>d</mi>
</msubsup>
<mo>&CenterDot;</mo>
<msubsup>
<mi>&alpha;</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
<mi>d</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>21</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula,Represent y, the active volume of the u classes unit in a regions in typical day d accounts for Pu,y,a,bRatio;
Represent unit capacity minimum output ratio of the u classes unit in a regions in typical day d;
1-4-9) power generating facilities and power grids investment construction is constrained
This is constrained to the limitation of the power supply gone into operation every year each department and interzone interconnection capacity, should meet as shown in formula (22):
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
<mi>N</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mi>max</mi>
</mrow>
<mi>N</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
<mo>,</mo>
<mi>y</mi>
</mrow>
<mi>I</mi>
</msubsup>
<mo>&le;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>y</mi>
<mi>max</mi>
</mrow>
<mi>I</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>22</mn>
<mo>)</mo>
</mrow>
</mrow>
1-4-10) carbon emission transaction constraint
The constraint is used as carbon emission measurement criteria using electricity consumption rather than generating;The quota that power industry can purchase in one's respective area is deposited
In the upper limit, there is the upper limit in total electricity consumption carbon emission, i.e.,:
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>a</mi>
<mo>&Element;</mo>
<mi>A</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>u</mi>
<mo>&Element;</mo>
<mi>U</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>d</mi>
<mo>&Element;</mo>
<mi>D</mi>
</mrow>
</munder>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>&Element;</mo>
<msub>
<mi>T</mi>
<mi>D</mi>
</msub>
</mrow>
</munder>
<msubsup>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
<mrow>
<mi>d</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>T</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>e</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>max</mi>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>23</mn>
<mo>)</mo>
</mrow>
</mrow>
E in formulamax y,aThe electricity consumption carbon emission maximum permissible value of y after consideration carbon transaction in a of region is represented, can be by regional
Government independently formulates, when the government department of system region has clearly towards the carbon emission control targe of power industry,
Emaxy,aDirectly give;
1-4-11) newly-built power supply nonnegativity restriction
I.e. all kinds of power supply operation total amounts in region need to meet nonnegativity restriction, as shown in formula (24):
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>a</mi>
<mo>&Element;</mo>
<mi>A</mi>
</mrow>
</munder>
<msubsup>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
<mi>N</mi>
</msubsup>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>24</mn>
<mo>)</mo>
</mrow>
</mrow>
1-4-12) electricity and interconnection decision variable nonnegativity restriction, as shown in formula (25):
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>u</mi>
<mo>,</mo>
<mi>y</mi>
<mo>,</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
<mrow>
<mi>d</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msubsup>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>y</mi>
</mrow>
<mi>I</mi>
</msubsup>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>25</mn>
<mo>)</mo>
</mrow>
</mrow>
1-4-13) constraint of contact line justification, is stated as shown in formula (26):
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
<mo>&Element;</mo>
<mover>
<mi>L</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</munder>
<msubsup>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>,</mo>
<mi>b</mi>
<mo>)</mo>
<mo>,</mo>
<mi>y</mi>
</mrow>
<mi>I</mi>
</msubsup>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>26</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula,Represent the set in region pair exchanged in the absence of the energy;
2) to the solution of low-carbon power generating facilities and power grids Optimized model
The low-carbon power generating facilities and power grids Optimized model is expressed as linear model:
<mrow>
<mtable>
<mtr>
<mtd>
<mi>min</mi>
</mtd>
<mtd>
<mrow>
<msup>
<mi>c</mi>
<mi>T</mi>
</msup>
<mo>&CenterDot;</mo>
<mi>x</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>A</mi>
<mo>&CenterDot;</mo>
<mi>x</mi>
<mo>&le;</mo>
<mi>b</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<mi>x</mi>
<mo>=</mo>
<msub>
<mi>b</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<mi>x</mi>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>27</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein A and b characterizes the inequality constraints in optimization problem, AeqAnd beqThe equality constraint in optimization problem is characterized, x is decision-making
Variable, cTFor system cost term coefficient;
2-1) by step 1) each class cost is represented as in the object functionForm, by formula (2)~(9) construct
ck;Order:
<mrow>
<mi>c</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>28</mn>
<mo>)</mo>
</mrow>
</mrow>
K represents the number of cost in object function in formula, thus the c in constructive formula (27);
According to step 1-4) in all kinds of constraintss, each class equality constraint and inequality constraints are represented as Ai·x≤biOr
AeqjX=beqjForm;By step 1-4) construction Ai、biAnd Aeqi、beqi:Order:
<mrow>
<mi>A</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mi>I</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>b</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mi>I</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
<mi>J</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mi>b</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
<mi>J</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>29</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula (29), I and J represent the number of inequality constraints and equality constraint respectively, thus by the constraint bar in Optimized model
Part is converted to available matrix computations form;
Solver is worked out by MATLAB, calls CPLEX12.5 to solve the low-carbon power generating facilities and power grids Optimized model,
Try to achieve the power generating facilities and power grids optimization planning of the low-carbon of all kinds of power supplys and transmission line of electricity in power system so that power system is in low-carbon
There is optimal economic benefit under development model.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510041210.XA CN104573875B (en) | 2015-01-27 | 2015-01-27 | A kind of method of the power generating facilities and power grids optimization planning of low-carbon |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510041210.XA CN104573875B (en) | 2015-01-27 | 2015-01-27 | A kind of method of the power generating facilities and power grids optimization planning of low-carbon |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104573875A CN104573875A (en) | 2015-04-29 |
CN104573875B true CN104573875B (en) | 2017-09-26 |
Family
ID=53089890
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510041210.XA Active CN104573875B (en) | 2015-01-27 | 2015-01-27 | A kind of method of the power generating facilities and power grids optimization planning of low-carbon |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104573875B (en) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105976046B (en) * | 2016-04-27 | 2022-02-25 | 中国电力科学研究院 | Low-carbon power grid planning method considering demand side management |
CN107093029A (en) * | 2017-05-03 | 2017-08-25 | 国家电网公司 | A kind of power source planning method constrained based on Renewable energy policy control |
CN107317325A (en) * | 2017-06-28 | 2017-11-03 | 国网山东省电力公司经济技术研究院 | A kind of low-carbon generation Electric Power Network Planning model for considering carbon emission flow |
CN108062617B (en) * | 2017-11-28 | 2021-01-12 | 华电电力科学研究院有限公司 | Method for adjusting carbon emission quota of thermal power unit |
CN109840655A (en) * | 2017-11-29 | 2019-06-04 | 中国电力科学研究院有限公司 | A kind of low-carbon energy method and system for planning |
CN108847661B (en) * | 2018-06-11 | 2020-05-19 | 华中科技大学 | Annual production simulation operation method and system for regional power system |
CN109711728B (en) * | 2018-12-27 | 2023-04-18 | 陕西师范大学 | Double-layer multi-target power dispatching method based on power uncertainty and low-carbon appeal |
CN113344410B (en) * | 2021-06-21 | 2023-04-28 | 广东电网有限责任公司 | Multi-dimensional evaluation method, system, terminal and storage medium for energy conservation transformation of power distribution network |
CN113746144B (en) * | 2021-11-08 | 2022-03-22 | 南京东博智慧能源研究院有限公司 | Source-grid-load real-time interaction electric carbon control method and intelligent management system thereof |
CN114400716B (en) * | 2022-03-29 | 2022-07-22 | 华中科技大学 | Capacity selection method for hydropower station unit expansion |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102170129A (en) * | 2011-04-22 | 2011-08-31 | 武汉大学 | Low-carbon dispatching method and device for electric power system based on large-scale wind power grid connection |
CN103414213A (en) * | 2013-08-13 | 2013-11-27 | 中国能源建设集团广东省电力设计研究院 | Power supply optimal configuration method for sea island type micro-grid |
CN103544541A (en) * | 2013-10-15 | 2014-01-29 | 国家电网公司 | Carbon emission reduction evaluation and calculation method for intelligent power distribution and utilization system |
CN103762589A (en) * | 2014-01-08 | 2014-04-30 | 河海大学 | Method for optimizing new energy capacity ratio in layers in power grid |
CN104050609A (en) * | 2014-07-14 | 2014-09-17 | 天津大学 | Low-carbon comprehensive benefit analyzing system of grid connection photovoltaic power generation |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4203602B2 (en) * | 2006-09-06 | 2009-01-07 | 国立大学法人 琉球大学 | Operation support method and apparatus for power supply equipment |
-
2015
- 2015-01-27 CN CN201510041210.XA patent/CN104573875B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102170129A (en) * | 2011-04-22 | 2011-08-31 | 武汉大学 | Low-carbon dispatching method and device for electric power system based on large-scale wind power grid connection |
CN103414213A (en) * | 2013-08-13 | 2013-11-27 | 中国能源建设集团广东省电力设计研究院 | Power supply optimal configuration method for sea island type micro-grid |
CN103544541A (en) * | 2013-10-15 | 2014-01-29 | 国家电网公司 | Carbon emission reduction evaluation and calculation method for intelligent power distribution and utilization system |
CN103762589A (en) * | 2014-01-08 | 2014-04-30 | 河海大学 | Method for optimizing new energy capacity ratio in layers in power grid |
CN104050609A (en) * | 2014-07-14 | 2014-09-17 | 天津大学 | Low-carbon comprehensive benefit analyzing system of grid connection photovoltaic power generation |
Also Published As
Publication number | Publication date |
---|---|
CN104573875A (en) | 2015-04-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104573875B (en) | A kind of method of the power generating facilities and power grids optimization planning of low-carbon | |
Xu et al. | Blockchain-based trustworthy energy dispatching approach for high renewable energy penetrated power systems | |
Chen et al. | Increasing the flexibility of combined heat and power for wind power integration in China: Modeling and implications | |
Skarvelis-Kazakos et al. | Implementing agent-based emissions trading for controlling Virtual Power Plant emissions | |
Fan et al. | Review on coordinated planning of source-network-load-storage for integrated energy systems | |
CN106529737A (en) | Planning and distribution method for peak load regulation power source on supply side of power distribution network | |
CN101917024A (en) | Generating method of universality cost space in security-constrained dispatch | |
Yuan et al. | A new hybrid multi-criteria decision-making approach for developing integrated energy systems in industrial parks | |
Guo et al. | Multi-objective bi-level quantity regulation scheduling method for electric-thermal integrated energy system considering thermal and hydraulic transient characteristics | |
CN105243600A (en) | Grid power generation adjustment method | |
Yu et al. | An improved combined heat and power economic dispatch model for natural gas combined cycle power plants | |
CN105528668A (en) | Dynamic environment and economy scheduling method of grid-connected wind power system | |
Pan et al. | The state grid corporation of China's practice and outlook for promoting new energy development | |
Nozarian et al. | Reliability-oriented planning framework for smart cities: From interconnected micro energy hubs to macro energy hub scale | |
Li et al. | Energy supply reliability assessment of the integrated energy system considering complementary and optimal operation during failure | |
Yang et al. | Regional integrated energy system reliability and low carbon joint planning considering multiple uncertainties | |
Yang et al. | Optimal scheduling of electro-thermal system considering refined demand response and source-load-storage cooperative hydrogen production | |
Ding et al. | Long-term operation rules of a hydro–wind–photovoltaic hybrid system considering forecast information | |
Pang et al. | Drivers and key pathways of the household energy consumption in the Yangtze river economic belt | |
CN117933648A (en) | Comprehensive energy system scheduling method based on double-layer carbon transaction model | |
Yang et al. | Optimal planning of energy storage system under the business model of cloud energy storage considering system inertia support and the electricity-heat coordination | |
Pang et al. | Optimal scheduling of a cogeneration system via Q-learning-based memetic algorithm considering demand-side response | |
CN105633950A (en) | Multi-target random, fuzzy and dynamic optimal power flow considering wind power injection uncertainty | |
Li et al. | Energy system and scheduling strategies of electric–gas networks deeply coupled under dual electricity prices | |
Chen et al. | A Low‐Carbon Planning Model for Regional Power Systems with Generation‐Load‐Storage Coordination considering New Energy Resources’ Consumption |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |