CN104951648B - 用于估算核电站严重事故后氧气浓度的方法 - Google Patents
用于估算核电站严重事故后氧气浓度的方法 Download PDFInfo
- Publication number
- CN104951648B CN104951648B CN201510268862.7A CN201510268862A CN104951648B CN 104951648 B CN104951648 B CN 104951648B CN 201510268862 A CN201510268862 A CN 201510268862A CN 104951648 B CN104951648 B CN 104951648B
- Authority
- CN
- China
- Prior art keywords
- hydrogen
- containment
- oxygen
- oxygen concentration
- calculated
- 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
Landscapes
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明涉及核电站安全领域，特别涉及一种用于估算核电站严重事故后氧气浓度的方法。该方法基于目前压水堆核电站普遍采用的非能动氢气复合器系统和氢气浓度监测系统，采用简化的计算方法，用来估算严重事故后安全壳内的氧气浓度，为核电站发生严重事故后的氢气风险评价和严重事故管理提供依据。
Description
技术领域
本发明涉及核电站安全领域，特别涉及一种用于估算核电站严重事故后氧气浓度的快速计算方法。
背景技术
核电站发生严重事故后，堆芯燃料包壳的金属锆会和高温的水蒸气发生锆水反应产生大量氢气。压力容器失效之后，堆芯熔融物和安全壳底板发生堆芯熔融物-混凝土(MCCI)反应，也会持续的释放出氢气。氢气在安全壳内的积聚到一定程度，有可能发生燃烧、爆炸，氢气爆炸带来的冲击波会威胁安全壳内的完整性。为了防止由于氢气爆炸导致的安全壳失效，核电站普遍设置了氢气控制系统和氢气监测系统，分别来消除安全壳内的氢气和监测氢气浓度，为严重事故管理提供依据。
随着氢气控制系统的工作，安全壳内的氧气也会随着不断消耗，氧气浓度不断降低。虽然氢气不断产生，在没有氧气的情况下，仍不会产生氢气爆炸风险。因此，估算严重事故后安全壳内的氧气浓度对严重事故后氢气风险的判断有一定意义。但由于核电站发生严重事故后，安全壳内放射性水平较高，无法将安全壳内的气体取样测量其氧气含量，因此给严重事故管理带来了盲区。
目前国内二代、二代加核电站(如岭澳核电站、秦山二期核电站)和三代核电站(如EPR、ACP1000)普遍采用非能动氢气复合器(详见“一种无须外部提供能源的氢与氧复合装置”，专利号：200620116543.0)来控制安全壳内的氢气浓度，使安全壳内的氢气浓度维持在10％以下，防止发生氢气爆炸。非能动氢气复合器是一种非能动的装置，采用催化复合的方法来使氢气和氧气发生反应，从而达到消除氢气的目的。氢气复合器的外形结构见附图1，防护壳2设在底座5上，底座5上设有进风孔4，在防护壳2的上端侧壁上设有出风口1，防护壳外壁设有固定架3，催化床为无底的抽屉式结构，抽屉6内设有催化剂金属网板8，抽屉6外壁固定有把手7。氢气复合器的消氢速率与安全壳内的温度、压力、氢气氧气浓度有关。
为了监测严重事故后安全壳内的氢气浓度，通常采用氢气监测装置来测量氢气浓度。其中一种氢气监测装置采用氢气催化复合的方法(详见“一种测量核电厂安全壳内氢气浓度的方法”，发明申请号201210400679.4)，其原理与非能动氢气复合器一致，通过使氢气氧气发生催化复合反应，测量核电厂环境温度T1和化合反应温度T2，根据温升计算氢气浓度。该方法流程见图2。
发明内容
本发明的目的在与针对目前严重事故下安全壳内氧气浓度监测手段的不足，提供了一种用于估算核电站严重事故后氧气浓度的方法,该方法基于目前压水堆核电站普遍采用的非能动氢气复合器系统和氢气浓度监测系统，采用简化的计算方法，用来估算严重事故后安全壳内的氧气浓度，为核电站发生严重事故后的氢气风险评价和严重事故管理提供依据。
为达到上述目的，本发明提供了一种用于估算核电站严重事故后氧气浓度的方法，包括以下步骤：
(S101)根据安全壳体积、正常运行压力、温度计算氧气初始总摩尔数NO20；
(S102)根据氢气测量系统测量得到的氢气浓度XH2推算由于非能动氢气复合器的复合反应导致的氧气消耗速率N’O2；
(S103)根据氧气消耗速率N’O2计算安全壳剩余氧气摩尔数N02；
(S104)根据安全壳内温度T、压力P和安全壳剩余氧气摩尔数N02计算氧气浓度XO2。
进一步，如上所述的用于估算核电站严重事故后氧气浓度的方法，步骤(S101)中，在计算氧气初始总摩尔数NO20时，按照安全壳大气100％湿度和理想气体模型计算。
进一步，如上所述的用于估算核电站严重事故后氧气浓度的方法，每隔一段时间重复步骤(S102)～步骤(S104)，对氧气浓度进行一次计算，第n次计算与第n-1次计算的时间差为时间步长DT(n)。
进一步，如上所述的用于估算核电站严重事故后氧气浓度的方法，第n次计算时氧气消耗速率根据氢气测量系统测量得到的氢气浓度XH2(n)乘以系数A确定，即N’O2(n)＝A×XH2(n)，
其中，系数A根据安全壳内非能动氢气复合器的数目和消氢能力计算得到。若核电厂有氢气复合器台数为k，在氢气浓度为X0时，每台氢气复合器的消氢速率为M,则系数A＝0.5×k×M/X0。
进一步，如上所述的用于估算核电站严重事故后氧气浓度的方法，第n次计算时安全壳剩余氧气摩尔数通过以下公式计算：
进一步，如上所述的用于估算核电站严重事故后氧气浓度的方法，在计算第n步氧气浓度XO2(n)时，通过以下公式计算：
其中，P(n)为第n步计算时安全壳压力，T(n)为第n步计算时安全壳温度，P(n)、T(n)分别由安全壳压力、温度监测仪表得到，V为安全壳体积，R为气体常量。
本发明的有益效果如下：通过本发明所提供的计算方法，在核电站发生严重事故后，在无法获得安全壳内氧气浓度的情况下，可以方便的估算出核电站严重事故后安全壳内的氧气浓度，从而能够判断出安全壳内有无氢气燃烧、爆炸风险，为严重事故管理、氢气风险的判断提供依据。
附图说明
图1为目前国内外普遍采用的非能动氢气复合器结构图；
图2为目前国内采用的氢气测量装置原理流程图；
图3为本发明实施例中严重事故后氧气浓度估算的简易计算方法流程图；
图4为本发明实施例中计算氧气初始总摩尔数NO20流程图。
具体实施方式
下面详细描述本发明的实施例。所述实施例的流程图在附图中示出。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能解释为对本发明的限制。
如图3所示，为本发明示例的严重事故后氧气浓度估算的简易计算方法流程，包括以下步骤：
步骤(S101)，根据安全壳体积、正常运行压力、温度计算氧气初始总摩尔数NO20。
本发明的实施例中，提出了计算氧气初始总摩尔数的方法，包括以下几个步骤，如图4所示：
步骤(S201)，计算当前温度下的饱和蒸汽压Psteam，此为公知技术；
步骤(S202)，计算当前压力下的空气分压Pair＝P-Psteam，其中P为安全壳正常运行压力；
步骤(S203)，计算氧气初始总摩尔数其中V为安全壳体积，R为气体常量。
步骤(S102)，根据氢气浓度测量得到的氢气浓度XH2推算由于非能动氢气复合器的复合反应导致的氧气消耗速率N’O2；
在本发明的一个实施例中，在计算时每隔一段时间对氧气浓度进行一次估算，即每隔一段时间重复步骤(S102)～步骤(S104)，第n次计算与第n-1次计算的时间差为时间步长DT(n)。
在本发明的一个实施例中，第n次计算时氧气消耗速率根据氢气测量系统测量得到的氢气浓度XH2(n)乘以系数A确定，即N’O2(n)＝A×XH2(n)。其中系数A根据安全壳内氢气复合器的数目和消氢能力计算得到。假如某核电厂有氢气复合器台数为k，在氢气浓度为X0时，每台氢气复合器的消氢速率为M,则系数A＝0.5×k×M/X0。
步骤(S103)，根据氧气消耗速率计算安全壳剩余氧气摩尔数N02；
在本发明的一个实施例中，第n次计算时安全壳剩余氧气摩尔数通过以下公式计算：
步骤S104，根据安全壳内温度T、压力P和安全壳剩余氧气摩尔数计算结果计算氧气浓度XO2。
在本发明的一个实施例中，在计算第n步氧气浓度XO2(n)时，通过以下公式计算：
其中，P(n)为第n步计算时安全壳压力，T(n)为第n步计算时安全壳温度，P(n)、T(n)分别由安全壳压力、温度监测仪表得到，V为安全壳体积，R为气体常量。
通过本发明，在核电站发生严重事故后，在无法获得安全壳内氧气浓度的情况下，可以方便的估算出核电站严重事故后安全壳内的氧气浓度，从而能够判断出安全壳内有无氢气燃烧、爆炸风险，为严重事故管理、氢气风险的判断提供依据。
尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种多样变化、修改、替换和变型，本发明的范围由所附权利要求及其等同限定。
Claims (2)
1.一种用于估算核电站严重事故后氧气浓度的方法,其特征在于，包括以下步骤：
(S101)根据安全壳体积、正常运行压力和温度计算氧气初始总摩尔数NO20，包括：
计算当前温度下的饱和蒸汽压Psteam；
计算当前压力下的空气分压Pair＝P-Psteam，其中P为安全壳正常运行压力；
计算氧气初始总摩尔数其中V为安全壳体积，R为气体常量，T为温度；
(S102)根据氢气测量系统测量得到的氢气浓度XH2推算由于非能动氢气复合器的复合反应导致的氧气消耗速率N’O2；
(S103)根据氧气消耗速率N’O2计算安全壳剩余氧气摩尔数N02，第n次计算时安全壳剩余氧气摩尔数通过以下公式计算：
NO2(n)＝氧气初始总摩尔数
(S104)根据安全壳内温度T、压力P和安全壳剩余氧气摩尔数N02计算氧气浓度XO2，在计算第n步氧气浓度XO2(n)时，通过以下公式计算：
其中，P(n)为第n步计算时安全壳压力，T(n)为第n步计算时安全壳温度，P(n)、T(n)分别由安全壳压力、温度监测仪表得到，V为安全壳体积，R为气体常量；
(S105)每隔一段时间重复步骤(S102)～步骤(S104)，对氧气浓度进行一次计算，第n次计算与第n-1次计算的时间差为时间步长DT(n)；第n次计算时氧气消耗速率根据氢气测量系统测量得到的氢气浓度XH2(n)乘以系数A确定，即N’O2(n)＝A×XH2(n)，
其中，系数A根据安全壳内非能动氢气复合器的数目和消氢能力计算得到，若核电厂有氢气复合器台数为k，在氢气浓度为X0时，每台氢气复合器的消氢速率为M,则系数A＝0.5×k×M/X0。
2.如权利要求1所述的用于估算核电站严重事故后氧气浓度的方法，其特征在于，步骤(S101)中，在计算氧气初始总摩尔数NO20时，按照安全壳大气100％湿度和理想气体模型计算。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510268862.7A CN104951648B (zh) | 2015-05-25 | 2015-05-25 | 用于估算核电站严重事故后氧气浓度的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510268862.7A CN104951648B (zh) | 2015-05-25 | 2015-05-25 | 用于估算核电站严重事故后氧气浓度的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104951648A CN104951648A (zh) | 2015-09-30 |
| CN104951648B true CN104951648B (zh) | 2019-01-04 |
Family
ID=54166302
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510268862.7A Active CN104951648B (zh) | 2015-05-25 | 2015-05-25 | 用于估算核电站严重事故后氧气浓度的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104951648B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107993731B (zh) * | 2017-11-06 | 2021-08-17 | 中国核电工程有限公司 | 一种反应堆严重事故后安全壳内气体可燃性监测系统 |
| CN112699335B (zh) * | 2020-12-21 | 2023-08-29 | 中国辐射防护研究院 | 后处理厂高放废液贮槽氢气爆炸事故释放源项估算方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102928550A (zh) * | 2012-10-19 | 2013-02-13 | 中国船舶重工集团公司第七一八研究所 | 一种测量核电厂安全壳内氢气浓度的方法 |
| CN103033441A (zh) * | 2012-12-18 | 2013-04-10 | 重庆川仪分析仪器有限公司 | 安全壳内气体浓度分析方法、装置及气体分析仪 |
| CN104538070A (zh) * | 2014-12-08 | 2015-04-22 | 中广核工程有限公司 | 核电站安全壳内氢气浓度测量方法、装置及系统 |
-
2015
- 2015-05-25 CN CN201510268862.7A patent/CN104951648B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102928550A (zh) * | 2012-10-19 | 2013-02-13 | 中国船舶重工集团公司第七一八研究所 | 一种测量核电厂安全壳内氢气浓度的方法 |
| CN103033441A (zh) * | 2012-12-18 | 2013-04-10 | 重庆川仪分析仪器有限公司 | 安全壳内气体浓度分析方法、装置及气体分析仪 |
| CN104538070A (zh) * | 2014-12-08 | 2015-04-22 | 中广核工程有限公司 | 核电站安全壳内氢气浓度测量方法、装置及系统 |
Non-Patent Citations (4)
| Title |
|---|
| 先进核电站严重事故下氢浓度监测技术研究;王占元;《仪器仪表用户》;20130808;第20卷(第4期);第44-46,9页 |
| 核电厂事故后安全壳内氢气浓度测量方法分析;刘静;《核安全》;20140630;第13卷(第2期);第2.1.1节 |
| 秦山二期扩建工程严重事故下安全壳内消氢措施的研究;尹峰;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20080815(第8期);表7-1 |
| 秦山二期核电厂严重事故下安全壳内氢气浓度分布及风险初步分析;邓坚等;《核动力工程》;20080415;第29卷(第2期);第78-84页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104951648A (zh) | 2015-09-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6400685B2 (ja) | 原子炉運転停止時におけるホウ素希釈監視方法 | |
| RU2009137958A (ru) | Способ определения трехмерного распределения мощности в активной зоне ядерного реактора | |
| CN104951648B (zh) | 用于估算核电站严重事故后氧气浓度的方法 | |
| CN105895175B (zh) | 一种标定核反应堆堆芯出口热电偶的方法 | |
| JP5424526B2 (ja) | 核軽水炉においてドライアウト特性を推測する方法、コンピュータプログラム製品、原子力エネルギープラント及び運転方法 | |
| CN107464591A (zh) | 一种安全壳泄漏在线监测系统传感器组合优化方法 | |
| WO2024061196A1 (zh) | 一种核电站蒸汽发生器的智能监测方法及系统 | |
| CN109659587A (zh) | 液流电池容量衰减控制系统及方法 | |
| CN208283037U (zh) | 一种储罐泄漏闪蒸射流行为特性实验测试装置 | |
| JP2014055942A (ja) | 電力を用いないで使用済燃料プールの液位を測定する方法およびシステム | |
| CN104332193B (zh) | 一种基于数字化技术的蒸汽发生器水位测量方法 | |
| Ko et al. | Structural analysis and response measurement locations of inner barrel assembly top plate in APR1400 | |
| Hugo et al. | Predicting evaporation rates from spent nuclear fuel storage pools | |
| CN115563762A (zh) | 一种受热面管屏状态的预警方法、装置、设备及存储介质 | |
| CN113639939A (zh) | 一种建筑物内空气泄漏率异常变化的实时监测方法及系统 | |
| CN111028965A (zh) | 一种基于安全壳氢气浓度的燃料元件过热损伤评价方法 | |
| Auban et al. | Experimental investigation of natural-circulation flow behavior under low-power/low-pressure conditions in the large-scale PANDA facility | |
| US10460847B2 (en) | Gravity-based, non-invasive reactor system and method for coolant inventory monitoring | |
| Lee et al. | Prediction of SG tube support plate flow area blockage rate using SG wide range level measurements and hydrodynamic analysis | |
| Nema et al. | Dynamics of hydrogen in sodium in LMFBR secondary circuit | |
| Kim | A study for modeling electrochemistry in light water reactors | |
| Proskuriakov et al. | Identification of standing pressure waves sources in primary loops of NPP with WWER and PWR | |
| McKee et al. | Use of in-sodium hydrogen detectors during steam generator tests in the SCTI | |
| Gardner et al. | Development of a hydrogen management concept for the Canadian supercritical water-cooled reactor | |
| Ye et al. | Research on the Influence of Degassing on the Release of Fission Gases in PWRs |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |