环境空气颗粒物来源解析正定矩阵因子分解模型计算技术指南.docxVIP

1

HJ1353—2024

环境空气颗粒物来源解析正定矩阵因子分解模型计算技术指南

1适用范围

本标准规定了正定矩阵因子分解（PMF）模型计算的数据准备、基础计算、旋转计算、颗粒物源贡献计算、结果合理性判断等技术要求。

本标准适用于应用PMF模型对环境空气颗粒物进行来源解析相关的计算工作。

2规范性引用文件

本标准引用了下列文件或其中的条款。凡是注明日期的引用文件，仅注日期的版本适用于本标准。凡是未注日期的引用文件，其最新版本（包括所有的修改单）适用于本标准。

HJ630环境监测质量管理技术导则

HJ1352环境空气颗粒物来源解析基于手工监测的受体模型法监测数据处理与检验

技术规范

《环境空气颗粒物来源解析监测技术方法指南》（监测函〔2020〕8号）

3术语和定义

下列术语和定义适用于本标准。

3.1

环境受体ambientreceptor

受到污染源排放影响的环境空气，简称受体。

3.2

正定矩阵因子分解模型positivematrixfactorizationmodel（PMF）

将样本数据的X矩阵分解为因子贡献G矩阵和因子谱F矩阵，对因子谱F矩阵进行识别，并定量计算样本的因子贡献的一种多元因子分析类模型。

3.3

污染源化学成分谱sourcechemicalprofile

特定污染源类排放颗粒物的相对稳定的化学组分信息（g/g），简称源成分谱。源成分谱应包含该类污染源的标识组分。

3.4

污染源贡献sourcecontribution

通过受体模型计算得到的污染源对环境空气颗粒物的贡献，用贡献浓度(μg/m3）或贡献占比表示。3.5

污染源标识组分sourcechemicaltracer

污染源化学成分谱中对该源类有指示作用的一种或多种颗粒物化学组分，又称示踪组分。源标识组分是区别该源类与其他源类的重要标识物，每种污染源有各自的标识组分。

2

HJ1353—2024

3.6

因子factor

PMF模型计算提取的G矩阵和F矩阵的每列信息，其中F矩阵的每列可识别为特定污染源类，G矩阵的每列可计算得到特定污染源类的贡献。

3.7

Q值函数functionQvalue

反映PMF模型计算结果质量的一种目标函数。

3.8

残差residual

颗粒物化学组分的实际测试值与PMF模型计算值之间的差值。

3.9

基础计算baserun

在计算过程中不加入任何数学、物理意义的约束条件，只根据模型Q值函数的收敛条件（Q值最小）得到因子贡献G矩阵和因子谱F矩阵的一种PMF模型计算方法。

3.10

旋转计算rotationalrun

对于基础计算得到的因子贡献G矩阵和因子谱F矩阵，加入一定的约束条件再次进行计算，转换成另一对矩阵（G*和F*），使各因子谱的标识组分更加突出，从而解决因子谱难以识别问题的一种PMF模型计算方法。

4模型原理和计算流程

4.1模型原理

在环境空气颗粒物来源解析中，PMF模型计算是将受体颗粒物化学组分浓度矩阵Xn×m因子化，分解为两个因子矩阵，G(n×p)和F(p×m)，以及一个残差矩阵E(n×m)（见公式1），在G矩阵和F矩阵中的数值都是正值（即非负限制）的限定下，通过定义一个目标Q值函数（见公式2），选择Q值函数值最小时的G矩阵（因子贡献矩阵）和F矩阵（因子谱矩阵）计算结果，再利用污染源标识组分将F矩阵识别为不同的源类，对G矩阵进行多元线性回归计算，从而得出不同源类对环境样品的贡献。

(n题m)(n题p)(p题

(n题m)(n题p)(p题m)(n题m)

式中：X(n×m)——受体颗粒物化学组分浓度X矩阵；

G(n×p)——因子贡献G矩阵；

F(p×m)——因子谱F矩阵；

E(n×m)——残差E矩阵；

n——样品个数；

m——化学成分种数；

p——解析的因子（污染源）数目。

Q=(Eij/σij）2

式中：Q——目标Q值函数；

Eij——第i个样品中第j个化学成分的残差；

σij——第i个样品中第j个化学成分的不确定度；

（1）

（2）

3

反复优化计算基础计算调整参数结果误差评估旋转计算调整参数结果误差评估数据处理数据检验基础计算运行结果判断分析旋转计算运行结果判断分析HJ

反复优化计算

基础