CN101119095A - 基于非线性鲁棒电力系统稳定器的大扰动实时仿真系统 - Google Patents

Abstract

基于非线性鲁棒电力系统稳定器的大扰动实时仿真系统属于电力系统稳定性仿真实验技术。其特征在于，该仿真系统由功率放大器、非线性鲁棒电力系统稳定器NR－PSS和电力系统实时数字仿真器RTDS依闭环串接而成，输向RTDS中发电机励磁控制同可控硅的合成信号是由NR－PSS产生的，其中会有，AVR的控制信号以及考虑了由发电机不确定性引起的大扰动信号的影响的辅助励磁控制信号，从而解决了用动态模拟大扰动实验由于机端参数和动模实验、现场实验时系统差异较大的问题，使得本发明能真实反应NR－PSS在世纪运行的系统中发生大扰动时其各方面的性能。

Description

基于非线性鲁棒电力系统稳定器的大扰动实时仿真系统

技术领域

本属于电力系统稳定控制技术和仿真实验方法领域。

背景技术

非线性电力系统稳定器(NR-PSS)装置的各项性能可以通过动模实验和现场试验加以检验。但是由于现场机组设备和系统运行安全的条件限制，只能检验系统发生小扰动的工况下NR-PSS的性能；而动模实验虽然可以模拟系统发生大扰动的工况，但是由于机组的参数和动模实验系统结构与现场实际运行系统差异较大，因而无法真实反应NR-PSS在实际运行系统中发生大扰动时其各方面性能。而上述种种局限可以通过采用RTDS大扰动仿真实验加以克服。

RTDS(Real Time Digital Simulator)实时数字仿真器是由加拿大RTDS技术公司研制的数字实时仿真装置，有计算软件(RSCAD)、计算处理和接口等硬件设备(Rack)组成，包括配套的工作站或微机。RTDS基于多处理器并行计算和高速DSP芯片计算电力系统各元件的数学模型，因而能够连续和实时地模拟电力系统的电磁暂态和机电暂态现象。由于RTDS能够维持实时条件下的连续运行，实际的控制保护设备就可以连接到RTDS进行闭环实验以分析和研究控制保护设备的性能。

发明内容

本发明专利的基本特征在于利用RTDS进行了NR-PSS装置的数字一物理混合的大扰动仿真实验。其中实物部分为NR-PSS装置；数学模型包括水轮机、调速器、发电机、变压器、线路、断路器、各种量测元件以及无穷大电力系统等部分。实验项目主要有三相短路故障暂态稳定极限传输功率测试实验、三相短路故障极限切除时间、稳定极限传输功率实验、无功功率进相极限测试、抑制各种频率的低频振荡等。

本发明的特征在于，依次含有以下步骤：

1、本系统含有一个由功率放大器、非线性鲁棒电力系统稳定器NR-PSS以及电力系统实时数字仿真器RTDS依次串连组成的闭环系统，其中：

功率放大器向所选NR-PSS输出发电机i的三相电压和三相电流，分别用u_a，u_b，u_c和i_a，i_b，i_c表示；

NR-PSS有数字集成电路构成，含有：模数转换电路、状态参数计算电路，所述NR-PSS输出的辅助励磁控制信号与V_fiN ^R _-PSS的计算电路、自动调压器AVR输出的主励磁控制信号V_fiAVR的计算电路以及输往发电机励磁的可控硅控制端的合成励磁控制信号V_fi的计算电路，其中：

模数转换电路，对发电机三相电压u_a，u_b，u_c以及三相电流i_a，i_b，i_c作模数转换；

状态参数计算电路，输入信号是所选魔术转换电路的输出信号，还有以下设定的发电机参数：给定频率值ω₀＝314.16，D轴电抗X_di、Q轴电抗X_qi、D轴暂态电抗X_di′、输出信号是标幺值下的Q轴电势E_qi、Q轴暂态电势E_qi′、D轴电流i_qi、Q轴电流i_qi、频率偏差Δω_i、功角偏差Δδ_i、有功功率偏差ΔP_ei以及线电压偏差ΔV_ei、微分量

和

，所述标幺制下的各状态参数：线电压和线电流的额定值、发电机组视在功率、稳态频率值ω₀以及额定线电压、额定线电流下的阻抗值的标幺值在[0，1]选取，其中各状态参数分别按照下式计算：

$E_{\mathrm{qi}}=\sqrt{{(V_{\mathrm{ti}}+\frac{Q_{\mathrm{ei}}{x}_{\mathrm{di}}}{V_{\mathrm{ti}}})}^2+{\left(\frac{P_{\mathrm{ei}}{x}_{\mathrm{di}}}{V_{\mathrm{ti}}}\right)}^2},$

$E_{\mathrm{qi}}^{\′}=\sqrt{{(V_{\mathrm{ti}}+\frac{Q_{\mathrm{ei}}{x}_{\mathrm{di}}^{\′}}{V_{\mathrm{ti}}})}^2+{\left(\frac{P_{\mathrm{ei}}{x}_{\mathrm{di}}^{\′}}{V_{\mathrm{ti}}}\right)}^2},$

$i_{\mathrm{di}}=\frac{c_i}{d_i}{E}_{\mathrm{qi}}^{\′}{I}_{\mathrm{ti}}^2,$

$i_{\mathrm{qi}}=\frac{b_i}{\mathrm{di}}{E}_{\mathrm{qi}}^{\′}{I}_{\mathrm{ti}}^2,$

其中：

$\left\{\begin{array}{c}{a}_i=Q_{\mathrm{ei}}+x_{\mathrm{di}}^{\′}{I}_{\mathrm{ti}}^2\\ b_i=P_{\mathrm{ei}}\\ c_i=Q_{\mathrm{ei}}+x_{\mathrm{qi}}{I}_{\mathrm{ti}}^2\\ d_i=a_i{c}_i+b_i^2\end{array}\right.,$

Δω_i＝(ω_i-ω_0i)，ω_i为发电机转速反馈量，

$\mathrm{\Δ}{\mathrm{\δ}}_i={\∫}_0^t\mathrm{\Δ}{\mathrm{\ω}}_i\mathrm{d\τ},$ t为积分周期，dτ为积分间隔，均为设定值，

ΔP_ei＝(P_ei-P_ei0)，P_ei由额定功率下相电压、相电流有效值算出，P_ei0为视在功率，

ΔV_ti＝(V_ti-V_ti0)，V_ti0为额定线电压，V_ti为线电压，

和

分别由如下所述线性拟合求导法求出：选取当前时刻的微分量i_qi或i_di，连同前面七个时刻的数值组成的向量I_q或I_d得到，

V_fiNR-PSS计算电路：

输入信号是：E_qi、i_qi、E_qi′、

i_qi、Δδ_i、Δω_i、ΔP_ei还有设定量：X_di、X_qi、X_di′、机组转动惯量H_i、反馈系数k_1i，k_2i，k_3i，各自在[0，300]间取值，阻尼调节系数C_1i，C_1i＝1/m_i，m_i在[0.1，10]间取值，轻载修正系数k_qi，以及空载修正系数C_qi，

当i_qi＝i_qimin，或i_qi＝-i_qimin时，k_qi＝0，C_qi＝0.1，

当i_qi＞i_qimin，或i_qi＜-i_qimin时，k_qi＝1，C_qi＝1，

输出信号V_fiNR-PSS按照下式求出：

$\begin{array}{c}{V}_{\mathrm{fiNR}-\mathrm{PSS}}=E_{\mathrm{qi}}-\frac{T_{d0i}}{i_{\mathrm{qi}}}[K_{\mathrm{qi}}{E}_{\mathrm{qi}}^{\′}{\stackrel{\·}{i}}_{\mathrm{qi}}+(x_{\mathrm{qi}}-x_{\mathrm{di}}^{\′})(i_{\mathrm{qi}}{\stackrel{\·}{i}}_{\mathrm{di}}+K_{\mathrm{qi}}{i}_{\mathrm{di}}{\stackrel{\·}{i}}_{\mathrm{qi}})]\\ +C_{\mathrm{qi}}{C}_{\mathrm{li}}\frac{{\mathrm{HT}}_{d0t}}{{\mathrm{\ω}}_0{i}_{\mathrm{qi}}}(k_{1i}\mathrm{\Δ}{\mathrm{\δ}}_i+k_{2i}\mathrm{\Δ}{\mathrm{\ω}}_i-k_{3i}\frac{{\mathrm{\ω}}_0}{H_i}\mathrm{\Δ}{P}_{\mathrm{ei}}),\end{array}$

V_fiAVR计算电路：

输入信号是ΔV_ti以及比例系数k_pi、积分系数K_Ii及微分系数k_di，k_pi在[1，300]间取值，k_Ii、k_di、在[0，100]间取值；

V_fiAVR按下式求出：

$V_{\mathrm{fiAVR}}=-k_{\mathrm{pi}}{\mathrm{\ΔV}}_{\mathrm{ti}}-k_{\mathrm{ii}}{\∫}_0^t\mathrm{\Δ}{V}_{\mathrm{ti}}\mathrm{d\τ}-k_{\mathrm{di}}\mathrm{\Δ}{\stackrel{\·}{V}}_{\mathrm{ti}},$ 其中：

$\mathrm{\Δ}{\stackrel{\·}{V}}_{\mathrm{ti}}=({\stackrel{\·}{V}}_{\mathrm{ti}}-{\stackrel{\·}{V}}_{t0}),$

V_fi计算电路，V_fi按下式算出：

V_fi＝C_3i·V_fiAVR+C_2i·V_fiNR-pss(C₁)，

其中：C_2i为NR-PSS的增益系数，取值范围为[0.1，1]，C_3i为AVR的增益系数，取值范围为[0.5，1]；

电力系统实时数值仿真器RTDS，含有：自并励励磁发电机、位于发电机输出母线上的电压互感器、电流互感器、变压器以及断路器，还含有位于输电线路上的多个断路器；其中：电压互感器和电流互感器输出分别向所述功率放大器输出所述三相电压和三相电流，自并励励磁发电机励磁线圈上串联的可控硅的控制信号是V_fiNR-PSS。

2、所选RTDS是加拿大RTDS技术公司的产品。

3、所得到的V_fiNR-PSS信号需用一个隔直环节

处理，以使Vfi_NR-PSS的变化量减小，T_wi一般取2～16秒。

4、所述T_wi为8秒。

5、所得到的V_fi需同一个一阶惯性环节处理，以滤除掉随机噪声引起的V_fiNR-PSS的输出抖动，T_di为0.03秒。

本发明基于RTDS提出了NR-PSS大扰动工况下的实时仿真方法，它具有以下优点：

1)所提出的方法与物理动模仿真相比，能够准确地模拟实际系统的元件设备和网络参数以及各种运行方式，从而检验NR-PSS各种工况下的性能。

2)所提出的方法能够简便有效地考核系统发生各种大扰动的工况下NR-PSS对于改善系统暂态安全稳定性的性能。

3)所提出的方法能够检验在各种极限工况下NR-PSS提高系统的静态稳定极限传输功率、暂态稳定极限传输功率和各种短路故障的极限切除时间。

4)所提出的方法能够检验当系统发生各种频率的低频振荡时，NR-PSS抑制振荡的性能。

5)所提出的方法能够检验NR-PSS对于系统中参数误差较大或者网络结构发生变化的鲁棒性，并对其各方面性能进行比较分析。

附图说明

图1NR-PSS实现方案

图2NR-PSS的RTDS大扰动实验结构图

图3NR-PSS的RTDS大扰动实验流程图

图4NR-PSS无补偿频率特性测量方法1

图5NR-PSS无补偿频率特性测量方法2

图6NR-PSS无补偿频率特性测量方法3

图7NR-PSS无补偿频率特性测量方法4

具体实施方式

NR-PSS的基本思路是建立电力系统励磁控制的鲁棒非线性动态模型的，采用反馈线性化方法得到系统的精确线性化模型，然后应用线性H_∞控制理论设计其鲁棒控制律，最后代回到设计的非线性反馈律中得到原系统的第i台发电机的NR-PSS控制律为

$V_{\mathrm{fiNR}-\mathrm{PSS}}=E_{\mathrm{qi}}-\frac{T_{d0i}^{\′}}{i_{\mathrm{qi}}}[E_{\mathrm{qi}}^{\′}{\stackrel{\·}{i}}_{\mathrm{qi}}+(x_{\mathrm{qi}}-x_{\mathrm{di}}^{\′})(i_{\mathrm{qi}}{\stackrel{\·}{i}}_{\mathrm{di}}+i_{\mathrm{di}}{\stackrel{\·}{i}}_{\mathrm{qi}})]C_1\frac{T_{\mathrm{ji}}{T}_{d0i}^{\′}}{{\mathrm{\ω}}_0{i}_{\mathrm{qi}}}(k_{1i}:\mathrm{\Δ\δ}+k_{2i}\mathrm{\Δ\ω}-k_{3i}\frac{{\mathrm{\ω}}_0}{T_j}{\mathrm{\ΔP}}_e)---\left(1\right)$

其中，下标i代表第i台发电机的参数(以下同)，i_di和i_qi分别为电枢电流的d轴和q轴分量；ω₀＝314.1529为稳态频率值；Δδ是转子运行角偏差(弧度)；Δω是角速度偏差值(弧度/秒)；ΔP_e是电磁功率偏差(标幺值)；E_q′、E_q为同步机暂态电势和空载电势(标幺值)；V_fiNR-PSS是控制器输出(标幺值)；x_d，x_x，x_d′是同步电抗和暂态电抗(标幺值)；T_d0′为定子开路时励磁绕组时间常数(秒)；T_j是转动惯量(秒)；C₁为阻尼调节系数。

考虑到励磁调节器还需要保证发电机机端电压调节精度，因而将控制律(1)与自动电压器(AVR)采用并联方式进行配合，即NR-PSS的输出与AVR的输出叠加，具体方案见图2。最终的励磁控制规律V_fi为：

V_fi＝C₃·V_fiAVR+C²·V_fiNR-PSS(C₁)    (2)

其中

$V_{\mathrm{fiNR}-\mathrm{PSS}}=E_{\mathrm{qi}}-\frac{T_{d0i}^{\′}}{i_{\mathrm{qi}}}[E_{\mathrm{qi}}^{\′}+(x_{\mathrm{qi}}-x_{\mathrm{di}}^{\′})(i_{\mathrm{qi}}{\stackrel{\·}{i}}_{\mathrm{di}}+i_{\mathrm{di}}{\stackrel{\·}{i}}_{\mathrm{qi}})]+C_i\frac{H_i{T}_{d0i}^{\′}}{{\mathrm{\ω}}_0{i}_{\mathrm{qi}}}(k_{1i}{\mathrm{\Δ\δ}}_i+k_{2i}{\mathrm{\Δ\ω}}_i-k_{3i}\frac{{\mathrm{\ω}}_0}{H_i}{\mathrm{\ΔP}}_{\mathrm{ei}})$

$V_{\mathrm{fiAVR}}=(k_p+k_{D}s)\frac{1}{1+k_{I}s}\mathrm{\Δ}{V}_{\mathrm{ti}}$

ΔV_ti为机端电压的偏差量，k_p、k_I和k_D分别为比例、积分和微分系数。另外，根据实际情况需要，具体工程实现中AVR也可以采用其他形式。NR-PSS控制律可以由DSP芯片实现，具体方案见图2所示。

本发明专利的目的是利用RTDS进行NR-PSS装置的数字-物理混合的大扰动仿真实验，并检验各种极限工况下NR-PSS的性能。它依次含有以下步骤：

步骤1：建立RTDS仿真实验系统。

利用RTDS、NR-PSS装置和功率放大器建立如图1所示的仿真实验系统，其中RTDS建立包括自并励励磁系统、发电机、调速系统、变压器、断路器、继电保护装置、故障控制模型、负荷模型以及无穷大等值电源的电力系统仿真环境，并向NR-PSS装置提供所需要的发电机机端电压、电流和主开关位置等模拟信号和数字信号输出；其中发电机机端电压和电流的信号经过功率放大器放大为与实际系统相符的信号输出到NR-PSS装置；NR-PSS装置采集到这些信号并进行控制律的计算，并输出励磁控制电压V_f1，经过自并励励磁系统的模型后得到发电机的励磁电压。RTDS将励磁电压代入仿真计算中，构成实时的闭环测试系统。实验结构图见附图1。

步骤2：仿真实验系统的初始化。

需要初始化的参数包括RTDS中的发电机组、变压器、线路、调速器、负荷参数，开关量以及继电保护动作时间等。参数可以根据现场实际运行工况，也可基于实际工况设置仿真所需要的工况。

步骤3：NR-PSS参数整定。

基于频域测试的方法，整定NR-PSS的各项参数，使得NR-PSS具有较为优良的性能。包括以下4个步骤：

步骤(3.1)：试验准备以及调整试验工况。

进行NR-PSS的参数整定前需要进行如下的准备：试验通道退出运行，对数字式励磁调节器须设置A/D变换器参数、设置加入点和外部信号的投切控制，检查外部信号的正确输入。设置的A/D变换器参数应有利于提高信噪比。

试验的工况为(以下各项试验均在此工况下进行)：被试验机组有功功率大于0.9额定功率，无功工况为0～0.2额定无功；调频和AGC(自动发电机控制)功能暂时退出；试验机组在机端有并列运行机组时，宜在并列机组停运下进行，也可在其运行下进行试验，但要求并列机组的无功调差率设为运行值。记录试验机组运行数据(包括机端电压、励磁电压、有功功率、无功功率以及原动机出力)、同厂其他机组的出力、各高压母线电压和各出线有功功率。

步骤(3.2)：测量NR-PSS输出点到机端电压无补偿频率特性曲线。

无补偿频率特性是指在不投入NR-PSS的情况下，输入外加干扰信号，测量发电机机端电压Vg对于外加信号的相频特性，又称励磁系统的滞后特性。

测量NR-PSS输出点到机端电压无补偿频率特性曲线要求励磁调节器具有外加信号入口，迭加到AVR的输入或者输出点上，对数字式励磁调节器需要有A/D变换器。实验具备频谱分析仪或低频正弦信号发生器，具备测量时间常数小于20ms的发电机电压变送器。

根据机组励磁调节器所具有的外加信号入口的不同位置和输入的干扰信号不同类型，测量无补偿频率特性曲线试验方法可分为以下4种方法：

·测量方法1

若励磁调节器干扰信号输入点可以叠加在AVR输出点，输入信号选择频谱分析仪输出的噪声信号，则采用图4所示的测量方法。

图4中干扰信号N是频谱分析仪输出的噪声信号，可选为随机信号或者周期性调频信号，前者使用均匀(uniform)滤波窗，后者使用汉宁(hanning)滤波窗。V_g为发电机机端电压，V_ref为电压给定值。

采用频谱分析仪测量无补偿频率特性的试验步骤如下：

(1)设置频谱分析仪的信号源种类、滤波窗、频率范围为0.2～2.0Hz、测量频率特性、输入端量程等。

(2)试验接线如图1所示。逐步增大噪声信号输出直至发电机机端电压V_g有微小摆动。

(3)测量0.2～2.0Hz的频率特性，观察频率特性曲线形状是否光滑；如果曲线不光滑，调整信号大小再进行测量。

(4)记录频率特性曲线的数据。

·测量方法2

如果励磁调节器干扰信号输入点可以叠加在AVR输出点，输入信号可以选择低频正弦信号采用图2所示的测量方法。

图2中干扰信号S为低频正弦信号发生器输出的正弦信号，V_g为发电机机端电压，V_ref为电压给定值。

采用低频正弦信号发生器和波形记录仪器测量无补偿频率特性的试验步骤如下：

(1)在0.2～2.0Hz之间每隔0.1Hz取一个频点，逐个频率点地进行测量。

(2)试验接线如图5所示。逐步增大信号输出直至发电机机端电压V_g有微小摆动，波形稳定后用波形记录仪记录波形，记录的波形应当基本正弦、光滑。

(3)每个频率点测量结束应将信号调到零，另一点频率信号从零起调。频率不同，信号大小不同。

(4)根据波形计算各个频率下的相位，并记录数据。

·测试方法3

若励磁调节器的外加信号入口叠加在AVR的输入点上，采用频谱分析仪噪声信号为输入干扰信号，则可按照图6所示的接线图测量无补偿频率特性曲线。图6中干扰信号N是频谱分析仪输出的噪声信号，可选为随机信号或者周期性调频信号，前者使用均匀(uniform)滤波窗，后者使用汉宁(hanning)滤波窗。V_g为发电机机端电压，V_ref为电压给定值。

由于需要测量的是NR-PSS输出点到机端电压的机组无补偿滞后特性曲线，因而试验前需要将AVR中的微分环节和积分环节退出，并将AVR中电压比例系数降低为10左右(避免引起振荡)。测量步骤如下：

(1)设置频谱分析仪的信号源种类、滤波窗、频率范围为0.2～3.0Hz、测量频率特性、输入端量程等。

(2)试验接线如图6所示。逐步增大噪声信号输出直至发电机机端电压Vg有微小摆动。

(3)测量0.2～3Hz的频率特性，观察频率特性曲线形状是否光滑；如果曲线不光滑，调整信号大小再进行测量。

(4)记录频率特性曲线的数据。

·测试方法4

若励磁调节器的外加信号入口叠加在AVR的输入点上，采用低频正弦信号为输入干扰信号，则可按照图7所示的接线图测量无补偿频率特性曲线。图7中干扰信号S为低频正弦信号发生器输出的正弦信号，V_g为发电机机端电压，V_ref为电压给定值。

与测试方法3类似，试验前需要将AVR中的微分环节和积分环节退出，并将AVR中电压比例系数降低为10左右(避免引起振荡)。测量步骤如下：

采用低频正弦信号发生器和波形记录仪器测量无补偿频率特性的试验步骤如下：

(1)在0.2～2.0Hz之间每隔0.1Hz取一个频点，逐个频率点地进行测量。

(2)试验接线如图7所示。逐步增大信号输出直至发电机机端电压V_g有微小摆动，波形稳定后用波形记录仪记录波形，记录的波形应当基本正弦、光滑。

(3)每个频率点测量结束应将信号调到零，另一点频率信号从零起调。频率不同，信号大小不同。

(4)根据波形计算各个频率下的相位，并记录数据。

步骤3：计算NR-PSS线性部分的反馈系数。

NR-PSS线性部分可以提取出来为V_l＝L₁Δδ+L₂Δω+L₃ΔP_e。按照《DL/T650-1998大型汽轮发电机自并励静止励磁系统技术条件》附录B2.1和《DL/T843-2003大型汽轮发电机自交流励磁机静止励磁系统技术条件》附录A2.1中要求：在0.2～2Hz低频振荡区应使PSS输出的力矩向量在Δω轴在超前10度～滞后45度以内，并使本机振荡频率力矩向量在Δω轴超前10度～滞后30度间。因此，根据步骤一测量得到的无补偿滞后特性曲线可以计算NR-PSS需要补偿的角度。

根据发电机的经典二阶模型，考虑原动机机械功率变化量ΔP_m等于0，可以计算出对于不同频率的输入信号，NR-PSS的线性部分相对于Δω轴的超前滞后角度如式(2)

${\mathrm{\α}}_l=\arctan (-\frac{L_1{\mathrm{\ω}}_0-4{\mathrm{\π}}^2{f}^2{T}_j{L}_3}{2\mathrm{\πf}{L}_2+2\mathrm{\πfD}{L}_3})---\left(2\right)$

其中f为输入信号的频率，ω₀为额定转速，T_j为发电机转动惯量，D是机组阻尼系数。

将式(2)计算得到的NR-PSS每个频率点的超前滞后角度α_l加上步骤1测量的无补偿频率特性的角度，可以得到每个频率点下NR-PSS输出的力矩相对于Δω轴的角度。通过调整参数L₁、L₂和L₃，可使得机组在NR-PSS线性部分补偿后的超前滞后角度满足上述标准要求。参数L₁、L₂和L₃的取值范围为[0，100]。

在验证NR-PSS的线性部分提供的力矩满足要求后，根据参数L₁、Lx和L₃以及机组参数，进一步求解相应的NR-PSS的参数K₁、K₂和K₃；即 $K_1=\frac{L_1{\mathrm{\ω}}_0{I}_q}{T_j{T_{d0}}^{\′}},$ $K_2=\frac{L_1{\mathrm{\ω}}_0{i}_q}{T_j{T_{d0}}^{\′}}$ $K_3=\frac{-L_3{i}_q}{{T_{d0}}^{\′}},$ 式中各符号的意义如前所述。

所求得到的K₁、K₂和K₃即可构成Ricatti方程的解P矩阵；通过Ricatti方程即可得到相应的Q矩阵和R矩阵；可以证明参数K₁、K₂和K₃即为Q矩阵和R矩阵所对应的二次型性能指标下的最优解。

步骤4：确定NR-PSS的临界增益。

确定NR-PSS的参数K₁、K₂和K₃之后，需要确定进一步NR-PSS的阻尼系数C₁和NR-PSS增益系数C₂，使闭环系统拥有优良的阻尼特性同时保证电压的稳定运行。确定C₁和C₂两个系数的方法包括临界放大倍数法和负载阶跃实验法。

临界放大倍数法是在投入NR-PSS后，逐步增大NR-PSS的阻尼调节系数C₁，仔细观察发电机励磁电压，当励磁电压开始出现连续不断的振荡时立即退出NR-PSS。此时的放大倍数即为临界增益；而使用增益可取临界增益的1/2～1/3；从而确定NR-PSS阻尼调节系数C₁的值。

负载阶跃法是输入端设置阶跃量为1％～4％额定发电机电压，阶跃后6s左右即反向阶跃回原值。逐步增大NR-PSS的增益系数C₂，并进行有无NR-PSS的负载阶跃实验，观察其电压动态响应和有功的稳定性。当有NR-PSS进行阶跃后有功振荡的次数大于5次，并且进入稳态运行后机端电压和有功功率出现连续不断的振荡，则减少C₂的数值至刚产生持续振荡时的放大倍数，即临界增益。使用增益可取临界增益的1/2～1/3；从而确定NR-PSS反馈C₂的值。

在完成上述NR-PSS整定步骤以后，即可确定NR-PSS的各项参数。

步骤4：进行常规小扰动实验验证NR-PSS参数的正确性。

进行2～5％电压给定值阶跃实验、无功反调实验、切机切线路切负荷等小扰动实验验证NR-PSS性能，并可对NR-PSS的参数进行微调。

步骤5：进行极端运行工况的打扰动实验下考验NR-PSS性能。

可以在系统各种运行工况下进行的大扰动实验包括以下类型：

①给定值±20％的电压阶跃实验：利用NR-PSS装置的电压给定值，在发电机额定负载下进行电压给定值±20％左右的阶跃实验，观察电压波形和有功功率振荡波形，检验NR-PSS能否使有功功率的振荡迅速衰减，恢复系统稳定；

②短路实验：设置RTDS的故障类型以及故障持续时间，主变压器高压侧或线路上进行单相接地短路实验、两相短路实验、两相接地短路实验、三相短路实验以及三相接地短路实验等考核NR-PSS的性能，观察有功功率的波形是否能够迅速恢复稳定；

③静稳极限传输功率测试：通过不断增加系统中所有发电机组的有功功率，观察有功功率波形直到其发生振荡并且振荡幅度逐渐增大至失步，记录下此时系统的总有功功率，即为该工况下的NR-PSS极限传输功率；

④无功功率极限测试：机组有功保持在额定状态，不断减小机组无功功率的输出，观察有功功率波形直到其发生振荡并且振荡幅度逐渐增大至失步，记录下此时系统的总无功功率即为该工况下的NR-PSS极限无功功率；

⑤暂稳极限功率测试：利用RTDS软件的故障控制模块可以设置在主变压器高压侧或者线路上发生0.1s三相接地短路故障，通过不断增加系统中发电机组的有功功率，观察有功功率波形，直到故障切除后有功功率振荡幅度逐渐增大至系统发生失步，记录下此时系统中发电机组的总有功功率，即为该工况下的NR-PSS暂稳极限传输功率；

⑥极限切除时间测试：可以在主变压器高压侧或者线路上发生三相短路故障，通过RTDS软件的故障控制模块不断增加短路持续时间，观察有功功率波形，直到故障切除后有功功率振荡幅度逐渐增大至系统发生失步，记录下此时的故障持续时间即为该工况下的NR-PSS极限切除时间。

另外，可以设置RTDS断路器开关，可以控制其他发电机和线路是否并网，使得仿真系统结构改变，即上述实验可以在系统单机单线或者多机多线路的工况下，进行以上各种极限传输功率实验和三相短路故障的极限切除时间等测试。

进行上述实验时，都要必须记录有功功率、机端电压、无功功率以及励磁电压等波形数据。

步骤6：NR-PSS抑制频率为0.1～3Hz的低频振荡实验

将发电机组的转动惯量常数增大或减小，使得进行电压阶跃时机组的本机振荡频率减小或增大；也可以调整线路参数使得系统振荡频率在0.1～3Hz之间变化。通过电压阶跃或者切机切线路等方法使得系统发生低频振荡，投入NR-PSS后观察其抑制各种频率的情况，并记录有功功率、机端电压、无功功率和励磁电压等波形数据。

步骤7：NR-PSS鲁棒性能检测实验。

改变发电机机组参数以及变压器或线路等电网的参数，使得参数发生了较大的变化，而输入NR-PSS装置中的机组及系统参数不变，使得控制器参数有了较大的误差。在高压母线或者线路上进行三相接地短路等大扰动实验，考验NR-PSS对于发电机组和系统网络参数发生变化时其性能，记录有功功率、机端电压、无功功率和励磁电压等波形数据。比较NR-PSS在参数正确和参数错误的情况下有功功率振荡的次数、幅值以及振荡平息时间，从而验证NR-PSS在系统各种扰动情况下对参数和网络结构变化的鲁棒性。

步骤8：整理和分析记录数据。

将上述仿真实验得到的数据进行分析整理，并得出NR-PSS装置的RTDS大扰动仿真实验的结论。

Claims (5)

1.基于非线性鲁棒电力系统稳定器的大扰动实时仿真系统，其特征在于，含有一个由功率放大器、非线性鲁棒电力系统稳定器NR-PSS以及电力系统实时数字仿真器RTDS依次串连组成的闭环系统，其中：

功率放大器向所选NR-PSS输出发电机i的三相电压和三相电流，分别用u_a，u_b，u_c和i_a，i_b，i_c表示；

NR-PSS有数字集成电路构成，含有：模数转换电路、状态参数计算电路，所述NR-PSS输出的辅助励磁控制信号与V_fiNR-PSS的计算电路、自动调压器AVR输出的主励磁控制信号V_fiAVR的计算电路以及输往发电机励磁的可控硅控制端的合成励磁控制信号V_fi的计算电路，其中：

模数转换电路，对发电机三相电压u_a，u_b，u_c以及三相电流i_a，i_b，i_c作模数转换；

状态参数计算电路，输入信号是所选魔术转换电路的输出信号，还有以下设定的发电机参数：给定频率值ω₀＝314.16，D轴电抗X_di、Q轴电抗X_qi、D轴暂态电抗X_di′、输出信号是标幺值下的Q轴电势E_qi、Q轴暂态电势E_qi′、D轴电流i_di、Q轴电流i_qi、频率偏差Δω_i、功角偏差△δ_i、有功功率偏差ΔP_ei以及线电压偏差ΔV_ei、微分量

和

，所述标幺制下的各状态参数：线电压和线电流的额定值、发电机组视在功率、稳态频率值ω₀。以及额定线电压、额定线电流下的阻抗值的标幺值在[0，1]选取，其中各状态参数分别按照下式计算：

$E_{\mathrm{qi}}=\sqrt{{(V_{\mathrm{ti}}+\frac{Q_{\mathrm{ei}}{x}_{\mathrm{di}}}{V_{\mathrm{ti}}})}^2+{\left(\frac{P_{\mathrm{ei}}{x}_{\mathrm{di}}}{V_{\mathrm{ti}}}\right)}^2},$

$E_{\mathrm{qi}}^{\′}=\sqrt{{(V_{\mathrm{ti}}+\frac{Q_{\mathrm{ei}}{x}_{\mathrm{di}}^{\′}}{V_{\mathrm{ti}}})}^2+{\left(\frac{P_{\mathrm{ei}}{x}_{\mathrm{di}}^{\′}}{V_{\mathrm{ti}}}\right)}^2},$

$i_{\mathrm{di}}=\frac{c_i}{d_i}{E}_{\mathrm{qi}}^{\′}{I}_{\mathrm{ti}}^2,$

$i_{\mathrm{qi}}=\frac{b_i}{\mathrm{di}}{E}_{\mathrm{qi}}^{\′}{I}_{\mathrm{ti}}^2,$

其中：

$\left\{\begin{array}{c}{a}_i=Q_{\mathrm{ei}}+x_{\mathrm{di}}^{\′}{I}_{\mathrm{ti}}^2\\ b_i=P_{\mathrm{ei}}\\ c_i=Q_{\mathrm{ei}}+x_{\mathrm{qi}}{I}_{\mathrm{ti}}^2\\ d_i=a_i{c}_i+b_i^2\end{array}\right.,$

Δω_i＝(ω_i-ω_0i)，ω_i为发电机转速反馈量，

${\mathrm{\Δ\δ}}_i={\∫}_0^t\mathrm{\Δ}{\mathrm{\ω}}_i\mathrm{d\τ},$ t为积分周期，dτ为积分间隔，均为设定值，

ΔP_ei＝(P_ei-P_ei0)，P_ei由额定功率下相电压、相电流有效值算出，P_ei0为视在功率，

ΔV_ti＝(V_ti-V_ti0)，V_ti0为额定线电压，V_ti为线电压，

和

分别由如下所述线性拟合求导法求出：选取当前时刻的微分量i_qi或i_di连同前面七个时刻的数值组成的向量I_q或I_d得到，

V_fiNR-PSS计算电路：

输入信号是：E_qi、i_qi、E_qi′、

i_qi、Δδ_i、Δω_i、ΔP_ei，还有设定量：X_di、X_qi、X_di′、机组转动惯量H_i、反馈系数k_1i，k_2i，k_3i，各自在[0，300]间取值，阻尼调节系数C_1i，C_1i＝1/m_i，m_i，在[0.1，10]间取值，轻载修正系数k_qi，以及空载修正系数C_qi，

当i_qi＝i_qimin，或i_qi＝-i_qimin时，k_qi＝0，C_qi=0.1，

当i_qi＞i_qimin，或i_qi＜-i_qimin时，k_qi＝1，C_qi＝1；

输出信号V_fiNR-PSS按照下式求出：

$V_{\mathrm{fiNR}-\mathrm{PSS}}=E_{\mathrm{qi}}-\frac{T_{d0i}}{i_{\mathrm{qi}}}[K_{\mathrm{qi}}{E}_{\mathrm{qi}}^{\′}{\stackrel{.}{i}}_{\mathrm{qi}}+(x_{\mathrm{qi}}-x_{\mathrm{di}}^{\′})(i_{\mathrm{qi}}{\stackrel{.}{i}}_{\mathrm{di}}+K_{\mathrm{qi}}{i}_{\mathrm{di}}{\stackrel{.}{i}}_{\mathrm{qi}})],$

$+C_{\mathrm{qi}}{C}_{1i}\frac{{\mathrm{HT}}_{d0i}}{{\mathrm{\ω}}_0{i}_{\mathrm{qi}}}(k_{1i}\mathrm{\Δ}{\mathrm{\δ}}_i+k_{2i}\mathrm{\Δ}{\mathrm{\ω}}_i-k_{3i}\frac{{\mathrm{\ω}}_0}{H_i}\mathrm{\Δ}{P}_{\mathrm{ei}})$

V_fiAVR计算电路：

输入信号是ΔV_ti以及比例系数k_pi，、积分系数k_li，及微分系数k_di，k_pi，在[1，300]间取值，k_Ii、k_di在[0，100]间取值；

V_fiAVR按下式求出：

$V_{\mathrm{fiAVR}}=-k_{\mathrm{pi}}\mathrm{\Δ}{V}_{\mathrm{ti}}-k_{\mathrm{ii}}{\∫}_0^t\mathrm{\Δ}{V}_{\mathrm{ti}}\mathrm{d\τ}-k_{\mathrm{di}}\mathrm{\Δ}{\stackrel{.}{V}}_{\mathrm{ti}},$ 其中：

Δ_ti＝(_ti-_ti0)，

V_fi计算电路，V_fi按下式算出：

V_fi＝C_3i·V_fiAVR+C_2i·V_fiNR-PSS(C₁)，

其中：C_2i为NR-PSS的增益系数，取值范围为[0.1，1]，C_3i为AVR的增益系数，取值范围为[0.5，1]；

电力系统实时数值仿真器RTDS，含有：自并励励磁发电机、位于发电机输出母线上的电压互感器、电流互感器、变压器以及断路器，还含有位于输电线路上的多个断路器；其中：电压互感器和电流互感器输出分别向所述功率放大器输出所述三相电压和三相电流，自并励励磁发电机励磁线圈上串联的可控硅的控制信号是V_fiNR-PSS。

2.根据权利要求1所述的基于非线性鲁棒电力系统稳定器的大扰动实时仿真系统，其特征在于，所选RTDS是加拿大RTDS技术公司的产品。

3.根据权利要求1所述的基于非线性鲁棒电力系统稳定器的大扰动实时仿真系统，其特征在于，所得到的V_fiNR-PSS信号需用一个隔直环节

处理，以使V_fiNR-PSS的变化量减小，T_wi一般取2～16秒。

4.根据权利要求1所述的基于非线性鲁棒电力系统稳定器的大扰动实时仿真系统，其特征在于，所述T_wi为8秒。

5.根据权利要求1所述的基于非线性鲁棒电力系统稳定器的大扰动实时仿真系统，其特征在于，所得到的V_fi需同一个一阶惯性环节

处理，以滤除掉随机噪声引起的V_fiNR-PSS的输出抖动，T_di为0.03秒。

Images