1.一种核设施事故场外后果实时在线评价方法，该方法首先以核设施为中心确定核设施事故的模拟区范围，所述模拟区范围包括以核设施为中心20km×20km的近场区域和以核设施为中心100km×100km的远场区域；
在确定出模拟区范围后，该在线评价方法包括模拟区的风场预测方法、风场诊断方法、拉格朗日烟团大气扩散模拟方法、气态放射性核素所致辐射剂量与防护措施建议的估算方法、以及气载途径食物链放射性活度估算方法；
其中，所述风场预测方法包括：
1)以核设施为中心，获取模拟区范围的不同来源的气象数据；所述气象数据为全球或中尺度数值天气预报数据或者为待预测中尺度区域的观测站测得的观测数据；
2)根据所述气象数据来源的不同，进行模拟区的风场预测，进行风场预测的具体方式为：
21)当所述气象数据为全球或中尺度数值天气预报数据时，采用非静力模式进行风场预测；
22)当所述气象数据为待预测中尺度区域的观测站测得的观测数据时，采用准静力模式进行风场预测；
所述风场诊断方法包括：
(1)获取模拟区的气象数据和地理环境数据，并根据获取的气象数据和地理位置数据建立模拟区的三维风场；所述气象数据包括模拟区的气象观测站的观测气象数据和气象部门的数值天气预报数据；所述地理环境数据包括模拟区的经纬度信息、下垫面高程、植被、水体分布特征和土地利用；三维风场为多层网格风场；
(2)根据所述气象部门的数值天气预报数据和地理环境数据对所述三维初始风场进行初步调整，得到初步调整后的风场；所述初步调整包括根据地形运动学效应对三维初始风场进行调整、根据坡流效应对三维初始风场进行调整和根据地形热动力学对风场的阻塞效应对三维初始风场进行调整；
(3)根据所述观测气象数据对初步调整后的风场进行处理，得到最终的三维风场；所述处理包括对初步调整后的风场依次进行插值处理、平滑处理和垂直风速分量的调整；
所述拉格朗日烟团大气扩散模拟方法通过依次释放一系列烟团模拟放射性核素的连续释放来计算气载放射性核素扩散产生的浓度，包括以下步骤：
A、建立模拟区范围的二维网格系统；
B、计算二维网格系统中每个网格点的各核素的浓度，计算方式为：
B1.确定计算的时间步长ΔT即浓度结果输出时间间隔和烟团的释放时间间隔Δt，根据所述烟团的释放时间间隔Δt顺序释放一系列烟团模拟放射性核素的连续释放；
B2.计算第M个时间步长的第W个时段释放的第i个烟团中的j核素对网格点(xg,yg,zg)空气中的j核素的浓度贡献 计算公式为：


其中，(xg,yg,zg)中(xg,yg)为网格点的二维坐标，zg为烟团中核素取对人体有影响的高度， 为为第M个时间步长的第W个时段释放的第i个烟团中j核素的源强，QWjo第M个时间步长所释放的所有烟团中j核素的源强，σxy(i)、σz(i)分别为所述第i个烟团水平方向和垂直方向的有效扩散参数， 为所述第i个烟团在第M个时间步长结束时刻的
质心坐标，zinv为逆温层层顶的垂直高度，λj为j核素的衰变常数， 为第M个时间步长结束时刻所述第i个烟团所经的迁移时间；
B3.计算从事故发生时刻起至第M个时间步长结束时刻网格点(xg,yg,zg)空气中的j核素的时间积分浓度χj(xg,yg,zg；M)，计算公式为：

其中，n表示每个时刻烟团的释放数；
C、计算核素的地表沉积量，计算方式为：
C1.计算从事故发生时刻起至第M个时间步长结束时刻因干沉积造成的(x,y)网格中j核素的地表干沉积总量WDj(x,y；M)，计算公式为：

其中，Vdj为j核素的干沉积速度，χj(x,y；M)为χj(x,y,z；M)在z＝0至z＝∞高度的积分；
C2.计算第M1个时间步长到第M2个时间步长期间的降雨造成的(x,y)网格中j核素的地表湿沉积总量Wwj(x,y,M1→M2)，计算公式为：

其中，∧j为冲洗因子，I为降雨强度，A为冲洗因子系数，A的取值范围为[3×10-5,3×10-3]，a的取值范围为[0.5,1]；
D、计算一个烟团中的所有核素的γ射线对指定受照点(x,y,z)造成的γ辐射剂量率dγ(Q,Eγ,σy,σz,H,Rxy)，计算公式为：


其中，K＝1.6×10-13，单位为Gy/s/MeV/Kg，σen为空气的能量吸收系数，单位为m2/Kg，Eγ-1
为γ射线的辐射能量，单位为MeV，B(μr)为累积因子，μ为空气的线性减弱因子，单位为m ，r为从离烟团中心(x＝y＝0,z＝-H)距离为Rxy的受照点到体积元dxdydz的距离，H为烟团中心高度，χ(x,y,z)为在指定受照点(x,y,z)处的瞬时空气浓度，(x,y,z)中z为指定受照点的高度，(x,y)为指定受照点多对应的网格坐标，Q为一个烟团中的放射性核素的活度，σxy、σz分别烟团水平方向和垂直方向的扩散参数；
所述气态放射性核素所致辐射剂量与防护措施建议的估算方法包括：
a、根据预设的网格间距对模拟区进行网格划分；
b、计算核设施事故释放期间，气态放射性核素在环境中迁移扩散对模拟区居民及居民生活环境所造成的影响，包括：
b1.对于每一个网格，计算不同照射途径下气态放射性核素对网格中居民所造成潜在剂量，并根据潜在剂量计算不同照射途径下气态放射性核素对网格中居民所造成的预期剂量；所述照射途径包括烟云γ外照射途径、地面沉积外照射途径和吸入内照射途径；
b2.对于每一个网格，计算气态放射性核素释放开始后两天放射性核素对网格中的居民所造成的预期剂量，根据计算结果判断是否有确定性效应发生的区域；根据计算结果判断是否有确定性效应发生的区域的方式为：
当气态放射性核素释放开始后两天内放射性核素对网格中的居民的各器官所造成的与其剂量是否大于预设的各器官的确定性效应阈值时，则判断有确定性效应发生。
b3.对于每一个网格，计算采取各种应急干预行动后，网格中居民的可防止剂量和剩余剂量；所述应急干预行动包括隐蔽、撤离、避迁、服用碘片；在采取撤离与避迁行动时，计算撤离与避迁过程中各撤离路线和避迁路线引起的附加剂量；
b4.计算事故释放期间，牛奶中的关键核素浓度和水库中的关键核素浓度，根据计算结果判定是否达到了禁止牛奶或地表水饮用的干预水平；所述关键核素包括34Cs、137Cs、103Ru、106Ru、89Sr、131I、90Sr；
所述气载途径食物链放射性活度估算方法包括：
I、对于模拟区，计算核设施事故期间所释放的污染物烟羽经过时所造成的核素沉积量；对于某一区域，核设施事故期间所释放的污染物烟羽经过时所造成的核素沉积量包括污染物烟羽经过该区域上方所造成的植物上的核素的干沉积量和湿沉积量、以及土壤表面核素的总沉积量；
II、根据所述核素沉积量，计算植物收割时植物可食部分中核素的浓度以及动物产品中核素的浓度；
III、根据植物收割时植物可食部分中核素的浓度以及动物产品中核素的浓度，计算人体通过食入途径对核素的摄入速率。
2.根据权利要求1所述的一种核设施事故场外后果实时在线评价方法，其特征在于：所述风场预测方法的步骤21)中，根据动力学方程、热力学方程、水汽方程、连续性方程进行非静力模式下的风场预测，其中，所述动力学方程、热力学方程、水汽方程、连续性方程如下：
动力学方程如下：




P0＝1000hPa


其中，x、y、z为三维笛卡尔坐标系下的东西向、南北向和垂直向的坐标，为地形追随坐标系中的垂直坐标， 为 坐标系中的三维速度分量，t为时间，θ为位
温，π为扰动艾克纳函数，f为科氏参数，Km为湍流扩散系数，g为重力加速度，θv为虚温，cp为干空气比定压热容，P为大气压，P0＝1000hPa为参考气压，R为干空气比气体常数，H为设定的顶部边界高度，zg为地形高度，z为笛卡尔坐标系下的高度；
热力学方程如下：

其中，kh为热量和湿气的湍流粘性系数， 表示以弧度表示
水汽方程如下：
其中，rn表示水物质混合比；
连续方程如下：

其中，ρ为流体密度。
步骤22)中，根据两个水平动力学方程、准静力方程、连续性方程、热力学方程和湍流能量方程进行静力模式下的风场预测，其中：两个水平动力学方程、准静力方程、连续性方程、热力学方程和湍流能量方程如下：
动力学方程如下：

准静力方程如下：


其中，x、y、z为三维笛卡尔坐标系下的东西向、南北向和垂直向的坐标，为地形追随坐标系中的垂直坐标， 为 坐标系中的三维速度分量，t为时间，zg为
地形高度，H为设定的顶部边界高度，θ为位温，π为扰动艾克纳函数，g为重力加速度，Fu和Fv分别为东西向风速的湍流扩散项和南北向风速的湍流扩散项；
热力学方程如下：

其中，Fθ为位温的湍流扩散项；
连续性方程如下：


其中，π为表示大气压的扰动艾克纳Exner函数，cp为干空气比定压热容，P为大气压，P0＝1000hPa为参考气压，R为干空气的比气体常数；
湍流能量方程如下：


其中，q2为湍流动能， 分别为x、y、z三个方向上的单位质量湍流动
能，Kw表示垂直动量湍流扩散系数，α为比容，Kθ为温度的热量湍流扩散系数，l为湍流尺度，B1l取经验值10.46，Fq2为湍流动能的湍流扩散项；
上述各公式中， 其中，为u、v、θ或者q2，KH、
KZ为对应湍流扩散项的水平扩散系数和垂直扩散系数。
其中，分别采用非静力或静力模式进行风场预测时，首先进行远场区域的风场预测，完成远场区域的风场预测后，以远场区域的风场预测结果作为输入数据获取近场区域的风场。
3.根据权利要求1所述的一种核设施事故场外后果实时在线评价方法，其特征在于：所述风场诊断方法的步骤(2)中，根据地形运动学效应对三维初始风场进行调整包括：根据地形运动学效应对垂直方向风速分量进行调整和根据地形运动学效应对水平方向风速分量进行调整；根据地形运动学效应对垂直方向风速分量进行调整包括以下步骤：
(2-1)、计算笛卡尔坐标系中风速的垂直分量w，计算公式为：



其中，V为模拟区域平均风速，ht为地形高度， 为ht的梯度，表示随地高度差，k为与大气稳定度有关的指数衰减系数，z为笛卡尔坐标系下的垂直坐标，N为布伦特-维萨拉频率，|V|为区域平均风速的绝对值，g为重力加速度，θ为环境位温；
(2-2)、计算地形跟踪坐标系下风速的垂直分量W，计算公式为：

其中，u、v分别为笛卡尔坐标系下风速在x轴水平方向的分量和在y轴水平方向的分量，为地形高度在x轴方向的梯度， 为地形高度在y轴方向的梯度；x轴水平方向为东西走向，y轴水平方向为南北走向；
根据地形运动学效应对水平方向风场分量进行修正的方式为：采用偏差最小化方法对水平风场分量采用迭代法进行调整，直至调整后的三维偏差小于设定值；
根据坡流效应对三维初始风场进行调整的方式包括以下步骤：
(2.1)计算坡流效应产生的坡流风速；坡流风速S的计算公式为：
S＝{[Qh g x sinα/[(ρcp T)(CD+k)]}1/3[1-exp(-x/Le)]1/3
Le＝h/(CD+k)
其中，Qh为显热通量，g为重力加速度，x为坡流距离坡顶的距离，α为坡流水平角度，ρ为空气密度，cp为空气热容，T为空气温度，CD为地面阻力系数，k为坡流层顶部的夹带系数，Le为平衡长度尺度，h为坡流高度；
(2.2)根据坡流风速和坡流水平角度α计算出x轴水平方向的坡流风速us和y轴水平方向的坡流风速vs；
当坡流为下坡流时，
CD＝K＝4×10-2
h＝0.05ΔZ
sinα＝minimum(sinα，ΔZ/x)
当坡流为上坡流时，(CD+k)～1，x sinα＝ΔZ；ΔZ为坡流高度；
(2.3)将坡流风速加入到三维初始风场，对三维初始风场的水平风速进行调整，调整的公式为：
u1′＝u1+us
v1′＝v1+vs
其中，u1′和v1′分别为根据坡流效应对三维初始风场进行修正后的x轴水平方向的风速分量和y轴水平方向的风速分量，u1和v1分别为三维初始风场中x轴水平方向的风速分量和y轴水平方向的风速分量。
根据地形热动力学对风场的阻塞效应对三维初始风场进行调整的方式为：
(21)计算三维初始风场中每个网格点的局地弗洛德数Fr，计算公式为：

Δht＝(hmax)ij-(z)ijk
其中，V为网格点的风速，N为布伦特-维萨拉频率，Δht为有效障碍高度，(hmax)ij为网格点(i,j)影响半径内的最大地形高度，(z)ijk为高空层k中网格点(i,j)的高度；
(22)判断网格点局地弗洛德数Fr是否小于设定的临界弗洛德数，若否，则不对三维初始风场进行调整，若是，则判断网格点的风场是否有向上的分量，若是，则将网格点的风向调整为与地形相切的方向，若否，则不对三维初始风场进行调整；
步骤(3)中，根据观测气象数据对初步调整后的风场进行插值处理的插值公式为：

其中，(u,v)′2为插值处理后的风场中网格点的水平方向的初始风速分量，(u,v)1为所述初步调整后的风场中网格点的水平方向的风速分量，k为气象观测站的标识序号，(uobs,vobs)k为第k个气象观测站的水平方向的观测风速分量，R为预设的初步调整后的风场的权重数；Rk为第k个气象观测站到网格点的距离；
根据观测气象数据对插值处理后风场进行平滑处理，平滑处理公式为：
(ui,j)2″＝0.5ui,j+0.125(ui-1,j+ui+1,j+ui,j-1+ui,j+1)
(vi,j)2″＝0.5vi,j+0.125(vi-1,j+vi+1,j+vi,j-1+vi,j+1)
其，(ui,j)2″、(vi,j)2″分别为平滑处理后在网格点(i,j)处的水平方向中x轴方向的风速分量和y轴方向的风速分量；ui,j、vi,j分别为平滑处理前插值处理后在网格点(i,j)处的水平方向中x轴方向的风速分量和y轴方向的风速分量；
根据观测气象数据对平滑处理后的风场进行垂直风速调整的方式包括：根据不压缩质量守恒方程从平滑处理后的水平风场分量中计算出垂直风速分量，计算出的垂直风速分量为最终的垂直风速分量，平滑处理后的水平风场为最终的水平风场，或者，
在根据不压缩质量守恒方程从平滑处理后的水平风场分量中计算出垂直风速分量后，对计算出的垂直风速进行调整，调整后的垂直风速分量为最终的垂直风速分量；
其中，根据不压缩质量守恒方程从平滑处理后的水平风场分量中计算出垂直风速分量W1(z)的公式为：

其中，u″、v″分别为平滑处理后的水平方向中x轴方向的风速分量和y轴方向的风速分量；
在根据不压缩质量守恒方程从平滑处理后的水平风场分量中计算出垂直风速分量后，对计算出的垂直风速进行调整，得到调整后的垂直风速分量W2(z)的计算公式为：
W2(z)＝W1(z)-(z/ztop)W1(z＝ztop)
其中，z为笛卡尔坐标系下的垂直坐标，ztop为模拟区的顶部高度,W1(z＝ztop)为根据不压缩质量守恒方程从平滑处理后的水平风场分量中计算出的模拟区顶部的垂直速度分量；
计算出调整后的垂直风速分量W2(z)后，采用偏差最小化方法调整两个水平方向的风速分量，偏差最小化方法在保持垂直分量为W2(z)不变的情况下，利用迭代方法将两个水平方向的风速分量的偏差调整到小于设定的偏差阈值，采用偏差化方法调整后的两个水平方向的水平分量为最终的水平风场分量；
计算出调整后的垂直风速分量W2(z)后，采用偏差最小化方法调整两个水平方向的风速分量的调整的方式为：
计算网格单元中每个网格点(i,j,k)处的水平风向的风速分量偏差值Dijk，计算公式为：

其中，zk+1/2-zk-1/2表示垂直方向上的层高，Δx和Δy分别表示在x轴和y轴方向上网格单元的边长；
调整网格点(i,j,k)的周围网格单元的水平风速分量，使网格点的偏差值Dijk小于设定的偏差阈值ε，调整公式为：
(unew)i+1,j,k＝ui+1,j,k+uadj
(unew)i-1,j,k＝ui-1,j,k-uadj
(vnew)i+1,j,k＝vi,j+1,k+vadj
(vnew)i,j-1,k＝vi,j-1,k-vadj


其中，uadj和vadj分别为在x轴和y轴方向上风速分量的调整量；
迭代使用上述调整方式，直至所有网格单元的偏差值均小于设定的偏差阈值ε。
4.根据权利要求1所述的一种核设施事故场外后果实时在线评价方法，其特征在于：所述拉格朗日烟团大气扩散模拟方法的步骤B2中，在计算第i个烟团中的j核素对网格点(xg,yg,zg)空气中的j核素的浓度之前，还包括判断烟团是否发生分裂的步骤，在烟团发生分裂，分别计算分裂后的每个子烟团中的j核素对网格点(xg,yg,zg)空气中的j核素的浓度贡献；
烟团分裂的具体方式为：一个分裂前半径为σp的原始烟团在水平方向上分裂为半径为σp/2的一个中心烟团和四个卫星烟团，中心烟团与原始烟团的圆心重合，四个卫星烟团沿中心烟团的圆周方向均匀分布且与原始烟团圆心的距离为0.89σp，中心烟团的质量为原始烟团的5.88％，卫星烟团的质量为原始烟团的23.5％。
5.根据权利要求1所述的一种核设施事故场外后果实时在线评价方法，其特征在于：气态放射性核素所致辐射剂量与防护措施建议的估算方法的步骤b1中，烟云γ外照射途径下气态放射性核素对居民所造成的潜在剂量的计算公式为：

其中，L表示网格，O表示人体器官，O＝1，2分别表示甲状腺、全身；NU表示气态放射性核素；M表示计算的时间步长的个数，ΔT为一个时间步长，ΔT＝0.5h，2T为计算的时间大于24小时的实时评价小时数； 为烟云γ外照射途径下所有气态放射性核素在48+2T个步长时间对网格L中居民的O器官所造成的潜在剂量；Dγ(L,NU,M)是网格L中NU核素在第M个步长时间的烟云γ剂量率；
烟云γ外照射途径下气态放射性核素对每个网格中的居民全身所造成的预期剂量的计算公式为：

其中， 为48+2T个步长时间内烟云γ外照射途径下气态放射性核素对L
网格中居民全身所造成的预期剂量，FNCL(L)me为L网格对应于烟云γ外照射途径的平均居址因子；
地面沉积外照射途径下气态放射性核素对居民所造成的潜在剂量的计算公式为：

其中，I表示地面沉积外照射累积照射时间，I＝1，2，3，4分别对应2天、7天、1个月和50年； 为所有气态放射性核素在地面沉积48+2T个步长时间内对网格L中居民
的O器官所造成的潜在剂量；WD(L,NU,M)为NU核素在第M个步长时间在网格L的地面沉积浓度；DRFGR(NU,O)是一个时间步长内NU核素对O器官的地面沉积剂量率转换因子；WD(L,NU,48+2T)为NU核素在第48+2T个步长时间在网格L的地面沉积浓度；DFGR(NU,O，I)是I时间内NU核素对O器官的地面沉积剂量率转换因子；α为地面沉积结束后，所剩的计算时段I-48-2T占总计算时段I的时间份额；
地面沉积外照射途径下气态放射性核素对每个网格中的居民全身所造成的预期剂量的计算公式为：

为48+2T个步长时间内地面沉积外照射途径下气态放射性核素对L
网格中居民全身所造成的预期剂量；FNGR(L)me为L网格对应地面沉积外照射途径的平均居址因子；
吸入内照射途径下气态放射性核素对居民所造成的潜在剂量的计算公式为：

其中， 表示吸入内照射途径下所有气态放射性核素在48+2T个步长时间
内对L网格居民的O器官所造成的潜在剂量；χ(L,NU,M)表示在第M个步长时间L网格中NU核素在空气中的浓度；DFIH(NU,O,I＝4)为居民的O器官对NU核素的吸入内照射途径剂量转换因子，RIH为居民的呼吸速率；
吸入内照射途径下气态放射性核素对每个网格中的居民全身所造成的预期剂量的计算公式为：

其中， 为48+2T个步长时间内L网格对应于吸入内照射途径的预期剂
量；FNIH(L)me(t)为当前时间t时刻L网格对应吸入内照射途径的平均居址因子；
步骤b2中，气态放射性核素释放开始后两天放射性核素对网格中的居民所造成的预期剂量的计算公式为：

其中，DSUMPOT(L,O)为气态放射性核素释放开始后两天放射性核素对L网格内居民的O器官所造成的预期剂量，
分别为释放开始后两天烟云γ外照射途径、地面沉积外照途径和吸入内照射途径下气态放射性核素对L网格内居民各器官的潜在剂量；
步骤b3中，当采取隐蔽的应急干预行动时，计算各网格内居民在事故释放开始后自第l小时起隐蔽2天的可防止剂量和剩余剂量，l为0.5h的整倍数，包括：
计算L网格居民事故释放后自第l小时起隐蔽2天的烟云γ外照射途径的可防止剂量计算公式为：

其中，FNCL(L)me为L网格对应于烟云γ外照射途径的平均居址因子，FSCL(L)me为L网格对应烟云γ外照射途径的平均建筑物屏蔽因子，当48+2T-2l>96时，取48+2T-2l＝96；
计算L网格居民事故释放后自第l小时起隐蔽2天的地面沉积外照射途径的可防止剂量计算公式为：

其中，FNGR(L)me为L网格对应地面沉积外照射途径的平均居址因子,FSGR(L)me表示L网格对应于地面沉积外照射途径的建筑物屏蔽因子，其中当48+2T-2l>96时，取48+2T-2l＝96；
计算L网格居民事故释放后自第l小时起隐蔽2天的吸入内照射途径的可防止剂量计算公式为：

其中，FNIH(L,M)me与FSIH(L,M)me分别为L网格第M个时间步长对应于吸入内照射途径的居址因子与建筑物屏蔽因子，其中当48+2T-2l>96时，取48+2T-2l＝96；
计算L网格居民事故释放后自第l小时起隐蔽2天的总的可防止剂量DAVE(L,O＝2,l→l+
2d)，计算公式为：

计算L网格居民事故释放后自第l小时起隐蔽2天的烟云γ外照射途径的剩余剂量计算公式为：

计算L网格居民事故释放后自第l小时起隐蔽2天的地面沉积外照射途径的剩余剂量计算公式为：

计算L网格居民事故释放后自第l小时起隐蔽2天的吸入内照射途径的剩余剂量计算公式为：

计算L网格居民事故释放后自第l小时起隐蔽2天的总的剩余剂量Dle(L,O＝2,l→l+
2d)，计算公式为：

当隐蔽2天的总的可防止剂量是否大于第一设定值时，则对对应网格的居民采取隐蔽措施；
当采取撤离的应急干预行动时，计算各网格内居民在事故释放开始后自第l小时起撤离7天的可防止剂量和剩余剂量，以及撤离过程中所造成的附加剂量；包括：
计算L网格居民事故释放后自第l小时起撤离7天的烟云γ外照射途径的可防止剂量计算公式为：

其中，N为L网格的撤离路线所经过的网格总数；Ln(x,y)是L网格居民撤离过程所经过的第n个网格的坐标，Δt′为撤离时间，取48+2T≤48×7+2l；
计算L网格居民事故释放后自第l小时起撤离7天的地面沉积外照射途径的可防止剂量计算公式为：

其中，当T>6×24+l时，取T＝6×24+l；
计算L网格居民事故释放后自第l小时起撤离7天的吸入内照射途径的可防止剂量计算公式为：

计算L网格居民事故释放后自第l小时起撤离7天的总的可防止剂量DAVE(L,O＝2,l→l+
7d)，计算公式为：

计算L网格居民事故释放后自第l小时起撤离7天的烟云γ外照射途径的剩余剂量计算公式为：

计算L网格居民事故释放后自第l小时起撤离7天的地面沉积外照射途径的剩余剂量计算公式为：

计算L网格居民事故释放后自第l小时起撤离7天的吸入内照射途径的剩余剂量计算公式为：

le
计算第l网格居民事故释放后自第l小时起撤离7天的总的剩余剂量D (L,O＝2,l→l+
7d)，计算公式为：

当撤离7天的总的可防止剂量是否大于第二设定值时，则对应网格的居民采取撤离措施；
计算采取撤离措施时，撤离过程中所造成的附加剂量，包括：
计算第m条撤离路线在撤离过程中烟云γ外照射途径的附加剂量 计算公式为：

其中，nm是第m条撤离路线所经过的网格总数，Lnm(x,y)是第m条撤离路线所经过的第n个网格的坐标，Δt′m是第m条撤离路线所需的撤离时间；
计算第m条撤离路线在撤离过程中地面沉积外照射途径的附加剂量 计算公式为：

计算第m条撤离路线在撤离过程中吸入内照射途径的附加剂量 计算公式为

计算第m条撤离路线在撤离过程中所受总的附加剂量 的计算公式为：

当采取复用碘片的应急干预行动时，计算各网格内居民事故释放后第l小时服碘片的甲状腺可防止剂量与剩余剂量，包括：
计算L网格居民所受的甲状腺预期剂量 计算公式为：

DFIH(NU,O＝1,I＝4)·RIH·ΔT·FNIH(L)me
计算L网格居民事故释放后第l小时服碘片的剩余剂量 计算公式
为：

其中，tIH为核素碘的吸入时间；a(M)为减少因子；
tIH-l·60分≤15分时，a(M)＝1-exp{-λ·[(l·60-tIH)+15分]}；
tIH-l·60分>15分时，a(M)＝0；
λ为时间常数，λ＝0.00289分-1；
计算L网格居民事故释放后第l小时服碘片的可防止剂量 计算
公式为：

若服碘片的可防止剂量大于第三设定值，则对对应网格的居民采取服碘片措施；
当采取避迁的应急干预行动时，计算各网格点居民事故释放开始后自第l小时起避迁的可防止剂量和剩余剂量，以及避迁过程中所造成的附加剂量，包括：
计算L网格居民事故释放后自第l小时起避迁的烟云γ外照射途径的可防止剂量计算公式为：

其中，N为L网格撤离路线所经过的网格总数，Ln(x,y)是L网格居民避迁过程所经过的第n个网格的坐标，Δt′为撤离时间，取48+2T≤48×30+2l；
计算L网格居民事故释放后自第l小时起避迁的地面沉积外照射途径的可防止剂量计算公式为：

其中，当T>29×24+l时，取T＝29×24+l；取48+2T≤48×30+2l；
计算L网格居民事故释放后第l小时起避迁的吸入内照射途径的可防止剂量
计算公式为：

其中，取48+2T≤48×30+2l；
计算L网格居民事故释放后第l小时起避迁的总的可防止剂量DAVE(L,O＝2,l→l+30d)，计算公式为：

当避迁的总的可防止剂量大于第四设定值时，则对应网格内居民采取撤离措施；
计算L网格居民事故释放后自第l小时起避迁的烟云γ外照射途径的剩余剂量
计算公式为：

计算L网格居民事故释放后自第l小时起避迁的地面沉积外照射途径的剩余剂量计算公式为：

计算L网格居民事故释放后自第l小时起避迁的吸入内照射途径的剩余剂量
计算公式为：

计算L网格居民事故释放后自第l小时起避迁的总的剩余剂量Dle(L,O＝2,l→l+30d)，计算公式为：

计算采取避迁措施时，避迁过程中所造成的附加剂量，包括：
计算第m条避迁路线在避迁过程中烟云γ外照射途径的附加剂量

其中，N为L网格避迁路线所经过的网格总数，Ln(x,y)是L网格居民避迁过程所经过的第n个网格的坐标，Δt′为避迁时间，取48+2T≤48×30+2l；
计算第m条避迁路线在避迁过程中地面沉积外照射途径的附加剂量 计算公式为：

计算第m条避迁路线在避迁过程中吸入内照射途径的附加剂量 计算公式为：

计算第m条避迁路线在避迁过程中总的附加剂量 计算公式为：

步骤B4中，计算牛奶中j核素的浓度Cmj(L)，计算公式为：


Qp＝16.1kF
其中，Fmj为动物每天摄入的放射性j核素出现在每公斤牛奶中的平均份额， 为饲料p
作物中的j核素的浓度，Q为动物每天消耗的饲料量；
WDj(x,y,M＝48+2T)表示事故停止时，坐标为(x,y)的网格的地面上j核素的干沉积浓度；Wwj(x,y,M＝48+2T)表示事故停止时由于事故期间降水造成坐标为(x,y)的网格的地面的j核素的湿沉积浓度；R为滞留即沉积放射性滞留在叶面上的份额；Tiv是易位因子，即滞留在植物外表面的放射性核素向植物食用部分的转移系数，对于饲料作物，偏保守地取R＝
0.5，Tiv＝1.0；Yva为饲料作物的单位面积产量； 是风化产生的清除速率常数；
λj为j核素的放射性衰变常数，λw为环境衰减常数，对于碘，取λw＝2.97×10-3d-1，对于其它-3 -1
粒子，取λw＝1.94×10 d ；te是农作物在生长季节受污染的时间，te取事故释放持续时间d，由事故源项给出；
F表示饲料作物的鲜重与干重之比，k表示本地饲料所占份额；
计算水库中的关键核素的浓度，包括：
假若事故释放期间，从M1至M2时段发生降水，降水强度为I，事故释放结束时水库水中j核素的浓度Cdj(x,y)的计算公式为：

Vl＝Sl·I·(M2-M1)·0.5×10-3
其中，WDj(x,y,M＝48+2T)表示水库所在位置(x,y)处事故结束时刻地面上(或水面上)因干沉积导致的j核素沉积浓度；Wwj(x,y,M1→M2)表示事故释放期间若从M1至M2时段发生降水，导致水库所在位置地面上或水面上的j核素湿沉积浓度；Sl为集雨面积，VC为水库库容，f为水库中的实际的储水体积与库容之比；Vl为集水量；
I为降水强度，M1和M2表示降水从M1时段开始，M2时段结束，因为每个时段步长为0.5h，故有因子0.5；
若牛奶中的各核素浓度值高于预置各核素的规定值，则应对牛奶的供应采取行动，或禁止饮用，或供应替代食品；若水库中各核素的浓度超过预置的关于饮水的各核素的通用行动水平，则应采取禁止饮用或其它水处理行动。
6.根据权利要求1所述的一种核设施事故场外后果实时在线评价方法，其特征在于：气载途径食物链放射性活度估算方法的步骤I中，对于所述某一区域，污染物烟羽经过该区域上方所造成的j类植物上i核素的干沉积量和湿沉积量、以及土壤表面的i核素的总沉积量的计算方式为：
污染物烟羽经过该区域上方所造成的j类植物上i核素的干沉积量Adij(te)的计算公式为：

Vdij＝Vdij,maxLAIj(te)/LAIj,max
其中，te为污染物烟羽尾部离开该区域的时刻，Vdij为i核素向j类植物的干沉积速度，xi为核设施事故发生后该区域地面空气中i核素的浓度，λi为i核素的放射性衰减常数，Δt为污染物烟羽经过该区域的时间长度；
Vdij,max为i核素向j类植物的最大沉积速度，LAIj(te)为污染物烟羽尾部离开该区域的时刻j类植物的植物叶面积指数，LAIj,max为j类植物的最大叶面积指数；植物的叶面积指数是指单位面积土壤上的植物叶子面积，最大沉积速度是指植物叶面积指数最大时的沉积速度；
污染物烟羽经过该区域上方期间由于降雨过程的湿沉积所造成的j类植物上i核素的湿沉积量Awij(te)的计算公式为：
Awij(te)＝fwjAwi(te)


R＝IkΔtk/3600
其中，fwj为j类植物的截获份额，Awi(te)污染物烟羽经过该区域上方期间由于降雨过程的湿沉积所造成的i核素的湿沉积量；
k为污染物烟羽经过该区域上方期间发生的降雨过程的次数，Λik为第k次降雨过程所对应的i核素的冲洗系数，Qi为i核素源强，u为该区域地面的平均风速，x为该区域距事故源项的下风距离，λi为i核素的放射性衰减常数，Δtk为第k次降雨过程的持续时间，te为污染物烟羽尾部离开该区域的时刻，tk为第k次降雨过程的结束时间；
LAIj为j类植物叶面积指数，Sj为j类植物有效贮水能力，R为污染物烟羽经过该区域上方期间的降雨总量，Ik为第k次降雨的降雨强度，Δtk为第k次降雨过程的持续时间；
污染物烟羽经过所述某一区域时造成的土壤表面的i核素的总沉积量Asij(te)的计算公式为：
或

其中，Vd表示i核素的沉降速度,fa表示沉积核素中未被叶面截获而到达土壤的份额；
步骤II中，对于所述某一区域，收割时j类植物可食部位中的i核素的浓度Aij(tp)的计算方式为：
仅考虑核素易位原因时，浓度Aij(tp)的计算公式为：

Tr(δt)＝Tr，max·exp[-b(δt-tmax)2]
其中，tp为j类植物的收割时间，δt为核素沉积结束至收割的时间间隔，δt＝tp-te Tr(δt)为i核素收割前δt时间发生的沉积相应的易位因子，Tr，max为易位因子最大值，Yj为j类植物的单产量，tmax为具有最大易位因子的时刻距收割的时间，b为表征实测曲线的斜率，实测曲线指的是距离收割的不同时间点测得的易位因子值与对应时间形成的曲线；
同时考虑核素易位及核素衰减作用时，浓度Aij(tp)的计算公式为：

其中，λw为环境衰减常数；
对于所述某一区域，T时刻m类动物的K类动物产品中的i核素的浓度CmKi(T)的计算公式为：


其中，FmKi表示m类动物食入含i核素的食物后转移到K类动物产品的转移系数，aKn为K类动物产品中发生n类生物转移的份额，Aami(t)表示t时刻m类动物对i核素的的摄入速度，λbkn为相应于K类动物产品发生n类生物转移时对应的转移速率常数，λi为i核素的放射性衰减常数；
AijT、AigT、AiH分别表示t时刻j类植物、牧草与干草中i核素的浓度，Ijm(t)、Igm(t)、IHm(t)分别表示m类动物t时刻对j类植物、牧草与干草三类的干重摄入率。
步骤III中，对于所述某一区域，t时刻人体通过食入途径对i核素的摄入速率的计算公式为：

其中，Aij(tp)为收割时第j类植物中i核素的浓度；Vj(t)和Vk(t)为不同年龄组居民对j类植物制品与k类动物产品的日消费量，Frj和Frk分别为j类植物与K类动物产品中的食品加工滞留因子，Pej和PeK分别为j类植物与K类动物产品的加工效率，CmKi(ts)为m类动物宰割时K类动物产品中i核素的浓度，ts为m类动物的宰割时刻。
7.一种核设施事故场外后果实时在线评价系统，包括：
模拟区范围确定模块，用于以核设施为中心确定核设施事故的模拟区范围，所述模拟区范围包括以核设施为中心20km×20km的近场区域和以核设施为中心100km×100km的远场区域；
该在线评价系统包括风场预测子系统、风场诊断子系统、大气扩散子系统、剂量与干预措施子系统、以及食物链子系统；
所述风场预测子系统包括气象数据获取模块和风场预测模块；其中：
气象数据获取模块，用于以核设施为中心，获取模拟区范围的不同来源的气象数据；所述气象数据为全球或中尺度数值天气预报数据或者为待预测中尺度区域的观测站测得的观测数据；
风场预测模块，用于根据所述气象数据来源的不同，进行模拟区的风场预测；该风场预测模块包括：
第一预测模块，用于当所述气象数据为全球或中尺度数值天气预报数据时，采用非静力模式进行风场预测；
第二预测模块，用于当所述气象数据为待预测中尺度区域的观测站测得的观测数据时，采用准静力模式进行风场预测；
所述风场诊断子系统包括三维初始风场建立模块、初始风场初步调整模块和风场再次调整模块；其中：
三维初始风场建立模块，用于获取模拟区的气象数据和地理环境数据，并根据获取的气象数据和地理位置数据建立模拟区的三维初始风场；所述气象数据包括模拟区的气象观测站的观测气象数据和气象部门的数值天气预报数据；所述地理环境数据包括模拟区的经纬度信息、下垫面高程、植被、水体分布特征和土地利用；三维风场为多层网格风场；
初始风场初步调整模块，用于根据所述气象部门的数值天气预报数据和地理环境数据对所述三维初始风场进行初步调整，得到初步调整后的风场；所述初步调整包括根据地形运动学效应对三维初始风场进行调整、根据坡流效应对三维初始风场进行调整和根据地形热动力学对风场的阻塞效应对三维初始风场进行调整；
风场再次调整模块，用于根据所述观测气象数据对初步调整后的风场进行处理，得到最终的三维风场；所述处理包括对初步调整后的风场依次进行插值处理、平滑处理和垂直风速分量的调整；
所述大气扩散子系统通过依次释放一系列烟团模拟放射性核素的连续释放来计算气载放射性核素扩散产生的浓度，该子系统包括二维网格系统建立模块、网格点核素浓度计算模块、地表沉积量计算模块和γ辐射剂量率计算模块；其中：
二维网格系统建立模块，用于建立模拟区范围的二维网格系统；
网格点核素浓度计算模块，用于计算二维网格系统中每个网格点的各核素的浓度，网格点核素浓度计算模块计算二维网格系统中每个网格点的各核素的浓度的方式为：
B1.确定计算的时间步长ΔT即浓度结果输出时间间隔和烟团的释放时间间隔Δt，根据所述烟团的释放时间间隔Δt顺序释放一系列烟团模拟放射性核素的连续释放；
B2.计算第M个时间步长的第W个时段释放的第i个烟团中的j核素对网格点(xg,yg,zg)空气中的j核素的浓度贡献 计算公式为：


其中，(xg,yg,zg)中(xg,yg)为网格点的二维坐标，zg为烟团中核素取对人体有影响的高度， 为为第M个时间步长的第W个时段释放的第i个烟团中j核素的源强，QWjo第M个时间步长所释放的所有烟团中j核素的源强，σxy(i)、σz(i)分别为所述第i个烟团水平方向和垂直方向的有效扩散参数， 为所述第i个烟团在第M个时间步长结束时刻的
质心坐标，zinv为逆温层层顶的垂直高度，λj为j核素的衰变常数， 为第M个时间步长结束时刻所述第i个烟团所经的迁移时间；
B3.计算从事故发生时刻起至第M个时间步长结束时刻网格点(xg,yg,zg)空气中的j核素的时间积分浓度χj(xg,yg,zg；M)，计算公式为：

其中，n表示每个时刻烟团的释放数；
地表沉积量计算模块，用于计算核素的地表沉积量，计算方式为：
C1.计算从事故发生时刻起至第M个时间步长结束时刻因干沉积造成的(x,y)网格中j核素的地表干沉积总量WDj(x,y；M)，计算公式为：

其中，Vdj为j核素的干沉积速度，χj(x,y；M)为χj(x,y,z；M)在z＝0至z＝∞高度的积分；
C2.计算第M1个时间步长到第M2个时间步长期间的降雨造成的(x,y)网格中j核素的地表湿沉积总量Wwj(x,y,M1→M2)，计算公式为：

其中，∧j为冲洗因子，I为降雨强度，A为冲洗因子系数，A的取值范围为[3×10-5,3×10-3]，a的取值范围为[0.5,1]；
γ辐射剂量率计算模块，用于计算一个烟团中的所有核素的γ射线对指定受照点(x,y,z)造成的γ辐射剂量率dγ(Q,Eγ,σy,σz,H,Rxy)，计算公式为：


其中，K＝1.6×10-13，单位为Gy/s/MeV/Kg，σen为空气的能量吸收系数，单位为m2/Kg，Eγ-1
为γ射线的辐射能量，单位为MeV，B(μr)为累积因子，μ为空气的线性减弱因子，单位为m ，r为从离烟团中心(x＝y＝0,z＝-H)距离为Rxy的受照点到体积元dxdydz的距离，H为烟团中心高度，χ(x,y,z)为在指定受照点(x,y,z)处的瞬时空气浓度，(x,y,z)中z为指定受照点的高度，(x,y)为指定受照点多对应的网格坐标，Q为一个烟团中的放射性核素的活度，σxy、σz分别烟团水平方向和垂直方向的扩散参数；
所述剂量与干预措施子系统包括网格划分模块和气态核素辐射剂量计算模块；其中：
网格划分模块，根据预设的网格间距对模拟区进行网格划分；
气态核素辐射剂量计算模块，用于计算核设施事故释放期间，气态放射性核素在环境中迁移扩散对模拟区居民及居民生活环境所造成的影响，包括：
b1.对于每一个网格，计算不同照射途径下气态放射性核素对网格中居民所造成潜在剂量，并根据潜在剂量计算不同照射途径下气态放射性核素对网格中居民所造成的预期剂量；所述照射途径包括烟云γ外照射途径、地面沉积外照射途径和吸入内照射途径；
b2.对于每一个网格，计算气态放射性核素释放开始后两天放射性核素对网格中的居民所造成的预期剂量，根据计算结果判断是否有确定性效应发生的区域；根据计算结果判断是否有确定性效应发生的区域的方式为：
当气态放射性核素释放开始后两天内放射性核素对网格中的居民的各器官所造成的与其剂量是否大于预设的各器官的确定性效应阈值时，则判断有确定性效应发生。
b3.对于每一个网格，计算采取各种应急干预行动后，网格中居民的可防止剂量和剩余剂量；所述应急干预行动包括隐蔽、撤离、避迁、服用碘片；在采取撤离与避迁行动时，计算撤离与避迁过程中各撤离路线和避迁路线引起的附加剂量；
b4.计算事故释放期间，牛奶中的关键核素浓度和水库中的关键核素浓度，根据计算结
34 137 103
果判定是否达到了禁止牛奶或地表水饮用的干预水平；所述关键核素包括 Cs、 Cs、 Ru、106Ru、89Sr、131I、90Sr；
所述食物链模块包括核素沉积量计算模块、动植物中核素浓度计算模块和人体核素摄入速率计算模块；其中：
核素沉积量计算模块，用于计算模拟区核设施事故期间所释放的污染物烟羽经过时所造成的核素沉积量；对于某一区域，核设施事故期间所释放的污染物烟羽经过时所造成的核素沉积量包括污染物烟羽经过该区域上方所造成的植物上的核素的干沉积量和湿沉积量、以及土壤表面核素的总沉积量；
动植物中核素浓度计算模块，用于根据所述核素沉积量，计算植物收割时植物可食部分中核素的浓度以及动物产品中核素的浓度；
人体核素摄入速率计算模块，用于根据植物收割时植物可食部分中核素的浓度以及动物产品中核素的浓度，计算人体通过食入途径对核素的摄入速率。
8.根据权利要求7所述的一种核设施事故场外后果实时在线评价系统，其特征在于：步骤B1中，在计算第i个烟团中的j核素对网格点(xg,yg,zg)空气中的j核素的浓度之前，还包括判断烟团是否发生分裂的步骤，在烟团发生分裂，分别计算分裂后的每个子烟团中的j核素对网格点(xg,yg,zg)空气中的j核素的浓度贡献；烟团分裂的具体方式为：一个分裂前半径为σp的原始烟团在水平方向上分裂为半径为σp/2的一个中心烟团和四个卫星烟团，中心烟团与原始烟团的圆心重合，四个卫星烟团沿中心烟团的圆周方向均匀分布且与原始烟团圆心的距离为0.89σp，中心烟团的质量为原始烟团的5.88％，卫星烟团的质量为原始烟团的23.5％。