Evaluation of Nuclear Accident Atmospheric Radionuclide Transport Model ARTM

早期的示踪实验主要有三种类型：（1）大气示踪剂可行性研究实验；（2）定性大气扩散实验；（3）大气扩散参数确定实验，B.Vanderborght等在1983年对这些实验进行了详细的综述^［1］，令人遗憾的是大量的实验都缺乏定量描述，而此后的示踪实验在吸取前人经验的基础上对实验条件和实验结果进行了相对详细的定量描述，使大气扩散模型的统计验证成为可能。按照实验规模，大气扩散实验分为大尺度大气示踪实验（200 km以上）和中小尺度（200 km以下）大气示踪实验。

大尺度扩散实验（200 km以上）由于其经费昂贵、实验设计和协调的复杂性而很少开展。现在能够应用于验证大尺度扩散模型的示踪实验数据集还不多，其中受到广泛关注的数据集有：CAPTEX、ANATEX、ETEX^［2］。

相对于大尺度示踪实验而言，中小尺度的示踪实验数据集数量就比较多了。Kincaid野外实验是EPRI烟羽模型验证与开发项目的一个重要试验^［3］，于1980年至1981年在美国伊利诺斯州Kincaid电厂进行，该核电厂的地理坐标为西经89.49°，北纬39.59°。海拔约180 m，周围是平坦的农场和少量湖泊，地表粗糙度约10 cm。示踪气体实验持续约350 h，示踪气体采用SF6，实验期间SF6从高187 m，直径9 m的电厂烟囱释放，监测器分别布置在距源0.5至50 km范围内的12条弧线上。数据集中包括了171 h，1 284组弧线-小时数据，每组数据都设置有质量指数表示数值的可信任程度。

丹麦国家环境研究所（NERI）开发的模型验证比对工具MVK中Kincaid数据集是包含数据最全面、应用最广泛的数据集，此外还有Indianapolis、Copenhagen、Lillestrøm三个大气扩散实验数据集。放射性核素的大气扩散可借助常规的实踪实验进行验证，如Ka-Hing Yau等应用Kincaid数据集验证了AUSTAL2000和CALPUFF扩散模型，EPA将该数据集应用于CALPUFF和其他大气扩散模型的验证，Steven R. Hanna将其用于ISC3、ADMS和AERMOD模型的验证比对，邹敬等也应用此数据集验证比对了实时在线核事故决策支持系统系统中RIMPUFF模型的验证比对^［4］。

目前国际上已经开展了一些关于空气质量模型一般评价方法的讨论，但是到目前为止还没有形成统一的评价方法和模型性能标准。应用统计评价方法分析验证模型中，MVK工具分析较为常用。

丹麦国家环境研究所（NERI）基于Hanna^［5-7］等的定量统计参数方法开发了模型验证工具MVK（Model Validation Kit）。该工具包包括Kincaid、Indianapolis、Copenhagen、Lillestrom四个大气扩散实验数据集，定量分析程序（BOOT）、残差分析程序（Residual Analysis）及二维绘图程序（Sigplot）等。实验观测值与模型预测值通过MVK定量统计分析、残差分析及定性图形分析的全面比对，可以有效的评价模型的准确性、适用性，揭示模型的错误与数据中的坏点。

目前BOOT程序已经涵盖了ASTM（American Society for Testing and Materials）的模型评价方法，诊断模型预测值的采用的统计量除平均值（Mean）、标准偏差（Sigma）外，还包括平均偏差（Bias）、比例偏差（FB）、几何平均偏差（MG）、标准化均方误差（NMSE）、几何平均方差（VG）、相关系数（r）及比例方差（FS），这些参数的计算公式如表1所示：

式中：CO为观测浓度；CP为模型预测浓度；σCo为实验观测浓度的标准差；σCp为模型预测浓度的标准差，单位均为μg∙m-3。

平均偏差Bias定义为预测平均值和观测平均值的差值，其数值可能是正数也可能是负数。当Bias为正数时，可能表示每一组数据的预测值均大于观测值，即模型具有高估的趋势，也可能表示预测值大于观测值数据对的偏差之和，大于预测值小于观测值数据对的偏差之和，因此Bias不具备表征模型高低估趋势的性能。其置信区间可以表示为：

式中：t_α取决于t分布的自由度和置信区间的振幅。

FB和MG都是反映系统误差的统计参数。其中FB更为常用，主要原因有2个：（1）FB是对称而有界的，可接受范围在－2（极端低估）到2（极端高估）之间；（2）FB是一个无量纲数，使用于涉及不同浓度量级比较的研究中比较方便。但对于大部分的空气污染物来说，其浓度接近于对数正态分布，因此个别的高观测浓度或预测浓度会对比例偏差FB产生显著的影响，而几何平均误差MG由于采用对数的处理方法，对很高的数值就不那么敏感了，但接近于0的浓度值会对MG影响显著。

NMSE是表征偏差的统计参量，对观测值与预测值的差值非常敏感，尤其峰值的偏差对NMSE的影响更大。一般来说，模型的评价与验证都是用模型预测值直接与实验观测值比对的，但是由于观测值的不确定性主要来自于大气的随机扰动，而模型预测值的不确定性主要来自于模型模拟的物理过程，两者来源不同，直接比较的方法有可能会得到误导性的结果，因此有的时候将系统误差和非系统（随机）误差分开考虑是必要的。总的NMSE可以表示为由系统误差引起的标准化均方误差NMSE_S与由非系统误差引起的标准化均方误差NMSE_U之和。当非系统误差为0时，NMSE取到最小值NMSE_S。

类似地，VG等于由系统误差引起的几何方差VG_S与由非系统误差引起的几何方差VG_U之和。

相关系数r是目前模型评价中使用最广泛的统计量。r为1时，表示预测值与测量值正相关，即数据对（CO，CP）分布在斜率为正数的直线上；r为-1时，则直线的斜率为负数；r为0时，表示预测值和测量值之间没有线性关系。

α因子表示满足1α≤CPCO≤α的数据所占的百分比（不包括（0，0）数据对）。α常取2和5。FAC对因为不受数据奇点的影响，是模型评价中最可靠的方法，可以作为消除极端值的工具和提高统计分析结果的方法。

MVK虽然在模型的验证与评价中得到了广泛的应用，但其应用仍具有一定的局限性，主要体现在以下几个方面：

（1）数量有限，只考虑了4个实验数据集；（2）仅适用于短程扩散模型；（3）只能处理一些相对简单的情景，主要适用于在均匀地表环境下单个点源释放无反应性气体的情景；（4）没有明确的解决由观测浓度的随机特性所产生的问题，程序只是直接将弧最大浓度和烟羽中心线浓度进行比较；（5）数据集主要关注的是弧最大浓度，只有部分数据集包括可信的侧风向积分浓度^［8］；（6）BOOT程序采用的统计量中并没有一个综合性的统计指标，验证人员只能根据统计结果得到模型在某一方面的统计评价，而没有总体的适用性评价。

本文应用MVK工具，将ARTM模拟结果与AERMOD模型、CALPUFF模型及Kincaid数据集进行比对。AERMOD模型和CALPUFF模型是《环境影响评价技术导则大气环境》（HJ2.2—2008）推荐的大气扩散模型，已经广泛应用放射性核素大气扩散分析^［9-12］，此外CALPUFF模型还能应用于核事故水体影响分析中^［13］。Kincaid数据集根据数据质量将数据进行分类，本文采用质量指示数分别为2和3的552组数据进行比对验证，数据分布在3~40 km范围内，模拟区域地形平坦。

ARTM是以德国空气质量标准推荐的扩散模式AUSTAL2000扩散模型为基础来计算大气中放射性污染物的扩散的，该模型系统主要包括诊断风场模型、雾羽抬升模型、浓度计算模型和核素沉积模型。模型的基本原理是拉格朗日烟团粒子模式，即把每个污染质点当成有标志的质点，粒子的运行轨迹可用下式表述：

式中：X为粒子的三维坐标分量；V为平均风速(u,v,w)；V'为湍流脉动速度分量(u',v',w')；t为时间序列；Δt为时间步长。

通过释放大量粒子，计算粒子的轨迹，这些粒子描述气载污染物在大气中的迁移扩散。粒子在流场中按平均风输送，同时又用一系列随机位移来模拟湍流扩散，这样就表达了平流和湍流扩散两种作用，最后由这些粒子在空间和时间上的总体分布估算出污染物的分布^［14］。

如图1所示，ARTM模型系统包括诊断风场模型、雾羽抬升模型、浓度计算模型和核素沉积模型。诊断风场模型结合地形和下垫面状况处理输入的气象数据，计算风场。浓度计算模型经过雾羽抬升模型的修正后输出每个网格的空气浓度。最后核素沉积模型输出每个网格的干、湿沉积浓度。

Figure 1.  Simulation process of ARTM model

实验的统计结果列于表2中，模型预测值与观测值越接近，则表示模型的预测性能越好。

表2中，ARTM模型的均值偏保守，其他两个模型均值偏小，这样的结果不仅与模型本身的属性有关，与实验条件也有很大关系，如ARTM因软件局限采用了平坦地形模拟，而AERMOD和CALPUFF均输入实际地形数据进行模拟，因此ARTM的预测结果在数值上会比实际地形略大。一般说来如果模型的统计结果满足－0.3<FB<0.3，0.7<MG<1.3，NMSE<1.5，VG<4，FAC2>0.5，则可以认为模型的模拟结果是有效的^［15］，ARTM的统计结果基本达到了这个标准。相关系数Cor是在检验模式预测值有效性中广泛应用的统计量，描述了两组变量相对于各自的均值以统一比例增加或减少的程度^［16］。ATRM模型的相关系数值略低，表明观测值与预测值相对于各自的均值变化的统一程度较差。此外ARTM模型比例方差数值上与理想值也是比较接近的。

模型预测值和观测值的对应的散点分布图如图2所示。数据点在图中呈大致的椭圆形带状散布，但一部分点偏离较远，造成相关系数r值偏低。图3为模型分位数图，分位数图是将模型预测值与观测值分别从高到低排序，然后将对应的点打印在图上，这个图的意义在于能反映模型的预测趋势^［17］。由图4柱状分析图可以看出，ARTM模型柱状图0 km~3 km距离段75%处约指向1.0，表明在该距离段有75%的数据的比值小于1.0。结合图1反映的预测趋势，说明ARTM模型在0 km~3 km距离段的预测有低估的趋势，在稍远处（5 km以外）ARTM模型的预测值有高估的趋势。

Figure 2.  Scatter plot of ARTM

Figure 3.  Quantile plot of ARTM

Figure 4.  Residual plot of ARTM

在大气扩散模型的验证评价中，需要结合模型的特点选择适合的实验数据集和验证评价方法。目前国际上应用广泛的示踪实验数据集和模型验证评价方法中，模型验证工具MVK（Model Validation Kit）采用了比较完善的统计分析方法和实用的图形工具，是众多验证评价方法中较为规范和成熟的。应用MVK验证评价放射性核素大气扩散模型ARTM的统计结果表明：ARTM的模拟结果基本上是可信赖的，3 km内的预测有低估的趋势，5 km外的预测有高估的趋势。

单一数据集只能代表有限的气象条件，对于模型验证评价应选取更多的数据集、综合采用多种验证评价方法做更广泛的比对。现有可靠性较高的数据集所能代表的源项、气象及地形条件仍十分有限，但现场扩散实验较为昂贵，为大气扩散模型适用范围和局限性的深入研究增加了困难，因此开展更多的高质量的示踪实验或者风洞实验用于模型验证比对也是十分必要的。