SignalCalc在励磁系统静态模型辨识中的应用的制作方法

本发明涉及励磁系统实验领域，具体是signalcalc在励磁系统静态模型辨识中的应用。
背景技术：
：励磁系统静态模型辨识主要是对励磁调节器pid和pss的各个传递函数环节进行频谱辨识，确定各个环节的准确性，具体辨识示意图如图1所示。由图1可以看出，辨识需要对pid及pss（这里为2a模型）共计9~14个传递函数环节进行辨识，每个环节的辨识需要读取记录每个频率点对应的相频和幅频特性，共计64组数据，则整个试验下来，仅读取记录的数据就多达576~896组，而目前的试验仪器安捷伦35670a由于功能的限制，试验数据的读取记录只能通过手动完成，使得原始数据的读取记录就占了50%的试验时间，如图2所示，因此要提高辨识的效率，就必须找到快捷的数据读取记录方法。表1两种试验仪器的优缺点分析试验仪器优点缺点安捷伦35670a已得到较好的应用。便携性差，参数设置复杂；抗干扰能力较差，容易产生坏数据；数据的读取记录需要通过人工手动完成。signalcalc外观小巧，易携带；通过计算机软件控制，可方便地导出扫频数据。主要应用于力学振动试验，无实践于励磁试验的先例；参数设置与试验效果均为未知。表1对安捷伦35670a和signalcalc的优缺点进行了分析，由表1可知，signalcalc在辨识完成后，即可直接导出数据，有效弥补了安捷伦35670a的缺点，进而提高试验效率，但设备无应用辨识的先例，其参数设置与效果均为未知，需要通过摸索尝试来进一步验证确认。为了更加高效准确的完成励磁试验，在国内首次完成三机无刷励磁系统mec7000的仿真建模，三菱mec7000励磁系统现场实测模型如图3所示，在进行bpa仿真时，根据对比验证，采用fm型无刷励磁系统模型作为仿真计算用的等效模型，模型图如图4所示。在模型的等效转换中，图3中的主励磁机及不可控功率整流器模块b可与图4中主励磁机及不可控功率整流器模型d分别对应等效，其中图3中的k2/rfe可由图4的二级调节器增益kb来等效，但由于k2/rfe的现场实测值为14104，超出了二级调节器增益kb的最大取值范围（最大限值为9999），因而无法在bpa数据卡中输入该值；图3中的励磁机电流反馈增益参数kf可由图4中的励磁机电流反馈增益kh进行等效，励磁机电流反馈增益参数kf取值为0.0015，而励磁机电流反馈增益kh输入的数值精度达不到要求（最小限值为0.001），如采用约数，则可能影响仿真效果。因此需要通过进一步的等效计算，在保证仿真效果的情况下，降低二级调节器增益kb，增大励磁机电流反馈增益kh。然而一直以来，在电力系统中，发电机励磁仿真建模校核计算多采用人工读数，手动计算的方式。由于此次是首次对三菱mec7000进行仿真建模，需要通过不断的尝试修正参数，使仿真阶跃响应曲线与现场实测曲线一致，每修正一次，均需要重新进行校核计算，如果采用上述方法，将消耗大量时间在读数与计算上面；另外采用人工读数的方法也有可能对计算精度产生影响，降低仿真校核计算的效率和精度。技术实现要素：为了将signalcalc动态信号分析仪应用于励磁系统静态模型辨识中，发挥其优点，本发明提供了以下技术方案：signalcalc在励磁系统静态模型辨识中的应用，包括以下步骤：（1）采用等效转换的方法，搭建仿真模型；（2）采用bpa仿真工具进行仿真校核，检验搭建的仿真模型的功能效果；（3）将signalcalc应用于励磁系统静态模型辨识中，完成模型环节的辨识工作。进一步，所述步骤（1）包含以下步骤：1）励磁系统现场实测模型进行仿真时采用无刷励磁系统作为等效模型，在等效转换时，分别对励磁系统现场实测模型的主励磁机及不可控功率整流器模型和无刷励磁系统模型的主励磁机及不可控功率整流器模型进行转换；2）建立励磁系统模型转换的等效传递函数，令kb1=k2/rfe，励磁系统现场实测模型的主励磁机及不可控功率整流器模型转换后得到等效传递函数为：（1）无刷励磁系统模型的主励磁机及不可控功率整流器模型转换后得到等效传递函数为：（2）其中：kb为第二级调节器增益，kh为励磁机电流反馈增益，se为励磁机饱和系数，va为励磁调节器输出电压，ve为不可控三相全波整流桥的输出电压，其中k2为励磁调节器二级放大增益，kf为励磁机电流反馈增益，ke为空载励磁电流的倒数，te为励磁机时间常数，rfe为励磁机磁场电阻,kb1为励磁调节器二级放大增益k2与励磁机磁场电阻rfe的比值。进一步，所述步骤2)中等效传递函数（1）和（2），采用解析求解的方法求解，求得第二级调节器增益kb、励磁机电流反馈增益kh的数值；励磁调节器二级放大增益k2取1890，励磁机电流反馈增益kf取0.0015，励磁机时间常数te取2.13,励磁机磁场电阻rfe取0.134，空载励磁电流的倒数ke取0.0121。进一步，所述采用bpa仿真工具进行仿真校核包括以下步骤：1）修正bpa数据卡片，调用bpa程序进行仿真计算；2）将bpa仿真计算结果导入“仿真.xls”；3）启动程序；4）点击“导入仿真数据”选择“仿真.xls”；5）点击“导入实测数据”选择现场实测数据；6）点击“校核计算”，同时计算仿真和实测数据的超调量、上升、峰值、调节时间，并进行误差计算；7）校核仿真结果是否符合要求；8）若仿真结果不符合要求，则返回步骤（1）重新开始进行以上步骤；若仿真结果符合要求，则仿真校核结束。进一步，所述程序为应用编程软件编写的可直接读取数据文件中数据的程序；所述程序开发相应的图形界面，并且读取数据时直接导入到相应的图形界面进行校核计算；所述的数据文件中的数据是通过电力系统分析软件bpa对修正bpa数据卡片进行仿真计算，得出计算结果导出至数据文件得到的；所述的数据文件是excel文件；所述的bpa数据卡片包括励磁系统fm卡、励磁系统f+数据卡、发电机mf卡。进一步，所述的励磁系统fm卡参数：调差系数xc为0、调节器输入滤波器时间常数tr为0.004s、调节器pid增益k为23.55、积分调节选择因子kv为0、电压调节器超前时间常数t1为20s、电压调节器滞后时间常数t2为20s、电压调节器超前时间常数t3为0.95s、电压调节器滞后时间常数t4为3.309s、电压调节器放大器增益ka为1、电压调节器放大器时间常数ta为0.004s、软负反馈放大倍数kf为0、励磁机电流反馈增益kh（标幺值)为0.107；其中，调节器pid增益k、积分调节选择因子kv、电压调节器放大器增益ka、软负反馈放大倍数kf、励磁机电流反馈增益kh均为标幺值。进一步，所述的励磁系统f+数据卡参数：调节器最大内部电压vamax为100、调节器最小内部电压vamin为-100、二级调节器增益kb为169.98、二级调节器时间常数t5为0.01、励磁机自励系数ke为1、励磁机时间常数te为2.13s、最大励磁电压处的励磁机饱和系数se1为0.154、75%最大励磁电压处的励磁机饱和系数se2为0.048、电压调节器最大输出电压vrmax为61.6、电压调节器最小输出电压vrmin为-31.7、换相电抗的整流器负载因子kc为0.279、去磁因子kd为0.22、励磁机励磁电流限制增益kl1为40、励磁机电流限制vl1r为5.6、最大励磁电压efdmax为6；其中，调节器最大内部电压vamax、调节器最小内部电压vamin、励磁机自励系数ke、电压调节器最大输出电压vrmax、电压调节器最小输出电压vrmin、换相电抗的整流器负载因子kc、去磁因子kd、励磁机励磁电流限制增益kl1、励磁机电流限制vl1r、最大励磁电压efdmax均为标幺值。进一步，所述的发电机mf卡参数：发电机母线电压为24.0kv、发电机标幺参数的基准容量为1222.2mva、直轴暂态电抗xd′为0.403、交轴暂态电抗xq′为0.571、直轴不饱和同步电抗xd为1.81、交轴不饱和同步电抗xq为1.74、直轴瞬态开路时间常数tdo′为8.208s、交轴瞬态开路时间常数tqo′为0.8s、定子漏抗为0.234、额定电压时电机饱和系数为0.0784、1.2倍额定电压时电机饱和系数为0.4941；其中直轴暂态电抗xd′、交轴暂态电抗xq′、直轴不饱和同步电抗xd、交轴不饱和同步电抗xq均为标幺值。进一步，所述步骤（3）将signalcalc应用于励磁系统静态模型辨识中，完成模型环节的辨识工作包含以下步骤：1）确定并设置signalcalc控制软件的参数，如以下表6所示：表6signalcalc控制软件的参数参数设置值参数说明fspan10分析频宽lines100谱线数fortf频域分析方式windowhanning窗函数triggerfreerun触发方式activech1,ch2激活的输入通道couplingdcse输入信号的耦合方式waveformrandom输出信号类型level300mv输出信号幅值2）基于signalcalc控制软件，右键软件主界面的单击幅频特性图或相频特性图；3）在弹出的菜单中选择“graph”—“copyactivetrace”；4）拷贝试验数据至excel，直接得到每个频率点对应的相位、幅值。本发明的有益效果为：在未进行模型转换前，由于实测数据超过了bpa数据文件输入的范围，使得仿真工作无法开展。通过模型等效转换以及采用仿真校核计算程序的不断修正计算，确定了所有励磁系统模型参数。signalcalc可采用拷贝的方式直接导出试验数据，避免了安捷伦35670a人工读数计数的弊端，另外其测试数据的可靠性较高，坏数据较少，减少了在数据修正上所耗费的时间。仿真计算效率和精度的提高，可为电网的安全稳定分析提供准确的励磁系统模型，以此可计算电网内部存在的安全风险隐患，防范电网发生低频振荡的风险，从而保障社会居民的用电安全。附图说明图1pid和pss模型频谱辨识示意图；图2采用安捷伦35670a进行辨识试验所需时间统计图；图3三菱mec7000励磁系统现场实测模型图；图4bpa无刷励磁系统模型图；图5三菱mec7000励磁系统现场实测模型中b模块等效图；图6bpa无刷励磁系统模型中d模块等效图；图7是仿真校核计算流程图；图8是仿真校核程序图形界面截图；图9是程序的图形界面截图；图10signalcacl控制软件的参数设置主界面截图；图11成果实施前后对比统计图。具体实施方式下面结合附图对本发明进一步说明：三菱mec7000励磁系统现场实测模型进行bpa仿真时采用fm型无刷励磁系统作为等效模型，在等效转换时，将ke计入模型中，分别对三菱mec7000励磁系统现场实测模型图3的主励磁机及不可控功率整流器模块b和bpa无刷励磁系统模型图4的主励磁机及不可控功率整流器模块d进行转换，分别得到的转换模型图如图5和图6所示。建立励磁系统模型转换图的等效传递函数，令kb1=k2/rfe，励磁系统现场实测模型的主励磁机及不可控功率整流器模型转换后的模型如图3，它的等效传递函数为：（1）无刷励磁系统模型的主励磁机及不可控功率整流器模型转换后的模型如图4，结合传递函数计算的一般方法，可以求得式（2），它的等效传递函数为：（2）其中：kb为第二级调节器增益，kh为励磁机电流反馈增益；se为励磁机饱和系数，va为励磁调节器输出电压，ve为不可控三相全波整流桥的输出电压，其中k2为励磁调节器二级放大增益，kf为励磁机电流反馈增益，te为励磁机时间常数，rfe为励磁机磁场电阻，kb1为励磁调节器二级放大增益k2与励磁机磁场电阻rfe的比值。等效传递函数（1）和（2），采用解析求解的方法求解，可求得第二级调节器增益kb、励磁机电流反馈增益kh的数值；励磁调节器二级放大增益参数取值k2取1890，励磁机电流反馈增益参数取值kf取0.0015,励磁机时间常数参数取值te取2.13,励磁机磁场电阻均为参数取值rfe取0.134。采用解析求解的方法求解等效传递函数，将参数代入式（1），令式（1）等于式（2），即：可求得第二级调节器增益kb、励磁机电流反馈增益kh的数值，如表2所示，表2kb、kh转换前后参数值对比bpa参数转换前的值转换后的值限值kb141041709999（最大值）kh0.00150.10650.001（最小值）由表1可得，模型转换后kb、kh的数值在其取值的限值内，另外由于式（1）等于式（2），可说明参数转换后，两个模型的仿真效果仍然一致的。采用bpa仿真工具进行仿真校核以下步骤：1）在进行校核计算前，先修正bpa数据卡片，调用bpa程序进行仿真计算，将bpa仿真计算结果导出至excel文件；这些数据卡片如表3、表4、表5所示。表3励磁系统fm卡参数表参数名称参数调差系数xc0调节器输入滤波器时间常数tr（秒）0.004调节器pid增益k（标幺值）23.55积分调节选择因子kv（标幺值）0电压调节器超前时间常数t1（秒）20电压调节器滞后时间常数t2（秒）20电压调节器超前时间常数t3（秒）0.95电压调节器滞后时间常数t4（秒）3.309电压调节器放大器增益ka（标幺值）1电压调节器放大器时间常数ta（秒）0.004软负反馈放大倍数kf(标幺值)0励磁机电流反馈增益kh（标幺值）0.107表4励磁系统f+数据卡参数表参数名称参数调节器最大内部电压vamax（标幺值）100调节器最小内部电压vamin（标幺值）-100二级调节器增益kb169.98二级调节器时间常数t50.01励磁机自励系数ke（标幺值）1励磁机时间常数te（秒）2.13最大励磁电压处的励磁机饱和系数se10.15475%最大励磁电压处的励磁机饱和系数se20.048电压调节器最大输出电压vrmax（标幺值）61.6电压调节器最小输出电压vrmin（标幺值）-31.7换相电抗的整流器负载因子kc（标幺值）0.279去磁因子kd（标幺值）0.22励磁机励磁电流限制增益kl1（标幺值）40励磁机电流限制vl1r（标幺值）5.6最大励磁电压efdmax（标幺值）6表5发电机mf卡参数名称参数发电机母线电压（kv）24.0发电机标幺参数的基准容量（mva）1222.2直轴暂态电抗xd′（标幺值）0.403交轴暂态电抗xq′（标幺值）0.571直轴不饱和同步电抗xd（标幺值）1.81交轴不饱和同步电抗xq（标幺值）1.74直轴瞬态开路时间常数tdo′（秒）8.208交轴瞬态开路时间常数tqo′（秒）0.8定子漏抗0.234额定电压时电机饱和系数0.07841.2倍额定电压时电机饱和系数0.49412）启动已经编写好的计算程序，对应着仿真校核程序图形界面如图9点击“导入仿真数据”和“导入实测数据”，最后点击“校核计算”，计算机程序同时计算仿真和实测数据的超调量、上升时间、峰值时间、调节时间并进行误差计算。3）通过仿真校核计算程序的不断修正计算，确定了所有励磁系统模型参数，得到图8为5%发电机电压阶跃仿真曲线与现场实测曲线的对比。4）其中excel文件包括仿真数据文件和实测数据文件。本发明经过不断的尝试验证，确定了signalcalc控制软件的参数，如表6所示，并据此对signalcalc控制软件的参数进行设置。表6signalcalc控制软件的参数参数设置值参数说明fspan10分析频宽lines100谱线数fortf频域分析方式windowhanning窗函数triggerfreerun触发方式activech1,ch2激活的输入通道couplingdcse输入信号的耦合方式waveformrandom输出信号类型level300mv输出信号幅值另外软件主界面右侧的四个显示框分别为输入输出信号的监视窗口、频谱分析的幅频特性图和相频特性图，如图10所示。试验结束后，对试验结果进行分析，不需要逐个记录每个频率点的幅值或相位，右键单击幅频特性图或相频特性图，在弹出的菜单中选择“graph”—“copyactivetrace”，拷贝试验数据至excel，直接得到每个频率点对应的相位、幅值，避免了安捷伦35670a手动逐个记录数据的缺陷，极大提高了工作效率。分别将安捷伦35670a和signalcalc应用于励磁系统的静态模型辨识试验中，完成了对10个模型环节的辨识工作，具体效果如表7所示。表7安捷伦35670a和signalcalc的应用对比注：①幅频特性的计算公式：(p2-p1)/p1。其中p2为signalcalc或安捷伦的实测幅值，p1为对应环节的理论计算幅值。②相频特性的计算公式：q2-q1。其中q2为signalcalc或安捷伦的实测相位值，q1为对应环节的理论计算相位值。③每个环节的辨识时间不包括接线及参数设置的时间。从表7可以看到，signalcalc在相频或幅频特性的辨识精度较安捷伦35670a均有明显提高，而每个环节的辨识时间则缩短了10至15分钟。主要原因是signalcalc可采用拷贝的方式直接导出试验数据，避免了安捷伦35670a人工读数计数的弊端，另外其测试数据的可靠性较高，坏数据较少，减少了在数据修正上所耗费的时间。本发明实施后，总的辨识试验工作时间缩短至45分钟（不包括接线及参数设置的时间），试验效率提高72.7%，提升对比如图11所示。另外本成果可提高仿真校核计算的效率和精度，具体的提升如表8：表8单次仿真和完成建模仿真在校核计算上所耗费时间单次仿真耗时完成建模仿真耗时成果实施前3~5分钟最多75分钟成果实施后不超过30秒最多5分钟效率提升90%93%由表8可知，采用人工读数进行校核计算，每次需耗时3~5分钟，而应用本成果后，每次仅需不到30秒，而一般情况下，完成一次建模仿真需进行8~15次校核计算，则应用本成果最多可节约93%的时间。经济效益的计算依据：试验和仿真效率的提升，可缩短整个励磁试验的时间，整个试验周期由原来的8天缩短至5天，相当于可多承接3/8个电厂试验项目，按每个电厂试验项目30万元计，即可增加潜在收益12万元（算上节省3天的差旅费）。signalcalc应用于核电励磁系统静态模型辨识试验后，试验时间缩短了3个小时，即电厂可以提前3个小时并网发电，由核电上网电价0.215元/千瓦时计算（上网电价0.43元/千瓦时，商业运行之前，按半价计算），可为电厂增加收益69万元；本发明不局限于以上所述的具体实施方式，以上所述仅为本发明的较佳实施案例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。当前第1页1 2 3