CN107846041B - Differential optimization control method for direct-drive permanent magnet synchronous wind power generation system - Google Patents

Differential optimization control method for direct-drive permanent magnet synchronous wind power generation system Download PDF

Info

Publication number
CN107846041B
CN107846041B CN201711122161.8A CN201711122161A CN107846041B CN 107846041 B CN107846041 B CN 107846041B CN 201711122161 A CN201711122161 A CN 201711122161A CN 107846041 B CN107846041 B CN 107846041B
Authority
CN
China
Prior art keywords
phase
permanent magnet
direct
power generation
grid
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
CN201711122161.8A
Other languages
Chinese (zh)
Other versions
CN107846041A (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.)
Wenzhou University
Original Assignee
Wenzhou University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wenzhou University filed Critical Wenzhou University
Priority to CN201711122161.8A priority Critical patent/CN107846041B/en
Publication of CN107846041A publication Critical patent/CN107846041A/en
Application granted granted Critical
Publication of CN107846041B publication Critical patent/CN107846041B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H02J3/386—
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/76—Power conversion electric or electronic aspects
    • 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
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
    • 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
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Economics (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Geometry (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

本发明公开一种直驱永磁同步风力发电系统的差分优化控制方法，采用背靠背全功率变流的电能控制结构，分别建立机侧变流器控制系统和网侧变流器控制系统的状态空间模型，并将直驱永磁同步风力发电系统中的机侧变流控制模块、网侧变流控制模块的跟踪误差的绝对值与时间的乘积和系统输出电压波形总谐波畸变率的加权叠加作为评估控制性能的适应度函数，设计基于实数编码的差分优化控制方法，实现直驱永磁同步风力发电系统中PI控制器参数的优化。本发明的控制方法能有效提高直驱永磁同步风力发电系统工作效率和发电质量，风电机组的能量转换效率、运行可靠性，并在电网发生扰动时确保直驱永磁同步风力发电系统稳定时间更短，稳态误差更小，鲁棒性更强。

Figure 201711122161

The invention discloses a differential optimization control method of a direct-drive permanent magnet synchronous wind power generation system, which adopts a back-to-back full-power conversion electric energy control structure to establish the state spaces of a machine-side converter control system and a grid-side converter control system respectively. model, and the product of the absolute value of the tracking error and the time of the machine-side converter control module and the grid-side converter control module in the direct-drive permanent magnet synchronous wind power generation system and the weighted superposition of the total harmonic distortion rate of the system output voltage waveform As a fitness function to evaluate the control performance, a differential optimization control method based on real number coding is designed to optimize the parameters of the PI controller in the direct-drive permanent magnet synchronous wind power generation system. The control method of the invention can effectively improve the working efficiency and power generation quality of the direct-drive permanent magnet synchronous wind power generation system, the energy conversion efficiency and operation reliability of the wind turbine, and ensure the stability time of the direct-drive permanent magnet synchronous wind power generation system when the power grid is disturbed Shorter, smaller steady-state error, more robust.

Figure 201711122161

Description

Differential optimization control method for direct-drive permanent magnet synchronous wind power generation system
Technical Field
The invention relates to the field of intelligent control technology of new energy power generation systems, in particular to a differential optimization control method of a direct-drive permanent magnet synchronous wind power generation system.
Background
A direct-drive Permanent Magnet synchronous generator (D-PMSG) omits a complex mechanical transmission mechanism, reduces the loss of a wind generating set by adopting full-power variable-current electric energy control, and improves the energy conversion efficiency and the operation reliability of the wind generating set. The power control strategy of the full-power converter has very critical influence on the operating characteristics and the power output quality of the D-PMSG. On the background that the capacity of a wind power access power grid is rapidly increased, a grid-connected direct-drive wind turbine generator set has strong fault resistance, and has strong robustness and certain anti-interference capability when the power grid is disturbed. Therefore, how to design an effective full-power grid-connected converter control system of the D-PMSG has important engineering application value.
The full-power grid-connected current conversion control strategy of the D-PMSG is a common organic side uncontrollable rectifier grid-connected side thyristor inversion control strategy, a machine side uncontrollable rectifier grid-connected side PWM voltage source type inversion control strategy, a machine side uncontrollable rectifier grid-connected Boost voltage source inversion control strategy, a back-to-back full-power current conversion control strategy and the like. The back-to-back full-power variable flow control strategy regulates and controls the rotating speed and the reactive power of the permanent magnet synchronous engine by adopting a machine side PWM rectification system to realize the tracking of the optimal power of wind energy; the intermediate direct-current voltage is stabilized by adopting a grid-side PWM inversion system, and the reactive power injected into the power grid is regulated and controlled. The back-to-back full-power variable flow control strategy is more flexible in control because both ends are PWM variable flow systems, but complexity and uncertainty of system optimization control are increased. The differential optimization control is used as a novel optimization control algorithm based on a cluster intelligent theory, the application research in the power system is late, and the method mainly focuses on the aspects of power grid planning, load economic distribution, optimal load flow calculation and the like. The research on how to design the DE optimization control method of the D-PMSG back-to-back full-power converter electric energy control system by setting a multi-objective weighting function so as to improve the working efficiency and the power generation quality of the wind power system under complex working conditions and uncertain factors is a well-recognized problem in the academic and scientific fields and the engineering application fields at home and abroad.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a differential optimization control method of a direct-drive permanent magnet synchronous wind power generation system, which can effectively improve the working efficiency and the power generation quality of a D-PMSG (direct-drive permanent magnet synchronous wind power generation system) and the energy conversion efficiency and the operation reliability of a wind turbine generator set, and ensure that the D-PMSG has shorter stabilization time, smaller steady-state error and stronger robustness when a power grid is disturbed. The specific technical scheme is as follows:
a differential optimization control method of a direct-drive permanent magnet synchronous wind power generation system is characterized in that the system adopts a back-to-back full-power converter electric energy control structure, wherein a machine side converter controls the rotating speed and the output power of a generator, a grid side converter stabilizes the voltage of a direct current bus and controls the active power and the reactive power output by the power generation system; the method comprises the following steps:
(1) establishing a state space model of the grid-side converter system by adopting a modeling method described by a duty ratio function:
Figure BDA0001467595750000021
wherein, L is filter inductance, R is equivalent resistance of the filter inductance, C filter capacitance, ekIs the three-phase grid potential, ikFor three-phase output current, dkIs the duty ratio, udcIs a DC bus voltage, edIs a direct current electromotive force at the side of the direct current bus, RdThe equivalent resistance is the equivalent resistance of the direct current bus side;
(2) d and q synchronous rotating coordinate transformation is carried out on the formula (1), so that the state space model of the grid-side converter system is adjusted as follows:
Figure BDA0001467595750000022
wherein e isd、eqIs d-axis component and q-axis component i under two-phase synchronous rotating coordinate system of three-phase power grid electromotive forced、iqD and q axis components under a two-phase synchronous rotating coordinate system of output current of a grid-side converter systemd、dqD and q axis components under an equivalent switching function two-phase synchronous rotating coordinate system; omega is the three-phase grid voltage angular frequency.
(3) By modeling the machine side converter system under the two-phase synchronous rotating coordinate system, a voltage equation of the three-phase stator winding can be obtained:
Figure BDA0001467595750000023
wherein u issa、usb、usc-stator three-phase voltage，isa、isb、isc-stator three phase currents. RsWinding resistance of stator,. psia、ψb、ψc-stator winding flux per phase expressed as follows:
Figure BDA0001467595750000024
wherein L isaa、Lbb、LccFor the stator winding inductance per phase, Mab＝Mba、Mbc＝Mcb、Mac＝McaFor mutual inductance, ψ, between windings of each phasefa、ψfb、ψfcPermanent magnet flux linkage of each pole, represented as a three-phase winding:
Figure BDA0001467595750000031
wherein psifExciting a flux linkage for a permanent magnet;
(4) the state space equation of the machine side converter system is obtained according to the formulas (3), (4) and (5) as follows:
Figure BDA0001467595750000032
where p is a differential operator, LsIs the winding inductance of the stator. In the case of a three-phase system,
isa+isb+isc＝0 (7)
combining the equations (4) - (7), and further transforming the machine side converter system state space equation into:
Figure BDA0001467595750000033
wherein, ω iseIs the electrical angular velocity of the rotor;
(5) d and q synchronous rotating coordinate transformation is carried out on the formula (8), and finally, the state space model of the machine side converter system is adjusted as follows:
Figure BDA0001467595750000034
Figure BDA0001467595750000035
wherein u issd、usqIs d and q axis components i of a stator three-phase voltage two-phase synchronous rotating coordinate systemsd、isqD and q axis components psi under stator three-phase current two-phase synchronous rotating coordinate systemd、ψqD and q axis components of generator flux linkage, Ld、LqD, q-axis component, omega, of stator winding inductancesIs the generator angular velocity.
(6) Combining formulas (2), (9) and (10) with a PI controller to obtain an electric energy control system of a back-to-back full-power converter in the direct-drive permanent magnet synchronous wind power generation system;
(7) and (5) carrying out differential optimization on the system obtained in the step (6), and setting the parameter values of the differential optimization: a variation factor F, a cross factor CR, a population scale M and a maximum iteration number G;
(8) randomly generating an initial population P ═ x of the variable of a PI controller of an electric energy control system for back-to-back full-power conversion in the real number coded direct-drive permanent magnet synchronous wind power generation meeting the constraint condition of a formula (11)1,x2,…,xpWherein the ith individual xiRepresenting a sequence of control deltas K to be optimizedp1i,Ki1i,Kp2i,Ki2i,Kp3i,Ki3i,Kp4i,Ki4iThe specific production process is as follows:
xij＝Δumin+randij(Δumax-Δumin),i＝1,2,...,p；j＝1,2,...,8 (11)
wherein, Δ uminAnd Δ umaxRespectively representing the lower and upper limits, rand, of the sequence of control incrementsijRepresents a set of random numbers generated between 0 and 1;
(9) for each individual x in the population P according to equation (11)i,i＝1,2,…,pPerforming an objective function fiAnd (4) calculating and evaluating, specifically as shown in formula (13), and setting the current minimum objective function value in the population as fbestSetting the corresponding individual as the current best solution Sbest；
Figure BDA0001467595750000041
Wherein e is1And e2Respectively representing the tracking errors of the machine side variable flow control module and the network side variable flow control module, T representing the system running time value, TminAnd TmaxRespectively representing the initial time and the end time of the operation of the three-phase inverter system, THD representing the harmonic distortion rate of the voltage waveform output by the direct-drive permanent magnet synchronous wind power generation system, w1And w2Represents a weight coefficient and satisfies w1+w2＝1；
(10) Performing mutation operation on the population, and randomly selecting 3 individuals x from the populationp1，xp2，xp3And i, p1, p2 and p3 are different from each other, and the specific mutation operations are as follows:
hij(t+1)＝xp1j(t)+F·(xp2j(t)-xp3j(t)) (13)
wherein h isij(t +1) temporarily storing the intermediate variable of the variant individual;
(11) in order to increase the diversity of interference parameter vectors, the population after mutation operation is subjected to cross operation, which specifically comprises the following steps:
Figure BDA0001467595750000042
wherein, randlijRepresenting a set of random numbers, v, generated between 0 and 1ij(t +1) is an intermediate variable of the individual after temporary storage of the crossover;
(12) to determine xi(t) whether it can be a member of the next generation, and obtaining an experimental vector v using the formula (12)i(t +1) and target vector xi(t) and comparing the values of the evaluation functions as follows:
Figure BDA0001467595750000051
(13) and (5) repeating the steps (9) to (12) until the system runs to the maximum iteration number G to obtain the optimal variable of the PI controller, so that the back-to-back full-power current transformation electric energy control system in the direct-drive permanent magnet synchronous wind power generation is controlled.
Compared with the prior art, the invention has the following beneficial effects:
by adopting the differential optimization control method of the D-PMSG back-to-back full-power conversion system, the working efficiency and the power generation quality of the D-PMSG, the energy conversion efficiency and the operation reliability of the wind turbine generator set can be effectively improved, and the D-PMSG is ensured to have shorter stabilization time, smaller steady-state error and stronger robustness when the power grid is disturbed.
Drawings
FIG. 1 is a structural block diagram of a direct-drive permanent magnet synchronous wind power generation system and a schematic diagram of DE optimization control principle;
FIG. 2 is a flow chart of an implementation process of the control method of the present invention;
FIG. 3 is a graph of the objective function value optimization process after the control method of the present invention is implemented;
FIG. 4 is a three-phase power output voltage waveform diagram of the D-PMSG implemented in the RT-LAB power real-time simulation platform.
Detailed Description
The following representative examples are intended to illustrate the invention, but are not intended to limit the scope of the invention described herein.
As shown in fig. 1 and 2, a differential optimization control method for a direct-drive permanent magnet synchronous wind power generation system adopts a back-to-back full-power converter electric energy control structure, wherein a machine-side converter controls the rotation speed and output power of a generator, a grid-side converter stabilizes the voltage of a direct-current bus, and controls the active power and reactive power output by the power generation system; the method comprises the following steps:
a differential optimization control method of a direct-drive permanent magnet synchronous wind power generation system is characterized in that the system adopts a back-to-back full-power converter electric energy control structure, wherein a machine side converter controls the rotating speed and the output power of a generator, a grid side converter stabilizes the voltage of a direct current bus and controls the active power and the reactive power output by the power generation system; the method comprises the following steps:
(1) establishing a state space model of the grid-side converter system by adopting a modeling method described by a duty ratio function:
Figure BDA0001467595750000061
wherein, L is filter inductance, R is equivalent resistance of the filter inductance, C filter capacitance, ekIs the three-phase grid potential, ikFor three-phase output current, dkIs the duty ratio, udcIs a DC bus voltage, edIs a direct current electromotive force at the side of the direct current bus, RdThe equivalent resistance is the equivalent resistance of the direct current bus side;
(2) d and q synchronous rotating coordinate transformation is carried out on the formula (1), so that the state space model of the grid-side converter system is adjusted as follows:
Figure BDA0001467595750000062
wherein e isd、eqIs d-axis component and q-axis component i under two-phase synchronous rotating coordinate system of three-phase power grid electromotive forced、iqD and q axis components under a two-phase synchronous rotating coordinate system of output current of a grid-side converter systemd、dqD and q axis components under an equivalent switching function two-phase synchronous rotating coordinate system; omega is the three-phase grid voltage angular frequency.
(3) By modeling the machine side converter system under the two-phase synchronous rotating coordinate system, a voltage equation of the three-phase stator winding can be obtained:
Figure BDA0001467595750000063
wherein u issa、usb、uscStator three-phase voltage isa、isb、isc-stator three phase currents. RsWinding resistance of stator,. psia、ψb、ψc-stator winding flux per phase expressed as follows:
Figure BDA0001467595750000064
wherein L isaa、Lbb、LccFor the stator winding inductance per phase, Mab＝Mba、Mbc＝Mcb、Mac＝McaFor mutual inductance, ψ, between windings of each phasefa、ψfb、ψfcPermanent magnet flux linkage of each pole, represented as a three-phase winding:
Figure BDA0001467595750000071
wherein psifExciting a flux linkage for a permanent magnet;
(4) the state space equation of the machine side converter system is obtained according to the formulas (3), (4) and (5) as follows:
Figure BDA0001467595750000072
where p is a differential operator, LsIs the winding inductance of the stator. In the case of a three-phase system,
isa+isb+isc＝0 (7)
combining the equations (4) - (7), and further transforming the machine side converter system state space equation into:
Figure BDA0001467595750000073
wherein, ω iseIs the electrical angular velocity of the rotor;
(5) d and q synchronous rotating coordinate transformation is carried out on the formula (8), and finally, the state space model of the machine side converter system is adjusted as follows:
Figure BDA0001467595750000074
Figure BDA0001467595750000075
wherein u issd、usqIs d and q axis components i of a stator three-phase voltage two-phase synchronous rotating coordinate systemsd、isqD and q axis components psi under stator three-phase current two-phase synchronous rotating coordinate systemd、ψqD and q axis components of generator flux linkage, Ld、LqD, q-axis component, omega, of stator winding inductancesIs the generator angular velocity.
(6) Combining formulas (2), (9) and (10) with a PI controller to obtain an electric energy control system of a back-to-back full-power converter in the direct-drive permanent magnet synchronous wind power generation system;
(7) and (5) carrying out differential optimization on the system obtained in the step (6), and setting the parameter values of the differential optimization: a variation factor F, a cross factor CR, a population scale M and a maximum iteration number G;
(8) randomly generating an initial population P ═ x of the variable of a PI controller of an electric energy control system for back-to-back full-power conversion in the real number coded direct-drive permanent magnet synchronous wind power generation meeting the constraint condition of a formula (11)1,x2,…,xpWherein the ith individual xiRepresenting a sequence of control deltas K to be optimizedp1i,Ki1i,Kp2i,Ki2i,Kp3i,Ki3i,Kp4i,Ki4iThe specific production process is as follows:
xij＝Δumin+randij(Δumax-Δumin),i＝1,2,...,p；j＝1,2,...,8 (11)
wherein, Δ uminAnd Δ umaxRespectively representing the lower and upper limits, rand, of the sequence of control incrementsijRepresents a set of random numbers generated between 0 and 1;
(9) for each individual in the population P according to equation (11)xiI 1,2, …, p performs the objective function fiAnd (4) calculating and evaluating, specifically as shown in formula (13), and setting the current minimum objective function value in the population as fbestSetting the corresponding individual as the current best solution Sbest；
Figure BDA0001467595750000081
Wherein e is1And e2Respectively representing the tracking errors of the machine side variable flow control module and the network side variable flow control module, T representing the system running time value, TminAnd TmaxRespectively representing the initial time and the end time of the operation of the three-phase inverter system, THD representing the harmonic distortion rate of the voltage waveform output by the direct-drive permanent magnet synchronous wind power generation system, w1And w2Represents a weight coefficient and satisfies w1+w2＝1；
(10) Performing mutation operation on the population, and randomly selecting 3 individuals x from the populationp1，xp2，xp3And i, p1, p2 and p3 are different from each other, and the specific mutation operations are as follows:
hij(t+1)＝xp1j(t)+F·(xp2j(t)-xp3j(t)) (13)
wherein h isij(t +1) temporarily storing the intermediate variable of the variant individual;
(11) in order to increase the diversity of interference parameter vectors, the population after mutation operation is subjected to cross operation, which specifically comprises the following steps:
Figure BDA0001467595750000082
wherein, randlijRepresenting a set of random numbers, v, generated between 0 and 1ij(t +1) is an intermediate variable of the individual after temporary storage of the crossover;
(12) to determine xi(t) whether it can be a member of the next generation, and obtaining an experimental vector v using the formula (12)i(t +1) and target vector xi(t) evaluating the values of the functions, comparing them, and performing the specific operationThe following were used:
Figure BDA0001467595750000091
(13) and (5) repeating the steps (9) to (12) until the system runs to the maximum iteration number G to obtain the optimal variable of the PI controller, so that the back-to-back full-power current transformation electric energy control system in the direct-drive permanent magnet synchronous wind power generation is controlled.
The optimization process curve of the objective function value after the control method is implemented is shown in fig. 3, the three-phase power output voltage waveform of the real-time simulation platform of the RT-LAB power of the D-PMSG after the control method is implemented is shown in fig. 4, and as can be seen from the graph, the differential optimization method has strong global convergence capability and robustness, and after the differential optimization calculation is adopted to carry out optimization control on the electric energy control system of the back-to-back full-power converter in the D-PMSG, the three-phase power output of the system is stable, and the electric energy quality is high.

Claims (1)

1. A differential optimization control method of a direct-drive permanent magnet synchronous wind power generation system is characterized in that the system adopts a back-to-back full-power converter electric energy control structure, wherein a machine side converter controls the rotating speed and the output power of a generator, a grid side converter stabilizes the voltage of a direct current bus and controls the active power and the reactive power output by the power generation system; the method comprises the following steps:
(1) establishing a state space model of the grid-side converter system by adopting a modeling method described by a duty ratio function:
Figure FDA0002386290260000011
wherein, L is filter inductance, R is equivalent resistance of filter inductance, C filter capacitance, epsilonkIs the three-phase grid potential, ikFor three-phase output current, dkIs the duty ratio, udcIs a DC bus voltage, εdIs a direct current electromotive force at the side of the direct current bus, RdThe equivalent resistance is the equivalent resistance of the direct current bus side;
(2) d and q synchronous rotating coordinate transformation is carried out on the formula (1), so that the state space model of the grid-side converter system is adjusted as follows:
Figure FDA0002386290260000012
wherein e isd、eqIs d-axis component and q-axis component i under two-phase synchronous rotating coordinate system of three-phase power grid electromotive forced、iqD and q axis components under a two-phase synchronous rotating coordinate system of output current of a grid-side converter systemd、dqD and q axis components under an equivalent switching function two-phase synchronous rotating coordinate system; omega is the angular frequency of the three-phase power grid voltage;
(3) by modeling the machine side converter system under the two-phase synchronous rotating coordinate system, a voltage equation of the three-phase stator winding can be obtained:
Figure FDA0002386290260000013
wherein u issa、usb、uscStator three-phase voltage isa、isb、isc-stator three phase currents; rsWinding resistance of stator,. psia、ψb、ψc-stator winding flux per phase expressed as follows:
Figure FDA0002386290260000021
wherein L isaa、Lbb、LccFor the stator winding inductance per phase, Mab＝Mba、Mbc＝Mcb、Mac＝McaFor mutual inductance, ψ, between windings of each phasefa、ψfb、ψfcPermanent magnet flux linkage of each pole, represented as a three-phase winding:
Figure FDA0002386290260000022
wherein psifExciting a flux linkage for a permanent magnet;
(4) the state space equation of the machine side converter system is obtained according to the formulas (3), (4) and (5) as follows:
Figure FDA0002386290260000023
wherein D is a differential operator, LsIs the winding inductance of the stator; in the case of a three-phase system,
isa+isb+isc＝0 (7)
combining the equations (4) - (7), and further transforming the machine side converter system state space equation into:
Figure FDA0002386290260000024
wherein, ω iseIs the electrical angular velocity of the rotor;
(5) d and q synchronous rotating coordinate transformation is carried out on the formula (8), and finally, the state space model of the machine side converter system is adjusted as follows:
Figure FDA0002386290260000025
Figure FDA0002386290260000031
wherein u issd、usqIs d and q axis components i of a stator three-phase voltage two-phase synchronous rotating coordinate systemsd、isqIs d-axis component and q-axis component psi of stator three-phase current two-phase synchronous rotating coordinate systemd、ψqIs the d and q axis components of the generator flux linkage, Ld、LqD, q-axis components, omega, of stator winding inductancesIs the generator angular velocity;
(6) combining formulas (2), (9) and (10) with a PI controller to obtain an electric energy control system of a back-to-back full-power converter in the direct-drive permanent magnet synchronous wind power generation system;
(7) and (5) carrying out differential optimization on the system obtained in the step (6), and setting the parameter values of the differential optimization: a variation factor F, a cross factor CR, a population scale M and a maximum iteration number G;
(8) randomly generating an initial population P ═ x of the variable of the PI controller of the electric energy control system of the back-to-back full-power conversion in the real number coded direct-drive permanent magnet synchronous wind power generation system which meets the constraint condition of the formula (11)1,x2,…,xpWherein the ith individual xiRepresenting a sequence of control deltas K to be optimizedp1i,Ki1i,Kp2i,Ki2i,Kp3i,Ki3i,Kp4i,Ki4iThe specific production process is as follows:
xij＝Δumin+randij(Δumax-Δumin),i＝1,2,...,p；j＝1,2,...,8 (11)
wherein, Δ uminAnd Δ umaxRespectively representing the lower and upper limits, rand, of the sequence of control incrementsijRepresents a set of random numbers generated between 0 and 1;
(9) for each individual x in the population P according to equation (12)iI 1,2, …, p performs the objective function fiCalculating and evaluating, and setting the current minimum objective function value in the population as fbestSetting the corresponding individual as the current best solution Sbest；
Figure FDA0002386290260000032
Wherein e is1And e2Respectively representing the tracking errors of a machine side variable flow control module and a network side variable flow control module, tau representing the running time value of the system, TminAnd TmaxRespectively representing the initial time and the termination time of the operation of a grid-side converter system, THD representing the harmonic distortion rate of the voltage waveform output by the direct-drive permanent magnet synchronous wind power generation system, w1And w2Represents a weight coefficient and satisfies w1+w2＝1；
(10) Performing mutation operation on the population, and randomly selecting 3 individuals x from the populationp1，xp2，xp3And i, p1, p2 and p3 are different from each other, and the specific mutation operations are as follows:
hij(t+1)＝xp1j(t)+F·(xp2j(t)-xp3j(t)) (13)
wherein h isij(t +1) temporarily storing the intermediate variable of the variant individual;
(11) in order to increase the diversity of interference parameter vectors, the population after mutation operation is subjected to cross operation, which specifically comprises the following steps:
Figure FDA0002386290260000041
wherein, randlijRepresenting a set of random numbers, v, generated between 0 and 1ij(t +1) is an intermediate variable of the individual after temporary storage of the crossover;
(12) to determine xi(t) whether it can be a member of the next generation, and obtaining an experimental vector v using the formula (12)i(t +1) and target vector xi(t) and comparing the values of the evaluation functions as follows:
Figure FDA0002386290260000042
(13) and (5) repeating the steps (9) to (12) until the system runs to the maximum iteration number G to obtain the optimal variable of the PI controller, so that the back-to-back full-power current transformation electric energy control system in the direct-drive permanent magnet synchronous wind power generation system is controlled.
CN201711122161.8A 2017-11-14 2017-11-14 Differential optimization control method for direct-drive permanent magnet synchronous wind power generation system Active CN107846041B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711122161.8A CN107846041B (en) 2017-11-14 2017-11-14 Differential optimization control method for direct-drive permanent magnet synchronous wind power generation system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711122161.8A CN107846041B (en) 2017-11-14 2017-11-14 Differential optimization control method for direct-drive permanent magnet synchronous wind power generation system

Publications (2)

Publication Number Publication Date
CN107846041A CN107846041A (en) 2018-03-27
CN107846041B true CN107846041B (en) 2020-04-24

Family

ID=61678831

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711122161.8A Active CN107846041B (en) 2017-11-14 2017-11-14 Differential optimization control method for direct-drive permanent magnet synchronous wind power generation system

Country Status (1)

Country Link
CN (1) CN107846041B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109193684B (en) * 2018-08-14 2021-09-07 河海大学 A real-time reactive power optimization method for power system based on two-stage optimization
CN109347382B (en) * 2018-11-24 2022-05-13 沈阳工业大学 A Rotor Position Estimation Method for Permanent Magnet Direct Drive Wind Turbine
CN113258843B (en) * 2021-06-11 2022-12-13 盛东如东海上风力发电有限责任公司 Direct-drive wind turbine generator motor rotating speed control method, control system and grid-connected system
CN113422385A (en) * 2021-07-06 2021-09-21 上海工程技术大学 Grid-connected wind power generation system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007143313A (en) * 2005-11-18 2007-06-07 Hokuriku Electric Power Co Inc:The Method for controlling power distribution line voltage
CN103023021A (en) * 2012-11-27 2013-04-03 上海电气集团股份有限公司 Decoupling control method for nonlinear power of double-fed wind power generation system
CN103138672A (en) * 2013-03-13 2013-06-05 华北电力大学（保定） Active disturbance rejection control method of direct-driven permanent magnet synchronization wind power system
CN105700353A (en) * 2016-01-30 2016-06-22 河南城建学院 A PID controller parameter optimal setting method based on a differential evolution method
CN105867126A (en) * 2016-04-12 2016-08-17 温州大学 Fractional order PI optimization control method of three-phase voltage source type inverter system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007143313A (en) * 2005-11-18 2007-06-07 Hokuriku Electric Power Co Inc:The Method for controlling power distribution line voltage
CN103023021A (en) * 2012-11-27 2013-04-03 上海电气集团股份有限公司 Decoupling control method for nonlinear power of double-fed wind power generation system
CN103138672A (en) * 2013-03-13 2013-06-05 华北电力大学（保定） Active disturbance rejection control method of direct-driven permanent magnet synchronization wind power system
CN105700353A (en) * 2016-01-30 2016-06-22 河南城建学院 A PID controller parameter optimal setting method based on a differential evolution method
CN105867126A (en) * 2016-04-12 2016-08-17 温州大学 Fractional order PI optimization control method of three-phase voltage source type inverter system

Also Published As

Publication number Publication date
CN107846041A (en) 2018-03-27

Similar Documents

Publication Publication Date Title
CN107846041B (en) Differential optimization control method for direct-drive permanent magnet synchronous wind power generation system
CN105429183A (en) Permanent magnetic direct-drive type offshore wind power plant grid-connected system topology structure and control method thereof
CN106992514B (en) A small-disturbance stability analysis method for wind-storage isolated grid system
CN106602916A (en) Hybrid level three-phase four-bridge arm converter device and control method
CN104579060A (en) Indirect power control methods of cage-type rotor brushless double-fed wind generator
CN112117947A (en) SRM torque ripple suppression control system and method based on current injection method
Hamitouche et al. A new control strategy of dual stator induction generator with power regulation
CN102420566B (en) Control device of brushless doubly-fed machine and control method thereof
Mishra et al. Six-phase DFIG-MPPT synergy: pioneering approaches for maximising energy yield in wind energy generation system
CN113659905B (en) Three-level power generation system model prediction control method based on time-varying disturbance compensation
Vattuone et al. Open-end-winding PMSG for wind energy conversion system with dual boost NPC converter
CN112039393A (en) Control method with multiphase permanent magnet synchronous generator system and MMC
Jbarah et al. Improved dfig dftc by using a fractional-order super twisting algorithms in wind power application
Hwang et al. Vector control of multiple-module transverse flux PM generator for large-scale direct-drive wind turbines
Camurça et al. High efficiency wind energy conversion system based on the Three-Level Delta-Switch T-Type Converter and PMSG Model-Based loss minimization
CN114553014A (en) Pulse width modulation method for synthesizing multiphase full-parameter controllable alternating voltage by three-phase alternating current
Li et al. Sequential predictive control with dynamic priority of three-level NPC back-to-back power converters in PMSG wind turbine systems
Aguilar et al. Adaptive controller for PMSG wind turbine systems with back-to-back PWM converters
Naguru et al. Comparative study of power control of DFIG using PI control and FeedBack Linearization Control
Hamidia et al. Fuzzy direct torque controlled permanent magnet synchronous, doubly fed and induction motors
Sekhar et al. An Interconnected Wind Driven SEIG System Using SVPWM Controlled TL Z-Source Inverter Strategy for Off-Shore WECS
Zhang et al. Direct power and torque control for three-level NPC based PWM ac/dc/ac converter
Krishna et al. Direct predictive current control of grid connected neutral point clamped inverter for wave power extraction
ALMAKKI et al. MODERN TRANSPORTATION SYSTEMS AND TECHNOLOGIES
Al-Gabalawy et al. Controlling PMSG Wind Turbine Connected to Three Level Neutral Point Clamped Inverter Based on MPC

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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20180327

Assignee: Zhejiang Jiuhong Power Engineering Co.,Ltd.

Assignor: Wenzhou University

Contract record no.: X2021330000833

Denomination of invention: A differential optimal control method for direct drive permanent magnet synchronous wind power generation system

Granted publication date: 20200424

License type: Common License

Record date: 20211222

EE01 Entry into force of recordation of patent licensing contract