CN113536537B - 一种大破口失水事故分析方法及系统 - Google Patents
一种大破口失水事故分析方法及系统 Download PDFInfo
- Publication number
- CN113536537B CN113536537B CN202110647075.9A CN202110647075A CN113536537B CN 113536537 B CN113536537 B CN 113536537B CN 202110647075 A CN202110647075 A CN 202110647075A CN 113536537 B CN113536537 B CN 113536537B
- Authority
- CN
- China
- Prior art keywords
- steady
- transient
- calculation
- state calculation
- parameter
- 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 111
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000004364 calculation method Methods 0.000 claims abstract description 203
- 230000001052 transient effect Effects 0.000 claims abstract description 92
- 238000004458 analytical method Methods 0.000 claims abstract description 37
- 238000002347 injection Methods 0.000 claims abstract description 35
- 239000007924 injection Substances 0.000 claims abstract description 35
- 239000007788 liquid Substances 0.000 claims abstract description 21
- 230000005611 electricity Effects 0.000 claims abstract description 9
- 238000012795 verification Methods 0.000 claims abstract description 5
- 238000005070 sampling Methods 0.000 claims description 60
- 238000009826 distribution Methods 0.000 claims description 32
- ZOXJGFHDIHLPTG-UHFFFAOYSA-N Boron Chemical compound [B] ZOXJGFHDIHLPTG-UHFFFAOYSA-N 0.000 claims description 27
- 229910052796 boron Inorganic materials 0.000 claims description 27
- 238000012544 monitoring process Methods 0.000 claims description 13
- 230000004907 flux Effects 0.000 claims description 9
- 238000012952 Resampling Methods 0.000 claims description 7
- 238000009827 uniform distribution Methods 0.000 claims description 5
- 238000012986 modification Methods 0.000 claims description 4
- 230000004048 modification Effects 0.000 claims description 4
- 238000013524 data verification Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 3
- 230000014509 gene expression Effects 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 5
- 238000013210 evaluation model Methods 0.000 description 8
- 238000004590 computer program Methods 0.000 description 3
- 230000033228 biological regulation Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000005253 cladding Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 239000003381 stabilizer Substances 0.000 description 1
- 238000013076 uncertainty analysis Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
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
- 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
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)
- Testing And Monitoring For Control Systems (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明涉及反应堆热工水力设计及安全分析技术领域，具体公开了一种大破口失水事故分析方法及系统。选取核电站大破口失水事故相关的指标参数；建立稳态计算模型，对大破口失水事故相关参数进行稳态计算，并进行稳态计算后的参数值校验；构建瞬态计算模型，并对大破口失水事故相关参数进行瞬态计算，并根据计算结果与实际的破口、核电厂外电情况进行对比，并在出现偏离时，重新构建瞬态计算模型并进行瞬态计算；进行安注水扣除，获得液位参数值，并对瞬态模型进行更新后，在稳态计算结果基础上进行再计算，并在完成计算后，进行参数显示及分析。该方法和系统解决了压水堆大破口失水事故分析的工况多、流程繁琐、人因失误率高的难题。
Description
技术领域
本发明属于反应堆热工水力设计及安全分析技术领域，具体涉及一种大破口失水事故分析方法及系统。
背景技术
大破口事故决定核电站运行功率的最要设计基准事故之一，其分析评价模型可以从保守评价模型、组合评价模型、最佳估算模型和风险指引评价模型中选择。这些模型主要包括计算机程序、系统部件可用性假设、初始条件和边界条件四部分。当计算机程序采用满足法规要求的附录K模型的系统程序、系统部件可用性采用保守假设、初始条件和边界条件采用实际数据+不确定性时，形成失水事故组合评价模型DMRM，该组合评价模型具有以下优点：1)保守的计算机程序开发成本低，保守性满足法规要求；2)实际数据+不确定性的初始条件和边界条件，释放了过度的保守裕量，在相同的系统设计下，能够提升系统容许的功率上限。
然而采用DMRM组合评价模型，不确定性分析需要开展多工况分析，每个工况的计算流程包括参数抽样、稳态输入卡生成、稳态工况计算、瞬态输入卡修改、瞬态计算、安注水扣除、再启动计算和结果处理等多个过程，依靠技术人员开展完整失水事故分析，不仅耗费分析时间，拖延设计进度，而且容易造成人因失误，无法确保技术质量。
发明内容
本发明的目的在于提供一种大破口失水事故分析方法及系统，解决大坡口失水事故分析效率低、分析过程人因失误的问题。
本发明的技术方案如下：一种大破口失水事故分析方法，该方法包括：
选取核电站大破口失水事故相关的指标参数，并对各参数进行样本抽样后进行参数分布检查，以验证参数样本值的有效性；
建立稳态计算模型，对大破口失水事故相关参数进行稳态计算，根据稳态计算值与参数抽样值的偏差值是否超过设定阈值，来进行稳态计算后的参数值校验；
构建瞬态计算模型，并对大破口失水事故相关参数进行瞬态计算，并根据计算结果与实际的破口、核电厂外电情况进行对比，并在出现偏离时，重新构建瞬态计算模型并进行瞬态计算；
进行安注水扣除，获得液位参数值，并对瞬态模型进行更新后，在稳态计算结果基础上进行再计算，并在完成计算后，提取大破口失水事故相关参数，并进行参数显示及分析。
所述核电站大破口失水事故相关的指标参数包括核电站大破口失水事故相关的保守性假设、初始条件和边界条件的参数、参数范围和分布形式构成的指标参数集；
所述对抽样的参数样本值有效性进行验证是通过对抽样后的参数分布进行检查从而确定抽样参数样本的有效性，其中，热流密度热通道因子FQ和核焓升热通道因子FΔH两参数抽样值需要满足正态分布；反应堆功率、反应堆平均温度、系统压力、安注箱水体积、安注箱压力、安注箱水温、安注水温度、堆芯上部功率分布、堆芯中部功率分布、破口面积、厂外电是否丧失，上述参数需要满足均匀分布。
所述构建稳态计算模型是通过获得的各工况下的参数抽样样本，根据各个参数的物理含义，利用python的正则表达式对关键词进行定位，将每个抽样工况获得的参数值写入输入文本卡对应的位置，完成稳态计算模型构建。
所述对稳态计算输出文件校验是通过对稳态计算输出文件进行识别，并判断稳态计算结束时间与预设时间是否相同，从而判断稳态计算是否完成而进行校验，若未完成稳态计算，则返回修改抽样形成的文本卡信息，并重新进行稳态计算。
所述瞬态计算模型是在完成稳态计算后，通过将稳态计算获得的二回路蒸汽流量参数值及抽样获得的破口面积、厂外电丧失情况信息写入瞬态卡中，形成瞬态计算模型。
所述液位参数值获取是在完成瞬态计算后，提取计算结果中的硼浓度值，并不断检测核电厂模型中下腔室内的硼浓度，当硼浓度超过所设定的阈值时，获得此时刻的水容积差值，从而计算此时安注箱内剩余水容积，并根据球形安注箱容积及水位的关系，将水容积转化获得液位值。
所述瞬态计算或稳态计算是利用python的os模块调用包含附录K模型的系统程序，进行瞬态计算或稳态计算。
所述更新瞬态模型后，在稳态计算结果基础上进行再计算时，是将安注箱液位值更新在已建立模型的瞬态卡信息中，并在稳态计算结果的基础上进行再计算；所述对大破口失水事故相关参数监控，并进行显示和分析是通过提取大破口失水事故相关参数值后，利用matplotlib库完成结果图的自动绘制，基于文档模板，自动生成分析报告。
一种大破口失水事故分析系统，该系统参数抽样模块、稳态计算模块、瞬态计算模块以及硼浓度监测模块，其中，所述参数抽样模块与所述稳态计算模块相互连接，并通过所述参数抽样模块可将采集的核电站大破口失水事故相关参数值形成输入文本后，输入所述稳态计算模块进行稳态计算；所述稳态计算模块将稳态计算值传输至相互连接的所述瞬态计算模块进行瞬态计算；所述硼浓度监测模块与所述瞬态计算模块相互连接，并通过所述硼浓度监测模块从瞬态计算模型提取并获得安注箱液位值后，对所述瞬态计算模型中的瞬态计算进行更新再计算，并利用与所述瞬态计算模块相连接的数据显示模块对相关参数进行显示和分析。
所述参数抽样模块可采集对失水事故下PCT影响较大的参数，包括但不限于反应堆功率、反应堆平均温度、系统压力、安注箱水体积、安注箱压力、安注箱水温、安注水温度、热流密度热通道因子FQ、核焓升热通道因子FΔH、堆芯上部功率分布、堆芯中部功率分布、破口面积、厂外电是否丧失；所述参数抽样模块可对采集的抽样参数进行数据验证，其中，采集的参数除热流密度热通道因子FQ、核焓升热通道因子FΔH两参数为正态分布外，其它参数均需满足均匀分布，如某项参数不满足该项参数的分布要求，则需要进行重新抽样。
所述参数抽样模块可采集核电厂工况下的各参数样本值，并形成输入文本卡；所述稳态计算模块根据所述输入文本卡，将每个抽样工况获得的参数值写入文本卡对应的位置，形成稳态计算模型；对所述文本卡中的参数进行稳态计算，对稳态计算后的参数进行核实，并在获得参数稳态计算值与抽样获得的参数值的偏差大于设定的阈值时，向参数抽样模块发送相应参数重新抽样请求信号。
所述稳态计算模块对稳态计算输出文件进行识别，并根据稳态计算结束时间与预设时间的一致性来判断稳态计算是否完成，并在未完成稳态计算时，则向所述参数抽样模块发生修改输入文本卡信息的请求信息，并利用所述稳态计算模块进行稳态计算。
所述瞬态计算模块将稳态计算模块稳态计算后的二回路蒸汽流量参数值及稳态计算时的破口面积和厂外电丧失情况信息输入瞬态卡，并进行瞬态计算，并根据瞬态计算结果与实际抽样获得的破口、厂外点情况信息进行校验，若出现偏离，瞬态计算模块则向稳态计算模块发送请求，在进行参数修改后，继续进行瞬态计算。
所述硼浓度监测模块可提取所述瞬态计算模块计算后的硼浓度值，并当硼浓度值超过所设定的阈值时，获得此刻水容积差值，并计算相应时刻安注箱内剩余水容积值，并根据球形安注箱容积及水位的关系，获得液位值；将该液位值传输至瞬态计算模块，并更新瞬态计算模块中的瞬态信息后进行再计算。
本发明的显著效果在于：本发明所述的一种大破口失水事故分析方法和系统，解决了压水堆大破口失水事故分析的工况多、流程繁琐、人因失误率高的难题。利用该发明，能够显著提升压水堆核电厂大破口失水事故分析的效率，缩短大破口失水事故评价周期，降低人因失误率，确保评价结果的合理性。该方法可用于我国自主化三代核电厂华龙一号及后续电厂的大破口失水事故分析，同时也可用于小破口失水事故分析以及一些需要开展大量敏感性研究的事故分析，有助于提升安全分析的效率，确保分析过程中参数处理的传递的正确性。
附图说明
图1为本发明所述的一种大破口失水事故分析系统结构示意图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细说明。
一种大破口失水事故分析方法，该方法具体包括如下步骤：
S1、选取核电站大破口失水事故相关的指标参数，并对各参数进行样本抽样后进行参数分布检查，以验证参数样本值的有效性；
分析核电站大破口失水事故相关的保守性假设、初始条件和边界条件的参数、参数范围和分布形式，并将上述信息以JSON格式输入到制定的输入文本中，其中，所制定的输入文本中抽样参数包含对失水事故下PCT影响较大的参数，包括但不限于反应堆功率、反应堆平均温度、系统压力、安注箱水体积、安注箱压力、安注箱水温、安注水温度、热流密度热通道因子FQ、核焓升热通道因子FΔH、堆芯上部功率分布、堆芯中部功率分布、破口面积、厂外电是否丧失，上述参数分布除热流密度热通道因子FQ、核焓升热通道因子FΔH两参数为正态分布外，其它参数均需满足均匀分布；
对上述参数进行抽样并获得各参数的样本值，并按工况进行输出；对每个工况下的参数分布进行检查，若任一参数不满足分布要求，则需要进行重新抽样；
S2、建立稳态计算模型，对大破口失水事故相关参数进行稳态计算，根据稳态计算值与参数抽样值的偏差值是否超过设定阈值，来进行稳态计算后的参数值校验；
S2.1、构建稳态计算模型；
针对获得的各工况下的参数抽样样本，根据各个参数的物理含义，利用python的正则表达式对关键词进行定位，将每个抽样工况获得的参数值写入输入文本卡对应的位置，完成稳态计算模型构建；
S2.2、进行各个工况下的稳态计算，并在计算获得的某项参数稳态计算值与抽样获得的参数值的偏差大于所设定的阈值后，转入步骤S1进行重新抽样；
利用python的os模块调用包含附录K模型的系统程序，对输入文本卡中的参数进行稳态计算，完成稳态计算后，主要对上述参数进行核实，若计算获得的参数稳态计算值与抽样获得的参数值的偏差大于设定的阈值，则返回步骤S1；
S2.3、对稳态计算输出文件进行校验，实时判定稳态计算是否完成；
对稳态计算输出文件进行识别，并判断稳态计算结束时间是否与预设时间相等，从而判断稳态计算是否成功或完成，若未完成稳态计算，则返回修改输入文本卡信息，并重新进行稳态计算；
S3、构建瞬态计算模型，并对大破口失水事故相关参数进行瞬态计算，并根据计算结果与实际的破口、核电厂外电情况进行对比，并在出现偏离时，重新构建瞬态计算模型并进行瞬态计算；
S3.1、在完成稳态计算后，将稳态计算获得的参数写入瞬态卡中，建立瞬态计算模型；
在判断稳态计算完成后，将稳态计算获得的二回路蒸汽流量参数值写入瞬态模型中，并将抽样获得的破口面积和核电厂外电丧失情况写入瞬态卡中，形成瞬态计算模型；
S3.2、进行瞬态计算，并对瞬态计算结果进行校验；
利用python的os模块调用包含附录K模型的系统程序，进行瞬态计算，并将瞬态计算结果与实际抽样的破口、核电厂外电情况进行校验，若出现偏离，则继续进行瞬态模型构建及瞬态计算；
S4、进行安注水扣除，获得液位参数值，并对瞬态模型进行更新后，在稳态计算结果基础上进行再计算，并在完成计算后，提取大破口失水事故相关参数，并进行参数显示及分析；
S4.1、完成瞬态计算后，利用所获得的硼浓度值，获得安注箱液位值；
在完成瞬态计算后，提取计算结果中的硼浓度值，并不断检测核电厂模型中下腔室内的硼浓度，当硼浓度超过所设定的阈值时，获得此时刻的水容积差值，从而计算此时安注箱内剩余水容积，并根据球形安注箱容积及水位的关系，将水容积转化获得液位值；
S4.2、更新瞬态模型，并在稳态计算结果基础上进行再计算，并监控大破口失水事故相关参数，并进行显示和分析；
将获得的安注箱液位值，在已建立模型的瞬态卡信息中，并对已有瞬态模型进行更新，并在稳态计算结果的基础上，进行再计算；
实时监控大破口失水事故相关参数计算值，提取如反应堆功率、稳压器压力、堆芯水位、安注流量、破口流量、包壳峰值温度等结果参数，利用matplotlib库完成结果图的自动绘制，基于文档模板，自动生成分析报告。
如图1所述，一种大破口失水事故分析系统，该系统包括参数抽样模块、稳态计算模块、瞬态计算模块以及硼浓度监测模块，其中，参数抽样模块与稳态计算模块相互连接，并通过参数抽样模块可将采集的核电站大破口失水事故相关参数值形成输入文本后，输入稳态计算模块进行稳态计算，其中，参数抽样模块可采集对失水事故下PCT影响较大的参数，包括但不限于反应堆功率、反应堆平均温度、系统压力、安注箱水体积、安注箱压力、安注箱水温、安注水温度、热流密度热通道因子FQ、核焓升热通道因子FΔH、堆芯上部功率分布、堆芯中部功率分布、破口面积、厂外电是否丧失等相关参数；同时，参数抽样模块可对上述采集的参数进行数据验证，其中，采集的参数除热流密度热通道因子FQ、核焓升热通道因子FΔH两参数为正态分布外，其它参数均需满足均匀分布，如某项参数不满足该项参数的分布要求，则需要进行重新抽样；针对上述参数进行抽样并获得各参数的样本值，并按照工况进行输出形成输入文本；
稳态计算模块根据参数抽样模块传输的包含参数值的输入文本，根据每个参数的物理含义，将每个抽样工况获得的参数值写入文本卡对应的位置，形成稳态计算模型；对文本卡中的参数进行稳态计算，对稳态计算后的参数进行核实，并在获得参数稳态计算值与抽样获得的参数值的偏差大于设定的阈值时，向参数抽样模块发送相应参数重新抽样请求信号；稳态计算模块对稳态计算输出文件进行识别，并根据稳态计算结束时间与预设时间的一致性来判断稳态计算是否完成，并在未完成稳态计算时，则向参数抽样模块发生修改输入文本卡信息的请求信息，并利用稳态计算模块进行稳态计算；
稳态计算模块还与瞬态计算模块相互连接，稳态计算模块在完成稳定计算后，将稳态计算获得的二回路蒸汽流量参数值，以及稳态计算时的破口面积和厂外电丧失情况信息传输至瞬态计算模块，并进行瞬态计算，并根据瞬态计算结果与实际抽样获得的破口、厂外点情况信息进行校验，若出现偏离，瞬态计算模块则向稳态计算模块发送请求，在进行参数修改后，继续进行瞬态计算；
硼浓度监测模块与瞬态计算模块相互连接，利用硼浓度监测模块提取瞬态计算模块计算后的硼浓度值，并当硼浓度值超过所设定的阈值时，获得此刻水容积差值，并计算相应时刻安注箱内剩余水容积值，并根据球形安注箱容积及水位的关系，获得液位值；将该液位值传输至瞬态计算模块，并更新瞬态计算模块中的瞬态信息后进行再计算；数据显示模块与瞬态计算模块相连接，通过数据显示模块可将传输至数据显示模块的核电站大破口失水事故相关参数值进行显示，并利用文档模板库，直接将获得的参数引入生成分析报告。
Claims (15)
1.一种大破口失水事故分析方法，其特征在于，该方法包括：
选取核电站大破口失水事故相关的指标参数，并对各参数进行样本抽样后进行参数分布检查，以验证参数样本值的有效性；
建立稳态计算模型，对大破口失水事故相关参数进行稳态计算，根据稳态计算值与参数抽样值的偏差值是否超过设定阈值，来进行稳态计算后的参数值校验；
构建瞬态计算模型，并对大破口失水事故相关参数进行瞬态计算，并根据计算结果与实际的破口、核电厂外电情况进行对比，并在出现偏离时，重新构建瞬态计算模型并进行瞬态计算；
进行安注水扣除，获得液位参数值，并对瞬态模型进行更新后，在稳态计算结果基础上进行再计算，并在完成计算后，提取大破口失水事故相关参数，并进行参数显示及分析。
2.根据权利要求1所述的一种大破口失水事故分析方法，其特征在于，所述核电站大破口失水事故相关的指标参数包括核电站大破口失水事故相关的保守性假设、初始条件和边界条件的参数、参数范围和分布形式构成的指标参数集。
3.根据权利要求1所述的一种大破口失水事故分析方法，其特征在于，对抽样的参数样本值有效性进行验证是通过对抽样后的参数分布进行检查从而确定抽样参数样本的有效性，其中，热流密度热通道因子FQ和核焓升热通道因子FΔH两参数抽样值需要满足正态分布；反应堆功率、反应堆平均温度、系统压力、安注箱水体积、安注箱压力、安注箱水温、安注水温度、堆芯上部功率分布、堆芯中部功率分布、破口面积、厂外电是否丧失，上述参数需要满足均匀分布。
4.根据权利要求1所述的一种大破口失水事故分析方法，其特征在于，构建稳态计算模型是通过获得的各工况下的参数抽样样本，根据各个参数的物理含义，利用python的正则表达式对关键词进行定位，将每个抽样工况获得的参数值写入输入文本卡对应的位置，完成稳态计算模型构建。
5.根据权利要求1所述的一种大破口失水事故分析方法，其特征在于，对稳态计算输出文件校验是通过对稳态计算输出文件进行识别，并判断稳态计算结束时间与预设时间是否相同，从而判断稳态计算是否完成而进行校验，若未完成稳态计算，则返回修改抽样形成的文本卡信息，并重新进行稳态计算。
6.根据权利要求1所述的一种大破口失水事故分析方法，其特征在于，所述瞬态计算模型是在完成稳态计算后，通过将稳态计算获得的二回路蒸汽流量参数值及抽样获得的破口面积、厂外电丧失情况信息写入瞬态卡中，形成瞬态计算模型。
7.根据权利要求1所述的一种大破口失水事故分析方法，其特征在于，所述液位参数值获取是在完成瞬态计算后，提取计算结果中的硼浓度值，并不断检测核电厂模型中下腔室内的硼浓度，当硼浓度超过所设定的阈值时，获得此时刻的水容积差值，从而计算此时安注箱内剩余水容积，并根据球形安注箱容积及水位的关系，将水容积转化获得液位值。
8.根据权利要求1所述的一种大破口失水事故分析方法，其特征在于，所述瞬态计算或稳态计算是利用python的os模块调用包含附录K模型的系统程序，进行瞬态计算或稳态计算。
9.根据权利要求1所述的一种大破口失水事故分析方法，其特征在于，更新瞬态模型后，在稳态计算结果基础上进行再计算时，是将安注箱液位值更新在已建立模型的瞬态卡信息中，并在稳态计算结果的基础上进行再计算；所述对大破口失水事故相关参数监控，并进行显示和分析是通过提取大破口失水事故相关参数值后，利用matplotlib库完成结果图的自动绘制，基于文档模板，自动生成分析报告。
10.一种大破口失水事故分析系统，其特征在于，该系统包括参数抽样模块、稳态计算模块、瞬态计算模块以及硼浓度监测模块，其中，所述参数抽样模块与所述稳态计算模块相互连接，并通过所述参数抽样模块可将采集的核电站大破口失水事故相关参数值形成输入文本后，输入所述稳态计算模块进行稳态计算；所述稳态计算模块将稳态计算值传输至相互连接的所述瞬态计算模块进行瞬态计算；所述硼浓度监测模块与所述瞬态计算模块相互连接，并通过所述硼浓度监测模块从瞬态计算模型提取并获得安注箱液位值后，对所述瞬态计算模型中的瞬态计算进行更新再计算，并利用与所述瞬态计算模块相连接的数据显示模块对相关参数进行显示和分析。
11.根据权利要求10所述的一种大破口失水事故分析系统，其特征在于，所述参数抽样模块可采集对失水事故下PCT影响较大的参数，包括反应堆功率、反应堆平均温度、系统压力、安注箱水体积、安注箱压力、安注箱水温、安注水温度、热流密度热通道因子FQ、核焓升热通道因子FΔH、堆芯上部功率分布、堆芯中部功率分布、破口面积、厂外电是否丧失；所述参数抽样模块可对采集的抽样参数进行数据验证，其中，采集的参数除热流密度热通道因子FQ、核焓升热通道因子FΔH两参数为正态分布外，其它参数均需满足均匀分布，如某项参数不满足该项参数的分布要求，则需要进行重新抽样。
12.根据权利要求10所述的一种大破口失水事故分析系统，其特征在于，所述参数抽样模块可采集核电厂工况下的各参数样本值，并形成输入文本卡；所述稳态计算模块根据所述输入文本卡，将每个抽样工况获得的参数值写入文本卡对应的位置，形成稳态计算模型；对所述文本卡中的参数进行稳态计算，对稳态计算后的参数进行核实，并在获得参数稳态计算值与抽样获得的参数值的偏差大于设定的阈值时，向参数抽样模块发送相应参数重新抽样请求信号。
13.根据权利要求12所述的一种大破口失水事故分析系统，其特征在于，所述稳态计算模块对稳态计算输出文件进行识别，并根据稳态计算结束时间与预设时间的一致性来判断稳态计算是否完成，并在未完成稳态计算时，则向所述参数抽样模块发生修改输入文本卡信息的请求信息，并利用所述稳态计算模块进行稳态计算。
14.根据权利要求10所述的一种大破口失水事故分析系统，其特征在于，所述瞬态计算模块将稳态计算模块稳态计算后的二回路蒸汽流量参数值及稳态计算时的破口面积和厂外电丧失情况信息输入瞬态卡，并进行瞬态计算，并根据瞬态计算结果与实际抽样获得的破口、厂外点情况信息进行校验，若出现偏离，瞬态计算模块则向稳态计算模块发送请求，在进行参数修改后，继续进行瞬态计算。
15.根据权利要求10所述的一种大破口失水事故分析系统，其特征在于，所述硼浓度监测模块可提取所述瞬态计算模块计算后的硼浓度值，并当硼浓度值超过所设定的阈值时，获得此刻水容积差值，并计算相应时刻安注箱内剩余水容积值，并根据球形安注箱容积及水位的关系，获得液位值；将该液位值传输至瞬态计算模块，并更新瞬态计算模块中的瞬态信息后进行再计算。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110647075.9A CN113536537B (zh) | 2021-06-10 | 2021-06-10 | 一种大破口失水事故分析方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110647075.9A CN113536537B (zh) | 2021-06-10 | 2021-06-10 | 一种大破口失水事故分析方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113536537A CN113536537A (zh) | 2021-10-22 |
| CN113536537B true CN113536537B (zh) | 2024-01-12 |
Family
ID=78124793
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110647075.9A Active CN113536537B (zh) | 2021-06-10 | 2021-06-10 | 一种大破口失水事故分析方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113536537B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114117792B (zh) * | 2021-11-26 | 2023-09-26 | 中国核动力研究设计院 | 一种失水事故分析方法 |
Citations (5)
| 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 |
| CA2818881A1 (en) * | 2010-11-23 | 2012-05-31 | Westinghouse Electric Company Llc | Full spectrum loca evaluation model and analysis methodology |
| CN104700222A (zh) * | 2015-03-18 | 2015-06-10 | 中科华核电技术研究院有限公司 | 核电厂大破口事故不确定性分析方法 |
| CN105653869A (zh) * | 2016-01-05 | 2016-06-08 | 中国核动力研究设计院 | 一种超临界水堆堆芯三维瞬态性能分析方法 |
| CN107644694A (zh) * | 2017-09-20 | 2018-01-30 | 岭东核电有限公司 | 一种核电厂大破口事故分析方法 |
-
2021
- 2021-06-10 CN CN202110647075.9A patent/CN113536537B/zh active Active
Patent Citations (5)
| 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 |
| CA2818881A1 (en) * | 2010-11-23 | 2012-05-31 | Westinghouse Electric Company Llc | Full spectrum loca evaluation model and analysis methodology |
| CN104700222A (zh) * | 2015-03-18 | 2015-06-10 | 中科华核电技术研究院有限公司 | 核电厂大破口事故不确定性分析方法 |
| CN105653869A (zh) * | 2016-01-05 | 2016-06-08 | 中国核动力研究设计院 | 一种超临界水堆堆芯三维瞬态性能分析方法 |
| CN107644694A (zh) * | 2017-09-20 | 2018-01-30 | 岭东核电有限公司 | 一种核电厂大破口事故分析方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113536537A (zh) | 2021-10-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5357389B2 (ja) | 原子炉運転に関する作動限界までのマージンを判定する方法 | |
| CN103973513B (zh) | 一种核电站全范围模拟机的远程监测方法和系统 | |
| US20220291654A1 (en) | Online monitoring method of nuclear power plant system based on isolation forest method and sliding window method | |
| CN107767025A (zh) | 一种风险指引的核电厂大破口失水事故分析方法 | |
| CN105006259A (zh) | 一种核电厂堆芯损伤评价方法 | |
| CN113536537B (zh) | 一种大破口失水事故分析方法及系统 | |
| CN112133460B (zh) | 一种快堆堆芯在线监督方法及系统 | |
| CN110739091A (zh) | 一种核电厂次临界度直接测量系统及硬件平台和测量方法 | |
| Franzén | Evaluation of Fuel Assembly Bow Penalty Peaking Factors for Ringhals 3: Based on a Cycle Specific Core Water Gap Distribution | |
| Mandelli et al. | Overview of new tools to perform safety analysis: BWR station black out test case | |
| CN112926686A (zh) | 基于brb和lstm模型的电力大数据用电异常检测方法及装置 | |
| CN116680905A (zh) | 针对反应堆严重事故不确定性与敏感性分析的方法及系统 | |
| Kim | Rigorous derivation of interfacing system LOCA frequency formulas for probabilistic safety assessment of nuclear power plants | |
| CN115564247A (zh) | 一种核电厂事故智能识别及决策的方法及系统 | |
| Wiltbank | Dynamic probabilistic risk assessment using ADAPT/RELAP5-3D for cybersecurity threat analysis | |
| Tian et al. | Application of Best-Estimate Plus Uncertainty Analysis Method in Nuclear Safety Evaluation | |
| CN113868959B (zh) | 一种基于自适应抽样与det结合的核电厂事故失效概率计算方法 | |
| KR102675124B1 (ko) | 중수로 신규연료 시범장전 중 연료특성 추적계산 방법 및 시스템 | |
| Fisher et al. | Application of Boron Feedback From BOA CIPS/CILC Analysis on Neutronic Models | |
| Karnaukhov et al. | Analysis of the sensitivity and uncertainty of calculation using the KORSAR/GP code for accidents with a failure of storage pools’ cooling system | |
| Yu et al. | Simulation Approach for Safety Evaluation of Smart Transmitters in NPPs | |
| CN115906267A (zh) | 故障分析、检测方法、系统、设备和介质 | |
| Cao et al. | Digital Modeling and Simulation Technology of Pumping Infrastructure Engineering | |
| CN113837535A (zh) | 一种用于核电厂严重事故进程回溯的方法 | |
| Iooss et al. | Uncertainty assessments in severe nuclear accident scenarios |
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 |