CN115632401B - SAPF parameter design method considering load and power grid impedance influence - Google Patents

SAPF parameter design method considering load and power grid impedance influence Download PDF

Info

Publication number
CN115632401B
CN115632401B CN202211645210.7A CN202211645210A CN115632401B CN 115632401 B CN115632401 B CN 115632401B CN 202211645210 A CN202211645210 A CN 202211645210A CN 115632401 B CN115632401 B CN 115632401B
Authority
CN
China
Prior art keywords
power grid
current loop
sapf
current
pcc
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202211645210.7A
Other languages
Chinese (zh)
Other versions
CN115632401A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202211645210.7A priority Critical patent/CN115632401B/en
Publication of CN115632401A publication Critical patent/CN115632401A/en
Application granted granted Critical
Publication of CN115632401B publication Critical patent/CN115632401B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • 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/24Arrangements for preventing or reducing oscillations of power in networks
    • 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/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02MAPPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
    • H02M1/00Details of apparatus for conversion
    • H02M1/12Arrangements for reducing harmonics from ac input or output
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02MAPPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
    • H02M7/00Conversion of ac power input into dc power output; Conversion of dc power input into ac power output
    • H02M7/42Conversion of dc power input into ac power output without possibility of reversal
    • H02M7/44Conversion of dc power input into ac power output without possibility of reversal by static converters
    • H02M7/48Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode
    • H02M7/53Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal
    • H02M7/537Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters
    • H02M7/539Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters with automatic control of output wave form or frequency
    • H02M7/5395Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters with automatic control of output wave form or frequency by pulse-width modulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/04Power grid 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]
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20The dispersed energy generation being of renewable origin
    • 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/20Active power filtering [APF]
    • 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)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Feedback Control In General (AREA)

Abstract

The invention discloses a SAPF parameter design method considering load and power grid impedance influences, and belongs to the technical field of alternating current-direct current conversion under a new energy power grid. Firstly, obtaining an optimal proportional coefficient expression of proportional link gains of a fundamental current loop and a power grid current loop controller under passive guidance according to a power grid current response equation of an SAPF double-loop current control system and correcting the optimal proportional coefficient expression; and selecting LCL filter parameters and calculating the gain of a proportional link of a power grid current loop/fundamental current loop controller according to a power grid current response equation, the corrected optimal proportional coefficient and the given minimum control bandwidth of a power grid current loop, and finally designing a proportional gain coefficient of a PCC power grid voltage feedforward term meeting the passivity requirement to ensure that the output admittance of the SAPF system is completely passive in the Nyquist range. The invention can ensure that the SAPF system is not influenced by the impedance fluctuation of a power grid and switching load, and has strong robustness.

Description

SAPF parameter design method considering load and power grid impedance influence
Technical Field
The invention belongs to the technical field of alternating current-direct current conversion under a new energy power grid, and particularly relates to an SAPF parameter design method considering load and power grid impedance influences.
Background
With the development of new power systems, more and more power electronics devices are used on the generation and demand sides, making the grid characteristics increasingly complex. The rectifier, especially the uncontrolled rectifier, brings serious power quality problems to the power grid, and the introduced harmonic current can cause power grid voltage distortion, generate extra line loss, accelerate the aging of power transmission and transformation equipment, and even cause the misoperation of a relay protection device. The parallel active power filter is used as a typical electric energy quality control device, has the characteristics of flexible compensation of harmonic current, quick response and the like, and can effectively relieve harmonic pollution in a power grid.
In order to suppress the inverter switching harmonics, an L-type or LCL-type filter is required to be provided at the output terminal thereof. Of these, LCL type filters are widely used due to their better harmonic attenuation capability. However, the inherent resonance characteristics of LCL type filters, coupled with the inevitable control delays in discrete control systems, pose a risk of instability in SAPF systems. In addition, the power grid environment of the SAPF grid-connected point can be changed due to the fluctuation of the power grid impedance and the switching of the load, and the system oscillation can be caused by the change of the power grid environment when the LCL filter and the controller are improperly designed or have poor robustness.
Disclosure of Invention
In view of the above, the present invention provides an SAPF parameter design method considering load and grid impedance influences, which takes passivity as a guide, obtains proportional link gain parameters of a fundamental current loop and a grid current loop controller and LCL filter parameters on the premise of simultaneously satisfying current loop control bandwidth and stability margin of an SAPF dual-loop current control system, and finally makes output admittance of an SAPF system completely passive by adjusting a proportional gain coefficient of a common coupling Point (PCC) grid voltage feedforward. The invention can ensure that the SAPF system still ensures the passivity of the LCL filter when the parameters of the LCL filter fluctuate, so that the three-phase LCL type SAPF system can stably operate when the impedance and the load of a power grid fluctuate.
The invention adopts the following technical scheme:
an SAPF parameter design method considering load and power grid impedance influences is provided, wherein SAPF parameters comprise proportional link gain of a fundamental current loop and a power grid current loop controller, LCL filter parameters and proportional gain coefficients of a PCC power grid voltage feedforward term;
the parameter design method comprises the following steps:
(1) Obtaining the optimal proportional coefficient expression of the proportional link gain of a fundamental current loop controller and the proportional link gain of a power grid current loop controller under passive guidance according to the real part information of the system output admittance when PCC voltage feedforward is not applied in the power grid current response equation of the SAPF double-loop current control system;
(2) Calculating a correction factor of the optimal proportionality coefficient expression in the step (1) according to the maximum negative deviation of the inductance value and the capacitance value of the inverter side in the LCL filter to obtain a corrected optimal proportionality coefficient expression;
(3) Obtaining an inequality equation set of LCL filter parameter values according to the corrected optimal proportionality coefficient expression, the stable condition of the power grid current response equation of the SAPF double-loop current control system and the given minimum control bandwidth of the power grid current loop, and selecting the LCL filter parameters in a feasible area surrounded by the inequality equation set;
(4) Substituting the LCL filter parameter value selected in the step (3) into the optimal proportional coefficient expression modified in the step (2) to obtain a modified optimal proportional coefficient, and obtaining the proportional link gain of the grid current loop controller and the proportional link gain of the fundamental current loop controller according to the modified optimal proportional coefficient, the selected minimum grid current loop control bandwidth and the grid side inductance value of the LCL filter;
(5) And (5) designing a proportional gain coefficient of a voltage feedforward term of the PCC power grid, which meets the passivity requirement, under the condition of the proportional link gain of the power grid current loop controller and the proportional link gain of the fundamental current loop controller obtained in the step (4).
Further, the SAPF dual-loop current control system comprises:
the bus voltage loop takes direct-current bus voltage and a direct-current bus voltage reference value as input, and obtains the amplitude of the fundamental current loop reference value after being adjusted by a bus voltage controller; obtaining a phase angle of a fundamental wave current loop reference value after the phase-locked loop is used for carrying out phase locking on the PCC network voltage; obtaining a fundamental current loop reference value according to the amplitude and the phase angle;
the fundamental current loop takes a fundamental current loop reference value and inverter side feedback current as input, and obtains output quantity after passing through the fundamental current loop controller;
the power grid current loop takes power grid current as input and obtains output quantity after passing through a power grid current loop controller;
the PCC power grid voltage feedforward item takes the PCC power grid voltage as input, and obtains output quantity after low-pass filtering and multiplying by a proportional gain coefficient;
and adding the output quantities of the fundamental wave current loop, the grid current loop and the PCC grid voltage feedforward item, modulating by SPWM to obtain a switching signal, and acting the switching signal on the three-phase voltage source inverter to obtain the output voltage of the inverter.
Further, the grid current response equation of the SAPF double-loop current control system is as follows:
Figure 711773DEST_PATH_IMAGE001
wherein the content of the first and second substances,i g (s)、u g (s)、i Lh (s) Harmonic current sources equivalent to the power grid current, the power grid voltage and the nonlinear load respectively;G L (s) Is a closed loop transfer function from the non-linear load current to the grid current;Y oa (s) For the total output admittance of the system,Y mL (s) The admittance generated for the nonlinear load input admittance coupling to the system,Y oc (s) For the system output admittance when applying the PCC voltage feed-forward,Y oca (s) To the output admittance when no PCC voltage feed-forward is applied,Y ocf (s) For the PCC voltage feed-forward effect equivalent admittance,Z g (s) Is the grid impedance.
Further, the expression of the corrected optimal proportionality coefficient in step (2) is as follows:
Figure 222389DEST_PATH_IMAGE002
wherein, the first and the second end of the pipe are connected with each other,K OPT in order to optimize the scaling factor before the correction,K MOPT in order to obtain the optimum scaling factor after the modification,
Figure 736547DEST_PATH_IMAGE003
in order to correct the factor(s),ω s in order to sample the angular frequency of the signal,ω r being the resonance angular frequency of the LCL filter,
Figure 226696DEST_PATH_IMAGE004
is the ratio of the inverter side inductance value and the grid side inductance value in the LCL filter.
Further, the correction factor is:
Figure 775489DEST_PATH_IMAGE005
wherein the content of the first and second substances,
Figure 66793DEST_PATH_IMAGE006
the inductance and the capacitance of the inverter side in the LCL filter are respectively the maximum negative deviation from the factory parameters.
Further, the inequality equation set of the LCL filter parameter values in step (3) includes the following constraint conditions:
the first constraint condition is:
Figure 192881DEST_PATH_IMAGE007
the second constraint condition is as follows:
Figure 860623DEST_PATH_IMAGE009
the third constraint condition is as follows:
Figure 263922DEST_PATH_IMAGE010
wherein the content of the first and second substances,ω cmin for a given minimum control bandwidth of the grid current loop,T s is the sampling time.
Further, in the step (4), a calculation formula of the proportional link gain of the grid current loop controller and the proportional link gain of the fundamental current loop controller is as follows:
wherein the content of the first and second substances,ω cmin for a given minimum control bandwidth of the grid current loop,K pf is the proportional link gain of the fundamental current loop controller,K ph the proportional link gain of the current loop controller of the power grid is obtained,L 2 is the grid side inductance value in the LCL filter.
Further, the step (5) is specifically as follows:
and (3) according to the value range of the proportional gain coefficient of the PCC power grid voltage feedforward term, the gain of the proportional link of the power grid current loop controller obtained in the step (5) and the gain of the proportional link of the fundamental current loop controller, obtaining the frequency response of the system output admittance when the PCC voltage feedforward is applied, and selecting the proportional gain coefficient of the PCC power grid voltage feedforward term as a design result in the range that the frequency response meets the passivity requirement.
Compared with the prior art, the invention has the following beneficial effects:
the invention provides a three-phase SAPF controller considering load fluctuation and power grid impedance influence and an LCL filter parameter design method. And then, proportional link gain parameters of fundamental waves and a power grid current loop controller are designed by utilizing an admittance model and taking passivity as guidance, and parameters of an LCL filter are obtained. And finally, adjusting a proportional gain coefficient of a PCC (point-to-point capacitor) power grid voltage feedforward link to ensure that the output admittance of the SAPF system is completely passive within the Nyquist range. The parameters designed by the invention can enable the SAPF system to simultaneously meet the requirements of control bandwidth and stability margin, so that the SAPF system not only has strong robustness on the impedance of a power grid on the basis of ensuring good harmonic wave treatment capability, but also can resist the influence brought by load fluctuation, and the application scenes of the SAPF are greatly expanded.
Drawings
FIG. 1 is a block diagram of a three-phase LCL type SAPF dual-loop current control system;
FIG. 2 is an equivalent transfer function control block diagram of a three-phase LCL model SAPF dual-loop current control system;
FIG. 3 is an equivalent circuit diagram of the SAPF dual-loop current control system, the grid, and the non-linear load;
FIG. 4 shows the current bandwidthω cmin Defined by three inequality constraints when chosen as 2000rad/s
Figure 850762DEST_PATH_IMAGE011
Andω r a feasible region map of (a);
FIG. 5 is a graph of current bandwidthω cmin 1000rad/s and 3000rad/s respectively, and is formed by three inequality constraint conditions
Figure 339512DEST_PATH_IMAGE012
Andω r a feasible region map of (a);
in FIG. 6, (a) to (h) are as
Figure 669999DEST_PATH_IMAGE013
And
Figure 927805DEST_PATH_IMAGE014
when wave motion occurs, adoptK OPT Parameter ofL 1 AndC f fluctuating admittance imaginary part
Figure 826491DEST_PATH_IMAGE015
A feasible domain variation graph of (1);
in FIG. 7, (a) to (d) are
Figure 160127DEST_PATH_IMAGE013
And
Figure 169671DEST_PATH_IMAGE014
when fluctuation occurs, adoptK MOPT Parameter ofL 1 AndC f fluctuating admittance imaginary part
Figure 406617DEST_PATH_IMAGE016
A feasible domain variation graph of (1);
fig. 8 shows the values of the LCL parameters of filter I and filter II;
fig. 9 (a) to (c) are graphs of SAPF grid current loop open-loop transfer functions Bode under three sets of parameters corresponding to the filter I and the filter II, respectively;
in fig. 10, (a) to (c) are three sets of parameters corresponding to the filter I and the filter II, respectively, whenL 1 AndC f SAPF output admittance when parameter fluctuation occursY oca (s) A frequency response characteristic map of (a);
FIG. 11 shows the case of applying the corresponding parameters of filter II
Figure 476205DEST_PATH_IMAGE017
Bode diagram of (1);
FIG. 12 (a) is a graph showing the absence of PCC voltage feed-forward term adjustmentY oca () AndY mL () Phasor diagram of interaction, fig. 12 (b) is after adjustment by the PCC voltage feed forward termY oc () And withY mL () A phasor diagram of the interactions;
FIG. 13 is a schematic view of
Figure 673968DEST_PATH_IMAGE018
Has a frequency response ofK ff A varying three-dimensional surface map;
FIG. 14 is a drawing showingY oc (s) Frequency response plot of (a);
fig. 15 (a) to (b) are respectively the power grid current waveform and Total Harmonic Distortion (THD) of the SAPF system using the parameters designed according to the present invention, when the power grid impedance is not applied;
FIG. 16 shows the system resonating after the grid impedance is applied without applying PCC voltage feed forward using the design parameters of the present invention; after the load is put into operation, the resonance is inhibited, and the waveform of the a-phase PCC voltage and the waveform of the power grid current are obtained;
fig. 17 (a) to (b) show the waveform and Total Harmonic Distortion (THD) of the current in the power grid, which are still stable after the system is stabilized after the impedance of the power grid is applied and the load is removed, respectively, in the SAPF system using the parameters designed according to the present invention.
Detailed Description
In order to more specifically describe the present invention, the following detailed description is provided for the technical solution of the present invention with reference to the accompanying drawings and the specific embodiments.
As shown in FIG. 1, in the SAPF dual-loop current control strategy, the grid currenti g Side current of inverteri inv Dc bus voltageu dc And PCC point voltageu pcc Is sampled as a control variable. Wherein for the network current loop, the goal is to convert the network currenti g Harmonic component ofi gh The control is zero. Considering that a Proportional Resonant (PR) controller employed by a grid current loop controller has a high gain at a specific subharmonic frequency, the grid currenti g Harmonic waves do not need to be extracted, and the harmonic waves can be directly used as the input of a power grid current loop controller. For the fundamental current loop, the control objective is to keep the bus voltage stable. For this purpose, the inverter side currenti inv Tracking a reference current by a Proportional Resonant (PR) controller of a fundamental current loopi f, ref . Reference current of fundamental current loopi f, ref By means of a Phase Locked Loop (PLL) with a PCC voltageu pcc The phase of (2) remains synchronized. In this embodiment, in order not to introduce the influence of the PLL, the bandwidth thereof may be set to be low, for example, 25Hz. PCC grid voltageu pcc Passing through a low pass filterG ff (s) Post-multiplying by a proportional gain factorK ff As a feedforward term, the feedforward term is superposed with the output of a power grid current loop controller and a fundamental current loop controller, a switching signal is obtained after SPWM modulation, the switching signal acts on a three-phase Voltage Source Inverter (VSI), and finally the output voltage of the inverter is obtainedu invu inv Obtaining output current after filtering by LCLi outi out Harmonic component in and non-linear load currenti L And the harmonic content in the power grid current is reduced by offsetting.
The principles and process flow of the present invention are described below.
(1) Using the rest shown in figure 2
Figure 611837DEST_PATH_IMAGE019
The transfer function control block diagram under the coordinate system establishes an admittance model of the SAPF system, and the detailed discussion is as follows:
the dual loop current control block diagram of FIG. 1 may be equivalent to the quiescent state of FIG. 2
Figure 578656DEST_PATH_IMAGE019
And controlling the block diagram of the transfer function under the coordinate system. As shown in fig. 2, which is a sampled-data system, there are continuous and discrete time signals, and the sampled discrete variables are represented by the corresponding continuous variables with marks. By means of sampling switches between continuous and discrete signals (with sampling period ofT s ) And zero order keeperG h (s) And (6) performing conversion. For a continuous signaliCorresponding to a sampled discrete signal ofi * The relationship between the two can be expressed in the frequency domain as:
Figure 84723DEST_PATH_IMAGE020
wherein the content of the first and second substances,
Figure 628837DEST_PATH_IMAGE021
in order to sample the angular frequency of the signal,sin order to be the laplacian operator,jis a symbol of an imaginary number,kthe values are integers. The above expression shows that the sampled signal is introduced into the image andω s is a period.
The zero order keeper (ZOH) in fig. 2 is used to equate to a pulse width modulation Process (PWM), and its corresponding transfer function can be expressed as follows:
Figure 245763DEST_PATH_IMAGE022
through ZOH, discrete voltage signalsu * inv Can be converted into actual voltage signalsu inv . In addition, the introduction of the discrete control is inevitableThe calculation delay, which is dependent on the configuration of the control system, is typically one sampling period, i.e.D(z) = z -1z = exp(sT s ). In the present embodiment, for the fundamental current reference valuei f, ref And grid current harmonic component reference valuei gh, ref Controller for realizing zero-error tracking and power grid current loopG ch (s) And fundamental current loop controllerG cf (s) All adopt PR controller, the resonance (R) unit of PR controller issThe domain can be expressed as:
Figure 693187DEST_PATH_IMAGE023
wherein the content of the first and second substances,K rn is as followsnThe resonant gain of the sub-harmonics,
Figure 370156DEST_PATH_IMAGE024
and are respectively the secondnThe corresponding compensation angle of the sub-R unit and the resonance angular frequency. The proportional gains of the fundamental current loop and the grid current loop controller are respectively controlled byK pf AndK ph and (4) showing.G ch (z) AndG cf (z) Are respectively composed ofG ch (s) AndG cf (s) Dispersing to obtain the product.G ff (s) Is a low-pass filter which is discrete toG ff (z)。K ff The value of the proportional gain coefficient is 0~1. Due to the DC bus voltageu dc Is controlled to be a constant value, and the corresponding control link is omitted in fig. 2.
For an LCL filter, its frequency characteristics can be expressed as the following admittance:
Figure 276933DEST_PATH_IMAGE025
a mixed output admittance model can be obtained by utilizing the inteswample modeling method, and the method can be used for solving the problem that the traditional method is not suitable for the traditional methodz = exp(sT s ) And realizing unification. From fig. 2, the relationship between all variables in the sapf control system can be expressed as:
Figure 556604DEST_PATH_IMAGE026
according to the above equation and ignoring the effect of the image introduced during sampling, byu pcc (s) The closed loop response under the SAPF dual loop current control framework as input can be derived as follows:
Figure 498015DEST_PATH_IMAGE027
wherein the content of the first and second substances,G i (s),G L (s) Respectively corresponding input-to-output closed loop transfer functions,Y oc (s) To output admittance, it can be expressed as:
Figure 345886DEST_PATH_IMAGE028
wherein the content of the first and second substances,Y c1 (z),Y c2 (z),Y o1 (z),Y o2 (z) Can be respectively obtained by the dispersion of the zero-order keeper by the corresponding continuous domain admittance.
Output admittanceY oc (s) Phasors can be decomposed into output admittance without applying PCC voltage feed-forwardY oca (s) Admittance equivalent to PCC voltage feedforward actionY ocf (s) Namely:
Figure 864592DEST_PATH_IMAGE029
wherein, the first and the second end of the pipe are connected with each other,
Figure 823320DEST_PATH_IMAGE030
the total output admittance of the SAPF control system can be affected by adjusting the PCC voltage feed-forward term, whereinY ocf (s) For the PCC voltage feed-forward effect equivalent admittance,Y oca (s) For system output admittance when PCC voltage feed-forward is not applied,Y oc (s) The system output admittance is fed forward for the application of the PCC voltage.
When the PCC voltage exists at a frequency off x In harmonic disturbance, neglecting the coupling effect between different subharmonics, the relationship between the harmonic voltage and the current on the ac side of the nonlinear load connected to the PCC point can be expressed as:
Figure 478293DEST_PATH_IMAGE031
wherein the content of the first and second substances,Y L the amplitude-frequency characteristic of the input admittance of the nonlinear load can be obtained by sweeping frequency to a load port, and the input admittance of the load presents inductance for a three-phase uncontrolled rectifying load with an incoming reactor, which is common in the industry. Thus, the load currenti L (s) Can be regarded as harmonic current sourcei Lh (s) Superposition of input admittance model with nonlinear loading:
Figure 497064DEST_PATH_IMAGE032
substituting the above formula into a power grid current response equation of the SAPF closed-loop system to obtain:
Figure 378433DEST_PATH_IMAGE033
wherein the content of the first and second substances,Y mL (s) = G L (s)Y L (s) The admittance generated by the nonlinear load input admittance coupling with the system;G L (s)=1- G i (s),G L (s) For closed-loop transfer function from non-linear load current to grid currentThe number of the first and second groups is,G i (s) Is a closed loop transfer function of a power grid current loop. Input admittance of non-linear load can be foundY L (s) The total output admittance of the PCC port is influenced by the coupling of the SAPF control loop, and the stability of a grid-connected system is further influenced.
The above equation can be expressed as a davinin equivalent circuit, as shown in fig. 3. Will be provided withu pcc (s)= u g (s)- i g (s)Z g (s) After the formula is substituted, an ideal power grid can be deducedu g (s) The closed loop system response equation for the input is as follows:
Figure 764021DEST_PATH_IMAGE034
wherein the content of the first and second substances,Y oa (s) = Y mL (s) + Y oca (s) + Y ocf (s) The total output admittance for the SAPF to the PCC port.
(2) Utilizing SAPF system output admittance when PCC voltage feedforward is not applied as obtained in (1)Y oca (s) Design of proportional gain of fundamental current loop controllerK pf Proportional gain to grid current loop controllerK ph Meanwhile, value constraint conditions of LCL filter element parameters are given, and the specific solving method is as follows:
the stability condition of the SAPF system can be known from the closed-loop system response equation of (1): 1)G L (s) Gradual stabilization; 2) Loop gainY oa (s)Z g (s) The Nyquist stability criterion is satisfied. The condition 1) can be realized by designing parameters of a fundamental current loop and a grid current loop controller, and the sufficient condition for the condition 2) to be satisfied is total output admittanceY oa (s) From 0 toω N At a passive origin, i.e. Re, in the range (Nyquist angular frequency) Y oa (s) Is > 0 or
Figure 414445DEST_PATH_IMAGE035
Since the resonant unit of the current controller only has great influence on the amplitude-frequency characteristic near the resonant frequency, the solution is carried outY oca (s) In the real part, neglecting it, i.e.G cf (z)≈K pfG ch (z)≈K ph . Substitution intos = z = exp(jωT s ) The following can be obtained:
Figure 338539DEST_PATH_IMAGE036
wherein the content of the first and second substances,
Figure 831837DEST_PATH_IMAGE037
wherein the content of the first and second substances,ωis the angular frequency. When the temperature is higher than the set temperatureωIs less thanω N Time, always positive, analysisY oca (s) The real part of (A) is not difficult to find, in the Nyquist frequency range, as long asK pf AndK ph satisfies the following conditions:
Figure 132369DEST_PATH_IMAGE038
namely can makeY oca (s) The real part is greater than zero (only inω = ω s Zero at/6). Wherein the content of the first and second substances,
Figure 27512DEST_PATH_IMAGE039
a ratio of an inverter-side inductance value and a grid-side inductance value in the LCL filter;ω r the LCL filter natural resonant angular frequency.
The parameters of the current loop controller are designed under the condition of satisfying the above equation relationG L (s) And (4) gradually stabilizing. As shown in fig. 2, ignoreu pcc Andi L the power grid current loop gain can be obtained as follows:
Figure 122507DEST_PATH_IMAGE040
in general, SAPF grid current loop control bandwidthω c Much less than the resonant frequency of the LCL filterω r . Thus at frequencyω c Nearby, the LCL filter can be approximated as a first order L filter. Likewise, delayD(z) Can be approximately replaced by a first order Taylor expansion, i.e.
Figure 244047DEST_PATH_IMAGE041
. Substitution intoK pf = K OPT K ph Loop gainT i (s) Can be simplified as follows:
Figure 207324DEST_PATH_IMAGE042
wherein the content of the first and second substances,L = L 1 +L 2 the inductance value represents the sum of the inverter-side inductance value and the grid-side inductance value in the LCL filter.
According to the above formula, inω c ToT i (s) The amplitude of (d) is:
Figure 566761DEST_PATH_IMAGE043
according to the formula, the proportional gain of the power grid current loop controllerK ph Determining the control bandwidth of the current loop of the power networkω c When the network current loop controls the bandwidthω c After the setting, the user can set the setting,K ph can be expressed as:
Figure 98236DEST_PATH_IMAGE044
loop gainT i (s) In thatω c The phase of (a) is:
Figure 802013DEST_PATH_IMAGE045
wherein the content of the first and second substances,
Figure 709926DEST_PATH_IMAGE046
is the phase margin PM of the system. Will be provided withK ph After the substitution into the formula (I) is carried out,
Figure 455028DEST_PATH_IMAGE047
can be rewritten as:
Figure 16459DEST_PATH_IMAGE048
as can be seen from the above formula, the PM of the system is the bandwidth set by the systemω c AndK OPT and (6) determining. ByK OPT The expression shows that after LCL filter parameters are designedK OPT Is uniquely determined. Bandwidth at that timeω c The size of (b) determines the stability margin of the SAPF control system. Therefore, it is required toω c And
Figure 847012DEST_PATH_IMAGE049
a compromise is made on the choice of (a) to obtain a larger stability margin with a smaller system bandwidth. Fundamentally, the passivity design goal reduces one degree of freedom in the SAPF dual current loop control, and under the guidance of the design, in order to simultaneously meet the requirements of high bandwidth and larger PM, the LCL filter parameters can be reasonably designed and adjustedK OPT Is achieved by the value of (c). For the
Figure 152091DEST_PATH_IMAGE049
In other words, since the first-order taylor expansion is used to replace the delay element with a certain conservative property in the phase lag, it can be set to 45 °, and the first constraint on stability can be derived from the above equation:
Figure 751700DEST_PATH_IMAGE050
wherein the content of the first and second substances,ω cmin for minimum control bandwidth of the grid current loop, because of the grid current loop controllerG ch (z) Including multiple R units, to avoid an additional 0-dB ride-through,ω cmin should be greater than the harmonic frequency of the highest order compensation.
For SAPF, the resonant frequency of the LCL filter is typically higher due to the higher harmonic frequencies to be compensatedω r Is generally set higher thanω s And 6, meeting the requirement of larger current loop bandwidth of the power grid. At this time, for the dual current loop control structure without delay compensation, a pair of open loop instability poles exist. In order for the system to ensure stability, the phase-frequency curve should be only in accordance with the Nyquist stability criterionω r There is a-180 positive crossing. For this purpose, the loop gainT i (s) In thatω s The amplitude at/6 is less than 0dB atω r The amplitude at (a) should be greater than 0dB. Considering disturbance influence such as parameter fluctuation, the amplitude margin (GM) is generally set to 3dB to improve the system robustness, so the corresponding stability condition can be expressed as
Figure 359399DEST_PATH_IMAGE051
And
Figure 801881DEST_PATH_IMAGE052
the second and third restrictions on stability, respectively, are as follows:
Figure 520439DEST_PATH_IMAGE054
can find that whenω cmin After determination, the variables of the first, second and third stability constraints are allω r And
Figure 240133DEST_PATH_IMAGE055
. When LCL filter parametersω r And
Figure 641902DEST_PATH_IMAGE055
when the three conditions are met, the SAPF system can meet the requirements of bandwidth and stability margin on the premise of ensuring that the output admittance is passive.
The arithmetic example of the present invention shows that the highest harmonic compensation time is 19 th Since the power frequency voltage is 50Hz, correspondingly,ω cmin the value can be 2000rad/s. The boundaries of the three above-mentioned stability constraints are shown in fig. 4 when the sampling frequency is 15kHz, provided that the selection is made in a common region enclosed in fig. 4ω r And
Figure 181468DEST_PATH_IMAGE056
therefore, the SAPF control system designed by taking passivity as a guide has sufficient stability margin and strong robustness to power grid impedance fluctuation. It can be found that there are countless selection schemes for the LCL filter parameters in the feasible domain defined by the three constraint conditions, and in practice, suitable LCL filter parameters should be designed in the feasible domain according to the power level, reactive capacity constraint, output current ripple limitation, voltage sag limitation, and the like of the SAPF system.
When in useω cmin When 1000rad/s is taken, a feasible region surrounded by three constraint condition boundaries is shown in (a) in fig. 5, and after the bandwidth of the power grid current loop is increased to 3000rad/s, as shown in (b) in fig. 5, the area of the feasible region is continuously reduced, that is, the selectable range of the LCL filter parameters is smaller and smaller. This means that a very high bandwidth, while improving the response of the system, sacrifices the flexibility of LCL parameter selection, so that the bandwidth is selected to be equal to or greater thanω cmin That is, when the bandwidth is selected asω cmin Proportional gain of time-of-flight, grid current loop controllerK ph The corresponding representation is:
Figure 93929DEST_PATH_IMAGE057
(3) Analysis ofL 1 AndC f after the parameters are changed, the SAPF outputs the change rule of admittance frequency response,and pairK OPT Make a correction to obtainK MOPT The specific discussion process is as follows:
in practice, the amount of the liquid to be used,L 1 andC f the parameters of (a) change after the operation of the SAPF,K OPT will be changed along with the change of the temperature,K pf = K OPT K ph this condition will no longer be satisfied. At this time, admittanceY oca (s) Will be in the Nyquist frequency rangeω s Non-passive regions occur near/6. For the case of the passive region, reference may be made to the document "PassionEnhance for LCL-Filtered observer With Grid Current Control and Capacitor Current Active pointing, xuehua Wang", published in the journal IEEE Transactions on Power Electronics, 2022. In order to compensate the non-passive region, it is necessary to analyze the characteristics of the non-passive region to try to make the region exhibit a consistent inductance or capacitance.
In the invention, all parameters after fluctuation are represented by adding a left-falling part to the upper right corner of the parameter before corresponding fluctuation, for example, by representingL 1 After the parameters fluctuate. Assuming that the parameters after fluctuation satisfy:
Figure 402551DEST_PATH_IMAGE058
at this time, admittance
Figure 617631DEST_PATH_IMAGE059
The real part of (d) can be expressed as:
Figure 769127DEST_PATH_IMAGE060
Figure 95066DEST_PATH_IMAGE061
after the parameter fluctuates
Figure 523773DEST_PATH_IMAGE062
If the value is still greater than zero, the above formula shows
Figure 768810DEST_PATH_IMAGE063
In thatω s 6 and
Figure 548547DEST_PATH_IMAGE064
with spacing less than zero, i.e. admittance
Figure 38697DEST_PATH_IMAGE065
Non-passive regions are present.
Further, the admittance in this region was investigated
Figure 321910DEST_PATH_IMAGE066
Positive and negative of imaginary part of will
Figure 878794DEST_PATH_IMAGE067
Substitution into
Figure 4881DEST_PATH_IMAGE068
Taking its imaginary part:
Figure 672623DEST_PATH_IMAGE069
from the above formula, becauseK OPT > 0,
Figure 810343DEST_PATH_IMAGE070
Always on, therefore, it is true that
Figure 928341DEST_PATH_IMAGE071
Positive or negative of imaginary part depends on
Figure 151512DEST_PATH_IMAGE072
And
Figure 622945DEST_PATH_IMAGE073
relationships between them, as can be seen in FIG. 4ω r Far greater thanω s 6, therefore
Figure 739805DEST_PATH_IMAGE074
Should still be greater than
Figure 904070DEST_PATH_IMAGE075
At this time
Figure 972127DEST_PATH_IMAGE076
Is positive, and the admittance exhibits a capacitive property at that frequency.
However,
Figure 981671DEST_PATH_IMAGE077
the imaginary part of (A) cannot be directly deduced intuitively, and the positive and negative parts of the imaginary part deviate along with the fluctuation of the parameter
Figure 359563DEST_PATH_IMAGE078
And
Figure 553784DEST_PATH_IMAGE079
but is changed. To compensate for the passive region of the output admittance,
Figure 751547DEST_PATH_IMAGE080
should be at
Figure 564782DEST_PATH_IMAGE078
And
Figure 656235DEST_PATH_IMAGE079
remains positive throughout the variation, thus admittance
Figure 896724DEST_PATH_IMAGE081
In thatω s 6 and
Figure 706417DEST_PATH_IMAGE082
the frequency band will exhibit consistent capacitance, requiring only a lagging phase compensation for this region.
Solving for
Figure 323343DEST_PATH_IMAGE083
And an imaginary part ofThe sign is positive, which results in:
Figure 144668DEST_PATH_IMAGE084
wherein the content of the first and second substances,
Figure 447736DEST_PATH_IMAGE085
deviation when parameter fluctuates
Figure 354512DEST_PATH_IMAGE086
And
Figure 509550DEST_PATH_IMAGE087
after determination, the above formula is aboutω r And
Figure 575595DEST_PATH_IMAGE088
can satisfy the inequality constraint condition to ensure that the current situation is satisfied
Figure 157886DEST_PATH_IMAGE089
And
Figure 676592DEST_PATH_IMAGE087
in the following, the first step is to put the paper into the bag,
Figure 635321DEST_PATH_IMAGE090
the sign of (1) is positive. To obtain following
Figure 696818DEST_PATH_IMAGE078
And
Figure 574644DEST_PATH_IMAGE087
the invention adopts perturbation method to draw the variable trend of the feasible region
Figure 721591DEST_PATH_IMAGE089
And
Figure 952852DEST_PATH_IMAGE087
changing from-15% to 0, and increasing to 15% again
Figure 515462DEST_PATH_IMAGE091
Aboutω r And
Figure 439556DEST_PATH_IMAGE088
is shown in fig. 6. Fig. 6 is a diagram of the internal stability feasible region of fig. 4, and if there is an intersection between two feasible regions in the fluctuation range of the parameters, it indicates that there is a set of LCL filter parameters that can make the SAPF system in the middle of the time periodL 1 AndC f the stability constraint is still satisfied when the change occurs, andω s 6 and
Figure 73800DEST_PATH_IMAGE092
capacitive non-passive regions occur within the frequency band.
As shown in fig. 6, when the parameter isL 1 AndC f when the value of (A) is reduced, with
Figure 233386DEST_PATH_IMAGE078
And
Figure 738316DEST_PATH_IMAGE087
the intersection area is continuously reduced until disappearing when the concentration is changed from-1% to-15%, which shows that when the concentration is changed to be zeroL 1 AndC f when large negative fluctuation occurs, the non-passive area of the output admittance cannot present consistent characteristics, and the characteristics of the non-passive area change suddenly at a certain frequency, so that the inductance is converted into the capacitance, and the design of a phase compensation link is not facilitated. When inL 1 AndC f as the value of the parameter increases, with
Figure 223524DEST_PATH_IMAGE089
And
Figure 345064DEST_PATH_IMAGE087
increasing from 1% to 15%, the two regions always coincide, indicating an output admittance factor of the SAPFL 1 AndC f the non-passive area generated by the forward fluctuation is always capacitive.
From the above rule, the invention is aimed atK OPT Making a correction by a factor ofK MOPT
Figure 918128DEST_PATH_IMAGE093
Wherein the content of the first and second substances,
Figure 667778DEST_PATH_IMAGE094
Figure 199253DEST_PATH_IMAGE095
and
Figure 903030DEST_PATH_IMAGE096
is composed ofL 1 AndC f the maximum negative deviation of the parameter values is, in this embodiment,
Figure 810943DEST_PATH_IMAGE095
=
Figure 556045DEST_PATH_IMAGE097
= -0.15。
after applying the modified coefficients, the output admittance isω s 6 and
Figure 117476DEST_PATH_IMAGE098
non-passive regions will occur between the frequency bands. As in the analysis described above, it is preferred that,
Figure 948029DEST_PATH_IMAGE099
is always positive, so that
Figure 862895DEST_PATH_IMAGE100
The feasible region of positive sign of (1) is also by perturbation
Figure 587138DEST_PATH_IMAGE078
And
Figure 460416DEST_PATH_IMAGE087
is plotted in fig. 7. As can be seen from FIG. 7, the corrected coefficients are usedK MOPT After, no matter the parametersL 1 AndC f how to fluctuate within a deviation of 15% such that
Figure 778265DEST_PATH_IMAGE101
For the positive parameter robustness feasible region always including the stability feasible region, it should be noted that when the modified coefficient is adoptedK MOPT Then, the internal stability feasible region in fig. 4 needs to be corrected at the same time, and the internal stability feasible region in fig. 7 is the corrected result.
To verify the correctness of the above passivity-oriented parameter design method, the present embodiment takes a 30kVA SAPF prototype as an example, and the parameters thereof are listed in table 1 below:
TABLE 1 SAPF operating parameters
Figure 621456DEST_PATH_IMAGE102
The two sets of LCL filter parameters are listed in table 2 below:
TABLE 2 LCL Filter parameters
Figure 341150DEST_PATH_IMAGE103
Both sets of parameters were selected within the feasible region obtained from the above analysis, as shown in FIG. 8, where filter I corresponds toω r Is 24407rad/s,
Figure 742919DEST_PATH_IMAGE088
2.54, located in the corrected feasible region, whichK MOPT Is 0.0589 whenω cmin When the amount is set to 2000rad/s, the calculation is performedK ph Is a high-frequency signal of 2.89 omega,K pf is 0.17 omega. In order to highlight the effect of the correction method provided by the invention, the filterII corresponds toω r Is 38730rad/s,
Figure 548064DEST_PATH_IMAGE104
8, located in the overlapping area of the two feasible fields: when not correctedK OPT Is 0.4804 whenω cmin When the value is set to 2000rad/s, the calculation can be madeK ph Is a high-frequency wave of 3.22 omega,K pf is 1.55 omega; when the correction coefficient is adoptedK MOPT Is 0.0696, same asω cmin Then, calculate to obtainK ph Is a gas flow rate of 2.83 omega,K pf and is 0.20 omega.
Firstly, the stability of the SAPF closed-loop system is verified, and when the three groups of parameters are adopted, the SAPF power grid current loop open-loop transfer functions Bode diagrams are respectively shown as (a), (b) and (c) in FIG. 9. In accordance with the nyquist stability criterion,N + - N - should be equal to half the number of open loop instability poles. When in useω r > ω s At/6, there are two open-loop instability poles in the system, shown in (a), (b), (c) of FIG. 9N + - N - All are 1, which indicates that the SAPF system closed loop stability can be ensured by adopting the parameter design method. Further, from Bode fig. 9, the system has PM, GM, and bandwidth that are all large enough. Note that the PM of (b) in fig. 9 is slightly less than 45 °, which is caused by an error introduced by the linearization process of the delay element.
Then analyzing SAPF output admittance under the condition of adopting passive-oriented parameter design methodY oca (s) The frequency response characteristic of (2) is, as shown in (a), (b) and (c) of fig. 10, as the parameter fluctuates, it is not difficult to find from (a) and (c) of fig. 10, and the corrected passivity coefficient is usedK MOPT When the parameter fluctuates between-15% and +15%, the output admittance isω s The non-passive area around/6 presents consistent capacitance; and adopts the original passivity coefficientK OPT When is coming into contact withL 1 AndC f after the parameters have larger negative fluctuation at the same time, the output admittance isω s The non-passive area near/6 has the capacityThe sex and the sensitivity are difficult to compensate for in a targeted manner, as shown in fig. 10 (b).
(4) Analyzing the influence of the uncontrolled rectifying load on the stability of the SAPF system, wherein the specific analysis process is as follows:
the uncontrolled rectifying load is widely applied in the industry, when the uncontrolled rectifying load is put into operation, harmonic waves can be introduced into a power grid, and meanwhile, the input admittance of the uncontrolled rectifying load is coupled with an SAPF control loop to generate harmonic wavesY mL (s) From the foregoing analysis, it can be seen that:
Figure 70312DEST_PATH_IMAGE105
after the sampling process has been linearized,Y mL (s) Can be simplified into a single frequency admittance model
Figure 769147DEST_PATH_IMAGE106
Namely:
Figure 718648DEST_PATH_IMAGE107
wherein the content of the first and second substances,
Figure 745510DEST_PATH_IMAGE108
is a power grid current loop open-loop transfer function under a single-frequency model,G D (s)=G h (s)D(z). According to the formula, the compound has the advantages of,
Figure 196083DEST_PATH_IMAGE109
respectively amplitude and phase angle ofY L (s) And
Figure 624790DEST_PATH_IMAGE110
the decision, can be expressed as follows:
Figure 869827DEST_PATH_IMAGE111
Y L (s) The expression of input impedance of (a) is:
Figure 649564DEST_PATH_IMAGE112
wherein, the first and the second end of the pipe are connected with each other,Z dc for dc-side load impedance, it is common in industrial applications to exhibit a resistive characteristic,f 1 is the fundamental frequency. The expression showsY L (s) At medium and high frequencies withZ dc Exhibit the same resistance to inductance. For the
Figure 513615DEST_PATH_IMAGE113
When the LCL parameter of filter II is used, andK ph is 2.83 the number of the omega is zero,K pf at 0.20 Ω, the Bode diagram is shown in fig. 11. The equivalent admittance of the SAPF system for the PCC port is as follows, regardless of the PCC voltage feed forward effectY mL ()+Y oca (). At the resonant frequencyf r Nearby, though
Figure 422927DEST_PATH_IMAGE114
Has a large phase change, but has a large harmonic amplitude value, so that
Figure 979811DEST_PATH_IMAGE115
Has a small amplitude and a relative phasorY oca () The change of (a) is slight; during the parameter design process, the requirements are always met
Figure 981265DEST_PATH_IMAGE116
And is
Figure 508061DEST_PATH_IMAGE117
Therefore, it is
Figure 911360DEST_PATH_IMAGE118
The value of (A) is positive and real, then always satisfy
Figure 763779DEST_PATH_IMAGE119
Thus, is atω s Near/6
Figure 986950DEST_PATH_IMAGE120
Is mainly determined byY L (s). As shown in (a) of FIG. 12, inω s Near/6 frequency, equivalent phasorY mL () When the fan-shaped shadow area is positioned, the fan-shaped shadow area can be effectively reduced or even eliminatedY oca () In thatω s And a non-passive area near the frequency of 6 improves the robustness of the SAPF system to the power grid impedance fluctuation. However, due to the fluctuation of the load, the passivity of the SAPF output admittance cannot be realized only by depending on the regulation of the load when the load is emptyY mL () =0 orY mL () In the non-shaded region of the fourth quadrant,Y oca () In thatω s There will still be non-passive regions around the/6 frequency.
(5) Designing PCC voltage feedforward proportional gainK ff The robustness of the SAPF system is improved, and the specific design process is as follows:
in order to eliminate the influence of the load as much as possible, the load can be reducedY oca () In thatω s Phase around/6 frequency, in the present invention byY ocf (s) To pairY oca (s) Is adjusted. ByY oc (s) = Y oca (s) + Y ocf (s) Can find out
Figure 458382DEST_PATH_IMAGE121
The expression is as follows:
Figure 840822DEST_PATH_IMAGE122
as can be seen from the above equation, when the first order low pass filterG ff (z) When a proper cut-off frequency is selected, the frequency can be effectively reducedY oca (s) Phase in the mid-high frequency band.G ff (z) Correspond tosOf domainsG ff (s) The transfer function of (c) can be expressed as follows:
Figure 473929DEST_PATH_IMAGE123
the above equation is a typical first-order low-pass filter,ω f is the cut-off angular frequency. To reduceY oca (s) In thatω s A phase around a/6 frequency, which can be selectedω f Is (rad/s).
In the present embodiment, drawing is adopted
Figure 449975DEST_PATH_IMAGE124
Is dependent on the phaseK ff Design parameters with varying frequency response curvesK ff As shown in fig. 13. As can be seen from figure 13 of the drawings,K ff main influence of
Figure 82688DEST_PATH_IMAGE125
Frequency characteristics in a low frequency band. In order to ensure that certain feedforward is used for improving the dynamic performance of APF during grid connection and load switching, the selection is carried outK ff Is 0.8, at this timeY oc (s) Is plotted in fig. 14, from fig. 14 it can be found thatY oc (s) The phase of the medium and high frequency band is obviously reduced, and when the 280uH power grid impedance is switched, the system cannot be in resonance instability with the power grid impedance. At this time from (b) in figure 12,Y mL () The feasible shadow sector area is obviously widened, and at the moment, the load is switched,Y mL () When the change occurs in the shadow sector, the total output admittance can be effectively ensuredY oa (s) The passivity of (A).
(1) The experiment verifies as follows:
the experiment adopts the parameters of the LCL filter II in the table 2, and the designed parameters are as follows according to the steps:K MOPT is a content of 0.0696 (g),ω c at 2000 (rad/s),K ph the content of the acid was 2.83,K pf is 0.20.PCC voltage feedforward linkω f Is (rad/s),K ff is 0.8。
With the above parameters, when the grid impedance is not switched, the system steady state performance is as shown in (a) and (b) in fig. 15, where (a) in 15 is abc three-phase grid current waveform, and (b) in 15 is corresponding Total Harmonic Distortion (THD) value. As can be seen from fig. 15, under the control of the SAPF dual-loop current, the THD of the grid current can be reduced to less than 5% by using the parameters in the present embodiment, and the industry standard is satisfied.
In order to verify the effect of the PCC voltage feedforward link, as shown in fig. 16, when the PCC voltage feedforward link is not applied, the grid impedance of 600 uH is applied, the system resonates, and after a load is applied, the resonance is obviously suppressed. The nonlinear load can effectively improve the robustness of the SAPF system to power grid impedance fluctuation to a certain extent, but in practice, the load is switched frequently, and the SAPF system is unstable again after the load is cut off.
After the PCC voltage feed-forward is applied, as shown in (a) of fig. 17, after the grid impedance is applied, the system is stable, and the corresponding grid current THD is 4.3% as shown in (b) of fig. 17, which meets the industry standard; after the load is removed, the system remains stable. The robustness of the SAPF system to the power grid impedance fluctuation is greatly improved after the PCC voltage feedforward is applied, and the whole process is not influenced by load switching.
Experimental results show that the parameter design method can ensure the stable operation of the three-phase LCL type SAPF when the power grid impedance and the load fluctuate, and achieve the ideal harmonic compensation effect.
The embodiments described above are intended to facilitate one of ordinary skill in the art in understanding and using the invention. It will be readily apparent to those skilled in the art that various modifications to the above-described embodiments may be made, and the generic principles defined herein may be applied to other embodiments without the use of inventive faculty. Therefore, the present invention is not limited to the above embodiments, and those skilled in the art should make improvements and modifications to the present invention based on the disclosure of the present invention within the protection scope of the present invention.

Claims (8)

1. An SAPF parameter design method considering load and power grid impedance influence, wherein the SAPF parameters comprise proportional link gain of a fundamental current loop and a power grid current loop controller, LCL filter parameters and a proportional gain coefficient of a PCC power grid voltage feedforward term;
the parameter design method is characterized by comprising the following steps:
(1) Obtaining the optimal proportional coefficient expression of the proportional link gain of a fundamental current loop controller and the proportional link gain of a power grid current loop controller under passive guidance according to the real part information of the system output admittance when PCC voltage feedforward is not applied in the power grid current response equation of the SAPF double-loop current control system;
(2) Calculating a correction factor of the optimal proportionality coefficient expression in the step (1) according to the maximum negative deviation of the inductance value and the capacitance value at the side of the inverter in the LCL filter to obtain a corrected optimal proportionality coefficient expression;
(3) Obtaining an inequality equation set of LCL filter parameter values according to the corrected optimal proportionality coefficient expression, the stable condition of the power grid current response equation of the SAPF double-loop current control system and the given minimum control bandwidth of the power grid current loop, and selecting the LCL filter parameters in a feasible area surrounded by the inequality equation set;
(4) Substituting the LCL filter parameter value selected in the step (3) into the optimal proportionality coefficient expression corrected in the step (2) to obtain a corrected optimal proportionality coefficient, and obtaining the proportional link gain of the power grid current loop controller and the proportional link gain of the fundamental current loop controller according to the corrected optimal proportionality coefficient, the selected minimum power grid current loop control bandwidth and the power grid side inductance value of the LCL filter;
(5) And (5) designing a proportional gain coefficient of a voltage feedforward term of the PCC power grid, which meets the passivity requirement, under the condition of the proportional link gain of the power grid current loop controller and the proportional link gain of the fundamental current loop controller obtained in the step (4).
2. The method of claim 1, wherein the SAPF double-loop current control system comprises:
the bus voltage loop takes direct-current bus voltage and a direct-current bus voltage reference value as input, and obtains the amplitude of the fundamental current loop reference value after being adjusted by a bus voltage controller; obtaining a phase angle of a fundamental wave current loop reference value after the phase-locked loop is used for carrying out phase locking on the PCC network voltage; obtaining a fundamental current loop reference value according to the amplitude and the phase angle;
the fundamental current loop takes a fundamental current loop reference value and inverter side feedback current as input, and obtains output quantity after passing through the fundamental current loop controller;
the power grid current loop takes power grid current as input and obtains output quantity after passing through a power grid current loop controller;
the PCC power grid voltage feedforward item takes the PCC power grid voltage as input, and obtains output quantity after low-pass filtering and multiplication by a proportional gain coefficient;
and the output quantities of the fundamental wave current loop, the grid current loop and the PCC grid voltage feedforward term are added and then are modulated by SPWM to obtain a switching signal, and the switching signal acts on the three-phase voltage source inverter to obtain the output voltage of the inverter.
3. The SAPF parameter design method considering load and grid impedance influences of claim 1, wherein a grid current response equation of the SAPF double-loop current control system is as follows:
Figure 722485DEST_PATH_IMAGE001
wherein the content of the first and second substances,i g (s)、u g (s)、i Lh (s) Harmonic current sources equivalent to the power grid current, the power grid voltage and the nonlinear load respectively;G L (s) Is a closed loop transfer function from the non-linear load current to the grid current;Y oa (s) For the total output admittance of the system,Y mL (s) The admittance generated for the nonlinear load input admittance coupling to the system,Y oc (s) For the system output admittance when applying the PCC voltage feed-forward,Y oca (s) To the output admittance when no PCC voltage feed-forward is applied,Y ocf (s) For the PCC voltage feed-forward effect equivalent admittance,Z g (s) Is the grid impedance.
4. The method for designing SAPF parameters in consideration of load and grid impedance influences according to claim 1, wherein the modified optimal scaling factor in step (2) is expressed as follows:
Figure 108467DEST_PATH_IMAGE002
wherein the content of the first and second substances,K OPT in order to optimize the scaling factor before the correction,K MOPT in order to obtain the optimum scaling factor after the modification,
Figure 980214DEST_PATH_IMAGE003
in order to correct the factor(s),ω s in order to sample the angular frequency of the signal,ω r being the resonance angular frequency of the LCL filter,
Figure 844265DEST_PATH_IMAGE004
is the ratio of the inverter side inductance value and the grid side inductance value in the LCL filter.
5. The SAPF parameter design method considering the influence of load and grid impedance of claim 4, wherein the correction factors are:
Figure 252113DEST_PATH_IMAGE005
wherein the content of the first and second substances,
Figure 543417DEST_PATH_IMAGE006
inverter side inductance values in LCL filter andthe maximum negative deviation between the capacitance value and the factory parameters.
6. The SAPF parameter design method considering load and grid impedance influence according to claim 4, wherein the inequality equation set of LCL filter parameter values in step (3) includes the following constraints:
the first constraint condition is:
Figure 403925DEST_PATH_IMAGE007
the second constraint condition is as follows:
Figure 71667DEST_PATH_IMAGE009
the third constraint condition is as follows:
Figure 334021DEST_PATH_IMAGE010
wherein the content of the first and second substances,ω cmin for a given minimum control bandwidth of the grid current loop,T s is the sampling time.
7. The SAPF parameter design method considering load and grid impedance influence according to claim 4, wherein in step (4), the calculation formulas of the grid current loop controller proportional link gain and the fundamental current loop controller proportional link gain are as follows:
Figure DEST_PATH_IMAGE011
wherein the content of the first and second substances,ω cmin for a given minimum control bandwidth of the grid current loop,K pf is the proportional link gain of the fundamental current loop controller,K ph the proportional link gain of the current loop controller of the power grid is obtained,L 2 is LCLGrid side inductance values in the filter.
8. The SAPF parameter design method considering load and grid impedance influence as claimed in claim 1, wherein the step (5) is specifically as follows:
and (3) according to the value range of the proportional gain coefficient of the PCC power grid voltage feedforward term, the gain of the proportional link of the power grid current loop controller obtained in the step (5) and the gain of the proportional link of the fundamental current loop controller, obtaining the frequency response of the system output admittance when the PCC voltage feedforward is applied, and selecting the proportional gain coefficient of the PCC power grid voltage feedforward term as a design result in the range that the frequency response meets the passivity requirement.
CN202211645210.7A 2022-12-21 2022-12-21 SAPF parameter design method considering load and power grid impedance influence Active CN115632401B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211645210.7A CN115632401B (en) 2022-12-21 2022-12-21 SAPF parameter design method considering load and power grid impedance influence

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211645210.7A CN115632401B (en) 2022-12-21 2022-12-21 SAPF parameter design method considering load and power grid impedance influence

Publications (2)

Publication Number Publication Date
CN115632401A CN115632401A (en) 2023-01-20
CN115632401B true CN115632401B (en) 2023-03-07

Family

ID=84910806

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211645210.7A Active CN115632401B (en) 2022-12-21 2022-12-21 SAPF parameter design method considering load and power grid impedance influence

Country Status (1)

Country Link
CN (1) CN115632401B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116760108B (en) * 2023-08-21 2024-01-05 浙江浙能技术研究院有限公司 LCL-SAPF stability control method based on unified stability constraint
CN117914172B (en) * 2024-03-20 2024-05-17 河海大学 Method, equipment and medium for designing control parameters of voltage loop of grid-formed inverter

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102510120A (en) * 2011-11-23 2012-06-20 中国科学院电工研究所 Micro-grid inverter voltage and current double-ring hanging control method based on virtual impedance
CN110768299A (en) * 2019-12-02 2020-02-07 湖南大学 Sequence impedance modeling and stability analysis method of load virtual synchronous machine
CN113629984A (en) * 2021-07-27 2021-11-09 浙江大学 Three-phase LCL type SAPF parameter design method based on double-loop current control strategy
CN114513006A (en) * 2022-02-18 2022-05-17 浙江浙能技术研究院有限公司 LCL type SAPF robust current control method based on adaptive IIR filtering link

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11177732B1 (en) * 2020-09-04 2021-11-16 Rockwell Automation Technologies, Inc. Dynamic resonance control for grid-tied power converters

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102510120A (en) * 2011-11-23 2012-06-20 中国科学院电工研究所 Micro-grid inverter voltage and current double-ring hanging control method based on virtual impedance
CN110768299A (en) * 2019-12-02 2020-02-07 湖南大学 Sequence impedance modeling and stability analysis method of load virtual synchronous machine
CN113629984A (en) * 2021-07-27 2021-11-09 浙江大学 Three-phase LCL type SAPF parameter design method based on double-loop current control strategy
CN114513006A (en) * 2022-02-18 2022-05-17 浙江浙能技术研究院有限公司 LCL type SAPF robust current control method based on adaptive IIR filtering link

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Guanghui Yang.A Novel Sequential Direct Torque Control Scheme for Seven-phase Induction Motors.2020,全文. *
孟德智.基于ESO的永磁同步电机转子动能反馈控制策略.2021,全文. *
齐虹.带LCL 滤波的SAPF 电流控制点选取的研究.2016,全文. *

Also Published As

Publication number Publication date
CN115632401A (en) 2023-01-20

Similar Documents

Publication Publication Date Title
CN115632401B (en) SAPF parameter design method considering load and power grid impedance influence
JP4230648B2 (en) Controller for reactive power series compensator
CN110112940B (en) PWM rectifier self-adaptive sliding mode QPIR control method based on αβ coordinate system
CN108667068B (en) Method for realizing hybrid damping of LCL grid-connected inverter based on PC-QPCI
CN107346946A (en) A kind of inverter is discrete to repeat sliding-mode control
CN108390413B (en) Method for improving adaptability of distributed photovoltaic to alternating current-direct current hybrid micro-grid
CN113629984B (en) Three-phase LCL type SAPF parameter design method based on double-loop current control strategy
CN111313474B (en) Improved micro-grid-connected presynchronization control method
CN110729752B (en) Output impedance remolding method of grid-connected inverter parallel system
CN107895966A (en) The light current electric voltage feed forward lag compensation control method off the net based on impedance self-adaptive
CN109193792A (en) The SOCVF feedforward of gird-connected inverter and phase compensating factor method for handover control under weak grid
CN109524992B (en) Grid-connected converter self-adaptive control method and control system under weak grid
CN114865633A (en) Self-adaptive quasi-PR active damping low-frequency harmonic suppression method
CN109327048B (en) Robust phase locking system and method for grid-connected converter
CN113489292B (en) Control method for negative feedback virtual impedance of LCL type grid-connected converter
CN109950926B (en) Grid-connected inverter stability control method based on q-axis voltage integral feedforward under weak network
CN117335677A (en) Voltage stabilization control method for energy storage converter
CN117318079A (en) Three-level parallel static var generator based on LCLLC filtering
Xu et al. The stability of LCL-type grid-tied inverter based on repetitive control and grid voltage feed-forward
CN115207968B (en) Method for improving stability margin of photovoltaic grid-connected inverter under weak current network
CN114285081B (en) Wide area power system stabilizing method based on self-adaptive virtual resistor
CN115498657A (en) Virtual admittance method for improving stability of grid-connected inverter
CN115378040A (en) Grid-connected system based on LCL type photovoltaic inverter and QVR control method
CN114552666A (en) Phase-locked loop compensation control circuit based on first-order complex vector filter under weak power grid
CN116599066B (en) Mixed APF flexibility improving method based on voltage loop gain locking

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