CN107220414A - 一种核电机组冷却剂屏蔽主泵动态建模及分析方法 - Google Patents
一种核电机组冷却剂屏蔽主泵动态建模及分析方法 Download PDFInfo
- Publication number
- CN107220414A CN107220414A CN201710311080.6A CN201710311080A CN107220414A CN 107220414 A CN107220414 A CN 107220414A CN 201710311080 A CN201710311080 A CN 201710311080A CN 107220414 A CN107220414 A CN 107220414A
- Authority
- CN
- China
- Prior art keywords
- mrow
- main pump
- cooling agent
- msub
- msubsup
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
- Control Of Positive-Displacement Pumps (AREA)
Abstract
本发明涉及一种核电机组冷却剂屏蔽主泵动态建模及分析方法。针对电网侧电压、频率对主泵的影响以及屏蔽电机本身的特性，建立冷却剂主泵流量计算模型，基于主泵电机电磁转矩公式，考虑冷却剂主泵的摩擦转矩，计算主泵电机整个系统阻力矩，并基于异步电动机运行特性，求解冷却剂主泵电磁转矩，建立用于电力系统分析的冷却剂主泵流量计算模型；根据已有设计及测试数据，进行参数辨识，辨识模型中相应参数。本发明方法考虑了冷却剂主泵的摩擦转矩，能够准确实时反映电网侧电压、频率波动对冷却剂流量影响，进而准确分析电网波动对核电安全影响，弥补目前大多数主泵冷却剂分析围绕主泵惰转和主泵卡死分析，对接入电网的核电分析具有更高的实时性。
Description
技术领域
本发明属于核电接入电网的安全稳定分析技术领域，涉及一种核电机组冷却剂屏蔽主泵动态建模及分析方法。
背景技术
反应堆冷却剂主循环泵是核电站最重要的设备之一，也是整个冷却剂系统中唯一的旋转设备，具有功率大、工作条件苛刻的特点。
正常运行情况下，核电厂的厂用电系统由电网供电，若电网侧出现事故，导致核电厂不能从电网侧获取厂用电时，反应堆会自动安全停堆，但在反应堆安全停堆过程中，需要反应堆冷却系统继续工作，以将堆芯中残余的热量传递出去。以AP1000为代表的第三代先进压水堆核电机组冷却剂主泵采用屏蔽电动机，其飞轮和电动机转子一起包容在RCS的承压部件内，使其径向尺寸受到极大限制，因此，相比于其他堆型，具有更小的转动惯量，冷却剂流量对电网侧波动反映更为明显，因此研究电网侧波动对冷却剂流量造成的影响，针对根据分析结果，进行继电保护整定，停堆等措施提供系统安全可靠性。
目前，常用的冷却剂主泵流量建模方法主要有，流体网格计算方法相关模型、基于四象限特性曲线的相关计算模型等。基于流体网格计算方法的冷却剂主泵流量计算模型主要应用于电站仿真系统中，将冷却剂主泵的工作过程简化为由管路、阀门等连接而成的流体网格系统，更多的是从泵的压力特性入手分析建模，忽略了泵本身的电学特性，并不适用于电力系统分析，也无法反映电网侧波动对冷却剂流量的影响。而基于四象限特性曲线的冷却剂主泵流量计算模型虽然考虑了冷却剂主泵的电学特性，但由于计算繁琐，参数不易获得，一般只用于独立分析主泵特性采用，难以与其他设备模型及电力系统仿真进行结合。
基于此，为了通过仿真研究电网波动对核电机组安全运行的影响，一方面要求冷却剂屏蔽主泵仿真模型能够精确快速地表征冷却剂流量的波动，以便于制定相应控制措施；另一方面，仿真模型应能够接入电力系统仿真软件和其他核电设备模型中进行仿真运算。
发明内容
本发明的目的在于提供一种核电机组冷却剂屏蔽主泵动态建模及分析方法，该方法考虑了冷却剂主泵的摩擦转矩，能够准确实时反映电网侧电压、频率波动对冷却剂流量影响，进而准确分析电网波动对核电安全影响，弥补目前已公开方法对主泵惰转和主泵卡死分析的不足，对接入电网的核电分析具有更高的实时性；同时，模型结构清晰、参数意义明确、易获取，具有很高的准确性和实用性。
为实现上述目的，本发明的技术方案是：一种核电机组冷却剂屏蔽主泵动态建模及分析方法，包括如下步骤，
S1、基于核电机组冷却剂屏蔽主泵特性，建立冷却剂主泵流量计算模型，包括转矩平衡方程、质量流量方程；
S2、基于冷却剂主泵设计参数及测试数据，计算或辨识得到主泵模型参数，并计算主泵电磁转矩、水力转矩和摩擦转矩；
S3、基于步骤S1所建立的冷却剂主泵流量计算模型，仿真计算主泵惰转工况，电压、频率波动工况时冷却剂流量实时变化，与测试曲线作对比，校验模型及参数；
S4、将步骤S1所建立的冷却剂主泵流量计算模型接入核电机组系统模型，并接入电力系统仿真软件，实现核电机组整体动态仿真分析。
在本发明一实施例中，在步骤S1中，所述冷却剂主泵流量计算模型的转矩平衡方程、质量流量方程，是采用以下方法建立的：
冷却剂屏蔽主泵转矩平衡方程为：
式中，Ip为冷却剂主泵的转动惯量，ωp为冷却剂主泵转子角速度，Mpe为冷却剂主泵电磁转矩，Mpm为冷却剂主泵受到的总阻力矩；
主泵电磁转矩Mpe计算表达式为：
式中，p为主泵电动机磁极对数，R2为主泵电动机转子每相绕组电阻，s为主泵电动机转差率，m1为主泵电动机磁极对数，U1n为电网侧系统基准电压，f1n为电网侧系统基准频率，上述各参数均为设计常数；Xσ为主泵电动机转子漏电抗，U1*、f1*分别为实际电网电压标幺值和频率标幺值；K1、K2为与冷却剂主泵相关的常系数；
屏蔽主泵因其飞轮和转子固封于承压部件内，使得其转子转动惯量小于常规泵；转速对转矩变化灵敏度高，摩擦转矩不可忽略，因此屏蔽主泵的阻力矩应同时计及冷却剂流量产生的阻力矩Mp1和由于机械摩擦造成的摩擦阻力矩Mp2，即有
Mpm＝Mp1+Mp2 (3)
冷却剂流量产生的阻力矩Mp1表达式为：
式中，D为冷却剂质量流量，ωp为为冷却剂主泵转子角速度，η为冷却剂主泵效率，H为冷却剂主泵扬程；
Mp2表示主泵机械摩擦引起的阻力矩，采用经验公式表示为
式中，K3为常数，由泵的实验值确定；n取1.7-2.4；
冷却剂质量流量正比于主泵转子角速度，但由于系统惯性，冷却剂流量变化滞后于转子角速度变化，将冷却剂质量流量方程表示为
式中，为转子角速度标幺值，D*为主泵质量流量标幺值，τ为常数，取0.2s。
在本发明一实施例中，在步骤S2中，基于冷却剂主泵设计参数及测试数据，计算或辨识得到主泵模型参数的具体实现如下：
S21、基于冷却剂主泵惰转工况测试数据，采用MATLAB优化设计工具箱进行参数辨识计算，辨识出主泵摩擦阻力矩相关系数K3，n；
S22、基于电网侧频率或电压波动情况下，冷却剂流量测试数据，并将主泵摩擦阻力矩相关系数K3，n作为已知参数，采用MATLAB优化设计工具箱进行参数辨识，得到冷却剂主泵电磁转矩相关系数K1、K2；
S23、基于辨识的参数，建立相应的电磁转矩，摩擦转矩模型方程。
在本发明一实施例中，所述步骤S3的具体实现过程如下：
S31、仿真冷却剂主泵惰转工况下，根据冷却剂流量变化曲线，求解惰转时间常数，与相关模型及参数进行比较分析，验证所述冷却剂主泵流量计算模型的有效性；
S32、仿真电网侧极端事故情况下，根据冷却剂流量与转速变化趋势，验证主泵低转速保护优先级高于冷却剂低流量保护，验证所述冷却剂主泵流量计算模型的合理性；
S33、仿真电网电压、频率阶跃波动工况下，根据冷却剂流量变化，验证冷却剂主泵流量计算模型对电网侧电压、频率扰动响应仿真的准确性。
在本发明一实施例中，在步骤S4中，实现核电机组整体动态仿真分析具体为分析电网侧波动情况下，核电侧主泵流量，燃料温度，蒸汽压力相关参数动态响应过程。
相较于现有技术，本发明具有以下有益效果：
本发明方法考虑了冷却剂主泵的的摩擦转矩，能够准确实时反映电网侧电压、频率波动对冷却剂流量影响，进而准确分析电网波动对核电安全影响，弥补目前已公开方法对主泵惰转和主泵卡死分析的不足，对接入电网的核电分析具有更高的实时性。同时，模型结构清晰、参数意义明确、易获取，具有很高的准确性和实用性。
附图说明
图1是本发明的冷却剂屏蔽主泵动态建模分析流程图。
图2是本发明的核电机组冷却剂屏蔽主泵电动机等效示意图。
图3是摩擦阻力矩相关参数K3、n辨识结果图。
图4是摩擦转矩参数拟合曲线图。
图5是电磁转矩相关参数K1、K2拟合结果图。
图6是冷却剂主泵流量计算模型传输框图。
图7是t＝10s时主泵失去外电源冷却剂流量变化曲线图。
图8是冷却剂主泵电动机供电频率以4Hz/s下降事故时堆芯流量与转速归一化曲线图。
图9是电压/频率阶跃扰动下冷却剂流量变化曲线图。
图10是电网频率下降时，核电机组动态响应图。
具体实施方式
下面结合附图，对本发明的技术方案进行具体说明。
本发明的一种核电机组冷却剂屏蔽主泵动态建模及分析方法，包括如下步骤，
S1、基于核电机组冷却剂屏蔽主泵特性，建立冷却剂主泵流量计算模型，包括转矩平衡方程、质量流量方程；
S2、基于冷却剂主泵设计参数及测试数据，计算或辨识得到主泵模型参数，并计算主泵电磁转矩、水力转矩和摩擦转矩；
S3、基于步骤S1所建立的冷却剂主泵流量计算模型，仿真计算主泵惰转工况，电压、频率波动工况时冷却剂流量实时变化，与测试曲线作对比，校验模型及参数；
S4、将步骤S1所建立的冷却剂主泵流量计算模型接入核电机组系统模型，并接入电力系统仿真软件，实现核电机组整体动态仿真分析。
在步骤S1中，所述冷却剂主泵流量计算模型的转矩平衡方程、质量流量方程，是采用以下方法建立的：
冷却剂屏蔽主泵转矩平衡方程为：
式中，Ip为冷却剂主泵的转动惯量，ωp为冷却剂主泵转子角速度，Mpe为冷却剂主泵电磁转矩，Mpm为冷却剂主泵受到的总阻力矩；
主泵电磁转矩Mpe计算表达式为：
式中，p为主泵电动机磁极对数，R2为主泵电动机转子每相绕组电阻，s为主泵电动机转差率，m1为主泵电动机磁极对数，U1n为电网侧系统基准电压，f1n为电网侧系统基准频率，上述各参数均为设计常数；Xσ为主泵电动机转子漏电抗，U1*、f1*分别为实际电网电压标幺值和频率标幺值；K1、K2为与冷却剂主泵相关的常系数；
屏蔽主泵因其飞轮和转子固封于承压部件内，使得其转子转动惯量小于常规泵；转速对转矩变化灵敏度高，摩擦转矩不可忽略，因此屏蔽主泵的阻力矩应同时计及冷却剂流量产生的阻力矩Mp1和由于机械摩擦造成的摩擦阻力矩Mp2，即有
Mpm＝Mp1+Mp2 (3)
冷却剂流量产生的阻力矩Mp1表达式为：
式中，D为冷却剂质量流量，ωp为冷却剂主泵转子角速度，η为冷却剂主泵效率，H为冷却剂主泵扬程；
Mp2表示主泵机械摩擦引起的阻力矩，采用经验公式表示为
式中，K3为常数，由泵的实验值确定；n取1.7-2.4；
冷却剂质量流量正比于主泵转子角速度，但由于系统惯性，冷却剂流量变化滞后于转子角速度变化，将冷却剂质量流量方程表示为
式中，为转子角速度标幺值，D*为主泵质量流量标幺值，τ为常数，取0.2s。
在步骤S2中，基于冷却剂主泵设计参数及测试数据，计算或辨识得到主泵模型参数的具体实现如下：
S21、基于冷却剂主泵惰转工况测试数据，采用MATLAB优化设计工具箱进行参数辨识计算，辨识出主泵摩擦阻力矩相关系数K3，n；
S22、基于电网侧频率或电压波动情况下，冷却剂流量测试数据，并将主泵摩擦阻力矩相关系数K3，n作为已知参数，采用MATLAB优化设计工具箱进行参数辨识，得到冷却剂主泵电磁转矩相关系数K1、K2；
S23、基于辨识的参数，建立相应的电磁转矩，摩擦转矩模型方程。
所述步骤S3的具体实现过程如下：
S31、仿真冷却剂主泵惰转工况下，根据冷却剂流量变化曲线，求解惰转时间常数，与相关模型及参数进行比较分析，验证所述冷却剂主泵流量计算模型的有效性；
S32、仿真电网侧极端事故情况下，根据冷却剂流量与转速变化趋势，验证主泵低转速保护优先级高于冷却剂低流量保护，验证所述冷却剂主泵流量计算模型的合理性；
S33、仿真电网电压、频率阶跃波动工况下，根据冷却剂流量变化，验证冷却剂主泵流量计算模型对电网侧电压、频率扰动响应仿真的准确性。
在步骤S4中，实现核电机组整体动态仿真分析具体为分析电网侧波动情况下，核电侧主泵流量，燃料温度，蒸汽压力相关参数动态响应过程。
以下为本发明的具体实现过程。
本实施例提供了一种核电机组冷却剂屏蔽主泵动态建模及分析方法，如图1所示，包括以下步骤：
步骤S1：基于核电机组冷却剂屏蔽主泵特性，建立冷却剂主泵流量计算模型，包括转矩平衡方程、质量流量方程。
冷却剂屏蔽主泵转矩平衡方程为：
在(1)式中，Ip为冷却剂主泵的转动惯量，ωp为冷却剂主泵转子角速度，Mpe为冷却剂主泵电磁转矩，Mpm为冷却剂主泵受到的总阻力矩，包含两部分，一部分为冷却剂流量产生的阻力矩Mp1，另一部分为主泵电动机摩擦阻力矩Mp2。
本实施例中，冷却剂屏蔽主泵等效电路如图2所示，主泵电磁转矩Mpe计算表达式为：
式中，p为主泵电动机磁极对数，R2为主泵电动机转子每相绕组电阻，s为主泵电动机转差率，m1为主泵电动机磁极对数，U1n为电网侧系统基准电压，f1n为电网侧系统基准频率，均为设计常数。Xσ为主泵电动机转子漏电抗，U1*、f1*分别为实际电网电压标幺值和频率标幺值。K1、K2为与冷却剂主泵相关的常系数，在本发明中采用已有实验数据进行辨识。
冷却剂流量产生的阻力矩Mp1表达式为：
式中D为冷却剂质量流量，ωp为冷却剂主泵转子角速度，η为主泵效率，H为主泵扬程。
主泵摩擦阻力矩Mp2表达式为：
式中，K3为常数，由泵的实验值确定；n为经验值，一般取1.7-2.4，由实验值确定。
冷却剂质量流量正比于主泵转子角速度，但由于系统惯性，冷却剂流量变化滞后于转子角速度变化，将冷却剂质量流量方程表示为
式中，为转子角速度标幺值，D*为主泵质量流量标幺值，τ为常数，可取值0.2s。
因此冷却剂主泵流量计算模型为：
步骤S2：基于冷却剂主泵设计参数及测试数据，辨识得到主泵模型参数，并计算主泵电磁转矩、水力转矩和摩擦转矩。
冷却剂屏蔽主泵模型需辨识的参数为：主泵电磁转矩相关参数K1，K2，主泵摩擦阻力矩相关参数K3、n；并由此获得主泵电磁转矩，主泵水力阻矩和主泵摩擦转矩方程。本实施例中，对AP1000核电机组冷却剂屏蔽主泵建立的动态模型进行参数识别，并运用辨识出的参数进行仿真分析运算。
S2.1基于冷却剂屏蔽主泵惰转工况，对(4)式进行参数识别。采用采用MATLAB优化设计工具箱，进行参数辨识，辨识方程为：
式(7)中Ip＝931kg·m2,ωp初始值为ωp0＝183.26rad/s，本实施例对摩擦阻力矩参数K3、n拟合结果如图3、4所示，可以看出K3在第2次拟合即开始收敛，收敛值为0.0617，n在第12次拟合开始收敛，收敛值为1.86。
S2.2将步骤S2.1获得的摩擦阻力矩参数K3、n代入公式(1)，得到
本实施例基于已有的冷却剂屏蔽主泵测试数据进行参数辨识，辨识结果如图5所示，在第30次拟合后，K1，K2即开始收敛，收敛值分别为K1＝3059.4，K2＝245.3。
步骤S3：模型仿真与校验
本实施例基于所建立的数学模型，本实施例在Simulink工具箱中搭建仿真模型，如图6所示。通过不同工况下冷却剂流量特性曲线分析验证本模型的有效性和准确性。
S3.1惰转曲线分析。本实施例给出了t＝10s时刻冷却剂屏蔽主泵失去外电源后惰转工况下本发明模型仿真曲线与主泵设计曲线的对比验证，如图7所示。对比验证结果表明本发明提出的模型仿真得到的主泵惰转曲线与设计惰转曲线基本一致，验证了本发明提出模型的有效性。
S3.2电网侧极端事故下的仿真验证。本实施例给出了冷却剂屏蔽主泵供电频率以4Hz/s下降时，根据本发明提出模型仿真获得的标幺化冷却剂流量和主泵转速变化曲线，如图8所示，本实施例仿真结果表明，主泵转速下降和环路流量下降是同步的，同时主泵下降转速超前于冷却剂流量下降，可以根据仿真结果制定主泵低转速停堆整定值的动作优先级高于冷却剂低流量停堆整定。
S3.3频率、电压波动下的仿真验证。本实施例给出了冷却剂屏蔽主泵供电电压阶跃下降到0.95pu时冷却剂流量变化，如图9所示。其中图9(a)为供电电压阶跃下降到0.95pu时冷却剂流量变化曲线，图9(b)为电网频率阶跃下降到0.95pu时冷却剂流量变化曲线。由本实施例仿真结果可以看出频率波动对冷却剂流量影响更为明显。
步骤S4：接入整体模型的仿真。
将搭建的模型接入已有的核电机组，并接入电力系统仿真软件进行全时域仿真，可以实时得到在电网侧故障时，核电侧的冷却剂流量，堆芯燃料温度等输出，进而根据结果作用于堆芯停堆等保护措施制定。
本实施例给出了在t＝2s时电网侧频率阶跃下降到0.95pu后，核电机组相关参数变化曲线，如图10所示。图10(a)为冷却剂流量变化，图10(b)为中子通量变化，图10(c)为蒸汽发生器蒸汽压力变化，图10(d)为冷却剂堆芯燃料温度变化。
从图10中可以看出，当t＝2s时电网侧频率阶跃下降到0.95pu时，冷却剂流量迅速降低，使得冷却剂带着堆芯内部温度减小，堆芯燃料温度上升，同时蒸汽压力降低，反应堆功率增加，中子通量上升。
以上是本发明的较佳实施例，凡依本发明技术方案所作的改变，所产生的功能作用未超出本发明技术方案的范围时，均属于本发明的保护范围。
Claims (5)
1.一种核电机组冷却剂屏蔽主泵动态建模及分析方法，其特征在于：包括如下步骤，
S1、基于核电机组冷却剂屏蔽主泵特性，建立冷却剂主泵流量计算模型，包括转矩平衡方程、质量流量方程；
S2、基于冷却剂主泵设计参数及测试数据，计算或辨识得到主泵模型参数，并计算主泵电磁转矩、水力转矩和摩擦转矩；
S3、基于步骤S1所建立的冷却剂主泵流量计算模型，仿真计算主泵惰转工况，电压、频率波动工况时冷却剂流量实时变化，与测试曲线作对比，校验模型及参数；
S4、将步骤S1所建立的冷却剂主泵流量计算模型接入核电机组系统模型，并接入电力系统仿真软件，实现核电机组整体动态仿真分析。
2.根据权利要求1所述的一种核电机组冷却剂屏蔽主泵动态建模及分析方法，其特征在于：在步骤S1中，所述冷却剂主泵流量计算模型的转矩平衡方程、质量流量方程，是采用以下方法建立的：
冷却剂屏蔽主泵转矩平衡方程为：
<mrow>
<msub>
<mi>I</mi>
<mi>p</mi>
</msub>
<mfrac>
<mrow>
<msub>
<mi>d&omega;</mi>
<mi>p</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>M</mi>
<mrow>
<mi>p</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>M</mi>
<mrow>
<mi>p</mi>
<mi>m</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，Ip为冷却剂主泵的转动惯量，ωp为冷却剂主泵转子角速度，Mpe为冷却剂主泵电磁转矩，Mpm为冷却剂主泵受到的总阻力矩；
主泵电磁转矩Mpe计算表达式为：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>p</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mi>p</mi>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;R</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mfrac>
<msubsup>
<mi>U</mi>
<mrow>
<mn>1</mn>
<mi>n</mi>
</mrow>
<mn>2</mn>
</msubsup>
<msub>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mi>n</mi>
</mrow>
</msub>
</mfrac>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msubsup>
<mi>sU</mi>
<mrow>
<mn>1</mn>
<mo>*</mo>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msubsup>
<mi>X</mi>
<mrow>
<mi>&sigma;</mi>
<mi>n</mi>
</mrow>
<mn>2</mn>
</msubsup>
<msubsup>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mo>*</mo>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<msubsup>
<mi>R</mi>
<mn>2</mn>
<mrow>
<mo>&prime;</mo>
<mn>2</mn>
</mrow>
</msubsup>
</mfrac>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<msub>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mo>*</mo>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msubsup>
<mi>sU</mi>
<mrow>
<mn>1</mn>
<mo>*</mo>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mrow>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<msubsup>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mo>*</mo>
</mrow>
<mn>2</mn>
</msubsup>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<msub>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mo>*</mo>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mi>p</mi>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;R</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mfrac>
<msubsup>
<mi>U</mi>
<mrow>
<mn>1</mn>
<mi>n</mi>
</mrow>
<mn>2</mn>
</msubsup>
<msub>
<mi>f</mi>
<mrow>
<mn>1</mn>
<mi>n</mi>
</mrow>
</msub>
</mfrac>
<mo>,</mo>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<msubsup>
<mi>X</mi>
<mrow>
<mi>&sigma;</mi>
<mi>n</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mrow>
<msup>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
式中，p为主泵电动机磁极对数，R2为主泵电动机转子每相绕组电阻，s为主泵电动机转差率，m1为主泵电动机磁极对数，U1n为电网侧系统基准电压，f1n为电网侧系统基准频率，上述各参数均为设计常数；Xσ为主泵电动机转子漏电抗，U1*、f1*分别为实际电网电压标幺值和频率标幺值；K1、K2为与冷却剂主泵相关的常系数；
屏蔽主泵因其飞轮和转子固封于承压部件内，使得其转子转动惯量小于常规泵；转速对转矩变化灵敏度高，摩擦转矩不可忽略，因此屏蔽主泵的阻力矩应同时计及冷却剂流量产生的阻力矩Mp1和由于机械摩擦造成的摩擦阻力矩Mp2，即有
Mpm＝Mp1+Mp2 (3)
冷却剂流量产生的阻力矩Mp1表达式为：
<mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>p</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>D</mi>
<mi>H</mi>
<mi>g</mi>
</mrow>
<mrow>
<msub>
<mi>&omega;</mi>
<mi>p</mi>
</msub>
<mi>&eta;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，D为冷却剂质量流量，ωp为冷却剂主泵转子角速度，η为冷却剂主泵效率，H为冷却剂主泵扬程；
Mp2表示主泵机械摩擦引起的阻力矩，采用经验公式表示为
<mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>p</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>K</mi>
<mn>3</mn>
</msub>
<msubsup>
<mi>&omega;</mi>
<mi>p</mi>
<mi>n</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，K3为常数，由泵的实验值确定；n取1.7-2.4；
冷却剂质量流量正比于主泵转子角速度，但由于系统惯性，冷却剂流量变化滞后于转子角速度变化，将冷却剂质量流量方程表示为
<mrow>
<msup>
<mi>D</mi>
<mo>*</mo>
</msup>
<mo>=</mo>
<msubsup>
<mi>&omega;</mi>
<mi>p</mi>
<mo>*</mo>
</msubsup>
<mo>+</mo>
<mi>&tau;</mi>
<mfrac>
<mrow>
<msubsup>
<mi>d&omega;</mi>
<mi>p</mi>
<mo>*</mo>
</msubsup>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，为转子角速度标幺值，D*为主泵质量流量标幺值，τ为常数，取0.2s。
3.根据权利要求2所述的一种核电机组冷却剂屏蔽主泵动态建模及分析方法，其特征在于：在步骤S2中，基于冷却剂主泵设计参数及测试数据，计算或辨识得到主泵模型参数的具体实现如下：
S21、基于冷却剂主泵惰转工况测试数据，采用MATLAB优化设计工具箱进行参数辨识计算，辨识出主泵摩擦阻力矩相关系数K3，n；
S22、基于电网侧频率或电压波动情况下，冷却剂流量测试数据，并将主泵摩擦阻力矩相关系数K3，n作为已知参数，采用MATLAB优化设计工具箱进行参数辨识，得到冷却剂主泵电磁转矩相关系数K1、K2；
S23、基于辨识的参数，建立相应的电磁转矩，摩擦转矩模型方程。
4.根据权利要求1所述的一种核电机组冷却剂屏蔽主泵动态建模及分析方法，其特征在于：所述步骤S3的具体实现过程如下：
S31、仿真冷却剂主泵惰转工况下，根据冷却剂流量变化曲线，求解惰转时间常数，与相关模型及参数进行比较分析，验证所述冷却剂主泵流量计算模型的有效性；
S32、仿真电网侧极端事故情况下，根据冷却剂流量与转速变化趋势，验证主泵低转速保护优先级高于冷却剂低流量保护，验证所述冷却剂主泵流量计算模型的合理性；
S33、仿真电网电压、频率阶跃波动工况下，根据冷却剂流量变化，验证冷却剂主泵流量计算模型对电网侧电压、频率扰动响应仿真的准确性。
5.根据权利要求1所述的一种核电机组冷却剂屏蔽主泵动态建模及分析方法，其特征在于：在步骤S4中，实现核电机组整体动态仿真分析具体为分析电网侧波动情况下，核电侧主泵流量，燃料温度，蒸汽压力相关参数动态响应过程。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710311080.6A CN107220414B (zh) | 2017-05-05 | 2017-05-05 | 一种核电机组冷却剂屏蔽主泵动态建模及分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710311080.6A CN107220414B (zh) | 2017-05-05 | 2017-05-05 | 一种核电机组冷却剂屏蔽主泵动态建模及分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107220414A true CN107220414A (zh) | 2017-09-29 |
| CN107220414B CN107220414B (zh) | 2021-01-29 |
Family
ID=59943778
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710311080.6A Active CN107220414B (zh) | 2017-05-05 | 2017-05-05 | 一种核电机组冷却剂屏蔽主泵动态建模及分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107220414B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109255176A (zh) * | 2018-08-31 | 2019-01-22 | 岭澳核电有限公司 | 动排气后一回路剩余空气体积标准值提升分析方法及系统 |
| CN109543223A (zh) * | 2018-10-19 | 2019-03-29 | 广西电网有限责任公司电力科学研究院 | 一种核电机组外特性分析方法 |
| CN110929407A (zh) * | 2019-11-29 | 2020-03-27 | 清华大学 | 一种供暖用屏蔽电泵的能效评价方法 |
| CN111247442A (zh) * | 2017-10-24 | 2020-06-05 | 三菱电机株式会社 | 异常诊断装置、异常诊断方法及异常诊断系统 |
| CN112016218A (zh) * | 2020-09-24 | 2020-12-01 | 哈尔滨电气动力装备有限公司 | 核主泵惰转模型的建模方法及应用 |
| CN117911195A (zh) * | 2024-03-18 | 2024-04-19 | 国网山西省电力公司电力科学研究院 | 一种异步电动机机群等效模型参数辨识方法及装置 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060120501A1 (en) * | 2004-08-31 | 2006-06-08 | Kabushiki Kaisha Toshiba | Power source for re-circulation pump and method of controlling the same |
| CN102279901A (zh) * | 2011-05-17 | 2011-12-14 | 湖北省电力公司电力试验研究院 | 一种针对第三代压水堆核电机组的建模方法 |
| CN103671228A (zh) * | 2013-12-10 | 2014-03-26 | 中广核工程有限公司 | 对反应堆冷却剂泵进行应急轴封的方法及轴封组件 |
| CN105545799A (zh) * | 2016-01-07 | 2016-05-04 | 江苏大学 | 一种核主泵的惰转模型叶轮水力设计方法 |
| CN105719092A (zh) * | 2016-01-26 | 2016-06-29 | 云南电网有限责任公司电力科学研究院 | 一种核电机组整体稳定性分析方法 |
| CN105760985A (zh) * | 2016-01-25 | 2016-07-13 | 武汉大学 | 一种针对核电接入电网的低频风险综合评估方法 |
-
2017
- 2017-05-05 CN CN201710311080.6A patent/CN107220414B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060120501A1 (en) * | 2004-08-31 | 2006-06-08 | Kabushiki Kaisha Toshiba | Power source for re-circulation pump and method of controlling the same |
| CN102279901A (zh) * | 2011-05-17 | 2011-12-14 | 湖北省电力公司电力试验研究院 | 一种针对第三代压水堆核电机组的建模方法 |
| CN103671228A (zh) * | 2013-12-10 | 2014-03-26 | 中广核工程有限公司 | 对反应堆冷却剂泵进行应急轴封的方法及轴封组件 |
| CN105545799A (zh) * | 2016-01-07 | 2016-05-04 | 江苏大学 | 一种核主泵的惰转模型叶轮水力设计方法 |
| CN105760985A (zh) * | 2016-01-25 | 2016-07-13 | 武汉大学 | 一种针对核电接入电网的低频风险综合评估方法 |
| CN105719092A (zh) * | 2016-01-26 | 2016-06-29 | 云南电网有限责任公司电力科学研究院 | 一种核电机组整体稳定性分析方法 |
Non-Patent Citations (3)
| Title |
|---|
| XU YI-MING等: ""Numerical Simulation of Interior Flow Field of Reactor Coolant Pump under Station Blackout Accident"", 《ADVANCED MATERIALS RESEARCH》 * |
| 徐一鸣等: ""核主泵惰转转速计算模型的比较"", 《发电设备》 * |
| 赵洁等: ""大型压水堆核电机组与电网相互影响机制的研究"", 《中国电机工程学报》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111247442A (zh) * | 2017-10-24 | 2020-06-05 | 三菱电机株式会社 | 异常诊断装置、异常诊断方法及异常诊断系统 |
| CN111247442B (zh) * | 2017-10-24 | 2022-04-26 | 三菱电机株式会社 | 异常诊断装置、异常诊断方法及异常诊断系统 |
| CN109255176A (zh) * | 2018-08-31 | 2019-01-22 | 岭澳核电有限公司 | 动排气后一回路剩余空气体积标准值提升分析方法及系统 |
| CN109255176B (zh) * | 2018-08-31 | 2023-05-09 | 岭澳核电有限公司 | 动排气后一回路剩余空气体积标准值提升分析方法及系统 |
| CN109543223A (zh) * | 2018-10-19 | 2019-03-29 | 广西电网有限责任公司电力科学研究院 | 一种核电机组外特性分析方法 |
| CN110929407A (zh) * | 2019-11-29 | 2020-03-27 | 清华大学 | 一种供暖用屏蔽电泵的能效评价方法 |
| CN110929407B (zh) * | 2019-11-29 | 2021-06-08 | 清华大学 | 一种供暖用屏蔽电泵的能效评价方法 |
| CN112016218A (zh) * | 2020-09-24 | 2020-12-01 | 哈尔滨电气动力装备有限公司 | 核主泵惰转模型的建模方法及应用 |
| CN112016218B (zh) * | 2020-09-24 | 2022-07-22 | 哈尔滨电气动力装备有限公司 | 核主泵惰转模型的建模方法及应用 |
| CN117911195A (zh) * | 2024-03-18 | 2024-04-19 | 国网山西省电力公司电力科学研究院 | 一种异步电动机机群等效模型参数辨识方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107220414B (zh) | 2021-01-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107220414A (zh) | 一种核电机组冷却剂屏蔽主泵动态建模及分析方法 | |
| CN109902433B (zh) | 压水堆非能动安全壳余热排出系统跨维度耦合方法 | |
| CN103246764B (zh) | 基于Ansys Workbench的潜油电机温度场模拟方法 | |
| CN102913464B (zh) | 离心泵转子瞬态流固耦合特性预测方法 | |
| Liu et al. | Pressure fluctuation prediction of a model Kaplan turbine by unsteady turbulent flow simulation | |
| CN1105950C (zh) | 汽轮机叶片温度预测方法 | |
| CN103336911B (zh) | 适用于多机电力系统调频分析的数学模型机理建模方法 | |
| CN103927416B (zh) | 一种余热排出泵抗震分析方法 | |
| CN104993482A (zh) | 电力系统暂态稳定性在线分析的分区合成方法及其装置 | |
| CN104978625A (zh) | 基于聚合rlc电路模型的电力系统次同步谐振分析方法 | |
| CN110516275A (zh) | 基于扰动信息的仿真模型参数校核方法 | |
| CN106778041A (zh) | 一种双馈电机三相短路电流最大值的简化计算方法 | |
| CN104008227A (zh) | 一种离心泵动静部件间隙优化设计方法 | |
| CN104484528A (zh) | 一种基于pscad的核电机组动力、电气混合仿真平台 | |
| An et al. | The correlation analysis of RCPs impeller geometrical parameters and optimization in coast-down process | |
| CN111130135A (zh) | 一种适用于高比例新能源接入下的电力系统惯量计算方法 | |
| Wang et al. | Study on the transient evolution law of internal flow field and dynamic stress of reactor coolant pump under rotor seizure accident | |
| CN104820727A (zh) | 发电机氢气冷却系统的仿真方法与装置 | |
| Lu et al. | Research on the influence of nonuniform inflow on steady-state hydraulic performance of reactor coolant pump | |
| CN105529736A (zh) | 撬棒投入双馈风电机群故障电流等值模型计算方法 | |
| Jayanthi et al. | Simulation and integrated testing of process models of PFBR operator training simulator | |
| Zhang et al. | Research of energy characteristics and fault position detection for cracked rotor system | |
| CN104809512B (zh) | 一种电力系统不同紧急控制措施性能指标的快速评估方法 | |
| CN104779613B (zh) | 基于试验的含变流器电力元件等效建模方法 | |
| CN103838961A (zh) | 超超临界汽轮机转子三维温度和热应力实时监测方法 |
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 |