CN114117792A - 一种失水事故分析方法 - Google Patents
一种失水事故分析方法 Download PDFInfo
- Publication number
- CN114117792A CN114117792A CN202111424450.XA CN202111424450A CN114117792A CN 114117792 A CN114117792 A CN 114117792A CN 202111424450 A CN202111424450 A CN 202111424450A CN 114117792 A CN114117792 A CN 114117792A
- Authority
- CN
- China
- Prior art keywords
- loss
- model
- analysis
- coolant accident
- cladding
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 53
- 239000002826 coolant Substances 0.000 title claims abstract description 45
- 238000010206 sensitivity analysis Methods 0.000 claims abstract description 44
- 238000005253 cladding Methods 0.000 claims abstract description 42
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 24
- 238000000034 method Methods 0.000 claims description 55
- 238000009835 boiling Methods 0.000 claims description 32
- 238000004364 calculation method Methods 0.000 claims description 30
- 238000002347 injection Methods 0.000 claims description 29
- 239000007924 injection Substances 0.000 claims description 29
- 230000008569 process Effects 0.000 claims description 24
- 238000012546 transfer Methods 0.000 claims description 24
- 239000007788 liquid Substances 0.000 claims description 13
- 238000005070 sampling Methods 0.000 claims description 10
- 239000003381 stabilizer Substances 0.000 claims description 10
- 230000008961 swelling Effects 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 9
- 238000013461 design Methods 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 8
- 230000007704 transition Effects 0.000 claims description 8
- 230000001105 regulatory effect Effects 0.000 claims description 6
- ZOXJGFHDIHLPTG-UHFFFAOYSA-N Boron Chemical compound [B] ZOXJGFHDIHLPTG-UHFFFAOYSA-N 0.000 claims description 5
- 229910052796 boron Inorganic materials 0.000 claims description 5
- 238000001816 cooling Methods 0.000 claims description 5
- 238000009826 distribution Methods 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000002474 experimental method Methods 0.000 description 5
- 239000000446 fuel Substances 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 239000000243 solution Substances 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 3
- 239000008188 pellet Substances 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000007921 spray Substances 0.000 description 2
- 238000005507 spraying Methods 0.000 description 2
- 208000034320 Autosomal recessive spastic ataxia of Charlevoix-Saguenay Diseases 0.000 description 1
- 201000003679 Charlevoix-Saguenay spastic ataxia Diseases 0.000 description 1
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 238000007664 blowing Methods 0.000 description 1
- 239000000460 chlorine Substances 0.000 description 1
- 229910052801 chlorine Inorganic materials 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000001595 flow curve Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 125000006850 spacer group Chemical group 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Images
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
-
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Geometry (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明实施例提供一种失水事故分析方法,以实现在满足法规的基础上进行失水事故分析，包括：S1.对一回路系统和二回路系统进行水力学几何建模；S2.以包壳峰值温度为主要指标进行敏感性分析；S3.结合失水事故程序的特点使用关键模型进行分析；S4.建立满足双95％要求的分析方法。本发明实施例通过水力学几何建模、敏感性分析、关键模型分析和双95％要求分析，实现了在满足法规的基础上对失水事故的分析。
Description
技术领域
本发明涉及一种失水事故分析方法。
背景技术
国际上常用的失水事故分析方法大致分为两种：使用最佳估算程序，同时准确界定模型以及分析方法的不确定性；采用保守模型，同时也配合采用一定的分析方法来挖掘一部分的裕量。第一种方法由于需要对所有模型进行验证并且要界定其不确定性，因此耗时长且工作量巨大，因此在进行三代改进核电项目以及海外出口项目时，拟采用保守模型，同时也配合采用一定的分析方法来挖掘一部分的裕量的方法来进行失水事故分析。
美国联邦法规10CFR 50.46附录K模型(简称“附录K模型”)是一套在核电站设计中被广泛应用的LOCA事故应急堆芯冷却系统(ECCS)评价标准，是获得NRC认可的。目前，随着自主化三代核电以及核电出口的需求，需要研发出一套具有自主知识产权的分析程序和分析方法用于最复杂的设计基准事故(即失水事故)分析。
发明内容
本发明实施例提供一种失水事故分析方法,以实现在满足法规的基础上进行失水事故分析。
本发明实施例通过下述技术方案实现：
第一方面，本发明实施例提供一种失水事故分析方法，包括：
S1.对一回路系统和二回路系统进行水力学几何建模；
S2.以包壳峰值温度为主要指标进行敏感性分析；
S3.结合失水事故程序的特点使用关键模型进行分析；
S4.建立满足双95％要求的分析方法。
进一步的，所述步骤S1包括：
对安全系统进行建模；
对热构件进行模拟；
对于失水事故相关的保护信号进行模拟；
对系统关键参量进行稳态调试。
进一步的，所述步骤S1还包括：
增加控制变量来追踪事故过程中下封头硼浓度的变化用以标定旁通结束时间，标定堆芯入口平均流速的变化用以判断再淹没阶段是否需要启用蒸汽冷却惩罚模型；
启用和水力学计算相关的满足法规要求的模型；
启用和热构件相关的满足法规要求的模型。
进一步的，所述步骤S2包括：
对喷放系数进行敏感性分析，确定能够获得最高包壳峰值温度的喷放系数值；
对破口节块划分进行敏感性分析，确定能够获得最高包壳峰值温度的破口节块划分方法；
对高、中、低压安注注入点以及安注箱注入点进行节块划分敏感性分析，分别进行失水事故计算，确定能够获得最高包壳峰值温度的注入点节块划分方法；
对堆芯节块划分进行敏感性分析；
对再淹没过程二维移动精细热传导模型中的最大网格数进行敏感性分析；
对下降段、下封头、下腔室、上腔室、上封头节块划分依次进行敏感性分析；
对破损环路热管段、SG传热管节块划分进行敏感性分析，分别进行失水事故计算，确定能够获得最高包壳峰值温度的热管段、SG传热管节块划分方法；
对时间步长进行敏感性分析，确定能够获得最高包壳峰值温度的时间步长。
进一步的，所述步骤S3包括：
使用Moody模型计算两相临界流；
使用间隙热传导模型进行分析；
使用二维移动精细热传导模型来计算再淹没过程径向以及轴向过程的热传导；
使用包壳肿胀破裂模型模拟肿胀破裂之后造成的流道阻塞对失水事故后续过程的影响；
使用液位跟踪模型进行分析；
使用1971ANS模型来计算衰变热；
使用经过UPTF test6验证的CCFL模型进行分析；
使用最小膜态沸腾温度模型进行分析；
使用壁面与液体、壁面与蒸汽之间的固-液辐射换热模型进行分析。
进一步的，步骤S4包括：
选择一系列与电厂初始状态相关的参数作为抽样对象；
选择再淹没阶段与临界后换热相关的重要模型；
通过临界后试验数据界定膜态沸腾换热模型、过渡沸腾换热模型、最小膜态沸腾温度模型的不确定度范围；
进行随机抽样获得124组参数，生成计算矩阵；
开展124组大破口失水事故分析，获得所有计算工况的包壳峰值温度；选取第三高的值作为满足双95％要求的结果。
进一步的，对系统关键参量进行稳态调试包括：对堆芯压降、旁流份额、回路压降、稳压器初始压力、稳压器初始水位、反应堆平均温度、蒸汽发生器初始水位以及热工设计流量进行稳态调试。
进一步的，对热构件进行模拟包括：
在进行堆芯热构件建模时，通过FQ、Fdh和FAP参量确定各热构件的功率份额，并且分别构造热棒、热组件、平均组件的轴向功率分布形状。
进一步的，对于失水事故相关的保护信号进行模拟，包括：
对反应堆停堆信号、安注注入信号、汽轮机停机信号和主泵停运信号进行模拟。
进一步的，一回路系统包括堆芯、下降段、下腔室、下封头、上腔室、上封头、稳压器、冷管段、热管段、主泵和蒸汽发生器一次侧；二回路系统包括：蒸汽发生器、主蒸汽管道和汽轮机。
本发明实施例与现有技术相比，具有如下的优点和有益效果：
本发明实施例的一种失水事故分析方法，通过水力学几何建模、敏感性分析、关键模型分析和双95％要求分析，实现了在满足法规的基础上对失水事故的分析。
附图说明
为了更清楚地说明本发明示例性实施方式的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1为失水事故分析方法流程示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。
在以下描述中，为了提供对本发明的透彻理解阐述了大量特定细节。然而，对于本领域普通技术人员显而易见的是：不必采用这些特定细节来实行本发明。在其他实施例中，为了避免混淆本发明，未具体描述公知的结构、电路、材料或方法。
在整个说明书中，对“一个实施例”、“实施例”、“一个示例”或“示例”的提及意味着：结合该实施例或示例描述的特定特征、结构或特性被包含在本发明至少一个实施例中。因此，在整个说明书的各个地方出现的短语“一个实施例”、“实施例”、“一个示例”或“示例”不一定都指同一实施例或示例。此外，可以以任何适当的组合和、或子组合将特定的特征、结构或特性组合在一个或多个实施例或示例中。此外，本领域普通技术人员应当理解，在此提供的示图都是为了说明的目的，并且示图不一定是按比例绘制的。这里使用的术语“和/或”包括一个或多个相关列出的项目的任何和所有组合。
在本发明的描述中，术语“前”、“后”、“左”、“右”、“上”、“下”、“竖直”、“水平”、“高”、“低”“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明保护范围的限制。
实施例
第一方面，本发明实施例提供一种失水事故分析方法，包括：
S1.对一回路系统和二回路系统进行水力学几何建模；
S2.以包壳峰值温度为主要指标进行敏感性分析；
S3.结合失水事故程序的特点使用关键模型进行分析；
S4.建立满足双95％要求的分析方法。
进一步的，所述步骤S1包括：
对安全系统进行建模；
对热构件进行模拟；
对于失水事故相关的保护信号进行模拟；
对系统关键参量进行稳态调试。
进一步的，所述步骤S1还包括：
增加控制变量来追踪事故过程中下封头硼浓度的变化用以标定旁通结束时间，标定堆芯入口平均流速的变化用以判断再淹没阶段是否需要启用蒸汽冷却惩罚模型；
启用和水力学计算相关的满足法规要求的模型；
启用和热构件相关的满足法规要求的模型。
进一步的，所述步骤S2包括：
对喷放系数进行敏感性分析，确定能够获得最高包壳峰值温度的喷放系数值；
对破口节块划分进行敏感性分析，确定能够获得最高包壳峰值温度的破口节块划分方法；
对高、中、低压安注注入点以及安注箱注入点进行节块划分敏感性分析，分别进行失水事故计算，确定能够获得最高包壳峰值温度的注入点节块划分方法；
对堆芯节块划分进行敏感性分析；
对再淹没过程二维移动精细热传导模型中的最大网格数进行敏感性分析；
对下降段、下封头、下腔室、上腔室、上封头节块划分依次进行敏感性分析；
对破损环路热管段、SG传热管节块划分进行敏感性分析，分别进行失水事故计算，确定能够获得最高包壳峰值温度的热管段、SG传热管节块划分方法；
对时间步长进行敏感性分析，确定能够获得最高包壳峰值温度的时间步长。
进一步的，所述步骤S3包括：
使用Moody模型计算两相临界流；
使用间隙热传导模型进行分析；
使用二维移动精细热传导模型来计算再淹没过程径向以及轴向过程的热传导；
使用包壳肿胀破裂模型模拟肿胀破裂之后造成的流道阻塞对失水事故后续过程的影响；
使用液位跟踪模型进行分析；
使用1971ANS模型来计算衰变热；
使用经过UPTF test6验证的CCFL模型进行分析；
使用最小膜态沸腾温度模型进行分析；
使用壁面与液体、壁面与蒸汽之间的固-液辐射换热模型进行分析。
进一步的，步骤S4包括：
选择一系列与电厂初始状态相关的参数作为抽样对象；
选择再淹没阶段与临界后换热相关的重要模型；
通过临界后试验数据界定膜态沸腾换热模型、过渡沸腾换热模型、最小膜态沸腾温度模型的不确定度范围；
进行随机抽样获得124组参数，生成计算矩阵；
开展124组大破口失水事故分析，获得所有计算工况的包壳峰值温度；选取第三高的值作为满足双95％要求的结果。
进一步的，对系统关键参量进行稳态调试包括：对堆芯压降、旁流份额、回路压降、稳压器初始压力、稳压器初始水位、反应堆平均温度、蒸汽发生器初始水位以及热工设计流量进行稳态调试。
进一步的，对热构件进行模拟包括：
在进行堆芯热构件建模时，通过FQ、Fdh和FAP参量确定各热构件的功率份额，并且分别构造热棒、热组件、平均组件的轴向功率分布形状。
进一步的，对于失水事故相关的保护信号进行模拟，包括：
对反应堆停堆信号、安注注入信号、汽轮机停机信号和主泵停运信号进行模拟。
进一步的，一回路系统包括堆芯、下降段、下腔室、下封头、上腔室、上封头、稳压器、冷管段、热管段、主泵和蒸汽发生器一次侧；二回路系统包括：蒸汽发生器、主蒸汽管道和汽轮机。
具体失水事故分析方法如下：
(1)分析程序的选择
选择一个能够模拟失水事故的分析软件，例如：ARSAC、RELAP、CATHARE、WCOBRA/TRAC等。推荐使用ARSAC程序，因为该程序是中国核动力研究设计院自主开发的系统分析程序，和国外程序相比，在很多方面，例如：矩阵求解、再淹没模型等，都有一定的优势。
(2)满足法规要求的建模方法的确定
这里对建模过程中的关键流程和需要重点关注的内容进行说明：
①按自主化三代核电厂的设计特点，对一回路系统(包括堆芯、下降段、下腔室、下封头、上腔室、上封头、稳压器、冷管段、热管段、主泵、蒸汽发生器一次侧等)、二回路系统(包括：蒸汽发生器、主蒸汽管道、汽轮机等)进行水力学几何建模。确保节块划分合理、各部分连接关系正确。水力学建模中关键的问题在于将堆芯划分为1个热组件和若干个平均组件，并要模拟热组件和平均组件之间的横向流动(考虑截面突变的影响)。对于环形下降段的模拟需要根据回路数目来确定将其划分为几个部分(例如：对于三环路电厂，需将下降环腔分为三部分)，各个部分之间应该用横向接管连接。要输入正确的主泵特性曲线，并真实模拟燃料组件上大、小定位格架造成的局部阻力损失。为了能够满足“旁通阶段结束之前将注入至系统内部的安注水扣除”这一法规要求，建议在下降段上部区域模拟一个虚拟的接管和阀门，该接管和阀门执行的功能为在下封头出现一定量的含硼水的时刻之前强制性的将安注水抽出至系统之外。
②对专设安全系统进行建模，包括：安注箱、高(或中)压安注系统、低压安注系统等进行建模。对于安注箱的模拟要确保相对高度、初始总容积、气容积、水容积、初始压力以及其与冷管段的接管建模正确；对于高、中、低压安注，选取保守的压力-流量曲线。此外，对于破损环路，不模拟安注注入。
③对系统所涉及的热构件进行模拟，热构件包含两种类型：带热源的热构件、不带热源的热构件。带热源的热构件只存在堆芯内部，包括：热棒热构件、热组件热构件以及平均组件热构件。在进行堆芯热构件建模时，需要通过FQ(堆芯最高功率点功率与平均功率的比值)、Fdh(热组件焓升与平均组件焓升的比值)、FAP(热棒平均功率和热组件中除热棒之外的单根燃料棒平均功率的比值)等参量确定各热构件的功率份额，并且通过合理的方式分别构造热棒、热组件、平均组件的轴向功率分布形状。此外，需要对芯块及包壳的材料物性、气隙组成及物性、边界条件等进行合理模拟。对于不带热源的热构件，存在于很多关键部件中，包括：稳压器、蒸汽发生器U型管、上封头等。对于这些热构件需要合理模拟材料物性以及壁面厚度。
④对于失水事故相关的保护信号进行模拟，包括：反应堆停堆信号、安注注入信号、汽轮机停机信号、主泵停运信号等。为了保证分析结果的保守性，在考虑信号定值时需要添加一定的偏差，这些偏差将会导致更晚停堆、更晚启动安注。
⑤对系统关键参量进行稳态调试，包括：堆芯压降、旁流份额、回路压降、稳压器初始压力、稳压器初始水位、反应堆平均温度、蒸汽发生器初始水位、热工设计流量等。尽量保证以上参量与目标值偏差在1％范围以内。此外，为了满足法规的要求，在处理气隙组分以及芯块物性时，要满足获得储能保守这一条件。
⑥在以上建模的基础上需要进行更加细化的处理，其中包括：
a)增加一些控制变量来追踪事故过程中下封头硼浓度的变化用以标定旁通结束时间，标定堆芯入口平均流速的变化用以判断再淹没阶段是否需要启用蒸汽冷却惩罚模型。此外还需要标定在堆芯完全裸露之后再次有一定量流体流入的时刻(代表再淹没开始)等。
b)要结合程序的特征在建模时启用和水力学计算相关的满足法规要求的模型，包括：合理的临界流模型，水力学构件内部的液位跟踪模型，由于燃料棒肿胀破裂造成的流道阻塞模型等。
c)要结合程序的特征在建模时启用和热构件相关的满足法规要求的模型，包括：合理的包壳氧化模型，燃料肿胀破裂模型等。
(3)规范的敏感性分析流程的确定
为了保证分析结果的合理性，应该以获得较高包壳峰值温度为目标，按照如下顺序进行敏感性分析：
①对喷放系数进行敏感性分析，确定能够获得最高包壳峰值温度的喷放系数值；
②对破口节块划分进行敏感性分析，确定能够获得最高包壳峰值温度的破口节块划分方法；
③对高、中、低压安注注入点以及安注箱注入点进行节块划分敏感性分析，分别进行失水事故计算，确定能够获得最高包壳峰值温度的注入点节块划分方法。不同安注注入点敏感性分析应该分开进行；
④对堆芯节块划分进行敏感性分析。由于堆芯轴向布置有大、小定位格架，在保证每一个定位格架处在一个水力学节块的前提下，对其余部分进行加粗或细分。在获得多种堆芯节块划分的基础上分别进行失水事故计算，确定能够获得最高包壳峰值温度的堆芯节块划分方法；
⑤对再淹没过程二维移动精细热传导模型中的最大网格数进行敏感性分析，，确定能够获得最高包壳峰值温度的二维移动精细热传导模型中的最大网格数划分方法。
⑥对下降段、下封头、下腔室、上腔室、上封头节块划分依次进行敏感性分析，在进行敏感性分析时，要保证细分的每个节块体积一致。通过进行多组失水事故计算，确定能够获得最高包壳峰值温度的下降段、下封头、下腔室、上腔室、上封头节块划分方法；
⑦对破损环路热管段、SG传热管节块划分进行敏感性分析，分别进行失水事故计算，确定能够获得最高包壳峰值温度的热管段、SG传热管节块划分方法。
⑧对时间步长进行敏感性分析。确定能够获得最高包壳峰值温度的时间步长。由于大破口失水事故现象非常复杂，在计算初期如果选择较大的时间步长可能会引起较大的数值不稳定性导致计算无法继续开展，因此可以在破口发生后的2秒的时间内统一选择用较小的时间步长0.0005秒，对后续计算过程进行4组时间敏感性分析。
(4)关键模型使用方法的规定
为了满足附录K的要求，并且结合失水事故程序的特点，要按照如下要求启用相关模型或选择合适的模型(如果所选用的程序没有如下模型，则需要添加相关模型)：
a)两相临界流计算选择Moody模型；
b)要启用间隙热传导模型；
c)要使用二维移动精细热传导模型来计算再淹没过程径向以及轴向过程的热传导；
d)要使用合理保守的包壳肿胀破裂模型，并且要模拟肿胀破裂之后造成的流道阻塞对失水事故后续过程的影响；
e)在喷放阶段，对于临界热流密度的计算，高压区域推荐使用B&W-2关系式；低压区域推荐使用Barnett关系式(或在该关系式的基础上进行一定的修正)。在再淹没阶段，对于临界热流密度的计算推荐使用Zuber关系式(或在该关系式的基础上进行修正)。
f)在喷放阶段，要使用满足附录K要求的换热逻辑，一旦出现偏离泡核沸腾情形，则之后的换热模式只能是过渡沸腾、膜态沸腾、单相蒸汽对流换热中的一种；一旦进入过渡沸腾且壁面过热度超过300℉之后，则之后的换热模式只能是膜态沸腾、单相蒸汽对流换热中的一种。在再淹没阶段，若如果流速小于1.0inch/s，则临界后换热区域的换热模式需要强制修改为单相蒸汽对流换热模型。
g)对于喷放阶段临界后换热，推荐使用McDonough-Milich-King的关系式来计算过渡沸腾换热；推荐使用Groeneveld 5.7关系式来计算膜态沸腾换热。
h)为了较好的模拟CCFL现象，推荐使用经过UPTF test6验证的CCFL模型。
i)要使用1971ANS模型来计算衰变热。
j)在堆芯内部应启用液位跟踪模型。
k)应使用能够体现包壳材料表面特性以及氧化层厚度影响的最小膜态沸腾温度模型。
l)对于中、高流速(>1.0inch/s)下的再淹没过程，推荐使用FLECHT-SEASET再淹没实验、RBHT再淹没实验、NEPTUN再淹没实验(这三个实验是国际上公认的能够作为再淹没模型验证的实验)各两组中、高流速工况作为标定，判断程序是否会低估包壳峰值温度。如果程序没有低估这6个工况的包壳峰值温度，则不需要修正再淹没阶段临界后换热系数，如果程序低估了包壳峰值温度，则应该使用小于1.0的因子来修正再淹没阶段临界后换热系数以保证程序不会低估这6个工况的包壳峰值温度。
m)在计算再淹没阶段的膜态沸腾换热时，需要根据空泡份额的大小区分是反环状流膜态沸腾还是弥散流膜态沸腾。根据不同的模式选取相应的模型。对于反环状流膜态沸腾，推荐选取Bromley模型，对于弥散流膜态沸腾，推荐使用Forslund-Rohsenow关系式。
n)需要启用壁面与液体，壁面与蒸汽之间的固-液辐射换热模型。
(5)满足双95％要求的分析方法的建立
在对基准工况确定了建模方法、开展了敏感性分析工作以及确定了模型使用方法的基础上，为了挖掘一定的分析裕量，以利于设计优化(例如：提升功率，减小安注箱尺寸等)，将采用在一定的边界条件的范围内通过非参数顺序统计以及简单抽样的方法开展多组计算，选取计算最高值。
为了获得满足双95％的结果，首先需要确定计算组数。在我们所建立的方法中，拟计算的组数为59组(只关注包壳峰值温度这一参量)或124组(同时关注包壳峰值温度、包壳氧化层厚度、产氢氯这三个参量)。挑选计算结果中最高包壳峰值温度或第三高包壳峰值温度(如果开展59组计算，则选择最高包壳峰值温度；若开展124组计算，则选择第三高包壳峰值温度)作为最终的分析结果。
具体方法如下：
(1)选择一系列与电厂初始状态相关的参数作为抽样对象。这些参数包括：初始功率、稳压器初始压力、稳压器初始水位、FQ、Fdh、功率分布形状、初始平均温度、燃耗相关参数(包括：与燃耗相关的内压、芯块热导率、气隙种类与份额等)、破口类型和尺寸等。
(2)选择再淹没阶段与临界后换热相关的重要模型，包括：膜态沸腾换热模型、过渡沸腾换热模型、最小膜态沸腾温度模型，将这三个模型的不确定度范围作为抽样对象。
(3)根据自主化三代核电厂的设计特点和运行特点，界定与电厂初始状态相关的参数的取值范围。
(4)通过临界后试验数据界定膜态沸腾换热模型、过渡沸腾换热模型、最小膜态沸腾温度模型的不确定度范围。
(5)进行随机抽样获得124组参数(这些参数中包括：与电厂初始状态相关的参数、模型不确定性相关参数)，生成计算矩阵。
(6)开展124组大破口失水事故分析，获得所有计算工况的包壳峰值温度，并同时关注包壳氧化层厚度和产氢率。选取第三高的值作为满足双95％要求的结果。
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种失水事故分析方法，其特征在于，包括：
S1.对一回路系统和二回路系统进行水力学几何建模；
S2.以包壳峰值温度为主要指标进行敏感性分析；
S3.结合失水事故程序的特点使用关键模型进行分析；
S4.建立满足双95％要求的分析方法。
2.如权利要求1所述失水事故分析方法，其特征在于，所述步骤S1包括：
对安全系统进行建模；
对热构件进行模拟；
对于失水事故相关的保护信号进行模拟；
对系统关键参量进行稳态调试。
3.如权利要求2所述失水事故分析方法，其特征在于，所述步骤S1还包括：
增加控制变量来追踪事故过程中下封头硼浓度的变化用以标定旁通结束时间，标定堆芯入口平均流速的变化用以判断再淹没阶段是否需要启用蒸汽冷却惩罚模型；
启用和水力学计算相关的满足法规要求的模型；
启用和热构件相关的满足法规要求的模型。
4.如权利要求1所述失水事故分析方法，其特征在于，所述步骤S2包括：
对喷放系数进行敏感性分析，确定能够获得最高包壳峰值温度的喷放系数值；
对破口节块划分进行敏感性分析，确定能够获得最高包壳峰值温度的破口节块划分方法；
对高、中、低压安注注入点以及安注箱注入点进行节块划分敏感性分析，分别进行失水事故计算，确定能够获得最高包壳峰值温度的注入点节块划分方法；
对堆芯节块划分进行敏感性分析；
对再淹没过程二维移动精细热传导模型中的最大网格数进行敏感性分析；
对下降段、下封头、下腔室、上腔室、上封头节块划分依次进行敏感性分析；
对破损环路热管段、SG传热管节块划分进行敏感性分析，分别进行失水事故计算，确定能够获得最高包壳峰值温度的热管段、SG传热管节块划分方法；
对时间步长进行敏感性分析，确定能够获得最高包壳峰值温度的时间步长。
5.如权利要求1所述失水事故分析方法，其特征在于，所述步骤S3包括：
使用Moody模型计算两相临界流；
使用间隙热传导模型进行分析；
使用二维移动精细热传导模型来计算再淹没过程径向以及轴向过程的热传导；
使用包壳肿胀破裂模型模拟肿胀破裂之后造成的流道阻塞对失水事故后续过程的影响；
使用液位跟踪模型进行分析；
使用1971ANS模型来计算衰变热；
使用经过UPTF test6验证的CCFL模型进行分析；
使用最小膜态沸腾温度模型进行分析；
使用壁面与液体、壁面与蒸汽之间的固-液辐射换热模型进行分析。
6.如权利要求1所述失水事故分析方法，其特征在于，步骤S4包括：
选择一系列与电厂初始状态相关的参数作为抽样对象；
选择再淹没阶段与临界后换热相关的重要模型；
通过临界后试验数据界定膜态沸腾换热模型、过渡沸腾换热模型、最小膜态沸腾温度模型的不确定度范围；
进行随机抽样获得124组参数，生成计算矩阵；
开展124组大破口失水事故分析，获得所有计算工况的包壳峰值温度；选取第三高的值作为满足双95％要求的结果。
7.如权利要求2所述失水事故分析方法，其特征在于，对系统关键参量进行稳态调试包括：对堆芯压降、旁流份额、回路压降、稳压器初始压力、稳压器初始水位、反应堆平均温度、蒸汽发生器初始水位以及热工设计流量进行稳态调试。
8.如权利要求2所述失水事故分析方法，其特征在于，对热构件进行模拟包括：
在进行堆芯热构件建模时，通过FQ、Fdh和FAP参量确定各热构件的功率份额，并且分别构造热棒、热组件、平均组件的轴向功率分布形状。
9.如权利要求2所述失水事故分析方法，其特征在于，对于失水事故相关的保护信号进行模拟，包括：
对反应堆停堆信号、安注注入信号、汽轮机停机信号和主泵停运信号进行模拟。
10.如权利要求1所述失水事故分析方法，其特征在于，一回路系统包括堆芯、下降段、下腔室、下封头、上腔室、上封头、稳压器、冷管段、热管段、主泵和蒸汽发生器一次侧；二回路系统包括：蒸汽发生器、主蒸汽管道和汽轮机。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111424450.XA CN114117792B (zh) | 2021-11-26 | 2021-11-26 | 一种失水事故分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111424450.XA CN114117792B (zh) | 2021-11-26 | 2021-11-26 | 一种失水事故分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114117792A true CN114117792A (zh) | 2022-03-01 |
| CN114117792B CN114117792B (zh) | 2023-09-26 |
Family
ID=80370407
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111424450.XA Active CN114117792B (zh) | 2021-11-26 | 2021-11-26 | 一种失水事故分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114117792B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115186561A (zh) * | 2022-07-28 | 2022-10-14 | 中国核动力研究设计院 | 一种loca事故下燃料元件性能模拟方法、系统及存储介质 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2077561A1 (en) * | 2008-01-01 | 2009-07-08 | Institute of Nuclear Energy Research Atomic Energy Council, Executive Yuan | Thermal limit analysis with hot-channel model for boiling water reactors |
| CN107609313A (zh) * | 2017-10-18 | 2018-01-19 | 哈尔滨工程大学 | 一种船用非能动安全系统可靠性分析方法 |
| CN113536537A (zh) * | 2021-06-10 | 2021-10-22 | 中国核动力研究设计院 | 一种大破口失水事故分析方法及系统 |
-
2021
- 2021-11-26 CN CN202111424450.XA patent/CN114117792B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2077561A1 (en) * | 2008-01-01 | 2009-07-08 | Institute of Nuclear Energy Research Atomic Energy Council, Executive Yuan | Thermal limit analysis with hot-channel model for boiling water reactors |
| CN107609313A (zh) * | 2017-10-18 | 2018-01-19 | 哈尔滨工程大学 | 一种船用非能动安全系统可靠性分析方法 |
| CN113536537A (zh) * | 2021-06-10 | 2021-10-22 | 中国核动力研究设计院 | 一种大破口失水事故分析方法及系统 |
Non-Patent Citations (4)
| Title |
|---|
| 吴丹等: "基于ARSAC-K程序研究...模方法对失水事故分析的影响", 中国核学会2019年学术年会论文集第10册（核安全分卷、核安保分卷）, pages 148 - 154 * |
| 吴丹等: "基于ARSAC-K程序研究建模方法对失水事故分析的影响", 中国核学会2019年学术年会论文集第10册（核安全分卷、核安保分卷）, pages 148 - 154 * |
| 黄涛等: "适用于DMRM 评价模型的评估方法研究", 中国核科学技术进展报告（第七卷）——中国核学会2021年学术年会论文集第10册（核安全分卷、核安保分卷、核环保分卷）, pages 1 - 8 * |
| 黄涛等: "适用于DMRM评价模型的评估方法研究", 中国核科学技术进展报告（第七卷）--中国核学会2021年学术年会论文集第10册（核安全分册、核安保分卷、核环保分卷）, pages 1 - 8 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115186561A (zh) * | 2022-07-28 | 2022-10-14 | 中国核动力研究设计院 | 一种loca事故下燃料元件性能模拟方法、系统及存储介质 |
| CN115186561B (zh) * | 2022-07-28 | 2023-09-12 | 中国核动力研究设计院 | 一种loca事故下燃料元件性能模拟方法、系统及存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114117792B (zh) | 2023-09-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zheng et al. | Study on the DLOFC and PLOFC accidents of the 200 MWe pebble-bed modular high temperature gas-cooled reactor with TINTE and SPECTRA codes | |
| Gharari et al. | Study the effects of various parameters on hydrogen production in the WWER1000/V446 | |
| CN114117792B (zh) | 一种失水事故分析方法 | |
| Park et al. | Counterpart test for LSTF 1% cold-leg break LOCA (SB-CL-32) utilizing ATLAS test facility | |
| Zhou et al. | Research on analysis and modeling methods for LBLOCA in nuclear power plants based on autonomous LOCA analysis platform ARSAC | |
| Kotsarev et al. | Qualification of coupled 3D neutron kinetic/thermal hydraulic code systems by the calculation of a VVER-440 benchmark–Re-connection of an isolated loop | |
| Lizorkin et al. | Development and application of the coupled thermal-hydraulics and neutron-kinetics code ATHLET/BIPR-VVER for Safety Analysis | |
| Choi et al. | Plant-specific assessment of the natural circulation-induced creep rupture characteristics in the RCS pressure boundary during the SBO accident | |
| Safarzadeh et al. | Full scope simulation of VVER-1000 blowdown source and containment pressurization in a LBLOCA by parallel coupling of TRACE and CONTAIN | |
| Safarzadeh et al. | The large-break LOCA uncertainty analysis in a VVER-1000 reactor using TRACE and DAKOTA | |
| Bae et al. | Experimental investigation and validation of TASS/SMR-S code for single-phase and two-phase natural circulation tests with SMART-ITL facility | |
| Kral et al. | Thermal Hydraulic Analyses for PTS Evaluation: Comparison of Temperature Fields at RPV Predicted by System TH Code and CFD Code | |
| Porhemmat et al. | Primary loop analysis for a PWR contains Passive Core Cooling System; LOCA and clad rising temperature | |
| Uspuras et al. | Validation of coupled neutronic/thermal-hydraulic code RELAP5-3D for RBMK-1500 reactor analysis application | |
| Hadjam et al. | Simulation of a LBLOCA in the CALLISTO test facility using the best estimate computer code RELAP5/SCDAP3. 2 | |
| Li et al. | A Chinese 300MWe Two-Loop PWR NPP LBLOCA Analysis Based on the Deterministic Realistic Hybrid Methodology | |
| Feng et al. | Integral validation for two-fluid thermal hydraulic system analysis code LOCUST 1.2 based on LOFT L9-3 experiment | |
| Wang et al. | VALIDATION of TRANSIENT SYSTEM ANALYSIS CODE GINKGO BASED ON a SLB TEST | |
| Vihavainen | VVER-440 Thermal Hydraulics as a Computer Code Validation Challenge | |
| Langenbuch et al. | ATHLET/BIPR-VVER–an Advanced Coupled Code System for VVER Safety Analysis | |
| Grgic et al. | Coupled Code Analysis of the Rod Withdrawal at Power Accident including Uncertainty Evaluation Using CIAU-TN Method | |
| Kim et al. | Enhancement of safety analysis capability for a CANDU-6 reactor using RELAP-CANDU/SCAN coupled code system | |
| Grush et al. | Results and predictions of scaled, nuclear large break loss-of-coolant experiments | |
| Austregesilo et al. | Analysis of the TMI-2 accident with the code ATHLET-CD | |
| Zhang et al. | Research and Development on Simulator of Fast Reactor in China |
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 |