CN101950972A - SVC composite control method based on rapid equivalent susceptance calculation - Google Patents

SVC composite control method based on rapid equivalent susceptance calculation Download PDF

Info

Publication number
CN101950972A
CN101950972A CN 201010516681 CN201010516681A CN101950972A CN 101950972 A CN101950972 A CN 101950972A CN 201010516681 CN201010516681 CN 201010516681 CN 201010516681 A CN201010516681 A CN 201010516681A CN 101950972 A CN101950972 A CN 101950972A
Authority
CN
China
Prior art keywords
phase
current
equivalent susceptance
susceptance
sequence
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
Application number
CN 201010516681
Other languages
Chinese (zh)
Other versions
CN101950972B (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.)
Baiyun Power Group Co., Ltd.
Original Assignee
Hunan University
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 Hunan University filed Critical Hunan University
Priority to CN2010105166819A priority Critical patent/CN101950972B/en
Publication of CN101950972A publication Critical patent/CN101950972A/en
Application granted granted Critical
Publication of CN101950972B publication Critical patent/CN101950972B/en
Expired - Fee Related 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
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/30Reactive power compensation

Landscapes

  • Control Of Electrical Variables (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses an SVC composite control method based on rapid equivalent susceptance calculation, wherein, an SVC comprises a thyristor controlled reactor (TCR) and a passive filter bank. The SVC composite control method comprises the following steps: firstly, calculating the corresponding equivalent susceptance by using an equivalent susceptance computing method after detecting a voltage signal and a current signal; adding up an open-loop equivalent susceptance output value BLabc and a closed-loop PI regulator output value to obtain a total equivalent susceptance compensation value BTabc; calculating the actual input group number n of PPFs and the actual equivalent susceptance value of the TCR by a power divider according to the compensation value; and searching an Hash table by the TCR according to the actual equivalent susceptance value to obtain a corresponding triggering angle and then triggering opening of a corresponding thyristor to send out corresponding inductive reactive current for compensating capacitive reactive power and negative-sequence current, and finally realizing complete compensation of the reactive power and the negative-sequence current of a system. The SVC composite control method effectively reduces the calculated quantity of the system and improves the response time of the system.

Description

A kind of SVC composite control method that calculates based on quick equivalent susceptance
Technical field
The invention belongs to negative phase-sequence and reactive power compensation control technology field, particularly a kind of SVC composite control method that calculates based on quick equivalent susceptance.
Background technology
Along with The development in society and economy, the power consumption of industrial circles such as metallurgy, mine and traffic is increasing.A large amount of utilizations as high-power perception such as gold metallurgy rectifying installation, arc furnace and electric locomotive and non-linear equipment, the negative phase-sequence of bringing thus, power quality problems such as idle and harmonic wave are also more and more serious, a large amount of negative phase-sequence and idlely not only caused a large amount of electric energy losses has also reduced the safety and stability performance of system simultaneously.
At present, SVC is acknowledged as the best equipment of reactive power compensation with its cost performance height, response characteristics such as fast, carry out the equilibrating compensation to eliminate negative sequence component but also can regulate by phase-splitting to the three-phase asymmetric load.Most SVC are used for realizing that reactive power compensation improves the power factor or the stable line voltage of system, and it mostly is the PI control method, and its control performance and precision also are difficult to reach requirement, and can not take into account imbalance compensation.Existing domestic and international list of references utilizes the phase splitting compensation principle of Steinmetz, and the method for regulating by phase-splitting is used for compensating negative-sequence current and reactive current with SVC, realizes the balance of system's three phase power.The voltage and current signal positive sequence component that the scholar generally adopts synchronous coordinate conversion or instantaneous reactive power theory to obtain, phase splitting compensation principle by Steinmetz can calculate three-phase compensation susceptance value then, because this control method is open loop control, compensation precision is not high, be subjected to the influence of the system's detection and the error of calculation easily, can not compensate by the real-time tracking power network current.
For compensation precision and the performance of dynamic tracking that improves SVC, Taiwan's scholars S.-Y.Lee professor proposes a kind of SVC Mathematical Modeling based on negative-sequence current compensation, and the The optimal compensation amount that calculates SVC by dynamic iteration and search realizes the dynamic compensation to unbalanced load.But this control method needs complicated iteration and search to calculate.
At present, Chang Yong detection algorithm is based on the p-q algorithm and the ip-iq algorithm of instantaneous reactive power theory.Wherein, the accuracy in detection of p-q method is influenced by grid voltage waveform, and testing result has error when there is certain distortion in line voltage; And ip-iq algorithm detection method has been owing to comprised space coordinate transformation 2 times, and algorithm is complicated, is unfavorable for that software realizes.
The FBD method of the new a kind of time domain detection that proposes of German scholar in recent years.Its explicit physical meaning, real-time be good, be applicable to single-phase or polyphase system.The basic ideas of FBD method are: is the load equivalent in the side circuit desired electrical guiding element, thinks that the power in the circuit all consumes on this equivalent conductance.According to equivalent conductance electric current is decomposed, the character of each current component is discussed.It is compared with the real-time detection method based on the three-phase instantaneous reactive power theory, improved the rapidity of detection method, Park conversion and 3/2 conversion have been saved, need not complicated matrixing, algorithm obtains simplifying, greatly reduced the expense of microprocessor, time-delay is further shortened, and also lays the first stone for further improving sample frequency and compensation precision.
Summary of the invention
The defective that exists at the open-loop control method of SVC negative phase-sequence and reactive power compensation, the present invention aims to provide a kind of SVC composite control method that calculates based on quick equivalent susceptance, can calculate the equivalent susceptance of SVC fast by the FBD algorithm, its explicit physical meaning, directly fast, effectively improved the compensation performance of SVC.Utilize the open loop feedforward action, realize idle and quick compensation negative-sequence current, improve system response time; Utilize the indifference tracking characteristics of close-loop feedback control, remedy the deficiency of open loop control, realize the stable state indifference, like this with open loop control and the effective combination of closed-loop control, realize negative phase-sequence in the centering high-voltage fence electric current and idle full remuneration, improved the quality of power supply of electric power system.
To achieve the above object of the invention, the technical solution adopted in the present invention is:
A kind of SVC composite control method that calculates based on quick equivalent susceptance comprises the steps:
A), at first by detecting three phase network voltage (u a, u b, u c), utilize the positive sequence voltage detection method can calculate the three-phase positive sequence voltage (
Figure BDA0000029189870000031
With );
B), detect the threephase load current i Labc, utilize the equivalent susceptance value B that tries to achieve threephase load electric current correspondence based on the quick equivalent susceptance computational methods of FBD algorithm LabcBy detecting the threephase load current i Labc, three-phase TCR current i TCRabcWith passive filter group current i PPFabc, their additions can be calculated the three phase network current i Sabc, again by the three phase network current i that obtains SabcCalculate the equivalent susceptance value B of three phase network electric current correspondence according to quick equivalent susceptance computational methods based on the FBD algorithm Sabc
C), the equivalent susceptance value B of the three phase network electric current correspondence that obtains of above-mentioned steps Sabc, after the adjusting of PI controller, can calculate the closed loop output valve
Figure BDA0000029189870000033
Again with the closed loop output valve that obtains
Figure BDA0000029189870000034
The corresponding equivalent susceptance value B of threephase load electric current that obtains with above-mentioned steps LabcAddition, the three-phase that obtains the total output of compound control compensates susceptance B Tabc(B Tab, B Tbc, B Tca);
D), the three-phase according to the total output of compound control compensates susceptance B Tabc, the actual needs that draws SVC by power divider drops into the actual equivalent susceptance value that PPF organizes number n and TCR
Figure BDA0000029189870000035
Controller sends control signal corresponding and cuts n group PPF like this, for system provides corresponding reactive capability; Simultaneously, TCR is according to actual equivalent susceptance value Show to obtain relative trigger angle α by searching Hash, pass through the driving and the amplification of drive circuit then, trigger opening of corresponding thyristor and send corresponding perceptual reactive current and come compensating reactive power and negative-sequence current, the system that finally realizes is idle and the full remuneration of negative-sequence current.
Described quick equivalent susceptance computational methods based on the FBD algorithm may further comprise the steps:
A), establishing three phase network voltage and three-phase current is:
u a = U sin ( wt ) u b = U sin ( wt - 2 π / 3 ) u c = U sin ( wt + 2 π / 3 ) - - - ( 1 )
Figure BDA0000029189870000041
Wherein, U is the effective value of three-phase voltage;
Figure BDA0000029189870000042
With
Figure BDA0000029189870000043
Effective value and phase angle for n forward-order current;
Figure BDA0000029189870000044
With
Figure BDA0000029189870000045
Effective value and phase angle for n negative-sequence current; Make u 1=(u a, u b, u c), i=(i a, i b, i c), with voltage u 1Each component obtain new vector behind the phase shift pi/2 forward and be designated as
Figure BDA0000029189870000046
According to the power definition of FBD principle, have that fundamental positive sequence equivalence is meritorious to be respectively ∑ p with reactive power +=<u 1, i 〉,
Figure BDA0000029189870000047
So can push away below fundamental positive sequence equivalence meritorious with idle electricity lead into:
G p + ( t ) = Σ p + | | u 1 | | 2 = 2 [ u a i a + u b i b + u c i c ] 3 U 2 G q + ( t ) = Σ q + | | u i 1 | | 2 = 2 [ u i a i a + u i b i b + u i c i c ] 3 U 2 - - - ( 3 )
If in like manner with u 1In u bAnd u cExchange one upper/lower positions then can obtain other one group of voltage vector u mutually 2(u 2=(u a, u c, u b)) and
Figure BDA0000029189870000049
Then can push away to such an extent that the meritorious electricity of first-harmonic negative phase-sequence equivalence is led the expression formula of leading with idle electricity:
G p - ( t ) = Σ p - | | u 2 | | 2 = 2 ( u a i a + u c i b + u b i c ) 3 U 2 G q - ( t ) = Σ q - | | u i 2 | | 2 = 2 ( u i a i a + u i c i b + u i b i c ) 3 U 2 - - - ( 4 )
Wherein have
Figure BDA00000291898700000411
U is the effective value of phase voltage;
B), with in the formula (3)
Figure BDA00000291898700000412
Behind low pass filter LPF filtering Alternating Component, obtain the direct current composition that the meritorious electricity of fundamental positive sequence equivalence is led
Figure BDA00000291898700000413
The direct current composition of leading with idle electricity
Figure BDA00000291898700000414
With in the formula (4)
Figure BDA00000291898700000415
Behind the low pass filter filters out Alternating Component, obtain the direct current composition that the meritorious electricity of first-harmonic negative phase-sequence equivalence is led
Figure BDA00000291898700000416
The direct current composition of leading with idle electricity Then according to the definition of positive sequence and negative-sequence current, promptly can push away following formula:
Wherein, U is the effective value of phase voltage; I +And I -Be fundamental positive sequence and negative-sequence current vector;
Figure BDA0000029189870000052
With
Figure BDA0000029189870000053
Be fundamental positive sequence and negative-sequence current effective value;
Figure BDA0000029189870000054
With
Figure BDA0000029189870000055
Power-factor angle for fundamental positive sequence and negative-sequence current;
C), realize that equilibrating compensation and raising power factor of electric network approach 1, just must satisfy following 2 conditions: the negative sequence component of 1) eliminating load current; 2) making the positive sequence component imaginary part of load current is 0; According to above-mentioned condition, can obtain the equivalent susceptance value that TCR need drop into electrical network and be:
B ab = - ( G q + - - 3 G p - - + G q - - ) / 3 B bc = - ( G q + - - 2 G q - - ) / 3 B ca = - ( G q + - + 3 G p - - + G q - - ) / 3 - - - ( 6 )
Utilize in the step a) positive sequence voltage detection method calculate the three-phase positive sequence voltage (
Figure BDA0000029189870000057
With
Figure BDA0000029189870000058
) step as follows:
Three phase network voltage (u a, u b, u c) through three or two coordinate transforms, be converted into dq reference axis direct current signal u dAnd u q, the three phase network voltage transitions is that dq axial coordinate computing formula is as follows:
u d u q = 2 3 sin wt sin ( wt - 2 π / 3 ) sin ( wt + 2 π / 3 ) - cos wt - cos ( wt - 2 π / 3 ) - cos ( wt + 2 π / 3 ) u a u b u c - - - ( 7 )
Wherein w is meant the first-harmonic angular frequency, the direct current signal u under the above-mentioned dq reference axis that obtains dAnd u q, after low-pass filtering, inverse transformation transforms, and is reduced to the three-phase positive sequence voltage
Figure BDA00000291898700000510
With Computing formula is as follows:
u a + u b + u c + = 2 3 sin wt - cos wt sin ( wt - 2 π / 3 ) - cos ( wt - 2 π / 3 ) sin ( wt + 2 π / 3 ) - cos ( wt + 2 π / 3 ) u d + u q + - - - ( 8 )
Further, in the step d), by total output B of compound control AbcThe actual needs that draws SVC drops into the actual equivalent susceptance value that PPF organizes number n and TCR
Figure BDA00000291898700000513
Power distribution method be:
If compensation capacity of each group passive filter is Q, then the passive filter group of input compensation is counted n and should be satisfied following condition:
max { B Tab , B Tbc , B Tca } < min | n Q 3 U 2 | - - - ( 9 )
Wherein U is the effective value of phase voltage, B Tab, B Tbc, B TcaThe total three-phase of exporting of compound control that calculates for step c) compensates susceptance.After dropping into n group passive filter, because three-phase power supply system is in capacitive and crosses the benefit state, idle to reach idle and full remuneration negative phase-sequence by the corresponding perception of TCR compensation, then the perceptual idle susceptance of TCR actual needs input is:
B ab * = B Tab - nQ 3 U 2 B bc * = B Tbc - nQ 3 U 2 B ca * = B Tca - nQ 3 U 2 - - - ( 10 )
Further, in the step d), TCR is according to actual equivalent susceptance value By searching the specific implementation step that Hash shows to obtain relative trigger angle α be:
Because between the equivalent susceptance that trigger angle α and TCR need compensate is a kind of non-linear relation, following is not:
B = 2 &pi; - 2 &alpha; - sin ( 2 &pi; - 2 &alpha; ) &pi;X - - - ( 11 )
Wherein X is the reactance value of connecting with thyristor, and α is the trigger angle of TCR; Therefore at first the following formula linearity is changed into the Hash question blank, because the excursion of the trigger angle α of TCR is [90 °, 180 °], trigger angle α is since 90 °, with 0.5 ° as incremental change, each trigger angle α can calculate an equivalent susceptance value, (α B), can calculate 180 groups of vectors altogether so just to constitute a vector, deposit in the corresponding Hash question blank, so just finished the linearization process of equivalent susceptance B to trigger angle α; And then according to actual equivalent susceptance value
Figure BDA0000029189870000065
Obtain relative trigger angle α by searching described Hash question blank.
Based on above-mentioned steps, technique effect of the present invention is:
1. owing to need not carry out complicated Park conversion, DQ conversion and inverse transformation thereof in the FBD algorithm that the present invention adopts, and end product is not subjected to the influence of zero-sequence current, reduced the load of controller, described TCR equivalent susceptance acquiring method based on the FBD algorithm can directly detect SVC instruction references signal fast, thereby improved the quality of the idle and negative sequence compensation of system.
2. by in conjunction with open loop and closed loop control method, remedy the shortcoming that open loop control compensation precision is not high and examined and error of calculation influence is big, improved the compensation precision of SVC.
In sum, the SVC composite control method that calculates based on quick equivalent susceptance of the present invention can improve the compensation precision and the dynamic property of SVC device.After adopting composite control method proposed by the invention, the negative phase-sequence and the reactive power compensation of electric power system are all improved greatly.
Description of drawings
Fig. 1 is imbalance compensation entire system block diagram among the present invention;
Fig. 2 among the present invention based on the quick equivalent susceptance computing block diagram of FBD algorithm;
Fig. 3 detects block diagram for three-phase positive sequence voltage among the present invention;
Fig. 4 is SVC composite control method block diagram described in the present invention;
In above-mentioned accompanying drawing,
1-is a thyristor-controlled reactor; The 2-passive filter; The 3-upper computer detection system;
The 4-DSP control system; The 5-Fiber Optical Communication System; 6-drives triggering system;
7-static reactive system.
Embodiment
As shown in Figure 1, include in the imbalance compensation system among the present invention by being connected in the Static Var Compensator SVC7 that thyristor-controlled reactor 1 in the electrical network and passive filter 2 are formed.Wherein thyristor-controlled reactor 1 is taked delta connection, each branch road is taked two equal-sized reactor series connection, the antiparallel two-way thyristor mode of connection placed in the middle, can make the voltage to earth of thyristor valve when normal operation lower like this, and when a reactor is short-circuited fault, another reactor still can play the effect of limiting short-circuit current, has strengthened reliability of system operation.Wherein PPF passive filter 2 not only can provide jumbo capacitive reactive power, but also has certain harmonic filtration function, can reduce the electrical network aberration rate.
The negative phase-sequence that causes at load, the idle power quality problem that waits, present embodiment takes the imbalance compensation control method to realize the balance of each phase power, the system configuration of this control method of specific implementation comprises the upper computer detection system 3, DSP control system 4, Fiber Optical Communication System 5 and the trigger board 6 that connect successively as shown in Figure 1.Wherein upper computer detection system 3 is used to gather voltage and current signal, and finishes the calculating of three-phase compensation equivalent susceptance, and the result is passed to DSP control system 4.DSP control system 4 is mainly finished the ripple of sending out of signal, and triggering signal is transferred to trigger board 6 through Fiber Optical Communication System 5, and signal triggers thyristor by the driving and the amplification of trigger board 6.This photoelectricity triggering mode has realized that not only the isolation of high electronegative potential also has comparatively ideal anti-electromagnetic interference capability, and response speed is fast, helps system and moves reliably and with long-term.
As shown in Figure 2, among the present invention based on the equivalent susceptance computing block diagram of FBD algorithm.
If three-phase voltage and electric current are:
u a = U sin ( wt ) u b = U sin ( wt - 2 &pi; / 3 ) u c = U sin ( wt + 2 &pi; / 3 ) - - - ( 12 )
Wherein, U is the effective value of three-phase voltage;
Figure BDA0000029189870000083
With
Figure BDA0000029189870000084
Effective value and phase angle for n forward-order current;
Figure BDA0000029189870000085
With
Figure BDA0000029189870000086
Effective value and phase angle for n negative-sequence current.Make u 1=(u a, u b, u c), i=(i a, i b, i c), with voltage u 1Each component obtain new vector behind the phase shift pi/2 forward and be designated as
Figure BDA0000029189870000087
According to the power definition of FBD principle, have that fundamental positive sequence equivalence is meritorious to be respectively ∑ p with reactive power +=<u 1, i 〉,
Figure BDA0000029189870000088
So can push away below fundamental positive sequence equivalence meritorious with idle electricity lead into:
G p + ( t ) = &Sigma; p + | | u 1 | | 2 = 2 [ u a i a + u b i b + u c i c ] 3 U 2 G q + ( t ) = &Sigma; q + | | u i 1 | | 2 = 2 [ u i a i a + u i b i b + u i c i c ] 3 U 2 - - - ( 14 )
To obtain in the following formula
Figure BDA00000291898700000810
Behind low pass filter LPF filtering Alternating Component, obtain the direct current composition that the meritorious electricity of fundamental positive sequence equivalence is led
Figure BDA00000291898700000811
The direct current composition of leading with idle electricity
Figure BDA00000291898700000812
If in like manner with u 1In u bAnd u cExchange one upper/lower positions then can obtain other one group of voltage vector u mutually 2(u 2=(u a, u c, u b)) and
Figure BDA0000029189870000091
Then can push away to such an extent that the meritorious electricity of first-harmonic negative phase-sequence equivalence is led the expression formula of leading with idle electricity:
G p - ( t ) = &Sigma; p - | | u 2 | | 2 = 2 ( u a i a + u c i b + u b i c ) 3 U 2 G q - ( t ) = &Sigma; q - | | u i 2 | | 2 = 2 ( u i a i a + u i c i b + u i b i c ) 3 U 2 - - - ( 15 )
Wherein have
Figure BDA0000029189870000093
U is the effective value of phase voltage.To obtain in the following formula
Figure BDA0000029189870000094
Behind the low pass filter filters out Alternating Component, obtain the direct current composition that the meritorious electricity of first-harmonic negative phase-sequence equivalence is led The direct current composition of leading with idle electricity
Figure BDA0000029189870000096
Then according to the definition of positive sequence and negative-sequence current, promptly can push away following formula:
Figure BDA0000029189870000097
Wherein, U is the effective value of phase voltage; I +And I -Be fundamental positive sequence and negative-sequence current vector;
Figure BDA0000029189870000098
With
Figure BDA0000029189870000099
Be fundamental positive sequence and negative-sequence current effective value; With
Figure BDA00000291898700000911
Power-factor angle for fundamental positive sequence and negative-sequence current.Realize that equilibrating compensation and raising power factor of electric network approach 1, just must satisfy following 2 conditions: the negative sequence component of 1) eliminating load current; 2) making the positive sequence component imaginary part of load current is 0.According to above-mentioned condition, can obtain the equivalent susceptance value that TCR need drop into electrical network and be:
B ab = - ( G q + - - 3 G p - - + G q - - ) / 3 B bc = - ( G q + - - 2 G q - - ) / 3 B ca = - ( G q + - + 3 G p - - + G q - - ) / 3 - - - ( 17 )
By above-mentioned analysis as can be known, zero-sequence current does not influence testing result, so that the FBD algorithm is applicable to is single-phase, the current detecting of three-phase three-wire system and three-phase four-wire system.And the FBD algorithm do not carry out complicated Park conversion, DQ conversion and inverse transformation thereof, reduced the calculated load amount of controller.
Since when equivalent susceptance less than 0 the time, then to drop into the idle susceptance of corresponding perception, otherwise, then need to drop into corresponding capacitive reactive power susceptance.According to the size of load reactive power demand, design N group passive filter carries out the compensation of capacitive reactive power, and the compensation capacity of each group passive filter is Q.Then the passive filter group of input compensation is counted n and should be satisfied following condition:
max { B Tab , B Tbc , B Tca } < min | n Q 3 U 2 | - - - ( 18 )
At first drop into after the n group passive filter, three-phase power supply system is in capacitive and crosses the benefit state, and is idle to reach idle and full remuneration negative phase-sequence by the corresponding perception of TCR compensation then.Then the perceptual idle susceptance of TCR actual needs input is:
B ab * = B Tab - nQ 3 U 2 B bc * = B Tbc - nQ 3 U 2 B ca * = B Tca - nQ 3 U 2 - - - ( 19 )
Wherein U is the effective value of phase voltage, B Tab, B Tbc, B TcaThree-phase compensation susceptance for the total output of compound control.By electrical network three-phase voltage and current signal are detected, utilize the FBD algorithm computation to obtain accordingly
Figure BDA0000029189870000103
Figure BDA0000029189870000104
And, utilize the composite control method of open loop and closed loop can calculate the total equivalent susceptance B of compound control output by corresponding processing and calculating Tabc, calculate the actual equivalent susceptance value that the actual input group of PPF is counted n and TCR by above-mentioned power division principle Controller sends control signal corresponding and cuts n group PPF like this, for system provides corresponding reactive capability, obtains the full remuneration that corresponding perceptual reactive current realizes negative-sequence current by triggering TCR again.Because be a kind of non-linear relation between the equivalent susceptance that trigger angle α and TCR need compensate, as follows:
B = 2 &pi; - 2 &alpha; - sin ( 2 &pi; - 2 &alpha; ) &pi;X - - - ( 20 )
Wherein X is the reactance value of connecting with thyristor, and α is the trigger angle of TCR.The control of TCR generally changes into the Hash question blank with the following formula linearity for convenience, like this can be according to the equivalent susceptance value, and table look-at obtains corresponding TCR trigger angle.Because the excursion of the trigger angle of TCR is [90 °, 180 °], α is since 90 °, with 0.5 ° as incremental change, each α can calculate an equivalent susceptance value, so just constitutes a vector (α, B), can calculate 180 groups of vectors altogether, deposit in the corresponding Hash question blank, so just finish the linearization process of equivalent susceptance B to trigger angle α.
As shown in Figure 3, positive sequence voltage detects block diagram among the present invention.
Because the imbalance and the harmonic current phenomenon of enterprise's electric power system load are serious, make to exist a large amount of negative phase-sequences and harmonic wave composition in the supply power voltage, make result of calculation that very big error will be arranged.Therefore present embodiment adopts instantaneous reactive power theory to extract the positive sequence voltage detection method, shown in figure (3).Three-phase voltage signal (u a, u b, u c) through three or two coordinate transforms, three-phase voltage signal is converted into dq reference axis direct current signal u dAnd u q, come filtering to convert corresponding AC signal to by low pass filter then by other negative phase-sequence and harmonic voltage, can obtain the direct current signal that forms by three-phase fundamental positive sequence voltage transitions like this With
Figure BDA0000029189870000112
Like this
Figure BDA0000029189870000113
With
Figure BDA0000029189870000114
To change through inverse transformation, be reduced to three-phase fundamental positive sequence voltage
Figure BDA0000029189870000115
With
Figure BDA0000029189870000116
It is as follows that three-phase voltage is converted to dq axial coordinate computing formula:
u d u q = T abc / dq ( w ) u a u b u c
u d u q = 2 3 sin wt sin ( wt - 2 &pi; / 3 ) sin ( wt + 2 &pi; / 3 ) - cos wt - cos ( wt - 2 &pi; / 3 ) - cos ( wt + 2 &pi; / 3 ) u a u b u c - - - ( 21 )
U under the above-mentioned dq reference axis that obtains dAnd u q, after low-pass filtering, inverse transformation is converted into the component of voltage under three phase coordinates, and computing formula is as follows:
u a + u b + u c + = T dq / abc ( w ) u d + u q +
u a + u b + u c + = 2 3 sin wt - cos wt sin ( wt - 2 &pi; / 3 ) - cos ( wt - 2 &pi; / 3 ) sin ( wt + 2 &pi; / 3 ) - cos ( wt + 2 &pi; / 3 ) u d + u q + - - - ( 22 )
Like this by above-mentioned positive sequence voltage detection method, effective negative phase-sequence and the harmonic component in the filtering three-phase voltage component improves the counting accuracy and the compensation performance of SVC bucking-out system.
As shown in Figure 4, SVC composite control method of the present invention, as described below:
Because the nonlinear characteristic of SVC, its control performance will inevitably be affected.In order to improve the dynamic property of SVC, the present invention combines traditional closed-loop and open loop control, utilizes the fast-response energy of feedforward open loop control, utilizes the performance of dynamic tracking of closed-loop control to carry out online self-adjusting again, makes the dynamic property of SVC more superior.After compensation is stable, power network current three-phase power balance, negative phase-sequence and reactive current are 0, thereby should be 0 according to the equivalent susceptance that power network current calculates, so can realize the dynamic tracking compensation of negative phase-sequence and reactive current with the equivalent susceptance of power network current correspondence as the closed-loop control target.
Concrete implementation step is as follows:
1), at first by detecting three phase network voltage (u a, u b, u c), utilize the positive sequence voltage detection method can calculate the three-phase positive sequence voltage (
Figure BDA0000029189870000121
With
Figure BDA0000029189870000122
);
2), detect the threephase load current i Labc, utilize the equivalent susceptance value B that tries to achieve threephase load electric current correspondence based on the quick equivalent susceptance computational methods of FBD algorithm LabcBy detecting the threephase load current i Labc, three-phase TCR current i TCRabcWith passive filter group current i PPFabc, their additions can be calculated the three phase network current i Sabc, the quick equivalent susceptance computational methods of basis can calculate the equivalent susceptance value B of three phase network electric current correspondence thus Sabc
3), the equivalent susceptance B of the three phase network electric current correspondence that obtains such as previous step Sabc, after the adjusting of PI controller, can calculate the closed loop output valve
Figure BDA0000029189870000123
(
Figure BDA0000029189870000124
k p, k iBe expressed as ratio and integral coefficient respectively), the corresponding equivalent susceptance value B of threephase load electric current that calculates suddenly with previous step like this LabcAddition can obtain total output B of compound control Tabc
Total output B according to compound control Tabc, the actual needs that calculates SVC by power divider drops into the actual equivalent susceptance value that PPF organizes number n and TCR
Figure BDA0000029189870000125
Controller sends control signal corresponding and cuts n group PPF like this, for system provides corresponding reactive capability.Yet TCR is according to actual equivalent susceptance value Show to obtain the relative trigger angle by searching Hash, pass through the driving and the amplification of drive circuit then, trigger opening of corresponding thyristor and send corresponding perceptual reactive current and come compensating reactive power and negative-sequence current, the system that finally realizes is idle and the full remuneration of negative-sequence current.Utilize the dynamic adjustments ability of closed-loop control to realize steady-state error like this, effectively remedied the shortcoming of open loop control, improve dynamic compensation performance and the compensation precision of SVC, realize negative phase-sequence and idle compensation and the balance of three phase power.

Claims (5)

1. a SVC composite control method that calculates based on quick equivalent susceptance is characterized in that, comprises the steps:
A), at first by detecting three phase network voltage (u a, u b, u c), utilize the positive sequence voltage detection method can calculate the three-phase positive sequence voltage (
Figure FDA0000029189850000011
With
Figure FDA0000029189850000012
);
B), detect the threephase load current i Labc, utilize the equivalent susceptance value B that tries to achieve threephase load electric current correspondence based on the quick equivalent susceptance computational methods of FBD algorithm LabcBy detecting the threephase load current i Labc, three-phase TCR current i TCRabcWith passive filter group current i PPFabc, their additions can be calculated the three phase network current i Sabc, thus according to the equivalent susceptance value B that can calculate three phase network electric current correspondence based on the quick equivalent susceptance computational methods of FBD algorithm Sabc
C), the equivalent susceptance B of the three phase network electric current correspondence that obtains such as previous step Sabc, after the adjusting of PI controller, can calculate the closed loop output valve
Figure FDA0000029189850000013
Again with the closed loop output valve that obtains The corresponding equivalent susceptance value B of threephase load electric current that obtains with above-mentioned steps LabcAddition, the three-phase that obtains the total output of compound control compensates susceptance B Tabc(B Tab, B Tbc, B Tca);
D), the three-phase according to the total output of compound control compensates susceptance B Tabc, the actual needs that calculates SVC by power divider drops into the actual equivalent susceptance value that PPF organizes number n and TCR
Figure FDA0000029189850000015
Controller sends control signal corresponding and cuts n group PPF like this, for system provides corresponding reactive capability; Simultaneously, TCR is according to actual equivalent susceptance value
Figure FDA0000029189850000016
Show to obtain the relative trigger angle by searching Hash, pass through the driving and the amplification of drive circuit then, trigger opening of corresponding thyristor and send corresponding perceptual reactive current and come compensating reactive power and negative-sequence current, the system that finally realizes is idle and the full remuneration of negative-sequence current.
2. according to the described SVC composite control method that calculates based on quick equivalent susceptance of claim 1, it is characterized in that described quick equivalent susceptance computational methods based on the FBD algorithm may further comprise the steps:
A), establishing three phase network voltage and three-phase current is:
u a = U sin ( wt ) u b = U sin ( wt - 2 &pi; / 3 ) u c = U sin ( wt + 2 &pi; / 3 ) - - - ( 1 )
Figure FDA0000029189850000022
Wherein, U is the effective value of three-phase voltage;
Figure FDA0000029189850000023
With Effective value and phase angle for n forward-order current;
Figure FDA0000029189850000025
With
Figure FDA0000029189850000026
Effective value and phase angle for n negative-sequence current; Make u 1=(u a, u b, u c), i=(i a, i b, i c), with voltage u 1Each component obtain new vector behind the phase shift pi/2 forward and be designated as
Figure FDA0000029189850000027
According to the power definition of FBD principle, have that fundamental positive sequence equivalence is meritorious to be respectively ∑ p with reactive power +=<u 1, i 〉,
Figure FDA0000029189850000028
So can push away below fundamental positive sequence equivalence meritorious with idle electricity lead into:
G p + ( t ) = &Sigma; p + | | u 1 | | 2 = 2 [ u a i a + u b i b + u c i c ] 3 U 2 G q + ( t ) = &Sigma; q + | | u i 1 | | 2 = 2 [ u i a i a + u i b i b + u i c i c ] 3 U 2 - - - ( 3 )
If in like manner with u 1In u bAnd u cExchange one upper/lower positions then can obtain other one group of voltage vector u mutually 2(u 2=(u a, u c, u b)) and
Figure FDA00000291898500000210
Then can push away to such an extent that the meritorious electricity of first-harmonic negative phase-sequence equivalence is led the expression formula of leading with idle electricity:
G p - ( t ) = &Sigma; p - | | u 2 | | 2 = 2 ( u a i a + u c i b + u b i c ) 3 U 2 G q - ( t ) = &Sigma; q - | | u i 2 | | 2 = 2 ( u i a i a + u i c i b + u i b i c ) 3 U 2 - - - ( 4 )
Wherein have
Figure FDA00000291898500000212
U is the effective value of phase voltage;
B), with in the formula (3)
Figure FDA00000291898500000213
Behind low pass filter LPF filtering Alternating Component, obtain the direct current composition that the meritorious electricity of fundamental positive sequence equivalence is led
Figure FDA00000291898500000214
The direct current composition of leading with idle electricity
Figure FDA00000291898500000215
With in the formula (4)
Figure FDA00000291898500000216
Behind the low pass filter filters out Alternating Component, obtain the direct current composition that the meritorious electricity of first-harmonic negative phase-sequence equivalence is led
Figure FDA00000291898500000217
The direct current composition of leading with idle electricity
Figure FDA00000291898500000218
Then according to the definition of positive sequence and negative-sequence current, promptly can push away following formula:
Wherein, U is the effective value of phase voltage; I +And I -Be fundamental positive sequence and negative-sequence current vector; With Be fundamental positive sequence and negative-sequence current effective value;
Figure FDA0000029189850000034
With Power-factor angle for fundamental positive sequence and negative-sequence current;
C), realize that equilibrating compensation and raising power factor of electric network approach 1, just must satisfy following 2 conditions: the negative sequence component of 1) eliminating load current; 2) making the positive sequence component imaginary part of load current is 0; According to above-mentioned condition, can obtain the equivalent susceptance value that TCR need drop into electrical network and be:
B ab = - ( G q + - - 3 G p - - + G q - - ) / 3 B bc = - ( G q + - - 2 G q - - ) / 3 B ca = - ( G q + - + 3 G p - - + G q - - ) / 3 . - - - ( 6 )
3. according to the described SVC composite control method that calculates based on quick equivalent susceptance of claim 1, it is characterized in that, utilize in the step a) positive sequence voltage detection method calculate the three-phase positive sequence voltage (
Figure FDA0000029189850000037
With
Figure FDA0000029189850000038
) step as follows:
Three phase network voltage (u a, u b, u c) through three or two coordinate transforms, be converted into dq reference axis direct current signal u dAnd u q, the three phase network voltage transitions is that dq axial coordinate computing formula is as follows:
u d u q = 2 3 sin wt sin ( wt - 2 &pi; / 3 ) sin ( wt + 2 &pi; / 3 ) - cos wt - cos ( wt - 2 &pi; / 3 ) - cos ( wt + 2 &pi; / 3 ) u a u b u c - - - ( 7 )
Wherein w is meant the first-harmonic angular frequency, the direct current signal u under the above-mentioned dq reference axis that obtains dAnd u q, after low-pass filtering, inverse transformation transforms, and is reduced to the three-phase positive sequence voltage
Figure FDA00000291898500000310
With
Figure FDA00000291898500000311
Computing formula is as follows:
u a + u b + u c + = 2 3 sin wt - cos wt sin ( wt - 2 &pi; / 3 ) - cos ( wt - 2 &pi; / 3 ) sin ( wt + 2 &pi; / 3 ) - cos ( wt + 2 &pi; / 3 ) u d + u q + . - - - ( 8 )
4. according to the described SVC composite control method that calculates based on quick equivalent susceptance of claim 1, it is characterized in that, in the step d), by total output B of compound control AbcThe actual needs that draws SVC drops into the actual equivalent susceptance value that PPF organizes number n and TCR Power distribution method be:
If compensation capacity of each group passive filter is Q, then the passive filter group of input compensation is counted n and should be satisfied following condition:
max { B Tab , B Tbc , B Tca } < min | n Q 3 U 2 | - - - ( 9 )
Wherein U is the effective value of phase voltage, B Tab, B Tbc, B TcaThe total three-phase of exporting of compound control that calculates for step c) in the claim 1 compensates susceptance; After dropping into n group passive filter, because three-phase power supply system is in capacitive and crosses the benefit state, idle to reach idle and full remuneration negative phase-sequence by the corresponding perception of TCR compensation, then the perceptual idle susceptance of TCR actual needs input is:
B ab * = B Tab - nQ 3 U 2 B bc * = B Tbc - nQ 3 U 2 B ca * = B Tca - nQ 3 U 2 . - - - ( 10 )
5. according to the described SVC composite control method that calculates based on quick equivalent susceptance of claim 1, it is characterized in that in the step d), TCR is according to actual equivalent susceptance value By searching the specific implementation step that Hash shows to obtain relative trigger angle α be:
Because be a kind of non-linear relation between the equivalent susceptance that trigger angle α and TCR need compensate, as follows:
B = 2 &pi; - 2 &alpha; - sin ( 2 &pi; - 2 &alpha; ) &pi;X - - - ( 11 )
Wherein X is the reactance value of connecting with thyristor, and α is the trigger angle of TCR; Therefore at first the following formula linearity is changed into the Hash question blank, because the excursion of the trigger angle α of TCR is [90 °, 180 °], trigger angle α is since 90 °, with 0.5 ° as incremental change, each trigger angle α can calculate an equivalent susceptance value, (α B), can calculate 180 groups of vectors altogether so just to constitute a vector, deposit in the corresponding Hash question blank, so just finished the linearization process of equivalent susceptance B to trigger angle α; And then according to actual equivalent susceptance value
Figure FDA0000029189850000046
Obtain relative trigger angle α by searching described Hash question blank.
CN2010105166819A 2010-10-22 2010-10-22 SVC composite control method based on rapid equivalent susceptance calculation Expired - Fee Related CN101950972B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010105166819A CN101950972B (en) 2010-10-22 2010-10-22 SVC composite control method based on rapid equivalent susceptance calculation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010105166819A CN101950972B (en) 2010-10-22 2010-10-22 SVC composite control method based on rapid equivalent susceptance calculation

Publications (2)

Publication Number Publication Date
CN101950972A true CN101950972A (en) 2011-01-19
CN101950972B CN101950972B (en) 2012-05-30

Family

ID=43454365

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010105166819A Expired - Fee Related CN101950972B (en) 2010-10-22 2010-10-22 SVC composite control method based on rapid equivalent susceptance calculation

Country Status (1)

Country Link
CN (1) CN101950972B (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102570503A (en) * 2011-12-30 2012-07-11 东方电气集团东方汽轮机有限公司 Double-fed wind power generation system
CN103138275A (en) * 2013-02-01 2013-06-05 中国科学院等离子体物理研究所 Switching virtual circuit (SVC) feedforward feedback controlling method based on high-power four-quadrant converter load
CN103472282A (en) * 2013-09-23 2013-12-25 东南大学 Improved FBD harmonic current detection method based on adaptive principle
CN104541428A (en) * 2012-06-07 2015-04-22 伽柏·法卡斯 Control method, control device and mobile electric power storage apparatus
CN104810830A (en) * 2015-05-04 2015-07-29 上海电力学院 Three-phase four-wire system four-leg static synchronous compensator and control method
CN105162138A (en) * 2015-09-11 2015-12-16 合肥工业大学 Reactive and harmonic current rapid detection method based on voltage sequence decomposition
CN106291170A (en) * 2016-07-27 2017-01-04 中国科学院等离子体物理研究所 TCR type SVC dynamic response performance test method towards Large Copacity impact load
CN106602576A (en) * 2017-03-03 2017-04-26 国家电网公司 Voltage reactive control system and control method based on static reactive voltage compensation
CN106771507A (en) * 2017-01-20 2017-05-31 泰州学院 Based on the synchronous Fast reactive current of voltage reference split-phase
CN110661247A (en) * 2019-11-12 2020-01-07 湖南大学 Power coefficient compensation-based power equalization control method and system for direct-current micro-grid
CN112072648A (en) * 2020-08-28 2020-12-11 武汉大学 Method for judging optimal access point of electric energy quality control device for inhibiting inter-harmonic resonance

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1766774A (en) * 2005-09-13 2006-05-03 中国电力科学研究院 Method for regulating static var compensator of electricity transmission system
CN101162841A (en) * 2007-10-18 2008-04-16 湖南大学 Non-linearity PI voltage control method of stillness reactive compensator
CN101183791A (en) * 2007-12-19 2008-05-21 湖南大学 Static reactive compensator and active power filter combined operation system and control method thereof
CN101202448A (en) * 2007-12-19 2008-06-18 湖南大学 Harmonic based on APF and SVC and reactive-load dynamic state compensating system and frequency-dividing phase-splitting electrical current control method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1766774A (en) * 2005-09-13 2006-05-03 中国电力科学研究院 Method for regulating static var compensator of electricity transmission system
CN101162841A (en) * 2007-10-18 2008-04-16 湖南大学 Non-linearity PI voltage control method of stillness reactive compensator
CN101183791A (en) * 2007-12-19 2008-05-21 湖南大学 Static reactive compensator and active power filter combined operation system and control method thereof
CN101202448A (en) * 2007-12-19 2008-06-18 湖南大学 Harmonic based on APF and SVC and reactive-load dynamic state compensating system and frequency-dividing phase-splitting electrical current control method

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102570503A (en) * 2011-12-30 2012-07-11 东方电气集团东方汽轮机有限公司 Double-fed wind power generation system
CN104541428A (en) * 2012-06-07 2015-04-22 伽柏·法卡斯 Control method, control device and mobile electric power storage apparatus
CN103138275A (en) * 2013-02-01 2013-06-05 中国科学院等离子体物理研究所 Switching virtual circuit (SVC) feedforward feedback controlling method based on high-power four-quadrant converter load
CN103138275B (en) * 2013-02-01 2015-04-15 中国科学院等离子体物理研究所 Switching virtual circuit (SVC) feedforward feedback controlling method based on high-power four-quadrant converter load
CN103472282A (en) * 2013-09-23 2013-12-25 东南大学 Improved FBD harmonic current detection method based on adaptive principle
CN103472282B (en) * 2013-09-23 2015-09-16 东南大学 A kind of FBD harmonic current detecting method based on adaptive principle
CN104810830A (en) * 2015-05-04 2015-07-29 上海电力学院 Three-phase four-wire system four-leg static synchronous compensator and control method
CN105162138A (en) * 2015-09-11 2015-12-16 合肥工业大学 Reactive and harmonic current rapid detection method based on voltage sequence decomposition
CN106291170A (en) * 2016-07-27 2017-01-04 中国科学院等离子体物理研究所 TCR type SVC dynamic response performance test method towards Large Copacity impact load
CN106291170B (en) * 2016-07-27 2020-09-29 中国科学院等离子体物理研究所 TCR type SVC dynamic response performance test method for high-capacity impact load
CN106771507A (en) * 2017-01-20 2017-05-31 泰州学院 Based on the synchronous Fast reactive current of voltage reference split-phase
CN106771507B (en) * 2017-01-20 2020-02-14 泰州学院 Reactive current rapid detection method based on voltage reference split-phase synchronization
CN106602576A (en) * 2017-03-03 2017-04-26 国家电网公司 Voltage reactive control system and control method based on static reactive voltage compensation
CN106602576B (en) * 2017-03-03 2019-09-13 国家电网公司 Voltage powerless control system and control method based on static reactive voltage compensation
CN110661247A (en) * 2019-11-12 2020-01-07 湖南大学 Power coefficient compensation-based power equalization control method and system for direct-current micro-grid
CN110661247B (en) * 2019-11-12 2021-03-23 湖南大学 Power coefficient compensation-based power equalization control method and system for direct-current micro-grid
CN112072648A (en) * 2020-08-28 2020-12-11 武汉大学 Method for judging optimal access point of electric energy quality control device for inhibiting inter-harmonic resonance

Also Published As

Publication number Publication date
CN101950972B (en) 2012-05-30

Similar Documents

Publication Publication Date Title
CN101950972B (en) SVC composite control method based on rapid equivalent susceptance calculation
CN101924371B (en) Mixed power quality treatment method
US20070014132A1 (en) Comprehensive power quality controller for substation in power system
CN102157948B (en) Static var compensation device current detection method based on transient power balance
CN102638049A (en) Direct-current bus inter-phase voltage balancing control method for chained type triangular connection STATCOM (Static Synchronous Compensator)
CN103956734A (en) Multi-target electric energy quality comprehensive control and optimization device
CN107769230A (en) A kind of active and passive combination power load three-phase imbalance adjusting means
CN104466968B (en) A kind of DSTATCOM negative-sequence currents cross coupling compensation control method
CN207625296U (en) A kind of active and passive combination power load three-phase imbalance regulating device
CN102570474B (en) Triggering and conducting method used for switching of capacitor of thyristor AC non-contact switch
CN101969204A (en) Vector control-based static volt-ampere reactive compensation (SVC) method and device
Kumbha et al. Power quality improvement of distribution lines using DSTATCOM under various loading conditions
Hou et al. Adaptive active voltage-type arc suppression strategy considering the influence of line parameters in active distribution network
Li et al. Research on current-limiting control strategy suitable for ground faults in AC microgrid
CN111917132A (en) Method for improving robustness of multi-inverter parallel low-voltage microgrid sag control system
CN204465011U (en) Distribution power consumption dynamic optimization bascule
Wang et al. Research on application of TCR+ FC typed SVC in power quality integrated management for power traction system
CN106291170A (en) TCR type SVC dynamic response performance test method towards Large Copacity impact load
Zheng et al. The method of harmonic compensation of UPQC based on Euler-Lagrange model
Dexin et al. Research on the method of using apf to improve the power quality of microgrid based on the control of energy storage
CN107017643A (en) A kind of static reacance generator of dead-beat current control
Ma et al. Research on reactive power compensation technology in distribution grid
Jin et al. Model predictive control for three-phase unbalance compensation under three-phase four-wire power system
Xu An improved current-limiting control strategy for shunt active power filter
Xiao et al. Research of power quality management technology according to distribution network involving electric arc furnace

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
ASS Succession or assignment of patent right

Owner name: BAIYUN ELECTRIC GROUP CO., LTD.

Effective date: 20130624

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20130624

Address after: 410082 Hunan province Changsha Lushan Road No. 2

Patentee after: Hunan University

Patentee after: Baiyun Power Group Co., Ltd.

Address before: 410082 Hunan province Changsha Lushan Road No. 2

Patentee before: Hunan University

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120530

Termination date: 20201022