CN115642597B - Distributed photovoltaic bearing capacity calculation method and device - Google Patents
Distributed photovoltaic bearing capacity calculation method and device Download PDFInfo
- Publication number
- CN115642597B CN115642597B CN202211661000.7A CN202211661000A CN115642597B CN 115642597 B CN115642597 B CN 115642597B CN 202211661000 A CN202211661000 A CN 202211661000A CN 115642597 B CN115642597 B CN 115642597B
- Authority
- CN
- China
- Prior art keywords
- distributed photovoltaic
- access
- capacity
- access scheme
- under
- 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
Images
Classifications
-
- 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
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
Abstract
The invention provides a distributed photovoltaic bearing capacity calculation method and a distributed photovoltaic bearing capacity calculation device, which relate to the technical field of electricity, and the method comprises the following steps: setting simulation initial values of a plurality of access nodes under each access scheme based on a plurality of access schemes of the distribution network access distributed photovoltaic; the access scheme comprises the number of access nodes and the positions of the access nodes; determining a plurality of simulated actual values of each access scheme after the total distributed photovoltaic capacity is increased by multiple times according to the simulated initial values of a plurality of access nodes under each access scheme; and calculating the distributed photovoltaic maximum bearing capacity of each access scheme under the constraint condition of the line carrying capacity of the power distribution network based on a plurality of simulated actual values of each access scheme. According to the distributed photovoltaic bearing capacity calculation method and device, the distributed photovoltaic maximum capacity corresponding to the distribution network node voltage reaching the upper limit is calculated through mechanism analysis of distribution network node voltage distribution, and the calculation efficiency of the distributed photovoltaic bearing capacity is effectively improved.
Description
Technical Field
The invention relates to the technical field of electricity, in particular to a distributed photovoltaic bearing capacity calculation method and device.
Background
With the increasing installed capacity of distributed photovoltaic systems, power distribution networks face multiple challenges: the excessive distributed photovoltaic output can cause the power flow return, the voltage lifting and even the out-of-limit of the power distribution network, and the safe operation of the power distribution network and the reliable power supply of users are seriously influenced. In addition, the phenomenon of light abandonment is serious due to overhigh permeability of the distributed photovoltaic, and great resource waste is caused. Therefore, it is necessary to evaluate the distributed photovoltaic carrying capacity of the power distribution network, so as to reasonably and effectively guide the distributed photovoltaic grid-connected scheme of the actual power distribution system.
At present, research methods for distributed photovoltaic bearing capacity evaluation of a power distribution network mainly comprise an analytical method, a mathematical optimization method and a random scene simulation method. Compared with the former two methods, the simulation scene of the random scene simulation method is not limited and has universality, the distributed photovoltaic carrying capacity of the actual power distribution network can be reflected more truly and accurately, and the problems of large calculation workload and low calculation efficiency exist at the same time.
In order to solve the problem of low computational efficiency of the random scene simulation method, the existing research mainly focuses on improving the sampling method so as to reduce the computational workload. For example, a sampling method such as normal distribution sampling, image Ellang distribution sampling, latin hypercube distribution sampling and the like is used for replacing the traditional uniform sampling, thereby reducing the invalid calculation amount. But the current research does not fully combine the mechanism of the overvoltage of the power distribution network to carry out optimization and improvement of the algorithm.
In view of the above, the present inventors have designed a distributed photovoltaic bearing capacity calculation method and apparatus through repeated experiments according to production design experiences in this field and related fields for many years, so as to solve the problems in the prior art.
Disclosure of Invention
The invention aims to provide a distributed photovoltaic bearing capacity calculation method and a distributed photovoltaic bearing capacity calculation device.
In order to achieve the above object, the present invention provides a distributed photovoltaic bearing capacity calculation method, wherein the method includes:
setting simulation initial values of a plurality of access nodes under each access scheme based on a plurality of access schemes for accessing a distribution network to distributed photovoltaic; wherein the access scheme comprises the number of access nodes and the position of the access nodes;
determining a plurality of simulated actual values of each access scheme after the total distributed photovoltaic capacity is increased for a plurality of times according to simulated initial values of a plurality of access nodes under each access scheme;
and calculating the distributed photovoltaic maximum bearing capacity of each access scheme under the constraint condition of the line carrying capacity of the power distribution network based on a plurality of simulated actual values of each access scheme.
The invention also provides a distributed photovoltaic bearing capacity calculation device, which comprises:
the system comprises a setting unit, a distribution network and a plurality of access schemes, wherein the setting unit is used for setting simulation initial values of a plurality of access nodes under each access scheme based on the plurality of access schemes for accessing the distribution network to the distributed photovoltaic; wherein the access scheme comprises the number of access nodes and the position of the access nodes;
the simulation calculation unit is used for determining a plurality of simulation actual values of each access scheme after the distributed photovoltaic total capacity is increased for a plurality of times according to the simulation initial values of the plurality of access nodes under each access scheme;
and the bearing capacity calculation unit is used for calculating the distributed photovoltaic maximum bearing capacity of each access scheme under the constraint condition of the line current-carrying capacity of the power distribution network on the basis of a plurality of simulated actual values of each access scheme.
The invention also proposes a computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method as described above when executing the computer program.
The invention also proposes a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program which, when executed by a processor, implements the method as described above.
The invention also proposes a computer program product, wherein the computer program product comprises a computer program which, when executed by a processor, implements the method as described above.
Compared with the prior art, the invention has the following characteristics and advantages:
compared with the traditional random scene simulation method, the method and the device for improving the calculation efficiency of the distributed photovoltaic bearing capacity do not perform load flow calculation on all scenes generated randomly, excessive iterative calculation is avoided, the calculation amount is effectively reduced on the premise of ensuring the accuracy, the calculation time is shortened, and the calculation efficiency is improved.
Drawings
The drawings described herein are for illustration purposes only and are not intended to limit the scope of the present disclosure in any way. In addition, the shapes, the proportional sizes, and the like of the respective members in the drawings are merely schematic for facilitating the understanding of the present invention, and do not specifically limit the shapes, the proportional sizes, and the like of the respective members of the present invention. Those skilled in the art, having the benefit of the teachings of this invention, may choose from the various possible shapes and proportional sizes to implement the invention as a matter of case.
Fig. 1 is a step diagram of a distributed photovoltaic bearing capacity calculation method according to the present invention.
FIG. 2 is a diagram of the steps of the present invention for determining a plurality of simulated actual values of an access scheme;
fig. 3 is a diagram illustrating a step of calculating the maximum load capacity of the distributed photovoltaic system according to the present invention.
Detailed Description
The details of the present invention can be more clearly understood in conjunction with the accompanying drawings and the description of the embodiments of the present invention. However, the specific embodiments of the present invention described herein are for the purpose of illustration only and should not be construed as limiting the invention in any way. Any possible variations based on the present invention may be conceived by the skilled person in the light of the teachings of the present invention, and these should be considered to fall within the scope of the present invention.
The invention provides a distributed photovoltaic bearing capacity calculation method, as shown in fig. 1 to 3, the method includes:
setting simulation initial values of a plurality of access nodes under each access scheme based on a plurality of access schemes for accessing the distribution network to the distributed photovoltaic; the access scheme comprises the number of access nodes and the positions of the access nodes;
determining a plurality of simulated actual values of each access scheme after the total distributed photovoltaic capacity is increased by multiple times according to the simulated initial values of a plurality of access nodes under each access scheme;
and calculating the distributed photovoltaic maximum bearing capacity of each access scheme under the constraint condition of the line carrying capacity of the power distribution network based on a plurality of simulated actual values of each access scheme.
Compared with the traditional random scene simulation method, the distributed photovoltaic bearing capacity calculation method provided by the invention has the advantages that load flow calculation is not carried out on all scenes generated randomly, excessive iterative calculation is avoided, the calculation amount is effectively reduced on the premise of ensuring the accuracy, the calculation time is shortened, and the calculation efficiency is improved.
In the invention, the maximum bearing capacity of the distributed photovoltaic means the maximum total capacity of the distributed photovoltaic which is accessed under the condition of meeting the safe operation of the power distribution network.
In an alternative embodiment of the invention, defining the access scheme comprises: a typical time profile is selected.
Specifically, annual time sequence power data of all loads in the power distribution network and distributed photovoltaic annual output time sequence data are obtained, and the time with the minimum ratio of total load power to total distributed photovoltaic output is selected as a typical time section.
In an alternative embodiment of the invention, defining the access scheme comprises: determining the number of access nodes of the distributed photovoltaic and determining the position of the access nodes of the distributed photovoltaic.
An alternative to this embodimentIn addition, the determining the number of access nodes of the distributed photovoltaic system includes: the total number of nodes of the power grid is set asNAll nodes except the feeder root node can be accessed to the distributed photovoltaic, and the probability of accessing to the distributed photovoltaic is the same. Sampling number of access nodes for generating distributed photovoltaic by using uniform distribution as principleN PV Is 1 to
Any integer in between.
In an optional example of this embodiment, determining the location of the access node of the distributed photovoltaic includes: the set of access nodes for setting and distributing power grid to access distributed photovoltaic
Wherein, in the step (A),Brepresenting a set of access nodes other than the feeder root node,B j is shown asjThe number principle of the access nodes is ' increasing from the head end to the tail end ' of the feeder line ', and the root node is not connected to the distributed photovoltaic cells, so that the distributed photovoltaic cells are integratedBThe middle node index is incremented from 2 toN. From the collectionBMedium random sample generationN pv Elements, resulting subsetsB PV I.e. a set, subset, of nodes accessing distributed photovoltaicsB PV The subscript of each element in the index table reflects the access node position of the distributed photovoltaic.
In an alternative example of this embodiment, the number of times of increasing the total capacity of the distributed photovoltaic system is set toMNext, the process is carried out.
In an optional embodiment of the present invention, if the access node of the distributed photovoltaic is a load node, the initial simulation value of the distributed photovoltaic access node is set as a load peak value; otherwise, setting the simulation initial value of the access node of the distributed photovoltaic as the mean value of all loads of the power distribution network.
Further, let usiUnder the individual access schemejThe initial simulation value of the distributed photovoltaic access capacity is
After increasing the number of times, the number is the sameiThe simulated initial value of each access scheme is
。
In an optional embodiment of the present invention, determining, according to initial simulation values of a plurality of access nodes under each access scheme, a plurality of actual simulation values of each access scheme after the total distributed photovoltaic capacity is increased for a plurality of times includes:
determining simulation actual values of a plurality of access nodes under each access scheme according to the simulation initial values of the access nodes under each access scheme and the incremental times of the total capacity of the distributed photovoltaic;
and summing the simulated actual values of the plurality of access nodes under each access scheme to obtain the simulated actual value of the distributed photovoltaic total capacity under each access scheme.
In an optional embodiment of the present invention, determining simulated actual values of a plurality of access nodes under each access scheme according to simulated initial values of the plurality of access nodes under each access scheme and the number of times of increasing the total distributed photovoltaic capacity includes:
setting a simulation incremental step length of a simulation initial value of each access node under each access scheme;
and determining the simulated actual value of each access node under each access scheme according to the incremental times and the simulated incremental step length of the total distributed photovoltaic capacity so as to obtain the simulated actual value sequence of the total distributed photovoltaic capacity of all the access schemes.
In an alternative example of this embodiment, the firstjThe simulated incremental step size of each access node isd j ， 
。
In an alternative example of this embodiment, the firstjDistributed photovoltaic is inxThe simulated actual value after the secondary increment is
， 
Then to the firstxAfter increasing by one time thisjThe simulated actual value of the access node is
。
In an optional embodiment of the present invention, determining, according to initial simulation values of a plurality of access nodes under each access scheme, a plurality of actual simulation values of each access scheme after the total distributed photovoltaic capacity is increased by multiple times, further includes:
determining the maximum actual voltage value of the access node under each access scheme through load flow calculation according to the increasing times of the total capacity of the distributed photovoltaic system so as to obtain the maximum actual voltage value sequence of the nodes of all the access schemes;
and setting the voltage data when the maximum voltage of the access node under each access scheme exceeds the voltage of the root node of the feeder line as the first element of the maximum actual voltage sequence of the access node.
Specifically, the first step is determined by load flow calculationxMaximum voltage of each access node of power distribution network during secondary increasing
Until the maximum voltage of each access node of the power distribution network exceeds the voltage of a feeder root node, and the simulation actual value of each access scheme of the distributed photovoltaic total capacity at the moment
Setting to access scheme simulated actual value sequenceP PV,i First element ofP PV,i(1) Setting the maximum voltage of the node under the scene as the maximum voltage sequenceU max,i First element ofU max,i(1) 。
In an optional example of this embodiment, before determining, by load flow calculation, a maximum voltage value of an access node under each access scheme according to the number of times of increasing the total capacity of the distributed photovoltaic system, the method further includes:
dividing the maximum voltage of each access node under each access scheme into the maximum voltage of the node at the first stage and the maximum voltage of the node at the second stage; and the maximum voltage of the access node of the first stage is equal to the voltage of the feeder root node of the distributed photovoltaic.
Specifically, as the total capacity of the distributed photovoltaic system increases, the maximum voltage of each access node of the power distribution network generally shows a two-stage variation trend. The maximum voltage of the access node of the first stage is equal to the voltage of the root node of the feederU 0 The maximum voltage of the access node at the second stage is increased along with the increase of the total capacity of the distributed photovoltaic system, and the maximum voltage of the access nodeU max,i With distributed total photovoltaic capacityP PV,i Approximately linear.
It should be noted that the existing mature technology of load flow calculation belongs to one of three basic calculations of an electric power system, and has an effect of calculating a voltage value of each node under the condition that a power grid topological structure and power of each node are known, and the load flow calculation can be calculated through electric power system simulation software such as matpower, openDSS and the like.
In an optional example, determining, by load flow calculation, a maximum voltage value of an access node under each access scheme according to the number of times of increasing the total capacity of the distributed photovoltaic system includes:
determining a node maximum voltage value under each access scheme based on a simulation actual value under each access scheme according to a linear relation between the maximum voltage of the access node and the total capacity of the distributed photovoltaic;
and setting the maximum voltage of the access node under each access scheme which does not exceed the voltage of the root node of the feeder line as other elements except the first element in the maximum actual voltage sequence of the access node.
Maximum voltage of access nodeU max,i With distributed total photovoltaic capacityP PV,i Approximately linear, can be approximately expressed as
Wherein, in the step (A),k i the maximum voltage variation quantity of the distribution network node caused by the increase of the unit capacity distributed photovoltaic is defined as the voltage increase rate,U b,i is defined as the virtual voltage intercept (vj) of the voltage,P PV,i for the simulated actual value sequence of each access scheme of the second-stage distributed photovoltaic,
， U max,i for the second stage maximum voltage sequence of each access node
Wherein, in the step (A),Tis the number of elements of the maximum voltage sequence of the node.
Is provided withP PV,i(t) To representP PV,i To middletThe number of the elements is one,U max,i(t) to representU max,i To middletThe number of the elements is one,
(Tsuggested 8).
When in usetWhen =1, can obtainP PV,i(t) =P PV,i(1) ， U max,i(t) =U max,i(1) . Fromt=2 starting, the total distributed photovoltaic capacity is in incremental stepsd j And (3) increasing the number:P x PV,i(t) =P PV,i(t-1) +0.1P 0 PV,i calculating to obtain the maximum voltage of the access node under the current scene by using power distribution network load flow calculation software and recording as a sequenceU max,i Of (2) element(s)U max,i(t) . If it is
Terminating the increment, sequenceP PV,i 、U max,i Completing the calculation; if not, then,tincrement by 1 (i.e.tIs updated tot+ 1) and repeatingP PV,i(t) 、U max,i(t) Is calculated until
。
Further, according to a linear relationship between the maximum node voltage and the total distributed photovoltaic capacity, the method for determining the maximum node voltage value under each access scheme based on the simulated actual value set under each access scheme includes:
and determining a voltage increase rate and a virtual voltage intercept according to the total distributed photovoltaic capacity and the maximum node voltage under each access scheme.
Further, the voltage increase rate is
(ii) a And the number of the first and second groups is,
wherein, in the process,P PV,i(t) sequence of distributed photovoltaic total capacities representing a second phaseP PV,i To middletThe number of the elements is one,
wherein, in the process,U max,i(t) node maximum voltage sequence representing the second phaseU max,i To middletAn element; when in uset≤TWhen the utility model is used, the water is discharged,Tis the number of elements of the maximum voltage sequence of the node.
Preferably, the virtual voltage intercept is
。
In an optional embodiment of the present invention, the constraint condition of the line ampacity of the power distribution network includes:
if maximum bearing capacity of distributed photovoltaic
If so, the maximum bearing capacity of the distributed photovoltaic can not meet the requirement of the line carrying capacity;
if the maximum bearing capacity of the distributed photovoltaic system
If so, the maximum bearing capacity of the distributed photovoltaic meets the current-carrying capacity requirement of the line;
wherein the content of the first and second substances,I SC indicating the maximum current that the line is allowed to pass,U Z is the reference voltage value of the feeder line.
In an optional example of the embodiment, calculating the distributed photovoltaic maximum bearing capacity of each access scheme under the constraint condition of the line ampacity of the power distribution network based on a plurality of simulated actual values of each access scheme includes:
and calculating the distributed photovoltaic bearing capacity under each access scheme under the constraint condition of the line current-carrying capacity of the power distribution network according to the voltage increase rate and the virtual voltage intercept.
Specifically, the firstiDistributed photovoltaic maximum bearing capacity under distributed photovoltaic access schemeP PV,i-max ， 
WhereinU up For the upper limit of the distribution network voltage allowed, for a 10kV distribution network,U up it is recommended to take 1.07p.u.s..
Further, repeating the calculating steps until the number of the access schemes of the distributed photovoltaic system reaches a preset valueM. The distributed photovoltaic maximum bearing capacity can be expressed as
WhereinP PV,i-max Is shown asiAnd the distributed photovoltaic maximum bearing capacity under the distributed photovoltaic access scheme.
In an optional embodiment of the present invention, the distributed photovoltaic bearing capacity calculation method further comprises: and drawing a scatter diagram of the maximum bearing capacity of the distributed photovoltaic system. And drawing a distributed photovoltaic maximum bearing capacity scatter diagram of the power distribution network by taking the distributed photovoltaic total capacity as a horizontal axis and the maximum voltage of the access node as a vertical axis so as to perform more visual analysis on the distributed photovoltaic bearing capacity of the power distribution network.
The distributed photovoltaic bearing capacity calculation method provided by the invention is based on the analysis of the maximum voltage of the node and the sensitivity of the distributed photovoltaic total capacity, load flow calculation is not carried out on all randomly generated scenes, the linear relation between the maximum voltage of the access node and the distributed photovoltaic total capacity in the second stage is fitted through a least square method, and the corresponding distributed photovoltaic total capacity is directly calculated when the voltage reaches the upper limit. The method effectively reduces the calculation amount on the premise of ensuring the accuracy, so that the calculation efficiency is improved.
The invention also provides a distributed photovoltaic bearing capacity calculation device, which comprises:
the setting unit is used for setting simulation initial values of a plurality of access nodes under each access scheme based on a plurality of access schemes for accessing the distribution network to the distributed photovoltaic; the access scheme comprises the number of access nodes and the positions of the access nodes;
the simulation calculation unit is used for determining a plurality of simulation actual values of each access scheme after the distributed photovoltaic total capacity is increased for a plurality of times according to the simulation initial values of a plurality of access nodes under each access scheme;
and the bearing capacity calculation unit is used for calculating the distributed photovoltaic maximum bearing capacity of each access scheme under the constraint condition of the line current-carrying capacity of the power distribution network on the basis of a plurality of simulated actual values of each access scheme.
The invention also provides a computer device, which includes a memory, a processor and a computer program stored in the memory and executable on the processor, and the processor implements the distributed photovoltaic bearing capacity calculation method as described above when executing the computer program.
The invention also proposes a computer-readable storage medium, in which a computer program is stored which, when being executed by a processor, implements the distributed photovoltaic load calculation method as described above.
The invention also proposes a computer program product comprising a computer program which, when executed by a processor, implements the distributed photovoltaic load calculation method as described above.
The specific implementation process of the distributed photovoltaic bearing capacity calculation method provided by the present invention is described in detail with reference to an embodiment:
step 1: a typical time profile is selected. Acquiring annual time sequence power data of all loads in the power distribution network and annual distributed photovoltaic output time sequence data, and selecting the moment with the minimum ratio of total load power to total distributed photovoltaic output as a typical time section.
And 2, step: determining the number of access nodes of the distributed photovoltaic. The total number of nodes of the power grid is set asNAll access nodes except the feeder root node can be accessed to the distributed photovoltaic, and the probability of accessing to the distributed photovoltaic is the same. Further, the number of access nodes for generating distributed photovoltaic is sampled on the basis of uniform distributionN PV ， N PV Is 1 to
Any integer therebetween.
And 3, step 3: and determining the position of each access node of the distributed photovoltaic. Node set of equipment power grid capable of accessing distributed photovoltaic
Wherein, in the step (A),Brepresenting a set of nodes other than the feeder root node,B j denotes the firstjThe numbering principle of each access node is that the number of each access node is increased from the head end to the tail end of the feeder line, and the root node is not accessed to the distributed photovoltaic system in the step 2, so that the distributed photovoltaic system is integratedBThe middle node subscript increments from 2 toN. From the setBMedium random sample generationN PV Element, resulting subsetB PV I.e. a set, subset, of access nodes to the distributed photovoltaicB PV The subscript of each element in the series reflects the distributed photovoltaicEach access node location.
And defining a distributed photovoltaic access scheme as the number of distributed photovoltaic access nodes and the positions of the access nodes generated in the step 2 and the step 3. The total simulation times of the distributed photovoltaic access scheme is set asMNext, the process is repeated. For the firstiAnd (4) executing the following steps 4-6.3 according to the distributed photovoltaic access scheme.
And 4, step 4: is provided with the firstiAnd (4) simulating an initial value of each distributed photovoltaic access node under each distributed photovoltaic access scheme. If the distributed photovoltaic grid-connected point is a load node, setting the simulation initial value of the distributed photovoltaic access node as a load peak value; otherwise, setting the simulation initial value of the access node of the distributed photovoltaic as the average value of all loads of the power distribution system. Wherein, is provided withiUnder the individual access schemejThe initial simulation value of the distributed photovoltaic access capacity is
And the simulation initial value of the total capacity of the distributed photovoltaic is
。
With the increase of the total capacity of the distributed photovoltaic system, the maximum voltage of each access node of the power distribution network generally shows a 2-stage change trend. The first stage is that the maximum voltage of each access node of the power distribution network is equal to the voltage of a root node of a feeder lineU 0 And in the second stage, the maximum voltage of the distribution network node is increased along with the increase of the total capacity of the distributed photovoltaic. In the second stage, the maximum voltage of the nodeU max,i With distributed total photovoltaic capacityP PV,i Approximately linear, can be approximately expressed as
Wherein, in the step (A),k i the maximum voltage variation quantity of the distribution network node caused by the increase of the unit capacity distributed photovoltaic is defined as the voltage increase rate,U b,i defined as the virtual voltage intercept.
And 5: calculating the distributed photovoltaic total capacity and the node maximum voltage corresponding to the node voltage of the feeder line root for the first time;
carrying out incremental increase and load flow calculation on the basis of the initial values of the distributed photovoltaic simulation, and setting the first stepjThe simulation increment step length of the distributed photovoltaic is
， 
. First, thejAn access node is inxThe actual value of the simulation after the secondary increment isP x PV,i(t) ， P x PV,i(t) =P PV,i(t-1) +0.1P 0 PV,i Then to the firstxThe simulated actual value of the distributed photovoltaic total capacity after the secondary increment is
. Determining through load flow calculationxMaximum voltage of distribution network node during secondary increment
Until the maximum voltage of the distribution network node exceeds the voltage of the root node of the feeder line, and the simulation actual value of the total capacity of the distributed photovoltaic at the moment
Is a sequence ofP PV,i First element ofP PV,i(1) Setting the maximum voltage of the node in the scene as a sequenceU max,i First element of (2)U max,i(1) Wherein, in the step (A),P PV,i for the second stage distributed photovoltaic total capacity sequence,
， U max,i for the second stage maximum voltage sequence of each access node
And T is the number of elements of the maximum voltage sequence of the node.
Step 6: calculate the firstiDistributed photovoltaic maximum bearing capacity under distributed photovoltaic access schemeP PV,i-max 。
Step 6.1: and calculating a distributed photovoltaic total capacity sequence and a node maximum voltage sequence of the second stage. Is provided withP PV,i(t) To representP PV,i( To middletThe number of the elements is one,U max,i(t) to representU max,i To middletThe number of the elements is one,
(8 is suggested).
When the temperature is higher than the set temperaturetIf =1, the result can be obtained from step 5P PV,i(t) =P PV,i(1) ， U max,i(t) =U max,i(1) FromtAnd =2, the total distributed photovoltaic capacity is increased in an incremental manner according to the step 5:P x PV,i(t) =P PV,i(t-1) +0.1P 0 PV,i calculating to obtain the maximum voltage of the node under the current scene by using power distribution network load flow calculation software and recording as a sequenceU max,i Of (2) element(s)U max,i(t) . If it is
Terminating the increment, sequence
、
Completing the calculation; if not, then,
increment by 1 (i.e. 1)tIs updated to
) And repeatingP PV,i(t) 、U max,i(t) Until calculating
。
Step 6.2: a voltage growth rate and a virtual voltage intercept are determined. The voltage increase rate is calculated by
Virtual voltage interceptU b,i Is calculated by the formula
Wherein
、
。
Step 6.3: calculating the maximum bearing capacity of the distributed photovoltaic under the ith distributed photovoltaic access schemeP PV,i-max ， 
WhereinU up For the upper limit of the distribution network voltage allowed, for a 10kV distribution network,U up it is recommended to take 1.07p.u..
And 7: repeating the step 2 to the step 6.3 until the number of the access schemes of the distributed photovoltaic reaches a preset valueM. The maximum bearing capacity of the distributed photovoltaic system can be expressed as
WhereinP PV,i-max Denotes the firstiAnd the distributed photovoltaic maximum bearing capacity under the distributed photovoltaic access scheme.
And 8: and checking whether the current-carrying capacity constraint condition of the line is met. If it is
If so, it indicates that the current maximum distributed photovoltaic bearing capacity does not meet the requirement of the line carrying capacity, and the maximum distributed photovoltaic bearing capacity is taken as the maximum line carrying capacity, that is, C =
(ii) a If it is
And if so, indicating that the current maximum bearing capacity of the distributed photovoltaic meets the current-carrying capacity requirement of the line. Wherein the content of the first and second substances,I SC indicating the maximum current that the line is allowed to pass,U Z is the reference voltage value of the feeder line.
And step 9: and drawing a scatter diagram of the maximum bearing capacity of the distributed photovoltaic system. And drawing a distributed photovoltaic maximum bearing capacity scatter diagram of the power distribution network by taking the distributed photovoltaic total capacity as a horizontal axis and the node maximum voltage as a vertical axis so as to carry out more visual analysis on the distributed photovoltaic bearing capacity of the power distribution network.
The present invention is not limited to the above embodiments, and in particular, various features described in different embodiments can be arbitrarily combined with each other to form other embodiments, and the features are understood to be applicable to any embodiment except the explicitly opposite descriptions, and are not limited to the described embodiments.
Claims (13)
1. A method of distributed photovoltaic load bearing calculation, the method comprising:
setting simulation initial values of a plurality of access nodes under each access scheme based on a plurality of access schemes for accessing a distribution network to distributed photovoltaic; wherein the access scheme comprises the number of access nodes and the position of the access nodes;
determining a plurality of simulated actual values of each access scheme after the distributed photovoltaic total capacity is increased for a plurality of times according to the simulated initial values of a plurality of access nodes under each access scheme; based on a plurality of simulated actual values of each access scheme, calculating the distributed photovoltaic maximum bearing capacity of each access scheme under the constraint condition of the line ampacity of the power distribution network;
the constraint condition of the line current-carrying capacity of the power distribution network comprises the following steps:
if the maximum bearing capacity of the distributed photovoltaic
If so, the maximum bearing capacity of the distributed photovoltaic can not meet the requirement of the line carrying capacity;
if the maximum bearing capacity of the distributed photovoltaic
If so, the maximum bearing capacity of the distributed photovoltaic meets the current-carrying capacity requirement of the line;
wherein the content of the first and second substances,I SC indicating the maximum current that the line is allowed to pass,U Z is the reference voltage value of the feeder line.
2. The method for calculating distributed photovoltaic bearing capacity according to claim 1, wherein determining a plurality of simulated actual values of each access scheme after a plurality of increments of total distributed photovoltaic capacity according to simulated initial values of a plurality of access nodes under each access scheme comprises:
determining simulation actual values of the plurality of access nodes under each access scheme according to the simulation initial values of the plurality of access nodes under each access scheme and the incremental times of the total distributed photovoltaic capacity;
and summing the simulated actual values of the plurality of access nodes under each access scheme to obtain the simulated actual value of the distributed photovoltaic total capacity under each access scheme.
3. The method for calculating distributed photovoltaic bearing capacity according to claim 2, wherein determining simulated actual values of the plurality of access nodes under each access scheme according to simulated initial values of the plurality of access nodes under each access scheme and the number of times of increasing the total capacity of the distributed photovoltaic system comprises:
setting the simulation incremental step length of the simulation initial value of each access node under each access scheme;
and determining the simulated actual value of each access node under each access scheme according to the incremental times of the distributed photovoltaic total capacity and the simulated incremental step length to obtain the simulated actual value sequence of the distributed photovoltaic total capacity of all the access schemes.
4. The method according to claim 2 or 3, wherein a plurality of simulated actual values of each access scheme after a plurality of increments of the total distributed photovoltaic capacity are determined according to simulated initial values of a plurality of access nodes under each access scheme, further comprising:
determining the maximum actual voltage value of the access node under each access scheme through load flow calculation according to the increasing times of the total distributed photovoltaic capacity so as to obtain the maximum actual voltage value sequence of the access node of all the access schemes;
and setting the voltage data of the maximum voltage of the access node under each access scheme to exceed the voltage of the root node of the feeder line as the first element of the maximum actual voltage sequence of the access node.
5. The distributed photovoltaic bearing capacity calculation method according to claim 4, wherein before determining the maximum voltage value of the access node under each access scheme through load flow calculation according to the number of times of increasing the total capacity of the distributed photovoltaic, the method further comprises:
dividing the maximum voltage of each access node under each access scheme into the maximum voltage of the node in the first stage and the maximum voltage of the node in the second stage; wherein the maximum voltage of the junction of the first stage is equal to the feeder root node voltage of the distributed photovoltaic.
6. The distributed photovoltaic bearing capacity calculation method according to claim 5, wherein determining the maximum voltage value of the access node under each access scheme through load flow calculation according to the number of times of increasing the total capacity of the distributed photovoltaic, comprises:
determining the maximum voltage value of the access node under each access scheme based on the simulated actual value under each access scheme according to the linear relation between the maximum voltage of the node and the total capacity of the distributed photovoltaic;
and setting the maximum voltage of the node under each access scheme which does not exceed the voltage of the root node of the feeder line as other elements except the head element in the maximum actual voltage sequence of the access node.
7. The distributed photovoltaic bearer calculation method of claim 6, wherein determining the access node maximum voltage value under each of the access schemes based on a set of simulated actual values under each of the access schemes according to a linear relationship between node maximum voltage and total distributed photovoltaic capacity comprises:
and determining a voltage increase rate and a virtual voltage intercept according to the total distributed photovoltaic capacity and the maximum voltage of the access node under each access scheme.
8. The distributed photovoltaic bearing capacity calculation method of claim 7, wherein the voltage increase rate is
(ii) a And the number of the first and second electrodes,
wherein, in the process, P PV,i(t) sequence of distributed photovoltaic total capacities representing a second phaseP PV,i To middletThe number of the elements is one,
PV,i for distributed photovoltaic total capacity sequencesP PV,i Average value of (2)；
Wherein, in the process,U max,i(t) node maximum voltage sequence representing the second phaseU max,i To middletThe number of the elements is one,
max,i is a sequence of maximum voltagesU max,i Average value of (d); when the temperature is higher than the set temperature
When the temperature of the water is higher than the set temperature,Tis the number of elements of the maximum voltage sequence of the node.
9. The distributed photovoltaic bearing capacity calculation method of claim 8, wherein the virtual voltage intercept is
。
10. The method for calculating distributed photovoltaic bearing capacity according to claim 7, wherein calculating the distributed photovoltaic maximum bearing capacity of each of the access schemes under a constraint condition of a line ampacity of a power distribution network based on a plurality of simulated actual values of each of the access schemes comprises:
and calculating the distributed photovoltaic bearing capacity under each access scheme under the constraint condition of the line current-carrying capacity of the power distribution network according to the voltage increase rate and the virtual voltage intercept.
11. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the distributed photovoltaic load calculation method of any of claims 1 to 10 when executing the computer program.
12. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program which, when executed by a processor, implements the distributed photovoltaic load calculation method of any of claims 1 to 10.
13. A computer program product, characterized in that it comprises a computer program which, when being executed by a processor, implements the distributed photovoltaic load bearing calculation method of any one of claims 1 to 10.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211661000.7A CN115642597B (en) | 2022-12-23 | 2022-12-23 | Distributed photovoltaic bearing capacity calculation method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211661000.7A CN115642597B (en) | 2022-12-23 | 2022-12-23 | Distributed photovoltaic bearing capacity calculation method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115642597A CN115642597A (en) | 2023-01-24 |
| CN115642597B true CN115642597B (en) | 2023-03-10 |
Family
ID=84949871
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211661000.7A Active CN115642597B (en) | 2022-12-23 | 2022-12-23 | Distributed photovoltaic bearing capacity calculation method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115642597B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112242703A (en) * | 2020-09-14 | 2021-01-19 | 河海大学 | Power distribution network photovoltaic consumption evaluation method based on PSO (particle swarm optimization) optimization Monte Carlo algorithm |
| CN113507120A (en) * | 2021-07-16 | 2021-10-15 | 国网山西省电力公司晋城供电公司 | Method and device for calculating bearing capacity and electronic equipment |
| CN113922411A (en) * | 2021-08-05 | 2022-01-11 | 国网冀北电力有限公司承德供电公司 | Method and device for determining distributed photovoltaic absorption capacity |
| CN114759615A (en) * | 2022-06-10 | 2022-07-15 | 南京邮电大学 | Distribution network distributed photovoltaic bearing capacity analysis and reduction method based on hybrid simulation |
| CN115455691A (en) * | 2022-09-09 | 2022-12-09 | 内蒙古电力（集团）有限责任公司内蒙古电力经济技术研究院分公司 | Power grid simulation method considering new energy power fluctuation and line dynamic current-carrying capacity |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20190356131A1 (en) * | 2018-05-21 | 2019-11-21 | Commonwealth Edison Company | Methods and systems for determining distribution grid hosting capacity |
-
2022
- 2022-12-23 CN CN202211661000.7A patent/CN115642597B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112242703A (en) * | 2020-09-14 | 2021-01-19 | 河海大学 | Power distribution network photovoltaic consumption evaluation method based on PSO (particle swarm optimization) optimization Monte Carlo algorithm |
| CN113507120A (en) * | 2021-07-16 | 2021-10-15 | 国网山西省电力公司晋城供电公司 | Method and device for calculating bearing capacity and electronic equipment |
| CN113922411A (en) * | 2021-08-05 | 2022-01-11 | 国网冀北电力有限公司承德供电公司 | Method and device for determining distributed photovoltaic absorption capacity |
| CN114759615A (en) * | 2022-06-10 | 2022-07-15 | 南京邮电大学 | Distribution network distributed photovoltaic bearing capacity analysis and reduction method based on hybrid simulation |
| CN115455691A (en) * | 2022-09-09 | 2022-12-09 | 内蒙古电力（集团）有限责任公司内蒙古电力经济技术研究院分公司 | Power grid simulation method considering new energy power fluctuation and line dynamic current-carrying capacity |
Non-Patent Citations (1)
| Title |
|---|
| 陈璨 等.分布式光伏边界渗透率快速定位及消纳方案择优.2022,第55卷(第55卷),摘要，正文第1节，图1-2. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115642597A (en) | 2023-01-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106532778B (en) | Method for calculating maximum access capacity of distributed photovoltaic grid connection | |
| CN109638870B (en) | Phase modulator configuration method of extra-high voltage direct current transmission end power grid | |
| CN109636009B (en) | Method and system for establishing neural network model for determining line loss of power grid | |
| CN110247397B (en) | Energy storage configuration method, system and device and readable storage medium | |
| CN107332252A (en) | A kind of power distribution network low-voltage comprehensive processing method of meter and broad sense reactive source | |
| Monteiro et al. | Long-term sizing of lead–acid batteries in order to reduce technical losses on distribution networks: A distributed generation approach | |
| CN115017854A (en) | Method for calculating maximum allowable capacity of DG (distributed generation) of power distribution network based on multidimensional evaluation index system | |
| CN105279297B (en) | A kind of PSS/E electric system equivalence methods based on engineer application | |
| CN110336322B (en) | Photovoltaic power generation access capacity determination method based on daily minimum load confidence interval | |
| CN107301499B (en) | Distribution feeder statistical line loss rate data cleaning method based on AMI data | |
| CN115642597B (en) | Distributed photovoltaic bearing capacity calculation method and device | |
| CN111342501B (en) | Reactive power control method for microgrid with distributed power supply | |
| Liu et al. | Hierarchical overvoltage predictive control scheme for a DFIG-based wind farm | |
| CN110323779B (en) | Method and system for dynamically aggregating power of distributed power generation and energy storage device | |
| CN105207255B (en) | A kind of power system peak regulation computational methods suitable for wind power output | |
| CN114285086B (en) | Electrochemical energy storage power station energy scheduling method and system based on improved ant colony algorithm | |
| CN111697607A (en) | Multi-terminal flexible direct-current transmission receiving-end power grid access method and system | |
| CN110829440B (en) | Three-phase medium-low voltage integrated power distribution network voltage control method and system | |
| Seepromting et al. | Optimal grid-connected with multi-solar PV placement and sizing for power loss reduction and voltage profile improvement | |
| CN108564249B (en) | Power distribution network confidence peak clipping benefit evaluation method considering distributed photovoltaic randomness | |
| Lamprianidou et al. | Impact of Data-Driven Modelling Approaches on the Analysis of Active Distribution Networks | |
| Zhou et al. | Research on transmission congestion of power system containing wind farms based on genetic algorithm | |
| Li et al. | Wind farms reactive power optimization using genetic/tabu hybrid algorithm | |
| CN112952869B (en) | Method and system for expanding and planning AC-DC hybrid system considering wind power access | |
| CN110445128B (en) | Sensitivity-based real-time evaluation method for new energy consumption capability of regional 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 |