关于GIS数据转SMOKE源清单求助大家

lm8005507771

目前用GIS做了某市的交通源清单，但现在只是网格数据，想要把它用在CMAQ中，结合MIX源清单模拟某市的交通源影响情况，查资料找不太到从网格数据转到CMAQ可以识别的SMOKE数据的方法，求助大家，多谢多谢！

aslamic

空气质量模型应用的核心难点就是清单的准备，两个途径：
1.        攻破smoke模型，毕竟这是cmaq模型体系的默认清单处理模型，但是由于SMOKE适应的是美国的清单数据体系，我国没有这样的体系，所以你需要根据他的文件特征解决IDA或者其他格式的排放清单文件编制（是个很简单的fortran问题，不同字段对应不同的格式，需要注意一些单位的转换），并完成相关的支撑文件，如COSTCY文件、时间索引、时间分配文件、物种分配文件、物种分配索引文件等，通常这些东西弄完你的smoke就算搞清楚三分之一了，这是个硬骨头，不要指望别的课题组会给你多少帮助，但是不是不可能完成的事情，毕竟SMOKE手册写得很清楚；
2.        跳过smoke，直接写成ioapi格式的数据，所谓的ioapi其实是一个特殊约定的netcdf格式文件，有一些特殊的attribute用来给ioapi读取用，相关的格式可以从cmaq示例文件里面找到，但是因为你跳过了smoke但是又不能不管smoke实现的功能，这些东西都需要你自己做，包括但不限于（默认你不是gis高手，因为gis高手会怎么做我不知道，我不懂gis）：
a)        空间分配，你需要从griddot2d和gridcro2d文件里面读取经纬度信息，获得每个网格的中心经纬度和对应的边界经纬度，将你的网格化清单通过resample也好怎么也好放到对应的网格里面；
b)        物种分配，根据你选择的机制（如cb05、cb06、ae5、ae6）将VOCs、NOx、PM2.5分为必要的组分，涉及单位换算（气态污染物为mol/s，颗粒态为g/s，但是清单是t/a）
c)        垂直分配，针对点源你需要利用布里吉斯公式或者参考薛文博等的因子法将点源排放放到不同的高度里面，模型层对应的高度在metcro3d文件中获取；
d)        时间分配，你需要根据排放特征（从活动水平、在线监测等处获取），将年排放分成具有明显时间属性的月排放、日排放（结合周几来做）和小时排放；
e)        文件生成，上述内容可以说主要是编程技巧和简单数学的范畴，那么现在你需要自己编程生成cmaq所需的ioapi数据了，你可以用matlab、ncl、R等提供netcdf读写能力的程序手动实现ioapi的约定，也可以用c、fortran直接调用ioapi来写，到了这里的话就可以测试了，如果测试通过，这件事就成了。
回到开头，不管你用1还是2，都是很难的一件事情，除开技术本身，更多的是需要你专心研究这件事，从其他课题组或者人那里能够得到的帮助是很有限的，，每个课题组都有自己的处理方法。
至于我，当然是第一种和第二种都会咯，其中第二种打算最近开源出来，但是实在是没空写文档，估计只有来年了。

L0725

你好，我现在也到了这一步，能用GIS做出网格清单，但是怎么转成可以被CMAQ读取的清单，可否给些指点。

lm8005507771

aslamic 发表于 2019-11-21 19:21
空气质量模型应用的核心难点就是清单的准备，两个途径：
1.        攻破smoke模型，毕竟这是cmaq模型体系的默认清 ...

高手！！！

lm8005507771

L0725 发表于 2019-11-21 14:55
你好，我现在也到了这一步，能用GIS做出网格清单，但是怎么转成可以被CMAQ读取的清单，可否给些指点。

目前应该就是高手说的两种途径，我是用了第一种方法，把smoke弄懂了三分之一的样子（参加了一个smoke培训，然后又看了看用户手册），然后用smoke处理了源清单，cmaq就可以识别了，如果不是算法大神的话第一种方法相对简单一点，毕竟有现成的手册给你看，网上都能找到（没记错的话就是CMAQ官网），从安装到使用都有。

L0725

aslamic 发表于 2019-11-21 19:21
空气质量模型应用的核心难点就是清单的准备，两个途径：
1.        攻破smoke模型，毕竟这是cmaq模型体系的默认清 ...

果然高手，期待大作。

L0725

lm8005507771 发表于 2019-11-22 08:47
目前应该就是高手说的两种途径，我是用了第一种方法，把smoke弄懂了三分之一的样子（参加了一个smoke培训 ...

也就说GIS网格化的清单和cmaq能够识别的清单还是有差距的？归根结底还是从smoke入手的

L0725

lm8005507771 发表于 2019-11-22 08:47
目前应该就是高手说的两种途径，我是用了第一种方法，把smoke弄懂了三分之一的样子（参加了一个smoke培训 ...

也就说GIS网格化的清单和cmaq能够识别的清单还是有差距的？归根结底还是从smoke入手的

lm8005507771

L0725 发表于 2019-11-24 19:25
也就说GIS网格化的清单和cmaq能够识别的清单还是有差距的？归根结底还是从smoke入手的

GIS网格化的清单，就是只是一个二维分布的清单，应该是一个按照某个规则把污染物总量分配到格点的矩阵（我做的交通源就按路网分配，不过这个格点分配比例是smoke里要用的，很重要的一个参数），但是一个文本格式的东西cmaq是不认识的，最起码要转成ioapi格式的，另外还要考虑物种分配、时间分配等等，我理解的其中物种分配应该是smoke最厉害的地方，里面考虑了几种不同机制，每种机制又涉及几十上百个化学方程式。所以从smoke入手的话比较方便

L0725

lm8005507771 发表于 2019-11-25 08:53
GIS网格化的清单，就是只是一个二维分布的清单，应该是一个按照某个规则把污染物总量分配到格点的矩阵（ ...

好的，非常感谢！