CN114552669B - Flexibility-considered partitioning method for distributed power supply distribution network containing high permeability - Google Patents
Flexibility-considered partitioning method for distributed power supply distribution network containing high permeability Download PDFInfo
- Publication number
- CN114552669B CN114552669B CN202210193972.1A CN202210193972A CN114552669B CN 114552669 B CN114552669 B CN 114552669B CN 202210193972 A CN202210193972 A CN 202210193972A CN 114552669 B CN114552669 B CN 114552669B
- Authority
- CN
- China
- Prior art keywords
- cluster
- index
- node
- distribution network
- flexibility
- 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
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000000638 solvent extraction Methods 0.000 title claims abstract description 18
- 230000035699 permeability Effects 0.000 title abstract description 8
- 238000005192 partition Methods 0.000 claims abstract description 32
- 230000006978 adaptation Effects 0.000 claims abstract description 30
- 230000005540 biological transmission Effects 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 21
- 238000004146 energy storage Methods 0.000 claims description 17
- 230000008859 change Effects 0.000 claims description 9
- 230000004044 response Effects 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 101000912561 Bos taurus Fibrinogen gamma-B chain Proteins 0.000 claims description 3
- 238000007599 discharging Methods 0.000 claims description 3
- 230000002222 downregulating effect Effects 0.000 claims description 3
- 238000009434 installation Methods 0.000 claims description 3
- 238000013139 quantization Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 2
- QCAWEPFNJXQPAN-UHFFFAOYSA-N methoxyfenozide Chemical compound COC1=CC=CC(C(=O)NN(C(=O)C=2C=C(C)C=C(C)C=2)C(C)(C)C)=C1C QCAWEPFNJXQPAN-UHFFFAOYSA-N 0.000 claims 1
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 201000006549 dyspepsia Diseases 0.000 abstract description 2
- 230000009194 climbing Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
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
- 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
- H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- 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/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
- H02J3/14—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by switching loads on to, or off from, network, e.g. progressively balanced loading
- H02J3/144—Demand-response operation of the power transmission or distribution network
-
- 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/28—Arrangements for balancing of the load in a network by storage of energy
- H02J3/32—Arrangements for balancing of the load in a network by storage of energy using batteries with converting means
-
- 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/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/22—The renewable source being solar energy
- H02J2300/24—The renewable source being solar energy of photovoltaic origin
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/28—The renewable source being wind energy
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E70/00—Other energy conversion or management systems reducing GHG emissions
- Y02E70/30—Systems combining energy storage with energy generation of non-fossil origin
Abstract
The invention discloses a power distribution network partitioning method with high permeability distributed power supply considering flexibility, comprising the following steps: 1. according to the regulation characteristics of existing flexible resources in the power distribution network and the fluctuation of original DGs and loads, respectively establishing a corresponding cluster flexibility supply model and a cluster flexibility demand model; 2. a payload adaptation rate index describing the flexibility supply in the cluster relative to the flexibility demand matching degree is proposed; 3. providing a branch load margin index for describing the spatial transmission characteristic of flexible resource output in a cluster; 4. providing a modularity value index describing cluster structure characteristics; 5. and giving the index to a certain weight coefficient to form a comprehensive index of the network partition. The invention can improve the autonomous characteristic of the cluster, and exert the adjustment capability of flexible resources in the cluster to the greatest extent, thereby being beneficial to solving the problem of difficult digestion caused by insufficient operation flexibility in the later stage of power distribution network planning.
Description
Technical Field
The invention relates to the field of high-permeability distributed power supply access distribution network planning, in particular to a method for partitioning a distribution network of a distributed power supply with high permeability, which takes flexibility into consideration.
Background
The large amount of penetration of intermittent distributed power sources such as wind power and photovoltaic, and the increase of novel loads with uncertain empty height such as electric automobiles, provide great challenges for the operation of a power distribution network, and cause problems such as out-of-limit voltage fluctuation, power flow pouring, increased system loss, reduced DG consumption level, unbalanced power among feeder lines and the like. Meanwhile, the distributed power supply is largely connected, and the characteristics of small single-machine capacity, large number and scattered geographic positions greatly increase the complexity of the power grid, so that the traditional centralized management structure of the power distribution network is difficult to meet the requirement of controlling time scale in the operation stage, and the power distribution network has the problem of serious insufficient flexibility in operation; on the other hand, the power supply planning problem and the operation control problem of the power distribution network are mutually influenced, so that the operation management of the power distribution network after the large-scale distributed power supply is accessed is necessary to be performed in a partition mode.
The power distribution network is partitioned by taking the clusters as basic units, and in a power system, the clusters have the advantages of coupling and cooperation of nodes in the clusters, and loose and labor division among the clusters. The application of the cluster in the power system mainly comprises two fields of dispatching control and power grid planning, and most of work is concentrated on dispatching control at present, wherein the fields comprise reactive voltage control, power grid partition, active power control and the like. The method specifically relates to the following partitioning method: the space electric distance is used as an index, and an immunity-central point clustering algorithm is adopted to carry out reactive voltage control partition on the power distribution network system; guiding the division of the power distribution network based on cluster modularity performance indexes of the electric distance and the regional voltage regulation capability; cluster division is carried out based on the modularity index of the electrical distance so as to facilitate reasonable calling of energy storage regulation resources, and partition control of the voltage of the power distribution network is realized; and constructing comprehensive performance indexes based on the electric distance, the reactive power balance degree and the active power balance degree, and carrying out network partitioning by taking the distribution network planning as an application scene.
The traditional cluster-oriented planning network partition comprehensive index generally takes the structure and the function as principles, namely, the structure meets the conditions that the intra-cluster connection is tight, the inter-cluster connection is sparse, and the intra-cluster power exchange is facilitated; functionally, the group should have source complementarity to reduce power exchange between the groups and promote the intra-group consumption of DG; the network partitioning by using the indexes determined by the structural and functional principles is helpful to solve the power balance under the static condition, but with the superposition of high-permeability DG grid connection and multi-type loads, the running environment of the cluster has more uncertainty and volatility, so that higher requirements are put on the dynamic power balance capability of the cluster, namely the cluster needs to have stronger flexibility. At present, related researches are provided for providing corresponding flexibility supply and demand balance indexes from the aspect of climbing flexibility shortage when network partitioning is carried out, so that the capacity of power real-time balance in a cluster is improved to a certain extent, but the flexibility of the flexibility resources in the cluster in response to the net load demand in space is ignored from the aspect of network side flexibility, and the method has certain limitation.
In summary, how to solve the problems of mismatching of supply and demand and unbalanced power between feeder lines in each cluster by considering the combination between each flexible resource node and the payload node in the coordination system and the combination between different lines, and to set flexibility indexes from different angles, fully exert the adjustment capability of flexible resources to improve the autonomous characteristics of the clusters is a problem to be solved by those skilled in the art.
Disclosure of Invention
The invention can overcome the defects of the existing network partitioning method, and provides the partitioning method for the distributed power supply power distribution network with high permeability, which considers flexibility, so that the power distribution network can be reasonably partitioned, and the adjustment capability of flexible resources in each cluster can be exerted to the greatest extent, thereby being beneficial to solving the problem of insufficient operation flexibility in the later stage of power distribution network planning under the background of high-permeability DG grid connection and multi-type load occurrence.
In order to achieve the aim of the invention, the invention adopts the following technical scheme:
the invention relates to a power distribution network partitioning method with high permeability distributed power supply considering flexibility, which is characterized by comprising the following steps:
a: establishing a corresponding cluster flexibility supply model and a cluster flexibility demand model;
A1：let it be assumed that the group G is adjustable at the ith node i And energy storage element ES i The output, charge and discharge power at time t are respectivelyThen the capability of up-regulating all flexibility resources in the c-th cluster at the time t total flexibility is obtained by using the formula (1-1) and the formula (1-2) respectively>And Down-regulating ability->
In the formula (1-1) and the formula (1-2),respectively represent the adjustable units G at the ith node i Active maximum, minimum force and ramp up and down rate limits; />Representing a node set accessed with an adjustable unit G in a c-th cluster; τ represents the response time scale; />Respectively representing the maximum charge and discharge power of the stored energy at the ith node; />A node set which represents the energy storage element ES accessed in the c-th cluster; />In the positive direction, the energy storage element ES representing the ith node i Discharging, negative, the energy storage element ES representing the ith node i Charging; wherein,energy storage element ES for the ith node i Maximum capacity of installation;
a2: obtaining quantitative indexes of the net load in the c-th cluster at the t moment by using the formulas (1-3) and (1-4) respectivelyQuantization index of flexibility requirement +.>
In the formulas (1-3) and (1-4),and->Respectively representing photovoltaic and wind power original output and load active data of an ith node at the moment t; />The net load of the c-th cluster at the time t and the time t+1 respectively; />The number of load nodes in the c-th cluster;
b: constructing a net load adaptation rate index, a branch load margin index and a modularity value index, and endowing a certain weight coefficient to form a comprehensive index of a network partition;
b1: constructing a payload adaptation rate index of the c-th cluster in a t period by using the formulas (1-5)
In the formula (1-5); sigma represents the upper limit of the ratio of the total capacity of flexible resources to the payload demand within the cluster;
representing the comprehensive adjustment capability of the c-th cluster in the t period, which corresponds to the actual change direction of the flexibility requirement, and the comprehensive adjustment capability is represented by the following formula (1-6):
in the formula (1-6), N c Dividing the number of clusters;
during the planning period T, the payload adaptation rate index is calculated by using the formulas (1-7)Performing per unit value to obtain a comprehensive payload adaptation rate index L of the whole power distribution network AR ：
B2: constructing branch load margin index of the c-th cluster at t moment by using (1-8)
In the formula (1-8), I max For the branch to transmit the maximum value of the current, I ij,t The transmission current of a line between the ith node and the jth node at the t moment is used as the transmission current;the branch number in the c-th cluster;
during the planning period T, the branch load margin index is calculated by using the formulas (1-7)Performing per unit value to obtain a comprehensive branch load margin index H of the whole power distribution network BM ：
B3: constructing a modularity index ρ by using the formula (1-10):
in the formula (1-10), A i,j The weight of the edge between the ith node and the jth node is obtained by the formula (1-11); n is a node set in the power distribution network; k (k) i ＝∑ j∈N A i,j The sum of the weights of all edges connected with the ith node; k (k) j ＝∑ i∈N A i,j The sum of the weights of all edges connected with the jth node; m= (Σ) i∈N ∑ j∈N A i,j ) 2 represents the sum of the weights of all the edges connected to the nodes in the distribution network; sigma (i, j) is an optimization variable of the partitioning problem, if sigma (i, j) =1, the i node and the j node are located in the same area, otherwise, the i node and the j node are representedThe points are not in the same area;
A i,j ＝1-L ij /max(L) (1-11)
in the formula (1-11), d ij The ratio between the voltage change value of the jth node and the voltage change value of the ith node after the jth node injects the unit reactive power is represented; l (L) ij The spatial electric distance between the ith node and the jth node for considering the influence of all nodes is obtained by the formulas (1-12); max (L) represents the maximum value of the elements in the electrical distance matrix L;
b4: constructing a comprehensive index gamma of the network partition by using the formula (1-13):
max γ＝λ 1 ρ+λ 2 L AR +λ 3 H BM (1-13)
in the formula (1-13), lambda 1 、λ 2 、λ 3 Modular value index rho and net load adaptation rate index L respectively corresponding to network partitions AR And branch load margin index H BM Weights of (2);
c: taking the comprehensive dividing index gamma as an optimization target of an improved FN algorithm and calculating to obtain an optimal partition result;
c1: initializing network partition, regarding each node in the power distribution network to be partitioned as an independent cluster, and calculating initial value rho of modularity index 0 Let the real-time value of the modularity of the distribution network be ρ new And initialize ρ new ＝ρ 0 The method comprises the steps of carrying out a first treatment on the surface of the Let the real-time value of the comprehensive index of the power distribution network partition be gamma new And initialize gamma new ＝0；
C2: according to the cluster division condition in the current power distribution network, calculating a module degree increment matrix delta rho 'under various cluster combination conditions before each combination to obtain a corresponding module degree value matrix rho' =rho new +Δρ'；
And C3: before any ith independent cluster and jth independent cluster in the power distribution network are combined, judging the ith independent cluster and the jth independent clusterIf there is a branch formed by nodes between j independent clusters, combining the ith independent cluster with the jth independent cluster, and calculating the combined payload adaptation rate indexAnd branch load margin index>Otherwise, the payload adaptation rate index after combining the ith independent cluster and the jth independent cluster is +.>And branch load margin indexMarking as zero, thereby completing the combination judgment and calculation of all clusters and obtaining a comprehensive index value judgment matrix gamma' =lambda after all combination conditions are considered 1 ρ'+λ 2 L AR '+λ 3 H BM' The method comprises the steps of carrying out a first treatment on the surface of the Wherein L is AR’ 、H BM' Respectively representing a matrix formed by each payload adaptation rate index and a matrix formed by each branch load margin index under the condition that all nodes in the power distribution network are combined;
and C4: selecting two clusters corresponding to the maximum value in the comprehensive index value judgment matrix gamma' in the step C3 for merging, and calculating the net load adaptation rate index of the merged clustersBranch load margin index>And a modularity index value ρ' max Thereby obtaining the comprehensive index gamma 'of the combined clusters' max ；
C5: will ρ' max Assigning a value to ρ new Will be gamma' max Assignment to gamma new After that, the step C2 is returned to be sequentially executed untilγ new And the partition is not increased any more, so that the optimal partition result is obtained.
Compared with the prior art, the invention has the following advantages:
1. according to the invention, network partition comprehensive indexes considering flexibility and structural characteristics are established, on one hand, the matching degree of the flexibility resource adjustment margin in the cluster relative to flexibility requirements is improved, on the other hand, the space flexibility is provided for the response of the flexibility resource to the net load fluctuation output, the adjustment capability of the flexibility resource in the cluster is brought into full play to the greatest extent, the autonomous characteristic of the cluster is improved, and the problem of difficult digestion caused by insufficient operation flexibility in the later stage of power distribution network planning under the complex background of high-permeability DG grid connection and multi-type load occurrence is solved.
2. The space dimension of the flexibility of the invention provides the branch load margin index, provides a flexible transmission channel for the response of the flexible resource to the net load output, ensures the balance of the distribution of the flexible resource power in the cluster among all feeders, and further solves the problem of insufficient flexibility of the later operation of the cluster planning.
3. Aiming at the defects of the traditional FN algorithm, the invention provides an improved FN algorithm, a judging module for judging the connectivity condition of nodes among the clusters to be combined is added in the algorithm, the unnecessary combined calculation process is reduced, the searching efficiency during algorithm optimizing is improved, and the accurate and efficient partition process is ensured when the comprehensive partition index of the method is used as the optimizing target of the improved FN algorithm.
Drawings
Figure 1 is a flow chart of the steps of the method of the present invention for partitioning a distributed power supply distribution network with high permeability, taking into account flexibility.
Fig. 2 is a flow chart of the present invention for power distribution network partitioning using the improved FN algorithm.
Detailed Description
The invention is further illustrated by the following figures and examples, which are not intended to be limiting.
Examples: the partitioning method of the distributed power supply distribution network with high permeability considering flexibility is based on angles of a source side and a network side, sets related flexibility indexes, combines coupling characteristics in a cluster to realize reasonable partitioning of the distribution network, and specifically, as shown in fig. 1, comprises the following steps:
a: and respectively establishing a corresponding cluster flexibility supply model and a cluster flexibility demand model according to the adjustment characteristics of the existing flexibility resources in the network and the fluctuation of the original DGs and loads.
Cluster flexibility supply and demand models are as follows:
a1: based on existing studies of the supply characteristics of each flexible resource, the main flexible resource in the network before planning is considered as an adjustable conventional unit and an energy storage element. Let it be assumed that the group G is adjustable at the ith node i And energy storage element ES i The output, charge and discharge power at time t are respectivelyThen the capability of up-regulating all flexibility resources in the c-th cluster at the time t total flexibility is obtained by using the formula (1-1) and the formula (1-2) respectively>And Down-regulating ability->
In the formula (1-1) and the formula (1-2),respectively represent the adjustable units G at the ith node i Maximum and minimum active force, and up and downA hill climbing rate limit; />Representing a node set accessed with an adjustable unit G in a c-th cluster; τ represents the response time scale; />Respectively representing the maximum charge and discharge power of the stored energy at the ith node; />A node set which represents the energy storage element ES accessed in the c-th cluster; />In the positive direction, the energy storage element ES representing the ith node i Discharging, negative, the energy storage element ES representing the ith node i Charging; wherein,energy storage element ES for the ith node i Maximum capacity of installation;
a2: flexibility requirements within the cluster. Cluster flexibility requirements are derived from the volatility and randomness of the original DG and load, and the quantized index of the payload in the c-th cluster at time t is obtained by using equations (1-3) and (1-4), respectivelyQuantization index of flexibility requirement +.>
In the formulas (1-3) and (1-4),and->Respectively representing photovoltaic and wind power original output and load active data of an ith node at the moment t; />The net load of the c-th cluster at the time t and the time t+1 respectively; />The number of load nodes in the c-th cluster;
b: based on two aspects of source side flexibility and network side flexibility, a net load adaptation rate index describing flexibility supply in a cluster relative to flexibility demand matching degree, a branch load margin index describing spatial transmission characteristics of flexible resource output in the cluster, a modularity value index describing cluster structural characteristics are given, and a certain weight coefficient is given to form a comprehensive index of network partition.
B1: constructing a payload adaptation rate index of the c-th cluster in a t period by using the formulas (1-5)
In the formulas (1-5), sigma represents the upper limit of the total adjustment capacity of flexible resources in the cluster relative to the ratio of the net load demand;
comprehensive representing actual change direction of flexibility requirement of c-th cluster in t periodAnd the ability to be regulated, and is represented by the formula (1-6):
in the formula (1-6), N c Dividing the number of clusters;
during the planning period T, the payload adaptation rate index is calculated by using the formulas (1-7)Performing per unit value to obtain a comprehensive payload adaptation rate index L of the whole power distribution network AR ：
In the formula (1-7), L AR The method comprises the steps of (1) indicating the payload adaptation rate of the whole system; t is a planning period;representing the maximum value of the payload adaptation rate in all clusters throughout the period.
B2: constructing branch load margin index of the c-th cluster at t moment by using (1-8)
In the formula (1-8), I max For the branch to transmit the maximum value of the current, I ij,t The transmission current of a line between the ith node and the jth node at the t moment is used as the transmission current;the branch number in the c-th cluster;
during the planning period T, the branch load margin index is calculated by using the formulas (1-9)Performing per unit value to obtain a comprehensive branch load margin index H of the whole power distribution network BM ：
B3: constructing a modularity index ρ by using the formula (1-10):
in the formula (1-10), A i,j The weight of the edge between the ith node and the jth node is obtained by the formula (1-11); n is a node set in the power distribution network; k (k) i ＝∑ j∈N A i,j The sum of the weights of all edges connected with the ith node; k (k) j ＝∑ i∈N A i,j The sum of the weights of all edges connected with the jth node; m= (Σ) i∈N ∑ j∈N A i,j ) 2 represents the sum of the weights of all the edges connected to the nodes in the distribution network; sigma (i, j) is an optimization variable of the division problem, if sigma (i, j) =1, the i node and the j node are located in the same area, otherwise, the i node and the j node are not located in the same area;
A i,j ＝1-L ij /max(L) (1-11)
network edge weight A i,j Based on the spatial electrical distance representation, the spatial electrical distance is used for measuring the tightness degree of electrical connection between two nodes which are influenced by other nodes in an n-dimensional space, and is generally obtained by a reactive sensitivity relation, and the expression is as follows:
ΔV＝S QV *ΔQ (1-12)
in the formula (1-12): s is S QV In order to consider the reactive sensitivity matrix of the line active effect, Δv, Δq represent the voltage and reactive variation values, respectively.
In the formula (1-11), d ij The ratio between the voltage change value of the jth node and the voltage change value of the ith node after the jth node is injected with the unit reactive power is represented, and the ratio is obtained by the formula (1-14); l (L) ij The spatial electric distance between the ith node and the jth node for considering the influence of all nodes is obtained by the formulas (1-13); max (L) represents the maximum value of the elements in the electrical distance matrix L;
b4: constructing a comprehensive index gamma of the network partition by using the formula (1-15):
max γ＝λ 1 ρ+λ 2 L AR +λ 3 H BM (1-15)
in the formula (1-15), lambda 1 、λ 2 、λ 3 Modular value index rho and net load adaptation rate index L respectively corresponding to network partitions AR And branch load margin index H BM Weights of (2);
c: providing an improved FN algorithm, adding a judging module for the connectivity condition of nodes among the clusters to be combined, taking the comprehensive dividing index as an optimizing target of the algorithm, calculating to obtain an optimal partitioning result, specifically, as shown in fig. 2, and performing the following steps:
c1: initializing network partition, regarding each node in the power distribution network to be partitioned as an independent cluster, and calculating initial value rho of modularity index 0 Let the real-time value of the modularity of the distribution network be ρ new And initialize ρ new ＝ρ 0 The method comprises the steps of carrying out a first treatment on the surface of the Let the real-time value of the comprehensive index of the power distribution network partition be gamma new And initialize gamma new ＝0；
C2: according to the cluster division condition in the current power distribution network, calculating a module degree increment matrix delta rho' under various cluster combination conditions before each combination to obtain a corresponding moduleDegree matrix ρ' =ρ new +Δρ'；
And C3: before any i independent cluster and j independent cluster in the power distribution network are combined, judging whether a node between the i independent cluster and the j independent cluster forms a branch, if so, combining the i independent cluster and the j independent cluster, and calculating a combined payload adaptation rate indexAnd branch load margin index>Otherwise, the payload adaptation rate index after combining the ith independent cluster and the jth independent cluster is +.>And branch load margin indexMarking as zero, thereby completing the combination judgment and calculation of all clusters and obtaining a comprehensive index value judgment matrix gamma' =lambda after all combination conditions are considered 1 ρ'+λ 2 L AR '+λ 3 H BM' The method comprises the steps of carrying out a first treatment on the surface of the Wherein L is AR’ 、H BM' Respectively representing a matrix formed by each net load adaptation rate index value and each branch load margin index value under the condition that all nodes in the power distribution network are combined;
and C4: selecting two clusters corresponding to the maximum value in the comprehensive index value judgment matrix gamma' in the step C3 for merging, and calculating the net load adaptation rate index of the merged clustersBranch load margin index>And a modularity index value ρ' max Thereby obtaining the comprehensive index gamma 'of the combined clusters' max ；
C5: will ρ' max Assigning a value to ρ new Will be gamma' max Assignment to gamma new After that, the step C2 is returned to be sequentially executed until gamma new And the partition is not increased any more, so that the optimal partition result is obtained.
Claims (1)
1. A power distribution network partitioning method with high-permeability distributed power supply considering flexibility is characterized by comprising the following steps:
a: establishing a corresponding cluster flexibility supply model and a cluster flexibility demand model;
a1: let it be assumed that the group G is adjustable at the ith node i And energy storage element ES i The output, charge and discharge power at time t are respectivelyThen the capability of up-regulating all flexibility resources in the c-th cluster at the time t total flexibility is obtained by using the formula (1-1) and the formula (1-2) respectively>And Down-regulating ability->
In the formula (1-1) and the formula (1-2),respectively represent the adjustable units G at the ith node i Maximum and minimum active force and the following components,A downhill climb rate limit; />Representing a node set accessed with an adjustable unit G in a c-th cluster; τ represents the response time scale; />Respectively representing the maximum charge and discharge power of the stored energy at the ith node; />A node set which represents the energy storage element ES accessed in the c-th cluster; />In the positive direction, the energy storage element ES representing the ith node i Discharging, negative, the energy storage element ES representing the ith node i Charging; wherein, energy storage element ES for the ith node i Maximum capacity of installation;
a2: obtaining quantitative indexes of the net load in the c-th cluster at the t moment by using the formulas (1-3) and (1-4) respectivelyQuantization index of flexibility requirement +.>
In the formulas (1-3) and (1-4),and->Respectively representing photovoltaic and wind power original output and load active data of an ith node at the moment t; />The net load of the c-th cluster at the time t and the time t+1 respectively;the number of load nodes in the c-th cluster;
b: constructing a net load adaptation rate index, a branch load margin index and a modularity value index, and endowing a certain weight coefficient to form a comprehensive index of a network partition;
b1: constructing a payload adaptation rate index of the c-th cluster in a t period by using the formulas (1-5)
In the formula (1-5); sigma represents the upper limit of the ratio of the total capacity of flexible resources to the payload demand within the cluster;
representing that the c-th cluster corresponds to the prodigy in period tThe comprehensive regulation capability of the actual change direction of the activity requirement is represented by the following formula (1-6):
in the formula (1-6), N c Dividing the number of clusters;
during the planning period T, the payload adaptation rate index is calculated by using the formulas (1-7)Performing per unit value to obtain a comprehensive payload adaptation rate index L of the whole power distribution network AR ：
B2: constructing branch load margin index of the c-th cluster at t moment by using (1-8)
In the formula (1-8), I max For the branch to transmit the maximum value of the current, I ij,t The transmission current of a line between the ith node and the jth node at the t moment is used as the transmission current;the branch number in the c-th cluster;
during the planning period T, the branch load margin index is calculated by using the formulas (1-7)Performing per unit value to obtain a comprehensive branch load margin index H of the whole power distribution network BM ：
B3: constructing a modularity index ρ by using the formula (1-10):
in the formula (1-10), A i,j The weight of the edge between the ith node and the jth node is obtained by the formula (1-11); n is a node set in the power distribution network; k (k) i ＝∑ j∈N A i，j The sum of the weights of all edges connected with the ith node; k (k) j ＝∑ i∈N A i，j The sum of the weights of all edges connected with the jth node; m= (Σ) i∈N ∑ j∈N A i，j ) 2 represents the sum of the weights of all the edges connected to the nodes in the distribution network; sigma (i, j) is an optimization variable of the division problem, if sigma (i, j) =1, the i node and the j node are located in the same area, otherwise, the i node and the j node are not located in the same area;
A i，j ＝l-L ij /max(L) (1-11)
in the formula (1-11), d ij The ratio between the voltage change value of the jth node and the voltage change value of the ith node after the jth node injects the unit reactive power is represented; l (L) ij The spatial electric distance between the ith node and the jth node for considering the influence of all nodes is obtained by the formulas (1-12); max (L) represents the maximum value of the elements in the electrical distance matrix L;
b4: constructing a comprehensive index gamma of the network partition by using the formula (1-13):
max γ＝λ 1 ρ+λ 2 L AR +λ 3 H BM (1-13)
in the formula (1-13), lambda 1 、λ 2 、λ 3 Modular value index rho and net load adaptation rate index L respectively corresponding to network partitions AR And branch load margin index H BM Weights of (2);
c: taking the comprehensive dividing index gamma as an optimization target of an improved FN algorithm and calculating to obtain an optimal partition result;
c1: initializing network partition, regarding each node in the power distribution network to be partitioned as an independent cluster, and calculating initial value rho of modularity index 0 Let the real-time value of the modularity of the distribution network be ρ new And initialize ρ new ＝ρ 0 The method comprises the steps of carrying out a first treatment on the surface of the Let the real-time value of the comprehensive index of the power distribution network partition be gamma new And initialize gamma new ＝0；
C2: according to the cluster division condition in the current power distribution network, calculating a module degree increment matrix delta rho 'under various cluster combination conditions before each combination to obtain a corresponding module degree value matrix rho' =rho new +△ρ′；
And C3: before any i independent cluster and j independent cluster in the power distribution network are combined, judging whether a node between the i independent cluster and the j independent cluster forms a branch, if so, combining the i independent cluster and the j independent cluster, and calculating a combined payload adaptation rate indexAnd branch load margin index>Otherwise, the payload adaptation rate index after combining the ith independent cluster and the jth independent cluster is +.>And branch load margin index>Marking as zero, thereby completing the combination judgment and calculation of all clusters and obtaining a comprehensive index value judgment matrix gamma' =lambda after all combination conditions are considered 1 ρ′+λ 2 L AR′ +λ 3 H BM′ The method comprises the steps of carrying out a first treatment on the surface of the Wherein L is AR′ 、H BM′ Respectively representing a matrix formed by each payload adaptation rate index and a matrix formed by each branch load margin index under the condition that all nodes in the power distribution network are combined;
and C4: selecting two clusters corresponding to the maximum value in the comprehensive index value judgment matrix gamma' in the step C3 for merging, and calculating the net load adaptation rate index of the merged clustersBranch load margin index>And a modularity index value ρ' max Thereby obtaining the comprehensive index gamma 'of the combined clusters' max ；
C5: will ρ' max Assigning a value to ρ new Will be gamma' max Assignment to gamma new After that, the step C2 is returned to be sequentially executed until gamma new And the partition is not increased any more, so that the optimal partition result is obtained.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210193972.1A CN114552669B (en) | 2022-03-01 | 2022-03-01 | Flexibility-considered partitioning method for distributed power supply distribution network containing high permeability |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210193972.1A CN114552669B (en) | 2022-03-01 | 2022-03-01 | Flexibility-considered partitioning method for distributed power supply distribution network containing high permeability |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114552669A CN114552669A (en) | 2022-05-27 |
| CN114552669B true CN114552669B (en) | 2024-03-12 |
Family
ID=81662338
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210193972.1A Active CN114552669B (en) | 2022-03-01 | 2022-03-01 | Flexibility-considered partitioning method for distributed power supply distribution network containing high permeability |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114552669B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117439090B (en) * | 2023-12-19 | 2024-04-02 | 浙江大学 | Flexible resource allocation or scheduling method taking flexible adjustment coefficient as index |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108448620A (en) * | 2018-04-04 | 2018-08-24 | 合肥工业大学 | High permeability distributed generation resource assemblage classification method based on integrated performance index |
| CN110429649A (en) * | 2019-08-13 | 2019-11-08 | 合肥工业大学 | Consider the high permeability renewable energy assemblage classification method of flexibility |
| CN110518575A (en) * | 2019-08-02 | 2019-11-29 | 南京理工大学 | Multiple Time Scales active distribution network voltage optimization control method based on region division |
| CN113364058A (en) * | 2020-03-05 | 2021-09-07 | 中国电力科学研究院有限公司 | Reactive power control method and system for power distribution network |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103577892A (en) * | 2013-10-30 | 2014-02-12 | 河海大学 | Progressive intelligent power distribution system scheduling method |
-
2022
- 2022-03-01 CN CN202210193972.1A patent/CN114552669B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108448620A (en) * | 2018-04-04 | 2018-08-24 | 合肥工业大学 | High permeability distributed generation resource assemblage classification method based on integrated performance index |
| CN110518575A (en) * | 2019-08-02 | 2019-11-29 | 南京理工大学 | Multiple Time Scales active distribution network voltage optimization control method based on region division |
| CN110429649A (en) * | 2019-08-13 | 2019-11-08 | 合肥工业大学 | Consider the high permeability renewable energy assemblage classification method of flexibility |
| CN113364058A (en) * | 2020-03-05 | 2021-09-07 | 中国电力科学研究院有限公司 | Reactive power control method and system for power distribution network |
Non-Patent Citations (2)
| Title |
|---|
| 含高渗透分布式电源配电网灵活性提升优化调度方法;王洪坤;王守相;潘志新;王建明;;电力系统自动化;20180702(第15期);全文 * |
| 采用综合性能指标的高渗透率分布式电源集群划分方法;丁明;刘先放;毕锐;胡迪;叶彬;张晶晶;;电力系统自动化;20180606(第15期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114552669A (en) | 2022-05-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106410861B (en) | A kind of micro-capacitance sensor optimization operation real-time control method based on schedulable ability | |
| CN110429649B (en) | High-permeability renewable energy cluster division method considering flexibility | |
| CN108092324B (en) | AGC control system and control method for wind power participating in peak shaving frequency modulation | |
| CN109657993B (en) | Non-cooperative game-based automatic demand response method for energy local area network energy storage system | |
| CN107994595A (en) | A kind of system of peak load shifting control method and system and the application control method | |
| CN107453381B (en) | Electric car cluster power regulating method and system based on two stages cross-over control | |
| CN113013929A (en) | Load curve adjustment-oriented active power distribution network simulation optimization operation method | |
| CN110676849B (en) | Method for constructing islanding micro-grid group energy scheduling model | |
| CN105896596B (en) | A kind of the wind power layering smoothing system and its method of consideration Demand Side Response | |
| CN112994097A (en) | High-proportion distributed photovoltaic cooperative control method based on intelligent distribution transformer terminal system | |
| CN114552669B (en) | Flexibility-considered partitioning method for distributed power supply distribution network containing high permeability | |
| CN114336785B (en) | Distributed power supply group control and group dispatching control method and device based on grid clustering | |
| CN109711706A (en) | Consider the active distribution network substation planning method of distributed generation resource and demand response | |
| CN111682536B (en) | Random-robust optimization operation method for virtual power plant participating in dual market before day | |
| CN110852495A (en) | Site selection method for distributed energy storage power station | |
| CN110323779B (en) | Method and system for dynamically aggregating power of distributed power generation and energy storage device | |
| CN110808600A (en) | Method and system for calculating access of battery energy storage station in transformer substation | |
| CN107591794A (en) | Active distribution network source storage capacity configuration optimizing method based on load classification | |
| CN113837449B (en) | Centralized optimization scheduling method for power grid system participated by virtual power plant | |
| CN113673141B (en) | Energy router modeling and optimization control method based on data driving | |
| CN110970939B (en) | Distributed energy cluster optimization method and system | |
| CN115207966A (en) | Active power distribution network distributed energy multi-scale consumption method | |
| CN115021327A (en) | Collaborative day-ahead optimization scheduling method of interconnected micro-grid system based on virtual energy storage | |
| CN113393085B (en) | Cluster dividing method considering flexibility supply and demand balance and response speed | |
| CN112103987B (en) | Reactive voltage two-stage partition and control method for wind power access power grid |
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 |