CN109766573B - Electromagnetic transient real-time simulation method and device for new energy station - Google Patents

Electromagnetic transient real-time simulation method and device for new energy station Download PDF

Info

Publication number
CN109766573B
CN109766573B CN201811449863.1A CN201811449863A CN109766573B CN 109766573 B CN109766573 B CN 109766573B CN 201811449863 A CN201811449863 A CN 201811449863A CN 109766573 B CN109766573 B CN 109766573B
Authority
CN
China
Prior art keywords
new energy
power generation
generation unit
energy power
access point
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.)
Active
Application number
CN201811449863.1A
Other languages
Chinese (zh)
Other versions
CN109766573A (en
Inventor
刘纯
刘可可
王伟胜
张剑云
李光辉
何国庆
金一丁
孙艳霞
何飞
李丽娜
居秀丽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Gansu Electric Power Co Ltd
Electric Power Research Institute of State Grid Ningxia Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Gansu Electric Power Co Ltd
Electric Power Research Institute of State Grid Ningxia Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI, State Grid Gansu Electric Power Co Ltd, Electric Power Research Institute of State Grid Ningxia Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201811449863.1A priority Critical patent/CN109766573B/en
Publication of CN109766573A publication Critical patent/CN109766573A/en
Application granted granted Critical
Publication of CN109766573B publication Critical patent/CN109766573B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention provides an electromagnetic transient real-time simulation method and device for a new energy station, which are used for calculating the reactance and resistance of a collecting line segment between adjacent new energy power generation units based on an acquired new energy station serial actual model; constructing a new energy station parallel equivalent model based on the new energy power generation unit access point and the reactance and the resistance; and carrying out electromagnetic transient simulation on the new energy station based on the new energy station parallel equivalent model. The invention has the advantages of high simulation speed and small occupied simulation resources, and can realize large-scale simulation of the new energy station. The invention builds the new energy station parallel equivalent model based on the new energy station serial actual model, greatly improves the simulation speed of the new energy station real-time simulation, and achieves the simulation precision and accuracy of the same original serial collection mode.

Description

Electromagnetic transient real-time simulation method and device for new energy station
Technical Field
The invention relates to the technical field of electromagnetic transient simulation, in particular to an electromagnetic transient real-time simulation method and device for a new energy station.
Background
The new energy station is an important component of the power grid, and the electromagnetic transient characteristic simulation of the new energy station is an important research direction in the field of new energy, and has great significance in researching the characteristics of the new energy station and the sending-out characteristics of the new energy station. However, because the new energy station has larger scale and larger difficulty in simulation technology, the new energy station is subjected to electromagnetic transient real-time simulation by adopting real-time simulation. The current common real-time simulation method of the new energy station simply equivalent a new energy station to a power generation unit, and the simple equivalent method ignores the difference of the power generation units in the new energy station and cannot accurately reflect the characteristics of the station. In order to accurately reflect the characteristics of the new energy station, a new energy station actual value model with a certain scale needs to be built.
The actual topological structure of the new energy station is in a series connection mode, a certain number of new energy power generation units are gathered into the system in a series connection mode, if modeling simulation is carried out in a real-time simulation system according to the actual series topological structure, the series connection mode needs to carry out cyclic iteration, limited simulation resources are occupied by iterative computation in a large amount, the simulation speed is very slow, the scale of few power generation units can be simulated only, and the simulation scale is severely limited.
Disclosure of Invention
In order to overcome the defects of low simulation speed, large occupied simulation resources and limited simulation scale in the prior art, the invention provides an electromagnetic transient real-time simulation method and device for a new energy station, and the reactance and resistance of a collecting line segment between adjacent new energy power generation units are calculated based on an acquired new energy station serial actual model; constructing a new energy station parallel equivalent model based on the new energy power generation unit access point and the reactance and the resistance; electromagnetic transient simulation is carried out on the new energy station based on the new energy station parallel equivalent model, so that the simulation speed is high, the occupied simulation resources are small, and the large-scale simulation of the new energy station can be realized.
In order to achieve the above purpose, the present invention adopts the following technical scheme:
in one aspect, the invention provides an electromagnetic transient real-time simulation method of a new energy station, which comprises the following steps:
calculating the reactance and resistance of a collecting line segment between adjacent new energy power generation units based on the obtained new energy station serial actual model;
constructing a new energy station parallel equivalent model based on the new energy power generation unit access point and the reactance and the resistance;
and carrying out electromagnetic transient simulation on the new energy station based on the new energy station parallel equivalent model.
The method for calculating the impedance of the collecting line segment between each new energy power generation unit access point and the grid-connected point based on the obtained new energy station serial actual model comprises the following steps:
and calculating the reactance and resistance of the converging line segments between adjacent new energy power generation units based on the number of segments of the converging line and the resistance and the reactance of the converging line in the new energy station serial actual model.
The construction of the new energy station parallel equivalent model comprises the following steps:
grid-connecting all new energy power generation units based on the access point to obtain a new energy station parallel topological structure;
calculating equivalent longitudinal pressure difference and equivalent transverse pressure difference from each new energy power generation unit access point to the grid-connected point;
calculating equivalent impedance of a collecting line segment between an access point and a grid-connected point of the new energy power generation unit based on the equivalent longitudinal differential pressure and the equivalent transverse differential pressure;
and constructing a new energy station parallel equivalent model based on the parallel topological structure and the equivalent impedance.
The equivalent impedance of the collecting line segment between the access point and the grid-connected point of the new energy power generation unit is as follows:
Z i _B=R i _B+jX i _B
wherein Z is i B is equivalent impedance of a collecting line segment between an i access point and a grid-connected point of the new energy power generation unit, R i B is the equivalent resistance of a collecting line segment between an i access point and a grid-connected point of the new energy power generation unit, X i And B is equivalent reactance of a collecting line segment between an access point and a grid-connected point of the new energy power generation unit i, and j is an imaginary unit.
The R is i _B、X i B is of the formula:
Figure BDA0001886420060000021
Figure BDA0001886420060000022
in Pt (Pt) i For the active power, qt, generated by the new energy power generation unit i i Reactive power generated by the new energy power generation unit i; u (U) 0 The voltage amplitude of the new energy power generation unit grid connection point is set; deltaU iZ_B Is equivalent longitudinal pressure difference delta U of a collecting line segment between an access point and a grid connection point of a new energy power generation unit i iH_B And the equivalent transverse pressure difference of the collecting line segment between the access point of the new energy power generation unit i and the grid-connected point is obtained.
The DeltaU is iZ_B 、ΔU iH_B The formula is as follows:
ΔU iZ_B =ΔU 1Z +ΔU 2Z +···+ΔU iZ
ΔU iH_B =ΔU 1H +ΔU 2H +···+ΔU iH
in the formula DeltaU iZ Is the longitudinal pressure difference delta U of the collecting line segment between the new energy power generation unit i and the new energy power generation unit i-1 iH The transverse pressure difference of the line segment is collected between the new energy power generation unit i and the new energy power generation unit i-1.
The DeltaU is iZ 、ΔU iH The formula is as follows:
Figure BDA0001886420060000031
Figure BDA0001886420060000032
in U 3i-3 Is the voltage of the node behind the access point of the new energy power generation unit i-1, X i Reactance of converging line segment between new energy power generation unit i and new energy power generation unit i-1, R i The resistance of the collecting line segment between the new energy power generation unit i and the new energy power generation unit i-1 is shown; p (P) 3i-3 The active power of the node behind the access point of the new energy power generation unit i-1 is Q 3i-3 And the reactive power of the node behind the access point of the new energy power generation unit i-1.
The P is 3i-3 、Q 3i-3 The formula is as follows:
P 3i-3 =P 3i-2 -I i 2 *R i
Q 3i-3 =Q 3i-2 -I i 2 *X i
wherein I is i The current of the new energy power generation unit i; p (P) 3i-2 Active power of node before access point of new energy power generation unit i, Q 3i-2 The reactive power of the node before the access point of the new energy power generation unit i is used; the P is 3i-2 、Q 3i-2 Determined as follows:
P 3i-2 =P 3i +Pt i
Q 3i-2 =Q 3i +Qt i
wherein P is 3i Active power of node after access point of new energy power generation unit i, Q 3i And the reactive power of the node after the new energy power generation unit i is accessed to the point is determined according to the following formula:
P 3i =P 3i+1 -I i+1 2 *R i+1
Q 3i =Q 3i+1 -I i+1 2 *R i+1
wherein I is i+1 The current of the new energy power generation unit i+1; r is R i+1 The resistor of the collecting line segment between the new energy power generation unit i+1 and the new energy power generation unit i is set; p (P) 3i+1 Active power of the node before the access point of the new energy power generation unit i+1, Q 3i+1 Reactive power of a node before the access point of the new energy power generation unit i+1; the P is 3i+1 、Q 3i+1 Active power P of the previous node of the access point of the new energy power generation unit n-1 respectively 3n-5 And reactive power Q 3n-5 Determining the P 3n-5 、Q 3n-5 Determined as follows:
P 3n-5 =P 3n-3 +Pt n-1
Q 3n-5 =Q 3n-3 +Qt n-1
in Pt (Pt) n-1 Active power Qt generated by the new energy power generation unit n-1 n-1 Reactive power generated by the new energy power generation unit n-1; p (P) 3n-3 Active power of a node behind an access point of the new energy power generation unit n-1 is Q 3n-3 Reactive power of a node behind the access point of the new energy power generation unit n-1; the P is 3n-3 、Q 3n-3 Determined as follows:
P 3n-3 =P 3n-2 -I n 2 *R n
Q 3n-3 =Q 3n-2 -I n 2 *R n
wherein I is n The current of the new energy power generation unit n; r is R n The resistance of the collecting line segment between the new energy power generation unit n and the new energy power generation unit n-1 is shown as nTotal number of new energy power generation units, and n=m; p (P) 3n-2 Active power of node before access point of new energy power generation unit n, Q 3n-2 The reactive power of the node before the access point of the new energy power generation unit n is determined according to the following formula:
P 3n-2 =P 3n-1 =Pt n
Q 3n-2 =Q 3n-1 =Qt n
wherein P is 3n-1 Active power of n access points of new energy power generation unit, Q 3n-1 Reactive power, pt, of n access points of new energy power generation units n For the active power, qt, generated by the new energy power generation unit n n Reactive power generated by the new energy power generation unit n.
After electromagnetic transient simulation is carried out on the new energy station based on the equivalent result, the method comprises the following steps:
when the parallel simulation result is inconsistent with the serial simulation result, correcting the parallel equivalent model of the new energy station;
and the serial simulation result is obtained by performing electromagnetic transient off-line simulation on the new energy station serial actual model.
Correcting the new energy station parallel equivalent model, which comprises the following steps:
determining the reactance and the reactance of the collecting line in unit length based on the model of the collecting line, and calculating the reactance and the resistance of a collecting line segment between adjacent new energy power generation units based on the reactance and the reactance of the collecting line in unit length;
and correcting a new energy station parallel equivalent model based on the new energy power generation unit access point and the reactance and resistance.
On the other hand, the invention provides an electromagnetic transient real-time simulation device of a new energy station, which comprises the following components:
the calculation module is used for calculating the reactance and the resistance of the collecting line segment between the adjacent new energy power generation units based on the obtained new energy station serial actual model;
the modeling module is used for constructing a new energy station parallel equivalent model based on the new energy power generation unit access point and the reactance and resistance;
and the simulation module is used for carrying out electromagnetic transient simulation on the new energy station based on the new energy station parallel equivalent model.
Compared with the closest prior art, the technical scheme provided by the invention has the following beneficial effects:
according to the electromagnetic transient real-time simulation method of the new energy station, the reactance and the resistance of the collecting line segments between adjacent new energy power generation units are calculated based on the obtained new energy station serial actual model; constructing a new energy station parallel equivalent model based on the new energy power generation unit access point and the reactance and the resistance; electromagnetic transient simulation is carried out on the new energy station based on the new energy station parallel equivalent model, so that the simulation speed is high, the occupied simulation resources are small, and the large-scale simulation of the new energy station can be realized;
the electromagnetic transient real-time simulation device of the new energy station comprises a calculation module, a modeling module and a simulation module, wherein the calculation module is used for calculating the reactance and the resistance of a collecting line segment between adjacent new energy power generation units based on an acquired new energy station serial actual model; the modeling module is used for constructing a new energy station parallel equivalent model based on the new energy power generation unit access point and the reactance and the resistance; the simulation module is used for carrying out electromagnetic transient simulation on the new energy station based on the new energy station parallel equivalent model, so that the simulation speed is high, the occupied simulation resources are small, and the large-scale simulation of the new energy station can be realized;
according to the technical scheme provided by the invention, the original serial collection mode is equivalent to the parallel collection mode, namely, the new energy station parallel equivalent model is built based on the new energy station serial actual model, so that the simulation speed of the new energy station real-time simulation is greatly improved, and the simulation precision and accuracy which are the same as those of the original serial collection mode are achieved;
the specific process of constructing the new energy station parallel equivalent model based on the new energy power generation unit access point, the reactance and the resistance is parallel calculation, iterative calculation is not needed, calculation time and resources are greatly saved, and meanwhile simulation precision and accuracy which are the same as those of the new energy station serial actual model are achieved.
Drawings
FIG. 1 is a flow chart of an electromagnetic transient real-time simulation method of a new energy station in an embodiment of the invention;
FIG. 2 is a diagram of a prior art series connection actual model of a new energy station;
fig. 3 is a diagram of a parallel equivalent model structure of a new energy station in an embodiment of the present invention.
Detailed Description
The invention is described in further detail below with reference to the accompanying drawings.
Example 1
The embodiment 1 of the invention provides an electromagnetic transient real-time simulation method of a new energy station, wherein a specific flow chart is shown in fig. 1, and the specific process is as follows:
s101: calculating the reactance and resistance of a collecting line segment between adjacent new energy power generation units based on the obtained new energy station serial actual model;
s102: constructing a new energy station parallel equivalent model (shown in figure 3) based on the new energy power generation unit access point, reactance and resistance;
s103: and carrying out electromagnetic transient simulation on the new energy station based on the new energy station parallel equivalent model.
As shown in fig. 2, the new energy station serial actual model is shown in fig. 2, any new energy power generation unit in the new energy station is represented by i, n new energy power generation units are total, any section of collection line segment between the parallel point of the new energy power generation units and the tail end of the collection line is represented by k, m sections are total, namely, 1.ltoreq.k.ltoreq.m, k=i (namely, line segment k between the i-access point of the new energy power generation unit and the i-1 access point of the new energy power generation unit corresponds to the new energy power generation unit i), n=m (namely, the number of the new energy units is equal to the number of the line segments on the collection line), the new energy power generation unit parallel point is used as a forward direction in the direction of the power grid, the new energy power generation unit parallel point is used as a backward direction in the tail end of the collection line, if i=1, the access point of the new energy power generation unit 1 is used as a node 2 on the collection line, and a node 0 is a grid-connected point, then a node 1 is the access point of the new energy power generation unit 1 is the node (namely, a node 2) is the node.
In S101, the impedance of the collecting line segment between each new energy power generation unit access point and the grid-connected point is calculated based on the obtained new energy station serial actual model, specifically, the reactance and the resistance of the collecting line segment between adjacent new energy power generation units are calculated based on the number of the collecting lines in the new energy station serial actual model and the resistance and the reactance of the collecting lines in unit length. The formula can be expressed as: r is R i =R d *l k ,X i =L d *l k ,X i Reactance of converging line segment between new energy power generation unit i and new energy power generation unit i-1, R i The resistance l of the collecting line segment between the new energy power generation unit i and the new energy power generation unit i-1 k For the length of the kth segment in the collecting line, k=i, and k is more than or equal to 1 and less than or equal to m, m is the number of segments of the collecting line, R d Resistance of the aggregate line per unit length, X d The reactance of the line is assembled for a unit length.
In S102, the construction of the new energy station parallel equivalent model includes:
all new energy power generation units are connected in parallel based on the access point to obtain a new energy station parallel topological structure;
calculating equivalent longitudinal pressure difference and equivalent transverse pressure difference from each new energy power generation unit access point to the grid-connected point;
calculating equivalent impedance of a collecting line segment between an access point and a grid-connected point of the new energy power generation unit based on the equivalent longitudinal differential pressure and the equivalent transverse differential pressure;
and constructing a new energy station parallel equivalent model based on the parallel topology structure and the equivalent impedance.
The equivalent impedance of a collecting line segment between a new energy power generation unit access point and a grid-connected point is as follows:
Z i _B=R i _B+jX i _B
wherein Z is i B is equivalent impedance of a collecting line segment between an i access point and a grid-connected point of the new energy power generation unit, R i B is the connection point of the new energy power generation unit i to the parallel connectionEquivalent resistance, X, of the pooling line segments between dots i And B is equivalent reactance of a collecting line segment between an access point and a grid-connected point of the new energy power generation unit i, and j is an imaginary unit.
R i _B、X i B is of the formula:
Figure BDA0001886420060000071
Figure BDA0001886420060000072
in Pt (Pt) i For the active power, qt, generated by the new energy power generation unit i i Reactive power generated by the new energy power generation unit i; u (U) 0 The voltage amplitude of the new energy power generation unit grid connection point is set; deltaU iZ_B Is equivalent longitudinal pressure difference delta U of a collecting line segment between an access point and a grid connection point of a new energy power generation unit i iH_B And the equivalent transverse pressure difference of the collecting line segment between the access point of the new energy power generation unit i and the grid-connected point is obtained.
ΔU iZ_B 、ΔU iH_B The formula is as follows:
ΔU iZ_B =ΔU 1Z +ΔU 2Z +···+ΔU iZ
ΔU iH_B =ΔU 1H +ΔU 2H +···+ΔU iH
in the formula DeltaU iZ Is the longitudinal pressure difference delta U of the collecting line segment between the new energy power generation unit i and the new energy power generation unit i-1 iH The transverse pressure difference of the line segment is collected between the new energy power generation unit i and the new energy power generation unit i-1.
ΔU iZ 、ΔU iH The formula is as follows:
Figure BDA0001886420060000073
Figure BDA0001886420060000074
in U 3i-3 Is the voltage of the node behind the access point of the new energy power generation unit i-1, P 3i-3 The active power of the node behind the access point of the new energy power generation unit i-1 is Q 3i-3 Reactive power P of the node behind the access point of the new energy power generation unit i-1 3i-3 、Q 3i-3 The formula is as follows:
P 3i-3 =P 3i-2 -I i 2 *R i
Q 3i-3 =Q 3i-2 -I i 2 *X i
wherein I is i The current of the new energy power generation unit i; p (P) 3i-2 Active power of node before access point of new energy power generation unit i, Q 3i-2 The reactive power of the node before the access point of the new energy power generation unit i is used; p (P) 3i-2 、Q 3i-2 Determined as follows:
P 3i-2 =P 3i +Pt i
Q 3i-2 =Q 3i +Qt i
wherein P is 3i Active power of node after access point of new energy power generation unit i, Q 3i And the reactive power of the node after the new energy power generation unit i is accessed to the point is determined according to the following formula:
P 3i =P 3i+1 -I i+1 2 *R i+1
Q 3i =Q 3i+1 -I i+1 2 *R i+1
wherein I is i+1 The current of the new energy power generation unit i+1; r is R i+1 The resistor of the collecting line segment between the new energy power generation unit i+1 and the new energy power generation unit i is set; p (P) 3i+1 Active power of the node before the access point of the new energy power generation unit i+1, Q 3i+1 Reactive power of a node before the access point of the new energy power generation unit i+1; by analogy, P 3i+1 、Q 3i+1 Active power P of the previous node of the access point of the new energy power generation unit n-1 respectively 3n-5 And reactive power Q 3n-5 Determining P 3n-5 、Q 3n-5 Determined as follows:
P 3n-5 =P 3n-3 +Pt n-1
Q 3n-5 =Q 3n-3 +Qt n-1
in Pt (Pt) n-1 Active power Qt generated by the new energy power generation unit n-1 n-1 Reactive power generated by the new energy power generation unit n-1; p (P) 3n-3 Active power of a node behind an access point of the new energy power generation unit n-1 is Q 3n-3 Reactive power of a node behind the access point of the new energy power generation unit n-1; p (P) 3n-3 、Q 3n-3 Determined as follows:
P 3n-3 =P 3n-2 -I n 2 *R n
Q 3n-3 =Q 3n-2 -I n 2 *R n
wherein I is n The current of the new energy power generation unit n; r is R n The resistance of a collecting line segment between the new energy power generation unit n and the new energy power generation unit n-1 is shown, n is the total number of the new energy power generation units, and n=m; p (P) 3n-2 Active power of node before access point of new energy power generation unit n, Q 3n-2 The reactive power of the node before the access point of the new energy power generation unit n is determined according to the following formula:
P 3n-2 =P 3n-1 =Pt n
Q 3n-2 =Q 3n-1 =Qt n
wherein P is 3n-1 Active power of n access points of new energy power generation unit, Q 3n-1 Reactive power, pt, of n access points of new energy power generation units n For the active power, qt, generated by the new energy power generation unit n n Reactive power generated by the new energy power generation unit n.
After the electromagnetic transient simulation is performed on the new energy station based on the equivalent result in S103, the method includes:
when the parallel simulation result is inconsistent with the serial simulation result, correcting the parallel equivalent model of the new energy station;
the series simulation result is obtained by performing electromagnetic transient off-line simulation on the new energy station series actual model.
Correcting the new energy station parallel equivalent model, which comprises the following steps:
determining the reactance and the reactance of the converging line in unit length based on the model of the converging line, and calculating the reactance and the resistance of the converging line segment between adjacent new energy power generation units based on the reactance and the reactance of the converging line in unit length;
and re-determining the equivalent impedance of the collecting line segment between the access point of the new energy power generation unit and the grid-connected point based on the reactance and the resistance.
Example 2
Based on the same inventive concept, embodiment 2 of the present invention further provides an electromagnetic transient real-time simulation device for a new energy station, which includes a calculation module, a modeling module and a simulation module, and the functions of the above modules are described in detail below:
the calculation module is used for calculating the reactance and the resistance of the collecting line segment between the adjacent new energy power generation units based on the obtained new energy station serial actual model;
the modeling module is used for constructing a new energy station parallel equivalent model based on the new energy power generation unit access point, the reactance and the resistance;
and the simulation module is used for carrying out electromagnetic transient simulation on the new energy station based on the new energy station parallel equivalent model.
The calculation module calculates the impedance of a collecting line segment between each new energy power generation unit access point and a grid-connected point based on the obtained new energy station serial actual model, and the calculating module comprises the following steps:
and calculating the reactance and resistance of the converging line segments between adjacent new energy power generation units based on the number of segments of the converging line and the resistance and the reactance of the converging line in the new energy station serial actual model. The formula can be expressed as: r is R i =R d *l k ,X i =L d *l k ,X i Reactance of converging line segment between new energy power generation unit i and new energy power generation unit i-1, R i Generating power sheet for new energyResistance of collecting line segment between element i and new energy power generation unit i-1, l k For the length of the kth segment in the collecting line, k=i, and k is more than or equal to 1 and less than or equal to m, m is the number of segments of the collecting line, R d Resistance of the aggregate line per unit length, X d The reactance of the line is assembled for a unit length.
The modeling module is used for constructing a new energy station parallel equivalent model and comprises the following steps:
the determining unit is used for connecting all the new energy power generation units in a grid mode based on the access point to obtain a new energy station parallel topological structure;
the first calculation unit is used for calculating equivalent longitudinal pressure difference and equivalent transverse pressure difference from each new energy power generation unit access point to the grid-connected point;
the second calculation unit is used for calculating the equivalent impedance of the collecting line segment between the access point and the grid-connected point of the new energy power generation unit based on the equivalent longitudinal differential pressure and the equivalent transverse differential pressure;
the construction unit is used for constructing a new energy station parallel equivalent model based on the parallel topology structure and the equivalent impedance.
The second calculation unit calculates equivalent impedance of a collecting line segment between the access point of the new energy power generation unit and the grid-connected point as follows:
Z i _B=R i _B+jX i _B
wherein Z is i B is equivalent impedance of a collecting line segment between an i access point and a grid-connected point of the new energy power generation unit, R i B is the equivalent resistance of a collecting line segment between an i access point and a grid-connected point of the new energy power generation unit, X i B is equivalent reactance of a collecting line segment between an access point and a grid-connected point of the new energy power generation unit i, and j is an imaginary unit; r is R i _B、X i B is of the formula:
Figure BDA0001886420060000101
Figure BDA0001886420060000102
in Pt (Pt) i For the active power, qt, generated by the new energy power generation unit i i Reactive power generated by the new energy power generation unit i; u (U) 0 The voltage amplitude of the new energy power generation unit grid connection point is set; deltaU iZ_B Is equivalent longitudinal pressure difference delta U of a collecting line segment between an access point and a grid connection point of a new energy power generation unit i iH_B And the equivalent transverse pressure difference of the collecting line segment between the access point of the new energy power generation unit i and the grid-connected point is obtained. DeltaU iZ_B 、ΔU iH_B The formula is as follows:
ΔU iZ_B =ΔU 1Z +ΔU 2Z +···+ΔU iZ
ΔU iH_B =ΔU 1H +ΔU 2H +···+ΔU iH
in the formula DeltaU iZ Is the longitudinal pressure difference delta U of the collecting line segment between the new energy power generation unit i and the new energy power generation unit i-1 iH The transverse pressure difference of the line segment is collected between the new energy power generation unit i and the new energy power generation unit i-1. DeltaU iZ 、ΔU iH The formula is as follows:
Figure BDA0001886420060000111
Figure BDA0001886420060000112
in U 3i-3 Is the voltage of the node behind the access point of the new energy power generation unit i-1, X i Reactance of converging line segment between new energy power generation unit i and new energy power generation unit i-1, R i The resistance of the collecting line segment between the new energy power generation unit i and the new energy power generation unit i-1 is shown; p (P) 3i-3 The active power of the node behind the access point of the new energy power generation unit i-1 is Q 3i-3 And the reactive power of the node behind the access point of the new energy power generation unit i-1.
P 3i-3 、Q 3i-3 The formula is as follows:
P 3i-3 =P 3i-2 -I i 2 *R i
Q 3i-3 =Q 3i-2 -I i 2 *X i
wherein I is i The current of the new energy power generation unit i; p (P) 3i-2 Active power of node before access point of new energy power generation unit i, Q 3i-2 The reactive power of the node before the access point of the new energy power generation unit i is used; p (P) 3i-2 、Q 3i-2 Determined as follows:
P 3i-2 =P 3i +Pt i
Q 3i-2 =Q 3i +Qt i
wherein P is 3i Active power of node after access point of new energy power generation unit i, Q 3i And the reactive power of the node after the new energy power generation unit i is accessed to the point is determined according to the following formula:
P 3i =P 3i+1 -I i+1 2 *R i+1
Q 3i =Q 3i+1 -I i+1 2 *R i+1
wherein I is i+1 The current of the new energy power generation unit i+1; r is R i+1 The resistor of the collecting line segment between the new energy power generation unit i+1 and the new energy power generation unit i is set; p (P) 3i+1 Active power of the node before the access point of the new energy power generation unit i+1, Q 3i+1 Reactive power of a node before the access point of the new energy power generation unit i+1; p (P) 3i+1 、Q 3i+1 Active power P of the previous node of the access point of the new energy power generation unit n-1 respectively 3n-5 And reactive power Q 3n-5 Determining P 3n-5 、Q 3n-5 Determined as follows:
P 3n-5 =P 3n-3 +Pt n-1
Q 3n-5 =Q 3n-3 +Qt n-1
in Pt (Pt) n-1 Active power Qt generated by the new energy power generation unit n-1 n-1 Reactive power generated by the new energy power generation unit n-1; p (P) 3n-3 The node behind the access point of the new energy power generation unit n-1 is provided withPower, Q 3n-3 Reactive power of a node behind the access point of the new energy power generation unit n-1; p (P) 3n-3 、Q 3n-3 Determined as follows:
P 3n-3 =P 3n-2 -I n 2 *R n
Q 3n-3 =Q 3n-2 -I n 2 *R n
wherein I is n The current of the new energy power generation unit n; r is R n The resistance of a collecting line segment between the new energy power generation unit n and the new energy power generation unit n-1 is shown, n is the total number of the new energy power generation units, and n=m; p (P) 3n-2 Active power of node before access point of new energy power generation unit n, Q 3n-2 The reactive power of the node before the access point of the new energy power generation unit n is determined according to the following formula:
P 3n-2 =P 3n-1 =Pt n
Q 3n-2 =Q 3n-1 =Qt n
wherein P is 3n-1 Active power of n access points of new energy power generation unit, Q 3n-1 Reactive power, pt, of n access points of new energy power generation units n For the active power, qt, generated by the new energy power generation unit n n Reactive power generated by the new energy power generation unit n.
The device provided by the embodiment 2 of the invention further comprises a judging module, wherein after the electromagnetic transient simulation is carried out on the new energy station based on the equivalent result, the judging module judges whether the parallel simulation result obtained by the electromagnetic transient simulation is consistent with the serial simulation result, when the parallel simulation result is inconsistent with the serial simulation result, the reactance and the reactance of the collecting line in unit length are determined again based on the model of the collecting line, and the reactance and the resistance of the collecting line segment between adjacent new energy power generation units are calculated based on the reactance and the reactance of the collecting line in unit length; and correcting the parallel equivalent model of the new energy station based on the access point, reactance and resistance of the new energy power generation unit.
For convenience of description, the parts of the above apparatus are described as being functionally divided into various modules or units, respectively. Of course, the functions of each module or unit may be implemented in the same piece or pieces of software or hardware when implementing the present application.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention and not for limiting the same, and a person skilled in the art may still make modifications and equivalents to the specific embodiments of the present invention with reference to the above embodiments, and any modifications and equivalents not departing from the spirit and scope of the present invention are within the scope of the claims of the present invention as filed herewith.

Claims (10)

1. The electromagnetic transient real-time simulation method of the new energy station is characterized by comprising the following steps of:
calculating the reactance and resistance of a collecting line segment between adjacent new energy power generation units based on the obtained new energy station serial actual model;
constructing a new energy station parallel equivalent model based on the new energy power generation unit access point and the reactance and the resistance;
calculating the reactance and resistance of a converging line segment between adjacent new energy power generation units based on the number of segments of the converging line and the resistance and the reactance of the converging line in the new energy station serial actual model; grid-connecting all new energy power generation units based on the access point to obtain a new energy station parallel topological structure; calculating equivalent longitudinal pressure difference and equivalent transverse pressure difference from each new energy power generation unit access point to the grid-connected point; calculating equivalent impedance of a collecting line segment between an access point and a grid-connected point of the new energy power generation unit based on the equivalent longitudinal differential pressure and the equivalent transverse differential pressure;
performing electromagnetic transient simulation on the new energy station based on the new energy station parallel equivalent model;
the method for calculating the impedance of the collecting line segment between each new energy power generation unit access point and the grid-connected point based on the obtained new energy station serial actual model comprises the following steps:
and calculating the reactance and resistance of the converging line segments between adjacent new energy power generation units based on the number of segments of the converging line and the resistance and the reactance of the converging line in the new energy station serial actual model.
2. The electromagnetic transient real-time simulation method of the new energy station according to claim 1, wherein the construction of the new energy station parallel equivalent model comprises the following steps:
grid-connecting all new energy power generation units based on the access point to obtain a new energy station parallel topological structure;
calculating equivalent longitudinal pressure difference and equivalent transverse pressure difference from each new energy power generation unit access point to the grid-connected point;
calculating equivalent impedance of a collecting line segment between an access point and a grid-connected point of the new energy power generation unit based on the equivalent longitudinal differential pressure and the equivalent transverse differential pressure;
and constructing a new energy station parallel equivalent model based on the parallel topological structure and the equivalent impedance.
3. The electromagnetic transient real-time simulation method of the new energy station according to claim 2, wherein the equivalent impedance of a collecting line segment between an access point and a grid-connected point of the new energy power generation unit is as follows:
Z i _B=R i _B+jX i _B
wherein Z is i B is equivalent impedance of a collecting line segment between an i access point and a grid-connected point of the new energy power generation unit, R i B is the equivalent resistance of a collecting line segment between an i access point and a grid-connected point of the new energy power generation unit, X i And B is equivalent reactance of a collecting line segment between an access point and a grid-connected point of the new energy power generation unit i, and j is an imaginary unit.
4. The electromagnetic transient real-time simulation method of the new energy station according to claim 3, wherein the R is as follows i _B、X i B is of the formula:
Figure FDA0004022974540000021
Figure FDA0004022974540000022
in Pt (Pt) i For the active power, qt, generated by the new energy power generation unit i i Reactive power generated by the new energy power generation unit i; u (U) 0 The voltage amplitude of the new energy power generation unit grid connection point is set; deltaU iZ_B Is equivalent longitudinal pressure difference delta U of a collecting line segment between an access point and a grid connection point of a new energy power generation unit i iH_B And the equivalent transverse pressure difference of the collecting line segment between the access point of the new energy power generation unit i and the grid-connected point is obtained.
5. The electromagnetic transient real-time simulation method of the new energy station according to claim 4, wherein the Δu is iZ_B 、ΔU iH_B The formula is as follows:
ΔU iZ_B =ΔU 1Z +ΔU 2Z +…+ΔU iZ
ΔU iH_B =ΔU 1H +ΔU 2H +…+ΔU iH
in the formula DeltaU iZ Is the longitudinal pressure difference delta U of the collecting line segment between the new energy power generation unit i and the new energy power generation unit i-1 iH The transverse pressure difference of the line segment is collected between the new energy power generation unit i and the new energy power generation unit i-1.
6. The electromagnetic transient real-time simulation method of the new energy station according to claim 5, wherein the Δu is iZ 、ΔU iH The formula is as follows:
Figure FDA0004022974540000023
Figure FDA0004022974540000024
in U 3i-3 Is the voltage of the node behind the access point of the new energy power generation unit i-1, X i Reactance of converging line segment between new energy power generation unit i and new energy power generation unit i-1, R i The resistance of the collecting line segment between the new energy power generation unit i and the new energy power generation unit i-1 is shown; p (P) 3i-3 The active power of the node behind the access point of the new energy power generation unit i-1 is Q 3i-3 And the reactive power of the node behind the access point of the new energy power generation unit i-1.
7. The electromagnetic transient real-time simulation method of the new energy station according to claim 6, wherein the P is as follows 3i-3 、Q 3i-3 The formula is as follows:
P 3i-3 =P 3i-2 -I i 2 *R i
Q 3i-3 =Q 3i-2 -I i 2 *X i
wherein I is i The current of the new energy power generation unit i; p (P) 3i-2 Active power of node before access point of new energy power generation unit i, Q 3i-2 The reactive power of the node before the access point of the new energy power generation unit i is used; the P is 3i-2 、Q 3i-2 Determined as follows:
P 3i-2 =P 3i +Pt i
Q 3i-2 =Q 3i +Qt i
wherein P is 3i Active power of node after access point of new energy power generation unit i, Q 3i And the reactive power of the node after the new energy power generation unit i is accessed to the point is determined according to the following formula:
P 3i =P 3i+1 -I i+1 2 *R i+1
Q 3i =Q 3i+1 -I i+1 2 *R i+1
wherein I is i+1 The current of the new energy power generation unit i+1; r is R i+1 The resistor of the collecting line segment between the new energy power generation unit i+1 and the new energy power generation unit i is set; p (P) 3i+1 Before the new energy power generation unit i+1 is connected with the pointActive power of a node, Q 3i+1 Reactive power of a node before the access point of the new energy power generation unit i+1; the P is 3i+1 、Q 3i+1 Active power P of the previous node of the access point of the new energy power generation unit n-1 respectively 3n-5 And reactive power Q 3n-5 Determining the P 3n-5 、Q 3n-5 Determined as follows:
P 3n-5 =P 3n-3 +Pt n-1
Q 3n-5 =Q 3n-3 +Qt n-1
in Pt (Pt) n-1 Active power Qt generated by the new energy power generation unit n-1 n-1 Reactive power generated by the new energy power generation unit n-1; p (P) 3n-3 Active power of a node behind an access point of the new energy power generation unit n-1 is Q 3n-3 Reactive power of a node behind the access point of the new energy power generation unit n-1; the P is 3n-3 、Q 3n-3 Determined as follows:
P 3n-3 =P 3n-2 -I n 2 *R n
Q 3n-3 =Q 3n-2 -I n 2 *R n
wherein I is n The current of the new energy power generation unit n; r is R n The resistance of a collecting line segment between the new energy power generation unit n and the new energy power generation unit n-1 is shown, n is the total number of the new energy power generation units, and n=m; p (P) 3n-2 Active power of node before access point of new energy power generation unit n, Q 3n-2 The reactive power of the node before the access point of the new energy power generation unit n is determined according to the following formula:
P 3n-2 =P 3n-1 =Pt n
Q 3n-2 =Q 3n-1 =Qt n
wherein P is 3n-1 Active power of n access points of new energy power generation unit, Q 3n-1 Reactive power, pt, of n access points of new energy power generation units n For the active power, qt, generated by the new energy power generation unit n n Reactive power generated by the new energy power generation unit n.
8. The electromagnetic transient real-time simulation method of the new energy station according to claim 1, wherein after electromagnetic transient simulation of the new energy station based on the equivalent result, comprising:
when the parallel simulation result is inconsistent with the serial simulation result, correcting the parallel equivalent model of the new energy station;
and the serial simulation result is obtained by performing electromagnetic transient off-line simulation on the new energy station serial actual model.
9. The electromagnetic transient real-time simulation method of the new energy station according to claim 8, wherein the correcting the new energy station parallel equivalent model comprises:
determining the reactance and the resistance of the collecting line in unit length based on the model of the collecting line, and calculating the reactance and the resistance of a collecting line segment between adjacent new energy power generation units based on the reactance and the resistance of the collecting line in unit length;
and correcting a new energy station parallel equivalent model based on the new energy power generation unit access point and the reactance and resistance.
10. An electromagnetic transient real-time simulation device of a new energy station is characterized by comprising:
the calculation module is used for calculating the reactance and the resistance of the collecting line segment between the adjacent new energy power generation units based on the obtained new energy station serial actual model;
the modeling module is used for constructing a new energy station parallel equivalent model based on the new energy power generation unit access point and the reactance and resistance;
calculating the reactance and resistance of a converging line segment between adjacent new energy power generation units based on the number of segments of the converging line and the resistance and the reactance of the converging line in the new energy station serial actual model; grid-connecting all new energy power generation units based on the access point to obtain a new energy station parallel topological structure; calculating equivalent longitudinal pressure difference and equivalent transverse pressure difference from each new energy power generation unit access point to the grid-connected point; calculating equivalent impedance of a collecting line segment between an access point and a grid-connected point of the new energy power generation unit based on the equivalent longitudinal differential pressure and the equivalent transverse differential pressure;
and the simulation module is used for carrying out electromagnetic transient simulation on the new energy station based on the new energy station parallel equivalent model.
CN201811449863.1A 2018-11-30 2018-11-30 Electromagnetic transient real-time simulation method and device for new energy station Active CN109766573B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811449863.1A CN109766573B (en) 2018-11-30 2018-11-30 Electromagnetic transient real-time simulation method and device for new energy station

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811449863.1A CN109766573B (en) 2018-11-30 2018-11-30 Electromagnetic transient real-time simulation method and device for new energy station

Publications (2)

Publication Number Publication Date
CN109766573A CN109766573A (en) 2019-05-17
CN109766573B true CN109766573B (en) 2023-04-21

Family

ID=66451179

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811449863.1A Active CN109766573B (en) 2018-11-30 2018-11-30 Electromagnetic transient real-time simulation method and device for new energy station

Country Status (1)

Country Link
CN (1) CN109766573B (en)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103730910B (en) * 2013-12-20 2016-02-24 国电南瑞科技股份有限公司 The Dynamic Equivalence that a kind of large-scale photovoltaic power station is grid-connected
CN104318088B (en) * 2014-10-11 2017-07-14 清华大学 A kind of electromagnetic transient in power system emulation mode containing many electronic power switches
CN106886616B (en) * 2015-12-15 2020-06-16 中国电力科学研究院 Automatic network division method for large-scale electromagnetic transient power grid simulation
CN106099891B (en) * 2016-07-26 2018-06-22 华南理工大学 The analysis method that shunt reactor is distributed rationally is carried out to marine wind electric field leading-out terminal extra large cable

Also Published As

Publication number Publication date
CN109766573A (en) 2019-05-17

Similar Documents

Publication Publication Date Title
CN103810646B (en) Improved projection integral algorithm based active power distribution system dynamic simulation method
CN104079003B (en) Containing the probability load flow calculation method of photo-voltaic power supply distribution network
CN107529644B (en) Linear approximation method for static voltage stability domain boundary of power system
CN110596530A (en) Low-current ground fault line selection method
CN106886839B (en) Hybrid integer programming-based water-fire-electricity generator set combination optimization scheduling method
CN103336866A (en) Processing method for negative resistance containing branch in electromagnetic transient simulation
CN104766142B (en) Transient Stability Constraints optimal load flow calculation method based on EEAC and trace sensitivity
CN103678900A (en) Network decoupling calculation method used for regional power distribution network real-time simulation
CN111900718B (en) Active power distribution network dynamic simulation method based on multi-stage optimization catch-up variational iteration method
CN103632046A (en) Power grid load flow calculation method
CN110765594B (en) Method and device for controlling internal interaction interface of wind power plant real-time simulator
CN115079592A (en) Pipe network simulation method for thermodynamic system of ship nuclear power device
CN109766573B (en) Electromagnetic transient real-time simulation method and device for new energy station
CN112287605B (en) Power flow checking method based on graph convolution network acceleration
CN106951089A (en) Gesture interaction method and system
CN112234598B (en) Electromagnetic transient simulation initialization method
CN114188945B (en) Method and device for calculating short-circuit current of power distribution network containing photovoltaic power supply
CN103199524A (en) Load flow calculation method adapting to multiple distributed power access
CN106951646A (en) Multi tate interface method and device in a kind of power system real-time simulation
CN116306236A (en) Wind power plant dynamic equivalent modeling method and system based on GRU-LSTM-FC combined network
CN107425519B (en) Method for calculating maximum power supply capacity of three-phase power distribution network containing distributed power supply
CN107657071A (en) Based on the power system uncertainty time-domain simulation method for improving sparse probability assignments method
CN111082414B (en) Transient voltage calculation method and system
CN113659559A (en) Method for transient stability analysis of DFIG grid-connected power system
CN117973096B (en) Direct-drive wind farm parallel simulation method, system, terminal and medium based on double-layer network division

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20200628

Address after: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15

Applicant after: CHINA ELECTRIC POWER RESEARCH INSTITUTE Co.,Ltd.

Applicant after: STATE GRID CORPORATION OF CHINA

Applicant after: ELECTRIC POWER RESEARCH INSTITUTE OF STATE GRID NINGXIA ELECTRIC POWER Co.,Ltd.

Applicant after: STATE GRID GANSU ELECTRIC POWER Co.

Address before: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15

Applicant before: CHINA ELECTRIC POWER RESEARCH INSTITUTE Co.,Ltd.

Applicant before: STATE GRID CORPORATION OF CHINA

Applicant before: ELECTRIC POWER RESEARCH INSTITUTE OF STATE GRID NINGXIA ELECTRIC POWER Co.,Ltd.

SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant