CN111898234A - Gas engine power generation system simulation modeling method suitable for FRTDS - Google Patents
Gas engine power generation system simulation modeling method suitable for FRTDS Download PDFInfo
- Publication number
- CN111898234A CN111898234A CN202010422076.9A CN202010422076A CN111898234A CN 111898234 A CN111898234 A CN 111898234A CN 202010422076 A CN202010422076 A CN 202010422076A CN 111898234 A CN111898234 A CN 111898234A
- Authority
- CN
- China
- Prior art keywords
- simulation
- network
- power generation
- generation system
- gas turbine
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 123
- 238000010248 power generation Methods 0.000 title claims abstract description 51
- 238000005094 computer simulation Methods 0.000 title claims abstract description 14
- 238000004088 simulation Methods 0.000 claims abstract description 80
- 238000004364 calculation method Methods 0.000 claims abstract description 19
- 230000008569 process Effects 0.000 claims abstract description 16
- 230000006855 networking Effects 0.000 claims abstract description 15
- 238000012545 processing Methods 0.000 claims description 14
- 230000003993 interaction Effects 0.000 claims description 13
- 230000004044 response Effects 0.000 claims description 7
- 230000001052 transient effect Effects 0.000 claims description 7
- 238000013178 mathematical model Methods 0.000 claims description 4
- 230000008030 elimination Effects 0.000 claims description 3
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 2
- 238000011160 research Methods 0.000 abstract description 6
- 238000004458 analytical method Methods 0.000 abstract description 2
- 239000007789 gas Substances 0.000 description 48
- 239000003990 capacitor Substances 0.000 description 9
- 239000000243 solution Substances 0.000 description 6
- 230000008901 benefit Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000014509 gene expression Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000001939 inductive effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 239000002918 waste heat Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
The invention discloses a gas engine power generation system simulation modeling method suitable for FRTDS, which adopts a real-time digital resolver (FRTDS) based on FPGA to comprehensively perform modeling analysis on models of main components such as a gas engine system, a thermodynamic system, an electric power system and the like and an inverter control system; aiming at the network characteristics of a power system comprising a gas turbine, the simulation method provides a network splitting method and a solving process suitable for the system from the aim of reducing calculated amount and realizing real-time simulation; the multi-rate simulation method provides an implicit difference method based on display prediction and an implicit difference method based on interpolation prediction on the basis of a multi-rate networking method by using a display Euler method, and researches the actual application and simulation process of the multi-rate networking method in a gas turbine power system. The method provided by the invention takes a gas turbine power generation system as an example, solves the real-time simulation problem of a large-scale system with a plurality of simulation step lengths, reduces the simulation calculation amount and enlarges the simulation scale.
Description
Technical Field
The invention belongs to the technical field of power automation, relates to a simulation model establishing method for real-time simulation of a power system, and particularly relates to a simulation modeling method for a gas turbine power generation system suitable for FRTDS.
Background
In the world, energy is regarded as the fundamental condition for human survival and the main driving force for social progress, and the problems of energy shortage, environmental pollution and the like caused by fossil fuels are highly regarded by various countries. Due to the vigorous development of the world economy, under the heavy pressure of great increase of energy demand and increasingly serious contradiction between supply and demand, how to reasonably develop and utilize energy, improve the utilization efficiency of energy, explore more novel energy and improve the economic benefit and ecological benefit of energy utilization becomes a problem of general attention.
The widespread use of micro gas turbines today results in increasingly complex grid operation, making it more difficult for power technicians to monitor, analyze, and control it. Therefore, the research on the digital real-time simulation of the gas turbine power generation system has very important significance for researching the actual operation state of the micro-grid containing the gas turbine power generation system. The digital real-time simulation of the gas turbine power generation system can simulate the on-line early warning, analysis and scheduling of a micro-grid containing a gas turbine in operation in practical application, and can perform hardware-in-loop experiments on secondary equipment and technical training on power system workers.
The real-time solver of the FRTDS is composed of a plurality of micro-processing cores for completing operation tasks, the micro-processing cores complete data interaction through a data exchange station to form a large-scale real-time solver structure, and the whole structure of the real-time solver is shown in figure 1. In order to realize multi-rate simulation and fully utilize the micro-processing core, the micro-processing core and the external equipment perform data interaction in a ping-pong operation mode, and a step length controller masters an interaction time point. The data interaction mode between the microprocessing cores comprises a hand-pulling mode and a data pipeline mode. Data interaction between two adjacent microprocessing cores adopts a hand-in-hand form, and data interaction between two microprocessing cores far away from each other adopts a data pipeline form. Meanwhile, an SFP/SFP + interface and a PCIe interface are equipped for the real-time resolver, the SFP/SFP + interface is connected with the actual equipment through a signal conversion device, and the PCIe interface is connected with the industrial controller.
The micro-processing core is composed of an arithmetic component, a data storage unit, an instruction unit, an input data controller and an output data controller, as shown in fig. 2. The arithmetic component is used for completing the execution of complex operation expressions including arithmetic operation expressions, logic operation expressions and comparison operation expressions. The instruction unit comprises an instruction storage unit and an instruction processing unit, and the instruction processing unit is responsible for obtaining instructions from the instruction storage unit and sorting the instructions into data structures required by the control unit through decomposition and indexing. The control unit is composed of a selection controller, an input data controller and an output data controller.
When multi-rate simulation is carried out in the FRTDS platform, the interaction frequency of large-step simulation operation and small-step simulation operation is different from that of external data. At this time, the information interaction between the employed arithmetic element and the input/output element is as shown in fig. 3. The modules 1-4 are all ping-pong circuits, wherein the module 1 and the module 3 perform ping-pong operations according to large step long time beats, and the module 2 and the module 4 perform ping-pong operations according to small step time beats. And the RAMA and the RAMB are sequentially connected with the input-output assembly and the operation assembly.
However, the existing research on the gas turbine simulation method has the following defects:
1. in general, the research on cogeneration considered from the integrated energy system ignores the establishment process and the real operation state of the actual model of the system, and is not suitable for the simulation aiming at the actual operation, such as real-time simulation.
2. In the research, the working principle and the dynamic characteristic of the single gas turbine power generation system or the single link of single waste heat utilization are mostly modeled and analyzed, and the research on the simulation method of the micro gas turbine thermoelectric combined operation characteristic is still insufficient.
Disclosure of Invention
The invention aims to overcome the defects in the prior art, provides a gas turbine power generation system simulation modeling method suitable for FRTDS, solves the problem of real-time simulation of a large-scale system with a plurality of simulation step lengths, reduces the simulation calculation amount, and enlarges the simulation scale.
The purpose of the invention is realized by the following technical scheme:
a gas engine power generation system simulation modeling method suitable for FRTDS is characterized in that a real-time digital solver (FRTDS) based on FPGA is composed of a plurality of micro-processing cores for completing operation tasks, the micro-processing cores complete data interaction through a data exchange station and can be applied to multi-rate simulation; the method comprises the following steps:
(1) establishing a mathematical model of a gas engine power generation system; the gas engine power generation system consists of an electric power system, a mechanical system and a thermal system;
(2) investigating the dynamic characteristics of each element in the gas engine power generation system, describing a single element by adopting an ordinary differential equation or a partial differential equation, discretizing, carrying out simulation calculation by using a discrete network method, and writing a network equation according to a connection relation series; establishing a network computing method and a solving process suitable for a gas turbine power generation system by combining the network characteristics of a power system containing a gas turbine;
(3) the method takes the transient response speed and the difference of calculated amount of different elements in the gas turbine power generation system into consideration, carries out multi-rate networking on the whole gas turbine power generation system during real-time simulation, and reduces the calculated amount of the gas turbine power generation system on the basis of meeting the response speed.
Further, the network computing method and the solving process adopted in the step (2) are as follows: network division is carried out on a power system part in the gas engine power generation system to obtain a multi-port Norton equivalent circuit of a sub-network, namely an equivalent sub-network, and then a constraint equation of voltage and current is written according to the connection relation of the equivalent sub-network; during real-time simulation, a node voltage method is adopted only for a gas engine power generation system formed by connecting standard Noton equivalent circuits, a node voltage elimination method is adopted to solve port input variables of each equivalent sub-network, and port node voltages of each equivalent sub-network at any moment are obtained; when the equivalent sub-networks are internally solved, the norton equivalent currents, the injection currents of the port nodes and the internal electrical quantities of the equivalent sub-networks of each equivalent sub-network can be calculated by using a linear combination method for updating.
Furthermore, a rectifier consisting of diodes, and an inverter consisting of an IGBT and a freewheeling diode which are connected in parallel in the gas turbine power generation system belong to nonlinear elements; describing the state of the power switch by adopting a binary resistance model, and repeatedly modifying a network equation and solving the network equation in the process of determining the state of the switch; the matrix dimension of the network equation before and after the power switch action is kept constant.
Furthermore, an implicit difference method based on display prediction and an implicit difference method based on interpolation prediction are provided on the basis of a common multi-rate networking method, namely a display Eulerian method, and the simulation process of the actual application and multi-rate simulation of the implicit difference method based on display prediction and the implicit difference method based on interpolation prediction in the multi-rate networking of the power system of the gas turbine is simulated.
Further, solving the equivalent current source of the equivalent sub-network through linear combination, wherein the equivalent current source is formed by linear combination of the previous time value of the state variable of each dynamic element and the current time value of the independent current source; and for the internal state variable of the equivalent sub-network, the voltage value of the network port is replaced by the independent voltage source value, and the current time value of the state variable of the output dynamic element is expressed as a linear combination of the previous time value of the state variable, the current time value of the independent current source and the current time voltage value of the port node of the equivalent sub-network.
Further, for first order differential equations
At tnTo tn+1Within the integration step length, x can be obtained by adopting display Euler method differencen+1＝xn+Δtf(xn,tn) (ii) a When the implicit trapezoidal method is used for multi-rate networking, the first order differential equation is processed at tnTo tn+1The difference can be obtained by adopting an implicit trapezoidal method in the integral step length
Implicit difference method based on display prediction, namely, firstly solving unknown variable f (x) in implicit trapezoidal method by explicit Euler methodn+1,tn+1) Predicting, and solving by using predicted data through an implicit trapezoidal method; implicit difference method based on interpolation prediction, namely, firstly solving unknown variable f (x) in subsequent implicit trapezoidal method by Lagrange interpolation polynomialn+1,tn+1) And predicting, and solving by using the predicted data through an implicit trapezoidal method.
Compared with the prior art, the technical scheme of the invention has the following beneficial effects:
1. a simulation model of a gas turbine power generation system containing combined heat and power is established, and a network calculation method and a solution method are designed for the system type from the purposes of reducing calculated amount and realizing real-time simulation aiming at the network characteristics of a power system containing a gas turbine.
2. The method aims at processing the difference of element transient response speeds in multi-system simulation, adopts multi-rate simulation on a simulation method, optimizes the traditional electromagnetic transient simulation, improves the defects that the traditional electromagnetic transient simulation cannot reflect the change trend of the whole system along with time, has small step length and slow simulation speed, and is not suitable for real-time simulation, and realizes the real-time simulation of a micro gas turbine simulation model.
3. The multi-rate networking method provides a multi-rate interface interactive data prediction method, and improves the multi-rate real-time simulation precision.
4. And a network equation is established by adopting a network block solving method, so that the simulation calculation amount is reduced, and the simulation scale is conveniently enlarged.
5. The problems of large calculation amount and long serial length of the traditional electromagnetic transient algorithm can be effectively solved.
Drawings
Fig. 1 is a schematic diagram of the overall structure of a real-time resolver.
FIG. 2 is a schematic diagram of a microprocessor core architecture.
FIG. 3 refers to the interaction between the compute component and the input-output component.
FIG. 4 is an overall model of the combined heat and power system of a micro gas turbine according to the present invention.
Fig. 5 is a network block solving form of the present invention.
FIG. 6 is a micro-grid simulation flow chart of the gas turbine-containing power generation system of the invention.
Fig. 7(a) and 7(b) are schematic diagrams of an implicit differential method decoupling LCL branch of explicit prediction according to the present invention.
FIG. 8 is a gas turbine power generation system subsystem of the present invention split.
FIG. 9 is a simulation time node between two adjacent systems of the present invention.
Detailed Description
The invention is described in further detail below with reference to the figures and specific examples. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
A simulation modeling method for a gas turbine power generation system suitable for FRTDS adopts a real-time digital solver (FRTDS) based on FPGA, as shown in figure 1, and comprises a plurality of micro-processing cores for completing operation tasks, the micro-processing cores complete data interaction through a data exchange station and can be applied to multi-rate simulation,
in this embodiment, the implementation of the gas turbine power generation system simulation modeling suitable for FRTDS requires the following steps:
(1) the establishment of a mathematical model of the gas engine power generation system is the basis of real-time simulation. In the FRTDS real-time solver, it is necessary to complete the simulation operation of the power system, the mechanical system, and the thermal system, which are composed of the gas engine system, the permanent magnet synchronous generator, the rectifier inverter, the circuit, the thermocouple unit, the heat supply loop, and the like, and it is necessary to determine the whole mathematical model, as shown in fig. 4.
(2) And (3) investigating the dynamic characteristics of the elements, describing a single element by adopting a normal differential equation or a partial differential equation, discretizing by adopting an implicit trapezoidal method, carrying out simulation calculation by using a discrete network method, and writing a network equation according to a connection relation series. The network computing method and the solving process suitable for the gas turbine power generation system are provided by combining the network characteristics of the power system with the gas turbine and the network characteristics of the power system with the gas turbine, and from the aim of reducing the computation amount and realizing real-time simulation.
(3) Considering the difference between the transient response speed and the calculated amount of different elements, the multi-rate networking is carried out on the whole gas turbine power generation system during real-time simulation, and the calculated amount of system operation is reduced on the basis of meeting the response speed.
In this embodiment, the network computing method and solving process adopted by the gas turbine power generation system are as follows: when a large-scale network is solved, the network division is carried out on the part of the power system in the gas engine power generation system to obtain the multi-port Nuoton equivalent circuits of the sub-networks, and then the constraint equation of the voltage and the current is written according to the connection relation of each equivalent sub-network. During real-time simulation, a node voltage method is only needed to be adopted for a power system formed by connecting standard Noton equivalent circuits, and a node voltage elimination method is adopted to solve port input variables of each sub-network, so that the port node voltage of each equivalent sub-network at any moment is obtained. When the equivalent sub-networks are solved internally, the norton equivalent currents, the injection currents of the port nodes and the internal electrical quantities of the equivalent sub-networks of each equivalent sub-network can be calculated by using a linear combination method to update. The formation of sub-networks of norton equivalent circuits and the resulting connected network blocks is shown in fig. 5.
In the embodiment, a rectifier composed of diodes, an inverter composed of an IGBT and a freewheeling diode in parallel in a gas turbine power generation system belong to nonlinear elements, the circuit topology is affected when a high-frequency power switch acts, and a binary resistance model is used for describing the state of the power switch, so that the matrix dimension of a network equation before and after the power switch acts is kept constant. The state of the switching element affects the external characteristic equations of the sub-network and thus the input variables of the port. Therefore, in determining the switch state, the network equations need to be iteratively modified and solved. This method takes a lot of computation time. When the port input variable of the sub-network is the current in the inductor or the voltage at the capacitor, the change of the state of the switch element inside the sub-network has little influence on the external characteristic equation because the inductor current and the capacitor voltage cannot change suddenly, thereby avoiding the whole network equation to be solved again. Under this assumption, the local iteration method shown in fig. 6 can be used to perform the microgrid simulation calculation.
In this embodiment, the multi-rate networking method is to provide an implicit difference method based on display prediction and an implicit difference method based on interpolation prediction on the basis of a common multi-rate networking method (display euler method), and study the actual application and the simulation process of multi-rate simulation of the multi-rate networking method in the gas turbine power system.
The functions and effects of the present invention are further explained below:
(1) implicit difference method based on interpolation prediction.
For a complex gas engine power generation system, certain single elements are described by ordinary differential or partial differential equations, and the state equation is
When the simulation step length is delta t and the implicit trapezoidal method is adopted to carry out difference on the above formula, the method can be obtained
x(t)＝H0u(t)+H1x(t-Δt)+H0u(t-Δt)
Wherein,
according to the formula, the value of x at the time t is not only determined by u and x of t-delta t but also related to u at the time t, so that if u (t) is obtained first to obtain x (t), the system can be decoupled into two subsystem circuits, and simulation is carried out by using different step sizes through a multi-rate technology. However, the value of u at the current time is often an unknown variable and needs to be predicted by an interpolation method.
When Lagrange interpolation polynomial is adopted to predict u (t), n times t before t time are taken0＝t-kΔt， t1＝t-(k-1)Δt，…，tn-1Data u (t) of t- Δ t0)，u(t1)，…，u(tn-1) Substituting the interpolation basis function to obtain
Wherein i is 0, 1.
the value of u (t) at time t can be obtained
And (3) substituting the u (t) predicted by the formula (4-7) for the formula (4-5) to calculate the value of x at the current moment, so as to realize the decoupling operation of the system.
(2) Implicit difference method based on explicit prediction
According to the implicit difference method of interpolation prediction, data needs to be predicted firstly, for some special circuits, an unknown variable in subsequent solution can be predicted through an explicit Euler method, then the predicted data is used for solution through the implicit difference method, and the actual value of the predicted variable is obtained and then enters the subsequent solution process.
Taking the LCL filter circuit as an example, when the capacitance C in fig. 7(a) is differentiated by the implicit difference method, the value of the capacitance voltage at time t is related not only to the capacitance voltage current at the previous time but also to the capacitance current at time t.
Inductor L1The voltage-current relationship between the two ends is obtained by differentiating by adopting an explicit Euler method and delta t as a simulation step length
Wherein i1(t- Δ t) is inductance L1The inductor current at the last moment in time,
is an inductance L1The inductor voltage at the last instant.
The inductance L can be obtained in the same way2Predicted current of
Wherein i2(t- Δ t) is inductance L2The inductor current at the last moment in time,
is an inductance L2The inductor voltage at the last instant.
At this time, the inductance L1And L2The current value at the time t is determined by the inductive voltage and the inductive current of t-delta t, and is independent of the voltage value at the time t, so that the current i flowing through the capacitor can be predicted by adopting an explicit Euler methodC *(t)＝i1 *(t)-i2 *(t)。
Then, the voltage-current relation between the two ends of the capacitor C is differentiated by using an implicit trapezoidal method to obtain
The voltage value of the capacitor differentiated by the above formula at the time t is not only related to the voltage and current at the previous time, but also related to the capacitance current value at the current time obtained by the explicit euler method prediction. Therefore, the solution of the LCL network is divided into five steps: firstly, predicting current i flowing through a capacitor at the current moment by an explicit Euler methodC *(t); ② calculating the current capacitor voltage u by using invisible trapezoid methodC(t); solving the size u of the node voltage value of the resistance-capacitance branch circuit0(t)＝uC(t)+Rd(i1 *(t)-i2 *(t)); fourthly, u is0(t) substitution into the sub-systems 1 and 2 in FIG. 7(b) to obtain the actual inductance current value i at the present time1(t)、i2(t); calculating the actual current i of the current RC branchCAnd (t) solving the capacitor voltage at the next moment. The circuit in fig. 7(a) is equivalent to the circuit in fig. 7(b) in which the decoupling of the two subsystems is completed, and the decoupled subsystem is simulated by using different step sizes through a multi-rate technology.
(3) The partial simulation of the gas turbine power generation system can be split into three subsystems, as shown in FIG. 8. The rectification inversion part comprises a plurality of power electronic switches, the switching frequency of the power electronic switches is very high, and a small step length delta t needs to be selected for simulation; the gas turbine power generation device, the permanent magnet synchronous generator and the power grid connected with the gas turbine can select large step length delta T for simulation, and the calculation amount of simulation is reduced.
Different step lengths are adopted for simulation between two adjacent subsystems, the subsystem 1 and the subsystem 3 adopt larger step length delta T for simulation, the subsystem 2 adopts smaller step length delta T for simulation, and the ratio of delta T to delta T is an integer m. In the time of delta T, the subsystem 1 and the subsystem 3 perform one calculation, the subsystem 2 performs m calculations, the three subsystems adopt parallel calculation time sequences, and simulation time nodes of the three subsystems are shown in FIG. 9.
The large step data need to be sent to the small step side for simulation, the small step side data need to be sent to the large step side for simulation, and interaction is performed once within one delta T.
The mechanical system and the electric system in the micro-grid are separated and are simulated by using a multi-rate technology, the mechanical system, namely the subsystem 1, adopts large-step simulation, the rectification inverter system, namely the subsystem 2, adopts small-step simulation, and the delta T is equal to m delta T.
In many power systems, the rotor of the prime mover directly drags the rotor of the generator to rotate, both having the same rotational angular velocity. And, the rotor equation of motion of the undamped generator is
Where ω is the angular velocity of the rotor, J is the moment of inertia of the rotor, TmMechanical torque, T, output for mechanical systemseIs the electromagnetic torque of the electrical power system. The prime motor and the generator rotor adopt an implicit difference method decoupling circuit based on interpolation prediction.
Differentiating the above formula by an implicit trapezoidal method, and defining the difference value of the mechanical torque and the electromagnetic torque as TdIs obtained by
It follows that the rotational speed of the subsystem 1 is
The rotation speed in the k-th period of the subsystem 2 is
Wherein,
the mechanical torque and the electromagnetic torque at the current moment are obtained through Lagrange interpolation polynomial prediction. In this way, the mechanical and electrical systems of the system can be separated at the rotors of the prime mover and the generator. In the simulation process, the rotating speed used by the subsystem 1 at the current moment T is obtained by the torque difference between the current moment and the moment T-delta T and the rotating speed at the moment T-delta T, and the rotating speed at each moment T-delta T + k delta T of the subsystem 2 is obtained by the mechanical torque at the moments T and T-delta T, the electromagnetic torque at the moment T-delta T + (k-1) delta T and the rotating speed at the moment T-delta T + (k-1) delta T.
Determining the rotational speed omegar1And ωr2Then, the value of omega is adjustedr1(T) substituting into the mechanical system, and finally obtaining the mechanical torque T at the current moment through the gas turbine and the gas compressor systemmAnd (t) solving the rotating speed at the next moment. Will omegar2And (t) substituting into the subsystem 2 to obtain the output current value of the permanent magnet synchronous generator, substituting the output current value into a node voltage equation of the rectification inverter system, and solving the voltage of each node by eliminating.
And (3) splitting an inverter system and a large power grid system in the microgrid and simulating by using a multi-rate technology. The inverter system and the large power grid system can be separated by using the LCL filter circuit, the inverter system adopts small step simulation, the power grid system adopts large step simulation, and the delta T is h delta T.
Due to the property of the LCL circuit, when the filter circuit is decoupled by adopting an implicit difference method based on explicit prediction, the LCL circuit has the advantages that
u02(t)＝uC2(t)+Rd(i1 *(t)-i2 *(t))
u01(t-ΔT+kΔt)＝uC1(t-ΔT+kΔt)+Rd(i1 *(t-ΔT+kΔt)-i2 *(t))
Wherein u is01Is an equivalent voltage source of the inverter system, u02Is an equivalent voltage source of the power grid system.
The capacitance voltage used by the inverter system for solving is obtained from the capacitance differential formula in (2)
Wherein,
the capacitor voltage used for solving the power grid system is
The present invention is not limited to the above-described embodiments. The foregoing description of the specific embodiments is intended to describe and illustrate the technical solutions of the present invention, and the specific embodiments described above are merely illustrative and not restrictive. Those skilled in the art, having the benefit of this disclosure, may effect numerous modifications thereto and changes may be made without departing from the scope of the invention as defined by the claims and their equivalents.
Claims (6)
1. A gas engine power generation system simulation modeling method suitable for FRTDS is characterized in that a real-time digital solver (FRTDS) based on FPGA is composed of a plurality of micro-processing cores for completing operation tasks, the micro-processing cores complete data interaction through a data exchange station and can be applied to multi-rate simulation; the method is characterized by comprising the following steps:
(1) establishing a mathematical model of a gas engine power generation system; the gas engine power generation system consists of an electric power system, a mechanical system and a thermal system;
(2) investigating the dynamic characteristics of each element in the gas engine power generation system, describing a single element by adopting an ordinary differential equation or a partial differential equation, discretizing, carrying out simulation calculation by using a discrete network method, and writing a network equation according to a connection relation series; establishing a network computing method and a solving process suitable for a gas turbine power generation system by combining the network characteristics of a power system containing a gas turbine;
(3) the method takes the transient response speed and the difference of calculated amount of different elements in the gas turbine power generation system into consideration, carries out multi-rate networking on the whole gas turbine power generation system during real-time simulation, and reduces the calculated amount of the gas turbine power generation system on the basis of meeting the response speed.
2. The simulation modeling method for gas turbine power generation system suitable for FRTDS as claimed in claim 1, wherein the network computing method and solving process adopted in step (2) is as follows: network division is carried out on a power system part in the gas engine power generation system to obtain a multi-port Norton equivalent circuit of a sub-network, namely an equivalent sub-network, and then a constraint equation of voltage and current is written according to the connection relation of the equivalent sub-network; during real-time simulation, a node voltage method is adopted only for a gas engine power generation system formed by connecting standard Noton equivalent circuits, a node voltage elimination method is adopted to solve port input variables of each equivalent sub-network, and port node voltages of each equivalent sub-network at any moment are obtained; when the equivalent sub-networks are internally solved, the norton equivalent currents, the injection currents of the port nodes and the internal electrical quantities of the equivalent sub-networks of each equivalent sub-network can be calculated by using a linear combination method for updating.
3. The gas turbine power generation system simulation modeling method suitable for FRTDS as claimed in claim 1, wherein the gas turbine power generation system comprises a rectifier composed of diodes, an inverter composed of IGBT and a freewheeling diode in parallel, and belongs to a nonlinear element; describing the state of the power switch by adopting a binary resistance model, and repeatedly modifying a network equation and solving the network equation in the process of determining the state of the switch; the matrix dimension of the network equation before and after the power switch action is kept constant.
4. The simulation modeling method for the gas turbine power generation system suitable for the FRTDS as claimed in claim 1, wherein an implicit difference method based on display prediction and an implicit difference method based on interpolation prediction are provided on the basis of a common multi-rate networking method, namely a display Euler method, and the simulation process of actual application and multi-rate simulation of the implicit difference method based on display prediction and the implicit difference method based on interpolation prediction in the multi-rate networking of the gas turbine power system is simulated.
5. The gas turbine power generation system simulation modeling method applicable to FRTDS of claim 2, characterized in that the equivalent current sources of the equivalent sub-network are solved by linear combination, and the equivalent current sources are composed of linear combination of previous time values of state variables of each dynamic element and current time values of independent current sources; for the internal state variable of the equivalent sub-network, the voltage value of the network port is replaced by the independent voltage source value, and the current time value of the state variable of the output dynamic element is expressed as the linear combination of the previous time value of the state variable, the current time value of the independent current source and the current time voltage value of the port node of the equivalent sub-network.
6. The simulation modeling method for gas turbine power generation system suitable for FRTDS of claim 4, wherein the first order differential equation is applied
At tnTo tn+1X can be obtained by adopting display Euler method difference in integration step lengthn+1＝xn+Δtf(xn,tn) (ii) a When the implicit trapezoidal method is used for multi-rate networking, the first order differential equation is processed at tnTo tn+1The difference can be obtained by adopting an implicit trapezoidal method in the integration step length
Implicit difference method based on display prediction, namely, firstly solving unknown variable f (x) in implicit trapezoidal method by explicit Euler methodn+1,tn+1) Predicting, and solving by using predicted data through an implicit trapezoidal method; implicit difference method based on interpolation prediction, namely, firstly solving unknown variable f (x) in subsequent implicit trapezoidal method by Lagrange interpolation polynomialn+1,tn+1) And predicting, and solving by using predicted data through an implicit trapezoidal method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010422076.9A CN111898234A (en) | 2020-05-18 | 2020-05-18 | Gas engine power generation system simulation modeling method suitable for FRTDS |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010422076.9A CN111898234A (en) | 2020-05-18 | 2020-05-18 | Gas engine power generation system simulation modeling method suitable for FRTDS |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111898234A true CN111898234A (en) | 2020-11-06 |
Family
ID=73206519
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010422076.9A Pending CN111898234A (en) | 2020-05-18 | 2020-05-18 | Gas engine power generation system simulation modeling method suitable for FRTDS |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111898234A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113155888A (en) * | 2021-04-25 | 2021-07-23 | 山东大学 | Dynamic equivalent circuit of combined heat and power system and working method thereof |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101478159A (en) * | 2009-01-12 | 2009-07-08 | 浙江大学 | Transient stabilized constraint tide optimization process based on differential equality constraint conversion order reduction |
| CN106372339A (en) * | 2016-09-05 | 2017-02-01 | 清华大学 | Multi-rate simulation method and device for power-electronizing power system |
| CN110059354A (en) * | 2019-03-14 | 2019-07-26 | 天津大学 | A kind of electric system real-time electromechanical transient emulation method based on FPGA |
| CN110135031A (en) * | 2019-04-30 | 2019-08-16 | 东南大学 | Electric power system transient stability calculation method based on half implicit runge kutta method |
-
2020
- 2020-05-18 CN CN202010422076.9A patent/CN111898234A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101478159A (en) * | 2009-01-12 | 2009-07-08 | 浙江大学 | Transient stabilized constraint tide optimization process based on differential equality constraint conversion order reduction |
| CN106372339A (en) * | 2016-09-05 | 2017-02-01 | 清华大学 | Multi-rate simulation method and device for power-electronizing power system |
| CN110059354A (en) * | 2019-03-14 | 2019-07-26 | 天津大学 | A kind of electric system real-time electromechanical transient emulation method based on FPGA |
| CN110135031A (en) * | 2019-04-30 | 2019-08-16 | 东南大学 | Electric power system transient stability calculation method based on half implicit runge kutta method |
Non-Patent Citations (4)
| Title |
|---|
| TOSHIJI KATO ET.AL.: "Multirate Analysis Method for a Power Electronic System by Circuit Partitioning", 《IEEE》 * |
| 易桂平，胡仁杰: "微型燃气轮机发电建模与仿真研究", 《江苏电机工程》 * |
| 杨俊鹏: "微电网的暂态建模与仿真", 《工程科技Ⅱ辑》 * |
| 赵丹: "微电网电磁暂态实时仿真并行计算方法研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113155888A (en) * | 2021-04-25 | 2021-07-23 | 山东大学 | Dynamic equivalent circuit of combined heat and power system and working method thereof |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109066688A (en) | Based on the Probabilistic Load Flow data capture method under renewable energy uncertainty | |
| CN103700036B (en) | A kind of transient stability projecting integral method being suitable to power system Multiple Time Scales | |
| CN104392135A (en) | Probabilistic optimal power flow calculation method for alternating-current and direct-current systems of offshore wind power plants subjected to VSC-HVDC (voltage source converter-high voltage direct current) grid connection | |
| CN104765282B (en) | A kind of digital Synthetic Simulation Platform in power station | |
| CN111598376A (en) | Method and system for carrying out auxiliary decision-making on power grid big data based on information driving | |
| Xiong et al. | Real-time FPGA-digital twin monitoring and diagnostics for PET applications | |
| CN103207924B (en) | Wind power generating set machinery transient state and electric transient hybrid simulation system and method | |
| CN111898234A (en) | Gas engine power generation system simulation modeling method suitable for FRTDS | |
| Bai et al. | Automatic modeling and optimization for the digital twin of a regional multi-energy system | |
| Gong et al. | A model review for controller-hardware-in-the-loop simulation in EV powertrain application | |
| Benhamadouche et al. | FPGA-based hardware-in-the-loop for multi-domain simulation | |
| Zhang et al. | Research on digital twin model of energy equipment | |
| CN117272110A (en) | Running state fault diagnosis and prediction method for fan generator and gearbox | |
| Zhou et al. | Digital Twin Modeling and Simulation of Distributed and Centralized Integrated Energy System on CloudPSS-IESLab | |
| Poppenborg et al. | Energy hub gas: A multi-domain system modelling and co-simulation approach | |
| Wang et al. | Architecture and application of real-time co-simulation platform for cyber-physical power system | |
| CN112307699A (en) | VSC-HVDC electromechanical transient simulation method and system based on phase component | |
| CN104933255A (en) | Combined simulation method of calling MATLAB engine by PSASP from background | |
| Sun et al. | Research on Carbon Traceability of Power System Based on Blockchain and Power Flow Calculation under Carbon Peaking and Carbon Neutrality Goals | |
| Dong et al. | Control configuration based on AOE and its application in the power allocation of BESSs | |
| CN113706094B (en) | Comprehensive energy real-time collaborative simulation system and method based on message bus | |
| CN112560245B (en) | Simulation calculation method for configuring switching device and heat dissipation equipment of electric traction system | |
| CN114386238B (en) | Asynchronous motor semi-physical modeling method considering temperature change | |
| Wang et al. | Implementation of Hardware-in-the-Loop Real-time Simulation of Per-Unit Doubly-fed Induction Generator | |
| CN118211884B (en) | Feature sample construction method for power grid transient stability evaluation and application thereof |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| AD01 | Patent right deemed abandoned | ||
| AD01 | Patent right deemed abandoned |
Effective date of abandoning: 20231117 |