CN105162159B - Photovoltaic grid inverter controller parameter identification method based on differential evolution algorithm - Google Patents
Photovoltaic grid inverter controller parameter identification method based on differential evolution algorithm Download PDFInfo
- Publication number
- CN105162159B CN105162159B CN201510472410.0A CN201510472410A CN105162159B CN 105162159 B CN105162159 B CN 105162159B CN 201510472410 A CN201510472410 A CN 201510472410A CN 105162159 B CN105162159 B CN 105162159B
- Authority
- CN
- China
- Prior art keywords
- population
- vector
- value
- formula
- parameter
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 21
- 230000008859 change Effects 0.000 claims abstract description 23
- 230000008569 process Effects 0.000 claims description 9
- 230000001360 synchronised effect Effects 0.000 claims description 5
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000010355 oscillation Effects 0.000 claims 1
- 230000001052 transient effect Effects 0.000 abstract description 7
- 230000007812 deficiency Effects 0.000 abstract description 4
- 238000004364 calculation method Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 6
- 230000004224 protection Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000012850 discrimination method Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
- Y02E10/56—Power conversion systems, e.g. maximum power point trackers
Abstract
The invention discloses a kind of photovoltaic grid inverter controller parameter identification method based on differential evolution algorithm, using differential evolution algorithm to photovoltaic combining inverter current inner loop d axles control system and the proportionality coefficient and integral coefficient of q axle control system proportional and integral controllers, the parameter such as integrator output violent change value and controller output violent change value limit is recognized, the initial value of each solution vector in initial population is determined by theoretical calculation method, effectively cover and entirely search prime field, reduce the hunting zone of differential evolution algorithm, accelerate the speed that parameter to be identified converges to actual value；The identification flow stage by stage of pi regulator parameter and amplitude limit link parameter is established, to reduce parameter identification dimension, improves controller parameter identification efficiency.Accurate recognition goes out the parameter of controller amplitude limit link, to accurately reflect the fault transient characteristic of photo-voltaic power supply, makes up the deficiency of existing controller parameter identification method, has a good application prospect.
Description
Technical field
The present invention relates to power system modeling technical field, and in particular to a kind of based on the grid-connected of differential evolution algorithm
Circuit control device parameter identification method.
Background technology
Photovoltaic power generation technology has the advantages that not limited by regional condition, scaleable, safe and reliable, pollution-free, thus
Rapid development and popularization are obtained in China, grid connection capacity is growing.But the fault current characteristics that photo-voltaic power supply is complicated
Huge challenge also is brought to traditional protection philosophy and fault detection method with intermittent, randomness output feature, is led
Cause its possibly can not failure judgement exactly position, cause relay fail or malfunction.
It is known that the basic task of relay protection is become according to caused electrical quantity (or non-electric quantity) after electric network fault
Change feature, discriminating and isolated fault equipment or element, therefore, electric network fault specificity analysis is the basis of relay protection research field
Content.
Because photo-voltaic power supply is sound of the combining inverter during electric network fault by combining inverter and grid interface
Characteristic is answered to determine the fault current characteristics of photo-voltaic power supply, and the fault response characteristics of photovoltaic combining inverter then depend on its control
The selection of device parameter processed.Therefore, the controller parameter of accurate recognition photovoltaic combining inverter (includes the ratio of proportional and integral controller
Example coefficient, integral coefficient and amplitude limit link parameter etc.), there is important meaning for the photo-voltaic power supply equivalent model for building high confidence level
Justice, solid foundation is established to study the electric network fault current characteristics of the access containing photo-voltaic power supply and protection philosophy research.
At present, the controller parameter discrimination method on photovoltaic combining inverter or other kinds of transverter, which is intended to, distinguishes
Know the proportionality coefficient and integral coefficient of controller inner and outer rings proportional and integral controller, do not recognize the output violent change of integrator
The parameter of the amplitude limit link such as value and controller output violent change value, is specifically described as follows,
As shown in figure 1, grid-connected photovoltaic power supply by photovoltaic array, DC bus capacitor, combining inverter, LC wave filter groups into,
The mathematical equation of photovoltaic combining inverter is as follows under two-phase synchronous rotary dq coordinate systems,
Wherein, ud、uqFor d, q axis component of bridge arm output voltage, ugd、ugqFor d, q axis component of line voltage, id、iqFor
D, q axis component of inverter output current, ω1For synchronous angular velocity, L1For AC filter inductance, R1For filter inductance etc.
Imitate resistance；
Using grid voltage orientation vector controlled, d axles are oriented to line voltage vector, and take feedforward compensation and PI to adjust
During the control mode of section, the double -loop control block diagram of photovoltaic combining inverter is as shown in Fig. 2 figure it is seen that system employs
Double closed-loop control system, outer shroud are power voltage ring, and inner ring is electric current loop, in order to realize that the unity power factor of photo-voltaic power supply is transported
OK, Q is typically taken*=0.PI controller of the outer shroud instruction compared with value of feedback, through overpower-outer voltage is obtained in electric current
Fourth finger makesWithWithWith the PI controllers compared with feedback current, by current inner loop, and feedforward benefit is carried out
Repay, obtain the three-phase voltage control instruction needed for inverter after coordinate transformWith
When a grid fault occurs, because trouble duration is general shorter, and photo-voltaic power supply is typically arranged with energy storage dress
Put, thus, it is believed that photovoltaic combining inverter DC side electric current capacitance voltage UdcIt is held essentially constant during failure.Meanwhile
For the response speed of photo-voltaic power supply during improving electric network fault, power voltage outer shroud locking can be carried out straight to current inner loop
Connect control.Therefore, in fault transient analysis, the influence of current inner loop controller operation characteristic need to only be considered, shown in Fig. 2
In the double -loop control block diagram of photovoltaic combining inverter, the output current limitation of photovoltaic combining inverter is not considered.However, by
Influenceed in by factors such as inverter itself capacity and overcurrent protections, the output current of photovoltaic combining inverter must be limited
System.Therefore, it is necessary to output violent change to controller during control strategy specific implementation.Amplitude limit is carried out being exported to controller
In the case of, if grid collapses cause the command value of current inner loop a larger step value occur, controller output
It can be limited by saturation, this will cause the hydraulic performance decline of actual closed-loop system, and such as overshoot increase, stabilization time increase, sternly
Stability of control system is even destroyed in the case of weight.It is this due to controller output violent change cause system closed loop response be deteriorated show
As referred to as windup phenomenons.In order to avoid the appearance of controller windup phenomenons, effective method most directly perceived is exactly to design
Anti-windup controllers, suppress the saturated phenomenon of controller, and allow system to exit saturation region as early as possible.Wherein, it is most simple
Single anti-windup controller design methods are that the output to the integral element in PI controllers carries out amplitude limit.
By taking d axles as an example, after considering controller output violent change and anti-windup designs (integrator output violent change), electricity
The specific implementation link of inner ring control system is flowed as shown in figure 3, due to the presence of amplitude limit link, and photovoltaic is simultaneously after grid collapses
Relation between the input of the controller of net inverter, output can not be described with analytical expression.
By foregoing description, when a grid fault occurs, the disturbance of line voltage may cause photovoltaic combining inverter
Controller enters nonlinear element, influences integrator output and controller output violent change, and this will largely effect on the event of photo-voltaic power supply
Hinder transient characterisitics.Therefore, the photo-voltaic power supply constructed using the controller parameter discrimination method of conventional photovoltaic combining inverter
Equivalent model, the fault transient characteristic of photo-voltaic power supply can not be accurately reflected, how to make up existing controller parameter identification method
Deficiency, accurate recognition go out the parameter of controller amplitude limit link, are current urgent problems.
The content of the invention
Technical problem solved by the invention is to overcome to recognize using the controller parameter of traditional photovoltaic combining inverter
The equivalent model for the photo-voltaic power supply that method constructs, the problem of fault transient characteristic of photo-voltaic power supply can not be accurately reflected.
In order to achieve the above object, the technical solution adopted in the present invention is：
A kind of photovoltaic grid inverter controller parameter identification method based on differential evolution algorithm, it is characterised in that：Bag
Include following steps,
Step (1), population scale NP=50 is set, entering for moment t1, t2 is respectively obtained according to formula (1) and formula (2)
Change colonyWith
Wherein,The proportionality coefficient of d, q shaft current inner ring pi regulator in i-th of population of t1 moment is corresponded to respectively
kP, integral coefficient kI, i=1,2 ..., NP；WithD, q shaft current in i-th of population of t2 moment are corresponded to respectively
Inner ring integrator output violent change value Lint_up、Lint_low, and current inner loop controller output violent change value Lout_up、Lout_low, i=
1,2 ..., NP；
T1=0, t2=0 are made, initial population is setIn section [kP_1/ 5,5kP_2] interior value,
In section [kI_1/ 5,5kI_2] interior value,Section (0,2.0] interior value,Section [- 2.0,
0] interior value, according to formula (3), k is calculatedP_1、kP_2、kI_1And kI_2Value,
ωcFor cut-off angular frequency, value is 10 ωs, ωsFor synchronous rotary angular speed；ξ is damping ratio, ωnTo shake naturally
Swing frequency, ξ=0.707, ωn=10 ωs；L1For the filter inductance of photovoltaic combining inverter AC, R1For filter inductance L1's
Internal resistance；
According to formula (4), correspondence function is listed
Wherein, idAnd i (k)q(k) it is the photovoltaic combining inverter current inner loop controller in the case of current inner loop command value
Reality output；WithFor in the case of identical current inner loop command value, photovoltaic combining inverter current inner loop control
The output of device trace model；For solution vector；N=48, it is
The sampling number of one power frequency period；
Grid entry point voltage is set to drop to 0.85pu from 1.0pu so that current inner loop controller is operated in the range of linearity, and
Perform step (2)；
Step (2), if t1>NG, NG are the identification of maximum evolutionary generation, the then proportionality coefficient and integral coefficient of pi regulator
Process terminates, populationThe minimum solution vector of middle correspondence function value is the proportionality coefficient and integration of the pi regulator picked out
Coefficient, and perform step (7)；Otherwise, step (3) is performed；
Step (3), utilizes correspondence functionTo populationIn each solution vector evaluated, if kind
GroupThe middle matching condition that solution vector satisfaction be present as shown in formula (5), then parameter identification success,
Wherein, ε values are 0.01；Corresponding solution vector is the proportionality coefficient and integral coefficient of the pi regulator picked out,
And perform step (7)；If populationIn each solution vector be unsatisfactory for matching condition as shown in formula (5), then perform step
Suddenly (4)；
Step (4), to populationIn solution vectorMutation operation is carried out according to formula (6), obtains variation vector
Wherein, r1, r2, r3 ∈ { 1,2,3 ..., NP }, not equal to i, and are not mutually equal two-by-two；F is zoom factor, is taken
It is [0,1] to be worth scope；
Step (5), by populationIn solution vectorWith variation vectorCarry out the intersection as shown in formula (7)
Operation, produce trial vector
Wherein, rand (j) ∈ [0,1], it is equally distributed random number；J ∈ { 1,2 }, represent j-th of variable；randn
(i) ∈ { 1,2 }, indexed for randomly selected dimension variable；CR is crossover probability constant, and span is [0,1]；
Step (6), by making a variation and caused trial vector after crossover operationWith populationIn solution vectorEnter
Competition of the row as shown in formula (8), and select the less vector of correspondence function value to be used as filial generation, obtain the population at t1+1 moment
Make t1=t1+1, return to step (2)；
Step (7), grid entry point voltage is set to drop to 0.4pu from 1.0pu, so that current inner loop controller is operated in non-
The range of linearity, and perform step (8)；
Step (8), if t2>NG, then amplitude limit link parameter identification process terminate, populationMiddle correspondence function value is minimum
Solution vector be the amplitude limit link parameter picked out；Otherwise, step (9) is performed；
Step (9), if populationMiddle have solution vector and meet matching condition as shown in formula (5), then amplitude limit link ginseng
Number recognizes successfully, and identification process terminates, and corresponding solution vector is the amplitude limit link parameter picked out；If populationIn it is each
Solution vector is unsatisfactory for the matching condition as shown in formula (5), then performs step (10)；
Step (10), to populationIn solution vectorMutation operation is carried out according to formula (9), obtains becoming incorgruous
Amount
Wherein, r1, r2, r3 ∈ { 1,2,3 ..., NP }, not equal to i, and are not mutually equal two-by-two；F is zoom factor, is taken
It is [0,1] to be worth scope；
Step (11), by populationIn solution vectorWith variation vectorCarry out as shown in formula (10)
Crossover operation, produce trial vector
Wherein, rand (j) ∈ [0,1], it is equally distributed random number, j ∈ { 3,4,5,6 }, represents j-th of variable；
Randn (i) ∈ { 3,4,5,6 }, indexed for randomly selected dimension variable；CR is crossover probability constant, span for [0,
1]；
Step (12), by making a variation and caused trial vector after crossover operationWith populationIn solution vectorThe competition as shown in formula (11) is carried out, and selects the less vector of correspondence function value to be used as filial generation, when obtaining t2+1
The population at quarter
Make t2=t2+1, return to step (8).
The foregoing photovoltaic grid inverter controller parameter identification method based on differential evolution algorithm, it is characterised in that：
Step (2), the maximum evolutionary generation NG of step (8) are 1000.
The foregoing photovoltaic grid inverter controller parameter identification method based on differential evolution algorithm, it is characterised in that：
Step (4), step (10) zoom factor F preferred value are 0.5.
The foregoing photovoltaic grid inverter controller parameter identification method based on differential evolution algorithm, it is characterised in that：
Step (5), step (11) crossover probability constant CR preferred value are 0.75.
The beneficial effects of the invention are as follows：Photovoltaic combining inverter current inner loop d axles are controlled using differential evolution algorithm and are
The proportionality coefficient and integral coefficient, integrator output violent change value and controller of system and q axle control system proportional and integral controllers are defeated
Go out the parameters such as amplitude limit value to be recognized, the initial value of each solution vector in initial population is determined by theoretical calculation method, is effectively covered
Cover and entirely searched prime field, reduced the hunting zone of differential evolution algorithm, accelerated the speed that parameter to be identified converges to actual value
Degree；The identification flow stage by stage of pi regulator parameter and amplitude limit link parameter is established, to reduce parameter identification dimension, is improved
Controller parameter recognizes efficiency.Accurate recognition goes out the parameter of controller amplitude limit link, to accurately reflect the failure of photo-voltaic power supply
Transient characterisitics, the deficiency of existing controller parameter identification method is made up, is had a good application prospect.
Brief description of the drawings
Fig. 1 is the system block diagram of grid-connected photovoltaic power supply.
Fig. 2 is the control block diagram of photovoltaic combining inverter double -loop control.
Fig. 3 is the control block diagram of current inner loop d axle control systems.
Fig. 4 is the principle of the photovoltaic grid inverter controller parameter identification method based on differential evolution algorithm of the present invention
Figure.
Embodiment
Below in conjunction with Figure of description, the present invention is further illustrated.
The photovoltaic grid inverter controller parameter identification method based on differential evolution algorithm of the present invention, is entered using difference
Change algorithm to photovoltaic combining inverter current inner loop d axles control system and the ratio system of q axle control system proportional and integral controllers
The parameter such as number and integral coefficient, integrator output violent change value and controller output violent change value is recognized, by theoretical calculation side
Method determines the initial value of each solution vector in initial population, effectively covers and entirely searches prime field, reduces searching for differential evolution algorithm
Rope scope, accelerate the speed that parameter to be identified converges to actual value；Establish pi regulator parameter and amplitude limit link parameter
Flow is recognized stage by stage, to reduce parameter identification dimension, improves controller parameter identification efficiency.Accurate recognition goes out controller limit
The parameter of width link, to accurately reflect the fault transient characteristic of photo-voltaic power supply, make up existing controller parameter identification method
Deficiency, schematic diagram is as shown in figure 4, according to the output of actual photovoltaic combining inverter and photovoltaic grid-connected inversion under grid fault conditions
The difference of the output of device current inner loop control system trace model, trace model is constantly repaiied by differential evolution algorithm
Just, so as to picking out controller parameter, following steps are specifically included,
Step (1), population scale NP=50 is set, entering for moment t1, t2 is respectively obtained according to formula (1) and formula (2)
Change colonyWith
Wherein,(i=1,2 ..., NP) corresponds to d, q shaft current inner ring PI in i-th of population of t1 moment and adjusted respectively
The proportionality coefficient k of deviceP, integral coefficient kI；With(i=1,2 ..., NP) corresponds to i-th of population of t2 moment respectively
Middle d, q shaft current inner ring integrator output violent change value Lint_up、Lint_low, and current inner loop controller output violent change value
Lout_up、Lout_low；
T1=0, t2=0 are made, initial population is setIn section [kP_1/ 5,5kP_2] interior value,
In section [kI_1/ 5,5kI_2] interior value,Section (0,2.0] interior value,Section [- 2.0,
0] interior value, according to formula (3), k is calculatedP_1、kP_2、kI_1And kI_2Value,
ωcFor cut-off angular frequency, value is 10 ωs, ωsFor synchronous rotary angular speed；ξ is damping ratio, ωnTo shake naturally
Swing frequency, ξ=0.707, ωn=10 ωs；L1For the filter inductance of photovoltaic combining inverter AC, R1For filter inductance L1's
Internal resistance；
According to formula (4), correspondence function is listed
Wherein, idAnd i (k)q(k) it is the photovoltaic combining inverter current inner loop controller in the case of current inner loop command value
Reality output；WithFor in the case of identical current inner loop command value, photovoltaic combining inverter current inner loop control
The output of device trace model processed；For solution vector；N=48,
For the sampling number of a power frequency period；
Grid entry point voltage is set to drop to 0.85pu from 1.0pu so that current inner loop controller is operated in the range of linearity, and
Perform step (2)；
Step (2), if t1>NG, NG are the identification of maximum evolutionary generation, the then proportionality coefficient and integral coefficient of pi regulator
Process terminates, populationThe minimum solution vector of middle correspondence function value is the proportionality coefficient and integration of the pi regulator picked out
Coefficient, and perform step (7)；Otherwise, step (3) is performed；
Step (3), utilizes correspondence functionTo populationIn each solution vector evaluated, if
PopulationThe middle matching condition that solution vector satisfaction be present as shown in formula (5), then parameter identification success,
Wherein, ε values are 0.01；Corresponding solution vector is the proportionality coefficient and integral coefficient of the pi regulator picked out,
And perform step (7)；If populationIn each solution vector be unsatisfactory for matching condition as shown in formula (5), then perform step
Suddenly (4)；
Step (4), to populationIn solution vectorMutation operation is carried out according to formula (6), obtains variation vector
Wherein, r1, r2, r3 ∈ { 1,2,3 ..., NP }, not equal to i, and are not mutually equal two-by-two；F is zoom factor, is taken
It is [0,1] to be worth scope；
Step (5), by populationIn solution vectorWith variation vectorCarry out the intersection as shown in formula (7)
Operation, produce trial vector
Wherein, rand (j) ∈ [0,1], it is equally distributed random number；J ∈ { 1,2 }, represent j-th of variable；randn
(i) ∈ { 1,2 }, indexed for randomly selected dimension variable；CR is crossover probability constant, and span is [0,1]；
Step (6), by making a variation and caused trial vector after crossover operationWith populationIn solution vectorEnter
Competition of the row as shown in formula (8), and select the less vector of correspondence function value to be used as filial generation, obtain the population at t1+1 moment
Make t1=t1+1, return to step (2)；
Step (7), grid entry point voltage is set to drop to 0.4pu from 1.0pu, so that current inner loop controller is operated in non-
The range of linearity, and perform step (8)；
Step (8), if t2>NG, then amplitude limit link parameter identification process terminate, populationMiddle correspondence function value is minimum
Solution vector be the amplitude limit link parameter picked out；Otherwise, step (9) is performed；
Step (9), if populationMiddle have solution vector and meet matching condition as shown in formula (5), then amplitude limit link ginseng
Number recognizes successfully, and identification process terminates, and corresponding solution vector is the amplitude limit link parameter picked out；If populationIn it is each
Solution vector is unsatisfactory for the matching condition as shown in formula (5), then performs step (10)；
Step (10), to populationIn solution vectorMutation operation is carried out according to formula (9), obtains becoming incorgruous
Amount
Wherein, r1, r2, r3 ∈ { 1,2,3 ..., NP }, not equal to i, and are not mutually equal two-by-two；F is zoom factor, is taken
It is [0,1] to be worth scope；
Step (11), by populationIn solution vectorWith variation vectorCarry out as shown in formula (10)
Crossover operation, produce trial vector
Wherein, rand (j) ∈ [0,1], it is equally distributed random number；J ∈ { 3,4,5,6 }, represent j-th of variable；
Randn (i) ∈ { 3,4,5,6 }, indexed for randomly selected dimension variable；CR is crossover probability constant, span for [0,
1]；
Step (12), by making a variation and caused trial vector after crossover operationWith populationIn solution vectorThe competition as shown in formula (11) is carried out, and selects the less vector of correspondence function value to be used as filial generation, when obtaining t2+1
The population at quarter
Make t2=t2+1, return to step (8).
The general principle and principal character and advantages of the present invention of the present invention has been shown and described above.The technology of the industry
Personnel are it should be appreciated that the present invention is not limited to the above embodiments, and the simply explanation described in above-described embodiment and specification is originally
The principle of invention, without departing from the spirit and scope of the present invention, various changes and modifications of the present invention are possible, these changes
Change and improvement all fall within the protetion scope of the claimed invention.The claimed scope of the invention by appended claims and its
Equivalent thereof.
Claims (4)
1. the photovoltaic grid inverter controller parameter identification method based on differential evolution algorithm, it is characterised in that：Including following
Step,
Step (1), population scale NP=50 is set, moment t1, t2 glade are respectively obtained according to formula (1) and formula (2)
BodyWith
Wherein,The proportionality coefficient k of d, q shaft current inner ring pi regulator in i-th of population of t1 moment is corresponded to respectivelyP, product
Divide coefficient kI, i=1,2 ..., NP；WithD, q shaft current inner ring in i-th of population of t2 moment are corresponded to respectively
Integrator output violent change value Lint_up、Lint_low, and current inner loop controller output violent change value Lout_up、Lout_low, i=1,
2 ..., NP；
T1=0, t2=0 are made, initial population is setIn section [kP_1/ 5,5kP_2] interior value,In area
Between [kI_1/ 5,5kI_2] interior value,Section (0,2.0] interior value,In section [- 2.0,0]
Value, according to formula (3), k is calculatedP_1、kP_2、kI_1And kI_2Value,
ωcFor cut-off angular frequency, value is 10 ωs, ωsFor synchronous rotary angular speed；ξ is damping ratio, ωnFor natural oscillation frequency
Rate, ξ=0.707, ωn=10 ωs；L1For the filter inductance of photovoltaic combining inverter AC, R1For filter inductance L1Internal resistance；
According to formula (4), correspondence function is listed
Wherein, idAnd i (k)q(k) it is the reality of photovoltaic combining inverter current inner loop controller in the case of current inner loop command value
Border exports；WithFor in the case of identical current inner loop command value, photovoltaic combining inverter current inner loop controller with
The output of track model；For solution vector；N=48, it is one
The sampling number of individual power frequency period；
Grid entry point voltage is set to drop to 0.85pu from 1.0pu so that current inner loop controller is operated in the range of linearity, and performs
Step (2)；
Step (2), if t1>NG, NG are maximum evolutionary generation, then the identification process of the proportionality coefficient of pi regulator and integral coefficient
Terminate, populationThe minimum solution vector of middle correspondence function value is the proportionality coefficient and integral coefficient of the pi regulator picked out,
And perform step (7)；Otherwise, step (3) is performed；
Step (3), utilizes correspondence functionTo populationIn each solution vector evaluated, if populationThe middle matching condition that solution vector satisfaction be present as shown in formula (5), then parameter identification success,
Wherein, ε values are 0.01；Corresponding solution vector is the proportionality coefficient and integral coefficient of the pi regulator picked out, and holds
Row step (7)；If populationIn each solution vector be unsatisfactory for matching condition as shown in formula (5), then perform step
(4)；
Step (4), to populationIn solution vectorMutation operation is carried out according to formula (6), obtains variation vector
Wherein, r1, r2, r3 ∈ { 1,2,3 ..., NP }, not equal to i, and are not mutually equal two-by-two；F is zoom factor, value model
Enclose for [0,1]；
Step (5), by populationIn solution vectorWith variation vectorThe crossover operation as shown in formula (7) is carried out,
Produce trial vector
Wherein, rand (j) ∈ [0,1], it is equally distributed random number；J ∈ { 1,2 }, represent j-th of variable；randn(i)∈
{ 1,2 }, indexed for randomly selected dimension variable；CR is crossover probability constant, and span is [0,1]；
Step (6), by making a variation and caused trial vector after crossover operationWith populationIn solution vectorCarry out
Competition as shown in formula (8), and select the less vector of correspondence function value to be used as filial generation, obtain the population at t1+1 moment
Make t1=t1+1, return to step (2)；
Step (7), grid entry point voltage is set to drop to 0.4pu from 1.0pu so that current inner loop controller is operated in inelastic region
Domain, and perform step (8)；
Step (8), if t2>NG, then amplitude limit link parameter identification process terminate, populationThe minimum solution of middle correspondence function value
Vector is the amplitude limit link parameter picked out；Otherwise, step (9) is performed；
Step (9), if populationMiddle have solution vector and meet matching condition as shown in formula (5), then amplitude limit link parameter is distinguished
Know successfully, identification process terminates, and corresponding solution vector is the amplitude limit link parameter picked out；If populationIn each solution to
Amount is unsatisfactory for the matching condition as shown in formula (5), then performs step (10)；
Step (10), to populationIn solution vectorMutation operation is carried out according to formula (9), obtains variation vector
Wherein, r1, r2, r3 ∈ { 1,2,3 ..., NP }, not equal to i, and are not mutually equal two-by-two；F is zoom factor, value model
Enclose for [0,1]；
Step (11), by populationIn solution vectorWith variation vectorCarry out the intersection as shown in formula (10)
Operation, produce trial vector
Wherein, rand (j) ∈ [0,1], it is equally distributed random number, j ∈ { 3,4,5,6 }, represents j-th of variable；randn
(i) ∈ { 3,4,5,6 }, indexed for randomly selected dimension variable；CR is crossover probability constant, and span is [0,1]；
Step (12), by making a variation and caused trial vector after crossover operationWith populationIn solution vector
The competition as shown in formula (11) is carried out, and selects the less vector of correspondence function value to be used as filial generation, obtains the t2+1 moment
Population
Make t2=t2+1, return to step (8).
2. the photovoltaic grid inverter controller parameter identification method according to claim 1 based on differential evolution algorithm,
It is characterized in that：Step (2), the maximum evolutionary generation NG of step (8) are 1000.
3. the photovoltaic grid inverter controller parameter identification method according to claim 1 based on differential evolution algorithm,
It is characterized in that：Step (4), step (10) zoom factor F preferred value are 0.5.
4. the photovoltaic grid inverter controller parameter identification method according to claim 1 based on differential evolution algorithm,
It is characterized in that：Step (5), step (11) crossover probability constant CR preferred value are 0.75.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510472410.0A CN105162159B (en) | 2015-08-05 | 2015-08-05 | Photovoltaic grid inverter controller parameter identification method based on differential evolution algorithm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510472410.0A CN105162159B (en) | 2015-08-05 | 2015-08-05 | Photovoltaic grid inverter controller parameter identification method based on differential evolution algorithm |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105162159A CN105162159A (en) | 2015-12-16 |
| CN105162159B true CN105162159B (en) | 2017-11-21 |
Family
ID=54802954
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510472410.0A Active CN105162159B (en) | 2015-08-05 | 2015-08-05 | Photovoltaic grid inverter controller parameter identification method based on differential evolution algorithm |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105162159B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106992551B (en) * | 2017-06-05 | 2019-06-18 | 合肥工业大学 | Photovoltaic inversion controller parameter discrimination method based on fuzzy C-mean algorithm and differential evolution hybrid algorithm |
| CN110504887B (en) * | 2019-08-29 | 2021-06-11 | 阳光电源股份有限公司 | Motor controller and control method thereof |
| CN113497457B (en) * | 2021-08-13 | 2022-06-10 | 国网天津市电力公司 | Photovoltaic power station grid-connected control method based on state parameter identification |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9331487B2 (en) * | 2013-03-14 | 2016-05-03 | Rockwell Automation Technologies, Inc. | Method and apparatus for islanding detection for grid tie converters |
| CN103166247B (en) * | 2013-03-15 | 2015-03-11 | 东方电气集团东方汽轮机有限公司 | System and method for controlling doubly-fed wind power generation grid-side converter |
| CN103208815B (en) * | 2013-04-02 | 2014-11-26 | 清华大学 | d-q axis parameter identification method for grid-connected inverter of photovoltaic power generation system |
| CN103956769B (en) * | 2014-04-08 | 2016-01-13 | 华中科技大学 | A kind of control method of the micro-grid connection inverter based on fuzzyPI |
| CN104319814A (en) * | 2014-11-10 | 2015-01-28 | 国家电网公司 | Method and system for determining control parameters of photovoltaic grid-connected inverter |
-
2015
- 2015-08-05 CN CN201510472410.0A patent/CN105162159B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105162159A (en) | 2015-12-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109672221B (en) | Direct-drive wind power plant dynamic equivalence method for subsynchronous oscillation analysis | |
| CN105703382A (en) | Subsynchronous oscillation suppression method and system for wind power plant | |
| CN105162159B (en) | Photovoltaic grid inverter controller parameter identification method based on differential evolution algorithm | |
| CN103078324B (en) | Based on the active filter current automatic adaptation control method of Fuzzy~+ PI | |
| CN101917150A (en) | Robust controller of permanent magnet synchronous motor based on fuzzy-neural network generalized inverse and construction method thereof | |
| CN109962495B (en) | Method for positioning and restraining ultralow-frequency oscillation disturbance source | |
| CN104836498A (en) | Auto-tuning PID generator control system based on artificial neural network | |
| CN106877363A (en) | A kind of SSSC suppression system sub-synchronous oscillation method and device | |
| CN110492479A (en) | A kind of rotary inertia of distributed grid-connected equipment and damping discrimination method | |
| CN103023418A (en) | Online parameter identification method of synchronous generator based on wide-area measurement information | |
| Aissaoui et al. | Application of fuzzy logic in control of electrical machines | |
| CN105337290A (en) | Reactive adjustment method applicable to low-frequency oscillation aid decision of electric system | |
| CN102790581A (en) | Constructing method for robust controller for radial position of bearingless asynchronous motor | |
| Ray et al. | Swarm and bacterial foraging based optimal power system stabilizer for stability improvement | |
| Singhal et al. | Advanced adaptive particle swarm optimization based SVC controller for power system stability | |
| Dehghani et al. | An improved PSO-GA based fuzzy controller of STATCOM for cascaded induction motors of drilling system | |
| CN103036252A (en) | Wind power plant model checking system and method based on trace sensitivity method | |
| CN111737919A (en) | Direct-drive wind power plant grouping method suitable for subsynchronous oscillation analysis | |
| CN103956767A (en) | Wake-effect-considering wind power plant grid-connection stability analysis method | |
| CN110838725A (en) | Parameter setting method and device for wind power plant closed-loop PI controller | |
| Peng et al. | Power system stabilizer parameters optimization using immune genetic algorithm | |
| Ma et al. | Research on fault location in DC distribution network based on adaptive artificial bee colony slime mould algorithm | |
| CN111525612B (en) | Real-time stable control device suitable for new forms of energy system that is incorporated into power networks | |
| Yin et al. | Reduced-dimensional skip-inception feature-aggregated classified proportional-integral-derivative for suppression of mixed-mode oscillations in hydropower units | |
| CN113721461A (en) | New energy unit parameter identification method and system based on multiple test scenes |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |