CN116316887A - Distribution line distributed photovoltaic bearing capacity assessment method based on N-1 safety criterion - Google Patents
Distribution line distributed photovoltaic bearing capacity assessment method based on N-1 safety criterion Download PDFInfo
- Publication number
- CN116316887A CN116316887A CN202310285711.7A CN202310285711A CN116316887A CN 116316887 A CN116316887 A CN 116316887A CN 202310285711 A CN202310285711 A CN 202310285711A CN 116316887 A CN116316887 A CN 116316887A
- Authority
- CN
- China
- Prior art keywords
- distribution line
- scene
- node
- moment
- distributed photovoltaic
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000011156 evaluation Methods 0.000 claims abstract description 8
- 238000013210 evaluation model Methods 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims abstract description 5
- 230000009467 reduction Effects 0.000 claims abstract description 4
- 230000005540 biological transmission Effects 0.000 claims description 23
- 238000004146 energy storage Methods 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 6
- 238000009434 installation Methods 0.000 claims description 6
- 241000243621 Vandenboschia maxima Species 0.000 claims 1
- 238000003064 k means clustering Methods 0.000 abstract description 6
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 238000001228 spectrum Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
Images
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
-
- 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
- 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/48—Controlling the sharing of the in-phase component
-
- 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/50—Controlling the sharing of the out-of-phase component
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention provides a distribution line distributed photovoltaic bearing capacity assessment method based on an N-1 safety criterion, which relates to the technical field of distribution line bearing capacity assessment and specifically comprises the following steps: collecting data of the distributed photovoltaic historical output data, and obtaining a typical output scene through a K-means clustering algorithm; and combining uncertainty characteristics of source load, and taking the acceptance capacity of the distributed photovoltaic power supply in a fault scene into consideration to establish a distribution line bearing capacity evaluation model based on an N-1 safety criterion. Based on a decoupling strategy, performing decoupling dimension reduction calculation on a feasible domain considering an N-1 fault scene; the feasible region considering the N-1 security constraint can be characterized as the intersection of the feasible regions in all fault scenes, and the evaluation result of the distributed photovoltaic power supply carrying capacity of the distribution line in the scene is obtained. And checking whether all the typical output scenes obtained through the K-means algorithm are traversed. The method provided by the invention is beneficial to providing reference for reasonable configuration of distribution lines and distributed photovoltaic power sources.
Description
Technical Field
The patent belongs to the technical field of distribution line bearing capacity assessment and relates to a distribution line distributed photovoltaic bearing capacity assessment method based on N-1 safety criteria.
Background
In order to reduce carbon emission, the generation proportion of new energy is greatly increased in recent years, and particularly, distributed generation is applied in a large scale. The capacity and permeability of a distributed power source represented by photovoltaic and wind power on a power source side in a power system are continuously increased, and the fluctuation and uncertainty of the output of the distributed power source provide challenges for the running stability and reliability of a power grid to the safe running of a distribution line. Therefore, the assessment of the bearing capacity of the distribution line has important significance for the safe operation of the distribution line.
The existing distributed photovoltaic access capacity is evaluated under the condition of safe operation of a distribution line, and in actual conditions, when the distribution line faces to faults or overhauls and the like, in order to ensure that a terminal user does not have the power failure and the like, the reliability of the distribution line is ensured, and the load rate of the distribution line is required to meet N-1 safety standards. Therefore, there is a need to propose a distribution line distributed photovoltaic load capacity assessment method based on N-1 safety guidelines.
Disclosure of Invention
In order to solve the defects and the shortcomings in the prior art, the patent adopts a distributed photovoltaic bearing capacity evaluation method of the distribution line based on the N-1 safety criterion, and the bearing capacity of the distribution line is evaluated under the condition of considering the receiving capacity of the distributed photovoltaic power supply under the fault scene.
Specifically, the power distribution line distributed photovoltaic bearing capacity evaluation method based on the N-1 safety rule is characterized by comprising the following steps of:
step 1: combining uncertainty characteristics of source load, and considering acceptance capacity of a distributed photovoltaic power supply in a fault scene to establish a distribution line bearing capacity assessment method based on an N-1 safety criterion;
step 2: and carrying out decoupling dimension reduction calculation on the feasible domain in the N-1 safety criterion fault scene based on a decoupling strategy.
In the step 1, based on the N-1 safety criterion, the installation capacity of the distributed photovoltaic power supply is taken as an objective function, and the line capacity, the node voltage, the tide and the short-circuit current are taken as constraint conditions:
the objective function is that the distributed photovoltaic power supply of the line has the largest installation capacity, and the expression is as follows:
Moreover, before the occurrence of the fault scene, the constraint condition of the distributed photovoltaic power supply bearing capacity evaluation model is as follows:
1) Line capacity constraints
Wherein:the active power flowing through the ij branch t of the distribution line at the moment of the scene s is obtained; />The reactive power flowing through the ij branch t of the distribution line at the moment of the scene s is obtained; />Line capacity for ij branch of distribution line;
2) Tidal current constraint
V s,j,t =V s,i,t -(r ij P s,ij,t +x ij Q s,ij,t )
Wherein: p (P) s,ij,t The active power flowing through the branch ij at the moment t of the distribution line under the scene s is obtained; q (Q) s,ij,t The reactive power flowing through the branch ij at the moment t of the distribution line in the scene s is obtained; p (P) s,jk,t Active power flowing out of the node j and flowing to the node k at the moment of the distribution line t in the scene s; q (Q) s,jk,t The reactive power flows out of the node j and flows to the node k at the moment of the distribution line t in the scene s;the active power is the load active power of a node j at the moment t of the distribution line under the scene s; />The reactive power is the load reactive power of a node j at the moment t of the distribution line under the scene s; />DG power of a node j at t moment of a distribution line in a scene s; />The energy storage and charging power of a node j at the moment t of the distribution line in the scene s; />The energy storage and discharge power of a node j at the moment t of the distribution line in the scene s is obtained; r is (r) ij The resistance of branch ij; x is x ij Is the reactance of branch ij.
Moreover, after the fault scene occurs, the constraint conditions of the distributed photovoltaic power supply bearing capacity evaluation model are as follows:
1) Line capacity constraints
For the transmission line, after a fault occurs or an overhaul scene is met, the condition that the N-1 check is met means that: when one power transmission line fails or overhauls and stops running, the load of the power transmission line can be transferred to other power transmission lines which normally run, and the short-time allowable overload capacity of the normal running line is not more than the rated capacity of a single line, and the power transmission line is based on the N-1 safety constraint as shown in the formula:
|P(t)≤(N-1)S|
wherein: n is the number of transmission lines; s is the rated capacity of a single power transmission line; p (t) is the net load of the transmission line at the moment t;
2) Tidal current constraint
V f,j,t =V f,i,t -(r ij P f,ij,t +x ij Q f,ij,t )
Wherein: p (P) f,ij,t The active power flowing through the branch ij at the moment t of the distribution line in the fault scene f; q (Q) f,ij,t Reactive power flowing through a branch ij at the moment t of the distribution line in the fault scene f; p (P) f,jk,t Active power flowing out of the node j and flowing to the node k at the moment of the distribution line t in the fault scene f; q (Q) f,jk,t The reactive power flows out of the node j and flows to the node k at the moment of the distribution line t in the fault scene f;the active power is the load active power of a node j at the moment t of the distribution line under the fault scene f; />The reactive power is the load reactive power of a node j at the moment t of the distribution line in the fault scene f; />The DG power of a node j at the moment t of the distribution line in the fault scene f; />The energy storage and charging power of a node j at the moment t of the distribution line in the fault scene f; />The energy storage and discharge power of a node j at the moment t of the distribution line in the fault scene f; r is (r) ij The resistance of branch ij; x is x ij Is the reactance of branch ij.
And in the step 2, the feasible domains meeting the constraint in the ith line fault scene and the constraint in the non-fault scene are obtained in parallel, the feasible domain considering the N-1 security constraint can be characterized as the intersection of the feasible domains in all fault scenes, and the specific flow is as follows:
(1) According to the maximum value and the minimum value in each dimension, constructing an initialized distributed photovoltaic access capacity feasible domain;
the minimum solution problem for the ith dimension:
S.t.
maximum solving problem for the i-th dimension:
S.t.
the obtained information under the non-fault scene and the fault sceneIs used as vertex V to construct a feasible domain omega c ;
(2) New vertex search: by shifting out the feasible region omega c Searching for new vertices; let the inequality expression of the kth boundary plane beObtaining a new vertex by solving the following linear programming problem:
wherein: v (V) p Is the feasible domain omega c Any point in the above;
the optimal solution of the linear programming is the top point of the constraint. Therefore, the optimal solution is the new vertex, and after all sides are translated, all the vertices are marked as V new ;
(3) Algorithm termination condition checking: calculating the difference between the feasible region volumes before and after the new vertex in the consideration (2) to be DeltaV; the smaller the delta V change, the smaller the feasible region change before and after adding the new vertex, thus, when delta V is smaller than the given threshold, the algorithm is terminated and the feasible region omega under the c-th line fault is recorded c To construct a feasible region Ω from the new vertices otherwise c Then return to (2);
when obtaining the feasible region Ω in all fault conditions c The feasible region Ω taking into account the N-1 security constraints is characterized as the feasible region Ω under all fault conditions c And (3) obtaining the carrying capacity evaluation result of the distributed photovoltaic power distribution line.
And step 2, the confidence degrees are ranked from large to small, the confidence degrees are opposite to the obtained bearing capacity ranking order, and the corresponding bearing capacity of the distributed photovoltaic power distribution line is determined after the confidence degrees are selected.
The application has the advantages and positive effects that:
the method for evaluating the distributed photovoltaic bearing capacity of the distribution line based on the N-1 safety criterion is based on the N-1 safety criterion, aims at the maximum bearing capacity of the distribution line, comprehensively considers the admittance capacity of the distributed photovoltaic power supply under a fault scene, and is favorable for providing references for reasonable configuration of the distribution line and the distributed photovoltaic power supply.
Drawings
FIG. 1 is a graph of historical output of a distributed photovoltaic power source.
Fig. 2 is a graph of typical output curves of clustered distributed photovoltaic power supplies.
Fig. 3 is a diagram of the structure of the line a.
Fig. 4 is a load time sequence characteristic curve.
Fig. 5 is a flow chart of the method.
Detailed Description
The present application is further described in detail below with reference to the attached drawings by way of specific examples, which are intended to be illustrative, not limiting, and not to limit the scope of the present application.
The distribution line distributed photovoltaic bearing capacity assessment method based on the N-1 safety rule comprises the following steps of:
step 1: and combining uncertainty characteristics of source load, and taking the acceptance capacity of the distributed photovoltaic power supply in a fault scene into consideration to establish a distribution line bearing capacity assessment method based on an N-1 safety criterion.
Step 2: and carrying out decoupling dimension reduction calculation on the feasible domain considering the N-1 fault scene based on a decoupling strategy.
In the step 1, a K-means clustering algorithm is adopted to obtain a distributed photovoltaic typical output scene, and the specific flow is as follows:
(1) Normalizing the collected historical distributed photovoltaic power output data;
(2) Randomly selecting K points from the normalized data to serve as mass centers;
(3) Respectively calculating the distance from the remaining data points to each centroid;
(4) Classifying the non-centroid points into clusters represented by centroids with shortest distances to the non-centroid points to form K clusters;
(5) According to the classified clusters, re-calculating mass centers in each cluster, and determining actual mass centers in the group;
(6) Repeating the operation of the step (3) until the iteration times or errors are smaller than the specified value, and obtaining the typical output scene of the distributed photovoltaic.
The specific calculation process is as follows:
note that sample dataset d= { X 1 ,X 2 ,...,X n N represents the number of samples, sample X of the ith i ={x i1 ,x i2 ,...,x ip P represents the attribute number or index category of the sample, and K represents the number of clusters. In K-means clustering, any two samples X i And X is f The euclidean distance between them is in the form:
wherein: d (x) i ,x f ) Is the distance between element i and the initial centroid f.
At the beginning of clustering, K samples are randomly selected as initial clustering centers, and then non-clustering center samples are classified into a category C represented by a clustering center h with the shortest distance to the non-clustering center samples according to the above formula h Is a kind of medium.
The cluster center is then updated by:
wherein: phi (X) i ) Represents the h th class C h Sample contained in (C) h I represents the h th category C h The number of samples contained in the sample, m h Represents the hClass C h Mean of the samples.
The goal of K-means clustering is a minimized square error criterion function, of the form:
the clustering center m is continuously updated through the formula h Samples phi (X) i ) Until the square error criterion function value no longer changes or the set number of iterations is reached.
After iteration is finished, the clustering center m at the moment is read h I.e. the typical sample we need, phi (X i ) I.e., the class to which each representative sample corresponds.
Based on the N-1 safety criterion, the step 1 is to use the maximum installation capacity of the distributed photovoltaic power supply of the distribution line as an objective function, and use the line capacity, the node voltage, the tide, the short-circuit current and the harmonic voltage distortion rate as constraint condition processes, and provides a distributed photovoltaic distribution line bearing capacity assessment method based on the N-1 criterion, which comprises the following steps:
the objective function is to maximize the distribution line distributed photovoltaic power installation capacity, expressed as follows:
Before a fault scene occurs, the constraint conditions of the distributed photovoltaic power supply bearing capacity assessment model are as follows:
1) Node voltage constraint
V i,min ≤V s,i,t ≤V i,max
Wherein: v (V) i,min A lower voltage limit for safe operation of the distribution line node i;V i,max the upper voltage limit for safe operation of the distribution line node i; v (V) s,i,t The voltage of the distribution line node i at the moment t under the scene s.
2) Line capacity constraints
Wherein:the active power flowing through the ij branch t of the 10kV medium-voltage distribution line at the moment of the scene s; />The reactive power flowing through the ij branch t of the 10kV medium-voltage distribution line at the moment of the scene s; />The line capacity of the ij branch of the 10kV medium voltage distribution line.
3) Tidal current constraint
V s,j,t =V s,i,t -(r ij P s,ij,t +x ij Q s,ij,t )
Wherein: p (P) s,ij,t The active power flowing through the branch ij at the moment t of the distribution line under the scene s is obtained; q (Q) s,ij,t The reactive power flowing through the branch ij at the moment t of the distribution line in the scene s is obtained; p (P) s,jk,t Active power flowing out of the node j and flowing to the node k at the moment of the distribution line t in the scene s; q (Q) s,jk,t The reactive power flows out of the node j and flows to the node k at the moment of the distribution line t in the scene s;the active power is the load active power of a node j at the moment t of the distribution line under the scene s; />The reactive power is the load reactive power of a node j at the moment t of the distribution line under the scene s; />DG power of a node j at t moment of a distribution line in a scene s; />The energy storage and charging power of a node j at the moment t of the distribution line in the scene s; />The energy storage and discharge power of a node j at the moment t of the distribution line in the scene s is obtained; r is (r) ij The resistance of branch ij; x is x ij Is the reactance of branch ij.
4) Harmonic voltage distortion rate constraint
The limitation of harmonic voltage distortion rate specified by national standards is expressed as follows:
THD V ≤THD V,max
wherein: THD (total heat transfer) V Is harmonic voltage distortion rate; THD (total heat transfer) V,max Maximum harmonic voltage distortion rate allowed for ensuring safe operation of the network; u (U) h-spectrum Is the harmonic voltage amplitude of the h-th typical spectrum; u (U) h The harmonic voltage amplitude of the h harmonic source; u (U) 1 Is the grid-connected point fundamental voltage.
5) Short circuit current restraint
Wherein:short-circuit current contributing to the original power supply; />Short-circuit current contributing to DG; i sc,max The maximum short-circuit current allowed to ensure safe operation of the network.
After the fault scene occurs, the constraint conditions of the distributed photovoltaic power supply bearing capacity evaluation model are as follows:
1) Node voltage constraint
V i,min ≤V s,i,t ≤V i,max
Wherein: v (V) i,min A lower voltage limit for safe operation of the distribution line node i; v (V) i,max The upper voltage limit for safe operation of the distribution line node i; v (V) s,i,t The voltage of the distribution line node i at the moment t under the scene s.
2) Line capacity constraints
For the transmission line, after a fault occurs or an overhaul scene is met, the condition that the N-1 check is met means that: when one transmission line fails or overhauls and stops running, the load of the transmission line can be transferred to other transmission lines which normally run, and the overload capacity of the normal running line is allowed to be smaller than the rated capacity of a single line. The N-1 based security constraint is shown in the formula.
|P 1 (t)≤(N-1)S|
Wherein: n is the number of transmission lines; s is the rated capacity of a single power transmission line; p (P) l And (t) is the net load of the transmission line at the moment t.
3) Tidal current constraint
V f,j,t =V f,i,t -(r ij P f,ij,t +x ij Q f,ij,t )
Wherein: p (P) f,ij,t The active power flowing through the branch ij at the moment t of the distribution line in the fault scene f; q (Q) f,ij,t Reactive power flowing through a branch ij at the moment t of the distribution line in the fault scene f; p (P) f,jk,t Active power flowing out of the node j and flowing to the node k at the moment of the distribution line t in the fault scene f; q (Q) f,jk,t The reactive power flows out of the node j and flows to the node k at the moment of the distribution line t in the fault scene f;the active power is the load active power of a node j at the moment t of the distribution line under the fault scene f; />The reactive power is the load reactive power of a node j at the moment t of the distribution line in the fault scene f; />The DG power of a node j at the moment t of the distribution line in the fault scene f; />The energy storage and charging power of a node j at the moment t of the distribution line in the fault scene f; />The energy storage and discharge power of a node j at the moment t of the distribution line in the fault scene f; r is (r) ij The resistance of branch ij; x is x ij Is the reactance of branch ij.
4) Harmonic voltage distortion rate constraint
The limitation of harmonic voltage distortion rate specified by national standards is expressed as follows:
THD V ≤THD V,max
wherein: THD (total heat transfer) V Is harmonic voltage distortion rate; THD (total heat transfer) V,max Maximum harmonic voltage distortion rate allowed for ensuring safe operation of the network; u (U) f,h-spectrum Is the harmonic voltage amplitude of the h-th typical frequency spectrum after the fault; u (U) f,h The harmonic voltage amplitude of the h harmonic source after the fault; u (U) 1 Is the grid-connected point fundamental voltage.
5) Short circuit current restraint
Wherein:short-circuit current contributing to the original power supply after failure; />Short-circuit current contributing to DG after failure; i sc,max The maximum short-circuit current allowed to ensure safe operation of the network.
Further, in step 2, the algorithm is initialized
(1) According to the maximum value and the minimum value in each dimension, constructing an initialized distributed photovoltaic access capacity feasible domain;
the minimum solution problem for the ith dimension:
S.t.
maximum solving problem for the i-th dimension:
S.t.
the obtained information under the non-fault scene and the fault sceneIs used as vertex V to construct a feasible domain omega c 。
(2) New vertex search: by shifting out the feasible region omega c Searching for new vertices; let the inequality expression of the kth boundary plane beObtaining a new vertex by solving the following linear programming problem:
wherein: v (V) p Is the feasible domain omega c Any point in the above;
the optimal solution of the linear programming is the top point of the constraint. Therefore, the optimal solution is the new vertex, and after all sides are translated, all the vertices are marked as V new ;
(3) Algorithm termination condition checking: calculating the difference between the feasible region volumes before and after the new vertex in the consideration (2) to be DeltaV; the smaller the delta V change, the smaller the feasible region change before and after adding the new vertex, thus, when delta V is smaller than the given threshold, the algorithm is terminated and the feasible region omega under the c-th line fault is recorded c To construct a feasible region Ω from the new vertices otherwise c Then return to (2);
when obtaining the feasible region Ω in all fault conditions c The feasible region Ω taking into account the N-1 security constraints is characterized as the feasible region Ω under all fault conditions c Is a complex of the two. By solving forAnd solving the optimization model to obtain a distributed photovoltaic power distribution line bearing capacity evaluation result based on an N-1 criterion.
Further, when the bearing capacity is smaller, the condition of line capacity out-of-limit cannot occur in different output scenes, so that the more the number of scenes for bearing capacity adaptation is, the higher the confidence is. And according to the confidence coefficient from large to small, namely the obtained bearing capacity is from small to large, the smaller the bearing capacity is, namely the larger the confidence coefficient is, and after the confidence coefficient is selected, the corresponding distributed photovoltaic power bearing capacity of the distribution line can be determined.
In the embodiment, six distribution lines in a certain area are selected, the numbers of the six distribution lines are a, b, c, d, e, f respectively, when the bearing capacity of the distribution lines is evaluated, K-means clustering is adopted to cut down a large number of distribution photovoltaic power output scenes obtained by the collection so as to obtain typical scenes and probability thereof, and the bearing capacity evaluation of the distribution lines of the distributed photovoltaic power is continuously carried out based on the obtained typical scenes.
(1) Historical scene of distributed photovoltaic power output
Drawing a large number of distributed photovoltaic power output curves according to the funding data, wherein the result is shown in the following graph:
the k-means clustering is adopted to cut down a large number of scenes, in the example, all scenes are divided into 10 types, the clustering center of each type is used as a cut-down scene, and the obtained cut-down typical output curve is shown in the following graph:
the probability of the scene data within each class and the various typical curves derived therefrom are shown in the following table:
TABLE 1 typical scenario and probability thereof
Category(s) | Number of scenes in class | Probability of occurrence | Output coefficient at |
1 | 40 | 0.11 | 0.67 |
2 | 47 | 0.13 | 0.42 |
3 | 19 | 0.05 | 0.21 |
4 | 20 | 0.05 | 0.82 |
5 | 46 | 0.13 | 0.68 |
6 | 67 | 0.18 | 0.56 |
7 | 44 | 0.12 | 0.07 |
8 | 12 | 0.04 | 0.55 |
9 | 43 | 0.12 | 0.43 |
10 | 27 | 0.07 | 0.39 |
(2) Bearing capacity assessment of distributed photovoltaic power distribution circuit
Taking a line a as an example to explain the carrying capacity evaluation process of the distributed photovoltaic power distribution line, the structure of the line a is shown in fig. 3, the load information of each node is shown in table 3, and the time sequence characteristic curve of the load is shown in fig. 4.
The following table shows impedance information of the line.
The following table shows the load information of each node of the line.
Table 3 load table for each node of line a
Node numbering | Maximum active power (kW) | Maximum reactive power (kVar) |
1 | 18 | 9 |
2 | 0 | 0 |
3 | 375 | 182 |
4 | 480 | 465 |
5 | 500 | 242 |
6 | 0 | 0 |
The lower graph shows the timing characteristics of the load.
And accessing the distributed photovoltaic power supply into a No. 3 node to obtain bearing capacity evaluation results under different confidence degrees, wherein the results are shown in the following table.
Table 4 results of bearing capacity at different confidence levels
Confidence (%) | 100% | 95% | 82% | 70% |
Bearing capacity (MW) | 10.395 | 12.385 | 12.847 | 16.084 |
It can be seen that the higher the confidence, the less load bearing. This is because, when the load bearing capacity is small, no voltage out-of-limit or line capacity out-of-limit occurs in different output scenarios, and the number of scenarios to which the load bearing capacity is adapted is large, and the confidence is high.
In this embodiment, the confidence level was chosen to be 95% and thus the load bearing capacity was 12.385MW.
Claims (6)
1. A method for evaluating distributed photovoltaic carrying capacity of a distribution line based on N-1 safety criteria, the method comprising the steps of:
step 1: combining uncertainty characteristics of source load, and considering acceptance capacity of a distributed photovoltaic power supply in a fault scene to establish a distribution line bearing capacity assessment method based on an N-1 safety criterion;
step 2: and carrying out decoupling dimension reduction calculation on the feasible domain in the N-1 safety criterion fault scene based on a decoupling strategy.
2. The method for evaluating the distributed photovoltaic carrying capacity of the distribution line based on the N-1 safety criterion according to claim 1, wherein in the step 1, the installation capacity of the distributed photovoltaic power supply is taken as an objective function, and the line capacity, the node voltage, the tide and the short-circuit current are taken as constraint conditions:
the objective function is that the distributed photovoltaic power supply of the line has the largest installation capacity, and the expression is as follows:
3. The method for evaluating the distributed photovoltaic load-bearing capacity of the distribution line based on the N-1 safety criterion according to claim 2, wherein the constraint condition of the distributed photovoltaic power load-bearing capacity evaluation model is as follows before a fault scene occurs:
1) Line capacity constraints
Wherein:the active power flowing through the ij branch t of the distribution line at the moment of the scene s is obtained; />The reactive power flowing through the ij branch t of the distribution line at the moment of the scene s is obtained; />Line capacity for ij branch of distribution line;
2) Tidal current constraint
V s,j,t =V s,i,t -(r ij P s,ij,t +x ij Q s,ij,t )
Wherein: p (P) s,ij,t The active power flowing through the branch ij at the moment t of the distribution line under the scene s is obtained; q (Q) s,ij,t The reactive power flowing through the branch ij at the moment t of the distribution line in the scene s is obtained; p (P) s,jk,t Active power flowing out of the node j and flowing to the node k at the moment of the distribution line t in the scene s; q (Q) s,jk,t The reactive power flows out of the node j and flows to the node k at the moment of the distribution line t in the scene s;the active power is the load active power of a node j at the moment t of the distribution line under the scene s; />The reactive power is the load reactive power of a node j at the moment t of the distribution line under the scene s; />DG power of a node j at t moment of a distribution line in a scene s; />The energy storage and charging power of a node j at the moment t of the distribution line in the scene s; />The energy storage and discharge power of a node j at the moment t of the distribution line in the scene s is obtained; r is (r) ij The resistance of branch ij; x is x ij Is the reactance of branch ij.
4. The method for evaluating the distributed photovoltaic load-bearing capacity of the distribution line based on the N-1 safety criterion according to claim 2, wherein after a fault scene occurs, the constraint condition of the distributed photovoltaic power load-bearing capacity evaluation model is as follows:
1) Line capacity constraints
For the transmission line, after a fault occurs or an overhaul scene is met, the condition that the N-1 check is met means that: when one power transmission line fails or overhauls and stops running, the load of the power transmission line can be transferred to other power transmission lines which normally run, and the short-time allowable overload capacity of the normal running line is not more than the rated capacity of a single line, and the power transmission line is based on the N-1 safety constraint as shown in the formula:
|P(t)≤(N-1)S|
wherein: n is the number of transmission lines; s is the rated capacity of a single power transmission line; p (t) is the net load of the transmission line at the moment t;
2) Tidal current constraint
V f,j,t =V f,i,t -(r ij P f,ij,t +x ij Q f,ij,t )
Wherein: p (P) f,ij,t The active power flowing through the branch ij at the moment t of the distribution line in the fault scene f; q (Q) f,ij,t Reactive power flowing through a branch ij at the moment t of the distribution line in the fault scene f; p (P) f,jk,t Active power flowing out of the node j and flowing to the node k at the moment of the distribution line t in the fault scene f; q (Q) f,jk,t The reactive power flows out of the node j and flows to the node k at the moment of the distribution line t in the fault scene f;the active power is the load active power of a node j at the moment t of the distribution line under the fault scene f; />The reactive power is the load reactive power of a node j at the moment t of the distribution line in the fault scene f; />The DG power of a node j at the moment t of the distribution line in the fault scene f; />The energy storage and charging power of a node j at the moment t of the distribution line in the fault scene f; />The energy storage and discharge power of a node j at the moment t of the distribution line in the fault scene f; r is (r) ij The resistance of branch ij; x is x ij Is the reactance of branch ij.
5. The method for evaluating the distributed photovoltaic bearing capacity of the distribution line based on the N-1 safety criterion according to claim 1, wherein in the step 2, the feasible regions meeting the constraint in the ith line fault scene and the constraint in the non-fault scene are obtained in parallel, and the feasible regions considering the N-1 safety constraint can be characterized as the intersection of the feasible regions in all fault scenes, and the specific flow is as follows:
(1) According to the maximum value and the minimum value in each dimension, constructing an initialized distributed photovoltaic access capacity feasible domain;
the minimum solution problem for the ith dimension:
S.t.
maximum solving problem for the i-th dimension:
S.t.
the obtained information under the non-fault scene and the fault sceneIs used as vertex V to construct a feasible domain omega c ;
(2) New vertex search: by shifting out the feasible region omega c Searching for new vertices; let the inequality expression of the kth boundary plane beObtaining a new vertex by solving the following linear programming problem:
wherein: v (V) p Is the feasible domain omega c Any point in the above;
the optimal solution of the linear programming is the top point of the constraint. Therefore, the optimal solution is the new vertex, and after all sides are translated, all the vertices are marked as V new ;
(3) Algorithm termination condition checking: calculating the difference between the feasible region volumes before and after the new vertex in the consideration (2) to be DeltaV; the smaller the delta V change, the smaller the feasible region change before and after adding the new vertex, thus, when delta V is smaller than the given threshold, the algorithm is terminated and the feasible region omega under the c-th line fault is recorded c To construct a feasible region Ω from the new vertices otherwise c Then return to (2);
when obtaining the feasible region Ω in all fault conditions c The feasible region Ω taking into account the N-1 security constraints is characterized as the feasible region Ω under all fault conditions c And (3) obtaining the carrying capacity evaluation result of the distributed photovoltaic power distribution line.
6. The method for evaluating distributed photovoltaic load-bearing capacity of power distribution lines based on N-1 safety criteria according to claim 5, wherein in step 2, confidence levels are sorted from large to small, the confidence levels are opposite to the obtained load-bearing capacity arrangement order, and after the confidence levels are selected, the corresponding distributed photovoltaic power distribution line load-bearing capacity is determined.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310285711.7A CN116316887A (en) | 2023-03-22 | 2023-03-22 | Distribution line distributed photovoltaic bearing capacity assessment method based on N-1 safety criterion |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310285711.7A CN116316887A (en) | 2023-03-22 | 2023-03-22 | Distribution line distributed photovoltaic bearing capacity assessment method based on N-1 safety criterion |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116316887A true CN116316887A (en) | 2023-06-23 |
Family
ID=86828478
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310285711.7A Pending CN116316887A (en) | 2023-03-22 | 2023-03-22 | Distribution line distributed photovoltaic bearing capacity assessment method based on N-1 safety criterion |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116316887A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118074240A (en) * | 2024-04-22 | 2024-05-24 | 国网浙江省电力有限公司宁波供电公司 | Maximum admittance capacity assessment method and equipment for distributed power supply |
CN118157219A (en) * | 2024-01-31 | 2024-06-07 | 广东电网有限责任公司 | Photovoltaic access method considering feeder line absorption space and carbon reduction benefits |
-
2023
- 2023-03-22 CN CN202310285711.7A patent/CN116316887A/en active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118157219A (en) * | 2024-01-31 | 2024-06-07 | 广东电网有限责任公司 | Photovoltaic access method considering feeder line absorption space and carbon reduction benefits |
CN118074240A (en) * | 2024-04-22 | 2024-05-24 | 国网浙江省电力有限公司宁波供电公司 | Maximum admittance capacity assessment method and equipment for distributed power supply |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116316887A (en) | Distribution line distributed photovoltaic bearing capacity assessment method based on N-1 safety criterion | |
Ye et al. | Optimization of hybrid energy storage system control strategy for pure electric vehicle based on typical driving cycle | |
CN109818369B (en) | Distributed power supply planning method considering output fuzzy randomness | |
CN107316113A (en) | A kind of Transmission Expansion Planning in Electric method and system | |
CN109034587B (en) | Active power distribution system optimal scheduling method for coordinating multiple controllable units | |
CN107947183B (en) | Power distribution network self-adaptive optimization method containing three-terminal SNOP (single-input single-output) based on differential evolution | |
Timmers et al. | A systematic review of DC wind farm collector cost-effectiveness | |
Penangsang et al. | Optimal placement and sizing of distributed generation in radial distribution system using K-means clustering method | |
Rezaeian‐Marjani et al. | Probabilistic assessment of D‐STATCOM operation considering correlated uncertain variables | |
Ebrahimi et al. | Probabilistic optimal planning in active distribution networks considering non‐linear loads based on data clustering method | |
Ali et al. | Optimal planning of uncertain renewable energy sources in unbalanced distribution systems by a multi‐objective hybrid PSO–SCO algorithm | |
CN116683518A (en) | Comprehensive evaluation method for power distribution network considering low-carbon benefits | |
CN116757406A (en) | Distributed power supply access distribution network site selection and volume determination method considering uncertainty | |
CN116522800A (en) | Multi-target wind-power storage station site selection and volume determination method based on scene probability | |
CN114243766B (en) | Regional multi-energy system optimal configuration method and system | |
CN111697607A (en) | Multi-terminal flexible direct-current transmission receiving-end power grid access method and system | |
CN115587763A (en) | Distributed power supply bearing capacity evaluation method considering multi-voltage-level coupling | |
CN106549412B (en) | For wind farm grid-connected power supply performance evaluation method | |
Su et al. | Cluster division and optimal scheduling of offshore wind power with energy storage based on fast unfolding algorithm | |
CN111092453A (en) | Power grid key line identification method for multipoint access photovoltaic power supply | |
CN118333474B (en) | Method for solving capacity planning problem of island comprehensive energy system | |
Wang et al. | Equivalent Modeling of PV Cluster Dynamic Process Based on K-medoids Clustering | |
Wu et al. | A distribution network flexible resource capacity configuration method with large renewable energy sources access | |
CN114039350B (en) | Energy storage point selection method and system for improving new energy station group digestion capacity | |
CN116599162B (en) | Method for determining new energy permeability under N-1 |
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 |