[发明专利]增加气象模式接口的中尺度大气光化学污染模拟预测算法

权利要求书

1.一种增加气象模式接口的中尺度大气光化学污染模拟预测算法，包括以下步骤：

步骤1：采用CALGRID化学模式，考虑大气化学反应、大气输送与扩散、沉降、地面面源和高架排放源的影响，并对化学物种原始浓度中物理量的平均量分解为脉动量，得到化学物种浓度变化方程如式(1)，

式中，C是化学物种平均浓度，V是三维风矢量的平均量，K是湍流扩散系数，E是污染源排放，是由于沉降引起的物种浓度变化，PCHEM是化学产生率，LCHEM是化学损失率；式(1)中二阶湍流扩散项▽·(K▽C)通过湍流扩散系数K理论闭合转化得到；式(1)中右侧多项式依次分别为平流项、垂直扩散项、源项、沉降项、化学变化项；

步骤2：由于平流项、垂直扩散项、源项、沉降项、化学变化项各过程的特征时间不一致，通过算子分列时间积分对各过程进行积分，

Cn+1＝AxAyAzAcAzAyAxCn (2)

式中，Ax、Ay是水平输送扩散算子；Az是垂直输送扩散、物质输入和物理损耗算子，Ac是化学反应算子；

步骤3：将气象模式MM5或气象模式WRF的水平差分跳点网格ArakawaB和C格点气象数据插值到式(1)中CALGRID化学模式的水平差分非跳点网格ArakawaA的格点上；

步骤4：将气象模式MM5或气象模式WRF的垂直σ地形跟随坐标系统的气象数据插值到式(1)中CALGRID化学模式的高度Z坐标系统；

步骤5：通过气象模式MM5或WRF输出下垫面类型，并从外部资料集中获取下垫面的包括地表粗糙度和叶面积指数LAI的特征参数，所述外部资料集包括：地表粗糙度、NASA/GSFC空间分辨率为1°的LAI资料；

步骤6：通过气象模式MM5或WRF的预报场，诊断获得式(1)中CALGRID化学模式所需的边界层气象参数，所述边界层气象参数包括：PGT稳定度分类、摩擦速度u*、对流速率尺度w*、莫宁－奥布霍夫长度L；

步骤7：输入步骤6得到的边界层气象参数，通过改进后CALGRID化学模式(1)计算并输出污染物的输送扩散数据。

2.根据权利要求1所述增加气象模式接口的中尺度大气光化学污染模拟预测算法，其特征在于：步骤4中，对于水平风速u、v和温度t的插值通过以下方法求取：

首先，根据σ值，σ＝(p0-pt)/(ps-pt)，以及模式地面气压ps、模式顶气压pt计算水平风速u、v、温度t和垂直风速w变量在σ坐标系中所处的等压面；

然后，利用压高公式计算等压面对应的高度；

最后，利用线性插值，得到CALGRID坐标系统各个层次对应的数值。