[发明专利]一种核设施事故场外后果评价中的风场诊断方法

说明书

技术领域

本发明涉及核设施事故后果分析技术领域，具体涉及一种核设施事故场外后果评价中数据预处理方法。

背景技术

核设施或其他核设施在严重偏离正常运行工况、发生事故时，放射性物质的释放可能或已经失去应有的控制，达到不可接受的水平，将会造成严重的后果。因此，对核设施事故危害的防范控制及分析，具有重要的战略意义。

核设施事故后果评价系统是通过将区域气象数据和放射性污染物迁移模式和后果评价集成为一体的评价系统，用于估计、评价和显示在特定范围内的事故后果，该系统是核设施或其他核设施应急准备和应急防护行动决策的重要组成部分。核设施事故后果评价系统，首先需要把类别与内容各异的源项、气象、环境信息等数据进行预处理，能够对上述各数据实现有效的处理对后续事故后果的分析结果有着至关重要的影响，本发明正是针对该问题而提出了一种核设施事故场外后果评价中的风场诊断方法。

发明内容

针对现有技术中存在的缺陷，本发明的目的在于提供一种能够构造出局地气象场的核设施事故场外后果评价中的风场诊断方法。

为实现上述目的，本发明采用的技术方案如下：

一种核设施事故场外后果评价中的风场诊断方法，包括以下步骤：

步骤一、以核设施为中心，确定核设施事故的模拟区范围，获取模拟区的气象数据和地理环境数据，并根据获取的气象数据和地理位置数据建立模拟区的三维初始风场；所述气象数据包括模拟区的气象观测站的观测气象数据和气象部门的数值天气预报数据；所述地理环境数据包括模拟区的经纬度信息、下垫面高程、植被、水体分布特征和土地利用；三维初始风场为多层网格风场；

步骤二、根据所述气象部门的数值天气预报数据和地理环境数据对所述三维三维初始风场进行初步调整，得到初步调整后的风场；所述初步调整包括根据地形运动学效应对三维初始风场三维初始风场进行调整、根据坡流效应对三维初始风场三维初始风场进行调整和根据地形热动力学对风场的阻塞效应对三维初始风场三维初始风场进行调整；

步骤三、根据所述观测气象数据对初步调整后的风场进行处理，得到最终风场；所述处理包括对初步调整后的风场依次进行插值处理、平滑处理和垂直风速分量的调整。

进一步，如上所述的一种核设施事故场外后果评价中的风场诊断方法，步骤二中，所述根据地形运动学效应对三维初始风场三维初始风场进行调整包括：根据地形运动学效应对垂直方向风速分量进行调整和根据地形运动学效应对水平方向风速分量进行调整；

根据地形运动学效应对垂直方向风速分量进行调整的方式为：

2.1)计算笛卡尔坐标系中风速的垂直分量w，计算公式为：

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

2.2)计算地形跟踪坐标系下风速的垂直分量W，计算公式为：

其中，u、v分别为笛卡尔坐标系下风速在x轴水平方向的分量和在y轴水平方向的分量，为地形高度在x轴方向的梯度，为地形高度在y轴方向的梯度；x轴水平方向为东西走向，y轴水平方向为南北走向；

根据地形运动学效应对水平方向风场分量进行修正的方式为：采用偏差最小化方法对水平风场分量采用迭代法进行调整，直至调整后的三维偏差小于设定值。

进一步，如上所述的一种核设施事故场外后果评价中的风场诊断方法，步骤二中，根据坡流效应对三维初始风场进行调整的方式为：

1)计算坡流效应产生的坡流风速；坡流风速S的计算公式为：

S＝{[Q_hgx sinα/[(ρc_pT)(C_D+k)]}^1/3[1-exp(-x/L_e)]^1/3

Le＝h/(C_D+k)

其中，Q_h为显热通量，g为重力加速度，x为坡流距离坡顶的距离，α为坡流水平角度，ρ为空气密度，c_p为空气热容，T为空气温度，C_D为地面阻力系数，k为坡流层顶部的夹带系数，L_e为平衡长度尺度，h为坡流高度；