CN109470667A - 一种结合水质参数和三维荧光光谱进行污染物溯源的方法 - Google Patents
一种结合水质参数和三维荧光光谱进行污染物溯源的方法 Download PDFInfo
- Publication number
- CN109470667A CN109470667A CN201811355655.5A CN201811355655A CN109470667A CN 109470667 A CN109470667 A CN 109470667A CN 201811355655 A CN201811355655 A CN 201811355655A CN 109470667 A CN109470667 A CN 109470667A
- Authority
- CN
- China
- Prior art keywords
- source
- water quality
- pollutant
- quality parameter
- fluorescence spectrum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 103
- 238000000034 method Methods 0.000 title claims abstract description 36
- 239000003344 environmental pollutant Substances 0.000 title claims abstract description 32
- 231100000719 pollutant Toxicity 0.000 title claims abstract description 32
- 238000002189 fluorescence spectrum Methods 0.000 title claims abstract description 31
- 238000001514 detection method Methods 0.000 claims abstract description 14
- 238000012544 monitoring process Methods 0.000 claims abstract description 13
- 238000005070 sampling Methods 0.000 claims description 16
- 238000004458 analytical method Methods 0.000 claims description 15
- QIVBCDIJIAJPQS-VIFPVBQESA-N L-tryptophane Chemical compound C1=CC=C2C(C[C@H](N)C(O)=O)=CNC2=C1 QIVBCDIJIAJPQS-VIFPVBQESA-N 0.000 claims description 12
- QIVBCDIJIAJPQS-UHFFFAOYSA-N Tryptophan Natural products C1=CC=C2C(CC(N)C(O)=O)=CNC2=C1 QIVBCDIJIAJPQS-UHFFFAOYSA-N 0.000 claims description 12
- 238000013459 approach Methods 0.000 claims description 10
- 238000009434 installation Methods 0.000 claims description 6
- PUKLDDOGISCFCP-JSQCKWNTSA-N 21-Deoxycortisone Chemical compound C1CC2=CC(=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@@](C(=O)C)(O)[C@@]1(C)CC2=O PUKLDDOGISCFCP-JSQCKWNTSA-N 0.000 claims description 5
- FCYKAQOGGFGCMD-UHFFFAOYSA-N Fulvic acid Natural products O1C2=CC(O)=C(O)C(C(O)=O)=C2C(=O)C2=C1CC(C)(O)OC2 FCYKAQOGGFGCMD-UHFFFAOYSA-N 0.000 claims description 5
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 5
- 239000002509 fulvic acid Substances 0.000 claims description 5
- 229940095100 fulvic acid Drugs 0.000 claims description 5
- 229910052698 phosphorus Inorganic materials 0.000 claims description 5
- 239000011574 phosphorus Substances 0.000 claims description 5
- XKMRRTOUMJRJIA-UHFFFAOYSA-N ammonia nh3 Chemical compound N.N XKMRRTOUMJRJIA-UHFFFAOYSA-N 0.000 claims description 4
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 3
- 229930002868 chlorophyll a Natural products 0.000 claims description 3
- ATNHDLDRLWWWCB-AENOIHSZSA-M chlorophyll a Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 ATNHDLDRLWWWCB-AENOIHSZSA-M 0.000 claims description 3
- 239000012535 impurity Substances 0.000 claims description 3
- 229910052760 oxygen Inorganic materials 0.000 claims description 3
- 239000001301 oxygen Substances 0.000 claims description 3
- 229910021642 ultra pure water Inorganic materials 0.000 claims description 3
- 239000012498 ultrapure water Substances 0.000 claims description 3
- MMDJDBSEMBIJBB-UHFFFAOYSA-N [O-][N+]([O-])=O.[O-][N+]([O-])=O.[O-][N+]([O-])=O.[NH6+3] Chemical compound [O-][N+]([O-])=O.[O-][N+]([O-])=O.[O-][N+]([O-])=O.[NH6+3] MMDJDBSEMBIJBB-UHFFFAOYSA-N 0.000 claims description 2
- DWNBOPVKNPVNQG-LURJTMIESA-N (2s)-4-hydroxy-2-(propylamino)butanoic acid Chemical compound CCCN[C@H](C(O)=O)CCO DWNBOPVKNPVNQG-LURJTMIESA-N 0.000 claims 1
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 239000004744 fabric Substances 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 239000003365 glass fiber Substances 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 3
- 238000000556 factor analysis Methods 0.000 abstract description 2
- 238000001917 fluorescence detection Methods 0.000 abstract 1
- 239000011159 matrix material Substances 0.000 description 15
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 6
- YBJHBAHKTGYVGT-ZKWXMUAHSA-N (+)-Biotin Chemical compound N1C(=O)N[C@@H]2[C@H](CCCCC(=O)O)SC[C@@H]21 YBJHBAHKTGYVGT-ZKWXMUAHSA-N 0.000 description 5
- 108090000623 proteins and genes Proteins 0.000 description 5
- FEPMHVLSLDOMQC-UHFFFAOYSA-N virginiamycin-S1 Natural products CC1OC(=O)C(C=2C=CC=CC=2)NC(=O)C2CC(=O)CCN2C(=O)C(CC=2C=CC=CC=2)N(C)C(=O)C2CCCN2C(=O)C(CC)NC(=O)C1NC(=O)C1=NC=CC=C1O FEPMHVLSLDOMQC-UHFFFAOYSA-N 0.000 description 5
- QGZKDVFQNNGYKY-UHFFFAOYSA-N Ammonia Chemical compound N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 description 4
- 230000005284 excitation Effects 0.000 description 4
- 239000011521 glass Substances 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 229910052757 nitrogen Inorganic materials 0.000 description 3
- 229910021529 ammonia Inorganic materials 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000012528 membrane Substances 0.000 description 2
- 238000012821 model calculation Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 239000003403 water pollutant Substances 0.000 description 2
- 238000003911 water pollution Methods 0.000 description 2
- 239000002028 Biomass Substances 0.000 description 1
- 229910002651 NO3 Inorganic materials 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000001506 fluorescence spectroscopy Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000010842 industrial wastewater Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- QJGQUHMNIGDVPM-UHFFFAOYSA-N nitrogen group Chemical group [N] QJGQUHMNIGDVPM-UHFFFAOYSA-N 0.000 description 1
- 229910000069 nitrogen hydride Inorganic materials 0.000 description 1
- 230000035764 nutrition Effects 0.000 description 1
- 235000016709 nutrition Nutrition 0.000 description 1
- 239000005416 organic matter Substances 0.000 description 1
- 239000002957 persistent organic pollutant Substances 0.000 description 1
- FGIUAXJPYTZDNR-UHFFFAOYSA-N potassium nitrate Chemical compound [K+].[O-][N+]([O-])=O FGIUAXJPYTZDNR-UHFFFAOYSA-N 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 239000010865 sewage Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/64—Fluorescence; Phosphorescence
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/18—Water
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/64—Fluorescence; Phosphorescence
- G01N2021/6417—Spectrofluorimetric devices
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/20—Controlling water pollution; Waste water treatment
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Engineering & Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
本发明提供了一种结合水质参数和三维荧光光谱进行污染物溯源的方法，步骤包括：水质检测步骤、荧光检测步骤、污染源类型识别步骤以及污染源确定步骤。该污染物溯源的方法结合PMF模型因子解析结果与特征荧光组分荧光强度比值，可快速判断水体环境是否受到人为活动影响，并同时判断其污染来源类型，再结合相关区域水质监测布点情况及水文参数定向筛查污染源，最终通过对比相应污染源检测结果明确水体主要污染来源，具有操作简便、准确度高、重现性好等优点，可以快速、有效识别水体环境中的典型污染源，为相关污染的管控与治理提供理论指导。
Description
技术领域
本发明涉及一种污染物溯源的方法，尤其是一种结合水质参数和三维荧光光谱进行污染物溯源的方法。
背景技术
随着人类社会高速发展，水体环境受到严重污染，其污染主要分为自然因素影响导致的水质变化和人为因素影响导致的水质恶化，如生活源的排放、农业面源污染、工业废水的排放等。为有效控制水体污染，需明确水体所受污染源类型及其特性，并采取相应管控和治理措施，做到快速排查并控制污染源，确保水体水质安全。
通常，水质监测中依靠化学需氧量(COD)、生物需氧量(BOD)、氨氮(NH3-N)等水质参数来反映水质优劣，但这些参数所反映的是水体有机物整体的综合信息，很难判定其单一或特定的污染源。正定矩阵因子分解模型(PMF)是一种多元统计分析工具，通常将样本数据矩阵分解为因子贡献和因子分布两个矩阵。与一般因子分析模型相同的是，PMF模型根据受体参数的数据集分配源和源的贡献率；不同的是该方法可以得到非负约束的因子得分，并对数据的不确定性进行加权，当出现缺失值时可用数据的平均值或中位数替代，降低缺失数据对分析结果的影响。在PMF模型运行时，需人为确定因子数，不同因子数的选择会导致源解析过程中各参数的载荷分配具有不确定性，可能会对源解析的结果产生较大的干扰，以往的研究中往往依据区域特征及不同因子解析结果选择合适的因子数，但具有一定的主观性。同时，水体中溶解性荧光类有机物(CDOM)保存了水体中有机污染物的详细特征，如三维荧光光谱能根据污染物与含量不同，具有与水样一一对应的特点，如人的指纹具有唯一性，被称为污水的“荧光指纹”。荧光图谱解析一般利用识峰法，利用最大荧光强度所在位置，确定激发与发射波长，继而确定荧光峰类型，操作较为麻烦，并且因为荧光峰位置的变化，使样品荧光峰强度无法在同一维度进行比较。若样品过多会造成工作量巨大，不利于高效快捷判断水体污染来源。PARAFAC法基于三线分解理论，采用交替最小二乘算法可实现荧光物质的“数学分离”，但其只能识别出特征荧光组分，无法定义各荧光组分准确来源。
发明内容
本发明的目的在于：提供一种结合水质参数和三维荧光光谱进行污染物溯源的方法，利用常规水质参数的解析以及三维荧光光谱的分析对水体污染物进行溯源。
为了实现上述发明目的，本发明提供了一种结合水质参数和三维荧光光谱进行污染物溯源的方法，包括如下步骤：
水质检测步骤：对水样进行常规水质参数检测，再利用PCA/FA分析法确定常规水质参数的主因子数，再由PMF模型以主因子数计算各个常规水质参数的因子贡献率，再根据因子贡献率确定主要污染源；
荧光检测步骤：对水样进行三维荧光光谱检测，再对检测获得的三维荧光光谱行PARAFAC法解析，获得各个特征荧光组分，并计算各个特征荧光组分间的荧光强度比值；
污染源类型识别步骤：根据主要污染源以及荧光强度比值判断该水样采集点的污染源类型；
污染源确定步骤：获取监测区域的各个水样采集点布设位置以及水文参数，结合各个水样采集点的布设位置、水文参数以及污染源类型确定污染物源头。
进一步地，在进行水样采集点布设时，需要根据监测区域的形状和范围大小确定采样点的布设数量以及采样点的位置，在进行各个采集点的水样采集时，需要用0.45μm玻璃纤维滤膜进行过滤，去除水样中的杂质。
进一步地常规水质参数包括水体浊度、溶解氧、化学需氧量、总磷、氨氮以及硝氮；若进行原水测定，则常规水质参数还包括叶绿素a。
进一步地，利用PCA/FA分析法确定常规水质参数的主因子数时，若特征值大于1，则确定为主因子。
进一步地，PMF模型根据因子贡献率确定主要污染源时，选取存在因子贡献率大于50％的主因子作为主要污染因子，根据主要污染因子确定主要污染源。
进一步地，在对水样进行三维荧光光谱检测时，以超纯水作为空白消除三维荧光光谱的瑞丽散射。
进一步地，在根据主要污染源以及荧光强度比值判断污染源类型时，若类色氨酸峰T1与类色氨酸峰T2的比值介于1.5～3，则判定为生活源污染；若类色氨酸峰T1与类色氨酸峰T2的比值介于4～5，则判定为石油化工类污染；若类色氨酸峰T1与类富里酸峰M1的比值大于3，则判定为农业面源污染。
进一步地，在利用PCA/FA分析法确定常规水质参数的主因子数前，需要对因子分析的适用性进行验证，利用Bartlett球度检验及KMO标准进行检验，若KMO值在(0.5,1)的范围内表明可以进行因子分析。
进一步地，PARAFAC法解析是根据各个荧光组分峰的位置来确定各个特征荧光组分类型的。
本发明的有益效果在于：本发明结合PMF因子解析结果与特征荧光组分荧光强度比值，可快速判断水体环境是否受到人为活动影响，并同时判断其污染来源类型，再结合相关区域水质监测布点情况及水文参数定向筛查污染源，最终通过对比相应污染源检测结果明确水体主要污染来源，具有操作简便、准确度高、重现性好等优点，可以快速、有效识别水体环境中的典型污染源，为相关污染的管控与治理提供理论指导。
附图说明
图1为本发明的方法流程图；
图2为本发明的PCA/FA碎石图；
图3为本发明的某河流水质PMF分析常规水质参数的因子贡献率示意图。
具体实施方式
如图1所示，本发明提供了一种结合水质参数和三维荧光光谱进行污染物溯源的方法，包括如下步骤：
水质检测步骤：对水样进行常规水质参数检测，再利用PCA/FA分析法确定常规水质参数的主因子数，再由PMF模型以主因子数计算各个常规水质参数的因子贡献率，再根据因子贡献率确定主要污染源；
荧光检测步骤：将水样取4mL于10mm比色皿中，采用三维荧光光谱分析仪，设置相关参数，对水样进行三维荧光光谱检测，再对检测获得的三维荧光光谱进行PARAFAC法解析，PARAFAC法解析由MATLAB软件加载DOMflour工具包进行，获得特征荧光组分，同时根据各荧光组分峰的位置判断荧光组分类型，如类蛋白质、类富里酸、类腐殖质等。若荧光组分中类蛋白质比例超过70％，则认定该位点受到人为影响较严重，并计算各个特征荧光组分类型间的荧光强度比值；
污染源类型识别步骤：根据主要污染源以及荧光强度比值判断该水样采集点的污染源类型；
污染源确定步骤：获取监测区域的各个水样采集点布设位置以及水文参数，结合各个水样采集点的布设位置、水文参数以及污染源类型确定污染物源头。
进一步地，在进行水样采集点布设时，需要根据监测区域的形状和范围大小确定采样点的布设数量以及采样点的位置，在进行各个采集点的水样采集时，需要用0.45μm玻璃纤维滤膜进行过滤，去除水样中的杂质。
进一步地，常规水质参数包括水体浊度、溶解氧、化学需氧量、总磷、氨氮以及硝氮；若进行原水测定，则常规水质参数还包括叶绿素a。
进一步地，利用PCA/FA分析法确定常规水质参数的主因子数时，若特征值大于1，则确定为主因子。
进一步地，PMF模型根据因子贡献率确定主要污染源时，选取存在因子贡献率大于50％的主因子作为主要污染因子，根据主要污染因子确定主要污染源。
进一步地，在对水样进行三维荧光光谱检测时，以超纯水作为空白消除三维荧光光谱的瑞丽散射。
进一步地，在根据主要污染源以及荧光强度比值判断污染源类型时，若类色氨酸峰T1与类色氨酸峰T2的比值介于1.5～3，则判定为生活源污染；若类色氨酸峰T1与类色氨酸峰T2的比值介于4～5，则判定为石油化工类污染；若类色氨酸峰T1与类富里酸峰M1的比值大于3，则判定为农业面源污染。
进一步地，在利用PCA/FA分析法确定常规水质参数的主因子数前，需要对因子分析的适用性进行验证，利用Bartlett球度检验及KMO标准进行检验，若KMO值在(0.5,1)的范围内表明可以进行因子分析。
进一步地，PARAFAC法解析是根据各个荧光组分峰的位置来确定各个特征荧光组分类型的。
进一步地，PCA/FA分析法是将单个变量的数据变成均值为0，方差为1的数据集，具体公式为：
式(1)中，式中，Zi表示标准化数据，Ci表示样本变量值，Ai表示样本变量的平均浓度，Si表示变量i的标准偏差。PCA/FA在统计软件SPSS 22.0中完成，利用最大方差法(Varimax)旋转PCs，使得每个公因子中高载荷的变量数目达到最小，并保留特征值大于1的VFs作为主因子用于确定主因子数。
进一步地，PMF模型为：
式(2)中，Xnm表示原始矩阵(n×m)，代表n个样本，m表示监测参数；分解成两个矩阵Gnp(n×p)和Fpm(p×m)，p表示计算出的源的数量(提取因子的个数)；G表示源贡献率矩阵；F表示源成分谱矩阵(因子负荷)；Enm(n×m)为残差矩阵，代表解析结果与实测值的差值，表示为：
式(3)中，eij、xij、gik和fkj分别是矩阵E、X、G、F中对应的元素，PMF5.0对每个数据进行不确定性加权。定义一个目标函数Q(E)：
式(4)中，sij表示第i个样本中第j个监测参数的不确定性，其他字母定义同前。约束条件是矩阵G和F中所以元素非负限制。Q(E)作为PMF模型的判据之一，利用最小二乘法进行迭代计算，Q(E)收敛后才能进一步分析，多次运行不断地将原始矩阵X分解为非负矩阵G、F，选择Q(E)值较小的继续分析。不确定性的计算方法如下：
当xij≤MDL，则
当xij＞MDL，则
式(5)和(6)中，MDL表示监测参数的方法检测限；RSD表示对应监测参数的相对标准偏差。另外，sij可直接定义为RSD。PMF模型通过不断更改不确定度定义来确定源数目、最优成分谱和源贡献率。
进一步地，PARAFAC算法是基于三线分解理论采用交替最小二乘算法实现的一种数学模型，可用PARAFAC法实现样品中荧光物质的“数学分离”。测量得到的荧光光谱数据可组成一个I×J×K的矩阵，其中为样品数I，J，K分别为样品的发射波长和激发波长的数目。用平行因子法将其分解为得分矩阵A和负载矩阵B和C，分解模型表示为：
式(7)中，Xijk第i个样品、发射波长为j、激发波长为k处的荧光强度值；ain，bjn和ckn分别为负载矩阵A，B，C中的元素；e为残差矩阵E(I×J×K)的组成元素；N代表因子数即成分数。利用平行因子分析法对所有样品进行模拟计算，可解析出主要荧光组分类型。
以类蛋白质为例，荧光组分类蛋白质比例及荧光强度比值利用相对荧光强度进行计算，相对荧光强度计算方法为：
I＝Fmax×Ex(λmax)×Em(λmax) (8)
式(8)中，Fmax代表第n种成份的相对荧光强度，Ex(λmax)代表第n种成份激发负载的最大值，Em(λmax)代表第n种成份发射负载的最大值。荧光组分类蛋白质比例的计算方法为：
P％＝∑IP/∑I (9)
式(9)中，∑Ip为荧光组分类蛋白质含量，∑I为全部荧光组分。
本发明在具体实施验证时，在某特征河流共设有5个采样位点，流经区域基本为城乡结合区域或农业面源区域，易受到生活源及农业面源影响，同时，该河流附近基本无工业区，因此受石油化工等工业源影响应较小。
利用本发明专利的方法对该条特征河流污染来源进行溯源解析：
采用PCA/FA对进行标准化后的常规水质参数进行因子分析。为验证数据对因子分析的适用性，采用Bartlett球度检验及Kaiser-Mryer-Olkin(KMO)标准对数据进行检验。KMO值越接近于1，意味着变量间的相关性越强，原有变量间的相关性越强，一般KMO>0.5表明适合做因子分析，在本实例中KMO＝0.608。由图2可以看到，特征值大于1的因子数共有5个，因此选择PMF模型运算的因子数为5。
采用PMF模型对该特征河流水质进行各常规水质参数的因子贡献率解析，共解析得到5个因子(S1-S5)，如图3所示。
其中，因子S1对氮营养盐的贡献率较高，对NO3 --N和NH3-N的贡献率分别为66.4％和75.1％，即表示受到农业面源污染。因子S2则主要为含氮、磷有机物的污染，对NO2 --N和PO4 3-的贡献率较高，分别为76.4％和57.2％。因子S3主要来源于SS的贡献率，贡献率达到了70.4％，可能为区域降雨导致大量SS进入河流。T、DO、EC、Cl-和SO4 2-在S4表现出的源贡献率相对集中，介于40.8-47.5％，由于贡献率均未超过50％，说明该源对河流水体的贡献不显著。因子S5对Chla的贡献率最大，高达81.0％，表明浮游植物的生物量较高，故因子S5可代表浮游植物的影响。结果表明，该河流污染源分为：农业面源；氮、磷所带来的有机物污染，主要来源基本为生活源；河流降雨的影响，可视为自然来源；浮游植物的影响，可视为生活源。
将校正后的三维荧光光谱导入MATLAB通过DOMflour工具包进行PARAFAC分析，对各荧光组分的识别及相对荧光强度进行计算。经残差及折半分析检验最终确定为3个组分，其特征荧光峰位置分别为组分1：类富里酸荧光峰M1(245/411)，组分2：类蛋白质荧光峰T1(230/349)以及组分3：类蛋白质荧光峰T2(225，275/321)。将模型计算结果导出，通过发明内容中的计算步骤得到各点位的平均相对荧光强度，如表1所示。
表1某河流不同点位各荧光组分相对荧光强度及比例
如表1所示，在该特征河流各点位的T1+T2占总相对荧光强度的比例均达到70％以上，说明受人为影响严重。PMF模型分析显示，除自然来源外，该特征河流受到农业面源及生活源的影响显著，而无明显工业废水污染源因子。同时，T1/T2均在1.5～3，说明受生活源影响显著，而受石油化工类别等工业污染不显著。该结果与其地理位置特征相对应，表明溯源解析结果可靠。同时，采样点1、2、5的T1/M1>3，说明其受到生活源影响的同时，还受到一定程度农业面源的影响。结合该河流所处地理位置特征，发现采样点1、2、5均位于农田及乡间附近，易受到农业面源影响，而采样点3、4位于城区，受到农业面源污染的影响较小。
综合PMF模型及CDOM三维荧光光谱分析，得出该河流采样点1、2和5河段主要受到生活源和农业面源影响，而采样点3和4则主要来源于生活源的影响。综合水质参数及相关标准，若区域位点水质较差，且相关指标未能达标，则应根据污染源类型进一步查询相应污染源。从而，可以通过本专利方法的实施，实现在区域水体污染物溯源解析中有效缩小核查范围，快速有效确定污染源，并控制其造成的环境污染。
Claims (9)
1.一种结合水质参数和三维荧光光谱进行污染物溯源的方法，其特征在于，包括如下步骤：
水质检测步骤：对水样进行常规水质参数检测，再利用PCA/FA分析法确定常规水质参数的主因子数，再由PMF模型以主因子数计算各个常规水质参数的因子贡献率，再根据因子贡献率确定主要污染源；
荧光检测步骤：对水样进行三维荧光光谱检测，再对检测获得的三维荧光光谱行PARAFAC法解析，获得各个特征荧光组分，并计算各个特征荧光组分间的荧光强度比值；
污染源类型识别步骤：根据主要污染源以及荧光强度比值判断该水样采集点的污染源类型；
污染源确定步骤：获取监测区域的各个水样采集点布设位置以及水文参数，结合各个水样采集点的布设位置、水文参数以及污染源类型确定污染物源头。
2.根据权利要求1所述的结合水质参数和三维荧光光谱进行污染物溯源的方法，其特征在于，在进行水样采集点布设时，需要根据监测区域的形状和范围大小确定采样点的布设数量以及采样点的位置，在进行各个采集点的水样采集时，需要用0.45μm玻璃纤维滤膜进行过滤，去除水样中的杂质。
3.根据权利要求1所述的结合水质参数和三维荧光光谱进行污染物溯源的方法，其特征在于，常规水质参数包括水体浊度、溶解氧、化学需氧量、总磷、氨氮以及硝氮；若进行原水测定，则常规水质参数还包括叶绿素a。
4.根据权利要求1所述的结合水质参数和三维荧光光谱进行污染物溯源的方法，其特征在于，利用PCA/FA分析法确定常规水质参数的主因子数时，若特征值大于1，则确定为主因子。
5.根据权利要求1所述的结合水质参数和三维荧光光谱进行污染物溯源的方法，其特征在于，PMF模型根据因子贡献率确定主要污染源时，选取存在因子贡献率大于50％的主因子作为主要污染因子，根据主要污染因子确定主要污染源。
6.根据权利要求1所述的结合水质参数和三维荧光光谱进行污染物溯源的方法，其特征在于，在对水样进行三维荧光光谱检测时，以超纯水作为空白消除三维荧光光谱的瑞丽散射。
7.根据权利要求1所述的结合水质参数和三维荧光光谱进行污染物溯源的方法，其特征在于，在根据主要污染源以及荧光强度比值判断污染源类型时，若类色氨酸峰T1与类色氨酸峰T2的比值介于1.5～3，则判定为生活源污染；若类色氨酸峰T1与类色氨酸峰T2的比值介于4～5，则判定为石油化工类污染；若类色氨酸峰T1与类富里酸峰M1的比值大于3，则判定为农业面源污染。
8.根据权利要求1所述的结合水质参数和三维荧光光谱进行污染物溯源的方法，其特征在于，在利用PCA/FA分析法确定常规水质参数的主因子数前，需要对因子分析的适用性进行验证，利用Bartlett球度检验及KMO标准进行检验，若KMO值在(0.5,1)的范围内表明可以进行因子分析。
9.根据权利要求1所述的结合水质参数和三维荧光光谱进行污染物溯源的方法，其特征在于，PARAFAC法是根据各荧光组分峰的位置来确定各个特征荧光组分类型的。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811355655.5A CN109470667A (zh) | 2018-11-14 | 2018-11-14 | 一种结合水质参数和三维荧光光谱进行污染物溯源的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811355655.5A CN109470667A (zh) | 2018-11-14 | 2018-11-14 | 一种结合水质参数和三维荧光光谱进行污染物溯源的方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109470667A true CN109470667A (zh) | 2019-03-15 |
Family
ID=65672777
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811355655.5A Pending CN109470667A (zh) | 2018-11-14 | 2018-11-14 | 一种结合水质参数和三维荧光光谱进行污染物溯源的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109470667A (zh) |
Cited By (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110390494A (zh) * | 2019-08-13 | 2019-10-29 | 成都理工大学 | 农业区生活垃圾填埋场地下水中“三氮”的溯源方法 |
| CN110531043A (zh) * | 2019-08-29 | 2019-12-03 | 中水三立数据技术股份有限公司 | 一种流域水环境污染智能管理系统及其方法 |
| CN111289486A (zh) * | 2020-04-01 | 2020-06-16 | 西北农林科技大学 | 一种利用溶解性有机物荧光特征定量测定点源污染和面源污染贡献率的方法 |
| CN111426668A (zh) * | 2020-04-28 | 2020-07-17 | 华夏安健物联科技(青岛)有限公司 | 一种利用三维荧光光谱特征信息对污染水体溯源分类识别的方法 |
| CN111537490A (zh) * | 2020-07-08 | 2020-08-14 | 佛山求识光谱数据科技有限公司 | 一种基于冰箱的肉类抗生素检测方法 |
| CN111861421A (zh) * | 2020-08-01 | 2020-10-30 | 生态环境部南京环境科学研究所 | 一种流域突发水污染快速溯源方法 |
| CN112198144A (zh) * | 2020-09-16 | 2021-01-08 | 安徽泽众安全科技有限公司 | 一种快速污水溯源的方法及系统 |
| CN112630202A (zh) * | 2020-12-10 | 2021-04-09 | 湖南大学 | 一种城市排水系统的雨天溢流污水来源的识别方法 |
| CN112903660A (zh) * | 2021-03-11 | 2021-06-04 | 广西大学 | 一种判别流域水体污染现状和污染来源的方法 |
| CN113011478A (zh) * | 2021-03-05 | 2021-06-22 | 清华大学 | 基于数据融合的污染源识别方法及系统 |
| CN113567401A (zh) * | 2020-04-28 | 2021-10-29 | 中国环境科学研究院 | 一种垃圾填埋场渗滤液污染地下水情况的快速检测方法和应用 |
| CN113916847A (zh) * | 2021-07-20 | 2022-01-11 | 江苏省扬州环境监测中心 | 一种基于光谱技术和线性支持向量算法的水质检测方法 |
| CN114965348A (zh) * | 2022-07-27 | 2022-08-30 | 浙江数翰科技有限公司 | 基于污水检测的光谱解析方法和系统 |
| CN115219472A (zh) * | 2022-08-12 | 2022-10-21 | 生态环境部华南环境科学研究所(生态环境部生态环境应急研究所) | 一种定量识别混合水体多污染源的方法及系统 |
| CN115266667A (zh) * | 2022-07-21 | 2022-11-01 | 杭州石炭纪环保科技有限公司 | 一种基于三维荧光光谱的江河流域污染来源确定方法及系统 |
| WO2022241883A1 (zh) * | 2021-05-17 | 2022-11-24 | 清华大学 | 基于三维液相色谱指纹的污染源识别方法及装置 |
| CN115825024A (zh) * | 2022-10-10 | 2023-03-21 | 中国农业大学烟台研究院 | 一种水质现场快速检测方法 |
| CN115983666A (zh) * | 2022-11-23 | 2023-04-18 | 生态环境部华南环境科学研究所(生态环境部生态环境应急研究所) | 一种评估水体来源和贡献的荧光指标 |
| CN116628598A (zh) * | 2023-05-15 | 2023-08-22 | 生态环境部华南环境科学研究所(生态环境部生态环境应急研究所) | 一种基于大数据和nmf模型的二噁英来源解析方法及系统 |
| CN117634178A (zh) * | 2023-11-24 | 2024-03-01 | 广东省环境科学研究院 | 一种高效溯源水质污染源的系统及方法 |
| CN118071370A (zh) * | 2024-02-07 | 2024-05-24 | 长江勘测规划设计研究有限责任公司 | 排水管网排口出流水体中溶解性有机物的溯源方法及系统 |
| CN118883880A (zh) * | 2024-10-08 | 2024-11-01 | 安徽新宇环保科技股份有限公司 | 一种基于多维数据融合分析的水污染溯源分析方法及系统 |
| CN119539994A (zh) * | 2025-01-23 | 2025-02-28 | 北京大学 | 地下水污染源解析方法、系统、设备、介质及产品 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101576485A (zh) * | 2009-06-04 | 2009-11-11 | 浙江大学 | 一种多源光谱融合水质分析方法 |
| CN102621118A (zh) * | 2012-03-18 | 2012-08-01 | 吉林大学 | 温室蔬菜病虫害的预警方法 |
| CN105532524A (zh) * | 2015-12-10 | 2016-05-04 | 浙江海洋学院 | 一种日本囊对虾虾苗耐氨氮性能的等级划分与优选方法 |
| CN105699345A (zh) * | 2016-01-25 | 2016-06-22 | 耿春茂 | 一种三维荧光光谱结合parafac算法测定污染物的方法 |
| CN106018359A (zh) * | 2016-05-12 | 2016-10-12 | 常州市排水管理处 | 一种污水厂水质监测预警方法及系统 |
| CN108052486A (zh) * | 2017-12-01 | 2018-05-18 | 南开大学 | 一种基于颗粒物无机组分和有机标识物的精细化源解析方法 |
-
2018
- 2018-11-14 CN CN201811355655.5A patent/CN109470667A/zh active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101576485A (zh) * | 2009-06-04 | 2009-11-11 | 浙江大学 | 一种多源光谱融合水质分析方法 |
| CN102621118A (zh) * | 2012-03-18 | 2012-08-01 | 吉林大学 | 温室蔬菜病虫害的预警方法 |
| CN105532524A (zh) * | 2015-12-10 | 2016-05-04 | 浙江海洋学院 | 一种日本囊对虾虾苗耐氨氮性能的等级划分与优选方法 |
| CN105699345A (zh) * | 2016-01-25 | 2016-06-22 | 耿春茂 | 一种三维荧光光谱结合parafac算法测定污染物的方法 |
| CN106018359A (zh) * | 2016-05-12 | 2016-10-12 | 常州市排水管理处 | 一种污水厂水质监测预警方法及系统 |
| CN108052486A (zh) * | 2017-12-01 | 2018-05-18 | 南开大学 | 一种基于颗粒物无机组分和有机标识物的精细化源解析方法 |
Non-Patent Citations (5)
| Title |
|---|
| COLIN A. STEDMON，ET.AL.: ""Tracing dissolved organic matter in aquatic environments using a new approach to fluorescence spectroscopy"", 《MARINE CHEMISTRY》 * |
| 祝鹏 等: ""平行因子分析法在太湖水体三维荧光峰比值分析中的应用"", 《光谱学与光谱分析》 * |
| 祝鹏 等: ""湖泊水体三维荧光光谱的PARAFAC法在污染源解析中的应用"", 《光谱学与光谱分析》 * |
| 郭田田: ""流域尺度水质时空变异特征及污染源解析研究"", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
| 陈茂福 等: ""城市污水的三维荧光指纹特征"", 《光学学报》 * |
Cited By (33)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110390494A (zh) * | 2019-08-13 | 2019-10-29 | 成都理工大学 | 农业区生活垃圾填埋场地下水中“三氮”的溯源方法 |
| CN110390494B (zh) * | 2019-08-13 | 2022-04-26 | 成都理工大学 | 农业区生活垃圾填埋场地下水中“三氮”的溯源方法 |
| CN110531043A (zh) * | 2019-08-29 | 2019-12-03 | 中水三立数据技术股份有限公司 | 一种流域水环境污染智能管理系统及其方法 |
| CN111289486A (zh) * | 2020-04-01 | 2020-06-16 | 西北农林科技大学 | 一种利用溶解性有机物荧光特征定量测定点源污染和面源污染贡献率的方法 |
| CN111289486B (zh) * | 2020-04-01 | 2022-11-01 | 西北农林科技大学 | 一种利用溶解性有机物荧光特征定量测定点源污染和面源污染贡献率的方法 |
| CN111426668B (zh) * | 2020-04-28 | 2023-09-08 | 华夏安健物联科技(青岛)有限公司 | 一种利用三维荧光光谱特征信息对污染水体溯源分类识别的方法 |
| CN111426668A (zh) * | 2020-04-28 | 2020-07-17 | 华夏安健物联科技(青岛)有限公司 | 一种利用三维荧光光谱特征信息对污染水体溯源分类识别的方法 |
| CN113567401B (zh) * | 2020-04-28 | 2022-09-30 | 中国环境科学研究院 | 一种垃圾填埋场渗滤液污染地下水情况的快速检测方法和应用 |
| CN113567401A (zh) * | 2020-04-28 | 2021-10-29 | 中国环境科学研究院 | 一种垃圾填埋场渗滤液污染地下水情况的快速检测方法和应用 |
| CN111537490A (zh) * | 2020-07-08 | 2020-08-14 | 佛山求识光谱数据科技有限公司 | 一种基于冰箱的肉类抗生素检测方法 |
| CN111861421A (zh) * | 2020-08-01 | 2020-10-30 | 生态环境部南京环境科学研究所 | 一种流域突发水污染快速溯源方法 |
| CN112198144A (zh) * | 2020-09-16 | 2021-01-08 | 安徽泽众安全科技有限公司 | 一种快速污水溯源的方法及系统 |
| CN112630202A (zh) * | 2020-12-10 | 2021-04-09 | 湖南大学 | 一种城市排水系统的雨天溢流污水来源的识别方法 |
| CN113011478A (zh) * | 2021-03-05 | 2021-06-22 | 清华大学 | 基于数据融合的污染源识别方法及系统 |
| CN112903660A (zh) * | 2021-03-11 | 2021-06-04 | 广西大学 | 一种判别流域水体污染现状和污染来源的方法 |
| WO2022241883A1 (zh) * | 2021-05-17 | 2022-11-24 | 清华大学 | 基于三维液相色谱指纹的污染源识别方法及装置 |
| CN113916847A (zh) * | 2021-07-20 | 2022-01-11 | 江苏省扬州环境监测中心 | 一种基于光谱技术和线性支持向量算法的水质检测方法 |
| CN113916847B (zh) * | 2021-07-20 | 2024-04-09 | 江苏省扬州环境监测中心 | 一种基于光谱技术和线性支持向量算法的水质检测方法 |
| CN115266667A (zh) * | 2022-07-21 | 2022-11-01 | 杭州石炭纪环保科技有限公司 | 一种基于三维荧光光谱的江河流域污染来源确定方法及系统 |
| CN114965348A (zh) * | 2022-07-27 | 2022-08-30 | 浙江数翰科技有限公司 | 基于污水检测的光谱解析方法和系统 |
| CN115219472A (zh) * | 2022-08-12 | 2022-10-21 | 生态环境部华南环境科学研究所(生态环境部生态环境应急研究所) | 一种定量识别混合水体多污染源的方法及系统 |
| CN115825024B (zh) * | 2022-10-10 | 2024-05-24 | 中国农业大学烟台研究院 | 一种水质现场快速检测方法 |
| CN115825024A (zh) * | 2022-10-10 | 2023-03-21 | 中国农业大学烟台研究院 | 一种水质现场快速检测方法 |
| CN115983666B (zh) * | 2022-11-23 | 2023-09-22 | 生态环境部华南环境科学研究所(生态环境部生态环境应急研究所) | 一种评估水体来源和贡献的荧光指标 |
| CN115983666A (zh) * | 2022-11-23 | 2023-04-18 | 生态环境部华南环境科学研究所(生态环境部生态环境应急研究所) | 一种评估水体来源和贡献的荧光指标 |
| CN116628598A (zh) * | 2023-05-15 | 2023-08-22 | 生态环境部华南环境科学研究所(生态环境部生态环境应急研究所) | 一种基于大数据和nmf模型的二噁英来源解析方法及系统 |
| CN116628598B (zh) * | 2023-05-15 | 2024-03-12 | 生态环境部华南环境科学研究所(生态环境部生态环境应急研究所) | 一种基于大数据和nmf模型的二噁英来源解析方法及系统 |
| CN117634178A (zh) * | 2023-11-24 | 2024-03-01 | 广东省环境科学研究院 | 一种高效溯源水质污染源的系统及方法 |
| CN117634178B (zh) * | 2023-11-24 | 2024-05-28 | 广东省环境科学研究院 | 一种高效溯源水质污染源的系统及方法 |
| CN118071370A (zh) * | 2024-02-07 | 2024-05-24 | 长江勘测规划设计研究有限责任公司 | 排水管网排口出流水体中溶解性有机物的溯源方法及系统 |
| CN118883880A (zh) * | 2024-10-08 | 2024-11-01 | 安徽新宇环保科技股份有限公司 | 一种基于多维数据融合分析的水污染溯源分析方法及系统 |
| CN118883880B (zh) * | 2024-10-08 | 2025-01-07 | 安徽新宇环保科技股份有限公司 | 一种基于多维数据融合分析的水污染溯源分析方法及系统 |
| CN119539994A (zh) * | 2025-01-23 | 2025-02-28 | 北京大学 | 地下水污染源解析方法、系统、设备、介质及产品 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109470667A (zh) | 一种结合水质参数和三维荧光光谱进行污染物溯源的方法 | |
| Demattê et al. | Soil analytical quality control by traditional and spectroscopy techniques: Constructing the future of a hybrid laboratory for low environmental impact | |
| CN114295749B (zh) | 一种水体有机污染智能化溯源方法及系统 | |
| CN105158353B (zh) | 土壤中多环芳烃污染源解析的方法 | |
| CN113916847B (zh) | 一种基于光谱技术和线性支持向量算法的水质检测方法 | |
| Su et al. | Temporal trend and source apportionment of water pollution in different functional zones of Qiantang River, China | |
| CN112735539B (zh) | 基于复合受体模型的污染场地特征污染物的源解析方法 | |
| CN117195135B (zh) | 一种水污染异常溯源检测方法及系统 | |
| CN106650020A (zh) | 一种复合受体模型污染源解析方法 | |
| CN109470668A (zh) | 一种结合氮氧同位素和三维荧光进行污染物溯源的方法 | |
| CN118275636A (zh) | 利用水系沉积物解析地下水重金属污染源及污染边界的方法 | |
| CN116297936A (zh) | 一种污染物的溯源方法、系统、装置及介质 | |
| CN113971989A (zh) | 一种基于opls的森林土壤有机碳含量高光谱建模方法 | |
| CN117974404A (zh) | 一种水陆协同的陆域污染源分析方法及系统 | |
| CN113283072A (zh) | 一种适用于多场景条件的水体cod检测方法 | |
| CN115389730A (zh) | 污染水域智能识别治理方法、系统及存储介质 | |
| CN113533672A (zh) | 一种水质在线监测方法及装置 | |
| CN118858564A (zh) | 一种地下水污染溯源方法 | |
| CN108663334A (zh) | 基于多分类器融合寻找土壤养分光谱特征波长的方法 | |
| NL2034211B1 (en) | Method and system for quantitatively identifying multi-pollution sources of mixed water body | |
| Lancellotti et al. | Accuracy of rapid tests used for analysis of advanced onsite wastewater treatment system effluent | |
| CN118111958A (zh) | 基于三维荧光指纹的排水管网清水混接点识别方法及系统 | |
| CN113376114A (zh) | 一种基于紫外-可见光谱数据的水污染溯源方法 | |
| CN112966783A (zh) | 一种矿井水常量离子特征标识及水源判别的方法 | |
| CN108491995B (zh) | 一种用于饮用水风险因子识别的关键控制因素筛选方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190315 |