CN104765282A - All-digital comprehensive simulation platform of power station - Google Patents

Abstract

An all-digital comprehensive simulation platform of a power station achieves graphical modeling of a power station power grid part, a flow net system and a control system. The established simulation models can perform independent simulation calculation of subsystems or devices or the like, united closed-loop computation of all system models of the whole power station, all-digital real-time simulation and simulation analysis experiments. The platform comprises a simulator module, real-time calculated data can be automatically linked with simulator virtual power station monitoring images or equipment variations in virtual equipment or physical subsystems, and then functions of the real-time digital simulator can be achieved. The platform directly calls the simulation models in analysis and computation, and conveniently meets simulation demands of a power station transient process, a middle-and-long term process, a flow network, the control system and the like.

Description

The digital Synthetic Simulation Platform in a kind of power station

Technical field

The digital Synthetic Simulation Platform in a kind of power station of the present invention, belongs to power system digital simulation field.

Background technology

Existing real time data display generally adopts ripe electromechanical transient simulation program or electromagnetic transients program, real-timedigital simulation lays particular emphasis on the interface realized with physical system, thus realization Digital Simulation System is to the test and check of physical system, and adopt special hardware platform, expensive, such as RTDS, ADPSS etc., the outstanding accuracy of simulation result and the real-time of calculating.For ADPSS, employing be the electromechanical transient simulation program of PSASP maturation be core calculations program; For RTDS, employing be the electromagnetic transients program of EMTD maturation be core calculations program.Existing electric system realtime digital simulation platform can carry out closed loop calculating by external inspecting device material object, thus the object that realization emulation electric system is tested supervising device in kind, there is very important using value.For complete electric power station system, include the systems such as electrical network, auxiliary system, prime mover energy conversion, control, if desired the research and learning is carried out to this holonomic system, need to realize the system modeling and simulations interdisciplinary such as electrical network, drift net, control, existing realtime digital simulation platform does not have such function.

The needs of transient state process research when existing Power System Stability Analysis program is to experience disturbance under realizing certain original state of electric system, such as PSASP, EMTP etc. are by emulating the instrument of analyzing and researching to certain concrete operating mode of electric system, the accuracy of outstanding simulation result.But some running status and initial operating mode thereof are difficult to obtain sometimes in electric system.Real-time calculation procedure, compared with the dynamic behaviour of the easier observing system of analytical calculation program, can provide initial operating mode for the analytical calculation of physical system dynamic perfromance.

Existing emulation tool reaches the object that professional studies electric system to a great extent, such as PSASP, EMTP etc., adopts strict mathematical model when their common feature is modeling.But the professional of most beginner or just intermediate level, such as many operations staffs, the student of power specialty, they lack one and adopt the model that represent real complete electric system of strict mathematical model modeling to carry out studying, learn and testing, also the proper object that all kinds of actual power device (as hydro-generating Unit, aut.eq. etc.) is studied, learns and tested is lacked, real-time simulation instrument may solve the problem, but 1): RTDS, ADPSS adopts special hardware platform, expensive; 2): the simulation problems that these emulation tools can only solve single subject usually does not possess the copying etc. of other system as systems such as drift nets; 3) needs that users carry out imparting knowledge to students, test and learning: the professional softwares such as RTDS, ADPSS lay particular emphasis on the experiment to actual physics device, are not positioned.Current emulation tool can only solve the simulation problems of single subject usually; electric power station system is an Iarge-scale system comprising drift net, electrical network, protection and control, and people are difficult to find suitable emulation tool to set up this complete perfect power station model and carry out the research and learning.

Existing power station simulation training system contains electrical network, auxiliary system, energy conversion system, Control protection system etc.; therefore more as a real virtual electric power system; but the general employing mathematical model of non-critical or the strict calculating etc. of logical model or non-closed loop; be difficult to that there is the meaning of research or the effect of help comprehension theory, be usually suitable only for the training of equipment operating.Existingly similar comprise electrical network, the replicating machine platform of drift net is developed as object mainly with replicating machine, what the mathematics library of grid parts generally adopted is power flow algorithm, and electrical network calculates, genset calculates, Protection control system calculates the strict closed loop being difficult to be associated as power station electric part model and calculates; The mathematical model of what the mathematics library of drift net generally adopted is single-phase flow or homogeneous flow.Association during whole system combined calculation between subsystem computation model calculates not strict association, closed loop calculates, some subsystem realistic model even logic-based model etc. and mathematical model of nonphysical system.Traditional electric system simulation machine feature: what Simulation Calculation adopted is non-strict mathematical model or logical model, thus the data deficiency Research Significance emulated or the reference significance to production, rudimentary training in being applicable to carrying out, Simulation Calculation adopts high-speed servers to calculate usually, training terminal machine is the operating terminal of Data receiver and student, even if the form of one machine with multi-mode, also be difficult to change multiple student the uses same simulation training system fact when exercise operates, training efficiency is low, and the startup of server time is long, cause that again to recover the initial operating mode time long, affect training efficiency.

Summary of the invention

The digital Synthetic Simulation Platform in a kind of power station of the present invention, adopts on the basis of strict mathematical model, forms complete, perfect power station simulation model by mathematics library mode.General commercial PC can be adopted as operation platform, realize the all-digital real-time simulation of power station total system overall process; Not only meet the needs that the operation and control of user to power station holonomic system and each subsystem, subsystem is studied, learns and tested, also can meet the needs that user studies Blast Furnace Top Gas Recovery Turbine Unit (TRT) various types of in system, aut.eq. etc., learn and tests.Real-time simulation data can the equipment variables of automatic corresponding monitoring power station picture, user is by the screen operation the same with the supervisory system of reality, the operation monitoring whole virtual power station, reach the object of training, teaching, examination, realize the function of power station real time simulator, owing to adopting strict mathematical model, also solve the problem that traditional simulation machine is generally only suitable for just middle rank training.

The digital Synthetic Simulation Platform in a kind of power station, comprise power station main system mathematics library, control system mathematics library, power station main system mathematics library is used for realizing primary electrical system, the mathematics library of fluid network, form primary electrical system realistic model respectively, auxiliary system realistic model, energy conversion system realistic model, in calculating in real time or analytical calculation, realize power station total system closed loop calculate, control system mathematics library presses motor, excitation system, prime mover and governing system thereof, aut.eq., protective device, supervising device, the physics subsystem that commutation inversion device etc. are actual or device form realistic model, this realistic model by the main system element calling of being correlated with or can be enabled before power station main system forms realistic model, participate in the total system closed loop calculating of power station, picture at different levels during power station total system real-time closed-loop calculates, as: each main system, the graphic picture of physics subsystem or device all can be opened and monitor,

Comprise real-time computation and analysis computing function, analytical calculation realizes by calling each realistic model, its initial operating mode can set in mathematics library process, or import the profile data preserved by the realistic model such as each subsystem or subsystem in calculating in real time, primary electrical system calculates in real time or analytical calculation realizes calculating or the debugging of primary electrical system and Control System Imitation model thereof, auxiliary system calculates in real time or analytical calculation realizes calculating or the debugging of auxiliary system and Control System Imitation model thereof, energy conversion system calculates in real time or analytical calculation realizes calculating or the debugging of energy conversion system and Control System Imitation model thereof, control system calculates in real time or analytical calculation realizes calculating or the debugging of Control System Imitation model, power station total system calculates in real time or analytical calculation realizes the associating closed loop calculating of power station total system realistic model or debugs,

Comprise replicating machine module, replicating machine module can control the real-time computing function of electric power station system realistic model, real-time calculating data can the equipment variables automatically with replicating machine in virtual monitoring power station picture or virtual unit or physics subsystem link, and realize real time simulator function.

Described total system closed loop calculates: comprise the iterative computation between primary electrical system realistic model, auxiliary system realistic model, energy conversion system realistic model and calculate with the closed loop of combining of related Control System realistic model.

Described power station main system mathematics library, be the realistic model of transient stability realistic model in order to form electric system and drift net system such as auxiliary system, prime mover side energy conversion system, drift net system simulation model is processed into the form consistent with network system realistic model, composed component topological relation after graphic elements connects, except node graphic elements, the element in physical system is all represented in component library, bus or node form electrical network/drift net node jointly, other elements can only be connected with bus or node, the type of node or bus need be indicated, as the PQ node in electrical network, PV node, SLACK node, drift net similar with it, initialization for realistic model calculates, component parameters is arranged all can on electrical network/flow-net diagram, graph of a correspondence element be revised or preserve, component parameters not only comprises the parameter of element itself, also the mark associated of this graphic elements and certain element on other figure or other figure is comprised, the selection of component models type, fault verification etc., real-time calculating data on figure at different levels, graph of a correspondence element can be monitored.

Described control system mathematics library, it is the realistic model in order to realize the actual physics subsystem such as genset, motor, excitation system, prime mover and governing system thereof, aut.eq., protective device, supervising device, commutation inversion device or device, calculate in real time in control system or realize calculating or the debugging of Control System Imitation model in analytical calculation, the realistic model formed can use in the modeling of supply station main system; Graphic elements in component library is various general purpose function modules, according to certain system transter block diagram, rear composed component topological relation is connected to graphic elements, component parameters is arranged all can on control system figure, graph of a correspondence element be revised or preserve, component parameters not only comprises the parameter of element itself, also comprise the mark that associate of this graphic elements with certain element on other figure or other figure, the selection, fault verification etc. of component models type, control system calculates data in real time and can monitor on the corresponding graphic elements of this figure; Control system mathematics library provides open component module; user also also passes through the realistic model of user oneself programming realization element by the parameter of the component module definable element opened by the topological relation of the required system of control system mathematics library acquisition, and carries out debugging, analyze or calling under control system mathematics library environment.

The digital Synthetic Simulation Platform in a kind of power station of the present invention, technique effect is as follows:

1), the modeling method of complete, the sophisticated model of power station total system and mathematics library platform, general commercial PC can be adopted as operation platform, realize the all-digital real-time simulation of power station total system overall process, not only meet the needs that the operation and control of user to power station holonomic system is studied, learns and tested, also can meet the needs that user studies Blast Furnace Top Gas Recovery Turbine Unit (TRT) various types of in system, aut.eq. or subsystem etc., learn and tests;

2) simulation analysis function and all-digital real-time simulation function, is included, profile data in real-time calculating is stored in realistic model by each subsystem or subsystem, in analytical calculation, directly call realistic model, realize the needs that power station transient state process, medium-term and long-term process, fluid network, control system etc. are emulated easily.

3), real-time calculating data can the equipment variables automatically with replicating machine in virtual monitoring power station picture or virtual unit or physics subsystem link, realize real time simulator function, user is by the virtual monitoring power station screen operation consistent with actual power station or equipment, monitor the operation in whole virtual power station, reach training, teaching, the object of examination, realize the function of the real time simulator of power station or equipment or physics subsystem, real time simulator is owing to adopting strict mathematical model, therefore the problem that traditional simulation machine is generally only suitable for training on operation is solved.

Accompanying drawing explanation

Fig. 1 is computation structure figure of the present invention.

Fig. 2 is functional module relation schematic diagram of the present invention.

Fig. 3 is total system closed loop Computing Principle block diagram of the present invention.

Fig. 4 is the realistic model figure of fluid network branch road m of the present invention.

Fig. 5 is fluid network calculation flow chart of the present invention.

Fig. 6 is that electric part of the present invention replaces derivation algorithm process flow diagram.

Embodiment

The digital Synthetic Simulation Platform in a kind of power station, comprises main system mathematics library, control system mathematics library.Main system mathematics library is completed by primary electrical system, fluid network mathematics library, form primary electrical system realistic model, auxiliary system realistic model, energy conversion system realistic model respectively, in calculating in real time or analytical calculation, realize total system associating closed loop calculate, Fig. 3 is shown in by total system closed loop Computing Principle block diagram.

Power station primary electrical system comprises the middle-size and small-size electric system of genset, station service electrical system, minitype high voltage electrical network, equivalent electric system formation, realistic model is made up of power system simulation model and genset realistic model, adopts the calculating of the Pattern completion primary electrical system alternately solved.Genset realistic model is made up of synchronous motor realistic model, excitation system realistic model, prime mover and governing system realistic model thereof.The each physics subsystem realistic model of genset can be generated by control system mathematics library, and called in primary electrical system mathematics library.Primary electrical system realistic model can calculate in real time at this platform primary electrical system or realize in analytical calculation calculating and debugging.Fig. 6 is shown in by power station primary electrical system realistic model Computing Principle block diagram.

Power plant auxiliary machine system comprises wind, water, oily three digest journals, the realistic model of power plant auxiliary machine system can be obtained by the mathematics library of this platform flow volume grid, Fig. 5 is shown in by the Computing Principle block diagram of auxiliary system realistic model, can this platform auxiliary system calculate in real time or analytical calculation realize calculate and debugging.Auxiliary system and monitoring/protection system simulation model thereof participate in total system simulation calculation.Theory diagram in total system simulation calculation is shown in Fig. 3.

Power station energy conversion system is for the therrmodynamic system of thermal power plant, the realistic model of Power Plant Thermal System can be obtained by the mathematics library of this platform flow volume grid, Fig. 5 is shown in by the Computing Principle block diagram of Module in Thermodynamic System Simulation model, can this platform energy conversion system calculate in real time or analytical calculation realize calculate and debugging.Therrmodynamic system and monitoring/protection system simulation model thereof participate in total system simulation calculation.Theory diagram in total system simulation calculation is shown in Fig. 3.

Power station control system mathematics library comprises the realistic model of the actual physics subsystem such as synchronous motor, excitation system, prime mover and governing system thereof, aut.eq., protective device, supervising device, commutation inversion device or device etc., can to calculate in real time or analytical calculation realizes calculating and debugging at this platform control system.These realistic models to be formed before realistic model by the main system element calling of being correlated with in power station main system or enable, and participate in during power station total system closed loop calculates, the theory diagram in total system simulation calculation is shown in Fig. 3.

Power station total system realistic model and each module calculated relationship, the as shown in Figure 1 digital Synthetic Simulation Platform total system in power station computation structure figure.

The formation of power station total system realistic model and the relation of functional module, as shown in Figure 2 power station all-digital real-time simulation system functional module relation schematic diagram.Power station total system realistic model comprises power station primary electrical system realistic model, complete comprise genset, the transient stability simulation calculation of middle-size and small-size electric system that station service electrical system, minitype high voltage electrical network, equivalent electric system are formed; Auxiliary system realistic model, completes the simulation calculation of power station wind, water, oil system; Energy conversion system realistic model, completes the simulation calculation etc. of power station prime mover side therrmodynamic system or Hydraulic Power System.The functional module of power station full digital trigger technique platform comprises: operating mode section imports, generate source file, electrical network calculates in real time, main system calculates in real time, subsidiary engine/energy conversion system calculates in real time, control system calculates in real time, analytical calculation, curve module, the preservation of operating mode section, replicating machine module, monitoring/protection picture module etc.As shown in Figure 2, real-time analytic method calculates and realizes by calling relevant realistic model; Main system calculates the combined calculation or debugging that realize main system and Control System Imitation model thereof in real time; Primary electrical system calculates the combined calculation or debugging that realize primary electrical system and Control System Imitation model thereof in real time; Subsidiary engine/energy conversion system calculates the combined calculation or debugging that realize subsidiary engine/energy conversion system and Control System Imitation model thereof in real time.Realize by calling corresponding realistic model the analytical calculation of each system, the profile data in calculating in real time is stored in each subsystem or subsystem realistic model, and corresponding profile data can call as initial operating mode in analytical calculation.

The graphic elements of mathematics library can make by C language, or makes graphic elements with existing graphics development software.The effect of graphic elements is the different physical component in representative system, and graphic elements can identify and connects and form the topological relation determined upon connection, can eject parameter interface to go forward side by side the preservation of line parameter or display by graphic elements.Graphical components library comprises power system component storehouse, compressible fluid component library, incompressible fluid component library, control system component storehouse etc.

Primary electrical system mathematics library, except node graphic elements, all represent the element in physical system in component library, bus or node form grid nodes jointly, and other elements can only be connected with bus or node; Component parameters is arranged all can on electrical network figure, graph of a correspondence element be revised or preserve, component parameters not only comprises the parameter of element itself, also the correlation parameter of the supervising device relevant with element and the selection, fault verification etc. of component models type is comprised, Load flow calculation data all on electrical network figure, graph of a correspondence element can show, and calculate data in real time and all on electrical network figure, graph of a correspondence element can be monitored.The type of node or bus need be indicated, as the PQ node in electrical network, PV node, SLACK node, drift net similar with it, the initialization for realistic model calculates.

Primary electrical system mathematical model comprises electric network model, genset model, wherein genset model is by generator model, excitation system model, prime mover and governing system model-composing, the realistic model that each element mathematical model is formed can stored in module library, dissimilar realistic model can be called in model bank in the initialization procedure of grid graph modeling, such as synchronous motor can select the utility model etc. such as 3 rank, 4 rank, 5 rank be applicable in networking operation, also can select the Synchronous Machine Models etc. be applicable in start-up course.

Power station electric part transient stability Digital Simulation computation model algorithm flow chart is shown in Fig. 6.Power station fluid network calculation flow chart is shown in Fig. 5.

Power station electric part Simulation Calculation is formed by adopting the power system simulation model that alternately solves and genset realistic model.The core of electric power system transient stability simulation calculation is simultaneous solution or alternately solves network equation group and genset system of equations.

One, power system simulation model:

Nodal voltage equation: $Y\stackrel{\·}{U}=\stackrel{\·}{I}---\left(1\right)$

Wherein Y is admittance matrix, for node voltage vector, for node Injection Current vector.

Generator unit stator winding voltage equation is made to be:

$\left[\begin{array}{c}{U}_d\\ U_q\end{array}\right]=\left[\begin{array}{c}{E}_d^{\′}\\ E_q^{\′}\end{array}\right]-\left[\begin{array}{cc}{r}_a& -X_q^{\′}\\ X_d^{\′}& r_a\end{array}\right]\left[\begin{array}{c}{I}_d\\ I_q\end{array}\right]---\left(2\right)$

As X ' _d≠ X ' _qtime, order $\left[\begin{array}{cc}{r}_a& -X_q^{\′}\\ X_d^{\′}& r_a\end{array}\right]=\left[\begin{array}{cc}{r}_a& -X_d^{\′}\\ X_d^{\′}& r_a\end{array}\right]+\left[\begin{array}{cc}0& X_d^{\′}-X_q^{\′}\\ 0& 0\end{array}\right],$

Voltage equation then on generator unit stator winding dq axle can arrange and be:

$\left[\begin{array}{c}{U}_d\\ U_q\end{array}\right]=\left[\begin{array}{c}{E}_d^{\′}\\ E_q^{\′}\end{array}\right]-\left[\begin{array}{cc}{r}_a& -X_d^{\′}\\ X_d^{\′}& r_a\end{array}\right]\left[\begin{array}{c}{I}_d\\ I_q\end{array}\right]+\left[\begin{array}{cc}0& X_d^{\′}-X_q^{\′}\\ 0& 0\end{array}\right]\left[\begin{array}{c}{I}_d\\ I_q\end{array}\right]---\left(3\right)$

Order, E ' _dx=E ' _d-(X ' _d-X ' _q) I _q(4)

Then (3) formula turns to:

$\left[\begin{array}{c}{U}_d\\ U_q\end{array}\right]=\left[\begin{array}{c}{E}_{\mathrm{dx}}^{\′}\\ E_q^{\′}\end{array}\right]-\left[\begin{array}{cc}{r}_a& -X_d^{\′}\\ X_d^{\′}& r_a\end{array}\right]\left[\begin{array}{c}{I}_d\\ I_q\end{array}\right]---\left(5\right)$

Obviously, (5) formula can be write as the form of vector, that is:

Here, can be write as vector form and nodal voltage equation after (5) formula carries out dq conversion simultaneous solution, and be expressed as following form:

$Y_E\stackrel{\·}{E}=\stackrel{\·}{I}---\left(6\right)$

Here for the virtual potential after the virtual dq axle reactance that the reactance of d, q axle is equal, matrix of coefficients Y _efor permanent matrix, in matrix, the value of each element can not change along with the change of δ.(6) namely formula also becomes power system simulation model.

Two, genset realistic model:

Here implicit trapezoid method is adopted to be solved to example with the differential equation, generator adopts 4 rank utility models, excitation system as Suo Shi (11) formula, prime mover and governing system thereof are depicted as example as (14) formula, and the calculating of genset model is described.

For without loss of generality, (7) to (17) formula is made to represent the calculating of i-th genset model.The compute classes of other genset models seemingly.

Generator amature winding mathematical model:

$\left\{\begin{array}{c}{T}_{d0}^{\′}p{E}_q^{\′}=E_f-E_q^{\′}-(X_d-X_d^{\′})i_d\\ T_{q0}^{\′}p{E}_d^{\′}=-E_d^{\′}+(X_q-X_q^{\′})i_q\end{array}\right.---\left(7\right)$

Adopt hiding-trapezium integral method differencing:

$\begin{array}{c}\left[\begin{array}{cc}{T}_{d0}^{\′}+\mathrm{Dt}& 0\\ 0& T_{q0}^{\′}+\mathrm{Dt}\end{array}\right]\left[\begin{array}{c}{E}_q^{\′}(t+\mathrm{\Δt})\\ E_d^{\′}(t+\mathrm{\Δt})\end{array}\right]\\ =\left[\begin{array}{cc}{T}_{d0}^{\′}+\mathrm{Dt}& 0\\ 0& T_{q0}^{\′}+\mathrm{Dt}\end{array}\right]\left[\begin{array}{c}{E}_q^{\′}\left(t\right)\\ E_d^{\′}\left(t\right)\end{array}\right]+\left[\begin{array}{cc}2\mathrm{Dt}(X_d^{\′}-X_q)& 0\\ 0& 2\mathrm{Dt}(X_q-X_q^{\′})\end{array}\right]\left[\begin{array}{c}{I}_d\left(t\right)\\ I_q\left(t\right)\end{array}\right]+\left[\begin{array}{c}2\mathrm{Dt}\\ 0\end{array}\right]E_f\left(t\right)\end{array}---\left(8\right)$

Above formula turns to:

$\begin{array}{c}\left[\begin{array}{c}{E}_d^{\′}(t+\mathrm{\Δt})\\ E_q^{\′}(t+\mathrm{\Δt})\end{array}\right]\\ ={\left[\begin{array}{cc}2T_{q0}^{\′}+\mathrm{Dt}& 0\\ 0& 2T_{d0}^{\′}+\mathrm{Dt}\end{array}\right]}^{-1}\left[\begin{array}{cc}2T_{q0}^{\′}-\mathrm{Dt}& 0\\ 0& 2T_{d0}^{\′}-\mathrm{Dt}\end{array}\right]\left[\begin{array}{c}{E}_d^{\′}\left(t\right)\\ E_q^{\′}\left(t\right)\end{array}\right]\\ +{\left[\begin{array}{cc}2T_{q0}^{\′}+\mathrm{Dt}& 0\\ 0& 2T_{d0}^{\′}+\mathrm{Dt}\end{array}\right]}^{-1}\left[\begin{array}{cc}0& 2\mathrm{Dt}(X_q-X_q^{\′})\\ 2\mathrm{Dt}(X_d^{\′}-X_d)& 0\end{array}\right]\left[\begin{array}{c}{I}_d\left(t\right)\\ I_q\left(t\right)\end{array}\right]\\ +{\left[\begin{array}{cc}2T_{q0}^{\′}+\mathrm{Dt}& 0\\ 0& 2T_{d0}^{\′}+\mathrm{Dt}\end{array}\right]}^{-1}\left[\begin{array}{c}0\\ 2\mathrm{Dt}\end{array}\right]E_f\left(t\right)\end{array}---\left(9\right)$

Here, perform (4) formula and (5) formula, then because (5) formula is vector form, thus (5) formula after the conversion of dq axle with nodal voltage equation simultaneous solution, obtains (6) formula.

In like manner form excitation system realistic model and obtain E _f(t+ △ t), and prime mover and governing system realistic model thereof obtain δ (t+ △ t).

Three, alternately calculating is solved:

Power system simulation model, genset realistic model have alternately solved the calculating of every step for the moment; Synchronous Machine Models type, excitation system types of models, governing system types of models are not limited to above special case, can select dissimilar model composition grid simulation module as required, genset emulation module alternately solves calculating.

Fig. 6 is shown in by power station electric part Computing Principle block diagram.

Drift net system mathematics library, composed component topological relation after graphic elements connects, all represents the element in physical system except node graphic elements in component library, node forms drift net node, and other elements can only be connected with node; Component library is divided into compressible fluid component library and incompressible fluid component library; Component parameters is arranged all can on flow-net diagram, graph of a correspondence element be modified or preserve, component parameters not only comprises the parameter of element itself, also comprise the correlation parameter of the supervising device relevant with element and the selection, fault verification etc. of component models type, calculate data in real time and all on flow-net diagram, graph of a correspondence element can be monitored.

The formation of drift net realistic model:

1), each for drift net element transient state equational numerical solution equation arranged or be processed into following form:

G _m(t+ △ t)=P _m(t+ △ t) * BX _m(t)+I _m(t), wherein, G _m(t+ △ t) flow for walking lower a period of time; P _m(t+ △ t) is for walking branch road two ends pressure reduction lower a period of time, BX _mt () is the equivalent conductance of branch road m in the step calculating of this time, I _mt () is the traffic sources of branch road m in the step calculating of this time;

2) the drift net node pressure equation, formed according to above formula has following form: GY*P=IM+H ρ+WG, this formula meets the principle of mass conservation, principle of conservation of momentum, conservation of energy principle, wherein P is node pressure vector, GY is node equivalent conductance matrix, and IM+H ρ+WG is node implantation quality flow; Conductance matrix GY is square formation, its exponent number just equals the nodes n of network, conductance matrix is sparse matrix, its each row non-zero nondiagonal element number just equals the series arm number that node corresponding with this row connects, the diagonal element GY (i, i) of conductance matrix just equals the summation of all serial or parallel connection branch road equivalence conductances that this node is connected, the nondiagonal element GY (i of conductance matrix, j) connected node i is equaled, the negative value of the branch road equivalence conductance of j; The mass rate that IM (i) injects node i for traffic sources in each branch road numerical solution of being connected with node i with, WG (i) for be parallel to load equipment in respective nodes or as pressure source equipment as press in the pressure head relevant with voltage, frequency of pump or gas bag the mass rate of injection node i that formed, the mass rate of the injection node i that H ρ (i) is formed for the difference of pressing in exert pressure outside each branch node of being connected with node i strong and branch road and, the density of working medium in this mass rate change pipe network branch road.

Bernoulli equation for the unsteady flo w viscous fluid of branch road m:

${\mathrm{\ρ}}_m{L}_m\frac{{\mathrm{dU}}_m}{\mathrm{dt}}=P_m+H_m^{\′}-h_w---\left(11\right)$

Wherein: H ' _m=ρ _mgZ _m(12)

$h_w=\mathrm{\Σ}(1/2\mathrm{\λ}\frac{L_m}{d_m}+\mathrm{\ξ}){\mathrm{\ρ}}_m{U}_m\left|U_m\right|---\left(13\right)$

And suppose that in branch road, the flow velocity of each point is equal.

By mass rate G _m=ρ _ma _mu _mwith branch road inertia and (13) substitution (11) formula obtains:

$I_m\frac{{\mathrm{dG}}_m}{\mathrm{dt}}=P_m+{\mathrm{\ρ}}_{m}g{Z}_m-R_m{G}_m{G}_m---\left(14\right)$

Wherein:

$R_m=\frac{R_c}{{\mathrm{\ρ}}_m}---\left(15\right)$

$\frac{R_c}{{\mathrm{\ρ}}_m}=\mathrm{\Σ}(1/2\mathrm{\λ}\frac{L_m}{d_m}+\mathrm{\ξ})/A_m^2---\left(16\right)$

Above and in later formula:

I, j are the endpoint number of branch road m, and m is branch road number.

ρ _m: working medium density in branch road m.

H _pump: lift of pump.

L _m: the length of branch road m.

A _m: the sectional area of branch road m.

R _m: the resistance coefficient (relevant with ρ) of branch road m.

H _w: the flow losses that in branch road m, working medium is total.

D _m: the equivalent diameter of branch road m.

Z _m: branch road m first and last end difference of elevation.

I _m: branch road inertia.

λ: frictional resistant coefficient.

ξ: coefficient of shock resistance.

(4) formula obtains after adopting hiding-trapezium integral method differencing:

$\begin{array}{c}{I}_m\frac{G_m(t+\mathrm{\Δt})-G_m\left(t\right)}{\mathrm{\Δt}}\\ =0.5*(P_m(t+\mathrm{\Δt})+P_m\left(t\right))-0.25*R_m*{(G_m(t+\mathrm{\Δt})+G_m\left(t\right))}^2+{\mathrm{\ρ}}_{m}g{Z}_m\end{array}---\left(17\right)$

Above formula can obtain after arranging:

0.25R _m△tG _m(t+△t)G _m(t+△t)+(I _m+0.5R _m△tG _m(t))G _m(t+△t)

＝0.5*P _m(t+△t)+0.5*P _m(t)-0.25R _m△tG _m ²(t)+I _mG _m(t)+ρ _mgZ _m△t (18)

Each for drift net element transient state equational numerical solution equation is arranged or is processed into following form:

G _m(t+△t)＝P _m(t+△t)*BX _m(t)+IM _m(t) (19)

Wherein, G _m(t+ △ t) flow for walking lower a period of time; P _m(t+ △ t) is for walking branch road two ends pressure reduction lower a period of time, BX _mt () is the equivalent conductance of branch road m in the step calculating of this time, IM _mt () is the traffic sources of branch road m in the step calculating of this time;

Here for Forecasting Methodology as follows, (14) formula is variable to be changed to:

DG _m=(P _m+ ρ _mgZ _m-R _mg _mg _m) dt/I _m, due to, G _m(t+ △ t)=G _m(t)+dG _mt (), can make,

G _m(t+△t)＝G _m(t)+(P _m(t)+ρ _mgZ _m-R _m*G _m(t)*G _m(t))△t/I _m

(18) formula can be expressed as following form:

(RX _m(t))G _m(t+△t)

＝P _m(t+△t)+P _m(t)-0.5R _m△tG _m ²(t)+2I _mG _m(t)+2ρ _mgZ _m△t (20)

Wherein,

RX _m(t)＝0.25R _m△t((P _m(t)-R _m*G _m(t)+ρ _mgZ _m)△t/I _m+G _m(t))+I _m+0.5R _m△tG _m(t)

Be organized into the form as shown in (19) further:

Wherein, BX _m(t)=1/RX _m(t) (21)

IM _m(t)＝(P _m(t)-0.5R _m△tG _m ²(t)+2I _mG _m(t)+2ρ _mgZ _m△t)/RX _m(t) (22)

This formula can be expressed as the realistic model figure shown in Fig. 4.IM _mt () represents the injection traffic sources of t branch road m, BX _mt () represents the equivalent conductance of t branch road m.

Fluid network calculation flow chart as shown in Figure 5.

The graphic elements of control system mathematics library can make by C language equally, or makes graphic elements with existing graphics development software.The representative of graphic elements be general purpose function block diagram in system, graphic elements can identify the topological relation connecting and also form the determination between each function block diagram upon connection, parameter interface can be ejected by graphic elements to go forward side by side the preservation of line parameter or display, the simulation algorithms such as hiding-trapezium integral method can be adopted equally by all equation difference differentiation in a width figure, adopt the mode of simultaneous solution to solve this DIFFERENCE EQUATIONS.Control system mathematics library, it is the realistic model in order to realize the actual physics subsystem such as genset, motor, excitation system, prime mover and governing system thereof, aut.eq., protective device, supervising device, commutation inversion device or device, calculate in real time in control system or realize calculating or the debugging of Control System Imitation model in analytical calculation, the realistic model formed can use in the modeling of supply station main system; Graphic elements in component library is various general purpose function modules, according to certain system transter block diagram, rear composed component topological relation is connected to graphic elements, component parameters is arranged all can on control system figure, graph of a correspondence element be revised or preserve, component parameters not only comprises the parameter of element itself, also comprise the mark that associate of this graphic elements with certain element on other figure or other figure, the selection, fault verification etc. of component models type, control system calculates data in real time and can monitor on the corresponding graphic elements of this figure; Control system mathematics library provides open component module; user also also passes through the realistic model of user oneself programming realization element by the parameter of the component module definable element opened by the topological relation of the required system of control system mathematics library acquisition, and carries out debugging, analyze or calling under control system mathematics library environment.

Replicating machine module comprises the management of monitoring power station picture, and realize calling of various kinds of equipment monitored picture by monitoring power station image menu or catalogue, real time data links with the equipment variables in monitoring power station picture automatically; For reducing the calculating scale of main system, the calculating of power station main system and control/protection system thereof does not comprise the calculating of disconnector, and the effect of disconnector in main system is realized by the disconnector module in replicating machine.The major function of disconnector module is the lock operation of disconnector and relative breaker.

Claims (4)

1. the digital Synthetic Simulation Platform in power station, comprise power station main system mathematics library, control system mathematics library, it is characterized in that, power station main system mathematics library is used for realizing primary electrical system, the mathematics library of fluid network, form primary electrical system realistic model respectively, auxiliary system realistic model, energy conversion system realistic model, in calculating in real time or analytical calculation, realize power station total system closed loop calculate, control system mathematics library presses motor, excitation system, prime mover and governing system thereof, aut.eq., protective device, supervising device, the physics subsystem that commutation inversion device etc. are actual or device form realistic model, this realistic model by the main system element calling of being correlated with or can be enabled before power station main system forms realistic model, participate in the total system closed loop calculating of power station, picture at different levels during power station total system real-time closed-loop calculates, as: each main system, the graphic picture of physics subsystem or device all can be opened and monitor,

Comprise real-time computation and analysis computing function, analytical calculation realizes by calling each realistic model, its initial operating mode can set in mathematics library process, or import the profile data preserved by the realistic model such as each subsystem or subsystem in calculating in real time, primary electrical system calculates in real time or analytical calculation realizes calculating or the debugging of primary electrical system and Control System Imitation model thereof, auxiliary system calculates in real time or analytical calculation realizes calculating or the debugging of auxiliary system and Control System Imitation model thereof, energy conversion system calculates in real time or analytical calculation realizes calculating or the debugging of energy conversion system and Control System Imitation model thereof, control system calculates in real time or analytical calculation realizes calculating or the debugging of Control System Imitation model, power station total system calculates in real time or analytical calculation realizes the associating closed loop calculating of power station total system realistic model or debugs,

Comprise replicating machine module, replicating machine module can control the real-time computing function of electric power station system realistic model, real-time calculating data can the equipment variables automatically with replicating machine in virtual monitoring power station picture or virtual unit or physics subsystem link, and realize real time simulator function.

2. the digital Synthetic Simulation Platform in a kind of power station according to claim 1, it is characterized in that, described total system closed loop calculates: comprise the iterative computation between primary electrical system realistic model, auxiliary system realistic model, energy conversion system realistic model and calculate with the closed loop of combining of related Control System realistic model.

3. the digital Synthetic Simulation Platform in a kind of power station according to claim 1, it is characterized in that, described power station main system mathematics library, be the realistic model of transient stability realistic model in order to form electric system and drift net system such as auxiliary system, prime mover side energy conversion system, drift net system simulation model is processed into the form consistent with network system realistic model, composed component topological relation after graphic elements connects, except node graphic elements, the element in physical system is all represented in component library, bus or node form electrical network/drift net node jointly, other elements can only be connected with bus or node, the type of node or bus need be indicated, as the PQ node in electrical network, PV node, SLACK node, drift net similar with it, initialization for realistic model calculates, component parameters is arranged all can on electrical network/flow-net diagram, graph of a correspondence element be revised or preserve, component parameters not only comprises the parameter of element itself, also the mark associated of this graphic elements and certain element on other figure or other figure is comprised, the selection of component models type, fault verification etc., real-time calculating data on figure at different levels, graph of a correspondence element can be monitored.

4. the digital Synthetic Simulation Platform in a kind of power station according to claim 1, it is characterized in that, described control system mathematics library, it is the realistic model in order to realize the actual physics subsystem such as genset, motor, excitation system, prime mover and governing system thereof, aut.eq., protective device, supervising device, commutation inversion device or device, calculate in real time in control system or realize calculating or the debugging of Control System Imitation model in analytical calculation, the realistic model formed can use in the modeling of supply station main system; Graphic elements in component library is various general purpose function modules, according to certain system transter block diagram, rear composed component topological relation is connected to graphic elements, component parameters is arranged all can on control system figure, graph of a correspondence element be revised or preserve, component parameters not only comprises the parameter of element itself, also comprise the mark that associate of this graphic elements with certain element on other figure or other figure, the selection, fault verification etc. of component models type, control system calculates data in real time and can monitor on the corresponding graphic elements of this figure; Control system mathematics library provides open component module; user also also passes through the realistic model of user oneself programming realization element by the parameter of the component module definable element opened by the topological relation of the required system of control system mathematics library acquisition, and carries out debugging, analyze or calling under control system mathematics library environment.