CN111259565A - Dynamic simulation method and system for voltage source type current converter - Google Patents
Dynamic simulation method and system for voltage source type current converter Download PDFInfo
- Publication number
- CN111259565A CN111259565A CN202010084659.5A CN202010084659A CN111259565A CN 111259565 A CN111259565 A CN 111259565A CN 202010084659 A CN202010084659 A CN 202010084659A CN 111259565 A CN111259565 A CN 111259565A
- Authority
- CN
- China
- Prior art keywords
- voltage source
- simulation
- source type
- converter
- dynamic
- 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
- 238000000034 method Methods 0.000 title claims abstract description 66
- 238000005094 computer simulation Methods 0.000 title claims abstract description 32
- 238000004088 simulation Methods 0.000 claims abstract description 107
- 230000001052 transient effect Effects 0.000 claims abstract description 29
- 238000000354 decomposition reaction Methods 0.000 claims description 29
- 230000003044 adaptive effect Effects 0.000 claims description 12
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000006798 recombination Effects 0.000 description 3
- 238000005215 recombination Methods 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention provides a dynamic simulation method and a dynamic simulation system for a voltage source type current converter. The method is based on the principle of a multi-band dynamic vector method, and the dynamic vector model of the voltage source type current converter is decomposed, recombined and frequency-shifted to obtain the multi-band dynamic vector model; the simulation time is decomposed into a plurality of time subintervals according to the operating condition, and the simulation model and the window width of the windowed Fourier are selected for self-adaptive simulation according to the operating condition of the voltage source type converter in each time subinterval, so that different simulation models are adopted for self-adaptation of the steady-state operating condition and the fault operating condition. The invention combines the multi-band dynamic vector model with the self-adaptive variable simulation step length, solves the problems that the simulation scale and the simulation speed of the traditional VSC electromagnetic transient model are limited and the simulation precision cannot be ensured by the traditional VSC dynamic phasor model and a frequency offset method, and realizes the high-efficiency simulation which considers the steady-state working condition and the fault working condition and adopts the multi-band dynamic phasor model with different window widths.
Description
Technical Field
The invention relates to the technical field of converter performance research, in particular to a dynamic simulation method and system of a voltage source type converter.
Background
The large-scale access of new energy power generation and the wide construction of medium and low Voltage direct current transmission/distribution networks cause the operation mode of a power system to be deeply changed, and the transient characteristics of a grid-connected Voltage Source Converter (VSC) directly influence the dynamic characteristics of the power system containing the large-scale new energy power generation. The VSC has complex action characteristics, and the converter simulation model which can establish both precision and efficiency is always the focus content of academic research. Nowadays, a converter simulation model suitable for various occasions is established, such as a VSC detailed model, and the sampling theorem is used as a basis to capture the precise change of a high-frequency switch through an extremely small step length, but on one hand, the calculated amount is very large due to the excessively small step length, and on the other hand, the admittance array of the simulated power system must be updated during each switching action, so that the simulation burden is increased. The VSC dynamic vector model changes a high-frequency simulation object into a low-frequency simulation object through frequency shift, can perform simulation with a large step, but cannot consider harmonic waves, so that the precision of the model is reduced to a certain extent. VSC multifrequency section dynamic vector model, it is fixed to decompose the window width in windowing Fourier decomposition process, is power frequency cycle usually, consequently can not carry out step length self-adaptation according to the operating mode demand of difference, so how to realize giving consideration to the dynamic simulation of voltage source transverter of simulation speed and simulation accuracy becomes a technical problem that awaits the solution urgently. The method is used for solving the problem that the simulation scale and the simulation speed of the traditional VSC electromagnetic transient model are limited.
Disclosure of Invention
The invention aims to provide a dynamic simulation method and a dynamic simulation system for a voltage source type converter, which aim to solve the problems that the simulation scale and the simulation speed of the traditional VSC electromagnetic transient model are limited and the simulation precision cannot be ensured by the traditional VSC dynamic phasor model and a frequency offset method.
In order to achieve the purpose, the invention provides the following scheme:
a dynamic simulation method of a voltage source type converter comprises the following steps:
establishing a dynamic vector model of the voltage source type current converter;
decomposing, recombining and frequency shifting the dynamic vector model based on a multi-band dynamic vector method principle to obtain a multi-band dynamic vector model;
decomposing the simulation time into a plurality of time subintervals according to the operating condition;
based on the self-adaptive strategy principle, the window width of the multi-band dynamic phasor model is self-adaptively switched according to the working condition of the voltage type converter in each time subinterval, wherein the working condition is a steady-state working condition or a fault working condition;
and simulating each time subinterval by adopting the window width after each time subinterval is switched to realize the dynamic simulation of the voltage source type current converter.
Optionally, the establishing a dynamic vector model of the voltage source converter specifically includes:
the dynamic vector model of the voltage source type current converter is established as follows:
wherein t represents time, ua、ub、ucRespectively representing three-phase voltages at an outlet side of the voltage source type converter; i.e. ia、ib、icRespectively representing three-phase currents on the inlet side of the voltage source type converter; u. ofdThe voltage of the direct current side of the voltage source type converter is represented, and R and L respectively represent the equivalent resistance and the equivalent inductance of the coupling transformer at the alternating current side of the voltage source type converter; < > represents a dynamic vector; sa、Sb、ScRespectively representing the switching functions of a three-phase full bridge arm of the voltage source type converter;
optionally, the decomposing, recombining, and frequency shifting the dynamic vector model based on the principle of the multi-band dynamic vector method to obtain the multi-band dynamic vector model specifically includes:
based on the principle of a multi-band dynamic vector method, decomposing, recombining and frequency shifting the dynamic vector model to obtain the multi-band dynamic vector model is as follows:
wherein,
respectively representing the sub-band dynamic phasors of the three-phase current signals at the inlet side of the voltage source type converter,
respectively representing the sub-band dynamic phasors of the three-phase voltage signal at the outlet side of the voltage source type converter,
and
respectively representing sub-frequency band dynamic phasors of switching function signals of a three-phase full bridge arm of the voltage source type converter;
wherein x (t) represents ia、ib、ic、ua、ub、uc、Sa、SbOr Sc；fnFrequency of the nth frequency band signal, BnRepresenting the nth frequency band signal, Hn-1And HnRespectively representing the last harmonic order, X, of the (n-1) th and nth frequency band signalsh(t) denotes the h-Fourier coefficient of x (t), ωhDenotes the angular frequency, ω, of the h-th harmonic of x (t)sAt the fundamental angular frequency, f (T)i) Indicating the window when performing a windowed Fourier decomposition of the state vector for the ith time subintervalWide, x (τ) denotes the signal of x (t) over a period τ.
Optionally, decomposing the simulation time into a plurality of time subintervals according to the operating condition specifically includes:
decomposing the simulation time T into a plurality of time subintervals according to the operation condition as follows: t ═ T1∪T2∪T3∪……∪TI；T1、T2、T3And TIRespectively, the 1 st, 2 nd, 3 rd and I th time subintervals are indicated.
Optionally, based on the adaptive strategy principle, the window width of the multi-band dynamic phasor model is adaptively switched according to the working condition of the voltage-type converter in each time subinterval, where the working condition is a steady-state working condition or a fault working condition, and the method specifically includes:
when the voltage source type converter is in a steady-state working condition in the ith time subinterval, switching the window width of the windowed Fourier to a power frequency period; wherein, I is 1,2,3 …, and I represents the number of sub-time intervals;
and when the voltage source type converter is in a transient working condition in the ith time subinterval, switching the window width of the windowed Fourier to a switching period.
The invention also provides a dynamic simulation system of the voltage source type converter, which specifically comprises:
the dynamic vector model establishing module is used for establishing a dynamic vector model of the voltage source type current converter;
the multi-band dynamic vector model acquisition module is used for decomposing, recombining and frequency shifting the dynamic vector model based on a multi-band dynamic vector method principle to obtain a multi-band dynamic vector model;
the simulation time decomposition module is used for decomposing the simulation time into a plurality of time subintervals according to the operation working condition;
the adaptive window width switching module is used for adaptively switching the window width of the multi-band dynamic phasor model according to the working condition of the voltage type converter in each time subinterval based on an adaptive strategy principle, wherein the working condition is a steady-state working condition or a fault working condition;
and the simulation module is used for simulating each time subinterval by adopting the window width after each time subinterval is switched, so as to realize the dynamic simulation of the voltage source type current converter.
Optionally, the dynamic vector model establishing module specifically includes:
the dynamic vector model establishing submodule is used for establishing a dynamic vector model of the voltage source type current converter as follows:
wherein t represents time, ua、ub、ucRespectively representing three-phase voltages at an outlet side of the voltage source type converter; i.e. ia、ib、icRespectively representing three-phase currents on the inlet side of the voltage source type converter; u. ofdThe voltage of the direct current side of the voltage source type converter is represented, and R and L respectively represent the equivalent resistance and the equivalent inductance of the coupling transformer at the alternating current side of the voltage source type converter; < > represents a dynamic vector; sa、Sb、ScRespectively representing the switching functions of a three-phase full bridge arm of the voltage source type converter;
optionally, the multi-band dynamic vector model obtaining module specifically includes:
the multi-band dynamic vector model obtaining submodule is used for decomposing, recombining and frequency shifting the dynamic vector model based on the principle of a multi-band dynamic vector method to obtain the multi-band dynamic vector model, and comprises the following steps:
wherein,
respectively representing the sub-band dynamic phasors of the three-phase current signals at the inlet side of the voltage source type converter,
respectively representing the sub-band dynamic phasors of the three-phase voltage signal at the outlet side of the voltage source type converter,
respectively representing sub-frequency band dynamic phasors of switching function signals of a three-phase full bridge arm of the voltage source type converter;
wherein x (t) represents ia、ib、ic、ua、ub、uc、Sa、SbOr Sc；fnFrequency of the nth frequency band signal, BnRepresenting the nth frequency band signal, Hn-1And HnRespectively representing the last harmonic order, X, of the (n-1) th and nth frequency band signalsh(t) denotes the h-Fourier coefficient of x (t), ωhDenotes the angular frequency, ω, of the h-th harmonic of x (t)sAt the fundamental angular frequency, f (T)i) The window width when performing windowed fourier decomposition on the state vector of the ith time subinterval is shown, and x (τ) represents the signal over a period τ of x (t).
Optionally, the simulation time decomposition module specifically includes:
the simulation time decomposition submodule is used for decomposing the simulation time T into a plurality of time subintervals according to the operation condition, and the time subintervals are as follows: t ═ T1∪T2∪T3∪……∪TI；T1、T2、T3And TIRespectively, the 1 st, 2 nd, 3 rd and I th time subintervals are indicated.
Optionally, the adaptive window width switching module specifically includes:
the first switching submodule is used for switching the window width of the windowed Fourier to a power frequency period when the voltage source type converter is in a steady-state working condition in the ith time subinterval; wherein, I is 1,2,3 …, and I represents the number of sub-time intervals;
and the second switching submodule is used for switching the window width of the windowed Fourier to a switching period when the voltage source type converter is in a transient working condition in the ith time subinterval.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the invention provides a dynamic simulation method and a dynamic simulation system for a voltage source type current converter. Firstly, establishing a dynamic vector model of a voltage source type converter; decomposing, recombining and frequency shifting the dynamic vector model based on a multi-band dynamic vector method principle to obtain a multi-band dynamic vector model; then, the simulation time is decomposed into a plurality of time subintervals according to the operation working condition, and the windowing Fourier decomposition window width of the multi-band dynamic phasor is adaptively switched to perform simulation according to the working condition of the voltage source type current converter in each time subinterval based on the adaptive strategy principle. The invention combines a multi-band dynamic vector model with a self-adaptive variable simulation step length, solves the problems that the simulation scale and the simulation speed of the traditional VSC electromagnetic transient model are limited and the simulation precision cannot be ensured by the traditional VSC dynamic vector model and a frequency offset method, and realizes the dynamic simulation of the voltage source type converter which gives consideration to both the simulation speed and the simulation precision.
The method is based on the self-adaptive strategy principle, and the windowed Fourier window width is self-adaptively switched according to the steady-state working condition or the fault working condition of the voltage source type converter in each time subinterval, so that a simulation model giving consideration to the steady-state working condition and the fault working condition is adopted, and the dynamic simulation of the voltage source type converter giving consideration to the simulation speed and the simulation precision is realized; and the simulation is started after the corresponding window width is adopted after the switching is carried out on the specific working condition, so that the simulation precision is ensured and the simulation speed is improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
Fig. 1 is a flowchart of a dynamic simulation method of a voltage source converter according to the present invention;
fig. 2 is a topology structure diagram of the voltage source converter provided by the present invention;
FIG. 3 is a schematic diagram of a multi-band dynamic vector method provided by the present invention;
FIG. 4 is a schematic diagram of the band segmentation principle provided by the present invention;
fig. 5 is a schematic diagram illustrating the principle of the adaptive strategy provided by the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide a dynamic simulation method and a dynamic simulation system for a voltage source type converter, which aim to solve the problems that the simulation scale and the simulation speed of the traditional VSC electromagnetic transient model are limited and the simulation precision cannot be ensured by the traditional VSC dynamic vector model and a frequency migration method.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
In order to achieve the above object, the present invention provides a dynamic simulation method of a voltage source converter, as shown in fig. 1, the simulation method includes the following steps:
the topological structure of the voltage source type converter is shown in fig. 2, and the dynamic vector model of the voltage source type converter established by the invention is as follows:
wherein t represents time, ua、ub、ucRespectively representing three-phase voltages at an outlet side of the voltage source type converter; i.e. ia、ib、icRespectively representing three-phase currents on the inlet side of the voltage source type converter; u. ofdThe voltage of the direct current side of the voltage source type converter is represented, and R and L respectively represent the equivalent resistance and the equivalent inductance of the coupling transformer at the alternating current side of the voltage source type converter; < > represents a dynamic vector; sa、Sb、ScRespectively representing the switching functions of a three-phase full bridge arm of the voltage source type converter;
and 102, decomposing, recombining and frequency shifting the dynamic vector model based on the principle of a multi-band dynamic vector method to obtain the multi-band dynamic vector model.
Based on the principle of a multi-band dynamic vector method, the invention decomposes, recombines and shifts the frequency of the dynamic vector model to obtain the multi-band dynamic vector model, which comprises the following steps:
wherein,
respectively representing the sub-band dynamic phasors of the three-phase current signals at the inlet side of the voltage source type converter,
individual watchShows the sub-band dynamic phasor of the three-phase voltage signal at the outlet side of the voltage source type converter,
and
respectively representing sub-frequency band dynamic phasors of switching function signals of a three-phase full bridge arm of the voltage source type converter;
wherein x (t) represents ia、ib、ic、ua、ub、uc、Sa、SbOr Sc；fnFrequency of the nth frequency band signal, BnRepresenting the nth frequency band signal, Hn-1And HnRespectively representing the last harmonic order, X, of the (n-1) th and nth frequency band signalsh(t) denotes the h-Fourier coefficient of x (t), ωhDenotes the angular frequency, ω, of the h-th harmonic of x (t)sAt the fundamental angular frequency, f (T)i) The window width when performing windowed fourier decomposition on the state vector of the ith time subinterval is shown, and x (τ) represents the signal over a period τ of x (t).
The principle of the multi-band dynamic vector method is shown in fig. 3, and includes the steps of decomposition, recombination and frequency shift.
The specific steps of signal decomposition and recombination are as follows: voltage, current, etc. in an electrical power system can be seen as periodically varying electrical signals. For period of T0In a period tau e (T-T) of a period x (tau)0,t]The fourier decomposition of its complex form is:
in the formula: x (t) is a complex form of the electrical signal x (τ), ωs＝2π/T0At fundamental angular frequency, Xh(t) is the h-th order Fourier coefficient, i.e. the h-order "dynamic phasor"。
The fourier coefficient h in the equation (1) is theoretically infinite, but in the electromagnetic transient simulation, the frequency corresponding to the simulation step length is generally 10 times of the signal frequency according to the requirements of precision and sampling theorem. Therefore, in electromagnetic transient simulation, the harmonic number h is generally finite, and the maximum value can be determined according to the simulation step size (for example, the maximum harmonic number h corresponding to the 50 μ s simulation step size is 40). At this time, formula (1) can be written as:
the mathematical meaning of the fourier transform is: any periodic signal satisfying the dirichlet condition can be represented as a group by ejhωstIs a linear combination of orthogonal bases. If the signal frequency is divided into N sub-bands according to fig. 4, by the linear combination law, equation (2) can be grouped and recombined according to the sub-bands, and the signals before and after combination are equal, which means the same signal x (t).
The result of each subband recombination is: a plurality of sub-signals with different frequencies obtained by Fourier decomposition in each frequency band are recombined into 1 frequency band signal. At this time, x (t) can be regarded as the sum of the signals of the frequency bands, that is:
wherein, Bn (t) is the recombined sub-band signal of the nth frequency band.
Comparing the terms on the right of the equations (2) and (3), it can be found that the term of the equation (2) is greatly reduced from the harmonic number 2M to the frequency band number 2N.
The signal frequency shift comprises the following steps: the complex signal X (t) is recombined in a segmented manner according to the frequency from small to large in the formula (3). For each sub-band signal bn (t), a frequency within a band (e.g., the center frequency of each band) may be selected to be shifted, respectively, that is:
wherein,
for the signals after frequency shift of the sub-bands, it is abbreviated as
ωrnFor a frequency-shifted angular frequency of the frequency band n,
is the upper and lower frequency limits, f, of the frequency band nhIs the frequency of the signal in the frequency band and satisfies
Balance
Is the sub-frequency band dynamic phasor (SFB-DP) of the signal. Obviously, compared with the conventional dynamic phasor Xh(t) dynamic phasors for sub-bands differing by only one frequency component
Is a signal having a bandwidth. If the bandwidth satisfies the narrow-band condition, then
Is also a low-frequency signal, and large-step simulation can be adopted for the low-frequency signal, thereby improving the simulation speed.
When the traditional dynamic phasor method is used for electromagnetic transient simulation, signals are firstly decomposed into dynamic phasors of each order according to the formula (2), and then the dynamic phasors of each order are placed in different CPU cores for parallel simulation by utilizing the characteristics of CPU multi-core cores. However, because the number of CPU cores is limited, the traditional dynamic phasor method can only adopt a few orders of dynamic phasors to form an equation set for parallel simulation. Therefore, the total bandwidth of the signals simulated by the traditional dynamic phasor method is far smaller than the actual bandwidth of the signals, so that the harmonic truncation error is large and the simulation precision is low. Different from the traditional dynamic phasor method, if the frequency-shifted sub-band dynamic phasor is decomposed according to the formula (3) for parallel simulation, because the sub-band signal has a certain bandwidth, under the condition that the frequency-shifted sub-band dynamic phasor has the same equation set number as the traditional dynamic phasor method, the bandwidth of the simulatable signal is far greater than that of the traditional dynamic phasor method, and the harmonic truncation error is far less than that of the traditional dynamic phasor method, so that the simulatable signal has extremely high simulation precision.
103, decomposing the simulation time into a plurality of time subintervals according to the operation condition; the method specifically comprises the following steps:
decomposing the simulation time T into a plurality of time subintervals according to the operation condition as follows: t ═ T1∪T2∪T3∪……∪TI；T1、T2、T3And TIRespectively, the 1 st, 2 nd, 3 rd and I th time subintervals are indicated.
104, based on the self-adaptive strategy principle, self-adaptively switching the window width of the multi-band dynamic phasor model according to the working condition of the voltage type current converter in each time subinterval, wherein the working condition is a steady-state working condition or a fault working condition;
and 105, simulating each time subinterval by using the window width switched by each time subinterval to realize the dynamic simulation of the voltage source type current converter.
Specifically, when simulation is performed, if the converter is in a steady-state working condition, the VSC multi-frequency is applied to the power frequency cycle time scaleAnd (3) performing simulation by using a large step length by using a power frequency period as a decomposition window width when the section dynamic phasor model, namely the VSC multi-band dynamic phasor model, is subjected to windowed Fourier decomposition. If the converter is disturbed (if a fault occurs), when the converter is in a transient working condition, the disturbance moment is advanced by a certain time quantum (the advanced time quantum is determined according to the actual simulation condition, and the advanced time quantum is 2T in the patentS，TSSwitching period), starting at a time after the time is advanced, applying a VSC multi-band dynamic phasor model on a switching period time scale, namely adopting the switching period as a decomposition window width when carrying out windowed Fourier decomposition on the VSC multi-band dynamic phasor model, and carrying out simulation by using small step length. If the simulation precision of the transient working condition of the converter is further improved, the rotating frequency of the multi-band dynamic vector can be zero, namely, a detailed model is adopted for simulation, and therefore extremely high simulation precision is obtained. FIG. 5 is a schematic diagram of an adaptive strategy for recovering a converter after the converter is disturbed and the operation condition changes.
The invention also provides a dynamic simulation system of the voltage source type converter, which specifically comprises:
and the dynamic vector model establishing module is used for establishing a dynamic vector model of the voltage source type current converter.
The dynamic vector model building module specifically comprises:
the dynamic vector model establishing submodule is used for establishing a dynamic vector model of the voltage source type current converter as follows:
wherein t represents time, ua、ub、ucRespectively representing three-phase voltages at an outlet side of the voltage source type converter; i.e. ia、ib、icRespectively representing three-phase currents on the inlet side of the voltage source type converter; u. ofdThe voltage of the direct current side of the voltage source type converter is represented, and R and L respectively represent the equivalent resistance and the equivalent inductance of the coupling transformer at the alternating current side of the voltage source type converter; < > represents a dynamic vector; sa、Sb、ScRespectively representing the switching functions of a three-phase full bridge arm of the voltage source type converter;
and the multi-band dynamic vector model acquisition module is used for decomposing, recombining and frequency shifting the dynamic vector model based on a multi-band dynamic vector method principle to obtain the multi-band dynamic vector model.
The multi-band dynamic vector model obtaining module specifically comprises:
the multi-band dynamic vector model obtaining submodule is used for decomposing, recombining and frequency shifting the dynamic vector model based on the principle of a multi-band dynamic vector method to obtain the multi-band dynamic vector model, and comprises the following steps:
wherein,
respectively representing the sub-band dynamic phasors of the three-phase current signals at the inlet side of the voltage source type converter,
respectively representing the sub-band dynamic phasors of the three-phase voltage signal at the outlet side of the voltage source type converter,
respectively representing sub-frequency band dynamic phasors of switching function signals of a three-phase full bridge arm of the voltage source type converter;
wherein x (t) represents ia、ib、ic、ua、ub、uc、Sa、SbOr Sc；fnFrequency of the nth frequency band signal, BnRepresenting the nth frequency band signal, Hn-1And HnRespectively representing the last harmonic order, X, of the (n-1) th and nth frequency band signalsh(t) denotes the h-Fourier coefficient of x (t), ωhDenotes the angular frequency, ω, of the h-th harmonic of x (t)sAt the fundamental angular frequency, f (T)i) The window width when performing windowed fourier decomposition on the state vector of the ith time subinterval is shown, and x (τ) represents the signal over a period τ of x (t).
And the simulation time decomposition module is used for decomposing the simulation time into a plurality of time subintervals according to the operating condition.
The simulation time decomposition module specifically comprises: the simulation time decomposition submodule is used for decomposing the simulation time T into a plurality of time subintervals according to the operation condition, and the time subintervals are as follows: t ═ T1∪T2∪T3∪……∪TI； T1、T2、T3And TIRespectively, the 1 st, 2 nd, 3 rd and I th time subintervals are indicated.
And the self-adaptive window width switching module is used for self-adaptively switching the window width of the multi-band dynamic phasor model according to the working condition of the voltage type current converter in each time subinterval based on the self-adaptive strategy principle, wherein the working condition is a steady-state working condition or a fault working condition.
The self-adaptive window width switching module comprises a first switching submodule and is used for switching the window width of the windowed Fourier to a power frequency period when the voltage source type converter is in a steady-state working condition in the ith time subinterval; wherein, I is 1,2,3 …, and I represents the number of sub-time intervals; and the second switching submodule is used for switching the window width of the windowed Fourier to a switching period when the voltage source type converter is in a transient working condition in the ith time subinterval.
And the simulation module is used for simulating each time subinterval by adopting the window width after each time subinterval is switched, so as to realize the dynamic simulation of the voltage source type current converter.
The method and the system of the invention have the following advantages: the self-adaptive strategy idea of the invention can self-adaptively switch the window width of each variable windowing Fourier decomposition according to the requirements of different operation working conditions, and uses large step length simulation under the steady state working condition and small step length simulation under the transient state working condition. The self-adaptive strategy specific method adopts a large window width to carry out Fourier decomposition under the steady-state working condition, the window width can be amplified to a power frequency period and is far larger than the window widths of other methods, and the simulation speed is improved. And under the transient working condition, the switching period is used as the window width, so that the simulation precision is ensured. The VSC self-adaptive multi-band dynamic phasor electromagnetic transient simulation model organically combines a self-adaptive strategy and a multi-band dynamic phasor method, so that the window width of the multi-band dynamic phasor method is variable and is not fixed any more when the multi-band dynamic phasor method carries out windowed Fourier decomposition on variables, and the simulation speed and the simulation precision are considered at the same time.
Specifically, the VSC adaptive multi-band dynamic phasor electromagnetic transient simulation model is compared with a VSC traditional dynamic phasor model.
1) The VSC self-adaptive multi-band dynamic phasor electromagnetic transient simulation model is different from the VSC traditional dynamic phasor model, and when the system frequency fluctuates (such as low-frequency oscillation), the problem of frequency spectrum leakage of a dynamic phasor method cannot occur.
2) Because the frequency range of the VSC self-adaptive multi-band dynamic phasor electromagnetic transient simulation model research is very large, the problem of harmonic wave does not need to be considered independently, and the harmonic wave truncation error is very small.
2. The VSC self-adaptive multi-band dynamic phasor electromagnetic transient simulation model is compared with a VSC frequency offset method model.
1) Compared with a VSC frequency offset method model, the VSC self-adaptive multi-band dynamic phasor electromagnetic transient simulation model has the advantages that the dominant frequency is the frequency band center frequency, the degree of freedom of selection is high, the defect that the VSC frequency offset method model can only shift the frequency of signals nearby a power frequency fundamental wave is overcome, and the simulation precision is high.
2) The number of sets of differential equations established by the VSC self-adaptive multi-band dynamic phasor electromagnetic transient simulation model is equal to the number N of sub-bands, compared with a VSC frequency offset method model, the number of equations is slightly increased, and due to mutual independence among the sub-bands, the CPU multi-core technology can be utilized for parallel solution, so that the simulation calculation speed can be greatly increased.
3. The VSC self-adaptive multi-band dynamic phasor electromagnetic transient simulation model is compared with the VSC multi-band dynamic phasor model.
The windowed Fourier decomposition window width of the VSC self-adaptive multi-band dynamic phasor electromagnetic transient simulation model can be changed, the defect that the window width of the VSC multi-band dynamic phasor model is not changed (usually power frequency) is overcome, large-step simulation can be adopted under the steady-state working condition, small-step simulation can be adopted under the transient working condition, and both the simulation precision and the simulation speed are considered.
The equivalent embodiments in the present specification are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts between the equivalent embodiments can be referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.
The principle and the implementation manner of the present invention are explained by applying specific examples, the above description of the embodiments is only used to help understanding the method of the present invention and the core idea thereof, the described embodiments are only a part of the embodiments of the present invention, not all embodiments, and all other embodiments obtained by a person of ordinary skill in the art based on the embodiments of the present invention without creative efforts belong to the protection scope of the present invention.
Claims (10)
1. A dynamic simulation method of a voltage source converter is characterized by comprising the following steps:
establishing a dynamic vector model of the voltage source type current converter;
decomposing, recombining and frequency shifting the dynamic vector model based on a multi-band dynamic vector method principle to obtain a multi-band dynamic vector model;
decomposing the simulation time into a plurality of time subintervals according to the operating condition;
based on the self-adaptive strategy principle, the window width of the multi-band dynamic phasor model is self-adaptively switched according to the working condition of the voltage type converter in each time subinterval, wherein the working condition is a steady-state working condition or a fault working condition;
and simulating each time subinterval by adopting the window width after each time subinterval is switched to realize the dynamic simulation of the voltage source type current converter.
2. The method according to claim 1, wherein the establishing a dynamic vector model of the voltage source converter specifically comprises:
the dynamic vector model of the voltage source type current converter is established as follows:
wherein t represents time, ua、ub、ucRespectively representing three-phase voltages at an outlet side of the voltage source type converter; i.e. ia、ib、icRespectively representing three-phase currents on the inlet side of the voltage source type converter; u. ofdThe voltage of the direct current side of the voltage source type converter is represented, and R and L respectively represent the equivalent resistance and the equivalent inductance of the coupling transformer at the alternating current side of the voltage source type converter; < > represents a dynamic vector; sa、Sb、ScRespectively representing the switching functions of a three-phase full bridge arm of the voltage source type converter;
3. the method according to claim 2, wherein the decomposing, recombining and frequency shifting the dynamic vector model based on the multi-band dynamic vector method principle to obtain the multi-band dynamic vector model specifically comprises:
based on the principle of a multi-band dynamic vector method, decomposing, recombining and frequency shifting the dynamic vector model to obtain the multi-band dynamic vector model is as follows:
wherein,
respectively representing the sub-band dynamic phasors of the three-phase current signals at the inlet side of the voltage source type converter,
respectively representing the sub-band dynamic phasors of the three-phase voltage signal at the outlet side of the voltage source type converter,
and
respectively representing sub-frequency band dynamic phasors of switching function signals of a three-phase full bridge arm of the voltage source type converter;
wherein x (t) represents ia、ib、ic、ua、ub、uc、Sa、SbOr Sc；fnFrequency of the nth frequency band signal, BnRepresenting the nth frequency band signal, Hn-1And HnRespectively representing the last harmonic order, X, of the (n-1) th and nth frequency band signalsh(t) denotes the h-Fourier coefficient of x (t), ωhDenotes the angular frequency, ω, of the h-th harmonic of x (t)sAt the fundamental angular frequency, f (T)i) The window width when performing windowed fourier decomposition on the state vector of the ith time subinterval is shown, and x (τ) represents the signal over a period τ of x (t).
4. The method according to claim 1, wherein the decomposing the simulation time into a plurality of time subintervals according to the operating condition specifically comprises:
decomposing the simulation time T into a plurality of time subintervals according to the operation condition as follows: t ═ T1∪T2∪T3∪……∪TI；T1、T2、T3And TIRespectively, the 1 st, 2 nd, 3 rd and I th time subintervals are indicated.
5. The dynamic simulation method of the voltage source converter according to claim 1, wherein the adaptively switching the window width of the multi-band dynamic phasor model according to the operating condition of the voltage source converter at each time subinterval based on the adaptive strategy principle, where the operating condition is a steady-state operating condition or a fault operating condition, specifically comprises:
when the voltage source type converter is in a steady-state working condition in the ith time subinterval, switching the window width of the windowed Fourier to a power frequency period; wherein, I is 1,2,3 …, and I represents the number of sub-time intervals;
and when the voltage source type converter is in a transient working condition in the ith time subinterval, switching the window width of the windowed Fourier to a switching period.
6. A dynamic simulation system of a voltage source converter is characterized by specifically comprising:
the dynamic vector model establishing module is used for establishing a dynamic vector model of the voltage source type current converter;
the multi-band dynamic vector model acquisition module is used for decomposing, recombining and frequency shifting the dynamic vector model based on a multi-band dynamic vector method principle to obtain a multi-band dynamic vector model;
the simulation time decomposition module is used for decomposing the simulation time into a plurality of time subintervals according to the operation working condition;
the adaptive window width switching module is used for adaptively switching the window width of the multi-band dynamic phasor model according to the working condition of the voltage type converter in each time subinterval based on an adaptive strategy principle, wherein the working condition is a steady-state working condition or a fault working condition;
and the simulation module is used for simulating each time subinterval by adopting the window width after each time subinterval is switched, so as to realize the dynamic simulation of the voltage source type current converter.
7. The dynamic simulation system of a voltage source converter according to claim 6, wherein the dynamic vector model building module specifically comprises:
the dynamic vector model establishing submodule is used for establishing a dynamic vector model of the voltage source type current converter as follows:
wherein t represents time, ua、ub、ucRespectively representing three-phase voltages at an outlet side of the voltage source type converter; i.e. ia、ib、icRespectively representing three-phase currents on the inlet side of the voltage source type converter; u. ofdThe voltage of the direct current side of the voltage source type converter is represented, and R and L respectively represent the equivalent resistance and the equivalent inductance of the coupling transformer at the alternating current side of the voltage source type converter; < > represents a dynamic vector; sa、Sb、ScRespectively representing the switching functions of a three-phase full bridge arm of the voltage source type converter;
8. the dynamic simulation system of a voltage source converter according to claim 7, wherein the multi-band dynamic vector model obtaining module specifically comprises:
the multi-band dynamic vector model obtaining submodule is used for decomposing, recombining and frequency shifting the dynamic vector model based on the principle of a multi-band dynamic vector method to obtain the multi-band dynamic vector model, and comprises the following steps:
wherein,
respectively representing the sub-band dynamic phasors of the three-phase current signals at the inlet side of the voltage source type converter,
respectively representing the sub-band dynamic phasors of the three-phase voltage signal at the outlet side of the voltage source type converter,
and
respectively representing sub-frequency band dynamic phasors of switching function signals of a three-phase full bridge arm of the voltage source type converter;
wherein x (t) represents ia、ib、ic、ua、ub、uc、Sa、SbOr Sc；fnFrequency of the nth frequency band signal, BnRepresenting the nth frequency band signal, Hn-1And HnRespectively representing the last harmonic order, X, of the (n-1) th and nth frequency band signalsh(t) denotes the h-Fourier coefficient of x (t), ωhDenotes the angular frequency, ω, of the h-th harmonic of x (t)sAt the fundamental angular frequency, f (T)i) The window width when performing windowed fourier decomposition on the state vector of the ith time subinterval is shown, and x (τ) represents the signal over a period τ of x (t).
9. The dynamic simulation system of the voltage source converter according to claim 6, wherein the simulation time decomposition module specifically comprises:
the simulation time decomposition submodule is used for decomposing the simulation time T into a plurality of time subintervals according to the operation condition, and the time subintervals are as follows: t ═ T1∪T2∪T3∪……∪TI；T1、T2、T3And TIRespectively, the 1 st, 2 nd, 3 rd and I th time subintervals are indicated.
10. The dynamic simulation system of a voltage source converter according to claim 6, wherein the adaptive window width switching module specifically comprises:
the first switching submodule is used for switching the window width of the windowed Fourier to a power frequency period when the voltage source type converter is in a steady-state working condition in the ith time subinterval; wherein, I is 1,2,3 …, and I represents the number of sub-time intervals;
and the second switching submodule is used for switching the window width of the windowed Fourier to a switching period when the voltage source type converter is in a transient working condition in the ith time subinterval.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010084659.5A CN111259565B (en) | 2020-02-10 | 2020-02-10 | Dynamic simulation method and system for voltage source type current converter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010084659.5A CN111259565B (en) | 2020-02-10 | 2020-02-10 | Dynamic simulation method and system for voltage source type current converter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111259565A true CN111259565A (en) | 2020-06-09 |
| CN111259565B CN111259565B (en) | 2021-12-14 |
Family
ID=70949206
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010084659.5A Active CN111259565B (en) | 2020-02-10 | 2020-02-10 | Dynamic simulation method and system for voltage source type current converter |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111259565B (en) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101908088A (en) * | 2010-07-22 | 2010-12-08 | 北京航空航天大学 | Time domain bidirectional iteration-based turbine vane flutter stress forecasting method |
| CN101968525A (en) * | 2010-10-09 | 2011-02-09 | 杭州市电力局 | Fault positioning method for power distribution network by combining simulation calculation and real-time monitoring |
| EP2319170A1 (en) * | 2008-08-26 | 2011-05-11 | ABB Technology AG | Control of a voltage source converter using synchronous machine emulation |
| US20110109280A1 (en) * | 2009-11-06 | 2011-05-12 | Abdel-Aty Edris | Automatic Self-Adjusting Voltage Controller |
| CN102542105A (en) * | 2011-12-24 | 2012-07-04 | 中国煤炭科工集团太原研究院 | Gear load wireless monitoring system and interactive multi-stage gear physical simulation method finished based on same |
| EP2549616A1 (en) * | 2011-07-22 | 2013-01-23 | ABB Research Ltd. | An arrangement and a method for supplying electric power |
| CN103326383A (en) * | 2013-05-31 | 2013-09-25 | 华北电力大学 | Wide area damping self-adaptation control system and control method thereof |
| CN105787149A (en) * | 2015-12-27 | 2016-07-20 | 北京航空航天大学 | Method for accurately converting power spectrum to tooth surface stress spectrum on drive system shaft of spiral bevel gear |
| CN108021719A (en) * | 2016-10-29 | 2018-05-11 | 南京理工大学 | A kind of wind farm grid-connected passive control method |
| CN110334935A (en) * | 2019-06-27 | 2019-10-15 | 南方电网科学研究院有限责任公司 | Method and device for evaluating transient stability of grid-connected exchanger and storage medium |
-
2020
- 2020-02-10 CN CN202010084659.5A patent/CN111259565B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2319170A1 (en) * | 2008-08-26 | 2011-05-11 | ABB Technology AG | Control of a voltage source converter using synchronous machine emulation |
| US20110109280A1 (en) * | 2009-11-06 | 2011-05-12 | Abdel-Aty Edris | Automatic Self-Adjusting Voltage Controller |
| CN101908088A (en) * | 2010-07-22 | 2010-12-08 | 北京航空航天大学 | Time domain bidirectional iteration-based turbine vane flutter stress forecasting method |
| CN101968525A (en) * | 2010-10-09 | 2011-02-09 | 杭州市电力局 | Fault positioning method for power distribution network by combining simulation calculation and real-time monitoring |
| EP2549616A1 (en) * | 2011-07-22 | 2013-01-23 | ABB Research Ltd. | An arrangement and a method for supplying electric power |
| CN102542105A (en) * | 2011-12-24 | 2012-07-04 | 中国煤炭科工集团太原研究院 | Gear load wireless monitoring system and interactive multi-stage gear physical simulation method finished based on same |
| CN103326383A (en) * | 2013-05-31 | 2013-09-25 | 华北电力大学 | Wide area damping self-adaptation control system and control method thereof |
| CN105787149A (en) * | 2015-12-27 | 2016-07-20 | 北京航空航天大学 | Method for accurately converting power spectrum to tooth surface stress spectrum on drive system shaft of spiral bevel gear |
| CN108021719A (en) * | 2016-10-29 | 2018-05-11 | 南京理工大学 | A kind of wind farm grid-connected passive control method |
| CN110334935A (en) * | 2019-06-27 | 2019-10-15 | 南方电网科学研究院有限责任公司 | Method and device for evaluating transient stability of grid-connected exchanger and storage medium |
Non-Patent Citations (2)
| Title |
|---|
| M. DARYABAK 等: "Modeling of LCC-HVDC Systems Using Dynamic Phasors", 《IEEE TRANSACTIONS ON POWER DELIVERY》 * |
| 孙栩 等: "VSC-HVDC系统的动态相量法建模仿真分析", 《电力系统自动化》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111259565B (en) | 2021-12-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kamwa et al. | Adaptive phasor and frequency-tracking schemes for wide-area protection and control | |
| CN110362937B (en) | Electromagnetic transient simulation method and system for modular multilevel converter | |
| CN110348161B (en) | Multi-frequency-band dynamic phasor electromagnetic transient simulation method and system for voltage source type converter | |
| CN102621382B (en) | Method for detecting frequency, phase and amplitude of electric signal in electric power system | |
| CN101572408B (en) | Method for improving simulation capacity of HVDC transmission system | |
| Jin et al. | A method for the identification of low frequency oscillation modes in power systems subjected to noise | |
| Dash et al. | Dynamic phasor and frequency estimation of time-varying power system signals | |
| CN110378022B (en) | Multi-frequency-band dynamic phasor electromagnetic transient simulation method and system | |
| CN110889210B (en) | Frequency shift transient simulation method, system, medium and device based on root matching method | |
| Usama et al. | Design and implementation of a wavelet analysis‐based shunt fault detection and identification module for transmission lines application | |
| Rivas et al. | Adaptive current harmonic estimation under fault conditions for smart grid systems | |
| CN104362628A (en) | Control method and device for avoiding large-wind-field resonance | |
| CN111310325B (en) | Dynamic simulation method and system of modular multilevel converter | |
| Samal et al. | New signal subspace approach to estimate the inter‐area oscillatory modes in power system using TLS‐ESPRIT algorithm | |
| Cassano et al. | Modeling of harmonic propagation of fast DC EV charging station in a Low Voltage network | |
| CN115133527A (en) | AC/DC coordination control method based on cascade hybrid DC multi-band additional damping control | |
| CN110378020B (en) | Multi-frequency-band dynamic phasor electromagnetic transient simulation method and system for power grid commutation converter | |
| CN111259565B (en) | Dynamic simulation method and system for voltage source type current converter | |
| CN112993988B (en) | Power grid line loss analysis method | |
| CN102854422B (en) | Three-phase unsymmetrical fault analyzing method for transformer branch | |
| CN110378021B (en) | Power transmission line simulation method and system | |
| CN103729502A (en) | Method for increasing electromagnetic transient simulation speed of power system | |
| CN103904653A (en) | Observability analyzing method in power grid harmonic state estimation | |
| CN117613881A (en) | New energy station harmonic responsibility division method based on principal component analysis method | |
| CN105140911A (en) | Power flow-based multi-direct-current feed-in receiving-end power grid voltage supporting capability evaluation method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Yao Yifan Inventor after: Wu Mengtong Inventor after: Ma Jiahao Inventor after: Yao Shujun Inventor after: Song Wenda Inventor after: Wang Yan Inventor after: Han Minxiao Inventor before: Yao Shujun Inventor before: Song Wenda Inventor before: Wang Yan Inventor before: Han Minxiao |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |