CN106127342B - Cross-regional tie line transaction optimization method based on SCUC - Google Patents
Cross-regional tie line transaction optimization method based on SCUC Download PDFInfo
- Publication number
- CN106127342B CN106127342B CN201610473659.8A CN201610473659A CN106127342B CN 106127342 B CN106127342 B CN 106127342B CN 201610473659 A CN201610473659 A CN 201610473659A CN 106127342 B CN106127342 B CN 106127342B
- Authority
- CN
- China
- Prior art keywords
- tie
- power
- line
- time
- peak
- 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
- 238000005457 optimization Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000005611 electricity Effects 0.000 claims description 18
- 230000005540 biological transmission Effects 0.000 claims description 8
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 238000012821 model calculation Methods 0.000 abstract 1
- 238000011161 development Methods 0.000 description 8
- 238000005265 energy consumption Methods 0.000 description 4
- 230000014509 gene expression Effects 0.000 description 4
- 230000001737 promoting effect Effects 0.000 description 4
- 230000010354 integration Effects 0.000 description 3
- 238000010248 power generation Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000013468 resource allocation Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 102100029469 WD repeat and HMG-box DNA-binding protein 1 Human genes 0.000 description 1
- 101710097421 WD repeat and HMG-box DNA-binding protein 1 Proteins 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000009194 climbing Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention provides a trans-regional junctor trading optimization method based on SCUC, which comprises the steps of setting parameters of a power grid, a unit and a trans-regional junctor of an input and output area; selecting a cross-regional power transaction plan compiling method; and constructing an input and output combined safety constraint unit combination model and optimizing model calculation. The technical scheme provided by the invention improves the pertinence of the electric power transaction plan compilation in different regions and different power supply types, can evaluate whether the capacity of the cross-regional connecting line meets the capacity expansion requirement, and provides technical support for the capacity expansion of the connecting line.
Description
Technical Field
The invention relates to an optimization method in the fields of electric power markets and electric power automation, in particular to a cross-regional tie line transaction optimization method based on SCUC.
Background
The government of China highly pays attention to the development and utilization of clean energy, and even promotes the power generation of the clean energy and the optimized configuration of the large-range clean energy to the national strategic altitude. With the development of a series of new energy development policies in China, large-scale clean energy development and utilization become important measures for improving energy structures, building low-carbon society and maintaining economic and social sustainable development.
Accompanying the large-scale integration of clean energy is a series of consumption difficulties of clean energy: on one hand, as large-scale clean energy is far away from a load center and is limited by the load level of a regional power grid and a power supply structure, the dispatching operation and the receiving pressure of the clean energy are huge, and the wind and water abandoning conditions are serious; on the other hand, the construction of the main power plant in the end system is not matched with the increasing speed of the power demand, so that large-area power limitation appears in multi-province cities, a large amount of electric energy needs to be transmitted remotely for relieving the power utilization pressure, and the dependence degree of the input end system on external power is continuously improved. The cross-regional transaction and large-scale consumption of the clean energy are effective ways for solving the problem that the clean energy is difficult to be consumed on site, and the dependence of an input end system on external power is met.
The global energy Internet thought is the key of clean energy substitution and electric energy substitution, and is a platform for guaranteeing the high-efficiency development of global clean energy; the development planning of national synchronous power grids of national power grids is also an effective way for promoting resource optimization configuration and peak shaving resource complementation. With the gradual construction of the ultra-high voltage backbone network of the three-horizontal three-vertical one-ring network, the transmission capacity between networks is rapidly enhanced, and a wide space is provided for large-scale cross-regional resource optimization configuration.
Aiming at the urgent need of developing clean energy in China and the contradiction between the current situation of a power grid and the consumption of the clean energy in China, the lean decision of cross-regional transaction power is realized by optimizing cross-regional clean energy transaction power, fully utilizing the existing power transmission channel and creating a transaction mode aiming at different regions and different power types from the perspective of large-scale resource allocation, the condition that the existing cross-regional transaction power depends on artificial decision is improved, and the scheme and the plan which promote the reliable consumption of the clean energy through a marketization way and have practical significance are researched.
In order to meet the needs of the prior art, the invention provides a cross-regional tie line transaction optimization method based on the SCUC.
Disclosure of Invention
To meet the needs of the prior art, the invention provides a tie line power trading scheme which optimally trades power across a tie line and promotes maximum consumption of clean energy.
The improvement of the cross-regional junctor trading optimization method provided by the invention is that the optimization method comprises the following steps:
(1) setting parameters of input and output areas;
(2) selecting a cross-regional power transaction plan compiling method;
(3) constructing an input and output combined safety constraint unit combination model;
(4) and optimizing the model.
Further, in the step (1),
(1-1) determining an input area, an output area, a trans-regional direct current tie line between the input area and the output area and other tie lines;
and (1-2) determining physical and economic parameters of a power grid, a unit and a cross-region tie line of an input region and an output region.
Further, the step (2) comprises,
(2-1) determining the peak-valley time period of the tie line power, and optimizing the peak-valley proportion according to the following formula;
the valley power is represented by the following formula (1):
PtieD(itie)＝RD*PMax(itie) (1)
the peak power is shown in the following equation (2):
PtieU(itie)＝RU*PMax(itie) (2)
wherein, the relationship between the peak proportion and the valley proportion is shown as the following formula (3):
RU≥RD (3)
in the above formula, PtieD(itie): tie line itieA power value at a trough period; ptieU(itie): tie line itieThe power value during peak hours; pMax(itie): tie line itieThe maximum delivery capacity of; rD: the ratio of the valley power to the maximum transmission power of the tie line; rU: the ratio of peak power to maximum power delivered by the tie; pMax(itie): tie line itieThe maximum delivery capacity of;
(2-2) determining the peak-to-valley proportion of the tie line power, and optimizing the peak-to-valley time period according to the following formula;
the full optimization time period only has one time of changing the valley into the peak and one time of changing the peak into the valley;
junctor state Utie(itieT) is represented by the following formula (4):
tie line power Ptie(itieT) is represented by the following formula (5):
wherein, UstartUp(itieT): tie line itieThe change from trough to peak at time t, represented by 0 or 1; u shapestartDown(itieT): tie line itieThe change from peak to valley during time t, represented by 0 or 1; rU(itie) And RD(itie): respectively representing the ratio of the low-valley power to the high-peak power on the contact line; u shapetie(itieT): if U istie(itieT) 1, representing a tie itieIs in a peak period during the period t; if U istie(itieT) is 0 and represents a tie line itieIn a valley period at a time t; pTieMax(itie): tie line itieMaximum power limit of;
(2-3) the peak-to-valley ratio is not limited, and the optimization is free.
Further, the step (3) includes a combined model of the purchasing, selling and transporting combined safety constraint unit established based on the input and output combined model;
the objective function of the combined model is shown in the following equation (6):
wherein max CleanEnergy: the maximum consumption of clean energy; i.e. ice: a clean energy unit; p is a radical ofi(iceT, b): section i of zone bceThe force at time t; t isPrdMin: a time interval comprises a number of times in units of minutes.
Further, the step (3) comprises:
(3-1) thermal power minimum operation mode constraint;
the minimum number of the regional thermal power generating units is shown as the following formula (7):
the minimum starting capacity of the regional thermal power generating unit is shown as the following formula (8):
in the above formula, NumMinOn(g, b): the minimum number of the starting units of the unit group g in the area b; capMinOn(g, b): minimum boot capacity of cluster g of zone b; i issGU(iceG, b): is represented by 0 or 1, IsGU(iceG, b) 1 indicates the unit i of the area bceIn cluster b; u shapei(iceT, b): section i of zone bceA start-stop state at time t;
(3-2) constraint to guarantee input area acceptance intent:
the input area abandoned wind power does not increase as shown in the following formula (9):
Wherein, PFW(t, b): inputting the predicted wind power of the region b at the moment t; pi(iwdT, b): wind turbine i of input region bwdThe force at time t; rAW(b) The method comprises the following steps The air abandoning proportion of the area b; setBrchbuy: collecting electricity purchasing areas;
electricity purchase cost C of input areaTB(b) Without increasing as shown in the following formula (10):
wherein, Ci(i, t, b) the cost of electricity generated by the wind turbine i in the area b at the time t; si(i, t, b) represents the start-stop cost of the wind turbine generator i in the area b at the moment t; dtie(itieT, b) a tie line itieIn the direction of the region b at the time t, dtie(itieT, b) indicates acceptance, dtie(itieT, b) — 1 represents sending; pr (Pr) oftie(itieT, b) denotes a tie itieThe price in region b at time t; i issTieOp(itie) Represented by 0 or 1, IsTieOp(itie) 1 represents itieIs an optimized tie line variable; cBThe original electricity purchasing cost;
further, the operation constraint of the dc link in step (3) includes:
i. the tie line capacity is shown by the following formula (11):
0≤Ptie(itie,t)≤PTieMax(itie) (11)
in the formula, PTieMax(itie): maximum power limit of the tie line;
ii. The adjustment rate of the dc line in the adjacent time period is shown in the following formula (12):
in the formula, RampUp(itieT): tie line itieA ramp rate at time t; rampDown(itieT): tie line itieA landslide rate at time t; Δ t: the interval length of the time period;
iii, the increase and decrease change of the power of the tie line in the adjacent time period cannot be adjusted in different directions;
the link power adjustment state is represented by the following expression (13):
x+(itie,t)+x-(itie,t)＝x(itie,t)≤1 (13)
the change in the tie line power value is shown in the following equation (14):
in the above formula, x+(itieT): whether the direct current sending power changes in the positive direction or not at each time interval; x- (i)tieT): whether the direct current sending power changes reversely at each time interval; z is a radical of1(itieT) and z2(itieT) is an auxiliary variable represented by 0 or 1; m1And M2Is an auxiliary positive value parameter;
iv, the interval of the tie line dc adjustment is represented by the following formula (15):
tie line itieThe state at time t is determined according to the following equation (16):
in the above formula, NT: DC line itieThe minimum number of adjustment interval periods; y (i)tieAnd t) is an auxiliary variable represented by 0 or 1.
Further, in the step (4), the nonlinear factors in the input and output combined safety constraint unit combination model are linearly expressed, and the electric power of the cross-region connecting line and the unit start-stop and output of the input and output ends are calculated by adopting a mixed integer programming method.
Moreover, compared with the closest prior art, the invention has the following excellent effects:
1) the technical scheme provided by the invention takes the maximum consumption of the clean energy as an optimization target, promotes the maximum consumption of the clean energy by optimizing the cross-regional connecting line trading power and the unit starting, stopping and outputting of the input and output regions, and improves the receiving capacity of the clean energy; the method comprises the steps of establishing a purchase, sale and transmission combined optimization model which comprehensively considers an input area, an output area and a power transmission channel and covers different types of constraints such as physical operation, economic operation and the like, realizing quantitative calculation of the cross-regional consumption of the clean energy, and solving the problem of evaluation of the cross-regional consumption capacity of the clean energy in a long-time scale.
2) The technical scheme provided by the invention aims at different regions and different power types, utilizes the peak-valley difference, the time difference and the load difference of power resources among the regions, fully utilizes the existing power transmission channel and optimizes the power curve of the interconnection line among the regions; developing a clean energy trans-regional trans-provincial mid-long term power transaction optimization algorithm, and realizing inter-regional tie line power curve optimization by using peak-valley difference, time difference and load difference of inter-regional power resources, so as to improve the consumption level of clean energy; the technical means is provided for the arrangement of the cross-regional clean energy trading mode, and the lean operation level of the clean energy trading is improved.
3) The technical scheme provided by the invention is based on the safety constraint unit combination, can simulate the core link of generating plan compilation, can improve the operation efficiency of the high-energy-efficiency large unit, enables the unit to operate at the optimal working point as far as possible, improves the economy and energy conservation of the system, and provides a basis for economic dispatching and safety check of the power system; the wind power consumption capability of an output end power grid and the wind power receiving capability of an input end power grid are optimized and evaluated, the risk brought by large-scale wind power integration to the safe operation of the power grid can be prevented in advance, the wind power receiving capability of the power grid under the current operation environment can be deeply sensed, effective technical support can be provided for scheduling personnel to make day-ahead power generation plans and control the real-time operation, and the safe operation level and the wind power consumption capability of the power grid after large-scale wind power integration are greatly improved.
4) The technical scheme of the invention provides a method for compiling the cross-regional and cross-provincial power trading plan of clean energy with fixed peak-valley ratio, fixed peak-valley time period, unlimited curve shape and the like, realizes multi-scene comparison of different consumption schemes of the clean energy, improves the pertinence of the power trading plan compilation in different regions and different power types, can optimize and evaluate whether the capacity of the cross-regional connecting line meets the requirement of capacity expansion, and provides technical support for the capacity expansion of the connecting line.
Drawings
FIG. 1 is a flow chart of an optimization method provided by the present invention.
Detailed Description
The technical scheme provided by the invention is clearly described in detail in the following with reference to the attached drawings of the specification.
The invention relates to a cross-regional tie line transaction optimization method for promoting clean energy consumption based on SCUC, which can provide an optimal cross-regional tie line transaction electric power, uniformly and coordinately consider the complementarity of load difference, peak-valley difference and time difference of a transmitting end and a receiving end of the clean energy, and select a tie line electric power transaction scheme and a unit combination scheme for promoting the maximum consumption of the clean energy.
Based on an input and output area interconnection system, the physical economic operation constraint of power grid operation is fully considered, the minimum operation mode constraint of a thermal power generating unit, the tie line electric power transaction scheme constraint and the direct current tie line operation constraint are considered, and the tie line electric power transaction capable of promoting clean energy consumption is optimized through a safety constraint unit combination method.
Aiming at the urgent need of developing clean energy in China and the contradiction between the current situation of a power grid and the consumption of the clean energy in China, the lean decision of cross-regional power trading is realized by optimizing cross-regional clean energy trading power, fully utilizing the existing power transmission channel aiming at different regions and different power types, further innovating a trading mode from the perspective of large-scale resource allocation, and selecting a scheme and a plan which promote the reliable consumption of the clean energy and have practical significance through a marketization approach. The invention provides technical support for lean formulation of a trading scheme, and effectively solves the problem that most trans-regional trading power curves of a power grid are decided by artificial experience and the power supply load difference of an input and output region cannot be considered comprehensively in a refined manner.
The optimization method provided by the invention comprises the following steps:
firstly), setting an input area and an output area of clean energy, and determining physical and economic parameters of a power grid, a unit and a cross-regional connecting line of the output and input areas;
the method is characterized in that the direction of the trans-regional transaction of the clean energy is determined, namely a trans-regional direct current link and other links between the output region and the input region of the clean energy are determined.
Specific physical and economic parameters are required including:
(1) inputting information of each unit in an output area, wherein the information comprises unit name, installed capacity, unit type, plant power rate, peak regulation capacity, whether the unit is a heat supply unit or not, minimum output of the unit, climbing rate, landslide rate, minimum outage time and minimum running time;
(2) inputting recent power supply planning information of an output area, wherein the recent power supply planning information comprises newly-added machine capacity and newly-added machine types;
(3) inputting the unit electricity price or the marking post electricity price of the output area;
(4) generating cost information of the thermal power generating unit comprises a micro-increment cost curve, cold start cost, warm start cost and hot start cost;
(5) inputting and outputting minimum operation mode information of the thermal power generating units in the area, wherein the minimum operation mode information comprises the minimum starting number and the minimum starting capacity of each power plant;
(6) inputting system load and rotation standby information of an output area;
(7) inputting wind and light predicted output of an output area;
(8) inputting the predicted output of the radial flow hydropower of the output area and the generated electric quantity range of the reservoir type hydropower;
(9) basic information of the trans-provincial/trans-regional junctor comprises a junctor name, a power upper limit and a power lower limit;
(10) the inter-provincial/inter-regional transaction information comprises inter-provincial/inter-regional link names, transaction buyers, transaction sellers, transaction components, transaction prices, transaction electric quantity and transaction electric power;
(11) the intra-provincial/regional framework structure information comprises a section name, a section limit value, units related to the section, sensitivity coefficients of the units corresponding to the section, and sensitivity coefficients of the links corresponding to the section;
(12) inputting wind abandoning information, light abandoning information and water abandoning information of an area; the cost information of purchasing electricity from outside the area is input.
Secondly), selecting a clean energy cross-region and cross-provincial trading plan compiling method according to requirements:
fixing the peak-valley time period of the electric power of the tie line, and optimizing the peak-valley proportion; fixing the peak-valley proportion of the electric power of the tie line, and optimizing the peak-valley time period; the peak-to-valley ratio is not limited, and the optimization is free;
in the actual power grid operation and transaction development process, according to the operation experience of the tie line or the negotiation result of both the purchasing and selling parties, the shape of the direct current tie line power is constrained, such as: fixing the peak-valley time period of the electric power of the tie line, and optimizing the peak-valley proportion; the peak-valley proportion of the electric power of the tie line is fixed, the peak-valley period is optimized, and the like, and the mathematical expressions corresponding to the above shape constraints are as follows:
1) fixing the peak-valley time period of the electric power of the tie line, and optimizing the peak-valley proportion;
the relationship between the valley ratio and the valley power is shown in the following formula (1):
PtieD(itie)＝RD·PMax(itie) (1)
the relationship between the peak proportion and the peak power is shown in the following formula (2):
PtieU(itie)＝RU·PMax(itie) (2)
the peak ratio is equal to or greater than the trough ratio, as shown in the following formula (3):
RU＞＝RD (3)
assigning peak power to the tie-line power during peak hours is shown in equation (4) below:
Ptie(itie,t)＝PtieU(itie) Wherein t satisfies IsTieU(t)＝1 (4)
The tie-line power for assigning the valley power to the valley period is shown in the following equation (5):
Ptie(itie,t)＝PtieD(itie) Wherein t satisfiesIsTieU(t)＝0 (5)
Wherein, PtieD(itie) Represents a tie line itiePower value in the valley period, PtieU(itie) Represents a tie line itieElectric power value during peak hours, PMax(itie) Represents a tie line itieMaximum transport capacity, RDRepresenting the ratio of the valley power to the maximum power delivered by the tie, RURepresenting the ratio of the valley power to the maximum power delivered by the tie, Ptie(itieT) represents the tie line power, IsTieU(t) '1' means that the period t is in the peak period, IsTieUWhen (t) is 0, the time period t is in the valley period;
2) fixing the peak-valley proportion of the electric power of the tie line, and optimizing the peak-valley time period;
the time when the full optimization period has only one valley to peak is shown as the following formula (6):
the time when the peak becomes low is only once in the full optimization period as shown in the following formula (7):
the relationship between the tie line state and the peak-to-valley state is shown in the following equation (8):
the relationship between the tie line power, the tie line state, and the peak-to-valley ratio is as shown in the following equation (9):
Ptie(itie，t)＝Utie(itie，t)·PTieMax(itie)·RU(itie)+(1-Utie(itie,t))·PTieMax(itie)·RD(itie) (9)
wherein, UstartUp(itie,t)、UstartDown(itieT) represents a tie line i by 0 or 1, respectivelytieNo change from valley to peak, tie-line i during time ttieNo change from peak to valley during time t; rUAnd RDRepresenting the ratio of valley to peak power on the tie line; u shapetie(itieT): if U istie(itieT) 1, representing a tie itieIs in a peak period during the period t; if U istie(itieT) is 0 and represents a tie line itieIn the valley period at the t period.
Thirdly), modeling a clean energy output area and an input area into a transmitting and receiving combined system through a cross-region connecting line, and establishing a safety constraint unit combined model;
according to the actual power grid model of the transmitting and receiving ends and the actual physical characteristics of the connecting lines, practical constraints such as thermal power minimum operation mode constraint, clean energy prediction constraint, hydroelectric power generation quantity constraint, direct-current connecting line power shape constraint, input region wind-abandoning light-abandoning water-abandoning power quantity constraint, input region electricity purchasing cost constraint and the like, and basic constraints such as system balance constraint, direct-current connecting line operation constraint, unit operation constraint, power grid safety constraint and the like are considered, the maximum clean energy consumption of the input end and the output end is taken as a target, and the optimized objects are the transaction power of the cross-regional connecting line and the unit start-stop and output of the input end and the output end;
establishing a combined model of a purchase, sale and transportation combined safety constraint unit based on an input and output combined model, wherein the target function is the maximum clean energy consumption of an output area, and the expression is as follows:
wherein iceRepresenting a clean energy unit; p is a radical ofi(iceT, b) stands for a unit i of the area bceForce applied at time T, TPrdMinIs to indicate a timeThe interval includes minutes.
(1) Thermal power minimum operation mode constraint:
the constraint of the minimum starting number of the regional thermal power generating units is shown as the following formula (11): :
the constraint of the minimum startup capacity of the regional thermal power generating unit is shown as the following formula (12):
wherein, formula (11) represents that the minimum number of the thermal power generator groups g in the region b is not less than NumMinOn(g, b); equation (12) represents that the minimum startup capacity of the thermal power unit group g of the region b is not less than CapMinOn(g,b)；IsGU(ituG, b) represents 0 or 1: i issGU(ituG, b) 1 indicates the unit i of the area bceIn cluster g; u shapei(ituT, b) stands for a unit i of the area bceA start-stop state at time t; n is a radical ofumMinOn(g, b) represents the minimum number of the sets of the set group g of the area b; capMinOn(g, b) represents the minimum boot capacity of the cluster g of the area b.
(2) Constraints that guarantee the acceptance willingness of the input area:
firstly, the input region abandoned wind power is not restricted, and a restriction expression that the input region abandoned wind power is not greater than the original abandoned wind power is shown as the following formula (13):
wherein, PFW(t, b) represents the predicted wind power of input region b at time t, Pi(iwdT, b) wind turbines i representing the input area bwdOutput at time t, RAW(b) Showing the wind curtailment ratio, S, of the area betBrchbuyTo representAnd (4) collecting electricity purchasing areas.
② the electricity purchasing cost of income area does not increase the restriction, and the electricity purchasing cost C of area b is inputtedTB(b) The electricity purchasing cost of the area b is not more than the original electricity purchasing cost CBThe constraint is shown in the following equation (14),
wherein itieSatisfy IsTieOp(itie)≠0；IsTieOp(itie) Represented by 0 or 1, IsTieOp(itie) 1 represents itieIs an optimized tie line variable; ci(i, t, b) represents the cost of electricity generated by the wind turbine i in zone b at time t, Si(i, t, b) represents the start-stop cost of the wind turbine generator i of the area b at the moment t, dtie(itieT, b) denotes a tie itieIn the direction of the area b at the time t, 1 indicates incoming and-1 indicates outgoing; pr (Pr) oftie(itieT, b) denotes a tie itieThe price in region b at time t; . (3) The operation constraint of the direct current tie line is as follows:
1) the tie capacity is shown by the following formula (15):
0≤Ptie(itie,t)≤PTieMax(itie) (15)
the above equation constrains the power of the tie-line to be within its maximum power limit.
2) Tie line power rate of change constraint
Ptie(itie,t)-Ptie(itie,t-1)≤RampUp(itie,t)Δt (15)
Ptie(itie,t-1)-Ptie(itie,t)≤RampDown(itie,t)Δt (16)
Equations (15) and (16) constrain that the adjustment rate of the dc line in the adjacent time period cannot exceed the limit of the dc operation mode; rampUp(itieAnd t) represents a tie line itieA ramp rate at time t; rampDown(itieT) represents a slip rate; Δ t represents the interval length of the period.
3) Link continuous time interval power adjustment direction constraint
The increase and decrease changes of adjacent time intervals cannot be adjusted in different directions as shown in the following formulas (18) and (19):
x+(itie,t)+x-(itie,t+1)≤1 (17)
x+(itie,t+1)+x-(itie,t)≤1 (18)
junctor power adjustment state x (i)tieAnd t) is related to the forward and reverse changes of power as shown in the following equation (20):
x+(itie,t)+x-(itie,t)＝x(itie,t)≤1 (19)
the relationship between the change in the tie line power value and the variables of the forward change and the direction change is expressed by the following equations (21), (22) and (23):
Ptie(itie,t)-Ptie(itie,t-1)≤M1z1(itie,t) (20)
Ptie(itie,t-1)-Ptie(itie,t)≤M2z2(itie,t) (21)
x+(itie,t)≥z1(itie,t) (22)
wherein, x (i)tie,t),x+(itie,t),x-(itieT), respectively adopting 0 or 1 to respectively indicate whether the direct current sending power changes in each time interval, whether the direct current sending power changes in a forward direction (power is increased), and whether the direct current sending power changes in a reverse direction (power is reduced); z is a radical of1(itieT) and z2(itieT) is an auxiliary variable, represented by 1 or 0; m1And M2Is an auxiliary positive value parameter.
(4) Adjusting interval constraint by the direct current of the tie line;
after the dc power is adjusted once, the dc power is operated steadily for at least a minimum time interval as shown in the following equation (24):
wherein, IsTieStart(itieT) and IsTieEnd(itieT) represents a tie line i by 0 or 1, respectivelytieWhether power adjustment is started or not and whether power adjustment is finished or not at the moment t; n is a radical ofTFor a direct current line itieThe minimum number of adjustment interval periods; y (i)tieT) is an auxiliary variable, represented by 0 or 1;
IsTieStart(itie,t)、IsTieEnd(itiet) and x (i)tie,t)、y(itieAnd t) are shown in the following formulas (25) to (29):
IsTieStart(itie,t)≥x(itie,t+1)-y(itie,t) (24)
IsTieEnd(itie,t)≥x(itie,t)-y(itie,t) (25)
y(itie,t)≤x(itie,t) (26)
y(itie,t)≤x(itie,t+1) (27)
y(itie,t)≥x(itie,t)+x(itie,t+1)-1 (28)
and fourthly) carrying out optimization calculation aiming at the input and output combined safety constraint unit combination model.
And (3) linearly expressing the nonlinear factors in the model, and calculating the power of the cross-region connecting line and the unit start-stop and output of the input end and the output end by adopting a mixed integer programming method.
Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art can make modifications and equivalents to the embodiments of the present invention without departing from the spirit and scope of the present invention, which is set forth in the claims of the present application.
Claims (2)
1. A cross-regional junctor transaction optimization method based on SCUC is characterized by comprising the following steps:
(1) setting parameters of input and output areas;
(2) selecting a cross-regional power transaction plan compiling method;
(3) constructing an input and output combined safety constraint unit combination model;
(4) calculating an optimization model;
in the step (1), the step (c),
(1-1) determining an input area, an output area, a trans-regional direct current tie line between the input area and the output area and other tie lines;
(1-2) determining physical and economic parameters of a power grid, a unit and a cross-region tie line of an input region and an output region;
the step (2) comprises the steps of,
(2-1) determining the peak-valley time period of the tie line power, and optimizing the peak-valley proportion according to the following formula;
the valley power is represented by the following formula (1):
PtieD(itie)＝RD*PMax(itie) (1)
the peak power is shown in the following equation (2):
PtieU(itie)＝RU*PMax(itie) (2)
wherein, the relationship between the peak proportion and the valley proportion is shown as the following formula (3):
RU≥RD (3)
in the above formula, PtieD(itie): tie line itieA power value at a trough period; ptieU(itie): tie line itieThe power value during peak hours; pMax(itie): tie line itieThe maximum delivery capacity of; rD: the ratio of the valley power to the maximum transmission power of the tie line; rU: the ratio of peak power to maximum power delivered by the tie; pMax(itie): tie line itieMaximum transport ofCapacity;
(2-2) determining the peak-to-valley proportion of the tie line power, and optimizing the peak-to-valley time period according to the following formula;
the full optimization time period only has one time of changing the valley into the peak and one time of changing the peak into the valley;
junctor state Utie(itieT) is represented by the following formula (4):
tie line power Ptie(itieT) is represented by the following formula (5):
Ptie(itie,t)＝Utie(itie,t)*PTieMax(itie)*RU(itie)+(1-Utie(itie,t))*PTieMax(itie)*RD(itie) (5)
wherein, UstartUp(itieT): tie line itieThe change from trough to peak at time t, represented by 0 or 1; u shapestartDown(itieT): tie line itieThe change from peak to valley during time t, represented by 0 or 1; rU(itie) And RD(itie): respectively representing the ratio of the low-valley power to the high-peak power on the contact line; u shapetie(itieT): if U istie(itieT) 1, representing a tie itieIs in a peak period during the period t; if U istie(itieT) is 0 and represents a tie line itieIn a valley period at a time t; pTieMax(itie): tie line itieMaximum power limit of;
(2-3) the peak-to-valley ratio is not limited, and the optimization is free;
the step (3) comprises a purchase, sale and transportation combined safety constraint unit combined model established based on the input and output combined model;
the objective function of the unit combination model is shown as the following formula (6):
wherein max CleanEnergy: the maximum consumption of clean energy; i.e. ice: a clean energy unit; p is a radical ofi(iceT, b): section i of zone bceThe force at time t; t isPrdMin: a time interval comprises a number of times in units of minutes;
the step (3) comprises the following steps:
(3-1) thermal power minimum operation mode constraint;
the minimum number of the regional thermal power generating units is shown as the following formula (7):
the minimum starting capacity of the regional thermal power generating unit is shown as the following formula (8):
in the above formula, NumMinOn(g, b): the minimum number of the starting units of the unit group g in the area b; capMinOn(g, b): minimum boot capacity of cluster g of zone b; i issGU(iceG, b): is represented by 0 or 1, IsGU(iceG, b) 1 indicates the unit i of the area bceIn cluster b; u shapei(iceT, b): section i of zone bceA start-stop state at time t;
(3-2) constraint to guarantee input area acceptance intent:
the input area abandoned wind power does not increase as shown in the following formula (9):
wherein, PFW(t, b): inputting the predicted wind power of the region b at the moment t; pi(iwdT, b): wind turbine i of input region bwdThe force at time t; rAW(b) The method comprises the following steps The air abandoning proportion of the area b; setBrchbuy: collecting electricity purchasing areas;
electricity purchase cost C of input areaTB(b) Without increasing as shown in the following formula (10):
wherein itieSatisfy IsTieOp(itie)≠0；IsTieOp(itie) Represented by 0 or 1, IsTieOp(itie) 1 represents itieIs an optimized tie line variable; ci(i, t, b) the cost of electricity generated by the wind turbine i in the area b at the time t; si(i, t, b) represents the start-stop cost of the wind turbine generator i in the area b at the moment t; dtie(itieT, b) a tie line itieIn the direction of the region b at the time t, dtie(itieT, b) indicates acceptance, dtie(itieT, b) — 1 represents sending; pr (Pr) oftie(itieT, b) denotes a tie itieThe price in region b at time t; cBThe original electricity purchasing cost;
the operation constraint of the direct current tie line in the step (3) comprises the following steps:
i. the tie line capacity is shown by the following formula (11):
0≤Ptie(itie,t)≤PTieMax(itie) (11)
in the formula, PTieMax(itie): maximum power limit of the tie line;
ii. The adjustment rate of the dc line in the adjacent time period is shown in the following formula (12):
in the formula, RampUp(itieT): tie line itieA ramp rate at time t; rampDown(itieT): tie line itieA landslide rate at time t; Δ t: the interval length of the time period;
iii, the increase and decrease change of the power of the tie line in the adjacent time period cannot be adjusted in different directions;
junctor power adjustment state x (i)tieT) is represented by the following formula (13):
x+(itie,t)+x-(itie,t)＝x(itie,t)≤1 (13)
the change in the tie line power value is shown in the following equation (14):
in the above formula, x+(itieT): whether the direct current sending power changes in the positive direction or not at each time interval; x is the number of-(itieT): whether the direct current sending power changes reversely at each time interval; z is a radical of1(itieT) and z2(itieT) is an auxiliary variable represented by 0 or 1; m1And M2Is an auxiliary positive value parameter;
iv, the interval of the tie line dc adjustment is represented by the following formula (15):
wherein, IsTieStart: tie line itieWhether the active power adjustment is started at the time t is represented by 0 or 1;
IsTieEnd: tie line itieWhether the active power adjustment is finished at the time t is represented by 0 or 1;
tie line itieThe state at time t is determined according to the following equation (16):
in the above formula, NT: DC line itieThe minimum number of adjustment interval periods; y (i)tieAnd t) is an auxiliary variable represented by 0 or 1.
2. The optimization method according to claim 1, wherein in the step (4), the nonlinear factors in the input and output combined safety constraint unit combination model are linearly expressed, and a mixed integer programming method is adopted to calculate the power of the cross-regional connecting line and the unit start-stop and output of the input and output ends.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610473659.8A CN106127342B (en) | 2016-06-24 | 2016-06-24 | Cross-regional tie line transaction optimization method based on SCUC |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610473659.8A CN106127342B (en) | 2016-06-24 | 2016-06-24 | Cross-regional tie line transaction optimization method based on SCUC |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106127342A CN106127342A (en) | 2016-11-16 |
| CN106127342B true CN106127342B (en) | 2021-03-02 |
Family
ID=57266695
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610473659.8A Active CN106127342B (en) | 2016-06-24 | 2016-06-24 | Cross-regional tie line transaction optimization method based on SCUC |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106127342B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106846172A (en) * | 2016-12-26 | 2017-06-13 | 国网山东省电力公司泰安供电公司 | Tie line plan power acquisition methods and device based on electricity transaction |
| CN108347048B (en) * | 2017-01-22 | 2024-03-19 | 中国电力科学研究院 | Planning method adapting to transregional and transnational scheduling modes |
| CN107798430B (en) * | 2017-10-25 | 2021-02-26 | 国网新疆电力公司电力科学研究院 | Bidding optimization method considering renewable energy cross-region consumption |
| CN107968408B (en) * | 2017-11-17 | 2021-05-04 | 中国南方电网有限责任公司 | Asynchronous networking direct-current power plan optimization method, system and device |
| CN108182506B (en) * | 2017-12-05 | 2021-08-27 | 中国电力科学研究院有限公司 | Step-by-step optimization method and device for monthly power generation plan |
| CN111815176B (en) * | 2020-07-10 | 2022-04-15 | 国网四川省电力公司电力科学研究院 | Long-term electric quantity multi-channel complementary coordination sending method and system in hydropower enrichment power grid |
| CN112200408A (en) * | 2020-09-03 | 2021-01-08 | 中国南方电网有限责任公司 | Method, system, device and medium for clearing cross-regional electric power spot market |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102280878A (en) * | 2011-07-26 | 2011-12-14 | 国电南瑞科技股份有限公司 | Wind power penetration optimization evaluation method based on SCED |
| CN104578176A (en) * | 2014-12-11 | 2015-04-29 | 国电南瑞科技股份有限公司 | Method for making power generation plan in consideration of direct current interaction |
-
2016
- 2016-06-24 CN CN201610473659.8A patent/CN106127342B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102280878A (en) * | 2011-07-26 | 2011-12-14 | 国电南瑞科技股份有限公司 | Wind power penetration optimization evaluation method based on SCED |
| CN104578176A (en) * | 2014-12-11 | 2015-04-29 | 国电南瑞科技股份有限公司 | Method for making power generation plan in consideration of direct current interaction |
Non-Patent Citations (4)
| Title |
|---|
| 以直流联络线运行方式优化提升新能源消纳能力的新模式;钟海旺等;《电力系统自动化》;20150210;第39卷(第3期);第36-41页 * |
| 直流跨区互联电网发输电计划模型与方法;王斌等;《电力系统自动化》;20160210;第40卷(第3期);第8-13页 * |
| 考虑跨区直流调峰的日前发电计划优化方法及分析;韩红卫等;《电力系统自动化》;20150825;第39卷(第16期);第138-142页 * |
| 跨区电网中风电消纳影响因素分析及综合评估方法研究;牛东晓;《电网技术》;20160430;第40卷(第4期);第1087-1092页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106127342A (en) | 2016-11-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106127342B (en) | Cross-regional tie line transaction optimization method based on SCUC | |
| Zhang et al. | Optimal operation of integrated electricity and heat system: A review of modeling and solution methods | |
| Wu et al. | Bi-level optimization model for integrated energy system considering the thermal comfort of heat customers | |
| Majidi et al. | Optimal robust operation of combined heat and power systems with demand response programs | |
| Chen et al. | Energy trading and market equilibrium in integrated heat-power distribution systems | |
| Liu et al. | Optimal coordinated operation of a multi-energy community considering interactions between energy storage and conversion devices | |
| Lu et al. | Coordinated dispatch of multi-energy system with district heating network: Modeling and solution strategy | |
| Awad et al. | Impact of energy storage systems on electricity market equilibrium | |
| Farahmand et al. | Balancing market integration in the Northern European continent: A 2030 case study | |
| Oskouei et al. | Scenario-based stochastic optimal operation of wind, photovoltaic, pump-storage hybrid system in frequency-based pricing | |
| Liu et al. | Heat-electricity coupled peak load shifting for multi-energy industrial parks: A Stackelberg game approach | |
| Tan et al. | The optimization model for multi-type customers assisting wind power consumptive considering uncertainty and demand response based on robust stochastic theory | |
| Wang et al. | Transactive energy sharing in a microgrid via an enhanced distributed adaptive robust optimization approach | |
| US20140257584A1 (en) | Energy management system, energy management method, medium, and server | |
| Agamah et al. | A heuristic combinatorial optimization algorithm for load-leveling and peak demand reduction using energy storage systems | |
| JP7181350B2 (en) | Microgrid operation planning device and method, and regional energy management device and energy management device used in microgrid operation planning device | |
| EP2966748A1 (en) | Energy management system, energy management method, program, and server | |
| Hou et al. | Multi-time scale optimization scheduling of microgrid considering source and load uncertainty | |
| CN108596442A (en) | The integrated energy system economic load dispatching method of meter and Conditional Lyapunov ExponentP | |
| Simoglou et al. | Optimal operation of insular electricity grids under high RES penetration | |
| Jones et al. | Evaluation of distributed building thermal energy storage in conjunction with wind and solar electric power generation | |
| Ren et al. | Optimal design and management of distributed energy network considering both efficiency and fairness | |
| CN111667136A (en) | Clearing method and device for regional power market and storage medium | |
| Zhao et al. | Multi-time scale dispatch approach for an AC/DC hybrid distribution system considering the response uncertainty of flexible loads | |
| Perković et al. | Coupling of cleaner production with a day-ahead electricity market: A hypothetical case study |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Address after: Nan Shui Road Gulou District of Nanjing city of Jiangsu Province, No. 8 210003 Applicant after: CHINA ELECTRIC POWER RESEARCH INSTITUTE Co.,Ltd. Applicant after: STATE GRID CORPORATION OF CHINA Applicant after: STATE GRID JIANGSU ELECTRIC POWER Co. Address before: Nan Shui Road Gulou District of Nanjing city of Jiangsu Province, No. 8 210003 Applicant before: China Electric Power Research Institute Applicant before: State Grid Corporation of China Applicant before: STATE GRID JIANGSU ELECTRIC POWER Co. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |