CN112350364B - 基于双同步旋转坐标系的全功率风力发电系统的建模方法 - Google Patents

基于双同步旋转坐标系的全功率风力发电系统的建模方法 Download PDF

Info

Publication number
CN112350364B
CN112350364B CN202011056354.XA CN202011056354A CN112350364B CN 112350364 B CN112350364 B CN 112350364B CN 202011056354 A CN202011056354 A CN 202011056354A CN 112350364 B CN112350364 B CN 112350364B
Authority
CN
China
Prior art keywords
current
module
axis component
sequence
bridge arm
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
CN202011056354.XA
Other languages
English (en)
Other versions
CN112350364A (zh
Inventor
刘芳
刘世权
李研
徐韫钰
王佳宁
刘威
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN202011056354.XA priority Critical patent/CN112350364B/zh
Publication of CN112350364A publication Critical patent/CN112350364A/zh
Application granted granted Critical
Publication of CN112350364B publication Critical patent/CN112350364B/zh
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/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • 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
    • H02J2300/28The renewable source being wind energy
    • 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
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects

Landscapes

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

Abstract

本发明公开了一种基于基于双同步旋转坐标系的全功率风力发电系统的建模方法,属于电力电子稳定性控制领域。尤其适用于变流器模组输入输出并联的直驱型全功率风力发电系统在双dq坐标系下建模方法。该建模方法包括采样,并建立了双dq坐标系下的7个模型。即通过改变传统的单dq建模方法,采用双dq控制方法建立直驱型全功率风力发电系统等效数学模型。而且建模过程是基于变流器模组并联得到的,不仅符合当下大功率风电变流器提高功率容量的方式,而且也能适用单模组变流器的情形,建模过程具有极强的适应性,更加贴近实际,大大简化了运算过程和系统维数,具有创新意义。

Description

基于双同步旋转坐标系的全功率风力发电系统的建模方法
技术领域
本发明涉及一种全功率风力发电系统建模方法,尤其是一种基于变流器模组输入输出并联的直驱型全功率风力发电系统在双dq坐标系下建模方法。
背景技术
近年来,随着能源的短缺与环境的恶化,风电、光伏等新能源得以迅猛发展,新能源在电网中的比例日益上升,逐渐形成了高比例新能源并网发电的电力系统格局。目前,作为变速恒频风力发电系统主流机型的直驱式全功率型风力发电系统,相对于双馈机组,永磁直驱式风力发电系统定子侧直接与背靠背PWM变流器相连,实现了发电机与电网的隔离,提高了直驱式永磁同步发电系统的低电压穿越能力,同时没有齿轮箱、电刷和滑环,增加了系统的可靠性。同时它采用背靠背PWM变流器,不仅易于实现变速恒频运行,并且易于实现功率解耦控制,使其成为当下研究的热点。随着风电并网在电网中的占比增多,特别是大区互联电网的出现和逐渐扩大,小干扰稳定性日益突出。这在高比例新能源并网的电力系统中会弱化系统的控制性能,给系统的安全稳定运行带来巨大压力和挑战。为此,很多专家学者提出了基于特征根的小信号分析方法来分析风力发电系统的小扰动稳定性问题。其中建立整个系统的小信号模型并进行仿真分析是重要的研究手段之一。
题为《基于小信号模型的永磁直驱风电系统特征值灵敏度分析》(向如意,杨金明,王林,余情,杨小斌.基于小信号模型的永磁直驱风电系统特征值灵敏度分析[J].电网与清洁能源,2015,31(05):94-101.)的文章建立了控制器理想输出时和基于背靠背PWM变流器的双闭环控制下的永磁直驱风电系统的小信号模型,并进行特征值灵敏度分析,为控制器的优化配置提供了理论支持。但并未指出控制器参数详细的优化方向。
题为《永磁风力发电机组小扰动稳定性分析》(冯闯.永磁风力发电机组小扰动稳定性分析[D].湖南大学,2013.)的文章推导了永磁风力发电机组完整的小信号模型,通过特征根分析法分析了永磁风力发电机组的振荡模式及其产生的机理,并给出了相应的解决方案提高永磁同步风力发电机组的阻尼,以此提高系统的稳定性。但建模过程中并未考虑锁相环的影响。
题为《直驱永磁型风电系统的小扰动稳定性分析》(李军军,吴政球,谭勋琼.直驱永磁型风电系统的小扰动稳定性分析[J].电力系统及其自动化学报,2012,24(05):49-57.)的文章建立了基于双PWM变换器的直驱永磁型风电系统的小扰动稳定性分析数学模型,研究了风电系统并网后不同参数变化对电力系统稳定性的影响。得出了若要保证系统具有较好的小扰动稳定性,线路电感不宜过大;直流电容要取值合适;母线电压可适当提高;增加浆距角调节机构,以保证系统在额定风速以上能安全地运行,为直驱永磁型风电并网安全稳定运行提供了可借鉴的理论依据。但是不仅在建模过程中没有考虑锁相环对系统的影响,而且其提高系统稳定性的方法主要是更改主电路的参数,这使得系统成本大幅增加,不利于大规模推广,具有一定的局限性。
题为《基于小扰动稳定分析的并网直驱永磁风电系统建模优化和仿真》(杨昆,杨秀.基于小扰动稳定分析的并网直驱永磁风电系统建模优化和仿真[J].华东电力,2014,42(10):2028-2033.)的文章建立了完整的风力发电系统小扰动模型,包括风轮机、永磁同步发电机、AC/DC控制器、L型滤波器、DC/AC控制器和直流侧电容,用特征值分析法分析研究了永磁直驱风力发电机接入无穷大电网的稳定性问题。但是仅在强电网下的一个稳态工作点进行了小干扰稳定性研究,且并未给出影响稳定性的因素。
综上所述,目前文献大多对直驱型全功率风力发电系统进行了各种不同近似程度的建模,但大多数都是采用单dq控制的建模方法,而实际中为应对端口电压的不平衡,常采用双dq的控制方法。同时常采用锁相环来跟踪电网频率,使变流器输出电压、电流频率与电网频率一致,而有的文献在建模过程中也没有考虑锁相环的影响。因此,理论分析与实际存在偏差,不能更好的反映实际系统小干扰稳定运行的情况,即直驱型全功率风力发电系统的双dq建模研究还有待于进一步深入和探索。
发明内容
本发明要解决的技术问题为克服上述各种技术方案的局限性。本发明针对一种直驱型全功率风力发电系统,尤其是一种基于模组输入输出并联、且采用永磁同步发电机的全功率型风力发电系统在双dq坐标系下的建模方法。
本发明的目的是这样实现的。本发明提供了一种基于双同步旋转坐标系的全功率风力发电系统的建模方法,所述全功率风力发电系统为变流器模组输入输出并联的直驱型全功率风力发电系统;所述变流器模组输入输出并联的直驱型全功率风力发电系统主电路的拓扑结构包括永磁同步电机M、机侧滤波电容C1、直流侧电容Cdc、机侧变流器K1、网侧变流器K2、网侧滤波电容C、网侧电感Lg和电网;所述机侧变流器K1包括模组K11、模组K12、模组K11的桥臂侧电感L3、模组K12的桥臂侧电感L4;所述网侧变流器K2包括模组K21、模组K22、模组K21的桥臂侧电感L1、模组K22的桥臂侧电感L2;所述桥臂侧电感L3与永磁同步电机定子侧相连,桥臂侧电感L4并接在桥臂侧电感L3与永磁同步电机M之间,机侧滤波电容C1也并接在臂侧电感L3与永磁同步电机M之间;桥臂侧电感L1串联网侧电感Lg后接入电网,桥臂侧电感L2并接在桥臂侧电感L1与网侧电感Lg之间,滤波电容C也并接在桥臂侧电感L1与网侧电感Lg之间;所述直流侧电容Cdc并联在机侧变流器K1和网侧变流器K2之间;
所述网侧滤波电容C与桥臂侧电感L1、桥臂侧电感L2构成了网侧滤波器;
将变流器模组K11向直流侧输出的电流记为模组K11机侧输出电流、变流器模组K12向直流侧输出的电流记为模组K12机侧输出电流、模组K21向电网侧输出的电流记为模组K21桥臂侧电感电流、模组K22向电网侧输出的电流记为模组K22桥臂侧电感电流;
所述建模方法包括建立以下7个模型:双dq坐标系下的网侧变流器控制模型、双dq坐标系下的网侧滤波器及电网阻抗模型、双dq坐标系下的网侧滤波器双dq变换时的解耦运算及滤波器模型、双dq坐标系下的电网接口模型、双dq坐标系下锁相环模型、双dq坐标系下各模组中间直流侧模型、双dq坐标系下机侧变流器控制模型;所述dq坐标系采用d轴定向、q轴超前于d轴90°;
所述建模方法的具体步骤如下:
步骤1,建立双dq坐标系下的网侧变流器控制模型,包括:建立正序dq坐标系下电压外环方程、负序dq坐标系下电压外环方程、正序dq坐标系下电流内环方程和负序dq坐标系下电流内环方程,具体过程及表达式如下:
步骤1.1,采样,即通过采集得到双dq坐标系下与网侧变流器控制相关的如下参数:直流侧电容Cdc两侧的直流侧电压Udc、模组K21桥臂侧正序电流d轴分量iLd1 p、模组K21桥臂侧正序电流q轴分量iLq1 p、模组K21桥臂侧负序电流d轴分量iLd1 n、模组K21桥臂侧负序电流q轴分量iLq1 n、模组K22桥臂侧正序电流d轴分量iLd2 p、模组K22桥臂侧正序电流q轴分量iLq2 p、模组K22桥臂侧负序电流d轴分量iLd2 n、模组K22桥臂侧负序电流q轴分量iLq2 n
步骤1.2,建立正序dq坐标系下电压外环方程,其表达式为:
Figure BDA0002706458700000051
式中:idref p为正序电流内环d轴参考值、Kvp为网侧变流器电压外环比例系数、Udcref为直流侧电压参考值、Udc为直流侧电压采样值、Ubase为电压基准值、intUdcerr为电压环积分项、Kvi为网侧变流器电压外环积分系数、iqref p为正序电流内环q轴参考值;
步骤1.3,建立负序dq坐标系下电压外环方程,其表达式为:
Figure BDA0002706458700000052
式中:idref n为负序电流内环d轴参考值、iqref n为负序电流内环q轴参考值;
步骤1.4,建立正序dq坐标系下电流内环方程,其表达式为:
Figure BDA0002706458700000053
Figure BDA0002706458700000054
其中:
Figure BDA0002706458700000055
式中:iod p为桥臂侧输出电流正序d轴分量、ioq p为桥臂侧输出电流正序q轴分量、iLd1 p为模组K21桥臂侧输出电流正序d轴分量、iLd2 p为模组K22桥臂侧输出电流正序d轴分量、iLq1 p为模组K21桥臂侧输出电流正序q轴分量、iLq2 p为模组K22桥臂侧输出电流正序q轴分量、ud p为电流内环输出正序电压d轴分量、uq p为电流内环输出正序电压q轴分量、Kpip为正序坐标系下电流内环PI调节器的比例系数、Ibase为电流基准值、intiderr p为网侧变流器电流环控制的积分项输出d轴正序分量、intiqerr p为网侧变流器电流环控制的积分项输出q轴正序分量、Kiip为正序坐标系下电流内环PI调节器的积分环节;
步骤1.5,建立负序dq坐标系下电流内环方程,其表达上为:
Figure BDA0002706458700000061
Figure BDA0002706458700000062
其中:
Figure BDA0002706458700000063
式中:iod n为桥臂侧输出电流负序d轴分量、ioq n为桥臂侧输出电流负序q轴分量、ud n为电流内环输出负序电压d轴分量、uq n为电流内环输出负序电压q轴分量、Kpin为负序坐标系下电流内环PI调节器的比例系数、intiderr n为网侧变流器电流环控制的积分项输出d轴负序分量、intiqerr n为网侧变流器电流环控制的积分项输出q轴负序分量、Kiin为负序坐标系下电流内环PI调节器的积分系数;
步骤2,建立双dq坐标系下的网侧滤波器及电网阻抗模型,包括:建立正序dq坐标系下的网侧滤波器及电网阻抗模型、负序dq坐标系下的网侧滤波器及电网阻抗模型,具体过程及表达式:
步骤2.1,采样,即通过采集获得双dq坐标系下与网侧滤波器及电网阻抗相关的如下参数:模组K21桥臂侧电感电流经正序解耦运算单元后获得的模组K21桥臂侧正序电流d轴分量iLd1 p_L、模组K21桥臂侧电感电流经正序解耦运算单元后获得的模组K21桥臂侧正序电流q轴分量iLq1 p_L、模组K21桥臂侧电感电流经负序解耦运算单元后获得的模组K21桥臂侧负序电流d轴分量iLd1 n_L、模组K21桥臂侧电感电流经负序解耦运算单元后获得的模组K21桥臂侧负序电流q轴分量iLq1 n_L、模组K22桥臂侧电感电流经正序解耦运算单元后获得的模组K22桥臂侧正序电流d轴分量iLd2 p_L、模组K22桥臂侧电感电流经正序解耦运算单元后获得的模组K22桥臂侧正序电流q轴分量iLq2 p_L、模组K22桥臂电感侧电流经负序解耦运算单元后获得的模组K22桥臂侧负序电流d轴分量iLd2 n_L、模组K22桥臂侧电感电流经负序解耦运算单元后获得的模组K22桥臂侧负序电流q轴分量iLq2 n_L、网侧滤波电容电压经正序解耦运算单元后获得的正序电容电压d轴分量ucd p_L、网侧滤波电容电压经正序解耦运算单元后获得的正序电容电压q轴分量ucq p_L、网侧滤波电容电压经负序解耦运算单元后获得的负序电容电压d轴分量ucd n_L、网侧滤波电容电压经负序解耦运算单元后获得的负序电容电压q轴分量ucq n _L、网侧电感电流经正序解耦运算单元后获得的正序网侧电感电流d轴分量igd p_L、网侧电感电流经正序解耦运算单元后获得的正序网侧电感电流q轴分量igq p_L、网侧电感电流经负序解耦运算单元后获得的负序网侧电感电流d轴分量igd n_L、网侧电感电流经负序解耦运算单元后获得的负序网侧电感电流q轴分量igq n_L、电网电压经正序解耦运算单元后获得的正序电网电压d轴分量ed p_L、电网电压经正序解耦运算单元后获得的正序电网电压q轴分量eq p_L、电网电压经负序解耦运算单元后获得的负序电网电压d轴分量ed n_L、电网电压经负序解耦运算单元后获得的负序电网电压q轴分量eq n_L
步骤2.2,建立正序dq坐标系下的网侧滤波器及电网阻抗模型,其表达式为:
Figure BDA0002706458700000081
Figure BDA0002706458700000082
其中,
Figure BDA0002706458700000083
式中:ω为锁相环输出的角频率、ud1 p为模组K21正序输出电压d轴分量、uq1 p为模组K21正序输出电压q轴分量、ud2 p为模组K22正序输出电压d轴分量、uq2 p为模组K22正序输出电压q轴分量、iod p_L为模组K21与模组K22桥臂侧电感电流之和经正序解耦运算单元后获得的正序电流d轴分量、ioq p_L为模组K21与模组K22桥臂侧电感电流之和经正序解耦运算单元后获得的正序电流q轴分量、R1为桥臂电侧感L1的寄生电阻、R2为桥臂侧电感L2的寄生电阻、Rg网为侧电感Lg的寄生电阻;
步骤2.3,建立负序dq坐标系下的网侧滤波器及电网阻抗模型,其表达式为:
Figure BDA0002706458700000091
Figure BDA0002706458700000092
其中:
Figure BDA0002706458700000093
式中:ud1 n为模组K21负序输出电压d轴分量、uq1 n为模组K21负序输出电压q轴分量、ud2 n为模组K22负序输出电压d轴分量、uq2 n为模组K22负序输出电压q轴分量、iod n_L为模组K21与模组K22桥臂侧电流之和经负序解耦运算单元后获得的负序电流d轴分量、ioq n_L为模组K21与模组K22桥臂侧电流之和经负序解耦运算单元后获得的负序电流q轴分量;
步骤3,建立双dq坐标系下的网侧滤波器双dq变换时的解耦运算及滤波器模型,包括:建立正负序dq坐标系下的桥臂侧电感电流双dq变换时的解耦运算及滤波器模型、建立正负序dq坐标系下的网侧滤波电容电压双dq变换时的解耦运算及滤波器模型、建立正负序dq坐标系下的网侧电感电流双dq变换时的解耦运算及滤波器模型,具体过程及表达式:
步骤3.1,建立正负序dq坐标系下的桥臂侧电感电流双dq变换时的解耦运算及滤波器模型,其表达式为:
模组K21:
Figure BDA0002706458700000101
Figure BDA0002706458700000102
式中:iLd1 p_1为模组K21桥臂侧电感电流经正序dq分解后的d轴分量、iLq1 p_1为模组K21桥臂侧电感电流经正序dq分解后的q轴分量、iLd1 n_1为模组K21桥臂侧电感电流经负序dq分解后的d轴分量、iLq1 n_1为模组K21桥臂侧电感电流经负序dq分解后的q轴分量、θPLL为锁相环输出的相角、wf为低通滤波器LPF截止频率;
模组K22:
Figure BDA0002706458700000111
Figure BDA0002706458700000112
式中:iLd2 p_1为模组K22桥臂侧电感电流经正序dq分解后的d轴分量、iLq2 p_1为模组K22桥臂侧电感电流经正序dq分解后的q轴分量、iLd2 n_1为模组K22桥臂侧电感电流经负序dq分解后的d轴分量、iLq2 n_1为模组K22桥臂侧电感电流经负序dq分解后的q轴分量;
步骤3.2,建立正负序dq坐标系下的网侧滤波电容电压双dq变换时的解耦运算及滤波器模型,其表达式为:
Figure BDA0002706458700000113
Figure BDA0002706458700000121
式中:ucd p_1为网侧滤波电容电压经正序dq分解后的d轴分量、ucq p_1为网侧滤波电容电压经正序dq分解后的q轴分量、ucd n_1为网侧滤波电容电压经负序dq分解后的d轴分量、ucq n_1为网侧滤波电容电压经负序dq分解后的q轴分量、ucd p为网侧滤波电容正序电容电压d轴分量、ucq p为网侧滤波电容正序电容电压q轴分量、ucd n为网侧滤波电容负序电容电压d轴分量、ucq n为网侧滤波电容负序电容电压q轴分量;
步骤3.3,建立正负序dq坐标系下的网侧电感电流双dq变换时的解耦运算及滤波器模型,其表达式为:
Figure BDA0002706458700000122
Figure BDA0002706458700000123
式中:igd p_1为网侧电感电流经正序dq分解后的d轴分量、igq p_1为网侧电感电流经正序dq分解后的q轴分量、igd n_1为网侧电感电流经负序dq分解后的d轴分量、igq n_1为网侧电感电流经负序dq分解后的q轴分量、igd p为网侧电感的正序电感电流d轴分量、igq p为网侧电感的正序电感电流q轴分量、igd n为网侧电感的负序电感电流d轴分量、igq n_1为网侧电感的负序电感电流q轴分量;
步骤4,建立双dq坐标系下的电网接口模型,包括:建立正序dq坐标系下的电网接口模型、负序dq坐标系下的电网接口模型、双dq坐标系下电网接口正序和负序dq变换时的解耦运算及滤波器模型,具体过程及表达式如下:
步骤4.1,建立正序dq坐标系下的电网接口模型,其表达式为:
ed p_1=Upmcos(θ+θp)+Unmcos(θ+2θPLLn);
eq p_1=Upmsin(θ+θp)-Unmsin(θ+2θPLLn);
式中:ed p_1为电网电压经正序dq变换得到的d轴分量、eq p_1为电网电压经正序dq变换得到的q轴分量、Upm为正序电压幅值、Unm为负序电压幅值、θ为通过锁相环得到的角度与电网角度的偏差角度、θp为正序电网的初始角度、θn为负序电网的初始角度;
步骤4.2,建立负序dq坐标系下的电网接口模型,其表达式为:
ed n_1=Unmcos(θ+θn)+Upmcos(θ+2θPLLp);
eq n_1=Unmsin(θ-θn)+Upmsin(θ+2θPLLp);
式中:ed n_1为电网电压经负序dq变换得到的d轴分量、eq n_1为电网电压经负序dq变换得到的q轴分量;
步骤4.3,建立双dq坐标系下电网接口解耦运算及滤波器模型,其表达式为:
Figure BDA0002706458700000141
其中:
Figure BDA0002706458700000142
式中:ed p_1为电网电压经正序dq变换得到的d轴分量、eq p_1为电网电压经正序dq变换得到的q轴分量、ed n_1为电网电压经负序dq变换得到的d轴分量、eq n_1为电网电压经负序dq变换得到的q轴分量、ed p为电网电压正序电压d轴分量、eq p为电网电压正序电压q轴分量、ed n为电网电压负序电压d轴分量、eq n为电网电压负序电压q轴分量;
步骤5,建立双dq坐标系下锁相环模型,包括:
Figure BDA0002706458700000143
式中:ω为锁相环输出的同步角频率、ω0为额定角频率、Kppll为锁相环比例系数、Kipll为锁相环积分系数、intucq p_L为锁相环控制的积分项输出;
步骤6,建立双dq坐标系下各模组中间直流侧模型,包括建立双dq坐标系下中间直流侧模型方程、网侧输出功率方程,具体过程及表达式如下:
步骤6.1,建立双dq坐标系下中间直流侧模型方程,其表达式为:
Figure BDA0002706458700000151
Figure BDA0002706458700000152
式中:Pdc11为模组K21向网侧输出的功率、Pdc12为模组K11向中间直流侧输入的功率、idc1d为模组K11向直流侧输入的电流d轴分量、Pdc21为模组K22向网侧输出的功率、Pdc22为模组K12向中间直流侧输入的功率、idc2d为模组K12向中间直流侧输入电流d轴分量;
步骤6.2,建立网侧输出功率方程,其表达式为:
模组K21:
Figure BDA0002706458700000153
其中:
Figure BDA0002706458700000154
式中:Pdc11o为模组K21输出功率的平均值、Pdc11c2为模组K21输出功率余弦量2次脉动的幅值、Pdc11s2为模组K21输出功率正弦量2次脉动的幅值;
模组K22:
Figure BDA0002706458700000155
其中:
Figure BDA0002706458700000161
式中:Pdc21o为模组K22输出功率的平均值、Pdc21c2为模组K22输出功率余弦量2次脉动的幅值、Pdc21s2为模组K22输出功率正弦量2次脉动的幅值;
步骤7,建立单dq坐标系下机侧变流器、永磁同步电机及风力机等效模型:
将单dq坐标系下的机侧变流器、永磁同步电机及风力机等效成一个已知电流源模型,则:
Figure BDA0002706458700000162
式中:idc1d为模组K11机侧输出电流的d轴分量,idc1q为模组K11机侧输出电流的q轴分量,idc2d为模组K12机侧输出电流的d轴分量,idc2q为模组K12机侧输出电流的q轴分量。
与现有技术相比,本发明的有益效果如下:
1、现今对直驱型全功率风力发电系统,尤其是采用双dq控制的直驱型全功率风力发电系统模型的研究甚少,本发明提出的基于双dq坐标系的全功率风力发电系统建模方法有别于传统的建模方法,小信号建模过程是基于变流器模组并联得到的,不仅符合当下大功率风电变流器提高功率容量的方式,而且也能适用单模组变流器的情形,建模过程具有极强的适应性,非常符合现场实际并具有创新意义。
2、由于直驱型全功率风力发电系统,机侧永磁同步电机经过背靠背PWM变流器同电网相连,实现了机侧与网侧的隔离,因此,本发明中将机侧等效为一个电流源,大大简化了运算过程和系统的阶数。
附图说明
图1是本发明实施例的直驱型全功率风力发电系统电路结构拓扑图。
图2是本发明实施例中的直驱型全功率变流器网侧控制框图。
图3是本发明的在双同步旋转坐标系下正负序分量获取原理框图。
图4是本发明的在双同步旋转坐标系下的锁相环原理框图。
具体实施方式
图1是本发明实施例的直驱型全功率风力发电系统电路结构拓扑图。由该图可见,本发明所述全功率风力发电系统为变流器模组输入输出并联的直驱型全功率风力发电系统。所述变流器模组输入输出并联的直驱型全功率风力发电系统主电路的拓扑结构包括永磁同步电机M、机侧滤波电容C1、机侧滤波电容C1的寄生电阻R、直流侧电容Cdc、机侧变流器K1、网侧变流器K2、网侧滤波电容C、网侧电感Lg、网侧电感Lg的寄生电阻Rg和电网。
所述机侧变流器K1包括模组K11、模组K12、模组K11的桥臂侧电感L3、模组K12的桥臂侧电感L4。所述网侧变流器K2包括模组K21、模组K22、模组K21的桥臂侧电感L1、桥臂电侧感L1的寄生电阻R1、模组K22的桥臂侧电感L2、桥臂侧电感L2的寄生电阻R2。所述桥臂侧电感L3与永磁同步电机定子侧相连,桥臂侧电感L4并接在桥臂侧电感L3与永磁同步电机M之间,机侧滤波电容C1也并接在臂侧电感L3与永磁同步电机M之间。桥臂侧电感L1串联网侧电感Lg后接入电网,桥臂侧电感L2并接在桥臂侧电感L1与网侧电感Lg之间,滤波电容C也并接在桥臂侧电感L1与网侧电感Lg之间。所述直流侧电容Cdc并联在机侧变流器K1和网侧变流器K2之间。
所述网侧滤波电容C与桥臂侧电感L1、桥臂侧电感L2构成了网侧滤波器。
将变流器模组K11向直流侧输出的电流记为模组K11机侧输出电流、变流器模组K12向直流侧输出的电流记为模组K12机侧输出电流、模组K21向电网侧输出的电流记为模组K21桥臂侧电感电流iL1、模组K22向电网侧输出的电流记为模组K22桥臂侧电感电流iL2
本发明所述建模方法包括建立以下7个模型:双dq坐标系下的网侧变流器控制模型、双dq坐标系下的网侧滤波器及电网阻抗模型、双dq坐标系下的网侧滤波器双dq变换时的解耦运算及滤波器模型、双dq坐标系下的电网接口模型、双dq坐标系下锁相环模型、双dq坐标系下各模组中间直流侧模型、双dq坐标系下机侧变流器控制模型;所述dq坐标系采用d轴定向、q轴超前于d轴90°。
所述建模方法的具体步骤如下:
步骤1,建立双dq坐标系下的网侧变流器控制模型,包括:建立正序dq坐标系下电压外环方程、负序dq坐标系下电压外环方程、正序dq坐标系下电流内环方程和负序dq坐标系下电流内环方程。双dq坐标系下的网侧变流器控制如图2所示。
具体过程及表达式如下:
步骤1.1,采样,即通过采集得到双dq坐标系下与网侧变流器控制相关的如下参数:直流侧电容Cdc两侧的直流侧电压Udc、模组K21桥臂侧正序电流d轴分量iLd1 p、模组K21桥臂侧正序电流q轴分量iLq1 p、模组K21桥臂侧负序电流d轴分量iLd1 n、模组K21桥臂侧负序电流q轴分量iLq1 n、模组K22桥臂侧正序电流d轴分量iLd2 p、模组K22桥臂侧正序电流q轴分量iLq2 p、模组K22桥臂侧负序电流d轴分量iLd2 n、模组K22桥臂侧负序电流q轴分量iLq2 n
步骤1.2,建立正序dq坐标系下电压外环方程,其表达式为:
Figure BDA0002706458700000191
式中:idref p为正序电流内环d轴参考值、Kvp为网侧变流器电压外环比例系数、Udcref为直流侧电压参考值、Udc为直流侧电压采样值、Ubase为电压基准值、intUdcerr为电压环积分项、Kvi为网侧变流器电压外环积分系数、iqref p为正序电流内环q轴参考值;
步骤1.3,建立负序dq坐标系下电压外环方程,其表达式为:
Figure BDA0002706458700000192
式中:idref n为负序电流内环d轴参考值、iqref n为负序电流内环q轴参考值;
步骤1.4,建立正序dq坐标系下电流内环方程,其表达式为:
Figure BDA0002706458700000193
Figure BDA0002706458700000194
其中:
Figure BDA0002706458700000195
式中:iod p为桥臂侧输出电流正序d轴分量、ioq p为桥臂侧输出电流正序q轴分量、iLd1 p为模组K21桥臂侧输出电流正序d轴分量、iLd2 p为模组K22桥臂侧输出电流正序d轴分量、iLq1 p为模组K21桥臂侧输出电流正序q轴分量、iLq2 p为模组K22桥臂侧输出电流正序q轴分量、ud p为电流内环输出正序电压d轴分量、uq p为电流内环输出正序电压q轴分量、Kpip为正序坐标系下电流内环PI调节器的比例系数、Ibase为电流基准值、intiderr p为网侧变流器电流环控制的积分项输出d轴正序分量、intiqerr p为网侧变流器电流环控制的积分项输出q轴正序分量、Kiip为正序坐标系下电流内环PI调节器的积分环节;
步骤1.5,建立负序dq坐标系下电流内环方程,其表达上为:
Figure BDA0002706458700000201
Figure BDA0002706458700000202
其中:
Figure BDA0002706458700000203
式中:iod n为桥臂侧输出电流负序d轴分量、ioq n为桥臂侧输出电流负序q轴分量、ud n为电流内环输出负序电压d轴分量、uq n为电流内环输出负序电压q轴分量、Kpin为负序坐标系下电流内环PI调节器的比例系数、intiderr n为网侧变流器电流环控制的积分项输出d轴负序分量、intiqerr n为网侧变流器电流环控制的积分项输出q轴负序分量、Kiin为负序坐标系下电流内环PI调节器的积分系数;
步骤2,建立双dq坐标系下的网侧滤波器及电网阻抗模型,包括:建立正序dq坐标系下的网侧滤波器及电网阻抗模型、负序dq坐标系下的网侧滤波器及电网阻抗模型。网侧滤波器及电网阻抗连接关系如图2所示。
具体过程及表达式:
步骤2.1,采样,即通过采集获得双dq坐标系下与网侧滤波器及电网阻抗相关的如下参数:模组K21桥臂侧电感电流经正序解耦运算单元后获得的模组K21桥臂侧正序电流d轴分量iLd1 p_L、模组K21桥臂侧电感电流经正序解耦运算单元后获得的模组K21桥臂侧正序电流q轴分量iLq1 p_L、模组K21桥臂侧电感电流经负序解耦运算单元后获得的模组K21桥臂侧负序电流d轴分量iLd1 n_L、模组K21桥臂侧电感电流经负序解耦运算单元后获得的模组K21桥臂侧负序电流q轴分量iLq1 n_L、模组K22桥臂侧电感电流经正序解耦运算单元后获得的模组K22桥臂侧正序电流d轴分量iLd2 p_L、模组K22桥臂侧电感电流经正序解耦运算单元后获得的模组K22桥臂侧正序电流q轴分量iLq2 p_L、模组K22桥臂电感侧电流经负序解耦运算单元后获得的模组K22桥臂侧负序电流d轴分量iLd2 n_L、模组K22桥臂侧电感电流经负序解耦运算单元后获得的模组K22桥臂侧负序电流q轴分量iLq2 n_L、网侧滤波电容电压经正序解耦运算单元后获得的正序电容电压d轴分量ucd p_L、网侧滤波电容电压经正序解耦运算单元后获得的正序电容电压q轴分量ucq p_L、网侧滤波电容电压经负序解耦运算单元后获得的负序电容电压d轴分量ucd n_L、网侧滤波电容电压经负序解耦运算单元后获得的负序电容电压q轴分量ucq n _L、网侧电感电流经正序解耦运算单元后获得的正序网侧电感电流d轴分量igd p_L、网侧电感电流经正序解耦运算单元后获得的正序网侧电感电流q轴分量igq p_L、网侧电感电流经负序解耦运算单元后获得的负序网侧电感电流d轴分量igd n_L、网侧电感电流经负序解耦运算单元后获得的负序网侧电感电流q轴分量igq n_L、电网电压经正序解耦运算单元后获得的正序电网电压d轴分量ed p_L、电网电压经正序解耦运算单元后获得的正序电网电压q轴分量eq p_L、电网电压经负序解耦运算单元后获得的负序电网电压d轴分量ed n_L、电网电压经负序解耦运算单元后获得的负序电网电压q轴分量eq n_L
步骤2.2,建立正序dq坐标系下的网侧滤波器及电网阻抗模型,其表达式为:
Figure BDA0002706458700000221
Figure BDA0002706458700000222
其中,
Figure BDA0002706458700000223
式中:ω为锁相环输出的角频率、ud1 p为模组K21正序输出电压d轴分量、uq1 p为模组K21正序输出电压q轴分量、ud2 p为模组K22正序输出电压d轴分量、uq2 p为模组K22正序输出电压q轴分量、iod p_L为模组K21与模组K22桥臂侧电感电流之和经正序解耦运算单元后获得的正序电流d轴分量、ioq p_L为模组K21与模组K22桥臂侧电感电流之和经正序解耦运算单元后获得的正序电流q轴分量、R1为桥臂电侧感L1的寄生电阻、R2为桥臂侧电感L2的寄生电阻、Rg网为侧电感Lg的寄生电阻;
步骤2.3,建立负序dq坐标系下的网侧滤波器及电网阻抗模型,其表达式为:
Figure BDA0002706458700000231
Figure BDA0002706458700000232
其中:
Figure BDA0002706458700000233
式中:ud1 n为模组K21负序输出电压d轴分量、uq1 n为模组K21负序输出电压q轴分量、ud2 n为模组K22负序输出电压d轴分量、uq2 n为模组K22负序输出电压q轴分量、iod n_L为模组K21与模组K22桥臂侧电流之和经负序解耦运算单元后获得的负序电流d轴分量、ioq n_L为模组K21与模组K22桥臂侧电流之和经负序解耦运算单元后获得的负序电流q轴分量;
步骤3,建立双dq坐标系下的网侧滤波器双dq变换时的解耦运算及滤波器模型,包括:建立正负序dq坐标系下的桥臂侧电感电流双dq变换时的解耦运算及滤波器模型、建立正负序dq坐标系下的网侧滤波电容电压双dq变换时的解耦运算及滤波器模型、建立正负序dq坐标系下的网侧电感电流双dq变换时的解耦运算及滤波器模型。双dq坐标系下变量通过正负序dq变换时的解耦运算及滤波器框图如图3所示。
具体过程及表达式:
步骤3.1,建立正负序dq坐标系下的桥臂侧电感电流双dq变换时的解耦运算及滤波器模型,其表达式为:
模组K21:
Figure BDA0002706458700000241
Figure BDA0002706458700000242
式中:iLd1 p_1为模组K21桥臂侧电感电流经正序dq分解后的d轴分量、iLq1 p_1为模组K21桥臂侧电感电流经正序dq分解后的q轴分量、iLd1 n_1为模组K21桥臂侧电感电流经负序dq分解后的d轴分量、iLq1 n_1为模组K21桥臂侧电感电流经负序dq分解后的q轴分量、θPLL为锁相环输出的相角、wf为低通滤波器LPF截止频率;
模组K22:
Figure BDA0002706458700000251
Figure BDA0002706458700000252
式中:iLd2 p_1为模组K22桥臂侧电感电流经正序dq分解后的d轴分量、iLq2 p_1为模组K22桥臂侧电感电流经正序dq分解后的q轴分量、iLd2 n_1为模组K22桥臂侧电感电流经负序dq分解后的d轴分量、iLq2 n_1为模组K22桥臂侧电感电流经负序dq分解后的q轴分量;
步骤3.2,建立正负序dq坐标系下的网侧滤波电容电压双dq变换时的解耦运算及滤波器模型,其表达式为:
Figure BDA0002706458700000253
Figure BDA0002706458700000261
式中:ucd p_1为网侧滤波电容电压经正序dq分解后的d轴分量、ucq p_1为网侧滤波电容电压经正序dq分解后的q轴分量、ucd n_1为网侧滤波电容电压经负序dq分解后的d轴分量、ucq n_1为网侧滤波电容电压经负序dq分解后的q轴分量、ucd p为网侧滤波电容正序电容电压d轴分量、ucq p为网侧滤波电容正序电容电压q轴分量、ucd n为网侧滤波电容负序电容电压d轴分量、ucq n为网侧滤波电容负序电容电压q轴分量;
步骤3.3,建立正负序dq坐标系下的网侧电感电流双dq变换时的解耦运算及滤波器模型,其表达式为:
Figure BDA0002706458700000262
Figure BDA0002706458700000263
式中:igd p_1为网侧电感电流经正序dq分解后的d轴分量、igq p_1为网侧电感电流经正序dq分解后的q轴分量、igd n_1为网侧电感电流经负序dq分解后的d轴分量、igq n_1为网侧电感电流经负序dq分解后的q轴分量、igd p为网侧电感的正序电感电流d轴分量、igq p为网侧电感的正序电感电流q轴分量、igd n为网侧电感的负序电感电流d轴分量、igq n_1为网侧电感的负序电感电流q轴分量;
步骤4,建立双dq坐标系下的电网接口模型,包括:建立正序dq坐标系下的电网接口模型、负序dq坐标系下的电网接口模型、双dq坐标系下电网接口正序和负序dq变换时的解耦运算及滤波器模型。双dq坐标系下变量通过正负序dq变换时的解耦运算及滤波器框图如图3所示。
具体过程及表达式如下:
步骤4.1,建立正序dq坐标系下的电网接口模型,其表达式为:
ed p_1=Upmcos(θ+θp)+Unmcos(θ+2θPLLn);
eq p_1=Upmsin(θ+θp)-Unmsin(θ+2θPLLn);
式中:ed p_1为电网电压经正序dq变换得到的d轴分量、eq p_1为电网电压经正序dq变换得到的q轴分量、Upm为正序电压幅值、Unm为负序电压幅值、θ为通过锁相环得到的角度与电网角度的偏差角度、θp为正序电网的初始角度、θn为负序电网的初始角度;
步骤4.2,建立负序dq坐标系下的电网接口模型,其表达式为:
ed n_1=Unmcos(θ+θn)+Upmcos(θ+2θPLLp);
eq n_1=Unmsin(θ-θn)+Upmsin(θ+2θPLLp);
式中:ed n_1为电网电压经负序dq变换得到的d轴分量、eq n_1为电网电压经负序dq变换得到的q轴分量;
步骤4.3,建立双dq坐标系下电网接口解耦运算及滤波器模型,其表达式为:
Figure BDA0002706458700000281
其中:
Figure BDA0002706458700000282
式中:ed p_1为电网电压经正序dq变换得到的d轴分量、eq p_1为电网电压经正序dq变换得到的q轴分量、ed n_1为电网电压经负序dq变换得到的d轴分量、eq n_1为电网电压经负序dq变换得到的q轴分量、ed p为电网电压正序电压d轴分量、eq p为电网电压正序电压q轴分量、ed n为电网电压负序电压d轴分量、eq n为电网电压负序电压q轴分量;
步骤5,建立双dq坐标系下锁相环模型,双dq坐标系下锁相环框图如图4所示,包括:
Figure BDA0002706458700000283
式中:ω为锁相环输出的同步角频率、ω0为额定角频率、Kppll为锁相环比例系数、Kipll为锁相环积分系数、intucq p_L为锁相环控制的积分项输出;
步骤6,建立双dq坐标系下各模组中间直流侧模型,包括建立双dq坐标系下中间直流侧模型方程、网侧输出功率方程。
具体过程及表达式如下:
步骤6.1,建立双dq坐标系下中间直流侧模型方程,其表达式为:
Figure BDA0002706458700000291
Figure BDA0002706458700000292
式中:Pdc11为模组K21向网侧输出的功率、Pdc12为模组K11向中间直流侧输入的功率、idc1d为模组K11向直流侧输入的电流d轴分量、Pdc21为模组K22向网侧输出的功率、Pdc22为模组K12向中间直流侧输入的功率、idc2d为模组K12向中间直流侧输入电流d轴分量;
步骤6.2,建立网侧输出功率方程,其表达式为:
模组K21:
Figure BDA0002706458700000293
其中:
Figure BDA0002706458700000294
式中:Pdc11o为模组K21输出功率的平均值、Pdc11c2为模组K21输出功率余弦量2次脉动的幅值、Pdc11s2为模组K21输出功率正弦量2次脉动的幅值;
模组K22:
Figure BDA0002706458700000301
其中:
Figure BDA0002706458700000302
式中:Pdc21o为模组K22输出功率的平均值、Pdc21c2为模组K22输出功率余弦量2次脉动的幅值、Pdc21s2为模组K22输出功率正弦量2次脉动的幅值;
步骤7,建立单dq坐标系下机侧变流器、永磁同步电机及风力机等效模型:
将单dq坐标系下的机侧变流器、永磁同步电机及风力机等效成一个已知电流源模型,则:
Figure BDA0002706458700000303
式中:idc1d为模组K11机侧输出电流的d轴分量,idc1q为模组K11机侧输出电流的q轴分量,idc2d为模组K12机侧输出电流的d轴分量,idc2q为模组K12机侧输出电流的q轴分量。

Claims (1)

1.一种基于双同步旋转坐标系的全功率风力发电系统的建模方法,所述全功率风力发电系统为变流器模组输入输出并联的直驱型全功率风力发电系统;所述变流器模组输入输出并联的直驱型全功率风力发电系统主电路的拓扑结构包括永磁同步电机M、机侧滤波电容C1、直流侧电容Cdc、机侧变流器K1、网侧变流器K2、网侧滤波电容C、网侧电感Lg和电网;所述机侧变流器K1包括模组K11、模组K12、模组K11的桥臂侧电感L3、模组K12的桥臂侧电感L4;所述网侧变流器K2包括模组K21、模组K22、模组K21的桥臂侧电感L1、模组K22的桥臂侧电感L2;所述桥臂侧电感L3与永磁同步电机定子侧相连,桥臂侧电感L4并接在桥臂侧电感L3与永磁同步电机M之间,机侧滤波电容C1也并接在桥臂侧电感L3与永磁同步电机M之间;桥臂侧电感L1串联网侧电感Lg后接入电网,桥臂侧电感L2并接在桥臂侧电感L1与网侧电感Lg之间,滤波电容C也并接在桥臂侧电感L1与网侧电感Lg之间;所述直流侧电容Cdc并联在机侧变流器K1和网侧变流器K2之间;
所述网侧滤波电容C与桥臂侧电感L1、桥臂侧电感L2构成了网侧滤波器;
将变流器模组K11向直流侧输出的电流记为模组K11机侧输出电流、变流器模组K12向直流侧输出的电流记为模组K12机侧输出电流、模组K21向电网侧输出的电流记为模组K21桥臂侧电感电流、模组K22向电网侧输出的电流记为模组K22桥臂侧电感电流;
其特征在于,所述建模方法包括建立以下7个模型:双dq坐标系下的网侧变流器控制模型、双dq坐标系下的网侧滤波器及电网阻抗模型、双dq坐标系下的网侧滤波器双dq变换时的解耦运算及滤波器模型、双dq坐标系下的电网接口模型、双dq坐标系下锁相环模型、双dq坐标系下各模组中间直流侧模型、双dq坐标系下机侧变流器控制模型;所述dq坐标系采用d轴定向、q轴超前于d轴90°;
所述建模方法的具体步骤如下:
步骤1,建立双dq坐标系下的网侧变流器控制模型,包括:建立正序dq坐标系下电压外环方程、负序dq坐标系下电压外环方程、正序dq坐标系下电流内环方程和负序dq坐标系下电流内环方程,具体过程及表达式如下:
步骤1.1,采样,即通过采集得到双dq坐标系下与网侧变流器控制相关的如下参数:直流侧电容Cdc两侧的直流侧电压Udc、模组K21桥臂侧正序电流d轴分量iLd1 p、模组K21桥臂侧正序电流q轴分量iLq1 p、模组K21桥臂侧负序电流d轴分量iLd1 n、模组K21桥臂侧负序电流q轴分量iLq1 n、模组K22桥臂侧正序电流d轴分量iLd2 p、模组K22桥臂侧正序电流q轴分量iLq2 p、模组K22桥臂侧负序电流d轴分量iLd2 n、模组K22桥臂侧负序电流q轴分量iLq2 n
步骤1.2,建立正序dq坐标系下电压外环方程,其表达式为:
Figure FDA0003629000170000021
式中:idref p为正序电流内环d轴参考值、Kvp为网侧变流器电压外环比例系数、Udcref为直流侧电压参考值、Udc为直流侧电压采样值、Ubase为电压基准值、intUdcerr为电压环积分项、Kvi为网侧变流器电压外环积分系数、iqref p为正序电流内环q轴参考值;
步骤1.3,建立负序dq坐标系下电压外环方程,其表达式为:
Figure FDA0003629000170000022
式中:idref n为负序电流内环d轴参考值、iqref n为负序电流内环q轴参考值;
步骤1.4,建立正序dq坐标系下电流内环方程,其表达式为:
Figure FDA0003629000170000031
Figure FDA0003629000170000032
其中:
Figure FDA0003629000170000033
式中:iod p为桥臂侧输出电流正序d轴分量、ioq p为桥臂侧输出电流正序q轴分量、iLd1 p为模组K21桥臂侧输出电流正序d轴分量、iLd2 p为模组K22桥臂侧输出电流正序d轴分量、iLq1 p为模组K21桥臂侧输出电流正序q轴分量、iLq2 p为模组K22桥臂侧输出电流正序q轴分量、ud p为电流内环输出正序电压d轴分量、uq p为电流内环输出正序电压q轴分量、Kpip为正序坐标系下电流内环PI调节器的比例系数、Ibase为电流基准值、intiderr p为网侧变流器电流环控制的积分项输出d轴正序分量、intiqerr p为网侧变流器电流环控制的积分项输出q轴正序分量、Kiip为正序坐标系下电流内环PI调节器的积分环节;
步骤1.5,建立负序dq坐标系下电流内环方程,其表达式为:
Figure FDA0003629000170000034
Figure FDA0003629000170000035
其中:
Figure FDA0003629000170000041
式中:iod n为桥臂侧输出电流负序d轴分量、ioq n为桥臂侧输出电流负序q轴分量、ud n为电流内环输出负序电压d轴分量、uq n为电流内环输出负序电压q轴分量、Kpin为负序坐标系下电流内环PI调节器的比例系数、intiderr n为网侧变流器电流环控制的积分项输出d轴负序分量、intiqerr n为网侧变流器电流环控制的积分项输出q轴负序分量、Kiin为负序坐标系下电流内环PI调节器的积分系数;
步骤2,建立双dq坐标系下的网侧滤波器及电网阻抗模型,包括:建立正序dq坐标系下的网侧滤波器及电网阻抗模型、负序dq坐标系下的网侧滤波器及电网阻抗模型,具体过程及表达式:
步骤2.1,采样,即通过采集获得双dq坐标系下与网侧滤波器及电网阻抗相关的如下参数:模组K21桥臂侧电感电流经正序解耦运算单元后获得的模组K21桥臂侧正序电流d轴分量iLd1 p_L、模组K21桥臂侧电感电流经正序解耦运算单元后获得的模组K21桥臂侧正序电流q轴分量iLq1 p_L、模组K21桥臂侧电感电流经负序解耦运算单元后获得的模组K21桥臂侧负序电流d轴分量iLd1 n_L、模组K21桥臂侧电感电流经负序解耦运算单元后获得的模组K21桥臂侧负序电流q轴分量iLq1 n_L、模组K22桥臂侧电感电流经正序解耦运算单元后获得的模组K22桥臂侧正序电流d轴分量iLd2 p_L、模组K22桥臂侧电感电流经正序解耦运算单元后获得的模组K22桥臂侧正序电流q轴分量iLq2 p_L、模组K22桥臂电感侧电流经负序解耦运算单元后获得的模组K22桥臂侧负序电流d轴分量iLd2 n_L、模组K22桥臂侧电感电流经负序解耦运算单元后获得的模组K22桥臂侧负序电流q轴分量iLq2 n_L、网侧滤波电容电压经正序解耦运算单元后获得的正序电容电压d轴分量ucd p_L、网侧滤波电容电压经正序解耦运算单元后获得的正序电容电压q轴分量ucq p_L、网侧滤波电容电压经负序解耦运算单元后获得的负序电容电压d轴分量ucd n_L、网侧滤波电容电压经负序解耦运算单元后获得的负序电容电压q轴分量ucq n_L、网侧电感电流经正序解耦运算单元后获得的正序网侧电感电流d轴分量igd p_L、网侧电感电流经正序解耦运算单元后获得的正序网侧电感电流q轴分量igq p_L、网侧电感电流经负序解耦运算单元后获得的负序网侧电感电流d轴分量igd n_L、网侧电感电流经负序解耦运算单元后获得的负序网侧电感电流q轴分量igq n_L、电网电压经正序解耦运算单元后获得的正序电网电压d轴分量ed p_L、电网电压经正序解耦运算单元后获得的正序电网电压q轴分量eq p_L、电网电压经负序解耦运算单元后获得的负序电网电压d轴分量ed n_L、电网电压经负序解耦运算单元后获得的负序电网电压q轴分量eq n_L
步骤2.2,建立正序dq坐标系下的网侧滤波器及电网阻抗模型,其表达式为:
Figure FDA0003629000170000051
Figure FDA0003629000170000061
其中,
Figure FDA0003629000170000062
式中:ω为锁相环输出的角频率、ud1 p为模组K21正序输出电压d轴分量、uq1 p为模组K21正序输出电压q轴分量、ud2 p为模组K22正序输出电压d轴分量、uq2 p为模组K22正序输出电压q轴分量、iod p_L为模组K21与模组K22桥臂侧电感电流之和经正序解耦运算单元后获得的正序电流d轴分量、ioq p_L为模组K21与模组K22桥臂侧电感电流之和经正序解耦运算单元后获得的正序电流q轴分量、R1为桥臂电侧感L1的寄生电阻、R2为桥臂侧电感L2的寄生电阻、Rg为网侧电感Lg的寄生电阻;
步骤2.3,建立负序dq坐标系下的网侧滤波器及电网阻抗模型,其表达式为:
Figure FDA0003629000170000063
Figure FDA0003629000170000071
其中:
Figure FDA0003629000170000072
式中:ud1 n为模组K21负序输出电压d轴分量、uq1 n为模组K21负序输出电压q轴分量、ud2 n为模组K22负序输出电压d轴分量、uq2 n为模组K22负序输出电压q轴分量、iod n_L为模组K21与模组K22桥臂侧电流之和经负序解耦运算单元后获得的负序电流d轴分量、ioq n_L为模组K21与模组K22桥臂侧电流之和经负序解耦运算单元后获得的负序电流q轴分量;
步骤3,建立双dq坐标系下的网侧滤波器双dq变换时的解耦运算及滤波器模型,包括:建立正负序dq坐标系下的桥臂侧电感电流双dq变换时的解耦运算及滤波器模型、建立正负序dq坐标系下的网侧滤波电容电压双dq变换时的解耦运算及滤波器模型、建立正负序dq坐标系下的网侧电感电流双dq变换时的解耦运算及滤波器模型,具体过程及表达式:
步骤3.1,建立正负序dq坐标系下的桥臂侧电感电流双dq变换时的解耦运算及滤波器模型,其表达式为:
模组K21:
Figure FDA0003629000170000081
Figure FDA0003629000170000082
式中:iLd1 p_1为模组K21桥臂侧电感电流经正序dq分解后的d轴分量、iLq1 p_1为模组K21桥臂侧电感电流经正序dq分解后的q轴分量、iLd1 n_1为模组K21桥臂侧电感电流经负序dq分解后的d轴分量、iLq1 n_1为模组K21桥臂侧电感电流经负序dq分解后的q轴分量、θPLL为锁相环输出的相角、wf为低通滤波器LPF截止频率;
模组K22:
Figure FDA0003629000170000083
Figure FDA0003629000170000084
式中:iLd2 p_1为模组K22桥臂侧电感电流经正序dq分解后的d轴分量、iLq2 p_1为模组K22桥臂侧电感电流经正序dq分解后的q轴分量、iLd2 n_1为模组K22桥臂侧电感电流经负序dq分解后的d轴分量、iLq2 n_1为模组K22桥臂侧电感电流经负序dq分解后的q轴分量;
步骤3.2,建立正负序dq坐标系下的网侧滤波电容电压双dq变换时的解耦运算及滤波器模型,其表达式为:
Figure FDA0003629000170000091
Figure FDA0003629000170000092
式中:ucd p_1为网侧滤波电容电压经正序dq分解后的d轴分量、ucq p_1为网侧滤波电容电压经正序dq分解后的q轴分量、ucd n_1为网侧滤波电容电压经负序dq分解后的d轴分量、ucq n_1为网侧滤波电容电压经负序dq分解后的q轴分量、ucd p为网侧滤波电容正序电容电压d轴分量、ucq p为网侧滤波电容正序电容电压q轴分量、ucd n为网侧滤波电容负序电容电压d轴分量、ucq n为网侧滤波电容负序电容电压q轴分量;
步骤3.3,建立正负序dq坐标系下的网侧电感电流双dq变换时的解耦运算及滤波器模型,其表达式为:
Figure FDA0003629000170000101
Figure FDA0003629000170000102
式中:igd p_1为网侧电感电流经正序dq分解后的d轴分量、igq p_1为网侧电感电流经正序dq分解后的q轴分量、igd n_1为网侧电感电流经负序dq分解后的d轴分量、igq n_1为网侧电感电流经负序dq分解后的q轴分量、igd p为网侧电感的正序电感电流d轴分量、igq p为网侧电感的正序电感电流q轴分量、igd n为网侧电感的负序电感电流d轴分量、igq n为网侧电感的负序电感电流q轴分量;
步骤4,建立双dq坐标系下的电网接口模型,包括:建立正序dq坐标系下的电网接口模型、负序dq坐标系下的电网接口模型、双dq坐标系下电网接口正序和负序dq变换时的解耦运算及滤波器模型,具体过程及表达式如下:
步骤4.1,建立正序dq坐标系下的电网接口模型,其表达式为:
ed p_1=Upmcos(θ+θp)+Unmcos(θ+2θPLLn);
eq p_1=Upmsin(θ+θp)-Unmsin(θ+2θPLLn);
式中:ed p_1为电网电压经正序dq变换得到的d轴分量、eq p_1为电网电压经正序dq变换得到的q轴分量、Upm为正序电压幅值、Unm为负序电压幅值、θ为通过锁相环得到的角度与电网角度的偏差角度、θp为正序电网的初始角度、θn为负序电网的初始角度;
步骤4.2,建立负序dq坐标系下的电网接口模型,其表达式为:
ed n_1=Unmcos(θ+θn)+Upmcos(θ+2θPLLp);
eq n_1=Unmsin(θ-θn)+Upmsin(θ+2θPLLp);
式中:ed n_1为电网电压经负序dq变换得到的d轴分量、eq n_1为电网电压经负序dq变换得到的q轴分量;
步骤4.3,建立双dq坐标系下电网接口解耦运算及滤波器模型,其表达式为:
Figure FDA0003629000170000111
其中:
Figure FDA0003629000170000112
式中:ed p_1为电网电压经正序dq变换得到的d轴分量、eq p_1为电网电压经正序dq变换得到的q轴分量、ed n_1为电网电压经负序dq变换得到的d轴分量、eq n_1为电网电压经负序dq变换得到的q轴分量、ed p为电网电压正序电压d轴分量、eq p为电网电压正序电压q轴分量、ed n为电网电压负序电压d轴分量、eq n为电网电压负序电压q轴分量;
步骤5,建立双dq坐标系下锁相环模型,包括:
Figure FDA0003629000170000121
式中:ω为锁相环输出的同步角频率、ω0为额定角频率、Kppll为锁相环比例系数、Kipll为锁相环积分系数、int ucq p_L为锁相环控制的积分项输出;
步骤6,建立双dq坐标系下各模组中间直流侧模型,包括建立双dq坐标系下中间直流侧模型方程、网侧输出功率方程,具体过程及表达式如下:
步骤6.1,建立双dq坐标系下中间直流侧模型方程,其表达式为:
Figure FDA0003629000170000122
Figure FDA0003629000170000123
式中:Pdc11为模组K21向网侧输出的功率、Pdc12为模组K11向中间直流侧输入的功率、idc1d为模组K11向直流侧输入的电流d轴分量、Pdc21为模组K22向网侧输出的功率、Pdc22为模组K12向中间直流侧输入的功率、idc2d为模组K12向中间直流侧输入电流d轴分量;
步骤6.2,建立网侧输出功率方程,其表达式为:
模组K21:
Figure FDA0003629000170000124
其中:
Figure FDA0003629000170000131
式中:Pdc11o为模组K21输出功率的平均值、Pdc11c2为模组K21输出功率余弦量2次脉动的幅值、Pdc11s2为模组K21输出功率正弦量2次脉动的幅值;
模组K22:
Figure FDA0003629000170000132
其中:
Figure FDA0003629000170000133
式中:Pdc21o为模组K22输出功率的平均值、Pdc21c2为模组K22输出功率余弦量2次脉动的幅值、Pdc21s2为模组K22输出功率正弦量2次脉动的幅值;
步骤7,建立单dq坐标系下机侧变流器、永磁同步电机及风力机等效模型:
将单dq坐标系下的机侧变流器、永磁同步电机及风力机等效成一个已知电流源模型,则:
Figure FDA0003629000170000134
式中:idc1d为模组K11机侧输出电流的d轴分量,idc1q为模组K11机侧输出电流的q轴分量,idc2d为模组K12机侧输出电流的d轴分量,idc2q为模组K12机侧输出电流的q轴分量。
CN202011056354.XA 2020-09-28 2020-09-28 基于双同步旋转坐标系的全功率风力发电系统的建模方法 Active CN112350364B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011056354.XA CN112350364B (zh) 2020-09-28 2020-09-28 基于双同步旋转坐标系的全功率风力发电系统的建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011056354.XA CN112350364B (zh) 2020-09-28 2020-09-28 基于双同步旋转坐标系的全功率风力发电系统的建模方法

Publications (2)

Publication Number Publication Date
CN112350364A CN112350364A (zh) 2021-02-09
CN112350364B true CN112350364B (zh) 2022-07-01

Family

ID=74361428

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011056354.XA Active CN112350364B (zh) 2020-09-28 2020-09-28 基于双同步旋转坐标系的全功率风力发电系统的建模方法

Country Status (1)

Country Link
CN (1) CN112350364B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114421517B (zh) * 2021-11-25 2022-12-13 广州鼎汉轨道交通装备有限公司 锁相环系统

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104253435A (zh) * 2013-06-27 2014-12-31 无锡乐华自动化科技有限公司 有源电力滤波器基于双dq变换的控制方法
DK3116087T3 (da) * 2015-07-07 2019-09-09 Siemens Gamesa Renewable Energy As Vindmølle, der er forbundet med et forsyningsnet via en HVDC-strømforbindelse gennem en netværksbrostyreenhed med effekt- og spændingsstyring
CN106877403A (zh) * 2015-12-13 2017-06-20 姚秋丽 一种永磁风力发电系统电网同步化方法
CN110994668A (zh) * 2019-10-29 2020-04-10 国网甘肃省电力公司电力科学研究院 基于双馈风电场并网系统输出阻抗模型的稳定性分析方法
CN110912186B (zh) * 2019-11-27 2021-04-06 国网冀北电力有限公司电力科学研究院 双馈型风力发电虚拟同步发电机建模方法

Also Published As

Publication number Publication date
CN112350364A (zh) 2021-02-09

Similar Documents

Publication Publication Date Title
CN109449958B (zh) 一种双馈风机并网系统稳定性分析方法
CN106451532B (zh) 一种多相直驱永磁风力发电变流一体化系统及其控制方法
CN105897013B (zh) 一种双向ac/dc变换器的虚拟惯性控制方法
CN108418253B (zh) 电流控制型虚拟同步发电机的阻抗建模与稳定性分析方法
CN102611138B (zh) 一种无延时的单相光伏并网功率调节方法
CN108429284B (zh) 基于电压源输出双馈风电机组的谐波电压比例前馈补偿方法
CN110912186B (zh) 双馈型风力发电虚拟同步发电机建模方法
CN101345423B (zh) 用于风力发电系统的五电平h桥级联背靠背变流器
CN110601268B (zh) 一种双馈风机并网端口输出阻抗建模及稳定性分析方法
CN101141110A (zh) 变速恒频双馈异步风力发电机转子电流无延时控制方法
CN112165115B (zh) 一种直驱风电机组灰箱模型的参数辨识方法及装置
CN108281986B (zh) 电压控制型虚拟同步发电机的阻抗建模与稳定性分析方法
CN105337312B (zh) 一种小功率风电变流器组网时的实时监控方法
CN105743091B (zh) 一种有源电力滤波器的双环解耦控制方法
CN103580055A (zh) 变速恒频双馈风力发电机组开放式并网实验系统及方法
CN112234654B (zh) 基于电容电压微分反馈的虚拟同步机的lc谐振抑制方法
CN109444541A (zh) 无变压器耦合注入式宽频带阻抗测量装置及其控制方法
CN104579060A (zh) 笼型转子无刷双馈风力发电机的间接功率控制方法
CN110336318A (zh) 一种单相光伏并网发电系统及控制方法
CN110350547A (zh) 一种适用于不平衡负载的多电平无功补偿系统及其控制方法
CN113346559A (zh) 极弱电网下直驱风电系统低电压穿越功率切换控制方法
CN112350364B (zh) 基于双同步旋转坐标系的全功率风力发电系统的建模方法
CN115036948A (zh) 一种双馈风机和牵引供电系统的负序协同补偿方法
CN113378347B (zh) 一种基于模块化多端口的风电机组频域阻抗建模方法
CN103366053B (zh) 一种电压定向矢量控制策略的改进及数学建模方法

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