CN112488465B - 一种八氧化三铀临界事故应急时总裂变次数估算方法 - Google Patents

一种八氧化三铀临界事故应急时总裂变次数估算方法 Download PDF

Info

Publication number
CN112488465B
CN112488465B CN202011270525.9A CN202011270525A CN112488465B CN 112488465 B CN112488465 B CN 112488465B CN 202011270525 A CN202011270525 A CN 202011270525A CN 112488465 B CN112488465 B CN 112488465B
Authority
CN
China
Prior art keywords
critical
triuranium octoxide
estimating
accident
gamma
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
Application number
CN202011270525.9A
Other languages
English (en)
Other versions
CN112488465A (zh
Inventor
贾林胜
杨亚鹏
王任泽
冯宗洋
王宁
梁博宁
张建岗
李国强
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Institute for Radiation Protection
Original Assignee
China Institute for Radiation Protection
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Institute for Radiation Protection filed Critical China Institute for Radiation Protection
Priority to CN202011270525.9A priority Critical patent/CN112488465B/zh
Publication of CN112488465A publication Critical patent/CN112488465A/zh
Application granted granted Critical
Publication of CN112488465B publication Critical patent/CN112488465B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0637Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00Energy generation of nuclear origin
    • Y02E30/30Nuclear fission reactors

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Educational Administration (AREA)
  • Mathematical Physics (AREA)
  • Operations Research (AREA)
  • Development Economics (AREA)
  • General Business, Economics & Management (AREA)
  • Game Theory and Decision Science (AREA)
  • Quality & Reliability (AREA)
  • Marketing (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Tourism & Hospitality (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Measurement Of Radiation (AREA)

Abstract

本发明涉及一种八氧化三铀临界事故应急时总裂变次数估算方法,所述方法包括以下步骤:(1)确定八氧化三铀系统是否发生临界事故;(2)确定γ临界报警数据能否获取且八氧化三铀为高富集度,则采用NUREG‑CR/6504推荐的方法;(3)若γ临界报警数据未能获取,或八氧化三铀为非高富集度,则采用RASCAL推荐的方法;(4)判断临界是否结束,若未结束,当前获取的临界信息一旦有更新,再从步骤(1)开始依次执行,估算临界裂变次数。本发明提供的方法与事故进行相结合,为应急决策者提供更准确的技术支持。

Description

一种八氧化三铀临界事故应急时总裂变次数估算方法
技术领域
本发明属于核临界事故裂变次数估算技术领域,具体涉及一种八氧化三铀临界事故应急时总裂变次数估算方法。
背景技术
核临界事故应急评价技术是核燃料循环设施应急准备与响应的重要组成部分,且我国应急管理也要求,对存在潜在核临界事故的设施需要具备该事故的应急评价能力。八氧化三铀临界事故是后处理厂可能发生的事故,需要做好相应的应急工作。临界总裂变次数是核临界事故应急评价的重要内容,也是技术难点之一。它反映了核临界事故的大小和规模,直接影响事故应急防护行动决策。
发生临界事故时,往往开始获取的数据有限,随着时间的推移,获取的数据越来越丰富,因此总裂变次数估算非常有必要与事故进程相结合,建立一套基于进程信息的估算方法,为应急决策者提供更有力的技术支持。
发明内容
针对现有技术中所存在的问题,本发明的目的在于提供一种与事故进程相结合的八氧化三铀临界事故应急时总裂变次数估算方法,为应急决策者提供更有力的技术支持。
为达到以上目的,本发明采用的一种技术方案是:一种八氧化三铀临界事故应急时总裂变次数估算方法,所述方法与事故进程相结合,所述方法包括以下步骤:
(1)确定八氧化三铀系统是否发生临界事故;
(2)如果发生临界事故,确定能否获取临界γ报警仪读数,如果能获取读数,且八氧化三铀为高富集度,则采用NUREG-CR/6504中基于U(93.2)3O8@H/235U=10情景下γ剂量率与时间和距离的关系曲线图的估算方法;
(3)若γ临界报警数据未能获取,或八氧化三铀为非高富集度,则采用RASCAL基于系统情景的估算方法;
(4)判断临界是否结束,若未结束,当前获取的临界信息是否有更新,一旦有更新,再从步骤(1)开始依次执行,估算临界裂变次数。
进一步的,所述步骤(2)中NUREG-CR/6504中基于
U(93.2)3O8@H/235U=10情景下γ剂量率与时间和距离的关系曲线图的估算方法具体步骤为:
a)根据NUREG-CR/6504中基于U(93.2)3O8@H/235U=10情景下γ剂量率与时间和距离的关系曲线图建立八氧化三铀临界裂变次数为
1017时,γ剂量率随时间、距离变化的拟合公式;
b)根据所述拟合公式估算出裂变次数为1017时,不同时间、不同临界报警仪与事故点距离对应的γ剂量率;
c)计算实际临界事故总裂变次数,计算公式为:
Figure BDA0002777566920000021
式中:
D实际-实际临界γ报警仪读数;
D估算-估算的γ剂量率。
进一步的,所述步骤(3)中RASCAL基于系统情景的估算方法为:当八氧化三铀为粉末时,估算首次裂变次数为3×1020,总裂变次数为3×1020;当八氧化三铀为块状铀时,估算首次裂变次数为3×1019,总裂变次数为3×1019
本发明的有益效果在于:1、本方法与事故进程相结合,根据不同的进程和所获得的数据,选择合适的估算方法,为应急决策者提供更准确的技术支持。2、建立了八氧化三铀临界裂变次数为1017下,γ剂量率随时间、距离变化的拟合公式,较NUREG-CR/6504中从曲线图中读取剂量率的方法更方便、准确。
附图说明
图1为本发明所述一种八氧化三铀临界事故应急时总裂变次数估算方法的流程图;
图2为NUREG/CR-6504给出的U(93.2)3O8@H/235U=10情景下γ剂量率与时间距离的关系。
具体实施方式
为使本发明解决的技术问题、采用的技术方案和达到的技术效果更加清楚,下面将结合附图对本发明实施例的技术方案作进一步的详细描述。显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,均属于本发明保护的范围。
参阅图1,图1为本发明所述方法一实施例的流程示意图,所述方法包括以下步骤:
步骤(1):确定八氧化三铀系统是否发生临界事故。γ临界报警仪会设置阈值,当读数超过该值会发出声光报警,表示八氧化三铀系统已发生临界事故。
步骤(2):如果发生临界事故,确定能否获取临界γ报警仪读数,如果能获取读数,且八氧化三铀为高富集度,则采用NUREG-CR/6504中基于U(93.2)3O8@H/235U=10情景下γ剂量率与时间和距离的关系曲线图的估算方法,具体计算步骤为:
a)参阅图2,图2为U(93.2)3O8@H/235U=10情景下γ剂量率与时间距离的关系曲线图,根据图2建立八氧化三铀临界裂变次数为1017时,γ剂量率随时间、距离变化的拟合公式,拟合公式见表1;
表1八氧化三铀γ剂量率随时间、距离变化的拟合公式表
Figure BDA0002777566920000041
式中:D-估算的剂量率,x-时间
b)根据所述拟合公式估算出裂变次数为1017时,不同时间、不同临界报警仪与事故点距离对应的γ剂量率。例如,当临界发生后10s,临界报警γ剂量率监测仪表距事故点的距离为1.52m,则采用拟合公式D=280000/x0.7068估算γ剂量率;
c)根据估算的γ剂量率与实际临界γ报警仪读数,计算实际临界事故总裂变次数,计算公式为:
Figure BDA0002777566920000042
式中:
D实际-实际临界γ报警仪读数;
D估算-估算的γ剂量率。
步骤(3):若γ临界报警数据未能获取,或八氧化三铀为非高富集度,则采用RASCAL基于系统情景的估算方法。
一般一次临界发生的瞬间,临界报警仪表会报警,但可能获取不到数据,因为剂量太高仪表堵死,但临界发生后数十秒或几分钟以后仪表会显示数据,这与仪器本身有关。当不能获取临界γ报警仪读数,或不满足高富集条件时,则采用RASCAL基于情景的估算方法,具体方法参阅表2。
表2不同系统情景对应的裂变次数
系统情景 首次脉冲裂变次数 总裂变次数
粉末 3×10<sup>20</sup> 3×10<sup>20</sup>
块状铀 3×10<sup>19</sup> 3×10<sup>19</sup>
步骤(4):判断临界是否结束,若未结束,当前获取的临界信息是否有更新,一旦有更新,再从步骤(1)开始依次执行,估算临界裂变次数。
本领域技术人员应该明白,本发明所述装置及方法并不限于具体实施方式中所述的实施例,上面的具体描述只是为了解释本发明的目的,并非用于限制本发明。本领域技术人员根据本发明的技术方案得出其他的实施方式,同样属于本发明的技术创新范围,本发明的保护范围由权利要求及其等同物限定。

Claims (5)

1.一种八氧化三铀临界事故应急时总裂变次数估算方法,其特征在于,所述方法与事故进程相结合,所述方法包括以下步骤:
(1)确定八氧化三铀系统是否发生临界事故;
(2)如果发生临界事故,确定能否获取临界γ报警仪读数,如果能获取读数,且八氧化三铀为高富集度,则采用NUREG-CR/6504中基于U(93.2)3O8@H/235U=10情景下γ剂量率与时间和距离的关系曲线图的估算方法;
(3)若γ临界报警数据未能获取,或八氧化三铀为非高富集度,则采用RASCAL基于系统情景的估算方法;
(4)判断临界是否结束,若未结束,当前获取的临界信息是否有更新,一旦有更新,再从步骤(1)开始依次执行,估算临界裂变次数;
所述NUREG-CR/6504中基于U(93.2)3O8@H/235U=10情景下γ剂量率与时间和距离的关系曲线图的估算方法具体步骤为:
a)根据NUREG-CR/6504中基于U(93.2)3O8@H/235U=10情景下γ剂量率与时间和距离的关系曲线图建立八氧化三铀临界裂变次数为1017时,γ剂量率随时间、距离变化的拟合公式;
b)根据所述拟合公式估算出裂变次数为1017时,不同时间、不同临界报警仪与事故点距离对应的γ剂量率;
c)计算实际临界总裂变次数,计算公式为:
Figure FDA0003774938900000011
式中:
D实际-实际临界γ报警仪读数;
D估算-估算的γ剂量率。
2.根据权利要求1所述的一种八氧化三铀临界事故应急时总裂变次数估算方法,其特征在于,所述RASCAL基于系统情景的估算方法为:当八氧化三铀为粉末时,估算首次裂变次数为3×1020
3.根据权利要求1所述的一种八氧化三铀临界事故应急时总裂变次数估算方法,其特征在于,所述RASCAL基于系统情景的估算方法为:当八氧化三铀为粉末时,估算总裂变次数为3×1020
4.根据权利要求1所述的一种八氧化三铀临界事故应急时总裂变次数估算方法,其特征在于,所述RASCAL基于系统情景的估算方法为:当八氧化三铀为块状铀时,估算首次裂变次数为3×1019
5.根据权利要求1所述的一种八氧化三铀临界事故应急时总裂变次数估算方法,其特征在于,所述RASCAL基于系统情景的估算方法为:当八氧化三铀为块状铀时,估算总裂变次数为3×1019
CN202011270525.9A 2020-11-13 2020-11-13 一种八氧化三铀临界事故应急时总裂变次数估算方法 Active CN112488465B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011270525.9A CN112488465B (zh) 2020-11-13 2020-11-13 一种八氧化三铀临界事故应急时总裂变次数估算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011270525.9A CN112488465B (zh) 2020-11-13 2020-11-13 一种八氧化三铀临界事故应急时总裂变次数估算方法

Publications (2)

Publication Number Publication Date
CN112488465A CN112488465A (zh) 2021-03-12
CN112488465B true CN112488465B (zh) 2022-11-18

Family

ID=74930725

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011270525.9A Active CN112488465B (zh) 2020-11-13 2020-11-13 一种八氧化三铀临界事故应急时总裂变次数估算方法

Country Status (1)

Country Link
CN (1) CN112488465B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102306508A (zh) * 2011-09-07 2012-01-04 中国核动力研究设计院 可用于反应堆临界监督和临界外推的反应性仪
JP2013003104A (ja) * 2011-06-21 2013-01-07 Toshiba Corp 未臨界度測定装置および未臨界度測定方法
CN109657905A (zh) * 2018-11-12 2019-04-19 中国辐射防护研究院 一种mox燃料的临界事故环境释放源项的估算方法
CN110390444A (zh) * 2018-04-16 2019-10-29 中国辐射防护研究院 一种核燃料循环设施uf6事故泄漏后果评估计算方法
CN111798643A (zh) * 2020-05-25 2020-10-20 中国原子能科学研究院 一种新型全局临界事故报警方法及系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013003104A (ja) * 2011-06-21 2013-01-07 Toshiba Corp 未臨界度測定装置および未臨界度測定方法
CN102306508A (zh) * 2011-09-07 2012-01-04 中国核动力研究设计院 可用于反应堆临界监督和临界外推的反应性仪
CN110390444A (zh) * 2018-04-16 2019-10-29 中国辐射防护研究院 一种核燃料循环设施uf6事故泄漏后果评估计算方法
CN109657905A (zh) * 2018-11-12 2019-04-19 中国辐射防护研究院 一种mox燃料的临界事故环境释放源项的估算方法
CN111798643A (zh) * 2020-05-25 2020-10-20 中国原子能科学研究院 一种新型全局临界事故报警方法及系统

Also Published As

Publication number Publication date
CN112488465A (zh) 2021-03-12

Similar Documents

Publication Publication Date Title
CN107301570B (zh) 业务量预测方法、异常业务量检测方法、装置和电子设备
EP3293737B1 (en) Correlation tolerance limit setting system using repetitive cross-validation and method therefor
CN110739091B (zh) 一种核电厂次临界度直接测量系统及硬件平台和测量方法
CN110334406A (zh) 一种考虑风速特大值的极值风速重现期确定方法和装置
CN112488465B (zh) 一种八氧化三铀临界事故应急时总裂变次数估算方法
CN113095579B (zh) 一种耦合伯努利-伽马-高斯分布的日尺度降水预报校正方法
CN112488461B (zh) 一种金属铀临界事故应急时总裂变次数估算方法
CN112488460B (zh) 一种二氧化铀临界事故应急时总裂变次数估算方法
CN108801320B (zh) 一种天然气测量系统的诊断方法
CN116499609A (zh) 核电厂波动瞬态的自动识别方法、系统、计算机设备和存储介质
CN112488459B (zh) 一种硝酸铀酰溶液临界事故应急时总裂变次数估算方法
Egozcue et al. Wave-height hazard analysis in Eastern Coast of Spain-Bayesian approach using generalized Pareto distribution
JP2008157669A (ja) 未臨界度測定方法及び未臨界度測定用プログラム、並びに未臨界度測定装置
CN107562695B (zh) 一种配电变压器的负载数据处理方法及装置
JPH04232500A (ja) 原子炉出力モニタリング方法及び装置
CN113435464B (zh) 异常数据检测方法、装置、电子设备及计算机存储介质
CN110244096B (zh) 一种电能量计量系统中对电表满码自动发现与处理的方法
CN108268694B (zh) 核事故辐射数据的模拟方法、装置、终端设备和存储介质
CN110020000B (zh) 判定异常风速数据的方法和装置
JP2004294193A (ja) 放射線出力監視装置
CN114334194B (zh) 高温气冷堆氦气泄漏预警方法、装置、设备及存储介质
CN206095136U (zh) 一种混凝土膨胀率/收缩率自动测试装置
CN111176931A (zh) 一种运营监控方法、装置、服务器及存储介质
CN112488464B (zh) 一种钚溶液临界事故应急时总裂变次数估算方法
CN111079267B (zh) 一种基于微观残余应力突度的装载机底盘可靠性检测方法

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