发明名称 一种针对第三代压水堆核电机组的建模方法
摘要 本发明公开了一种针对第三代压水堆核电机组的建模方法,包括步骤:步骤1.将核电机组系统模型分解为若干子系统模型;步骤2.根据热工学和能量传递及转换规律,建立步骤1中所述子系统模型;步骤3.将步骤2所得的各子系统模型组合成核电机组全系统模型,并连接核电机组全系统模型与电力系统模型,得到核电机组与电力系统联合模型;步骤4.基于步骤3所得的核电机组与电力系统联合模型,建立第三代压水堆核电机组自定义模型,并根据上述自定义模型模拟核电机组性能及机网相互影响。本发明方法有效解决了核电机组与电力机组的联合仿真,可用于核电厂和电网的机网协调分析,实用性强。
申请公布号 CN102279901B 申请公布日期 2013.01.09
申请号 CN201110127164.7 申请日期 2011.05.17
申请人 湖北省电力公司电力科学研究院;武汉大学;国家电网公司 发明人 刘涤尘;赵洁;杜治;雷庆生;王静
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 武汉科皓知识产权代理事务所(特殊普通合伙) 42222 代理人 张火春
主权项 1.一种针对第三代压水堆核电机组的建模方法,其特征在于,包括以下步骤:步骤1,基于核电厂内各系统设备与电力系统的联系,将核电机组系统模型分解为若干子系统模型,所述子系统模型包括:堆芯中子动态模型、堆芯燃料及冷却剂温度模型、热线温度模型、冷线温度模型、一回路平均温度模型、蒸汽发生器模型、反应堆控制系统模型、汽轮发电机及其调速系统模型、汽轮机旁路调节系统模型,所述蒸汽发生器模型包括蒸汽发生器中的一回路冷却剂温度模型、U型管温度模型、二回路蒸汽压力模型,所述反应堆控制系统模型包括反应堆平均温度调节系统模型和反应堆功率控制棒调节系统模型;步骤2,根据热工学、能量传递及转换规律,建立步骤1中所述子系统模型,其中,a、采用能量平衡方程式建立堆芯燃料及冷却剂温度模型,所述堆芯燃料及冷却剂温度模型为:<maths num="0001"><![CDATA[<math><mrow><mfrac><msub><mi>dT</mi><mi>F</mi></msub><mi>dt</mi></mfrac><mo>=</mo><mfrac><mrow><mi>f</mi><msub><mi>P</mi><mn>0</mn></msub></mrow><mrow><msub><mi>m</mi><mi>F</mi></msub><msub><mi>c</mi><mi>PF</mi></msub></mrow></mfrac><mi>N</mi><mo>+</mo><mfrac><mi>hA</mi><mrow><msub><mi>m</mi><mi>F</mi></msub><msub><mi>c</mi><mi>PF</mi></msub></mrow></mfrac><mrow><mo>(</mo><msub><mi>T</mi><mi>av</mi></msub><mo>-</mo><msub><mi>T</mi><mi>F</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mfrac><msub><mi>dT</mi><mi>av</mi></msub><mi>dt</mi></mfrac><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>f</mi><mo>)</mo></mrow><msub><mi>P</mi><mn>0</mn></msub></mrow><mrow><msub><mi>m</mi><mi>C</mi></msub><msub><mi>c</mi><mi>PC</mi></msub></mrow></mfrac><mi>N</mi><mo>+</mo><mfrac><mi>hA</mi><mrow><msub><mi>m</mi><mi>C</mi></msub><msub><mi>c</mi><mi>PC</mi></msub></mrow></mfrac><mrow><mo>(</mo><msub><mi>T</mi><mi>F</mi></msub><mo>-</mo><mfrac><mrow><msub><mi>T</mi><mrow><mi>&theta;</mi><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>T</mi><mi>av</mi></msub></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mo>+</mo><mn>2</mn><mfrac><msub><mover><mi>m</mi><mo>&CenterDot;</mo></mover><mi>C</mi></msub><msub><mi>m</mi><mi>C</mi></msub></mfrac><mrow><mo>(</mo><msub><mi>T</mi><mrow><mi>&theta;</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>T</mi><mi>av</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><mfrac><msub><mi>dT</mi><mrow><mi>&theta;</mi><mn>2</mn></mrow></msub><mi>dt</mi></mfrac><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>f</mi><mo>)</mo></mrow><msub><mi>P</mi><mn>0</mn></msub></mrow><mrow><msub><mi>m</mi><mi>C</mi></msub><msub><mi>c</mi><mi>PC</mi></msub></mrow></mfrac><mi>N</mi><mo>+</mo><mfrac><mi>hA</mi><mrow><msub><mi>m</mi><mi>C</mi></msub><msub><mi>c</mi><mi>PC</mi></msub></mrow></mfrac><mrow><mo>(</mo><msub><mi>T</mi><mi>F</mi></msub><mo>-</mo><mfrac><mrow><msub><mi>T</mi><mrow><mi>&theta;</mi><mn>2</mn></mrow></msub><mo>+</mo><msub><mi>T</mi><mi>av</mi></msub></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mo>+</mo><mn>2</mn><mfrac><msub><mover><mi>m</mi><mo>&CenterDot;</mo></mover><mi>C</mi></msub><msub><mi>m</mi><mi>C</mi></msub></mfrac><mrow><mo>(</mo><msub><mi>T</mi><mi>av</mi></msub><mo>-</mo><msub><mi>T</mi><mrow><mi>&theta;</mi><mn>2</mn></mrow></msub><mo>)</mo></mrow></mrow></math>]]></maths>其中:T<sub>F</sub>为堆芯燃料温度;N为堆芯内中子通量密度;T<sub>θ1</sub>为反应堆冷却剂入口温度;T<sub>av</sub>为反应堆冷却剂平均温度;T<sub>θ2</sub>为反应堆冷却剂出口温度;P<sub>0</sub>为堆芯初始功率;f为燃料温升所占堆芯功率的百分比;h为堆芯中燃料与冷却剂之间的传热系数;A为堆芯中燃料与冷却剂之间的传热面积;m<sub>F</sub>、c<sub>PF</sub>为燃料的质量和比热;m<sub>C</sub>、c<sub>PC</sub>为堆芯中冷却剂的质量与比热;<img file="FDA00002139469300014.GIF" wi="50" he="48" />为冷却剂流过堆芯时的质量流量;b、采用热平衡方程式和一阶惯性环节建立热线温度模型、冷线温度模型和一回路平均温度模型,所述热线温度模型为:<maths num="0004"><![CDATA[<math><mrow><mfrac><msub><mi>dT</mi><mi>HL</mi></msub><mi>dt</mi></mfrac><mo>=</mo><mfrac><msub><mover><mi>m</mi><mo>&CenterDot;</mo></mover><mi>HL</mi></msub><msub><mi>m</mi><mi>HL</mi></msub></mfrac><mrow><mo>(</mo><msub><mi>T</mi><mrow><mi>&theta;</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mi>T</mi><mi>HL</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths>其中:m<sub>HL</sub>为热线中冷却剂的质量;<img file="FDA00002139469300016.GIF" wi="62" he="49" />为冷却剂流过热线时的质量流量;T<sub>θ2</sub>为反应堆冷却剂出口温度;T<sub>HL</sub>为热线温度;所述冷线温度模型为:<maths num="0005"><![CDATA[<math><mrow><mfrac><mrow><mi>d</mi><msub><mi>T</mi><mrow><mi>&theta;</mi><mn>1</mn></mrow></msub></mrow><mi>dt</mi></mfrac><mo>=</mo><mfrac><msub><mover><mi>m</mi><mo>&CenterDot;</mo></mover><mi>CL</mi></msub><msub><mi>m</mi><mi>CL</mi></msub></mfrac><mrow><mo>(</mo><msub><mi>T</mi><mi>CL</mi></msub><mo>-</mo><msub><mi>T</mi><mrow><mi>&theta;</mi><mn>1</mn></mrow></msub><mo>)</mo></mrow></mrow></math>]]></maths>其中:m<sub>CL</sub>为冷线中冷却剂的质量;<img file="FDA00002139469300021.GIF" wi="60" he="48" />为冷却剂流过冷线时的质量流量;T<sub>θ1</sub>为反应堆冷却剂入口温度;T<sub>CL</sub>为冷线温度;所述一回路平均温度模型为:<maths num="0006"><![CDATA[<math><mrow><msub><mi>T</mi><mi>avg</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><msub><mi>T</mi><mi>HL</mi></msub><mo>+</mo><msub><mi>T</mi><mi>CL</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>&tau;</mi><mi>c</mi></msub><mfrac><mrow><mi>d</mi><msub><mi>T</mi><mi>avg</mi></msub></mrow><mi>dt</mi></mfrac></mrow></math>]]></maths>其中:T<sub>avg</sub>为一回路平均温度测量值;T<sub>HL</sub>为热线温度;T<sub>CL</sub>为冷线温度;τ<sub>c</sub>为阻温探测器的时延;c、根据一回路和U型管之间的热阻参数、能量守恒定理和热传递规律,建立蒸汽发生器的一回路冷却剂温度模型;根据一回路到U型管之间的热阻参数和U型管与二回路之间的热阻参数,建立U型管温度模型;根据质量平衡、体积平衡及能量平衡方程式,建立蒸汽发生器的二回路蒸汽压力模型:所述一回路冷却剂温度模型为:<maths num="0007"><![CDATA[<math><mrow><msub><mi>M</mi><mi>sp</mi></msub><msub><mi>C</mi><mi>pc</mi></msub><mfrac><mrow><mi>d</mi><msub><mi>T</mi><mi>P</mi></msub></mrow><mi>dt</mi></mfrac><mo>=</mo><mn>2</mn><msub><mi>D</mi><mi>sp</mi></msub><msub><mi>C</mi><mi>pc</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>HL</mi></msub><mo>-</mo><msub><mi>T</mi><mi>P</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>K</mi><mi>sp</mi></msub><msub><mi>S</mi><mi>s</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>P</mi></msub><mo>-</mo><msub><mi>T</mi><mi>m</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><msub><mi>M</mi><mi>sp</mi></msub><msub><mi>C</mi><mi>pc</mi></msub><mfrac><mrow><mi>d</mi><msub><mi>T</mi><mi>CL</mi></msub></mrow><mi>dt</mi></mfrac><mo>=</mo><mn>2</mn><msub><mi>D</mi><mi>sp</mi></msub><msub><mi>C</mi><mi>pc</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>P</mi></msub><mo>-</mo><msub><mi>T</mi><mi>CL</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>K</mi><mi>sp</mi></msub><msub><mi>S</mi><mi>s</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>P</mi></msub><mo>-</mo><msub><mi>T</mi><mi>m</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths>其中,T<sub>P</sub>为蒸汽发生器中一回路冷却剂的平均温度;T<sub>m</sub>为U型管平均温度;T<sub>HL</sub>为蒸汽发生器中一回路冷却剂的入口温度;T<sub>CL</sub>为蒸汽发生器中一回路冷却剂的出口温度;M<sub>SP</sub>为蒸汽发生器中一回路冷却剂的质量;C<sub>pc</sub>为一回路冷却剂的比热;D<sub>sp</sub>为蒸汽发生器中一回路冷却剂的质量流量;K<sub>sp</sub>为一回路冷却剂至U型管的传热系数;S<sub>s</sub>为U型管的传热面积;所述U型管温度模型为:<maths num="0009"><![CDATA[<math><mrow><msub><mi>M</mi><mi>m</mi></msub><msub><mi>C</mi><mi>m</mi></msub><mfrac><msub><mi>dT</mi><mi>m</mi></msub><mi>dt</mi></mfrac><mo>=</mo><msub><mi>S</mi><mi>s</mi></msub><msub><mi>K</mi><mi>sp</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>p</mi></msub><mo>-</mo><msub><mi>T</mi><mi>m</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>K</mi><mi>ss</mi></msub><msub><mi>S</mi><mi>s</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mi>m</mi></msub><mo>-</mo><msub><mi>T</mi><mi>s</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths>其中:T<sub>m</sub>为U型管平均温度;T<sub>s</sub>为二回路侧的蒸汽饱和温度;M<sub>m</sub>、C<sub>m</sub>分别为U型管的质量和比热;K<sub>sp</sub>为一回路冷却剂至U型管的传热系数;S<sub>s</sub>为U型管的传热面积;T<sub>P</sub>为一回路冷却剂的平均温度;K<sub>ss</sub>为U型管至二回路侧冷却水的传热系数;所述蒸汽发生器的二回路蒸汽压力模型为:<maths num="0010"><![CDATA[<math><mrow><mfrac><msub><mi>dP</mi><mi>S</mi></msub><mi>dt</mi></mfrac><mo>=</mo><mfrac><mn>1</mn><mi>K</mi></mfrac><mo>[</mo><msub><mi>K</mi><mi>ss</mi></msub><msub><mi>S</mi><mi>s</mi></msub><msub><mi>T</mi><mi>m</mi></msub><mo>-</mo><msub><mi>K</mi><mi>ss</mi></msub><msub><mi>S</mi><mi>s</mi></msub><mfrac><msub><mrow><mo>&PartialD;</mo><mi>T</mi></mrow><mi>s</mi></msub><mrow><mo>&PartialD;</mo><msub><mi>P</mi><mi>s</mi></msub></mrow></mfrac><msub><mi>P</mi><mi>s</mi></msub><mo>+</mo><mi>&mu;</mi><mrow><mo>(</mo><msub><mi>h</mi><mi>fw</mi></msub><mo>-</mo><msub><mi>h</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow></math>]]></maths>其中:P<sub>S</sub>为蒸汽发生器出口蒸汽压力;K<sub>ss</sub>为U型管至二回路侧冷却水的传热系数;S<sub>s</sub>为U型管的传热面积;T<sub>m</sub>为U型管平均温度;μ为汽轮机汽门开度;h<sub>fw</sub>为给水焓;h<sub>s</sub>蒸汽发生器中二回路出口蒸汽的焓;<img file="FDA00002139469300031.GIF" wi="114" he="89" />为饱和蒸汽温度与压力之间的系数;<img file="FDA00002139469300032.GIF" wi="443" he="101" />M<sub>sf</sub>和h<sub>f</sub>分别为蒸汽发生器中二回路侧冷却水的质量和焓,M<sub>sg</sub>和h<sub>g</sub>分别为蒸汽发生器中二回路侧蒸汽的质量和焓;步骤3,将步骤2所得的各子系统模型组合成核电机组全系统模型,并通过汽轮发电机调速系统模型和发电机系统模型连接核电机组全系统模型与电力系统模型,得到核电机组与电力系统联合模型;步骤4,基于步骤3所得的核电机组与电力系统联合模型,建立第三代压水堆核电机组的自定义模型,并采用上述第三代压水堆核电机组自定义模型模拟核电机组性能及机网相互影响。
地址 430077 湖北省武汉市武昌区徐东大街361号