CN102592003A - Data exchange method for electromechanical and electromagnetic transient hybrid simulation - Google Patents
Data exchange method for electromechanical and electromagnetic transient hybrid simulation Download PDFInfo
- Publication number
- CN102592003A CN102592003A CN2011104194286A CN201110419428A CN102592003A CN 102592003 A CN102592003 A CN 102592003A CN 2011104194286 A CN2011104194286 A CN 2011104194286A CN 201110419428 A CN201110419428 A CN 201110419428A CN 102592003 A CN102592003 A CN 102592003A
- Authority
- CN
- China
- Prior art keywords
- simulation
- power grid
- electromagnetic transient
- electromechanical
- mrow
- 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.)
- Granted
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 258
- 230000001052 transient effect Effects 0.000 title claims abstract description 142
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000012545 processing Methods 0.000 claims abstract description 13
- 230000001360 synchronised effect Effects 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 7
- 238000002347 injection Methods 0.000 claims description 5
- 239000007924 injection Substances 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 4
- 238000011160 research Methods 0.000 abstract description 4
- 230000008569 process Effects 0.000 description 13
- 238000010586 diagram Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000000284 extract Substances 0.000 description 2
- 101100499229 Mus musculus Dhrsx gene Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000012549 training Methods 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a data exchange method for electromechanical and electromagnetic transient hybrid simulation and relates to the technical field of transient stability research and digital simulation of power systems. The method comprises the following steps of: simulating an electromagnetic transient simulation power grid to obtain the electromagnetic transient simulation data, simulating an electromechanical simulation power grid to obtain the electromechanical simulation data, processing the electromechanical simulation data to obtain a three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid, correcting a power supply frequency and transmitting the parameters of the three-phase instantaneous value Thevenin equivalent circuit model to the electromagnetic transient simulation power grid through an interface bus; and processing the electromagnetic transient simulation data to obtain the active power and the reactive power injected into the electromechanical simulation power grid through the interface bus and form an electromechanical simulation interface bus model, so that the active power and the reactive power can be transmitted to the electromechanical simulation power grid. According to the data exchange method for electromechanical and electromagnetic transient hybrid simulation disclosed by the invention, the real-time simulation of the power grid in a relatively large scale is realized; and the hybrid simulation is not limited by the scale of the power grid.
Description
Technical Field
The invention relates to the technical field of transient stability research and digital simulation of a power system, in particular to a data exchange method for electromechanical and electromagnetic transient hybrid simulation.
Background
The digital simulation of the power system is widely applied to the fields of transient stability analysis, setting and control of a protection device, real-time training simulation of a dispatcher and the like. The Transient process of the power system is very complex, and can be divided into an Electromagnetic Transient process and an electromechanical Transient process according to different time scales, and the two processes are usually researched by adopting two methods, namely an EMTP (Electromagnetic Transient Program) method and a TSP (Transient simulation Program).
The EMTP calculates the instantaneous response of the voltage current in the grid dynamics with a small simulation step (typically in microsecond steps (h), e.g. 50 μ s). The TSP analysis system has the advantages that the simulation step size of the swing process of the generator power angle and the bus voltage in the dynamic state is large (millisecond step size (H) is generally adopted, for example, 2ms), so that the large-scale power grid can be simulated in real time.
In the process of implementing the invention, the inventor finds that at least the following disadvantages and shortcomings exist in the prior art:
1. EMTP simulation needs a mathematical model for describing the dynamic behavior of elements in a network in detail, so that the real-time simulation of a large-scale power grid is difficult to perform;
2. with the continuous development of power electronic technology and high-voltage direct-current transmission technology, the electromagnetic dynamic process of the novel elements is frequently considered for the safety analysis of a large power grid, but TSP simulation cannot perform the simulation.
Disclosure of Invention
The invention provides a data exchange method for electromechanical and electromagnetic transient hybrid simulation, which combines electromechanical transient simulation and electromagnetic transient simulation together to perform hybrid simulation by the interface method, so as to realize real-time simulation of a large-scale power grid, which is described in detail in the following:
a method of data exchange for hybrid electromechanical and electromagnetic transient simulation, the method comprising the steps of:
(1) regarding an area where an electrical element is located as an electromagnetic transient simulation power grid, performing electromagnetic transient simulation on the electromagnetic transient simulation power grid to obtain electromagnetic transient simulation data, regarding a network except the electromagnetic transient simulation power grid as an electromechanical simulation power grid, performing electromechanical transient simulation on the electromechanical simulation power grid to obtain electromechanical simulation data, and regarding a bus between the electromechanical simulation power grid and the electromagnetic transient simulation power grid as an interface bus of hybrid simulation;
(2) processing the obtained electromechanical simulation data to obtain a three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid, correcting the power frequency, and transmitting parameters of the three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid to the electromagnetic transient simulation power grid through the interface bus;
(3) and processing the obtained electromagnetic transient simulation data, obtaining active power P and reactive power Q of the interface bus injected into an electromagnetic transient simulation power grid, and forming an electromechanical simulation model by the active power P and the reactive power Q obtained by each interface to transmit to the electromechanical simulation power grid.
The step (2) of processing the obtained electromechanical simulation data to obtain a three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid, modifying the power frequency, and transmitting parameters of the three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid to the electromagnetic transient simulation power grid through the interface bus comprises the following specific contents:
1) obtaining the voltage (V) of each synchronous generator in the d-q coordinate system from the electromechanical simulation datad、Vq) Current value (I)d、Iq) Network node admittance matrix Y, node voltage
And injecting current
2) According to the voltage (V)d、Vq) And the current value (I)d、Iq) Obtaining the equivalent generator admittance Y in the unified circuit model of each generatorGAnd injecting current
And the equivalent generator admittance YGAnd the injection current
Adding the power grid into the electromechanical simulation power grid;
3) according to the network node admittance matrix Y and the node voltage
And the injection current
Acquiring a single-phase multi-port Thevenin phasor equivalent circuit of the electromechanical simulation power grid by a network equivalent method;
4) converting the single-phase multi-port Thevenin phasor equivalent circuit of the electromechanical simulation power grid into a three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid;
5) correcting the power supply frequency of the three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid;
6) and transmitting the parameters of the three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid to the electromagnetic transient simulation power grid through the interface bus.
The method for correcting the power frequency of the three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid specifically comprises the following steps:
fjthe frequency of the Thevenin equivalent power supply of the interface time boundary bus j; f. of0Is a rated frequency; h is an electromechanical simulation step length; thetaj，a(t + H) and θj，aAnd (t) is the phase angle of the jth interface bus voltage at the moment t + H and t respectively.
The unified circuit model specifically comprises:
wherein,
and
equivalent potentials and reactances representing the d-axis and q-axis, respectively; r represents the stator resistance.
The specific contents of processing the obtained electromagnetic transient simulation data in the step (3), obtaining active power P and reactive power Q injected into an electromagnetic transient simulation power grid through the interface bus, and forming an electromechanical simulation model by the active power P and the reactive power Q obtained by each interface to be transmitted to the electromechanical simulation power grid include:
1) obtaining three-phase instantaneous value voltage (u) of each interface bus from the electromagnetic transient simulation dataa，ub，uc) Injecting a three-phase instantaneous value current vector sequence of the electromagnetic transient simulation power grid;
2) obtaining the fundamental frequency amplitude and the phase angle corresponding to each three-phase instantaneous value current vector sequence by adopting a least square method, and converting the fundamental frequency amplitude and the phase angle into three-phase phasors
And
3) applying a symmetrical component method to the three-phase phasor
Andperforming coordinate transformation to extract positive sequence component of node voltage at the interface bus
And injecting a current positive sequence component of the electromagnetic transient simulation grid
4) According to the positive sequence component of the node voltage of each interface bus
And a current positive sequence component of the electromagnetic transient simulation power grid
Obtaining the interface bus bar noteEntering active power P and reactive power Q of the electromagnetic transient simulation power grid;
5) and forming an electromechanical simulation model by the active power P and the reactive power Q obtained by the interfaces, and transmitting the electromechanical simulation model to the electromechanical simulation power grid.
The technical scheme provided by the invention has the beneficial effects that:
the invention provides a data exchange method for electromechanical and electromagnetic transient hybrid simulation, which processes and transmits electromechanical transient data and electromagnetic transient hybrid simulation data, realizes hybrid simulation combining electromechanical transient simulation and electromagnetic transient simulation, and realizes real-time simulation of a large-scale power grid; the hybrid simulation is not limited by the scale of a power grid, electromagnetic transient hybrid simulation can be performed on the power grid inside the hybrid simulation, electromechanical transient simulation is performed on the external power grid, and meanwhile, the simulation efficiency and the electromagnetic transient result of the internal power grid research are considered; the hybrid simulation can better reflect the influence of the power electronic equipment and the high-voltage direct-current transmission equipment on the transient power angle stability of the power system after being connected to the power grid, and the electromagnetic simulation result can reflect the influence of the power angle swing of the system in return; moreover, the unified circuit model of the synchronous generator provided by the invention expands the application range of the hybrid simulation model; the three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid considering the frequency deviation can correctly reflect the physical process that the voltage and the frequency of a generator rotor and each boundary bus swing after the system is disturbed; the invention extracts the voltage and current fundamental frequency components by using the least square technique and the symmetrical component method, and can carry out mixed simulation on the asymmetric faults.
Drawings
FIG. 1 is a schematic diagram of a unified circuit model of a synchronous generator provided by the present invention;
FIG. 2 is a schematic diagram of a hybrid simulation data exchange process provided by the present invention;
FIG. 3 is a wiring diagram of the 3-machine 9-node checking system with the transformer substation provided by the invention;
FIG. 4 is a schematic diagram of an EMTP simulation circuit model provided in the present invention;
fig. 5 is a schematic diagram of a TSP simulation circuit model provided by the present invention;
FIG. 6 is a schematic diagram of a comparison curve of a mixed simulation result of the three-phase current of B1-B and a PSCAD full-network electromagnetic simulation;
FIG. 7 is a flow chart of a data exchange method for electromechanical and electromagnetic transient hybrid simulation according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention will be described in detail with reference to the accompanying drawings.
In order to combine the electromechanical transient simulation and the electromagnetic transient simulation for hybrid simulation, an embodiment of the present invention provides a data exchange method for electromechanical and electromagnetic transient hybrid simulation, which is described in detail in fig. 2 and 7:
101: regarding the area where the electrical element is located as an electromagnetic transient simulation power grid, performing electromagnetic transient simulation on the electromagnetic transient simulation power grid to obtain electromagnetic transient simulation data, regarding a network except the electromagnetic transient simulation power grid as an electromechanical simulation power grid, performing electromechanical transient simulation on the electromechanical simulation power grid to obtain electromechanical simulation data, and regarding a bus between the electromechanical simulation power grid and the electromagnetic transient simulation power grid as an interface bus for hybrid simulation;
102: processing the obtained electromechanical simulation data to obtain a three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid, correcting the power frequency, and transmitting parameters of the three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid to the electromagnetic transient simulation power grid through an interface bus;
wherein, this step specifically includes:
1) obtaining the voltage (V) of each synchronous generator in d-q coordinate system from the electromechanical simulation datad、Vq) Current value (I)d、Iq) Network node admittance matrix Y, node voltageAnd injecting current
Wherein the d-axis represents the direct axis of the synchronous generator and the q-axis represents the quadrature axis of the synchronous generator.
2) According to voltage (V)d、Vq) Sum current value (I)d、Iq) Obtaining the equivalent generator admittance Y in the unified circuit model of each generatorGAnd injecting currentAnd the equivalent generator is admitted to YGAnd injecting current
Adding the power grid into an electromechanical simulation power grid;
the unified circuit model in the embodiment of the present invention is the model shown in fig. 1. For different kinds of generators, the unified stator voltage equation in the d-q coordinate system can be expressed as:
in the formula of Ud、Uq、IdAnd IqRepresenting the voltage and current of d-axis and q-axis of the synchronous motor respectively;
and
equivalent (transient or sub-transient) potentials and reactances representing the d-and q-axes, respectively; r represents the stator resistance.
In the embodiment of the invention, the synchronous motor is connected with the electromechanical simulation power grid through a coordinate transformation equation shown in a formula (2), wherein delta is a power angle of the generator, and U isx+jUyAnd Ix+jIyThe terminal voltage and current phasor of the generator under the system x-y coordinate system are shown in the figure 1.
Substituting the formula (2) into the formula (1) to derive a model of the synchronous generator in an x-y coordinate system, which is as follows:
in the formula:
after the TSP subsystem performs the one-step simulation, the hybrid simulation component receives the synchronous generator data and calculates parameters of the synchronous generator circuit model according to equation (3), and adds the parameters to the network node admittance matrix and injection current data received by the component.
3) According to network node admittance matrix Y, node voltage
And injecting current
Acquiring a single-phase multi-port Thevenin phasor equivalent circuit of the electromechanical simulation power grid by a network equivalent method;
4) converting a single-phase multi-port Thevenin phasor equivalent circuit of the electromechanical simulation power grid into a three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid;
because the voltage frequency of each boundary bus in the power system dynamic state can swing near a rated value, the power frequency of a three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid needs to be corrected.
5) Correcting the power frequency of a three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid;
wherein the steps are as follows: the power supply frequency is revised according to the formula (4),
fjthe frequency of the Thevenin equivalent power supply of the interface time boundary bus j; f. of0Is a rated frequency; h is an electromechanical simulation step length; thetaj，a(t + H) and θj，aAnd (t) is the phase angle of the jth interface bus voltage at the moment t + H and t respectively.
6) And transmitting the parameters of the three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid to the electromagnetic transient simulation power grid through an interface bus.
Wherein, the parameters specifically include: voltage amplitude, phase angle, corrected power frequency and impedance, etc.
103: and processing the obtained electromagnetic transient simulation data, obtaining active power P and reactive power Q which are injected into the electromagnetic transient simulation power grid through an interface bus, and forming an electromechanical simulation model by the active power P and the reactive power Q which are obtained by each interface to be transmitted to the electromechanical simulation power grid.
Wherein, this step specifically includes:
1) obtaining three-phase instantaneous value voltage (u) of each interface bus from electromagnetic transient simulation dataa，ub，uc) Injecting a three-phase instantaneous value current vector sequence of the electromagnetic transient simulation power grid;
2) obtaining fundamental frequency amplitude and phase angle corresponding to each three-phase instantaneous value current vector sequence by adopting a least square method, and converting the fundamental frequency amplitude and the phase angle into three-phase phasors
And
in an electromechanical simulation step length H period, electromagnetic transient simulation data is a group of discrete time domain signals y (t) with the corresponding step length Hi). In general, y (t) can be represented by the following Fourier series:
according to the embodiment of the invention, the least square method is adopted to process electromagnetic transient simulation data to obtain the voltage and current discrete signals at the interface bus, and the corresponding fundamental wave voltage and current phasor is solved. Aiming at the characteristics of actual operation of the power system, the formula (5) is simplified, higher harmonics are omitted, and a direct current component, a fundamental frequency and a frequency multiplication component are reserved, so that the formula (6) is obtained.
The above formula is expanded according to a trigonometric formula:
y(t)＝C1F1(t)+C2F2(t) (7)
in the formula:
F1＝cosωt，F2sin ω t. Wherein C is1And C2To await a quantity, F1And F2Is a known time-sequential phasor.
The following error function is established using the least squares method.
Wherein, y (t)i) And the instantaneous three-phase voltage and current transmitted to the interface bus of the hybrid simulation component by the electromagnetic transient simulation power grid are represented.
E＝[[Y]-[F][C]]T[[Y]-[F][C]] (9)
Wherein: [ Y ]]＝[y(t1)，y(t2)，L，y(tn)]T；
[C]＝[C1，C2]T。
The derivation of the above formula is made equal to zero, the unknown quantity C can be solved, and the amplitude A of the base frequency cosine function can be extracted1Angle of sum
By the method, the fundamental phasor of the voltage and the current at the interface bus can be obtainedAnd
3) using symmetrical component method to measure three-phase phasor
And
coordinate transformation is carried out, and the positive sequence component of the node voltage at the bus of the interface is extractedAnd injecting the current positive sequence component of the electromagnetic transient simulation power grid
4) According to the positive sequence component of the node voltage of each interface busAnd current positive sequence component of electromagnetic transient simulation power grid
Obtaining active power P and reactive power Q of an interface bus injected into an electromagnetic transient simulation power grid;
5) and forming an electromechanical simulation model by using the active power P and the reactive power Q obtained by the interfaces, and transmitting the electromechanical simulation model to an electromechanical simulation power grid.
The method comprises the following steps: the active power P and the reactive power Q calculated by each interface are regarded as the constant load of the interface bus TSP simulation, i.e. the electromechanical simulation model, see fig. 5, for example: p1+jQ1And P2+jQ2。
The feasibility of a data exchange method for electromechanical and electromagnetic transient hybrid simulation provided by the embodiment of the present invention is verified with reference to an example, which is described in detail below:
the embodiment of the invention is modified on a 3-machine 9-node experimental system, and a simulation transformer substation is added, wherein the simulation transformer substation comprises two distribution transformers and 4 lines, as shown in figure 3. Performing electromechanical transient simulation on an external network (an electromechanical simulation power grid) of the transformer substation, wherein a generator of the transformer substation selects a double-shaft secondary transient model; an internal network (an electromagnetic transient simulation power grid) of the transformer substation adopts electromagnetic transient simulation; the B1 and B2 buses are interface buses. TSP simulation step size is 0.02s, EMTP simulation step size is 50us, single-phase short-circuit fault occurs at BL2 bus, the fault starts from 1.5 seconds, and 1.7 seconds of fault disappears. And carrying out full-network load flow calculation on the whole experimental system, solving initial information of each state variable in the system, and respectively transmitting the initial information to TSP simulation and EMTP simulation. In order to ensure the stability of the system interface in the actual simulation, the two systems are connected at the moment of 0.5 second. The embodiment of the invention takes the control step of the system simulation by the 1.0-1.02 second hybrid simulation interface component as an example.
Step 1: firstly, receiving and processing data information of electromechanical simulation data in 1.0 second;
the hybrid simulation in the embodiment of the invention comprises two interfaces, and the simulation forms a three-phase instantaneous value Thevenin equivalent circuit model of a corresponding 2-port electromagnetic transient simulation power grid, as shown in FIG. 4, and corresponding data information is shown in Table 1.
Table 1 list of data information transmitted from electromechanical simulation grid to electromagnetic transient simulation grid in 1.0 second
Step 2: the electromechanical simulation power grid transmits the parameter information of the table 1 to the electromagnetic transient simulation power grid through an interface bus;
and step 3: performing electromagnetic simulation on the electromagnetic transient simulation power grid for H seconds (note that H is the step length of TSP simulation); after three-phase voltage and current waveform data obtained by EMTP simulation of the electromagnetic transient simulation power grid are processed, corresponding positive sequence components of voltage and current are extracted, active power P and reactive power Q of each interface are further solved to form an electromechanical simulation model as shown in FIG. 5, and power data information received by the electromechanical simulation power grid is shown in Table 2.
Table 2 data information of electromagnetic transient simulation power grid transferred to electromechanical simulation power grid at 1.0+ H second
And 4, step 4: and the electromagnetic transient simulation power grid transmits the parameter information of the table 2 to the electromechanical simulation power grid and carries out one-step electromechanical transient simulation.
And performing mixed simulation according to the steps until the simulation end time is 3 seconds. Fig. 6 shows three-phase currents of a line B1-B, where a solid line is a hybrid simulation result and a dotted line is a commercial software PSCAD full-network electromagnetic simulation result, and a comparison shows that two curves can be well fitted, which indicates that the hybrid simulation data exchange method is accurate and effective, and verifies the feasibility of the embodiment of the present invention.
In summary, the embodiment of the present invention provides a data exchange method for electromechanical and electromagnetic transient hybrid simulation, and the embodiment of the present invention processes and transfers data of electromechanical transient and data of electromagnetic transient hybrid simulation, thereby implementing hybrid simulation combining electromechanical transient simulation and electromagnetic transient simulation, and implementing real-time simulation of a large-scale power grid; the hybrid simulation is not limited by the scale of a power grid, electromagnetic transient hybrid simulation can be performed on the power grid inside the hybrid simulation, electromechanical transient simulation is performed on the external power grid, and meanwhile, the simulation efficiency and the electromagnetic transient result of the internal power grid research are considered; the hybrid simulation can better reflect the influence of the power electronic equipment and the high-voltage direct-current transmission equipment on the transient power angle stability of the power system after being connected to the power grid, and the electromagnetic simulation result can reflect the influence of the power angle swing of the system in return; in addition, the unified circuit model of the synchronous generator provided by the embodiment of the invention expands the application range of the hybrid simulation model; the electromagnetic transient simulation circuit model considering the frequency deviation can correctly reflect the physical process that the voltage and the frequency of a generator rotor and each boundary bus swing after the system is disturbed; the embodiment of the invention extracts the fundamental frequency components of the voltage and the current by using a least square technology and a symmetric component method, and can perform hybrid simulation on asymmetric faults.
Those skilled in the art will appreciate that the drawings are only schematic illustrations of preferred embodiments, and the above-described embodiments of the present invention are merely provided for description and do not represent the merits of the embodiments.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents, improvements and the like that fall within the spirit and principle of the present invention are intended to be included therein.
Claims (5)
1. A method of data exchange for hybrid electro-mechanical and electromagnetic transient simulation, the method comprising the steps of:
(1) regarding an area where an electrical element is located as an electromagnetic transient simulation power grid, performing electromagnetic transient simulation on the electromagnetic transient simulation power grid to obtain electromagnetic transient simulation data, regarding a network except the electromagnetic transient simulation power grid as an electromechanical simulation power grid, performing electromechanical transient simulation on the electromechanical simulation power grid to obtain electromechanical simulation data, and regarding a bus between the electromechanical simulation power grid and the electromagnetic transient simulation power grid as an interface bus of hybrid simulation;
(2) processing the obtained electromechanical simulation data to obtain a three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid, correcting the power frequency, and transmitting parameters of the three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid to the electromagnetic transient simulation power grid through the interface bus;
(3) and processing the obtained electromagnetic transient simulation data, obtaining active power P and reactive power Q of the interface bus injected into an electromagnetic transient simulation power grid, and forming an electromechanical simulation model by the active power P and the reactive power Q obtained by each interface to transmit to the electromechanical simulation power grid.
2. The data exchange method for hybrid electromechanical and electromagnetic transient simulation according to claim 1, wherein the processing of the obtained electromechanical simulation data in step (2) to obtain a three-phase instantaneous value thevenin equivalent circuit model of an electromagnetic transient simulation power grid, modifying a power frequency, and transferring parameters of the three-phase instantaneous value thevenin equivalent circuit model of the electromagnetic transient simulation power grid to the specific contents of the electromagnetic transient simulation power grid through the interface bus comprises:
1) obtaining the voltage (V) of each synchronous generator in the d-q coordinate system from the electromechanical simulation datad、Vq) Current value (I)d、Iq) Network node admittance matrix Y, node voltage
And injecting current
2) According to the voltage (V)d、Vq) And the current value (I)d、Iq) Obtaining the equivalent generator admittance Y in the unified circuit model of each generatorGAnd injecting current
And the equivalent generator admittance YGAnd the injection currentAdding the power grid into the electromechanical simulation power grid;
3) according to the network node admittance matrix Y and the node voltage
And the injection current
Acquiring a single-phase multi-port Thevenin phasor equivalent circuit of the electromechanical simulation power grid by a network equivalent method;
4) converting the single-phase multi-port Thevenin phasor equivalent circuit of the electromechanical simulation power grid into a three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid;
5) correcting the power supply frequency of the three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid;
6) and transmitting the parameters of the three-phase instantaneous value Thevenin equivalent circuit model of the electromagnetic transient simulation power grid to the electromagnetic transient simulation power grid through the interface bus.
3. The data exchange method for the hybrid electromechanical and electromagnetic transient simulation according to claim 2, wherein the modifying the power supply frequency of the three-phase instantaneous value thevenin equivalent circuit model of the electromagnetic transient simulation power grid is specifically:
fjthe frequency of the Thevenin equivalent power supply of the interface time boundary bus j; f. of0Is a rated frequency; h is an electromechanical simulation step length; thetaj，a(t + H) and θj，aAnd (t) is the phase angle of the jth interface bus voltage at the moment t + H and t respectively.
4. The method according to claim 2, wherein the unified circuit model is specifically:
wherein,
and
equivalent potentials and reactances representing the d-axis and q-axis, respectively; r represents the stator resistance.
5. The method according to claim 1, wherein the processing the electromagnetic transient simulation data obtained in step (3), obtaining active power P and reactive power Q injected into an electromagnetic transient simulation grid by the interface bus, and transferring the active power P and the reactive power Q solved by each interface to the electromechanical simulation grid to form an electromechanical simulation model specifically includes:
1) obtaining three-phase instantaneous value voltage (u) of each interface bus from the electromagnetic transient simulation dataa，ub，uc) Injecting a three-phase instantaneous value current vector sequence of the electromagnetic transient simulation power grid;
2) obtaining the fundamental frequency amplitude and the phase angle corresponding to each three-phase instantaneous value current vector sequence by adopting a least square method, and converting the fundamental frequency amplitude and the phase angle into three-phase phasors
And
3) applying a symmetrical component method to the three-phase phasor
Andperforming coordinate transformation to extract positive sequence component of node voltage at the interface bus
And injecting a current positive sequence component of the electromagnetic transient simulation grid
4) According to the positive sequence component of the node voltage of each interface busAnd a current positive sequence component of the electromagnetic transient simulation power grid
Obtaining active power P and reactive power Q of the electromagnetic transient simulation power grid injected by the interface bus;
5) and forming the active power P and the reactive power Q obtained by each interface into the electromechanical simulation model, and transmitting the electromechanical simulation model to the electromechanical simulation power grid.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2011104194286A CN102592003B (en) | 2011-12-15 | 2011-12-15 | Data exchange method for electromechanical and electromagnetic transient hybrid simulation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2011104194286A CN102592003B (en) | 2011-12-15 | 2011-12-15 | Data exchange method for electromechanical and electromagnetic transient hybrid simulation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102592003A true CN102592003A (en) | 2012-07-18 |
| CN102592003B CN102592003B (en) | 2013-11-13 |
Family
ID=46480637
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2011104194286A Expired - Fee Related CN102592003B (en) | 2011-12-15 | 2011-12-15 | Data exchange method for electromechanical and electromagnetic transient hybrid simulation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102592003B (en) |
Cited By (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102841962A (en) * | 2012-07-23 | 2012-12-26 | 南方电网科学研究院有限责任公司 | Electromagnetic-electromechanical transient hybrid real-time simulation interface interaction method based on advanced calculation |
| CN102891483A (en) * | 2012-09-07 | 2013-01-23 | 中国南方电网有限责任公司电网技术研究中心 | Electromagnetic and electromechanical transient digitalized hybrid real-time simulation method for power system |
| CN103279590A (en) * | 2013-04-22 | 2013-09-04 | 中国南方电网有限责任公司电网技术研究中心 | Initial self-correction computation method of interface power in electrical power system hybrid real-time simulation |
| CN104732033A (en) * | 2015-03-31 | 2015-06-24 | 河南行知专利服务有限公司 | Simulation analysis method of electromechanical transient process of electrical power system |
| CN105160137A (en) * | 2015-09-30 | 2015-12-16 | 中国南方电网有限责任公司电网技术研究中心 | PSCAD based electromagnetic-electromechanical transient mixed simulation realization method |
| CN105259781A (en) * | 2015-09-25 | 2016-01-20 | 上海交通大学 | Electric power system mixed simulation system for different types of real-time digital simulation apparatuses |
| CN106204331A (en) * | 2016-07-19 | 2016-12-07 | 北京智中能源互联网研究院有限公司 | Embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation method |
| WO2016188503A3 (en) * | 2015-05-22 | 2016-12-29 | 中国电力科学研究院 | Implementation method suitable for electromagnetic transient multi-time scale real-time simulation interface |
| CN107066732A (en) * | 2017-04-13 | 2017-08-18 | 清华大学 | Towards the distributed hybrid simulation analysis method and device of alternating current-direct current hybrid network |
| CN107622152A (en) * | 2017-09-01 | 2018-01-23 | 南方电网科学研究院有限责任公司 | Electromechanical-electromagnetic transient hybrid simulation method for large power grid |
| CN107947155A (en) * | 2017-11-14 | 2018-04-20 | 南方电网科学研究院有限责任公司 | Hybrid simulation method based on electromechanical transient and electromagnetic transient |
| CN108009338A (en) * | 2017-11-27 | 2018-05-08 | 国网北京市电力公司 | Transient emulation processing method and processing device |
| CN109066612A (en) * | 2018-08-27 | 2018-12-21 | 国网湖北省电力有限公司电力科学研究院 | The switching overvoltage calculation method combined based on electromechanical transient and electro-magnetic transient |
| CN109116183A (en) * | 2018-06-28 | 2019-01-01 | 全球能源互联网研究院有限公司 | Harmonic-model parameter identification method, device, storage medium and electronic equipment |
| CN109460575A (en) * | 2018-10-09 | 2019-03-12 | 武汉大学 | A kind of hybrid simulation interface system and interpolation algorithm based on Bei Jielong equivalence route |
| CN109783895A (en) * | 2018-12-27 | 2019-05-21 | 中国电力科学研究院有限公司 | A kind of electric system hybrid simulation method and system containing distributed generation resource |
| CN109815607A (en) * | 2019-01-31 | 2019-05-28 | 上海交通大学 | Full electromagnetical transient emulation method towards extensive AC network based on time-frequency coordinate transform |
| CN112036002A (en) * | 2020-07-02 | 2020-12-04 | 中国电力科学研究院有限公司 | Electromechanical transient-electromagnetic transient hybrid simulation method and system |
| CN112231981A (en) * | 2020-10-20 | 2021-01-15 | 国网宁夏电力有限公司电力科学研究院 | Method for establishing large-scale electromagnetic transient simulation example |
| CN113541189A (en) * | 2021-07-19 | 2021-10-22 | 国网上海市电力公司 | Interface positioning method suitable for large-scale alternating current-direct current power grid electromechanical-electromagnetic transient hybrid simulation |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6807520B1 (en) * | 2000-12-11 | 2004-10-19 | Synopsys, Inc. | System and method for simulation of an integrated circuit design using a hierarchical input netlist and divisions along hierarchical boundaries thereof |
| CN101246505A (en) * | 2007-12-14 | 2008-08-20 | 南方电网技术研究中心 | Electric network electromagnet transient and electromechanical transient hybrid simulation system and simulation method thereof |
| CN101957872A (en) * | 2010-07-09 | 2011-01-26 | 南方电网技术研究中心 | Hybrid real-time simulation method for alternating current-direct current power system |
-
2011
- 2011-12-15 CN CN2011104194286A patent/CN102592003B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6807520B1 (en) * | 2000-12-11 | 2004-10-19 | Synopsys, Inc. | System and method for simulation of an integrated circuit design using a hierarchical input netlist and divisions along hierarchical boundaries thereof |
| CN101246505A (en) * | 2007-12-14 | 2008-08-20 | 南方电网技术研究中心 | Electric network electromagnet transient and electromechanical transient hybrid simulation system and simulation method thereof |
| CN101957872A (en) * | 2010-07-09 | 2011-01-26 | 南方电网技术研究中心 | Hybrid real-time simulation method for alternating current-direct current power system |
Cited By (28)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102841962A (en) * | 2012-07-23 | 2012-12-26 | 南方电网科学研究院有限责任公司 | Electromagnetic-electromechanical transient hybrid real-time simulation interface interaction method based on advanced calculation |
| CN102891483A (en) * | 2012-09-07 | 2013-01-23 | 中国南方电网有限责任公司电网技术研究中心 | Electromagnetic and electromechanical transient digitalized hybrid real-time simulation method for power system |
| CN103279590A (en) * | 2013-04-22 | 2013-09-04 | 中国南方电网有限责任公司电网技术研究中心 | Initial self-correction computation method of interface power in electrical power system hybrid real-time simulation |
| CN103279590B (en) * | 2013-04-22 | 2016-05-11 | 中国南方电网有限责任公司电网技术研究中心 | The initial self-correcting computational methods of interface power in power system hybrid real-time simulation |
| CN104732033A (en) * | 2015-03-31 | 2015-06-24 | 河南行知专利服务有限公司 | Simulation analysis method of electromechanical transient process of electrical power system |
| WO2016188503A3 (en) * | 2015-05-22 | 2016-12-29 | 中国电力科学研究院 | Implementation method suitable for electromagnetic transient multi-time scale real-time simulation interface |
| CN105259781A (en) * | 2015-09-25 | 2016-01-20 | 上海交通大学 | Electric power system mixed simulation system for different types of real-time digital simulation apparatuses |
| CN105160137A (en) * | 2015-09-30 | 2015-12-16 | 中国南方电网有限责任公司电网技术研究中心 | PSCAD based electromagnetic-electromechanical transient mixed simulation realization method |
| CN106204331B (en) * | 2016-07-19 | 2019-07-12 | 北京智中能源互联网研究院有限公司 | Embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation method |
| CN106204331A (en) * | 2016-07-19 | 2016-12-07 | 北京智中能源互联网研究院有限公司 | Embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation method |
| CN107066732A (en) * | 2017-04-13 | 2017-08-18 | 清华大学 | Towards the distributed hybrid simulation analysis method and device of alternating current-direct current hybrid network |
| CN107622152B (en) * | 2017-09-01 | 2020-10-02 | 南方电网科学研究院有限责任公司 | Electromechanical-electromagnetic transient hybrid simulation method for large power grid |
| CN107622152A (en) * | 2017-09-01 | 2018-01-23 | 南方电网科学研究院有限责任公司 | Electromechanical-electromagnetic transient hybrid simulation method for large power grid |
| CN107947155A (en) * | 2017-11-14 | 2018-04-20 | 南方电网科学研究院有限责任公司 | Hybrid simulation method based on electromechanical transient and electromagnetic transient |
| CN108009338A (en) * | 2017-11-27 | 2018-05-08 | 国网北京市电力公司 | Transient emulation processing method and processing device |
| CN109116183A (en) * | 2018-06-28 | 2019-01-01 | 全球能源互联网研究院有限公司 | Harmonic-model parameter identification method, device, storage medium and electronic equipment |
| CN109116183B (en) * | 2018-06-28 | 2020-07-31 | 全球能源互联网研究院有限公司 | Harmonic model parameter identification method and device, storage medium and electronic equipment |
| CN109066612A (en) * | 2018-08-27 | 2018-12-21 | 国网湖北省电力有限公司电力科学研究院 | The switching overvoltage calculation method combined based on electromechanical transient and electro-magnetic transient |
| CN109066612B (en) * | 2018-08-27 | 2019-09-27 | 国网湖北省电力有限公司电力科学研究院 | The switching overvoltage calculation method combined based on electromechanical transient and electro-magnetic transient |
| CN109460575A (en) * | 2018-10-09 | 2019-03-12 | 武汉大学 | A kind of hybrid simulation interface system and interpolation algorithm based on Bei Jielong equivalence route |
| CN109460575B (en) * | 2018-10-09 | 2023-10-31 | 武汉大学 | Hybrid simulation interface system based on Bergeron equivalent circuit and interpolation algorithm |
| CN109783895A (en) * | 2018-12-27 | 2019-05-21 | 中国电力科学研究院有限公司 | A kind of electric system hybrid simulation method and system containing distributed generation resource |
| CN109815607A (en) * | 2019-01-31 | 2019-05-28 | 上海交通大学 | Full electromagnetical transient emulation method towards extensive AC network based on time-frequency coordinate transform |
| CN112036002A (en) * | 2020-07-02 | 2020-12-04 | 中国电力科学研究院有限公司 | Electromechanical transient-electromagnetic transient hybrid simulation method and system |
| CN112231981A (en) * | 2020-10-20 | 2021-01-15 | 国网宁夏电力有限公司电力科学研究院 | Method for establishing large-scale electromagnetic transient simulation example |
| CN112231981B (en) * | 2020-10-20 | 2023-04-07 | 国网宁夏电力有限公司电力科学研究院 | Method for establishing large-scale electromagnetic transient simulation example |
| CN113541189A (en) * | 2021-07-19 | 2021-10-22 | 国网上海市电力公司 | Interface positioning method suitable for large-scale alternating current-direct current power grid electromechanical-electromagnetic transient hybrid simulation |
| CN113541189B (en) * | 2021-07-19 | 2024-02-02 | 国网上海市电力公司 | Interface positioning method for electromechanical-electromagnetic transient hybrid simulation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102592003B (en) | 2013-11-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102592003A (en) | Data exchange method for electromechanical and electromagnetic transient hybrid simulation | |
| CN104866665B (en) | Based on the equivalent hybrid simulation method containing power electronic equipment with interacting of interface | |
| Huang et al. | Integrated transmission and distribution system power flow and dynamic simulation using mixed three-sequence/three-phase modeling | |
| CN103559347B (en) | A kind of construction method of extensive AC and DC power system electromagnetic transient simulation model | |
| CN103345162B (en) | A kind of power level digital-analog hybrid simulation system | |
| Huang et al. | A constant-parameter voltage-behind-reactance synchronous machine model based on shifted-frequency analysis | |
| CN102521452B (en) | Computing system of large power grid closed loop | |
| Meliopoulos et al. | Wide area dynamic monitoring and stability controls | |
| CN107093901A (en) | The machine-electricity transient model and emulation mode of a kind of Distributed Power Flow controller | |
| Jia et al. | Fault characteristics and protection adaptability analysis in VSC‐HVDC‐connected offshore wind farm integration system | |
| CN103593566B (en) | The power system comprehensive state method of estimation of mixing quadratic programming form | |
| Farantatos et al. | A predictive out of step protection scheme based on PMU enabled dynamic state estimation | |
| CN104252567B (en) | A kind of modeling method of power transformer three-phase leakage field admittance matrix universal model | |
| CN110516314A (en) | The multi-modal joint simulation method of multi tate for large-scale wind power field | |
| CN102708225A (en) | Fragmentation debugging method for electromagnetic transient simulation of AC/DC large power grid | |
| Sauer | Time-scale features and their applications in electric power system dynamic modeling and analysis | |
| Zhijun et al. | Hybrid simulation of power systems with SVC dynamic phasor model | |
| Zhu et al. | Hybrid-model transient stability simulation using dynamic phasors based HVDC system model | |
| CN104615479A (en) | Equivalent method for electromagnetic and electromechanical transient hybrid simulation electromagnetic side system | |
| CN105914789A (en) | Inverter-interfaced distributed generation device simplified modeling method | |
| Kotian et al. | Dynamic phasor modelling and simulation | |
| Vargas et al. | Two-port network equivalent of VSC-HVDC for power flow studies | |
| Tian et al. | Electromagnetic and electromechanical hybrid simulation based on mosaik framework | |
| Dufour et al. | Testing 750 node distribution grids and devices | |
| Liang et al. | Comparisons of impact on the modeling detail on real time simulation of large power systems with HVDC |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20131113 Termination date: 20201215 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |