CN117316326A - 一种基于监测数据及气象要素的VOCs污染来源解析方法 - Google Patents
一种基于监测数据及气象要素的VOCs污染来源解析方法 Download PDFInfo
- Publication number
- CN117316326A CN117316326A CN202311291677.0A CN202311291677A CN117316326A CN 117316326 A CN117316326 A CN 117316326A CN 202311291677 A CN202311291677 A CN 202311291677A CN 117316326 A CN117316326 A CN 117316326A
- Authority
- CN
- China
- Prior art keywords
- vocs
- substances
- concentration
- pollution
- analysis
- 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
- 239000012855 volatile organic compound Substances 0.000 title claims abstract description 122
- 238000012544 monitoring process Methods 0.000 title claims abstract description 51
- 238000004458 analytical method Methods 0.000 title claims abstract description 40
- CBENFWSGALASAD-UHFFFAOYSA-N Ozone Chemical compound [O-][O+]=O CBENFWSGALASAD-UHFFFAOYSA-N 0.000 claims abstract description 38
- 238000001514 detection method Methods 0.000 claims abstract description 35
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 19
- 231100000719 pollutant Toxicity 0.000 claims abstract description 19
- 238000010219 correlation analysis Methods 0.000 claims abstract description 13
- 239000003570 air Substances 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims abstract description 10
- 239000000809 air pollutant Substances 0.000 claims abstract description 9
- 239000012080 ambient air Substances 0.000 claims abstract description 7
- 238000009826 distribution Methods 0.000 claims abstract description 4
- 239000000126 substance Substances 0.000 claims description 99
- -1 aldehyde ketone Chemical class 0.000 claims description 39
- 238000000034 method Methods 0.000 claims description 27
- 239000013076 target substance Substances 0.000 claims description 25
- 238000011160 research Methods 0.000 claims description 16
- 150000001335 aliphatic alkanes Chemical class 0.000 claims description 10
- 150000001336 alkenes Chemical class 0.000 claims description 10
- 150000001345 alkine derivatives Chemical class 0.000 claims description 10
- 150000004945 aromatic hydrocarbons Chemical class 0.000 claims description 10
- UHOVQNZJYSORNB-UHFFFAOYSA-N Benzene Chemical compound C1=CC=CC=C1 UHOVQNZJYSORNB-UHFFFAOYSA-N 0.000 claims description 9
- YXFVVABEGXRONW-UHFFFAOYSA-N Toluene Chemical compound CC1=CC=CC=C1 YXFVVABEGXRONW-UHFFFAOYSA-N 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 9
- 238000002290 gas chromatography-mass spectrometry Methods 0.000 claims description 8
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Natural products C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 8
- 238000004811 liquid chromatography Methods 0.000 claims description 6
- 238000004519 manufacturing process Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 239000007789 gas Substances 0.000 claims description 5
- MGWGWNFMUOTEHG-UHFFFAOYSA-N 4-(3,5-dimethylphenyl)-1,3-thiazol-2-amine Chemical compound CC1=CC(C)=CC(C=2N=C(N)SC=2)=C1 MGWGWNFMUOTEHG-UHFFFAOYSA-N 0.000 claims description 4
- 238000011088 calibration curve Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000002485 combustion reaction Methods 0.000 claims description 4
- JCXJVPUVTGWSNB-UHFFFAOYSA-N nitrogen dioxide Inorganic materials O=[N]=O JCXJVPUVTGWSNB-UHFFFAOYSA-N 0.000 claims description 4
- 238000004445 quantitative analysis Methods 0.000 claims description 4
- 230000009257 reactivity Effects 0.000 claims description 4
- 239000002803 fossil fuel Substances 0.000 claims description 3
- 239000003960 organic solvent Substances 0.000 claims description 3
- 239000003973 paint Substances 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000011109 contamination Methods 0.000 claims 2
- 230000007613 environmental effect Effects 0.000 abstract description 3
- 238000011954 pollution control method Methods 0.000 abstract description 2
- WEVYAHXRMPXWCK-UHFFFAOYSA-N Acetonitrile Chemical compound CC#N WEVYAHXRMPXWCK-UHFFFAOYSA-N 0.000 description 6
- 239000000243 solution Substances 0.000 description 5
- 150000001299 aldehydes Chemical class 0.000 description 3
- MWUXSHHQAYIFBG-UHFFFAOYSA-N nitrogen oxide Inorganic materials O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 description 3
- 239000002243 precursor Substances 0.000 description 3
- 238000002076 thermal analysis method Methods 0.000 description 3
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- HSFWRNGVRCDJHI-UHFFFAOYSA-N alpha-acetylene Natural products C#C HSFWRNGVRCDJHI-UHFFFAOYSA-N 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 125000002534 ethynyl group Chemical group [H]C#C* 0.000 description 2
- 238000004128 high performance liquid chromatography Methods 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 238000012163 sequencing technique Methods 0.000 description 2
- HORQAOAYAYGIBM-UHFFFAOYSA-N 2,4-dinitrophenylhydrazine Chemical compound NNC1=CC=C([N+]([O-])=O)C=C1[N+]([O-])=O HORQAOAYAYGIBM-UHFFFAOYSA-N 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 239000000443 aerosol Substances 0.000 description 1
- 239000000926 atmospheric chemistry Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007664 blowing Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000013329 compounding Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 239000003480 eluent Substances 0.000 description 1
- 238000010828 elution Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000008014 freezing Effects 0.000 description 1
- 238000007710 freezing Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 150000002576 ketones Chemical class 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000004949 mass spectrometry Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 150000002894 organic compounds Chemical class 0.000 description 1
- 239000007800 oxidant agent Substances 0.000 description 1
- 239000013618 particulate matter Substances 0.000 description 1
- 238000006552 photochemical reaction Methods 0.000 description 1
- 238000004094 preconcentration Methods 0.000 description 1
- 238000005057 refrigeration Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000002904 solvent Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
- G01N30/02—Column chromatography
- G01N30/86—Signal analysis
- G01N30/8696—Details of Software
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/20—Identification of molecular entities, parts thereof or of chemical compositions
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Tourism & Hospitality (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computing Systems (AREA)
- Crystallography & Structural Chemistry (AREA)
- Strategic Management (AREA)
- Primary Health Care (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Economics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明涉及环境空气挥发性有机物监测领域，具体为一种基于监测数据及气象要素的VOCs污染来源解析方法，首先确定监测范围和监测因子；对环境空气中的VOCs进行定量分析；进行VOCs及组分浓度分布分析；评估各组分对臭氧生成的贡献；对VOCs与环境空气污染物和气象要素进行相关性分析；识别区域污染物排放源；采用正定矩阵因子分解模型进行VOCs污染来源解析。本发明可确定全检测区域以及各个检测点位全年VOCs污染浓度占比、时间变化趋势、各组分对臭氧生成贡献、VOCs的组分与气象要素间线性关系以及污染来源等数据，全面直观地分析污染来源；便于准确制定污染防治方法，建立VOCs污染防治体系。
Description
技术领域
本发明涉及环境空气挥发性有机物监测领域，具体为一种基于监测数据及气象要素的VOCs污染来源解析方法。
背景技术
挥发性有机物(Volatile Organic Compounds，VOCs)是一类有机物化合物的统称，对人体健康有较大危害。VOCs可以和氮氧化物发生光化学反应，形成光化学烟雾；也能与大气中的·OH、NO3 -、O3等氧化剂发生多途径发应，生成二次有机气溶胶(SOA)。挥发性有机物是形成臭氧污染的重要前体物，近年来以京津冀周边地区等为代表的城市大气污染问题日益突出，尤其臭氧污染态势未形成有效管控，成为仅次于颗粒物污染的第二大污染物。由VOCs引起的大气环境复合污染问题已经成为大气化学领域的研究热点。
挥发性有机物是形成细颗粒物(PM2.5)、臭氧(O3)等二次污染物的重要前体物，进而引发灰霾、光化学烟雾等大气环境问题。随着工业化和城市化的快速发展以及能源消费的持续增长，以PM2.5和臭氧为主要特征的区域性复合型大气污染日益突出，区域内空气重污染现象大范围同时出现的频次日益增多，由于VOCs排放来源复杂、排放形式多样、物质种类繁多，建立VOCs污染防治体系难度较大。
因此，亟需提供一种基于监测数据及气象要素的VOCs污染来源解析方法。
发明内容
为避免现有技术存在的上述问题，本发明的目的在于提供一种基于监测数据及气象要素的VOCs污染来源解析方法。
为实现上述目的，本发明提供了如下技术方案：一种基于监测数据及气象要素的VOCs污染来源解析方法，具体包括一下步骤：
S1：确定监测范围和监测因子；
S2：采用气质联用法或液相色谱法对环境空气中的VOCs进行定量分析；
S3：VOCs及组分浓度分布计算；
S4：利用最大增量反应，分析与评估各组分对臭氧生成的贡献；
S5：利用VOCs监测数据，对VOCs与主要污染物和主要气象要素进行相关性分析；
S6：利用不同挥发性有机物组分之间的相关性及特定物质的特征比值来识别区域污染物排放源，使用甲苯/苯比值，即T/B比值的大小来判断物污染物的来源；
S7：采用正定矩阵因子分解模型开展大气挥发性有机物VOCs污染来源解析研究。
本发明进一步设置为，步骤S1具体为，所述监测因子为VOCs中的70种目标物，即57种非甲烷烃类物质、13种醛酮类物质；
根据区域的主导风向在被检测区域内设置多个检测点位，在监测时段内，每隔6天进行一次连续采样，其中非甲烷烃类物质连续监测24小时，醛酮类在一天连续采样监测3小时。
本发明进一步设置为，步骤S2中，对非甲烷烃类物质用气质联用法作1.25nmol/mol、2.5nmol/mol、5nmol/mol、10nmol/mol、15nmol/mol、20nmol/mol不同浓度梯度的标准曲线；
对醛酮类物质采用液相色谱法，参考《环境空气醛、酮类化合物的测定高效液相色谱法(HJ 683-2014)》方法标准，配制标准系列分别为0.1μg/ml、0.2μg/ml、0.5μg/ml、1.0μg/ml、2.0μg/ml，校准曲线的相关系数≥0.999。
本发明进一步设置为，步骤S3具体为，根据监测结果，对VOCs的浓度水平、VOCs化学组分及VOCs关键物质进行污染趋势研究和特征分析；
所述化学成分即57种非甲烷烃类物质包含的烷烃、烯烃、炔烃和芳香烃4种化学组分共计57种目标物质，其中烷烃包含29种目标物质，烯烃包含11种目标物质，炔烃包含1种目标物质即乙炔，芳香烃包含16种目标物质。
本发明进一步设置为，所述污染趋势研究和特征分析，即确定全检测区域以及每个点位的非甲烷烃类物质、醛酮类物质的日浓度范围、平均浓度范围和时间变化；
所述VOCs化学组分的污染趋势研究和特征分析，即确定全检测区域以及每个点位芳香烃、炔烃、烷烃、烯烃分别的日浓度范围、平均浓度范围、平均浓度占比和时间变化；
所述VOCs关键物质的污染趋势研究和特征分析，即根据监测结果，将57种非甲烷烃类物质和13种醛酮类物质进行物质浓度大小排序，分析主要目标物质浓度水平占比和时间变化。
本发明进一步设置为，步骤S4具体为，各种污染组分的大气浓度与其最大增量反应性乘积为臭氧生产潜势OFP值，分别计算全检测区域以及每个点位的非甲烷烃类物质和醛酮类物质的OFP浓度范围和平均OFP浓度；根据计算结果，对VOCs的OFP特征、VOCs化学组分的OFP特征及VOCs关键物质的OFP特征进行分析，评估各组分对臭氧生成的贡献，组分的OFP值大，则说明该组分对臭氧生成潜势大，反之则小。
所述VOCs的OFP特征分析，即分别计算全检测区域以及每个点位的非甲烷烃类物质和醛酮类物质的OFP浓度范围和平均OFP浓度；确定每个点位的非甲烷烃类物质和醛酮类物质对全检测区域的臭氧生成贡献表现，以及全检测区域和各点位全年OFP浓度变化。
所述VOCs化学组分的OFP特征分析，即确定全检测区域以及每个点位芳香烃、炔烃、烷烃、烯烃分别的OFP浓度范围、OFP平均浓度、各组分平均占比，以及各组分全年OFP浓度变化。
所述VOCs关键物质的OFP特征分析，即将57种非甲烷烃类物质和13种醛酮类物质进行OFP浓度大小排序，分析主要目标物质OFP浓度水平占比即臭氧生成贡献和时间变化。
本发明进一步设置为，步骤S5具体为，步骤S6具体为，分别将非甲烷烃类物质和醛酮类物质的浓度水平与步骤S3中确定的环境空气主要污染物的浓度水平进行相关性分析，判断非甲烷烃类物质和醛酮类物质的浓度水平是否与环境空气主要污染物的浓度水平存在明显线性相关关系；所述环境空气污染物为臭氧、二氧化氮和PM2.5；
选取在VOCs的70种目标物质中选取浓度占比高的组分或关键物质与气象要素进行相关性分析；所述气象要素包括气压、气温、湿度、风速和风向；判断VOCs的目标物质是否与气象要素存在明显线性相关关系。
本发明进一步设置为，步骤S6具体为，分别计算全检测区域以及每个点位的T/B比值；若T/B比值大于2，说明污染物来源主要和涂料等有机溶剂有关；若T/B比值为2，说明污染物来源和交通机动车尾气排放有关；若T/B比值小于2，说明污染来源和石油化工生产、化石燃料燃烧有关。
本发明进一步设置为，步骤S7中，所述正定矩阵因子分解模型如下：
(CONC>MDL)
(CONC<MDL)
其中，xij为i样品中j组分的浓度；gik为第k个源对第i个样品的贡献；fkj为第k个排放源中j组分的含量；eij为残差；u为不确定度，uij为i样品中j组分的不确定度；Q为模型输入样品浓度及与对应物种相关的不确定度，通过迭代方法获取因子贡献和成分谱贡献的目标值Q，通过多次运行最终得到最稳健的Q值，MDL为检测限，EF为Error Fraction，可根据浓度与大小关系设为5％-20％，CONC为浓度。
依据正定矩阵因子分解PMF模型运算结果的信噪比S/N，对物质是否参与运算进行设置，其中S/N低于0.5的物质设置为bad进行剔除，S/N大于0.5小于1的，设置为weak以降低其权重，S/N大于1的则设置成strong；通过输入模型的数据进行分析，结合信噪比、物种相关性及物种时间序列浓度对物种及样本的可用性进行分析，运算给出模型结果。
需要注意的是，输入进模型的数据缺失值依据实际情况可指定为“-999”，可设定“-999”为物种特异性中位数或平均值代替。
本发明通过评价环境空气质量状况，统计分析濮阳市大气挥发性有机物的浓度水平、化学组分及关键物质污染特征，继而建立模型对影响臭氧生成的大气挥发性有机物开展潜势量计算分析，指出挥发性有机物中对臭氧生成影响较大的组分、种类及关键物质，结合气象因素进行相关性分析，从而明确对臭氧污染影响较大的挥发性有机物，分析与评估VOCs各组分对臭氧生成的贡献。
利用挥发性有机物监测数据，选择与之对应的臭氧、二氧化氮和PM2.5监测数据以及气象要素进行相关性分析，利用不同挥发性有机物组分之间的相关性及特定物质的特征比值来识别区域污染物排放源，利用PMF模型进行挥发性有机物污染来源解析，准确识别大气挥发性有机物的主要污染来源。
综上，本发明的上述技术方案的有益效果如下：
本发明根据检测区域内大气挥发性有机物监测结果，对VOCs的浓度水平、化学组分及关键物质进行污染特征分析，利用最大增量反应性，分析与评估VOCs各组分对臭氧生成的贡献，利用VOCs监测数据，选择与之对应的臭氧、二氧化氮和PM2.5监测数据以及气象要素进行相关性分析，利用不同VOCs组分之间的相关性及特定物质的特征比值来识别区域污染物排放源，利用PMF模型进行VOCs污染来源解析。
根据本发明所述方法，可得到全检测区域以及各个检测点位全年的VOCs污染浓度占比、时间变化趋势、各组分对臭氧生成的贡献、VOCs的主要组分与气象要素间线性关系以及污染来源等数据，全面、直观地分析VOCs污染来源；便于准确制定污染防治方法，建立VOCs污染防治体系。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明所述基于监测数据及气象要素的VOCs污染来源解析方法流程图。
图2为本发明所述13种醛酮类物质清单。
图3为本发明实施例中某市四点位均值各因子浓度及贡献率分布图。
具体实施方式
为了使本领域的人员更好地理解本发明的技术方案，下面结合本发明的附图，对本发明的技术方案进行清楚、完整的描述，基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的其它类同实施例，都应当属于本申请保护的范围。
下面结合附图和较佳的实施例对本发明作进一步说明。
实施例：
如图1所示，为本发明较佳实施例，一种基于监测数据及气象要素的VOCs污染来源解析方法，具体包括以下步骤：
S1：确定监测范围和监测因子；
所述监测因子为VOCs中的70种目标物，即57种非甲烷烃类物质、13种醛酮类物质；其中57种非甲烷烃类物质即原PAMS清单中物质，13种醛酮类物质见图2所示清单。
根据区域的主导风向在被检测区域内设置多个检测点位，在监测时段内，每隔6天进行一次连续采样，其中非甲烷烃类物质连续监测24小时，醛酮类在一天连续采样监测3小时。
S2：采用气质联用法或液相色谱法对环境空气中的VOCs进行定量分析；
对非甲烷烃类物质用气质联用法作1.25nmol/mol、2.5nmol/mol、5nmol/mol、10nmol/mol、15nmol/mol、20nmol/mol不同浓度梯度的标准曲线；
对醛酮类物质采用液相色谱法，参考《环境空气醛、酮类化合物的测定高效液相色谱法(HJ 683-2014)》方法标准，配制标准系列分别为0.1μg/ml、0.2μg/ml、0.5μg/ml、1.0μg/ml、2.0μg/ml，校准曲线的相关系数≥0.999。
其中57种PAMS使用一台1900自动采样器系统和3台RS1200+苏玛罐智能采样器10:00整点开始采样，每个样品采样时间为24小时，采样流量均为3.6mL/min，每10个样品采集一个全程序空白和一个样品平行。采集的样品通过大气预浓缩仪前处理，进入气相-质谱仪分析，根据校准曲线获得每个样品的VOCs浓度。而13种醛酮类物质用空气采样器通过DNPH+臭氧小柱串联采集，采样流量为1.0L/min，采样时间从中午12点到下午3点(共3小时)，每10个样品采集一个全程序空白和一个样品平行。
测量57种非甲烷烃类采用博赛德1900自动采样设备和7200大气预浓缩仪，该系统采用超低温预浓缩与气相色谱-质谱联用(GC-MS)检测技术，应用于痕量挥发性有机物测量，测量范围涵盖美国EPA光化学评估监测体系规定的57种臭氧前体物，13种含氧挥发性有机物。7200大气预浓缩仪通过为三级冷阱浓缩，大气样品被抽入仪器中，其中的痕量挥发性有机物分别被冷冻捕集下来，热解析后进入色谱柱，其中一路(C2-C5的碳氢化合物)由FID检测器检出，外一路(C5-C10的碳氢化合物和含氧挥发性有机物)由MS检测器检出。其中，冷冻富集装置采取电式制冷的方式，在-180℃低温下将目标化合物捕集，并在热解析的时候快速加热到80℃，随后样品进入GC-MS中分析。分离和检测采用安捷伦7890B/5977B气质联用仪。按照每次分析的时间顺序，一次完整的采样主要有4个步骤：样品采集、热解析、GC-MS分析和加热反吹。测量13种醛酮物质用乙腈按采样时气流相反的方向洗脱采样管，洗脱液收集至5mL容量瓶中，定容。用液相色谱仪进行分析，流动相为乙腈和水，采用梯度洗脱方法，流速1.0mL/min，柱温40℃，进样量20μL。
S3：VOCs及组分浓度分布计算；
根据监测结果，对VOCs的浓度水平、VOCs化学组分及VOCs关键物质进行污染趋势研究和特征分析；
所述化学成分即57种非甲烷烃类物质包含的烷烃、烯烃、炔烃和芳香烃4种化学组分共计57种目标物质，其中烷烃包含29种目标物质，烯烃包含11种目标物质，炔烃包含1种目标物质即乙炔，芳香烃包含16种目标物质。
所述VOCs的浓度水平的污染趋势研究和特征分析，即确定全检测区域以及每个点位非甲烷烃类物质、醛酮类物质的日浓度范围、平均浓度范围和时间变化；
所述VOCs化学组分的污染趋势研究和特征分析，即确定全检测区域以及每个点位芳香烃、炔烃、烷烃、烯烃分别的日浓度范围、平均浓度范围、平均浓度占比和时间变化；
所述VOCs关键物质的污染趋势研究和特征分析，即根据监测结果，将57种非甲烷烃类物质和13种醛酮类物质进行物质浓度大小排序，分析主要目标物质浓度水平占比和时间变化。
S4：利用最大增量反应，分析与评估各组分对臭氧生成的贡献；
各种污染组分的大气浓度与其最大增量反应性乘积为臭氧生产潜势OFP值，分别计算全检测区域以及每个点位的非甲烷烃类物质和醛酮类物质的OFP浓度范围和平均OFP浓度；根据计算结果，对VOCs的OFP特征、VOCs化学组分的OFP特征及VOCs关键物质的OFP特征进行分析，评估各组分对臭氧生成的贡献，组分的OFP值大，则说明该组分对臭氧生成潜势大，反之则小。
所述VOCs的OFP特征分析，即分别计算全检测区域以及每个点位的非甲烷烃类物质和醛酮类物质的OFP浓度范围和平均OFP浓度；确定每个点位的非甲烷烃类物质和醛酮类物质对全检测区域的臭氧生成贡献表现，以及全检测区域和各点位全年OFP浓度变化。
所述VOCs化学组分的OFP特征分析，即确定全检测区域以及每个点位芳香烃、炔烃、烷烃、烯烃分别的OFP浓度范围、OFP平均浓度、各组分平均占比，以及各组分全年OFP浓度变化。
所述VOCs关键物质的OFP特征分析，即将57种非甲烷烃类物质和13种醛酮类物质进行OFP浓度大小排序，分析主要目标物质OFP浓度水平占比即臭氧生成贡献和时间变化。
S5：利用VOCs监测数据，对VOCs与主要污染物和主要气象要素进行相关性分析；
分别将非甲烷烃类物质和醛酮类物质的浓度水平与步骤S3中确定的环境空气主要污染物的浓度水平进行相关性分析，判断非甲烷烃类物质和醛酮类物质的浓度水平是否与环境空气主要污染物的浓度水平存在明显线性相关关系；所述环境空气污染物为臭氧、二氧化氮和PM2.5；
选取在VOCs的70种目标物质中选取浓度占比高的组分或关键物质与气象要素进行相关性分析；所述气象要素包括气压、气温、湿度、风速和风向；判断VOCs的目标物质是否与气象要素存在明显线性相关关系。
S6：利用不同挥发性有机物组分之间的相关性及特定物质的特征比值来识别区域污染物排放源，使用甲苯/苯比值，即T/B比值的大小来判断物污染物的来源；
分别计算全检测区域以及每个点位的T/B比值；若T/B比值大于2，说明污染物来源主要和涂料等有机溶剂有关；若T/B比值为2，说明污染物来源和交通机动车尾气排放有关；若T/B比值小于2，说明污染来源和石油化工生产、化石燃料燃烧有关。
S7：采用正定矩阵因子分解模型开展大气挥发性有机物VOCs污染来源解析研究。
所述正定矩阵因子分解模型如下：
(CONC>MDL)
(CONC<MDL)
其中，xij为i样品中j组分的浓度；gik为第k个源对第i个样品的贡献；fkj为第k个排放源中j组分的含量；eij为残差；u为不确定度，uij为i样品中j组分的不确定度；Q为模型输入样品浓度及与对应物种相关的不确定度，通过迭代方法获取因子贡献和成分谱贡献的目标值Q，通过多次运行最终得到最稳健的Q值，MDL为检测限，EF为Error Fraction，可根据浓度与大小关系设为5％-20％，CONC为浓度。
依据正定矩阵因子分解PMF模型运算结果的信噪比S/N，对物质是否参与运算进行设置，其中S/N低于0.5的物质设置为bad进行剔除，S/N大于0.5小于1的，设置为weak以降低其权重，S/N大于1的则设置成strong；通过输入模型的数据进行分析，结合信噪比、物种相关性及物种时间序列浓度对物种及样本的可用性进行分析，运算给出模型结果。
需要注意的是，输入进模型的数据缺失值依据实际情况可指定为“-999”，可设定“-999”为物种特异性中位数或平均值代替。
需要说明的是，源解析数据共分为两组，分别为在线监测数据及手工监测数据，本实施例以某市检测数据为例，在该市设立四个检测点位，模型根据所选物质数量及多次运算尝试，共解析出四个点位均值大气挥发性有机物的主要源类5种，结果如图3所示。可见，该市VOCs来源解析共解析源5类，分别为机动车、天然源、燃烧源、溶剂使用及石化，其占比分别为41％、22.9％、17.4％、15.7％及3.1％。
以上所述仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (9)
1.一种基于监测数据及气象要素的VOCs污染来源解析方法，其特征在于，包括以下步骤：
S1：确定监测范围和监测因子；
S2：采用气质联用法或液相色谱法对环境空气中的VOCs进行定量分析；
S3：VOCs及组分浓度分布计算分析；
S4：利用最大增量反应，分析与评估各组分对臭氧生成的贡献；
S5：利用VOCs监测数据，对VOCs与环境空气污染物和气象要素进行相关性分析；
S6：利用不同VOCs组分之间的相关性及特定物质的特征比值来识别区域污染物排放源；
S7：采用正定矩阵因子分解模型进行VOCs污染来源解析研究。
2.根据权利要求1所述的一种基于监测数据及气象要素的VOCs污染来源解析方法，其特征在于，步骤S1具体为，所述监测因子为VOCs中的70种目标物，即57种非甲烷烃类物质和13种醛酮类物质；
根据区域的主导风向在被检测区域内设置多个检测点，在监测时段内，每隔6天进行一次连续采样，其中非甲烷烃类物质连续监测24小时，醛酮类在一天连续采样监测3小时。
3.根据权利要求1所述的一种基于监测数据及气象要素的VOCs污染来源解析方法，其特征在于，步骤S2中，对非甲烷烃类物质用气质联用法进行定量分析，作1.25nmol/mol、2.5nmol/mol、5nmol/mol、10nmol/mol、15nmol/mol、20nmol/mol不同浓度梯度的标准曲线；
对醛酮类物质采用液相色谱法进行定量分析，配制标准系列分别为0.1μg/ml、0.2μg/ml、0.5μg/ml、1.0μg/ml、2.0μg/ml，校准曲线的相关系数≥0.999。
4.根据权利要求1所述的一种基于监测数据及气象要素的VOCs污染来源解析方法，其特征在于，步骤S3具体为，根据监测结果，对VOCs的浓度水平、VOCs化学组分及VOCs关键物质进行污染趋势研究和特征分析；
所述化学成分即57种非甲烷烃类物质包含的烷烃、烯烃、炔烃和芳香烃4种化学组分共计57种目标物质。
5.根据权利要求4所述的一种基于监测数据及气象要素的VOCs污染来源解析方法，其特征在于，所述污染趋势研究和特征分析，即确定全检测区域以及每个点位的非甲烷烃类物质、醛酮类物质的日浓度范围、平均浓度范围和时间变化；
所述VOCs化学组分的污染趋势研究和特征分析，即确定全检测区域以及每个点位芳香烃、炔烃、烷烃、烯烃分别的日浓度范围、平均浓度范围、平均浓度占比和时间变化；
所述VOCs关键物质的污染趋势研究和特征分析，即根据监测结果，将57种非甲烷烃类物质和13种醛酮类物质进行物质浓度大小排序，分析主要目标物质浓度水平占比和时间变化。
6.根据权利要求1所述的一种基于监测数据及气象要素的VOCs污染来源解析方法，其特征在于，步骤S4具体为，各种污染组分的大气浓度与其最大增量反应性乘积为臭氧生产潜势OFP值，分别计算全检测区域以及每个点位的非甲烷烃类物质和醛酮类物质的OFP浓度范围和平均OFP浓度；根据计算结果，对VOCs的OFP特征、VOCs化学组分的OFP特征及VOCs关键物质的OFP特征进行分析，评估各组分对臭氧生成的贡献，组分的OFP值大，则说明该组分对臭氧生成潜势大，反之则小。
7.根据权利要求1所述的一种基于监测数据及气象要素的VOCs污染来源解析方法，其特征在于，步骤S5具体为，分别将非甲烷烃类物质和醛酮类物质的浓度水平与步骤S3中确定的环境空气污染物的浓度水平进行相关性分析，判断非甲烷烃类物质和醛酮类物质的浓度水平是否与环境空气污染物的浓度水平存在明显线性相关关系；所述环境空气污染物为臭氧、二氧化氮和PM2.5；
选取在VOCs的70种目标物质中选取浓度占比高的组分或关键物质与气象要素进行相关性分析；所述气象要素包括气压、气温、湿度、风速和风向；判断VOCs的目标物质是否与气象要素存在明显线性相关关系。
8.根据权利要求1所述的一种基于监测数据及气象要素的VOCs污染来源解析方法，其特征在于，步骤S6具体为，使用甲苯/苯比值，即T/B比值的大小来判断物污染物的来源；分别计算全检测区域以及每个点位的T/B比值；若T/B比值大于2，说明污染物来源和涂料等有机溶剂有关；若T/B比值为2，说明污染物来源和交通机动车尾气排放有关；若T/B比值小于2，说明污染来源和石油化工生产、化石燃料燃烧有关。
9.根据权利要求1所述的一种基于监测数据及气象要素的VOCs污染来源解析方法，其特征在于，步骤S7中，所述正定矩阵因子分解模型如下：
其中，xij为i样品中j组分的浓度；gik为第k个源对第i个样品的贡献；fkj为第k个排放源中j组分的含量；eij为残差；u为不确定度，uij为i样品中j组分的不确定度；Q为因子贡献和成分谱贡献的目标值Q；MDL为检测限，EF为Error Fraction，可根据浓度与大小关系设为5％-20％，CONC为浓度；
依据正定矩阵因子分解PMF模型运算结果的信噪比S/N，对物质是否参与运算进行设置，其中S/N低于0.5的物质设置为bad进行剔除，S/N大于0.5小于1的，设置为weak以降低其权重，S/N大于1的则设置成strong；通过输入模型数据进行分析，结合信噪比、物种相关性及物种时间序列浓度对物种及样本的可用性进行分析，运算给出模型结果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311291677.0A CN117316326A (zh) | 2023-10-08 | 2023-10-08 | 一种基于监测数据及气象要素的VOCs污染来源解析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311291677.0A CN117316326A (zh) | 2023-10-08 | 2023-10-08 | 一种基于监测数据及气象要素的VOCs污染来源解析方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117316326A true CN117316326A (zh) | 2023-12-29 |
Family
ID=89261771
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311291677.0A Pending CN117316326A (zh) | 2023-10-08 | 2023-10-08 | 一种基于监测数据及气象要素的VOCs污染来源解析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117316326A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117741016A (zh) * | 2024-02-20 | 2024-03-22 | 中汽研汽车检验中心（天津）有限公司 | 挥发性有机物逸出排放测试方法及装置 |
-
2023
- 2023-10-08 CN CN202311291677.0A patent/CN117316326A/zh active Pending
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117741016A (zh) * | 2024-02-20 | 2024-03-22 | 中汽研汽车检验中心（天津）有限公司 | 挥发性有机物逸出排放测试方法及装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Turpin et al. | Measuring and simulating particulate organics in the atmosphere: problems and prospects | |
| Materić et al. | Methods in plant foliar volatile organic compounds research | |
| Jiang et al. | Seasonal and diurnal variations of BTEX compounds in the semi-urban environment of Orleans, France | |
| Williams et al. | An in-situ instrument for speciated organic composition of atmospheric aerosols: T hermal desorption a erosol GC/MS-FID (TAG) | |
| Gros et al. | Volatile organic compounds sources in Paris in spring 2007. Part I: qualitative analysis | |
| Hopkins et al. | A two-column method for long-term monitoring of non-methane hydrocarbons (NMHCs) and oxygenated volatile organic compounds (o-VOCs) | |
| CN117316326A (zh) | 一种基于监测数据及气象要素的VOCs污染来源解析方法 | |
| Gensch et al. | Isotope ratio studies of atmospheric organic compounds: Principles, methods, applications and potential | |
| van Drooge et al. | Thermal desorption gas chromatography–mass spectrometry as an enhanced method for the quantification of polycyclic aromatic hydrocarbons from ambient air particulate matter | |
| Klemp et al. | High-speed GC analysis of VOCs: sample collection and inlet systems. 1. | |
| Liu et al. | Application of thermal desorption methods for airborne polycyclic aromatic hydrocarbon measurement: A critical review | |
| CN111665311A (zh) | 一种氢气中硫化物含量的检测方法 | |
| US8088627B2 (en) | On-line gas chromatographic analysis of airborne particles | |
| Rhee et al. | Continuous‐flow isotope analysis of the deuterium/hydrogen ratio in atmospheric hydrogen | |
| Liu et al. | Separation of PAHs in aerosol by thin layer chromatography for compound-specific stable carbon isotope analysis | |
| CN113466369A (zh) | 基于双冷阱预浓缩和飞行时间质谱的在线VOCs快速监测系统及其方法 | |
| Bandy et al. | Determination of part-per-trillion by volume levels of atmospheric carbon disulfide by gas chromatography/mass spectrometry | |
| Peng et al. | Real-time measurement of phase partitioning of organic compounds using a proton-transfer-reaction time-of-flight mass spectrometer coupled to a CHARON inlet | |
| CN218099031U (zh) | 一种环境大气中VOCs在线和离线快速监测装置 | |
| CN204479552U (zh) | 环境在线监测系统 | |
| Wernis et al. | Development of an in situ dual-channel thermal desorption gas chromatography instrument for consistent quantification of volatile, intermediate-volatility and semivolatile organic compounds | |
| Pierotti | Analysis of trace oxygenated hydrocarbons in the environment | |
| Whalen et al. | Detection of aromatic hydrocarbons in the atmosphere at ppt levels | |
| Cheng et al. | Analytical method development of long-chain ketones in PM2. 5 aerosols using accelerated solvent extraction and GC/FID/MSD | |
| Leinster et al. | Detection and measurement of volatile hydrocarbons at ambient concentrations in the atmosphere |
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 |