地下水型饮用水水源地污染源解析方法技术

本发明专利技术提供一种地下水型饮用水水源地污染源解析方法，包括：对采样数据进行预处理；计算预处理后的采样数据的特征值和特征向量；利用所述特征值矩阵和特征向量矩阵求得因子载荷矩阵和因子得分矩阵；对所述因子载荷矩阵和因子得分矩阵进行非负约束的旋转；使用最小二乘方法进行迭代运算重复所述对所述因子载荷矩阵和因子得分矩阵进行非负约束的旋转，确定因子载荷矩阵和因子得分矩阵。对所述因子载荷矩阵和因子得分矩阵进行非负约束的旋转，避免矩阵分解的结果中出现负值，使得因子载荷和因子得分具有可解释性和明确的物理意义，同时确定污染源和贡献，而不需要事先知道源成分谱，运用简便，从而有效消除测量检测的不确定性。

【技术实现步骤摘要】

本专利技术涉及水体检测
，特别是指一种地下水型饮用水水源地污染源解析方法。
技术介绍
近年来，就全国范围而言，我国地下水质量下降趋势明显，不宜作为饮用水水源的比例逐渐增大，即使存在地下水水源地饮用水指标并不一定受到直接污染而超标现象，但其水质可能已经不同程度承受污染风险，且水质同时受地下水运动演化影响，因此考虑选择反映水质演化的地下水水化学类型常规指标，进行污染来源解析。地下水水质指标在取样与分析过程中，由于取样和样品处理、试剂和水纯度、仪器量度和仪器洁净、采用的分析方法、测定过程以及数据处理等过程均会产生测量误差(系统误差，随机误差，过失误差)。从取样到分析结果计算误差都绝对存在，虽然在各个过程中进行质量控制，但无法完全消除不确定性的影响。
技术实现思路
本专利技术要解决的技术问题是提供一种地下水型饮用水水源地污染源解析方法，能够有效消除测量检测的消除不确定性。为解决上述技术问题，本专利技术的实施例提供一种地下水型饮用水水源地污染源解析方法，所述地下水型饮用水水源地污染源解析方法包括：对采样数据进行预处理，所述预处理包括缺失值处理、不确定性处理和数据合理性分析中的至少一种；计算预处理后的采样数据的特征值和特征向量，得到特征值矩阵和特征向量矩阵；利用所述特征值矩阵和特征向量矩阵求得因子载荷矩阵和因子得分矩阵；对所述因子载荷矩阵和因子得分矩阵进行非负约束的旋转；使用最小二乘方法进行迭代运算重复所述对所述因子载荷矩阵和因子得分矩阵进行非负约束的旋转，确定因子载荷矩阵和因子得分矩阵，所述因子载荷矩阵对应污染源指纹图谱，所述因子得分矩阵对应主要污染源贡献。优选的，对采样数据进行缺失值处理包括：丢弃采样数据中含缺失数据的记录；用全局变量或属性的平均值代替缺失数据；先根据欧式距离或相关分析来确定距离缺失数据样本最近的K个样本，将这K个值加权平均来估计该样本的缺失数据，其中k为大于零的整数；采用神经网络模型预测每一个缺失数据；对低于数据检测限的数据用数据检测限值或1/2检测限代替；优选的，对采样数据进行不确定性处理包括：利用下式计算采样数据的不确定性：U=(s×c)2+l2]]>其中，s为误差百分数；c为指标浓度值；l为因子数据检出限。优选的，对采样数据进行数据合理性分析，包括：采用下式对采样数据进行筛选：S/N=Σi=1n(xij-sij)2Σi=1nsij2]]>其中，S/N为信噪比，xij表示第i采样点第j个样品的浓度，sij表示第i采样点第j个样品的标准偏差。优选的，所述计算预处理后采样数据的特征值和特征向量，得到特征值矩阵和特征向量矩阵，包括：对预处理后的采样数据进行无量纲化；计算无量纲化后采样数据的协方差矩阵；利用雅各布方法计算协方差矩阵的特征值和特征向量，得到特征值矩阵和特征向量矩阵；优选的，所述计算无量纲化后采样数据的协方差矩阵，包括：利用下式计算无量纲化后采样数据的协方差矩阵：Z＝DD′其中，Z为协方差矩阵，D为无量纲化后的样品数据矩阵，D′为D的转置；利用雅各布方法计算协方差矩阵的特征值和特征向量，得到特征值矩阵和特征向量矩阵，协方差矩阵表示为：Z＝QEQ′其中，E为特征值矩阵，Q为特征向量矩阵，Q′为Q的转置。优选的，所述利用所述特征值矩阵和特征向量矩阵求得因子载荷矩阵和因子得分矩阵之前，包括：利用累计方差贡献率提取显著性因子；所述利用累计方差贡献率提取显著性因子的提取条件为：CV=Σj=1nλjΣj=1mλj×100%]]>＞第一限定值其中：n为显著性因子个数；m为污染物个数；λ为特征值。优选的，所述利用所述特征值矩阵和特征向量矩阵求得因子载荷矩阵和因子得分矩阵，所述因子载荷矩阵表示为：S＝QE1/2因子得分矩阵可表示为：C＝(S′S)-1S′D其中,S为因子载荷矩阵，C为因子得分矩阵，S’为S的转置。优选的，所述对所述因子载荷矩阵和因子得分矩阵进行非负约束的旋转，包括：利用下式对所述因子载荷矩阵和因子得分矩阵进行非负约束的旋转：S1＝ST1-1C1＝T1C其中，S1为旋转后的因子载荷矩阵；C1为旋转后的因子得分矩阵；T1为转换矩阵，且T1＝(C*C′)(C′C*)-1，C*为把C中的负值替换为零后的因子得分矩阵；利用因子载荷矩阵S1构建转换矩阵T2对因子载荷矩阵S1和因子得分矩阵C1继续旋转：S2＝S1T2C2＝T2-1C1其中：S2为二次旋转后的因子载荷矩阵；C2为二次旋转后的因子得分矩阵；T2为二次转换矩阵，且T2＝(S1′S1)-1(S1′S1*)，S1*为S1中的负值换为零后的因子载荷矩阵。优选的，所述使用最小二乘方法进行迭代运算重复所述对所述因子载荷矩阵和因子得分矩阵进行非负约束的旋转，确定因子载荷矩阵和因子得分矩阵，包括：重复所述对所述因子载荷矩阵和因子得分矩阵进行非负约束的旋转，因子载荷中负值的平方和小于预先设定的第二限定值时，输出因子载荷矩阵和因子得分矩阵。本专利技术的上述技术方案的有益效果如下：上述方案中，对所述因子载荷矩阵和因子得分矩阵进行非负约束的旋转，避免矩阵分解的结果中出现负值，使得因子载荷和因子得分具有可解释性和明确的物理意义，使用最小二乘方法进行迭代运算重复所述对所述因子载荷矩阵和因子得分矩阵进行非负约束的旋转，确定因子载荷矩阵和因子得分矩阵，能够同时确定污染源和贡献，而不需要事先知道源成分谱，运用简便，从而有效消除测量检测的不确定性。附图说明图1为本专利技术实施例的地下水型饮用水水源地污染源解析方法流程图。具体实施方式为使本专利技术要解决的技术问题、技术方案和优点更加清楚，下面将结合附图及具体实施例进行详细描述。如图1所示，本专利技术的实施例一种地下水型饮用水水源地污染源解析方法，所述地下水型饮用水水源地污染源解析方法包括：步骤101：对采样数据进行预处理，所述预处理包括缺失值处理、不确定性处理和数据合理性分析中的至少一种。步骤102：计算预处理后的采样数据的特征值和特征向量，得到特征值矩阵和特征向量矩阵。步骤103：利用所述特征值矩阵和特征向量矩阵求得因子载荷矩阵和因子得分矩阵。步骤104：对所述因子载荷矩阵和因子得分矩阵进行非负约束的旋转。步骤105：使用最小二乘方法进行迭代运算重复所述对本文档来自技高网...

【技术保护点】
一种地下水型饮用水水源地污染源解析方法，其特征在于，所述地下水型饮用水水源地污染源解析方法包括：对采样数据进行预处理，所述预处理包括缺失值处理、不确定性处理和数据合理性分析中的至少一种；计算预处理后的采样数据的特征值和特征向量，得到特征值矩阵和特征向量矩阵；利用所述特征值矩阵和特征向量矩阵求得因子载荷矩阵和因子得分矩阵；对所述因子载荷矩阵和因子得分矩阵进行非负约束的旋转；使用最小二乘方法进行迭代运算重复所述对所述因子载荷矩阵和因子得分矩阵进行非负约束的旋转，确定因子载荷矩阵和因子得分矩阵，所述因子载荷矩阵对应污染源指纹图谱，所述因子得分矩阵对应主要污染源贡献。

【技术特征摘要】
1.一种地下水型饮用水水源地污染源解析方法，其特征在于，所述地下水型饮用水水
源地污染源解析方法包括：
对采样数据进行预处理，所述预处理包括缺失值处理、不确定性处理和数据合理性分
析中的至少一种；
计算预处理后的采样数据的特征值和特征向量，得到特征值矩阵和特征向量矩阵；
利用所述特征值矩阵和特征向量矩阵求得因子载荷矩阵和因子得分矩阵；
对所述因子载荷矩阵和因子得分矩阵进行非负约束的旋转；
使用最小二乘方法进行迭代运算重复所述对所述因子载荷矩阵和因子得分矩阵进行
非负约束的旋转，确定因子载荷矩阵和因子得分矩阵，所述因子载荷矩阵对应污染源指纹
图谱，所述因子得分矩阵对应主要污染源贡献。
2.根据权利要求1所述的地下水型饮用水水源地污染源解析方法，其特征在于，对采样
数据进行缺失值处理包括：
丢弃采样数据中含缺失数据的记录；
用全局变量或属性的平均值代替缺失数据；
先根据欧式距离或相关分析来确定距离缺失数据样本最近的K个样本，将这K个值加权
平均来估计该样本的缺失数据，其中k为大于零的整数；
采用神经网络模型预测每一个缺失数据；
对低于数据检测限的数据用数据检测限值或1/2检测限代替。
3.根据权利要求1所述的地下水型饮用水水源地污染源解析方法，其特征在于，对采样
数据进行不确定性处理包括：
利用下式计算采样数据的不确定性：
其中，s为误差百分数；c为指标浓度值；l为因子数据检出限。
4.根据权利要求1所述的地下水型饮用水水源地污染源解析方法，其特征在于，对采样
数据进行数据合理性分析，包括：
采用下式对采样数据进行筛选：
其中，S/N为信噪比，xij表示第i采样点第j个样品的浓度，sij表示第i采样点第j个样品
的标准偏差。
5.根据权利要求1所述的地下水型饮用水水源地污染源解析方法，其特征在于，所述计
算预处理后采样数据的特征值和特征向量，得到特征值矩阵和特征向量矩阵，包括：
对预处理后的采样数据进行无量纲化；
计算无量纲化后采样数据的协方差矩阵；
利用雅各布方法计算协方差矩阵的特征值和特征向量，得到特征值矩阵和特征向量矩
阵。
6.根据权利要求5所述的地下水型饮用水水源地污染源解析方法，其特征在于，所述计
算无量纲化后采样数据的协方差矩阵，包括：
利用下式计算无量纲化...

【专利技术属性】
技术研发人员：左锐，张宏凯，王金生，滕彦国，杨洁，翟远征，
申请(专利权)人：北京师范大学，
类型：发明
国别省市：北京;11