CN110970142B - 一种压水堆大破口失水事故始发应急工况预测方法 - Google Patents
一种压水堆大破口失水事故始发应急工况预测方法 Download PDFInfo
- Publication number
- CN110970142B CN110970142B CN201911147446.6A CN201911147446A CN110970142B CN 110970142 B CN110970142 B CN 110970142B CN 201911147446 A CN201911147446 A CN 201911147446A CN 110970142 B CN110970142 B CN 110970142B
- Authority
- CN
- China
- Prior art keywords
- temperature
- accident
- formula
- reactor
- coolant
- 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.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 75
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000000977 initiatory effect Effects 0.000 title claims description 6
- 238000004364 calculation method Methods 0.000 claims abstract description 34
- 239000000446 fuel Substances 0.000 claims abstract description 21
- 238000005253 cladding Methods 0.000 claims abstract description 18
- ATYZRBBOXUWECY-UHFFFAOYSA-N zirconium;hydrate Chemical compound O.[Zr] ATYZRBBOXUWECY-UHFFFAOYSA-N 0.000 claims abstract description 15
- 230000008859 change Effects 0.000 claims abstract description 13
- 238000001816 cooling Methods 0.000 claims abstract description 12
- 230000007774 longterm Effects 0.000 claims abstract description 12
- 238000005507 spraying Methods 0.000 claims abstract description 5
- 239000002826 coolant Substances 0.000 claims description 57
- 238000012546 transfer Methods 0.000 claims description 38
- 238000006243 chemical reaction Methods 0.000 claims description 16
- 238000002347 injection Methods 0.000 claims description 12
- 239000007924 injection Substances 0.000 claims description 12
- 230000004992 fission Effects 0.000 claims description 9
- QCWXUUIWCKQGHC-UHFFFAOYSA-N Zirconium Chemical compound [Zr] QCWXUUIWCKQGHC-UHFFFAOYSA-N 0.000 claims description 6
- 239000000376 reactant Substances 0.000 claims description 6
- 238000007664 blowing Methods 0.000 claims description 5
- 238000004519 manufacturing process Methods 0.000 claims description 5
- 229910052726 zirconium Inorganic materials 0.000 claims description 5
- 230000001052 transient effect Effects 0.000 abstract description 18
- 230000008569 process Effects 0.000 abstract description 8
- 238000004458 analytical method Methods 0.000 abstract description 7
- 238000011156 evaluation Methods 0.000 abstract description 6
- 239000011162 core material Substances 0.000 description 67
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 8
- 239000001257 hydrogen Substances 0.000 description 8
- 229910052739 hydrogen Inorganic materials 0.000 description 8
- 239000007788 liquid Substances 0.000 description 6
- 230000008018 melting Effects 0.000 description 6
- 238000002844 melting Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 238000013528 artificial neural network Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 239000003381 stabilizer Substances 0.000 description 4
- MCMNRKCIXSYSNV-UHFFFAOYSA-N Zirconium dioxide Chemical compound O=[Zr]=O MCMNRKCIXSYSNV-UHFFFAOYSA-N 0.000 description 3
- 229910001093 Zr alloy Inorganic materials 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 239000010865 sewage Substances 0.000 description 3
- 239000007921 spray Substances 0.000 description 3
- 229910045601 alloy Inorganic materials 0.000 description 2
- 239000000956 alloy Substances 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000004880 explosion Methods 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 239000012071 phase Substances 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 229910017076 Fe Zr Inorganic materials 0.000 description 1
- 229910052770 Uranium Inorganic materials 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000000306 component Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 230000003090 exacerbative effect Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- XLYOFNOQVPJJNP-ZSJDYOACSA-N heavy water Substances [2H]O[2H] XLYOFNOQVPJJNP-ZSJDYOACSA-N 0.000 description 1
- 229910001026 inconel Inorganic materials 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000007791 liquid phase Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 239000000155 melt Substances 0.000 description 1
- 229910044991 metal oxide Inorganic materials 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 239000008188 pellet Substances 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 125000006850 spacer group Chemical group 0.000 description 1
- 239000010935 stainless steel Substances 0.000 description 1
- 229910001220 stainless steel Inorganic materials 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000005514 two-phase flow Effects 0.000 description 1
- DNYWZCXLKNTFFI-UHFFFAOYSA-N uranium Chemical compound [U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U][U] DNYWZCXLKNTFFI-UHFFFAOYSA-N 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 230000008016 vaporization Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21C—NUCLEAR REACTORS
- G21C17/00—Monitoring; Testing ; Maintaining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
- G06Q50/265—Personal security, identity or safety
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- Development Economics (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Educational Administration (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Computer Security & Cryptography (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Plasma & Fusion (AREA)
- General Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明提供一种压水堆大破口失水事故始发应急工况预测方法，所述方法包括以下步骤：根据大破口失水事故特点将事故分为喷放、再灌水/再淹没和长期冷却三个阶段；计算包壳温度，包壳温度变化值＝(堆芯产热量+锆水反应放热量‑堆芯换热量)/燃料棒比热；根据堆芯温度评价事故应急工况。本发明提供的一种压水堆大破口失水事故始发应急工况预测方法，将复杂的一回路进行合理简化，对一回路进行瞬态分析，通过将事故划分为三个阶段，根据不同的阶段特点使用不同的经典公式近似估算应急工况进程，使得计算速度将大大加快，提升评价效率。
Description
技术领域
本发明属于核辐射安全技术领域，具体涉及一种压水堆大破口失水事故始发应急工况预测方法。
背景技术
福岛事故再一次为人类敲响了警钟，再一次表明：虽然反应堆严重事故发生概率极低，但是它仍然可能发生。IAEA提出，应直接根据堆芯损伤评价的结果指导防护行动，而无需依据源项估算结果。在应急情况下依据几个关键的反应堆工况参数实时、快速地预测未来的反应堆工况与事故进程，可以为应急决策提供支持，为场内和场外防护行动的实施赢得更多的时间。
由于神经网络的方法具有计算速度快这一显著优点，所以国内外学者使用该方法进行了一些事故工况预测的研究。
韩国学者Sangjun Park,Jinkyun Park等开发了二回路瞬态的监测系统，用以在事故早期阶段探查到瞬态的发生，确定瞬态的场景类型，并告知操作员距离汽轮机停机的剩余时间。在1000MW压水堆模拟机上的瞬态示例的数据库中，建立了所有的瞬态场景。瞬态示例展示了在没有操作员干预的情况下汽轮机停机之前核电站的行为。与此同时，从核电站的计算机上周期性地获取运行数据，并将之与瞬态示例数据库进行对比，找到与之高度相似的瞬态示例，从而确定瞬态的类型和预计的距离汽轮机停机剩余时间。对运行数据的预处理包括标准化、降低维数、确定权重因子。数据库包含上百个变量，所以很难迅速比较各个示例、及时给出结论，故使用基本组件分析方法降低变量的维度。本方法集中在预测汽轮机停机时间的计算和瞬态类型的确定上，然而本方法也可以用于在核电站中探查系统失效。主要流程为：建立二回路系统瞬态实例的数据库，与实时数据进行相似性分析，最终给出瞬态类型和距离汽轮机停机的剩余时间。详见图1。
韩国学者Seung Geun Kim等使用支持向量计算器，利用核电站状态变量输入在短时间内的变化，预测严重事故中可能发生的事件的发生时间。因为核电站的尺寸和复杂性一直在增加，操作员在理解核电站系统问题和缓解问题中面临着越来越大的挑战。在瞬态发生时，操作员要检查大量的仪表数据以确定异常事件的类型。由于数据量大、诊断时间短等复杂因素，操作员可能会做出错误的操作，从而加剧事故的严重性。如果在核电站发生瞬态事件，操作员会根据推荐操作规程估计异常事件的类型、尝试缓解异常事件的后果。对于操作员而言，严重事故相关事件的发生时间是非常重要的。因此，在严重事故发生时，能够提供此类信息的支持系统是大有作为的。首先，确定破口位置和破口尺寸。使用MAAP程序获得演练和测试数据。结果表明所开发的算法能够在LOCA事故中正确地确定破口位置，准确确定破口尺寸，也能够在合理的误差范围内预测主要事件的发生时间。根据这些结果，该文献提出的算法可以用于真实的核电站，因为该算法仅使用紧急停堆后短时间内的数据。如今有许多计算机代码可以进行事故相关的安全分析，但是它们的计算时间偏长，不适合作为实时支持系统。有研究使用人工智能技术处理监测数据，成功预测了事故相关事件的发生时间，包括堆芯裸露时间、堆芯温度超过1200℃的时间、压力容器失效时间。然而该研究只分析没有认为干预的情况，这种情况不符合实际。该研究将LOCA事故作为始发事件，分析在LOCA事故下的8种不同的事故序列，分析的堆型为APR1400。为了区分破口发生在热管段还是冷管段，支持向量分级(SVC)使用13种核电站状态参数：破口段蒸汽发生器(SG)的压力、水位、温度；未破口段SG的压力、水位、温度；稳压器的压力和水位；堆芯水位和水温；污水坑的水位；CTMT的压力和温度。确定破口尺寸也使用上述13种电站状态参数。当衰变热大于一回路向二回路的传热量时，堆芯最大温度才可能超过1200℃，所以选择与热损失相关的参数。压力容器失效可以看作是堆芯温度上升的延伸，所以选择的参数是一样的。预测堆芯最大温度超过1200℃使用的参数有：未破损段S/G的温度，安全壳压力、温度，压力容器水位；压力容器失效时间使用的参数：未破损段S/G的温度，安全壳压力、温度，压力容器水位；安全壳失效时间使用的参数(压力超过4个大气压时，认为安全壳失效)：稳压器压力，安全壳压力和温度，压力容器水位。得出的结果主要为不同破口尺寸与主要事件发生时间的对应关系曲线；未来可以继续开展的工作主要是降低预测的误差。
神经网络方法已经广泛地应用于核工程领域中，有许多学者使用神经网络的方法进行反应堆事故进程预测的研究。韩国学者Soon Ho Park等使用模糊神经网络预测严重事故中的压力容器水位，因为在严重事故中压力容器的水位无法测量到。该研究使用培训数据开发模型，并使用与之相互独立的测试数据进行验证；所用数据由MAAP4程序模拟OAP1000获得。
近来，核电站的严重事故已经成为全世界关心的问题。在事故初期，核电站的主要安全参数急剧变化，使得操作员没有充足的时间做出响应。事故期间，氢气从一回路泄漏至安全壳内，氢气浓度达到一定水平后会有发生氢爆的危险，严重威胁安全壳的完整性。氢气来源于锆水反应、安全壳内的蒸汽受到辐照活化和安全壳喷淋造成内部材料的腐蚀。为了防止发生氢气爆炸，安全壳内的氢气浓度应低于4％。韩国学者Dong Yeong Kim等使用人工智能的方法预测安全壳内氢气浓度的变化。该研究将破口位置分为热管段、冷管段以及蒸汽发生器传热管三种。破口尺寸分为210个步长；对于LOCA事故，破口尺寸从双端断裂尺寸的1/10000～1；对于蒸汽发生器，传热管破裂数量从1根变化到210根。破口尺寸的误差为0.4％。预测安全壳氢气浓度所用的两个输入为破口尺寸、停堆时长；使用MAAP4的模拟数据进行模型的开发和验证。
韩国学者Man Gyun Na,Sun Ho Shin等使用概率神经网络方法预测何时堆芯发生裸露、CET温度何时超过648.9℃(一般认为此时开始严重事故管理)、压力容器何时失效。热管段LOCA为始发事件时，预测堆芯何时裸露使用的输入有3个：安全壳温度、稳压器水位、完好段S/G温度。预测CET温度超过648.9℃所使用的输入有3个：污水池水位、压力容器水位、破口段S/G水位。模糊约束的数目为4个。预测压力容器失效的输入有3个：安全壳温度、稳压器水位、稳压器压力。冷管段LOCA为始发事件时，预测堆芯何时裸露使用的输入有3个：压力容器水位、污水池水位、堆芯水位。预测CET温度超过648.9℃所使用的输入有4个：压力容器水位、污水池水位、完好段S/G水位、破口段S/G水位。预测压力容器失效的输入有3个：安全壳压力、稳压器压力、破口段S/G水位。
清华大学核能与新能源技术研究院针对位于山东荣成的高温气冷堆示范电站，开发了一套核应急响应支持系统。该系统有两个任务：一是基于测量仪器的异常进行故障诊断；二是基于测量仪器的读数和操作员的行动预测事故进程。为了更好地执行这两项任务，开发了动态贝叶斯网络，在该网络中使用贝叶斯机器学习(LBP)算法进行事故诊断和预测。
综上所述，反应堆严重事故预测研究汇总，可总结为表1。
表1
针对压水堆大破口事故始发应急，现有技术中主要采用的方法是，以流体的质量守恒、和能量守恒方程为模型的基础，守恒方程的形式取决于所采用的两相流模型。但判断为气相或液相后，再采用质量守恒、和能量守恒方程计算工况评价所需参数的过程较为繁琐，需要的时间较长。
发明内容
针对现有技术中存在的缺陷，本发明的目的是提供一种压水堆大破口失水事故始发应急工况预测方法，将复杂的一回路进行合理简化，对一回路进行瞬态分析，通过将事故划分为三个阶段，根据不同的阶段特点使用不同的经典公式近似估算应急工况进程，使得计算速度将大大加快，提升评价效率。
为达到以上目的，本发明采用的技术方案是：一种压水堆大破口失水事故始发应急工况预测方法，所述应用包括以下步骤：根据大破口失水事故特点将事故分为喷放、再灌水/再淹没和长期冷却三个阶段；
计算包壳温度，包壳温度变化值＝(堆芯产热量+锆水反应放热量-堆芯换热量)/燃料棒比热；
根据包壳温度评价事故应急工况。
进一步的，堆芯换热量的计算公式为：
Qtransfer＝hAΔT
式中，h-换热系数，W/(m2·K)；A-燃料棒外表面积，m2；ΔT-燃料棒温度与冷却剂温度差，K。
进一步的，在喷放和再灌水/再淹没阶段，当堆芯冷却剂与燃料棒的对流方式属于强迫对流时，所述换热系数h的计算方法包括：
计算层流强迫对流换热系数包括以下公式：
Nu＝C(n)gdev
式中，C(n)-对于棒束和单管均取4.36；(z-z0)-距离流动入口的长度，1000m；Dh-水力学直径；Re-雷诺数；Pr-普朗特数；gdev-充分发展的Nusselt数；F(z)-无量纲入口长度；
计算紊流强迫对流换热系数公式为：
Nu＝0.023Re0.8Pr0.4
式中，Re-雷诺数；Pr-普朗特数；
比较层流强迫对流换热系数和紊流强迫对流换热系数的大小，取两者最大值作为强迫对流换热的对流换热系数h。
进一步的，在长期冷却阶段，堆芯冷却剂与燃料棒的对流方式属自然对流时，所述换热系数h的计算方法公式包括以下两种，取两者最大值作为自然对流换热的对流换热系数：
Nu＝0.18Ra1/4(L/Dh)-1/9
Nu＝0.065Ra1/3(L/Dh)-1/9
式中：L-流道长度；Ra-瑞利数；Dh-表征流动的特征长度。
进一步的，所述锆水反应放热量计算方式为：
Qox(T)＝Qox(T0)+Hrp(T)-Hrp(T0)
Hrp(T)＝Hr(T)-Hp(T)
式中：T0-参考温度，298.15K；Qox(T)-温度T下锆水反应放热量，J；Hr(T)-温度T下反应物焓值之和，J；HP(T)-温度T下反应物焓值之和，J；Qox(T0)-锆水反应对应值为5.797×106J/kgZr。
进一步的，所述堆芯产热量计算公式为：
Pdecay＝N(τ)+Ns1(τ)+Ns2(τ)
式中：N(τ)-τ时刻剩余中子引起的裂变功率，W；Ns1(τ)-τ时刻裂变产物的衰变功率，W；Ns2(τ)-τ时刻中子俘获产物的衰变功率，W；
N(τ)/N(0)＝0.15exp(-0.1τ)
式中：τ-停堆时间，s；N(0)-停堆时反应堆功率，W；
式中：Ns1(τ)-停堆τ秒后的裂变产物的衰变功率；A和a为系数；τ0-连续运行的时间；
式中：c和α为常数。
进一步的，当10-1s≤τ＜10s，A＝12.05，a＝0.0639；
当10s≤τ≤1.5×102s，A＝15.31，a＝0.1807；
当1.5×102s≤τ＜4×106s，A＝26.02，a＝0.2834；
当4×106s≤τ＜2×108s，A＝53.18，a＝0.3350。
进一步的，当为喷放阶段时，冷却剂温度的计算中需要的参数为堆芯压力、泄漏流量、堆芯水位；
当为再灌水/再淹没和长期冷却时，冷却剂温度的计算中需要的参数为中压安注、堆芯压力、堆芯水位。
本发明的效果在于，将复杂的一回路进行合理简化，对一回路进行瞬态分析，通过将事故划分为三个阶段，根据不同的阶段特点使用不同的经典公式近似估算应急工况进程，使得计算速度将大大加快，提升评价效率。
附图说明
图1为现有技术中诊断与预测流程示意图；
图2为压水堆大破口失水事故阶段划分原理示意图；
图3为堆芯活性区节点划分示意图；
图4为应急工况预测程序的计算流程示意图。
具体实施方式
为使本发明解决的技术问题、采用的技术方案和达到的技术效果更加清楚，下面将结合附图对本发明实施例的技术方案作进一步的详细描述。显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，均属于本发明保护的范围。
本发明提供一种压水堆大破口失水事故始发应急工况预测方法，所述方法包括以下步骤：
根据大破口失水事故特点将事故分为喷放、再灌水/再淹没和长期冷却三个阶段。参阅图2，各阶段的特点如下：
喷放阶段的特点包括：冷却剂快速泄露、一回路迅速泄压、堆芯泡份额增加以及传热恶化，堆芯温度上升；
再灌水/在淹没阶段的特点包括：中压安注开始注水、堆芯水位上升并完全淹没堆芯以及传热有所改善，堆芯温度暂时下降；
长期冷却阶段的特点包括：中压安注箱排空、衰变热加热堆芯剩余冷却剂、冷却剂汽化，从破口逸出、水位下降，堆芯逐渐裸露以及堆芯温度上升，直至熔化。
还需要说明的是，参阅图3，在本方法中将堆芯活性区沿径向划分为4个径向环，沿轴向划分为10个轴向层，共40个节点。在计算堆芯温度时，认为每个节点的温度是相同的。
还需要说明的是，参阅图4，根据堆芯压力是否大于安全壳压力来判断是否已由喷放阶段转换为再灌水/在淹没阶段；根据堆芯是否满水位判断是否已由再灌水/在淹没阶段转换为长期冷却阶段。
计算包壳温度，包壳温度变化值＝(堆芯产热量+锆水反应放热量-堆芯换热量)/燃料棒比热。
(1)堆芯换热量计算
堆芯换热量的计算公式为：
Qtransfer＝hAΔT
式中，h—换热系数，W/(m2·K)；A—燃料棒外表面积，m2；ΔT—燃料棒温度与冷却剂温度差，K。
其中，在喷放和再灌水/再淹没阶段，当堆芯冷却剂与燃料棒的对流方式属于强迫对流时，所述换热系数h的计算方法包括：
计算层流强迫对流换热系数以下公式：
Nu＝C(n)gdev
式中，C(n)-对于棒束和单管均取4.36；(z-z0)-距离流动入口的长度，1000m；Dh-水力学直径；Re-雷诺数；Pr-普朗特数；gdev-充分发展的Nusselt数；F(z)-无量纲入口长度；
计算紊流强迫对流换热系数公式为：
Nu＝0.023Re0.8Pr0.4
式中，Re-雷诺数；Pr-普朗特数；
比较层流强迫对流换热系数和紊流强迫对流换热系数的大小，取两者最大值作为强迫对流换热的对流换热系数h。
在长期冷却阶段，堆芯冷却剂与燃料棒的对流方式属自然对流时，所述换热系数h的计算方法包括以下两种，取两者最大值作为自然对流换热的对流换热系数：
Nu＝0.18Ra1/4(L/Dh)-1/9
Nu＝0.065Ra1/3(L/Dh)-1/9
式中：L-流道长度；Ra-瑞利数；Dh-表征流动的特征长度。
(2)锆水反应放热量计算
若在上述阶段中发生过锆水反应，则计算锆水反应放热量。锆水反应化学方程式为：
Zr+2H2O→ZrO2+2H2+Qox
式中，Qox表示反应放出的热量，J。
则，锆水反应放热量计算方式为：
Qox(T)＝Qox(T0)+Hrp(T)-Hrp(T0)
Hrp(T)＝Hr(T)-Hp(T)
式中：T0-参考温度，298.15K；Qox(T)-温度T下锆水反应放热量，J；Hr(T)-温度T下反应物焓值之和，J；HP(T)-温度T下反应物焓值之和，J；Qox(T0)-锆水反应对应值为5.797×106J/kgZr。
(3)堆芯产热量计算
堆芯产热量计算公式为：
Pdecay＝N(τ)+Ns1(τ)+Ns2(τ)
式中：N(τ)-τ时刻剩余中子引起的裂变功率，W；Ns1(τ)-τ时刻裂变产物的衰变功率，W；Ns2(τ)-τ时刻中子俘获产物的衰变功率，W；
N(τ)/N(0)＝0.15exp(-0.1τ)
式中：τ-停堆时间，s；N(0)-停堆时反应堆功率，W；
式中：Ns1(τ)-停堆τ秒后的裂变产物的衰变功率；A和a为系数；τ0-连续运行的时间。根据不同的堆停时间，A和a取不同的值，具体查下表2可得。
表2
中子俘获产物的衰变功率Ns2(τ)可由下式计算：
对于低富集度铀作燃料的压水堆，可取c＝0.6，a＝0.2。在其他实施例中，中子俘获产物的衰变功率Ns2(τ)在上式所得数值基础上再乘以系数1.1。
根据包壳温度评价事故应急工况。
根据包壳温度可判断堆芯工况，判断依据见表3。原始堆芯几何形状可在数分钟至数小时的时间内逐渐丧失，对于轻水堆，温度范围在1500K～3000K。对于低压事故序列，堆芯温度达到1000K～1200K时，锆包壳开始肿胀、破裂；对于高压事故序列，包壳在低温下坍塌到燃料上，锆合金包壳失效时间可能推迟，堆芯温度达到1500K时才会开始肿胀、破裂。
就严重事故而言，因堆芯材料氧化而放出的热量可能超出衰变热。金属锆氧化释热率在1500K时与衰变热相当，在1800K以上的温度下可能达到衰变热的10倍。
堆芯材料的整体熔化取决于每种堆芯材料的熔化温度，但由于堆芯内不同材料在相互紧密接触的情况下发生化学反应，其整体熔化会受到所形成低熔点合金的影响。由于因科镍合金、不锈钢和包壳的剧烈相互作用，定位格架将因低熔点合金的形成二在相对低的温度下而受损。对于典型压水堆设计，在1500K～1700K温度范围内，Fe-Zr、Ag-Zr等之间的化学反应会导致堆芯构件的早期液化与迁移。
锆合金的熔点一般介于2000K～2250K之间，包壳温度一旦超过2000K，锆合金包壳就会熔化，而且在某些情况下会流动到堆芯下部区域。熔融的锆会部分溶解UO2芯块和包壳外的ZrO2，使其在远低于自身熔点的温度下液化。
表3
在堆芯换热量的计算中运用到了ΔT(燃料棒温度与冷却剂温度差)。其中关于冷却剂温度的计算，根据不同的阶段特点，参阅图4计算所需的参数也不相同。
当为喷放阶段时，冷却剂温度的计算中需要的参数为堆芯压力、泄漏流量、堆芯水位；当为再灌水/再淹没和长期冷却时，冷却剂温度的计算中需要的参数为中压安注、堆芯压力、堆芯水位。
(1)中压安注
利用伯努利方程计算中压安注流量，当中压安注箱中的水耗尽或中压安注箱的压力小于等于一回路压力时，流量为0。
(2)堆芯压力
堆芯压力为冷却剂温度对应的饱和压力，使用下式计算：
式中：Ps-水蒸汽的饱和压力，kPa；T-水蒸汽温度，K；A、B、C为常数，查下表4可得。
表4
(3)泄漏流量计算
通过流入和流出一回路的水质量平衡来计算泄漏流量。流入流出一回路的各流量、一回路压力温度等，使用核电厂的实时测量数据。
泄漏流量使用下式计算：
式中，WB-泄漏流量，kg/s；Win-冷却剂流入一回路的流量，kg/s；Wout-冷却剂流出一回路的流量kg/s；
-计算时间间隔内，一回路冷却剂因温度压力变化而产生的质量变化，kg/s。
冷却剂流入一回路的流量使用下式计算：
Win＝WSI+Wfeed+Wseal+Wspray
式中，WSI-安注流量，kg/s；Wfeed-化容补水流量，kg/s；Wseal-主泵轴封注入流量，kg/s；Wspray-稳压器辅助喷淋流量，kg/s。
冷却剂流出一回路的流量使用下式计算：
Wout＝Wletdown+WSL
式中，Wletdown-一回路下泄流量，kg/s；WSL-主泵轴封回水流量，kg/s。
式中，
-一回路冷却剂除稳压器外的平均密度，kg/m3；V-一回路冷却剂系统除稳压器外的体积，m3；ρl_PRE ρl_PRE-稳压器液态冷却剂密度，kg/m3；S-稳压器横截面积，m2；hPRE-稳压器液位，m；ρg_PRE-稳压器气态冷却剂密度，kg/m3；HPRE-稳压器总高度，m。
(4)堆芯水位
初始堆芯冷却剂总质量减去冷却剂泄漏总质量即为堆芯冷却剂剩余质量(与泄漏量相比蒸发量很小可以忽略)，剩余冷却剂除以冷却剂密度得到剩余冷却剂体积，剩余冷却剂体积除以压力容器横截面积即可得到堆芯水位，即：
式中：Mtotal-初始堆芯冷却剂总质量，kg；Mleak-冷却剂泄漏总质量，kg；ρcoolant-冷却剂密度，kg/m3；SRPV-压力容器横截面积，m2。
(5)冷却剂温度计算
气态冷却剂与裸露的燃料元件换热，液态冷却剂与未裸露的燃料元件换热，温度上升均取决于换热量及其比热。由冷却剂温度初始值、每时间步长的温度变化，可知未来冷却剂随时间变化。
ΔTcoolant_g＝Qtran_g/cp_g
式中：ΔTcoolant_g-气态冷却剂温度变化；Qtran_g-气态冷却剂与堆芯换热量；cp_g-气态冷却剂比热。
ΔTcoolant_l＝Qtran_l/cp_l
式中：ΔTcoolant_l-液态冷却剂温度变化；Qtran_l-液态冷却剂与堆芯换热量；cp_l-液态冷却剂比热。
区别于现有技术，本发明提供的一种压水堆大破口失水事故始发应急工况预测方法，将复杂的一回路进行合理简化，对一回路进行瞬态分析，通过将事故划分为三个阶段，根据不同的阶段特点使用不同的经典公式近似估算应急工况进程，使得计算速度将大大加快，提升评价效率。
本领域技术人员应该明白，本发明所述的方法并不限于具体实施方式中所述的实施例，上面的具体描述只是为了解释本发明的目的，并非用于限制本发明。本领域技术人员根据本发明的技术方案得出其他的实施方式，同样属于本发明的技术创新范围，本发明的保护范围由权利要求及其等同物限定。
Claims (4)
1.一种压水堆大破口失水事故始发应急工况预测方法，其特征在于，所述方法包括以下步骤：
根据大破口失水事故特点将事故分为喷放、再灌水/再淹没和长期冷却三个阶段；
计算包壳温度，包壳温度变化值＝(堆芯产热量+锆水反应放热量-堆芯换热量)/燃料棒比热；
根据包壳温度评价事故应急工况；
堆芯换热量的计算公式为：
Qtransfer＝hAΔT
式中，h-换热系数，W/(m2·K)；A-燃料棒外表面积，m2；ΔT-燃料棒温度与冷却剂温度差，K；
在喷放和再灌水/再淹没阶段，当堆芯冷却剂与燃料棒的对流方式属于强迫对流时，所述换热系数h的计算方法包括：
计算层流强迫对流换热系数包括以下公式：
Nu＝C(n)gdev
式中，C(n)-对于棒束和单管均取4.36；(z-z0)-距离流动入口的长度，1000m；Dh-水力学直径；Re-雷诺数；Pr-普朗特数；gdev-充分发展的Nusselt数；F(z)-无量纲入口长度；
计算紊流强迫对流换热系数公式为：
Nu＝0.023Re0.8Pr0.4
式中，Re-雷诺数；Pr-普朗特数；
比较层流强迫对流换热系数和紊流强迫对流换热系数的大小，取两者最大值作为强迫对流换热的对流换热系数h；
在长期冷却阶段，堆芯冷却剂与燃料棒的对流方式属自然对流时，所述换热系数h的计算方法公式包括以下两种，取两者最大值作为自然对流换热的对流换热系数：
Nu＝0.18Ra1/4(L/Dh)-1/9
Nu＝0.065Ra1/3(L/Dh)-1/9
式中：L-流道长度；Ra-瑞利数；Dh-表征流动的特征长度；
当为喷放阶段时，冷却剂温度的计算中需要的参数为堆芯压力、泄漏流量、堆芯水位；
当为再灌水/再淹没和长期冷却时，冷却剂温度的计算中需要的参数为中压安注、堆芯压力、堆芯水位。
2.根据权利要求1所述的一种压水堆大破口失水事故始发应急工况预测方法，其特征在于，所述锆水反应放热量计算方式为：
Qox(T)＝Qox(T0)+Hrp(T)-Hrp(T0)
Hrp(T)＝Hr(T)-Hp(T)
式中：T0-参考温度，298.15K；Qox(T)-温度T下锆水反应放热量，J；Hr(T)-温度T下反应物焓值之和，J；HP(T)-温度T下反应物焓值之和，J；Qox(T0)-锆水反应对应值为5.797×106J/kgZr。
3.根据权利要求1所述的一种压水堆大破口失水事故始发应急工况预测方法，其特征在于，所述堆芯产热量计算公式为：
Pdecay＝N(τ)+Ns1(τ)+Ns2(τ)
式中：N(τ)-τ时刻剩余中子引起的裂变功率，W；Ns1(τ)-τ时刻裂变产物的衰变功率，W；Ns2(τ)-τ时刻中子俘获产物的衰变功率，W；
N(τ)/N(0)＝0.15exp(-0.1τ)
式中：τ-停堆时间，s；N(0)-停堆时反应堆功率，W；
式中：Ns1(τ)-停堆τ秒后的裂变产物的衰变功率；A和a为系数；τ0-连续运行的时间；
式中：c和α为常数。
4.根据权利要求3所述的一种压水堆大破口失水事故始发应急工况预测方法，其特征在于，
当10-1s≤τ＜10s，A＝12.05，a＝0.0639；
当10s≤τ≤1.5×102s，A＝15.31，a＝0.1807；
当1.5×102s≤τ＜4×106s，A＝26.02，a＝0.2834；
当4×106s≤τ＜2×108s，A＝53.18，a＝0.3350。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911147446.6A CN110970142B (zh) | 2019-11-21 | 2019-11-21 | 一种压水堆大破口失水事故始发应急工况预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911147446.6A CN110970142B (zh) | 2019-11-21 | 2019-11-21 | 一种压水堆大破口失水事故始发应急工况预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110970142A CN110970142A (zh) | 2020-04-07 |
| CN110970142B true CN110970142B (zh) | 2022-04-19 |
Family
ID=70031157
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911147446.6A Active CN110970142B (zh) | 2019-11-21 | 2019-11-21 | 一种压水堆大破口失水事故始发应急工况预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110970142B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111666655B (zh) * | 2020-05-08 | 2023-07-14 | 中国辐射防护研究院 | 一种六氟化铀泄漏源项的计算方法 |
| CN111540491B (zh) * | 2020-05-14 | 2022-04-01 | 中国核动力研究设计院 | 一种棒状燃料源项释放特性研究实验装置及其使用方法 |
| CN112418642B (zh) * | 2020-11-18 | 2022-03-11 | 中国核动力研究设计院 | 核电厂正常余热排出系统完全丧失事故的分析方法及系统 |
| CN112862148B (zh) * | 2020-12-21 | 2022-06-07 | 中国辐射防护研究院 | 后处理厂高放废液蒸发器红油爆炸事故释放源项估算方法 |
| CN113972016B (zh) * | 2021-10-26 | 2024-01-26 | 中国核动力研究设计院 | 核电厂安全壳外失水事故应对方法、装置、设备及介质 |
| CN115597010A (zh) * | 2022-10-12 | 2023-01-13 | 中广核工程有限公司（Cn） | 化容系统破口位置诊断方法、系统、装置以及存储介质 |
| CN116415414A (zh) * | 2023-02-16 | 2023-07-11 | 华中科技大学 | 压水堆失水事故数据适用性评价方法、系统、设备及终端 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107451398A (zh) * | 2017-07-07 | 2017-12-08 | 西安交通大学 | 压水堆核电厂严重事故分析方法 |
| CN110991006A (zh) * | 2019-11-06 | 2020-04-10 | 中国辐射防护研究院 | 一种基于裸露时间的压水堆大loca事故堆芯损伤评价方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140376678A1 (en) * | 2013-06-25 | 2014-12-25 | Robert H. Leyse | Method of and Apparatus for Monitoring a Nuclear Reactor Core Under Normal and Accident Conditions |
-
2019
- 2019-11-21 CN CN201911147446.6A patent/CN110970142B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107451398A (zh) * | 2017-07-07 | 2017-12-08 | 西安交通大学 | 压水堆核电厂严重事故分析方法 |
| CN110991006A (zh) * | 2019-11-06 | 2020-04-10 | 中国辐射防护研究院 | 一种基于裸露时间的压水堆大loca事故堆芯损伤评价方法 |
Non-Patent Citations (2)
| Title |
|---|
| Emergency condition analysis for MBLOCA along with SBO initiated severe accident using MELCOR;王宁等;《IOP Conference Series-Earth and Environmental Science》;20190626;全文 * |
| 大破口始发应急工况预测方法研究;王宁等;《辐射防护》;20201120;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110970142A (zh) | 2020-04-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110970142B (zh) | 一种压水堆大破口失水事故始发应急工况预测方法 | |
| US9997263B2 (en) | Method of determining the value of a parameter representative of the operability of a nuclear reactor, determining system, computer and corresponding medium | |
| Wang et al. | The development of a zirconium oxidation calculating program module for Module In-vessel Degraded Analysis Code MIDAC | |
| Saghafi et al. | Determination of PAR configuration for PWR containment design: A hydrogen mitigation strategy | |
| Wu et al. | System code evaluation of near-term accident tolerant claddings during boiling water reactor short-term and long-term station blackout accidents | |
| Hosseini et al. | Analysis of the natural circulation flow map uncertainties in an integral small modular reactor | |
| Liu et al. | Estimation of the depressurization process of Fukushima Daiichi NPP Unit 1 with SAMPSON | |
| Gauntt et al. | MELCOR simulations of the severe accident at the Fukushima Daiichi Unit 1 reactor | |
| Choi et al. | Plant-specific assessment of the natural circulation-induced creep rupture characteristics in the RCS pressure boundary during the SBO accident | |
| Ross et al. | Terry Turbopump Analytical Modeling Efforts in Fiscal Year 2016-Progress Report | |
| US20220051824A1 (en) | Method for protecting a nuclear reactor and corresponding nuclear reactor | |
| Saez et al. | Analysis of the turbine deblading in an HTGR with the CATHARE code | |
| Ahn | Creep Rupture Analysis of the RCS Pressure Boundary for a PWR SBO Accident | |
| Kim | Dynamic Event Tree Analysis of Accident Tolerant Fuel Safety Benefits Considering Steam Generator Tube Degradation | |
| Yu et al. | Accident and Safety Analysis | |
| Ma et al. | RESEARCH ON THE UTILITY OF MAAP AT THE ACCIDENT TRANSIENTS | |
| Pistora et al. | Evaluation of pressurised thermal shocks for VVER 440/213 reactor pressure vessel in NPP Dukovany | |
| Likhanskii et al. | Development of expert system for fuel monitoring and analysis in WWER-1000 units | |
| Cunning et al. | An Examination of CANDU Fuel Performance Margins Derived from a Statistical Assessment of Industrial Manufacturing Data | |
| Novog et al. | A statistical methodology for determination of safety systems actuation setpoints based on extreme value statistics | |
| Papini et al. | DEVELOPMENT OF A HOT-CHANNEL MODEL FOR THE VERIFICATION OF THE LOCA CRITERIA USING THE FRAPTRAN CODE | |
| Burgazzi et al. | Valutazioni di rischio e scenari incidentali-PAR 2017 | |
| Ning et al. | Method for Prognosis of LBLOCA Initiated Emergency Condition | |
| Xiao et al. | Analysis of hydrogen and source term under severe accident conditions for thousand megawatt PWR | |
| Скалозубов et al. | Modernized loss-of-coolant & blackout accident management strategy at nuclear power plants with wwer |
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 |