CN112165115B - 一种直驱风电机组灰箱模型的参数辨识方法及装置 - Google Patents
一种直驱风电机组灰箱模型的参数辨识方法及装置 Download PDFInfo
- Publication number
- CN112165115B CN112165115B CN202010971256.2A CN202010971256A CN112165115B CN 112165115 B CN112165115 B CN 112165115B CN 202010971256 A CN202010971256 A CN 202010971256A CN 112165115 B CN112165115 B CN 112165115B
- Authority
- CN
- China
- Prior art keywords
- side converter
- impedance
- direct
- grid
- port
- 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 78
- 238000012546 transfer Methods 0.000 claims description 15
- 238000010408 sweeping Methods 0.000 claims description 12
- 230000001360 synchronised effect Effects 0.000 claims description 9
- 238000013461 design Methods 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 230000006855 networking Effects 0.000 abstract description 2
- 238000011160 research Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 239000003990 capacitor Substances 0.000 description 2
- 238000011217 control strategy Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006880 cross-coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- 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
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- 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]
-
- 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
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/28—The renewable source being wind energy
-
- 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
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B10/00—Integration of renewable energy sources in buildings
- Y02B10/30—Wind power
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
本发明提供了一种直驱风电机组灰箱模型的参数辨识方法及装置,首先,通过建立机侧系统(包括发电机和机侧变流器)直流端口的理论阻抗和测量阻抗,利用最小二乘法辨识得到机侧变流器的控制器参数。然后,建立直驱发电机组网侧交流端口的整体理论阻抗,并将辨识得到的机侧变流器控制器参数代入,通过与测量阻抗的比较,利用最小二乘法辨识得到网侧变流器的控制器参数。本发明利用交、直流端口阻抗频率特性实现直驱风电机组机侧变流器和网侧变流器不同带宽控制环节控制参数的精确辨识,且能够保证直驱风电机组端口宽频阻抗特性的一致性,适用于具有多带宽控制环节的风电机组参数辨识,具有实用性强、实现方便、辨识精确等优点。
Description
技术领域
本发明涉及风力发电技术领域,具体地说,涉及一种直驱风电机组灰箱模型的参数辨识方法及装置
背景技术
在能源资源短缺以及环境污染不断严重的背景下,风能作为可再生能源重要的一部分越来越引起人们的重视。目前,我国风电总装机容量世界第一,已成为我国第三大主力电源。随着陆上风电开发日趋饱和,海上风电将成为未来风力发电的主要发展方向。业内普遍认为,陆上风电机组将以双馈机型为主,海上风电机组将以直驱全功率机型为主。
直驱风电机组采用背靠背全功率变流器作为发电机与电网的接口装备,对风电机组的动态特性起决定性影响。直驱风电机组包含机侧变流器和网侧变流器,机、网侧变流器均含有复杂的控制结构。为研究风电机组/风电场并网的动态特性,往往需要建立风电机组的详细模型,包括仿真模型和数学模型,因此需要知道风电机组所采用的控制策略和控制参数。然而,风机制造商因商业保密,一般不会公开其所采用的具体控制策略和控制参数,即存在黑/灰箱问题,这给风电机组的准确建模提出了挑战。为解决该问题,相关学者开展了风电机组的参数辨识研究。目前,已有研究多集中在电气参数辨识方面,如发电机的参数辨识等,部分文献研究了风电变流器控制参数的辨识,但多采用时域辨识方法,仅能较准确地辨识单一参数,而对具有多带宽控制环节的控制系统参数辨识精度不高。此外,少数文献提出控制器参数的频域辨识方法,但所提方法基于控制器的解耦模型,导致控制器内环参数的辨识精度偏低。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种直驱风电机组灰箱模型的参数辨识方法及装置。
为实现上述目的,本发明采用以下技术方案。
根据本发明的一个方面,提供了一种直驱风电机组灰箱模型的参数辨识方法,包括以下步骤:
基于直驱风电机组的发电机与机侧变流器及其控制,建立机侧变流器直流端口理论阻抗;
在机侧变流器直流端口注入扰动电压,测量机侧变流器直流端口阻抗;
给定机侧变流器的控制器参数初始值,构建损失函数,使得测量得到的直流端口阻抗与直流端口理论阻抗在各个频率点处的阻抗差值的平方和最小,通过非线性最小二乘算法辨识获得机侧变流器的控制器参数;
建立直驱风电机组的网侧变流器交流端口整体理论正序阻抗,并将辨识得到的机侧变流器控制器参数代入此正序阻抗中;
在网侧变流器交流端口注入三相正序扰动电压,测量网侧变流器交流端口正序阻抗;
给定网侧变流器的控制器参数初始值,构建损失函数,使得测量得到的交流端口正序阻抗与代入了机侧变流器控制器参数的交流端口整体理论正序阻抗在各个频率点处的阻抗差值的平方和最小,通过非线性最小二乘算法辨识获得网侧变流器的控制器参数。
优选地,所述机侧变流器的控制器参数包括:转矩外环的比例和积分系数以及电流内环的比例和积分系数。
优选地,所述网侧变流器的控制器参数包括:直流电压外环的比例和积分系数、无功功率外环的比例和积分系数、电流内环的比例和积分系数以及锁相环的比例和积分系数。
优选地,所述建立机侧变流器直流端口理论阻抗模型的方法,包括:
式中,Rs为定子电阻;Lsd,Lsq分别为同步发电机定子d轴和q轴自感;ωr0为发电机稳态时的电气转速;
建立机侧变流器的交、直流三端口模型为:
式中,Htor(s),Hc(s)分别为转矩外环和电流内环的传递函数,ψm为体磁链,np为极对数,Ucd0,Ucq0,Icd0,Icq0分别为逆变器出口dq轴稳态电压与电流,Pc0为发电机输出的有功功率,Vdc0为直流端口稳态电压。
根据同步发电机和机侧变流器端口互联的电压电流约束条件,得到机侧变流器直流端口阻抗Zdc(s)为:
Zdc(s)=1/Ydc(s) (4)
优选地,所述在机侧变流器直流端口注入扰动电压,测量机侧变流器直流端口阻抗的方法,包括:
在机侧变流器直流端口注入谐波扰动电压,通过FFT分析提取谐波电压Δudc和谐波电流Δidc,计算得到机侧变流器直流端口阻抗Zdc为:
优选地,所述获得机侧变流器的控制器参数的方法,包括:
在直驱风电机组的系统电气参数已知的情况下,代入控制器参数初始值,计算得到直驱风电机组直流侧阻抗Zdc(f1,f2,...,fn),其中,f1,f2,...,fn代表n个不同的频率;
如果损失函数大于给定的极小值ξ,则根据设定的步长更新初始值为求得此更新初始值下的直驱风电机组机侧变流器直流侧阻抗Z′dc(f1,f2,...,fn),以此为基础进一步更新此初始值下损失函数的大小,重复此步骤,得到使损失函数最小下的控制器参数即为最终获得的机侧变流器的控制器参数。
优选地,所述建立网侧变流器交流端口整体理论正序阻抗模型的方法,包括:
建立网侧变流器交、直流三端口模型为:
其中:
式中,Ucd0,Ucq0,Icd0,Icq0分别为逆变器出口dq轴稳态电压与电流,Hpq(s)是无功功率外环的传递函数,Hdc(s)是直流电压外环的传递函数,Hc为电流内环的传递函数,Vdc0是网侧变流器直流端口稳态电压Hpll(s)是PLL的传递函数,Upccd0为PCC点的d轴稳态电压,Rf,Lf是滤波器的电阻和电感;
式中,Ccap是直流侧电容值;
式中:
其中,其中j为虚数单位。
将辨识得到的机侧变流器控制器参数代入上述得到的正序阻抗,其方法为:
优选地,所述在网侧变流器交流端口注入三相正序扰动电压,测量网侧变流器交流端口正序阻抗的方法,包括:
在网侧变流器交流端口注入三相扰动电压,通过FFT分析提取正序扰动电压分量Δup和正序扰动电流分量Δip,计算得到网侧变流器交流端口正序阻抗Zp为:
优选地,所述获得网侧变流器的控制器参数的方法,包括:
根据风电机组控制器参数一般设计方法,结合获得的机侧变流器控制器参数,计算得到网侧变流器控制器参数,作为网侧变流器控制器参数的初始值其中,分别为电流内环的比例、积分常数,为直流电压外环的比例、积分常数,分别为功率外环的比例、积分常数,为PLL的比例、积分常数;
在直驱风电机组的系统电气参数已知的情况下,代入控制器参数初始值,计算得到直驱风电机组网侧变流器交流侧正序阻抗Zp(f1,f2,...,fn),其中,f1,f2,...,fn代表n个不同的频率;
如果损失函数大于给定的极小值ξ,则根据设定的步长更新初始值为求得此更新初始值下的直驱风电机组网侧变流器交流侧正序阻抗Z′p(f1,f2,...,fn),以此为基础进一步更新此初始值下损失函数的大小,重复此步骤,得到使损失函数最小下的控制器参数即为最终获得的网侧变流器的控制器参数。
根据本发明的另一个方面,提供了一种直驱风电机组灰箱模型的参数辨识装置,包括:直流端口阻抗建立模块、直流端口阻抗扫频模块、机侧变流器控制器参数辨识模块、交流端口阻抗建立模块、交流端口扫频模块、网侧变流器控制器参数辨识模块;其中:
所述直流端口阻抗建立模块,基于直驱风电机组的发电机与机侧变流器及其控制,用于建立机侧变流器直流端口的理论阻抗;
所述直流端口阻抗扫频模块,在机侧变流器直流端口注入扰动电压,用于测量机侧变流器直流端口阻抗;
所述机侧变流器控制器参数辨识模块,采用非线性最小二乘法,将直流端口理论阻抗与测量阻抗做差,使得差的平方和最小,用于辨识获得机侧变流器控制器参数;
所述交流端口阻抗建立模块,代入用于建立网侧变流器交流端口的理论正序阻抗,将辨识得到的机侧变流器控制器参数代入此正序阻抗;
所述交流端口扫频模块,在网侧变流器交流端口注入三相正序扰动电压,用于测量网侧变流器交流端口的正序阻抗;
所述网侧变流器控制器参数辨识模块,采用非线性最小二乘法,将交流端口理论正序阻抗与测量正序阻抗做差,使得差的平方和最小,用于辨识获得网侧变流器控制器参数。
根据本发明的第三个方面,提供了一种直驱风电机组,采用上述任一项所述的参数辨识方法对其灰箱参数进行辨识,获得机侧变流器的控制器参数和网侧变流器的控制器参数。
由于采用了上述技术方案,本发明与现有技术相比,具有如下的有益效果:
本发明提供的直驱风电机组灰箱模型的参数辨识方法及装置,基于风电机组端口测量阻抗数据进行辨识,具有较强的工程实用价值,且算法简单,易于实现。
本发明提供的直驱风电机组灰箱模型的参数辨识方法及装置,所提辨识方法适用于具有不同带宽控制环节的控制系统参数辨识,且辨识精度高,所辨识参数能够保证系统端口宽频阻抗特性的一致性,为风电并网的宽频振荡问题研究提供准确的模型基础。
实施本发明的任一产品并不一定需要同时达到以上所述的所有优点。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明一优选实施例中永磁直驱风电机组的拓扑结构示意图;
图2为本发明一优选实施例中机侧变流器控制系统框图;
图3为本发明一优选实施例中网侧变流器控制系统框图;
图4为本发明一优选实施例中直驱风电机组灰箱模型的参数辨识方法的流程图;
图5为本发明一优选实施例中直驱风电机组灰箱模型的参数辨识装置结构示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进。这些都属于本发明的保护范围。
本发明一实施例提供了一种直驱风电机组灰箱模型的参数辨识方法,该方法虑发电机与机侧变流器及其控制,建立机侧变流器的直流端口理论阻抗模型;然后,在永磁直驱风电机组机侧变流器直流侧端口注入小扰动电压,测得直流端口的阻抗;根据风电机组控制器参数一般设计方法,给定机侧变流器控制器参数合适的初始值,并构建损失函数,使直流端口测量阻抗和理论阻抗在各频率点的阻抗差值的平方和最小,使用非线性最小二乘法寻优,最终辨识出机侧变流器控制器参数。接着,将辨识得到的机侧变流器控制器参数代入,建立永磁直驱风电机组交流端口整体dq阻抗理论模型,并利用dq阻抗与正序阻抗之间的转换关系,将dq阻抗转换成正序阻抗;在网侧变流器交流测端口注入三相小扰动电压,测得交流端口的正序阻抗;根据风电机组控制器参数一般设计方法,给定网侧变流器控制器参数合适的初值,并构建损失函数,使交流端口测量正序阻抗与理论正序阻抗在各个频率点的阻抗差值的平方和最小,使用非线性最小二乘法寻优,最终辨识出网侧变流器控制器参数。
具体地,本实施例提供的方法,包括以下步骤:
步骤S1:考虑发电机与机侧变流器及其控制(即转矩外环和电流内环控制),建立机侧变流器的直流端口理论阻抗模型。
步骤S2:在直驱风电机组机侧变流器直流端口注入小扰动电压,测量其直流端口阻抗。
步骤S3:给定机侧变流器控制器参数初始值,构建损失函数,使得机侧变流器直流端口测量阻抗与理论阻抗在各个频率点处的阻抗差值的平方和最小,通过非线性最小二乘算法辨识获得机侧变流器的控制器参数。
步骤S4:建立直驱风电机组的网侧变流器交流端口整体理论正序阻抗,并将辨识得到的机侧变流器控制器参数代入该正序阻抗中;
步骤S5:在直驱风电机组网侧变流器交流端口注入三相正序小扰动电压,测量其交流端口正序阻抗。
步骤S6:给定网侧变流器控制器参数初始值,构建损失函数,使得交流端口测量正序阻抗与代入了机侧变流器控制器参数的理论正序阻抗在各个频率点处的阻抗差值的平方和最小,通过非线性最小二乘算法辨识获得网侧变流器的控制器参数。
作为一优选实施例,所述机侧变流器的控制器参数包括:转矩外环的比例和积分系数以及电流内环的比例和积分系数。
作为一优选实施例,所述网侧变流器的控制器参数包括:直流电压外环的比例和积分系数、无功功率外环的比例和积分系数、电流内环的比例和积分系数以及锁相环的比例和积分系数。
作为一优选实施例,所述建立机侧变流器直流端口理论阻抗模型的方法,包括:
式中Rs为定子电阻;Lsd,Lsq分别为同步发电机定子d轴和q轴自感;ωr0为发电机稳态时的电气转速;
建立机侧变流器的交、直流三端口模型为:
式中,Htor(s),Hc(s)分别为转矩外环和电流内环的传递函数,ψm为体磁链,np为极对数,Ucd0,Ucq0,Icd0,Icq0分别为逆变器出口dq轴稳态电压与电流,Pc0为发电机输出的有功功率,Vdc0为直流端口稳态电压;
根据同步发电机和机侧变流器端口互联的电压电流约束条件,得到机侧变流器直流端口阻抗Zdc(s)为:
Zdc(s)=1/Ydc(s) (4)
作为一优选实施例,所述在机侧变流器直流端口注入扰动电压,测量机侧变流器直流端口阻抗的方法,包括:
在机侧变流器直流端口注入谐波扰动电压,通过FFT分析提取谐波电压Δudc和谐波电流Δidc,计算得到机侧变流器直流端口阻抗Zdc为:
作为一优选实施例,所述获得机侧变流器的控制器参数的方法,包括:
根据风电机组控制器参数一般设计方法计算得到机侧变流器控制器参数,作为机侧变流器控制器参数的初始值其中分别为机侧变流器电流内环的比例、积分常数和转矩外环的比例积分常数;在直驱风电机组的系统电气参数已知的情况下,代入控制器参数初始值,计算得到直驱风电机组直流侧阻抗Zdc(f1,f2,...,fn),其中f1,f2,...,fn代表n个不同的频率;
如果损失函数大于给定的极小值ξ(即为一设定的阈值),则根据设定的步长更新初始值为求得此更新初始值下的直驱风电机组机侧变流器直流侧阻抗Z′dc(f1,f2,...,fn),以此为基础进一步更新此初始值下损失函数的大小,重复此步骤,得到使损失函数最小下的控制器参数即为最终获得的机侧变流器的控制器参数。
作为一优选实施例,所述建立网侧变流器交流端口整体理论正序阻抗模型的方法,包括:
建立网侧变流器交、直流三端口模型为:
其中:
式中,Ucd0,Ucq0,Icd0,Icq0分别为逆变器出口dq轴稳态电压与电流,Hpq(s)是无功功率外环的传递函数,Hdc(s)是直流电压外环的传递函数,Hc为电流内环的传递函数,Vdc0为直流端口稳态电压,Hpll(s)是PLL的传递函数,Upccd0为并网点的d轴稳态电压,Rf,Lf是滤波器的电阻和电感;
式中,Ccap是直流侧电容值;
式中:
其中,Vdc0为直流端口稳态电压;
将辨识得到的机侧变流器控制器参数代入上述得到的正序阻抗,其方法为:
作为一优选实施例,所述在网侧变流器交流端口注入三相正序扰动电压,测量网侧变流器交流端口正序阻抗的方法,包括:
在网侧变流器交流端口注入三相扰动电压,通过FFT分析提取正序扰动电压分量Δup和正序扰动电流分量Δip,计算得到网侧变流器交流端口正序阻抗Zp为:
作为一优选实施例,所述获得网侧变流器的控制器参数的方法,包括:
根据风电机组控制器参数一般设计方法,结合获得的机侧变流器控制器参数,计算得到网侧变流器控制器参数,作为网侧变流器控制器参数的初始值其中分别为电流内环的比例、积分常数,为直流电压外环的比例、积分常数,分别为功率外环的比例、积分常数,为PLL的比例、积分常数;在直驱风电机组的系统电气参数已知的情况下,代入控制器参数初始值,计算得到直驱风电机组网侧变流器交流侧正序阻抗Zp(f1,f2,...,fn),其中f1,f2,...,fn代表n个不同的频率;
如果损失函数大于给定的极小值ξ(即为一设定的阈值),则根据设定的步长更新初始值为求得此更新初始值下的直驱风电机组网侧变流器交流侧正序阻抗Z′p(f1,f2,...,fn),以此为基础进一步更新此初始值下损失函数的大小,重复此步骤,得到使损失函数最小下的控制器参数即为最终获得的网侧变流器的控制器参数。
本发明另一实施例提供了一种直驱风电机组灰箱模型的参数辨识装置,用于实现本发明上述实施例所提供的直驱风电机组灰箱模型的参数辨识方法,包括:直流端口阻抗建立模块、直流端口阻抗扫频模块、机侧变流器控制器参数辨识模块、交流端口阻抗建立模块、交流端口扫频模块、网侧变流器控制器参数辨识模块;其中,
所述直流端口阻抗建立模块,基于直驱风电机组的发电机与机侧变流器及其控制(即转矩外环和电流内环控制),用于建立机侧变流器直流端口的理论阻抗,需考虑其具体的控制结构;
所述直流端口阻抗扫频模块,在机侧变流器直流端口注入扰动电压,用于测量机侧变流器直流端口阻抗;
所述机侧变流器控制器参数辨识模块采用非线性最小二乘法,将直流端口理论阻抗与测量阻抗做差,使得差的平方和最小,可辨识出机侧变流器控制器参数;
所述交流端口阻抗建立模块,用于建立网侧变流器交流端口的理论正序阻抗,并将辨识得到的机侧变流器控制器参数代入此正序阻抗,需考虑其具体的控制结构;
所述交流端口扫频模块,在网侧变流器交流端口注入三相正序扰动电压,用于测量网侧变流器交流端口的正序阻抗。
所述网侧变流器控制器参数辨识模块采用非线性最小二乘法,将交流端口理论正序阻抗与测量正序阻抗做差,使得差的平方和最小,可辨识出网侧变流器控制器参数。
本发明第三个实施例提供了一种直驱风电机组,采用本发明上述实施例中任一项所述的参数辨识方法对其灰箱参数进行辨识,获得机侧变流器的控制器参数和网侧变流器的控制器参数。
下面结合附图以及具体应用实例,对本发明上述实施例所提供的技术方案进一步详细描述如下。
具体实例1:
本具体实例中,基于永磁直驱风电机组机,交流电网电压为690V/50Hz,风电场PCC点电压为690V/50Hz,永磁同步发电机,额定功率为2MW。机侧变流器控制系统转矩给定值为-1.6e6N,网侧变流器控制系统直流电压给定值为1120V。
如图1所示,本具体实例中,永磁直驱风场由以下几个部分构成,包括:永磁直驱风力发电机,机侧变流器,直流侧电容,网侧变流器,滤波电感以及交流电网。
如图2所示,控制结构为转矩外环、电流内环的双闭环控制结构,外环控制转矩,Tor*是转矩参考值,Tor是转矩反馈值,两者之差经外环PI调节器(PI_Tor模块)送入机侧变流器电流内环,并作为电流内环q轴的给定值电流内环q轴给定值与反馈值之差并加上发电机定子q轴电压作为SPWM的q轴驱动信号;机侧变流器电流内环d轴参考值设置为零,与d轴电流反馈值之差并加上发电机定子d轴电压作为SPWM的d轴驱动信号。SPWM的d轴与q轴驱动信号通过转子位置角θr的dq/abc的坐标变换为三相驱动信号,驱动IGBT的开通。
如图3(a)和(b)所示,(a)中控制结构为直流电压外环/无功功率外环、电流内环的双闭环控制结构,d轴外环控制直流端口的电压,Udc*是直流电压参考值,Udc是直流电压反馈值,两者之差经外环PI调节器(PI_Udc模块)送入网侧变流器电流内环,并作为电流内环d轴的给定值电流内环d轴给定值与反馈值之差并加上电网电压前馈量和交叉耦合分量作为SPWM的d轴驱动信号,其中为PCC点的d轴电压分量,ω为同步角频率,Lf为滤波电感;q轴外环控制逆变器输出的的无功功率,Q*是无功功率参考值,Q是无功功率反馈值,两者之差经外环PI调节器(PI_PQ模块)送入网侧变流器电流内环,并作为电流内环q轴的给定值电流内环q轴给定值与反馈值之差并加上电压前馈量和交叉耦合分量作为SPWM的q轴驱动信号,其中为PCC点的q轴电压分量;(b)中给出了电网电压定向的原理,是电网三相电压,经过abc/dq变换得到电网电压的dq轴分量和经过锁相环PI调节器PI_PLL后加上电网角频率w并经过积分环节即可得到电网位置角θg。SPWM的dq轴驱动信号即可经过电网位置角θg的dq/abc的坐标变换变换为三相驱动信号,驱动IGBT的开通。如图4所示,本具体实例中,所述的一种直驱风电机组灰箱模型的参数辨识方法,包括如下步骤:
S1:考虑发电机与机侧变流器及其控制,建立机侧变流器的直流端口理论阻抗模型。
S2:在永磁直驱风电机组机侧变流器直流端口注入小扰动电压,测量其直流端口阻抗。
S3:给定机侧变流器控制器参数初始值,构建损失函数,使得机侧变流器直流端口测量阻抗与理论阻抗在各个频率点处的阻抗差值的平方和最小,通过非线性最小二乘算法辨识获得机侧变流器的控制器参数。
S4:建立永磁直驱风电机组交流端口整体理论阻抗模型,并将辨识得到的机侧变流器控制器参数代入该阻抗模型。
S5:在永磁直驱风电机组交流端口注入三相正序小扰动电压,测量其交流端口正序阻抗。
S6:给定网侧变流器控制器参数初始值,构建损失函数,使得交流端口测量正序阻抗与理论正序阻抗在各个频率点处的阻抗差值的平方和最小,通过非线性最小二乘算法辨识获得网侧变流器的控制器参数。
如图5所示,本具体实例中,直驱风电机组灰箱模型的参数辨识装置,包括如下模块:
直流端口阻抗建立模块:用于建立机侧变流器直流端口的理论阻抗,需考虑其具体的控制结构;
直流端口阻抗扫频模块:用于测量机侧变流器直流端口阻抗;
机侧变流器控制器参数辨识模块:采用非线性最小二乘法,将直流端口理论阻抗与测量阻抗做差,使得差的平方和最小,可辨识出机侧变流器控制器参数;
交流端口阻抗建立模块:用于建立网侧变流器交流端口的理论正序阻抗,需考虑其具体的控制结构;
交流端口扫频模块:用于测量网侧变流器交流端口的正序阻抗;
网侧变流器控制器参数辨识模块:采用非线性最小二乘法,将交流端口理论正序阻抗与测量正序阻抗做差,使得差的平方和最小,可辨识出网侧变流器控制器参数。
如表1所示,在本具体实例中,对永磁直驱风电机组机侧变流器控制器参数进行参数辨识,给出了机侧变流器控制器参数的辨识结果,可以看到,电流内环控制器比例和积分系数的误差均在1%以内,而转矩外环控制器比例和积分系数的误差较大,原因主要有两个:(1)转矩外环控制器的比例和积分系数太小,导致辨识出的结果相对误差较大;(2)转矩外环对阻抗特性的灵敏度较低,导致参数的可辨识性较差。但辨识出的结果最大偏差为23%,辨识精度在可接受范围内。
表1机侧变流器控制器参数辨识结果
如表2所示,在本具体实例中,对永磁直驱风电机组网侧变流器控制器参数进行参数辨识,给出了网侧变流器控制器参数的辨识结果,可以看到,辨识的8个控制器参数中,有一个最大误差为4.7%,其余均小于2%,辨识精度较高。
表2网侧变流器控制器参数辨识结果
由此可见,本发明上述实施例提供的直驱风电机组灰箱模型的参数辨识方法及装置,通过端口阻抗特性分步辨识机侧变流器和网侧变流器的控制器参数,使得端口理论阻抗与测量阻抗的频域特性完全一致,实现了直驱风电机组不同控制环节控制器参数的精确辨识,辨识过程操作简便、算法简单。
本发明上述实施例提供的直驱风电机组灰箱模型的参数辨识方法及装置,考虑了风电机组内部详细拓扑及各个组成元件,包括永磁同步发电机、机侧变流器及其控制、直流侧电容、网侧变流器及其控制、网侧滤波器等。首先,通过建立机侧系统(包括发电机和机侧变流器)直流端口的理论阻抗和测量阻抗,利用最小二乘法辨识得到机侧变流器的控制器参数。然后,将辨识得到的机侧变流器控制器参数代入,通过建立直驱发电机组网侧交流端口的整体理论阻抗和测量阻抗,利用最小二乘法辨识得到网侧变流器的控制器参数。本发明上述实施例提供的方法及装置,利用交、直流端口阻抗频率特性实现直驱风电机组机侧变流器和网侧变流器不同带宽控制环节控制参数的精确辨识,且能够保证直驱风电机组端口宽频阻抗特性的一致性,适用于具有多带宽控制环节的风电机组参数辨识,具有实用性强、实现方便、辨识精确等优点。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变形或修改,这并不影响本发明的实质内容。
Claims (10)
1.一种直驱风电机组灰箱模型的参数辨识方法,其特征在于,包括以下步骤:
基于直驱风电机组的发电机与机侧变流器及其控制,建立机侧变流器直流端口理论阻抗;
在机侧变流器直流端口注入扰动电压,测量机侧变流器直流端口阻抗;
给定机侧变流器的控制器参数初始值,构建损失函数,使得测量得到的直流端口阻抗与直流端口理论阻抗在各个频率点处的阻抗差值的平方和最小,通过非线性最小二乘算法辨识获得机侧变流器的控制器参数;
建立直驱风电机组的网侧变流器交流端口整体理论正序阻抗,并将辨识得到的机侧变流器控制器参数代入此正序阻抗中;
在网侧变流器交流端口注入三相正序扰动电压,测量网侧变流器交流端口正序阻抗;
给定网侧变流器的控制器参数初始值,构建损失函数,使得测量得到的交流端口正序阻抗与代入了机侧变流器控制器参数的交流端口整体理论正序阻抗在各个频率点处的阻抗差值的平方和最小,通过非线性最小二乘算法辨识获得网侧变流器的控制器参数。
2.根据权利要求1所述的一种直驱风电机组灰箱模型的参数辨识方法,其特征在于:所述机侧变流器的控制器参数包括:
转矩外环的比例和积分系数以及电流内环的比例和积分系数;
所述网侧变流器的控制器参数包括:
直流电压外环的比例和积分系数、无功功率外环的比例和积分系数、电流内环的比例和积分系数以及锁相环的比例和积分系数。
3.根据权利要求1所述的一种直驱风电机组灰箱模型的参数辨识方法,其特征在于,所述建立机侧变流器直流端口理论阻抗的方法,包括:
式中,Rs为定子电阻;Lsd,Lsq分别为同步发电机定子d轴和q轴自感;ωr0为发电机稳态时的电气转速;
建立机侧变流器的交、直流三端口模型为:
式中,Htor(s),Hc(s)分别为转矩外环和电流内环的传递函数,ψm为体磁链,np为极对数,Ucd0,Ucq0,Icd0,Icq0分别为逆变器出口dq轴稳态电压与电流,Pc0为发电机输出的有功功率,Vdc0为直流端口稳态电压;
根据同步发电机和机侧变流器端口互联的电压电流约束条件,得到机侧变流器直流端口阻抗Zdc(s)为:
Zdc(s)=1/Ydc(s) (4)
5.根据权利要求1所述的一种直驱风电机组灰箱模型的参数辨识方法,其特征在于,所述获得机侧变流器的控制器参数的方法,包括:
在直驱风电机组的系统电气参数已知的情况下,代入控制器参数初始值,计算得到直驱风电机组直流侧阻抗Zdc(f1,f2,...,fn),其中,f1,f2,...,fn代表n个不同的频率;
6.根据权利要求1所述的一种直驱风电机组灰箱模型的参数辨识方法,其特征在于,所述建立网侧变流器交流端口整体理论正序阻抗的方法,包括:
建立网侧变流器交、直流三端口模型为:
其中:
式中,Ucd0,Ucq0,Icd0,Icq0分别为逆变器出口dq轴稳态电压与电流,Hpq(s)是无功功率外环的传递函数,Hdc(s)是直流电压外环的传递函数,Hc为电流内环的传递函数,Vdc0是网侧变流器直流端口稳态电压,Hpll(s)是PLL的传递函数,Upccd0为并网点的d轴稳态电压,Rf,Lf是滤波器的电阻和电感;Pc0为发电机发出的有功功率;
式中,Ccap是直流侧电容值,Ydc(s)为机侧变流器直流端口导纳;
式中:
其中,j为虚数单位;
将辨识得到的机侧变流器控制器参数代入上述得到的正序阻抗,其方法为:
8.根据权利要求1所述的一种直驱风电机组灰箱模型的参数辨识方法,其特征在于,所述获得网侧变流器的控制器参数的方法,包括:
根据风电机组控制器参数设计方法,结合获得的机侧变流器控制器参数,计算得到网侧变流器控制器参数,作为网侧变流器控制器参数的初始值其中,分别为电流内环的比例、积分常数,为直流电压外环的比例、积分常数,分别为功率外环的比例、积分常数,为PLL的比例、积分常数;
在直驱风电机组的系统电气参数已知的情况下,代入控制器参数初始值,计算得到直驱风电机组网侧变流器交流侧正序阻抗Zp(f1,f2,...,fn),f1,f2,...,fn代表n个不同的频率;
9.一种直驱风电机组灰箱模型的参数辨识装置,其特征在于,包括:直流端口阻抗建立模块、直流端口阻抗扫频模块、机侧变流器控制器参数辨识模块、交流端口阻抗建立模块、交流端口扫频模块、网侧变流器控制器参数辨识模块;其中:
所述直流端口阻抗建立模块,基于直驱风电机组的发电机与机侧变流器及其控制,用于建立机侧变流器直流端口的理论阻抗;
所述直流端口阻抗扫频模块,在机侧变流器直流端口注入扰动电压,用于测量机侧变流器直流端口阻抗;
所述机侧变流器控制器参数辨识模块,采用非线性最小二乘法,将直流端口理论阻抗与测量阻抗做差,使得差的平方和最小,用于辨识获得机侧变流器控制器参数;
所述交流端口阻抗建立模块,代入用于建立网侧变流器交流端口的理论正序阻抗,并将辨识得到的机侧变流器控制器参数代入此正序阻抗;
所述交流端口扫频模块,在网侧变流器交流端口注入三相正序扰动电压,用于测量网侧变流器交流端口的正序阻抗;
所述网侧变流器控制器参数辨识模块,采用非线性最小二乘法,将交流端口理论正序阻抗与测量正序阻抗做差,使得差的平方和最小,用于辨识获得网侧变流器控制器参数。
10.一种直驱风电机组,其特征在于,采用权利要求1-8任一项所述的参数辨识方法对其灰箱参数进行辨识,获得机侧变流器的控制器参数和网侧变流器的控制器参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010971256.2A CN112165115B (zh) | 2020-09-16 | 2020-09-16 | 一种直驱风电机组灰箱模型的参数辨识方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010971256.2A CN112165115B (zh) | 2020-09-16 | 2020-09-16 | 一种直驱风电机组灰箱模型的参数辨识方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112165115A CN112165115A (zh) | 2021-01-01 |
CN112165115B true CN112165115B (zh) | 2022-08-02 |
Family
ID=73857954
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010971256.2A Active CN112165115B (zh) | 2020-09-16 | 2020-09-16 | 一种直驱风电机组灰箱模型的参数辨识方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112165115B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113659588A (zh) * | 2021-02-09 | 2021-11-16 | 内蒙古能源发电投资集团有限公司 | 一种弱电网下永磁直驱风电机组振荡抑制方法 |
CN112865181B (zh) * | 2021-03-02 | 2022-07-29 | 国网冀北电力有限公司电力科学研究院 | 基于端口阻抗特性的光伏逆变器参数辨识方法及装置 |
CN114123336B (zh) * | 2021-11-30 | 2023-06-16 | 国网山东省电力公司德州供电公司 | 一种直驱风机变流器的控制参数辨识方法及系统 |
CN115249980B (zh) * | 2022-08-18 | 2024-04-19 | 上海交通大学 | 数据和知识联合驱动的新能源场站阻抗辨识方法及系统 |
CN117277412B (zh) * | 2023-09-26 | 2024-02-27 | 中国电力科学研究院有限公司 | 风电机组变流器的控制参数识别方法、系统、设备和介质 |
CN117578596B (zh) * | 2024-01-16 | 2024-03-29 | 湖南大学 | 一种直驱永磁风力发电机组连续故障穿越控制方法及系统 |
CN118244184B (zh) * | 2024-05-28 | 2024-08-13 | 中国电力科学研究院有限公司 | 一种并网变流器的序阻抗校核方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103887815A (zh) * | 2014-02-21 | 2014-06-25 | 华南理工大学 | 基于运行数据的风电场参数辨识及动态等值方法 |
CN107328996A (zh) * | 2017-06-12 | 2017-11-07 | 华北电力大学 | 抗差输电线路正序参数辨识方法 |
CN108631358A (zh) * | 2018-03-30 | 2018-10-09 | 中国电力科学研究院有限公司 | 基于控制硬件在环的直驱风电机组阻抗确定方法和装置 |
CN109672221A (zh) * | 2019-02-26 | 2019-04-23 | 西南交通大学 | 一种用于次同步振荡分析的直驱风电场动态等值方法 |
CN110165956A (zh) * | 2019-05-14 | 2019-08-23 | 北京理工大学 | 一种开绕组永磁同步电机零序电感在线辨识方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8138631B2 (en) * | 2007-12-21 | 2012-03-20 | Eiq Energy, Inc. | Advanced renewable energy harvesting |
-
2020
- 2020-09-16 CN CN202010971256.2A patent/CN112165115B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103887815A (zh) * | 2014-02-21 | 2014-06-25 | 华南理工大学 | 基于运行数据的风电场参数辨识及动态等值方法 |
CN107328996A (zh) * | 2017-06-12 | 2017-11-07 | 华北电力大学 | 抗差输电线路正序参数辨识方法 |
CN108631358A (zh) * | 2018-03-30 | 2018-10-09 | 中国电力科学研究院有限公司 | 基于控制硬件在环的直驱风电机组阻抗确定方法和装置 |
CN109672221A (zh) * | 2019-02-26 | 2019-04-23 | 西南交通大学 | 一种用于次同步振荡分析的直驱风电场动态等值方法 |
CN110165956A (zh) * | 2019-05-14 | 2019-08-23 | 北京理工大学 | 一种开绕组永磁同步电机零序电感在线辨识方法 |
Non-Patent Citations (4)
Title |
---|
Online parameter estimator of the DC bus capacitor bank for Doubly-Fed Induction Generators;Cristian Blanco et al.;《2019 21st European Conference on Power Electronics and Applications (EPE "19 ECCE Europe)》;20191128;第1-9页 * |
基于控制硬件在环的风电机组阻抗测量及影响因素分析;李光辉等;《电网技术》;20190531;第43卷(第5期);第1624-1631页 * |
直驱永磁同步风力发电机的电气参数辨识;张仰飞等;《电力系统自动化》;20120725;第36卷(第14期);第150-153页 * |
风电机组电磁暂态建模及验证;张兴等;《中国电力》;20200731;第53卷(第7期);第106-112页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112165115A (zh) | 2021-01-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112165115B (zh) | 一种直驱风电机组灰箱模型的参数辨识方法及装置 | |
Xue et al. | A complete impedance model of a PMSG-based wind energy conversion system and its effect on the stability analysis of MMC-HVDC connected offshore wind farms | |
CN107732939B (zh) | 基于电压源型变换器解耦控制的抑制次同步振荡控制方法 | |
CN111525611A (zh) | 计及频率耦合效应的双馈并网系统次同步振荡分析方法 | |
Molina et al. | Dynamic modeling of wind farms with variable-speed direct-driven PMSG wind turbines | |
CN112615393B (zh) | 基于矢量拟合的直驱风电机组控制器参数辨识方法及装置 | |
CN110676874B (zh) | 计及频率耦合效应的直驱式风机次同步振荡电气量分析方法 | |
CN108281986B (zh) | 电压控制型虚拟同步发电机的阻抗建模与稳定性分析方法 | |
CN109921421B (zh) | 双馈风电机组输出谐波电流模型建立方法 | |
CN112234654B (zh) | 基于电容电压微分反馈的虚拟同步机的lc谐振抑制方法 | |
CN112800381B (zh) | 一种弱电网下pmsg并网系统的阻抗测量方法 | |
CN106452234B (zh) | 一种双馈风力发电机定子匝间短路负序电流抑制方法 | |
CN109586337A (zh) | 基于频域建模的vsc并网系统次同步振荡风险评估方法 | |
CN106451539B (zh) | 一种计及永磁直驱风电机组动态特性的风电场并网稳定分析方法 | |
CN111769597B (zh) | 一种双馈风力发电机的降维建模分析方法 | |
CN110212575A (zh) | 双馈风机单机无穷大系统的小干扰稳定性模态分析方法 | |
CN113378347B (zh) | 一种基于模块化多端口的风电机组频域阻抗建模方法 | |
CN106655247A (zh) | 一种LCCL型光伏并网逆变器dq轴解耦控制设计方法 | |
CN105552951B (zh) | 一种基于重复滑模的dfig系统控制方法 | |
CN110210170A (zh) | 一种大规模风电机群等值小信号模型建模方法 | |
CN112994113B (zh) | 一种双馈风力发电机序阻抗建模方法及系统 | |
CN112653184B (zh) | 风力发电装备黑箱模型的辨识方法、装置、终端及介质 | |
CN112260263B (zh) | 陆上双馈风场经传统直流送出系统的阻抗测量方法及装置 | |
CN114744640A (zh) | 一种分布式调相机序阻抗建模方法 | |
CN108448989A (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 |