并联换流器锁相环参数计算方法与流程

1.本发明涉及电子换流器控制技术领域，尤其涉及一种基于相平面法的并联换流器锁相环参数计算方法。


背景技术：

2.换流器可以实现交直流的相互转供，提高电网的供配电效率及运行可靠性；但换流器与电网连接时需要解决换流器与电网电压或电流同步的问题；电网同步是一个自适应过程，一般通过控制算法生成一个内部参考信号并使该参考信号与特定的电网变量步调一致来实现。并联换流器通常采用电压作为同步输入信号，即基于电压的同步方式；其中，锁相环(phase-locked loop，pll)是现阶段最为流行的同步技术。在现有电力电子换流器控制技术中，锁相环参数在计算时大多基于锁相环本身的线性化模型，当电网较强时，可以认为锁相环的参数和电网特性是解耦的，所得到的整定参数是能让锁相环的输出快速准确地跟踪系统频率和相位角变化的，即实现电网同步。
3.随着新型电力系统建设的不断推进，风电、光伏发电的占比将快速提升，同步发电机将逐步退出，电力系统强度将变弱，从而导致风电、光伏发电等机组必须适应弱系统运行这一新场景。然而，电网较弱时并联换流器向电网注入的功率会对换流器电网接入点的电压产生比较大的影响，导致锁相环的参数和电网的特性耦合在一起，无法使锁相环的输出快速准确地跟踪系统频率和相位角变化，致使风电、光伏发电等机组的并联换流器容易与电网失去同步，从而不能正常发挥作用，甚至脱网，进而威胁到电网安全稳定运行。若考虑电网影响，现有锁相环参数的计算变得很复杂、繁重，且单纯的时域仿真分析，不仅计算量非常大、计算时间很长，还不能对换流器的暂态稳定性有一个较为全面的了解，也无法确定暂态稳定的临界点。因此，有必要采用新的方法来确定锁相环的参数，从而提升弱电网下换流器的同步稳定性。


技术实现要素：

4.本发明的目的在于提供一种并联换流器锁相环参数计算方法，基于相位偏差二阶系统模型，采用相平面法获取锁相环参数的终值及并联换流器的同步稳定临界点，能够对并联换流器的同步稳定性有一个较为全面地了解。
5.为了达到上述目的，本发明通过以下技术方案实现：
6.一种并联换流器锁相环参数计算方法，包括：
7.计算锁相环输出相位与电网电压相位的偏差；
8.根据所述锁相环输出相位与所述电网电压相位的偏差，构建相位偏差二阶系统模型；
9.根据阻尼参数和调整时间参数计算所述锁相环参数的初始值；以及
10.基于所述相位偏差二阶系统模型，采用相平面法和所述锁相环参数的初始值获取相位偏差的相轨迹，以得到所述锁相环参数的终值及并联换流器的同步稳定临界点。
11.优选地，所述锁相环参数包括比例参数和积分参数；且所述锁相环输出相位采用如下公式表示：
12.θ
pll
＝∫(k
ppllvpccq
+k
ipll
∫v
pccq
dt+ωg)dt
13.其中，θ
pll
为所述锁相环输出相位；k
ppll
为所述锁相环的比例参数；k
ipll
为所述锁相环的积分参数；v
pcc
为并网点电压；v
pccq
为并网点电压的q轴分量；ωg为所述电网角频率基准；t为参变量；
14.所述电网电压相位采用如下公式表示：
15.θg＝∫ωgdt
16.其中，θg为所述电网电压相对于参考坐标系的相位；
17.所述锁相环输出相位与所述电网电压相位的偏差采用如下公式进行计算：
18.δ＝θ
pll-θg19.＝∫(k
ppllvpccq
+k
ipll
∫v
pccq
dt)dt
20.其中，δ为所述锁相环输出相位与所述电网电压相位的偏差。
21.优选地，所述并网点电压v
pcc
的q轴分量采用如下公式进行计算：
[0022]vpccq
＝v
g sin(θ
g-θ
pll
)+lg·icm
·
ωg[0023]
其中，vg为所述电网电压幅值；lg为所述电网电感；i
cm
为并联换流器内部电流指令幅值。
[0024]
优选地，所述相位偏差二阶系统模型采用如下公式表示：
[0025][0026]
优选地，所述锁相环的比例参数和积分参数采用如下公式进行计算：
[0027][0028][0029]
其中，ts表示调整时间参数；ζ表示阻尼参数。
[0030]
优选地，根据阻尼参数和调整时间参数计算所述锁相环参数的初始值的步骤包括：
[0031]
保持所述阻尼参数和所述调整时间参数中一个值不变，另一个设置不同值，对应计算所述锁相环的比例参数和积分参数的第一初始值；以及
[0032]
将所述锁相环的比例参数和积分参数的第一初始值构成第一初始条件。
[0033]
优选地，基于所述相位偏差二阶系统模型，采用相平面法和所述锁相环参数的初始值获取相位偏差的相轨迹的步骤包括：
[0034]
以相变量δ为横坐标、相变量δ
′
为纵坐标构成相平面，且δ
′
为δ的一阶导数；
[0035]
每一所述第一初始条件下，相变量δ和δ
′
随着参变量t的推移在相平面上运动形成一个第一相轨迹；以及
[0036]
将收敛的第一相轨迹对应的所述锁相环的比例参数和积分参数的第一初始值作为所述锁相环参数的终值。
[0037]
优选地，根据阻尼参数和调整时间参数计算所述锁相环参数的初始值的步骤还包括：
[0038]
根据所述阻尼参数和所述调整时间参数计算所述锁相环的比例参数和积分参数的基准值；
[0039]
将所述锁相环积分参数的基准值设置为积分参数的第二初始值且将比例参数的基准值的0.1～2倍对应设置为比例参数的第二初始值；或者将所述锁相环比例参数的基准值设置为比例参数的第二初始值且将积分参数的基准值的0.1～2倍对应设置为积分参数的第二初始值；以及
[0040]
将所述锁相环的比例参数和积分参数的第二初始值作为第二初始条件。
[0041]
优选地，基于所述相位偏差二阶系统模型，采用相平面法和所述锁相环参数的初始值获取相位偏差的相轨迹的步骤还包括：
[0042]
每一所述第二初始条件下，相变量δ和δ
′
随着参变量t的推移在相平面上运动形成一个第二相轨迹；
[0043]
将收敛的第二相轨迹和未收敛的第二相轨迹的分界点所对应的所述锁相环的比例参数和积分参数的第二初始值作为并联换流器的同步稳定临界点。
[0044]
优选地，所述并联换流器锁相环参数计算方法，还包括：采用时域仿真分析对所述锁相环参数的终值进行验证；若时域仿真结果满足预设要求，则应用所述锁相环参数的终值。
[0045]
本发明与现有技术相比至少具有以下优点之一：
[0046]
本发明提供的并联换流器锁相环参数计算方法，可以根据锁相环输出相位与电网电压相位的偏差构建相位偏差二阶系统模型，以及根据阻尼参数和调整时间参数计算锁相环参数的初始值，并基于相位偏差二阶系统模型，采用相平面法和锁相环参数的初始值获取相位偏差的相轨迹，从而得到锁相环参数的终值及并联换流器的同步稳定临界点。
[0047]
本发明采用相平面法绘制相位偏差二阶系统模型的相轨迹，可以直观、准确地反映锁相环的比例参数和积分参数对并联换流器的同步稳定性的影响，并可以根据相轨迹的收敛情况确定锁相环的比例参数和积分参数的终值。
[0048]
本发明还可以根据相轨迹的收敛情况确定并联换流器的暂态同步稳定的临界点，并能对暂态同步稳定性有一个较为全面的了解。
[0049]
本发明还能显著降低锁相环参数确定过程中的计算量，且所确定的锁相环的比例参数和积分参数的终值能保证并联换流器的同步稳定性。
附图说明
[0050]
图1是本发明一实施例提供的一种并联换流器锁相环参数计算方法的流程图；
[0051]
图2是本发明一实施例提供的一种并联换流器锁相环参数计算方法中锁相环的结构示意图；
[0052]
图3是本发明一实施例提供的一种并联换流器锁相环参数计算方法中并联换流器接入电网的单线图；
[0053]
图4是本发明一实施例提供的一种并联换流器锁相环参数计算方法中设置不同阻尼参数的相平面仿真结果；
[0054]
图5是本发明一实施例提供的一种并联换流器锁相环参数计算方法中设置不同调整时间参数的相平面仿真结果；
[0055]
图6是本发明一实施例提供的一种并联换流器锁相环参数计算方法中调整时间参数为0.2s且阻尼参数为0.1的时域仿真结果；
[0056]
图7是本发明一实施例提供的一种并联换流器锁相环参数计算方法中调整时间参数为0.2s且阻尼参数为0.5的时域仿真结果；
[0057]
图8是本发明一实施例提供的一种并联换流器锁相环参数计算方法中调整时间参数为0.2s且阻尼参数为1.0的时域仿真结果；
[0058]
图9是本发明一实施例提供的一种并联换流器锁相环参数计算方法中阻尼参数为0.4且调整时间参数为0.1s的时域仿真结果；
[0059]
图10是本发明一实施例提供的一种并联换流器锁相环参数计算方法中阻尼参数为0.4且调整时间参数为0.5s的时域仿真结果；
[0060]
图11是本发明一实施例提供的一种并联换流器锁相环参数计算方法中阻尼参数为0.4且调整时间参数为1.0s的时域仿真结果；
[0061]
图12是本发明一实施例提供的一种并联换流器锁相环参数计算方法中改变锁相环比例参数的相平面仿真结果；
[0062]
图13是本发明一实施例提供的一种并联换流器锁相环参数计算方法中改变锁相环积分参数的相平面仿真结果。
具体实施方式
[0063]
以下结合附图和具体实施方式对本发明提出的一种并联换流器锁相环参数计算方法做进一步详细说明。根据下面说明，本发明的优点和特征将更清楚。需要说明的是，附图采用非常简化的形式且均使用非精准的比例，仅用以方便、明晰地辅助说明本发明实施方式的目的。为了使本发明的目的、特征和优点能够更加明显易懂，请参阅附图。须知，本说明书所附图式所绘示的结构、比例、大小等，均仅用以配合说明书所揭示的内容，以供熟悉此技术的人士了解与阅读，并非用以限定本发明实施的限定条件，故不具技术上的实质意义，任何结构的修饰、比例关系的改变或大小的调整，在不影响本发明所能产生的功效及所能达成的目的下，均应仍落在本发明所揭示的技术内容能涵盖的范围内。
[0064]
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者还包括这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
[0065]
结合附图1～13所示，本实施例提供一种并联换流器锁相环参数计算方法，包括：步骤s110、计算锁相环输出相位与电网电压相位的偏差；步骤s120、根据所述锁相环输出相位与所述电网电压相位的偏差，构建相位偏差二阶系统模型；步骤s130、根据阻尼参数和调整时间参数计算所述锁相环参数的初始值；以及步骤s140、基于所述相位偏差二阶系统模
型，采用相平面法和所述锁相环参数的初始值获取相位偏差的相轨迹，以得到所述锁相环参数的终值及并联换流器的同步稳定临界点。
[0066]
可以理解的是，在一些其他的实施例中，所述锁相环参数包括比例参数和积分参数。
[0067]
具体的，在本实施例中，所述并联换流器锁相环是一种同步旋转轴锁相环(synchronous reference frame phase-locked loop，srf-pll)，其典型结构如图2所示。从图2中可以看出，并网点(即接入点)的三相电压v
pcc
通过abc到dq的坐标变换，可得到v
pccd
和v
pccq
，其中v
pccd
为并网点电压v
pcc
的d轴分量，v
pccq
为并网点电压v
pcc
的q轴分量。所述锁相环pll的参数有三个，分别为比例参数k
ppll
、积分参数k
ipll
和电网角频率基准ωg；其中电网角频率基准ωg通常是固定不变的，对于工频为50hz的电网而言，电网角频率基准ωg为2
×
50
×
π，近似为314.159265。因此，所述锁相环参数中需要设计和整定的只有比例参数k
ppll
和积分参数k
ipll
，但本发明不以此为限。
[0068]
具体的，在本实施例中，所述并联换流器接入所述电网的单线图如图3所示，其中所述并联换流器可以用受控电流源来等效，所述电网可以用无穷大电源等效，电网输电线路可以用电感等效；θ
ic
为所述并联换流器输出电流相对于参考坐标系的相位，且θ
ic
＝θi+θ
pll
；θi为所述并联换流器输出电流相位与所述锁相环pll输出相位θ
pll
的偏差，对于新能源机组中的并联换流器，一般都按单位功率因数运行，此时θi＝0；zg代表所述电网阻抗，在本实施例中只考虑所述电网阻抗zg的电感部分，即仅考虑所述电网电感lg，但本发明不以此为限。
[0069]
请同时参考图1、图2和图3，所述步骤s110中，所述锁相环输出相位采用如下公式表示：
[0070]
θ
pll
＝∫(k
ppllvpccq
+k
ipll
∫v
pccq
dt+ωg)dt
ꢀꢀꢀꢀꢀꢀꢀ
(1)
[0071]
其中，θ
pll
为所述锁相环输出相位；k
ppll
为所述锁相环的比例参数；k
ipll
为所述锁相环的积分参数；v
pcc
为并网点电压；v
pccq
为并网点电压的q轴分量；ωg为所述电网角频率基准；t为参变量；
[0072]
所述电网电压相位采用如下公式表示：
[0073]
θg＝∫ωgdt
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0074]
其中，θg为所述电网电压相对于参考坐标系的相位；
[0075]
所述锁相环输出相位与所述电网电压相位的偏差采用如下公式进行计算：
[0076]
δ＝θ
pll-θg[0077]
＝∫(k
ppllvpccq
+k
ipll
∫v
pccq
dt)dt
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0078]
其中，δ为所述锁相环输出相位与所述电网电压相位的偏差。
[0079]
可以理解的是，在一些其他的实施例中，所述并网点电压v
pcc
的q轴分量采用如下公式进行计算：
[0080]vpccq
＝vgsin(θ
g-θ
pll
)+lg·icm
·
ωgꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)
[0081]
其中，vg为所述电网电压幅值；lg为所述电网电感；i
cm
为并联换流器内部电流指令幅值。
[0082]
请同时参考图1、图2和图3，所述步骤s120中，所述相位偏差二阶系统模型采用如下公式表示：
[0083][0084]
具体的，在本实施例中，对公式(3)两边求取二阶导数，可得：
[0085][0086]
将v
pccq
的表达式即公式(4)代入公式(6)中，然后通过一系列的化简和整理，则可得到所述相位偏差二阶系统模型的表达式即公式(5)。
[0087]
在其他的实施例中，还可以取状态变量(相变量)δ＝x1，δ
′
＝x2，且δ
′
为δ的一阶导数，则所述相位偏差二阶系统模型的状态方程如下：
[0088][0089]
公式(7)则与公式(5)等价，皆可用于绘制相轨迹，但本发明不以此为限。
[0090]
请同时参考图1、图2和图3，所述步骤s130中，所述锁相环的比例参数和积分参数采用如下公式进行计算：
[0091][0092][0093]
其中，ts表示调整时间参数；ζ表示阻尼参数。
[0094]
具体的，在本实施例中，根据公式(5)或(7)所表示的所述相位偏差二阶系统模型，可以利用所述阻尼参数和所述调整时间参数来描述动态过程；且所述阻尼参数和所述调整时间参数的计算公式如下：
[0095][0096][0097]
通过对所述调整时间参数的计算公式(11)变形即可得到所述锁相环的比例参数的计算公式(8)；通过对所述阻尼参数的计算公式(10)变形即可得到所述锁相环的积分参数的计算公式(9)，但本发明不以此为限。
[0098]
请同时参考图1、图4和图5，所述步骤s130包括：保持所述阻尼参数和所述调整时间参数中一个值不变，另一个设置不同值，对应计算所述锁相环的比例参数和积分参数的第一初始值；以及将所述锁相环的比例参数和积分参数的第一初始值构成第一初始条件。
[0099]
可以理解的是，在一些其他的实施例中，所述步骤s140包括：以相变量δ为横坐标、相变量δ
′
为纵坐标构成相平面，且δ
′
为δ的一阶导数；每一所述第一初始条件下，相变量δ和
δ
′
随着参变量t的推移在相平面上运动形成一个第一相轨迹；以及将收敛的第一相轨迹对应的所述锁相环的比例参数和积分参数的第一初始值作为所述锁相环参数的终值。
[0100]
具体的，相平面法(phase plane method)是一种基于时域的分析方法，该方法通过图解法将一阶和二阶系统的运动过程转化为位置和速度平面上的相轨迹，从而比较直观、准确地反映系统的稳定性、平衡状态和稳态精度，以及初始条件及参数对系统运动的影响。由于在相平面上只能表示两个独立的变量，故相平面法只适用于一二阶线性和非线性系统。更具体的，所述相位偏差二阶系统模型为二阶线性方程，故本实施例中可以利用相平面法分析所述锁相环的比例参数和积分参数对所述并联换流器的同步稳定性的影响，但本发明不以此为限。
[0101]
在本实施例中，根据微分方程解的存在与唯一性定理，对于任一所述第一初始条件，公式(5)或公式(7)在相平面上都有一条相轨迹与之对应。多个所述初始条件下的运动则对应多条相轨迹，从而形成所述相轨迹簇。通过这些相轨迹或相轨迹簇可以优化所述锁相环的比例参数和积分参数，并判断所述并联换流器的同步稳定性。
[0102]
具体的，在本实施例中，可以保持所述调整时间参数不变(例如ts＝0.2s)，设置不同的所述阻尼参数，并得到相应的相平面分析结果(如图4所示)。从图4中可以看出，所述阻尼参数ζ＝0.1对应的第一相轨迹是不收敛的，即所述并联换流器的暂态同步是不稳定的；所述阻尼参数ζ＝0.5和ζ＝1.0对应的第一相轨迹是收敛的，即所述并联换流器的暂态同步是稳定的；可知，所述阻尼参数值略大应该是有利于所述并联换流器同步稳定性的。此时，可以选取所述调整时间参数ts＝0.2s，所述阻尼参数ζ＝0.5或ζ＝1.0时计算得到的所述锁相环的比例参数值和积分参数值为所述锁相环参数的终值，但本发明不以此为限。
[0103]
此外，还可以保持所述阻尼参数不变(ζ＝0.4)，设置不同所述调整时间参数，并得到相应的相平面分析结果(如图5所示)。从图5中可以看出，所述调整时间参数ts＝0.1s对应的第一相轨迹是不收敛的，即所述并联换流器的暂态同步是不稳定的；所述调整时间参数ts＝0.5s和ts＝1.0s时对应的第一相轨迹是收敛的，即所述并联换流器的暂态同步是稳定的；可知，所述调整时间参数值略大应该是有利于所述并联换流器同步稳定性的。此时，可以选取所述阻尼参数ζ＝0.4，所述调整时间参数ts＝0.5s或ts＝1.0s时计算得到的所述锁相环的比例参数值和积分参数值为所述锁相环参数的终值，但本发明不以此为限。
[0104]
请同时参考图1、图6至图11，所述并联换流器锁相环参数计算方法，还包括：采用时域仿真分析对所述锁相环参数的终值进行验证；若时域仿真结果满足预设要求，则应用所述锁相环参数的终值。
[0105]
具体的，在本实施例中，图6、图7和图8分别给出了保持所述调整时间参数ts＝0.2s不变，所述阻尼参数ζ＝0.1、ζ＝0.5和ζ＝1.0时的暂态同步稳定性时域仿真分析中所述锁相环输出相位θ
pll
的仿真结果。可以看出，时域仿真结果和相平面法的仿真结果是相互符合的，即所述阻尼参数ζ＝0.1时所述并联换流器的暂态同步是不稳定的，所述阻尼参数ζ＝0.5和ζ＝1.0时所述并联换流器的暂态同步是稳定的，进一步验证了所述阻尼参数值较大是有利于所述并联换流器同步稳定性的，同时可以将所述调整时间参数ts＝0.2s，所述阻尼参数ζ＝0.5或ζ＝1.0时计算得到的所述锁相环的比例参数值和积分参数值应用于所述锁相环中。优选地，可以将一组、二组或三组所述锁相环的比例参数值和积分参数值应用于所述锁相环中，以保证所述并联换流器的同步稳定性；且不宜采用更多组的所述锁相环
的比例参数值和积分参数值，但本发明不以此为限。
[0106]
图9、图10和图11分别给出了保持所述阻尼参数ζ＝0.4不变，所述调整时间参数ts＝0.1s、ts＝0.5s和ts＝1.0s时的暂态同步稳定性时域仿真分析中所述锁相环输相位θ
pll
的仿真结果。可以看出，时域仿真结果和相平面法的仿真结果是相互符合的，即所述调整时间参数ts＝0.1s时所述并联换流器的暂态同步是不稳定的，所述调整时间参数ts＝0.5s和ts＝1.0s时所述并联换流器的暂态同步是稳定的，进一步验证了所述调整时间参数较大是有利于所述并联换流器同步稳定性的，同时可以将所述阻尼参数ζ＝0.4，所述调整时间参数ts＝0.5s或ts＝1.0s时计算得到的所述锁相环的比例参数值和积分参数值应用于所述锁相环中，但本发明不以此为限。
[0107]
请同时参考图1、图12和图13，所述步骤s130还包括：根据所述阻尼参数和所述调整时间参数计算所述锁相环的比例参数和积分参数的基准值；将所述锁相环积分参数的基准值设置为积分参数的第二初始值且将比例参数的基准值的0.1～2倍对应设置为比例参数的第二初始值；或者将所述锁相环比例参数的基准值设置为比例参数的第二初始值且将积分参数的基准值的0.1～2倍对应设置为积分参数的第二初始值；以及将所述锁相环的比例参数和积分参数的第二初始值作为第二初始条件。
[0108]
所述步骤s140还包括：每一所述第二初始条件下，相变量δ和δ
′
随着参变量t的推移在相平面上运动形成一个第二相轨迹；将收敛的第二相轨迹和未收敛的第二相轨迹的分界点所对应的所述锁相环的比例参数和积分参数的第二初始值作为所述并联换流器的同步稳定临界点。
[0109]
具体的，在本实施例中，可以保持所述锁相环的积分参数k
ipll
不变而改变比例参数k
ppll
，或者保持所述锁相环的比例参数k
ppll
不变而改变积分参数k
ipll
，并做出一系列的渐变线，从而明确所述锁相环的比例参数k
ppll
和积分参数k
ipll
对所述并联换流器暂态同步稳定性的影响程度，进而寻找临界点。更具体的，可以以所述调整时间参数ts＝0.2s和所述阻尼参数ζ＝0.5为基准，计算得到所述锁相环的比例参数基准值和积分参数基准值
[0110]
若保持所述锁相环的积分参数基准值不变，则可以绘制所述锁相环的比例参数为对应的第二相轨迹，并得到相应的第二相轨迹簇(如图12所示)；从图12中可以看出，根据所述第二相轨迹的收敛情况，所述并联换流器的同步稳定临界点为所述调整时间参数ts＝0.286s和所述阻尼参数ζ＝0.35时计算得到的所述锁相环的比例参数值和积分参数值。若保持所述锁相环的比例参数基准值不变，则可以绘制所述锁相环的积分参数为对应的第二相轨迹，并得到相应的第二相轨迹簇(如图13所示)；从图13中可以看出，根据所述第二相轨迹的收敛情况，所述并联换流器的同步稳定临界点为所述调整时间参数ts＝0.2s和所述阻尼参数ζ＝0.42时计算得到的所述锁相环的比例参数值和积分参数值。此外，从图12和图13的相轨迹簇可以看出，增大所述锁相环的比例参数k
ppll
值，所述阻尼参数值增大，有利于所述并联换流器暂态同步稳定；减少所述锁相环的积分参数k
ipll
，所述阻尼参数值增大，有利于所述并联换流器暂态同步稳定，但本发明不以此为限。
[0111]
综上所述，本实施例提供一种并联换流器锁相环参数计算方法，可以根据锁相环输出相位与电网电压相位的偏差，构建相位偏差二阶系统模型；根据阻尼参数和调整时间参数计算锁相环参数的初始值；以及基于相位偏差二阶系统模型，采用相平面法和锁相环参数的初始值获取相位偏差的相轨迹，以得到锁相环参数的终值及并联换流器的同步稳定临界点。本实施例采用相平面法绘制相位偏差二阶系统模型的相轨迹，可以直观、准确地反映锁相环的比例参数和积分参数对并联换流器的同步稳定性的影响，并可以根据相轨迹的收敛情况确定锁相环的比例参数和积分参数的终值；还可以根据相轨迹的收敛情况确定并联换流器的暂态同步稳定的临界点，并能对暂态同步稳定性有一个较为全面的了解。此外，还能显著降低锁相环参数确定过程中的计算量，且所确定的锁相环的比例参数和积分参数的终值能保证并联换流器的同步稳定性。
[0112]
尽管本发明的内容已经通过上述优选实施例作了详细介绍，但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后，对于本发明的多种修改和替代都将是显而易见的。因此，本发明的保护范围应由所附的权利要求来限定。