[发明专利]一种估算细颗粒物及其前体物环境容量的新方法

摘要

本发明提供一种估算细颗粒物及其前体物环境容量的新方法，包括得到大气静态总容量；得到污染物水平方向输送总量；得到干沉降总量和湿沉降总量，再求和得到沉降总量；得到水平扩散总量和垂直扩散总量，再求和得到扩散清除总量；将水平方向输送总量、沉降总量和扩散清除总量求和，得到大气动态容量；求和大气静态总容量和大气动态容量，再进行修正，得到污染物订正后的环境容量。优点为从污染物生成、转化、消亡等物理化学过程量化大气对污染物的容纳能力，并且细颗粒物及其前体物的环境容量具有时空动态特征，采用本发明提供的方法估算大气污染物环境容量时，结果更为逼近真实监测值。还具有估算过程简单快速、计算量小的优点。

主权项

一种估算细颗粒物及其前体物环境容量的方法，其特征在于，包括以下步骤：步骤1，对于被估算的三维评价区域，建立空间投影坐标系XYZ；其中，X方向和Y方向位于水平面；Z方向为垂直方向；步骤2，在水平方向，按设定的网格分辨率，将所述三维评价区域进行网格化处理，将其格网化为x*y个网格，其中，x为X方向网格数量，y为Y方向网格数量；在垂直方向，按照设定的垂直方向层数，再将已水平方向网格化的所述三维评价区域划分为z层，其中，z为Z方向的层数；步骤3，利用中尺度气象模式模拟计算所述三维评价区域评价年的气象要素逐小时均值，获得第三代空气质量模式模拟计算必需输入的气象参数，包括：每个网格在x方向的风速ui,j,k以及每个网格在y方向的风速vi,j,k；步骤4，利用第三代空气质量模式模拟计算所述三维评价区域评价年的大气污染物浓度逐小时均值，获得步骤2得到的每个网格的空气质量参数，包括：污染物成分、污染物的浓度ci,j,k、干沉降过程清除的污染物浓度湿沉降过程清除的污染物浓度水平扩散过程清除的污染物浓度和垂直扩散过程清除的污染物浓度步骤5，读取基于国家环境空气质量标准等级的污染物标准浓度限值cstd以及各层评价柱体的体积Vk，基于公式1计算得到三维评价空间的大气静态总容量Qstatic：Qstatic=Σk=1zcstd×Vk---(1)]]>步骤6，对于任意的第i,j,k网格，读取该网格的污染物浓度ci,j,k、在x方向的风速ui,j,k、在y方向的风速vi,j,k、在x方向网格分辨率dxi,j、在y方向的网格分辨率dyi,j和网格高度dzi,j,k；再读取到第i‑1,j,k网格的在x方向的风速ui‑1,j,k、在y方向的网格分辨率dyi‑1,j和网格高度dzi‑1,j,k；再读取到第i,j‑1,k网格的在y方向的风速vi,j‑1,k、在x方向网格分辨率dxi,j‑1和网格高度dzi,j‑1,k；然后，基于公式2计算得到第i,j,k网格的网格水平输送量Fi,j,k：Fi,j,k=ci,j,k×ui,j,k×dyi,j×dzi,j,k+ci,j,k×ui-1,j,k×dyi-1,j×dzi-1,j,k+ci,j,k×vi,j,k×dxi,j×dzi,j,k+ci,j,k×vi,j-1,k×dxi,j-1×dzi,j-1,k---(2)]]>在分别计算得到每个网格的网格水平输送量后，基于公式3，对每个网格的网格水平输送量进行求和计算，得到三维评价空间的污染物水平方向输送总量F：F=Σi=1xΣj=1yΣk=1zFi,j,k---(3)]]>步骤7，对于任意的第i,j,k网格，读取干沉降过程中该网格清除的污染物浓度和该网格保留的污染物浓度ci,j,k；然后，基于公式4计算得到该网格的干沉降量Ddryi,j,k=Vdryi,j,k×cstd=cdryi,j,kcdryi,j,k+ci,j,k×cstd---(4)]]>其中，代表干沉降速率；在分别计算得到每个网格的干沉降量后，基于公式5，对每个网格的干沉降量进行求和计算，得到三维评价空间的干沉降总量Ddry：Ddry=Σi=1xΣj=1yΣk=1zDdryi,j,k---(5)]]>步骤8，对于任意的第i,j,k网格，读取湿沉降过程中该网格清除的污染物浓度以及该网格保留的污染物浓度ci,j,k；然后，基于公式6计算得到该网格的湿沉降量Dweti,j,k=Vweti,j,k×cstd=cweti,j,kcweti,j,k+ci,j,k×cstd---(6)]]>其中，代表湿沉降速率；在分别计算得到每个网格的湿沉降量后，基于公式7，对每个网格的湿沉降量进行求和计算，得到三维评价空间的湿沉降总量Dwet：Dwet=Σi=1xΣj=1yΣk=1zDweti,j,k---(7)]]>步骤9，基于公式8，将步骤7得到的干沉降总量Ddry和步骤8得到的湿沉降总量Dwet进行求和计算，得到三维评价空间的沉降总量Dep：Dep＝Ddry+Dwet                      (8)步骤10，对于任意的第i,j,k网格，读取水平扩散过程中该网格清除的污染物浓度和该网格保留的污染物浓度ci,j,k；然后，基于公式9计算得到该网格的水平扩散量Dhori,j,k=Vhori,j,k×cstd=chori,j,kchori,j,k+ci,j,k×cstd---(9)]]>其中，代表水平扩散速率；在分别计算得到每个网格的水平扩散量后，基于公式10，对每个网格的水平扩散量进行求和计算，得到三维评价空间的水平扩散总量Dhor：Dhor=Σi=1xΣj=1yΣk=1zDhori,j,k---(10)]]>步骤11，对于任意的第i,j,k网格，读取垂直扩散过程中该网格清除的污染物浓度和该网格保留的污染物浓度ci,j,k；然后，基于公式11计算得到该网格的垂直扩散量Dveri,j,k=Vveri,j,k×cstd=cveri,j,kcveri,j,k+ci,j,k×cstd---(11)]]>其中，代表垂直扩散速率；在分别计算得到每个网格的垂直扩散量后，基于公式12，对每个网格的垂直扩散量进行求和计算，得到三维评价空间的垂直扩散总量Dver：Dver=Σi=1xΣj=1yΣk=1zDveri,j,k---(12)]]>步骤12，基于公式13，将步骤10得到的水平扩散总量Dhor和步骤11得到的垂直扩散总量Dver进行求和计算，得到三维评价空间的扩散清除总量Dif：Dif＝Dhor+Dver                      (13)步骤13，基于公式14，将步骤6得到的污染物水平方向输送总量F、步骤9得到的沉降总量Dep和步骤12得到的扩散清除总量Dif进行求和计算，得到三维评价空间的大气动态容量Qdynamic；Qdynamic＝F+Dep+Dif                  (14)步骤14，基于公式15，将步骤5得到的大气静态总容量Qstatic和步骤13得到的大气动态容量Qdynamic进行求和计算，得到三维评价空间的污染物初始环境容量Q；Q＝Qstatic+Qdynamic                     (15)步骤15，读取到污染物转化率α，基于公式16，得到污染物订正后的环境容量Q’：Q'＝Q×1/α                         (16)；步骤15中，当污染物初始环境容量Q为一次细颗粒物初始环境容量时，按公式16得到的环境容量Q’即为一次细颗粒物环境容量，并且，此时α取值为1；当污染物初始环境容量Q为硫酸盐时，按公式16得到的环境容量Q’为二氧化硫环境容量，并且，此时α取值为：0＜α＜1；当污染物初始环境容量Q为硝酸盐时，按公式16得到的环境容量Q’为氮氧化物环境容量，并且，此时α取值为：0＜α＜1。