CN107294095A - Occur the system and control method with ZIP virtual loads with harmonic voltage - Google Patents

Occur the system and control method with ZIP virtual loads with harmonic voltage Download PDF

Info

Publication number
CN107294095A
CN107294095A CN201710468735.0A CN201710468735A CN107294095A CN 107294095 A CN107294095 A CN 107294095A CN 201710468735 A CN201710468735 A CN 201710468735A CN 107294095 A CN107294095 A CN 107294095A
Authority
CN
China
Prior art keywords
msub
mrow
msup
phase
mfrac
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
CN201710468735.0A
Other languages
Chinese (zh)
Other versions
CN107294095B (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.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
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 Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN201710468735.0A priority Critical patent/CN107294095B/en
Publication of CN107294095A publication Critical patent/CN107294095A/en
Application granted granted Critical
Publication of CN107294095B publication Critical patent/CN107294095B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/01Arrangements for reducing harmonics or ripples
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • 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/40Arrangements for reducing harmonics

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Inverter Devices (AREA)

Abstract

The invention discloses a kind of system and control method occurred with harmonic voltage with ZIP virtual loads, including power network, three-phase STS static state switching switches, module, three-phase series manifold type transformer group and ZIP virtual load modules occur for harmonic voltage.Wherein described ZIP virtual loads module is by acquisition system output voltage electric current and gives outside active power and reactive power, according to ZIP virtual load algorithms, realizes the foundation of ZIP virtual load models.A kind of system and control method occurred with harmonic voltage with ZIP virtual loads disclosed by the invention, both the local load in power system can have been simulated, the function of load on-site elimination can also be realized, it can be applied to draw water energy storage, photovoltaic from a variety of occasions such as cleer and peaceful electric automobile, simultaneity factor is also equipped with the function of harmonic wave generation, improves the completeness and reliability of system.

Description

Occur the system and control method with ZIP virtual loads with harmonic voltage
Technical field
The invention belongs to power quality analysis and control field, it is related to a kind of with harmonic voltage generation and ZIP virtual negatives The system and control method of lotus, can both simulate the local load in power system, and can also realize the work(of load on-site elimination Can, can be applied to draw water energy storage, photovoltaic are from a variety of occasions such as cleer and peaceful electric automobile, and simultaneity factor is also equipped with the work(of harmonic wave generation Can, improve the completeness and novelty of system.
Background technology
In power system, the characteristic such as voltage magnitude and frequency change is common existing in power system caused by load change As being that the control of this off-grid operation needs to set about from these fundamental problems, to propose reliable and effective solution.
At present, the research on load model in power system has become focus, such as Chinese invention patent document《It is a kind of The modeling method of power system load》(publication number 101789598A) and《The Dynamic Equivalence of middle small hydropower station cluster》It is (open Number 106651641A), wherein:
Chinese invention patent application file (101789598A) is in disclosed in 28 days July in 2010《A kind of power system is born The modeling method of lotus》, the present invention from load point classification, load model parameters polymerization, solve total body examination debate method identifiability and The problem of many solutions, disclose a kind of modeling method of power system load;
Chinese invention patent application file (106651641A) is in disclosed in 10 days Mays in 2017《Middle small hydropower station cluster Dynamic Equivalence》, present invention selection Equivalent Model is three rank parallel operation of generator static state ZIP load models, small power station in selection Active power, reactive power, voltage, the frequency of cluster power transmission interconnection Phasor Measurements data, to input, are adopted as input data Equivalent Model parameter is recognized with Self Adaptive Control parameter improvement differential evolution algorithm.
The existing granted patent on load model is made a general survey of, the problem of it is present is as follows:
1st, the foundation of load model can not simulate the local load in power system, cannot more realize load on-site elimination Function;
2nd, load model system sets up not comprehensive enough on the basis of fundamental wave mostly, does not account for the feelings comprising harmonic voltage Condition, and the presence of Harmonic in Power System voltage is extensive, it is necessary to be included within research range.
The content of the invention
In order to solve the problem of prior art is present, occur the invention discloses one kind with harmonic voltage and ZIP is virtual The system and control method of load, pass through acquisition system output voltage electric current and given outside active power and idle work( Rate, according to ZIP virtual load algorithms, realizes the foundation of ZIP virtual load models, the control method can be simulated in power system Local load, and realize the function of load on-site elimination.
The object of the present invention is achieved like this.
The invention provides a kind of system occurred with harmonic voltage with ZIP virtual loads, including power network, three-phase STS Module, three-phase series manifold type transformer group and ZIP virtual load modules occur for static state switching switch, harmonic wave;The harmonic wave hair Raw module includes DC source, three-phase tri-level PWM inverter, three-phase LC wave filters and three-phase STS1Static state switching switch;It is described ZIP virtual loads module includes electric capacity C, three-phase tri-level PWM converter and three-phase L-type wave filter;
The power network three-phase output end connects one to one with three-phase STS static state switching switch input terminals, and described three Phase STS static state switching output switching terminals connect one to one with the three-phase series manifold type transformer group secondary input;Institute State direct current source output terminal to connect one to one with the three-phase tri-level PWM inverter input, the three-phase tri-level PWM is inverse Become device output end and the inductance input of the three-phase LC wave filters connects one to one, the electric capacity of the three-phase LC wave filters is defeated Go out end and the three-phase STS1Static state switching switch input terminal connects one to one, the three-phase STS1Static state switching switch output End connects one to one with the three-phase series manifold type transformer group primary side input, the three-phase series manifold type transformer Group primary side output end and the electric capacity neutral point N short circuits of the three-phase LC wave filters;The three-phase series manifold type transformer group pair Side output end connects one to one with the three-phase L-type filter input end, the three-phase L-type filter output and described three Phase three-level pwm current transformer input connects one to one, and the three-phase tri-level PWM converter output end and electric capacity C are one by one Correspondence connection.
Present invention also offers a kind of control method occurred with harmonic voltage with the system of ZIP virtual loads, including Following steps:
Step 1, the three-phase series manifold type transformer group secondary output end three-phase phase in a switch periods is first gathered Voltage ua、ub、ucWith three-phase current ia、ib、ic, then first pass through coordinate transform and respectively obtain three-phase phase voltage α beta -axis components uα、uβ With three-phase current α beta -axis components iα、iβ, then lock phase angle γ obtained by phase-angle detection;
The coordinate transform formula is respectively:
The phase-angle detection formula is:
Step 2, according to the three-phase phase voltage α beta -axis components u obtained in step 1α、uβWith three-phase current α beta -axis components iα、iβ And lock phase angle γ, three-phase phase voltage dq axis components u is obtained by coordinate transformd、uqWith three-phase current dq axis components id、iq
The coordinate transform formula is respectively:
ud=uβcosγ+uαsinγ
uq=-uβsinγ+uαcosγ
id=iβcosγ+iαsinγ
iq=-iβsinγ+iαcosγ
Step 3, according to the three-phase phase voltage dq axis components u obtained in step 2d、uq, give active-power Po, give idle Power Qo, obtained three-phase phase voltage frequency departure amount △ f are detected, three-phase current dq axles are obtained by ZIP virtual load algorithms Instruct component idref、iqref
The ZIP virtual loads algorithmic formula is:
Wherein, kpf,kqfRespectively the frequency variation coefficient of d axles and q axles, p1,p2,p3The ZIP virtual negatives respectively given Lotus middle impedance, electric current, power component are in the ratio shared by q axles, q1,q2,q3ZIP virtual loads middle impedance, the electricity respectively given Stream, power component are in d axle proportions, if p1,p2,p3Summation, q1,q2,q3Summation is 1;
Step 4, according to the three-phase current dq axis components i obtained in step 2d、iqWith the three-phase current dq obtained in step 3 Axle instruction component idref、iqref, the drive signal dq axles point for obtaining three-phase tri-level PWM converter are calculated by controller equation Measure dd、dq, then obtain by coordinate transform the three-phase driving signal d of three-phase tri-level PWM convertera、dbAnd dc
Controller equation is:
Wherein, KpFor the proportionality coefficient of controller, KiFor the integral coefficient of controller, s is Laplace operator;
Coordinate transform formula is:
Relative to existing technology, beneficial effects of the present invention are:
1st, a kind of system and control method with harmonic voltage generation and ZIP virtual loads disclosed by the invention, both may be used To simulate the local load in power system, the function of load on-site elimination can also be realized, draw water energy storage, photovoltaic is can be applied to From a variety of occasions such as cleer and peaceful electric automobile;
2nd, system possesses the function of harmonic wave generation, improves the completeness and reliability of system.
Brief description of the drawings
Fig. 1 is that the system topology figure with ZIP virtual loads occurs with harmonic voltage;
Fig. 2 is that the System Control Figure with ZIP virtual loads occurs with harmonic voltage.
Each several part label and part corresponding relation are as follows in figure:
Module occurs for 20- harmonic voltages
21- DC sources
22- three-phase tri-level PWM inverters
30- three-phase series manifold type transformer groups
40-ZIP virtual load modules
41- three-phase tri-level PWM converters
Embodiment
Utilize Matlab/Simulink emulation platform building system models.
Reference picture 1, it is of the present invention to occur to include power network, three-phase STS with the system of ZIP virtual loads with harmonic voltage Module 20, three-phase series manifold type transformer group 30 and ZIP virtual loads module 40 occur for static state switching switch, harmonic wave;It is described Module 20, which occurs, for harmonic wave includes DC source 21, three-phase tri-level PWM inverter 22, three-phase LC wave filters, three-phase STS1Static state is cut Change switch;The ZIP virtual loads module 40 includes electric capacity C, three-phase tri-level PWM converter 41, three-phase L-type wave filter.
The power network three-phase output end connects one to one with three-phase STS static state switching switch input terminals, and described three Phase STS static state switching output switching terminals connect one to one with the secondary input of three-phase series manifold type transformer group 30; The direct current source output terminal connects one to one with the input of three-phase tri-level PWM inverter 22, the three-phase tri-level The output end of PWM inverter 22 and the inductance input of the three-phase LC wave filters connect one to one, the three-phase LC wave filters Electric capacity output end and the three-phase STS1Static state switching switch input terminal connects one to one, the three-phase STS1Static state switching Output switching terminal connects one to one with the primary side input of three-phase series manifold type transformer group 30, the three-phase series coupling The box-like primary side output end of transformer group 30 and the electric capacity neutral point N short circuits of the three-phase LC wave filters;The three-phase series coupling The secondary output end of formula transformer group 30 connects one to one with the three-phase L-type filter input end, the three-phase L-type wave filter Output end connects one to one with the input of three-phase tri-level PWM converter 41, the three-phase tri-level PWM converter 41 Output end connects one to one with electric capacity C.
In the present embodiment, power network is that line voltage is 380V three-phase alternating currents;The voltage of DC source 21 is 800V;Three-phase tri-level PWM inverter 22 power grade 20KW, rated voltage 380V, switching frequency 16KHz;Three-phase LC filter inductances L1= 0.12mH, electric capacity C1=30uF;The no-load voltage ratio of three-phase series manifold type transformer group 30 is 1;The power of three-phase tri-level PWM converter 41 Grade 20KW, rated voltage 380V, switching frequency 16KHz, three-phase L-type filter inductance L2=0.12mH.
Described occurs the system control method with ZIP virtual loads referring to Fig. 2, including following step with harmonic voltage Suddenly:
Step 1, the three-phase series manifold type transformer group secondary output end three-phase phase in a switch periods is first gathered Voltage ua、ub、ucWith three-phase current ia、ib、ic, then first pass through abc/ α β coordinate transforms and respectively obtain three-phase phase voltage α β axles point Measure uα、uβWith three-phase current α beta -axis components iα、iβ, then lock phase angle γ obtained by phase-angle detection.
Abc/ α β coordinate transform formula are respectively:
Phase-angle detection formula is:
Step 2, according to the three-phase phase voltage α beta -axis components u obtained in step 1α、uβWith three-phase current α beta -axis components iα、iβ And lock phase angle γ, three-phase phase voltage dq axis components u is obtained by α β/dq coordinate transformsd、uqWith three-phase current dq axis components id、 iq
α β/dq coordinate transform formula are respectively:
ud=uβcosγ+uαsinγ
uq=-uβsinγ+uαcosγ
id=iβcosγ+iαsinγ
iq=-iβsinγ+iαcosγ
Step 3, according to the three-phase phase voltage dq axis components u obtained in step 2d、uq, give active-power Po, give idle Power Qo, obtained three-phase phase voltage frequency departure amount △ f are detected, is calculated by ZIP virtual loads algorithm and obtains three-phase current Dq axles instruction component idref、iqref
The ZIP virtual loads algorithmic formula is:
Wherein, kpf,kqfRespectively the frequency variation coefficient of d axles and q axles, p1,p2,p3The ZIP virtual negatives respectively given Lotus middle impedance, electric current, power component are in the ratio shared by q axles, q1,q2,q3ZIP virtual loads middle impedance, the electricity respectively given Stream, power component are in d axle proportions, if p1,p2,p3Summation, q1,q2,q3Summation is 1.
Step 4, according to the three-phase current dq axis components i obtained in step 2d、iqWith the three-phase current dq obtained in step 3 Axle instruction component idref、iqref, the drive signal dq axles point for obtaining three-phase tri-level PWM converter are calculated by controller equation Measure dd、dq, then obtain by dq/abc coordinate transforms the three-phase driving signal d of three-phase tri-level PWM convertera、dbAnd dc
Controller equation is:
Wherein, KpFor the proportionality coefficient of controller, KiFor the integral coefficient of controller, s is Laplace operator.This implementation K is taken in examplep=1, Ki=10.
Dq/abc coordinate transform formula are:

Claims (2)

1. a kind of system occurred with harmonic voltage with ZIP virtual loads, it is characterised in that static including power network, three-phase STS Module (20), three-phase series manifold type transformer group (30) and ZIP virtual load modules occur for switching switch, harmonic wave;It is described humorous Module (20), which occurs, for ripple includes DC source (21), three-phase tri-level PWM inverter (22), three-phase LC wave filters and three-phase STS1It is quiet State switching switch;The ZIP virtual loads module (40) includes electric capacity C, three-phase tri-level PWM converter (41) and three-phase L-type Wave filter;
The power network three-phase output end connects one to one with three-phase STS static state switching switch input terminals, the three-phase STS Static state switching output switching terminal connects one to one with three-phase series manifold type transformer group (30) the secondary input;It is described DC source (21) output end connects one to one with three-phase tri-level PWM inverter (22) input, the electricity of three-phase three Flat PWM inverter (22) output end and the inductance input of the three-phase LC wave filters connect one to one, the three-phase LC filters The electric capacity output end of ripple device and the three-phase STS1Static state switching switch input terminal connects one to one, the three-phase STS1It is static Switching output switching terminal connects one to one with three-phase series manifold type transformer group (30) the primary side input, the three-phase Series coupled formula transformer group (30) primary side output end and the electric capacity neutral point N short circuits of the three-phase LC wave filters;The three-phase Series coupled formula transformer group (30) secondary output end connects one to one with the three-phase L-type filter input end, and described three Phase L-type filter output connects one to one with three-phase tri-level PWM converter (41) input, the electricity of three-phase three Flat PWM converter (41) output end connects one to one with electric capacity C.
2. a kind of control method as claimed in claim 1 occurred with harmonic voltage with the system of ZIP virtual loads, it is special Levy and be, comprise the following steps:
Step 1, three-phase series manifold type transformer group (30) the secondary output end three-phase phase in a switch periods is first gathered Voltage ua、ub、ucWith three-phase current ia、ib、ic, then first pass through coordinate transform and respectively obtain three-phase phase voltage α beta -axis components uα、uβ With three-phase current α beta -axis components iα、iβ, then lock phase angle γ obtained by phase-angle detection;
The coordinate transform formula is respectively:
<mrow> <msub> <mi>u</mi> <mi>&amp;alpha;</mi> </msub> <mo>=</mo> <mfrac> <msqrt> <mn>3</mn> </msqrt> <mn>3</mn> </mfrac> <msub> <mi>u</mi> <mi>b</mi> </msub> <mo>-</mo> <mfrac> <msqrt> <mn>3</mn> </msqrt> <mn>3</mn> </mfrac> <msub> <mi>u</mi> <mi>c</mi> </msub> </mrow>
<mrow> <msub> <mi>u</mi> <mi>&amp;beta;</mi> </msub> <mo>=</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <msub> <mi>u</mi> <mi>a</mi> </msub> <mo>-</mo> <mfrac> <mn>1</mn> <mn>3</mn> </mfrac> <msub> <mi>u</mi> <mi>b</mi> </msub> <mo>-</mo> <mfrac> <mn>1</mn> <mn>3</mn> </mfrac> <msub> <mi>u</mi> <mi>c</mi> </msub> </mrow>
<mrow> <msub> <mi>i</mi> <mi>&amp;alpha;</mi> </msub> <mo>=</mo> <mfrac> <msqrt> <mn>3</mn> </msqrt> <mn>3</mn> </mfrac> <msub> <mi>i</mi> <mi>b</mi> </msub> <mo>-</mo> <mfrac> <msqrt> <mn>3</mn> </msqrt> <mn>3</mn> </mfrac> <msub> <mi>i</mi> <mi>c</mi> </msub> </mrow>
<mrow> <msub> <mi>i</mi> <mi>&amp;beta;</mi> </msub> <mo>=</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <msub> <mi>i</mi> <mi>a</mi> </msub> <mo>-</mo> <mfrac> <mn>1</mn> <mn>3</mn> </mfrac> <msub> <mi>i</mi> <mi>b</mi> </msub> <mo>-</mo> <mfrac> <mn>1</mn> <mn>3</mn> </mfrac> <msub> <mi>i</mi> <mi>c</mi> </msub> </mrow>
The phase-angle detection formula is:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;gamma;</mi> <mo>=</mo> <mo>=</mo> <mfrac> <msub> <mi>u</mi> <mi>&amp;beta;</mi> </msub> <msqrt> <mrow> <msup> <msub> <mi>u</mi> <mi>&amp;alpha;</mi> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>u</mi> <mi>&amp;beta;</mi> </msub> <mn>2</mn> </msup> </mrow> </msqrt> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;gamma;</mi> <mo>=</mo> <mo>=</mo> <mfrac> <msub> <mi>u</mi> <mi>&amp;alpha;</mi> </msub> <msqrt> <mrow> <msup> <msub> <mi>u</mi> <mi>&amp;alpha;</mi> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>u</mi> <mi>&amp;beta;</mi> </msub> <mn>2</mn> </msup> </mrow> </msqrt> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced>
Step 2, according to the three-phase phase voltage α beta -axis components u obtained in step 1α、uβWith three-phase current α beta -axis components iα、iβAnd lock Phase angle γ, three-phase phase voltage dq axis components u is obtained by coordinate transformd、uqWith three-phase current dq axis components id、iq
The coordinate transform formula is respectively:
ud=uβcosγ+uαsinγ
uq=-uβsinγ+uαcosγ
id=iβcosγ+iαsinγ
iq=-iβsinγ+iαcosγ
Step 3, according to the three-phase phase voltage dq axis components u obtained in step 2d、uq, give active-power Po, give reactive power Qo, obtained three-phase phase voltage frequency departure amount △ f are detected, obtaining three-phase current dq axles by ZIP virtual load algorithms instructs Component idref、iqref
The ZIP virtual loads algorithmic formula is:
<mrow> <msub> <mi>i</mi> <mi>d</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>u</mi> <mi>q</mi> </msub> <msub> <mi>P</mi> <mi>o</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>(</mo> <mfrac> <mrow> <msup> <msub> <mi>u</mi> <mi>q</mi> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>u</mi> <mi>d</mi> </msub> <mn>2</mn> </msup> </mrow> <mrow> <msup> <msub> <mi>u</mi> <mi>o</mi> </msub> <mn>2</mn> </msup> </mrow> </mfrac> <mo>)</mo> <mo>+</mo> <msub> <mi>p</mi> <mn>2</mn> </msub> <msqrt> <mfrac> <mrow> <msup> <msub> <mi>u</mi> <mi>q</mi> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>u</mi> <mi>d</mi> </msub> <mn>2</mn> </msup> </mrow> <mrow> <msup> <msub> <mi>u</mi> <mi>o</mi> </msub> <mn>2</mn> </msup> </mrow> </mfrac> </msqrt> <mo>+</mo> <msub> <mi>p</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <msub> <mi>k</mi> <mrow> <mi>p</mi> <mi>f</mi> </mrow> </msub> <mi>&amp;Delta;</mi> <mi>f</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>u</mi> <mi>d</mi> </msub> <msub> <mi>Q</mi> <mi>o</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <mo>(</mo> <mfrac> <mrow> <msup> <msub> <mi>u</mi> <mi>q</mi> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>u</mi> <mi>d</mi> </msub> <mn>2</mn> </msup> </mrow> <mrow> <msup> <msub> <mi>u</mi> <mi>o</mi> </msub> <mn>2</mn> </msup> </mrow> </mfrac> <mo>)</mo> <mo>+</mo> <msub> <mi>q</mi> <mn>2</mn> </msub> <msqrt> <mfrac> <mrow> <msup> <msub> <mi>u</mi> <mi>q</mi> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>u</mi> <mi>d</mi> </msub> <mn>2</mn> </msup> </mrow> <mrow> <msup> <msub> <mi>u</mi> <mi>o</mi> </msub> <mn>2</mn> </msup> </mrow> </mfrac> </msqrt> <mo>+</mo> <msub> <mi>q</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <msub> <mi>k</mi> <mrow> <mi>q</mi> <mi>f</mi> </mrow> </msub> <mi>&amp;Delta;</mi> <mi>f</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msup> <msub> <mi>u</mi> <mi>q</mi> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>u</mi> <mi>d</mi> </msub> <mn>2</mn> </msup> </mrow> </mfrac> </mrow>
<mrow> <msub> <mi>d</mi> <mi>q</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>u</mi> <mi>d</mi> </msub> <msub> <mi>P</mi> <mi>o</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>p</mi> <mn>1</mn> </msub> <mo>(</mo> <mfrac> <mrow> <msup> <msub> <mi>u</mi> <mi>q</mi> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>u</mi> <mi>d</mi> </msub> <mn>2</mn> </msup> </mrow> <mrow> <msup> <msub> <mi>u</mi> <mi>o</mi> </msub> <mn>2</mn> </msup> </mrow> </mfrac> <mo>)</mo> <mo>+</mo> <msub> <mi>p</mi> <mn>2</mn> </msub> <msqrt> <mfrac> <mrow> <msup> <msub> <mi>u</mi> <mi>q</mi> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>u</mi> <mi>d</mi> </msub> <mn>2</mn> </msup> </mrow> <mrow> <msup> <msub> <mi>u</mi> <mi>o</mi> </msub> <mn>2</mn> </msup> </mrow> </mfrac> </msqrt> <mo>+</mo> <msub> <mi>p</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <msub> <mi>k</mi> <mrow> <mi>p</mi> <mi>f</mi> </mrow> </msub> <mi>&amp;Delta;</mi> <mi>f</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>u</mi> <mi>q</mi> </msub> <msub> <mi>Q</mi> <mi>o</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mn>1</mn> </msub> <mo>(</mo> <mfrac> <mrow> <msup> <msub> <mi>u</mi> <mi>q</mi> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>u</mi> <mi>d</mi> </msub> <mn>2</mn> </msup> </mrow> <mrow> <msup> <msub> <mi>u</mi> <mi>o</mi> </msub> <mn>2</mn> </msup> </mrow> </mfrac> <mo>)</mo> <mo>+</mo> <msub> <mi>q</mi> <mn>2</mn> </msub> <msqrt> <mfrac> <mrow> <msup> <msub> <mi>u</mi> <mi>q</mi> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>u</mi> <mi>d</mi> </msub> <mn>2</mn> </msup> </mrow> <mrow> <msup> <msub> <mi>u</mi> <mi>o</mi> </msub> <mn>2</mn> </msup> </mrow> </mfrac> </msqrt> <mo>+</mo> <msub> <mi>q</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <msub> <mi>k</mi> <mrow> <mi>q</mi> <mi>f</mi> </mrow> </msub> <mi>&amp;Delta;</mi> <mi>f</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msup> <msub> <mi>u</mi> <mi>q</mi> </msub> <mn>2</mn> </msup> <mo>+</mo> <msup> <msub> <mi>u</mi> <mi>d</mi> </msub> <mn>2</mn> </msup> </mrow> </mfrac> </mrow>
Wherein, kpf,kqfRespectively the frequency variation coefficient of d axles and q axles, p1,p2,p3Respectively hindered in given ZIP virtual loads Anti-, electric current, power component are in the ratio shared by q axles, q1,q2,q3ZIP virtual loads middle impedance, electric current, the work(respectively given Rate component is in d axle proportions, if p1,p2,p3Summation, q1,q2,q3Summation is 1;
Step 4, according to the three-phase current dq axis components i obtained in step 2d、iqRefer to the three-phase current dq axles obtained in step 3 Make component idref、iqref, the drive signal dq axles point for obtaining three-phase tri-level PWM converter (41) are calculated by controller equation Measure dd、dq, then obtain by coordinate transform the three-phase driving signal d of three-phase tri-level PWM converter (41)a、dbAnd dc
Controller equation is:
<mrow> <msub> <mi>d</mi> <mi>d</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>K</mi> <mi>p</mi> </msub> <mo>+</mo> <mfrac> <msub> <mi>K</mi> <mi>i</mi> </msub> <mi>s</mi> </mfrac> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msub> <mi>i</mi> <mrow> <mi>d</mi> <mi>r</mi> <mi>e</mi> <mi>f</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>i</mi> <mi>d</mi> </msub> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>d</mi> <mi>q</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>K</mi> <mi>p</mi> </msub> <mo>+</mo> <mfrac> <msub> <mi>K</mi> <mi>i</mi> </msub> <mi>s</mi> </mfrac> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msub> <mi>i</mi> <mrow> <mi>q</mi> <mi>r</mi> <mi>e</mi> <mi>f</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>i</mi> <mi>q</mi> </msub> <mo>)</mo> </mrow> </mrow>
Wherein, KpFor the proportionality coefficient of controller, KiFor the integral coefficient of controller, s is Laplace operator;
Coordinate transform formula is:
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>d</mi> <mi>a</mi> </msub> <mo>=</mo> <msub> <mi>d</mi> <mi>d</mi> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;gamma;</mi> <mo>+</mo> <msub> <mi>d</mi> <mi>q</mi> </msub> <mi>cos</mi> <mi>&amp;gamma;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>d</mi> <mi>b</mi> </msub> <mo>=</mo> <msub> <mi>d</mi> <mi>d</mi> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mo>(</mo> <mrow> <mi>&amp;gamma;</mi> <mo>-</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&amp;pi;</mi> </mrow> <mo>)</mo> <mo>+</mo> <msub> <mi>d</mi> <mi>q</mi> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mi>&amp;gamma;</mi> <mo>-</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&amp;pi;</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>d</mi> <mi>c</mi> </msub> <mo>=</mo> <msub> <mi>d</mi> <mi>d</mi> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mi>&amp;gamma;</mi> <mo>+</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&amp;pi;</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>d</mi> <mi>q</mi> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mi>&amp;gamma;</mi> <mo>+</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&amp;pi;</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>.</mo> </mrow> 2
CN201710468735.0A 2017-06-20 2017-06-20 System and control method with harmonic voltage generation and ZIP virtual load Active CN107294095B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710468735.0A CN107294095B (en) 2017-06-20 2017-06-20 System and control method with harmonic voltage generation and ZIP virtual load

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710468735.0A CN107294095B (en) 2017-06-20 2017-06-20 System and control method with harmonic voltage generation and ZIP virtual load

Publications (2)

Publication Number Publication Date
CN107294095A true CN107294095A (en) 2017-10-24
CN107294095B CN107294095B (en) 2019-09-10

Family

ID=60096800

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710468735.0A Active CN107294095B (en) 2017-06-20 2017-06-20 System and control method with harmonic voltage generation and ZIP virtual load

Country Status (1)

Country Link
CN (1) CN107294095B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109193749A (en) * 2018-07-23 2019-01-11 合肥工业大学 Harmonic wave generating device and control method based on recursive discrete Fourier transformation
CN117783757A (en) * 2024-02-23 2024-03-29 山东华天电气有限公司 Modularized simulation device and control method thereof
CN117783757B (en) * 2024-02-23 2024-05-14 山东华天电气有限公司 Modularized simulation device and control method thereof

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102185337A (en) * 2011-05-06 2011-09-14 深圳市禾望电气有限公司 Grid stimulating device
CN204228843U (en) * 2014-09-25 2015-03-25 上海市质量监督检验技术研究院 A kind of AC network simulator topological structure
CN106487240A (en) * 2015-08-28 2017-03-08 中国电力科学研究院 A kind of grid simulator with accurate harmonic voltage and virtual impedance control

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102185337A (en) * 2011-05-06 2011-09-14 深圳市禾望电气有限公司 Grid stimulating device
CN204228843U (en) * 2014-09-25 2015-03-25 上海市质量监督检验技术研究院 A kind of AC network simulator topological structure
CN106487240A (en) * 2015-08-28 2017-03-08 中国电力科学研究院 A kind of grid simulator with accurate harmonic voltage and virtual impedance control

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张炳达 等: "基于PWM技术的功率负荷模拟器", 《电力电子技术》 *
朱虹 等: "大功率电网模拟器系统研究", 《电源学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109193749A (en) * 2018-07-23 2019-01-11 合肥工业大学 Harmonic wave generating device and control method based on recursive discrete Fourier transformation
CN117783757A (en) * 2024-02-23 2024-03-29 山东华天电气有限公司 Modularized simulation device and control method thereof
CN117783757B (en) * 2024-02-23 2024-05-14 山东华天电气有限公司 Modularized simulation device and control method thereof

Also Published As

Publication number Publication date
CN107294095B (en) 2019-09-10

Similar Documents

Publication Publication Date Title
CN103107559B (en) A kind of method determining Distributed Power Flow controller system parameter
CN103236800B (en) Novel topological structure voltage source type inverter and adjusting method
CN103401459B (en) Triangularly connected chain H bridge suspended type inverter interphase DC side voltage balancing control method
CN106655805A (en) MMC-based multi-port hybrid power electronic transformer and control method thereof
CN106549417A (en) A kind of virtual synchronous generator control method and device of photovoltaic energy storage system
Song et al. Current harmonic minimization of a grid-connected photovoltaic 500kW three-phase inverter using PR control
CN106410839A (en) Active and reactive current coordination control-based photovoltaic grid-connected inverter control method
Changizian et al. Three-phase multistage system (DC-AC-DC-AC) for connecting solar cells to the grid
CN107332231A (en) Suitable for MW grades of high-power grid simulator topological structure
Sahoo et al. Hybrid generalised power theory for power quality enhancement
Naghizadeh et al. New single–phase three-wire interlinking converter and hybrid AC/LVDC microgrid
CN114884126A (en) Hybrid distribution transformer of integrated photovoltaic power generation system and control method
Sun et al. The application of modular multilevel matrix converter in fractional frequency offshore wind power system
CN107579666A (en) Multi-functional hybrid power electronic transformer and control method based on MMC matrix converters
CN107046297A (en) DC series permanent magnetism wind power plant topological structure and its control method
CN201797326U (en) Multifunctional power network simulator
Abdelkader et al. Integration of renewable energy sources in the dynamic voltage restorer for improving power quality using ANFIS controller
CN105514972A (en) PSCAD modeling and simulation method for grid-connected inverter during unsymmetrical failure of power grid
CN107294095A (en) Occur the system and control method with ZIP virtual loads with harmonic voltage
Dorantes et al. A filter design methodology of a VSC-HVDC system
Chen et al. Dual-AC-port modular multilevel converter for hybrid frequency integration of offshore wind power: Modeling and control
Parreiras et al. The true unit power factor converter with flexible arrangement for battery energy systems: Series and separate dc bus connection
Reddy et al. Hybrid renewable energy sources based four leg inverter for power quality improvement
Sule et al. Modeling and Performance assessment of an efficient photovoltaic based grid-connected inverter system with integration of LCL filter
Bina et al. High-accuracy power sharing in parallel inverters in an islanded microgrid using modified sliding mode control approach

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant