CN104935013A - Method for calculating voltage distribution of feed line containing DG - Google Patents

Method for calculating voltage distribution of feed line containing DG Download PDF

Info

Publication number
CN104935013A
CN104935013A CN201510416929.7A CN201510416929A CN104935013A CN 104935013 A CN104935013 A CN 104935013A CN 201510416929 A CN201510416929 A CN 201510416929A CN 104935013 A CN104935013 A CN 104935013A
Authority
CN
China
Prior art keywords
sigma
node
interval
voltage
power source
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.)
Granted
Application number
CN201510416929.7A
Other languages
Chinese (zh)
Other versions
CN104935013B (en
Inventor
张新慧
梁永省
陈羽
徐丙垠
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University of Technology
Original Assignee
Shandong University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong University of Technology filed Critical Shandong University of Technology
Priority to CN201510416929.7A priority Critical patent/CN104935013B/en
Publication of CN104935013A publication Critical patent/CN104935013A/en
Application granted granted Critical
Publication of CN104935013B publication Critical patent/CN104935013B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The present invention relates to a method for calculating voltage distribution of a feed line containing a DG (Distributed Generation), characterized in that calculation of voltage distribution is performed through a superimposed principle, i.e., when a system power supply is in effect and when the DG is in effect, line loss is respectively calculated, then voltage of each node is calculated, and voltage distribution of the feed line is derived. When the plurality of arbitrary DGs are connected, the method also can derive the rule of voltage changes of each node of the feed line. Nowadays, grid connection capacity of the DG is bigger and bigger, the number of the DGs connected in the feed line is correspondingly increased. Therefore, the method can derive the voltage distribution case after the DGs with certain capacity are connected into the feed line, and can derive the optimum position and capacity of the DG when the plurality of DGs are connected into the feed line, thereby solving the problem of node voltage off-limit after connection of the DGs.

Description

Containing the computational methods that the feeder voltage of DG distributes
Technical field
The present invention relates to the computational methods that a kind of feeder voltage containing DG distributes, belong to distribution network technology field.
Background technology
Distributed power source is at energy field in occupation of consequence, and its development alleviates the deterioration process of scarcity of resources and problem of environmental pollution, becomes one of approach that nature regenerative resource effectively utilizes.Its impact on power distribution network is not considered in the access of early stage of development distributed power source, namely " namely connect and namely forget ", but along with the development of distributed generation technology, grid connection capacity is also increasing, original distribution network also becomes complicated, become the active power distribution network of two-way flow from the distribution of power one-way flow, the impact of power distribution network also can not be ignored.Power distribution network can be made to flow into harmonic wave after DG access, cause the quality of power supply defective, all can be able to impact the planning of power distribution network, relaying protection etc. in addition.
The access of following distributed power source will be not only list and power supply duplicate supply, but towards the future development of many plant-grid connection, circuit therefore will be made more complicated.Node voltage can be caused to raise after distributed power source access feeder line, even out-of-limit, rational capacity and on-position configure the voltage condition can improving feeder line, but irrational configuration easily causes feeder voltage change unbalanced and stablize, and even occurs the situation that some node voltage is out-of-limit.Therefore a kind of method is needed to the change in voltage rule of feeder line of deriving with derivation feeder voltage distribution situation.These computational methods can derive the access of feeder line any number of distributed power source time circuit voltage's distribiuting situation, can lay the foundation for the reasonable disposition of further analysis distribution formula position of source and capacity.
Summary of the invention
According to above deficiency of the prior art, the problem to be solved in the present invention is: provide a kind of simplicity of design, easy to use, touch manner can be used to control digital control system and data display, be convenient to secondary development upgrading, enhance robustness and the fail safe of touch control mode, improve the aobvious control equipment integrating based on Intelligent flat and the data communications method thereof of the data transmission efficiency in system.
The technical solution adopted for the present invention to solve the technical problems is:
The computational methods that a kind of feeder voltage containing DG distributes are provided, adopt the principle of superposition to carry out Voltage distribution computation, line loss when namely computing system power supply and distributed power source act on respectively, and then the voltage calculating each node, concrete step is:
After a, distributed power source access feeder line, circuit is divided into some sections, and top node and first distributed power source access point are first paragraph, and the n-th distributed power source access point and endpoint node are one section, and all the other are between every adjacent node be one section;
B, be continuous print function node voltage formula fitting;
C, respectively to the node voltage function differentiate between adjacent distributions formula power supply;
D, the rule of feeder line voltage's distribiuting when analyzing any number of distributed power source access situation.
The described following formula of step b interior joint voltage equation represents:
As k ∈ [0, d 1] time,
U K = U 0 - k Σ i = k N P i R + Σ i = k N Q i X U N + k Σ j = 1 n P d j R U N - - - ( 1 ) ;
As k ∈ [d x, d x+1] time, wherein d xrepresent the node at an xth DG place.
U K = U 0 - k Σ i = k N P i R + Σ i = k N Q i X U N + Σ j = 1 x P d j Rd i U N + k Σ j = x + 1 n P d j R U N - - - ( 2 ) ;
As k ∈ [d n, N] time,
U K = U 0 - k Σ i = k N P i R + Σ i = k N Q i X U N + Σ j = 1 n P d j Rd i U N - - - ( 3 ) ;
Suppose k value be (0, N] on any real number, then U kthe continuous function on interval can be defined as, and at interval (0, d 1), (d x, d x+1), (d n, N) namely except the whole story node and distributed power source access point except all can lead,
Now suppose d 0=0, represent top node, corresponding access capacity P d0=0; d n+1=N, represents endpoint node, corresponding access capacity P d (n+1)=0;
Then formula (4) can be merged in formula (1), (2), (3)
U K = U 0 - k Σ i = k N P i R + Σ i = k N Q i X U N + Σ j = 1 x P d j Rd i U N + k Σ j = x + 1 n P d j R U N , k ∈ [ d x , d x + 1 ] - - - ( 4 )
Wherein d xrepresent the node at an xth DG place, x is the arbitrary integer in [0, n] interval;
Now the span of k in formula (4) is expanded, make it get arbitrary real number, then U in [0, N] kbecome a continuous print piecewise function, and except arbitrary d xpoint function all can be led, then by solving the rule of feeder voltage distribution to the mode of piecewise function differentiate piecemeal, step is as follows:
1. first, input d 1, d 2d n, n, N, P d1, P d2p dnvalue and load value size, assignment x=0;
2. ask interval [d x, d x+1] i.e. interval [d 0, d 1] interior U kderived function
3. solve d respectively x, d x+1the right-hand derivative at place and left derivative, and whether the long-pending of both judgements is greater than 0;
If 3.1 both amassing are greater than 0, then judge whether the right-hand derivative of dx is greater than 0, if be greater than 0, export k at (d x, d x+1) increase progressively in interval, if be less than 0, export k at (d x, d x+1) successively decrease in interval;
If 3.2 both amassing are less than 0, then judge that k is at [d x+ 1, d x+1-1] whether exist in interval and a bit make be 0, if so, then export k at (d x, d x+1) first successively decreasing in interval increases progressively afterwards, if do not exist, judges U dx-U d (x+1)whether be greater than 0;
If be 3.2.1 greater than 0, then export k at (d x, d x+1) successively decrease in interval;
If be 3.2.2 less than 0, then export k at (d x, d x+1) increase progressively in interval;
4. judge whether x is greater than n, if be greater than n, terminate, if be less than n, then make x=x+1, repeat step and play 2, until x>n;
5. sum up the feeder voltage regularity of distribution.
The beneficial effect that the present invention has is:
When the access of any number of distributed power source, this algorithm can be derived the rule of each node voltage change of feeder line.Now, the grid connection capacity of distributed power source is increasing, and the number accessing distributed power source in feeder line also increases accordingly.Therefore voltage's distribiuting situation after the distributed power source access feeder line utilizing this algorithm can derive under certain capacity, and then the optimum position of DG and capacity during many plant-grid connection in feeder line of can deriving, solve the out-of-limit problem of node voltage after distributed power source access.
The method can Computation distribution formula plant-grid connection distribution feeder time stable state node voltage, and then derivation feed connection node change in voltage trend.When distributed power source accesses the quantity of feeder line and position determines, the rule of voltage's distribiuting during different situations of can deriving.Simultaneously because distributed power source access capacity gets more and more, access quantity is also different, therefore can by analyzing the Voltage Distribution of feeder line and then releasing best on-position and the capacity limit of distributed power source by this algorithm, to avoid node voltage out-of-limit, make line voltage distribution improve effect and reach best.
Accompanying drawing explanation
Fig. 1 is that n DG accesses feeder line schematic diagram;
Fig. 2 is voltage's distribiuting reasoning flow figure;
Embodiment
Below in conjunction with accompanying drawing, embodiments of the invention are described further:
As depicted in figs. 1 and 2, suppose have in n distributed power source access feeder line, the node of access is respectively d 1, d 2d nrepresent the 1st, 2 respectively ... the node of n DG access, access capacity is respectively P d1, P d2p dn.After deriving by analysis, the voltage equation of each node k can point situation represent with formula (1)-(3) respectively.
As k ∈ [0, d 1] time, U K = U 0 - k Σ i = k N P i R + Σ i = k N Q i X U N + k Σ j = 1 n P d j R U N - - - ( 1 )
As k ∈ [d x, d x+1] time, wherein d xrepresent the node at an xth DG place.
U K = U 0 - k Σ i = k N P i R + Σ i = k N Q i X U N + Σ j = 1 x P d j Rd i U N + k Σ j = x + 1 n P d j R U N - - - ( 2 )
As k ∈ [d n, N] time, U K = U 0 - k Σ i = k N P i R + Σ i = k N Q i X U N + Σ j = 1 n P d j Rd i U N - - - ( 3 )
Suppose k value be (0, N] on any real number, then U kthe continuous function on interval can be defined as, and at interval (0, d 1), (d x, d x+1), (d n, N) namely except the whole story node and distributed power source access point except all can lead.Therefore the derivative that finds a function of by stages, and then solve the change of feeder voltage in respective bins.
Now suppose d 0=0, represent top node, corresponding access capacity P d0=0; d n+1=N, represents endpoint node, corresponding access capacity P d (n+1)=0.
Then formula (4) can be merged in formula (1), (2), (3)
U K = U 0 - k Σ i = k N P i R + Σ i = k N Q i X U N + Σ j = 1 x P d j Rd i U N + k Σ j = x + 1 n P d j R U N , k ∈ [ d x , d x + 1 ] - - - ( 4 )
Wherein d xrepresent the node at an xth DG place, x is the arbitrary integer in [0, n] interval.
Now the span of k in formula (4) is expanded, make it get arbitrary real number, then U in [0, N] kbecome a continuous print piecewise function, and except arbitrary d xpoint function all can be led.Then by solving the rule of feeder voltage distribution to the mode of piecewise function differentiate piecemeal.Step is as follows.
1. first, input d 1, d 2d n, n, N, P d1, P d2p dnvalue and load value size, assignment x=0;
2. ask interval [d x, d x+1] i.e. interval [d 0, d 1] interior U kderived function
3. solve d respectively x, d x+1right-hand derivative and left derivative, and whether long-pending both judging is greater than 0;
If 3.1 both amassing are greater than 0, then judge d xright-hand derivative whether be greater than 0, if be greater than 0, export k at (d x, d x+1) increase progressively in interval, if be less than 0, export k at (d x, d x+1) successively decrease in interval;
If 3.2 both amassing are less than 0, then judge that k is at [d x+ 1, d x+1-1] whether exist in interval and a bit make be 0, if so, then export k at (d x, d x+1) first successively decreasing in interval increases progressively afterwards, if do not exist, judges U dx-U d (x+1)whether be greater than 0;
If be 3.2.1 greater than 0, then export k at (d x, d x+1) successively decrease in interval;
If be 3.2.2 less than 0, then export k at (d x, d x+1) increase progressively in interval;
4. judge whether x is greater than n, if be greater than n, terminate, if be less than n, then make x=x+1, repeat step and play 2, until x>n;
5. sum up the feeder voltage regularity of distribution.

Claims (2)

1. containing the computational methods that the feeder voltage of DG distributes, it is characterized in that, adopt the principle of superposition to carry out Voltage distribution computation, line loss when namely computing system power supply and distributed power source act on respectively, and then the voltage calculating each node, concrete step is:
After a, distributed power source access feeder line, circuit is divided into some sections, and top node and first distributed power source access point are first paragraph, and the n-th distributed power source access point and endpoint node are one section, and all the other are between every adjacent node be one section;
B, be continuous print function node voltage formula fitting;
C, respectively to the node voltage function differentiate between adjacent distributions formula power supply;
D, the rule of feeder line voltage's distribiuting when analyzing any number of distributed power source access situation.
2. the computational methods that distribute of the feeder voltage containing DG according to claim 1, it is characterized in that, the described following formula of step b interior joint voltage equation represents:
As k ∈ [0, d 1] time,
U K = U 0 - k Σ i = k N P i R + Σ i = k N Q i X U N + k Σ j = k N P d j R U N - - - ( 1 ) ;
As k ∈ [d x, d x+1] time, wherein d xrepresent the node at an xth DG place.
U K = U 0 - k Σ i = 1 N P i R + Σ i = 1 N Q i X U N + Σ j = 1 x P d j Rd i U N + k Σ j = x + 1 n P d j R U N - - - ( 2 ) ;
As k ∈ [d n, N] time,
U K = U 0 - k Σ i = k N P i R + Σ i = k N Q i X U N + Σ j = 1 n P d j Rd i U N - - - ( 3 ) ;
Suppose k value be (0, N] on any real number, then U kthe continuous function on interval can be defined as, and at interval (0, d 1), (d x, d x+1), (d n, N) namely except the whole story node and distributed power source access point except all can lead,
Now suppose d 0=0, represent top node, corresponding access capacity P d0=0; d n+1=N, represents endpoint node, corresponding access capacity P d (n+1)=0;
Then formula (4) can be merged in formula (1), (2), (3)
U K = U 0 - k Σ i = k N P i R + Σ i = k N Q i X U N + Σ j = 1 n P d j Rd i U N + k Σ j = x + 1 n P d j R U N , k ∈ [ d x , d x + 1 ] - - - ( 4 )
Wherein d xrepresent the node at an xth DG place, x is the arbitrary integer in [0, n] interval;
Now the span of k in formula (4) is expanded, make it get arbitrary real number, then U in [0, N] kbecome a continuous print piecewise function, and except arbitrary d xpoint function all can be led, then by solving the rule of feeder voltage distribution to the mode of piecewise function differentiate piecemeal, step is as follows:
(1). first, input d 1, d 2d n, n, N, P d1, P d2p dnvalue and load value size, assignment x=0;
(2). ask interval [d x, d x+1] i.e. interval [d 0, d 1] interior U kderived function
(3). solve respectively at d x, d x+1the right-hand derivative at place and left derivative, and whether the long-pending of both judgements is greater than 0;
If 3.1 both amassing are greater than 0, then judge d xright-hand derivative whether be greater than 0, if be greater than 0, export k at (d x, d x+1) increase progressively in interval, if be less than 0, export k at (d x, d x+1) successively decrease in interval;
If 3.2 both amassing are less than 0, then judge that k is at [d x+ 1, d x+1-1] whether exist in interval and a bit make be 0, if so, then export k at (d x, d x+1) first successively decreasing in interval increases progressively afterwards, if do not exist, judges U dx-U d (x+1)whether be greater than 0;
If be 3.2.1 greater than 0, then export k at (d x, d x+1) successively decrease in interval;
If be 3.2.2 less than 0, then export k at (d x, d x+1) increase progressively in interval;
(4). judge whether x is greater than n, if be greater than n, terminate, if be less than n, then make x=x+1, repeat step and play 2, until x>n;
(5). sum up the feeder voltage regularity of distribution.
CN201510416929.7A 2015-07-16 2015-07-16 The computational methods of feeder voltage distribution containing DG Expired - Fee Related CN104935013B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510416929.7A CN104935013B (en) 2015-07-16 2015-07-16 The computational methods of feeder voltage distribution containing DG

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510416929.7A CN104935013B (en) 2015-07-16 2015-07-16 The computational methods of feeder voltage distribution containing DG

Publications (2)

Publication Number Publication Date
CN104935013A true CN104935013A (en) 2015-09-23
CN104935013B CN104935013B (en) 2018-06-19

Family

ID=54122047

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510416929.7A Expired - Fee Related CN104935013B (en) 2015-07-16 2015-07-16 The computational methods of feeder voltage distribution containing DG

Country Status (1)

Country Link
CN (1) CN104935013B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107104432A (en) * 2017-05-12 2017-08-29 山东理工大学 A kind of computational methods of distributed power access point voltage
CN107196293A (en) * 2017-06-07 2017-09-22 国网上海市电力公司 The application of the steady state voltage distributed model of distributed photovoltaic multiple spot access distribution
CN112994102A (en) * 2021-03-19 2021-06-18 国网湖北省电力有限公司电力科学研究院 Voltage calculation optimization method for power distribution network feeder line access distributed power supply

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103971026A (en) * 2014-05-28 2014-08-06 南京工程学院 General method for calculating tide of positive power distribution networks
CN104730420A (en) * 2015-03-26 2015-06-24 曾江 Node voltage estimation and voltage sag rapid detection method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103971026A (en) * 2014-05-28 2014-08-06 南京工程学院 General method for calculating tide of positive power distribution networks
CN104730420A (en) * 2015-03-26 2015-06-24 曾江 Node voltage estimation and voltage sag rapid detection method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
G.X.LUO & A.SEMLYEN: "Efficient load flow for large weakly meshed networks", 《IEEE TRANSACTIONS ON POWER SYSTEMS》 *
钟浩等: "基于曲线拟合的泰勒级数电网故障后电压快速准确计算", 《湖南大学学报》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107104432A (en) * 2017-05-12 2017-08-29 山东理工大学 A kind of computational methods of distributed power access point voltage
CN107104432B (en) * 2017-05-12 2019-04-09 山东理工大学 A kind of calculation method of distributed power access point voltage
CN107196293A (en) * 2017-06-07 2017-09-22 国网上海市电力公司 The application of the steady state voltage distributed model of distributed photovoltaic multiple spot access distribution
CN107196293B (en) * 2017-06-07 2020-05-12 国网上海市电力公司 Application method of steady-state voltage distribution model of distributed photovoltaic multi-point access distribution network
CN112994102A (en) * 2021-03-19 2021-06-18 国网湖北省电力有限公司电力科学研究院 Voltage calculation optimization method for power distribution network feeder line access distributed power supply

Also Published As

Publication number Publication date
CN104935013B (en) 2018-06-19

Similar Documents

Publication Publication Date Title
Li et al. Distributed adaptive robust voltage/var control with network partition in active distribution networks
Gupta et al. Countrywide PV hosting capacity and energy storage requirements for distribution networks: The case of Switzerland
CN106451556A (en) Method and apparatus for determining capability of receiving distributed power supply by distributed network feeder line
CN103337864B (en) Many direct currents concentrate the dynamic passive compensation mounting points optimization method of feed-in receiving end electrical network
CN102593829B (en) Method and device for determining transmission section in consideration of key branch circuit
CN105449713A (en) Active power distribution network soft normally open point (SNOP) planning method considering distributed generator characteristics
CN106503839A (en) A kind of marine wind electric field annular current collection network hierarchy planing method
CN103544541A (en) Carbon emission reduction evaluation and calculation method for intelligent power distribution and utilization system
CN102570450B (en) Static reliability assessment method for complex power grid
CN105337275A (en) Medium-voltage power distribution network power supply capacity evaluation method based on reliability of power distribution system
CN105244889A (en) Power distribution network reactive power compensation planning method
CN109039766B (en) Power CPS network risk propagation threshold determination method based on seepage probability
CN108023364A (en) Power distribution network distributed generation resource maximum access capability computational methods based on convex difference planning
CN104077664B (en) Confidence capacity assessment method of energy storage and generation system of wind power
Peng et al. Reliability analysis in interdependent smart grid systems
CN103093391A (en) City power distribution network arranging method
CN105375461A (en) Active power distribution network power supply capacity real-time assessment method based on prediction technology
CN102163845B (en) Optimal configuration method of distributed generations (DG) based on power moment algorithm
CN104915724A (en) AC-DC parallel power transmission channel power optimization distribution method and system
CN103887792B (en) A kind of low-voltage distribution network modeling method containing distributed power source
CN108462210A (en) Photovoltaic based on data mining can open the computational methods of capacity
CN104821578A (en) Large-scale wind power-containing power transmission system planning method taking available transmission capacity into account
CN104935013A (en) Method for calculating voltage distribution of feed line containing DG
CN106845848A (en) The appraisal procedure and device of reciprocal effect between Multi-infeed HVDC transmission system current conversion station
CN107622360A (en) A kind of critical circuits recognition methods for considering subjective and objective factor

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180619

Termination date: 20190716

CF01 Termination of patent right due to non-payment of annual fee