CN104281978A - Available transmission capacity calculation method based on probabilistic power flow - Google Patents

Available transmission capacity calculation method based on probabilistic power flow Download PDF

Info

Publication number
CN104281978A
CN104281978A CN201410475723.7A CN201410475723A CN104281978A CN 104281978 A CN104281978 A CN 104281978A CN 201410475723 A CN201410475723 A CN 201410475723A CN 104281978 A CN104281978 A CN 104281978A
Authority
CN
China
Prior art keywords
electricity
atc
transmission
batc
transmission line
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201410475723.7A
Other languages
Chinese (zh)
Other versions
CN104281978B (en
Inventor
李卫东
刘奭昕
王海霞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dalian University of Technology
State Grid Liaoning Electric Power Co Ltd
Original Assignee
Dalian University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN201410475723.7A priority Critical patent/CN104281978B/en
Publication of CN104281978A publication Critical patent/CN104281978A/en
Application granted granted Critical
Publication of CN104281978B publication Critical patent/CN104281978B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • General Physics & Mathematics (AREA)
  • Development Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Theoretical Computer Science (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention relates to an available transmission capacity calculation method based on the probabilistic power flow. According to the method, the probabilistic ATC of a probabilistic power flow calculating system is used, and the calculated amount of the method is greatly reduced compared with that of a simulation method, so that the calculating speed of the system probabilistic ATC is increased, and the calculating speed of AQTC is increased. The concepts of transmission element ATC and system ATC influence factors of the transmission element ATC are raised, transmission elements are sequentially added to a power system according to the sequence of the influence factors from large to small, more influences of the transmission elements are considered in a system probabilistic ATC result, the deviation between the result and actual system probabilistic ATC expectation is gradually reduced, the requirement for the probabilistic ATC calculation accuracy of the system is met, and the requirement for the AQTC calculation accuracy is met. Thus, the available transmission capacity calculation method can be widely applied to the fields of probabilistic available electric quantity transmission capacity calculation and probabilistic available transmission capacity calculation.

Description

A kind of Calculation of Available Transfer Capability method based on Probabilistic Load Flow
Technical field
The present invention relates to a kind of available conveying electricity computing method, particularly about a kind of Calculation of Available Transfer Capability method based on Probabilistic Load Flow.
Background technology
Wind-powered electricity generation, as clean regenerative resource, in today that energy crisis is increasingly serious, obtains development energetically.But, along with the continuous expansion of wind-electricity integration scale, wind energy turbine set abandons the universal phenomenon that wind becomes operation of power networks just gradually.Because the interregional resource complementation of China is with the obvious advantage, by coordinating on a large scale to give full play to the complementary effect that interregional load avoids the peak hour between wrong paddy effect, power supply architecture, thus the regenerative resources such as wind-powered electricity generation of farthest dissolving, so wind-powered electricity generation needs to carry out long-distance sand transport across multi-region electric network.Generally speaking, the production schedule of electricity power enterprise is in the mode of annual electricity, is provided by higher authority; And wind power is very strong in short-term uncertainty, but there is certain rule on long terms, can more adequately be predicted.Therefore, wind-powered electricity generation transaction is more suitable for signing contract with long period yardstick.The quantity of wind-powered electricity generation transaction, not only relevant with received quantity with the wind energy quantity that both parties can provide, also depend on the transport capacity of approach electrical network.Therefore, how to determine the quantity of the available conveying wind-powered electricity generation of electrical network between both parties, namely determining the available power transport capacity (Available Quantity Transfer Capability, AQTC) between both parties, is one of key issue urgently to be resolved hurrily before both sides conclude the business.
So-called available power transport capacity (AQTC), refers to actual physics power transmission network at the appointed time residue and can be used for the electricity transport capacity that business uses in section.The concept of existing available transmission capacity (Available Transfer Capability, ATC) refers to remaining in actual physics power transmission network, to can be used for business use transmission capacity.Both compare, although the two is all available transport capacitys of reflection electrical network, what ATC reflected is power, and AQTC refers to electricity.
The influence factor of AQTC is a lot, the generation schedule in the such as cycle, load prediction, Unit Combination, turnaround plan etc., therefore its computing method are comparatively complicated, and is difficult to accurate Calculation, does not also have the calculating AQTC method of relative maturity at present.And use for reference existing ATC algorithm, a thinking can be provided for the calculating of AQTC: within the minor cycle that changes of operating modes is little, choose a typical way, the uncertainty of discontinuity surface when utilizing each in this cycle of probability simulation of ATC.Once larger change occurs the method for operation, just again choose typical way, calculate the probability ATC under new typical way.Finally, the probability ATC under each typical way and time scale are weighted, obtain macrocyclic available conveying electricity.
Choose typical way obtained good application in Practical Project, staff rule of thumb, finds out and comparatively ripe, easy chooses mode, such as summer maximum mode, summer minimum mode, winter maximum mode, winter minimum mode.But the calculating of probability ATC still also exists some shortcomings.At present, the computing method for probability available transmission capacity (probability ATC) mainly comprise following a few class:
(1) based on the algorithm of stochastic programming: in the process calculating ATC, first use SPR (two-stage Stochastic Programming with Recourse) method by discrete variable serialization, subsequently according to the result of calculation of SPR, use CCP (Chance Constrained Programming) method process continuous variable, try to achieve the ATC value under probability meaning.Because the method has used the process of discrete variable and continuous variable in computation process, so computing velocity is not very desirable.
(2) enumerative technique; The enumerative technique that later stage proposes is that enumerating of system state is combined calculating ATC with optimized algorithm.When calculating ATC, if the catastrophic failure of this system is less, then the method is relatively effective, but for the electric system of reality, the exponential time characteristic of enumerative technique still limits its application, and the medium-term and long-term ATC that really effectively cannot process Iarge-scale system calculates.
(3) based on the algorithm of Monte-Carlo Simulation: application Monte Carlo method is sampled to system state.The method can process the uncertain factor of huge number in actual electric network easily, and its computing time does not connect the increase of complexity with system scale or network and sharply increases, and is apply a kind of computing method very widely.But, in order to ensure precision, frequency in sampling during emulation is thousands of, and institute is very consuming time in this way, is difficult to ensure counting yield.
(4) based on the algorithm of Bootstrap: this method is applied to during long-term ATC calculates, the market information (the exerting oneself of such as generator, node load level etc.) that can make full use of nearest a period of time removes the sample as estimating.But, as a kind of emerging algorithm, the uncertainty (random fault etc. of such as power transmission line) of some network parameter can't be processed well, need further improvement.
(5) based on the ATC computing method of probability: document Lei Dong, Saifeng Li, Yihan Yang, Hai Bao.Asia-Pacific Power and Energy Engineering Conference (APPEEC2010), utilize Probabilistic Load Flow in 2010 " The Calculation of Transfer Reliability Margin Based on the Probabilistic Load Flow ", after taking into full account TRM and CBM two kinds of nargin, calculate to obtain ATC.Although speed significantly improves, only considered a bottleneck circuit, have ignored the probability impact of other transmission of electricity element, make the computational accuracy of this algorithm be difficult to be met.
In sum, in existing probability ATC computing method, method (1)-(4) model is too complicated, calculates too consuming time; (5) although computing velocity ensured, consider bottleneck condition very few, precision is inadequate.If probability ATC computing velocity is too slow, then the counting yield of AQTC can be made to be difficult to be protected, also just be difficult in engineering be used; If the precision that probability ATC calculates is inadequate, then can make the follow-up AQTC that calculates and actual value deviation excessive, do not have the directive function to both parties, cause waste or the shortage of electric power resource.
Summary of the invention
For the problems referred to above, the object of this invention is to provide a kind of Calculation of Available Transfer Capability method based on Probabilistic Load Flow taking into account accuracy and runtime.
For achieving the above object, the present invention takes following technical scheme: a kind of Calculation of Available Transfer Capability method based on Probabilistic Load Flow, it comprises the following steps: 1) to needing the electric system calculating AQTC, several typical way are chosen within the cycle, this typical way comprises typical way 1, typical way 2,, typical way n, and set the duration of each typical way within the cycle; 2) successively using the method for operation of each typical way as this electric system, calculate the probability ATC of this electric system under each method of operation to expect, it comprises the following steps: 2-1) read in the data of electric system under a certain method of operation that will calculate, this electric power system data comprises the parameters of each transmission of electricity element; And set this electric system accuracy requirement required under this method of operation; 2-2) respectively to transmit electricity according to gained the parameter of element, calculate the determinacy ATC of this electric system, be designated as ATC c; 2-3) make the ATC of gained cfor power is flowed out in feeder section injecting power, powered region, calculate the Probabilistic Load Flow distribution of now this electric system, it comprises the power distribution of transmission line of electricity and the voltage's distribiuting of transmission of electricity bus; 2-4) distribute according to the Probabilistic Load Flow of this electric system, calculate each transmission of electricity element to the factor of influence of system ATC, and from big to small transmission of electricity element is sorted according to factor of influence; 2-5) choose the initial transmission of electricity element of the sets of elements that the maximum transmission of electricity element of factor of influence is formed for interpolation transmission of electricity element as this electric system, calculate its transmission of electricity element ATC using the initial distribution curve of the probability distribution curve of the transmission of electricity element ATC of gained as this Probabilistic ATC, using the initial expectation of the expectation of initial distribution curve as this Probabilistic ATC; 2-6) according to factor of influence ranking results from big to small, add transmission of electricity element successively in sets of elements; 2-7) judge that whether the transmission of electricity element for adding is separate with each transmission of electricity element added; If independent, then carry out next step; If not independent, then reject, turn back to step 2-6) add next transmission of electricity element; 2-8) calculate institute and add the probability ATC of element of transmitting electricity, its distribution curve and this electric system original probability ATC distribution curve are carried out matching, obtains the probability ATC distribution curve that this electric system is new, and calculate the expectation of new probability ATC distribution curve; If 2-9) difference of the new expectation of this electric system gained and the former expectation of last gained meets the accuracy requirement preset, then using new probability ATC distribution curve as the final distribution curve of this electric system; If do not meet, then return step 2-6) continue to add new transmission of electricity element, until meet the accuracy requirement preset; 3) the probability ATC under each typical way is expected that correspondence was multiplied with the duration of each typical way, and the AQTC obtained in the cycle that each product of gained is sued for peace.
Described step 2) 2-4) comprise following content: transmission of electricity element factor of influence μ on system probability ATC affect as follows: establish transmission line of electricity i to the factor of influence μ of system probability ATC ifor: wherein, P imaxfor the peak power of the permission of transmission line of electricity i; P ifor the probability power of transmission line of electricity i; If bus j is to the factor of influence μ of system probability ATC jfor: wherein, V jminfor the lower limit of busbar voltage; V jmaxfor the upper limit of busbar voltage; V jfor the voltage of bus j; By the size of factor of influence μ value, can sort to the influence degree of each transmission of electricity element to electric system, it represents μ iand μ jbe referred to as factor of influence μ, μ value is larger, and influence degree is larger, and vice versa.
Described step 2) 2-5) comprise following content: transmission of electricity element ATC, comprise transmission line of electricity available transmission capacity LATC and transmission of electricity bus available transmission capacity BATC; 1. the available transmission capacity LATC of transmission line of electricity: set with any transmission line of electricity i as bottleneck factor, it is the available transmission capacity LATC of this transmission line of electricity i that electric system can be used further to the transmission capacity that business uses i; If Δ LATC ifor under the system state of the determination of known all parameters, during with transmission line of electricity i for bottleneck factor, the power that feed, powered two regions also can increase or reduce; If P i, Δ LATC iboth relation with increases are positive correlation, then have: LATC i = ATC c + &Delta; LATC i = ATC c + P i max - P i S Ai - S Bi ; Wherein, P imaxfor the peak power of the permission of transmission line of electricity i; P ifor the probability power of transmission line of electricity i; S aifor feeder section power increases the sensitivity to transmission line of electricity i power; S bifor powered area power increases the sensitivity to transmission line of electricity i power; If P i, Δ LATC iboth relation with increases are negative correlation, then think LATC now ifor infinite, do not consider that this transmission line of electricity i is on the impact of probability results; 2. transmit electricity the available transmission capacity BATC of bus: set with any one bus j as bottleneck factor, and electric system can be used further to the available transmission capacity BATC that transmission capacity that business uses is bus j j; If Δ BATC jfor under the system state of the determination of known all parameters, during with bus j for bottleneck factor, the power that feed, powered two regions also can increase or reduce; If V j, Δ BATC jboth relation with increases are positive correlation, then have: work as V j< V jmintime: ; Work as V j> V jmintime: BATC j = ATC c + &Delta; BATC j = ATC c + V j max - V j S Aj - S Bj ; Wherein, V jminfor the lower limit of busbar voltage; V jmaxfor the upper limit of busbar voltage; S ajfor feeder section power increases the voltage sensitivity to bus j; S bjfor powered area power increases the voltage sensitivity to bus j; S during positive correlation aj-S bj> 0; If V j, Δ BATC jboth relation with increases are negative correlation, then have: work as V j> V jmaxtime: ; Work as V j< V jmaxtime: BATC j = ATC c + &Delta; BATC j = ATC c + V j - V j min S Aj - S Bj .
Described step 2) 2-7) in the relevant condition of any two transmission lines of electricity as follows: if two transmission lines of electricity can meet following three requirements numerically simultaneously, then these two transmission lines of electricity are relevant: dominate degree of correlation β dbe greater than maximum relation degree S max75%: ; Total relevance β is greater than 75% of the leading degree of correlation: jack per line number λ is greater than 75% of total node number N or is less than 25% of total node number N: or
The present invention is owing to taking above technical scheme, it has the following advantages: 1, the present invention utilizes the probability ATC of probabilistic load flow system, compare simulation and greatly reduce calculated amount, thus improve the computing velocity of system probability ATC, and then improve the computing velocity of AQTC.2, transmission of electricity element ATC and the concept to system ATC factor of influence thereof are proposed, in electric system, transmission of electricity element is added successively according to factor of influence order from big to small, in the middle of the probability ATC result impact of more transmission of electricity elements being considered system, reduce the deviation expected with system actual probabilities ATC step by step, and then the probability ATC computational accuracy met required by system, meet the computational accuracy of AQTC simultaneously.In view of above reason, the present invention can be widely used in probability available power transport capacity and calculate and probability Calculation of Available Transfer Capability field.
Accompanying drawing explanation
Fig. 1 is that typical way chooses schematic diagram
Fig. 2 is the schematic flow sheet that in the present invention, probability ATC calculates
Fig. 3 is the probability distribution schematic diagram of transmission line of electricity i power
Fig. 4 is the power P of transmission line of electricity i iwith Δ LATC ibetween qualitative relationships schematic diagram
Fig. 5 is the probability distribution schematic diagram of transmission of electricity bus j voltage
Fig. 6 is the voltage V of transmission of electricity bus j jwith Δ BATC jbetween qualitative relationships schematic diagram
Embodiment
Below in conjunction with drawings and Examples, the present invention is described in detail.
1) as shown in Figure 1, first needs are calculated to the electric system of AQTC (Available Quantity Transfer Capability, available power transport capacity), several typical way are chosen within the cycle, i.e. typical way 1, typical way 2 ..., typical way n.Set the duration of each typical way within the cycle, selection principle is as follows: within the minor cycle that changes of operating modes is little, choose a typical way representatively, once larger change occurs the method for operation, just again chooses typical way.In practical power systems, staff has found out comparatively ripe typical way choosing method according to the experience run for many years, can directly apply into.
2) successively using the method for operation of each typical way as this electric system, and calculate the expectation of the probability ATC of this electric system under each method of operation, it comprises the following steps:
It should be noted that, the step that probability ATC expectation is asked in this electric system under different running method is consistent, and therefore no longer repeat specification is only representatively described with the process of asking for of this kind of method of operation lower probability ATC.
2-1) as shown in Figure 2, read in the data of electric system under this method of operation that will calculate, this electric power system data comprises the parameters of the transmission of electricity such as generator, transmission line of electricity, transformer, the reactive power compensator element commonly used this area, to calculate determinacy ATC and the Probabilistic Load Flow of this electric system according to these existing parameters.Utilize the result of Probabilistic Load Flow can obtain the probability ATC of transmission of electricity element, be the available transmission capacity of transmission of electricity element, it comprises the available transmission capacity LATC (Line Available Transfer Capability) of transmission line of electricity and the available transmission capacity BATC (Bus Available Transfer Capability) of transmission of electricity bus.The transmission of electricity element that refering in particular in the present invention affects Probabilistic ATC size is transmission line of electricity and generator, and the ATC of transmission line of electricity is LATC, and the ATC of generator is BATC.
2-2) respectively transmit electricity according to gained the parameter of element, the ATC algorithm based on continuous tide utilizing this area to commonly use calculates the determinacy ATC of electric system, is designated as ATC c;
2-3) make the ATC of gained cfor power is flowed out in feeder section injecting power, powered region, adopt the Probabilistic Load Flow algorithm in conjunction with cumulant and Gram-Charlier series expansion that this area is commonly used, calculate the Probabilistic Load Flow distribution of now this electric system, it comprises the power distribution of transmission line of electricity and the voltage's distribiuting of transmission of electricity bus;
2-4) distribute according to the Probabilistic Load Flow of this electric system, calculate each transmission of electricity element to the factor of influence of system ATC, and from big to small transmission of electricity element is sorted according to factor of influence;
System probability ATC is limited by the transmission of electricity element of out-of-limit (having exceeded operational limit) at first, so operate in from the transmission of electricity element close to ultimate limit state, may be larger on the impact of system probability ATC, therefore propose transmission of electricity element factor of influence μ in order to respectively transmit electricity element to the size of system probability ATC influence degree.
If transmission line of electricity i is to the factor of influence μ of system probability ATC ifor:
&mu; i = 1 - &Integral; 0 P i max P i ( x ) dx - - - ( 1 )
If bus j is to the factor of influence μ of system probability ATC jfor:
&mu; j = 1 - &Integral; V i min V i max V j ( x ) dx - - - ( 2 )
By factor of influence μ (μ iand μ jbe referred to as factor of influence μ) size of value, can sort to the influence degree of each transmission of electricity element to electric system: μ value is larger, and influence degree is larger, and vice versa.So, when interpolation transmission of electricity element carries out probability simulation iteration, the less transmission of electricity element of influence degree can just be ignored when meeting precision.For the real system that, bottleneck transmission of electricity element comparatively large for system scale is comparatively concentrated, greatly reduce calculated amount.
2-5) choose the initial transmission of electricity element of the sets of elements that the maximum transmission of electricity element of factor of influence is formed for interpolation transmission of electricity element as this electric system, calculate its transmission of electricity element ATC using the initial distribution curve of its probability distribution curve as this Probabilistic ATC, using the initial expectation of the expectation of initial distribution curve as this Probabilistic ATC;
1. the available transmission capacity LATC of transmission line of electricity
If with any transmission line of electricity i for bottleneck factor, it is the available transmission capacity LATC of this transmission line of electricity i that electric system can be used further to the transmission capacity that business uses i.
As shown in Figure 3, Figure 4, if Δ LATC ithe power that under the system state of the determination of known all parameters, during with transmission line of electricity i for bottleneck factor, feed, powered two regions also can increase (or minimizing), wherein, Δ LATC i>0 illustrates under this system state, and transmission line of electricity i does not also reach itself operational limit; Δ LATC i=0 illustrates under this system state, and transmission line of electricity i just reaches operational limit; Δ LATC i<0 illustrates under this system state, and transmission line of electricity i has exceeded operational limit already.If P ifor the probability power of transmission line of electricity i.
If along with Δ LATC iincrease, power (absolute value) P of transmission line of electricity i ialso increase, i.e. Δ LATC iwith P irelation with increase be positive correlation, then both sides relation is as follows:
&Delta; LATC i = P i max - P i S Ai - S Bi S Ai - S Bi > 0 - - - ( 3 )
Wherein, P imaxfor the maximum transmission power that transmission line of electricity i allows; S aifor feeder section power increases the sensitivity to transmission line of electricity i power; S bifor powered area power increases the sensitivity to transmission line of electricity i power, owing to being load growth, so band negative sign; S ai-S bi> 0 is equivalent to Δ LATC iwith P ibetween relation with increase be positive correlation.
If along with Δ LATC iincrease, the power (absolute value) of transmission line of electricity i reduces, i.e. Δ LATC iwith P irelation with increase be negative correlation (S ai-S bi< 0), so can think Δ LATC now ifor just infinite.Under normal circumstances, if along with Δ LATC iincrease, the power (absolute value) of transmission line of electricity i reduces, and is so constantly increasing Δ LATC iprocess in, transmission line of electricity i does not reach oneself the limit forever, because the limit of circuit is an inaccessiable maximal value, does not have lower limit, so can think Δ LATC ifor just infinite, in other words when negative correlation, do not consider that this transmission line of electricity i is on the impact of probability results.
Formula (4) is utilized to calculate the available transmission capacity LATC of transmission line of electricity i i:
LATC i = ATC c + &Delta; LATC i = ATC c + P i max - P i S Ai - S Bi - - - ( 4 )
2. transmit electricity the available transmission capacity BATC of bus
If with any one bus j for bottleneck factor, electric system can be used further to transmission capacity that business the uses available transmission capacity BATC that is bus j j.
As shown in Figure 5, Figure 6, if Δ BATC jthe power that under the system state of the determination of known all parameters, during with bus j for bottleneck factor, feed, powered two regions also can increase (or minimizing), wherein, Δ BATC j>0 illustrates under this system state, and bus j does not also reach operational limit; Δ BATC j=0 illustrates under this system state, and bus j just reaches operational limit; Δ BATC j<0 illustrates under this system state, and bus j has exceeded operational limit already.
If V j, Δ BATC jboth relation with increases are positive correlation, then have:
Work as V j< V jmintime: BATC j=0 (5)
Work as V j> V jmintime: BATC j = ATC c + V j max - V j S Aj - S Bj - - - ( 6 )
Wherein, V jfor the voltage of bus j; V jminfor the lower limit of the voltage of bus j; V jmaxfor the upper limit of the voltage of bus j; S ajfor feeder section power increases the voltage sensitivity to bus j; S bjfor powered area power increases the voltage sensitivity to bus j; S during positive correlation aj-S bj> 0.
If V j, Δ BATC jboth relation with increases are negative correlation, then have:
Work as V j> V jmaxtime: BATC j=0 (7)
Work as V j< V jmaxtime: BATC j = ATC c + V j - V j min S Aj - S Bj - - - ( 8 )
It should be noted that, due to the power P of transmission line of electricity i i, transmission of electricity bus j voltage V jfor probable value, so the available transmission capacity LATC of transmission line of electricity i iwith the available transmission capacity BATC of transmission of electricity bus j jbe probabilistic ATC.
2-6) according to factor of influence ranking results from big to small, add transmission of electricity element successively in sets of elements, owing to multiple transmission of electricity element can be considered in the middle of probability results successively according to factor of influence size, so can think, add transmission of electricity element more, from exact value more close to, precision is higher, while meeting computing velocity, accomplish the control in precision;
2-7) judge that whether the transmission of electricity element for adding is separate with each transmission of electricity element added; If independent, then carry out next step; If relevant (not independent), then reject, turn back to step 2-6) add next transmission of electricity element, judge that independently process is as follows:
If be bus for the transmission of electricity element added, because be independently between bus and transmission line of electricity, then directly add and carry out next step, not needing to judge the mutual independence with the transmission of electricity element added; If the transmission of electricity element added is transmission line of electricity, then carries out the calculating of the degree of correlation with the transmission line of electricity that added in sets of elements, judge mutual independence with this, do not need to carry out independence judgement with the transmission of electricity bus that added in sets of elements.
Correlativity determination methods between transmission line of electricity is as follows:
If any two transmission lines of electricity are m-n, p-q, define the total relevance (Line Relevancy, LR) of these two transmission lines of electricity for they are to the sum of products of each node injection sensitivity, then the total relevance β between any two transmission lines of electricity:
&beta; = &Sigma; j = 1 N S m - n , j &times; S p - q , j - - - ( 9 )
Wherein, β is the total relevance between transmission line of electricity m-n, p-q; N is calculated electric system interior joint number; S m-n, ipower for calculated electric system median generatrix j injects the sensitivity to transmission line of electricity m-n; S p-q, ipower for calculated electric system median generatrix j injects the sensitivity to transmission line of electricity p-q.
In view of generally, transmission line of electricity power has the greatest impact by transmission line of electricity two-end-point injecting power, be embodied in the middle of the l-w sensitivity matrix between transmission line of electricity power l and node injecting power w, show as end points injecting power to the sensitivity number of transmission line of electricity at all node powers to maximum among its sensitivity.
So the injecting power sensitivity of transmission line of electricity end points is considered separately, the sensitivity sum of products of definition transmission line of electricity end points to two transmission lines of electricity takes the degree of correlation (Dominant Relevancy as the leading factor, DR), if two transmission lines of electricity have same endpoints, then coincidence node is only calculated once:
&beta; D = &Sigma; j = m , n , p , q S m - n , j &times; S p - q , j - - - ( 10 )
Wherein, β dfor the leading degree of correlation of transmission line of electricity m-n and transmission line of electricity p-q.
Defining the sensitivity sum of products of all the other nodes to two transmission lines of electricity is the non-dominant degree of correlation (Non-Dominant Relevancy, NDR):
&beta; N = &Sigma; j &NotEqual; m , n , p , q S m - n , j &times; S p - q , j - - - ( 11 )
Wherein, β nfor the non-dominant degree of correlation of transmission line of electricity m-n and transmission line of electricity p-q.
Add up the jack per line number λ of two transmission lines of electricity to each node sensitivity, and calculate the maximal value of sensitivity product, i.e. maximum relation degree S max:
S max＝max{S m-n,j×S p-q,j} (12)
The leading degree of correlation embodies each Pilot bus (transmission line of electricity end points) random fluctuation affects sum to two transmission line of electricity power.If two transmission lines of electricity are correlated with, then the product of each Pilot bus sensitivity keeps jack per line (or contrary sign nodal values is very little) substantially, and the leading degree of correlation is larger.
Total relevance embodies each node random fluctuation to two transmission line of electricity power influences sums.If two transmission lines of electricity are correlated with, then the Pilot bus degree of correlation can be comparatively large, the non-dominant degree of correlation should with leading degree of correlation jack per line, or with its contrary sign but numerical value is very little, the total relevance that both add up can be larger.
Jack per line number embodies each node random fluctuation affects direction to two transmission line of electricity power.If two transmission line of electricity positive correlations, then jack per line nodes will much larger than contrary sign nodes; If two transmission line of electricity negative correlation, then jack per line nodes will much smaller than contrary sign nodes.
Three requirements numerically below if two transmission lines of electricity can meet simultaneously, then these two transmission lines of electricity are relative to each other, namely independent: leading degree of correlation β dbe greater than maximum relation degree S max75%: total relevance β is greater than leading degree of correlation β d75%: jack per line number λ is greater than 75% of total node number N or is less than 25% of total node number N: or
2-8) calculate institute and add the probability ATC of element of transmitting electricity, its distribution curve and this electric system original probability ATC distribution curve are carried out matching, obtain the probability ATC distribution curve that this electric system is new, and calculate the expectation of new probability ATC distribution curve, fit procedure is as follows:
P AB(x)＝P A(x)×P B(B≥x)+P A(A≥x)×P B(x)-P A(x)×P B(x) (13)
Wherein, AB is the system be made up of any two transmission of electricity elements A, B; P aBfor the distribution of the probability ATC of AB system; P afor the probability ATC of transmission of electricity element A distributes; P bfor the probability ATC of transmission of electricity element B distributes; X ∈ (0 ,+∞).
Introduce cumulative function F A ( x ) = &Integral; - &infin; x P A ( x ) dx , F B ( x ) = &Integral; - &infin; x P B ( x ) dx , Then have:
P AB(x)＝P A(x)×(1-F B(x))+(1-F A(x))×P B(x)+P A(x)×P B(x) (14)
If 2-9) difference of the new expectation of this electric system gained and the former expectation of last gained meets the accuracy requirement preset, then using new probability ATC distribution curve as the final distribution curve of this electric system; If do not meet, then return step 2-6) continue to add new transmission of electricity element, until meet the accuracy requirement preset.
3) the probability ATC under each typical way is expected that correspondence was multiplied with the duration of each typical way, and the AQTC obtained in the cycle that each product of gained is sued for peace.
The above; be only the present invention's preferably embodiment; but protection scope of the present invention is not limited thereto; anyly be familiar with those skilled in the art in the technical scope that the present invention discloses; be equal to according to technical scheme of the present invention and inventive concept thereof and replace or change, all should be encompassed within protection scope of the present invention.

Claims (5)

1., based on a Calculation of Available Transfer Capability method for Probabilistic Load Flow, it comprises the following steps:
1) to needing the electric system calculating AQTC, within the cycle, choose several typical way, this typical way comprises typical way 1, typical way 2 ..., typical way n, and set the duration of each typical way within the cycle;
2) successively using the method for operation of each typical way as this electric system, calculate the probability ATC of this electric system under each method of operation and expect, it comprises the following steps:
2-1) read in the data of electric system under a certain method of operation that will calculate, this electric power system data comprises the parameters of each transmission of electricity element; And set this electric system accuracy requirement required under this method of operation;
2-2) respectively to transmit electricity according to gained the parameter of element, calculate the determinacy ATC of this electric system, be designated as ATC c;
2-3) make the ATC of gained cfor power is flowed out in feeder section injecting power, powered region, calculate the Probabilistic Load Flow distribution of now this electric system, it comprises the power distribution of transmission line of electricity and the voltage's distribiuting of transmission of electricity bus;
2-4) distribute according to the Probabilistic Load Flow of this electric system, calculate each transmission of electricity element to the factor of influence of system ATC, and from big to small transmission of electricity element is sorted according to factor of influence;
2-5) choose the initial transmission of electricity element of the sets of elements that the maximum transmission of electricity element of factor of influence is formed for interpolation transmission of electricity element as this electric system, calculate its transmission of electricity element ATC using the initial distribution curve of the probability distribution curve of the transmission of electricity element ATC of gained as this Probabilistic ATC, using the initial expectation of the expectation of initial distribution curve as this Probabilistic ATC;
2-6) according to factor of influence ranking results from big to small, add transmission of electricity element successively in sets of elements;
2-7) judge that whether the transmission of electricity element for adding is separate with each transmission of electricity element added; If independent, then carry out next step; If not independent, then reject, turn back to step 2-6) add next transmission of electricity element;
2-8) calculate institute and add the probability ATC of element of transmitting electricity, its distribution curve and this electric system original probability ATC distribution curve are carried out matching, obtains the probability ATC distribution curve that this electric system is new, and calculate the expectation of new probability ATC distribution curve;
If 2-9) difference of the new expectation of this electric system gained and the former expectation of last gained meets the accuracy requirement preset, then using new probability ATC distribution curve as the final distribution curve of this electric system; If do not meet, then return step 2-6) continue to add new transmission of electricity element, until meet the accuracy requirement preset;
3) the probability ATC under each typical way is expected that correspondence was multiplied with the duration of each typical way, and the AQTC obtained in the cycle that each product of gained is sued for peace.
2. a kind of Calculation of Available Transfer Capability method based on Probabilistic Load Flow as claimed in claim 1, is characterized in that: described step 2) 2-4) comprise following content:
The factor of influence μ of transmission of electricity element affects as follows on system probability ATC:
If transmission line of electricity i is to the factor of influence μ of system probability ATC ifor:
&mu; i = 1 - &Integral; 0 P i max P i ( x ) dx
Wherein, P imaxfor the peak power of the permission of transmission line of electricity i; P ifor the probability power of transmission line of electricity i;
If bus is to the factor of influence μ of system probability ATC jfor:
&mu; j = 1 - &Integral; V i min V i max V j ( x ) dx
Wherein, V jminfor the lower limit of busbar voltage; V jmaxfor the upper limit of busbar voltage; V jfor the voltage of bus j;
By the size of factor of influence μ value, can sort to the influence degree of each transmission of electricity element to electric system, it represents μ iand μ jbe referred to as factor of influence μ, μ value is larger, and influence degree is larger, and vice versa.
3. a kind of Calculation of Available Transfer Capability method based on Probabilistic Load Flow as claimed in claim 1, is characterized in that: described step 2) 2-5) comprise following content:
The ATC of transmission of electricity element, comprises the available transmission capacity LATC of transmission line of electricity and the available transmission capacity BATC of transmission of electricity bus;
1. the available transmission capacity LATC of transmission line of electricity
If with any transmission line of electricity i for bottleneck factor, it is the available transmission capacity LATC of this transmission line of electricity i that electric system can be used further to the transmission capacity that business uses i; If Δ LATC ifor under the system state of the determination of known all parameters, during with transmission line of electricity i for bottleneck factor, the power that feed, powered two regions also can increase or reduce;
If P i, Δ LATC iboth relation with increases are positive correlation, then have:
LATC i = ATC c + &Delta; LATC i = ATC c + P i max - P i S Ai - S Bi
Wherein, P imaxfor the peak power of the permission of transmission line of electricity i; P ifor the probability power of transmission line of electricity i; S aifor feeder section power increases the sensitivity to transmission line of electricity i power; S bifor powered area power increases the sensitivity to transmission line of electricity i power;
If P i, Δ LATC iboth relation with increases are negative correlation, then think LATC now ifor infinite, do not consider that this transmission line of electricity i is on the impact of probability results;
2. transmit electricity the available transmission capacity BATC of bus
If with any one bus j for bottleneck factor, electric system can be used further to transmission capacity that business the uses available transmission capacity BATC that is bus j j; If Δ BATC jfor under the system state of the determination of known all parameters, during with bus j for bottleneck factor, the power that feed, powered two regions also can increase or reduce;
If V j, Δ BATC jboth relation with increases are positive correlation, then have:
Work as V j< V jmintime: BATC j=0
Work as V j> V jmintime: BATC j = ATC c + &Delta; BATC j = ATC c + V j max - V j S Aj - S Bj
Wherein, V jminfor the lower limit of busbar voltage; V jmaxfor the upper limit of busbar voltage; S ajfor feeder section power increases the voltage sensitivity to bus j; S bjfor powered area power increases the voltage sensitivity to bus j; S during positive correlation aj-S bj> 0;
If V j, Δ BATC jboth relation with increases are negative correlation, then have:
Work as V j> V jmaxtime: BATC j=0
Work as V j< V jmaxtime: BATC j = ATC c + &Delta; BATC j = ATC c + V j - V j min S Aj - S Bj .
4. a kind of Calculation of Available Transfer Capability method based on Probabilistic Load Flow as claimed in claim 2, is characterized in that: described step 2) 2-5) comprise following content:
The ATC of transmission of electricity element, comprises the available transmission capacity LATC of transmission line of electricity and the available transmission capacity BATC of transmission of electricity bus;
1. the available transmission capacity LATC of transmission line of electricity
If with any transmission line of electricity i for bottleneck factor, it is the available transmission capacity LATC of this transmission line of electricity i that electric system can be used further to the transmission capacity that business uses i; If Δ LATC ifor under the system state of the determination of known all parameters, during with transmission line of electricity i for bottleneck factor, the power that feed, powered two regions also can increase or reduce;
If P i, Δ LATC iboth relation with increases are positive correlation, then have:
LATC i = ATC c + &Delta; LATC i = ATC c + P i max - P i S Ai - S Bi
Wherein, P imaxfor the peak power of the permission of transmission line of electricity i; P ifor the probability power of transmission line of electricity i; S aifor feeder section power increases the sensitivity to transmission line of electricity i power; S bifor powered area power increases the sensitivity to transmission line of electricity i power;
If P i, Δ LATC iboth relation with increases are negative correlation, then think LATC now ifor infinite, do not consider that this transmission line of electricity i is on the impact of probability results;
2. transmit electricity the available transmission capacity BATC of bus
If with any one bus j for bottleneck factor, electric system can be used further to transmission capacity that business the uses available transmission capacity BATC that is bus j j; If Δ BATC jfor under the system state of the determination of known all parameters, during with bus j for bottleneck factor, the power that feed, powered two regions also can increase or reduce;
If V j, Δ BATC jboth relation with increases are positive correlation, then have:
Work as V j< V jmintime: BATC j=0
Work as V j> V jmintime: BATC j = ATC c + &Delta; BATC j = ATC c + V j max - V j S Aj - S Bj
Wherein, V jminfor the lower limit of busbar voltage; V jmaxfor the upper limit of busbar voltage; S ajfor feeder section power increases the voltage sensitivity to bus j; S bjfor powered area power increases the voltage sensitivity to bus j; S during positive correlation aj-S bj> 0;
If V j, Δ BATC jboth relation with increases are negative correlation, then have:
Work as V j> V jmaxtime: BATC j=0
Work as V j< V jmaxtime: BATC j = ATC c + &Delta; BATC j = ATC c + V j - V j min S Aj - S Bj .
5. a kind of Calculation of Available Transfer Capability method based on Probabilistic Load Flow as claimed in claim 1 or 2 or 3 or 4, is characterized in that: described step 2) 2-7) in the relevant condition of any two transmission lines of electricity as follows:
Three requirements numerically below if two transmission lines of electricity can meet simultaneously, then these two transmission lines of electricity are relevant: leading degree of correlation β dbe greater than maximum relation degree S max75%: total relevance β is greater than 75% of the leading degree of correlation: jack per line number λ is greater than 75% of total node number N or is less than 25% of total node number N: &lambda; > 3 4 N Or &lambda; < 1 4 N .
CN201410475723.7A 2014-09-12 2014-09-17 A kind of Calculation of Available Transfer Capability method based on Probabilistic Load Flow Expired - Fee Related CN104281978B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410475723.7A CN104281978B (en) 2014-09-12 2014-09-17 A kind of Calculation of Available Transfer Capability method based on Probabilistic Load Flow

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
CN201410467880 2014-09-12
CN2014104678803 2014-09-12
CN201410467880.3 2014-09-12
CN201410475723.7A CN104281978B (en) 2014-09-12 2014-09-17 A kind of Calculation of Available Transfer Capability method based on Probabilistic Load Flow

Publications (2)

Publication Number Publication Date
CN104281978A true CN104281978A (en) 2015-01-14
CN104281978B CN104281978B (en) 2017-07-28

Family

ID=52256829

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410475723.7A Expired - Fee Related CN104281978B (en) 2014-09-12 2014-09-17 A kind of Calculation of Available Transfer Capability method based on Probabilistic Load Flow

Country Status (1)

Country Link
CN (1) CN104281978B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104617578A (en) * 2015-02-26 2015-05-13 河海大学 Method for acquiring available power transmission capability of power system with wind power plant
CN104734155A (en) * 2015-04-14 2015-06-24 河海大学 Method of acquiring available transfer capacity of electricity-gas interconnected energy system
CN107370157A (en) * 2017-09-12 2017-11-21 四川大学 A kind of available transfer capacity of transmission network risk benefit decision method based on trend entropy
CN112736914A (en) * 2020-12-29 2021-04-30 国网吉林省电力有限公司 Available transmission capacity probability calculation method considering wind power correlation
CN113363981A (en) * 2021-07-12 2021-09-07 中国电力科学研究院有限公司 Available transmission capacity calculation method, system, device and storage medium

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101895116A (en) * 2010-07-16 2010-11-24 国电南瑞科技股份有限公司 Online available transmission capacity computing method based on distribution factor
CN103985058A (en) * 2014-04-04 2014-08-13 河海大学 Available transfer capability calculation method based on improved multiple centrality-correction interior point method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101895116A (en) * 2010-07-16 2010-11-24 国电南瑞科技股份有限公司 Online available transmission capacity computing method based on distribution factor
CN103985058A (en) * 2014-04-04 2014-08-13 河海大学 Available transfer capability calculation method based on improved multiple centrality-correction interior point method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李赛峰: "基于概率模型的可用输电能力计算", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104617578A (en) * 2015-02-26 2015-05-13 河海大学 Method for acquiring available power transmission capability of power system with wind power plant
CN104617578B (en) * 2015-02-26 2017-05-03 河海大学 Method for acquiring available power transmission capability of power system with wind power plant
CN104734155A (en) * 2015-04-14 2015-06-24 河海大学 Method of acquiring available transfer capacity of electricity-gas interconnected energy system
CN107370157A (en) * 2017-09-12 2017-11-21 四川大学 A kind of available transfer capacity of transmission network risk benefit decision method based on trend entropy
CN112736914A (en) * 2020-12-29 2021-04-30 国网吉林省电力有限公司 Available transmission capacity probability calculation method considering wind power correlation
CN112736914B (en) * 2020-12-29 2022-11-11 国网吉林省电力有限公司 Available transmission capacity probability calculation method considering wind power correlation
CN113363981A (en) * 2021-07-12 2021-09-07 中国电力科学研究院有限公司 Available transmission capacity calculation method, system, device and storage medium

Also Published As

Publication number Publication date
CN104281978B (en) 2017-07-28

Similar Documents

Publication Publication Date Title
CN104281978B (en) A kind of Calculation of Available Transfer Capability method based on Probabilistic Load Flow
Xu et al. Trajectory sensitivity analysis on the equivalent one‐machine‐infinite‐bus of multi‐machine systems for preventive transient stability control
CN102928738A (en) Dynamic layer diagnostic device and method of smart power grid fault
CN104779611A (en) Economic dispatch method for micro grid based on centralized and distributed double-layer optimization strategy
Chen et al. Review of restoration technology for renewable‐dominated electric power systems
CN103825272A (en) Reliability determination method for power distribution network with distributed wind power based on analytical method
CN104318397A (en) Risk assessment and analysis method based on power grid short-time run-time behaviors
CN115456304A (en) Offshore wind farm reliability index calculation method and device considering typhoon influence
CN103971175A (en) Short-term load prediction method of multistage substations
CN107171365A (en) Multiple target stochastic and dynamic economic load dispatching method with asynchronous iteration is decoupled based on scene
CN102570451A (en) Static reliability assessment method for power transmission network
Shabbir et al. Short-term wind energy forecasting using deep learning-based predictive analytics
CN104021315A (en) Method for calculating station service power consumption rate of power station on basis of BP neutral network
Sharma et al. Particle swarm optimization approach for economic load dispatch: a review
CN105787604A (en) Transmission-distribution-cooperative distribution network situation quick awareness method
Wang et al. Reliability evaluation method for distribution network with distributed generations considering feeder fault recovery and network reconfiguration
CN117473384A (en) Power grid line safety constraint identification method, device, equipment and storage medium
Shan et al. A Day-Ahead Economic Dispatch Scheme for Transmission System with High Penetration of Renewable Energy
Wang et al. Optimal recovery strategy of DERs integrated distribution network based on scheduling rationality
CN110097292A (en) A kind of Demand-side interconnection effect reliability estimation method based on Monte Carlo
CN116613762A (en) Method and system for identifying key nodes in power grid area
Ofosu et al. Electricity consumption forecast for Tarkwa using autoregressive integrated moving average and adaptive neuro fuzzy inference system
Bienstock et al. Robust linear control of nonconvex battery operation in transmission systems
CN113393074A (en) Model construction method for power generation extension planning considering flexibility of external network
CN105610160A (en) Power supply security assessment method and system for power transmission 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
CB03 Change of inventor or designer information

Inventor after: Liu Shixin

Inventor after: Jin Shijun

Inventor after: Zhao Shiyi

Inventor after: Li Zhengwen

Inventor after: Wang Haixia

Inventor after: Li Weidong

Inventor before: Li Weidong

Inventor before: Liu Shixin

Inventor before: Wang Haixia

CB03 Change of inventor or designer information
TA01 Transfer of patent application right

Effective date of registration: 20170615

Address after: 116024 Dalian high tech park, Liaoning Ling Road, No. 2

Applicant after: DALIAN University OF TECHNOLOGY

Applicant after: STATE GRID LIAONING ELECTRIC POWER SUPPLY Co.,Ltd.

Address before: 116024 Dalian high tech park, Liaoning Ling Road, No. 2

Applicant before: Dalian University of Technology

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170728

Termination date: 20210917

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