CN108268694B - 核事故辐射数据的模拟方法、装置、终端设备和存储介质 - Google Patents
核事故辐射数据的模拟方法、装置、终端设备和存储介质 Download PDFInfo
- Publication number
- CN108268694B CN108268694B CN201711248992.XA CN201711248992A CN108268694B CN 108268694 B CN108268694 B CN 108268694B CN 201711248992 A CN201711248992 A CN 201711248992A CN 108268694 B CN108268694 B CN 108268694B
- Authority
- CN
- China
- Prior art keywords
- radiation
- reactor
- nuclear
- atmospheric region
- region
- 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
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
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21C—NUCLEAR REACTORS
- G21C17/00—Monitoring; Testing ; Maintaining
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Plasma & Fusion (AREA)
- High Energy & Nuclear Physics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明涉及核电站事故应急的技术领域,提出一种核事故辐射数据的模拟方法、装置、终端设备和计算机存储介质。该方法包括:在模拟核事故发生后,确定核反应堆的辐射量;在核辐射从核反应堆向预设的多个目标区域扩散的过程中,按照核扩散路径的先后顺序依次计算各个目标区域的辐射量,所述核扩散路径由核反应堆和各个目标区域组成,且核反应堆处于所述核扩散路径的首端;各个目标区域的辐射量根据前级区域的辐射量、沉积于前级区域的辐射量以及泄露至后级区域的辐射量计算得到。与传统的依靠经验值设置各个区域的辐射量相比,本方法能够更准确地得到各个区域的模拟辐射量,提高核事故演习的效果。
Description
技术领域
本发明涉及核电站事故应急的技术领域,尤其涉及一种核事故辐射数据的模拟方法、装置、终端设备和计算机存储介质。
背景技术
核事故是指大型核设施(例如核燃料生产厂、核反应堆、核电厂、核动力舰船及后处理厂等)发生的意外事件,可能造成人员受到放射损伤和放射性污染等严重后果。
核事故的演习对于提高核电站相关人员的事故应急处理能力至关重要,在进行核事故演习时通常需要模拟设置各种场景下的辐射数据。目前人们主要依靠经验值设置这些辐射数据,然而仅仅依靠经验值设置的辐射数据往往不够准确,造成演习效果差。
发明内容
有鉴于此,本发明实施例提供了一种核事故辐射数据的模拟方法、装置、终端设备和计算机存储介质,旨在解决核事故演习中辐射数据模拟不准确,演习效果差的问题。
本发明实施例的第一方面提供了一种核事故辐射数据的模拟方法,包括:
在模拟核事故发生后,确定核反应堆的辐射量;
在核辐射从所述核反应堆向预设的多个目标区域扩散的过程中,按照核扩散路径的先后顺序依次计算各个所述目标区域的辐射量,所述核扩散路径由所述核反应堆和各个所述目标区域组成,且所述核反应堆处于所述核扩散路径的首端;
其中,所述目标区域的辐射量通过以下步骤计算:
根据前级区域的辐射量、沉积于所述前级区域的辐射量以及泄露至后级区域的辐射量计算所述目标区域的辐射量,所述前级区域为所述核扩散路径中处于所述目标区域前一级的区域,所述后级区域为所述核扩散路径中处于所述目标区域后一级的区域,沉积于所述核扩散路径中的每个区域的辐射量均为预设的已知量,所述泄露至后级区域的辐射量与所述目标区域的辐射量之间呈预设的比例关系。
本发明实施例的第二方面提供了一种核事故辐射数据的模拟装置,包括:
反应堆辐射量确定模块,用于在模拟核事故发生后,确定核反应堆的辐射量;
目标区域辐射量计算模块,用于在核辐射从所述核反应堆向预设的多个目标区域扩散的过程中,按照核扩散路径的先后顺序依次计算各个所述目标区域的辐射量,所述核扩散路径由所述核反应堆和各个所述目标区域组成,且所述核反应堆处于所述核扩散路径的首端;
其中,所述目标区域的辐射量通过以下步骤计算:
根据前级区域的辐射量、沉积于所述前级区域的辐射量以及泄露至后级区域的辐射量计算所述目标区域的辐射量,所述前级区域为所述核扩散路径中处于所述目标区域前一级的区域,所述后级区域为所述核扩散路径中处于所述目标区域后一级的区域,沉积于所述核扩散路径中的每个区域的辐射量均为预设的已知量,所述泄露至后级区域的辐射量与所述目标区域的辐射量之间呈预设的比例关系。
本发明实施例的第三方面提供了一种终端设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如本发明实施例的第一方面提供的核事故辐射数据的模拟方法的步骤。
本发明实施例的第四方面提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现如本发明实施例的第一方面提供的核事故辐射数据的模拟方法的步骤。
在本发明实施例中,在模拟核事故发生后,确定核反应堆的辐射量;在核辐射从所述核反应堆向预设的多个目标区域扩散的过程中,按照核扩散路径的先后顺序依次计算各个所述目标区域的辐射量,所述核扩散路径由所述核反应堆和各个所述目标区域组成,且所述核反应堆处于所述核扩散路径的首端;各个目标区域的辐射量根据前级区域的辐射量、沉积于所述前级区域的辐射量以及泄露至后级区域的辐射量计算得到。与传统的依靠经验值设置各个区域的辐射量相比,本方法通过预置核反应堆的辐射量,在核辐射的扩散过程中按照核扩散路径的先后顺序依次计算各个所述目标区域的辐射量,能够更准确地得到各个区域的模拟辐射量,提高核事故演习的效果。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种核事故辐射数据的模拟方法的第一个实施例的流程图;
图2是本发明实施例提供的一种核事故辐射数据的模拟方法的第二个实施例的流程图;
图3是本发明实施例提供的一种核事故辐射数据的模拟方法的第三个实施例的流程图;
图4是本发明实施例提供的一种核事故辐射数据的模拟装置的一个实施例的结构图;
图5是本发明实施例提供的一种终端设备的示意图。
具体实施方式
以下描述中,为了说明而不是为了限定,提出了诸如特定***结构、技术之类的具体细节,以便透彻理解本发明实施例。然而,本领域的技术人员应当清楚,在没有这些具体细节的其它实施例中也可以实现本发明。在其它情况中,省略对众所周知的***、装置、电路以及方法的详细说明,以免不必要的细节妨碍本发明的描述。
本发明实施例提供了一种核事故辐射数据的模拟方法、装置、终端设备和计算机存储介质,旨在解决核事故演习中辐射数据模拟不准确,演习效果差的问题。
请参阅图1,本发明实施例中一种核事故辐射数据的模拟方法的第一个实施例包括:
101、在模拟核事故发生后,确定核反应堆的辐射量;
核反应堆的辐射量为核扩散过程的初始源项,可以根据核反应堆的能量大小预先设置。
可选的,为了进一步提高数据模拟的真实性,步骤101可以包括:
(1)获取核反应堆的堆芯温度;
(2)根据所述堆芯温度确定当前所处的核事故阶段;
(3)将当前所处的核事故阶段对应的初始辐射量作为所述核反应堆的辐射量。
对于上述步骤(1)至(3),获取核反应堆的堆芯温度,根据所述堆芯温度确定当前所处的核事故阶段。可以将核事故辐射的迁移过程划分为多个核事故阶段,比如初始阶段、破损释放气体阶段和芯块熔化阶段,各阶段之间跳转时其源项呈数量级的增长。当堆芯温度小于某个阈值(比如650℃)时,判定当前处于初始阶段,将初始阶段对应的初始辐射量作为所述核反应堆的辐射量;当堆芯温度持续上升超过某个阈值时,判定进入破损释放气体阶段,将破损释放气体阶段对应的初始辐射量作为所述核反应堆的辐射量;当堆芯温度进一步上升超过某个阈值(比如2800℃)时,假设无任何注入手段挽救堆芯,则判定进入芯块熔化阶段,将芯块熔化阶段对应的初始辐射量作为所述核反应堆的辐射量。对于每个阶段来说,可以用该阶段最大的源项作为初始辐射量,以包络各种事故的可能。比如,在初始阶段可以用所有燃料棒破损0.25%产生的辐射量作为对应的初始辐射量,在破损释放气体阶段可以用所有燃料棒破损100%产生的辐射量作为对应的初始辐射量,在芯块熔化阶段可以用所有燃料棒100%熔化产生的辐射量作为对应的初始辐射量。
102、在核辐射从所述核反应堆向预设的多个目标区域扩散的过程中,按照核扩散路径的先后顺序依次计算各个所述目标区域的辐射量。
在确定核反应堆的辐射量之后,在核辐射从所述核反应堆向预设的多个目标区域扩散的过程中,按照核扩散路径的先后顺序依次计算各个所述目标区域的辐射量。所述核扩散路径由所述核反应堆和各个所述目标区域组成,且所述核反应堆处于所述核扩散路径的首端。比如某个核扩散路径为:核反应堆—目标区域A—目标区域B—目标区域C,在核反应堆的辐射量已知的情况下,首先计算目标区域A的辐射量,然后计算目标区域B的辐射量,最后计算目标区域C的辐射量。
某个目标区域的辐射量通过以下步骤计算:
根据前级区域的辐射量、沉积于所述前级区域的辐射量以及泄露至后级区域的辐射量计算所述目标区域的辐射量。
所述前级区域为所述核扩散路径中处于所述目标区域前一级的区域,所述后级区域为所述核扩散路径中处于所述目标区域后一级的区域,沉积于所述核扩散路径中的每个区域的辐射量均为预设的已知量,所述泄露至后级区域的辐射量与所述目标区域的辐射量之间呈预设的比例关系。同样以核扩散路径:核反应堆—目标区域A—目标区域B—目标区域C作为例子,在计算得到目标区域A的辐射量之后,目标区域B的辐射量根据前级区域的辐射量(即目标区域A的辐射量)、沉积于所述前级区域的辐射量(即沉积于目标区域A的辐射量)以及泄露至后级区域的辐射量(即泄露至目标区域C的辐射量)计算。其中目标区域A的辐射量和沉积于目标区域A的辐射量都为已知量,而泄露至目标区域C的辐射量与目标区域B的辐射量之间呈预设的比例关系(比如60%),因此可以通过建立方程式计算得到目标区域B的辐射量以及泄露至目标区域C的辐射量。在计算得到目标区域B的辐射量之后,可以采用相同的方法计算目标区域C的辐射量。
进一步的,在步骤102之后,还可以包括:
(1)确定各个所述目标区域中的关键测量点;
(2)获取各个关键测量点的辐射参数,所述辐射参数为用于计算辐射浓度的预设参数;
(3)根据计算得到的各个所述目标区域的辐射量,结合各个所述关键测量点的辐射参数,分别计算得到各个所述关键测量点的辐射浓度。
对于上述步骤(1),确定各个所述目标区域中的关键测量点,关键测量点可以是预设的任意已知测量点,比如模拟的核电站场景中各个虚拟测量仪器(比如辐射监测仪表KRT)所处的位置。
对于上述步骤(2),获取各个关键测量点的辐射参数,所述辐射参数为用于计算辐射浓度的预设参数,比如可以是测量点的温度、压力、封闭区域的容积、压差以及墙壁厚度等各类与辐射浓度计算相关的参数。在某个已知的虚拟场景中,这些辐射参数均为预设的已知量。
对于上述步骤(3),在计算得到的各个所述目标区域的辐射量之后,结合各个所述关键测量点的辐射参数,分别计算得到各个所述关键测量点的辐射浓度。在得到辐射浓度之后,可以通过虚拟的仪表将辐射浓度显示出来,使得相关人员可以直观地看到各个关键测量点的辐射浓度。
在本发明实施例中,在模拟核事故发生后,确定核反应堆的辐射量;在核辐射从所述核反应堆向预设的多个目标区域扩散的过程中,按照核扩散路径的先后顺序依次计算各个所述目标区域的辐射量,所述核扩散路径由所述核反应堆和各个所述目标区域组成,且所述核反应堆处于所述核扩散路径的首端;各个目标区域的辐射量根据前级区域的辐射量、沉积于所述前级区域的辐射量以及泄露至后级区域的辐射量计算得到。与传统的依靠经验值设置各个区域的辐射量相比,本方法通过预置核反应堆的辐射量,在核辐射的扩散过程中按照核扩散路径的先后顺序依次计算各个所述目标区域的辐射量,能够更准确地得到各个区域的模拟辐射量,提高核事故演习的效果。
请参阅图2,本发明实施例中一种核事故辐射数据的模拟方法的第二个实施例包括:
201、在模拟核事故发生后,确定核反应堆的辐射量;
步骤201与步骤101相同,具体可参照步骤101的相关说明。
202、在核辐射从所述核反应堆向预设的多个目标区域扩散的过程中,按照核扩散路径的先后顺序依次计算核反应堆的一回路、反应堆厂房、辅助厂房以及大气区域的辐射量;
在本实施例中,所述目标区域包括核反应堆的一回路、反应堆厂房、辅助厂房以及大气区域,所述反应堆厂房为放置所述核反应堆以及所述一回路的封闭式区域,所述辅助厂房为与所述反应堆厂房相邻的封闭式区域。核扩散路径可分为两路,一路为:核反应堆—一回路—反应堆厂房—辅助厂房—大气区域,另一路为:核反应堆—一回路—反应堆厂房—大气区域。通过这样设置,更贴近于实际的核电站场景,模拟效果更好。
具体的,某个目标区域的辐射量可以通过以下公式(2-1)计算:
C(t+dt)=C(t)+[S(t)-D(t)-F(t)]*dt(2-1)
其中C(t+dt)表示当前时刻该目标区域的辐射量,C(t)表示上一时刻该目标区域的辐射量,S(t)表示该目标区域的前级区域的辐射量,D(t)表示沉积于该目标区域的前级区域的辐射量,F(t)表示泄露至该目标区域的后级区域的辐射量,t表示时间。C(t+dt)和F(t)之间呈预设的比例关系,即F(t)=C(t+dt)*a,a为已知常数。显然,以上各个辐射量均为时间的函数,随时间的改变而改变。
进一步的,还可以对核裂变产物的族进行区分,进而将公式(2-1)修改为以下的公式(2-2)。
C(t+dt)=C(t)+[S(t)-D(t)-F(t)]*dt(2-2)
iiiii
其中,下标i表示裂变产物的族。若i=x,则公式(2-2)计算的是x族裂变产物的辐射量,若i=y,则公式(2-2)计算的是y族裂变产物的辐射量,以此类推。
根据该核扩散路径,在核反应堆的辐射量已知的情况下,首先计算一回路的辐射量。具体的,所述一回路的辐射量通过以下步骤计算:
将所述核反应堆的辐射量减去沉积于所述核反应堆的辐射量,得到从所述核反应堆扩散至所述一回路的辐射总量;
将从所述核反应堆扩散至所述一回路的辐射总量减去从所述一回路泄露至所述反应堆厂房的辐射量,得到所述一回路的辐射量,从所述一回路泄露至所述反应堆厂房的辐射量等于所述一回路的辐射量乘以预设比例系数。
假设核反应堆的辐射量为S1(t),沉积于核反应堆的辐射量为D1(t),一回路当前时刻的辐射量为C2(t+dt),一回路上一时刻的辐射量为C2(t),从一回路泄露至反应堆厂房的辐射量为F2(t)=C2(t+dt)*a(a为预设的常数,比如60%),则C2(t+dt)=C2(t)+[S1(t)-D1(t)-F2(t)],对于初始时刻t=0来说,C2(0)=0,在S1(t)和D1(t)已知的情况下即可计算得到C2(t+dt)。需要注意的是,这里的下标1和2用于区分各个目标区域,比如1代表核反应堆,2代表一回路,与公式(2-2)中的下标i完全不同。
在计算得到一回路的辐射量之后,开始计算核扩散路径的下一个目标区域,即反应堆厂房的辐射量。具体的,反应堆厂房的辐射量通过以下步骤计算:
将所述一回路的辐射量减去沉积于所述一回路的辐射量,得到从所述一回路扩散至所述反应堆厂房的辐射总量;
将从所述一回路扩散至所述反应堆厂房的辐射总量减去从所述反应堆厂房泄露至所述辅助厂房和所述大气区域的辐射量,得到所述反应堆厂房的辐射量,从所述反应堆厂房泄露至所述辅助厂房和所述大气区域的辐射量等于所述反应堆厂房的辐射量乘以预设比例系数。
假设一回路的辐射量为C2(t+dt),沉积于一回路的辐射量为D2(t),反应堆厂房当前时刻的辐射量为C3(t+dt),反应堆厂房上一时刻的辐射量为C3(t),从反应堆厂房泄露至辅助厂房和大气区域的辐射量为F3(t)=C3(t+dt)*b(b为预设的常数,可以通过墙壁厚度等参数计算得到,比如0.3%),则C3(t+dt)=C3(t)+[C2(t+dt)-D2(t)-F3(t)],对于初始时刻t=0来说,C3(0)=0,在C2(t+dt)和D2(t)已知的情况下即可计算得到C3(t+dt)。
在计算得到反应堆厂房的辐射量后,开始计算辅助厂房的辐射量。具体的,辅助厂房的辐射量通过以下步骤计算:
将所述反应堆厂房的辐射量减去沉积于所述反应堆厂房的辐射量,得到从所述反应堆厂房扩散至所述辅助厂房和所述大气区域的辐射总量;
将从所述反应堆厂房扩散至所述辅助厂房和所述大气区域的辐射总量乘以第一比例系数,得到从所述反应堆厂房扩散至所述辅助厂房的辐射总量;
将从所述反应堆厂房扩散至所述辅助厂房的辐射总量减去从所述辅助厂房泄露至所述大气区域的辐射量,得到所述辅助厂房的辐射量,从所述辅助厂房泄露至所述大气区域的辐射量等于所述辅助厂房的辐射量乘以预设比例系数。
假设反应堆厂房的辐射量为C3(t+dt),沉积于反应堆厂房的辐射量为D3(t),辅助厂房当前时刻的辐射量为C4(t+dt),辅助厂房上一时刻的辐射量为C4(t),从辅助厂房泄露至大气区域的辐射量为F4(t)=C4(t+dt)*c(c为预设的常数,比如60%),第一比例系数为x1(小于1的正数,比如97%),则C4(t+dt)=C4(t)+[(C3(t+dt)-D3(t))*x1-F4(t)],对于初始时刻t=0来说,C4(0)=0,在C3(t+dt)和D3(t)已知的情况下即可计算得到C4(t+dt)。
在计算得到辅助厂房的辐射量后,开始计算大气区域的辐射量,大气区域的前级区域有两个,分别是反应堆厂房和辅助厂房,因此扩散至所述大气区域的辐射总量等于从所述辅助厂房扩散至所述大气区域的辐射总量和从所述反应堆厂房扩散至所述大气区域的辐射总量之和。具体的,大气区域的辐射量通过以下步骤计算:
将所述辅助厂房的辐射量减去沉积于所述辅助厂房的辐射量,得到从所述辅助厂房扩散至所述大气区域的辐射总量;
将从所述反应堆厂房扩散至所述辅助厂房和所述大气区域的辐射总量乘以第二比例系数,得到从所述反应堆厂房扩散至所述大气区域的辐射总量,所述第一比例系数和所述第二比例系数均为小于1的正数,且所述第一比例系数和所述第二比例系数的和为1;
将从所述辅助厂房扩散至所述大气区域的辐射总量和从所述反应堆厂房扩散至所述大气区域的辐射总量相加,得到扩散至所述大气区域的辐射总量;
将扩散至所述大气区域的辐射总量减去从所述大气区域泄露至所述大气区域的后级区域的辐射量,得到所述大气区域的辐射量,从所述大气区域泄露至所述大气区域的后级区域的辐射量等于所述大气区域的辐射量乘以预设比例系数。
假设反应堆厂房的辐射量为C3(t+dt),沉积于反应堆厂房的辐射量为D3(t),辅助厂房的辐射量为C4(t+dt),沉积于辅助厂房的辐射量为D4(t),大气区域当前时刻的辐射量为C5(t+dt),大气区域上一时刻的辐射量为C5(t),从大气区域泄露至后级区域的辐射量为F5(t)=C5(t+dt)*d(d为预设的常数,比如60%),第二比例系数为x2(小于1的正数,且和第一比例系数的和为1,即若x1为97%,则x2为3%),则C5(t+dt)=C5(t)+[(C3(t+dt)-D3(t))*x2+C4(t+dt)-D4(t)-F5(t)],对于初始时刻t=0来说,C5(0)=0,在C3(t+dt)、D3(t)、C4(t+dt)和D4(t)均已知的情况下即可计算得到C5(t+dt)。
进一步的,在计算得到核反应堆的一回路、反应堆厂房、辅助厂房以及大气区域的辐射量之后,还可以采用相同的方法计算大气区域的后级区域(比如核电站外部某个建筑物)的辐射量。
203、获取核反应堆的一回路、反应堆厂房、辅助厂房以及大气区域的关键测量点的辐射参数,结合各区域的辐射量和辐射参数计算这些关键测量点的实时辐射浓度。
具体的,步骤203可以包括但不限于:
(1)根据反应堆厂房的辐射量以及机组状态、事故序列、芯块、包壳、RCP破口、温度、压力、容积等辐射参数计算反应堆厂房内各个指定测量点的辐射浓度(比如厂房底部测量点、厂房安全壳高度测量点的辐射浓度)。
(2)根据核反应堆的辐射量、滞留比例、以及一回路中的各类管道参数(如管径、材料、弯曲角度等)计算一回路管道的实时辐射浓度(比如水放射性活度及变化等),得到实时辐射浓度后可通过设置于对应位置的模拟KRT设备显示实时的数值。
(3)根据反应堆厂房的放射性浓度,在正常泄漏情况下,根据安全壳打压实验中反应堆厂房与辅助厂房贯穿件的比例等比设置泄漏比例,根据反应堆厂房与辅助厂房、大气区域的泄漏截面积和压差确定泄漏比率,从而计算辅助厂房中各个指定测量点的实时辐射浓度。
(4)根据通风情况、外部区域某建筑的墙壁厚度以及正负压差计算该建筑内的实时辐射浓度。
(5)根据反应堆厂房正常泄漏或异常泄漏情况下其迁移路径上各位置的放射性浓度计算应急抢修位置的气体放射性浓度和剂量率的数值及变化。
(6)根据反应堆厂房正常泄漏或异常泄漏情况下迁移到大气中的放射性浓度及气象条件模拟厂区主要通道及重点区域(比如UA/UD、EU、主控室、辐射防护值班室等)的气体放射性浓度和剂量率的数值和变化。
(7)根据反应堆厂房正常泄漏或异常泄漏情况下迁移到大气中的放射性浓度及气象条件模拟电厂固定伽马辐射监测站的监测数据。
(8)根据反应堆厂房正常泄漏或异常泄漏情况下迁移到大气中的放射性浓度及气象条件模拟核电站基地内及周边任意指定地点的监测数据。
在本发明实施例中,在模拟核事故发生后,确定核反应堆的辐射量;在核辐射从所述核反应堆向预设的多个目标区域扩散的过程中,按照核扩散路径的先后顺序依次计算核反应堆的一回路、反应堆厂房、辅助厂房以及大气区域的辐射量;获取核反应堆的一回路、反应堆厂房、辅助厂房以及大气区域的关键测量点的辐射参数,结合各区域的辐射量和辐射参数计算这些关键测量点的实时辐射浓度。本实施例通过核扩散路径的先后顺序依次计算一回路、反应堆厂房、辅助厂房以及大气区域的辐射量,然后结合各区域的辐射量和辐射参数计算各个关键测量点的实时辐射浓度,能够更准确地得到各个区域的模拟辐射量以及各个关键测量点的实时辐射浓度,提高核事故演习的效果。
请参阅图3,本发明实施例中一种核事故辐射数据的模拟方法的第三个实施例包括:
301、在模拟核事故发生后,确定核反应堆的辐射量;
步骤301与步骤101相同,具体可参照步骤101的相关说明。
302、在核辐射从所述核反应堆向预设的多个目标区域扩散的过程中,按照核扩散路径的先后顺序依次计算核反应堆的一回路、核反应堆的二回路以及大气区域的辐射量;
在本实施例中,所述目标区域包括核反应堆的一回路、核反应堆的二回路以及大气区域,核扩散路径为:核反应堆—一回路—二回路—大气区域。在进行计算时,同样采用公式(2-1)。
所述一回路的辐射量通过以下步骤计算:
将所述核反应堆的辐射量减去沉积于所述核反应堆的辐射量,得到从所述核反应堆扩散至所述一回路的辐射总量;
将从所述核反应堆扩散至所述一回路的辐射总量减去从所述一回路泄露至所述二回路的辐射量,得到所述一回路的辐射量,从所述一回路泄露至所述二回路的辐射量等于所述一回路的辐射量乘以预设比例系数。
假设核反应堆的辐射量为S1(t),沉积于核反应堆的辐射量为D1(t),一回路当前时刻的辐射量为C2(t+dt),一回路上一时刻的辐射量为C2(t),从一回路泄露至反应堆厂房的辐射量为F2(t)=C2(t+dt)*a(a为预设的常数,比如40%),则C2(t+dt)=C2(t)+[S1(t)-D1(t)-F2(t)],对于初始时刻t=0来说,C2(0)=0,在S1(t)和D1(t)已知的情况下即可计算得到C2(t+dt)。
在计算得到一回路的辐射量之后,开始计算二回路的辐射量,二回路的辐射量通过以下步骤计算:
将所述一回路的辐射量减去沉积于所述一回路的辐射量,得到从所述一回路扩散至所述二回路的辐射总量;
将从所述一回路扩散至所述二回路的辐射总量减去从所述二回路泄露至所述大气区域的辐射量,得到所述二回路的辐射量,从所述二回路泄露至所述大气区域的辐射量等于所述二回路的辐射量乘以预设比例系数。
假设一回路当前时刻的辐射量为C2(t+dt),沉积于一回路的辐射量为D2(t),二回路当前时刻的辐射量为C3(t+dt),二回路上一时刻的辐射量为C3(t),从二回路泄露至大气区域的辐射量为F3(t)=C2(t+dt)*b(b为预设的常数,比如60%),则C3(t+dt)=C3(t)+[C2(t+dt)-D2(t)-F3(t)],对于初始时刻t=0来说,C3(0)=0,在C2(t+dt)和D2(t)已知的情况下即可计算得到C3(t+dt)。
在计算得到二回路的辐射量之后,可以计算大气区域的辐射量,大气区域的辐射量通过以下步骤计算:
将所述二回路的辐射量减去沉积于所述二回路的辐射量,得到从所述二回路扩散至所述大气区域的辐射总量;
将从所述二回路扩散至所述大气区域的辐射总量减去从所述大气区域泄露至所述大气区域的后级区域的辐射量,得到所述大气区域的辐射量,从所述大气区域泄露至所述大气区域的后级区域的辐射量等于所述大气区域的辐射量乘以预设比例系数。
假设二回路当前时刻的辐射量为C3(t+dt),沉积于二回路的辐射量为D3(t),大气区域当前时刻的辐射量为C4(t+dt),大气区域上一时刻的辐射量为C4(t),从大气区域泄露至后级区域的辐射量为F4(t)=C3(t+dt)*c(c为预设的常数,比如30%),则C4(t+dt)=C4(t)+[C3(t+dt)-D3(t)-F4(t)],对于初始时刻t=0来说,C4(0)=0,在C3(t+dt)和D3(t)已知的情况下即可计算得到C4(t+dt)。
303、获取核反应堆的一回路、二回路以及大气区域的关键测量点的辐射参数,结合各区域的辐射量和辐射参数计算这些关键测量点的实时辐射浓度。
在计算得到核反应堆的一回路、二回路以及大气区域的辐射量之后,可以获取各区域的关键测量点的辐射参数,并结合各区域的辐射量和辐射参数计算这些关键测量点的实时辐射浓度。
具体的,可以包括但不限于:
(1)根据核反应堆的辐射量、滞留比例、以及一回路中的各类管道参数(如管径、材料、弯曲角度等)计算一回路管道的实时辐射浓度(比如水放射性活度及变化等),得到实时辐射浓度后可通过设置于对应位置的模拟KRT设备显示实时的数值。
(2)根据一回路的辐射量、N16的浓度、SG等效破口管子的数量、一二回路之间的压差、二回路蒸汽管道参数(如管径、材料、弯曲角度等)计算二回路管道的实时辐射浓度(包括辐射液浓度以及辐射气体浓度),得到实时辐射浓度后可通过设置于对应位置的模拟KRT设备显示实时的数值。
在本发明实施例中,在模拟核事故发生后,确定核反应堆的辐射量;在核辐射从所述核反应堆向预设的多个目标区域扩散的过程中,按照核扩散路径的先后顺序依次计算核反应堆的一回路、核反应堆的二回路以及大气区域的辐射量;获取核反应堆的一回路、二回路以及大气区域的关键测量点的辐射参数,结合各区域的辐射量和辐射参数计算这些关键测量点的实时辐射浓度。本实施例通过核扩散路径的先后顺序依次计算核反应堆的一回路、二回路以及大气区域的辐射量,然后结合各区域的辐射量和辐射参数计算各个关键测量点的实时辐射浓度,能够更准确地得到各个区域的模拟辐射量以及各个关键测量点的实时辐射浓度,提高核事故演习的效果。
应理解,上述各个实施例中各步骤的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本发明实施例的实施过程构成任何限定。
上面主要描述了一种核事故辐射数据的模拟方法,下面将对一种核事故辐射数据的模拟装置进行描述。
请参阅图4,本发明实施例中一种核事故辐射数据的模拟装置的一个实施例包括:
反应堆辐射量确定模块401,用于在模拟核事故发生后,确定核反应堆的辐射量;
目标区域辐射量计算模块402,用于在核辐射从所述核反应堆向预设的多个目标区域扩散的过程中,按照核扩散路径的先后顺序依次计算各个所述目标区域的辐射量,所述核扩散路径由所述核反应堆和各个所述目标区域组成,且所述核反应堆处于所述核扩散路径的首端;
其中,所述目标区域的辐射量通过以下步骤计算:
根据前级区域的辐射量、沉积于所述前级区域的辐射量以及泄露至后级区域的辐射量计算所述目标区域的辐射量,所述前级区域为所述核扩散路径中处于所述目标区域前一级的区域,所述后级区域为所述核扩散路径中处于所述目标区域后一级的区域,沉积于所述核扩散路径中的每个区域的辐射量均为预设的已知量,所述泄露至后级区域的辐射量与所述目标区域的辐射量之间呈预设的比例关系。
进一步的,所述反应堆辐射量确定模块401可以包括:
堆芯温度获取单元,用于获取核反应堆的堆芯温度;
核事故阶段确定单元,用于根据所述堆芯温度确定当前所处的核事故阶段;
反应堆辐射量确定单元,用于将当前所处的核事故阶段对应的初始辐射量作为所述核反应堆的辐射量。
进一步的,所述核事故辐射数据的模拟装置还可以包括:
关键测量点确定模块,用于确定各个所述目标区域中的关键测量点;
辐射参数获取模块,用于获取各个关键测量点的辐射参数,所述辐射参数为用于计算辐射浓度的预设参数;
辐射浓度计算模块,用于根据计算得到的各个所述目标区域的辐射量,结合各个所述关键测量点的辐射参数,分别计算得到各个所述关键测量点的辐射浓度。
本发明实施例还提供一种终端设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如图1至图3表示的任意一种核事故辐射数据的模拟方法的步骤。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现如图1至图3表示的任意一种核事故辐射数据的模拟方法的步骤。
图5是本发明一实施例提供的终端设备的示意图。如图5所示,该实施例的终端设备5包括:处理器50、存储器51以及存储在所述存储器51中并可在所述处理器50上运行的计算机程序52。所述处理器50执行所述计算机程序52时实现上述各个核事故辐射数据的模拟方法的实施例中的步骤,例如图1所示的步骤101至102。或者,所述处理器50执行所述计算机程序52时实现上述各装置实施例中各模块/单元的功能,例如图4所示模块401至402的功能。
所述计算机程序52可以被分割成一个或多个模块/单元,所述一个或者多个模块/单元被存储在所述存储器51中,并由所述处理器50执行,以完成本发明。所述一个或多个模块/单元可以是能够完成特定功能的一系列计算机程序指令段,该指令段用于描述所述计算机程序52在所述终端设备5中的执行过程。
所述终端设备5可以是各种类型的手机、桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。所述终端设备可包括,但不仅限于,处理器50、存储器51。本领域技术人员可以理解,图5仅仅是终端设备5的示例,并不构成对终端设备5的限定,可以包括比图示更多或更少的部件,或者组合某些部件,或者不同的部件,例如所述终端设备5还可以包括输入输出设备、网络接入设备、总线等。
所称处理器50可以是中央处理单元(Central Processing Unit,CPU),还可以是其他通用处理器、数字信号处理器(Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现成可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
所述存储器51可以是所述终端设备5的内部存储单元,例如终端设备5的硬盘或内存。所述存储器51也可以是所述终端设备5的外部存储设备,例如所述终端设备5上配备的插接式硬盘,智能存储卡(Smart Media Card,SMC),安全数字(Secure Digital,SD)卡,闪存卡(Flash Card)等。进一步地,所述存储器51还可以既包括所述终端设备5的内部存储单元也包括外部存储设备。所述存储器51用于存储所述计算机程序以及所述终端设备所需的其他程序和数据。所述存储器51还可以用于暂时地存储已经输出或者将要输出的数据。
所属领域的技术人员可以清楚地了解到,为了描述的方便和简洁,仅以上述各功能单元、模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能单元、模块完成,即将所述装置的内部结构划分成不同的功能单元或模块,以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中,上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。另外,各功能单元、模块的具体名称也只是为了便于相互区分,并不用于限制本申请的保护范围。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的***,装置和单元的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述或记载的部分,可以参见其它实施例的相关描述。
本领域普通技术人员可以意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
在本发明所提供的实施例中,应该理解到,所揭露的装置和方法,可以通过其它的方式实现。例如,以上所描述的***实施例仅仅是示意性的,例如,所述模块或单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个***,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通讯连接可以是通过一些接口,装置或单元的间接耦合或通讯连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明实现上述实施例方法中的全部或部分流程,也可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一计算机可读存储介质中,该计算机程序在被处理器执行时,可实现上述各个方法实施例的步骤。其中,所述计算机程序包括计算机程序代码,所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括:能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,RandomAccess Memory)、电载波信号、电信信号以及软件分发介质等。需要说明的是,所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减,例如在某些司法管辖区,根据立法和专利实践,计算机可读介质不包括是电载波信号和电信信号。
以上所述实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围,均应包含在本发明的保护范围之内。
Claims (8)
1.一种核事故辐射数据的模拟方法,其特征在于,包括:
在模拟核事故发生后,确定核反应堆的辐射量;
在核辐射从所述核反应堆向预设的多个目标区域扩散的过程中,按照核扩散路径的先后顺序依次计算各个所述目标区域的辐射量,所述核扩散路径由所述核反应堆和各个所述目标区域组成,且所述核反应堆处于所述核扩散路径的首端;
其中,所述目标区域包括核反应堆的一回路、反应堆厂房、辅助厂房以及大气区域,所述反应堆厂房为放置所述核反应堆以及所述一回路的封闭式区域,所述辅助厂房为与所述反应堆厂房相邻的封闭式区域,所述核扩散路径包括两路,其中一路依次为:核反应堆、核反应堆的一回路、反应堆厂房、辅助厂房、大气区域,另一路依次为:核反应堆、核反应堆的一回路、反应堆厂房、大气区域;
所述一回路的辐射量通过以下步骤计算:
将所述核反应堆的辐射量减去沉积于所述核反应堆的辐射量,得到从所述核反应堆扩散至所述一回路的辐射总量;
将从所述核反应堆扩散至所述一回路的辐射总量减去从所述一回路泄露至所述反应堆厂房的辐射量,得到所述一回路的辐射量,从所述一回路泄露至所述反应堆厂房的辐射量等于所述一回路的辐射量乘以预设比例系数;
所述反应堆厂房的辐射量通过以下步骤计算:
将所述一回路的辐射量减去沉积于所述一回路的辐射量,得到从所述一回路扩散至所述反应堆厂房的辐射总量;
将从所述一回路扩散至所述反应堆厂房的辐射总量减去从所述反应堆厂房泄露至所述辅助厂房和所述大气区域的辐射量,得到所述反应堆厂房的辐射量,从所述反应堆厂房泄露至所述辅助厂房和所述大气区域的辐射量等于所述反应堆厂房的辐射量乘以预设比例系数;
所述辅助厂房的辐射量通过以下步骤计算:
将所述反应堆厂房的辐射量减去沉积于所述反应堆厂房的辐射量,得到从所述反应堆厂房扩散至所述辅助厂房和所述大气区域的辐射总量;
将从所述反应堆厂房扩散至所述辅助厂房和所述大气区域的辐射总量乘以第一比例系数,得到从所述反应堆厂房扩散至所述辅助厂房的辐射总量;
将从所述反应堆厂房扩散至所述辅助厂房的辐射总量减去从所述辅助厂房泄露至所述大气区域的辐射量,得到所述辅助厂房的辐射量,从所述辅助厂房泄露至所述大气区域的辐射量等于所述辅助厂房的辐射量乘以预设比例系数;
所述大气区域的辐射量通过以下步骤计算:
将所述辅助厂房的辐射量减去沉积于所述辅助厂房的辐射量,得到从所述辅助厂房扩散至所述大气区域的辐射总量;
将从所述反应堆厂房扩散至所述辅助厂房和所述大气区域的辐射总量乘以第二比例系数,得到从所述反应堆厂房扩散至所述大气区域的辐射总量,所述第一比例系数和所述第二比例系数均为小于1的正数,且所述第一比例系数和所述第二比例系数的和为1;
将从所述辅助厂房扩散至所述大气区域的辐射总量和从所述反应堆厂房扩散至所述大气区域的辐射总量相加,得到扩散至所述大气区域的辐射总量;
将扩散至所述大气区域的辐射总量减去从所述大气区域泄露至所述大气区域的后级区域的辐射量,得到所述大气区域的辐射量,从所述大气区域泄露至所述大气区域的后级区域的辐射量等于所述大气区域的辐射量乘以预设比例系数;
或者
其中,所述目标区域包括核反应堆的一回路、核反应堆的二回路以及大气区域,所述核扩散路径依次为:核反应堆、核反应堆的一回路、核反应堆的二回路、大气区域;
所述一回路的辐射量通过以下步骤计算:
将所述核反应堆的辐射量减去沉积于所述核反应堆的辐射量,得到从所述核反应堆扩散至所述一回路的辐射总量;
将从所述核反应堆扩散至所述一回路的辐射总量减去从所述一回路泄露至所述二回路的辐射量,得到所述一回路的辐射量,从所述一回路泄露至所述二回路的辐射量等于所述一回路的辐射量乘以预设比例系数;
所述二回路的辐射量通过以下步骤计算:
将所述一回路的辐射量减去沉积于所述一回路的辐射量,得到从所述一回路扩散至所述二回路的辐射总量;
将从所述一回路扩散至所述二回路的辐射总量减去从所述二回路泄露至所述大气区域的辐射量,得到所述二回路的辐射量,从所述二回路泄露至所述大气区域的辐射量等于所述二回路的辐射量乘以预设比例系数;
所述大气区域的辐射量通过以下步骤计算:
将所述二回路的辐射量减去沉积于所述二回路的辐射量,得到从所述二回路扩散至所述大气区域的辐射总量;
将从所述二回路扩散至所述大气区域的辐射总量减去从所述大气区域泄露至所述大气区域的后级区域的辐射量,得到所述大气区域的辐射量,从所述大气区域泄露至所述大气区域的后级区域的辐射量等于所述大气区域的辐射量乘以预设比例系数。
2.如权利要求1所述的核事故辐射数据的模拟方法,其特征在于,所述确定核反应堆的辐射量包括:
获取核反应堆的堆芯温度;
根据所述堆芯温度确定当前所处的核事故阶段;
将当前所处的核事故阶段对应的初始辐射量作为所述核反应堆的辐射量。
3.如权利要求1或2所述的核事故辐射数据的模拟方法,其特征在于,在按照核扩散路径的先后顺序依次计算各个所述目标区域的辐射量之后,还包括:
确定各个所述目标区域中的关键测量点;
获取各个关键测量点的辐射参数,所述辐射参数为用于计算辐射浓度的预设参数;
根据计算得到的各个所述目标区域的辐射量,结合各个所述关键测量点的辐射参数,分别计算得到各个所述关键测量点的辐射浓度。
4.一种核事故辐射数据的模拟装置,其特征在于,包括:
反应堆辐射量确定模块,用于在模拟核事故发生后,确定核反应堆的辐射量;
目标区域辐射量计算模块,用于在核辐射从所述核反应堆向预设的多个目标区域扩散的过程中,按照核扩散路径的先后顺序依次计算各个所述目标区域的辐射量,所述核扩散路径由所述核反应堆和各个所述目标区域组成,且所述核反应堆处于所述核扩散路径的首端;
其中,所述目标区域包括核反应堆的一回路、反应堆厂房、辅助厂房以及大气区域,所述反应堆厂房为放置所述核反应堆以及所述一回路的封闭式区域,所述辅助厂房为与所述反应堆厂房相邻的封闭式区域,所述核扩散路径包括两路,其中一路依次为:核反应堆、核反应堆的一回路、反应堆厂房、辅助厂房、大气区域,另一路依次为:核反应堆、核反应堆的一回路、反应堆厂房、大气区域;
所述目标区域辐射量计算模块包括:
一回路辐射量计算单元,用于将所述核反应堆的辐射量减去沉积于所述核反应堆的辐射量,得到从所述核反应堆扩散至所述一回路的辐射总量;将从所述核反应堆扩散至所述一回路的辐射总量减去从所述一回路泄露至所述反应堆厂房的辐射量,得到所述一回路的辐射量,从所述一回路泄露至所述反应堆厂房的辐射量等于所述一回路的辐射量乘以预设比例系数;
反应堆厂房辐射量计算单元,用于将所述一回路的辐射量减去沉积于所述一回路的辐射量,得到从所述一回路扩散至所述反应堆厂房的辐射总量;将从所述一回路扩散至所述反应堆厂房的辐射总量减去从所述反应堆厂房泄露至所述辅助厂房和所述大气区域的辐射量,得到所述反应堆厂房的辐射量,从所述反应堆厂房泄露至所述辅助厂房和所述大气区域的辐射量等于所述反应堆厂房的辐射量乘以预设比例系数;
辅助厂房辐射量计算单元,用于将所述反应堆厂房的辐射量减去沉积于所述反应堆厂房的辐射量,得到从所述反应堆厂房扩散至所述辅助厂房和所述大气区域的辐射总量;将从所述反应堆厂房扩散至所述辅助厂房和所述大气区域的辐射总量乘以第一比例系数,得到从所述反应堆厂房扩散至所述辅助厂房的辐射总量;将从所述反应堆厂房扩散至所述辅助厂房的辐射总量减去从所述辅助厂房泄露至所述大气区域的辐射量,得到所述辅助厂房的辐射量,从所述辅助厂房泄露至所述大气区域的辐射量等于所述辅助厂房的辐射量乘以预设比例系数;
大气区域辐射量计算单元,用于将所述辅助厂房的辐射量减去沉积于所述辅助厂房的辐射量,得到从所述辅助厂房扩散至所述大气区域的辐射总量;将从所述反应堆厂房扩散至所述辅助厂房和所述大气区域的辐射总量乘以第二比例系数,得到从所述反应堆厂房扩散至所述大气区域的辐射总量,所述第一比例系数和所述第二比例系数均为小于1的正数,且所述第一比例系数和所述第二比例系数的和为1;将从所述辅助厂房扩散至所述大气区域的辐射总量和从所述反应堆厂房扩散至所述大气区域的辐射总量相加,得到扩散至所述大气区域的辐射总量;将扩散至所述大气区域的辐射总量减去从所述大气区域泄露至所述大气区域的后级区域的辐射量,得到所述大气区域的辐射量,从所述大气区域泄露至所述大气区域的后级区域的辐射量等于所述大气区域的辐射量乘以预设比例系数;
或者
其中,所述目标区域包括核反应堆的一回路、核反应堆的二回路以及大气区域,所述核扩散路径依次为:核反应堆、核反应堆的一回路、核反应堆的二回路、大气区域;
所述目标区域辐射量计算模块包括:
一回路辐射量计算单元,用于将所述核反应堆的辐射量减去沉积于所述核反应堆的辐射量,得到从所述核反应堆扩散至所述一回路的辐射总量;将从所述核反应堆扩散至所述一回路的辐射总量减去从所述一回路泄露至所述二回路的辐射量,得到所述一回路的辐射量,从所述一回路泄露至所述二回路的辐射量等于所述一回路的辐射量乘以预设比例系数;
二回路辐射量计算单元,用于将所述一回路的辐射量减去沉积于所述一回路的辐射量,得到从所述一回路扩散至所述二回路的辐射总量;将从所述一回路扩散至所述二回路的辐射总量减去从所述二回路泄露至所述大气区域的辐射量,得到所述二回路的辐射量,从所述二回路泄露至所述大气区域的辐射量等于所述二回路的辐射量乘以预设比例系数;
大气区域辐射量计算单元,用于将所述二回路的辐射量减去沉积于所述二回路的辐射量,得到从所述二回路扩散至所述大气区域的辐射总量;将从所述二回路扩散至所述大气区域的辐射总量减去从所述大气区域泄露至所述大气区域的后级区域的辐射量,得到所述大气区域的辐射量,从所述大气区域泄露至所述大气区域的后级区域的辐射量等于所述大气区域的辐射量乘以预设比例系数。
5.如权利要求4所述的核事故辐射数据的模拟装置,其特征在于,所述反应堆辐射量确定模块包括:
堆芯温度获取单元,用于获取核反应堆的堆芯温度;
核事故阶段确定单元,用于根据所述堆芯温度确定当前所处的核事故阶段;
反应堆辐射量确定单元,用于将当前所处的核事故阶段对应的初始辐射量作为所述核反应堆的辐射量。
6.如权利要求4或5所述的核事故辐射数据的模拟装置,其特征在于,还包括:
关键测量点确定模块,用于确定各个所述目标区域中的关键测量点;
辐射参数获取模块,用于获取各个关键测量点的辐射参数,所述辐射参数为用于计算辐射浓度的预设参数;
辐射浓度计算模块,用于根据计算得到的各个所述目标区域的辐射量,结合各个所述关键测量点的辐射参数,分别计算得到各个所述关键测量点的辐射浓度。
7.一种终端设备,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如权利要求1至3中任一项所述的核事故辐射数据的模拟方法的步骤。
8.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至3中任一项所述的核事故辐射数据的模拟方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711248992.XA CN108268694B (zh) | 2017-12-01 | 2017-12-01 | 核事故辐射数据的模拟方法、装置、终端设备和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711248992.XA CN108268694B (zh) | 2017-12-01 | 2017-12-01 | 核事故辐射数据的模拟方法、装置、终端设备和存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108268694A CN108268694A (zh) | 2018-07-10 |
CN108268694B true CN108268694B (zh) | 2021-03-26 |
Family
ID=62771821
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711248992.XA Active CN108268694B (zh) | 2017-12-01 | 2017-12-01 | 核事故辐射数据的模拟方法、装置、终端设备和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108268694B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110458461B (zh) * | 2019-08-14 | 2021-11-09 | 苏州热工研究院有限公司 | 核动力厂多堆事故放射性后果评价方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288174A (zh) * | 2011-06-22 | 2011-12-21 | 中核能源科技有限公司 | 一种基于核岛辐射剂量作业路径的确定方法及*** |
CN105629288A (zh) * | 2015-12-28 | 2016-06-01 | 深圳先进技术研究院 | 核电站辐射分布的预测*** |
CN106548023A (zh) * | 2016-10-27 | 2017-03-29 | 长江勘测规划设计研究有限责任公司 | 一种地下核电站严重事故气载放射性源项精确评估方法 |
CN107194103A (zh) * | 2017-06-07 | 2017-09-22 | 哈尔滨工程大学 | 一种基于3dsMax的核设施模型辐射场剂量仿真方法 |
CN107330186A (zh) * | 2017-06-30 | 2017-11-07 | 哈尔滨工程大学 | 一种3dsMax核设施模型辐射场剂量分布仿真方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9569583B2 (en) * | 2014-04-07 | 2017-02-14 | TallannQuest LLC | Method and system for computer-aided design of radiation-hardened integrated circuits |
-
2017
- 2017-12-01 CN CN201711248992.XA patent/CN108268694B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102288174A (zh) * | 2011-06-22 | 2011-12-21 | 中核能源科技有限公司 | 一种基于核岛辐射剂量作业路径的确定方法及*** |
CN105629288A (zh) * | 2015-12-28 | 2016-06-01 | 深圳先进技术研究院 | 核电站辐射分布的预测*** |
CN106548023A (zh) * | 2016-10-27 | 2017-03-29 | 长江勘测规划设计研究有限责任公司 | 一种地下核电站严重事故气载放射性源项精确评估方法 |
CN107194103A (zh) * | 2017-06-07 | 2017-09-22 | 哈尔滨工程大学 | 一种基于3dsMax的核设施模型辐射场剂量仿真方法 |
CN107330186A (zh) * | 2017-06-30 | 2017-11-07 | 哈尔滨工程大学 | 一种3dsMax核设施模型辐射场剂量分布仿真方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108268694A (zh) | 2018-07-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sanchez-Saez et al. | Uncertainty analysis of a large break loss of coolant accident in a pressurized water reactor using non-parametric methods | |
Di Maio et al. | A dynamic probabilistic safety margin characterization approach in support of Integrated Deterministic and Probabilistic Safety Analysis | |
Arkoma et al. | Statistical analysis of fuel failures in large break loss-of-coolant accident (LBLOCA) in EPR type nuclear power plant | |
Ahn et al. | FONESYS: the FOrum & NEtwork of SYStem thermal-hydraulic codes in nuclear reactor thermal-hydraulics | |
CN108268694B (zh) | 核事故辐射数据的模拟方法、装置、终端设备和存储介质 | |
CN106844208A (zh) | 用于反应堆物理计算程序适用性验证的方法 | |
Alcaro et al. | BEPU analysis of a passive decay heat removal system with RELAP5/MOD3. 3 and RELAP5-3D | |
Freixa et al. | On the validation of BEPU methodologies through the simulation of integral experiments: Application to the PKL test facility | |
Xu et al. | ATHLET simulation of PKL I2. 2 IB-loca benchmark test and quantitative assessment | |
Mochizuki | Evaluation of spent fuel pool temperature and water level during SBO | |
CN111157193B (zh) | 核电站密封性试验方法、密封性试验装置及电子设备 | |
CN104898155A (zh) | 一种基于中子平均能量的次临界***次临界度测量方法 | |
CN101834001A (zh) | 针对核事故损失及应急行动代价的优化控制方法 | |
Park et al. | Smart sensing of the RPV water level in NPP severe accidents using a GMDH algorithm | |
Ivanov et al. | Best estimate plus uncertainty (BEPU): why it is still not widely used | |
CN112488465B (zh) | 一种八氧化三铀临界事故应急时总裂变次数估算方法 | |
Zio et al. | Bootstrap and Order Statistics for Quantifying Thermal‐Hydraulic Code Uncertainties in the Estimation of Safety Margins | |
Park et al. | Numerical study on the thermal-hydraulic behavior in the ultimate heat sink of passive residual heat removal system in the SMART | |
CN112488461B (zh) | 一种金属铀临界事故应急时总裂变次数估算方法 | |
Suescún-Díaz et al. | Novel approach to solving the inverse equation of point kinetics by the Bernoulli number generalisation method | |
CN105373665B (zh) | 获取用于核电站辐射仿真***的多核素等效参数的方法 | |
CN112819053A (zh) | 模型库的建立方法及装置、诊断方法及装置、预测方法 | |
CN105468831A (zh) | 一种反应堆工程仿真机辐射值仿真方法 | |
Yum et al. | Uncertainty assessment for loss of flow accident of a 5 MW pool-type research reactor | |
Wang et al. | Preliminary reliability allocation method for fusion facility based on probabilistic safety assessment |
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 |