基于主成分分析法的浑河抚顺段多环芳烃来源分

1、基于主成分分析法的浑河抚顺段多环芳烃来源分析摘要：随着城镇化程度的提高，给生态环境带来负面影响，多环芳烃（PAHs）由于难降解等特性存在于河流、土壤环境中，对居民健康造成一定威胁。本文选取浑河抚顺段与抚新城为研究区域，利用同分异构体比值法与主成分分析法对城镇化进程中多环芳烃来源进行解析，发现河流水体干湿两季中PAHs皆来源于石油燃烧与煤、柴等生物质燃烧，湿季中PAHs来源存在部分石油输入；沈抚新城干湿两季表层土壤中PAHs来源皆包括石油泄露与大气沉降混合源、石油泄露与柴油燃烧混合源以及煤燃烧与机动车排放混合源，差异在于湿季PAHs来源包括柴油燃烧源；对比两种方法发现主成分分析法得到的结果更具体

2、，为相关研究提供了理论依据。关键词：PAHs；同分异构体比值法；主成分分析法；干湿两季Source Analysis of Polycyclic Aromatic Hydrocarbons in Fushun Section of Hunhe River Based on Principal Component AnalysisAbstract: As the degree of urbanization increases, it has a negative impact on the ecological environment. Polycyclic aromatic hydrocar

3、bons (PAHs), which are hard to degrade, exist in the rivers and soil environment, posing a certain threat to the residents health. In this paper, the Fushun section of the Fushun River and the new city of Fushun were selected as the study area. The PAHs were analyzed by the method of isomer ratio an

4、d principal component analysis. The results showed that PAHs in both wet and dry river water sources were all from The PAHs sources in the wet season were partly oil input. The sources of PAHs in surface soil of both wet and dry seasons in Shen-Fu new city included both the mixture source of oil lea

5、kage and atmospheric deposition, the mixture source of oil leakage and diesel combustion and the mixed sources of coal combustion and vehicle emissions. The difference lies in the sources of PAHs in wet season including diesel combustion sources. Comparing the two methods, the results obtained by pr

6、incipal component analysis are more specific and provide the theoretical basis for the related research.Key words: PAHs; isomer ratio method; principal component analysis; dry and wet seasons0 引言随着我国城镇化程度的提高，一些污染物不断累积在河流、土壤等环境中，给居民健康带来威胁。由于多环芳烃污染物具有难挥发、难降解的特性引起社会的广泛关注1。环境中PAHs主要来源包括人为源与天然源，人为源包括某些化工

7、产品的生成与使用、未经燃烧的石油类产品与各种燃料的燃烧，天然源包括火山爆发、天然火灾以及生物、微生物合成等自然活动2，而天然源所占比例很小，环境中PAHs主要来源于人为源，因此展开城镇化进程中多环芳烃来源解析研究具有十分重要的意义。本文选取浑河抚顺段与抚新城为研究区域，利用同分异构体比值法与主成分分析法对城镇化进程中多环芳烃来源进行解析，为相关研究提供了理论依据。1 PAHs来源分析方法1.1 同分异构体比值法互为同分异构体的多环芳烃通常具有相似的动力学质量转移系数与热力学分配系数，因此能够根据其含量比值来区分不同热源与石油源等，PAHs来源可以通过lnP/(lnP+BghiP)、lnP/(l

8、nP+BaP)、BaA/(BaA+Chry) 和Fla/(Fla+Pyr)等比值法进行分析，其标准为：lnP/(lnP+BghiP)0.2、lnP/(lnP+BaP)0.2、BaA/(BaA+Chry)0.2、Fla/(Fla+Pyr)0.5、lnP/(lnP+BaP)0.5、BaA/(BaA+Chry)0.35、Fla/(Fla+Pyr)0.5，PAHs来源为草、木、煤燃烧源；0.2lnP/(lnP+BghiP)0.5、0.2lnP/(lnP+BaP)0.5、0.2BaA/(BaA+Chry)0.35、0.4Fla/(Fla+Pyr)0.5，PAHs来源为石油燃烧源3。1.2 主成分分析法尽

9、管同分异构体比值法能够对PAHs来源进行初步判断，但是容易受大气沉降、生物降解等因素的影响，得到的结论准确性需要进一步完善3.73,88。利用主成分分析“降维”的思路，根据SPSS软件进行主成分分析等统计学方法分析PAHs来源，能够得到反映原始数据的综合变量，依据方差进行排列，方差最大代表最多的数据信息，依次类推，最后根据各个主成分PAHs的载荷量来分析污染物的来源89。通常，NaP、Ace和Acy等PAHs代表石油挥发或者石油泄露；Phe、Fla、Ant、Chry、BaA、Pyr等表示煤燃烧的特征指示物；BaP、dBaAnt和BghiP代表汽油燃烧；BkF、BbF和InP表示柴油燃烧的特征指

10、示物4,5。2 城镇化进程中多环芳烃来源解析2.1 利用同分异构体比值法进行解析2.1.1河流多环芳烃来源分析由于河流中Chry与BaA不易检测，因此利用lnP/(lnP+BghiP)、lnP/(lnP+BaP)、Fla/(Fla+Pyr)对河流样品进行同分异构体来源分析，结果见图1。图1 河流水体干湿两季PAHs同分异构体比值由图1可知，河流水体干季中PAHs主要来源于石油燃烧与煤、柴等生物质燃烧，湿季中PAHs同样来源于石油燃烧与草、木、煤燃烧，同时存在部分石油输入源，两季来源差异体现在湿季来源多了石油输入。浑河抚顺段处于我国严寒地区，河流结冰现象普遍，环境中PAHs很难迁移至水中，因此干

11、季水体中累积的PAHs来自燃烧源；由于湿季温度高，PAHs能够伴随大气沉降进行迁移进入水体。我国城市依河而建，河流中PAHs多数来源于人类活动，因此两季来源差异与季节变化相关。2.1.2 表层土壤多环芳烃来源分析由于研究区干、湿两季中Fla、Chry、BaA、Pyr、BaP、BghiP、Inp等PAHs单体检出率较高，所以使用lnP/(lnP+BaP)、lnP/(lnP+BghiP)、Fla/(Fla+Pyr)和BaA/(BaA+Chry)四种同分异构体比值对表层土壤PAHs进行分析，结果见图2。图2沈抚新城干湿两季PAHs同分异构体比值由图2可知，沈抚新城干季土壤中PAHs主要来源于石油燃烧

12、与草、木、煤等生物质燃烧，部分来源于石油输入，湿季土壤中PAHs皆来源于石油燃烧与草、木、煤等生物质燃烧。干季采样时沈抚新城处于寒冬供暖期，草、木、煤等生物质燃烧向空气中排放大量PAHs，另外，沈抚新城处于城市转变过渡时期，一些现代化器械的使用必然会带来石油燃烧PAHs；湿季中不存在石油输入源PAHs，这是由于季节更迭，温度变化，使得石油输入源PAHs发生迁移，因此没有检测到石油输入源，而具体的PAHs来源与当地气候、居民生活方式有很大的关系。2.2利用主成分分析法进行解析2.2.1河流多环芳烃来源分析通过软件SPSS19.0结合方差极大正交旋转法进行河流多环芳烃成分分析，提取特征值根大于的主

13、因子，得到河流干湿两季水体主成分因子载荷，结果见表1。表1 河流干湿两季水体样品方差极大旋转后的主成分因子载荷河流水体PAHs干季湿季PC1PC2PC3PC1PC2PC3MNaP10.972-0.1360.092-0.3900.2960.831MNaP20.970-0.1610.062-0.3660.3350.804NaP0.1780.6140.740-0.3900.1670.796Ace0.980-0.1360.0080.3430.824-0.260Acy0.9660.0650.091-0.4110.3300.458Phe0.9320.2540.0870.8410.435-0.209Dibt

14、0.980-0.1070.0680.0680.639-0.224FI0.984-0.0110.0680.4160.82-0.055Fla0.2720.918-0.2040.9530.1420.156Ret0.2470.8060.5190.916-0.2260.165Ant0.970-0.1830.0450.7460.456-0.336Chry0.7570.566-0.1960.936-0.1510.199BaA0.957-0.2470.0270.9440.0720.163Pyr0.1920.866-0.4120.897-0.1480.269BaP0.9610.032-0.1770.890-0.

15、3280.259BkF0.978-0.145-0.0340.331-0.7440.020BbF0.6140.602-0.3250.943-0.038-0.071dBaAnt0.958-0.2440.0250.859-0.3350.124Pery0.969-0.198-0.0160.7670.1840.276BghiP0.981-0.086-0.0810.921-0.1740.210Inp0.981-0.127-0.0340.935-0.1830.210解释方差变量（%）73.38517.7186.03456.66616.76513.937累计方差贡献率（%）73.38591.10397.137

16、56.66673.52087.458由表1可知，河流干季水体PAHs的三个成分累积解释了方差变量的97.137%，其中成分解释了方差变量的73.385%，成分1的载荷较大指示物有FI、Ant、Phe、Chry、BaA、BaP、BkF、BbF、dBaAnt、BghiP和Inp，分别代表煤炭燃烧和机动车排放；成分2解释了方差变量的17.818%，贡献较大的为Fla与Pyr，划分为煤炭燃烧；成分解释了方差变量的6.034%，载荷最大为NaP，代表石油泄露。河流湿季水体PAHs的三个成分累积解释了方差变量的87.458%，其中成分解释了方差变量的56.666%，成分1的载荷较大指示物包括Phe、Fla

17、、Ret、Chry、BaA、BbF等，同样代表煤炭燃烧和机动车排放；成分2解释了方差变量的16.854%，载荷较大为Ace、Dibt和FI，划分为石油泄露；成分解释了方差变量的13.938%，载荷最大为MNaP1、MNaP2、NaP和Acy，划分为石油泄露与大气沉降。以上可知，河流水体多环芳烃来源分析结果与同分异构体分析结果一致，并且比同分异构体比值法更加透彻。相对浑河抚顺段来说，由于湿季中增加了大气沉降源，表明季节变化能够对水体中PAHs产生一定的影响；干季中普遍的煤炭供暖必然会成为多环芳烃的主要来源；干湿两季都包含石油泄露输入与机动车排放，这与浑河船舶航运以及两岸机动车行驶密不可分。2.2

18、.2 表层土壤多环芳烃来源分析利用同种方法对沈抚新城干、湿两季表层土壤PAHs进行主成分析，结果见表2。表2 沈抚新城干湿两季土壤样品方差极大旋转后的主成分因子载荷河流水体PAHs干季湿季PC1PC2PC3PC1PC2PC3MNaP10.3160.932-0.0560.5060.8450.054MNaP20.2410.9460.1410.4700.8570.103NaP-0.0030.8220.5360.0670.5540.468Ace0.9540.1210.0450.900-0.1900.308Acy0.931-0.025-0.0170.5820.743-0.080Phe0.9640.049

19、-0.0460.9760.1360.017Dibt0.9440.150-0.1550.904-0.2900.230F10.9570.146-0.0220.8840.3210.202Fla0.8670.098-0.2450.5660.668-0.281Ret0.3970.455-0.7510.4080.557-0.495Ant0.9780.020-0.0920.9450.1810.153Chry0.982-0.0650.0300.959-0.205-0.044BaA0.975-0.1160.0520.976-0.1570.031Pyr0.980-0.0770.0130.885-0.0830.29

20、1BaP0.972-0.1400.0100.951-0.242-0.050BkF0.873-0.2370.2220.941-0.256-0.101BbF0.981-0.0970.0330.947-0.249-0.100dBaAnt0.974-0.1290.0870.829-0.146-0.402Pery0.967-0.1440.1070.932-0.3140.010BghiP0.972-0.1460.0880.928-0.282-0.125Inp0.974-0.1310.0820.937-0.270-0.099解释方差变量（%）76.31013.9595.14268.32518.6235.15

21、8累计方差贡献率（%）76.31090.26994.42168.32586.94892.109由表2可知，沈抚新城干季表层土壤中PAHs的三个成分累积解释了方差变量的94.421%，其中成分解释了方差变量的76.310%，成分1的载荷较大指示物有Ace、Acy、Dibt、F1、Fla、Ant、Chry、BaA、Pyr、BaP、BbF、dBaAnt、Pery、BghiP、Inp，划分为机动车排放与煤炭燃烧；成分2解释 了方差变量的13.959%，贡献较大的为MNaP1、MNaP2、NaP，代表大气沉降与石油泄露；成分解释了方差变量的4.152%，载荷最大为NaP，代表石油泄露。抚新城湿季表层土壤

22、中PAHs的三个成分累积解释了方差变量的92.109%，其中成分解释了方差变量的68.325%，成分1的载荷较大指示物包括Ace、Phe、Dibt、Ant、Chry、BaA等，表明成分1代表机动车排放与煤炭燃烧；成分2解释了方差变量的18.623%，载荷较大为MNaP1、MNaP2、NaP与Fla，划分为柴油燃烧与石油泄露；成分解释了方差变量的5.161%，载荷最大为NaP、Ace、Dibt和F1，划分为石油泄露与大气沉降。综上可知，与同分异构体法相比，主成分分析法的分析结果更加具体化。沈抚新城干湿季表层土壤中PAHs来源存在一定的差异，湿季主要来源为石油泄露与大气沉降混合源、石油泄露与柴油燃

23、烧混合源以及煤燃烧与机动车排放混合源，干季主要来源为石油泄露输入源、石油泄露输入与大气沉降混合源和煤燃烧与机动车排放混合源，两季差异在于柴油燃烧源。3结论本文首先介绍了同分异构体比值法与主成分分析法，然后对浑河抚顺段与沈抚新城干湿两季河流与土壤中多环芳烃进行来源解析，得到以下主要结论：（1）利用同分异构体比值法进行城镇化进程中多环芳烃来源解析，发现河流水体干湿两季中PAHs皆来源于石油燃烧与煤、柴等生物质燃烧，而湿季中PAHs来源存在部分石油输入源；沈抚新城干湿两季土壤中PAHs来源于石油燃烧与草、木、煤等生物质燃烧，湿季并未发现石油输入源，原因是季节变化影响了石油输入源PAHs的迁移。（2）利用主成分分析法进行城镇化进程中多环芳烃来源解析，得到河流与土壤中主成分因子荷载。结果表明，浑河抚顺段湿季