CN107086568B - A kind of optimal support unit localization method of electric system decomposed based on forecast failure - Google Patents

A kind of optimal support unit localization method of electric system decomposed based on forecast failure Download PDF

Info

Publication number
CN107086568B
CN107086568B CN201710347116.6A CN201710347116A CN107086568B CN 107086568 B CN107086568 B CN 107086568B CN 201710347116 A CN201710347116 A CN 201710347116A CN 107086568 B CN107086568 B CN 107086568B
Authority
CN
China
Prior art keywords
bus
forecast failure
constraint
generating set
climbing capacity
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710347116.6A
Other languages
Chinese (zh)
Other versions
CN107086568A (en
Inventor
罗志浩
尹峰
陈波
丁宁
苏烨
丁俊宏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
YINENG ELECTRIC TECHNOLOGY Co Ltd HANGZHOU
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Zhejiang Electric Power Co Ltd
Original Assignee
YINENG ELECTRIC TECHNOLOGY Co Ltd HANGZHOU
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Zhejiang Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by YINENG ELECTRIC TECHNOLOGY Co Ltd HANGZHOU, State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Zhejiang Electric Power Co Ltd filed Critical YINENG ELECTRIC TECHNOLOGY Co Ltd HANGZHOU
Priority to CN201710347116.6A priority Critical patent/CN107086568B/en
Publication of CN107086568A publication Critical patent/CN107086568A/en
Application granted granted Critical
Publication of CN107086568B publication Critical patent/CN107086568B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/001—Methods to deal with contingencies, e.g. abnormalities, faults or failures

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)

Abstract

The invention discloses a kind of optimal support unit localization methods of electric system decomposed based on forecast failure.For electric power grid frequency stabilization, the unit that several pairs of electricity net safety stables play decisive role need to be chosen.The present invention considers power grid difference forecast failure, establishes to minimize power system security constrained optimum tide model of the generating set climbing capacity as target；Use forecast failure decomposition method by the model decomposition for normal operating mode under subproblem under primal problem and each forecast failure method of operation；The solution of above-mentioned security constraint optimal load flow model is realized by the subproblem under the primal problem and each forecast failure method of operation under iterative solution normal operating mode；The minimum generating set climbing capacity revamping workload of gained is solved according to step, the generating set of revamping workload non-zero is subjected to sort descending accordingly, the optimal support unit of decreasing priority is successively denoted as in order, exports as a result.The present invention has preferable applicability, has better meet actual demand.

Description

A kind of optimal support unit localization method of electric system decomposed based on forecast failure
Technical field
The invention belongs to technical field of power systems, specifically a kind of electric system decomposed based on forecast failure is most Excellent support unit localization method.
Background technique
Current power construction is fast-developing, one side various regions power source construction projects still rapid growth, the specific gravity of wind-powered electricity generation, nuclear power It rises year by year；On the other hand a plurality of extra high voltage line is built up, it is transregional come capacitance accounting be continuously improved.While national economy is Structural adjustment is carried out, biggish variation has occurred in society's electricity consumption structure therewith.The increasingly increase of power grid daily load curve, thermoelectricity Unit is declined to a great extent using hourage, when unit operating load reduce when, unit heat consumption is substantially increased, unit efficiency substantially under Drop.Therefore when fired power generating unit is declined to a great extent using hourage, Large-scale fire-electricity unit long-time underrun becomes electricity market It is new under normal conditions, it is necessary to consider the fm capacity of large-scale thermal power machine group.
Under the demand background, thermal power generation unit requires to guarantee the safety of oneself first, while different in major network appearance Ability with certain aid in treatment accident when reason condition.Primary frequency function is exactly that thermal power generation unit is different in power grid appearance The spinning reserve capacity quick response frequency of steam turbine and boiler heat storage and unit is made full use of to change in normal situation, to make up electricity Net power generation and load difference are away from electric power grid frequency stabilization.For this purpose, the machine that several pairs of electricity net safety stables play decisive role need to be chosen Group, the subject of implementation as above-mentioned support function.
Summary of the invention
In consideration of it, it is an object of the present invention to according to the different power grids under the fault condition of extra-high voltage each in power grid pick-up point Response characteristic, the position of the optimal power supply supporting point in each region when providing electric network fault according to region, to determine depth frequency modulation The optimal unit optimum selection supported is the transformation target of frequency modulation function as having height plus going out.
The present invention is realized using following scheme: a kind of optimal support unit positioning of electric system decomposed based on forecast failure Method, it is characterised in that the following steps are included:
Step 1): the steady-state load flow data of electric system, the security constraint data of generator, bus and route, hair are loaded Motor climbing capacity data and power grid forecast failure set data；
Step 2): the Load flow calculation under power grid normal mode obtains initial launch point；
Step 3): it is based on step 1) -2) the data obtained, it establishes to minimize generating set climbing capacity as the electric power of target System security constraint optimal load flow model；
Step 4): based on model obtained by step 3), use forecast failure decomposition method by the model decomposition for normal operation Subproblem under mode under primal problem and each forecast failure method of operation solves the primal problem under normal operating mode；
Step 5): step 4) acquired results are based on, successively solve the subproblem under each forecast failure method of operation, and generate Accordingly additional inequality constraints, is added in primal problem described in step 4)；
Step 6): calculating the increment size of the optimum results of the solved primal problem of step 4) in adjacent two-wheeled iteration, if its Less than given threshold value, then to step 7)；Otherwise return step 5) continue to iteratively solve；
Step 7): according to minimum generating set climbing capacity revamping workload obtained by step 4), by the generator of revamping workload non-zero Group carries out sort descending accordingly, is successively denoted as the optimal support unit of decreasing priority in order, exports as a result.
Further, the particular content of the step 3) are as follows: to minimize generating set climbing capacity revamping workload as optimization Power flow equation after target and the transformation of generating set climbing capacity under power grid normal operating mode and all forecast failure modes For equality constraint, with generator Climing constant, the generator under power grid normal operating mode and all forecast failure methods of operation Active power output constraint, the constraint of generator reactive units limits, node voltage amplitude and route thermostabilization are constrained to inequality constraints, Establish one group of Non-linear Optimal Model.
Further, the Non-linear Optimal Model is specific as follows:
1) objective function
In formula,The respectively revamping workload of the positive and negative climbing capacity of generating set i；To operate normally The active power output of balancing generator group under mode；wi, p is respectively the weight coefficient of climbing capacity revamping workload and network loss；SGFor power generation The set of unit；
2) equality constraint
Equality constraint is the power flow equation under different running method:
In formula, variable subscript represents the method for operation, wherein 0 represents normal operating mode, other represent forecast failure operation Mode；SkFor the set of the forecast failure method of operation；SBFor the set of bus；
UiFor the voltage magnitude of bus i, i ∈ SB；UjFor the voltage magnitude of bus j；PGi, QGiRespectively generating set i's Active and idle power output, i ∈ SG；PLi, QLiThe respectively active and load or burden without work of bus i, i ∈ SG；δijBetween bus ij Phase angle difference；λkFor the load growth factor under method of operation k, for describing the fluctuation of load and its with the growth of national economy；The real and imaginary parts of corresponding element respectively in node admittance matrix；
3) inequality constraints
Inequality constraints includes generating set Climing constant
In formula,The original positive and negative climbing capacity of respectively generating set i；It is issued for normal operating mode The active power output of motor group i,For the active power output of generating set i under method of operation k；
And other operation constraints, including node voltage constraint, generated power and idle units limits and Line Flow Constraint
In formula, FijTrend between bus i to bus j on route；The upper and lower scribing line of each variable is respectively indicated to strain The upper and lower bound of amount.
Further, the particular content of the step 4) are as follows: the Non-linear Optimal Model for being established step 3) is decomposed into 1 The subproblem of the primal problem of a corresponding normal operating mode and corresponding several forecast failure methods of operation.
Further, the particular content of primal problem are as follows:
1) objective function
In formula,The respectively revamping workload of the positive and negative climbing capacity of generating set i；To operate normally The active power output of balancing generator group under mode；wi, p is respectively the weight coefficient of climbing capacity revamping workload and network loss；SGFor power generation The set of unit；
2) equality constraint
Equality constraint is the power flow equation under normal operating mode
In formula, variable subscript 0 represents normal operating mode, SBFor the set of bus；
UiFor the voltage magnitude of bus i, i ∈ SB；UjFor the voltage magnitude of bus j；PGi, QGiRespectively generating set i's Active and idle power output, i ∈ SG；PLi, QLiThe respectively active and load or burden without work of bus i, i ∈ SG；δijBetween bus ij Phase angle difference；λ0For the load growth factor under normal operating mode, for describing the fluctuation of load and its with the increasing of national economy It is long；Gij, BijThe real and imaginary parts of corresponding element respectively in node admittance matrix；
3) inequality constraints
Inequality constraints includes node voltage constraint, generated power and idle units limits, Line Flow constraint:
In formula, FijTrend between bus i to bus j on route；The upper and lower scribing line of each variable is respectively indicated to strain The upper and lower bound of amount.
In the step 5), the particular content of subproblem corresponding to k-th of forecast failure are as follows:
1) objective function
In formula,The positive and negative climbing capacity slack variable of respectively generating set i；
2) equality constraint
Equality constraint is the power flow equation under forecast failure method of operation k
In formula, SBFor the set of bus；UiFor the voltage magnitude of bus i, i ∈ SB；UjFor the voltage magnitude of bus j；PGi, QGiThe active and idle power output of respectively generating set i, i ∈ SG；PLi, QLiThe respectively active and load or burden without work of bus i, i ∈ SG；δijFor the phase angle difference between bus ij；λkFor the load growth factor under method of operation k, for describe load fluctuation and Its with national economy growth； The real and imaginary parts of corresponding element respectively in node admittance matrix；
3) inequality constraints
Inequality constraints includes the unit ramp loss after relaxation
In formula,The original positive and negative climbing capacity of respectively generating set i；
Unit climbing capacity slack variable constraint:
And other operation constraints, including the constraint of forecast failure method of operation k lower node voltage, generated power and idle Units limits and Line Flow constraint:
In formula, FijTrend between bus i to bus j on route；The upper and lower scribing line of each variable is respectively indicated to strain The upper and lower bound of amount.
Further, optimization problem corresponding to above-mentioned forecast failure method of operation k is solved, if optimization aim J non-zero, gives birth to At following inequality constraints, it is added in step 4) in the inequality constraints of primal problem:
In formula, u0,It is asked to be controlled in variable (including bus voltage amplitude etc.) and step 4) under normal operating mode The numerical value for the corresponding control variable that solution primal problem obtains；T indicates transposition；
Changing for the positive and negative climbing capacity of generating set i that primal problem obtains respectively is solved in step 4) The numerical value for the amount of making；
μ,The inequality of the respectively dual variable of equality constraint, the constraint of climbing capacity lower limit and upper limit constraint is about The dual variable of beam.
By implementing above-mentioned steps, the unit for playing decisive role to electricity net safety stable can be accurately positioned, establish phase The priority orders for answering unit located the optimal support unit in the case of power grid different faults, and give unit climbing energy The optimal modification scheme of power, the security feature after electric network fault is met with minimum cost.
Detailed description of the invention
Fig. 1 is flow diagram of the invention.
Specific embodiment
The present invention will be further described with reference to the accompanying drawings and embodiments.
The present embodiment provides a kind of methods of electrical power system transient catastrophe failure screening, as shown in Figure 1, including following step It is rapid:
Step (1): the steady-state load flow data of electric system, the security constraint number of the elements such as generator, bus, route are loaded According to existing generator climbing capacity data and power grid forecast failure set data；
Step (2): the Load flow calculation under power grid normal mode obtains initial launch point；
Step (3): being based on step (1)-(2) the data obtained, establishes to minimize generating set climbing capacity as target Power system security constrained optimum tide model；
Step (4): based on model obtained by step (3), use forecast failure decomposition method by the model decomposition for normal fortune Subproblem under line mode under primal problem and each forecast failure method of operation；Solve the primal problem under normal operating mode；
Step (5): being based on step (4) acquired results, successively solves the subproblem under each forecast failure method of operation, and raw At corresponding additional inequality constraints, it is added in step (4) described primal problem；
Step (6): calculating the increment size of the optimum results of step (4) solved primal problem in adjacent 2 wheel iteration, if its Less than given threshold value, then to step (7)；Otherwise return step (5) continues to iteratively solve；
Step (7): according to minimum generating set climbing capacity revamping workload obtained by step (4), by the power generation of revamping workload non-zero Unit carries out sort descending accordingly, is successively denoted as the optimal support unit of decreasing priority in order, exports as a result.
In the present embodiment, the step (3) is specifically, to minimize generating set climbing capacity revamping workload as optimization mesh Mark, the power flow equation in terms of and after the transformation of generating set climbing capacity under power grid normal operating mode and all forecast failure modes For equality constraint, with generator Climing constant, the generator under power grid normal operating mode and all forecast failure methods of operation The security constraints such as active power output constraint, the constraint of generator reactive units limits, node voltage amplitude, route thermostabilization constraint are not Equality constraint establishes one group of Non-linear Optimal Model, specific as follows:
1) objective function
In formula,The respectively revamping workload of the positive and negative climbing capacity of generating set i；To operate normally The active power output of balancing generator group under mode；wi, p is respectively the weight coefficient of climbing capacity revamping workload and network loss；SGFor power generation The set of unit.
2) equality constraint
Equality constraint is the power flow equation under different running method
In formula, variable subscript represents the method for operation, wherein 0 represents normal operating mode, other represent forecast failure operation Mode；SkFor the set of the forecast failure method of operation；SBFor the set of bus；
UiFor the voltage magnitude of bus i, i ∈ SB；UjFor the voltage magnitude of bus j；PGi, QGiRespectively generating set i's Active and idle power output, i ∈ SG；PLi, QLiThe respectively active and load or burden without work of bus i, i ∈ SG；δijBetween bus ij Phase angle difference；λkFor the load growth factor under method of operation k, for describing the fluctuation of load and its with the growth of national economy；The real and imaginary parts of corresponding element respectively in node admittance matrix.
3) inequality constraints
Inequality constraints includes unit ramp loss
In formula,The original positive and negative climbing capacity of respectively generating set i；It is issued for normal operating mode The active power output of motor group i,For the active power output of generating set i under method of operation k；
And other operation constraints, including node voltage constraint, generated power and idle units limits, Line Flow are about Beam etc.:
In formula, FijTrend between bus i to bus j on route；The upper and lower scribing line of each variable is respectively indicated to strain The upper and lower bound of amount.
In the present embodiment, the step (4) specifically: the Non-linear Optimal Model that step (3) is established is decomposed into 1 The subproblem of the primal problem of a corresponding normal operating mode and corresponding several forecast failure methods of operation.Wherein, primal problem is specific Are as follows:
1) objective function
In formula, variable meaning is same as described above.
2) equality constraint
Equality constraint is the power flow equation under normal operating mode
In formula, variable subscript 0 represents normal operating mode, SBFor the set of bus；
UiFor the voltage magnitude of bus i, i ∈ SB；UjFor the voltage magnitude of bus j；PGi, QGiRespectively generating set i's Active and idle power output, i ∈ SG；PLi, QLiThe respectively active and load or burden without work of bus i, i ∈ SG；δijBetween bus ij Phase angle difference；λ0For the load growth factor under normal operating mode, for describing the fluctuation of load and its with the increasing of national economy It is long；Gij, BijThe real and imaginary parts of corresponding element respectively in node admittance matrix.
3) inequality constraints
Inequality constraints includes node voltage constraint, generated power and idle units limits, Line Flow constraint etc.:
In formula, FijTrend between bus i to bus j on route；The upper and lower scribing line of each variable is respectively indicated to strain The upper and lower bound of amount；
And by a series of step (5) inequality constraints generated.
In the present embodiment, the step (5) specifically: subproblem corresponding to k-th of forecast failure specifically:
1) objective function
In formula,The positive and negative climbing capacity slack variable of respectively generating set i；
2) equality constraint
Equality constraint is the power flow equation under forecast failure method of operation k.
In formula, SBFor the set of bus；UiFor the voltage magnitude of bus i, i ∈ SB；UjFor the voltage magnitude of bus j；PGi, QGiThe active and idle power output of respectively generating set i, i ∈ SG；PLi, QLiThe respectively active and load or burden without work of bus i, i ∈ SG；δijFor the phase angle difference between bus ij；λkFor the load growth factor under method of operation k, for describe load fluctuation and Its with national economy growth； The real and imaginary parts of corresponding element respectively in node admittance matrix.
3) inequality constraints
Inequality constraints includes the unit ramp loss after relaxation
In formula,The original positive and negative climbing capacity of respectively generating set i；
Unit climbing capacity slack variable constraint:
And other operation constraints, including the constraint of forecast failure method of operation k lower node voltage, generated power and idle Units limits, Line Flow constraint etc.:
In formula, FijTrend between bus i to bus j on route；The upper and lower scribing line of each variable is respectively indicated to strain The upper and lower bound of amount.
Optimization problem corresponding to above-mentioned forecast failure method of operation k is solved, if optimization aim J non-zero, is generated as follows not Equality constraint is added in step (4) in the inequality constraints set of primal problem:
In formula, u0,It is obtained for solution primal problem in control variable under normal operating mode and step 4) corresponding Control the numerical value of variable；T indicates transposition；
Changing for the positive and negative climbing capacity of generating set i that primal problem obtains respectively is solved in step 4) The numerical value for the amount of making；
μ,The inequality of the respectively dual variable of equality constraint, the constraint of climbing capacity lower limit and upper limit constraint is about The dual variable of beam.
By implementing above-mentioned steps, the optimal modification scheme of unit climbing capacity is given, and full with minimum cost accordingly Security feature after foot electric network fault.As a result, above-mentioned specific implementation step give it is optimal under power grid different faults Support unit.
The foregoing is merely presently preferred embodiments of the present invention, all equalizations done according to claims of the present invention protection scope Variation and modification, are all covered by the present invention.

Claims (7)

1. a kind of optimal support unit localization method of electric system decomposed based on forecast failure, it is characterised in that including following step It is rapid:
Step 1): the steady-state load flow data of electric system, the security constraint data of generator, bus and route, generator are loaded Climbing capacity data and power grid forecast failure set data；
Step 2): the Load flow calculation under power grid normal mode obtains initial launch point；
Step 3): it is based on step 1) -2) the data obtained, it establishes to minimize generating set climbing capacity as the electric system of target Security constraint optimal load flow model；
Step 4): based on model obtained by step 3), use forecast failure decomposition method by the model decomposition for normal operating mode Subproblem under lower primal problem and each forecast failure method of operation solves the primal problem under normal operating mode；
Step 5): being based on step 4) acquired results, successively solves the subproblem under each forecast failure method of operation, and generates corresponding Additional inequality constraints, is added in primal problem described in step 4)；
Step 6): the increment size of the optimum results of the solved primal problem of step 4) in adjacent two-wheeled iteration is calculated, if it is less than Given threshold value, then to step 7)；Otherwise return step 5) continue to iteratively solve；
Step 7): according to minimum generating set climbing capacity revamping workload obtained by step 4), by the generating set evidence of revamping workload non-zero This carries out sort descending, is successively denoted as the optimal support unit of decreasing priority in order, exports as a result.
2. the electric system optimal support unit localization method according to claim 1 decomposed based on forecast failure, special Sign is,
The particular content of the step 3) are as follows: to minimize generating set climbing capacity revamping workload as optimization aim, and power generation Power flow equation after the transformation of unit climbing capacity under power grid normal operating mode and all forecast failure modes is equality constraint, with Power grid normal operating mode and generator Climing constant under all forecast failure methods of operation, generated power units limits, Generator reactive units limits, node voltage amplitude constraint and route thermostabilization be constrained to inequality constraints, establish one group it is non-thread Property Optimized model.
3. the electric system optimal support unit localization method according to claim 2 decomposed based on forecast failure, special Sign is,
The Non-linear Optimal Model is specific as follows:
1) objective function
In formula,The respectively revamping workload of the positive and negative climbing capacity of generating set i；For normal operating mode The active power output of lower balancing generator group；wi, p is respectively the weight coefficient of climbing capacity revamping workload and network loss；SGFor generating set Set；
2) equality constraint
Equality constraint is the power flow equation under different running method:
In formula, variable subscript represents the method for operation, wherein 0 represents normal operating mode, other represent the forecast failure method of operation； SkFor the set of the forecast failure method of operation；SBFor the set of bus；
UiFor the voltage magnitude of bus i, i ∈ SB；UjFor the voltage magnitude of bus j；PGi, QGiRespectively generating set i active and Idle power output, i ∈ SG；PLi, QLiThe respectively active and load or burden without work of bus i, i ∈ SG；δijFor the phase angle difference between bus ij； λkFor the load growth factor under method of operation k, for describing the fluctuation of load and its with the growth of national economy； The real and imaginary parts of corresponding element respectively in node admittance matrix；
3) inequality constraints
Inequality constraints includes generating set Climing constant
In formula,The original positive and negative climbing capacity of respectively generating set i；For generator under normal operating mode The active power output of group i,For the active power output of generating set i under method of operation k；
And other operation constraints, including node voltage constraint, generated power and idle units limits and Line Flow constraint
In formula, FijTrend between bus i to bus j on route；The upper and lower scribing line of each variable is respectively indicated to dependent variable Upper and lower bound.
4. the electric system optimal support unit localization method according to claim 3 decomposed based on forecast failure, special Sign is,
The particular content of the step 4) are as follows: the Non-linear Optimal Model for being established step 3) is decomposed into 1 correspondence and normally transports The subproblem of the primal problem of line mode and corresponding several forecast failure methods of operation.
5. the electric system optimal support unit localization method according to claim 4 decomposed based on forecast failure, special Sign is,
The particular content of primal problem are as follows:
1) objective function
In formula,The respectively revamping workload of the positive and negative climbing capacity of generating set i；For under normal operating mode The active power output of balancing generator group；wi, p is respectively the weight coefficient of climbing capacity revamping workload and network loss；SGFor generating set Set；
2) equality constraint
Equality constraint is the power flow equation under normal operating mode
In formula, variable subscript 0 represents normal operating mode, SBFor the set of bus；
UiFor the voltage magnitude of bus i, i ∈ SB；UjFor the voltage magnitude of bus j；PGi, QGiRespectively generating set i active and Idle power output, i ∈ SG；PLi, QLiThe respectively active and load or burden without work of bus i, i ∈ SG；δijFor the phase angle difference between bus ij； λ0For the load growth factor under normal operating mode, for describing the fluctuation of load and its with the growth of national economy；Gij, Bij The real and imaginary parts of corresponding element respectively in node admittance matrix；
3) inequality constraints
Inequality constraints includes node voltage constraint, generated power and idle units limits, Line Flow constraint:
In formula, FijTrend between bus i to bus j on route；The upper and lower scribing line of each variable is respectively indicated to dependent variable Upper and lower bound.
6. the electric system optimal support unit localization method according to claim 5 decomposed based on forecast failure, special Sign is,
In the step 5), the particular content of subproblem corresponding to k-th of forecast failure are as follows:
1) objective function
In formula,The positive and negative climbing capacity slack variable of respectively generating set i；
2) equality constraint
Equality constraint is the power flow equation under forecast failure method of operation k
In formula, SBFor the set of bus；UiFor the voltage magnitude of bus i, i ∈ SB；UjFor the voltage magnitude of bus j；PGi, QGiPoint Not Wei generating set i active and idle power output, i ∈ SG；PLi, QLiThe respectively active and load or burden without work of bus i, i ∈ SG；δij For the phase angle difference between bus ij；λkFor the load growth factor under method of operation k, for describing the fluctuation of load and its with state The growth of people's economy； The real and imaginary parts of corresponding element respectively in node admittance matrix；
3) inequality constraints
Inequality constraints includes the unit ramp loss after relaxation
In formula,The original positive and negative climbing capacity of respectively generating set i；
Unit climbing capacity slack variable constraint:
And other operation constraints, including the constraint of forecast failure method of operation k lower node voltage, generated power and idle power output Constraint and Line Flow constraint:
In formula, FijTrend between bus i to bus j on route；The upper and lower scribing line of each variable is respectively indicated to dependent variable Upper and lower bound.
7. the electric system optimal support unit localization method according to claim 6 decomposed based on forecast failure, special Sign is,
Optimization problem corresponding to above-mentioned forecast failure method of operation k is solved, if optimization aim J non-zero, generates such as lower inequality Constraint is added in step 4) in the inequality constraints of primal problem:
In formula, u0,Control variable and the middle solution primal problem of step 4) obtain corresponding respectively under normal operating mode Control the numerical value of variable；T indicates transposition；
The revamping workload for the positive and negative climbing capacity of generating set i that primal problem obtains respectively is solved in step 4) Numerical value；
μ,The inequality constraints of the respectively dual variable of equality constraint, the constraint of climbing capacity lower limit and upper limit constraint Dual variable.
CN201710347116.6A 2017-05-17 2017-05-17 A kind of optimal support unit localization method of electric system decomposed based on forecast failure Active CN107086568B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710347116.6A CN107086568B (en) 2017-05-17 2017-05-17 A kind of optimal support unit localization method of electric system decomposed based on forecast failure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710347116.6A CN107086568B (en) 2017-05-17 2017-05-17 A kind of optimal support unit localization method of electric system decomposed based on forecast failure

Publications (2)

Publication Number Publication Date
CN107086568A CN107086568A (en) 2017-08-22
CN107086568B true CN107086568B (en) 2019-11-12

Family

ID=59609037

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710347116.6A Active CN107086568B (en) 2017-05-17 2017-05-17 A kind of optimal support unit localization method of electric system decomposed based on forecast failure

Country Status (1)

Country Link
CN (1) CN107086568B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108808738B (en) * 2018-05-29 2021-07-13 国电南瑞科技股份有限公司 Power grid safety correction control method considering constraint priority
CN109586275B (en) * 2018-10-18 2020-07-07 清华大学 Quick starting method and device for electromagnetic transient simulation of alternating current-direct current hybrid system
CN110336590A (en) * 2019-07-23 2019-10-15 广东电网有限责任公司 A kind of Fault Locating Method of power telecom network, device and equipment
CN110556824A (en) * 2019-08-23 2019-12-10 广西电网有限责任公司 Power transmission capacity improving method based on bus splitting

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102684224A (en) * 2012-05-25 2012-09-19 浙江大学 Unit combination method for resolving and considering wind power volatility
CN106096751A (en) * 2016-05-15 2016-11-09 国电南瑞科技股份有限公司 Consider that new forms of energy access and participate in Short Term Generation Schedules arrangement and standby Optimal Configuration Method with Demand Side Response

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102684224A (en) * 2012-05-25 2012-09-19 浙江大学 Unit combination method for resolving and considering wind power volatility
CN106096751A (en) * 2016-05-15 2016-11-09 国电南瑞科技股份有限公司 Consider that new forms of energy access and participate in Short Term Generation Schedules arrangement and standby Optimal Configuration Method with Demand Side Response

Also Published As

Publication number Publication date
CN107086568A (en) 2017-08-22

Similar Documents

Publication Publication Date Title
Abazari et al. Coordination strategies of distributed energy resources including FESS, DEG, FC and WTG in load frequency control (LFC) scheme of hybrid isolated micro-grid
Ramírez et al. Placement and sizing of battery energy storage for primary frequency control in an isolated section of the Mexican power system
Ruttledge et al. Emulated inertial response from wind turbines: gain scheduling and resource coordination
CN107086568B (en) A kind of optimal support unit localization method of electric system decomposed based on forecast failure
Kim et al. Supercapacitor to provide ancillary services with control coordination
Cingoz et al. Optimized resource management for PV–fuel-cell-based microgrids using load characterizations
Dhundhara et al. Grid frequency enhancement using coordinated action of wind unit with redox flow battery in a deregulated electricity market
Ding et al. Operation optimization for microgrids under centralized control
Akram et al. Hybrid energy stoarage system for frequency regulation in microgrids with source and load uncertainties
Saxena et al. Enhancement of system performance using STATCOM as dynamic compensator with squirrel cage induction generator (SCIG) based microgrid
CN108390415A (en) A kind of regional power grid new energy consumption capacity calculation methods and system
Varghese et al. Energy storage management of hybrid solar/wind standalone system using adaptive neuro‐fuzzy inference system
Lazaar et al. Modeling and control of a hydrogen-based green data center
Murali et al. Exploration of novel optimal fuzzy‐based controller for enhancement of frequency regulation of deregulated hybrid power system with modified HVDC tie‐line
Obradović et al. Frequency dynamics of the Northern European AC/DC power system: a look-ahead study
Somefun et al. Review of different methods for siting and sizing distributed generator
Chowdhury et al. Optimal placement and sizing of renewable distributed generation in electricity networks considering different load models
Kim et al. Modelling and control coordination scheme of a wind‐to‐hydrogen set for future renewable‐based power systems
Ji et al. Coordinated control and power management of diesel‐PV‐battery in hybrid stand‐alone microgrid system
Wadi et al. Optimization of Distributed Generation Using Homer Software and Fuzzy Logic Control
Illindala Flexible distribution of energy and storage resources
Koutroulis et al. Optimal design and economic evaluation of a battery energy storage system for the maximization of the energy generated by wind farms in isolated electric grids
Anusha et al. Energy Management of an Off-Grid and Grid Connected Hybrid Renewable Energy Source Micro-Grid System for Commercial Load
Ramírez‐González et al. Assessment of inertial and primary frequency control from wind power plants in the Mexican electric power grid
Xu et al. An ABC algorithm for optimization of restoration path in a power grid with HVDC connection

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant