CN105787604B - A kind of distribution situation method for quickly sensing of transmission & distribution collaboration - Google Patents
A kind of distribution situation method for quickly sensing of transmission & distribution collaboration Download PDFInfo
- Publication number
- CN105787604B CN105787604B CN201610165392.6A CN201610165392A CN105787604B CN 105787604 B CN105787604 B CN 105787604B CN 201610165392 A CN201610165392 A CN 201610165392A CN 105787604 B CN105787604 B CN 105787604B
- Authority
- CN
- China
- Prior art keywords
- distribution
- power
- transmission
- node
- thevenin
- 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.)
- Expired - Fee Related
Links
- 230000005540 biological transmission Effects 0.000 title claims abstract description 106
- 238000000034 method Methods 0.000 title claims abstract description 103
- 230000005611 electricity Effects 0.000 claims abstract description 16
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 23
- 230000002452 interceptive effect Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 3
- 230000002265 prevention Effects 0.000 abstract description 8
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000011160 research Methods 0.000 abstract description 3
- 230000019771 cognition Effects 0.000 abstract description 2
- 238000011217 control strategy Methods 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 12
- 230000003993 interaction Effects 0.000 description 9
- 230000008569 process Effects 0.000 description 9
- 238000011161 development Methods 0.000 description 7
- 230000008447 perception Effects 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 230000000052 comparative effect Effects 0.000 description 6
- 238000005259 measurement Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 235000013399 edible fruits Nutrition 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 230000010181 polygamy Effects 0.000 description 2
- 238000004088 simulation Methods 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
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 235000000332 black box Nutrition 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
- 238000007796 conventional method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000012850 discrimination method 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
- 230000006870 function Effects 0.000 description 1
- 235000008434 ginseng Nutrition 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
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification 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
- 238000004064 recycling Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000006467 substitution reaction Methods 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 kind of distribution situation method for quickly sensing of transmission & distribution collaboration, the following steps are included: by the whole equivalent PQ load by being connect on transmission and distribution network boundary node of distribution, the transmission of electricity net state estimated is calculated with algorithm quicksort, calculates power transmission network Thevenin's equivalence parameter；Thevenin's equivalence parameter is accessed into distribution root node, calculates amendment loss using the method that rewinds is pushed forward；The PQ node power information for updating distribution equivalence is lost using amendment, updates and obtains Thevenin's equivalence parameter；Updated Thevenin's equivalence parameter is transmitted to the root node of distribution, utilizes the state being pushed forward after rewinding method calculating distribution power swing.The present invention has stepped essential step to the realization of real-time situation cognition technology, can be used as the basis of all kinds of problems in research transmission & distribution collaborative framework, influence of the transmission & distribution collaboration prevention and control to transmission & distribution collaboration power grid situation track is such as studied, to find optimal Prevention and control strategy.
Description
Technical field
The present invention relates to a kind of distribution situation method for quickly sensing of transmission & distribution collaboration.
Background technique
Intelligent network distribution is as smart grid inalienable part, just towards having flexible, reliable, efficient power distribution network
Grid structure, the communication network of high reliability and high security, high permeability distributed generation resource access, distribution system it is quick
General direction and the overall goal development of emulation and self-healing control.Distributed power generation and generation of electricity by new energy are persistently incorporated to, intelligence
Can power distribution network more flexible changeable trend is presented so that contacting for power transmission network and power distribution network is increasingly close, therefore transmission & distribution cooperate with
Analysis is to guarantee the important means of economic power system and safe operation.To realize that transmission & distribution cooperate with operation of power networks state development trend
Comprehensive control of accurate anticipation and development process, the situational awareness of distribution under background is cooperateed with there is an urgent need to promote transmission & distribution, is
Situation security evaluation and prevention and control under all kinds of uncertain scenes establish technical foundation.
Power grid Situation Awareness is the important technical for grasping operation of power networks track, right by wide area space-time unique
It is related to acquisition, understanding and the prediction of all kinds of elements of operation of power networks variation, makes every effort to the security postures for accurately and effectively grasping power grid,
So that the safety management of power grid becomes actively from passive.Existing pertinent literature is situational awareness techniques in the following intelligent power transmission network
The application aspect key technology having sketched the contours detailed blueprint, but be related to do not provide and be discussed in detail, need further deeply to be ground
Study carefully.
Research at present about transmission & distribution collaboration power grid situational awareness techniques is still few, and existing method solves transmission & distribution collaboration back
The problem of bound power mismatch and voltage mismatch easily occur for transmission side and distribution side under scape, however since calculating process needs instead
Multiple principal and subordinate's iteration, calculation amount are larger.With the development and construction of intelligent network distribution, the uncertain factor of system operation increases, and needs
The future for carrying out safety situation evaluation does not know number of scenes and is exponentially increased yet, and divides theoretical transmission & distribution based on principal and subordinate and cooperates with tide
Flow calculation methodologies control for meeting intelligent network distribution situation and the requirement of real time of control decision there is also difficulties, there is an urgent need to study
The method that transmission & distribution cooperate with distribution situation quick sensing under background.
Summary of the invention
The present invention to solve the above-mentioned problems, proposes a kind of distribution situation method for quickly sensing of transmission & distribution collaboration, we
Method is the transmission & distribution collaboration distribution Situation Awareness new method based on wide area measurement Thevenin's equivalence, and this method avoid based on principal and subordinate point
The multiple principal and subordinate's iterative process for splitting theoretical method, has good calculating speed, lays a good foundation for the perception of distribution real-time situation,
Simulation example demonstrates the feasibility and validity of proposed method.
To achieve the goals above, the present invention adopts the following technical scheme:
A kind of distribution situation method for quickly sensing of transmission & distribution collaboration, comprising the following steps:
(1) it by the whole equivalent PQ load by being connect on transmission and distribution network boundary node of distribution, is calculated with algorithm quicksort pre-
The transmission of electricity net state estimated calculates power transmission network Thevenin's equivalence parameter；
(2) Thevenin's equivalence parameter is accessed into distribution root node, calculates amendment loss using the method that rewinds is pushed forward；
(3) the PQ node power information that distribution equivalence is updated using amendment loss, is updated and obtains Thevenin's equivalence parameter；
(4) updated Thevenin's equivalence parameter is transmitted to the root node of distribution, calculates distribution using the method that rewinds is pushed forward
State after power swing.
In the step (1), specifically include:
(1-1) from any load bus to from system side, by complicated electric system be converted into equivalent potential by etc.
The equivalent network that value impedance is powered to the load；
(1-2) calculates the load equivalent impedance of the equivalent network, is by being modified to origin node admittance matrix
The node admittance matrix of system；
(1-3) calculates open-circuit voltage when disconnecting the load impedance that each node goes out, and obtains Thevenin's equivalence potential, again
The node admittance matrix of update the system constructs the mathematical relationship between power transmission network flow state and Thevenin's equivalence parameter.
In the step (1), voltage vector, current vector, active power, reactive power and the load resistance of each node are acquired
Anti- calculating equivalent impedance.
In the step (1-3), node admittance matrix is hindered load by being modified to obtain to origin node admittance matrix
It is anti-to be integrated into system admittance matrix, the diagonal element of admittance battle array corresponding with load bus is modified, revised node is obtained
Admittance matrix.
In the step (2), the power of power distribution network includes two parts, and a part is the load consumption of each bus nodes connection
Power be to send the PQ node of load by distributed generation resource node processing to the distribution containing distributed generation resource；Another part is
Distribution power loss.
In the step (2), specifically include:
(2-1) definition amendment loss, to substitute the accurate distribution loss estimated under scene；
The distribution that the injecting power for estimating the connected PQ node of power distribution network is set as current operating conditions is lost and is estimated by (2-2)
The sum of Distribution Network Load Data net power calculates power transmission network Thevenin's equivalence parameter；
Thevenin's equivalence parameter is accessed distribution root node by (2-3), then calculates distribution loss with the method that rewinds is pushed forward, and is repaired
Positive loss.
In the step (2), interactive can be obtained of primary information is carried out to transmission & distribution collaboration power grid and meets engineer application precision
Result is lost in distribution.
In the step (3), the power that boundary node is injected to power distribution network is constituted with amendment loss and Distribution Network Load Data, is updated
The PQ node power information of distribution equivalence, regains Thevenin's equivalence parameter.
In the step (4), the Thevenin's equivalence parameter access distribution root node of update is equivalent to and is increased in distribution
One new node.
In the step (4), using be pushed forward rewind method when, need to increase newly and calculate branch between two nodes and be pushed forward power and two
Node rewinds voltage.
The invention has the benefit that
(1) present invention has stepped essential step to the realization of real-time situation cognition technology, can be used as research transmission & distribution collaboration frame
Influence of the transmission & distribution collaboration prevention and control to transmission & distribution collaboration power grid situation track is such as studied, to seek in the basis of all kinds of problems in frame
Look for optimal Prevention and control strategy；
(2) present invention analyzes new energy diverse access point and connects in combination with the new energy such as wind-powered electricity generation power output probability forecasting method
Enter the influence under mode to transmission & distribution collaboration power grid static voltage stability risk, obtains reasonable new energy access scheme；It can also be with
Situation by the way that transmission & distribution are cooperateed with power grid N-1 scene carries out quickly calculating etc.；
(3) present invention is that transmission & distribution cooperate with power grid security analysis and prevention and control etc. to lay a good foundation, and is had wide
Application prospect.
Detailed description of the invention
Fig. 1 is of the invention based on Thevenin's equivalence distribution Situation Awareness schematic diagram；
Fig. 2 is system schematic after Thevenin's equivalence of the invention；
Fig. 3 is distribution Situation Awareness schematic diagram of the invention；
Fig. 4 is 9 node system topological diagram of IEEE of the invention；
Fig. 5 is that distribution network voltage amplitude of the invention calculates comparative result figure；
Fig. 6 is distribution network voltage phase calculation comparative result figure of the invention；
Fig. 7 is that 6 distribution network voltage amplitude of node of the invention calculates comparative result figure；
Fig. 8 is 6 distribution network voltage phase calculation comparative result figure of node of the invention；
Fig. 9 is that 8 distribution network voltage amplitude of node of the invention calculates comparative result figure；
Figure 10 is 8 distribution network voltage phase calculation comparative result figure of node of the invention.
Specific embodiment:
The invention will be further described with embodiment with reference to the accompanying drawing.
The basic thought of 1 transmission & distribution collaboration distribution situation quick sensing
Power grid Situation Awareness process is divided into the acquisition of situation element, real-time situation understands, further trend predicts 3 stages,
The prediction of middle further trend is to summarize to the development and change rule of power grid situation and reasoning, and then predict the development of further trend
Variation tendency, and prevention and control countermeasure is provided to the potential safety problem under the following operation of power networks state, be in Situation Awareness most
High-level requirement.Generation of electricity by new energy power output has the characteristics that randomness and intermittent, when perceiving the further trend of distribution, need
The probability for doing multiple scenes to generation of electricity by new energy power output state is estimated, in conjunction with the regulation scene such as distribution rack reconstruct, to distribution
The number of scenes that the prediction of further trend and security evaluation are faced is huge, more compels to the requirement of real-time in line computation
It cuts.Therefore, the Situation Awareness speed for improving transmission & distribution collaboration distribution is significant.
Perception to transmission & distribution collaboration distribution situation, primary analysis method is that theoretical method is divided based on principal and subordinate at present.It should
Power distribution network equivalence is PQ node when carrying out power transmission network state computation by method, and the value of PQ node is that power distribution network is interactive to power transmission network
Power information, initial value is generally the net load power of power distribution network, then calculates transmission of electricity net state with algorithm quicksort and obtains
The voltage of boundary node, then information of voltage is transmitted to the root node of power distribution network, state of electric distribution network is calculated using the method that rewinds is pushed forward
And power distribution network is updated to the power information of power transmission network interaction, iteration is repeatedly until convergence repeatedly.This method computational accuracy is high,
The problem of overcoming conventional method transmission side and distribution side bound power mismatch and voltage mismatch, but need principal and subordinate's iteration repeatedly
Convergence can be reached, speed owes excellent when causing to calculate larger power grid, and principal and subordinate's splitting method is applied to transmission & distribution collaboration and is matched
There are still certain speed bottlenecks for the real-time situation perception of net.
Coupled single-ended mouth equivalence can be handled complicated electric power system as black-box model, recognized to parameter, accurately distinguished
Single port after knowing the certifiable equivalence of result is in the sound of response and complicated electric power system under same disturbance under external microvariations
It should be consistent.Thus it is extended to transmission & distribution collaboration power grid, power distribution network is with the continuous interactive iteration information of power transmission network in principal and subordinate's splitting method
The main reason for generating speed bottle-neck if being coupled single-ended mouth the transmission of electricity mesh portions equivalence in transmission & distribution collaboration power grid, then is used
It is pushed forward the situation for rewinding method perception power distribution network, the interactive iteration of transmission side and distribution side boundary node information can be saved, so as to
Improve distribution Situation Awareness speed.The present invention proposes that the transmission & distribution based on transmission system coupled single-ended mouth equivalence cooperate with distribution as a result,
Situation Awareness method.
For the outer net static equivalence method of single port, Thevenin's equivalence parameter identification be study more mature, achievement compared with
More a kind of methods can totally be divided into and be measured and two class of wide area measurement based on local.System is assumed based on the method that local measures
Operating status is constant in repeatedly measuring the time, this hypothesis is difficult to set up in practical applications, and easily parameter drift occurs and ask
Topic influences identification precision.With the development of PMU, wide area measurement system reaches its maturity, the Thevenin's equivalence ginseng based on wide area measurement
Number discrimination method proposes that such method is based on single section tidal current data identification Thevenin's equivalence parameter in succession, avoids and is based on
The inadequate natural endowment of local method for measurement considers the coupled relation of equivalent internal system, and the single port after equivalence is to external response
Complication system before analog is equivalent is to external response.
Match for this purpose, the present invention uses to access after carrying out power transmission network with power transmission network boundary node Thevenin's equivalence from power distribution network
Net, then situation calculating is carried out to distribution by distribution traditional method of rewinding that is pushed forward, as shown in Figure 1:
In Fig. 1, power transmission network and power distribution network are indicated with M and S respectively, and B is the boundary node and power distribution network of transmission and distribution network
Root node, EthAnd ZthThevenin's equivalence potential and impedance after carrying out Thevenin's equivalence for transmission system.
When carrying out transmission & distribution collaboration distribution Situation Awareness using Thevenin's equivalence, it is first PQ node by power distribution network equivalence, is
Save the process that transmission & distribution information repeatedly interacts, it is desirable that the final result that the PQ node power information is calculated with principal and subordinate's splitting method connects
Closely, it therefore needs accurately to estimate power information of the power distribution network to power transmission network interaction first.Then Thevenin's equivalence is carried out to power transmission network
Afterwards, the root node of Thevenin's equivalence potential equivalence obtained and impedance access distribution, is pushed forward as shown in Figure 1, directlying adopt again
The distribution situation for rewinding method perception, eliminates 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 transmission & distribution collaboration background.
2 Thevenin's equivalence parameter Fast Identification Methods
It is to be equivalent to the load under a single state section based on wide area measurement that Thevenin's equivalence parameter, which calculates core ideas,
Load impedance, the admittance battle array of update the system seek the open-circuit voltage of node based on nodal voltage equation, and then obtain all equivalent
Parameter.
At a certain moment from any load bus to from system side, complicated electric system can be regarded as equivalent potential Eth
By equivalent impedance ZthThe equivalent network powered to the load, as shown in Figure 2.
Under the moment state section, it is calculate by the following formula and load is equivalent to impedance:
In formula:Pi、Qi、ZLiRespectively the voltage vector, current vector of node i, active power, reactive power,
Load impedance.
System after being equivalent to impedance to load, node admittance matrix can be by repairing origin node admittance matrix
Just obtain.If the original admittance matrix of system is Y0, Y0Middle element is Y0ij, load impedance, which is integrated into system admittance matrix only, to be repaired
Change and load
The diagonal element of the corresponding admittance battle array of node, if revised admittance matrix is Y', diagonal element is Y 'ii, repair
Positive formula is as follows:
R is the number of PQ node.When solving the Thevenin's equivalence parameter at load bus i, open-circuit voltage is calculated first
Obtain equivalent electrical potential parameters.It calculates the load impedance only needed at disconnected node i, calculates its open-circuit voltage, as Dai Weinan etc.
It is worth potential.
After disconnecting the load impedance at i, update the system node admittance matrix again is needed, if amendment posterior nodal point admittance matrix is
Y, element 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 nodes
VoltageFor known quantity, PQ node voltage phasor is amount to be asked.After node load equivalent impedance, node injection electricity
Flow vectorInIt is 0, according to nodal voltage equation,
The available equation group being made of r equation, as shown in formula (4):
It will contain in above formulaContinuous item moves on on the right side of equation, and formula (4) can be rewritten into following matrix form:
AU=-B (5)
Equation (5) are solved using Gaussian elimination method, open-circuit voltage at i-node can be solvedAs Dai Weinan etc.
It is worth potentialThen its equivalent impedance Z is sought according to formula (6)thi。
In formula:Indicate that the complex power of node i, subscript * indicate conjugation.So far, using the whole network state, sight can be quickly obtained
Examine the Thevenin's equivalence parameter outside node.
To sum up, the mathematical relationship between power transmission network flow state and Thevenin's equivalence parameter can be constructed, is indicated with following formula:
(Eth,Zth)=f (P, Q, V, θ, Y) (7)
In formula: P, Q, V, θ, Y respectively indicate node active power, reactive power, voltage, phase angle and power grid admittance matrix,
Eth,ZthIndicate Thevenin's equivalence potential and impedance, f () indicates the amount between Thevenin's equivalence parameter and power grid known parameters
Change relationship.Formula (7) illustrates, can the corresponding Thevenin's equivalence parameter of the quick obtaining state to the real-time tracking of transmission of electricity net state.
3 transmission & distribution cooperate with distribution situation method for quickly sensing under background
The accurate predictor method of 3.1 transmission and distribution network power interaction
The uncertainty and intermittence that new energy is contributed in distribution are brought huge to the safe operation of transmission & distribution collaboration power grid
Challenge need to accurately perceive distribution for the distribution operation conditions under the following uncertain scene of assessment and decision Optimal regulation and control measure
Situation.Mentioned basic thought through the invention, the Situation Awareness result precision of distribution are reliable on condition that Thevenin's equivalence parameter
Accurate identification, and Thevenin's equivalence parameter accurate recognition require have accurately transmission of electricity net state section, that is, require transmission and distribution network
Between accurately interact power.Therefore, the interaction power of transmission and distribution network is accurately estimated to transmission & distribution collaboration distribution Situation Awareness result
Precision is most important.
To power distribution network, power is mainly made of two parts, and a part is the function of the load consumption of each bus nodes connection
Distributed generation resource node processing can be to send the PQ node of load to the distribution containing distributed generation resource by rate；Another part be with
Net power loss is formulated as follows:
In formula:For the power that boundary node is injected to power distribution network,It is known for the net load power in power distribution network
Amount,It accounts for leading；It is unknown quantity for the power loss in power distribution network,Middle accounting example is smaller.Therefore, accurate
Power information of the distribution to power transmission network interaction is estimated, is needed in accurate predictor (8)It is lesser to power swing to estimate field
Scape, distribution loss is close with distribution loss under current state, can directly be substituted with the distribution loss under current state, to knot
The influence of fruit precision is smaller, but is directed to the biggish scene of distribution power swing, then carries out substitution with current distribution loss and will lead to meter
Calculation error is larger, and the present invention herein proposes the interaction power predictor method of suitable distribution power swing scene.
Definition amendment loss first, is usedIt indicating, amendment loss is used to substitute the accurate distribution loss estimated under scene,
Therefore it is required that its value and active lossDeviation is little, and needs to calculate simple and quick.Based on above-mentioned consideration, to the field estimated
The process that scape calculates amendment loss is as follows:
(1) the distribution loss that the injecting power for estimating the connected PQ node of power distribution network is set as current operating conditions is matched with estimating
The sum of net load net power can carry out Load flow calculation with algorithm quicksort and estimate transmission of electricity net state, then be calculated and transmitted electricity with formula (7)
Net Thevenin's equivalence parameter；
(2) Thevenin's equivalence parameter is accessed into distribution root node, as shown in Figure 1, being calculated again with network loss with being pushed forward the method that rewinds
Consumption, the loss are to correct loss
Since distribution loss accounting example in Distribution Network Load Data is smaller, carry out amendment loss that the interaction of transmission & distribution obtains with it is real
Ratio of the difference of border loss in Distribution Network Load Data is smaller, therefore carries out primary information interaction to transmission & distribution collaboration power grid and can be obtained
Meet engineer application precision distribution loss as a result, such as formula (9):
It is lost with amendmentAnd Distribution Network Load DataIt constitutesIt replaces formula (8) to be calculated with formula (9), simplifies
The multiple interactive process of principal and subordinate's splitting method facilitates the perception speed that transmission & distribution collaboration distribution situation is improved under the premise of guaranteeing precision
Degree.
3.2 transmission & distribution cooperate with distribution situation method for quickly sensing step under background
The following distribution network systems will continually introduce new energy, distributed generation resource and microgrid, bring actively to distribution operation control
While property and flexibility, the uncertainty with Running State is significantly enhanced, and frequent fluctuation is caused to run by distribution power
State constantly changes.By taking distribution power swing as an example, distribution situation quick sensing is cooperateed with to the transmission & distribution after distribution power swing, in detail
Thin method is as follows:
(1) by the whole equivalent PQ load by being connect on transmission and distribution network boundary node of distribution, if its performance number is current fortune
The distribution of row state is lost and estimates the sum of Distribution Network Load Data net power.The transmission of electricity net state estimated is calculated with algorithm quicksort, then
Power transmission network Thevenin's equivalence parameter E ' is calculated with formula (7)th、Z′th；
(2) by Thevenin's equivalence parameter E 'th、Z′thDistribution root node is accessed, as shown in Figure 1, rewinding method meter with being pushed forward
Calculate amendment loss
(3) the PQ node power information that distribution equivalence is updated with formula (9) repeats step (1) and obtains Thevenin's equivalence parameter
E″th、Z″th；
(4) by E "th、Z″thIt is transmitted to the root node of distribution, recycling is pushed forward after the method that rewinds calculates distribution power swing
State.
In the above process, Thevenin's equivalence parameter access distribution root node is equivalent to increased in distribution one it is new
Node, as shown in Figure 3:
In Fig. 3, j is first wife's net root node, and i is the bus that access distribution is newly-increased after power transmission network is equivalent.It is pushed back before utilization
When band method, need to increase newly branch between node i and node j be pushed forward horse-power formula (10) and two nodes rewind voltage equation
(11):
In formula, Pj、Qj、VjThe injection of respectively node j is active and reactive and voltage magnitude, ViFor the voltage magnitude of node i,
Subscript k and k+1 respectively indicate the kth time and k+1 iteration for being pushed forward the process of rewinding, and Δ V, δ V are respectively the longitudinal direction point of voltage landing
Amount and cross stream component.
4, simulated example
The uncertainty and intermittence of new energy power output cause transmission & distribution collaboration constantly to change with Running State, and at present also
Can not Accurate Prediction future time instance new energy power output, a variety of possible scenes of future time instance new energy power output can only be provided, it is false
If subsequent time uncertainty power supply is contributed, possible number of scenes is K, next if the number for not knowing power supply in distribution is N
The number of scenes of the uncertain power output of moment power grid is KN, it can be seen that, it is right if traversal calculates the net state of matching of each possible scene
Calculating speed is more demanding.For this purpose, by taking single scene as an example, carrying out calculating speed with principal and subordinate's splitting method during simulating, verifying
The comparison of degree and precision can then illustrate the mentioned side of the present invention as can significantly improving calculating speed under the premise of precision is met the requirements
Method can be obviously improved Situation Awareness speed, be based on the fact that the optimization prevention and control of Situation Awareness lay the foundation.First to include list
The result precision and calculating speed that the present invention proposed method and principal and subordinate's splitting method is compared for the transmission & distribution cooperative system of distribution, are tested
The validity of proposed method is demonstrate,proved, then verify the precision of proposed method by taking the transmission & distribution cooperative system comprising polygamy net as an example again and is counted
Calculate speed.
This simulation calculation hardware platform is THINKPAD W530 work station, and CPU i7-3740QM, dominant frequency 2.7GHz are interior
8G is deposited, software platform is the tool box MATLAB and PSAT.
4.1 single distribution examples
It is emulated with 3 machine of IEEE, 9 node electrical transmission network systems collocation, 33 node power distribution net system of IEEE, 9 node of IEEE
System topological is shown in Fig. 4.Do following change to system: 9 node system trunk rack voltage class of IEEE is set as 121kV, uses IEEE
33 node distribution network systems substitute node, 6 load, distribution network systems voltage class are 35kV.
The load for changing distribution carrys out simulated power fluctuation, and the Distribution Network Load Data power after fluctuation is shown in appendix A.Compare two kinds of calculating
Method obtains the Running State of matching of transmission & distribution collaboration, and method one is that theoretical calculation method is divided based on principal and subordinate, ties in this way
Fruit is standard value；Method two is proposed method by the present invention.
1 distribution voltage magnitude of table calculates error and velocity contrast
By table 1 and Fig. 5, Fig. 6 it is found that method two and one calculated result of method are almost the same, error is very small, wherein voltage
Only 1.5507e-4%, the maximum absolute error of voltage phase angle are 0.0024% to amplitude maximum relative error.But in identical calculations ring
Under border, only used time 0.53 second of method two, and 1.25 seconds one used times of method.Wherein method one principal and subordinate's iteration 4 times.Work as load level
When being promoted to 2 times of initial level, principal and subordinate's the number of iterations is 10 times, and the calculating used time of method one is about 3.75s at this time, and method
Two due to transmission and distribution network information interaction times be it is fixed, calculate the time be still 0.53s.Therefore, method two and one knot of method
Fruit keeps same precision level, and calculating speed is obviously improved, if the calculating of estimating to transmission of electricity net state uses sensitivity method,
It can further improve the calculating speed of method two.
4.2 polygamy net examples
On the basis of upper example system, with 33 node distribution network systems substitute node of IEEE, 8 load.Change two to match
The load power of net goes out fluctuation to simulate new energy, and the Distribution Network Load Data power after fluctuation is shown in appendix A, compares two methods and is matching
State of electric distribution network after the fluctuation of net load power, is shown in Fig. 7.
2 node of table, 6 distribution network voltage amplitude calculates error and velocity contrast
3 node of table, 8 distribution network voltage amplitude calculates error and velocity contrast
Appendix A matches network parameters and each node power
By Fig. 7~10 and table 2,3 it is found that method two and one calculated result of method are almost the same, error is smaller, two of them
The voltage magnitude maximum relative error of distribution is respectively 0.002% and 0.0014%, the maximum absolute error difference of voltage phase angle
For 0.0073% and 0.0047%, but under identical calculations environment, only used time 0.8 second of method two, and 1.83 seconds one used times of method,
Calculating speed is obviously improved.
Two examples are compared, the calculating time of single game scape are increased with the expansion of power grid scale, it is contemplated that work as needle
When to actual electric network system comprising thousands of nodes, the numerous uncertain scene compositions of value volume and range of product will lead to calculation scale
It is huge, the mentioned method of the present invention compared with advantage of principal and subordinate's splitting method in calculating speed by the more obvious of embodiment, to transmission & distribution association
Realization with distribution situation real-time perception has stepped essential step.
Above-mentioned, although the foregoing specific embodiments of the present invention is described with reference to the accompanying drawings, not protects model to the present invention
The limitation enclosed, those skilled in the art should understand that, based on the technical solutions of the present invention, those skilled in the art are not
Need to make the creative labor the various modifications or changes that can be made still within protection scope of the present invention.
Claims (10)
1. a kind of distribution situation method for quickly sensing of transmission & distribution collaboration, it is characterized in that: the following steps are included:
(1) it by the whole equivalent PQ load by being connect on transmission and distribution network boundary node of distribution, is estimated with algorithm quicksort calculating
Transmission of electricity net state, calculates power transmission network Thevenin's equivalence parameter；
(2) Thevenin's equivalence parameter is accessed into distribution root node, calculates amendment loss using the method that rewinds is pushed forward；
(3) the PQ node power information that distribution equivalence is updated using amendment loss, is updated and obtains Thevenin's equivalence parameter；
(4) updated Thevenin's equivalence parameter is transmitted to the root node of distribution, calculates distribution power using the method that rewinds is pushed forward
State after fluctuation.
2. the distribution situation method for quickly sensing of a kind of transmission & distribution collaboration as described in claim 1, it is characterized in that: the step
(1) it in, specifically includes:
(1-1), to from system side, converts equivalent potential by equivalent resistance for complicated electric system from any load bus
The anti-equivalent network powered to the load；
(1-2) calculates the load equivalent impedance of the equivalent network, by being modified to obtain system to origin node admittance matrix
Node admittance matrix；
(1-3) calculates open-circuit voltage when disconnecting the load impedance at each node, obtains Thevenin's equivalence potential, corrects again
The node admittance matrix of system constructs the mathematical relationship between power transmission network flow state and Thevenin's equivalence parameter.
3. the distribution situation method for quickly sensing of a kind of transmission & distribution collaboration as described in claim 1, it is characterized in that: the step
(1) in, voltage vector, current vector, active power, reactive power and the load impedance for acquiring each node calculate equivalent impedance.
4. the distribution situation method for quickly sensing of a kind of transmission & distribution collaboration as claimed in claim 2, it is characterized in that: the step
In (1-3), node admittance matrix is integrated into system admittance by being modified to obtain to origin node admittance matrix, by load impedance
Matrix modifies the diagonal element of admittance battle array corresponding with load bus, obtains revised node admittance matrix.
5. the distribution situation method for quickly sensing of a kind of transmission & distribution collaboration as described in claim 1, it is characterized in that: the step
(2) in, the power of power distribution network includes two parts, and a part is the power of the load consumption of each bus nodes connection, to containing distribution
Distributed generation resource node processing is to send the PQ node of load by the distribution of formula power supply；Another part is distribution power loss.
6. the distribution situation method for quickly sensing of a kind of transmission & distribution collaboration as described in claim 1, it is characterized in that: the step
(2) it in, specifically includes:
(2-1) definition amendment loss, to substitute the accurate distribution loss estimated under scene；
The injecting power for estimating the connected PQ node of power distribution network is set as the distribution loss of current operating conditions and estimates distribution by (2-2)
The sum of load net power calculates power transmission network Thevenin's equivalence parameter；
Thevenin's equivalence parameter is accessed distribution root node by (2-3), then calculates distribution loss with the method that rewinds is pushed forward, and obtains amendment damage
Consumption.
7. the distribution situation method for quickly sensing of a kind of transmission & distribution collaboration as described in claim 1, it is characterized in that: the step
(2) in, interactive can be obtained of primary information is carried out to transmission & distribution collaboration power grid and meets engineer application precision distribution loss result.
8. the distribution situation method for quickly sensing of a kind of transmission & distribution collaboration as described in claim 1, it is characterized in that: the step
(3) in, the power that boundary node is injected to power distribution network is constituted with amendment loss and Distribution Network Load Data, updates the PQ node of distribution equivalence
Power information regains Thevenin's equivalence parameter.
9. the distribution situation method for quickly sensing of a kind of transmission & distribution collaboration as described in claim 1, it is characterized in that: the step
(4) in, the Thevenin's equivalence parameter access distribution root node of update is equivalent to and increases a new node in distribution.
10. the distribution situation method for quickly sensing of a kind of transmission & distribution collaboration as described in claim 1, it is characterized in that: the step
(4) in, using be pushed forward rewind method when, need to increase newly calculate branch between two nodes be pushed forward power and two nodes 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 CN105787604A (en) | 2016-07-20 |
| CN105787604B true 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) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| CN111463794B (en) * | 2020-04-29 | 2021-07-30 | 山东大学 | Transmission and distribution network cooperative power flow control method and system with additional earth impedance |
| CN111563623B (en) * | 2020-04-30 | 2022-05-10 | 国网山东省电力公司威海供电公司 | Typical scene extraction method and system for wind power system planning |
| CN114113898B (en) * | 2021-11-29 | 2023-11-14 | 大连海事大学 | 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 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105787604A (en) | 2016-07-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107453357B (en) | Power distribution network state estimation method based on layered solution | |
| CN105787604B (en) | A kind of distribution situation method for quickly sensing of transmission & distribution collaboration | |
| CN106443246B (en) | The on-line identification method of small interference stability parameter based on PMU metric data | |
| CN101958543B (en) | Method for estimating three-phase impedance-free nonlinear multi-source state of transformer substation | |
| CN106383296B (en) | A kind of improvement impedance type active power distribution network Fault Location Algorithm based on phasor analysis | |
| CN108448568B (en) | Power distribution network hybrid state estimation method based on multiple time period measurement data | |
| CN101527455B (en) | Interconnected electric network distributed current calculating method on the basis of alternation and iteration of current module | |
| CN106921156A (en) | A kind of active distribution network method for estimating state based on many sampling period hybrid measurements | |
| CN103700036B (en) | A kind of transient stability projecting integral method being suitable to power system Multiple Time Scales | |
| CN107843810A (en) | A kind of active power distribution network fault section tuning on-line method based on state estimation | |
| CN104778367A (en) | On-line wide-area Thevenin equivalent parameter calculation method based on single state section | |
| CN104037759B (en) | Electric power system injects the robust least-square estimation method of transfer distribution factor | |
| CN109100614A (en) | A kind of transmission open acess system and method based on PMU device | |
| CN108448585A (en) | A kind of electric network swim equation solution method of linearization based on data-driven | |
| CN104156542A (en) | Implicit-projection-based method for simulating stability of active power distribution system | |
| CN114091816A (en) | Power distribution network state estimation method based on gated graph neural network of data fusion | |
| CN103956735B (en) | The harmonic trend of a kind of distributed generation system analyzes method | |
| CN108388688A (en) | A kind of active distribution network reliability estimation method | |
| CN110299762A (en) | Active distribution network Robust filter method based on PMU near-realtime data | |
| CN103532137A (en) | Method for estimating state of three-phase four-wire low-voltage distribution network | |
| CN106208099A (en) | A kind of Method for Reactive Power Optimization in Power based on bi-level programming and application thereof | |
| CN106372440B (en) | A kind of adaptive robust state estimation method of the power distribution network of parallel computation and device | |
| CN116316637A (en) | Dynamic topology identification method, system, equipment and storage medium for power distribution network | |
| CN105021952A (en) | A large power grid fault locating method based on wide-area phase angle measurement information | |
| Hu et al. | Frequency prediction model combining ISFR model and LSTM network |
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: 20190514 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |