CN114552652B - 基于参数敏感性因子的并网逆变器锁相环参数辨识方法 - Google Patents
基于参数敏感性因子的并网逆变器锁相环参数辨识方法 Download PDFInfo
- Publication number
- CN114552652B CN114552652B CN202210250369.2A CN202210250369A CN114552652B CN 114552652 B CN114552652 B CN 114552652B CN 202210250369 A CN202210250369 A CN 202210250369A CN 114552652 B CN114552652 B CN 114552652B
- Authority
- CN
- China
- Prior art keywords
- phase
- locked loop
- parameter
- coefficient
- frequency band
- 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 66
- 230000035945 sensitivity Effects 0.000 title claims abstract description 66
- 239000002245 particle Substances 0.000 claims abstract description 156
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 66
- 238000005457 optimization Methods 0.000 claims abstract description 53
- 238000002847 impedance measurement Methods 0.000 claims abstract description 48
- 238000005070 sampling Methods 0.000 claims description 37
- 230000008707 rearrangement Effects 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 17
- 230000006870 function Effects 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 5
- 238000010168 coupling process Methods 0.000 claims description 5
- 238000005859 coupling reaction Methods 0.000 claims description 5
- 230000010354 integration Effects 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 8
- 238000004088 simulation Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000003337 fertilizer Substances 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
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
-
- 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/24—Arrangements for preventing or reducing oscillations of power in networks
- H02J3/241—The oscillation concerning frequency
-
- 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
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/40—Arrangements for reducing harmonics
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Measurement Of Resistance Or Impedance (AREA)
- Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
Abstract
本发明提供了一种基于参数敏感性因子的并网逆变器锁相环参数辨识方法，属于电能质量和控制领域。所述锁相环参数包括锁相环比例系数和锁相环积分系数。该方法建立含有未知锁相环参数的正序阻抗模型，提出用参数敏感性因子得到参数敏感性因子满足要求的频率段，使用粒子群优化算法结合在参数敏感性因子满足要求的频率段内测量得到的阻抗进行参数估计，获得多组锁相环参数估计值，确定锁相环参数估计值的平滑频段，获得参数辨识结果。本发明综合了并网逆变器序阻抗建模、参数敏感性因子、阻抗测量、粒子群优化算法和平滑频段确定方法的优点，缩短了阻抗测量的频率范围，能够准确地辨识出并网逆变器锁相环参数。
Description
技术领域
本发明属于电能质量和控制领域，具体涉及一种并网条件下在并网逆变器锁相环参数未知的情况时基于参数敏感性因子进行锁相环参数辨识的方法。
背景技术
基于阻抗模型的并网系统稳定性分析，需要对并网逆变器进行阻抗建模获得其阻抗模型。而阻抗模型的建立需要获取并网逆变器的控制参数，如电流环参数、锁相环参数等。参数辨识方法被广泛用于在逆变器控制参数未知的情况下辨识出逆变器控制参数。
参数辨识方法需要并网逆变器的测量阻抗数据，然而由于噪声和谐波损耗等的影响，测量阻抗与理论阻抗之间存在误差，这导致参数辨识结果的不准确。例如：
文献1“宋依桐，胡海涛，肖冬华，潘鹏宇，何正友.牵引供电系统注入式谐波阻抗测量的扰动限值选定方法[J].中国电机工程学报，2020，40(15)：4770-4780.”考虑噪声、背景谐波等因素进行指标形式的推导，能够在面对已知对象、工况或部分基本参数的条件下，在实施阻抗测量前进行明确的注入扰动限值判定，使阻抗测量准确。
在并网逆变器阻抗建模中，锁相环参数的影响十分重要不可忽略，而现有的并网逆变器参数辨识方法忽略了锁相环参数，仅针对除锁相环参数以外的参数进行辨识。例如：
文献2“孙黎霞，林雪，金宇清，潘磊，施浩波.基于粒子群优化算法的并网光伏发电单元建模[J].电网技术，2015，39(05)：1213-1218.”提出在忽略前馈和耦合的影响下，分别建立电流内环和电压外环的控制模型，然后采用粒子群优化算法分步辨识模型中的未知参数，但此种辨识方法是在dq坐标系下建立起电流内环和电压外环控制模型，因而忽略了锁相环参数，所以在辨识过程中也没有考虑锁相环参数的辨识。
文献3“孙丽玲.并网逆变器惯量和阻尼参数辨识研究[D].合肥工业大学，2020.”对虚拟同步机控制的并网逆变器建立了考虑锁相环的统一状态空间模型并对其进行参数辨识，但是在实际的辨识过程中锁相环带宽远远大于系统带宽，所以算法不能准确找到锁相环参数信息，也无法对锁相环参数进行辨识。
综上所述，逆变器参数辨识方法的现有技术中还存在以下不足：
1、阻抗测量环节存在误差，导致辨识误差存在不确定性，并且通过提高测量设备精度的方案所需成本高。
2、对逆变器参数进行辨识时是在dq坐标系下建模并施加扰动，或者在辨识过程中直接忽略了锁相环参数的影响，这些原因导致了锁相环参数无法被辨识出来。
3、锁相环与电流环等其他控制环耦合紧密，对锁相环施加的扰动也会对其他控制环产生影响，所以无法通过传统的施加扰动观察响应的方法来对锁相环参数进行辨识。
发明内容
本发明为了解决传统并网逆变器参数辨识过程中由于建模和锁相环与其他控制环耦合而导致锁相环参数无法辨识以及阻抗测量误差的存在而导致锁相环参数辨识不准确的问题，提供了一种基于并网逆变器序阻抗模型、参数敏感性因子、粒子群优化算法与平滑频段确定方法的并网逆变器锁相环参数辨识方法，所述锁相环参数包括锁相环比例系数和锁相环积分系数。
为解决本发明的技术问题，本发明提供了一种基于参数敏感性因子的并网逆变器锁相环参数辨识方法，所述锁相环参数包括锁相环比例系数和锁相环积分系数；所述并网逆变器锁相环参数辨识方法所涉及的拓扑结构包括并网逆变器、电网和锁相环参数辨识装置；所述锁相环参数辨识装置包括阻抗测量及计算单元和参数估计及数据处理单元；所述阻抗测量及计算单元包括采样单元和计算单元，所述采样单元接入电网与并网逆变器相连接的公共耦合点PCC处，实现对公共耦合点PCC处电压和电流的采样，所述计算单元与采样单元的输出端相连接，计算单元与参数估计及数据处理单元在算法内部实现数据交互；
所述参数辨识方法建立含有未知并网逆变器锁相环参数的待估计阻抗，提出用参数敏感性因子得到锁相环带宽段下参数敏感性因子满足要求的频率段，然后使用粒子群优化算法结合在锁相环带宽段下参数敏感性因子满足要求的频率段内测量得到的阻抗进行参数估计，获得多组并网逆变器锁相环估计参数，对多组并网逆变器锁相环估计参数进行处理，确定锁相环参数估计值的平滑频段，最终获得参数辨识结果，具体的，包括以下步骤：
步骤1，参数设定
包括以下参数：锁相环比例系数辨识下限KpPLLmin，锁相环比例系数辨识上限KpPLLmax，锁相环积分系数辨识下限KiPLLmin，锁相环积分系数辨识上限KiPLLmax；
步骤2，记并网逆变器的正序阻抗为Zip(KpPLL，KiPLL，f)，其中，KpPLL为待辨识锁相环比例系数，KiPLL为待辨识锁相环积分系数，f为频率，即正序阻抗Zip(KpPLL，KiPLL，f)仅包含待辨识锁相环比例系数KpPLL和待辨识锁相环积分系数KiPLL这2个未知参数；
步骤3，求解参数敏感性因子满足要求的频率段FSFS
步骤3.1，记锁相环比例系数参数敏感性因子为Fp(f)，锁相环积分系数参数敏感性因子为Fi(f)，其表达式如下：
式中的Re()表示对()中的参数取实部，Im()表示对()中的参数取虚部，KpPLLt为实际锁相环比例系数，KiPLLt为实际锁相环积分系数，其表达式分别如下：
其中，V1是基波相电压的峰值，fBW是锁相环带宽；
设锁相环带宽fBW的取值范围为[fBWmin，fBWmax]，并记为锁相环带宽段fBWd，其中，fBWmin是锁相环带宽段fBWd的取值下限，fBWmax是锁相环带宽段fBWd的取值上限；
步骤3.2，将锁相环带宽fBW下锁相环比例系数参数敏感性因子Fp(f)的最大值记为Fpmax，将锁相环带宽fBW下，由Fp(f)＞μ％·Fpmax确定的频率范围认定为锁相环带宽fBW下锁相环比例系数参数敏感性因子满足要求的频率段，并记为第一频率段FSFSps(fBW)，其中的μ是常数；
令锁相环带宽段fBWd下锁相环比例系数参数敏感性因子满足要求的频率段FSFSp是锁相环带宽段fBWd内每个第一频率段FSFSps(fBW)的并集频率段，其表达式如下：
FSFSp＝∪FSFSps(fBW)
步骤3.3，将锁相环带宽fBW下锁相环积分系数参数敏感性因子Fi(f)的最大值记为Fimax，将锁相环带宽fBW下，由Fi(f)＞μ％·Fimax确定的频率范围记为锁相环带宽fBW下锁相环积分系数参数敏感性因子满足要求的频率段，并记为第二频率段FSFSis(fBW)；
令锁相环带宽段fBWd下锁相环积分系数参数敏感性因子满足要求的频率段FSFSi是锁相环带宽段fBWd内每个第二频率段FSFSis(fBW)的并集频率段，其表达式如下：
FSFSi＝∪FSFSis(fBW)
步骤3.4，定义参数敏感性因子满足要求的频率段FSFS为锁相环带宽段fBWd下锁相环比例系数参数敏感性因子满足要求的频率段FSFSp和锁相环带宽段fBWd下锁相环积分系数参数敏感性因子满足要求的频率段FSFSi两者的并集频率段，其表达式如下：
FSFS＝FSFSpUFSFSi
步骤4，将锁相环参数辨识装置接入公共耦合点PCC，并通过采样单元对公共耦合点PCC处的三相线电压和三相电流进行N次采样，采样的频率范围是步骤3得到的参数敏感性因子满足要求的频率段FSFS，将N次采样中的任意一次采样记为第y次采样，与第y次采样对应的频率记为阻抗测量频率fy，与阻抗测量频率fy对应的三相线电压采样值记为三相线电压Uy，与阻抗测量频率fy对应的三相电流采样值记为三相电流Iy，y＝1，2，...，N，N为正整数；
步骤5，将步骤4采样得到的N个三相线电压Uy和N个三相电流Iy送入计算单元，得到N个测量正序阻抗Zrpy，y＝1，2，...，N；
步骤6，将N个阻抗测量频率fy和N个测量正序阻抗Zrpy送入参数识别及数据处理单元，使用N个测量正序阻抗Zrpy进行参数识别；
步骤6.1，定义待估计正序阻抗Zipdy(KpPLLdy，KiPLLdy，fddy)，其中的KpPLLdy为待估计锁相环比例系数，KiPLLdy为待估计锁相环积分系数，fddy为待估计正序阻抗Zipdy中的频率，y＝1，2，...，N；
令fddy＝fy，则待估计正序阻抗Zipdy表达为：Zipdy(KpPLLdy，KiPLLdy，fy)；
步骤6.2，使用粒子群优化算法使待估计正序阻抗Zipdy的实部和虚部与测量正序阻抗Zrpy的实部和虚部一一逼近，直至满足估计条件，得到锁相环比例系数估计值KpPLLFy和锁相环积分系数估计值KiPLLFy，具体的，建立粒子群优化算法函数fo(KpPLLdy，KiPLLdy)，通过计算得到粒子群群体最佳位置gbest，gbest＝(gbest1 gbest2)，其中gbest1为粒子群群体最佳位置第一列，gbest2为粒子群群体最佳位置第二列；
则锁相环比例系数估计值KpPLLFy和锁相环积分系数估计值KiPLLFy分别为：
KpPLLFy＝gbest1
KiPLLFy＝gbest2
步骤6.3，按照步骤6.1-步骤6.2方法得到N个阻抗测量频率fy对应的N个锁相环比例系数估计值KpPLLFy和N个锁相环积分系数估计值KiPLLFy，y＝1，2，...，N；
步骤7，参数估计及数据处理单元根据预存的算法，对步骤6得到的N个锁相环比例系数估计值KpPLLFy和锁相环积分系数估计值KiPLLFy进行处理，得到锁相环比例系数辨识结果KpPLLF和锁相环积分系数辨识结果KiPLLF；
步骤7.1，参数估计及数据处理单元根据预存的算法，对步骤6得到的N个锁相环比例系数估计值KpPLLFy进行处理，得到锁相环比例系数辨识结果KpPLLF，具体如下：
将N个锁相环比例系数估计值KpPLLFy按阻抗测量频率fy的值从小到大依次进行重新排序，并将其中任意一个重新排序后的锁相环比例系数估计值记为重排锁相环比例系数估计值其中Ap为重排锁相环比例系数估计值/>的序号，并记为第一重排序号Ap，Ap＝1，2，...，N；
考察N个重排锁相环比例系数估计值的分布状态，设在N个重排锁相环比例系数估计值/>中有Γp个满足条件E的第一待选频段，所述条件E为：在该第一待选频段中，所有的重排锁相环比例系数估计值/>均围绕其中的一个重排锁相环比例系数估计值/>上下波动，且上下波动的值小于该重排锁相环比例系数估计值/>的σ％，σ为设定的常数；Γp为正整数且Γp≥1；
若Γp＝1，则选择第一待选频段为第一平滑频段SFSp；若Γp＞1，则在Γp个满足条件E的第一待选频段中选择重排锁相环比例系数估计值个数最多的第一待选频段记为第一平滑频段SFSp；
记录第一平滑频段SFSp的以下参数：第一平滑频段SFSp中第一重排序号Ap的最小值Apmin和第一重排序号Ap的最大值Apmax，第一平滑频段SFSp中阻抗测量频率fy的最小值fpmin和阻抗测量频率fy的最大值fpmax；
锁相环比例系数辨识结果KpPLLF的表达式如下：
步骤7.2，参数估计及数据处理单元根据预存的算法，对步骤6得到的N个锁相环积分系数估计值KiPLLFy进行处理，得到锁相环积分系数辨识结果KiPLLF，具体如下：
将N个锁相环积分系数估计值KiPLLFy按阻抗测量频率fy的值从小到大依次进行重新排序，并将其中任意一个重新排序后的锁相环积分系数估计值记为重排锁相环积分系数估计值其中Ai为重排锁相环积分系数估计值/>的序号，并记为第二重排序号Ai，Ai＝1，2，...，N；
考察N个重排锁相环积分系数估计值的分布状态，设在N个重排锁相环积分系数估计值/>中有Γi个满足条件F的第二待选频段，所述条件F为：在该第二待选频段中，所有的重排锁相环积分系数估计值/>均围绕其中的一个重排锁相环积分系数估计值/>上下波动，且上下波动的值小于该重排锁相环积分系数估计值/>的σ％，σ为设定的常数；Γi为正整数且Γi≥1；
若Γi＝1，则选择第二待选频段为第二平滑频段SFSi；若Γi＞1，则在Γi个满足条件F的第二待选频段中选择重排锁相环积分系数估计值个数最多的第二待选频段记为第二平滑频段SFSi；
记录第二平滑频段SFSi的以下参数：第二平滑频段SFSi中第二重排序号Ai的最小值Aimin和第二重排序号Ai的最大值Aimax，第二平滑频段SFSi中阻抗测量频率fy的最小值fimin和阻抗测量频率fy的最大值fimax；
锁相环积分系数辨识结果KiPLLF的表达式如下：
本发明基于参数敏感性因子的并网逆变器锁相环参数辨识方法，综合了并网逆变器序阻抗建模、参数敏感性因子、阻抗测量、粒子群优化算法和平滑频段确定方法的优点。其有益效果具体体现如下：
1)建立起的并网逆变器序阻抗模型能够克服传统并网逆变器建模的缺点，解决电流环与锁相环两者耦合紧密的问题，可以用于辨识并网逆变器锁相环参数。
2)提出的参数敏感性因子可以对阻抗测量频率段的选取进行预处理，缩短了阻抗测量的频率范围，降低了阻抗测量误差导致的参数辨识不准确，同时避免了为提高测量设备精度而带来的高成本。
3)阻抗测量等操作对并网逆变器系统影响很小，不影响并网系统的稳定运行，能够在线辨识出并网逆变器的锁相环参数。
4)使用粒子群优化算法对未知参数进行参数估计，克服了传统最小二乘法等算法不能辨识非线性模型的缺点。
5)提出锁相环参数估计值的平滑频段的选定方法，减小了辨识过程中偶然误差对辨识结果造成的影响。
附图说明
图1为本发明所涉及的拓扑图。
图2为本发明实施例中建立的仿真电路结构图。
图3为本发明方法的流程图。
图4为锁相环比例系数辨识结果的辨识示意图。
图5为锁相环积分系数辨识结果的辨识示意图。
具体实施方式
以下结合附图，对本发明进行进一步详细说明。
图1所示为本发明基于参数敏感性因子的并网逆变器锁相环参数辨识方法所涉及的拓扑图。由图1可见，所述并网逆变器锁相环参数辨识方法所涉及的拓扑结构包括并网逆变器10、电网40和锁相环参数辨识装置；所述锁相环参数辨识装置包括阻抗测量及计算单元20和参数估计及数据处理单元30；所述阻抗测量及计算单元20包括采样单元201和计算单元202，所述采样单元201接入电网40与并网逆变器10相连接的公共耦合点PCC处，实现对公共耦合点PCC处电压和电流的采样，所述计算单元202与采样单元201的输出端相连接，计算单元202与参数估计及数据处理单元30在算法内部实现数据交互。
图2为本实施例中建立的基于参数敏感性因子的并网逆变器锁相环参数辨识方法仿真电路结构图，并网逆变器10直流电压Vdc＝1500V，额定功率PN＝300kW，滤波器电感L＝1.4mH，电网40频率f1＝60Hz，电网线电压有效值Vg＝690V，电流控制器输入为PCC点处三相电流Iy，锁相环输入为PCC点处三相线电压Uy，输出为相角θPLL。
图3为本发明方法的流程图。由该图可见，本发明是一种基于参数敏感性因子的并网逆变器锁相环参数辨识方法，所述锁相环参数包括锁相环比例系数和锁相环积分系数。所述参数辨识方法建立含有未知并网逆变器锁相环参数的待估计阻抗，提出用参数敏感性因子得到锁相环带宽段下参数敏感性因子满足要求的频率段，然后使用粒子群优化算法结合在锁相环带宽段下参数敏感性因子满足要求的频率段内测量得到的阻抗进行参数估计，获得多组并网逆变器锁相环估计参数，对多组并网逆变器锁相环估计参数进行处理，确定锁相环参数估计值的平滑频段，最终获得参数辨识结果，具体的，包括以下步骤：
步骤1，参数设定
包括以下参数：锁相环比例系数辨识下限KpPLLmin，锁相环比例系数辨识上限KpPLLmax，锁相环积分系数辨识下限KiPLLmin，锁相环积分系数辨识上限KiPLLmax。
在本实施例中，KpPLLmin＝0，KpPLLmax＝10，KiPLLmin＝0，KiPLLmax＝1000。
步骤2，记并网逆变器10的正序阻抗为Zip(KpPLL，KiPLL，f)，其中，KpPLL为待辨识锁相环比例系数，KiPLL为待辨识锁相环积分系数，f为频率，即正序阻抗Zip(KpPLL，KiPLL，f)仅包含待辨识锁相环比例系数KpPLL和待辨识锁相环积分系数KiPLL这2个未知参数。
步骤3，求解参数敏感性因子满足要求的频率段FSFS
步骤3.1，记锁相环比例系数参数敏感性因子为Fp(f)，锁相环积分系数参数敏感性因子为Fi(f)，其表达式如下：
式中的Re()表示对()中的参数取实部，Im()表示对()中的参数取虚部，KpPLLt为实际锁相环比例系数，KiPLLt为实际锁相环积分系数，其表达式分别如下：
其中，V1是基波相电压的峰值，fBW是锁相环带宽；
设锁相环带宽fBW的取值范围为[fBWmin，fBWmax]，并记为锁相环带宽段fBWd，其中，fBWmin是锁相环带宽段fBWd的取值下限，fBWmax是锁相环带宽段fBWd的取值上限；
步骤3.2，将锁相环带宽fBW下锁相环比例系数参数敏感性因子Fp(f)的最大值记为Fpmax，将锁相环带宽fBW下，由Fp(f)＞μ％·Fpmax确定的频率范围认定为锁相环带宽fBW下锁相环比例系数参数敏感性因子满足要求的频率段，并记为第一频率段FSFSps(fBW)，其中的μ是常数；
令锁相环带宽段fBWd下锁相环比例系数参数敏感性因子满足要求的频率段FSFSp是锁相环带宽段fBWd内每个第一频率段FSFSps(fBW)的并集频率段，其表达式如下：
FSFSp＝∪FSFSps(fBW)
步骤3.3，将锁相环带宽fBW下锁相环积分系数参数敏感性因子Fi(f)的最大值记为Fimax，将锁相环带宽fBW下，由Fi(f)＞μ％·Fimax确定的频率范围记为锁相环带宽fBW下锁相环积分系数参数敏感性因子满足要求的频率段，并记为第二频率段FSFSis(fBW)；
令锁相环带宽段fBWd下锁相环积分系数参数敏感性因子满足要求的频率段FSFSi是锁相环带宽段fBWd内每个第二频率段FSFSis(fBW)的并集频率段，其表达式如下：
FSFSi＝∪FSFSis(fBW)
步骤3.4，定义参数敏感性因子满足要求的频率段FSFS为锁相环带宽段fBWd下锁相环比例系数参数敏感性因子满足要求的频率段FSFSp和锁相环带宽段fBWd下锁相环积分系数参数敏感性因子满足要求的频率段FSFSi两者的并集频率段，其表达式如下：
FSFS＝FSFSp∪FSFSi。
在本实施例中，锁相环带宽段fBWd为[10，30]，由此得到的参数敏感性因子满足要求的频率段FSFS为[55，99]，μ＝80。
步骤4，将锁相环参数辨识装置接入公共耦合点PCC，并通过采样单元201对公共耦合点PCC处的三相线电压和三相电流进行N次采样，采样的频率范围是步骤3得到的参数敏感性因子满足要求的频率段FSFS，将N次采样中的任意一次采样记为第y次采样，与第y次采样对应的频率记为阻抗测量频率fy，与阻抗测量频率fy对应的三相线电压采样值记为三相线电压Uy，与阻抗测量频率fy对应的三相电流采样值记为三相电流Iy，y＝1，2，...，N，N为正整数。
在本实施例中，N＝45。
步骤5，将步骤4采样得到的N个三相线电压Uy和N个三相电流Iy送入计算单元202，得到N个测量正序阻抗Zrpy，y＝1，2，...，N。
在本实施例中，共得到45个测量正序阻抗Zrpy
步骤6，将N个阻抗测量频率fy和N个测量正序阻抗Zrpy送入参数识别及数据处理单元30，使用N个测量正序阻抗Zrpy进行参数识别；
步骤6.1，定义待估计正序阻抗Zipdy(KpPLLdy，KiPLLdy，fddy)，其中的KpPLLdy为待估计锁相环比例系数，KiPLLdy为待估计锁相环积分系数，fddy为待估计正序阻抗Zipdy中的频率，y＝1，2，...，N；
令fddy＝fy，则待估计正序阻抗Zipdy表达为：Zipdy(KpPLLdy，KiPLLdy，fy)；
步骤6.2，使用粒子群优化算法使待估计正序阻抗Zipdy的实部和虚部与测量正序阻抗Zrpy的实部和虚部一一逼近，直至满足估计条件，得到锁相环比例系数估计值KpPLLFy和锁相环积分系数估计值KiPLLFy，具体的，建立粒子群优化算法函数fo(KpPLLdy，KiPLLdy)，通过计算得到粒子群群体最佳位置gbest，gbest＝(gbest1 gbest2)，其中gbest1为粒子群群体最佳位置第一列，gbest2为粒子群群体最佳位置第二列；
则锁相环比例系数估计值KpPLLFy和锁相环积分系数估计值KiPLLFy分别为：
KpPLLFy＝gbest1
KiPLLFy＝gbest2
步骤6.3，按照步骤6.1-步骤6.2方法得到N个阻抗测量频率fy对应的N个锁相环比例系数估计值KpPLLFy和N个锁相环积分系数估计值KiPLLFy，y＝1，2，...，N。
步骤7，参数估计及数据处理单元30根据预存的算法，对步骤6得到的N个锁相环比例系数估计值KpPLLFy和锁相环积分系数估计值KiPLLFy进行处理，得到锁相环比例系数辨识结果KpPLLF和锁相环积分系数辨识结果KiPLLF。
步骤7.1，参数估计及数据处理单元30根据预存的算法，对步骤6得到的N个锁相环比例系数估计值KpPLLFy进行处理，得到锁相环比例系数辨识结果KpPLLF，具体如下：
将N个锁相环比例系数估计值KpPLLFy按阻抗测量频率fy的值从小到大依次进行重新排序，并将其中任意一个重新排序后的锁相环比例系数估计值记为重排锁相环比例系数估计值其中Ap为重排锁相环比例系数估计值/>的序号，并记为第一重排序号Ap，Ap＝1，2，...，N；
考察N个重排锁相环比例系数估计值的分布状态，设在N个重排锁相环比例系数估计值/>中有Γp个满足条件E的第一待选频段，所述条件E为：在该第一待选频段中，所有的重排锁相环比例系数估计值/>均围绕其中的一个重排锁相环比例系数估计值/>上下波动，且上下波动的值小于该重排锁相环比例系数估计值/>的σ％，σ为设定的常数；Γp为正整数且Γp≥1；
若Γp＝1，则选择第一待选频段为第一平滑频段SFSp；若Γp＞1，则在Γp个满足条件E的第一待选频段中选择重排锁相环比例系数估计值个数最多的第一待选频段记为第一平滑频段SFSp；
记录第一平滑频段SFSp的以下参数：第一平滑频段SFSp中第一重排序号Ap的最小值Apmin和第一重排序号Ap的最大值Apmax，第一平滑频段SFSp中阻抗测量频率fy的最小值fpmin和阻抗测量频率fy的最大值fpmax；
锁相环比例系数辨识结果KpPLLF的表达式如下：
在本实施例中，σ＝0.5。
图4为锁相环比例系数辨识结果的辨识示意图，横坐标是阻抗测量频率fy，纵坐标是重排锁相环比例系数估计值由图4可见，第一平滑频段SFSp为[76，89]，即fpmin＝76，fpmax＝89，第一重排序号Ap的最小值Apmin＝22，第一重排序号Ap的最大值Apmax＝35，锁相环比例系数辨识结果KpPLLF为0.192。为了佐证，对结果进行了仿真，仿真中待辨识锁相环比例系数KpPLL设置值为0.193。
步骤7.2，参数估计及数据处理单元30根据预存的算法，对步骤6得到的N个锁相环积分系数估计值KiPLLFy进行处理，得到锁相环积分系数辨识结果KiPLLF，具体如下：
将N个锁相环积分系数估计值KiPLLFy按阻抗测量频率fy的值从小到大依次进行重新排序，并将其中任意一个重新排序后的锁相环积分系数估计值记为重排锁相环积分系数估计值其中Ai为重排锁相环积分系数估计值/>的序号，并记为第二重排序号Ai，Ai＝1，2，...，N；
考察N个重排锁相环积分系数估计值的分布状态，设在N个重排锁相环积分系数估计值/>中有Γi个满足条件F的第二待选频段，所述条件F为：在该第二待选频段中，所有的重排锁相环积分系数估计值/>均围绕其中的一个重排锁相环积分系数估计值/>上下波动，且上下波动的值小于该重排锁相环积分系数估计值/>的σ％，σ为设定的常数；Γi为正整数且Γi≥1；
若Γi＝1，则选择第二待选频段为第二平滑频段SFSi；若Γi＞1，则在Γi个满足条件F的第二待选频段中选择重排锁相环积分系数估计值个数最多的第二待选频段记为第二平滑频段SFSi；
记录第二平滑频段SFSi的以下参数：第二平滑频段SFSi中第二重排序号Ai的最小值Aimin和第二重排序号Ai的最大值Aimax，第二平滑频段SFSi中阻抗测量频率fy的最小值fimin和阻抗测量频率fy的最大值fimax；
锁相环积分系数辨识结果KiPLLF的表达式如下：
图5为锁相环积分系数辨识结果的辨识示意图，横坐标是阻抗测量频率fy，纵坐标是重排锁相环积分系数估计值由图5可见，第二平滑频段SFSi为[63，77]，即fimin＝63，fimax＝77，第二重排序号Ai的最小值Aimin＝9，第二重排序号Ai的最大值Aimax＝23，锁相环积分系数辨识结果KiPLLF为36.361。为了佐证，对结果进行了仿真，仿真中待辨识锁相环积分系数KiPLL设置值为36.400。
至此，并网逆变器锁相环参数辨识结束。
在以上步骤中，步骤6.2所述的粒子群群体最佳位置gbest获取方式如下：
建立粒子群优化算法函数fo(KpPLLdy，KiPLLdy)，其表达式如下：
fo(KpPLLdy，KiPLLdy)＝||Re(zrpy)-Re[Zipdy(KpPLLdy，KiPLLdy，fy)]|+|Im(Zrpy)-Im[Zipdy(KpPLLdy，KiPLLdy，fy)]||
使用粒子群优化算法寻找满足要求的待估计锁相环比例系数KpPLLdy和待估计锁相环积分系数KiPLLdy，使粒子群优化算法函数fo(KpPLLdy，KiPLLdy)的值最小，具体步骤如下：
步骤6.2.1，设置粒子群规模为n，粒子群优化函数变量个数d为2，粒子群优化算法学习因子1为c1，粒子群优化算法学习因子2为c2，粒子群优化算法速度权重W1，粒子群优化算法位置权重W2，粒子群优化迭代次数K，粒子速度最小值vmin，粒子速度最大值vmax；
步骤6.2.2，随机生成n行2列的二维矩阵记为粒子群优化算法粒子种群xn×2，且第一列元素的随机数范围在锁相环比例系数辨识下限KpPLLmin到锁相环比例系数辨识上限KpPLLmax之间，粒子群优化算法粒子种群xn×2第二列元素的随机数范围在锁相环积分系数辨识下限KiPLLmin到锁相环积分系数辨识上限KiPLLmax之间，其中，记n行中任意一行为第a行，a＝1，2，...，n，a为正整数，2列中任意一列为第b列，b＝1，2，b为正整数，粒子群优化算法粒子种群xn×2中的任意一个粒子记为Xab，a＝1，2，...，n，b＝1，2，粒子群优化算法粒子种群xn×2中的任意一行粒子记为Xa，a＝1，2，...，n，任意一行粒子Xa包含任意一行第一列的粒子Xa1，a＝1，2，...，n和任意一行第二列的粒子Xa2，a＝1，2，...，n；
粒子群优化算法粒子种群xn×2和粒子群优化算法粒子种群xn×2中的任意一行粒子Xa如下：
Xa＝(Xa1 Xa2)；
随机生成n行2列的二维矩阵记为粒子群优化算法粒子速度vn×2，且每个元素的随机数范围在粒子速度最小值vmin和粒子速度最大值vmax之间，其中，粒子群优化算法粒子速度vn×2中的任意一个粒子速度记为Vab，a＝1，2，...，n，b＝1，2，粒子群优化算法粒子速度vn×2中的任意一行粒子速度记为Va，a＝1，2，...，n，任意一行粒子速度Va包含任意一行第一列的粒子速度Va1，a＝1，2，...，n和任意一行第二列的粒子Va2，a＝1，2，...，n；
粒子群优化算法粒子速度vn×2和粒子群优化算法粒子速度vn×2中的任意一行粒子速度Va如下：
Va＝(Va1 Va2)；
步骤6.2.3，计算粒子群优化算法粒子种群xn×2中的任意一行粒子Xa的种群适应度，记为粒子群适应度Fita，a＝1，2，...，n，粒子群适应度Fita的值如下：
Fita＝fo(Xa1，Xa2)；
记粒子群任意一行最佳位置为pbesta，a＝1，2，...，n，pbesta＝(pbesta1pbesta2)，其中，粒子群任意一行最佳位置第一列为pbesta1，a＝1，2，...，n，粒子群任意一行最佳位置第二列为pbesta2，a＝1，2，...，n，粒子群任意一行最佳位置为pbesta如下：
pbesta1＝Xa1
pbesta2＝Xa2
pbesta＝Xa；
寻找Fita中的最小值记为粒子群最小适应度Fitamin，粒子群最小适应度Fitamin所对应的那一行粒子，记为粒子群最小适应度行粒子Xamin，Xamin＝(Xa1minXa2min)，其中，Xa1min为粒子群最小适应度行第一列粒子，Xa2min为粒子群最小适应度行第二列粒子；
记粒子群群体最佳位置为gbest，gbest＝(gbest1 gbest2)，其中gbest1为粒子群群体最佳位置第一列，gbest2为粒子群群体最佳位置第二列，粒子群群体最佳位置为gbest如下：
gbest1＝Xa1min
gbest2＝Xa2min
gbest＝Xamin；
步骤6.2.4，更新粒子群优化算法粒子种群xn×2和粒子群优化算法粒子速度vn×2，并进行K次更新，记K次更新中任意一次更新为第k次更新，k＝1，2，...，K，k为正整数；
记第k次更新时的粒子群优化算法粒子种群为xn×2 k，第k次更新时的粒子群优化算法粒子种群xn×2 k中的任意一个粒子记为Xab k，a＝1，2，...，n，b＝1，2，第k次更新时的粒子群优化算法粒子种群xn×2 k中的任意一行粒子记为Xa k，a＝1，2，...，n，任意一行粒子Xa k包含任意一行第一列的粒子Xa1 k，a＝1，2，...，n和任意一行第二列的粒子Xa2 k，a＝1，2，...，n，第k次更新时的粒子群优化算法粒子速度为vn×2 k，第k次更新时的粒子群优化算法粒子速度vn×2 k中的任意一个粒子速度记为Vab k，a＝1，2，...，n，b＝1，2，第k次更新时的粒子群优化算法粒子速度vn×2 k中的任意一行粒子速度记为Va k，a＝1，2，...，n，任意一行粒子速度Va k包含任意一行第一列的粒子速度Va1 k，a＝1，2，...，n和任意一行第二列的粒子Va2 k，a＝1，2，...，n，k＝1，2，...，K；
第k次更新时的粒子群优化算法粒子种群xn×2 k和第k次更新时的粒子群优化算法粒子种群xn×2 k中的任意一行粒子Xa k如下：
Xa k＝(Xa1 k Xa2 k)；
第k次更新时的粒子群优化算法粒子速度vn×2 k和第k次更新时的粒子群优化算法粒子速度vn×2 k中的任意一行粒子Va k如下：
Va k＝(Va1 k Va2 k)；
第k次更新的步骤如下：
(1)记第k次更新时的粒子群优化算法粒子速度为vn×2 k，更新公式如下：
Va k＝W1×Va (k-1)+c1×rand×[pbest-Xa (k-1)]+c2×rand×[gbest-Xa (k-1)]
其中，rand为0到1之间的随机数，k＝1，2，...，K；
当第一次更新即k＝1时，Va (k-1)＝Va 0＝Va，Xa (k-1)＝Xa 0＝Xa；
若Va1 k和Va2 k中有一个大于vmax则其值替换成vmax，若Va1 k和Va2 k中有一个小于vmin则其值替换成vmin；
(2)第k次更新时的粒子群优化算法粒子种群xn×2 k，更新公式如下：
Xa k＝Xa (k-1)+W2×Va (k-1)
其中，k＝1，2，...，K；
当第一次更新即k＝1时，Xa (k-1)＝Xa 0＝Xa；
若Xa1 k和Xa2 k中有一个大于KpPLLmax则其值替换成KpPLLmax，Xa1 k和Xa2 k中有一个小于KpPLLmin则其值替换成KpPLLmin；
(3)重新计算适应度，第k次更新时的粒子群优化算法粒子种群xn×2 k任意一行粒子Xa k的种群适应度，记为粒子群适应度Fita k，a＝1，2，...，n，粒子群适应度Fita k的值如下：
Fita k＝fo(Xa1 k，Xa2 k)
其中，k＝1，2，...，K；
若Fita k＜fo(pbesta1，pbesta2)，则pbesta1＝Xa1 k，pbesta2＝Xa2 k，pbesta＝Xa k；反之，粒子群任意一行最佳位置pbesta保持不变；
寻找Fita k中的最小值记为第k次更新时粒子群最小适应度Fitamin k，粒子群最小适应度Fitamin k所对应的那一行粒子，记为第k次更新时粒子群最小适应度行粒子Xamin k，Xamin k＝(Xa1min kXa2min k)，其中，Xa1min k为第k次更新时粒子群最小适应度行第一列粒子，Xa2min k为第k次更新时粒子群最小适应度行第二列粒子；
若Fitamin k＜fo(gbest1，gbest2)，则gbest1＝Xa1min k，gbest2＝Xa2min k，gbest＝Xamin k；反之，粒子群群体最佳位置gbest保持不变；
步骤6.2.5，输出粒子群群体最佳位置gbest。
相对于传统的参数辨识方法因建模和扰动选择等方面的问题，无法对锁相环参数进行辨识以及因阻抗测量存在误差的问题，锁相环参数辨识不准确。而本发明所提出的基于参数敏感性因子的并网逆变器锁相环参数辨识方法，首先，建立了并网逆变器序阻抗模型，在建模上考虑了锁相环参数的影响；其次，提出用参数敏感性因子得到锁相环带宽段下参数敏感性因子满足要求的频率段，缩短了阻抗测量的频率范围，同时减小了阻抗测量的误差对参数辨识结果的影响；然后，使用阻抗测量的方法测量锁相环带宽段下参数敏感性因子满足要求的频率段内的序阻抗；接着，选用粒子群优化算法进行参数估计；最后处理多组并网逆变器锁相环估计参数，确定锁相环参数估计值的平滑频段，得出锁相环参数的辨识结果。本方法能够在线辨识锁相环参数，不影响系统的稳定运行，辨识出的锁相环参数为系统稳定性分析奠定了基础。
综上所述，该方法实现简单，只需建立含有未知锁相环参数的正序阻抗模型，利用提出的参数敏感性因子预处理得到阻抗测量的频率范围，使用粒子群优化算法结合测量阻抗进行参数估计，获得多组并网逆变器锁相环估计参数，处理多组并网逆变器锁相环估计参数，确定锁相环参数估计值的平滑频段，最终获得参数辨识结果，具有一定的可行性。
Claims (1)
1.一种基于参数敏感性因子的并网逆变器锁相环参数辨识方法，所述锁相环参数包括锁相环比例系数和锁相环积分系数；所述并网逆变器锁相环参数辨识方法所涉及的拓扑结构包括并网逆变器(10)、电网(40)和锁相环参数辨识装置；所述锁相环参数辨识装置包括阻抗测量及计算单元(20)和参数估计及数据处理单元(30)；所述阻抗测量及计算单元(20)包括采样单元(201)和计算单元(202)，所述采样单元(201)接入电网(40)与并网逆变器(10)相连接的公共耦合点PCC处，实现对公共耦合点PCC处电压和电流的采样，所述计算单元(202)与采样单元(201)的输出端相连接，计算单元(202)与参数估计及数据处理单元(30)在算法内部实现数据交互；
其特征在于，所述参数辨识方法建立含有未知并网逆变器锁相环参数的待估计阻抗，提出用参数敏感性因子得到锁相环带宽段下参数敏感性因子满足要求的频率段，然后使用粒子群优化算法结合在锁相环带宽段下参数敏感性因子满足要求的频率段内测量得到的阻抗进行参数估计，获得多组并网逆变器锁相环估计参数，对多组并网逆变器锁相环估计参数进行处理，确定锁相环参数估计值的平滑频段，最终获得参数辨识结果，具体的，包括以下步骤：
步骤1，参数设定
包括以下参数：锁相环比例系数辨识下限KpPLLmin，锁相环比例系数辨识上限KpPLLmax，锁相环积分系数辨识下限KiPLLmin，锁相环积分系数辨识上限KiPLLmax；
步骤2，记并网逆变器(10)的正序阻抗为Zip(KpPLL，KiPLL，f)，其中，KpPLL为待辨识锁相环比例系数，KiPLL为待辨识锁相环积分系数，f为频率，即正序阻抗Zip(KpPLL，KiPLL，f)仅包含待辨识锁相环比例系数KpPLL和待辨识锁相环积分系数KiPLL这2个未知参数；
步骤3，求解参数敏感性因子满足要求的频率段FSFS
步骤3.1，记锁相环比例系数参数敏感性因子为Fp(f)，锁相环积分系数参数敏感性因子为Fi(f)，其表达式如下：
式中的Re()表示对()中的参数取实部，Im()表示对()中的参数取虚部，KpPLLt为实际锁相环比例系数，KiPLLt为实际锁相环积分系数，其表达式分别如下：
其中，V1是基波相电压的峰值，fBW是锁相环带宽；
设锁相环带宽fBW的取值范围为[fBWmin，fBWmax]，并记为锁相环带宽段fBWd，其中，fBWmin是锁相环带宽段fBWd的取值下限，fBWmax是锁相环带宽段fBWd的取值上限；
步骤3.2，将锁相环带宽fBW下锁相环比例系数参数敏感性因子Fp(f)的最大值记为Fpmax，将锁相环带宽fBW下，由Fp(f)＞μ％·Fpmax确定的频率范围认定为锁相环带宽fBW下锁相环比例系数参数敏感性因子满足要求的频率段，并记为第一频率段FSFSps(fBW)，其中的μ是常数；
令锁相环带宽段fBWd下锁相环比例系数参数敏感性因子满足要求的频率段FSFSp是锁相环带宽段fBWd内每个第一频率段FSFSps(fBW)的并集频率段，其表达式如下：
FSFSp＝∪FSFSps(fBW)
步骤3.3，将锁相环带宽fBW下锁相环积分系数参数敏感性因子Fi(f)的最大值记为Fimax，将锁相环带宽fBW下，由Fi(f)＞μ％·Fimax确定的频率范围记为锁相环带宽fBW下锁相环积分系数参数敏感性因子满足要求的频率段，并记为第二频率段FSFSis(fBW)；
令锁相环带宽段fBWd下锁相环积分系数参数敏感性因子满足要求的频率段FSFSi是锁相环带宽段fBWd内每个第二频率段FSFSis(fBW)的并集频率段，其表达式如下：
FSFSi＝∪FSFSis(fBW)
步骤3.4，定义参数敏感性因子满足要求的频率段FSFS为锁相环带宽段fBWd下锁相环比例系数参数敏感性因子满足要求的频率段FSFSp和锁相环带宽段fBWd下锁相环积分系数参数敏感性因子满足要求的频率段FSFSi两者的并集频率段，其表达式如下：
FSFS＝FSFSp∪FSFSi
步骤4，将锁相环参数辨识装置接入公共耦合点PCC，并通过采样单元(201)对公共耦合点PCC处的三相线电压和三相电流进行N次采样，采样的频率范围是步骤3得到的参数敏感性因子满足要求的频率段FSFS，将N次采样中的任意一次采样记为第y次采样，与第y次采样对应的频率记为阻抗测量频率fy，与阻抗测量频率fy对应的三相线电压采样值记为三相线电压Uy，与阻抗测量频率fy对应的三相电流采样值记为三相电流Iy，y＝1，2，...，N，N为正整数；
步骤5，将步骤4采样得到的N个三相线电压Uy和N个三相电流Iy送入计算单元(202)，得到N个测量正序阻抗Zrpy，y＝1，2，...，N；
步骤6，将N个阻抗测量频率fy和N个测量正序阻抗Zrpy送入参数识别及数据处理单元(30)，使用N个测量正序阻抗Zrpy进行参数识别；
步骤6.1，定义待估计正序阻抗Zipdy(KpPLLdy，KiPLLdy，fddy)，其中的KpPLLdy为待估计锁相环比例系数，KiPLLdy为待估计锁相环积分系数，fddy为待估计正序阻抗Zipdy中的频率，y＝1，2，...，N；
令fddy＝fy，则待估计正序阻抗Zipdy表达为：Zipdy(KpPLLdy，KiPLLdy，fy)；
步骤6.2，使用粒子群优化算法使待估计正序阻抗Zipdy的实部和虚部与测量正序阻抗Zrpy的实部和虚部一一逼近，直至满足估计条件，得到锁相环比例系数估计值KpPLLFy和锁相环积分系数估计值KiPLLFy，具体的，建立粒子群优化算法函数fo(KpPLLdy，KiPLLdy)，通过计算得到粒子群群体最佳位置gbest，gbest＝(gbest1 gbest2)，其中gbest1为粒子群群体最佳位置第一列，gbest2为粒子群群体最佳位置第二列；
则锁相环比例系数估计值KpPLLFy和锁相环积分系数估计值KiPLLFy分别为：
KpPLLFy＝gbest1
KiPLLFy＝gbest2
步骤6.3，按照步骤6.1-步骤6.2方法得到N个阻抗测量频率fy对应的N个锁相环比例系数估计值KpPLLFy和N个锁相环积分系数估计值KiPLLFy，y＝1，2，...，N；
步骤7，参数估计及数据处理单元(30)根据预存的算法，对步骤6得到的N个锁相环比例系数估计值KpPLLFy和锁相环积分系数估计值KiPLLFy进行处理，得到锁相环比例系数辨识结果KpPLLF和锁相环积分系数辨识结果KiPLLF；
步骤7.1，参数估计及数据处理单元(30)根据预存的算法，对步骤6得到的N个锁相环比例系数估计值KpPLLFy进行处理，得到锁相环比例系数辨识结果KpPLLF，具体如下：
将N个锁相环比例系数估计值KpPLLFy按阻抗测量频率fy的值从小到大依次进行重新排序，并将其中任意一个重新排序后的锁相环比例系数估计值记为重排锁相环比例系数估计值其中Ap为重排锁相环比例系数估计值/>的序号，并记为第一重排序号Ap，Ap＝1，2，...，N；
考察N个重排锁相环比例系数估计值的分布状态，设在N个重排锁相环比例系数估计值/>中有Γp个满足条件E的第一待选频段，所述条件E为：在该第一待选频段中，所有的重排锁相环比例系数估计值/>均围绕其中的一个重排锁相环比例系数估计值上下波动，且上下波动的值小于该重排锁相环比例系数估计值/>的σ％，σ为设定的常数；Γp为正整数且Γp≥1；
若Γp＝1，则选择第一待选频段为第一平滑频段SFSp；若Γp＞1，则在Γp个满足条件E的第一待选频段中选择重排锁相环比例系数估计值个数最多的第一待选频段记为第一平滑频段SFSp；
记录第一平滑频段SFSp的以下参数：第一平滑频段SFSp中第一重排序号Ap的最小值Apmin和第一重排序号Ap的最大值Apmax，第一平滑频段SFSp中阻抗测量频率fy的最小值fpmin和阻抗测量频率fy的最大值fpmax；
锁相环比例系数辨识结果KpPLLF的表达式如下：
步骤7.2，参数估计及数据处理单元(30)根据预存的算法，对步骤6得到的N个锁相环积分系数估计值KiPLLFy进行处理，得到锁相环积分系数辨识结果KiPLLF，具体如下：
将N个锁相环积分系数估计值KiPLLFy按阻抗测量频率fy的值从小到大依次进行重新排序，并将其中任意一个重新排序后的锁相环积分系数估计值记为重排锁相环积分系数估计值其中Ai为重排锁相环积分系数估计值/>的序号，并记为第二重排序号Ai，Ai＝1，2，...，N；
考察N个重排锁相环积分系数估计值的分布状态，设在N个重排锁相环积分系数估计值/>中有Γi个满足条件F的第二待选频段，所述条件F为：在该第二待选频段中，所有的重排锁相环积分系数估计值/>均围绕其中的一个重排锁相环积分系数估计值上下波动，且上下波动的值小于该重排锁相环积分系数估计值/>的σ％，σ为设定的常数；Γi为正整数且Γi≥1；
若Γi＝1，则选择第二待选频段为第二平滑频段SFSi；若Γi＞1，则在Γi个满足条件F的第二待选频段中选择重排锁相环积分系数估计值个数最多的第二待选频段记为第二平滑频段SFSi；
记录第二平滑频段SFSi的以下参数：第二平滑频段SFSi中第二重排序号Ai的最小值Aimin和第二重排序号Ai的最大值Aimax，第二平滑频段SFSi中阻抗测量频率fy的最小值fimin和阻抗测量频率fy的最大值fimax；
锁相环积分系数辨识结果KiPLLF的表达式如下：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210250369.2A CN114552652B (zh) | 2022-03-10 | 2022-03-10 | 基于参数敏感性因子的并网逆变器锁相环参数辨识方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210250369.2A CN114552652B (zh) | 2022-03-10 | 2022-03-10 | 基于参数敏感性因子的并网逆变器锁相环参数辨识方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114552652A CN114552652A (zh) | 2022-05-27 |
| CN114552652B true CN114552652B (zh) | 2024-03-26 |
Family
ID=81663304
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210250369.2A Active CN114552652B (zh) | 2022-03-10 | 2022-03-10 | 基于参数敏感性因子的并网逆变器锁相环参数辨识方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114552652B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110518581A (zh) * | 2019-08-19 | 2019-11-29 | 上海交通大学 | 计及采样滤波锁相环的逆变器阻抗优化方法 |
| CN111506869A (zh) * | 2020-04-15 | 2020-08-07 | 浙江大学 | 基于多正弦波信号注入的并网逆变器频率耦合阻抗特性提取方法 |
| CN112104275A (zh) * | 2020-11-13 | 2020-12-18 | 国网浙江省电力有限公司电力科学研究院 | 一种无锁相环直接功率控制的dfig阻抗重塑控制方法 |
| CN112421683A (zh) * | 2021-01-22 | 2021-02-26 | 中国电力科学研究院有限公司 | 一种并网逆变器的多回路控制参数辨识方法及系统 |
| CN113690933A (zh) * | 2021-08-06 | 2021-11-23 | 合肥工业大学 | 并网逆变器锁相环参数辨识方法 |
| WO2022037581A1 (zh) * | 2020-08-17 | 2022-02-24 | 中国电力科学研究院有限公司 | 一种光伏逆变器控制参数的辨识方法和装置 |
-
2022
- 2022-03-10 CN CN202210250369.2A patent/CN114552652B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110518581A (zh) * | 2019-08-19 | 2019-11-29 | 上海交通大学 | 计及采样滤波锁相环的逆变器阻抗优化方法 |
| CN111506869A (zh) * | 2020-04-15 | 2020-08-07 | 浙江大学 | 基于多正弦波信号注入的并网逆变器频率耦合阻抗特性提取方法 |
| WO2022037581A1 (zh) * | 2020-08-17 | 2022-02-24 | 中国电力科学研究院有限公司 | 一种光伏逆变器控制参数的辨识方法和装置 |
| CN112104275A (zh) * | 2020-11-13 | 2020-12-18 | 国网浙江省电力有限公司电力科学研究院 | 一种无锁相环直接功率控制的dfig阻抗重塑控制方法 |
| CN112421683A (zh) * | 2021-01-22 | 2021-02-26 | 中国电力科学研究院有限公司 | 一种并网逆变器的多回路控制参数辨识方法及系统 |
| CN113690933A (zh) * | 2021-08-06 | 2021-11-23 | 合肥工业大学 | 并网逆变器锁相环参数辨识方法 |
Non-Patent Citations (2)
| Title |
|---|
| 三相并网逆变器锁相环频率特性分析及其稳定性研究;王程;陈新;张旸;陈杰;龚春英;;中国电机工程学报;20170705(13);183-193 * |
| 无锁相环直接功率控制双馈感应发电机系统阻抗建模和稳定性分析;年珩;胡彬;陈亮;徐韵扬;;中国电机工程学报;20201231(03);277-288 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114552652A (zh) | 2022-05-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108667048B (zh) | 新能源并网系统振荡稳定性的频域判稳方法及装置 | |
| CN109004652B (zh) | 一种含异步电机的直流多馈入系统受端电网强度获得方法 | |
| CN109713971B (zh) | 一种永磁同步电机的扰动抑制方法 | |
| CN108809273A (zh) | 基于lms自适应滤波的复数直接频率估计方法 | |
| CN110417050A (zh) | 一种并联结构大型风电场等值计算方法 | |
| CN112685969A (zh) | 一种大型火力发电厂冷端的优化方法 | |
| CN109194225A (zh) | 一种双馈电机参数在线辨识方法 | |
| WO2024114279A1 (zh) | 功率模块老化测试的漏电流校正方法、装置及测试系统 | |
| CN111308207B (zh) | 一种针对单相交流系统的dq阻抗测量方法 | |
| CN114552652B (zh) | 基于参数敏感性因子的并网逆变器锁相环参数辨识方法 | |
| CN107104430A (zh) | 一种配电网双解耦状态估计方法 | |
| CN113690933B (zh) | 并网逆变器锁相环参数辨识方法 | |
| CN111737919B (zh) | 一种适用于次同步振荡分析的直驱式风电场分群方法 | |
| CN116520095B (zh) | 故障测距方法、系统以及计算机可读存储介质 | |
| CN111553080B (zh) | 一种配电台区负荷动态等值非机理模型参数闭环辨识方法 | |
| CN116776525B (zh) | 一种数据驱动型igbt器件的驱动优化方法 | |
| CN115343570B (zh) | 一种电网阻抗在线辨识方法及装置 | |
| CN115117884A (zh) | 电力系统暂态稳定域边界的计算方法 | |
| CN112649746A (zh) | 一种结合电路等效和递推迭代的荷电状态估计方法 | |
| TWI803916B (zh) | 一種晶片測試機測量晶片內阻的方法 | |
| CN105092980B (zh) | 一种输入输出阻抗智能化测试方法 | |
| CN112332734B (zh) | 一种提高大范围调速响应能力的超高速电动空压机变电压扩稳控制系统及方法 | |
| CN115293090A (zh) | 基于重构数据处理的多谐波源责任量化方法 | |
| JP3620930B2 (ja) | 電力系統の特性推定装置および特性推定方法 | |
| CN116167232B (zh) | 一种dfig序阻抗模型辨识方法及系统 |
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 |