PMF5.0计算

6.3 EPA PMF 5.0 计算步骤

6.3.1 软件下载

可在美 国 EPA 网站免 费下载 ，下载地址 ： https://www.epa.gov/air-research/positive-matrix-factorization-model-environmental-data-analyses。

计算机（Win7以上系统，3.1 GHz 以上处理器，4 GB 以上 RAM，需安装 framework4.0 以上版本）。

6.3.3 数据准备

3 种格式的数据文件可输入至 EPA PMF5.0 模型，包括 tab 间隔的 txt 文件、逗号间隔的csv 文件和 excel 文件，推荐使用 excel 文件。

PMF 模型需要输入两个数据文件：（1）颗粒物及组分质量浓度的数据文件，（2）颗粒物及组分质量浓度不确定度的数据文件。两个数据文件内容需要对应。第一行为组分名称，组分名称必须是唯一，中间不能含有逗号。第一列为样品编号，如样品 ID、日期或时间、样品 ID 加上日期时间等。两个数据文件均不允许有数据缺失、0 值和负值的情况。对于样品部分组分浓度未检出的情况，可用1/2 检出限来代替；对于样品数据中部分组分数据缺失的情况，建议输入-999 来标识，对应的不确定度可暂用5/6 检出限代替（模型计算时会自动替代）。

6.3.4 数据导入和初步检验

6.3.4.1  数据导入

启动 PMF 模型软件，进入“Model Data”界面（见图 15）。在“Data Files”的“Input Files”下选择样品组分浓度文件和不确定度文件的输入路径，点击“Load”（导入）选择相应的数据文件。在“Output Files”中，设置确定输出文件的路径、文件类型和输出文件前缀（Output File Profix），输出文件前缀（Output File Profix）只能是数字或字母。如对数据文件进行了修改，则需要重新加载输入文件。

图15  数据导入的设置界面

PMF 模型中可以保存新的配置文件或者载入之前的配置文件。配置文件中的信息包括“Data Files” 中的各类设置（例如，输入文件，输出文件位置和输出文件类型），“Concentration/Uncertainty”（浓度/不确定度）的设置以及“Base Model Run”（模型基础计算）、“Fpeak Rotation”（F 矩阵峰值旋转）、“Constrained Model Runs”（模型约束计算）中的所有运行参数。

6.3.4.2  初步检验

（1）在“Concentration/Uncertainty”中对受体化学组分数据的权重、异常值等进行设置或处理（界面见图 16）。根据 S/N（信噪比）值确定 PMF 中组分的计算权重，权重包括Bad（差）、Weak（弱）和 Strong（强）三个水平。对于被选择 Strong 的组分，按输入的浓度和不确定度参与计算；对于被选择 Weak 的组分，模型在计算时会自动增加该组分的不确定度，进而而降低该组分浓度的权重；对于被选择 bad 的组分，则不参与计算。建议

将颗粒物的关键标识组分（如 OC、EC、SO42-、NO3-、NH4+、Al、Si、Ca、K、Fe 等组分）设为 Strong。其他非关键标识组分（如：Cl-、Na、Mg、Cu、Cd、Cr 等组分）建议设置Strong，但是在分析不确定性较高时，可设置成 Weak 参与计算。当某一组分的缺失数据较

多时，可通过采用设定较大的不确定度来降低该组分的计算权重。此外，计算各因子对颗粒物的贡献时，需要将颗粒物质量浓度数据输入模型，并将颗粒物质量浓度设置为 Total Variable（总变量），其计算权重默认为 Weak。

（2）在“Concentration Scatter