CN105787604A - Transmission-distribution-cooperative distribution network situation quick awareness method - Google Patents
Transmission-distribution-cooperative distribution network situation quick awareness method Download PDFInfo
- Publication number
- CN105787604A CN105787604A CN201610165392.6A CN201610165392A CN105787604A CN 105787604 A CN105787604 A CN 105787604A CN 201610165392 A CN201610165392 A CN 201610165392A CN 105787604 A CN105787604 A CN 105787604A
- Authority
- CN
- China
- Prior art keywords
- distribution
- power
- transmission
- network
- node
- 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
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 251
- 230000005540 biological transmission Effects 0.000 title claims abstract description 107
- 238000000034 method Methods 0.000 title claims abstract description 105
- 238000012937 correction Methods 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 23
- 238000004422 calculation algorithm Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 3
- 230000002265 prevention Effects 0.000 abstract description 8
- 238000005516 engineering process Methods 0.000 abstract description 4
- 238000011217 control strategy Methods 0.000 abstract description 2
- 238000000354 decomposition reaction Methods 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 11
- 230000008447 perception Effects 0.000 description 9
- 230000008569 process Effects 0.000 description 8
- 230000000052 comparative effect Effects 0.000 description 6
- 238000011161 development Methods 0.000 description 6
- 238000005259 measurement Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 4
- 230000005611 electricity Effects 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000000205 computational method Methods 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 208000032170 Congenital Abnormalities Diseases 0.000 description 1
- 235000000332 black box Nutrition 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000013277 forecasting method Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000004899 motility Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 102000030592 phosphoserine aminotransferase Human genes 0.000 description 1
- 108010088694 phosphoserine aminotransferase Proteins 0.000 description 1
- 230000036544 posture Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
Abstract
The invention discloses a transmission-distribution-cooperative distribution network situation quick awareness method. The method comprises the following steps: enabling a distribution network whole body to be equivalent to PQ load connected to boundary nodes of a transmission and distribution network; calculating pre-estimated power transmission network state through a quick decomposition method; calculating Thevenin equivalent parameters of a power transmission network; inputting the Thevenin equivalent parameters to root nodes of the distribution network, and calculating correction loss through a push-forward back-drawing method; updating distribution network equivalent PQ node power information through the correction loss to obtain Thevenin equivalent parameters; and transmitting the updated Thevenin equivalent parameters to the root nodes of the distribution network, and calculating the state of the distribution network after power fluctuation through the push-forward back-drawing method. The method makes an important step for realization of real-time situation awareness technology, and can serve as a basis for studying various problems in a transmission and distribution collaborative framework, such as studying influence of transmission and distribution collaborative prevention and control on situation trend of a transmission and distribution collaborative power grid, thereby finding an optimum prevention and control strategy.
Description
Technical field
The present invention relates to the distribution situation method for quickly sensing that a kind of transmission & distribution are collaborative.
Background technology
Intelligent network distribution is as intelligent grid inalienable part, just towards possessing distribution network shelf structure flexible, reliable, efficient, high reliability and the communication network of high security, the distributed power source access of high permeability, the high-speed simulation of distribution system and the general direction of self-healing control and overall goal development.The continuing of distributed power generation and generation of electricity by new energy is incorporated to, and intelligent distribution network presents more flexible and changeable trend so that power transmission network is day by day tight with contacting of power distribution network, and therefore transmission & distribution Cooperative Analysis is to ensure that the important means of economic power system and safe operation.Comprehensive control of accurate anticipation with evolution for realizing the collaborative operation of power networks state development trend of transmission & distribution, in the urgent need to promoting the situational awareness of distribution under the collaborative background of transmission & distribution, establish technical foundation for the situation security evaluation under all kinds of uncertain scenes and prevention and control.
Electrical network Situation Awareness is the important technical grasping operation of power networks track, by in wide area space-time unique, to the collection of all kinds of key elements relating to operation of power networks change, understanding and prediction, make every effort to the security postures grasping electrical network accurately and effectively so that the safety management of electrical network is from passively becoming actively.Existing pertinent literature is that situational awareness techniques has sketched the contours detailed blueprint in the application aspect of following intelligence power transmission network, but the key technology related to does not provide and is discussed in detail, and needs to be furtherd investigate further.
Research about the collaborative electrical network situational awareness techniques of transmission & distribution is still few at present, existing method solves the problem that under the collaborative background of transmission & distribution, bound power mismatch and voltage mismatch very easily occur for transmission side and distribution side, needing principal and subordinate's iteration repeatedly yet with calculating process, amount of calculation is bigger.Development and construction along with intelligent network distribution, the uncertain factor that system is run increases, the uncertain number of scenes in future that need to carry out safety situation evaluation is also exponentially increased, the collaborative tidal current computing method of transmission & distribution dividing theory based on principal and subordinate there is also difficulty for the requirement of real time meeting the control of intelligent network distribution situation and control decision, in the urgent need to the method for distribution situation quick sensing under the collaborative background of research transmission & distribution.
Summary of the invention
The present invention is to solve the problems referred to above, propose the distribution situation method for quickly sensing that a kind of transmission & distribution are collaborative, this method is based on the collaborative distribution Situation Awareness new method of transmission & distribution of wide area measurement Thevenin's equivalence, this method avoid the repeatedly principal and subordinate's iterative process dividing theoretical method based on principal and subordinate, possesses good calculating speed, laying a good foundation for distribution real-time situation perception, simulation example demonstrates feasibility and the effectiveness of institute's extracting method.
To achieve these goals, the present invention adopts the following technical scheme that
The distribution situation method for quickly sensing that a kind of transmission & distribution are collaborative, comprises the following steps:
(1) by distribution entirety equivalence it is the PQ load connect on transmission and distribution network boundary node, calculates the power transmission network state estimated with algorithm quicksort, calculate power transmission network Thevenin's equivalence parameter；
(2) Thevenin's equivalence parameter is accessed distribution root node, push back band method before utilization and calculate correction loss；
(3) utilize the PQ node power information revising loss renewal distribution equivalence, update and obtain Thevenin's equivalence parameter；
(4) the Thevenin's equivalence parameter after renewal is delivered to the root node of distribution, pushes back before utilization and calculate the state after distribution power swing with method.
In described step (1), specifically include:
(1-1) observe from arbitrary load bus to system side, complicated power system is converted into the equivalent network that equivalent electromotive force is powered to this load through equivalent impedance；
(1-2) the load equiva lent impedance of this equivalent network is calculated, by origin node admittance matrix is modified obtaining the bus admittance matrix of system；
(1-3) open-circuit voltage during load impedance that each node of calculating disconnection goes out, obtains Thevenin's equivalence electromotive force, again the bus admittance matrix of update the system, the mathematical relationship between structure power transmission network flow state and Thevenin's equivalence parameter.
In described step (1), gather the voltage vector of each node, current vector, active power, reactive power and load impedance and calculate equiva lent impedance.
In described step (1-3), bus admittance matrix is by being modified obtaining to origin node admittance matrix, load impedance is integrated into system admittance matrix, revises the diagonal element of the admittance battle array corresponding with load bus, obtain revised bus admittance matrix.
In described step (2), the power of power distribution network includes two parts, and a part is the power that the load that each bus nodes connects consumes, and to the distribution containing distributed power source, is send the PQ node of load by distributed power source node processing；Another part is distribution power attenuation.
In described step (2), specifically include:
(2-1) loss is revised in definition, in order to substitute the accurate distribution loss estimating under scene；
(2-2) injecting power estimating the connected PQ node of power distribution network it is set to the distribution loss of current operating conditions and estimates Distribution Network Load Data net power sum, calculating power transmission network Thevenin's equivalence parameter；
(2-3) Thevenin's equivalence parameter is accessed distribution root node, then calculate distribution loss with front pushing back band method, obtain revising loss.
In described step (2), the collaborative electrical network of transmission & distribution is carried out primary information and can be met engineer applied precision distribution loss result alternately.
In described step (3), constitute, with correction loss and Distribution Network Load Data, the power that boundary node injects to power distribution network, update the PQ node power information that distribution is equivalent, regain Thevenin's equivalence parameter.
In described step (4), the Thevenin's equivalence parameter of renewal is accessed distribution root node and is equivalent in distribution, add a new node.
In described step (4), when pushing back band method before utilization, need to increase newly calculate push away power and two nodes between two nodes before branch road rewind voltage.
The invention have the benefit that
(1) present invention has stepped essential step to the realization of real-time situation cognition technology, can as the basis of all kinds of problems in research transmission & distribution collaborative framework, as studied the transmission & distribution collaborative prevention and control impact on the collaborative electrical network situation track of transmission & distribution, thus finding optimum Prevention and control strategy；
(2) present invention exerts oneself probability forecasting method in combinations with new forms of energy such as wind-powered electricity generations, analyzes the impact on the collaborative electrical network static voltage stability risk of transmission & distribution under new forms of energy diverse access point and access way, obtains rational new forms of energy access scheme；Can also by the situation of the collaborative electrical network N-1 scene of transmission & distribution be carried out quick calculating etc.；
(3) present invention is that the collaborative aspect such as power grid security analysis and prevention and control of transmission & distribution is laid a good foundation, and has broad application prospects.
Accompanying drawing explanation
Fig. 1 be the present invention based on Thevenin's equivalence distribution Situation Awareness schematic diagram；
Fig. 2 be the present invention Thevenin's equivalence after system schematic；
Fig. 3 is the distribution Situation Awareness schematic diagram of the present invention；
Fig. 4 is the IEEE9 node system topological diagram of the present invention；
The distribution network voltage amplitude that Fig. 5 is the present invention calculates Comparative result figure；
Fig. 6 is the distribution network voltage phase calculation Comparative result figure of the present invention；
The node 6 distribution network voltage amplitude that Fig. 7 is the present invention calculates Comparative result figure；
Fig. 8 is the node 6 distribution network voltage phase calculation Comparative result figure of the present invention；
The node 8 distribution network voltage amplitude that Fig. 9 is the present invention calculates Comparative result figure；
Figure 10 is the node 8 distribution network voltage phase calculation Comparative result figure of the present invention.
Detailed description of the invention:
Below in conjunction with accompanying drawing, the invention will be further described with embodiment.
The basic thought of the collaborative distribution situation quick sensing of 1 transmission & distribution
Electrical network Situation Awareness process is divided into situation key element to gather, real-time situation is understood, further trend 3 stages of prediction, wherein further trend prediction is that the development and change rule to electrical network situation is summarized and reasoning, and then the development tendency of prediction further trend, and the potential safety problem under following operation of power networks state is provided prevention and control countermeasure, it is the requirement of highest level in Situation Awareness.Generation of electricity by new energy is exerted oneself and is had randomness and intermittent feature, when the further trend of perception distribution, the probability that the generation of electricity by new energy state of exerting oneself need to be done multiple scene is estimated, in conjunction with regulation and control scenes such as distribution rack reconstruct, the number of scenes that the prediction of distribution further trend and security evaluation is encountered is huge, and the requirement of real-time in line computation is comparatively urgent.Therefore, the Situation Awareness speed improving the collaborative distribution of transmission & distribution is significant.
Perception to the collaborative distribution situation of transmission & distribution, method of currently mainly analyzing is based on the method that principal and subordinate divides theory.Power distribution network equivalence is PQ node when carrying out power transmission network state computation by the method, the value of PQ node is power distribution network to the mutual power information of power transmission network, initial value is generally the net load power of power distribution network, then calculate power transmission network state with algorithm quicksort and obtain the voltage of boundary node, information of voltage is delivered to the root node of power distribution network again, pushing back before utilization and calculate state of electric distribution network with method and update power distribution network to the mutual power information of power transmission network, so iterating repeatedly until restraining.The method computational accuracy is high, overcome traditional method transmission side and the problem of distribution side bound power mismatch and voltage mismatch, but need principal and subordinate's iteration repeatedly to can be only achieved convergence, causing that larger electrical network is calculated hourly velocity owes excellent, the real-time situation perception that principal and subordinate's splitting method is applied to the collaborative distribution of transmission & distribution still suffers from certain speed bottleneck.
It is black-box model that complicated electric power system can be processed by coupled single-ended mouth equivalence, and parameter carries out identification, and accurate identification result can ensure that the response under outside microvariations of the single port after equivalence is consistent with complicated electric power system response under same disturbance.Thus it is extended to the collaborative electrical network of transmission & distribution, in principal and subordinate's splitting method, power distribution network and the continuous interactive iteration information of power transmission network are the main causes producing speed bottle-neck, if being coupled single-ended mouth the power transmission network part equivalence in collaborative for transmission & distribution electrical network, again by the situation front pushed back with method perception power distribution network, transmission side and the interactive iteration of distribution side boundary node information can be saved, thus distribution Situation Awareness speed can be improved.Thus, the present invention proposes based on the collaborative distribution Situation Awareness method of the transmission & distribution that transmission system coupled single-ended mouth is equivalent.
For the outer net static equivalence method of single port, Thevenin's equivalence parameter identification is to study the class method comparatively ripe, achievement is more, totally can be divided into and measuring and wide area measurement two class based on local.The method measured based on local assumes that system is repeatedly to measure running status in the time constant, and this hypothesis is difficult to set up in actual applications, and easily occurs that parameter drift problem affects identification precision.Development along with PMU, wide area measurement system reaches its maturity, Thevenin's equivalence parameter identification method based on wide area measurement proposes in succession, such method is based on single section tidal current data identification Thevenin's equivalence parameter, avoid the congenital defect based on local method for measurement, considering the coupled relation of equivalent internal system, outside response can be simulated the complication system before equivalence to outside response by the single port after equivalence.
For this, the present invention adopts and accesses distribution after power transmission network carrying out Thevenin's equivalence from power distribution network and power transmission network boundary node, then by distribution traditional before push back band method distribution carried out situation calculating, as shown in Figure 1:
In Fig. 1, power transmission network and power distribution network represent with M and S respectively, and B is the boundary node of transmission and distribution network, are also the root nodes of power distribution network, EthAnd ZthThe Thevenin's equivalence electromotive force after Thevenin's equivalence and impedance is carried out for transmission system.
When utilizing Thevenin's equivalence to carry out the collaborative distribution Situation Awareness of transmission & distribution, first it is PQ node by power distribution network equivalence, the process repeatedly mutual for saving transmission & distribution information, require that the final result that this PQ node power information and principal and subordinate's splitting method calculate is close, therefore need first accurately to estimate power distribution network to the mutual power information of power transmission network.Then after power transmission network being carried out Thevenin's equivalence, the Thevenin's equivalence electromotive force obtain equivalence and impedance access the root node of distribution, as shown in Figure 1, the distribution situation with method perception is pushed back before directly adopting again, eliminate the process of the continuous interface boundary information of transmission and distribution network in principal and subordinate's splitting method, significantly improve the speed of distribution Situation Awareness under the collaborative background of transmission & distribution.
2 Thevenin's equivalence parameter Fast Identification Methods
Thevenin's equivalence parameter calculates core thinking and is based on wide area measurement, and the load under a single state section is equivalent to load impedance, and the admittance battle array of update the system is asked for the open-circuit voltage of node, and then obtained whole equivalent parameters based on nodal voltage equation.
Observing from arbitrary load bus to system side in certain moment, complicated power system can be regarded as equivalent electromotive force EthThrough equivalent impedance ZthTo the equivalent network that this load is powered, as shown in Figure 2.
Under this moment state section, calculated by following formula and load be equivalent to impedance:
In formula:Pi、Qi、ZLiRespectively the voltage vector of node i, current vector, active power, reactive power, load impedance.
Load is equivalent to the system after impedance, and its bus admittance matrix can by being modified obtaining to origin node admittance matrix.If the original admittance matrix of system is Y0, Y0Middle element is Y0ij, load impedance is integrated into system admittance matrix and only need to revise and load
The diagonal element of the admittance battle array that node is corresponding, if revised admittance matrix is Y', diagonal element is Y 'ii, correction formula is as follows:
R is the number of PQ node.When solving the Thevenin's equivalence parameter at load bus i place, first calculate open-circuit voltage and obtain equivalent electrical potential parameters.It calculates the load impedance only needing disconnected node i place, calculates its open-circuit voltage, is Thevenin's equivalence electromotive force.
After disconnecting the load impedance at i place, needing update the system bus admittance matrix again, if revising posterior nodal point admittance matrix is Y, element is Yij.Correction formula is as follows:
Except YiiThe element of outer Y and Y' is all identical.
When solving i-node open-circuit voltage, node voltage vectorIn, PV node and balance node voltageFor known quantity, PQ node voltage phasor is amount to be asked.After adopting node load equivalent impedance, node injection current vectorInIt is 0, according to nodal voltage equation,
The equation group being made up of r equation can be obtained, as shown in formula (4):
By in above formula containingContinuous item moves on on the right side of equation, and formula (4) can be rewritten into following matrix form:
AU=-B (5)
Utilize Gaussian elimination method solving equation (5), i-node place open-circuit voltage can be solvedIt is Thevenin's equivalence electromotive forceThen its equivalent impedance Z is asked for according to formula (6)thi。
In formula:Representing the complex power of node i, subscript * represents conjugation.So far, utilize the whole network state, can quickly obtain the Thevenin's equivalence parameter observed outside node.
To sum up, the mathematical relationship between power transmission network flow state and Thevenin's equivalence parameter can be constructed, represent with following formula:
(Eth,Zth)=f (P, Q, V, θ, Y) (7)
In formula: P, Q, V, θ, Y represent node active power, reactive power, voltage, phase angle and electrical network admittance matrix, E respectivelyth,ZthRepresenting Thevenin's equivalence electromotive force and impedance, f () represents the quantitative relationship between Thevenin's equivalence parameter and electrical network known parameters.Formula (7) illustrates, can Thevenin's equivalence parameter corresponding to this state of quick obtaining to the real-time tracking of power transmission network state.
Distribution situation method for quickly sensing under the collaborative background of 3 transmission & distribution
The accurate predictor method that 3.1 transmission and distribution network power are mutual
Uncertainty that in distribution, new forms of energy are exerted oneself and intermittence, bring huge challenge to the safe operation of the collaborative electrical network of transmission & distribution, for assessing distribution operation conditions under following uncertain scene decision-making Optimal regulation and control measure, needs the situation of accurate perception distribution.By the carried basic thought of the present invention, the Situation Awareness reliable premise of result precision of distribution is the accurate identification of Thevenin's equivalence parameter, and Thevenin's equivalence parameter accurate recognition requires there is power transmission network state section accurately, namely require between transmission and distribution network mutual power accurately.Therefore, the mutual power accurately estimating transmission and distribution network is most important to the precision of the collaborative distribution Situation Awareness result of transmission & distribution.
To power distribution network, its power is mainly made up of two parts, and a part is the power that the load that each bus nodes connects consumes, and to the distribution containing distributed power source, can be send the PQ node of load by distributed power source node processing；Another part is distribution power attenuation, is formulated as follows:
In formula:For the power that boundary node injects to power distribution network,For the net load power in power distribution network, for known quantity,Account for leading；For the power attenuation in power distribution network, it is unknown quantity,Middle accounting example is less.Therefore, to accurately estimate distribution to the mutual power information of power transmission network, need in accurate predictor (8)To power swing less estimate scene, distribution loss is close with distribution loss under current state, can the distribution loss under direct current state substitute, the impact of result precision is less, but for the scene that distribution power swing is bigger, carrying out substituting with current distribution loss again and will cause calculating error relatively greatly, the present invention proposes to be suitable for the mutual power predictor method of distribution power swing scene at this.
First loss is revised in definition, usesRepresent, revise loss and be used for substituting the accurate distribution loss estimating under scene, therefore it is required that its value and active lossDeviation is little, and needs calculating simple and quick.Based on above-mentioned consideration, the scene estimated is calculated the process revising loss as follows:
(1) injecting power estimating the connected PQ node of power distribution network is set to the distribution loss of current operating conditions and estimates Distribution Network Load Data net power sum, available algorithm quicksort carries out Load flow calculation and estimates power transmission network state, then calculates power transmission network Thevenin's equivalence parameter with formula (7)；
(2) Thevenin's equivalence parameter is accessed distribution root node, as it is shown in figure 1, again with front pushing back the calculating distribution loss of band method, this loss is correction loss
Owing to distribution loss accounting example in Distribution Network Load Data is less, carry out that transmission & distribution obtain alternately to revise loss less with the difference of active loss ratio in Distribution Network Load Data, therefore the collaborative electrical network of transmission & distribution is carried out primary information and can be met engineer applied precision distribution loss result alternately, such as formula (9):
With revising lossAnd Distribution Network Load DataConstituteNamely replace formula (8) to be calculated by formula (9), simplified principal and subordinate's splitting method repeatedly interaction, contribute to improving the perception velocities of the collaborative distribution situation of transmission & distribution under the premise ensureing precision.
Distribution situation method for quickly sensing step under the collaborative background of 3.2 transmission & distribution
Following distribution network systems will continually introduce new forms of energy, distributed power source and microgrid, run to distribution while controlling to bring initiative and motility, and the uncertainty of distribution running status is significantly enhanced, and frequent fluctuation is caused that running status is continually changing by distribution power.For distribution power swing, to the collaborative distribution situation quick sensing of the transmission & distribution after distribution power swing, method detailed is as follows:
(1) be the PQ load connect on transmission and distribution network boundary node by distribution entirety equivalence, if the distribution loss that its performance number is current operating conditions with estimate Distribution Network Load Data net power sum.Calculate the power transmission network state estimated with algorithm quicksort, then calculate power transmission network Thevenin's equivalence parameter E ' with formula (7)th、Z′th；
(2) by Thevenin's equivalence parameter E 'th、Z′thAccess distribution root node, as it is shown in figure 1, calculate correction loss with front pushing back band method
(3) update, by formula (9), the PQ node power information that distribution is equivalent, repeat step (1) and obtain Thevenin's equivalence parameter E "th、Z″th；
(4) by E "th、Z″thIt is delivered to the root node of distribution, pushes back before recycling and calculate the state after distribution power swing with method.
In said process, Thevenin's equivalence parameter is accessed distribution root node and is equivalent in distribution, add a new node, as shown in Figure 3:
In Fig. 3, j is first wife's net root node, and i accesses the bus that distribution is newly-increased after power transmission network equivalence.When pushing back band method before utilization, need to increase newly push away horse-power formula (10) and two nodes between node i and node j before branch road rewind voltage equation (11):
In formula, Pj、Qj、VjRespectively the injection of node j is meritorious, idle and voltage magnitude, ViFor the voltage magnitude of node i, subscript k and k+1 pushes back the kth with process time and k+1 iteration before representing respectively, longitudinal component that Δ V, δ V respectively voltage lands and cross stream component.
4, simulated example
Uncertainty and intermittence that new forms of energy are exerted oneself cause that the collaborative distribution running status of transmission & distribution is continually changing, and at present also cannot the exerting oneself of Accurate Prediction future time instance new forms of energy, the multiple possible scene that future time instance new forms of energy are exerted oneself can only be provided, assume that subsequent time uncertainty power supply possible number of scenes of exerting oneself is K, if the number of uncertain power supply is N in distribution, then the uncertain number of scenes exerted oneself of subsequent time electrical network is KN, it can be seen that, if traversal calculates the distribution state of each possible scene, to calculating, rate request is higher.For this, in simulating, verifying process, for single scene, with the contrast that principal and subordinate's splitting method is calculated speed and precision, require to significantly improve under premise calculating speed if meeting in precision, then can illustrate that institute of the present invention extracting method can be obviously improved Situation Awareness speed, for based on the fact that the optimization prevention and control of Situation Awareness lay the foundation.First contrast the result precision of institute's extracting method of the present invention and principal and subordinate's splitting method for the transmission & distribution cooperative system comprising single distribution and calculate speed, the effectiveness of checking institute extracting method, then verifies the precision of institute's extracting method for the transmission & distribution cooperative system comprising many distributions again and calculates speed.
This simulation calculation hardware platform is THINKPADW530 work station, and CPU is i7-3740QM, dominant frequency 2.7GHz, internal memory 8G, and software platform is MATLAB and PSAT workbox.
4.1 single distribution examples
Emulating with IEEE3 machine 9 node electrical transmission network systems collocation IEEE33 node power distribution net system, IEEE9 node system topology is shown in Fig. 4.System being done following change: IEEE9 node system trunk rack electric pressure is set to 121kV, with IEEE33 node distribution network systems substitute node 6 load, distribution network systems electric pressure is 35kV.
The load changing distribution carrys out simulated power fluctuation, and the Distribution Network Load Data power after fluctuation is shown in appendix A.Contrasting two kinds of computational methods and obtain the distribution running status that transmission & distribution are collaborative, method one is divide the computational methods of theory based on principal and subordinate, and result is standard value in this way；Method two is institute of the present invention extracting method.
Table 1 distribution voltage magnitude calculates error and velocity contrast
By table 1 and Fig. 5, Fig. 6 it can be seen that method two is basically identical with method one result of calculation, error is very little, wherein voltage magnitude maximum relative error only 1.5507e-4%, and the maximum absolute error of voltage phase angle is 0.0024%.But under identical calculations environment, method two is used time 0.53 second only, and 1.25 seconds method one used times.Wherein method one principal and subordinate's iteration 4 times.When load level is promoted to 2 times of initial level, principal and subordinate's iterations is 10 times, and now the calculating used time of method one is about 3.75s, and method two is fixing due to transmission and distribution network information interaction times, and the calculating time remains 0.53s.Therefore, method two keeps same precision level with method one result, calculates speed and is obviously improved, if the calculating of estimating of power transmission network state is adopted sensitivity method, can further improve the calculating speed of method two.
Distribution example more than 4.2
On the basis of a upper example system, with IEEE33 node distribution network systems substitute node 8 load.The load power of two distributions of change is simulated new forms of energy and is gone out fluctuation, and the Distribution Network Load Data power after fluctuation is shown in appendix A, compares two kinds of methods state of electric distribution network after Distribution Network Load Data power swing, sees Fig. 7.
Table 2 node 6 distribution network voltage amplitude calculates error and velocity contrast
Table 3 node 8 distribution network voltage amplitude calculates error and velocity contrast
Appendix A distribution parameter and each node power
From Fig. 7～10 and table 2,3, method two is basically identical with method one result of calculation, error is less, the voltage magnitude maximum relative error of two of which distribution respectively 0.002% and 0.0014%, the maximum absolute error of voltage phase angle respectively 0.0073% and 0.0047%, but under identical calculations environment, method two is used time 0.8 second only, and 1.83 seconds method one used times, calculate speed and be obviously improved.
Contrast two examples, the calculating time of single game scape is increased along with the expansion of electrical network scale, it is contemplated that, when for when comprising the actual electric network system of thousands of node, the various uncertain scene composition of value volume and range of product will cause calculating huge, institute of the present invention extracting method is compared with principal and subordinate's splitting method the advantage calculated in speed becoming apparent from embodiment, and the realization to the collaborative distribution situation real-time perception of transmission & distribution has stepped essential step.
The specific embodiment of the present invention is described in conjunction with accompanying drawing although above-mentioned; but not limiting the scope of the invention; one of ordinary skill in the art should be understood that; on the basis of technical scheme, those skilled in the art need not pay various amendments or deformation that creative work can make still within protection scope of the present invention.
Claims (10)
1. the distribution situation method for quickly sensing that transmission & distribution are collaborative, is characterized in that: comprise the following steps:
(1) by distribution entirety equivalence it is the PQ load connect on transmission and distribution network boundary node, calculates the power transmission network state estimated with algorithm quicksort, calculate power transmission network Thevenin's equivalence parameter；
(2) Thevenin's equivalence parameter is accessed distribution root node, push back band method before utilization and calculate correction loss；
(3) utilize the PQ node power information revising loss renewal distribution equivalence, update and obtain Thevenin's equivalence parameter；
(4) the Thevenin's equivalence parameter after renewal is delivered to the root node of distribution, pushes back before utilization and calculate the state after distribution power swing with method.
2. the distribution situation method for quickly sensing that a kind of transmission & distribution as claimed in claim 1 are collaborative, is characterized in that: in described step (1), specifically include:
(1-1) observe from arbitrary load bus to system side, complicated power system is converted into the equivalent network that equivalent electromotive force is powered to this load through equivalent impedance；
(1-2) the load equiva lent impedance of this equivalent network is calculated, by origin node admittance matrix is modified obtaining the bus admittance matrix of system；
(1-3) open-circuit voltage during load impedance that each node of calculating disconnection goes out, obtains Thevenin's equivalence electromotive force, again the bus admittance matrix of update the system, the mathematical relationship between structure power transmission network flow state and Thevenin's equivalence parameter.
3. the distribution situation method for quickly sensing that a kind of transmission & distribution as claimed in claim 1 are collaborative, is characterized in that: in described step (1), gathers the voltage vector of each node, current vector, active power, reactive power and load impedance and calculates equiva lent impedance.
4. the distribution situation method for quickly sensing that a kind of transmission & distribution as claimed in claim 2 are collaborative, it is characterized in that: in described step (1-3), bus admittance matrix is by being modified obtaining to origin node admittance matrix, load impedance is integrated into system admittance matrix, revise the diagonal element of the admittance battle array corresponding with load bus, obtain revised bus admittance matrix.
5. the distribution situation method for quickly sensing that a kind of transmission & distribution as claimed in claim 1 are collaborative, it is characterized in that: in described step (2), the power of power distribution network includes two parts, a part is the power that the load that each bus nodes connects consumes, to the distribution containing distributed power source, it is send the PQ node of load by distributed power source node processing；Another part is distribution power attenuation.
6. the distribution situation method for quickly sensing that a kind of transmission & distribution as claimed in claim 1 are collaborative, is characterized in that: in described step (2), specifically include:
(2-1) loss is revised in definition, in order to substitute the accurate distribution loss estimating under scene；
(2-2) injecting power estimating the connected PQ node of power distribution network it is set to the distribution loss of current operating conditions and estimates Distribution Network Load Data net power sum, calculating power transmission network Thevenin's equivalence parameter；
(2-3) Thevenin's equivalence parameter is accessed distribution root node, then calculate distribution loss with front pushing back band method, obtain revising loss.
7. the distribution situation method for quickly sensing that a kind of transmission & distribution as claimed in claim 1 are collaborative, is characterized in that: in described step (2), the collaborative electrical network of transmission & distribution is carried out primary information and can be met engineer applied precision distribution loss result alternately.
8. the distribution situation method for quickly sensing that a kind of transmission & distribution as claimed in claim 1 are collaborative, it is characterized in that: in described step (3), the power that boundary node injects is constituted to power distribution network with correction loss and Distribution Network Load Data, update the PQ node power information that distribution is equivalent, regain Thevenin's equivalence parameter.
9. the distribution situation method for quickly sensing that a kind of transmission & distribution as claimed in claim 1 are collaborative, is characterized in that: in described step (4), the Thevenin's equivalence parameter of renewal is accessed distribution root node and is equivalent to add a new node in distribution.
10. the distribution situation method for quickly sensing that a kind of transmission & distribution as claimed in claim 1 are collaborative, is characterized in that: in described step (4), when pushing back band method before utilization, need to increase newly calculate push away power and two nodes between two nodes before branch road rewind voltage.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610165392.6A CN105787604B (en) | 2016-03-22 | 2016-03-22 | A kind of distribution situation method for quickly sensing of transmission & distribution collaboration |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610165392.6A CN105787604B (en) | 2016-03-22 | 2016-03-22 | A kind of distribution situation method for quickly sensing of transmission & distribution collaboration |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105787604A true CN105787604A (en) | 2016-07-20 |
| CN105787604B CN105787604B (en) | 2019-05-14 |
Family
ID=56390533
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610165392.6A Expired - Fee Related CN105787604B (en) | 2016-03-22 | 2016-03-22 | A kind of distribution situation method for quickly sensing of transmission & distribution collaboration |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105787604B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110120673A (en) * | 2019-05-09 | 2019-08-13 | 山东大学 | Distributed transmission & distribution collaboration idle work optimization method and system based on Thevenin's equivalence parameter identification |
| CN111463794A (en) * | 2020-04-29 | 2020-07-28 | 山东大学 | Transmission and distribution network cooperative power flow control method and system with additional earth impedance |
| CN111563623A (en) * | 2020-04-30 | 2020-08-21 | 国网山东省电力公司威海供电公司 | Typical scene extraction method and system for wind power system planning |
| CN114113898A (en) * | 2021-11-29 | 2022-03-01 | 大连海事大学 | Power distribution network loss analysis method and system based on multi-source measurement data |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101505061A (en) * | 2009-03-17 | 2009-08-12 | 中国电力科学研究院 | Computation method capable of tracking Davinan equivalence parameter base on time domain simulation |
| CN102622488A (en) * | 2012-03-21 | 2012-08-01 | 江西省电力科学研究院 | Distributed power capacity planning method for distribution network |
| CN103632046A (en) * | 2013-11-20 | 2014-03-12 | 国网河南省电力公司南阳供电公司 | Power grid load flow calculation method |
| CN104778367A (en) * | 2015-04-20 | 2015-07-15 | 山东大学 | On-line wide-area Thevenin equivalent parameter calculation method based on single state section |
-
2016
- 2016-03-22 CN CN201610165392.6A patent/CN105787604B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101505061A (en) * | 2009-03-17 | 2009-08-12 | 中国电力科学研究院 | Computation method capable of tracking Davinan equivalence parameter base on time domain simulation |
| CN102622488A (en) * | 2012-03-21 | 2012-08-01 | 江西省电力科学研究院 | Distributed power capacity planning method for distribution network |
| CN103632046A (en) * | 2013-11-20 | 2014-03-12 | 国网河南省电力公司南阳供电公司 | Power grid load flow calculation method |
| CN104778367A (en) * | 2015-04-20 | 2015-07-15 | 山东大学 | On-line wide-area Thevenin equivalent parameter calculation method based on single state section |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110120673A (en) * | 2019-05-09 | 2019-08-13 | 山东大学 | Distributed transmission & distribution collaboration idle work optimization method and system based on Thevenin's equivalence parameter identification |
| CN110120673B (en) * | 2019-05-09 | 2020-10-02 | 山东大学 | Distributed transmission and distribution cooperative reactive power optimization method and system based on Thevenin equivalent parameter identification |
| CN111463794A (en) * | 2020-04-29 | 2020-07-28 | 山东大学 | Transmission and distribution network cooperative power flow control method and system with additional earth impedance |
| CN111463794B (en) * | 2020-04-29 | 2021-07-30 | 山东大学 | Transmission and distribution network cooperative power flow control method and system with additional earth impedance |
| CN111563623A (en) * | 2020-04-30 | 2020-08-21 | 国网山东省电力公司威海供电公司 | Typical scene extraction method and system for wind power system planning |
| CN111563623B (en) * | 2020-04-30 | 2022-05-10 | 国网山东省电力公司威海供电公司 | Typical scene extraction method and system for wind power system planning |
| CN114113898A (en) * | 2021-11-29 | 2022-03-01 | 大连海事大学 | Power distribution network loss analysis method and system based on multi-source measurement data |
| CN114113898B (en) * | 2021-11-29 | 2023-11-14 | 大连海事大学 | Power distribution network loss analysis method and system based on multi-source measurement data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105787604B (en) | 2019-05-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10158229B2 (en) | Distributed gradient descent for solving optimal power flow in radial networks | |
| CN101976835B (en) | Parallel computation method for Newton power flow of large-scale electric power system | |
| CN103700036B (en) | A kind of transient stability projecting integral method being suitable to power system Multiple Time Scales | |
| CN102354981B (en) | Distributed computation based voltage stability assessment method of sub-networks in interconnected power network | |
| CN103928925B (en) | Power distribution network load flow calculation method based on forward-backward sweep | |
| CN103810646B (en) | Improved projection integral algorithm based active power distribution system dynamic simulation method | |
| CN101527455B (en) | Interconnected electric network distributed current calculating method on the basis of alternation and iteration of current module | |
| CN105787604A (en) | Transmission-distribution-cooperative distribution network situation quick awareness method | |
| CN106383296A (en) | Improved impedance type active distribution network fault location algorithm based on phasor analysis | |
| CN105656031A (en) | Security risk assessment method of wind-power-included electric power system based on Gaussian mixture distribution characteristics | |
| CN105140919A (en) | Method for rapid topology analysis and establishment of topological island based on graph theory | |
| CN108388688A (en) | A kind of active distribution network reliability estimation method | |
| CN114091816A (en) | Power distribution network state estimation method based on gated graph neural network of data fusion | |
| CN104156542A (en) | Implicit-projection-based method for simulating stability of active power distribution system | |
| CN104113061B (en) | A kind of distribution network three-phase power flow method containing distributed power source | |
| CN109100614A (en) | A kind of transmission open acess system and method based on PMU device | |
| CN103474989A (en) | Network reconstruction method based on sensitivity analysis | |
| CN105205244A (en) | Closed loop operation simulation system based on electromechanical-electromagnetic hybrid simulation technology | |
| CN104037759B (en) | Electric power system injects the robust least-square estimation method of transfer distribution factor | |
| CN111079972A (en) | Method, device and medium for planning reliability of active power distribution network | |
| CN105811404B (en) | A kind of distribution static voltage stability situation monitoring method of transmission ＆ distribution collaboration | |
| CN106897942A (en) | A kind of power distribution network distributed parallel method for estimating state and device | |
| CN106208099A (en) | A kind of Method for Reactive Power Optimization in Power based on bi-level programming and application thereof | |
| CN107257130A (en) | The low-voltage network loss computing method of decoupling is measured based on region | |
| CN106372440B (en) | A kind of adaptive robust state estimation method of the power distribution network of parallel computation and device |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190514 |