压水堆MSLB-SGTR事故模拟与分析:反应堆安全性验证

需积分: 10 0 下载量 164 浏览量 更新于2024-08-19 收藏 191KB PDF 举报
本文主要探讨了一种复杂且关键的核能系统事故——压水堆(Pressurized Water Reactor,PWR)中的主蒸汽管道断裂(Main Steam Line Break,MSLB)与蒸汽发生器传热管破裂(Steam Generator Tube Rupture,SGTR)叠加事故的热工水力学分析。该研究使用了高级仿真软件Relap 5/mod3.2进行事故过程模拟,以评估事故的演变情况和安全应对策略。 首先,研究人员对比了模拟计算得出的一回路压力、一次侧向二次侧泄漏流量等关键参数与实际参考电站的数据,结果显示两者变化趋势高度吻合,安全注射(Safety Injection,即在事故时向反应堆注入冷却水)等安全措施的响应也基本一致,这验证了所采用的模拟方法的准确性与可靠性。 接着,作者深入分析了事故过程中主要参数的变化机制,例如压力下降、泄漏流量增加的原因,这些变化通常涉及到蒸汽的快速泄漏、热量的迅速释放以及冷却系统的响应等因素。通过这种分析,作者能够构建出事故的完整序列,并提出操作员可能采取的干预措施,如控制蒸汽排放、启动应急冷却系统等。 文章的重点部分着重讨论了堆芯水位在事故干预下的变化。研究发现,在操作员在10分钟内有效干预的情况下,堆芯水位保持在安全范围内,避免了堆芯裸露,确保了反应堆的安全运行。这一结果对于理解和预防类似事故的发生具有重要意义。 最后,文章将研究成果归结到核反应堆事故分析的领域,强调了Relap 5/mod3.2程序在模拟这类复杂事故中的应用价值,以及堆芯保护措施对防止严重后果的关键作用。关键词包括反应堆事故分析、Relap 5软件、MSLB和SGTR破裂,这些都表明了研究的实用性和理论价值。 这篇论文提供了一种有效的工具和技术来理解和模拟压水堆中的多重故障情况,有助于核能行业的事故预防和应急响应策略的制定。

clc,clear,close all; b=dir('E:\Work\工作事务文件夹\潮汐课程设计\数据\B站水位数据.txt'); for i=1:length(b) fid=fopen(b(i).name); dataB{i}=textscan(fid,'%f-%f-%f %f:%f %f'); fclose(fid); end a=dir('E:\Work\工作事务文件夹\潮汐课程设计\数据\A站水位数据.txt'); for i=1:length(a) fid=fopen(a(i).name); dataA{i}=textscan(fid,'%f-%f-%f %f:%f %f'); fclose(fid); end %% %求日月距平 yearA=dataA{1,1}{1,1} yearB=dataB{1,1}{1,1} monthA=dataA{1,1}{1,2} monthB=dataB{1,1}{1,2} heightB=dataB{1,1}{1,6} heightA=dataA{1,1}{1,6} MSLB_L=mean(heightB); MSLA_L=mean(heightA); for i=1:1096%%三年的日距平 MSLA_S_day(i)=sum(heightA(i*24-23:i*24))/24;% A站日距平 MSLB_S_day(i)=sum(heightB(i*24-23:i*24))/24;% B站日距平 end MSLA_S_day=MSLA_S_day-MSLA_L; MSLB_S_day=MSLB_S_day-MSLB_L; T=1991; for i=1:3%%三年的月距平 m=find(yearA==T); T=T+1 hA=0; k=0; for j=1:12 k=monthA(m,1); hA=heightA(m,1); n=find(k==j); MSLA_S_month(1,(i-1)*12+j)=sum(hA(n,1))/length(n); end end T=1991; for i=1:3 m=find(yearB==T); T=T+1 hB=0; k=0; for j=1:12 k=monthB(m,1); hB=heightB(m,1); n=find(k==j); MSLB_S_month(1,(i-1)*12+j)=sum(hB(n,1))/length(n); end end MSLA_S_month=MSLA_S_month-MSLA_L; MSLB_S_month=MSLB_S_month-MSLB_L; figure(1) plot(MSLA_S_day) hold on plot(MSLB_S_day) legend('A站三年日距平','B站三年日距平') set(gca,'XTickLabel',{'1月','6月','12月','6月','12月','6月','12月'}); figure(2) plot(MSLA_S_month(1,:)) hold on plot(MSLB_S_month(1,:)) legend('A站三年月距平','B站三年月距平') set(gca,'XTickLabel',{'1月','6月','12月','6月','12月','6月','12月'}); %%保存数据 save MSLA_S_month ; save MSLB_S_month; save MSLA_S_day; save MSLB_S_day; save MSLA_L; save MSLB_L; 分析一下代码每一步

2023-07-10 上传