CN105914755A - Electric power system dynamic reactive power optimization configuration method and system - Google Patents

Abstract

The invention relates to an electric power system dynamic reactive power optimization configuration method and system. Time domain simulation analysis is performed on each node based on power grid basic data. The transient voltage stability evaluation index of an electric power system is calculated according to the simulation result so that the voltage qualified rate of the electric power system can be obtained. When the voltage qualified rate is less than a preset qualified threshold, candidate nodes are extracted and reactive power compensation sensitivity is calculated. Dynamic reactive power compensation capacity of the corresponding candidate nodes is increased according to the reactive power compensation sensitivity until the voltage qualified rate reaches the preset qualified threshold, and finally an optimal reactive power compensation capacity scheme of the electric power system is obtained and outputted. The dynamic characteristics of the system under the condition of different faults are analyzed through the sensitivity, and dynamic reactive power optimization configuration is realized according to the reasonable criterion and index so that the dynamic reactive power installing position and the capacity size of the electric power system can be automatically calculated, voltage coordination control under the situation of complex reactive power configuration of the power grid can be systematically guided, and voltage stability of the electric power system can be enhanced.

Description

Electrical power system dynamic reactive power Optimal Configuration Method and system

Technical field

The present invention relates to electric power network technical field, particularly relate to a kind of electrical power system dynamic reactive power optimization and join Put method and system.

Background technology

Power distribution network is the end link of electric energy configuration, has that overlay area is wide, a feature such as circuit and load intensive, Its quality of power supply, safety in operation are not only directly connected to the economic benefit of users and electric power enterprise, also It it is the important embodiment of social benefit.

The reactive power source distribution of power system directly affects the safety of power system and stablizes.On the one hand, nothing Merit underpower will cause system voltage to reduce, electrical equipment can not reasonable employment, voltage even can be caused to collapse Burst etc. a series of accidents；On the other hand, reactive power surplus can cause unnecessary investment waste, even leads Cause system voltage deterioration, harm system and the safety of equipment.Large-scale receiving end urban distribution network load is heavier, The load proportions such as motor are higher, and when suffering the large disturbances such as short circuit, Enhancement of Transient Voltage Stability problem is particularly Prominent.

Therefore, how distribution network is carried out idle work optimization and improve system voltage stability, to power system Safety and economic operation is significant.

Summary of the invention

Based on this, it is necessary to for the problems referred to above, it is provided that a kind of electricity improving Power System Voltage Stability Force system Dynamic reactive power optimization collocation method and system.

A kind of electrical power system dynamic reactive power Optimal Configuration Method, comprises the following steps:

Receiving the electrical network basic data of power system, described electrical network basic data includes type and the joint of fault The parameter of point；

Type according to described fault and the parameter of described node, carry out time-domain-simulation analysis to each node, Node voltage data to corresponding node；

The Enhancement of Transient Voltage Stability assessment being calculated described power system according to described node voltage data refers to Mark, and the rate of qualified voltage of described power system it is calculated according to described Enhancement of Transient Voltage Stability evaluation index；

In described rate of qualified voltage less than when presetting qualified threshold value, obtain voltage and do not meet the node of preset standard Obtain both candidate nodes, and calculate the reactive-load compensation sensitivity of each both candidate nodes；

The dynamic passive compensation capacity of corresponding both candidate nodes is increased according to described reactive-load compensation sensitivity, until institute Stating rate of qualified voltage is described default qualified threshold value；

Dynamic passive compensation capacity according to each both candidate nodes obtains the optimal reactive compensation of described power system to be held Amount scheme also exports.

A kind of electrical power system dynamic reactive power Optimizing Configuration System, including:

Data acquisition module, for receiving the electrical network basic data of power system, described electrical network basic data bag Include the type of fault and the parameter of node；

Time-domain-simulation module, for the type according to described fault and the parameter of described node, enters each node Row time-domain-simulation is analyzed, and obtains the node voltage data of corresponding node；

Qualification rate computing module, for being calculated the temporary of described power system according to described node voltage data State voltage stability evaluation index, and it is calculated described electricity according to described Enhancement of Transient Voltage Stability evaluation index The rate of qualified voltage of Force system；

Node screening module, for when described rate of qualified voltage is less than default qualified threshold value, obtaining voltage not The node meeting preset standard obtains both candidate nodes, and calculates the reactive-load compensation sensitivity of each both candidate nodes；

Reactive compensation module, for increasing the dynamic nothing of corresponding both candidate nodes according to described reactive-load compensation sensitivity Merit compensation capacity, until described rate of qualified voltage is described default qualified threshold value；

Optimize output module, for obtaining described power train according to the dynamic passive compensation capacity of each both candidate nodes System optimal reactive compensation volume solutions and export.

Above-mentioned electrical power system dynamic reactive power Optimal Configuration Method and system, based on electrical network basic data to each node Carrying out time-domain-simulation analysis, the Enhancement of Transient Voltage Stability evaluation index calculating power system according to simulation result enters And obtain the rate of qualified voltage of power system.In rate of qualified voltage less than extracting candidate's joint when presetting qualified threshold value Point calculates reactive-load compensation sensitivity, and the dynamic reactive increasing corresponding both candidate nodes according to reactive-load compensation sensitivity is mended Repay capacity, until rate of qualified voltage is for presetting qualified threshold value, finally obtain the optimal reactive compensation of power system Volume solutions also exports.By the dynamic characteristic of system in the case of sensitive analysis different faults, according to rationally Criterion and index realize distributing rationally of dynamic reactive, can automatically calculate the installation of electrical power system dynamic reactive power Position and amount of capacity, it is possible to systematically instruct the voltage under the complicated idle configuration situation of electrical network to coordinate to control, Improve Power System Voltage Stability.

Accompanying drawing explanation

Fig. 1 is the flow chart of electrical power system dynamic reactive power Optimal Configuration Method in an embodiment；

Fig. 2 is the structure chart of electrical power system dynamic reactive power Optimizing Configuration System in an embodiment.

Detailed description of the invention

A kind of electrical power system dynamic reactive power Optimal Configuration Method, as it is shown in figure 1, comprise the following steps:

Step S110: receive the electrical network basic data of power system.Electrical network basic data includes the type of fault And the parameter of node.Node can be specifically transformer station etc., the parameter of node can include each node type, Number and title etc..The selection form of node is indefinite, specifically can need to select according to user to be seen Key node, it is also possible to be the node obtaining all power systems.Additionally, electrical network basic data may also include The place of all types of faults.Being appreciated that different according to the actual requirements, electrical network basic data may also include it His electric network data.Collect and arrange electrical network basic data, clearly requiring the input data of adjustment, in order to after Continuous data process.

Step S120: each node is carried out time-domain-simulation analysis according to the type of fault and the parameter of node, Node voltage data to corresponding node.The present embodiment can use electromechanical transient simulation software to power system Carry out the time-domain-simulation analysis of fault scanning.In conjunction with electrical network basic data, analyze the typical fault drawing electrical network, Write fault set and output collection file.Wherein, contingent multiple typical cases during fault set comprises power system Fault；Output collection is to need to select key node to be seen according to user, and these nodes have and can characterize The characteristic of system transient modelling voltage stability.

Specifically, by fault set and dynamo-electric Transient State Simulation Software PSD-BPA (the Power System of output collection write Department-Bonneville Power Administration, power system institute-Bang Naweier power office, Refer to a set of widely used power system analysis software instrument) input file, write typical fault set, institute The node output collection being concerned about and the txt file of the dynamic passive compensation corresponding with output collection node.At PSD-BPA Program is originally inputted in swi file and adds specific character string, is written fault set, output collection by judgement Position with the txt file of dynamic passive compensation.Consider the practical operating experiences of power system, electrical network peace Full property requires and the many factors such as user's request, and read failure is concentrated and specified fault, and with output card and In swi file, the most steady parameter card is spliced into interim swi file, and interim swi filename is set to fault name.

In programming software platform, call electromechanical transient simulation software selected fault set is carried out time-domain-simulation divide Analysis, specifically can be entered selected fault set by Matlab routine call power system analysis software PSD-BPA Row time-domain-simulation is analyzed, and the voltage condition after a failure of the node selected by record, obtain emulation is defeated Go out file to process, obtain node voltage data.Specifically emulation is obtained out file and carries out data process, Out file includes explanatory note and data two parts, obtains one group of effective clear data by programming.Out file The data of node voltage are that the setting according to cycle and step-length obtains, finally according to cycle, step-length and defeated The requirement of card release, obtains each node voltage data.

The present embodiment calls electromechanical transient simulation software each node is carried out time-domain-simulation analysis, simplify electricity The optimal reactive power allocation analysis of Force system, also assures that the reliability of analysis result simultaneously.Combine programming soft Part platform enhances coverage rate and the fineness of computational analysis, makes quickly to analyze, finds, plans that electrical network is weak Link and safety and stability leak are possibly realized, and economic and social benefit is notable.

Step S130: the Enhancement of Transient Voltage Stability assessment being calculated power system according to node voltage data refers to Mark, and the rate of qualified voltage of power system it is calculated according to Enhancement of Transient Voltage Stability evaluation index.

Enhancement of Transient Voltage Stability assessment is carried out, specifically after obtaining node voltage data by time-domain-simulation analysis The voltage dip severity index of each node and system can be obtained by processing each node voltage data. Before carrying out idle configuration, network is analyzed, under the given method of operation and fault set calculate system and Severity index drops in node voltage temporarily, finds out the weak link of electrical network, in order to follow-up to carry out reactive power excellent Change.

Wherein in an embodiment, Enhancement of Transient Voltage Stability evaluation index includes that the voltage dip of each node is tight Principal characteristic index and voltage qualification index, and the voltage dip severity index of power system and voltage qualified Property index.Step S130 includes step 132 and step 134.

Step 132: be calculated voltage dip severity index and the voltage of each node according to node voltage data Qualification index.Single failure lower node voltage dip severity index and electricity is carried out according to node voltage data The calculating of pressure qualification index.

Specifically, the voltage stability criterion of power system is, system centring point bus electricity in transient state and dynamic process Pressure fall temporarily is consistently less than transient voltage threshold value 1 (typically can be taken as 0.75pu, on the basis of rating average voltage) Time less than setting (typically taking 1 second), and dynamic process terminates rear 220kV and above Backbone point busbar voltage is not less than steady state voltage threshold value 2 and (typically can be taken as 0.9pu, with rating average voltage be Benchmark).Wherein in an embodiment, node voltage data include that node is when occurring different types of fault Dynamic electric voltage.Set two corresponding voltage thresholds, be respectively designated as the first voltage threshold and the second voltage Threshold value, uses V_L1And V_L2Representing, value is 0.9pu and 0.75pu.Step 132 specifically includes:

${\mathrm{Se}}_{i,j}=D_{i,j}+{\∫}_{t_0}^{t_0+{\mathrm{\Δt}}_1}(V_{L1}-V_{i,j}(t\left)\right)dt$

Wherein, V_i,j(t) for when break down i time node j the perunit value of dynamic electric voltage, Se_i,jRepresent node Voltage dip severity index, Fp_i,jRepresent the voltage qualification index of node.V_L1And V_L2Represent first respectively Voltage threshold and the second voltage threshold, and V_L1More than V_L2。D_i,jFor penalty function, p is penalty function D_i,jTo discontented The penalty value of afc voltage stability criterion, is set as 10 in the present embodiment.s_vFor node voltage less than the second voltage Threshold value V_L2Persistent period limit value, the concrete value that sets can be 1 second, t₀The time started for fault, Δ t₁、 Δt₂When being respectively fault i, node j dynamic electric voltage is less than the first voltage threshold V_L1, the second voltage threshold V_L2Hold The continuous time, V_i,j(∞) it is the voltage after dynamic process terminates.The present embodiment judges according to voltage stability criterion Fault condition lower node voltage is the most qualified, if voltage is qualified, value is 1, and voltage is defective, and value is 0, Judgment accuracy is high, easy and simple to handle reliably.

The dynamic V of node voltage_i,jT () can be obtained by time-domain-simulation, therefore, the voltage dip seriousness of node refers to Mark Se_i,jImplement and use the data that obtain of time-domain-simulation, wherein integral and calculating also can change sum formula into, I.e. corresponding computing formula can be equivalent to:

Se_i,j=D_i,j+Σk_i,j(t)(V_L1-V_i,j(t))Δt

$k_{i,j}\left(t\right)=\left\{\begin{array}{c}0,V_{i,j}\left(t\right)>V_{L1}\\ 1,V_{i,j}\left(t\right)\≤V_{L1}\end{array}\right.$

Wherein, Δ t is time-domain-simulation material calculation, k_i,jT under the conditions of ()=0 expression fault i, node j is in t Voltage is more than the first voltage threshold, k_i,jT ()=1 represents that voltage is less than the first voltage threshold.By the electricity to node The computing formula of pressure fall severity index temporarily deforms, and can improve data and calculate simplicity.

Step 134: voltage dip severity index and voltage qualification index according to each node are calculated electricity The Enhancement of Transient Voltage Stability evaluation index of Force system and voltage qualification index, and according to the voltage of power system Qualification index is calculated the rate of qualified voltage of power system.

The voltage qualification index expression of power system in institute's attention location system set of node, total joint that voltage is qualified Count.The rate of qualified voltage of power system represents voltage, and qualified node accounts for the percentage ratio of paid close attention to node total number. The joint that severity index needs to consider system backbone node and needs are paid close attention to drops in system voltage temporarily Point.Owing to each node significance level is different, different weights can be set.For each node in system, Need to consider multiple typical fault posterior nodal point voltage recovery situation.Wherein in an embodiment, step 134 specifically include:

$Fp=\underset{j=1}{\overset{N}{\Σ}}\underset{i=1}{\overset{M}{\Π}}{\mathrm{Fp}}_{i,j}$

$\mathrm{Pr}=\frac{Fp}{N}\×100\%$

Wherein, Se represents the voltage dip severity index of power system, and Fp represents that the voltage of power system closes Lattice index, Pr represents the rate of qualified voltage of power system；M is fault number, and N is node number, M Individual node voltage qualification index Fp under individual typical fault_Σ,jIt is that the faulty situation lower node j voltage of institute closes Lattice index Fp_i,jCompany's multiplication, its implication is that and if only if node j transient voltage under all typical faults Fp is met time the most qualified_Σ,j=1, if defective at certain fault lower node transient voltage, then it is assumed that and this node is temporary State voltage is defective, now Fp_Σ,j=0.w_iFor the weighted value of fault i, meet formulaFor node j Weighted value, weight value mainly considers following two aspects: the weight of (1) different electric pressure transformer stations node Want degree difference very big, transformer station's importance is divided by electric pressure；(2) power transformation of same electric pressure Its importance is weighed by substation capacity size between tiny node.

Step S140: in rate of qualified voltage less than when presetting qualified threshold value, obtain voltage and do not meet preset standard Node obtain both candidate nodes, and calculate the reactive-load compensation sensitivity of each both candidate nodes.

Preset the concrete value of qualified threshold value, and preset standard all can adjust according to practical situation.This enforcement Presetting qualified threshold value in example is 100%, and corresponding with step S130, preset standard is the voltage qualification of node Index is 1, and i.e. obtaining voltage qualification index is that the node of 1 is as both candidate nodes.Complete node and electric power After the calculating of the voltage dip severity index of system, calculate the rate of qualified voltage of power system.If voltage Qualification rate reaches 100%, then explanation is made without reactive-load compensation, can the stable joint of direct output power system Point is to check.

If rate of qualified voltage can not reach 100%, then choose the node that the defective situation of transient voltage occurs, make For needing to carry out the both candidate nodes of reactive-load compensation, in the input file of the dynamo-electric Transient State Simulation Software of write, again Carrying out the time-domain-simulation analysis of fault scanning, calculating each node increases after dynamic passive compensation power train The reactive-load compensation sensitivity of the voltage dip severity index of system, the reactive-load compensation sensitivity of each node constitutes nothing Merit compensates sensitivity matrix, as the dynamic reactive amount of capacity of follow-up selected node.The node that sensitivity is big The reactive compensation capacity increased is needed to be greater than the node that sensitivity is little.

Wherein in an embodiment, step S140 calculates the reactive-load compensation sensitivity of each both candidate nodes, bag Include step 142 and step 144.

Step 142: successively each both candidate nodes is increased the dynamic passive compensation value of identical capacity, and again calculates The voltage dip severity index of power system.The concrete value of dynamic passive compensation value is not unique, respectively Each both candidate nodes is increased dynamic passive compensation value, and recalculates corresponding both candidate nodes increase dynamic reactive benefit Repay the voltage dip severity index of power system after value.

Step 144: calculate corresponding both candidate nodes according to the variable quantity of the voltage dip severity index of power system Reactive-load compensation sensitivity.Particularly as follows:

${\mathrm{\λ}}_i=\frac{\mathrm{\Δ}Se}{\mathrm{\Δ}cap}=\frac{Se(k-1)-Se{\left(k\right)}_i}{\mathrm{\Δ}cap},(i=1,2,3......)$

Wherein, λ_iRepresenting the reactive-load compensation sensitivity of both candidate nodes i, Se (k-1) is for increasing dynamic passive compensation The voltage dip severity index of power system, Se (k) before value_iThe reactive-load compensation of Δ cap is individually increased for node i The voltage dip severity index of power system after capacity, Δ cap represents the dynamic passive compensation value of increase, Concrete value can according to circumstances be set as arbitrarily being more than the numerical value of 0, and representative value is the integer between 1-10.

Step S150: increase the dynamic passive compensation capacity of corresponding both candidate nodes according to reactive-load compensation sensitivity, Until rate of qualified voltage is for presetting qualified threshold value.By being calculated the candidate needing to install reactive-load compensation equipment After the reactive-load compensation sensitivity matrix of node, carry out idle according to the size of reactive-load compensation sensitivity further Distribute rationally.Wherein in an embodiment, step S150 includes step 152 and step 154.

Step 152: according to the size of reactive-load compensation sensitivity, increase corresponding candidate's joint with identical interval capacity The dynamic passive compensation capacity of point.Specifically, according to the result of calculation of reactive-load compensation sensitivity, candidate is saved Point is ranked up from small to large according to sensitivity, the dynamic reactive of the both candidate nodes installation that its medium sensitivity is minimum Compensation capacity is Δ cap Mvar, and the capacity (Δ cap) identical according to each node separation is sequentially overlapped each The dynamic passive compensation capacity of node.

Step 154: after detection increase dynamic passive compensation capacity, whether the rate of qualified voltage of power system is less than pre- If qualified threshold value；The most then return step S140, again obtain voltage and do not meet the node of preset standard and obtain To both candidate nodes, and calculate the reactive-load compensation sensitivity of each both candidate nodes.If it is not, then carry out step S160.

Specifically, the dynamic passive compensation of corresponding both candidate nodes is increased in the size according to reactive-load compensation sensitivity After capacity, again calculate the Enhancement of Transient Voltage Stability evaluation index of power system, and again by power system Rate of qualified voltage with preset qualified threshold ratio relatively.Equally as a example by default qualified threshold value is 100%, if voltage Qualification rate reaches 100%, illustrates that the setting of this reactive capability is rational, can carry out step S160.

If rate of qualified voltage does not still reach 100%, then the reactive-load compensation recalculating each both candidate nodes is sensitive Degree, is updated reactive-load compensation sensitivity matrix.Utilize the reactive-load compensation sensitivity matrix after updating again Increase the dynamic passive compensation capacity of corresponding both candidate nodes, calculate again after increasing dynamic passive compensation capacity The Enhancement of Transient Voltage Stability evaluation index of power system, and the rate of qualified voltage calculating power system compares. So circulation is until the rate of qualified voltage of power system is 100% so that system voltage reaches stable.

Additionally, for follow each time gone bad time, exportable result of calculation is to check, the present embodiment is fallen into a trap Calculation result includes: Enhancement of Transient Voltage Stability evaluation index (includes that severity index and voltage drop in system voltage temporarily Qualification rate) and current dynamic passive compensation scheme.

Step S160: obtain the OPTIMAL REACTIVE POWER of power system according to the dynamic passive compensation capacity of each both candidate nodes Compensation capacity scheme also exports.After the rate of qualified voltage of power system reaches to preset qualified threshold value, to each time The dynamic passive compensation capacity selecting node carries out sorting-out in statistics and obtains optimal reactive compensation volume solutions, is used for marking Show the dynamic passive compensation capacity to each both candidate nodes, for carrying out Dynamic reactive power optimization configuration.Export excellent nothing Merit compensation capacity scheme can be specifically to be shown by display screen so that staff watches, it is also possible to is output Being stored for data call to memorizer, concrete form does not limits.In output optimal reactive compensation capacity side After case, the stable node of the most exportable output power system is to check.

Above-mentioned electrical power system dynamic reactive power Optimal Configuration Method, introduces power system by the concept of sensitive analysis Dynamic reactive power optimization configuration among, by the dynamic characteristic of system in the case of sensitive analysis different faults, Realize distributing rationally of dynamic reactive according to rational criterion and index, can automatically calculate Electrical Power System Dynamic without The installation site of merit and amount of capacity, it is possible to systematically instruct the voltage association under the complicated idle configuration situation of electrical network Regulation and control system, effectively promotes receiving end electrical network and tackles the ability of various potential disturbances, improve Network Voltage Stability And management level.

Present invention also offers a kind of electrical power system dynamic reactive power Optimizing Configuration System, as in figure 2 it is shown, include Data acquisition module 110, time-domain-simulation module 120, qualification rate computing module 130, node screening module 140, Reactive compensation module 150 and optimization output module 160.

Data acquisition module 110 is for receiving the electrical network basic data of power system.Electrical network basic data includes The type of fault and the parameter of node.Node can be specifically transformer station etc., and the parameter of node can include respectively The type of node, number and title etc..The selection form of node is indefinite, can be specifically according to user's needs Select key node to be seen, it is also possible to be the node obtaining all power systems.Additionally, electrical network is basic Data may also include the place of all types of fault.It is appreciated that different according to the actual requirements, electrical network basis number According to may also include other electric network datas.Collect and arrange electrical network basic data, clearly requiring the input number of adjustment According to, in order to follow-up data process.

Time-domain-simulation module 120, for the type according to fault and the parameter of node, carries out time domain to each node Simulation analysis, obtains the node voltage data of corresponding node.The present embodiment can use electromechanical transient simulation soft Part carries out the time-domain-simulation analysis of fault scanning to power system.In conjunction with electrical network basic data, analyze and draw electricity The typical fault of net, writes fault set and output collection file.Wherein, may during fault set comprises power system The multiple typical faults occurred；Output collection is to need to select key node to be seen according to user, these joints Point has the characteristic that can characterize system transient modelling voltage stability.

Call electromechanical transient simulation software and each node is carried out time-domain-simulation analysis, simplify the nothing of power system Merit distributes analysis rationally, also assures that the reliability of analysis result simultaneously.Combine programming software platform to strengthen The coverage rate of computational analysis and fineness, make quickly to analyze, find, plan electrical network weak link and safety Stablizing leak to be possibly realized, economic and social benefit is notable.

Qualification rate computing module 130 for being calculated the transient voltage of power system according to node voltage data Stability assessment index, and the voltage conjunction of power system it is calculated according to Enhancement of Transient Voltage Stability evaluation index Lattice rate.

Enhancement of Transient Voltage Stability assessment is carried out, specifically after obtaining node voltage data by time-domain-simulation analysis The voltage dip severity index of each node and system can be obtained by processing each node voltage data. Before carrying out idle configuration, network is analyzed, under the given method of operation and fault set calculate system and Severity index drops in node voltage temporarily, finds out the weak link of electrical network, in order to follow-up to carry out reactive power excellent Change.

Wherein in an embodiment, Enhancement of Transient Voltage Stability evaluation index includes that the voltage dip of each node is tight Principal characteristic index and voltage qualification index, and the voltage dip severity index of power system and voltage qualified Property index.Qualification rate computing module 130 includes node indicator calculating unit and system index computing unit.

Node indicator calculating unit is serious for the voltage dip being calculated each node according to node voltage data Property index and voltage qualification index.Single failure lower node voltage dip is carried out tight according to node voltage data Principal characteristic index and the calculating of voltage qualification index.Wherein in an embodiment, node voltage data include The node dynamic electric voltage when there is different types of fault.Set two corresponding voltage thresholds, order respectively Entitled first voltage threshold and the second voltage threshold, use V_L1And V_L2Representing, value is 0.9pu and 0.75pu. Node indicator calculating unit according to node voltage data be calculated each node voltage dip severity index and Voltage qualification index specifically includes:

${\mathrm{Se}}_{i,j}=D_{i,j}+{\∫}_{t_0}^{t_0+{\mathrm{\Δt}}_1}(V_{L1}-V_{i,j}(t\left)\right)dt$

Wherein, V_i,j(t) for when break down i time node j the perunit value of dynamic electric voltage, Se_i,jRepresent node Voltage dip severity index, Fp_i,jRepresent the voltage qualification index of node.V_L1And V_L2Represent first respectively Voltage threshold and the second voltage threshold, and V_L1More than V_L2。D_i,jFor penalty function, p is penalty function D_i,jTo discontented The penalty value of afc voltage stability criterion, is set as 10 in the present embodiment.s_vFor node voltage less than the second voltage Threshold value V_L2Persistent period limit value, the concrete value that sets can be 1 second, t₀The time started for fault, Δ t₁、 Δt₂When being respectively fault i, node j dynamic electric voltage is less than the first voltage threshold V_L1, the second voltage threshold V_L2Hold The continuous time, V_i,j(∞) it is the voltage after dynamic process terminates.The present embodiment judges according to voltage stability criterion Fault condition lower node voltage is the most qualified, if voltage is qualified, value is 1, and voltage is defective, and value is 0, Judgment accuracy is high, easy and simple to handle reliably.

The dynamic V of node voltage_i,jT () can be obtained by time-domain-simulation, therefore, the voltage dip seriousness of node refers to Mark Se_i,jImplement and use the data that obtain of time-domain-simulation, wherein integral and calculating also can change sum formula into, I.e. corresponding computing formula can be equivalent to:

Se_i,j=D_i,j+Σk_i,j(t)(V_L1-V_i,j(t))Δt

$k_{i,j}\left(t\right)=\left\{\begin{array}{c}0,V_{i,j}\left(t\right)>V_{L1}\\ 1,V_{i,j}\left(t\right)\≤V_{L1}\end{array}\right.$

Wherein, Δ t is time-domain-simulation material calculation, k_i,jT under the conditions of ()=0 expression fault i, node j is in t Voltage is more than the first voltage threshold, k_i,jT ()=1 represents that voltage is less than the first voltage threshold.By the electricity to node The computing formula of pressure fall severity index temporarily deforms, and can improve data and calculate simplicity.

System index computing unit is for referring to according to voltage dip severity index and the voltage qualification of each node Mark is calculated Enhancement of Transient Voltage Stability evaluation index and the voltage qualification index of power system, and according to electricity The voltage qualification index of Force system is calculated the rate of qualified voltage of power system.

The voltage qualification index expression of power system in institute's attention location system set of node, total joint that voltage is qualified Count.The rate of qualified voltage of power system represents voltage, and qualified node accounts for the percentage ratio of paid close attention to node total number. The joint that severity index needs to consider system backbone node and needs are paid close attention to drops in system voltage temporarily Point.Owing to each node significance level is different, different weights can be set.For each node in system, Need to consider multiple typical fault posterior nodal point voltage recovery situation.Wherein in an embodiment, system Indicator calculating unit is calculated electricity according to voltage dip severity index and the voltage qualification index of each node The Enhancement of Transient Voltage Stability evaluation index of Force system and voltage qualification index, and according to the voltage of power system Qualification index is calculated the rate of qualified voltage of power system and specifically includes:

$Fp=\underset{j=1}{\overset{N}{\Σ}}\underset{i=1}{\overset{M}{\Π}}{\mathrm{Fp}}_{i,j}$

$\mathrm{Pr}=\frac{Fp}{N}\×100\%$

Wherein, Se represents the voltage dip severity index of power system, and Fp represents that the voltage of power system closes Lattice index, Pr represents the rate of qualified voltage of power system；M is fault number, and N is node number, M Individual node voltage qualification index Fp under individual typical fault_Σ,jIt is that the faulty situation lower node j voltage of institute closes Lattice index Fp_i,jCompany's multiplication, its implication is that and if only if node j transient voltage under all typical faults Fp is met time the most qualified_∑,j=1, if defective at certain fault lower node transient voltage, then it is assumed that and this node is temporary State voltage is defective, now Fp_∑,j=0.w_iFor the weighted value of fault i, meet formulaFor node j Weighted value, weight value mainly considers following two aspects: the weight of (1) different electric pressure transformer stations node Want degree difference very big, transformer station's importance is divided by electric pressure；(2) power transformation of same electric pressure Its importance is weighed by substation capacity size between tiny node.

Node screening module 140 is for when rate of qualified voltage is less than default qualified threshold value, obtaining voltage and be not inconsistent The node closing preset standard obtains both candidate nodes, and calculates the reactive-load compensation sensitivity of each both candidate nodes.

Preset the concrete value of qualified threshold value, and preset standard all can adjust according to practical situation.This enforcement In example preset qualified threshold value be 100%, accordingly, preset standard be the voltage qualification index of node be 1, I.e. obtaining voltage qualification index is that the node of 1 is as both candidate nodes.At the electricity completing node and power system After the calculating of pressure fall severity index temporarily, calculate the rate of qualified voltage of power system.If rate of qualified voltage reaches To 100%, then explanation is made without reactive-load compensation, can the stable node of direct output power system to look into See.

If rate of qualified voltage can not reach 100%, then choose the node that the defective situation of transient voltage occurs, make For needing to carry out the both candidate nodes of reactive-load compensation, in the input file of the dynamo-electric Transient State Simulation Software of write, again Carrying out the time-domain-simulation analysis of fault scanning, calculating each node increases after dynamic passive compensation power train The reactive-load compensation sensitivity of the voltage dip severity index of system, the reactive-load compensation sensitivity of each node constitutes nothing Merit compensates sensitivity matrix, as the dynamic reactive amount of capacity of follow-up selected node.The node that sensitivity is big The reactive compensation capacity increased is needed to be greater than the node that sensitivity is little.

Wherein in an embodiment, node screening module 140 includes node screening unit, node processing list Unit and Calculation of Sensitivity unit.

Node screening unit is for when rate of qualified voltage is less than default qualified threshold value, obtaining voltage and do not meet pre- The node of bidding standard obtains both candidate nodes.If rate of qualified voltage can not reach 100%, then choose and transient state electricity occurs Press the node of defective situation, as needing to carry out the both candidate nodes of reactive-load compensation.

Endpoint processing unit is used for increasing each both candidate nodes successively the dynamic passive compensation value of identical capacity, and Again calculate the voltage dip severity index of power system.The concrete value of dynamic passive compensation value is the most only One, respectively each both candidate nodes is increased dynamic passive compensation value, and it is dynamic to recalculate corresponding both candidate nodes increase The voltage dip severity index of power system after state reactive compensation values.

Calculation of Sensitivity unit is right for the variable quantity calculating of the voltage dip severity index according to power system Answer the reactive-load compensation sensitivity of both candidate nodes.Particularly as follows:

${\mathrm{\λ}}_i=\frac{\mathrm{\Δ}Se}{\mathrm{\Δ}cap}=\frac{Se(k-1)-Se{\left(k\right)}_i}{\mathrm{\Δ}cap},(i=1,2,3......)$

Wherein, λ_iRepresenting the reactive-load compensation sensitivity of both candidate nodes i, Se (k-1) is for increasing dynamic passive compensation The voltage dip severity index of power system, Se (k) before value_iThe reactive-load compensation of Δ cap is individually increased for node i The voltage dip severity index of power system after capacity, Δ cap represents the dynamic passive compensation value of increase, Concrete value can according to circumstances be set as arbitrarily being more than the numerical value of 0, and representative value is the integer between 1-10.

Reactive compensation module 150 for increasing the dynamic reactive of corresponding both candidate nodes according to reactive-load compensation sensitivity Compensation capacity, until rate of qualified voltage is for presetting qualified threshold value.Need reactive-load compensation is installed by being calculated After the reactive-load compensation sensitivity matrix of the both candidate nodes of equipment, further according to the size of reactive-load compensation sensitivity Carry out optimal reactive power allocation.Wherein in an embodiment, reactive compensation module 150 includes reactive-load compensation Unit and qualification rate detector unit.

Reactive compensation unit is for the size according to reactive-load compensation sensitivity, and it is right to increase with identical interval capacity Answer the dynamic passive compensation capacity of both candidate nodes.Specifically, according to the result of calculation of reactive-load compensation sensitivity, Both candidate nodes is ranked up from small to large according to sensitivity, the both candidate nodes installation that its medium sensitivity is minimum Dynamic passive compensation capacity is Δ cap Mvar, and the capacity (Δ cap) identical according to each node separation comes successively The dynamic passive compensation capacity of each node of superposition.

Qualification rate detector unit is the rate of qualified voltage of power system after detecting increase dynamic passive compensation capacity Whether less than presetting qualified threshold value；The most then control node screening module 140 again to obtain voltage and do not meet The node of preset standard obtains both candidate nodes, and calculates the reactive-load compensation sensitivity of each both candidate nodes；If it is not, Then control to optimize output module 160 and obtain power system according to the dynamic passive compensation capacity of each both candidate nodes Optimal reactive compensation volume solutions also exports.

Specifically, the dynamic passive compensation of corresponding both candidate nodes is increased in the size according to reactive-load compensation sensitivity After capacity, again calculate the Enhancement of Transient Voltage Stability evaluation index of power system, and again by power system Rate of qualified voltage with preset qualified threshold ratio relatively.Equally as a example by default qualified threshold value is 100%, if voltage Qualification rate reaches 100%, illustrates that the setting of this reactive capability is rational, the optimum of exportable power system Reactive compensation capacity scheme.

If rate of qualified voltage does not still reach 100%, then the reactive-load compensation recalculating each both candidate nodes is sensitive Degree, is updated reactive-load compensation sensitivity matrix.Utilize the reactive-load compensation sensitivity matrix after updating again Increase the dynamic passive compensation capacity of corresponding both candidate nodes, calculate again after increasing dynamic passive compensation capacity The Enhancement of Transient Voltage Stability evaluation index of power system, and the rate of qualified voltage calculating power system compares. So circulation is until the rate of qualified voltage of power system is 100% so that system voltage reaches stable.

Additionally, for follow each time gone bad time, exportable result of calculation is to check, the present embodiment is fallen into a trap Calculation result includes: Enhancement of Transient Voltage Stability evaluation index (includes that severity index and voltage drop in system voltage temporarily Qualification rate) and current dynamic passive compensation scheme.

Optimize output module 160 for obtaining power system according to the dynamic passive compensation capacity of each both candidate nodes Optimal reactive compensation volume solutions and export.When the rate of qualified voltage of power system reaches to preset qualified threshold value After, the dynamic passive compensation capacity of each both candidate nodes is carried out sorting-out in statistics and obtains optimal reactive compensation capacity side Case, for indicating the dynamic passive compensation capacity to each both candidate nodes, for carrying out Dynamic reactive power optimization configuration. Exporting excellent reactive compensation capacity scheme can be specifically to be shown so that staff watches by display screen, it is possible to To be that output is stored for data call to memorizer, concrete form does not limits.Mend in output OPTIMAL REACTIVE POWER After repaying volume solutions, the stable node of the most exportable output power system is to check.

Above-mentioned electrical power system dynamic reactive power Optimizing Configuration System, introduces power system by the concept of sensitive analysis Dynamic reactive power optimization configuration among, by the dynamic characteristic of system in the case of sensitive analysis different faults, Realize distributing rationally of dynamic reactive according to rational criterion and index, can automatically calculate Electrical Power System Dynamic without The installation site of merit and amount of capacity, it is possible to systematically instruct the voltage association under the complicated idle configuration situation of electrical network Regulation and control system, effectively promotes receiving end electrical network and tackles the ability of various potential disturbances, improve Network Voltage Stability And management level.

Embodiment described above only have expressed the several embodiments of the present invention, and it describes more concrete and detailed, But can not therefore be construed as limiting the scope of the patent.It should be pointed out that, for this area For those of ordinary skill, without departing from the inventive concept of the premise, it is also possible to make some deformation and change Entering, these broadly fall into protection scope of the present invention.Therefore, the protection domain of patent of the present invention should be with appended power Profit requires to be as the criterion.

Claims (10)

1. an electrical power system dynamic reactive power Optimal Configuration Method, it is characterised in that comprise the following steps:

Receiving the electrical network basic data of power system, described electrical network basic data includes type and the joint of fault The parameter of point；

Type and the parameter of described node according to described fault carry out time-domain-simulation analysis to each node, obtain The node voltage data of corresponding node；

The Enhancement of Transient Voltage Stability assessment being calculated described power system according to described node voltage data refers to Mark, and the rate of qualified voltage of described power system it is calculated according to described Enhancement of Transient Voltage Stability evaluation index；

In described rate of qualified voltage less than when presetting qualified threshold value, obtain voltage and do not meet the node of preset standard Obtain both candidate nodes, and calculate the reactive-load compensation sensitivity of each both candidate nodes；

The dynamic passive compensation capacity of corresponding both candidate nodes is increased according to described reactive-load compensation sensitivity, until institute Stating rate of qualified voltage is described default qualified threshold value；

Dynamic passive compensation capacity according to each both candidate nodes obtains the optimal reactive compensation of described power system to be held Amount scheme also exports.

Electrical power system dynamic reactive power Optimal Configuration Method the most according to claim 1, it is characterised in that Described Enhancement of Transient Voltage Stability evaluation index includes voltage dip severity index and the voltage qualification of each node Index, and the voltage dip severity index of described power system and voltage qualification index；Described basis Described node voltage data are calculated the Enhancement of Transient Voltage Stability evaluation index of described power system, and according to Described Enhancement of Transient Voltage Stability evaluation index is calculated the step of the rate of qualified voltage of described power system, bag Include following steps:

Voltage dip severity index and the voltage of each described node it is calculated according to described node voltage data Qualification index；

Voltage dip severity index and voltage qualification index according to each described node are calculated described electricity The Enhancement of Transient Voltage Stability evaluation index of Force system and voltage qualification index, and according to described power system Voltage qualification index is calculated the rate of qualified voltage of described power system.

Electrical power system dynamic reactive power Optimal Configuration Method the most according to claim 2, it is characterised in that Described node voltage data include the node dynamic electric voltage when there is different types of fault；According to described joint Point voltage data are calculated voltage dip severity index and the voltage qualification index of each described node, bag Include:

${\mathrm{Se}}_{i,j}=D_{i,j}+{\∫}_{t_0}^{t_0+{\mathrm{\Δt}}_1}(V_{L1}-V_{i,j}(t\left)\right)dt$

Wherein, V_i,j(t) for when break down i time node j the perunit value of dynamic electric voltage, Se_i,jRepresent node Voltage dip severity index, Fp_i,jRepresent the voltage qualification index of node；V_L1And V_L2Represent first respectively Voltage threshold and the second voltage threshold, and V_L1More than V_L2；D_i,jFor penalty function, p is penalty function D_i,jTo discontented The penalty value of afc voltage stability criterion, s_vFor node voltage less than the second voltage threshold V_L2Persistent period limit value, t₀The time started for fault, Δ t₁、Δt₂When being respectively fault i, node j dynamic electric voltage is less than the first voltage threshold Value V_L1, the second voltage threshold V_L2Persistent period, V_i,j(∞) it is the voltage after dynamic process terminates.

Electrical power system dynamic reactive power Optimal Configuration Method the most according to claim 3, it is characterised in that The described voltage dip severity index according to each described node and voltage qualification index are calculated described electricity The Enhancement of Transient Voltage Stability evaluation index of Force system and voltage qualification index, and according to described power system Voltage qualification index is calculated the rate of qualified voltage of described power system, including:

$Fp=\underset{j=1}{\overset{N}{\Σ}}\underset{i=1}{\overset{M}{\Π}}{\mathrm{Fp}}_{i,j}$

$\mathrm{Pr}=\frac{Fp}{N}\×100\%$

Wherein, Se represents the voltage dip severity index of power system, and Fp represents that the voltage of power system closes Lattice index, Pr represents the rate of qualified voltage of power system；M is fault number, and N is node number, w_iFor The weighted value of fault i, meets formula Weighted value for node j.

Electrical power system dynamic reactive power Optimal Configuration Method the most according to claim 2, it is characterised in that The step of the described reactive-load compensation sensitivity being calculated each both candidate nodes, comprises the following steps:

Successively each both candidate nodes is increased the dynamic passive compensation value of identical capacity, and again calculates described electric power The voltage dip severity index of system；

The variable quantity of the voltage dip severity index according to described power system calculates the nothing of corresponding both candidate nodes Merit compensates sensitivity, particularly as follows:

${\mathrm{\λ}}_i=\frac{\mathrm{\Δ}Se}{\mathrm{\Δ}cap}=\frac{Se(k-1)-Se{\left(k\right)}_i}{\mathrm{\Δ}cap},(i=1,2,3......)$

Wherein, λ_iRepresenting the reactive-load compensation sensitivity of both candidate nodes i, Se (k-1) is for increasing dynamic passive compensation The voltage dip severity index of power system, Se (k) before value_iThe reactive-load compensation of Δ cap is individually increased for node i The voltage dip severity index of power system after capacity, Δ cap represents the dynamic passive compensation value of increase.

Electrical power system dynamic reactive power Optimal Configuration Method the most according to claim 1, it is characterised in that The described dynamic passive compensation capacity increasing corresponding both candidate nodes according to described reactive-load compensation sensitivity, until institute Stating rate of qualified voltage is the described step presetting qualified threshold value, comprises the following steps:

According to the size of described reactive-load compensation sensitivity, increase corresponding both candidate nodes with identical interval capacity Dynamic passive compensation capacity；

After detection increase dynamic passive compensation capacity, whether the rate of qualified voltage of power system is less than presetting qualified threshold Value；

The most then return the node that described acquisition voltage do not meets preset standard and obtain both candidate nodes, and calculate The step of the reactive-load compensation sensitivity of each both candidate nodes；

If it is not, then carry out the described dynamic passive compensation capacity according to each both candidate nodes to obtain described power system Optimal reactive compensation volume solutions and the step that exports.

7. an electrical power system dynamic reactive power Optimizing Configuration System, it is characterised in that including:

Data acquisition module, for receiving the electrical network basic data of power system, described electrical network basic data bag Include the type of fault and the parameter of node；

Time-domain-simulation module, for the type according to described fault and the parameter of described node, enters each node Row time-domain-simulation is analyzed, and obtains the node voltage data of corresponding node；

Qualification rate computing module, for being calculated the temporary of described power system according to described node voltage data State voltage stability evaluation index, and it is calculated described electricity according to described Enhancement of Transient Voltage Stability evaluation index The rate of qualified voltage of Force system；

Node screening module, for when described rate of qualified voltage is less than default qualified threshold value, obtaining voltage not The node meeting preset standard obtains both candidate nodes, and calculates the reactive-load compensation sensitivity of each both candidate nodes；

Reactive compensation module, for increasing the dynamic nothing of corresponding both candidate nodes according to described reactive-load compensation sensitivity Merit compensation capacity, until described rate of qualified voltage is described default qualified threshold value；

Optimize output module, for obtaining described power train according to the dynamic passive compensation capacity of each both candidate nodes System optimal reactive compensation volume solutions and export.

Electrical power system dynamic reactive power Optimizing Configuration System the most according to claim 7, it is characterised in that Described Enhancement of Transient Voltage Stability evaluation index includes voltage dip severity index and the voltage qualification of each node Index, and the voltage dip severity index of described power system and voltage qualification index；Described qualified Rate computing module includes:

Node indicator calculating unit, for being calculated the electricity of each described node according to described node voltage data Pressure fall severity index and voltage qualification index temporarily；

System index computing unit, closes for the voltage dip severity index according to each described node and voltage Lattice index is calculated the Enhancement of Transient Voltage Stability evaluation index of described power system and voltage qualification refers to Mark, and it is qualified to be calculated the voltage of described power system according to the voltage qualification index of described power system Rate.

Electrical power system dynamic reactive power Optimizing Configuration System the most according to claim 8, it is characterised in that Described node screening module includes:

Node screening unit, for when described rate of qualified voltage is less than default qualified threshold value, obtaining voltage not The node meeting preset standard obtains both candidate nodes；

Endpoint processing unit, for successively each both candidate nodes being increased the dynamic passive compensation value of identical capacity, And again calculate the voltage dip severity index of described power system；

Calculation of Sensitivity unit, for the variable quantity of the voltage dip severity index according to described power system Calculate the reactive-load compensation sensitivity of corresponding both candidate nodes, particularly as follows:

${\mathrm{\λ}}_i=\frac{\mathrm{\Δ}Se}{\mathrm{\Δ}cap}=\frac{Se(k-1)-Se{\left(k\right)}_i}{\mathrm{\Δ}cap},(i=1,2,3......)$

Wherein, λ_iRepresenting the reactive-load compensation sensitivity of both candidate nodes i, Se (k-1) is for increasing dynamic passive compensation The voltage dip severity index of power system, Se (k) before value_iThe reactive-load compensation of Δ cap is individually increased for node i The voltage dip severity index of power system after capacity, Δ cap represents the dynamic passive compensation value of increase.

Electrical power system dynamic reactive power Optimizing Configuration System the most according to claim 7, it is characterised in that Described reactive compensation module includes:

Reactive compensation unit, for the size according to described reactive-load compensation sensitivity, with identical interval capacity Increase the dynamic passive compensation capacity of corresponding both candidate nodes；

Qualification rate detector unit, after detecting increase dynamic passive compensation capacity, the voltage of power system is qualified Whether rate is less than presetting qualified threshold value；The most then control described node screening module again to obtain voltage and be not inconsistent The node closing preset standard obtains both candidate nodes, and calculates the reactive-load compensation sensitivity of each both candidate nodes；If it is not, Then control described optimization output module and obtain described power train according to the dynamic passive compensation capacity of each both candidate nodes System optimal reactive compensation volume solutions and export.