CN104779635A - Controller device suitable for VSC-MTDC (voltage source converter-multi-terminal direct current) system - Google Patents
Controller device suitable for VSC-MTDC (voltage source converter-multi-terminal direct current) system Download PDFInfo
- Publication number
- CN104779635A CN104779635A CN201510137472.6A CN201510137472A CN104779635A CN 104779635 A CN104779635 A CN 104779635A CN 201510137472 A CN201510137472 A CN 201510137472A CN 104779635 A CN104779635 A CN 104779635A
- Authority
- CN
- China
- Prior art keywords
- phase
- current
- power
- axle
- represent
- 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
- 230000005540 biological transmission Effects 0.000 claims abstract description 19
- 230000003068 static effect Effects 0.000 claims description 20
- 230000001105 regulatory effect Effects 0.000 claims description 11
- 230000009466 transformation Effects 0.000 claims description 11
- 230000001052 transient effect Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 5
- 230000001276 controlling effect Effects 0.000 claims description 4
- WCUXLLCKKVVCTQ-UHFFFAOYSA-M Potassium chloride Chemical compound [Cl-].[K+] WCUXLLCKKVVCTQ-UHFFFAOYSA-M 0.000 claims description 3
- 230000007935 neutral effect Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 3
- 238000013178 mathematical model Methods 0.000 abstract description 3
- 230000010354 integration Effects 0.000 abstract 1
- 238000013461 design Methods 0.000 description 8
- 230000008859 change Effects 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000011217 control strategy Methods 0.000 description 3
- 238000000034 method Methods 0.000 description 3
- 230000033228 biological regulation Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000004870 electrical engineering Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/36—Arrangements for transfer of electric power between ac networks via a high-tension dc link
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/60—Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Inverter Devices (AREA)
Abstract
The invention provides a controller device suitable for a VSC-MTDC (voltage source converter-multi-terminal direct current) system. The controller device comprises a rectifier-side local controller used for controlling a current transformer station with a rectification effect in an MTDC power transmission system, an inverter-side local controller used for controlling a current transformer station with an inversion effect in the MTDC power transmission system, and a coordination controller used for coordinating electric power physical parameters such as power, current and the like of each current transformer station in the MTDC power transmission system. The controller device effectively avoids PI (proportional integration) parameter setting and calculation of conventional double closed loop control and greatly lowers the complexity of control on rectification and inversion of current transformers without a linear controller and PWM (pulse width modulation); at the same time, a discrete mathematical model is fully utilized based on a P-DPC (predictive-direct power control) strategy, calculation is simple, and digitization is easy to realize.
Description
Technical field
The present invention relates to electrical engineering field, specifically one is applicable to the controller of VSC-MTDC (multi-terminal HVDC transmission based on voltage source converter) system, is especially applicable to the control device of VSC-MTDC system based on prediction-Direct Power (P-DPC).
Background technology
Along with the approach exhaustion of traditional energy, new forms of energy more and more receive the concern of people, and offshore wind farm then gains great popularity because of its plurality of advantages.At present, the long-distance transmissions of Large Scale Offshore Wind Farm adopts the system based on VSC-HVDC usually.On the basis of both-end high voltage direct current transmission, multi-terminal HVDC transmission MTDC technology is improved gradually and is applied.Multi-terminal direct current transmission system at least comprises the rotary substation of more than 3 or 3.Due to developing rapidly of electric power transfer, traditional two ends HVDC can not meet the demands gradually, and multi-terminal HVDC transmission is more and more subject to people and payes attention to.Compared with both-end high-voltage direct current, multi-terminal direct current transmission system is except needing to consider the control of each rotary substation self, also to consider the cooperation control between each rotary substation, adopt the stable operation of effective Controller gain variations to VSC-MTDC transmission system most important.Generally speaking, tradition has following shortcoming based on PI controller:
1, the local control design case of current transformer is complicated, comprises voltage, current controller;
2, have multiple PI controller, PI parameter regulates complicated;
3, cooperation control cannot carry out operational mode handoff functionality, and system reliability is low.
Summary of the invention
For defect of the prior art, the object of the invention is to overcome above proposed shortcoming and defect, propose a kind of control device being applicable to VSC-MTDC system.The present invention is directed to VSC-MTDC system, this locality of current transformer is controlled, devise rectification side local controller and inverter side local controller respectively; Devise the tuning controller between rotary substation simultaneously.
Rectification side local controller, for the control of rotary substation playing rectified action in multi-terminal direct current transmission system.
Inverter side local controller, for the control of rotary substation playing reversion reaction in multi-terminal direct current transmission system.
Tuning controller, for the controller coordinated electric power physical quantitys such as the power in multi-terminal direct current transmission system between each rotary substation, electric currents.
According to a kind of control device being applicable to VSC-MTDC system provided by the invention, comprise model apparatus for establishing;
Described model apparatus for establishing, for setting up three-phase VSC model, specific as follows:
Three-phase VSC comprises IGBT switch S 1, S2, S3, S4, S5, S6, also comprises DC power supply;
Be parallel with A phase brachium pontis, B phase brachium pontis, C phase brachium pontis between DC power supply, wherein, IGBT switch S 1, S4 form A phase brachium pontis, and IGBT switch S 2, S5 form B phase brachium pontis, and IGBT switch S 3, S6 form C phase brachium pontis;
The IGBT switch function of IGBT switch S 1, S2, S3, S4, S5, S6 is respectively S
1, S
2, S
3, S
4, S
5, S
6;
Brachium pontis on off state is defined as follows:
S
arepresent A phase brachium pontis switch function, S
brepresent B phase brachium pontis switch function, S
crepresent C phase brachium pontis switch function; Then the resultant vector S of A phase, B phase, C phase brachium pontis switch function is:
S=S
a+ α S
b+ α
2s
c(4) wherein, α=e
j2 π/3, e is natural constant;
Calculate three-phase VSC output voltage vector U
i, calculating formula is:
U
i=SU
dci=0,1,2,…,7 (5)
U
irepresent i-th three-phase VSC output voltage vector, U
dcrepresent the direct voltage of DC power supply;
Obtain through three phase static abc coordinate system and two-phase static α β coordinate system transformation:
U
αrepresent the three-phase VSC output voltage along α axle in two-phase static α β coordinate system, U
βrepresent the three-phase VSC output voltage along β axle in two-phase static α β coordinate system;
Following transient current equation is obtained by Kirchhoff's current law (KCL):
I
a, i
b, i
crepresent A phase, B phase, C phase current respectively, t represents the time, U
aN, U
bN, U
cNrepresent the voltage of A phase, B phase, the relative neutral point N of C respectively, e
a, e
b, e
crepresent A phase, B phase, the C phase voltage of three-phase alternating current system respectively, R represents resistance, and L represents reactance;
Suppose three-phase grid balance, above-mentioned transient current equation obtains following equation through three phase static abc coordinate system and two-phase static α β coordinate system transformation:
E
αrepresent the three-phase alternating current system voltage component along α axle in two-phase static α β coordinate system, e
βrepresent the three-phase alternating current system voltage component along β axle in two-phase static α β coordinate system, i
αrepresent the electric current along α axle in two-phase static α β coordinate system, i
βrepresent the electric current along β axle in two-phase static α β coordinate system.
Preferably, inverter side controller is also comprised;
Described inverter side controller, for realizing following function:
By formula (8) discretization, obtain following equation:
K represents current time, and k+1 represents subsequent time, i
α(k+1) electric current of subsequent time along α axle is represented, i
β(k+1) electric current of subsequent time along β axle is represented, T
srepresent the sampling time, i
αk () represents the electric current of current time along α axle, i
βk () represents the electric current of current time along β axle, U
αk () represents the output voltage of current time along the three-phase VSC of α axle, U
βk () represents the output voltage of current time along the three-phase VSC of β axle, e
αk () represents the three-phase alternating current system voltage component of current time along α axle, e
βk () represents the three-phase alternating current system voltage component of current time along β axle;
The equation under dq two-phase rotating coordinate system can be obtained through rotation transformation to formula (9):
I
d(k+1) the three-phase grid-connected inverter output current d axle component that the k moment dopes is represented, i
q(k+1) represent the three-phase grid-connected inverter output current q axle component that the k moment dopes, θ represents the space angle of electrical network;
Active power and the reactive power in (k+1) moment is drawn by following prediction equation:
P (k+1) represents the active power in (k+1) moment, and Q (k+1) represents the reactive power in (k+1) moment, e
drepresent the three-phase alternating current system voltage component of d axle, e
qrepresent the three-phase alternating current system voltage component of q axle;
Calculate the first cost function g, calculating formula is:
g=|(Q
*-Q(k+1)|+|(P
*-P(k+1)| (12)
P
*for given active power reference value, Q
*for given reactive power reference qref;
Calculate respectively at 8 three-phase VSC output voltage vector U
i8 the first cost function g under state, obtain the minimum value g of the first cost function g
min; Obtain the minimum value g corresponding to the first cost function g
minclosest to expecting power, thus obtain corresponding U
α(k), U
βk the size of (), by U
α(k), U
βk () substitutes into the U in formula (6) respectively
α, U
β, obtain the switch function that this k moment expects, and the switch function of this expectation passed to the brachium pontis on off state of subsequent time, and then control inverter side current transformer; Wherein, described closest to expecting that power is the minimum performance number Q (k+1) of the value of instruction first cost function g and P (k+1), i=0,1,2 ... 7.
Preferably, rectification side controller is also comprised;
Described rectification side controller, for realizing following function:
Set up rectification side three-phase VSC model:
Discretization formula (13) obtains predicted current function:
K represents current time, and k+1 represents subsequent time, i
α(k+1) electric current of subsequent time along α axle is represented, i
β(k+1) electric current of subsequent time along β axle is represented, T
srepresent the sampling time, i
αk () represents the electric current of current time along α axle, i
βk () represents the electric current of current time along β axle, U
αk () represents the output voltage of current time along the three-phase VSC of α axle, U
βk () represents the output voltage of current time along the three-phase VSC of β axle, e
αk () represents the three-phase alternating current system voltage component of current time along α axle, e
βk () represents the three-phase alternating current system voltage component of current time along β axle;
The equation under dq two-phase rotating coordinate system can be obtained through rotation transformation to formula (14):
Active power and the reactive power in (k+1) moment is drawn by following prediction equation:
P (k+1) represents the active power in (k+1) moment, and Q (k+1) represents the reactive power in (k+1) moment, e
drepresent the three-phase alternating current system voltage component of d axle, e
qrepresent the three-phase alternating current system voltage component of q axle; i
d(k+1) the three-phase grid-connected inverter output current d axle component that the k moment dopes is represented, i
q(k+1) the three-phase grid-connected inverter output current q axle component that the k moment dopes is represented; θ represents the space angle of electrical network;
Calculate the second cost function g, calculating formula is:
g=|(Q
*-Q(k+1)|+|(P
dc+P
i-P(k+1))| (15)
Q
*for given reactive power reference qref, P
dcrepresent the direct current power after regulating, P
irepresent the active power that i-th other rotary substation consumes except rectification side current transformer in multi-terminal direct current transmission system, i>0;
Calculate respectively at the active-power P that i other rotary substation consumes
ii under state the second cost function g, obtains the minimum value g of the second cost function g
min; Obtain the minimum value g corresponding to the second cost function g
minclosest to expecting power, thus obtain corresponding U
α(k), U
βk the size of (), by U
α(k), U
βk () substitutes into the U in formula (6) respectively
α, U
β, obtain the switch function that this k moment expects, and the switch function of this expectation passed to the brachium pontis on off state of subsequent time, and then control rectification side current transformer; Wherein, described closest to expecting that power is the minimum performance number Q (k+1) of the value of instruction second cost function g and P (k+1), i>0.
Preferably, also tuning controller is comprised;
Described tuning controller, for controlling main convertor and there is the current transformer of power adjustments, specific as follows:
Work as P
1(k+1) allow in adjustable range at power, i.e. P
1min<P
1(k+1) <P
1max, then under making VSC2, VSC3 be operated in power mode; When VSC1 power adjustments is not enough or out of service, P
1(k+1) adjustable range [P is not allowed at power
1min, P
1max] in, then make the operational mode of VSC2 switch to direct voltage pattern to bear power shortage, under VSC3 is operated in power mode by power mode;
Wherein: VSC1 represents main convertor, VSC2 is the current transformer with power adjustments, and VSC3 is the current transformer without regulating power; P
1minrepresent that rotary substation is exerted oneself lower limit, i.e. minimum load, P
1(k+1) value of current predictive power is represented, P
1maxrepresent that rotary substation is exerted oneself the upper limit, i.e. maximum output.
Compared with prior art, the present invention has following beneficial effect:
(1) local controller method for designing is simple, and principle is also uncomplicated, without the need to voltage and current controller, does not also need PWM.
(2) there is no PI controller in local controller, without the need to PI calculation of parameter and the adjustment of complexity, greatly improve control efficiency.
(3) tuning controller designed by has the function switching current transformer mode of operation, greatly improves the stability of a system.
(4) take full advantage of discrete mathematical model, calculate simple, easy Digital Realization.
Accompanying drawing explanation
By reading the detailed description done non-limiting example with reference to the following drawings, other features, objects and advantages of the present invention will become more obvious:
Fig. 1 is three-phase VSC combining inverter.
Fig. 2 is the design of inverter side local controller.
Fig. 3 is the design of rectification side local controller.
Fig. 4 is the tuning controller based on prediction-direct Power Control.
Fig. 5 is VSC-MTDC total system simulation model.
Fig. 6 is normal operating conditions DC voltage waveform.
Fig. 7 is each current transformer active power waveform under normal operating conditions.
Fig. 8 is the DC voltage waveform under active power sudden change exceeds VSC1 adjustable range.
Fig. 9 is each current transformer active power waveform under active power sudden change exceeds VSC1 adjustable range.
Figure 10 is the DC voltage waveform of VSC3 rotary substation fault when exiting in short-term.
Figure 11 is each current transformer active power waveform of VSC3 rotary substation fault when exiting in short-term.
Embodiment
Below in conjunction with specific embodiment, the present invention is described in detail.Following examples will contribute to those skilled in the art and understand the present invention further, but not limit the present invention in any form.It should be pointed out that to those skilled in the art, without departing from the inventive concept of the premise, some distortion and improvement can also be made.These all belong to protection scope of the present invention.
The control device being applicable to VSC-MTDC system provided by the invention, comprising: model apparatus for establishing, rectification side local controller, inverter side local controller, tuning controller;
Described model apparatus for establishing, for setting up three-phase VSC Mathematical Modeling, as shown in Figure 1, comprises 6 IGBT switches.Wherein, have three brachium pontis, upper and lower each 1 the IGBT switch of each brachium pontis.
Three-phase VSC comprises IGBT switch S 1, S2, S3, S4, S5, S6, also comprises DC power supply;
Be parallel with A phase brachium pontis, B phase brachium pontis, C phase brachium pontis between DC power supply, wherein, IGBT switch S 1, S4 form A phase brachium pontis, and IGBT switch S 2, S5 form B phase brachium pontis, and IGBT switch S 3, S6 form C phase brachium pontis;
The IGBT switch function of IGBT switch S 1, S2, S3, S4, S5, S6 is respectively S
1, S
2, S
3, S
4, S
5, S
6;
In Fig. 1, U
dcrepresent direct voltage, a represents A phase, and b represents B phase, and c represents C phase, i
arepresent A phase current, i
brepresent B phase current, i
crepresent C phase current, e
a, e
b, e
crepresent each phase voltage of three-phase alternating current system, N represents neutral point;
Brachium pontis on off state is defined as follows:
S
arepresent A phase brachium pontis switch function, S
brepresent B phase brachium pontis switch function, S
crepresent C phase brachium pontis switch function;
Then the resultant vector S of above-mentioned switch function is:
S=S
a+ α S
b+ α
2s
c(4) wherein, α=e
j2 π/3, e is constant, and be a nonterminating and non-recurring decimal, its value approximates 2.718281828459.
Because switch combination number has 8 kinds, the output voltage vector of three-phase VSC can calculate:
U
i=SU
dci=0,1,2,…,7 (5)
U
irepresent i-th three-phase VSC output voltage vector, U
dcrepresent direct voltage;
Calculate the output voltage of three-phase VSC, obtain through abc/ α β coordinate transform:
U
αrepresent the output voltage along the VSC of α axle, U
βrepresent the output voltage along the VSC of β axle;
Following transient current equation is obtained by Kirchhoff's current law (KCL):
I
a, i
b, i
crepresent three-phase current, t represents the time, U
aN, U
bN, U
cNrepresent the voltage of three-phase centering point N, e
a, e
b, e
crepresent each phase voltage of three-phase alternating current system, R represents resistance, and L represents reactance;
Suppose three-phase grid balance, above-mentioned transient current equation obtains following equation through coordinate transform:
E
αrepresent the AC system component of voltage along α axle, e
βrepresent the AC system component of voltage along β axle, i
αrepresent the electric current along α axle, i
βrepresent the electric current along β axle;
Described inverter side controller, for realizing following function:
By discretization three-phase VSC Mathematical Modeling, introduce prediction-Direct Power (P-DPC) control strategy.By formula (8) discretization, obtain following equation:
K represents current time, i
α(k+1) electric current of subsequent time along α axle is represented, i
β(k+1) electric current of subsequent time β axle is represented, T
srepresent the sampling time, i
αk () represents the electric current of current time along α axle, i
βk () represents the electric current of current time β axle, U
αk () represents the output voltage of the VSC of current time α axle, U
βk () represents the output voltage of the VSC of current time β axle, e
αk () represents the AC system component of voltage of current time α axle, e
βk () represents the AC system component of voltage of current time β axle;
From formula (9), after discretization inverter Mathematical Modeling, contacting between the electric current of subsequent time and a upper moment electric current can be obtained, namely can go out the electric current in (k+1) moment from the current forecasting in k moment.
Equation under rotation transformation can obtain d, q rotational coordinates:
I
d(k+1) the three-phase grid-connected inverter output current d axle component that the k moment dopes is represented, i
q(k+1) represent the three-phase grid-connected inverter output current q axle component that the k moment dopes, θ represents the space angle of electrical network;
Suppose the stabilization of power grids, then three-phase voltage can be regarded as substantially constant.Active power and the reactive power in (k+1) moment so just can be drawn by following prediction equation.
P (k+1) represents the active power in (k+1) moment, and Q (k+1) represents the reactive power in (k+1) moment, e
drepresent the AC system component of voltage of d axle, e
qrepresent the AC system component of voltage of q axle.
Next is selected cost function.Valency value function minimum principle and power-balance, adopt the first cost function g:
g=|(Q
*-Q(k+1)|+|(P
*-P(k+1)| (12)
Wherein P
*with Q
*for given active power and reactive power reference qref.Utilize the first cost function formula (12), pass through global optimizing, instead can release the size closest to expecting power, closest to expecting that power is the minimum performance number Q (k+1) of instruction type (12) value and P (k+1), thus the voltage U under this state can be obtained
α(k), U
β(k) size, more most suitable on off state can be obtained by formula (6).
Be specially, calculate respectively at 8 three-phase VSC output voltage vector U
i8 the first cost function g under state, obtain the minimum value g of the first cost function g
min; Obtain the minimum value g corresponding to the first cost function g
minclosest to expecting power, thus obtain corresponding U
α(k), U
βthe size of (k), then through type (6) provides the brachium pontis on off state of subsequent time, and then control inverter side current transformer; Wherein, described closest to expecting that power is the minimum performance number Q (k+1) of the value of instruction first cost function g and P (k+1), i=0,1,2 ... 7.
Whole process is and provides switching vector selector by prediction-direct Power Control.Whole controller architecture design as shown in Figure 2.
Described rectification side controller, for realizing following function:
Rectification side current transformer assume responsibility for stable DC voltage and provides the function of system power, and therefore, rectification side current transformer must adopt constant DC voltage control pattern, controls direct voltage and reactive power.Due to the structural similarity with inverter side current transformer, be not difficult to obtain rectification side VSC Mathematical Modeling:
Owing to thinking DC voltage stability, U
α, U
βstill can be calculated by formula (6).Discretization formula (13) can obtain predicted current function:
Again through d, q rotation transformation and mathematical computations, the active power and the reactive power that input VSC can be obtained, avatar still can be represented by formula (11), is specially: can to obtain the equation under dq two-phase rotating coordinate system through rotation transformation to formula (14):
Active power and the reactive power in (k+1) moment is drawn by following prediction equation:
Because rectification side current transformer must bear the function of regulating power, so direct current power governing loop need be added at controller.Choosing for cost function, whole system must could reliability service under power-balance, adopts following cost function:
g=|(Q
*-Q(k+1)|+|(P
dc+P
i-P(k+1))| (15)
In formula (15), P
dcfor the direct current power after adjustment, P
ifor the active power that other rotary substations except rectification side current transformer consume, i represents number, and the span of i is i>0.Each P
isize-dependent in communication system faster.Whole process is still the vector control based on prediction-Direct Power, and control design case block diagram as shown in Figure 3.
Particularly, calculate respectively at the active-power P that i other rotary substation consumes
ii under state the second cost function g, obtains the minimum value g of the second cost function g
min; Obtain the minimum value g corresponding to the second cost function g
minclosest to expecting power, thus obtain corresponding U
α(k), U
βthe size of (k), then through type (6) provides the brachium pontis on off state of subsequent time, and then control rectification side current transformer; Wherein, described closest to expecting that power is the minimum performance number Q (k+1) of the value of instruction second cost function g and P (k+1), i>0.
In Fig. 3,4, U
* dcrepresent direct voltage reference value, Pref represents active power reference value, P
minrepresent rotary substation power output lower limit, i.e. minimum load, P
maxrepresent the rotary substation power output upper limit, i.e. maximum output.
Described tuning controller, for controlling main convertor and there is the current transformer of power adjustments, specific as follows:
Work as P
1(k+1) allow in adjustable range at power, i.e. P
1min<P
1(k+1) <P
1max, then under making VSC2, VSC3 be operated in power mode; When VSC1 power adjustments is not enough or out of service, P
1(k+1) adjustable range [P is not allowed at power
1min, P
1max] in, then make the operational mode of VSC2 switch to direct voltage pattern to bear power shortage, under VSC3 is operated in power mode by power mode;
Wherein: VSC1 represents main convertor, VSC2 is the current transformer with power adjustments, and VSC3 is the current transformer without regulating power; P
1minrepresent that rotary substation is exerted oneself lower limit, i.e. minimum load, P
1(k+1) value of current predictive power is represented, P
1maxrepresent that rotary substation is exerted oneself the upper limit, i.e. maximum output.
Particularly, design at tuning controller, if only specify a current transformer to adopt constant DC voltage control in multi-terminal system, system operation reliability difference and active power can be caused accurately not to control, therefore need to adopt pattern switching controls to some current transformers, when main convertor cannot meet power adjustments time, the current transformer with power regulation functions can bear power shortage automatically.Be implemented as follows, the power that each current transformer is born is transferred to main convertor by communication system and has the current transformer of power adjustments, change main convertor by prediction-direct Power Control or there is the control signal of current transformer of power adjustments, the operational mode changing current transformer and then the power controlling current transformer are born, and ensure system power balance.As shown in Figure 3, be designed tuning controller, this controller be applied to main convertor and there is the current transformer of power adjustments.As shown in Figure 4, when being operated in normal condition, P
1min<P
1(k+1) <P
1max, switch is in 1 state, under VSC2 is operated in power mode; When VSC1 power adjustments is not enough or out of service, P
1(k+1) do not allow in adjustable range at power, switch is in 2 states, and VSC2 operational mode switches to direct voltage pattern.P
1minrepresent that rotary substation is exerted oneself lower limit, i.e. minimum load, P
1(k+1) value of current predictive power is represented, P
1maxrepresent that rotary substation is exerted oneself the upper limit, i.e. maximum output; Wherein: VSC1 represents main convertor, VSC2 is the current transformer with power adjustments of specifying.
Tuning controller is operated in following three kinds of operating states:
(1) normal operating conditions, the predicted power of each current transformer is transferred to main convertor by communication system, utilize formula (6) and cost function formula (15), provide the control gate signal of main convertor under normal operating conditions, realize prediction-direct Power Control, ensure that each current transformer Direct Power balances; Now VSC1 is operated in and determines direct voltage pattern.VSC2, VSC3 are operated in power mode.Normal conditions, system works is in this state.VSC3 represents the current transformer without regulating power.
(2) VSC3 chugging, the power output of VSC1 current transformer exceeds its freely regulated allowed band [P
1min, P
1max], exceed power section and bear by current transformer VSC2 with power adjustments ability all in direct current system.Now, VSC1 keeps maximum power output or minimum power, is operated in constant DC voltage control pattern.Based on P-DPC cooperation control, VSC2 is switched to direct voltage pattern by power mode, automatically bears power shortage, makes whole system again keep power-balance.VSC3 is still operated in power mode.
(3) VSC1 main convertor is out of service, and VSC2 bears power shortage completely, and be switched to direct voltage pattern by VSC2 by power mode, VSC3 is still operated in power mode.Wherein: VSC1 represents main convertor, VSC2 is the current transformer with power adjustments of specifying, and VSC3 is the current transformer without power regulation functions.
Further, as shown in Figure 2, net side three-phase voltage and electric current is gathered in the k moment, through type (6), (9), (10), (11) calculate the meritorious of (k+1) moment and reactive power, compare with reactive power reference qref with given active power, by cost function formula (12) and formula (6), global optimizing, select the switching signal vector of the most applicable subsequent time, and then control inverter side current transformer.Wherein, with active power predicted value and the absolute value of difference of active power reference value and the absolute value sum of the difference of reactive power predicted value and reactive power reference qref for being worth reference index, its concrete form is as follows:
g=|(Q
*-Q(k+1)|+|(P
*-P(k+1)|
As shown in Figure 3, rectification side local controller designs, net side three-phase voltage and electric current is gathered in the k moment, through type (6), (14), (10), (11) calculate the meritorious of (K+1) moment and reactive power, with by communication system transmission come each current transformer power sum and modulation direct current power compared with, by cost function formula (15) and formula (6), global optimizing, select the switching signal vector of the most applicable subsequent time, and then control rectification side current transformer.Wherein, with active power predicted value and the absolute value of difference of all current transformer active power and the direct current power summation after modulating and the absolute value sum of the difference of reactive power predicted value and system reactive power reference value for being worth reference index, concrete form is as follows:
g=|(Q
*-Q(k+1)|+|(P
dc+P
i-P(k+1))|
Wherein, P
dcfor the direct current power after adjustment, P
ifor the active power that other rotary substations except rectification side current transformer consume.Each P
isize-dependent in communication system faster.
As shown in Figure 4, tuning controller designs, and is designed tuning controller, this controller is applied to main convertor and has the current transformer of power adjustments.When being operated in normal condition, P
1min<P
1(k+1) <P
1max, switch is in 1 state, under VSC2 is operated in power mode; When VSC1 power adjustments is not enough or out of service, P
1(k+1) do not allow in adjustable range at power, switch is in 2 states, and VSC2 operational mode switches to direct voltage pattern.Wherein: VSC1 represents main convertor, VSC2 is the current transformer with power adjustments of specifying.
Further, the application of controller provided by the invention in multi-terminal direct current transmission system.Fig. 5 is three end DC transmission system simulation models.Wherein VSC1 is main convertor, adopts constant DC voltage control; VSC2 introduces P-DPC cooperation control, and operating state can switch between power mode and direct voltage pattern, namely has power adjustments ability; VSC3 is operated in power mode, not changeable.Concrete simulation parameter is as following table:
Table 1 simulation parameter
The present embodiment will verify the validity of designed controller from following three kinds of operating states.
(1) normal operating conditions, VSC3 rotary substation side active power is suddenlyd change, and VSC1 power output is within adjustable range;
(2) VSC3 rotary substation side active power sudden change, VSC1 power output exceeds adjustable range;
(3) VSC1 rotary substation because of fault out of service in short-term.
In Fig. 5: Pulse1 ~ 3 represent control impuls; V-I 1 ~ 4 represents voltage and current measurement module; RL1 ~ 3 represent line impedance.
System responses is as shown in Figure of description.
Wherein, Fig. 6, Fig. 7 are the system responses under (1) state, can find out when being in normal operating conditions, and be that power or direct voltage are all very stable, wherein, direct voltage is the same with preset value, remains on 750V.When 0.1s ~ 0.2s, current transformer VSC3 power produces fluctuation, and when changing to-6KW by 6KW, the automatic regulating power balance of main variable flow station VSC1, system run all right, illustrates that the control strategy based on P-DPC is effective.In the case, VSC1 is operated in and determines direct voltage pattern, VSC2 and VSC3 is operated in power mode.
Fig. 8, system responses under Fig. 9 correspondence (2) state, can find out, when 0.1s ~ 0.2s, the active power sudden change of VSC3, changes to-1.3MW by original 6kW, cause direct current system unbalanced power, the power injecting DC network is less than its power output, and direct voltage declines, and system works point changes.For checking the validity of controller, VSC2 initial power is set to 2MW here, and now power stage has exceeded the adjustable range [-3MW, 3MW] of main convertor VSC1.VSC1 keeps maximum power 3MW to export, and VSC2 mode of operation automatically switches to direct voltage pattern by power mode, regulating system power-balance.
System responses under Figure 10, Figure 11 correspondence (3) state, 0.1s ~ 0.2s, VSC1 rotary substation out of service in short-term, cause direct current system unbalanced power, direct voltage fluctuates.VSC2 switches to direct voltage pattern by power mode, regulating system power-balance, stable DC voltage.After 0.2s, circuit breaker reclosing, system is recovered normal and is run in 0.05.
The present invention devises current transformer local controller based on P-DPC and tuning controller, effectively prevent PI parameter tuning and the calculating of traditional double closed-loop control, without the need to linear controller and PWM, greatly reduces the control complexity of current transformer rectification and inversion.Meanwhile, the control strategy based on P-DPC takes full advantage of discrete mathematical model, calculates simple, easy Digital Realization.
Above specific embodiments of the invention are described.It is to be appreciated that the present invention is not limited to above-mentioned particular implementation, those skilled in the art can make various distortion or amendment within the scope of the claims, and this does not affect flesh and blood of the present invention.
Claims (4)
1. be applicable to a control device for VSC-MTDC system, it is characterized in that, comprise model apparatus for establishing;
Described model apparatus for establishing, for setting up three-phase VSC model, specific as follows:
Three-phase VSC comprises IGBT switch S 1, S2, S3, S4, S5, S6, also comprises DC power supply;
Be parallel with A phase brachium pontis, B phase brachium pontis, C phase brachium pontis between DC power supply, wherein, IGBT switch S 1, S4 form A phase brachium pontis, and IGBT switch S 2, S5 form B phase brachium pontis, and IGBT switch S 3, S6 form C phase brachium pontis;
The IGBT switch function of IGBT switch S 1, S2, S3, S4, S5, S6 is respectively S
1, S
2, S
3, S
4, S
5, S
6;
Brachium pontis on off state is defined as follows:
S
arepresent A phase brachium pontis switch function, S
brepresent B phase brachium pontis switch function, S
crepresent C phase brachium pontis switch function; Then the resultant vector S of A phase, B phase, C phase brachium pontis switch function is:
S=S
a+αS
b+α
2S
c(4)
Wherein, α=e
j2 π/3, e is natural constant;
Calculate three-phase VSC output voltage vector U
i, calculating formula is:
U
i=SU
dci=0,1,2,…,7 (5)
U
irepresent i-th three-phase VSC output voltage vector, U
dcrepresent the direct voltage of DC power supply;
Obtain through three phase static abc coordinate system and two-phase static α β coordinate system transformation:
U
αrepresent the three-phase VSC output voltage along α axle in two-phase static α β coordinate system, U
βrepresent the three-phase VSC output voltage along β axle in two-phase static α β coordinate system;
Following transient current equation is obtained by Kirchhoff's current law (KCL):
I
a, i
b, i
crepresent A phase, B phase, C phase current respectively, t represents the time, U
aN, U
bN, U
cNrepresent the voltage of A phase, B phase, the relative neutral point N of C respectively, e
a, e
b, e
crepresent A phase, B phase, the C phase voltage of three-phase alternating current system respectively, R represents resistance, and L represents reactance;
Suppose three-phase grid balance, above-mentioned transient current equation obtains following equation through three phase static abc coordinate system and two-phase static α β coordinate system transformation:
E
αrepresent the three-phase alternating current system voltage component along α axle in two-phase static α β coordinate system, e
βrepresent the three-phase alternating current system voltage component along β axle in two-phase static α β coordinate system, i
αrepresent the electric current along α axle in two-phase static α β coordinate system, i
βrepresent the electric current along β axle in two-phase static α β coordinate system.
2. the control device being applicable to VSC-MTDC system according to claim 1, is characterized in that, also comprises inverter side controller;
Described inverter side controller, for realizing following function:
By formula (8) discretization, obtain following equation:
K represents current time, and k+1 represents subsequent time, i
α(k+1) electric current of subsequent time along α axle is represented, i
β(k+1) electric current of subsequent time along β axle is represented, T
srepresent the sampling time, i
αk () represents the electric current of current time along α axle, i
βk () represents the electric current of current time along β axle, U
αk () represents the output voltage of current time along the three-phase VSC of α axle, U
βk () represents the output voltage of current time along the three-phase VSC of β axle, e
αk () represents the three-phase alternating current system voltage component of current time along α axle, e
βk () represents the three-phase alternating current system voltage component of current time along β axle;
The equation under dq two-phase rotating coordinate system can be obtained through rotation transformation to formula (9):
I
d(k+1) the three-phase grid-connected inverter output current d axle component that the k moment dopes is represented, i
q(k+1) represent the three-phase grid-connected inverter output current q axle component that the k moment dopes, θ represents the space angle of electrical network;
Active power and the reactive power in (k+1) moment is drawn by following prediction equation:
P (k+1) represents the active power in (k+1) moment, and Q (k+1) represents the reactive power in (k+1) moment, e
drepresent the three-phase alternating current system voltage component of d axle, e
qrepresent the three-phase alternating current system voltage component of q axle;
Calculate the first cost function g, calculating formula is:
g=|(Q
*-Q(k+1)|+|(P
*-P(k+1)| (12)
P
*for given active power reference value, Q
*for given reactive power reference qref;
Calculate respectively at 8 three-phase VSC output voltage vector U
i8 the first cost function g under state, obtain the minimum value g of the first cost function g
min; Obtain the minimum value g corresponding to the first cost function g
minclosest to expecting power, thus obtain corresponding U
α(k), U
βk the size of (), by U
α(k), U
βk () substitutes into the U in formula (6) respectively
α, U
β, obtain the switch function that this k moment expects, and the switch function of this expectation passed to the brachium pontis on off state of subsequent time, and then control inverter side current transformer; Wherein, described closest to expecting that power is the minimum performance number Q (k+1) of the value of instruction first cost function g and P (k+1), i=0,1,2 ... 7.
3. the control device being applicable to VSC-MTDC system according to claim 1, is characterized in that, also comprises rectification side controller;
Described rectification side controller, for realizing following function:
Set up rectification side three-phase VSC model:
Discretization formula (13) obtains predicted current function:
K represents current time, and k+1 represents subsequent time, i
α(k+1) electric current of subsequent time along α axle is represented, i
β(k+1) electric current of subsequent time along β axle is represented, T
srepresent the sampling time, i
αk () represents the electric current of current time along α axle, i
βk () represents the electric current of current time along β axle, U
αk () represents the output voltage of current time along the three-phase VSC of α axle, U
βk () represents the output voltage of current time along the three-phase VSC of β axle, e
αk () represents the three-phase alternating current system voltage component of current time along α axle, e
βk () represents the three-phase alternating current system voltage component of current time along β axle;
The equation under dq two-phase rotating coordinate system can be obtained through rotation transformation to formula (14):
Active power and the reactive power in (k+1) moment is drawn by following prediction equation:
P (k+1) represents the active power in (k+1) moment, and Q (k+1) represents the reactive power in (k+1) moment, e
drepresent the three-phase alternating current system voltage component of d axle, e
qrepresent the three-phase alternating current system voltage component of q axle; i
d(k+1) the three-phase grid-connected inverter output current d axle component that the k moment dopes is represented, i
q(k+1) the three-phase grid-connected inverter output current q axle component that the k moment dopes is represented; θ represents the space angle of electrical network;
Calculate the second cost function g, calculating formula is:
g=|(Q
*-Q(k+1)|+|(P
dc+P
i-P(k+1))| (15)
Q
*for given reactive power reference qref, P
dcrepresent the direct current power after regulating, P
irepresent the active power that i-th other rotary substation consumes except rectification side current transformer in multi-terminal direct current transmission system, i>0;
Calculate respectively at the active-power P that i other rotary substation consumes
ii under state the second cost function g, obtains the minimum value g of the second cost function g
min; Obtain the minimum value g corresponding to the second cost function g
minclosest to expecting power, thus obtain corresponding U
α(k), U
βk the size of (), by U
α(k), U
βk () substitutes into the U in formula (6) respectively
α, U
β, obtain the switch function that this k moment expects, and the switch function of this expectation passed to the brachium pontis on off state of subsequent time, and then control rectification side current transformer; Wherein, described closest to expecting that power is the minimum performance number Q (k+1) of the value of instruction second cost function g and P (k+1), i>0.
4. the control device being applicable to VSC-MTDC system according to claim 1, is characterized in that, also comprise tuning controller;
Described tuning controller, for controlling main convertor and there is the current transformer of power adjustments, specific as follows:
Work as P
1(k+1) allow in adjustable range at power, i.e. P
1min<P
1(k+1) <P
1max, then under making VSC2, VSC3 be operated in power mode; When VSC1 power adjustments is not enough or out of service, P
1(k+1) adjustable range [P is not allowed at power
1min, P
1max] in, then make the operational mode of VSC2 switch to direct voltage pattern to bear power shortage, under VSC3 is operated in power mode by power mode;
Wherein: VSC1 represents main convertor, VSC2 is the current transformer with power adjustments, and VSC3 is the current transformer without regulating power; P
1minrepresent that rotary substation is exerted oneself lower limit, i.e. minimum load, P
1(k+1) value of current predictive power is represented, P
1maxrepresent that rotary substation is exerted oneself the upper limit, i.e. maximum output.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510137472.6A CN104779635B (en) | 2015-03-26 | 2015-03-26 | Suitable for the control device of VSC MTDC systems |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510137472.6A CN104779635B (en) | 2015-03-26 | 2015-03-26 | Suitable for the control device of VSC MTDC systems |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104779635A true CN104779635A (en) | 2015-07-15 |
CN104779635B CN104779635B (en) | 2017-07-11 |
Family
ID=53620944
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510137472.6A Expired - Fee Related CN104779635B (en) | 2015-03-26 | 2015-03-26 | Suitable for the control device of VSC MTDC systems |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104779635B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106803683A (en) * | 2017-03-15 | 2017-06-06 | 郑州轻工业学院 | A kind of two-way AC/DC convertor model prediction current control method |
CN114070049A (en) * | 2021-10-26 | 2022-02-18 | 西安理工大学 | BOOST converter power switching control method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2068416A1 (en) * | 2007-12-05 | 2009-06-10 | Siemens Aktiengesellschaft | Advance control for a voltage source converter |
CN103050967A (en) * | 2013-01-06 | 2013-04-17 | 华北电力大学(保定) | Active disturbance control method of flexible direct current power transmission system |
CN103701148A (en) * | 2013-12-16 | 2014-04-02 | 南方电网科学研究院有限责任公司 | Starting control method for accessing VSC-MTDC system of large wind farm |
-
2015
- 2015-03-26 CN CN201510137472.6A patent/CN104779635B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2068416A1 (en) * | 2007-12-05 | 2009-06-10 | Siemens Aktiengesellschaft | Advance control for a voltage source converter |
CN103050967A (en) * | 2013-01-06 | 2013-04-17 | 华北电力大学(保定) | Active disturbance control method of flexible direct current power transmission system |
CN103701148A (en) * | 2013-12-16 | 2014-04-02 | 南方电网科学研究院有限责任公司 | Starting control method for accessing VSC-MTDC system of large wind farm |
Non-Patent Citations (5)
Title |
---|
SANGSHIN KWAK: "Model-Predictive Direct Power Control With Vector Preselection Technique for Highly Efficient Active Rectifiers", 《IEEE TRANSACTIONS ON INDUSTRIAL INFORMATICS》 * |
任敬国: "基于改进定有功功率控制特性的VSC-MTDC***仿真", 《电力***自动化》 * |
蔡长青等: "光伏并网逆变***无差拍预测直接功率控制方法", 《电力电子技术》 * |
赵方平等: "三相并网逆变器直接功率控制和直接功率预测控制的对比", 《电工技术学报》 * |
阮思烨等: "多端电压源型直流输电***的控制策略", 《电力***自动化》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106803683A (en) * | 2017-03-15 | 2017-06-06 | 郑州轻工业学院 | A kind of two-way AC/DC convertor model prediction current control method |
CN106803683B (en) * | 2017-03-15 | 2019-08-02 | 郑州轻工业学院 | A kind of two-way AC/DC convertor model prediction current control method |
CN114070049A (en) * | 2021-10-26 | 2022-02-18 | 西安理工大学 | BOOST converter power switching control method |
CN114070049B (en) * | 2021-10-26 | 2024-01-16 | 西安理工大学 | BOOST converter power switching control method |
Also Published As
Publication number | Publication date |
---|---|
CN104779635B (en) | 2017-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11791632B2 (en) | High-frequency uncontrolled rectifier-based DC transmission system for offshore wind farm | |
CN104659802B (en) | A kind of VSC-HVDC alternating voltage-frequency coordination control method improving AC system transient stability | |
CN105262113A (en) | Photovoltaic power generation system reactive power control method based on probabilistic fuzzy neural network | |
Taheri et al. | Application of Synchronous Static Series Compensator (SSSC) on enhancement of voltage stability and power oscillation damping | |
CN104810843A (en) | Active power control method for MW-level battery energy storage system based on battery charge state | |
Jurado | Power supply quality improvement with a SOFC plant by neural-network-based control | |
CN104779635A (en) | Controller device suitable for VSC-MTDC (voltage source converter-multi-terminal direct current) system | |
CN103972921A (en) | Power balance based low-voltage riding through and coordinative controlling method of permanent magnet synchronous wind power generator | |
Xi et al. | STATCOM operation strategy under power system faults | |
CN114583706A (en) | Direct-current power distribution system for treating low voltage at tail end | |
Tyagi et al. | Small hydro based grid forming converter having power sharing and synchronization capability with DFIG based WECS | |
Magueed et al. | Design of robust converter interface for wind power applications | |
CN204290321U (en) | Micro-capacitance sensor voltage perturbation control system | |
Djehaf et al. | Modeling of a multi-level converter based VSC HVDC supplying a dead load | |
Chao et al. | Approach on nonlinear control theory for designing STATCOM controller | |
Liu et al. | A novel nonlinear decoupled controller for VSC-HVDC system | |
Vignali | The impact of renewable energy on voltage stability and fault level in power systems | |
Aluko | Modelling and performance analysis of doubly fed induction generator wind farm | |
Theisen et al. | Stability analysis of an offshore grid supplied by a HVDC-VSC | |
CN104007349A (en) | Fuzzy control island detection method based on wavelet transformation | |
Wang et al. | Two-stage voltage control strategy for PV plants based on variable droop control | |
Karkevandi et al. | A novel adaptive hierarchical control structure for inverter based DGs under microgrid dynamics | |
Zhang et al. | Research of MMC-UPFC Control Strategy | |
Jing | Control and operation of MMC-HVDC system for connecting offshore wind farm | |
Xue et al. | Hierarchical Voltage/VAR security boundaries and margin assessments in Multi-HVDC infeed systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170711 |
|
CF01 | Termination of patent right due to non-payment of annual fee |