CN107394775A - A kind of active power distribution network sequential islet operation method of meter and load priority - Google Patents

A kind of active power distribution network sequential islet operation method of meter and load priority Download PDF

Info

Publication number
CN107394775A
CN107394775A CN201710772825.9A CN201710772825A CN107394775A CN 107394775 A CN107394775 A CN 107394775A CN 201710772825 A CN201710772825 A CN 201710772825A CN 107394775 A CN107394775 A CN 107394775A
Authority
CN
China
Prior art keywords
mrow
msubsup
msub
node
load
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.)
Pending
Application number
CN201710772825.9A
Other languages
Chinese (zh)
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.)
Tianjin University
State Grid Corp of China SGCC
State Grid Economic and Technological Research Institute
Original Assignee
Tianjin University
State Grid Corp of China SGCC
State Grid Economic and Technological Research Institute
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 Tianjin University, State Grid Corp of China SGCC, State Grid Economic and Technological Research Institute filed Critical Tianjin University
Priority to CN201710772825.9A priority Critical patent/CN107394775A/en
Publication of CN107394775A publication Critical patent/CN107394775A/en
Pending legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/067Enterprise or organisation modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/388Islanding, i.e. disconnection of local power supply from the network

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Physics & Mathematics (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Power Engineering (AREA)
  • Health & Medical Sciences (AREA)
  • Educational Administration (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The present invention relates to a kind of meter and the active power distribution network sequential islet operation method of load priority, its step:According to selected distribution system, distribution system structure and initial parameter value are inputted;According to distribution system structure and initial parameter value, the mathematical modeling of the active power distribution network sequential isolated island partition problem of foundation meter and load priority;The mathematical modeling of active power distribution network sequential isolated island partition problem is solved using prim al- dual interior point m ethod;Solving result is exported, including target function value, each node voltage amplitude, each node recover load coefficient, operation reserve, block switch and the interconnection switch of controlled distribution formula power supply cut-off state.The present invention can establish a kind of meter and the active power distribution network isolated island partitioning model of load priority, consider a variety of safe operations constraints, finally determine optimal isolated island division.

Description

A kind of active power distribution network sequential islet operation method of meter and load priority
Technical field
The present invention relates to a kind of active power distribution network islet operation method, counted and load priority has more particularly to a kind of Source power distribution network sequential islet operation method.
Background technology
Self-healing function construction in intelligent grid by be intelligent grid system construction the ring of key one.Wherein, distribution is utilized Formula power supply realizes the division of electric power isolated island with operating to one of key technology in self-healing control.Power distribution network is direct as power system A user oriented ring, the closest with contacting for user, the influence to user is the most direct.For power distribution network, due to The permeability of distributed power source (distributed generator, DG) is more and more stronger, occurs in transmission system or distribution system , can be by being formed with can because part DG has the ability of actively regulation in the case that failure causes the comprehensive dead electricity of distribution system Control DG recovers the power supply of a part of important load for the power distribution network electric power islet operation of power supply, substantially increases the confession of power distribution network Electric reliability.
A kind of emergency mode that active power distribution network islet operation is handled as electric network fault, system confession can be obviously improved The reliability and security of electricity.Power distribution network is directly connected with user, the target of its islet operation be distribution network failure safety every From rear, the preferential recovery of important load guarantee in dead electricity region, and recover other insignificant loads as far as possible, recover so as to realize Power supply Income Maximum, is a problem urgently to be resolved hurrily for intelligent distribution network.
Generally in large scale containing the grid-connected active power distribution network of a large amount of regenerative resources, connection type is extremely complex, While realizing optimal energy source and distribution, also contain the hidden danger of major accident generation.When major accident occurs, how Optimal electric power isolated island is fast and effeciently determined, realizes Fault Isolation, reduces power failure range as far as possible, play has independently as far as possible The effect of the distributed power source of generating capacity is to improve power supply reliability.
After power distribution network breaks down, the priority of load is considered, be the important load in recovery systems as more as possible, lead to Cross and interconnection switch and block switch are cut-off, change the original topology of network, improve the utilization rate of distributed power source, make system The weighting burden with power amount of recovery is maximum.
The content of the invention
For above-mentioned active power distribution network isolated island partition problem, counted and load priority it is an object of the invention to provide a kind of Active power distribution network sequential islet operation method, it can establish a kind of meter and the active power distribution network isolated island division mould of load priority Type, consider a variety of safe operation constraints, finally determine optimal isolated island division.
To achieve the above object, the present invention takes following technical scheme:A kind of active power distribution network of meter and load priority Sequential islet operation method, it is characterised in that following steps:1) according to selected distribution system, distribution system structure and ginseng are inputted Number initial value;2) according to distribution system structure and initial parameter value, the active power distribution network sequential for establishing meter and load priority is lonely The mathematical modeling of island partition problem;3) mathematical modeling using prim al- dual interior point m ethod to active power distribution network sequential isolated island partition problem Solved;4) solving result is exported, including target function value, each node voltage amplitude, each node recover load coefficient, controllable Operation reserve, block switch and the interconnection switch of distributed power source cut-off state.
Further, in the step 1), initial value includes line parameter circuit value, load level, network topology annexation, distribution System operation voltage level and branch current limitation, controllable and uncontrollable distributed power source on-position, capacity and peak power Factor angle, load and uncontrollable distributed power source operation curve, the priority factor of load, the system failure moment and it is lasting when Between, reference voltage and reference power.
Further, in the step 2), the mathematical modeling of active power distribution network sequential isolated island partition problem, which is specifically established, to be included: The weighting burden with power amount of setting distribution system service restoration is up to object function, considers the radial constraint of system, system tide respectively Stream constraint, working voltage constraint, branch current constraint, uncontrollable distributed power source operation constraint and the operation of controlled distribution formula power supply Constraint.
Further, the weighting burden with power amount maximum target function of the distribution system service restoration is:
In formula, NTFor the time of system islet operation;tFFor system jam at the time of;NNFor all nodes in system Set;For the burden with power in t node i;λiFor the recovery coefficient of load in node i;wiFor in node i Load priority factor.
Further, the system is radial is constrained to:
αijijji,ij∈ΩB
αij∈{0,1}
0≤βij≤1,0≤βji≤1
In formula, ΩBRepresent the set of all branch roads equipped with block switch or interconnection switch;NNFor the collection of system interior joint Close, NSThe set of source node in expression system;αijRepresent branch road ij upper switch cut-offs state, αij=1 represents switch closure, αij=0 represents to switch off;βijRepresent node i and node j relation, βij=1 expression node j be node i father node, βij =0 expression node j is not the father node of node i.
Further, the trend constraint is:
-Mαij≤Pt,ij≤Mαij
-Mαij≤Qt,ij≤Mαij
In formula, rijFor branch road ij resistance, xijFor branch road ij reactance;Pt,ijIt is active to be flowed through on t branch road ij Power, Qt,ijFor the reactive power flowed through on t branch road ij, k ∈ j;Pt,iFor the active power injected in t node i it With, Qt,iFor the reactive power sum injected in t node i;Load respectively in t node i Active power and reactive power;For in t node i controlled distribution formula power supply inject active power,For t section On point i controlled distribution formula power supply inject reactive power,Injected for uncontrollable distributed power source in t node i active Power,For the reactive power of uncontrollable distributed power source injection in t node i;Parameter M expressions one are greatly normal Amount;Ut,iRepresent the voltage magnitude of t node i;λiRepresent the recovery coefficient of load in node i;It,ijRepresent t branch road ij On current amplitude.
Further, the working voltage is constrained to:
In formula,WithThe respectively voltage upper and lower limit of node i.
Further, the branch current is constrained to:
In formula,For branch road ij upper current limit.
Further, the uncontrollable distributed power source operation is constrained to:
In formula,The upper limit of the active power output of uncontrollable distributed power source is connect by t node i;Represent Node i connects the capacity of uncontrollable distributed power source;For the maximum power factor of the operation of distributed power source in node i Angle.
Further, the controlled distribution formula power supply operation is constrained to:
In formula,Represent that node i connects the capacity of controlled distribution formula power supply.
For the present invention due to taking above technical scheme, it has advantages below:1st, the present invention is based on the active distribution of solution Isolated island partition problem in net, in the case where considering network reconfiguration and load priority, take into full account that distributed power source is contributed The influence of moment and trouble duration to active power distribution network isolated island partition problem occurs with load timing, and failure, builds The active power distribution network sequential isolated island partitioning model of a kind of vertical meter and load priority, its mathematics essence is non-linear programming problem, Solved using prim al- dual interior point m ethod, obtain service restoration scheme.2nd, the present invention is on the basis of network reconfiguration is considered, simultaneously Consider that distributed power source is contributed and the sequential operation characteristic of load, ensures that important load is held in active power distribution network to greatest extent Continuous stable operation.3rd, it is different at the time of the present invention breaks down according to power distribution network, and trouble duration difference, obtain difference Isolated island partition strategy, to ensure that it is maximum that system recovers load.
Brief description of the drawings
Fig. 1 is overall flow schematic diagram of the present invention;
Fig. 2 is the improved node example structure charts of IEEE 33;
Fig. 3 is load characteristic curve;
Fig. 4 is photovoltaic power curve;
Fig. 5 is the system jam at 0, the schematic diagram of the isolated island division result of 4 hours fault times;
Fig. 6 is the system jam at 8, the schematic diagram of the isolated island division result of 4 hours fault times.
Embodiment
The present invention is described in detail with reference to the accompanying drawings and examples.
As shown in figure 1, the present invention provides a kind of active power distribution network sequential islet operation method of meter and load priority, its Comprise the following steps:
1) according to selected distribution system, distribution system structure and initial parameter value are inputted;
Wherein, initial value includes line parameter circuit value, load level, network topology annexation, distribution system working voltage water The limitation of gentle branch current, controllable and uncontrollable distributed power source on-position, capacity and maximum power factor angle, load and not Controlled distribution formula power supply operation curve, the priority factor of load, system failure moment and duration, reference voltage and benchmark Power etc..
In the present embodiment, first input the node systems of IEEE 33 in circuit element impedance value, load cell it is active Power, reactive power, load operation curve, distributed electrical source dates, distributed power source operation curve, the priority system of load Number, network topology annexation, example structure as shown in Fig. 2 detail parameters as shown in table 1, table 2, table 3, table 4, Fig. 3 and Fig. 4; Permanent three-phase fault occurs respectively when 0 and 8 between setting branch road 1-2, fault time continues 4 hours;The base of initialization system Quasi- voltage is 12.66kV, reference power 1MVA.
The IEEE33 nodes example load on-position of table 1 and power
The IEEE33 node example line parameter circuit values of table 2
The distributed power source configuring condition of table 3
The load priority of table 4
2) the distribution system structure and initial parameter value provided according to step 1), establishes meter and the active of load priority is matched somebody with somebody The mathematical modeling of power network sequential isolated island partition problem;
The mathematical modeling of active power distribution network sequential isolated island partition problem, which is specifically established, to be included:Set adding for distribution system service restoration Power burden with power amount be up to object function, respectively consider system it is radial constraint, system load flow constraint, working voltage constraint, Branch current constraint, uncontrollable distributed power source operation constraint and the operation constraint of controlled distribution formula power supply;Specific method for building up is such as Under:
2.1) the weighting burden with power amount maximum target function of distribution system service restoration is:
In formula, NTFor the time of system islet operation;tFFor system jam at the time of;NNFor all nodes in system Set;For the burden with power in t node i;λiIt is 0-1 variables for the recovery coefficient of load in node i, λi= 1 represents that the node load recovers, λi=0 represents that the node load does not recover;wiFor the load priority factor in node i.
2.2) system is radial is constrained to:
αijijji,ij∈ΩB (2)
αij∈{0,1} (5)
0≤βij≤1,0≤βji≤1 (6)
In formula, ΩBRepresent the set of all branch roads equipped with block switch or interconnection switch;NNFor the collection of system interior joint Close, NSThe set of source node in expression system;αijRepresent branch road ij upper switch cut-offs state, αij=1 represents switch closure, αij=0 represents to switch off;βijRepresent node i and node j relation, βij=1 expression node j be node i father node, βij =0 expression node j is not the father node of node i.
Network reconfiguration not only can be carried out to system by the radial constraint of system, can also be to controlled distribution formula power supply Operation method is selected, when the controlled distribution formula power supply in node i meets formula ∑ij∈ΩBβijWhen=1, distributed power source choosing Take PQ control modes;When the controlled distribution formula power supply in node i meets formula ∑ij∈ΩBβijWhen=0, the distributed power source chooses V/f Control mode.
2.3) trend constraint is:
-Mαij≤Pt,ij≤Mαij (14)
-Mαij≤Qt,ij≤Mαij (15)
In formula, rijFor branch road ij resistance, xijFor branch road ij reactance;Pt,ijIt is active to be flowed through on t branch road ij Power, Qt,ijFor the reactive power flowed through on t branch road ij, k ∈ j;Pt,iFor the active power injected in t node i it With, Qt,iFor the reactive power sum injected in t node i;Load respectively in t node i Active power and reactive power;For in t node i controlled distribution formula power supply inject active power,For t section On point i controlled distribution formula power supply inject reactive power,Injected for uncontrollable distributed power source in t node i active Power,For the reactive power of uncontrollable distributed power source injection in t node i;Parameter M expressions one are greatly normal Amount, generally takes 9999;Ut,iRepresent the voltage magnitude of t node i;λiRepresent the recovery coefficient of load in node i;It,ijRepresent Current amplitude on t branch road ij.
2.4) working voltage is constrained to:
In formula,WithThe respectively voltage upper and lower limit of node i.
2.5) branch current is constrained to:
In formula,For branch road ij upper current limit.
2.6) uncontrollable distributed power source operation is constrained to:
In formula,The upper limit of the active power output of uncontrollable distributed power source is connect by t node i;Represent Node i connects the capacity of uncontrollable distributed power source;For the maximum power factor of the operation of distributed power source in node i Angle.
2.7) controlled distribution formula power supply operation is constrained to:
In formula,Represent that node i connects the capacity of controlled distribution formula power supply.
3) mathematical modeling of active power distribution network sequential isolated island partition problem is solved using prim al- dual interior point m ethod;
4) solving result of step 3) is exported, including target function value, each node voltage amplitude, each node recover load system Number, operation reserve, block switch and the interconnection switch of controlled distribution formula power supply cut-off state etc..
Embodiment:The meter of the present invention and the active power distribution network sequential islet operation method of load priority, using interior point method Solved, obtain optimal islet operation strategy, controlled distribution formula power supply operation method is shown in Table 5, recovers load condition and is shown in Table 6。
The operation method of the controlled distribution formula power supply of table 5
Node 13 20 23 30
When 0 V/f is controlled PQ is controlled V/f is controlled PQ is controlled
When 8 PQ is controlled PQ is controlled V/f is controlled PQ is controlled
The load restoration situation of table 6
It is Intel (R) Core (TM) i5-3470CPU to perform the computer hardware environment that optimization calculates, and dominant frequency is 3.20GHz, inside save as 4GB;Software environment is the operating systems of Windows 7.
Distributed power source is considered to contribute and load timing, to greatest extent in guarantee system important load continue it is steady Fixed operation.Load in system is divided into three classes, i.e. 1 type load, 2 type loads and 3 type loads, weight takes 100,10 and 1 respectively. It is different at the time of being broken down according to power distribution network, and trouble duration difference, network reconfiguration is carried out to system, obtains difference Isolated island partition strategy, to ensure that the active income amount that system obtains is maximum.When with the system jam moment being respectively 0 and 8 When, the system islet operation time is 4 hours, optimal isolated island result as shown in Figure 5 and Figure 6, wherein solid node represent load Recover, hollow node represents that load does not recover.
The various embodiments described above are merely to illustrate the present invention, and each step can be all varied from, in the technology of the present invention On the basis of scheme, all improvement carried out according to the principle of the invention to separate step and equivalents, it should not exclude in this hair Outside bright protection domain.

Claims (10)

1. a kind of active power distribution network sequential islet operation method of meter and load priority, it is characterised in that following steps:
1) according to selected distribution system, distribution system structure and initial parameter value are inputted;
2) according to distribution system structure and initial parameter value, the active power distribution network sequential isolated island for establishing meter and load priority divides The mathematical modeling of problem;
3) mathematical modeling of active power distribution network sequential isolated island partition problem is solved using prim al- dual interior point m ethod;
4) solving result is exported, including target function value, each node voltage amplitude, each node recover load coefficient, controlled distribution Operation reserve, block switch and the interconnection switch of formula power supply cut-off state.
2. the active power distribution network sequential islet operation method of a kind of meter as claimed in claim 1 and load priority, its feature It is:In the step 1), initial value includes line parameter circuit value, load level, network topology annexation, distribution system operation electricity Voltage levels and branch current limitation, controllable and uncontrollable distributed power source on-position, capacity and maximum power factor angle, load With uncontrollable distributed power source operation curve, the priority factor of load, system failure moment and duration, reference voltage and Reference power.
3. the active power distribution network sequential islet operation method of a kind of meter as claimed in claim 1 and load priority, its feature It is:In the step 2), the mathematical modeling of active power distribution network sequential isolated island partition problem, which is specifically established, to be included:Set power distribution system The weighting burden with power amount that system recovers is up to object function, considers the radial constraint of system, system load flow constraint, operation respectively Voltage constraint, branch current constraint, uncontrollable distributed power source operation constraint and the operation constraint of controlled distribution formula power supply.
4. the active power distribution network sequential islet operation method of a kind of meter as claimed in claim 3 and load priority, its feature It is:The weighting burden with power amount maximum target function of the distribution system service restoration is:
<mrow> <mi>f</mi> <mo>=</mo> <msubsup> <mi>max&amp;Sigma;</mi> <mrow> <mi>t</mi> <mo>=</mo> <msup> <mi>t</mi> <mi>F</mi> </msup> </mrow> <msub> <mi>N</mi> <mi>T</mi> </msub> </msubsup> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <msub> <mi>N</mi> <mi>N</mi> </msub> </msubsup> <msub> <mi>&amp;lambda;</mi> <mi>i</mi> </msub> <msub> <mi>w</mi> <mi>i</mi> </msub> <msubsup> <mi>P</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mrow> <mi>L</mi> <mi>O</mi> <mi>A</mi> <mi>D</mi> </mrow> </msubsup> </mrow>
In formula, NTFor the time of system islet operation;tFFor system jam at the time of;NNFor the collection of all nodes in system Close;For the burden with power in t node i;λiFor the recovery coefficient of load in node i;wiFor the load in node i Priority factor.
5. the active power distribution network sequential islet operation method of a kind of meter as claimed in claim 4 and load priority, its feature It is:The system is radial to be constrained to:
αijijji,ij∈ΩB
<mrow> <msub> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mi>j</mi> <mo>&amp;Element;</mo> <msub> <mi>&amp;Omega;</mi> <mi>B</mi> </msub> </mrow> </msub> <msub> <mi>&amp;beta;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>i</mi> <mo>&amp;Element;</mo> <msub> <mi>N</mi> <mi>N</mi> </msub> <mo>\</mo> <msub> <mi>N</mi> <mi>S</mi> </msub> </mrow>
<mrow> <msub> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mi>j</mi> <mo>&amp;Element;</mo> <msub> <mi>&amp;Omega;</mi> <mi>B</mi> </msub> </mrow> </msub> <msub> <mi>&amp;beta;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>i</mi> <mo>&amp;Element;</mo> <msub> <mi>N</mi> <mi>S</mi> </msub> </mrow>
αij∈{0,1}
0≤βij≤1,0≤βji≤1
In formula, ΩBRepresent the set of all branch roads equipped with block switch or interconnection switch;NNFor the set of system interior joint, NS The set of source node in expression system;αijRepresent branch road ij upper switch cut-offs state, αij=1 represents switch closure, αij=0 Expression switches off;βijRepresent node i and node j relation, βij=1 expression node j be node i father node, βij=0 table It is not the father node of node i to show node j.
A kind of 6. active power distribution network sequential islet operation side of meter and load priority as described in any one of claim 3 to 5 Method, it is characterised in that:The trend constraint is:
<mrow> <msub> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mi>k</mi> <mo>&amp;Element;</mo> <msub> <mi>&amp;Omega;</mi> <mi>B</mi> </msub> </mrow> </msub> <msub> <mi>P</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> <mi>k</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>&amp;Sigma;</mi> <mrow> <mi>j</mi> <mi>i</mi> <mo>&amp;Element;</mo> <msub> <mi>&amp;Omega;</mi> <mi>B</mi> </msub> </mrow> </msub> <mrow> <mo>(</mo> <msub> <mi>P</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>j</mi> <mi>i</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>r</mi> <mrow> <mi>j</mi> <mi>i</mi> </mrow> </msub> <msubsup> <mi>I</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>j</mi> <mi>i</mi> </mrow> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>P</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> </mrow>
<mrow> <msub> <mi>P</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>=</mo> <msubsup> <mi>P</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>C</mi> </msubsup> <mo>+</mo> <msubsup> <mi>P</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mrow> <mi>D</mi> <mi>G</mi> </mrow> </msubsup> <mo>-</mo> <msub> <mi>&amp;lambda;</mi> <mi>i</mi> </msub> <msubsup> <mi>P</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mrow> <mi>L</mi> <mi>O</mi> <mi>A</mi> <mi>D</mi> </mrow> </msubsup> </mrow>
<mrow> <msub> <mi>Q</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>=</mo> <msubsup> <mi>Q</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>C</mi> </msubsup> <mo>+</mo> <msubsup> <mi>Q</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mrow> <mi>D</mi> <mi>G</mi> </mrow> </msubsup> <mo>-</mo> <msub> <mi>&amp;lambda;</mi> <mi>i</mi> </msub> <msubsup> <mi>Q</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mrow> <mi>L</mi> <mi>O</mi> <mi>A</mi> <mi>D</mi> </mrow> </msubsup> </mrow>
<mrow> <msubsup> <mi>I</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> <mi>j</mi> </mrow> <mn>2</mn> </msubsup> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>P</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> <mi>j</mi> </mrow> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>Q</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> <mi>j</mi> </mrow> <mn>2</mn> </msubsup> </mrow> <msubsup> <mi>U</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mn>2</mn> </msubsup> </mfrac> </mrow>
<mrow> <msubsup> <mi>U</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>U</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>j</mi> </mrow> <mn>2</mn> </msubsup> <mo>-</mo> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <msub> <mi>P</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <msub> <mi>Q</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <msubsup> <mi>r</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> <msubsup> <mi>I</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> <mi>j</mi> </mrow> <mn>2</mn> </msubsup> <mo>+</mo> <mi>M</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>&amp;GreaterEqual;</mo> <mn>0</mn> </mrow>
<mrow> <msubsup> <mi>U</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>U</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>j</mi> </mrow> <mn>2</mn> </msubsup> <mo>-</mo> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <msub> <mi>P</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>x</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <msub> <mi>Q</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <msubsup> <mi>r</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> <msubsup> <mi>I</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> <mi>j</mi> </mrow> <mn>2</mn> </msubsup> <mo>-</mo> <mi>M</mi> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>&amp;le;</mo> <mn>0</mn> </mrow>
-Mαij≤Pt,ij≤Mαij
-Mαij≤Qt,ij≤Mαij
<mrow> <mn>0</mn> <mo>&amp;le;</mo> <msubsup> <mi>I</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> <mi>j</mi> </mrow> <mn>2</mn> </msubsup> <mo>&amp;le;</mo> <msub> <mi>M&amp;alpha;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mrow>
In formula, rijFor branch road ij resistance, xijFor branch road ij reactance;Pt,ijFor flowed through on t branch road ij active power, Qt,ijFor the reactive power flowed through on t branch road ij, k ∈ j;Pt,iFor injected in t node i active power sum, Qt,iFor the reactive power sum injected in t node i;Load is active respectively in t node i Power and reactive power;For in t node i controlled distribution formula power supply inject active power,For t node i The reactive power of upper controlled distribution formula power supply injection,For the wattful power of uncontrollable distributed power source injection in t node i Rate,For the reactive power of uncontrollable distributed power source injection in t node i;Parameter M represents a great constant; Ut,iRepresent the voltage magnitude of t node i;λiRepresent the recovery coefficient of load in node i;It,ijRepresent on t branch road ij Current amplitude.
7. the active power distribution network sequential islet operation method of a kind of meter as claimed in claim 6 and load priority, its feature It is:The working voltage is constrained to:
<mrow> <msubsup> <mi>U</mi> <mi>i</mi> <mi>min</mi> </msubsup> <mo>&amp;le;</mo> <msub> <mi>U</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>&amp;le;</mo> <msubsup> <mi>U</mi> <mi>i</mi> <mi>max</mi> </msubsup> </mrow>
In formula,WithThe respectively voltage upper and lower limit of node i.
8. the active power distribution network sequential islet operation method of a kind of meter as claimed in claim 7 and load priority, its feature It is:The branch current is constrained to:
<mrow> <mo>-</mo> <msubsup> <mi>I</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>max</mi> </msubsup> <mo>&amp;le;</mo> <msub> <mi>I</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>&amp;le;</mo> <msubsup> <mi>I</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>max</mi> </msubsup> </mrow>
In formula,For branch road ij upper current limit.
9. the active power distribution network sequential islet operation method of a kind of meter as claimed in claim 8 and load priority, its feature It is:The uncontrollable distributed power source operation is constrained to:
<mrow> <msubsup> <mi>P</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mrow> <mi>D</mi> <mi>G</mi> </mrow> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mrow> <mi>D</mi> <mi>G</mi> <mo>,</mo> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msubsup> </mrow>
<mrow> <mo>|</mo> <msubsup> <mi>Q</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mrow> <mi>D</mi> <mi>G</mi> </mrow> </msubsup> <mo>|</mo> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mrow> <mi>D</mi> <mi>G</mi> </mrow> </msubsup> <mi>t</mi> <mi>a</mi> <mi>n</mi> <mrow> <mo>(</mo> <msup> <mi>cos</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msubsup> <mi>&amp;theta;</mi> <mi>i</mi> <mi>max</mi> </msubsup> <mo>)</mo> </mrow> </mrow>
<mrow> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msubsup> <mi>P</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mrow> <mi>D</mi> <mi>G</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mi>Q</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mrow> <mi>D</mi> <mi>G</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>&amp;le;</mo> <msubsup> <mi>S</mi> <mi>i</mi> <mrow> <mi>D</mi> <mi>G</mi> </mrow> </msubsup> </mrow>
In formula,The upper limit of the active power output of uncontrollable distributed power source is connect by t node i;Represent node i Connect the capacity of uncontrollable distributed power source;For the maximum power factor angle of the operation of distributed power source in node i.
10. the active power distribution network sequential islet operation method of a kind of meter as claimed in claim 9 and load priority, its feature It is:The controlled distribution formula power supply operation is constrained to:
<mrow> <mo>|</mo> <msubsup> <mi>Q</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>C</mi> </msubsup> <mo>|</mo> <mo>&amp;le;</mo> <msubsup> <mi>P</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>C</mi> </msubsup> <mi>t</mi> <mi>a</mi> <mi>n</mi> <mrow> <mo>(</mo> <msup> <mi>cos</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msubsup> <mi>&amp;theta;</mi> <mi>i</mi> <mi>max</mi> </msubsup> <mo>)</mo> </mrow> </mrow>
<mrow> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msubsup> <mi>P</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>C</mi> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mi>Q</mi> <mrow> <mi>t</mi> <mo>,</mo> <mi>i</mi> </mrow> <mi>C</mi> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>&amp;le;</mo> <msubsup> <mi>S</mi> <mi>i</mi> <mi>C</mi> </msubsup> </mrow>
In formula,Represent that node i connects the capacity of controlled distribution formula power supply.
CN201710772825.9A 2017-08-31 2017-08-31 A kind of active power distribution network sequential islet operation method of meter and load priority Pending CN107394775A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710772825.9A CN107394775A (en) 2017-08-31 2017-08-31 A kind of active power distribution network sequential islet operation method of meter and load priority

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710772825.9A CN107394775A (en) 2017-08-31 2017-08-31 A kind of active power distribution network sequential islet operation method of meter and load priority

Publications (1)

Publication Number Publication Date
CN107394775A true CN107394775A (en) 2017-11-24

Family

ID=60349048

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710772825.9A Pending CN107394775A (en) 2017-08-31 2017-08-31 A kind of active power distribution network sequential islet operation method of meter and load priority

Country Status (1)

Country Link
CN (1) CN107394775A (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109038560A (en) * 2018-08-03 2018-12-18 国家电网有限公司 Power distribution network distributed energy storage Economic Analysis Method and system based on operation reserve
CN110112736A (en) * 2019-06-06 2019-08-09 南京工程学院 A kind of power distribution network self-healing method based on mixed integer linear programming
CN110263981A (en) * 2019-05-30 2019-09-20 天津大学 Consider that the gas-of flexible scheduling strategy is electrically coupled integrated energy system planing method
CN110417007A (en) * 2019-07-25 2019-11-05 天津大学 A kind of timing quantitative analysis method of active distribution system operational flexibility
CN110492526A (en) * 2019-08-29 2019-11-22 国网上海市电力公司 A kind of isolated island division methods based on deferrable load
CN113129167A (en) * 2021-05-17 2021-07-16 国网河北省电力有限公司电力科学研究院 Method and device for recovering power supply of transformer area and terminal equipment

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104934964A (en) * 2015-03-18 2015-09-23 华南理工大学 Power distribution network reconstruction and island division method containing distributed power supply
CN105023058A (en) * 2015-07-07 2015-11-04 天津大学 Power distribution network intelligent soft switch operation optimization method with simultaneous consideration of switch motion
CN105977934A (en) * 2016-06-21 2016-09-28 天津大学 Method for recovering power supply of soft open point of power distribution network in consideration of load importance

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104934964A (en) * 2015-03-18 2015-09-23 华南理工大学 Power distribution network reconstruction and island division method containing distributed power supply
CN105023058A (en) * 2015-07-07 2015-11-04 天津大学 Power distribution network intelligent soft switch operation optimization method with simultaneous consideration of switch motion
CN105977934A (en) * 2016-06-21 2016-09-28 天津大学 Method for recovering power supply of soft open point of power distribution network in consideration of load importance

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109038560A (en) * 2018-08-03 2018-12-18 国家电网有限公司 Power distribution network distributed energy storage Economic Analysis Method and system based on operation reserve
CN109038560B (en) * 2018-08-03 2020-08-28 国家电网有限公司 Power distribution network distributed energy storage economy evaluation method and system based on operation strategy
CN110263981A (en) * 2019-05-30 2019-09-20 天津大学 Consider that the gas-of flexible scheduling strategy is electrically coupled integrated energy system planing method
CN110263981B (en) * 2019-05-30 2023-04-07 天津大学 Gas-electric coupling comprehensive energy system planning method considering flexible scheduling strategy
CN110112736A (en) * 2019-06-06 2019-08-09 南京工程学院 A kind of power distribution network self-healing method based on mixed integer linear programming
CN110417007A (en) * 2019-07-25 2019-11-05 天津大学 A kind of timing quantitative analysis method of active distribution system operational flexibility
CN110417007B (en) * 2019-07-25 2023-03-10 天津大学 Time sequence quantitative analysis method for operation flexibility of active power distribution system
CN110492526A (en) * 2019-08-29 2019-11-22 国网上海市电力公司 A kind of isolated island division methods based on deferrable load
CN110492526B (en) * 2019-08-29 2023-04-07 国网上海市电力公司 Island division method based on adjustable load
CN113129167A (en) * 2021-05-17 2021-07-16 国网河北省电力有限公司电力科学研究院 Method and device for recovering power supply of transformer area and terminal equipment

Similar Documents

Publication Publication Date Title
CN107394775A (en) A kind of active power distribution network sequential islet operation method of meter and load priority
CN110350508B (en) Method for constructing unified fault recovery model of active power distribution network
Luo et al. Optimal sizing and control strategy of isolated grid with wind power and energy storage system
Maknouninejad et al. Realizing unified microgrid voltage profile and loss minimization: A cooperative distributed optimization and control approach
CN105119280B (en) Alternating current-direct current mixed structure active power distribution network running optimizatin method based on cone optimization
CN108988322A (en) The microgrid optimization of operation strategy method of consideration system time variation
CN105186500B (en) A kind of power distribution network power dissipation coordination optimizing method based on weighting acceleration Lagrangian again
CN106026169B (en) A kind of composition decomposition optimization method that power distribution network is incorporated to based on more micro-capacitance sensors
CN107341623A (en) A kind of active power distribution network source storage lotus islet operation method of meter and network reconfiguration
CN103904644B (en) A kind of Automatic load distribution method for intelligent transformer substation accessed based on distributed power source
CN104300538B (en) Plant stand coordinates update method and the system of coordination constraint limit value in automatism voltage control
CN105790255B (en) Multi-microgrid power supply system and control method
CN107221930A (en) A kind of intelligent Sofe Switch service restoration method of active power distribution network
Debbarma et al. Grid frequency support from V2G aggregators and HVdc links in presence of nonsynchronous units
Mohammadi et al. Adaptive voltage-based load shedding scheme for the DC microgrid
CN107451743A (en) A kind of active power distribution network islet operation method of consideration source lotus operation characteristic
CN107332290B (en) Regional load transfer method based on direct current circuit
Dozein et al. Frequency response capabilities of utility-scale battery energy storage systems, with application to the august 2018 separation event in australia
CN107332277A (en) Store up the active power distribution network islet operation method of lotus operation characteristic in a kind of consideration source
CN104124700A (en) Method and system for generating power distribution network black start scheme
CN112598159A (en) Island fusion recovery method based on multi-source network reconstruction
CN108206540A (en) Consider the active power distribution network isolated island division methods in important load continuous working period
Omari et al. An online control strategy for DC coupled hybrid power systems
CN108599974B (en) Micro-grid distributed communication topology design method based on graph theory connectivity
Valarmathi et al. Simulation and analysis of wind energy and photo voltaic hybrid system

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20171124