CN106293720A - A kind of can graphic programming debugging Power System Simulator - Google Patents

Abstract

The invention discloses a kind of can the Power System Simulator of graphic programming debugging, represent the element storing unit of the internal each assembly graphic element of electrical network including: storage, set and the parameter setting unit of record component graphic element design parameter value and editing interface unit；By transferring the assembly graphic element in element storing unit, combine in described editing interface unit, completed the making of electric power mixing system by parameter setting unit；Also there is parameter setting unit and the processing unit of assignment component graphic element design parameter；During work, the electrical network basic data at least including machine unit characteristic, load information and electric power system tide that described processing unit sets according to parameter setting unit, the electrical network networking prediction scheme obtaining editor carries out simulating, verifying, if checking not up to system requirements, then updates networking plan.

Description

A kind of can graphic programming debugging Power System Simulator

Technical field

The invention belongs to visual programming field, particularly relate to a kind of can graphic programming debugging Power System Simulator.

Background technology

Along with Urbanization in China is constantly accelerated, the regional hub city area of hot zones is constantly increasing, and leads Electricity of applying increases sharply, old network system, effectively allocates model owing to lacking, thus causes generating mixing system to expire The most growing generated energy demand of foot, new dilatation area, the present city of excrescence needs subregion to ration the power supply to solve electric power allotment and need Ask, and these local produce often, logistics type enterprise, ration the power supply and can cause bigger economic loss.

And electrical network upgrading needs to put into great amount of cost, the especially design during electrical network upgrading and checking In the stage, owing to existing program means is single, causing checking amendment long in time limit, cost is high.

Summary of the invention

For the proposition of problem above, the present invention propose a kind of can the Power System Simulator of graphic programming debugging, bag Include:

Storage represents the element storing unit of the internal each assembly graphic element of electrical network, sets and record component graphic element tool The parameter setting unit of body parameter value and editing interface unit；

By transferring the assembly graphic element in element storing unit, combine in described editing interface unit, by ginseng Number setup unit completes the making of electric power mixing system；

Also there is parameter setting unit and the processing unit of assignment component graphic element design parameter；

During work, described processing unit at least includes machine unit characteristic, load information according to what parameter setting unit set With the electrical network basic data of electric power system tide, the electrical network networking prediction scheme obtaining editor carries out simulating, verifying, if checking does not reaches To system requirements, then update networking plan.

As preferred embodiment, described processing unit workflow is as follows: current system trend is included in n node The network loss total amount that system is current is calculated on the basis of power system ground state trendThe unit of node i is meritorious exerts oneself as P_Gi ⁰, negative The meritorious amount of lotus is P_Li ⁰；

Add-on system disturbance: the i.e. generated power of node i is exerted oneself or load meritorious amount sudden change, wherein, generated power Sudden Changing Rate is Δ P_Gr-i, load Sudden Changing Rate of gaining merit is Δ P_Lr-i, now in system, network loss total amount is unknown, utilizes the net that step 1 obtains Harm is current system network loss approximation, and then solves current system imbalance power approximation, then current system imbalance merit Rate P_accFor:

${P_{acc}}^0=(\underset{i=1}{\overset{n}{\Σ}}{P_{Gi}}^2+\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\ΔP}}_{Gr-i})-(\underset{i=1}{\overset{n}{\Σ}}{P_{Li}}^0+\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\ΔP}}_{Lr-i})-P_{\Σ\mathrm{lo}ss}^0$

Understand according to the electromechanical properties of generating set:

$P_{acc}=T_{J\mathrm{\Σ}}\frac{d\mathrm{\ω}}{dt}+(K_{G\mathrm{\Σ}}+K_{L\mathrm{\Σ}})\frac{\mathrm{\Δ}\mathrm{\ω}}{2\mathrm{\π}}$

Then system currently approximates frequency difference Δ ω:

$\mathrm{\Δ}\mathrm{\ω}=\frac{{P_{acc}}^0}{(\frac{T_{J\mathrm{\Σ}}}{\mathrm{\Δ}t}+\frac{K_{G\mathrm{\Σ}}+K_{L\mathrm{\Σ}}}{2\mathrm{\π}})}$

According to generating set electromechanical properties, Primary regulation characteristic and load quiet frequency characteristic, elapsed time Δ t, power train In system, generating set inertia is exerted oneself variation delta P_GTi, Primary regulation is exerted oneself variation delta P_Gi, and load is in quiet frequency characteristic effect Under variation delta P_Li:

${\mathrm{\ΔP}}_{GTi}=T_{Ji}\frac{\mathrm{\Δ}\mathrm{\ω}}{\mathrm{\Δ}t}$

${\mathrm{\ΔP}}_{Gi}=K_{Gi}\frac{\mathrm{\Δ}\mathrm{\ω}}{2\mathrm{\π}}$

${\mathrm{\ΔP}}_{Li}=K_{Li}\frac{\mathrm{\Δ}\mathrm{\ω}}{2\mathrm{\π}}$

The unit of corresponding current power system node i is meritorious exerts oneself as P_Gi ¹, the meritorious amount of load be P_Li ¹As follows:

P_Gi ¹=P_Gi ⁰+ΔP_Gr-i+ΔP_Gi+ΔP_GTi

P_Li ¹=P_Li ⁰+ΔP_Lr-i+ΔP_Li

Wherein: Δ t is simulation time step-length, unit: s, T_JiFor the inertia time constant of system generator group i, unit: s, T_J∑For the inertia time constant sum of all generating sets of system, unit: s, K_GiFor the unit power regulation of generating set i, list Position: MW/Hz, K_GΣFor the unit power regulation sum of all generating sets of system, unit: MW/Hz, K_LiUnit for load is adjusted Joint power, unit: MW/Hz, K_LiFor the unit power regulation sum of all loads of system, unit: MW/Hz；

Utilize that each unit of obtaining is meritorious exerts oneself and the load amount of gaining merit carries out Load flow calculation, obtain new power system net Damage total amountThen relative to the ground state trend before the Δ t time, the total variation of current power system lossesFor:

${\mathrm{\ΔP}}_{\mathrm{\Σ}loss}^k=P_{\mathrm{\Σ}loss}^k-P_{\mathrm{\Σ}loss}^{k-1}(k\≥1),$

Wherein: k is Load flow calculation iterations after disturbance occurs, ground state trend correspondence iterations k=0；

Arranging iteration precision value ε, accuracy value ε is the least, and result of calculation degree of accuracy is the highest；IfIllustrate to pass through After time Δ t, system losses change, when the power system imbalance power approximation that obtains is big with actual phase ratio error, This error approximation and the total variation of power system network lossEqual, then order:

${P_{acc}}^k={\mathrm{\ΔP}}_{\mathrm{\Σ}loss}^k+{\mathrm{\ΔP}}_{GTi}+{\mathrm{\ΔP}}_{Gi}-{\mathrm{\ΔP}}_{Li}$

ΔP_Gr-i=0,

ΔP_Lr-i=0

Load flow calculation iterations k adds 1, returns step 3 and is iterated calculating；

IfIterative computation stops.

Further, loss allocation COEFFICIENT K_{loss_i}Value is the unit power regulation K of each unit_Gi, unit: MW/ Hz, the system that accounts for all units unit power regulation sum Σ K_GiPercentage ratio, it may be assumed that

$K_{loss\_i}=\frac{K_{Gi}}{\underset{i=1}{\overset{N}{\Σ}}{K}_{Gi}}\×100\%(1\≤i\≤N);$

When the regulation allowance that unit i participates in frequency regulation is 0, make K_Gi=0, it may be assumed that the corresponding current loss allocation system of unit i Number K_{loss_i}=0.

Further, nodal information change refers mainly to the generated power of power system interior joint i and exerts oneself change or have Workload changes, and meets:

Power system ground state trend in the method refer to this power system disturbance occur before, whole system supply and demand put down System load flow under weighing apparatus, frequency steady state.

Accompanying drawing explanation

In order to be illustrated more clearly that the embodiment of the present application or technical scheme of the prior art, below will be to embodiment or existing In having technology to describe, the required accompanying drawing used is briefly described, it should be apparent that, the accompanying drawing in describing below is only this Some embodiments described in application, for those of ordinary skill in the art, on the premise of not paying creative work, Other accompanying drawing can also be obtained according to these accompanying drawings.

Fig. 1 is the system flow chart of the present invention

Detailed description of the invention

For making technical scheme and advantage clearer, below in conjunction with the accompanying drawing in the embodiment of the present invention, to this Technical scheme in inventive embodiments carries out the most complete description:

As shown in Figure 1 a kind of can graphic programming debugging Power System Simulator, including:

Storage represents the element storing unit of the internal each assembly graphic element of electrical network, sets and record component graphic element tool The parameter setting unit of body parameter value and editing interface unit；

By transferring the assembly graphic element in element storing unit, combine in described editing interface unit, by ginseng Number setup unit completes the making of electric power mixing system；

Also there is parameter setting unit and the processing unit of assignment component graphic element design parameter；

During work, described processing unit at least includes machine unit characteristic, load information according to what parameter setting unit set With the electrical network basic data of electric power system tide, the electrical network networking prediction scheme obtaining editor carries out simulating, verifying, if checking does not reaches To system requirements, then update networking plan.

As preferred embodiment, described processing unit workflow is as follows: current system trend is included in n node The network loss total amount that system is current is calculated on the basis of power system ground state trendThe unit of node i is meritorious exerts oneself as P_Gi ⁰, negative The meritorious amount of lotus is P_Li ⁰；

Add-on system disturbance: the i.e. generated power of node i is exerted oneself or load meritorious amount sudden change, wherein, generated power Sudden Changing Rate is Δ P_Gr-i, load Sudden Changing Rate of gaining merit is Δ P_Lr-i, now in system, network loss total amount is unknown, utilizes the net that step 1 obtains Harm is current system network loss approximation, and then solves current system imbalance power approximation, then current system imbalance merit Rate P_accFor:

${P_{acc}}^0=(\underset{i=1}{\overset{n}{\Σ}}{P_{Gi}}^0+\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\ΔP}}_{Gr-i})-(\underset{i=1}{\overset{n}{\Σ}}{P_{Li}}^0+\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\ΔP}}_{Lr-i})-P_{\Σ\mathrm{lo}ss}^0$

Understand according to the electromechanical properties of generating set:

$P_{acc}=T_{J\mathrm{\Σ}}\frac{d\mathrm{\ω}}{dt}+(K_{G\mathrm{\Σ}}+K_{L\mathrm{\Σ}})\frac{\mathrm{\Δ}\mathrm{\ω}}{2\mathrm{\π}}$

Then system currently approximates frequency difference Δ ω:

$\mathrm{\Δ}\mathrm{\ω}=\frac{{P_{acc}}^0}{(\frac{T_{J\mathrm{\Σ}}}{\mathrm{\Δ}t}+\frac{K_{G\mathrm{\Σ}}+K_{L\mathrm{\Σ}}}{2\mathrm{\π}})}$

According to generating set electromechanical properties, Primary regulation characteristic and load quiet frequency characteristic, elapsed time Δ t, power train In system, generating set inertia is exerted oneself variation delta P_GTi, Primary regulation is exerted oneself variation delta P_Gi, and load is in quiet frequency characteristic effect Under variation delta P_Li:

${\mathrm{\ΔP}}_{GTi}=T_{Ji}\frac{\mathrm{\Δ}\mathrm{\ω}}{\mathrm{\Δ}t}$

${\mathrm{\ΔP}}_{Gi}=K_{Gi}\frac{\mathrm{\Δ}\mathrm{\ω}}{2\mathrm{\π}}$

${\mathrm{\ΔP}}_{Li}=K_{Li}\frac{\mathrm{\Δ}\mathrm{\ω}}{2\mathrm{\π}}$

The unit of corresponding current power system node i is meritorious exerts oneself as P_Gi ¹, the meritorious amount of load be P_Li ¹As follows:

P_Gi ¹=P_Gi ⁰+ΔP_Gr-i+ΔP_Gi+ΔP_GTi

P_Li ¹=P_Li ⁰+ΔP_Lr-i+ΔP_Li

Wherein: Δ t is simulation time step-length, unit: s, T_JiFor the inertia time constant of system generator group i, unit: s, T_J∑For the inertia time constant sum of all generating sets of system, unit: s, K_GiFor the unit power regulation of generating set i, list Position: MW/Hz, K_GΣFor the unit power regulation sum of all generating sets of system, unit: MW/Hz, K_LiUnit for load is adjusted Joint power, unit: MW/Hz, K_LiFor the unit power regulation sum of all loads of system, unit: MW/Hz；

Utilize that each unit of obtaining is meritorious exerts oneself and the load amount of gaining merit carries out Load flow calculation, obtain new power system net Damage total amountThen relative to the ground state trend before the Δ t time, the total variation of current power system lossesFor:

${\mathrm{\ΔP}}_{\mathrm{\Σ}loss}^k=P_{\mathrm{\Σ}loss}^k-P_{\mathrm{\Σ}loss}^{k-1}(k\≥1),$

Wherein: k is Load flow calculation iterations after disturbance occurs, ground state trend correspondence iterations k=0；

Arranging iteration precision value ε, accuracy value ε is the least, and result of calculation degree of accuracy is the highest；IfIllustrate to pass through After time Δ t, system losses change, when the power system imbalance power approximation that obtains is big with actual phase ratio error, This error approximation and the total variation of power system network lossEqual, then order:

${P_{acc}}^k={\mathrm{\ΔP}}_{\mathrm{\Σ}loss}^k+{\mathrm{\ΔP}}_{GTi}+{\mathrm{\ΔP}}_{Gi}-{\mathrm{\ΔP}}_{Li}$

ΔP_Gr-i=0,

ΔP_Lr-i=0

Load flow calculation iterations k adds 1, returns step 3 and is iterated calculating；

IfIterative computation stops.

Further, loss allocation COEFFICIENT K_{loss_i}Value is the unit power regulation K of each unit_Gi, unit: MW/ Hz, the system that accounts for all units unit power regulation sum Σ K_GiPercentage ratio, it may be assumed that

$K_{loss\_i}=\frac{K_{Gi}}{\underset{i=1}{\overset{N}{\Σ}}{K}_{Gi}}\×100\%(1\≤i\≤N);$

When the regulation allowance that unit i participates in frequency regulation is 0, make K_Gi=0, it may be assumed that the corresponding current loss allocation system of unit i Number K_{loss_i}=0.

Further, nodal information change refers mainly to the generated power of power system interior joint i and exerts oneself change or have Workload changes, and meets:

Power system ground state trend in the method refer to this power system disturbance occur before, whole system supply and demand put down System load flow under weighing apparatus, frequency steady state.

The above, the only present invention preferably detailed description of the invention, but protection scope of the present invention is not limited thereto, Any those familiar with the art in the technical scope that the invention discloses, according to technical scheme and Inventive concept equivalent or change in addition, all should contain within protection scope of the present invention.

Claims (4)

1. one kind can graphic programming debugging Power System Simulator, it is characterised in that including:

Storage represents that the element storing unit of the internal each assembly graphic element of electrical network, setting record component graphic element are specifically joined The parameter setting unit of numerical value and editing interface unit；

By transferring the assembly graphic element in element storing unit, combine in described editing interface unit, set by parameter Cell completes the making of electric power mixing system；

Also there is parameter setting unit and the processing unit of assignment component graphic element design parameter；

During work, described processing unit at least includes machine unit characteristic, load information and electricity according to what parameter setting unit set The electrical network basic data of Force system trend, the electrical network networking prediction scheme obtaining editor carries out simulating, verifying, if checking is not up to System requirement, then update networking plan；

Now, in described processing unit rule of thumb storehouse, networking is detected, to not meeting experience by the empirical data of storage The content of data, provides amending advice；

When being verified, after reaching system requirements, described processing unit stores the prediction scheme reaching requirement, simultaneously to remote service Device end sends this prediction scheme, and then sends to the terminal specified or mobile terminal.

The visual programming system of consideration electrical network dynamic power flow the most according to claim 1, is further characterized in that described Processing unit workflow is as follows: current system trend calculates system on the basis of being included in n node power Ground trend and works as Front network loss total amountThe unit of node i is meritorious exerts oneself as P_Gi ⁰, the meritorious amount of load is P_Li ⁰；

Add-on system disturbance: the i.e. generated power of node i is exerted oneself or load meritorious amount sudden change, and wherein, generated power suddenlys change Amount is Δ P_Gr-i, load Sudden Changing Rate of gaining merit is Δ P_Lr-i, now in system, network loss total amount is unknown, and the network loss utilizing step 1 to obtain is made For current system network loss approximation, and then solve current system imbalance power approximation, then current system imbalance power P_acc For:

${P_{acc}}^0=(\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{P_{Gi}}^0+\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\ΔP}}_{Gr-i})-(\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{P_{Li}}^0+\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\ΔP}}_{Lr-i})-P_{\mathrm{\Σ}loss}^0$

Understand according to the electromechanical properties of generating set:

$P_{acc}=T_{J\mathrm{\Σ}}\frac{d\mathrm{\ω}}{dt}+(K_{G\mathrm{\Σ}}+K_{L\mathrm{\Σ}})\frac{\mathrm{\Δ}\mathrm{\ω}}{2\mathrm{\π}}$

Then system currently approximates frequency difference Δ ω:

$\mathrm{\Δ}\mathrm{\ω}=\frac{{P_{acc}}^0}{(\frac{T_{J\mathrm{\Σ}}}{\mathrm{\Δ}t}+\frac{K_{G\mathrm{\Σ}}+K_{L\mathrm{\Σ}}}{2\mathrm{\π}})}$

According to generating set electromechanical properties, Primary regulation characteristic and load quiet frequency characteristic, elapsed time Δ t, in power system Generating set inertia is exerted oneself variation delta P_GTi, Primary regulation is exerted oneself variation delta P_Gi, and load is under quiet frequency characteristic effect Variation delta P_Li:

${\mathrm{\ΔP}}_{GTi}=T_{Ji}\frac{\mathrm{\Δ}\mathrm{\ω}}{\mathrm{\Δ}t}$

${\mathrm{\ΔP}}_{Gi}=K_{Gi}\frac{\mathrm{\Δ}\mathrm{\ω}}{2\mathrm{\π}}$

${\mathrm{\ΔP}}_{Li}=K_{Li}\frac{\mathrm{\Δ}\mathrm{\ω}}{2\mathrm{\π}}$

The unit of corresponding current power system node i is meritorious exerts oneself as P_Gi ¹, the meritorious amount of load be P_Li ¹As follows:

P_Gi ¹=P_Gi ⁰+ΔP_Gr-i+ΔP_Gi+ΔP_GTi

P_Li ¹=P_Li ⁰+ΔP_Lr-i+ΔP_Li

Wherein: Δ t is simulation time step-length, unit: s, T_JiFor the inertia time constant of system generator group i, unit: s, T_J∑ For the inertia time constant sum of all generating sets of system, unit: s, K_GiFor the unit power regulation of generating set i, unit: MW/Hz, K_GΣFor the unit power regulation sum of all generating sets of system, unit: MW/Hz, K_LiUnit for load regulates merit Rate, unit: MW/Hz, K_LiFor the unit power regulation sum of all loads of system, unit: MW/Hz；

Utilize that each unit of obtaining is meritorious exerts oneself and the load amount of gaining merit carries out Load flow calculation, obtain new power system network loss total AmountThen relative to the ground state trend before the Δ t time, the total variation of current power system lossesFor:

${\mathrm{\ΔP}}_{\mathrm{\Σ}loss}^k=P_{\mathrm{\Σ}loss}^k-P_{\mathrm{\Σ}loss}^{k-1},(k\≥1),$

Wherein: k is Load flow calculation iterations after disturbance occurs, ground state trend correspondence iterations k=0；

Arranging iteration precision value ε, accuracy value ε is the least, and result of calculation degree of accuracy is the highest；IfElapsed time is described After Δ t, system losses change, and when the power system imbalance power approximation that obtains is big with actual phase ratio error, this is by mistake Difference approximation and the total variation of power system network lossEqual, then order:

${P_{\mathrm{acc}}}^k={\mathrm{\ΔP}}_{\mathrm{\Σloss}}^k+{\mathrm{\ΔP}}_{\mathrm{GTi}}+{\mathrm{\ΔP}}_{\mathrm{Gi}}-{\mathrm{\ΔP}}_{\mathrm{Li}}$

ΔP_Gr-i=0,

ΔP_Lr-i=0

Load flow calculation iterations k adds 1, returns step 3 and is iterated calculating；

IfIterative computation stops.

The most according to claim 2 a kind of can graphic programming debugging Power System Simulator, be further characterized in that: network loss Partition coefficient K_{loss_i}Value is the unit power regulation K of each unit_Gi, unit: MW/Hz, the system that accounts for all units unit regulates Power sum Σ K_GiPercentage ratio, it may be assumed that

$K_{loss\_i}=\frac{K_{Gi}}{\underset{i=1}{\overset{N}{\Σ}}{K}_{Gi}}\×100\%,(1\≤i\≤N);$

When the regulation allowance that unit i participates in frequency regulation is 0, make K_Gi=0, it may be assumed that the corresponding current loss allocation coefficient of unit i K_{loss_i}=0.

The most according to claim 3 a kind of can graphic programming debugging Power System Simulator, be further characterized in that:

Nodal information change refers mainly to the generated power of power system interior joint i and exerts oneself change or burden with power change, and full Foot:Or

Power system ground state trend in the method refer to this power system before disturbance occurs, the whole system equilibrium of supply and demand, frequently System load flow under rate steady statue.