CN106712055A - 一种与低励限制功能相协调的电力系统稳定器配置方法 - Google Patents

Abstract

本发明涉及一种与低励限制功能相协调的电力系统稳定器配置方法，其特征在于，包括：设置基于低励限制边界线的发电机典型进相运行工况点集合；确定低励限制未动作时传统电力系统稳定器PSS₁经过励磁环节产生的电磁转矩增量信号相对于PSS₁输入信号的开环传递函数G_X(s)；确定低励限制动作时新增电力系统稳定器PSS₂经过励磁环节产生的电磁转矩增量信号相对于PSS₂输入信号的开环传递函数G′_X(s)；利用传递函数G_X(s)和G′_X(s)，构建PSS₂的超前‐滞后环节时间常数的优化数学模型并求解。该配置方法具有考虑因素全面，操作简单方便、实用性强等优点，可有效解决低励限制动作时引起的电力系统动态稳性定恶化的问题。

Description

一种与低励限制功能相协调的电力系统稳定器配置方法

技术领域

本发明涉及电力系统领域，具体涉及一种与低励限制功能相协调的电力系统稳定器配置方法。

背景技术

低励限制(UEL)作为励磁系统中一项重要的辅助功能，其作用在于到防止发电机因进相过深而失稳，防止定子铁芯端部过热和降低发电机失磁保护误动的风险。UEL按接入励磁系统方式的不同可划分为叠加型和选择型2种类型，前者目前应用更为广泛。UEL以机组实时有功P_t，无功Q_t和机端电压U_t为输入变量。根据P_t和U_t计算获得的实时的UEL无功边界值Q_VR，将Q_VR与Q_t的差值信号经过一个超前—滞后环节输出至励磁系统电压叠加点。

电力系统低频振荡是振荡频率在0.1～2Hz范围内的电气量的持续振荡现象，在弱联系、远距离、重负荷输电线路上更容易出现。低频振荡产生的主要原因是现代电力系统采用的快速且高放大倍数励磁系统使电力系统原有的正阻尼减弱甚至出现负阻尼。电力系统稳定器(PSS)已经被证实为一种抑制低频振荡的有效手段，并且已经广泛应用于现代电力系统中。PSS的基本原理是引入母线频率变化信号、加速功率信号、轴滑差值信号、电磁功率差值信号中的一种或几种，通过控制环节补偿励磁系统的相位滞后特性，产生一个附加的阻尼转矩以抑制低频振荡。《电力系统稳定器整定实验导则》(DL/T1231‐2013)中规定了工程上基于相位补偿法确定PSS参数的整定方法。

因为UEL和PSS的输出信号均输出至励磁系统电压叠加点，二者在机组受扰动态过程中均参与励磁调节而相互影响，传统的PSS参数无法起到理想的效果，现场已发生多起因UEL作用使得电力系统动态稳定性恶化的案例。

发明内容

本发明的目的是提供一种与低励限制功能相协调的电力系统稳定器配置方法，构建并求解考虑低励限制动作时新增电力系统稳定器的超前‐滞后环节时间常数的优化数学模型，该优化数学模型的构建避免了低励限制投入对系统动态稳定性的恶化作用。

为实现上述目的，本发明采用以下技术方案：

一种与低励限制功能相协调的电力系统稳定器配置方法，所述方法包括以下步骤：

(1)设置基于低励限制边界线的发电机典型进相运行工况点集合；

(2)考虑低励限制未动作时传统电力系统稳定器PSS₁投入的情况，基于单机‐无穷大系统Heffron‐Phillips模型，确定PSS₁经过励磁环节产生的电磁转矩增量信号相对于PSS输入信号的开环传递函数G_X(s)；

(3)考虑低励限制动作时新增电力系统稳定器PSS₂投入的情况，基于单机‐无穷大系统扩展Heffron‐Phillips模型，确定PSS₂经过励磁环节产生的电磁转矩增量信号相对于PSS₂输入信号的开环传递函数G′_X(s)；

(4)利用传递函数G_X(s)和G′_X(s)，构建PSS₂的超前‐滞后环节时间常数的优化数学模型并求解。

优选的，所述步骤(1)中设置基于低励限制边界线的发电机典型进相运行工况点集合包括2个工况点；

针对发电机通过联系阻抗jX_e接入无穷大容量的母线构成的单机‐无穷大系统，以发电机额定视在功率S_N和额定电压U_N为基准值，定义工况点1的有功功率P₁、无功功率Q₁、机端电压U₁和联系电抗X_e1为：

P₁＝0.5pu，Q₁＝‐0.1+Q_VR1pu，U₁＝1.0pu，X_e1＝0.4pu；

其中，Q_VR1为根据已知的低励限制边界线函数，Q_VR＝f(P,U)，将P₁和U₁带入后计算的值；

定义工况点2的有功功率P₂、无功功率Q₂、机端电压U₂和联系电抗X_e2为：P₂＝1pu，Q₂＝‐0.1+Q_VR2pu，U₂＝1.0pu，X_e2＝0.2pu；

其中，Q_VR2为根据已知的低励限制边界线函数，Q_VR＝f(P,U)，将P₂和U₂带入后计算的值。

优选的，所述步骤(2)中考虑低励限制未动作时传统电力系统稳定器PSS₁投入的情况，基于单机‐无穷大系统Heffron‐Phillips模型，确定PSS₁经过励磁环节产生的电磁转矩增量信号相对于PSS输入信号的开环传递函数G_X(s)为：

G_X(s)＝K₂G(s)G_PSS1(s) (1)

式(1)中，G(s)＝K₃K_A/(1+K₃K₆K_A+sK₃KT′_d0)，G_PSS1(s)为PSS₁的传递函数，K_A为励磁系统主环放大倍数，T′_d0为直轴暂态开路时间常数，K、K₂、K₃和K₆皆为Heffron‐Phillips模型中根据系统工况和发电机参数计算获得的系数。

优选的，所述步骤(3)中考虑低励限制动作时新增电力系统稳定器PSS₂投入的情况，基于单机‐无穷大系统扩展Heffron‐Phillips模型，确定PSS₂经过励磁环节产生的电磁转矩增量信号相对于PSS₂输入信号的开环传递函数G′_X(s)为：

G′_X(s)＝K₂G(s)G_PSS2(s)/(1-G(s)G_UEL(s)K₈) (2)

式(2)中，G(s)＝K₃K_A/(1+K₃K₆K_A+sK₃KT′_d0)；G_UEL(s)为UEL的传递函数，G_PSS2(s)为PSS₂的传递函数，K、K₂、K₃、K₆和K₈皆为为扩展Heffron‐Phillips模型中根据系统工况和发电机参数计算获得的系数。

优选的，所述步骤(4)中利用传递函数G_X(s)和G′_X(s)，构建PSS₂的超前‐滞后环节时间常数的优化数学模型和求解方法包括如下步骤：

4‐1确定PSS₂的超前滞后环节的传递函数；

4‐2构造计算时间常数T₁′、T₂′、T₃′和T₄′的优化数学模型；

4‐3采用粒子群优化算法求解优化数学模型。

进一步地，所述步骤(4‐1)包括：确定PPSS₂的超前滞后环节的传递函数为(1+T₁'s)/(1+T₂'s)*(1+T₃'s)/(1+T₄'s)，其中，T₁′、T₂′、T₃′、T₄′为时间常数。

进一步地，所述步骤(4‐2)中构建计算时间常数T₁′、T₂′、T₃′和T₄′的优化数学模型如下式所示：

式(3)中，T_j′为PSS₂的超前滞后环节时间常数，f_i为振荡频率点，取值范围为0.1～2Hz；为典型进相运行工况点1时，考虑低励限制未启动时传统电力系统稳定器PSS₁作用下，开环传递函数G_X(s)在0.1～2Hz频率范围内对应的相位值序列；为典型进相运行工况点1时，考虑低励限制动作时新增电力系统稳定器PSS₂作用下，开环传递函数G′_X(s)在0.1～2Hz频率范围内对应的相位值序列；为典型进相运行工况点2时，考虑低励限制未启动时传统电力系统稳定器PSS₁作用下，开环传递函数G_X(s)在0.1～2Hz频率范围内对应的相位值序列；为典型进相运行工况点2时，考虑低励限制启动时新增电力系统稳定器PSS₂作用下，开环传递函数G′_X(s)在0.1～2Hz频率范围内对应的相位值序列；

f_mk为PSS超前滞后环节的中心频率，PSS₁对应的2个中心频率分别为：和PSS₂对应的2个中心频率分别为：和

和最接近的现有技术比，本发明提供技术方案具有以下优异效果：

1、本发明提供的技术方案通过发电机典型进相运行工况点集合的设置，考虑了低励限制动作时新增电力系统稳定器对发电机大负荷方式和小负荷方式不同运行工况的适应性，保证在运行方式发生变化时仍发挥理想的效果；

2、本发明提供的技术方案考虑了新增电力系统稳定器对整个低频振荡频段的有效覆盖，对不同振荡模式均可提供有效阻尼；

3、本发明提供的技术方案具有操作简单方便，实用性强的优点。

附图说明

图1为本发明中所提出技术方案的流程图；

图2为不包含UEL的单机无穷大系统Heffron‐Phillips模型；

图3为本发明中包含UEL的单机无穷大系统扩展Heffron‐Phillips模型；

图4为本发明中工况1条件下UEL动作前后PSS₁提供电磁转矩相频特性对比图；

图5为本发明中工况2条件下UEL动作前后PSS₁提供电磁转矩相频特性对比图；

图6为本发明中工况1条件下PSS₁和PSS₂提供电磁转矩的相频特性对比图；

图7为本发明中工况2条件下PSS₁和PSS₂提供电磁转矩的相频特性对比图。

具体实施方式

下面结合具体实例对本发明进行详细说明。

如图1所示，本发明提出一种与低励限制功能相协调的电力系统稳定器配置方法，包括：

(1)设置基于低励限制边界线的发电机典型进相运行工况点集合；

(2)考虑低励限制未动作时传统电力系统稳定器PSS₁投入的情况，基于单机‐无穷大系统Heffron‐Phillips模型，确定PSS₁经过励磁环节产生的电磁转矩增量信号相对于PSS输入信号的开环传递函数G_X(s)；

(3)考虑低励限制动作时新增电力系统稳定器PSS₂投入的情况，基于单机‐无穷大系统扩展Heffron‐Phillips模型，确定PSS₂经过励磁环节产生的电磁转矩增量信号相对于PSS₂输入信号的开环传递函数G′_X(s)；

(4)利用传递函数G_X(s)和G′_X(s)，构建PSS₂的超前‐滞后环节时间常数的优化数学模型并求解。

具体通过以下具体实施例中的计算均按照该流程执行：

图2所示为基于K₁～K₆系数的单机无穷大系统Heffron-Phillips模型，该模型的推导过程及其模型系数K₁～K₆表达式见文献《电力系统稳定性及发电机励磁控制》(刘取著，《中国电力出版社》出版，2007年)相关章节；

图3所示为基于K₁～K₈系数且含UEL的单机无穷大系统Heffron-Phillips模型，该模型是在原Heffron-Phillips模型的基础上考虑UEL动作时的调节作用推导获得的，这里给出图2中扩展Heffron-Phillips模型新增系数K₇和K₈的表达式如下：

式中U_d0和U_q0为稳态运行点的机端电压直轴和交轴分量；I_d0和I_q0为扰动前机端电流直轴和交轴分量；x′_d为发电机直轴暂态电抗；x_q为发电机交轴电抗；x_e为系统联系电抗；K_u和C分别为为低励限制线斜率和截距。

令G₃(s)＝K₃/(1+sK₃T′_d0)，G_AVR(s)＝K_A/(1+sT_E)，并记图1中虚线框内传递函数为G(s)，则：

电力系统中广泛采用以静止励磁系统为代表的快速励磁系统，该类励磁系统时间常数一般小于0.1秒，可取T_E＝0，则：

在基于K₁～K₆系数的不包含UEL的单机无穷大系统Heffron-Phillips模型中电磁转矩信号相对于PSS₁输入信号的开环传递函数为：

G_X(s)＝K₂G(s)G_PSS1(s) (8)

式中G_PSS1(s)为PSS₁的传递函数，G(s)＝K₃K_A/(1+K₃K₆K_A+sK₃KT′_d0)，K_A为励磁系统主环放大倍数，T′_d0为直轴暂态开路时间常数，其他参数为Heffron-Phillips模型中根据系统工况和发电机参数计算获得的系数。

在基于K₁～K₈系数的包含UEL的单机无穷大系统扩展Heffron-Phillips模型中电磁转矩信号相对于PSS输入信号的开环传递函数为：

G′_X(s)＝K₂G(s)G_PSS2(s)/(1-G(s)G_UEL(s)K₈) (9)式中G(s)＝K₃K_A/(1+K₃K₆K_A+sK₃KT′_d0)；G_UEL(s)为UEL的传递函数，G_PSS2(s)为PSS₂的传递函数，K_A为励磁系统主环放大倍数，T′_d0为直轴暂态开路时间常数，其他参数为扩展Heffron-Phillips模型中根据系统工况和发电机参数计算获得的系数。

以某典型单机无穷大系统作为算例，该系统基本参数如表1所示。

表1发电机参数表

项目	值
		发电机额定视在功率SN(MVA)	639.5
发电机额定有功功率PN(MW)	550
		发电机额定无功功率QN(MVar)	326.35
发电机直轴同步电抗xd(pu.)	1.627
		发电机直轴暂态电抗x′d(pu.)	0.523
发电机交轴同步电抗xq(pu.)	1.16
		发电机转子惯性时间常数TJ(s)	8.4
发电机直轴暂态开路时间常数T′d0(s)	10
		发电机阻尼系数D(pu.)	2.0

发电机励磁系统参数为：K_A＝39.575；T_E＝0.01s；UEL控制环节参数为：K_H＝10；T_H1＝0.2s；T_H2＝2s；设UEL不动作时的PSS₁和UEL动作时的PSS₂均采用为以电磁功率信号为输入的PSS1A型，其传递函数表达式为：G_PSS(s)＝sK_QST_Q(1+sT₁)(1+sT₃)/(1+sT_Q)/(1+sT₂)/(1+sT₄)；已知PSS₁的参数为：K_QS＝10；T_QS＝6s；T_Q＝6s；T₁＝0.15s；T₂＝0.04s；T₃＝0.5s；T₄＝0.04s。PSS₂的参数为：K′_QS＝10；T′_QS＝6s；T′_Q＝6s；超前滞后环节时间常数T′₁、T′₂、T′₃、T′₄为待求量。

计算步骤(1)中两种工况的扩展Heffron-Phillips模型参数分别为：

工况1：P＝0.5S_N(319.75MW)，Q＝-0.1S_N+Q_VR|_P＝0.5pu＝-0.3S_N(-191.8Mvar)，U＝1.0pu，X_e＝0.2pu；K₁＝0.7092；K₂＝0.7217；K₃＝0.3957；K₄＝0.7967；K₅＝0.0039；K₆＝0.2460；K₇＝0.0054；K₈＝-1.1963；(U_d0＝0.4450；U_q0＝0.8635；I_d0＝-0.0252；I_q0＝0.3836；x_e＝0.4pu)；

工况2：P＝1.0S_N(639.5MW)，Q＝-0.1S_N+Q_VR|_P＝1.0pu＝-0.1S_N(-63.95Mvar)，U＝1.0pu，X_e＝0.4pu；K₁＝0.7891；K₂＝0.9135；K₃＝0.4554；K₄＝1.0085；K₅＝-0.0761；K₆＝0.3360；K₇＝0.0602；K₈＝-0.9395；(U_d0＝0.6270；U_q0＝0.7687；I_d0＝0.3568；I_q0＝0.0.5405；x_e＝0.2pu)。

将已求得模型参数代入式(8)和式(9)(以G_PSS1(s)代替式中G_PSS2(s))可得工况1和工况2下UEL是否动作时PSS₁提供电磁转矩的相频特性分别如图3和图4所示。

根据图3和图4可见，PSS₁在UEL不动作时能够为系统提供合适的阻尼，但UEL动作时PSS₁提供电磁转矩的相位较UEL不动作时超前很多，需整定PSS₂；为使PSS₂在UEL动作后尽可能达到PSS₁在UEL动作之前的控制效果，首先应使PSS₂和PSS₁所提供转矩的相位尽可能接近。

利用传递函数G_X(s)和G′_X(s)，构造确定PSS₂的超前-滞后环节时间常数的优化数学模型和求解方法，包含如下步骤：

1)PSS₁超前滞后环节传递函数为：(1+T₁s)/(1+T₂s)*(1+T₃s)/(1+T₄s)；PSS₂的超前滞后环节的传递函数为：(1+T₁'s)/(1+T₂'s)*(1+T₃'s)/(1+T₄'s)，时间常数T₁、T₂、T₃、T₄为已知量，时间常数T₁′、T₂′、T₃′、T₄′为待求量。

2)构造计算时间常数T₁′、T₂′、T₃′、T₄′的优化数学模型如下：

式(10)中T_j′为PSS₂的超前滞后环节时间常数。f_i为振荡频率点，取值范围为0.1～2Hz。为典型进相运行工况点1时，考虑低励限制未动作时传统电力系统稳定器PSS₁作用下，开环传递函数G_X(s)在0.1～2Hz频率范围内对应的相位值序列；为典型进相运行工况点1时，考虑低励限制动作时新增电力系统稳定器PSS₂作用下，开环传递函数G′_X(s)在0.1～2Hz频率范围内对应的相位值序列。为典型进相运行工况点2时，考虑低励限制未动作时传统电力系统稳定器PSS₁作用下，开环传递函数G_X(s)在0.1～2Hz频率范围内对应的相位值序列；为典型进相运行工况点2时，考虑低励限制动作时新增电力系统稳定器PSS₂作用下，开环传递函数G′_X(s)在0.1～2Hz频率范围内对应的相位值序列。

式(3)中f_mk为PSS超前滞后环节的中心频率，PSS₁对应有2个中心频率为：PSS₂对应有2个中心频率为：

3)采用粒子群优化算法，求解式(10)中的优化问题，设置粒子群规模为PSOsize＝50，微粒的维数为PSOxvs＝4，加速因子为c₁＝c₂＝2.05；惯性常数w＝0.5；设置粒子的初始位置为PSS₁的参数，最大迭代次数为MAXiter＝50；最大迭代误差为MAXerr＝0.001；交叉概率为Pcros＝0.2。

计算结果为：T′₁＝0.3525；T′₂＝0.01s；T′₃＝0.2825；T′₄＝9.0399；将所得参数代入式(8)和(9)中计算可得工况1和工况2下到PSS₂在UEL动作时提供电磁转矩和PSS₁在UEL不动作时提供电磁转矩的相频特性对比图如图5和图6所示，可见PSS₂相位补偿环节时间常数满足相位补偿要求。

最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，所属领域的普通技术人员尽管参照上述实施例应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，这些未脱离本发明精神和范围的任何修改或者等同替换，均在申请待批的本发明的权利要求保护范围之内。

Claims (7)

1.一种与低励限制功能相协调的电力系统稳定器配置方法，其特征在于，所述方法包括以下步骤：

(1)设置基于低励限制边界线的发电机典型进相运行工况点集合；

(2)考虑低励限制未动作时传统电力系统稳定器PSS₁投入的情况，基于单机‐无穷大系统Heffron‐Phillips模型，确定PSS₁经过励磁环节产生的电磁转矩增量信号相对于PSS输入信号的开环传递函数G_X(s)；

(3)考虑低励限制动作时新增电力系统稳定器PSS₂投入的情况，基于单机‐无穷大系统扩展Heffron‐Phillips模型，确定PSS₂经过励磁环节产生的电磁转矩增量信号相对于PSS₂输入信号的开环传递函数G′_X(s)；

(4)利用传递函数G_X(s)和G′_X(s)，构建PSS₂的超前‐滞后环节时间常数的优化数学模型并求解。

2.如权利要求1所述的一种与低励限制功能相协调的电力系统稳定器配置方法，其特征在于，所述步骤(1)中设置基于低励限制边界线的发电机典型进相运行工况点集合包括2个工况点；

针对发电机通过联系阻抗jX_e接入无穷大容量的母线构成的单机‐无穷大系统，以发电机额定视在功率S_N和额定电压U_N为基准值，定义工况点1的有功功率P₁、无功功率Q₁、机端电压U₁和联系电抗X_e1为：

P₁＝0.5pu，Q₁＝‐0.1+Q_VR1pu，U₁＝1.0pu，X_e1＝0.4pu；

其中，Q_VR1为根据已知的低励限制边界线函数，Q_VR＝f(P,U)，将P₁和U₁带入后计算的值；

定义工况点2的有功功率P₂、无功功率Q₂、机端电压U₂和联系电抗X_e2为：P₂＝1pu，Q₂＝‐0.1+Q_VR2pu，U₂＝1.0pu，X_e2＝0.2pu；

其中，Q_VR2为根据已知的低励限制边界线函数，Q_VR＝f(P,U)，将P₂和U₂带入后计算的值。

3.如权利要求1所述的一种与低励限制功能相协调的电力系统稳定器配置方法，其特征在于：所述步骤(2)中考虑低励限制未动作时传统电力系统稳定器PSS₁投入的情况，基于单机‐无穷大系统Heffron‐Phillips模型，确定PSS₁经过励磁环节产生的电磁转矩增量信号相对于PSS输入信号的开环传递函数G_X(s)为：

$G_X\left(s\right)=K_{2}G\left(s\right)G_{{\mathrm{PSS}}_1}\left(s\right)---\left(1\right)$

式(1)中，G(s)＝K₃K_A/(1+K₃K₆K_A+sK₃KT′_d0)，G_PSS1(s)为PSS₁的传递函数，K_A为励磁系统主环放大倍数，T′_d0为直轴暂态开路时间常数，K、K₂、K₃和K₆皆为Heffron‐Phillips模型中根据系统工况和发电机参数计算获得的系数。

4.如权利要求1所述的一种与低励限制功能相协调的电力系统稳定器配置方法，其特征在于，所述步骤(3)中考虑低励限制动作时新增电力系统稳定器PSS₂投入的情况，基于单机‐无穷大系统扩展Heffron‐Phillips模型，确定PSS₂经过励磁环节产生的电磁转矩增量信号相对于PSS₂输入信号的开环传递函数G′_X(s)为：

$G_X^{\′}\left(s\right)=K_{2}G\left(s\right)G_{{\mathrm{PSS}}_2}\left(s\right)/(1-G(s\left)G_{UEL}\right(s\left)K_8\right)---\left(2\right)$

式(2)中，G(s)＝K₃K_A/(1+K₃K₆K_A+sK₃KT′_d0)；G_UEL(s)为UEL的传递函数，G_PSS2(s)为PSS₂的传递函数，K、K₂、K₃、K₆和K₈皆为扩展Heffron‐Phillips模型中根据系统工况和发电机参数计算获得的系数。

5.如权利要求1所述的一种与低励限制功能相协调的电力系统稳定器配置方法，其特征在于，所述步骤(4)中利用传递函数G_X(s)和G′_X(s)，构建PSS₂的超前‐滞后环节时间常数的优化数学模型和求解方法包括如下步骤：

4‐1确定PSS₂的超前滞后环节的传递函数；

4‐2构造计算时间常数T₁′、T₂′、T₃′和T₄′的优化数学模型；

4‐3采用粒子群优化算法求解优化数学模型。

6.如权利要求5所述的一种与低励限制功能相协调的电力系统稳定器配置方法，其特征在于，所述步骤(4‐1)包括：确定PSS₂的超前滞后环节的传递函数为(1+T₁'s)/(1+T₂'s)*(1+T₃'s)/(1+T₄'s),其中，T₁′、T₂′、T₃′、T₄′为时间常数。

7.如权利要求5所述的一种与低励限制功能相协调的电力系统稳定器配置方法，其特征在于，所述步骤(4‐2)中构建计算时间常数T₁′、T₂′、T₃′和T₄′的优化数学模型如下式所示：

式(3)中，T_j′为PSS₂的超前滞后环节时间常数，f_i为振荡频率点，取值范围为0.1～2Hz； 为典型进相运行工况点1时，考虑低励限制未启动时传统电力系统稳定器PSS₁作用下，开环传递函数G_X(s)在0.1～2Hz频率范围内对应的相位值序列；为典型进相运行工况点1时，考虑低励限制动作时新增电力系统稳定器PSS₂作用下，开环传递函数G′_X(s)在0.1～2Hz频率范围内对应的相位值序列； 为典型进相运行工况点2时，考虑低励限制未启动时传统电力系统稳定器PSS₁作用下，开环传递函数G_X(s)在0.1～2Hz频率范围内对应的相位值序列；为典型进相运行工况点2时，考虑低励限制启动时新增电力系统稳定器PSS₂作用下，开环传递函数G′_X(s)在0.1～2Hz频率范围内对应的相位值序列；

f_mk为PSS超前滞后环节的中心频率，PSS₁对应的2个中心频率分别为： $f_{m1}=1/\left(2\mathrm{\π}\sqrt{T_1{T}_2}\right)$ 和 $f_{m2}=1/\left(2\mathrm{\π}\sqrt{T_3{T}_4}\right);$ PSS₂对应的2个中心频率分别为： $f_{m3}=1/\left(2\mathrm{\π}\sqrt{T_1^{\′}{T}_2^{\′}}\right)$ 和 $f_{m4}=1/\left(2\mathrm{\π}\sqrt{T_3^{\′}{T}_4^{\′}}\right).$