CN105811401A - Auxiliary decision making method for provincial power grid accident emergency control - Google Patents

Abstract

The invention discloses an auxiliary decision making method for provincial power grid accident emergency control. The method comprises the following steps: automatically scanning and judging the property of a provincial power grid fault by a power system, obtaining provincial power grid operation data in the fault and building a power grid real-time power flow database, a security restriction database, a control measures database and a power grid sensitivity database; making systems of 220kV and above for a provincial power grid equivalent to a node-branch topology structure; with the minimum total control cost after a power grid accident as a target function and load balance, power flow security and a control variable as constraint conditions, building an auxiliary decision making model for provincial power grid accident emergency control; and determining a model, and obtaining an auxiliary decision making scheme for provincial power grid accident emergency control. The decision making method is provided from the provincial power grid particularity by combining the provincial power grid dispatching principle, the control requirements and the practical operating experience; the performability and the operability of the decision making scheme are improved; the workload when a dispatcher treats the power grid accident is reduced; the treatment strategy is optimized; and the timeliness of power grid accident treatment is ensured.

Description

Provincial power network accident emergency controls aid decision-making method

Technical field

The invention belongs to power system automatic field, be specifically related to a kind of provincial power network accident emergency and control aid decision Method.

Background technology

In recent years, along with national grid big operation architecture construction is pushed forward comprehensively, the whole nation extra-high voltage interconnected network add run-up If tradition relies on the drawback of dispatching of power netwoks pattern of power dispatcher's experience and the most progressively highlights, state, point, save three Levels Scheduling electrical network On-line security and stability analysis application is thus arisen at the historic moment.At present, on-line security and stability analysis system embedment a new generation D5000 scheduling Platform, based on the whole network real-time data synchronization and shared, periodical safety scanner uni anticipation mode analysis mechanisms, monitoring current electric grid Operation risk also provides practical aid decision, can find the weak spot that current electric grid runs in time, it is achieved electrical network preventive control.

After provincial power network accident, emergency flight control is the emergent control means under power grid accident or abnormal conditions, and it is the most auxiliary Decision-making is helped to directly affect the promptness and correctness controlled after power grid accident.Current on-line system is built from top to bottom, and accident should The anxious Decision Model that controls is set up based on extra-high voltage interconnected power grid, its control generation controlling provincial power network accident emergency Considering the most subtly and optimizing of the aspects such as valency, control device, stability control strategy, Hydro-Thermal Systems control characteristic, may cause Decision scheme enforceability is strong or the problem that controls high cost.

Summary of the invention

It is an object of the invention to provide a kind of for provincial power network accident, consider all types of unit output of electrical network and The provincial power network accident emergency of electrical network power reguirements controls aid decision-making method.

This provincial power network accident emergency that the present invention provides controls aid decision-making method, comprises the steps:

S1., when provincial power network has an accident, power system automatically scans and judges the character of provincial power network fault, and obtains Take provincial power network service data during fault, set up electrical network Real-time Power Flow, security restriction, control measure and electrical network sensitivity data Storehouse；

S2. provincial power network 220kV and system above are equivalent to node-branch road topological structure；

S3. to control total Least-cost after power grid accident as object function, by each Branch Power Flow safety, save between interconnection merit Rate balance, as primarily controlling requirement, with balancing the load, trend safety and control variable as constraints, specifically considers water, fire Electricity exert oneself adjustment, Hydropower Unit start and stop, thermoelectricity emergency shutdown, load fall for, small power supply exert oneself control, element legs switching, draw The throwing of limited load, unit Reactive-power control, capacitive reactance device is moved back, main shift adjusts, unit trip surely arranges the machine of cutting, cutting load factor, extreme misery Electricity adjustable extent and water power regimen factor, set up provincial power network accident emergency and control Decision Model；

S4. the provincial power network accident emergency set up step S3 controls Decision Model and solves, and obtains provincial electricity Online story event emergency flight control aid decision scheme.

Node described in step S2-branch road topological structure, for transmission line of electricity, contact become main transformer as branch road, with power plant, Between province, this side bus of interconnection is power supply node, with all transformer stations 220kV bus as load bus, by off line load, online Small power supply exerts oneself all equivalences to 220kV bus, formation node-branch road topological structure.

Object function described in step S3, is specially and uses following formula as object function:

$\begin{array}{cc}\min & F=\underset{v=1}{\overset{11}{\Σ}}{\mathrm{\β}}_v{C}_v\\ \left\{\begin{array}{c}{C}_1=\underset{i=1}{\overset{a}{\Σ}}\left|{\mathrm{\Δp}}_{s-i}\right|\\ C_2=\underset{j=1}{\overset{b}{\Σ}}\left|{\mathrm{\Δp}}_{h-j}\right|\\ C_3=\underset{i=1}{\overset{a}{\Σ}}\left|{\mathrm{\Δp}}_{s-i}\right|t_{s-i}\\ C_4=\underset{k=1}{\overset{m}{\Σ}}\left|{\mathrm{\Δp}}_{d-k}\right|\\ C_5=\underset{k=1}{\overset{m}{\Σ}}\left|{\mathrm{\Δp}}_{x-k}\right|\end{array}\right.& \left\{\begin{array}{c}{C}_6=\underset{e=1}{\overset{n}{\Σ}}{r}_e\\ C_7=\underset{j=1}{\overset{b}{\Σ}}\left|{\mathrm{\Δp}}_{h-j}\right|t_{h-j}\\ C_8=\underset{k=1}{\overset{m}{\Σ}}\left|{\mathrm{\Δp}}_{l-k}\right|\\ C_9=\underset{i=1}{\overset{a}{\Σ}}\left|{\mathrm{\Δq}}_{s-i}\right|+\underset{j=1}{\overset{b}{\Σ}}\left|{\mathrm{\Δq}}_{h-j}\right|\\ C_{10}=\underset{k=1}{\overset{m}{\Σ}}\left|{\mathrm{\Δq}}_{l-k}\right|\\ C_{11}=\underset{k=1}{\overset{m}{\Σ}}{\mathrm{\Δ\τ}}_k\end{array}\right.\end{array}$

C in formula₁～C₈Represent that exert oneself adjustments, thermoelectricity of water power is exerted oneself adjustments, Hydropower Unit start and stop, load confession, little electricity respectively Source exert oneself control, element legs switching, thermoelectricity emergency shutdown, draw the control cost of limited load eight class real power control strategy, C₉~ C₁₁Represent that the idle adjustment of exerting oneself of water, thermoelectricity, Substation Reactive-power Compensation adjust, transforming plant main transformer on-load transfers the files three class voltages respectively The control cost of control strategy；β₁～β₁₁Then represent the corresponding cost coefficient controlling cost；Δp_s-i、Δq_s-i(1≤i≤a)、 Δp_h-j、Δq_h-jHydroelectric power plant i, thermal power plant j are gained merit by (1≤j≤b) expression, idle adjustment amount of exerting oneself, Δ p_d-k、Δp_x-k、Δ p_l-k、Δq_l-k、τ_kLoad bus k is poured out burden with power, the meritorious regulation of exerting oneself of small power supply, draws limit to have by (1≤k≤m) expression Workload, reactive-load compensation adjustment amount and main transformer are transferred the files number；A, b, m are respectively the whole network water, thermoelectricity power supply node and load bus Quantity；t_s-iFor the start and stop unit number of units of hydroelectric power plant i, t_h-jEmergency shutdown number of units for j thermal power plant；r_eRepresent the whole network branch road, if Branch road e participates in throwing/cut then r_e=1, otherwise r_e=0.

Described cost coefficient, for using following formula calculation cost coefficient:

β_v=Tim_v×Eco_v

Tim in formula_vIt is the long-run cost rate of v regulation measure, Eco_vIt is the unit-economy one-tenth of v regulation measure This；Tim_vAnd Eco_vValue be to combine the data statistics run in reality Unit Commitment and regulation, scheduling operation to obtain.

Described Tim_vAnd Eco_vValue, shown in table 1 specific as follows:

The Tim of each regulation measure of table 1_v、Eco_v, reference value

Described in step S3 with balancing the load, trend safety and control variable as constraints, be specially use following formula make For constraints:

$(P_w^{*}-p_w)+\underset{i=1}{\overset{a}{\Σ}}{\mathrm{\Δp}}_{s-i}+\underset{j=1}{\overset{b}{\Σ}}{\mathrm{\Δp}}_{h-j}+\underset{k=1}{\overset{m}{\Σ}}{\mathrm{\Δp}}_{x-k}+\underset{k=1}{\overset{m}{\Σ}}{\mathrm{\Δp}}_{l-k}=0$

$D+{\Π}_G\·{\mathrm{\Δp}}_G+{\Π}_L\·{\mathrm{\Δp}}_L+{\Π}_R\·R\≤\stackrel{\‾}{D}$

$\underset{\‾}{U}\≤U+\mathrm{\Φ}\·\mathrm{\Δ}q+{\mathrm{\Φ}}_T\·{\mathrm{\Δ\τ}}_T\≤\stackrel{\‾}{U}$

$\underset{\‾}{p_e}\≤p_e\≤\stackrel{\‾}{p_e}$

$\left\{\begin{array}{c}\underset{\‾}{{\mathrm{\Δp}}_{s-i}}\≤{\mathrm{\Δp}}_{s-i}\≤\stackrel{\‾}{{\mathrm{\Δp}}_{s-i}}\\ -p_{h-j}\≤{\mathrm{\Δp}}_{h-j}\≤\stackrel{\‾}{{\mathrm{\Δp}}_{h-j}}\\ 0\≤{\mathrm{\Δp}}_{d-k}\≤\stackrel{\‾}{{\mathrm{\Δp}}_{d-k}}\\ \underset{\‾}{{\mathrm{\Δp}}_{x-k}}\≤{\mathrm{\Δp}}_{x-k}\≤\stackrel{\‾}{{\mathrm{\Δp}}_{x-k}}\\ 0\≤{\mathrm{\Δp}}_{l-k}\≤\stackrel{\‾}{{\mathrm{\Δp}}_{l-k}}\end{array}\right.$

$\left\{\begin{array}{c}\underset{\‾}{{\mathrm{\Δq}}_{s-i}}\≤{\mathrm{\Δq}}_{s-i}\≤\stackrel{\‾}{{\mathrm{\Δq}}_{s-i}}\\ \underset{\‾}{{\mathrm{\Δq}}_{h-j}}\≤{\mathrm{\Δq}}_{h-j}\≤\stackrel{\‾}{{\mathrm{\Δq}}_{h-j}}\\ \underset{\‾}{{\mathrm{\Δq}}_{l-k}}\≤{\mathrm{\Δq}}_{l-k}\≤\stackrel{\‾}{{\mathrm{\Δq}}_{l-k}}\end{array}\right.$

In formula, p^* _w、p_wRepresent respectively and send into meritorious planned value and instantaneous value when proparea outreaches winding thread；D=[d₁…d_n]、Representing the whole network effective power flow section currency respectively and control limit list, each limit that controls is by linear system Unite and electrical network after accident is determined through real-time N-1 scanning；ΔP_G、ΔP_L, R represent whole power node, load bus and branch road respectively List, wherein:

ΔP_G=Δ P_G=[Δ p_s-1…Δp_s-a…Δp_h-1…Δ_ph-b…]^T,

ΔP_L=[... (Δ p_d-k+Δp_x-k+Δp_l-k)…(-Δp_d-k+Δp_x-p+Δp_l-p)…]^T(k, p load each other falls Go out, pour node into, 1≤k ≠ p≤m), R=[r₁…r_e…r_n]^T；∏_G、∏_L、∏_RRepresent that power supply node-section is sensitive accordingly Degree, load bus-section sensitivity and branch road-section sensitivity matrix；U、 URepresent the whole network node voltage currency respectively And control upper and lower limit list, wherein U=[u_i…u_j…u_k…]^T；Φ、Φ_TIt is respectively each node of the whole network and injects idle-voltage Sensitivity, main transformer are transferred the files-voltage sensibility matrix, matrix element Φ_ij、Φ_T-ijRepresent that node j regulation unit is idle or main respectively The shift position Voltage Cortrol amount to node i；Δ q is node Reactive-power control amount list, Δ q=[Δ q_i…Δq_j…Δq_k…]^T, Mainly power supply node idle exert oneself, load bus reactive-load compensation etc.；Δτ_TTransfer the files list for load bus on-load voltage regulation main transformer, Δτ_T=[... Δ τ_k…]T。p_e、 p _e Represent that each branch road is gained merit and the steady upper and lower limit of heat； Δp _s-i 、 Δ p _H-j 、 Δp _x-k Represent the adjustable meritorious scope of exerting oneself of water, thermoelectricity and small power supply respectively,Represent load joint Point is pourable, can draw the meritorious limit of limit； Δq _s-i 、 Δq _h-j 、 Δq _l-k For corresponding water, thermoelectricity is idle Exert oneself and Substation Reactive-power Compensation adjustable extent； Δp _s-i ′For the real-time upper and lower spinning reserve amount of i factory, Avoiding the adjustable upper and lower limit of unit after vibrating area for i factory Hydropower Unit, λ is natural number；p_h-j、For having in real time of thermal power plant j Merit is exerted oneself and the combustion of unit minimum steady is exerted oneself, t_h-jFor its emergency shutdown number of units, η is natural number；If branch road e is to all out-of-limit The throwing of road f/cut sensitivity is negative value and then considers that branch road e participates in throwing/cutting, r_ePut 1, otherwise r_eSet to 0.

Solving described in step S4, comprises the steps:

Step 1) to step S3 set up provincial power network accident emergency control Decision Model parameter be optimized；

Step 2) to step 1) provincial power network accident emergency after the parameters optimization that obtains controls Decision Model, use Heuritic approach optimizing solves.

Step 1) described in optimization, including sensitivity matrix rarefaction, decoupling matrix depression of order and branch road merge depression of order.

This provincial power network accident emergency that the present invention provides controls aid decision-making method, former in conjunction with provincial power network scheduling Then, control to require and practical operating experiences, from the particularity of provincial power network, propose practical provincial power network accident emergency control Make online aid decision modeling method, compensate for the deficiency of existing model, improve enforceability and the operability of aid decision, Minimizing dispatcher is processed workload during power grid accident, optimizes accident treatment strategy, it is ensured that electric network accident management promptness, right Control the great practical value of operation of power networks risk.

Accompanying drawing explanation

Fig. 1 is the method flow diagram of the present invention.

Fig. 2 is the electric network composition schematic diagram of the specific embodiment of the present invention.

Detailed description of the invention

It is illustrated in figure 1 the method flow diagram of the present invention: this provincial power network accident emergency that the present invention provides controls auxiliary Help decision method, comprise the steps:

S1., when provincial power network has an accident, power system automatically scans and judges the character of provincial power network fault, and obtains Take provincial power network service data during fault, set up electrical network Real-time Power Flow, security restriction, control measure and electrical network sensitivity data Storehouse；

S2. provincial power network 220kV and system above are equivalent to node-branch road topological structure: become with transmission line of electricity, contact Main transformer is branch road, so that between power plant, province, this side bus of interconnection is as power supply node, with all transformer stations 220kV bus as load Node, exerts oneself all equivalences to 220kV bus, formation node-branch road topological structure by off line load, online small power supply.

S3. to control total Least-cost after power grid accident as object function, by each Branch Power Flow safety, save between interconnection merit Rate balance, as primarily controlling requirement, with balancing the load, trend safety and control variable as constraints, specifically considers water, fire Electricity exert oneself adjustment, Hydropower Unit start and stop, thermoelectricity emergency shutdown, load fall for, small power supply exert oneself control, element legs switching, draw The throwing of limited load, unit Reactive-power control, capacitive reactance device is moved back, main shift adjusts, unit trip surely arranges the machine of cutting, cutting load factor, extreme misery Electricity adjustable extent and water power regimen factor, set up provincial power network accident emergency and control Decision Model；

Object function, is specially and uses following formula as object function:

$\begin{array}{cc}\min & F=\underset{v=1}{\overset{11}{\Σ}}{\mathrm{\β}}_v{C}_v\\ \left\{\begin{array}{c}{C}_1=\underset{i=1}{\overset{a}{\Σ}}\left|{\mathrm{\Δp}}_{s-i}\right|\\ C_2=\underset{j=1}{\overset{b}{\Σ}}\left|{\mathrm{\Δp}}_{h-j}\right|\\ C_3=\underset{i=1}{\overset{a}{\Σ}}\left|{\mathrm{\Δp}}_{s-i}\right|t_{s-i}\\ C_4=\underset{k=1}{\overset{m}{\Σ}}\left|{\mathrm{\Δp}}_{d-k}\right|\\ C_5=\underset{k=1}{\overset{m}{\Σ}}\left|{\mathrm{\Δp}}_{x-k}\right|\end{array}\right.& \left\{\begin{array}{c}{C}_6=\underset{e=1}{\overset{n}{\Σ}}{r}_e\\ C_7=\underset{j=1}{\overset{b}{\Σ}}\left|{\mathrm{\Δp}}_{h-j}\right|t_{h-j}\\ C_8=\underset{k=1}{\overset{m}{\Σ}}\left|{\mathrm{\Δp}}_{l-k}\right|\\ C_9=\underset{i=1}{\overset{a}{\Σ}}\left|{\mathrm{\Δq}}_{s-i}\right|+\underset{j=1}{\overset{b}{\Σ}}\left|{\mathrm{\Δq}}_{h-j}\right|\\ C_{10}=\underset{k=1}{\overset{m}{\Σ}}\left|{\mathrm{\Δq}}_{l-k}\right|\\ C_{11}=\underset{k=1}{\overset{m}{\Σ}}{\mathrm{\Δ\τ}}_k\end{array}\right.\end{array}$

C in formula₁～C₈Represent that exert oneself adjustments, thermoelectricity of water power is exerted oneself adjustments, Hydropower Unit start and stop, load confession, little electricity respectively Source exert oneself control, element legs switching, thermoelectricity emergency shutdown, draw the control cost of limited load eight class real power control strategy, C₉~ C₁₁Represent that the idle adjustment of exerting oneself of water, thermoelectricity, Substation Reactive-power Compensation adjust, transforming plant main transformer on-load transfers the files three class voltages respectively The control cost of control strategy；β₁～β₁₁Then represent the corresponding cost coefficient controlling cost；Δp_s-i、Δq_s-i(1≤i≤a)、 Δp_h-j、Δq_h-jHydroelectric power plant i, thermal power plant j are gained merit by (1≤j≤b) expression, idle adjustment amount of exerting oneself, Δ p_d-k、Δp_x-k、Δ p_l-k、Δq_l-k、τ_kLoad bus k is poured out burden with power, the meritorious regulation of exerting oneself of small power supply, draws limit to have by (1≤k≤m) expression Workload, reactive-load compensation adjustment amount and main transformer are transferred the files number；A, b, m are respectively the whole network water, thermoelectricity power supply node and load bus Quantity；t_s-iFor the start and stop unit number of units of hydroelectric power plant i, t_h-jEmergency shutdown number of units for j thermal power plant；r_eRepresent the whole network branch road, if Branch road e participates in throwing/cut then r_e=1, otherwise r_e=0.

Wherein, cost coefficient can use following formula to calculate:

β_v=Tim_v×Eco_v

Tim in formula_vIt is the long-run cost rate of v regulation measure, Eco_vIt is the unit-economy one-tenth of v regulation measure This；Tim_vAnd Eco_vValue be to combine the data statistics run in reality Unit Commitment and regulation, scheduling operation to obtain.

Concrete, Tim_vAnd Eco_vValue as shown in table 1 below:

The Tim of each regulation measure of table 1_v、Eco_v, reference value

And described in this step with balancing the load, trend safety and control variable as constraints, be specially use under Formula is as constraints:

$(P_w^{*}-p_w)+\underset{i=1}{\overset{a}{\Σ}}{\mathrm{\Δp}}_{s-i}+\underset{j=1}{\overset{b}{\Σ}}{\mathrm{\Δp}}_{h-j}+\underset{k=1}{\overset{m}{\Σ}}{\mathrm{\Δp}}_{x-k}+\underset{k=1}{\overset{m}{\Σ}}{\mathrm{\Δp}}_{l-k}=0$

$D+{\Π}_G\·{\mathrm{\Δp}}_G+{\Π}_L\·{\mathrm{\Δp}}_L+{\Π}_R\·R\≤\stackrel{\‾}{D}$

$\underset{\‾}{U}\≤U+\mathrm{\Φ}\·\mathrm{\Δ}q+{\mathrm{\Φ}}_T\·{\mathrm{\Δ\τ}}_T\≤\stackrel{\‾}{U}$

$\underset{\‾}{p_e}\≤p_e\≤\stackrel{\‾}{p_e}$

$\left\{\begin{array}{c}\underset{\‾}{{\mathrm{\Δp}}_{s-i}}\≤{\mathrm{\Δp}}_{s-i}\≤\stackrel{\‾}{{\mathrm{\Δp}}_{s-i}}\\ -p_{h-j}\≤{\mathrm{\Δp}}_{h-j}\≤\stackrel{\‾}{{\mathrm{\Δp}}_{h-j}}\\ 0\≤{\mathrm{\Δp}}_{d-k}\≤\stackrel{\‾}{{\mathrm{\Δp}}_{d-k}}\\ \underset{\‾}{{\mathrm{\Δp}}_{x-k}}\≤{\mathrm{\Δp}}_{x-k}\≤\stackrel{\‾}{{\mathrm{\Δp}}_{x-k}}\\ 0\≤{\mathrm{\Δp}}_{l-k}\≤\stackrel{\‾}{{\mathrm{\Δp}}_{l-k}}\end{array}\right.$

$\left\{\begin{array}{c}\underset{\‾}{{\mathrm{\Δq}}_{s-i}}\≤{\mathrm{\Δq}}_{s-i}\≤\stackrel{\‾}{{\mathrm{\Δq}}_{s-i}}\\ \underset{\‾}{{\mathrm{\Δq}}_{h-j}}\≤{\mathrm{\Δq}}_{h-j}\≤\stackrel{\‾}{{\mathrm{\Δq}}_{h-j}}\\ \underset{\‾}{{\mathrm{\Δq}}_{l-k}}\≤{\mathrm{\Δq}}_{l-k}\≤\stackrel{\‾}{{\mathrm{\Δq}}_{l-k}}\end{array}\right.$

In formula, p^* _w、p_wRepresent respectively and send into meritorious planned value and instantaneous value when proparea outreaches winding thread；D=[d₁…d_n]、Representing the whole network effective power flow section currency respectively and control limit list, each limit that controls is by linear system Unite and electrical network after accident is determined through real-time N-1 scanning；ΔP_G、ΔP_L, R represent whole power node, load bus and branch road respectively List, wherein:

ΔP_G=Δ P_G=[Δ_ps-1…Δ_ps-a…Δp_h-1…Δ_ph-b…]^T,

ΔP_L=[... (Δ p_d-k+Δp_x-k+Δp_l-k)…(-Δp_d-k+Δp_x-p+Δp_l-p)…]^T(k, p load each other falls Go out, pour node into, 1≤k ≠ p≤m), R=[r₁…r_e…r_n]^T；∏_G、∏_L、∏_RRepresent that power supply node-section is sensitive accordingly Degree, load bus-section sensitivity and branch road-section sensitivity matrix；U、 URepresent the whole network node voltage currency respectively And control upper and lower limit list, wherein U=[u_i…u_j…u_k…]^T；φ、φ_TIt is respectively each node of the whole network and injects idle-voltage Sensitivity, main transformer are transferred the files-voltage sensibility matrix, matrix element Φ_ij、Φ_T-ijRepresent that node j regulation unit is idle or main respectively The shift position Voltage Cortrol amount to node i；Δ q is node Reactive-power control amount list, Δ q=[Δ q_i…Δq_j…Δq_k…]^T, Mainly power supply node idle exert oneself, load bus reactive-load compensation etc.；Δτ_TTransfer the files list for load bus on-load voltage regulation main transformer, Δτ_T=[... Δ τ_k…]T。p_e、 p _e Represent that each branch road is gained merit and the steady upper and lower limit of heat； Δp _s-i 、 Δ p _h-j 、 Δp _x-k Represent the adjustable meritorious scope of exerting oneself of water, thermoelectricity and small power supply respectively,Represent load joint Point is pourable, can draw the meritorious limit of limit； Δq _s-i 、 Δq _h-j 、 Δq _l-k For corresponding water, thermoelectricity is idle Exert oneself and Substation Reactive-power Compensation adjustable extent； Δp _s-i ' for the real-time upper and lower spinning reserve amount of i factory, Avoiding the adjustable upper and lower limit of unit after vibrating area for i factory Hydropower Unit, λ is natural number；p_h-j、For having in real time of thermal power plant j Merit is exerted oneself and the combustion of unit minimum steady is exerted oneself, t_h-jFor its emergency shutdown number of units, η is natural number；If branch road e is to all out-of-limit The throwing of road f/cut sensitivity is negative value and then considers that branch road e participates in throwing/cutting, r_ePut 1, otherwise r_eSet to 0.

S4. the provincial power network accident emergency set up step S3 controls Decision Model and solves, and obtains provincial electricity Online story event emergency flight control aid decision scheme.

Owing to the inventive method enriches control measure type, but also improve model ginseng, variable dimension simultaneously, cause asking Solve difficulty, the time that calculates increases；For ensureing its online real-time, need when Practical Calculation to carry out following parameter optimization.

(1) sensitivity matrix rarefaction: trend security constraint formula is the important step of model solution, wherein sensitivity matrix According to electrical network scale up to the most thousands of hundreds of rank, serious restriction calculates speed.In actual electric network, to a certain out-of-limit section, only The minority degree of coupling is high, electrical distance is short node, branch road are higher to its sensitivity, and remaining great deal of nodes, branch road are to its sensitivity The lowest even zero and keep relative stability, this category node, branch road are not involved in controlling substantially.To this, at initial sensitivity matrix In sorted by sensitivity or set sensitivity threshold, be zero by a large number number being put in relatively low sensitivity, improved by rarefaction Solution efficiency.Sensitivity threshold is set to 5% by suggestion, and this measure can improve the matrix degree of rarefication of about 70%, and does not affect control Effect.

(2) decoupling matrix depression of order: owing to node-section and branch road-section sensitivity are closely related with electrical distance, will system Each by stages relatively low for sensitivity, according to the subregion such as electric pressure, electric network composition, is solved by system interior joint, branch road and sensitivity matrix Coupling, forms the node of multiple low order, branch road sensitivity matrix calculates, and can be effectively improved calculating speed.According to this province electrical network Structure, will contact more weak, power supply point less 220kV section electrical network as decoupling subregion with major network, and set up collaborative with major network solves The sensitivity matrix of coupling depression of order.

(3) branch road merges: for double-circuit lines on the same pole, if circuit is included into N-1 category, and this line parameter circuit value phase with jumping Closely, branch road can be modeled merge, it is achieved model reduction.

After model optimization completes, use heuritic approach that the model after optimizing is carried out optimizing and solve, draw final Provincial power network accident emergency control aid decision scheme.

This provincial power network accident emergency provided the present invention below in conjunction with a concrete example controls aid decision side Method is further described:

As in figure 2 it is shown, be 500kV transformer station for certain provincial power network maintenance mode partial structurtes schematic diagram: A～D, a～n is 220kV transformer station, cd line interruption maintenance, s1, s2, h1, h2 are respectively four water, thermal power plants.Take peak load tide one day in summer Flow section, Aa-I, Aa-II line three fault forever of simulation parallel lines on same tower is jumped together, and accident trend is directed respectively into on-line system and basis Model carries out simulation calculation, and system core index see table 2.

Key index value after table 2 accident

According to upper table, after accident, Be line, eb line effective power flow seriously exceed the steady limit value of heat, and A becomes off line meritorious slight overload, But A becomes the+B off line section of change not to be transshipped；In addition 220kV minimum voltage gets over lower limit, and low-limit frequency is the most out-of-limit.

To this, on-line system aid decision (decision-making 1) and this model-aided decision-making (decision-making 2) see table 3,4.If after accident Each control measure start simultaneously at execution immediately, and according to table 1 above, two class aid decision implementation effect contrasts see table 5.

Table 3 on-line system state of emergency aid decision (decision-making 1)

4 model accident emergency aid decisions (decision-making 2) of table

Table 5 decision-making 1, the contrast of decision-making 2 key index

Knowable to table 3～table 5, two class aid decisions all can reach finally to control target.Decision-making 1 is gained merit by adjusting power plant Exerting oneself, draw limited load to control effective power flow, adjustment power plant is idle exerts oneself and Substation Reactive-power Compensation control voltage, master control cost F reaches 331360, wherein draws limited load 219MW, heavy overload circuit more limit to reach 17.4min；Decision-making 2 is except using above-mentioned arranging Execute, additionally use water power start, thermoelectricity emergency shutdown, fall load, adjust small power supply exert oneself, put into the measures such as extension wire, always Control cost F and reach 99508.8, wherein draw limited load 35MW, heavy overload circuit more limit to reach 0.1min.Contrast shows, certainly The master control cost of plan 2, draw limited load amount, heavy overload circuit to get over limit all far below decision-making 1, electrical network thing can be substantially reduced Therefore risk and impact.Visible, when electrical network major accident, complicated abnormal trend, this model accident emergency aid decision can reach To controlling requirement, and control efficiency, quality are higher.

Claims (8)

1. provincial power network accident emergency controls an aid decision-making method, comprises the steps:

S1. when provincial power network has an accident, power system scans and judges the character of provincial power network fault automatically, obtains fault Time provincial power network service data, set up electrical network Real-time Power Flow, security restriction, control measure and electrical network sensitivity data storehouse；

S2. provincial power network 220kV and system above are equivalent to node-branch road topological structure；

S3. to control total Least-cost after power grid accident as object function, by each Branch Power Flow safety, save between dominant eigenvalues put down Weighing apparatus, as primarily controlling requirement, with balancing the load, trend safety and control variable as constraints, specifically considers that water, thermoelectricity go out Power adjustment, Hydropower Unit start and stop, thermoelectricity emergency shutdown, load fall for, small power supply exert oneself controls, element legs switching, draw limit to bear Lotus, unit Reactive-power control, capacitive reactance device throw move back, main shift adjust, unit trip surely arrange the machine of cutting, cutting load factor, Hydro-Thermal Systems can Tune scope and water power regimen factor, set up provincial power network accident emergency and control Decision Model；

S4. the provincial power network accident emergency set up step S3 controls Decision Model and solves, and obtains provincial power network thing Therefore emergency flight control aid decision scheme.

Provincial power network accident emergency the most according to claim 1 controls aid decision-making method, it is characterised in that step S2 institute The node stated-branch road topological structure, for becoming main transformer as branch road with transmission line of electricity, contact, with power plant, save between this side of interconnection female Line is power supply node, with all transformer stations 220kV bus as load bus, exerts oneself the most equivalent by off line load, online small power supply To 220kV bus, form node-branch road topological structure.

Provincial power network accident emergency the most according to claim 1 and 2 controls aid decision-making method, it is characterised in that step S3 Described object function is:

$\begin{array}{cc}\min & F=\underset{v=1}{\overset{11}{\Σ}}{\mathrm{\β}}_v{C}_v\end{array}$

$\left\{\begin{array}{c}{C}_1=\underset{i=1}{\overset{a}{\Σ}}\left|{\mathrm{\Δp}}_{s-i}\right|\\ C_2=\underset{j=1}{\overset{b}{\Σ}}\left|{\mathrm{\Δp}}_{h-j}\right|\\ C_3=\underset{i=1}{\overset{a}{\Σ}}\left|{\mathrm{\Δp}}_{s-i}\right|t_{s-i}\\ C_4=\underset{k=1}{\overset{m}{\Σ}}\left|{\mathrm{\Δp}}_{d-k}\right|\\ C_5=\underset{k=1}{\overset{m}{\Σ}}\left|{\mathrm{\Δp}}_{x-k}\right|\end{array}\right.\left\{\begin{array}{c}{C}_6=\underset{e=1}{\overset{n}{\Σ}}{r}_e\\ C_7=\underset{j=1}{\overset{b}{\Σ}}\left|{\mathrm{\Δp}}_{h-j}\right|t_{h-j}\\ C_8=\underset{k=1}{\overset{m}{\Σ}}\left|{\mathrm{\Δp}}_{l-k}\right|\\ C_9=\underset{i=1}{\overset{a}{\Σ}}\left|{\mathrm{\Δq}}_{s-i}\right|+\underset{j=1}{\overset{b}{\Σ}}\left|{\mathrm{\Δq}}_{h-j}\right|\\ C_{10}=\underset{k=1}{\overset{m}{\Σ}}\left|{\mathrm{\Δq}}_{l-k}\right|\\ C_{11}=\underset{k=1}{\overset{m}{\Σ}}{\mathrm{\Δ\τ}}_k\end{array}\right.$

C in formula₁～C₈Represent that water power exert oneself adjustments, Hydropower Unit start and stop, load confession, small power supply of adjustments, thermoelectricity of exerting oneself goes out respectively Power control, element legs switching, thermoelectricity emergency shutdown, draw the control cost of limited load eight class real power control strategy, C₉～C₁₁Point Do not represent that the idle adjustment of exerting oneself of water, thermoelectricity, Substation Reactive-power Compensation adjust, transforming plant main transformer on-load transfers the files three class Control of Voltage The control cost of strategy；β₁～β₁₁Then represent the corresponding cost coefficient controlling cost；Δp_s-i、Δq_s-i(1≤i≤a)、Δ p_h-j、Δq_h-jHydroelectric power plant i, thermal power plant j are gained merit by (1≤j≤b) expression, idle adjustment amount of exerting oneself, Δ p_d-k、Δp_x-k、Δ p_l-k、Δq_l-k、τ_kLoad bus k is poured out burden with power, the meritorious regulation of exerting oneself of small power supply, draws limit to have by (1≤k≤m) expression Workload, reactive-load compensation adjustment amount and main transformer are transferred the files number；A, b, m are respectively the whole network water, thermoelectricity power supply node and load bus Quantity；t_s-iFor the start and stop unit number of units of hydroelectric power plant i, t_h-jEmergency shutdown number of units for j thermal power plant；r_eRepresent the whole network branch road, if Branch road e participates in throwing/cut then r_e=1, otherwise r_e=0。

Provincial power network accident emergency the most according to claim 3 controls aid decision-making method, it is characterised in that described generation Valency coefficient, for using following formula calculation cost coefficient:

β_v=Tim_v×Eco_v

Tim in formula_vIt is the long-run cost rate of v regulation measure, Eco_vIt it is the unit-economy cost of v regulation measure； Tim_vAnd Eco_vValue be to combine the data statistics run in reality Unit Commitment and regulation, scheduling operation to obtain.

Provincial power network accident emergency the most according to claim 4 controls aid decision-making method, it is characterised in that described Tim_v And Eco_vValue, shown in table 1 specific as follows:

Provincial power network accident emergency the most according to claim 1 and 2 controls aid decision-making method, it is characterised in that step S3 Described with balancing the load, trend safety and control variable as constraints, it is specially and uses following formula as constraints:

$(P_w^{*}-p_w)+\underset{i=1}{\overset{a}{\Σ}}{\mathrm{\Δp}}_{s-i}+\underset{j=1}{\overset{b}{\Σ}}{\mathrm{\Δp}}_{h-j}+\underset{k=1}{\overset{m}{\Σ}}{\mathrm{\Δp}}_{x-k}+\underset{k=1}{\overset{m}{\Σ}}{\mathrm{\Δp}}_{l-k}=0$

$D+{\Π}_G\·{\mathrm{\Δp}}_G+{\Π}_L\·{\mathrm{\Δp}}_L+{\Π}_R\·R\≤\stackrel{\‾}{D}$

$\underset{\‾}{U}\≤U+\mathrm{\Φ}\·\mathrm{\Δ}q+{\mathrm{\Φ}}_T\·{\mathrm{\Δ\τ}}_T\≤\stackrel{\‾}{U}$

$\underset{\‾}{p_e}\≤p_e\≤\stackrel{\‾}{p_e}$

$\left\{\begin{array}{c}\underset{\‾}{{\mathrm{\Δp}}_{s-i}}\≤{\mathrm{\Δp}}_{s-i}\≤\stackrel{\‾}{{\mathrm{\Δp}}_{s-i}}\\ -p_{h-j}\≤{\mathrm{\Δp}}_{h-j}\≤\stackrel{\‾}{{\mathrm{\Δp}}_{h-j}}\\ 0\≤{\mathrm{\Δp}}_{d-k}\≤\stackrel{\‾}{{\mathrm{\Δp}}_{d-k}}\\ \underset{\‾}{{\mathrm{\Δp}}_{x-k}}\≤{\mathrm{\Δp}}_{x-k}\≤\stackrel{\‾}{{\mathrm{\Δp}}_{x-k}}\\ 0\≤{\mathrm{\Δp}}_{l-k}\≤\stackrel{\‾}{{\mathrm{\Δp}}_{l-k}}\end{array}\right.$

$\left\{\begin{array}{c}\underset{\‾}{{\mathrm{\Δq}}_{s-i}}\≤{\mathrm{\Δq}}_{s-i}\≤\stackrel{\‾}{{\mathrm{\Δq}}_{s-i}}\\ \underset{\‾}{{\mathrm{\Δq}}_{h-j}}\≤{\mathrm{\Δq}}_{h-j}\≤\stackrel{\‾}{{\mathrm{\Δq}}_{h-j}}\\ \underset{\‾}{{\mathrm{\Δq}}_{l-k}}\≤{\mathrm{\Δq}}_{l-k}\≤\stackrel{\‾}{{\mathrm{\Δq}}_{l-k}}\end{array}\right.$

In formula, p^* _w、p_wRepresent respectively and send into meritorious planned value and instantaneous value when proparea outreaches winding thread；D=[d₁…d_n]、Representing the whole network effective power flow section currency respectively and control limit list, each limit that controls is by linear system Unite and electrical network after accident is determined through real-time N-1 scanning；ΔP_G、ΔP_L, R represent whole power node, load bus and branch road respectively List, wherein:

ΔP_G=ΔP_G=[Δp_s-1…Δp_s-a…Δp_h-1…Δp_h-b…]^T,

ΔP_L=[…(Δp_d-k+Δp_x-k+Δp_l-k)…(-Δp_d-k+Δp_x-p+Δp_l-p)…]^T(k, p load each other is poured out, is poured into Node, 1≤k ≠ p≤m), R=[r₁…r_e…r_n]^T；∏_G、∏_L、∏_RRepresent power supply node-section sensitivity, load joint accordingly Point-section sensitivity and branch road-section sensitivity matrix；U、 URepresent the whole network node voltage currency and control thereof respectively Upper and lower limit list, wherein U=[u_i…u_j…u_k…]^T；Φ、Φ_TIt is respectively each node of the whole network and injects idle-voltage sensibility, master Modified tone shelves-voltage sensibility matrix, matrix element Φ_ij、Φ_T-ijUnit is idle or main shift pair to represent node j regulation respectively The Voltage Cortrol amount of node i；Δ q is node Reactive-power control amount list, Δ q=[Δ q_i…Δq_j…Δq_k…]^T, mainly electricity Source node is idle exerts oneself, load bus reactive-load compensation etc.；Δτ_TTransfer the files list for load bus on-load voltage regulation main transformer, Δ τ_T=[... Δτ_k…]T。p_e、 p _e Represent that each branch road is gained merit and the steady upper and lower limit of heat； Δp _s-i 、 Δp _h-j 、 Δp _x-k Represent the adjustable meritorious scope of exerting oneself of water, thermoelectricity and small power supply respectively,Represent load bus pourable, The meritorious limit of limit can be drawn； Δq _s-i 、Δq_h-j、 Δq _l-k For corresponding water, thermoelectricity is idle exerts oneself and power transformation Stand reactive-load compensation adjustable extent； Δp _s-i ′For the real-time upper and lower spinning reserve amount of i factory,For i factory water power Unit avoids the adjustable upper and lower limit of the unit after vibrating area, and λ is natural number；p_h-j、Exert oneself for thermal power plant j in real time meritorious and The combustion of its unit minimum steady is exerted oneself, t_h-jFor its emergency shutdown number of units, η is natural number；If branch road e to the throwing of all out-of-limit branch road f/ Cut sensitivity be negative value then consider branch road e participate in throw/cut, r_ePut 1, otherwise r_eSet to 0.

Provincial power network accident emergency the most according to claim 1 and 2 controls aid decision-making method, it is characterised in that step S4 Described solves, and comprises the steps:

Step 1) to step S3 set up provincial power network accident emergency control Decision Model parameter be optimized；

Step 2) to step 1) provincial power network accident emergency after the parameters optimization that obtains controls Decision Model, use and inspire Formula algorithm optimizing solves.

Provincial power network accident emergency the most according to claim 7 controls aid decision-making method, it is characterised in that step 1) institute The optimization stated, merges depression of order including sensitivity matrix rarefaction, decoupling matrix depression of order and branch road.